Showing 1–50 of 416 results for author: Brandenburg, A

Intergalactic medium rotation measure of primordial magnetic fields

Authors: Salome Mtchedlidze, Paola Domínguez-Fernández, Xiaolong Du, Ettore Carretti, Franco Vazza, Shane Patrick O'Sullivan, Axel Brandenburg, Tina Kahniashvili

Abstract: The Faraday rotation effect, quantified by the Rotation Measure (RM), is a powerful probe of the large-scale magnetization of the Universe - tracing magnetic fields not only on galaxy and galaxy cluster scales but also in the intergalactic Medium (IGM; referred to as $\mathrm{RM}_{\text{IGM}}$). The redshift dependence of the latter has extensively been explored with observations. It has also been… ▽ More The Faraday rotation effect, quantified by the Rotation Measure (RM), is a powerful probe of the large-scale magnetization of the Universe - tracing magnetic fields not only on galaxy and galaxy cluster scales but also in the intergalactic Medium (IGM; referred to as $\mathrm{RM}_{\text{IGM}}$). The redshift dependence of the latter has extensively been explored with observations. It has also been shown that this relation can help to distinguish between different large-scale magnetization scenarios. We study the evolution of this $\mathrm{RM}_{\text{IGM}}$ for different primordial magnetogenesis scenarios to search for the imprints of primordial magnetic fields (PMFs; magnetic fields originating in the early Universe) on the redshift-dependence of $\mathrm{RM}_{\text{IGM}}$. We use cosmological magnetohydrodynamic (MHD) simulations for evolving PMFs during large-scale structure formation, coupled to the light cone analysis to produce a realistic statistical sample of mock $\mathrm{RM}_{\text{IGM}}$ images. We study the predicted behavior for the cosmic evolution of $\mathrm{RM}_{\text{IGM}}$ for different correlation lengths of PMFs, and provide fitting functions for their dependence on redshifts. We compare these mock RM trends with the recent analysis of the the LOw-Frequency ARray (LOFAR) RM Grid and find that large-scale-correlated PMFs should have (comoving) strengths $\lesssim 0.75$ nanoGauss, if originated during inflation with the scale invariant spectrum and (comoving) correlation length $\sim 19$ cMpc/h or $ \lesssim 30$ nanoGauss if they originated during phase-transition epochs with the comoving correlation length $\sim 1$ cMpc/h. Our findings agree with previous observations and confirm the results of semi-analytical studies, showing that upper limits on the PMF strength decrease as their coherence scales increase. △ Less

Submitted 23 June, 2024; originally announced June 2024.

Comments: 16 pages, 9 figures, comments welcome!

Report number: NORDITA-2024-023

Efficiency of dynamos from autonomous generation of chiral asymmetry

Authors: Jennifer Schober, Igor Rogachevskii, Axel Brandenburg

Abstract: At high energies, the dynamics of a plasma with charged fermions can be described in terms of chiral magnetohydrodynamics. Using direct numerical simulations, we demonstrate that chiral magnetic waves (CMWs) can produce a chiral asymmetry $μ_5 = μ_\mathrm{L} - μ_\mathrm{R}$ from a spatially fluctuating (inhomogeneous) chemical potential $μ= μ_\mathrm{L} + μ_\mathrm{R}$, where $μ_\mathrm{L}$ and… ▽ More At high energies, the dynamics of a plasma with charged fermions can be described in terms of chiral magnetohydrodynamics. Using direct numerical simulations, we demonstrate that chiral magnetic waves (CMWs) can produce a chiral asymmetry $μ_5 = μ_\mathrm{L} - μ_\mathrm{R}$ from a spatially fluctuating (inhomogeneous) chemical potential $μ= μ_\mathrm{L} + μ_\mathrm{R}$, where $μ_\mathrm{L}$ and $μ_\mathrm{R}$ are the chemical potentials of left- and right-handed electrically charged fermions, respectively. If the frequency of the CMW is less than or comparable to the characteristic growth rate of the chiral dynamo instability, the magnetic field can be amplified on small spatial scales. The growth rate of this small-scale chiral dynamo instability is determined by the spatial maximum value of $μ_5$ fluctuations. Therefore, the magnetic field amplification occurs during periods when $μ_5$ reaches temporal maxima during the CMW. If the small-scale chiral dynamo instability leads to a magnetic field strength that exceeds a critical value, which depends on the resistivity and the initial value of $μ$, magnetically-dominated turbulence is produced. Turbulence gives rise to a large-scale dynamo instability, which we find to be caused by the magnetic alpha effect. Our results have consequences for the dynamics of certain high-energy plasmas, such as the early Universe or proto-neutron stars. △ Less

Submitted 11 April, 2024; originally announced April 2024.

Comments: 15 pages, 11 figures, submitted to PRD

Report number: NORDITA-2024-008

Resistively controlled primordial magnetic turbulence decay

Authors: A. Brandenburg, A. Neronov, F. Vazza

Abstract: Magnetic fields generated in the early Universe undergo turbulent decay during the radiation-dominated era. The decay is governed by a decay exponent and a decay time. It has been argued that the latter is prolonged by magnetic reconnection, which depends on the microphysical resistivity and viscosity. Turbulence, on the other hand, is not usually expected to be sensitive to microphysical dissipat… ▽ More Magnetic fields generated in the early Universe undergo turbulent decay during the radiation-dominated era. The decay is governed by a decay exponent and a decay time. It has been argued that the latter is prolonged by magnetic reconnection, which depends on the microphysical resistivity and viscosity. Turbulence, on the other hand, is not usually expected to be sensitive to microphysical dissipation, which affects only very small scales. We want to test and quantify the reconnection hypothesis in decaying hydromagnetic turbulence. We performed high-resolution numerical simulations with zero net magnetic helicity using the Pencil Code with up to $2048^3$ mesh points and relate the decay time to the Alfvén time for different resistivities and viscosities. The decay time is found to be longer than the Alfvén time by a factor that increases with increasing Lundquist number to the 1/4 power. The decay exponent is as expected from the conservation of the Hosking integral, but a timescale dependence on resistivity is unusual for developed turbulence and not found for hydrodynamic turbulence. In two dimensions, the Lundquist number dependence is shown to be leveling off above values of $\approx25,000$, independently of the value of the viscosity. Our numerical results suggest that resistivity effects have been overestimated in earlier work. Instead of reconnection, it may be the magnetic helicity density in smaller patches that is responsible for the resistively slow decay. The leveling off at large Lundquist number cannot currently be confirmed in three dimensions. △ Less

Submitted 11 April, 2024; v1 submitted 16 January, 2024; originally announced January 2024.

Comments: 15 pages, 17 figures, 3 tables; A&A (in press)

Report number: NORDITA-2024-001

Journal ref: Astron. Astrophys. 687, A186 (2024)

Backreaction of axion-SU(2) dynamics during inflation

Authors: Oksana Iarygina, Evangelos I. Sfakianakis, Ramkishor Sharma, Axel Brandenburg

Abstract: We consider the effects of backreaction on axion-SU(2) dynamics during inflation. We use the linear evolution equations for the gauge field modes and compute their backreaction on the background quantities numerically using the Hartree approximation. We show that the spectator chromo-natural inflation attractor is unstable when back-reaction becomes important. Working within the constraints of the… ▽ More We consider the effects of backreaction on axion-SU(2) dynamics during inflation. We use the linear evolution equations for the gauge field modes and compute their backreaction on the background quantities numerically using the Hartree approximation. We show that the spectator chromo-natural inflation attractor is unstable when back-reaction becomes important. Working within the constraints of the linear mode equations, we find a new dynamical attractor solution for the axion field and the vacuum expectation value of the gauge field, where the latter has an opposite sign with respect to the chromo-natural inflation solution. Our findings are of particular interest to the phenomenology of axion-SU(2) inflation, as they demonstrate the instability of the usual trajectory due to large backreaction effects. The viable parameter space of the model becomes significantly altered, provided future non-Abelian lattice simulations confirm the existence of the new dynamical attractor. In addition, the backreaction effects lead to characteristic oscillatory features in the primordial gravitational wave background that are potentially detectable with upcoming gravitational wave detectors. △ Less

Submitted 19 June, 2024; v1 submitted 13 November, 2023; originally announced November 2023.

Comments: 25 pages, 13 figures, 2 tables

Report number: NORDITA-2023-070

Journal ref: JCAP04(2024)018

Shallow relic gravitational wave spectrum with acoustic peak

Authors: Ramkishor Sharma, Jani Dahl, Axel Brandenburg, Mark Hindmarsh

Abstract: We study the gravitational wave (GW) spectrum produced by acoustic waves in the early universe, such as would be produced by a first order phase transition, focusing on the low-frequency side of the peak. We confirm with numerical simulations the Sound Shell model prediction of a steep rise with wave number $k$ of $k^9$ to a peak whose magnitude grows at a rate $(H/k_\text{p})H$, where $H$ is the… ▽ More We study the gravitational wave (GW) spectrum produced by acoustic waves in the early universe, such as would be produced by a first order phase transition, focusing on the low-frequency side of the peak. We confirm with numerical simulations the Sound Shell model prediction of a steep rise with wave number $k$ of $k^9$ to a peak whose magnitude grows at a rate $(H/k_\text{p})H$, where $H$ is the Hubble rate and $k_\text{p}$ the peak wave number, set by the peak wave number of the fluid velocity power spectrum. We also show that hitherto neglected terms give a shallower part with amplitude $(H/k_\text{p})^2$ in the range $H \lesssim k \lesssim k_\text{p}$, which in the limit of small $H/k$ rises as $k$. This linear rise has been seen in other modelling and also in direct numerical simulations. The relative amplitude between the linearly rising part and the peak therefore depends on the peak wave number of the velocity spectrum and the lifetime of the source, which in an expanding background is bounded above by the Hubble time $H^{-1}$. For slow phase transitions, which have the lowest peak wave number and the loudest signals, the acoustic GW peak appears as a localized enhancement of the spectrum, with a rise to the peak less steep than $k^9$. The shape of the peak, absent in vortical turbulence, may help to lift degeneracies in phase transition parameter estimation at future GW observatories. △ Less

Submitted 10 November, 2023; v1 submitted 24 August, 2023; originally announced August 2023.

Comments: 21 pages, 10 figures, Accepted for publication in JCAP

Report number: NORDITA-2023-051, HIP-2023-13/TH

Journal ref: J. Cosmol. Astropart. Phys. 12 (2023) 042

Electromagnetic conversion into kinetic and thermal energies

Authors: Axel Brandenburg, Nousaba Nasrin Protiti

Abstract: Conversion of electromagnetic energy into magnetohydrodynamic energy occurs when the electric conductivity changes from negligible to finite values. This process is relevant during the epoch of reheating of the early Universe at the end of inflation and before the emergence of the radiation-dominated era. We find that conversion into kinetic and thermal energies is primarily the result of electric… ▽ More Conversion of electromagnetic energy into magnetohydrodynamic energy occurs when the electric conductivity changes from negligible to finite values. This process is relevant during the epoch of reheating of the early Universe at the end of inflation and before the emergence of the radiation-dominated era. We find that conversion into kinetic and thermal energies is primarily the result of electric energy dissipation and that the magnetic energy plays only a secondary role in this process. This means that, since electric energy dominates over magnetic energy during inflation and reheating, significant amounts of electric energy can be converted into magnetohydrodynamic energy when conductivity emerges early enough, before the relevant length scales become stable. △ Less

Submitted 29 August, 2023; v1 submitted 1 August, 2023; originally announced August 2023.

Comments: 22 pages, 16 figures, 3 tables, published in Entropy, special issue on Energy Transfer and Dissipation in Plasma Turbulence

Report number: NORDITA-2023-043

Journal ref: Entropy 25, 1270 (2023)

Chiral anomaly and dynamos from inhomogeneous chemical potential fluctuations

Authors: Jennifer Schober, Igor Rogachevskii, Axel Brandenburg

Abstract: In the standard model of particle physics, the chiral anomaly can occur in relativistic plasmas and plays a role in the early Universe, protoneutron stars, heavy-ion collisions, and quantum materials. It gives rise to a magnetic instability if the number densities of left- and right-handed electrically charged fermions are unequal. Using direct numerical simulations, we show this can result just f… ▽ More In the standard model of particle physics, the chiral anomaly can occur in relativistic plasmas and plays a role in the early Universe, protoneutron stars, heavy-ion collisions, and quantum materials. It gives rise to a magnetic instability if the number densities of left- and right-handed electrically charged fermions are unequal. Using direct numerical simulations, we show this can result just from spatial fluctuations of the chemical potential, causing a chiral dynamo instability, magnetically driven turbulence, and ultimately a large-scale magnetic field through the magnetic alpha effect. △ Less

Submitted 5 February, 2024; v1 submitted 27 July, 2023; originally announced July 2023.

Comments: 5 pages, 4 figures

Report number: NORDITA-2023-041

Journal ref: Phys. Rev. Lett. 132, 065101 (2024)

Relic Gravitational Waves from the Chiral Plasma Instability in the Standard Cosmological Model

Authors: Axel Brandenburg, Emma Clarke, Tina Kahniashvili, Andrew J. Long, Guotong Sun

Abstract: In the primordial plasma, at temperatures above the scale of electroweak symmetry breaking, the presence of chiral asymmetries is expected to induce the development of helical hypermagnetic fields through the phenomenon of chiral plasma instability. It results in magnetohydrodynamic turbulence due to the high conductivity and low viscosity and sources gravitational waves that survive in the univer… ▽ More In the primordial plasma, at temperatures above the scale of electroweak symmetry breaking, the presence of chiral asymmetries is expected to induce the development of helical hypermagnetic fields through the phenomenon of chiral plasma instability. It results in magnetohydrodynamic turbulence due to the high conductivity and low viscosity and sources gravitational waves that survive in the universe today as a stochastic polarized gravitational wave background. In this article, we show that this scenario only relies on Standard Model physics, and therefore the observable signatures, namely the relic magnetic field and gravitational background, are linked to a single parameter controlling the initial chiral asymmetry. We estimate the magnetic field and gravitational wave spectra, and validate these estimates with 3D numerical simulations. △ Less

Submitted 8 May, 2024; v1 submitted 18 July, 2023; originally announced July 2023.

Comments: 8 pages, 4 figures, 1 table; v2 updated to match published version, results unchanged and PRD accepted; v3: fixed broken reference link

Report number: NORDITA-2023-034

Journal ref: Phys. Rev. D 109, 043534 (2024)

Inverse cascading for initial MHD turbulence spectra between Saffman and Batchelor

Authors: Axel Brandenburg, Ramkishor Sharma, Tanmay Vachaspati

Abstract: In decaying magnetohydrodynamic (MHD) turbulence with a strong magnetic field, the spectral magnetic energy density is known to increase with time at small wavenumbers $k$, provided the spectrum at low $k$ is sufficiently steep. This process is called inverse cascading and occurs for an initial Batchelor spectrum, where the magnetic energy per linear wavenumber interval increases like $k^4$. For a… ▽ More In decaying magnetohydrodynamic (MHD) turbulence with a strong magnetic field, the spectral magnetic energy density is known to increase with time at small wavenumbers $k$, provided the spectrum at low $k$ is sufficiently steep. This process is called inverse cascading and occurs for an initial Batchelor spectrum, where the magnetic energy per linear wavenumber interval increases like $k^4$. For an initial Saffman spectrum that is proportional to $k^2$, however, inverse cascading has not been found in the past. We study here the case of an intermediate $k^3$ spectrum, which may be relevant for magnetogenesis in the early Universe during the electroweak epoch. This case is not well understood in view of the standard Taylor expansion of the magnetic energy spectrum for small $k$. Using high resolution MHD simulations, we show that also in this case there is inverse cascading with a strength just as expected from the conservation of the Hosking integral, which governs the decay of an initial Batchelor spectrum. Even for shallower $k^α$ spectra with spectral index $α>3/2$, our simulations suggest a spectral increase at small $k$ with time $t$ proportional to $t^{4α/9-2/3}$. The critical spectral index of $α=3/2$ is related to the slope of the spectral envelope in the Hosking phenomenology. Our simulations with $2048^3$ mesh points now suggest inverse cascading even for an initial Saffman spectrum. △ Less

Submitted 26 October, 2023; v1 submitted 10 July, 2023; originally announced July 2023.

Comments: 18 pages, 9 figures, 3 tables, submitted to J. Plasma Physics, improved evidence for inverse cascading for nonhelical $k^2$ spectrum

Report number: NORDITA-2023-035

Journal ref: J. Plasma Phys. 89, 905890606 (2023)

Confronting the neutron star population with inverse cascades

Authors: Nikhil Sarin, Axel Brandenburg, Brynmor Haskell

Abstract: The origin and evolution of magnetic fields of neutron stars from birth has long been a source of debate. Here, motivated by recent simulations of the Hall cascade with magnetic helicity, we invoke a model where the large-scale magnetic field of neutron stars grows as a product of small-scale turbulence through an inverse cascade. We apply this model to a simulated population of neutron stars at b… ▽ More The origin and evolution of magnetic fields of neutron stars from birth has long been a source of debate. Here, motivated by recent simulations of the Hall cascade with magnetic helicity, we invoke a model where the large-scale magnetic field of neutron stars grows as a product of small-scale turbulence through an inverse cascade. We apply this model to a simulated population of neutron stars at birth and show how this model can account for the evolution of such objects across the $P\dot{P}$ diagram, explaining both pulsar and magnetar observations. Under the assumption that small-scale turbulence is responsible for large-scale magnetic fields, we place a lower limit on the spherical harmonic degree of the energy-carrying magnetic eddies of $\approx 40$. Our results favor the presence of a highly resistive pasta layer at the base of the neutron star crust. We further discuss the implications of this paradigm on direct observables, such as the nominal age and braking index of pulsars. △ Less

Submitted 26 July, 2023; v1 submitted 23 May, 2023; originally announced May 2023.

Comments: Published in ApJL. 8 pages, 3 figures

Report number: NORDITA-2023-023

Journal ref: Astrophys. J. Lett. 952, L21 (2023)

Chiral magnetohydrodynamics with zero total chirality

Authors: Axel Brandenburg, Kohei Kamada, Kyohei Mukaida, Kai Schmitz, Jennifer Schober

Abstract: We study the evolution of magnetic fields coupled with chiral fermion asymmetry in the framework of chiral magnetohydrodynamics with zero initial total chirality. The initial magnetic field has a turbulent spectrum peaking at a certain characteristic scale and is fully helical with positive helicity. The initial chiral chemical potential is spatially uniform and negative. We consider two opposite… ▽ More We study the evolution of magnetic fields coupled with chiral fermion asymmetry in the framework of chiral magnetohydrodynamics with zero initial total chirality. The initial magnetic field has a turbulent spectrum peaking at a certain characteristic scale and is fully helical with positive helicity. The initial chiral chemical potential is spatially uniform and negative. We consider two opposite cases where the ratio of the length scale of the chiral plasma instability (CPI) to the characteristic scale of the turbulence is smaller and larger than unity. These initial conditions might be realized in cosmological models such as certain types of axion inflation. The magnetic field and chiral chemical potential evolve with inverse cascading in such a way that the magnetic helicity and chirality cancel each other at all times. The CPI time scale is found to determine mainly the time when the magnetic helicity spectrum attains negative values at high wave numbers. The turnover time of the energy-carrying eddies, on the other hand, determines the time when the peak of the spectrum starts to shift to smaller wave numbers via an inverse cascade. The onset of helicity decay is determined by the time when the chiral magnetic effect becomes efficient at the peak of the initial magnetic energy spectrum. When spin flipping is important, the chiral chemical potential vanishes and the magnetic helicity becomes constant, which leads to a faster increase of the correlation length, as expected from magnetic helicity conservation. This also happens when the initial total chirality is imbalanced. Our findings have important implications for baryogenesis after axion inflation. △ Less

Submitted 22 August, 2023; v1 submitted 13 April, 2023; originally announced April 2023.

Comments: 24 pages, 21 figures, 6 tables, corrected Fig.17, Phys. Rev. D, in press

Report number: NORDITA-2023-014, RESCEU-5/23, KEK-TH-2504, MS-TP-23-13

Journal ref: Phys. Rev. D 108, 063529 (2023)

Turbulent processes and mean-field dynamo

Authors: Axel Brandenburg, Detlef Elstner, Youhei Masada, Valery Pipin

Abstract: Mean-field dynamo theory has important applications in solar physics and galactic magnetism. We discuss some of the many turbulence effects relevant to the generation of large-scale magnetic fields in the solar convection zone. The mean-field description is then used to illustrate the physics of the $α$ effect, turbulent pumping, turbulent magnetic diffusivity, and other effects on a modern solar… ▽ More Mean-field dynamo theory has important applications in solar physics and galactic magnetism. We discuss some of the many turbulence effects relevant to the generation of large-scale magnetic fields in the solar convection zone. The mean-field description is then used to illustrate the physics of the $α$ effect, turbulent pumping, turbulent magnetic diffusivity, and other effects on a modern solar dynamo model. We also discuss how turbulence transport coefficients are derived from local simulations of convection and then used in mean-field models. △ Less

Submitted 21 November, 2023; v1 submitted 22 March, 2023; originally announced March 2023.

Comments: 74 pages, 20 figures, 2 tables, Space Science Reviews, accepted, special issue "Solar and stellar dynamos: a new era" (corrected fig17)

Report number: NORDITA-2023-011

Journal ref: Spa. Sci. Rev. 219, 55 (2023)

Decay law of magnetic turbulence with helicity balanced by chiral fermions

Authors: Axel Brandenburg, Kohei Kamada, Jennifer Schober

Abstract: In plasmas composed of massless electrically charged fermions, chirality can be interchanged with magnetic helicity while preserving the total chirality through the quantum chiral anomaly. The decay of turbulent energy in plasmas such as those in the early Universe and compact stars is usually controlled by certain conservation laws. In the case of zero total chirality, when the magnetic helicity… ▽ More In plasmas composed of massless electrically charged fermions, chirality can be interchanged with magnetic helicity while preserving the total chirality through the quantum chiral anomaly. The decay of turbulent energy in plasmas such as those in the early Universe and compact stars is usually controlled by certain conservation laws. In the case of zero total chirality, when the magnetic helicity density balances with the appropriately scaled chiral chemical potential to zero, the total chirality no longer determines the decay. We propose that in such a case, an adaptation to the Hosking integral, which is conserved in nonhelical magnetically dominated turbulence, controls the decay in turbulence with helicity balanced by chiral fermions. We show, using a high resolution numerical simulation, that this is indeed the case. The magnetic energy density decays and the correlation length increases with time just like in nonhelical turbulence with vanishing chiral chemical potential. But here, the magnetic helicity density is nearly maximum and shows a scaling with time $t$ proportional to $t^{-2/3}$. This is unrelated to the $t^{-2/3}$ decay of magnetic {\it energy} in fully helical magnetic turbulence. The modulus of the chiral chemical potential decays in the same fashion. This is much slower than the exponential decay previously expected in theories of asymmetric baryon production from the hypermagnetic helicity decay after axion inflation. △ Less

Submitted 12 May, 2023; v1 submitted 1 February, 2023; originally announced February 2023.

Comments: 6 pages, 4 figures, corrected Fig.2

Report number: NORDITA-2023-001, RESCEU-1/23

Journal ref: Phys. Rev. Res. 5, L022028 (2023)

Modified propagation of gravitational waves from the early radiation era

Authors: Yutong He, Alberto Roper Pol, Axel Brandenburg

Abstract: We study the propagation of cosmological gravitational wave (GW) backgrounds from the early radiation era until the present day in modified theories of gravity. Comparing to general relativity (GR), we study the effects that modified gravity parameters, such as the GW friction $α_{\rm M}$ and the tensor speed excess $α_{\rm T}$, have on the present-day GW spectrum. We use both the WKB estimate, wh… ▽ More We study the propagation of cosmological gravitational wave (GW) backgrounds from the early radiation era until the present day in modified theories of gravity. Comparing to general relativity (GR), we study the effects that modified gravity parameters, such as the GW friction $α_{\rm M}$ and the tensor speed excess $α_{\rm T}$, have on the present-day GW spectrum. We use both the WKB estimate, which provides an analytical description but fails at superhorizon scales, and numerical simulations that allow us to go beyond the WKB approximation. We show that a constant $α_{\rm T}$ makes relatively insignificant changes to the GR solution, especially taking into account the constraints on its value from GW observations by the LIGO--Virgo collaboration, while $α_{\rm M}$ can introduce modifications to the spectral slopes of the GW energy spectrum in the low-frequency regime depending on the considered time evolution of $α_{\rm M}$. The latter effect is additional to the damping or growth occurring equally at all scales that can be predicted by the WKB approximation. In light of the recent observations by pulsar timing array (PTA) collaborations, and the potential observations by future detectors such as SKA, LISA, DECIGO, BBO, or ET, we show that, in most of the cases, constraints cannot be placed on the effects of $α_{\rm M}$ and the initial GW energy density $\mathcal{E}_{\rm GW}^*$ separately, but only on the combined effects of the two, unless the signal is observed at different frequency ranges. In particular, we provide some constraints on the combined effects from the reported PTA observations. △ Less

Submitted 20 June, 2023; v1 submitted 12 December, 2022; originally announced December 2022.

Comments: 36 pages, 11 figures, 2 tables

Report number: NORDITA-2022-080

Journal ref: JCAP06(2023)025

Hosking integral in nonhelical Hall cascade

Authors: Axel Brandenburg

Abstract: The Hosking integral, which characterizes magnetic helicity fluctuations in subvolumes, is known to govern the decay of magnetically dominated turbulence. Here we show that, when the evolution of the magnetic field is controlled by the motion of electrons only, as in neutron star crusts, the decay of the magnetic field is still controlled by the Hosking integral, but now it has effectively differe… ▽ More The Hosking integral, which characterizes magnetic helicity fluctuations in subvolumes, is known to govern the decay of magnetically dominated turbulence. Here we show that, when the evolution of the magnetic field is controlled by the motion of electrons only, as in neutron star crusts, the decay of the magnetic field is still controlled by the Hosking integral, but now it has effectively different dimensions than in ordinary magnetohydrodynamic (MHD) turbulence. This causes the correlation length to increase with time $t$ like $t^{4/13}$ instead of $t^{4/9}$ in MHD. The magnetic energy density decreases like $t^{-10/13}$, which is slower than in MHD, where it decays like $t^{-10/9}$. These new analytic results agree with earlier numerical simulations for the nonhelical Hall cascade. △ Less

Submitted 27 December, 2022; v1 submitted 25 November, 2022; originally announced November 2022.

Comments: 9 pages, 4 figures, 2 tables

Report number: NORDITA-2022-083

Journal ref: J. Plasma Phys. 89, 175890101 (2023)

Galactic Dynamos

Authors: Axel Brandenburg, Evangelia Ntormousi

Abstract: 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 to explain either particular aspects of the problem or to reproduce galactic magnetic fields globally. Our main con… ▽ More 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 to explain either 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 only be explained if helical magnetic fields of small or moderate length scales can rapidly be 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. △ Less

Submitted 30 August, 2023; v1 submitted 7 November, 2022; originally announced November 2022.

Comments: 47 pages, 15 figures, 4 tables

Report number: NORDITA-2022-077

Journal ref: Annu. Rev. Astron. Astrophys. 61, 561-606 (2023)

Inflationary and phase-transitional primordial magnetic fields in galaxy clusters

Authors: Salome Mtchedlidze, Paola Domínguez-Fernández, Xiaolong Du, Wolfram Schmidt, Axel Brandenburg, Jens Niemeyer, Tina Kahniashvili

Abstract: Primordial magnetic fields (PMFs) are possible candidates for explaining the observed magnetic fields in galaxy clusters. Two competing scenarios of primordial magnetogenesis have been discussed in the literature: inflationary and phase-transitional. We study the amplification of both large- and small-scale correlated magnetic fields, corresponding to inflation- and phase transition-generated PMFs… ▽ More Primordial magnetic fields (PMFs) are possible candidates for explaining the observed magnetic fields in galaxy clusters. Two competing scenarios of primordial magnetogenesis have been discussed in the literature: inflationary and phase-transitional. We study the amplification of both large- and small-scale correlated magnetic fields, corresponding to inflation- and phase transition-generated PMFs, in a massive galaxy cluster. We employ high-resolution magnetohydrodynamic cosmological zoom-in simulations to resolve the turbulent motions in the intracluster medium. We find that the turbulent amplification is more efficient for the large-scale inflationary models, while the phase transition-generated seed fields show moderate growth. The differences between the models are imprinted on the spectral characteristics of the field (such as the amplitude and the shape of the magnetic power spectrum) and therefore on the final correlation length. We find a one order of magnitude difference between the final strengths of the inflation- and phase transition-generated magnetic fields, and a factor of 1.5 difference between their final coherence scales. Thus, the final configuration of the magnetic field retains information about the PMF generation scenarios. Our findings have implications for future extragalactic Faraday rotation surveys with the possibility of distinguishing between different magnetogenesis scenarios. △ Less

Submitted 16 February, 2023; v1 submitted 18 October, 2022; originally announced October 2022.

Comments: 21 pages, 13 figures, comments welcome

Report number: NORDITA-2022-072

Dissipative magnetic structures and scales in small-scale dynamos

Authors: A. Brandenburg, I. Rogachevskii, J. Schober

Abstract: Small-scale dynamos play important roles in modern astrophysics, especially on Galactic and extragalactic scales. Owing to dynamo action, purely hydrodynamic Kolmogorov turbulence hardly exists and is often replaced by hydromagnetic turbulence. Understanding the size of dissipative magnetic structures is important in estimating the time scale of Galactic scintillation and other observational and t… ▽ More Small-scale dynamos play important roles in modern astrophysics, especially on Galactic and extragalactic scales. Owing to dynamo action, purely hydrodynamic Kolmogorov turbulence hardly exists and is often replaced by hydromagnetic turbulence. Understanding the size of dissipative magnetic structures is important in estimating the time scale of Galactic scintillation and other observational and theoretical aspects of interstellar and intergalactic small-scale dynamos. Here we show that, during the kinematic phase of the small-scale dynamo, the cutoff wavenumber of the magnetic energy spectra scales as expected for large magnetic Prandtl numbers, but continues in the same way also for moderately small values - contrary to what is expected. For a critical magnetic Prandtl number of about 0.3, the dissipative and resistive cutoffs are found to occur at the same wavenumber. In the nonlinearly saturated regime, the critical magnetic Prandtl number becomes unity. The cutoff scale now has a shallower scaling with magnetic Prandtl number below a value of about three, and a steeper one otherwise compared to the kinematic regime. △ Less

Submitted 14 November, 2022; v1 submitted 18 September, 2022; originally announced September 2022.

Comments: 9 pages, 12 figures, 4 tables, submitted to MNRAS, revised version

Report number: NORDITA 2022-064

Journal ref: Mon. Not. Roy. Astron. Soc. 518, 6367-6375 (2023)

Magnetohydrodynamics predicts heavy-tailed distributions of axion-photon conversion

Authors: Pierluca Carenza, Ramkishor Sharma, M. C. David Marsh, Axel Brandenburg, Eike Ravensburg

Abstract: The interconversion of axionlike particles (ALPs) and photons in magnetised astrophysical environments provides a promising route to search for ALPs. The strongest limits to date on light ALPs use galaxy clusters as ALP-photon converters. However, such studies traditionally rely on simple models of the cluster magnetic fields, with the state-of-the-art being Gaussian random fields (GRFs). We prese… ▽ More The interconversion of axionlike particles (ALPs) and photons in magnetised astrophysical environments provides a promising route to search for ALPs. The strongest limits to date on light ALPs use galaxy clusters as ALP-photon converters. However, such studies traditionally rely on simple models of the cluster magnetic fields, with the state-of-the-art being Gaussian random fields (GRFs). We present the first systematic study of ALP-photon conversion in more realistic, turbulent fields from dedicated magnetohydrodynamic (MHD) simulations, which we compare with GRF models. For GRFs, we analytically derive the distribution of conversion ratios at fixed energy and find that it follows an exponential law. We find that the MHD models agree with the exponential law for typical, small-amplitude mixings but exhibit distinctly heavy tails for rare and large mixings. We explain how non-Gaussian features, e.g.~coherent structures and local spikes in the MHD magnetic field, are responsible for the heavy tail. Our results suggest that limits placed on ALPs using GRFs are robust. △ Less

Submitted 16 November, 2023; v1 submitted 8 August, 2022; originally announced August 2022.

Comments: 18 pages, 9 figures. v2: major changes to match the published version. Extended appendix to include details on the statistical analysis

Report number: NORDITA 2022-060

Journal ref: Phys.Rev.D 108 (2023) 10, 103029

Batchelor, Saffman, and Kazantsev spectra in galactic small-scale dynamos

Authors: Axel Brandenburg, Hongzhe Zhou, Ramkishor Sharma

Abstract: The magnetic fields in galaxy clusters and probably also in the interstellar medium are believed to be generated by a small-scale dynamo. Theoretically, during its kinematic stage, it is characterized by a Kazantsev spectrum, which peaks at the resistive scale. It is only slightly shallower than the Saffman spectrum that is expected for random and causally connected magnetic fields. Causally disco… ▽ More The magnetic fields in galaxy clusters and probably also in the interstellar medium are believed to be generated by a small-scale dynamo. Theoretically, during its kinematic stage, it is characterized by a Kazantsev spectrum, which peaks at the resistive scale. It is only slightly shallower than the Saffman spectrum that is expected for random and causally connected magnetic fields. Causally disconnected fields have the even steeper Batchelor spectrum. Here we show that all three spectra are present in the small-scale dynamo. During the kinematic stage, the Batchelor spectrum occurs on scales larger than the energy-carrying scale of the turbulence, and the Kazantsev spectrum on smaller scales within the inertial range of the turbulence -- even for a magnetic Prandtl number of unity. In the saturated state, the dynamo develops a Saffman spectrum on large scales, suggestive of the build-up of long-range correlations. At large magnetic Prandtl numbers, elongated structures are seen in synthetic synchrotron emission maps showing the parity-even E polarization. We also observe a significant excess in the E polarization over the parity-odd B polarization at subresistive scales, and a deficiency at larger scales. This finding is at odds with the observed excess in the Galactic microwave foreground emission, which is believed to be associated with larger scales. The E and B polarizations may be highly non-Gaussian and skewed in the kinematic regime of the dynamo. For dust emission, however, the polarized emission is always nearly Gaussian, and the excess in the E polarization is much weaker. △ Less

Submitted 3 November, 2022; v1 submitted 19 July, 2022; originally announced July 2022.

Comments: 14 pages, 25 figures, 5 tables, MNRAS, in press

Report number: NORDITA 2022-051

Journal ref: Mon. Not. Roy. Astron. Soc. 518, 3312-3325 (2023)

Scaling of the Hosking integral in decaying magnetically-dominated turbulence

Authors: Hongzhe Zhou, Ramkishor Sharma, Axel Brandenburg

Abstract: The Saffman helicity invariant of Hosking and Schekochihin (2021, PRX 11, 041005), which we here call the Hosking integral, has emerged as an important quantity that may govern the decay properties of magnetically dominated nonhelical turbulence. Using a range of different computational methods, we confirm that this quantity is indeed gauge-invariant and nearly perfectly conserved in the limit of… ▽ More The Saffman helicity invariant of Hosking and Schekochihin (2021, PRX 11, 041005), which we here call the Hosking integral, has emerged as an important quantity that may govern the decay properties of magnetically dominated nonhelical turbulence. Using a range of different computational methods, we confirm that this quantity is indeed gauge-invariant and nearly perfectly conserved in the limit of large Lundquist numbers. For direct numerical simulations with ordinary viscosity and magnetic diffusivity operators, we find that the solution develops in a nearly self-similar fashion. In a diagram quantifying the instantaneous decay coefficients of magnetic energy and integral scale, we find that the solution evolves along a line that is indeed suggestive of the governing role of the Hosking integral. The solution settles near a line in this diagram that is expected for a self-similar evolution of the magnetic energy spectrum. The solution will settle in a slightly different position when the magnetic diffusivity decreases with time, which would be compatible with the decay being governed by the reconnection time scale rather than the Alfvén time. △ Less

Submitted 26 September, 2022; v1 submitted 14 June, 2022; originally announced June 2022.

Comments: 19 pages, 7 figures, 2 tables; accepted by JPP

Report number: NORDITA 2022-040

Low frequency tail of gravitational wave spectra from hydromagnetic turbulence

Authors: Ramkishor Sharma, Axel Brandenburg

Abstract: Hydrodynamic and magnetohydrodynamic turbulence in the early Universe can drive gravitational waves (GWs) and imprint their spectrum onto that of GWs, which might still be observable today. We study the production of the GW background from freely decaying magnetohydrodynamic turbulence from helical and nonhelical initial magnetic fields. To understand the produced GW spectra, we develop a simple m… ▽ More Hydrodynamic and magnetohydrodynamic turbulence in the early Universe can drive gravitational waves (GWs) and imprint their spectrum onto that of GWs, which might still be observable today. We study the production of the GW background from freely decaying magnetohydrodynamic turbulence from helical and nonhelical initial magnetic fields. To understand the produced GW spectra, we develop a simple model on the basis of the evolution of the magnetic stress tensor. We find that the GW spectra obtained in this model reproduce those obtained in numerical simulations if we consider the detailed time evolution of the low-frequency tail of the stress spectrum from numerical simulations. We also show that the shapes of the produced GW frequency spectra are different for helical and nonhelical cases for the same initial magnetic energy spectra. Such differences can help distinguish helical and nonhelical initial magnetic fields from a polarized background of GWs -- especially when the expected circular polarization cannot be detected directly. △ Less

Submitted 28 November, 2022; v1 submitted 31 May, 2022; originally announced June 2022.

Comments: 12 pages, 12 figures, 3 tables, Accepted for publication in Phys. Rev. D

Report number: NORDITA 2022-036

Journal ref: Phys. Rev. D 106, 103536 (2022)

A Comparative Analysis of Machine-learning Models for Solar Flare Forecasting: Identifying High-performing Active Region Flare Indicators

Authors: Suvadip Sinha, Om Gupta, Vishal Singh, B. Lekshmi, Dibyendu Nandy, Dhrubaditya Mitra, Saikat Chatterjee, Sourangshu Bhattacharya, Saptarshi Chatterjee, Nandita Srivastava, Axel Brandenburg, Sanchita Pal

Abstract: Solar flares create adverse space weather impacting space and Earth-based technologies. However, the difficulty of forecasting flares, and by extension severe space weather, is accentuated by the lack of any unique flare trigger or a single physical pathway. Studies indicate that multiple physical properties contribute to active region flare potential, compounding the challenge. Recent development… ▽ More Solar flares create adverse space weather impacting space and Earth-based technologies. However, the difficulty of forecasting flares, and by extension severe space weather, is accentuated by the lack of any unique flare trigger or a single physical pathway. Studies indicate that multiple physical properties contribute to active region flare potential, compounding the challenge. Recent developments in machine learning (ML) have enabled analysis of higher-dimensional data leading to increasingly better flare forecasting techniques. However, consensus on high-performing flare predictors remains elusive. In the most comprehensive study to date, we conduct a comparative analysis of four popular ML techniques (k-nearest neighbor, logistic regression, random forest classifier, and support vector machine) by training these on magnetic parameters obtained from the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO) for the entirety of solar cycle 24. We demonstrate that the logistic regression and support vector machine algorithms perform extremely well in forecasting active region flaring potential. The logistic regression algorithm returns the highest true skill score of $0.967 \pm 0.018$, possibly the highest classification performance achieved with any strictly parametric study. From a comparative assessment, we establish that the magnetic properties like total current helicity, total vertical current density, total unsigned flux, R_VALUE, and total absolute twist are the top-performing flare indicators. We also introduce and analyze two new performance metrics, namely, severe and clear space weather indicators. Our analysis constrains the most successful ML algorithms and identifies physical parameters that contribute most to active region flare productivity. △ Less

Submitted 2 November, 2022; v1 submitted 12 April, 2022; originally announced April 2022.

Report number: NORDITA-2022-025

Journal ref: Astrophys. J. 935, 45 (2022)

Dynamo effect in unstirred self-gravitating turbulence

Authors: Axel Brandenburg, Evangelia Ntormousi

Abstract: In many astrophysical environments, self-gravity can generate kinetic energy, which, in principle, is available for driving dynamo action. Using direct numerical simulations, we show that in unstirred self-gravitating subsonic turbulence with helicity and a magnetic Prandtl number of unity, there is a critical magnetic Reynolds number of about 25 above which the work done against the Lorentz force… ▽ More In many astrophysical environments, self-gravity can generate kinetic energy, which, in principle, is available for driving dynamo action. Using direct numerical simulations, we show that in unstirred self-gravitating subsonic turbulence with helicity and a magnetic Prandtl number of unity, there is a critical magnetic Reynolds number of about 25 above which the work done against the Lorentz force exceeds the Ohmic dissipation. The collapse itself drives predominantly irrotational motions that cannot be responsible for dynamo action. We find that, with a weak magnetic field, one-third of the work done by the gravitational force goes into compressional heating and the remaining two-thirds go first into kinetic energy of the turbulence before a fraction of it is converted further into magnetic and finally thermal energies. Close to the collapse, however, these fractions change toward 1/4 and 3/4 for compressional heating and kinetic energy, respectively. When the magnetic field is strong, the compressional heating fraction is unchanged. Out of the remaining kinetic energy, one quarter goes directly into magnetic energy via work against the Lorentz force. The fraction of vortical motions diminishes in favor of compressive motions that are almost exclusively driven by the Jeans instability. For an initially uniform magnetic field, field amplification at scales larger than those of the initial turbulence are driven by tangling. △ Less

Submitted 14 March, 2022; v1 submitted 7 December, 2021; originally announced December 2021.

Comments: 17 pages, 19 figures, 1 table, resubmitted to MNRAS

Report number: NORDITA 2021-107

Big Bang Nucleosynthesis Limits and Relic Gravitational Waves Detection Prospects

Authors: Tina Kahniashvili, Emma Clarke, Jonathan Stepp, Axel Brandenburg

Abstract: We revisit the big bang nucleosynthesis (BBN) limits on primordial magnetic fields and/or turbulent motions accounting for the decaying nature of turbulent sources between the time of generation and BBN. This leads to larger estimates for the gravitational wave (GW) signal than previously expected. We address the detection prospects through space-based interferometers (for GWs generated around the… ▽ More We revisit the big bang nucleosynthesis (BBN) limits on primordial magnetic fields and/or turbulent motions accounting for the decaying nature of turbulent sources between the time of generation and BBN. This leads to larger estimates for the gravitational wave (GW) signal than previously expected. We address the detection prospects through space-based interferometers (for GWs generated around the electroweak energy scale) as well as pulsar timing arrays and astrometric missions (for GWs generated around the quantum chromodynamics energy scale). △ Less

Submitted 18 June, 2022; v1 submitted 18 November, 2021; originally announced November 2021.

Comments: 9 pages; 6 figures; discussions and references added; typos fixed; conclusions unchanged

Report number: NORDITA-2021-089

Journal ref: Phys.Rev.Lett. 128 (2022) 22, 221301

Leading-order nonlinear gravitational waves from reheating magnetogeneses

Authors: Yutong He, Alberto Roper Pol, Axel Brandenburg

Abstract: We study the leading-order nonlinear gravitational waves (GWs) produced by an electromagnetic (EM) stress in reheating magnetogenesis scenarios. Both nonhelical and helical magnetic fields are considered. By numerically solving the linear and leading-order nonlinear GW equations, we find that the GW energy from the latter is usually larger. We compare their differences in terms of the GW spectrum… ▽ More We study the leading-order nonlinear gravitational waves (GWs) produced by an electromagnetic (EM) stress in reheating magnetogenesis scenarios. Both nonhelical and helical magnetic fields are considered. By numerically solving the linear and leading-order nonlinear GW equations, we find that the GW energy from the latter is usually larger. We compare their differences in terms of the GW spectrum and parameterize the GW energy difference due to the nonlinear term, $Δ\mathcal{E}_{\rm GW}$, in terms of EM energy $\mathcal{E}_{\rm EM}$ as $Δ\mathcal{E}_{\rm GW}=(\tilde p\mathcal{E}_{\rm EM}/k_*)^3$, where $k_*$ is the characteristic wave number, $\tilde p=0.84$ and $0.88$ are found in the nonhelical and helical cases, respectively, with reheating around the QCD energy scale, while $\tilde p=0.45$ is found at the electroweak energy scale. We also compare the polarization spectrum of the linear and nonlinear cases and find that adding the nonlinear term usually yields a decrease in the polarization that is proportional to the EM energy density. We parameterize the fractional polarization suppression as $|Δ\mathcal{P}_{\rm GW}/\mathcal{P}_{\rm GW}|=\tilde r \mathcal{E}_{\rm EM}/k_*$ and find $\tilde r = 1.2 \times 10^{-1}$, $7.2 \times 10^{-4}$, and $3.2 \times 10^{-2}$ for the helical cases with reheating temperatures $T_{\rm r} = 300 {\rm TeV}$, $8 {\rm GeV}$, and $120 {\rm MeV}$, respectively. Prospects of observation by pulsar timing arrays, space-based interferometers, and other novel detection proposals are also discussed. △ Less

Submitted 27 October, 2021; originally announced October 2021.

Comments: 17 pages, 7 figures, 3 tables, to be submitted to PRD

Report number: NORDITA-2021-080

Chirality in Astrophysics

Authors: Axel Brandenburg

Abstract: Chirality, or handedness, enters astrophysics in three distinct ways. Magnetic field and vortex lines tend to be helical and have a systematic twist in the northern and southern hemispheres of a star or a galaxy. Helicity is here driven by external factors. Chirality can also enter at the microphysical level and can then be traced back to the parity-breaking weak force. Finally, chirality can aris… ▽ More Chirality, or handedness, enters astrophysics in three distinct ways. Magnetic field and vortex lines tend to be helical and have a systematic twist in the northern and southern hemispheres of a star or a galaxy. Helicity is here driven by external factors. Chirality can also enter at the microphysical level and can then be traced back to the parity-breaking weak force. Finally, chirality can arise spontaneously, but this requires not only the presence of an instability, but also the action of nonlinearity. Examples can be found both in magnetohydrodynamics and in astrobiology, where homochirality among biomolecules probably got established at the origin of life. In this review, all three types of chirality production will be explored and compared. △ Less

Submitted 5 May, 2022; v1 submitted 14 October, 2021; originally announced October 2021.

Comments: 20 pages, 6 figures, 1 table, to appear in the proceedings to Nobel Symposium 167: Chiral Matter, eds. E. Babaev, D. Kharzeev, M. Larsson, A. Molochkov, & V. Zhaunerchyk, World Scientific, in press

Report number: NORDITA-2021-088

Evolution of primordial magnetic fields during large-scale structure formation

Authors: Salome Mtchedlidze, Paola Domínguez-Fernández, Xiaolong Du, Axel Brandenburg, Tina Kahniashvili, Shane O'Sullivan, Wolfram Schmidt, Marcus Brüggen

Abstract: Primordial magnetic fields could explain the large-scale magnetic fields present in the Universe. Inflation and phase transitions in the early Universe could give rise to such fields with unique characteristics. We investigate the magneto-hydrodynamic evolution of these magnetogenesis scenarios with cosmological simulations. We evolve inflation-generated magnetic fields either as (i) uniform (homo… ▽ More Primordial magnetic fields could explain the large-scale magnetic fields present in the Universe. Inflation and phase transitions in the early Universe could give rise to such fields with unique characteristics. We investigate the magneto-hydrodynamic evolution of these magnetogenesis scenarios with cosmological simulations. We evolve inflation-generated magnetic fields either as (i) uniform (homogeneous) or as (ii) scale-invariant stochastic fields, and phase transition-generated ones either as (iii) helical or as (iv) non-helical fields from the radiation-dominated epoch. We find that the final distribution of magnetic fields in the simulated cosmic web shows a dependence on the initial strength and the topology of the seed field. Thus, the observed field configuration retains information on the initial conditions at the moment of the field generation. If detected, primordial magnetic field observations would open a new window for indirect probes of the early universe. The differences between the competing models are revealed on the scale of galaxy clusters, bridges, as well as filaments and voids. The distinctive spectral evolution of different seed fields produces imprints on the correlation length today. We discuss how the differences between rotation measures from highly ionized regions can potentially be probed with forthcoming surveys. △ Less

Submitted 20 April, 2022; v1 submitted 28 September, 2021; originally announced September 2021.

Comments: 28 pages, 17 figures, accepted for publication in ApJ

Report number: NORDITA 2021-082

Journal ref: Astrophys. J. 929, 127 (2022)

Dynamo instabilities in plasmas with inhomogeneous chiral chemical potential

Authors: Jennifer Schober, Igor Rogachevskii, Axel Brandenburg

Abstract: We study the dynamics of magnetic fields in chiral magnetohydrodynamics, which takes into account the effects of an additional electric current related to the chiral magnetic effect in high energy plasmas. We perform direct numerical simulations, considering weak seed magnetic fields and inhomogeneities of the chiral chemical potential mu_5 with a zero mean. We demonstrate that a small-scale chira… ▽ More We study the dynamics of magnetic fields in chiral magnetohydrodynamics, which takes into account the effects of an additional electric current related to the chiral magnetic effect in high energy plasmas. We perform direct numerical simulations, considering weak seed magnetic fields and inhomogeneities of the chiral chemical potential mu_5 with a zero mean. We demonstrate that a small-scale chiral dynamo can occur in such plasmas if fluctuations of mu_5 are correlated on length scales that are much larger than the scale on which the dynamo growth rate reaches its maximum. Magnetic fluctuations grow by many orders of magnitude due to the small-scale chiral dynamo instability. Once the nonlinear backreaction of the generated magnetic field on fluctuations of mu_5 sets in, the ratio of these scales decreases and the dynamo saturates. When magnetic fluctuations grow sufficiently to drive turbulence via the Lorentz force before reaching maximum field strength, an additional mean-field dynamo phase is identified. The mean magnetic field grows on a scale that is larger than the integral scale of turbulence after the amplification of the fluctuating component saturates. The growth rate of the mean magnetic field is caused by a magnetic alpha effect that is proportional to the current helicity. With the onset of turbulence, the power spectrum of mu_5 develops a universal k^(-1) scaling independently of its initial shape, while the magnetic energy spectrum approaches a k^(-3) scaling. △ Less

Submitted 12 February, 2022; v1 submitted 27 July, 2021; originally announced July 2021.

Comments: 22 pages, 22 figures

Report number: NORDITA-2021-067

Journal ref: Phys. Rev. D 105, 043507 (2022)

Production of a chiral magnetic anomaly with emerging turbulence and mean-field dynamo action

Authors: Jennifer Schober, Igor Rogachevskii, Axel Brandenburg

Abstract: In relativistic magnetized plasmas, asymmetry in the number densities of left- and right-handed fermions, i.e., a non-zero chiral chemical potential mu_5, leads to an electric current along the magnetic field. This causes a chiral dynamo instability for a uniform mu_5, but our simulations reveal a dynamo even for fluctuating mu_5 with zero mean. It produces magnetically-dominated turbulence and ge… ▽ More In relativistic magnetized plasmas, asymmetry in the number densities of left- and right-handed fermions, i.e., a non-zero chiral chemical potential mu_5, leads to an electric current along the magnetic field. This causes a chiral dynamo instability for a uniform mu_5, but our simulations reveal a dynamo even for fluctuating mu_5 with zero mean. It produces magnetically-dominated turbulence and generates mean magnetic fields via the magnetic alpha effect. Eventually, a universal scale-invariant k^(-1) spectrum of mu_5 and a k^(-3) magnetic spectrum are formed independently of the initial condition. △ Less

Submitted 12 February, 2022; v1 submitted 27 July, 2021; originally announced July 2021.

Comments: 6 pages, 4 figures, published in Phys. Rev. Lett

Report number: NORDITA-2021-066

Journal ref: Phys. Rev. Lett. 128, 065002 (2022)

Simulations of helical inflationary magnetogenesis and gravitational waves

Authors: Axel Brandenburg, Yutong He, Ramkishor Sharma

Abstract: Using numerical simulations of helical inflationary magnetogenesis in a low reheating temperature scenario, we show that the magnetic energy spectrum is strongly peaked at a particular wavenumber that depends on the reheating temperature. Gravitational waves (GWs) are produced at frequencies between 3 nHz and 50 mHz for reheating temperatures between 150 MeV and 3x10^5 GeV, respectively. At and be… ▽ More Using numerical simulations of helical inflationary magnetogenesis in a low reheating temperature scenario, we show that the magnetic energy spectrum is strongly peaked at a particular wavenumber that depends on the reheating temperature. Gravitational waves (GWs) are produced at frequencies between 3 nHz and 50 mHz for reheating temperatures between 150 MeV and 3x10^5 GeV, respectively. At and below the peak frequency, the stress spectrum is always found to be that of white noise. This implies a linear increase of GW energy per logarithmic wavenumber interval, instead of a cubic one, as previously thought. Both in the helical and nonhelical cases, the GW spectrum is followed by a sharp drop for frequencies above the respective peak frequency. In this magnetogenesis scenario, the presence of a helical term extends the peak of the GW spectrum and therefore also the position of the aforementioned drop toward larger frequencies compared to the case without helicity. This might make a difference in it being detectable with space interferometers. The efficiency of GW production is found to be almost the same as in the nonhelical case, and independent of the reheating temperature, provided the electromagnetic energy at the end of reheating is fixed to be a certain fraction of the radiation energy density. Also, contrary to the case without helicity, the electric energy is now less than the magnetic energy during reheating. The fractional circular polarization is found to be nearly hundred per cent in a certain range below the peak frequency range. △ Less

Submitted 26 July, 2021; originally announced July 2021.

Comments: 16 pages, 8 figures, 4 tables, submitted to ApJ

Report number: NORDITA-2021-063

Journal ref: Astrophys. J. 922, 192 (2021)

Polarization of gravitational waves from helical MHD turbulent sources

Authors: Alberto Roper Pol, Sayan Mandal, Axel Brandenburg, Tina Kahniashvili

Abstract: We use direct numerical simulations of decaying primordial hydromagnetic turbulence with helicity to compute the resulting gravitational wave (GW) production and its degree of circular polarization. We find a clear dependence of the polarization of the resulting GWs on the fractional helicity of the turbulent source and we show that driven magnetic fields produce GWs more efficiently than magnetic… ▽ More We use direct numerical simulations of decaying primordial hydromagnetic turbulence with helicity to compute the resulting gravitational wave (GW) production and its degree of circular polarization. We find a clear dependence of the polarization of the resulting GWs on the fractional helicity of the turbulent source and we show that driven magnetic fields produce GWs more efficiently than magnetic fields that are initially present, leading to larger spectral amplitudes. The helicity does not have a huge impact on the maximum spectral amplitude in any of the two types of turbulence considered. However, the GW spectrum at wave numbers away from the peak becomes smaller for larger values of the magnetic fractional helicity. The degree of circular polarization approaches zero at frequencies below the peak, and reaches its maximum at the peak. At higher frequencies, it stays finite if the magnetic field is initially present, and it approaches zero if it is driven. We predict that the spectral peak of the GW signal can be detected by LISA if the turbulent energy density is at least $\sim 3\%$ of the radiation energy density, and the characteristic scale is a hundredth of the horizon at the electroweak scale. We show that the resulting GW polarization is unlikely to be detectable by the anisotropies induced by our proper motion in the dipole response function of LISA. Such signals can, however, be detectable by cross-correlating data from the LISA-Taiji network for turbulent energy densities of $\sim 5\%$, and fractional helicity of 0.5 to 1. Second-generation space-based GW detectors, such as BBO and DECIGO, would allow for the detection of a larger range of the GW spectrum and smaller amplitudes of the magnetic field. △ Less

Submitted 4 April, 2022; v1 submitted 12 July, 2021; originally announced July 2021.

Comments: 49 pages, 18 figures, 1 table, accepted in JCAP

Report number: NORDITA-2021-062

Simulating relic gravitational waves from inflationary magnetogenesis

Authors: Axel Brandenburg, Ramkishor Sharma

Abstract: We present three-dimensional direct numerical simulations of the production of magnetic fields and gravitational waves (GWs) in the early Universe during a low energy scale matter-dominated post-inflationary reheating era, and during the early subsequent radiative era, which is strongly turbulent. The parameters of the model are determined such that it avoids a number of known physical problems an… ▽ More We present three-dimensional direct numerical simulations of the production of magnetic fields and gravitational waves (GWs) in the early Universe during a low energy scale matter-dominated post-inflationary reheating era, and during the early subsequent radiative era, which is strongly turbulent. The parameters of the model are determined such that it avoids a number of known physical problems and produces magnetic energy densities between 0.2% and 2% of the critical energy density at the end of reheating. During the subsequent development of a turbulent magnetohydrodynamic cascade, magnetic fields and GWs develop a spectrum that extends to higher frequencies in the millihertz (nanohertz) range for models with reheating temperatures of around 100 GeV (150 MeV) at the beginning of the radiation-dominated era. However, even though the turbulent cascade is fully developed, the GW spectrum shows a sharp drop for frequencies above the peak value. This suggests that the turbulence is less efficient in driving GWs than previously thought. The peaks of the resulting GW spectra may well be in the range accessible to space interferometers, pulsar timing arrays, and other facilities. △ Less

Submitted 12 October, 2021; v1 submitted 7 June, 2021; originally announced June 2021.

Comments: 20 pages, 15 figures, 4 tables, Astrophys. J. 920, 26 (2021)

Report number: NORDITA-2021-041

Journal ref: Astrophys. J. 920, 26 (2021)

Compressible test-field method and its application to shear dynamos

Authors: Maarit J. Käpylä, Matthias Rheinhardt, Axel Brandenburg

Abstract: In this study we present a compressible test-field method (CTFM) for computing $α$ effect and turbulent magnetic diffusivity tensors, as well as those relevant for mean ponderomotive force and mass source, applied to the full MHD equations. We describe the theoretical background of the method, and compare it to the quasi-kinematic test-field method, and to the previously studied variant working in… ▽ More In this study we present a compressible test-field method (CTFM) for computing $α$ effect and turbulent magnetic diffusivity tensors, as well as those relevant for mean ponderomotive force and mass source, applied to the full MHD equations. We describe the theoretical background of the method, and compare it to the quasi-kinematic test-field method, and to the previously studied variant working in simplified MHD (SMHD). We present several test cases using velocity and magnetic fields of the Roberts geometry, and also compare with the imposed-field method. We show that, for moderate imposed field strengths, the nonlinear CTFM (nCTFM) gives results in agreement with the imposed-field method. Comparison of different flavors of the nCTFM in the shear dynamo case also agree up to equipartition field strengths. Some deviations between the CTFM and SMHD variants exist. As a relevant physical application, we study non-helically forced shear flows, which exhibit large-scale dynamo action, and present a re-analysis of low Reynolds number, moderate shear systems, where we previously neglected the pressure gradient in the momentum equation, and found no coherent shear-current effect. Another key difference is that in the earlier study we used magnetic forcing to mimic small-scale dynamo action, while here it is self-consistently driven by purely kinetic forcing. The kinematic CTFM with general validity forms the core of our analysis. We still find no coherent shear-current effect, but do recover strong large-scale dynamo action that, according to our analysis, is driven through the incoherent effects. △ Less

Submitted 3 February, 2022; v1 submitted 2 June, 2021; originally announced June 2021.

Comments: 23 pages, 11 figures, submitted to the Astrophysical Journal

Report number: NORDITA-2021-040

Spectral characterisation of inertial particle clustering in turbulence

Authors: N. E. L. Haugen, A. Brandenburg, C. Sandin, L. Mattsson

Abstract: Clustering of inertial particles is important for many types of astrophysical and geophysical turbulence, but it has been studied predominately for incompressible flows. Here we study compressible flows and compare clustering in both compressively (irrotationally) and vortically (solenoidally) forced turbulence. Vortically and compressively forced flows are driven stochastically either by solenoid… ▽ More Clustering of inertial particles is important for many types of astrophysical and geophysical turbulence, but it has been studied predominately for incompressible flows. Here we study compressible flows and compare clustering in both compressively (irrotationally) and vortically (solenoidally) forced turbulence. Vortically and compressively forced flows are driven stochastically either by solenoidal waves or by circular expansion waves, respectively. For compressively forced flows, the power spectrum of the density of inertial particles is a useful tool for displaying particle clustering relative to the fluid density enhancement. Power spectra are shown to be particularly sensitive for studying large-scale particle clustering, while conventional tools such as radial distribution functions are more suitable for studying small-scale clustering. Our primary finding is that particle clustering through shock interaction is particularly prominent in turbulence driven by spherical expansion waves. It manifests itself through a double-peaked distribution of spectral power as a function of Stokes number. The two peaks are associated with two distinct clustering mechanisms; shock interaction for smaller Stokes numbers and the centrifugal sling effect for larger values. The clustering of inertial particles is associated with the formation of caustics. Such caustics can only be captured in the Lagrangian description, which allows us to assess the relative importance of caustics in vortically and compressively forced turbulence. We show that the statistical noise resulting from the limited number of particles in the Lagrangian description can be removed from the particle power spectra, allowing us a more detailed comparison of the residual spectra. We focus on the Epstein drag law relevant for rarefied gases, but show that our findings apply also to the usual Stokes drag. △ Less

Submitted 16 December, 2021; v1 submitted 3 May, 2021; originally announced May 2021.

Comments: 33 pages, 25 figures, 5 tables, J. Fluid Mech., in press

Report number: NORDITA-2021-033

Journal ref: J. Fluid Mech. 934, A37 (2022)

Tensor spectrum of turbulence-sourced gravitational waves as a constraint on graviton mass

Authors: Yutong He, Axel Brandenburg, Aditya Sinha

Abstract: We consider a generic dispersive massive gravity theory and numerically study its resulting modified energy and strain spectra of tensor gravitational waves (GWs) sourced by (i) fully developed turbulence during the electroweak phase transition (EWPT) and (ii) forced hydromagnetic turbulence during the QCD phase transition (QCDPT). The GW spectra are then computed in both spatial and temporal Four… ▽ More We consider a generic dispersive massive gravity theory and numerically study its resulting modified energy and strain spectra of tensor gravitational waves (GWs) sourced by (i) fully developed turbulence during the electroweak phase transition (EWPT) and (ii) forced hydromagnetic turbulence during the QCD phase transition (QCDPT). The GW spectra are then computed in both spatial and temporal Fourier domains. We find, from the spatial spectra, that the slope modifications are weakly dependent on the eddy size at QCDPT, and, from the temporal spectra, that the modifications are pronounced in the $1$--$10{\rm nHz}$ range -- the sensitivity range of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) -- for a graviton mass $m_{\rm g}$ in the range $2\times10^{-23}{\rm eV}\lesssim m_{\rm g}c^2\lesssim7\times10^{-22}{\rm eV}$. △ Less

Submitted 6 July, 2021; v1 submitted 7 April, 2021; originally announced April 2021.

Comments: 15 pages, 3 figures, 2 tables, comments welcome

Report number: NORDITA 2021-029

Journal ref: JCAP07(2021)015

The scalar, vector, and tensor modes in gravitational wave turbulence simulations

Authors: Axel Brandenburg, Grigol Gogoberidze, Tina Kahniashvili, Sayan Mandal, Alberto Roper Pol, Nakul Shenoy

Abstract: We study the gravitational wave (GW) signal sourced by primordial turbulence that is assumed to be present at cosmological phase transitions like the electroweak and quantum chromodynamics phase transitions. We consider various models of primordial turbulence, such as those with and without helicity, purely hydrodynamical turbulence induced by fluid motions, and magnetohydrodynamic turbulence whos… ▽ More We study the gravitational wave (GW) signal sourced by primordial turbulence that is assumed to be present at cosmological phase transitions like the electroweak and quantum chromodynamics phase transitions. We consider various models of primordial turbulence, such as those with and without helicity, purely hydrodynamical turbulence induced by fluid motions, and magnetohydrodynamic turbulence whose energy can be dominated either by kinetic or magnetic energy, depending on the nature of the turbulence. We also study circularly polarized GWs generated by parity violating sources such as helical turbulence. Our ultimate goal is to determine the efficiency of GW production through different classes of turbulence. We find that the GW energy and strain tend to be large for acoustic or irrotational turbulence, even though its tensor mode amplitude is relatively small at most wave numbers. Only at very small wave numbers is the spectral tensor mode significant, which might explain the efficient GW production in that case. △ Less

Submitted 24 June, 2021; v1 submitted 1 March, 2021; originally announced March 2021.

Comments: 25 pages, 8 figures

Report number: NORDITA-2021-019

Journal ref: Class. Quantum Grav. 38, 145002 (2021)

Can we observe the QCD phase transition-generated gravitational waves through pulsar timing arrays?

Authors: Axel Brandenburg, Emma Clarke, Yutong He, Tina Kahniashvili

Abstract: We perform numerical simulations of gravitational waves (GWs) induced by hydrodynamic and hydromagnetic turbulent sources that might have been present at cosmological quantum chromodynamic (QCD) phase transitions. For turbulent energies of about 4% of the radiation energy density, the typical scale of such motions may have been a sizable fraction of the Hubble scale at that time. The resulting GWs… ▽ More We perform numerical simulations of gravitational waves (GWs) induced by hydrodynamic and hydromagnetic turbulent sources that might have been present at cosmological quantum chromodynamic (QCD) phase transitions. For turbulent energies of about 4% of the radiation energy density, the typical scale of such motions may have been a sizable fraction of the Hubble scale at that time. The resulting GWs are found to have an energy fraction of about $10^{-9}$ of the critical energy density in the nHz range today and may already have been observed by the NANOGrav collaboration. This is further made possible by our findings of shallower spectra proportional to the square root of the frequency for nonhelical hydromagnetic turbulence. This implies more power at low frequencies than for the steeper spectra previously anticipated. The behavior toward higher frequencies depends strongly on the nature of the turbulence. For vortical hydrodynamic and hydromagnetic turbulence, there is a sharp drop of spectral GW energy by up to five orders of magnitude in the presence of helicity, and somewhat less in the absence of helicity. For acoustic hydrodynamic turbulence, the sharp drop is replaced by a power law decay, albeit with a rather steep slope. Our study supports earlier findings of a quadratic scaling of the GW energy with the magnetic energy of the turbulence and inverse quadratic scaling with the peak frequency, which leads to larger GW energies under QCD conditions. △ Less

Submitted 26 July, 2021; v1 submitted 24 February, 2021; originally announced February 2021.

Comments: 14 pages, 8 figures, 3 tables; page numbers, figure captions, references updated and added; discussions expanded; results unchanged and PRD accepted

Report number: NORDITA-2021-016

Journal ref: Phys. Rev. D 104, 043513 (2021)

Relic gravitational waves from the chiral magnetic effect

Authors: Axel Brandenburg, Yutong He, Tina Kahniashvili, Matthias Rheinhardt, Jennifer Schober

Abstract: Relic gravitational waves (GWs) can be produced by primordial magnetic fields. However, not much is known about the resulting GW amplitudes and their dependence on the details of the generation mechanism. Here we treat magnetic field generation through the chiral magnetic effect (CME) as a generic mechanism and explore its dependence on the speed of generation (the product of magnetic diffusivity… ▽ More Relic gravitational waves (GWs) can be produced by primordial magnetic fields. However, not much is known about the resulting GW amplitudes and their dependence on the details of the generation mechanism. Here we treat magnetic field generation through the chiral magnetic effect (CME) as a generic mechanism and explore its dependence on the speed of generation (the product of magnetic diffusivity and characteristic wavenumber) and the speed characterizing the maximum magnetic field strength expected from the CME. When the latter exceeds the former (regime I), the regime applicable to the early universe, we obtain an inverse cascade with moderate GW energy that scales with the third power of the magnetic energy. When the generation speed exceeds the CME limit (regime II), the GW energy continues to increase without a corresponding increase of magnetic energy. In the early kinematic phase, the GW energy spectrum (per linear wavenumber interval) has opposite slopes in both regimes and is characterized by an inertial range spectrum in regime I and a white noise spectrum in regime II. The occurrence of these two slopes is shown to be a generic consequence of a nearly monochromatic exponential growth of the magnetic field. The resulting GW energy is found to be proportional to the fifth power of the limiting CME speed and the first power of the generation speed. △ Less

Submitted 22 April, 2021; v1 submitted 20 January, 2021; originally announced January 2021.

Comments: 15 pages, 15 figures, 4 tables, ApJ

Report number: NORDITA-2021-005

Journal ref: Astrophys. J. 911, 110 (2021)

Circular Polarization of Gravitational Waves from Early-Universe Helical Turbulence

Authors: Tina Kahniashvili, Axel Brandenburg, Grigol Gogoberidze, Sayan Mandal, Alberto Roper Pol

Abstract: We perform direct numerical simulations to compute the net circular polarization of gravitational waves from helical (chiral) turbulent sources in the early Universe for a variety of initial conditions, including driven (stationary) and decaying turbulence. We investigate the resulting gravitational wave signal assuming different turbulent geneses such as magnetically or kinetically driven cases.… ▽ More We perform direct numerical simulations to compute the net circular polarization of gravitational waves from helical (chiral) turbulent sources in the early Universe for a variety of initial conditions, including driven (stationary) and decaying turbulence. We investigate the resulting gravitational wave signal assuming different turbulent geneses such as magnetically or kinetically driven cases. Under realistic physical conditions in the early Universe we compute numerically the wave number-dependent polarization degree of the gravitational waves. We find that the spectral polarization degree strongly depends on the initial conditions. The peak of the spectral polarization degree occurs at twice the typical wavenumber of the source, as expected, and for fully helical decaying turbulence, it reaches its maximum of nearly 100\% {\it only} at the peak. We determine the temporal evolution of the turbulent sources as well as the resulting gravitational waves, showing that the dominant contribution to their spectral energy density happens shortly after the activation of the source. Only through an artificially prolonged decay of the turbulence can further increase of the gravitational wave amplitude be achieved. We estimate the detection prospects for the net polarization, arguing that its detection contains {\it clean} information (including the generation mechanisms, time, and strength) about the sources of possible parity violations in the early Universe. △ Less

Submitted 8 February, 2021; v1 submitted 11 November, 2020; originally announced November 2020.

Comments: 9 pages, 5 figures; discussions added; references updated and added; conclusions unchanged; Phys. Rev. Research in press

Report number: NORDITA-2020-102

Journal ref: Phys. Rev. Research 3, 013193 (2021)

Turbulent radiative diffusion and turbulent Newtonian cooling

Authors: Axel Brandenburg, Upasana Das

Abstract: Radiation transport plays important roles in stellar atmospheres, but the effects of turbulence are being obscured by other effects such as stratification. Using radiative hydrodynamic simulations of forced turbulence, we determine the decay rates of sinusoidal large-scale temperature perturbations of different wavenumbers in the optically thick and thin regimes. Increasing the wavenumber increase… ▽ More Radiation transport plays important roles in stellar atmospheres, but the effects of turbulence are being obscured by other effects such as stratification. Using radiative hydrodynamic simulations of forced turbulence, we determine the decay rates of sinusoidal large-scale temperature perturbations of different wavenumbers in the optically thick and thin regimes. Increasing the wavenumber increases the rate of decay in both regimes, but this effect is much weaker than for the usual turbulent diffusion of passive scalars, where the increase is quadratic for small wavenumbers. The turbulent decay is well described by an enhanced Newtonian cooling process in the optically thin limit, which is found to show a weak increase proportional to the square root of the wavenumber. In the optically thick limit, the increase in turbulent decay is somewhat steeper for wavenumbers below the energy-carrying wavenumber of the turbulence, but levels off toward larger wavenumbers. In the presence of turbulence, the typical cooling time is comparable to the turbulent turnover time. We observe that the temperature takes a long time to reach equilibrium in both the optically thin and thick cases, but in the former, the temperature retains smaller scale structures for longer. △ Less

Submitted 28 September, 2021; v1 submitted 13 October, 2020; originally announced October 2020.

Comments: 10 pages, 7 figures, 1 table, Physics of Fluids

Report number: NORDITA-2020-095

Journal ref: Phys. Fluids 33, 095125 (2021)

The Pencil Code, a modular MPI code for partial differential equations and particles: multipurpose and multiuser-maintained

Authors: A. Brandenburg, A. Johansen, P. A. Bourdin, W. Dobler, W. Lyra, M. Rheinhardt, S. Bingert, N. E. L. Haugen, A. Mee, F. Gent, N. Babkovskaia, C. -C. Yang, T. Heinemann, B. Dintrans, D. Mitra, S. Candelaresi, J. Warnecke, P. J. Käpylä, A. Schreiber, P. Chatterjee, M. J. Käpylä, X. -Y. Li, J. Krüger, J. R. Aarnes, G. R. Sarson , et al. (12 additional authors not shown)

Abstract: The Pencil Code is a highly modular physics-oriented simulation code that can be adapted to a wide range of applications. It is primarily designed to solve partial differential equations (PDEs) of compressible hydrodynamics and has lots of add-ons ranging from astrophysical magnetohydrodynamics (MHD) to meteorological cloud microphysics and engineering applications in combustion. Nevertheless, the… ▽ More The Pencil Code is a highly modular physics-oriented simulation code that can be adapted to a wide range of applications. It is primarily designed to solve partial differential equations (PDEs) of compressible hydrodynamics and has lots of add-ons ranging from astrophysical magnetohydrodynamics (MHD) to meteorological cloud microphysics and engineering applications in combustion. Nevertheless, the framework is general and can also be applied to situations not related to hydrodynamics or even PDEs, for example when just the message passing interface or input/output strategies of the code are to be used. The code can also evolve Lagrangian (inertial and noninertial) particles, their coagulation and condensation, as well as their interaction with the fluid. △ Less

Submitted 17 September, 2020; originally announced September 2020.

Comments: 7 pages, submitted to the Journal for Open Source Software (JOSS)

Report number: NORDITA-2020-087

Journal ref: Journal of Open Source Software 6, 2807 (2021)

Hall cascade with fractional magnetic helicity in neutron star crusts

Authors: Axel Brandenburg

Abstract: The ohmic decay of magnetic fields in the crusts of neutron stars is generally believed to be governed by Hall drift which leads to what is known as a Hall cascade. Here we show that helical and fractionally helical magnetic fields undergo strong inverse cascading like in magnetohydrodynamics (MHD), but the magnetic energy decays more slowly with time $t$: $\propto t^{-2/5}$ instead of… ▽ More The ohmic decay of magnetic fields in the crusts of neutron stars is generally believed to be governed by Hall drift which leads to what is known as a Hall cascade. Here we show that helical and fractionally helical magnetic fields undergo strong inverse cascading like in magnetohydrodynamics (MHD), but the magnetic energy decays more slowly with time $t$: $\propto t^{-2/5}$ instead of $\propto t^{-2/3}$ in MHD. Even for a nonhelical magnetic field there is a certain degree of inverse cascading for sufficiently strong magnetic fields. The inertial range scaling with wavenumber $k$ is compatible with earlier findings for the forced Hall cascade, i.e., proportional to $k^{-7/3}$, but in the decaying cases, the subinertial range spectrum steepens to a novel $k^5$ slope instead of the $k^4$ slope in MHD. The energy of the large-scale magnetic field can increase quadratically in time through inverse cascading. For helical fields, the energy dissipation is found to be inversely proportional to the large-scale magnetic field and proportional to the fifth power of the root-mean square (rms) magnetic field. For neutron star conditions with an rms magnetic field of a few times $10^{14}\,$G, the large-scale magnetic field might only be $10^{11}\,$G, while still producing magnetic dissipation of $10^{33}\,$erg$\,$s$^{-1}$ for thousands of years, which could manifest itself through X-ray emission. Finally, it is shown that the conclusions from local unstratified models agree rather well with those from stratified models with boundaries. △ Less

Submitted 3 August, 2020; v1 submitted 23 June, 2020; originally announced June 2020.

Comments: 14 pages, 17 figures, 5 tables, ApJ, in press

Report number: NORDITA-2020-063

Journal ref: Astrophys. J. 901, 18 (2020)

On the existence of shear-current effects in magnetized burgulence

Authors: Maarit J. Käpylä, Javier Álvarez Vizoso, Matthias Rheinhardt, Axel Brandenburg, Nishant K. Singh

Abstract: The possibility of explaining shear flow dynamos by a magnetic shear-current (MSC) effect is examined via numerical simulations. Our primary diagnostics is the determination of the turbulent magnetic diffusivity tensor $\boldsymbolη$. In our setup, a negative sign of its component $η_{yx}$ is necessary for coherent dynamo action by the SC effect. To be able to measure turbulent transport coefficie… ▽ More The possibility of explaining shear flow dynamos by a magnetic shear-current (MSC) effect is examined via numerical simulations. Our primary diagnostics is the determination of the turbulent magnetic diffusivity tensor $\boldsymbolη$. In our setup, a negative sign of its component $η_{yx}$ is necessary for coherent dynamo action by the SC effect. To be able to measure turbulent transport coefficients from systems with magnetic background turbulence, we present an extension of the test-field method (TFM), applicable to our setup where the pressure gradient is dropped from the momentum equation: the nonlinear TFM (NLTFM). Our momentum equation is related to Burgers' equation and the resulting flows are referred to as magnetized burgulence. We use both stochastic kinetic and magnetic forcings to mimic cases without and with simultaneous small-scale dynamo action (SSD). When we force only kinetically, negative $η_{yx}$ are obtained with exponential growth in both the radial and azimuthal mean magnetic field components. Using isotropic magnetic forcing, the field growth is no longer exponential, while NLTFM yields positive $η_{yx}$. By employing an alternative forcing from which wavevectors having small components are removed, the exponential growth is recovered, but the NLTFM results do not change significantly. Analyzing the dynamo excitation conditions for the coherent SC and incoherent $α$ and SC effects shows that the incoherent effects are the main drivers of the dynamo in the majority of cases. We find no evidence for MSC-effect-driven dynamos in our simulations. △ Less

Submitted 12 January, 2021; v1 submitted 10 June, 2020; originally announced June 2020.

Comments: 21 pages, 11 figures, published in Astrophysical Journal

Report number: NORDITA-2020-052

Journal ref: Astrophys. J. 905, 179 (2020)

The effect of a dynamo-generated field on the Parker wind

Authors: P. Jakab, A. Brandenburg

Abstract: Stellar winds are an integral part of the underlying dynamo, the motor of stellar activity. The wind controls the star's angular momentum loss, which depends on the magnetic field geometry which varies significantly in time and latitude. Here we study basic properties of a self-consistent model that includes simple representations of both the global stellar dynamo in a spherical shell and the exte… ▽ More Stellar winds are an integral part of the underlying dynamo, the motor of stellar activity. The wind controls the star's angular momentum loss, which depends on the magnetic field geometry which varies significantly in time and latitude. Here we study basic properties of a self-consistent model that includes simple representations of both the global stellar dynamo in a spherical shell and the exterior in which the wind accelerates and becomes supersonic. We numerically solve an axisymmetric mean-field model for the induction, momentum, and continuity equations using an isothermal equation of state. The model allows for the simultaneous generation of a mean magnetic field and the development of a Parker wind. The resulting flow is transonic at the critical point, which we arrange to be between the inner and outer radii of the model. The boundary conditions are assumed to be such that the magnetic field is antisymmetric about the equator, i.e., dipolar. At the solar rotation rate, the dynamo is oscillatory and of $α^2$ type. In most of the domain, the magnetic field corresponds to that of a split monopole. The magnetic energy flux is largest between the stellar surface and the critical point. The angular momentum flux is highly variable in time and can reach negative values, especially at midlatitudes. At rapid rotation of up to 50 times the solar value, most of the magnetic field is lost along the axis within the inner tangential cylinder of the model. The model reveals unexpected features that are not generally anticipated from models that are designed to reproduce the solar wind: highly variable angular momentum fluxes even from just an $α^2$ dynamo in the star. A major caveat of our isothermal models with a magnetic field produced by a dynamo is the difficulty to reach small enough plasma betas without the dynamo itself becoming unrealistically strong inside the star. △ Less

Submitted 15 March, 2021; v1 submitted 3 June, 2020; originally announced June 2020.

Comments: 15 pages, 18 figures, 1 table, A&A, published

Report number: NORDITA-2020-050

Journal ref: Astron. Astrophys. 647, A18 (2021)

On the measurement of handedness in Fermi Large Area Telescope data

Authors: Julia Asplund, Guðlaugur Jóhannesson, Axel Brandenburg

Abstract: A handedness in the arrival directions of high-energy photons from outside our Galaxy can be related to the helicity of an intergalactic magnetic field. Previous estimates by arXiv:1310.4826 and arXiv:1412.3171 showed a hint of a signal present in the photons observed by the Fermi Large Area Telescope (LAT). An update on the measurement of handedness in Fermi-LAT data is presented using more than… ▽ More A handedness in the arrival directions of high-energy photons from outside our Galaxy can be related to the helicity of an intergalactic magnetic field. Previous estimates by arXiv:1310.4826 and arXiv:1412.3171 showed a hint of a signal present in the photons observed by the Fermi Large Area Telescope (LAT). An update on the measurement of handedness in Fermi-LAT data is presented using more than 10 years of observations. Simulations are performed to study the uncertainty of the measurements, taking into account the structure of the exposure caused by the energy-dependent instrument response and its observing profile, as well as the background from the interstellar medium. The simulations are required to accurately estimate the uncertainty and to show that previously the uncertainty was significantly underestimated. The apparent signal in the earlier analysis of {\em Fermi}-LAT data is rendered non-significant. △ Less

Submitted 26 May, 2020; originally announced May 2020.

Comments: 11 pages, 7 figures. Accepted for publication in ApJ

Report number: NORDITA-2020-048

Journal ref: Astrophys. J. 898, 124 (2020)

Primordial magnetic helicity evolution with a homogeneous magnetic field from inflation

Authors: Axel Brandenburg, Ruth Durrer, Yiwen Huang, Tina Kahniashvili, Sayan Mandal, Shinji Mukohyama

Abstract: Motivated by a scenario of magnetogenesis in which a homogeneous magnetic field is generated during inflation, we study the magnetohydrodynamic evolution of the primordial plasma motions for two kinds of initial conditions -- (i) a spatially homogeneous field with an unlimited correlation length, and (ii) a zero flux scale-invariant statistically homogeneous magnetic field. In both cases, we apply… ▽ More Motivated by a scenario of magnetogenesis in which a homogeneous magnetic field is generated during inflation, we study the magnetohydrodynamic evolution of the primordial plasma motions for two kinds of initial conditions -- (i) a spatially homogeneous field with an unlimited correlation length, and (ii) a zero flux scale-invariant statistically homogeneous magnetic field. In both cases, we apply, for a short initial time interval, monochromatic forcing at a certain wave number so that the correlation length is finite, but much smaller than the typical length scale of turbulence. In particular, we investigate the decay of nonhelical and helical hydromagnetic turbulence. We show that, in the presence of a homogeneous magnetic field, the decay of helical and nonhelical small-scale fields can occur rapidly. This is a special property of a system with a perfectly homogeneous magnetic field, which is sometimes considered as a local approximation to a slowly varying background field. It can never change and acts as an imposed magnetic field. This is in a sharp contrast to the case of a statistically homogeneous magnetic field, where we recover familiar decay properties: a much slower decay of magnetic energy and a faster growth of the correlation length, especially in the case with magnetic felicity. The result suggests that a homogeneous magnetic field, if generated during inflation, should persist under the influence of small-scale fields and could be the origin of the large-scale magnetic field in the Universe. △ Less

Submitted 25 July, 2020; v1 submitted 13 May, 2020; originally announced May 2020.

Comments: 16 pages, 7 figures, 1 table. Journal version

Report number: NORDITA-2020-044, YITP-20-47, IPMU20-0042

Journal ref: Phys. Rev. D 102, 023536 (2020)

Magnetic field evolution in solar-type stars

Authors: Axel Brandenburg

Abstract: We discuss selected aspects regarding the magnetic field evolution of solar-type stars. Most of the stars with activity cycles are in the range where the normalized chromospheric Calcium emission increases linearly with the inverse Rossby number. For Rossby numbers below about a quarter of the solar value, the activity saturates and no cycles have been found. For Rossby numbers above the solar val… ▽ More We discuss selected aspects regarding the magnetic field evolution of solar-type stars. Most of the stars with activity cycles are in the range where the normalized chromospheric Calcium emission increases linearly with the inverse Rossby number. For Rossby numbers below about a quarter of the solar value, the activity saturates and no cycles have been found. For Rossby numbers above the solar value, again no activity cycles have been found, but now the activity goes up again for a major fraction of the stars. Rapidly rotating stars show nonaxisymmetric large-scale magnetic fields, but there is disagreement between models and observations regarding the actual value of the Rossby number where this happens. We also discuss the prospects of detecting the sign of magnetic helicity using various linear polarization techniques both at the stellar surface using the parity-odd contribution to linear polarization and above the surface using Faraday rotation. △ Less

Submitted 31 March, 2020; originally announced April 2020.

Comments: 13 pages, 3 figures, 3 tables, proceedings of IAU symposium 354, Solar and Stellar Magnetic Fields: Origins and Manifestations

Report number: NORDITA-2020-029

Journal ref: Proc. IAU 15 (2019) 169-180

Hemispheric handedness in the Galactic synchrotron polarization foreground

Authors: Axel Brandenburg, Marcus Brüggen

Abstract: The large-scale magnetic field of the Milky Way is thought to be created by an alpha-Omega dynamo, which implies that it should have opposite handedness North and South of the Galactic midplane. Here we attempt to detect a variation in handedness using polarization data from the Wilkinson Microwave Anisotropy Probe. Previous analyzes of the parity-even and parity-odd parts of linear polarization o… ▽ More The large-scale magnetic field of the Milky Way is thought to be created by an alpha-Omega dynamo, which implies that it should have opposite handedness North and South of the Galactic midplane. Here we attempt to detect a variation in handedness using polarization data from the Wilkinson Microwave Anisotropy Probe. Previous analyzes of the parity-even and parity-odd parts of linear polarization of the global dust and synchrotron emission have focused on quadratic correlations in spectral space of, and between, these two components. Here, by contrast, we analyze the parity-odd polarization itself and show that it has, on average, opposite signs in Northern and Southern Galactic hemispheres. Comparison with a Galactic mean-field dynamo model shows broad qualitative agreement and reveals that the sign of the observed hemispheric dependence of the azimuthally averaged parity-odd polarization is not determined by the sign of alpha, but by the sense of differential rotation. △ Less

Submitted 24 May, 2020; v1 submitted 29 March, 2020; originally announced March 2020.

Comments: 7 pages, 4 figures, 1 table, submitted to ApJL

Report number: NORDITA-2020-027

Journal ref: Astrophys. J. Lett. 896, L14 (2020)

Application of a helicity proxy to edge-on galaxies

Authors: Axel Brandenburg, Ray S. Furuya

Abstract: We study the prospects of detecting magnetic helicity in galaxies by observing the dust polarization of the edge-on galaxy NGC 891. Our numerical results of mean-field dynamo calculations show that there should be a large-scale component of the rotationally invariant parity-odd B polarization that we predict to be negative in the first and third quadrants, and positive in the second and fourth qua… ▽ More We study the prospects of detecting magnetic helicity in galaxies by observing the dust polarization of the edge-on galaxy NGC 891. Our numerical results of mean-field dynamo calculations show that there should be a large-scale component of the rotationally invariant parity-odd B polarization that we predict to be negative in the first and third quadrants, and positive in the second and fourth quadrants. The large-scale parity-even E polarization is predicted to be negative near the axis and positive further away in the outskirts. These properties are shown to be mostly a consequence of the magnetic field being azimuthal and the polarized intensity being maximum at the center of the galaxy and are not a signature of magnetic helicity. △ Less

Submitted 17 June, 2020; v1 submitted 16 March, 2020; originally announced March 2020.

Comments: 12 pages, 11 figures, 2 tables, MNRAS, in press

Report number: NORDITA-2020-025

Journal ref: Mon. Not. Roy. Astron. Soc. 496, 4749-4759 (2020)