We characterize Galactic dust filaments by correlating BICEP/Keck and Planck data with polarization templates based on neutral hydrogen (H
Observing in six frequency bands from 27 to 280 GHz over a large sky area, the Simons Observatory (SO) is poised to address many questions in Galactic astrophysics in addition to its principal cosmological goals. In this work, we provide quantitative forecasts on astrophysical parameters of interest for a range of Galactic science cases. We find that SO can: constrain the frequency spectrum of polarized dust emission at a level of Δ A supplement describing author contributions to this paper can be found at
- NSF-PAR ID:
- 10479408
- Author(s) / Creator(s):
- ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more »
- Publisher / Repository:
- DOI PREFIX: 10.3847
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 929
- Issue:
- 2
- ISSN:
- 0004-637X
- Format(s):
- Medium: X Size: Article No. 166
- Size(s):
- ["Article No. 166"]
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract i ) 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, Hi is 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 Hi emission into distinct velocity components and detect dust polarization correlated with the local Galactic Hi but not with the Hi associated with Magellanic Streami . We present a robust, multifrequency detection of polarized dust emission correlated with the filamentary Hi morphology template down to 95 GHz. For assessing its utility for foreground cleaning, we report that the Hi morphology template correlates inB modes 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
Water vapor (H2O) is one of the brightest molecular emitters after carbon monoxide (CO) in galaxies with high infrared (IR) luminosity, allowing us to investigate the warm and dense phase of the interstellar medium (ISM) where star formation occurs. However, due to the complexity of its radiative spectrum, H2O is not frequently exploited as an ISM tracer in distant galaxies. Therefore, H2O studies of the warm and dense gas at high-
z remain largely unexplored. In this work, we present observations conducted with the Northern Extended Millimeter Array (NOEMA) toward threez > 6 IR-bright quasarsJ2310+1855 ,J1148+5251 , andJ0439+1634 targeted in their multiple para- and ortho-H2O transitions (312 − 303, 111 − 000, 220 − 211, and 422 − 413), as well as their far-IR (FIR) dust continuum. By combining our data with previous measurements from the literature, we estimated the dust masses and temperatures, continuum optical depths, IR luminosities, and star formation rates (SFR) from the FIR continuum. We modeled the H2O lines using the MOLPOP-CEP radiative transfer code, finding that water vapor lines in our quasar host galaxies are primarily excited in the warm, dense (with a gas kinetic temperature and density ofT kin = 50 K,n H2 ∼ 104.5 − 105 cm−3) molecular medium with a water vapor column density ofN H2O ∼ 2 × 1017 − 3 × 1018 cm−3. High-J H2O lines are mainly radiatively pumped by the intense optically-thin far-IR radiation field associated with a warm dust component at temperatures ofT dust ∼ 80 − 190 K that account for < 5 − 10% of the total dust mass. In the case of J2310+1855, our analysis points to a relatively high value of the continuum optical depth at 100 μm (τ 100 ∼ 1). Our results are in agreement with expectations based on the H2O spectral line energy distribution of local and high-z ultra-luminous IR galaxies and active galactic nuclei (AGN). The analysis of the Boltzmann diagrams highlights the interplay between collisions and IR pumping in populating the high H2O energy levels and it allows us to directly compare the excitation conditions in the targeted quasar host galaxies. In addition, the observations enable us to sample the high-luminosity part of the H2O–total-IR (TIR) luminosity relations (L H2O −L TIR). Overall, our results point to supralinear trends that suggest H2O–TIR relations are likely driven by IR pumping, rather than the mere co-spatiality between the FIR continuum- and line-emitting regions. The observedL H2O/L TIRratios in ourz > 6 quasars do not show any strong deviations with respect to those measured in star-forming galaxies and AGN at lower redshifts. This supports the notion that H2O can be likely used to trace the star formation activity buried deep within the dense molecular clouds. -
more » « less
Abstract High-fidelity simulators that connect theoretical models with observations are indispensable tools in many sciences. If the likelihood is known, inference can proceed using standard techniques. However, when the likelihood is intractable or unknown, a simulator makes it possible to infer the parameters of a theoretical model directly from real and simulated observations when coupled with machine learning. We introduce an extension of the recently proposed likelihood-free frequentist inference (
LF2I ) approach that makes it possible to construct confidence sets with thep -value function and to use the same function to check the coverage explicitly at any given parameter point. LikeLF2I , this extension yields provably valid confidence sets in parameter inference problems for which a high-fidelity simulator is available. The utility of our algorithm is illustrated by applying it to three pedagogically interesting examples: the first is from cosmology, the second from high-energy physics and astronomy, both with tractable likelihoods, while the third, with an intractable likelihood, is from epidemiology3Code to reproduce all of our results is available on
https://github.com/AliAlkadhim/ALFFI . -
more » « less
Abstract We present the first linear polarization measurements from the 2015 long-duration balloon flight of
Spider , which is an experiment that is designed to map the polarization of the cosmic microwave background (CMB) on degree angular scales. The results from these measurements include maps and angular power spectra from observations of 4.8% of the sky at 95 and 150 GHz, along with the results of internal consistency tests on these data. While the polarized CMB anisotropy from primordial density perturbations is the dominant signal in this region of sky, Galactic dust emission is also detected with high significance. Galactic synchrotron emission is found to be negligible in theSpider bands. We employ two independent foreground-removal techniques to explore the sensitivity of the cosmological result to the assumptions made by each. The primary method uses a dust template derived fromPlanck data to subtract the Galactic dust signal. A second approach, which constitutes a joint analysis ofSpider andPlanck data in the harmonic domain, assumes a modified-blackbody model for the spectral energy distribution of the dust with no constraint on its spatial morphology. Using a likelihood that jointly samples the template amplitude andr parameter space, we derive 95% upper limits on the primordial tensor-to-scalar ratio from Feldman–Cousins and Bayesian constructions, findingr < 0.11 andr < 0.19, respectively. Roughly half the uncertainty inr derives from noise associated with the template subtraction. New data at 280 GHz fromSpider ’s second flight will complement thePlanck polarization maps, providing powerful measurements of the polarized Galactic dust emission. -
more » « less
Abstract Reliable studies of the long-term dynamics of planetary systems require numerical integrators that are accurate and fast. The challenge is often formidable because the chaotic nature of many systems requires relative numerical error bounds at or close to machine precision (∼10−16, double-precision arithmetic); otherwise, numerical chaos may dominate over physical chaos. Currently, the speed/accuracy demands are usually only met by symplectic integrators. For example, the most up-to-date long-term astronomical solutions for the solar system in the past (widely used in, e.g., astrochronology and high-precision geological dating) have been obtained using symplectic integrators. However, the source codes of these integrators are unavailable. Here I present the symplectic integrator
orbitN (lean version 1.0) with the primary goal of generating accurate and reproducible long-term orbital solutions for near-Keplerian planetary systems (here the solar system) with a dominant massM 0. Among other features,orbitN-1.0 includesM 0’s quadrupole moment, a lunar contribution, and post-Newtonian corrections (1PN) due toM 0(fast symplectic implementation). To reduce numerical round-off errors, Kahan compensated summation was implemented. I useorbitN to provide insight into the effect of various processes on the long-term chaos in the solar system. Notably, 1PN corrections have the opposite effect on chaoticity/stability on a 100 Myr versus Gyr timescale. For the current application,orbitN is about as fast as or faster (factor 1.15–2.6) than comparable integrators, depending on hardware.1The orbitN source code (C) is available at
http://github.com/rezeebe/orbitN .
