No Arabic abstract
Scattering maps from strained or disordered nano-structures around a Bragg reflection can either be computed quickly using approximations and a (Fast) Fourier transform, or using individual atomic positions. In this article we show that it is possible to compute up to 4.10^10 $reflections.atoms/s using a single graphic card, and we evaluate how this speed depends on number of atoms and points in reciprocal space. An open-source software library (PyNX) allowing easy scattering computations (including grazing incidence conditions) in the Python language is described, with examples of scattering from non-ideal nanostructures.
The structural investigations of nanomaterials motivated by their large variety and diverse set of applications have attracted considerable attention. In particular, the ever-improving machinery, both in laboratory and at large scale facilities, together with the methodical improvements available for studying nanostructures ranging from epitaxial nanomaterials, nanocrystalline thin films and coatings, to nanoparticles and colloidal nanocrystals allows us to gain a more detailed understanding of their structural properties. As the structure essentially determines the physical properties of the materials, this advances the possibilities of structural studies and also enables a deeper understanding of the structure to property relationships. In this special issue entitled Investigation of Nanostructures with X-ray Scattering Techniques five contributions show the recent progress in various research fields. Contributions cover topics as diverse as neutron scattering on magnetic multilayer films, epitaxial orientation of organic thin films, nanoparticle ordering and chemical composition analysis, and the combination of nanofocused X-ray beams with electrical measurements.
Hidden Markov models (HMMs) are general purpose models for time-series data widely used across the sciences because of their flexibility and elegance. However fitting HMMs can often be computationally demanding and time consuming, particularly when the the number of hidden states is large or the Markov chain itself is long. Here we introduce a new Graphical Processing Unit (GPU) based algorithm designed to fit long chain HMMs, applying our approach to an HMM for nonvolcanic tremor events developed by Wang et al.(2018). Even on a modest GPU, our implementation resulted in a 1000-fold increase in speed over the standard single processor algorithm, allowing a full Bayesian inference of uncertainty related to model parameters. Similar improvements would be expected for HMM models given large number of observations and moderate state spaces (<80 states with current hardware). We discuss the model, general GPU architecture and algorithms and report performance of the method on a tremor dataset from the Shikoku region, Japan.
In this paper, we develop a highly efficient molecular dynamics code fully implemented on graphics processing units for thermal conductivity calculations using the Green-Kubo formula. We compare two different schemes for force evaluation, a previously used thread-scheme where a single thread is used for one particle and each thread calculates the total force for the corresponding particle, and a new block-scheme where a whole block is used for one particle and each thread in the block calculates one or several pair forces between the particle associated with the given block and its neighbor particle(s) associated with the given thread. For both schemes, two different classical potentials, namely, the Lennard-Jones potential and the rigid-ion potential are implemented. While the thread-scheme performs a little better for relatively large systems, the block-scheme performs much better for relatively small systems. The relative performance of the block-scheme over the thread-scheme also increases with the increasing of the cutoff radius. We validate the implementation by calculating lattice thermal conductivities of solid argon and lead telluride.
A computational fluid dynamics (CFD) simulation framework for predicting complex flows is developed on the Tensor Processing Unit (TPU) platform. The TPU architecture is featured with accelerated performance of dense matrix multiplication, large high bandwidth memory, and a fast inter-chip interconnect, which makes it attractive for high-performance scientific computing. The CFD framework solves the variable-density Navier-Stokes equation using a Low-Mach approximation, and the governing equations are discretized by a finite difference method on a collocated structured mesh. It uses the graph-based TensorFlow as the programming paradigm. The accuracy and performance of this framework is studied both numerically and analytically, specifically focusing on effects of TPU-native single precision floating point arithmetic on solution accuracy. The algorithm and implementation are validated with canonical 2D and 3D Taylor Green vortex simulations. To demonstrate the capability for simulating turbulent flows, simulations are conducted for two configurations, namely the decaying homogeneous isotropic turbulence and a turbulent planar jet. Both simulations show good statistical agreement with reference solutions. The performance analysis shows a linear weak scaling and a super-linear strong scaling up to a full TPU v3 pod with 2048 cores.
The correlations between the sequence of monomers in a polymer and its three-dimensional structure is a grand challenge in polymer science and biology. The properties and functions of macromolecules depend on their 3D shape that has appeared to be dictated by their monomer sequence. However, the progress towards understanding the sequence-structure-property correlations and their utilization in materials engineering are slow because it is almost impossible to characterize astronomically large number of possible sequences of a copolymer using traditional experimental and simulation methods. To address this problem, here, we combine evolutionary computing and coarse-grained molecular dynamics simulation and study the sequence-structure correlations of a model AB type copolymer system. The CGMD based evolutionary algorithm screens the sequence space of the copolymer efficiently and identifies wide range of single molecule structures including extremal radius of gyrations. The data provide new insights on the sequence-Rg correlations of the copolymer system and their impact on the structure and functionality of polymeric materials. The work highlights the opportunities of sequence specific control of macromolecular structure for designing materials with exceptional properties.