Quantum diffusion in liquid water from ring polymer molecular dynamics

We have used the ring polymer molecular-dynamics method to study the translational and orientational motions in an extended simple point charge model of liquid water under ambient conditions. We find, in agreement with previous studies, that quantum-mechanical effects increase the self-diffusion coefficient D and decrease the relaxation times around the principal axes of the water molecule by a factor of around 1.5. These results are consistent with a simple Stokes-Einstein picture of the molecular motion and suggest that the main effect of the quantum fluctuations is to decrease the viscosity of the liquid by about a third. We then go on to consider the system-size scaling of the calculated self-diffusion coefficient and show that an appropriate extrapolation to the limit of infinite system size increases D by a further factor of around 1.3 over the value obtained from a simulation of a system containing 216 water molecules. These findings are discussed in light of the widespread use of classical molecular-dynamics simulations of this sort of size to model the dynamics of aqueous systems

Torsional path integral Monte Carlo method for calculating the absolute quantum free energy of large molecules

A new technique for evaluating the absolute free energy of large molecules is presented. Quantum-mechanical contributions to the intramolecular torsions are included via the torsional path integral Monte Carlo (TPIMC) technique. Importance sampling schemes based on uncoupled free rotors and harmonic oscillators facilitate the use of the TPIMC technique for the direct evaluation of quantum partition functions. Absolute free energies are calculated for the molecules ethane, n-butane, n-octane, and enkephalin, and quantum contributions are found to be significant. Comparison of the TPIMC technique with the harmonic oscillator approximation and a variational technique is performed for the ethane molecule. For all molecules, the quantum contributions to free energy are found to be significant but slightly smaller than the quantum contributions to internal energy

Quantum free energies of the conformers of glycine on an ab initio potential energy surface

The torsional path integral Monte Carlo (TPIMC) technique is used to study the five lowest-energy conformers of glycine. The theoretical method provides an anharmonic and quantum-mechanical description of conformational free energy and is used for the first time with an ab initio potential energy surface. The 3-dimensional torsional potential energy surface of glycine was obtained at the MP2/6-311++G** level of theory and is optimized with respect to the non-torsional degrees of freedom. Calculated conformer populations compare well with those reported in recent matrix-isolation infrared spectroscopy experiments. An additional conformer, not yet observed, is predicted to be heavily populated in the thermal equilibria probed by experiment, and a new explanation for its elusiveness is provided. Quantum effects, such as zero point energy, are found to substantially alter conformer populations, and an algorithm for estimating the role of non-torsional vibrations in the conformational thermodynamics of a molecule is introduced

Quantum diffusion in liquid para-hydrogen from ring-polymer molecular dynamics

We have used the ring-polymer molecular dynamics method to calculate approximate Kubo-transformed velocity autocorrelation functions and self-diffusion coefficients for low-pressure liquid para-hydrogen at temperatures of 25 and 14 K. The resulting diffusion coefficients are shown to be consistent with experimental shear viscosities and the established finite-size relation D(L)~=D([infinity])–2.837kBT/6pietaL, where kB is the Boltzmann constant, T the absolute temperature, eta the shear viscosity, and L the length of the (cubic) simulation cell. The diffusion coefficients D(L) obtained in simulations with finite system sizes are therefore too small. However, the extrapolation to infinite system size corrects this deficiency and leads to excellent agreement with experimental results. This both demonstrates the influence of system-size effects on quantum mechanical diffusion coefficients and provides further evidence that ring-polymer molecular dynamics is an accurate as well as practical way of including quantum effects in condensed phase molecular dynamics

Torsional path integral Monte Carlo method for the quantum simulation of large molecules

A molecular application is introduced for calculating quantum statistical mechanical expectation values of large molecules at nonzero temperatures. The Torsional Path Integral Monte Carlo (TPIMC) technique applies an uncoupled winding number formalism to the torsional degrees of freedom in molecular systems. The internal energy of the molecules ethane, n-butane, n-octane, and enkephalin are calculated at standard temperature using the TPIMC technique and compared to the expectation values obtained using the harmonic oscillator approximation and a variational technique. All studied molecules exhibited significant quantum mechanical contributions to their internal energy expectation values according to the TPIMC technique. The harmonic oscillator approximation approach to calculating the internal energy performs well for the molecules presented in this study but is limited by its neglect of both anharmonicity effects and the potential coupling of intramolecular torsion

Torsional anharmonicity in the conformational thermodynamics of flexible molecules

We present an algorithm for calculating the conformational thermodynamics of large, flexible molecules that combines ab initio electronic structure theory calculations with a torsional path integral Monte Carlo (TPIMC) simulation. The new algorithm overcomes the previous limitations of the TPIMC method by including the thermodynamic contributions of non-torsional vibrational modes and by affordably incorporating the ab initio calculation of conformer electronic energies, and it improves the conventional ab initio treatment of conformational thermodynamics by accounting for the anharmonicity of the torsional modes. Using previously published ab initio results and new TPIMC calculations, we apply the algorithm to the conformers of the adrenaline molecule

MIT Space Engineering Research Center

The Space Engineering Research Center (SERC) at MIT, started in Jul. 1988, has completed two years of research. The Center is approaching the operational phase of its first testbed, is midway through the construction of a second testbed, and is in the design phase of a third. We presently have seven participating faculty, four participating staff members, ten graduate students, and numerous undergraduates. This report reviews the testbed programs, individual graduate research, other SERC activities not funded by the Center, interaction with non-MIT organizations, and SERC milestones. Published papers made possible by SERC funding are included at the end of the report

On the Areas of Cyclic and Semicyclic Polygons

We investigate the ``generalized Heron polynomial'' that relates the squared area of an n-gon inscribed in a circle to the squares of its side lengths. For a (2m+1)-gon or (2m+2)-gon, we express it as the defining polynomial of a certain variety derived from the variety of binary (2m-1)-forms having m-1 double roots. Thus we obtain explicit formulas for the areas of cyclic heptagons and octagons, and illuminate some mysterious features of Robbins' formulas for the areas of cyclic pentagons and hexagons. We also introduce a companion family of polynomials that relate the squared area of an n-gon inscribed in a circle, one of whose sides is a diameter, to the squared lengths of the other sides. By similar algebraic techniques we obtain explicit formulas for these polynomials for all n <= 7.Comment: 22 page

What Effect has Population Growth on People?

Author Institution: Department of Zoology and Entomology, The Ohio State University, Columbus 1

Middeck Active Control Experiment (MACE), phase A

A rationale to determine which structural experiments are sufficient to verify the design of structures employing Controlled Structures Technology was derived. A survey of proposed NASA missions was undertaken to identify candidate test articles for use in the Middeck Active Control Experiment (MACE). The survey revealed that potential test articles could be classified into one of three roles: development, demonstration, and qualification, depending on the maturity of the technology and the mission the structure must fulfill. A set of criteria was derived that allowed determination of which role a potential test article must fulfill. A review of the capabilities and limitations of the STS middeck was conducted. A reference design for the MACE test article was presented. Computing requirements for running typical closed-loop controllers was determined, and various computer configurations were studied. The various components required to manufacture the structure were identified. A management plan was established for the remainder of the program experiment development, flight and ground systems development, and integration to the carrier. Procedures for configuration control, fiscal control, and safety, reliabilty, and quality assurance were developed