Hill Regime Accretion in Astrophysical Disks

• Hill regime accretion is the process where a body’s gravitational influence, defined by its Hill sphere, dominates over stellar tides in shear-dominated disks.
• The accretion rate is regulated by encounter kinematics and local transport physics, influencing pebble capture, gas accretion, and planetesimal growth.
• Applications span diverse systems, from protoplanetary and circumplanetary disks to ring moons and AGN environments, each governed by unique dissipative processes.

Hill regime accretion is the class of accretion processes in shearing astrophysical disks for which the relevant gravitational sphere of influence is the Hill sphere, $R_{\rm H}=a\left(M/(3M_\ast)\right)^{1/3}$, rather than a purely Bondi-type sphere set only by pressure support. In this regime, encounter kinematics are set primarily by local Keplerian shear, the characteristic relative speed is the Hill velocity $v_{\rm H}=\Omega R_{\rm H}$, and mass delivery is controlled by the structure and transport physics inside or across the Hill sphere. The phrase is used in several related literatures: pebble capture by planetary cores, gas accretion by superthermal protoplanets, trans-Hill planetesimal growth, collision-dominated inner-disk accretion, ring-moon accretion in Saturn’s rings, and Hill-limited gas retention by compact objects embedded in AGN disks (Lambrechts et al., 2012, Xu et al., 2017, Li et al., 2020, McKernan et al., 2019).

1. Common dynamical framework

The basic scales are the Hill radius,

$R_{\rm H}=a\left(\frac{M}{3M_\ast}\right)^{1/3},$

the Hill velocity,

$v_{\rm H}=\Omega R_{\rm H},$

and the Bondi radius,

$R_{\rm B}=\frac{GM}{c_s^2}$

or, in pebble-accretion notation with a headwind $\Delta v$,

$R_{\rm B}=\frac{GM}{\Delta v^2}.$

The Hill sphere is the region in which the accretor’s gravity dominates over stellar tides in the rotating frame; the Bondi radius is the region in which gravity dominates over pressure or headwind-driven drift. Hill-regime accretion begins when the gravitational sphere of influence is effectively truncated by stellar tides, so that capture and supply are regulated by the Hill geometry and by shear-driven encounters rather than by purely embedded Bondi flow (Xu et al., 2017, Lambrechts et al., 2012).

In gas-accretion studies, the thermal limit is the condition $r_{\rm H}\approx H$. One formulation writes

$M_{\rm th}\simeq 3h^3M_\ast,$

with $h\equiv H/a$, while another introduces

$v_{\rm H}=\Omega R_{\rm H}$0

The 3D global simulations of runaway gas inflow distinguish three ranges: $v_{\rm H}=\Omega R_{\rm H}$1 for Bondi-embedded accretion, $v_{\rm H}=\Omega R_{\rm H}$2 for 3D Hill accretion, and $v_{\rm H}=\Omega R_{\rm H}$3 for a strongly superthermal regime in which the Hill sphere protrudes vertically out of the disk and inflow becomes effectively 2D (Choksi et al., 2023). In pebble accretion, the analogous transition is set by the “transition mass”

$v_{\rm H}=\Omega R_{\rm H}$4

or, in the Lambrechts–Johansen notation,

$v_{\rm H}=\Omega R_{\rm H}$5

defined by $v_{\rm H}=\Omega R_{\rm H}$6 or $v_{\rm H}=\Omega R_{\rm H}$7 (Xu et al., 2017, Lambrechts et al., 2012).

A second organizing distinction is geometric. Pebble accretion is “2D” when the pebble layer thickness is smaller than the capture radius or Hill radius, and “3D” when it is thicker. Gas accretion onto giant planets is likewise 3D when the Hill sphere remains embedded in the disk thickness, but becomes vertically limited when $v_{\rm H}=\Omega R_{\rm H}$8. These dimensional transitions alter the mass scaling even when the Hill sphere remains the governing gravitational scale (Xu et al., 2017, Choksi et al., 2023).

2. Pebble accretion onto planetary cores

In the pebble-accretion literature, the Hill regime is the high-efficiency mode entered once $v_{\rm H}=\Omega R_{\rm H}$9. For marginally coupled pebbles, $R_{\rm H}=a\left(\frac{M}{3M_\ast}\right)^{1/3},$0, particles entering the Hill sphere are captured with near-unit efficiency in laminar flow because gas drag dissipates energy during the $R_{\rm H}=a\left(\frac{M}{3M_\ast}\right)^{1/3},$1 encounter. The canonical 2D Hill rate is

$R_{\rm H}=a\left(\frac{M}{3M_\ast}\right)^{1/3},$2

while the 3D rate is reduced by the factor $R_{\rm H}=a\left(\frac{M}{3M_\ast}\right)^{1/3},$3 when $R_{\rm H}=a\left(\frac{M}{3M_\ast}\right)^{1/3},$4. For tightly coupled particles the effective Hill cross section shrinks, with

$R_{\rm H}=a\left(\frac{M}{3M_\ast}\right)^{1/3},$5

so that $R_{\rm H}=a\left(\frac{M}{3M_\ast}\right)^{1/3},$6 (Lambrechts et al., 2012).

MRI turbulence does not erase the Hill-regime capture mechanism. In shearing-box simulations, accretion for $R_{\rm H}=a\left(\frac{M}{3M_\ast}\right)^{1/3},$7 remains intrinsically efficient even in strong ideal-MHD turbulence once the pebble layer thickness is normalized out. The paper defines

$R_{\rm H}=a\left(\frac{M}{3M_\ast}\right)^{1/3},$8

and a modified 2D/3D normalization

$R_{\rm H}=a\left(\frac{M}{3M_\ast}\right)^{1/3},$9

and finds comparable $v_{\rm H}=\Omega R_{\rm H},$0 in hydro, ambipolar-diffusion, and ideal-MHD runs. Turbulence broadens the feeding zone and lowers the capture probability within the Hill sphere, but these effects largely cancel; the remaining change in the absolute rate is primarily geometric, through the turbulence-set pebble layer thickness

$v_{\rm H}=\Omega R_{\rm H},$1

(Xu et al., 2017).

Analytical generalizations retain the same Hill control but replace monodisperse capture by size-integrated rates. Lyra et al. derive an exact monodisperse bridging formula between 3D and 2D,

$v_{\rm H}=\Omega R_{\rm H},$2

with the 2D Hill limit

$v_{\rm H}=\Omega R_{\rm H},$3

For a polydisperse MRN distribution with $v_{\rm H}=\Omega R_{\rm H},$4, $v_{\rm H}=\Omega R_{\rm H},$5, and $v_{\rm H}=\Omega R_{\rm H},$6, the exact Hill-regime correction is

$v_{\rm H}=\Omega R_{\rm H},$7

so the Hill rate is reduced by an exact factor $v_{\rm H}=\Omega R_{\rm H},$8 relative to the monodisperse case (Lyra et al., 2023).

Radiative thermodynamics does not necessarily alter Hill-regime pebble capture. For a $v_{\rm H}=\Omega R_{\rm H},$9 embryo at $R_{\rm B}=\frac{GM}{c_s^2}$0, $R_{\rm B}=\frac{GM}{c_s^2}$1 under the disk conditions studied, so the embryo resides in the Hill regime. Radiative transfer lowers temperatures in the outer Hill sphere by $R_{\rm B}=\frac{GM}{c_s^2}$2 and raises inner densities by up to a factor of two, while convective velocity amplitudes increase by roughly an order of magnitude near the embryo; nevertheless the measured pebble accretion rates remain essentially identical to the purely convective runs and preserve the linear scaling $R_{\rm B}=\frac{GM}{c_s^2}$3 over $R_{\rm B}=\frac{GM}{c_s^2}$4 to $R_{\rm B}=\frac{GM}{c_s^2}$5 (Popovas et al., 2018).

3. Gas accretion by giant planets in the Hill regime

For gas giants, Hill-regime accretion is not simply geometric capture from the full Hill cross section. After partial gap opening, residual gas diffuses into the co-orbital region, settles onto horseshoe streamlines, and enters the Hill sphere near the Lagrange channels. Inside $R_{\rm B}=\frac{GM}{c_s^2}$6, the incoming horseshoe flow encounters the outer circumplanetary disk (CPD), where gas is on nearly circular, shear-dominated orbits around the planet. Above the thermal mass, stellar tidal and Coriolis forces deflect the converging streamlines, and the interface between horseshoe flow and outer CPD develops discontinuities in vortensity and Bernoulli constant. In the inviscid limit these quantities are conserved along streamlines, so the interface acts as a “tidal barrier” that impedes free mixing and makes accretion diffusion-limited rather than purely supply-limited (Li et al., 2020).

The corresponding semi-analytic rate is written as

$R_{\rm B}=\frac{GM}{c_s^2}$7

with partial-gap depletion

$R_{\rm B}=\frac{GM}{c_s^2}$8

and a Hill-regime inflow coefficient

$R_{\rm B}=\frac{GM}{c_s^2}$9

The exponential term is the tidal barrier. It becomes severe when $\Delta v$0, so low-viscosity superthermal planets experience an exponential throttling,

$\Delta v$1

This is the principal departure from the classical supply-limited estimate $\Delta v$2, which overestimates accretion above the thermal limit because it ignores CPD coupling and the barrier (Li et al., 2020).

The mass dependence is therefore sensitive to viscosity and thickness. In weakly viscous, modest-thickness disks with $\Delta v$3 and $\Delta v$4, growth stalls at $\Delta v$5 because $\Delta v$6; in strongly viscous and/or thick disks with $\Delta v$7 and $\Delta v$8, planets can reach several $\Delta v$9. Two-dimensional calculations show eccentric horseshoe instabilities and accretion surges for very massive planets, whereas eccentric streamlines remain stable in 3D over $R_{\rm B}=\frac{GM}{\Delta v^2}.$0 orbits, keeping the barrier effective (Li et al., 2020).

A complementary 3D global calibration measures maximum one-way inflow into a sphere around the planet. It finds

$R_{\rm B}=\frac{GM}{\Delta v^2}.$1

for $R_{\rm B}=\frac{GM}{\Delta v^2}.$2, with $R_{\rm B}=\frac{GM}{\Delta v^2}.$3 without sink cells and $R_{\rm B}=\frac{GM}{\Delta v^2}.$4 with sink cells, and

$R_{\rm B}=\frac{GM}{\Delta v^2}.$5

for $R_{\rm B}=\frac{GM}{\Delta v^2}.$6, with $R_{\rm B}=\frac{GM}{\Delta v^2}.$7. These are strict upper limits to the true accretion because they quantify inflow into the Hill sphere rather than permanent accretion through the CPD and envelope (Choksi et al., 2023).

4. Circumplanetary delivery, envelope structure, and reduced-order prescriptions

Local 3D hydrodynamic simulations of gas delivery to circumplanetary disks separate a Bondi-dominated low-mass regime from a Hill-dominated high-mass regime at roughly $R_{\rm B}=\frac{GM}{\Delta v^2}.$8 at $R_{\rm B}=\frac{GM}{\Delta v^2}.$9. In the high-mass regime the CPD radius scales as

$r_{\rm H}\approx H$0

gas accretes vertically from above and below the midplane, and the radially averaged accretion-band width obeys

$r_{\rm H}\approx H$1

for $r_{\rm H}\approx H$2. A robust result is that

$r_{\rm H}\approx H$3

over $r_{\rm H}\approx H$4. When combined with the Kanagawa gap formula,

$r_{\rm H}\approx H$5

this yields a semi-analytical 3D CPD accretion law with the same $r_{\rm H}\approx H$6 dependence as the earlier 2D case but a normalization lower by $r_{\rm H}\approx H$7 (Maeda et al., 2022).

Non-isothermal 3D simulations of a $r_{\rm H}\approx H$8 planet at $r_{\rm H}\approx H$9 add a thermodynamic boundary inside the Hill sphere. In that model, $M_{\rm th}\simeq 3h^3M_\ast,$0 and the flow is a 3D, vertically dominated inflow that becomes supersonic near the Hill sphere. The inflow halts at an ionization surface defined by $M_{\rm th}\simeq 3h^3M_\ast,$1, with quasi-steady values $M_{\rm th}\simeq 3h^3M_\ast,$2 without feedback and $M_{\rm th}\simeq 3h^3M_\ast,$3 with radiative feedback. The mean accretion rate into the ionized envelope rises from $M_{\rm th}\simeq 3h^3M_\ast,$4 to $M_{\rm th}\simeq 3h^3M_\ast,$5, a gain factor of $M_{\rm th}\simeq 3h^3M_\ast,$6, while only $M_{\rm th}\simeq 3h^3M_\ast,$7 of gas crossing $M_{\rm th}\simeq 3h^3M_\ast,$8 reaches $M_{\rm th}\simeq 3h^3M_\ast,$9 without feedback and $h\equiv H/a$0 with feedback. This redefines the effective accretion boundary: line emission and accretion luminosity should be anchored to $h\equiv H/a$1 rather than to $h\equiv H/a$2 (Montesinos et al., 19 Mar 2025).

One-dimensional evolutionary models calibrated to 3D hydrodynamics implement disc-limited Hill accretion through an inverse growth time

$h\equiv H/a$3

with recommended parameters $h\equiv H/a$4, $h\equiv H/a$5, and $h\equiv H/a$6. The gas sink is applied cell-by-cell across the feeding zone,

$h\equiv H/a$7

so mass is removed from the same radial location from which it is accreted. In the calibrated Hill regime this prescription reproduces 3D hydro growth and migration tracks at the $h\equiv H/a$8 level across the parameter sets tested (Schib et al., 2022).

5. Trans-Hill growth, oligarchy, and collision-dominated inner disks

In coagulation theory, “trans-Hill” denotes the stage at which the bodies dominating viscous stirring satisfy $h\equiv H/a$9, where $v_{\rm H}=\Omega R_{\rm H}$00 is the random speed of the small bodies. The trans-Hill radius is

$v_{\rm H}=\Omega R_{\rm H}$01

with $v_{\rm H}=\Omega R_{\rm H}$02. After onset, $v_{\rm H}=\Omega R_{\rm H}$03 and the characteristic stirrer size evolve in lockstep,

$v_{\rm H}=\Omega R_{\rm H}$04

and the size spectrum below $v_{\rm H}=\Omega R_{\rm H}$05 becomes self-similar. In the collisionless trans-Hill substage the efficiency remains very low, with $v_{\rm H}=\Omega R_{\rm H}$06 at the peak and differential slope $v_{\rm H}=\Omega R_{\rm H}$07; in the collisional trans-Hill substage, entered at

$v_{\rm H}=\Omega R_{\rm H}$08

the efficiency rises as $v_{\rm H}=\Omega R_{\rm H}$09 and the differential slope becomes $v_{\rm H}=\Omega R_{\rm H}$10. Termination occurs either by mutual accretion of big bodies, at

$v_{\rm H}=\Omega R_{\rm H}$11

or by oligarchy when the oligarchy parameter

$v_{\rm H}=\Omega R_{\rm H}$12

approaches unity (Lithwick, 2013).

The approach to Hill-regime pebble accretion is a distinct bottleneck. In a streaming-instability filament model, the transition mass is

$v_{\rm H}=\Omega R_{\rm H}$13

and in the irradiated disk used there it scales as $v_{\rm H}=\Omega R_{\rm H}$14. The paper finds that planetesimal accretion reaches the Hill-transition mass only inside $v_{\rm H}=\Omega R_{\rm H}$15 around a solar-mass star within $v_{\rm H}=\Omega R_{\rm H}$16, while beyond $v_{\rm H}=\Omega R_{\rm H}$17 there is little or no collisional growth. This implies that, in cold giant-planet formation zones, pebble accretion must act directly on the most massive planetesimals already present in the initial mass function, with little or no help from mutual planetesimal collisions (Lorek et al., 2022).

At very short orbital periods, Hill-regime accretion enters a different limit because the physical body occupies a substantial fraction of its Hill sphere. For a rocky body,

$v_{\rm H}=\Omega R_{\rm H}$18

so $v_{\rm H}=\Omega R_{\rm H}$19 is mass-independent and scales as $v_{\rm H}=\Omega R_{\rm H}$20. Direct N-body simulations show that when $v_{\rm H}=\Omega R_{\rm H}$21 the system transitions to a collision-dominated Hill regime in which most close encounters end in mergers rather than scattering, accretion becomes nearly $v_{\rm H}=\Omega R_{\rm H}$22 percent efficient, embryos grow beyond the classical isolation mass, and the familiar oligarch/planetesimal bimodality does not develop. For realistic rocky densities the inferred boundary is $v_{\rm H}=\Omega R_{\rm H}$23 days around a $v_{\rm H}=\Omega R_{\rm H}$24 star (Wallace et al., 2023).

6. Ring moons, AGN disks, and observational consequences

In Saturn’s rings, Hill-regime accretion governs the growth of equatorial ridges on embedded moons when ring material approaches through the Hill sphere with shear-driven speeds $v_{\rm H}=\Omega R_{\rm H}$25. The relevant eccentricity scale is the Hill eccentricity,

$v_{\rm H}=\Omega R_{\rm H}$26

and synchronous rotation imposes a torque limit. The accretion torque is

$v_{\rm H}=\Omega R_{\rm H}$27

while the Saturnian tidal torque is

$v_{\rm H}=\Omega R_{\rm H}$28

Equating them gives a critical accretion rate

$v_{\rm H}=\Omega R_{\rm H}$29

For Pan, Atlas, and Daphnis the paper quotes $v_{\rm H}=\Omega R_{\rm H}$30, $v_{\rm H}=\Omega R_{\rm H}$31, and $v_{\rm H}=\Omega R_{\rm H}$32, with minimum ridge-growth times of $v_{\rm H}=\Omega R_{\rm H}$33, $v_{\rm H}=\Omega R_{\rm H}$34, and $v_{\rm H}=\Omega R_{\rm H}$35, and argues that the observed non-axisymmetric ridges imply actual durations $v_{\rm H}=\Omega R_{\rm H}$36 (Quillen et al., 2020).

In AGN disks, Hill-regime accretion appears around stellar-mass black holes embedded in the disk of a supermassive black hole. The pre-merger gas reservoir is limited by the Hill sphere rather than by the Bondi sphere because shear truncates capture at

$v_{\rm H}=\Omega R_{\rm H}$37

After a gravitational-wave kick, ram-pressure stripping of this Hill-sphere gas produces a prompt electromagnetic transient if $v_{\rm H}=\Omega R_{\rm H}$38. The onset time scales as

$v_{\rm H}=\Omega R_{\rm H}$39

the duration is $v_{\rm H}=\Omega R_{\rm H}$40, and the luminosity is

$v_{\rm H}=\Omega R_{\rm H}$41

If $v_{\rm H}=\Omega R_{\rm H}$42, the flare is delayed and weakened by photon diffusion through the disk atmosphere (McKernan et al., 2019).

Across these applications, Hill-regime accretion has a common structural meaning: stellar tides and disk shear set the relevant sphere of influence, but the mass-growth law depends on the transport physics specific to the system. In pebble accretion, that dependence is set by stopping time and pebble-layer thickness; in giant-planet gas accretion, by gap depletion, CPD coupling, viscosity, and, in some models, a thermodynamic boundary inside the Hill sphere; in planetesimal growth, by the relation between $v_{\rm H}=\Omega R_{\rm H}$43 and $v_{\rm H}=\Omega R_{\rm H}$44; and in ring or AGN environments, by torque constraints or ram-pressure stripping. This suggests that “Hill regime accretion” is not a single rate law but a family of shear-dominated accretion problems linked by a common dynamical scale, $v_{\rm H}=\Omega R_{\rm H}$45, and differentiated by the dissipative and supply processes that operate within it.