Molecular dynamic simulations of atom–cluster collision processes

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.

Exploring the importance of quantum effects in nucleation: the archetypical Ne(n) case

The effect of quantum mechanics (QM) on the details of the nucleation process is explored employing Ne clusters as test cases due to their semi-quantal nature. In particular, we investigate the impact of quantum mechanics on both condensation and dissociation rates in the framework of the microcanonical ensemble. Using both classical trajectories and two semi-quantal approaches (zero point averaged dynamics, ZPAD, and Gaussian-based time dependent Hartree, G-TDH) to model cluster and collision dynamics, we simulate the dissociation and monomer capture for Ne(8) as a function of the cluster internal energy, impact parameter and collision speed. The results for the capture probability P(s)(b) as a function of the impact parameter suggest that classical trajectories always underestimate capture probabilities with respect to ZPAD, albeit at most by 15%-20% in the cases we studied. They also do so in some important situations when using G-TDH. More interestingly, dissociation rates k(diss) are grossly overestimated by classical mechanics, at least by one order of magnitude. We interpret both behaviours as mainly due to the reduced amount of kinetic energy available to a quantum cluster for a chosen total internal energy. We also find that the decrease in monomer dissociation energy due to zero point energy effects plays a key role in defining dissociation rates. In fact, semi-quantal and classical results for k(diss) seem to follow a common "corresponding states" behaviour when the proper definition of internal and dissociation energies are used in a transition state model estimation of the evaporation rate constants.

A comparison of rigid and flexible water models in collisions of monomers and small clusters

In this study we have investigated the dynamics of small water clusters using microcanonical molecular dynamics simulations. The clusters are formed by colliding vapor monomers with target clusters of two and five molecules. The monomers are sampled from a thermal ensemble at T=300 K and target clusters with several total energies are considered. We compare rigid extended simple point charge water with flexible counterparts having intramolecular harmonic bonds with force constants 10(3) and 10(5) kcal(mol A2). We show that the lifetimes of the clusters formed via collision process are similar for the rigid model and the flexible model with the bigger force constant, if the translational temperatures of the target cluster molecules are equal. The model with the smaller force constant results in much longer lifetimes due to the stabilizing effect caused by the kinetic energy transfer into internal vibration of the molecules. This process may take several hundreds of picoseconds, giving rise to time-dependent decay rates of constant-energy clusters. A study of binary collisions of water molecules shows that the introduction of flexibility to the molecules increases the possibility of dimer formation and thus offers an alternative route for dimer production in vapors. Our results imply that allowing for internal degrees of freedom is likely to enhance gas-liquid nucleation rates in water simulations.

Postcollision relaxation of small atomic clusters

Molecular-dynamics simulations are performed to investigate the effects caused by the lack of internal equilibration on the dynamics and properties of atomic clusters. The studied systems consist of Lennard-Jones clusters of five to ten atoms and a colliding vapor monomer. Cluster radius and potential energy are shown to reach a time-independent value within 30 ps after a collision with a vapor monomer. The relaxation in terms of rotational energy takes at least 200 ps. During the first couple of picoseconds after the collision time-dependent cluster decay rates are observed. The unrelaxed cluster states are expected to have minimal effect on gas-liquid nucleation rates.

Molecular dynamics simulations of cluster nucleation during inert gas condensation

Molecular dynamics simulations of vapor-phase nucleation of germanium in an argon atmosphere were performed and a unexpected channel of nucleation was observed. This channel, vapor-induced cluster splitting, is important for more refractory materials since the critical nucleus size can fall below the size of a dimer. As opposed to conventional direct vapor nucleation of the dimer, which occurs by three-body collisions, cluster-splitting nucleation is a second-order reaction. The most important cluster-splitting reaction is the collision of a vapor atom and a trimer that leads to the formation of two dimers. The importance of the cluster-splitting nucleation channel relative to the direct vapor nucleation channel is observed to increase with decreasing vapor density and increasing ratio of vapor to carrier gas atoms.

The role of dimers in evaporation of small argon clusters

Evaporation of small Lennard-Jones argon clusters has been studied using molecular dynamic simulations. An extensive library of clusters with 4, 5, 6, 11, and 21 atoms has been obtained from an earlier study. Analysis of the evaporation properties of the clusters indicate, that the fraction of dimer evaporations of all evaporation events increases with the total energy of the cluster. The fraction of evaporated dimers from clusters with a constant lifetime is independent of the cluster size for short-lived clusters and increases with cluster size for long-lived clusters. Only a few percent of the clusters which are long lived enough to participate in vapor-liquid nucleation decay by emitting dimers. The mean cluster lifetime as a function of total energy shows the same exponentially decreasing trend for monomer and dimer evaporation channels. The fraction of trimer evaporations is found to be vanishingly small.