Showing 1–50 of 94 results for author: Ysard, N

JWST observations of the Horsehead photon-dominated region I. First results from multi-band near- and mid-infrared imaging

Authors: A. Abergel, K. Misselt, K. D. Gordon, A. Noriega-Crespo, P. Guillard, D. Van De Putte, A. N. Witt, N. Ysard, M. Baes, H. Beuther, P. Bouchet, B. R. Brandl, M. Elyajouri, O. Kannavou, S. Kendrew, P. Klassen, B. Trahin

Abstract: The JWST has captured the sharpest IR images ever taken of the Horsehead nebula, a prototypical moderately irradiated PDR that is fully representative of most of the UV-illuminated molecular gas in the Milky Way and star-forming galaxies. We investigate the impact of FUV radiation of a molecular cloud and constrain the structure of the edge of the PDR and its illumination conditions. We used NIRCa… ▽ More The JWST has captured the sharpest IR images ever taken of the Horsehead nebula, a prototypical moderately irradiated PDR that is fully representative of most of the UV-illuminated molecular gas in the Milky Way and star-forming galaxies. We investigate the impact of FUV radiation of a molecular cloud and constrain the structure of the edge of the PDR and its illumination conditions. We used NIRCam and MIRI to obtain 17 broadband and 6 narrowband maps from 0.7 to 28 $μみゅー$m. We mapped the dust emission, scattered light, and several gas phase lines. We also used HST-WFC3 maps at 1.1 and 1. 6 $μみゅー$m, along with HST-STIS spectroscopic observations of the H$αあるふぁ$ line. We probed the structure of the edge of the Horsehead and resolved its spatial complexity. We detected a network of faint striated features extending perpendicularly to the PDR front into the H\,II region in filters sensitive to nano-grain emission and light scattered by larger grains. This may indeed figure as the first detection of the entrainment of dust particles in the evaporative flow. The map of the 1-0 S(1) line of H$_2$ presents sharp sub-structures on scales as small as 1.5 arcsec. The ionization and dissociation fronts appear at distances 1-2 arcsec behind the edge of the PDR and seem to spatially coincide, indicating a thickness of the neutral atomic layer below 100 au. All broadband maps present strong color variations which can be explained by dust attenuation. Deviations of the emissions in the H$αあるふぁ$, Pa$αあるふぁ,$ and Br$αあるふぁ$ lines also indicate dust attenuation. With a very simple model, we derive the main features of the extinction curve. A small excess of extinction at 3 $μみゅー$m may be attributed to icy H$_2$O mantles onto grains. In all lines of sight crossing the inner regions of the Horsehead, it appears that dust attenuation is non-negligible over the entire spectral range. △ Less

Submitted 24 April, 2024; originally announced April 2024.

Comments: 25 pages, 17 figures

PDRs4All VIII: Mid-IR emission line inventory of the Orion Bar

Authors: Dries Van De Putte, Raphael Meshaka, Boris Trahin, Emilie Habart, Els Peeters, Olivier Berné, Felipe Alarcón, Amélie Canin, Ryan Chown, Ilane Schroetter, Ameek Sidhu, Christiaan Boersma, Emeric Bron, Emmanuel Dartois, Javier R. Goicoechea, Karl D. Gordon, Takashi Onaka, Alexander G. G. M. Tielens, Laurent Verstraete, Mark G. Wolfire, Alain Abergel, Edwin A. Bergin, Jeronimo Bernard-Salas, Jan Cami, Sara Cuadrado , et al. (113 additional authors not shown)

Abstract: 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 observat… ▽ More 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. △ Less

Submitted 3 April, 2024; originally announced April 2024.

Comments: 26 pages, 12 figures, 3 tables. Submitted to A&A, under review (1st revision)

Journal ref: A&A 687, A86 (2024)

A far-ultraviolet-driven photoevaporation flow observed in a protoplanetary disk

Authors: Olivier Berné, Emilie Habart, Els Peeters, Ilane Schroetter, Amélie Canin, Ameek Sidhu, Ryan Chown, Emeric Bron, Thomas J. Haworth, Pamela Klaassen, Boris Trahin, Dries Van De Putte, Felipe Alarcón, Marion Zannese, Alain Abergel, Edwin A. Bergin, Jeronimo Bernard-Salas, Christiaan Boersma, Jan Cami, Sara Cuadrado, Emmanuel Dartois, Daniel Dicken, Meriem Elyajouri, Asunción Fuente, Javier R. Goicoechea , et al. (121 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 29 February, 2024; originally announced March 2024.

Journal ref: Science, 383, 6686, 2024

The magnetic field in the Flame nebula

Authors: Ivana Bešlić, Simon Coudé, Dariusz C. Lis, Maryvonne Gerin, Paul F. Goldsmith, Jerome Pety, Antoine Roueff, Karine Demyk, Charles D. Dowell, Lucas Einig, Javier R. Goicoechea, Francois Levrier, Jan Orkisz, Nicolas Peretto, Miriam G. Santa-Maria, Nathalie Ysard, Antoine Zakardjian

Abstract: Star formation is essential in galaxy evolution and the cycling of matter. The support of interstellar clouds against gravitational collapse by magnetic (B-) fields has been proposed to explain the low observed star formation efficiency in galaxies and the Milky Way. Despite the Planck satellite providing a 5-15' all-sky map of the B-field geometry in the diffuse interstellar medium, higher spatia… ▽ More Star formation is essential in galaxy evolution and the cycling of matter. The support of interstellar clouds against gravitational collapse by magnetic (B-) fields has been proposed to explain the low observed star formation efficiency in galaxies and the Milky Way. Despite the Planck satellite providing a 5-15' all-sky map of the B-field geometry in the diffuse interstellar medium, higher spatial resolution observations are required to understand the transition from diffuse gas to gravitationally unstable filaments. NGC 2024, the Flame Nebula, in the nearby Orion B molecular cloud, contains a young, expanding HII region and a dense filament that harbors embedded protostellar objects. Therefore, NGC 2024 is an excellent opportunity to study the role of B-fields in the formation, evolution, and collapse of filaments, as well as the dynamics and effects of young HII regions on the surrounding molecular gas. We combine new 154 and 216 micron dust polarization measurements carried out using the HAWC+ instrument aboard SOFIA with molecular line observations of 12CN(1-0) and HCO+(1-0) from the IRAM 30-meter telescope to determine the B-field geometry and to estimate the plane of the sky magnetic field strength across the NGC 2024. The HAWC+ observations show an ordered B-field geometry in NGC 2024 that follows the morphology of the expanding HII region and the direction of the main filament. The derived plane of the sky B-field strength is moderate, ranging from 30 to 80 micro G. The strongest B-field is found at the northern-west edge of the HII region, characterized by the highest gas densities and molecular line widths. In contrast, the weakest field is found toward the filament in NGC 2024. The B-field has a non-negligible influence on the gas stability at the edges of the expanding HII shell (gas impacted by the stellar feedback) and the filament (site of the current star formation). △ Less

Submitted 7 February, 2024; v1 submitted 30 January, 2024; originally announced January 2024.

Comments: 36 pages, 26 figures Accepted for publication in Astronomy & Astrophysics

THEMIS 2.0: A self-consistent model for dust extinction, emission, and polarisation

Authors: Nathalie Ysard, Anthony Peter Jones, Vincent Guillet, Karine Demyk, Marjorie Decleir, Laurent Verstraete, Ilyes Choubani, Marc-Antoine Miville-Deschênes, Lapo Fanciullo

Abstract: Recent observations in emission, extinction, and polarisation have at least partially invalidated most of the astronomical standard grain models for the diffuse ISM. Moreover, lab measurements on interstellar silicate analogues have shown differences with the optical properties used in these standard models. To address these issues, our objective is twofold: (i) to update the optical properties of… ▽ More Recent observations in emission, extinction, and polarisation have at least partially invalidated most of the astronomical standard grain models for the diffuse ISM. Moreover, lab measurements on interstellar silicate analogues have shown differences with the optical properties used in these standard models. To address these issues, our objective is twofold: (i) to update the optical properties of silicates and (ii) to develop the THEMIS dust model to allow the calculation of polarised extinction and emission. Based on optical constants measured in the lab for amorphous silicates and on observational constraints in mid-IR extinction and X-ray scattering, we defined new optical constants for the THEMIS silicates. Absorption and scattering efficiencies for spheroidal grains were then derived with the discrete dipole approximation. These new optical properties make it possible to explain the dust emission and extinction, both total and polarised. The model is not yet pushed to its limits since it does not require the perfect alignment of all grains to explain the observations and it therefore has the potential to accommodate the highest polarisation levels inferred from extinction measures. Moreover, the dispersion of the optical properties of the different lab silicates naturally explain the variations in both the total and polarised emission and extinction observed in the diffuse ISM. A single, invariant model calibrated on one single set of observations is obsolete for explaining contemporary observations. We are proposing a completely flexible dust model based entirely on lab measurements that has the potential to make major advances in understanding the nature of ISM grains and how they evolve as a function of their environment. Even if challenging, this is also relevant for future missions that will aim to perform precise measurements of the CMB spectral distortions and polarisation. △ Less

Submitted 30 May, 2024; v1 submitted 15 January, 2024; originally announced January 2024.

Comments: Accepted for publication by A&A ; abstract modified to match the arxiv standard ; typo corrected on 19/01/2024 in Tab. 3 ; units corrected on 07/02/2024 in Figs. 9 and C.2

PDRs4All. V. Modelling the dust evolution across the illuminated edge of the Orion Bar

Authors: M. Elyajouri, N. Ysard, A. Abergel, E. Habart, L. Verstraete, A. Jones, M. Juvela, T. Schirmer, R. Meshaka, E. Dartois, J. Lebourlot, G. Rouille, T. Onaka, E. Peeters, O. Berne, F. Alarcon, J. Bernard-Salas, M. Buragohain, J. Cami, A. Canin, R. Chown, K. Demyk, K. Gordon, O. Kannavou, M. Kirsanova , et al. (9 additional authors not shown)

Abstract: We study the emission of dust grains within the Orion Bar - a well-known, highly far-UV (FUV)-irradiated PDR. The Orion Bar because of its edge-on geometry provides an exceptional benchmark for characterizing dust evolution and the associated driving processes under varying physical conditions. Our goal is to constrain the local properties of dust by comparing its emission to models. Taking advant… ▽ More We study the emission of dust grains within the Orion Bar - a well-known, highly far-UV (FUV)-irradiated PDR. The Orion Bar because of its edge-on geometry provides an exceptional benchmark for characterizing dust evolution and the associated driving processes under varying physical conditions. Our goal is to constrain the local properties of dust by comparing its emission to models. Taking advantage of the recent JWST PDRs4All data, we follow the dust emission as traced by JWST NIRCam (at 3.35 and 4.8 micron) and MIRI (at 7.7, 11.3, 15.0, and 25.5 micron), along with NIRSpec and MRS spectroscopic observations. First, we constrain the minimum size and hydrogen content of carbon nano-grains from a comparison between the observed dust emission spectra and the predictions of the THEMIS dust model coupled to the numerical code DustEM. Using this dust model, we then perform 3D radiative transfer simulations of dust emission with the SOC code and compare to data obtained along well chosen profiles across the Orion Bar. The JWST data allows us, for the first time, to spatially resolve the steep variation of dust emission at the illuminated edge of the Orion Bar PDR. By considering a dust model with carbonaceous nano-grains and submicronic coated silicate grains, we derive unprecedented constraints on the properties of across the Orion Bar. To explain the observed emission profiles with our simulations, we find that the nano-grains must be strongly depleted with an abundance (relative to the gas) 15 times less than in the diffuse ISM. The NIRSpec and MRS spectroscopic observations reveal variations in the hydrogenation of the carbon nano-grains. The lowest hydrogenation levels are found in the vicinity of the illuminating stars suggesting photo-processing while more hydrogenated nano-grains are found in the cold and dense molecular region, potentially indicative of larger grains. △ Less

Submitted 2 January, 2024; originally announced January 2024.

Comments: 18 pages, 11 figures, accepted for publication in A&A

Journal ref: A&A 685, A76 (2024)

PDRs4All III: JWST's NIR spectroscopic view of the Orion Bar

Authors: Els Peeters, Emilie Habart, Olivier Berne, Ameek Sidhu, Ryan Chown, Dries Van De Putte, Boris Trahin, Ilane Schroetter, Amelie Canin, Felipe Alarcon, Bethany Schefter, Baria Khan, Sofia Pasquini, Alexander G. G. M. Tielens, Mark G. Wolfire, Emmanuel Dartois, Javier R. Goicoechea, Alexandros Maragkoudakis, Takashi Onaka, Marc W. Pound, Silvia Vicente, Alain Abergel, Edwin A. Bergin, Jeronimo Bernard-Salas, Christiaan Boersma , et al. (113 additional authors not shown)

Abstract: (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… ▽ More (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. △ Less

Submitted 12 October, 2023; originally announced October 2023.

Comments: 52 pages, 30 figures, submitted to A&A

Journal ref: A&A 685, A74 (2024)

Exploring the interstellar medium of NGC 891 at millimeter wavelengths using the NIKA2 camera

Authors: S. Katsioli, R. Adam, P. Ade, H. Ajeddig, P. André, E. Artis, H. Aussel, M. Baes, A. Beelen, A. Benoît, S. Berta, L. Bing, O. Bourrion, M. Calvo, A. Catalano, C. J. R. Clark, I. De Looze, M. De Petris, F. -X. Désert, S. Doyle, E. F. C. Driessen, G. Ejlali, M. Galametz, F. Galliano, A. Gomez , et al. (39 additional authors not shown)

Abstract: In the framework of the IMEGIN Large Program, we used the NIKA2 camera on the IRAM 30-m telescope to observe the edge-on galaxy NGC 891 at 1.15 mm and 2 mm and at a FWHM of 11.1" and 17.6", respectively. Multiwavelength data enriched with the new NIKA2 observations fitted by the HerBIE SED code (coupled with the THEMIS dust model) were used to constrain the physical properties of the ISM. Emission… ▽ More In the framework of the IMEGIN Large Program, we used the NIKA2 camera on the IRAM 30-m telescope to observe the edge-on galaxy NGC 891 at 1.15 mm and 2 mm and at a FWHM of 11.1" and 17.6", respectively. Multiwavelength data enriched with the new NIKA2 observations fitted by the HerBIE SED code (coupled with the THEMIS dust model) were used to constrain the physical properties of the ISM. Emission originating from the diffuse dust disk is detected at all wavelengths from mid-IR to mm, while mid-IR observations reveal warm dust emission from compact HII regions. Indications of mm excess emission have also been found in the outer parts of the galactic disk. Furthermore, our SED fitting analysis constrained the mass fraction of the small (< 15 Angstrom) dust grains. We found that small grains constitute 9.5% of the total dust mass in the galactic plane, but this fraction increases up to ~ 20% at large distances (|z| > 3 kpc) from the galactic plane. △ Less

Submitted 6 October, 2023; originally announced October 2023.

Comments: To appear in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Sciences

Journal ref: EPJ Web of conferences 293 (2024) 00026

The stratification of ISM properties in the edge-on galaxy NGC 891 revealed by NIKA2

Authors: S. Katsioli, E. M. Xilouris, C. Kramer, R. Adam, P. Ade, H. Ajeddig, P. André, E. Artis, H. Aussel, M. Baes, A. Beelen, A. Benoît, S. Berta, L. Bing, O. Bourrion, M. Calvo, A. Catalano, C. J. R. Clark, I. De Looze, M. De Petris, F. -X. Désert, S. Doyle, E. F. C. Driessen, G. Ejlali, M. Galametz , et al. (38 additional authors not shown)

Abstract: As the millimeter wavelength range remains a largely unexplored spectral region for galaxies, the IMEGIN large program aims to map the millimeter continuum emission of 22 nearby galaxies at 1.15 and 2 mm. Using the high-resolution maps produced by the NIKA2 camera, we explore the existence of very cold dust and take possible contamination by free-free and synchrotron emission into account. We stud… ▽ More As the millimeter wavelength range remains a largely unexplored spectral region for galaxies, the IMEGIN large program aims to map the millimeter continuum emission of 22 nearby galaxies at 1.15 and 2 mm. Using the high-resolution maps produced by the NIKA2 camera, we explore the existence of very cold dust and take possible contamination by free-free and synchrotron emission into account. We study the IR-to-radio emission coming from different regions along the galactic plane and at large vertical distances. New observations of NGC 891, using the NIKA2 camera on the IRAM 30m telescope, along with a suite of observations at other wavelengths were used to perform a multiwavelength study of the spectral energy distribution in the interstellar medium in this galaxy. This analysis was performed globally and locally, using the advanced hierarchical Bayesian fitting code, HerBIE, coupled with the THEMIS dust model. Our dust modeling is able to reproduce the near-IR to millimeter emission of NGC 891, with the exception of an excess at a level of 25% obtained by the NIKA2 observations in the outermost parts of the disk. The radio continuum and thermal dust emission are distributed differently in the disk and galaxy halo. Different dusty environments are also revealed by a multiwavelength investigation of the emission features. Our detailed decomposition at millimeter and centimeter wavelengths shows that emission at 1 mm is purely originated by dust. Radio components become progressively important with increasing wavelengths. Finally, we find that emission arising from small dust grains accounts for ~ 9.5% of the total dust mass, reaching up to 20% at large galactic latitudes. Shock waves in the outflows that shatter the dust grains might explain this higher fraction of small grains in the halo. △ Less

Submitted 15 September, 2023; originally announced September 2023.

Comments: 15 pages, 10 figures, Accepted for publication in A&A

Journal ref: A&A 679, A7 (2023)

PDRs4All II: JWST's NIR and MIR imaging view of the Orion Nebula

Authors: Emilie Habart, Els Peeters, Olivier Berné, Boris Trahin, Amélie Canin, Ryan Chown, Ameek Sidhu, Dries Van De Putte, Felipe Alarcón, Ilane Schroetter, Emmanuel Dartois, Sílvia Vicente, Alain Abergel, Edwin A. Bergin, Jeronimo Bernard-Salas, Christiaan Boersma, Emeric Bron, Jan Cami, Sara Cuadrado, Daniel Dicken, Meriem Elyajouri, Asunción Fuente, Javier R. Goicoechea, Karl D. Gordon, Lina Issa , et al. (117 additional authors not shown)

Abstract: 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 fron… ▽ More 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. △ Less

Submitted 2 September, 2023; v1 submitted 31 August, 2023; originally announced August 2023.

Journal ref: A&A 685, A73 (2024)

Low-temperature optical constants of amorphous silicate dust analogues

Authors: K. Demyk, V. Gromov, C. Meny, N. Ysard, D. Paradis, A. P. Jones, D. Petitprez, P. Hubert, H. Leroux, C. Nayral, F. Delpech

Abstract: Cosmic dust models are key ingredients in advancing our understanding of astronomical environments as diverse as interstellar clouds in galaxies, circumstellar envelopes around evolved and young stars, and protoplanetary disks. Such models consist of several dust populations, each with different compositions and size distributions. They may also consider different grain shapes, although most model… ▽ More Cosmic dust models are key ingredients in advancing our understanding of astronomical environments as diverse as interstellar clouds in galaxies, circumstellar envelopes around evolved and young stars, and protoplanetary disks. Such models consist of several dust populations, each with different compositions and size distributions. They may also consider different grain shapes, although most models assume spherical grains. All include a component of silicate dust. The absorption and emission properties of these dust components are calculated from the optical constants of each dust material which have various experimental, phenomenological, and theoretical origins depending on the model. We aim to provide the community with new sets of optical constants for amorphous silicate dust analogues at low temperatures. The analogues consist of four Mg-rich silicate samples of stoichiometry ranging from enstatite to olivine, and of eight samples of Mg and Fe rich silicates with a pyroxene stoichiometry and differing magnesium and iron content. We calculated the optical constants from transmission measurements using the Kramers-Kronig relations, assuming that the grains are small compared to the wavelength and prolate in shape with axis ratios of 1.5 and 2 for the Mg and Fe rich samples, respectively. New optical constants for silicate dust analogues were calculated over the wavelength range from 5 to 800-1000 microns, depending on the sample, and at temperatures of 10, 30, 100, 200, and 300 K.We determined the uncertainties on the derived optical constants based on the assumptions used to calculate them. To facilitate the use of these data in cosmic dust models, we provide optical constants extrapolated outside the measured spectral range into the UV-NIR and mm-cm wavelength ranges, as well as formulae that can be used to interpolate them at any temperature in the range 10 - 300 K. △ Less

Submitted 14 September, 2022; originally announced September 2022.

Journal ref: A&A 666, A192 (2022)

Nano-grain depletion in photon-dominated regions

Authors: T. Schirmer, N. Ysard, E. Habart, A. P. Jones, A. Abergel, L. Verstraete

Abstract: Context. Carbonaceous nano-grains play a fundamental role in the physico-chemistry of the interstellar medium (ISM) and especially of photon-dominated regions (PDRs). Their properties vary with the local physical conditions and affect the local chemistry and dynamics. Aims. We aim to highlight the evolution of carbonaceous nano-grains in three different PDRs and propose a scenario of dust evolutio… ▽ More Context. Carbonaceous nano-grains play a fundamental role in the physico-chemistry of the interstellar medium (ISM) and especially of photon-dominated regions (PDRs). Their properties vary with the local physical conditions and affect the local chemistry and dynamics. Aims. We aim to highlight the evolution of carbonaceous nano-grains in three different PDRs and propose a scenario of dust evolution as a response to the physical conditions. Methods. We used Spitzer/IRAC (3.6, 4.5, 5.8, and 8 $μみゅー$m) and Spitzer/MIPS (24 $μみゅー$m) together with Herschel/PACS (70 $μみゅー$m) to map dust emission in IC63 and the Orion Bar. To assess the dust properties, we modelled the dust emission in these regions using the radiative transfer code SOC together with the THEMIS dust model. Results. Regardless of the PDR, we find that nano-grains are depleted and that their minimum size is larger than in the diffuse ISM (DISM), which suggests that the mechanisms that lead nano-grains to be photo-destroyed are very efficient below a given critical size limit. The evolution of the nano-grain dust-to-gas mass ratio with both G0 and the effective temperature of the illuminating star indicates a competition between the nano-grain formation through the fragmentation of larger grains and nano-grain photo-destruction. We modelled dust collisions driven by radiative pressure with a classical 1D approach to show that this is a viable scenario for explaining nano-grain formation through fragmentation and, thus, the variations observed in nano-grain dust-to-gas mass ratios from one PDR to another. Conclusions. We find a broad variation in the nano-grain dust properties from one PDR to another, along with a general trend of nano-grain depletion in these regions. We propose a viable scenario of nano-grain formation through fragmentation of large grains due to radiative pressure-induced collisions. △ Less

Submitted 9 September, 2022; originally announced September 2022.

Journal ref: A&A 666, A49 (2022)

Influence of grain growth on CO2 ice spectroscopic profiles : Modelling for dense cores and disks

Authors: Emmanuel Dartois, Jennifer A. Noble, Nathalie Ysard, Karine Demyk, Marin Chabot

Abstract: Interstellar dust grain growth in dense clouds and protoplanetary disks, even moderate, affects the observed interstellar ice profiles as soon as a significant fraction of dust grains is in the size range close to the wave vector at the considered wavelength. The continuum baseline correction made prior to analysing ice profiles influences the subsequent analysis and hence the estimated ice compos… ▽ More Interstellar dust grain growth in dense clouds and protoplanetary disks, even moderate, affects the observed interstellar ice profiles as soon as a significant fraction of dust grains is in the size range close to the wave vector at the considered wavelength. The continuum baseline correction made prior to analysing ice profiles influences the subsequent analysis and hence the estimated ice composition, typically obtained by band fitting using thin film ice mixture spectra. We model the effect of grain growth on ice mantle spectroscopic profiles, focusing on CO2 to see how it can affect interstellar ice mantle spectral analysis and interpretation. Using the Discrete Dipole Approximation for Scattering and Absorption of Light, the mass absorption coefficients of distributions of grains composed of ellipsoidal silicate cores with water and carbon dioxide ice mantles are calculated. A few other ice mantle compositions are also calculated. We explore the size distribution evolution from dense clouds to simulate the first steps of grain growth up to three microns in size. The results are injected into RADMC-3D full scattering radiative transfer models of spherical clouds and protoplanetary disk templates to retrieve observable spectral energy distributions. We focus on calculated profile of the CO2 antisymmetric stretching mode ice band profile at 4.27 microns, a meaningful indicator of grain growth. The observed profiles toward dense cores with the Infrared space observatory and Akari satellites already showed profiles possibly indicative of moderate grain growth.The observation of protoplanetary disks at high inclination with the JWST should present distorted profiles that will put constraints on the extent of dust growth. The more evolved the dust size distribution, the more the extraction of the ice mantle composition will require both understanding and taking into account grain growth. △ Less

Submitted 19 July, 2022; originally announced July 2022.

Comments: 19 pages, 16 figures

Journal ref: A&A 666, A153 (2022)

Dust emissivity in resolved spiral galaxies

Authors: S. Bianchi, V. Casasola, E. Corbelli, F. Galliano, L. Magrini, A. Nersesian, F. Salvestrini, M. Baes, L. P. Cassara, C. J. R. Clark, I. De Looze, A. P. Jones, S. C. Madden, A. Mosenkov, N. Ysard

Abstract: Context: The far-infrared (FIR) and sub-millimeter (submm) emissivity of the Milky Way (MW) cirrus is an important benchmark for dust grain models. Dust masses in other galaxies are generally derived from the FIR/submm using the emission properties of these MW-calibrated models. Aims: We seek to derive the FIR/submm emissivity in nine nearby spiral galaxies to check its compatibility with MW cirru… ▽ More Context: The far-infrared (FIR) and sub-millimeter (submm) emissivity of the Milky Way (MW) cirrus is an important benchmark for dust grain models. Dust masses in other galaxies are generally derived from the FIR/submm using the emission properties of these MW-calibrated models. Aims: We seek to derive the FIR/submm emissivity in nine nearby spiral galaxies to check its compatibility with MW cirrus measurements. Methods: We obtained values of the emissivity at 70 to 500 um, using maps of dust emission from the Herschel satellite and of gas surface density from the THINGS and HERACLES surveys on a scale generally corresponding to 440 pc. We studied the variation of the emissivity with the surface brightness ratio I(250um)/I(500um), a proxy for the intensity of the interstellar radiation field heating the dust. Results: We find that the average value of the emissivity agrees with MW estimates for pixels sharing the same color as the cirrus, namely, for I(250um)/I(500um) = 4.5. For I(250um)/I(500um) > 5, the measured emissivity is instead up to a factor ~2 lower than predicted from MW dust models heated by stronger radiation fields. Regions with higher I(250um)/I(500um) are preferentially closer to the galactic center and have a higher overall (stellar+gas) surface density and molecular fraction. The results do not depend strongly on the adopted CO-to-molecular conversion factor and do not appear to be affected by the mixing of heating conditions. Conclusions: Our results confirm the validity of MW dust models at low density, but are at odds with predictions for grain evolution in higher density environments. If the lower-than-expected emissivity at high I(250um)/I(500um) is the result of intrinsic variations in the dust properties, it would imply an underestimation of the dust mass surface density of up to a factor ~2 when using current dust models. △ Less

Submitted 27 June, 2022; originally announced June 2022.

Comments: 16 pages, 11 figures, A&A accepted

Journal ref: A&A 664, A187 (2022)

Performance of the polarization leakage correction in the PILOT data

Authors: J-Ph. Bernard, A. Bernard, H. Roussel, I. Choubani, D. Alina, J. Aumont, A. Hughes, I. Ristorcelli, S. Stever, T. Matsumura S. Sugiyama, K. Komatsu, G. de Gasperis, K. Ferriere, V. Guillet, N. Ysard, P. Ade, P. de Bernardis, N. Bray, B. Crane, J. P. Dubois, M. Griffin, P. Hargrave, Y. Longval, S. Louvel, B. Maffei , et al. (11 additional authors not shown)

Abstract: The Polarized Instrument for Long-wavelength Observation of the Tenuous interstellar medium (PILOT) is a balloon-borne experiment that aims to measure the polarized emission of thermal dust at a wavelength of 240 um (1.2 THz). The PILOT experiment flew from Timmins, Ontario, Canada in 2015 and 2019 and from Alice Springs, Australia in April 2017. The in-flight performance of the instrument during… ▽ More The Polarized Instrument for Long-wavelength Observation of the Tenuous interstellar medium (PILOT) is a balloon-borne experiment that aims to measure the polarized emission of thermal dust at a wavelength of 240 um (1.2 THz). The PILOT experiment flew from Timmins, Ontario, Canada in 2015 and 2019 and from Alice Springs, Australia in April 2017. The in-flight performance of the instrument during the second flight was described in Mangilli et al. 2019. In this paper, we present data processing steps that were not presented in Mangilli et al. 2019 and that we have recently implemented to correct for several remaining instrumental effects. The additional data processing concerns corrections related to detector cross-talk and readout circuit memory effects, and leakage from total intensity to polarization. We illustrate the above effects and the performance of our corrections using data obtained during the third flight of PILOT, but the methods used to assess the impact of these effects on the final science-ready data, and our strategies for correcting them will be applied to all PILOT data. We show that the above corrections, and in particular that for the intensity to polarization leakage, which is most critical for accurate polarization measurements with PILOT, are accurate to better than 0.4 % as measured on Jupiter during flight#3. △ Less

Submitted 7 May, 2022; originally announced May 2022.

Spinning nano-carbon grains: Viable origin for anomalous microwave emission

Authors: Nathalie Ysard, Marc-Antoine Miville-Deschênes, Laurent Verstraete, Anthony Peter Jones

Abstract: Context. Excess microwave emission, commonly known as anomalous microwave emission (AME), is now routinely detected in the Milky Way. Although its link with the rotation of interstellar (carbonaceous) nano-grains seems to be relatively well established at cloud scales, large-scale observations show a lack of correlation between the different tracers of nano-carbons and AME, which has led the commu… ▽ More Context. Excess microwave emission, commonly known as anomalous microwave emission (AME), is now routinely detected in the Milky Way. Although its link with the rotation of interstellar (carbonaceous) nano-grains seems to be relatively well established at cloud scales, large-scale observations show a lack of correlation between the different tracers of nano-carbons and AME, which has led the community to question the viability of this link. Aims. Using ancillary data and spinning dust models for nano-carbons and nano-silicates, we explore the extent to which the AME that come out of the Galactic Plane might originate with one or another carrier. Methods. In contrast to previous large-scale studies, our method is not built on comparing the correlations of the different dust tracers with each other, but rather on comparing the poor correlations predicted by the models with observed correlations. This is based on estimates that are as realistic as possible of the gas ionisation state and grain charge as a function of the local radiation field and gas density. Results. First, nano-carbon dust can explain all the observations for medium properties, in agreement with the latest findings about the separation of cold and warm neutral medium in the diffuse interstellar medium. The dispersion in the observations can be accounted for with little variations in the dust size distribution, abundance, or electric dipole moment. Second, regardless of the properties and abundance of the nano-silicate dust we considered, spinning nano-silicates are excluded as the sole source of the AME. Third, the best agreement with the observations is obtained when the emission of spinning nano-carbons alone is taken into account. However, a marginal participation of nano-silicates in AME production cannot be excluded as long as their abundance does not exceed $Y_{\rm Si} \sim 1\%$. △ Less

Submitted 11 May, 2022; v1 submitted 3 May, 2022; originally announced May 2022.

Comments: accepted for publication by A&A

Journal ref: A&A 663, A65 (2022)

Radial distribution of the carbonaceous nano-grains in the protoplanetary disk around HD 169142

Authors: Marie Devinat, Émilie Habart, Éric Pantin, Nathalie Ysard, Anthony Jones, Lucas Labadie, Emmanuel Di Folco

Abstract: HD 169142 is part of the class of (pre-)transitional protoplanetary disks showing multiple carbon nanodust spectroscopic signatures (aromatic, aliphatic) dominating the infrared spectrum. Precise constraints on the spatial distribution and properties of carbonaceous dust particles are essential to understanding the physics of the disk. The HD 169142 disk is seen almost face-on and thus offers a un… ▽ More HD 169142 is part of the class of (pre-)transitional protoplanetary disks showing multiple carbon nanodust spectroscopic signatures (aromatic, aliphatic) dominating the infrared spectrum. Precise constraints on the spatial distribution and properties of carbonaceous dust particles are essential to understanding the physics of the disk. The HD 169142 disk is seen almost face-on and thus offers a unique opportunity to study the dust radial evolution. We investigate the spatial distribution and properties of the carriers of several dust aromatic emission features in the disk across a broad spatial range (10-200 AUえーゆー). We analysed imaging and spectroscopic observations in the 8-12 microns range from VLT/VISIR, as well as adaptive optics spectroscopic observations in the 3-4 microns range from VLT/NACO. The data probes the spatial evolution of the 3.3, 8.6, and 11.3 microns aromatic bands. To constrain the radial distribution of carbonaceous nano-grains, the observations were compared to models using The Heterogeneous dust Evolution Model for Interstellar Solids (THEMIS), integrated into the POLARIS radiative transfer code by calculating the thermal and stochastic heating of sub-micrometer dust grains. Our data show predominant nano-particle emission at all radii (resolution of about 0.1", 12 AUえーゆー at 3 microns and 0.3", 35 AUえーゆー at 10 microns) in the HD 169142 disk. This unambiguously shows that carbonaceous nano-grains dominate radiatively the infrared spectrum in most of the disk, as suggested by previous studies. In order to account for both VISIR and NACO emission maps, we show the need for aromatic particles distributed within the disk from the outermost regions to a radius of 20 AUえーゆー, corresponding to the outer limit of the inner cavity derived from previous observations. In the inner cavity, these aromatic particles might be present but their abundance would then be significantly decreased. △ Less

Submitted 30 April, 2022; originally announced May 2022.

Comments: 11 pages, 8 figures

Journal ref: A&A 663, A151 (2022)

Searching for Anomalous Microwave Emission in nearby galaxies. K-band observations with the Sardinia Radio Telescope

Authors: S. Bianchi, M. Murgia, A. Melis, V. Casasola, F. Galliano, F. Govoni, A. P. Jones, S. C. Madden, R. Paladino, F. Salvestrini, E. M. Xilouris, N. Ysard

Abstract: We observed four nearby spiral galaxies (NGC 3627, NGC 4254, NGC 4736 and NGC 5055) in the K band with the 64-m Sardinia Radio Telescope, with the aim of detecting the Anomalous Microwave Emission (AME), a radiation component presumably due to spinning dust grains, observed so far in the Milky Way and in a handful of other galaxies only (most notably, M 31). We mapped the galaxies at 18.6 and 24.6… ▽ More We observed four nearby spiral galaxies (NGC 3627, NGC 4254, NGC 4736 and NGC 5055) in the K band with the 64-m Sardinia Radio Telescope, with the aim of detecting the Anomalous Microwave Emission (AME), a radiation component presumably due to spinning dust grains, observed so far in the Milky Way and in a handful of other galaxies only (most notably, M 31). We mapped the galaxies at 18.6 and 24.6 GHz and studied their global photometry together with other radio-continuum data from the literature, in order to find AME as emission in excess of the synchrotron and thermal components. We only find upper limits for AME. These non-detections, and other upper limits in the literature, are nevertheless consistent with the average AME emissivity from the few detections: it is $εいぷしろん^\mathrm{AME}_{\mathrm{30~GHz}} = 2.4\pm0.4 \times 10^{-2}$ MJy sr$^{-1}$ (M$_\odot$ pc$^{-2}$)$^{-1}$ in units of dust surface density (equivalently, $1.4\pm0.2 \times 10^{-18}$ Jy sr$^{-1}$ (H cm$^{-2}$)$^{-1}$ in units of H column density). We finally suggest to search for AME in quiescent spirals with relatively low radio luminosity, such as M~31. △ Less

Submitted 24 January, 2022; originally announced January 2022.

Comments: A&A Letter, accepted

Journal ref: A&A 658, L8 (2022)

PDRs4All: A JWST Early Release Science Program on radiative feedback from massive stars

Authors: Olivier Berné, Émilie Habart, Els Peeters, Alain Abergel, Edwin A. Bergin, Jeronimo Bernard-Salas, Emeric Bron, Jan Cami, Stéphanie Cazaux, Emmanuel Dartois, Asunción Fuente, Javier R. Goicoechea, Karl D. Gordon, Yoko Okada, Takashi Onaka, Massimo Robberto, Markus Röllig, Alexander G. G. M. Tielens, Silvia Vicente, Mark G. Wolfire, Felipe Alarcon, C. Boersma, Ameélie Canin, Ryan Chown, Daniel Dicken , et al. (112 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 13 January, 2022; originally announced January 2022.

Comments: Submitted to PASP

A derivation of nano-diamond optical constants: Here be nano-diamonds

Authors: A. P. Jones, Nathalie Ysard

Abstract: % context {Nano-diamonds are an enticing and enigmatic dust component yet their origin is still unclear. They have been unequivocally detected in only a few astronomical objects, yet they are the most abundant of the pre-solar grains, both in terms of mass and number.} %aims {Our goal is to derive a viable set of nano-diamond optical constants and optical properties to enable their modelling in an… ▽ More % context {Nano-diamonds are an enticing and enigmatic dust component yet their origin is still unclear. They have been unequivocally detected in only a few astronomical objects, yet they are the most abundant of the pre-solar grains, both in terms of mass and number.} %aims {Our goal is to derive a viable set of nano-diamond optical constants and optical properties to enable their modelling in any type of astrophysical object where, primarily, the local (inter)stellar radiation field is well-determined.} % methods {The complex indices of refraction, $m(n,k)$, of nano-diamonds, constrained by available laboratory measurements, were calculated as a function of size, surface hydrogenation, and internal (dis)order, using the THEMIS a-C(:H) methodology optEC$_{\rm (s)}$(a).} % results {To demonstrate the utility of the optical properties (the efficiency factors $Q_{\rm ext}$, $Q_{\rm sca}$, and $Q_{\rm abs}$), calculated using the derived $m(n,k)$ data, we show that nano-diamonds could be abundant in the interstellar medium (ISM) and yet remain undetectable there.} % conclusions {The derived optical constants provide a means to explore the existence and viability of nano-diamonds in a wide range of astronomical sources. Here we show that up to a few percent of the available carbon budget could be hidden in the form of nano-diamonds in the diffuse ISM, in abundances comparable to the pre-solar nano-diamond abundances in primitive meteorites.} △ Less

Submitted 8 November, 2021; originally announced November 2021.

Comments: 17 pages, 8 figures, 4 tables

Journal ref: A&A 657, A128 (2022)

The mm-to-cm SED of spiral galaxies. Synergies between NIKA2 and SRT instruments

Authors: S. Bianchi, M. Murgia, A. Melis, V. Casasola, M. Galametz, F. Galliano, F. Govoni, A. Jones, S. Madden, R. Paladino, E. Xilouris, N. Ysard

Abstract: The mm-to-cm range of the Spectral Energy Distribution of spiral galaxies remains largely unexplored. Its coverage is required to disentangle the contribution of dust emission, free-free and synchrotron radiation and can provide constraints on dust models, star-formation rates and ISM properties. We present the case for a synergy between NIKA2 observations of nearby spirals and those from planned… ▽ More The mm-to-cm range of the Spectral Energy Distribution of spiral galaxies remains largely unexplored. Its coverage is required to disentangle the contribution of dust emission, free-free and synchrotron radiation and can provide constraints on dust models, star-formation rates and ISM properties. We present the case for a synergy between NIKA2 observations of nearby spirals and those from planned and current instrumentation at the Sardinia Radio Telescope, and report on a pilot K-band program to search for Anomalous Microwave Emission, an elusive emission component which is presumably related to dust. △ Less

Submitted 2 November, 2021; originally announced November 2021.

Comments: To appear in the Proceedings of the International Conference entitled "mm Universe @ NIKA2", Rome (Italy), June 2021, EPJ Web of conferences

The formation of planetary systems with SPICA

Authors: I. Kamp, M. Honda, H. Nomura, M. Audard, D. Fedele, L. B. F. M. Waters, Y. Aikawa, A. Banzatti, J. E. Bowey, M. Bradford, C. Dominik, K. Furuya, E. Habart, D. Ishihara, D. Johnstone, G. Kennedy, M. Kim, Q. Kral, S. P. Lai, B. Larsson, M. McClure, A. Miotello, M. Momose, T. Nakagawa, D. Naylor , et al. (16 additional authors not shown)

Abstract: In this era of spatially resolved observations of planet forming disks with ALMA and large ground-based telescopes such as the VLT, Keck and Subaru, we still lack statistically relevant information on the quantity and composition of the material that is building the planets, such as the total disk gas mass, the ice content of dust, and the state of water in planetesimals. SPICA is an infrared spac… ▽ More In this era of spatially resolved observations of planet forming disks with ALMA and large ground-based telescopes such as the VLT, Keck and Subaru, we still lack statistically relevant information on the quantity and composition of the material that is building the planets, such as the total disk gas mass, the ice content of dust, and the state of water in planetesimals. SPICA is an infrared space mission concept developed jointly by JAXA and ESA to address these questions. The key unique capabilities of SPICA that enable this research are (1) the wide spectral coverage 10-220 micron, (2) the high line detection sensitivity of (1-2) 10-19 W m-2 with R~2000-5000 in the far-IR (SAFARI) and 10-20 W m-2 with R~29000 in the mid-IR (SMI, spectrally resolving line profiles), (3) the high far-IR continuum sensitivity of 0.45 mJy (SAFARI), and (4) the observing efficiency for point source surveys. This paper details how mid- to far-IR infrared spectra will be unique in measuring the gas masses and water/ice content of disks and how these quantities evolve during the planet forming period. These observations will clarify the crucial transition when disks exhaust their primordial gas and further planet formation requires secondary gas produced from planetesimals. The high spectral resolution mid-IR is also unique for determining the location of the snowline dividing the rocky and icy mass reservoirs within the disk and how the divide evolves during the build-up of planetary systems. Infrared spectroscopy (mid- to far-IR) of key solid state bands is crucial for assessing whether extensive radial mixing, which is part of our Solar System history, is a general process occurring in most planetary systems and whether extrasolar planetesimals are similar to our Solar System comets/asteroids. ... (abbreviated) △ Less

Submitted 25 June, 2021; originally announced June 2021.

Comments: accepted for publication in PASA

First MATISSE L-band observations of HD 179218. Is the inner 10 au region rich in carbon dust particles?

Authors: E. Kokoulina, A. Matter, B. Lopez, E. Pantin, N. Ysard, G. Weigelt, E. Habart, J. Varga, A. Jones, A. Meilland, E. Dartois, L. Klarmann, J. -C. Augereau, R. van Boekel, M. Hogerheijde, G. Yoffe, L. B. F. M. Waters, C. Dominik, W. Jaffe, F. Millour, Th. Henning, K. -H. Hofmann, D. Schertl, S. Lagarde, R. G. Petrov , et al. (36 additional authors not shown)

Abstract: Carbon is one of the most abundant components in the Universe. While silicates have been the main focus of solid phase studies in protoplanetary discs (PPDs), little is known about the solid carbon content especially in the planet-forming regions ($\sim $0.1 to 10 au). Fortunately, several refractory carbonaceous species present C-H bonds (such as hydrogenated nano-diamond and amorphous carbon as… ▽ More Carbon is one of the most abundant components in the Universe. While silicates have been the main focus of solid phase studies in protoplanetary discs (PPDs), little is known about the solid carbon content especially in the planet-forming regions ($\sim $0.1 to 10 au). Fortunately, several refractory carbonaceous species present C-H bonds (such as hydrogenated nano-diamond and amorphous carbon as well as polycyclic aromatic hydrocarbons (PAHs)), which generate infrared (IR) features that can be used to trace the solid carbon reservoirs. The new mid-IR instrument MATISSE, installed at the Very Large Telescope Interferometer (VLTI), can spatially resolve the inner regions ($\sim$ 1 to 10 au) of PPDs and locate, down to the au-scale, the emission coming from carbon grains. Our aim is to provide a consistent view on the radial structure, down to the au-scale, as well as basic physical properties and the nature of the material responsible for the IR continuum emission in the inner disk region around HD 179218. We implemented a temperature-gradient model to interpret the disk IR continuum emission, based on a multiwavelength dataset comprising a broadband spectral energy distribution (SED) and VLTI H-, L-, and N-bands interferometric data obtained in low spectral resolution. Then, we added a ring-like component, representing the carbonaceous L-band features-emitting region, to assess its detectability in future higher spectral resolution observations employing mid-IR interferometry. △ Less

Submitted 29 July, 2021; v1 submitted 24 June, 2021; originally announced June 2021.

Journal ref: A&A 652, A61 (2021)

Influence of the nano-grain depletion in photon-dominated regions: Application to the gas physics and chemistry in the Horsehead

Authors: T. Schirmer, E. Habart, N. Ysard, E. Bron, J. Le Bourlot, L. Verstraete, A. Abergel, A. P. Jones, E. Roueff, F. Le Petit

Abstract: The large disparity in physical conditions from the diffuse interstellar medium (ISM) to denser clouds such as photon-dominated regions (PDRs) triggers an evolution of the dust properties (i.e. composition, size, and shape). The gas physics and chemistry are tightly connected to these dust properties and are therefore affected by dust evolution and especially the nano-grain depletion in the outer… ▽ More The large disparity in physical conditions from the diffuse interstellar medium (ISM) to denser clouds such as photon-dominated regions (PDRs) triggers an evolution of the dust properties (i.e. composition, size, and shape). The gas physics and chemistry are tightly connected to these dust properties and are therefore affected by dust evolution and especially the nano-grain depletion in the outer irradiated part of PDRs. We highlight the influence of nano-grain depletion on the gas physics and chemistry in the Horsehead nebula, a prototypical PDR. We used a model for atomic and molecular gas in PDRs, the Meudon PDR code, using diffuse ISM-like dust and Horsehead-like dust to study the influence of nano-grain depletion on the gas physics and chemistry, focusing on the impact on photoelectric heating and H2 formation and, therefore, on the H2 gas lines. We find that nano-grain depletion in the Horsehead strongly affects gas heating through the photoelectric effect and thus the gas temperature and the H2 formation, hence the H -> H2 position. Consequently, the first four pure rotational lines of H2 (e.g. 0-0 S(0), S(1), S(2), and S(3)) vary by a factor of 2 to 14. The 0-0 S(3) line that is often underestimated in models is underestimated even more when taking nano-grain depletion into account due to the decrease in gas heating through the photoelectric effect. This strongly suggests that our understanding of the excitation of H2 and/or of heating processes in the Horsehead, and more generally in PDRs, is still incomplete. Nano-grain depletion in the outer part of the Horsehead has a strong influence on several gas tracers that will be prominent in JWST observations of irradiated clouds. We therefore need to take this depletion into account in order to improve our understanding of the Horsehead, and more generally PDRs, and to contribute to the optimal scientific return of the mission. △ Less

Submitted 16 April, 2021; originally announced April 2021.

Comments: Accepted in A&A, comments welcome

Journal ref: A&A 649, A148 (2021)

Multi-wavelength observations and modelling of a quiescent cloud LDN1512

Authors: M. Saajasto, M. Juvela, C. Lefevre, L. Pagani, N. Ysard

Abstract: Light scattering at near-infrared wavelengths has been used to study the optical properties of the interstellar dust grains, but these studies are limited by the assumptions on the strength of the radiation field. On the other hand, thermal dust emission can be used to constrain the properties of the radiation field, although this is hampered by uncertainty about the dust emissivity. We test if cu… ▽ More Light scattering at near-infrared wavelengths has been used to study the optical properties of the interstellar dust grains, but these studies are limited by the assumptions on the strength of the radiation field. On the other hand, thermal dust emission can be used to constrain the properties of the radiation field, although this is hampered by uncertainty about the dust emissivity. We test if current dust models allow us to model a molecular cloud simultaneously in the near infrared (NIR) and far infrared (FIR) wavelengths and compare the results with observations. Our aim is to place constraints on the properties of the dust grains and the strength of the radiation field. We present computations of dust emission and scattered light of a quiescent molecular cloud LDN1512. We construct radiative transfer models for LDN1512 that include an anisotropic radiation field and a three-dimensional cloud model. We are able to reproduce the observed FIR observations, with a radiation field derived from the DIRBE observations, with all of the tested dust models. However, with the same density distribution and the assumed radiation field, the models fail to reproduce the observed NIR scattering in all cases except for models that take into account dust evolution via coagulation and mantle formation. We find that the column densities derived from our radiative transfer modelling can differ by a factor of up to two, compared to the column densities derived from the observations with modified blackbody fits. The discrepancy in the column densities is likely caused because of temperature difference between a modified blackbody fit and the real spectra. We show that the observed dust emission can be reproduced with several different assumptions about the properties of the dust grains. However, in order to reproduce the observed scattered surface brightness dust evolution must be taken into account. △ Less

Submitted 28 January, 2021; originally announced January 2021.

Comments: 16 pages, 16 figures. Recommended for publication to A&A

A Nearby Galaxy Perspective on Dust Evolution. Scaling relations and constraints on the dust build-up in galaxies with the DustPedia and DGS samples

Authors: Frédéric Galliano, Angelos Nersesian, Simone Bianchi, Ilse De Looze, Sambit Roychowdhury, Maarten Baes, Viviana Casasola, Letizia, P. Cassará, Wouter Dobbels, Jacopo Fritz, Maud Galametz, Anthony P. Jones, Suzanne C. Madden, Aleksandr Mosenkov, Emmanuel M. Xilouris, Nathalie Ysard

Abstract: Methods. We have modelled a sample of ~800 nearby galaxies, spanning a wide range of metallicity, gas fraction, specific star formation rate and Hubble stage. We have derived the dust properties of each object from its spectral energy distribution. Through an additional level of analysis, we have inferred the timescales of dust condensation in core-collapse supernova ejecta, grain growth in cold c… ▽ More Methods. We have modelled a sample of ~800 nearby galaxies, spanning a wide range of metallicity, gas fraction, specific star formation rate and Hubble stage. We have derived the dust properties of each object from its spectral energy distribution. Through an additional level of analysis, we have inferred the timescales of dust condensation in core-collapse supernova ejecta, grain growth in cold clouds and dust destruction by shock waves. Throughout this paper, we have adopted a hierarchical Bayesian approach, resulting in a single large probability distribution of all the parameters of all the galaxies, to ensure the most rigorous interpretation of our data. Results. We confirm the drastic evolution with metallicity of the dust-to-metal mass ratio (by two orders of magnitude), found by previous studies. We show that dust production by core-collapse supernovae is efficient only at very low-metallicity, a single supernova producing on average less than ~0.03 Msun/SN of dust. Our data indicate that grain growth is the dominant formation mechanism at metallicity above ~1/5 solar, with a grain growth timescale shorter than ~50 Myr at solar metallicity. Shock destruction is relatively efficient, a single supernova clearing dust on average in at least ~1200 Msun/SN of gas. These results are robust when assuming different stellar initial mass functions. In addition, we show that early-type galaxies are outliers in several scaling relations. This feature could result from grain thermal sputtering in hot X-ray emitting gas, an hypothesis supported by a negative correlation between the dust-to-stellar mass ratio and the X-ray photon rate per grain. Finally, we confirm the well-known evolution of the aromatic-feature-emitting grain mass fraction as a function of metallicity and interstellar radiation field intensity. Our data indicate the relation with metallicity is significantly stronger. △ Less

Submitted 2 January, 2021; originally announced January 2021.

Comments: 64 pages, 36 figures, accepted for publication in Astronomy & Astrophysics

Spatial distribution of the aromatic and aliphatic carbonaceous nano-grain features in the protoplanetary disk around HD 100546

Authors: Emilie Habart, Thomas Boutéraon, Robert Brauer, Nathalie Ysard, Eric Pantin, Antoine Marchal, Anthony P. Jones

Abstract: Carbonaceous nano-grains are present at the surface of protoplanetary disks around Herbig Ae/Be stars, where most of the central star UV energy is dissipated. Efficiently coupled to the gas, nano-grains are able to trace the disk outer flaring parts, and possibly the gaps from which the larger grains are missing. We examine the spatial distribution and evolution of the nano-dust emission in the (p… ▽ More Carbonaceous nano-grains are present at the surface of protoplanetary disks around Herbig Ae/Be stars, where most of the central star UV energy is dissipated. Efficiently coupled to the gas, nano-grains are able to trace the disk outer flaring parts, and possibly the gaps from which the larger grains are missing. We examine the spatial distribution and evolution of the nano-dust emission in the (pre-)transitional disk HD100546 that shows annular gaps, rings, and spirals, and reveals rich carbon nano-dust spectroscopic signatures (aromatic, aliphatic) in a wide spatial range (~20-200au). We analyse adaptive optics spectroscopic observations from 3 to 4um and imaging and spectroscopic observations from 8 to 12um. We compare the data to model predictions using the THEMIS model with the radiative transfer code POLARIS calculating heating of micro- and nanometric dust grains for a given disk structure. The aromatic features at 3.3, 8.6 and 11.3um, as well as, the aliphatic ones from 3.4 to 3.5um are spatially extended with band morphologies dependong on local physical conditions. The aliphatic-to-aromatic band ratio 3.4/3.3 increases with the distance from the star suggesting UV processing. In the 8-12um observed spectra, features characteristic of aromatic particles and crystalline silicates are detected with their relative contribution changing with distance to the star. The model predicts that the features and adjacent continuum are due to different combinations of grain sub-populations, with a dependence on the UV field intensity. Shorter wavelength features are dominated by the smallest grains (< 0.7nm) throughout the disk, while at longer wavelengths what dominates the emission close to the star is a mix of several grain populations, and far away from the star is the largest nano-grain population. △ Less

Submitted 27 November, 2020; originally announced November 2020.

Comments: accepted for publication in A&A, abridged abstract to comply with arxiv requirements

Journal ref: A&A 649, A84 (2021)

Dust evolution: going beyond the empirical

Authors: Nathalie Ysard

Abstract: A key element when modeling dust in any astrophysical environment is a self-consistent treatment of the evolution of the dust material properties (size distribution, chemical composition and structure) as they react to and adjust to the local radiation field intensity and hardness and to the gas density and dynamics. The best way to achieve this goal is to anchore as many model parameters as possi… ▽ More A key element when modeling dust in any astrophysical environment is a self-consistent treatment of the evolution of the dust material properties (size distribution, chemical composition and structure) as they react to and adjust to the local radiation field intensity and hardness and to the gas density and dynamics. The best way to achieve this goal is to anchore as many model parameters as possible to laboratory data. In this paper, I present two examples to illustrate how outstanding questions in dust modeling have been/are being moved forward by recent advances in laboratory astrophysics and what laboratory data are still needed to further advance dust evolution models. △ Less

Submitted 27 November, 2020; originally announced November 2020.

Dust emission, extinction, and scattering in LDN 1642

Authors: Mika Juvela, Sharma Neha, Emma Mannfors, Mika Saajasto, Nathalie Ysard, Veli-Matti Pelkonen

Abstract: We study the near-infrared (NIR) scattering in LDN 1642, its correlation with the cloud structure, and the ability of dust models to simultaneously explain sub-millimetre emission, NIR extinction, and NIR scattering. We use observations from the HAWK-I instrument to measure the NIR surface brightness and extinction. These are compared with Herschel data on dust emission and, with radiative trans… ▽ More We study the near-infrared (NIR) scattering in LDN 1642, its correlation with the cloud structure, and the ability of dust models to simultaneously explain sub-millimetre emission, NIR extinction, and NIR scattering. We use observations from the HAWK-I instrument to measure the NIR surface brightness and extinction. These are compared with Herschel data on dust emission and, with radiative transfer modelling, with predictions calculated for different dust models. We find an optical depth ratio $τたう(250\,μみゅー{\rm m})/τたう(J)\approx 10^{-3}$, confirming earlier findings of high sub-millimetre emissivity. The relationships between the column density derived from dust emission and the NIR colour excesses is linear and consistent with the standard NIR extinction curve. The extinction peaks at $A_J=2.6\,$mag, the NIR surface brightness remaining correlated with $N({\rm H}_2)$ without saturation. Radiative transfer models can fit the sub-millimetre data with any of the tested dust models. However, these predict a NIR extinction that is higher and a NIR surface brightness that is lower than in observations. If the dust sub-millimetre emissivity is rescaled to the observed value of $τたう(250\,μみゅー{\rm m})/τたう(J)$, dust models with high NIR albedo can reach the observed level of NIR surface brightness. The NIR extinction of the models tends to be higher than directly measured, which is reflected in the shape of the NIR surface brightness spectra. The combination of emission, extinction, and scattering measurements provides strong constraints on dust models. The observations of LDN 1642 indicate clear dust evolution, including a strong increase in the sub-millimetre emissivity, not yet fully explained by the current dust models. △ Less

Submitted 26 September, 2020; originally announced September 2020.

Comments: 16 pages, 16 figures, submitted to A&A

High-resolution, 3D radiative transfer modelling IV. AGN-powered dust heating in NGC 1068

Authors: S. Viaene, A. Nersesian, J. Fritz, S. Verstocken, M. Baes, S. Bianchi, V. Casasola, L. Cassarà, C. Clark, J. Davies, I. De Looze, P. De Vis, W. Dobbels, M. Galametz, F. Galliano, A. Jones, S. Madden, A. Mosenkov, A. Trcka, E. M. Xilouris, N. Ysard

Abstract: Dust emission, an important diagnostic of star formation and ISM mass throughout the Universe, can be powered by sources unrelated to ongoing star formation. In the framework of the DustPedia project we have set out to disentangle the radiation of the ongoing star formation from that of the older stellar populations. This is done through detailed, 3D radiative transfer simulations of face-on spira… ▽ More Dust emission, an important diagnostic of star formation and ISM mass throughout the Universe, can be powered by sources unrelated to ongoing star formation. In the framework of the DustPedia project we have set out to disentangle the radiation of the ongoing star formation from that of the older stellar populations. This is done through detailed, 3D radiative transfer simulations of face-on spiral galaxies. In this particular study, we focus on NGC 1068, which contains an active galactic nucleus (AGN). The effect of diffuse dust heating by AGN (beyond the torus) was so far only investigated for quasars. This additional dust heating source further contaminates the broadband fluxes on which classic galaxy modelling tools rely to derive physical properties. We aim to fit a realistic model to the observations of NGC 1068 and quantify the contribution of the several dust heating sources. Our model is able to reproduce the global spectral energy distribution of the galaxy. It matches the resolved optical and infrared images fairly well, but deviates in the UV and the submm. We find a strong wavelength dependency of AGN contamination to the broadband fluxes. It peaks in the MIR, drops in the FIR, but rises again at submm wavelengths. We quantify the contribution of the dust heating sources in each 3D dust cell and find a median value of 83% for the star formation component. The AGN contribution is measurable at the percentage level in the disc, but quickly increases in the inner few 100 pc, peaking above 90%. This is the first time the phenomenon of an AGN heating the diffuse dust beyond its torus is quantified in a nearby star-forming galaxy. NGC 1068 only contains a weak AGN, meaning this effect can be stronger in galaxies with a more luminous AGN. This could significantly impact the derived star formation rates and ISM masses for such systems. △ Less

Submitted 4 May, 2020; originally announced May 2020.

Comments: 18 pages, 11 figures, accepted for publication in A&A

High-resolution, 3D radiative transfer modelling III. The DustPedia barred galaxies

Authors: Angelos Nersesian, Sam Verstocken, Sebastien Viaene, Maarten Baes, Emmanuel M. Xilouris, Simone Bianchi, Viviana Casasola, Christopher J. R. Clark, Jonathan I. Davies, Ilse De Looze, Pieter De Vis, Wouter Dobbels, Jacopo Fritz, Maud Galametz, Frederic Galliano, Anthony P. Jones, Suzanne C. Madden, Aleksandr V. Mosenkov, Ana Trcka, Nathalie Ysard

Abstract: Context: Dust in late-type galaxies in the local Universe is responsible for absorbing approximately one third of the energy emitted by stars. It is often assumed that dust heating is mainly attributable to the absorption of UV and optical photons emitted by the youngest (<= 100 Myr) stars. Consequently, thermal re-emission by dust at FIR wavelengths is often linked to the star-formation activity… ▽ More Context: Dust in late-type galaxies in the local Universe is responsible for absorbing approximately one third of the energy emitted by stars. It is often assumed that dust heating is mainly attributable to the absorption of UV and optical photons emitted by the youngest (<= 100 Myr) stars. Consequently, thermal re-emission by dust at FIR wavelengths is often linked to the star-formation activity of a galaxy. However, several studies argue that the contribution to dust heating by much older stars might be more significant. Advances in radiation transfer (RT) simulations finally allow us to actually quantify the heating mechanisms of diffuse dust by the stellar radiation field. Aims: As one of the main goals in the DustPedia project, we have constructed detailed 3D stellar and dust RT models for nearby galaxies. We analyse the contribution of the different stellar populations to the dust heating in four face-on barred galaxies: NGC1365, M83, M95, and M100. We aim to quantify the fraction directly related to young stars, both globally and on local scales, and to assess the influence of the bar on the heating fraction. Results: We derive global attenuation laws for each galaxy and confirm that galaxies of high sSFR have shallower attenuation curves and weaker UV bumps. On average, 36.5% of the bolometric luminosity is absorbed by dust. We report a clear effect of the bar structure on the radial profiles of the dust-heating fraction by the young stars, and the dust temperature. We find that the young stars are the main contributors to the dust heating, donating, on average ~59% of their luminosity to this purpose throughout the galaxy. This dust-heating fraction drops to ~53% in the bar region and ~38% in the bulge region where the old stars are the dominant contributors to the dust heating. We also find a strong link between the heating fraction by the young stars and the sSFR. △ Less

Submitted 7 April, 2020; originally announced April 2020.

Dust evolution across the Horsehead Nebula

Authors: T. Schirmer, A. Abergel, L. Verstraete, N. Ysard, M. Juvela, A. P. Jones, E. Habart

Abstract: Micro-physical processes on interstellar dust surfaces are tightly connected to dust properties (i.e. dust composition, size and shape) and play a key role in numerous phenomena in the interstellar medium (ISM). The large disparity in physical conditions (i.e. density, gas temperature) in the ISM triggers an evolution of dust properties. The analysis of how dust evolves with the physical condition… ▽ More Micro-physical processes on interstellar dust surfaces are tightly connected to dust properties (i.e. dust composition, size and shape) and play a key role in numerous phenomena in the interstellar medium (ISM). The large disparity in physical conditions (i.e. density, gas temperature) in the ISM triggers an evolution of dust properties. The analysis of how dust evolves with the physical conditions is a stepping-stone towards a more thorough understanding of interstellar dust. The aim of this paper is to highlight dust evolution in the Horsehead Nebula PDR region. We use Spitzer/IRAC (3.6, 4.5, 5.8 and 8 μみゅーm), Spitzer/MIPS (24 μみゅーm) together with Herschel/PACS (70 and 160 μみゅーm) and Herschel/SPIRE (250, 350 and 500 μみゅーm) to map the spatial distribution of dust in the Horsehead over the entire emission spectral range. We model dust emission and scattering using the THEMIS interstellar dust model together with the 3D radiative transfer code SOC. We find that the nano-grains dust-to-gas ratio in the irradiated outer part of the Horsehead is 6 to 10 times lower than in the diffuse ISM. Their minimum size is 2 to 2.25 times larger than in the diffuse ISM and the power-law exponent of their size distribution, 1.1 to 1.4 times lower than in the diffuse ISM. Regarding the denser part of the Horsehead, it is necessary to use evolved grains (i.e. aggregates, with or without an ice mantle). It is not possible to explain the observations using grains from the diffuse medium. We therefore propose the following scenario to explain our results. In the outer part of the Horsehead, all the nano-grains have not yet had time to re-form completely through photo-fragmentation of aggregates and the smallest of the nano-grains that are sensitive to the radiation field are photo-destroyed. In the inner part of the Horsehead, grains most likely consist of multi-compositional, mantled aggregates. △ Less

Submitted 12 March, 2020; originally announced March 2020.

Journal ref: A&A 639, A144 (2020)

Dust emissivity and absorption cross section in DustPedia late-type galaxies

Authors: S. Bianchi, V. Casasola, M. Baes, C. J. R. Clark, E. Corbelli, J. I. Davies, I. De Looze, P. De Vis, W. Dobbels, M. Galametz, F. Galliano, A. P. Jones, S. C. Madden, L. Magrini, A. Mosenkov, A. Nersesian, S. Viaene, E. M. Xilouris, N. Ysard

Abstract: Aims: We compare the far-infrared to sub-millimetre dust emission properties measured in high Galactic latitude cirrus with those determined in a sample of 204 late-type DustPedia galaxies. The aim is to verify if it is appropriate to use Milky Way dust properties to derive dust masses in external galaxies. Methods: We used Herschel observations and atomic and molecular gas masses to estimate the… ▽ More Aims: We compare the far-infrared to sub-millimetre dust emission properties measured in high Galactic latitude cirrus with those determined in a sample of 204 late-type DustPedia galaxies. The aim is to verify if it is appropriate to use Milky Way dust properties to derive dust masses in external galaxies. Methods: We used Herschel observations and atomic and molecular gas masses to estimate the disc-averaged dust emissivity at 250 micrometres, and from this, the absorption cross section per H atom and per dust mass. The emissivity requires one assumption, which is the CO-to-H_2 conversion factor, and the dust temperature is additionally required for the absorption cross section per H atom; yet another constraint on the dust-to-hydrogen ratio D/H, depending on metallicity, is required for the absorption cross section dust mass. Results: We find epsilon(250) = 0.82 +/- 0.07 MJy sr^-1 (1E20 H cm^-2)^-1 for galaxies with 4 < F(250)/F(500) < 5. This depends only weakly on the adopted CO-to-H_2 conversion factor. The value is almost the same as that for the Milky Way at the same colour ratio. Instead, for F(250)/F(500) > 6, epsilon(250) is lower than predicted by its dependence on the heating conditions. The reduction suggests a variation in dust emission properties for spirals of earlier type, higher metallicity, and with a higher fraction of molecular gas. When the standard emission properties of Galactic cirrus are used for these galaxies, their dust masses might be underestimated by up to a factor of two. Values for the absorption cross sections at the Milky Way metallicity are also close to those of the cirrus. Mild trends of the absorption cross sections with metallicity are found, although the results depend on the assumptions made. △ Less

Submitted 27 September, 2019; originally announced September 2019.

Comments: 13 pages, 6 figures, 2 tables, A&A accepted

Journal ref: A&A 631, A102 (2019)

From grains to pebbles: the influence of size distribution and chemical composition on dust emission properties

Authors: N. Ysard, M. Koehler, I. Jimenez-Serra, A. P. Jones, L. Verstraete

Abstract: The size and composition of dust grains are critical in setting the dynamical, physical and chemical evolution of the media in which they are present. Thanks to facilities such as ALMA and in the future the SKA, their thermal emission in the (sub)mm to cm has become a convenient way to trace grain properties. Our aim is to understand the influence of the grain composition and size distribution on… ▽ More The size and composition of dust grains are critical in setting the dynamical, physical and chemical evolution of the media in which they are present. Thanks to facilities such as ALMA and in the future the SKA, their thermal emission in the (sub)mm to cm has become a convenient way to trace grain properties. Our aim is to understand the influence of the grain composition and size distribution on the shape of their SED in dense ISM regions such as molecular clouds, prestellar cores, YSOs and protoplanetary discs. Starting from the optical constants defined in the THEMIS model for amorphous hydrogenated carbon and silicate grains in addition to water ice, we define 6 material mixtures representative of the expected dust composition in dense ISM regions. The optical properties of 0.01 micron to 10 cm grains are then calculated with effective medium and Mie theories. The corresponding SEDs are calculated for isolated clouds either externally heated by the ISRF alone or in addition to an internal source. The 3 main outcomes of this study are as follows. First, the dust mass absorption coefficient strongly depends on both its composition and size distribution potentially leading to errors in dust mass estimates by factors up to 3 and 20, respectively. Second, it appears almost impossible to retrieve the grain composition from the (sub)mm to cm SED shape alone as its spectral index for lambda > 3 mm does not depend on composition. Third, using the true dust opacity spectral index to estimate grain sizes may lead to erroneous findings as the observed spectral index can be highly modified by the dust temperature distribution along the line-of-sight, which depends on the specific heating source and on the geometry of the studied region. Based on the interpretation of only the spectral shape of (sub)mm to cm SEDs, the determination of the dust masses, compositions and sizes are highly uncertain. △ Less

Submitted 5 November, 2019; v1 submitted 11 September, 2019; originally announced September 2019.

Comments: accepted by A&A, abstract modified to match arXiv standards

Journal ref: A&A 631, A88 (2019)

The First Maps of $κかっぱ_{d}$ -- the Dust Mass Absorption Coefficient -- in Nearby Galaxies, with DustPedia

Authors: Christopher J. R. Clark, Pieter De Vis, Maarten Baes, Simone Bianchi, Viviana Casasola, Letizia P. Cassarà, Jonathan I. Davies, Wouter Dobbels, Sofia Lianou, Ilse De Looze, Ruth Evans, Maud Galametz, Frederic Galliano, Anthony P. Jones, Suzanne C. Madden, Alexander V. Mosenkov, Sam Verstocken, Sébastien Viaene, E. Manolis Xilouris, Nathalie Ysard

Abstract: The dust mass absorption coefficient, $κかっぱ_{d}$, is the conversion function used to infer physical dust masses from observations of dust emission. However, it is notoriously poorly constrained, and it is highly uncertain how it varies, either between or within galaxies. Here we present the results of a proof-of concept study, using the DustPedia data for two nearby face-on spiral galaxies M74 (NGC 6… ▽ More The dust mass absorption coefficient, $κかっぱ_{d}$, is the conversion function used to infer physical dust masses from observations of dust emission. However, it is notoriously poorly constrained, and it is highly uncertain how it varies, either between or within galaxies. Here we present the results of a proof-of concept study, using the DustPedia data for two nearby face-on spiral galaxies M74 (NGC 628) and M83 (NGC 5236), to create the first ever maps of $κかっぱ_{d}$ in galaxies. We determine $κかっぱ_{d}$ using an empirical method that exploits the fact that the dust-to-metals ratio of the interstellar medium is constrained by direct measurements of the depletion of gas-phase metals. We apply this method pixel-by-pixel within M74 and M83, to create maps of $κかっぱ_{d}$. We also demonstrate a novel method of producing metallicity maps for galaxies with irregularly-sampled measurements, using the machine learning technique of Gaussian process regression. We find strong evidence for significant variation in $κかっぱ_{d}$. We find values of $κかっぱ_{d}$ at 500 $μみゅー$m spanning the range 0.11-0.25 ${\rm m^{2}\,kg^{-1}}$ in M74, and 0.15-0.80 ${\rm m^{2}\,kg^{-1}}$ in M83. Surprisingly, we find that $κかっぱ_{d}$ shows a distinct inverse correlation with the local density of the interstellar medium. This inverse correlation is the opposite of what is predicted by standard dust models. However, we find this relationship to be robust against a large range of changes to our method - only the adoption of unphysical or highly unusual assumptions would be able to suppress it. △ Less

Submitted 22 April, 2022; v1 submitted 12 August, 2019; originally announced August 2019.

Comments: Corrected typographical error in Equation A1, as per erratum accepted by MNRAS on 20th April 2022. No results or conclusions effected by the error, or by this correction

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 489, Issue 4, November 2019, Pages 5256-5283

Probing the cold magnetized Universe with SPICA-POL (B-BOP)

Authors: Ph. André, A. Hughes, V. Guillet, F. Boulanger, A. Bracco, E. Ntormousi, D. Arzoumanian, A. J. Maury, J. -Ph. Bernard, S. Bontemps, I. Ristorcelli, J. M. Girart, F. Motte, K. Tassis, E. Pantin, T. Montmerle, D. Johnstone, S. Gabici, A. Efstathiou, Shantanu Basu, M. Béthermin, H. Beuther, J. Braine, J. Di Francesco, E. Falgarone , et al. (31 additional authors not shown)

Abstract: SPICA, the cryogenic infrared space telescope recently pre-selected for a `Phase A' concept study as one of the three remaining candidates for ESA's fifth medium class (M5) mission, is foreseen to include a far-infrared polarimetric imager (SPICA-POL, now called B-BOP), which would offer a unique opportunity to resolve major issues in our understanding of the nearby, cold magnetized Universe. This… ▽ More SPICA, the cryogenic infrared space telescope recently pre-selected for a `Phase A' concept study as one of the three remaining candidates for ESA's fifth medium class (M5) mission, is foreseen to include a far-infrared polarimetric imager (SPICA-POL, now called B-BOP), which would offer a unique opportunity to resolve major issues in our understanding of the nearby, cold magnetized Universe. This paper presents an overview of the main science drivers for B-BOP, including high dynamic range polarimetric imaging of the cold interstellar medium (ISM) in both our Milky Way and nearby galaxies. Thanks to a cooled telescope, B-BOP will deliver wide-field 100-350 micron images of linearly polarized dust emission in Stokes Q and U with a resolution, signal-to-noise ratio, and both intensity and spatial dynamic ranges comparable to those achieved by Herschel images of the cold ISM in total intensity (Stokes I). The B-BOP 200 micron images will also have a factor ~30 higher resolution than Planck polarization data. This will make B-BOP a unique tool for characterizing the statistical properties of the magnetized interstellar medium and probing the role of magnetic fields in the formation and evolution of the interstellar web of dusty molecular filaments giving birth to most stars in our Galaxy. B-BOP will also be a powerful instrument for studying the magnetism of nearby galaxies and testing galactic dynamo models, constraining the physics of dust grain alignment, informing the problem of the interaction of cosmic rays with molecular clouds, tracing magnetic fields in the inner layers of protoplanetary disks, and monitoring accretion bursts in embedded protostars. △ Less

Submitted 9 May, 2019; originally announced May 2019.

Comments: Accepted for publication in PASA; 37 pages, 16 figures

Journal ref: PASA, 36, e029 (2019)

DustPedia - the relationships between stars, gas and dust for galaxies residing in different environments

Authors: J. I. Davies, A. Nersesian, M. Baes, S. Bianchi, V. Casasola, L. P. Cassara, C. J. R. Clark, I. De Looze, P. De Vis, R. Evans, J. Fritz, M. Galametz, F. Galliano, A. P. Jones, S. Lianou, S. C. Madden, A. V. Mosenkov, M. W. L. Smith, S. Verstocken, S. Viaene, M. Vika, E. Xilouris, N. Ysard

Abstract: We use a sub-set of the DustPedia galaxy sample (461 galaxies) to investigate the effect the environment has had on galaxies. We consider Virgo cluster and field samples and also assign a density contrast parameter to each galaxy, as defined by the local density of SDSS galaxies. We consider their chemical evolution (using M_{Dust}/M_{Baryon} and M_{Gas}/M_{Baryon}), their specific star formation… ▽ More We use a sub-set of the DustPedia galaxy sample (461 galaxies) to investigate the effect the environment has had on galaxies. We consider Virgo cluster and field samples and also assign a density contrast parameter to each galaxy, as defined by the local density of SDSS galaxies. We consider their chemical evolution (using M_{Dust}/M_{Baryon} and M_{Gas}/M_{Baryon}), their specific star formation rate (SFR/M_{Stars}), star formation efficiency (SFR/M_{Gas}), stars-to-dust mass ratio (M_{Stars}/M_{Dust}), gas-to-dust mass ratio (M_{Gas}/M_{Dust}) and the relationship between star formation rate per unit mass of dust and dust temperature (SFR/M_{Dust} and T_{Dust}). Late type galaxies (later than Sc) in all of the environments can be modelled using simple closed box chemical evolution and a simple star formation history (SFR(t) \propto t\exp{-t/τたう}). For earlier type galaxies the physical mechanisms that give rise to their properties are clearly much more varied and require a more complicated model (mergers, gas in or outflow). However, we find little or no difference in the properties of galaxies of the same morphological type within the cluster, field or with different density contrasts. It appears that it is morphology, how and whenever this is laid down, and consistent internal physical processes that primarily determine the derived properties of galaxies in the DustPedia sample and not processes related to differences in the local environment. △ Less

Submitted 30 April, 2019; originally announced April 2019.

Comments: Accepted for publication in Astronomy and Astrophysics

The Life Cycle of Dust

Authors: Sarah Sadavoy, Mikako Matsuura, Lee Armus, Cara Battersby, Caitlin Casey, Christopher Clark, Asantha Cooray, Karine Demyk, Neal Evans, Karl Gordon, Frederic Galliano, Maryvonne Gerin, Benne Holwerda, Nia Imara, Doug Johnstone, Alvaro Labiano, David Leisawitz, Wanggi Lim, Leslie Looney, Margaret Meixner, Eric Murphy, Roberta Paladini, Julia Roman-Duval, Karin Sandstrom, John-David Smith , et al. (2 additional authors not shown)

Abstract: Dust offers a unique probe of the interstellar medium (ISM) across multiple size, density, and temperature scales. Dust is detected in outflows of evolved stars, star-forming molecular clouds, planet-forming disks, and even in galaxies at the dawn of the Universe. These grains also have a profound effect on various astrophysical phenomena from thermal balance and extinction in galaxies to the buil… ▽ More Dust offers a unique probe of the interstellar medium (ISM) across multiple size, density, and temperature scales. Dust is detected in outflows of evolved stars, star-forming molecular clouds, planet-forming disks, and even in galaxies at the dawn of the Universe. These grains also have a profound effect on various astrophysical phenomena from thermal balance and extinction in galaxies to the building blocks for planets, and changes in dust grain properties will affect all of these phenomena. A full understanding of dust in all of its forms and stages requires a multi-disciplinary investigation of the dust life cycle. Such an investigation can be achieved with a statistical study of dust properties across stellar evolution, star and planet formation, and redshift. Current and future instrumentation will enable this investigation through fast and sensitive observations in dust continuum, polarization, and spectroscopy from near-infrared to millimeter wavelengths. △ Less

Submitted 24 April, 2019; originally announced April 2019.

Comments: Science white paper submitted to the Astro2020 Decadal Survey

A systematic metallicity study of DustPedia galaxies reveals evolution in the dust-to-metal ratios

Authors: P. De Vis, A. Jones, S. Viaene, V. Casasola, C. J. R. Clark, M. Baes, S. Bianchi, L. P. Cassara, J. I. Davies, I. De Looze, M. Galametz, F. Galliano, S. Lianou, S. Madden, A. Manilla-Robles, A. V. Mosenkov, A. Nersesian, S. Roychowdhury, E. M. Xilouris, N. Ysard

Abstract: Observations of evolution in the dust-to-metal ratio allow us to constrain the dominant dust processing mechanisms. In this work, we present a study of the dust-to-metal and dust-to-gas ratios in a sub-sample of ~500 DustPedia galaxies. Using literature and MUSE emission line fluxes, we derived gas-phase metallicities (oxygen abundances) for over 10000 individual regions and determine characterist… ▽ More Observations of evolution in the dust-to-metal ratio allow us to constrain the dominant dust processing mechanisms. In this work, we present a study of the dust-to-metal and dust-to-gas ratios in a sub-sample of ~500 DustPedia galaxies. Using literature and MUSE emission line fluxes, we derived gas-phase metallicities (oxygen abundances) for over 10000 individual regions and determine characteristic metallicities for each galaxy. We study how the relative dust, gas, and metal contents of galaxies evolve by using metallicity and gas fraction as proxies for evolutionary state. The global oxygen abundance and nitrogen-to-oxygen ratio are found to increase monotonically as galaxies evolve. Additionally, unevolved galaxies (gas fraction > 60%, metallicity 12 + log(O/H) < 8.2) have dust-to-metal ratios that are about a factor of 2.1 lower (a factor of six lower for galaxies with gas fraction > 80%) than the typical dust-to-metal ratio (Md/MZ ~ 0.214) for more evolved sources. However, for high gas fractions, the scatter is larger due to larger observational uncertainties as well as a potential dependence of the dust grain growth timescale and supernova dust yield on local conditions and star formation histories. We find chemical evolution models with a strong contribution from dust grain growth describe these observations reasonably well. The dust-to-metal ratio is also found to be lower for low stellar masses and high specific star formation rates (with the exception of some sources undergoing a starburst). Finally, the metallicity gradient correlates weakly with the HI-to-stellar mass ratio, the effective radius and the dust-to-stellar mass ratio, but not with stellar mass. △ Less

Submitted 25 October, 2019; v1 submitted 25 January, 2019; originally announced January 2019.

Comments: 24 pages, 15 figures, published in A&A, bibref: 2019A&A...623A...5D, This Version 2 has increased accuracy in Table B1 in appendix B

Journal ref: A&A 623, A5 (2019)

Nano carbon dust emission in proto-planetary disks: the aliphatic-aromatic components

Authors: T. Boutéraon, E. Habart, N. Ysard, A. P. Jones, E. Dartois, T. Pino

Abstract: In the interstellar medium, carbon (nano-)grains are a major component of interstellar dust. This solid phase is more vulnerable to processing and destruction than its silicate counterpart. It exhibits a complex, size-dependent evolution due to interactions within different radiative and dynamical environments. Infrared signatures of these nanocarbon grains are seen in a large number of disks arou… ▽ More In the interstellar medium, carbon (nano-)grains are a major component of interstellar dust. This solid phase is more vulnerable to processing and destruction than its silicate counterpart. It exhibits a complex, size-dependent evolution due to interactions within different radiative and dynamical environments. Infrared signatures of these nanocarbon grains are seen in a large number of disks around Herbig HAeBe stars. We probe the composition and evolution of carbon nano-grains at the surface of (pre-)transitional protoplanetary disks around Herbig stars. We present spatially resolved infrared emission spectra obtained with NAOS CONICA at the VLT in the 3-4 $μみゅー$m range with a spatial resolution of 0.1", which allow us to trace aromatic, olefinic and aliphatic bands which are attributed to sub-nanometer hydrocarbon grains. We apply a gaussian fitting to analyse the observed spectral signatures. Finally, we propose an interpretation in the framework of the The Heterogeneous dust Evolution Model of Interstellar Solids (THEMIS). We show the presence of several spatially extended spectral features, related to aromatic and aliphatic hydrocarbon material in disks around Herbig stars, from ~ 10 to 50-100 au, and even in inner gaps devoided of large grains. The correlation and constant intensity ratios between aliphatic and aromatic CH stretching bands suggest a common nature of the carriers. Given their expected high destruction rates due to UV photons, our observations suggest that they are continuously replenished at the disk surfaces. △ Less

Submitted 22 January, 2019; originally announced January 2019.

Comments: accepted for publication by A&A

Journal ref: A&A 623, A135 (2019)

Dust spectrum and polarisation at 850 um in the massive IRDC G035.39-00.33

Authors: Mika Juvela, Vincent Guillet, Tie Liu, Isabelle Ristorcelli, Veli-Matti Pelkonen, Dana Alina, Leonardo Bronfman, David J. Eden, Kee Tae Kim, Patrick M. Koch, Woojin Kwon, Chang Won Lee, Johanna Malinen, Elisabetta Micelotta, Julien Montillaud, Mark G. Rawlings, Patricio Sanhueza, Archana Soam, Alessio Traficante, Nathalie Ysard, Chuan-Peng Zhang

Abstract: The dust sub-millimetre polarisation of star-forming clouds carries information on dust and the role of magnetic fields in cloud evolution. With observations of a dense filamentary cloud G035.39-00.33, we aim to characterise the dust emission properties and the variations of the polarisation fraction. JCMT SCUBA-2/POL-2 data at 850um are combined with Planck 850um (353GHz) data to map polarisation… ▽ More The dust sub-millimetre polarisation of star-forming clouds carries information on dust and the role of magnetic fields in cloud evolution. With observations of a dense filamentary cloud G035.39-00.33, we aim to characterise the dust emission properties and the variations of the polarisation fraction. JCMT SCUBA-2/POL-2 data at 850um are combined with Planck 850um (353GHz) data to map polarisation fractions. With previous SCUBA-2 observations (450um and 850um) and Herschel data, the column densities are determined via modified blackbody fits and via radiative transfer modelling. Models are constructed to examine how the polarisation angles and fractions depend on potential magnetic field geometries and grain alignment. POL-2 data show clear changes in the magnetic field orientation. The filament has a peak column density of N(H2)~7 10^22 cm-2, a minimum dust temperature of T~12 K, and a mass of some 4300Msun for the area N(H2)> 5 10^21 cm-2. The estimated average value of the dust opacity spectral index is beta ~ 1.9. The ratio of sub-millimetre and J band optical depths is tau(250 um)/tau(J) ~ 2.5 10^-3, more than four times the typical values for diffuse medium. The polarisation fraction decreases as a function of column density to p ~ 1% in the central filament. Because of noise, the observed decrease of p(N) is significant only at N(H2)>2 10^22 cm-2. The observations suggest that the grain alignment is not constant. Although the data can be explained with a complete loss of alignment at densities above ~ 10^4 cm-3 or using the predictions of radiative torques alignment, the uncertainty of the field geometry and the spatial filtering of the SCUBA-2 data prevent strong conclusions. G035.39-00.33 shows strong signs of dust evolution and the low polarisation fraction is suggestive of a loss of polarised emission from its densest parts. △ Less

Submitted 6 September, 2018; v1 submitted 4 September, 2018; originally announced September 2018.

Comments: accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 620, A26 (2018)

The optical properties of dust: the effects of composition, size, and structure

Authors: N. Ysard, A. P. Jones, K. Demyk, T. Boutéraon, M. Koehler

Abstract: Grains are determinant for setting the chemical, physical, dynamical, and radiative properties of all the media in which they are present. Their influence depends on the grain composition, size, and geometrical structure which vary throughout the dust lifecycle. Grain growth arises in dense media as traced by an enhancement of the FIR emissivity and by cloudshine and coreshine. We investigate the… ▽ More Grains are determinant for setting the chemical, physical, dynamical, and radiative properties of all the media in which they are present. Their influence depends on the grain composition, size, and geometrical structure which vary throughout the dust lifecycle. Grain growth arises in dense media as traced by an enhancement of the FIR emissivity and by cloudshine and coreshine. We investigate the imprint of the grain characteristics on the dust unpolarised optical properties from the visible to the FIR for isolated grains and aggregates. Using optical constants for both carbonaceous and silicate materials, we derive the absorption and scattering efficiencies, the asymmetry factor of the phase function, the albedo, and the mass opacity, using either Mie theory or DDA. We study the effects of the size, porosity and shape of the grains, and of the monomers constituting the aggregates, on the optical properties. For aggregates we study the influence of the number of monomers and of mixing monomer sizes. Grain structure changes result in variations at all wavelengths. Porosity, grain elongation, and aggregation all produce an increase in the FIR opacity. The spectral dependence of this increase depends on the nature of the material composing the grain: it is independent of the wavelength for insulators but not for conductors. For aggregates, the FIR increase does not depend on the monomer size and saturates for aggregates containing 6 or more monomers. For silicates, the MIR spectral feature at 18 microns is more sensitive to the details of the grain structure than the 10 microns feature. This study provides a basis for understanding the range of variations achievable as a result of varying the grain characteristics. It emphasises the importance of considering the detailed grain structure in determining the dust optical properties and of using exact methods. △ Less

Submitted 14 June, 2018; originally announced June 2018.

Comments: accepted for publication by A&A

Journal ref: A&A 617, A124 (2018)

French SKA White Book - The French Community towards the Square Kilometre Array

Authors: F. Acero, J. -T. Acquaviva, R. Adam, N. Aghanim, M. Allen, M. Alves, R. Ammanouil, R. Ansari, A. Araudo, E. Armengaud, B. Ascaso, E. Athanassoula, D. Aubert, S. Babak, A. Bacmann, A. Banday, K. Barriere, F. Bellossi, J. -P. Bernard, M. G. Bernardini, M. Béthermin, E. Blanc, L. Blanchet, J. Bobin, S. Boissier , et al. (153 additional authors not shown)

Abstract: The "Square Kilometre Array" (SKA) is a large international radio telescope project characterised, as suggested by its name, by a total collecting area of approximately one square kilometre, and consisting of several interferometric arrays to observe at metric and centimetric wavelengths. The deployment of the SKA will take place in two sites, in South Africa and Australia, and in two successive p… ▽ More The "Square Kilometre Array" (SKA) is a large international radio telescope project characterised, as suggested by its name, by a total collecting area of approximately one square kilometre, and consisting of several interferometric arrays to observe at metric and centimetric wavelengths. The deployment of the SKA will take place in two sites, in South Africa and Australia, and in two successive phases. From its Phase 1, the SKA will be one of the most formidable scientific machines ever deployed by mankind, and by far the most impressive in terms of data throughput and required computing power. With the participation of almost 200 authors from forty research institutes and six private companies, the publication of this French SKA white paper illustrates the strong involvement in the SKA project of the French astronomical community and of a rapidly growing number of major scientific and technological players in the fields of Big Data, high performance computing, energy production and storage, as well as system integration. △ Less

Submitted 28 March, 2018; v1 submitted 19 December, 2017; originally announced December 2017.

Comments: Editor in chief: C. Ferrari; Editors: M. Alves, S. Bosse, S. Corbel, A. Ferrari, K. Ferrière, S. Gauffre, E. Josselin, G. Lagache, S. Lambert, G. Marquette, J.-M. Martin, M.-A. Miville-Deschênes, L. Montier, B. Semelin, G. Theureau, S. Vergani, N. Vilmer, P. Zarka; Original file with high resolution figures at SKA-France link: https://ska-france.oca.eu/images/SKA-France-Media/FWB_051017.pdf

Dust models compatible with Planck intensity and polarization data in translucent lines of sight

Authors: V. Guillet, L. Fanciullo, L. Verstraete, F. Boulanger, A. P. Jones, M. -A. Miville-Deschênes, N. Ysard, F. Levrier, M. Alves

Abstract: The dust properties inferred from the analysis of Planck observations in total and polarized emission challenge current dust models. We propose new dust models compatible with polarized and unpolarized data in extinction and emission for translucent lines of sight ($0.5 < A_V < 2.5$). We amended the DustEM tool to model polarized extinction and emission. We fit the spectral dependence of the mean… ▽ More The dust properties inferred from the analysis of Planck observations in total and polarized emission challenge current dust models. We propose new dust models compatible with polarized and unpolarized data in extinction and emission for translucent lines of sight ($0.5 < A_V < 2.5$). We amended the DustEM tool to model polarized extinction and emission. We fit the spectral dependence of the mean extinction, polarized extinction, SED and polarized SED with PAHs, astrosilicates and amorphous carbon (a-C). The astrosilicate population is aligned along the magnetic field lines, while the a-C population may be aligned or not. With their current optical properties, oblate astrosilicate grains are not emissive enough to reproduce the emission to extinction polarization ratio $P_{353}/p_V$ derived with Planck data. Models using prolate astrosilicate grains with an elongation $a/b=3$ and an inclusion of 20% of porosity succeed. The spectral dependence of the polarized SED is steeper in our models than in the data. Models perform slightly better when a-C grains are aligned. A small (6%) volume inclusion of a-C in the astrosilicate matrix removes the need for porosity and perfect grain alignment, and improves the fit to the polarized SED. Dust models based on astrosilicates can be reconciled with Planck data by adapting the shape of grains and adding inclusions of porosity or a-C in the astrosilicate matrix. △ Less

Submitted 12 October, 2017; originally announced October 2017.

Journal ref: A&A 610, A16 (2018)

Radial distribution of dust, stars, gas, and star-formation rate in DustPedia face-on galaxies

Authors: V. Casasola, L. P. Cassara, S. Bianchi, S. Verstocken, E. Xilouris, L. Magrini, M. W. L. Smith, I. De Looze, M. Galametz, S. C. Madden, M. Baes, C. Clark, J. Davies, P. De Vis, R. Evans, J. Fritz, F. Galliano, A. P. Jones, A. V. Mosenkov, S. Viaene, N. Ysard

Abstract: The purpose of this work is the characterization of the radial distribution of dust, stars, gas, and star-formation rate (SFR) in a sub-sample of 18 face-on spiral galaxies extracted from the DustPedia sample. This study is performed by exploiting the multi-wavelength, from UV to sub-mm bands, DustPedia database, in addition to molecular (12CO) and atomic (HI) gas maps and metallicity abundance in… ▽ More The purpose of this work is the characterization of the radial distribution of dust, stars, gas, and star-formation rate (SFR) in a sub-sample of 18 face-on spiral galaxies extracted from the DustPedia sample. This study is performed by exploiting the multi-wavelength, from UV to sub-mm bands, DustPedia database, in addition to molecular (12CO) and atomic (HI) gas maps and metallicity abundance information available in the literature. We fitted the surface brightness profiles of the tracers of dust and stars, the mass surface density profiles of dust, stars, molecular gas, and total gas, and the SFR surface density profiles with an exponential curve and derived their scale-lengths. We also developed a method to solve for the CO-to-H2 conversion factor (alpha_CO) per galaxy by using dust and gas mass profiles. Although each galaxy has its own peculiar behaviour, we identified a common trend of the exponential scale-lengths vs. wavelength. On average, the scale-lengths normalized to the B-band 25 mag/arcsec^2 radius decrease from UV to 70 micron, from 0.4 to 0.2, and then increase back up to 0.3 at 500 microns. The main result is that, on average, the dust mass surface density scale-length is about 1.8 times the stellar one derived from IRAC data and the 3.6 micron surface brightness, and close to that in the UV. We found a mild dependence of the scale-lengths on the Hubble stage T: the scale-lengths of the Herschel bands and the 3.6 micron scale-length tend to increase from earlier to later types, the scale-length at 70 micron tends to be smaller than that at longer sub-mm wavelength with ratios between longer sub-mm wavelengths and 70 micron that decrease with increasing T. The scale-length ratio of SFR and stars shows a weak increasing trend towards later types. △ Less

Submitted 16 June, 2017; originally announced June 2017.

Comments: 33 pages, 9 figures (at low resolution here), 8 tables, Accepted for publication in A&A

Journal ref: A&A 605, A18 (2017)

The global dust modelling framework THEMIS (The Heterogeneous dust Evolution Model for Interstellar Solids)

Authors: A. P. Jones, M. Koehler, N. Ysard, M. Bocchio, L. Verstraete

Abstract: Here we introduce the interstellar dust modelling framework THEMIS (The Heterogeneous dust Evolution Model for Interstellar Solids), which takes a global view of dust and its evolution in response to the local conditions in interstellar media. This approach is built upon a core model that was developed to explain the dust extinction and emission in the diffuse interstellar medium. The model was th… ▽ More Here we introduce the interstellar dust modelling framework THEMIS (The Heterogeneous dust Evolution Model for Interstellar Solids), which takes a global view of dust and its evolution in response to the local conditions in interstellar media. This approach is built upon a core model that was developed to explain the dust extinction and emission in the diffuse interstellar medium. The model was then further developed to self-consistently include the effects of dust evolution in the transition to denser regions. The THEMIS approach is under continuous development and currently we are extending the framework to explore the implications of dust evolution in HII regions and the photon-dominated regions associated with star formation. We provide links to the THEMIS, DustEM and DustPedia websites where more information about the model, its input data and applications can be found. △ Less

Submitted 2 March, 2017; originally announced March 2017.

Comments: 9 pages, 1 figure

Journal ref: A&A 602, A46 (2017)

Galactic Cold Cores. VIII. Filament formation and evolution: Filament properties in context with evolutionary models

Authors: A. Rivera-Ingraham, I. Ristorcelli, M. Juvela, J. Montillaud, A. Men'shchikov, J. Malinen, V. -M. Pelkonen, A. Marston, P. G. Martin, L. Pagani, R. Paladini, D. Paradis, N. Ysard, D. Ward-Thompson, J. -P. Bernard, D. J. Marshall, L. Montier, L. V. Toth

Abstract: Filaments are key for star formation models. As part of the study carried out by the Herschel GCC Programme, here we study the filament properties presented in GCC.VII in context with theoretical models of filament formation and evolution. A conservative sample of filaments at a distance D<500pc was extracted with the Getfilaments algorithm. Their physical structure was quantified according to two… ▽ More Filaments are key for star formation models. As part of the study carried out by the Herschel GCC Programme, here we study the filament properties presented in GCC.VII in context with theoretical models of filament formation and evolution. A conservative sample of filaments at a distance D<500pc was extracted with the Getfilaments algorithm. Their physical structure was quantified according to two main components: the central (Gaussian) region (core component), and the power-law like region dominating the filament column density profile at larger radii (wing component). The properties and behaviour of these components relative to the total linear mass density of the filament and its environmental column density were compared with theoretical models describing the evolution of filaments under gravity-dominated conditions. The feasibility of a transition to supercritical state by accretion is dependent on the combined effect of filament intrinsic properties and environmental conditions. Reasonably self-gravitating (high Mline-core) filaments in dense environments (av\sim3mag) can become supercritical in timescales of t\sim1Myr by accreting mass at constant or decreasing width. The trend of increasing Mline-tot (Mline-core and Mline-wing), and ridge Av with background also indicates that the precursors of star-forming filaments evolve coevally with their environment. The simultaneous increase of environment and filament Av explains the association between dense environments and high Mline-core values, and argues against filaments remaining in constant single-pressure equilibrium states. The simultaneous growth of filament and background in locations with efficient mass assembly, predicted in numerical models of collapsing clouds, presents a suitable scenario for the fulfillment of the combined filament mass-environment criterium that is in quantitative agreement with Herschel observations. △ Less

Submitted 24 February, 2017; originally announced February 2017.

Comments: Accepted in A&A

Journal ref: A&A 601, A94 (2017)

Modeling dust emission in the Magellanic Clouds with Spitzer and Herschel

Authors: Jeremy Chastenet, Caroline Bot, Karl D. Gordon, Marco Bocchio, Julian Roman-Duval, Anthony P. Jones, Nathalie Ysard

Abstract: Dust modeling is crucial to infer dust properties and budget for galaxy studies. However, there are systematic disparities between dust grain models that result in corresponding systematic differences in the inferred dust properties of galaxies. Quantifying these systematics requires a consistent fitting analysis. We compare the output dust parameters and assess the differences between two dust gr… ▽ More Dust modeling is crucial to infer dust properties and budget for galaxy studies. However, there are systematic disparities between dust grain models that result in corresponding systematic differences in the inferred dust properties of galaxies. Quantifying these systematics requires a consistent fitting analysis. We compare the output dust parameters and assess the differences between two dust grain models, Compiegne et al (2011), and THEMIS (Jones et al 2013, Kohler et al 2015). In this study, we use a single fitting method applied to all the models to extract a coherent and unique statistical analysis. We fit the models to the dust emission seen by Spitzer and Herschel in the Small and Large Magellanic Clouds (SMC and LMC). The observations cover the infrared (IR) spectrum from a few microns to the sub-millimeter range. For each fitted pixel, we calculate the full n-D likelihood, based on the method described in Gordon et al (2014). The free parameters are both environmental ($U$, the interstellar radiation field strength; $αあるふぁ_\mathrm{ISRF}$, power-law coefficient for a multi-U environment; $Ωおめが^*$, the starlight strength) and intrinsic to the model ($Y_\mathrm{i}$: abundances of the grain species $i$; $αあるふぁ_\mathrm{sCM20}$, coefficient in the small carbon grain size distribution). Fractional residuals of 5 different sets of parameters show that fitting THEMIS brings a more accurate reproduction of the observations than the Compiegne model. However, independent variations of the dust species show strong model-dependencies. We find that the abundance of silicates can only be constrained to an upper-limit and the silicate/carbon ratio is different than that seen in our Galaxy. In the LMC, our fits result in dust masses slightly lower than those found in literature , by a factor lower than 2. In the SMC, we find dust masses in agreement with previous studies. △ Less

Submitted 22 December, 2016; originally announced December 2016.

Comments: 17 pages, 17 figures; accepted for publication in A&A

Journal ref: A&A 601, A55 (2017)

DustPedia - A Definitive Study of Cosmic Dust in the Local Universe

Authors: J. I. Davies, M. Baes, S. Bianchi, A. Jones, S. Madden, M. Xilouris, M. Bocchio, V. Casasola, L. Cassara, C. Clark, I. De Looze, R. Evans, J. Fritz, M. Galametz, F. Galliano, S. Lianou, A. V. Mosenkov, M. Smith, S. Verstocken, S. Viaene, M. Vika, G. Wagle, N. Ysard

Abstract: The European Space Agency has invested heavily in two cornerstones missions; Herschel and Planck. The legacy data from these missions provides us with an unprecedented opportunity to study cosmic dust in galaxies so that we can answer fundamental questions about, for example: the origin of the chemical elements, physical processes in the interstellar medium (ISM), its effect on stellar radiation,… ▽ More The European Space Agency has invested heavily in two cornerstones missions; Herschel and Planck. The legacy data from these missions provides us with an unprecedented opportunity to study cosmic dust in galaxies so that we can answer fundamental questions about, for example: the origin of the chemical elements, physical processes in the interstellar medium (ISM), its effect on stellar radiation, its relation to star formation and how this relates to the cosmic far infrared background. In this paper we describe the DustPedia project, which is enabling us to develop tools and computer models that will help us relate observed cosmic dust emission to its physical properties (chemical composition, size distribution, temperature), to its origins (evolved stars, super novae, growth in the ISM) and the processes that destroy it (high energy collisions and shock heated gas). To carry out this research we will combine the Herschel/Planck data with that from other sources of data, providing observations at numerous wavelengths (< 41) across the spectral energy distribution, thus creating the DustPedia database. To maximise our spatial resolution and sensitivity to cosmic dust we limit our analysis to 4231 local galaxies (v < 3000 km/s) selected via their near infrared luminosity (stellar mass). To help us interpret the data we have developed a new physical model for dust (THEMIS), a new Bayesian method of fitting and interpreting spectral energy distributions (HerBIE) and a state-of-the-art Monte Carlo photon tracing radiative transfer model (SKIRT). In this the first of the DustPedia papers we describe the project objectives, data sets used and provide an insight into the new scientific methods we plan to implement. △ Less

Submitted 20 September, 2016; originally announced September 2016.

Comments: 66 pages, accepted for publication in PASP

The Herschel Virgo Cluster Survey. XX. Dust and gas in the foreground Galactic cirrus

Authors: S. Bianchi, C. Giovanardi, M. W. L. Smith, J. Fritz, J. I. Davies, M. P. Haynes, R. Giovanelli, M. Baes, M. Bocchio, S. Boissier, M. Boquien, A. Boselli, V. Casasola, C. J. R. Clark, I. De Looze, S. di Serego Alighieri, M. Grossi, A. P. Jones, T. M. Hughes, L. K. Hunt, S. Madden, L. Magrini, C. Pappalardo, N. Ysard, S. Zibetti

Abstract: We study the correlation between far-infared/submm dust emission and atomic gas column density in order to derive the properties of the high Galactic latitude, low density, Milky Way cirrus in the foreground of the Virgo cluster of galaxies. Dust emission maps from 60 to 850 um are obtained from SPIRE observations carried out within the Herschel Virgo Cluster Survey, complemented by IRAS-IRIS and… ▽ More We study the correlation between far-infared/submm dust emission and atomic gas column density in order to derive the properties of the high Galactic latitude, low density, Milky Way cirrus in the foreground of the Virgo cluster of galaxies. Dust emission maps from 60 to 850 um are obtained from SPIRE observations carried out within the Herschel Virgo Cluster Survey, complemented by IRAS-IRIS and Planck-HFI maps. Data from the Arecibo legacy Fast ALFA Survey is used to derive atomic gas column densities for two broad velocity components, low and intermediate velocity clouds. Dust emissivities are derived for each gas component and each far-infared/submm band. For the low velocity clouds, we measure an average emissivity 0.79 +/- 0.08 times 1E-20 MJy sr^-1 cm^2 at 250 um. After fitting a modified blackbody to the available bands, we estimated a dust absorption cross-section 0.49 +/- 0.13 times 1E-25 cm^2 H^-1 at 250 um (with dust temperature T = 20.4 +/- 1.5 K and spectral index beta = 1.53 +/- 0.17). The results are in excellent agreement with those obtained by Planck over a much larger coverage of the high Galactic latitude cirrus (50% of the sky vs 0.2% in our work). For dust associated with intermediate velocity gas, we confirm earlier Planck results and find a higher temperature and lower emissivity and cross-section. After subtracting the modelled components, we find regions at scales smaller than 20' where the residuals deviate significantly from the average, cosmic-infrared-background dominated, scatter. These large residuals are most likely due to local variations in the cirrus dust properties (and/or the dust/atomic-gas correlation) or to high-latitude molecular clouds with average N_H2 <~ 1E20 cm^-2. We find no conclusive evidence for intracluster dust emission in Virgo. △ Less

Submitted 19 September, 2016; originally announced September 2016.

Comments: 25 pages, 19 figures, A&A accepted

Journal ref: A&A 597, A130 (2017)