Identification and Characterization of the HCl–DMS Gas Phase Molecular Complex via Infrared Spectroscopy and Electronic Structure Calculations

Thermodynamic and optical properties of HCOOH(H2O)n and HCOOH(NH3)(H2O)(n-1) clusters at various temperatures and pressures: a computational study

Density functional theory has been used to compute the gas-phase geometries, binding energies, ZPE-corrected binding energies, BSSE-corrected binding energies, binding enthalpies, and binding free energies of HCOOH(H₂O)_n and HCOOH(NH₃)(H₂O)_(n-1) clusters with n = 1-8, 10, 12, 14, 16, 18, and 20. Enthalpies and free energies are calculated for a range of atmospherically relevant temperatures (T) and pressures (P) (from T = 298.15 K, P = 1013.25 hPa to T = 216.65 K, P = 226.32 hPa). The optical properties of those clusters have been studied at the CAM-B3LYP/aug-cc-pVDZ level of theory.

Thermodynamic properties of forming methanol-water and ethanol-water clusters at various temperatures and pressures and implications for atmospheric chemistry: A DFT study

The gas-phase geometries, binding energies, enthalpies, and free energies of methanol-(water)_n and ethanol-(water)_n clusters containing n=1-10,20,30,40, and 50 water molecules have been calculated using density functional theory. The binding energies are calculated at 0 K. The enthalpies are calculated at a temperature of 298.15 K and pressure of 1013.25 hPa (1 atm). The free energies are calculated at a wide range of temperature (T) and pressure (P) (from T = 298.15 K, P = 1013.25 hPa to T = 216.65 K, P = 226.32 hPa). The results show that the free energy of the formation of a specific cluster from its free molecules is negative (i.e., favorable) only below some critical temperature and pressure, which depends on the cluster's size. One of the most common volatile organic compounds (VOCs) in the troposphere is methanol, ethanol, and atmospheric aerosols containing methanol and ethanol. The Rayleigh scattering properties of methanol-water and ethanol-water clusters have been investigated. The scattering intensities were computed at static (∞ nm) and different wavelengths (700, 600, 500, and 400 nm) of naturally polarized light. Rayleigh scattering intensities increase about 9%-10% at 400 nm compared to the static limit (∞ nm) for both methanol-water and ethanol-water clusters.

Role of glycine on sulfuric acid-ammonia clusters formation: Transporter or participator

Glycine (Gly) is ubiquitous in the atmosphere and plays a vital role in new particle formation (NPF). However, the potential mechanism of its on sulfuric acid (SA) - ammonia (A) clusters formation under various atmospheric conditions is still ambiguous. Herein, a (Gly)_x·(SA)_y·(A)_z (z ≤ x + y ≤ 3) multicomponent system was investigated by using density functional theory (DFT) combined with Atmospheric Cluster Dynamics Code (ACDC) at different temperatures and precursor concentrations. The results show that Gly, with one carboxyl (-COOH) and one amine (-NH₂) group, can interact strongly with SA and A in two directions through hydrogen bonds or proton transfer. Within the relevant range of atmospheric concentrations, Gly can enhance the formation rate of SA-A-based clusters, especially at low temperature, low [SA], and median [A]. The enhancement (R) of Gly on NPF can be up to 340 at T = 218.15 K, [SA] = 10⁴, [A] = 10⁹, and [Gly] = 10⁷ molecules/cm³. In addition, the main growth paths of clusters show that Gly molecules participate into cluster formation in the initial stage and eventually leave the cluster by evaporation in subsequent cluster growth at low [Gly], it acts as an important "transporter" to connect the smaller and larger cluster. With the increase of [Gly], it acts as a "participator" directly participating in NPF.

Can nitrous acid contribute to atmospheric new particle formation from nitric acid and water?

The properties of (HNO₃)(HONO)(H₂O)_n (n = 1–6) clusters are reported including thermodynamics, structures, temperature-dependence, intermolecular forces, optical properties, and evaporation rates.

Theoretical study of the cis-pinonic acid and its atmospheric hydrolysate participation in the atmospheric nucleation

cis-Pinonic acid (CPA), one of the major photooxidation products of α-pinene, is believed to contribute to the formation of aerosols formed over forested areas. In the current study, we implement quantum chemical calculation to investigate the interaction between sulfuric acid (SA) and CPA as well as the hydrolysate of CPA (HCPA) in the presence of water or ammonia in the atmosphere. The lowest free energy configurations, reactants, transition states, intermediates, and products were optimized at 298/278K and 1atm at the M06-2X/6-311+G(3df,3pd) level. Our results show that one CPA molecule might initially nucleate with SA molecules and subsequently participate in the formation and growth of the new particle in the form of HCPA. More than one HCPA molecule may be involved in the critical nuclei. Furthermore, the hydrolysis reaction of CPA can be effectively catalyzed by SA and nitric acid (NA) in presence of water, which significantly increases the HCPA content in the atmosphere and subsequently promotes the particle nucleation. Overall, the current study elucidates a new mechanism of atmospheric nucleation driven by CPA and its hydrolysate.

A density functional theory study of aldehydes and their atmospheric products participating in nucleation

The products of aldehydes from aldol condensation, hydration, and polymerization reactions can promote new particle formation by stabilizing sulfuric acid.

Formation of atmospheric molecular clusters consisting of sulfuric acid and C8H12O6 tricarboxylic acid

We present the structures and thermochemical properties of (MBTCA)₁₋₃(H₂SO₄)₁₋₄ atmospheric molecular clusters.

Basis set convergence of the binding energies of strongly hydrogen-bonded atmospheric clusters

We present the first binding energy benchmark set at the CBS limit of strongly hydrogen bonded atmospheric molecular clusters.

Interaction of oxalic acid with dimethylamine and its atmospheric implications

Oxalic acid and dimethylamine are the most common organic acid and base in the atmosphere, and are recognized as significant precursor species in atmospheric new particle formation.

Acid-Base Formalism Extended to Excited State for O-H···S Hydrogen Bonding Interaction

Hydrogen bond can be regarded as an interaction between a base and a proton covalently bound to another base. In this context the strength of hydrogen bond scales with the proton affinity of the acceptor base and the pKa of the donor, i.e., it follows the acid-base formalism. This has been amply demonstrated in conventional hydrogen bonds. Is this also true for the unconventional hydrogen bonds involving lesser electronegative elements such as sulfur atom? In our previous work, we had established that the strength of O-H···S hydrogen bonding (HB) interaction scales with the proton affinity (PA) of the acceptor. In this work, we have investigated the other counterpart, i.e., the H-bonding interaction between the photoacids with different pKa values with a common base such as the H2O and H2S. The 1:1 complexes of five para substituted phenols p-aminophenol, p-cresol, p-fluorophenol, p-chlorophenol, and p-cyanophenol with H2O and H2S were investigated experimentally and computationally. The investigations were also extended to the excited states. The experimental observations of the spectral shifts in the O-H stretching frequency and the S1-S0 band origins were correlated with the pKa of the donors. Ab initio calculations at the MP2 and various dispersion corrected density functional levels of theory were performed to compute the dissociation energy (D0) of the complexes. The quantum theory of atoms in molecules (QTAIM), noncovalent interaction (NCI) method, natural bonding orbital (NBO) analysis, and natural decomposition analysis (NEDA) were carried out for further characterization of HB interaction. The O-H stretching frequency red shifts and the dissociation energies were found to be lower for the O-H···S hydrogen bonded systems compared to those for the O-H···O H-bound systems. Despite being dominated by the dispersion interaction the O-H···S interaction in the H2S complexes also conformed to the acid-base formalism, i.e., the D0 and the O-H red shift scaled with the pKa of the donor, similar to that observed in the O-H···O interaction. However, the two classes of H-bonds follow different correlations. In addition we also discuss the nuances associated with the similarity and differences in the hydrogen bonding properties of the two classes in the ground electronic state as well as in the excited state.

Spectroscopy and dynamics of double proton transfer in formic acid dimer

We present the isolated gas phase infrared spectra of formic acid dimer, (HCOOH)2, and its deuterated counterpart formic-d acid, (DCOOH)2, at room temperature. The formic acid dimer spectrum was obtained by spectral subtraction of a spectrum of formic acid vapor recorded at low pressure from that recorded at a higher pressure. The spectra of formic acid vapor contain features from both formic acid monomer and formic acid dimer, but at low and high pressures of formic acid, the equilibrium is pushed towards the monomer and dimer, respectively. A similar approach was used for the formic-d acid dimer. Building on the previous development of the Molecular Mechanics with Proton Transfer (MMPT) force field for simulating proton transfer reactions, molecular dynamics (MD) simulations were carried out to interpret the experimental spectra in the OH-stretching region. Within the framework of MMPT, a combination of symmetric single and double minimum potential energy surfaces (PESs) provides a good description of the double proton transfer PES. In a next step, potential morphing together with electronic structure calculations at the B3LYP and MP2 level of theory was used to align the computed and experimentally observed spectral features in the OH-stretching region. From this analysis, a barrier for double proton transfer between 5 and 7 kcal mol(-1) was derived, which compares with a CCSD(T)/aug-cc-pVTZ calculated barrier of 7.9 kcal mol(-1). Such a combination of experimental and computational techniques for estimating barriers for proton transfer in gas phase systems is generic and holds promise for further improved PESs and energetics of these important systems. Additional MD simulations at the semi-empirical DFTB level of theory agree quite well for the center band position but underestimate the width of the OH-stretching band.

Computational (DFT and MP2) and spectral interpretations, normal coordinate analysis, force constants and barriers to internal rotations of Trimethylacetonitrile

The Raman (3500[Formula: see text]cm[Formula: see text]100[Formula: see text]cm[Formula: see text] and IR spectra (4000[Formula: see text]cm[Formula: see text]400[Formula: see text]cm[Formula: see text] of liquid trimethylacetonitrile (C5H9N, TMA) have been obtained. In addition, the 1H and [Formula: see text]C NMR spectra of TMA were obtained in DMSO-[Formula: see text] and CDCl3. The staggered conformer (C[Formula: see text] was favored using MP2 and DFT(B3LYP/[Formula: see text]B97XD) quantum mechanical calculations utilizing a 6-311[Formula: see text]G(d,p) basis set. High energy difference estimates of 4534[Formula: see text]cm[Formula: see text]5338[Formula: see text]cm[Formula: see text] (12.96[Formula: see text]kcal/mol[Formula: see text]15.26[Formula: see text]kcal/mol) were predicted, along with three imaginary torsion frequencies for the eclipsed conformer, therefore considered a transition state. The 1H and [Formula: see text]C NMR chemical shifts were predicted with B3LYP and [Formula: see text]B97XD methods using the GIAO approximation and 6-311[Formula: see text]G(d,p) basis set. B3LYP frequencies calculation is favored herein owing to the relatively good compilation with the experimental measurements. The computed structural parameters are well correlated to those reported from electron diffraction and microwave studies. Moreover, the [Formula: see text]C[Formula: see text]H coupling constant was estimated and found consistent with that observed for the sample dissolved in DMSO-[Formula: see text]/CDCl3 solvents. Using the observed methyl torsion at 266[Formula: see text]cm[Formula: see text] in gas phase and the kinetic parameter F number, a potential function (V[Formula: see text] of [Formula: see text][Formula: see text]cm[Formula: see text]([Formula: see text][Formula: see text]kcal/mol) was obtained, this barrier to internal rotation is well correlated to 1527[Formula: see text]cm[Formula: see text] (4.37[Formula: see text]kcal/mol) predicted from MP2/6-311[Formula: see text]G(d,p) potential surface scan. Aided by the predicted wavenumbers and their IR intensity/Raman activity, the observed IR/Raman bands were intensively discussed and therefore assigned to their corresponding fundamentals, in agreement with novel normal coordinate analysis and potential energy distributions (PEDs).

Strong Hydrogen Bonded Molecular Interactions between Atmospheric Diamines and Sulfuric Acid

We investigate the molecular interaction between methyl-substituted N,N,N',N'-ethylenediamines, propane-1,3-diamine, butane-1,4-diamine, and sulfuric acid using computational methods. Molecular structure of the diamines and their dimer clusters with sulfuric acid is studied using three density functional theory methods (PW91, M06-2X, and ωB97X-D) with the 6-31++G(d,p) basis set. A high level explicitly correlated CCSD(T)-F12a/VDZ-F12 method is used to obtain accurate binding energies. The reaction Gibbs free energies are evaluated and compared with values for reactions involving ammonia and atmospherically relevant monoamines (methylamine, dimethylamine, and trimethylamine). We find that the complex formation between sulfuric acid and the studied diamines provides similar or more favorable reaction free energies than dimethylamine. Diamines that contain one or more secondary amino groups are found to stabilize sulfuric acid complexes more efficiently. Elongating the carbon backbone from ethylenediamine to propane-1,3-diamine or butane-1,4-diamine further stabilizes the complex formation with sulfuric acid by up to 4.3 kcal/mol. Dimethyl-substituted butane-1,4-diamine yields a staggering formation free energy of -19.1 kcal/mol for the clustering with sulfuric acid, indicating that such diamines could potentially be a key species in the initial step in the formation of new particles. For studying larger clusters consisting of a diamine molecule with up to four sulfuric acid molecules, we benchmark and utilize a domain local pair natural orbital coupled cluster (DLPNO-CCSD(T)) method. We find that a single diamine is capable of efficiently stabilizing sulfuric acid clusters with up to four acid molecules, whereas monoamines such as dimethylamine are capable of stabilizing at most 2-3 sulfuric acid molecules.

The Effect of Water and Bases on the Clustering of a Cyclohexene Autoxidation Product C6H8O7 with Sulfuric Acid

We investigate the molecular interaction between sulfuric acid and a C6H8O7 ketodiperoxy acid compound (a proxy for highly oxidized products of, e.g., monoterpene autoxidation) in the presence of water, ammonia, or dimethylamine. The molecular geometries are obtained using density functional theory (M06-2X, PW91, and ωB97X-D) with the 6-31++G(d,p) basis set, and the binding energy is corrected utilizing a high-level DLPNO-CCSD(T)/def2-QZVPP calculation. The formation free energies were calculated (ΔG at 298 K and 1 atm), and the stability of the molecular clusters was evaluated. The presence of bases is found to enhance the interaction between ketodiperoxy acid compounds and sulfuric acid. The addition of C6H8O7 compounds to (H2SO4)(NH3) or (H2SO4)((CH3)2NH) clusters is, however, not able to compete with the corresponding uptake of another sulfuric acid molecule, even at a high loading of organic compounds. We furthermore investigate the origin of the weak binding of peroxyacid compounds using atoms in molecules and natural bonding orbital analysis. The weak binding is caused by an internal hydrogen bond and the lack of a strong hydrogen bond acceptor in the peroxyacid group. These findings indicate that autoxidation products containing solely or mainly hydroperoxide and carbonyl functional groups do not participate in the initial step of new particle formation and thereby only contribute to the growth of atmospheric particles.

Predicting accurate absolute binding energies in aqueous solution: thermodynamic considerations for electronic structure methods

Recent predictions of absolute binding free energies of host-guest complexes in aqueous solution using electronic structure theory have been encouraging for some systems, while other systems remain problematic. In this paper I summarize some of the many factors that could easily contribute 1-3 kcal mol(-1) errors at 298 K: three-body dispersion effects, molecular symmetry, anharmonicity, spurious imaginary frequencies, insufficient conformational sampling, wrong or changing ionization states, errors in the solvation free energy of ions, and explicit solvent (and ion) effects that are not well-represented by continuum models. While I focus on binding free energies in aqueous solution the approach also applies (with minor adjustments) to any free energy difference such as conformational or reaction free energy differences or activation free energies in any solvent.

Standard free energy of the equilibrium between the trans-monomer and the cyclic-dimer of acetic acid in the gas phase from infrared spectroscopy

Survey jet-cooled spectra of acetic acid have been recorded in the infrared region (200-4000 cm(-1)) over a wide range of expansion conditions. From the variations of the relative intensities of the signals, vibrational transitions have been assigned unambiguously to the trans-monomer and cyclic-dimer. The IR-active fundamental frequencies have been determined at the instrumental accuracy of 0.5 cm(-1). This analysis of the jet-cooled spectra supported by electronic structure calculations permitted us to characterize the trans-monomer/cyclic-dimer equilibrium. From static cell spectra at 298 K, variations of the molar fractions ratio as a function of the total pressure were used to estimate the equilibrium constant and the Gibbs free energy of dimerization at 298 K. The very good agreement with the literature data shows that the present method is able to produce, from a single study, a free energy value as reliable as the one obtained from a large collection of data. In addition, the semi-empirical free energy value was used to estimate the accuracy of electronic structure calculations and in turn the accuracy of the derived useful information such as the dissociation energy of the complex (i.e. the strength of the hydrogen bonds) or the relative energies within the conformational landscape.

Computational Study of the Clustering of a Cyclohexene Autoxidation Product C6H8O7 with Itself and Sulfuric Acid

We investigate the molecular interactions between sulfuric acid and a recently reported C6H8O7 ketodiperoxy acid formed through autoxidation from cyclohexene and ozone. Structurally similar but larger ELVOC (extremely low volatility organic compound) products formed from autoxidation of monoterpenes are believed to play a major role in the formation and early growth of atmospheric aerosol particles. Utilizing density functional theory geometries, with a DLPNO-CCSD(T)/def2-QZVPP single point energy correction, the stepwise Gibbs free energies of formation have been calculated, and the stabilities of the molecular clusters have been evaluated. C6H8O7 interacts weakly with both itself and sulfuric acid, with standard free energies of formation (ΔG at 298 K and 1 atm) around or above 0 kcal/mol. This is due to the presence of strong intramolecular hydrogen bonds in the peroxyacid groups of C6H8O7. These stabilize the isolated molecule with respect to its clusters, and lead to unfavorable interaction energies. The addition of sulfuric acid to clusters containing C6H8O7 is somewhat more favorable, but the formed clusters are still far more likely to evaporate than to grow further in atmospheric conditions. These findings indicate that the O/C ratio cannot exclusively be used as a proxy for volatility in atmospheric new particle formation involving organic compounds. The specific molecular structure, and especially the number of strong hydrogen binding moieties, are equally important. The interaction between the C6H8O7 compound and aqueous phase sulfate ions indicates that ELVOC-type compounds can contribute to aerosol mass by effectively partitioning into the condensed phase.

Structure, stability and intramolecular interaction of M(N5)2(M = Mg, Ca, Sr and Ba) : a theoretical study

The potential energetic materials, alkaline earth metal complexes of the pentazole anion (M(N₅)₂, M = Mg²⁺, Ca²⁺, Sr²⁺and Ba²⁺), were studied using the density functional theory.

Rayleigh light scattering properties of atmospheric molecular clusters consisting of sulfuric acid and bases

The Rayleigh light scattering properties of (H₂SO₄)_a(NH₃)_b and (H₂SO₄)_a((CH₃)₂NH)_b atmospheric molecular clusters have been investigated using a response theory approach.

Benchmarking ab initio binding energies of hydrogen-bonded molecular clusters based on FTIR spectroscopy

Models of formation and growth of atmospheric aerosols are highly dependent on accurate cluster binding energies. These are most often calculated by ab initio electronic structure methods but remain associated with significant uncertainties. We present a computational benchmarking study of the Gibbs free binding energies in molecular complexes and clusters based on gas phase FTIR spectroscopy. The acetonitrile-HCl molecular complex is identified via its redshifted H-Cl stretching vibrational mode. We determine the Gibbs free binding energy, ΔG°295 K, to between 4.8 and 7.9 kJ mol(-1) and compare this range to predictions from several widely used electronic structure methods, including five density functionals, Møller-Plesset perturbation theory, and five coupled cluster methods up to CCSDT quality, considering also the D3 dispersion correctional scheme. With some exceptions, we find that most electronic structure methods overestimate ΔG°295 K. The effects of vibrational anharmonicity is approximated using scaling factors, reducing ΔG°295 K by ca. 1.8 kJ mol(-1), whereby ΔG°295 K predictions well within the experimental range can be obtained.

Computational study of the Rayleigh light scattering properties of atmospheric pre-nucleation clusters

The Rayleigh and hyper Rayleigh scattering properties of the binary (H2SO4)(H2O)n and ternary (H2SO4)(NH3)(H2O)n clusters are investigated using a quantum mechanical response theory approach. The molecular Rayleigh scattering intensities are expressed using the dipole polarizability α and hyperpolarizability β tensors. Using density functional theory, we elucidate the effect of cluster morphology on the scattering properties using a combinatorial sampling approach. We find that the Rayleigh scattering intensity depends quadratically on the number of water molecules in the cluster and that a single ammonia molecule is able to induce a high anisotropy, which further increases the scattering intensity. The hyper Rayleigh scattering activities are found to be extremely low. This study presents the first attempt to map the scattering of atmospheric molecular clusters using a bottom-up approach.

The effect of fluorine substitution in alcohol–amine complexes

The effect of fluorine substitution on the hydrogen bond strength in alcohol–amine molecular complexes was investigated, with a combination of vapour phase infrared and near infrared spectroscopy and theoretical calculations.

Resolving the anomalous infrared spectrum of the MeCN–HCl molecular cluster using ab Initio molecular dynamics

Molecular dynamics simulations reveal that the asymmetric peak seen in IR spectra of acetonitrile–HCl molecular complexes are due to high population of complexes with partially broken hydrogen bonds.