Pub Date : 2026-01-06DOI: 10.1021/acsearthspacechem.5c00248
Matteo Michielan, , , Kim Steenbakkers, , , Daniela Ascenzi, , , Jake A. Diprose, , , Miroslav Polášek, , , Sandra Brünken, , , Claire Romanzin, , , Brianna R. Heazlewood, , , Cristina Puzzarini, , and , Vincent Richardson*,
The thioketenyl cation (HCCS+) has been recently detected in the dark cloud TMC-1 by radioastronomical observations within the QUIJOTE survey. However, the infrared (IR) spectrum of this ion is yet to be reported in the literature. Spectroscopic reference data are essential for the search of HCCS+ using the James Webb Space Telescope, not only in molecular clouds and star-forming regions, but also in the ionospheres and upper atmospheres of exoplanets. In this work, we demonstrate a method for the selective generation of the HCCS+ ion in its triplet ground state (3Σ–) and use this method to obtain IR band positions for HCCS+. The IR-action spectrum of H2-tagged HCCS+ has been measured in a cryogenic 22-pole ion trap via IR photodissociation (IR-PD) spectroscopy with the FELIX light source in the wavenumber regions 450–1850 and 3000–3350 cm–1. Spectral information is complemented by theoretical calculations on the fragmentation mechanisms leading to the formation of HCCS+ from dissociative ionization of 2,5-dibromothiophene. The assignment of the experimental HCCS+ vibrational bands is aided by comparison with ab initio computed values from literature and from calculations at the UB3LYP/cc-pVQZ level of theory, for both the triplet (3Σ–) and singlet (1Σ+) states of HCCS+. The experimental HCCS+ spectra show an overall good agreement with the scaled theoretical values (to account for anharmonicity effects), facilitating assignment of the IR spectral features. These findings will enable new reactivity investigations and spectroscopic measurements to be conducted, and for HCCS+ to be included in astrochemical models and databases.
{"title":"IR-Action Spectroscopy of the Astrochemically Relevant HCCS+ Cation","authors":"Matteo Michielan, , , Kim Steenbakkers, , , Daniela Ascenzi, , , Jake A. Diprose, , , Miroslav Polášek, , , Sandra Brünken, , , Claire Romanzin, , , Brianna R. Heazlewood, , , Cristina Puzzarini, , and , Vincent Richardson*, ","doi":"10.1021/acsearthspacechem.5c00248","DOIUrl":"https://doi.org/10.1021/acsearthspacechem.5c00248","url":null,"abstract":"<p >The thioketenyl cation (HCCS<sup>+</sup>) has been recently detected in the dark cloud TMC-1 by radioastronomical observations within the QUIJOTE survey. However, the infrared (IR) spectrum of this ion is yet to be reported in the literature. Spectroscopic reference data are essential for the search of HCCS<sup>+</sup> using the James Webb Space Telescope, not only in molecular clouds and star-forming regions, but also in the ionospheres and upper atmospheres of exoplanets. In this work, we demonstrate a method for the selective generation of the HCCS<sup>+</sup> ion in its triplet ground state (<sup>3</sup>Σ<sup>–</sup>) and use this method to obtain IR band positions for HCCS<sup>+</sup>. The IR-action spectrum of H<sub>2</sub>-tagged HCCS<sup>+</sup> has been measured in a cryogenic 22-pole ion trap via IR photodissociation (IR-PD) spectroscopy with the FELIX light source in the wavenumber regions 450–1850 and 3000–3350 cm<sup>–1</sup>. Spectral information is complemented by theoretical calculations on the fragmentation mechanisms leading to the formation of HCCS<sup>+</sup> from dissociative ionization of 2,5-dibromothiophene. The assignment of the experimental HCCS<sup>+</sup> vibrational bands is aided by comparison with <i>ab initio</i> computed values from literature and from calculations at the UB3LYP/cc-pVQZ level of theory, for both the triplet (<sup>3</sup>Σ<sup>–</sup>) and singlet (<sup>1</sup>Σ<sup>+</sup>) states of HCCS<sup>+</sup>. The experimental HCCS<sup>+</sup> spectra show an overall good agreement with the scaled theoretical values (to account for anharmonicity effects), facilitating assignment of the IR spectral features. These findings will enable new reactivity investigations and spectroscopic measurements to be conducted, and for HCCS<sup>+</sup> to be included in astrochemical models and databases.</p>","PeriodicalId":15,"journal":{"name":"ACS Earth and Space Chemistry","volume":"10 1","pages":"148–156"},"PeriodicalIF":2.9,"publicationDate":"2026-01-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://pubs.acs.org/doi/pdf/10.1021/acsearthspacechem.5c00248","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145972669","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-02DOI: 10.1021/acsearthspacechem.5c00243
Wallace C. H. Hui, and , Kono H. Lemke*,
<p >Using <i>ab initio</i> molecular dynamics and thermodynamic integration calculations, we have examined the speciation and stability of H<sub>2</sub>SO<sub>3</sub> from 25 to 350 °C in bulk liquid and water nanodroplets. H<sub>2</sub>SO<sub>3</sub> ionization constants have been determined by using constraint force and thermodynamic integration techniques. The temperature dependence of the first ionization (yielding HSO<sub>3</sub><sup>–</sup>) and second ionization (yielding SO<sub>3</sub><sup>2–</sup>) constants up to 350 °C is given by the following van’t Hoff equations: log<i>K</i><sub>1</sub> = −32.5686 – 0.0297·<i>T</i> + 6.8894·ln(<i>T</i>) and log<i>K</i><sub>2</sub> = −86.3360 – 0.0608·<i>T</i> + 16.8911·ln(<i>T</i>). Our calculated ionization constants are in good agreement with experimental data at 25 °C, with p<i>K</i><sub>1</sub> = 2.08 ± 0.13 (Exp. = 1.88, [<contrib-group><span>Beyad, Y.</span></contrib-group> <cite><i>Dalton Trans.</i></cite> <span>2014</span>, <em>43</em>(5), 2147–2152]) and p<i>K</i><sub>2</sub> = 8.14 ± 0.10 (Exp. = 7.12, [<contrib-group><span>Millero, F. J.</span></contrib-group> <cite><i>J. Atmos. Chem.</i></cite> <span>1989</span>, <em>8</em>, 377–389]). With increasing temperature, values of p<i>K</i><sub>1</sub> and p<i>K</i><sub>2</sub> increase to 6.93 ± 0.13 and 15.66 ± 0.16 at 350 °C, respectively, highlighting the reduced acidity of H<sub>2</sub>SO<sub>3</sub> and HSO<sub>3</sub><sup>–</sup>. The derived enthalpy changes (Δ<i>H</i>°) shift toward more exothermic values with increasing temperature from −2.67 kcal/mol at 25 °C to −33.08 kcal/mol at 350 °C, for the first ionization step, and from −1.69 kcal/mol at 25 °C to −59.90 kcal/mol at 350 °C kcal/mol for the second ionization step, the temperature-weighted entropy (<i>T</i>Δ<i>S</i>°) shift from −5.50 kcal/mol (25 °C) to −52.85 kcal/mol (350 °C) and from −12.80 kcal/mol (25 °C) to −104.56 kcal/mol (350 °C) for the first and second ionization steps, respectively. Calculated values of Δ<i>H</i>° and Δ<i>S</i>° demonstrate that with increasing temperature entropic effects outweigh enthalpic contributions, giving rise to a decreased acidity of H<sub>2</sub>SO<sub>3</sub>. Vapor-phase simulations of H<sub>2</sub>SO<sub>3</sub> at 350 °C reveal a ∼30% larger ionization free energy Δ<i>G</i> compared to solution phase ionization, indicating an enhanced energetic penalty for deprotonation in low-density aqueous media. Simulations in H<sub>2</sub>SO<sub>3</sub>(H<sub>2</sub>O)<sub>100</sub> nanodroplets at 25 and 100 °C demonstrate surface-enhanced stabilization for H<sub>2</sub>SO<sub>3</sub>, with the deprotonation free energy Δ<i>G</i> 2.6–4.6 kcal/mol larger at surface solvation sites versus droplet interior. Finally, these molecular dynamics results reveal that at neutral pH in bulk liquid, sulfurous acid exists as the ionized HSO<sub>3</sub><sup>–</sup> form, while aqueous hydrothermal nanodroplet surfaces and interiors favor the fully associat
{"title":"Ab Initio Molecular Dynamics Study of H2SO3 Ionization in Hydrothermal Fluids and Nanodroplets","authors":"Wallace C. H. Hui, and , Kono H. Lemke*, ","doi":"10.1021/acsearthspacechem.5c00243","DOIUrl":"https://doi.org/10.1021/acsearthspacechem.5c00243","url":null,"abstract":"<p >Using <i>ab initio</i> molecular dynamics and thermodynamic integration calculations, we have examined the speciation and stability of H<sub>2</sub>SO<sub>3</sub> from 25 to 350 °C in bulk liquid and water nanodroplets. H<sub>2</sub>SO<sub>3</sub> ionization constants have been determined by using constraint force and thermodynamic integration techniques. The temperature dependence of the first ionization (yielding HSO<sub>3</sub><sup>–</sup>) and second ionization (yielding SO<sub>3</sub><sup>2–</sup>) constants up to 350 °C is given by the following van’t Hoff equations: log<i>K</i><sub>1</sub> = −32.5686 – 0.0297·<i>T</i> + 6.8894·ln(<i>T</i>) and log<i>K</i><sub>2</sub> = −86.3360 – 0.0608·<i>T</i> + 16.8911·ln(<i>T</i>). Our calculated ionization constants are in good agreement with experimental data at 25 °C, with p<i>K</i><sub>1</sub> = 2.08 ± 0.13 (Exp. = 1.88, [<contrib-group><span>Beyad, Y.</span></contrib-group> <cite><i>Dalton Trans.</i></cite> <span>2014</span>, <em>43</em>(5), 2147–2152]) and p<i>K</i><sub>2</sub> = 8.14 ± 0.10 (Exp. = 7.12, [<contrib-group><span>Millero, F. J.</span></contrib-group> <cite><i>J. Atmos. Chem.</i></cite> <span>1989</span>, <em>8</em>, 377–389]). With increasing temperature, values of p<i>K</i><sub>1</sub> and p<i>K</i><sub>2</sub> increase to 6.93 ± 0.13 and 15.66 ± 0.16 at 350 °C, respectively, highlighting the reduced acidity of H<sub>2</sub>SO<sub>3</sub> and HSO<sub>3</sub><sup>–</sup>. The derived enthalpy changes (Δ<i>H</i>°) shift toward more exothermic values with increasing temperature from −2.67 kcal/mol at 25 °C to −33.08 kcal/mol at 350 °C, for the first ionization step, and from −1.69 kcal/mol at 25 °C to −59.90 kcal/mol at 350 °C kcal/mol for the second ionization step, the temperature-weighted entropy (<i>T</i>Δ<i>S</i>°) shift from −5.50 kcal/mol (25 °C) to −52.85 kcal/mol (350 °C) and from −12.80 kcal/mol (25 °C) to −104.56 kcal/mol (350 °C) for the first and second ionization steps, respectively. Calculated values of Δ<i>H</i>° and Δ<i>S</i>° demonstrate that with increasing temperature entropic effects outweigh enthalpic contributions, giving rise to a decreased acidity of H<sub>2</sub>SO<sub>3</sub>. Vapor-phase simulations of H<sub>2</sub>SO<sub>3</sub> at 350 °C reveal a ∼30% larger ionization free energy Δ<i>G</i> compared to solution phase ionization, indicating an enhanced energetic penalty for deprotonation in low-density aqueous media. Simulations in H<sub>2</sub>SO<sub>3</sub>(H<sub>2</sub>O)<sub>100</sub> nanodroplets at 25 and 100 °C demonstrate surface-enhanced stabilization for H<sub>2</sub>SO<sub>3</sub>, with the deprotonation free energy Δ<i>G</i> 2.6–4.6 kcal/mol larger at surface solvation sites versus droplet interior. Finally, these molecular dynamics results reveal that at neutral pH in bulk liquid, sulfurous acid exists as the ionized HSO<sub>3</sub><sup>–</sup> form, while aqueous hydrothermal nanodroplet surfaces and interiors favor the fully associat","PeriodicalId":15,"journal":{"name":"ACS Earth and Space Chemistry","volume":"10 1","pages":"120–130"},"PeriodicalIF":2.9,"publicationDate":"2026-01-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145969072","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-12-29DOI: 10.1021/acsearthspacechem.5c00311
Francesco Benedetti, , , Mauro Satta, , , Tommaso Grassi, , , Stefan Vogt-Geisse, , and , Stefano Bovino*,
Surface chemistry on interstellar dust grains is recognized as a central component in astrochemical models, representing a plausible formation route for many of the observed complex molecular species. However, key parameters governing interstellar surface chemistry, such as diffusion energy barriers, remain poorly constrained. In particular, surface diffusion constitutes a fundamental step in the synthesis of complex organic molecules and plays a crucial role in understanding the desorption process. In this paper, the diffusion dynamics of carbon monoxide (CO) on amorphous solid water (ASW) surfaces, representative of interstellar ices, is modeled with quantum-chemical methods. Employing a representative ensemble of water clusters, each made of 22 molecules, diffusion energy barriers between the binding sites are computed using density functional theory. Diffusion rate coefficients are then determined by applying the harmonic approximation of the transition state theory. The results, in agreement with experimental studies, revealed a wide distribution of diffusion energies. This reflects the intrinsic topological heterogeneity of ASW surfaces and highlights how surface mobility significantly influences CO’s desorption dynamics and, as a consequence, surface-mediated reactivity in interstellar environments. We show that key parameters commonly employed in astrochemical models, such as the ratio between binding and diffusion energy, should be carefully revised.
{"title":"CO Diffusion on Interstellar Amorphous Solid Water: A Computational Study","authors":"Francesco Benedetti, , , Mauro Satta, , , Tommaso Grassi, , , Stefan Vogt-Geisse, , and , Stefano Bovino*, ","doi":"10.1021/acsearthspacechem.5c00311","DOIUrl":"https://doi.org/10.1021/acsearthspacechem.5c00311","url":null,"abstract":"<p >Surface chemistry on interstellar dust grains is recognized as a central component in astrochemical models, representing a plausible formation route for many of the observed complex molecular species. However, key parameters governing interstellar surface chemistry, such as diffusion energy barriers, remain poorly constrained. In particular, surface diffusion constitutes a fundamental step in the synthesis of complex organic molecules and plays a crucial role in understanding the desorption process. In this paper, the diffusion dynamics of carbon monoxide (CO) on amorphous solid water (ASW) surfaces, representative of interstellar ices, is modeled with quantum-chemical methods. Employing a representative ensemble of water clusters, each made of 22 molecules, diffusion energy barriers between the binding sites are computed using density functional theory. Diffusion rate coefficients are then determined by applying the harmonic approximation of the transition state theory. The results, in agreement with experimental studies, revealed a wide distribution of diffusion energies. This reflects the intrinsic topological heterogeneity of ASW surfaces and highlights how surface mobility significantly influences CO’s desorption dynamics and, as a consequence, surface-mediated reactivity in interstellar environments. We show that key parameters commonly employed in astrochemical models, such as the ratio between binding and diffusion energy, should be carefully revised.</p>","PeriodicalId":15,"journal":{"name":"ACS Earth and Space Chemistry","volume":"10 1","pages":"224–234"},"PeriodicalIF":2.9,"publicationDate":"2025-12-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://pubs.acs.org/doi/pdf/10.1021/acsearthspacechem.5c00311","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145969057","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Fluorinated esters (FESs) are emerging atmospheric species formed as secondary products from the degradation of hydrofluoroethers (HFEs) and other fluorinated precursors. Fluoroalkylfluoroformates (FAFs) are a subgroup of fluorinated esters, where a fluoroalkyl group is bonded to a fluoroformate group (C(O)F). In this work, we have done a thorough kinetic and mechanistic analysis of the atmospheric degradation of specific FAFs, including FC(O)OCHF2, FC(O)OCH2F, FC(O)OCH3, and FC(O)OCH2CF3, with hydroxyl (OH) radical and Chlorine (Cl) atoms. Geometry optimizations and frequency analyses were performed at the M06–2X/aug-cc-pVTZ and MP2/aug-cc-pVTZ levels. The single-point (SP) energies were refined using the CCSD(T)/aug-cc-pVXZ (D and T) basis sets, and complete basis set (CBS) extrapolation was performed to obtain more reliable rate coefficient values for the reactions. Eight distinct H-abstraction pathways were identified, and the corresponding potential energy surfaces were constructed. Rate constants were computed using canonical variational transition state theory with small-curvature tunneling (CVT/SCT) over the 250–500 K range. We also examined the impact of these substances on the atmosphere.
{"title":"A Theoretical Investigation on the Hydrogen Abstraction Reactions of FC(O)OX (X = CHF2, CH2F, CH3, CH2CF3) with OH Radicals and Cl Atoms and Their Atmospheric Impacts","authors":"Suresh Tiwari, , , Arvind Kumar, , and , Ranga Subramanian*, ","doi":"10.1021/acsearthspacechem.5c00268","DOIUrl":"https://doi.org/10.1021/acsearthspacechem.5c00268","url":null,"abstract":"<p >Fluorinated esters (FESs) are emerging atmospheric species formed as secondary products from the degradation of hydrofluoroethers (HFEs) and other fluorinated precursors. Fluoroalkylfluoroformates (FAFs) are a subgroup of fluorinated esters, where a fluoroalkyl group is bonded to a fluoroformate group (C(O)F). In this work, we have done a thorough kinetic and mechanistic analysis of the atmospheric degradation of specific FAFs, including FC(O)OCHF<sub>2</sub>, FC(O)OCH<sub>2</sub>F, FC(O)OCH<sub>3</sub>, and FC(O)OCH<sub>2</sub>CF<sub>3</sub>, with hydroxyl (OH) radical and Chlorine (Cl) atoms. Geometry optimizations and frequency analyses were performed at the M06–2X/aug-cc-pVTZ and MP2/aug-cc-pVTZ levels. The single-point (SP) energies were refined using the CCSD(T)/aug-cc-pVXZ (D and T) basis sets, and complete basis set (CBS) extrapolation was performed to obtain more reliable rate coefficient values for the reactions. Eight distinct H-abstraction pathways were identified, and the corresponding potential energy surfaces were constructed. Rate constants were computed using canonical variational transition state theory with small-curvature tunneling (CVT/SCT) over the 250–500 K range. We also examined the impact of these substances on the atmosphere.</p>","PeriodicalId":15,"journal":{"name":"ACS Earth and Space Chemistry","volume":"10 1","pages":"167–178"},"PeriodicalIF":2.9,"publicationDate":"2025-12-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145969060","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-12-24DOI: 10.1021/acsearthspacechem.5c00249
Ryan C. Fortenberry*, and , Vincent J. Esposito,
The long theorized polycyclic aromatic hydrocarbons (PAHs) crashed onto the astrochemical scene with the radioastronomical discovery of cyanobenzene (benzonitrile) in 2018, but the astronomical observation of the CN stretch for this or any other known CN-PAHs in the IR has yet to be reported. With the wealth of James Webb Space Telescope (JWST) data being actively returned, the need for reference data for IR features of known interstellar molecules like CN-PAHs has never been greater. Due to complexities in working with PAHs in the laboratory, quantum chemistry provides the best means of high-throughput vibrational reference data generation. This work gives an overview for the state-of-the-art in quantum chemical approaches and results for computing IR spectroscopic data for CN-PAHs, specifically for the CN stretch. The use of the hybrid rDSD/junTZ + B3LYP/N07D method through quartic force fields (QFFs) and second-order vibrational perturbation theory (VPT2) has computed the CN stretch of C2v 9-cyanoanthracene to be 2207 cm–1 (4.531 μm), exactly the same as free electron laser experiments report. Additionally, due to the structural regularity of the attached PAHs, the CN stretch does not vary greatly (less than 20 cm–1) from CN-PAH molecule to different CN-PAH molecule. However, this spectroscopic region is home to Mg IV and Ar VI lines that may hinder observations of the CN stretch in astronomical environments where vibrational emission is possible.
{"title":"Quantum Chemical Computation of Vibrational Spectra for CN-PAHs","authors":"Ryan C. Fortenberry*, and , Vincent J. Esposito, ","doi":"10.1021/acsearthspacechem.5c00249","DOIUrl":"https://doi.org/10.1021/acsearthspacechem.5c00249","url":null,"abstract":"<p >The long theorized polycyclic aromatic hydrocarbons (PAHs) crashed onto the astrochemical scene with the radioastronomical discovery of cyanobenzene (benzonitrile) in 2018, but the astronomical observation of the CN stretch for this or any other known CN-PAHs in the IR has yet to be reported. With the wealth of James Webb Space Telescope (JWST) data being actively returned, the need for reference data for IR features of known interstellar molecules like CN-PAHs has never been greater. Due to complexities in working with PAHs in the laboratory, quantum chemistry provides the best means of high-throughput vibrational reference data generation. This work gives an overview for the state-of-the-art in quantum chemical approaches and results for computing IR spectroscopic data for CN-PAHs, specifically for the CN stretch. The use of the hybrid rDSD/junTZ + B3LYP/N07D method through quartic force fields (QFFs) and second-order vibrational perturbation theory (VPT2) has computed the CN stretch of <i>C</i><sub>2<i>v</i></sub> 9-cyanoanthracene to be 2207 cm<sup>–1</sup> (4.531 μm), exactly the same as free electron laser experiments report. Additionally, due to the structural regularity of the attached PAHs, the CN stretch does not vary greatly (less than 20 cm<sup>–1</sup>) from CN-PAH molecule to different CN-PAH molecule. However, this spectroscopic region is home to Mg IV and Ar VI lines that may hinder observations of the CN stretch in astronomical environments where vibrational emission is possible.</p>","PeriodicalId":15,"journal":{"name":"ACS Earth and Space Chemistry","volume":"10 1","pages":"1–10"},"PeriodicalIF":2.9,"publicationDate":"2025-12-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145969102","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-12-24DOI: 10.1021/acsearthspacechem.5c00214
Christina Breitenstein, , , Manuela van Pinxteren, , , Kenneth Wagner, , , Lorenz Reschke, , , Sebastian Zeppenfeld, , and , Hartmut Herrmann*,
A sensitive analytical method with a run time of 17 min is presented to detect and quantify free and combined amino acids (FAA/CAA) and nucleobases in complex marine matrices. The method combines hydrophilic interaction chromatography (HILIC) with quadrupole time-of-flight mass spectrometry (Q-TOF-MS). Retention time interday relative standard deviation (RSD) was 0.09–7.1%; peak area RSD was <6%; and individual limits of quantification (LOQs) ranged between 0.5 and 5 μg/L. The optimized solid-phase extraction (SPE) desalting protocol for salt water achieved a recovery of more than 50%, except for threonine (30%) and uracil (15%). Despite lower recoveries, compounds were consistently and reliably detected in ambient samples. Recovery RSD was ≤10% for most analytes, with uracil (23%) and cystine (13%) as exceptions. Hydrolysis in polypropylene vials to quantify CAA yielded recoveries between 61 and 102%. Matrix effects were evaluated via standard addition and isotope-labeled standards, revealing ion suppression/enhancement in desalted seawater samples, thus confirming the need for matrix-matched calibration. For aerosol particle extracts, no desalting procedure was necessary. The method was applied to polar-region samples, yielding atmospheric concentrations that ranged between 0.5 pg/m3 for adenine and 3.5 ng/m3 for glutamine. In seawater samples, concentrations ranged between 1.3 μg/L for glycine and 48 μg/L for uracil. Altogether, the optimized method enables reliable analysis, without the need for a derivatization step. Thus, it is allowing for the sensitive determination of amino acids and nucleobases that are challenging to analyze and rarely reported in environmental studies to date.
{"title":"Protocol for the Analysis of Combined and Free Amino Acids in Seawater and Marine Aerosol Particles Using Hydrophilic Interaction LC–TOF-MS","authors":"Christina Breitenstein, , , Manuela van Pinxteren, , , Kenneth Wagner, , , Lorenz Reschke, , , Sebastian Zeppenfeld, , and , Hartmut Herrmann*, ","doi":"10.1021/acsearthspacechem.5c00214","DOIUrl":"https://doi.org/10.1021/acsearthspacechem.5c00214","url":null,"abstract":"<p >A sensitive analytical method with a run time of 17 min is presented to detect and quantify free and combined amino acids (FAA/CAA) and nucleobases in complex marine matrices. The method combines hydrophilic interaction chromatography (HILIC) with quadrupole time-of-flight mass spectrometry (Q-TOF-MS). Retention time interday relative standard deviation (RSD) was 0.09–7.1%; peak area RSD was <6%; and individual limits of quantification (LOQs) ranged between 0.5 and 5 μg/L. The optimized solid-phase extraction (SPE) desalting protocol for salt water achieved a recovery of more than 50%, except for threonine (30%) and uracil (15%). Despite lower recoveries, compounds were consistently and reliably detected in ambient samples. Recovery RSD was ≤10% for most analytes, with uracil (23%) and cystine (13%) as exceptions. Hydrolysis in polypropylene vials to quantify CAA yielded recoveries between 61 and 102%. Matrix effects were evaluated via standard addition and isotope-labeled standards, revealing ion suppression/enhancement in desalted seawater samples, thus confirming the need for matrix-matched calibration. For aerosol particle extracts, no desalting procedure was necessary. The method was applied to polar-region samples, yielding atmospheric concentrations that ranged between 0.5 pg/m<sup>3</sup> for adenine and 3.5 ng/m<sup>3</sup> for glutamine. In seawater samples, concentrations ranged between 1.3 μg/L for glycine and 48 μg/L for uracil. Altogether, the optimized method enables reliable analysis, without the need for a derivatization step. Thus, it is allowing for the sensitive determination of amino acids and nucleobases that are challenging to analyze and rarely reported in environmental studies to date.</p>","PeriodicalId":15,"journal":{"name":"ACS Earth and Space Chemistry","volume":"10 1","pages":"80–93"},"PeriodicalIF":2.9,"publicationDate":"2025-12-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://pubs.acs.org/doi/pdf/10.1021/acsearthspacechem.5c00214","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145969089","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-12-23DOI: 10.1021/acsearthspacechem.5c00306
Qihui Shi, , , Zhaodong Xi*, , , Yidan Ding, , , Gary G. Lash, , and , Shuheng Tang,
Shale organic matter (OM) enrichment plays a critical role in evaluating the future prospects of shale gas worldwide. The Ordovician-Silurian transition represents a critical interval notable for extensive deposits of organic-rich shale formations. Despite this, the underlying factors controlling OM enrichment during this interval, particularly in the Wufeng-Longmaxi (WF-LMX) shale, need comprehensive reevaluation. Our study compiles major and trace element concentrations for 1376 samples and TOC values for 3143 samples from the WF-LMX shale in South China to define quantitative thresholds for redox environment and paleo-productivity. Paleoenvironmental characteristics associated with the deposition of shales containing varying OM abundances (TOC <2%, 2–4%, and >4%) were clarified. Random Forest (RF) and Artificial Neural Network (ANN) models were applied to identify the primary factors controlling OM enrichment of the WF-LMX succession. Results of this work show that MoEF, UEF, VEF, and CrEF values of 10, 10, 1.5, and 1.5, respectively, represent the threshold for anoxic conditions during shale deposition. SiO2(exc), P/Ti, and Babio values of approximately 30%, 0.2, and 1000 ppm, respectively, signify high paleo-productivity. RF and ANN results suggest that redox environment and paleo-productivity as the foremost influences on OM accumulation, followed by paleoclimate. Organic-rich shales (Type-a, TOC > 4%) were primarily deposited under warm and humid climate conditions. Enhanced nutrient supply promoted elevated productivity and established anoxic bottom conditions that facilitated OM preservation. Type-c shales (TOC < 2%) accumulated under low-moderate productivity and oxygenated bottom conditions, leading to extensive degradation and dilution of organic matter. Type-b shales (TOC = 2–4%) exhibit intermediate characteristics between Type-a and Type-c. This study offers new perspectives on OM enrichment mechanisms in the WF-LMX shale succession and provides a theoretical basis for shale gas exploration within this stratigraphic interval.
{"title":"Insights of Big Data and Machine Learning on the Depositional Environment and Organic Matter Enrichment of the Wufeng-Longmaxi Shales","authors":"Qihui Shi, , , Zhaodong Xi*, , , Yidan Ding, , , Gary G. Lash, , and , Shuheng Tang, ","doi":"10.1021/acsearthspacechem.5c00306","DOIUrl":"https://doi.org/10.1021/acsearthspacechem.5c00306","url":null,"abstract":"<p >Shale organic matter (OM) enrichment plays a critical role in evaluating the future prospects of shale gas worldwide. The Ordovician-Silurian transition represents a critical interval notable for extensive deposits of organic-rich shale formations. Despite this, the underlying factors controlling OM enrichment during this interval, particularly in the Wufeng-Longmaxi (WF-LMX) shale, need comprehensive reevaluation. Our study compiles major and trace element concentrations for 1376 samples and TOC values for 3143 samples from the WF-LMX shale in South China to define quantitative thresholds for redox environment and paleo-productivity. Paleoenvironmental characteristics associated with the deposition of shales containing varying OM abundances (TOC <2%, 2–4%, and >4%) were clarified. Random Forest (RF) and Artificial Neural Network (ANN) models were applied to identify the primary factors controlling OM enrichment of the WF-LMX succession. Results of this work show that Mo<sub>EF</sub>, <i>U</i><sub>EF</sub>, <i>V</i><sub>EF</sub>, and Cr<sub>EF</sub> values of 10, 10, 1.5, and 1.5, respectively, represent the threshold for anoxic conditions during shale deposition. SiO<sub>2(exc)</sub>, P/Ti, and Ba<sub>bio</sub> values of approximately 30%, 0.2, and 1000 ppm, respectively, signify high paleo-productivity. RF and ANN results suggest that redox environment and paleo-productivity as the foremost influences on OM accumulation, followed by paleoclimate. Organic-rich shales (Type-a, TOC > 4%) were primarily deposited under warm and humid climate conditions. Enhanced nutrient supply promoted elevated productivity and established anoxic bottom conditions that facilitated OM preservation. Type-c shales (TOC < 2%) accumulated under low-moderate productivity and oxygenated bottom conditions, leading to extensive degradation and dilution of organic matter. Type-b shales (TOC = 2–4%) exhibit intermediate characteristics between Type-a and Type-c. This study offers new perspectives on OM enrichment mechanisms in the WF-LMX shale succession and provides a theoretical basis for shale gas exploration within this stratigraphic interval.</p>","PeriodicalId":15,"journal":{"name":"ACS Earth and Space Chemistry","volume":"10 1","pages":"210–223"},"PeriodicalIF":2.9,"publicationDate":"2025-12-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145972681","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-12-23DOI: 10.1021/acsearthspacechem.5c00223
Rafael Flores-Larrañaga, , , Luis Armando Gonzalez-Ortiz, , , Filiberto Ortíz-Chi, , , María Eugenia Castro, , , Francisco J. Melendez*, , , Lisset Noriega*, , and , Gabriel Merino*,
The sulfur depletion problem in the interstellar medium motivates a detailed characterization of sulfur-bearing molecules such as CH4S2, whose isomers may act as potential reservoirs of interstellar sulfur. Using high-level ab initio methods, we identified eight structural isomers of CH4S2, with 1-methylhydrodisulfide (CH3SSH) as the global minimum and methanedithiol lying 3.1 kcal·mol–1 higher in energy. These low-energy isomers show both kinetic and thermodynamic stability against decomposition into H2CS and H2S, even in the presence of water. In contrast, hypervalent sulfur isomers are more than 60 kcal·mol–1 higher in energy. Rotational constants obtained with the jun-Cheap Scheme (jun-ChS), including anharmonic corrections reproduce experimental values within 0.17%, ensuring reliable spectral predictions. The stability of CH4S2 isomers under interstellar conditions and their dipole moments (up to 1.9 D) suggest that they are promising targets for radioastronomical detection and may help clarify the chemical pathways of sulfur in cold molecular clouds.
{"title":"Computational Characterization of CH4S2 Isomers as Key Candidates in Interstellar Sulfur Chemistry","authors":"Rafael Flores-Larrañaga, , , Luis Armando Gonzalez-Ortiz, , , Filiberto Ortíz-Chi, , , María Eugenia Castro, , , Francisco J. Melendez*, , , Lisset Noriega*, , and , Gabriel Merino*, ","doi":"10.1021/acsearthspacechem.5c00223","DOIUrl":"https://doi.org/10.1021/acsearthspacechem.5c00223","url":null,"abstract":"<p >The sulfur depletion problem in the interstellar medium motivates a detailed characterization of sulfur-bearing molecules such as CH<sub>4</sub>S<sub>2</sub>, whose isomers may act as potential reservoirs of interstellar sulfur. Using high-level <i>ab initio</i> methods, we identified eight structural isomers of CH<sub>4</sub>S<sub>2</sub>, with 1-methylhydrodisulfide (CH<sub>3</sub>SSH) as the global minimum and methanedithiol lying 3.1 kcal·mol<sup>–1</sup> higher in energy. These low-energy isomers show both kinetic and thermodynamic stability against decomposition into H<sub>2</sub>CS and H<sub>2</sub>S, even in the presence of water. In contrast, hypervalent sulfur isomers are more than 60 kcal·mol<sup>–1</sup> higher in energy. Rotational constants obtained with the jun-Cheap Scheme (jun-ChS), including anharmonic corrections reproduce experimental values within 0.17%, ensuring reliable spectral predictions. The stability of CH<sub>4</sub>S<sub>2</sub> isomers under interstellar conditions and their dipole moments (up to 1.9 D) suggest that they are promising targets for radioastronomical detection and may help clarify the chemical pathways of sulfur in cold molecular clouds.</p>","PeriodicalId":15,"journal":{"name":"ACS Earth and Space Chemistry","volume":"10 1","pages":"110–119"},"PeriodicalIF":2.9,"publicationDate":"2025-12-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://pubs.acs.org/doi/pdf/10.1021/acsearthspacechem.5c00223","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145969085","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-12-22DOI: 10.1021/acsearthspacechem.5c00187
Silvia Spezzano*, , , Wiebke Riedel, , , Paola Caselli, , , Olli Sipilä, , , Yuxin Lin, , , Hayley A. Bunn, , , Elena Redaelli, , , Laurent H. Coudert, , , Andrés Megías, , and , Izaskun Jimenez-Serra,
Isotopic fractionation is a very powerful tool to follow the evolution of material from one stage to the next in the star-formation process. Prestellar cores exhibit some of the highest levels of deuteration because their physical conditions (T ≤ 10 K and n(H2) ≥ 105 cm–3) greatly favor deuteration processes. Deuteration maps are a measure of the effectiveness of the deuteration across the core, and they are useful to study both the deuteration and the formation mechanism (either in the gas-phase or on grain surfaces) of the main species. Methanol is the simplest complex organic molecule (COM) that is O-bearing and detected in the interstellar medium (ISM). It represents the beginning of molecular complexity in star-forming regions; thus, a complete understanding of its formation and deuteration is a necessary step to understand the development of further chemical complexity. In this paper, we use single-dish observations with the IRAM 30 m telescope and state-of-the-art chemical models to investigate the deuteration of methanol toward the prototypical prestellar core L1544. We also compare the results of the chemical models with previous observations of deuterated methanol toward the presttellar cores HMM1 and L694-2. The spectra extracted from the CHD2OH map show that the emission is concentrated in the center and toward the northwest of the core. Using deep observations toward the dust and the methanol peaks of the core, we derive a very large deuterium fraction for methanol (∼20%) toward both peaks. The comparison of our observational results with chemical models has highlighted the importance of H-abstraction processes in the formation and deuteration of methanol. Deep observations combined with state-of-the-art chemical models are of fundamental importance in understanding the development of molecular complexity in the ISM. Our analysis also shows the importance of non-LTE effects when measuring the D/H ratios in methanol.
{"title":"High Deuteration of Methanol in L1544","authors":"Silvia Spezzano*, , , Wiebke Riedel, , , Paola Caselli, , , Olli Sipilä, , , Yuxin Lin, , , Hayley A. Bunn, , , Elena Redaelli, , , Laurent H. Coudert, , , Andrés Megías, , and , Izaskun Jimenez-Serra, ","doi":"10.1021/acsearthspacechem.5c00187","DOIUrl":"https://doi.org/10.1021/acsearthspacechem.5c00187","url":null,"abstract":"<p >Isotopic fractionation is a very powerful tool to follow the evolution of material from one stage to the next in the star-formation process. Prestellar cores exhibit some of the highest levels of deuteration because their physical conditions (<i>T</i> ≤ 10 K and <i>n</i>(H<sub>2</sub>) ≥ 10<sup>5</sup> cm<sup>–3</sup>) greatly favor deuteration processes. Deuteration maps are a measure of the effectiveness of the deuteration across the core, and they are useful to study both the deuteration and the formation mechanism (either in the gas-phase or on grain surfaces) of the main species. Methanol is the simplest complex organic molecule (COM) that is O-bearing and detected in the interstellar medium (ISM). It represents the beginning of molecular complexity in star-forming regions; thus, a complete understanding of its formation and deuteration is a necessary step to understand the development of further chemical complexity. In this paper, we use single-dish observations with the IRAM 30 m telescope and state-of-the-art chemical models to investigate the deuteration of methanol toward the prototypical prestellar core L1544. We also compare the results of the chemical models with previous observations of deuterated methanol toward the presttellar cores HMM1 and L694-2. The spectra extracted from the CHD<sub>2</sub>OH map show that the emission is concentrated in the center and toward the northwest of the core. Using deep observations toward the dust and the methanol peaks of the core, we derive a very large deuterium fraction for methanol (∼20%) toward both peaks. The comparison of our observational results with chemical models has highlighted the importance of H-abstraction processes in the formation and deuteration of methanol. Deep observations combined with state-of-the-art chemical models are of fundamental importance in understanding the development of molecular complexity in the ISM. Our analysis also shows the importance of non-LTE effects when measuring the D/H ratios in methanol.</p>","PeriodicalId":15,"journal":{"name":"ACS Earth and Space Chemistry","volume":"10 1","pages":"57–65"},"PeriodicalIF":2.9,"publicationDate":"2025-12-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://pubs.acs.org/doi/pdf/10.1021/acsearthspacechem.5c00187","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145969069","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-12-18DOI: 10.1021/acsearthspacechem.5c00208
Vittorio Bariosco*, , , Stefano Pantaleone, , , Cecilia Ceccarelli, , , Piero Ugliengo, , and , Albert Rimola*,
Understanding diffusion on interstellar ices is key to modeling the chemical evolution of cold molecular clouds, where low temperatures severely limit molecular mobility. In this study, we introduce a robust and fully automated multiscale computational framework to quantify diffusion processes of adsorbates at the surface of amorphous solid water (ASW). Using H2S as a test case, whose binding sites were previously studied at the ab initio level, we constructed a detailed network of 141 adsorption sites connected by over 270 transition states. All density functional energetics were benchmarked against DLPNO–CCSD(T), achieving chemical accuracy in the description of diffusion barriers, which span from 0.1 to 27 kJ mol–1 with a median value of 5.4 kJ mol–1. An off-lattice kinetic Monte Carlo (kMC) model adopting both the ab initio diffusion barriers and binding energies for the desorption processes was carried out to compute temperature-dependent diffusion coefficients and to reconstruct the temperature-programmed desorption (TPD) curve. Our simulations reveal that thermal diffusion of H2S is negligible below 20 K, with diffusion coefficients as low as 10–48 cm2 s–1 at 10 K, thus excluding Langmuir–Hinshelwood surface encounters under typical dense cloud conditions. Moreover, under submonolayer conditions, diffusion was found to have negligible influence on the reconstructed TPD peak position. Furthermore, our results demonstrate that a universal scaling factor f to guess the diffusion barriers (ΔEdiff) from the sole knowledge of BE: f = ΔEdiff/BE does not apply as it exhibits wide variability across the sampled configurations. These findings highlight the need for incorporating statistically meaningful distributions of binding energies and diffusion barriers in astrochemical models to more accurately capture diffusion and surface reactivity on interstellar ices.
{"title":"A Fully ab Initio Kinetic Monte Carlo Approach for Modeling Adsorption and Diffusion in Interstellar Icy Grain Mantles: The Case of H2S","authors":"Vittorio Bariosco*, , , Stefano Pantaleone, , , Cecilia Ceccarelli, , , Piero Ugliengo, , and , Albert Rimola*, ","doi":"10.1021/acsearthspacechem.5c00208","DOIUrl":"https://doi.org/10.1021/acsearthspacechem.5c00208","url":null,"abstract":"<p >Understanding diffusion on interstellar ices is key to modeling the chemical evolution of cold molecular clouds, where low temperatures severely limit molecular mobility. In this study, we introduce a robust and fully automated multiscale computational framework to quantify diffusion processes of adsorbates at the surface of amorphous solid water (ASW). Using H<sub>2</sub>S as a test case, whose binding sites were previously studied at the ab initio level, we constructed a detailed network of 141 adsorption sites connected by over 270 transition states. All density functional energetics were benchmarked against DLPNO–CCSD(T), achieving chemical accuracy in the description of diffusion barriers, which span from 0.1 to 27 kJ mol<sup>–1</sup> with a median value of 5.4 kJ mol<sup>–1</sup>. An off-lattice kinetic Monte Carlo (kMC) model adopting both the ab initio diffusion barriers and binding energies for the desorption processes was carried out to compute temperature-dependent diffusion coefficients and to reconstruct the temperature-programmed desorption (TPD) curve. Our simulations reveal that thermal diffusion of H<sub>2</sub>S is negligible below 20 K, with diffusion coefficients as low as 10<sup>–48</sup> cm<sup>2</sup> s<sup>–1</sup> at 10 K, thus excluding Langmuir–Hinshelwood surface encounters under typical dense cloud conditions. Moreover, under submonolayer conditions, diffusion was found to have negligible influence on the reconstructed TPD peak position. Furthermore, our results demonstrate that a universal scaling factor <i>f</i> to guess the diffusion barriers (Δ<i>E</i><sub>diff</sub>) from the sole knowledge of BE: <i>f</i> = Δ<i>E</i><sub>diff</sub>/BE does not apply as it exhibits wide variability across the sampled configurations. These findings highlight the need for incorporating statistically meaningful distributions of binding energies and diffusion barriers in astrochemical models to more accurately capture diffusion and surface reactivity on interstellar ices.</p>","PeriodicalId":15,"journal":{"name":"ACS Earth and Space Chemistry","volume":"10 1","pages":"66–79"},"PeriodicalIF":2.9,"publicationDate":"2025-12-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://pubs.acs.org/doi/pdf/10.1021/acsearthspacechem.5c00208","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145969068","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}