Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reproducibility”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

A Model of Solar Equilibrium: The Hydrodynamic Limit

Helioseismology has revealed the internal density and rotation profiles of the Sun. Yet, knowledge of its magnetic fields and meridional circulation is confined much closer to the surface, and latitudinal entropy gradients are below detectable limits. While numerical simulations can offer insight into the interior dynamics and help identify which ingredients are necessary to reproduce particular observations, some features of the Sun can be understood analytically from an equilibrium perspective. Examples of such features include: the 1D density profile arising from steady-state energy transport from the core to the surface, and the tilting of isorotation contours in the convection zone (CZ) due to baroclinic forcing. To help answer the question of which features can be explained by equilibrium, we propose analyzing stationary axisymmetric ideal magnetohydrodynamic flows in the solar regime. By prescribing an appropriate entropy at the surface, we recover a rotation profile that reasonably matches observations in the bulk of the CZ. Additionally, by including the effects of poloidal flow, we reproduce a feature that is reminiscent of the near surface shear layer. Furthermore, no tachocline-like feature is seen in hydrodynamic equilibrium, suggesting the importance of either dynamics or magnetic fields in its description.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Radiation Hydrodynamical Simulations of the First Quasars

Supermassive black holes (SMBHs) are the central engines of luminous quasars and are found in most massive galaxies today. But the recent discoveries of ULAS J1120+0641, a 2 × 10 9 ${M}_{\odot }$ black hole (BH) at z = 7.1, and ULAS J1342+0928, a 8.0 × 10 8 ${M}_{\odot }$ BH at z = 7.5, now push the era of quasar formation up to just 690 Myr after the Big Bang. Here we report new cosmological simulations of SMBHs with X-rays fully coupled to primordial chemistry and hydrodynamics which show that J1120 and J1342 can form from direct collapse black holes if their growth is fed by cold, dense accretion streams, like those thought to fuel rapid star formation in some galaxies at later epochs. Our models reproduce all of the observed properties of J1120: its mass, luminosity, and H ii region as well as star formation rates and metallicities in its host galaxy. They also reproduce the dynamical mass of the innermost 1.5 kpc of its emission region recently measured by ALMA and J-band magnitudes that are in good agreement with those found by the VISTA Hemisphere Survey.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Constraining the Evolution of the Proton Distribution Function in the Heliotail

We use Interstellar Boundary Explorer (IBEX) measurements of energetic neutral atoms (ENAs) to constrain the proton (mostly pickup ion, PUI) distribution in the heliotail. In our previous study, we solved the Parker transport equation and found that the velocity diffusion coefficient $D(v)$ for PUIs is approximately $D(v)$ ~ 1.1 × 10 -8 km 2 s -3 ($v/v_0$) 1.3 , assuming the initial proton distribution processed by the termination shock (TS), $f_{p,0}$, is a kappa distribution with kappa index $ₖ_{p,0}$ = 1.63. In this study, we test different forms for $f_{p,0}$. We find that if $f_{p,0}$ is kappa-distributed and $D(v)$ = $D_0(v/v_0)^{1.3}$, any kappa index in the range 1.5 < $κ_{p,0}$ < 10 is consistent with IBEX data if $D_0 ~$ 0.8–1.3 × 10 -8 km 2 s -3 . While the case where $D(v)$ ∝ $v^{1.3}$ yields ENA fluxes that appear to best reproduce IBEX data for any $κ_{p,0}$, it is possible for $D(v)$ to scale close to $~v^{2/3}$ or $~v^2$ within our uncertainties by changing $D_0$. We also show that an upstream PUI filled-shell distribution that is heated by a quasi-stationary TS, generating a downstream filled-shell with large cutoff speed, yields an excess of ENAs > 2 keV compared to IBEX. However, using a fully kinetic particle-in-cell simulation to process a PUI filled-shell across the TS yields ENA spectra consistent with IBEX , reinforcing the significance of self-consistent, preferential PUI heating and diffusion at the TS. Interestingly, an upstream PUI distribution inferred from the particle-in-cell simulation to reproduce Voyager 2 observations of the nose-ward TS is inconsistent with IBEX observations from the heliotail, suggesting differences in the upstream PUI distribution or TS properties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

PUSHing Core-collapse Supernovae to Explosions in Spherical Symmetry. II. Explodability and Remnant Properties

In a previously presented proof-of-principle study, we established a parameterized spherically symmetric explosion method (PUSH) that can reproduce many features of core-collapse supernovae. The present paper goes beyond a specific application that is able to reproduce observational properties of SN 1987A and performs a systematic study of an extensive set of nonrotating, solar metallicity stellar progenitor models in the mass range from 10.8 to 120 M ⊙ . This includes the transition from neutron stars to black holes as the final result of the collapse of massive stars, and the relation of the latter to supernovae, possibly faint supernovae, and failed supernovae. We discuss the explosion properties of all models and predict remnant mass distributions within this approach. The present paper provides the basis for extended nucleosynthesis predictions in a forthcoming paper to be employed in galactic evolution models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Actinide Production in the Neutron-rich Ejecta of a Neutron Star Merger

The rapid neutron-capture ("r-") process is responsible for synthesizing many of the heavy elements observed in both the solar system and Galactic metal-poor halo stars. Simulations of r-process nucleosynthesis can reproduce abundances derived from observations with varying success, but so far they fail to account for the observed overenhancement of actinides, present in about 30% of r-process-enhanced stars. In this work, we investigate actinide production in the dynamical ejecta of a neutron star merger (NSM) and explore whether varying levels of neutron-richness can reproduce the actinide boost. We also investigate the sensitivity of actinide production on nuclear physics properties: fission distribution, β-decay, and mass model. For most cases, the actinides are overproduced in our models if the initial conditions are sufficiently neutron-rich for fission cycling. Here, we find that actinide production can be so robust in the dynamical ejecta that an additional lanthanide-rich, actinide-poor component is necessary in order to match observations of actinide-boost stars. We present a simple actinide-dilution model that folds in estimated contributions from two nucleosynthetic sites within a merger event. Our study suggests that while the dynamical ejecta of an NSM are likely production sites for the formation of actinides, a significant contribution from another site or sites (e.g., the NSM accretion disk wind) is required to explain abundances of r-process-enhanced, metal-poor stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Neutron Star Mergers Might Not Be the Only Source of r -process Elements in the Milky Way

Probing the origin of r -process elements in the universe represents a multidisciplinary challenge. We review the observational evidence that probes the properties of r -process sites, and address them using galactic chemical evolution simulations, binary population synthesis models, and nucleosynthesis calculations. Our motivation is to define which astrophysical sites have significantly contributed to the total mass of r -process elements present in our Galaxy. We found discrepancies with the neutron star (NS–NS) merger scenario. When we assume that they are the only site, the decreasing trend of [Eu/Fe] at [Fe/H] > -1 in the disk of the Milky Way cannot be reproduced while accounting for the delay-time distribution (DTD) of coalescence times ( ∝t -1 ) derived from short gamma-ray bursts (GRBs) and population synthesis models. Steeper DTD functions ( ∝t -1.5 ) or power laws combined with a strong burst of mergers before the onset of supernovae (SNe) Ia can reproduce the [Eu/Fe] trend, but this scenario is inconsistent with the similar fraction of short GRBs and SNe Ia occurring in early-type galaxies, and it reduces the probability of detecting GW170817 in an early-type galaxy. One solution is to assume an additional production site of Eu that would be active in the early universe, but would fade away with increasing metallicity. If this is correct, this additional site could be responsible for roughly 50% of the Eu production in the early universe before the onset of SNe Ia. Rare classes of supernovae could be this additional r -process source, but hydrodynamic simulations still need to ensure the conditions for a robust r -process pattern.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

A New Model for Electron-capture Supernovae in Galactic Chemical Evolution

We examine the contribution of electron-capture supernovae (ECSNe), low-mass SNe from collapsing Fe cores (FeCCSNe), and rotating massive stars to the chemical composition of the Galaxy. Our model includes contributions to chemical evolution from both thermonuclear ECSNe (tECSNe) and gravitational collapse ECSNe (cECSNe). We show that if ECSNe are predominantly gravitational collapse SNe but about 15% are partial thermonuclear explosions, the model is able to reproduce the solar abundances of several important and problematic isotopes including ${}^{48}\mathrm{Ca}$, ${}^{50}\mathrm{Ti}$, and 54 Cr together with 58 Fe, 64 Ni, 82 Se, and 86 Kr and several of the Zn–Zr isotopes. A model in which no cECSNe occur, only tECSNe with low-mass FeCCSNe or rotating massive stars, proves also very successful at reproducing the solar abundances for these isotopes. Despite the small mass range for the progenitors of ECSNe and low-mass FeCCSNe, the large production factors suffice for the solar inventory of the above isotopes. Our model is compelling because it introduces no new tensions with the solar abundance distribution for a Milky Way model—only tending to improve the model predictions for several isotopes. Here, the proposed astrophysical production model thus provides a natural and elegant way to explain one of the last uncharted territories on the periodic table of astrophysical element production.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The Polarized Image of a Synchrotron-emitting Ring of Gas Orbiting a Black Hole

Synchrotron radiation from hot gas near a black hole results in a polarized image. The image polarization is determined by effects including the orientation of the magnetic field in the emitting region, relativistic motion of the gas, strong gravitational lensing by the black hole, and parallel transport in the curved spacetime. We explore these effects using a simple model of an axisymmetric, equatorial accretion disk around a Schwarzschild black hole. By using an approximate expression for the null geodesics derived by Beloborodov and conservation of the Walker–Penrose constant, we provide analytic estimates for the image polarization. We test this model using currently favored general relativistic magnetohydrodynamic simulations of M87*, using ring parameters given by the simulations. For a subset of these with modest Faraday effects, we show that the ring model broadly reproduces the polarimetric image morphology. Our model also predicts the polarization evolution for compact flaring regions, such as those observed from Sgr A* with GRAVITY. With suitably chosen parameters, our simple model can reproduce the EVPA pattern and relative polarized intensity in Event Horizon Telescope images of M87*. Under the physically motivated assumption that the magnetic field trails the fluid velocity, this comparison is consistent with the clockwise rotation inferred from total intensity images.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Ultracool Dwarfs Observed with the Spitzer Infrared Spectrograph. I. An Accurate Look at the L-to-T Transition at ~300 Myr from Optical Through Mid-infrared Spectrophotometry

Here we present Spitzer IRS 5–14 μm spectra and 16 μm and 22 μm photometry of the T2.5 companion to the ~300 Myr old G0V star HN Peg. We incorporate previous 0.8–5 μm observations to obtain the most comprehensive spectral energy distribution (SED) of an intermediate-gravity L/T-transition dwarf that, together with an accurate Gaia EDR3 parallax of the primary, enables us to derive precise fundamental parameters. We find that young (≈0.1–0.3 Gyr) early-T dwarfs on average have ≈140 K lower effective temperatures, ≈20% larger radii, and similar bolometric luminosities compared to ≳1 Gyr old field dwarfs with similar spectral types. Our accurate infrared spectrophotometry offers new detail at wavelengths where the dominant carbon-bearing molecules have their strongest transitions: at 3.4 μm for methane and at 4.6 μm for carbon monoxide. We assess the performance of various widely available photospheric models and find that models with condensates and/or clouds better reproduce the full SED of this moderately young early-T dwarf. However, cloud-free models incorporating a more general convective instability treatment reproduce at least the low-resolution near-infrared spectrum similarly well. Our analysis of R ≈ 2300 J-band spectra shows that the near-infrared potassium absorption lines in HN Peg B have similar strengths to those seen in both younger and older T2–T3 dwarfs. We conclude that while alkali lines are well established as surface gravity indicators for L-type or warmer stars, they are insensitive to surface gravity in early-T dwarfs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

ADDGALS: Simulated Sky Catalogs for Wide Field Galaxy Surveys

Abstract We present a method for creating simulated galaxy catalogs with realistic galaxy luminosities, broadband colors, and projected clustering over large cosmic volumes. The technique, denoted Addgals (Adding Density Dependent GAlaxies to Lightcone Simulations), uses an empirical approach to place galaxies within lightcone outputs of cosmological simulations. It can be applied to significantly lower-resolution simulations than those required for commonly used methods such as halo occupation distributions, subhalo abundance matching, and semi-analytic models, while still accurately reproducing projected galaxy clustering statistics down to scales of r ∼ 100 h −1 kpc . We show that Addgals catalogs reproduce several statistical properties of the galaxy distribution as measured by the Sloan Digital Sky Survey (SDSS) main galaxy sample, including galaxy number densities, observed magnitude and color distributions, as well as luminosity- and color-dependent clustering. We also compare to cluster–galaxy cross correlations, where we find significant discrepancies with measurements from SDSS that are likely linked to artificial subhalo disruption in the simulations. Applications of this model to simulations of deep wide-area photometric surveys, including modeling weak-lensing statistics, photometric redshifts, and galaxy cluster finding, are presented in DeRose et al., and an application to a full cosmology analysis of Dark Energy Survey (DES) Year 3 like data is presented in DeRose et al. We plan to publicly release a 10,313 square degree catalog constructed using Addgals with magnitudes appropriate for several existing and planned surveys, including SDSS, DES, VISTA, Wide-field Infrared Survey Explorer, and Rubin Observatory’s Legacy Survey of Space and Time.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

BASS. XXXVI. Constraining the Local Supermassive Black Hole–Halo Connection with BASS DR2 AGNs

We investigate the connection between supermassive black holes (SMBHs) and their host dark matter halos in the local universe using the clustering statistics and luminosity function of active galactic nuclei (AGNs) from the Swift/BAT AGN Spectroscopic Survey (BASS DR2). By forward-modeling AGN activity into snapshot halo catalogs from N -body simulations, we test a scenario in which SMBH mass correlates with dark matter (sub)halo mass for fixed stellar mass. We compare this to a model absent of this correlation, where stellar mass alone determines the SMBH mass. We find that while both simple models are able to largely reproduce the abundance and overall clustering of AGNs, the model in which black hole mass is tightly correlated with halo mass is preferred by the data by 1.8 σ . When including an independent measurement on the black hole mass–halo mass correlation, this model is preferred by 4.6 σ . We show that the clustering trends with black hole mass can further break the degeneracies between the two scenarios and that our preferred model reproduces the measured clustering differences on one-halo scales between large and small black hole masses. These results indicate that the halo binding energy is fundamentally connected to the growth of SMBHs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The Evolution and Mass Dependence of Galaxy Cluster Pressure Profiles at 0.05 ≤ z ≤ 0.60 and 4 × 10 14 M ⊙ ≤ M 500 ≤ 30 × 10 14 M ⊙

Abstract We have combined X-ray observations from Chandra with Sunyaev–Zel’dovich effect data from Planck and Bolocam to measure intracluster medium pressure profiles from 0.03 R 500 ≤ R ≤ 5 R 500 for a sample of 21 low- z galaxy clusters with a median redshift of 〈 z 〉 = 0.08 and a median mass of 〈 M 500 〉 = 6.1 × 10 14 M ⊙ and a sample of 19 mid- z galaxy clusters with 〈 z 〉 = 0.50 and 〈 M 500 〉 = 10.6 × 10 14 M ⊙ . The mean scaled pressure in the low- z sample is lower at small radii and higher at large radii, a trend that is accurately reproduced in similarly selected samples from The Three Hundred simulations. This difference appears to be primarily due to dynamical state at small radii, evolution at intermediate radii, and a combination of evolution and mass dependence at large radii. Furthermore, the overall flattening of the mean scaled pressure profile in the low- z sample compared to the mid- z sample is consistent with expectations due to differences in the mass accretion rate and the fractional impact of feedback mechanisms. In agreement with previous studies, the fractional scatter about the mean scaled pressure profile reaches a minimum of ≃20% near 0.5 R 500 . This scatter is consistent between the low- z and mid- z samples at all radii, suggesting it is not strongly impacted by sample selection, and this general behavior is reproduced in The Three Hundred simulations. Finally, analytic functions that approximately describe the mass and redshift trends in mean pressure profile shape are provided.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The Progenitors of Superluminous Type Ia Supernovae

Abstract Recent observations of type Ia supernovae (SNe Ia) have discovered a subclass of “super-Chandrasekhar” SNe Ia (SC SNe Ia) whose high luminosities and low ejecta velocities suggest that they originate from the explosions of white dwarfs (WDs) with masses that exceed the Chandrasekhar mass limit. Different models have been proposed to explain the progenitors of these explosions, including a “magnetized WD” model and a “WD merger” model. To test the robustness of these models, we conduct a 1D numerical parameter survey of WD explosions using these models as initial conditions. We follow the explosions using the hydrodynamics code Castro and then use the radiation transport code SuperNu to create light curves and spectra for the models. We find that while both classes of models fall within the range of SC SNe Ia observations on the light-curve width–luminosity relation, only the WD merger models reproduce the observed low ejecta velocities. The light curves of our merger models are more similar photometrically to observations than our magnetized models. Given this, we discuss possible explanations for the brightest SC SNe Ia observations that cannot be reproduced with our WD merger models. This study provides the basis for future SC SNe Ia observations and higher-dimensional numerical models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Collapsar Gamma-Ray Bursts Grind Their Black Hole Spins to a Halt

Abstract The spin of a newly formed black hole (BH) at the center of a massive star evolves from its natal value due to two competing processes: accretion of gas angular momentum that increases the spin and extraction of BH angular momentum by outflows that decreases the spin. Ultimately, the final, equilibrium spin is set by a balance between both processes. In order for the BH to launch relativistic jets and power aγ-ray burst (GRB), the BH magnetic field needs to be dynamically important. Thus, we consider the case of a magnetically arrested disk (MAD) driving the spin evolution of the BH. By applying the semianalytic MAD BH spin evolution model of Lowell et al. to collapsars, we show that if the BH accretes ∼20% of its initial mass, its dimensionless spin inevitably reaches small values,a≲ 0.2. For such spins, and for mass accretion rates inferred from collapsar simulations, we show that our semianalytic model reproduces the energetics of typical GRB jets,L jet ∼ 10 50 erg s −1 . We show that our semianalytic model reproduces the nearly constant power of typical GRB jets. If the MAD onset is delayed, this allows powerful jets at the high end of the GRB luminosity distribution,L jet ∼ 10 52 erg s −1 , but the final spin remains low,a≲ 0.3. These results are consistent with the low spins inferred from gravitational wave detections of binary BH mergers. In a companion paper by Gottlieb et al., we use GRB observations to constrain the natal BH spin to bea≃ 0.2.

Astronomy & Astrophysics↗
DOE OSTI2024

Neutrino-driven Outflows and the Elemental Abundance Patterns of Very Metal-poor Stars

Abstract The elemental abundances between strontium and silver ( Z = 38–47) observed in the atmospheres of very metal-poor stars in the Galaxy may contain the fingerprint of the weak r -process and ν p -process occurring in early core-collapse supernovae explosions. In this work, we combine various astrophysical conditions based on a steady-state model to cover the richness of the supernova ejecta in terms of entropy, expansion timescale, and electron fraction. The calculated abundances based on different combinations of conditions are compared with stellar observations, with the aim of constraining supernova ejecta conditions. We find that some conditions of the neutrino-driven outflows consistently reproduce the observed abundances of our sample. In addition, from the successful combinations, the neutron-rich trajectories better reproduce the observed abundances of Sr–Zr ( Z = 38–40), while the proton-rich ones, Mo–Pd ( Z = 42–47).

Astronomy & Astrophysics↗
DOE OSTI2024

Mercury’s Chaotic Secular Evolution as a Subdiffusive Process

Abstract Mercury’s orbit can destabilize, generally resulting in a collision with either Venus or the Sun. Chaotic evolution can causeg 1 to decrease to the approximately constant value ofg 5 and create a resonance. Previous work has approximated the variation ing 1 as stochastic diffusion, which leads to a phenomological model that can reproduce the Mercury instability statistics of secular andN-body models on timescales longer than 10 Gyr. Here we show that the diffusive model significantly underpredicts the Mercury instability probability on timescales less than 5 Gyr, the remaining lifespan of the solar system. This is becauseg 1 exhibits larger variations on short timescales than the diffusive model would suggest. To better model the variations on short timescales, we build a new subdiffusive phenomological model forg 1 . Subdiffusion is similar to diffusion but exhibits larger displacements on short timescales and smaller displacements on long timescales. We choose model parameters based on the behavior of theg 1 trajectories in theN-body simulations, leading to a tuned model that can reproduce Mercury instability statistics from 1–40 Gyr. This work motivates fundamental questions in solar system dynamics: why does subdiffusion better approximate the variation ing 1 than standard diffusion? Why is there an upper bound ong 1 , but not a lower bound that would prevent it from reachingg 5 ?

Astronomy & Astrophysics↗
DOE OSTI2017

DISCERNING THE GAMMA-RAY-EMITTING REGION IN THE FLAT SPECTRUM RADIO QUASARS

A model-dependent method is proposed to determine the location of the γ -ray-emitting region for a given flat spectrum radio quasar (FSRQ). In the model, the extra-relativistic electrons are injected at the base of the jet and non-thermal photons are produced by both synchrotron radiation and inverse-Compton (IC) scattering in the energy dissipation region. The target photons dominating IC scattering originate from both synchrotron photons and external ambient photon fields, and the energy density of external radiation field is a function of the distance between the position of the dissipation region and a central supermassive black hole, and their spectra are seen in the comoving frame. Moreover, the energy dissipation region could be determined by the model parameter through reproducing the γ -ray spectra. Such a model is applied to reproduce the quasi-simultaneous multi-wavelength observed data for 36 FSRQs. In order to define the width of the broad-line region (BLR) shell and dusty molecular torus (MT) shell, a simple numerical constraint is used to determine the outer boundary of the BLR and dusty MT. Our results show that (1) the γ -ray-emitting regions are located at the range from 0.1 to 10 pc; (2) the γ -ray-emitting regions are located outside the BLRs and within the dusty molecular tori; and (3) the γ -ray-emitting regions are located closer to the dusty MT ranges than the BLRs. Therefore, it may be concluded that direct evidence for the far site scenario could be obtained on the basis of the model results.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

VULCAN: An Open-source, Validated Chemical Kinetics Python Code for Exoplanetary Atmospheres

We present an open-source and validated chemical kinetics code for studying hot exoplanetary atmospheres, which we name VULCAN. It is constructed for gaseous chemistry from 500 to 2500 K, using a reduced C–H–O chemical network with about 300 reactions. It uses eddy diffusion to mimic atmospheric dynamics and excludes photochemistry. We have provided a full description of the rate coefficients and thermodynamic data used. We validate VULCAN by reproducing chemical equilibrium and by comparing its output versus the disequilibrium-chemistry calculations of Moses et al. and Rimmer and Helling. It reproduces the models of HD 189733b and HD 209458b by Moses et al., which employ a network with nearly 1600 reactions. We also use VULCAN to examine the theoretical trends produced when the temperature–pressure profile and carbon-to-oxygen ratio are varied. Assisted by a sensitivity test designed to identify the key reactions responsible for producing a specific molecule, we revisit the quenching approximation and find that it is accurate for methane but breaks down for acetylene, because the disequilibrium abundance of acetylene is not directly determined by transport-induced quenching, but is rather indirectly controlled by the disequilibrium abundance of methane. Therefore we suggest that the quenching approximation should be used with caution and must always be checked against a chemical kinetics calculation. A one-dimensional model atmosphere with 100 layers, computed using VULCAN, typically takes several minutes to complete. VULCAN is part of the Exoclimes Simulation Platform (ESP; exoclime.net) and publicly available at https://github.com/exoclime/VULCAN.

79 ASTRONOMY AND ASTROPHYSICS↗