Search

(Abridged) JWST continues to deliver incredibly detailed infrared (IR) images of star forming regions in the Milky Way and beyond. IR emission from star-forming regions is very spectrally rich due to emission from gas-phase atoms, ions, and polycyclic aromatic hydrocarbons (PAHs). Physically interpreting IR images of these regions relies on assumptions about the underlying spectral energy distribution in the imaging bandpasses. We aim to provide empirical prescriptions linking line, PAH, and continuum intensities from JWST images, to facilitate the interpretation of JWST images in a wide variety of contexts. We use JWST PDRs4All Near-Infrared Camera (NIRCam) and Mid-Infrared Instrument (MIRI) imaging and Near-Infrared Spectrograph (NIRSpec) integral field unit (IFU) and MIRI Medium Resolution Spectrograph (MRS) spectroscopic observations of the Orion Bar photodissociation region (PDR), to directly compare and cross-calibrate imaging and IFU data at ~100 AU resolution over a region where the radiation field and ISM environment evolves from the hot ionized gas to the cold molecular gas. We measure the relative contributions of line, PAH, and continuum emission to the NIRCam and MIRI filters as functions of local physical conditions. We provide empirical prescriptions based on NIRCam and MIRI images to derive intensities of emission lines and PAH features. Within the range of the environments probed in this study, these prescriptions accurately predict Pa-alpha, Br-alpha, PAH 3.3 um and 11.2 um intensities, while those for FeII 1.644 um, H_2 1--0 S(1) 2.12 um and 1--0 S(9) 4.96 um, and PAH 7.7 um show more complicated environmental dependencies. Linear combinations of JWST NIRCam and MIRI images provide effective tracers of ionized gas, H_2, and PAH emission in PDRs. We expect these recipes to be useful for both the Galactic and extragalactic communities., Comment: 31 pages, 24 figures, submitted to A&A

Infrared emission features toward interstellar gas of the IC 348 star cluster in Perseus have been recently proposed to originate from the amino acid tryptophan. The assignment was based on laboratory infrared spectra of tryptophan pressed into pellets, a method which is known to cause large frequency shifts compared to the gas phase. We assess the validity of the assignment based on the original Spitzer data as well as new data from JWST. In addition, we report new spectra of tryptophan condensed in para-hydrogen matrices to compare with the observed spectra. The JWST MIRI data do not show evidence for tryptophan, despite deeper integration toward IC 348. In addition, we show that several of the lines attributed to tryptophan are likely due to instrumental artifacts. This, combined with the new laboratory data, allows us to conclude that there is no compelling evidence for the tryptophan assignment.

One of the main problems in astrochemistry is determining the amount of sulfur in volatiles and refractories in the interstellar medium. The detection of the main sulfur reservoirs (icy H$_2$S and atomic gas) has been challenging, and estimates are based on the reliability of models to account for the abundances of species containing less than 1% of the total sulfur. The high sensitivity of the James Webb Space Telescope provides an unprecedented opportunity to estimate the sulfur abundance through the observation of the [S I] 25.249 $\mu$m line. We used the [S III] 18.7 $\mu$m, [S IV] 10.5 $\mu$m, and [S l] 25.249 $\mu$m lines to estimate the amount of sulfur in the ionized and molecular gas along the Orion Bar. For the theoretical part, we used an upgraded version of the Meudon photodissociation region (PDR) code to model the observations. New inelastic collision rates of neutral atomic sulfur with ortho- and para- molecular hydrogen were calculated to predict the line intensities. The [S III] 18.7 $\mu$m and [S IV] 10.5 $\mu$m lines are detected over the imaged region with a shallow increase (by a factor of 4) toward the HII region. We estimate a moderate sulfur depletion, by a factor of $\sim$2, in the ionized gas. The corrugated interface between the molecular and atomic phases gives rise to several edge-on dissociation fronts we refer to as DF1, DF2, and DF3. The [S l] 25.249 $\mu$m line is only detected toward DF2 and DF3, the dissociation fronts located farthest from the HII region. The detailed modeling of DF3 using the Meudon PDR code shows that the emission of the [S l] 25.249 $\mu$m line is coming from warm ($>$ 40 K) molecular gas located at A$_{\rm V}$ $\sim$ 1$-$5 mag from the ionization front. Moreover, the intensity of the [S l] 25.249 $\mu$m line is only accounted for if we assume the presence of undepleted sulfur., Comment: 16 pages, 6 figures. Accepted for publication in Astronomy and Astrophysics

Mid-infrared emission features probe the properties of ionized gas, and hot or warm molecular gas. The Orion Bar is a frequently studied photodissociation region (PDR) containing large amounts of gas under these conditions, and was observed with the MIRI IFU aboard JWST as part of the "PDRs4All" program. The resulting IR spectroscopic images of high angular resolution (0.2") reveal a rich observational inventory of mid-IR emission lines, and spatially resolve the substructure of the PDR, with a mosaic cutting perpendicularly across the ionization front and three dissociation fronts. We extracted five spectra that represent the ionized, atomic, and molecular gas layers, and measured the most prominent gas emission lines. An initial analysis summarizes the physical conditions of the gas and the potential of these data. We identified around 100 lines, report an additional 18 lines that remain unidentified, and measured the line intensities and central wavelengths. The H I recombination lines originating from the ionized gas layer bordering the PDR, have intensity ratios that are well matched by emissivity coefficients from H recombination theory, but deviate up to 10% due contamination by He I lines. We report the observed emission lines of various ionization stages of Ne, P, S, Cl, Ar, Fe, and Ni, and show how certain line ratios vary between the five regions. We observe the pure-rotational H$_2$ lines in the vibrational ground state from 0-0 S(1) to 0-0 S(8), and in the first vibrationally excited state from 1-1 S(5) to 1-1 S(9). We derive H$_2$ excitation diagrams, and approximate the excitation with one thermal (~700 K) component representative of an average gas temperature, and one non-thermal component (~2700 K) probing the effect of UV pumping. We compare these results to an existing model for the Orion Bar PDR and highlight the differences with the observations., Comment: 26 pages, 12 figures, 3 tables. Submitted to A&A, under review (1st revision)

Most low-mass stars form in stellar clusters that also contain massive stars, which are sources of far-ultraviolet (FUV) radiation. Theoretical models predict that this FUV radiation produces photo-dissociation regions (PDRs) on the surfaces of protoplanetary disks around low-mass stars, impacting planet formation within the disks. We report JWST and Atacama Large Millimetere Array observations of a FUV-irradiated protoplanetary disk in the Orion Nebula. Emission lines are detected from the PDR; modelling their kinematics and excitation allows us to constrain the physical conditions within the gas. We quantify the mass-loss rate induced by the FUV irradiation, finding it is sufficient to remove gas from the disk in less than a million years. This is rapid enough to affect giant planet formation in the disk.

Aromatic infrared bands (AIBs) are a set of broad emission bands at 3.3, 6.2, 7.7, 8.6, 11.2, and 12.7 $\mu$m, seen in the infrared spectra of most galaxies. With JWST, the 3.3 $\mu$m AIB can in principle be detected up to a redshift of $\sim$ 7. Relating the evolution of the 3.3 $\mu$m AIB to local physical properties of the ISM is thus of paramount importance. By applying a dedicated machine learning algorithm to JWST NIRSpec observations of the Orion Bar photodissociation region obtained as part of the PDRs4All Early Release Science (ERS) program, we extracted two template spectra capturing the evolution of the AIB-related emission in the 3.2-3.6 $\mu$m range, which includes the AIB at 3.3 $\mu$m and its main satellite band at 3.4 $\mu$m. In the Orion Bar, we analyze the spatial distribution of the templates and their relationship with the fluorescent emission of H$_2$ in the near infrared. We find that one template ("AIB$_{\rm Irrad}$") traces regions of neutral atomic gas with strong far-UV fields, while the other template ("AIB$_{\rm Shielded}$") corresponds to shielded regions with lower FUV fields and a higher molecular gas fraction. We then show that these two templates can be used to fit the NIRSpec AIB-related spectra of nearby galaxies. The relative weight of the two templates (AIB$_{\rm Irrad/Shielded}$) is a tracer of the radiative feedback from massive stars on the ISM. We derive an estimate of AIB$_{\rm Irrad/Shielded}$ in a $z$ = 4.22 lensed galaxy, and find that it has a lower value than for local galaxies. This pilot study illustrates how a detailed analysis of AIB emission in nearby regions can be used to probe the physical conditions of the extragalactic ISM., Comment: 10 pages, 7 figures, accepted for publication in A&A

Forty years ago it was proposed that gas phase organic chemistry in the interstellar medium was initiated by the methyl cation CH3+, but hitherto it has not been observed outside the Solar System. Alternative routes involving processes on grain surfaces have been invoked. Here we report JWST observations of CH3+ in a protoplanetary disk in the Orion star forming region. We find that gas-phase organic chemistry is activated by UV irradiation., Comment: Published in Nature

Water is a key ingredient for the emergence of life as we know it. Yet, its destruction and reformation in space remains unprobed in warm gas. Here, we detect the hydroxyl radical (OH) emission from a planet-forming disk exposed to external far-ultraviolet (FUV) radiation with the James Webb Space Telescope. The observations are confronted with the results of quantum dynamical calculations. The highly excited OH infrared rotational lines are the tell-tale signs of H2O destruction by FUV. The OH infrared ro-vibrational lines are attributed to chemical excitation via the key reaction O+H=OH+H which seeds the formation of water in the gas-phase. We infer that the equivalent of the Earth ocean's worth of water is destroyed per month and replenished. These results show that under warm and irradiated conditions water is destroyed and efficiently reformed via gas-phase reactions. This process, assisted by diffusive transport, could reduce the HDO/H2O ratio in the warm regions of planet-forming disks., Comment: Version submitted to Nature Astronomy

To keep current global warming below 1.5{\deg}C compared with the pre-industrial era, measures must be taken as quickly as possible in all spheres of society. Astronomy must also make its contribution. In this proceeding, and during the workshop to which it refers, different levers of actions are discussed through various examples: individual efforts, laboratory-level actions, impact evaluation and mitigation in major projects, institutional level, and involvement through collectives., Comment: 5 pages, 0 figures, SF2A conference 2023

(Abridged) We investigate the impact of radiative feedback from massive stars on their natal cloud and focus on the transition from the HII region to the atomic PDR (crossing the ionisation front (IF)), and the subsequent transition to the molecular PDR (crossing the dissociation front (DF)). We use high-resolution near-IR integral field spectroscopic data from NIRSpec on JWST to observe the Orion Bar PDR as part of the PDRs4All JWST Early Release Science Program. The NIRSpec data reveal a forest of lines including, but not limited to, HeI, HI, and CI recombination lines, ionic lines, OI and NI fluorescence lines, Aromatic Infrared Bands (AIBs including aromatic CH, aliphatic CH, and their CD counterparts), CO2 ice, pure rotational and ro-vibrational lines from H2, and ro-vibrational lines HD, CO, and CH+, most of them detected for the first time towards a PDR. Their spatial distribution resolves the H and He ionisation structure in the Huygens region, gives insight into the geometry of the Bar, and confirms the large-scale stratification of PDRs. We observe numerous smaller scale structures whose typical size decreases with distance from Ori C and IR lines from CI, if solely arising from radiative recombination and cascade, reveal very high gas temperatures consistent with the hot irradiated surface of small-scale dense clumps deep inside the PDR. The H2 lines reveal multiple, prominent filaments which exhibit different characteristics. This leaves the impression of a "terraced" transition from the predominantly atomic surface region to the CO-rich molecular zone deeper in. This study showcases the discovery space created by JWST to further our understanding of the impact radiation from young stars has on their natal molecular cloud and proto-planetary disk, which touches on star- and planet formation as well as galaxy evolution., Comment: 52 pages, 30 figures, submitted to A&A

The JWST has captured the most detailed and sharpest infrared images ever taken of the inner region of the Orion Nebula, the nearest massive star formation region, and a prototypical highly irradiated dense photo-dissociation region (PDR). We investigate the fundamental interaction of far-ultraviolet photons with molecular clouds. The transitions across the ionization front (IF), dissociation front (DF), and the molecular cloud are studied at high-angular resolution. These transitions are relevant to understanding the effects of radiative feedback from massive stars and the dominant physical and chemical processes that lead to the IR emission that JWST will detect in many Galactic and extragalactic environments. Due to the proximity of the Orion Nebula and the unprecedented angular resolution of JWST, these data reveal that the molecular cloud borders are hyper structured at small angular scales of 0.1-1" (0.0002-0.002 pc or 40-400 au at 414 pc). A diverse set of features are observed such as ridges, waves, globules and photoevaporated protoplanetary disks. At the PDR atomic to molecular transition, several bright features are detected that are associated with the highly irradiated surroundings of the dense molecular condensations and embedded young star. Toward the Orion Bar PDR, a highly sculpted interface is detected with sharp edges and density increases near the IF and DF. This was predicted by previous modeling studies, but the fronts were unresolved in most tracers. A complex, structured, and folded DF surface was traced by the H2 lines. This dataset was used to revisit the commonly adopted 2D PDR structure of the Orion Bar. JWST provides us with a complete view of the PDR, all the way from the PDR edge to the substructured dense region, and this allowed us to determine, in detail, where the emission of the atomic and molecular lines, aromatic bands, and dust originate.

(Abridged) Mid-infrared observations of photodissociation regions (PDRs) are dominated by strong emission features called aromatic infrared bands (AIBs). The most prominent AIBs are found at 3.3, 6.2, 7.7, 8.6, and 11.2 $\mu$m. The most sensitive, highest-resolution infrared spectral imaging data ever taken of the prototypical PDR, the Orion Bar, have been captured by JWST. We provide an inventory of the AIBs found in the Orion Bar, along with mid-IR template spectra from five distinct regions in the Bar: the molecular PDR, the atomic PDR, and the HII region. We use JWST NIRSpec IFU and MIRI MRS observations of the Orion Bar from the JWST Early Release Science Program, PDRs4All (ID: 1288). We extract five template spectra to represent the morphology and environment of the Orion Bar PDR. The superb sensitivity and the spectral and spatial resolution of these JWST observations reveal many details of the AIB emission and enable an improved characterization of their detailed profile shapes and sub-components. While the spectra are dominated by the well-known AIBs at 3.3, 6.2, 7.7, 8.6, 11.2, and 12.7 $\mu$m, a wealth of weaker features and sub-components are present. We report trends in the widths and relative strengths of AIBs across the five template spectra. These trends yield valuable insight into the photochemical evolution of PAHs, such as the evolution responsible for the shift of 11.2 $\mu$m AIB emission from class B$_{11.2}$ in the molecular PDR to class A$_{11.2}$ in the PDR surface layers. This photochemical evolution is driven by the increased importance of FUV processing in the PDR surface layers, resulting in a "weeding out" of the weakest links of the PAH family in these layers. For now, these JWST observations are consistent with a model in which the underlying PAH family is composed of a few species: the so-called 'grandPAHs'., Comment: 25 pages, 10 figures, to appear in A&A

We present VLT/MUSE Narrow Field Mode (NFM) observations of a pair of disc-bearing young stellar objects towards the Orion Bar: 203-504 and 203-506. Both of these discs are subject to external photoevaporation, where winds are launched from their outer regions due to environmental irradiation. Intriguingly, despite having projected separation from one another of only 1.65{\arcsec} (660au at 400pc), 203-504 has a classic teardrop shaped ``proplyd'' morphology pointing towards $\theta^2$Ori A (indicating irradiation by the EUV of that star, rather than $\theta^1$ Ori C) but 203-506 has no ionisation front, indicating it is not irradiated by stellar EUV at all. However, 203-506 does show [CI] 8727{\AA} and [OI] 6300{\AA} in emission, indicating irradiation by stellar FUV. This explicitly demonstrates the importance of FUV irradiation in driving mass loss from discs. We conclude that shielding of 203-506 from EUV is most likely due to its position on the observers side of an ionized layer lying in the foreground of the Huygens Region. We demonstrate that the outflow HH 519, previously thought to be emanating from 203-504 is actually an irradiated cloud edge and identify a new compact outflow from that object approximately along our line of sight with a velocity $\sim130$\,km\,s$^{-1}$., Comment: 14 pages, 8 figures. Accepted for publication in MNRAS

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$_2$ 1-0 S(1) at 2.12~$\mu$m, tracing the dissociation front, and the [FeII] and Br$\gamma$ lines, at 1.64 and 2.16~$\mu$m respectively, tracing the ionization front. We obtained narrow-band filter images in these key gas line diagnostic over $\sim 40''$ at spatial scales of $\sim$0.1$''$ ($\sim$0.0002~pc or $\sim$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$_2$ 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$_2$ 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.

Massive stars disrupt their natal molecular cloud material through radiative and mechanical feedback processes. These processes have profound effects on the evolution of interstellar matter in our Galaxy and throughout the Universe, from the era of vigorous star formation at redshifts of 1-3 to the present day. The dominant feedback processes can be probed by observations of the Photo-Dissociation Regions (PDRs) where the far-ultraviolet photons of massive stars create warm regions of gas and dust in the neutral atomic and molecular gas. PDR emission provides a unique tool to study in detail the physical and chemical processes that are relevant for most of the mass in inter- and circumstellar media including diffuse clouds, proto-planetary disks and molecular cloud surfaces, globules, planetary nebulae, and star-forming regions. PDR emission dominates the infrared (IR) spectra of star-forming galaxies. Most of the Galactic and extragalactic observations obtained with the James Webb Space Telescope (JWST) will therefore arise in PDR emission. In this paper we present an Early Release Science program using the MIRI, NIRSpec, and NIRCam instruments dedicated to the observations of an emblematic and nearby PDR: the Orion Bar. These early JWST observations will provide template datasets designed to identify key PDR characteristics in JWST observations. These data will serve to benchmark PDR models and extend them into the JWST era. We also present the Science-Enabling products that we will provide to the community. These template datasets and Science-Enabling products will guide the preparation of future proposals on star-forming regions in our Galaxy and beyond and will facilitate data analysis and interpretation of forthcoming JWST observations., Comment: Submitted to PASP

The James Webb Space Telescope (JWST) was launched on December 25, 2021. This document presents a simulation of the Near Infrared Spectrograph (NIRSpec) observations of the Orion Bar which will be performed as part of the Early Release Sciences (ERS) program "PDRs4all". The methodology to produce this data relies on the use of a direct forward model of the instrument applied to a synthetic scene of the Orion Bar, coupled to format matching in order to deliver data in JWST-pipeline data format. The resulting 3D cube for one order is provided publicly, and is compatible with tools developed by the STScI (e.g. Cubeviz) and with the science enabling products developed by thePDRs4all team. This cube can be used as a template observation for proposers who would like to apply for NIRSpec observations of extended sources with JWST.

Research institutions are bound to contribute to greenhouse gas emission (GHG) reduction efforts for several reasons. First, part of the scientific community's research deals with climate change issues. Second, scientists contribute to students' education: they must be consistent and role models. Third the literature on the carbon footprint of researchers points to the high level of some individual footprints. In a quest for consistency and role models, scientists, teams of scientists or universities have started to quantify their carbon footprints and debate on reduction options. Indeed, measuring the carbon footprint of research activities requires tools designed to tackle its specific features. In this paper, we present an open-source web application, GES 1point5, developed by an interdisciplinary team of scientists from several research labs in France. GES 1point5 is specifically designed to estimate the carbon footprint of research activities in France. It operates at the scale of research labs, i.e. laboratoires, which are the social structures around which research is organized in France and the smallest decision making entities in the French research system. The application allows French research labs to compute their own carbon footprint along a standardized, open protocol. The data collected in a rapidly growing network of labs will be used as part of the Labos 1point5 project to estimate France's research carbon footprint. At the time of submitting this manuscript, 89 research labs had engaged with GES 1point5 to estimate their greenhouse gas emissions. We expect that an international adoption of GES 1point5 (adapted to fit domestic specifics) could contribute to establishing a global understanding of the drivers of the research carbon footprint worldwide and the levers to decrease it.

[Abridged] An increasing number of astronomical instruments (on Earth and space-based) provide hyperspectral images, that is three-dimensional data cubes with two spatial dimensions and one spectral dimension. The intrinsic limitation in spatial resolution of these instruments implies that the spectra associated with pixels of such images are most often mixtures of the spectra of the "pure" components that exist in the considered region. In order to estimate the spectra and spatial abundances of these pure components, we here propose an original blind signal separation (BSS), that is to say an unsupervised unmixing method. Our approach is based on extensions and combinations of linear BSS methods that belong to two major classes of methods, namely nonnegative matrix factorization (NMF) and Sparse Component Analysis (SCA). The former performs the decomposition of hyperspectral images, as a set of pure spectra and abundance maps, by using nonnegativity constraints, but the estimated solution is not unique: It highly depends on the initialization of the algorithm. The considered SCA methods are based on the assumption of the existence of points or tiny spatial zones where only one source is active (i.e., one pure component is present). In real conditions, the assumption of perfect single-source points or zones is not always realistic. In such conditions, SCA yields approximate versions of the unknown sources and mixing coefficients. We propose to use part of these preliminary estimates from the SCA to initialize several runs of the NMF to constrain the convergence of the NMF algorithm. Detailed tests with synthetic data show that the decomposition achieved with such hybrid methods is nearly unique and provides good performance, illustrating the potential of applications to real data., Comment: Accepted for publication in Astronomy and Astrophysics

We investigate inequalities in access to permanent positions in professional astronomy in France, focusing on the hiring stage. We use results from a national survey conducted on behalf of the French society of astronomy and astrophysics (SF2A) aimed at young astronomers holding a PhD obtained in France, and answered by over 300 researchers. We find that women are nearly two times less likely than men to be selected by the (national or local) committees attributing permanent positions ($p=0.06$). We also find that applicants who did their undergraduate studies in an elite school ("Grande \'Ecole"), where women are largely under-represented, rather than in a university, are nearly three times more likely to succeed in obtaining a position ($p=0.0026$). Our analysis suggests the existence of two biases in committees attributing permanent positions in astronomy in France: a gender bias, and a form of elitism. These biases against women in their professional life impacts their personal life as our survey shows that a larger fraction of them declare that having children can have a negative effect on their careers. They are half as many as men having children in the sample. National committees (such as the CNRS) have acknowledged this issue for several years now, hence one can hope that changes will be seen in the next decade.

This paper aims at providing a comprehensive framework to generate an astrophysical scene and to simulate realistic hyperspectral and multispectral data acquired by two JWST instruments, namely NIRCam Imager and NIRSpec IFU. We want to show that this simulation framework can be resorted to assess the benefits of fusing these images to recover an image of high spatial and spectral resolutions. To do so, we create a synthetic scene associated with a canonical infrared source, the Orion Bar. This scene combines pre-existing modelled spectra provided by the JWST Early Release Science Program 1288 and real high resolution spatial maps from the Hubble space and ALMA telescopes. We develop forward models including corresponding noises for the two JWST instruments based on their technical designs and physical features. JWST observations are then simulated by applying the forward models to the aforementioned synthetic scene. We test a dedicated fusion algorithm we developed on these simulated observations. We show the fusion process reconstructs the high spatio-spectral resolution scene with a good accuracy on most areas, and we identify some limitations of the method to be tackled in future works. The synthetic scene and observations presented in the paper are made publicly available and can be used for instance to evaluate instrument models (aboard the JWST or on the ground), pipelines, or more sophisticated algorithms dedicated to JWST data analysis. Besides, fusion methods such as the one presented in this paper are shown to be promising tools to fully exploit the unprecedented capabilities of the JWST.

Hyperspectral imaging has become a significant source of valuable data for astronomers over the past decades. Current instrumental and observing time constraints allow direct acquisition of multispectral images, with high spatial but low spectral resolution, and hyperspectral images, with low spatial but high spectral resolution. To enhance scientific interpretation of the data, we propose a data fusion method which combines the benefits of each image to recover a high spatio-spectral resolution datacube. The proposed inverse problem accounts for the specificities of astronomical instruments, such as spectrally variant blurs. We provide a fast implementation by solving the problem in the frequency domain and in a low-dimensional subspace to efficiently handle the convolution operators as well as the high dimensionality of the data. We conduct experiments on a realistic synthetic dataset of simulated observation of the upcoming James Webb Space Telescope, and we show that our fusion algorithm outperforms state-of-the-art methods commonly used in remote sensing for Earth observation.

The JWST will deliver large data sets of high-quality spectral data over the 0.6-28 $\mu$m range. It will combine sensitivity, spectral and spatial resolution. Specific tools are required to provide efficient scientific analysis of such large data sets. Our aim is to illustrate the potential of unsupervised learning methods to get insights into chemical variations in the populations that carry the aromatic infrared bands (AIBs), more specifically PAH species and carbonaceous very small grains (VSGs). We present a method based on linear fitting and blind signal separation (BSS) for extracting representative spectra for a spectral data set. The method is fast and robust, which ensures its applicability to JWST spectral cubes. We tested this method on a sample of ISO-SWS data, which resemble most closely the JWST spectra in terms of spectral resolution and coverage. Four representative spectra were extracted. Their main characteristics appear consistent with previous studies with populations dominated by cationic PAHs, neutral PAHs, evaporating VSGs, and large ionized PAHs, known as the PAH$^x$ population. In addition, the 3 $\mu$m range, which is considered here for the first time in a BSS method, reveals the presence of aliphatics connected to neutral PAHs. Each representative spectrum is found to carry second-order spectral signatures (e.g. small bands), which are connected with the underlying chemical diversity of populations. However, the precise attribution of theses signatures remains limited by the combined small size and heterogeneity of the sample of astronomical spectra available in this study. The upcoming JWST data will allow us to overcome this limitation. The large data sets of hyperspectral images provided by JWST analysed with the proposed method, which is fast and robust, will open promising perspectives for our understanding of the chemical evolution of the AIB carriers., Comment: 31 pages, 23 figures, accepted for A&A

Dense, fast-moving ejecta knots in supernova remnants are prime sites for molecule and dust formation. We present SOFIA far-IR spectrometer FIFI-LS observations of CO-rich knots in Cas A which cover a ~1 square arc minute area of the northern shell, in the [O III] 52 and 88 micron and [O I] 63 micron lines. The FIFI-LS spectra reveal that the line profiles of [O III] and [O I] are similar to those of the Herschel PACS [O III] and CO lines. We find that the [O III] maps show very different morphology than the [O I] map. The [O III] maps reveal diffuse, large-scale structures and the ratio of the two [O III] lines imply the presence of gas with a range of density 500 - 10,000 per cm^3 within the mapped region. In contrast, the [O I] map shows bright emission associated with the dense CO-rich knots. The 63 micron [O I] line traces cooled, dense post-shocked gas of ejecta. We find that IR-dominated [O III] emission is from post-shocked gas based on its morphology, high column density, and velocity profile. We describe multi-phase ejecta knots, a lifetime of clumps, and survival of dust in the young supernova remnants., Comment: submitted to ApJ

One serious challenge for planet formation is the rapid inward drift of pebble-sized dust particles in protoplanetary disks. Dust trapping at local maxima in the disk gas pressure has received much theoretical attention but still lacks observational support. The cold dust emission in the AB Aur disk forms an asymmetric ring at a radius of about 120 au, which is suggestive of dust trapping in a gas vortex. We present high spatial resolution (0".58x0".78 ~ 80x110 au) NOEMA observations of the 1.12 mm and 2.22 mm dust continuum emission from the AB Aur disk. Significant azimuthal variations of the flux ratio at both wavelengths indicate a size segregation of the large dust particles along the ring. Our continuum images also show that the intensity variations along the ring are smaller at 2.22 mm than at 1.12 mm, contrary to what dust trapping models with a gas vortex have predicted. Our two-fluid (gas+dust) hydrodynamical simulations demonstrate that this feature is well explained if the gas vortex has started to decay due to turbulent diffusion, and dust particles are thus losing the azimuthal trapping on different timescales depending on their size. The comparison between our observations and simulations allows us to constrain the size distribution and the total mass of solid particles in the ring, which we find to be of the order of 30 Earth masses, enough to form future rocky planets., Comment: 9 pags, 4 Figures

Protoplanetary disks undergo substantial mass-loss by photoevaporation, a mechanism which is crucial to their dynamical evolution. However, the processes regulating the gas energetics have not been well constrained by observations so far. We aim at studying the processes involved in disk photoevaporation when it is driven by far-UV photons. We present a unique Herschel survey and new ALMA observations of four externally-illuminated photoevaporating disks (a.k.a. proplyds). For the analysis of these data, we developed a 1D model of the photodissociation region (PDR) of a proplyd, based on the Meudon PDR code and computed the far infrared line emission. We successfully reproduce most of the observations and derive key physical parameters, i.e. densities at the disk surface of about $10^{6}$ cm$^{-3}$ and local gas temperatures of about 1000 K. Our modelling suggests that all studied disks are found in a transitional regime resulting from the interplay between several heating and cooling processes that we identify. These differ from those dominating in classical PDRs, i.e. grain photo-electric effect and cooling by [OI] and [CII] FIR lines. This energetic regime is associated to an equilibrium dynamical point of the photoevaporation flow: the mass-loss rate is self-regulated to set the envelope column density at a value that maintains the temperature at the disk surface around 1000 K. From our best-fit models, we estimate mass-loss rates - of the order of $10^{-7}$ $\mathrm{M}_\odot$/yr - that are in agreement with earlier spectroscopic observation of ionised gas tracers. This holds only if we assume an evaporation flow launched from the disk surface at sound speed (supercritical regime). We have identified the energetic regime regulating FUV-photoevaporation in proplyds. This regime could be implemented into models of the dynamical evolution of protoplanetary disks., Comment: 17 pages (6 more with appendices), 11 figures, accepted for publication in A&A

The near- and mid-IR spectrum of many astronomical objects is dominated by emission bands due to UV-excited polycyclic aromatic hydrocarbons (PAH) and evaporating very small grains (eVSG). Previous studies with the ISO, Spitzer and AKARI space telescopes have shown that the spectral variations of these features are directly related to the local physical conditions that induce a photo-chemical evolution of the band carriers. Because of the limited sensitivity and spatial resolution, these studies have focused mainly on galactic star-forming regions. We discuss how the advent of JWST will allow to extend these studies to previously unresolved sources such as near-by galaxies, and how the analysis of the infrared signatures of PAHs and eVSGs can be used to determine their physical conditions and chemical composition., Comment: To appear in the Proceedings of the annual meeting of the French society of astronomy and astrophysics (SF2A 2015)

Recent studies have confirmed the presence of buckminsterfullerene (C$_{60}$) in different interstellar and circumstellar environments. However, several aspects regarding C$_{60}$ in space are not well understood yet, such as the formation and excitation processes, and the connection between C$_{60}$ and other carbonaceous compounds in the interstellar medium, in particular polycyclic aromatic hydrocarbons (PAHs). In this paper we study several photodissociation regions (PDRs) where C$_{60}$ and PAHs are detected and the local physical conditions are reasonably well constrained, to provide observational insights into these questions. C$_{60}$ is found to emit in PDRs where the dust is cool ($T_d = 20-40$ K) and even in PDRs with cool stars. These results exclude the possibility for C$_{60}$ to be locked in grains at thermal equilibrium in these environments. We observe that PAH and C$_{60}$ emission are spatially uncorrelated and that C$_{60}$ is present in PDRs where the physical conditions (in terms of radiation field and hydrogen density) allow for full dehydrogenation of PAHs, with the exception of Ced 201. We also find trends indicative of an increase in C$_{60}$ abundance within individual PDRs, but these trends are not universal. These results support models where the dehydrogenation of carbonaceous species is the first step towards C$_{60}$ formation. However, this is not the only parameter involved and C$_{60}$ formation is likely affected by shocks and PDR age.

Using the IRAM 30m telescope we have surveyed a $1\times0.8^{\circ}$ part of the Orion molecular cloud in the $^{12}$CO and $^{13}$CO (2-1) lines with a maximal spatial resolution of $\sim$11" and spectral resolution of $\sim$ 0.4 km~s$^{-1}$. The cloud appears filamentary, clumpy and with a complex kinematical structure. We derive an estimated mass of the cloud of 7700 M$_{\text{Sun}}$ (half of which is found in regions with visual extinctions $A_V$ below $\sim$10) and a dynamical age for the nebula of the order of 0.2 Myrs. The energy balance suggests that magnetic fields play an important role in supporting the cloud, at large and small scales. According to our analysis, the turbulent kinetic energy in the molecular gas due to outflows is comparable to turbulent kinetic energy resulting from the interaction of the cloud with the HII region. This latter feedback appears negative, i.e. the triggering of star formation by the HII region is inefficient in Orion. The reduced data as well as additional products such as the column density map are made available online at http://userpages.irap.omp.eu/~oberne/Olivier_Berne/Data.

The mid-infrared (IR; 5-15~$\mu$m) spectrum of a wide variety of astronomical objects exhibits a set of broad emission features at 6.2, 7.7, 8.6, 11.3 and 12.7 $\mu$m. About 30 years ago it was proposed that these signatures are due to emission from a family of UV heated nanometer-sized carbonaceous molecules known as polycyclic aromatic hydrocarbons (PAHs), causing them to be referred to as aromatic IR bands (AIBs). Today, the acceptance of the PAH model is far from settled, as the identification of a single PAH in space has not yet been successful and physically relevant theoretical models involving ``true'' PAH cross sections do not reproduce the AIBs in detail. In this paper, we use the NASA Ames PAH IR Spectroscopic Database, which contains over 500 quantum-computed spectra, in conjunction with a simple emission model, to show that the spectrum produced by any random mixture of at least 30 PAHs converges to the same 'kernel'-spectrum. This kernel-spectrum captures the essence of the PAH emission spectrum and is highly correlated with observations of AIBs, strongly supporting PAHs as their source. Also, the fact that a large number of molecules are required implies that spectroscopic signatures of the individual PAHs contributing to the AIBs spanning the visible, near-infrared, and far infrared spectral regions are weak, explaining why they have not yet been detected. An improved effort, joining laboratory, theoretical, and observational studies of the PAH emission process, will support the use of PAH features as a probe of physical and chemical conditions in the nearby and distant Universe.

We report a Herschel detection of high-J rotational CO lines from a dense knot in the supernova remnant Cas A. Based on a combined analysis of these rotational lines, and previously observed ro-vibrational CO lines, we find the gas to be warm (two components at 400 and 2000 K) and dense (1e6-7 cm-3), with a CO column density of 5e17 cm-2. This, along with the broad line widths (400 kms-1), suggests that the CO emission originates in the post-shock region of the reverse shock. As the passage of the reverse shock dissociates any existing molecules, the CO has most likely reformed in the last few years, in the post-shock gas. The CO cooling time is comparable to the CO formation time, so possible heating sources (UV photons from the shock front, X-rays, electron conduction) to maintain the large column density of warm CO are discussed., Comment: Accepted to A&A Letters, 4 pages, 4 figures

In the early 90s, C60+ was proposed as the carrier of two diffuse interstellar bands (DIBs) at 957.7 and 963.2 nm, but a firm identification still awaits gas-phase spectroscopic data. Neutral C60, on the other hand, was recently detected through its infrared emission bands in the interstellar medium and evolved stars. In this contribution, we present the detection of C60+ through its infrared vibrational bands in the NGC 7023 nebula, based on spectroscopic observations with the Spitzer space telescope, quantum chemistry calculation, and laboratory data from the literature. This detection supports the idea that C60+ could be a DIB carrier, and provides robust evidence that fullerenes exist in the gas-phase in the interstellar medium. Modeling efforts to design specific observations, combined with new gas-phase data, will be essential to confirm this proposal. A definitive attribution of the 957.7 and 963.2 nm DIBs to C60+ would represent a significant step forward in the field., Comment: To appear in "Proceedings of IAU 297 symposium on the Diffuse Interstellar Bands", eds. J. Cami and N. Cox (5 pages)

We observationally investigate the relation between the photoelectric heating efficiency in PDRs and the charge of PAHs, which are considered to play a key role in photoelectric heating. Using PACS onboard Herschel, we observed six PDRs spanning a wide range of FUV radiation fields (G_0=100-10^5). To measure the photoelectric heating efficiency, we obtained the intensities of the main cooling lines, i.e., the [OI]63um, 145um, and [CII]158um, as well as the FIR continuum intensity. We used Spitzer/IRS spectroscopic mapping observations to investigate the MIR PAH features in the same regions. We decomposed the MIR PAH emission into that of neutral (PAH^0) and positively ionized (PAH^+) species to derive the fraction of the positively charged PAHs, and compare it to the photoelectric heating efficiency. The heating efficiency traced by ([OI]63um+[OI]145um+[CII]158um) / TIR, ranges between 0.1% and 0.9% in different sources, and the fraction of PAH^+ relative to (PAH^0 + PAH^+) spans from 0(+11)% to 87(+/-10)%. All positions with a high PAH^+ fraction show a low heating efficiency, and all positions with a high heating efficiency have a low PAH^+ fraction, supporting the scenario in which a positive grain charge results in a decreased heating efficiency. Theoretical estimates of the photoelectric heating efficiency show a stronger dependence on the charging parameter gamma=G_0 T^{1/2}/n_e than the observed efficiency reported in this study, and the discrepancy is significant at low gamma. The photoelectric heating efficiency on PAHs, traced by ([OI]63um+[OI]145um+[CII]158um) / (PAH+[OI]63um+[OI]145um+[CII]158um), shows a much better match between the observations and the theoretical estimates. The good agreement of the photoelectric heating efficiency on PAHs with a theoretical model indicates the dominant contribution of PAHs to the photoelectric heating. (abridged for arXiv), Comment: 14 pages with 10 figures, accepted for publication in A&A

Proplyds are photodissociation region (PDR)-like cometary cocoons around young stars which are thought to originate through photo-evaporation of the central protoplanetary disk by external UV radiation from the nearby OB stars. This letter presents spatially resolved mid-infrared imaging and spectroscopy of the proplyd HST10 obtained with the VLT/VISIR instrument. These observations allow us to detect Polycyclic Aromatic Hydrocarbons (PAH) emission in the proplyd photodissociation region and to study the general properties of PAHs in proplyds for the first time. We find that PAHs in HST10 are mostly neutral and at least 50 times less abundant than typical values found for the diffuse ISM or the nearby Orion Bar. With such a low PAH abundance, photoelectric heating is significantly reduced. If this low abundance pertains also to the original disk material, gas heating rates could be too low to efficiently drive photoevaporation unless other processes can be identified. Alternatively, the model behind the formation of proplyds as evaporating disks may have to be revised., Comment: 5 pages, 3 figures, 1 table

Buckminsterfullerene (C60) was recently detected through its infrared emission bands in the interstellar medium (ISM), including in the proximity of massive stars, where physical conditions could favor the formation of the cationic form, C60+. In addition, C60+ was proposed as the carrier of two diffuse interstellar bands in the near-IR, although a firm identification still awaits for gas-phase spectroscopic data. We examined in details the Spitzer IRS spectra of the NGC 7023 reflection nebula, at a position close (7.5") to the illuminating B star HD 200775, and found four previously unreported bands at 6.4, 7.1, 8.2 and 10.5 \mu m in addition to the classical bands attributed to Polycylic Aromatic Hydrocarbons (PAHs) and neutral C60. These 4 bands are observed only in this region of the nebula, while C60 emission is still present slightly further away from the star, and PAH emission even further away. Based on this observation, on theoretical calculations we perform, and on laboratory studies, we attribute these bands to C60+. The detection of C60+ confirms the idea that large carbon molecules exist in the gas-phase in these environments. In addition, the relative variation of the C60, and C60+, band intensities constitutes a potentially powerful probe of the physical conditions in highly UV-irradiated regions., Comment: Accepted for publication in A&A, v2: minor corrections of typos and language and additional reference included

Hydrodynamical instabilities are believed to power some of the small scale (0.1-10 pc) turbulence and chemical mixing in the interstellar medium. Identifying such instabilities has always been difficult but recent observations of a wavelike structure (the Ripples) in the Orion nebula have been interpreted as a signature of the Kelvin-Helmholtz instability (KHI), occurring at the interface between the HII region and the molecular cloud. However, this has not been verified theoretically. In this letter, we investigate theoretically the stability of this interface using observational constraints for the local physical conditions. A linear analysis shows that the HII/molecular cloud interface is indeed KH unstable for a certain range of magnetic field orientation. We find that the maximal growth-rates correspond to typical timescales of a few 1e4 years and instability wavelengths of 0.06 to 0.6 pc. We predict that after 2e5 years the KHI saturates and forms a turbulent layer of thickness ~0.5 pc. The KHI can remain in linear phase over a maximum distance of 0.75 pc. These spatial and time scales are compatible with the Ripples representing the linear phase of the KHI. These results suggest that the KHI may be crucial to generate turbulence and to bring heavy elements injected by the winds of massive stars in HII regions to colder regions where planetary systems around low mass stars are being formed. This could apply to the transport of 26Al injected by a massive star in an HII region to the nascent solar-system., Comment: Accepted for publication in ApJ Letters

Large spatial-spectral surveys are more and more common in astronomy. This calls for the need of new methods to analyze such mega- to giga-pixel data-cubes. In this paper we present a method to decompose such observations into a limited and comprehensive set of components. The original data can then be interpreted in terms of linear combinations of these components. The method uses non-negative matrix factorization (NMF) to extract latent spectral end-members in the data. The number of needed end-members is estimated based on the level of noise in the data. A Monte-Carlo scheme is adopted to estimate the optimal end-members, and their standard deviations. Finally, the maps of linear coefficients are reconstructed using non-negative least squares. We apply this method to a set of hyperspectral data of the NGC 7023 nebula, obtained recently with the HIFI instrument onboard the Herschel space observatory, and provide a first interpretation of the results in terms of 3-dimensional dynamical structure of the region., Comment: Proceedings of the 2012 meeting of the french astronomical society (SF2A) in Nice

Buckminsterfullerene (C60) was recently confirmed to be the largest molecule identified in space. However, it remains unclear how, and where this molecule is formed. It is generally believed that C60 is formed from the build up of small carbonaceous compounds, in the hot and dense envelopes of evolved stars. Analyzing infrared observations, obtained by Spitzer and Herschel, we found that C60 is efficiently formed in the tenuous and cold environment of an interstellar cloud illuminated by strong ultraviolet (UV) radiation fields. This implies that another formation pathway, efficient at low densities, must exist. Based on recent laboratory and theoretical studies, we argue that Polycyclic Aromatic Hydrocarbons are converted into graphene, and subsequently C60, under UV irradiation from massive stars. This shows that alternative - top-down - routes are key to understanding the organic inventory in space., Comment: 22 pages, 4 figures, 3 sup. figures, 1 sup. video

Herbig Ae stars (HAe) are the precursors of Vega-type systems and, therefore, crucial objects in planet formation studies. Thus far, only a few disks associated with HAe stars have been studied using millimetre interferometers. Our aim is to determine the dust evolution and the lifetime of the disks associated with Herbig Ae stars. We imaged the continuum emission at 3 mm and 1.3 mm of the Herbig Ae/Be stars BD+61154, RR Tau, VY Mon and LkHa 198 using the Plateau de Bure Interferometer (PdBI). These stars are in the upper end of the stellar mass range of the Herbig Ae stars (stellar mass greater than 3 solar masses). Our measurements were used to complete the Spectral Energy Distribution (SED). The modelling of the SED, in particular the FIR-mm part, allow us to determine the masses and dust properties of these disks. We detected the disks associated with BD+61154, RR Tau and VY Mon with disk masses of 0.35 Msun, 0.05 Msun and 0.40 Msun respectively. The disk around LkHa 198 was not detected with an upper limit to the disk mass of 0.004 Msun. We detected, however, the disks associated with the younger stellar objects LkHa 198--IR and LkHa 198-mm that are located in the vicinity of LkHa 198. The fitting of the mm part of the SED reveal that the grains in the mid-plane of the disks around BD+61154, RR Tau and VY Mon have sizes of 1--1000 microns. Therefore, grains have not grown to centimetre sizes in these disks yet. These massive (M>3 Msun) and young (about 1 Myr) HAe stars are surrounded by massive (>0.04 Msun) disks with grains of micron-millimetre sizes. Although grain growth is proceeding in these disks, their evolutionary stage is prior to the formation of planetesimals. These disks are less evolved than those detected around T Tauri and Herbig Be stars.

Massive stars influence their parental molecular cloud, and it has long been suspected that the development of hydrodynamical instabilities can compress or fragment the cloud. Identifying such instabilities has proved difficult. It has been suggested that elongated structures (such as the `pillars of creation') and other shapes arise because of instabilities, but alternative explanations are available. One key signature of an instability is a wave-like structure in the gas, which has hitherto not been seen. Here we report the presence of `waves' at the surface of the Orion molecular cloud near where massive stars are forming. The waves seem to be a Kelvin-Helmholtz instability that arises during the expansion of the nebula as gas heated and ionized by massive stars is blown over pre-existing molecular gas., Comment: Preprint of publication in Nature

We present the detection and characterization of a peculiar low-mass protostar (IRAS 22129+7000) located ~0.4 pc from Ced 201 Photodissociation Region (PDR) and ~0.2 pc from the HH450 jet. The cold circumstellar envelope surrounding the object has been mapped through its 1.2 mm dust continuum emission with IRAM-30m/MAMBO. The deeply embedded protostar is clearly detected with Spitzer/MIPS (70 um), IRS (20-35 um) and IRAC (4.5, 5.8, and 8 um) but also in the K_s band (2.15 um). Given the large "near- and mid-IR excess" in its spectral energy distribution, but large submillimeter-to-bolometric luminosity ratio (~2%), IRAS 22129+7000 must be a transition Class 0/I source and/or a multiple stellar system. Targeted observations of several molecular lines from CO, 13CO, C18O, HCO+ and DCO+ have been obtained. The presence of a collimated molecular outflow mapped with the CSO telescope in the CO J=3-2 line suggests that the protostar/disk system is still accreting material from its natal envelope. Indeed, optically thick line profiles from high density tracers such as HCO+ J=1-0 show a red-shifted-absorption asymmetry reminiscent of inward motions. We construct a preliminary physical model of the circumstellar envelope (including radial density and temperature gradients, velocity field and turbulence) that reproduces the observed line profiles and estimates the ionization fraction. The presence of both mechanical and (non-ionizing) FUV-radiative input makes the region an interesting case to study triggered star formation.

The use of Blind Signal Separation methods (ICA and other approaches) for the analysis of astrophysical data remains quite unexplored. In this paper, we present a new approach for analyzing the infrared emission spectra of interstellar dust, obtained with NASA's Spitzer Space Telescope, using FastICA and Non-negative Matrix Factorization (NMF). Using these two methods, we were able to unveil the source spectra of three different types of carbonaceous nanoparticles present in interstellar space. These spectra can then constitute a basis for the interpretation of the mid-infrared emission spectra of interstellar dust in the Milky Way and nearby galaxies. We also show how to use these extracted spectra to derive the spatial distribution of these nanoparticles.