The First Stars: formation under cosmic ray feedback


Jacob 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


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




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 i