Next Article in Journal
The Effect of Long-Term Administration of Fatty Acid Amide Hydrolase Inhibitor URB597 on Oxidative Metabolism in the Heart of Rats with Primary and Secondary Hypertension
Next Article in Special Issue
Ab Initio Molecular Dynamics Simulation Study on the Stereo Reactions between Atomic Oxygen Anion and Methane
Previous Article in Journal
QSAR Study of N-Myristoyltransferase Inhibitors of Antimalarial Agents
Previous Article in Special Issue
One-Electron Reduction Potentials: Calibration of Theoretical Protocols for Morita–Baylis–Hillman Nitroaromatic Compounds in Aprotic Media
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 23
Issue 9
10.3390/molecules23092349
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Theoretical Insights into the Electron Capture Behavior of H2SO4···N2O Complex: A DFT and Molecular Dynamics Study

by
Wei-Hua Wang
*,
Wen-Ling Feng
,
Wen-Liang Wang
and
Ping Li
*
Key Laboratory of Life-Organic Analysis, School of Chemistry and Chemical Engineering, Qufu Normal University, Qufu 273165, China
*
Authors to whom correspondence should be addressed.
Molecules 2018, 23(9), 2349; https://doi.org/10.3390/molecules23092349
Submission received: 18 August 2018 / Revised: 4 September 2018 / Accepted: 5 September 2018 / Published: 13 September 2018
(This article belongs to the Special Issue Theoretical Investigations of Reaction Mechanisms)
Browse Figures
Versions Notes

Abstract

:
Both sulfuric acid (H2SO4) and nitrous oxide (N2O) play a central role in the atmospheric chemistry in regulating the global environment and climate changes. In this study, the interaction behavior between H2SO4 and N2O before and after electron capture has been explored using the density functional theory (DFT) method as well as molecular dynamics simulation. The intermolecular interactions have been characterized by atoms in molecules (AIM), natural bond orbital (NBO), and reduced density gradient (RDG) analyses, respectively. It was found that H2SO4 and N2O can form two transient molecular complexes via intermolecular H-bonds within a certain timescale. However, two molecular complexes can be transformed into OH radical, N2, and HSO4− species upon electron capture, providing an alternative formation source of OH radical in the atmosphere. Expectedly, the present findings not only can provide new insights into the transformation behavior of H2SO4 and N2O, but also can enable us to better understand the potential role of the free electron in driving the proceeding of the relevant reactions in the atmosphere.

1. Introduction

As the primary component in the production of secondary aerosols, sulfuric acid (H2SO4) plays a central role in atmospheric aerosol nucleation [1,2,3,4]. It is mainly produced from sulfur dioxide via oxidation by hydroxyl radical. Generally, when H2SO4 collides with other molecules, it can form small clusters of molecules that can grow into new stable aerosols [1]. As a result, many interaction models of H2SO4 with water, ammonia, and organic molecules have been constructed to elucidate new atmospheric particle formation and nucleation mechanisms [5,6,7,8,9,10]. Meanwhile, H2SO4 also plays an important role in perturbing atmospheric hydroperoxyl radical levels via the formation of strong H-bonding complexes [11,12]. In addition, H2SO4 has been involved in the formation of the polar stratospheric clouds (PSCs) associated with the stratospheric ozone depletion [13]. Therefore, it is necessary to acquire the detailed interactions of H2SO4 with other atmospheric species so as to better understand the potential role of H2SO4 in the atmosphere.
As one of the important nitrogen-containing species in the atmosphere and a main greenhouse gas, nitrous oxide (N2O) has a long atmospheric lifetime and radiative forcing of about 200 times higher than that of CO2 [14,15]. It has been demonstrated that the emission of N2O is mainly from both natural (about 60%) and anthropogenic sources (about 40%) [16]. In 2016, the globally averaged N2O concentration reached 328.9 ± 0.1 ppb, which is 122% of the pre-industrial level (270 ppb) [16]. In addition, N2O has also been implicated in the ozone depletion. Ravishankara’s [17] studies suggested that N2O emission is the single most important ozone-depleting emission and is expected to remain the largest throughout the 21st century. Therefore, it is also necessary to investigate the transformation of N2O to reduce its negative effects on environments.
Given the fact that H2SO4 and N2O coexist in the atmosphere, it is expected that both of them could interact with each other in some way. Can they form stable complexes? If so, what are the interaction modes adopted between them? Especially, what happens to these complexes if an additional electron is introduced since a large number of electrons can be generated by cosmic ray radiation in the atmosphere? Here, it is well known that the excess electrons can be captured by the stratospheric molecules [18,19,20] and the cosmic-ray-driven electron-induced reactions have been involved in the ozone depletion and global climate change [21,22,23]. Unfortunately, despite the atmospheric importance of H2SO4 and N2O, the interaction behavior between them in the absence and presence of an additional electron remains unclear. No relevant studies have been reported so far to the best of our knowledge. In this case, theoretical studies are highly desirable to carry out to provide useful information for the relevant experiments.
Therefore, to address those questions mentioned above, the interaction behavior between H2SO4 and N2O before and after the electron capture has been systematically explored by using the density functional theory (DFT) method and ab initio molecular dynamics simulation. Hopefully, the present study can provide new insights into the potential role of H2SO4 and the excess electrons in the atmosphere as well as the transformation of N2O.

2. Computational Details

In this study, full geometric optimizations were performed employing the B3LYP and M06-2X methods as well as the 6-311++G(3df,3pd) basis set, where the reliability of the DFT method has been confirmed by many systems [24,25,26,27,28,29]. As a result, it was found that both methods can give consistent results for the geometry and energy parameters overall. Given the compromise between computational accuracy and cost, and to keep the consistency of the results for the optimization and molecular dynamics at the same level of theory, the results at the B3LYP/6-311++G(3df,3pd) level of theory have been mainly discussed below unless otherwise noted. To explore the nature of the optimized structure, we performed vibration frequency analysis at the same level of theory. Moreover, intrinsic reaction coordinate (IRC) calculations have been performed to further confirm the connectivity between the given transition state and the corresponding reactant and product. In addition, to further improve the energy parameters in the formation process of the neutral complexes, single-point energy calculations have also been carried out at the CCSD(T)/AUG-cc-pVTZ level of theory based on the optimized geometries mentioned above.
To explore the implicit solvent effects on the title reaction, the solvent model density (SMD) solvation model has been employed in aqueous solution at the B3LYP/6-311++G(3df,3pd) level of theory. As a result, similar results have been obtained as those of the results in the gas phase.
To characterize the intermolecular H-bonds between H2SO4 and N2O before and after electron capture, atoms in molecules (AIM) [30], natural bond orbital (NBO) [31], and the reduced density gradient (RDG) analyses have been performed based on the optimized geometries [32,33]. In the AIM theory, the presence and nature of the intermolecular H-bonds can be characterized by the bond critical point (BCP). As for the formation mechanism of the intermolecular H-bonds, it can be elucidated from the orbital interactions between the lone pairs of the proton acceptor and the anti-bonding σ orbital (σ*) of the proton donor within the NBO method, accompanying with the electron transfer from the former to the latter. Moreover, the orbital interaction strength can be assessed from the second-order stabilization energy E(2). In addition, the intermolecular H-bonds can also be indicated by the presence of the spikes within the RDG analyses.
To obtain the interaction strength between H2SO4 and N2O, interaction energy between them has been calculated as the energy differences between the formed complexes and the reactants. Meanwhile, zero-point vibrational energy (ZPVE) corrections and basis set superposition errors (BSSEs) have also been considered [34].
To investigate the dynamic stability of the formed complexes and the microscopic details of the electron capture process for those complexes, we performed ab initio molecular dynamics simulation employing the Atom Centered Density Matrix Propagation (ADMP) method at the B3LYP/6-311++G(3df,3pd) level of theory on the basis of their optimized geometries [35,36,37]. The whole simulation was performed at 298 K for the neutral complexes and the electron captured products with a canonical ensemble. For the electron capture processes of the neutral complexes, a microcanonical ensemble has been employed. The total simulation time is 1 ps with a time step of 0.5 fs.
To evaluate the ability to accept an electron for those neutral complexes, the adiabatic electron affinity (AEA), vertical electron affinity (VEA), and vertical electron detachment energy (VEDE) have been calculated. Here, AEA is calculated as the enthalpy difference between the optimized neutral and anionic complexes and VEA is calculated as the electronic energy difference between the optimized neutral and anionic states on the basis of the optimized neutral complexes. As for the VEDE, it is calculated as the electronic energy difference between the neutral and anionic states at the same geometry of the anionic state.
All the calculations were completed employing Gaussian 09 program (Gaussian, Inc., Wallingford, CT, USA) [38].

3. Results and Discussion

3.1. Formation of Neutral Molecular Complexes

3.1.1. Structural Features

As displayed in Figure 1, two complexes named as A and B have been formed between H2SO4 and N2O through intermolecular H-bonds. Namely, one of the protons in H2SO4 interacts with the terminal N and O atoms of N2O in complexes A and B, respectively. This point can be confirmed by the location of the corresponding BCPs between H2SO4 and N2O as shown in Figure 1. The intermolecular H-bond distances of 2.030 (2.171) and 1.935 (1.895) A suggest that the H-bonding interactions should be weak, where the data in parentheses refer to the results at the M06-2X/6-311++G(3df,3pd) level of theory. Actually, as shown in Table 1, the positive Laplacians and energy densities at the BCPs of the H-bonds indicate that the above H-bonds are predominated by electrostatic interaction.
To further verify the formation of intermolecular H-bonds mentioned above, we have performed RDG analyses for the complexes A and B. As shown in Figure 2, the intermolecular H-bonds can be verified by the presence of the corresponding spikes between H2SO4 and N2O.
To explore the formation mechanisms for the above intermolecular H-bonds, NBO analyses have been carried out based on the optimized complexes. As shown in Figure 3, the formation of the intermolecular H-bonds in complexes A and B is mainly attributed to the electron transfer from the lone pairs of the terminal N and O atoms of N2O to the anti-bonding σ orbital (σ*) of the O4-H5 bond of H2SO4, respectively. Moreover, as given in Table 2, the larger total second-order stabilization energies of the orbital interactions in complex B imply that the H-bond in complex B is stronger than that of complex A, which can be further confirmed by the interaction energy mentioned below.
According to the NBO analyses, the O4-H5 bond of H2SO4 participating in the H-bonds should be elongated due to the electron transfer from the lone pairs of the proton acceptor to the σ* orbital of the proton donor. As a result, compared with the isolated state, the O4-H5 bond of H2SO4 has been elongated by about 0.006 and 0.008 A, respectively. Correspondingly, as shown in Figure 4, red-shifts of 115 and 152 cm−1 occur for the vibrational frequencies of the O4-H5 bond in complexes A and B, respectively. Moreover, the characteristic absorption peaks at 3650.0 and 3612.8 cm−1 are originated from the stretching vibration of the O4-H5 bond of H2SO4 fragment in complexes A and B, providing the direct spectra proof for the existence of intermolecular H-bonds. For the sake of reference, the detailed IR results including peak positions and intensities have been given in Table S1 of the Supplementary Material (SM).
To investigate the electron redistribution details for H2SO4 and N2O upon complexation, the electron density difference maps of them have been constructed based on the optimized complexes A and B. As shown in Figure 5, electron redistribution mainly occurs along the intermolecular H-bond axis in the complexes, which is consistent with the features of the H-bonding complexes [39,40,41]. In more details, in complex A, the loss of the electron is mainly concentrated on the positions of the proton involved in the intermolecular H-bond in H2SO4 and the other two atoms of N2O not participating in the H-bonds. As for the gain of the electron, it is mainly concentrated on the terminal N atom of N2O and the O atom of H2SO4 involving in the H-bonds. Similarly, the same is also true for complex B. As a result, it was found that N2O and H2SO4 act as an electron donor and electron acceptor in both complexes, where the net electron transfer is 0.014 and 0.017 in complexes A and B based on the NBO analyses, respectively.

3.1.2. Energy Analyses

As displayed in Table 3, the calculated interaction energies between H2SO4 and N2O are −2.14(−2.72) and −2.77(−4.22) kcal/mol in complexes A and B, respectively, where the data in parentheses refer to the results at the CCSD(T)/AUG-cc-pVTZ level of theory. As a result, complex B is more stable than complex A by 0.76 (1.72) kcal/mol, which is in agreement with the larger H-bond in complex B.
Thermodynamically, the calculated enthalpy changes (ΔH) for complexes A and B are all negative, indicating that the formation processes of them are exothermic reactions. On the other hand, the small positive values for the Gibbs free energy changes (ΔG) suggest that the formation processes of them are non-spontaneous at 298 K. As shown in Figure 6, the ΔG values become more negative with the decreasing of the temperature. Therefore, low temperature should be favorable for the formation of complexes thermodynamically.

3.1.3. Molecular Dynamics Analyses

To further explore the dynamic stability of the complexes A and B, molecular dynamics simulations have been performed on the basis of their optimized geometries. The time evolution of the potential energy and the H-bond in both complexes have been given in Figure 7. On the one hand, both intermolecular H-bonds in complexes can still be observed during the initial stage of the dynamics simulation, indicating that the complexes can exist in a short period of time. On the other hand, the intermolecular H-bonds disappear with the evolution of the simulation time, implying the instability of the complexes. Therefore, complexes A and B can exist transiently within a certain period of time.

3.2. Electron Capture Process

On the basis of the neutral complexes A and B mentioned above, the electron capture properties of them have been investigated.

3.2.1. Structural Features for the Electron Capture Products

As shown in Figure 8, complex A′ has been located for the complex A upon electron capture. Obviously, the original H5 atom of H2SO4 fragment has been transferred to the terminal N10 atom of N2O, resulting in the formation of two fragments HSO4− and HNNO as well as a new intermolecular H-bond N10-H5···O4. At the same time, an additional intermolecular H-bond O6-H7···O8 has also formed between the other proton of H2SO4 and the terminal O8 atom of N2O. As shown in Figure 9, complex A′ is characterized by the double intermolecular H-bonds from the presence of the BCPs and the spikes between two fragments in RDG maps as displayed in Figure S1 of Supplementary Materials. Especially, for the intermolecular H-bond N10-H5···O4, its larger strength should be highlighted from its short H-bonding distance of 1.491 A. Actually, as presented in Table 1, the positive and negative values of the Laplacian and energy density at the BCP of H-bond suggest that it has partial covalent properties. In this case, can H atom transfer occur for the H-bond N10-H5···O4? To clarify this point, the H atom transfer behavior has been studied for complex A′. As a result, it was found that another complex A′′ can be obtained via the double H atom transfer, i.e., the H5 atom has been transferred from N10 of N2O to O4 of H2SO4 accompanying with the transfer of H7 from O6 of H2SO4 to O8 of N2O, where the corresponding transition state A′-TS has been confirmed by the IRC calculations. Similar to complex A′, as shown in Figure 9 and Figure S1 of Supplementary Materials, complex A′′ is also characterized by the double intermolecular H-bonds. Moreover, analyses of the energy parameters suggest that the forward and reverse barrier heights are 0.31(−1.44) and 9.05(7.62) kcal/mol, where the data in parentheses are the results including the ZPVE corrections. Clearly, the disappearance of the forward barrier height suggests that the reaction can take place readily. Meanwhile, the calculated ∆H and ∆G for the conversion process from complex A′ to A′′ are −8.74 and −9.39 kcal/mol, implying that the conversion process is also favorable thermodynamically. Therefore, complex A should be converted to complex A′′ upon electron capture, where the dynamic stability of the complex A′′ will be explored below. In addition, as displayed in Figure 10, the excess electron has been mainly distributed on the HNNO and NNOH fragments, exhibiting the nature of them as a radical.
For complex B, as displayed in Figure 8, its original intermolecular H-bond has disappeared upon electron capture and the whole complex has been collapsed into three fragments including OH radical, N2, and HSO4−. As displayed in Figure 10, the excess electron is mainly localized on the OH fragment. As a result, OH radical can be produced in the electron capture process of complex B. Therefore, besides the main source of the OH production (e.g., solar radiation), the present study can also provide an alternative source of the OH radical in the atmosphere although its contribution to the total OH concentration is unknown at present.
In addition, to verify the fact that the above electron capture products are free from the overestimated stabilization of the delocalized state for B3LYP method, relevant calculations have also been carried out employing the long range corrected functional CAM-B3LYP method and the 6-311++G(3df,3pd) basis set. As a result, the calculated results are well consistent with the B3LYP results.

3.2.2. Electron Affinity

To explore the ability of the neutral complexes A and B to accept an electron thermodynamically, the electron affinities of them have been calculated. As shown in Table 4, the AEAs are 1.80 and 3.18 eV for complexes A and B, respectively. Here, the small AEA of complex A relative to complex B should be attributed to the instability of the electron capture product A′ mentioned above. The positive values of AEA suggest that the electron capture products are more stable than those of the neutral complexes. At the same time, the large difference between AEA and VEA suggests the significant geometrical changes for complexes A and B before and after electron capture. Obviously, this point is consistent with the optimized geometries for complexes A′ and B′ mentioned above. In addition, the calculated VEDEs of the complexes A′ and B′ are 5.22 and 7.44 eV, respectively. Therefore, the electron capture products are stable relative to the detachment of an electron.
In addition, the ∆G for the electron capture process of complexes A and B are −38.9 and −80.2 kcal/mol, respectively. Therefore, the electron capture processes for the neutral complexes are exothermic and spontaneous processes and can occur easily thermodynamically.

3.2.3. Electron Capture Dynamics

To further explore the microscopic details in the electron capture process, the molecular dynamics calculations have been carried out based on the optimized complexes A and B.
As shown in Figure 11 and Figure S2 of Supplementary Materials, for complex A upon electron capture, H2SO4 and N2O begin to approach to each other, accompanying with the strengthening of the intermolecular H-bond. Meanwhile, the O4-H5 bond participating in the H-bond in H2SO4 fragment increases gradually. On the contrary, the distance between the terminal N10 atom of N2O and the H5 atom of H2SO4 fragment decreases simultaneously. As a result, the H5 atom of H2SO4 fragment has been transferred to the terminal N10 atom of N2O fragment. For example, at 88.5 fs, the distances N10···H5 and O4···H5 have been changed from the initial 2.030 and 0.974 A to 0.981 and 1.699 A, respectively. At 134.5 fs, the H5 atom of H2SO4 has been transferred to the N10 atom of N2O, producing the HNNO and HSO4− fragments, which is consistent with the optimized complex A′.
Moreover, as mentioned above, complex A′ can be transformed into complex A′′ easily. To further explore the dynamic stability of complex A′′, molecular dynamics calculations have been carried out based on its optimized geometry. As displayed in Figure 12 and Figure S3 of Supplementary Materials, the N-O bond of HONN fragment has been elongated gradually with the evolution of simulation time. At 60 and 90 fs, the N-O bond has been elongated to 2.128 and 2.950 A, respectively. Therefore, complex A′′ is unstable and can be further transformed into OH radical, N2, and HSO4− fragments.
As for the complex B upon electron capture, as shown in Figure 13 and Figure S4 of Supplementary Materials, the H5 atom involved in the intermolecular H-bond in H2SO4 fragment approaches the terminal O8 atom of N2O fragment gradually, accompanying with the elongation of the N9-O8 bond of N2O fragment simultaneously. At 50 fs, the O4-H5 bond of H2SO4 and N9-O8 bond of N2O have been changed from 0.976 and 1.191 to 1.473 and 2.047 A, respectively. With the evolution of the simulation time, complex B has been transformed into OH radical, N2, and HSO4− fragments, which is consistent with the optimized geometry of complex B′.
In addition, as shown in Figure 11 and Figure 13, the energies of the neutral complexes have been decreased significantly upon electron capture, reflecting the easy proceeding of the electron capture process once again.
As a result, the neutral complexes A and B can be eventually transformed into OH radical, N2, and HSO4− fragments upon electron capture. Certainly, further experiments are needed to further confirm the present findings.

4. Conclusions

In the present study, the electron capture properties of the formed complex of H2SO4 with N2O in the gas phase have been systematically explored theoretically. It was found that H2SO4 and N2O can form two transient complexes via intermolecular H-bonds in the atmosphere, which has been characterized by using AIM, NBO, and RDG analyses. Upon electron capture, two neutral complexes have been eventually transformed into OH radical, N2, and HSO4− species. Therefore, the present study not only provides an alternative formation source of OH radical in the atmosphere as well as new insights into the degradation and transformation behavior of N2O, but also enables us to better understand the potential role of H2SO4 and the excess electron in the atmosphere. Certainly, more extensive experiments are highly desirable to verify the present findings.

Supplementary Materials

RDG maps of the electron capture products, the time evolution of the selected bond distances for complexes A, A′′, and B in the molecular dynamics calculations, and the IR data for complexes A and B. The following are available online, Figure S1: RDG maps of the electron capture products, Figure S2: The evolution of the selected bond distances for complex A in the electron capture process as a function of time, Figure S3: The evolution of the selected bond distances for complex A′′ as a function of time, Figure S4: The evolution of the selected bond distances for complex B in the electron capture process as a function of time, Table S1: The calculated IR results for complexes A and B.

Author Contributions

Data curation, W.-H.W. and W.-L.W.; Investigation, W.-H.W.; Methodology, W.-L.F.; Software, P.L.; Supervision, P.L.; Writing—Original draft, P.L.

Funding

This work is supported by National Natural Science Foundation of China (Nos. 21577076, 21303093, 21805163), Natural Science Foundation of Shandong Province (No. ZR2018MB020), and the Doctoral Foundation of Shandong Province (Nos. ZR2016BB20, ZR2017BB055).

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Tomicic, M.; Enghoff, M.B.; Svensmark, H. Experimental study of H2SO4 aerosol nucleation at high ionization levels. Atmos. Chem. Phys. 2018, 18, 5921–5930. [Google Scholar] [CrossRef]
  2. Brus, D.; Škrabalová, L.; Herrmann, E.; Olenius, T.; Trávničková, T.; Makkonen, U.; Merikanto, J. Temperature-dependent diffusion of H2SO4 in air at atmospherically relevant conditions: Laboratory measurements using laminar flow technique. Atmosphere 2017, 8, 132. [Google Scholar] [CrossRef]
  3. Sipilä, M.; Berndt, T.; Petäjä, T.; Brus, D.; Vanhanen, J.; Stratmann, F.; Patokoski, J.; Mauldin, R.L., III; Hyvärinen, A.; Lihavainen, H.; et al. The role of sulfuric acid in atmospheric nucleation. Science 2010, 327, 1243–1246. [Google Scholar]
  4. Anglada, J.M.; Olivella, S.; Sole, A. Hydrogen transfer between sulfuric acid and hydroxyl radical in the gas phase: Competition among hydrogen atom transfer, proton-coupled electron-transfer, and double proton transfer. J. Phys. Chem. A 2006, 110, 1982–1990. [Google Scholar] [CrossRef] [PubMed]
  5. Almeida, J.; Schobesberger, S.; Kürten, A.; Ortega, I.K.; Kupiainen-Määttä, O.; Praplan, A.P.; Adamov, A.; Amorim, A.; Bianchi, F.; Breitenlechner, M.; et al. Molecular understanding of sulphuric acid–amine particle nucleation in the atmosphere. Nature 2013, 502, 359–363. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Zhang, R.; Suh, I.; Zhao, J.; Zhang, D.; Fortner, E.C.; Tie, X.; Molina, L.T.; Molina, M.J. Atmospheric new particle formation enhanced by organic acids. Science 2004, 304, 1487–1490. [Google Scholar] [CrossRef] [PubMed]
  7. Kirkby, J.; Curtius, J.; Almeida, J.; Dunne, E.; Duplissy, J.; Ehrhart, S.; Franchin, A.; Gagné, S.; Ickes, L.; Kürten, A.; et al. Role of sulphuric acid, ammonia and galactic cosmic rays in atmospheric aerosol nucleation. Nature 2011, 476, 429–433. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Ehn, M.; Thornton, J.A.; Kleist, E.; Sipilä, M.; Junninen, H.; Pullinen, I.; Springer, M.; Rubach, F.; Tillmann, R.; Lee, B.; et al. A large source of low-volatility secondary organic aerosol. Nature 2014, 506, 476–479. [Google Scholar] [CrossRef] [PubMed]
  9. Beichert, P.; Schrems, O. Complexes of sulfuric acid with hydrogen chloride, water, nitric acid, chlorine nitrate, and hydrogen peroxide: An ab initio investigation. J. Phys. Chem. A 1998, 102, 10540–10544. [Google Scholar] [CrossRef]
  10. Metzger, A.; Verheggen, B.; Dommen, J.; Duplissy, J.; Prevot, A.S.H.; Weingartner, E.; Riipinen, I.; Kulmala, M.; Spracklen, D.V.; Carslaw, K.S.; et al. Evidence for the role of organics in aerosol particle formation under atmospheric conditions. Proc. Natl. Acad. Sci. USA 2010, 107, 6646–6651. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Miller, C.E.; Francisco, J.S. The formation of a surprisingly stable HO2-H2SO4 complex. J. Am. Chem. Soc. 2001, 123, 10387–10388. [Google Scholar] [CrossRef] [PubMed]
  12. Li, P.; Ma, Z.; Wang, W.; Zhai, Y.; Sun, H.; Bi, S.; Bu, Y. Theoretical studies on the coupling interactions in H2SO4···HOO●···(H2O)n(n = 0–2) clusters: Toward understanding the role of water molecules in the uptake of HOO● radical by sulfuric acid aerosols. Phys. Chem. Chem. Phys. 2011, 13, 941–953. [Google Scholar] [CrossRef] [PubMed]
  13. Molina, M.J.; Zhang, R.; Wooldridge, P.J.; McMahon, J.R.; Kim, J.E.; Chang, H.Y.; Beyer, K.D. Physical chemistry of the H2SO4/HNO3/H2O system: Implications for polar stratospheric clouds. Science 1993, 261, 1418–1423. [Google Scholar] [CrossRef] [PubMed]
  14. Bouwman, A.F.; Van der Hoek, K.W.; Olivier, J.G.J. Uncertainties in the global source distribution of nitrous oxide. J. Geophys. Res. 1995, 100, 2785–2800. [Google Scholar] [CrossRef]
  15. Crutzen, P.J.; Mosier, A.R.; Smith, K.A.; Winiwarter, W. N2O release from agro-biofuel production negates global warming reduction by replacing fossil fuels. Atmos. Chem. Phys. 2008, 8, 389–395. [Google Scholar] [CrossRef] [Green Version]
  16. World Meteorological Organization. The state of greenhouse gases in the atmosphere based on global observations through 2016. WMO Greenh. Gas Bull. 2017, 13, 1–8. [Google Scholar]
  17. Ravishankara, A.R.; Daniel, J.S.; Portmann, R.W. Nitrous oxide (N2O): The dominant ozone-depleting substance emitted in the 21st century. Science 2009, 326, 123–125. [Google Scholar] [CrossRef] [PubMed]
  18. Tachikawa, H.; Abe, S. Reaction dynamics following electron capture of chlorofluorocarbon adsorbed on water cluster: A direct density functional theory molecular dynamics study. J. Chem. Phys. 2007, 126, 194310. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Tachikawa, H. Dissociative electron capture of halocarbon caused by the internal electron transfer from water trimer anion. Phys. Chem. Chem. Phys. 2008, 10, 2200–2206. [Google Scholar] [CrossRef] [PubMed]
  20. Tachikawa, H. Direct ab initio MD study on the electron capture dynamics of hydroperoxy radical (HOO)-water complexes. J. Phys. Chem. A 2010, 114, 4951–4956. [Google Scholar] [CrossRef] [PubMed]
  21. Lu, Q.B.; Sanche, L. Effects of cosmic rays on atmospheric chlorofluorocarbon dissociation and ozone depletion. Phys. Rev. Lett. 2001, 87, 078501. [Google Scholar] [CrossRef] [PubMed]
  22. Lu, Q.B. Correlation between cosmic rays and ozone depletion. Phys. Rev. Lett. 2009, 102, 118501. [Google Scholar] [CrossRef] [PubMed]
  23. Lu, Q.B. Cosmic-ray-driven electron-induced reactions of halogenated molecules adsorbed on ice surfaces: Implications for atmospheric ozone depletion and global climate change. Phys. Rep. 2010, 487, 141–167. [Google Scholar] [CrossRef] [Green Version]
  24. Xu, K.; Wang, W.; Wei, W.; Feng, W.; Sun, Q.; Li, P. Insights into the reaction mechanism of Criegee intermediate CH2OO with methane and implications for the formation of methanol. J. Phys. Chem. A 2017, 121, 7236–7245. [Google Scholar] [CrossRef] [PubMed]
  25. Guo, C.; Wang, W.; Feng, W.; Li, P. Insights into the one-electron reduction behavior of tetrachloro-o-benzoquinone: A DFT and molecular dynamics study. RSC Adv. 2017, 7, 12775–12782. [Google Scholar] [CrossRef]
  26. Li, P.; Guo, C.; Feng, W.; Sun, Q.; Wang, W. A DFT study on the reaction mechanism between tetrachloro-o-benzoquinone and H2O2 and an alternative reaction approach to produce the hydroxyl radical. RSC Adv. 2017, 7, 22919–22926. [Google Scholar] [CrossRef]
  27. Wang, W.; Zhang, X.; Li, P.; Sun, Q.; Li, Z.; Ren, C.; Guo, C. CO2 capture and separation from N2/CH4 mixtures by Co@B8/Co@B8– and M@B9/M@B9– (M = Ir, Rh, Ru) clusters: A theoretical study. J. Phys. Chem. A 2015, 119, 796–805. [Google Scholar] [CrossRef] [PubMed]
  28. Li, P.; Wang, W.; Sun, Q.; Li, Z.; Du, A.; Bi, S.; Zhao, Y. Insights into the mechanism of the reaction between tetrachloro-p-benzoquinone and hydrogen peroxide and their implications in the catalytic role of water molecules in producing the hydroxyl radical. Chem. Phys. Chem. 2013, 14, 2737–2743. [Google Scholar] [CrossRef] [PubMed]
  29. Li, P.; Ma, Z.; Wang, W.; Song, R.; Zhai, Y.; Bi, S.; Sun, H.; Bu, Y. Theoretical studies on the electron capture properties of the H2SO4…HOO● complex and its implications as an alternative source of HOOH. Phys. Chem. Phys. Chem. 2011, 13, 5931–5939. [Google Scholar] [CrossRef] [PubMed]
  30. Bader, R.F.W. Atoms in Molecules: A Quantum Theory; Oxford University Press: Oxford, UK, 1990. [Google Scholar]
  31. Reed, A.E.; Curtiss, L.A.; Weinhold, F. Intermolecular interactions from a natural bond orbital, donor-acceptor viewpoint. Chem. Rev. 1988, 88, 899–926. [Google Scholar] [CrossRef]
  32. Johnson, E.R.; Keinan, S.; Mori-Sánchez, P.; Contreras-García, J.; Cohen, A.J.; Yang, W. Revealing noncovalent interactions. J. Am. Chem. Soc. 2010, 132, 6498–6506. [Google Scholar] [CrossRef] [PubMed]
  33. Lu, T.; Chen, F. Multiwfn: A multifunctional wavefunction analyzer. J. Comput. Chem. 2012, 33, 580–592. [Google Scholar] [CrossRef] [PubMed]
  34. Boys, S.F.; Bernardi, F. The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errors. Mol. Phys. 1970, 19, 553–566. [Google Scholar] [CrossRef]
  35. Schlegel, H.B.; Millam, J.M.; Iyengar, S.S.; Voth, G.A.; Daniels, A.D.; Scuseria, G.E.; Frisch, M.J. Ab initio molecular dynamics: Propagating the density matrix with Gaussian orbitals. J. Chem. Phys. 2001, 114, 9758–9763. [Google Scholar] [CrossRef]
  36. Iyengar, S.S.; Schlegel, H.B.; Millam, J.M.; Voth, G.A.; Scuseria, G.E.; Frisch, M.J. Ab initio molecular dynamics: Propagating the density matrix with Gaussian orbitals. II. Generalizations based on mass-weighting, idempotency, energy conservation and choice of initial conditions. J. Chem. Phys. 2001, 115, 10291–10302. [Google Scholar] [CrossRef]
  37. Schlegel, H.B.; Iyengar, S.S.; Li, X.; Millam, J.M.; Voth, G.A.; Scuseria, G.E.; Frisch, M.J. Ab initio molecular dynamics: Propagating the density matrix with Gaussian orbitals. III. Comparison with Born-Oppenheimer dynamics. J. Chem. Phys. 2002, 117, 8694–8704. [Google Scholar] [CrossRef]
  38. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G.A. Gaussian, 9th ed.; Gaussian Inc.: Wallingford, CT, USA, 2009. [Google Scholar]
  39. Li, P.; Ma, Z.; Wang, W.; Shen, Z.; Bi, S.; Sun, H.; Bu, Y. Coupling interactions between sulfurous acid and the hydroperoxyl radical and implications for the formation of highly stable intermediates. ChemPhysChem 2010, 11, 696–705. [Google Scholar] [CrossRef] [PubMed]
  40. Torrent-Sucarrat, M.; Anglada, J.M. On the gas phase hydrogen bond complexes between formic acid and hydroperoxyl radical. A theoretical study. J. Phys. Chem. A 2006, 110, 9718–9726. [Google Scholar] [CrossRef] [PubMed]
  41. Fonseca Guerra, C.; Bickelhaupt, F.M.; Snijders, J.G.; Baerends, E.J. The nature of the hydrogen bond in DNA base pairs: The role of charge transfer and resonance assistance. Chem. Eur. J. 1999, 5, 3581–3594. [Google Scholar] [CrossRef]
Sample Availability: Samples of the compounds are not available from the authors.
Molecules 23 02349 g001 550
Figure 1. The optimized neutral complexes (A) and (B) and their molecular graphs, where the units of the bond distances and angles in parentheses are in angstroms and degrees, respectively. The bond critical point (BCP) is denoted as a small red dot.
Figure 1. The optimized neutral complexes (A) and (B) and their molecular graphs, where the units of the bond distances and angles in parentheses are in angstroms and degrees, respectively. The bond critical point (BCP) is denoted as a small red dot.
Molecules 23 02349 g001
Molecules 23 02349 g002 550
Figure 2. The scatter diagrams (left) and isosurfaces (right) of RDG analyses for complexes (A) and (B). The blue, green, and red colors represent the H-bond interactions, Van der Waals interactions, and steric effects, respectively.
Figure 2. The scatter diagrams (left) and isosurfaces (right) of RDG analyses for complexes (A) and (B). The blue, green, and red colors represent the H-bond interactions, Van der Waals interactions, and steric effects, respectively.
Molecules 23 02349 g002
Molecules 23 02349 g003 550
Figure 3. The NBO orbital interaction diagrams for the H-bonds in complexes A and B, where the isodensity contours are 0.001 electron/bohr3.
Figure 3. The NBO orbital interaction diagrams for the H-bonds in complexes A and B, where the isodensity contours are 0.001 electron/bohr3.
Molecules 23 02349 g003
Molecules 23 02349 g004 550
Figure 4. Calculated IR spectra of complexes A and B as well as H2SO4.
Figure 4. Calculated IR spectra of complexes A and B as well as H2SO4.
Molecules 23 02349 g004
Molecules 23 02349 g005 550
Figure 5. Electron density difference maps for complexes (A) and (B). Blue and violet areas denote the loss and gain of electron density, respectively. The isodensity contours are 0.001 electron/bohr3.
Figure 5. Electron density difference maps for complexes (A) and (B). Blue and violet areas denote the loss and gain of electron density, respectively. The isodensity contours are 0.001 electron/bohr3.
Molecules 23 02349 g005
Molecules 23 02349 g006 550
Figure 6. The dependence of ΔG on temperatures in the formation process of complexes A and B.
Figure 6. The dependence of ΔG on temperatures in the formation process of complexes A and B.
Molecules 23 02349 g006
Molecules 23 02349 g007 550
Figure 7. The time evolution of the potential energy (top) and the intermolecular H-bonds in complexes A and B (bottom) in molecular dynamics calculations.
Figure 7. The time evolution of the potential energy (top) and the intermolecular H-bonds in complexes A and B (bottom) in molecular dynamics calculations.
Molecules 23 02349 g007
Molecules 23 02349 g008 550
Figure 8. The optimized electron capture products and the related transition state.
Figure 8. The optimized electron capture products and the related transition state.
Molecules 23 02349 g008
Molecules 23 02349 g009 550
Figure 9. The molecular graphs for the electron capture products, where the BCP and ring critical point are denoted as small red and yellow dots, respectively.
Figure 9. The molecular graphs for the electron capture products, where the BCP and ring critical point are denoted as small red and yellow dots, respectively.
Molecules 23 02349 g009
Molecules 23 02349 g010 550
Figure 10. Spin density distribution diagram of the electron capture products, where the isodensity contours are 0.03 electron/bohr3.
Figure 10. Spin density distribution diagram of the electron capture products, where the isodensity contours are 0.03 electron/bohr3.
Molecules 23 02349 g010
Molecules 23 02349 g011 550
Figure 11. The changes of potential energy and main geometries in the electron capture process of complex A along with the simulation time.
Figure 11. The changes of potential energy and main geometries in the electron capture process of complex A along with the simulation time.
Molecules 23 02349 g011
Molecules 23 02349 g012 550
Figure 12. The changes of potential energy and main geometries for complex A′′ along with the simulation time.
Figure 12. The changes of potential energy and main geometries for complex A′′ along with the simulation time.
Molecules 23 02349 g012
Molecules 23 02349 g013 550
Figure 13. The changes of potential energy and main geometries in the electron capture process of complex B along with the simulation time.
Figure 13. The changes of potential energy and main geometries in the electron capture process of complex B along with the simulation time.
Molecules 23 02349 g013
Table
Table 1. The topological parameters at the BCPs for the intermolecular H-bonds a.
Table 1. The topological parameters at the BCPs for the intermolecular H-bonds a.
ComplexesBCPρbcp∇2ρbcpGbcpVbcpHbcp
AO4-H5···N100.01980.06940.0155−0.01370.0018
BO4-H5···O80.02300.08130.0190−0.01770.0013
A′O6-H7···O80.02010.06810.0156−0.01420.0014
N10-H5···O40.08160.08970.0550−0.0876−0.0326
A″O8-H7···O60.07270.10290.0523−0.0788−0.0266
O4-H5···N100.01470.04940.0103−0.00830.0020
B′O8-H5···O40.04240.10900.0340−0.0408−0.0068
a Atomic numbering refers to Figure 1 and Figure 8.
Table
Table 2. NBO analyses for the H-bonds in complexes A and B a.
Table 2. NBO analyses for the H-bonds in complexes A and B a.
ComplexesDonor NBOAcceptor NBOE(2)
ALP(N10)σ*(H5-O4)6.38
BLP1(O8)σ*(H5-O4)2.91
LP2(O8)σ*(H5-O4)4.79
a The unit of E(2) is in kcal/mol.
Table
Table 3. The calculated interaction energies (ΔEInt), relative energies (ΔERel), enthalpy and Gibbs free energy changes (ΔH and ΔG) during the formation process of complexes A and B a.
Table 3. The calculated interaction energies (ΔEInt), relative energies (ΔERel), enthalpy and Gibbs free energy changes (ΔH and ΔG) during the formation process of complexes A and B a.
ComplexesΔEIntΔERelΔHΔG
A−2.14(−2.72)0.76(1.72)−2.07(−2.34)4.19(3.92)
B−2.77(−4.22)0.00(0.00)−2.92(−3.93)4.58(3.57)
a The data in parentheses refer to the results at the CCSD(T)/AUG-cc-pVTZ level of theory plus the thermodynamic corrections at the B3LYP/6-311++G(3df,3pd) level of theory. All the units are in kcal/mol.
Table
Table 4. AEA, VEA, and VEDE results for complexes A and B a.
Table 4. AEA, VEA, and VEDE results for complexes A and B a.
ComplexesAEAVEAVEDE
A1.80−0.275.22
B3.18−0.277.44
a All the units are in eV.

Share and Cite

MDPI and ACS Style

Wang, W.-H.; Feng, W.-L.; Wang, W.-L.; Li, P. Theoretical Insights into the Electron Capture Behavior of H2SO4···N2O Complex: A DFT and Molecular Dynamics Study. Molecules 2018, 23, 2349. https://doi.org/10.3390/molecules23092349

AMA Style

Wang W-H, Feng W-L, Wang W-L, Li P. Theoretical Insights into the Electron Capture Behavior of H2SO4···N2O Complex: A DFT and Molecular Dynamics Study. Molecules. 2018; 23(9):2349. https://doi.org/10.3390/molecules23092349

Chicago/Turabian Style

Wang, Wei-Hua, Wen-Ling Feng, Wen-Liang Wang, and Ping Li. 2018. "Theoretical Insights into the Electron Capture Behavior of H2SO4···N2O Complex: A DFT and Molecular Dynamics Study" Molecules 23, no. 9: 2349. https://doi.org/10.3390/molecules23092349

APA Style

Wang, W. -H., Feng, W. -L., Wang, W. -L., & Li, P. (2018). Theoretical Insights into the Electron Capture Behavior of H2SO4···N2O Complex: A DFT and Molecular Dynamics Study. Molecules, 23(9), 2349. https://doi.org/10.3390/molecules23092349

Article Metrics

Back to TopTop