We develop an efficient numerical algorithm for the identification of a large number of saddle points of the potential energy function of Lennard- Jones clusters. Knowledge of the saddle points allows us to find many thousand adjacent minima of clusters containing up to 80 argon atoms and to locate many pairs of minima with the right characteristics to form two-level systems (TLS). The true TLS are singled out by calculating the ground-state tunneling splitting. The entropic contribution to all barriers is evaluated and discussed.
We report on numerical procedures for, and preliminary results on the search for, tunnelling centres in Lennard-Jones clusters, seen as simple model systems of glasses. Several of the double-well potentials identified are good candidates to give rise to two-level systems. The role of boundary effects, and the application of the semiclassical WKB approximation in multidimensional spaces for the calculation of the ground state splitting are discussed.
The phase diagram of the prototypical two-dimensional Lennard-Jones system, while extensively investigated, is still debated. In particular, there are controversial results in the literature as concern the existence of the hexatic phase and the melting scenario. Here, we study the phase behaviour of 2D LJ particles via large-scale numerical simulations. We demonstrate that at high temperature, when the attraction in the potential plays a minor role, melting occurs via a continuous solid-hexatic transition followed by a first-order hexatic-fluid transition. As the temperature decreases, the density range where the hexatic phase occurs shrinks so that at low-temperature melting occurs via a first-order liquid-solid transition. The temperature where the hexatic phase disappears is well above the liquid-gas critical temperature. The evolution of the density of topological defects confirms this scenario.
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.
We present a systematic study of the thermodynamics of two and three-dimensional generalized Lennard-Jones ($LJ$) systems focusing on the relationship between the range of the potential, the system density and its dimension. We found that the existence of negative specific heats depends on these three factors and not only on the potential range and the density of the system as stated in recent contributions.
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.