Far-Ultraviolet to Mid-Infrared Dust Extinction Curves in the Milky Way, LMC, and SMC

Slide 1

Slide 1 text

Karl D. Gordon, Astronomer STScI, Baltimore, MD, USA Visiting Prof, Ghent University Ghent University 29 Apr 2024 [email protected] karllark@github “Have Dust – Will Study” Slides on speakerdeck Far-Ultraviolet to Mid-Infrared Dust Extinction Curves in the Milky Way, LMC, and SMC

Slide 2

Slide 2 text

Extinction ● Critical clues to dust grain abundance, size, composition, and shape ● Clues in the continuum and features ● Straightforward measurement ● Focus on spectroscopic measurements ● Biased by my views, not comprehensive

Slide 3

Slide 3 text

Extinction (not Attenuation) ● Extinction – Absorption and scattering out of line- of-sight – Specific to a point source dust – Proportional to dust grain properties ● Attenuation – scattering into the line-of-sight – Varying extinction to stars – Applies to galaxies, circumstellar dust, etc.

Slide 4

Slide 4 text

Extinction vs Attenuation Github: karllark/dust_attenuation Github: karllark/dust_extinction

Slide 5

Slide 5 text

Pair Method Github: karllark/measure_extinction OB stars ideal targets Luminous at all wavelengths Fairly simple atmospheres

Slide 6

Slide 6 text

Basic measurement is the color “excess” versus wavelength (in magnitudes as we are astronomers) Normalize so measurements with different dust columns can be compared A(V) determined by extrapolating to infinite wavelength e.g.,

Slide 7

Slide 7 text

Wavelength axis variations Allow full FUV-MIR extinction to be seen Wavelength scale is proportional to energy Emphasizes UV wavelengths Versus λ Versus 1/λ UV Opt NIR MIR

Slide 8

Slide 8 text

Diffuse vs Dense Dust ● Diffuse sightlines – Most detailed extinction studies – Generally A(V) values less than a few – No 3.0 μm H2 0 ice feature ● Dense sightlines – Generally studied in the NIR/MIR only – Have 3.0 μm H2 0 ice feature Decleir et al. (2022, ApJ, 930, 15) H20 Ice

Slide 9

Slide 9 text

Our Galaxy

Slide 10

Slide 10 text

Stecher 1969, ApJ, 157, 125

Slide 11

Slide 11 text

International Ultraviolet Explorer 18 years! 100,000 UV spectra!

Slide 12

Slide 12 text

No content

Slide 13

Slide 13 text

Witt, Bohlin, & Stecher 1984, 279, 698

Slide 14

Slide 14 text

2175 A Bump width → strong variation center → almost no variation Fitzpatrick & Massa (1986, ApJ, 307, 286)

Slide 15

Slide 15 text

Fitzpatrick & Massa 1988, ApJ, 328, 734

Slide 16

Slide 16 text

No content

Slide 17

Slide 17 text

Cardelli, Clayton, & Mathis 1988, ApJ, 329L, 33; 1989, ApJ, 345, 245

Slide 18

Slide 18 text

Valencic, Clayton, & Gordon 2004, ApJ, 616, 912

Slide 19

Slide 19 text

Far Ultraviolet Spectroscopic Explorer 905 to 1195 Å

Slide 20

Slide 20 text

Pre-FUSE ORFEUS: Sasseen et al. 2002, ApJ, 566, 267 Voyager: Snow et al. 1990, ApJ, 399, L23 See also: Buss et al. 1994, France et al. 2004; Lewis et al. 2005

Slide 21

Slide 21 text

Gordon, Cartledge, & Clayton 2009, ApJ, 705, 1320 See also: Sofia et al. 2005

Slide 22

Slide 22 text

FUV (& all UV) extinction smooth Residuals due to stellar or ISM HI FUV rise consistent with feature peaking at ~800 Å

Slide 23

Slide 23 text

ApJ, 944, 33 Dries Van De Putte

Slide 24

Slide 24 text

FUV rise extinction component & H2 column Very strong correlation & goes through zero!!! Grain responsible for FUV rise and H2 co-spatial Not the case for the 2175 Å bump

Slide 25

Slide 25 text

Hubble Space Telescope 912 Å to 2.5 μm

Slide 26

Slide 26 text

Pre-Hubble ● Continuum measured mainly via photometry (until 2019) ● Except for Orion dust – Cadelli & Clayton (1988, AJ, 95, 516) ● Very Broad Structure with low resolution spectra – Next slide

Slide 27

Slide 27 text

Very Broad Structure = Deviations from linear w/ 1/lambda Whiteoak 1966, ApJ, 144, 305 Hayes 1973, IAUS, 53, 83

Slide 28

Slide 28 text

No content

Slide 29

Slide 29 text

Broad Optical Features (origin unknown) Long known Very Broad Structure Explained! Massa, Fitzpatrick, & Gordon 2020, ApJ, 891, 67 Centers: 4370, 4870, & 6300 Å Widths: ~10% Two blue correlate w/ 2175 Å

Slide 30

Slide 30 text

ISS features as strong as Si features

Slide 31

Slide 31 text

NIR/MIR Extinction ● Continuum measured via photometry (until 2021) ● Features measured spectroscopically – Narrow wavelength ranges or towards “complicated” stars – 3.0 μm ice – 3.4 μm hydrocarbon – 10/20 μm silicate features

Slide 32

Slide 32 text

● Many ice features ● 3.0 μm H2 0 the strongest ● Requires A(V) > 3 Boogert, Gerakines, & Whittet

Slide 33

Slide 33 text

3.4 μm Hydrocarbon grains Pendelton & Allamandola (2002, ApJS, 138, 75)

Slide 34

Slide 34 text

10 & 20 μm Silicate Features Chiar & Tielens (2006, ApJ, 637, 774)

Slide 35

Slide 35 text

Crystalline Silicate Fraction ~0.2% (or zero) Kemper, Viernd, & Tielens (2004, ApJ, 609, 826) Amorphous Crystalline

Slide 36

Slide 36 text

NASA’s Infrared Telescope Facility 0.8 to 20 μm

Slide 37

Slide 37 text

Marjorie Decleir

Slide 38

Slide 38 text

3.0 μm ice feature measured in the diffuse average at A(ice)/A(V) = 0.0019 +/- 0.007 Not a significant detection, but intriguing at the level predicted by Potapov et al. (2021) if ice is present in shadowed pits in silicate grains

Slide 39

Slide 39 text

Spitzer Space Telescope 3 to 160 μm

Slide 40

Slide 40 text

No content

Slide 41

Slide 41 text

No content

Slide 42

Slide 42 text

Variations and 1st Direct Comparison of UV and MIR Strong variation Not correlated with grain size Silicates not correlated with 2175 A bump

Slide 43

Slide 43 text

Putting all the Wavelengths Together

Slide 44

Slide 44 text

Cardelli, Clayton, & Mathis 1988, ApJ, 329L, 33; 1989, ApJ, 345, 245

Slide 45

Slide 45 text

ApJ, 2023, 950, 86 Based on spectroscopy at all wavelengths

Slide 46

Slide 46 text

Average Variation with R(V)/grain size

Slide 47

Slide 47 text

Average Variation with R(V)/grain size Intermediate Scale Structure (ISS) features (Massa et al. 2020) size dependent

Slide 48

Slide 48 text

Significant Deviations from Literature Curves at Select R(V)

Slide 49

Slide 49 text

MW Extinction Summary ● Spectroscopic from FUV to MIR – FUV smooth & consistent with feature peaking ~800 Å – H2 correlates with FUV rise → co-spatial! – New optical features → unknown origin – NIR powerlaw & may contain H2 0 ice – MIR versus UV features → 2175 Å not due to silicates – R(V) dependent extinction relationship → one relationship for all wavelengths

Slide 50

Slide 50 text

Beyond Milky Way

Slide 51

Slide 51 text

Clayton & Martin 1985, ApJ, 288, 558 Add figure showing 30 Dor/LMC2

Slide 52

Slide 52 text

Misselt et al. (1999, ApJ, 515, 128)

Slide 53

Slide 53 text

Spatial Distribution Misselt et al. (1999, ApJ, 515, 128)

Slide 54

Slide 54 text

Prevot et al. 1984, A&A, 132, 389

Slide 55

Slide 55 text

SMC Extinction Curves Milky Way-like! (2175 Å bump) 4 similar curves are found in the star forming bar of the SMC! Ha image of the SMC Gordon & Clayton (1998, ApJ, 500, 816) Gordon et al. (2003, ApJ, 594, 279) STIS Small region Lots of variation! Maiz Apellaniz & Rubio (2012, A&A, 541, 54)

Slide 56

Slide 56 text

ApJ, submitted (favorable referee report received)

Slide 57

Slide 57 text

Bump and Bumpless sightlines throughout galaxy

Slide 58

Slide 58 text

No content

Slide 59

Slide 59 text

SMC AzV 456 SMC Bar LMC2 LMC General MW Quiescent Processed Continuum of Properties Gordon et al. (2003, ApJ, 594, 279) Gordon et al. (2016, ApJ, 826, 104)

Slide 60

Slide 60 text

Family of Curves = Order in the Chaos

Slide 61

Slide 61 text

Hints at Origin

Slide 62

Slide 62 text

N(HI)/A(V) traces metallicity and dust formation & destruction

Slide 63

Slide 63 text

UV vs MIR Carbonaceous Features

Slide 64

Slide 64 text

Future ● HST programs – M31/M33 UV extinction (PI: Clayton) ● Extinction beyond the MW/LMC/SMC – MW expanded high/low R(V) sample (PI: Decleir) ● Dust at the extremes ● JWST programs – WISCI (PI: Zeegers) & MEAD (PI: Decleir) ● NIR/MIR MW extinction continuum and features – LMC/SMC MIR spectra (PI: Gordon) ● 1st measure of 10 um silicate feature

Slide 65

Slide 65 text

Summary ● Extinction measurements are fundamental to understanding dust ● Spectroscopic Extinction from FUV to MIR in Milky Way – Overall quite smooth – Few broad features including some new ones – One intriguing correlation, comforting non-correlations – One R(V) extinction relationship for all wavelengths ● Beyond Milky Way – Larger variation than MW – Intriguing correlations including with gas-to-dust [N(HI)/A(V)] ● Future bright w/ HST and JWST

Slide 66

Slide 66 text

Thanks ● Hiring in Fall – postdoc position, PDRs, extinction curves, dust, ... ● Interested in joint PhD between UGhent and STScI? – Talk to me ● Interested in jobs at STScI – Talk to me

Slide 67

Slide 67 text

Not Covered ● Dust scattering (e.g., albedo, scattering phase function) ● Dust emission ● Atomic composition of dust (e.g., depletions) ● Dust Polarization ● Dust grain models ● All important, but not the focus of this talk

Slide 68

Slide 68 text

Slide from Dries Van De Puttte

Slide 69

Slide 69 text

No content

Slide 70

Slide 70 text

Fitting a Line New penultimate technique* ● Fully accounts for correlated x & y uncertainties ● Allows for different correlations for different points ● Line integral for each point *The authors await knowledge of what paper to cite where this was first shown (likely from decades ago)

Slide 71

Slide 71 text

Clayton et al. 2015, ApJ, 815, 14

Slide 72

Slide 72 text

M31/M33 Extinction ● All w/ HST photometry ● Similar to MW ● Radial, metallicity variations ● 1st measurements in M33 Petia Yanchulov Merica-Jones Clayton, G et al., in prep Yanchulov M-J, P. et al., in prep

Slide 73

Slide 73 text

Dust Grain Modeling: DGFit ● Calculate how new observations affect our understanding of dust grains – Bigger/smaller grains, new materials, ... ● Quantitatively compare different grain materials – Bayesian statistics may help (Bayes Factors) ● Add additional observational constraints – E.g., albedo and scattering phase function asymmetry

Slide 74

Slide 74 text

MW Average

Slide 75

Slide 75 text

SMC Average

Slide 76

Slide 76 text

MEAD Measuring Extinction and Abundances of Dust ● Cycle 1 JWST program ● PI: Marjorie Decleir (+3 Co-Is) ● 9 sightlines towards OB stars – A(V) from 1.2 – 2.5 mag – R(V) from 2.5 – 5, f(H2) from 0.1 – 0.65 – Preexisting HST program, PI: M. Decleir ● Measure atomic dust abundances via UV spectra of ISM gas absorption lines (depletions) – Existing IUE spectra, optical/NIR photometry ● MIRI MRS + NIRCam grism observations – Spectra from 2.5 – 4 µm & 5 – 28 µm ● Data status – NIRCam grism: all data taken – MIRI MRS: 6/9 sightlines (just!) taken

Slide 77

Slide 77 text

WISCI Webb Investigation of Silicates, Carbon, and Ices ● Cycle 1 JWST program ● PI: Sascha Zeegers (+20 Co-Is) ● 12 sightlines towards OB stars – A(V) from 4.6 – 8.1 mag – Followup HST program, PI: Zeegers ● Optical for all, UV for half ● MIRI MRS + NIRCam grism observations – Complete spectrum from 2.5 – 28 µm ● Data status – NIRCam grism: all data taken – MIRI MRS: 8/12 sightlines taken Sascha Zeegers ESA

Slide 78

Slide 78 text

NIRCam+MIRI spectra WISCI MEAD RJ = Rayleigh-Jeans Units = λ2 F(ν)

Slide 79

Slide 79 text

JWST 10 µm Silicate Feature WISCI MEAD Models A(V) = 5 Models A(V) = 1.5 Si-O Si-O Huα Pfα/Huβ Huα Pfα/Huβ

Slide 80

Slide 80 text

JWST NIR Carbonaceous Features WISCI MEAD Models A(V) = 5 Models A(V) = 1.5 C-H (aliphatic) C-H (aromatic) C-H (aliphatic) C-H (aromatic) Pfδ Pfε Pfγ Pfδ Pfε Pfγ

Slide 81

Slide 81 text

JWST MIR Carbonaceous Features I WISCI MEAD Models A(V) = 5 Models A(V) = 1.5 C=C C=O C=C=C? C-H? C=C C=O C=C=C? C-H? C=C=C C=O: carbonyl C=C: aromatic/olefinic C-H: aromatic/aliphatic Huγ Huγ