Search

Crucial for star formation is the interplay between gravity and turbulence. The observed cloud virial parameter, $\alpha_{\mathrm{vir}}$, which is the ratio of twice the turbulent kinetic energy to the gravitational energy, is found to vary significantly in different environments, where the scatter among individual star-forming clouds can exceed an order of magnitude. Therefore, a strong dependence of the initial mass function (IMF) on $\alpha_{\mathrm{vir}}$ may challenge the notion of a universal IMF. To determine the role of $\alpha_{\mathrm{vir}}$ on the IMF, we compare the star-particle mass functions obtained in high-resolution magnetohydrodynamical simulations including jet and heating feedback, with $\alpha_{\mathrm{vir}}=0.0625$, $0.125$, and $0.5$. We find that varying $\alpha_{\mathrm{vir}}$ from $\alpha_{\mathrm{vir}}\sim0.5$ to $\alpha_{\mathrm{vir}}<0.1$ shifts the peak of the IMF to lower masses by a factor of $\sim2$ and increases the star formation rate by a similar factor. The dependence of the IMF and star formation rate on $\alpha_{\mathrm{vir}}$ is non-linear, with the dependence subsiding at $\alpha_{\mathrm{vir}}<0.1$. Our study shows a systematic dependence of the IMF on $\alpha_{\mathrm{vir}}$. Yet, it may not be measurable easily in observations, considering the uncertainties, and the relatively weak dependence found in this study., Comment: 13 pages, 12 figures, 2 tables (MNRAS submitted)

Even though the interstellar medium (ISM) of star-forming galaxies has been known to have a multiphase structure (broadly hot, warm, and cold phases) since the 1970s, how magnetic fields differ between the ISM phases is still unknown. Using results from numerical simulations, this work explores how the multiphase nature of the ISM shapes magnetic fields and then discusses possible implications of those results for polarisation observations of the Milky Way and high-redshift galaxies. These findings will enhance our understanding of the role of magnetic fields in galaxy evolution and prepare us to harness the upcoming wealth of radio polarisation data from the Square Kilometre Array and its pathfinders., Comment: 8 pages, 3 figures, Conference Proceedings for "Cosmic Magnetism in the Pre-SKA Era", May 27 - 31, 2024, Kagoshima, Japan

We present the largest Galactic neutral hydrogen HI absorption survey to date, utilizing the Australian SKA Pathfinder Telescope at an unprecedented spatial resolution of 30''. This survey, GASKAP-HI, unbiasedly targets 2,714 continuum background sources over 250 square degrees in the direction of the Magellanic Clouds, a significant increase compared to a total of 373 sources observed by previous Galactic absorption surveys across the entire Milky Way. We aim to investigate the physical properties of cold (CNM) and warm (WNM) neutral atomic gas in the Milky Way foreground, characterized by two prominent filaments at high Galactic latitudes (between $-45^{\circ}$ and $-25^{\circ}$). We detected strong HI absorption along 462 lines of sight above the 3$\sigma$ threshold, achieving an absorption detection rate of 17%. GASKAP-HI's unprecedented angular resolution allows for simultaneous absorption and emission measurements to sample almost the same gas clouds along a line of sight. A joint Gaussian decomposition is then applied to absorption-emission spectra to provide direct estimates of HI optical depths, temperatures, and column densities for the CNM and WNM components. The thermal properties of CNM components are consistent with those previously observed along a wide range of Solar neighborhood environments, indicating that cold HI properties are widely prevalent throughout the local interstellar medium. Across our region of interest, CNM accounts for ~30% of the total HI gas, with the CNM fraction increasing with column density toward the two filaments. Our analysis reveals an anti-correlation between CNM temperature and its optical depth, which implies that CNM with lower optical depth leads to a higher temperature., Comment: Largest Galactic HI Absorption Survey To Date (GASKAP-HI): Cold Atomic Gas in the Magellanic Cloud foreground using Australian SKA Pathfinder. This paper has 19 pages, 17 figures. This paper has been accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal

Radio observations have revealed magnetic fields in the intracluster medium (ICM) of galaxy clusters nearly in equipartition with turbulence. This suggests magnetic field amplification by dynamo processes. However, observations are limited to redshifts z <~ 0.7, and the weakly collisional nature of the ICM complicates studying magnetic field evolution at higher redshifts through theoretical models and simulations. Using a model of the weakly collisional dynamo, we modelled the evolution of the Faraday rotation measure (RM) in clusters of different masses, up to z <~ 1.5, and investigated its properties. We compared our results with radio observations of various galaxy clusters. We used merger trees generated by the modified GALFORM algorithm to track the evolution of plasma quantities during galaxy cluster formation. Assuming the magnetic field remains in equipartition with the turbulent velocity field, we generated RM maps to study their properties. We find that both the standard deviation of RM, sigma_RM, and the absolute average |mu_RM| increase with cluster mass. Due to redshift dilution, RM values for a fixed cluster mass remain nearly constant between z=0 and z=1.5. For r/r200 >~ 0.4, sigma_RM does not vary significantly with L/r200, with L being the size of the observed RM patch. Below this limit, sigma_RM increases as L decreases. We find that radial RM profiles have a consistent shape, proportional to 10^{-1.2(r/r200)}, and are nearly independent of redshift. Our z~0 profiles for M_clust = 10^15 Msol match RM observations in the Coma cluster but show discrepancies with Perseus, possibly due to high gas mixing. Models for clusters with M_clust = 10^13 and M_clust = 10^15 Msol at z=0 and z = 0.174 align well with Fornax and A2345 data for r/r200 <~ 0.4. Our model can be useful for generating mock polarization observations for radio telescopes., Comment: 14 pages, 11 figures

We introduce $\texttt{A}$strophysical $\texttt{H}$ybrid-$\texttt{K}$inetic simulations with the $\texttt{flASH}$ code ($\texttt{AHKASH}$) -- a new Hybrid particle-in-cell (PIC) code developed within the framework of the multi-physics code $\texttt{FLASH}$. The new code uses a second-order accurate Boris integrator and a predictor-predictor-corrector algorithm for advancing the Hybrid-kinetic equations, using the constraint transport method to ensure that magnetic fields are divergence-free. The code supports various interpolation schemes between the particles and grid cells, with post-interpolation smoothing to reduce finite particle noise. We further implement a $\delta f$ method to study instabilities in weakly collisional plasmas. The new code is tested on standard physical problems such as the motion of charged particles in uniform and spatially varying magnetic fields, the propagation of Alfv\'en and whistler waves, and Landau damping of ion acoustic waves. We test different interpolation kernels and demonstrate the necessity of performing post-interpolation smoothing. We couple the $\texttt{TurbGen}$ turbulence driving module to the new Hybrid PIC code, allowing us to test the code on the highly complex physical problem of the turbulent dynamo. To investigate steady-state turbulence with a fixed sonic Mach number, it is important to maintain isothermal plasma conditions. Therefore, we introduce a novel cooling method for Hybrid PIC codes and provide tests and calibrations of this method to keep the plasma isothermal. We describe and test the `hybrid precision' method, which significantly reduces (by a factor $\sim1.5$) the computational cost, without compromising the accuracy of the numerical solutions. Finally, we test the parallel scalability of the new code, showing excellent scaling up to 10,000~cores., Comment: 23 pages, 20 figures, Accepted for publication in Monthly Notices of the Royal Astronomical Society

Magnetic fields are an energetically important component of star-formation galaxies, but it is often difficult to measure their properties from observations. One of the complexities stems from the fact that the magnetic fields, especially in spiral galaxies, have a two-scale nature: a large-scale field, coherent over ${\rm kpc}$ scales and a small-scale, random field with a scale of $\lesssim$ $100~{\rm pc}$. Moreover, it is known that the strength of small- and large-scale fields are comparable and this makes it even harder to find their imprints in radio polarisation observations such as the Faraday rotation measure, ${\rm RM}$, which is the integral over the path length of the product of the thermal electron density and the parallel component of the magnetic field to the line of sight. Here, we propose and demonstrate the use of second-order structure functions of ${\rm RM}$ computed with multiple higher-order stencils as a powerful analysis to separate the small- and large-scale magnetic field components. In particular, we provide new methods and calibrations to compute the scale and the strength of the large-scale magnetic field in the presence of small-scale magnetic fluctuations. We then apply the method to find the scale of large-scale magnetic fields in the nearby galaxies M51 and NGC 6946, using archival data and further discuss the need for computing the ${\rm RM}$ structure functions with higher-order stencils. With multiple modern radio polarisation observatories and eventually the Square Kilometre Array, ${\rm RM}$ observations will significantly improve in quantity and quality, and the higher-order stencil structure function techniques developed here can be used to extract information about multiscale magnetic fields in galaxies., Comment: 12 pages, 7 figures, 1 table; accepted for publication in MNRAS

We study the wide-binary eccentricity ($e$) distribution in young star clusters and the role of turbulence in setting the form of the $e$ distribution using magnetohydrodynamical (MHD) simulations of star cluster formation. The simulations incorporate gravity, turbulence, magnetic fields, protostellar heating, and jets/outflows. We find that (1) simulations that employ purely compressive turbulence driving produce binaries with a superthermal $e$ distribution ($\alpha>1$ in $p(e) \propto e^\alpha$), while simulations with purely solenoidal driving or natural mixture of driving modes produce subthermal/thermal distributions ($\alpha \leq$ 1), (2) the $e$ distribution over the full range of binary separations in our simulations is set at the early stages of the star cluster formation process, (3) while binaries (separation of $r_{\mathrm{pair}} \leq 1000\, \mathrm{AU}$) have subthermal to thermal $e$ distributions ($\alpha \sim 0.8$), wide binaries ($r_{\mathrm{pair}} > 1000\, \mathrm{AU}$) have a superthermal distribution ($\alpha \sim 1.8$), and (4) low-mass binary systems (system masses of $M_{\mathrm{sys}} \leq 0.8\, \mathrm{M_\odot}$) have a highly superthermal distribution ($\alpha \sim 2.4$), whereas high-mass systems ($M_{\mathrm{sys}} > 0.8\, \mathrm{M_\odot}$) exhibit a subthermal/thermal distribution ($\alpha \sim 0.8$). The binary eccentricity distribution is often modelled as a thermal distribution. However, our results suggest that the $e$ distribution depends on the range of separation of the sampled binaries, which agrees with the findings from recent Gaia observations. We conclude that the dependence of the $e$ distribution on the binary separation and mass is linked to the binary formation mechanism governed by the turbulent properties of the parent cloud., Comment: 15 pages, 9 figures, 1 table (added additional citation, published in MNRAS)

Molecular clouds are prime locations to study the process of star formation. These clouds contain filamentary structures and cores, which are crucial sites for the formation of young stars. The star-formation process has been investigated using various techniques, including polarimetry for tracing magnetic fields. In this small review-cum-short report, we put together the efforts (mainly from the Indian community) to understand the roles of turbulence and magnetic fields in star formation. These are two components of the ISM competing against gravity, which is primarily responsible for the collapse of gas to form stars. We also include attempts made using simulations of molecular clouds to study this competition. Studies on feedback and magnetic fields are combined and listed to understand the importance of the interaction between two energies in setting the current observed star formation efficiency. We have listed available and upcoming facilities with the polarization capabilities needed to trace magnetic fields. We have also stated the importance of ongoing and desired collaborations between Indian communities and facilities abroad to shed more light on the roles of turbulence and magnetic fields in the process of star formation., Comment: 12 pages, 2 figures, 1 table, accepted in Journal of Astrophysics and Astronomy

Several different methods are regularly used to infer the properties of the neutral interstellar medium (ISM) using atomic hydrogen (H I) 21cm absorption and emission spectra. In this work, we study various techniques used for inferring ISM gas phase properties, namely the correlation between brightness temperature and optical depth $(T_B(v)$, $\tau(v))$ at each channel velocity ($v$), and decomposition into Gaussian components, by creating mock spectra from a 3D magnetohydrodynamic simulation of a two-phase, turbulent ISM. We propose a physically motivated model to explain the $T_B(v)-\tau(v)$ distribution and relate the model parameters to properties like warm gas spin temperature and cold cloud length scales. Two methods based on Gaussian decomposition -- using only absorption spectra and both absorption and emission spectra -- are used to infer the column density distribution as a function of temperature. In observations, such analysis reveals the puzzle of large amounts (significantly higher than in simulations) of gas with temperature in the thermally unstable range of $\sim$200 K to $\sim$2000 K and a lack of the expected bimodal (two-phase) temperature distribution. We show that, in simulation, both methods are able to recover the actual gas distribution in the simulation till temperatures $\lesssim2500$~K (and the two-phase distribution in general) reasonably well. We find our results to be robust to a range of effects such as noise, varying emission beam size, and simulation resolution. This shows that the observational inferences are unlikely to be artifacts, thus highlighting a tension between observations and simulations. We discuss possible reasons for this tension and ways to resolve it., Comment: 22 pages (including appendixes), 16 figures, 3 tables, Accepted for publication in MNRAS

The Australian SKA Pathfinder (ASKAP) radio telescope has carried out a survey of the entire Southern Sky at 887.5MHz. The wide area, high angular resolution, and broad bandwidth provided by the low-band Rapid ASKAP Continuum Survey (RACS-low) allow the production of a next-generation rotation measure (RM) grid across the entire Southern Sky. Here we introduce this project as Spectral and Polarisation in Cutouts of Extragalactic sources from RACS (SPICE-RACS). In our first data release, we image 30 RACS-low fields in Stokes $I$, $Q$, $U$ at 25'' angular resolution, across 744 to 1032MHz with 1MHz spectral resolution. Using a bespoke, highly parallelised, software pipeline we are able to rapidly process wide-area spectro-polarimetric ASKAP observations. Notably, we use 'postage stamp' cutouts to assess the polarisation properties of \ncomponents\ radio components detected in total intensity. We find that our Stokes $Q$ and $U$ images have an rms noise of ~80$\mu$Jy/PSF, and our correction for instrumental polarisation leakage allows us to characterise components with >1% polarisation fraction over most of the field of view. We produce a broadband polarised radio component catalogue that contains \nrms\ RM measurements over an area of ~1300deg^2 with an average error in RM of 1.6+1.1-1.0rad/m^2, and an average linear polarisation fraction 3.4+3.0-1.6%. We determine this subset of components using the conditions that the polarised signal-to-noise ratio is $>8$, the polarisation fraction is above our estimated polarised leakage, and the Stokes $I$ spectrum has a reliable model. Our catalogue provides an areal density of $4\pm2$ RMs/deg^2; an increase of $\sim4$ times over the previous state-of-the-art (Taylor et al. 2009). Meaning that, having used just 3% of the RACS-low sky area, we have produced the 3rd largest RM catalogue to date. This catalogue has broad applications for studying..., Comment: 42 pages, 24 figures, 6 tables. Accepted for publications in PASA

The turbulent small-scale dynamo (SSD) is likely to be responsible for the magnetisation of the interstellar medium (ISM) that we observe in the Universe today. The SSD efficiently converts kinetic energy $E_{\rm kin}$ into magnetic energy $E_{\rm mag} $, and is often used to explain how an initially weak magnetic field with $E_{\rm mag} \ll E_{\rm kin}$ is amplified, and then maintained at a level $E_{\rm mag} \lesssim E_{\rm kin}$. Usually, this process is studied by initialising a weak seed magnetic field and letting the turbulence grow it to saturation. However, in this Part I of the Growth or Decay series, using three-dimensional, visco-resistive magnetohydrodynamical turbulence simulations up to magnetic Reynolds numbers of 2000, we show that the same final state in the integral quantities, energy spectra, and characteristic scales of the magnetic field can also be achieved if initially $E_{\rm mag} \sim E_{\rm kin}$ or even if initially $E_{\rm mag} \gg E_{\rm kin}$. This suggests that the final saturated state of the turbulent dynamo is set by the turbulence and the material properties of the plasma, independent of the initial structure or amplitude of the magnetic field. We discuss the implications this has for the maintenance of magnetic fields in turbulent plasmas and future studies exploring the dynamo saturation., Comment: 16 pages. 11 figures. Accepted to MNRAS. Part I focuses on the journey towards saturation. Part II focuses on the decay into saturation

Turbulence is a critical ingredient for star formation, yet its role for the initial mass function (IMF) is not fully understood. Here we perform magnetohydrodynamical (MHD) simulations of star cluster formation including gravity, turbulence, magnetic fields, stellar heating and outflow feedback to study the influence of the mode of turbulence driving on IMF. We find that simulations that employ purely compressive turbulence driving (COMP) produce a higher fraction of low-mass stars as compared to simulations that use purely solenoidal driving (SOL). The characteristic (median) mass of the sink particle (protostellar) distribution for COMP is shifted to lower masses by a factor of ~ 1.5 compared to SOL. Our simulation IMFs capture the important features of the observed IMF form. We find that turbulence-regulated theories of the IMF match our simulation IMFs reasonably well in the high-mass and low-mass range, but underestimate the number of very low-mass stars, which form towards the later stages of our simulations and stop accreting due to dynamical interactions. Our simulations show that for both COMP and SOL, the multiplicity fraction is an increasing function of the primary mass, although the multiplicity fraction in COMP is higher than that of SOL for any primary mass range. We find that binary mass ratio distribution is independent of the turbulence driving mode. The average specific angular momentum of the sink particles in SOL is a factor of 2 higher than that for COMP. Overall, we conclude that the turbulence driving mode plays a significant role in shaping the IMF., Comment: 26 pages, 16 figures, 4 tables, accepted for publication in MNRAS, minor changes, part of methodology is same as in Mathew & Federrath 2021. arXiv admin note: text overlap with arXiv:2106.06521

Magnetic fields and turbulence are important components of the interstellar medium (ISM) of star-forming galaxies. It is challenging to measure the properties of the small-scale ISM magnetic fields (magnetic fields at scales smaller than the turbulence driving scale). Using numerical simulations, we demonstrate how the second-order rotation measure (RM, which depends on thermal electron density, $n_{\rm e}$, and magnetic field, $b$) structure function can probe the properties of small-scale $b$. We then apply our results to observations of the Small and Large Magellanic Clouds (SMC and LMC). First, using Gaussian random $b$, we show that the characteristic scale where the RM structure function flattens is approximately equal to the correlation length of $b$. We also show that computing the RM structure function with a higher-order stencil (more than the commonly-used two-point stencil) is necessary to accurately estimate the slope of the structure function. Then, using Gaussian random $b$ and lognormal $n_{\rm e}$ with known power spectra, we derive an empirical relationship between the slope of the power spectrum of $b$, $n_{\rm e}$, and RM. We apply these results to the SMC and LMC and estimate the following properties of small-scale $b$: correlation length ($160~\pm 21~{\rm pc}$ for the SMC and $87~\pm~17~{\rm pc}$ for the LMC), strength ($14~\pm 2~\mu{\rm G}$ for the SMC and $15~\pm 3~\mu{\rm G}$ for the LMC), and slope of the magnetic power spectrum ($-1.3~\pm~0.4$ for the SMC and $-1.6~\pm~0.1$ for the LMC). We also find that $n_{\rm e}$ is practically constant over the estimated $b$ correlation scales., Comment: 27 pages (including 6 appendices), 29 figures (16 in the main text and 13 in the appendices), and 4 tables (2 in the main text and 2 in the appendices); MNRAS accepted

Cosmic rays (CRs) are a dynamically important component of the interstellar medium (ISM) of galaxies. The $\sim$GeV CRs that carry most CR energy and pressure are likely confined by self-generated turbulence, leading them to stream along magnetic field lines at the ion Alfv\'en speed. However, the consequences of self-confinement for CR propagation on galaxy scales remain highly uncertain. In this paper, we use a large ensemble of magnetohydrodynamical turbulence simulations to quantify how the basic parameters describing ISM turbulence -- the sonic Mach number, $\mathcal{M}$ (plasma compressibility), Alfv\'en Mach number, $\mathcal{M}_{A0}$ (strength of the large-scale field with respect to the turbulence), and ionisation fraction by mass, $\chi$ -- affect the transport of streaming CRs. We show that the large-scale transport of CRs whose small-scale motion consists of streaming along field lines is well described as a combination of streaming along the mean field and superdiffusion both along (parallel to) and across (perpendicular to) it; $\mathcal{M}_{A0}$ drives the level of anisotropy between parallel and perpendicular diffusion and $\chi$ modulates the magnitude of the diffusion coefficients, while in our choice of units, $\mathcal{M}$ is unimportant except in the sub-Alfv\'enic ($\mathcal{M}_{A0} \lesssim 0.5$) regime. Our finding that superdiffusion is ubiquitous potentially explains the apparent discrepancy between CR diffusion coefficients inferred from measurements close to individual sources compared to those measured on larger, Galactic scales. Finally, we present empirical fits for the diffusion coefficients as a function of plasma parameters that may be used as sub-grid recipes for global interstellar medium, galaxy or cosmological simulations., Comment: Accepted in MNRAS

The turbulent dynamo is a powerful mechanism that converts turbulent kinetic energy to magnetic energy. A key question regarding the magnetic field amplification by turbulence, is, on what scale, $k_{\rm p}$, do magnetic fields become most concentrated? There has been some disagreement about whether $k_{\rm p}$ is controlled by the viscous scale, $k_\nu$ (where turbulent kinetic energy dissipates), or the resistive scale, $k_\eta$ (where magnetic fields dissipate). Here we use direct numerical simulations of magnetohydrodynamic turbulence to measure characteristic scales in the kinematic phase of the turbulent dynamo. We run $104$-simulations with hydrodynamic Reynolds numbers of $10 \leq {\rm Re} \leq 3600$, and magnetic Reynolds numbers of $270 \leq {\rm Rm} \leq 4000$, to explore the dependence of $k_{\rm p}$ on $k_\nu$ and $k_\eta$. Using physically motivated models for the kinetic and magnetic energy spectra, we measure $k_\nu$, $k_\eta$ and $k_{\rm p}$, making sure that the obtained scales are numerically converged. We determine the overall dissipation scale relations $k_\nu = (0.025^{+0.005}_{-0.006})\, k_{\rm turb}\, {\rm Re}^{3/4}$ and $k_\eta = (0.88^{+0.21}_{-0.23})\, k_\nu\, {\rm Pm}^{1/2}$, where $k_{\rm turb}$ is the turbulence driving wavenumber and ${\rm Pm}={\rm Rm}/{\rm Re}$ is the magnetic Prandtl number. We demonstrate that the principle dependence of $k_{\rm p}$ is on $k_\eta$. For plasmas where ${\rm Re} \gtrsim 100$, we find that $k_{\rm p} = (1.2_{-0.2}^{+0.2})\, k_\eta$, with the proportionality constant related to the power-law `Kazantsev' exponent of the magnetic power spectrum. Throughout this study, we find a dichotomy in the fundamental properties of the dynamo where ${\rm Re} > 100$, compared to ${\rm Re} < 100$. We report a minimum critical hydrodynamic Reynolds number, ${\rm Re}_{\rm crit} = 100$ for bonafide turbulent dynamo action., Comment: 15 pages, 10 figures; accepted for publication in MNRAS, 2022 April 4

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\mathbf{B}$ and large-scale, $\mathbf{B}_0$ fields, $\left< (\delta\mathbf{B}\cdot\mathbf{B}_0)^{2} \right>^{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\'enic 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 MHD turbulence. Furthermore, we explore the relation between the turbulent, mean-field and total Alfv\'en Mach numbers, and demonstrate that sub-Alfv\'enic turbulence can only be developed through a strong, large-scale magnetic field, which supports an extremely super-Alfv\'enic turbulent magnetic field. This means that the magnetic field fluctuations are significantly subdominant to the velocity fluctuations in the sub-Alfv\'enic large-scale field regime. Throughout our study, we broadly discuss the implications for observations of magnetic fields and understanding the dynamics in the magnetised ISM., Comment: Accepted in MNRAS. 20 pages. 16 figures. 1 table

Magnetic fields are a dynamically important component of the turbulent interstellar medium (ISM) of star-forming galaxies. These magnetic fields are due to a dynamo action, which is a process of converting turbulent kinetic energy to magnetic energy. A dynamo that acts at scales less than the turbulent driving scale is known as the turbulent dynamo. The ISM is a multiphase medium and observations suggest that the properties of magnetic fields differ with the phase. Here, we aim to study how the properties of the turbulent dynamo depend on the phase. We simulate the non-isothermal turbulent dynamo in a two-phase medium (most previous work assumes an isothermal gas). We show that the warm phase ($T\ge10^3~{\rm K}$) is transsonic and the cold phase ($T<10^3~{\rm K}$) is supersonic. We find that the growth rate of magnetic fields in the exponentially growing stage is similar in both phases. We compute the terms responsible for amplification and destruction of vorticity and show that in both phases vorticity is amplified due to turbulent motions, further amplified by the baroclinic term in the warm phase, and destroyed by the term for viscous interactions in the presence of logarithmic density gradients in the cold phase. We find that the final ratio of magnetic to turbulent kinetic energy is lower in the cold phase due to a stronger Lorentz force. We show that the non-isothermal turbulent dynamo is significantly different from its isothermal counterpart and this demonstrates the need for studying the turbulent dynamo in a multiphase medium., Comment: 18 pages (including 3 appendices), 17 figures (12 in the main text and 5 in the appendices), and 3 tables; accepted in MNRAS

We present results from a search for the HI 21-cm line in absorption towards 16 bright radio sources with the 6-antenna commissioning array of the Australian Square Kilometre Array Pathfinder (ASKAP). Our targets were selected from the 2-Jy sample, a flux-limited survey of the southern radio sky with extensive multi-wavelength follow-up. Two sources were detected in HI absorption including a new detection towards the bright FRII radio galaxy PKS 0409-75 at a redshift of z=0.674. The HI absorption line is blueshifted by ~3300 km/s compared to the optical redshift of the host galaxy of PKS 0409-75 at z=0.693. Deep optical imaging and spectroscopic follow-up with the GMOS instrument on the Gemini-South telescope reveal that the HI absorption is associated with a galaxy in front of the southern radio lobe with a stellar mass of $3.2 - 6.8 \times 10^{11}M_\odot$, a star-formation rate of $\sim 1.24 M_\odot$ yr$^{-1}$ and an estimated HI column density of $2.16\times10^{21}$ cm$^{-2}$, assuming a spin temperature of $T_{\rm spin}=500$ K and source covering factor of $C_{\rm f}=0.3$. Using polarisation measurements of PKS 0409-75 from the literature we estimate the magnetic field of the absorbing galaxy to be ~14.5$\mu$G, consistent with field strengths observed in nearby spiral galaxies, but larger than expected for an elliptical galaxy. Results from this pilot study can inform future surveys as new wide-field telescopes allow us to search for 21-cm HI absorption towards all bright radio sources as opposed to smaller targeted samples., Comment: 14 pages, 6 figures, accepted for publication in MNRAS

Magnetic fields in several astrophysical objects are amplified and maintained by a dynamo mechanism, which is the conversion of the turbulent kinetic energy to magnetic energy. A dynamo that amplifies magnetic fields at scales $<$ the driving scale of turbulence is known as the fluctuation dynamo. We study the properties of the fluctuation dynamo in supersonic turbulent plasmas, which is of relevance to the ISM, structure formation, and lab experiments of laser-plasma turbulence. Using simulations, we explore the properties of the exponentially growing and saturated state of the fluctuation dynamo for subsonic and supersonic turbulence. We confirm that the fluctuation dynamo efficiency decreases with compressibility. We show that the fluctuation dynamo generated magnetic fields are spatially intermittent and the level of intermittency decreases as the field saturates. We find a stronger back reaction of the magnetic field on the velocity for the subsonic case as compared to the supersonic case. Locally, we find that the level of alignment between vorticity and velocity, velocity and magnetic field, and current density and magnetic field in the saturated stage is enhanced in comparison to the exponentially growing phase for the subsonic case, but only the current density and magnetic field alignment is enhanced for the supersonic case. We show that both the magnetic field amplification (due to weaker stretching of field lines) and diffusion decreases when the field saturates, but the diffusion is enhanced relative to amplification. This occurs throughout the volume in the subsonic turbulence, but primarily in the strong-field regions for the supersonic case. This leads to the saturation of the fluctuation dynamo. Overall, both the amplification and diffusion of magnetic fields are affected and thus a drastic change in either of them is not required for the saturation. [Abstract abridged], Comment: 24 pages, 17 figures, 2 tables; Accepted for publication in Physical Review Fluids

Studying the driving modes of turbulence is important for characterizing the impact of turbulence in various astrophysical environments. The driving mode of turbulence is parameterized by $b$, which relates the width of the gas density PDF to the turbulent Mach number; $b\approx 1/3$, $1$, and $0.4$ correspond to driving that is solenoidal, compressive, and a natural mixture of the two, respectively. In this work, we use high-resolution (sub-pc) ALMA $^{12}$CO ($J$ = $2-1$), $^{13}$CO ($J$ = $2-1$), and C$^{18}$O ($J$ = $2-1$) observations of filamentary molecular clouds in the star-forming region N159E (the Papillon Nebula) in the Large Magellanic Cloud (LMC) to provide the first measurement of turbulence driving parameter in an extragalactic region. We use a non-local thermodynamic equilibrium (NLTE) analysis of the CO isotopologues to construct a gas density PDF, which we find to be largely log-normal in shape with some intermittent features indicating deviations from lognormality. We find that the width of the log-normal part of the density PDF is comparable to the supersonic turbulent Mach number, resulting in $b \approx 0.9$. This implies that the driving mode of turbulence in N159E is primarily compressive. We speculate that the compressive turbulence could have been powered by gravo-turbulent fragmentation of the molecular gas, or due to compression powered by H I flows that led to the development of the molecular filaments observed by ALMA in the region. Our analysis can be easily applied to study the nature of turbulence driving in resolved star-forming regions in the local as well as the high-redshift Universe., Comment: 14 pages, 7 figures. MNRAS in press

Magnetic fields in elliptical galaxies are poorly constrained due to a lack of significant synchrotron emission from them. This paper explores properties of magnetic fields in ellipticals via two methods. First, we exploit the Laing-Garrington effect (asymmetry in the observed polarization fraction between radio galaxy jets) for 57 galaxies with redshifts up to $0.5$. We use the differences in polarization fraction and rotation measure between the jet and counterjet to estimate the small- and large-scale magnetic fields in and around ellipticals (including their circumgalactic medium). We find that the small-scale field (at scales smaller than the driving scale of turbulence, approximately $300~{\rm pc}$) lies in the range $0.1~\text{--}~2.75~\mu{\rm G}$. The large-scale field (at scales of $100~{\rm kpc}$) is an order of magnitude smaller than the small-scale field. In the second method, we cross-match the Faraday rotation measures (RM) of a few hundred (out of $3098$) extragalactic radio sources with galaxy catalogs to explore the effect of the number and morphology of intervening galaxies on the observed RM distribution. We use both Gaussian and non-Gaussian functions to describe the RM distribution and derive its statistical properties. Finally, using the difference in the observed polarization fraction between the intervening spirals and ellipticals, we estimate the small-scale magnetic fields at the center of ellipticals to be $\sim6~\mu{\rm G}$. Both methods with different observations and analysis techniques give magnetic field strengths consistent with previous studies ($\leq10\mu{\rm G}$), and the results can be used to constrain dynamo theories and galaxy evolution simulations., Comment: 18 pages (incl. 1 appendix), total 12 figures (incl. 1 in appendix), published in MNRAS

Pulsars can act as an excellent probe of the Milky Way magnetic field. The average strength of the Galactic magnetic field component parallel to the line of sight can be estimated as $\langle B_\parallel \rangle = 1.232 \, \text{RM}/\text{DM}$, where $\text{RM}$ and $\text{DM}$ are the rotation and dispersion measure of the pulsar. However, this assumes that the thermal electron density and magnetic field of the interstellar medium are uncorrelated. Using numerical simulations and observations, we test the validity of this assumption. Based on magnetohydrodynamical simulations of driven turbulence, we show that the correlation between the thermal electron density and the small-scale magnetic field increases with increasing Mach number of the turbulence. We find that the assumption of uncorrelated thermal electron density and magnetic fields is valid only for subsonic and transsonic flows, but for supersonic turbulence, the field strength can be severely overestimated by using $1.232 \, \text{RM}/\text{DM}$. We then correlate existing pulsar observations from the Australia Telescope National Facility with regions of enhanced thermal electron density and magnetic fields probed by ${^{12} \mathrm {CO}}$ data of molecular clouds, magnetic fields from the Zeeman splitting of the 21 cm line, neutral hydrogen column density, and H$\alpha$ observations. Using these observational data, we show that the thermal electron density and magnetic fields are largely uncorrelated over kpc scales. Thus, we conclude that the relation $\langle B_\parallel \rangle = 1.232 \, \text{RM}/\text{DM}$ provides a good estimate of the magnetic field on Galactic scales but might break down on sub - kpc scales., Comment: Accepted for publication in MNRAS; 18 pages (including one appendix), 13 figures (12 in the main text and 1 in the appendix), 6 tables in the main text

Fluctuation dynamos are thought to play an essential role in magnetized galaxy evolution, saturating within $\sim0.01~$Gyr and thus potentially acting as seeds for large-scale dynamos. However, unambiguous observational confirmation of the fluctuation dynamo action in a galactic environment is still missing. This is because, in spiral galaxies, it is difficult to differentiate between small-scale magnetic fields generated by a fluctuation dynamo and those due to the tangling of the large-scale field. We propose that observations of magnetic fields in elliptical galaxies would directly probe the fluctuation dynamo action. This is motivated by the fact that in ellipticals, due to their lack of significant rotation, the conventional large-scale dynamo is absent and the fluctuation dynamo is responsible for controlling the magnetic field strength and structure. By considering turbulence injected by Type Ia supernova explosions and possible magnetic field amplification by cooling flows, we estimate expected magnetic field strengths of $0.2~-~1 \,\mu{\rm G}$ in quiescent elliptical galaxies. We use a semi-analytic model of galaxy formation to estimate the distribution and redshift evolution of field strengths, tentatively finding a decrease in magnetic field strength with decreasing redshift. We analyse a sample of radio sources that exhibit the Laing-Garrington (LG) effect (radio polarization asymmetry in jets) and infer magnetic field strengths between $0.14~-~1.33 \,\mu{\rm G}$ for a uniform thermal electron density and between $1.36~-~6.21\,\mu{\rm G}$ for the thermal electron density following the King profile. We examine observational techniques for measuring the magnetic field saturation state in elliptical galaxies, focusing on Faraday RM grids, the LG effect, synchrotron emission, and gravitational lensing, finding appealing prospects for future empirical analysis., Comment: 16 pages, 7 figures, 1 table, accepted for publication in The Astrophysical Journal

Magnetic fields in galaxies and galaxy clusters are amplified from a very weak seed value to the observed $\mu{\rm G}$ strengths by the turbulent dynamo. The seed magnetic field can be of primordial or astrophysical origin. The strength and structure of the seed field, on the galaxy or galaxy cluster scale, can be very different, depending on the seed-field generation mechanism. The seed field first encounters the small-scale dynamo, thus we investigate the effects of the strength and structure of the seed field on the small-scale dynamo action. Using numerical simulations of driven turbulence and considering three different seed-field configurations: 1) uniform field, 2) random field with a power-law spectrum, and 3) random field with a parabolic spectrum, we show that the strength and statistical properties of the dynamo-generated magnetic fields are independent of the details of the seed field. We demonstrate that, even when the small-scale dynamo is not active, small-scale magnetic fields can be generated and amplified linearly due to the tangling of the large-scale field. In the presence of the small-scale dynamo action, we find that any memory of the seed field for the non-linear small-scale dynamo generated magnetic fields is lost and thus, it is not possible to trace back seed-field information from the evolved magnetic fields in a turbulent medium., Comment: 12 pages, 9 figures, 1 table, accepted for publication in MNRAS

The rich structure that we observe in molecular clouds is due to the interplay between strong magnetic fields and supersonic (turbulent) velocity fluctuations. The velocity fluctuations interact with the magnetic field, causing it too to fluctuate. Using numerical simulations, we explore the nature of such magnetic field fluctuations, $\vec{\delta B}$, over a wide range of turbulent Mach numbers, $\mathcal{M} = 2 - 20$ (i.e., from weak to strong compressibility), and Alfv\'en Mach numbers, $\mathcal{M}_{\text{A}0} = 0.1 - 100$ (i.e., from strong to weak magnetic mean fields, $B_0$). We derive a compressible quasi-static fluctuation model from the magnetohydrodynamical (MHD) equations and show that velocity gradients parallel to the mean magnetic field give rise to compressible modes in sub-Alfv\'enic flows, which prevents the flow from becoming two-dimensional, as is the case in incompressible MHD turbulence. We then generalise an analytical model for the magnitude of the magnetic fluctuations to include $\mathcal{M}$, and find $|\vec{\delta B}| = \delta B = c_s\sqrt{\pi\rho_0}\mathcal{M}\mathcal{M}_{\text{A}0}$, where $c_s$ is the sound speed and $\rho_0$ is the mean density of gas. This new relation fits well in the strong $B$-field regime. We go on to study the anisotropy between the perpendicular ($ B_{\perp}$) and parallel ($ B_{\parallel}$) fluctuations and the mean-normalised fluctuations, which we find follow universal scaling relations, invariant of $\mathcal{M}$. We provide a detailed analysis of the morphology for the $\delta B_{\perp}$ and $\delta B_{\parallel}$ probability density functions and find that eddies aligned with $B_0$ cause parallel fluctuations that reduce $B_{\parallel}$ in the most anisotropic simulations. We discuss broadly the implications of our fluctuation models for magnetised gases in the interstellar medium., Comment: 19 pages, 9 figures. Accepted for publication in MNRAS

The presence of magnetic fields in many astrophysical objects is due to dynamo action, whereby a part of the kinetic energy is converted into magnetic energy. A turbulent dynamo that produces magnetic field structures on the same scale as the turbulent flow is known as the fluctuation dynamo. We use numerical simulations to explore the nonlinear, statistically steady state of the fluctuation dynamo in driven turbulence. We demonstrate that as the magnetic field growth saturates, its amplification and diffusion are both affected by the back-reaction of the Lorentz force upon the flow. The amplification of the magnetic field is reduced due to stronger alignment between the velocity field, magnetic field, and electric current density. Furthermore, we confirm that the amplification decreases due to a weaker stretching of the magnetic field lines. The enhancement in diffusion relative to the field line stretching is quantified by a decrease in the computed local value of the magnetic Reynolds number. Using the Minkowski functionals, we quantify the shape of the magnetic structures produced by the dynamo as magnetic filaments and ribbons in both kinematic and saturated dynamos and derive the scalings of the typical length, width, and thickness of the magnetic structures with the magnetic Reynolds number. We show that all three of these magnetic length scales increase as the dynamo saturates. The magnetic intermittency, strong in the kinematic dynamo (where the magnetic field strength grows exponentially) persists in the statistically steady state, but intense magnetic filaments and ribbons are more volume-filling., Comment: 17 pages, 19 figures, accepted for publication in Physical Review Fluids

Energy equipartition between cosmic rays and magnetic fields is often assumed to infer magnetic field properties from the synchrotron observations of star-forming galaxies. However, there is no compelling physical reason to expect the same. We aim to explore the validity of the energy equipartition assumption. After describing popular arguments in favour of the assumption, we first discuss observational results which support it at large scales and how certain observations show significant deviations from equipartition at scales smaller than $\approx 1 \, {\rm kpc}$, probably related to the propagation length of the cosmic rays. Then we test the energy equipartition assumption using test-particle and MHD simulations. From the results of the simulations, we find that the energy equipartition assumption is not valid at scales smaller than the driving scale of the ISM turbulence ($\approx 100 \, {\rm pc}$ in spiral galaxies), which can be regarded as the lower limit for the scale beyond which equipartition is valid. We suggest that one must be aware of the dynamical scales in the system before assuming energy equipartition to extract magnetic field information from synchrotron observations. Finally, we present ideas for future observations and simulations to investigate in more detail under which conditions the equipartition assumption is valid or not., Comment: Invited review article for the special issue "New Perspectives on Galactic Magnetism" of the journal "Galaxies", accepted for publication

Interpretations of synchrotron observations often assume a tight correlation between magnetic and cosmic ray energy densities. We examine this assumption using both test-particle simulations of cosmic rays and MHD simulations which include cosmic rays as a diffusive fluid. We find no spatial correlation between the cosmic rays and magnetic field energy densities at turbulent scales. Moreover, the cosmic ray number density and magnetic field energy density are statistically independent. Nevertheless, the cosmic ray spatial distribution is highly inhomogeneous, especially at low energies because the particles are trapped between random magnetic mirrors. These results can significantly change the interpretation of synchrotron observations and thus our understanding of the strength and structure of magnetic fields in the Milky Way and nearby spiral galaxies., Comment: 4 pages, 2 figures, Proceedings of Focus Meeting FM8: New Insights in Extragalactic Magnetic Fields, XXXth IAU General Assembly, Vienna, August 29-31, 2018

In this white paper we introduce the IMAGINE Consortium and its scientific background, goals and structure. Our purpose is to coordinate and facilitate the efforts of a diverse group of researchers in the broad areas of the interstellar medium, Galactic magnetic fields and cosmic rays, and our goal is to develop more comprehensive insights into the structures and roles of interstellar magnetic fields and their interactions with cosmic rays. To achieve a higher level of self-consistency, depth and rigour can only be achieved by the coordinated efforts of experts in diverse areas of astrophysics involved in observational, theoretical and numerical work. We present our view of the present status of this topic, identify its key unsolved problems and suggest a strategy that will underpin our work. The backbone of the consortium is the Interstellar MAGnetic field INference Engine, a publicly available Bayesian platform that employs robust statistical methods to explore the multi-dimensional likelihood space using any number of modular inputs. It provides an interpretation and modelling framework that has the power and flexibility to include a variety of observational, theoretical and numerical lines of evidence into a self-consistent and comprehensive picture of the thermal and non-thermal interstellar media. An important innovation is that a consistent understanding of the phenomena that are directly or indirectly influenced by the Galactic magnetic field, such as the deflection of ultra-high energy cosmic rays or extragalactic backgrounds, is made an integral part of the modelling. The IMAGINE Consortium, which is informal by nature and open to new participants, hereby presents a methodological framework for the modelling and understanding of Galactic magnetic fields that is available to all communities whose research relies on a state-of-the-art solution to this problem. (Abridged.), Comment: 56 pages, 15 figures, minor updates to match published JCAP version

Synchrotron radiation from cosmic rays is a key observational probe of the galactic magnetic field. Interpreting synchrotron emission data requires knowledge of the cosmic ray number density, which is often assumed to be in energy equipartition (or otherwise tightly correlated) with the magnetic field energy. However, there is no compelling observational or theoretical reason to expect such tight correlation to hold across all scales. We use test particle simulations, tracing the propagation of charged particles (protons) through a random magnetic field, to study the cosmic ray distribution at scales comparable to the correlation scale of the turbulent flow in the interstellar medium ($\simeq 100\,{\rm pc}$ in spiral galaxies). In these simulations, we find that there is no spatial correlation between the cosmic ray number density and the magnetic field energy density. In fact, their distributions are approximately statistically independent. We find that low-energy cosmic rays can become trapped between magnetic mirrors, whose location depends more on the structure of the field lines than on the field strength., Comment: 14 pages, 15 figures, 1 table, Accepted for publication in MNRAS, Video showing trajectory of a cosmic ray particle in a magnetic trap- https://youtu.be/EzabLlio3FE

The propagation of cosmic rays in turbulent magnetic fields is a diffusive process driven by the scattering of the charged particles by random magnetic fluctuations. Such fields are usually highly intermittent, consisting of intense magnetic filaments and ribbons surrounded by weaker, unstructured fluctuations. Studies of cosmic ray propagation have largely overlooked intermittency, instead relying on Gaussian random magnetic fields. Using test particle simulations, we investigate cosmic ray diffusivity in intermittent, dynamo-generated magnetic fields. The results are compared with those obtained from non-intermittent magnetic fields having identical power spectra. The presence of magnetic intermittency significantly enhances cosmic ray diffusion over a wide range of particle energies. We demonstrate that the results can be interpreted in terms of a correlated random walk., Comment: 5 pages, 3 figures, Accepted for publication in ApJL

Zeldovich's stretch-twist fold (STF) dynamo provided a breakthrough in conceptual understanding of fast dynamos, including fluctuation or small scale dynamos. We study the evolution and saturation behaviour of two types of Baker's map dynamos, which have been used to model Zeldovich's STF dynamo process. Using such maps allows one to analyze dynamos at much higher magnetic Reynolds numbers $R_M$ as compared to direct numerical simulations. In the 2-strip map dynamo there is constant constructive folding while the 4-strip map dynamo also allows the possibility of field reversal. Incorporating a diffusive step parameterised by $R_M$, we find that the magnetic field $B(x)$ is amplified only above a critical $R_M=R_{crit} \sim 4$ for both types of dynamos. We explore the saturation of these dynamos in 3 ways; by a renormalized decrease of the effective $R_M$ (Case I) or due to a decrease in the efficiency of field amplification by stretching (Case II), or a combination of both effects (Case III). For Case I, we show that $B(x)$ in the saturated state, for both types of maps, goes back to the marginal eigenfunction, which is obtained for the critical $R_M=R_{crit}$. This is independent of the initial $R_M=R_{M0}$. On the other hand in Case II, for the 2-strip map, we show that $B(x)$ now saturates preserving the structure of the kinematic eigenfunction. Thus the energy is transferred to larger scales in Case I but remains at the smallest resistive scales in Case II. For the 4-strip map, the $B(x)$ oscillates with time, although with a structure similar to the kinematic eigenfunction. Interestingly, the saturated state for Case III shows an intermediate behaviour, with $B(x)$ now similar to the kinematic eigenfunction for an intermediate $R_M=R_{sat}$, with $R_{M0}>R_{sat}>R_{crit}$. These saturation properties are akin to the ones discussed in the context of fluctuation dynamos., Comment: 17 pages, 14 figures, Accepted in JPP (Special Issue - Zeldovich)

We investigate the anomalies in the Earth - Moon system using ancient eclipse data. We identify nine groups of anomalous eclipses between 400 and 1800 AD recorded in parts of India that should have completely missed the subcontinent as per NASA simulations (Espenak and Meeus, 2011). We show that the typical correction to the lunar location required to reconcile the anomalous eclipses is relatively small and consistent with the fluctuations in the length of day that are observed in recent periods. We then investigate the change in Earth's moment of inertia due to differential acceleration of land and water that can account for this discrepancy. We show that 80 percent of these discrepancies occur when the Moon is at declinations greater than 10 deg and closer to its major standstill of 28 deg while it spends 46 percent of the time in this region. We simulate the differential interaction of the Moon's gravity with landmass and water using finite element method to account for landmass and water mass. We show that the results of eclipse error are consistent with the estimate of a small differential acceleration when the Moon is over land at high latitudes. However, we encounter some examples where the results from simulations studies cannot explain the phenomena. Hence we propose that dT corrections have to be coupled with some other mechanism possibly a small vertical oscillation in the Moon's rotational plane with a period of the order of a few hundred years to achieve the required adjustment in the Eclipse maps., Comment: 11 pages, 2 tables, 15 figures in the text and 2 in appendix