Refine
Has Fulltext
- no (31)
Year of publication
Document Type
- Article (31) (remove)
Language
- English (31)
Is part of the Bibliography
- yes (31)
Keywords
- methods: numerical (31) (remove)
Institute
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.
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.
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.
Context. For magnetically driven events, the magnetic energy of the system is the prime energy reservoir that fuels the dynamical evolution. In the solar context, the free energy (i.e., the energy in excess of the potential field energy) is one of the main indicators used in space weather forecasts to predict the eruptivity of active regions. A trustworthy estimation of the magnetic energy is therefore needed in three-dimensional (3D) models of the solar atmosphere, e. g., in coronal fields reconstructions or numerical simulations.
Aims. The expression of the energy of a system as the sum of its potential energy and its free energy (Thomson's theorem) is strictly valid when the magnetic field is exactly solenoidal. For numerical realizations on a discrete grid, this property may be only approximately fulfilled. We show that the imperfect solenoidality induces terms in the energy that can lead to misinterpreting the amount of free energy present in a magnetic configuration.
Methods. We consider a decomposition of the energy in solenoidal and nonsolenoidal parts which allows the unambiguous estimation of the nonsolenoidal contribution to the energy. We apply this decomposition to six typical cases broadly used in solar physics. We quantify to what extent the Thomson theorem is not satisfied when approximately solenoidal fields are used.
Results. The quantified errors on energy vary from negligible to significant errors, depending on the extent of the nonsolenoidal component of the field. We identify the main source of errors and analyze the implications of adding a variable amount of divergence to various solenoidal fields. Finally, we present pathological unphysical situations where the estimated free energy would appear to be negative, as found in some previous works, and we identify the source of this error to be the presence of a finite divergence.
Conclusions. We provide a method of quantifying the effect of a finite divergence in numerical fields, together with detailed diagnostics of its sources. We also compare the efficiency of two divergence-cleaning techniques. These results are applicable to a broad range of numerical realizations of magnetic fields.
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.
Cosmic-ray hydrodynamics
(2019)
Star formation in galaxies appears to be self-regulated by energetic feedback processes. Among the most promising agents of feedback are cosmic rays (CRs), the relativistic ion population of interstellar and intergalactic plasmas. In these environments, energetic CRs are virtually collisionless and interact via collective phenomena mediated by kinetic-scale plasma waves and large-scale magnetic fields. The enormous separation of kinetic and global astrophysical scales requires a hydrodynamic description. Here, we develop a new macroscopic theory for CR transport in the self-confinement picture, which includes CR diffusion and streaming. The interaction between CRs and electromagnetic fields of Alfvenic turbulence provides the main source of CR scattering, and causes CRs to stream along the magnetic field with the Alfven velocity if resonant waves are sufficiently energetic. However, numerical simulations struggle to capture this effect with current transport formalisms and adopt regularization schemes to ensure numerical stability. We extent the theory by deriving an equation for the CRmomentum density along the mean magnetic field and include a transport equation for the Alfven-wave energy. We account for energy exchange of CRs and Alfven waves via the gyroresonant instability and include other wave damping mechanisms. Using numerical simulations, we demonstrate that our new theory enables stable, self-regulated CR transport. The theory is coupled to magnetohydrodynamics, conserves the total energy and momentum, and correctly recovers previous macroscopic CR transport formalisms in the steady-state flux limit. Because it is free of tunable parameters, it holds the promise to provide predictable simulations of CR feedback in galaxy formation.
We investigate the origin and physical properties of O vi absorbers at low redshift (z = 0.25) using a subset of cosmological, hydrodynamical simulations from the OverWhelmingly Large Simulations (OWLS) project. Intervening O vi absorbers are believed to trace shock-heated gas in the warm-hot intergalactic medium (WHIM) and may thus play a key role in the search for the missing baryons in the present-day Universe. When compared to observations, the predicted distributions of the different O vi line parameters (column density, Doppler parameter, rest equivalent width W-r) from our simulations exhibit a lack of strong O vi absorbers, a discrepancy that has also been found by Oppenheimer & Dave. This suggests that physical processes on subgrid scales (e.g. turbulence) may strongly influence the observed properties of O vi systems. We find that the intervening O vi absorption arises mainly in highly metal enriched (10-1 < Z/Z(circle dot) less than or similar to 1) gas at typical overdensities of 1 < /<<>> less than or similar to 102. One-third of the O vi absorbers in our simulation are found to trace gas at temperatures T < 105 K, while the rest arises in gas at higher temperatures, most of them around T = 105.3 +/- 0.5 K. These temperatures are much higher than inferred by Oppenheimer & Dave, probably because that work did not take the suppression of metal-line cooling by the photoionizing background radiation into account. While the O vi resides in a similar region of (, T)-space as much of the shock-heated baryonic matter, the vast majority of this gas has a lower metal content and does not give rise to detectable O vi absorption. As a consequence of the patchy metal distribution, O vi absorbers in our simulations trace only a very small fraction of the cosmic baryons (< 2 per cent) and the cosmic metals. Instead, these systems presumably trace previously shock-heated, metal-rich material from galactic winds that is now mixing with the ambient gas and cooling. The common approach of comparing O vi and H i column densities to estimate the physical conditions in intervening absorbers from QSO observations may be misleading, as most of the H i (and most of the gas mass) is not physically connected with the high-metallicity patches that give rise to the O vi absorption.
We investigate the physical state of H?i absorbing gas at low redshift (z = 0.25) using a subset of cosmological, hydrodynamic simulations from the OverWhelmingly Large Simulations project, focusing in particular on broad (bHI=40 km s-1) H?i Lya absorbers (BLAs), which are believed to originate in shock-heated gas in the warm-hot intergalactic medium (WHIM). Our fiducial model, which includes radiative cooling by heavy elements and feedback by supernovae and active galactic nuclei, predicts that by z = 0.25 nearly 60?per cent of the gas mass ends up at densities and temperatures characteristic of the WHIM and we find that half of this fraction is due to outflows. The standard H?i observables (distribution of H?i column densities NH?I, distribution of Doppler parameters bHI, bHINH?I correlation) and the BLA line number density predicted by our simulations are in remarkably good agreement with observations. BLAs arise in gas that is hotter, more highly ionized and more enriched than the gas giving rise to typical Lya forest absorbers. The majority of the BLAs arise in warm-hot [log?(T/?K) similar to 5] gas at low (log?? < 1.5) overdensities. On average, thermal broadening accounts for at least 60?per cent of the BLA linewidth, which in turn can be used as a rough indicator of the thermal state of the gas. Detectable BLAs account for only a small fraction of the true baryon content of the WHIM at low redshift. In order to detect the bulk of the mass in this gas phase, a sensitivity at least one order of magnitude better than achieved by current ultraviolet spectrographs is required. We argue that BLAs mostly trace gas that has been shock heated and enriched by outflows and that they therefore provide an important window on a poorly understood feedback process.
At z < 1 a large fraction of the baryons is thought to reside in diffuse gas that has been shock-heated to high temperatures (10 (5)-10 (6) K). Absorption by the 770.41, 780.32 A doublet of Ne viii in quasar spectra represents a unique tool to study this elusive warm-hot phase. We have developed an analytic model for the properties of Ne viii absorbers that allows for an inhomogeneous metal distribution. Our model agrees with the predictions of a simulation from the OverWhelmingly Large Simulations project indicating that the average line-of-sight metal-filling fraction within the absorbing gas is low (c(L) similar to 0.1). Most of the Ne viii in our model is produced in low-density, collisionally ionized gas (n(H) = 10(-6)-10(-4) cm(-3), T = 10 (5)-10 (6) K). Strong Ne viii absorbers (log(10)(N-NeVIII/cm(-2))14), like those recently detected by Hubble Space Telescope/Cosmic Origins Spectrograph, are found to arise in higher density gas (n(H) greater than or similar to 10(-4) cm(-3), T approximate to 5 x 10 (5) K). Ne viii cloudlets harbour only 1 per cent of the cosmic baryon budget. The baryon content of the surrounding gas (which has similar densities and temperatures as the Ne viii cloudlets) is a factor c(-1)L higher. We conclude that Ne viii absorbers are robust probes of shock-heated diffuse gas, but that spectra with signal-to-noise ratios S/N > 100 would be required to detect the bulk of the baryons in warm-hot gas.
Context. The High Energy Stereoscopic System Galactic plane survey (HGPS) is to date the most comprehensive census of Galactic gamma -ray sources at very high energies (VHE; 100 GeV <= E <= 100 TeV). As a consequence of the limited sensitivity of this survey, the 78 detected gamma -ray sources comprise only a small and biased subsample of the overall population. The larger part consists of currently unresolved sources, which contribute to large-scale diffuse emission to a still uncertain amount.Aims. We study the VHE gamma -ray source population in the Milky Way. For this purpose population-synthesis models are derived based on the distributions of source positions, extents, and luminosities.Methods. Several azimuth-symmetric and spiral-arm models are compared for spatial source distribution. The luminosity and radius function of the population are derived from the source properties of the HGPS data set and are corrected for the sensitivity bias of the HGPS. Based on these models, VHE source populations are simulated and the subsets of sources detectable according to the HGPS are compared with HGPS sources.Results. The power-law indices of luminosity and radius functions are determined to range between -1.6 and -1.9 for luminosity and -1.1 and -1.6 for radius. A two-arm spiral structure with central bar is discarded as spatial distribution of VHE sources, while azimuth-symmetric distributions and a distribution following a four-arm spiral structure without bar describe the HGPS data reasonably well. The total number of Galactic VHE sources is predicted to be in the range from 800 to 7000 with a total luminosity and flux of (1.6-6.3) x 10(36) ph s(-1) and (3-15) x 10(-10) ph cm(-2) s(-1), respectively.Conclusions. Depending on the model, the HGPS sample accounts for (68-87)% of the emission of the population in the scanned region. This suggests that unresolved sources represent a critical component of the diffuse emission measurable in the HGPS. With the foreseen jump in sensitivity of the Cherenkov Telescope Array, the number of detectable sources is predicted to increase by a factor between 5 and 9.