ABSTRACT
We present the first simulations evolving resolved spectra of cosmic rays (CRs) from MeV–TeV energies (including electrons, positrons, (anti)protons, and heavier nuclei), in live ...kinetic-magnetohydrodynamics galaxy simulations with star formation and feedback. We utilize new numerical methods including terms often neglected in historical models, comparing Milky Way analogues with phenomenological scattering coefficients ν to Solar-neighbourhood Local interstellar medium (LISM) observations (spectra, B/C, e+/e−, $\mathrm{\bar{p}}/\mathrm{p}$, 10Be/9Be, ionization, and γ-rays). We show it is possible to reproduce observations with simple single-power-law injection and scattering coefficients (scaling with rigidity R), similar to previous (non-dynamical) calculations. We also find: (1) The circumgalactic medium in realistic galaxies necessarily imposes an $\sim 10\,$ kpc CR scattering halo, influencing the required ν(R). (2) Increasing the normalization of ν(R) re-normalizes CR secondary spectra but also changes primary spectral slopes, owing to source distribution and loss effects. (3) Diffusive/turbulent reacceleration is unimportant and generally sub-dominant to gyroresonant/streaming losses, which are sub-dominant to adiabatic/convective terms dominated by $\sim 0.1-1\,$ kpc turbulent/fountain motions. (4) CR spectra vary considerably across galaxies; certain features can arise from local structure rather than transport physics. (5) Systematic variation in CR ionization rates between LISM and molecular clouds (or Galactic position) arises naturally without invoking alternative sources. (6) Abundances of CNO nuclei require most CR acceleration occurs around when reverse shocks form in SNe, not in OB wind bubbles or later Sedov–Taylor stages of SNe remnants.
Abstract
Interstellar dust at high Galactic latitudes can influence astronomical foreground subtraction, produce diffuse scattered light, and soften the UV spectra of quasars. In a sample of 94 sight ...lines toward quasars at high latitude and low extinction, we evaluate the interstellar “gas-to-dust ratio”
N
H
/
E
(
B
−
V
), using hydrogen column densities (H
i
and H
2
) and far-IR (FIR) estimates of dust reddening. In the Galactic plane, this ratio is 6.0 ± 0.2 (in units of 10
21
cm
−2
mag
−1
). On average, recent Planck estimates of
E
(
B
−
V
) in low-reddening sight lines are 12% higher than those from Schlafly & Finkbeiner, and
N
H I
exhibits significant variations when measured at different radio telescopes. In a sample of 51 quasars with measurements of both H
i
and H
2
and 0.01 ≤
E
(
B
−
V
) ≲ 0.1, we find mean ratios 10.3 ± 0.4 (gas at all velocities) and 9.2 ± 0.3 (low-velocity only) using Planck
E
(
B
−
V
) data. High-latitude H
2
fractions are generally small (2%–3% on average), although nine of 39 sight lines at ∣
b
∣ ≥ 40° have
f
H2
of 1%–17%. Because FIR-inferred
E
(
B
−
V
) is sensitive to modeled dust temperature
T
d
and emissivity index
β
, gas-to-dust ratios have large, asymmetric errors at low
E
(
B
−
V
). The ratios are elevated in sight lines with high-velocity clouds, which contribute
N
H
but little reddening. In Complex C, the ratio decreases by 40% when high-velocity gas is excluded. Decreases in dust content are expected in low-metallicity gas above the Galactic plane, resulting from grain destruction in shocks, settling to the disk, and thermal sputtering in hot halo gas.
Abstract
We present optical and infrared imaging and spectroscopy of the R Coronae Borealis-type (R Cor Bor) star IRAS 00450+7401. Optical spectra further confirm its classification as a cool R Cor ...Bor system, having a hydrogen-deficient carbon star spectral subclass of HdC5 or later. Mid-infrared spectroscopy reveals the typical ∼8
μ
m “hump” seen in other R Cor Bor stars and no other features. A modern-epoch spectral energy distribution shows bright emission from hot dust having
T
dust
> 600 K. Historical infrared data reveal generally cooler dust color temperatures combined with long-term fading trends, but provide no discernible correlation between flux level and temperature. Investigating the most mid-infrared variable R Cor Bor stars found in IRAS, AKARI, and WISE data reveals similar fading trends, bursts that can show a factor of up to 10 change in flux density between epochs, and blackbody-fit dust color temperatures that span 400–1300 K. While some R Cor Bor stars such as IRAS 00450+7401 appear to undergo fade/burst cycles in the mid-infrared, significant gaps in temporal coverage prevent conclusively identifying any preferred timescale for their mid-infrared variability and circumstellar dust temperature changes.
ABSTRACT
Characterizing the sub-mm Galactic emission has become increasingly critical especially in identifying and removing its polarized contribution from the one emitted by the cosmic microwave ...background (CMB). In this work, we present a parametric foreground removal performed on to sub-patches identified in the celestial sphere by means of spectral clustering. Our approach takes into account efficiently both the geometrical affinity and the similarity induced by the measurements and the accompanying errors. The optimal partition is then used to parametrically separate the Galactic emission encoding thermal dust and synchrotron from the CMB one applied on two nominal observations of forthcoming experiments from the ground and from the space. Moreover, the clustering is performed on tracers that are different from the data used for component separation, e.g. the spectral index maps of dust and synchrotron. Performing the parametric fit singularly on each of the clustering derived regions results in an overall improvement: both controlling the bias and the uncertainties in the CMB B-mode recovered maps. We finally apply this technique using the map of the number of clouds along the line of sight, $\mathcal {N}_c$, as estimated from H i emission data and perform parametric fitting on to patches derived by clustering on this map. We show that adopting the $\mathcal {N}_c$ map as a tracer for the patches related to the thermal dust emission, results in reducing the B-mode residuals post-component separation. The code is made publicly available https://github.com/giuspugl/fgcluster.
ABSTRACT
The physics of magnetic fields (B) and cosmic rays (CRs) have recently been included in simulations of galaxy formation. However, significant uncertainties remain in how these components ...affect galaxy evolution. To understand their common observational tracers, we analyse the magnetic fields in a set of high-resolution, magnetohydrodynamic, cosmological simulations of Milky-Way-like galaxies from the FIRE-2 project. We compare mock observables of magnetic field tracers for simulations with and without CRs to observations of Zeeman splitting and rotation/dispersion measures. We find reasonable agreement between simulations and observations in both the neutral and the ionized interstellar medium (ISM). We find that the simulated galaxies with CRs show weaker ISM |B| fields on average compared to their magnetic-field-only counterparts. This is a manifestation of the effects of CRs in the diffuse, low density inner circumgalactic medium (CGM). We find that equipartition between magnetic and cosmic ray energy densities may be valid at large (> 1 kpc) scales for typical ISM densities of Milky-Way-like galaxies, but not in their haloes. Within the ISM, the magnetic fields in our simulated galaxies follow a power-law scaling with gas density. The scaling extends down to neutral hydrogen number densities < 300 cm−3, in contrast to observationally derived models, but consistent with the observational measurements. Finally, we generate synthetic rotation measure (RM) profiles for projections of the simulated galaxies and compare to observational constraints in the CGM. While consistent with upper limits, improved data are needed to detect the predicted CGM RMs at 10–200 kpc and better constrain theoretical predictions.
If a single line of sight (LOS) intercepts multiple dust clouds with different spectral energy distributions and magnetic field orientations, then the frequency scaling of each of the Stokes
Q
and
U
...parameters of the thermal dust emission may be different, a phenomenon we refer to as LOS frequency decorrelation. We present first evidence for LOS frequency decorrelation in
Planck
data using independent measurements of neutral-hydrogen (H
I
) emission to probe the 3D structure of the magnetized interstellar medium (ISM). We use H
I
-based measurements of the number of clouds per LOS and the magnetic field orientation in each cloud to select two sets of sightlines: (i) a target sample of pixels that are likely to exhibit LOS frequency decorrelation and (ii) a control sample of pixels that lack complex LOS structure. We test the null hypothesis that LOS frequency decorrelation is not detectable in
Planck
353 and 217 GHz polarization data at high Galactic latitudes. We reject the null hypothesis at high significance based on data that show that the combined effect of polarization angle variation with frequency and depolarization are detected in the target sample. This detection is robust against the choice of cosmic microwave background (CMB) map and map-making pipeline. The observed change in polarization angle due to LOS frequency decorrelation is detectable above the
Planck
noise level. The probability that the detected effect is due to noise alone ranges from 5 × 10
−2
to 4 × 10
−7
, depending on the CMB subtraction algorithm and treatment of residual systematic errors; correcting for residual systematic errors consistently increases the significance of the effect. Within the target sample, the LOS decorrelation effect is stronger for sightlines with more misaligned magnetic fields, as expected. With our sample, we estimate that an intrinsic variation of ~15% in the ratio of 353 to 217 GHz polarized emission between clouds is sufficient to reproduce the measured effect. Our finding underlines the importance of ongoing studies to map the three-dimensional structure of the magnetized and dusty ISM that could ultimately help component separation methods to account for frequency decorrelation effects in CMB polarization studies.
The availability of large data sets with stellar distance and polarization information will enable a tomographic reconstruction of the (plane-of-the-sky-projected) interstellar magnetic field in the ...near future. We demonstrate the feasibility of such a decomposition within a small region of the diffuse interstellar medium (ISM). We combine measurements of starlight (R-band) linear polarization obtained using the RoboPol polarimeter with stellar distances from the second Gaia data release. The stellar sample is brighter than 17 mag in the R-band and reaches out to several kiloparsecs from the Sun. H i emission spectra reveal the existence of two distinct clouds along the line of sight. We decompose the line-of-sight-integrated stellar polarizations to obtain the mean polarization properties of the two clouds. The two clouds exhibit significant differences in terms of column density and polarization properties. Their mean plane-of-the-sky magnetic field orientation differs by 60°. We show how our tomographic decomposition can be used to constrain our estimates of the polarizing efficiency of the clouds as well as the frequency dependence of the polarization angle of polarized dust emission. We also demonstrate a new method to constrain cloud distances based on this decomposition. Our results represent a preview of the wealth of information that can be obtained from a tomographic map of the ISM magnetic field.
ABSTRACT
Radiation-dust driven outflows, where radiation pressure on dust grains accelerates gas, occur in many astrophysical environments. Almost all previous numerical studies of these systems have ...assumed that the dust was perfectly coupled to the gas. However, it has recently been shown that the dust in these systems is unstable to a large class of ‘resonant drag instabilities’ (RDIs) which de-couple the dust and gas dynamics and could qualitatively change the non-linear outcome of these outflows. We present the first simulations of radiation-dust driven outflows in stratified, inhomogeneous media, including explicit grain dynamics and a realistic spectrum of grain sizes and charge, magnetic fields and Lorentz forces on grains (which dramatically enhance the RDIs), Coulomb and Epstein drag forces, and explicit radiation transport allowing for different grain absorption and scattering properties. In this paper, we consider conditions resembling giant molecular clouds (GMCs), H ii regions, and distributed starbursts, where optical depths are modest (≲1), single-scattering effects dominate radiation-dust coupling, Lorentz forces dominate over drag on grains, and the fastest-growing RDIs are similar, such as magnetosonic and fast-gyro RDIs. These RDIs generically produce strong size-dependent dust clustering, growing non-linear on time-scales that are much shorter than the characteristic times of the outflow. The instabilities produce filamentary and plume-like or ‘horsehead’ nebular morphologies that are remarkably similar to observed dust structures in GMCs and H ii regions. Additionally, in some cases they strongly alter the magnetic field structure and topology relative to filaments. Despite driving strong micro-scale dust clumping which leaves some gas ‘behind,’ an order-unity fraction of the gas is always efficiently entrained by dust.
ABSTRACT
Synchrotron emission is one of few observable tracers of galactic magnetic fields (B) and cosmic rays (CRs). Much of our understanding of B in galaxies comes from utilizing synchrotron ...observations in conjunction with several simplifying assumptions of equipartition models, however, it remains unclear how well these assumptions hold, and what B these estimates physically represent. Using Feedback in Realistic Environments project simulations which self-consistently evolve CR proton, electron, and positron spectra from MeV to TeV energies, we present the first synthetic synchrotron emission predictions from simulated L* galaxies with ‘live’ spectrally resolved CR-magnetohydrodynamic. We find that synchrotron emission can be dominated by relatively cool and dense gas, resulting in equipartition estimates of B with fiducial assumptions underestimating the ‘true’ B in the gas that contributes the most emission by factors of 2–3 due to small volume-filling factors. Motivated by our results, we present an analytical framework that expands upon equipartition models for estimating B in a multiphase medium. Comparing our spectrally resolved synchrotron predictions to simpler spectral assumptions used in galaxy simulations with CRs, we find that spectral evolution can be crucial for accurate synchrotron calculations towards galactic centres, where loss terms are large.
Context. Galactic dust emission is polarized at unexpectedly high levels, as revealed by Planck. Aims. The origin of the observed ≃20% polarization fractions can be identified by characterizing the ...properties of optical starlight polarization in a region with maximally polarized dust emission. Methods. We measure the R-band linear polarization of 22 stars in a region with a submillimeter polarization fraction of ≃20%. A subset of 6 stars is also measured in the B, V, and I bands to investigate the wavelength dependence of polarization. Results. We find that starlight is polarized at correspondingly high levels. Through multiband polarimetry we find that the high polarization fractions are unlikely to arise from unusual dust properties, such as enhanced grain alignment. Instead, a favorable magnetic field geometry is the most likely explanation, and is supported by observational probes of the magnetic field morphology. The observed starlight polarization exceeds the classical upper limit of pV/E(B−V)max = 9% mag−1 and is at least as high as 13% mag−1, as inferred from a joint analysis of Planck data, starlight polarization, and reddening measurements. Thus, we confirm that the intrinsic polarizing ability of dust grains at optical wavelengths has long been underestimated.