We present a major update to the open-source atmospheric modeling package \texttt{PICASO}, designed for simulating the thermal structure and spectra of hydrogen-rich atmospheres of brown dwarfs and exoplanets. This release, \texttt{PICASO 4.0}, expands upon the existing radiative-convective equilibrium model framework by incorporating several new capabilities. Key additions include the integration of \texttt{Virga} for self-consistent cloud modeling, new flexible treatments for rainout and cold trapping of volatile species, and support for photochemistry. We also introduce a parameterized energy injection scheme to simulate additional external or internal heating processes. These features are motivated by lessons from recent JWST observations that reveal the prevalence of non-equilibrium chemistry and clouds. We benchmark the new functionalities against previously published results in the literature, including the Sonora Diamondback grid, energy injected atmospheres, patchy cloud models, and other photochemical models of WASP-39b. \texttt{PICASO} continues to be actively developed as an open-source package aimed at enabling reproducible, community-driven atmospheric modeling of all
The mechanisms by which light interacts with ice and the impact of photo-induced reactions are central to our understanding of environmental, atmospheric and astrophysical processes. However, a microscopic description of the photoproducts originating from UV absorption and emission processes has remained elusive. Here we explore the photochemistry of ice using time-dependent hybrid density functional theory on various models of pristine and defective ice Ih. Our investigation of the excited state potential energy surface of the crystal shows that UV absorption can lead to the formation of hydronium ions, hydroxyl radicals and excess electrons. One of the dominant mechanisms of decay from the excited to the ground-state involves the recombination of the electron with the hydroxyl radical yielding hydronium-hydroxide ion-pairs. We find that the details of this charge recombination process sensitively depend on the presence of defects in the lattice, such as vacancies and pre-existing photoproducts. We also observe the formation of Bjerrum defects following UV absorption; we suggest that, together with hydroxide anions, they are likely responsible for prominent features experimentally
We present an efficient heating/cooling method coupled with chemistry and ultraviolet (UV) radiative transfer, which can be applied to numerical simulations of the interstellar medium (ISM). We follow the time-dependent evolution of hydrogen species (H$_2$, H, H$^+$), assume carbon/oxygen species (C, C$^+$, CO, O, and O$^+$) are in formation-destruction balance given the non-steady hydrogen abundances, and include essential heating/cooling processes needed to capture thermodynamics of all ISM phases. UV radiation from discrete point sources and the diffuse background is followed through adaptive ray tracing and a six-ray approximation, respectively, allowing for H$_2$ self-shielding; cosmic ray (CR) heating and ionization are also included. To validate our methods and demonstrate their application for a range of density, metallicity, and radiation field, we conduct a series of tests, including the equilibrium curves of thermal pressure vs. density, the chemical and thermal structure in photo-dissociation regions, H I-to-H$_2$ transitions, and the expansion of H II regions and radiative supernova remnants. Careful treatment of photochemistry and CR ionization is essential for many a
A detailed modeling of simultaneous UV-photochemical and thermochemical processes in exoplanet atmosphere-like conditions is essential for the analysis and interpretation of a vast amount of current and future spectral data from exoplanets. However, a detailed reaction kinetic model that incorporates both UV photochemistry and thermal chemistry is challenging due to the massive size of the chemical system as well as to the lack of understanding of photochemistry compared to thermal-only chemistry. Here, we utilize an automatic chemical reaction mechanism generator to build a high-fidelity thermochemical reaction kinetic model later then incorporated with UV-photochemistry enhanced by metastable triplet-state carbon monoxide (a3Pi). Our model results show that two different photochemical reactions driven by Lyman-a photons (i.e. H2 + CO(a3Pi) -> H + HCO and CO(X1Sig+) + CO(a3Pi) -> C(3P) + CO2) can enhance thermal chemistry resulting in significant increases in the formation of CH4, H2O, and CO2 in H2-dominated systems with trace amounts of CO, which qualitatively matches with the observations from previous experimental studies. Our model also suggests that at temperatures abo
Some super Earths and mini Neptunes will likely have thick atmospheres that are not H2-dominated. We have developed a photochemistry-thermochemistry kinetic-transport model for exploring the compositions of thick atmospheres on super Earths and mini Neptunes, applicable for both H2-dominated atmospheres and non-H2-dominated atmospheres. Using this model to study thick atmospheres for wide ranges of temperatures and elemental abundances, we classify them into hydrogen-rich atmospheres, water-rich atmospheres, oxygen-rich atmospheres, and hydrocarbon-rich atmospheres. We find that carbon has to be in the form of CO2 rather than CH4 or CO in a H2-depleted water-dominated thick atmosphere, and that the preferred loss of light elements from an oxygen-poor carbon-rich atmosphere leads to formation of unsaturated hydrocarbons (C2H2 and C2H4). We apply our self-consistent atmosphere models to compute spectra and diagnostic features for known transiting low-mass exoplanets GJ 1214 b, HD 97658 b, and 55 Cnc e. For GJ 1214 b like planets we find that (1) C2H2 features at 1.0 and 1.5 micron in transmission and C2H2 and C2H4 features at 9-14 micron in thermal emission are diagnostic for hydroca
HCl was detected in the Martian atmosphere by the NOMAD and ACS spectrometers aboard the ExoMars TGO. Photochemical models show that using gas-phase chemistry alone is insufficient to reproduce these data. Recent work has developed a heterogeneous chemical network within a 1D photochemistry model, guided by the seasonal variability in HCl. The aim of this work is to show that incorporating heterogeneous chlorine chemistry into a global 3D model of Martian photochemistry with conventional gas-phase chemistry can reproduce spatial and temporal changes in hydrogen chloride on Mars. We incorporated this heterogeneous chlorine scheme into the MPCM to model chlorine photochemistry during MYs 34 and 35. These two years provide contrasting dust scenarios, with MY 34 featuring a global dust storm. We also examined correlations in the model results between HCl and other key atmospheric quantities, as well as production and loss processes, to understand the impact of different factors driving changes in HCl. We find that this 3D model of Martian is consistent with the changes in HCl observed by ACS in MY 34 and MY 35, including detections and 70% of non-detections. For the remaining 30%, mode
The interstellar medium is characterized by a rich and diverse chemistry. Many of its complex organic molecules are proposed to form through radical chemistry in icy grain mantles. Radicals form readily when interstellar ices (composed of water and other volatiles) are exposed to UV photons and other sources of dissociative radiation, and, if sufficiently mobile, the radicals can react to form larger, more complex molecules. The resulting complex organic molecules (COMs) accompany star and planet formation, and may eventually seed the origins of life on nascent planets. Experiments of increasing sophistication have demonstrated that known interstellar COMs as well as the prebiotically interesting amino acids can form through ice photochemistry. We review these experiments and discuss the qualitative and quantitative kinetic and mechanistic constraints they have provided. We finally compare the effects of UV radiation with those of three other potential sources of radical production and chemistry in interstellar ices: electrons, ions and X-rays.
Astrochemistry is a well-established multidisciplinary field devoted to study molecules in space. While most astrochemists are oriented to observe molecules in the gas phase and reproduce their abundances by modeling the physical conditions of the medium, the microscopic dust particles wandering in the interstellar medium deserve the attention of a smaller community. Radiation and thermally-driven processes taking place in the bare dust, and particularly in dust particles covered by ice mantles, are mimicked in the laboratory. In addition to water, interstellar ice contains other simple species. In this Review we present our current knowledge on ice photochemistry and thermal processing that ultimately leads to formation of complex organic molecules (COMs). Numerous COMs are of astrobiological interest and match those present in comets and asteroids. Upon impact of these minor bodies, water and COMs were delivered to the earth and might have intervened in the first prebiotic reactions.
This Perspective is part of a Special Topic that explored the maturity of nonadiabatic molecular dynamics for predicting photochemical processes. In 2023, a prediction challenge was issued to the community of computational photochemists to simulate the photochemistry of cyclobutanone, photoexcited at 200 nm, and the resulting time-resolved MeV-UED signal. The challenge attracted 15 theoretical predictions from more than 70 researchers, employing a wide range of strategies for electronic structure and nonadiabatic molecular dynamics to predict the time-resolved MeV-UED signal before the experiment had been conducted at SLAC (Stanford, USA). The MeV-UED instrument at Shanghai Jiao Tong University was also used to provide a second independent time-resolved MeV-UED signal for the photochemistry of cyclobutanone. This Perspective discusses the various approaches and strategies used by the participants to predict the photochemistry of cyclobutanone. This work also summarizes the strengths and weaknesses of various methods used for photoexcitation, electronic structure, nonadiabatic dynamics, and calculation of observables, as agreed by the participants during a CECAM workshop dedicated t
We studied several representative circumstellar disks surrounding the Herbig Ae star MWC 480 and the T Tauri stars LkCa 15 and DM Tau at (sub-)millimeter wavelengths in lines of CCH. Our aim is to characterize photochemistry in the heavily UV-irradiated MWC 480 disk and compare the results to the disks around cooler T Tauri stars. We detected and mapped CCH in these disks with the IRAM Plateau de Bure Interferome- ter in the C- and D-configurations in the (1-0) and (2-1) transitions. Using an iterative minimization technique, the CCH column densities and excitation conditions are con- strained. Very low excitation temperatures are derived for the T Tauri stars. These values are compared with the results of advanced chemical modeling, which is based on a steady-state flared disk structure with a vertical temperature gradient, and a gas- grain chemical network with surface reactions. Both model and observations suggest that CCH is a sensitive tracer of the X-ray and UV irradiation. The predicted radial dependency and source to source variations of CCH column densities qualitatively agree with the observed values, but the predicted column densities are too low by a factor of several.
We investigate the impact of photochemistry and X-ray ionization on the molecular composition of, and ionization fraction in, a protoplanetary disk surrounding a typical T Tauri star. We use a sophisticated physical model, which includes a robust treatment of the radiative transfer of UV and X-ray radiation, and calculate the time-dependent chemical structure using a comprehensive chemical network. In previous work, we approximated the photochemistry and X-ray ionization, here, we recalculate the photoreaction rates using the explicit UV wavelength spectrum and wavelength-dependent reaction cross sections. We recalculate the X-ray ionization rate using our explicit elemental composition and X-ray energy spectrum. We find photochemistry has a larger influence on the molecular composition than X-ray ionization. Observable molecules sensitive to the photorates include OH, HCO+, N2H+, H2O, CO2 and CH3OH. The only molecule significantly affected by the X-ray ionization is N2H+ indicating it is safe to adopt existing approximations of the X-ray ionization rate in typical T Tauri star-disk systems. The recalculation of the photorates increases the abundances of neutral molecules in the ou
In comets the atomic oxygen green to red-doublet emission intensity ratio (G/R ratio) of 0.1 has been used to confirm H$_2$O as the parent species producing oxygen emission lines. The larger ($>$0.1) value of G/R ratio observed in a few comets is ascribed to the presence of higher CO$_2$ and CO relative abundances in the cometary coma. We aim to study the effect of CO$_2$ and CO relative abundances on the observed G/R ratio in comets observed at large ($>$2 au) heliocentric distances by accounting for important production and loss processes of O($^1$S) and O($^1$D) in the cometary coma. Recently we have developed a coupled chemistry-emission model to study photochemistry of O($^1$S) and O($^1$D) atoms and the production of green and red-doublet emissions in comets Hyakutake and Hale-Bopp. In the present work we applied the model to six comets where green and red-doublet emissions are observed when they are beyond 2 au from the Sun. In a water-dominated cometary coma and with significant ($>$10%) CO$_2$ relative abundance, photodissociation of H$_2$O mainly governs the red-doublet emission, whereas CO$_2$ controls the green line emission. If a comet has equal composition of
We present a comprehensive photochemistry model for exploration of the chemical composition of terrestrial exoplanet atmospheres. The photochemistry model is designed from the ground up to have the capacity to treat all types of terrestrial planet atmospheres, ranging from oxidizing through reducing, which makes the code suitable for applications for the wide range of anticipated terrestrial exoplanet compositions. The one-dimensional chemical transport model treats up to 800 chemical reactions, photochemical processes, dry and wet deposition, surface emission and thermal escape of O, H, C, N and S bearing species, as well as formation and deposition of elemental sulfur and sulfuric acid aerosols. We validate the model by computing the atmospheric composition of current Earth and Mars and find agreement with observations of major trace gases in Earth's and Mars' atmospheres. We simulate several plausible atmospheric scenarios of terrestrial exoplanets, and choose three benchmark cases for atmospheres from reducing to oxidizing. The most interesting finding is that atomic hydrogen is always a more abundant reactive radical than the hydroxyl radical in anoxic atmospheres. Whether ato
Energetic processing of interstellar ice mantles and planetary atmospheres via photochemistry is a critical mechanism in the extraterrestrial synthesis of prebiotic molecules. Photochemistry is defined as chemical processes initiated by photon-induced electronic excitation, not involving ionization. In contrast, photons with energies above the ionization threshold initiate radiation chemistry (radiolysis). Vacuum-ultraviolet (6.2-12.4 eV) light may initiate photochemistry and radiation chemistry because the threshold for producing secondary electrons is lower in the condensed phase than in the gas phase. Approximately half of cosmic-ray induced photons incident on interstellar ices in star-forming regions initiate photochemistry while the rest initiate radiation chemistry. While experimental techniques such as velocity map imaging may be used to extract exquisite details about gas-phase photochemistry, such detailed information cannot be obtained for condensed-phase photochemistry, which involves greater complexity, including the production of excitons, excimers, and exciplexes. Because a primary objective of chemistry is to provide molecular-level mechanistic explanations for macr
Sulfur gases are common components in the volcanic and biological emission on Earth, and are expected to be important input gases for atmospheres on terrestrial exoplanets. We study the atmospheric composition and the spectra of terrestrial exoplanets with sulfur compounds (i.e., H2S and SO2) emitted from their surfaces. We use a comprehensive one-dimensional photochemistry model and radiative transfer model to investigate the sulfur chemistry in atmospheres ranging from reducing to oxidizing. The most important finding is that both H2S and SO2 are chemically short-lived in virtually all types of atmospheres on terrestrial exoplanets, based on models of H2, N2, and CO2 atmospheres. This implies that direct detection of surface sulfur emission is unlikely, as their surface emission rates need to be extremely high (>1000 times Earth's volcanic sulfur emission) for these gases to build up to a detectable level. We also find that sulfur compounds emitted from the surface lead to photochemical formation of elemental sulfur and sulfuric acid in the atmosphere, which would condense to form aerosols if saturated. For terrestrial exoplanets in the habitable zone of Sun-like stars or M st
Mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT) has recently emerged as an attractive electronic-structure method for studying photochemical processes, given that it bridges the computational efficiency of single-reference approaches with the versatility of multireference methods. In the following, we critically assess the general applicability of MRSF-TDDFT to photochemistry and identify two important limitations. First, the doubly-excited configurations included in MRSF-TDDFT come at the cost of missing some singly-excited configurations. Second, MRSF-TDDFT provides unreliable excited-state energies when its triplet reference - a cornerstone of the method - abruptly changes its nature, e.g., when the T$_1$ and T$_2$ triplet states become nearly degenerate and exchange electronic character. This change of character of the triplet reference can induce discontinuities or sharp distortions in electronic potential energy curves of the response states in unsuspected regions of the nuclear configuration space. We propose strategies and diagnostics to detect these limitations in the exploration of potential energy surfaces and nonadiabatic molecular dynami
Atmospheric nitrogen may provide important constraints on giant planet formation. Following our semi-analytical work (Ohno & Fortney 2022), we further pursue the relation between observable NH3 and an atmosphere's bulk nitrogen abundance by applying the photochemical kinetics model VULCAN across planetary equilibrium temperature, mass, age, eddy diffusion coefficient, atmospheric composition, and stellar spectral type. We confirm that the quenched NH3 abundance coincides with the bulk nitrogen abundance only at sub-Jupiter mass (< 1MJ) planets and old ages (> 1 Gyr) for solar composition atmospheres, highlighting important caveats for inferring atmospheric nitrogen abundances. Our semi-analytical model reproduces the quenched NH3 abundance computed by VULCAN and thus helps to infer the bulk nitrogen abundance from a retrieved NH3 abundance. By computing transmission and emission spectra, we predict that the equilibrium temperature range of 400--1000 K is optimal for detecting NH3 because NH3 depletion by thermochemistry and photochemistry is significant at hotter planets whereas entire spectral features become weak at colder planets. For Jupiter-mass planets around Sun-li
Graphitic carbon nitride ($g$-CN) has attracted vast interest as a promising inexpensive metal-free photocatalyst for water splitting with solar photons. The heptazine (Hz) molecule is the building block of graphitic carbon nitride. The photochemistry of the Hz molecule and derivatives thereof in protic environments has been the subject of several recent experimental and computational studies. In the present work, the hydrogen-bonded Hz$\cdots$H$_2$O complex was adopted as a model system for the exploration of photoinduced electron and proton transfer processes in this complex with quasi-classical nonadiabatic trajectory simulations, using the $ab$ $initio$ ADC(2) electronic-structure method and a computationally efficient surface-hopping algorithm. The population of the optically excited bright $^1ππ^*$ state of the Hz chromophore relaxes through three $^1nπ^*$ states and a low-lying charge-transfer state, which drives proton transfer from H$_2$O to Hz, to the long-lived optically dark S$_1$($ππ^*$) state of Hz. The imaging of this ultrafast and complex dynamics with femtosecond time-resolved transient absorption (TA) pump-probe (PP) spectroscopy and two-dimensional (2D) electroni
Compared to the diversity seen in exoplanets, Venus is a veritable astrophysical twin of the Earth, however its global cloud layer truncates features in transmission spectroscopy, masking its non-Earth-like nature. Observational indicators that can distinguish an exo-Venus from an exo-Earth must therefore survive above the cloud layer. The above-cloud atmosphere is dominated by photochemistry, which depends on the spectrum of the host star and therefore changes between stellar systems. We explore the systematic changes in photochemistry above the clouds of Venus-like exoplanets orbiting K-Dwarf or M-Dwarf host stars, using a recently validated model of the full Venus atmosphere (0-115 km) and stellar spectra from the MUSCLES Treasury survey. SO2, OCS and H2S are key gas species in Venus-like planets that are not present in Earth-like planets, and could therefore act as observational discriminants if their atmospheric abundances are high enough to be detected. We find that SO2, OCS and H2S all survive above the cloud layer when irradiated by the coolest K-Dwarf and all seven M-Dwarfs, whereas these species are heavily photochemically depleted above the clouds of Venus. The productio
Ice mantles that formed on top of dust grains are photoprocessed by the secondary ultraviolet (UV) field in cold and dense molecular clouds. UV photons induce photochemistry and desorption of ice molecules. Experimental simulations dedicated to ice analogs under astrophysically relevant conditions are needed to understand these processes. We present UV-irradiation experiments of a pure CO2 ice analog. Calibration of the QMS allowed us to quantify the photodesorption of molecules to the gas phase. This information was added to the data provided by the FTIR on the solid phase to obtain a complete quantitative study of the UV photoprocessing of an ice analog. Experimental simulations were performed in an ultra-high vacuum chamber. Ice samples were deposited onto an infrared transparent window at 8K and were subsequently irradiated with a microwave-discharged hydrogen flow lamp. After irradiation, ice samples were warmed up until complete sublimation was attained. Photolysis of CO2 molecules initiates a network of photon-induced chemical reactions leading to the formation of CO, CO3 ,O2 , and O3 . During irradiation, photon-induced desorption of CO and, to a lesser extent, O2 and CO2 t