We present the first Bayesian method for tomographic decomposition of the plane-of-sky orientation of the magnetic field with the use of stellar polarimetry and distance. This standalone tomographic inversion method presents an important step forward in reconstructing the magnetized interstellar medium (ISM) in three dimensions within dusty regions. We develop a model in which the polarization signal from the magnetized and dusty ISM is described by thin layers at various distances, a working assumption which should be satisfied in small-angular circular apertures. Our modeling makes it possible to infer the mean polarization (amplitude and orientation) induced by individual dusty clouds and to account for the turbulence-induced scatter in a generic way. We present a likelihood function that explicitly accounts for uncertainties in polarization and parallax. We develop a framework for reconstructing the magnetized ISM through the maximization of the log-likelihood using a nested sampling method. We test our Bayesian inversion method on mock data, representative of the high Galactic latitude sky, taking into account realistic uncertainties from Gaia and as expected for the optical polarization survey P ASIPHAE according to the currently planned observing strategy. We demonstrate that our method is effective at recovering the cloud properties as soon as the polarization induced by a cloud to its background stars is higher than ~0.1% for the adopted survey exposure time and level of systematic uncertainty. The larger the induced polarization is, the better the method’s performance, and the lower the number of required stars. Our method makes it possible to recover not only the mean polarization properties but also to characterize the intrinsic scatter, thus creating new ways to characterize ISM turbulence and the magnetic field strength. Finally, we apply our method to an existing data set of starlight polarization with known line-of-sight decomposition, demonstrating agreement with previous results and an improved quantification of uncertainties in cloud properties.
more »
« less
The first degree-scale starlight-polarization-based tomography map of the magnetized interstellar medium
We present the first degree-scale tomography map of the dusty magnetized interstellar medium (ISM) from stellar polarimetry and distance measurements. We used the RoboPol polarimeter at Skinakas Observatory to conduct a survey of the polarization of starlight in a region of the sky of about four square degrees. We propose a Bayesian method to decompose the stellar-polarization source field along the distance to invert the three-dimensional (3D) volume occupied by the observed stars. We used this method to obtain the first 3D map of the dusty magnetized ISM. Specifically, we produced a tomography map of the orientation of the plane-of-sky component of the magnetic field threading the diffuse, dusty regions responsible for the stellar polarization. For the targeted region centered on Galactic coordinates (l,b) ≈ (103.3°, 22.3°), we identified several ISM clouds. Most of the lines of sight intersect more than one cloud. A very nearby component was detected in the foreground of a dominant component from which most of the polarization signal comes and which we identified as being an intersection of the wall of the Local Bubble and the Cepheus Flare. Farther clouds, with a distance of up to 2 kpc, were similarly detected. Some of them likely correspond to intermediate-velocity clouds seen in HIspectra in this region of the sky. We found that the orientation of the plane-of-sky component of the magnetic field changes along distance for most of the lines of sight. Our study demonstrates that starlight polarization data coupled to distance measures have the power to reveal the great complexity of the dusty magnetized ISM in 3D and, in particular, to provide local measurements of the plane-of-sky component of the magnetic field in dusty regions. This demonstrates that the inversion of large data volumes, as expected from the PASIPHAEsurvey, will provide the necessary means to move forward in the modeling of the Galactic magnetic field and of the dusty magnetized ISM as a contaminant in observations of the cosmic microwave background polarization.
more »
« less
- Award ID(s):
- 2109127
- PAR ID:
- 10547710
- Author(s) / Creator(s):
- ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more »
- Publisher / Repository:
- EDP Sciences
- Date Published:
- Journal Name:
- Astronomy & Astrophysics
- Volume:
- 684
- ISSN:
- 0004-6361
- Page Range / eLocation ID:
- A162
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract We characterize Galactic dust filaments by correlating BICEP/Keck and Planck data with polarization templates based on neutral hydrogen (Hi) observations. Dust polarization is important for both our understanding of astrophysical processes in the interstellar medium (ISM) and the search for primordial gravitational waves in the cosmic microwave background (CMB). In the diffuse ISM, Hiis strongly correlated with the dust and partly organized into filaments that are aligned with the local magnetic field. We analyze the deep BICEP/Keck data at 95, 150, and 220 GHz, over the low-column-density region of sky where BICEP/Keck has set the best limits on primordial gravitational waves. We separate the Hiemission into distinct velocity components and detect dust polarization correlated with the local Galactic Hibut not with the Hiassociated with Magellanic Streami. We present a robust, multifrequency detection of polarized dust emission correlated with the filamentary Himorphology template down to 95 GHz. For assessing its utility for foreground cleaning, we report that the Himorphology template correlates inBmodes at a ∼10%–65% level over the multipole range 20 <ℓ< 200 with the BICEP/Keck maps, which contain contributions from dust, CMB, and noise components. We measure the spectral index of the filamentary dust component spectral energy distribution to beβ= 1.54 ± 0.13. We find no evidence for decorrelation in this region between the filaments and the rest of the dust field or from the inclusion of dust associated with the intermediate velocity Hi. Finally, we explore the morphological parameter space in the Hi-based filamentary model.more » « less
-
Abstract We adopt the deep learning methodcasi-3d(convolutional approach to structure identification-3D) to infer the orientation of magnetic fields in sub-/trans-Alfvénic turbulent clouds from molecular line emission. We carry out magnetohydrodynamic simulations with different magnetic field strengths and use these to generate synthetic observations. We apply the 3D radiation transfer coderadmc-3dto model12CO and13CO (J = 1−0) line emission from the simulated clouds and then train acasi-3dmodel on these line emission data cubes to predict magnetic field morphology at the pixel level. The trainedcasi-3dmodel is able to infer magnetic field directions with a low error (≲10° for sub-Alfvénic samples and ≲30° for trans-Alfvénic samples). We further test the performance ofcasi-3don a real sub-/trans- Alfvénic region in Taurus. Thecasi-3dprediction is consistent with the magnetic field direction inferred from Planck dust polarization measurements. We use our developed methods to produce a new magnetic field map of Taurus that has a three times higher angular resolution than the Planck map.more » « less
-
Abstract Cool, dusty interstellar material plays an important role in the chemical evolution of galaxies. We present an analysis of this material across galaxy type through a spatially resolved spectral stacking analysis of galaxies from the MaNGA survey. With stellar population synthesis, we isolate neutral gas signals from resonance lines, comparing outcomes across model types, galactic geometry, and host stellar mass and age. We find that both synthetic and empirical models fail to capture the range of galactic chemical abundances. There is also notable Naicontamination from the Galaxy’s interstellar medium (ISM) in the MILES empirical stellar library. We are unable to reliably determine the column density of the gas due to the accuracy of absorption measurements, but differential analysis across radius and inclination reveals consistent and significant path-length dependent absorption in the equivalent width of Nai. We note similar but lesser trends in a narrow Caiiindex. We find no trends in Caior in a broad Caiiindex, indicating its ISM insensitivity and providing evidence in favor of its utility in determining the age and chemical content of stellar populations. Our data shows there is a cool ISM component in most external galaxies withDn(4000) < 1.7 that can be traced by Nai. Lastly, we caution that the characterization of gas kinematics traced by Naiin such low-resolution spectra is subject to systematic effects due to the chosen approach to stellar population modeling.more » « less
-
Abstract We assess the impact of Caiiλλ3934, 3969 and Naiλλ5891, 5897 absorption arising in the interstellar medium (ISM) on the Sloan Digital Sky Survey-IV MaNGA Stellar Library (MaStar) and produce corrected spectroscopy for 80% of the 24,162-star catalog. We model the absorption strength of these transitions as a function of the stellar distance, Galactic latitude, and dust reddening based on high-spectral resolution studies. With this model, we identify 6342 MaStar stars that have negligible ISM absorption (WISM(CaiiK) < 0.07 Å andWISM(Nai5891) < 0.05 Å). For 12,110 of the remaining stars, we replace their NaiD profile (and their Caiiprofile for effective temperaturesTeff> 9000 K) with a coadded spectrum of low-ISM stars with similarTeff, surface gravity, and metallicity. For 738 additional stars withTeff> 9000 K, we replace these spectral regions with a matching ATLAS9-based BOSZ model. This results in a mean reduction inW(CaiiK) (W(NaiD)) of 0.4–0.7 Å (0.6–1.1 Å) for hot stars (Teff> 7610 K), and a mean reduction inW(NaiD) of 0.1–0.2 Å for cooler stars. We show that interstellar absorption in the simple stellar population (SSP) model spectra constructed from the original library artificially enhancesW(CaiiK) by ≳20% at young ages (<400 Myr); dramatically enhances the strength of stellar NaiD in starbursting systems (by ≳50%); and enhances stellar NaiD in older stellar populations (≳10 Gyr) by ≳10%. We provide SSP spectra constructed from the cleaned library and discuss the implications of these effects for stellar population synthesis analyses constraining the stellar age, [Na/Fe] abundance, and initial mass function.more » « less
An official website of the United States government

