skip to main content


Title: Stars Crushed by Black Holes. II. A Physical Model of Adiabatic Compression and Shock Formation in Tidal Disruption Events
Abstract

We develop a Newtonian model of a deep tidal disruption event (TDE), for which the pericenter distance of the star,rp, is well within the tidal radius of the black hole,rt, i.e., whenβrt/rp≫ 1. We find that shocks form forβ≳ 3, but they are weak (with Mach numbers ∼1) for allβ, and that they reach the center of the star prior to the time of maximum adiabatic compression forβ≳ 10. The maximum density and temperature reached during the TDE follow much shallower relations withβthan the previously predictedρmaxβ3andTmaxβ2scalings. Belowβ≃ 10, this shallower dependence occurs because the pressure gradient is dynamically significant before the pressure is comparable to the ram pressure of the free-falling gas, while aboveβ≃ 10, we find that shocks prematurely halt the compression and yield the scalingsρmaxβ1.62andTmaxβ1.12. We find excellent agreement between our results and high-resolution simulations. Our results demonstrate that, in the Newtonian limit, the compression experienced by the star is completely independent of the mass of the black hole. We discuss our results in the context of existing (affine) models, polytropic versus non-polytropic stars, and general relativistic effects, which become important when the pericenter of the star nears the direct capture radius, atβ∼ 12.5 (2.7) for a solar-like star disrupted by a 106M(107M) supermassive black hole.

 
more » « less
Award ID(s):
2006684
NSF-PAR ID:
10362619
Author(s) / Creator(s):
;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
926
Issue:
1
ISSN:
0004-637X
Format(s):
Medium: X Size: Article No. 47
Size(s):
["Article No. 47"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    State transitions in black hole X-ray binaries are likely caused by gas evaporation from a thin accretion disk into a hot corona. We present a height-integrated version of this process, which is suitable for analytical and numerical studies. With radiusrscaled to Schwarzschild units and coronal mass accretion rateṁcto Eddington units, the results of the model are independent of black hole mass. State transitions should thus be similar in X-ray binaries and an active galactic nucleus. The corona solution consists of two power-law segments separated at a break radiusrb∼ 103(α/0.3)−2, whereαis the viscosity parameter. Gas evaporates from the disk to the corona forr>rb, and condenses back forr<rb. Atrb,ṁcreaches its maximum,ṁc,max0.02(α/0.3)3. If atrrbthe thin disk accretes withṁd<ṁc,max, then the disk evaporates fully before reachingrb, giving the hard state. Otherwise, the disk survives at all radii, giving the thermal state. While the basic model considers only bremsstrahlung cooling and viscous heating, we also discuss a more realistic model that includes Compton cooling and direct coronal heating by energy transport from the disk. Solutions are again independent of black hole mass, andrbremains unchanged. This model predicts strong coronal winds forr>rb, and aT∼ 5 × 108K Compton-cooled corona forr<rb. Two-temperature effects are ignored, but may be important at small radii.

     
    more » « less
  2. Abstract

    M dwarfs are common host stars to exoplanets but often lack atmospheric abundance measurements. Late-M dwarfs are also good analogs to the youngest substellar companions, which share similarTeff∼ 2300–2800 K. We present atmospheric analyses for the M7.5 companion HIP 55507 B and its K6V primary star with Keck/KPIC high-resolution (R∼ 35,000)K-band spectroscopy. First, by including KPIC relative radial velocities between the primary and secondary in the orbit fit, we improve the dynamical mass precision by 60% and findMB=88.03.2+3.4MJup, putting HIP 55507 B above the stellar–substellar boundary. We also find that HIP 55507 B orbits its K6V primary star witha=383+4au ande= 0.40 ± 0.04. From atmospheric retrievals of HIP 55507 B, we measure [C/H] = 0.24 ± 0.13, [O/H] = 0.15 ± 0.13, and C/O = 0.67 ± 0.04. Moreover, we strongly detect13CO (7.8σsignificance) and tentatively detectH218O(3.7σsignificance) in the companion’s atmosphere and measure12CO/13CO=9822+28andH216O/H218O=24080+145after accounting for systematic errors. From a simplified retrieval analysis of HIP 55507 A, we measure12CO/13CO=7916+21andC16O/C18O=28870+125for the primary star. These results demonstrate that HIP 55507 A and B have consistent12C/13C and16O/18O to the <1σlevel, as expected for a chemically homogeneous binary system. Given the similar flux ratios and separations between HIP 55507 AB and systems with young substellar companions, our results open the door to systematically measuring13CO andH218Oabundances in the atmospheres of substellar or even planetary-mass companions with similar spectral types.

     
    more » « less
  3. Abstract

    We present measurements of the rest-frame UV spectral slope,β, for a sample of 36 faint star-forming galaxies atz∼ 9–16 discovered in one of the deepest JWST NIRCam surveys to date, the Next Generation Deep Extragalactic Exploratory Public Survey. We use robust photometric measurements for UV-faint galaxies (down toMUV∼ −16), originally published in Leung et al., and measure values of the UV spectral slope via photometric power-law fitting to both the observed photometry and stellar population models obtained through spectral energy distribution (SED) fitting withBagpipes. We obtain a median and 68% confidence interval forβfrom photometric power-law fitting ofβPL=2.70.5+0.5and from SED fitting,βSED=2.30.1+0.2for the full sample. We show that when only two to three photometric detections are available, SED fitting has a lower scatter and reduced biases than photometric power-law fitting. We quantify this bias and find that after correction the medianβSED,corr=2.50.2+0.2. We measure physical properties for our galaxies withBagpipesand find that our faint (MUV=18.10.9+0.7) sample is low in mass (log[M*/M]=7.70.5+0.5), fairly dust-poor (Av=0.10.1+0.2mag), and modestly young (log[age]=7.80.8+0.2yr) with a median star formation rate oflog(SFR)=0.30.4+0.4Myr1. We find no strong evidence for ultrablue UV spectral slopes (β∼ −3) within our sample, as would be expected for exotically metal-poor (Z/Z< 10−3) stellar populations with very high Lyman continuum escape fractions. Our observations are consistent with model predictions that galaxies of these stellar masses atz∼ 9–16 should have only modestly low metallicities (Z/Z∼ 0.1–0.2).

     
    more » « less
  4. Abstract

    We present a Keck/MOSFIRE rest-optical composite spectrum of 16 typical gravitationally lensed star-forming dwarf galaxies at 1.7 ≲z≲ 2.6 (zmean= 2.30), all chosen independent of emission-line strength. These galaxies have a median stellar mass oflog(M*/M)med=8.290.43+0.51and a median star formation rate ofSFRHαmed=2.251.26+2.15Myr1. We measure the faint electron-temperature-sensitive [Oiii]λ4363 emission line at 2.5σ(4.1σ) significance when considering a bootstrapped (statistical-only) uncertainty spectrum. This yields a direct-method oxygen abundance of12+log(O/H)direct=7.880.22+0.25(0.150.06+0.12Z). We investigate the applicability at highzof locally calibrated oxygen-based strong-line metallicity relations, finding that the local reference calibrations of Bian et al. best reproduce (≲0.12 dex) our composite metallicity at fixed strong-line ratio. At fixedM*, our composite is well represented by thez∼ 2.3 direct-method stellar mass—gas-phase metallicity relation (MZR) of Sanders et al. When comparing to predicted MZRs from the IllustrisTNG and FIRE simulations, having recalculated our stellar masses with more realistic nonparametric star formation histories(log(M*/M)med=8.920.22+0.31), we find excellent agreement with the FIRE MZR. Our composite is consistent with no metallicity evolution, at fixedM*and SFR, of the locally defined fundamental metallicity relation. We measure the doublet ratio [Oii]λ3729/[Oii]λ3726 = 1.56 ± 0.32 (1.51 ± 0.12) and a corresponding electron density ofne=10+215cm3(ne=10+74cm3) when considering the bootstrapped (statistical-only) error spectrum. This result suggests that lower-mass galaxies have lower densities than higher-mass galaxies atz∼ 2.

     
    more » « less
  5. Abstract

    We present13CO(J= 1 → 0) observations for the EDGE-CALIFA survey, which is a mapping survey of 126 nearby galaxies at a typical spatial resolution of 1.5 kpc. Using detected12CO emission as a prior, we detect13CO in 41 galaxies via integrated line flux over the entire galaxy and in 30 galaxies via integrated line intensity in resolved synthesized beams. Incorporating our CO observations and optical IFU spectroscopy, we perform a systematic comparison between the line ratio12/13I[12CO(J=10)]/I[13CO(J=10)]and the properties of the stars and ionized gas. Higher12/13values are found in interacting galaxies compared to those in noninteracting galaxies. The global12/13slightly increases with infrared colorF60/F100but appears insensitive to other host-galaxy properties such as morphology, stellar mass, or galaxy size. We also present azimuthally averaged12/13profiles for our sample up to a galactocentric radius of 0.4r25(∼6 kpc), taking into account the13CO nondetections by spectral stacking. The radial profiles of12/13are quite flat across our sample. Within galactocentric distances of 0.2r25, the azimuthally averaged12/13increases with the star formation rate. However, Spearman rank correlation tests show the azimuthally averaged12/13does not strongly correlate with any other gas or stellar properties in general, especially beyond 0.2r25from the galaxy centers. Our findings suggest that in the complex environments in galaxy disks,12/13is not a sensitive tracer for ISM properties. Dynamical disturbances, like galaxy interactions or the presence of a bar, also have an overall impact on12/13, which further complicates the interpretations of12/13variations.

     
    more » « less