5,195 research outputs found
Density of states of a binary Lennard-Jones Glass
We calculate the density of states of a binary Lennard-Jones glass using a
recently proposed Monte Carlo algorithm. Unlike traditional molecular
simulation approaches, the algorithm samples distinct configurations according
to self-consistent estimates of the density of states, thereby giving rise to
uniform internal-energy histograms. The method is applied to simulate the
equilibrium, low-temperature thermodynamic properties of a widely studied glass
former consisting of a binary mixture of Lennard-Jones particles. We show how a
density-of-states algorithm can be combined with particle identity swaps and
configurational bias techniques to study that system. Results are presented for
the energy and entropy below the mode coupling temperature.Comment: 6 pages, 3 figures, accepted by J Chem Phy
Constant Pressure Hybrid Molecular Dynamics-Monte Carlo Simulations
New hybrid Molecular Dynamics-Monte Carlo methods are proposed to increase
the efficiency of constant-pressure simulations. Two variations of the isobaric
Molecular Dynamics component of the algorithms are considered. In the first, we
use the extended-ensemble method of Andersen [H. C. Andersen J. Chem. Phys.
{\bf 72},2384 (1980)]. In the second, we arrive at a new constant-pressure
Monte Carlo technique based on the reversible generalization of the
weak-coupling barostat [H. J. C. Berendsen et. al J. Chem. Phys. {\bf 81},
3684(1984)]. This latter technique turns out to be highly effective in
equilibrating and maintaining a target pressure. It is superior to the
extended-ensemble method, which in turn is superior to simple volume-rescaling
algorithms. The efficiency of the proposed methods is demonstrated by studying
two systems. The first is a simple Lennard-Jones fluid. The second is a mixture
of polyethylene chains of 200 monomers.Comment: 10 pages, 4 figure
Density of States Monte Carlo Method for Simulation of Fluids
A Monte Carlo method based on a density-of-states sampling is proposed for
study of arbitrary statistical mechanical ensembles in a continuum. A random
walk in the two-dimensional space of particle number and energy is used to
estimate the density of states of the system; this density of states is
continuously updated as the random walk visits individual states. The validity
and usefulness of the method are demonstrated by applying it to the simulation
of a Lennard-Jones fluid. Results for its thermodynamic properties, including
the vapor-liquid phase coexistence curve, are shown to be in good agreement
with high-accuracy literature data.Comment: 5 pages, 3 figures, accepted by J Chem Phy
Multicanonical Parallel Tempering
We present a novel implementation of the parallel tempering Monte Carlo
method in a multicanonical ensemble. Multicanonical weights are derived by a
self-consistent iterative process using a Boltzmann inversion of global energy
histograms. This procedure gives rise to a much broader overlap of
thermodynamic-property histograms; fewer replicas are necessary in parallel
tempering simulations, and the acceptance of trial swap moves can be made
arbitrarily high. We demonstrate the usefulness of the method in the context of
a grand-multicanonical ensemble, where we use multicanonical simulations in
energy space with the addition of an unmodified chemical potential term in
particle-number space. Several possible implementations are discussed, and the
best choice is presented in the context of the liquid-gas phase transition of
the Lennard-Jones fluid. A substantial decrease in the necessary number of
replicas can be achieved through the proposed method, thereby providing a
higher efficiency and the possibility of parallelization.Comment: 8 pages, 3 figure, accepted by J Chem Phy
Low-temperature anomalies of a vapor deposited glass
We investigate the low temperature properties of two-dimensional
Lennard-Jones glass films, prepared in silico both by liquid cooling and by
physical vapor deposition. We identify deep in the solid phase a crossover
temperature , at which slow dynamics and enhanced heterogeneity emerge.
Around , localized defects become visible, leading to vibrational
anomalies as compared to standard solids. We find that on average,
decreases in samples with lower inherent structure energy, suggesting that such
anomalies will be suppressed in ultra-stable glass films, prepared both by very
slow liquid cooling and vapor deposition.Comment: 10 pages including appendices, 8 figures. Version accepted for
Physical Review Material
- …