skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: The maximum accretion rate of a protoplanet: how fast can runaway be?
ABSTRACT The hunt is on for dozens of protoplanets hypothesized to reside in protoplanetary discs with imaged gaps. How bright these planets are, and what they will grow to become, depend on their accretion rates, which may be in the runaway regime. Using 3D global simulations, we calculate maximum gas accretion rates for planet masses Mp from 1$$\, \mathrm{ M}_{{\oplus }}$$ to $$10\, \mathrm{ M}_{\rm J}$$. When the planet is small enough that its sphere of influence is fully embedded in the disc, with a Bondi radius rBondi smaller than the disc’s scale height Hp – such planets have thermal mass parameters qth ≡ (Mp/M⋆)/(Hp/Rp)3 ≲ 0.3, for host stellar mass M⋆ and orbital radius Rp – the maximum accretion rate follows a Bondi scaling, with $$\max \dot{M}_{\rm p} \propto \rho _{\rm g}M_{\rm p}^2 / (H_{\rm p}/R_{\rm p})^3$$ for ambient disc density ρg. For more massive planets with 0.3 ≲ qth ≲ 10, the Hill sphere replaces the Bondi sphere as the gravitational sphere of influence, and $$\max \dot{M}_{\rm p} \propto \rho _{\rm g}M_{\rm p}^1$$, with no dependence on Hp/Rp. In the strongly superthermal limit when qth ≳ 10, the Hill sphere pops well out of the disc, and $$\max \dot{M}_{\rm p} \propto \rho _{\rm g}M_{\rm p}^{2/3} (H_{\rm p}/R_{\rm p})^1$$. Applied to the two confirmed protoplanets PDS 70b and c, our numerically calibrated maximum accretion rates imply that their Jupiter-like masses may increase by up to a factor of ∼2 before their parent disc dissipates.  more » « less
Award ID(s):
2205500
PAR ID:
10626946
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
MNRAS
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
525
Issue:
2
ISSN:
0035-8711
Page Range / eLocation ID:
2806 to 2819
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; ; ; ; ; ; ; et al (, Monthly Notices of the Royal Astronomical Society)
    null (Ed.)
    ABSTRACT We report the discovery of a warm sub-Saturn, TOI-257b (HD 19916b), based on data from NASA’s Transiting Exoplanet Survey Satellite (TESS). The transit signal was detected by TESS and confirmed to be of planetary origin based on radial velocity observations. An analysis of the TESS photometry, the Minerva-Australis, FEROS, and HARPS radial velocities, and the asteroseismic data of the stellar oscillations reveals that TOI-257b has a mass of MP = 0.138 ± 0.023 $$\rm {M_J}$$ (43.9 ± 7.3 $$\, M_{\rm \oplus}$$), a radius of RP = 0.639 ± 0.013 $$\rm {R_J}$$ (7.16 ± 0.15 $$\, \mathrm{ R}_{\rm \oplus}$$), bulk density of $$0.65^{+0.12}_{-0.11}$$ (cgs), and period $$18.38818^{+0.00085}_{-0.00084}$$ $$\rm {days}$$. TOI-257b orbits a bright (V = 7.612 mag) somewhat evolved late F-type star with M* = 1.390 ± 0.046 $$\rm {M_{sun}}$$, R* = 1.888 ± 0.033 $$\rm {R_{sun}}$$, Teff = 6075 ± 90 $$\rm {K}$$, and vsin i = 11.3 ± 0.5 km s−1. Additionally, we find hints for a second non-transiting sub-Saturn mass planet on a ∼71 day orbit using the radial velocity data. This system joins the ranks of a small number of exoplanet host stars (∼100) that have been characterized with asteroseismology. Warm sub-Saturns are rare in the known sample of exoplanets, and thus the discovery of TOI-257b is important in the context of future work studying the formation and migration history of similar planetary systems. 
    more » « less
  2. ; ; ; ; (, Monthly Notices of the Royal Astronomical Society)
    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
  3. ; ; ; (, Monthly Notices of the Royal Astronomical Society)
    ABSTRACT We revisit the question of ‘hot mode’ versus ‘cold mode’ accretion on to galaxies using steady-state cooling flow solutions and idealized 3D hydrodynamic simulations. We demonstrate that for the hot accretion mode to exist, the cooling time is required to be longer than the free-fall time near the radius where the gas is rotationally supported, Rcirc, i.e. the existence of the hot mode depends on physical conditions at the galaxy scale rather than on physical conditions at the halo scale. When allowing for the depletion of the halo baryon fraction relative to the cosmic mean, the longer cooling times imply that a virialized gaseous halo may form in halo masses below the threshold of $$\sim 10^{12}\, {\rm M_{\odot }}$$ derived for baryon-complete haloes. We show that for any halo mass there is a maximum accretion rate for which the gas is virialized throughout the halo and can accrete via the hot mode of $${\dot{M}}_{\rm crit}\approx 0.7(v_{\rm c}/100\, \rm km\ s^{-1})^{5.4}(R_{\rm circ}/10\, {\rm kpc})(Z/\, {\rm Z_{\odot }})^{-0.9}\, {\rm M_{\odot }}\, {\rm yr}^{-1}$$, where Z and vc are the metallicity and circular velocity measured at Rcirc. For accretion rates $$\gtrsim {\dot{M}}_{\rm crit}$$ the volume-filling gas phase can in principle be ‘transonic’ – virialized in the outer halo but cool and free-falling near the galaxy. We compare $${\dot{M}}_{\rm crit}$$ to the average star formation rate (SFR) in haloes at 0 < z < 10 implied by the stellar-mass–halo-mass relation. For a plausible metallicity evolution with redshift, we find that $${\rm SFR}\lesssim {\dot{M}}_{\rm crit}$$ at most masses and redshifts, suggesting that the SFR of galaxies could be primarily sustained by the hot mode in halo masses well below the classic threshold of $$\sim 10^{12}\, {\rm M_{\odot }}$$. 
    more » « less
  4. ; ; ; ; (, Monthly Notices of the Royal Astronomical Society)
    ABSTRACT Despite many methods developed to find young massive planets in protoplanetary discs, it is challenging to directly detect low-mass planets that are embedded in discs. On the other hand, the core-accretion theory suggests that there could be a large population of embedded low-mass young planets at the Kelvin-Helmholtz (KH) contraction phase. We adopt both 1D models and 3D simulations to calculate the envelopes around low-mass cores (several to tens of M⊕) with different luminosities, and derive their thermal fluxes at radio wavelengths. We find that, when the background disc is optically thin at radio wavelengths, radio observations can see through the disc and probe the denser envelope within the planet’s Hill sphere. When the optically thin disc is observed with the resolution reaching one disc scale height, the radio thermal flux from the planetary envelope around a 10 M⊕ core is more than 10 per cent higher than the flux from the background disc. The emitting region can be extended and elongated. Finally, our model suggests that the au-scale clump at 52 au in the TW Hydrae disc revealed by ALMA is consistent with the envelope of an embedded 10–20 M⊕ planet, which can explain the detected flux, the spectral index dip, and the tentative spirals. The observation is also consistent with the planet undergoing pebble accretion. Future ALMA and ngVLA observations may directly reveal more such low-mass planets, enabling us to study core growth and even reconstruct the planet formation history using the embedded ‘protoplanet’ population. 
    more » « less
  5. ; ; ; ; ; ; ; ; ; ; et al (, Monthly Notices of the Royal Astronomical Society)
    ABSTRACT We apply the barred Schwarzschild method developed by Tahmasebzadeh et al. (2022) to a barred S0 galaxy, NGC 4371, observed by IFU instruments from the TIMER and ATLAS3D projects. We construct the gravitational potential by combining a fixed black hole mass, a spherical dark matter halo, and stellar mass distribution deprojected from 3.6 μm S$^4$G image considering an axisymmetric disc and a triaxial bar. We independently modelled kinematic data from TIMER and ATLAS3D. Both models fit the data remarkably well. We find a consistent bar pattern speed from the two sets of models with $$\Omega _{\rm p} = 23.6 \pm 2.8 \, \mathrm{km \, s^{-1} \, kpc^{-1} }$$ and $$\Omega _{\rm p} = 22.4 \pm 3.5 \, \mathrm{km \, s^{-1} \, kpc^{-1} }$$, respectively. The dimensionless bar rotation parameter is determined to be $$\mathcal {R} \equiv R_{\rm cor}/R_{\rm bar}=1.88 \pm 0.37$$, indicating a likely slow bar in NGC 4371. Additionally, our model predicts a high amount of dark matter within the bar region ($$M_{\rm DM}/ M_{\rm total}$$\sim 0.51 \pm 0.06$$), which, aligned with the predictions of cosmological simulations, indicates that fast bars are generally found in baryon-dominated discs. Based on the best-fitting model, we further decompose the galaxy into multiple 3D orbital structures, including a BP/X bar, a classical bulge, a nuclear disc, and a main disc. The BP/X bar is not perfectly included in the input 3D density model, but BP/X-supporting orbits are picked through the fitting to the kinematic data. This is the first time a real barred galaxy has been modelled utilizing the Schwarzschild method including a 3D bar. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Text Encoded Conversion
  • XML
  • Save / Share this Record
  • Send to Email