Introduction

Atmospheric aerosols affect air quality, human health, and the global climate1,2,3. Aerosols act as cloud condensation nuclei (CCN) and ice nuclei, and impact frequency of occurrence, and lifetime of clouds4,5,6,7. Despite their broad impact, the exact formation pathways for the atmospheric particles remain largely unknown8. Sulfuric acid (H2SO4) is recognized as perhaps the most important nucleating species in the atmosphere9,10,11. However, binary homogenous nucleation of H2SO4 and water (H2O) is insufficient to explain the nucleation events under actual atmospheric conditions12, suggesting that other species may also participate in nucleation events13,14,15.

Amines and carboxylic acids are an important class of compounds that have recently been found to contribute to the particle formation under certain conditions8,16. It is noted that molecules with high vapor pressures, which includes certain carboxylic acids, would be unlikely to contribute to new particle formation. However, recent field studies from a variety of locations around the world have found evidence of carboxylic acids contributing to new particle formation. For example, the mass spectra of nanoparticles obtained during new particle formation events in Hyytiala, Finland is suggestive of the presence of alkylaminium and carboxylate ions in the particles17. Organic species are also predicted to play an important role in the particle growth in Tecamac, Mexico13. Zhang et al. observed a marked increase in the particle concentration when benzoic, p-toluic, m-toluic, or cis-pinonic acid vapor was added to the H2SO4–H2O vapor system18,19. Additional measurements in Shanghai, China observed increased rates of particle formation containing markers for the sulfuric acid dimer, H2SO4-dimethylamine (DMA) clusters and H2SO4-DMA-H2O nucleation20. The role of oxalic acid in particle formation from vapor phase methanesulfonic acid, methylamine, and H2O has been recently examined experimentally21. The addition of water to the mixture of oxalic acid and methylamine enhanced the rate of particle formation by an order of magnitude. Although the amine concentrations are 1–3 orders of magnitude lower than that of NH3 in the atmosphere22, amines are experimentally found to be more effective than NH3 in enhancing particle formation11,23,24,25,26. For example, the CLOUD chamber experiments at CERN have demonstrated the superior particle forming efficacy of DMA compared to NH311. The vapor pressure of monocarboxylic acids is higher, by a factor of 102 to 104 than that of the corresponding dicarboxylic acids27, suggesting that the latter is likely to play a bigger role in the new particle formation. Despite that, monocarboxylic acids have been detected in aerosol particles collected in various regions of the world28,29,30,31,32,33,34. Supplementary Table 1 summarizes the concentration of HCOOH and CH3COOH measured in air samples collected in a variety of environments. In the urban environment of New Mexico City, HCOOH and CH3COOH were predominantly detected in the particulate phase35. On average, 53% of HCOOH and 67% of CH3COOH were present in particulate matter. The presence of particulate organic salts was also raised as a possible explanation for the relatively high aerosol hygroscopicity and CCN concentrations observed in the Amazon basin36.

Despite being ubiquitous in diverse environments, the mechanism of organic particulate formation is yet to be fully established. Specifically, identifying key organic species and the underlying chemical mechanisms responsible for the nucleation and growth of atmospheric particles remains a significant challenge. Our current understanding about the role of organics in initial atmospheric nucleation events is based on quantum chemical calculations8,37,38,39,40,41,42, which are typically limited to few atom clusters. This hampers the generalization of gas-phase results to the water surfaces, such as fog, snow, clouds and water microdroplets, which are present in the troposphere and are believed to impact the chemistries occurring there43,44,45,46,47,48,49.

We have performed experiments and BOMD simulations to explore the particle formation from HCOOH, TMA ((CH3)3N), and water. The BOMD simulations and performed experiments provide a further understanding as to the mechanism behind high vapor pressure molecules contributing to new particle formation. The results of this experiment-theory study also help in understanding the role of organics in both haze and cloud droplet formation, and to nanoparticle growth in urban, rural and remote regions.

Results

Experimental results of particle formation enhanced by trimethylamine

Experiments were conducted by flowing gas mixtures of formic acid, water vapor, and TMA in a flow reaction cell. The instrumental set-up is shown in Supplementary Fig. 1. Three sets of conditions were analyzed for each precursor, which resulted in 27 conditions (Supplementary Table 2), each analyzed at 6 different reaction times (8, 16, 24, 32, 40, and 48 s). The analysis of the 8 minimal and maximal concentrations are discussed in the following figures. However, it should be noted that the additional 19 conditions agree with the minimal and maximal results. The size distribution and total number of particles at each reaction time under each reaction condition was measured. Figure 1a compares particle formation and particle size distributions as a result of flowing a mixture of formic acid (140 ppm) and water vapor (630 ppm) in the absence and presence of 200 ppb TMA. The total number of particles and particle size distribution is plotted for both 8 and 48 s reaction times. Figure 1b offers a closer look at particle concentrations vs particle diameters ranging from 300–500 nm. Figure 1c shows particle concentrations vs diameter from 0–100 nm. No particles are observed while flowing only N2, 630 ppm H2O, 140 ppm formic acid or 200 ppb TMA. Particles did form as a result of flowing a mixture of formic acid (140 ppm) and water vapor (630 ppm); however, with the addition 200 ppb of TMA to the mixture, there is a dramatic increase in the formation of smaller particles at both reactions. The mixture of formic acid (140 ppm) and water vapor (630 ppm) without any TMA, generated a maximum of 3.3 ± 0.3 × 106 particles cm−3 at a diameter of 1.29 nm and a total number of particles of 7.2 ± 0.3 × 106 particles ranging in diameter from 1.09–498 nm. 99.4 % of the particles were <40 nm in diameter and 0.6% of particles were >40 nm. Inclusion of 200 ppb TMA to this formic acid/water vapor mixture resulted in an increase in the total number of particles to 7.1 ± 0.05 × 107 particles between 1.09 and 498 nm in diameter with a maximum of 3.4 ± 0.3 × 106 particles cm−3 generated at a diameter of 1.09 nm at 8 s of reaction time. The percentage of particles formed at 8 s with a diameter of <40 nm decreased to 98.7 % while the percentage of particles with a diameter of >40 nm rose to 1.3%. With an increase of reaction time to 48 s, the total number of particles rose to 4.1 ± 0.11 × 108 particles with a maximum concentration of 7.3 ± 0.9 × 106 particles cm−3 generated at a diameter of 1.54 nm. The percentage of particles formed <40 nm in diameter dropped to 98.6 % while the percentage of particles >40 nm rose to 1.4%. By comparison, the formic acid/water vapor mixture containing 200 ppb TMA shows an increased formation of both smaller diameter particles, <40 nm, and larger diameter particles (>300 nm), compared to the formic acid/water vapor mixture without TMA. The difference in the number of total particles and formation of particles (>300 nm) increases with longer reaction times for mixtures that include TMA.

Fig. 1
figure 1

Particle formation. Particle formation initiated by the addition of 140 ppm formic acid, 630 ppm water vapor, 200 ppb trimethylamine (TMA) at different reaction times (green-no TMA, red-8 s and blue-48 s). a Overall comparison of particle concentration vs particle diameter (0.8–500 nm). b Comparison of larger diameter particle concentration vs particle diameter (300–500 nm) at different times (blue- 48 s, red- 8 s, green-no TMA). c Comparison of smaller diameter particle concentration vs particle diameter (0.8–100 nm). Numerical values are given in Supplementary Table 3

Table 1 lists the conditions probed for each trial and the measured rate of particle formation. Figure 2 compares particle formation as a function of initial formic acid, water vapor, and TMA concentrations. Additional comparisons of Table 1 conditions are found in the supplemental data. See Supplementary Fig. 2 and Supplementary Note 1. Figure 2a shows that the particle formation occurs even at the lowest concentration of formic acid (140 ppm), water vapor (630 ppm) and TMA (200 ppb) in the reaction cell shown in red compared to increased TMA concentration (400 ppb) shown in black. The rate of particle formation under these conditions is measured to be 1.0 ± .02 × 106 particles cm−3 s−1. The rate of particle formation under higher concentrations could not be determined due to the aggregation of smaller sized particles leading to an increased concentration of larger particles, which skewed the calculations needed to determine the rate of particle formation. Figure 2a compares trials A and C (shown in black) with lower formic acid concentrations and an increase concentration of TMA from 200 ppb to 400 ppb, which shows a definite increase in the smaller particles at the shorter reaction times. Increasing the TMA concentration by a factor of 2 results in a 6.6% increase in the number of particles formed with diameters <20 nm at 8 s. A shift in the size distribution of smaller particles (<20 nm) and an increase in the number of particles >40 nm is also observed. In Fig. 2b, trials B (red) and D (black) are shown. The formic acid (140 ppm) and higher water vapor (1550 ppm) concentrations were kept constant but the concentration of TMA was double in trial D compared to B. The total number of particles grown in trial D at 8 s was 36% larger compared to that in B. In addition, evaluation of the particle size distribution shows 3.9% of the total particles in B are >20 nm in diameter but upon doubling of the concentration of TMA the total number of particles >20 nm in trial D increased to 4.2%. Figure 2c, which shows trials E (red) and G (black), the formic acid (540 ppm) and water concentrations (630 ppm) were kept constant but the TMA concentration was double in trial G compared to E. The total number of particles in trial G at 8 s was 17% smaller compared to that in E. In addition, evaluation of the particle size distribution shows that 3.1% of the total particles in E are >20 nm in diameter but upon increasing the concentration of TMA, the total number of particles >20 nm in trial G increased to 3.3%. In Fig. 2d, which shows trials F (red) and H (black), the formic acid (540 ppm) and water concentrations (630 ppm) were kept constant but the TMA concentration was double in trial H compared to F. The total number of particles grown in trial F at 8 s was 4.8% larger compared to that in H. In addition, evaluation of the particle size distribution shows that 4.1% of the total particles in F are >20 nm in diameter but upon increasing the concentration of TMA, the total number of particles >20 nm in trial H increased to 4.5%. Comparison of Fig. 2a–d (Supplementary Table 4) indicate an increase in the concentration of TMA increases the rate of particle formation and perturbs the particle size distribution a table of the values is found in supplemental materials. However, due to the high concentrations of formic acid and TMA in these trials, nucleation was not the only mechanism contributing to particle formation as is discussed later.

Table 1 Conditions probed for each trial and the measured rate of particle formation
Fig. 2
figure 2

Comparison of particle formation under different conditions. Inclusion of TMA into the formic acid water vapor system is shown to increase particle formation. Comparison of particle size distribution measured by an SMPS with different concentrations of formic acid, H2O and TMA. a Concentrations A (red) vs C(black); b concentrations B (red) vs D (black); c concentrations E (red) vs G (black); d concentrations F (red) vs H (black)

The effect of TMA on the kinetics of particle formation was next measured by changing the TMA concentrations and comparing the changed rates of particle formation. To overcome a secondary mechanism via aggregation and to determine the rate of particle formation via nucleation, the formic acid and TMA concentrations were decreased minimizing aggregation in the flow cell. A gas mixture of 60 ppm formic acid and 630 ppm water vapor with varying concentrations of TMA was introduced into the reaction cell. The concentration of TMA was varied in the reaction mixture between 12 and 20 ppb and the total number of particles >2.5 nm in diameter was measured at reaction times ranging between 8 and 48 s. These experiments allowed for measurement of the reaction rate coefficient, k > 2.5 nm, for new particle formation. The rate of particle formation was determined by finding the linear slope between the reaction times (8–48 s) and the total number of particles. Figure 3 shows a comparison of k > 2.5 nm for three different formic acid, water vapor, TMA mixtures. A linear relationship between total number of particles formed >2.5 nm in diameter and increasing concentration of TMA is observed. The rate of particle formation is observed to increase with the TMA concentration. The only anomaly observed in these series of experiments occurred when 200–400 ppb of TMA was introduced into the reaction cell. Under these conditions, the total number of particles formed, and shift in the particle size distribution did not follow the linear trend (Fig. 4), but instead, a decrease in rate of particle formation was measured. This rate decrease can be associated to the inability to measure particles with a diameter >500 nm using the SMPS’s used in this project. Figure 4a compares formation of particles >300 nm in diameter as a function of the reaction time under conditions of flowing 140 ppm formic acid, 630 ppm water vapor and 200 ppb TMA. At 8 s, the shortest reaction time probed, particle concentration peaks at 1.9 ± 0.4 × 104 particles cm−3 at 332 nm in diameter. As the reaction time increases to 24 s, the maximum concentration increases to 4.7 ± 0.1 × 104 particles cm−3 at a particle diameter of 332 nm. A long, shallow tail is observed in the observed particle diameters for both reaction times, 8 and 24 s. As the reaction time increases to 32 s the previously observed maximum in particle concentration observed in reaction times 8 and 24 s at 332 nm disappears, and a strong, broad distribution of particles >375 nm forms. At 48 s reaction time, a large peak in concentration is measured at 375 nm and a second broad tail is observed rising from 440 nm in diameter. Fine particles are classified into three distinct ranges: the ultrafine particle range (<10 nm), the transient nuclei range (10–100 nm) and the accumulation range (100 nm–1.2 μm) 50. Figure 4b separates particle diameters for each reaction time into these three modes versus particle concentration. At the earliest reaction time, there is a high particle concentration for the ultrafine particle range and a constant particle concentration for the accumulation range with an increase in reaction time. However, at the next shortest reaction time, there is an unexpected decrease in particle concentration in the ultrafine particle range followed by a slow increase in the ultrafine particle concentration with an increase in reaction time. In contrast, there is still an increase in the total number of particles in the accumulation range. The sudden decrease and subsequent increase in smaller diameter particles (<10 nm) along with the continuous increase in larger diameter particles indicates a two-fold nature of particle growth occurring in the reaction cell via nucleation and aggregation. To better comprehend the influence of TMA and H2O on the two-fold mechanism of particle growth, Fig. 5 compares the change in particle diameter distribution versus either a change in TMA or H2O concentration. Figure 5a analyzes particle concentration versus particle diameter distribution with 140 ppm formic acid, 630 ppm water vapor, 200 ppb TMA. Figure 5b shows a change in particle distribution with an increase in TMA. While Fig. 5c shows the change in particle distribution with an increase in H2O. The concentrations of H2O and TMA were respectively increased from 630 ppm to 1550 ppm, and 200 ppb to 400 ppb. Figure 5b indicates that with increased TMA concentrations particles with a diameter between 2.5 and 10 nm decreased and increased the number of particles with diameters between 100 and 500 nm six-fold at 8 s. Compared with Fig. 5c, there is a decrease in the number of particles with in diameter range diameters between 2.5 and 10 nm and a three-fold increase in particles with diameters between 100 and 500 nm at 8 s. Increased reaction times show the expected decrease in the overall particle concentration followed by an increase in the larger sized particles. The overall decrease is attributed to the inability to measure particles large than 500 nm in diameter. The aggregation of smaller sized particles leading to the formation particles with diameters 100–500 nm may be enhanced by the dipole moment caused by the TMA complexing with formic acid which in turn forms a hydrophobic/hydrophilic interaction as stated in the computational section. A similar enhancement may occur with increased H2O concentration, however to a lesser extent due to the smaller dipole moment. Ultrafine particles are formed as formic acid/water/TMA clusters grow with increasing reaction time in the cell. As the concentration of these particles increases the probability that these particles will collide with one another and aggregate to form a larger diameter particle increases. Aggregation of particles in the ultrafine and accumulation size range results in a decrease in the total number of these particles but an increase in the formation of particles >100 nm in diameter. We see both particle growth as water adheres to the nucleating site but also observe aggregation as small particles collide into one another to produce larger diameter particles.

Fig. 3
figure 3

Rate of particle formation. Addition of TMA to formic acid and water vapor system increased the rate of particle formation. Comparative rates of particle formation at 60 ppm formic acid, 630 ppm H2O and different TMA concentrations

Fig. 4
figure 4

Indepth comparison of particle formation. Inclusion of TMA into the formic acid and water vapor system increased the formation of larger diameter particles. Comparison of particle formation with 140 ppm formic acid, 630 ppm water vapor, 200 ppb TMA at different reaction times. a particle concentration vs particle diameter (300–493 nm) b bar graph comparing concentration to particle diameter sized (2.5–10 nm, 10–100 nm, and 100–493 nm) numerical values are given in Supplementary Tables 5, 6

Fig. 5
figure 5

Bar graph comparison of particle formation. Increasing water vapor or TMA concentrations increased the formation of larger diameter particles. Bar graph comparing concentration to particle diameter sized (2.5–10 nm, 10–100 nm, and 100–493 nm) a particle formation with 140 ppm formic acid, 630 ppm water vapor, 200 ppb TMA; b particle formation with 140 ppm formic acid, 630 ppm water vapor, 400 ppb TMA; c particle formation with 140 ppm formic acid, 1550 ppm water vapor, 200 ppb TMA). Numerical values are given in Supplementary Table 5

Simulation evidence of the formic acid–TMA interaction at air–water interface

The BOMD simulations provide useful information into the time scale, molecular mechanism of the ion-pair formation as well as the dynamic behavior of the ion-pair formed on the aqueous surface. The BOMD simulations were performed on a water droplet of 191 H2O molecules probing the nature of the interaction between HCOOH and (CH3)3N. The details of the simulation procedure are provided in the Methods section. The initial configuration of the HCOOH··(CH3)3N complex adsorbed on the water droplet surface are given in the Supplementary Fig. 3. We initiated BOMD simulations from hydrogen-bonded and non-hydrogen-bonded complexes of HCOOH with (CH3)3N. Surprisingly, only the hydrogen-bonded configurations resulted in the ion-pair formation on the water surface, implying that the formation of hydrogen-bonded HCOOH··(CH3)3N complex in the gas-phase constitutes a crucial step in the particle formation on the water surface. In a recent experimental study48, Eugene et al. also observed interfacial proton transfer from simpler carboxylic acids. Though the role of gas-to-particle partitioning in the particle formation has been speculated before8, our simulations provide a mechanistic rationale why the gas-to-particle conversion is actually required for the particle formation.

The HCOOH··(CH3)3N interaction follows a typical trajectory of acid-base chemistry and involves a proton transfer between HCOOH and (CH3)3N without the direct involvement of surface water molecules. This results in the formation of HCOO· ·(CH3)3NH + ion-pair on a picosecond (ps) time scale (Fig. 6). As just described, Eugene et al. have recently also observed interfacial proton transfer from simpler carboxylic acids that can play an important role in aerosol formation48. The role of water droplet here is to stabilize the ion-pair particle by forming a hydration shell around it. These findings are also consistent with field measurements predicting the presence of aminium salts in aged organic carbon particles in certain regions of California51,52. The reactive uptake of NH3 onto slightly soluble organic acid particle is also found to enhance the CCN activity and the formation rates of these particles53. For the reaction between HCOOH and (CH3)3N, the transition state like complex is formed at 5.18 ps (Fig. 6a and Supplementary Movie 1). In this complex, the hydroxyl proton of HCOOH is partially dissociated and transferred towards (CH3)3N, i.e., the O1–H1 bond is 1.33 Å long whereas the H1–N1 bond is 1.29 Å long. This complex is converted into the HCOO··(CH3)3NH+ ion pair at 5.23 ps. The O1–H1 bond is now lengthened to 1.70 Å, indicative of a hydrogen bonding interaction whereas the H1–N1 bond has become a true covalent bond (H1–N = 1.06 Å).

Fig. 6
figure 6

The simulation results detailing the interfacial proton transfer between formic acid and TMA ((CH3)3N). a Snapshot structures taken from the BOMD simulations of the reaction of formic acid with (CH3)3N, which illustrates the formation of HCOO··(CH3)3NH+ ion-pair on the water droplet. b Time evolution of key bond distances, O1–H1 and H1–N1 involved in the HCOO··(CH3)3NH+ ion-pair forming reaction. c Combined radial/radial distribution function involving H1–O1 and H1–N1 bond distances in the HCOO··(CH3)3NH+ ion-pair. d Combined distribution function involving an angular distribution function between H1–O1 and H1–N1 vectors and a radical distribution function between O1 and H1 bond distance in the HCOO··(CH3)3NH+ ion-pair

Dynamic behavior of HCOO··(CH3)3NH+ ion-pair at the air–water interface

To gain deeper insights into the dynamic behavior of the HCOO··(CH3)3NH+ ion-pair, we next analyzed its locus on the water droplet. Supplementary Fig. 4 shows the distance between the center of mass of the ion-pair and that of the water droplet as a function of the simulation time. The HCOO ··(CH3)3NH+ ion-pair is situated at 10–13 Å distance from the center of the water droplet, implying that it preferentially resides at the air-water interface. There are two main interactions responsible for the interfacial locus of the HCOO··(CH3)3NH+ ion-pair: (i) intraparticle hydrogen bonding (O1··H1–N1), and (ii) interfacial hydration shell. The intraparticle hydrogen bond in HCOO··(CH3)3NH+ ion-pair is quite strong as is evidenced from the combined distribution functions shown in Fig. 6c, d, respectively. The calculated average number of hydrogen bonds formed by the HCOO··(CH3)3NH+ ion-pair with interfacial water molecules provides useful insight into its interfacial hydration shell. Our model specified a hydrogen bond between a formate oxygen and H2O, if the O1/O2··H-O or O··H1–N1 distance was <2.5 Å and the O1/O2··H-O or O··H1–N1 hydrogen bond angle was >150°. The COO group forms 1.9 hydrogen bonds whereas the aminium proton does not form any hydrogen bond with H2O molecules. The hydrophobicity of trimethyl group in the HCOO··(CH3)3NH+ ion-pair prevents any interaction between the aminium proton and interfacial water molecules. To deeply understand the solvation structure of HCOO··(CH3)3NH+ ion-pair at the air-water interface, we next identified key [m,n] configurations and calculated their probabilities (Fig. 7). Here m and n are the number of interfacial H2O molecules bound to the HCOO and (CH3)3NH+, respectively. The configurations, in which only the HCOO group of HCOO··(CH3)3NH+ ion-pair binds to one, two and three interfacial H2O molecules, are the most probable ones and account for 31%, 38, and 26% of the total configurations, respectively (Fig. 7). Since there remains significant uncertainty about the exact composition of the particle in air8,37,38,39,40,41,42, our results may play a crucial role in revealing a fundamental yet vital piece of information on the particle formation on the water surfaces.

Fig. 7
figure 7

Simulation results on the hydration behavior of the HCOO··(CH3)3NH+ ion-pair on the water surface. a Schematic showing the [m,n] interfacial waters forming hydrogen bonds with oxygens and aminium protons and of HCOO··(CH3)3NH+ where m and n are the number of interfacial water molecules bound to HCOO and (CH3)3NH+, respectively. b Histograms of probabilities of different [m,n] configurations for HCOO··(CH3)3NH+

Atmospheric implications

The growth of atmospheric nanoparticles was initially proposed to occur via direct condensation of organic acids54. However, the high saturation pressure of organic acids over nanoparticles55 negates the possibility of such a mechanism. On contrary, thermal desorption chemical ionization mass spectrometry (TDCIMS) analysis of the particles formed in Tecamac, Mexico shows the presence of carboxylic and hydroxycarboxylic organic acids13, suggesting that organics play a key role in the nanoparticle growth.

Our experimental and theoretical results may help in resolving the discrepancy in the role of organic acids in nanoparticle growth. The current BOMD simulations suggest a new gas-to-particle partitioning mechanism for the formation of organic particles that involves two steps: (i) formation of hydrogen-bonded complexes of organic acids with amines in the gas-phase, and their adsorption on the aqueous surface, and (ii) subsequent proton transfer between hydrogen-bonded entities. The mechanistic beauty of this new mechanism lies in the fact that it does not require direct condensation of organic acids, but rather involves hydrogen bonding between organic acids and amines. The organic-amine salts formed in this mechanism will have reduced volatility, as has been pointed out in a recent study56. This new mechanism is supported by combined laboratory investigation and field measurements using TDCIMS and ultrafine hygroscopicity tandem differential mobility analyzer (TDMA) confirming that carboxylate-alkylaminium salts contribute to the particle growth in atmosphere17. Additional support for this mechanism comes from the study of Dinar et al. showing that the reactive uptake of NH3 can enhance the CCN activity and hygroscopic growth of organic acid particles53. Our BOMD simulations suggest that the carboxylate-aminium ion-pair remained stable on the water surface over the simulated time scale of 20 ps. This is again consistent with high thermal stability of alkylaminium carboxylate salts from the organic acid-amine reactions57.

Discussion

Our experiments show the nanoparticle formation under conditions that are typically outside what is expected if you use classical nucleation theory to describe particle formation. That is to say, we see nanoparticle growth under conditions of <100% RH. Under the classical nucleation model, it is suggested that 100% RH is necessary for particle formation and growth to occur. We speculate that by adding trace amounts of amines and/or formic acid, we may “seed” a cloud to increase the rate of particle formation. Our experimental results are consistent with what is observed around the world in terms of measuring carboxylic acids in particles as shown in the Sao Paulo, Brazil and Tecama, Mexico studies. The conclusions from our work could be used to better understand the role of other monocarboxylic acids such as CH3COOH, and dicarboxylic acids such as oxalic acid and malonic acid in the particle formation under diverse environmental conditions, including semiurban sites in the northeastern United States58, urban environments59,60,61, and remote locations62.

Though considered unlikely due to their high vapor pressures atmospheric gaseous- and particle-phase, carboxylic acids have been measured in the urban environment of Sao Paulo, Brazil63. Carboxylic acids accounted for a fraction of 7% of the total organic carbon in the particle phase, with formate and oxalate being two of the most abundant carboxylates. Though the major source of these carboxylates is suggested to be traffic emissions, the comprehensive insights into their emission source are still lacking. Our results suggest that the acid-base chemistries between organic acids and organic nitrogens could be an overlooked source of organic particulate in urban air. The organic nitrogens could be transported from the rural agricultural sites to the urban regions where they react with organic acids and result in the particulate formation. These new findings may help in improving the accuracy of existing atmospheric models.

Methods

Experimental details

A flow cell set-up was used to experimentally verify the formation of particles initiated by formic acid, water vapor, and TMA. Supplementary Fig. 1 shows the experimental setup, a 180 cm long Pyrex flow cell (5.1 cm i.d.) was connected to two aluminum boxes (26.5 cm × 26.5 cm × 26.5 cm) that allowed for UV light to pass through the flow cell for detection of formic acid and TMA. Formic acid, water vapor, and N2 were separately introduced at the top of the cell. TMA was introduced at varying points in the flow cell by a Teflon coated shower ring. Particle size distribution with particle diameter sizes ranging from 1.09 nm to 493.95 nm were analyzed using two Grimm Aerosol Technik scanning mobility particle sizer (SMPS, model: 5.710) each consisting of a differential mobility analyzer (DMA) Electrostatic Classifier (Vienna/Reishl type, model 55-u; 55–100) and a faraday cup electrometer (FCE, model 5.705).

Prior to each experiment, the flow cell was cleaned out with ozone, N2 (2 SLM) and O2 (100 sccm) for 1 h between experiments and overnight. This allowed for the complete removal of any formic acid or TMA that had been introduced in previous experiments. The flow cell was kept at a constant temperature (22 ± 1 °C) and pressure (650 ± 5 Torr). In all experiments, formic acid and water vapor were introduced separately into the flow cell by bubbling N2 (100–300 sccm) through 88% proof formic acid and purified water. Additional N2 was introduced into the flow cell for a total flow rate of 3.125 SLM. As previously mentioned TMA (25–125 sccm) was introduced at varying points in the flow cell by a Teflon coated shower ring. The shower ring (i.d. 4.5 cm) with pin holes was attached to a stainless-steel rod which allowed for movement up and down the flow cell. This method of introducing TMA allows for reaction times varying from 8 to 48 s. A range of concentrations for formic acid and TMA were varied by varying the flow of gases and subsequently using UV absorption spectroscopy to determine formic acid and TMA concentrations. Using Beer’s Law and previously published cross sections from 220.16 to 250.9 nm the concentrations of TMA were determined to be 200 and 400 ppb. Using previously published formic acid cross sections from 210.88 to 230.86 nm the formic acid concentrations were determined to be 140 and 540 ppm. A water vapor calibration curve was determined by introducing known concentrations of water vapor via a syringe pump (KD Scientific KDS-100 Syringe Pump) at known flow rates. Water was introduced into a stainless-steel tee packed with glass wool heated to 150 °C. N2 (3.1 SLM) passed through the glass wool into the flow cell and the absorption from the water concentration was determined by dithering an IR diode over 1380.47 nm to 1384.6 nm and the area under the desired peak was integrated for the known concentration of water vapor. Water vapor concentrations used were 630 and 1550 ppm.

Computational details

Born–Oppenheimer molecular dynamics (BOMD) simulations were performed based on a density functional theory (DFT) method implemented in the CP2K64 code. In the BOMD simulation, the droplet system contained 191 water molecules, one HCOOH molecule and one N(CH3)3 molecule. We have explored a total number of 10 different trajectories, in which five of them started from non-hydrogen-bonded configurations of HCOOH and N(CH3)3 whereas the other five started from hydrogen-bonded configurations. The dimension of the simulation box is x = 35 Å, y = 35 Å, z = 35 Å, which is large enough to neglect interactions between adjacent periodic images of the water droplet. Prior to the BOMD simulation, the system was fully relaxed using a DFT method, in which the exchange and correlation interaction is treated with the Becke-Lee-Yang-Parr (BLYP) functional65,66. The Grimme’s dispersion correction method is applied to account for the weak dispersion interaction67,68. A double-ζ Gaussian basis set combined with an auxiliary basis set and the Goedecker-Teter-Hutter (GTH) norm-conserved pseudopotentials were adopted to treat the valence electrons and the core electrons, respectively69,70. An energy cutoff of 280 Rydberg was set for the plane-wave basis set and 40 Rydberg for the Gaussian basis set. The BOMD simulations were carried out in the constant volume and temperature (NVT) ensemble, with the Nose-Hoover chain method for controlling the temperature (300 K) of the system. The integration step was set as 1fs, which had been proven to achieve sufficient energy conservation for the water system.