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 - Tepper-Garcia, Thor A1 - Richter, Philipp A1 - Schaye, Joop T1 - Absorption signatures of warm-hot gas at low redshift - ne viii JF - Monthly notices of the Royal Astronomical Society N2 - 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. KW - methods: analytical KW - methods: numerical KW - galaxies: formation KW - intergalactic medium KW - quasars: absorption lines KW - cosmology: theory Y1 - 2013 U6 - https://doi.org/10.1093/mnras/stt1712 SN - 0035-8711 SN - 1365-2966 VL - 436 IS - 3 SP - 2063 EP - 2081 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Tepper-Garcia, Thorsten A1 - Richter, Philipp A1 - Schaye, Joop A1 - Booth, C. M. A1 - Vecchia, Claudio Dalla A1 - Theuns, Tom A1 - Wiersma, Robert P. C. T1 - Absorption signatures of warm-hot gas at low redshift o vi JF - Monthly notices of the Royal Astronomical Society N2 - 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. KW - methods: numerical KW - galaxies: formation KW - intergalactic medium KW - quasars: absorption lines KW - cosmology: theory Y1 - 2011 U6 - https://doi.org/10.1111/j.1365-2966.2010.18123.x SN - 0035-8711 VL - 413 IS - 1 SP - 190 EP - 212 PB - Wiley-Blackwell CY - Hoboken ER - TY - JOUR A1 - Tepper-Garcia, Thorsten A1 - Richter, Philipp A1 - Schaye, Joop A1 - Booth, C. M. A1 - Dalla Vecchia, Claudio A1 - Theuns, Tom T1 - Absorption signatures of warm-hot gas at low redshift: broad H?i Lya absorbers JF - Monthly notices of the Royal Astronomical Society N2 - 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. KW - methods: numerical KW - galaxies: formation KW - intergalactic medium KW - quasars: absorption lines KW - cosmology: theory Y1 - 2012 U6 - https://doi.org/10.1111/j.1365-2966.2012.21545.x SN - 0035-8711 VL - 425 IS - 3 SP - 1640 EP - 1663 PB - Wiley-Blackwell CY - Hoboken ER - TY - JOUR A1 - Valori, Gherarod A1 - Demoulin, Pascal A1 - Pariat, E. A1 - Masson, S. T1 - Accuracy of magnetic energy computations JF - Astronomy and astrophysics : an international weekly journal N2 - 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. KW - magnetic fields KW - methods: numerical KW - Sun: surface magnetism KW - Sun: corona Y1 - 2013 U6 - https://doi.org/10.1051/0004-6361/201220982 SN - 0004-6361 VL - 553 IS - 2 PB - EDP Sciences CY - Les Ulis 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 - 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 - Thomas, T. A1 - Pfrommer, Christoph T1 - Cosmic-ray hydrodynamics BT - alfvén-wave regulated transport of cosmic rays JF - Monthly notices of the Royal Astronomical Society N2 - 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. KW - hydrodynamics KW - radiative transfer KW - methods: analytical KW - methods: numerical KW - cosmic rays Y1 - 2019 U6 - https://doi.org/10.1093/mnras/stz263 SN - 0035-8711 SN - 1365-2966 VL - 485 IS - 3 SP - 2977 EP - 3008 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Sachse, Manuel A1 - Kappel, David A1 - Tirsch, Daniela A1 - Otto, Katharina A. T1 - Discrete element modeling of aeolian-like morphologies on comet 67P/Churyumov-Gerasimenko JF - Astronomy and astrophysics : an international weekly journal N2 - Context. Even after the Rosetta mission, some of the mechanical parameters of comet 67P/Churyumov-Gerasimenko's surface material are still not well constrained. They are needed to improve our understanding of cometary activity or for planning sample return procedures. Aims. We discuss the physical process dominating the formation of aeolian-like surface features in the form of moats and wind taillike bedforms around obstacles and investigate the mechanical and geometrical parameters involved. Methods. By applying the discrete element method (DEM) in a low-gravity environment, we numerically simulated the dynamics of the surface layer particles and the particle stream involved in the formation of aeolian-like morphological features. The material is composed of polydisperse spherical particles that consist of a mixture of dust and water ice, with interparticle forces given by the Hertz contact model, cohesion, friction, and rolling friction. We determined a working set of parameters that enables simulations to be reasonably realistic and investigated morphological changes when modifying these parameters. Results. The aeolian-like surface features are reasonably well reproduced using model materials with a tensile strength on the order of 0.1-1 Pa. Stronger materials and obstacles with round shapes impede the formation of a moat and a wind tail. The integrated dust flux required for the formation of moats and wind tails is on the order of 100 kg m(-2), which, based on the timescale of morphological changes inferred from Rosetta images, translates to a near-surface particle density on the order of 10(-6)-10(-4) kg m(-3). Conclusions. DEM modeling of the aeolian-like surface features reveals complex formation mechanisms that involve both deposition of ejected material and surface erosion. More numerical work and additional in situ measurements or sample return missions are needed to better investigate mechanical parameters of cometary surface material and to understand the mechanics of cometary activity. KW - comets: general KW - comets: individual: 67P/Churyumov-Gerasimenko KW - methods: numerical Y1 - 2022 U6 - https://doi.org/10.1051/0004-6361/202141296 SN - 0004-6361 SN - 1432-0746 VL - 662 PB - EDP Sciences CY - Les Ulis ER - TY - JOUR A1 - Kappel, David A1 - Sachse, Manuel A1 - Haack, David A1 - Otto, Katharina A. T1 - Discrete element modeling of boulder and cliff morphologies on comet 67P/Churyumov-Gerasimenko JF - Astronomy and astrophysics : an international weekly journal N2 - Context: Even after the Rosetta mission, some of the mechanical parameters of comet 67P/Churyumov-Gerasimenko's surface material are not yet well constrained. These parameters are needed to improve our understanding of cometary activity or for planning sample return missions. Aims: We study some of the physical processes involved in the formation of selected surface features and investigate the mechanical and geometrical parameters involved. Methods: Applying the discrete element method (DEM) in a low-gravity environment, we numerically simulated the surface layer particle dynamics involved in the formation of selected morphological features. The material considered is a mixture of polydisperse ice and dust spheres with inter-particle forces given by the Hertz contact model, translational friction, rolling friction, cohesion from unsintered contacts, and optionally due to bonds from ice sintering. We determined a working set of parameters that enables the simulations to be reasonably realistic and investigated morphological changes due to modifications thereof. Results: The selected morphological features are reasonably well reproduced using model materials with a tensile strength on the order of 1-10 Pa. Increasing the diameters of the spherical particles decreases the material strength, and increasing the friction leads to a more brittle but somewhat stronger material. High friction is required to make the material sufficiently brittle to match observations, which points to the presence of very rough, even angular particles. Reasonable seismic activity does not suffice to trigger the collapses of cliffs without material heterogeneities or structural defects. Conclusions: DEM modeling can be a powerful tool to investigate mechanical parameters of cometary surface material. However, many uncertainties arise from our limited understanding of particle shapes, spatial configurations, and size distributions, all on multiple length scales. Further numerical work, in situ measurements, and sample return missions are needed to better understand the mechanics of cometary material and cometary activity. KW - comets: general KW - comets: individual: 67P KW - Churyumov-Gerasimenko KW - methods: numerical Y1 - 2020 U6 - https://doi.org/10.1051/0004-6361/201937152 SN - 0004-6361 SN - 1432-0746 VL - 641 PB - EDP Sciences CY - Les Ulis ER - TY - JOUR A1 - Bohdan, Artem A1 - Niemiec, Jacek A1 - Kobzar, Oleh A1 - Pohl, Martin T1 - Electron Pre-acceleration at Nonrelativistic High-Mach-number Perpendicular Shocks JF - The astrophysical journal : an international review of spectroscopy and astronomical physics N2 - We perform particle-in-cell simulations of perpendicular nonrelativistic collisionless shocks to study electron heating and pre-acceleration for parameters that permit the extrapolation to the conditions at young supernova remnants. Our high-resolution large-scale numerical experiments sample a representative portion of the shock surface and demonstrate that the efficiency of electron injection is strongly modulated with the phase of the shock reformation. For plasmas with low and moderate temperature (plasma beta beta p =5.10(-4) and 0.5 beta p =), we explore the nonlinear shock structure and electron pre-acceleration for various orientations of the large-scale magnetic field with respect to the simulation plane, while keeping it at 90 degrees to the shock normal. Ion reflection off of the shock leads to the formation of magnetic filaments in the shock ramp, resulting from Weibel-type instabilities, and electrostatic Buneman modes in the shock foot. In all of the cases under study, the latter provides first-stage electron energization through the shock-surfing acceleration mechanism. The subsequent energization strongly depends on the field orientation and proceeds through adiabatic or second-order Fermi acceleration processes for configurations with the out-of-plane and in-plane field components, respectively. For strictly out-of-plane field, the fraction of suprathermal electrons is much higher than for other configurations, because only in this case are the Buneman modes fully captured by the 2D simulation grid. Shocks in plasma with moderate bp provide more efficient pre-acceleration. The relevance of our results to the physics of fully 3D systems is discussed. KW - acceleration of particles KW - instabilities KW - ISM: supernova remnants KW - methods: numerical KW - plasmas KW - shock Y1 - 2017 U6 - https://doi.org/10.3847/1538-4357/aa872a SN - 0004-637X SN - 1538-4357 VL - 847 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Hani, Maan H. A1 - Sparre, Martin A1 - Ellison, Sara L. A1 - Torrey, Paul A1 - Vogelsberger, Mark T1 - Galaxy mergers moulding the circum-galactic medium BT - I. The impact of a major merger JF - Monthly notices of the Royal Astronomical Society N2 - Galaxies are surrounded by sizeable gas reservoirs which host a significant amount of metals: the circum-galactic medium (CGM). The CGM acts as a mediator between the galaxy and the extragalactic medium. However, our understanding of how galaxy mergers, a major evolutionary transformation, impact the CGM remains deficient. We present a theoretical study of the effect of galaxy mergers on the CGM. We use hydrodynamical cosmological zoom-in simulations of a major merger selected from the Illustris project such that the z = 0 descendant has a halo mass and stellar mass comparable to the Milky Way. To study the CGM we then re-simulated this system at a 40 times better mass resolution, and included detailed post-processing ionization modelling. Our work demonstrates the effect the merger has on the characteristic size of the CGM, its metallicity, and the predicted covering fraction of various commonly observed gas-phase species, such as H I, C IV, and O VI. We show that merger-induced outflows can increase the CGM metallicity by 0.2-0.3 dex within 0.5 Gyr post-merger. These effects last up to 6 Gyr post-merger. While the merger increases the total metal covering fractions by factors of 2-3, the covering fractions of commonly observed UV ions decrease due to the hard ionizing radiation from the active galactic nucleus, which we model explicitly. Our study of the single simulated major merger presented in this work demonstrates the significant impact that a galaxy interaction can have on the size, metallicity, and observed column densities of the CGM. KW - methods: numerical KW - galaxies: evolution KW - galaxies: haloes KW - galaxies: interactions Y1 - 2017 U6 - https://doi.org/10.1093/mnras/stx3252 SN - 0035-8711 SN - 1365-2966 VL - 475 IS - 1 SP - 1160 EP - 1176 PB - Oxford Univ. 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 - Seiler, Michael A1 - Seiß, Martin A1 - Hoffmann, Holger A1 - Spahn, Frank T1 - Hydrodynamic Simulations of Asymmetric Propeller Structures in Saturn's Rings JF - The astrophysical journal : an international review of spectroscopy and astronomical physics ; Supplement series N2 - The observation of the non-Keplerian behavior of propeller structures in Saturn's outer A ring raises the question: how does the propeller respond to the wandering of the central embedded moonlet? Here, we study numerically how the structural imprint of the propeller changes for a libration of the moonlet. It turns out that the libration induces an asymmetry in the propeller, which depends on the libration period and amplitude of the moonlet. Further, we study the dependence of the asymmetry on the libration period and amplitude for a moonlet with a 400 m Hill radius, which is located in the outer A ring. This allows us to apply our findings to the largest known propeller Blériot, which is expected to be of a similar size. For Blériot, we can conclude that, supposing the moonlet is librating with the largest observed period of 11.1 yr and an azimuthal amplitude of about 1845 km, a small asymmetry should be measurable but depends on the moonlet's libration phase at the observation time. The longitude residuals of other trans-Encke propellers (e.g., Earhart) show amplitudes similar to Blériot, which might allow us to observe larger asymmetries due to their smaller azimuthal extent, allowing us to scan the whole gap structure for asymmetries in one observation. Although the librational model of the moonlet is a simplification, our results are a first step toward the development of a consistent model for the description of the formation of asymmetric propellers caused by a freely moving moonlet. KW - Hydrodynamics KW - methods: data analysis KW - methods: numerical KW - planets and satellites: dynamical evolution and stability KW - planets and satellites: individual (Saturn) KW - planets and satellites: rings Y1 - 2019 U6 - https://doi.org/10.3847/1538-4365/ab26b0 SN - 0067-0049 SN - 1538-4365 VL - 243 IS - 2 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Sparre, Martin A1 - Pfrommer, Christoph A1 - Ehlert, Kristian T1 - Interaction of a cold cloud with a hot wind BT - the regimes of cloud growth and destruction and the impact of magnetic fields JF - Monthly notices of the Royal Astronomical Society N2 - Multiphase galaxy winds, the accretion of cold gas through galaxy haloes, and gas stripping from jellyfish galaxies are examples of interactions between cold and hot gaseous phases. There are two important regimes in such systems. A sufficiently small cold cloud is destroyed by the hot wind as a result of Kelvin-Helmholtz instabilities, which shatter the cloud into small pieces that eventually mix and dissolve in the hot wind. In contrast, stripped cold gas from a large cloud mixes with the hot wind to intermediate temperatures, and then becomes thermally unstable and cools, causing a net accretion of hot gas to the cold tail. Using the magneto-hydrodynamical code AREPO, we perform cloud crushing simulations and test analytical criteria for the transition between the growth and destruction regimes to clarify a current debate in the literature. We find that the hot-wind cooling time sets the transition radius and not the cooling time of the mixed phase. Magnetic fields modify the wind-cloud interaction. Draping of wind magnetic field enhances the field upstream of the cloud, and fluid instabilities are suppressed by a turbulently magnetized wind beyond what is seen for a wind with a uniform magnetic field. We furthermore predict jellyfish galaxies to have ordered magnetic fields aligned with their tails. We finally discuss how the results of idealized simulations can be used to provide input to subgrid models in cosmological (magneto-)hydrodynamical simulations, which cannot resolve the detailed small-scale structure of cold gas clouds in the circumgalactic medium. KW - methods: numerical KW - ISM: jets and outflows KW - galaxies: formation Y1 - 2020 U6 - https://doi.org/10.1093/mnras/staa3177 SN - 0035-8711 SN - 1365-2966 VL - 499 IS - 3 SP - 4261 EP - 4281 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Bohdan, Artem A1 - Niemiec, Jacek A1 - Pohl, Martin A1 - Matsumoto, Yosuke A1 - Amano, Takanobu A1 - Hoshino, Masahiro T1 - Kinetic Simulations of Nonrelativistic Perpendicular Shocks of Young Supernova Remnants BT - I. Electron Shock-surfing Acceleration JF - The astrophysical journal : an international review of spectroscopy and astronomical physics N2 - Electron injection at high Mach number nonrelativistic perpendicular shocks is studied here for parameters that are applicable to young SNR shocks. Using high-resolution large-scale two-dimensional fully kinetic particle-in-cell simulations and tracing individual particles, we in detail analyze the shock-surfing acceleration (SSA) of electrons at the leading edge of the shock foot. The central question is to what degree the process can be captured in 2D3V simulations. We find that the energy gain in SSA always arises from the electrostatic field of a Buneman wave. Electron energization is more efficient in the out-of-plane orientation of the large-scale magnetic field because both the phase speed and the amplitude of the waves are higher than for the in-plane scenario. Also, a larger number of electrons is trapped by the waves compared to the in-plane configuration. We conclude that significant modifications of the simulation parameters are needed to reach the same level of SSA efficiency as in simulations with out-of-plane magnetic field or 3D simulations. KW - acceleration of particles KW - instabilities KW - ISM: supernova remnants KW - methods: numerical KW - plasmas KW - shock waves Y1 - 2019 U6 - https://doi.org/10.3847/1538-4357/ab1b6d SN - 0004-637X SN - 1538-4357 VL - 878 IS - 1 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Mizuno, Yosuke A1 - Pohl, Martin A1 - Niemiec, Jacek A1 - Zhang, Bing A1 - Nishikawa, Ken-Ichi A1 - Hardee, Philip E. T1 - Magnetic field amplification and saturation in turbulence behind a relativistic shock JF - Monthly notices of the Royal Astronomical Society N2 - We have investigated via 2D relativistic magnetohydrodynamic simulations the long-term evolution of turbulence created by a relativistic shock propagating through an inhomogeneous medium. In the post-shock region, magnetic field is strongly amplified by turbulent motions triggered by pre-shock density inhomogeneities. Using a long-simulation box we have followed the magnetic field amplification until it is fully developed and saturated. The turbulent velocity is subrelativistic even for a strong shock. Magnetic field amplification is controlled by the turbulent motion and saturation occurs when the magnetic energy is comparable to the turbulent kinetic energy. Magnetic field amplification and saturation depend on the initial strength and direction of the magnetic field in the pre-shock medium, and on the shock strength. If the initial magnetic field is perpendicular to the shock normal, the magnetic field is first compressed at the shock and then can be amplified by turbulent motion in the post-shock region. Saturation occurs when the magnetic energy becomes comparable to the turbulent kinetic energy in the post-shock region. If the initial magnetic field in the pre-shock medium is strong, the post-shock region becomes turbulent but significant field amplification does not occur. If the magnetic energy after shock compression is larger than the turbulent kinetic energy in the post-shock region, significant field amplification does not occur. We discuss possible applications of our results to gamma-ray bursts and active galactic nuclei. KW - MHD KW - relativistic processes KW - shock waves KW - turbulence KW - methods: numerical KW - gamma-ray burst: general Y1 - 2014 U6 - https://doi.org/10.1093/mnras/stu196 SN - 0035-8711 SN - 1365-2966 VL - 439 IS - 4 SP - 3490 EP - 3503 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Mizuno, Yosuke A1 - Pohl, Martin A1 - Niemiec, Jacek A1 - Zhang, Bing A1 - Nishikawa, Ken-Ichi A1 - Hardee, Philip E. T1 - Magnetic-field amplification by turbulence in a relativistic shockpropagating through an inhomogeneous medium JF - The astrophysical journal : an international review of spectroscopy and astronomical physics N2 - We perform two-dimensional relativistic magnetohydrodynamic simulations of a mildly relativistic shock propagating through an inhomogeneous medium. We show that the postshock region becomes turbulent owing to preshock density inhomogeneity, and the magnetic field is strongly amplified due to the stretching and folding of field lines in the turbulent velocity field. The amplified magnetic field evolves into a filamentary structure in two-dimensional simulations. The magnetic energy spectrum is flatter than the Kolmogorov spectrum and indicates that a so-called small-scale dynamo is occurring in the postshock region. We also find that the amount of magnetic-field amplification depends on the direction of the mean preshock magnetic field, and the timescale of magnetic-field growth depends on the shock strength. KW - gamma-ray burst: general KW - magnetohydrodynamics (MHD) KW - methods: numerical KW - relativistic processes KW - shock waves KW - turbulence Y1 - 2011 U6 - https://doi.org/10.1088/0004-637X/726/2/62 SN - 0004-637X VL - 726 IS - 2 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Bustamante, Sebastian A1 - Sparre, Martin A1 - Springel, Volker A1 - Grand, Robert J. J. T1 - Merger-induced metallicity dilution in cosmological galaxy formation simulations JF - Monthly notices of the Royal Astronomical Society N2 - Observational studies have revealed that galaxy pairs tend to have lower gas-phase metallicity than isolated galaxies. This metallicity deficiency can be caused by inflows of low-metallicity gas due to the tidal forces and gravitational torques associated with galaxy mergers, diluting the metal content of the central region. In this work we demonstrate that such metallicity dilution occurs in state-of-the-art cosmological simulations of galaxy formation. We find that the dilution is typically 0.1 dex for major mergers, and is noticeable at projected separations smaller than 40 kpc. For minor mergers the metallicity dilution is still present, even though the amplitude is significantly smaller. Consistent with previous analysis of observed galaxies we find that mergers are outliers from the fundamental metallicity relation, with deviations being larger than expected for a Gaussian distribution of residuals. Our large sample of mergers within full cosmological simulations also makes it possible to estimate how the star formation rate enhancement and gas consumption timescale behave as a function of the merger mass ratio. We confirm that strong starbursts are likely to occur in major mergers, but they can also arise in minor mergers if more than two galaxies are participating in the interaction, a scenario that has largely been ignored in previous work based on idealised isolated merger simulations. KW - methods: numerical KW - galaxies: interactions KW - galaxies: star formation KW - galaxies: evolution Y1 - 2018 U6 - https://doi.org/10.1093/mnras/sty1692 SN - 0035-8711 SN - 1365-2966 VL - 479 IS - 3 SP - 3381 EP - 3392 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Steppa, Constantin A1 - Egberts, Kathrin T1 - Modelling the Galactic very-high-energy gamma-ray source population JF - Astronomy and astrophysics : an international weekly journal N2 - 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. KW - astroparticle physics KW - gamma rays: general KW - gamma rays: diffuse KW - background KW - methods: observational KW - methods: numerical Y1 - 2020 U6 - https://doi.org/10.1051/0004-6361/202038172 SN - 0004-6361 SN - 1432-0746 VL - 643 PB - EDP Sciences CY - Les Ulis ER -