New Interstellar Extinction Maps Based on Gaia and Other Sky Surveys

We present a new two-dimensional (2D) map of total Galactic extinction, $A_\mathrm{V}$, across the entire dust half-layer from the Sun to extragalactic space for Galactic latitudes $|b| > 13\degr$, as well as a three-dimensional (3D) map of $A_\mathrm{V}$ within 2~kpc of the Sun.
These maps are based on $A_\mathrm{V}$ and distance estimates derived from a dataset, which utilizes {\it Gaia} Data Release 3 parallaxes and multi-band photometry for nearly 100 million dwarf stars. We apply our own corrections to account for significant systematics in this dataset. Our 2D map achieves an angular resolution of 6.1~arcmin, while the 3D map offers a transverse resolution of 3.56~pc --- corresponding to variable angular resolution depending on distance --- and a radial resolution of 50~pc. In constructing these maps, we pay particular attention to the solar neighborhood (within 200~pc) and to high Galactic latitudes.
The 3D map predicts $A_\mathrm{V}$ from the Sun to any extended object within the Galactic dust layer with an accuracy of $\sigma(A_\mathrm{V}) = 0.1$~mag. The 2D map provides $A_\mathrm{V}$ estimates for the entire dust half-layer up to extragalactic distances with an accuracy of $\sigma(A_\mathrm{V}) = 0.07$~mag. We provide $A_\mathrm{V}$ estimates from our maps for various classes of extended celestial objects \textbf{with angular sizes primarily in the range of 2--40~arcmin}, including 19,809 galaxies and quasars, 170 Galactic globular clusters, 458 open clusters, \textbf{and several hundreds molecular clouds from two lists}. We also present extinction values for 8,293 Type Ia supernovae. Comparison of our extinction estimates with those from previous maps and literature sources reveals systematic differences, indicating large-scale spatial variations in the extinction law and suggesting that earlier 2D reddening maps based on infrared dust emission tend to underestimate low extinction values.

Three dimensional interstellar extinction maps provide a powerful tool for stellar population analysis. We use data from the VISTA Variables in the Via Lactea survey together with the Besan\c{c}on stellar population synthesis model of the Galaxy to determine interstellar extinction as a function of distance in the Galactic bulge covering $ -10 < l < 10$ and $-10 < b <5$. We adopted a recently developed method to calculate the colour excess. First we constructed the H-Ks vs. Ks and J-Ks vs. Ks colour-magnitude diagrams based on the VVV catalogues that matched 2MASS. Then, based on the temperature-colour relation for M giants and the distance-colour relations, we derived the extinction as a function of distance. The observed colours were shifted to match the intrinsic colours in the Besan\c{c}on model as a function of distance iteratively. This created an extinction map with three dimensions: two spatial and one distance dimension along each line of sight towards the bulge. We present a 3D extinction map that covers the whole VVV area with a resolution of 6' x 6', using distance bins of 0.5 kpc. The high resolution and depth of the photometry allows us to derive extinction maps for a range of distances up to 10 kpc and up to 30 magnitudes of extinction in $A_{V}$. Integrated maps show the same dust features and consistent values as other 2D maps. We discuss the spatial distribution of dust features in the line of sight, which suggests that there is much material in front of the Galactic bar, specifically between 5-7 kpc. We compare our dust extinction map with high-resolution $\rm ^{12}CO$ maps towards the Galactic bulge, where we find a good correlation between $\rm ^{12}CO$ and $\rm A_{V}$. We determine the X factor by combining the CO map and our dust extinction map. Our derived average value is consistent with the canonical value of the Milky Way.

Aims: We used 12,530 photometrically-selected blue horizontal branch (BHB) stars from the Sloan Digital Sky Survey (SDSS) to estimate the total extinction of the Milky Way at the high Galactic latitudes, $R_V$ and $A_V$ in each line of sight. Methods: A Bayesian method was developed to estimate the reddening values in the given lines of sight. Based on the most likely values of reddening in multiple colors, we were able to derive the values of $R_V$ and $A_V$. Results: We selected 94 zero-reddened BHB stars from seven globular clusters as the template. The reddening in the four SDSS colors for the northern Galactic cap were estimated by comparing the field BHB stars with the template stars. The accuracy of this estimation is around 0.01\,mag for most lines of sight. We also obtained $<R_V>$ to be around 2.40$\pm1.05$ and $A_V$ map within an uncertainty of 0.1\,mag. The results, including reddening values in the four SDSS colors, $A_V$, and $R_V$ in each line of sight, are released on line. In this work, we employ an up-to-date parallel technique on GPU card to overcome time-consuming computations. We plan to release online the C++ CUDA code used for this analysis. Conclusions: The extinction map derived from BHB stars is highly consistent with that from Schlegel, Finkbeiner & Davis(1998). The derived $R_V$ is around 2.40$\pm1.05$. The contamination probably makes the $R_V$ be larger.

Fast radio bursts (FRBs) are an emerging class of bright, highly dispersed radio pulses. Recent work by Thornton et al. (2013) has revealed a population of FRBs in the High Time Resolution Universe (HTRU) survey at high Galactic latitudes. A variety of progenitors have been proposed including cataclysmic events at cosmological distances, Galactic flare stars, and terrestrial radio frequency interference. Here we report on a search for FRBs at intermediate Galactic latitudes ($-15^{\circ}$ $< b <$ 15$^{\circ}$) in data taken as part of the HTRU survey. No FRBs were discovered in this region. Several effects such as dispersion, scattering, sky temperature and scintillation decrease the sensitivity by more than 3$\sigma$ in $\sim$20\% of survey pointings. Including all of these effects, we exclude the hypothesis that FRBs are uniformly distributed on the sky with 99\% confidence. This low probability implies that additional factors -- not accounted for by standard Galactic models -- must be included to ease the discrepancy between the detection rates at high and low Galactic latitudes. A revised rate estimate or another strong and heretofore unknown selection effect in Galactic latitude would provide closer agreement between the surveys' detection rates. The dearth of detections at low Galactic latitude disfavors a Galactic origin for these bursts.

Observations of the submillimetre emission from Galactic dust, in both total intensityIand polarization, have received tremendous interest thanks to thePlanckfull-sky maps. In this paper we make use of such full-sky maps of dust polarized emission produced from the third public release ofPlanckdata. As the basis for expanding on astrophysical studies of the polarized thermal emission from Galactic dust, we present full-sky maps of the dust polarization fractionp, polarization angleψ, and dispersion function of polarization angles 𝒮. The joint distribution (one-point statistics) ofpandNHconfirms that the mean and maximum polarization fractions decrease with increasingNH. The uncertainty on the maximum observed polarization fraction,pmax= 22.0−1.4+3.5% at 353 GHz and 80′ resolution, is dominated by the uncertainty on the Galactic emission zero level in total intensity, in particular towards diffuse lines of sight at high Galactic latitudes. Furthermore, the inverse behaviour betweenpand 𝒮 found earlier is seen to be present at high latitudes. This follows the 𝒮 ∝ p−1relationship expected from models of the polarized sky (including numerical simulations of magnetohydrodynamical turbulence) that include effects from only the topology of the turbulent magnetic field, but otherwise have uniform alignment and dust properties. Thus, the statistical properties ofp,ψ, and 𝒮 for the most part reflect the structure of the Galactic magnetic field. Nevertheless, we search for potential signatures of varying grain alignment and dust properties. First, we analyse the product map 𝒮 × p, looking for residual trends. While the polarization fractionpdecreases by a factor of 3−4 betweenNH = 1020 cm−2andNH = 2 × 1022 cm−2, out of the Galactic plane, this product 𝒮 × ponly decreases by about 25%. Because 𝒮 is independent of the grain alignment efficiency, this demonstrates that the systematic decrease inpwithNHis determined mostly by the magnetic-field structure and not by a drop in grain alignment. This systematic trend is observed both in the diffuse interstellar medium (ISM) and in molecular clouds of the Gould Belt. Second, we look for a dependence of polarization properties on the dust temperature, as we would expect from the radiative alignment torque (RAT) theory. We find no systematic trend of 𝒮 × pwith the dust temperatureTd, whether in the diffuse ISM or in the molecular clouds of the Gould Belt. In the diffuse ISM, lines of sight with high polarization fractionpand low polarization angle dispersion 𝒮 tend, on the contrary, to have colder dust than lines of sight with lowpand high 𝒮. We also compare thePlanckthermal dust polarization with starlight polarization data in the visible at high Galactic latitudes. The agreement in polarization angles is remarkable, and is consistent with what we expect from the noise and the observed dispersion of polarization angles in the visible on the scale of thePlanckbeam. The two polarization emission-to-extinction ratios,RP/pandRS/V, which primarily characterize dust optical properties, have only a weak dependence on the column density, and converge towards the values previously determined for translucent lines of sight. We also determine an upper limit for the polarization fraction in extinction,pV/E(B − V), of 13% at high Galactic latitude, compatible with the polarization fractionp ≈ 20% observed at 353 GHz. Taken together, these results provide strong constraints for models of Galactic dust in diffuse gas.

We report the detection of two radio transients in the Nasu 1.4 GHz wide-field survey. In the survey, we use four pairs of the two-element interferometer aligned east-west to monitor the wide-field sky and simultaneously survey the region at +32° < δ < +42° in drift scanning. In 27 days of continuous observation at a declination between +41° and +42°, we have detected two radio transients of 1 Jy intensity. Since they appeared on only a single day during the 27 days, we consider these detections to be 1 Jy class bursts that brightened and faded within 2 days and have constant emission smaller than 200 mJy, the detection limit of the 27 days of integrated data. While one transient was in low Galactic latitudes, the other transient was detected in high Galactic latitudes and has counterparts only in γ-ray databases. Therefore, the high Galactic latitude transient might be one of the active galactic nuclei that are normally very faint in X-ray and quiet in radio wavelengths.

Studies of the properties of the inner Galactic Bulge depend strongly on the assumptions about the interstellar extinction. Most of the extinction maps available in the literature lack the information about the distance. We combine the observations with the Besancon model of the Galaxy to investigate the variations of extinction along different lines of sight towards the inner Galactic bulge as a function of distance. In addition we study the variations in the extinction law in the Bulge. We construct color-magnitude diagrams with the following sets of colors: H-Ks and J-Ks from the VVV catalogue as well as Ks-[3.6], Ks-[4.5], Ks-[5.8] and Ks-[8.0] from GLIMPSE-II catalogue matched with 2MASS. Using the newly derived temperature-color relation for M giants that match better the observed color-magnitude diagrams we then use the distance-color relations to derive the extinction as a function of distance. The observed colors are shifted to match the intrinsic colors in the Besan\c{c}on model, as a function of distance, iteratively thereby creating an extinction map with three dimensions: two spatial and one distance dimension along each line of sight towards the bulge. Colour excess maps are presented at a resolution of 15' x 15' for 6 different combinations of colors, in distance bins of 1 kpc. The high resolution and depth of the photometry allows to derive extinction maps to 10 kpc distance and up to 35 magnitudes of extinction in Av (3.5 mag in Aks). Integrated maps show the same dust features and consistent values with the other 2D maps. Starting from the color excess between the observations and the model we investigate the extinction law in near-infrared and its variation along different lines of sight.

Large-scale synchrotron loops are recognized as the main source of diffuse radio-continuum emission in the Galaxy at intermediate and high Galactic latitudes, yet their origin remain unexplained. For the first time, using a combination of multi-frequency data in the radio band of total and polarized intensities, we were able to associate one arc, hereafter, the Orion-Taurus ridge, with the wall of the most prominent stellar-feedback blown shell in the Solar neighborhood, namely, the Orion-Eridanus superbubble. We traced the Orion-Taurus ridge using 3D maps of interstellar dust extinction and column-density maps of molecular gas,NH2. We found the Orion-Taurus ridge at a distance of 400 pc, with a plane-of-the-sky extent of 180 pc. Its medianNH2value is (1.4−0.6+2.6) × 1021cm−2. Thanks to the broadband observations below 100 MHz of the Long Wavelength Array (LWA), we were also able to compute the low-frequency spectral-index map of synchrotron emissivity,β, in the Orion-Taurus ridge. We found a flat distribution ofβwith a median value of −2.24−0.02+0.03, which we interpreted in terms of depletion of low-energy (&lt; GeV) cosmic-ray electrons in recent supernova remnants (105–106yr). Our results are consistent with plane-of-the-sky magnetic-field strengths in the Orion-Taurus ridge greater than a few tens of μG (&gt; 30 − 40 μG). We report the first detection of diffuse synchrotron emission from cold-neutral, partly molecular gas in the surroundings of the Orion-Eridanus superbubble. This observation opens a new perspective for studies of the multiphase and magnetized interstellar medium with the advent of future high-sensitivity radio facilities, such as theC-Band All-Sky Survey (C-BASS) and the Square Kilometre Array (SKA).

Extinction in ultraviolet is much more significant than in optical or infrared, which can be very informative to precisely measure the extinction and understand the dust properties in the low-extinction areas. The high Galactic latitude sky is such an area, important for studying the extragalactic sky and the universe. Based on the stellar parameters measured by the LAMOST and GALAH spectroscopy and the ultraviolet photometry by the Galaxy Evolution Explorer (GALEX) space telescope, the extinction of 1,244,504 stars in the GALEX/near-UV band and 56,123 stars in the GALEX/far-UV band is calculated precisely. The error of color excess is 0.009, 0.128, and 0.454 mag for , , and , respectively. They delineate the GALEX/near-UV extinction map of about a third of the sky mainly at the high Galactic latitude area with an angular resolution of ∼ 0.°4. The mean color excess ratio in the entire sky area is derived to be 3.25, 2.95, and −0.37 for , , and , respectively, which is in general agreement with the previous works, and their changes with the Galactic latitude and the interstellar extinction are discussed.

Characterizing the structure of the Galactic interstellar medium (ISM) in three dimensions is of high importance for accurate modeling of dust emission as a foreground to the cosmic microwave background (CMB). At high Galactic latitude, where the total dust content is low, accurate maps of the 3D structure of the ISM are lacking. We develop a method to quantify the complexity of the distribution of dust along the line of sight with the use of H i line emission. The method relies on a Gaussian decomposition of the H i spectra to disentangle the emission from overlapping components in velocity. We use this information to create maps of the number of clouds along the line of sight. We apply the method to (a) the high Galactic latitude sky and (b) the region targeted by the BICEP/Keck experiment. In the north Galactic cap we find on average three clouds per 0.2 square degree pixel, while in the south the number falls to 2.5. The statistics of the number of clouds are affected by intermediate-velocity clouds (IVCs), primarily in the north. IVCs produce detectable features in the dust emission measured by Planck. We investigate the complexity of H i spectra in the BICEP/Keck region and find evidence for the existence of multiple components along the line of sight. The data (doi: 10.7910/DVN/8DA5LH) and software are made publicly available and can be used to inform CMB foreground modeling and 3D dust mapping.

Galactic interstellar extinction maps are powerful and necessary tools for Milky Way structure and stellar population analyses, particularly toward the heavily-reddened bulge and in the midplane. However, due to the difficulty of obtaining reliable extinction measures and distances for a large number of stars that are independent of these maps, tests of their accuracy and systematics have been limited. Our goal is to assess a variety of photometric stellar extinction estimates, including both 2D and 3D extinction maps, using independent extinction measures based on a large spectroscopic sample of stars towards the Milky Way bulge. We employ stellar atmospheric parameters derived from high-resolution $H$-band APOGEE spectra, combined with theoretical stellar isochrones, to calculate line-of-sight extinction and distances for a sample of more than 2400 giants towards the Milky Way bulge. We compare these extinction values to those predicted by individual near-IR and near+mid-IR stellar colors, 2D bulge extinction maps and 3D extinction maps. The long baseline, near+mid-IR stellar colors are, on average, the most accurate predictors of the APOGEE extinction estimates, and the 2D and 3D extinction maps derived from different stellar populations along different sightlines show varying degrees of reliability. We present the results of all of the comparisons and discuss reasons for the observed discrepancies. We also demonstrate how the particular stellar atmospheric models adopted can have a strong impact on this type of analysis, and discuss related caveats.

The contribution from the residuals of the foreground can have a significant impact on the temperature maps of the Cosmic Microwave Background (CMB). Mostly, the focus has been on the galactic plane, when foreground cleaning has taken place. However, in this paper, we will investigate the possible foreground contamination, from sources outside the galactic plane in the CMB maps. We will analyze the correlation between the Faraday rotation map and the CMB temperature map. The Faraday rotation map is dependent on the galactic magnetic field, as well as the thermal electron density, and both may contribute to the CMB temperature. We find that the standard deviation for the mean cross correlation deviate from that of simulations at the 99.9% level. Additionally, a comparison between the CMB temperature extrema and the extremum points of the Faraday rotation is also performed, showing a general overlap between the two. Also we find that the CMB Cold Spot is located at an area of strong negative cross correlation, meaning that it may be explained by a galactic origin. Further, we investigate nearby supernova remnants in the galaxy, traced by the galactic radio loops. These super nova remnants are located at high and low galactic latitude, and thus well outside the galactic plane. We find some correlation between the Faraday Rotation and the CMB temperature, at select radio loops. This indicate, that the galactic foregrounds may affect the CMB, at high galactic latitudes

The Milky Way dust extinction curve in the near-infrared (NIR) follows a power-law form, but the value of the slope, β NIR, is debated. Systematic variations in the slope of the Milky Way UV extinction curve are known to be correlated with variations in the optical slope (through R V ), but whether such a dependence extends to the NIR is unclear. Finally, because of low dust column densities, the NIR extinction law is poorly understood at high Galactic latitudes where most extragalactic work takes place. In this paper, we construct extinction curves from 56,649 stars with Sloan Digital Sky Survey (SDSS) and Two Micron All Sky Survey photometry, based on stellar parameters from SDSS spectra. We use dust maps to identify dust-free stars, from which we calibrate the relation between stellar parameters and intrinsic colors. Furthermore, to probe the low-dust regime at high latitudes, we use aggregate curves based on many stars. We find no systematic variation of β NIR across low-to-moderate dust columns (0.02 < E(B − V) ≲ 1), and report average β NIR = 1.85 ± 0.01, in agreement with the law in the 2019 Fitzpatrick et al. study, but steeper than the Cardelli et al. and 1999 Fitzpatrick laws. Star-to-star scatter in β NIR is relatively small (σ(β NIR) = 0.13). We also find no intrinsic correlation between β NIR and R V (there is an apparent correlation that is the result of the correlated uncertainties in the two values). These results hold for typical sightlines; we do not probe very dusty regions near the Galactic Center, nor rare sightlines with R V > 4. Finally, we find R H = 0.345 ± 0.007 and comment on its bearing on Cepheid calibrations and the determination of H 0.

The product of the previously constructed 3D maps of stellar reddening (Gontcharov 2010) and $R_V$ variations (Gontcharov 2012) has allowed us to produce a 3D interstellar extinction map within the nearest kiloparsec from the Sun with a spatial resolution of 50 pc and an accuracy of $0.2^m$. This map is compared with the 2D reddening map by Schlegel et al. (1998), the 3D extinction map at high latitudes by Jones et al. (2011), and the analytical extinction models by Arenou et al. (1992) and Gontcharov (2009). In all cases, we have found good agreement and show that there are no systematic errors in the new map everywhere except the direction toward the Galactic center. We have found that the map by Schlegel et al. (1998) reaches saturation near the Galactic equator at $E_{(B-V)}>0.8^m$, has a zero-point error and systematic errors gradually increasing with reddening, and among the analytical models those that take into account the extinction in the Gould Belt are more accurate. Our extinction map shows that it is determined by reddening variations at low latitudes and $R_V$ variations at high ones. This naturally explains the contradictory data on the correlation or anticorrelation between reddening and $R_V$ available in the literature. There is a correlation in a thin layer near the Galactic equator, because both reddening and $R_V$ here increase toward the Galactic center. There is an anticorrelation outside this layer, because higher values of $R_V$ correspond to lower reddening at high and middle latitudes. Systematic differences in sizes and other properties of the dust grains in different parts of the Galaxy manifest themselves in this way. The largest structures within the nearest kiloparsec, including the Local Bubble, the Gould Belt, the Great Tunnel, the Scorpius, Perseus, Orion, and other complexes, have manifested themselves in the constructed map.

We use SDSS photometry of 73 million stars to simultaneously obtain best-fit main-sequence stellar energy distribution (SED) and amount of dust extinction along the line of sight towards each star. Using a subsample of 23 million stars with 2MASS photometry, whose addition enables more robust results, we show that SDSS photometry alone is sufficient to break degeneracies between intrinsic stellar color and dust amount when the shape of extinction curve is fixed. When using both SDSS and 2MASS photometry, the ratio of the total to selective absorption, $R_V$, can be determined with an uncertainty of about 0.1 for most stars in high-extinction regions. These fits enable detailed studies of the dust properties and its spatial distribution, and of the stellar spatial distribution at low Galactic latitudes. Our results are in good agreement with the extinction normalization given by the Schlegel et al. (1998, SFD) dust maps at high northern Galactic latitudes, but indicate that the SFD extinction map appears to be consistently overestimated by about 20% in the southern sky, in agreement with Schlafly et al. (2010). The constraints on the shape of the dust extinction curve across the SDSS and 2MASS bandpasses support the models by Fitzpatrick (1999) and Cardelli et al. (1989). For the latter, we find an $R_V=3.0\pm0.1$(random) $\pm0.1$(systematic) over most of the high-latitude sky. At low Galactic latitudes (|b|<5), we demonstrate that the SFD map cannot be reliably used to correct for extinction as most stars are embedded in dust, rather than behind it. We introduce a method for efficient selection of candidate red giant stars in the disk, dubbed "dusty parallax relation", which utilizes a correlation between distance and the extinction along the line of sight. We make these best-fit parameters, as well as all the input SDSS and 2MASS data, publicly available in a user-friendly format.