The Role of Anharmonicity in Hydrogen-Bonded Systems: The Case of Water Clusters

The nature of vibrational anharmonicity has been examined for the case of small water clusters using second-order vibrational perturbation theory (VPT2) applied on second-order Møller–Plesset perturbation theory (MP2) potential energy surfaces. Using a training set of 16 water clusters (H2O)n=2–6,8,9 with a total of 723 vibrational modes, we determined scaling factors that map the harmonic frequencies onto anharmonic ones. The intermolecular modes were found to be substantially more anharmonic than intramolecular bending and stretching modes. Due to the varying levels of anharmonicity of the intermolecular and intramolecular modes, different frequency scaling factors for each region were necessary to achieve the highest accuracy. Furthermore, new scaling factors for zero-point vibrational energies (ZPVE) and vibrational corrections to the enthalpy (ΔHvib) and the entropy (Svib) have been determined. All the scaling factors reported in this study are different from previous works in that they are intended for hydrogen-bonded systems, while others were built using experimental frequencies of covalently bonded systems. An application of our scaling factors to the vibrational frequencies of water dimer and thermodynamic functions of 11 larger water clusters highlights the importance of anharmonic effects in hydrogen-bonded systems.

Thermodynamics of Forming Water Clusters at Various Temperatures and Pressures by Gaussian-2, Gaussian-3, Complete Basis Set-QB3, and Complete Basis Set-APNO Model Chemistries; Implications for Atmospheric Chemistry

The Gaussian-2, Gaussian-3, complete basis set- (CBS-) QB3, and CBS-APNO methods have been used to calculate ΔH° and ΔG° values for neutral clusters of water, (H2O)n, where n = 2−6. The structures are similar to those determined from experiment and from previous high-level calculations. The thermodynamic calculations by the G2, G3, and CBS-APNO methods compare well against the estimated MP2(CBS) limit. The cyclic pentamer and hexamer structures release the most heat per hydrogen bond formed of any of the clusters. While the cage and prism forms of the hexamer are the lowest energy structures at very low temperatures, as temperature is increased the cyclic structure is favored. The free energies of cluster formation at different temperatures reveal interesting insights, the most striking being that the cyclic trimer, cyclic tetramer, and cyclic pentamer, like the dimer, should be detectable in the lower troposphere. We predict water dimer concentrations of 9 × 1014 molecules/cm3, water trimer concentrations of 2.6 × 1012 molecules/cm3, tetramer concentrations of approximately 5.8 × 1011 molecules/cm3, and pentamer concentrations of approximately 3.5 × 1010 molecules/cm3 in saturated air at 298 K. These results have important implications for understanding the gas-phase chemistry of the lower troposphere.

Prediction of Accurate Anharmonic Experimental Vibrational Frequencies for Water Clusters, (H2O)n, n = 2−5

Accurate anharmonic experimental vibrational frequencies for water clusters consisting of 2−5 water molecules have been predicted on the basis of comparing different methods with MP2/aug-cc-pVTZ calculated and experimental anharmonic frequencies. The combination of using HF/6-31G* scaled frequencies for intramolecular modes and anharmonic frequencies for intermolecular modes gives excellent agreement with experiment for the water dimer and trimer and are as good as the expensive anharmonic MP2 calculations. The water trimer, the cyclic Ci and S4 tetramers, and the cyclic pentamer all have unique peaks in the infrared spectrum between 500 and 800 cm-1 and between 3400 and 3700 cm-1. Under the right experimental conditions these different clusters can be uniquely identified using high-resolution IR spectroscopy.

In search of CS2(H2O)n = 1–4 clusters

Gaussian-3 and MP2/aug-cc-pVnZ methods have been used to calculate geometries and thermochemistry of CS2(H2O)n, where n = 1–4. An extensive molecular dynamics search followed by optimization using these two methods located two dimers, six trimers, six tetramers, and two pentamers. The MP2/aug-cc-pVDZ structure matched best with the experimental result for the CS2(H2O) dimer, showing that diffuse functions are necessary to model the interactions found in this complex. For larger CS2(H2O)n clusters, the MP2/aug-cc-pVDZ minima are significantly different from the MP2(full)/6-31G* structures, revealing that the G3 model chemistry is not suitable for investigation of sulfur containing van der Waals complexes. Based on the MP2/aug-cc-pVTZ free energies, the concentration of saturated water in the atmosphere and the average amount of CS2 in the atmosphere, the concentrations of these clusters are predicted to be on the order of 105CS2(H2O) clusters∙cm−3 and 102 CS2(H2O)2 clusters∙cm−3 at 298.15 K. The MP2/aug-cc-pVDZ scaled harmonic and anharmonic frequencies of the most abundant dimer cluster at 298 K are presented, along with the MP2/aug-cc-pVDZ scaled harmonic frequencies for the CS2(H2O)n structures predicted to be present in a low-temperature molecular beam experiment.

The Role of Anharmonicity in Hydrogen-Bonded Systems: The Case of Water Clusters

The nature of vibrational anharmonicity has been examined for the case of small water clusters using second-order vibrational perturbation theory (VPT2) applied on second-order Møller–Plesset perturbation theory (MP2) potential energy surfaces. Using a training set of 16 water clusters (H2O)n=2–6,8,9 with a total of 723 vibrational modes, we determined scaling factors that map the harmonic frequencies onto anharmonic ones. The intermolecular modes were found to be substantially more anharmonic than intramolecular bending and stretching modes. Due to the varying levels of anharmonicity of the intermolecular and intramolecular modes, different frequency scaling factors for each region were necessary to achieve the highest accuracy. Furthermore, new scaling factors for zero-point vibrational energies (ZPVE) and vibrational corrections to the enthalpy (ΔHvib) and the entropy (Svib) have been determined. All the scaling factors reported in this study are different from previous works in that they are intended for hydrogen-bonded systems, while others were built using experimental frequencies of covalently bonded systems. An application of our scaling factors to the vibrational frequencies of water dimer and thermodynamic functions of 11 larger water clusters highlights the importance of anharmonic effects in hydrogen-bonded systems.

Computational Study of the Hydration of Sulfuric Acid Dimers: Implications for Acid Dissociation and Aerosol Formation

We have investigated the thermodynamics of sulfuric acid dimer hydration using ab initio quantum mechanical methods. For (H2SO4)2(H2O)n where n = 0−6, we employed high-level ab initio calculations to locate the most stable minima for each cluster size. The results presented herein yield a detailed understanding of the first deprotonation of sulfuric acid as a function of temperature for a system consisting of two sulfuric acid molecules and up to six waters. At 0 K, a cluster of two sulfuric acid molecules and one water remains undissociated. Addition of a second water begins the deprotonation of the first sulfuric acid leading to the di-ionic species (the bisulfate anion HSO4−, the hydronium cation H3O+, an undissociated sulfuric acid molecule, and a water). Upon the addition of a third water molecule, the second sulfuric acid molecule begins to dissociate. For the (H2SO4)2(H2O)3 cluster, the di-ionic cluster is a few kcal mol−1 more stable than the neutral cluster, which is just slightly more stable than the tetra-ionic cluster (two bisulfate anions, two hydronium cations, and one water). With four water molecules, the tetra-ionic cluster, (HSO4−)2(H3O+)2(H2O)2, becomes as favorable as the di-ionic cluster H2SO4(HSO4−)(H3O+)(H2O)3 at 0 K. Increasing the temperature favors the undissociated clusters, and at room temperature we predict that the di-ionic species is slightly more favorable than the neutral cluster once three waters have been added to the cluster. The tetra-ionic species competes with the di-ionic species once five waters have been added to the cluster. The thermodynamics of stepwise hydration of sulfuric acid dimer is similar to that of the monomer; it is favorable up to n = 4−5 at 298 K. A much more thermodynamically favorable pathway forming sulfuric acid dimer hydrates is through the combination of sulfuric acid monomer hydrates, but the low concentration of sulfuric acid relative to water vapor at ambient conditions limits that process.

5a-Butyl-1,3,8,10-tetra-chloro-7,13-bis-(4-nitro-benzo-yl)-5a,6a,12a,12b-tetra-hydro-7H,13H-thieno[2,3-b:4,5-b']bis-(1,4-benzoxazine)

The title compound, C(34)H(24)Cl(4)N(4)O(8)S, is a linear penta-cyclic system formed of two substituted benzoxazinyl groups fused to 2-n-butyl-tetra-hydro-thio-phene. The oxazine ring, which is fused to the n-butyl-substituted side of the thio-phene ring, is in a boat conformation. The other fused oxazine ring and the tetra-hydro-thiene ring are each in an envelope conformation. The bridgehead C atom alpha to both the S and N atoms forms the flap of each envelope. This results in a twist of the penta-cyclic system such that the dihedral angle between the terminal dichloro-benzene rings is 82.92 (8)°. In the crystal, inversion-related mol-ecules form a weakly hydrogen-bonded dimer, with two C-H⋯O inter-actions between an H atom on the oxazine ring and an amide O atom. Additionally, C-H⋯O inter-actions occur between an H atom on a screw-related nitro-benzene ring and an O atom on the nitro-benzene ring of one mol-ecule. One of the Cl atoms and the butyl group are disordered over two sets of sites with occupancy ratios of 0.94 (2):0.06 (2) and 0.624 (4):0.376 (4), respectively.

Purification and Preliminary Characterization of the TM0727 Protein from Thermotoga maritima

The TM0727 gene of Thermotoga maritima is responsible for encoding what has been reported to be a modulator of DNA gyrase (pmbA). Although the function of pmbA is still unknown, it is believedto be involved in cell division, carbon storage regulation, and the synthesis of the antibiotic peptide microcin B17. It is suggested that it serves together with tldD, a known zinc dependent protease, tomodulate DNA gyrase. TM0727 is believed to be a zinc dependent protease that binds zinc in the central active site of the molecule, located between two equivalent monomeric units. However, thecrystal structure determined by Wilson et al. (2005) did not contain zinc. It therefore remains to be seen if TM0727 requires zinc for activity, or regulation, and if the protein is indeed a protease. To begin studying this protein, the gene was expressed in BL21(DE3) pLysS cells and the induction time was optimized. Using affinity and ion exchange chromatography, the protein has been successfully purified. The purification procedure can be replicated to obtain sufficient protein for characterization. Purification results show that the protein loses stability after 24 hours and remains stable under an imidazole-free lysis workup. Preliminary characterization of TM0727 has focused on understanding the protein’s structuralproperties through tryptophan fluorescence anisotropy measurements. The four tryptophan residues located within the TM0727 dimer fluoresce at different maximum wavelengths and with differentintensities upon excitation with 295nm light. These emission properties are highly sensitive to the environment (solvent, surrounding residues) of each tryptophan residue. The low number oftryptophans allows for a specific monitoring of the protein’s structure as it denatures. As more denaturant is added to the protein, its tryptophan environments have clearly altered. This is indicative of unfolding and increased solvent exposure of the protein. This unfolding has been confirmed with the addition of a fluorescent quencher. Additionally, fluorescence anisotropy measurements have been carried out on the protein to gain a preliminary understanding of the rotational dynamics of the tryptophan residues. These experiments excite the tryptophan residues within the sample using a polarized light source. Polarized emission is then detected, the degree of which depends on the rotational dynamics and local environment of the tryptophan residues. The protein was denatured and the changes in emission were recorded to detect these structural changes. Results have shown a large change in quaternary structure, consistent with a dimer to monomer transition, occurs at 1.5M Guandidine HCl. There has also been an examination of the crystal structure for the location of a potential active site. The inner cavity of the protein was inspected visually to locate a potential location for a catalytic triad, specifically the amino acids found in the active sites of serine, cyteine, and aspartateproteases. It was found that a potential aspartic protease active site may be located between the Asparate286 and Aspartate287 residues. Further investigation is warranted to test this remotepossibility.

Synthesis of Mid-Chain Functionalized Polymers

Polymers with mid-chain alkoxyamine functionality were synthesized by activating monohalogenated polymers in the presence of nitroso or nitrone radical traps. The resulting polymers were either polystyrene (PSt) homopolymers with a mid-chain alkoxyamine or PSt-poly(methyl acrylate) (PMA) diblock copolymers with an alkoxyamine unit at the junction between the segments. Monohalogenated polymers where synthesized by atom transfer radical polymerization (ATRP) and were then reacted to form polymer radicals in the presence of a radical trap, nitrone or nitroso. When only polystyrene radicals were reacted with the radical trap a dimer was formed with an alkoxyamine functionality in the center of the polymer chain. This functionality allowed the polymer chain to be cleaved in order to visualize the extent of the alkoxyamine functionality incorporation into the polymer chains. It was found that near quantitative alkoxyamine mid-chain functionality could be achieved by activating the PStBr in the presence of 10 equivalents of nitrone, 5 equivalents of copper bromide, and 2 equivalents of copper metal. Further reducing the amount of copper metal led to incomplete coupling, while increasing the equivalents beyond 2 generated polymer dimers with less than quantitative mid-chain functionality. Monochlorinated polystyrene (PStCl) precursors gave much poorer coupling results compared to reactions with PStBr, which is consistent with the stronger C-Cl bond resisting activation and the formation of the polystyryl radicals. When poly (methyl acrylate) (PMABr) is reacted with PStBr in the presence of a nitroso group at reduced temperatures (30 oC) block copolymers were selectively formed with an alkoxyamine functionality in the center. This was done by first activating the PSt-Br to form a polymer radical that would react with the radical trap to form a persistent radical on the oxygen. The PMA-Br, once activated, reacted with the radical on the oxygen to form the block copolymer. To test the amount of functionality incorporated, a coupling reaction was performed with no nitroso present, and found that no reaction occurred. This showed that the radical trap is essential for the coupling to occur, and cleavage of the diblock indicated that the alkoxyamine functionality was indeed incorporated into the diblock.