Quantum effects in the solid–liquid phase diagram of Ne13 and (para-H2)13

Using a reduced dimensionality model to compute the thermodynamic properties of finite polypeptide aggregates

By implementing a simple reduced dimensionality model to describe the interactions in finite systems composed of two seven-amino-acid peptides, the thermodynamic properties of ordered and disordered aggregates were computed. Within this model, the hydrophobicity of each amino acid was varied, and the stability of the systems computed. Accurate averages in the canonical ensemble were obtained using various replica exchange Monte Carlo algorithms. Low and high temperature regions were encountered where the ordered and disordered aggregates were stabilized. It was observed that as the degree of hydrophobicity increased, the stability of the aggregates increased, with a significant energetic stabilization obtained for the ordered aggregates. Upon decreasing the concentration of the solution, the stability of the amorphous aggregates increased when compared to the ordered systems.

Parallel tempering: theory, applications, and new perspectives

We review the history of the parallel tempering simulation method. From its origins in data analysis, the parallel tempering method has become a standard workhorse of physicochemical simulations. We discuss the theory behind the method and its various generalizations. We mention a selected set of the many applications that have become possible with the introduction of parallel tempering, and we suggest several promising avenues for future research.

Thermodynamics and equilibrium structure of Ne38 cluster: Quantum mechanics versus classical

The equilibrium properties of classical Lennard-Jones (LJ38) versus quantum Ne38 Lennard-Jones clusters are investigated. The quantum simulations use both the path-integral Monte Carlo (PIMC) and the recently developed variational-Gaussian wave packet Monte Carlo (VGW-MC) methods. The PIMC and the classical MC simulations are implemented in the parallel tempering framework. The classical heat capacity Cv(T) curve agrees well with that of Neirotti et al. [J. Chem. Phys. 112, 10340 (2000)], although a much larger confining sphere is used in the present work. The classical Cv(T) shows a peak at about 6 K, interpreted as a solid-liquid transition, and a shoulder at approximately 4 K, attributed to a solid-solid transition involving structures from the global octahedral (Oh) minimum and the main icosahedral (C5v) minimum. The VGW method is used to locate and characterize the low energy states of Ne38, which are then further refined by PIMC calculations. Unlike the classical case, the ground state of Ne38 is a liquidlike structure. Among the several liquidlike states with energies below the two symmetric states (Oh and C5v), the lowest two exhibit strong delocalization over basins associated with at least two classical local minima. Because the symmetric structures do not play an essential role in the thermodynamics of Ne38, the quantum heat capacity is a featureless curve indicative of the absence of any structural transformations. Good agreement between the two methods, VGW and PIMC, is obtained. The present results are also consistent with the predictions by Calvo et al. [J. Chem. Phys. 114, 7312 (2001)] based on the quantum superposition method within the harmonic approximation. However, because of its approximate nature, the latter method leads to an incorrect assignment of the Ne38 ground state as well as to a significant underestimation of the heat capacity.

Pressure dependent study of the solid-solid phase change in 38-atom Lennard-Jones cluster

Phase change phenomena in clusters are often modeled by augmenting physical interaction potentials with an external constraining potential to handle evaporation processes in finite temperature simulations. These external constraining potentials exert a pressure on the cluster. The influence of this constraining pressure on phase change phenomena in 38-atom Lennard-Jones clusters is investigated, and it is demonstrated that modest changes in the parameters of the constraining potential can lead to an order of magnitude change in the constraining pressure. At sufficiently high pressures the solid to solidlike phase change region in the 38-atom Lennard-Jones cluster is completely eliminated.

Phase changes in selected Lennard-Jones X13−nYn clusters

Detailed studies of the thermodynamic properties of selected binary Lennard-Jones clusters of the type X13-nYn (where n=1, 2, 3) are presented. The total energy, heat capacity, and first derivative of the heat capacity as a function of temperature are calculated by using the classical and path integral Monte Carlo methods combined with the parallel tempering technique. A modification in the phase change phenomena from the presence of impurity atoms and quantum effects is investigated.