-
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.
-
Measuring Dust Attenuation Curves of SINGS/KINGFISH Galaxies Using Swift/UVOT Photometry
Authors:
Alexander Belles,
Marjorie Decleir,
William P. Bowman,
Lea M. Z. Hagen,
Caryl Gronwall,
Michael H. Siegel
Abstract:
We present Swift/Ultraviolet Optical Telescope (UVOT) integrated light photometry of the Spitzer Infrared Nearby Galaxies Survey (SINGS) and the Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel (KINGFISH) samples of nearby galaxies. Combining the Swift/UVOT data with archival photometry, we investigate a variety of dust attenuation curves derived using MCSED, a flexible spectra…
▽ More
We present Swift/Ultraviolet Optical Telescope (UVOT) integrated light photometry of the Spitzer Infrared Nearby Galaxies Survey (SINGS) and the Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel (KINGFISH) samples of nearby galaxies. Combining the Swift/UVOT data with archival photometry, we investigate a variety of dust attenuation curves derived using MCSED, a flexible spectral energy distribution fitting code. We fit the panchromatic data using three different star formation history (SFH) parameterizations: a decaying exponential, a double power law, and a piecewise function with breaks at physically motivated ages. We find that the average attenuation law of the sample changes slightly based on the SFH assumed. Specifically, the exponential SFH leads to the shallowest attenuation curves. Using simulated data, we also find the exponential SFH fails to outperform the more complex SFHs. Finally, we find a systematic offset in the derived bump strength between SED fits with and without UVOT data, where the inclusion of UVOT data leads to smaller bump strengths, highlighting the importance of the UVOT data. This discrepancy is not seen in fits to mock photometry. Understanding dust attenuation in the local universe is key to understanding high redshift objects where rest-frame far-infrared data is unavailable.
△ Less
Submitted 9 May, 2023;
originally announced May 2023.
-
One Relation for All Wavelengths: The Far-Ultraviolet to Mid-Infrared Milky Way Spectroscopic R(V) Dependent Dust Extinction Relationship
Authors:
Karl D. Gordon,
Geoffrey C. Clayton,
Marjorie Decleir,
E. L. Fitzpatrick,
Derck Massa,
Karl A. Misselt,
Erik J. Tollerud
Abstract:
Dust extinction is one of the fundamental measurements of dust grain sizes, compositions, and shapes. Most of the wavelength dependent variations seen in Milky Way extinction are strongly correlated with the single parameter R(V)=A(V)/E(B-V). Existing R(V) dependent extinction relationships use a mixture of spectroscopic and photometry observations, hence do not fully capture all the important dus…
▽ More
Dust extinction is one of the fundamental measurements of dust grain sizes, compositions, and shapes. Most of the wavelength dependent variations seen in Milky Way extinction are strongly correlated with the single parameter R(V)=A(V)/E(B-V). Existing R(V) dependent extinction relationships use a mixture of spectroscopic and photometry observations, hence do not fully capture all the important dust features nor continuum variations. Using four existing samples of spectroscopically measured dust extinction curves, we consistently measure the R(V) dependent extinction relationship spectroscopically from the far-ultraviolet to mid-infrared for the first time. Linear fits of A(lambda)/A(V) dependent on R(V) are done using a method that fully accounts for their significant and correlated uncertainties. These linear parameters are fit with analytic wavelength dependent functions to determine the smooth R(V) (2.3-5.6) and wavelength (912 A-32 micron) dependent extinction relationship. This relationship shows that the far-UV rise, 2175 A bump, and the three broad optical features are dependent on R(V), but the 10 and 20 micron features are not. Existing literature relationships show significant deviations compared to this relationship especially in the far-ultraviolet and infrared. Extinction curves that clearly deviate from this relationship illustrate that this relationship only describes the average behavior versus R(V). We find tentative evidence that the relationship may not be linear with 1/R(V) especially in the ultraviolet. For the first time, this relationship provides measurements of dust extinction that spectroscopically resolve the continuum and features in the ultraviolet, optical, and infrared as a function of R(V) enabling detailed studies of dust grains properties and full spectroscopic accounting for the effects of dust extinction on astrophysical objects.
△ Less
Submitted 6 April, 2023; v1 submitted 4 April, 2023;
originally announced April 2023.
-
The recent star formation history of NGC 628 on resolved scales
Authors:
Maria Lomaeva,
Ilse De Looze,
Amélie Saintonge,
Marjorie Decleir
Abstract:
Star formation histories (SFHs) are integral to our understanding of galaxy evolution. We can study recent SFHs by comparing the star formation rate (SFR) calculated using different tracers, as each probes a different timescale. We aim to calibrate a proxy for the present-day rate of change in SFR, dSFR/dt, which does not require full spectral energy distribution (SED) modeling and depends on as f…
▽ More
Star formation histories (SFHs) are integral to our understanding of galaxy evolution. We can study recent SFHs by comparing the star formation rate (SFR) calculated using different tracers, as each probes a different timescale. We aim to calibrate a proxy for the present-day rate of change in SFR, dSFR/dt, which does not require full spectral energy distribution (SED) modeling and depends on as few observables as possible, to guarantee its broad applicability. To achieve this, we create a set of models in CIGALE and define an SFR change diagnostic as the ratio of the SFR averaged over the past 5 and 200 Myr, <SFR5>/<SFR200>, probed by the H$α$-FUV colour. We apply <SFR5>/<SFR200> to the nearby spiral NGC 628 and find that its star formation activity has overall been declining in the recent past, with the spiral arms, however, maintaining a higher level of activity. The impact of the spiral arm structure is observed to be stronger on <SFR5>/<SFR200> than on the star formation efficiency (SFE$_\text{H$_2$}$). In addition, increasing disk pressure tends to increase recent star formation, and consequently <SFR5>/<SFR200>. We conclude that <SFR5>/<SFR200> is sensitive to the molecular gas content, spiral arm structure, and disk pressure. The <SFR5>/<SFR200> indicator is general and can be used to reconstruct the recent SFH of any star-forming galaxy for which H$α$, FUV, and either mid- or far-IR photometry is available, without the need of detailed modeling.
△ Less
Submitted 13 October, 2022;
originally announced October 2022.
-
SpeX near-infrared spectroscopic extinction curves in the Milky Way
Authors:
Marjorie Decleir,
Karl D. Gordon,
Jennifer E. Andrews,
Geoffrey C. Clayton,
Michael C. Cushing,
Karl A. Misselt,
Yvonne Pendleton,
John Rayner,
William D. Vacca,
D. C. B. Whittet
Abstract:
Interstellar dust extinction curves provide valuable information about dust properties, including the composition and size of the dust grains, and are essential to correct observations for the effects of interstellar dust. In this work, we measure a representative sample of near-infrared (NIR; 0.8-5.5 $μ$m) spectroscopic extinction curves for the first time, enabling us to investigate the extincti…
▽ More
Interstellar dust extinction curves provide valuable information about dust properties, including the composition and size of the dust grains, and are essential to correct observations for the effects of interstellar dust. In this work, we measure a representative sample of near-infrared (NIR; 0.8-5.5 $μ$m) spectroscopic extinction curves for the first time, enabling us to investigate the extinction at wavelengths where it is usually only measured in broad photometric bands. We use IRTF/SpeX spectra of a sample of reddened and comparison stars to measure 15 extinction curves with the pair method. Our sample spans A(V) values from 0.78 to 5.65 and R(V) values from 2.43 to 5.33. We confirm that the NIR extinction curves are well fit by a power law, with indices and amplitudes differing from sight line to sight line. Our average diffuse NIR extinction curve can be represented by a single power law with index $α= 1.7$, but because of the sight line-to-sight line variations, the shape of any average curve will depend on the parental sample. We find that most of the variation in our sample can be linked to the ratio of total-to-selective extinction R(V), a rough measurement of the average dust grain size. Two sight lines in our sample clearly show the ice extinction feature at 3 $μ$m, which can be fitted by a modified Drude profile. We find tentative ice detections with slightly over 3$σ$ significance in two other sight lines. In our average diffuse extinction curve, we measure a 3$σ$ upper limit of A(ice)/A(V) = 0.0021 for this ice feature.
△ Less
Submitted 10 May, 2022; v1 submitted 28 April, 2022;
originally announced April 2022.
-
Milky Way Mid-Infrared Spitzer Spectroscopic Extinction Curves: Continuum and Silicate Features
Authors:
Karl D. Gordon,
Karl A. Misselt,
Jeroen Bouwman,
Geoffrey C. Clayton,
Marjorie Decleir,
Dean C. Hines,
Yvonne Pendleton,
George Rieke,
J. D. T. Smith,
D. C. B. Whittet
Abstract:
We measured the mid-infrared (MIR) extinction using Spitzer photometry and spectroscopy (3.6--37 micron) for a sample of Milky Way sightlines (mostly) having measured ultraviolet extinction curves. We used the pair method to determine the MIR extinction that we then fit with a power law for the continuum and modified Drude profiles for the silicate features. We derived 16 extinction curves having…
▽ More
We measured the mid-infrared (MIR) extinction using Spitzer photometry and spectroscopy (3.6--37 micron) for a sample of Milky Way sightlines (mostly) having measured ultraviolet extinction curves. We used the pair method to determine the MIR extinction that we then fit with a power law for the continuum and modified Drude profiles for the silicate features. We derived 16 extinction curves having a range of A(V) (1.8-5.5) and R(V) values (2.4-4.3). Our sample includes two dense sightlines that have 3 micron ice feature detections and weak 2175 A bumps. The average A(lambda)/A(V) diffuse sightline extinction curve we calculate is lower than most previous literature measurements. This agrees better with literature diffuse dust grain models, though it is somewhat higher. The 10 micron silicate feature does not correlate with the 2175 A bump, for the first time providing direct observational confirmation that these two features arise from different grain populations. The strength of the 10 micron silicate feature varies by $\sim$2.5 and is not correlated with A(V) or R(V). It is well fit by a modified Drude profile with strong correlations seen between the central wavelength, width, and asymmetry. We do not detect other features with limits in A(lambda)/A(V) units of 0.0026 (5--10 micron), 0.004 (10--20 micron), and 0.008 (20-40 micron). We find that the standard prescription of estimating R(V) from C times E(K_s-V)/E(B-V) has C = -1.14 and a scatter of $\sim$7%. Using the IRAC 5.6 micron band instead of K_s gives C = -1.03 and the least scatter of $\sim$3\%.
△ Less
Submitted 12 May, 2021; v1 submitted 11 May, 2021;
originally announced May 2021.
-
JINGLE -- IV. Dust, HI gas and metal scaling laws in the local Universe
Authors:
I. De Looze,
I. Lamperti,
A. Saintonge,
M. Relano,
M. W. L. Smith,
C. J. R. Clark,
C. D. Wilson,
M. Decleir,
A. P. Jones,
R. C. Kennicutt,
G. Accurso,
E. Brinks,
M. Bureau,
P. Cigan,
D. L. Clements,
P. De Vis,
L Fanciullo,
Y. Gao,
W. K. Gear,
L. C. Ho,
H. S. Hwang,
M. J. Michalowski,
J. C. Lee,
C. Li,
L. Lin
, et al. (7 additional authors not shown)
Abstract:
Scaling laws of dust, HI gas and metal mass with stellar mass, specific star formation rate and metallicity are crucial to our understanding of the buildup of galaxies through their enrichment with metals and dust. In this work, we analyse how the dust and metal content varies with specific gas mass ($M_{\text{HI}}$/$M_{\star}$) across a diverse sample of 423 nearby galaxies. The observed trends a…
▽ More
Scaling laws of dust, HI gas and metal mass with stellar mass, specific star formation rate and metallicity are crucial to our understanding of the buildup of galaxies through their enrichment with metals and dust. In this work, we analyse how the dust and metal content varies with specific gas mass ($M_{\text{HI}}$/$M_{\star}$) across a diverse sample of 423 nearby galaxies. The observed trends are interpreted with a set of Dust and Element evolUtion modelS (DEUS) - incluidng stellar dust production, grain growth, and dust destruction - within a Bayesian framework to enable a rigorous search of the multi-dimensional parameter space. We find that these scaling laws for galaxies with $-1.0\lesssim \log M_{\text{HI}}$/$M_{\star}\lesssim0$ can be reproduced using closed-box models with high fractions (37-89$\%$) of supernova dust surviving a reverse shock, relatively low grain growth efficiencies ($ε$=30-40), and long dus lifetimes (1-2\,Gyr). The models have present-day dust masses with similar contributions from stellar sources (50-80\,$\%$) and grain growth (20-50\,$\%$). Over the entire lifetime of these galaxies, the contribution from stardust ($>$90\,$\%$) outweighs the fraction of dust grown in the interstellar medium ($<$10$\%$). Our results provide an alternative for the chemical evolution models that require extremely low supernova dust production efficiencies and short grain growth timescales to reproduce local scaling laws, and could help solving the conundrum on whether or not grains can grow efficiently in the interstellar medium.
△ Less
Submitted 2 June, 2020;
originally announced June 2020.
-
DustKING - the story continues: dust attenuation in NGC628
Authors:
Marjorie Decleir,
Ilse De Looze,
Médéric Boquien,
Maarten Baes
Abstract:
Dust attenuation is a crucial but highly uncertain parameter that hampers the determination of intrinsic galaxy properties, such as stellar masses, star formation rates and star formation histories. The shape of the dust attenuation law is not expected to be uniform between galaxies, nor within a galaxy. Our DustKING project was introduced at the first BINA workshop in 2016 and aims to study the v…
▽ More
Dust attenuation is a crucial but highly uncertain parameter that hampers the determination of intrinsic galaxy properties, such as stellar masses, star formation rates and star formation histories. The shape of the dust attenuation law is not expected to be uniform between galaxies, nor within a galaxy. Our DustKING project was introduced at the first BINA workshop in 2016 and aims to study the variations of dust attenuation curves in nearby galaxies. At the second BINA workshop in 2018, I presented the results of our pilot study for the spiral galaxy NGC628. We find that the average attenuation law of this galaxy is characterised by a MW-like bump and a steep UV slope. Furthermore, we observe intriguing variations within the galaxy, with regions of high $A_V$ exhibiting a shallower attenuation curve. Finally, we discuss how our work might benefit from data taken with the UVIT from the Indian AstroSat mission.
△ Less
Submitted 12 April, 2019;
originally announced April 2019.
-
Revealing the dust attenuation properties on resolved scales in NGC628 with SWIFT UVOT data
Authors:
Marjorie Decleir,
Ilse De Looze,
Médéric Boquien,
Maarten Baes,
Sam Verstocken,
Daniela Calzetti,
Laure Ciesla,
Jacopo Fritz,
Rob Kennicutt,
Angelos Nersesian,
Mathew Page
Abstract:
Understanding how dust attenuation laws vary between and within galaxies is a key question if we want to reliably measure the physical properties of galaxies at both global and local scales. To shed new light on this question, we present a detailed study of the slope and bump strength of the attenuation law in the nearby spiral galaxy NGC628 at the resolved spatial scale of 325 pc. To do so, we ha…
▽ More
Understanding how dust attenuation laws vary between and within galaxies is a key question if we want to reliably measure the physical properties of galaxies at both global and local scales. To shed new light on this question, we present a detailed study of the slope and bump strength of the attenuation law in the nearby spiral galaxy NGC628 at the resolved spatial scale of 325 pc. To do so, we have modelled a broad multi-wavelength dataset from the ultraviolet (UV) to the infrared (IR) with the state-of-the-art SED fitting code CIGALE, including SWIFT UVOT data for which we have developed a new optimized reduction pipeline. We find that the median dust attenuation curve of NGC628 is fairly steep, but not as steep as the SMC curve, and has a sub-MW-type UV bump. We observe intriguing variations within the galaxy, with regions of high $A_V$ exhibiting a shallower attenuation curve. We argue that the flattening of the curve is due to a dominance of absorption-over-scattering events at higher $A_V$. No trend between the bump strength and the IRAC 8.0 micron emission was found. However, this does not necessarily rule out PAHs as the main contributors to the UV bump.
△ Less
Submitted 15 March, 2019;
originally announced March 2019.
-
Old and young stellar populations in DustPedia galaxies and their role in dust heating
Authors:
A. Nersesian,
E. M. Xilouris,
S. Bianchi,
F. Galliano,
A. P. Jones,
M. Baes,
V. Casasola,
L. P. Cassara,
C. J. R. Clark,
J. I. Davies,
M. Decleir,
W. Dobbels,
I. De Looze,
P. De Vis,
J. Fritz,
M. Galametz,
S. C. Madden,
A. V. Mosenkov,
A. Trcka,
S. Verstocken,
S. Viaene,
S. Lianou
Abstract:
Within the framework of the DustPedia project we investigate the properties of cosmic dust and its interaction with the stellar radiation (originating from different stellar populations) for 814 galaxies in the nearby Universe, all observed by the Herschel Space Observatory. We take advantage of the widely used galaxy SED fitting code CIGALE, properly adapted to include the state-of-the-art dust m…
▽ More
Within the framework of the DustPedia project we investigate the properties of cosmic dust and its interaction with the stellar radiation (originating from different stellar populations) for 814 galaxies in the nearby Universe, all observed by the Herschel Space Observatory. We take advantage of the widely used galaxy SED fitting code CIGALE, properly adapted to include the state-of-the-art dust model THEMIS. Using the DustPedia photometry we determine the physical properties of the galaxies, such as, the dust and stellar mass, the star-formation rate, the bolometric luminosity as well as the unattenuated and the absorbed by dust stellar light, for both the old (> 200 Myr) and young (<= 200 Myr) stellar populations. We show how the mass of stars, dust, and atomic gas, as well as the star-formation rate and the dust temperature vary between galaxies of different morphologies and provide recipes to estimate these parameters given their Hubble stage (T). We find a mild correlation between the mass fraction of the small a-C(:H) grains with the specific star-formation rate. On average, young stars are very efficient in heating the dust, with absorption fractions reaching as high as ~77% of the total, unattenuated luminosity of this population. On the other hand, the maximum absorption fraction of old stars is ~24%. Dust heating in early-type galaxies is mainly due to old stars, up to a level of ~90%. Young stars progressively contribute more for `typical' spiral galaxies and they become the dominant source of dust heating for Sm type and irregular galaxies, donating up to ~60% of their luminosity to this purpose. Finally, we find a strong correlation of the dust heating fraction by young stars with morphology and the specific star-formation rate.
△ Less
Submitted 16 April, 2019; v1 submitted 14 March, 2019;
originally announced March 2019.
-
HERschel Observations of Edge-on Spirals (HEROES). IV. Dust energy balance problem
Authors:
Aleksandr V. Mosenkov,
Flor Allaert,
Maarten Baes,
Simone Bianchi,
Peter Camps,
Christopher J. R. Clark,
Marjorie Decleir,
Gert De Geyter,
Ilse De Looze,
Jacopo Fritz,
Gianfranco Gentile,
Benne W. Holwerda,
Thomas M. Hughes,
Fraser Lewis,
Matthew W. L. Smith,
Joris Verstappen,
Sam Verstocken,
Sébastien Viaene
Abstract:
We present results of the detailed dust energy balance study for the seven large edge-on galaxies in the HEROES sample using 3D radiative transfer (RT) modelling. Based on available optical and near-infrared observations of the HEROES galaxies, we derive the 3D distribution of stars and dust in these galaxies. For the sake of uniformity, we apply the same technique to retrieve galaxy properties fo…
▽ More
We present results of the detailed dust energy balance study for the seven large edge-on galaxies in the HEROES sample using 3D radiative transfer (RT) modelling. Based on available optical and near-infrared observations of the HEROES galaxies, we derive the 3D distribution of stars and dust in these galaxies. For the sake of uniformity, we apply the same technique to retrieve galaxy properties for the entire sample: we use a stellar model consisting of a Sérsic bulge and three double-exponential discs (a superthin disc for a young stellar population and thin and thick discs for old populations). For the dust component, we adopt a double-exponential disc with the new THEMIS dust-grain model. We fit oligochromatic radiative transfer (RT) models to the optical and near-infrared images with the fitting algorithm FitSKIRT and do panchromatic simulations with the SKIRT code at wavelengths ranging from ultraviolet to submillimeter. We confirm the previously stated dust energy balance problem in galaxies: for the HEROES galaxies, the dust emission derived from our RT calculations underestimates the real observations by a factor 1.5-4 for all galaxies except NGC 973 and NGC 5907 (apparently, the latter galaxy has a more complex geometry than we used). The comparison between our RT simulations and the observations at mid-infrared-submillimeter wavelengths shows that most of our galaxies exhibit complex dust morphologies (possible spiral arms, star-forming regions, more extended dust structure in the radial and vertical directions). We suggest that, in agreement with the results from Saftly et al. (2015), the large- and small-scale structure is the most probable explanation for the dust energy balance problem.
△ Less
Submitted 24 April, 2018;
originally announced April 2018.