Van der Waals type loop in microcanonical caloric curves of finite systems

Where macro meets micro

Reconciling or somehow linking the macroscopic and microscopic approaches to chemical and physical processes has been a challenge unaddressed for many years. One approach, presented here, treats the issue by examining individual phenomena well described by a macro approach that fails when applied to small systems. The key to the approach is determining the approximate system size below which the breakdown of the macro description is observable. The most developed example is the failure of the Gibbs phase rule for sufficiently small atomic clusters. Other examples, such as the onset, at sufficient size, of the insulator-to-metal transition, are discussed, as are some still more challenging phenomena.

Generalized simulated tempering realized on expanded ensembles of non-Boltzmann weights

A generalized version of the simulated tempering operated in the expanded ensembles of non-Boltzmann weights has been proposed to mitigate a quasiergodicity problem occurring in simulations of rough energy landscapes. In contrast to conventional simulated tempering employing the Boltzmann weight, our method utilizes a parametrized, generalized distribution as a workhorse for stochastic exchanges of configurations and subensembles transitions, which allows a considerable enhancement for the rate of convergence of Monte Carlo and molecular dynamics simulations using delocalized weights. A feature of our method is that the exploration of the parameter space encouraging subensembles transitions is greatly accelerated using the dynamic update scheme for the weight via the average guide specific to the energy distribution. The performance and characteristic feature of our method have been validated in the liquid-solid transition of Lennard-Jones clusters and the conformational sampling of alanine dipeptide by taking two types of Tsallis [C. Tsallis, J. Stat. Phys. 52, 479 (1988)] expanded ensembles associated with different parametrization schemes.

Dynamical origin of enhanced conformational searches of Tsallis statistics sampling

The characteristic sampling dynamics of importance samplings driven by the Tsallis weight [C. Tsallis, J. Stat. Phys. 52, 479 (1988)] has been analyzed in terms of recently developed Langevin stochastic model by considering the effects of the density of states and the potential smoothing of the Tsallis transformation. Our study reveals that the fixed points, which are determined by the crossing points of the statistical temperature and the Tsallis effective temperature, play a critical role in overall dynamics of the Tsallis statistics sampling. The dynamical origin of enhanced conformational searches of the Tsallis weight has been investigated by unveiling the intimate relationship between the sampling dynamics and the stability change of corresponding fixed points. Based on this stochastic analysis, we propose one effective method to realize a broad energy distribution in the Tsallis statistics sampling by determining optimal Tsallis parameters systematically based on preliminary canonical samplings. The effectiveness of our method has been validated in the folding simulation of Met-Enkephalin and liquid-solid transition simulation of Lennard-Jones cluster systems.

Determination of multicanonical weight based on a stochastic model of sampling dynamics

Based on the stochastic interpretation of the sampling process modeled by a Langevin equation, we present an effective iteration scheme to determine the weight in multicanonical molecular dynamics. Our method enables an automatic determination of the weight producing a uniform energy sampling via an iterative cancellation of the deterministic force in a Langevin equation. The deterministic force has been calculated from the energy trajectory by identifying the moments of the transition probability of a Fokker-Planck equation associated with a Langevin equation. The intimate relationship between the sampling process and the stochastic dynamics has been verified by applying the iteration scheme to a helix-coil transition of the 8-polyalanine system in a gas phase.

Dynamical origin of uniform sampling in multicanonical ensemble

The stochastic model describing the sampling process in multicanonical ensemble has been derived by considering the sampling process as an overdamped Brownian motion on the free energy surface. The essential dynamics of the multicanonical sampling has been characterized by a Langevin equation in a piecewise multivalleyed free energy landscape, modulated by a temperature-dependent curvature. Based on the stochastic model we showed that the multicanonical weight can be determined by interpolating maximum probability energy points of the canonical samplings at different temperatures.

Adiabatic invariance with first integrals of motion

The construction of a microthermodynamic formalism for isolated systems based on the concept of adiabatic invariance is an old but seldom appreciated effort in the literature, dating back at least to P. Hertz [Ann. Phys. (Leipzig) 33, 225 (1910)]. An apparently independent extension of such formalism for systems bearing additional first integrals of motion was recently proposed by Hans H. Rugh [Phys. Rev. E 64, 055101 (2001)], establishing the concept of adiabatic invariance even in such singular cases. After some remarks in connection with the formalism pioneered by Hertz, it will be suggested that such an extension can incidentally explain the success of a dynamical method for computing the entropy of classical interacting fluids, at least in some potential applications where the presence of additional first integrals cannot be ignored.