Search

The present paper explores a simple approach to the question of parallel tempering temperature selection. We argue that to optimize the performance of parallel tempering it is reasonable to require that the increase in entropy between successive temperatures be uniform over the entire ensemble. An estimate of the system’s heat capacity, obtained either from experiment, a preliminary simulation, or a suitable physical model, thus provides a means for generating the desired tempering ensemble. Applications to the two-dimensional Ising problem indicate that the resulting method is effective, simple to implement, and robust with respect to its sensitivity to the quality of the underlying heat capacity model. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

We present studies of the potential energy landscape of selected binary Lennard-Jones 13 atom clusters. The effect of adding selected impurity atoms to a homogeneous cluster is explored. We analyze the energy landscapes of the studied systems using disconnectivity graphs. The required inherent structures and transition states for the construction of disconnectivity graphs are found by combination of conjugate gradient and eigenvector-following methods. We show that it is possible to controllably induce new structures as well as reorder and stabilize existing structures that are characteristic of higher-lying minima. Moreover, it is shown that the selected structures can have experimentally relevant lifetimes. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

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. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

The structures and energetic effects of molecular nitrogen adsorbates on nickel clusters are investigated using an extended Hückel model coupled with two models of the adsorbate–nickel interaction. The potential parameters for the adsorbates are chosen to mimic experimental information about the binding strength of nitrogen on both cluster and bulk surface phases of nickel. The first model potential is a simple Lennard-Jones interaction that leads to binding sites in holes defined by sets of near-neighbor nickel atoms. The second model potential has a simple three-body form that forces the model nitrogen adsorbates to bind directly to single nickel atoms. Significant rearrangement of the core nickel structures are found in both models. A disconnectivity graph analysis of the potential energy surfaces implies that the rearrangements arise from low transition state barriers and the small differences between available isomers in the nickel core. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

Previous heat capacity estimators used in path integral simulations either have large variances that grow to infinity with the number of path variables or require the evaluation of first- and second-order derivatives of the potential. In the present paper, we show that the evaluation of the total energy by the T-method estimator and of the heat capacity by the TT-method estimator can be implemented by a finite difference scheme in a stable fashion. As such, the variances of the resulting estimators are finite and the evaluation of the estimators requires the potential function only. By comparison with the task of computing the partition function, the evaluation of the estimators requires k+1 times more calls to the potential, where k is the order of the difference scheme employed. Quantum Monte Carlo simulations for the Ne[sub 13] cluster demonstrate that a second order central-difference scheme should suffice for most applications. © 2003 American Institute of Physics. [ABSTRACT FROM AUTHOR]

We perform a thorough analysis on the choice of estimators for random series path integral methods. In particular, we show that both the thermodynamic (T-method) and the direct (H-method) energy estimators have finite variances and are straightforward to implement. It is demonstrated that the agreement between the T-method and the H-method estimators provides an important consistency check on the quality of the path integral simulations. We illustrate the behavior of the various estimators by computing the total, kinetic, and potential energies of a molecular hydrogen cluster using three different path integral techniques. Statistical tests are employed to validate the sampling strategy adopted as well as to measure the performance of the parallel random number generator utilized in the Monte Carlo simulation. Some issues raised by previous simulations of the hydrogen cluster are clarified. © 2003 American Institute of Physics. [ABSTRACT FROM AUTHOR]

The reweighted random series techniques provide finite-dimensional approximations to the quantum density matrix of a physical system that have fast asymptotic convergence. We study two special reweighted techniques that are based upon the Lévy–Ciesielski and Wiener–Fourier series, respectively. In agreement with the theoretical predictions, we demonstrate by numerical examples that the asymptotic convergence of the two reweighted methods is cubic for smooth enough potentials. For each reweighted technique, we propose some minimalist quadrature techniques for the computation of the path averages. These quadrature techniques are designed to preserve the asymptotic convergence of the original methods. © 2003 American Institute of Physics. [ABSTRACT FROM AUTHOR]

We report our studies of the potential energy surface (PES) of selected binary Lennard-Jones clusters. The effect of adding selected impurity atoms to a homogeneous cluster is explored. Inherent structures and transition states are found by combination of conjugate gradient and eigenvector-following methods while the topography of the PES is mapped with the help of a disconnectivity analysis. We show that we can controllably induce new structures as well as reorder and stabilize existing structures that are characteristic of higher-lying minima. [ABSTRACT FROM AUTHOR]

In the present paper we describe a stochastic quadrature method is designed for the evaluation of generalized, complex averages. Motivated by recent advances in sparse sampling techniques, this method is based on a combination of parallel tempering and stationary phase filtering methods. Numerical applications of the resulting "stationary tempering" approach are presented. We also examine inherent structure decomposition from a probabilistic clustering perspective. [ABSTRACT FROM AUTHOR]

Presents calculations of the torsional state dependence of A = B and C rotational constants of (D[sub 2]O)[sub 3], for torsional energies of up to 100 cm [sup -1] above the ground state. Changes of rms torsional angles of the 'free' O-D bonds; Changes in the oxygen-oxygen separation which accompany torsional excitation.

Reports that a four-dimensional intermolecular potential energy surfaces for the water trimer was calculated in the subspace of the torsional coordinates and the symmetric intermolecular stretch coordinate. Quantum three-dimensional calculation of torsional levels of water trimers; Torsional averaging of the rotational constants of water trimer.

Presents the quantum four-dimensional (4D) model calculations of coupled intermolecular torsional vibrations of cyclic homodromic water tetrameres. Context of a 4D discrete variable representation approach; Calculation of the lowest 50 torsional levels; Estimation of higher excited pseudorotational states.

Studies the intermolecular vibrations and van der Waals isomerism of the 2,3-dimethylnaphthalene...He van der Waals complex. Potential energy surface; Quantum calculations of the van der Waals vibrational levels; Reproduction of observed bands as intermolecular vibration excitations; Delocalization of van der Waals over the potential surface.

Studies effects of intermolecular torsional and symmetric stretching exciations of rotational constants water trimer. Decrease in the rotational constant with the increase of torsional energy; Calculated rotational constant and the root-mean-square torsional angles of water.

In the present paper we explore the use of generalized Gaussian quadrature methods in the context of equilibrium path integral applications. Using moment techniques, we devise a compact, self-adaptive approach for use in conjunction with selected classes of interaction potentials. We demonstrate that, when applicable, the resulting approach reduces appreciably the number of potential energy evaluations required in equilibrium path integral simulations. © 2000 American Institute of Physics. [ABSTRACT FROM AUTHOR]