Cosmic-ray Energy Density Research Articles

A source of gamma rays coincident with the shell of the supernova remnant CTB 80

ABSTRACT CTB 80 (G69.0+2.7) is a relatively old (50–80 kyr) supernova remnant (SNR) with a complex radio morphology showing three extended radio arms and a radio and X-ray nebula near the location of the pulsar PSR B1951+32. We report on a study of the GeV emission in the region of CTB 80 with Fermi-Large Area Telescope data. An extended source with a size of 1.3°, matching the size of the infrared shell associated to the SNR, was discovered. The GeV emission, detected up to an energy of ∼20 GeV, is more significant at the location of the northern radio arm where previous observations imply that the SNR shock is interacting with ambient material. Both hadronic and leptonic scenarios can reproduce the multiwavelength data reasonably well. The hadronic cosmic ray energy density required is considerably larger than the local Galactic value and the gamma-ray leptonic emission is mainly due to bremsstrahlung interactions. We conclude that GeV particles are still trapped or accelerated by the SNR producing the observed high-energy emission when interacting with ambient material.

Testing physical models for cosmic ray transport coefficients on galactic scales: self-confinement and extrinsic turbulence at ∼GeV energies

ABSTRACT The microphysics of ∼ GeV cosmic ray (CR) transport on galactic scales remain deeply uncertain, with almost all studies adopting simple prescriptions (e.g. constant diffusivity). We explore different physically motivated, anisotropic, dynamical CR transport scalings in high-resolution cosmological Feedback In Realistic Environment (FIRE) simulations of dwarf and ∼L* galaxies where scattering rates vary with local plasma properties motivated by extrinsic turbulence (ET) or self-confinement (SC) scenarios, with varying assumptions about e.g. turbulent power spectra on un-resolved scales, Alfvén-wave damping, etc. We self-consistently predict observables including γ-rays (Lγ), grammage, residence times, and CR energy densities to constrain the models. We demonstrate many non-linear dynamical effects (not captured in simpler models) tend to enhance confinement. For example, in multiphase media, even allowing arbitrary fast transport in neutral gas does not substantially reduce CR residence times (or Lγ), as transport is rate-limited by the ionized WIM and ‘inner CGM’ gaseous halo (104–106 K gas within $\lesssim 10\!-\!30\,$ kpc), and Lγ can be dominated by trapping in small ‘patches’. Most physical ET models contribute negligible scattering of ∼1–10 GeV CRs, but it is crucial to account for anisotropy and damping (especially of fast modes) or else scattering rates would violate observations. We show that the most widely assumed scalings for SC models produce excessive confinement by factors ≳100 in the warm ionized medium (WIM) and inner CGM, where turbulent and Landau damping dominate. This suggests either a breakdown of quasi-linear theory used to derive the CR transport parameters in SC, or that other novel damping mechanisms dominate in intermediate-density ionized gas.

Evidence of Cosmic-Ray Excess from Local Giant Molecular Clouds

Abstract We report the analysis of the Fermi-Large Area Telescope data from six nearby giant molecular clouds (MCs) belonging to the Gould Belt and the Aquila Rift regions. The high statistical γ-ray spectra above 3 GeV well described by power laws make it possible to derive precise estimates of the cosmic-ray (CR) distribution in the MCs. The comparison of γ-ray spectra of Taurus, Orion A, and Orion B clouds with the model expected from Alpha Magnetic Spectrometer (AMS-02) CR measurements confirms these clouds as passive clouds, immersed in an AMS-02-like CR spectrum. A similar comparison of Aquila Rift, Rho Oph, and Cepheus spectra yields significant deviation in both spectral indices and absolute fluxes, which can imply an additional acceleration of CRs throughout the entire clouds. Besides, the theoretical modeling of the excess γ-ray spectrum of these clouds, assuming π 0-decay interaction of CRs in the cloud, gives a considerable amount of an enhanced CR energy density and it shows a significant deviation in spectral shapes compared to the average AMS-02 CR spectrum between 30 GeV and 10 TeV. We suggest that this variation in the CR spectrum of Cepheus could be accounted for by an efficient acceleration in the shocks of winds of OB associations, while in Rho Oph, similar acceleration can be provided by multiple T-Tauri stars populated in the whole cloud. In the case of Aquila Rift, the excess in absolute CR flux could be related to an additional acceleration of CRs by supernova remnants or propagation effects in the cloud.

Constraints on the Emission of Gamma-Rays from M31 with HAWC

Abstract Cosmic rays, along with stellar radiation and magnetic fields, are known to make up a significant fraction of the energy density of galaxies such as the Milky Way. When cosmic rays interact in the interstellar medium, they produce gamma-ray emission which provides an important indication of how the cosmic rays propagate. Gamma-rays from the Andromeda galaxy (M31), located 785 kpc away, provide a unique opportunity to study cosmic-ray acceleration and diffusion in a galaxy with a structure and evolution very similar to the Milky Way. Using 33 months of data from the High Altitude Water Cherenkov Observatory, we search for teraelectronvolt gamma-rays from the galactic plane of M31. We also investigate past and present evidence of galactic activity in M31 by searching for Fermi bubble-like structures above and below the galactic nucleus. No significant gamma-ray emission is observed, so we use the null result to compute upper limits on the energy density of cosmic rays >10 TeV in M31.

Radio haloes of star-forming galaxies

ABSTRACT We study the synchrotron radio emission from extra-planar regions of star-forming galaxies. We use ideal magnetohydrodynamic simulations of a rotating Milky Way-type disc galaxy with distributed star formation sites for three star formation rates (0.3, 3, 30 M⊙ yr−1). From our simulations, we see emergence of galactic-scale magnetized outflows, carrying gas from the disc. We compare the morphology of the outflowing gas with hydrodynamic simulations. We look at the spatial distribution of magnetic field in the outflows. Assuming that a certain fraction of gas energy density is converted into cosmic ray energy density, and using information about the magnetic field, we obtain synchrotron emissivity throughout the simulation domain. We generate the surface brightness maps at 1.4 GHz. The outflows are more extended in the vertical direction than radial and hence have an oblate shape. We further find that the matter right behind the outer shock shines brighter in these maps than that above or below. To understand whether this feature can be observed, we produce vertical intensity profiles. We convolve the vertical intensity profile with the typical beam sizes of radio telescopes, for a galaxy located at 10 Mpc to estimate the radio scale height and compare with observations. The radio scale height is ∼300–1200 pc, depending on the resolution of the telescope. We relate the advection speed of the outer shock with the surface density of star formation as ${\rm v}_{\rm adv} \propto \Sigma _{\rm SFR}^{0.3}$, which is consistent with earlier observations and analytical estimates.

CHANG-ES

Context. Cosmic-ray electrons (CREs) originating from the star-forming discs of spiral galaxies frequently form extended radio haloes that are best observable in edge-on galaxies, where their properties can be directly investigated as a function of vertical height above the disc.Aims. For the present study, we selected two nearby edge-on galaxies from the Continuum Halos in Nearby Galaxies – an EVLA Survey (CHANG-ES), NGC 891 and 4565, which differ largely in their detectable halo extent and their star-formation rates (SFRs). Our aim is to figure out how such differences are related to the (advective and/or diffusive) CRE transport in the disc and in the halo.Methods. We use wide-band 1.5 and 6 GHz Very Large Array (VLA) observations obtained in the B, C, and D configurations, and combine the 6 GHz images with Effelsberg observations to correct for missing short spacings. After subtraction of the thermal emission, we investigate the spatially resolved synchrotron spectral index distribution in terms of CRE spectral ageing. We further compute total magnetic field strengths assuming equipartition between the cosmic-ray (CR) energy density and the magnetic field, and measure synchrotron scale heights at both frequencies. Based on the fitted vertical profiles of the synchrotron intensity and on the spectral index profile between 1.5 and 6 GHz, we create purely advective and purely diffusive CRE transport models by numerically solving the 1D diffusion–loss equation. In particular, we investigate for the first time the radial dependence of synchrotron and magnetic field scale heights, advection speeds, and diffusion coefficients, whereas previous studies of these two galaxies only determined global values of these quantities.Results. We find that the overall spectral index distribution of NGC 891 is mostly consistent with continuous CRE injection. In NGC 4565, many of the local synchrotron spectra (even in the disc) feature a break between 1.5 and 6 GHz and are thus more in line with discrete-epoch CRE injection (Jaffe–Perola (JP) or Kardashev–Pacholczyk (KP) models). This implies that CRE injection time-scales are lower than the synchrotron cooling time-scales. The synchrotron scale height of NGC 891 increases with radius, indicating that synchrotron losses are significant. NGC 891 is probably dominated by advective CRE transport at a velocity of ≳150 km s−1. In contrast, NGC 4565 is diffusion-dominated up toz = 1 kpc or higher, with a diffusion coefficient of ≥2 × 1028 cm2 s−1.

H3+, the ideal probe for in situ measurement of the Galactic cosmic rays.

The reason why [Formula: see text] is an ideal probe for in situ measurement of cosmic ray energy density in the Galaxy is discussed. The variation of measured values of cosmic ray ionization rate ζ of H2 since the theoretical prediction by Spitzer & Tomasko (Spitzer & Tomasko 1968 Astrophys. J. 152, 971-986) is reviewed. How the measurements by [Formula: see text] established ζ as of the order of 10-17 s-1 for dense clouds, 10-16 s-1 for diffuse clouds in the Galactic disc, and 10-14 s-1 for warm and diffuse gas in the Galactic centre is discussed. This article is part of a discussion meeting issue 'Advances in hydrogen molecular ions: H3+, H5+ and beyond'.

Fickian and non-Fickian diffusion of cosmic rays

Abstract Fluid approximations to cosmic ray (CR) transport are often preferred to kinetic descriptions in studies of the dynamics of the interstellar medium (ISM) of galaxies, because they allow simpler analytical and numerical treatments. Magnetohydrodynamic simulations of the ISM usually incorporate CR dynamics as an advection–diffusion equation for CR energy density, with anisotropic, magnetic-field-aligned diffusion with the diffusive flux assumed to obey Fick’s law. We compare test particle and fluid simulations of CRs in a random magnetic field. We demonstrate that a non-Fickian prescription of CR diffusion, which corresponds to the telegraph equation for the CR energy density, can be easily calibrated to match the test particle simulations with great accuracy. In particular, we consider a random magnetic field in the fluid simulation that has a lower spatial resolution than that used in the particle simulation to demonstrate that an appropriate choice of the diffusion tensor can account effectively for the unresolved (subgrid) scales of the magnetic field. We show that the characteristic time that appears in the telegraph equation can be physically interpreted as the time required for the particles to reach a diffusive regime and we stress that the Fickian description of the CR fluid is unable to describe complex boundary or initial conditions for the CR energy flux.

Cosmic-ray propagation in the bi-stable interstellar medium

Context. Cosmic rays propagate through the galactic scales down to the smaller scales at which stars form. Cosmic rays are close to energy equipartition with the other components of the interstellar medium and can provide a support against gravity if pressure gradients develop. Aims. We study the propagation of cosmic rays within the turbulent and magnetised bi-stable interstellar gas. The conditions necessary for cosmic-ray trapping and cosmic-ray pressure gradient development are investigated. Methods. We derived an analytical value of the critical diffusion coefficient for cosmic-ray trapping within a turbulent medium, which follows the observed scaling relations. We then presented a numerical study using 3D simulations of the evolution of a mixture of interstellar gas and cosmic rays, in which turbulence is driven at varying scales by stochastic forcing within a box of 40 pc. We explored a large parameter space in which the cosmic-ray diffusion coefficient, the magnetisation, the driving scale, and the amplitude of the turbulence forcing, as well as the initial cosmic-ray energy density, vary. Results. We identify a clear transition in the interstellar dynamics for cosmic-ray diffusion coefficients below a critical value deduced from observed scaling relations. This critical diffusion depends on the characteristic length scale L of Dcrit ≃ 3.1 × 1023 cm2 s−1(L/1 pc)q+1, where the exponent q relates the turbulent velocity dispersion σ to the length scale as σ ~ Lq. Hence, in our simulations this transition occurs around Dcrit ≃ 1024–1025 cm2 s−1. The transition is recovered in all cases of our parameter study and is in very good agreement with our simple analytical estimate. In the trapped cosmic-ray regime, the induced cosmic-ray pressure gradients can modify the gas flow and provide a support against the thermal instability development. We discuss possible mechanisms that can significantly reduce the cosmic-ray diffusion coefficients within the interstellar medium. Conclusions. Cosmic-ray pressure gradients can develop and modify the evolution of thermally bi-stable gas for diffusion coefficients D0 ≤ 1025 cm2 s−1 or in regions where the cosmic-ray pressure exceeds the thermal one by more than a factor of ten. This study provides the basis for further works including more realistic cosmic-ray diffusion coefficients, as well as local cosmic-ray sources.

ATCA 16 cm observation of CIZA J1358.9−4750: Implication of merger stage and constraint on non-thermal properties

Abstract We report the Australia Telescope Compact Array 16 cm observation of CIZA J1358.9−4750. Recent X-ray studies imply that this galaxy cluster is composed of merging, binary clusters. Using the EW367 configuration, we found no significant diffuse radio emission in and around the cluster. An upper limit of the total radio power at 1.4 GHz is ∼1.1 × 1022 W Hz−1 in 30 square arcminutes, which is a typical size for radio relics. It is known that an empirical relation holds between the total radio power and X-ray luminosity of the host cluster. The upper limit is about one order of magnitude lower than the power expected from the relation. Very young (∼70 Myr) shocks with low Mach numbers (∼1.3), which are often seen at an early stage of merger simulations, are suggested by the previous X-ray observation. The shocks may generate cosmic-ray electrons with a steep energy spectrum, which is consistent with non-detection of bright (>1023 W Hz−1) relic in this 16 cm band observation. Based on the assumption of energy equipartition, the upper limit gives a magnetic field strength of below 0.68f(Dlos/1 Mpc)−1(γmin/200)−1 μG, where f is the cosmic-ray total energy density over the cosmic-ray electron energy density, Dlos is the depth of the shock wave along the sightline, and γmin is the lower cutoff Lorentz factor of the cosmic-ray electron energy spectrum.

On the Foundation of Equipartition in Supernova Remnants

Abstract A widely accepted paradigm is that equipartition (eqp) between the energy density of cosmic rays (CRs) and the energy density of the magnetic field cannot be sustained in supernova remnants (SNRs). However, our 3D hydrodynamic supercomputer simulations, coupled with a nonlinear diffusive shock acceleration model, provide evidence that eqp may be established at the end of the Sedov phase of evolution in which most SNRs spend the longest portions of their lives. We introduce the term “constant partition” for any constant ratio between the CR energy density and the energy density of the magnetic field in an SNR, while the term “equipartition” should be reserved for the case of approximately the same values of the energy density (also, it is constant partition in the order of magnitude) of ultra-relativistic electrons only (or CRs in total) and the energy density of the magnetic field. Our simulations suggest that this approximate constant partition exists in all but the youngest SNRs. We speculate that since evolved SNRs at the end of the Sedov phase of evolution can reach eqp between CRs and magnetic fields, they may be responsible for initializing this type of eqp in the interstellar medium. Additionally, we show that eqp between the electron component of CRs and the magnetic field may be used for calculating the magnetic field strength directly from observations of synchrotron emission from SNRs. The values of magnetic field strengths in SNRs given here are approximately 2.5 times lower than values calculated by Arbutina et al.

Efficient cold outflows driven by cosmic rays in high-redshift galaxies and their global effects on the IGM

We present semi-analytical models of galactic outflows in high redshift galaxies driven by both hot thermal gas and non-thermal cosmic rays. Thermal pressure alone may not sustain a large scale outflow in low mass galaxies (i.e $M\sim 10^8$~M$_\odot$), in the presence of supernovae (SNe) feedback with large mass loading. We show that inclusion of cosmic ray pressure allows outflow solutions even in these galaxies. In massive galaxies for the same energy efficiency, cosmic ray driven winds can propagate to larger distances compared to pure thermally driven winds. On an average gas in the cosmic ray driven winds has a lower temperature which could aid detecting it through absorption lines in the spectra of background sources. Using our constrained semi-analytical models of galaxy formation (that explains the observed UV luminosity functions of galaxies) we study the influence of cosmic ray driven winds on the properties of the intergalactic medium (IGM) at different redshifts. In particular, we study the volume filling factor, average metallicity, cosmic ray and magnetic field energy densities for models invoking atomic cooled and molecular cooled halos. We show that the cosmic rays in the IGM could have enough energy that can be transferred to the thermal gas in presence of magnetic fields to influence the thermal history of the intergalactic medium. The significant volume filling and resulting strength of IGM magnetic fields can also account for recent $\gamma$-ray observations of blazars.

Magnetic field formation in the Milky Way like disc galaxies of the Auriga project

AbstractThe magnetic fields observed in the Milky Way and nearby galaxies appear to be in equipartition with the turbulent, thermal and cosmic ray energy densities, and hence are expected to be dynamically important. However, the origin of these strong magnetic fields is still unclear, and most previous attempts to simulate galaxy formation from cosmological initial conditions have ignored them altogether. Here, we analyse the magnetic fields predicted by the simulations of the Auriga Project, a set of 30 high-resolution cosmological zoom simulations of Milky Way like galaxies, carried out with a moving-mesh magnetohydrodynamics code and a detailed galaxy formation physics model. We find that the magnetic fields grow exponentially at early times owing to a small-scale dynamo with an e-folding time of roughly 100 Myr in the centre of haloes until saturation occurs around z = 2–3, when the magnetic energy density reaches about 10 per cent of the turbulent energy density with a typical strength of $10\text{--}50\,\rm {\mu G}$. In the galactic centres, the ratio between magnetic and turbulent energies remains nearly constant until z = 0. At larger radii, differential rotation in the discs leads to linear amplification that typically saturates around z = 0.5–0. The final radial and vertical variations of the magnetic field strength can be well described by two joint exponential profiles, and are in good agreement with observational constraints. Overall, the magnetic fields have only little effect on the global evolution of the galaxies as it takes too long to reach equipartition. We also demonstrate that our results are well converged with numerical resolution.

Low-frequency radio constraints on the synchrotron cosmic web

We present a search for the synchrotron emission from the synchrotron cosmic web by cross correlating 180MHz radio images from the Murchison Widefield Array with tracers of large scale structure (LSS). We use two versions of the radio image covering $21.76\times 21.76$ degrees with point sources brighter than 0.05 Jy subtracted, with and without filtering of Galactic emission. As tracers of the LSS we use the Two-Micron-All-Sky-Survey (2MASS) and the Widefield InfraRed Explorer (WISE) redshift catalogues to produce galaxy number density maps. The cross correlation functions all show peak amplitudes at zero degrees, decreasing with varying slopes towards zero correlation over a range of one degree. The cross correlation signals include components from point source, Galactic, and extragalactic diffuse emission. We use models of the diffuse emission from smoothing the density maps with Gaussians of sizes 1-4 Mpc to find limits on the cosmic web components. From these models we find surface brightness 99.7 per cent upper limits in the range of 0.09-2.20 mJy beam$^{-1}$ (average beam size of 2.6 arcmin), corresponding to 0.01-0.30 mJy arcmin$^{-2}$. Assuming equipartition between energy densities of cosmic rays and the magnetic field, the flux density limits translate to magnetic field strength limits of 0.03-1.98 $\mu$G, depending heavily on the spectral index. We conclude that for a 3$\sigma$ detection of 0.1 $\mu$G magnetic field strengths via cross correlations, image depths of sub-mJy to sub-$\mu$Jy are necessary. We include discussion on the treatment and effect of extragalactic point sources and Galactic emission, and next steps for building on this work.

Intensity of galactic cosmic rays in the early sun epoch

The process of heliospheric modulation of intensity of galactic cosmic rays is investigated by solving the transport equation. The spatial-energetic distribution of cosmic rays in the present epoch and in the past is analyzed. It is demonstrated that the particle density and the energy density of cosmic rays in the Solar System in the distant past were much lower than the corresponding current values. The cosmic ray intensity modulation in the early heliosphere was especially strong in the case of low-energy particles.

MODELING THE GAMMA-RAY EMISSION IN THE GALACTIC CENTER WITH A FADING COSMIC-RAY ACCELERATOR

ABSTRACT Recent HESS observations of the ∼200 pc scale diffuse gamma-ray emission from the central molecular zone (CMZ) suggest the presence of a PeV cosmic-ray accelerator (PeVatron) located in the inner 10 pc region of the Galactic center. Interestingly, the gamma-ray spectrum of the point-like source (HESS J1745-290) in the Galactic center shows a cutoff at ∼10 TeV, implying a cutoff around 100 TeV in the cosmic-ray proton spectrum. Here we propose that the gamma-ray emission from the inner and the outer regions may be explained self-consistently by run-away protons from a single yet fading accelerator. In this model, gamma-rays from the CMZ region are produced by protons injected in the past, while gamma-rays from the inner region are produced by protons injected more recently. We suggest that the blast wave formed in a tidal disruption event (TDE) caused by the supermassive black hole (Sgr A*) could serve as such a fading accelerator. With typical parameters of the TDE blast wave, gamma-ray spectra of both the CMZ region and HESS J1745-290 can be reproduced simultaneously. Meanwhile, we find that the cosmic-ray energy density profile in the CMZ region may also be reproduced in the fading accelerator model when appropriate combinations of the particle injection history and the diffusion coefficient of cosmic rays are adopted.

Equipartition and cosmic ray energy densities in central molecular zones of starbursts

Abstract The energy densities in magnetic fields and cosmic rays (CRs) in galaxies are often assumed to be in equipartition, allowing for an indirect estimate of the magnetic field strength from the observed radio synchrotron spectrum. However, both primary and secondary CRs contribute to the synchrotron spectrum, and the CR electrons also loose energy via bremsstrahlung and inverse Compton. While classical equipartition formulae avoid these intricacies, there have been recent revisions that account for the extreme conditions in starbursts. Yet, the application of the equipartition formula to starburst environments also presupposes that time-scales are long enough to reach equilibrium. Here, we test equipartition in the central molecular zones (CMZs) of nearby starburst galaxies by modelling the observed γ-ray spectra, which provide a direct measure of the CR energy density, and the radio spectra, which provide a probe of the magnetic field strength. We find that in starbursts, the magnetic field energy density is significantly larger than the CR energy density, demonstrating that the equipartition argument is frequently invalid for CMZs.

Turbulence in the diffuse magneto-ionized medium: observational aspects

AbstractTurbulence in the interstellar medium is ubiquitous. The turbulent energy density in the gas is significant, and comparable to energy densities of magnetic fields and cosmic rays. Studies of the turbulent interstellar gas in the Milky Way have mostly focused on the neutral gas component, since various spectral lines can give velocity information. Probing turbulent properties in the ionized gas, let alone in magnetic fields, is observationally more difficult. A number of observational methods are discussed below which provide estimates of the maximum scale of fluctuations, the Mach number and other turbulence characteristics.

EFFECTIVE DESTRUCTION OF CO BY COSMIC RAYS: IMPLICATIONS FOR TRACING H2GAS IN THE UNIVERSE

We report on the effects of cosmic rays (CRs) on the abundance of CO in $\rm H_2$ clouds under conditions typical for star-forming galaxies in the Universe. We discover that this most important molecule for tracing H$_2$ gas is very effectively destroyed in ISM environments with CR energy densities $\rm U_{CR}\sim(50-10^{3})\times U_{CR,Gal}$, a range expected in numerous star-forming systems throughout the Universe. This density-dependent effect operates volumetrically rather than only on molecular cloud surfaces (i.e. unlike FUV radiation that also destroys CO), and is facilitated by: a) the direct destruction of CO by CRs, and b) a reaction channel activated by CR-produced He$^{+}$. The effect we uncover is strong enough to render Milky-Way type Giant Molecular Clouds (GMCs) very CO-poor (and thus CO-untraceable), even in ISM environments with rather modestly enhanced average CR energy densities of $\rm U_{CR}\sim(10-50)\times\rm U_{CR,Gal}$. We conclude that the CR-induced destruction of CO in molecular clouds, unhindered by dust absorption, is perhaps the single most important factor controlling the CO-visibility of molecular gas in vigorously star-forming galaxies. We anticipate that a second order effect of this CO destruction mechanism will be to make the H$_2$ distribution in the gas-rich disks of such galaxies appear much clumpier in CO $J$=1--0, 2--1 line emission than it actually is. Finally we give an analytical approximation of the CO/H$_2$ abundance ratio as a function of gas density and CR energy density for use in galaxy-size or cosmological hydrodynamical simulations, and propose some key observational tests.

Erich Regener and the ionisation maximum of the atmosphere

Abstract. In the 1930s the German physicist Erich Regener (1881–1955) did important work on the measurement of the rate of production of ionisation deep under water and in the atmosphere. Along with one of his students, Georg Pfotzer, he discovered the altitude at which the production of ionisation in the atmosphere reaches a maximum, often, but misleadingly, called the Pfotzer maximum. Regener was one of the first to estimate the energy density of cosmic rays, an estimate that was used by Baade and Zwicky to bolster their postulate that supernovae might be their source. Yet Regener's name is less recognised by present-day cosmic ray physicists than it should be, largely because in 1937 he was forced to take early retirement by the National Socialists as his wife had Jewish ancestors. In this paper we briefly review his work on cosmic rays and recommend an alternative naming of the ionisation maximum. The influence that Regener had on the field through his son, his son-in-law, his grandsons and his students, and through his links with Rutherford's group in Cambridge, is discussed in an appendix. Regener was nominated for the Nobel Prize in Physics by Schrödinger in 1938. He died in 1955 at the age of 73.