Computerized polar plots by a cathode ray tube/grid overlay method

Overlay is aligned with four calibration dots so it is not affected by CRT drift or changes in vertical or horizontal gain when producing Nyquist /frequency-response phase/amplitude/ plots. Method produces over 50 plots per hour

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

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

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

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

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

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

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

Radiolysis of Macromolecular Organic Material in Mars-Relevant Mineral Matrices

The fate of organic material on Mars after deposition is crucial to interpreting the source of these molecules. Previous work has addressed how various organic compounds at millimeter depths in sediments respond to ultraviolet radiation. In contrast, this study addressed how highenergy particle radiation (200MeV protons, simulating the effect of galactic cosmic rays and solar wind at depths of <45 cm) inuences organic macromolecules in sediments. Specically, we report the generation of organicacid radiolysis products after exposure to radiation doses equivalent to geological time scales (17 Myr). We found that formate and oxalate were produced from a variety of organic starting materials and mineral matrices. Unlike ultravioletdriven reactions that can invoke Fenton chemistry to produce organic acids, our work suggests that irradiation of semiconductor surfaces, such as TiO2 or possible clay minerals found on Mars, forms oxygen and hydroxyl radical species, which can break down macromolecules into organic acids. We also investigated the metastability of benzoate in multiple mineral matrices. Benzoate was added to samples prior to irradiation and persisted up to 500 kGys of exposure. Our ndings suggest that organic acids are likely a major component of organic material buried at depth on Mars