Spectroscopy and Photochemistry of Aluminum-Bearing Species in the Universe

ConspectusMetal-bearing molecules impact the chemical and physical environment of many astronomical sources such as the circumstellar envelopes of large asymptotic giant branch and red supergiant stars, the interstellar medium, and planetary atmospheres (e.g., ablation of ∼20 tons per day into the Earth's upper atmosphere). In recent decades, the number of successfully detected metal-containing molecules has increased via rotational spectroscopic observations, which are driven by theoretical and experimental investigations. Following formation, the ultimate fate of each species (stabilization, dissociation, etc.) is determined by its electronic structure and electronic spectroscopic properties as it encounters the pervasive radiation fields in the vacuum of space. Studying these properties can evince the possibility of detection and predict the impact each molecule has on its surrounding environment. Aluminum, one of the most abundant elements and metals, is distributed throughout the universe as a constituent of gas-phase molecules (e.g., AlO, AlOH, AlCl, etc.) as well as condensed onto solid dust grains such as Al2O3. Free gas-phase aluminum-bearing molecules are synthesized by nonthermal equilibrium processes such as shocks and pulsations near the stellar photosphere or via the reaction of molecules on the surface of dust grains. Recent investigations in our research group utilizing quantum chemical methods, such as coupled cluster (CCSD(T) and CCSD(T)-F12) and multireference configuration interaction (MRCI) with large basis sets, have explored a wide breadth of spectroscopy and photochemistry of small (triatomic and tetratomic) aluminum-bearing molecules, including Al-H, Al-C, Al-N, Al-O, Al-Si, Al-P, and Al-S bonds, among others. The ground-state spectroscopy (rotational and vibrational) of various aluminum-bearing molecules is discussed in the context of experimental and observational detection potentials. These detection potentials depend on various factors, such as the magnitude of the permanent dipole moment (PDM) and the population of states yielding transition frequencies in detectable ranges. Many aluminum-bearing molecules possess large PDMs and may be prime candidates for astronomical and laboratory detection. Within this discussion, interesting aspects of the ground-state molecular orbital configuration of OAlNO are shown to lead to an uncommon triplet ground state. Additionally, the electronic absorption spectrum of the quasi-isoenergetic ground-state isomers of AlOSO is discussed as a sensitive method for detecting this species and differentiating between the two isomers. Finally, photochemical mechanisms key to the production of AlO and AlOH in low-density regions and the destruction of AlCO and AlOC are also discussed in order to understand the radiation-induced formation and destruction of these molecules.

Accelerated Photolysis of H2O2 at the Air-Water Interface of a Microdroplet

Photochemical homolysis of hydrogen peroxide (H2O2) occurs widely in nature and is a key source of hydroxyl radicals (·OH). The kinetics of H2O2 photolysis play a pivotal role in determining the efficiency of ·OH production, which is currently mainly investigated in bulk systems. Here, we report considerably accelerated H2O2 photolysis at the air-water interface of microdroplets, with a rate 1.9 × 103 times faster than that in bulk water. Our simulations show that due to the trans quasiplanar conformational preference of H2O2 at the air-water interface compared to the bulk or gas phase, the absorption peak in the spectrum of H2O2 is significantly redshifted by 45 nm, corresponding to greater absorbance of photons in the sunlight spectrum and faster photolysis of H2O2. This discovery has great potential to solve current problems associated with ·OH-centered heterogeneous photochemical processes in aerosols. For instance, we show that accelerated H2O2 photolysis in microdroplets could lead to markedly enhanced oxidation of SO2 and volatile organic compounds.

Distinctive Heterogeneous Reaction Mechanism of ClNO2 on the Air-Water Surface Containing Cl

The heterogeneous reaction of nitryl chloride (ClNO2) on the air-water surface plays a significant role in the chloride lifecycle. The air-water surface is ubiquitous on ice surfaces under supercooled conditions, affecting the uptake and heterogeneous reaction processes of trace gases. Previous studies suggest that ClNO2 is formed on Cl-doped ice surfaces following the N2O5 uptake. Herein, a distinctive heterogeneous reaction mechanism of ClNO2 is suggested on an air-water surface containing Cl under supercooled conditions using combined classic molecular dynamics (MD) and Born-Oppenheimer MD simulations. It is found that N2O5 dissociates into a NO2+ and NO3- ionic pair on the top air-water surface. In the top layer of the surface containing barely any Cl-, NO2+ proceeds through hydrolysis and produces H3O+ and HNO3. Thus, surface acidification appears because of H3O+ yields. With NO2+ diffusion to the deep layer of the surface, NO2+ reacts with Cl- and forms ClNO2. Note that ClNO2 formation competes with NO2+ hydrolysis, and the rate of ClNO2 formation is 27.7[Cl-] larger than that of NO2+ hydrolysis. Afterward, the reaction of ClNO2 with Cl- becomes barrierless with the catalysis by H3O+, which is not feasible on a neutral surface. Cl2 is thus generated and escapes into the atmosphere (low solubility of Cl2), contributing to the Cl radical. The proposed mechanism bolsters the current understanding of ClNO2's fate and its role in Cl chemistry in extremely cold environments like the Arctic and other high-latitude regions in wintertime.

Spectroscopy and Photochemistry of OAlNO and Implications for New Metal Chemistry in the Atmosphere

A new aluminum-bearing species, OAlNO, which has the potential to impact the chemistry of the Earth's upper atmosphere, is characterized via high-level, ab initio, spectroscopic methods. Meteor-ablated aluminum atoms are quickly oxidized to aluminum oxide (AlO) in the mesosphere and lower thermosphere (MLT), where a steady-state layer of AlO then builds up. Concurrent formation of nitric oxide (NO) in the same region of the atmosphere will lead to the bimolecular formation of the OAlNO molecule. Molecular orbital analysis provides fundamental insights into the chemical bonding and energetic arrangement of the triplet (1 3A″) ground state and singlet (1 1A') excited-state species of OAlNO. Additionally, unpaired electrons on the terminal oxygen atom of triplet (1 3A″) OAlNO cause it to be reactive to atmospheric species, potentially impacting climate science and high-altitude chemistry. The triplet (1 3A″) ground-state species exhibits a large permanent dipole moment useful for rotational spectroscopic detection; however, similar rotational constants to the singlet (1 1A') excited-state species will hamper differentiation in a spectrum. Strong infrared intensities will assist in detection and discrimination of the different spin states and isomers. Repulsive electronic excited states of OAlNO will lead to photolysis of the Al-N bond and formation of various electronic states of AlO + NO through nonadiabatic pathways. Reaction through the OAlNO intermediate represents a means for the production of electronically excited AlO, leading to new chemistry in the atmosphere. Excitation to higher-lying electronic states will lead to fluorescence with a minor Stokes shift, useful for laboratory investigation. Such physical properties of this molecule will allow for new, unexplored chemical pathways in the MLT to be considered.

Molecular Insights into the Spontaneous Generation of Cl2O in the Reaction of ClONO2 and HOCl at the Air-Water Interface

Chemical processes involving chlorine nitrate (ClONO2) at the surface of stratospheric aerosols are crucial to ozone depletion. Herein, we show a reaction route for the formation of Cl2O, which is a source of stratospheric chlorine, in the ClONO2 + HOCl reaction at the air-water interface. Our ab initio molecular dynamics (AIMD) simulations show that the (ClONO2)Cl···O(HOCl) halogen bond plays a key role in the reaction and is the main interaction between ClONO2 and HOCl both at the air-water interface and in the bulk liquid water. Furthermore, metadynamics-based AIMD simulations reveal two pathways: (i) The OCl fragment of HOCl binds to the Cl atom in ClONO2, resulting in the formation of Cl2O and NO3-. Simultaneously, the remaining hydrogen atom is transferred to a water molecule to form H3O+. (ii) HOCl acts as a bridge for Cl atom transfer from ClONO2 to the O atom of a water molecule, and this water molecule transfers one of its H atoms to another water molecule, forming two HOCl molecules, NO3-, and H3O+. Free-energy calculations show that the former is the energetically more favorable process. More importantly, the free-energy barrier for Cl2O formation at the air-water interface is only ∼0.8 kcal/mol, and the reaction is exothermic. These findings provide insights into the importance of fundamental chlorine chemistry and the broader implications of the aerosol air-water interface for atmospheric chemistry.

Resolving the Formation Mechanism of HONO via Ammonia-Promoted Photosensitized Conversion of Monomeric NO2 on Urban Glass Surfaces

Understanding the formation processes of nitrous acid (HONO) is crucial due to its role as a primary source of hydroxyl radicals (OH) in the urban atmosphere and its involvement in haze events. In this study, we propose a new pathway for HONO formation via the UVA-light-promoted photosensitized conversion of nitrogen dioxide (NO₂) in the presence of ammonia (NH₃) and polycyclic aromatic hydrocarbons (PAHs, common compounds in urban grime). This new mechanism differs from the traditional mechanism, as it does not require the formation of the NO₂ dimer. Instead, the enhanced electronic interaction between the UVA-light excited triplet state of PAHs and NO₂-H₂O/NO₂-NH₃-H₂O significantly reduces the energy barrier and facilitates the exothermic formation of HONO from monomeric NO₂. Furthermore, the performed experiments confirmed our theoretical findings and revealed that the synergistic effect from light-excited PAHs and NH₃ boosts the HONO formation with determined HONO fluxes of 3.6 × 10¹⁰ molecules cm^-2 s^-1 at 60% relative humidity (RH) higher than any previously reported HONO fluxes. Intriguingly, light-induced NO₂ to HONO conversion yield on authentic urban grime in presence of NH₃ is unprecedented 130% at 60% RH due to the role of NH₃ acting as a hydrogen carrier, facilitating the transfer of hydrogen from H₂O to NO₂. These results show that NH₃-assisted UVA-light-induced NO₂ to HONO conversion on urban surfaces can be a dominant source of HONO in the metropolitan area.

Imaging of pH distribution inside individual microdroplet by stimulated Raman microscopy

Aerosol microdroplets as microreactors for many important atmospheric reactions are ubiquitous in the atmosphere. pH largely regulates the chemical processes within them; however, how pH and chemical species spatially distribute within an atmospheric microdroplet is still under intense debate. The challenge is to measure pH distribution within a tiny volume without affecting the chemical species distribution. We demonstrate a method based on stimulated Raman scattering microscopy to visualize the three-dimensional pH distribution inside single microdroplets of varying sizes. We find that the surface of all microdroplets is more acidic, and a monotonic trend of pH decreasing is observed in the 2.9-μm aerosol microdroplet from center to edge, which is well supported by molecular dynamics simulation. However, bigger cloud microdroplet differs from small aerosol for pH distribution. This size-dependent pH distribution in microdroplets can be related to the surface-to-volume ratio. This work presents noncontact measurement and chemical imaging of pH distribution in microdroplets, filling the gap in our understanding of spatial pH in atmospheric aerosol.

Spectroscopic characterization of [H, Cl, S, O] molecular system: Potential candidate for detection in Venus atmosphere

Accurate spectroscopic parameters have been obtained for the identification of the [H, Cl, S, O] molecular system in the Venus atmosphere using computational methods. These calculations employed both standard and explicitly correlated coupled cluster techniques. All isomers possess C1 symmetry, with HOSCl being the most stable isomer. Only HOSCl and trigonal-HSOCl isomers are thermodynamically stable relative to the first dissociation limit HCl + SO. Fundamental modes of the lowest three isomers exhibit many anharmonic resonances, resulting in complex spectra. All isomers are found to be stable in the visible region as the calculation of vertical energy transition indicates. No electronic states were found to strongly absorb in the near UV-vis region.

Chemical Implications of Rapid Reactive Absorption of I2O4 at the Air-Water Interface

Marine aerosol formation involving iodine-bearing species significantly affects the global climate and radiation balance. Although recent studies outline the critical role of iodine oxide in nucleation, much less is known about its contribution to aerosol growth. This paper presents molecular-level evidence that the air-water interfacial reaction of I₂O₄ mediated by potent atmospheric chemicals, such as sulfuric acid (H₂SO₄) and amines [e.g., dimethylamine (DMA) and trimethylamine (TMA)], can occur rapidly on a picosecond time scale by Born-Oppenheimer molecular dynamics simulations. The interfacial water bridges the reactants while facilitating the DMA-mediated proton transfer and stabilizing the ionic products of H₂SO₄-involved reactions. The identified heterogeneous mechanisms exhibit the dual contribution to aerosol growth: (i) the ionic products (e.g., IO₃^-, DMAH⁺, TMAH⁺, and HSO₄^-) formed by reactive adsorption possess less volatility than the reactants and (ii) these ions, such as alkylammonium salts (e.g., DMAH⁺), are also highly hydrophilic, further facilitating hygroscopic growth. This investigation enhances not only our understanding of heterogeneous iodine chemistry but also the impact of iodine oxide on aerosol growth. Also, these findings can bridge the gap between the abundance of I₂O₄ in the laboratory and its absence in field-collected aerosols and provide an explanation for the missing source of IO₃^-, HSO₄^-, and DMAH⁺ in marine aerosols.

Spontaneous Oxidation of Thiols and Thioether at the Air-Water Interface of a Sea Spray Microdroplet

The transport of dissolved organic sulfur, including thiols and thioethers, from the ocean surface to the atmosphere through sea spray aerosol (SSA) is of great importance for the global sulfur cycle. Thiol/thioether in SSA undergoes rapid oxidation that is historically linked to photochemical processes. Here, we report the discovery of a non-photochemical, spontaneous path of thiol/thioether oxidation in SSA. Among 10 investigated naturally abundant thiol/thioether, seven species displayed rapid oxidation in SSA, with disulfide, sulfoxide, and sulfone comprising the major products. We suggest that such spontaneous oxidation of thiol/thioether was mainly fueled by thiol/thioether enrichment at the air-water interface and generation of highly reactive radicals by the loss of an electron from ions (e.g., glutathionyl radical produced from ionization of deprotonated glutathione) at or near the surface of the water microdroplet. Our work sheds light on a ubiquitous but previously overlooked pathway of thiol/thioether oxidation, which could contribute to an accelerated sulfur cycle as well as related metal transformation (e.g., mercury) at ocean-atmosphere interfaces.

Photocatalytic Oxidation of NO2 on TiO2: Evidence of a New Source of N2O5

N2O5 is an important intermediate in the atmospheric nitrogen cycle. Using a flow tube reactor, N2O5 was found to be released from the TiO2 surface during the photocatalytic oxidation of NO2, revealing a previously unreported source of N2O5. The rate of N2O5 release from TiO2 was dependent on the initial NO2 concentration, relative humidity, O2/N2 ratio, and irradiation intensity. Experimental evidences and quantum chemical calculations showed that NO2 can react with the surface hydroxyl groups and the generated electron holes on the TiO2, followed by combining with another NO2 molecule to form N2O5. The latter was physisorbed on TiO2 and had a low adsorption energy of -0.13 eV. Box model simulations indicated that the new source of N2O5 released from TiO2 can increase the daytime N2O5 concentration by up to 20% in urban areas if abundant TiO2-containing materials and high NOx concentrations were present. This joint experimental/theoretical study not only demonstrates a new chemical mechanism for N2O5 formation but also has important implications for air quality in urban areas.

Evidence of Formation of 1-10 nm Diameter Ice Nanotubes in Double-Walled Carbon Nanotube Capillaries

Water exhibits rich phase behaviors in nanoscale confinement. Since the simulation evidence of the formation of single-walled ice nanotubes (INTs) in single-walled carbon nanotubes was confirmed experimentally, INTs have been recognized as a form of low-dimensional hydrogen-bonding network. However, the single-walled INTs reported in the literature all possess subnanometer diameters (<1 nm). Herein, based on systematic and large-scale molecular dynamics simulations, we demonstrate the spontaneous freezing transition of liquid water to single-walled INTs with diameters reaching ∼10 nm when confined to capillaries of double-walled carbon nanotubes (DW-CNTs). Three distinct classes of INTs are observed, namely, INTs with flat square walls (INTs-FSW), INTs with puckered rhombic walls (INTs-PRW), and INTs with bilayer hexagonal walls (INTs-BHW). Surprisingly, when water is confined in DW-CNT (3, 3)@(13, 13), an INT-FSW freezing temperature of 380 K can be reached, which is even higher than the boiling temperature of bulk water at atmospheric pressure. The freezing temperatures of INTs-FSW decrease as their caliber increases, approaching to the freezing temperature of two-dimensional flat square ice at the large-diameter limit. In contrast, the freezing temperature of INTs-PRW is insensitive to their diameter. Ab initio molecular dynamics simulations are performed to examine the stability of the INT-FSW and INT-PRW. The highly stable INTs with diameters beyond subnanometer scale can be exploited for potential applications in nanofluidic technologies and for mass transport as bioinspired nanochannels.

Widespread detection of chlorine oxyacids in the Arctic atmosphere

Chlorine radicals are strong atmospheric oxidants known to play an important role in the depletion of surface ozone and the degradation of methane in the Arctic troposphere. Initial oxidation processes of chlorine produce chlorine oxides, and it has been speculated that the final oxidation steps lead to the formation of chloric (HClO₃) and perchloric (HClO₄) acids, although these two species have not been detected in the atmosphere. Here, we present atmospheric observations of gas-phase HClO₃ and HClO₄. Significant levels of HClO₃ were observed during springtime at Greenland (Villum Research Station), Ny-Ålesund research station and over the central Arctic Ocean, on-board research vessel Polarstern during the Multidisciplinary drifting Observatory for the Study of the Arctic Climate (MOSAiC) campaign, with estimated concentrations up to 7 × 10⁶ molecule cm^-3. The increase in HClO₃, concomitantly with that in HClO₄, was linked to the increase in bromine levels. These observations indicated that bromine chemistry enhances the formation of OClO, which is subsequently oxidized into HClO₃ and HClO₄ by hydroxyl radicals. HClO₃ and HClO₄ are not photoactive and therefore their loss through heterogeneous uptake on aerosol and snow surfaces can function as a previously missing atmospheric sink for reactive chlorine, thereby reducing the chlorine-driven oxidation capacity in the Arctic boundary layer. Our study reveals additional chlorine species in the atmosphere, providing further insights into atmospheric chlorine cycling in the polar environment.

Revealing Structural Insights of Solid Electrolyte Interphase in High-Concentrated Non-Flammable Electrolyte for Li Metal Batteries by Cryo-TEM

High-concentrated non-flammable electrolytes (HCNFE) in lithium metal batteries prevent thermal runaway accidents, but the microstructure of their solid electrolyte interphase (SEI) remains largely unexplored, due to the lack of direct imaging tools. Herein, cryo-HRTEM is applied to directly visualize the native state of SEI at the atomic scale. In HCNFE, SEI has a uniform laminated crystalline-amorphous structure that can prevent further reaction between the electrolyte and lithium. The inorganic SEI component, Li₂ S₂ O₇ , is precisely identified by cryo-HRTEM. Density functional theory (DFT) calculations demonstrate that the final Li₂ S₂ O₇ phase has suitable natural transmission channels for Li-ion diffusion and excellent ionic conductivity of 1.2 × 10^-5 S cm^-1 .

Fast Hydroxyl Radical Generation at the Air-Water Interface of Aerosols Mediated by Water-Soluble PM2.5 under Ultraviolet A Radiation

Due to the adverse health effects and the role in the formation of secondary organic aerosols, hydroxyl radical (OH) generation by atmospheric fine particulate matter (PM) has been of particular research interest in both bulk solutions and the gas phase. However, OH generation by PM at the air-water interface of atmospheric water droplets, a unique environment where reactions can be accelerated by orders of magnitude, has long been overlooked. Using the field-induced droplet ionization mass spectrometry methodology that selectively samples molecules at the air-water interface, here, we show significant oxidation of amphiphilic lipids and isoprene mediated by water-soluble PM_2.5 at the air-water interface under ultraviolet A irradiation, with the OH generation rate estimated to be 1.5 × 10¹⁶ molecule·s^-1·m^-2. Atomistic molecular dynamics simulations support the counter-intuitive affinity for the air-water interface of isoprene. We opine that it is the carboxylic chelators of the surface-active molecules in PM that enrich photocatalytic metals such as iron at the air-water interface and greatly enhance the OH generation therein. This work provides a potential new heterogeneous OH generation channel in the atmosphere.

Charge Density Evolution Governing Interfacial Friction

It is well-known that the electron nature of a solid in contact plays a predominant role in determining the many properties of the contact systems, but the general rules of electron coupling that govern interfacial friction remain an open issue for the surface/interface community. Here, density functional theory calculations were used to investigate the physical origins of friction of solid interfaces. It was found that interfacial friction can be inherently traced back to the electronic barrier to the change in the contact configuration of the joints in slip due to the resistance of energy level rearrangement leading to electron transfer, which applies for various interface types ranging from van der Waals, metallic, and ionic to covalent joints. The variation of the electron density accompanying contact conformation changes along the sliding pathways is defined to track the frictional energy dissipation process occurring in slip. The results demonstrate that the frictional energy landscapes evolve synchronously with responding charge density evolution along sliding pathways, yielding an explicitly linear dependence of frictional dissipation on electronic evolution. The correlation coefficient enables us to interpret the fundamental concept of shear strength. The present charge evolution model thereby provides insights into the classic hypothesis that the friction force scales with the real contact area. This may shed light on the intrinsic origin of friction at the electronic level, opening the way to the rational design of nanomechanical devices as well as the understanding of the natural faults.

Counterintuitive Oxidation of Alcohols at Air-Water Interfaces

This study shows that the oxidation of alcohols can rapidly occur at air-water interfaces. It was found that methanediols (HOCH₂OH) orient at air-water interfaces with a H atom of the -CH₂- group pointing toward the gaseous phase. Counterintuitively, gaseous hydroxyl radicals do not prefer to attack the exposed -CH₂- group but the -OH group that forms hydrogen bonds with water molecules at the surface via a water-promoted mechanism, leading to the formation of formic acids. Compared with gaseous oxidation, the water-promoted mechanism at the air-water interface significantly lowers free-energy barriers from ∼10.7 to ∼4.3 kcal·mol^-1 and therefore accelerates the formation of formic acids. The study unveils a previously overlooked source of environmental organic acids that are bound up with aerosol formation and water acidity.

Water Complex of Imidogen

Imidogen (NH) is the simplest nitrogen hydride that plays an important role in combustion and interstellar chemistry, and its combination with H₂O is the prototypical amidation reaction of O-H bonds involving a nitrene intermediate. Herein, we report the observation of the elusive water complex of NH, a prereaction complex associated with the amidation reaction in a solid N₂ matrix at 10 K. The hydrogen-bonded structure of NH···OH₂ (versus HN···HOH) is confirmed via IR spectroscopy with comprehensive isotope labeling (D, ¹⁸O, and ¹⁵N) and quantum chemical calculations at the UCCSD(T)/aug-cc-pVQZ level of theory. In line with the observed absorption at 350 nm, irradiation of the complex at 365 nm leads to O-H bond insertion, yielding hydroxylamine NH₂OH.

Mechanistic Insights into the Reactive Uptake of Chlorine Nitrate at the Air-Water Interface

It is well-known that the aqueous-phase processing of chlorine nitrate (ClONO₂) plays a crucial role in ozone depletion. However, many of the physical and chemical properties of ClONO₂ at the air-water interface or in bulk water are unknown or not understood on a microscopic scale. Here, the solvation and hydrolysis of ClONO₂ at the air-water interface and in bulk water at 300 K were investigated by classical and ab initio molecular dynamics (AIMD) simulations combined with free energy methods. Our results revealed that ClONO₂ prefers to accumulate at the air-water interface rather than in the bulk phase. Specifically, halogen bonding interactions (ClONO₂)Cl···O(H₂O) were found to be the predominant interactions between ClONO₂ and H₂O. Moreover, metadynamics-biased AIMD simulations revealed that ClONO₂ hydrolysis is catalyzed at the air-water interface with an activation barrier of only ∼0.2 kcal/mol; additionally, the difference in free energy between the product and reactant is only ∼0.1 kcal/mol. Surprisingly, the near-barrierless reaction and the comparable free energies of the reactant and product suggested that the ClONO₂ hydrolysis at the air-water interface is reversible. When the temperature is lowered from 300 to 200 K, the activation barrier for the ClONO₂ hydrolysis at the air-water interface is increased to ∼5.4 kcal/mol. These findings have important implications for the interpretation of experiments.

Spectroscopy and photochemistry of ClSSO

Sulfur-chlorine cycles play a role in the atmosphere of Venus. It is thought that many sulfur-chlorine bearing molecules could be present in Venus's atmosphere and play an important role in its chemical processes. The goal of this work is to provide new insight into the electronic structure and spectroscopy of the [Cl, S, S, O] molecular system. Eight isomers could be formed, but only three were found to be thermodynamically stable relative to the first dissociation limit. We spectroscopically characterized the two lowest energy stable isomers, C1-ClSSO and trans-ClSSO, using the accurate CCSD(T)-F12/aug-cc-pVTZ method. The dipole moments of the two lowest energy stable isomers are predicted to be 1.90 and 1.60 debye, respectively. The C1-ClSSO isomer is suitable for laser induced fluorescence detection since the lowest excited electronic states absorb in the visible, ∼610 nm, and near UV region, 330 nm. We mapped the evolution of the low-lying excited electronic states along the ClS, SS, and SO bond lengths to find that the production of ClS, SO, or S₂O is plausible, whereas the production of ClS₂ is not allowed.

Efficient exploration of complex free-energy landscapes by stepwise multi-subphase space metadynamics

We present an efficient method based on an extension of metadynamics for exploring complex free energy landscapes (FELs). The method employs two-step metadynamics simulations. In the first step, rapid metadynamics simulations using broad and tall Gaussians are performed to identify a free energy pathway (FEP) connecting two states of interest. The FEP is then divided into a series of independent subphase spaces that comprise selected discrete images of the system. Using appropriate collective variables (CVs) chosen according to the FEP, the accurate FEL of each subphase space is separately calculated in subsequent divide-and-conquer metadynamics simulations with narrow and low Gaussians. Finally, all FELs calculated in each subphase phase are merged to obtain the full FEL. We show that the method greatly improves the performance of the metadynamics approach. In particular, we were able to efficiently model chemical systems with complex FELs, such as chemical reactions at the air/water interface. We demonstrated the performance of this method on two model reactions: the hydrolysis of formaldehyde in the gas phase and at the air/water interface.

Evidence of Formation of Monolayer Hydrated Salts in Nanopores

Despite much effort being devoted to the study of ionic aqueous solutions at the nanoscale, our fundamental understanding of the microscopic kinetic and thermodynamic behaviors in these systems remains largely incomplete. Herein, we reported the first 10 μs molecular dynamics simulation, providing evidence of the spontaneous formation of monolayer hexagonal honeycomb hydrated salts of XCl₂·6H₂O (X = Ba, Sr, Ca, and Mg) from electrolyte aqueous solutions confined in an angstrom-scale slit under ambient conditions. By using both the classical molecular dynamics simulations and the first-principles Born-Oppenheimer molecular dynamics simulations, we further demonstrated that the hydrated salts were stable not only at ambient temperature but also at elevated temperatures. This phenomenon of formation of hydrated salt in water is contrary to the conventional view. The free energy calculations and dehydration analyses indicated that the spontaneous formation of hydrated salts can be attributed to the interplay between ion hydration and Coulombic attractions in the highly confined water. In addition to providing molecular-level insights into the novel behavior of ionic aqueous solutions at the nanoscale, our findings may have implications for the future exploration of potential existence of water molecules in the saline deposits on hot planets.