Molecular dynamics simulations of cluster nucleation during inert gas condensation

Study on liquid-like SiGe cluster growth during co-condensation from supersaturated vapor mixtures by molecular dynamics simulation

Based on the co-condensation processes in the Si-Ge system upon cooling, as determined by molecular dynamics (MD) simulation, we explored the mixed cluster growth dynamics and structural properties leading to the synthesis of liquid-like SiGe nanoclusters. The results indicated that the cluster size quickly increased to large clusters by the coalescence of transient small clusters in the growth stage during co-condensation. The transient clusters at different temperatures were verified to have slightly Si-rich compositions and liquid-like structures. The coalescence of such nanoclusters at high temperatures led to spherical clusters with homogeneous intermixing. However, irregularly shaped clusters with attached mixed parts were obtained owing to incomplete coalescence at low temperatures. Ge atoms tended to move to the cluster surface to exploit their energetically favorable state during the restructuring process, leading to slightly Ge-rich components on the cluster surface. The degree of intermixing for SiGe nanoclusters was related to cluster size. Generally, small clusters appeared to be more segregated during restructuring.

Alloy, Janus and core-shell nanoparticles: numerical modeling of their nucleation and growth in physical synthesis

While alloy, core-shell and Janus binary nanoclusters are found in more and more technological applications, their formation mechanisms are still poorly understood, especially during synthesis methods involving physical approaches. In this work, we employ a very simple model of such complex systems using Lennard-Jones interactions and inert gas quenching. After demonstrating the ability of the model to well reproduce the formation of alloy, core-shell or Janus nanoparticles, we studied their temporal evolution from the gas via droplets to nanocrystalline particles. In particular, we showed that the growth mechanisms exhibit qualitative differences between these three chemical orderings. Then, we determined how the quenching rate can be used to finely tune structural characteristics of the final nanoparticles, including size, shape and crystallinity.

Predicting stability limits for pure and doped dicationic noble gas clusters undergoing coulomb explosion: A parallel tempering based study

We have used a replica exchange Monte-Carlo procedure, popularly known as Parallel Tempering, to study the problem of Coulomb explosion in homogeneous Ar and Xe dicationic clusters as well as mixed Ar-Xe dicationic clusters of varying sizes with different degrees of relative composition. All the clusters studied have two units of positive charges. The simulations reveal that in all the cases there is a cutoff size below which the clusters fragment. It is seen that for the case of pure Ar, the value is around 95 while that for Xe it is 55. For the mixed clusters with increasing Xe content, the cutoff limit for suppression of Coulomb explosion gradually decreases from 95 for a pure Ar to 55 for a pure Xe cluster. The hallmark of this study is this smooth progression. All the clusters are simulated using the reliable potential energy surface developed by Gay and Berne (Gay and Berne, Phys. Rev. Lett. 1982, 49, 194). For the hetero clusters, we have also discussed two different ways of charge distribution, that is one in which both positive charges are on two Xe atoms and the other where the two charges are at a Xe atom and at an Ar atom. The fragmentation patterns observed by us are such that single ionic ejections are the favored dissociating pattern. © 2017 Wiley Periodicals, Inc.

Au–Ir nanoalloy nucleation during the gas-phase condensation: a comprehensive MD study including different effects

The number of formed clusters and their size increases with the increasing temperature and pressure, which is in good agreement with the experimental results.

A parallel tempering based study of Coulombic explosion and identification of dissociating fragments in charged noble gas clusters

In this communication, we would like to test the feasibility of a parallel tempering based study of dissociation in dicationic noble gas clusters, namely, Ar(n)(2+), Kr(n)(2+), and Xe(n)(2+), where "n" is the size of the cluster units. We would like to find out the correct limit for sizes of each of these systems, above which the clusters stay intact as a single unit and does not dissociate into fragments by the process of Coulomb explosion. Moreover, we would also like to, for a specific case, i.e., Ar(n)(2+), study in detail the fragmentation patterns and point out the switchover from the non-fission way to the fission mechanism of dissociation. In all these calculations, we would like to analyse, how close we are in our predictions with that of experimental results. As a further check on the dissociating patterns found out by parallel tempering, we also conduct basin hopping based study on representative sizes of the clusters and find that parallel tempering, as used for this present work as an optimizer, is able to predict correct features when compared with other celebrated methods like the basin hopping algorithm.

Comparative study on methodology in molecular dynamics simulation of nucleation

Gas-liquid nucleation of 1000 Lennard-Jones atoms is simulated to evaluate temperature regulation methods and methods to obtain nucleation rate. The Berendsen and the Andersen thermostats are compared. The Berendsen thermostat is unable to control the temperature of clusters larger than the critical size. Independent of the thermostating method the velocities of individual atoms and the translational velocities of clusters up to at least six atoms are accurately described by the Maxwell velocity distribution. Simulations with the Andersen thermostat yield about two times higher nucleation rates than those with the Berendsen thermostat. Nucleation rate is extracted from the simulations by direct observation of times of nucleation onset and by the method of Yasuoka and Matsumoto [J. Chem. Phys. 109, 8451 (1998)]. Compared to the direct observation, the nucleation rates obtained from the method of Yasuoka and Matsumoto are higher by a factor of 3.

Comparison between the classical theory predictions and molecular simulation results for heterogeneous nucleation of argon

We have performed Monte Carlo simulations of homogeneous and heterogeneous nucleations of Lennard-Jones argon clusters. The simulation results were interpreted using the major concept posing a difference between the homogeneous and heterogeneous classical nucleation theories-the contact parameter. Our results show that the multiplication concept of the classical heterogeneous nucleation theory describes the cluster-substrate interaction surprisingly well even for small molecular clusters. However, in the case of argon nucleating on a rigid monolayer of fcc(111) substrate at T=60 K, the argon-substrate atom interaction being approximately one-third as strong as the argon-argon interaction, the use of the classical theory concept results in an underestimation of the heterogeneous nucleation rate by two to three orders of magnitude even for large clusters. The main contribution to this discrepancy is induced by the failure of the classical theory of homogeneous nucleation to predict the energy involved in bringing one molecule from the vapor to the cluster for clusters containing less than approximately 15 molecules.

Molecular-dynamics study of the density scaling of inert gas condensation

The initial stages of vapor condensation of Ge in the presence of a cold Ar atmosphere were studied by molecular-dynamics simulations. The state variables of interest included the densities of condensing vapor and gas, the density of clusters, and the average cluster size, while the temperatures of the vapor and the clusters were separately monitored with time. Three condensation processes were explicitly identified: nucleation, monomeric growth, and cluster aggregation. Our principal finding is that both the average cluster size and the number of clusters scale with the linear dimension of the computation cell, L, and Ln, with the scaling parameter n approximately 4, corresponding to a reaction order of nu approximately 2.33. This small value of n is explained by an unexpected nucleation path involving the formation of Ge dimers via two-body collisions.