Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userFusion process of Lennard-Jones clusters: global minima and magic numbers formation
77
0
0.0
(
0
)
Added by
Andrey V. Solov'yov
Publication date
2003
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
We present a new theoretical framework for modelling the fusion process of Lennard-Jones (LJ) clusters. Starting from the initial tetrahedral cluster configuration, adding new atoms to the system and absorbing its energy at each step, we find cluster growing paths up to the cluster sizes of up to 150 atoms. We demonstrate that in this way all known global minima structures of the LJ-clusters can be found. Our method provides an efficient tool for the calculation and analysis of atomic cluster structure. With its use we justify the magic number sequence for the clusters of noble gas atoms and compare it with experimental observations. We report the striking correspondence of the peaks in the dependence on cluster size of the second derivative of the binding energy per atom calculated for the chain of LJ-clusters based on the icosahedral symmetry with the peaks in the abundance mass spectra experimentally measured for the clusters of noble gas atoms. Our method serves an efficient alternative to the global optimization techniques based on the Monte-Carlo simulations and it can be applied for the solution of a broad variety of problems in which atomic cluster structure is important.
rate research
Read More
A relation $mathcal{M}_{mathrm{SHS}tomathrm{LJ}}$ between the set of non-isomorphic sticky hard sphere clusters $mathcal{M}_mathrm{SHS}$ and the sets of local energy minima $mathcal{M}_{LJ}$ of the $(m,n)$-Lennard-Jones potential $V^mathrm{LJ}_{mn}(r) = frac{varepsilon}{n-m} [ m r^{-n} - n r^{-m} ]$ is established. The number of nonisomorphic stable clusters depends strongly and nontrivially on both $m$ and $n$, and increases exponentially with increasing cluster size $N$ for $N gtrsim 10$. While the map from $mathcal{M}_mathrm{SHS}to mathcal{M}_{mathrm{SHS}tomathrm{LJ}}$ is non-injective and non-surjective, the number of Lennard-Jones structures missing from the map is relatively small for cluster sizes up to $N=13$, and most of the missing structures correspond to energetically unfavourable minima even for fairly low $(m,n)$. Furthermore, even the softest Lennard-Jones potential predicts that the coordination of 13 spheres around a central sphere is problematic (the Gregory-Newton problem). A more realistic extended Lennard-Jones potential chosen from coupled-cluster calculations for a rare gas dimer leads to a substantial increase in the number of nonisomorphic clusters, even though the potential curve is very similar to a (6,12)-Lennard-Jones potential.
There has long been a discrepancy between the size distributions of Ar$_n^+$ clusters measured by different groups regarding whether or not magic numbers appear at sizes corresponding to the closure of icosahedral (sub-)shells. We show that the previously observed magic cluster size distributions are likely the result of an unresolved Ar$_n$H$^+$ component, that is, from protonated argon clusters. We find that the proton impurity gives cluster geometries that are much closer to those for neutral rare gas clusters, which are known to form icosahedral structures, than the pure cationic clusters, explaining why the mass spectra from protonated argon clusters better matches these structural models. Our results thus show that even small impurities, e.g. a single proton, can significantly influence the properties of clusters.
We propose a method for quantifying charge-driven instabilities in clusters, based on equilibrium simulations under confinement at constant external pressure. This approach makes no assumptions about the mode of decay and allows different clusters to be compared on an equal footing. A comprehensive survey of stability in model clusters of 309 Lennard-Jones particles augmented with Coulomb interactions is presented. We proceed to examine dynamic signatures of instability, finding that rate constants for ejection of charged particles increase smoothly as a function of total charge with no sudden changes. For clusters where many particles carry charge, ejection of individual charges competes with a fission process that leads to more symmetric division of the cluster into large fragments. The rate constants for fission depend much more sensitively on total charge than those for ejection of individual particles.
The water-graphite interaction potential proposed recently (Gonzalez et al.emph{J. Phys. Chem. C} textbf{2007}, emph{111}, 14862), the three TIP$N$P ($N=3,:4,:5$) water-water interaction models, and basin-hopping global optimization are used to find the likely candidates for the global potential energy minima of (H$_{2}$O)$_{n}$ clusters with $nleq21$ on the (0001)-surface of graphite and to perform a comparative study of these minima. We show that, except for the smaller clusters ($n<6$), for which ab-initio results are available, the three water-water potential models provide mostly inequivalent conformations. While TIP3P seems to favor monolayer water structures for $n<18$, TIP4P and TIP5P favor bilayer or volume structures for $n>6$. These $n$ values determine the threshold of dominance of the hydrophobic nature of the water-graphite interaction at the nanoscopic scale for these potential models.
The fission of highly charged sodium clusters with fissilities X>1 is studied by {em ab initio} molecular dynamics. Na_{24}^{4+} is found to undergo predominantly sequential Na_{3}^{+} emission on a time scale of 1 ps, while Na_{24}^{Q+} (5 leq Q leq 8) undergoes multifragmentation on a time scale geq 0.1 ps, with Na^{+} increasingly the dominant fragment as Q increases. All singly-charged fragments Na_{n}^{+} up to size n=6 are observed. The observed fragment spectrum is, within statistical error, independent of the temperature T of the parent cluster for T leq 1500 K. These findings are consistent with and explain recent trends observed experimentally.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
