ROUGH DRAFT authorea.com/7907

# -1cmJacob A. Hummel$$^{1}$$, Athena Stacy$$^{2}$$ and Volker Bromm$$^{1}$$ $$^1$$Department of Astronomy, The University of Texas at Austin, TX 78712, USA $$^2$$Department of Astronomy, University of California, Berkeley, CA 94720, USA

-1cm

stars: formation — stars: Population III — cosmology: theory — early Universe — dark ages, first stars

-1cm

# Introduction

\label{intro}

Cosmic rays (CRs) have long been known to play an important role in the complex physics and chemistry of the interstellar medium (ISM) in local galaxies (Spitzer et al. 1968; Spitzer et al. 1969; Glassgold et al. 1973; Goldsmith et al. 1978; Cravens et al. 1978; Mannheim et al. 1994; Tielens 2005; recently reviewed by Strong et al. 2007; Grenier et al. 2015). They are an effective source of ionisation and heating in various environments, from preheating the primordial intergalactic medium (IGM; Sazonov et al., 2015), to driving outflows in the diffuse ISM (e.g., Enßlin et al., 2007; Jubelgas et al., 2008; Salem et al., 2014; Hanasz et al., 2013; Booth et al., 2013; Salem et al., 2014a), to providing an important (often dominant) source of heating and ionisation in deeply embedded gas clouds and protostellar discs (Spitzer, 1978; Dalgarno et al., 1999; Indriolo et al., 2009; Padovani et al., 2009; Glassgold et al., 2012; Padovani et al., 2013; Padovani et al., 2015).

CRs are particularly interesting in the context of Population III (Pop III) star formation, as they provide a continually replenished source of free electrons, enhancing the formation of molecular hydrogen—the only coolant available in primordial gas (Lepp et al., 1984; Abel et al., 1997; Galli et al., 1998; Bromm et al., 2002; Glover et al., 2008). This enhances the ability of the gas to cool, modifying the characteristic density and temperature at which runaway gravitational collapse sets in, and possibly the characteristic mass of the stars thus formed. The characteristic mass of Pop III is critical, as it largely controls the extent to which the first stars influence their environment, determining both their total luminosity and ionising radiation output (Schaerer, 2002), in addition to the details of their eventual demise (Heger et al., 2003; Heger et al., 2010; Maeder et al., 2012). As such, a thorough understanding of how the very first stars impact subsequent episodes of metal-free star formation—sometimes referred to as Pop III.1 and Pop III.2, respectively (McKee et al., 2008)—is crucial to developing a comprehensive picture of cosmic evolution.

In the absence of any feedback, pioneering numerical studies suggested that the very first stars were quite massive—on the order of $$100{\,{\rm M}_{\odot}}$$ (e.g., Bromm et al., 1999; Bromm et al., 2002; Abel et al., 2002; Yoshida et al., 2003; Bromm et al., 2004; Yoshida et al., 2006; O’Shea et al., 2007). However, more recent simulations, aided by increased resolution, have found that significant fragmentation occurs during the star formation process (Turk et al., 2009; Stacy et al., 2010; Clark et al., 2011; Clark et al., 2011a; Greif et al., 2011; Greif et al., 2012; Stacy et al., 2013; Hirano et al., 2014; Hosokawa et al., 2015). While the precise functional form of the Pop III initial mass function (IMF) remains uncertain, the emerging consensus is that these stars were somewhat more massive than stars formed at the present epoch, with a characteristic mass of $$\sim$$ a few $$\times 10{\,{\rm M}_{\odot}}$$ (Bromm, 2013; Greif, 2015). Several studies have investigated the impact these stars have on chemical enrichment (Madau et al. 2001; Mori et al. 2002; Bromm et al. 2003; Mackey et al. 2003; Heger et al. 2003; Umeda et al. 2003; Bromm et al. 2004; Kitayama et al. 2005; Tornatore et al. 2007; Greif et al. 2007; Greif et al. 2010; Wise et al. 2008; Maio et al. 2011; recently reviewed in Whalen et al. 2008; Joggerst et al. 2010; Karlsson et al. 2013) and reionisation (Kitayama et al., 2004; Sokasian et al., 2004; Whalen et al., 2004; Alvarez et al., 2006; Johnson et al., 2007; Robertson et al., 2010), as well as the effects of a Pop III-generated UV background on star formation in neighboring minihalos (Shapiro et al., 2004; Iliev et al., 2005; O’Shea et al., 2005; Susa et al., 2006; Susa, 2007; Whalen et al., 2008; Hasegawa et al., 2009; Whalen et al., 2010). The impact of high-energy radiation on subsequent Pop III star formation however has been less thoroughly explored.

As Pop III stars form in a predominantly neutral medium, the majority of their ionising output is absorbed, allowing only radiation less energetic than the Lyman-$$\alpha$$ transition to escape the immediate vicinity of the star-forming halo. While far-ultraviolet radiation in the Lyman-Werner (LW) bands ($$11.2 - 13.6{\,{\rm eV}}$$) lacks sufficient energy to interact with atomic hydrogen, it can still effectively dissociate molecular hydrogen. Similarly, photons capable of dissociating $$H^-$$ limit the ability of the gas to cool by eliminating the $$H^-$$ channel for producing $${\mathrm{H}_2}$$ (Agarwal et al., 2012; Agarwal et al., 2016). However, studies have found that the expected mean value of such LW radiation is far below the critical flux required to suppress $${\mathrm{H}_2}$$ cooling ; moreover, not all LW photons produced escape their host halo to contribute to this background (Schauer et al., 2015). At the high energy end, the neutral hydrogen cross section for X-rays and cosmic rays is small, allowing them to easily escape their host minihaloes. We recently investigated the impact of a cosmic X-ray background generated by high-mass X-ray binaries on primordial star formation (Hummel et al., 2015); here we focus on the impact of a CR background consisting of particles accelerated in supernova shock waves via the first-order Fermi process (see Section \ref{sec:context}).

While the uncertainties involved in estimating the strength of the high-$$z$$ CR background are quite large, measurements of the $$^6{\rm Li}$$ abundance from metal-poor stars in the Galactic halo provide a useful constraint: the observed abundance of $$^6{\rm Li}$$ is approximately 1000 times higher than predicted by big bang nucleosynthesis (Asplund et al., 2006), strong evidence for the existence of a CR spallation channel to provide a $$^6{\rm Li}$$ bedrock abundance prior to the bulk of star formation (Rollinde et al., 2005; Rollinde et al., 2006). Production of this first pervasive CR background by shock-acceleration in Pop III supernovae dovetails nicely with the upper limits this places on the CR energy density at high redshifts (Rollinde et al., 2006).

Early studies of the impact of CRs on primordial star formation focused on the production of ultra-high-energy CRs (UHECRs) by the decay of ultra-heavy X particles (Shchekinov et al., 2004; Vasiliev et al., 2006; Ripamonti et al., 2007). With energies above the Greisen–Zatsepin–Kuzmin (GZK) cutoff (Greisen, 1966; Zatsepin et al., 1966), UHECRs interact with the cosmic microwave background (CMB) to produce ionising photons, which in turn enhance the free electron fraction of the gas. Other work investigated the direct collisional ionisation of neutral hydrogen by SN shock-generated CRs, using one-zone models to determine the impact of a CR background on the chemical and thermal evolution of the gas in a minihalo (Stacy et al., 2007; Jasche et al., 2007). These studies found that the presence of a CR background enhanced molecular hydrogen formation in the minihalo, cooling the gas and lowering the Jeans mass, and by extension, the characteristic mass of the stars formed. We expand upon these one-zone models, using three-dimensional ab initio cosmological simulations to investigate the impact of a CR background on Pop III stars forming in a minihalo.

This paper is organized as follows: In Section \ref{sec:context} we provide the cosmological context for this study, estimating the expected intensity of the CR background. Our numerical methodology is described in Section \ref{sec:methods}, while our results are presented in Section \ref{sec:results}. Finally, our conclusions are gathered in Section \ref{conclusions}. Throughout this paper we adopt a $$\Lambda$$CDM model of hierarchical structure formation, using the following cosmological parameters, consistent with the latest measurements from the Planck Collaboration (Ade et al., 2015): $$\Omega_{\Lambda} = 0.7$$, $$\Omega_{\rm m} = 0.3$$, $$\Omega_{\rm B} = 0.04$$, and $$H_0 = 70 {\,\mathrm{km}\,\mathrm{s}^{-1}}{\,{\rm Mpc}}^{-1}$$.