Adsorption of Polycyclic Aromatic Hydrocarbons and C60onto Forsterite: C-H Bond Activation by the Schottky Vacancy

Understanding how to catalytically break the C-H bond of aromatic molecules, such as polycyclic aromatic hydrocarbons (PAHs), is currently a big challenge and a subject of study in catalysis, astrochemistry, and planetary science. In the latter, the study of the breakdown reaction of PAHs on mineral surfaces is important to understand if PAHs are linked to prebiotic molecules in regions of star and planet formation. In this work, we employed a periodic density functional theory along with Grimme's D4 (DFT-D4) approach for studying the adsorption of a sample of PAHs (naphthalene, anthracene, fluoranthene, pyrene, coronene, and benzocoronene) and fullerene on the [010] forsterite surface and its defective surfaces (Fe-doped and Ni-doped surfaces and a MgO-Schottky vacancy) for their implications in catalysis and astrochemistry. On the basis of structural and binding energy analysis, large PAHs and fullerene present stronger adsorption on the pristine, Fe-doped, and Ni-doped forsterite surfaces than small PAHs. On a MgO-Schottky vacancy, parallel adsorption of the PAH leads to the chemisorption process (C-Si and/or C-O bonds), whereas perpendicular orientation of the PAH leads to the catalytic breaking of the aromatic C-H bond via a barrierless reaction. Spin density and charge analysis show that C-H dissociation is promoted by electron donation from the vacancy to the PAH. As a result of the undercoordinated Si and O atoms, the vacancy acts as a Frustrated Lewis Pair (FLP) catalyst. Therefore, a MgO-Schottky vacancy [010] forsterite surface proved to have potential catalytic activity for the activation of C-H bond in aromatic molecules

PDRs4All: A JWST Early Release Science Program on Radiative Feedback from Massive Stars

22 pags., 8 figs., 1 tab.Massive stars disrupt their natal molecular cloud material through radiative and mechanical feedback processes. These processes have profound effects on the evolution of interstellar matter in our Galaxy and throughout the universe, from the era of vigorous star formation at redshifts of 1-3 to the present day. The dominant feedback processes can be probed by observations of the Photo-Dissociation Regions (PDRs) where the far-ultraviolet photons of massive stars create warm regions of gas and dust in the neutral atomic and molecular gas. PDR emission provides a unique tool to study in detail the physical and chemical processes that are relevant for most of the mass in inter-and circumstellar media including diffuse clouds, proto-planetary disks, and molecular cloud surfaces, globules, planetary nebulae, and star-forming regions. PDR emission dominates the infrared (IR) spectra of star-forming galaxies. Most of the Galactic and extragalactic observations obtained with the James Webb Space Telescope (JWST) will therefore arise in PDR emission. In this paper we present an Early Release Science program using the MIRI, NIRSpec, and NIRCam instruments dedicated to the observations of an emblematic and nearby PDR: the Orion Bar. These early JWST observations will provide template data sets designed to identify key PDR characteristics in JWST observations. These data will serve to benchmark PDR models and extend them into the JWST era. We also present the Science-Enabling products that we will provide to the community. These template data sets and Science-Enabling products will guide the preparation of future proposals on star-forming regions in our Galaxy and beyond and will facilitate data analysis and interpretation of forthcoming JWST observations.Support for JWST-ERS program ID 1288 was provided through grants from the STScI under NASA contract NAS5-03127 to STScI (K.G., D.V.D.P., M.R.), Univ. of Maryland (M.W., M.P.), Univ. of Michigan (E.B., F.A.), and Univ. of Toledo (T.S.-Y.L.). O.B. and E.H. are supported by the Programme National “Physique et Chimie du Milieu Interstellaire” (PCMI) of CNRS/INSU with INC/INP co-funded by CEA and CNES, and through APR grants 6315 and 6410 provided by CNES. E. P. and J.C. acknowledge support from the National Science and Engineering Council of Canada (NSERC) Discovery Grant program (RGPIN-2020-06434 and RGPIN-2021-04197 respectively). E.P. acknowledges support from a Western Strategic Support Accelerator Grant (ROLA ID 0000050636). J.R.G. and S.C. thank the Spanish MCINN for funding support under grant PID2019-106110GB-I00. Work by M.R. and Y.O. is carried out within the Collaborative Research Centre 956, subproject C1, funded by the Deutsche Forschungsgemeinschaft (DFG)—project ID 184018867. T.O. acknowledges support from JSPS Bilateral Program, grant No. 120219939. M.P. and M.W. acknowledge support from NASA Astrophysics Data Analysis Program award #80NSSC19K0573. C.B. is grateful for an appointment at NASA Ames Research Center through the San José State University Research Foundation (NNX17AJ88A) and acknowledges support from the Internal Scientist Funding Model (ISFM) Directed Work Package at NASA Ames titled: “Laboratory Astrophysics—The NASA Ames PAH IR Spectroscopic Database.”Peer reviewe

Interaction of Aromatic Molecules with Forsterite: Accuracy of the Periodic DFT-D4 Method

Density functional theory (DFT) has provided deep atomic-level insights into the adsorption behavior of aromatic molecules on solid surfaces. However, modeling the surface phenomena of large molecules on mineral surfaces with accurate plane wave methods (PW) can be orders of magnitude more computationally expensive than localized atomic orbitals (LCAO) methods. In the present work, we propose a less costly approach based on the DFT-D4 method (PBE-D4), using LCAO, to study the interactions of aromatic molecules with the {010} forsterite (Mg2SiO4) surface for their relevance in astrochemistry. We studied the interaction of benzene with the pristine {010} forsterite surface and with transition-metal cations (Fe2+ and Ni2+) using PBE-D4 and a vdW-inclusive density functional (Dion, Rydberg, Schröder, Langreth, and Lundqvist (DRSLL)) with LCAO methods. PBE-D4 shows good agreement with coupled-cluster methods (CCSD(T)) for the binding energy trend of cation complexes and with PW methods for the binding energy of benzene on the forsterite surface with a difference of about 0.03 eV. The basis set superposition error (BSSE) correction is shown to be essential to ensure a correct estimation of the binding energies even when large basis sets are employed for single-point calculations of the optimized structures with smaller basis sets. We also studied the interaction of naphthalene and benzocoronene on pristine and transition-metal-doped {010} forsterite surfaces as a test case for PBE-D4. Yielding results that are in good agreement with the plane wave methods with a difference of about 0.02-0.17 eV, the PBE-D4 method is demonstrated to be effective in unraveling the binding structures and the energetic trends of aromatic molecules on pristine and transition-metal-doped forsterite mineral surfaces. Furthermore, PBE-D4 results are in good agreement with its predecessor PBE-D3(BJM) and with the vdW-inclusive density functionals, as long as transition metals are not involved. Hence, PBE-D4/CP-DZP has been proven to be a robust theory level to study the interaction of aromatic molecules on mineral surfaces

Adsorption of Polycyclic Aromatic Hydrocarbons and C60onto Forsterite: C-H Bond Activation by the Schottky Vacancy

Understanding how to catalytically break the C-H bond of aromatic molecules, such as polycyclic aromatic hydrocarbons (PAHs), is currently a big challenge and a subject of study in catalysis, astrochemistry, and planetary science. In the latter, the study of the breakdown reaction of PAHs on mineral surfaces is important to understand if PAHs are linked to prebiotic molecules in regions of star and planet formation. In this work, we employed a periodic density functional theory along with Grimme's D4 (DFT-D4) approach for studying the adsorption of a sample of PAHs (naphthalene, anthracene, fluoranthene, pyrene, coronene, and benzocoronene) and fullerene on the [010] forsterite surface and its defective surfaces (Fe-doped and Ni-doped surfaces and a MgO-Schottky vacancy) for their implications in catalysis and astrochemistry. On the basis of structural and binding energy analysis, large PAHs and fullerene present stronger adsorption on the pristine, Fe-doped, and Ni-doped forsterite surfaces than small PAHs. On a MgO-Schottky vacancy, parallel adsorption of the PAH leads to the chemisorption process (C-Si and/or C-O bonds), whereas perpendicular orientation of the PAH leads to the catalytic breaking of the aromatic C-H bond via a barrierless reaction. Spin density and charge analysis show that C-H dissociation is promoted by electron donation from the vacancy to the PAH. As a result of the undercoordinated Si and O atoms, the vacancy acts as a Frustrated Lewis Pair (FLP) catalyst. Therefore, a MgO-Schottky vacancy [010] forsterite surface proved to have potential catalytic activity for the activation of C-H bond in aromatic molecules

Interaction of Aromatic Molecules with Forsterite: Accuracy of the Periodic DFT-D4 Method

Density functional theory (DFT) has provided deep atomic-level insights into the adsorption behavior of aromatic molecules on solid surfaces. However, modeling the surface phenomena of large molecules on mineral surfaces with accurate plane wave methods (PW) can be orders of magnitude more computationally expensive than localized atomic orbitals (LCAO) methods. In the present work, we propose a less costly approach based on the DFT-D4 method (PBE-D4), using LCAO, to study the interactions of aromatic molecules with the {010} forsterite (Mg2SiO4) surface for their relevance in astrochemistry. We studied the interaction of benzene with the pristine {010} forsterite surface and with transition-metal cations (Fe2+ and Ni2+) using PBE-D4 and a vdW-inclusive density functional (Dion, Rydberg, Schröder, Langreth, and Lundqvist (DRSLL)) with LCAO methods. PBE-D4 shows good agreement with coupled-cluster methods (CCSD(T)) for the binding energy trend of cation complexes and with PW methods for the binding energy of benzene on the forsterite surface with a difference of about 0.03 eV. The basis set superposition error (BSSE) correction is shown to be essential to ensure a correct estimation of the binding energies even when large basis sets are employed for single-point calculations of the optimized structures with smaller basis sets. We also studied the interaction of naphthalene and benzocoronene on pristine and transition-metal-doped {010} forsterite surfaces as a test case for PBE-D4. Yielding results that are in good agreement with the plane wave methods with a difference of about 0.02-0.17 eV, the PBE-D4 method is demonstrated to be effective in unraveling the binding structures and the energetic trends of aromatic molecules on pristine and transition-metal-doped forsterite mineral surfaces. Furthermore, PBE-D4 results are in good agreement with its predecessor PBE-D3(BJM) and with the vdW-inclusive density functionals, as long as transition metals are not involved. Hence, PBE-D4/CP-DZP has been proven to be a robust theory level to study the interaction of aromatic molecules on mineral surfaces

SOFIA Far-infrared [O III] and [O I] Observations of Dense CO-knots in the Supernova Remnant Cassiopeia A: Multi-phase Ejecta

Dense, fast-moving ejecta knots in supernova remnants are prime sites for molecule and dust formation. We present SOFIA far-IR spectrometer FIFI-LS observations of CO-rich knots in Cas A which cover a ~1 square arc minute area of the northern shell, in the [O III] 52 and 88 micron and [O I] 63 micron lines. The FIFI-LS spectra reveal that the line profiles of [O III] and [O I] are similar to those of the Herschel PACS [O III] and CO lines. We find that the [O III] maps show very different morphology than the [O I] map. The [O III] maps reveal diffuse, large-scale structures and the ratio of the two [O III] lines imply the presence of gas with a range of density 500 - 10,000 per cm^3 within the mapped region. In contrast, the [O I] map shows bright emission associated with the dense CO-rich knots. The 63 micron [O I] line traces cooled, dense post-shocked gas of ejecta. We find that IR-dominated [O III] emission is from post-shocked gas based on its morphology, high column density, and velocity profile. We describe multi-phase ejecta knots, a lifetime of clumps, and survival of dust in the young supernova remnants