skip to main content


Title: Stability and pulsation of the first dark stars
ABSTRACT The first bright objects to form in the Universe might not have been ‘ordinary’ fusion-powered stars, but ‘dark stars’ (DSs) powered by the annihilation of dark matter (DM) in the form of weakly interacting massive particles (WIMPs). If discovered, DSs can provide a unique laboratory to test DM models. DSs are born with a mass of the order of M⊙ and may grow to a few million solar masses; in this work we investigate the properties of early DSs with masses up to $\sim \! 1000 \, \mathrm{ M}_\odot$, fueled by WIMPS weighing 100 GeV. We improve the previous implementation of the DM energy source into the stellar evolution code mesa. We show that the growth of DSs is not limited by astrophysical effects: DSs up to $\sim \!1000 \, \mathrm{ M}_{\odot }$ exhibit no dynamical instabilities; DSs are not subject to mass-loss driven by super-Eddington winds. We test the assumption of previous work that the injected energy per WIMP annihilation is constant throughout the star; relaxing this assumption does not change the properties of the DSs. Furthermore, we study DS pulsations, for the first time investigating non-adiabatic pulsation modes, using the linear pulsation code gyre. We find that acoustic modes in DSs of masses smaller than $\sim \! 200 \, \mathrm{ M}_\odot$ are excited by the κ − γ and γ mechanism in layers where hydrogen or helium is (partially) ionized. Moreover, we show that the mass-loss rates potentially induced by pulsations are negligible compared to the accretion rates.  more » « less
Award ID(s):
1716436 1663696
NSF-PAR ID:
10287877
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
503
Issue:
3
ISSN:
0035-8711
Page Range / eLocation ID:
3677 to 3691
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT

    Ultra-massive white dwarf stars are currently being discovered at a considerable rate, thanks to surveys such as the Gaia space mission. These dense and compact stellar remnants likely play a major role in Type Ia supernova explosions. It is possible to probe the interiors of ultra-massive white dwarfs through asteroseismology. In the case of the most massive white dwarfs, general relativity could affect their structure and pulsations substantially. In this work, we present results of relativistic pulsation calculations employing relativistic ultra-massive ONe-core white dwarf models with hydrogen-rich atmospheres and masses ranging from 1.29 to $1.369 \ \mathrm{M}_{\odot }$ with the aim of assessing the impact of general relativity on the adiabatic gravity (g)-mode period spectrum of very high mass ZZ Ceti stars. Employing the relativistic Cowling approximation for the pulsation analysis, we find that the critical buoyancy (Brunt–Väisälä) and acoustic (Lamb) frequencies are larger for the relativistic case, compared to the Newtonian case, due to the relativistic white dwarf models having smaller radii and higher gravities for a fixed stellar mass. In addition, the g-mode periods are shorter in the relativistic case than those in the Newtonian computations, with relative differences of up to ∼$50$ per cent for the highest mass models ($1.369 \ \mathrm{M}_{\odot }$) and for effective temperatures typical of the ZZ Ceti instability strip. Hence, the effects of general relativity on the structure, evolution, and pulsations of white dwarfs with masses larger than ∼$1.29 \ \mathrm{M}_{\odot }$ cannot be ignored in the asteroseismological analysis of ultra-massive ZZ Ceti stars.

     
    more » « less
  2. ABSTRACT We present and study a large suite of high-resolution cosmological zoom-in simulations, using the FIRE-2 treatment of mechanical and radiative feedback from massive stars, together with explicit treatment of magnetic fields, anisotropic conduction and viscosity (accounting for saturation and limitation by plasma instabilities at high β), and cosmic rays (CRs) injected in supernovae shocks (including anisotropic diffusion, streaming, adiabatic, hadronic and Coulomb losses). We survey systems from ultrafaint dwarf ($M_{\ast }\sim 10^{4}\, \mathrm{M}_{\odot }$, $M_{\rm halo}\sim 10^{9}\, \mathrm{M}_{\odot }$) through Milky Way/Local Group (MW/LG) masses, systematically vary uncertain CR parameters (e.g. the diffusion coefficient κ and streaming velocity), and study a broad ensemble of galaxy properties [masses, star formation (SF) histories, mass profiles, phase structure, morphologies, etc.]. We confirm previous conclusions that magnetic fields, conduction, and viscosity on resolved ($\gtrsim 1\,$ pc) scales have only small effects on bulk galaxy properties. CRs have relatively weak effects on all galaxy properties studied in dwarfs ($M_{\ast } \ll 10^{10}\, \mathrm{M}_{\odot }$, $M_{\rm halo} \lesssim 10^{11}\, \mathrm{M}_{\odot }$), or at high redshifts (z ≳ 1–2), for any physically reasonable parameters. However, at higher masses ($M_{\rm halo} \gtrsim 10^{11}\, \mathrm{M}_{\odot }$) and z ≲ 1–2, CRs can suppress SF and stellar masses by factors ∼2–4, given reasonable injection efficiencies and relatively high effective diffusion coefficients $\kappa \gtrsim 3\times 10^{29}\, {\rm cm^{2}\, s^{-1}}$. At lower κ, CRs take too long to escape dense star-forming gas and lose their energy to collisional hadronic losses, producing negligible effects on galaxies and violating empirical constraints from spallation and γ-ray emission. At much higher κ CRs escape too efficiently to have appreciable effects even in the CGM. But around $\kappa \sim 3\times 10^{29}\, {\rm cm^{2}\, s^{-1}}$, CRs escape the galaxy and build up a CR-pressure-dominated halo which maintains approximate virial equilibrium and supports relatively dense, cool (T ≪ 106 K) gas that would otherwise rain on to the galaxy. CR ‘heating’ (from collisional and streaming losses) is never dominant. 
    more » « less
  3. ABSTRACT

    Previous studies of fueling black holes in galactic nuclei have argued (on scales ${\sim}0.01{-}1000\,$pc) accretion is dynamical with inflow rates $\dot{M}\sim \eta \, M_{\rm gas}/t_{\rm dyn}$ in terms of gas mass Mgas, dynamical time tdyn, and some η. But these models generally neglected expulsion of gas by stellar feedback, or considered extremely high densities where expulsion is inefficient. Studies of star formation, however, have shown on sub-kpc scales the expulsion efficiency fwind = Mejected/Mtotal scales with the gravitational acceleration as $(1-f_{\rm wind})/f_{\rm wind}\sim \bar{a}_{\rm grav}/\langle \dot{p}/m_{\ast }\rangle \sim \Sigma _{\rm eff}/\Sigma _{\rm crit}$ where $\bar{a}_{\rm grav}\equiv G\, M_{\rm tot}(\lt r)/r^{2}$ and $\langle \dot{p}/m_{\ast }\rangle$ is the momentum injection rate from young stars. Adopting this as the simplest correction for stellar feedback, $\eta \rightarrow \eta \, (1-f_{\rm wind})$, we show this provides a more accurate description of simulations with stellar feedback at low densities. This has immediate consequences, predicting the slope and normalization of the MBH − σ and MBH − Mbulge relation, LAGN −SFR relations, and explanations for outliers in compact Es. Most strikingly, because star formation simulations show expulsion is efficient (fwind ∼ 1) below total-mass surface density $M_{\rm tot}/\pi \, r^{2}\lt \Sigma _{\rm crit}\sim 3\times 10^{9}\, \mathrm{M}_{\odot }\, {\rm kpc^{-2}}$ (where $\Sigma _{\rm crit}=\langle \dot{p}/m_{\ast }\rangle /(\pi \, G)$), BH mass is predicted to specifically trace host galaxy properties above a critical surface brightness Σcrit (B-band $\mu _{\rm B}^{\rm crit}\sim 19\, {\rm mag\, arcsec^{-2}}$). This naturally explains why BH masses preferentially reflect bulge properties or central surface densities (e.g. $\Sigma _{1\, {\rm kpc}}$), not ‘total’ galaxy properties.

     
    more » « less
  4. null (Ed.)
    ABSTRACT Understanding the evolution of self-gravitating, isothermal, magnetized gas is crucial for star formation, as these physical processes have been postulated to set the initial mass function (IMF). We present a suite of isothermal magnetohydrodynamic (MHD) simulations using the gizmo code that follow the formation of individual stars in giant molecular clouds (GMCs), spanning a range of Mach numbers found in observed GMCs ($\mathcal {M} \sim 10\!-\!50$). As in past works, the mean and median stellar masses are sensitive to numerical resolution, because they are sensitive to low-mass stars that contribute a vanishing fraction of the overall stellar mass. The mass-weighted median stellar mass M50 becomes insensitive to resolution once turbulent fragmentation is well resolved. Without imposing Larson-like scaling laws, our simulations find $M_\mathrm{50} \,\, \buildrel\propto \over \sim \,\,M_\mathrm{0} \mathcal {M}^{-3} \alpha _\mathrm{turb}\, \mathrm{SFE}^{1/3}$ for GMC mass M0, sonic Mach number $\mathcal {M}$, virial parameter αturb, and star formation efficiency SFE = M⋆/M0. This fit agrees well with previous IMF results from the ramses, orion2, and sphng codes. Although M50 has no significant dependence on the magnetic field strength at the cloud scale, MHD is necessary to prevent a fragmentation cascade that results in non-convergent stellar masses. For initial conditions and SFE similar to star-forming GMCs in our Galaxy, we predict M50 to be $\gt 20 \, \mathrm{M}_{\odot }$, an order of magnitude larger than observed ($\sim 2 \, \mathrm{M}_\odot$), together with an excess of brown dwarfs. Moreover, M50 is sensitive to initial cloud properties and evolves strongly in time within a given cloud, predicting much larger IMF variations than are observationally allowed. We conclude that physics beyond MHD turbulence and gravity are necessary ingredients for the IMF. 
    more » « less
  5. ABSTRACT Self-interacting dark matter (SIDM) cosmologies admit an enormous diversity of dark matter (DM) halo density profiles, from low-density cores to high-density core-collapsed cusps. The possibility of the growth of high central density in low-mass haloes, accelerated if haloes are subhaloes of larger systems, has intriguing consequences for small-halo searches with substructure lensing. However, following the evolution of ${\lesssim}10^8 \, \mathrm{M}_\odot$ subhaloes in lens-mass systems (${\sim}10^{13}\, \mathrm{M}_\odot$) is computationally expensive with traditional N-body simulations. In this work, we develop a new hybrid semi-analytical + N-body method to study the evolution of SIDM subhaloes with high fidelity, from core formation to core-collapse, in staged simulations. Our method works best for small subhaloes (≲1/1000 host mass), for which the error caused by dynamical friction is minimal. We are able to capture the evaporation of subhalo particles by interactions with host halo particles, an effect that has not yet been fully explored in the context of subhalo core-collapse. We find three main processes drive subhalo evolution: subhalo internal heat outflow, host-subhalo evaporation, and tidal effects. The subhalo central density grows only when the heat outflow outweighs the energy gain from evaporation and tidal heating. Thus, evaporation delays or even disrupts subhalo core-collapse. We map out the parameter space for subhaloes to core-collapse, finding that it is nearly impossible to drive core-collapse in subhaloes in SIDM models with constant cross-sections. Any discovery of ultracompact dark substructures with future substructure lensing observations favours additional degrees of freedom, such as velocity-dependence, in the cross-section. 
    more » « less