Reparameterization of GFN1-xTB for atmospheric molecular clusters: applications to multi-acid-multi-base systems

Atmospheric molecular clusters, the onset of secondary aerosol formation, are a major part of the current uncertainty in modern climate models. Quantum chemical (QC) methods are usually employed in a funneling approach to identify the lowest free energy cluster structures. However, the funneling approach highly depends on the accuracy of low-cost methods to ensure that important low-lying minima are not missed. Here we present a reparameterized GFN1-xTB model based on the clusteromics I-V datasets for studying atmospheric molecular clusters (AMC), denoted AMC-xTB. The AMC-xTB model reduces the mean of electronic binding energy errors from 7-11.8 kcal mol-1 to roughly 0 kcal mol-1 and the root mean square deviation from 7.6-12.3 kcal mol-1 to 0.81-1.45 kcal mol-1. In addition, the minimum structures obtained with AMC-xTB are closer to the ωB97X-D/6-31++G(d,p) level of theory compared to GFN1-xTB. We employ the new parameterization in two new configurational sampling workflows that include an additional meta-dynamics sampling step using CREST with the AMC-xTB model. The first workflow, denoted the "independent workflow", is a commonly used funneling approach with an additional CREST step, and the second, the "improvement workflow", is where the best configuration currently known in the literature is improved with a CREST + AMC-xTB step. Testing the new workflow we find configurations lower in free energy for all the literature clusters with the largest improvement being up to 21 kcal mol-1. Lastly, by employing the improvement workflow we massively screened 288 new multi-acid-multi-base clusters containing up to 8 different species. For these new multi-acid-multi-base cluster systems we observe that the improvement workflow finds configurations lower in free energy for 245 out of 288 (85.1%) cluster structures. Most of the improvements are within 2 kcal mol-1, but we see improvements up to 8.3 kcal mol-1. Hence, we can recommend this new workflow based on the AMC-xTB model for future studies on atmospheric molecular clusters.

Current and future machine learning approaches for modeling atmospheric cluster formation

The formation of strongly bound atmospheric molecular clusters is the first step towards forming new aerosol particles. Recent advances in the application of machine learning models open an enormous opportunity for complementing expensive quantum chemical calculations with efficient machine learning predictions. In this Perspective, we present how data-driven approaches can be applied to accelerate cluster configurational sampling, thereby greatly increasing the number of chemically relevant systems that can be covered.

Hydrogen-Bonding Interactions of Malic Acid with Common Atmospheric Bases

Malic acid (MA) (C₄H₆O₅) is one of the most important organic constituents of fruits that is widely used in food and beverage industries. It is also detected in the atmospheric aerosol samples collected in different parts of the world. Considering the fact that secondary organic aerosols have adverse impacts on the global atmosphere and climate and a molecular-level understanding of the compositions and formation mechanism of secondary organic aerosols is necessary, we have performed systematic density functional electronic structure calculations to investigate the hydrogen-bonding interactions between MA and several naturally occurring nitrogen-containing atmospheric bases such as ammonia and amines that are derived from ammonia by the substitution of hydrogens by a methyl group. The base molecules were allowed to interact with the carboxylic COOH and the hydroxyl-OH group of the MA separately. While at both sites, MA produces energetically stable binary complexes with bases with large negative values of binding energy, the thermodynamical stability, at an ambient temperature and pressure of 298.15 K and 1 atm, respectively, is favored only for the clusters formed at the COOH site. A much larger red shift of the carboxylic-OH stretch than that of the hydroxyl-OH reinforces the preference of this site for cluster formation. Both the binding electronic energy and binding free energy of MA-ammonia complexes are lower than those of MA-amine complexes, although the amines are derivatives of NH₃. The large increase in the Rayleigh activities upon cluster formation indicates that the MA-atmospheric base cluster may interact strongly with solar radiation. The detailed analysis of the structural, energetic, electrical, and spectroscopic properties of the binary complexes formed by MA with atmospheric bases shows that MA could participate in the atmospheric nucleation processes and subsequently contribute effectively to new particle formation in the atmosphere.

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.

Computational chemistry of cluster: Understanding the mechanism of atmospheric new particle formation at the molecular level

New particle formation (NPF), which exerts significant influence over human health and global climate, has been a hot topic and rapidly expands field of research in the environmental and atmospheric chemistry recent years. Generally, NPF contains two processes: formation of critical nucleus and further growth of the nucleus. However, due to the complexity of the atmospheric nucleation, which is a multicomponent process, formation of critical clusters as well as their growth is still connected to large uncertainties. Detection limits of instruments in measuring specific gaseous aerosol precursors and chemical compositions at the molecular level call for computational studies. Computational chemistry could effectively compensate the deficiency of laboratory experiments as well as observations and predict the nucleation mechanisms. We review the present theoretical literatures that discuss nucleation mechanism of atmospheric clusters. Focus of this review is on different nucleation systems involving sulfur-containing species, nitrogen-containing species and iodine-containing species. We hope this review will provide a deep insight for the molecular interaction of nucleation precursors and reveal nucleation mechanism at the molecular level.

Interaction of Glutamic Acid/Protonated Glutamic Acid with Amide and Water Molecules: A Theoretical Study

Amino acids are important nitrogen-containing compounds and organic carbon components that exist widely in the atmosphere. The formation of atmospheric aerosols is affected by their interactions with amides. The dimers formed by glutamic acid (Glu) or protonated glutamic acid (Glu⁺) with three kinds of amide molecules (formamide FA, acetamide AA, urea U) and the hydrated clusters formed by Glu or Glu⁺, U molecules along with one to six water molecules were systematically studied at the M06-2X/6-311++G(3df,3pd) level. U is predicted to form a more stable structure with Glu/Glu⁺ than FA and AA by thermodynamics. If the concentration ratio of FA to U is less than 10⁴, U will play a critical role in NPF. The degree of hydration in Glu⁺-mU-nW is higher than that of Glu-mU-nW (m = 0, 1; n = 0-6) clusters. Notably, Glu contributes more to the Rayleigh scattering properties than glutaric acid and sulfuric acid, and thus may lead to the destruction of atmospheric visibility. This study is helpful to better understand the properties of organic aerosols containing amino acids or amides.

Theoretical analysis of sulfuric acid–dimethylamine–oxalic acid–water clusters and implications for atmospheric cluster formation

In recent years, organic compounds potentially involved in atmospheric particle formation have received increased attention. However, the contributions of organic acids as precursors in nucleation remain ambiguous. In this study, the low-lying structures and thermodynamics of the sulfuric acid–dimethylamine–oxalic acid–water system are obtained at the M06-2X/6-311+G(2d,p) level, and the single point energy of the clusters has been calculated at the DF-LMP2-F12/VDZ-F12 level. The formations of the multicomponent clusters are predicted based on thermodynamics, involving proton transfer and hydrogen bonding interactions. Oxalic acid can synergistically promote the formation of the sulfuric acid–dimethylamine–oxalic acid–water system while inhibiting this with the addition of more sulfuric acid molecules. The results of hydrate distribution show that un-hydrate clusters play a dominant role during formation. Moreover, dimethylamine and oxalic acid have similar effects on Rayleigh scattering properties, and the clusters involving complex mixtures of compounds can have high optical activities.

The structure of SA₂.DMA.OA.W₄ cluster.

Characterization of dicarboxylic acids, oxoacids, and α-dicarbonyls in PM2.5 within the urban boundary layer in southern China: Sources and formation pathways

Low-molecular-weight dicarboxylic acids, which are important components of secondary organic aerosols, have been extensively studied in recent years. Many studies have focused on ground-level observations and literature reports on the vertical distribution of the organic aerosols within the urban boundary layer are limited. In this study, the vertical profiles of dicarboxylic acids and related organic compounds (DCRCs) in PM_2.5 were investigated at altitudinal levels (ground level and 488 m above the ground level) at the Canton Tower in Guangzhou, southern China, to elucidate their primary sources and secondary formation processes. The concentrations of DCRCs at ground level were generally higher than those at 488 m. Oxalic acid (C₂) was the most abundant species, followed by succinic acid (C₄) and malonic acid (C₃) at both heights. The higher ratio of DCRCs-bound carbon to organic carbon (i.e., DCRCs-C/OC) at 488 m (4.8 ± 1.2%) relative to that at ground level (2.7 ± 0.5%) indicated a higher degree of aerosol aging at 488 m. The abundance of C₂ was increased and the conversion of C₄ to C₃ was enhanced due to the photochemical oxidation of its homologues during long-range transport periods. The increase in C₂ was associated with in-cloud processes during pollution periods. Principal component analysis showed that DCRCs were mainly derived from atmospheric secondary processing and biomass burning was also an important source of long-chain carboxylic acids during autumn in Guangzhou. Our results illustrate that secondary processing and biomass burning play prominent roles in controlling the abundance of DCRCs. Furthermore, DCRCs are affected by air masses from regional areas, oxidation of their precursors via vertical transport and in-cloud processes.

Synergistic effect of glutaric acid and ammonia/amine/amide on their hydrates in the clustering: A theoretical study

The formation of molecular clusters makes influence on the atmosphere. The clusters of glutaric acid (GA) and common ammonia (A), amine (methylamine MA, dimethylamine DMA) and representative amide (urea U) along with water molecule were systematically studied theoretically. GA-A-nW (n = 1, 2), GA-MA-nW (n = 1, 2), GA-DMA-1W and GA-U-nW (n = 1-6) are predicted to be feasible thermodynamically with the hydrogen bonds as interaction force. GA and urea promote the clustering synergistically, and ammonia, methylamine, dimethylamine promote the clustering of small GA hydrates (n = 1-2), while inhibit that of large GA hydrates (n = 3-6). The results of humidity show that un-hydrate or mono-hydrate is the main form of GA-mbase-nW (m = 0, 1; n = 1-6) under relative humidity of 20%, 50% and 80%. The global minima remain dominant over the temperature range of 220-320 K. GA contributes more to the Rayleigh scattering properties than sulfuric acid. More importantly, the local minima can undergo isomerization to form the global minima crossing a free energy barrier ranging from 6.66 to 11.78 kcal mol^-1. This study indicates that GA and base molecules play a synergistic role to promote the formation of clusters. We hope it can provide more insights on interesting clustering in 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.

Low molecular weight dicarboxylic acids, oxocarboxylic acids and α-dicarbonyls as ozonolysis products of isoprene: Implication for the gaseous-phase formation of secondary organic aerosols

Oxidation of isoprene, a major biogenic volatile organic compound emitted from forest canopies, is a potential source of oxalic acid; the dominant species in organic aerosols. We evaluated here ozonolysis of isoprene in dry darkness as a source of oxalic (C₂), malonic (C₃) and succinic (C₄) acids. We found that oxalic acid and methylglyoxal are dominant products within 10 min of reaction followed by glyoxylic, malonic or succinic acids. Interestingly, molecular distributions of oxidation products from early reactions (9-29 min) were characterized by the predominance of methylglyoxal followed by C₂, which became dominant after 30 min. The isoprene-derived secondary organic aerosols (SOAs) showed chemical evolution with reaction time towards the molecular characteristics of dicarboxylic acids similar to those of ambient aerosols (C₂>C₃≥C₄). The carbon-based relative abundances of methylglyoxal decreased steadily (40%→30%), while those of C₂ increased with reaction time (15%→25%), but no such variations persisted for glyoxal (6-10%). This finding means that methylglyoxal is more important intermediate of oxalic acid than glyoxal. In contrast, smaller variability and lower concentrations of pyruvic and glyoxylic acids than other intermediates indicate that oxalic acid formation under dry conditions follows a different pathway than in aqueous-phase heterogeneous chemistry usually invoked for cloud/fog/atmospheric waters. Here, we propose new reaction schemes for high levels of methylglyoxal and oxalic acid via gas-phase chemical reactions with ozone and OH radicals to better interpret the ambient SOA composition. Furthermore, the relative abundances of C₂ exhibit small variability from 1 to 8 h, suggesting its stable character towards the oxidation by hydroxyl radicals.

Atmospherically Relevant Hydrogen-Bonded Interactions between Methanesulfonic Acid and H2SO4 Clusters: A Computational Study

A detailed and systematic quantum-chemical calculation has been performed with high-level density functional theory (DFT) to analyze the electrostatic interaction of methanesulfonic acid (CH₃SO₃H), also known as MSA, with pre-formed clusters of sulfuric acid (H₂SO₄) molecules in ambient conditions. Both MSA and H₂SO₄ are considered as atmospheric molecules that might play active roles in aerosol formation. The interactions between MSA and H₂SO₄ clusters lead to the formation of MSA···(H₂SO₄)_n (n = 2, 3) complexes stabilized by the formation of different types of intermolecular hydrogen bond networks. Analyses of cluster binding energies and free energy changes associated with their formation indicate that MSA could bring additional stability into the atmospheric molecular clusters responsible for aerosol formation. Variations of Gibbs free energy with temperature and pressure have been analyzed. The lower temperatures and pressures at the higher altitudes of the troposphere are found to play in favor of higher stability of the MSA···(H₂SO₄)_n clusters. Effects of hydrogen bond formation on dipole moment, mean polarizability, and anisotropy of polarizability of the clusters have been analyzed. Rayleigh scattering intensities are found to increase many-fold when light interacts with the MSA···(H₂SO₄)_n clusters.

Integrated experimental and theoretical approach to probe the synergistic effect of ammonia in methanesulfonic acid reactions with small alkylamines

While new particle formation events have been observed worldwide, our fundamental understanding of the precursors remains uncertain. It has been previously shown that small alkylamines and ammonia (NH3) are key actors in sub-3 nm particle formation through reactions with acids such as sulfuric acid (H2SO4) and methanesulfonic acid (CH3S(O)(O)OH, MSA), and that water also plays a role. Because NH3 and amines co-exist in air, we carried out combined experimental and theoretical studies examining the influence of the addition of NH3 on particle formation from the reactions of MSA with methylamine (MA) and trimethylamine (TMA). Experiments were performed in a 1 m flow reactor at 1 atm and 296 K. Measurements using an ultrafine condensation particle counter (CPC) and a scanning mobility particle sizer (SMPS) show that new particle formation was systematically enhanced upon simultaneous addition of NH3 to the MSA + amine binary system, with the magnitude depending on the amine investigated. For the MSA + TMA reaction system, the addition of NH3 at ppb concentrations produced a much greater effect (i.e. order of magnitude more particles) than the addition of ∼12 000 ppm water (corresponding to ∼45-50% relative humidity). The effect of NH3 on the MSA + MA system, which is already very efficient in forming particles on its own, was present but modest. Calculations of energies, partial charges and structures of small cluster models of the multi-component particles likewise suggest synergistic effects due to NH3 in the presence of MSA and amine. The local minimum structures and the interactions involved suggest mechanisms for this effect.

Atmospheric chemistry of thiourea: nucleation with urea and roles in NO2 hydrolysis

The nucleation with urea and roles in NO₂ hydrolysis in the presence of thiourea.

A molecular understanding of the interaction of typical aromatic acids with common aerosol nucleation precursors and their atmospheric implications

Aromatic acids, which are generated from numerous anthropogenic emissions and secondary transformations, have been considered to play a crucial role in new particle formation. In this study, we performed theoretical calculations at the PW91PW91/6-311++G(3df,3pd) level to investigate the interaction between typical aromatic acids namely benzoic acid (BA), phenylacetic acid (PAA), phthalic acid (PA), isophthalic acid (mPA), and terephthalic acid (PTA) and common atmospheric nucleation precursors namely sulfuric acid (SA), water (H2O), ammonia (NH3), methylamine (MA), dimethylamine (DMA), and trimethylamine (TMA). The geometric analysis, Gibbs free energy analysis, OH/NH-stretching vibrational frequency calculation, and atoms in molecules (AIM) analysis were conducted to determine the interactions in the complexes. The heterodimers formed a six to eight membered ring through four types of hydrogen bond, and the bond strength could be ranked in descending order: SO-H⋯O > O-H⋯O/N > N-H⋯O. The BA/PAA/mPA/PTA-SA complexes had the lowest Gibbs free energy values. PA was more likely to interact with NH3 or amines rather than SA due to an intra-molecular hydrogen bond. Additionally, the aromatic acids have similar ability to interact with SA and NH3 as monocarboxylic/dicarboxylic acid. The formation potential of the heterodimers from aromatic acids with common nucleation precursors in ambient atmosphere was investigated.

A molecular-scale study on the hydration of sulfuric acid-amide complexes and the atmospheric implication

Amides are ubiquitous in atmosphere. However, the role of amides in new particle formation (NPF) is poorly understood. Herein, the interaction of urea and formamide with sulfuric acid (SA) and up to four water (W) molecules has been studied at the M06-2X/6-311++G(3df,3pd) level of theory. The structures and properties of (Formamide)(SA)(W)_n (n = 0-4) and (Urea)(SA)(W)_n (n = 0-4) clusters were investigated. Results show that the interaction of SA with the CO group of amides plays a more important role in amide clusters compared with the NH₂ group. Proton transfer to water molecule become dominant in highly hydrated amide clusters at lower temperatures. There is no proton transfer to CO group in formamide clusters. The Rayleigh light scattering intensities of amide clusters are comparable to that of amine and oxalic acid clusters reported previously. Moreover, unhydrated (Amide)(SA) clusters have similar or even higher ability than hydrated SA clusters to participate in ion-induced nucleation. In comparison with formamide, urea has more interacting sites and its clusters have higher Rayleigh light scattering intensities, larger dipole moment, stronger interaction with SA and lower water affinity. The intermolecular interaction in (Formamide)(SA) is slightly weaker than that of SA dimer, which may be compensated by the high concentration of formamide, thus enabling formamide to participate in initial steps of NPF. This study may bring new insight into the role of amides in initial steps of NPF from molecular scale and could help better understand the properties of amide-containing organic aerosol.

Molecular understanding of the interaction of amino acids with sulfuric acid in the presence of water and the atmospheric implication

Amino acids are important components of atmospheric aerosols. Despite the diversity of amino acids structures, however, the role of amino acids with additional non-characteristic functional groups in new particle formation (NPF) has almost remained unexplored. Herein, the interaction of serine (Ser) and threonine (Thr), which feature a hydroxyl group and differ by a methyl-substitution, with sulfuric acid (SA) and up to three water (W) molecules has been investigated at the M06-2X/6-311++G (3df, 3pd) level of theory. The effects of structural differences of amino acids on the structure and properties of clusters were also pointed out. Results show that serine may play more important role in stabilizing sulfuric acid to promote NPF in initial steps compared with threonine, glycine and alanine. Meanwhile, threonine may participate in ion-induced nucleation due to the high dipole moment of (Thr) (SA) isomers. Moreover, the effects of structure differences of amino acids can be seen in several aspects. Firstly, methyl substitution and hydroxyl group of amino acids have great influence on the structure of clusters. Secondly, hydrated (Ser) (SA) and (Tur) (SA) clusters could retain water even at low relative humidity, which may due to the hydroxyl group in serine and threonine. In addition, the Rayleigh light scattering intensities of amino acid-containing clusters are higher than trimethylamine, monoethanolamine and oxalic acid-involved counterparts. The effect of carboxyl group and methyl substitution on optical properties of clusters is also discussed. This study may bring new insight into the role of amino acids with additional non-characteristic functional groups in initial steps of NPF.

Hydration of Atmospheric Molecular Clusters: A New Method for Systematic Configurational Sampling

We present a new systematic configurational sampling algorithm for investigating the potential energy surface of hydrated atmospheric molecular clusters. The algorithm is based on creating a Fibonacci sphere around each atom in the cluster and adding water molecules to each point in nine different orientations. For the sampling of water molecules to existing hydrogen bonds, the cluster is displaced along the hydrogen bond, and a water molecule is placed in between in three different orientations. Generated redundant structures are eliminated based on minimizing the root-mean-square distance of different conformers. Initially, the clusters are sampled using the semiempirical PM6 method and subsequently using density functional theory (M06-2X and ωB97X-D) with the 6-31++G(d,p) basis set. Applying the developed algorithm, we study the hydration of sulfuric acid with up to 15 water molecules. We find that the addition of the first four water molecules "saturate" the sulfuric acid molecule and that they are more thermodynamically favorable than the addition of water molecules 5-15. Using the large generated set of conformers, we assess the performance of approximate methods (ωB97X-D, M06-2X, PW91, and PW6B95-D3) in calculating the binding energies and assigning the global minimum conformation compared to high level CCSD(T)-F12a/VDZ-F12 reference calculations. The tested DFT functionals systematically overestimate the binding energies compared to coupled cluster calculations, and we find that this deficiency can be corrected by a simple scaling factor.

Closed-Shell Organic Compounds Might Form Dimers at the Surface of Molecular Clusters

The role of covalently bound dimer formation is studied using high-level quantum chemical methods. Reaction free energy profiles for dimer formation between common oxygen-containing functional groups are calculated, and based on the Gibbs free energy differences between transition states and reactants, we show that none of the studied two-component gas-phase reactions are kinetically feasible at 298.15 K and 1 atm. Therefore, the catalyzing effect of water, base, or acid molecules is calculated, and sulfuric acid is identified to lower the activation free energies significantly. We find that the reactions yielding hemiacetal, peroxyhemiacetal, α-hydroxyester, and geminal diol products occur with activation free energies of less than 10 kcal/mol with sulfuric acid as a catalyst, indicating that these reactions could potentially take place on the surface of sulfuric acid clusters. Additionally, the formed dimer products bind stronger onto the pre-existing cluster than the corresponding reagent monomers do. This implies that covalent dimerization reactions stabilize the existing cluster thermodynamically and make it less likely to evaporate. However, the studied small organic compounds, which contain only one functional group, are not able to form dimer products that are stable against evaporation at atmospheric conditions. Calculations of dimer formation onto a cluster surface and the clustering ability of dimer products should be extended to large terpene oxidation products in order to estimate the real atmospheric significance.

Interaction of oxalic acid with methylamine and its atmospheric implications

Oxalic acid, which is one of the most common dicarboxylic acids, is expected to be an important component of atmospheric aerosols. However, the contribution of oxalic acid to the generation of new particles is still poorly understood. In this study, the structural characteristics and thermodynamics of (C2H2O4)(CH3NH2) n (n = 1-4) were investigated at the PW91PW91/6-311++G(3df,3pd) level of theory. We found that clusters formed by oxalic acid and methylamine are relatively stable, and the more the atoms participating in the formation of a ring-like structure, the more stable is the cluster. In addition, via the analysis of atmospheric relevance, it can be revealed that clusters of (C2H2O4)(CH3NH2) n (n = 1-4) have a noteworthy concentration in the atmosphere, which indicates that these clusters could be participating in new particle formation. Moreover, by comparison with (H2C2O4)(NH3) n (n = 1-6) species, it can be seen that oxalic acid is more readily bound to methylamine than to ammonia, which promotes nucleation or new particle formation. Finally, the Rayleigh scattering properties of clusters of (C2H2O4)(CH3NH2) n (n = 1-4) were investigated for the first time to determine their atmospheric implications.

Hydration of 3-hydroxy-4,4-dimethylglutaric acid with dimethylamine complex and its atmospheric implications

Changes in temperature affects the distribution of isomers, which facilitates the understanding of new particle formation in the atmosphere.

Hydration of the methanesulfonate–ammonia/amine complex and its atmospheric implications

Methanesulfonate (MSA⁻), found in substantial concentrations in the atmosphere, is expected to enhance aerosol nucleation and the growth of nanoparticles, but the details of methanesulfonate clusters are poorly understood. In this study, MSA⁻ was chosen along with ammonia (NH₃) or three common amines and water (H₂O) to discuss the roles of ternary homogeneous nucleation and ion-induced nucleation in aerosol formation. We studied the structural characteristics and thermodynamics of the clusters using density functional theory at the PW91PW91/6-311++G(3df,3pd) level. The analysis of noncovalent interactions predicts that the amines can form more stable clusters with MSA⁻ than NH₃, in agreement with the results from structures and thermodynamics; however, the enhancement in stability for amines is not large enough to overcome the difference in the concentrations of NH₃ and amines under typical atmospheric conditions. In addition, the favorable free energies of formation for the (MSA⁻)(NH₃/amines)(H₂O)_n (n = 0–3) clusters at 298.15 K show that MSA⁻ could contribute to the aerosol nucleation process with binding NH₃/amines and H₂O up to n = 3. There are strong temperature and humidity dependences for the formation of complexes; higher humidity and temperature promote the formation of larger hydrates. Finally, for the (MSA⁻)(NH₃/amines)(H₂O)_n clusters, the evaporation rates were determined to further investigate the atmospheric implications.

Methanesulfonate (MSA⁻), found in substantial concentrations in the atmosphere, is expected to enhance aerosol nucleation and the growth of nanoparticles, but the details of methanesulfonate clusters are poorly understood.

Contribution of methane sulfonic acid to new particle formation in the atmosphere

Methane sulfonic acid (MSA) is present in substantial concentrations in the gas phase over oceans and coastal regions. We present an investigation into the contribution of MSA to new particle formation with the common atmospheric aerosol nucleation precursors including MSA, methanol, formic acid, acetone, dimethylether, formaldehyde, methyl formate, by making use a quantum chemical approach. Density functional theory calculations indicate that these bimolecular complexes are characterized by the presence of strong inter-molecular hydrogen bonds (SOH⋯O) with large binding energies and thermodynamic equilibrium constants. Topological analysis employing quantum theory of atoms in molecules shows that the charge density of the SOH⋯O hydrogen bonds of the MSA complexes falls in the range of hydrogen bonding criteria, but the Laplacian at bond critical points exceeds the range, which is due to the strong hydrogen bonding interactions. In all the studied complexes, the electrostatic interactions are found to be the main attractive force by localized molecular orbital energy decomposition analysis. All these indicate the environmental fate of MSA could play the role of nucleation centers in new particle formation. The effect of the atmospheric heights (0-12 km) was also considered. The Gibbs free energy of formation decreases with the increase of the atmospheric height owing to the decrease of the atmospheric temperature and pressure. The calculated Gibbs free energies of formation within the atmospheric temperature and pressure range could help to understand the atmospheric pollution.

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.

Particle formation and growth from oxalic acid, methanesulfonic acid, trimethylamine and water: a combined experimental and theoretical study

A combined experimental-theoretical study on the effect of oxalic acid on particle formation and growth from the reaction of MSA with trimethylamine in the absence and presence of water.

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.

Effects of Global and Local Anharmonicities on the Thermodynamic Properties of Sulfuric Acid Monohydrate

We use state-of-the-art electronic structure calculation methods and large basis sets to obtain reliable values for the thermodynamic properties of sulfuric acid monohydrate and study the effects of vibrational anharmonicity on these properties. We distinguish between two forms of vibrational anharmonicity: local anharmonicity, which refers to the anharmonicity of the vibrational modes of a given cluster conformer, and global anharmonicity, which originates from accounting for the presence of different conformers in the first place. In our most accurate approach, we solve the nuclear Schrödinger equation variationally for the intermolecular large-amplitude motions, thus quantum-mechanically accounting for the presence of higher-energy conformers for both reactants and products, while using the standard vibrational perturbational approach for the other vibrational modes. This results in a value of -11.0 kJ/mol for the reaction Gibbs free energy at 298.15 K. When standard vibrational perturbational approaches are employed, the effects of local anharmonicity depend heavily on the choice of the electronic structure calculation basis set. In fact, better results can often be achieved by combining a simple harmonic treatment for the vibrational partition function with a statistical mechanical accounting of global anharmonicity. Thus, we recommend that future studies that intend to include anharmonicity start by accounting for the presence of higher-energy conformers and only then consider whether local anharmonicity calculations are feasible and necessary.

Hydrogen bonding in cyclic complexes of carboxylic acid–sulfuric acid and their atmospheric implications

Carboxylic acids form cyclic ring structures with sulfuric acid and they could potentially be important in new particle formation.

Hydration of oxalic acid–ammonia complex: atmospheric implication and Rayleigh-scattering properties

A previous study of the binary system (H₂C₂O₄)(NH₃)_n (n = 1–6) suggested that an oxalic acid–ammonia complex may participate in atmospheric aerosol formations.

Hydration of a sulfuric acid–oxalic acid complex: acid dissociation and its atmospheric implication

We have investigated structural characteristics and thermodynamics of the hydration of a sulfuric acid–oxalic acid complex using density functional theory to gain insight into the ternary nucleation and its atmospheric implication.

On the properties and atmospheric implication of amine-hydrated clusters

Amines have been recognized as important precursor species in the formation of new atmospheric particles.