In this work we investigate the impact of calculating non-equilibrium chemical abundances consistently with the temperature structure for the atmospheres of highly-irradiated, close-in gas giant ...exoplanets. Chemical kinetics models have been widely used in the literature to investigate the chemical compositions of hot Jupiter atmospheres which are expected to be driven away from chemical equilibrium via processes such as vertical mixing and photochemistry. All of these models have so far used pressure–temperature (P–T) profiles as fixed model input. This results in a decoupling of the chemistry from the radiative and thermal properties of the atmosphere, despite the fact that in nature they are intricately linked. We use a one-dimensional radiative-convective equilibrium model, ATMO, which includes a sophisticated chemistry scheme to calculate P–T profiles which are fully consistent with non-equilibrium chemical abundances, including vertical mixing and photochemistry. Our primary conclusion is that, in cases of strong chemical disequilibrium, consistent calculations can lead to differences in the P–T profile of up to 100 K compared to the P–T profile derived assuming chemical equilibrium. This temperature change can, in turn, have important consequences for the chemical abundances themselves as well as for the simulated emission spectra. In particular, we find that performing the chemical kinetics calculation consistently can reduce the overall impact of non-equilibrium chemistry on the observable emission spectrum of hot Jupiters. Simulated observations derived from non-consistent models could thus yield the wrong interpretation. We show that this behaviour is due to the non-consistent models violating the energy budget balance of the atmosphere.
The anomalously large radii of strongly irradiated exoplanets have remained a major puzzle in astronomy. Based on a two-dimensional steady-state atmospheric circulation model, the validity of which ...is assessed by comparison to three-dimensional calculations, we reveal a new mechanism, namely the advection of the potential temperature due to mass and longitudinal momentum conservation, a process occurring in the Earth's atmosphere or oceans. In the deep atmosphere, the vanishing heating flux forces the atmospheric structure to converge to a hotter adiabat than the one obtained with 1D calculations, implying a larger radius for the planet. Not only do the calculations reproduce the observed radius of HD 209458b, but also reproduce the observed correlation between radius inflation and irradiation for transiting planets. Vertical advection of potential temperature induced by non-uniform atmospheric heating thus provides a robust mechanism to explain the inflated radii of irradiated hot Jupiters.
We present a new set of solar metallicity atmosphere and evolutionary models for very cool brown dwarfs and self-luminous giant exoplanets, which we term
ATMO
2020. Atmosphere models are generated ...with our state-of-the-art 1D radiative-convective equilibrium code
ATMO
, and are used as surface boundary conditions to calculate the interior structure and evolution of 0.001–0.075
M
⊙
objects. Our models include several key improvements to the input physics used in previous models available in the literature. Most notably, the use of a new H–He equation of state including ab initio quantum molecular dynamics calculations has raised the mass by ~1−2% at the stellar–substellar boundary and has altered the cooling tracks around the hydrogen and deuterium burning minimum masses. A second key improvement concerns updated molecular opacities in our atmosphere model
ATMO
, which now contains significantly more line transitions required to accurately capture the opacity in these hot atmospheres. This leads to warmer atmospheric temperature structures, further changing the cooling curves and predicted emission spectra of substellar objects. We present significant improvement for the treatment of the collisionally broadened potassium resonance doublet, and highlight the importance of these lines in shaping the red-optical and near-infrared spectrum of brown dwarfs. We generate three different grids of model simulations, one using equilibrium chemistry and two using non-equilibrium chemistry due to vertical mixing, all three computed self-consistently with the pressure-temperature structure of the atmosphere. We show the impact of vertical mixing on emission spectra and in colour-magnitude diagrams, highlighting how the 3.5−5.5
μ
m flux window can be used to calibrate vertical mixing in cool T–Y spectral type objects.
The survey of the mid-infrared sky by the Wide-field Infrared Survey Explorer (WISE) led to the discovery of extremely cold, low-mass brown dwarfs, classified as Y dwarfs, which extend the T class to ...lower temperatures. Twenty-four Y dwarfs are known at the time of writing. Here we present improved parallaxes for four of these, determined using Spitzer images. We give new photometry for four late-type T and three Y dwarfs and new spectra of three Y dwarfs, obtained at Gemini Observatory. We also present previously unpublished photometry taken from HST, ESO, Spitzer, and WISE archives of 11 late-type T and 9 Y dwarfs. The near-infrared data are put onto the same photometric system, forming a homogeneous data set for the coolest brown dwarfs. We compare recent models to our photometric and spectroscopic data set. We confirm that nonequilibrium atmospheric chemistry is important for these objects. Nonequilibrium cloud-free models reproduce well the near-infrared spectra and mid-infrared photometry for the warmer Y dwarfs with 425 ≤ Teff (K) ≤ 450. A small amount of cloud cover may improve the model fits in the near-infrared for the Y dwarfs with 325 ≤ Teff (K) ≤ 375. Neither cloudy nor cloud-free models reproduce the near-infrared photometry for the Teff = 250 K Y dwarf W0855. We use the mid-infrared region, where most of the flux originates, to constrain our models of W0855. We find that W0855 likely has a mass of 1.5-8 Jupiter masses and an age of 0.3-6 Gyr. The Y dwarfs with measured parallaxes are within 20 pc of the Sun and have tangential velocities typical of the thin disk. The metallicities and ages we derive for the sample are generally solar-like. We estimate that the known Y dwarfs are 3 to 20 Jupiter-mass objects with ages of 0.6-8.5 Gyr.
Atmospheric modeling of low-gravity (VL-G) young brown dwarfs remains challenging. The presence of very thick clouds is a possible source of this challenge, because of their extremely red ...near-infrared (NIR) spectra, but no cloud models provide a good fit to the data with a radius compatible with the evolutionary models for these objects. We show that cloudless atmospheres assuming a temperature gradient reduction caused by fingering convection provide a very good model to match the observed VL-G NIR spectra. The sequence of extremely red colors in the NIR for atmospheres with effective temperatures from ∼2000 K down to ∼1200 K is very well reproduced with predicted radii typical of young low-gravity objects. Future observations with NIRSPEC and MIRI on the James Webb Space Telescope (JWST) will provide more constraints in the mid-infrared, helping to confirm or refute whether or not the NIR reddening is caused by fingering convection. We suggest that the presence or absence of clouds will be directly determined by the silicate absorption features that can be observed with MIRI. JWST will therefore be able to better characterize the atmosphere of these hot young brown dwarfs and their low-gravity exoplanet analogs.
ABSTRACT This work aims to improve the current understanding of the atmospheres of brown dwarfs, especially cold ones with spectral types T and Y, whose modeling is a current challenge. Silicate and ...iron clouds are believed to disappear at the photosphere at the L/T transition, but cloudless models fail to reproduce correctly the spectra of T dwarfs, advocating for the addition of more physics, e.g., other types of clouds or internal energy transport mechanisms. We use a one-dimensional radiative/convective equilibrium code ATMO to investigate this issue. This code includes both equilibrium and out-of-equilibrium chemistry and solves consistently the PT structure. Included opacity sources are H2-H2, H2-He, H2O, CO, CO2, CH4, NH3, K, Na, and TiO, VO if they are present in the atmosphere. We show that the spectra of Y dwarfs can be accurately reproduced with a cloudless model if vertical mixing and NH3 quenching are taken into account. T dwarf spectra still have some reddening in, e.g., J-H, compared to cloudless models. This reddening can be reproduced by slightly reducing the temperature gradient in the atmosphere. We propose that this reduction of the stabilizing temperature gradient in these layers, leading to cooler structures, is due to the onset of fingering convection, triggered by the destabilizing impact of condensation of very thin dust.
Aims. For many years feedback processes generated by OB-stars in molecular clouds, including expanding ionization fronts, stellar winds, or UV-radiation, have been proposed to trigger subsequent star ...formation. However, hydrodynamic models including radiation and gravity show that UV-illumination has little or no impact on the global dynamical evolution of the cloud. Instead, gravitational collapse of filaments and/or merging of filamentary structures can lead to building up dense high-mass star-forming clumps. However, the overall density structure of the cloud has a large influence on this process, and requires a better understanding. Methods. The Rosette molecular cloud, irradiated by the NGC 2244 cluster, is a template region for triggered star-formation, and we investigated its spatial and density structure by applying a curvelet analysis, a filament-tracing algorithm (DisPerSE), and probability density functions (PDFs) on Herschel column density maps, obtained within the HOBYS key program. Results. The analysis reveals not only the filamentary structure of the cloud but also that all known infrared clusters except one lie at junctions of filaments, as predicted by turbulence simulations. The PDFs of sub-regions in the cloud show systematic differences. The two UV-exposed regions have a double-peaked PDF we interprete as caused by shock compression, while the PDFs of the center and other cloud parts are more complex, partly with a power-law tail. A deviation of the log-normal PDF form occurs at AV ≈ 9m for the center, and around 4m for the other regions. Only the part of the cloud farthest from the Rosette nebula shows a log-normal PDF. Conclusions. The deviations of the PDF from the log-normal shape typically associated with low- and high-mass star-forming regions at AV ≈ 3–4m and 8–10m, respectively, are found here within the very same cloud. This shows that there is no fundamental difference in the density structure of low- and high-mass star-forming regions. We conclude that star-formation in Rosette – and probably in high-mass star-forming clouds in general – is not globally triggered by the impact of UV-radiation. Moreover, star formation takes place in filaments that arose from the primordial turbulent structure built up during the formation of the cloud. Clusters form at filament mergers, but star formation can be locally induced in the direct interaction zone between an expanding H II-region and the molecular cloud.
High-mass stars and star clusters commonly form within hub-filament systems. Monoceros R2 (hereafter Mon R2), at a distance of 830 pc, harbors one of the closest such systems, making it an excellent ...target for case studies.
We investigate the morphology, stability and dynamical properties of the Mon R2 hub-filament system.
We employ observations of the
CO and C
O 1→0 and 2→1 lines obtained with the IRAM-30m telescope. We also use H
column density maps derived from
dust emission observations.
We identified the filamentary network in Mon R2 with the DisPerSE algorithm and characterized the individual filaments as either main (converging into the hub) or secondary (converging to a main filament) filaments. The main filaments have line masses of 30-100
pc
and show signs of fragmentation, while the secondary filaments have line masses of 12-60
pc
and show fragmentation only sporadically. In the context of Ostriker's hydrostatic filament model, the main filaments are thermally supercritical. If non-thermal motions are included, most of them are trans-critical. Most of the secondary filaments are roughly transcritical regardless of whether non-thermal motions are included or not. From the morphology and kinematics of the main filaments, we estimate a mass accretion rate of 10
-10
yr
into the central hub. The secondary filaments accrete into the main filaments with a rate of 0.1-0.4×10
yr
. The main filaments extend into the central hub. Their velocity gradients increase towards the hub, suggesting acceleration of the gas.We estimate that with the observed infall velocity, the mass-doubling time of the hub is ~ 2:5 Myr, ten times larger than the free-fall time, suggesting a dynamically old region. These timescales are comparable with the chemical age of the Hii region. Inside the hub, the main filaments show a ring- or a spiral-like morphology that exhibits rotation and infall motions. One possible explanation for the morphology is that gas is falling into the central cluster following a spiral-like pattern.
We report the novel detection of complex high column density tails in the probability distribution functions (PDFs) for three high-mass star-forming regions (CepOB3, MonR2, NGC 6334), obtained from ...dust emission observed with Herschel. The low column density range can be fitted with a lognormal distribution. A first power-law tail starts above an extinction (A
V) of ∼6–14. It has a slope of α = 1.3–2 for the ρ ∝ r
−α profile for an equivalent density distribution (spherical or cylindrical geometry), and is thus consistent with free-fall gravitational collapse. Above A
V ∼40, 60, and 140, we detect an excess that can be fitted by a flatter power-law tail with α > 2. It correlates with the central regions of the cloud (ridges/hubs) of size ∼1 pc and densities above 104 cm−3. This excess may be caused by physical processes that slow down collapse and reduce the flow of mass towards higher densities. Possible are: (1) rotation, which introduces an angular momentum barrier, (2) increasing optical depth and weaker cooling, (3) magnetic fields, (4) geometrical effects, and (5) protostellar feedback. The excess/second power-law tail is closely linked to high-mass star-formation though it does not imply a universal column density threshold for the formation of (high-mass) stars.
We present here a reanalysis of the Spitzer Space Telescope phase curves of the hot Jupiter WASP43 b, using the wavelet pixel-independent component analysis, a blind signal-source separation method. ...The data analyzed were recorded with the Infrared Array Camera and consisted of two visits at 3.6 m, and one visit at 4.5 m, each visit covering one transit and two eclipse events. To test the robustness of our technique we repeated the analysis on smaller portions of the phase curves, and by employing different instrument ramp models. Our reanalysis presents significant updates of the planetary parameters compared to those reported in the original phase curve study of WASP43 b. In particular, we found (1) higher nightside temperatures, (2) smaller hotspot offsets, (3) a greater consistency (∼1 ) between the two 3.6 m visits, and (4) a greater similarity with the predictions of the atmospheric circulation models. Our parameter results are consistent within 1 with those reported by a recent reanalysis of the same data sets. For each visit we studied the variation of the retrieved transit parameters as a function of various sets of stellar limb-darkening coefficients, finding significant degeneracy between the limb-darkening models and the analysis output. Furthermore, we performed the analysis of the single transit and eclipse events, and we examined the differences between these results with the ones obtained with the whole phase curve. Finally we provide a formula useful to optimize the trade-off between precision and duration of observations of transiting exoplanets.