Oppenheimer Molecular Dynamics Research Articles

Reaction of SO3 with H2SO4 and its implications for aerosol particle formation in the gas phase and at the air–water interface

Abstract. The reactions between SO3 and atmospheric acids are indispensable in improving the formation of aerosol particles. However, relative to those of SO3 with organic acids, the reaction of SO3 with inorganic acids has not received much attention. Here, we explore the atmospheric reaction between SO3 and H2SO4, a typical inorganic acid, in the gas phase and at the air–water interface using quantum chemical (QC) calculations and Born–Oppenheimer molecular dynamics simulations. We also report the effect of H2S2O7, the product of the reaction between SO3 and H2SO4, on new particle formation (NPF) in various environments using the Atmospheric Cluster Dynamics Code (ACDC) kinetic model and QC calculations. The present findings show that the gas-phase reactions of SO3 + H2SO4 without and with water molecules are both low-energy-barrier processes. With the involvement of interfacial water molecules, H2O induced the formation of the S2O72-⋯H3O+ ion pair, HSO4- mediated the formation of the HSO4-⋯H3O+ ion pair, and the deprotonation of H2S2O7 was observed and proceeded on the picosecond timescale. The present findings suggest the potential contribution of the SO3–H2SO4 reaction to NPF and aerosol particle growth, showing that (i) although H2S2O7 is easily hydrolyzed with water to form H2SO4, it can directly participate in H2SO4–NH3-based cluster formation and can present a more obvious enhancement effect on SA–A-based cluster formation, and (ii) the formed interfacial S2O72- can attract candidate species from the gas phase to the water surface and, thus, accelerate particle growth.

Surface electronic structure and photo activity of double perovskite Ba[formula omitted]PrBiO[formula omitted]: First-principles investigations

Double perovskite Ba2PrBiO6 exhibits photo catalytic activity in the visible light range with the observed optical band gap of ∼1.0 eV. On the other hand, density functional theory (DFT) calculations predict a larger bulk band gap of 2.0 eV even with the standard Perdew–Burke–Ernzerhof (PBE) functional, which is known to underestimate the band gap of semiconductors. Much larger bulk band gap of 3.13 eV is predicted with a hybrid functional treatment for the exchange–correlation, which generally gives better descriptions for them. This discrepancy between experiment and theory implies the presence of unknown mechanism which effectively reduces the band gap of Ba2PrBiO6 and thereby gains the sensitivity to the visible light. The first finding is that the surface band gaps of Ba2PrBiO6 calculated using an appropriate hybrid functional are smaller than the bulk band gaps and are 2.49 eV and 2.85 eV for the Pr–Bi and Ba–O polar surfaces, respectively. We also find that work function is quite sensitive to the surface structures and its accurate calculations are essential to show that the reduction and oxidation levels are fit to the gap region. When a substitutional defect PrBi, in which Bi is substituted by Pr, is introduced into the Pr–Bi polar surface, surface band gap is further reduced to 0.93 eV, being very close to the experimental value. For the Ba–O polar surface, similar reduction to 1.59 eV is found but not substantial as compared to the Pr–Bi polar surface. These narrow surface band gaps may be the origin of the photo activity of Ba2PrBiO6 in the visible light region. We also performed Born–Oppenheimer molecular dynamics simulations to clarify the reactivity of the polar surfaces to water. We find that two hydroxyl groups were created in the adsorption of a single H2O molecule. The decomposition of water molecule in the adsorption indicates that the polar surfaces provide a promising stage for the photo catalytic activity. Further adsorption of H2O molecules leads to their diffusion on the surface.

The reaction mechanism of SO3 with the multifunctional compound ethanolamine and its atmospheric implications.

SO3 is an important reactive species in sulfur cycle and sulfuric acid formation processes and its reactions with some functional group substances, such as H2O, NH3, CH3OH, and organic and inorganic acids, have been extensively studied. However, its loss mechanism with multifunctional species is still lacking in detail. Herein, the reaction mechanism between SO3 and monoethanolamide (MEA) was investigated in the gas phase and on water droplets. The quantum chemical calculations indicate that the gas-phase reactions of SO3 with the OH and NH2 moieties of MEA hardly occur as their reaction energy barriers are up to 21.9-29.4 kcal mol-1. When a single water molecule is added into the SO3 + MEA reaction, it not only decreases the reaction barrier by at least 15.0 kcal mol-1 and thus enhances the rate obviously, but can also lead to the main product changing from HOCH2CH2NHSO3H to NH2CH2CH2OSO3H. The Born Oppenheimer molecular dynamics simulations on a water droplet show that the routes of the NH2CH2CH2OSO3-⋯H3O+ ion pair, HSO4- and HOCH2CH2NH3+ ions and zwitterionic formations of HOCH2CH2NH2+-SO3- and SO3--OCH2CH2NH3+ occur through a loop-structure route or chain reaction process, and can be finished within several picoseconds. Interestingly, the nucleation simulations show that the products of HOCH2CH2NHSO3H and NH2CH2CH2OSO3H have a potential ability to participate in the formation of new particles as they can form larger clusters with H2SO4, NH3 and H2O molecules within 20 ns. Thus, the present study will not only give new insight into the reaction between SO3 and multifunctional compounds, but also provide a new potential formation mechanism for particles resulting from the loss of SO3.

Planar tetracoordinate beryllium compounds with a partially covalent Be-Ng bond.

An analysis of the thermodynamic and kinetic stability and the nature of the chemical bond in hypercoordinated compounds with the formula BeH3Ng+ (Ng = He-Rn) through high-level calculations is presented in this work. Thermochemical calculations show that, for the heavier noble gases (Ar-Rn), these systems are thermodynamically stable at room temperature; however, this stability decreases due to a weakening of the Be-H2 interaction, while the Be-Ng bond strengthens going down the periodic table. These results are complemented by Born Oppenheimer molecular dynamics simulations, in which the increasing tendency to dissociate the Be-H2 bond is evidenced. The nature of the chemical bonding depends on the analysis performed. On the one hand, the interacting quantum atoms method indicates that the covalent contribution is around 25 to 30%. On the other hand, the electron density topology indicates a covalent nature for compounds with Kr-Rn, while Hirshfeld population analysis in conjunction with Mayer's bond order establishes polar covalent behavior. The geometrical parameters and natural energy decomposition analysis (NEDA) indicate a covalent nature, allowing us to consider that the Be-Ng bond has a partially covalent character.

Combining stochastic density functional theory with deep potential molecular dynamics to study warm dense matter

In traditional finite-temperature Kohn–Sham density functional theory (KSDFT), the partial occupation of a large number of high-energy KS eigenstates restricts the use of first-principles molecular dynamics methods at extremely high temperatures. However, stochastic density functional theory (SDFT) can overcome this limitation. Recently, SDFT and the related mixed stochastic–deterministic density functional theory, based on a plane-wave basis set, have been implemented in the first-principles electronic structure software ABACUS [Q. Liu and M. Chen, Phys. Rev. B 106, 125132 (2022)]. In this study, we combine SDFT with the Born–Oppenheimer molecular dynamics method to investigate systems with temperatures ranging from a few tens of eV to 1000 eV. Importantly, we train machine-learning-based interatomic models using the SDFT data and employ these deep potential models to simulate large-scale systems with long trajectories. Subsequently, we compute and analyze the structural properties, dynamic properties, and transport coefficients of warm dense matter.

Electron binding energies of SO2 at the surface of a water cluster

The electronic properties of SO2 at the surface of a water cluster were investigated by employing a combination of Born–Oppenheimer molecular dynamics and electron propagator theory (EPT). In our work, we utilized a revised version of the Perdew–Burke–Ernzerhof (PBE) exchange-correlation functional, which incorporates empirical corrections for dispersion interactions in line with a recent study of the air–water interface conducted by Ohto et al. [J. Phys. Chem. Lett. 10(17), 4914–4919 (2019)]. Polarization effects induce a significant broadening of the electron binding energy distribution, as predicted by EPT. This broadening can result in a substantial increase in electron affinity, impacting the chemical reactivity of SO2 at the air–water interface, a topic of significant and recent research interest. We discuss the relationship between electron binding energies (EBEs) and the specific connections of SO2 to water. The results indicate that configurations involving an OS⋯H bond tend to yield higher electron affinities compared to complex formation through S⋯OW bonds. Surprisingly, SO2 molecules not bound to water molecules according to a specific criterion may also exhibit higher electron affinities. This feature can be explained by the role played by the polarization field from water molecules. Our best estimate for the HOMO–LUMO (H–L) gap of SO2 at the surface of a water cluster is 11.6 eV. Very similar H–L gaps are predicted for isolated and micro-solvated SO2. Fukui functions for the gas phase, and the micro-solvated SO2–H2O complex supports the view that the LUMO is predominantly localized on the SO2 moiety.

First results from the interaction between hydroxyl radicals and condensed water droplets: A spin-polarized DFT-MD study

The formation of clouds results from condensation of water vapor into liquid water droplets. These water droplets form on the surface of tiny particles that are suspended in the air, like dust. Atmospheric chemistry relies heavily on the hydroxyl radical (OH*), which is highly reactive and prone to combining with hydrogen atoms in other compounds. Despite previous studies concerning liquid, solid, and gas phase interactions, how hydroxyl radicals interact with condensed water droplets is still unknown. OH*-condensed droplet interactions show unique behavior, according to our studies. Several hydrogen transfers were detected during the interaction process. It takes approximately 3.5 kcal/mol of energy to initiate hydrogen transfer. Moreover, we revealed the dominant structures of hydrated radicals during the mobility of OH* into condensed droplets using Born–Oppenheimer molecular dynamics simulations. In contrast to our previous studies, the hemibond doesn't significantly affect the stability of OH*(H2O)n structures in condensed droplets. Our results showed that condensed droplets hinder the mobility of OH*, causing its diffusion coefficient to decrease. Moreover, several analyses were carried out to investigate condensed droplet surfaces, bonding patterns, and the position of OH* in conjunction with water molecules, including mobility analyses, radial distribution functions, coordination analyses, and population analyses.

Possible Roles of Transition Metal Cations in the Formation of Interstellar Benzene via Catalytic Acetylene Cyclotrimerization.

Polycyclic aromatic hydrocarbons (PAHs) are ubiquitous interstellar molecules. However, the formation mechanisms of PAHs and even the simplest cyclic aromatic hydrocarbon, benzene, are not yet fully understood. Recently, we reported the statistical and dynamical properties in the reaction mechanism of Fe+-catalyzed acetylene cyclotrimerization, whereby three acetylene molecules are directly converted to benzene. In this study, we extended our previous work and explored the possible role of the complex of other 3d transition metal cations, TM+ (TM = Sc, Ti, Mn, Co, and Ni), as a catalyst in acetylene cyclotrimerization. Potential energy profiles for bare TM+-catalyst (TM = Sc and Ti), for TM+NC--catalyst (TM = Sc, Ti, Mn, Co, and Ni), and for TM+-(H2O)8-catalyst (TM = Sc and Ti) systems were obtained using quantum chemistry calculations, including the density functional theory levels. The calculation results show that the scandium and titanium cations act as efficient catalysts in acetylene cyclotrimerization and that reactants, which contain an isolated acetylene and (C2H2)2 bound to a bare (ligated) TM cation (TM = Sc and Ti), can be converted into a benzene-metal-cation product complex without an entrance barrier. We found that the number of electrons in the 3d orbitals of the transition metal cation significantly contributes to the catalytic efficiency in the acetylene cyclotrimerization process. On-the-fly Born-Oppenheimer molecular dynamics (BOMD) simulations of the Ti+-NC- and Ti+-(H2O)8 complexes were also performed to comprehensively understand the nuclear dynamics of the reactions. The computational results suggest that interstellar benzene can be produced via acetylene cyclotrimerization reactions catalyzed by transition metal cation complexes.

Born–Oppenheimer molecular dynamics study on collective protein dynamics invoked by ultrafast photoisomerization of retinal chromophore in bacteriorhodopsin

Large-scale quantum-chemistry-based Born–Oppenheimer molecular dynamics simulations of collective protein dynamics of bacteriorhodopsin initiated by photoinduced all-trans→ 13-cis isomerization of retinal Schiff base chromophore compensate for missing information from time-resolved serial femtosecond X-ray crystallography measurements in terms of electronic states, energy transfer, and details of structural dynamics.

Reversible Hydrogen Storage by Planar Hypercoordinate Carbon Clusters

Planar hypercoordinate carbon clusters show a high capacity of hydrogen storage with better stability, as revealed by using density functional theory. The interaction energy of H2 with the cluster lies in between chemisorption and physisorption (0.11–0.19 eV). Topological analysis indicates closed shell-type interaction of H2 with the clusters. The clusters remain stable before and after adsorption. The gravimetric densities are found to be 25.6, 39.7, and 13.5 wt % for planar tetracoordinate, pentacoordinate, and hexacoordinate clusters, respectively. Furthermore, Born–Oppenheimer molecular dynamics simulations indicate the reversibility of adsorption. The clusters can accommodate a large number of H2 molecules and can be released at an elevated temperature (400 K). Our computation reveals that the planar clusters may be used as ideal 2D materials for energy storage.

Hydrogen Storage in a Ti-Functionalized Metal Carbyne Framework: Insights from a First-Principles Study

Hydrogen (H2) has the ability to be a future green-energy source for the automotive industry due to its versatility, clean nature, and efficiency as an energy carrier. The biggest obstacle to achieving hydrogen as a potential energy source is the nonavailability of effective hydrogen storage materials. In this proposed work, the hydrogen storage capacity and reversibility of a Ti-functionalized metal carbyne framework (MCF) having carbyne as an organic linker and Mg4O as an inorganic moiety are explored. All structural relaxation and simulation calculations are performed by employing the dispersion-corrected density functional theory computation. In this Ti-functionalized MCF, each Ti can adsorb up to five H2 molecules via the Kubas mechanism, with adsorption energy lying between 0.33 and 0.44 eV, which falls within the range outlined by the U.S. Department of Energy (DOE). For a detailed study of the thermal stability of framework, hydrogen reversibility, and mechanism of hydrogen sorption, Born–Oppenheimer molecular dynamics (BOMD) simulations are carried out at various temperatures. From the findings of BOMD simulations, the Ti-functionalized MCF can reversibly adsorb the H2 molecules with a high gravimetric density of about 9.4 wt %. This proposed Ti-functionalized MCF stands out as a potential hydrogen storage medium meeting the standards established by the DOE.

Shadow energy functionals and potentials in Born-Oppenheimer molecular dynamics.

In Born-Oppenheimer molecular dynamics (BOMD) simulations based on the density functional theory (DFT), the potential energy and the interatomic forces are calculated from an electronic ground state density that is determined by an iterative self-consistent field optimization procedure, which, in practice, never is fully converged. The calculated energies and forces are, therefore, only approximate, which may lead to an unphysical energy drift and instabilities. Here, we discuss an alternative shadow BOMD approach that is based on backward error analysis. Instead of calculating approximate solutions for an underlying exact regular Born-Oppenheimer potential, we do the opposite. Instead, we calculate the exact electron density, energies, and forces, but for an underlying approximate shadow Born-Oppenheimer potential energy surface. In this way, the calculated forces are conservative with respect to the approximate shadow potential and generate accurate molecular trajectories with long-term energy stabilities. We show how such shadow Born-Oppenheimer potentials can be constructed at different levels of accuracy as a function of the integration time step, δt, from the constrained minimization of a sequence of systematically improvable, but approximate, shadow energy density functionals. For each energy functional, there is a corresponding ground state Born-Oppenheimer potential. These pairs of shadow energy functionals and potentials are higher-level generalizations of the original "zeroth-level" shadow energy functionals and potentials used in extended Lagrangian BOMD [Niklasson, Eur. Phys. J. B 94, 164 (2021)]. The proposed shadow energy functionals and potentials are useful only within this extended dynamical framework, where also the electronic degrees of freedom are propagated as dynamical field variables together with the atomic positions and velocities. The theory is quite general and can be applied to MD simulations using approximate DFT, Hartree-Fock, or semi-empirical methods, as well as to coarse-grained flexible charge models.

Geometric structures and hydrogen storage properties of M2B7 (M=Be, Mg, Ca) clusters

Based on the density functional theory, we investigate the electronic properties of the clusters M2B7 (M = Be, Mg, Ca) and their hydrogen storage properties systematically in this paper. Extensive global search results show that the global minimal structures of the three systems (Be2B7, Mg2B7 and Ca2B7) are heptagonal biconical structure, and the two alkaline earth metals are located at the top of the biconical. Chemical bonding analyses show that M2B7 clusters have 6σ and 6π delocalized electrons, which are doubly aromatic. At the wB97XD level, the three systems have good hydrogen storage capabilities. The hydrogen storage density of Be2B7 is as high as 23.03 wt%, while Mg2B7 and Ca2B7 also far exceed the hydrogen storage target set by the U.S. Department of Energy in 2017. Their average adsorption energies of H2 molecules all ranged from 0.1 eV/H2 to 0.48 eV/H2, which is fall in between physisorption and chemisorption. Extensive Born Oppenheimer molecular dynamics (BOMD) simulations show that the H2 molecules of the three systems can be completely released at a certain temperature. Therefore, M2B7 systems can achieve reversible adsorption of H2 molecules at normal temperature and pressure. It can be seen that the B7 clusters modified by alkaline earth metals may become a promising new nano-hydrogen storage material.

CAl4X4 (X = Te, Po): Double Aromatic Molecular Stars Containing Planar Tetracoordinate Carbon Atoms

Planar tetracoordinate carbon (ptC) species are scarce and exotic. Introducing four peripheral Te/Po auxiliary atoms is an effective strategy to flatten the tetrahedral structure of CAl4 (Td, 1A1). Neutral CAl4X4 (X = Te, Po) clusters possess quadrangular star structures containing perfect ptC centers. Unbiased density functional theory (DFT) searches and high-level CCSD(T) calculations suggest that these ptC species are the global minima on the potential energy surfaces. Bonding analyses indicate that 40 valence-electron (VE) is ideal for the ptC CAl4X4 (X = Te, Po): one delocalized π and three σ bonds for the CAl4 core; four lone pairs (LPs) of four X atoms, eight localized Al-X σ bonds, and four delocalized Al-X-Al π bonds for the periphery. Thus, the ptC CAl4X4 (X = Te, Po) clusters possess the stable eight electron structures and 2π + 6σ double aromaticity. Born-Oppenheimer molecular dynamics (BOMD) simulations indicate that neutral ptC CAl4X4 (X = Te, Po) clusters are robust.

Explainable Artificial Intelligence Approach to Identify the Origin of Phonon‐Assisted Emission in WSe2 Monolayer

The application of explainable artificial intelligence in nanomaterial research has emerged in the past few years, which has facilitated the discovery of novel physical findings. However, a fundamental question arises concerning the physical insights presented by deep neural networks; the model interpretation results have not been carefully evaluated. Herein, explainable artificial intelligence and quantum mechanical calculations is bridged to investigate the correlation between light scattering and emission in a WSe2 monolayer. Convolutional neural networks using light scattering and emission data are first trained, while expecting the networks to determine the relationships between them. The trained models are interpreted and the specific phonon contribution during the exciton relaxation process is derived. Finally, the findings are independently evaluated through quantum mechanical calculations, such as the Born–Oppenheimer molecular dynamics simulation and density functional perturbation theory. The study provides reliable fundamental physical insight by evaluating the results of neural networks and suggests a novel methodology that can be applied in materials science.

Dehydroxylation of smectites from Nudged Elastic Band and Born-Oppenheimer Molecular Dynamics simulations

High–temperature treatment of smectite minerals leads to the dehydroxylation process. One of the consequences of this process is the loss of sorption ability, which is an interesting, but not fully understood, phenomenon. A combination of static and dynamic quantum–chemical methods was used to analyze structural changes associated with the dehydroxylation of smectite layers. The work aims in the determination of multi–path reaction mechanism. The Nudged Elastic Band (NEB) method allowed for the recognition of intermediate structures, transition states and activation energies in the studied reaction mechanisms. Born–Oppenheimer Molecular Dynamics (BOMD) was used to investigate the interatomic distances time–evolution of atoms involved in the dehydroxylation process at different temperatures. The simulations were carried out in the crystalline phase, for which the computational model was prepared based on the dioctahedral structure of smectites with the chemical formula: Ca2+Al4Si8O22(OH)2. It was found, that the process requires a proton transfer to the oxygen atom of the hydroxyl group located in the octahedral sheet, which results in the release of water molecule from the layer. In addition, it has been shown that the process of a water molecule formation on the surface is also possible. The process occurring inside the layer is more energetically favoured than that on the surface. However, high activation energies were obtained in both cases, what proves that the dehydroxylation reaction takes place only at a relatively high temperatures. It has also been demonstrated that this process can follow several different mechanisms.

Graph-based quantum response theory and shadow Born-Oppenheimer molecular dynamics.

Graph-based linear scaling electronic structure theory for quantum-mechanical molecular dynamics simulations [A. M. N. Niklasson et al., J. Chem. Phys. 144, 234101 (2016)] is adapted to the most recent shadow potential formulations of extended Lagrangian Born-Oppenheimer molecular dynamics, including fractional molecular-orbital occupation numbers [A. M. N. Niklasson, J. Chem. Phys. 152, 104103 (2020) and A. M. N. Niklasson, Eur. Phys. J. B 94, 164 (2021)], which enables stable simulations of sensitive complex chemical systems with unsteady charge solutions. The proposed formulation includes a preconditioned Krylov subspace approximation for the integration of the extended electronic degrees of freedom, which requires quantum response calculations for electronic states with fractional occupation numbers. For the response calculations, we introduce a graph-based canonical quantum perturbation theory that can be performed with the same natural parallelism and linear scaling complexity as the graph-based electronic structure calculations for the unperturbed ground state. The proposed techniques are particularly well-suited for semi-empirical electronic structure theory, and the methods are demonstrated using self-consistent charge density-functional tight-binding theory both for the acceleration of self-consistent field calculations and for quantum-mechanical molecular dynamics simulations. Graph-based techniques combined with the semi-empirical theory enable stable simulations of large, complex chemical systems, including tens-of-thousands of atoms.

Entropic Effects on the Aqueous Microsolvation of Protonated Glycine and Protonated β-Alanine. Hybrid Density Functional Theory Born-Oppenheimer Molecular Dynamics Studies.

Recent low-temperature infrared-based experimental studies provided information about the effects of aqueous microsolvation on the intramolecular hydrogen bond of protonated glycine and β-alanine [J. Phys. Chem. A 2019, 123, 3355]. Here we address the temperature-dependent entropic effects on the aqueous microsolvation patterns of these protonated amino acids using the AAH+(H2O)n (n = 1-8) cluster model at 50 K and room temperature with Born-Oppenheimer molecular dynamics using a calibrated hybrid density functional. The CCOOH-Ow, N-Ow, and center-of-mass-Ow radial distribution functions provide accurate structural data and temperature-dependent water coordination numbers vs. solvation degree. The solvation patterns for protonated glycine at 50 K show structural features in agreement with previous static optimizations. However, entropic effects at room temperature play a crucial role in the evolution of the intramolecular HB strength vs. solvation degree for both protonated amino acids. With increasing hydration entropic effects favor the making of solvent hydrogen bond networks over full solvation of protonated glycine. At room temperature four water molecules are needed to build the first solvation shell for protonated glycine while five are required for protonated β-alanine. A new statistical Cumulative Percentage of Structures (CPS) scheme is proposed; when the CPS data are analyzed in light of the empirical formula of Rozenberg et al. [Phys. Chem. Chem. Phys. 2000, 2, 2699] and the hydrogen bond relative strength (HBRS) criteria of Jeffrey [An Introduction to Hydrogen Bonding; Oxford University: 1997] we can provide a detailed molecular mechanism for the weakening of the intramolecular hydrogen bond based on the average dynamical structures, which clearly reveals the temperature dependence of this process. The new CPS-HBRS scheme proposed here can be utilized using any type of molecular dynamics trajectory (classical, BOMD, CPMD, etc.).

Indication of 310-Helix Structure in Gas-Phase Neutral Pentaalanine.

We investigate the gas-phase structure of the neutral pentaalanine peptide. The IR spectrum in the 340-1820 cm-1 frequency range is obtained by employing supersonic jet cooling, infrared multiphoton dissociation, and vacuum-ultraviolet action spectroscopy. Comparison with quantum chemical spectral calculations suggests that the molecule assumes multiple stable conformations, mainly of two structure types. In the most stable conformation theoretically found, the N-terminus forms a C5 ring and the backbone resembles that of an 310-helix with two β-turns. Additionally, the conformational preferences of pentaalanine have been evaluated using Born-Oppenheimer molecular dynamics, showing that a nonzero simulation time step causes a systematic frequency shift.

Significant influence of water molecules on the SO3 + HCl reaction in the gas phase and at the air-water interface.

The products resulting from the reactions between atmospheric acids and SO3 have a catalytic effect on the formation of new particles in aerosols. However, the SO3 + HCl reaction in the gas-phase and at the air-water interface has not been considered. Herein, this reaction was explored exhaustively by using high-level quantum chemical calculations and Born Oppenheimer molecular dynamics (BOMD) simulations. The quantum calculations show that the gas-phase reaction of SO3 + HCl is highly unlikely to occur under atmospheric conditions with a high energy barrier of 22.6 kcal mol-1. H2O and (H2O)2 play obvious catalytic roles in reducing the energy barrier of the SO3 + HCl reaction by over 18.2 kcal mol-1. The atmospheric lifetimes of SO3 show that the (H2O)2-assisted reaction dominates over the H2O-assisted reaction within the altitude range of 0-5 km, whereas the H2O-assisted reaction is more favorable within an altitude range of 10-50 km. BOMD simulations show that H2O-induced formation of the ClSO3-⋯H3O+ ion pair and HCl-assisted formation of the HSO4-⋯H3O+ ion pair were identified at the air-water interface. These routes followed a stepwise reaction mechanism and proceeded at a picosecond time scale. Interestingly, the formed ClSO3H in the gas phase has a tendency to aggregate with sulfuric acids, ammonias, and water molecules to form stable clusters within 40 ns simulation time, while the interfacial ClSO3- and H3O+ can attract H2SO4, NH3, and HNO3 for particle formation from the gas phase to the water surface. Thus, this work will not only help in understanding the SO3 + HCl reaction driven by water molecules in the gas-phase and at the air-water interface, but it will also provide some potential routes of aerosol formation from the reaction between SO3 and inorganic acids.