Accurate Predictions of Water Cluster Formation, (H2O)n=2−10

An efficient mixed molecular dynamics/quantum mechanics model has been applied to the water cluster system. The use of the MP2 method and correlation consistent basis sets, with appropriate correction for BSSE, allows for the accurate calculation of electronic and free energies for the formation of clusters of 2−10 water molecules. This approach reveals new low energy conformers for (H2O)n=7,9,10. The water heptamer conformers comprise five different structural motifs ranging from a three-dimensional prism to a quasi-planar book structure. A prism-like structure is favored energetically at low temperatures, but a chair-like structure is the global Gibbs free energy minimum past 200 K. The water nonamers exhibit less complexity with all the low energy structures shaped like a prism. The decamer has 30 conformers that are within 2 kcal/mol of the Gibbs free energy minimum structure at 298 K. These structures are categorized into four conformer classes, and a pentagonal prism is the most stable structure from 0 to 320 K. Results can be used as benchmark values for empirical water models and density functionals, and the method can be applied to larger water clusters.

Microscopic Picture of Aging in SiO2

We investigate the aging dynamics of amorphous SiO2 via molecular dynamics simulations of a quench from a high temperature Ti to a lower temperature Tf.We obtain a microscopic picture of aging dynamics by analyzing single particle trajectories, identifying jump events when a particle escapes the cage formed by its neighbors, and determining how these jumps depend on the waiting time tw, the time elapsed since the temperature quench to Tf. We find that the only tw-dependent microscopic quantity is the number of jumping particles per unit time, which decreases with age. Similar to previous studies for fragile glass formers, we show here for the strong glass former SiO2 that neither the distribution of jump lengths nor the distribution of times spent in the cage are tw dependent.We conclude that the microscopic aging dynamics is surprisingly similar for fragile and strong glass formers.

Quantum Mechanical Study of Sulfuric Acid Hydration: Atmospheric Implications

The role of the binary nucleation of sulfuric acid in aerosol formation and its implications for global warming is one of the fundamental unsettled questions in atmospheric chemistry. We have investigated the thermodynamics of sulfuric acid hydration using ab initio quantum mechanical methods. For H2SO4(H2O)n where n = 1–6, we used a scheme combining molecular dynamics configurational sampling with high-level ab initio calculations to locate the global and many low lying local minima for each cluster size. For each isomer, we extrapolated the Møller–Plesset perturbation theory (MP2) energies to their complete basis set (CBS) limit and added finite temperature corrections within the rigid-rotor-harmonic-oscillator (RRHO) model using scaled harmonic vibrational frequencies. We found that ionic pair (HSO4–·H3O+)(H2O)n−1clusters are competitive with the neutral (H2SO4)(H2O)n clusters for n ≥ 3 and are more stable than neutral clusters for n ≥ 4 depending on the temperature. The Boltzmann averaged Gibbs free energies for the formation of H2SO4(H2O)n clusters are favorable in colder regions of the troposphere (T = 216.65–273.15 K) for n = 1–6, but the formation of clusters with n ≥ 5 is not favorable at higher (T > 273.15 K) temperatures. Our results suggest the critical cluster of a binary H2SO4–H2O system must contain more than one H2SO4 and are in concert with recent findings(1) that the role of binary nucleation is small at ambient conditions, but significant at colder regions of the troposphere. Overall, the results support the idea that binary nucleation of sulfuric acid and water cannot account for nucleation of sulfuric acid in the lower troposphere.

Structural Analysis of alpha-Fetoprotein (AFP)-like Peptides with Anti-Breast-Cancer Properties

The abundance of alpha-fetoprotein (AFP), a natural protein produced by the fetal yolk sac during pregnancy, correlates with lower incidence of estrogen receptor positive (ER+) breast cancer. The pharmacophore region of AFP has been narrowed down to a four amino acid (AA) region in the third domain of the 591 AA peptide. Our computational study focuses on a 4-mer segment consisting of the amino acids threonine-proline-valine-asparagine (TPVN). We have run replica exchange molecular dynamics (REMD) simulations and used 120 configurational snapshots from the total trajectory as starting configurations for quantum chemical calculations. We optimized structures using semiempirical (PM3, PM6, PM6-D2, PM6-H2, PM6-DH+, PM6-DH2) and density functional methods (TPSS, PBE0, M06-2X). By comparing the accuracy of these methods against RI-MP2 benchmarks, we devised a protocol for calculating the lowest energy conformers of these peptides accurately and efficiently. This protocol screens out high-energy conformers using lower levels of theory and outlines a general method for predicting small peptide structures.

Hydration of the Sulfuric Acid−Methylamine Complex and Implications for Aerosol Formation

The binary H2SO4−H2O nucleation is one of the most important pathways by which aerosols form in the atmosphere, and the presence of ternary species like amines increases aerosol formation rates. In this study, we focus on the hydration of a ternary system of sulfuric acid (H2SO4), methylamine (NH2CH3), and up to six waters to evaluate its implications for aerosol formation. By combining molecular dynamics (MD) sampling with high-level ab initio calculations, we determine the thermodynamics of forming H2SO4(NH2CH3)(H2O)n, where n = 0−6. Because it is a strong acid−base system, H2SO4−NH2CH3 quickly forms a tightly bound HSO4−−NH3CH3+ complex that condenses water more readily than H2SO4 alone. The electronic binding energy of H2SO4−NH2CH3 is −21.8 kcal mol−1 compared with −16.8 kcal mol−1 for H2SO4−NH3 and −12.8 kcal mol−1 for H2SO4−H2O. Adding one to two water molecules to the H2SO4−NH2CH3 complex is more favorable than adding to H2SO4 alone, yet there is no systematic difference for n ≥ 3. However, the average number of water molecules around H2SO4−NH2CH3 is consistently higher than that of H2SO4, and it is fairly independent of temperature and relative humidity.

Hydration of the Sulfuric Acid-Methylamine Complex and Implications for Aerosol Formation

The binary H2SO4-H2O nucleation is one of the most important pathways by which aerosols form in the atmosphere, and the presence of ternary species like amines increases aerosol formation rates. In this study, we focus on the hydration of a ternary system of sulfuric acid (H2SO4), methylamine (NH2CH3), and up to six waters to evaluate its implications for aerosol formation. By combining molecular dynamics (MD) sampling with high-level ab initio calculations, we determine the thermodynamics of forming H2SO4(NH2CH3)(H2O)n, where n = 0-6. Because it is a strong acid-base system, H2SO4-NH2CH3 quickly forms a tightly bound HSO4(-)-NH3CH3(+) complex that condenses water more readily than H2SO4 alone. The electronic binding energy of H2SO4-NH2CH3 is -21.8 kcal mol(-1) compared with -16.8 kcal mol(-1) for H2SO4-NH3 and -12.8 kcal mol(-1) for H2SO4-H2O. Adding one to two water molecules to the H2SO4-NH2CH3 complex is more favorable than adding to H2SO4 alone, yet there is no systematic difference for n ≥ 3. However, the average number of water molecules around H2SO4-NH2CH3 is consistently higher than that of H2SO4, and it is fairly independent of temperature and relative humidity.