Experimental and Theoretical Study of the OH Vibrational Spectra and Overtone Chemistry of Gas-Phase Vinylacetic Acid

In this study we present the gas-phase vibrational spectrum of vinylacetic acid with a focus on the ν = 1−5 vibrational states of the OH stretching transitions. Cross sections for ν = 1, 2, 4 and 5 of the OH stretching vibrational transitions are derived on the basis of the vapor pressure data obtained for vinylacetic acid. Ab initio calculations are used to assist in the band assignments of the experimental spectra, and to determine the threshold for the decarboxylation of vinylacetic acid. When compared to the theoretical energy barrier to decarboxylation, it is found that the νOH = 4 transition with thermal excitation of low frequency modes or rotational motion and νOH = 5 transitions have sufficient energy for the reaction to proceed following overtone excitation.

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−1 clusters 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.

Ability of the PM3 Quantum Mechanical Method to Model Intermolecular Hydrogen Bonding between Neutral Molecules

The PM3 semiempirical quantum-mechanical method was found to systematically describe intermolecular hydrogen bonding in small polar molecules. PM3 shows charge transfer from the donor to acceptor molecules on the order of 0.02-0.06 units of charge when strong hydrogen bonds are formed. The PM3 method is predictive; calculated hydrogen bond energies with an absolute magnitude greater than 2 kcal mol-' suggest that the global minimum is a hydrogen bonded complex; absolute energies less than 2 kcal mol-' imply that other van der Waals complexes are more stable. The geometries of the PM3 hydrogen bonded complexes agree with high-resolution spectroscopic observations, gas electron diffraction data, and high-level ab initio calculations. The main limitations in the PM3 method are the underestimation of hydrogen bond lengths by 0.1-0.2 for some systems and the underestimation of reliable experimental hydrogen bond energies by approximately 1-2 kcal mol-l. The PM3 method predicts that ammonia is a good hydrogen bond acceptor and a poor hydrogen donor when interacting with neutral molecules. Electronegativity differences between F, N, and 0 predict that donor strength follows the order F > 0 > N and acceptor strength follows the order N > 0 > F. In the calculations presented in this article, the PM3 method mirrors these electronegativity differences, predicting the F-H- - -N bond to be the strongest and the N-H- - -F bond the weakest. It appears that the PM3 Hamiltonian is able to model hydrogen bonding because of the reduction of two-center repulsive forces brought about by the parameterization of the Gaussian core-core interactions. The ability of the PM3 method to model intermolecular hydrogen bonding means reasonably accurate quantum-mechanical calculations can be applied to small biologic systems.

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 (

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.

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.

Structure and thermodynamics of H3O+(H2O)8 clusters: A combined molecular dynamics and quantum mechanics approach

We have studied the structure and stability of (H3O+)(H2O)8 clusters using a combination of molecular dynamics sampling and high-level ab initio calculations. 20 distinct oxygen frameworks are found within 2 kcal/mol of the electronic or standard Gibbs free energy minimum. The impact of quantum zero-point vibrational corrections on the relative stability of these isomers is quite significant. The box-like isomers are favored in terms of electronic energy, but with the inclusion of zero-point vibrational corrections and entropic effects tree-like isomers are favored at higher temperatures. Under conditions from 0 to 298.15 K, the global minimum is predicted to be a tree-like structure with one dangling singly coordinated water molecule. Above 298.15 K, higher entropy tree-like isomers with two or more singly coordinated water molecules are favored. These assignments are generally consistent with experimental IR spectra of (H3O+)(H2O)8 obtained at 150 K.

Single-Molecule Junctions Based on Nitrile-Terminated Biphenyls: A Promising New Anchoring Group

We present a combined experimental and theoretical study of the electronic transport through single-molecule junctions based on nitrile-terminated biphenyl derivatives. Using a scanning tunneling microscope-based break-junction technique, we show that the nitrile-terminated compounds give rise to well-defined peaks in the conductance histograms resulting from the high selectivity of the N-Au binding. Ab initio calculations have revealed that the transport takes place through the tail of the LUMO. Furthermore, we have found both theoretically and experimentally that the conductance of the molecular junctions is roughly proportional to the square of the cosine of the torsion angle between the two benzene rings of the biphenyl core, which demonstrates the robustness of this structure-conductance relationship.

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.

The calculation of Afρ and mass loss rate for comets

Ab initio calculations of Afρ are presented using Mie scattering theory and a Direct Simulation Monte Carlo (DSMC) dust outflow model in support of the Rosetta mission and its target 67P/Churyumov-Gerasimenko (CG). These calculations are performed for particle sizes ranging from 0.010 μm to 1.0 cm. The present status of our knowledge of various differential particle size distributions is reviewed and a variety of particle size distributions is used to explore their effect on Afρ , and the dust mass production View the MathML sourcem˙. A new simple two parameter particle size distribution that curtails the effect of particles below 1 μm is developed. The contributions of all particle sizes are summed to get a resulting overall Afρ. The resultant Afρ could not easily be predicted a priori and turned out to be considerably more constraining regarding the mass loss rate than expected. It is found that a proper calculation of Afρ combined with a good Afρ measurement can constrain the dust/gas ratio in the coma of comets as well as other methods presently available. Phase curves of Afρ versus scattering angle are calculated and produce good agreement with observational data. The major conclusions of our calculations are: – The original definition of A in Afρ is problematical and Afρ should be: qsca(n,λ)×p(g)×f×ρqsca(n,λ)×p(g)×f×ρ. Nevertheless, we keep the present nomenclature of Afρ as a measured quantity for an ensemble of coma particles.– The ratio between Afρ and the dust mass loss rate View the MathML sourcem˙ is dominated by the particle size distribution. – For most particle size distributions presently in use, small particles in the range from 0.10 to 1.0 μm contribute a large fraction to Afρ. – Simplifying the calculation of Afρ by considering only large particles and approximating qsca does not represent a realistic model. Mie scattering theory or if necessary, more complex scattering calculations must be used. – For the commonly used particle size distribution, dn/da ∼ a−3.5 to a−4, there is a natural cut off in Afρ contribution for both small and large particles. – The scattering phase function must be taken into account for each particle size; otherwise the contribution of large particles can be over-estimated by a factor of 10. – Using an imaginary index of refraction of i = 0.10 does not produce sufficient backscattering to match observational data. – A mixture of dark particles with i ⩾ 0.10 and brighter silicate particles with i ⩽ 0.04 matches the observed phase curves quite well. – Using current observational constraints, we find the dust/gas mass-production ratio of CG at 1.3 AU is confined to a range of 0.03–0.5 with a reasonably likely value around 0.1.

Microwave Spectroscopy and Structures of Several Substituted Acetylenes and Fluorocarbons

The structures of two substituted acetylene compounds have been characterized from their microwave rotational spectra. In the first study, two structures of 6-methyl-3-heptyne have been determined. This compound can be thought of as an ethyl group separated from an isobutyl group by a C≡C spacer. Both structures have the ethyl and isobutyl groups eclipsed, consistent with the dominant interaction determining the orientation about the acetylene axis being the weak dispersion attraction between the end groups. One structure is with the isobutyl group in a symmetric conformation and the other with the isobutyl group asymmetric. In addition, the microwave spectrum of the butane analogue 3,5-octadiyne has been observed. This compound consists of two ethyl groups separated by two C≡C spacers. The study is still in progress, but it appears that the ethyl end groups are freely rotating. Therefore, it seems that the dispersion attractions between the end groups are too weak at this longer distance of about 7 Å. The structures of several fluorocarbons have also been studied by microwave spectroscopy. The structures of perfluoropentane and perfluorohexane have been shown to be helical, like the polymer polytetrafluoroethylene (Teflon©). The structure of perfluoropropane and two conformers of 1H-heptafluoropropane have been determined to be non-helical. It is apparent that the steric and dipole repulsions between fluorine atoms that have been attributed to the helical structure of longer fluorocarbon chains are not sufficient in a three carbon chain to cause a twist in the structures.

CCSD(T) Study of CD3-O-CD3 and CH3-O-CD3 Far-Infrared Spectra

From a vibrationally corrected 3D potential energy surface determined with highly correlated ab initio calculations (CCSD(T)), the lowest vibrational energies of two dimethyl-ether isotopologues, 12CH3–16O–12CD3 (DME-d3) and 12CD3–16O–12CD3 (DME-d6), are computed variationally. The levels that can be populated at very low temperatures correspond to the COC-bending and the two methyl torsional modes. Molecular symmetry groups are used for the classification of levels and torsional splittings. DME-d6 belongs to the G36 group, as the most abundant isotopologue 12CH3–16O–12CH3 (DME-h6), while DME-d3 is a G18 species. Previous assignments of experimental Raman and far-infrared spectra are discussed from an effective Hamiltonian obtained after refining the ab initio parameters. Because a good agreement between calculated and experimental transition frequencies is reached, new assignments are proposed for various combination bands corresponding to the two deuterated isotopologues and for the 020 → 030 transition of DME-d6. Vibrationally corrected potential energy barriers, structural parameters, and anharmonic spectroscopic parameters are provided. For the 3N – 9 neglected vibrational modes, harmonic and anharmonic fundamental frequencies are obtained using second-order perturbation theory by means of CCSD and MP2 force fields. Fermi resonances between the COC-bending and the torsional modes modify DME-d3 intensities and the band positions of the torsional overtones.

Band gap control via tuning of inversion degree in CdIn2S4 spinel

Based on theoretical arguments, we propose a possible route for controlling the band-gap in the promising photovoltaic material CdIn2S4. Our ab initio calculations show that the experimental degree of inversion in this spinel (fraction of tetrahedral sites occupied by In) corresponds approximately to the equilibrium value given by the minimum of the theoretical inversion free energy at a typical synthesis temperature. Modification of this temperature, or of the cooling rate after synthesis, is then expected to change the inversion degree, which in turn sensitively tunes the electronic band-gap of the solid, as shown here by screened hybrid functional calculations.

Nonhydrogenic exciton spectrum in perovskite CH3NH3PbI3

The excitons in the orthorhombic phase of the perovskite CH3NH3PbI3 are studied using the effective mass approximation. The electron–hole interaction is screened by a distance-dependent dielectric function, as described by the Haken potential or the Pollmann–Büttner potential. The energy spectrum and the eigenfunctions are calculated for both cases. The results show that the Pollmann–Büttner model, using the corresponding parameters obtained from ab initio calculations, provides better agreement with the experimental results.