-
BOWIE-ALIGN: How formation and migration histories of giant planets impact atmospheric compositions
Authors:
Anna B. T. Penzlin,
Richard A. Booth,
James Kirk,
James E. Owen,
Eva-Maria Ahrer,
Duncan A. Christie,
Alastair B. Claringbold,
Emma Esparza-Borges,
M. López-Morales,
N. J. Mayne,
Mason McCormack,
Annabella Meech,
Vatsal Panwar,
Diana Powell,
Denis E. Sergeev,
Jake Taylor,
Peter J. Wheatley,
Maria Zamyatina
Abstract:
Hot Jupiters present a unique opportunity for measuring how planet formation history shapes present-day atmospheric composition. However, due to the myriad pathways influencing composition, a well-constructed sample of planets is needed to determine whether formation history can be accurately traced back from atmospheric composition. To this end, the BOWIE-ALIGN survey will compare the composition…
▽ More
Hot Jupiters present a unique opportunity for measuring how planet formation history shapes present-day atmospheric composition. However, due to the myriad pathways influencing composition, a well-constructed sample of planets is needed to determine whether formation history can be accurately traced back from atmospheric composition. To this end, the BOWIE-ALIGN survey will compare the compositions of 8 hot Jupiters around F stars, 4 with orbits aligned with the stellar rotation axis and 4 misaligned. Using the alignment as an indicator for planets that underwent disc migration or high-eccentricity migration, one can determine whether migration history produces notable differences in composition between the two samples of planets. This paper describes the planet formation model that motivates our observing programme. Our model traces the accretion of chemical components from the gas and dust in the disc over a broad parameter space to create a full, unbiased model sample from which we can estimate the range of final atmospheric compositions. For high metallicity atmospheres (O/H > 10 times solar), the C/O ratios of aligned and misaligned planets diverge, with aligned planets having lower C/O (< 0.25) due to the accretion of oxygen-rich silicates from the inner disc. However, silicates may rain out instead of releasing their oxygen into the atmosphere. This would significantly increase the C/O of aligned planets (C/O > 0.6), inverting the trend between the aligned and misaligned planets. Nevertheless, by comparing statistically significant samples of aligned and misaligned planets, we expect atmospheric composition to constrain how planets form.
△ Less
Submitted 4 July, 2024; v1 submitted 3 July, 2024;
originally announced July 2024.
-
BOWIE-ALIGN: A JWST comparative survey of aligned vs misaligned hot Jupiters to test the dependence of atmospheric composition on migration history
Authors:
James Kirk,
Eva-Maria Ahrer,
Anna B. T. Penzlin,
James E. Owen,
Richard A. Booth,
Lili Alderson,
Duncan A. Christie,
Alastair B. Claringbold,
Emma Esparza-Borges,
Chloe E. Fisher,
Mercedes López-Morales,
N. J. Mayne,
Mason McCormack,
Annabella Meech,
Vatsal Panwar,
Diana Powell,
Jake Taylor,
Denis E. Sergeev,
Daniel Valentine,
Hannah R. Wakeford,
Peter J. Wheatley,
Maria Zamyatina
Abstract:
A primary objective of exoplanet atmosphere characterisation is to learn about planet formation and evolution, however, this is challenged by degeneracies. To determine whether differences in atmospheric composition can be reliably traced to differences in evolution, we are undertaking a new survey with JWST to compare the compositions of a sample of hot Jupiters that orbit F stars above the Kraft…
▽ More
A primary objective of exoplanet atmosphere characterisation is to learn about planet formation and evolution, however, this is challenged by degeneracies. To determine whether differences in atmospheric composition can be reliably traced to differences in evolution, we are undertaking a new survey with JWST to compare the compositions of a sample of hot Jupiters that orbit F stars above the Kraft break with different orbital alignments. Under the assumption that aligned planets migrate through the inner disc, while misaligned planets migrate after disc dispersal, the act of migrating through the inner disc should lead to a measurable difference in the C/O between aligned and misaligned planets. We expect the amplitude and sign of this difference to depend on the amount of planetesimal accretion and whether silicates accreted from the inner disc release their oxygen. Here, we identify all known exoplanets that are suitable for testing this hypothesis, describe our JWST survey, and use noise simulations and atmospheric retrievals to estimate our survey's sensitivity. With the selected sample of four aligned and four misaligned hot Jupiters, we will be sensitive to the predicted differences in C/O between aligned and misaligned hot Jupiters for a wide range of model scenarios.
△ Less
Submitted 3 July, 2024;
originally announced July 2024.
-
The CUISINES Framework for Conducting Exoplanet Model Intercomparison Projects, Version 1.0
Authors:
Linda E. Sohl,
Thomas J. Fauchez,
Shawn Domagal-Goldman,
Duncan A. Christie,
Russell Deitrick,
Jacob Haqq-Misra,
C. E. Harman,
Nicolas Iro,
Nathan J. Mayne,
Kostas Tsigaridis,
Geronimo L. Villanueva,
Amber V. Young,
Guillaume Chaverot
Abstract:
As JWST begins to return observations, it is more important than ever that exoplanet climate models can consistently and correctly predict the observability of exoplanets, retrieval of their data, and interpretation of planetary environments from that data. Model intercomparisons play a crucial role in this context, especially now when few data are available to validate model predictions. The CUIS…
▽ More
As JWST begins to return observations, it is more important than ever that exoplanet climate models can consistently and correctly predict the observability of exoplanets, retrieval of their data, and interpretation of planetary environments from that data. Model intercomparisons play a crucial role in this context, especially now when few data are available to validate model predictions. The CUISINES Working Group of NASA's Nexus for Exoplanet System Science (NExSS) supports a systematic approach to evaluating the performance of exoplanet models, and provides here a framework for conducting community-organized exoplanet Model Intercomparison Projects (exoMIPs). The CUISINES framework adapts Earth climate community practices specifically for the needs of exoplanet researchers, encompassing a range of model types, planetary targets, and parameter space studies. It is intended to help researchers to work collectively, equitably, and openly toward common goals. The CUISINES framework rests on five principles: 1) Define in advance what research question(s) the exoMIP is intended to address. 2) Create an experimental design that maximizes community participation, and advertise it widely. 3) Plan a project timeline that allows all exoMIP members to participate fully. 4) Generate data products from model output for direct comparison to observations. 5) Create a data management plan that is workable in the present and scalable for the future. Within the first years of its existence, CUISINES is already providing logistical support to 10 exoMIPs, and will continue to host annual workshops for further community feedback and presentation of new exoMIP ideas.
△ Less
Submitted 13 June, 2024;
originally announced June 2024.
-
Longitudinal Filtering, Sponge Layers, and Equatorial Jet Formation in a General Circulation Model of Gaseous Exoplanets
Authors:
D. A. Christie,
N. J. Mayne,
M. Zamyatina,
H. Baskett,
T. M. Evans-Soma,
N. Wood,
K. Kohary
Abstract:
General circulation models are a useful tool in understanding the three dimensional structure of hot Jupiter and sub-Neptune atmospheres; however, understanding the validity of the results from these simulations requires an understanding the artificial dissipation required for numerical stability. In this paper, we investigate the impact of the longitudinal filter and vertical ``sponge'' used in t…
▽ More
General circulation models are a useful tool in understanding the three dimensional structure of hot Jupiter and sub-Neptune atmospheres; however, understanding the validity of the results from these simulations requires an understanding the artificial dissipation required for numerical stability. In this paper, we investigate the impact of the longitudinal filter and vertical ``sponge'' used in the Met Office's {\sc Unified Model} when simulating gaseous exoplanets. We demonstrate that excessive dissipation can result in counter-rotating jets and a catastrophic failure to conserve angular momentum. Once the dissipation is reduced to a level where a super-rotating jet forms, however, the jet and thermal structure are relatively insensitive to the dissipation, except in the nightside gyres where temperatures can vary by $\sim 100\,\mathrm{K}$. We do find, however, that flattening the latitudinal profile of the longitudinal filtering alters the results more than a reduction in the strength of the filtering itself. We also show that even in situations where the temperatures are relatively insensitive to the dissipation, the vertical velocities can still vary with the dissipation, potentially impacting physical processes that depend on the local vertical transport.
△ Less
Submitted 4 June, 2024;
originally announced June 2024.
-
JWST Reveals CH$_4$, CO$_2$, and H$_2$O in a Metal-rich Miscible Atmosphere on a Two-Earth-Radius Exoplanet
Authors:
Björn Benneke,
Pierre-Alexis Roy,
Louis-Philippe Coulombe,
Michael Radica,
Caroline Piaulet,
Eva-Maria Ahrer,
Raymond Pierrehumbert,
Joshua Krissansen-Totton,
Hilke E. Schlichting,
Renyu Hu,
Jeehyun Yang,
Duncan Christie,
Daniel Thorngren,
Edward D. Young,
Stefan Pelletier,
Heather A. Knutson,
Yamila Miguel,
Thomas M. Evans-Soma,
Caroline Dorn,
Anna Gagnebin,
Jonathan J. Fortney,
Thaddeus Komacek,
Ryan MacDonald,
Eshan Raul,
Ryan Cloutier
, et al. (6 additional authors not shown)
Abstract:
Even though sub-Neptunes likely represent the most common outcome of planet formation, their natures remain poorly understood. In particular, planets near 1.5-2.5$\,R_\oplus$ often have bulk densities that can be explained equally well with widely different compositions and interior structures, resulting in grossly divergent implications for their formation. Here, we present the full 0.6-5.2…
▽ More
Even though sub-Neptunes likely represent the most common outcome of planet formation, their natures remain poorly understood. In particular, planets near 1.5-2.5$\,R_\oplus$ often have bulk densities that can be explained equally well with widely different compositions and interior structures, resulting in grossly divergent implications for their formation. Here, we present the full 0.6-5.2$\,μ\mathrm{m}$ JWST NIRISS/SOSS+NIRSpec/G395H transmission spectrum of the 2.2$\,R_\oplus$ TOI-270d ($4.78\,M_\oplus$, $T_\mathrm{eq}$=350-380 K), delivering unprecedented sensitivity for atmospheric characterization in the sub-Neptune regime. We detect five vibrational bands of CH$_4$ at 1.15, 1.4, 1.7, 2.3, and 3.3$\,μ$m (9.4$σ$), the signature of CO$_2$ at 4.3$\,μ$m (4.8$σ$), water vapor (2.5$σ$), and potential signatures of SO$_2$ at 4.0$\,μ\mathrm{m}$ and CS$_2$ at 4.6$\,μ\mathrm{m}$. Intriguingly, we find an overall highly metal-rich atmosphere, with a mean molecular weight of $5.47_{-1.14}^{+1.25}$. We infer an atmospheric metal mass fraction of $58_{-12}^{+8}\%$ and a C/O of $0.47_{-0.19}^{+0.16}$, indicating that approximately half the mass of the outer envelope is in high-molecular-weight volatiles (H$_2$O, CH$_4$, CO, CO$_2$) rather than H$_2$/He. We introduce a sub-Neptune classification scheme and identify TOI-270d as a "miscible-envelope sub-Neptune" in which H$_2$/He is well-mixed with the high-molecular-weight volatiles in a miscible supercritical metal-rich envelope. For a fully miscible envelope, we conclude that TOI-270d's interior is $90_{-4}^{+3}\,$wt$\,\%$ rock/iron, indicating that it formed as a rocky planet that accreted a few wt % of H$_2$/He, with the overall envelope metal content explained by magma-ocean/envelope reactions without the need for significant ice accretion. TOI-270d may well be an archetype of the overall population of sub-Neptunes.
△ Less
Submitted 5 March, 2024;
originally announced March 2024.
-
Quenching-driven equatorial depletion and limb asymmetries in hot Jupiter atmospheres: WASP-96b example
Authors:
Maria Zamyatina,
Duncan A. Christie,
Eric Hébrard,
Nathan J. Mayne,
Michael Radica,
Jake Taylor,
Harry Baskett,
Ben Moore,
Craig Lils,
Denis Sergeev,
Eva-Maria Ahrer,
James Manners,
Krisztian Kohary,
Adina D. Feinstein
Abstract:
Transport-induced quenching in hot Jupiter atmospheres is a process that determines the boundary between the part of the atmosphere at chemical equilibrium and the part of the atmosphere at thermochemical (but not photothermochemical) disequilibrium. The location of this boundary, the quench level, depends on the interplay between the dynamical and chemical timescales in the atmosphere, with quenc…
▽ More
Transport-induced quenching in hot Jupiter atmospheres is a process that determines the boundary between the part of the atmosphere at chemical equilibrium and the part of the atmosphere at thermochemical (but not photothermochemical) disequilibrium. The location of this boundary, the quench level, depends on the interplay between the dynamical and chemical timescales in the atmosphere, with quenching occurring when these timescales are equal. We explore the sensitivity of the quench level position to an increase in the planet's atmospheric metallicity using aerosol-free 3D GCM simulations of a hot Jupiter WASP-96b. We find that the temperature increase at pressures of $\sim$$10^{4}-10^{7}$ Pa that occurs when metallicity is increased could shift the position of the quench level to pressures dominated by the jet, and cause an equatorial depletion of $CH_4$, $NH_3$ and $HCN$. We discuss how such a depletion affects the planet's transmission spectrum, and how the analysis of the evening-morning limb asymmetries, especially within $\sim3-5 μm$, could help distinguish atmospheres of different metallicities that are at chemical equilibrium from those with the upper layers at thermochemical disequilibrium.
△ Less
Submitted 22 February, 2024;
originally announced February 2024.
-
GMC Collisions As Triggers of Star Formation. VIII. The Core Mass Function
Authors:
Chia-Jung Hsu,
Jonathan C. Tan,
Duncan Christie,
Yu Cheng,
Theo J. O'Neill
Abstract:
Compression in giant molecular cloud (GMC) collisions is a promising mechanism to trigger formation of massive star clusters and OB associations. We simulate colliding and non-colliding magnetised GMCs and examine the properties of prestellar cores, selected from projected mass surface density maps, including after synthetic {\it ALMA} observations. We then examine core properties, including mass,…
▽ More
Compression in giant molecular cloud (GMC) collisions is a promising mechanism to trigger formation of massive star clusters and OB associations. We simulate colliding and non-colliding magnetised GMCs and examine the properties of prestellar cores, selected from projected mass surface density maps, including after synthetic {\it ALMA} observations. We then examine core properties, including mass, size, density, velocity, velocity dispersion, temperature and magnetic field strength. After four Myr, $\sim1,000$ cores have formed in the GMC collision and the high-mass end of the core mass function (CMF) can be fit by a power law $dN/d{\rm{log}}M\propto{M}^{-α}$ with $α\simeq0.7$, i.e., relatively top-heavy compared to a Salpeter mass function. Depending on how cores are identified, a break in the power law can appear around a few $\times10\:M_\odot$. The non-colliding GMCs form fewer cores with a CMF with $α\simeq0.8$ to 1.2, i.e., closer to the Salpeter index. We compare the properties of these CMFs to those of several observed samples of cores. Considering other properties, cores formed from colliding clouds are typically warmer, have more disturbed internal kinematics and are more likely to be gravitational unbound, than cores formed from non-colliding GMCs. The dynamical state of the protocluster of cores formed in the GMC-GMC collision is intrinsically subvirial, but can appear to be supervirial if the total mass measurement is affected by observations that miss mass on large scales or at low densities.
△ Less
Submitted 3 March, 2023; v1 submitted 25 January, 2023;
originally announced January 2023.
-
A JWST NIRSpec Phase Curve for WASP-121b: Dayside Emission Strongest Eastward of the Substellar Point and Nightside Conditions Conducive to Cloud Formation
Authors:
Thomas Mikal-Evans,
David K. Sing,
Jiayin Dong,
Daniel Foreman-Mackey,
Tiffany Kataria,
Joanna K. Barstow,
Jayesh M. Goyal,
Nikole K. Lewis,
Joshua D. Lothringer,
Nathan J. Mayne,
Hannah R. Wakeford,
Duncan A. Christie,
Zafar Rustamkulov
Abstract:
We present the first exoplanet phase curve measurement made with the JWST NIRSpec instrument, highlighting the exceptional stability of this newly-commissioned observatory for exoplanet climate studies. The target, WASP-121b, is an ultrahot Jupiter with an orbital period of 30.6 hr. We analyze two broadband light curves generated for the NRS1 and NRS2 detectors, covering wavelength ranges of 2.70-…
▽ More
We present the first exoplanet phase curve measurement made with the JWST NIRSpec instrument, highlighting the exceptional stability of this newly-commissioned observatory for exoplanet climate studies. The target, WASP-121b, is an ultrahot Jupiter with an orbital period of 30.6 hr. We analyze two broadband light curves generated for the NRS1 and NRS2 detectors, covering wavelength ranges of 2.70-3.72 micron and 3.82-5.15 micron, respectively. Both light curves exhibit minimal systematics, with approximately linear drifts in the baseline flux level of 30 ppm/hr (NRS1) and 10 ppm/hr (NRS2). Assuming a simple brightness map for the planet described by a low-order spherical harmonic dipole, our light curve fits suggest that the phase curve peaks coincide with orbital phases $3.36 \pm 0.11$ deg (NRS1) and $2.66 \pm 0.12$ deg (NRS2) prior to mid-eclipse. This is consistent with the strongest dayside emission emanating from eastward of the substellar point. We measure planet-to-star emission ratios of $3,924 \pm 7$ ppm (NRS1) and $4,924 \pm 9$ ppm (NRS2) for the dayside hemisphere, and $136 \pm 8$ ppm (NRS1) and $630 \pm 10$ ppm (NRS2) for the nightside hemisphere. The latter nightside emission ratios translate to planetary brightness temperatures of $926 \pm 12$ K (NRS1) and $1,122 \pm 10$ K (NRS2), which are low enough for a wide range of refractory condensates to form, including enstatite and forsterite. A nightside cloud deck may be blocking emission from deeper, hotter layers of the atmosphere, potentially helping to explain why cloud-free 3D general circulation model simulations systematically over-predict the nightside emission for WASP-121b.
△ Less
Submitted 16 February, 2023; v1 submitted 9 January, 2023;
originally announced January 2023.
-
Observability of signatures of transport-induced chemistry in clear atmospheres of hot gas giant exoplanets
Authors:
Maria Zamyatina,
Éric Hébrard,
Benjamin Drummond,
Nathan J. Mayne,
James Manners,
Duncan A. Christie,
Pascal Tremblin,
David K. Sing,
Krisztian Kohary
Abstract:
Transport-induced quenching, i.e., the homogenisation of chemical abundances by atmospheric advection, is thought to occur in the atmospheres of hot gas giant exoplanets. While some numerical modelling of this process exists, the three-dimensional nature of transport-induced chemistry is underexplored. Here we present results of 3D cloud- and haze-free simulations of the atmospheres of HAT-P-11b,…
▽ More
Transport-induced quenching, i.e., the homogenisation of chemical abundances by atmospheric advection, is thought to occur in the atmospheres of hot gas giant exoplanets. While some numerical modelling of this process exists, the three-dimensional nature of transport-induced chemistry is underexplored. Here we present results of 3D cloud- and haze-free simulations of the atmospheres of HAT-P-11b, HD 189733b, HD 209458b, and WASP-17b including coupled hydrodynamics, radiative transfer and chemistry. Our simulations were performed with two chemical schemes: a chemical kinetics scheme, which is capable of capturing transport-induced quenching, and a simpler, more widely used chemical equilibrium scheme. We find that transport-induced quenching is predicted to occur in atmospheres of all planets in our sample; however, the extent to which it affects their synthetic spectra and phase curves varies from planet to planet. This implies that there is a "sweet spot" for the observability of signatures of transport-induced quenching, which is controlled by the interplay between the dynamics and chemistry.
△ Less
Submitted 21 November, 2022; v1 submitted 16 November, 2022;
originally announced November 2022.
-
CAMEMBERT: A Mini-Neptunes GCM Intercomparison, Protocol Version 1.0. A CUISINES Model Intercomparison Project
Authors:
Duncan A. Christie,
Elspeth K. H. Lee,
Hamish Innes,
Pascal A. Noti,
Benjamin Charnay,
Thomas J. Fauchez,
Nathan J. Mayne,
Russell Deitrick,
Feng Ding,
Jennifer J. Greco,
Mark Hammond,
Isaac Malsky,
Avi Mandell,
Emily Rauscher,
Michael T. Roman,
Denis E. Sergeev,
Linda Sohl,
Maria E. Steinrueck,
Martin Turbet,
Eric T. Wolf,
Maria Zamyatina,
Ludmila Carone
Abstract:
With an increased focus on the observing and modelling of mini-Neptunes, there comes a need to better understand the tools we use to model their atmospheres. In this paper, we present the protocol for the CAMEMBERT (Comparing Atmospheric Models of Extrasolar Mini-neptunes Building and Envisioning Retrievals and Transits) project, an intercomparison of general circulation models (GCMs) used by the…
▽ More
With an increased focus on the observing and modelling of mini-Neptunes, there comes a need to better understand the tools we use to model their atmospheres. In this paper, we present the protocol for the CAMEMBERT (Comparing Atmospheric Models of Extrasolar Mini-neptunes Building and Envisioning Retrievals and Transits) project, an intercomparison of general circulation models (GCMs) used by the exoplanetary science community to simulate the atmospheres of mini-Neptunes. We focus on two targets well studied both observationally and theoretically with planned JWST Cycle 1 observations: the warm GJ~1214b and the cooler K2-18b. For each target, we consider a temperature-forced case, a clear sky dual-grey radiative transfer case, and a clear sky multi band radiative transfer case, covering a range of complexities and configurations where we know differences exist between GCMs in the literature. This paper presents all the details necessary to participate in the intercomparison, with the intention of presenting the results in future papers. Currently, there are eight GCMs participating (ExoCAM, Exo-FMS, FMS PCM, Generic PCM, MITgcm, RM-GCM, THOR, and the UM), and membership in the project remains open. Those interested in participating are invited to contact the authors.
△ Less
Submitted 8 November, 2022;
originally announced November 2022.
-
The Impact of Phase Equilibrium Cloud Models on GCM Simulations of GJ~1214b
Authors:
D. A. Christie,
N. J. Mayne,
R. M. Gillard,
J. Manners,
E. Hébrard,
S. Lines,
K. Kohary
Abstract:
We investigate the impact of clouds on the atmosphere of GJ~1214b using the radiatively-coupled, phase-equilibrium cloud model {\sc EddySed} coupled to the {\sc Unified Model} general circulation model. We find that, consistent with previous investigations, high metallicity ($100\times$ solar) and clouds with large vertical extents (a sedimentation factor of $f_\mathrm{sed} = 0.1$) are required to…
▽ More
We investigate the impact of clouds on the atmosphere of GJ~1214b using the radiatively-coupled, phase-equilibrium cloud model {\sc EddySed} coupled to the {\sc Unified Model} general circulation model. We find that, consistent with previous investigations, high metallicity ($100\times$ solar) and clouds with large vertical extents (a sedimentation factor of $f_\mathrm{sed} = 0.1$) are required to best match the observations, although metallicities even higher than those investigated here may be required to improve agreement further. We additionally find that in our case which best matches the observations ($f_\mathrm{sed}=0.1$), the velocity structures change relative to the clear sky case with the formation of a superrotating jet being suppressed, although further investigation is required to understand the cause of the suppression. The increase in cloud extent with $f_\mathrm{sed}$ results in a cooler planet due to a higher albedo, causing the atmosphere to contract. This also results in a reduced day-night contrast seen in the phase curves, although the introduction of cloud still results in a reduction of the phase offset. We additionally investigate the impact the the {\sc Unified Model}'s pseudo-spherical irradiation scheme on the calculation of heating rates, finding that the introduction of nightside shortwave heating results in slower mid-latitude jets compared to the plane parallel irradiation scheme used in previous works. We also consider the impact of a gamma distribution, as opposed to a log-normal distribution, for the distribution of cloud particle radii and find the impact to be relatively minor.
△ Less
Submitted 25 September, 2022;
originally announced September 2022.
-
3D radiative-transfer for exoplanet atmospheres. gCMCRT: a GPU accelerated MCRT code
Authors:
Elspeth K. H. Lee,
Joost P. Wardenier,
Bibiana Prinoth,
Vivien Parmentier,
Simon L. Grimm,
Robin Baeyens,
Ludmila Carone,
Duncan Christie,
Russell Deitrick,
Daniel Kitzmann,
Nathan Mayne,
Michael Roman,
Brian Thorsbro
Abstract:
Radiative-transfer (RT) is a key component for investigating atmospheres of planetary bodies. With the 3D nature of exoplanet atmospheres being important in giving rise to their observable properties, accurate and fast 3D methods are required to be developed to meet future multi-dimensional and temporal data sets. We develop an open source GPU RT code, gCMCRT, a Monte Carlo RT forward model for ge…
▽ More
Radiative-transfer (RT) is a key component for investigating atmospheres of planetary bodies. With the 3D nature of exoplanet atmospheres being important in giving rise to their observable properties, accurate and fast 3D methods are required to be developed to meet future multi-dimensional and temporal data sets. We develop an open source GPU RT code, gCMCRT, a Monte Carlo RT forward model for general use in planetary atmosphere RT problems. We aim to automate the post-processing pipeline, starting from direct global circulation model (GCM) output to synthetic spectra. We develop albedo, emission and transmission spectra modes for 3D and 1D input structures. We include capability to use correlated-k and high-resolution opacity tables, the latter of which can be Doppler shifted inside the model. We post-process results from several GCM groups including ExoRad, SPARC/MITgcm THOR, UK Met Office UM, Exo-FMS and the Rauscher model. Users can therefore take advantage of desktop and HPC GPU computing solutions. gCMCRT is well suited for post-processing large GCM model grids produced by members of the community and for high-resolution 3D investigations.
△ Less
Submitted 18 March, 2022; v1 submitted 29 October, 2021;
originally announced October 2021.
-
The Impact of Mixing Treatments on Cloud Modelling in 3D Simulations of Hot Jupiters
Authors:
D. A. Christie,
N. J. Mayne,
S. Lines,
V. Parmentier,
J. Manners,
I. Boutle,
B. Drummond,
T. Mikal-Evans,
D. K. Sing,
K. Kohary
Abstract:
We present results of 3D hydrodynamical simulations of HD209458b including a coupled, radiatively-active cloud model ({\sc EddySed}). We investigate the role of the mixing by replacing the default convective treatment used in previous works with a more physically relevant mixing treatment ($K_{zz}$) based on global circulation. We find that uncertainty in the efficiency of sedimentation through th…
▽ More
We present results of 3D hydrodynamical simulations of HD209458b including a coupled, radiatively-active cloud model ({\sc EddySed}). We investigate the role of the mixing by replacing the default convective treatment used in previous works with a more physically relevant mixing treatment ($K_{zz}$) based on global circulation. We find that uncertainty in the efficiency of sedimentation through the sedimentation factor $f_\mathrm{sed}$ plays a larger role in shaping cloud thickness and its radiative feedback on the local gas temperatures -- e.g. hot spot shift and day-to-night side temperature gradient -- than the switch in mixing treatment. We demonstrate using our new mixing treatments that simulations with cloud scales which are a fraction of the pressure scale height improve agreement with the observed transmission spectra, the emission spectra, and the Spitzer 4.5 $\mathrm{μm}$ phase curve, although our models are still unable to reproduce the optical and UV transmission spectra. We also find that the inclusion of cloud increases the transit asymmetry in the optical between the east and west limbs, although the difference remains small ($\lesssim 1\%$).
△ Less
Submitted 12 July, 2021;
originally announced July 2021.
-
GMC Collisions as Triggers of Star Formation. VII. The Effect of Magnetic Field Strength on Star Formation
Authors:
Benjamin Wu,
Jonathan C. Tan,
Duncan Christie,
Fumitaka Nakamura
Abstract:
We investigate the formation of stars within giant molecular clouds (GMCs) evolving in environments of different global magnetic field strength and large-scale dynamics. Building upon a series of magnetohydrodynamic (MHD) simulations of non-colliding and colliding GMCs, we employ density- and magnetically-regulated star formation sub-grid models in clouds which range from moderately magnetically s…
▽ More
We investigate the formation of stars within giant molecular clouds (GMCs) evolving in environments of different global magnetic field strength and large-scale dynamics. Building upon a series of magnetohydrodynamic (MHD) simulations of non-colliding and colliding GMCs, we employ density- and magnetically-regulated star formation sub-grid models in clouds which range from moderately magnetically supercritical to near critical. We examine gas and star cluster morphologies, magnetic field strengths and relative orientations, pre-stellar core densities, temperatures, mass-to-flux ratios and velocities, star formation rates and efficiencies over time, spatial clustering of stars, and kinematics of the stars and natal gas. The large scale magnetic criticality of the region greatly affects the overall gas evolution and star formation properties. GMC collisions enhance star formation rates and efficiencies in magnetically supercritical conditions, but may actually inhibit them in the magnetically critical case. This may have implications for star formation in different Galactic environments such as the Galactic Center and the main Galactic disk.
△ Less
Submitted 28 February, 2020; v1 submitted 8 November, 2019;
originally announced November 2019.
-
A Model of the H$α$ and Na Transmission Spectrum of HD 189733b
Authors:
Chenliang Huang,
Phil Arras,
Duncan Christie,
Zhi-Yun Li
Abstract:
This paper presents a detailed hydrostatic model of the upper atmosphere of HD 189733b, with the goal of constraining its temperature, particle densities, and radiation field over the pressure range $10^{-4}-10\, μ\rm bar$, where the observed H$α$ transmission spectrum is produced. The atomic hydrogen level population is computed including both collisional and radiative transition rates. The Ly…
▽ More
This paper presents a detailed hydrostatic model of the upper atmosphere of HD 189733b, with the goal of constraining its temperature, particle densities, and radiation field over the pressure range $10^{-4}-10\, μ\rm bar$, where the observed H$α$ transmission spectrum is produced. The atomic hydrogen level population is computed including both collisional and radiative transition rates. The Ly$α$ resonant scattering is computed using a Monte-Carlo simulation. The model transmission spectra are in broad agreement with the data. Excitation of the H(2$\ell$) population is mainly by Ly$α$ radiative excitation due to the large Ly$α$ intensity. The density of H(2$\ell$) is nearly flat over two decades in pressure, and is optically thick to H$α$. Additional models computed for a range of the stellar Lyman continuum (LyC) flux suggest that the variability in H$α$ transit depth may be due to the variability in the stellar LyC. Since metal lines provide the dominant cooling of this part of the atmosphere, the atmosphere structure is sensitive to the density of species such as Mg and Na which may themselves be constrained by observations. Since the H$α$ and Na D lines have comparable absorption depths, we argue that the center of the Na D lines are also formed in the atomic layer where the H$α$ line is formed.
△ Less
Submitted 15 November, 2017;
originally announced November 2017.
-
GMC Collisions as Triggers of Star Formation. VI. Collision-Induced Turbulence
Authors:
Benjamin Wu,
Jonathan C. Tan,
Fumitaka Nakamura,
Duncan Christie,
Qi Li
Abstract:
We investigate collisions between giant molecular clouds (GMCs) as potential generators of their internal turbulence. Using magnetohydrodynamic (MHD) simulations of self-gravitating, magnetized, turbulent, GMCs, we compare kinematic and dynamic properties of dense gas structures formed when such clouds collide compared to those that form in non-colliding clouds as self-gravity overwhelms decaying…
▽ More
We investigate collisions between giant molecular clouds (GMCs) as potential generators of their internal turbulence. Using magnetohydrodynamic (MHD) simulations of self-gravitating, magnetized, turbulent, GMCs, we compare kinematic and dynamic properties of dense gas structures formed when such clouds collide compared to those that form in non-colliding clouds as self-gravity overwhelms decaying turbulence. We explore the nature of turbulence in these structures via distribution functions of density, velocity dispersions, virial parameters, and momentum injection. We find that the dense clumps formed from GMC collisions have higher effective Mach number, greater overall velocity dispersions, sustain near-virial equilibrium states for longer times, and are the conduit for injection of turbulent momentum into high density gas at high rates.
△ Less
Submitted 3 November, 2017;
originally announced November 2017.
-
GMC Collisions As Triggers of Star Formation. IV. The Role of Ambipolar Diffusion
Authors:
Duncan Christie,
Benjamin Wu,
Jonathan C. Tan
Abstract:
We investigate the role of ambipolar diffusion (AD) in collisions between magnetized giant molecular clouds (GMCs), which may be an important mechanism for triggering star cluster formation. Three dimensional simulations of GMC collisions are performed using a version of the Enzo magnetohydrodynamics code that has been extended to include AD. The resistivities are calculated using the 31-species c…
▽ More
We investigate the role of ambipolar diffusion (AD) in collisions between magnetized giant molecular clouds (GMCs), which may be an important mechanism for triggering star cluster formation. Three dimensional simulations of GMC collisions are performed using a version of the Enzo magnetohydrodynamics code that has been extended to include AD. The resistivities are calculated using the 31-species chemical model of Wu et al. (2015). We find that in the weak-field, $10\:{\rm μG}$ case, AD has only a modest effect on the dynamical evolution during the collision. However, for the stronger-field, $30\:{\rm μG}$ case involving near-critical clouds, AD results in formation of dense cores in regions where collapse is otherwise inhibited. The overall efficiency of formation of cores with $n_{\rm H}\geq10^{6}\:{\rm cm}^{-3}$ in these simulations is increases from about 0.2% to 2% once AD is included, comparable to observed values in star-forming GMCs. The gas around these cores typically has relatively slow infall at speeds that are a modest fraction of the free-fall speed.
△ Less
Submitted 19 September, 2017; v1 submitted 21 June, 2017;
originally announced June 2017.
-
The Interstellar Medium and Star Formation of Galactic Disks. I. ISM and GMC properties with Diffuse FUV and Cosmic Ray Backgrounds
Authors:
Qi Li,
Jonathan C. Tan,
Duncan Christie,
Thomas G. Bisbas,
Benjamin Wu
Abstract:
We present a series of adaptive mesh refinement (AMR) hydrodynamic simulations of flat rotation curve galactic gas disks with a detailed treatment of the interstellar medium (ISM) physics of the atomic to molecular phase transition under the influence of diffuse FUV radiation fields and cosmic ray backgrounds. We explore the effects of different FUV intensities, including a model with a radial gra…
▽ More
We present a series of adaptive mesh refinement (AMR) hydrodynamic simulations of flat rotation curve galactic gas disks with a detailed treatment of the interstellar medium (ISM) physics of the atomic to molecular phase transition under the influence of diffuse FUV radiation fields and cosmic ray backgrounds. We explore the effects of different FUV intensities, including a model with a radial gradient designed to mimic the Milky Way. The effects of cosmic rays, including radial gradients in their heating and ionization rates, are also explored. The final simulations in this series achieve $4\:$pc resolution across the $\sim20\:$kpc global disk diameter, with heating and cooling followed down to temperatures of $\sim10\:$K. The disks are evolved for $300\:$Myr, which is enough time for the ISM to achieve a quasi-statistical equilibrium. In particular, the mass fraction of molecular gas stabilizes by $\sim$200 Myr. Additional global ISM properties are analysed. Giant molecular clouds (GMCs) are also identified and the statistical properties of their populations examined. GMCs are tracked as the disks evolve. GMC collisions, which may be a means of triggering star cluster formation, are counted and the rates compared with analytic models. Relatively frequent GMC collision rates are seen in these simulations and their implications for understanding GMC properties, including the driving of internal turbulence, are discussed.
△ Less
Submitted 13 November, 2017; v1 submitted 12 June, 2017;
originally announced June 2017.
-
GMC Collisions as Triggers of Star Formation. III. Density and Magnetically Regulated Star Formation
Authors:
Benjamin Wu,
Jonathan C. Tan,
Duncan Christie,
Fumitaka Nakamura,
Sven Van Loo,
David Collins
Abstract:
We study giant molecular cloud (GMC) collisions and their ability to trigger star cluster formation. We further develop our three dimensional magnetized, turbulent, colliding GMC simulations by implementing star formation sub-grid models. Two such models are explored: (1) "Density-Regulated," i.e., fixed efficiency per free-fall time above a set density threshold; (2) "Magnetically-Regulated," i.e…
▽ More
We study giant molecular cloud (GMC) collisions and their ability to trigger star cluster formation. We further develop our three dimensional magnetized, turbulent, colliding GMC simulations by implementing star formation sub-grid models. Two such models are explored: (1) "Density-Regulated," i.e., fixed efficiency per free-fall time above a set density threshold; (2) "Magnetically-Regulated," i.e., fixed efficiency per free-fall time in regions that are magnetically supercritical. Variations of parameters associated with these models are also explored. In the non-colliding simulations, the overall level of star formation is sensitive to model parameter choices that relate to effective density thresholds. In the GMC collision simulations, the final star formation rates and efficiencies are relatively independent of these parameters. Between non-colliding and colliding cases, we compare the morphologies of the resulting star clusters, properties of star-forming gas, time evolution of the star formation rate (SFR), spatial clustering of the stars, and resulting kinematics of the stars in comparison to the natal gas. We find that typical collisions, by creating larger amounts of dense gas, trigger earlier and enhanced star formation, resulting in 10 times higher SFRs and efficiencies. The star clusters formed from GMC collisions show greater spatial sub-structure and more disturbed kinematics.
△ Less
Submitted 26 April, 2017; v1 submitted 26 February, 2017;
originally announced February 2017.
-
GMC Collisions as Triggers of Star Formation. II. 3D Turbulent, Magnetized Simulations
Authors:
Benjamin Wu,
Jonathan C. Tan,
Fumitaka Nakamura,
Sven Van Loo,
Duncan Christie,
David Collins
Abstract:
We investigate giant molecular cloud (GMCs) collisions and their ability to induce gravitational instability and thus star formation. This mechanism may be a major driver of star formation activity in galactic disks. We carry out a series of three dimensional, magnetohydrodynamics (MHD), adaptive mesh refinement (AMR) simulations to study how cloud collisions trigger formation of dense filaments a…
▽ More
We investigate giant molecular cloud (GMCs) collisions and their ability to induce gravitational instability and thus star formation. This mechanism may be a major driver of star formation activity in galactic disks. We carry out a series of three dimensional, magnetohydrodynamics (MHD), adaptive mesh refinement (AMR) simulations to study how cloud collisions trigger formation of dense filaments and clumps. Heating and cooling functions are implemented based on photo-dissociation region (PDR) models that span the atomic to molecular transition and can return detailed diagnostic information. The clouds are initialized with supersonic turbulence and a range of magnetic field strengths and orientations. Collisions at various velocities and impact parameters are investigated. Comparing and contrasting colliding and non-colliding cases, we characterize morphologies of dense gas, magnetic field structure, cloud kinematic signatures, and cloud dynamics. We present key observational diagnostics of cloud collisions, especially: relative orientations between magnetic fields and density structures, like filaments; $^{13}$CO($J$=2-1), $^{13}$CO($J$=3-2), and $^{12}$CO($J$=8-7) integrated intensity maps and spectra; and cloud virial parameters. We compare these results to observed Galactic clouds.
△ Less
Submitted 28 November, 2016; v1 submitted 3 June, 2016;
originally announced June 2016.
-
Investigation of the environment around close-in transiting exoplanets using CLOUDY
Authors:
Jake D. Turner,
Duncan Christie,
Phil Arras,
Robert E. Johnson,
Carl Schmidt
Abstract:
It has been suggested that hot stellar wind gas in a bow shock around an exoplanet is sufficiently opaque to absorb stellar photons and give rise to an observable transit depth at optical and UV wavelengths. In the first part of this paper, we use the CLOUDY plasma simulation code to model the absorption from X-ray to radio wavelengths by 1-D slabs of gas in coronal equilibrium with varying densit…
▽ More
It has been suggested that hot stellar wind gas in a bow shock around an exoplanet is sufficiently opaque to absorb stellar photons and give rise to an observable transit depth at optical and UV wavelengths. In the first part of this paper, we use the CLOUDY plasma simulation code to model the absorption from X-ray to radio wavelengths by 1-D slabs of gas in coronal equilibrium with varying densities ($10^{4}-10^{8} \, {\rm cm^{-3}}$) and temperatures ($2000-10^{6} \ {\rm K}$) illuminated by a solar spectrum. For slabs at coronal temperatures ($10^{6} \ {\rm K}$) and densities even orders of magnitude larger than expected for the compressed stellar wind ($10^{4}-10^{5} \, {\rm cm^{-3}}$), we find optical depths orders of magnitude too small ($> 3\times10^{-7}$) to explain the $\sim3\%$ UV transit depths seen with Hubble. Using this result and our modeling of slabs with lower temperatures ($2000-10^4 {\rm K}$), the conclusion is that the UV transits of WASP-12b and HD 189733b are likely due to atoms originating in the planet, as the stellar wind is too highly ionized. A corollary of this result is that transport of neutral atoms from the denser planetary atmosphere outward must be a primary consideration when constructing physical models. In the second part of this paper, additional calculations using CLOUDY are carried out to model a slab of planetary gas in radiative and thermal equilibrium with the stellar radiation field. Promising sources of opacity from the X-ray to radio wavelengths are discussed, some of which are not yet observed.
△ Less
Submitted 3 March, 2016;
originally announced March 2016.
-
Axisymmetric Simulations of Hot Jupiter-Stellar Wind Hydrodynamic Interaction
Authors:
Duncan Christie,
Phil Arras,
Zhi-Yun Li
Abstract:
Gas giant exoplanets orbiting at close distances to the parent star are subjected to large radiation and stellar wind fluxes. In this paper, hydrodynamic simulations of the planetary upper atmosphere and its interaction with the stellar wind are carried out to understand the possible flow regimes and how they affect the Lyman-alpha transmission spectrum. Following Tremblin and Chiang, charge excha…
▽ More
Gas giant exoplanets orbiting at close distances to the parent star are subjected to large radiation and stellar wind fluxes. In this paper, hydrodynamic simulations of the planetary upper atmosphere and its interaction with the stellar wind are carried out to understand the possible flow regimes and how they affect the Lyman-alpha transmission spectrum. Following Tremblin and Chiang, charge exchange reactions are included to explore the role of energetic atoms as compared to thermal particles. In order to understand the role of the tail as compared to the leading edge of the planetary gas, the simulations were carried out under axisymmetry, and photoionization and stellar wind electron impact ionization reactions were included to limit the extent of the neutrals away from the planet. By varying the planetary gas temperature, two regimes are found. At high temperature, a supersonic planetary wind is found, which is turned around by the stellar wind and forms a tail behind the planet. At lower temperatures, the planetary wind is shut off when the stellar wind penetrates inside where the sonic point would have been. In this regime mass is lost by viscous interaction at the boundary between planetary and stellar wind gases. Absorption by cold hydrogen atoms is large near the planetary surface, and decreases away from the planet as expected. The hot hydrogen absorption is in an annulus and typically dominated by the tail, at large impact parameter, rather than by the thin leading edge of the mixing layer near the substellar point.
△ Less
Submitted 20 January, 2016;
originally announced January 2016.
-
Hα Absorption in Transiting Exoplanet Atmospheres
Authors:
Duncan Christie,
Phil Arras,
Zhi-Yun Li
Abstract:
Absorption of stellar Hα by the upper atmosphere of the planet HD189733b has recently been detected by Jensen et al. Motivated by this observation, we have developed a model for atomic hydrogen in the n=2 state and compared the resulting Hα line profile to the observations. The model atmosphere is in hydrostatic balance, as well as thermal and photoionization equilibrium. Collisional and radiative…
▽ More
Absorption of stellar Hα by the upper atmosphere of the planet HD189733b has recently been detected by Jensen et al. Motivated by this observation, we have developed a model for atomic hydrogen in the n=2 state and compared the resulting Hα line profile to the observations. The model atmosphere is in hydrostatic balance, as well as thermal and photoionization equilibrium. Collisional and radiative transitions are included in the determination of the n=2 state level population. We find that Hα absorption is dominated by an optical depth τ ~ 1 shell, composed of hydrogen in the metastable 2s state that is located below the hydrogen ionization layer. The number density of the 2s state within the shell is found to vary slowly with radius, while that of the 1s state falls rapidly. Thus while the Lyαabsorption, for a certain wavelength, occurs inside a relatively well defined impact parameter, the contribution to Hα absorption is roughly uniform over the entire atomic hydrogen layer. The model can approximately reproduce the observed Lyα and Hα integrated transit depths for HD189733b by using an ionization rate enhanced over that expected for the star by an order of magnitude. For HD 209458b, we are unable to explain the asymmetric Hα line profile observed by Jensen et al., as the model produces a symmetric line profile with transit depth comparable to that of HD 189733b. In an appendix, we study the effect of the stellar Lyα absorption on the net cooling rate.
△ Less
Submitted 16 July, 2013;
originally announced July 2013.
-
Do Lognormal Column-Density Distributions in Molecular Clouds Imply Supersonic Turbulence?
Authors:
K. Tassis,
D. A. Christie,
A. Urban,
J. L. Pineda,
T. Ch. Mouschovias,
H. W. Yorke,
H. Martel
Abstract:
Recent observations of column densities in molecular clouds find lognormal distributions with power-law high-density tails. These results are often interpreted as indications that supersonic turbulence dominates the dynamics of the observed clouds. We calculate and present the column-density distributions of three clouds, modeled with very different techniques, none of which is dominated by supers…
▽ More
Recent observations of column densities in molecular clouds find lognormal distributions with power-law high-density tails. These results are often interpreted as indications that supersonic turbulence dominates the dynamics of the observed clouds. We calculate and present the column-density distributions of three clouds, modeled with very different techniques, none of which is dominated by supersonic turbulence. The first star-forming cloud is simulated using smoothed particle hydrodynamics (SPH); in this case gravity, opposed only by thermal-pressure forces, drives the evolution. The second cloud is magnetically subcritical with subsonic turbulence, simulated using nonideal MHD; in this case the evolution is due to gravitationally-driven ambipolar diffusion. The third cloud is isothermal, self-gravitating, and has a smooth density distribution analytically approximated with a uniform inner region and an r^-2 profile at larger radii. We show that in all three cases the column-density distributions are lognormal. Power-law tails develop only at late times (or, in the case of the smooth analytic profile, for strongly centrally concentrated configurations), when gravity dominates all opposing forces. It therefore follows that lognormal column-density distributions are generic features of diverse model clouds, and should not be interpreted as being a consequence of supersonic turbulence.
△ Less
Submitted 14 June, 2010;
originally announced June 2010.
-
Formation of Interstellar Clouds: Parker Instability with Phase Transitions
Authors:
Telemachos Ch. Mouschovias,
Matthew W. Kunz,
Duncan A. Christie
Abstract:
We follow numerically the nonlinear evolution of the Parker instability in the presence of phase transitions from a warm to a cold HI interstellar medium in two spatial dimensions. The nonlinear evolution of the system favors modes that allow the magnetic field lines to cross the galactic plane. Cold HI clouds form with typical masses ~= 10^5 M_sun, mean densities ~= 20 cm^-3, mean magnetic fiel…
▽ More
We follow numerically the nonlinear evolution of the Parker instability in the presence of phase transitions from a warm to a cold HI interstellar medium in two spatial dimensions. The nonlinear evolution of the system favors modes that allow the magnetic field lines to cross the galactic plane. Cold HI clouds form with typical masses ~= 10^5 M_sun, mean densities ~= 20 cm^-3, mean magnetic field strengths ~= 4.3 muG (rms field strengths ~= 6.4 muG), mass-to-flux ratios ~= 0.1 - 0.3 relative to critical, temperatures ~= 50 K, (two-dimensional) turbulent velocity dispersions ~= 1.6 km s^-1, and separations ~= 500 pc, in agreement with observations. The maximum density and magnetic field strength are ~= 10^3 cm^-3 and ~= 20 muG, respectively. Approximately 60% of all HI mass is in the warm neutral medium. The cold neutral medium is arranged into sheet-like structures both perpendicular and parallel to the galactic plane, but it is also found almost everywhere in the galactic plane, with the density being highest in valleys of the magnetic field lines. `Cloudlets' also form whose physical properties are in quantitative agreement with those observed for such objects by Heiles (1967). The nonlinear phase of the evolution takes ~< 30 Myr, so that, if the instability is triggered by a nonlinear perturbation such as a spiral density shock wave, interstellar clouds can form within a time suggested by observations.
△ Less
Submitted 7 January, 2009;
originally announced January 2009.