Magnetized Interstellar Medium Research Articles

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 H I spectra 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 PASIPHAE survey, 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.

Energetic particles in the central starburst, disc, and halo of NGC253

ABSTRACT Detailed modelling of the spectro-spatial distributions of energetic electrons and protons in galactic discs and haloes of starburst galaxies (SBGs) is needed in order to follow their interactions with the magnetized interstellar medium and radiation fields, determine their radiative yields, and for estimating their residual spectral densities in intergalactic environments. We have developed a semi-analytical approach for calculating the particle spectro-spatial distributions in the disc and halo based on a diffusion model for particle propagation from acceleration sites in the central SB and disc regions, including all their relevant interaction modes. Important overall normalization of our models is based on previous modelling of the Galactic disc (with the galprop code), scaled to the higher star-formations rate in NGC253, and on spatially resolved radio measurements of the central SB and disc. These provide the essential input for determining the particle distributions and their predicted radiative yields in the outer disc and inner halo for a range of values of the key parameters that affect diffusion rate and energy losses. Results of our work clearly indicate that quantitative description of non-thermal emission in SBGs has to be based on modelling of the particle distributions in the entire disc, not just the central SB region.

A new star formation recipe for magnetohydrodynamics simulations of galaxy formation

ABSTRACT Star formation has been observed to occur at globally low yet locally varying efficiencies. As such, accurate capture of star formation in numerical simulations requires mechanisms that can replicate both its smaller scale variations and larger scale properties. Magnetic fields are thought to play an essential role within the turbulent interstellar medium (ISM) and affect molecular cloud collapse. However, it remains to be fully explored how a magnetized model of star formation might influence galaxy evolution. We present a new model for a sub-grid star formation recipe that depends on the magnetic field. We run isolated disc galaxy simulations to assess its impact on the regulation of star formation using the code ramses. Building upon existing numerical methods, our model derives the star formation efficiency from local properties of the sub-grid magnetized ISM turbulence, assuming a constant Alfvén speed at sub-parsec scales. Compared to its non-magnetized counterpart, our star formation model suppresses the initial starburst by a factor of 2 while regulating star formation later on to a nearly constant rate of ∼1 M⊙ yr−1. Differences also arise in the local Schmidt law with a shallower power-law index for the magnetized star formation model. Our results encourage further examination into the notion that magnetic fields are likely to play a non-trivial role in our understanding of star and galaxy formation.

Frequency dependence of the thermal dust E/B ratio and EB correlation: Insights from the spin-moment expansion

The change of physical conditions across the turbulent and magnetized interstellar medium induces a 3D spatial variation of the properties of Galactic polarized emission. The observed signal results from the averaging of different spectral energy distributions (SEDs) and polarization angles along and between lines of sight. As a consequence, the total Stokes parameters Q and U will have different frequency dependencies, both departing from the canonical emission law, so that the polarization angle becomes frequency dependent. In the present work, we show how this phenomenon similarly induces a different, distorted SED for the three polarized angular power spectra 𝒟𝓁EE, 𝒟𝓁BB, and 𝒟𝓁EB, implying a variation of the 𝒟𝓁EE/𝒟𝓁BB ratio with frequency. We demonstrate how the previously introduced “spin-moment” formalism provides a natural framework to grasp these effects and enables us to derive analytical predictions for the spectral behaviors of the polarized spectra, focusing here on the example of thermal dust polarized emission. After a quantitative discussion based on a model combining emission from a filament with its background, we further reveal that the spectral complexity implemented in the dust models commonly used by the cosmic microwave background (CMB) community includes different distortions for the three polarized power-spectra. This new understanding is crucial for CMB component separation, in which extreme accuracy is required for the modeling of the dust signal to allow for the search of the primordial imprints of inflation or cosmic birefringence. For the latter, as long as the dust EB signal is not measured accurately, great caution is required regarding the assumptions made to model its spectral behavior, as it may not be inferred from the other dust angular power spectra.

Mixing of materials in magnetized core-collapse supernova remnants

ABSTRACTCore-collapse supernova remnants are structures of the interstellar medium (ISM) left behind the explosive death of most massive stars ($\lesssim\!40\, \rm {\rm M}_{\odot }$). Since they result in the expansion of the supernova shock wave into the gaseous environment shaped by the star’s wind history, their morphology constitutes an insight into the past evolution of their progenitor star. Particularly, fast-moving massive stars can produce asymmetric core-collapse supernova remnants. We investigate the mixing of materials in core-collapse supernova remnants generated by a moving massive $35\, \rm {\rm M}_{\odot }$ star, in a magnetized ISM. Stellar rotation and the wind magnetic field are time-dependently included into the models which follow the entire evolution of the stellar surroundings from the zero-age main-sequence to $80\, \rm kyr$ after the supernova explosion. It is found that very little main-sequence material is present in remnants from moving stars, that the Wolf-Rayet wind mixes very efficiently within the $10\, \rm kyr$ after the explosion, while the red supergiant material is still unmixed by 30 per cent within $50\, \rm kyr$ after the supernova. Our results indicate that the faster the stellar motion, the more complex the internal organization of the supernova remnant and the more effective the mixing of ejecta therein. In contrast, the mixing of stellar wind material is only weakly affected by progenitor motion, if at all.

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 PASIPHAE 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.

Ion alfvén velocity fluctuations and implications for the diffusion of streaming cosmic rays

The interstellar medium (ISM) of star-forming galaxies is magnetized and turbulent. Cosmic rays (CRs) propagate through it, and those with energies from ∼ GeV − TeV are likely subject to the streaming instability, whereby the wave damping processes balances excitation of resonant ionic Alfvén waves by the CRs, reaching an equilibrium in which the propagation speed of the CRs is very close to the local ion Alfvén velocity. The transport of streaming CRs is therefore sensitive to ionic Alfvén velocity fluctuations. In this paper we systematically study these fluctuations using a large ensemble of compressible MHD turbulence simulations. We show that for sub-Alfvénic turbulence, as applies for a strongly magnetized ISM, the ionic Alfvén velocity probability density function (PDF) is determined solely by the density fluctuations from shocked gas forming parallel to the magnetic field, and we develop analytical models for the ionic Alfvén velocity PDF up to second moments. For super-Alfvénic turbulence, magnetic and density fluctuations are correlated in complex ways, and these correlations as well as contributions from the magnetic fluctuations sets the ionic Alfvén velocity PDF. We discuss the implications of these findings for underlying “macroscopic” diffusion mechanisms in CRs undergoing the streaming instability, including modeling the macroscopic diffusion coefficient for the parallel transport in sub-Alfvénic plasmas. We also describe how, for highly-magnetized turbulent gas, the gas density PDF, and hence column density PDF, can be used to access information about ionic Alfvén velocity structure from observations of the magnetized ISM.

Shock-induced transformation of non-magnetic to magnetic ISM dust analogue

ABSTRACT The fate of organometallic dust subjected to extreme conditions, especially the shock fronts, in the interstellar medium (ISM) has not been explored to date. Iron and cyclopentadiene, which have been found to be present in the ISM, are known to react and produce the organometallic compound ferrocene under terrestrial conditions. In our experiment we subjected ferrocene, a possible proxy of such dust, to ∼5.6 Mach (M) shock commensurate with conditions encountered in the ISM, leading to a temperature rise up to 7300 K within 2 ms. Analysis of the post-shock residue showed the presence of a α-Fe and Fe3C composite that responded to an external magnetic field. These results show that a non-magnetic dust composed of molecules containing transition metals undergoing shock processing in the ISM can dissociate and synthesize dust that is then magnetic. Such drastic transformations from non-magnetic to magnetic dust induced by shocks might be of importance in interstellar polarization.

Energy balance and Alfvén Mach numbers in compressible magnetohydrodynamic turbulence with a large-scale magnetic field

ABSTRACT Energy equipartition is a powerful theoretical tool for understanding astrophysical plasmas. It is invoked, for example, to measure magnetic fields in the interstellar medium (ISM), as evidence for small-scale turbulent dynamo action, and, in general, to estimate the energy budget of star-forming molecular clouds. In this study, we motivate and explore the role of the volume-averaged root-mean-squared (rms) magnetic coupling term between the turbulent, $\delta {\boldsymbol{B}}$ , and large-scale, ${\boldsymbol{B}}_0$, fields, ${\left\langle (\delta \mathrm{{\boldsymbol {\mathit {B}}}}\cdot {\mathrm{{\boldsymbol {\mathit {B}}}}_0})^{2} \right\rangle ^{1/2}_{\mathcal {V}}}$. By considering the second moments of the energy balance equations we show that the rms coupling term is in energy equipartition with the volume-averaged turbulent kinetic energy for turbulence with a sub-Alfvénic large-scale field. Under the assumption of exact energy equipartition between these terms, we derive relations for the magnetic and coupling term fluctuations, which provide excellent, parameter-free agreement with time-averaged data from 280 numerical simulations of compressible magnetohydrodynamic (MHD) turbulence. Furthermore, we explore the relation between the turbulent mean field and total Alfvén Mach numbers, and demonstrate that sub-Alfvénic turbulence can only be developed through a strong, large-scale magnetic field, which supports an extremely super-Alfvénic turbulent magnetic field. This means that the magnetic field fluctuations are significantly subdominant to the velocity fluctuations in the sub-Alfvénic large-scale field regime. Throughout our study, we broadly discuss the implications for observations of magnetic fields and understanding the dynamics in the magnetized ISM.

First look at the multiphase interstellar medium using synthetic observations from low-frequency Faraday tomography

Faraday tomography of radio polarimetric data below 200 MHz from the LOw Frequency ARray (LOFAR) has been providing new perspectives on the diffuse and magnetized interstellar medium (ISM). One aspect of particular interest is the unexpected discovery of Faraday-rotated synchrotron polarization associated with structures of neutral gas, as traced by atomic hydrogen (HI) and dust. Here, we present the first in-depth numerical study of these LOFAR results. We produced and analyzed comprehensive synthetic observations of low-frequency synchrotron polarization from magneto-hydrodynamical (MHD) simulations of colliding super shells in the multiphase ISM from the literature. Using an analytical approach to derive the ionization state of the multiphase gas, we defined five distinct gas phases over more than four orders of magnitude in gas temperature and density, ranging from hot, and warm fully ionized gas to a cold neutral medium. We focused on establishing the contribution of each gas phase to synthetic observations of both rotation measure and synchrotron polarized intensity below 200 MHz. We also investigated the link between the latter and synthetic observations of optically thin HI gas. We find that it is not only the fully ionized gas, but also the warm partially ionized and neutral phases that strongly contribute to the total rotation measure and polarized intensity. However, the contribution of each phase to the observables strongly depends on the choice of the integration axis and the orientation of the mean magnetic field with respect to the shell collision axis. A strong correlation between the HI synthetic data and synchrotron polarized intensity, reminiscent of the LOFAR results, is obtained with lines of sight perpendicular to the mean magnetic field direction. Our study suggests that multiphase modeling of MHD processes is needed in order to interpret observations of the radio sky at low frequencies. This work is a first step toward understanding the complexity of low-frequency synchrotron emission that will be soon revolutionized thanks to large-scale surveys with LOFAR and the Square Kilometre Array.

Polarization power spectra and dust cloud morphology

Context In the framework of studying cosmic microwave background polarization and characterizing its Galactic foregrounds, the angular power spectrum analysis of the thermal dust polarization map has led to intriguing evidence of an E/B asymmetry and a positive TE correlation. The interpretation of these observations is the subject of theoretical and simulation-driven studies in which the correlation between the density structure of the interstellar medium (ISM) and the magnetic field appears to be a key aspect. In this context, and when the magnetized ISM structures are modeled in three dimensions, dust clouds are generally considered to be filamentary structures only, but both filamentary and sheet-like shapes are supported by observational and theoretical evidence. Aims. We aim to study the influence of the cloud shape and its connection to the local magnetic field, as well as the influence from the viewing angle, on the angular power spectra measured on thermal dust polarization maps; we specifically focus on the dependence of the E/B power asymmetry and TE correlation. Methods. To this end, we simulated realistic interstellar clouds with both filament-like and sheet-like shapes using the software Asterion, which also allowed us to generate synthetic maps of thermal dust polarized emission with an area of 400 square degrees. Then, we computed their polarization power spectra in the multipole range ℓ ϵ [100, 500] and focused on the E/B power asymmetry, quantified through the ℛEB ratio, and the correlation coefficient rTE between Τ and Ε modes. We quantified the dependence of ℛEB and rTE values on the offset angle (between the longest cloud axis and local magnetic field lines) and inclination angle (between the line of sight and the magnetic field) for both types of cloud shapes, either embedded in a regular magnetic field or coupled to a nonregular field to mimic turbulence. Results. We find that both types of cloud shapes cover the same regions of the (ℛEB, rTE) parameter space. The dependence on the inclination and offset angles is similar for both shapes, although sheet-like structures generally show larger scatter than filamentary structures. In addition to the known dependence on the offset angle, we find a strong dependence of ℛEB and rTE on the inclination angle. Conclusions. The very fact that filament-like and sheet-like structures may lead to polarization power spectra with similar (ℛEB,rTE) values complicates their interpretation. We argue that interpreting them solely in terms of filament characteristics is risky, and in future analyses, this degeneracy should be accounted for, as should the connection to the magnetic field geometry. Our results based on maps of 400 square degrees clarify that the overall geometrical arrangement of the magnetized ISM surrounding the observer leaves its marks on polarization power spectra.

A Comparison of Multiphase Magnetic Field Tracers in a High Galactic Latitude Region of the Filamentary Interstellar Medium

Abstract Understanding how the Galactic magnetic field threads the multiphase interstellar medium (ISM) remains a considerable challenge, as different magnetic field tracers probe dissimilar phases and field components. We search for evidence of a common magnetic field shared between the ionized and neutral ISM by comparing 1.4 GHz radio continuum polarization and H i line emission from the Galactic Arecibo L-Band Feed Array Continuum Transit Survey (GALFACTS) and Galactic Arecibo L-Band Feed Array H i (GALFA-H i) survey, respectively. We compute the polarization gradient of the continuum emission and search for associations with diffuse/translucent H i structures. The polarization gradient is sensitive to changes in the integrated product of the thermal electron density and line-of-sight field strength (B ∥) in warm ionized gas, while narrow H i structures highlight the plane-of-sky field orientation in cold neutral gas. We identified one region in the high Galactic latitude Arecibo sky, G216+26 centered on (ℓ, b) ∼ (216°, +26°), containing filaments in the polarization gradient that are aligned with narrow H i structures roughly parallel to the Galactic plane. We present a comparison of multiphase observations and magnetic field tracers of this region, demonstrating that the warm ionized and cold neutral media are connected likely via a common magnetic field. We quantify the physical properties of a polarization gradient filament associated with Hα emission, measuring a line-of-sight field strength B ∥ = 6 ± 4 μG and a plasma beta β = 2.1 − 2.1 + 3.1 . We discuss the lack of widespread multiphase magnetic field alignments and consider whether this region is associated with a short-timescale or physically rare phenomenon. This work highlights the utility of multitracer analyses for understanding the magnetized ISM.

A new approach for the statistical denoising of Planck interstellar dust polarization data

Dust emission is the main foreground for cosmic microwave background polarization. Its statistical characterization must be derived from the analysis of observational data because the precision required for a reliable component separation is far greater than what is currently achievable with physical models of the turbulent magnetized interstellar medium. This Letter takes a significant step toward this goal by proposing a method that retrieves non-Gaussian statistical characteristics of dust emission from noisy Planck polarization observations at 353 GHz. We devised a statistical denoising method based on wavelet phase harmonics (WPH) statistics, which characterize the coherent structures in non-Gaussian random fields and define a generative model of the data. The method was validated on mock data combining a dust map from a magnetohydrodynamic simulation and Planck noise maps. The denoised map reproduces the true power spectrum down to scales where the noise power is an order of magnitude larger than that of the signal. It remains highly correlated to the true emission and retrieves some of its non-Gaussian properties. Applied to Planck data, the method provides a new approach to building a generative model of dust polarization that will characterize the full complexity of the dust emission. We also release PyWPH, a public Python package, to perform GPU-accelerated WPH analyses on images.

Evidence for line-of-sight frequency decorrelation of polarized dust emission in Planck data

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 (HI) emission to probe the 3D structure of the magnetized interstellar medium (ISM). We use HI-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 multiphase and magnetized neutral hydrogen seen by LOFAR

Faraday tomography of polarimetric observations at low frequency in the radio is a unique tool for studying the structure of the magneto-ionic diffuse interstellar medium (ISM) based on Faraday depth. LOFAR data below 200 MHz have revealed a plethora of features in polarization, whose origin remains unknown. Previous studies have highlighted the remarkable association of such features with tracers of the magnetized-neutral ISM, such as interstellar dust and atomic hydrogen (HI). However, the physical conditions responsible for the correlation between magneto-ionic and neutral media have not been clarified yet. In this Letter we further investigate the correlation between LOFAR data and the HI spectroscopic observations at 21 cm from the Effelsberg-Bonn HI Survey (EBHIS). We focus on the multiphase properties of the HI gas. We present the first statistical study on the morphological correlation between LOFAR tomographic data and the cold (CNM), lukewarm (LNM), and warm (WNM) neutral medium HI phases. We use the Regularized Optimization for Hyper-Spectral Analysis approach to decompose the HI phases based on a Gaussian decomposition of the HI spectra. We study four fields of view – Fields 3C196, A, B, and C – and find, in at least the first two, a significant correlation between the LOFAR and EBHIS data using the histograms of oriented gradients (HOG) feature. The absence of a correlation in Fields B and C is caused by a low signal-to-noise ratio in polarization. The observed HOG correlation in Fields 3C196 and A is associated with all HI phases and it is surprisingly dominant in the CNM and LNM phases. We discuss possible mechanisms that would explain the correlation between CNM, LNM, and WNM with polarized emission at Faraday depths up to 10 rad m−2. Our results show how the complex structure of the ionic medium seen by the LOFAR tomographic data is tightly related to phase transition in the diffuse and magnetized neutral ISM traced by HI spectroscopic data.

Two Modes of Carbonaceous Dust Alignment

Abstract Radiative torques (RATs) or mechanical torques acting on irregular grains can induce the alignment of dust grains in respect to the alignment axis (AA), which can be either the direction of the magnetic field or the direction of the radiation. We show that carbonaceous grains can be aligned with their axes both parallel and perpendicular to the AA, and we explore the conditions where the particular mode of alignment takes place. We identify a new process of alignment of charged carbonaceous grains in a turbulent, magnetized interstellar medium with respect to an electric field. This field acts on grains accelerated in a turbulent medium and gyrorotating about a magnetic field. The electric field can also arise from the temporal variations of the magnetic field strength in turbulent, compressible media. The direction of the electric field is perpendicular to the magnetic field, and the carbonaceous grains precess in the electric field because of their electric moments. If this precession is faster than Larmor precession in the magnetic field, the alignment of such grains is with their long axes parallel to the magnetic field. We explore the parameter space for which the new mechanism aligns grains with long axes parallel to the magnetic field. We compare this mechanism with another process that provides the same type of alignment, namely, the RAT alignment of grains with insufficiently fast internal relaxation. We describe the conditions for which the particular mode of carbonaceous grain alignment is realized and discuss what information can be obtained by measuring the resulting polarization.

Statistical description of dust polarized emission from the diffuse interstellar medium

The statistical characterization of the diffuse magnetized interstellar medium (ISM) and Galactic foregrounds to the cosmic microwave background (CMB) poses a major challenge. To account for their non-Gaussian statistics, we need a data analysis approach capable of efficiently quantifying statistical couplings across scales. This information is encoded in the data, but most of it is lost when using conventional tools, such as one-point statistics and power spectra. The wavelet scattering transform (WST), a low-variance statistical descriptor of non-Gaussian processes introduced in data science, opens a path towards this goal. To establish the methodology, we applied the WST to noise-free maps of dust polarized thermal emission computed from a numerical simulation of magnetohydrodynamical turbulence in the diffuse ISM. We analyzed normalized complex Stokes maps and maps of the polarization fraction and polarization angle. The WST yields a few thousand coefficients; some of them measure the amplitude of the signal at a given scale, and the others characterize the couplings between scales and orientations. The dependence on orientation can be fitted with the reduced wavelet scattering transform (RWST), an angular model introduced in previous works for total intensity maps. The RWST provides a statistical description of the polarization maps, quantifying their multiscale properties in terms of isotropic and anisotropic contributions. It allowed us to exhibit the dependence of the map structure on the orientation of the mean magnetic field and to quantify the non-Gaussianity of the data. We also used RWST coefficients, complemented by additional constraints, to generate random synthetic maps with similar statistics. Their agreement with the original maps demonstrates the comprehensiveness of the statistical description provided by the RWST. This work is a step forward in the analysis of observational data and the modeling of CMB foregrounds. We also release PyWST, a public Python package to perform WST and RWST analyses of two-dimensional data.

Detecting pulsar polarization below 100 MHz with the Long Wavelength Array

ABSTRACT Using the first station of the Long Wavelength Array (LWA1), we examine polarized pulsar emission between 25 and 88 MHz. Polarized light from pulsars undergoes Faraday rotation as it passes through the magnetized interstellar medium. Observations from low-frequency telescopes are ideal for obtaining precise rotation measures (RMs) because the effect of Faraday rotation is proportional to the square of the observing wavelength. With these RMs, we obtained polarized pulse profiles to see how polarization changes in the 25–88 MHz range. The RMs were also used to derive values for the electron-density-weighted average Galactic magnetic field along the line of sight. We present RMs and polarization profiles of 15 pulsars acquired using data from LWA1. These results provide new insight into low-frequency polarization characteristics and pulsar emission heights, and complement measurements at higher frequencies.

Interstellar medium magnetic field in bow shock modelling

PurposeThis paper aims to demonstrate the impact of interstellar (IS) magnetic field on stellar shocks existence, shape and size in the stellar wind (SW) vs interstellar medium (ISM) numerical models.Design/methodology/approachComparison of hydrodynamics (HD) and magnetohydrodynamic (MHD) models results with or without ISM magnetic field, its intensity and ISM parameters.FindingsISM magnetic field facilitates formation and stabilises bow shocks around all astrophysical objects. ISM magnetic field may also be one of the reasons for a bow shock existence around the Sun.Practical implicationsISM magnetic field should be implemented in MHD and future kinetic numerical models of the SW interaction with ISM plasma.Originality/valueThis paper presents the results of HD and MHD models of bow shocks and the importance of ISM magnetic field implementation, according to astronomical bow shock observations. The study also presents a review of the most important papers showing the numerical results of bow shock formation.

Statistics of the polarized submillimetre emission maps from thermal dust in the turbulent, magnetized, diffuse ISM

Context. The interstellar medium (ISM) is now widely acknowledged to display features ascribable to magnetized turbulence. With the public release of Planck data and the current balloon-borne and ground-based experiments, the growing amount of data tracing the polarized thermal emission from Galactic dust in the submillimetre provides choice diagnostics to constrain the properties of this magnetized turbulence. Aims. We aim to constrain these properties in a statistical way, focussing in particular on the power spectral index βB of the turbulent component of the interstellar magnetic field in a diffuse molecular cloud, the Polaris Flare. Methods. We present an analysis framework based on simulating polarized thermal dust emission maps using model dust density (proportional to gas density nH) and magnetic field cubes, integrated along the line of sight (LOS), and comparing these statistically to actual data. The model fields are derived from fractional Brownian motion (fBm) processes, which allows a precise control of their one- and two-point statistics. The parameters controlling the model are (1)–(2) the spectral indices of the density and magnetic field cubes, (3)–(4) the RMS-to-mean ratios for both fields, (5) the mean gas density, (6) the orientation of the mean magnetic field in the plane of the sky (POS), (7) the dust temperature, (8) the dust polarization fraction, and (9) the depth of the simulated cubes. We explore the nine-dimensional parameter space through a Markov chain Monte Carlo analysis, which yields best-fitting parameters and associated uncertainties. Results. We find that the power spectrum of the turbulent component of the magnetic field in the Polaris Flare molecular cloud scales with wavenumber as k−βB with a spectral index βB = 2.8 ± 0.2. It complements a uniform field whose norm in the POS is approximately twice the norm of the fluctuations of the turbulent component, and whose position angle with respect to the north-south direction is χ0 ≈−69°. The density field nH is well represented by a log-normally distributed field with a mean gas density 〈nH〉≈40 cm−3, a fluctuation ratio σnH/〈nH〉≈1.6, and a power spectrum with an index βn=1.7−0.3+0.4. We also constrain the depth of the cloud to be d ≈ 13 pc, and the polarization fraction p0 ≈ 0.12. The agreement between the Planck data and the simulated maps for these best-fitting parameters is quantified by a χ2 value that is only slightly larger than unity. Conclusions. We conclude that our fBm-based model is a reasonable description of the diffuse, turbulent, magnetized ISM in the Polaris Flare molecular cloud, and that our analysis framework is able to yield quantitative estimates of the statistical properties of the dust density and magnetic field in this cloud.