-
A formation pathway for terrestrial planets with moderate water content involving atmospheric-volatile recycling
Authors:
Jonas Müller,
Bertram Bitsch,
Aaron David Schneider
Abstract:
Of the many recently discovered terrestrial exoplanets, some are expected to harbor moderate water mass fractions of a few percent. The formation pathways that can produce planets with these water mass fractions are not fully understood. Here, we use the code chemcomp, which consists of a semi-analytical 1D protoplanetary disk model harboring a migrating and accreting planet, to model the growth a…
▽ More
Of the many recently discovered terrestrial exoplanets, some are expected to harbor moderate water mass fractions of a few percent. The formation pathways that can produce planets with these water mass fractions are not fully understood. Here, we use the code chemcomp, which consists of a semi-analytical 1D protoplanetary disk model harboring a migrating and accreting planet, to model the growth and composition of planets with moderate water mass fractions by pebble accretion in a protoplanetary disk around a TRAPPIST-1 analog star. This star is accompanied by seven terrestrial planets, of which the outer four planets likely contain water mass fractions of between 1\% and 10\%. We adopt a published model that considers the evaporation of pebbles in the planetary envelope, from where recycling flows can transport the volatile vapor back into the disk. We find that with this model, the planetary water content depends on the influx rate of pebbles onto the planet. A decreasing pebble influx with time reduces the envelope temperature and consequently allows the formation of planets with moderate water mass fractions as inferred for the outer TRAPPIST-1 planets for a number of different simulation configurations. This is further evidence that the recycling of vapor is an important component of planet formation needed to explain the vast and diverse population of exoplanets.
△ Less
Submitted 30 June, 2024; v1 submitted 13 June, 2024;
originally announced June 2024.
-
Examining the detectability of ringing on highly eccentric exoplanets
Authors:
Mathijs Vanrespaille,
Robin Baeyens,
Aaron David Schneider,
Ludmila Carone,
Leen Decin
Abstract:
Eccentric exoplanets offer an opportunity to study the response of an atmosphere to changing thermal forcing and the robustness of the super-rotating equatorial jet seen on tidally locked hot Jupiters. However, the atmospheric dynamics on eccentric planets strongly depend on the planetary rotation period, which is difficult to constrain observationally. The ringing phenomenon, whereby the observed…
▽ More
Eccentric exoplanets offer an opportunity to study the response of an atmosphere to changing thermal forcing and the robustness of the super-rotating equatorial jet seen on tidally locked hot Jupiters. However, the atmospheric dynamics on eccentric planets strongly depend on the planetary rotation period, which is difficult to constrain observationally. The ringing phenomenon, whereby the observed emission increases and decreases after the periastron passage as the flash-heated hemisphere rotates into and out of view, can provide a tight constraint on rotation. We studied five highly eccentric transiting exoplanets HAT-P-2 b, HD 80606 b, TOI-3362 b, TOI-4127 b and HD 17156 b to find which displays strong ringing signals that are sufficiently strong for the James Webb Space Telescope (JWST) to detect. We implemented the treatment of eccentricity and non-synchronous rotation in the non-grey climate model expeRT/MITgcm and generated synthetic light curves. We find four detectable ringing peaks on HD 80606 b and some undetectable ringing on TOI-4127 b and HD 17156 b. The lack of clouds, photo-chemistry and obliquity in our models may have led us to overestimate the amplitude of the ringing however. The strength of the ringing signal is mostly determined by the eccentricity, planetary rotation period, planet-to-star radius ratio and apparent magnitude of the system. We searched for more exoplanets that could show ringing but found no candidates as promising as HD 80606 b. We recommend prioritising HD 80606 b as a target for ringing with JWST. A baseline of five days after the periastron passage would capture three ringing peaks, which is sufficient to tightly constrain the planetary rotation period. An extension to seven days would add a fourth peak, which would allow us to verify the rotation period.
△ Less
Submitted 26 March, 2024; v1 submitted 1 March, 2024;
originally announced March 2024.
-
chemcomp: Modeling the chemical composition of planets formed in protoplanetary disks
Authors:
Aaron David Schneider,
Bertram Bitsch
Abstract:
Future observations of exoplanets will hopefully reveal detailed constraints on planetary compositions. Recently, we have developed and introduced chemcomp (Schneider & Bitsch 2021a), which simulates the formation of planets in viscously evolving protoplanetary disks by the accretion of pebbles and gas. The chemical composition of planetary building blocks (pebbles and gas) is traced by including…
▽ More
Future observations of exoplanets will hopefully reveal detailed constraints on planetary compositions. Recently, we have developed and introduced chemcomp (Schneider & Bitsch 2021a), which simulates the formation of planets in viscously evolving protoplanetary disks by the accretion of pebbles and gas. The chemical composition of planetary building blocks (pebbles and gas) is traced by including a physical approach of the evaporation and condensation of volatiles at evaporation lines. We have now open-sourced the chemcomp code to enable comparisons between planet formation models and observational constraints by the community. The code can be found at https://github.com/AaronDavidSchneider/chemcomp, is easy to use (using configuration files) and comes with a detailed documentation and examples.
△ Less
Submitted 22 November, 2023;
originally announced January 2024.
-
Nightside clouds and disequilibrium chemistry on the hot Jupiter WASP-43b
Authors:
Taylor J. Bell,
Nicolas Crouzet,
Patricio E. Cubillos,
Laura Kreidberg,
Anjali A. A. Piette,
Michael T. Roman,
Joanna K. Barstow,
Jasmina Blecic,
Ludmila Carone,
Louis-Philippe Coulombe,
Elsa Ducrot,
Mark Hammond,
João M. Mendonça,
Julianne I. Moses,
Vivien Parmentier,
Kevin B. Stevenson,
Lucas Teinturier,
Michael Zhang,
Natalie M. Batalha,
Jacob L. Bean,
Björn Benneke,
Benjamin Charnay,
Katy L. Chubb,
Brice-Olivier Demory,
Peter Gao
, et al. (58 additional authors not shown)
Abstract:
Hot Jupiters are among the best-studied exoplanets, but it is still poorly understood how their chemical composition and cloud properties vary with longitude. Theoretical models predict that clouds may condense on the nightside and that molecular abundances can be driven out of equilibrium by zonal winds. Here we report a phase-resolved emission spectrum of the hot Jupiter WASP-43b measured from 5…
▽ More
Hot Jupiters are among the best-studied exoplanets, but it is still poorly understood how their chemical composition and cloud properties vary with longitude. Theoretical models predict that clouds may condense on the nightside and that molecular abundances can be driven out of equilibrium by zonal winds. Here we report a phase-resolved emission spectrum of the hot Jupiter WASP-43b measured from 5-12 $μ$m with JWST's Mid-Infrared Instrument (MIRI). The spectra reveal a large day-night temperature contrast (with average brightness temperatures of 1524$\pm$35 and 863$\pm$23 Kelvin, respectively) and evidence for water absorption at all orbital phases. Comparisons with three-dimensional atmospheric models show that both the phase curve shape and emission spectra strongly suggest the presence of nightside clouds which become optically thick to thermal emission at pressures greater than ~100 mbar. The dayside is consistent with a cloudless atmosphere above the mid-infrared photosphere. Contrary to expectations from equilibrium chemistry but consistent with disequilibrium kinetics models, methane is not detected on the nightside (2$σ$ upper limit of 1-6 parts per million, depending on model assumptions).
△ Less
Submitted 23 January, 2024;
originally announced January 2024.
-
Harnessing machine learning for accurate treatment of overlapping opacity species in general circulation models
Authors:
Aaron David Schneider,
Paul Mollière,
Gilles Louppe,
Ludmila Carone,
Uffe Gråe Jørgensen,
Leen Decin,
Christiane Helling
Abstract:
To understand high precision observations of exoplanets and brown dwarfs, we need detailed and complex general circulation models (GCMs) that incorporate hydrodynamics, chemistry, and radiation. For this study, we specifically examined the coupling between chemistry and radiation in GCMs and compared different methods for the mixing of opacities of different chemical species in the correlated-k as…
▽ More
To understand high precision observations of exoplanets and brown dwarfs, we need detailed and complex general circulation models (GCMs) that incorporate hydrodynamics, chemistry, and radiation. For this study, we specifically examined the coupling between chemistry and radiation in GCMs and compared different methods for the mixing of opacities of different chemical species in the correlated-k assumption, when equilibrium chemistry cannot be assumed. We propose a fast machine learning method based on DeepSets (DS), which effectively combines individual correlated-k opacities (k-tables). We evaluated the DS method alongside other published methods such as adaptive equivalent extinction (AEE) and random overlap with rebinning and resorting (RORR). We integrated these mixing methods into our GCM (expeRT/MITgcm) and assessed their accuracy and performance for the example of the hot Jupiter HD~209458 b. Our findings indicate that the DS method is both accurate and efficient for GCM usage, whereas RORR is too slow. Additionally, we observed that the accuracy of AEE depends on its specific implementation and may introduce numerical issues in achieving radiative transfer solution convergence. We then applied the DS mixing method in a simplified chemical disequilibrium situation, where we modeled the rainout of TiO and VO, and confirmed that the rainout of TiO and VO would hinder the formation of a stratosphere. To further expedite the development of consistent disequilibrium chemistry calculations in GCMs, we provide documentation and code for coupling the DS mixing method with correlated-k radiative transfer solvers. The DS method has been extensively tested to be accurate enough for GCMs; however, other methods might be needed for accelerating atmospheric retrievals.
△ Less
Submitted 6 December, 2023; v1 submitted 1 November, 2023;
originally announced November 2023.
-
Constraining the formation history of the HAT-P-11 system by atmospheric abundances
Authors:
Lena Chatziastros,
Bertram Bitsch,
Aaron David Schneider
Abstract:
The chemical fingerprint of a planet reveals information about its formation history regarding when and where it formed. The water content of a planet can help to constrain its formation pathway: If the planet formed in the outer regions of the disk and migrated inward, it will be water-rich due to the accretion of water-ice-rich solids. Conversely, formation in the inner disk, where water-ice is…
▽ More
The chemical fingerprint of a planet reveals information about its formation history regarding when and where it formed. The water content of a planet can help to constrain its formation pathway: If the planet formed in the outer regions of the disk and migrated inward, it will be water-rich due to the accretion of water-ice-rich solids. Conversely, formation in the inner disk, where water-ice is not available, will result in a smaller atmospheric water content due to the limited accretion of water vapor. This process complicates with the presence of gap-opening giant planets. A gas giant exerts a pressure bump exterior to its orbit, preventing the further influx of pebbles into the inner system, resulting in a water-poor environment. The different formation scenarios can help to constrain the formation of the HAT-P-11 system, which contains an inner sub-Neptune with a mass of 23.4 $\mathrm{M_{\oplus}}$ and substellar water abundances ($X_\mathrm{H_2O} \approx 0.11$), as well as an outer giant planet orbiting exterior to the water-ice line. Our planet formation model encompasses planetary growth through pebble and gas accretion, along with a pebble drift and evaporation module that enables us to track the chemical composition of the disk and the planets. We find that the presence of the gas giant is necessary to block water-ice-rich material, resulting in a substellar water content for the inner sub-Neptune. If the giant planet forms too early, not enough solid material can enter the inner disk regions, preventing the efficient growth of the inner planet. This highlights the importance of the timing of giant planet formation to explain the inner system structure. Our simulations predict a roughly stellar C/O ratio with superstellar C/H and O/H for HAT-P-11b, providing constraints for future observations that are essential for gaining a more detailed understanding of its formation.
△ Less
Submitted 25 October, 2023; v1 submitted 19 October, 2023;
originally announced October 2023.
-
Photodissociation and induced chemical asymmetries on ultra-hot gas giants. A case study of HCN on WASP-76 b
Authors:
Robin Baeyens,
Jean-Michel Désert,
Annemieke Petrignani,
Ludmila Carone,
Aaron David Schneider
Abstract:
Recent observations have resulted in the detection of chemical gradients on ultra-hot gas giants. Notwithstanding their high temperature, chemical reactions in ultra-hot atmospheres may occur in disequilibrium, due to vigorous day-night circulation and intense UV radiation from their stellar hosts. The goal of this work is to explore whether photochemistry is affecting the composition of ultra-hot…
▽ More
Recent observations have resulted in the detection of chemical gradients on ultra-hot gas giants. Notwithstanding their high temperature, chemical reactions in ultra-hot atmospheres may occur in disequilibrium, due to vigorous day-night circulation and intense UV radiation from their stellar hosts. The goal of this work is to explore whether photochemistry is affecting the composition of ultra-hot giant planets, and if it can introduce horizontal chemical gradients. In particular, we focus on hydrogen cyanide (HCN) on WASP-76 b, as it is a photochemically active molecule with a reported detection on only one side of this planet. We use a pseudo-2D chemical kinetics code to model the chemical composition of WASP-76 b along its equator. Our approach improves on chemical equilibrium models by computing vertical mixing, horizontal advection, and photochemistry. We find that production of HCN is initiated through thermal and photochemical dissociation of CO and N$_2$ on the day side of WASP-76 b. The resulting radicals are subsequently transported to the night side via the equatorial jet stream, where they recombine into different molecules. This process results in an HCN gradient with a maximal abundance on the planet's morning limb. We verified that photochemical dissociation is a necessary condition for this mechanism, as thermal dissociation alone proves insufficient. Other species produced via night-side disequilibrium chemistry are SO$_2$ and S$_2$. Our model acts as a proof of concept for chemical gradients on ultra-hot exoplanets. We demonstrate that even ultra-hot planets can exhibit disequilibrium chemistry and recommend that future studies do not neglect photochemistry in their analyses of ultra-hot planets.
△ Less
Submitted 11 March, 2024; v1 submitted 1 September, 2023;
originally announced September 2023.
-
Evidence of Radius Inflation in Radiative GCM Models of WASP-76b due to the Advection of Potential Temperature
Authors:
Felix Sainsbury-Martinez,
Pascal Tremblin,
Aaron David Schneider,
Ludmila Carone,
Isabelle Baraffe,
Gilles Chabrier,
Christiane Helling,
Leen Decin,
Uffe Gråe Jørgensen
Abstract:
Understanding the discrepancy between the radii of observed hot Jupiters and standard 'radiative-convective' models remains a hotly debated topic in the exoplanet community. One mechanism which has been proposed to bridge this gap, and which has recently come under scrutiny, is the vertical advection of potential temperature from the irradiated outer atmosphere deep into the interior, heating the…
▽ More
Understanding the discrepancy between the radii of observed hot Jupiters and standard 'radiative-convective' models remains a hotly debated topic in the exoplanet community. One mechanism which has been proposed to bridge this gap, and which has recently come under scrutiny, is the vertical advection of potential temperature from the irradiated outer atmosphere deep into the interior, heating the deep, unirradiated, atmosphere, warming the internal adiabat, and resulting in radius inflation. Specifically, a recent study which explored the atmosphere of WASP-76b using a 3D, non-grey, GCM suggested that their models lacked radius inflation, and hence any vertical enthalpy advection. Here we perform additional analysis of these, and related, models, focusing on an explicit analysis of vertical enthalpy transport and the resulting heating of the deep atmosphere compared with 1D models. Our results indicate that, after any evolution linked with initialisation, all the WASP-76b models considered here exhibit significant vertical enthalpy transport, heating the deep atmosphere significantly when compared with standard 1D models. Furthermore, comparison of a long time-scale (and hence near steady-state) model with a Jupiter-like internal-structure model suggests not only strong radius-inflation, but also that the model radius, $1.98 \mathrm{R_{J}}$, may be comparable with observations ($1.83\pm0.06 \mathrm{R_{J}}$). We thus conclude that the vertical advection of potential temperature alone is enough to explain the radius inflation of WASP-76b, and potentially other irradiated gas giants, albeit with the proviso that the exact strength of the vertical advection remains sensitive to model parameters, such as the inclusion of deep atmospheric drag.
△ Less
Submitted 21 June, 2023;
originally announced June 2023.
-
Thermal Effects in Binary Neutron Star Mergers
Authors:
Jacob Fields,
Aviral Prakash,
Matteo Breschi,
David Radice,
Sebastiano Bernuzzi,
André da Silva Schneider
Abstract:
We study the impact of finite-temperature effects in numerical-relativity simulations of binary neutron star mergers with microphysical equations of state and neutrino transport in which we vary the effective nucleon masses in a controlled way. We find that, as the specific heat is increased, the merger remnants become colder and more compact due to the reduced thermal pressure support. Using a fu…
▽ More
We study the impact of finite-temperature effects in numerical-relativity simulations of binary neutron star mergers with microphysical equations of state and neutrino transport in which we vary the effective nucleon masses in a controlled way. We find that, as the specific heat is increased, the merger remnants become colder and more compact due to the reduced thermal pressure support. Using a full Bayesian analysis, we demonstrate that this effect will be measurable in the postmerger gravitational wave signal with next-generation observatories at signal-to-noise ratios of 15.
△ Less
Submitted 10 July, 2023; v1 submitted 22 February, 2023;
originally announced February 2023.
-
WASP-39b: exo-Saturn with patchy cloud composition, moderate metallicity, and underdepleted S/O
Authors:
Ludmila Carone,
David A. Lewis,
Dominic Samra,
Aaron D. Schneider,
Christiane Helling
Abstract:
WASP-39b is one of the first extrasolar giant gas planets that has been observed within the JWST ERS program. Fundamental properties that may enable the link to exoplanet formation differ amongst retrieval methods, for example metallicity and mineral ratios.
In this work, the formation of clouds in the atmosphere of WASP-39b is explored to investigate how inhomogeneous cloud properties (particle…
▽ More
WASP-39b is one of the first extrasolar giant gas planets that has been observed within the JWST ERS program. Fundamental properties that may enable the link to exoplanet formation differ amongst retrieval methods, for example metallicity and mineral ratios.
In this work, the formation of clouds in the atmosphere of WASP-39b is explored to investigate how inhomogeneous cloud properties (particle sizes, material composition, opacity) may be for this intermediately warm gaseous exoplanet. WASP-39b's atmosphere has a comparable day-night temperature median with sufficiently low temperatures that clouds may form globally. The presence of clouds on WASP-39b can explain observations without resorting to a high (> 100x solar) metallicity atmosphere for a reduced vertical mixing efficiency. The assessment of mineral ratios shows an under-depletion of S/O due to condensation compared to C/O, Mg/O, Si/O, Fe/O ratios. Vertical patchiness due to heterogeneous cloud composition challenges simple cloud models. An equal mixture of silicates and metal oxides is expected to characterise the cloud top. Further, optical properties of Fe and Mg silicates in the mid-infrared differ significantly which will impact the interpretation of JWST observations. We conclude that WASP-39b's atmosphere contains clouds and the underdepletion of S/O by atmospheric condensation processes suggest the use of sulphur gas species as a possible link to primordial element abundances. Over-simplified cloud models do not capture the complex nature of mixed-condensate clouds in exoplanet atmospheres. The clouds in the observable upper atmosphere of WASP-39b are a mixture of different silicates and metal oxides. The use of constant particles sizes and/or one-material cloud particles alone to interpret spectra may not be sufficient to capture the full complexity available through JWST observations.
△ Less
Submitted 20 January, 2023;
originally announced January 2023.
-
No evidence for radius inflation in hot Jupiters from vertical advection of heat
Authors:
Aaron David Schneider,
Ludmila Carone,
Leen Decin,
Uffe Gråe Jørgensen,
Christiane Helling
Abstract:
Understanding the radiative-dynamical coupling between upper photosphere and deeper atmosphere is a key in understanding the abnormal large radii of hot Jupiters. One needs very long integration times of 3D general circulation models (GCMs) with self consistent radiative transfer to achieve a better understanding of the feedback process between dynamics and radiation. We here present the longest 3…
▽ More
Understanding the radiative-dynamical coupling between upper photosphere and deeper atmosphere is a key in understanding the abnormal large radii of hot Jupiters. One needs very long integration times of 3D general circulation models (GCMs) with self consistent radiative transfer to achieve a better understanding of the feedback process between dynamics and radiation. We here present the longest 3D non-gray GCM study (86000 d) of an ultra hot Jupiter (WASP-76 b) published to this date that reached a final converged state. Furthermore, we present a method that can be used to accelerate the path towards temperature convergence in the deep atmospheric layers. We find that the final converged temperature profile is cold in the deep atmospheric layers, lacking any sign of vertical transport of potential temperature by large scale atmospheric motions. We thus conclude that the coupling between radiation and dynamics alone is not sufficient to explain the abnormal large radii of inflated hot gas giants.
△ Less
Submitted 4 October, 2022;
originally announced October 2022.
-
How drifting and evaporating pebbles shape giant planets III: The formation of WASP-77A b and $τ$ Boötis b
Authors:
Bertram Bitsch,
Aaron David Schneider,
Laura Kreidberg
Abstract:
Atmospheric abundances are thought to constrain the planet formation pathway, because different species evaporate at different temperatures leaving distinct signatures in the accreted atmosphere. The planetary C/O ratio is thought to constrain the planet formation pathway, because of the condensation sequence of H$_2$O, CO$_2$, CH$_4$, and CO, resulting in an increase of the gas phase C/O ratio wi…
▽ More
Atmospheric abundances are thought to constrain the planet formation pathway, because different species evaporate at different temperatures leaving distinct signatures in the accreted atmosphere. The planetary C/O ratio is thought to constrain the planet formation pathway, because of the condensation sequence of H$_2$O, CO$_2$, CH$_4$, and CO, resulting in an increase of the gas phase C/O ratio with increasing distance. Here we use a disc evolution model including pebble growth, drift and evaporation coupled with a planet formation model that includes pebble and gas accretion as well as planet migration to compute the atmospheric compositions of giant planets. We compare our results to the recent observations of the hot Jupiters WASP-77A b and $τ$ Boötis b, which feature sub-solar and super-solar C/H and O/H values, respectively. Our simulations reproduce these measurements and show that giants like WASP-77A b should start to form beyond the CO$_2$ evaporation front, while giants like $τ$ Boötis b should originate from beyond the H$_2$O line. Our model allows the formation of sub- and super-solar atmospheric compositions. However simulations without pebble evaporation can not reproduce the super-solar C/H and O/H ratios of $τ$ Boötis b's atmosphere without additional accretion of solids. We identify the $α$ viscosity parameter of the disc as a key ingredient, because the viscosity drives the inward motion of volatile enriched vapor, responsible for the accretion of gaseous C and O. Depending on the planet's migration history order-of-magnitude differences in atmospheric C/H and O/H are expected. Our simulations also predict super-solar N/H for $τ$ Boötis b and solar N/H for WASP-77A b. We conclude that pebble evaporation is a key ingredient to explain the variety of exoplanet atmospheres, because it can explain both, sub- and super-solar atmospheric abundances.
△ Less
Submitted 14 July, 2022; v1 submitted 13 July, 2022;
originally announced July 2022.
-
Exploring the deep atmospheres of HD 209458b and WASP-43b using a non-gray general circulation model
Authors:
Aaron David Schneider,
Ludmila Carone,
Leen Decin,
Uffe Gråe Jørgensen,
Paul Mollière,
Robin Baeyens,
Sven Kiefer,
Christiane Helling
Abstract:
Simulations with a 3D general circulation model (GCM) suggest that one potential driver behind the observed radius inflation in hot Jupiters may be the downward advection of energy from the highly irradiated photosphere into the deeper layers. Here, we compare dynamical heat transport within the non-inflated hot Jupiter WASP-43b and the canonical inflated hot Jupiter HD 209458b, with similar effec…
▽ More
Simulations with a 3D general circulation model (GCM) suggest that one potential driver behind the observed radius inflation in hot Jupiters may be the downward advection of energy from the highly irradiated photosphere into the deeper layers. Here, we compare dynamical heat transport within the non-inflated hot Jupiter WASP-43b and the canonical inflated hot Jupiter HD 209458b, with similar effective temperatures. We investigate to what extent the radiatively driven heating and cooling in the photosphere (at pressures smaller than 1 bar) influence the deeper temperature profile (at pressures between 1 to 700 bar). Our simulations with the new non-gray 3D radiation-hydrodynamical model expeRT/MITgcm show that the deep temperature profile of WASP-43b is associated with a relatively cold adiabat. The deep layers of HD 209458b, however, do not converge and remain nearly unchanged regardless of whether a cold or a hot initial state is used. Furthermore, we show that different flow structures in the deep atmospheric layers arise. There, we find that WASP-43b exhibits a deep equatorial jet, driven by the relatively fast tidally locked rotation of this planet (0.81 days), as compared to HD 209458b (3.47 days). However, by comparing simulations with different rotation periods, we find that the resulting flow structures only marginally influence the temperature evolution in the deep atmosphere, which is almost completely dominated by radiative heating and cooling. Furthermore, we find that the evolution of deeper layers can influence the 3D temperature structure in the photosphere of WASP-43b. Thus, dayside emission spectra of WASP-43b may shed more light onto the dynamical processes occurring at greater depths.
△ Less
Submitted 3 October, 2022; v1 submitted 18 February, 2022;
originally announced February 2022.
-
Numerical relativity simulations of prompt collapse mergers: threshold mass and phenomenological constraints on neutron star properties after GW170817
Authors:
Rahul Kashyap,
Abhishek Das,
David Radice,
Surendra Padamata,
Aviral Prakash,
Domenico Logoteta,
Albino Perego,
Daniel A. Godzieba,
Sebastiano Bernuzzi,
Ignazio Bombaci,
Farrukh J. Fattoyev,
Brendan T. Reed,
André da Silva Schneider
Abstract:
We determine the threshold mass for prompt (no bounce) black hole formation in equal-mass neutron star (NS) mergers using a new set of 227 numerical relativity simulations. We consider 23 phenomenological and microphysical finite temperature equations of state (EOS), including models with hyperons and first-order phase transitions to deconfined quarks. We confirm the existence of EOS-insensitive r…
▽ More
We determine the threshold mass for prompt (no bounce) black hole formation in equal-mass neutron star (NS) mergers using a new set of 227 numerical relativity simulations. We consider 23 phenomenological and microphysical finite temperature equations of state (EOS), including models with hyperons and first-order phase transitions to deconfined quarks. We confirm the existence of EOS-insensitive relations between the threshold mass, the binary tidal parameter at the threshold ($Λ_{th}$), the maximum mass of nonrotating NSs, and the radii of reference mass NSs. We correct the systematic errors in previously reported fitting coefficients that were obtained with approximate general-relativity simulations. We combine the EOS-insensitive relations, phenomenological constraints on NS properties and observational data from GW170817 to derive an improved lower limit on radii of maximum mass and 1.6 M$_\odot$ NS of 9.81 km and 10.90 km, respectively. We also constrain the radius and quadrupolar tidal deformability ($Λ$) of a 1.4 $M_\odot$ NS to be larger than 10.74 km and 172, respectively. We consider uncertainties in all independent parameters -- fitting coefficients as well as GW170817 masses while reporting the range of radii constraints. We introduce new methods to constrain the upper as well as lower limit of NS maximum mass using future BNS detections and their identification as prompt or delayed collapse. With future observations it will be possible to derive even tighter constraints on the properties of matter at and above nuclear density using the method proposed in this work.
△ Less
Submitted 8 April, 2022; v1 submitted 9 November, 2021;
originally announced November 2021.
-
How drifting and evaporating pebbles shape giant planets II: Volatiles and refractories in atmospheres
Authors:
Aaron David Schneider,
Bertram Bitsch
Abstract:
Upcoming studies of extrasolar gas giants will give precise insights into the composition of planetary atmospheres with the ultimate goal to link it to the formation history of the planet. Here, we investigate how drifting and evaporating pebbles that enrich the gas phase of the disk influence the chemical composition of growing and migrating gas giants. To achieve this goal, we perform semi analy…
▽ More
Upcoming studies of extrasolar gas giants will give precise insights into the composition of planetary atmospheres with the ultimate goal to link it to the formation history of the planet. Here, we investigate how drifting and evaporating pebbles that enrich the gas phase of the disk influence the chemical composition of growing and migrating gas giants. To achieve this goal, we perform semi analytical 1D models of protoplanetary disks including viscous evolution, pebble drift and evaporation to simulate the growth of planets from planetary embryos to Jupiter mass objects by the accretion of pebbles and gas while they migrate through the disk. The gas phase of the protoplanetary disk is enriched due to the evaporation of inward drifting pebbles crossing evaporation lines, leading to the accretion of large amounts of volatiles into the planetary atmosphere. As a consequence, gas accreting planets are enriched in volatiles (C, O, N) compared to refractories (e.g., Mg, Si, Fe) by up to a factor of 100, depending on the chemical species, its exact abundance and volatility as well as the disk's viscosity. A simplified model for the formation of Jupiter reveals that its nitrogen content can be explained by inward diffusing nitrogen rich vapor, implying that Jupiter does not need to form close to the N2 evaporation front as indicated by previous simulations. However, our model predicts a too low oxygen abundance for Jupiter, implying either Jupiter's migration across the water ice line or an additional accretion of solids into the atmosphere, which can also increase Jupiter's carbon abundance. The accretion of solids will increase the refractory to volatile ratio in planetary atmospheres substantially. We thus conclude that the volatile to refractory ratio can place a strong constraint on planet formation theories making it an important target for future observations.
△ Less
Submitted 15 September, 2021; v1 submitted 8 September, 2021;
originally announced September 2021.
-
Equation of State Dependence of Gravitational Waves in Core-Collapse Supernovae
Authors:
Oliver Eggenberger Andersen,
Shuai Zha,
André da Silva Schneider,
Aurore Betranhandy,
Sean M. Couch,
Evan P. O'Connor
Abstract:
Gravitational waves (GWs) provide unobscured insight into the birthplace of neutron stars (NSs) and black holes in core-collapse supernovae (CCSNe). The nuclear equation of state (EOS) describing these dense environments is yet uncertain, and variations in its prescription affect the proto-neutron star (PNS) and the post-bounce dynamics in CCSNe simulations, subsequently impacting the GW emission.…
▽ More
Gravitational waves (GWs) provide unobscured insight into the birthplace of neutron stars (NSs) and black holes in core-collapse supernovae (CCSNe). The nuclear equation of state (EOS) describing these dense environments is yet uncertain, and variations in its prescription affect the proto-neutron star (PNS) and the post-bounce dynamics in CCSNe simulations, subsequently impacting the GW emission. We perform axisymmetric simulations of CCSNe with Skyrme-type EOSs to study how the GW signal and PNS convection zone are impacted by two experimentally accessible EOS parameters, (1) the effective mass of nucleons, $m^\star$, which is crucial in setting the thermal dependence of the EOS, and (2) the isoscalar incompressibility modulus, $K_{\rm{sat}}$. While $K_{\rm{sat}}$ shows little impact, the peak frequency of the GWs has a strong effective mass dependence due to faster contraction of the PNS for higher values of $m^\star$ owing to a decreased thermal pressure. These more compact PNSs also exhibit more neutrino heating which drives earlier explosions and correlates with the GW amplitude via accretion plumes striking the PNS, exciting the oscillations. We investigate the spatial origin of the GWs and show the agreement between a frequency-radial distribution of the GW emission and a perturbation analysis. We do not rule out overshoot from below via PNS convection as another moderately strong excitation mechanism in our simulations. We also study the combined effect of effective mass and rotation. In all our simulations we find evidence for a power gap near $\sim$1250 Hz, we investigate its origin and report its EOS dependence.
△ Less
Submitted 17 June, 2021;
originally announced June 2021.
-
How drifting and evaporating pebbles shape giant planets I: Heavy element content and atmospheric C/O
Authors:
Aaron David Schneider,
Bertram Bitsch
Abstract:
Recent observations of extrasolar gas giants suggest super-stellar C/O ratios in planetary atmospheres, while interior models of observed extrasolar giant planets additionally suggest high heavy element contents. Furthermore, recent observations of protoplanetary disks revealed super-solar C/H ratios, which are explained by inward drifting and evaporating pebbles, enhancing the volatile content of…
▽ More
Recent observations of extrasolar gas giants suggest super-stellar C/O ratios in planetary atmospheres, while interior models of observed extrasolar giant planets additionally suggest high heavy element contents. Furthermore, recent observations of protoplanetary disks revealed super-solar C/H ratios, which are explained by inward drifting and evaporating pebbles, enhancing the volatile content of the disk. We investigate how the inward drift and evaporation of volatile rich pebbles influences the atmospheric C/O ratio and heavy element content of giant planets growing by pebble and gas accretion. To achieve this goal, we perform semi analytical 1D models of protoplanetary disks including the treatment of viscous evolution and heating, pebble drift and simple chemistry to simulate the growth of planets from planetary embryos to Jupiter mass objects by accretion of pebbles and gas while they migrate through the disk. Our simulations show that the composition of the planetary gas atmosphere is dominated by the accretion of vapour, originating from inward drifting evaporating pebbles. This process allows the giant planets to harbour large heavy element contents. In addition, our model reveals that giant planets originating further away from the central star have a higher C/O ratio on average due to the evaporation of methane rich pebbles in the outer disk. These planets can then also harbour super-solar C/O ratios, in line with exoplanet observations. However, planets formed in the outer disk harbour a smaller heavy element content, due to a smaller vapour enrichment of the outer disk. Our model predicts that giant planets with low/large atmospheric C/O should harbour a large/low total heavy element content. We further conclude that the inclusion of pebble evaporation at evaporation lines is a key ingredient to determine the heavy element content and composition of giant planets.
△ Less
Submitted 9 September, 2021; v1 submitted 27 May, 2021;
originally announced May 2021.
-
Dry or water world? How the water contents of inner sub-Neptunes constrain giant planet formation and the location of the water ice line
Authors:
Bertram Bitsch,
Sean N. Raymond,
Lars A. Buchhave,
Aaron Bello-Arufe,
Alexander D. Rathcke,
Aaron David Schneider
Abstract:
In the pebble accretion scenario, the pebbles that form planets drift inward from the outer disk regions, carrying water ice with them. At the water ice line, the water ice on the inward drifting pebbles evaporates and is released into the gas phase, resulting in water-rich gas and dry pebbles that move into the inner disk regions. Large planetary cores can block the inward drifting pebbles by for…
▽ More
In the pebble accretion scenario, the pebbles that form planets drift inward from the outer disk regions, carrying water ice with them. At the water ice line, the water ice on the inward drifting pebbles evaporates and is released into the gas phase, resulting in water-rich gas and dry pebbles that move into the inner disk regions. Large planetary cores can block the inward drifting pebbles by forming a pressure bump outside their orbit in the protoplanetary disk. Depending on the relative position of a growing planetary core relative to the water ice line, water-rich pebbles might be blocked outside or inside the water ice line. Pebbles blocked outside the water ice line do not evaporate and thus do not release their water vapor into the gas phase, resulting in a dry inner disk, while pebbles blocked inside the water ice line release their water vapor into the gas phase, resulting in water vapor diffusing into the inner disk. As a consequence, close-in sub-Neptunes that accrete some gas from the disk should be dry or wet, respectively, if outer gas giants are outside or inside the water ice line, assuming that giant planets form fast, as has been suggested for Jupiter in our Solar System. Alternatively, a sub-Neptune could form outside the water ice line, accreting a large amount of icy pebbles and then migrating inward as a very wet sub-Neptune. We suggest that the water content of inner sub-Neptunes in systems with giant planets that can efficiently block the inward drifting pebbles could constrain the formation conditions of these systems, thus making these sub-Neptunes exciting targets for detailed characterization (e.g., with JWST, ELT, or ARIEL). In addition, the search for giant planets in systems with already characterized sub-Neptunes can be used to constrain the formation conditions of giant planets as well.
△ Less
Submitted 3 May, 2021; v1 submitted 23 April, 2021;
originally announced April 2021.
-
Progenitor Dependence of Hadron-quark Phase Transition in Failing Core-collapse Supernovae
Authors:
Shuai Zha,
Evan P. O'Connor,
André da Silva Schneider
Abstract:
We study the consequences of a hadron-quark phase transition (PT) in failing core-collapse supernovae (CCSNe) which give birth to stellar-mass black holes (BH). We perform a suite of neutrino-transport general-relativistic hydrodynamic simulations in spherical symmetry with 21 progenitor models and a hybrid equation of state (EoS) including hadrons and quarks. We find that the effect of the PT on…
▽ More
We study the consequences of a hadron-quark phase transition (PT) in failing core-collapse supernovae (CCSNe) which give birth to stellar-mass black holes (BH). We perform a suite of neutrino-transport general-relativistic hydrodynamic simulations in spherical symmetry with 21 progenitor models and a hybrid equation of state (EoS) including hadrons and quarks. We find that the effect of the PT on the CCSN postbounce dynamics is a function of the bounce compactness parameter $ξ_{2.2}$. For $ξ_{2.2}\gtrsim0.24$, the PT leads to a second dynamical collapse of the protocompact star (PCS). While BH formation starts immediately after this second collapse for models with $ξ_{2.2}\gtrsim0.51$, the PCS experiences a second bounce and oscillations for models with $0.24\lesssimξ_{2.2}\lesssim0.51$. These models emit potent oscillatory neutrino signals with a period of $\sim$ms for tens of ms after the second bounce, which can be a strong indicator of the PT in failing CCSNe if detected in the future. However, no shock revival occurs and BH formation inevitably takes place in our spherically-symmetric simulations. Furthermore, via a diagram of mass-specific entropy evolution of the PCS, the progenitor dependence can be understood through the appearance of third-family of compact stars emerging at large entropy induced by the PT.
△ Less
Submitted 3 March, 2021;
originally announced March 2021.
-
Surface waves in protoplanetary disks induced by outbursts: Concentric rings in scattered light
Authors:
A. D. Schneider,
C. P. Dullemond,
B. Bitsch
Abstract:
Context: Vertically hydrostatic protoplanetary disk models are based on the assumption that the main heating source, stellar irradiation, does not vary much with time. However, it is known that accreting young stars are variable sources of radiation. This is particularly evident for outbursting sources such as EX Lupi and FU Orionis stars. Aim: We investigate how such outbursts affect the vertical…
▽ More
Context: Vertically hydrostatic protoplanetary disk models are based on the assumption that the main heating source, stellar irradiation, does not vary much with time. However, it is known that accreting young stars are variable sources of radiation. This is particularly evident for outbursting sources such as EX Lupi and FU Orionis stars. Aim: We investigate how such outbursts affect the vertical structure of the outer regions of the protoplanetary disk, in particular their appearance in scattered light at optical and near-infrared wavelengths. Methods: We employ the 3D FARGOCA radiation-hydrodynamics code, in polar coordinates, to compute the time-dependent behavior of the axisymmetric disk structure. The outbursting inner disk region is not included explicitly. Instead, its luminosity is added to the stellar luminosity and is thus included in the irradiation of the outer disk regions. For time snapshots of interest we insert the density structure into the RADMC-3D radiative transfer code and compute the appearance of the disk at optical/near-infrared wavelengths. Results: We find that, depending on the amplitude of the outbursts, the vertical structure of the disk can become highly dynamic, featuring circular surface waves of considerable amplitude. These "hills" and "valleys" on the disk's surface show up in the scattered light images as bright and dark concentric rings. Initially these rings are small and act as standing waves, but they subsequently lead to outward propagating waves, like the waves produced by a stone thrown into a pond. These waves continue long after the actual outburst has died out. Conclusions: We propose that some of the multi-ringed structures seen in optical/infrared images of several protoplanetary disks may have their origin in outbursts that occurred decades or centuries ago.
△ Less
Submitted 8 September, 2018;
originally announced September 2018.
-
The Progenitor Dependence of Three-Dimensional Core-Collapse Supernovae
Authors:
C. D. Ott,
L. F. Roberts,
A. da Silva Schneider,
J. M. Fedrow,
R. Haas,
E. Schnetter
Abstract:
We present a first study of the progenitor star dependence of the three-dimensional (3D) neutrino mechanism of core-collapse supernovae. We employ full 3D general-relativistic multi-group neutrino radiation-hydrodynamics and simulate the post-bounce evolutions of progenitors with zero-age main sequence masses of $12$, $15$, $20$, $27$, and $40\,M_\odot$. All progenitors, with the exception of the…
▽ More
We present a first study of the progenitor star dependence of the three-dimensional (3D) neutrino mechanism of core-collapse supernovae. We employ full 3D general-relativistic multi-group neutrino radiation-hydrodynamics and simulate the post-bounce evolutions of progenitors with zero-age main sequence masses of $12$, $15$, $20$, $27$, and $40\,M_\odot$. All progenitors, with the exception of the $12\,M_\odot$ star, experience shock runaway by the end of their simulations. In most cases, a strongly asymmetric explosion will result. We find three qualitatively distinct evolutions that suggest a complex dependence of explosion dynamics on progenitor density structure, neutrino heating, and 3D flow. (1) Progenitors with massive cores, shallow density profiles, and high post-core-bounce accretion rates experience very strong neutrino heating and neutrino-driven turbulent convection, leading to early shock runaway. Accretion continues at a high rate, likely leading to black hole formation. (2) Intermediate progenitors experience neutrino-driven, turbulence-aided explosions triggered by the arrival of density discontinuities at the shock. These occur typically at the silicon/silicon-oxygen shell boundary. (3) Progenitors with small cores and density profiles without strong discontinuities experience shock recession and develop the 3D standing-accretion shock instability (SASI). Shock runaway ensues late, once declining accretion rate, SASI, and neutrino-driven convection create favorable conditions. These differences in explosion times and dynamics result in a non-monotonic relationship between progenitor and compact remnant mass.
△ Less
Submitted 4 December, 2017;
originally announced December 2017.
-
A New Open-Source Nuclear Equation of State Framework based on the Liquid-Drop Model with Skyrme Interaction
Authors:
Andre da Silva Schneider,
Luke F. Roberts,
Christian D. Ott
Abstract:
The equation of state (EOS) of dense matter is an essential ingredient for numerical simulations of core-collapse supernovae and neutron star mergers. The properties of matter near and above nuclear saturation density are uncertain, which translates into uncertainties in astrophysical simulations and their multi-messenger signatures. Therefore, a wide range of EOSs spanning the allowed range of nu…
▽ More
The equation of state (EOS) of dense matter is an essential ingredient for numerical simulations of core-collapse supernovae and neutron star mergers. The properties of matter near and above nuclear saturation density are uncertain, which translates into uncertainties in astrophysical simulations and their multi-messenger signatures. Therefore, a wide range of EOSs spanning the allowed range of nuclear interactions are necessary for determining the sensitivity of these astrophysical phenomena and their signatures to variations in input microphysics. We present a new set of finite temperature EOSs based on experimentally allowed Skyrme forces. We employ a liquid drop model of nuclei to capture the non-uniform phase of nuclear matter at sub-saturation density, which is blended into a nuclear statistical equilibrium EOS at lower densities. We also provide a new, open-source code for calculating EOSs for arbitrary Skyrme parametrizations. We then study the effects of different Skyrme parametrizations on thermodynamical properties of dense astrophysical matter, the neutron star mass-radius relationship, and the core collapse of 15 and 40 solar mass stars.
△ Less
Submitted 27 July, 2017; v1 submitted 5 July, 2017;
originally announced July 2017.