Nonisothermal nucleation in the gas phase is driven by cool subcritical clusters.

Nucleation of clusters from the gas phase is a widely encountered phenomenon, yet rather little is understood about the underlying out-of-equilibrium dynamics of this process. The classical view of nucleation assumes isothermal conditions where the nucleating clusters are in thermal equilibrium with their surroundings. However, in all first-order phase transitions, latent heat is released, potentially heating the clusters and suppressing the nucleation. The question of how the released energy affects cluster temperatures during nucleation as well as the growth rate remains controversial. To investigate the nonisothermal dynamics and energetics of homogeneous nucleation, we have performed molecular dynamics simulations of a supersaturated vapor in the presence of thermalizing carrier gas. The results obtained from these simulations are compared against kinetic modeling of isothermal nucleation and classical nonisothermal theory. For the studied systems, we find that nucleation rates are suppressed by two orders of magnitude at most, despite substantial release of latent heat. Our analyses further reveal that while the temperatures of the entire cluster size populations are elevated, the temperatures of the specific clusters driving the nucleation flux evolve from cold to hot when growing from subcritical to supercritical sizes and resolve the apparent contradictions regarding cluster temperatures. Our findings provide unprecedented insight into realistic nucleation events and allow us to directly assess earlier theoretical considerations of nonisothermal nucleation.

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 (H2O) 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.

Molecular-scale evidence of aerosol particle formation via sequential addition of HIO3.

Homogeneous nucleation and subsequent cluster growth leads to the formation of new aerosol particles in the atmosphere. The nucleation of sulfuric acid and organic vapours is thought to be responsible for the formation of new particles over continents, whereas iodine oxide vapours have been implicated in particle formation over coastal regions. The molecular clustering pathways that are involved in atmospheric particle formation have been elucidated in controlled laboratory studies of chemically simple systems, but direct molecular-level observations of nucleation in atmospheric field conditions that involve sulfuric acid, organic or iodine oxide vapours have yet to be reported. Here we present field data from Mace Head, Ireland, and supporting data from northern Greenland and Queen Maud Land, Antarctica, that enable us to identify the molecular steps involved in new particle formation in an iodine-rich, coastal atmospheric environment. We find that the formation and initial growth process is almost exclusively driven by iodine oxoacids and iodine oxide vapours, with average oxygen-to-iodine ratios of 2.4 found in the clusters. On the basis of this high ratio, together with the high concentrations of iodic acid (HIO3) observed, we suggest that cluster formation primarily proceeds by sequential addition of HIO3, followed by intracluster restructuring to I2O5 and recycling of water either in the atmosphere or on dehydration. Our study provides ambient atmospheric molecular-level observations of nucleation, supporting the previously suggested role of iodine-containing species in the formation of new aerosol particles, and identifies the key nucleating compound.

Unexpected quenching effect on new particle formation from the atmospheric reaction of methanol with SO3.

Despite the high abundance in the atmosphere, alcohols in general and methanol in particular are believed to play a small role in atmospheric new particle formation (NPF) largely due to the weak binding abilities of alcohols with the major nucleation precursors, e.g., sulfuric acid (SA) and dimethylamine (DMA). Herein, we identify a catalytic reaction that was previously overlooked, namely, the reaction between methanol and SO3, catalyzed by SA, DMA, or water. We found that alcohols can have unexpected quenching effects on the NPF process, particularly in dry and highly polluted regions with high concentrations of alcohols. Specifically, the catalytic reaction between methanol and SO3 can convert methanol into a less-volatile species--methyl hydrogen sulfate (MHS). The latter was initially thought to be a good nucleation agent for NPF. However, our simulation results suggest that the formation of MHS consumes an appreciable amount of atmospheric SO3, disfavoring further reactions of SO3 with H2O. Indeed, we found that MHS formation can cause a reduction of SA concentration up to 87%, whereas the nucleation ability of MHS toward new particles is not as good as that of SA. Hence, a high abundance of methanol in the atmosphere can lower the particle nucleation rate by as much as two orders of magnitude. Such a quenching effect suggests that the recently identified catalytic reactions between alcohols and SO3 need to be considered in atmospheric modeling in order to predict SA concentration from SO2, while also account for their potentially negative effect on NPF.

Homogeneous nucleation of carbon dioxide in supersonic nozzles II: molecular dynamics simulations and properties of nucleating clusters.

Large scale molecular dynamics simulations of the homogeneous nucleation of carbon dioxide in an argon atmosphere were carried out at temperatures between 75 and 105 K. Extensive analyses of the nucleating clusters' structural and energetic properties were performed to quantify these details for the supersonic nozzle experiments described in the first part of this series [Dingilian et al., Phys. Chem. Chem. Phys., 2020, 22, 19282-19298]. We studied ten different combinations of temperature and vapour pressure, leading to nucleation rates of 1023-1025 cm-3 s-1. Nucleating clusters possess significant excess energy from monomer capture, and the observed cluster temperatures during nucleation - on both sides of the critical cluster size - are higher than that of the carrier gas. Despite strong undercooling with respect to the triple point, most clusters are clearly liquid-like during the nucleation stage. Only at the lowest simulation temperatures and vapour densities, clusters containing over 100 molecules are able to undergo a second phase transition to a crystalline solid. The formation free energies retrieved from the molecular dynamics simulations were used to improve the classical nucleation theory by introducing a Tolman-like term into the classical liquid-drop model expression for the formation free energy. This simulation-based theory predicts the simulated nucleation rates perfectly, and improves the prediction of the experimental rates compared to self-consistent classical nucleation theory.

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.

Structural Effects of Amines in Enhancing Methanesulfonic Acid-Driven New Particle Formation.

Atmospheric amines can enhance methanesulfonic acid (MSA)-driven new particle formation (NPF), but the mechanism is fundamentally different compared to that of the extensively studied sulfuric acid (SA)-driven process. Generally, the enhancing potentials of amines in SA-driven NPF follow the basicity, while this is not the case for MSA-driven NPF, where structural effects dominate, making MSA-driven NPF more prominent for methylamine (MA) compared to dimethylamine (DMA). Therefore, probing structural factors determining the enhancing potentials of amines on MSA-driven NPF is key to fully understanding the contribution of MSA to NPF. Here, we performed a comparative study on DMA and MA enhancing MSA-driven NPF by examining cluster formation using computational methods. The results indicate that DMA-MSA clusters are more stable than the corresponding MA-MSA clusters for cluster sizes up to (DMA)2(MSA)2, indicating that the basicity of amines dominates the initial cluster formation. The methyl groups of DMA were found to present significant steric hindrance beyond the (DMA)2(MSA)2 cluster and this adds to the lower hydrogen bonding capacity of DMA, making the cluster growth less favorable compared to MA. This study implies that several amines could synergistically enhance MSA-driven NPF by maximizing the advantage of different amines in different amine-MSA cluster growth stages.

Homogeneous nucleation of carbon dioxide in supersonic nozzles I: experiments and classical theories.

We studied the homogeneous nucleation of carbon dioxide in the carrier gas argon for concentrations of CO2 ranging from 2 to 39 mole percent using three experimental methods. Position-resolved pressure trace measurements (PTM) determined that the onset of nucleation occurred at temperatures between 75 and 92 K with corresponding CO2 partial pressures of 39 to 793 Pa. Small angle X-ray scattering (SAXS) measurements provided particle size distributions and aerosol number densities. Number densities of approximately 1012 cm-3, and characteristic times ranging from 6 to 13 µs, resulted in measured nucleation rates on the order of 5 × 1017 cm-3 s-1, values that are consistent with other nucleation rate measurements in supersonic nozzles. Finally, we used Fourier transform infrared (FTIR) spectroscopy to identify that the condensed CO2 particles were crystalline cubic solids with either sharp or rounded corners. Molecular dynamics simulations, however, suggest that CO2 forms liquid-like critical clusters before transitioning to the solid phase. Furthermore, the critical clusters are not in thermal equilibrium with the carrier gas. Comparisons with nucleation theories were therefore made assuming liquid-like critical clusters and incorporating non-isothermal correction factors.

Seed-Adsorbate Interactions as the Key of Heterogeneous Butanol and Diethylene Glycol Nucleation on Ammonium Bisulfate and Tetramethylammonium Bromide.

Condensation particle counter (CPC) instruments are commonly used to detect atmospheric nanoparticles. They operate on the basis of condensing an organic working fluid on the nanoparticle seeds to grow the particles to a detectable size, and at the size of few nanometers, their efficiency depends on how well the working fluid interacts with the seeds under the measurement conditions. This study models the first steps of heterogeneous nucleation of two working fluids commonly used in CPCs (diethylene glycol (DEG) and n-butanol) onto two positively charged seeds, ammonium bisulfate and tetramethylammonium bromide. The nucleation process is modeled on a molecular level using a combination of systematic configurational sampling and density functional theory (DFT). We take into account the conformational flexibility of DEG and n-butanol and determine the key factors that can improve the efficiency of nanoparticle measurements by CPCs. The results show that hydrogen bonding between the seed and the working fluid molecules is central to the adsorption of the first DEG/n-butanol molecules onto the seeds. However, intermolecular hydrogen bonding between the adsorbed molecules can also enhance the nucleation process for the weakly adsorbing vapor molecules. Accordingly, the heterogeneous nucleation probability is higher for working fluid-nanoparticle combinations with a higher potential for hydrogen bonding; in this case, DEG and ammonium bisulfate. Moreover, conformational analysis and methodology evaluations indicate that the consideration of adsorbate conformers and step-wise addition of the vapor molecules to the seeds is not essential for qualitative modeling of heterogeneous nucleation systems, at least for systems where the adsorbate and seed chemical properties are clearly different. This is the first molecular-level modeling study reporting detailed chemical reasons for experimentally observed seed and working fluid preferences in CPCs and reproducing the experimental observations. Our presented approach can be likely used for predicting preferences in similar nucleating systems.

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.

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 ω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(T0)/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.

Methanesulfonic Acid-driven New Particle Formation Enhanced by Monoethanolamine: A Computational Study.

Amines are recognized as significant enhancing species on methanesulfonic acid (MSA)-driven new particle formation (NPF). Monoethanolamine (MEA) has been detected in the atmosphere, and its concentration could be significantly increased once MEA-based postcombustion CO2 capture technology is widely implemented. Here, we evaluated the enhancing potential of MEA on MSA-driven NPF by examining the formation of MEA-MSA clusters using a combination of quantum chemical calculations and kinetics modeling. The results indicate that the -OH group of MEA can form at least one hydrogen bond with MSA or MEA in all MEA-containing clusters. The enhancing potential of MEA is higher than that of the strongest enhancing agent known so far, methylamine (MA), for MSA-driven NPF. Such high enhancing potential can be ascribed to not only the higher gas-phase basicity but also the role of the additional -OH group of MEA in increasing the binding free energy by forming additional hydrogen bonds. This clarifies the importance of hydrogen-bonding capacity from the nonamino group of amines in enhancing MSA-driven NPF. The main growth pathway for MEA-MSA clusters proceeds via the initial formation of the (MEA)1(MSA)1 cluster, followed by alternately adding one MSA and one MEA molecule, differing from the case of MA-MSA clusters.

Piperazine Enhancing Sulfuric Acid-Based New Particle Formation: Implications for the Atmospheric Fate of Piperazine.

Piperazine (PZ), a cyclic diamine, is one of 160 detected atmospheric amines and an alternative solvent to the widely used monoethanolamine in post-combustion CO2 capture. Participating in H2SO4 (sulfuric acid, SA)-based new particle formation (NPF) could be an important removal pathway for PZ. Here, we employed quantum chemical calculations and kinetics modeling to evaluate the enhancing potential of PZ on SA-based NPF by examining the formation of PZ-SA clusters. The results indicate that PZ behaves more like a monoamine in stabilizing SA and can enhance SA-based NPF at the parts per trillion (ppt) level. The enhancing potential of PZ is less than that of the chainlike diamine putrescine and greater than that of dimethylamine, which is one of the strongest enhancing agents confirmed by ambient observations and experiments. After the initial formation of the (PZ)1(SA)1 cluster, the cluster mainly grows by gradual addition of SA or PZ monomer, followed by addition of (PZ)1(SA)1 cluster. We find that the ratio of PZ removal by NPF to that by the combination of NPF and oxidations is 0.5-0.97 at 278.15 K. As a result, we conclude that participation in the NPF pathway could significantly alter the environmental impact of PZ compared to only considering oxidation pathways.

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//ω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 concentration of tropospheric dust aerosol and its inefficiency in increasing the effective reaction rate coefficients.

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.

Modeling on Fragmentation of Clusters inside a Mass Spectrometer.

Atmospheric clusters are weakly bound and can fragment inside the measuring instruments, in particular, mass spectrometers. Since the clusters accelerate under electric fields, the fragmentation cannot be described in terms of rate constants under equilibrium conditions. Using basic statistical principles, we have developed a model for fragmentation of clusters moving under an external force. The model describes an energy transfer to the cluster internal modes caused by collisions with residual carrier gas molecules. As soon as enough energy is accumulated in the cluster internal modes, it can fragment. The model can be used for interpreting experimental measurements by atmospheric pressure interface mass spectrometers.

Self-Catalytic Reaction of SO3 and NH3 To Produce Sulfamic Acid and Its Implication to Atmospheric Particle Formation.

Sulfur trioxide (SO3) is one of the most active chemical species in the atmosphere, and its atmospheric fate has profound implications to air quality and human health. The dominant gas-phase loss pathway for SO3 is generally believed to be the reaction with water molecules, resulting in sulfuric acid. The latter is viewed as a critical component in the new particle formation (NPF). Herein, a new and competitive loss pathway for SO3 in the presence of abundant gas-phase ammonia (NH3) species is identified. Specifically, the reaction between SO3 and NH3, which produces sulfamic acid, can be self-catalyzed by the reactant (NH3). In dry and heavily polluted areas with relatively high concentrations of NH3, the effective rate constant for the bimolecular SO3-NH3 reaction can be sufficiently fast through this new loss pathway for SO3 to become competitive with the conventional loss pathway for SO3 with water. Furthermore, this study shows that the final product of the reaction, namely, sulfamic acid, can enhance the fastest possible rate of NPF from sulfuric acid and dimethylamine (DMA) by about a factor of 2. An alternative source of stabilizer for acid-base clustering in the atmosphere is suggested, and this new mechanism for NPF has potential to improve atmospheric modeling in highly polluted regions.

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.

Deviation from equilibrium conditions in molecular dynamic simulations of homogeneous nucleation.

We present a comparison between Monte Carlo (MC) results for homogeneous vapour-liquid nucleation of Lennard-Jones clusters and previously published values from molecular dynamics (MD) simulations. Both the MC and MD methods sample real cluster configuration distributions. In the MD simulations, the extent of the temperature fluctuation is usually controlled with an artificial thermostat rather than with more realistic carrier gas. In this study, not only a primarily velocity scaling thermostat is considered, but also Nosé-Hoover, Berendsen, and stochastic Langevin thermostat methods are covered. The nucleation rates based on a kinetic scheme and the canonical MC calculation serve as a point of reference since they by definition describe an equilibrated system. The studied temperature range is from T = 0.3 to 0.65 ϵ/k. The kinetic scheme reproduces well the isothermal nucleation rates obtained by Wedekind et al. [J. Chem. Phys. 127, 064501 (2007)] using MD simulations with carrier gas. The nucleation rates obtained by artificially thermostatted MD simulations are consistently lower than the reference nucleation rates based on MC calculations. The discrepancy increases up to several orders of magnitude when the density of the nucleating vapour decreases. At low temperatures, the difference to the MC-based reference nucleation rates in some cases exceeds the maximal nonisothermal effect predicted by classical theory of Feder et al. [Adv. Phys. 15, 111 (1966)].

Clustering mechanism of oxocarboxylic acids involving hydration reaction: Implications for the atmospheric models.

The formation of atmospheric aerosol particles from condensable gases is a dominant source of particulate matter in the boundary layer, but the mechanism is still ambiguous. During the clustering process, precursors with different reactivities can induce various chemical reactions in addition to the formation of hydrogen bonds. However, the clustering mechanism involving chemical reactions is rarely considered in most of the nucleation process models. Oxocarboxylic acids are common compositions of secondary organic aerosol, but the role of oxocarboxylic acids in secondary organic aerosol formation is still not fully understood. In this paper, glyoxylic acid, the simplest and the most abundant atmospheric oxocarboxylic acid, has been selected as a representative example of oxocarboxylic acids in order to study the clustering mechanism involving hydration reactions using density functional theory combined with the Atmospheric Clusters Dynamic Code. The hydration reaction of glyoxylic acid can occur either in the gas phase or during the clustering process. Under atmospheric conditions, the total conversion ratio of glyoxylic acid to its hydration reaction product (2,2-dihydroxyacetic acid) in both gas phase and clusters can be up to 85%, and the product can further participate in the clustering process. The differences in cluster structures and properties induced by the hydration reaction lead to significant differences in cluster formation rates and pathways at relatively low temperatures.