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.

Theoretical Calculations of Acid Dissociation Constants: A Review Article

Acid dissociation constants, or pKa values, are essential for understanding many fundamental reactions in chemistry. These values reveal the deprotonation state of a molecule in a particular solvent. There is great interest in using theoretical methods to calculate the pKa values for many different types of molecules. These include molecules that have not been synthesized, those for which experimental pKa determinations are difficult, and for larger molecules where the local environment changes the usual pKa values, such as for certain amino acids that are part of a larger polypeptide chain. Chemical accuracy in pKa calculations is difficult to achieve, because an error of 1.36 kcal/mol in the change of free energy for deprotonation in solvent results in an error of 1 pKa unit. In this review the most valuable methods for determining accurate pKa values in aqueous solution are presented for educators interested in explaining or using these methods for their students.

Benchmark Structures and Binding Energies of Small Water Clusters with Anharmonicity Corrections

For (H2O)n where n = 1–10, we used a scheme combining molecular dynamics sampling with high level ab initio calculations to locate the global and many low lying local minima for each cluster. For each isomer, we extrapolated the RI-MP2 energies to their complete basis set limit, included a CCSD(T) correction using a smaller basis set and added finite temperature corrections within the rigid-rotor-harmonic-oscillator (RRHO) model using scaled and unscaled harmonic vibrational frequencies. The vibrational scaling factors were determined specifically for water clusters by comparing harmonic frequencies with VPT2 fundamental frequencies. We find the CCSD(T) correction to the RI-MP2 binding energy to be small (<1%) but still important in determining accurate conformational energies. Anharmonic corrections are found to be non-negligble; they do not alter the energetic ordering of isomers, but they do lower the free energies of formation of the water clusters by as much as 4 kcal/mol at 298.15 K.

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.

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.

Broadband Fourier transform rotational spectroscopy for structure determination: The water heptamer

Over the recent years chirped-pulse, Fourier-transform microwave (CP-FTMW) spectrometers have chan- ged the scope of rotational spectroscopy. The broad frequency and large dynamic range make possible structural determinations in molecular systems of increasingly larger size from measurements of heavy atom (13C, 15N, 18O) isotopes recorded in natural abundance in the same spectrum as that of the parent isotopic species. The design of a broadband spectrometer operating in the 2–8 GHz frequency range with further improvements in sensitivity is presented. The current CP-FTMW spectrometer performance is benchmarked in the analyses of the rotational spectrum of the water heptamer, (H2O)7, in both 2– 8 GHz and 6–18 GHz frequency ranges. Two isomers of the water heptamer have been observed in a pulsed supersonic molecular expansion. High level ab initio structural searches were performed to pro- vide plausible low-energy candidates which were directly compared with accurate structures provided from broadband rotational spectra. The full substitution structure of the most stable species has been obtained through the analysis of all possible singly-substituted isotopologues (H218O and HDO), and a least-squares rm(1) geometry of the oxygen framework determined from 16 different isotopic species compares with the calculated O–O equilibrium distances at the 0.01 Å level.

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.

The Carrington event not observed in most ice core nitrate records

The Carrington Event of 1859 is considered to be among the largest space weather events of the last 150 years. We show that only one out of 14 well-resolved ice core records from Greenland and Antarctica has a nitrate spike dated to 1859. No sharp spikes are observed in the Antarctic cores studied here. In Greenland numerous spikes are observed in the 40 years surrounding 1859, but where other chemistry was measured, all large spikes have the unequivocal signal, including co-located spikes in ammonium, formate, black carbon and vanillic acid, of biomass burning plumes. It seems certain that most spikes in an earlier core, including that claimed for 1859, are also due to biomass burning plumes, and not to solar energetic particle (SEP) events. We conclude that an event as large as the Carrington Event did not leave an observable, widespread imprint in nitrate in polar ice. Nitrate spikes cannot be used to derive the statistics of SEPs.

EVOLUTION OF SELECTED ISOPRENE OXIDATION PRODUCTS IN DARK AQUEOUS AMMONIUM SULFATE

The aqueous phase processing of glyoxylic acid, pyruvic acid, oxalic acid and methylglyoxal was studied simulating dark and radical free atmospheric aqueous aerosol. A novel observation of the cleavage of a carbon-carbon bond in pyruvic acid and glyoxylic acid leading to their decarboxylation was made in the presence of ammonium salts but no decarboxylation was observed from oxalic acid. The empirical rate constants for decarboxylation were determined. The structure of the acid, ionic environment of solution and concentration of species found to affect the decarboxylation process. A tentative set of reaction mechanisms was proposed involving nucleophilic attack by ammonia on the carbonyl carbon leading to fragmentation of the carbon-carbon bond between the carbonyl and carboxyl carbons. Whereas, the formation of high molecular weight organic species was observed in the case of methylglyoxal. The elemental compositions of the species were determined. It was concluded that, additional pathways that are not currently known likely contribute to aqueous phase processing leading to high molecular weight organic species. Under similar conditions in atmospheric aerosol, the aqueous phase processing will markedly impact the physicochemical properties of aerosol.

Development of microsatellite markers in Cannabis Sativa for fingerprinting and genetic relatedness analyses

Microsatellite markers were developed for Cannabis sativa L. (marijuana) to estimate the level of polymorphism, usefulness for DNA typing (genotype identification), and to measure the genetic relationships between the different plants. Twelve different oligonucleotide probes were used to screen an enriched microsatellite library of Cannabis sativa in which 49% of the clones contained microsatellite sequences. Characterization of microsatellite loci in Cannabis revealed that GA/CT was the most abundant class of isolated microsatellites representing 50% overall. Eleven polymorphic SSR markers were developed, derived from dinucleotide motifs and eight from trinucleotide motifs. A total of 52 alleles were detected averaging 4.7 alleles/locus. The expected heterozygosity of the eleven loci ranged between 0.368 and 0.710 and the common probability of identical genotypes was 1.8 x 107. The loci identified 27 unique profiles of the 41 Cannabis samples. The eleven microsatellite markers developed in this study were found to be useful for DNA fingerprinting and for assessing genetic relationships in Cannabis.

Development of a Surface-Enhanced Raman Spectroscopy Method for the Detection of Benzodiazepines in Urine

Benzodiazepines are among the most prescribed compounds for anti-anxiety and are present in many toxicological screens. These drugs are also prominent in the commission of drug facilitated sexual assaults due their effects on the central nervous system. Due to their potency, a low dose of these compounds is often administered to victims; therefore, the target detection limit for these compounds in biological samples is 10 ng/mL. Currently these compounds are predominantly analyzed using immunoassay techniques; however more specific screening methods are needed. The goal of this dissertation was to develop a rapid, specific screening technique for benzodiazepines in urine samples utilizing surface-enhanced Raman spectroscopy (SERS), which has previously been shown be capable of to detect trace quantities of pharmaceutical compounds in aqueous solutions. Surface enhanced Raman spectroscopy has the advantage of overcoming the low sensitivity and fluorescence effects seen with conventional Raman spectroscopy. The spectra are obtained by applying an analyte onto a SERS-active metal substrate such as colloidal metal particles. SERS signals can be further increased with the addition of aggregate solutions. These agents cause the nanoparticles to amass and form hot-spots which increase the signal intensity. In this work, the colloidal particles are spherical gold nanoparticles in aqueous solution with an average size of approximately 30 nm. The optimum aggregating agent for the detection of benzodiazepines was determined to be 16.7 mM MgCl2, providing the highest signal intensities at the lowest drug concentrations with limits of detection between 0.5 and 127 ng/mL. A supported liquid extraction technique was utilized as a rapid clean extraction for benzodiazepines from urine at a pH of 5.0, allowing for clean extraction with limits of detection between 6 and 640 ng/mL. It was shown that at this pH other drugs that are prevalent in urine samples can be removed providing the selective detection of the benzodiazepine of interest. This technique has been shown to provide rapid (less than twenty minutes), sensitive, and specific detection of benzodiazepines at low concentrations in urine. It provides the forensic community with a sensitive and specific screening technique for the detection of benzodiazepines in drug facilitated assault cases.