Source author record

P. Tremblin

P. Tremblin appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

27works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

27 published item(s)

preprint2025arXiv

JWST/MIRI observations of the young TWA 27 system: hydrocarbon disk chemistry, silicate clouds, evidence for a CPD

The Mid-Infrared Instrument (MIRI) on the James Webb Space Telescope (JWST) enables the characterisation of young self-luminous gas giants at previously inaccessible wavelengths, revealing physical processes in gas, dust, and clouds. We characterise the young planetary system TWA 27 (2M1207) in the mid-infrared (MIR), studying the atmosphere and disk spectra of the M9 brown dwarf TWA 27A and its L6 planetary-mass companion TWA 27b. We obtained data with the MIRI Medium Resolution Spectrometer (MRS) from 4.9 to 20 um, and MIRI Imaging in the F1000W and F1500W filters. We applied high-contrast imaging methods to extract the companion's spectral energy distribution up to 15 um at 0.78 arcsec separation and a contrast of 60. Combining these with published JWST/NIRSpec spectra, we analysed the 1-20 um range using self-consistent atmospheric grids and 0D slab models for molecular disk emission. The atmosphere of TWA 27A is well fitted by a BT-SETTL model with Teff 2780 K, log g 4.3, plus a 740 K blackbody for the inner disk rim. The disk shows at least 11 organic molecules, with no water or silicate dust emission detected. The atmosphere of TWA 27b is matched by a Teff 1400 K low-gravity model with extinction, best fit by the ExoREM grid. MIRI spectra and photometry for TWA 27b reveal a silicate cloud absorption feature between 8-10 um and significant (>5 sigma) infrared excess at 15 um consistent with circumplanetary disk emission. These MIRI observations provide new insights into TWA 27, revealing diverse features to study the formation and evolution of circumplanetary disks and young dusty atmospheres.

preprint2025arXiv

Multi-modal atmospheric characterization of $β$ Pictoris b: Adding high-resolution continuum spectra from GRAVITY

We present the first VLTI/GRAVITY observations at R$_λ\sim 4000$ of $β$ Pic b. These four high S/N ($\sim$20) K-band spectra conserve both the pseudo-continuum and molecular absorption patterns. We analyze them with four self-consistent forward model grids (Exo-REM, ATMO, BT-Settl, Sonora) exploring $T_{\mathrm{eff}}$, log(g), metallicity, C/O, and $^{12}$CO/$^{13}$CO ratio. We also upgrade our forward modeling code ForMoSA to account for the data multi-modality and combine the GRAVITY epochs with published 1-5 $μ$m photometry, low- to medium-resolution spectra (0.9-7 $μ$m), and high-resolution echelle spectra (2.1-5.2 $μ$m). Sonora and Exo-REM are statistically preferred. Exo-REM yields $T_{\mathrm{eff}}$ $=1607.45^{+4.85}_{-6.20}$ K and log(g) $=4.46^{+0.02}_{-0.04}$ dex from GRAVITY alone, and $T_{\mathrm{eff}}$ $=1502.74^{+2.32}_{-2.14}$ K and log(g) $=4.00\pm0.01$ dex when including all datasets. Archival data significantly affect the retrieved parameters. C/O remains solar ($0.552^{+0.003}_{-0.002}$) while [M/H] reaches super-solar values (0.50$\pm$0.01). We report the first tentative constraint on log($^{12}$CO/$^{13}$CO) $\sim$1.12, though this remains inconclusive due to telluric residuals. Additionally, we estimate the luminosity to be log(L/L$_\odot$) $=-4.01^{+0.04}_{-0.05}$, implying a heavy-element content of up to $\sim$5% (20-80 M$_\oplus$) given the system age and dynamical mass measurements. Access to both continuum and molecular lines at K-band significantly impacts the metallicity, possibly owing to collision-induced absorption shaping the continuum. Echelle spectra do not dominate the final fit with respect to lower resolution data. Future multi-modal frameworks should include weighting schemes reflecting bandwidth and central wavelength coverage.

preprint2022arXiv

Nonideal self-gravity and cosmology: the importance of correlations in the dynamics of the large-scale structures of the Universe

Inspired by the role of correlations in the statistical mechanics of nonideal self-interacting fluids, we suggest that unresolved sub-structures (i.e. correlations) have to be taken into account in the Virial theorem of self-gravitating astrophysical systems. We demonstrate that their omission leads to a missing mass problem by using the semi-analytic polytropic solutions of the Lane-Emden equation. This problem suggests to extend the Friedmann equations to the nonideal regime by taking into account correlations in the dynamics of the expansion. The increase of correlations induced by the formation of the large-scale structures could explain naturally the accelerated expansion of the Universe in such a paradigm.

preprint2021arXiv

Non-ideal self-gravity and cosmology: the importance of correlations in the dynamics of the large-scale structures of the Universe

Inspired by the statistical mechanics of an ensemble of interacting particles (BBGKY hierarchy), we propose to account for small-scale inhomogeneities in self-gravitating astrophysical fluids by deriving a non-ideal Virial theorem and non-ideal Navier-Stokes equations using a decomposition of the gravitational force into a near- and far-field component. These equations involve the pair radial distribution function (similar to the two-point correlation function), similarly to the interaction energy and equation of state in liquids. Small-scale correlations lead to a non-ideal amplification of the gravitational interaction energy, whose omission leads to a missing mass problem, e.g., in galaxies and galaxy clusters. We also propose an extension of the Friedmann equations in the non-ideal regime. We estimate the non-ideal amplification factor of the gravitational interaction energy of the baryons to lie between 5 and 20, potentially explaining the observed value of the Hubble parameter. Within this framework, the acceleration of the expansion emerges naturally because of the increasing number of sub-structures induced by gravitational collapse, which increases their contribution to the total gravitational energy. A simple estimate predicts a non-ideal deceleration parameter qni~-1;this is potentially the first determination of the observed value based on an intuitively physical argument. We suggest that correlations and gravitational interactions could produce a transition to a viscous regime that can lead to flat rotation curves. This transition could also explain the dichotomy between (Keplerian) LSB elliptical galaxy and (non-Keplerian) spiral galaxy rotation profiles. Overall, our results demonstrate that non-ideal effects induced by inhomogeneities must be taken into account in order to properly determine the gravitational dynamics of galaxies and the larger scale universe.

preprint2020arXiv

A high-performance and portable asymptotic preserving radiation hydrodynamics code with the M1 model

Aims. We present a new radiation hydrodynamics code, called "ARK-RT" which uses a two-moment model with the M1 closure relation for radiative transfer. This code aims at being ready for high-performance computing, on exascale architectures. Methods. The two-moment model is solved using a finite volume scheme. The scheme is asymptotic preserving to capture accurately both optically thick and thin regimes. We also propose a well-balanced discretization of the radiative flux source term able to capture constant flux steady states with discontinuities in opacity. We use the library Trilinos for linear algebra and the package Kokkos allows us to reach high-performance computing and portability across different architectures, such as multi-core, many-core, and GP-GPU. Results. ARK-RT is able to reproduce standard tests in both free-streaming and diffusive limits, including purely radiative tests and radiation hydrodynamics ones. Using a time-implicit solver is profitable as soon as the time step given by the hydrodynamics is 50-100 times larger than the explicit time step for radiative transfer, depending on the preconditioner and the architecture. Albeit more work is needed to ensure stability in all circumstances. Using ARK-RT, we study the propagation of an ionization front in convective dense cores. We show that the ionization front is strongly stable against perturbations even with destabilizing convective motions. As a result, the presence of instabilities should be interpreted with caution. Overall, ARK-RT is well-suited to study many astrophysical problems involving convection and radiative transfer such as the dynamics of H ii regions in massive pre-stellar dense cores and future applications could include planetary atmospheres.

preprint2016arXiv

Cloudless atmospheres for L/T dwarfs and extra-solar giant planets

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.

preprint2016arXiv

Globules and Pillars in Cygnus X I. Herschel Far-infrared imaging of the Cyg OB2 environment

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.

preprint2016arXiv

Near-Infrared Spectroscopy of the Y0 WISEP J173835.52+273258.9 and the Y1 WISE J035000.32-565830.2: the Importance of Non-Equilibrium Chemistry

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.

preprint2016arXiv

Understanding star formation in molecular clouds III. Probability distribution functions of molecular lines in Cygnus X

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.

preprint2015arXiv

Detection of two power-law tails in the probability distribution functions of massive GMCs

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.

preprint2015arXiv

Filaments in the Lupus molecular clouds

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.

preprint2015arXiv

Fingering convection and cloudless models for cool brown dwarf atmospheres

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.

preprint2015arXiv

From forced collapse to H ii region expansion in Mon R2: Envelope density structure and age determination with Herschel

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.

preprint2015arXiv

Mopra CO Observations of the Bubble HII Region RCW120

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$α$ 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.

preprint2015arXiv

StarBench: The D-type expansion of an HII region

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.

preprint2015arXiv

Understanding star formation in molecular clouds I. Effects of line-of-sight contamination on the column density structure

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.

preprint2015arXiv

Understanding star formation in molecular clouds II. Signatures of gravitational collapse of IRDCs

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.

preprint2014arXiv

Age, size, and position of H ii regions in the Galaxy. Expansion of ionized gas in turbulent molecular clouds

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.

preprint2014arXiv

Ionization compression impact on dense gas distribution and star formation, Probability density functions around H ii regions as seen by Herschel

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.

preprint2013arXiv

Pillars and globules at the edges of H ii regions, Confronting Herschel observations and numerical simulations

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.

preprint2012arXiv

3D simulations of globules and pillars formation around HII regions: turbulence and shock curvature

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.

preprint2012arXiv

Cluster-formation in the Rosette molecular cloud at the junctions of filaments

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.

preprint2012arXiv

Globules and pillars seen in the [CII] 158 micron line with SOFIA

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.

preprint2012arXiv

Worldwide site comparison for submillimetre astronomy

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 μm and 3 mm for these values are shown. The Antarctic and South-American sites present very good conditions for submillimetre astronomy. The 350 μm and 450 μm atmospheric windows are open all year long whereas the 200 μ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.

preprint2011arXiv

3D simulations of pillars formation around HII regions: the importance of shock curvature

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.

preprint2011arXiv

Site testing for submillimetre astronomy at Dome C, Antarctica

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 μ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.

preprint2009arXiv

The Science Case for PILOT I: Summary and Overview

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).