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

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

Collision-induced conformational changes in glycine

We present quantum dynamical calculations on the conformational changes of glycine in collisions with the He, Ne, and Ar rare-gas atoms. For two conformer interconversion processes (III-->I and IV-->I), we find that the probability of interconversion is dependent on several factors, including the energy of the collision, the angle at which the colliding atom approaches the glycine molecule, and the strength of the glycine-atom interaction. Furthermore, we show that attractive interactions between the colliding atom and the glycine molecule catalyze conformer interconversion at low collision energies. In previous infrared spectroscopy studies of glycine trapped in rare-gas matrices and helium clusters, conformer III has been consistently observed, but conformer IV has yet to be conclusively detected. Because of the calculated thermodynamic stability of conformer IV, its elusiveness has been attributed to the IV-->I conformer interconversion process. However, our calculations present little indication that IV-->I interconversion occurs more readily than III-->I interconversion. Although we cannot determine whether conformer IV interconverts during experimental Ne- and Ar-matrix depositions, our evidence suggests that the conformer should be present in helium droplets. Anharmonic vibrational frequency calculations illustrate that previous efforts to detect conformer IV may have been hindered by the overlap of its IR-absorption bands with those of other conformers. We propose that the redshifted symmetric –CH2 stretch of conformer IV provides a means for its conclusive experimental detection

Case Reports : Rheumatoid pleural effusions and trapped lung : An uncommon complication of rheumatoid arthritis

Rheumatoid arthritis (RA) is a disease affecting approximately 1[percent] of the population. It is familiarly defined as "chronic, symmetric, debilitating and destructive inflammatory polyarthritis characterized by proliferative synovial tissue (pannus) formation in affected joints" (1). Pain symptoms are typically worse in the morning, thedisease affects females more than males, and ulnar deviation and swan neck deformities are common. Perhaps less well known is the extensive list of extra-articular manifestations (EAMs) that can occur at any time duringthe course of the disease. EAMs affect an estimated 18-41[percent] of patients with RA (1). Renal, pulmonary cardiovascular, nervous, and integumentary system manifestations have all been described. Our case is that of anRA patient with shortness of breath and pleural effusions

Quantum dynamics in the smallest water droplet

Water plays a central role in scientific disciplines ranging from geology to astronomy to biology. Yet it is an extraordinarily dif cult liquid to understand because of its complex, ever-changing patterns of hydrogen bonds. Studies of small water clusters have provided important insights into the concerted hydrogen-bond motions that can occur in water. These studies are also crucial for developing an accurate potential function for simulating the properties of liquid water and ice (1). On page 1310 of this issue, Richardson et al. (2) provide evidence for a concerted type of motion in which two hydrogen bonds in a water cluster are broken simultaneously (see the figure). The results have implications for many areas of scientific study, including the chemistry of polar solvents, the conformations of proteins, and the dissolution of ions in minerals