Search

The magnetic buoyancy (MBI) and Parker instabilities are strong and generic instabilities expected to occur in most astrophysical systems with sufficiently strong magnetic fields. In galactic and accretion discs, large-scale magnetic fields are thought to result from the mean-field dynamo action, in particular, the $\alpha^2\Omega$. Using non-ideal MHD equations, we model a section of the galactic disc in which the large-scale magnetic field is generated by an imposed $\alpha$-effect and differential rotation. We extend our earlier study of the interplay between magnetic buoyancy and the mean-field dynamo. We add differential rotation which enhances the dynamo and cosmic rays which enhance magnetic buoyancy. We construct a simple 1D model which replicates all significant features of the 3D simulations. We confirm that magnetic buoyancy can lead to oscillatory magnetic fields and discover that it can vary the magnetic field parity between quadrupolar and dipolar, and that inclusion of the differential rotation is responsible for the switch in field parity. Our results suggest that the large-scale magnetic field can have a dipolar parity within a few kiloparsecs of the galactic centre, provided the MBI is significantly stronger the the dynamo. Quadrupolar parity can remain predominant in the outer parts of a galactic disc. Cosmic rays accelerate both the dynamo and the MBI and support oscillatory non-linear states, a spatial magnetic field structure similar to the alternating magnetic field directions observed in some edge-on galaxies., Comment: 13 pages, 2 tables,, 9 figures

(abridged) Ultra-high energy cosmic rays (UHECRs) are extremely energetic charged particles with energies surpassing $10^{18}$ eV. Their sources remain elusive, obscured by deflections caused by the Galactic magnetic field (GMF). This challenge is further complicated by our limited understanding of the three-dimensional structure of the GMF, as current GMF observations consist primarily of quantities integrated along the line-of-sight (LOS). Nevertheless, data from upcoming stellar polarisation surveys along with Gaia's stellar parallax data are expected to yield local GMF measurements.. In this work, we employ methods of Bayesian statistical inference in order to sample the posterior distribution of the GMF within part of the Galaxy. By assuming a known rigidity and arrival direction of an UHECR, we backtrack its trajectory through various GMF configurations drawn from the posterior distribution. Our objective is to rigorously evaluate our algorithm's performance in scenarios that closely mirror the setting of expected future applications. In pursuit of this, we condition the posterior to synthetic integrated LOS measurements of the GMF, in addition to synthetic local POS-component measurements. In this proof of concept work, we assume the ground truth to be a magnetic field produced by a dynamo simulation of the Galactic ISM. Our results demonstrate that for all locations of the observed arrival direction on the POS, our algorithm is able to substantially update our knowledge on the original arrival direction of UHECRs with rigidity $E/Z = 5 \times 10^{19}$ eV, even in the case of complete absence of LOS information. If integrated data is included in the inference, then the regions of the celestial sphere where the maximum error occurs diminishes greatly. Even in those regions the maximum error is diminished by a factor of about $3$ in the specific setting studied., Comment: 11 pages, 7 figures

Stratified disks with strong horizontal magnetic fields are susceptible to magnetic buoyancy instability (MBI). Modifying the magnetic field and gas distributions, this can play an important role in galactic evolution. The MBI and the Parker instability, in which cosmic rays exacerbate MBI, are often studied using an imposed magnetic field. However, in galaxies and accretion discs, the magnetic field is continuously replenished by a large-scale dynamo action. Using non-ideal MHD equations, we model a section of the galactic disc (we neglect rotation and cosmic rays considered elsewhere), in which the large-scale field is generated by an imposed $\alpha$-effect of variable intensity to explore the interplay between dynamo instability and MBI. The system evolves through three distinct phases: the linear (kinematic) dynamo stage, the onset of linear MBI when the magnetic field becomes sufficiently strong and the nonlinear, statistically steady state. Nonlinear effects associated with the MBI introduce oscillations which do not occur when the field is produced by the dynamo alone. The MBI initially accelerates the magnetic field amplification but the growth is quenched by the vertical motions produced by MBI. We construct a 1D model, which replicates all significant features of 3D simulations to confirm that magnetic buoyancy alone can quench the dynamo and is responsible for the magnetic field oscillations. Unlike with an imposed magnetic field (arXiv:2305.03318,arXiv:2212.03215), the nonlinear interactions do not reduce the gas scale height, so the consequences of the magnetic buoyancy depend on how the magnetic field is maintained.

We describe a novel method to compute the components of dynamo tensors from direct magnetohydrodynamic (MHD) simulations. Our method relies upon an extension and generalisation of the standard H\"ogbom CLEAN algorithm widely used in radio astronomy to systematically remove the impact of the strongest beams onto the corresponding image. This generalisation, called the Iterative Removal of Sources (IROS) method, has been adopted here to model the turbulent electromotive force (EMF) in terms of the mean magnetic fields and currents. Analogous to the CLEAN algorithm, IROS treats the time series of the mean magnetic field and current as beams that convolve with the dynamo coefficients which are treated as (clean) images to produce the EMF time series (the dirty image). We apply this method to MHD simulations of galactic dynamos, to which we have previously employed other methods of computing dynamo coefficients such as the test-field method, the regression method, as well as local and non-local versions of the singular value decomposition (SVD) method. We show that our new method reliably recovers the dynamo coefficients from the MHD simulations. It also allows priors on the dynamo coefficients to be incorporated easily during the inversion, unlike in earlier methods. Moreover, using synthetic data, we demonstrate that it may serve as a viable post-processing tool in determining the dynamo coefficients, even when the power of additive noise to the EMF is twice as much the actual EMF., Comment: 11 pages, 9 figures

We analyse the observational signatures of galactic magnetic fields that are self-consistently generated in magnetohydrodynamic simulations of the interstellar medium through turbulence driven by supernova (SN) explosions and differential rotation. In particular, we study the time evolution of the Faraday rotation measure (RM), synchrotron radiation, and Stokes parameters by characterising the typical structures formed in the plane of observation. We do this by defining two distinct models for both thermal and cosmic ray (CR) electron distributions. Our results indicate that the maps of RM have structures which are sheared and rendered anisotropically by differential rotation and that they depend on the choice of thermal electrons model as well as the SN rate. Synchrotron maps are qualitatively similar to the maps of the mean magnetic field along the line of sight and structures are only marginally affected by the CR model. Stokes parameters and related quantities, such as the degree of linear polarisation, are highly dependent on both frequency and resolution of the observation., Comment: 15 pages, 19 figures

The generation of large-scale magnetic fields ($\overline{\mathbf{B}}$) in astrophysical systems is driven by the mean turbulent electromotive force ($\overline{{\cal E}}$), the cross correlation between local fluctuations of velocity and magnetic fields. This can depend non-locally on $\overline{\mathbf{B}}$ through a convolution kernel $K_{ij}$. In a new approach to find $K_{ij}$, we directly fit the time series data of $\overline{{\cal E}}$ versus $\overline{\mathbf{B}}$ from a galactic dynamo simulation using singular value decomposition. We calculate the usual turbulent transport coefficients as moments of $K_{ij}$, show the importance of including non-locality over eddy length scales to fully capture their amplitudes and that higher order corrections to the standard transport coefficients are small in the present case., Comment: 11 pages, 13 figures

Large-scale coherent magnetic fields observed in the nearby galaxies are thought to originate by a mean-field dynamo. This is governed via the turbulent electromotive force (EMF, $\overline{\mathcal{E}} $) generated by the helical turbulence driven by supernova (SN) explosions in the differentially rotating interstellar medium (ISM). In this paper we aim to investigate the possibility of dynamo action by the virtue of buoyancy due to a cosmic ray (CR) component injected through the SN explosions. We do this by analysing the magnetohydrodynamic simulations of local shearing box of ISM, in which the turbulence is driven via random SN explosions and the energy of the explosion is distributed in the CR and/or thermal energy components. We use the magnetic field aligned diffusion prescription for the propagation of CR. We compare the evolution of magnetic fields in the models with the CR component to our previous models that did not involve the CR. We demonstrate that the inclusion of CR component enhances the growth of dynamo slightly. We further compute the underlying dynamo coefficients using the test-fields method, and argue that the entire evolution of the large scale mean magnetic field can be reproduced with an $\alpha-\Omega$ dynamo model. We also show that the inclusion of CR component leads to an unbalanced turbulent pumping between magnetic field components and additional dynamo action by the R\"adler effect., Comment: 9 pages, 5 figures, 1 tables

Coherent magnetic fields in disc galaxies are thought to be generated by a large-scale (or mean-field) dynamo operating in their interstellar medium. A key driver of mean magnetic field growth is the turbulent electromotive force (EMF), which represents the influence of correlated small-scale (or fluctuating) velocity and magnetic fields on the mean field. The EMF is usually expressed as a linear expansion in the mean magnetic field and its derivatives, with the dynamo tensors as expansion coefficients. Here, we adopt the singular value decomposition (SVD) method to directly measure these turbulent transport coefficients in a simulation of the turbulent interstellar medium that realizes a large-scale dynamo. Specifically, the SVD is used to least-square fit the time series data of the EMF with that of the mean field and its derivatives, to determine these coefficients. We demonstrate that the spatial profiles of the EMF reconstructed from the SVD coefficients match well with that taken directly from the simulation. Also, as a direct test, we use the coefficients to simulate a 1-D mean-field dynamo model and find an overall similarity in the evolution of the mean magnetic field between the dynamo model and the direct simulation. We also compare the results with those which arise using simple regression and the ones obtained previously using the test-field (TF) method, to find reasonable qualitative agreement. Overall, the SVD method provides an effective post-processing tool to determine turbulent transport coefficients from simulations., Comment: 14 pages, 14 figures