Showing 1–50 of 74 results for author: Crocker, R M

TeV $γ$-ray emission near globular cluster Terzan 5 as a probe of cosmic ray transport

Authors: Mark R. Krumholz, Roland M. Crocker, Arash Bahramian, Pol Bordas

Abstract: Cosmic rays travelling through interstellar space have their propagation directions repeatedly scattered by fluctuating interstellar magnetic fields. The nature of this scattering is a major unsolved problem in astrophysics, one that has resisted solution largely due to a lack of direct observational constraints on the scattering rate. Here we show that very high-energy $γ$-ray emission from the g… ▽ More Cosmic rays travelling through interstellar space have their propagation directions repeatedly scattered by fluctuating interstellar magnetic fields. The nature of this scattering is a major unsolved problem in astrophysics, one that has resisted solution largely due to a lack of direct observational constraints on the scattering rate. Here we show that very high-energy $γ$-ray emission from the globular cluster Terzan 5, which has unexpectedly been found to be displaced from the cluster, presents a direct probe of this process. We show that this displacement is naturally explained by cosmic rays accelerated in the bow shock around the cluster propagating a finite distance before scattering processes re-orient enough of them towards Earth to produce a detectable $γ$-ray signal. The angular distance between the cluster and the signal places tight constraints on the scattering rate, which we show are consistent with a model whereby scattering is primarily due to excitation of magnetic waves by the cosmic rays themselves. The analysis method we develop here will make it possible to use sources with similarly displaced non-thermal X-ray and TeV $��$-ray signals as direct probes of cosmic ray scattering across a range of Galactic environments. △ Less

Submitted 26 June, 2024; originally announced June 2024.

Comments: Accepted to Nature Astronomy; 16 pages, 7 figures; this submission contains the main article, supplementary information, and extended data in a single document

CONGRuENTS (COsmic-ray, Neutrino, Gamma-ray and Radio Non-Thermal Spectra). II. Population-level correlations between galactic infrared, radio, and γ-ray emission

Authors: Matt A. Roth, Mark R. Krumholz, Roland M. Crocker, Todd A. Thompson

Abstract: Galaxies obey a number of empirical correlations between their radio, γ-ray, and infrared emission, but the physical origins of these correlations remain uncertain. Here we use the CONGRuENTS model for broadband non-thermal emission from star-forming galaxies, which self-consistently calculates energy-dependent transport and non-thermal emission from cosmic ray hadrons and leptons, to predict radi… ▽ More Galaxies obey a number of empirical correlations between their radio, γ-ray, and infrared emission, but the physical origins of these correlations remain uncertain. Here we use the CONGRuENTS model for broadband non-thermal emission from star-forming galaxies, which self-consistently calculates energy-dependent transport and non-thermal emission from cosmic ray hadrons and leptons, to predict radio and γ-ray emission for a synthetic galaxy population with properties drawn from a large deep-field survey. We show that our synthetic galaxies reproduce observed relations such as the FIR-radio correlation, the FIR-γ correlation, and the distribution of radio spectral indices, and we use the model to explain the physical origins of these relations. Our results show that the FIR-radio correlation arises because the amount of cosmic ray electron power ultimately radiated as synchrotron emission varies only weakly with galaxy star formation rate as a result of the constraints imposed on gas properties by hydrostatic balance and turbulent dynamo action; the same physics dictates the extent of proton calorimetry in different galaxies, and thus sets the FIR-γ-ray correlation. We further show that galactic radio spectral indices result primarily from competition between thermal free-free emission and energy-dependent loss of cosmic ray electrons to bremsstrahlung and escape into galactic halos, with shaping of the spectrum by inverse Compton, synchrotron, and ionisation processes typically playing a sub-dominant role. In addition to explaining existing observations, we use our analysis to predict a heretofore unseen correlation between the curvature of galaxies' radio spectra and their pion-driven γ-ray emission, a prediction that will be testable with upcoming facilities. △ Less

Submitted 15 June, 2024; v1 submitted 9 October, 2023; originally announced October 2023.

Comments: 17 pages, 14 figures

Journal ref: MNRAS, Volume 530, Issue 2, May 2024, Pages 1849-1865

Solar gamma ray probe of local cosmic ray electrons

Authors: Hong-Gang Yang, Yu Gao, Yin-Zhe Ma, Roland M. Crocker

Abstract: TeV-range cosmic ray electrons and positrons (CREs) have been directly measured in the search for new physics or unknown astrophysical sources. CREs can inverse-Compton scatter solar photons and boost their energies into gamma ray bands. Any potential CRE excess would enhance the resultant inverse Compton emission spectrum in the relevant energy range, offering a new window to verify the measured… ▽ More TeV-range cosmic ray electrons and positrons (CREs) have been directly measured in the search for new physics or unknown astrophysical sources. CREs can inverse-Compton scatter solar photons and boost their energies into gamma ray bands. Any potential CRE excess would enhance the resultant inverse Compton emission spectrum in the relevant energy range, offering a new window to verify the measured CRE spectrum. In this paper, we show that an excess in the TeV range of the CRE spectrum, such as the one indicated by the DAMPE experiment, can induce a characteristic solar gamma ray signal. Accounting for contamination from extragalactic gamma ray backgrounds (EGB), we forecast the DAMPE feature is testable ($\gtrsim 4 σ$) with a $\sim 10^{5}\,\mathrm{m}^2\,{\rm yr}$ exposure in the off-disk direction. This can be achieved by long-exposure observations of water Cherenkov telescopes, such as LHAASO (7.2 years) and HAWC (25.9 years). △ Less

Submitted 25 September, 2023; v1 submitted 9 September, 2023; originally announced September 2023.

Comments: 6 pages, 4 figures

Journal ref: Physical Review D (Letter) 108, L061304 (2023)

The spin axes of globular clusters and correlations with gamma-ray emission

Authors: Ciaran A. J. O'Hare, Alberto Krone-Martins, Celine Boehm, Roland M. Crocker

Abstract: A growing number of Milky Way globular clusters have been identified to possess a noticeable degree of solid-body rotation. For several clusters, the combination of stellar proper motions and radial velocities allows for 3-dimensional spin axes to be extracted. In this paper we consider the orientations of these spin axes, and ask whether they are correlated with any other properties of the cluste… ▽ More A growing number of Milky Way globular clusters have been identified to possess a noticeable degree of solid-body rotation. For several clusters, the combination of stellar proper motions and radial velocities allows for 3-dimensional spin axes to be extracted. In this paper we consider the orientations of these spin axes, and ask whether they are correlated with any other properties of the clusters -- either global properties to do with their orbits and origin, or internal properties related to the cluster composition. We discuss the possibility of alignments between the spin axes of globular clusters, chemodynamical groupings, and their orbital poles. We also point out a previously unidentified negative correlation between the measured gamma-ray emissivities and the inclination of the globular cluster spins with respect to the line of sight. Given that this correlation is not present in other wavelengths, we cannot conclusively attribute it solely to sampling bias. If the correlation holds up to scrutiny with more data, it may be indicative of sources of anisotropic gamma-ray emission in globular clusters. We discuss the plausibility of such an anisotropy arising from a population of dynamically formed millisecond pulsars with some degree of spin-orbit alignment. △ Less

Submitted 20 June, 2023; originally announced June 2023.

Comments: 22 pages, 14 figures

Prospective Dark Matter Annihilation Signals From the Sagittarius Dwarf Spheroidal

Authors: Thomas A. A. Venville, Alan R. Duffy, Roland M. Crocker, Oscar Macias, Thor Tepper-García

Abstract: The Sagittarius Dwarf Spheroidal galaxy (Sgr) is investigated as a target for DM annihilation searches utilising J-factor distributions calculated directly from a high-resolution hydrodynamic simulation of the infall and tidal disruption of Sgr around the Milky Way. In contrast to past studies, the simulation incorporates DM, stellar and gaseous components for both the Milky Way and the Sgr progen… ▽ More The Sagittarius Dwarf Spheroidal galaxy (Sgr) is investigated as a target for DM annihilation searches utilising J-factor distributions calculated directly from a high-resolution hydrodynamic simulation of the infall and tidal disruption of Sgr around the Milky Way. In contrast to past studies, the simulation incorporates DM, stellar and gaseous components for both the Milky Way and the Sgr progenitor galaxy. The simulated distributions account for significant tidal disruption affecting the DM density profile. Our estimate of the J-factor value for Sgr, $J_{\text{Sgr}}=1.48\times 10^{10}$ M$_\odot^2$ kpc$^{-5}$ ($6.46\times10^{16}\ \text{GeV}\ \text{cm}^{-5}$), is significantly lower than found in prior studies. This value, while formally a lower limit, is likely close to the true J-factor value for Sgr. It implies a DM cross-section incompatibly large in comparison with existing constraints would be required to attribute recently observed $γ$-ray emission from Sgr to DM annihilation. We also calculate a J-factor value using a NFW profile fitted to the simulated DM density distribution to facilitate comparison with past studies. This NFW J-factor value supports the conclusion that most past studies have overestimated the dark matter density of Sgr on small scales. This, together with the fact that the Sgr has recently been shown to emit $γ$-rays of astrophysical origin, complicate the use of Sgr in indirect DM detection searches. △ Less

Submitted 13 November, 2023; v1 submitted 23 March, 2023; originally announced March 2023.

Comments: 15 pages, 15 figures. Accepted for publication in MNRAS. Accepted manuscript uploaded

CONGRuENTS (COsmic-ray, Neutrino, Gamma-ray and Radio Non-Thermal Spectra). I. A predictive model for galactic non-thermal emission

Authors: Matt A. Roth, Mark R. Krumholz, Roland M. Crocker, Todd A. Thompson

Abstract: The total luminosity and spectral shape of the non-thermal emission produced by cosmic rays depends on their interstellar environment, a dependence that gives rise to correlations between galaxies' bulk properties -- star formation rate, stellar mass, and others -- and their non-thermal spectra. Understanding the physical mechanisms of cosmic ray transport, loss, and emission is key to understandi… ▽ More The total luminosity and spectral shape of the non-thermal emission produced by cosmic rays depends on their interstellar environment, a dependence that gives rise to correlations between galaxies' bulk properties -- star formation rate, stellar mass, and others -- and their non-thermal spectra. Understanding the physical mechanisms of cosmic ray transport, loss, and emission is key to understanding these correlations. Here, in the first paper of the series, we present a new method to compute the non-thermal spectra of star-forming galaxies, and describe an open-source software package -- COsmic-ray, Neutrino, Gamma-ray and Radio Non-Thermal Spectra (CONGRuENTS) -- that implements it. As a crucial innovation, our method requires as input only a galaxy's effective radius, star formation rate, stellar mass, and redshift, all quantities that are readily available for large samples of galaxies and do not require expensive, spatially resolved gas measurements. From these inputs we derive individual, galaxy-by-galaxy models for the background gas and radiation field through which cosmic rays propagate, from which we compute steady state cosmic ray spectra for hadronic and leptonic particles in both the galactic disc and halo by solving the full kinetic equation. We invoke modern models for cosmic ray transport and include all significant emission and loss mechanisms. In this paper we describe the model and validate it against non-thermal emission measured in nearby star-forming galaxies that span four orders of magnitude in star formation rate. △ Less

Submitted 16 May, 2023; v1 submitted 19 December, 2022; originally announced December 2022.

Comments: 23 pages, 14 figures, 1 table, accepted for publication in MNRAS

The cosmic ray ionisation and $γ$-ray budgets of star-forming galaxies

Authors: Mark R. Krumholz, Roland M. Crocker, Stella S. R. Offner

Abstract: Cosmic rays in star-forming galaxies are a dominant source of both diffuse $γ$-ray emission and ionisation in gas too deeply shielded for photons to penetrate. Though the cosmic rays responsible for $γ$-rays and ionisation are of different energies, they are produced by the same star formation-driven sources, and thus galaxies' star formation rates, $γ$-ray luminosities, and ionisation rates shoul… ▽ More Cosmic rays in star-forming galaxies are a dominant source of both diffuse $γ$-ray emission and ionisation in gas too deeply shielded for photons to penetrate. Though the cosmic rays responsible for $γ$-rays and ionisation are of different energies, they are produced by the same star formation-driven sources, and thus galaxies' star formation rates, $γ$-ray luminosities, and ionisation rates should all be linked. In this paper we use up-to-date cross-section data to determine this relationship, finding that cosmic rays in a galaxy of star formation rate $\dot{M}_*$ and gas depletion time $t_\mathrm{dep}$ produce a maximum primary ionisation rate $ζ\approx 1\times 10^{-16} (t_\mathrm{dep}/\mbox{Gyr})^{-1}$ s$^{-1}$ and a maximum $γ$-ray luminosity $L_γ\approx 4\times 10^{39} (\dot{M}_*/\mathrm{M}_\odot\mbox{ yr}^{-1})$ erg s$^{-1}$ in the 0.1 - 100 GeV band. These budgets imply either that the ionisation rates measured in Milky Way molecular clouds include a significant contribution from local sources that elevate them above the Galactic mean, or that CR-driven ionisation in the Milky Way is enhanced by sources not linked directly to star formation. Our results also imply that ionisation rates in starburst systems are only moderately enhanced compared to those in the Milky Way. Finally, we point out that measurements of $γ$-ray luminosities can be used to place constraints on galactic ionisation budgets in starburst galaxies that are nearly free of systematic uncertainties on the details of cosmic ray acceleration. △ Less

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

Comments: 19 pages, 8 figures, accepted for publication in MNRAS

Cosmic Ray Interstellar Propagation Tool using Itô Calculus (criptic): software for simultaneous calculation of cosmic ray transport and observational signatures

Authors: Mark R. Krumholz, Roland M. Crocker, Matt L. Sampson

Abstract: We present criptic, the Cosmic Ray Interstellar Propagation Tool using Itô Calculus, a new open-source software package to simulate the propagation of cosmic rays through the interstellar medium and to calculate the resulting observable non-thermal emission. Criptic solves the Fokker-Planck equation describing transport of cosmic rays on scales larger than that on which their pitch angles become a… ▽ More We present criptic, the Cosmic Ray Interstellar Propagation Tool using Itô Calculus, a new open-source software package to simulate the propagation of cosmic rays through the interstellar medium and to calculate the resulting observable non-thermal emission. Criptic solves the Fokker-Planck equation describing transport of cosmic rays on scales larger than that on which their pitch angles become approximately isotropic, and couples this to a rich and accurate treatment of the microphysical processes by which cosmic rays in the energy range $\sim$MeV to $\sim$PeV lose energy and produce emission. Criptic is deliberately agnostic as to both the cosmic ray transport model and the state of the background plasma through which cosmic rays travel. It can solve problems where cosmic rays stream, diffuse, or perform arbitrary combinations of both, and the coefficients describing these transport processes can be arbitrary functions of the background plasma state, the properties of the cosmic rays themselves, and local integrals of the cosmic ray field itself (e.g., the local cosmic ray pressure or pressure gradient). The code is parallelised using a hybrid OpenMP-MPI paradigm, allowing rapid calculations exploiting multiple cores and nodes on modern supercomputers. Here we describe the numerical methods used in the code, our treatment of the microphysical processes, and the set of code tests and validations we have performed. △ Less

Submitted 19 September, 2022; v1 submitted 27 July, 2022; originally announced July 2022.

Comments: 29 pages, 15 figures, accepted to MNRAS; very minor changes from first version. The software described is available from https://bitbucket.org/krumholz/criptic/src/master/, and a users' guide can be found at https://criptic.readthedocs.io/

Turbulent diffusion of streaming cosmic rays in compressible, partially ionised plasma

Authors: Matt L. Sampson, James R. Beattie, Mark R. Krumholz, Roland M. Crocker, Christoph Federrath, Amit Seta

Abstract: Cosmic rays (CRs) are a dynamically important component of the interstellar medium (ISM) of galaxies. The $\sim$GeV CRs that carry most CR energy and pressure are likely confined by self-generated turbulence, leading them to stream along magnetic field lines at the ion Alfvén speed. However, the consequences of self-confinement for CR propagation on galaxy scales remain highly uncertain. In this p… ▽ More Cosmic rays (CRs) are a dynamically important component of the interstellar medium (ISM) of galaxies. The $\sim$GeV CRs that carry most CR energy and pressure are likely confined by self-generated turbulence, leading them to stream along magnetic field lines at the ion Alfvén speed. However, the consequences of self-confinement for CR propagation on galaxy scales remain highly uncertain. In this paper, we use a large ensemble of magnetohydrodynamical turbulence simulations to quantify how the basic parameters describing ISM turbulence -- the sonic Mach number, $\mathcal{M}$ (plasma compressibility), Alfvén Mach number, $\mathcal{M}_{A0}$ (strength of the large-scale field with respect to the turbulence), and ionisation fraction by mass, $χ$ -- affect the transport of streaming CRs. We show that the large-scale transport of CRs whose small-scale motion consists of streaming along field lines is well described as a combination of streaming along the mean field and superdiffusion both along (parallel to) and across (perpendicular to) it; $\mathcal{M}_{A0}$ drives the level of anisotropy between parallel and perpendicular diffusion and $χ$ modulates the magnitude of the diffusion coefficients, while in our choice of units, $\mathcal{M}$ is unimportant except in the sub-Alfvénic ($\mathcal{M}_{A0} \lesssim 0.5$) regime. Our finding that superdiffusion is ubiquitous potentially explains the apparent discrepancy between CR diffusion coefficients inferred from measurements close to individual sources compared to those measured on larger, Galactic scales. Finally, we present empirical fits for the diffusion coefficients as a function of plasma parameters that may be used as sub-grid recipes for global interstellar medium, galaxy or cosmological simulations. △ Less

Submitted 3 November, 2022; v1 submitted 17 May, 2022; originally announced May 2022.

Comments: Accepted in MNRAS

Gamma-ray emission from the Sagittarius Dwarf Spheroidal galaxy due to millisecond pulsars

Authors: Roland M. Crocker, Oscar Macias, Dougal Mackey, Mark R. Krumholz, Shin'ichiro Ando, Shunsaku Horiuchi, Matthew G. Baring, Chris Gordon, Thomas Venville, Alan R. Duffy, Rui-Zhi Yang, Felix Aharonian, J. A. Hinton, Deheng Song, Ashley J. Ruiter, Miroslav D. Filipović

Abstract: The Fermi Bubbles are giant, gamma-ray emitting lobes emanating from the nucleus of the Milky Way discovered in ~1-100 GeV data collected by the Large Area Telescope on board the Fermi Gamma-Ray Space Telescope. Previous work has revealed substructure within the Fermi Bubbles that has been interpreted as a signature of collimated outflows from the Galaxy's super-massive black hole. Here we show vi… ▽ More The Fermi Bubbles are giant, gamma-ray emitting lobes emanating from the nucleus of the Milky Way discovered in ~1-100 GeV data collected by the Large Area Telescope on board the Fermi Gamma-Ray Space Telescope. Previous work has revealed substructure within the Fermi Bubbles that has been interpreted as a signature of collimated outflows from the Galaxy's super-massive black hole. Here we show via a spatial template analysis that much of the gamma-ray emission associated to the brightest region of substructure -- the so-called cocoon -- is likely due to the Sagittarius dwarf spheroidal (Sgr dSph) galaxy. This large Milky Way satellite is viewed through the Fermi Bubbles from the position of the Solar System. As a tidally and ram-pressure stripped remnant, the Sgr dSph has no on-going star formation, but we nevertheless demonstrate that the dwarf's millisecond pulsar (MSP) population can plausibly supply the gamma-ray signal that our analysis associates to its stellar template. The measured spectrum is naturally explained by inverse Compton scattering of cosmic microwave background photons by high-energy electron-positron pairs injected by MSPs belonging to the Sgr dSph, combined with these objects' magnetospheric emission. This finding plausibly suggests that MSPs produce significant gamma-ray emission amongst old stellar populations, potentially confounding indirect dark matter searches in regions such as the Galactic Centre, the Andromeda galaxy, and other massive Milky Way dwarf spheroidals. △ Less

Submitted 8 September, 2022; v1 submitted 25 April, 2022; originally announced April 2022.

Comments: Updated to match version accepted for publication in Nature Astronomy (2022). 14 pages main text, 3 main figures, 7 extended data figures. For published version of the paper, see https://rdcu.be/cUZ6X

Journal ref: Nature Astronomy (2022)

The Andromeda Gamma-Ray Excess: Background Systematics of the Millisecond Pulsars and Dark Matter Interpretations

Authors: Fabian Zimmer, Oscar Macias, Shin'ichiro Ando, Roland M. Crocker, Shunsaku Horiuchi

Abstract: Since the discovery of an excess in gamma rays in the direction of M31, its cause has been unclear. Published interpretations focus on a dark matter or stellar related origin. Studies of a similar excess in the Milky Way center motivate a correlation of the spatial morphology of the signal with the distribution of stellar mass in M31. However, a robust determination of the best theory for the obse… ▽ More Since the discovery of an excess in gamma rays in the direction of M31, its cause has been unclear. Published interpretations focus on a dark matter or stellar related origin. Studies of a similar excess in the Milky Way center motivate a correlation of the spatial morphology of the signal with the distribution of stellar mass in M31. However, a robust determination of the best theory for the observed excess emission is very challenging due to large uncertainties in the astrophysical gamma-ray foreground model. Here we perform a spectro-morphological analysis of the M31 gamma-ray excess using state-of-the-art templates for the distribution of stellar mass in M31 and novel astrophysical foreground models for its sky region. We construct maps for the old stellar populations of M31 based on observational data from the PAndAS survey and carefully remove the foreground stars. We also produce improved astrophysical foreground models by using novel image inpainting techniques based on machine learning methods. We find that our stellar maps, taken as a proxy for the location of a putative population of millisecond pulsars in the bulge of M31, reach a statistical significance of $5.4σ$, making them as strongly favoured as the simple phenomenological models usually considered in the literature, e.g., a disk-like template with uniform brightness. Our detection of the stellar templates is robust to generous variations of the astrophysical foreground model. Once the stellar templates are included in the astrophysical model, we show that the dark matter annihilation interpretation of the signal is unwarranted. Using the results of a binary population synthesis model we demonstrate that a population of about one million unresolved MSPs could naturally explain the observed gamma-ray luminosity per stellar mass, energy spectrum, and stellar bulge-to-disk flux ratio. △ Less

Submitted 1 April, 2022; originally announced April 2022.

Comments: 15 pages, 11 figures, comments are welcome

Ion Alfvén velocity fluctuations and implications for the diffusion of streaming cosmic rays

Authors: James R. Beattie, Mark R. Krumholz, Christoph Federrath, Matt Sampson, Roland M. Crocker

Abstract: The interstellar medium (ISM) of star-forming galaxies is magnetized and turbulent. Cosmic rays (CRs) propagate through it, and those with energies from $\sim\,\rm{GeV} - \rm{TeV}$ are likely subject to the streaming instability, whereby the wave damping processes balances excitation of resonant ionic Alfvén waves by the CRs, reaching an equilibrium in which the propagation speed of the CRs is ver… ▽ More The interstellar medium (ISM) of star-forming galaxies is magnetized and turbulent. Cosmic rays (CRs) propagate through it, and those with energies from $\sim\,\rm{GeV} - \rm{TeV}$ are likely subject to the streaming instability, whereby the wave damping processes balances excitation of resonant ionic Alfvén waves by the CRs, reaching an equilibrium in which the propagation speed of the CRs is very close to the local ion Alfvén velocity. The transport of streaming CRs is therefore sensitive to ionic Alfvén velocity fluctuations. In this paper we systematically study these fluctuations using a large ensemble of compressible MHD turbulence simulations. We show that for sub-Alfvénic turbulence, as applies for a strongly magnetized ISM, the ionic Alfvén velocity probability density function (PDF) is determined solely by the density fluctuations from shocked gas forming parallel to the magnetic field, and we develop analytical models for the ionic Alfvén velocity PDF up to second moments. For super-Alfvénic turbulence, magnetic and density fluctuations are correlated in complex ways, and these correlations as well as contributions from the magnetic fluctuations sets the ionic Alfvén velocity PDF. We discuss the implications of these findings for underlying "macroscopic" diffusion mechanisms in CRs undergoing the streaming instability, including modeling the macroscopic diffusion coefficient for the parallel transport in sub-Alfvénic plasmas. We also describe how, for highly-magnetized turbulent gas, the gas density PDF, and hence column density PDF, can be used to access information about ionic Alfvén velocity structure from observations of the magnetized ISM. △ Less

Submitted 7 August, 2022; v1 submitted 25 March, 2022; originally announced March 2022.

Comments: Invited submission accepted for publication in Frontiers in Astronomy and Space Sciences: Multi-scale Magnetic Field Measurements in the Multi-Phase Interstellar Medium; 46 pages; 12 Figures, 2 tables

MeerKAT uncovers the physics of an Odd Radio Circle

Authors: Ray P. Norris, J. D. Collier, Roland M. Crocker, Ian Heywood, Peter Macgregor, L. Rudnick, Stas Shabala, Heinz Andernach, Elisabete da Cunha, Jayanne English, Miroslav Filipovic, Baaerbel S. Koribalski, Kieran Luken, Aaron Robotham, Srikrishna Sekhar, Jessica E. Thorne, Tessa Vernstrom

Abstract: Odd Radio Circles (ORCs) are recently-discovered faint diffuse circles of radio emission, of unknown cause, surrounding galaxies at moderate redshift ($z ~ 0.2-0.6). Here we present detailed new MeerKAT radio images at 1284 MHz of the first ORC, originally discovered with the Australian Square Kilometre Array Pathfinder, with higher resolution (6 arcsec) and sensitivity (~ 2.4 uJy/bm). In additi… ▽ More Odd Radio Circles (ORCs) are recently-discovered faint diffuse circles of radio emission, of unknown cause, surrounding galaxies at moderate redshift ($z ~ 0.2-0.6). Here we present detailed new MeerKAT radio images at 1284 MHz of the first ORC, originally discovered with the Australian Square Kilometre Array Pathfinder, with higher resolution (6 arcsec) and sensitivity (~ 2.4 uJy/bm). In addition to the new images, which reveal a complex internal structure consisting of multiple arcs, we also present polarisation and spectral index maps. Based on these new data, we consider potential mechanisms that may generate the ORCs. △ Less

Submitted 20 March, 2022; originally announced March 2022.

Comments: Accepted by MNRAS

Energy balance and Alfvén Mach numbers in compressible magnetohydrodynamic turbulence with a large-scale magnetic field

Authors: James R. Beattie, Mark R. Krumholz, Raphael Skalidis, Christoph Federrath, Amit Seta, Roland M. Crocker, Philip Mocz, Neco Kriel

Abstract: Energy equipartition is a powerful theoretical tool for understanding astrophysical plasmas. It is invoked, for example, to measure magnetic fields in the interstellar medium (ISM), as evidence for small-scale turbulent dynamo action, and, in general, to estimate the energy budget of star-forming molecular clouds. In this study we motivate and explore the role of the volume-averaged root-mean-squa… ▽ More Energy equipartition is a powerful theoretical tool for understanding astrophysical plasmas. It is invoked, for example, to measure magnetic fields in the interstellar medium (ISM), as evidence for small-scale turbulent dynamo action, and, in general, to estimate the energy budget of star-forming molecular clouds. In this study we motivate and explore the role of the volume-averaged root-mean-squared (rms) magnetic coupling term between the turbulent, $δ\mathbf{B}$ and large-scale, $\mathbf{B}_0$ fields, $\left< (δ\mathbf{B}\cdot\mathbf{B}_0)^{2} \right>^{1/2}_{\mathcal{V}}$. By considering the second moments of the energy balance equations we show that the rms coupling term is in energy equipartition with the volume-averaged turbulent kinetic energy for turbulence with a sub-Alfvénic large-scale field. Under the assumption of exact energy equipartition between these terms, we derive relations for the magnetic and coupling term fluctuations, which provide excellent, parameter-free agreement with time-averaged data from 280 numerical simulations of compressible MHD turbulence. Furthermore, we explore the relation between the turbulent, mean-field and total Alfvén Mach numbers, and demonstrate that sub-Alfvénic turbulence can only be developed through a strong, large-scale magnetic field, which supports an extremely super-Alfvénic turbulent magnetic field. This means that the magnetic field fluctuations are significantly subdominant to the velocity fluctuations in the sub-Alfvénic large-scale field regime. Throughout our study, we broadly discuss the implications for observations of magnetic fields and understanding the dynamics in the magnetised ISM. △ Less

Submitted 24 July, 2022; v1 submitted 25 February, 2022; originally announced February 2022.

Comments: Accepted in MNRAS. 20 pages. 16 figures. 1 table

Mysterious Odd Radio Circle near the Large Magellanic Cloud -- An Intergalactic Supernova Remnant?

Authors: Miroslav D. Filipović, J. L. Payne, R. Z. E. Alsaberi, R. P. Norris, P. J. Macgregor, L. Rudnick, B. S. Koribalski, D. Leahy, L. Ducci, R. Kothes, H. Andernach, L. Barnes, I. S. Bojičić, L. M. Bozzetto, R. Brose, J. D. Collier, E. J. Crawford, R. M. Crocker, S. Dai, T. J. Galvin, F. Haberl, U. Heber, T. Hill, A. M. Hopkins, N. Hurley-Walker , et al. (26 additional authors not shown)

Abstract: We report the discovery of J0624-6948, a low-surface brightness radio ring, lying between the Galactic Plane and the Large Magellanic Cloud (LMC). It was first detected at 888 MHz with the Australian Square Kilometre Array Pathfinder (ASKAP), and with a diameter of ~196 arcsec. This source has phenomenological similarities to Odd Radio Circles (ORCs). Significant differences to the known ORCs - a… ▽ More We report the discovery of J0624-6948, a low-surface brightness radio ring, lying between the Galactic Plane and the Large Magellanic Cloud (LMC). It was first detected at 888 MHz with the Australian Square Kilometre Array Pathfinder (ASKAP), and with a diameter of ~196 arcsec. This source has phenomenological similarities to Odd Radio Circles (ORCs). Significant differences to the known ORCs - a flatter radio spectral index, the lack of a prominent central galaxy as a possible host, and larger apparent size - suggest that J0624-6948 may be a different type of object. We argue that the most plausible explanation for J0624-6948 is an intergalactic supernova remnant due to a star that resided in the LMC outskirts that had undergone a single-degenerate type Ia supernova, and we are seeing its remnant expand into a rarefied, intergalactic environment. We also examine if a massive star or a white dwarf binary ejected from either galaxy could be the supernova progenitor. Finally, we consider several other hypotheses for the nature of the object, including the jets of an active galactic nucleus (AGN) or the remnant of a nearby stellar super-flare. △ Less

Submitted 24 January, 2022; originally announced January 2022.

Comments: 20 pages accepted to MNRAS

The diffuse $γ$-ray background is dominated by star-forming galaxies

Authors: Matt A. Roth, Mark R. Krumholz, Roland M. Crocker, Silvia Celli

Abstract: The Fermi Gamma-ray Space Telescope has revealed a diffuse $γ$-ray background at energies from 0.1 GeV to 1 TeV, which can be separated into Galactic emission and an isotropic, extragalactic component. Previous efforts to understand the latter have been hampered by the lack of physical models capable of predicting the $γ$-ray emission produced by the many candidate sources, primarily active galact… ▽ More The Fermi Gamma-ray Space Telescope has revealed a diffuse $γ$-ray background at energies from 0.1 GeV to 1 TeV, which can be separated into Galactic emission and an isotropic, extragalactic component. Previous efforts to understand the latter have been hampered by the lack of physical models capable of predicting the $γ$-ray emission produced by the many candidate sources, primarily active galactic nuclei and star-forming galaxies, leaving their contributions poorly constrained. Here we present a calculation of the contribution of star-forming galaxies to the $γ$-ray background that does not rely on empirical scalings, and is instead based on a physical model for the $γ$-ray emission produced when cosmic rays accelerated in supernova remnants interact with the interstellar medium. After validating the model against local observations, we apply it to the observed cosmological star-forming galaxy population and recover an excellent match to both the total intensity and the spectral slope of the $γ$-ray background, demonstrating that star-forming galaxies alone can explain the full diffuse, isotropic $γ$-ray background. △ Less

Submitted 15 September, 2021; originally announced September 2021.

Comments: 18 pages, 10 figures. This work has been published in Nature. The version deposited here is the author's pre-print and may not reflect post-acceptance corrections or formatting related changes. The published version (Version of Record) of this manuscript is available at https://www.nature.com/articles/s41586-021-03802-x

Journal ref: Nature 597, 341-344 (2021)

First extragalactic measurement of the turbulence driving parameter: ALMA observations of the star-forming region N159E in the Large Magellanic Cloud

Authors: Piyush Sharda, Shyam H. Menon, Christoph Federrath, Mark R. Krumholz, James R. Beattie, Katherine E. Jameson, Kazuki Tokuda, Blakesley Burkhart, Roland M. Crocker, Charles J. Law, Amit Seta, Terrance J. Gaetz, Nickolas M. Pingel, Ivo R. Seitenzahl, Hidetoshi Sano, Yasuo Fukui

Abstract: Studying the driving modes of turbulence is important for characterizing the impact of turbulence in various astrophysical environments. The driving mode of turbulence is parameterized by $b$, which relates the width of the gas density PDF to the turbulent Mach number; $b\approx 1/3$, $1$, and $0.4$ correspond to driving that is solenoidal, compressive, and a natural mixture of the two, respective… ▽ More Studying the driving modes of turbulence is important for characterizing the impact of turbulence in various astrophysical environments. The driving mode of turbulence is parameterized by $b$, which relates the width of the gas density PDF to the turbulent Mach number; $b\approx 1/3$, $1$, and $0.4$ correspond to driving that is solenoidal, compressive, and a natural mixture of the two, respectively. In this work, we use high-resolution (sub-pc) ALMA $^{12}$CO ($J$ = $2-1$), $^{13}$CO ($J$ = $2-1$), and C$^{18}$O ($J$ = $2-1$) observations of filamentary molecular clouds in the star-forming region N159E (the Papillon Nebula) in the Large Magellanic Cloud (LMC) to provide the first measurement of turbulence driving parameter in an extragalactic region. We use a non-local thermodynamic equilibrium (NLTE) analysis of the CO isotopologues to construct a gas density PDF, which we find to be largely log-normal in shape with some intermittent features indicating deviations from lognormality. We find that the width of the log-normal part of the density PDF is comparable to the supersonic turbulent Mach number, resulting in $b \approx 0.9$. This implies that the driving mode of turbulence in N159E is primarily compressive. We speculate that the compressive turbulence could have been powered by gravo-turbulent fragmentation of the molecular gas, or due to compression powered by H I flows that led to the development of the molecular filaments observed by ALMA in the region. Our analysis can be easily applied to study the nature of turbulence driving in resolved star-forming regions in the local as well as the high-redshift Universe. △ Less

Submitted 19 October, 2021; v1 submitted 8 September, 2021; originally announced September 2021.

Comments: 14 pages, 7 figures. MNRAS in press

Measuring the smearing of the Galactic 511 keV signal: positron propagation or supernova kicks?

Authors: Thomas Siegert, Roland M. Crocker, Oscar Macias, Fiona H. Panther, Francesca Calore, Deheng Song, Shunsaku Horiuchi

Abstract: We use 15 years of $γ$-ray data from INTEGRAL/SPI in a refined investigation of the morphology of the Galactic bulge positron annihilation signal. Our spatial analysis confirms that the signal traces the old stellar population in the bulge and reveals for the first time that it traces the boxy bulge and nuclear stellar bulge. Using a 3D smoothing kernel, we find that the signal is smeared out over… ▽ More We use 15 years of $γ$-ray data from INTEGRAL/SPI in a refined investigation of the morphology of the Galactic bulge positron annihilation signal. Our spatial analysis confirms that the signal traces the old stellar population in the bulge and reveals for the first time that it traces the boxy bulge and nuclear stellar bulge. Using a 3D smoothing kernel, we find that the signal is smeared out over a characteristic length scale of $150 \pm 50\,$pc, suggesting either annihilation in situ at astrophysical sources kicked at formation or positron propagation away from sources. The former is disfavoured by its requiring kick velocities different between the Galactic nucleus ($\gtrsim 50\,\mathrm{km\,s^{-1}}$) and wider bulge ($\lesssim 15\,\mathrm{km\,s^{-1}}$) source. Positron propagation prior to annihilation can explain the overall phenomenology of the 511 keV signal for positrons injection energies $\lesssim 1.4\,$MeV, suggesting a nucleosynthesis origin. △ Less

Submitted 8 September, 2021; originally announced September 2021.

Comments: 5 pages, 4 figures, 2 tables; submitted with MNRAS Letters

Millisecond Pulsars from Accretion Induced Collapse as the Origin of the Galactic Centre Gamma-ray Excess Signal

Authors: Anuj Gautam, Roland M. Crocker, Lilia Ferrario, Ashley J. Ruiter, Harrison Ploeg, Chris Gordon, Oscar Macias

Abstract: Gamma-ray data from the Fermi-Large Area Telescope reveal an unexplained, apparently diffuse, signal from the Galactic bulge that peaks near 2 GeV with an approximately spherical intensity profile $\propto r^{-2.4}$ that extends to angular radial scales of at least 10 degrees, possibly to 20 degrees. The origin of this "Galactic Centre Excess" (GCE) has been debated with proposed sources prominent… ▽ More Gamma-ray data from the Fermi-Large Area Telescope reveal an unexplained, apparently diffuse, signal from the Galactic bulge that peaks near 2 GeV with an approximately spherical intensity profile $\propto r^{-2.4}$ that extends to angular radial scales of at least 10 degrees, possibly to 20 degrees. The origin of this "Galactic Centre Excess" (GCE) has been debated with proposed sources prominently including self-annihilating dark matter and a hitherto undetected population of millisecond pulsars (MSPs). However, the conventional channel for the generation of MSPs has been found to predict too many low mass X-ray binary (LMXB) systems and, because of the expected large natal kicks, may not accommodate the close spatial correspondence between the GCE signal and stars in the bulge. Here we report a binary population synthesis forward model that demonstrates that an MSP population arising from the accretion induced collapse (AIC) of O-Ne white dwarfs in Galactic bulge binaries can naturally reproduce the morphology, spectral shape, and intensity of the GCE signal while also obeying LMXB constraints. Synchrotron emission from MSP-launched cosmic ray electrons and positrons may simultaneously explain the mysterious, microwave "haze" from the inner Galaxy △ Less

Submitted 2 May, 2022; v1 submitted 1 June, 2021; originally announced June 2021.

Comments: 4 pages main text, 3 main figures, 6 pages Methods, 33 pages supplementary material. Updated to match version accepted for publication in Nature Astronomy (2022). For published version of the paper, see https://rdcu.be/cMjKY

Journal ref: Nature Astronomy (2022)

Prospects of direct detection of $^{48}$V gamma-rays from thermonuclear supernovae

Authors: Fiona H. Panther, Ivo R. Seitenzahl, Ashley J. Ruiter, Thomas Siegert, Stuart Sim, Roland M. Crocker

Abstract: Detection of gamma-rays emitted by radioactive isotopes synthesized in stellar explosions can give important insights into the processes that power transients such as supernovae, as well as providing a detailed census of the abundance of different isotope species relevant to the chemical evolution of the Universe. Observations of nearby supernovae have yielded observational proof that $^{57}$Co po… ▽ More Detection of gamma-rays emitted by radioactive isotopes synthesized in stellar explosions can give important insights into the processes that power transients such as supernovae, as well as providing a detailed census of the abundance of different isotope species relevant to the chemical evolution of the Universe. Observations of nearby supernovae have yielded observational proof that $^{57}$Co powered the late-time evolution of SN1987A's lightcurve, and conclusive evidence that $^{56}$Ni and its daughter nuclei power the light curves of Type Ia supernovae. In this paper we describe the prospects for detecting nuclear decay lines associated with the decay of $^{48}$V, the daughter nucleus of $^{48}$Cr, which is expected to be synthesised in large quantities - $M_{\mathrm{Cr}}\sim1.9\times10^{-2}\,\mathrm{M_\odot}$ - in transients initiated by explosive helium burning ($α$-capture) of a thick helium shell. We calculate emergent gamma-ray line fluxes for a simulated explosion model of a thermonuclear explosion of carbon-oxygen white dwarf core of mass $0.45\,M_{\odot}$ surrounded by a thick helium layer of mass $0.21\,M_{\odot}$. We present observational limits on the presence of $^{48}$V in nearby SNe Ia 2014J using the \textit{INTEGRAL} space telescope, excluding a $^{48}$Cr production on the surface of more than $0.1\,\mathrm{M_{\odot}}$. We find that the future gamma-ray mission AMEGO will have an approximately 5 per cent chance of observing $^{48}$V gamma-rays from such events during the currently-planned operational lifetime, based on our birthrate predictions of faint thermonuclear transients. We describe the conditions for a $3σ$ detection by the gamma-ray telescopes \textit{INTEGRAL}/SPI, COSI and AMEGO. △ Less

Submitted 8 September, 2021; v1 submitted 31 March, 2021; originally announced March 2021.

Comments: 9 pages, 3 figures, submitted to MNRAS, minor revisions Sept 2021

Cherenkov Telescope Array sensitivity to the putative millisecond pulsar population responsible for the Galactic center excess

Authors: Oscar Macias, Harm van Leijen, Deheng Song, Shin'ichiro Ando, Shunsaku Horiuchi, Roland M. Crocker

Abstract: The leading explanation of the $\textit{Fermi}$ Galactic center $γ$-ray excess is the extended emission from a unresolved population of millisecond pulsars (MSPs) in the Galactic bulge. Such a population would, along with the prompt $γ$ rays, also inject large quantities of electrons/positrons ($e^\pm$) into the interstellar medium. These $e^\pm$ could potentially inverse-Compton (IC) scatter ambi… ▽ More The leading explanation of the $\textit{Fermi}$ Galactic center $γ$-ray excess is the extended emission from a unresolved population of millisecond pulsars (MSPs) in the Galactic bulge. Such a population would, along with the prompt $γ$ rays, also inject large quantities of electrons/positrons ($e^\pm$) into the interstellar medium. These $e^\pm$ could potentially inverse-Compton (IC) scatter ambient photons into $γ$ rays that fall within the sensitivity range of the upcoming Cherenkov Telescope Array (CTA). In this article, we examine the detection potential of CTA to this signature by making a realistic estimation of the systematic uncertainties on the Galactic diffuse emission model at TeV-scale $γ$-ray energies. We forecast that, in the event that $e^\pm$ injection spectra are harder than $E^{-2}$, CTA has the potential to robustly discover the IC signature of a putative Galactic bulge MSP population sufficient to explain the GCE for $e^\pm$ injection efficiencies in the range $\approx 2.9-74.1\%$, or higher, depending on the level of mismodeling of the Galactic diffuse emission components. On the other hand, for spectra softer than $E^{-2.5}$, a reliable CTA detection would require an unphysically large $e^\pm$ injection efficiency of $\gtrsim 158\%$. However, even this pessimistic conclusion may be avoided in the plausible event that MSP observational and/or modeling uncertainties can be reduced. We further find that, in the event that an IC signal were detected, CTA can successfully discriminate between an MSP and a dark matter origin for the radiating $e^\pm$. △ Less

Submitted 26 June, 2021; v1 submitted 10 February, 2021; originally announced February 2021.

Comments: 21+6 pages, 12+5 figures. Updated to match version soon to be published by MNRAS

Journal ref: MNRAS 506, 1741-1760 (2021)

Evidence for a high-energy tail in the gamma-ray spectra of globular clusters

Authors: Deheng Song, Oscar Macias, Shunsaku Horiuchi, Roland M. Crocker, David M. Nataf

Abstract: Millisecond pulsars are very likely the main source of gamma-ray emission from globular clusters. However, the relative contributions of two separate emission processes--curvature radiation from millisecond pulsar magnetospheres vs. inverse Compton emission from relativistic pairs launched into the globular cluster environment by millisecond pulsars--have long been unclear. To address this, we sea… ▽ More Millisecond pulsars are very likely the main source of gamma-ray emission from globular clusters. However, the relative contributions of two separate emission processes--curvature radiation from millisecond pulsar magnetospheres vs. inverse Compton emission from relativistic pairs launched into the globular cluster environment by millisecond pulsars--have long been unclear. To address this, we search for evidence of inverse Compton emission in 8-year $\textit{Fermi}$-LAT data from the directions of 157 Milky Way globular clusters. We find a mildly statistically significant (3.8$σ$) correlation between the measured globular cluster gamma-ray luminosities and their photon field energy densities. However, this may also be explained by a hidden correlation between the photon field densities and the stellar encounter rates of globular clusters. Analysed $\textit{in toto}$, we demonstrate that the gamma-ray emission of globular clusters can be resolved spectrally into two components: i) an exponentially cut-off power law and ii) a pure power law. The latter component--which we uncover at a significance of 8.2$σ$--has a power index of 2.79 $\pm$ 0.25. It is most naturally interpreted as inverse Compton emission by cosmic-ray electrons and positrons injected by millisecond pulsars. We find the luminosity of this power-law component is comparable to, or slightly smaller than, the luminosity of the curved component, suggesting the fraction of millisecond pulsar spin-down luminosity into relativistic leptons is similar to the fraction of the spin-down luminosity into prompt magnetospheric radiation. △ Less

Submitted 26 September, 2021; v1 submitted 29 January, 2021; originally announced February 2021.

Comments: 16+9 pages, 13+4 figures, 4+1 tables. match the published version in MNRAS

Journal ref: MNRAS 507, 5161-5176 (2021)

A supernova remnant associated with a nascent black hole low-mass X-ray binary

Authors: N. I. Maxted, A. J. Ruiter, K. Belczynski, I. R. Seitenzahl, R. M. Crocker

Abstract: Energy released when the core of a high-mass star collapses into a black hole often powers an explosion that creates a supernova remnant. Black holes have limited windows of observability, and consequently are rarely identified in association with supernova remnants. Analysing multi-messenger data, we show that MAXI J1535-571 is the black hole produced in the stellar explosion that gave rise to th… ▽ More Energy released when the core of a high-mass star collapses into a black hole often powers an explosion that creates a supernova remnant. Black holes have limited windows of observability, and consequently are rarely identified in association with supernova remnants. Analysing multi-messenger data, we show that MAXI J1535-571 is the black hole produced in the stellar explosion that gave rise to the supernova remnant G323.7-1.0, making it the first case of an association between a black hole low-mass X-ray binary and a supernova remnant. Given this connection, we can infer from our modelling that the progenitor system was a close binary whose primary star had an initial mass of approx. 23-35 solar masses with a companion star about 10 times less massive. △ Less

Submitted 28 October, 2020; originally announced October 2020.

Comments: 4 pages text + 2 figures (8 pages total). To be submitted. Comments and feedback welcome!

Cosmic rays across the star-forming galaxy sequence. II: Stability limits and the onset of cosmic ray-driven outflows

Authors: Roland M. Crocker, Mark R. Krumholz, Todd A. Thompson

Abstract: Cosmic rays (CRs) are a plausible mechanism for launching winds of cool material from the discs of star-forming galaxies. However, there is no consensus on what types of galaxies likely host CR-driven winds, or what role these winds might play in regulating galaxies' star formation rates. Using a detailed treatment of the transport and losses of hadronic CRs developed in the previous paper in this… ▽ More Cosmic rays (CRs) are a plausible mechanism for launching winds of cool material from the discs of star-forming galaxies. However, there is no consensus on what types of galaxies likely host CR-driven winds, or what role these winds might play in regulating galaxies' star formation rates. Using a detailed treatment of the transport and losses of hadronic CRs developed in the previous paper in this series, here we develop a semi-analytic model that allows us to assess the viability of using CRs to launch cool winds from galactic discs. In particular, we determine the critical CR fluxes -- and corresponding star formation rate surface densities -- above which hydrostatic equilibrium within a given galaxy is precluded because CRs drive the gas off in a wind or otherwise render it unstable. We show that, for star-forming galaxies with lower gas surface densities typical of the Galaxy and local dwarfs, the locus of this CR stability curve patrols the high side of the observed distribution of galaxies in the Kennicutt-Schmidt parameter space of star formation rate versus gas surface density. However, hadronic losses render CRs unable to drive winds in galaxies with higher surface densities. Our results show that quiescent, low surface density galaxies like the Milky Way are poised on the cusp of instability, such that small changes to ISM parameters can lead to the launching of CR-driven outflows, and we suggest that, as a result, CR feedback sets an ultimate limit to the star formation efficiency of most modern galaxies. △ Less

Submitted 30 June, 2020; v1 submitted 29 June, 2020; originally announced June 2020.

Comments: 14 pages, 5 figures. Submitted to MNRAS. Comments welcome. (Legacy version of paper mistakenly appended to .pdf of original arXiv posting removed; main paper unchanged.)

Cosmic rays across the star-forming galaxy sequence. I: Cosmic ray pressures and calorimetry

Authors: Roland M. Crocker, Mark R. Krumholz, Todd A. Thompson

Abstract: In the Milky Way, cosmic rays (CRs) are dynamically important in the interstellar medium, contribute to hydrostatic balance, and may help regulate star formation. However, we know far less about the importance of CRs in galaxies whose gas content or star formation rate differ significantly from those of the Milky Way. Here we construct self-consistent models for hadronic CR transport, losses, and… ▽ More In the Milky Way, cosmic rays (CRs) are dynamically important in the interstellar medium, contribute to hydrostatic balance, and may help regulate star formation. However, we know far less about the importance of CRs in galaxies whose gas content or star formation rate differ significantly from those of the Milky Way. Here we construct self-consistent models for hadronic CR transport, losses, and contribution to pressure balance as a function of galaxy properties, covering a broad range of parameters from dwarfs to extreme starbursts. While the CR energy density increases from $\sim 1$ eV cm$^{-3}$ to $\sim 1$ keV cm$^{-3}$ over the range from sub-Milky Way dwarfs to bright starbursts, strong hadronic losses render CRs increasingly unimportant dynamically as the star formation rate surface density increases. In Milky Way-like systems, CR pressure is typically comparable to turbulent gas and magnetic pressure at the galactic midplane, but the ratio of CR pressure to gas pressure drops to $\sim 10^{-3}$ in dense starbursts. Galaxies also become increasingly CR calorimetric and gamma-ray bright in this limit. The degree of calorimetry at fixed galaxy properties is sensitive to the assumed model for CR transport, and in particular to the time CRs spend interacting with neutral ISM, where they undergo strong streaming losses. We also find that in some regimes of parameter space hydrostatic equilibrium discs cannot exist, and in Paper II of this series we use this result to derive a critical surface in the plane of star formation surface density and gas surface density beyond which CRs may drive large-scale galactic winds. △ Less

Submitted 18 January, 2021; v1 submitted 29 June, 2020; originally announced June 2020.

Comments: Accepted for publication in MNRAS. 23 pages, 10 figures. Changes from v1 include a new bibliography of symbols and a new appendix

Cosmic ray transport in starburst galaxies

Authors: Mark R. Krumholz, Roland M. Crocker, Siyao Xu, A. Lazarian, M. T. Rosevear, Jasper Bedwell-Wilson

Abstract: Starburst galaxies are efficient $γ$-ray producers, because their high supernova rates generate copious cosmic ray (CR) protons, and their high gas densities act as thick targets off which these protons can produce neutral pions and thence $γ$-rays. In this paper we present a first-principles calculation of the mechanisms by which CRs propagate through such environments, combining astrochemical mo… ▽ More Starburst galaxies are efficient $γ$-ray producers, because their high supernova rates generate copious cosmic ray (CR) protons, and their high gas densities act as thick targets off which these protons can produce neutral pions and thence $γ$-rays. In this paper we present a first-principles calculation of the mechanisms by which CRs propagate through such environments, combining astrochemical models with analysis of turbulence in weakly ionised plasma. We show that CRs cannot scatter off the strong large-scale turbulence found in starbursts, because efficient ion-neutral damping prevents such turbulence from cascading down to the scales of CR gyroradii. Instead, CRs stream along field lines at a rate determined by the competition between streaming instability and ion-neutral damping, leading to transport via a process of field line random walk. This results in an effective diffusion coefficient that is nearly energy-independent up to CR energies of $\sim 1$ TeV. We apply our computed diffusion coefficient to a simple model of CR escape and loss, and show that the resulting $γ$-ray spectra are in good agreement with the observed spectra of the starbursts NGC 253, M82, and Arp 220. In particular, our model reproduces these galaxies' relatively hard GeV $γ$-ray spectra and softer TeV spectra without the need for any fine-tuning of advective escape times or the shape of the CR injection spectrum. △ Less

Submitted 18 February, 2020; v1 submitted 21 November, 2019; originally announced November 2019.

Comments: 19 pages, 7 figures, MNRAS in press; minor text changes from previous version, no change in primary results

Constraints on positron annihilation kinematics in the inner Galaxy

Authors: Thomas Siegert, Roland M. Crocker, Roland Diehl, Martin G. H. Krause, Fiona H. Panther, Moritz M. M. Pleintinger, Christoph Weinberger

Abstract: The annihilation of cosmic positrons ($e^+$) with electrons in the interstellar medium (ISM) results in the strongest persistent gamma-ray line signal in the sky. For 50 years, this 511 keV emission has puzzled observers and theoreticians. A key issue for understanding $e^+$-astrophysics is found in cosmic-ray propagation, especially at low kinetic energies (< 10 MeV). We want to shed light on how… ▽ More The annihilation of cosmic positrons ($e^+$) with electrons in the interstellar medium (ISM) results in the strongest persistent gamma-ray line signal in the sky. For 50 years, this 511 keV emission has puzzled observers and theoreticians. A key issue for understanding $e^+$-astrophysics is found in cosmic-ray propagation, especially at low kinetic energies (< 10 MeV). We want to shed light on how $e^+$s propagate and the resulting morphology of the emission. We approach this "positron puzzle" by inferring kinematic information of the 511 keV line in the inner radian of the Galaxy. This constrains propagation scenarios and source populations. By dissecting the 511 keV emission as measured with INTEGRAL/SPI, we derive spectra for individual regions in the sky. The centroid energies are converted into Doppler-shifts, representing the line-of-sight velocity along different longitudes. This results in a longitude-velocity diagram of $e^+$-annihilation. We also determine Doppler-broadenings to study annihilation conditions as they vary across the Galaxy. We find line-of-sight velocities in the 511 keV line that are consistent with zero, as well as with galactic rotation from CO measurements, and measurements of radioactive Al-26. The velocity gradient in the inner 60 deg is determined to be $4\pm6$ km/s/deg. The 511 keV line width is constant as a function of longitude at $2.43\pm0.14$ keV. The positronium fraction is found to be 1.0 along the galactic plane. The weak signals in the disk leave open the question whether $e^+$-annihilation is associated with the high velocities seen in Al-26 or rather with ordinarily rotating components of the Galaxy's ISM. We confirm previous results that $e^+$s are slowed down to the 10 eV energy scale before annihilation, and constrain bulk Doppler-broadening contributions to <1.25 keV. Consequently, the true annihilation conditions remain unclear. △ Less

Submitted 2 June, 2019; originally announced June 2019.

Comments: 8 pages, 5 main figures, 7 pages appendix with 6 auxiliary figures, accepted for publication in A&A

Journal ref: A&A 627, A126 (2019)

SN1991bg-like supernovae are associated with old stellar populations

Authors: Fiona H. Panther, Ivo R. Seitenzahl, Ashley J. Ruiter, Roland M. Crocker, Chris Lidman, Ella Xi Wang, Brad E. Tucker, Brent Groves

Abstract: SN1991bg-like supernovae are a distinct subclass of thermonuclear supernovae (SNe Ia). Their spectral and photometric peculiarities indicate their progenitors and explosion mechanism differ from `normal' SNe Ia. One method of determining information about supernova progenitors we cannot directly observe is to observe the stellar population adjacent to the apparent supernova explosion site to infer… ▽ More SN1991bg-like supernovae are a distinct subclass of thermonuclear supernovae (SNe Ia). Their spectral and photometric peculiarities indicate their progenitors and explosion mechanism differ from `normal' SNe Ia. One method of determining information about supernova progenitors we cannot directly observe is to observe the stellar population adjacent to the apparent supernova explosion site to infer the distribution of stellar population ages and metallicities. We obtain integral field observations and analyse the spectra extracted from regions of projected radius $\sim\,\mathrm{kpc}$ about the apparent SN explosion site for 11 91bg-like SNe in both early- and late-type galaxies. We utilize full-spectrum spectral fitting to determine the ages and metallicities of the stellar population within the aperture. We find that the majority of the stellar populations that hosted 91bg-like supernovae have little recent star formation. The ages of the stellar populations suggest that that 91bg-like SN progenitors explode after delay times of $>6\,\mathrm{Gyr}$, much longer than the typical delay time of normal SNe Ia, which peaks at $\sim 1\,\mathrm{Gyr}$. △ Less

Submitted 23 April, 2019; originally announced April 2019.

Comments: 12 pages, 3 figures, 3 tables, submitted to Publications of the Astronomical Society of Australia

Journal ref: Publ. Astron. Soc. Aust. 36 (2019) e031

S-band Polarization All Sky Survey (S-PASS): survey description and maps

Authors: E. Carretti, M. Haverkorn, L. Staveley-Smith, G. Bernardi, B. M. Gaensler, M. J. Kesteven, S. Poppi, S. Brown, R. M. Crocker, C. Purcell, D. H. F. M. Schnitzler, X. Sun

Abstract: We present the S-Band Polarization All Sky Survey (S-PASS), a survey of polarized radio emission over the southern sky at Dec~$< -1^\circ$ taken with the Parkes radio telescope at 2.3~GHz. The main aim was to observe at a frequency high enough to avoid strong depolarization at intermediate Galactic latitudes (still present at 1.4 GHz) to study Galactic magnetism, but low enough to retain ample Sig… ▽ More We present the S-Band Polarization All Sky Survey (S-PASS), a survey of polarized radio emission over the southern sky at Dec~$< -1^\circ$ taken with the Parkes radio telescope at 2.3~GHz. The main aim was to observe at a frequency high enough to avoid strong depolarization at intermediate Galactic latitudes (still present at 1.4 GHz) to study Galactic magnetism, but low enough to retain ample Signal-to-Noise ratio (S/N) at high latitudes for extragalactic and cosmological science. We developed a new scanning strategy based on long azimuth scans, and a corresponding map-making procedure to make recovery of the overall mean signal of Stokes $Q$ and $U$ possible, a long-standing problem with polarization observations. We describe the scanning strategy, map-making procedure, and validation tests. The overall mean signal is recovered with a precision better than 0.5\%. The maps have a mean sensitivity of 0.81 mK on beam--size scales and show clear polarized signals, typically to within a few degrees of the Galactic plane, with ample S/N everywhere (the typical signal in low emission regions is 13 mK, and 98.6\% of the pixels have S/N $> 3$). The largest depolarization areas are in the inner Galaxy, associated with the Sagittarius Arm. We have also computed a Rotation Measure map combining S-PASS with archival data from the WMAP and Planck experiments. A Stokes $I$ map has been generated, with a sensitivity limited to the confusion level of 9 mK. △ Less

Submitted 22 March, 2019; originally announced March 2019.

Comments: Accepted for publication on MNRAS. Maps are available for download at the website indicated in the manuscript

Strong Evidence that the Galactic Bulge is Shining in Gamma Rays

Authors: Oscar Macias, Shunsaku Horiuchi, Manoj Kaplinghat, Chris Gordon, Roland M. Crocker, David M. Nataf

Abstract: There is growing evidence that the Galactic Center Excess identified in the $\textit{Fermi}$-LAT gamma-ray data arises from a population of faint astrophysical sources. We provide compelling supporting evidence by showing that the morphology of the excess traces the stellar over-density of the Galactic bulge. By adopting a template of the bulge stars obtained from a triaxial 3D fit to the diffuse… ▽ More There is growing evidence that the Galactic Center Excess identified in the $\textit{Fermi}$-LAT gamma-ray data arises from a population of faint astrophysical sources. We provide compelling supporting evidence by showing that the morphology of the excess traces the stellar over-density of the Galactic bulge. By adopting a template of the bulge stars obtained from a triaxial 3D fit to the diffuse near-infrared emission, we show that it is detected at high significance. The significance deteriorates when either the position or the orientation of the template is artificially shifted, supporting the correlation of the gamma-ray data with the Galactic bulge. In deriving these results, we have used more sophisticated templates at low-latitudes for the $\textit{Fermi}$ bubbles compared to previous work and the three-dimensional Inverse Compton (IC) maps recently released by the ${\tt GALPROP}$ team. Our results provide strong constraints on Millisecond Pulsar (MSP) formation scenarios proposed to explain the excess. We find that an $\textit{admixture formation}$ scenario, in which some of the relevant binaries are $\textit{primordial}$ and the rest are formed $\textit{dynamically}$, is preferred over a primordial-only formation scenario at $7.6σ$ confidence level. Our detailed morphological analysis also disfavors models of the disrupted globular clusters scenario that predict a spherically symmetric distribution of MSPs in the Galactic bulge. For the first time, we report evidence of a high energy tail in the nuclear bulge spectrum that could be the result of IC emission from electrons and positrons injected by a population of MSPs and star formation activity from the same site. △ Less

Submitted 20 September, 2019; v1 submitted 12 January, 2019; originally announced January 2019.

Comments: 21 pages, 13 figures, V2: Minor changes to match submitted version, V3: matches JCAP published version

Journal ref: JCAP09(2019)042

Radiation Pressure Limits on the Star Formation Efficiency and Surface Density of Compact Stellar Systems

Authors: Roland M. Crocker, Mark R. Krumholz, Todd A. Thompson, Holger Baumgardt, Dougal Mackey

Abstract: The large columns of dusty gas enshrouding and fuelling star-formation in young, massive stellar clusters may render such systems optically thick to radiation well into the infrared. This raises the prospect that both "direct" radiation pressure produced by absorption of photons leaving stellar surfaces and "indirect" radiation pressure from photons absorbed and then re-emitted by dust grains may… ▽ More The large columns of dusty gas enshrouding and fuelling star-formation in young, massive stellar clusters may render such systems optically thick to radiation well into the infrared. This raises the prospect that both "direct" radiation pressure produced by absorption of photons leaving stellar surfaces and "indirect" radiation pressure from photons absorbed and then re-emitted by dust grains may be important sources of feedback in such systems. Here we evaluate this possibility by deriving the conditions under which a spheroidal, self-gravitating, mixed gas-star cloud can avoid catastrophic disruption by the combined effects of direct and indirect radiation pressure. We show that radiation pressure sets a maximum star cluster formation efficiency of $ε_{\rm max} \sim 0.9$ at a (very large) gas surface density of $\sim 10^5 M_\odot$ pc$^{-2} (Z_\odot/Z) \simeq 20$ g cm$^{-2} (Z_\odot/Z)$, but that gas clouds above this limit undergo significant radiation-driven expansion during star formation, leading to a maximum stellar surface density very near this value for all star clusters. Data on the central surface mass density of compact stellar systems, while sparse and partly confused by dynamical effects, are broadly consistent with the existence of a metallicity-dependent upper-limit comparable to this value. Our results imply that this limit may preclude the formation of the progenitors of intermediate-mass black holes for systems with $Z \gtrsim 0.2 Z_\odot$. △ Less

Submitted 5 October, 2018; v1 submitted 6 August, 2018; originally announced August 2018.

Comments: 13 pages, 7 figures. Accepted for publication in MNRAS

Interstellar magnetic cannon targeting the Galactic halo : A young bubble at the origin of the Ophiuchus and Lupus molecular complexes

Authors: J. -F. Robitaille, A. M. M. Scaife, E. Carretti, M. Haverkorn, R. M. Crocker, M. J. Kesteven, S. Poppi, L. Staveley-Smith

Abstract: We report the detection of a new Galactic bubble at the interface between the halo and the Galactic disc. We suggest that the nearby Lupus complex and parts of the Ophiuchus complex constitute the denser parts of the structure. This young bubble, < 3 Myr old, could be the remnant of a supernova and expands inside a larger HI loop that has been created by the outflows of the Upper Scorpius OB assoc… ▽ More We report the detection of a new Galactic bubble at the interface between the halo and the Galactic disc. We suggest that the nearby Lupus complex and parts of the Ophiuchus complex constitute the denser parts of the structure. This young bubble, < 3 Myr old, could be the remnant of a supernova and expands inside a larger HI loop that has been created by the outflows of the Upper Scorpius OB association. An HI cavity filled with hot X-ray gas is associated with the structure, which is consistent with the Galactic chimney scenario. The X-ray emission extends beyond the west and north-west edges of the bubble, suggesting that hot gas outflows are breaching the cavity, possibly through the fragmented Lupus complex. Analyses of the polarised radio synchrotron and of the polarised dust emission of the region suggest the connection of the Galactic centre spur with the young Galactic bubble. A distribution of HI clumps that spatially corresponds well to the cavity boundaries was found at V_LSR ~ -100 km/s. Some of these HI clumps are forming jets, which may arise from the fragmented part of the bubble. We suggest that these clumps might be `dripping' cold clouds from the shell walls inside the cavity that is filled with hot ionised gas. It is possible that some of these clumps are magnetised and were then accelerated by the compressed magnetic field at the edge of the cavity. Such a mechanism would challenge the Galactic accretion and fountain model, where high-velocity clouds are considered to be formed at high Galactic latitude from hot gas flows from the Galactic plane. △ Less

Submitted 11 July, 2018; originally announced July 2018.

Comments: 16 pages, 16 figures

Journal ref: A&A 617, A101 (2018)

The effect of positron-alkali metal atom interactions in the diffuse ISM

Authors: Fiona H. Panther, Ivo R. Seitenzahl, Roland M. Crocker, Joshua R. Machacek, Dan J. Murtagh, Thomas Siegert, Roland Diehl

Abstract: In the Milky Way galaxy, positrons, which are responsible for the diffuse $511\,\mathrm{keV}$ gamma ray emission observed by space-based gamma ray observatories, are thought to annihilate predominantly through charge exchange interactions with neutral hydrogen. These charge exchange interactions can only take place if positrons have energies greater than $6.8\,\mathrm{eV}$, the minimum energy requ… ▽ More In the Milky Way galaxy, positrons, which are responsible for the diffuse $511\,\mathrm{keV}$ gamma ray emission observed by space-based gamma ray observatories, are thought to annihilate predominantly through charge exchange interactions with neutral hydrogen. These charge exchange interactions can only take place if positrons have energies greater than $6.8\,\mathrm{eV}$, the minimum energy required to liberate the electron bound to the hydrogen atom and then form positronium, a short-lived bound state composed of a positron-electron pair. Here we demonstrate the importance of positron interactions with neutral alkali metals in the warm interstellar medium (ISM). Positrons may undergo charge exchange with these atoms at any energy. In particular, we show that including positron interactions with sodium at solar abundance in the warm ISM can significantly reduce the annihilation timescale of positrons with energies below $6.8\,\mathrm{eV}$ by at least an order of magnitude. We show that including these interactions in our understanding of positron annihilation in the Milky Way rules out the idea that the number of positrons in the Galactic ISM could be maintained in steady state by injection events occurring at a typical periodicity $>\mathrm{Myr}$. △ Less

Submitted 5 July, 2018; originally announced July 2018.

Comments: 8 pages, 3 figures. Accepted for publication Physical Review D

The Maximum Flux of Star-Forming Galaxies

Authors: Roland M. Crocker, Mark R. Krumholz, Todd A. Thompson, Julie Clutterbuck

Abstract: The importance of radiation pressure feedback in galaxy formation has been extensively debated over the last decade. The regime of greatest uncertainty is in the most actively star-forming galaxies, where large dust columns can potentially produce a dust-reprocessed infrared radiation field with enough pressure to drive turbulence or eject material. Here we derive the conditions under which a self… ▽ More The importance of radiation pressure feedback in galaxy formation has been extensively debated over the last decade. The regime of greatest uncertainty is in the most actively star-forming galaxies, where large dust columns can potentially produce a dust-reprocessed infrared radiation field with enough pressure to drive turbulence or eject material. Here we derive the conditions under which a self-gravitating, mixed gas-star disc can remain hydrostatic despite trapped radiation pressure. Consistently taking into account the self-gravity of the medium, the star- and dust-to-gas ratios, and the effects of turbulent motions not driven by radiation, we show that galaxies can achieve a maximum Eddington-limited star formation rate per unit area $\dotΣ_{\rm *,crit} \sim 10^3 M_{\odot}$ pc$^{-2}$ Myr$^{-1}$, corresponding to a critical flux of $F_{\rm *,crit} \sim 10^{13} L_{\odot}$ kpc$^{-2}$ similar to previous estimates; higher fluxes eject mass in bulk, halting further star formation. Conversely, we show that in galaxies below this limit, our one-dimensional models imply simple vertical hydrostatic equilibrium and that radiation pressure is ineffective at driving turbulence or ejecting matter. Because the vast majority of star-forming galaxies lie below the maximum limit for typical dust-to-gas ratios, we conclude that infrared radiation pressure is likely unimportant for all but the most extreme systems on galaxy-wide scales. Thus, while radiation pressure does not explain the Kennicutt-Schmidt relation, it does impose an upper truncation on it. Our predicted truncation is in good agreement with the highest observed gas and star formation rate surface densities found both locally and at high redshift. △ Less

Submitted 21 April, 2018; v1 submitted 8 February, 2018; originally announced February 2018.

Comments: Version accepted for publication in MNRAS. 12 pages, 8 figures. New appendix on photon tiring

On the formation of neutron stars via accretion-induced collapse in binaries

Authors: A. J. Ruiter, L. Ferrario, K. Belczynski, I. R. Seitenzahl, R. M. Crocker, A. I. Karakas

Abstract: We investigate evolutionary pathways leading to neutron star formation through the collapse of oxygen-neon white dwarf (ONe WD) stars in interacting binaries. We consider (1) non-dynamical mass transfer where an ONe WD approaches the Chandrasekhar mass leading to accretion-induced collapse (AIC) and (2) dynamical timescale merger-induced collapse (MIC) between an ONe WD and another WD. We present… ▽ More We investigate evolutionary pathways leading to neutron star formation through the collapse of oxygen-neon white dwarf (ONe WD) stars in interacting binaries. We consider (1) non-dynamical mass transfer where an ONe WD approaches the Chandrasekhar mass leading to accretion-induced collapse (AIC) and (2) dynamical timescale merger-induced collapse (MIC) between an ONe WD and another WD. We present rates, delay times, and progenitor properties for two different treatments of common envelope evolution. We show that AIC neutron stars are formed via many different channels and the most dominant channel depends on the adopted common envelope physics. Most AIC and MIC neutron stars are born shortly after star formation, though some have delay times >10 Gyr. The shortest delay time (25-50 Myr) AIC neutron stars have stripped-envelope, compact, helium-burning star donors, though many prompt AIC neutron stars form via wind-accretion from an asymptotic giant branch star. The longest delay time AIC neutron stars, which may be observed as young milli-second pulsars among globular clusters, have a red giant or main sequence donor at the time of NS formation and will eventually evolve into NS + helium WD binaries. We discuss AIC & MIC binaries as potential gravitational wave sources for LISA. Neutron stars created via AIC undergo a LMXB phase, offering an electromagnetic counterpart for those shortest orbital period sources that LISA could identify. The formation of neutron stars from interacting WDs in binaries is likely to be a key mechanism for the production of LIGO/Virgo gravitational wave sources (NS-NS and BH-NS mergers) in globular clusters. △ Less

Submitted 8 January, 2019; v1 submitted 7 February, 2018; originally announced February 2018.

Comments: Accepted for publication in MNRAS

Positron Annihilation in the Nuclear Outflows of the Milky Way

Authors: Fiona H. Panther, Roland M. Crocker, Yuval Birnboim, Ivo R. Seitenzahl, Ashley J. Ruiter

Abstract: Observations of soft gamma rays emanating from the Milky Way from SPI/\textit{INTEGRAL} reveal the annihilation of $\sim2\times10^{43}$ positrons every second in the Galactic bulge. The origin of these positrons, which annihilate to produce a prominent emission line centered at 511 keV, has remained mysterious since their discovery almost 50 years ago. A plausible origin for the positrons is in as… ▽ More Observations of soft gamma rays emanating from the Milky Way from SPI/\textit{INTEGRAL} reveal the annihilation of $\sim2\times10^{43}$ positrons every second in the Galactic bulge. The origin of these positrons, which annihilate to produce a prominent emission line centered at 511 keV, has remained mysterious since their discovery almost 50 years ago. A plausible origin for the positrons is in association with the intense star formation ongoing in the Galactic center. Moreover, there is strong evidence for a nuclear outflow in the Milky Way. We find that advective transport and subsequent annihilation of positrons in such an outflow cannot simultaneously replicate the observed morphology of positron annihilation in the Galactic bulge and satisfy the requirement that $90$ per cent of positrons annihilate once the outflow has cooled to $10^4\,\mathrm{K}$. △ Less

Submitted 23 November, 2017; v1 submitted 6 October, 2017; originally announced October 2017.

Comments: 5 pages, 1 figure. Accepted for publication in Monthly Notices of the Royal Astronomical Society Letters

A Unified Model for Galactic Discs: Star Formation, Turbulence Driving, and Mass Transport

Authors: Mark R. Krumholz, Blakesley Burkhart, John C. Forbes, Roland M. Crocker

Abstract: We introduce a new model for the structure and evolution of the gas in galactic discs. In the model the gas is in vertical pressure and energy balance. Star formation feedback injects energy and momentum, and non-axisymmetric torques prevent the gas from becoming more than marginally gravitationally unstable. From these assumptions we derive the relationship between galaxies' bulk properties (gas… ▽ More We introduce a new model for the structure and evolution of the gas in galactic discs. In the model the gas is in vertical pressure and energy balance. Star formation feedback injects energy and momentum, and non-axisymmetric torques prevent the gas from becoming more than marginally gravitationally unstable. From these assumptions we derive the relationship between galaxies' bulk properties (gas surface density, stellar content, and rotation curve) and their star formation rates, gas velocity dispersions, and rates of radial inflow. We show that the turbulence in discs can be powered primarily by star formation feedback, radial transport, or a combination of the two. In contrast to models that omit either radial transport or star formation feedback, the predictions of this model yield excellent agreement with a wide range of observations, including the star formation law measured in both spatially resolved and unresolved data, the correlation between galaxies' star formation rates and velocity dispersions, and observed rates of radial inflow. The agreement holds across a wide range of galaxy mass and type, from local dwarfs to extreme starbursts to high-redshifts discs. We apply the model to galaxies on the star-forming main sequence, and show that it predicts a transition from mostly gravity-driven turbulence at high redshift to star formation-driven turbulence at low redshift. This transition, and the changes in mass transport rates that it produces, naturally explain why galaxy bulges tend to form at high redshift and discs at lower redshift, and why galaxies tend to quench inside-out. △ Less

Submitted 29 March, 2018; v1 submitted 31 May, 2017; originally announced June 2017.

Comments: 26 pages, 7 figures, MNRAS in press; minor changes from previous version; all plotting scripts and data available from https://bitbucket.org/krumholz/kbfc17

Galactic Bulge Preferred Over Dark Matter for the Galactic Center Gamma-Ray Excess

Authors: Oscar Macias, Chris Gordon, Roland M. Crocker, Brenna Coleman, Dylan Paterson, Shunsaku Horiuchi, Martin Pohl

Abstract: An anomalous gamma-ray excess emission has been found in Fermi Large Area Telescope data covering the centre of the Galaxy. Several theories have been proposed for this `Galactic Centre Excess'. They include self-annihilation of dark matter particles, an unresolved population of millisecond pulsars, an unresolved population of young pulsars, or a series of burst events. Here we report on a new a… ▽ More An anomalous gamma-ray excess emission has been found in Fermi Large Area Telescope data covering the centre of the Galaxy. Several theories have been proposed for this `Galactic Centre Excess'. They include self-annihilation of dark matter particles, an unresolved population of millisecond pulsars, an unresolved population of young pulsars, or a series of burst events. Here we report on a new analysis that exploits hydrodynamical modelling to register the position of interstellar gas associated with diffuse Galactic gamma-ray emission. We find evidence that the Galactic Centre Excess gamma rays are statistically better described by the stellar over-density in the Galactic bulge and the nuclear stellar bulge, rather than a spherical excess. Given its non-spherical nature, we argue that the Galactic Centre Excess is not a dark matter phenomenon but rather associated with the stellar population of the Galactic bulge and nuclear bulge. △ Less

Submitted 25 September, 2018; v1 submitted 20 November, 2016; originally announced November 2016.

Comments: 21 pages, 7 figures, V2: Minor modifications, main conclusions unchanged. V3: version submitted to Nature Astronomy. V4: matches Nature Astronomy published version. The hydrodynamical interstellar maps used in this work are available at https://github.com/chrisgordon1/galactic_bulge

Journal ref: Nature Astronomy (2018)

The link between solenoidal turbulence and slow star formation in G0.253+0.016

Authors: C. Federrath, J. M. Rathborne, S. N. Longmore, J. M. D. Kruijssen, J. Bally, Y. Contreras, R. M. Crocker, G. Garay, J. M. Jackson, L. Testi, A. J. Walsh

Abstract: Star formation in the Galactic disc is primarily controlled by gravity, turbulence, and magnetic fields. It is not clear that this also applies to star formation near the Galactic Centre. Here we determine the turbulence and star formation in the CMZ cloud G0.253+0.016. Using maps of 3mm dust emission and HNCO intensity-weighted velocity obtained with ALMA, we measure the volume-density variance… ▽ More Star formation in the Galactic disc is primarily controlled by gravity, turbulence, and magnetic fields. It is not clear that this also applies to star formation near the Galactic Centre. Here we determine the turbulence and star formation in the CMZ cloud G0.253+0.016. Using maps of 3mm dust emission and HNCO intensity-weighted velocity obtained with ALMA, we measure the volume-density variance $σ_{ρ/ρ_0} = 1.3 \pm 0.5$ and turbulent Mach number $\mathcal{M} = 11 \pm 3$. Combining these with turbulence simulations to constrain the plasma $β= 0.34 \pm 0.35$, we reconstruct the turbulence driving parameter $b = 0.22 \pm 0.12$ in G0.253+0.016. This low value of $b$ indicates solenoidal (divergence-free) driving of the turbulence in G0.253+0.016. By contrast, typical clouds in the Milky Way disc and spiral arms have a significant compressive (curl-free) driving component ($b > 0.4$). We speculate that shear causes the solenoidal driving in G0.253+0.016 and show that this may reduce the star formation rate by a factor of 7 compared to nearby clouds. △ Less

Submitted 27 September, 2016; originally announced September 2016.

Comments: IAU Symposium 322 "The Multi-Messenger Astrophysics of the Galactic Centre", eds. R. M. Crocker, S. N. Longmore, G. V. Bicknell

The link between turbulence, magnetic fields, filaments, and star formation in the Central Molecular Zone cloud G0.253+0.016

Authors: C. Federrath, J. M. Rathborne, S. N. Longmore, J. M. D. Kruijssen, J. Bally, Y. Contreras, R. M. Crocker, G. Garay, J. M. Jackson, L. Testi, A. J. Walsh

Abstract: Star formation is primarily controlled by the interplay between gravity, turbulence, and magnetic fields. However, the turbulence and magnetic fields in molecular clouds near the Galactic Center may differ substantially from spiral-arm clouds. Here we determine the physical parameters of the central molecular zone (CMZ) cloud G0.253+0.016, its turbulence, magnetic field and filamentary structure.… ▽ More Star formation is primarily controlled by the interplay between gravity, turbulence, and magnetic fields. However, the turbulence and magnetic fields in molecular clouds near the Galactic Center may differ substantially from spiral-arm clouds. Here we determine the physical parameters of the central molecular zone (CMZ) cloud G0.253+0.016, its turbulence, magnetic field and filamentary structure. Using column-density maps based on dust-continuum emission observations with ALMA+Herschel, we identify filaments and show that at least one dense core is located along them. We measure the filament width W_fil=0.17$\pm$0.08pc and the sonic scale λ_sonic=0.15$\pm$0.11pc of the turbulence, and find W_fil~λ_sonic. A strong velocity gradient is seen in the HNCO intensity-weighted velocity maps obtained with ALMA+Mopra, which is likely caused by large-scale shearing of G0.253+0.016, producing a wide double-peaked velocity PDF. After subtracting the gradient to isolate the turbulent motions, we find a nearly Gaussian velocity PDF typical for turbulence. We measure the total and turbulent velocity dispersion, 8.8$\pm$0.2km/s and 3.9$\pm$0.1km/s, respectively. Using magnetohydrodynamical simulations, we find that G0.253+0.016's turbulent magnetic field B_turb=130$\pm$50$μ$G is only ~1/10 of the ordered field component. Combining these measurements, we reconstruct the dominant turbulence driving mode in G0.253+0.016 and find a driving parameter b=0.22$\pm$0.12, indicating solenoidal (divergence-free) driving. We compare this to spiral-arm clouds, which typically have a significant compressive (curl-free) driving component (b>0.4). Motivated by previous reports of strong shearing motions in the CMZ, we speculate that shear causes the solenoidal driving in G0.253+0.016 and show that this reduces the star formation rate (SFR) by a factor of 6.9 compared to typical nearby clouds. △ Less

Submitted 19 September, 2016; originally announced September 2016.

Comments: 18 pages, 6 figures, 1 table containing new Brick parameters; accepted by ApJ

Journal ref: C. Federrath et al 2016 ApJ 832 143

SN1991bg-like supernovae are a compelling source of most Galactic antimatter

Authors: Fiona H. Panther, Roland M. Crocker, Ivo R. Seitenzahl, Ashley J. Ruiter

Abstract: The Milky Way Galaxy glows with the soft gamma ray emission resulting from the annihilation of $\sim 5 \times 10^{43}$ electron-positron pairs every second. The origin of this vast quantity of antimatter and the peculiar morphology of the 511keV gamma ray line resulting from this annihilation have been the subject of debate for almost half a century. Most obvious positron sources are associated wi… ▽ More The Milky Way Galaxy glows with the soft gamma ray emission resulting from the annihilation of $\sim 5 \times 10^{43}$ electron-positron pairs every second. The origin of this vast quantity of antimatter and the peculiar morphology of the 511keV gamma ray line resulting from this annihilation have been the subject of debate for almost half a century. Most obvious positron sources are associated with star forming regions and cannot explain the rate of positron annihilation in the Galactic bulge, which last saw star formation some $10\,\mathrm{Gyr}$ ago, or else violate stringent constraints on the positron injection energy. Radioactive decay of elements formed in core collapse supernovae (CCSNe) and normal Type Ia supernovae (SNe Ia) could supply positrons matching the injection energy constraints but the distribution of such potential sources does not replicate the required morphology. We show that a single class of peculiar thermonuclear supernova - SN1991bg-like supernovae (SNe 91bg) - can supply the number and distribution of positrons we see annihilating in the Galaxy through the decay of $^{44}$Ti synthesised in these events. Such $^{44}$Ti production simultaneously addresses the observed abundance of $^{44}$Ca, the $^{44}$Ti decay product, in solar system material. △ Less

Submitted 30 November, 2016; v1 submitted 15 September, 2016; originally announced September 2016.

Comments: Accepted for publication in Proceedings of IAU Symposium 322: The Multimessenger Astrophysics of the Galactic Center 4 pages

Diffuse Galactic antimatter from faint thermonuclear supernovae in old stellar populations

Authors: Roland M. Crocker, Ashley J. Ruiter, Ivo R. Seitenzahl, Fiona H. Panther, Stuart Sim, Holger Baumgardt, Anais Moller, David M. Nataf, Lilia Ferrario, J. J. Eldridge, Martin White, Brad E. Tucker, Felix Aharonian

Abstract: Our Galaxy hosts the annihilation of a few $\times 10^{43}$ low-energy positrons every second. Radioactive isotopes capable of supplying such positrons are synthesised in stars, stellar remnants, and supernovae. For decades, however, there has been no positive identification of a main stellar positron source leading to suggestions that many positrons originate from exotic sources like the Galaxy's… ▽ More Our Galaxy hosts the annihilation of a few $\times 10^{43}$ low-energy positrons every second. Radioactive isotopes capable of supplying such positrons are synthesised in stars, stellar remnants, and supernovae. For decades, however, there has been no positive identification of a main stellar positron source leading to suggestions that many positrons originate from exotic sources like the Galaxy's central super-massive black hole or dark matter annihilation. %, but such sources would not explain the recently-detected positron signal from the extended Galactic disk. Here we show that a single type of transient source, deriving from stellar populations of age 3-6 Gyr and yielding ~0.03 $M_\odot$ of the positron emitter $^{44}$Ti, can simultaneously explain the strength and morphology of the Galactic positron annihilation signal and the solar system abundance of the $^{44}$Ti decay product $^{44}$Ca. This transient is likely the merger of two low-mass white dwarfs, observed in external galaxies as the sub-luminous, thermonuclear supernova known as SN1991bg-like. △ Less

Submitted 23 May, 2017; v1 submitted 12 July, 2016; originally announced July 2016.

Comments: 28 pages main text with 4 figures in preprint style; 26 pages of Supplementary Information

Journal ref: Nature Astronomy 1, Article number: 0135 (2017)

A Dynamical Model for Gas Flows, Star Formation, and Nuclear Winds in Galactic Centres

Authors: Mark R. Krumholz., J. M. Diederik Kruijssen, Roland M. Crocker

Abstract: We present a dynamical model for gas transport, star formation, and winds in the nuclear regions of galaxies, focusing on the Milky Way's Central Molecular Zone (CMZ). In our model angular momentum and mass are transported by a combination of gravitational and bar-driven acoustic instabilities. In gravitationally-unstable regions the gas can form stars, and the resulting feedback drives both turbu… ▽ More We present a dynamical model for gas transport, star formation, and winds in the nuclear regions of galaxies, focusing on the Milky Way's Central Molecular Zone (CMZ). In our model angular momentum and mass are transported by a combination of gravitational and bar-driven acoustic instabilities. In gravitationally-unstable regions the gas can form stars, and the resulting feedback drives both turbulence and a wind that ejects mass from the CMZ. We show that the CMZ is in a quasi-steady state where mass deposited at large radii by the bar is transported inward to a star-forming, ring-shaped region at $\sim 100$ pc from the Galactic Centre, where the shear reaches a minimum. This ring undergoes episodic starbursts, with bursts lasting $\sim 5-10$ Myr occurring at $\sim 20-40$ Myr intervals. During quiescence the gas in the ring is not fully cleared, but is driven out of a self-gravitating state by the momentum injected by expanding supernova remnants. Starbursts also drive a wind off the star-forming ring, with a time-averaged mass flux comparable to the star formation rate. We show that our model agrees well with the observed properties of the CMZ, and places it near a star formation minimum within the evolutionary cycle. We argue that such cycles of bursty star formation and winds should be ubiquitous in the nuclei of barred spiral galaxies, and show that the resulting distribution of galactic nuclei on the Kennicutt-Schmidt relation is in good agreement with that observed in nearby galaxies. △ Less

Submitted 6 December, 2016; v1 submitted 10 May, 2016; originally announced May 2016.

Comments: 22 pages, 15 figures, MNRAS in press; code available at https://bitbucket.org/krumholz/cmzsf

Faraday Tomography of the North Polar Spur: Constraints on the distance to the Spur and on the Magnetic Field of the Galaxy

Authors: X. H. Sun, T. L. Landecker, B. M. Gaensler, E. Carretti, W. Reich, J. P. Leahy, N. M. McClure-Griffiths, R. M. Crocker, M. Wolleben, M. Haverkorn, K. A. Douglas, A. D. Gray

Abstract: We present radio continuum and polarization images of the North Polar Spur (NPS) from the Global Magneto-Ionic Medium Survey (GMIMS) conducted with the Dominion Radio Astrophysical Observatory 26-m Telescope. We fit polarization angle versus wavelength squared over 2048 frequency channels from 1280 to 1750 MHz to obtain a Faraday Rotation Measure (RM) map of the NPS. Combining this RM map with a p… ▽ More We present radio continuum and polarization images of the North Polar Spur (NPS) from the Global Magneto-Ionic Medium Survey (GMIMS) conducted with the Dominion Radio Astrophysical Observatory 26-m Telescope. We fit polarization angle versus wavelength squared over 2048 frequency channels from 1280 to 1750 MHz to obtain a Faraday Rotation Measure (RM) map of the NPS. Combining this RM map with a published Faraday depth map of the entire Galaxy in this direction, we derive the Faraday depth introduced by the NPS and the Galactic interstellar medium (ISM) in front of and behind the NPS. The Faraday depth contributed by the NPS is close to zero, indicating that the NPS is an emitting only feature. The Faraday depth caused by the ISM in front of the NPS is consistent with zero at b>50 degree, implying that this part of the NPS is local at a distance of approximately several hundred parsecs. The Faraday depth contributed by the ISM behind the NPS gradually increases with Galactic latitude up to b=44 degree, and decreases at higher Galactic latitudes. This implies that either the part of the NPS at b<44 degree is distant or the NPS is local but there is a sign change of the large-scale magnetic field. If the NPS is local, there is then no evidence for a large-scale anti-symmetry pattern in the Faraday depth of the Milky Way. The Faraday depth introduced by the ISM behind the NPS at latitudes b>50 degree can be explained by including a coherent vertical magnetic field. △ Less

Submitted 16 August, 2015; originally announced August 2015.

Comments: 9 pages, 8 figures, accepted for publication in ApJ. Some figures have been degraded to reduce sizes, for a high resolution version, see http://physics.usyd.edu.au/~xhsun/ms_nps.pdf

A Unified Model of the Fermi Bubbles, Microwave Haze, and Polarized Radio Lobes: Reverse Shocks in the Galactic Center's Giant Outflows

Authors: Roland M. Crocker, Geoffrey V. Bicknell, Andrew M. Taylor, Ettore Carretti

Abstract: The Galactic Center's giant outflows are manifest in three different, non-thermal phenomena: i) the hard-spectrum, γ-ray `Fermi Bubbles' emanating from the nucleus and extending to |b| ~ 50 degrees; ii) the hard-spectrum, total-intensity microwave (~ 20-40 GHz) `Haze' extending to |b| ~ 35 degrees in the lower reaches of the Fermi Bubbles; and iii) the steep spectrum, polarized, `S-PASS' radio (~… ▽ More The Galactic Center's giant outflows are manifest in three different, non-thermal phenomena: i) the hard-spectrum, γ-ray `Fermi Bubbles' emanating from the nucleus and extending to |b| ~ 50 degrees; ii) the hard-spectrum, total-intensity microwave (~ 20-40 GHz) `Haze' extending to |b| ~ 35 degrees in the lower reaches of the Fermi Bubbles; and iii) the steep spectrum, polarized, `S-PASS' radio (~ 2-20 GHz) Lobes that envelop the Bubbles and extend to |b| ~ 60 degrees. We find that the nuclear outflows inflate a genuine bubble in each Galactic hemisphere that has the classical structure, working outwards, of reverse shock, contact discontinuity, and forward shock. Expanding into the finite pressure of the halo and given appreciable cooling and gravitational losses, the contact discontinuity of each bubble is now expanding only very slowly. We find observational signatures in both hemispheres of giant, reverse shocks at heights of ~ 1 kpc above the nucleus; their presence ultimately explains all three of the non-thermal phenomena mentioned above. Synchrotron emission from shock-reaccelerated cosmic-ray electrons explains the spectrum, morphology, and vertical extent of the microwave Haze and the polarized radio Lobes. Collisions between shock-reaccelerated hadrons and denser gas in cooling condensations that form inside the contact discontinuity account for most of the Bubbles' γ-ray emissivity. △ Less

Submitted 17 June, 2015; v1 submitted 23 December, 2014; originally announced December 2014.

Comments: Accepted for publication in the Astrophysical Journal. 14 figures. Qualitative results unchanged from vers. 1 but calculations made more robust and text clarified

Journal ref: ApJ, 808, 107 (2015)

Non-Thermal Radiation from the Inner Galaxy

Authors: Roland M. Crocker

Abstract: I review our current state of knowledge about non-thermal radiation from the Galactic Centre (GC) and Inner Galaxy. Definitionally, the Galactic nucleus is at the bottom of the Galaxy's gravitational well, rendering it a promising region to seek the signatures of dark matter decay or annihilation. It also hosts, however, the Milky Way's resident supermassive black hole and up to 10% of current mas… ▽ More I review our current state of knowledge about non-thermal radiation from the Galactic Centre (GC) and Inner Galaxy. Definitionally, the Galactic nucleus is at the bottom of the Galaxy's gravitational well, rendering it a promising region to seek the signatures of dark matter decay or annihilation. It also hosts, however, the Milky Way's resident supermassive black hole and up to 10% of current massive star formation in the Galaxy. Thus the Galactic nucleus is a dynamic and highly-energized environment implying that extreme caution must be exercised in interpreting any unusual or unexpected signal from (or emerging from) the region as evidence for dark matter-related processes. One spectacular example of an `unexpected' signal is the discovery within the last few years of the `Fermi Bubbles' and, subsequently, their polarised radio counterparts. These giant lobes extend ~7 kpc from the nucleus into both north and south Galactic hemispheres. Hard-spectrum, microwave emission coincident with the lower reaches of the Bubbles has also been detected, first in WMAP, and more recently in Planck data. Debate continues as to the origin of the Bubbles and their multi-wavelength emissions: are they the signatures of relatively recent (in the last ~Myr) activity of the supermassive black hole or, alternatively, nuclear star formation? I will briefly review evidence that points to the latter interpretation. △ Less

Submitted 26 November, 2014; originally announced November 2014.

Comments: 9 pages, no figures; for proceedings of "Cosmic Rays and the InterStellar Medium" PoS(CRISM2014)022

Steady-State Hadronic Gamma-Ray Emission from 100-Myr-Old Fermi Bubbles

Authors: Roland M. Crocker, Geoffrey V. Bicknell, Ettore Carretti, Alex S. Hill, Ralph S. Sutherland

Abstract: The Fermi Bubbles are enigmatic γ-ray features of the Galactic bulge. Both putative activity (within $\sim$ few $\times$ Myr) connected to the Galactic center super-massive black hole and, alternatively, nuclear star formation have been claimed as the energising source of the Bubbles. Likewise, both inverse-Compton emission by non-thermal electrons (`leptonic' models) and collisions between non-th… ▽ More The Fermi Bubbles are enigmatic γ-ray features of the Galactic bulge. Both putative activity (within $\sim$ few $\times$ Myr) connected to the Galactic center super-massive black hole and, alternatively, nuclear star formation have been claimed as the energising source of the Bubbles. Likewise, both inverse-Compton emission by non-thermal electrons (`leptonic' models) and collisions between non-thermal protons and gas (`hadronic' models) have been advanced as the process supplying the Bubbles' γ-ray emission. An issue for any steady state hadronic model is that the very low density of the Bubbles' plasma seems to require that they accumulate protons over a multi-Gyr timescale, much longer than other natural timescales occurring in the problem. Here we present a hadronic model where the timescale for generating the Bubbles' hadronic γ-ray emission is $\sim$ few $\times 10^8$ years. Our model invokes collapse of the Bubbles' thermally-unstable plasma, leading to an accumulation of cosmic rays and magnetic field into localised, warm ($\sim 10^4$ K), and likely filamentary condensations of higher density gas. Under the condition that these filaments are supported by non-thermal pressure, we can predict the hadronic emission from the Bubbles to be $L_γ\simeq 2 \times 10^{37}$ erg/s $\ \dot{M}_\mathrm{in}/(0.1 \ M_{Sun}/$ year $) \ T_\mathrm{FB}^2/(3.5 \times 10^7 K) ^2 M_{fil}/M_{pls}$ ; precisely their observed luminosity (normalizing to the star-formation-driven mass flux into the Bubbles and their measured plasma temperature and adopting the further result that the mass in the filaments, $M_{fil}$ is approximately equal to that of the Bubbles' plasma, $M_{pls}$). △ Less

Submitted 17 August, 2014; v1 submitted 2 December, 2013; originally announced December 2013.

Comments: Accepted for publication in ApJ Letters. Shortened from ver.1 with an important modification to equation 2; qualitative results the same; vers 3: title fixed

Journal ref: ApJ, 791, L20 (2014)

Giant Magnetized Outflows from the Centre of the Milky Way

Authors: E. Carretti, R. M. Crocker, L. Staveley-Smith, M. Haverkorn, C. Purcell, B. M. Gaensler, G. Bernardi, M. J. Kesteven, S. Poppi

Abstract: The nucleus of the Milky Way is known to harbour regions of intense star formation activity as well as a super-massive black hole. Recent Fermi space telescope observations have revealed regions of γ-ray emission reaching far above and below the Galactic Centre, the so-called Fermi bubbles. It is uncertain whether these were generated by nuclear star formation or by quasar-like outbursts of the ce… ▽ More The nucleus of the Milky Way is known to harbour regions of intense star formation activity as well as a super-massive black hole. Recent Fermi space telescope observations have revealed regions of γ-ray emission reaching far above and below the Galactic Centre, the so-called Fermi bubbles. It is uncertain whether these were generated by nuclear star formation or by quasar-like outbursts of the central black hole and no information on the structures' magnetic field has been reported. Here we report on the detection of two giant, linearly-polarized radio Lobes, containing three ridge-like sub-structures, emanating from the Galactic Centre. The Lobes each extend ~60 deg, bear a close correspondence to the Fermi bubbles, are located in the Galactic bulge, and are permeated by strong magnetic fields of up to 15 μG. Our data signal that the radio Lobes originate in a bi-conical, star-formation (rather than black hole) driven outflow from the Galaxy's central 200 pc that transports a massive magnetic energy of ~10^55 erg into the Galactic halo. The ridges wind around this outflow and, we suggest, constitute a `phonographic' record of nuclear star formation activity over at least 10 Myr. △ Less

Submitted 3 January, 2013; originally announced January 2013.

Comments: Published in Nature, main article. 20 pages, 4 figures

Journal ref: Nature, 493, p66 (2013)

Magnetic substructure in the northern Fermi Bubble revealed by polarized WMAP emission

Authors: David I. Jones, Roland M. Crocker, Wolfgang Reich, Jürgen Ott, Felix A. Aharonian

Abstract: We report a correspondence between giant, polarized microwave structures emerging north from the Galactic plane near the Galactic center and a number of GeV gamma-ray features, including the eastern edge of the recently-discovered northern Fermi Bubble. The polarized microwave features also correspond to structures seen in the all-sky 408 MHz total intensity data, including the Galactic center spu… ▽ More We report a correspondence between giant, polarized microwave structures emerging north from the Galactic plane near the Galactic center and a number of GeV gamma-ray features, including the eastern edge of the recently-discovered northern Fermi Bubble. The polarized microwave features also correspond to structures seen in the all-sky 408 MHz total intensity data, including the Galactic center spur. The magnetic field structure revealed by the polarization data at 23 GHz suggests that neither the emission coincident with the Bubble edge nor the Galactic center spur are likely to be features of the local ISM. On the basis of the observed morphological correspondences, similar inferred spectra, and the similar energetics of all sources, we suggest a direct connection between the Galactic center spur and the northern Fermi Bubble. △ Less

Submitted 2 February, 2012; v1 submitted 21 January, 2012; originally announced January 2012.

Comments: Accepted for publication in The Astrophysical Journal Letters after minor changes

Journal ref: 2012 ApJ 747 L12

The Galactic Centre - A Laboratory for Starburst Galaxies (?)

Authors: Roland M. Crocker

Abstract: The Galactic centre - as the closest galactic nucleus - holds both intrinsic interest and possibly represents a useful analogue to star-burst nuclei which we can observe with orders of magnitude finer detail than these external systems. The environmental conditions in the GC - here taken to mean the inner 200 pc in diameter of the Milky Way - are extreme with respect to those typically encountered… ▽ More The Galactic centre - as the closest galactic nucleus - holds both intrinsic interest and possibly represents a useful analogue to star-burst nuclei which we can observe with orders of magnitude finer detail than these external systems. The environmental conditions in the GC - here taken to mean the inner 200 pc in diameter of the Milky Way - are extreme with respect to those typically encountered in the Galactic disk. The energy densities of the various GC ISM components are typically ~two orders of magnitude larger than those found locally and the star-formation rate density ~three orders of magnitude larger. Unusually within the Galaxy, the Galactic centre exhibits hard-spectrum, diffuse TeV (=10^12 eV) gamma-ray emission spatially coincident with the region's molecular gas. Recently the nuclei of local star-burst galaxies NGC 253 and M82 have also been detected in gamma-rays of such energies. We have embarked on an extended campaign of modelling the broadband (radio continuum to TeV gamma-ray), non- thermal signals received from the inner 200 pc of the Galaxy. On the basis of this modelling we find that star-formation and associated supernova activity is the ultimate driver of the region's non-thermal activity. This activity drives a large-scale wind of hot plasma and cosmic rays out of the GC. The wind advects the locally-accelerated cosmic rays quickly, before they can lose much energy in situ or penetrate into the densest molecular gas cores where star-formation occurs. The cosmic rays can, however, heat/ionize the lower density/warm H2 phase enveloping the cores. On very large scales (~10 kpc) the non-thermal signature of the escaping GC cosmic rays has probably been detected recently as the spectacular 'Fermi bubbles' and corresponding 'WMAP haze'. △ Less

Submitted 29 December, 2011; originally announced December 2011.

Comments: Invited talk to appear in Proceedings of IAU Symposium No. 284, 2011 (R.J. Tuffs & C.C. Popescu, eds.) `The Spectral Energy Distribution of Galaxies'