Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userSearch for tunnelling centres in Lennard-Jones clusters
65
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
We use molecular dynamics (MD) to simulate an unstable homogeneous mixture of binary fluids (AB), confined in a slit pore of width $D$. The pore walls are assumed to be flat and structureless, and attract one component of the mixture (A) with the same strength. The pair-wise interactions between the particles is modeled by the Lennard-Jones potential, with symmetric parameters that lead to a miscibility gap in the bulk. In the thin-film geometry, an interesting interplay occurs between surface enrichment and phase separation. We study the evolution of a mixture with equal amounts of A and B, which is rendered unstable by a temperature quench. We find that A-rich surface enrichment layers form quickly during the early stages of the evolution, causing a depletion of A in the inner regions of the film. These surface-directed concentration profiles propagate from the walls towards the center of the film, resulting in a transient layered structure. This layered state breaks up into a columnar state, which is characterized by the lateral coarsening of cylindrical domains. The qualitative features of this process resemble results from previous studies of diffusive Ginzburg-Landau-type models [S.~K. Das, S. Puri, J. Horbach, and K. Binder, Phys. Rev. E {bf 72}, 061603 (2005)], but quantitative aspects differ markedly. The relation to spinodal decomposition in a strictly 2-$d$ geometry is also discussed.
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 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.
The definitions of breaks and clusters in a one-dimensional chain in equilibrium are discussed. Analytical expressions are obtained for the expected cluster length, $langle K rangle$, as a function of temperature and pressure in a one-dimensional Lennard-Jones chain. These expressions are compared with results from molecular dynamics simulations. It is found that $langle K rangle$ increases exponentially with $beta = 1/k_BT$ and with pressure, $P$ in agreement with previous results in the literature. A method is illustrated for using $langle K rangle (beta, P)$ to generate a phase diagram for the Lennard-Jones chain. Some implications for the study of heat transport in Lennard-Jones chains are discussed.
This paper studies physical aging by computer simulations of a 2:1 Kob-Andersen binary Lennard-Jones mixture, a system that is less prone to crystallization than the standard 4:1 composition. Starting from thermal-equilibrium states, the time evolution of the following four quantities is monitored following up and down jumps in temperature: the potential energy, the virial, the average squared force, and the Laplacian of the potential energy. Despite the fact that significantly larger temperature jumps are studied here than in previous experiments, to a good approximation all four quantities conform to the single-parameter-aging scenario derived and validated for small jumps in experiments [Hecksher et al., J. Chem. Phys. 142, 241103 (2015)]. As a further confirmation of single-parameter aging with a common material time for the different quantities monitored, their relaxing parts are found to be almost identical for all temperature jumps.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
