Computational Study of the Structure and Thermodynamic Properties of Ammonium Chloride Clusters Using a Parallel J-Walking Approach

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 publicatio

Fourier Path Integral Monte Carlo Method for the Calculation of the Microcanonical Density of States

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 Physic

A computer model of solar panel-plasma interactions

High power solar arrays for satellite power systems are presently being planned with dimensions of kilometers, and with tens of kilovolts distributed over their surface. Such systems face many plasma interaction problems, such as power leakage to the plasma, particle focusing, and anomalous arcing. These effects cannot be adequately modeled without detailed knowledge of the plasma sheath structure and space charge effects. Laboratory studies of 1 by 10 meter solar array in a simulated low Earth orbit plasma are discussed. The plasma screening process is discussed, program theory is outlined, and a series of calibration models is presented. These models are designed to demonstrate that PANEL is capable of accurate self consistant space charge calculations. Such models include PANEL predictions for the Child-Langmuir diode problem

Locating transition states using double-ended classical trajectories

In this paper we present a method for locating transition states and higher-order saddles on potential energy surfaces using double-ended classical trajectories. We then apply this method to 7- and 8-atom Lennard-Jones clusters, finding one previously unreported transition state for the 7-atom cluster and two for the 8-atom cluster.Comment: Journal of Chemical Physics, 13 page

Taming the rugged landscape: production, reordering, and stabilization of selected cluster inherent structures in the X_(13-n)Y_n system

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 correctl

Fault-Tolerant Quantum Dynamical Decoupling

Dynamical decoupling pulse sequences have been used to extend coherence times in quantum systems ever since the discovery of the spin-echo effect. Here we introduce a method of recursively concatenated dynamical decoupling pulses, designed to overcome both decoherence and operational errors. This is important for coherent control of quantum systems such as quantum computers. For bounded-strength, non-Markovian environments, such as for the spin-bath that arises in electron- and nuclear-spin based solid-state quantum computer proposals, we show that it is strictly advantageous to use concatenated, as opposed to standard periodic dynamical decoupling pulse sequences. Namely, the concatenated scheme is both fault-tolerant and super-polynomially more efficient, at equal cost. We derive a condition on the pulse noise level below which concatenated is guaranteed to reduce decoherence.Comment: 5 pages, 4 color eps figures. v3: Minor changes. To appear in Phys. Rev. Let

Energy estimators for random series path-integral methods

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 table

Cumulant Methods and Short Time Propagators

The present paper clarifies a number of issues concerning the general problem of constructing improved short time quantum mechanical propagators. Cumulant methods are shown to be a particularly convenient tool for this task. Numerical results comparing methods based on partial averaging and on gradient approaches are presented for simple model problems and for many particle quantum fluids

Heat capacity estimators for random series path-integral methods by finite-difference schemes

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 figure

Phase changes in 38 atom Lennard-Jones clusters. II: A parallel tempering study of equilibrium and dynamic properties in the molecular dynamics and microcanonical

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, accepte