We present results from a set of simulations using a fully coupled three-dimensional (3D) chemistry-radiation-hydrodynamics model and investigate the effect of transport of chemical species by the large-scale atmospheric flow in hot Jupiter atmospheres. We couple a flexible chemical kinetics scheme to the Met Office Unified Model which enables the study of the interaction of chemistry, radiative transfer and fluid dynamics. We use a newly-released "reduced" chemical network comprising 30 chemical species that has been specifically developed for application in 3D atmosphere models. We simulate the atmospheres of the well-studied hot Jupiters HD~209458b and HD~189733b which both have dayside--nightside temperature contrasts of several hundred Kelvin and superrotating equatorial jets. We find qualitatively quite different chemical structures between the two planets, particularly for methane (CH$_4$), when advection of chemical species is included. Our results show that consideration of 3D chemical transport is vital in understanding the chemical composition of hot Jupiter atmospheres. 3D mixing leads to significant changes in the abundances of absorbing gas-phase species compared with what would be expected by assuming local chemical equilibrium, or from models including 1D - and even 2D - chemical mixing. We find that CH$_4$, carbon dioxide (CO$_2$) and ammonia (NH$_3$) are particularly interesting as 3D mixing of these species leads to prominent signatures of out-of-equilibrium chemistry in the transmission and emission spectra, detectable with near-future instruments.

We present new parallax measurements from the CFHT Infrared Parallax Program and the Pan-STARRS 3$\pi$ Steradian Survey for the young ($\approx150-300$ Myr) triple system VHS J125601.92$-$125723.9. This system is composed of a nearly equal-flux binary ("AB") and a wide, possibly planetary-mass companion ("b"). The system's published parallactic distance ($12.7\pm1.0$ pc) implies absolute magnitudes unusually faint compared to known young objects and is in tension with the spectrophotometric distance for the central binary ($17.2\pm2.6$ pc). Our CFHT and Pan-STARRS parallaxes are consistent, and the more precise CFHT result places VHS J1256-1257 at $22.2^{+1.1}_{-1.2}$ pc. Our new distance results in higher values for the companion's mass ($19\pm5$ M$_{\rm Jup}$) and temperature ($1240\pm50$ K), and also brings the absolute magnitudes of all three components into better agreement with known young objects.

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\,\mathrm{M_{\odot}}$ 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 $\sim1-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\,\mathrm{\mu m}$ flux window can be used to calibrate vertical mixing in cool T-Y spectral type objects.

The James Webb Space Telescope (JWST) is expected to revolutionize the field of exoplanets. The broad wavelength coverage and the high sensitivity of its instruments will allow characterization of exoplanetary atmospheres with unprecedented precision. Following the Call for the Cycle 1 Early Release Science Program, the Transiting Exoplanet Community was awarded time to observe several targets, including WASP-43b. The atmosphere of this hot Jupiter has been intensively observed but still harbors some mysteries, especially concerning the day-night temperature gradient, the efficiency of the atmospheric circulation, and the presence of nightside clouds. We will constrain these properties by observing a full orbit of the planet and extracting its spectroscopic phase curve in the 5--12 $\mu$m range with JWST/MIRI. To prepare for these observations, we performed an extensive modeling work with various codes: radiative transfer, chemical kinetics, cloud microphysics, global circulation models, JWST simulators, and spectral retrieval. Our JWST simulations show that we should achieve a precision of 210 ppm per 0.1 $\mu$m spectral bin on average, which will allow us to measure the variations of the spectrum in longitude and measure the night-side emission spectrum for the first time. If the atmosphere of WASP-43b is clear, our observations will permit us to determine if its atmosphere has an equilibrium or disequilibrium chemical composition, providing eventually the first conclusive evidence of chemical quenching in a hot Jupiter atmosphere. If the atmosphere is cloudy, a careful retrieval analysis will allow us to identify the cloud composition.

Several chemical networks have been developed to study warm (exo)planetary atmospheres. The kinetics of the reactions related to the methanol chemistry included in these schemes have been questioned. The goal of this paper is to update the methanol chemistry for such chemical networks thanks to recent publications in the combustion literature. We aim also at studying the consequences of this update on the atmospheric compositions of (exo)planetary atmospheres and brown dwarfs. We have performed an extensive review of combustion experimental studies and revisited the sub-mechanism describing methanol combustion in the scheme of Venot et al. (2012, A&A 624, A58). The updated scheme involves 108 species linked by a total of 1906 reactions. We have then applied our 1D kinetic model with this new scheme to several case studies (HD 209458b, HD 189733b, GJ 436b, GJ 1214b, ULAS J1335+11, Uranus, Neptune), and compared the results obtained with those obtained with the former scheme. The update of the scheme has a negligible impact on hot Jupiters atmospheres. However, the atmospheric composition of warm Neptunes and brown dwarfs is modified sufficiently to impact observational spectra in the wavelength range JWST will operate. Concerning Uranus and Neptune, the update of the chemical scheme modifies the abundance of CO and thus impacts the deep oxygen abundance required to reproduce the observational data. For future 3D kinetics models, we also derived a reduced scheme containing 44 species and 582 reactions. Chemical schemes should be regularly updated in order to maintain a high level of reliability on the results of kinetic models and be able to improve our knowledge on planetary formation.

Context: The anomalously large radii of hot Jupiters has long been a mystery. However, by combining both theoretical arguments and 2D models, a recent study has suggested that the vertical advection of potential temperature leads to an adiabatic temperature profile in the deep atmosphere hotter than the profile obtained with standard 1D models. Aims: In order to confirm the viability of that scenario, we extend this investigation to three dimensional, time-dependent, models. Methods: We use a 3D GCM, DYNAMICO to perform a series of calculations designed to explore the formation and structure of the driving atmospheric circulations, and detail how it responds to changes in both upper and deep atmospheric forcing. Results: In agreement with the previous, 2D, study, we find that a hot adiabat is the natural outcome of the long-term evolution of the deep atmosphere. Integration times of order $1500$ years are needed for that adiabat to emerge from an isothermal atmosphere, explaining why it has not been found in previous hot Jupiter studies. Models initialised from a hotter deep atmosphere tend to evolve faster toward the same final state. We also find that the deep adiabat is stable against low-levels of deep heating and cooling, as long as the Newtonian cooling time-scale is longer than $\sim 3000$ years at $200$ bar. Conclusions: We conclude that the steady-state vertical advection of potential temperature by deep atmospheric circulations constitutes a robust mechanism to explain hot Jupiter inflated radii. We suggest that future studies of hot Jupiters are evolved for a longer time than currently done, and, when possible, include models initialised with a hot deep adiabat. We stress that this mechanism stems from the advection of entropy by irradiation induced mass flows and does not require (finely tuned) dissipative process, in contrast with most previously suggested scenarios.

We present individual dynamical masses for the nearby M9.5+T5.5 binary WISE J072003.20$-$084651.2AB, a.k.a. Scholz's star. Combining high-precision CFHT/WIRCam photocenter astrometry and Keck adaptive optics resolved imaging, we measure the first high-quality parallactic distance ($6.80_{-0.06}^{+0.05}$ pc) and orbit ($8.06_{-0.25}^{+0.24}$ yr period) for this system composed of a low-mass star and brown dwarf. We find a moderately eccentric orbit ($e = 0.240_{-0.010}^{+0.009}$), incompatible with previous work based on less data, and dynamical masses of $99\pm6$ $M_{\rm Jup}$ and $66\pm4$ $M_{\rm Jup}$ for the two components. The primary mass is marginally inconsistent (2.1$\sigma$) with the empirical mass$-$magnitude$-$metallicity relation and models of main-sequence stars. The relatively high mass of the cold ($T_{\rm eff} = 1250\pm40$ K) brown dwarf companion indicates an age older than a few Gyr, in accord with age estimates for the primary star, and is consistent with our recent estimate of $\approx$70 $M_{\rm Jup}$ for the stellar/substellar boundary among the field population. Our improved parallax and proper motion, as well as an orbit-corrected system velocity, improve the accuracy of the system's close encounter with the solar system by an order of magnitude. WISE J0720$-$0846AB passed within $68.7\pm2.0$ kAU of the Sun $80.5\pm0.7$ kyr ago, passing through the outer Oort cloud where comets can have stable orbits.

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 $\mu$m, and one visit at 4.5 $\mu$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 ($\sim$1 $\sigma$) between the two 3.6~$\mu$m visits, and (4) a greater similarity with the predictions of the atmospheric circulation models. Our parameter results are consistent within 1 $\sigma$ 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.

WASP-121b is a transiting gas giant exoplanet orbiting close to its Roche limit, with an inflated radius nearly double that of Jupiter and a dayside temperature comparable to a late M dwarf photosphere. Secondary eclipse observations covering the 1.1-1.6 micron wavelength range have revealed an atmospheric thermal inversion on the dayside hemisphere, likely caused by high altitude absorption at optical wavelengths. Here we present secondary eclipse observations made with the Hubble Space Telescope Wide Field Camera 3 spectrograph that extend the wavelength coverage from 1.1 micron down to 0.8 micron. To determine the atmospheric properties from the measured eclipse spectrum, we performed a retrieval analysis assuming chemical equilibrium, with the effects of thermal dissociation and ionization included. Our best-fit model provides a good fit to the data with reduced chi^2=1.04. The data diverge from a blackbody spectrum and instead exhibit emission due to H- shortward of 1.1 micron. The best-fit model does not reproduce a previously reported bump in the spectrum at 1.25 micron, possibly indicating this feature is a statistical fluctuation in the data rather than a VO emission band as had been tentatively suggested. We estimate an atmospheric metallicity of [M/H]=1.09(-0.69,+0.57), and fit for the carbon and oxygen abundances separately, obtaining [C/H]=-0.29(-0.48,+0.61) and [O/H]=0.18(-0.60,+0.64). The corresponding carbon-to-oxygen ratio is C/O=0.49(-0.37,+0.65), which encompasses the solar value of 0.54, but has a large uncertainty.

High-mass stars and star clusters commonly form within hub-filament systems. Monoceros R2, 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 hub-filament system on basis of 13CO and C18O observations obtained with the IRAM-30m telescope and H2 column density maps derived from Herschel dust emission observations. We identified the filamentary network 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 Msun/pc and show signs of fragmentation. The secondary filaments have line masses of 12-60 Msun/pc and show fragmentation only sporadically. In the context of Ostriker's hydrostatic filament model, the main filaments are thermally super-critical. If non-thermal motions are included, most of them are trans-critical. Most of the secondary filaments are roughly trans-critical regardless of whether non-thermal motions are included or not. From the main filaments, we estimate a mass accretion rate of 10(-4)-10(-3) Msun/pc into the hub. The secondary filaments accrete into the main filaments with a rate of 0.1-0.4x10(-4) Msun/pc. The main filaments extend into the 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.

Although the 2P3/2-2P1/2 transition of [CII] at 158um is known to be an excellent tracer of active star formation, we still do not have a complete understanding of where within star formation regions the emission originates. Here, we use SOFIA upGREAT observations of [CII] emission toward the HII region complex Sh2-235 (S235) to better understand in detail the origin of [CII] emission. We complement these data with a fully-sampled Green Bank Telescope radio recombination line map tracing the ionized hydrogen gas. About half of the total [CII] emission associated with S235 is spatially coincident with ionized hydrogen gas, although spectroscopic analysis shows little evidence that this emission is coming from the ionized hydrogen volume. Velocity-integrated [CII] intensity is strongly correlated with WISE 12um intensity across the entire complex, indicating that both trace ultra-violet radiation fields. The 22um and radio continuum intensities are only correlated with [CII] intensity in the ionized hydrogen portion of the S235 region and the correlations between the [CII] and molecular gas tracers are poor across the region. We find similar results for emission averaged over a sample of external galaxies, although the strength of the correlations is weaker. Therefore, although many tracers are correlated with the strength of [CII] emission, only WISE 12um emission is correlated on small-scales of the individual HII region S235 and also has a decent correlation at the scale of entire galaxies. Future studies of a larger sample of Galactic HII regions would help to determine whether these results are truly representative.

Convection is an important physical process in astrophysics well-studied using numerical simulations under the Boussinesq and/or anelastic approximations. However these approaches reach their limits when compressible effects are important in the high Mach flow regime, e.g. in stellar atmospheres or in the presence of accretion shocks. In order to tackle these issues, we propose a new high performance and portable code, called "ARK" with a numerical solver well-suited for the stratified compressible Navier-Stokes equations. We take a finite volume approach with machine precision conservation of mass, transverse momentum and total energy. Based on previous works in applied mathematics we propose the use of a low Mach correction to achieve a good precision in both low and high Mach regimes. The gravity source term is discretized using a well-balanced scheme in order to reach machine precision hydrostatic balance. This new solver is implemented using the Kokkos library in order to achieve high performance computing and portability across different architectures (e.g. multi-core, many-core, and GP-GPU). We show that the low-Mach correction allows to reach the low-Mach regime with a much better accuracy than a standard Godunov-type approach. The combined well-balanced property and the low-Mach correction allowed us to trigger Rayleigh-B\'enard convective modes close to the critical Rayleigh number. Furthermore we present 3D turbulent Rayleigh-B\'enard convection with low diffusion using the low-Mach correction leading to a higher kinetic energy power spectrum. These results are very promising for future studies of high Mach and highly stratified convective problems in astrophysics.

The L/T transition is an important evolutionary phase in brown dwarf atmospheres, providing us with a unique opportunity to explore the effects of clouds, convection, winds, gravity and metallicity across a very narrow temperature range. Understanding these physical processes is critical for understanding ultracool atmospheres. In the next decade, we will answer three key questions regarding L/T transition atmospheres: 1. What is the physical mechanism behind the L/T transition? 2. What is the spatial extent of atmospheric structures at the L/T transition? 3. How do gravity and metallicity affect the L/T transition? The theory and methods developed for brown dwarfs will be used in the 2030s and beyond for solar-system age giant exoplanets and eventually habitable zone earth analogues. Developing these techniques now are crucial.

Filaments are a key step on the path that leads from molecular clouds to star formation. However, their characteristics are heavily debated, and the exact processes that lead to their formation and fragmentation into dense cores still remain to be fully understood. We aim at characterising the mass, kinematics, and stability against gravitational collapse of a statistically significant sample of filaments in the Orion B molecular cloud, which is renown for its very low star formation efficiency. We characterise the gas column densities and kinematics over a field of 1.9 deg$^2$, using C$^{18}$O(J=1-0) data from the IRAM-30m large programme ORION-B. Using two different Hessian-based filters, we extract and compare two filamentary networks, each containing over 100 filaments. Independent of the extraction method, the filaments have widths of 0.12$\pm$0.04 pc, and show a wide range of linear (1 - 100 $M_{\odot}$pc$^{-1}$) and volume densities (2.10$^3$ - 2.10$^5$ cm$^{-3}$). Compared to previous studies, the filament population is dominated by low-density, thermally sub-critical structures, suggesting that most of the identified filaments are not collapsing to form stars. In fact, only ~1% of the Orion B cloud mass covered by our observations can be found in super-critical, star-forming filaments, explaining the low star formation efficiency of the region. The velocity profiles observed across the filaments show quiescence in the centre, and coherency in the plane of the sky, despite being mostly supersonic. The filaments in Orion B apparently belong to a continuum which contains a few elements comparable to already studied star-forming filaments as well as many lower-density, gravitationally unbound structures. This comprehensive study of the Orion B filaments shows that the mass fraction in super-critical filaments is a key factor in determining star formation efficiency.

In the original manuscript arXiv:1710.10269, we presented a grid of forward model transmission spectra for hot Jupiter exoplanets. However, we recently identified an error in the treatment of rainout in our 1D atmosphere model ATMO. We explain the error, the correction, its validation and the changed conclusions in this erratum. The correction of this error led to changes in the equilibrium chemical abundances using rainout condensation and thereby the transmission spectra. We note that this error only affects the online library that includes rainout condensation, the library with local condensation (without rainout) is unaffected. The online library was updated with the correction in July 2018.

Three dimensional models that account for chemistry are useful tools to predict the chemical composition of (exo)planet and brown dwarf atmospheres and interpret observations of future telescopes. Recent Juno observations of the NH3 tropospheric distribution in Jupiter also indicate that 3D chemical modelling may be necessary to constrain the deep composition of the giant planets of the Solar System. However, due to the high computational cost of chemistry calculations, 3D chemical modelling has so far been limited. A solution to include chemical kinetics in a 3D model is to use a reduced chemical scheme. In this view, we have developed a reduced scheme containing 30 species and 181 reversible reactions, from the full chemical scheme of Venot et al. 2012. This scheme reproduces accurately the vertical profiles of the observable species (H2O, CH4, CO, CO2, NH3, and HCN) and has a large range of validity. It can be used to study all kind of warm atmospheres, except hot C-rich ones, which contains a high amount of C2H2. It can be used in 1D models, for fast computations, but also in 3D models for hot giant (exo)planet and brown dwarf atmospheres.

By generalizing the theory of convection to any type of thermal and compositional source terms (diabatic processes), we show that thermohaline convection in Earth oceans, fingering convection in stellar atmospheres, and moist convection in Earth atmosphere are deriving from the same general diabatic convective instability. We show also that "radiative convection" triggered by CO/CH4 transition with radiative transfer in the atmospheres of brown dwarfs is analog 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 atmosphere thus possibly explaining the reddening in brown-dwarf spectra. By using idealized two-dimensional 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 could be 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 to use brown dwarfs to better understand some aspects of the physics at play in the climate of our own planet.

We present an atmospheric transmission spectrum for the ultra-hot Jupiter WASP-121b, measured using the Space Telescope Imaging Spectrograph (STIS) onboard the Hubble Space Telescope (HST). Across the 0.47-1 micron wavelength range, the data imply an atmospheric opacity comparable to - and in some spectroscopic channels exceeding - that previously measured at near-infrared wavelengths (1.15-1.65 micron). Wavelength-dependent variations in the opacity rule out a gray cloud deck at a confidence level of 3.8-sigma and may instead be explained by VO spectral bands. We find a cloud-free model assuming chemical equilibrium for a temperature of 1500K and metal enrichment of 10-30x solar matches these data well. Using a free-chemistry retrieval analysis, we estimate a VO abundance of -6.6(-0.3,+0.2) dex. We find no evidence for TiO and place a 3-sigma upper limit of -7.9 dex on its abundance, suggesting TiO may have condensed from the gas phase at the day-night limb. The opacity rises steeply at the shortest wavelengths, increasing by approximately five pressure scale heights from 0.47 to 0.3 micron in wavelength. If this feature is caused by Rayleigh scattering due to uniformly-distributed aerosols, it would imply an unphysically high temperature of 6810+/-1530K. One alternative explanation for the short-wavelength rise is absorption due to SH (mercapto radical), which has been predicted as an important product of non-equilibrium chemistry in hot Jupiter atmospheres. Irrespective of the identity of the NUV absorber, it likely captures a significant amount of incident stellar radiation at low pressures, thus playing a significant role in the overall energy budget, thermal structure, and circulation of the atmosphere.

In this paper we present three-dimensional atmospheric simulations of the hot Jupiter HD~189733b under two different scenarios: local chemical equilibrium and including advection of the chemistry by the resolved wind. Our model consistently couples the treatment of dynamics, radiative transfer and chemistry, completing the feedback cycle between these three important processes. The effect of wind--driven advection on the chemical composition is qualitatively similar to our previous results for the warmer atmosphere of HD~209458b, found using the same model. However, we find more significant alterations to both the thermal and dynamical structure for the cooler atmosphere of HD~189733b, with changes in both the temperature and wind velocities reaching $\sim10\%$. We also present the contribution function, diagnosed from our simulations, and show that wind--driven chemistry has a significant impact on its three--dimensional structure, particularly for regions where methane is an important absorber. Finally, we present emission phase curves from our simulations and show the significant effect of wind--driven chemistry on the thermal emission, particularly within the 3.6 \textmu m Spitzer/IRAC channel.

The James Webb Space Telescope (JWST) presents the opportunity to transform our understanding of planets and the origins of life by revealing the atmospheric compositions, structures, and dynamics of transiting exoplanets in unprecedented detail. However, the high-precision, time-series observations required for such investigations have unique technical challenges, and prior experience with other facilities indicates that there will be a steep learning curve when JWST becomes operational. In this paper we describe the science objectives and detailed plans of the Transiting Exoplanet Community Early Release Science (ERS) Program, which is a recently approved program for JWST observations early in Cycle 1. The goal of this project, for which the obtained data will have no exclusive access period, is to accelerate the acquisition and diffusion of technical expertise for transiting exoplanet observations with JWST, while also providing a compelling set of representative datasets that will enable immediate scientific breakthroughs. The Transiting Exoplanet Community ERS Program will exercise the time-series modes of all four JWST instruments that have been identified as the consensus highest priorities, observe the full suite of transiting planet characterization geometries (transits, eclipses, and phase curves), and target planets with host stars that span an illustrative range of brightnesses. The observations in this program were defined through an inclusive and transparent process that had participation from JWST instrument experts and international leaders in transiting exoplanet studies. Community engagement in the project will be centered on a two-phase Data Challenge that culminates with the delivery of planetary spectra, time-series instrument performance reports, and open-source data analysis toolkits in time to inform the agenda for Cycle 2 of the JWST mission.

The G-type star GJ504A is known to host a 3 to 35 MJup companion whose temperature, mass, and projected separation all contribute to make it a test case for the planet formation theories and for atmospheric models of giant planets and light brown dwarfs. We collected data from the CHARA interferometer, SOPHIE spectrograph, and VLT/SPHERE high contrast imager to revisit the properties of the system. We measure a radius of 1.35+/- 0.04Rsun for GJ504A which yields isochronal ages of 21+/-2Myr or 4.0+/-1.8Gyr for the system and line-of-sight stellar rotation axis inclination of $162.4_{-4.3}^{+3.8}$ degrees or $18.6_{-3.8}^{+4.3}$ degrees. We re-detect the companion in the Y2, Y3, J3, H2, and K1 dual band SPHERE images. The complete 1-4 $\mu$m SED shape of GJ504b is best reproduced by T8-T9.5 objects with intermediate ages ($\leq1.5$Gyr), and/or unusual dusty atmospheres and/or super-solar metallicities. All six atmospheric models used yield $\mathrm{T_{eff}=550 \pm 50}$K for GJ504b and point toward a low surface gravity (3.5-4.0 dex). The accuracy on the metallicity value is limited by model-to-model systematics. It is not degenerate with the C/O ratio. We derive $\mathrm{log\:L/L_{\odot}=-6.15\pm0.15}$ dex for the companion compatible with masses of $\mathrm{M=1.3^{+0.6}_{-0.3}M_{Jup}}$ and $\mathrm{M=23^{+10}_{-9} M_{Jup}}$ for the young and old age ranges, respectively. The semi-major axis (sma) is above 27.8 au and the eccentricity lower than 0.55. The posterior on GJ~504b's orbital inclination suggests a misalignment with GJ~504A rotation axis. We combine the radial velocity and multi-epoch imaging data to exclude additional objects (90\% prob.) more massive than 2.5 and 30 $\mathrm{M_{Jup}}$ with sma in the range 0.01-80 au for the young and old system ages, respectively. The companion is in the envelope of the population of planets synthetized with our core-accretion model.

To understand and compare the 3D atmospheric structure of HD 209458 b and HD 189733 b, focusing on the formation and distribution of cloud particles, as well as their feedback on the dynamics and thermal profile. We couple the 3D Met Office Unified Model (UM), including detailed treatments of atmospheric radiative transfer and dynamics, to a kinetic cloud formation scheme. The resulting model self--consistently solves for the formation of condensation seeds, surface growth and evaporation, gravitational settling and advection, cloud radiative feedback via absorption and, crucially, scattering. Fluxes directly obtained from the UM are used to produce synthetic SEDs and phase curves. Our simulations show extensive cloud formation in both planets, however, cooler temperatures in the HD 189733 b result in higher cloud particle number densities. Sub-micron particles are suspended by vertical flows leading to extensive upper-atmosphere cloud cover. A combination of meridional advection and efficient cloud formation in cooler high latitude regions, result in enhanced cloud coverage for latitudes > 30 degrees and leads to a zonally banded structure for all our simulations. The cloud bands extend around the entire planet(s), as the temperatures, even on the day side, remain below the condensation temperature of silicates and oxides. Therefore, our simulated optical phase curve for HD 209458 b shows no `offset', in contrast to observations. Efficient scattering by cloud results an atmospheric cooling of up to 250K, and an advection-driven fluctuating cloud opacity causes temporal variability in the thermal emission. The inclusion of this fundamental cloud-atmosphere radiative feedback leads to significant differences with approaches neglecting these physical elements and suggests both a note of caution of interpretations neglecting such cloud feedback and scattering, and merits further study.

We present a study of the effect of wind-driven advection on the chemical composition of hot Jupiter atmospheres using a fully-consistent 3D hydrodynamics, chemistry and radiative transfer code, the Met Office Unified Model (UM). Chemical modelling of exoplanet atmospheres has primarily been restricted to 1D models that cannot account for 3D dynamical processes. In this work we couple a chemical relaxation scheme to the UM to account for the chemical interconversion of methane and carbon monoxide. This is done consistently with the radiative transfer meaning that departures from chemical equilibrium are included in the heating rates (and emission) and hence complete the feedback between the dynamics, thermal structure and chemical composition. In this letter we simulate the well studied atmosphere of HD~209458b. We find that the combined effect of horizontal and vertical advection leads to an increase in the methane abundance by several orders of magnitude; directly opposite to the trend found in previous works. Our results demonstrate the need to include 3D effects when considering the chemistry of hot Jupiter atmospheres. We calculate transmission and emission spectra, as well as the emission phase curve, from our simulations. We conclude that gas-phase non-equilibrium chemistry is unlikely to explain the model-observation discrepancy in the 4.5\,{\textmu m} {\it Spitzer}/IRAC channel. However, we highlight other spectral regions, observable with the James Webb Space Telescope, where signatures of wind-driven chemistry are more prominant.

In this work we have performed a series of simulations of the atmosphere of GJ~1214b assuming different metallicities using the Met Office Unified Model (UM). The UM is a general circulation model (GCM) that solves the deep, non-hydrostatic equations of motion and uses a flexible and accurate radiative transfer scheme, based on the two-stream and correlated-$k$ approximations, to calculate the heating rates. In this work we consistently couple a well-tested Gibbs energy minimisation scheme to solve for the chemical equilibrium abundances locally in each grid cell for a general set of elemental abundances, further improving the flexibility and accuracy of the model. As the metallicity of the atmosphere is increased we find significant changes in the dynamical and thermal structure, with subsequent implications for the simulated phase curve. The trends that we find are qualitatively consistent with previous works, though with quantitative differences. We investigate in detail the effect of increasing the metallicity by splitting the mechanism into constituents, involving the mean molecular weight, the heat capacity and the opacities. We find the opacity effect to be the dominant mechanism in altering the circulation and thermal structure. This result highlights the importance of accurately computing the opacities and radiative transfer in 3D GCMs.

Given the forthcoming launch of the James Webb Space Telescope (JWST) which will allow observ- ing exoplanet atmospheres with unprecedented signal-over-noise ratio, spectral coverage and spatial resolution, the uncertainties in the atmosphere modelling used to interpret the data need to be as- sessed. As the first step, we compare three independent 1D radiative-convective models: ATMO, Exo-REM and petitCODE. We identify differences in physical and chemical processes taken into ac- count thanks to a benchmark protocol we developed. We study the impact of these differences on the analysis of observable spectra. We show the importance of selecting carefully relevant molecular linelists to compute the atmospheric opacity. Indeed, differences between spectra calculated with Hitran and ExoMol exceed the expected uncertainties of future JWST observations. We also show the limitation in the precision of the models due to uncertainties on alkali and molecule lineshape, which induce spectral effects also larger than the expected JWST uncertainties. We compare two chemical models, Exo-REM and Venot Chemical Code, which do not lead to significant differences in the emission or transmission spectra. We discuss the observational consequences of using equilibrium or out-of- equilibrium chemistry and the major impact of phosphine, detectable with the JWST.Each of the models has benefited from the benchmarking activity and has been updated. The protocol developed in this paper and the online results can constitute a test case for other models.

Previous attempts at segmenting molecular line maps of molecular clouds have focused on using position-position-velocity data cubes of a single line to separate the spatial components of the cloud. In contrast, wide field spectral imaging with large spectral bandwidth in the (sub)mm domain now allows to combine multiple molecular tracers to understand the different physical and chemical phases that constitute giant molecular clouds. We aim at using multiple tracers (sensitive to different physical processes) to segment a molecular cloud into physically/chemically similar regions (rather than spatially connected components). We use a machine learning clustering method (the Meanshift algorithm) to cluster pixels with similar molecular emission, ignoring spatial information. Simple radiative transfer models are used to interpret the astrophysical information uncovered by the clustering. A clustering analysis based only on the J=1-0 lines of 12CO, 13CO and C18O reveals distinct density/column density regimes (nH~100, 500, and >1000 cm-3), closely related to the usual definitions of diffuse, translucent and high-column-density regions. Adding two UV-sensitive tracers, the (1-0) lines of HCO+ and CN, allows us to distinguish two clearly distinct chemical regimes, characteristic of UV-illuminated and UV-shielded gas. The UV-illuminated regime shows overbright HCO+ and CN emission, which we relate to photochemical enrichment. We also find a tail of high CN/HCO+ intensity ratio in UV-illuminated regions. Finer distinctions in density classes (nH~7E3, and 4E4 cm-3) for the densest regions are also identified, likely related to the higher critical density of the CN and HCO+ (1-0) lines. The association of simultaneous multi-line, wide-field mapping and powerful machine learning methods such as the Meanshift algorithm reveals how to decode the complex information available in molecular tracers.

Atmospheric modeling of low-gravity (VL-G) young brown dwarfs remains a challenge. The presence of very thick clouds has been suggested 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 evolutionary models for these objects. We show that cloudless atmospheres assuming a temperature gradient reduction caused by fingering convection provides 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 temperature 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 constrains in the mid-infrared, helping to confirm/refute whether or not the NIR reddening is caused by fingering convection. We suggest that the presence/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 analogues.

Infrared radiation emitted from a planet contains information about the chemical composition and vertical temperature profile of its atmosphere. If upper layers are cooler than lower layers, molecular gases will produce absorption features in the planetary thermal spectrum. Conversely, if there is a stratosphere - where temperature increases with altitude - these molecular features will be observed in emission. It has been suggested that stratospheres could form in highly irradiated exoplanets, but the extent to which this occurs is unresolved both theoretically and observationally. A previous claim for the presence of a stratosphere remains open to question, owing to the challenges posed by the highly variable host star and the low spectral resolution of the measurements. Here we report a near-infrared thermal spectrum for the ultrahot gas giant WASP-121b, which has an equilibrium temperature of approximately 2,500 kelvin. Water is resolved in emission, providing a detection of an exoplanet stratosphere at 5-sigma confidence. These observations imply that a substantial fraction of incident stellar radiation is retained at high altitudes in the atmosphere, possibly by absorbing chemical species such as gaseous vanadium oxide and titanium oxide.

We present a comprehensive statistical analysis of star-forming objects located in the vicinities of 1 360 bubble structures throughout the Galactic Plane and their local environments. The compilation of ~70 000 star-forming sources, found in the proximity of the ionized (Hii) regions and detected in both Hi-GAL and GLIMPSE surveys, provided a broad overview of the different evolutionary stages of star-formation in bubbles, from prestellar objects to more evolved young stellar objects (YSOs). Surface density maps of star-forming objects clearly reveal an evolutionary trend where more evolved star-forming objects are found spatially located near the center, while younger star-forming objects are found at the edge of the bubbles. We derived dynamic ages for a subsample of 182 Hii regions for which kinematic distances and radio continuum flux measurements were available. We detect ~80% more star-forming sources per unit area in the direction of bubbles than in the surrounding fields. We estimate ~10% clump formation efficiency (CFE) of Hi-GAL clumps in bubbles, twice the CFE in fields not affected by feedback. We find higher CFE of protostellar clumps in younger bubbles, whose density of the bubble shells is higher. We argue that the formation rate from prestellar to protostellar phase is probably higher during the early stages of the bubble expansion. Evaluation of the fragmentation time inside the shell of bubbles advocates the preexistence of clumps in the medium before the bubble, as supported by numerical simulations. Approximately 23% of the Hi-GAL clumps are found located in the direction of a bubble, with 15% for prestellar clumps and 41% for protostellar clumps. We argue that the high fraction of protostellar clumps may be due to the acceleration of the star-formation process cause by the feedback of the (Hii) bubbles.

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 microns, we conducted a detailed atmospheric study of the transiting Neptune-mass exoplanet HAT-P-26b. We detected prominent H2O absorption bands with a maximum base-to-peak amplitude of 525ppm 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.

The anomalously large radii of strongly irradiated exoplanets have remained a major puzzle in astronomy. Based on a 2D steady state atmospheric circulation model, the validity of which is assessed by comparison to 3D calculations, we reveal a new mechanism, namely the advection of the potential temperature due to mass and longitudinal momentum conservation, a process occuring in the Earth's atmosphere or oceans. At depth, 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 HD209458b, but also 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 explaining the inflated radii of irradiated hot Jupiters.

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 on to 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 non-equilibrium atmospheric chemistry is important for these objects. Non-equilibrium cloud-free models reproduce well the near-infrared spectra and mid-infrared photometry for the warmer Y dwarfs with 425 <= T_eff K <= 450. A small amount of cloud cover may improve the model fits in the near-infrared for the Y dwarfs with 325 <= T_eff K <= 375. Neither cloudy nor cloud-free models reproduce the near-infrared photometry for the T_eff = 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 to 8.5 Gyr.

We have followed up the three Y0 dwarfs WISEPA J041022.71+150248.5, WISEPA J173835.53+273258.9 and WISEPC J205628.90+145953.3 using the UKIRT/WFCAM telescope/instruments. We find parallaxes that are more consistent and accurate than previously published values. We estimate absolute magnitudes in photometric pass-bands from $Y$ to $W3$ and find them to be consistent between the three Y0 dwarfs indicating the inherent cosmic absolute magnitude spread of these objects is small. We examine the MKO $J$ magnitudes over the four year time line and find small but significant monotonic variations. Finally we estimate physical parameters from a comparison of spectra and parallax to equilibrium and non-equilibrium models finding values consistent with solar metallicity, an effective temperature of 450-475\,K and log~g of 4.0-4.5.

Context. The combination of wideband receivers and spectrometers currently available in (sub-)millimeter observatories deliver wide- field hyperspectral imaging of the interstellar medium. Tens of spectral lines can be observed over degree wide fields in about fifty hours. This wealth of data calls for restating the physical questions about the interstellar medium in statistical terms. Aims. We aim at gaining information on the physical structure of the interstellar medium from a statistical analysis of many lines from different species over a large field of view, without requiring detailed radiative transfer or astrochemical modeling. Methods. We coupled a nonlinear rescaling of the data with one of the simplest multivariate analysis methods, namely the Principal Component Analysis, to decompose the observed signal into components that we interpret first qualitatively and then quantitatively based on our deep knowledge of the observed region and of the astrochemistry at play. Results. We identify 3 principal components, linear compositions of line brightness temperatures, that are correlated at various levels with the column density, the volume density and the UV radiation field. Conclusions. When sampling a sufficiently diverse mixture of physical parameters, it is possible to decompose the molecular emission in order to gain physical insight on the observed interstellar medium. This opens a new avenue for future studies of the interstellar medium.

The nature of turbulence in molecular clouds is one of the key parameters that control star formation efficiency: compressive motions, as opposed to solenoidal motions, can trigger the collapse of cores, or mark the expansion of Hii regions. We try to observationally derive the fractions of momentum density ($\rho v$) contained in the solenoidal and compressive modes of turbulence in the Orion B molecular cloud and relate these fractions to the star formation efficiency in the cloud. The implementation of a statistical method developed by Brunt & Federrath (2014), applied to a $^{13}$CO(J=1-0) datacube obtained with the IRAM-30m telescope, allows us to retrieve 3-dimensional quantities from the projected quantities provided by the observations, yielding an estimate of the compressive versus solenoidal ratio in various regions of the cloud. Despite the Orion B molecular cloud being highly supersonic (mean Mach number $\sim$ 6), the fractions of motion in each mode diverge significantly from equipartition. The cloud's motions are on average mostly solenoidal (excess > 8 % with respect to equipartition), which is consistent with its low star formation rate. On the other hand, the motions around the main star-forming regions (NGC 2023 and NGC 2024) prove to be strongly compressive. We have successfully applied to observational data a method that was so far only tested on simulations, and have shown that there can be a strong intra-cloud variability of the compressive and solenoidal fractions, these fractions being in turn related to the star formation efficiency. This opens a new possibility for star-formation diagnostics in galactic molecular clouds.

The correlated-k method is frequently used to speed up radiation calculations in both one-dimensional and three-dimensional atmosphere models. An inherent difficulty with this method is how to treat overlapping absorption, i.e. absorption by more than one gas in a given spectral region. We have evaluated the applicability of three different methods in hot Jupiter and brown dwarf atmosphere models, all of which have been previously applied within models in the literature: (i) Random overlap, both with and without resorting and rebinning, (ii) equivalent extinction and (iii) pre-mixing of opacities, where (i) and (ii) combine k-coefficients for different gases to obtain k-coefficients for a mixture of gases, while (iii) calculates k-coefficients for a given mixture from the corresponding mixed line-by-line opacities. We find that the random overlap method is the most accurate and flexible of these treatments, and is fast enough to be used in one-dimensional models with resorting and rebinning. In three-dimensional models such as GCMs it is too slow, however, and equivalent extinction can provide a speed-up of at least a factor of three with only a minor loss of accuracy while at the same time retaining the flexibility gained by combining k-coefficients computed for each gas individually. Pre-mixed opacities are significantly less flexible, and we also find that particular care must be taken when using this method in order to to adequately resolve steep variations in composition at important chemical equilibrium boundaries. We use the random overlap method with resorting and rebinning in our one-dimensional atmosphere model and equivalent extinction in our GCM, which allows us to e.g. consistently treat the feedback of non-equilibrium chemistry on the total opacity and therefore the calculated P-T profiles in our models.

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 a agreement between observations and both our day side emission and hot spot offset, however, our night side emission is too large. Overall our results are qualitatively similar to those found by Showman et al. (2009) 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.

Molecular cloud structure is regulated by stellar feedback in various forms. Two of the most important feedback processes are UV photoionisation and supernovae from massive stars. However, the precise response of the cloud to these processes, and the interaction between them, remains an open question. In particular, we wish to know under which conditions the cloud can be dispersed by feedback, which in turn can give us hints as to how feedback regulates the star formation inside the cloud. We perform a suite of radiative magnetohydrodynamic simulations of a 10^5 solar mass cloud with embedded sources of ionising radiation and supernovae, including multiple supernovae and a hypernova model. A UV source corresponding to 10% of the mass of the cloud is required to disperse the cloud, suggesting that the star formation efficiency should be on the order of 10%. A single supernova is unable to significantly affect the evolution of the cloud. However, energetic hypernovae and multiple supernovae are able to add significant quantities of momentum to the cloud, approximately 10^{43} g cm/s of momentum per 10^{51} ergs of supernova energy. This is on the lower range of estimates in other works, since dense gas clumps that remain embedded inside the HII region cause rapid cooling in the supernova blast. We argue that supernovae alone are unable to regulate star formation in molecular clouds, and that strong pre-supernova feedback is required to allow supernova blastwaves to propagate efficiently into the interstellar medium

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.

We present a new set of analytic models for the expansion of HII regions powered by UV photoionisation from massive stars and compare them to a new suite of radiative magnetohydrodynamic simulations of turbulent, self-gravitating molecular clouds. To perform these simulations we use the Eulerian adaptive mesh magnetohydrodynamics code RAMSES-RT, including radiative transfer of UV photons. Our analytic models successfully predict the global behaviour of the HII region provided the density and velocity structure of the cloud is known. We give estimates for the HII region behaviour based on a power law fit to the density field assuming that the system is virialised. We give a radius at which the ionisation front should stop expanding ("stall"). If this radius is smaller than the distance to the edge of the cloud, the HII region will be trapped by the cloud. This effect is more severe in collapsing clouds than in virialised clouds, since the density in the former increases dramatically over time, with much larger photon emission rates needed for the HII region to escape a collapsing cloud. We also measure the response of Jeans unstable gas to the HII regions to predict the impact of UV radiation on star formation in the cloud. We find that the mass in unstable gas can be explained by a model in which the clouds are evaporated by UV photons, suggesting that the net feedback on star formation should be negative

The radiative feedback of massive stars on molecular clouds creates pillars, globules and other features at the interface between the HII region and molecular cloud. We present here Herschel observations between 70 and 500 micron of the immediate environment of the Cygnus OB2 association, performed within the HOBYS program. All structures were detected based on their appearance at 70 micron, and have been classified as pillars, globules, evaporating gasous globules (EGGs), proplyd-like objects, and condensations. From the 70 and 160 micron flux maps, we derive the local FUV field on the PDR surfaces. In parallel, we use a census of the O-stars to estimate the overall FUV-field, that is 10^3-10^4 G_0 close to the central OB cluster (within 10 pc) and decreases down to a few tens G_0, in a distance of 50 pc. From a SED fit to the four longest Herschel wavelengths, we determine column density and temperature maps and derive masses, volume densities and surface densities for these structures. We find that the morphological classification corresponds to distinct physical properties. Pillars and globules have the longest estimated photoevaporation lifetimes, a few 10^6 yr, while all other features should survive less than that. These lifetimes are consistent with that found in simulations of turbulent, UV-illuminated clouds. We propose a tentative evolutionary scheme in which pillars can evolve into globules, which in turn then evolve into EGGs, condensations and proplyd-like objects.

We present new near-infrared spectra, obtained at Gemini Observatory, for two Y dwarfs: WISE J035000.32-565830.2 (W0350) and WISEP J173835.52+273258.9 (W1738). A FLAMINGOS-2 R=540 spectrum was obtained for W0350, covering 1.0 < lambda um < 1.7, and a cross-dispersed GNIRS R=2800 spectrum was obtained for W1738, covering 0.993-1.087 um, 1.191-1.305 um, 1.589-1.631 um, and 1.985-2.175 um, in four orders. We also present revised YJH photometry for W1738, using new NIRI Y and J imaging, and a re-analysis of the previously published NIRI H band images. We compare these data, together with previously published data for late-T and Y dwarfs, to cloud-free models of solar metallicity, calculated both in chemical equilibrium and with disequilibrium driven by vertical transport. We find that for the Y dwarfs the non-equilibrium models reproduce the near-infrared data better than the equilibrium models. The remaining discrepancies suggest that fine-tuning the CH_4/CO and NH_3/N_2 balance is needed. Improved trigonometric parallaxes would improve the analysis. Despite the uncertainties and discrepancies, the models reproduce the observed near-infrared spectra well. We find that for the Y0, W1738, T_eff = 425 +/- 25 K and log g = 4.0 +/- 0.25, and for the Y1, W0350, T_eff = 350 +/- 25 K and log g = 4.0 +/- 0.25. W1738 may be metal-rich. Based on evolutionary models, these temperatures and gravities correspond to a mass range for both Y dwarfs of 3-9 Jupiter masses, with W0350 being a cooler, slightly older, version of W1738; the age of W0350 is 0.3-3 Gyr, and the age of W1738 is 0.15-1 Gyr.

Column density (N) PDFs serve as a powerful tool to characterize the physical processes that influence the structure of molecular clouds. Star-forming clouds can best be characterized by lognormal PDFs for the lower N range and a power-law tail for higher N, commonly attributed to turbulence and self-gravity and/or pressure, respectively. We report here on PDFs obtained from observations of 12CO, 13CO, C18O, CS, and N2H+ in the Cygnus X North region and compare to a PDF derived from dust observations with the Herschel satellite. The PDF of 12CO is lognormal for Av~1-30, but is cut for higher Av due to optical depth effects. The PDFs of C18O and 13CO are mostly lognormal up for Av~1-15, followed by excess up to Av~40. Above that value, all CO PDFs drop, most likely due to depletion. The high density tracers CS and N2H+ exhibit only a power law distribution between Av~15 and 400, respectively. The PDF from dust is lognormal for Av~3-15 and has a power-law tail up to Av~500. Absolute values for the molecular line column densities are, however, rather uncertain due to abundance and excitation temperature variations. Taken the dust PDF face value, we 'calibrate' the molecular line PDF of CS to the one of the dust and determined an abundance [CS]/[H2] of 10^-9. The slopes of the power-law tails of the CS, N2H+, and dust PDFs are consistent with free-fall collapse of filaments and clumps. A quasi static configuration of filaments and clumps can possibly also account for the observed N-PDFs, as long as they have a sufficiently condensed density structure and external ram pressure by gas accretion is provided. The somehow flatter slopes of N2H+ and CS can reflect an abundance change and/or subthermal excitation at low column densities.

The admitted, conventional scenario to explain the complex spectral evolution of brown dwarfs (BD) since their first detections twenty years ago, has always been the key role played by micron-size condensates, called "dust" or "clouds", in their atmosphere. This scenario, however, faces major problems, in particular the J-band brightening and the resurgence of FeH absorption at the L to T transition, and a physical first-principle understanding of this transition is lacking. In this paper, we propose a new, completely different explanation for BD and extrasolar giant planet (EGP) spectral evolution, without the need to invoke clouds. We show that, due to the slowness of the CO/CH4 and N2/NH3 chemical reactions, brown dwarf (L and T, respectively) and EGP atmospheres are subject to a thermo-chemical instability similar in nature to the fingering or chemical convective instability present in Earth oceans and at the Earth core/mantle boundary. The induced small-scale turbulent energy transport reduces the temperature gradient in the atmosphere, explaining the observed increase in near infrared J-H and J-K colors of L dwarfs and hot EGPs, while a warming up of the deep atmosphere along the L to T transition, as the CO/CH4 instability vanishes, naturally solves the two aforementioned puzzles, and provides a physical explanation of the L to T transition. This new picture leads to a drastic revision of our understanding of BD and EGP atmospheres and their evolution.

The surroundings of HII regions can have a profound influence on their development, morphology, and evolution. This paper explores the effect of the environment on H II regions in the MonR2 molecular cloud. We aim to investigate the density structure of envelopes surrounding HII regions and to determine their collapse and ionisation expansion ages. The Mon R2 molecular cloud is an ideal target since it hosts an H II region association. Column density and temperature images derived from Herschel data were used together to model the structure of HII bubbles and their surrounding envelopes. The resulting observational constraints were used to follow the development of the Mon R2 ionised regions with analytical calculations and numerical simulations. The four hot bubbles associated with H II regions are surrounded by dense, cold, and neutral gas envelopes. The radial density profiles are reminiscent of those of low-mass protostellar envelopes. The inner parts of envelopes of all four HII regions could be free-falling because they display shallow density profiles. As for their outer parts, the two compact HII regions show a density profile, which is typical of the equilibrium structure of an isothermal sphere. In contrast, the central UCHii region shows a steeper outer profile, that could be interpreted as material being forced to collapse. The size of the heated bubbles, the spectral type of the irradiating stars, and the mean initial neutral gas density are used to estimate the ionisation expansion time, texp, 0.1Myr,for the dense UCHII and compact HII regions and 0.35 Myr for the extended HII region. The envelope transition radii between the shallow and steeper density profiles are used to estimate the time elapsed since the formation of the first proto stellar embryo, Tinf : 1Myr, for the ultra-compact, 1.5 / 3Myr for the compact, and greater than 6Myr for the extended HII regions.

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, NGC6334), obtained from dust emission observed with Herschel. The low column density range can be fit with a lognormal distribution. A first power-law tail starts above an extinction (Av) of ~6-14. It has a slope of alpha=1.3-2 for the rho~r^-alpha profile for an equivalent density distribution (spherical or cylindrical geometry), and is thus consistent with free-fall gravitational collapse. Above Av~40, 60, and 140, we detect an excess that can be fitted by a flatter power law tail with alpha>2. It correlates with the central regions of the cloud (ridges/hubs) of size ~1 pc and densities above 10^4 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 have studied the filaments extracted from the column density maps of the nearby Lupus 1, 3, and 4 molecular clouds, derived from photometric maps observed with the Herschel satellite. Filaments in the Lupus clouds have quite low column densities, with a median value of $\sim$1.5$\times$10$^{21}$ cm$^{-2}$ and most have masses per unit length lower than the maximum critical value for radial gravitational collapse. Indeed, no evidence of filament contraction has been seen in the gas kinematics. We find that some filaments, that on average are thermally subcritical, contain dense cores that may eventually form stars. This is an indication that in the low column density regime, the critical condition for the formation of stars may be reached only locally and this condition is not a global property of the filament. Finally, in Lupus we find multiple observational evidences of the key role that the magnetic field plays in forming filaments, and determining their confinement and dynamical evolution.

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 HII region. The aim of this work is to understand the differences observed between the twelve participating numerical codes against the various analytical expressions examining the D-type phase of HII region expansion. To do this, we propose two well-defined tests which are tackled by 1D and 3D grid- and SPH- 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 twelve 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 $\lesssim2\%$. 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 dataset for calibrating the treatment of ionizing radiation hydrodynamics codes.

We present an analysis of Spitzer/IRAC primary transit and secondary eclipse lightcurves measured for HD209458b, using Gaussian process models to marginalise over the intrapixel sensitivity variations in the 3.6 micron and 4.5 micron channels and the ramp effect in the 5.8 micron and 8.0 micron channels. The main advantage of this approach is that we can account for a broad range of degeneracies between the planet signal and systematics without actually having to specify a deterministic functional form for the latter. Our results do not confirm a previous claim of water absorption in transmission. Instead, our results are more consistent with a featureless transmission spectrum, possibly due to a cloud deck obscuring molecular absorption bands. For the emission data, our values are not consistent with the thermal inversion in the dayside atmosphere that was originally inferred from these data. Instead, we agree with another re-analysis of these same data, which concluded a non-inverted atmosphere provides a better fit. We find that a solar-abundance clear-atmosphere model without a thermal inversion underpredicts the measured emission in the 4.5 micron channel, which may suggest the atmosphere is depleted in carbon monoxide. An acceptable fit to the emission data can be achieved by assuming that the planet radiates as an isothermal blackbody with a temperature of $1484\pm 18$ K.

This work aims to improve the current understanding of the atmospheres of brown dwarfs, especially cold ones with spectral type 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 (1D) 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.

We analyse column density and temperature maps derived from Herschel dust continuum observations of a sample of massive infrared dark clouds (G11.11-0.12, G18.82-0.28, G28.37+0.07, G28.53-0.25). We disentangle the velocity structure of the clouds using 13CO 1-0 and 12CO 3-2 data, showing that these IRDCs are the densest regions in massive giant molecular clouds and not isolated features. The probability distribution function (PDF) of column densities for all clouds have a power-law distribution over all (high) column densities, regardless of the evolutionary stage of the cloud: G11.11-0.12, G18.82-0.28, and G28.37+0.07 contain (proto)-stars, while G28.53-0.25 shows no signs of star formation. This is in contrast to the purely log-normal PDFs reported for near/mid-IR extinction maps. We only find a log-normal distribution for lower column densities, if we perform PDFs of the column density maps of the whole GMC in which the IRDCs are embedded. By comparing the PDF slope and the radial column density profile, we attribute the power law to the effect of large-scale gravitational collapse and to local free-fall collapse of pre- and protostellar cores. Independent from the PDF analysis, we find infall signatures in the spectral profiles of 12CO for G28.37+0.07 and G11.11-0.12, supporting the scenario of gravitational collapse. IRDCs are the densest regions within GMCs, which may be the progenitors of massive stars or clusters. At least some of the IRDCs are probably the same features as ridges (high column density regions with N>1e23 cm-2 over small areas), which were defined for nearby IR-bright GMCs. Because IRDCs are only confined to the densest (gravity dominated) cloud regions, the PDF constructed from this kind of a clipped image does not represent the (turbulence dominated) low column density regime of the cloud.

Column-density maps of molecular clouds are one of the most important observables in the context of molecular cloud- and star-formation (SF) studies. With the Herschel satellite it is now possible to determine the column density from dust emission. We use observations and simulations to demonstrate how LOS contamination affects the column density probability distribution function (PDF). We apply a first-order approximation (removing a constant level) to the molecular clouds of Auriga, Maddalena, Carina and NGC3603. In perfect agreement with the simulations, we find that the PDFs become broader, the peak shifts to lower column densities, and the power-law tail of the PDF flattens after correction. All PDFs have a lognormal part for low column densities with a peak at Av~2, a deviation point (DP) from the lognormal at Av(DP)~4-5, and a power-law tail for higher column densities. Assuming a density distribution rho~r^-alpha, the slopes of the power-law tails correspond to alpha(PDF)=1.8, 1.75, and 2.5 for Auriga, Carina, and NGC3603 (alpha~1.5-2 is consistent gravitational collapse). We find that low-mass and high-mass SF clouds display differences in the overall column density structure. Massive clouds assemble more gas in smaller cloud volumes than low-mass SF ones. However, for both cloud types, the transition of the PDF from lognormal shape into power-law tail is found at the same column density (at Av~4-5 mag). Low-mass and high-mass SF clouds then have the same low column density distribution, most likely dominated by supersonic turbulence. At higher column densities, collapse and external pressure can form the power-law tail. The relative importance of the two processes can vary between clouds and thus lead to the observed differences in PDF and column density structure.

We use the Mopra radio telescope to test for expansion of the molecular gas associated with the bubble HII region RCW120. A ring, or bubble, morphology is common for Galactic HII regions, but the three-dimensional geometry of such objects is still unclear. Detected near- and far-side expansion of the associated molecular material would be consistent with a three-dimensional spherical object. We map the $J = 1\rightarrow 0$ transitions of $^{12}$CO, $^{13}$CO, C$^{18}$O, and C$^{17}$O, and detect emission from all isotopologues. We do not detect the $0_0\rightarrow 1_{-1} E$ masing lines of CH$_3$OH at 108.8939 GHz. The strongest CO emission is from the photodissociation region (PDR), and there is a deficit of emission toward the bubble interior. We find no evidence for expansion of the molecular material associated with RCW120 and therefore can make no claims about its geometry. The lack of detected expansion is roughly in agreement with models for the time-evolution of an HII region like RCW120, and is consistent with an expansion speed of $< 1.5\, {\rm km\, s^{-1}}$. Single-position CO spectra show signatures of expansion, which underscores the importance of mapped spectra for such work. Dust temperature enhancements outside the PDR of RCW120 coincide with a deficit of emission in CO, confirming that these temperature enhancements are due to holes in the RCW120 PDR. H$\alpha$ emission shows that RCW120 is leaking $\sim5\%$ of the ionizing photons into the interstellar medium (ISM) through PDR holes at the locations of the temperature enhancements. H-alpha emission also shows a diffuse "halo" from leaked photons not associated with discrete holes in the PDR. Overall $25\pm10\%$ of all ionizing photons are leaking into the nearby ISM.

This work aims at improving the current understanding of the interaction between H ii regions and turbulent molecular clouds. We propose a new method to determine the age of a large sample of OB associations by investigating the development of their associated H ii regions in the surrounding turbulent medium. Using analytical solutions, one-dimensional (1D), and three-dimensional (3D) simulations, we constrained the expansion of the ionized bubble depending on the turbulent level of the parent molecular cloud. A grid of 1D simulations was then computed in order to build isochrone curves for H ii regions in a pressure-size diagram. This grid of models allowed to date large sample of OB associations and was used on the H ii Region Discovery Survey (HRDS). Analytical solutions and numerical simulations showed that the expansion of H ii regions is slowed down by the turbulence up to the point where the pressure of the ionized gas is in a quasi-equilibrium with the turbulent ram pressure. Based on this result, we built a grid of 1D models of the expansion of H ii regions in a profile based on Larson laws. The 3D turbulence is taken into account by an effective 1D temperature profile. The ages estimated by the isochrones of this grid agree well with literature values of well-known regions such as Rosette, RCW 36, RCW 79, and M16. We thus propose that this method can be used to give ages of young OB associations through the Galaxy such as the HRDS survey and also in nearby extra-galactic sources.

The treatment of radiation transport in global circulation models (GCMs) is crucial to correctly describe Earth and exoplanet atmospheric dynamics processes. The two-stream approximation and correlated-k method are currently state-of-the-art approximations applied in both Earth and hot Jupiter GCM radiation schemes to facilitate rapid calculation of fluxes and heating rates. Their accuracy have been tested extensively for Earth-like conditions, but verification of the methods' applicability to hot Jupiter-like conditions is lacking in the literature. We are adapting the UK Met Office GCM, the Unified Model (UM), for the study of hot Jupiters, and present in this work the adaptation of the Edwards-Slingo radiation scheme based on the two-stream approximation and the correlated-k method. We discuss the calculation of absorption coefficients from high temperature line lists and highlight the large uncertainty in the pressure-broadened line widths. We compare fluxes and heating rates obtained with our adapted scheme to more accurate discrete ordinate (DO) line-by-line (LbL) calculations ignoring scattering effects. We find that, in most cases, errors stay below 10 % for both heating rates and fluxes using ~ 10 k-coefficients in each band and a diffusivity factor D = 1.66. The two-stream approximation and the correlated-k method both contribute non-negligibly to the total error. We also find that using band-averaged absorption coefficients, which have previously been used in radiative-hydrodynamical simulations of a hot Jupiter, may yield errors of ~ 100 %, and should thus be used with caution.

Ionization feedback should impact the probability distribution function (PDF) of the column density around the ionized gas. We aim to quantify this effect and discuss its potential link to the Core and Initial Mass Function (CMF/IMF). We used in a systematic way Herschel column density maps of several regions observed within the HOBYS key program: M16, the Rosette and Vela C molecular cloud, and the RCW 120 H ii region. We fitted the column density PDFs of all clouds with two lognormal distributions, since they present a double-peak or enlarged shape in the PDF. Our interpretation is that the lowest part of the column density distribution describes the turbulent molecular gas while the second peak corresponds to a compression zone induced by the expansion of the ionized gas into the turbulent molecular cloud. The condensations at the edge of the ionized gas have a steep compressed radial profile, sometimes recognizable in the flattening of the power-law tail. This could lead to an unambiguous criterion able to disentangle triggered from pre-existing star formation. In the context of the gravo-turbulent scenario for the origin of the CMF/IMF, the double peaked/enlarged shape of the PDF may impact the formation of objects at both the low-mass and the high-mass end of the CMF/IMF. In particular a broader PDF is required by the gravo-turbulent scenario to fit properly the IMF with a reasonable initial Mach number for the molecular cloud. Since other physical processes (e.g. the equation of state and the variations among the core properties) have already been suggested to broaden the PDF, the relative importance of the different effects remains an open question.

Pillars and globules are present in many high-mass star-forming regions, such as the Eagle nebula (M16) and the Rosette molecular cloud, and understanding their origin will help characterize triggered star formation. The formation mechanisms of these structures are still being debated. Recent numerical simulations have shown how pillars can arise from the collapse of the shell in on itself and how globules can be formed from the interplay of the turbulent molecular cloud and the ionization from massive stars. The goal here is to test this scenario through recent observations of two massive star-forming regions, M16 and Rosette. The column density structure of the interface between molecular clouds and H ii regions was characterized using column density maps obtained from far-infrared imaging of the Herschel HOBYS key programme. Then, the DisPerSe algorithm was used on these maps to detect the compressed layers around the ionized gas and pillars in different evolutionary states. Finally, their velocity structure was investigated using CO data, and all observational signatures were tested against some distinct diagnostics established from simulations. The column density profiles have revealed the importance of compression at the edge of the ionized gas. The velocity properties of the structures, i.e. pillars and globules, are very close to what we predict from the numerical simulations. We have identified a good candidate of a nascent pillar in the Rosette molecular cloud that presents the velocity pattern of the shell collapsing on itself, induced by a high local curvature. Globules have a bulk velocity dispersion that indicates the importance of the initial turbulence in their formation, as proposed from numerical simulations. Altogether, this study re-enforces the picture of pillar formation by shell collapse and globule formation by the ionization of highly turbulent clouds.

When transiting their host stars, hot Jupiters absorb about 10% of the light in the wings of the stellar Lyman-alpha emission line. The absorption occurs at wavelengths Doppler-shifted from line center by +/- 100 km/s - larger than the thermal speeds with which partially neutral, 10E4 K hydrogen escapes from hot Jupiter atmospheres. It has been proposed that the absorption arises from 10E6 K hydrogen from the host stellar wind, made momentarily neutral by charge exchange with planetary H I. The +/-100 km/s velocities would then be attributed to the typical velocity dispersions of protons in the stellar wind - as inferred from spacecraft measurements of the Solar wind. To test this proposal, we perform 2D hydrodynamic simulations of colliding hot Jupiter and stellar winds, augmented by a chemistry module to compute the amount of hot neutral hydrogen produced by charge exchange. The Kelvin-Helmholtz instability mixes the two winds; in the mixing layer, charge exchange reactions establish, within tens of seconds. Enough hot neutral hydrogen is generated to reproduce the transit observations, and the amount of absorption converges with both spatial resolution and time. Our calculations support the idea that charge transfer between colliding winds correctly explains Lyman-alpha transit observations - modulo the effects of magnetic fields, which we do not model but which may suppress mixing. Other neglected effects include, in order of decreasing importance, rotational forces related to orbital motion, gravity, and stellar radiation pressure; we discuss quantitatively the errors introduced by our approximations. How hot stellar hydrogen cools when it collides with cold planetary hydrogen is also considered; a more careful treatment of how the mixing layer thermally equilibrates might explain the recent detection of Balmer H-alpha absorption in transiting hot Jupiters.

The most important limitation for ground-based submillimetre (submm) astronomy is the broad-band absorption of the total water vapour in the atmosphere above an observation site, often expressed as the Precipitable Water Vapour (PWV). A long-term statistic on the PWV is thus mandatory to characterize the quality of an existing or potential site for observational submm-astronomy. In this study we present a three-year statistic (2008-2010) of the PWV for ground-based telescope sites all around the world and for stratospheric altitudes relevant for SOFIA (Stratospheric Observatory for Far-infrared astronomy). The submm-transmission is calculated for typical PWVs using an atmospheric model. We present the absolute PWV values for each site sorted by year and time percentage. The PWV corresponding to the first decile (10%) and the quartiles (25%, 50%, 75%) are calculated and transmission curves between 150 {\mu}m and 3 mm for these values are shown. The Antarctic and South-American sites present very good conditions for submillimetre astronomy. The 350 {\mu}m and 450 {\mu}m atmospheric windows are open all year long whereas the 200 {\mu}m atmospheric window opens reasonably for 25 % of the time in Antarctica and the extremely high-altitude sites in Chile. Potential interesting new facilities are Macon in Argentinia and Summit in Greenland that show similar conditions as for example Mauna Kea (Hawaii). For SOFIA, we present in more detail transmission curves for different altitudes (11 to 14 km), PWV values, and higher frequencies (up to 5 THz). Though the atmosphere at these altitude is generally very transparent, the absorption at very high frequencies becomes more important, partly caused by minor species. In conclusion, the method presented in this paper could identify sites on Earth with a great potential for submillimetre astronomy, and guide future site testing campaigns in situ.

We investigate the interplay between the ionization radiation from massive stars and the turbulence inside the surrounding molecular gas thanks to 3D numerical simulations. We used the 3D hydrodynamical code HERACLES to model an initial turbulent medium that is ionized and heated by an ionizing source. Three different simulations are performed with different mean Mach numbers (1, 2 and 4). A non-equilibrium model for the ionization and the associated thermal processes was used. This revealed to be crucial when turbulent ram pressure is of the same order as the ionized-gas pressure. The density structures initiated by the turbulence cause local curvatures of the dense shell formed by the ionization compression. When the curvature of the shell is sufficient, the shell collapse on itself to form a pillar while a smaller curvature leads to the formation of dense clumps that are accelerated with the shell and therefore remain in the shell during the simulation. When the turbulent ram pressure of the cold gas is sufficient to balance the ionized-gas pressure, some dense-gas bubbles have enough kinetic energy to penetrate inside the ionized medium, forming cometary globules. This suggests a direct relation in the observations between the presence of globules and the relative importance of the turbulence compared to the ionized-gas pressure. The probability density functions present a double peak structure when the turbulence is low relative to the ionized-gas pressure. This could be used in observations as an indication of the turbulence inside molecular clouds.

Molecular globules and pillars are spectacular features, found only in the interface region between a molecular cloud and an HII-region. Impacting Far-ultraviolet (FUV) radiation creates photon dominated regions (PDRs) on their surfaces that can be traced by typical cooling lines. With the GREAT receiver onboard SOFIA we mapped and spectrally resolved the [CII] 158 micron atomic fine-structure line and the highly excited 12CO J=11-10 molecular line from three objects in Cygnus X (a pillar, a globule, and a strong IRAS source). We focus here on the globule and compare our data with existing Spitzer data and recent Herschel Open-Time PACS data. Extended [CII] emission and more compact CO-emission was found in the globule. We ascribe this emission mainly to an internal PDR, created by a possibly embedded star-cluster with at least one early B-star. However, external PDR emission caused by the excitation by the Cyg OB2 association cannot be fully excluded. The velocity-resolved [CII] emission traces the emission of PDR surfaces, possible rotation of the globule, and high-velocity outflowing gas. The globule shows a velocity shift of ~2 km/s with respect to the expanding HII-region, which can be understood as the residual turbulence of the molecular cloud from which the globule arose. This scenario is compatible with recent numerical simulations that emphazise the effect of turbulence. It is remarkable that an isolated globule shows these strong dynamical features traced by the [CII]-line, but it demands more observational studies to verify if there is indeed an embedded cluster of B-stars.

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. The Rosette molecular cloud, irradiated by the NGC2244 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. 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. 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 HII--region and the molecular cloud.

Radiative feedback from massive stars is a key process to understand how HII regions may enhance or inhibit star formation in pillars and globules at the interface with molecular clouds. We aim to contribute to model the interactions between ionization and gas clouds to better understand the processes at work. We study in detail the impact of modulations on the cloud-HII region interface and density modulations inside the cloud. We run three-dimensional hydrodynamical simulations based on Euler equations coupled with gravity using the HERACLES code. We implement a method to solve ionization/recombination equations and we take into account typical heating and cooling processes at work in the interstellar medium and due to ionization/recombination physics. UV radiation creates a dense shell compressed between an ionization front and a shock ahead. Interface modulations produce a curved shock that collapses on itself leading to stable growing pillar-like structures. The narrower the initial interface modulation, the longer the resulting pillar. We interpret pillars resulting from density modulations in terms of the ability of these density modula- tions to curve the shock ahead the ionization front. The shock curvature is a key process to understand the formation of structures at the edge of HII regions. Interface and density modulations at the edge of the cloud have a direct impact on the morphology of the dense shell during its formation. Deeper in the cloud, structures have less influence due to the high densities reached by the shell during its expansion.

Over the past few years a major effort has been put into the exploration of potential sites for the deployment of submillimetre astronomical facilities. Amongst the most important sites are Dome C and Dome A on the Antarctic Plateau, and the Chajnantor area in Chile. In this context, we report on measurements of the sky opacity at 200 um over a period of three years at the French-Italian station, Concordia, at Dome C, Antarctica. We also present some solutions to the challenges of operating in the harsh polar environ- ment. Dome C offers exceptional conditions in terms of absolute atmospheric transmission and stability for submillimetre astron- omy. Over the austral winter the PWV exhibits long periods during which it is stable and at a very low level (0.1 to 0.3 mm). Higher values (0.2 to 0.8 mm) of PWV are observed during the short summer period. Based on observations over three years, a transmission of around 50% at 350 um is achieved for 75% of the time. The 200-um window opens with a typical transmission of 10% to 15% for 25% of the time. Dome C is one of the best accessible sites on Earth for submillimetre astronomy. Observations at 350 or 450 {\mu}m are possible all year round, and the 200-um window opens long enough and with a sufficient transparency to be useful. Although the polar environment severely constrains hardware design, a permanent observatory with appropriate technical capabilities is feasible. Because of the very good astronomical conditions, high angular resolution and time series (multi-year) observations at Dome C with a medium size single dish telescope would enable unique studies to be conducted, some of which are not otherwise feasible even from space.

PILOT (the Pathfinder for an International Large Optical Telescope) is a proposed 2.5 m optical/infrared telescope to be located at Dome C on the Antarctic plateau. Conditions at Dome C are known to be exceptional for astronomy. The seeing (above ~30 m height), coherence time, and isoplanatic angle are all twice s good as at typical mid-latitude sites, while the water-vapour column, and the atmosphere and telescope thermal emission are all an order of magnitude better. These conditions enable a unique scientific capability for PILOT, which is addressed in this series of papers. The current paper presents an overview of the optical and instrumentation suite for PILO and its expected performance, a summary of the key science goals and observational approach for the facility, a discussion of the synergies between the science goals for PILOT and other telescopes, and a discussion of the future of Antarctic astronomy. Paper II and Paper III present details of the science projects divided, respectively, between the distant Universe (i.e., studies of first light, and the assembly and evolution of structure) and the nearby Universe (i.e., studies of Local Group galaxies, the Milky Way, and the Solar System).