89 research outputs found
Stellar Feedback in the ISM Revealed by Wide-Field Far-Infrared Spectral-Imaging
The radiative and mechanical interaction of stars with their environment
drives the evolution of the ISM and of galaxies as a whole. The far-IR emission
(lambda ~30 to 350 microns) from atoms and molecules dominates the cooling of
the warm gas in the neutral ISM, the material that ultimately forms stars.
Far-IR lines are thus the most sensitive probes of stellar feedback processes,
and allow us to quantify the deposition and cycling of energy in the ISM. While
ALMA (in the (sub)mm) and JWST (in the IR) provide astonishing sub-arcsecond
resolution images of point sources and their immediate environment, they cannot
access the main interstellar gas coolants, nor are they designed to image
entire star-forming regions (SFRs). Herschel far-IR photometric images of the
interstellar dust thermal emission revealed the ubiquitous large-scale
filamentary structure of SFRs, their mass content, and the location of
thousands of prestellar cores and protostars. These images, however, provide a
static view of the ISM: not only they dont constrain the cloud dynamics,
moreover they cannot reveal the chemical composition and energy transfer within
the cloud, thus giving little insight into the regulation process of star
formation by stellar feedback. In this white paper we emphasize the need of a
space telescope with wide-field spectral-imaging capabilities in the critical
far-IR domain.Comment: White Paper submitted to the Astro 2020 Decadal Survey on Astronomy
and Astrophysics (National Academies of Science, Engineering, and Medicine
Efficient ortho-para conversion of H2 on interstellar grain surfaces
Context: Fast surface conversion between ortho- and para-H2 has been observed
in laboratory studies, and this mechanism has been proposed to play a role in
the control of the ortho-para ratio in the interstellar medium. Observations of
rotational lines of H2 in Photo-Dissociation Regions (PDRs) have indeed found
significantly lower ortho-para ratios than expected at equilibrium. The
mechanisms controlling the balance of the ortho-para ratio in the interstellar
medium thus remain incompletely understood, while this ratio can affect the
thermodynamical properties of the gas (equation of state, cooling function).
Aims: We aim to build an accurate model of ortho-para conversion on dust
surfaces based on the most recent experimental and theoretical results, and to
validate it by comparison to observations of H2 rotational lines in PDRs.
Methods: We propose a statistical model of ortho-para conversion on dust grains
with fluctuating dust temperatures, based on a master equation approach. This
computation is then coupled to full PDR models and compared to PDR
observations. Results: We show that the observations of rotational H2 lines
indicate a high conversion efficiency on dust grains, and that this high
efficiency can be accounted for if taking dust temperature fluctuations into
account with our statistical model of surface conversion. Simpler models
neglecting the dust temperature fluctuations do not reach the high efficiency
deduced from the observations. Moreover, this high efficiency induced by dust
temperature fluctuations is quite insensitive to the values of microphysical
parameters of the model. Conclusions: Ortho-para conversion on grains is thus
an efficient mechanism in most astrophysical conditions that can play a
significant role in controlling the ortho-para ratio.Comment: Accepted in Astronomy & Astrophysic
Surface chemistry in the Interstellar Medium II. formation on dust with random temperature fluctuations
The formation on grains is known to be sensitive to dust
temperature, which is also known to fluctuate for small grain sizes due to
photon absorption. We aim at exploring the consequences of simultaneous
fluctuations of the dust temperature and the adsorbed H-atom population on the
formation rate under the full range of astrophysically relevant
UV intensities and gas conditions. The master equation approach is generalized
to coupled fluctuations in both the grain's temperature and its surface
population and solved numerically. The resolution can be simplified in the case
of the Eley-Rideal mechanism, allowing a fast computation. For the
Langmuir-Hinshelwood mechanism, it remains computationally expensive, and
accurate approximations are constructed. We find the Langmuir-Hinshelwood
mechanism to become an efficient formation mechanism in unshielded photon
dominated region (PDR) edge conditions when taking those fluctuations into
account, despite hot average dust temperatures. It reaches an importance
comparable to the Eley-Rideal mechanism. However, we show that a simpler rate
equation treatment gives qualitatively correct observable results in full cloud
simulations under most astrophysically relevant conditions. Typical differences
are a factor of 2-3 on the intensities of the lines. We
also find that rare fluctuations in cloud cores are sufficient to significantly
reduce the formation efficiency. Our detailed analysis confirms that the usual
approximations used in numerical models are adequate when interpreting
observations, but a more sophisticated statistical analysis is required if one
is interested in the details of surface processes.Comment: 21 pages, 28 figures, accepted in A&
Efficient sampling of non log-concave posterior distributions with mixture of noises
This paper focuses on a challenging class of inverse problems that is often
encountered in applications. The forward model is a complex non-linear
black-box, potentially non-injective, whose outputs cover multiple decades in
amplitude. Observations are supposed to be simultaneously damaged by additive
and multiplicative noises and censorship. As needed in many applications, the
aim of this work is to provide uncertainty quantification on top of parameter
estimates. The resulting log-likelihood is intractable and potentially
non-log-concave. An adapted Bayesian approach is proposed to provide
credibility intervals along with point estimates. An MCMC algorithm is proposed
to deal with the multimodal posterior distribution, even in a situation where
there is no global Lipschitz constant (or it is very large). It combines two
kernels, namely an improved version of (Preconditioned Metropolis Adjusted
Langevin) PMALA and a Multiple Try Metropolis (MTM) kernel. Whenever smooth,
its gradient admits a Lipschitz constant too large to be exploited in the
inference process. This sampler addresses all the challenges induced by the
complex form of the likelihood. The proposed method is illustrated on classical
test multimodal distributions as well as on a challenging and realistic inverse
problem in astronomy
High angular resolution near-IR view of the Orion Bar revealed by Keck/NIRC2
Nearby Photo-Dissociation Regions (PDRs), where the gas and dust are heated
by the far UV-irradiation emitted from stars, are ideal templates to study the
main stellar feedback processes. With this study we aim to probe the detailed
structures at the interfaces between ionized, atomic, and molecular gas in the
Orion Bar. This nearby prototypical strongly irradiated PDR will be among the
first targets of the James Webb Space Telescope (JWST) within the framework of
the PDRs4All Early Release Science program. We employed the sub-arcsec
resolution accessible with Keck-II NIRC2 and its adaptive optics system to
obtain the most detailed and complete images, ever performed, of the
vibrationally excited line H 1-0 S(1) at 2.12~m, tracing the
dissociation front, and the [FeII] and Br lines, at 1.64 and
2.16~m respectively, tracing the ionization front. We obtained narrow-band
filter images in these key gas line diagnostic over at spatial
scales of 0.1 (0.0002~pc or 40~AU at 414~pc). The
Keck/NIRC2 observations spatially resolve a plethora of irradiated
sub-structures such as ridges, filaments, globules and proplyds. A remarkable
spatial coincidence between the H 1-0 S(1) vibrational and HCO J=4-3
rotational emission previously obtained with ALMA is observed. This likely
indicates the intimate link between these two molecular species and highlights
that in high pressure PDR the H/H and C/C/CO transitions zones come
closer as compared to a typical layered structure of a constant density PDR.
This is in agreement with several previous studies that claimed that the Orion
Bar edge is composed of very small, dense, highly irradiated PDRs at high
thermal pressure immersed in a more diffuse environment
H2 formation on interstellar dust grains: the viewpoints of theory, experiments, models and observations
Molecular hydrogen is the most abundant molecule in the universe. It is the first one to form and survive photo-dissociation in tenuous environments. Its formation involves catalytic reactions on the surface of interstellar grains. The micro-physics of the formation process has been investigated intensively in the last 20 years, in parallel of new astrophysical observational and modeling progresses. In the perspectives of the probable revolution brought by the future satellite JWST, this article has been written to present what we think we know about the H formation in a variety of interstellar environments.VW’s research is funded by an ERC Starting Grant (3DICE, grant agreement 336474). GV acknowledges financial support from the National Science Foundation’s Astronomy & Astrophysics Division (Grants No. 1311958 and 1615897). LH acknowledges support from ERC Consolidator Grant GRANN (grant agreement no. 648551). GN acknowledges support from the Swedish Research Council. VW, FD and SM acknowledge the CNRS program ”Physique et Chimie du Milieu Interstellaire” (PCMI) co-funded bythe Centre National d’Etudes Spatiales (CNES). SDP acknowledges funding from STFC, UK. V.V acknowledges funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (MagneticYSOS project, grant agreement No 679937)
Radiative and mechanical feedback into the molecular gas in the Large Magellanic Cloud. I. N159W
We present Herschel SPIRE Fourier Transform Spectrometer (FTS) observations
of N159W, an active star-forming region in the Large Magellanic Cloud (LMC). In
our observations, a number of far-infrared cooling lines including CO(4-3) to
CO(12-11), [CI] 609 and 370 micron, and [NII] 205 micron are clearly detected.
With an aim of investigating the physical conditions and excitation processes
of molecular gas, we first construct CO spectral line energy distributions
(SLEDs) on 10 pc scales by combining the FTS CO transitions with ground-based
low-J CO data and analyze the observed CO SLEDs using non-LTE radiative
transfer models. We find that the CO-traced molecular gas in N159W is warm
(kinetic temperature of 153-754 K) and moderately dense (H2 number density of
(1.1-4.5)e3 cm-3). To assess the impact of the energetic processes in the
interstellar medium on the physical conditions of the CO-emitting gas, we then
compare the observed CO line intensities with the models of photodissociation
regions (PDRs) and shocks. We first constrain the properties of PDRs by
modelling Herschel observations of [OI] 145, [CII] 158, and [CI] 370 micron
fine-structure lines and find that the constrained PDR components emit very
weak CO emission. X-rays and cosmic-rays are also found to provide a negligible
contribution to the CO emission, essentially ruling out ionizing sources
(ultraviolet photons, X-rays, and cosmic-rays) as the dominant heating source
for CO in N159W. On the other hand, mechanical heating by low-velocity C-type
shocks with ~10 km/s appears sufficient enough to reproduce the observed warm
CO.Comment: accepted for publication in A&
Neural network-based emulation of interstellar medium models
The interpretation of observations of atomic and molecular tracers in the
galactic and extragalactic interstellar medium (ISM) requires comparisons with
state-of-the-art astrophysical models to infer some physical conditions.
Usually, ISM models are too time-consuming for such inference procedures, as
they call for numerous model evaluations. As a result, they are often replaced
by an interpolation of a grid of precomputed models.
We propose a new general method to derive faster, lighter, and more accurate
approximations of the model from a grid of precomputed models.
These emulators are defined with artificial neural networks (ANNs) designed
and trained to address the specificities inherent in ISM models. Indeed, such
models often predict many observables (e.g., line intensities) from just a few
input physical parameters and can yield outliers due to numerical instabilities
or physical bistabilities. We propose applying five strategies to address these
characteristics: 1) an outlier removal procedure; 2) a clustering method that
yields homogeneous subsets of lines that are simpler to predict with different
ANNs; 3) a dimension reduction technique that enables to adequately size the
network architecture; 4) the physical inputs are augmented with a polynomial
transform to ease the learning of nonlinearities; and 5) a dense architecture
to ease the learning of simple relations.
We compare the proposed ANNs with standard classes of interpolation methods
to emulate the Meudon PDR code, a representative ISM numerical model.
Combinations of the proposed strategies outperform all interpolation methods by
a factor of 2 on the average error, reaching 4.5% on the Meudon PDR code. These
networks are also 1000 times faster than accurate interpolation methods and
require ten to forty times less memory.
This work will enable efficient inferences on wide-field multiline
observations of the ISM
- …