Computational approaches for efficiently modelling of small atmospheric clusters

TD-DFT analysis of excitation of the closed-form spironaphthopyran in methanol solution: The contribution of vibronic transitions and intermolecular hydrogen bonds

A set of forty hybrid functionals, combined with a 6-31++G(d,p) basis set and an IEFPCM solvent description, was applied to calculate the S0→S1 excitation of the spironaphthopyran (SNP) photochromic compound. None of the hybrid functionals resulted in Cspiro-O bond cleavage upon optimization of the S1 excited state of the SNP, in contrast to the spiropyran BIPS which we studied earlier. At the same time, the deformations of the SNP during excitation turned out to be very small (in contrast to those for BIPS) which made it possible to calculate its vibronic absorption spectrum. The best agreement between the calculated and experimental absorption bands was demonstrated by the SOGGA11X functional. The main result of this work is confirmation of the vibration participation in the breaking of the Cspiro-O bond in the S1 excited state of the SNP. This assumption based on experimental data was previously put forward by other authors. Strong H-bonds between the SNP and two methanol molecules enhance the photoinduced distortions of the latter and are slightly weakened upon excitation.

Toward Modeling the Growth of Large Atmospheric Sulfuric Acid-Ammonia Clusters

Studying large atmospheric molecular clusters is needed to understand the transition between clusters and aerosol particles. In this work, we studied the (SA)n(AM)n clusters with n up to 30 and the (SA)m(AM)m±2 clusters, with m = 6-20. The cluster configurations are sampled using the ABCluster program, and the cluster geometries and thermochemical parameters are calculated using GFN1-xTB. The cluster binding energies are calculated using B97-3c. We find that the addition of sulfuric acid is preferred to the addition of ammonia. The addition free energies were found to have large uncertainties, which could potentially be attributed to errors in the applied level of theory. Based on DLPNO-CCSD(T0)/aug-cc-pVTZ benchmarks of the binding energies of the large (SA)8-9(AM)10 and (SA)10(AM)10-11 clusters, we find that ωB97X-D3BJ with a large basis set is required to yield accurate binding and addition energies. However, based on recalculations of the single-point energy at r2SCAN-3c and ωB97X-D3BJ/6-311++G(3df,3pd), we show that the single-point energy contribution is not the primary source of error. We hypothesize that a larger source of error might be present in the form of insufficient configurational sampling. Finally, we train Δ machine learning model on (SA)n(AM)n clusters with n up to 5 and show that we can predict the binding energies of clusters up to sizes of (SA)30(AM)30 with a binding energy error below 0.6 %. This is an encouraging approach for accurately modeling the binding energies of large acid-base clusters in the future.

Nitric Acid and Organic Acids Suppress the Role of Methanesulfonic Acid in Atmospheric New Particle Formation

Multicomponent atmospheric molecular clusters, typically comprising a combination of acids and bases, play a pivotal role in our climate system and contribute to the perplexing uncertainties embedded in modern climate models. Our understanding of cluster formation is limited by the lack of studies on complex mixed-acid-mixed-base systems. Here, we investigate multicomponent clusters consisting of mixtures of several acid and base molecules: sulfuric acid (SA), methanesulfonic acid (MSA), nitric acid (NA), formic acid (FA), along with methylamine (MA), dimethylamine (DMA), and trimethylamine (TMA). We calculated the binding free energies of a comprehensive set of 252 mixed-acid-mixed-base clusters at the DLPNO-CCSD(T0)/aug-cc-pVTZ//ωB97X-D/6-31++G(d,p) level of theory. Combined with the existing datasets, we simulated the new particle formation (NPF) rates using the Atmospheric Cluster Dynamics Code (ACDC). We find that the presence of NA and FA had a substantial impact, increasing the NPF rate by 60% at realistic conditions. Intriguingly, we find that NA and FA suppress the role of MSA in NPF. These findings suggest that even high concentration of MSA has a limited impact on NPF in polluted regions with high FA and NA. We outline a method for generating a lookup table that could potentially be used in climate models that sufficiently incorporates all the required chemistry. By unraveling the molecular mechanisms of mixed-acid-mixed-base clusters, we get one step closer to comprehending their implications for our global climate system.

Limited Role of Malonic Acid in Sulfuric Acid-Dimethylamine New Particle Formation

Aerosols play an important role in climate and air quality; however, the mechanisms behind aerosol particle formation in the atmosphere are poorly understood. Studies have identified sulfuric acid, water, oxidized organics, and ammonia/amines as key precursors for forming aerosol particles in the atmosphere. Theoretical and experimental investigations have indicated that other species, such as organic acids, may be involved in atmospheric nucleation and growth of freshly formed aerosol particles. Organic acids, such as dicarboxylic acids, which are abundant in the atmosphere, have been measured in ultrafine aerosol particles. These observations suggest that organic acids may contribute to new particle formation in the atmosphere but their role remains ambiguous. This study examines how malonic acid interacts with sulfuric acid and dimethylamine to form new particles at warm boundary layer conditions using experimental observations from a laminar flow reactor and quantum chemical calculations coupled with cluster dynamics simulations. Observations reveal that malonic acid does not contribute to the initial steps (formation of <1 nm diameter particle) of nucleation with sulfuric acid-dimethylamine. In addition, malonic acid was found to not participate in the subsequent growth of the freshly nucleated 1 nm particles from sulfuric acid-dimethylamine reactions to diameters of 2 nm.

Transition Metal Sensing with Nitrogenated Holey Graphene: A First-Principles Investigation

The toxicity of transition metals, including copper(II), manganese(II), iron(II), zinc(II), hexavalent chromium, and cobalt(II), at elevated concentrations presents a significant threat to living organisms. Thus, the development of efficient sensors capable of detecting these metals is of utmost importance. This study explores the utilization of two-dimensional nitrogenated holey graphene (C₂N) nanosheet as a sensor for toxic transition metals. The C₂N nanosheet's periodic shape and standard pore size render it well suited for adsorbing transition metals. The interaction energies between transition metals and C₂N nanosheets were calculated in both gas and solvent phases and were found to primarily result from physisorption, except for manganese and iron which exhibited chemisorption. To assess the interactions, we employed NCI, SAPT0, and QTAIM analyses, as well as FMO and NBO analysis, to examine the electronic properties of the TM@C₂N system. Our results indicated that the adsorption of copper and chromium significantly reduced the HOMO-LUMO energy gap of C₂N and significantly increased its electrical conductivity, confirming the high sensitivity of C₂N towards copper and chromium. The sensitivity test further confirmed the superior sensitivity and selectivity of C₂N towards copper. These findings offer valuable insight into the design and development of sensors for the detection of toxic transition metals.

Heterogeneous Ion-Induced Nucleation of Water and Butanol Vapors Studied via Computational Quantum Chemistry beyond Prenucleation and Critical Cluster Sizes

Water and butanol are used as working fluids in condensation particle counters, and condensation of a single vapor onto an ion can be used as a simple model system for the study of ion-induced nucleation in the atmosphere. Motivated by this, we examine heterogeneous nucleation of water (H₂O) and n-butanol (BuOH) vapors onto three positively (Li⁺, Na⁺, K⁺) and three negatively charged (F^-, Cl^-, Br^-) ions using classical nucleation theory and computational quantum chemistry methods. We study phenomena that cannot be captured by Kelvin-Thomson equation for small nucleation ion cores. Our quantum chemistry calculations reveal the molecular mechanism behind ion-induced nucleation for each studied system. Typically, ions become solvated from all sides after several vapor molecules condense onto the ion. However, we show that the clusters of water and large negatively charged ions (Cl^- and Br^-) thermodynamically prefer the ion being migrated to the cluster surface. Although our methods generally do not show clear sign-preference for ion-water nucleation, we identified positive sign-preference for ion-butanol nucleation caused by the possibility to form stabilizing hydrogen bonds between butanol molecules condensed onto a positively charged ion. These bonds cannot form when butanol condenses onto a negatively charged ion. Therefore, we show that ion charge, its sign, as well as vapor properties have effects on the prenucleation and critical cluster/droplet sizes and also on the molecular mechanism of ion-induced nucleation.

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.

Clusteromics V: Organic Enhanced Atmospheric Cluster Formation

Formic acid (FA) is a prominent candidate for organic enhanced nucleation due to its high abundance and stabilizing effect on smaller clusters. Its role in new particle formation is studied through the use of state-of-the-art quantum chemical methods on the cluster systems (acid)_1-2(FA)₁(base)_1-2 with the acids being sulfuric acid (SA)/methanesulfonic acid (MSA) and the bases consisting of ammonia (A), methylamine (MA), dimethylamine (DMA), trimethylamine (TMA), and ethylenediamine (EDA). A funneling approach is used to determine the cluster structures with initial configurations generated through the ABCluster program, followed by semiempirical PM7 and ωB97X-D/6-31++G(d,p) calculations. The final binding free energy is calculated at the DLPNO-CCSD(T₀)/aug-cc-pVTZ//ωB97X-D/6-31++G(d,p) level of theory using the quasi-harmonic approximation. Cluster dynamics simulations show that FA has a minuscule or negligible effect on the MSA-FA-base systems as well as most of the SA-FA-base systems. The SA-FA-DMA cluster system shows the highest influence from FA with an enhancement of 21%, compared to its non-FA counterpart.

Massive Assessment of the Binding Energies of Atmospheric Molecular Clusters

Quantum chemical studies of the formation and growth of atmospheric molecular clusters are important for understanding aerosol particle formation. However, the search for the lowest free-energy cluster configuration is extremely time consuming. This makes high-level benchmark data sets extremely valuable in the quest for the global minimum as it allows the identification of cost-efficient computational methodologies, as well as the development of high-level machine learning (ML) models. Herein, we present a highly versatile quantum chemical data set comprising a total of 11 749 (acid)_1-2(base)_1-2 cluster configurations, containing up to 44 atoms. Utilizing the LUMI supercomputer, we calculated highly accurate PNO-CCSD(F12*)(T)/cc-pVDZ-F12 binding energies of the full set of cluster configurations leading to an unprecedented data set both in regard to sheer size and with respect to the level of theory. We employ the constructed benchmark set to assess the performance of various semiempirical and density functional theory methods. In particular, we find that the r²-SCAN-3c method shows excellent performance across the data set related to both accuracy and CPU time, making it a promising method to employ during cluster configurational sampling. Furthermore, applying the data sets, we construct ML models based on Δ-learning and provide recommendations for future application of ML in cluster configurational sampling.

The driving effects of common atmospheric molecules for formation of prenucleation clusters: the case of sulfuric acid, formic acid, nitric acid, ammonia, and dimethyl amine

How secondary aerosols form is critical as aerosols' impact on Earth's climate is one of the main sources of uncertainty for understanding global warming. The beginning stages for formation of prenucleation complexes, that lead to larger aerosols, are difficult to decipher experimentally. We present a computational chemistry study of the interactions between three different acid molecules and two different bases. By combining a comprehensive search routine covering many thousands of configurations at the semiempirical level with high level quantum chemical calculations of approximately 1000 clusters for every possible combination of clusters containing a sulfuric acid molecule, a formic acid molecule, a nitric acid molecule, an ammonia molecule, a dimethylamine molecule, and 0-5 water molecules, we have completed an exhaustive search of the DLPNO-CCSD(T)/CBS//ωB97X-D/6-31++G** Gibbs free energy surface for this system. We find that the detailed geometries of each minimum free energy cluster are often more important than traditional acid or base strength. Addition of a water molecule to a dry cluster can enhance stabilization, and we find that the (SA)(NA)(A)(DMA)(W) cluster has special stability. Equilibrium calculations of SA, FA, NA, A, DMA, and water using our quantum chemical ΔG° values for cluster formation and realistic estimates of the concentrations of these monomers in the atmosphere reveals that nitric acid can drive early stages of particle formation just as efficiently as sulfuric acid. Our results lead us to believe that particle formation in the atmosphere results from the combination of many different molecules that are able to form highly stable complexes with acid molecules such as SA, NA, and FA.

Contribution of Methanesulfonic Acid to the Formation of Molecular Clusters in the Marine Atmosphere

Because of the lack of long-term measurements, new particle formation (NPF) in the marine atmosphere remains puzzling. Using quantum chemical methods, this study elucidates the cluster formation and further growth of sulfuric acid-methanesulfonic acid-dimethylamine (SA-MSA-DMA) clusters, relevant to NPF in the marine atmosphere. The cluster structures and thermochemical parameters of (SA)_n(MSA)_m(DMA)_l (n + m ≤ 4 and l ≤ 4) systems are calculated using density functional theory at the ωB97X-D/6-31++G(d,p) level of theory, and the single-point energies are calculated using high-level DLPNO-CCSD(T₀)/aug-cc-pVTZ calculations. The calculated thermochemistry is used as input to the Atmospheric Cluster Dynamics Code (ACDC) to gain insight into the cluster dynamics. At ambient conditions (298.15 K, 1 atm), we find that the distribution of outgrowing clusters primarily consists of SA and DMA, with a minor contribution from the mixed SA-MSA-DMA clusters. At lower temperature (278.15 K, 1 atm) the distribution broadens, and clusters containing one or more MSA molecules emerge. These findings show that in the cold marine atmosphere MSA likely participates in atmospheric NPF.

Clusteromics IV: The Role of Nitric Acid in Atmospheric Cluster Formation

Nitric acid (NA) has previously been shown to affect atmospheric new particle formation; however, its role still remains highly uncertain. Through the employment of state-of-the-art quantum chemical methods, we study the (acid)_1-2(base)_1-2 and (acid)₃(base)₂ clusters containing at least one nitric acid (NA) and sulfuric acid (SA) or methanesulfonic acid (MSA) with bases ammonia (A), methylamine (MA), dimethylamine (DMA), trimethylamine (TMA), and ethylenediamine (EDA). The initial cluster configurations are generated using the ABCluster program. PM7 and ωB97X-D/6-31++G(d,p) calculations are used to reduce the number of relevant configurations. The thermochemical parameters are calculated at the ωB97X-D/6-31++G(d,p) level of theory with the quasi-harmonic approximation, and the final single-point energies are calculated with high-level DLPNO-CCSD(T₀)/aug-cc-pVTZ calculations. The enhancing effect from the presence of nitric acid on cluster formation is studied using the calculated thermochemical data and cluster dynamics simulations. We find that when NA is in excess compared with the other acids, it has a substantial enhancing effect on the cluster formation potential.

Clusteromics III: Acid Synergy in Sulfuric Acid-Methanesulfonic Acid-Base Cluster Formation

Acid-base molecular clusters are an important stage in atmospheric new particle formation. While such clusters are most likely multicomponent in nature, there are very few reports on clusters consisting of multiple acid molecules and multiple base molecules. By applying state-of-the-art quantum chemical methods, we herein study electrically neutral (SA)₁(MSA)₁(base)_0-2 clusters with base = ammonia (A), methylamine (MA), dimethylamine (DMA), trimethylamine (TMA) and ethylenediamine (EDA). The cluster structures are obtained using a funneling approach employing the ABCluster program, semiempirical PM7 calculations and ωB97X-D/6-31++G(d,p) calculations. The final binding free energies are calculated at the DLPNO-CCSD(T₀)/aug-cc-pVTZ//ωB97X-D/6-31++G(d,p) level of theory using the quasi-harmonic approximation. Based on the calculated cluster geometries and thermochemistry (at 298.15 K and 1 atm), we find that the mixed (SA)₁(MSA)₁(base)_1-2 clusters more resemble the (SA)₂(base)_1-2 clusters compared to the (MSA)₂(base)_1-2 clusters. Hence, some of the steric hindrance and lack of hydrogen bond capacity previously observed in the (MSA)₂(base)_1-2 clusters is diminished in the corresponding (SA)₁(MSA)₁(base)_1-2 clusters. Cluster kinetics simulations reveal that the presence of an MSA molecule in the clusters enhances the cluster formation potential by up to a factor of 20. We find that the SA-MSA-DMA clusters have the highest cluster formation potential, and thus, this system should be further extended to larger sizes in future studies.

Hydrogen-Bond Topology Is More Important Than Acid/Base Strength in Atmospheric Prenucleation Clusters

We explored the hypothesis that on the nanoscale level, acids and bases might exhibit different behavior than in bulk solution. Our study system consisted of sulfuric acid, formic acid, ammonia, and water. We calculated highly accurate Domain-based Local pair-Natural Orbital- Coupled-Cluster/Complete Basis Set (DLPNO-CCSD(T)/CBS) energies on DFT geometries and used the resulting Gibbs free energies for cluster formation to compute the overall equilibrium constants for every possible cluster. The equilibrium constants combined with the initial monomer concentrations were used to predict the formation of clusters at the top and the bottom of the troposphere. Our results show that formic acid is as effective as ammonia at forming clusters with sulfuric acid and water. The structure of formic acid is uniquely suited to form hydrogen bonds with sulfuric acid. Additionally, it can partner with water to form bridges from one side of sulfuric acid to the other, hence demonstrating that hydrogen bonding topology is more important than acid/base strength in these atmospheric prenucleation clusters.

Monomers of Glycine and Serine Have a Limited Ability to Hydrate in the Atmosphere

The role of atmospheric aerosols on climate change is one of the biggest uncertainties in most global climate models. Organic aerosols have been identified as potential cloud condensation nuclei (CCN), and amino acids are organic molecules that could serve as CCN. Amino acids make up a significant portion of the total organic material in the atmosphere, and herein we present a systematic study of hydration for two of the most common atmospheric amino acids, glycine and serine. We compute DLPNO/CCSD(T)//M08-HX/MG3S thermodynamic properties and atmospheric concentrations of Gly(H₂O)_n and Ser(H₂O)_n, where n = 1-5. We predict that serine-water clusters have higher concentrations at n = 1 and 5, while glycine-water clusters have higher concentrations at n = 2-4. However, both glycine and serine are inferred to exist primarily in their nonhydrated monomer forms in the absence of other species such as sulfuric acid.

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.

New Particle Formation and Growth from Dimethyl Sulfide Oxidation by Hydroxyl Radicals

Dimethyl sulfide (DMS) is produced by plankton in oceans and constitutes the largest natural emission of sulfur to the atmosphere. In this work, we examine new particle formation from the primary pathway of oxidation of gas-phase DMS by OH radicals. We particularly focus on particle growth and mass yield as studied experimentally under dry conditions using the atmospheric simulation chamber AURA. Experimentally, we show that aerosol mass yields from oxidation of 50-200 ppb of DMS are low (2-7%) and that particle growth rates (8.2-24.4 nm/h) are comparable with ambient observations. An HR-ToF-AMS was calibrated using methanesulfonic acid (MSA) to account for fragments distributed across both the organic and sulfate fragmentation table. AMS-derived chemical compositions revealed that MSA was always more dominant than sulfate in the secondary aerosols formed. Modeling using the Aerosol Dynamics, gas- and particle-phase chemistry kinetic multilayer model for laboratory CHAMber studies (ADCHAM) indicates that the Master Chemical Mechanism gas-phase chemistry alone underestimates experimentally observed particle formation and that DMS multiphase and autoxidation chemistry is needed to explain observations. Based on quantum chemical calculations, we conclude that particle formation from DMS oxidation in the ambient atmosphere will most likely be driven by mixed sulfuric acid/MSA clusters clustering with both amines and ammonia.

Assessment of the DLPNO Binding Energies of Strongly Noncovalent Bonded Atmospheric Molecular Clusters

This work assesses the performance of DLPNO-CCSD(T0), DLPNO-MP2, and density functional theory methods in calculating the binding energies of a representative test set of 45 atmospheric acid-acid, acid-base, and acid-water dimer clusters. The performance of the approximate methods is compared to high level explicitly correlated CCSD(F12*)(T)/complete basis set (CBS) reference calculations. Out of the tested density functionals, ωB97X-D3(BJ) shows the best performance with a mean deviation of 0.09 kcal/mol and a maximum deviation of 0.83 kcal/mol. The RI-CC2/aug-cc-pV(T+d)Z level of theory severely overpredicts the cluster binding energies with a mean deviation of -1.31 kcal/mol and a maximum deviation up to -3.00 kcal/mol. Hence, RI-CC2/aug-cc-pV(T+d)Z should not be utilized for studying atmospheric molecular clusters. The DLPNO variants are tested both with and without the inclusion of explicit correlation (F12) in the wavefunction, with different pair natural orbital (PNO) settings (loosePNO, normalPNO, and tightPNO) and using both double and triple zeta basis sets. The performance of the DLPNO-MP2 methods is found to be independent of PNO settings and yield low mean deviations of -0.84 kcal/mol or below. However, DLPNO-MP2 requires explicitly correlated wavefunctions to yield maximum deviations below 1.40 kcal/mol. For obtaining high accuracy, with maximum deviation below ∼1.0 kcal/mol, either DLPNO-CCSD(T0)/aug-cc-pVTZ (normalPNO) calculations or DLPNO-CCSD(T0)-F12/cc-pVTZ-F12 (normalPNO) calculations are required. The most accurate level of theory is found to be DLPNO-CCSD(T0)-F12/cc-pVTZ-F12 using a tightPNO criterion which yields a mean deviation of 0.10 kcal/mol, with a maximum deviation of 0.20 kcal/mol, compared to the CCSD(F12*)(T)/CBS reference.

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.

Strong Even/Odd Pattern in the Computed Gas-Phase Stability of Dicarboxylic Acid Dimers: Implications for Condensation Thermodynamics

The physical properties of small straight-chain dicarboxylic acids are well known to exhibit even/odd alternations with respect to the carbon chain length. For example, odd numbered diacids have lower melting points and higher saturation vapor pressures than adjacent even numbered diacids. This alternation has previously been explained in terms of solid-state properties, such as higher torsional strain of odd number diacids. Using quantum chemical methods, we demonstrate an additional contribution to this alternation in properties resulting from gas-phase dimer formation. Due to a combination of hydrogen bond strength and torsional strain, dimer formation in the gas phase occurs efficiently for glutaric acid (C5) and pimelic acid (C7) but is unfavorable for succinic acid (C4) and adipic acid (C6). Our results indicate that a significant fraction of the total atmospheric gas-phase concentration of glutaric and pimelic acid may consist of dimers.

Configurational Sampling of Noncovalent (Atmospheric) Molecular Clusters: Sulfuric Acid and Guanidine

We studied the configurational sampling of noncovalently bonded molecular clusters relevant to the atmosphere. In this article, we discuss possible approaches to searching for optimal configurations and present one alternative based on systematic configurational sampling, which seems able to overcome the typical problems associated with searching for global minima on multidimensional potential energy surfaces. Since atmospheric molecular clusters are usually held together by intermolecular bonds, we also present a cost-effective strategy for treating hydrogen bonding and proton transferring by using rigid molecules and ions in different protonation states and illustrate its performance on clusters containing guanidine and sulfuric acid.

Benchmarking sampling methodology for calculations of Rayleigh light scattering properties of atmospheric molecular clusters

We present four different computational methods for benchmarking the sampling and Rayleigh light scattering of hydrogen bonded atmospheric molecular clusters.

Guanidine: A Highly Efficient Stabilizer in Atmospheric New-Particle Formation

The role of a strong organobase, guanidine, in sulfuric acid-driven new-particle formation is studied using state-of-the-art quantum chemical methods and molecular cluster formation simulations. Cluster formation mechanisms at the molecular level are resolved, and theoretical results on cluster stability are confirmed with mass spectrometer measurements. New-particle formation from guanidine and sulfuric acid molecules occurs without thermodynamic barriers under studied conditions, and clusters are growing close to a 1:1 composition of acid and base. Evaporation rates of the most stable clusters are extremely low, which can be explained by the proton transfers and symmetrical cluster structures. We compare the ability of guanidine and dimethylamine to enhance sulfuric acid-driven particle formation and show that more than 2000-fold concentration of dimethylamine is needed to yield as efficient particle formation as in the case of guanidine. At similar conditions, guanidine yields 8 orders of magnitude higher particle formation rates compared to dimethylamine. Highly basic compounds such as guanidine may explain experimentally observed particle formation events at low precursor vapor concentrations, whereas less basic and more abundant bases such as ammonia and amines are likely to explain measurements at high concentrations.

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.

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.

Ion-induced nucleation of pure biogenic particles

Atmospheric aerosols and their effect on clouds are thought to be important for anthropogenic radiative forcing of the climate, yet remain poorly understood. Globally, around half of cloud condensation nuclei originate from nucleation of atmospheric vapours. It is thought that sulfuric acid is essential to initiate most particle formation in the atmosphere, and that ions have a relatively minor role. Some laboratory studies, however, have reported organic particle formation without the intentional addition of sulfuric acid, although contamination could not be excluded. Here we present evidence for the formation of aerosol particles from highly oxidized biogenic vapours in the absence of sulfuric acid in a large chamber under atmospheric conditions. The highly oxygenated molecules (HOMs) are produced by ozonolysis of α-pinene. We find that ions from Galactic cosmic rays increase the nucleation rate by one to two orders of magnitude compared with neutral nucleation. Our experimental findings are supported by quantum chemical calculations of the cluster binding energies of representative HOMs. Ion-induced nucleation of pure organic particles constitutes a potentially widespread source of aerosol particles in terrestrial environments with low sulfuric acid pollution.

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.

Characterization of the nucleation precursor (H2SO4–(CH3)2NH) complex: intra-cluster interactions and atmospheric relevance

Amines have been proposed to participate in the nucleation process, but the electron density analysis and the determination of a temperature dependence of the clusters are still lacking.