skip to main content


Title: Starlight-polarization-based tomography of the magnetized ISM: PASIPHAE’s line-of-sight inversion method
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
Award ID(s):
2109127
NSF-PAR ID:
10465505
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; ; « less
Date Published:
Journal Name:
Astronomy & Astrophysics
Volume:
670
ISSN:
0004-6361
Page Range / eLocation ID:
A164
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Context. Atomic gas in the diffuse interstellar medium (ISM) is organized in filamentary structures. These structures usually host cold and dense molecular clumps. The Galactic magnetic field is considered to play an important role in the formation of these clumps. Aims. Our goal is to explore the role of the magnetic field in the H I -H 2 transition process. Methods. We targeted a diffuse ISM filamentary cloud toward the Ursa Major cirrus where gas transitions from atomic to molecular. We probed the magnetic field properties of the cloud with optical polarization observations. We performed multiwavelength spectroscopic observations of different species in order to probe the gas phase properties of the cloud. We observed the CO ( J = 1−0) and ( J = 2−1) lines in order to probe the molecular content of the cloud. We also obtained observations of the [C ii ] 157.6 µ m emission line in order to trace the CO-dark H 2 gas and estimate the mean volume density of the cloud. Results. We identified two distinct subregions within the cloud. One of the regions is mostly atomic, while the other is dominated by molecular gas, although most of it is CO-dark. The estimated plane-of-the-sky magnetic field strength between the two regions remains constant within uncertainties and lies in the range 13–30 µG. The total magnetic field strength does not scale with density. This implies that gas is compressed along the field lines. We also found that turbulence is trans-Alfvénic, with M A ≈ 1. In the molecular region, we detected an asymmetric CO clump whose minor axis is closer, with a 24° deviation, to the mean magnetic field orientation than the angle of its major axis. The H i velocity gradients are in general perpendicular to the mean magnetic field orientation except for the region close to the CO clump, where they tend to become parallel. This phenomenon is likely related to gas undergoing gravitational infall. The magnetic field morphology of the target cloud is parallel to the H i column density structure of the cloud in the atomic region, while it tends to become perpendicular to the H i structure in the molecular region. On the other hand, the magnetic field morphology seems to form a smaller offset angle with the total column density shape (including both atomic and molecular gas) of this transition cloud. Conclusions. In the target cloud where the H i –H 2 transition takes place, turbulence is trans-Alfvénic, and hence the magnetic field plays an important role in the cloud dynamics. Atomic gas probably accumulates preferentially along the magnetic field lines and creates overdensities where molecular gas can form. The magnetic field morphology is probed better by the total column density shape of the cloud, and not its H i column density shape. 
    more » « less
  2. null (Ed.)
    Abstract We present measurements of cosmic shear two-point correlation functions (TPCFs) from Hyper Suprime-Cam Subaru Strategic Program (HSC) first-year data, and derive cosmological constraints based on a blind analysis. The HSC first-year shape catalog is divided into four tomographic redshift bins ranging from $z=0.3$ to 1.5 with equal widths of $\Delta z =0.3$. The unweighted galaxy number densities in each tomographic bin are 5.9, 5.9, 4.3, and $2.4\:$arcmin$^{-2}$ from the lowest to highest redshifts, respectively. We adopt the standard TPCF estimators, $\xi _\pm$, for our cosmological analysis, given that we find no evidence of significant B-mode shear. The TPCFs are detected at high significance for all 10 combinations of auto- and cross-tomographic bins over a wide angular range, yielding a total signal-to-noise ratio of 19 in the angular ranges adopted in the cosmological analysis, $7^{\prime }<\theta <56^{\prime }$ for $\xi _+$ and $28^{\prime }<\theta <178^{\prime }$ for $\xi _-$. We perform the standard Bayesian likelihood analysis for cosmological inference from the measured cosmic shear TPCFs, including contributions from intrinsic alignment of galaxies as well as systematic effects from PSF model errors, shear calibration uncertainty, and source redshift distribution errors. We adopt a covariance matrix derived from realistic mock catalogs constructed from full-sky gravitational lensing simulations that fully account for survey geometry and measurement noise. For a flat $\Lambda$ cold dark matter model, we find $S\,_8 \equiv \sigma _8\sqrt{\Omega _{\rm m}/0.3}=0.804_{-0.029}^{+0.032}$, and $\Omega _{\rm m}=0.346_{-0.100}^{+0.052}$. We carefully check the robustness of the cosmological results against astrophysical modeling uncertainties and systematic uncertainties in measurements, and find that none of them has a significant impact on the cosmological constraints. 
    more » « less
  3. Abstract Modifications of the matter power spectrum due to baryonic physics are one of the major theoretical uncertainties in cosmological weak lensing measurements. Developing robust mitigation schemes for this source of systematic uncertainty increases the robustness of cosmological constraints, and may increase their precision if they enable the use of information from smaller scales. Here we explore the performance of two mitigation schemes for baryonic effects in weak lensing cosmic shear: the principal component analysis (PCA) method and the halo-model approach in hmcode. We construct mock tomographic shear power spectra from four hydrodynamical simulations, and run simulated likelihood analyses with cosmolike assuming LSST-like survey statistics. With an angular scale cut of ℓmax < 2000, both methods successfully remove the biases in cosmological parameters due to the various baryonic physics scenarios, with the PCA method causing less degradation in the parameter constraints than hmcode. For a more aggressive ℓmax = 5000, the PCA method performs well for all but one baryonic physics scenario, requiring additional training simulations to account for the extreme baryonic physics scenario of Illustris; hmcode exhibits tensions in the 2D posterior distributions of cosmological parameters due to lack of freedom in describing the power spectrum for $k \gt 10\ h^{-1}\, \mathrm{Mpc}$. We investigate variants of the PCA method and improve the bias mitigation through PCA by accounting for the noise properties in the data via Cholesky decomposition of the covariance matrix. Our improved PCA method allows us to retain more statistical constraining power while effectively mitigating baryonic uncertainties even for a broad range of baryonic physics scenarios. 
    more » « less
  4. ABSTRACT

    Observing 3D magnetic fields, including orientation and strength, within the interstellar medium is vital but notoriously difficult. However, recent advances in our understanding of anisotropic magnetohydrodynamic (MHD) turbulence demonstrate that MHD turbulence and 3D magnetic fields leave their imprints on the intensity features of spectroscopic observations. Leveraging these theoretical frameworks, we propose a novel Convolutional Neural Network (CNN) model to extract this embedded information, enabling the probe of 3D magnetic fields. This model examines the plane-of-the-sky magnetic field orientation (ϕ), the magnetic field’s inclination angle (γ) relative to the line-of-sight, and the total magnetization level (M$_{\rm A}^{-1}$) of the cloud. We train the model using synthetic emission lines of 13CO (J  = 1–0) and C18O (J  = 1–0), generated from 3D MHD simulations that span conditions from sub-Alfvénic to super-Alfvénic molecular clouds. Our tests confirm that the CNN model effectively reconstructs the 3D magnetic field topology and magnetization. The median uncertainties are under 5° for both ϕ and γ, and less than 0.2 for MA in sub-Alfvénic conditions (MA ≈ 0.5). In super-Alfvénic scenarios (MA ≈ 2.0), they are under 15° for ϕ and γ, and 1.5 for MA. We applied this trained CNN model to the L1478 molecular cloud. Results show a strong agreement between the CNN-predicted magnetic field orientation and that derived from Planck 353 GHz polarization. The CNN approach enabled us to construct the 3D magnetic field map for L1478, revealing a global inclination angle of ≈76° and a global MA of ≈1.07.

     
    more » « less
  5. Context. Molecular filaments and hubs have received special attention recently thanks to new studies showing their key role in star formation. While the (column) density and velocity structures of both filaments and hubs have been carefully studied, their magnetic field (B-field) properties have yet to be characterized. Consequently, the role of B-fields in the formation and evolution of hub-filament systems is not well constrained. Aims. We aim to understand the role of the B-field and its interplay with turbulence and gravity in the dynamical evolution of the NGC 6334 filament network that harbours cluster-forming hubs and high-mass star formation. Methods. We present new observations of the dust polarized emission at 850 μ m toward the 2 pc × 10 pc map of NGC 6334 at a spatial resolution of 0.09 pc obtained with the James Clerk Maxwell Telescope (JCMT) as part of the B-field In STar-forming Region Observations (BISTRO) survey. We study the distribution and dispersion of the polarized intensity ( PI ), the polarization fraction ( PF ), and the plane-of-the-sky B-field angle ( χ B_POS ) toward the whole region, along the 10 pc-long ridge and along the sub-filaments connected to the ridge and the hubs. We derived the power spectra of the intensity and χ B POS along the ridge crest and compared them with the results obtained from simulated filaments. Results. The observations span ~3 orders of magnitude in Stokes I and PI and ~2 orders of magnitude in PF (from ~0.2 to ~ 20%). A large scatter in PI and PF is observed for a given value of I . Our analyses show a complex B-field structure when observed over the whole region (~ 10 pc); however, at smaller scales (~1 pc), χ B POS varies coherently along the crests of the filament network. The observed power spectrum of χ B POS can be well represented with a power law function with a slope of − 1.33 ± 0.23, which is ~20% shallower than that of I . We find that this result is compatible with the properties of simulated filaments and may indicate the physical processes at play in the formation and evolution of star-forming filaments. Along the sub-filaments, χ B POS rotates frombeing mostly perpendicular or randomly oriented with respect to the crests to mostly parallel as the sub-filaments merge with the ridge and hubs. This variation of the B-field structure along the sub-filaments may be tracing local velocity flows of infalling matter in the ridge and hubs. Our analysis also suggests a variation in the energy balance along the crests of these sub-filaments, from magnetically critical or supercritical at their far ends to magnetically subcritical near the ridge and hubs. We also detect an increase in PF toward the high-column density ( N H 2 ≳ 10 23  cm −2 ) star cluster-forming hubs. These latter large PF values may be explained by the increase in grain alignment efficiency due to stellar radiation from the newborn stars, combined with an ordered B-field structure. Conclusions. These observational results reveal for the first time the characteristics of the small-scale (down to ~ 0.1 pc) B-field structure of a 10 pc-long hub-filament system. Our analyses show variations in the polarization properties along the sub-filaments that may be tracing the evolution of their physical properties during their interaction with the ridge and hubs. We also detect an impact of feedback from young high-mass stars on the local B-field structure and the polarization properties, which could put constraints on possible models for dust grain alignment and provide important hints as to the interplay between the star formation activity and interstellar B-fields. 
    more » « less