In-order to understand the results of recent observations of exoplanets, models have become increasingly complex. Unfortunately this increases both the computational cost and output size of said ...models. We intend to explore if AI-image-recognition can alleviate this burden. We used DYNAMICO to run a series of HD209458-like models with different orbital-radii. Training data for a number of features of interest was selected from the initial outputs of these models. This was used to train a pair of multi-categorisation convolutional-neural-networks (CNN), which we applied to our outer-atmosphere-equilibrated models. The features detected by our CNNs revealed that our models fall into two regimes: models with a shorter orbital-radii exhibit significant global mixing which shapes the entire atmospheres dynamics. Whereas, models with longer orbital-radii exhibit negligible mixing except at mid-pressures. Here, the initial non-detection of any trained features revealed a surprise: a night-side hot-spot. Analysis suggests that this occurs when rotational influence is sufficiently weak that divergent flows from the day-side to the night-side dominate over rotational-driven transport, such as the equatorial jet. We suggest that image-classification may play an important role in future, computational, atmospheric studies. However special care must be paid to the data feed into the model, from the colourmap, to training the CNN on features with enough breadth and complexity that the CNN can learn to detect them. However, by using preliminary-studies and prior-models, this should be more than achievable for future exascale calculations, allowing for a significant reduction in future workloads and computational resources.
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.
To study the complexity of hot Jupiter atmospheres revealed by observations of increasing quality, we have adapted the UK Met Office Global Circulation Model (GCM), the Unified Model (UM), to these ...exoplanets. The UM solves the full 3D Navier-Stokes equations with a height-varying gravity, avoiding the simplifications used in most GCMs currently applied to exoplanets. In this work we present the coupling of the UM dynamical core to an accurate radiation scheme based on the two-stream approximation and correlated-k method with state-of-the-art opacities from ExoMol. Our first application of this model is devoted to the extensively studied hot Jupiter HD 209458b. We have derived synthetic emission spectra and phase curves, and compare them to both previous models also based on state-of-the-art radiative transfer, and to observations. We find a reasonable agreement between observations and both our days side emission and hot spot offset, however, our night side emissions is too large. Overall our results are qualitatively similar to those found by Showman et al. (2009, ApJ, 699, 564) with the SPARC/MITgcm, however, we note several quantitative differences: Our simulations show significant variation in the position of the hottest part of the atmosphere with pressure, as expected from simple timescale arguments, and in contrast to the “vertical coherency” found by Showman et al. (2009). We also see significant quantitative differences in calculated synthetic observations. Our comparisons strengthen the need for detailed intercomparisons of dynamical cores, radiation schemes and post-processing tools to understand these differences. This effort is necessary in order to make robust conclusions about these atmospheres based on GCM results.
By generalizing the theory of convection to any type of thermal and compositional source terms (diabatic processes), we show that thermohaline convection in Earth's oceans, fingering convection in ...stellar atmospheres, and moist convection in Earth's atmosphere are derived from the same general diabatic convective instability. We also show that "radiative convection" triggered by the CO/CH4 transition with radiative transfer in the atmospheres of brown dwarfs is analogous to moist and thermohaline convection. We derive a generalization of the mixing-length theory to include the effect of source terms in 1D codes. We show that CO/CH4 "radiative" convection could significantly reduce the temperature gradient in the atmospheres of brown dwarfs similarly to moist convection in Earth's atmosphere, thus possibly explaining the reddening in brown dwarf spectra. By using idealized 2D hydrodynamic simulations in the Ledoux unstable regime, we show that compositional source terms can indeed provoke a reduction of the temperature gradient. The L/T transition could be explained by a bifurcation between the adiabatic and diabatic convective transports and seen as a giant cooling crisis: an analog of the boiling crisis in liquid/steam-water convective flows. This mechanism, with other chemical transitions, could be present in many giant and Earth-like exoplanets. The study of the impact of different parameters (effective temperature, compositional changes) on CO/CH4 radiative convection and the analogy with Earth moist and thermohaline convection is opening the possibility of using brown dwarfs to better understand some aspects of the physics at play in the climate of our own planet.
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.
ABSTRACT
Transport-induced quenching, i.e. the homogenization of chemical abundances by atmospheric advection, is thought to occur in the atmospheres of hot gas giant exoplanets. While some numerical ...modelling of this process exists, the three-dimensional nature of transport-induced quenching is underexplored. Here, we present results of 3D cloud- and haze-free simulations of the atmospheres of HAT-P-11b, HD 189733b, HD 209458b, and WASP-17b including coupled hydrodynamics, radiative transfer, and chemistry. Our simulations were performed with two chemical schemes: a chemical kinetics scheme, which is capable of capturing transport-induced quenching, and a simpler, more widely used chemical equilibrium scheme. We find that transport-induced quenching is predicted to occur in atmospheres of all planets in our sample; however, the extent to which it affects their synthetic spectra and phase curves varies from planet to planet. This implies that there is a ‘sweet spot’ for the observability of signatures of transport-induced quenching, which is controlled by the interplay between the dynamics and chemistry.
HAT-P-26b Wakeford, Hannah R.; Sing, David K.; Kataria, Tiffany ...
Science (American Association for the Advancement of Science),
05/2017, Volume:
356, Issue:
6338
Journal Article
Peer reviewed
Open access
A correlation between giant-planet mass and atmospheric heavy elemental abundance was first noted in the past century from observations of planets in our own Solar System and has served as a ...cornerstone of planet-formation theory. Using data from the Hubble and Spitzer Space Telescopes from 0.5 to 5 micrometers, we conducted a detailed atmospheric study of the transiting Neptune-mass exoplanet HAT-P-26b. We detected prominent H₂O absorption bands with a maximum base-to-peak amplitude of 525 parts per million in the transmission spectrum. Using the water abundance as a proxy for metallicity, we measured HAT-P-26b’s atmospheric heavy element content (
4.8
−
4.0
+
21.5
times solar). This likely indicates that HAT-P-26b’s atmosphere is primordial and obtained its gaseous envelope late in its disk lifetime, with little contamination from metal-rich planetesimals.
Ultra-hot rocky super-Earths are thought to be sufficiently irradiated by their host star to melt their surface and allow for long-lasting magma oceans as a result. A number of processes have been ...proposed to explain how such planets may have retained the primordial hydrogen captured during their formation, while moving inward in the planetary system. The new generation of space telescopes such as the James Webb Space Telescope may provide observations that are precise enough to characterize the atmospheres and perhaps the interiors of such exoplanets. We used a vaporization model that calculates the gas-liquid equilibrium between the atmosphere (including hydrogen) and the magma ocean to compute the elemental composition of a variety of atmospheres with different quantities of hydrogen. We then used the elemental composition in a steady-state atmospheric model (ATMO) to compute the atmospheric structure and generate synthetic emission spectra. With this method, we were able to confirm previous results showing that silicate atmospheres exhibit a thermal inversion, with a notable emission peak of SiO at 9 μm. We compared our method to the literature on the inclusion of hydrogen in the atmosphere to show that hydrogen reduces the thermal inversion because of the formation of H 2 O, which has a strong greenhouse potential. However, planets that are significantly irradiated by their host star are sufficiently hot to dissociate H 2 O, thus also allowing them to maintain a thermal inversion. The observational implications are twofold: (1) H 2 O is more likely to be detected in colder atmospheres and (2) detecting a thermal inversion in hotter atmospheres does not a priori exclude the presence of H (in its atomic form). Due to the impact of H on the overall chemistry and atmospheric structure (and, thus, observations), we emphasize the importance of including volatiles in the calculation of the gas-liquid equilibrium. Finally, we provide a criterion to determine potential targets for observation in light of these findings.
Abstract
We present a grid of forward model transmission spectra, adopting an isothermal temperature-pressure profile, alongside corresponding equilibrium chemical abundances for 117 observationally ...significant hot exoplanets (equilibrium temperatures of 547–2710 K). This model grid has been developed using a 1D radiative–convective–chemical equilibrium model termed ATMO, with up-to-date high-temperature opacities. We present an interpretation of observations of 10 exoplanets, including best-fitting parameters and χ2 maps. In agreement with previous works, we find a continuum from clear to hazy/cloudy atmospheres for this sample of hot Jupiters. The data for all the 10 planets are consistent with subsolar to solar C/O ratio, 0.005 to 10 times solar metallicity and water rather than methane-dominated infrared spectra. We then explore the range of simulated atmospheric spectra for different exoplanets, based on characteristics such as temperature, metallicity, C/O ratio, haziness and cloudiness. We find a transition value for the metallicity between 10 and 50 times solar, which leads to substantial changes in the transmission spectra. We also find a transition value of C/O ratio, from water to carbon species dominated infrared spectra, as found by previous works, revealing a temperature dependence of this transition point ranging from ∼0.56 to ∼1–1.3 for equilibrium temperatures from ∼900 to ∼2600 K. We highlight the potential of the spectral features of HCN and C2H2 to constrain the metallicities and C/O ratios of planets, using James Webb Space Telescope (JWST) observations. Finally, our entire grid (∼460 000 simulations) is publicly available and can be used directly with the JWST simulator PandExo for planning observations.
starbench: the D-type expansion of an H ii region Bisbas, T. G; Haworth, T. J; Williams, R. J. R ...
Monthly notices of the Royal Astronomical Society,
10/2015, Volume:
453, Issue:
2
Journal Article
Peer reviewed
Open access
starbench is a project focused on benchmarking and validating different star formation and stellar feedback codes. In this first starbench paper we perform a comparison study of the D-type expansion ...of an H ii region. The aim of this work is to understand the differences observed between the 12 participating numerical codes against the various analytical expressions examining the D-type phase of H ii region expansion. To do this, we propose two well-defined tests which are tackled by 1D and 3D grid- and smoothed particle hydrodynamics-based codes. The first test examines the ‘early phase’ D-type scenario during which the mechanical pressure driving the expansion is significantly larger than the thermal pressure of the neutral medium. The second test examines the ‘late phase’ D-type scenario during which the system relaxes to pressure equilibrium with the external medium. Although they are mutually in excellent agreement, all 12 participating codes follow a modified expansion law that deviates significantly from the classical Spitzer solution in both scenarios. We present a semi-empirical formula combining the two different solutions appropriate to both early and late phases that agrees with high-resolution simulations to ≲ 2 per cent. This formula provides a much better benchmark solution for code validation than the Spitzer solution. The present comparison has validated the participating codes and through this project we provide a data set for calibrating the treatment of ionizing radiation hydrodynamics codes.