Ionization energies in solution with the QM : QM approach

We discuss a fragment-based QM:QM scheme as a practical way to access the energetics of vertical electronic processes in the condensed phase. In the QM:QM scheme, we decompose the large molecular system into small fragments, which interact solely electrostatically. The energies of the fragments are calculated in a self-consistent field generated by the other fragments and the total energy of the system is calculated as a sum of the fragment energies. We show on two test cases (cytosine and a sodium cation) that the method allows one to accurately simulate the shift of vertical ionization energies (VIE) while going from the gas phase to the bulk. For both examples, the predicted solvent shifts and peak widths estimated at the DFT level agree well with the experimental observations. We argue that the QM:QM approach is more suitable than either an electrostatic embedding based QM/MM approach, a full quantum description at the DFT level with a generally used functional or a combination of both. We also discuss the potential scope of the applicability for other electronic processes such as Auger decay.Peer reviewe

Impact of Quantum Chemistry Parameter Choices and Cluster Distribution Model Settings on Modeled Atmospheric Particle Formation Rates

We tested the influence of various parameters on the new particle formation rate predicted for the sulfuric acid–ammonia system using quantum chemistry and cluster distribution dynamics simulations, in our case, Atmospheric Cluster Dynamics Code (ACDC). We found that consistent consideration of the rotational symmetry number of monomers (sulfuric acid and ammonia molecules, and bisulfate and ammonium ions) leads to a significant rise in the predicted particle formation rate, whereas inclusion of the rotational symmetry number of the clusters only changes the results slightly, and only in conditions where charged clusters dominate the particle formation rate. This is because most of the clusters stable enough to participate in new particle formation have a rotational symmetry number of 1, and few exceptions to this rule are positively charged clusters. In contrast, the application of the quasi-harmonic correction for low-frequency vibrational modes tends to generally decrease predicted new particle formation rates and also significantly alters the slope of the formation rate curve plotted against the sulfuric acid concentration, which is a typical convention in atmospheric aerosol science. The impact of the maximum size of the clusters explicitly included in the simulations depends on the simulated conditions. The errors arising from a limited set of clusters are higher for higher evaporation rates, and thus tend to increase with temperature. Similarly, the errors tend to be higher for lower vapor concentrations. The boundary conditions for outgrowing clusters (that are counted as formed particles) have only a small influence on the results, provided that the definition is chemically reasonable and that the set of simulated clusters is sufficiently large. A comparison with data from the Cosmics Leaving OUtdoor Droplets (CLOUD) chamber and a cluster distribution dynamics model using older quantum chemistry input data shows improved agreement when using our new input data and the proposed combination of symmetry and quasi-harmonic corrections.Peer reviewe

Computational Study of the Effect of Mineral Dust on Secondary Organic Aerosol Formation by Accretion Reactions of Closed-Shell Organic Compounds

The effect of dust aerosols on accretion reactions of water, formaldehyde, and formic acid was studied in the conditions of earth's troposphere at the DLPNO-CCSD(T)/aug-cc-pVTZ//omega B97X-D/6-31++G** level of theory. A detailed analysis of the reaction mechanisms in the gas phase and on the surface of mineral dust, represented by mono- and trisilicic acid, revealed that mineral dust has the potential of decreasing reaction barrier heights. Specifically, at 0 K, mineral dust can lower the apparent energy barrier of the reaction of formaldehyde with formic acid to zero. However, when the entropic contributions to the reaction free energies were accounted for, mineral dust was found to selectively enhance the reaction of water with formaldehyde, while inhibiting the reaction of formaldehyde and formic acid, in the lower parts of the troposphere (with temperatures around 298 K). In the upper troposphere (with temperatures closer to 198 K), mineral dust catalyzes both reactions and also the reaction of methanol with formic acid. Despite the intrinsic potential of mineral dust, calculation of the catalytic enhancement parameter for a likely range of dust aerosol concentrations suggested that dust aerosols will not contribute to secondary organic aerosol formation via dimerization of closed-shell organic compounds. The main reason for this is the relatively low absolute concentratign of tropospheric dust aerosol and its inefficiency in increasing the effective reaction rate coefficients.Peer reviewe

Intersystem Crossings Drive Atmospheric Gas-Phase Dimer Formation

High molecular weight "ROOR" dimers, likely formed in the gas phase through self- and cross-reactions of complex peroxy radicals (RO2), have been suggested to play a key role in forming ultrafine aerosol particles in the atmosphere. However, the molecular-level reaction mechanism producing these dimers remains unknown. Using multireference quantum chemical methods, we explore one potentially competitive pathway for ROOR' production, involving the initial formation of triplet alkoxy radical (RO) pairs, followed by extremely rapid intersystem crossings (ISC) to the singlet surface, permitting subsequent recombination to ROOR'. Using CH3OO + CH3OO as a model system, we show that the initial steps of this reaction mechanism are likely to be very fast, as the transition states for both the formation and the decomposition of the CH3O4CH3 tetroxide intermediate are far below the reactants in energy. Next, we compute ISC rates for seven different atmospherically relevant (3)(RO center dot center dot center dot R'O) complexes. The ISC rates vary significantly depending on the conformation of the complex and also exhibit strong stereoselectivity. Furthermore, the fastest ISC process is usually not between the lowest-energy triplet and singlet states but between the triplet ground state and an exited singlet state. For each studied (RO center dot center dot center dot R'O) system, at least one low-energy conformer with an ISC rate above 10(8) s(-1) can be found. This demonstrates that gas-phase dimer formation in the atmosphere very likely involves ISCs originating in relativistic quantum mechanics.Peer reviewe

Comparing Reaction Routes for (3)(RO center dot center dot center dot OR ') Intermediates Formed in Peroxy Radical Self- and Cross-Reactions

Organic peroxy radicals (RO2) are key intermediates in the chemistry of the atmosphere. One of the main sink reactions of RO2 is the recombination reaction RO2 + R'O-2, which has three main channels (all with O-2 as a coproduct): (1) R-H=O + R'OH, (2) RO + R'O, and (3) ROOR'. The RO + R'O "alkoxy" channel promotes radical and oxidant recycling, while the ROOR' "dimer" channel leads to low-volatility products relevant to aerosol processes. The ROOR' channel has only recently been discovered to play a role in the gas phase. Recent computational studies indicate that all of these channels first go through an intermediate complex( 1)(RO center dot center dot center dot O-3(2)center dot center dot center dot OR'). Here, O-3(2) is very weakly bound and will likely evaporate from the system, giving a triplet cluster of two alkoxy radicals: (3)(RO center dot center dot center dot OR'). In this study, we systematically investigate the three reaction channels for an atmospherically representative set of RO + R'O radicals formed in the corresponding RO2+ R'O-2 reaction. First, we systematically sample the possible conformations of the RO center dot center dot center dot OR' clusters on the triplet potential energy surface. Next, we compute energetic parameters and attempt to estimate reaction rate coefficients for the three channels: evaporation/dissociation to RO + R'O, a hydrogen shift leading to the formation of R'(-H)=O + ROH, and "spin-flip" (intersystem crossing) leading to, or at least allowing, the formation of ROOR' dimers. While large uncertainties in the computed energetics prevent a quantitative comparison of reaction rates, all three channels were found to be very fast (with typical rates greater than 10 6 s(-1)). This qualitatively demonstrates that the computationally proposed novel RO2 + R'O-2 reaction mechanism is compatible with experimental data showing non-negligible branching ratios for all three channels, at least for sufficiently complex RO2.Peer reviewe

Heterogeneous Nucleation of Butanol on NaCl : A Computational Study of Temperature, Humidity, Seed Charge, and Seed Size Effects

Using a combination of quantum chemistry and cluster size distribution dynamics, we study the heterogeneous nucleation of n-butanol and water onto sodium chloride (NaCl)(10) seeds at different butanol saturation ratios and relative humidities. We also investigate how the heterogeneous nucleation of butanol is affected by the seed size through comparing (NaCl)(5), (NaCl)(10), and ( NaCl)(25) seeds and by seed electrical charge through comparing (Na10Cl9)(+), (NaCl)(10), and (Na9Cl10)(-) seeds. Butanol is a common working fluid for condensation particle counters used in atmospheric aerosol studies, and NaCl seeds are frequently used for calibration purposes and as model systems, for example, sea spray aerosol. In general, our simulations reproduce the experimentally observed trends for the NaCl-BuOH-H2O system, such as the increase of nucleation rate with relative humidity and with temperature (at constant supersaturation of butanol). Our results also provide molecular-level insights into the vapor-seed interactions driving the first steps of the heterogeneous nucleation process. The main purpose of this work is to show that theoretical studies can provide molecular understanding of initial steps of heterogeneous nucleation and that it is possible to find cost-effective yet accurate-enough combinations of methods for configurational sampling and energy evaluation to successfully model heterogeneous nucleation of multicomponent systems. In the future, we anticipate that such simulations can also be extended to chemically more complex seeds.Peer reviewe

Modeling the formation and growth of atmospheric molecular clusters : A review

Molecular clusters are ubiquitous constituents of the ambient atmosphere, that can grow into larger sizes forming new aerosol particles. The formation and growth of small clusters into aerosol particles remain one of the largest uncertainties in global climate predictions. This has made the modeling of atmospheric molecular clustering into an active field of research, yielding direct molecular level information about the formation mechanism. We review the present state of-the-art quantum chemical methods and cluster distribution dynamics models that are applied to study the formation and growth of atmospheric molecular clusters. We outline the current challenges in applying theoretical methods and the future directions to move the field forward.Peer reviewe

Hydration of Atmospheric Molecular Clusters III : Procedure for Efficient Free Energy Surface Exploration of Large Hydrated Clusters

Sampling the shallow free energy surface of hydrated atmospheric molecular clusters is a significant challenge. Using computational methods, we present an efficient approach to obtain minimum free energy structures for large hydrated clusters of atmospheric relevance. We study clusters consisting of two to four sulfuric acid (sa) molecules and hydrate them with up to five water (w) molecules. The structures of the "dry" clusters are obtained using the ABCluster program to yield a large pool of low-lying conformer minima with respect to free energy. The conformers (up to ten) lowest in free energy are then hydrated using our recently developed systematic hydrate sampling technique. Using this approach, we identify a total of 1145 unique (sa)(2-4)(w)(1-5) cluster structures. The cluster geometries and thermochemical parameters are calculated at the omega B97X-D/6-31++G(d,p) level of theory, at 298.15 K and 1 atm. The single-point energy of the most stable clusters is calculated using a high-level DLPNO-CCSD(T-0)/aug-cc-pVTZ method. Using the thermochemical data, we calculate the equilibrium hydrate distribution of the clusters under atmospheric conditions and find that the larger (sa)(3) and (sa)(4) clusters are significantly more hydrated than the smaller (sa)(2) cluster or the sulfuric acid (sa)(1) molecule. These findings indicate that more than five water molecules might be required to fully saturate the sulfuric acid clusters with water under atmospheric conditions. The presented methodology gives modelers a tool to take the effect of water explicitly into account in atmospheric particle formation models based on quantum chemistry.Peer reviewe

A study on the fragmentation of sulfuric acid and dimethylamine clusters inside an atmospheric pressure interface time-of-flight mass spectrometer

Sulfuric acid and dimethylamine vapours in the atmosphere can form molecular clusters, which participate in new particle formation events. In this work, we have produced, measured, and identified clusters of sulfuric acid and dimethylamine using an electrospray ionizer coupled with a planar-differential mobility analyser, connected to an atmospheric pressure interface time-of-flight mass spectrometer (ESI–DMA–APi-TOF MS). This set-up is suitable for evaluating the extent of fragmentation of the charged clusters inside the instrument. We evaluated the fragmentation of 11 negatively charged clusters both experimentally and using a statistical model based on quantum chemical data. The results allowed us to quantify the fragmentation of the studied clusters and to reconstruct the mass spectrum by removing the artifacts due to the fragmentation.Peer reviewe

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.Peer reviewe