Three parallel algorithms for classical molecular dynamics are presented. The first assigns each processor a fixed subset of atoms; the second assigns each a fixed subset of inter-atomic forces to compute; the third assigns each a fixed spatial region. The algorithms are suitable for molecular dynamics models which can be difficult to parallelize efficiently—those with short-range forces where the neighbors of each atom change rapidly. They can be implemented on any distributed-memory parallel machine which allows for message-passing of data between independently executing processors. The algorithms are tested on a standard Lennard-Jones benchmark problem for system sizes ranging from 500 to 100,000,000 atoms on several parallel supercomputers--the nCUBE 2, Intel iPSC/860 and Paragon, and Cray T3D. Comparing the results to the fastest reported vectorized Cray Y-MP and C90 algorithm shows that the current generation of parallel machines is competitive with conventional vector supercomputers even for small problems. For large problems, the spatial algorithm achieves parallel efficiencies of 90% and a 1840-node Intel Paragon performs up to 165 faster than a single Cray C9O processor. Trade-offs between the three algorithms and guidelines for adapting them to more complex molecular dynamics simulations are also discussed.

Fast parallel algorithms for short-range molecular dynamics

The nucleation of crystals in liquids is one of nature’s most ubiquitous phenomena, playing an important role in areas such as climate change and the production of drugs. As the early stages of nucleation involve exceedingly small time and length scales, atomistic computer simulations can provide unique insights into the microscopic aspects of crystallization. In this review, we take stock of the numerous molecular dynamics simulations that, in the past few decades, have unraveled crucial aspects of crystal nucleation in liquids. We put into context the theoretical framework of classical nucleation theory and the state-of-the-art computational methods by reviewing simulations of such processes as ice nucleation and the crystallization of molecules in solutions. We shall see that molecular dynamics simulations have provided key insights into diverse nucleation scenarios, ranging from colloidal particles to natural gas hydrates, and that, as a result, the general applicability of classical nucleation theory...

Crystal Nucleation in Liquids: Open Questions and Future Challenges in Molecular Dynamics Simulations.

The aggregative growth and oriented attachment of nanocrystals and nanoparticles are reviewed, and they are contrasted to classical LaMer nucleation and growth, and to Ostwald ripening. Kinetic and mechanistic models are presented, and experiments directly observing aggregative growth and oriented attachment are summarized. Aggregative growth is described as a nonclassical nucleation and growth process. The concept of a nucleation function is introduced, and approximated with a Gaussian form. The height (Γmax) and width (Δtn) of the nucleation function are systematically varied by conditions that influence the colloidal stability of the small, primary nanocrystals participating in aggregative growth. The nucleation parameters Γmax and Δtn correlate with the final nanocrystal mean size and size distribution, affording a potential means of achieving nucleation control in nanocrystal synthesis.

Kinetics and Mechanisms of Aggregative Nanocrystal Growth

The nucleation of crystals in liquids is one of nature's most ubiquitous phenomena, playing an important role in areas such as climate change and the production of drugs. As the early stages of nucleation involve exceedingly small time and length scales, atomistic computer simulations can provide unique insight into the microscopic aspects of crystallization. In this review, we take stock of the numerous molecular dynamics simulations that in the last few decades have unraveled crucial aspects of crystal nucleation in liquids. We put into context the theoretical framework of classical nucleation theory and the state of the art computational methods, by reviewing simulations of e.g. ice nucleation or crystallization of molecules in solutions. We shall see that molecular dynamics simulations have provided key insight into diverse nucleation scenarios, ranging from colloidal particles to natural gas hydrates, and that in doing so the general applicability of classical nucleation theory has been repeatedly called into question. We have attempted to identify the most pressing open questions in the field. We believe that by improving (i.) existing interatomic potentials; and (ii.) currently available enhanced sampling methods, the community can move towards accurate investigations of realistic systems of practical interest, thus bringing simulations a step closer to experiments.

Crystal Nucleation in Liquids: Open Questions and Future Challenges in Molecular Dynamics Simulations

Can completely homogeneous nucleation occur? Large scale molecular dynamics simulations performed on a graphics-processing-unit rich supercomputer can shed light on this long-standing issue. Here, a billion-atom molecular dynamics simulation of homogeneous nucleation from an undercooled iron melt reveals that some satellite-like small grains surrounding previously formed large grains exist in the middle of the nucleation process, which are not distributed uniformly. At the same time, grains with a twin boundary are formed by heterogeneous nucleation from the surface of the previously formed grains. The local heterogeneity in the distribution of grains is caused by the local accumulation of the icosahedral structure in the undercooled melt near the previously formed grains. This insight is mainly attributable to the multi-graphics processing unit parallel computation combined with the rapid progress in high-performance computational environments.Nucleation is a fundamental physical process, however it is a long-standing issue whether completely homogeneous nucleation can occur. Here the authors reveal, via a billion-atom molecular dynamics simulation, that local heterogeneity exists during homogeneous nucleation in an undercooled iron melt.

/pdf/heterogeneity-in-homogeneous-nucleation-from-billion-atom-2thbquxnx8.pdf

Heterogeneity in homogeneous nucleation from billion-atom molecular dynamics simulation of solidification of pure metal

We present a novel method to identify local structures in disordered systems according to topological criteria. Its effectiveness is demonstrated in the analysis of the atomic structures in the rapid cooling of silver liquid. The method is parameter free and scale independent, and can generally be used for structural analysis of amorphous systems involving atoms or particles at different length scales.

https://iopscience.iop.org/article/10.1209/0295-5075/96/36001/pdf

A new method for analyzing the local structures of disordered systems

A molecular dynamics simulation study has been performed to investigate the solid microstructures of Ag at room temperature resulted from rapid solidification with six cooling rates by using Quantum Sutton–Chen (QSC) many body potential. Several structural analysis methods have been employed to characterize and analyze the atomic configurations in the system. The results reveal that as the cooling rate γ ⩾ 1.0 × 10 13  K/s, an amorphous structure can be obtained; as the cooling rates in the range of γ = 5.0 × 10 12  ∼ 1.38 × 10 11  K/s, the coexistence structures of the meta-stable hcp with stable fcc configurations can be formed with phase separating or layering structures. Moreover, by using the 3D graphics technology, the arrangement of atoms in the system reveals that the phase separating occurs before the layering structures can be formed following the cooling rate decreasing. And an interesting layering structures have been found, in which the central atoms of fcc and hcp basic clusters are arranged by alternate close-packed-layers with a new repeating sequence ABCB.

Molecular dynamics simulation for cooling rate dependence of solidification microstructures of silver

The effect of the cooling rate on the solidification process of liquid aluminium is studied using a large-scale molecular dynamics method. It is found that there are various types of short-range order (SRO) structures in the liquid, among which the icosahedral (ICO)-like structures are dominant. These SRO structures are in dynamic fluctuation and transform each other. The effect of the cooling rate on the microstructure is very weak at high temperatures and in supercooled liquids, and it appears only below the liquid–solid transition temperature. Fast cooling rates favour the formation of amorphous structures with ICO-like features, while slow cooling rates favour the formation of FCC crystalline structures. Furthermore, FCC and HCP structures can coexist in crystalline structures. It is also found that nanocrystalline aluminium can be achieved at appropriate cooling rates, and its formation mechanism is thoroughly investigated by tracing the evolution of nanoclusters. The arrangement of FCC and HCP atoms in the nanograins displays various twinned structures as observed using visualization analysis, which is different from the layering or phase separation structures observed in the solidification of Lennard-Jones fluids and some metal liquids.

Cooling rate dependence of solidification for liquid aluminium: a large-scale molecular dynamics simulation study

In order to understand the formation mechanism of nano-clusters during rapid solidification process in liquid metals, a molecular dynamics (MD) simulation study for a large-scale system consisting of 400 000 atoms of liquid metal Al has been performed. To describe the clusters, especially the nano-clusters formed in the system during rapid solidification processes, a cluster-type index method (CTIM) has been proposed. It is demonstrated that the nano-clusters (such as containing 150 atoms) are formed by combining some middle clusters and they sharply differ from those obtained by gaseous deposition, ionic spray methods and so on. Though the nano-clusters have different shapes, they all have the corners that can be starting points of dendrite growth in solidification processes in liquid metals.

Formation and description of nano-clusters formed during rapid solidification processes in liquid metals

The structural evolution of tantalum (Ta) during rapid cooling was investigated by molecular dynamics simulation, in terms of the system energy, the pair distribution function and the largest standard cluster analysis. It was found that the critical cooling rate for vitrification was about R ≥ 0.25 K ps−1, two orders lower than other metals (such as Au, Ag, Al, Zr and Zn) and that the meta-stable σ phase (β-Ta) not only appears on the pathway from liquid to the stable body-centred cubic crystal, but is also easily obtained at room temperature as a long-lived metastable phase with some probability. The most interesting point is that the liquid, amorphous and β-Ta phases share a nontrivial structural homology; the intrinsic topologically close-packed (TCP) structures in liquids are inherited and developed in different ways, resulting in amorphous or crystalline solids, respectively. With highly local packing fractions and geometrical incompatibility with the global close-packed (such as hcp, fcc and bcc) crystals, TCP structures inevitably result in structural heterogeneity and favour vitrification. As a superset of icosahedrons, TCP structures are ubiquitous in metallic melts, and just before the onset of crystallization reach their maximal number, which is much bigger in Ta than in other poor-GFA metals; so we argue that the strong forming ability of TCP local structures significantly enhances the glass forming ability of pure metals. These findings open up a new perspective that could have a profound impact on the research into metallic glasses.

Rang-su Liu

Papers

A new method for analyzing the local structures of disordered systems

Molecular dynamics simulation for cooling rate dependence of solidification microstructures of silver

Cooling rate dependence of solidification for liquid aluminium: a large-scale molecular dynamics simulation study

Formation and description of nano-clusters formed during rapid solidification processes in liquid metals

Topologically close-packed characteristic of amorphous tantalum.