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

Stability and Structure of Potentially Atmospherically Relevant Glycine Ammonium Bisulfate Clusters

New particle formation (NPF) is the process by which trace atmospheric acids and bases cluster and grow into particles that ultimately impact climate. Sulfuric acid concentration drives NPF, but nitrogen-containing bases promote the formation of more stable clusters via salt bridge formation. Recent computational efforts have suggested that amino acids can enhance NPF, predicting that they can stabilize new particles via multiple protonation sites, but there has yet to be experimental validation of these predictions. We used mass spectrometry and infrared spectroscopy to study the structure and stability of cationic clusters composed of glycine, sulfuric acid, and ammonia. When collisionally activated, clusters were significantly more likely to eliminate ammonia or sulfuric acid than glycine, while quantum chemical calculations predicted lower binding free energies for ammonia but similar binding free energies for glycine and sulfuric acid. These calculations predicted several low-energy structures, so we compared experimental and computed vibrational spectra to attempt to validate the computationally predicted minimum energy structure. Unambiguous identification of the experimental structure by comparison to these calculations was made difficult by the complexity of the experimental spectra and the fact that the identity of the computed lowest-energy structure depended strongly on temperature. If their vapors are present, amino acids are likely to be enriched in new particles by displacing more weakly bound ammonia, similar to the behavior of other atmospheric amines. The carboxylic acid groups were found to preferentially interact with other carboxylic acids, suggesting incipient organic/inorganic phase separation even at these small sizes.

Revealing the Impacts of Chemical Complexity on Submicrometer Sea Spray Aerosol Morphology

Sea spray aerosol (SSA) ejected through bursting bubbles at the ocean surface is a complex mixture of salts and organic species. Submicrometer SSA particles have long atmospheric lifetimes and play a critical role in the climate system. Composition impacts their ability to form marine clouds, yet their cloud-forming potential is difficult to study due to their small size. Here, we use large-scale molecular dynamics (MD) simulations as a "computational microscope" to provide never-before-seen views of 40 nm model aerosol particles and their molecular morphologies. We investigate how increasing chemical complexity impacts the distribution of organic material throughout individual particles for a range of organic constituents with varying chemical properties. Our simulations show that common organic marine surfactants readily partition between both the surface and interior of the aerosol, indicating that nascent SSA may be more heterogeneous than traditional morphological models suggest. We support our computational observations of SSA surface heterogeneity with Brewster angle microscopy on model interfaces. These observations indicate that increased chemical complexity in submicrometer SSA leads to a reduced surface coverage by marine organics, which may facilitate water uptake in the atmosphere. Our work thus establishes large-scale MD simulations as a novel technique for interrogating aerosols at the single-particle level.

Roles of Amides on the Formation of Atmospheric HONO and the Nucleation of Nitric Acid Hydrates

Nitrous acid (HONO) is hazardous to the human respiratory system, and the hydrolysis of NO₂ is the source of HONO. Hence, the investigation on the removal and transformation of HONO is urgently established. The effects of amide on the mechanism and kinetics of the formation of HONO with acetamide, formamide, methylformamide, urea, and its clusters of the catalyst were studied theoretically. The results show that amide and its small clusters reduce the energy barrier, the substituent improves the catalytic efficiency, and the catalytic effect order is dimer > monohydrate > monomer. Meanwhile, the clusters composed of nitric acid (HNO₃), amides, and 1-6 water molecules were investigated in the amide-assisted nitrogen dioxide (NO₂) hydrolysis reaction after HONO decomposes by combining the system sampling technique and density functional theory. The study on thermodynamics, intermolecular forces, optics properties of the clusters, as well as the influence of humidity, temperature, atmospheric pressure, and altitude shows that amide molecules promote the clustering and enhance the optical properties. The substituent facilitates the clustering of amide and nitric acid hydrate and lowers the humidity sensitivity of the clusters. The findings will help to control the atmospheric aerosol particle and then reduce the harm of poisonous organic chemicals on human health.

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.

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.

Amino Acids Compete with Ammonia in Sulfuric Acid-Based Atmospheric Aerosol Prenucleation: The Case of Glycine and Serine

We present a computational investigation of the sulfuric acid, glycine, serine, ammonia, and water system to understand if this system can form prenucleation clusters, which are precursors to larger aerosols in the atmosphere. We have performed a comprehensive configurational search of all possible clusters in this system, starting with the four different monomers and zero to five waters. Accurate Gibbs free energies of formation have been calculated with the DLPNO-CCSD(T)/complete basis set (CBS) method on ωb97xd/6-31++G** geometries. For the dry dimers of sulfuric acid, the weakest base, serine, is found to form the most stable complex, which is a consequence of the strong di-ionic complex formed between the bisulfate ion and the protonated serine cation. For the dry dimers without sulfuric acid, the glycine-serine complex is more stable than the glycine-ammonia or serine-ammonia complexes, stemming from the detailed structure and not related to base strength. For the larger complexes, sulfuric acid deprotonates and the proton is shifted to glycine, serine, or ammonia. The two amino acids and ammonia are almost interchangeable and there is no easy way to predict which molecule will be protonated without the calculated results. Assuming reasonable starting concentrations and a closed system of sulfuric acid, glycine, serine, ammonia, and five waters, we predict the concentrations of all possible complexes at two temperatures spanning the troposphere. The most negative ΔG° values are a function of the detailed molecular interactions of these clusters. These details are more important than the base strength of ammonia, glycine, and serine.

Gas-phase and aqueous-surface reaction mechanism of Criegee radicals with serine and nucleation of products: A theoretical study

Criegee intermediates (CIs) are short-lived carbonyl oxides, which can affect the budget of OH radicals, ozone, ammonia, organic/inorganic acids in the troposphere. This study investigated the reaction of CIs with serine (Ser) in the gas phase by using density functional theory (DFT) calculations and at the gas-liquid interface by using Born-Oppenheimer molecular dynamics (BOMD). The results reveal that the reactivity of the three functional groups of Ser can be ordered as follows: COOH > NH₂ > OH. Water-mediated reactions of CIs with NH₂ and OH groups of Ser on the droplet follow the proton exchange mechanism. The products, sulfuric acids, ammonia, and water molecules form stable clusters within 20 ns. This study shows that hydroperoxide products can contribute to new particle formation (NPF). The result deepens the understanding of the reaction of CIs with multifunctional pollutants and atmospheric behavior of CIs in polluted areas.

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.

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.

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

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

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

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

Enhancement of Atmospheric Nucleation by Highly Oxygenated Organic Molecules: A Density Functional Theory Study

New particle formation (NPF) by gas-particle conversion is the main source of atmospheric aerosols. Highly oxygenated organic molecules (HOMs) and sulfuric acid (SA) are important NPF participants. 2-Methylglyceric acid (MGA), a kind of HOMs, is a tracer of isoprene-derived secondary organic aerosols. The nucleation mechanisms of MGA with SA were studied using density functional theory and atmospheric cluster dynamics simulation in this study, along with that of MGA with methanesulfonic acid (MSA) as a comparison. Our theoretical works indicate that the (MGA)(SA) and (MGA)(MSA) clusters are the most stable ones in the (MGA) _i(SA) _j ( i = 1-2, j = 1-2) and (MGA) _i(MSA) _j ( i = 1-2, j = 1-2) clusters, respectively. Both the formation rates of (MGA)(SA) and (MGA)(MSA) clusters are quite large and could have significant contributions to NPF. The results imply that the homomolecular nucleation of MGA is unlikely to occur in the atmosphere, and MGA and SA can effectively contribute to heteromolecular nucleation mainly in the form of heterodimers. MSA exhibits properties similar to SA in its ability to form clusters with MGA but is slightly weaker than SA.

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.