Small-scale Dynamo Research Articles

Solar internetwork magnetic fields: Statistical comparison between observations and MHD simulations

Although the magnetic fields in the quiet Sun account for the majority of the magnetic energy in the solar photosphere, inferring their exact spatial distribution, origin, and evolution poses an important challenge because the signals lie at the limit of today's instrumental precision. This severely hinders and biases our interpretations, which are mostly made through nonlinear model-fitting approaches. Our goal is to directly compare simulated and observed polarization signals in the Fe I $ angstrom $ and $ angstrom $ spectral lines in the very quiet Sun, the so-called solar internetwork (IN). This way, we aim to constrain the mechanism responsible for the generation of the quiet Sun magnetism while avoiding the biases that plague other diagnostic methods. We used three different three-dimensional radiative magneto-hydrodynamic simulations representing different scenarios of magnetic field generation in the internetwork: small-scale dynamo, decay of active regions, and horizontal flux emergence. We synthesized Stokes profiles at different viewing angles and degraded them according to the instrumental specifications of the spectro-polarimeter (SP) on board the Hinode satellite. Finally, we statistically compared the simulated spectra to the Hinode/SOT/SP observations at the appropriate viewing angles. Of the three simulations, the small-scale dynamo best reproduced the statistical properties of the observed polarization signals. This is especially prominent for the disk center viewing geometry, where the agreement is excellent. Moving toward more inclined lines of sight, the agreement worsens slightly. The agreement between the small-scale dynamo simulation and observations at the disk center suggests that small-scale dynamo action plays an important role in the generation of quiet Sun magnetism. However, the magnetic field around 50 km above the continuum layer in this simulation does not reproduce observations as well as at the very base of the photosphere.

Large-scale ordered magnetic fields generated in mergers of helium white dwarfs

Stellar mergers are one important path to highly magnetised stars. Mergers of two low-mass white dwarfs may create up to every third hot subdwarf star. The merging process is usually assumed to dramatically amplify magnetic fields. However, so far only four highly magnetised hot subdwarf stars have been found, suggesting a fraction of less than $1<!PCT!>$. We present two high-resolution magnetohydrodynamical (MHD) simulations of the merger of two helium white dwarfs in a binary system with the same total mass of $0.6\,M_ We analysed an equal-mass merger with two $0.3\,M_ white dwarfs, and an unequal-mass merger with white dwarfs of $0.25\,M_ and $0.35\,M_ We simulated the inspiral, merger, and further evolution of the merger remnant for about $50$ rotations. We found efficient magnetic field amplification in both mergers via a small-scale dynamo, reproducing previous results of stellar merger simulations. The magnetic field saturates at a similar strength for both simulations. We then identified a second phase of magnetic field amplification in both merger remnants that happens on a timescale of several tens of rotational periods of the merger remnant. This phase generates a large-scale ordered azimuthal field via a large-scale dynamo driven by the magneto-rotational instability. Finally, we speculate that in the unequal-mass merger remnant, helium burning will initially start in a shell around a cold core, rather than in the centre. This forms a convection zone that coincides with the region that contains most of the magnetic energy, and likely destroys the strong, ordered field. Ohmic resistivity might then quickly erase the remaining small-scale field. Therefore, the mass ratio of the initial merger could be the selecting factor that decides if a merger remnant will stay highly magnetised long after the merger.

Vorticity and magnetic dynamo from subsonic expansion waves. II. Dependence on the magnetic Prandtl number, forcing scale, and cooling time

The amplification of astrophysical magnetic fields takes place via dynamo instability in turbulent environments. Vorticity is usually present in any dynamo, but its role is not yet fully understood. This work is an extension of previous research on the effect of an irrotational subsonic forcing on a magnetized medium in the presence of rotation or a differential velocity profile. We aim to explore a wider parameter space in terms of Reynolds numbers, the magnetic Prandtl number, the forcing scale, and the cooling timescale in a Newtonian cooling. We studied the effect of imposing that either the acceleration or the velocity forcing function be curl-free and evaluated the terms responsible for the evolution vorticity. We used direct numerical simulations to solve the fully compressible, resistive magnetohydrodynamic equations with the Pencil Code. We studied both isothermal and non-isothermal regimes and addressed the relative importance of different vorticity source terms. We report no small-scale dynamo for the models that do not include shear. We find a hydro instability, followed by a magnetic one, when a shearing velocity profile is applied. The vorticity production is found to be numerical in the purely irrotational acceleration case. Non-isothermality, rotation, shear, and density-dependent forcing, when included, contribute to increasing the vorticity. As in our previous study, we find that turbulence driven by subsonic expansion waves can amplify the vorticity and magnetic field only in the presence of a background shearing profile. The presence of a cooling function makes the instability occur on a shorter timescale. We estimate critical Reynolds and magnetic Reynolds numbers of 40 and 20, respectively.

Hysteresis Near the Transition of the Large-Scale Dynamo in the Presence of the Small-Scale Dynamo

Hysteresis Near the Transition of the Large-Scale Dynamo in the Presence of the Small-Scale Dynamo

Small-scale Dynamo with Nonzero Correlation Time

The small-scale dynamo is typically studied by assuming that the correlation time of the velocity field is zero. Some authors have used a smooth renovating flow model to study how the properties of the dynamo are affected by the correlation time being nonzero. Here, we assume the velocity is an incompressible Gaussian random field (which need not be smooth), and derive the lowest-order corrections to the evolution equation for the two-point correlation of the magnetic field in Fourier space. Using this, we obtain the evolution equation for the longitudinal correlation function of the magnetic field (M L ) in nonhelical turbulence, valid for arbitrary Prandtl number. The nonresistive terms of this equation do not contain spatial derivatives of M L of order greater than 2. We further simplify this equation in the limit of high Prandtl number, and find that the growth rate of the magnetic energy is much smaller than previously reported. Nevertheless, the magnetic power spectrum still retains the Kazantsev form at high Prandtl number.

Thermal conductivity with bells and whistlers: Suppression of the magnetothermal instability in galaxy clusters

In the hot and dilute intracluster medium (ICM) in galaxy clusters, kinetic plasma instabilities that are excited at the particle gyroradius may play an important role in the transport of heat and momentum, thus affecting the large-scale evolution of these systems. In this paper, we continue our investigation of the effect of whistler suppression of thermal conductivity on the magneto-thermal instability (MTI), which may be active in the periphery of galaxy clusters and may contribute to the observed levels of turbulence. We use a subgrid closure for the heat flux inspired from kinetic simulations and show that MTI turbulence with whistler suppression exhibits a critical transition as the suppression parameter is increased: for modest suppression of the conductivity, the turbulent velocities generated by the MTI decrease accordingly, in agreement with scaling laws found in previous studies of the MTI. However, for suppression above a critical threshold, the MTI loses its ability to maintain equipartition-level magnetic fields through a small-scale dynamo (SSD), and the system enters a ``death-spiral.'' We show that analogous levels of suppression of thermal conductivity with a simple model of flat uniform suppression would not inhibit the dynamo. We propose a model to explain this critical transition, and speculate that conditions in the hot ICM are such that in substantial portions of the galaxy cluster periphery the MTI might struggle to sustain its own dynamo. We then look at spatial correlations and energy transfers in spectral space and find that, with whistler suppression, most of the heat is transported along thin bundles of strong magnetic fields (the Autobahns of electrons), while high-beta regions are brought out of thermal equilibrium. We link this behavior to the intermittent nature of magnetic fields, and we observe an overall reduction of the efficiency of MTI turbulent driving at the largest turbulent scales. Finally, we show that external turbulence interferes with the MTI and leads to reduced levels of MTI turbulence. While individually both external turbulence and whistler suppression reduce MTI turbulence, we find that they can exhibit a complex interplay when acting in conjunction, with external turbulence boosting the whistler-suppressed thermal conductivity and even reviving a ``dead'' MTI. Our study illustrates how extending magnetohydrodynamics with a simple prescription for microscale plasma physics can lead to the formation of a complicated dynamical system and demonstrates that further work is needed in order to bridge the gap between micro- and macro scales in galaxy clusters.

Toward Cosmological Simulations of the Magnetized Intracluster Medium with Resolved Coulomb Collision Scale

We present the first results of one extremely high-resolution, nonradiative magnetohydrodynamical cosmological zoom-in simulation of a massive cluster with a virial mass of M vir = 2.0 × 1015 solar masses. We adopt a mass resolution of 4 × 105 M ⊙ with a maximum spatial resolution of around 250 pc in the central regions of the cluster. We follow the detailed amplification process in a resolved small-scale turbulent dynamo in the intracluster medium (ICM) with strong exponential growth until redshift 4, after which the field grows weakly in the adiabatic compression limit until redshift 2. The energy in the field is slightly reduced as the system approaches redshift zero in agreement with adiabatic decompression. The field structure is highly turbulent in the center and shows field reversals on a length scale of a few tens of kiloparsecs and an anticorrelation between the radial and angular field components in the central region that is ordered by small-scale turbulent dynamo action. The large-scale field on megaparsec scales is almost isotropic, indicating that the structure formation process in massive galaxy cluster formation suppresses any memory of both the initial field configuration and the amplified morphology via the turbulent dynamo. We demonstrate that extremely high-resolution simulations of the magnetized ICM are within reach that can simultaneously resolve the small-scale magnetic field structure, which is of major importance for the injection of and transport of cosmic rays in the ICM. This work is a major cornerstone for follow-up studies with an on-the-fly treatment of cosmic rays to model in detail electron-synchrotron and gamma-ray emissions.

Magnetic field amplification in massive primordial halos

Context. The potential importance of magnetic fields during structure formation and gravitational collapse in the early Universe has been shown in several studies. In particular, magnetic field amplification by the small-scale dynamo plays an important role in addition to the pure amplification expected from gravitational collapse. Aims. In this paper we study the small-scale dynamo for halos of ≳107 M⊙ collapsing at z ≳ 12, under different ambient conditions due to the strength of the Lyman-Werner background. Additionally, we estimate the approximate saturation level by varying the initial magnetic field strength. Methods. We performed cosmological magnetohydrodynamical simulations for three distinct halos of ∼107 M⊙ at z ≥ 12 by varying the Jeans resolution from 32 − 256 cells and employed Lyman Werner background flux of strengths 102 − 105 in units of J21, where J21 = 10−21 erg cm−2 sr−1 s−1Hz−1. To follow the chemical and thermal evolution of the gas, we made use of the KROME package. Results. In addition to the compression by collapse, we find magnetic field amplification via the dynamo in the regimes of atomic and molecular hydrogen cooling. Moreover, we find a lower saturation level in the molecular hydrogen cooling regime. This behaviour can be understood in terms of the generally reduced radial infall velocities and vorticities in this regime, as well as the higher Mach numbers of the gas, which give rise to a smaller saturation ratio. Conclusions. Our results overall suggest that the dynamo operates over a large range of conditions in the collapsing gas.

Simulating small-scale dynamo action in cool main-sequence stars

Context. The origin of the ubiquitous small-scale magnetic field observed on the solar surface can be attributed to the presence of a small-scale dynamo (SSD) operating in the sub-surface layers of the Sun. It is expected that a similar process could self-sustain a considerable amount of magnetic energy also in the near-surface layers of cool main-sequence stars other than the Sun. Aims. In this paper the properties of the magnetic field resulting from SSD action operating in the near-surface layers of four cool main-sequence stars and its self-organization into magnetic flux concentrations are investigated numerically. Methods. Three-dimensional radiative magnetohydrodynamic simulations of SSD action in the near-surface layers of four cool main-sequence stars of spectral types K8V, K2V, G2V, and F5V are carried out with the CO5BOLD code. The simulations are set up to have approximately the same Reynolds and magnetic Reynolds numbers, and to disentangle the impact of the effective temperature and the surface gravity on the SSD action from numerical effects. Results. It is found that the SSD growth rates in SI units differ for the four stellar models; the highest and lowest growth rate is for the K2V and F5V model, respectively. This is due to the different turnover times in the four simulations. Even so, the SSD field strengths reached in the saturation phases are similar in all models, with the same amount of kinetic energy converted into magnetic energy. If the magnetic energy that is pumped out from the computational domain across the bottom boundary is partially replenished from outside of the computational domain, we find that the SSD action leads to a sufficient reduction in the convective velocities to reduce the convective horizontal length scales in the convection zone by 5–10%, vanishing towards the optical depth unity level. In this case, strong kilogauss magnetic flux concentrations emerge at the surface, leading to magnetic bright features, which are more numerous and conspicuous for the K2V and G2V models than for the K8V and F5V models. Their vertical magnetic field component on the surface of optical depth unity increases from 1 kG to 1.6 kG with decreasing effective temperature from F5V to K8V. However, more than 90% of the magnetic flux through any of these stellar surfaces has a field strength of less than 1 kG.

The effect of pressure-anisotropy-driven kinetic instabilities on magnetic field amplification in galaxy clusters

Context. The intracluster medium (ICM) is the low-density diffuse gas that fills the space between galaxies within galaxy clusters. It is primarily composed of magnetized plasma, which reaches virial temperatures of up to 108 K, probably due to mergers of subhalos. Under these conditions, the plasma is weakly collisional and therefore has an anisotropic pressure tensor with respect to the local direction of the magnetic field. This triggers very fast, Larmor-scale, pressure-anisotropy-driven kinetic instabilities that alter magnetic field amplification. Aims. We aim to study magnetic field amplification through a turbulent, small-scale dynamo, including the effects of the kinetic instabilities, during the evolution of a typical massive galaxy cluster. A specific aim of this work is to establish a redshift limit from which a dynamo has to start to amplify the magnetic field up to equipartition with the turbulent velocity field at redshift z = 0. Methods. We implemented one-dimensional radial profiles for various plasma quantities for merger trees generated with the modified GALFORM algorithm. We assumed that turbulence is driven by successive mergers of dark matter halos and constructed effective models for the Reynolds number Reeff dependence on the magnetic field in three different magnetization regimes (unmagnetized, magnetized “kinetic”, and magnetized “fluid”), including the effects of kinetic instabilities. The magnetic field growth rate is calculated for the different Reeff models. Results. The model results in a higher magnetic field growth rate at higher redshift. For all scenarios considered in this study, to reach equipartition at z = 0, it is sufficient for the amplification of the magnetic field to start at redshift zstart ≈ 1.5 and above. The time to reach equipartition can be significantly shorter in cases with systematically smaller turbulent forcing scales and for the highest Reeff models. Conclusions. The origin of magnetic fields in the weakly collisional ICM can be explained by the small-scale turbulent dynamo, provided that the dynamo process starts beyond a given redshift. Merger trees are useful tools for studying the evolution of magnetic fields in weakly collisional plasmas, and could also be used to constrain the different stages of the dynamo that could potentially be observed by future radio telescopes.

A Floor in the Sun's Photospheric Magnetic Field: Implications for an Independent Small-scale Dynamo

Clette recently showed that F 10.7 systematically approaches a quiet Sun daily value of 67 solar flux units (sfu) at solar minima as the number of spotless days on the Sun increases. Previously, a floor of ∼2.8 nT had been proposed for the solar wind (SW) magnetic field strength (B). F 10.7, which closely tracks the Sun's unsigned photospheric magnetic flux, and SW B exhibit different relationships to their floors at 11 yr solar minima during the last ∼50 yr. While F 10.7 approaches 67 sfu at each minimum, the corresponding SW B is offset above ∼2.8 nT by an amount approximately proportional to the solar polar field strength—which varied by a factor of ∼2.5 during this interval. This difference is substantiated by ∼130 yr of reconstructed F 10.7 (via the range of the diurnal variation of the East-component (rY) of the geomagnetic field) and SW B (based on the interdiurnal variability geomagnetic activity index). For the last ∼60 yr, the contribution of the slow SW to SW B has exhibited a floor-like behavior at ∼2 nT, in contrast to the contributions of coronal mass ejections and high-speed streams that vary with the solar cycle. These observations, as well as recent SW studies based on Parker Solar Probe and Solar Dynamics Observatory data, suggest that (1) the Sun has a small-scale turbulent dynamo that is independent of the 11 yr sunspot cycle; and (2) the small-scale magnetic fields generated by this nonvarying turbulent dynamo maintain a constant open flux carried to the heliosphere by the Sun's floor-like slow SW.

Transition from Small-scale to Large-scale Dynamo in a Supernova-driven, Multiphase Medium

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.

Small-scale dynamo in cool stars

Context. Some of the quiet solar magnetic flux could be attributed to a small-scale dynamo (SSD) operating in the convection zone. An SSD operating in cool main-sequence stars is expected to affect the atmospheric structure, in particular, the convection, and should have observational signatures. Aims. We investigate the distribution of SSD magnetic fields and their effect on bolometric intensity characteristics, vertical velocity, and spatial distribution of the kinetic energy (KE) and magnetic energy (ME) in the lower photosphere of different spectral types. Methods. We analyzed the SSD and purely hydrodynamic simulations of the near surface layers of F3V, G2V, K0V, and M0V stars. We compared the time-averaged distributions and power spectra in SSD setups relative to the hydrodynamic setup. The properties of the individual magnetic fields are also considered. Results. The probability density functions with a field strength at the τ = 1 surface are quite similar for all cases. The M0V star displays the strongest fields, but relative to the gas pressure, the fields on the F3V star reach the highest values. In all stars, the horizontal field is stronger than the vertical field in the middle photosphere, and this excess becomes increasingly prominent toward later spectral types. These fields result in a decrease in the upflow velocities and a slight decrease in granule size, and also lead to formation of bright points in intergranular lanes. The spatial distribution of the KE and ME is also similar for all cases, implying that important scales are proportional to the pressure scale height. Conclusions. The SSD fields have rather similar effects on the photospheres of cool main-sequence stars: a significant reduction in convective velocities, as well as a slight reduction in granule size and a concentration of the field to kilogauss levels in intergranular lanes that is associated with the formation of bright points. The distribution of the field strengths and energies is also rather similar.

ANOMALOUS MAGNETIC REGIONS ON THE SUN

We studied the anomalous magnetic regions observed near the minima of solar cycles 24 and 25. The peculiarity of these areas was the deviation of their configuration from Hale's law of magnetic polarity and Joy's law about the inclination of the axes of bipolar groups to the latitudinal direction. Therefore, they belong to the class of so-called anti-Hale active regions. We paid special attention to the flare activity of anti-Hale regions, as this is important for forecasting space weather and magnetic storms in the Earth's atmosphere. The detected anomalies of the surface magnetism of the active regions studied by us may indicate the influence of the mechanisms of the deep small-scale dynamo on their evolution. In this regard we analyzed the possible mechanisms of the formation of anti-Hale magnetic regions. In particular, such mechanisms can be the mechanisms of a small-scale magnetic dynamo. In connection with this an urgent problem today is the search for observed evidence of the existence of the theoretically proposed by Brandenburg A. et al. (2012) of a new physical entity – a small-scale magnetic field hidden in the solar depths, excited by two qualitatively different mechanisms of a small-scale dynamo (SSD). The first mechanism is the SSD of macroscopic MHD (SSD1), while the second is the diffusion SSD of classical MHD (SSD2). However, the small contributions of these sources are very difficult to distinguish observationally. To solve this complication, Sokoloff, Khlystova and Abramenko (2015) proposed a test for separating the contributions of two sources based on a statistical probabilistic model. Such an important feature of the differences between of the two SSD is the behavior of the percentage of anti-Hail groups of sunspots (in relation to the total number of spots) in the minima of solar cycles. According to statistical studies of long series of observations Sokoloff, Khlystova and Abramenko (2015) found that the percentage of anti-Haile groups of spots increases during minima of the solar cycles, suggesting in favor of SSD2. We believe that the detected magnetic anomalies of the studied regions may be caused by the influence of a SSD2 in the depths of the convective zone of the Sun, since this source gives the most noticeable contribution to the surface magnetism near cycle minima.

Turbulent dynamo action and its effects on the mixing at the convective boundary of an idealized oxygen-burning shell

Convection is one of the most important mixing processes in stellar interiors. Hydrodynamic mass entrainment can bring fresh fuel from neighboring stable layers into a convection zone, modifying the structure and evolution of the star. Because flows in stellar convection zones are highly turbulent, multidimensional hydrodynamic simulations are fundamental to accurately capture the physics of mixing processes. Under some conditions, strong magnetic fields can be sustained by the action of a turbulent dynamo, adding another layer of complexity and possibly altering the dynamics in the convection zone and at its boundaries. In this study, we used our fully compressible SEVEN-LEAGUE HYDRO code to run detailed and highly resolved three-dimensional magnetohydrodynamic simulations of turbulent convection, dynamo amplification, and convective boundary mixing in a simplified setup whose stratification is similar to that of an oxygen-burning shell in a star with an initial mass of 25 M⊙. We find that the random stretching of magnetic field lines by fluid motions in the inertial range of the turbulent spectrum (i.e., a small-scale dynamo) naturally amplifies the seed field by several orders of magnitude in a few convective turnover timescales. During the subsequent saturated regime, the magnetic-to-kinetic energy ratio inside the convective shell reaches values as high as 0.33, and the average magnetic field strength is ∼1010 G. Such strong fields efficiently suppress shear instabilities, which feed the turbulent cascade of kinetic energy, on a wide range of spatial scales. The resulting convective flows are characterized by thread-like structures that extend over a large fraction of the convective shell. The reduced flow speeds and the presence of magnetic fields with strengths up to 60% of the equipartition value at the upper convective boundary diminish the rate of mass entrainment from the stable layer by ≈20% as compared to the purely hydrodynamic case.

Effect of flow shear on the onset of dynamos

Understanding the origin and structure of mean magnetic fields in astrophysical conditions is a major challenge. Shear flows often coexist in such astrophysical conditions, and the role of flow shear on the dynamo mechanism is of great interest. Here, we present a direct numerical simulation study of the effect of flow shear on dynamo instability for EPI2D flows [Yoshida et al., Phys. Rev. Lett. 119, 244501 (2017)] with controllable mirror symmetry (i.e., fluid helicity). Our numerical observations suggest that for helical base flows, the effect of shear is to reduce the small-scale dynamo (SSD) growth rate moderately. For non-helical base flows, flow shear has the opposite effect of amplifying the SSD action. The magnetic energy growth rate (γ) for non-helical base flows has been found to follow an algebraic nature of the form, γ=−aS+bS23, where a,b>0 are real constants, S is the shear flow strength, and γ is found to be independent of the scale of flow shear. Studies with different shear profiles and shear scale lengths for non-helical base flows have been performed to test the universality of our finding.

Azimuthal anisotropy of magnetic fields in the circumgalactic medium driven by galactic feedback processes

ABSTRACT We use the TNG50 cosmological magnetohydrodynamical simulation of the IllustrisTNG project to show that magnetic fields in the circumgalactic medium (CGM) have significant angular structure. This azimuthal anisotropy at fixed distance is driven by galactic feedback processes that launch strong outflows into the halo, preferentially along the minor axes of galaxies. These feedback-driven outflows entrain strong magnetic fields from the interstellar medium, dragging fields originally amplified by small-scale dynamos into the CGM. At the virial radius, z = 0 galaxies with M⋆ ∼ $10^{10}\, \rm {M_\odot }$ show the strongest anisotropy (∼0.35 dex). This signal weakens with decreasing impact parameter, and is also present but weaker for lower mass as well as higher mass galaxies. Creating mock Faraday rotation measure (RM) sightlines through the simulated volume, we find that the angular RM trend is qualitatively consistent with recent observational measurements. We show that rich structure is present in the circumgalactic magnetic fields of galaxies. However, TNG50 predicts small RM amplitudes in the CGM that make detection difficult as a result of other contributions along the line of sight.

Deciphering the radio–star formation correlation on kpc scales

The relation between the resolved star formation rate (SFR) per unit area and the nonthermal radio continuum emission is studied in 21 Virgo cluster galaxies and the two nearby spiral galaxies, NGC 6946 and M 51. For the interpretation and understanding of our results, we used a 3D model where star formation, 2D cosmic-ray (CR) propagation, and the physics of synchrotron emission are included. Based on the linear correlation between the SFR per unit area and the synchrotron emission and its scatter, radio-bright and radio-dim regions can be robustly defined for our sample of spiral galaxies. We identified CR diffusion or streaming as the physical causes of radio-bright regions of unperturbed symmetric spiral galaxies as NGC 6946. The enhanced magnetic field in the region of interstellar medium (ISM) compression via ram pressure is responsible for the southwestern radio-bright region in NGC 4501. We identified the probable causes of radio-bright regions in several galaxies as CR transport, via either gravitational tides (M 51) or galactic winds (NGC 4532) or ram pressure stripping (NGC 4330 and NGC 4522). Three galaxies are overall radio dim: NGC 4298, NGC 4535, and NGC 4567. Based on our model of synchrotron-emitting disks, we suggest that the overall radio-dim galaxies have a significantly lower magnetic field than expected by equipartition between the magnetic and turbulent energy densities. We suggest that this is linked to the difference between the timescales of the variation in the SFR and the small-scale dynamo. In NGC 4535, shear motions increase the total magnetic field strength via the induction equation, which leads to enhanced synchrotron emission with respect to the SFR in an otherwise radio-dim galactic disk. Radio-bright regions frequently coincide with asymmetric ridges of polarized radio continuum emission, and we found a clear albeit moderate correlation between the polarized radio continuum emission and the radio/SFR ratio. When compression or shear motions of the ISM are present in the galactic disk, the radio-bright regions are linked to the commonly observed asymmetric ridges of polarized radio continuum emission and represent a useful tool for the interaction diagnostics. The magnetic field is enhanced (as observed in NGC 4535 and NGC 4501) and ordered by these ISM compression and shear motions. Whereas the enhancement of the magnetic field is rather modest and does not significantly influence the radio-SFR correlation, the main effect of ISM compression and shear motions is the ordering of the magnetic field, which significantly affects the CR transport. Cosmic-ray energy losses and transport also affect the spectral index, which we measured between 4.85 and 1.4 GHz. The influence of CR losses and transport on the spectral index distribution with respect to the synchrotron/SFR ratio is discussed with the help of model calculations. Based on our results, we propose a scenario for the interplay between star formation, CR electrons, and magnetic fields in spiral galaxies.

Vorticity and magnetic dynamo from subsonic expansion waves

Context. The main driving forces supplying energy to the interstellar medium (ISM) are supernova explosions and stellar winds. Such localized sources are assimilable to curl-free velocity fields as a first approximation. They need to be combined with other physical processes to replicate real galactic environments, such as the presence of turbulence and a dynamo-sustained magnetic field in the ISM. Aims. This work is focused on the effect of an irrotational forcing on a magnetized flow in the presence of rotation, baroclinicity, shear, or a combination of any of the three. It follows an earlier analysis with a similar focus, namely, subsonic spherical expansion waves in hydrodynamic simulations. By including magnetic field in the model, we can evaluate the occurrence of dynamo on both small and large scales. We aim to identify the minimum ingredients needed to trigger a dynamo instability as well as the relation between dynamo and the growth of vorticity. Methods. We used the Pencil code to run resistive magnetohydrodynamic direct numerical simulations, exploring the ranges of values of several physical and numerical parameters of interest. We explored Reynolds numbers up to a few hundreds. We analyzed the temporal evolution of vorticity, kinetic, and magnetic energy, as well as their features in Fourier space. Results. We report the absence of a small-scale dynamo in all cases where only rotation is included, regardless of the given equation of state and rotation rate. Conversely, the inclusion of a background sinusoidal shearing profile leads to an hydrodynamic instability that produces an exponential growth of the vorticity at all scales, starting from small ones. This is know as vorticity dynamo. The onset of this instability occurs after a rather long temporal evolution of several thousand turbulent turnover times. The vorticity dynamo in turn drives an exponential growth of the magnetic field, first at small scales, followed by large ones. The instability is then saturated and the magnetic field approximately reaches equipartition with the turbulent kinetic energy. During the saturation phase, we can observe a winding of the magnetic field in the direction of the shearing flow. By varying the intensity of the shear, we see that the growth rates of this instability change. The inclusion of the baroclinic term has the main effect of delaying the onset of the vorticity dynamo, but then leads to a more rapid growth. Conclusions. Our work demonstrates how even purely irrotational forcing may lead to dynamo action in the presence of shear, thus amplifying the field to an equipartition level. At the same time, we confirm that purely irrotational forcing alone does not lead to any growth in terms of the vorticity, nor the magnetic field. This picture does not change in the presence of rotation or baroclinicity, at least up to a resolution of 2563 mesh points. To further generalize such a conclusion, we will need to explore how this setup works both at higher magnetic Reynolds numbers and with different prescriptions of the irrotational forcing.

Galactic Dynamos

Spiral galaxies, including the Milky Way, have large-scale magnetic fields with significant energy densities. The dominant theory attributes these magnetic fields to a large-scale dynamo. We review the current status of dynamo theory and discuss various numerical simulations designed either to explain particular aspects of the problem or to reproduce galactic magnetic fields globally. Our main conclusions can be summarized as follows: ▪Idealized direct numerical simulations produce mean magnetic fields, whose saturation energy density tends to decline with increasing magnetic Reynolds number. This is still an unsolved problem.▪Large-scale galactic magnetic fields of microgauss strengths can probably be explained only if helical magnetic fields of small or moderate length scales can be rapidly ejected or destroyed.▪Small-scale dynamos are important throughout a galaxy's life and probably provide strong seed fields at early stages.▪The circumgalactic medium (CGM) may play an important role in driving dynamo action at small and large length scales. These interactions between the galactic disk and the CGM may provide important insights into our understanding of galactic dynamos. We expect future research in galactic dynamos to focus on the cosmological history of galaxies and the interaction with the CGM as means of replacing the idealized boundary conditions used in earlier work.