Search

Cosmic rays (CRs) are ubiquitous in the interstellar medium (ISM) of nearby galaxies, but many of their properties are not well-constrained. Gamma-ray observations provide a powerful tool in this respect, allowing us to constrain both the interaction of CR protons with the ISM and their transport properties. To help better understand the link between observational signatures and CR physics, we use a series of magneto-hydrodynamical (MHD) AREPO simulations of isolated galaxies performed using spectrally-resolved CR transport in every computational cell, with subsequent gamma-ray emission calculated using the CRAYON+ (Cosmic RAY emissiON) code. In each of our simulated halos, modelling the energy-dependent spatial diffusion of CRs leads to a more extended distribution of high-energy (~100 GeV) gamma rays compared to that predicted by a 'grey' steady-state model, which is especially visible in the corresponding emission maps and radial profiles. Despite this, the total gamma-ray spectra can often be well approximated by the steady-state model, although recovering the same spectral index typically requires a minor variation of the energy dependence of the diffusion coefficient. Our simulations reproduce the observed spectral indices and gamma-ray spectra of nearby star-forming galaxies and also match recent observations of the far infrared--gamma-ray relation. We find, however, that the spectrally resolved model yields marginally smaller luminosities for lower star formation rates compared to grey simulations of CRs. Our work highlights the importance of modelling spectrally resolved CR transport for an accurate prediction of spatially resolved high-energy gamma-ray emission, as will be probed by the upcoming Cherenkov Telescope Array observatory., Comment: 18 pages, 11 figures, submitted to MNRAS, comments are welcome!

Cosmic rays (CRs) are dynamically important for the formation and evolution of galaxies by regulating star formation and by powering galactic outflows. However, to what extent CRs regulate galaxy formation depends on the coupling strength of CRs with the ambient plasma and the effective CR transport speed along the magnetic field. Moreover, both properties sensitively depend on the CR momentum, which is largely unexplored in three-dimensional hydrodynamical simulations. We perform magneto-hydrodynamical simulations of entire galaxies with masses ranging from $10^{10}$ to $10^{12}\,\mathrm{M}_\odot$ and compare dynamically coupled CRs in the grey approximation with a spectrally resolved model that includes CR momenta from $0.1\,\mathrm{GeV}~c^{-1}$ to $100\,\mathrm{TeV}~c^{-1}$. We find that hadronic cooling of CRs dominates over Alfv\'{e}n cooling, with the latter emulating CR losses as a result of streaming of CRs down their pressure gradient. While star formation rates and galaxy morphologies are only mildly affected by the spectral CR modelling, mass loading factors of galactic outflows can differ by up to a factor of four in dwarf galaxies. All simulated low-mass halos ($M=10^{10}$, $10^{11}$, and $3\times10^{11}\,\mathrm{M}_\odot$) drive strong outflows, where CR transport is temporally dominated by advection. In contrast, the Milky Way-mass galaxy with $M=10^{12}\,\mathrm{M}_\odot$ does not drive sustained outflows, so that CR transport is entirely dominated by diffusion. The effective energy weighted diffusion coefficients vary by two orders of magnitude from the canonical energy-weighted values of $\langle{D}\rangle_{e_\mathrm{cr}}\sim10^{28}\,\mathrm{cm^2\,s^{-1}}$ in the disc up to $3\times10^{29}\,\mathrm{cm^2\,s^{-1}}$ in the circumgalactic medium, where we observe substantial temperature and CR pressure differences between our grey and spectral CR models., Comment: 20 pages, submitted to MNRAS, comments welcome

The $��$-ray emission of star-forming (SF) galaxies is attributed to hadronic interactions of cosmic ray (CR) protons with the interstellar gas and contributions from CR electrons via bremsstrahlung and inverse Compton (IC) scattering. The relative importance of these processes in different galaxy types is still unclear. We model these processes in three-dimensional magneto-hydrodynamical (MHD) simulations of the formation of isolated galactic discs using the moving-mesh code AREPO, including dynamically coupled CR protons and adopting different CR transport models. We calculate steady-state CR spectra and also account for the emergence of secondary electrons and positrons. This allows us to produce detailed $��$-ray maps, luminosities and spectra of our simulated galaxies at different evolutionary stages. Our simulations with anisotropic CR diffusion and a low CR injection efficiency at supernovae (SNe, $��_{\mathrm{SN}}=0.05$) can successfully reproduce the observed far infrared (FIR)-$��$-ray relation. Starburst galaxies are close to the calorimetric limit, where CR protons lose most of their energy due to hadronic interactions and hence, their $��$-ray emission is dominated by neutral pion decay. However, in low SF galaxies, the increasing diffusive losses soften the CR proton spectra due to energy-dependent diffusion, and likewise steepen the pionic $��$-ray spectra. In turn, IC emission hardens the total spectra and can contribute up to $\sim40$ per cent of the total luminosity in low SF galaxies. Furthermore, in order to match the observed $��$-ray spectra of starburst galaxies, we require a weaker energy dependence of the CR diffusion coefficient, $D\propto{E}^{0.3}$, in comparison to Milky Way-like galaxies., 20 pages, 11 figures, accepted to MNRAS

Poster contribution for the conference "New Perspectives on Galactic Magnetism", {"references":["Caprioli, D and Spitkovsky, A , 2014 ApJ, 783, 91","Pfrommer, C. et al., 2017, MNRAS, 465, 4500","Riquelme, M. A. and Spitkovsky, A., 2011, ApJ, 733, 63","Sarazin, C., 1999, ApJ, 520, 529","Selig, M. et al., 2015, A&A, 581, A126","Springel, V., 2010, MNRAS, 401, 791","Winner, G. et al., arXiv:1903.01467"]}

We present a novel implementation of cosmic rays (CR) in the magneto-hydrodynamic code FLASH. CRs are described as separate fluids with different energies. CR advection, energy dependent anisotropic diffusion with respect to the magnetic field and adiabatic losses to follow the evolution of spectra are taken into account. We present a first study of the transport and immediate (~150 kyr) dynamical impact of CRs on the turbulent magnetised interstellar medium around supernova remnants on scales up to 80 pc. CR diffusion quickly leads to an efficient acceleration of low-density gas (mainly perpendicular to the magnetic field) with accelerations up to two orders of magnitude above the thermal values. Peaked (at 1 GeV) CR injection spectra have a stronger impact on the dynamics than power-law spectra. For realistic magnetic field configurations low energy CRs (with smaller diffusion coefficients) distribute anisotropically with large spatial variations of a factor of ten and more. Adiabatic losses can change the local spectra perceptibly but do not have an integral effect on the dynamics at the spatial and temporal scales considered here. We discuss the potential global impact of CRs and anisotropic transport on the dynamical structure of the ISM and also detail the limitations of the model., 22 pages, submitted to MNRAS, comments welcome

The time evolution of the probability density function (PDF) of the mass density is formulated and solved for systems in free-fall using a simple appoximate function for the collapse of a sphere. We demonstrate that a pressure-free collapse results in a power-law tail on the high-density side of the PDF. The slope quickly asymptotes to the functional form $\mathrm{P}_v(\rho)\propto\rho^{-1.54}$ for the (volume-weighted) PDF and $\mathrm{P}_m(\rho)\propto\rho^{-0.54}$ for the corresponding mass-weighted distribution. From the simple approximation of the PDF we derive analytic descriptions for mass accretion, finding that dynamically quiet systems with narrow density PDFs lead to retarded star formation and low star formation rates. Conversely, strong turbulent motions that broaden the PDF accelerate the collapse causing a bursting mode of star formation. Finally, we compare our theoretical work with observations. The measured star formation rates are consistent with our model during the early phases of the collapse. Comparison of observed column density PDFs with those derived from our model suggests that observed star-forming cores are roughly in free-fall., Comment: accepted for publication, 13 pages

This thesis investigates the impact of the initial conditions on present-day star formation. Using numerical simulations, we follow the gravitational collapse of dense molecular clouds under different initial turbulent motions and initial density distributions. Our analysis focuses on the morphology of the cloud, the time and location of the formation of stars, the energetics during the collapse, the formation of clusters including their internal structure, their accretion behaviour as well as their mass distribution. The morphology of the cloud and the total number of stars are strongly influenced by the initial type of turbulence and the initial density profile. The results range from almost unperturbed cores with a single star to strongly filamentary cores with hundreds of stars in disconnected clusters. The internal structure of protostellar clusters is systematically but not significantly influenced by the initial conditions. Concerning the accretion rates as well as the dynamical interactions of stars within the clusters, we observe a fairly uniform behaviour, not reflecting the large variations in the initial conditions. The simulations presented in this thesis were performed using the grid-based code FLASH, developed mainly at the University of Chicago.

We investigate the impact of different initial conditions for the initial density profile and the initial turbulence on the formation process of protostellar clusters. We study the collapse of dense molecula r cloud cores with three-dimensional adaptive mesh refinement simulations. We focus our discussion on the distribution of the gas among the protostellar objects in the turbulent dynamical cores. Despite the large variations in the initial configurations and the resulting gas and cluster morphology we find that all stellar clusters follow a very similar gas accretion behaviour. Once secondary protostars begin to form, the central region of a cluster is efficiently shielded from further accretion. Hence, objects located close to the centre are starved of material, as indicated by a strong decrease of the central acc retion rate. This Fragmentation Induced Starvation occurs not only in rotationally supported discs and filaments, but also in more spherically symmetric clusters with complex chaotic motions., 15 pages, MNRAS, in press, small changes (referee's suggestions and comments)

We investigate the formation of protostellar clusters during the collapse of dense molecular cloud cores with a focus on the evolution of potential and kinetic energy, the degree of substructure, and the early phase of mass segregation. Our study is based on a series of hydrodynamic simulations of dense cores, where we vary the initial density profile and the initial turbulent velocity. In the three-dimensional adaptive mesh refinement simulations, we follow the dynamical formation of filaments and protostars until a star formation efficiency of 20%. Despite the different initial configurations, the global ensemble of all protostars in a setup shows a similar energy evolution and forms sub-virial clusters with an energy ratio $E_\mathrm{kin}/|E_\mathrm{pot}|\sim0.2$. Concentrating on the innermost central region, the clusters show a roughly virialised energy balance. However, the region of virial balance only covers the innermost $\sim10-30%$ of all the protostars. In all simulations with multiple protostars, the total kinetic energy of the protostars is higher than the kinetic energy of the gas cloud, although the protostars only contain 20% of the total mass. The clusters vary significantly in size, mass, and number of protostars, and show different degrees of substructure and mass segregation. Flat density profiles and compressive turbulent modes produce more subclusters then centrally concentrated profiles and solenoidal turbulence. We find that dynamical relaxation and hence dynamical mass segregation is very efficient in all cases from the very beginning of the nascent cluster, i.e., during a phase when protostars are constantly forming and accreting., Comment: 19 pages, MNRAS accepted

We present a detailed parameter study of collapsing turbulent cloud cores, varying the initial density profile and the initial turbulent velocity field. We systematically investigate the influence of different initial conditions on the star formation process, mainly focusing on the fragmentation, the number of formed stars, and the resulting mass distributions. Our study compares four different density profiles (uniform, Bonnor-Ebert type, $\rho\propto r^{-1.5}$, and $\rho\propto r^{-2}$), combined with six different supersonic turbulent velocity fields (compressive, mixed, and solenoidal, initialised with two different random seeds each) in three-dimensional simulations using the adaptive-mesh refinement, hydrodynamics code FLASH. The simulations show that density profiles with flat cores produce hundreds of low-mass stars, either distributed throughout the entire cloud or found in subclusters, depending on the initial turbulence. Concentrated density profiles always lead to the formation of one high-mass star in the centre of the cloud and, if at all, low-mass stars surrounding the central one. In uniform and Bonnor-Ebert type density distributions, compressive initial turbulence leads to local collapse about 25% earlier than solenoidal turbulence. However, central collapse in the steep power-law profiles is too fast for the turbulence to have any significant influence. We conclude that (I) the initial density profile and turbulence mainly determine the cloud evolution and the formation of clusters, (II) the initial mass function (IMF) is not universal for all setups, and (III) that massive stars are much less likely to form in flat density distributions. The IMFs obtained in the uniform and Bonnor-Ebert type density profiles are more consistent with the observed IMF, but shifted to lower masses., Comment: 20 pages, MNRAS in press