Frederick Armstrong Gent

Magnetic fields are critical at many scales to galactic dynamics and structure, including multiphase pressure balance, dust processing, and star formation. Dynamo action determines their dynamical structure and strength. Simulations of combined large- and small-scale dynamos have successfully developed mean fields with strength and topology consistent with observations, but with turbulent fields much weaker than observed, while simulations of small-scale dynamos with parameters relevant to the interstellar medium yield turbulent fields 1 order of magnitude below the values observed or expected theoretically. We use the Pencil Code accelerated on GPUs with A staroth to perform high-resolution simulations of a supernova-driven galactic dynamo, including heating and cooling in a periodic domain. Our models show that the strength of the turbulent field produced by the small-scale dynamo approaches an asymptote at only modest magnetic Prandtl numbers. This allows us to use these models to suggest the essential characteristics of this constituent of the magnetic field for inclusion in global galactic models. The asymptotic limit occurs already at a magnetic Prandtl number of only a few hundred, many orders of magnitude below physical values in the interstellar medium and consistent with previous findings for isothermal compressible flows.

Magnetic buoyancy (MBI) and Parker instabilities are strong, 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 arise from mean-field dynamo action, particularly the $\alpha ^2\Omega$-dynamo. Using non-ideal MHD equations, we model a section of the galactic disc where the large-scale magnetic field is generated by an imposed $\alpha$-effect and differential rotation. We extend our previous study of the interplay between magnetic buoyancy and the mean-field dynamo by incorporating differential rotation, which enhances the dynamo, and cosmic rays, which amplify magnetic buoyancy. We construct a simple one-dimensional (1D) model which replicates all significant features of the three-dimensional (3D) simulations. Simulations confirm that magnetic buoyancy can lead to oscillatory magnetic fields and reveal that it can change the magnetic field parity between quadrupolar and dipolar states. Differential rotation facilitates this switch in parity, suggesting that the large-scale magnetic field can adopt a dipolar parity within a few kiloparsecs of the galactic centre. In contrast, quadrupolar parity may remain predominant in the outer parts of a galactic disc. Cosmic rays accelerate both the dynamo and the MBI, supporting oscillatory non-linear states and a spatial magnetic field structure similar to the alternating magnetic field directions observed in the haloes of some edge-on galaxies.

Dust in the interstellar medium (ISM) is critical to the absorption and intensity of emission profiles used widely in astronomical observations, and necessary for star and planet formation. Supernovae (SNe) both produce and destroy ISM dust. In particular the destruction rate is difficult to assess. Theory and prior simulations of dust processing by SNe in a uniform ISM predict quite high rates of dust destruction, potentially higher than the supernova dust production rate in some cases. Here we show simulations of supernova-induced dust processing with realistic ISM dynamics including magnetic field effects and demonstrate how ISM inhomogeneity and magnetic fields inhibit dust destruction. Compared to the non-magnetic homogeneous case, the dust mass destroyed within 1 Myr per SNe is reduced by more than a factor of two, which can have a great impact on the ISM dust budget.

Dust in the interstellar medium (ISM) is critical to the absorption and intensity of emission profiles used widely in astronomical observations, and necessary for star and planet formation. Supernovae (SNe) both produce and destroy ISM dust. In particular the destruction rate is difficult to assess. Theory and prior simulations of dust processing by SNe in a uniform ISM predict quite high rates of dust destruction, potentially higher than the supernova dust production rate in some cases. Here we show simulations of supernova-induced dust processing with realistic ISM dynamics including magnetic field effects and demonstrate how ISM inhomogeneity and magnetic fields inhibit dust destruction. Compared to the non-magnetic homogeneous case, the dust mass destroyed within 1 Myr per SNe is reduced by more than a factor of two, which can have a great impact on the ISM dust budget.

Galaxies are observed to host magnetic fields with a typical total strength of around 15microgauss. A coherent large-scale field constitutes up to a few microgauss of the total, while the rest is built from strong magnetic fluctuations over a wide range of spatial scales. This represents sufficient magnetic energy for it to be dynamically significant. Several questions immediately arise: What is the physical mechanism that gives rise to such magnetic fields? How do these magnetic fields affect the formation and evolution of galaxies? In which physical processes do magnetic fields play a role, and how can that role be characterized? Numerical modelling of magnetized flows in galaxies is playing an ever-increasing role in finding those answers. We review major techniques used for these models. Current results strongly support the conclusion that field growth occurs during the formation of the first galaxies on timescales shorter than their accretion timescales due to small-scale turbulent dynamos. The saturated small-scale dynamo maintains field strengths at a few percent of equipartition with turbulence. The subsequent action of large-scale dynamos in differentially rotating discs produces observed modern field strengths in equipartition with the turbulence and having power at large scales. The field structure resulting appears consistent with observations including Faraday rotation and polarisation from synchrotron and dust thermal emission. Major remaining challenges include scaling numerical models toward realistic scale separations and Prandtl and Reynolds numbers.

Magnetic fields are now widely recognized as critical at many scales to galactic dynamics and structure, including multiphase pressure balance, dust processing, and star formation. Using imposed magnetic fields cannot reliably model the interstellar medium's (ISM) dynamical structure nor phase interactions. Dynamos must be modeled. ISM models exist of turbulent magnetic fields using small-scale dynamo (SSD). Others model the large-scale dynamo (LSD) organizing magnetic fields at the scale of the disk or spiral arms. Separately, neither can fully describe the galactic magnetic field dynamics nor topology. We model the LSD and SSD together at a sufficient resolution to use the low explicit Lagrangian resistivity required. The galactic SSD saturates within 20 Myr. We show that the SSD is quite insensitive to the presence of an LSD and is even stronger in the presence of a large-scale shear flow. The LSD grows more slowly in the presence of SSD, saturating after 5 Gyr versus 1–2 Gyr in studies where the SSD is weak or absent. The LSD primarily grows in warm gas in the galactic midplane. Saturation of the LSD occurs due to α-quenching near the midplane as the growing mean-field produces a magnetic α that opposes the kinetic α. The magnetic energy in our models of the LSD shows a slightly sublinear response to increasing resolution, indicating that we are converging toward the physical solution at 1 pc resolution. Clustering supernovae in OB associations increases the growth rates for both the SSD and the LSD, compared to a horizontally uniform supernova distribution.

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 $α$-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.

Magnetic fields on small scales are ubiquitous in the Universe. Although they can often be observed in detail, their generation mechanisms are not fully understood. One possibility is the so-called small-scale dynamo (SSD). Prevailing numerical evidence, however, appears to indicate that an SSD is unlikely to exist at very low magnetic Prandtl numbers (Pr_M) such as those that are present in the Sun and other cool stars. Here we have performed high-resolution simulations of isothermal forced turbulence using the lowest Pr_M values achieved so far. Contrary to earlier findings, the SSD not only turns out to be possible for Pr_M down to 0.0031 but also becomes increasingly easier to excite for Pr_M below about 0.05. We relate this behaviour to the known hydrodynamic phenomenon referred to as the bottleneck effect. Extrapolating our results to solar values of Pr_M indicates that an SSD would be possible under such conditions.