DRAFT VERSION AUGUST 5, 2021 Typeset using LATEX twocolumn style in AASTeX63 A Comparative Study of Atmospheric Chemistry with VULCAN SHANG-MIN TSAI,1 MATEJ MALIK,2 DANIEL KITZMANN,3 JAMES R. LYONS,4 ALEXANDER FATEEV,5 ELSPETH LEE,6 AND KEVIN HENG6 1Atmospheric, Ocean, and Planetary Physics, Department of Physics, University of Oxford, UK 2Department of Astronomy, University of Maryland, College Park, MD 20742, USA 3University of Bern, Center for Space and Habitability 4Arizona State University, School of Earth and Space Exploration, Bateman Physical Sciences, Tempe, USA 5Technical University of Denmark, Department of Chemical and Biochemical, Denmark 6University of Bern, Center for Space and Habitability, Switzerland ABSTRACT We present an update of the open-source photochemical kinetics code VULCAN (Tsai et al.(2017); https://github.com/exoclime/VULCAN) to include C-H-N-O-S networks and photochemistry. Additional new features are advection transport, condensation, various boundary conditions, and temperature-dependent UV cross-sections. First, we validate our photochemical model for hot Jupiter atmospheres by performing an in- tercomparison of HD 189733b models between Moses et al.(2011), Venot et al.(2012), and VULCAN, to diagnose possible sources of discrepancy. Second, we set up a model of Jupiter extending from the deep tro- posphere to upper stratosphere to verify the kinetics for low temperature. Our model reproduces hydrocarbons consistent with observations, and the condensation scheme successfully predicts the locations of water and am- monia ice clouds. We show that vertical advection can regulate the local ammonia distribution in the deep atmosphere. Third, we validate the model for oxidizing atmospheres by simulating Earth and find agreement with observations. Last, VULCAN is applied to four representative cases of extrasolar giant planets: WASP- 33b, HD 189733b, GJ 436b, and 51 Eridani b. We look into the effects of the C/O ratio and chemistry of titanium/vanadium species for WASP-33b; we revisit HD 189733b for the effects of sulfur and carbon con- densation; the effects of internal heating and vertical mixing (Kzz) are explored for GJ 436b; we test updated planetary properties for 51 Eridani b with S8 condensates. We find sulfur can couple to carbon or nitrogen and impact other species such as hydrogen, methane, and ammonia. The observable features of the synthetic spectra and trends in the photochemical haze precursors are discussed for each case. 1. INTRODUCTION Birkby et al. 2013; Brogi et al. 2014; Hoeijmakers et al. Understanding the chemical compositions has been a cen- 2018). tral aspect in atmospheric characterization for planets within The majority of discovered exoplanets have sizes between and beyond the Solar System. Photochemical kinetics mod- Earth and Neptune. Their heavy elemental abundances (i.e. els establish the link between our knowledge of chemical metallicity) can vary considerably, as often inferred by the reactions and various planetary processes (e.g., atmospheric water detection (e.g., Wakeford et al. 2017; Chachan et al. dynamics, radiative transfer, outgassing process, etc.), pro- 2019). While CH4 is expected to be more abundant in cooler viding a theoretical basis for interpreting observations and (Teq . 1000 K) atmospheres, understanding how disequilib- addressing habitability. rium chemistry and other processes alter the CH4/CO abun- arXiv:2108.01790v1 [astro-ph.EP] 3 Aug 2021 Hot Jupiters are the first discovered and best characterized dance ratio remains an ongoing task. class of exoplanets. Transit and eclipse observations have The direct imaging technique provides a complementary made various initial detections of chemical species in their window to resolve young planets at a far orbit (e.g., see the reviews of Crossfield 2015; Pueyo 2018)). The new gener- atmospheres such as Na, K, H2O, CH4, CO, CO2 (e.g. see the review of Kreidberg 2018). An extreme class of exceed- ation of instruments like GPI and SPHERE (Chauvin 2018) ingly irradiated hot Jupiter around bright stars have equilib- have identified a number of interesting young Jupiter ana- rium temperature higher than 2000 K. They are prime tar- log. These young planets are self-luminous from their heat gets for emission observations, and recent high-resolution of formation and receive UV fluxes from the star at the same spectroscopic measurements reveal atomic and ionic features time, giving insights on the planet forming conditions outside that make their atmospheres resemble low-mass stars (e.g., the snow lines and the transition between planets and brown dwarfs. 2 Across the various types of aforementioned planetary at- cles are possibly ubiquitous, with diverse compositions (Gao mospheres, photochemical kinetics and atmospheric trans- et al. 2020) including cloud particles formed from conden- port are the dominant mechanisms that control the ma- sation or produced by photolysis at high altitudes. Micro- jor chemical abundances. Photodissociation occurs when physics models (Helling & Woitke 2006; Lavvas & Koski- molecules are split into reactive radicals by high-energy pho- nen 2017; Kawashima & Ikoma 2018; Gao & Benneke 2018; tons while atmospheric transport shapes the abundance distri- Ohno et al. 2020) have investigated trends and properties of bution. Disequilibrium processes can drive abundances con- aerosols for various environments. One particularly interest- siderably away from the chemical equilibrium state and are ing candidate of aerosols is the sulfur family, such as sulfuric best studied in chemical kinetics models. clouds (Hu et al. 2013; Misra et al. 2015; Loftus et al. 2019) Kinetics models stem from simulating the atmospheric in an oxidizing atmosphere or elemental sulfur in a reducing compositions in Solar System planets (e.g., Kasting et al. atmosphere (Hu et al. 2013; Gao et al. 2017). Photochem- 1979; Yung et al. 1984; Nair et al. 1994; Wilson & Atreya istry generally sets off the initial steps in the gas phase, then 2004; Lavvas et al. 2008; Hu et al. 2012; Krasnopolsky the condensable species can form particles when saturated in 2012), which focus on photochemistry and radical reactions. a broad range of altitudes (Gao et al. 2017). The relatively The low temperature regime makes thermochemistry less rel- simple sulfur particles in H2-dominated atmospheres allow a evant in most cases. Liang et al.(2003) first applied a pho- consistent photochemical-aerosol kinetics modeling, which tochemical kinetics model, Caltech/JPL KINETICS (Allen we will conduct in this work. Although the formation path- et al. 1981), to the hot Jupiter HD 209458b and identified ways of organic haze particles are highly complex, we will the photochemical source of water for producing atomic H. focus on a group of haze precursors and investigate their pho- However, some reaction rates in their study are extrapolated tochemical stability in the hope of providing complementary from measurements at low temperatures and not suitable for insights on the haze-forming conditions. hot Jupiter conditions. Line et al.(2010) adopt the high- The exclusive access to often proprietary chemical mod- temperature rate coefficients for the major molecules and els motivates us to develop an open-source, chemical kinet- use the lower boundary to mimic mixing from the thermo- ics code VULCAN (Tsai et al. 2017). The initial version chemical equilibrium region. A new group of models incor- of VULCAN includes a reduced-size C-H-O thermochemical porating kinetics data valid at high temperatures started to network and treats eddy diffusion. In Tsai et al.(2017), VUL- emerge since then. Zahnle et al.(2009) reverse the reactions CAN is validated by comparing the quench behavior with to ensure kinetics consistent with thermodynamic calcula- ARGO (Rimmer & Helling 2016) and Moses et al.(2011). tions and consider sulfur chemistry on hot Jupiters. Moses Since then, VULCAN has been continuously updated and ap- et al.(2011) implement high-temperature reactions in KI- plied to several studies such as Zilinskas et al.(2020) who NETICS to model hot Jupiters HD 189733b and HD 209458b identify key molecules of hot super-Earths with nitrogen- with detailed pathway analysis. Venot et al.(2012) adopt the dominated atmospheres, and Shulyak et al.(2020) who ex- combustion mechanisms validated for industrial applications plore the effects of XUV for different stellar types. to model the same canonical hot Jupiters but find different In this work, we present the new version of 1-D photo- quenching and photolysis profiles from Moses et al.(2011). chemical model VULCAN, with embedded chemical net- Hobbs et al.(2021) recently extend Zahnle et al.(2009) to works now including hydrogen, oxygen, carbon, nitrogen, include sulfur photochemistry and find the inclusion of sul- and sulfur. The chemical network is customizable and does fur can impact other non-sulfur species on HD 209458b and not require separating fast and slow species. The major up- 51 Eridani b. As the discovery of diverse exoplanets pro- dates of VULCAN from Tsai et al.(2017) are: gresses, more kinetics models have been applied to study a wide range of aspects, such as the compositional diversity • C-H-N-O-S chemical networks with about 100 within an atmospheric-grid framework (Moses et al. 2013; species, including a simplified benzene forming mech- Miguel & Kaltenegger 2014; Molaverdikhani et al. 2019), anism atmospheric evolution with loss and/or outgassing processes (Hu et al. 2015; Wordsworth et al. 2018; Lincowski