Search

In the present paper we identify a rigorous property of a number of tempering-based Monte Carlo sampling methods, including parallel tempering as well as partial and infinite swapping. Based on this property we develop a variety of performance measures for such rare-event sampling methods that are broadly applicable, informative, and straightforward to implement. We illustrate the use of these performance measures with a series of applications involving the equilibrium properties of simple Lennard-Jones clusters, applications for which the performance levels of partial and infinite swapping approaches are found to be higher than those of conventional parallel tempering., Comment: 18 figures

The cumulant representation of the Fourier path integral method is examined to determine the asymptotic convergence characteristics of the imaginary-time density matrix with respect to the number of path variables $N$ included. It is proved that when the cumulant expansion is truncated at order $p$, the asymptotic convergence rate of the density matrix behaves like $N^{-(2p+1)}$. The complex algebra associated with the proof is simplified by introducing a diagrammatic representation of the contributing terms along with an associated linked-cluster theorem. The cumulant terms at each order are expanded in a series such that the the asymptotic convergence rate is maintained without the need to calculate the full cumulant at order $p$. Using this truncated expansion of each cumulant at order $p$, the numerical cost in developing Fourier path integral expressions having convergence order $N^{-(2p+1)}$ is shown to be approximately linear in the number of required potential energy evaluations making the method promising for actual numerical implementation., Comment: 47 pages, 2 figures, submitted to PRE

We present studies of the potential energy landscape of selected binary Lennard-Jones thirteen 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., Comment: 12 pages, 14 figures, submitted to J. Chem. Phys. Reasons for replacing a paper: figures 2, 3, 7 and 11 did not show up correctly

Detailed studies of the thermodynamic properties of selected binary Lennard-Jones clusters of the type X_{13-n}Y_n (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., Comment: 14 pages, 13 figures. submitted to J. Chem. Phys

In the present work we investigate the adequacy of broken-symmetry unrestricted density functional theory (DFT) for constructing the potential energy curve of nickel dimer and nickel hydride, as a model for larger bare and hydrogenated nickel cluster calculations. We use three hybrid functionals: the popular B3LYP, Becke's newest optimized functional Becke98, and the simple FSLYP functional (50% Hartree-Fock and 50% Slater exchange and LYP gradient-corrected correlation functional) with two basis sets: all-electron (AE) Wachters+f basis set and Stuttgart RSC effective core potential (ECP) and basis set. We find that, overall, the best agreement with experiment, comparable to that of the high-level CASPT2, is obtained with B3LYP/AE, closely followed by Becke98/AE and Becke98/ECP. FSLYP/AE and B3LYP/ECP give slightly worse agreement with experiment, and FSLYP/ECP is the only method among the ones we studied that gives an unaceptably large error, underestimating the dissociation energy of nickel dimer by 28%, and being in the largest disagreement with the experiment and the other theoretical predictions., Comment: 17 pages, 7 tables, 7 figures; submitted to J. Chem. Phys.; Revtex4/LaTeX2e. v2 (8/5/04): New (and better) ECP results, without charge density fitting (which was found to give large errors). Subtracted the relativistic corrections from all experimental values

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_13 cluster demonstrate that a second order central-difference scheme should suffice for most applications., Comment: 11 pages, 4 figures

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., Comment: 15 pages, 1 figure, 3 tables

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., Comment: 9 pages, 9 figures, accepted for publication in J. Chem. Phys

We study the asymptotic convergence of the partial averaging method, a technique used in conjunction with the random series implementation of the Feynman-Kac formula. We prove asymptotic bounds valid for most series representations in the case when the potential has first order Sobolev derivatives. If the potential has also second order Sobolev derivatives, we prove a sharper theorem which gives the exact asymptotic behavior of the density matrices. The results are then specialized for the Wiener-Fourier series representation. It is found that the asymptotic behavior is o(1/n^2) if the potential has first order Sobolev derivatives. If the potential has second order Sobolev derivatives, the convergence is shown to be O(1/n^3) and we give the exact expressions for the convergence constants of the density matrices., Comment: 21 pages, 1 figure, submitted to J. Chem. Phys

The heat capacity and isomer distributions of the 38 atom Lennard-Jones cluster have been calculated in the canonical ensemble using parallel tempering Monte Carlo methods. A distinct region of temperature is identified that corresponds to equilibrium between the global minimum structure and the icosahedral basin of structures. This region of temperatures occurs below the melting peak of the heat capacity and is accompanied by a peak in the derivative of the heat capacity with temperature. Parallel tempering is shown to introduce correlations between results at different temperatures. A discussion is given that compares parallel tempering with other related approaches that ensure ergodic simulations., Comment: Journal of Chemical Physics, accepted

We study the 38-atom Lennard-Jones cluster with parallel tempering Monte Carlo methods in the microcanonical and molecular dynamics ensembles. A new Monte Carlo algorithm is presented that samples rigorously the molecular dynamics ensemble for a system at constant total energy, linear and angular momenta. By combining the parallel tempering technique with molecular dynamics methods, we develop a hybrid method to overcome quasi-ergodicity and to extract both equilibrium and dynamical properties from Monte Carlo and molecular dynamics simulations. Several thermodynamic, structural and dynamical properties are investigated for LJ$_{38}$, including the caloric curve, the diffusion constant and the largest Lyapunov exponent. The importance of insuring ergodicity in molecular dynamics simulations is illustrated by comparing the results of ergodic simulations with earlier molecular dynamics simulations., Comment: Journal of Chemical Physics, accepted

The thermodynamic and structural properties of (NH$_4$Cl)$_n$ clusters, n=3-10 are studied. Using the method of simulated annealing, the geometries of several isomers for each cluster size are examined. Jump-walking Monte Carlo simulations are then used to compute the constant-volume heat capacity for each cluster size over a wide temperature range. To carry out these simulations a new parallel algorithm is developed using the Parallel Virtual Machine (PVM) software package. Features of the cluster potential energy surfaces, such as energy differences among isomers and rotational barriers of the ammonium ions, are found to play important roles in determining the shape of the heat capacity curves., Comment: Journal of Chemical Physics, accepted for publication

Monte Carlo methods are used to study the phase transition in ammonium chloride from the orientationally ordered $\delta$ phase to the orientationally disordered $\gamma$ phase. An effective pair potential is used to model the interaction between ions. Thermodynamic properties are computed in the canonical and isothermal-isobaric ensembles. Each ammonium ion is treated as a rigidly rotating body and the lattice is fixed in the low-temperature CsCl geometry. A simple extension of the Metropolis Monte Carlo method is used to overcome quasiergodicity in the rotational sampling. In the constant-$NVT$ calculations the lattice is held rigid; in the constant-$NpT$ calculations the lattice parameter is allowed to fluctuate. In both ensembles the order parameter rapidly falls to zero in the range (200 - 250)K, suggesting that the model disorders at a temperature in fair agreement with the experimental disordering temperature (243K). Peaks in the heat capacity and thermal expansivity curves are also found in the same temperature range., Comment: submitted to Journal of Chemical Physics

Using a Hubbard-Stratonovich transformation coupled with Fourier path integral methods, expressions are derived for the numerical evaluation of the microcanonical density of states for quantum particles obeying Boltzmann statistics. A numerical algorithmis suggested to evaluate the quantum density of states and illustrated on a one-dimensional model system., Comment: Journal of Chemical Physics