TY - JOUR A1 - Thomas, Timon A1 - Pfrommer, Christoph A1 - Pakmor, Rüdiger T1 - A finite volume method for two-moment cosmic ray hydrodynamics on a moving mesh JF - Monthly notices of the Royal Astronomical Society N2 - We present a new numerical algorithm to solve the recently derived equations of two-moment cosmic ray hydrodynamics (CRHD). The algorithm is implemented as a module in the moving mesh AREPO code. Therein, the anisotropic transport of cosmic rays (CRs) along magnetic field lines is discretized using a path-conservative finite volume method on the unstructured time-dependent Voronoi mesh of AREPO. The interaction of CRs and gyroresonant Alfven waves is described by short time-scale source terms in the CRHD equations. We employ a custom-made semi-implicit adaptive time stepping source term integrator to accurately integrate this interaction on the small light-crossing time of the anisotropic transport step. Both the transport and the source term integration step are separated from the evolution of the magnetohydrodynamical equations using an operator split approach. The new algorithm is tested with a variety of test problems, including shock tubes, a perpendicular magnetized discontinuity, the hydrodynamic response to a CR overpressure, CR acceleration of a warm cloud, and a CR blast wave, which demonstrate that the coupling between CR and magnetohydrodynamics is robust and accurate. We demonstrate the numerical convergence of the presented scheme using new linear and non-linear analytic solutions. KW - hydrodynamics KW - MHD KW - methods: numerical KW - cosmic rays Y1 - 2021 U6 - https://doi.org/10.1093/mnras/stab397 SN - 0035-8711 SN - 1365-2966 VL - 503 IS - 2 SP - 2242 EP - 2264 PB - Oxford University Press CY - Oxford ER - TY - JOUR A1 - Sparre, Martin A1 - Whittingham, Joseph A1 - Damle, Mitali A1 - Hani, Maan H. A1 - Richter, Philipp A1 - Ellison, Sara L. A1 - Pfrommer, Christoph A1 - Vogelsberger, Mark T1 - Gas flows in galaxy mergers BT - supersonic turbulence in bridges, accretion from the circumgalactic medium, and metallicity dilution JF - Monthly notices of the Royal Astronomical Society N2 - In major galaxy mergers, the orbits of stars are violently perturbed, and gas is torqued to the centre, diluting the gas metallicity and igniting a starburst. In this paper, we study the gas dynamics in and around merging galaxies using a series of cosmological magnetohydrodynamical zoom-in simulations. We find that the gas bridge connecting the merging galaxies pre-coalescence is dominated by turbulent pressure, with turbulent Mach numbers peaking at values of 1.6-3.3. This implies that bridges are dominated by supersonic turbulence, and are thus ideal candidates for studying the impact of extreme environments on star formation. We also find that gas accreted from the circumgalactic medium (CGM) during the merger significantly contributes (27-51 percent) to the star formation rate (SFR) at the time of coalescence and drives the subsequent reignition of star formation in the merger remnant. Indeed, 19-53 percent of the SFR at z = 0 originates from gas belonging to the CGM prior the merger. Finally, we investigate the origin of the metallicity-diluted gas at the centre of merging galaxies. We show that this gas is rapidly accreted on to the Galactic Centre with a time-scale much shorter than that of normal star-forming galaxies. This explains why coalescing galaxies are not well-captured by the fundamental metallicity relation. KW - MHD KW - methods: numerical KW - galaxies: interactions KW - galaxies: starburst Y1 - 2021 U6 - https://doi.org/10.1093/mnras/stab3171 SN - 1365-2966 VL - 509 IS - 2 SP - 2720 EP - 2735 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Werhahn, Maria A1 - Pfrommer, Christoph A1 - Girichidis, Philipp T1 - Cosmic rays and non-thermal emission in simulated galaxies - III. Probing cosmic-ray calorimetry with radio spectra and the FIR-radio correlation JF - Monthly notices of the Royal Astronomical Society N2 - An extinction-free estimator of the star formation rate (SFR) of galaxies is critical for understanding the high-redshift universe. To this end, the nearly linear, tight correlation of far-infrared (FIR), and radio luminosity of star-forming galaxies is widely used. While the FIR is linked to massive star formation, which also generates shock-accelerated cosmic-ray (CR) electrons and radio synchrotron emission, a detailed understanding of the underlying physics is still lacking. Hence, we perform three-dimensional magnetohydrodynamical (MHD) simulations of isolated galaxies over a broad range of halo masses and SFRs using the moving-mesh code AREPO, and evolve the CR proton energy density self-consistently. In post-processing, we calculate the steady-state spectra of primary, shock-accelerated and secondary CR electrons, which result from hadronic CR proton interactions with the interstellar medium. The resulting total radio luminosities correlate with the FIR luminosities as observed and are dominated by primary CR electrons if we account for anisotropic CR diffusion. The increasing contribution of secondary emission up to 30 per cent in starbursts is compensated by the larger bremsstrahlung and Coulomb losses. CR electrons are in the calorimetric limit and lose most of their energy through inverse Compton interactions with star light and cosmic microwave background (CMB) photons while less energy is converted into synchrotron emission. This implies steep steady-state synchrotron spectra in starbursts. Interestingly, we find that thermal free-free emission flattens the total radio spectra at high radio frequencies and reconciles calorimetric theory with observations while free-free absorption explains the observed low-frequency flattening towards the central regions of starbursts. KW - MHD KW - methods: numerical KW - cosmic rays KW - galaxies: magnetic fields KW - galaxies: starburst KW - radio continuum: galaxies Y1 - 2021 U6 - https://doi.org/10.1093/mnras/stab2535 SN - 0035-8711 SN - 1365-2966 VL - 508 IS - 3 SP - 4072 EP - 4095 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Rüdiger, Günther A1 - Schultz, Manfred T1 - On the toroidal-velocity antidynamo theorem under the presence of nonuniform electric conductivity JF - Astronomische Nachrichten = Astronomical notes N2 - Laminar electrically conducting Couette flows with the hydrodynamically stable quasi-Keplerian rotation profile and nonuniform conductivity are probed for dynamo instability. In spherical geometry, the equations for the poloidal and the toroidal field components completely decouple, resulting in free decay, regardless of the spatial distribution of the electric conductivity. In cylindrical geometry the poloidal and toroidal components do not decouple, but here also we do not find dynamo excitations for the cases that the electric conductivity only depends on the radius or - much more complex- that it only depends on the azimuthal or the axial coordinate. The transformation of the plane-flow dynamo model of Busse and Wicht (1992) to cylindrical or spherical geometry therefore fails. It is also shown that even the inclusion of axial flows of both directions does not support the dynamo mechanism. The Elsasser toroidal-velocity antidynamo theorem, according to which dynamos without any radial velocity component cannot work, is thus not softened by nonuniform conductivity distributions. KW - antidynamo theorem KW - MHD KW - Taylor-Couette flow Y1 - 2022 U6 - https://doi.org/10.1002/asna.20224011 SN - 0004-6337 SN - 1521-3994 VL - 343 IS - 5 PB - Wiley-VCH CY - Weinheim ER - TY - JOUR A1 - Werhahn, Maria A1 - Pfrommer, Christoph A1 - Girichidis, Philipp A1 - Puchwein, Ewald A1 - Pakmor, Rüdiger T1 - Cosmic rays and non-thermal emission in simulated galaxies BT - I. Electron and proton spectra compared to Voyager-1 data JF - Monthly notices of the Royal Astronomical Society N2 - Current-day cosmic ray (CR) propagation studies use static Milky Way models and fit parametrized source distributions to data. Instead, we use three-dimensional magnetohydrodynamic (MHD) simulations of isolated galaxies with the moving-mesh code arepo that self-consistently accounts for hydrodynamic effects of CR protons. In post-processing, we calculate their steady-state spectra, taking into account all relevant loss processes. We show that this steady-state assumption is well justified in the disc and generally for regions that emit non-thermal radio and gamma rays. Additionally, we model the spectra of primary electrons, accelerated by supernova remnants, and secondary electrons and positrons produced in hadronic CR proton interactions with the gas. We find that proton spectra above 10 GeV only weakly depend on galactic radius, while they acquire a radial dependence at lower energies due to Coulomb interactions. Radiative losses steepen the spectra of primary CR electrons in the central galactic regions, while diffusive losses dominate in the outskirts. Secondary electrons exhibit a steeper spectrum than primaries because they originate from the transported steeper CR proton spectra. Consistent with Voyager-1 and AMS-02 data, our models (i) show a turnover of proton spectra below GeV energies due to Coulomb interactions so that electrons start to dominate the total particle spectra and (ii) match the shape of the positron fraction up to 10 GeV. We conclude that our steady-state CR modelling in MHD CR galaxy simulations is sufficiently realistic to capture the dominant transport effects shaping their spectra, arguing for a full MHD treatment to accurately model CR transport in the future. KW - astroparticle physics KW - MHD KW - methods: numerical KW - cosmic rays KW - local KW - interstellar matter Y1 - 2021 U6 - https://doi.org/10.1093/mnras/stab1324 SN - 0035-8711 SN - 1365-2966 VL - 505 IS - 3 SP - 3273 EP - 3294 PB - Oxford University Press CY - Oxford ER - TY - JOUR A1 - Whittingham, Joseph A1 - Sparre, Martin A1 - Pfrommer, Christoph A1 - Pakmor, Rüdiger T1 - The impact of magnetic fields on cosmological galaxy mergers BT - I. Reshaping gas and stellar discs JF - Monthly notices of the Royal Astronomical Society N2 - Mergers play an important role in galaxy evolution. In particular, major mergers are able to have a transformative effect on galaxy morphology. In this paper, we investigate the role of magnetic fields in gas-rich major mergers. To this end, we run a series of high-resolution magnetohydrodynamic (MHD) zoom-in simulations with the moving-mesh code arepo and compare the outcome with hydrodynamic simulations run from the same initial conditions. This is the first time that the effect of magnetic fields in major mergers has been investigated in a cosmologically consistent manner. In contrast to previous non-cosmological simulations, we find that the inclusion of magnetic fields has a substantial impact on the production of the merger remnant. Whilst magnetic fields do not strongly affect global properties, such as the star formation history, they are able to significantly influence structural properties. Indeed, MHD simulations consistently form remnants with extended discs and well-developed spiral structure, whilst hydrodynamic simulations form more compact remnants that display distinctive ring morphology. We support this work with a resolution study and show that whilst global properties are broadly converged across resolution and physics models, morphological differences only develop given sufficient resolution. We argue that this is due to the more efficient excitement of a small-scale dynamo in higher resolution simulations, resulting in a more strongly amplified field that is better able to influence gas dynamics. KW - MHD KW - methods: numerical KW - galaxies: interactions KW - galaxies: magnetic KW - fields Y1 - 2021 U6 - https://doi.org/10.1093/mnras/stab1425 SN - 0035-8711 SN - 1365-2966 VL - 506 IS - 1 SP - 229 EP - 255 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Rüdiger, Günther A1 - Schultz, Manfred T1 - Large-scale dynamo action of magnetized Taylor-Couette flows JF - Monthly notices of the Royal Astronomical Society N2 - A conducting Taylor-Couette flow with quasi-Keplerian rotation law containing a toroidal magnetic field serves as a mean-field dynamo model of the Tayler-Spruit type. The flows are unstable against non-axisymmetric perturbations which form electromotive forces defining a effect and eddy diffusivity. If both degenerated modes with m = +/- 1 are excited with the same power then the global a effect vanishes and a dynamo cannot work. It is shown, however, that the Tayler instability produces finite alpha effects if only an isolated mode is considered but this intrinsic helicity of the single-mode is too low for an alpha(2) dynamo. Moreover, an alpha Omega dynamo model with quasi-Keplerian rotation requires a minimum magnetic Reynolds number of rotation of Rm similar or equal to 2000 to work. Whether it really works depends on assumptions about the turbulence energy. For a steeper-than-quadratic dependence of the turbulence intensity on the magnetic field, however, dynamos are only excited if the resulting magnetic eddy diffusivity approximates its microscopic value, eta(T) similar or equal to eta. By basically lower or larger eddy diffusivities the dynamo instability is suppressed. KW - dynamo KW - instabilities KW - MHD KW - magnetic fields Y1 - 2020 U6 - https://doi.org/10.1093/mnras/staa293 SN - 0035-8711 SN - 1365-2966 VL - 493 IS - 1 SP - 1249 EP - 1260 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Meyer, Dominique M.-A. T1 - On the bipolarity of Wolf-Rayet nebulae JF - Monthly notices of the Royal Astronomical Society N2 - Wolf-Rayet stars are amongst the rarest but also most intriguing massive stars. Their extreme stellar winds induce famous multiwavelength circumstellar gas nebulae of various morphologies, spanning from circles and rings to bipolar shapes. This study is devoted to the investigation of the formation of young, asymmetric Wolf-Rayet gas nebulae and we present a 2.5-dimensional magneto-hydrodynamical toy model for the simulation of Wolf-Rayet gas nebulae generated by wind-wind interaction. Our method accounts for stellar wind asymmetries, rotation, magnetization, evolution, and mixing of materials. It is found that the morphology of the Wolf-Rayet nebulae of blue supergiant ancestors is tightly related to the wind geometry and to the stellar phase transition time interval, generating either a broadened peanut-like or a collimated jet-like gas nebula. Radiative transfer calculations of our Wolf-Rayet nebulae for dust infrared emission at 24 mu m show that the projected diffuse emission can appear as oblate, bipolar, ellipsoidal, or ring structures. Important projection effects are at work in shaping observed Wolf-Rayet nebulae. This might call a revision of the various classifications of Wolf-Rayet shells, which are mostly based on their observed shape. Particularly, our models question the possibility of producing pre-Wolf-Rayet wind asymmetries, responsible for bipolar nebulae like NGC 6888, within the single red supergiant evolution channel scenario. We propose that bipolar Wolf-Rayet nebulae can only be formed within the red supergiant scenario by multiple/merged massive stellar systems, or by single high-mass stars undergoing additional, e.g. blue supergiant, evolutionary stages prior to the Wolf-Rayet phase. KW - MHD KW - radiative transfer KW - circumstellar matter KW - stars: massive KW - stars: KW - Wolf-Rayet Y1 - 2021 U6 - https://doi.org/10.1093/mnras/stab2426 SN - 0035-8711 SN - 1365-2966 VL - 507 IS - 4 SP - 4697 EP - 4714 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Lebiga, O. A1 - Santos-Lima, Reinaldo A1 - Yan, Huirong T1 - Kinetic-MHD simulations of gyroresonance instability driven by CR pressure anisotropy JF - Monthly notices of the Royal Astronomical Society N2 - The transport of cosmic rays (CRs) is crucial for the understanding of almost all high-energy phenomena. Both pre-existing large-scale magnetohydrodynamic (MHD) turbulence and locally generated turbulence through plasma instabilities are important for the CR propagation in astrophysical media. The potential role of the resonant instability triggered by CR pressure anisotropy to regulate the parallel spatial diffusion of low-energy CRs (less than or similar to 100 GeV) in the interstellar and intracluster medium of galaxies has been shown in previous theoretical works. This work aims to study the gyroresonance instability via direct numerical simulations, in order to access quantitatively the wave-particle scattering rates. For this, we employ a 1D PIC-MHD code to follow the growth and saturation of the gyroresonance instability. We extract from the simulations the pitch-angle diffusion coefficient D-mu mu produced by the instability during the linear and saturation phases, and a very good agreement (within a factor of 3) is found with the values predicted by the quasi-linear theory (QLT). Our results support the applicability of the QLT for modelling the scattering of low-energy CRs by the gyroresonance instability in the complex interplay between this instability and the large-scale MHD turbulence. KW - MHD KW - plasmas KW - turbulence KW - cosmic rays Y1 - 2018 U6 - https://doi.org/10.1093/mnras/sty309 SN - 0035-8711 SN - 1365-2966 VL - 476 IS - 2 SP - 2779 EP - 2791 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Santos de Lima, Reinaldo A1 - Yan, Huirong A1 - de Gouveia Dal Pino, E. M. A1 - Lazarian, A. T1 - Limits on the ion temperature anisotropy in the turbulent intracluster medium JF - Monthly notices of the Royal Astronomical Society N2 - Turbulence in the weakly collisional intracluster medium (ICM) of galaxies is able to generate strong thermal velocity anisotropies in the ions (with respect to the local magnetic field direction), if the magnetic moment of the particles is conserved in the absence of Coulomb collisions. In this scenario, the anisotropic pressure magnetohydrodynamic (AMHD) turbulence shows a very different statistical behaviour from the standard MHD one and is unable to amplify seed magnetic fields. This is in contrast to previous cosmological MHD simulations that are successful in explaining the observed magnetic fields in the ICM. On the other hand, temperature anisotropies can also drive plasma instabilities that can relax the anisotropy. This work aims to compare the relaxation rate with the growth rate of the anisotropies driven by the turbulence. We employ quasi-linear theory to estimate the ion scattering rate resulting from the parallel firehose, mirror and ion-cyclotron instabilities, for a set of plasma parameters resulting from AMHD simulations of the turbulent ICM. We show that the ICM turbulence can sustain only anisotropy levels very close to the instability thresholds. We argue that the AMHD model that bounds the anisotropies at the marginal stability levels can describe the Alfv,nic turbulence cascade in the ICM. KW - MHD KW - plasmas KW - turbulence KW - galaxies: clusters: intracluster medium Y1 - 2016 U6 - https://doi.org/10.1093/mnras/stw1079 SN - 0035-8711 SN - 1365-2966 VL - 460 SP - 2492 EP - 2504 PB - Oxford Univ. Press CY - Oxford ER -