---
title: Inside-Out Planet Formation (IOPF)
url: https://www.emergentmind.com/topics/inside-out-planet-formation-iopf
type: topic
---

# Inside-Out Planet Formation (IOPF)

Inside-Out Planet Formation (IOPF) is an in situ framework for the origin of Systems with Tightly-packed Inner Planets (STIPs), in which close-in Earth- to super-Earth-scale planets are assembled sequentially at a pressure maximum associated with the dead-zone inner boundary of a protoplanetary disk. In the canonical picture, pebbles produced farther out drift inward by gas drag, accumulate at the innermost pressure trap, form the first planet, and then cease to feed that planet once gap opening displaces the trap and induces outward retreat of the dead-zone boundary; the same cycle then repeats at progressively larger radii. The model was introduced to account for compact, coplanar, mostly non-resonant Kepler multi-planet systems without requiring either large-scale resonant migration chains or extremely massive inner disks [1306.0576], [1510.06703].

## 1. Origins, explanatory target, and relation to competing models

IOPF was formulated to address the architecture of STIPs: compact systems with several planets of roughly Earth to super-Earth scale on sub-AU, well-aligned orbits. In the review formulation, these systems typically show period ratios around \(1.5\)–3 and only a weak tendency to occupy exact mean-motion resonances, a pattern that is awkward for migration scenarios that naturally yield resonant chains and for classical inner-disk in situ models that begin by assuming a strongly enhanced solid reservoir already concentrated inside \(\sim 1\) AU [1510.06703]. The original proposal therefore reframed radial drift from a loss channel into a delivery mechanism: solids formed over a much larger disk are brought inward as pebbles and concentrated at a dynamically selected inner trap rather than being required to exist there from the outset [1306.0576].

Within this literature, two broad alternatives are used as foils. The first is formation at larger radii followed by inward migration; the principal concern is that such models often overproduce resonant chains, although the caveat is noted that low-mass planets may evade trapping and later processes may break resonances. The second is classical in situ formation from a very massive local inner disk; the challenge here is that the required solid surface densities can reach \(\gtrsim 20\times\) the minimum-mass solar nebula, which the papers describe as potentially difficult to reconcile with standard viscous disk theory and liable to raise self-gravity issues [1411.2629]. IOPF instead posits that the inner solid reservoir is built dynamically by pebble drift from a supply zone that can extend to \(\gtrsim 10\) AU [1411.2629].

A terminological feature of the series is the designation of the first, innermost planet as the “Vulcan” planet. The term does not denote a separate class of exoplanets; it identifies the first planet to form at the innermost trap, closest to the star, in analogy with the historical hypothetical planet interior to Mercury [1411.2629].

## 2. Dead-zone physics, MRI activation, and the inner pressure maximum

The organizing structure in IOPF is the transition between an inner MRI-active region and an outer, weakly ionized dead zone. In the early analytic formulation, the key thermal criterion is that the disk midplane reaches approximately
\[
T \simeq 1200~\mathrm{K},
\]
the temperature at which thermal ionization of alkali metals is taken to permit MRI activation. Because the effective viscosity rises inward across this transition, a steady accretion flow requires a corresponding drop in surface density interior to the boundary, producing a radial pressure maximum where inward-drifting pebbles stall and accumulate [1411.2629]. In the CT14/CT15 parameterization used for the Vulcan-planet analysis, the radius of the \(1200\) K front is
\[
r_{1200\,\mathrm{K}} = 0.178\, \phi_{\rm DZIB}\, \gamma_{1.4}^{-2/9} \kappa_{10}^{2/9} \alpha_{-3}^{-2/9} m_{*,1}^{1/3} \left(f_r\dot m_{-9}\right)^{4/9}\ \mathrm{AU},
\]
with a fiducial correction \(\phi_{\rm DZIB}=0.5\) introduced to absorb uncertainties such as protostellar heating and wind losses [1411.2629].

Later work replaced the ad hoc active-zone viscosity with MRI criteria based on non-ideal MHD diffusivities and pressure-weighted vertical averaging. That revision changed the geometric interpretation of the trap. In the MRI-coupled steady-state solutions, the pressure maximum does not generally coincide with the first radius where the midplane becomes dead. Instead, because MRI-active layers persist above the dead midplane, the vertically averaged \(\bar\alpha\) continues to decrease outward beyond the actual midplane DZIB, so the pressure maximum forms farther out, within the dead zone itself [1712.07049]. For the fiducial \(M_*=1\,M_\odot\), \(\dot M=10^{-9}\,M_\odot\,{\rm yr}^{-1}\), \(\alpha_{\rm DZ}=10^{-4}\) model, the true midplane DZIB lies at \(\sim 0.09\) AU whereas the pressure maximum lies at \(\sim 0.25\) AU [1712.07049]. This does not remove the IOPF trap; it relocates it and makes its position depend explicitly on the layered MRI structure.

The same MRI-based calculations introduced additional dependencies absent from the original shorthand picture. Hall resistivity dominates near the midplane in the fiducial inner disk, X-ray ionization can be competitive with thermal ionization in the low-\(\Sigma\) inner active layers, and the resulting steady inner-disk solutions are viscously unstable to surface-density perturbations [1712.07049]. A plausible implication is that the IOPF trap is more structurally delicate than a simple \(\alpha\)-jump model suggests, even if the basic existence of a pressure maximum remains intact.

## 3. Pebble delivery, ring formation, and the sequential inside-out cycle

IOPF begins with a pebble flux. In the original formulation, cm–m solids drift inward because gas is pressure-supported and therefore slightly sub-Keplerian, so solids feel a headwind and lose angular momentum [1306.0576]. For a pressure profile \(P=P_0(r/r_0)^{-k_P}\), the radial drift speed is written
\[
v_{r,p}\simeq -k_P(c_s/v_K)^2 (\tau_{\rm fric}+\tau_{\rm fric}^{-1})^{-1} v_K,
\]
with the drift time \(t_{\rm drift}=r/|v_{r,p}|\) [1306.0576]. At the pressure maximum, the pressure gradient changes sign and the drift stalls; provided the pebble is large enough that outward drift relative to the gas exceeds inward gas advection, it is trapped rather than carried into the star [1306.0576].

The earliest IOPF papers allowed two routes from the trapped ring to a planet: direct gravitational instability of the pebble ring or planet formation by core accretion from a pebble-rich ring. In the gravitationally unstable version, the Toomre fragment mass is of order \(\sim 10^{-3}\,M_\oplus\), but the relevant annular “Toomre ring mass” is of order \(\sim 1\,M_\oplus\) under fiducial conditions, so the unstable ring can plausibly consolidate into an Earth-mass body before later growth to the gap-opening scale [1306.0576], [1510.06703]. Subsequent papers increasingly emphasized continued pebble accretion as the dominant growth channel after the initial seed forms [1510.06703].

The supply problem was addressed explicitly in “Pebble Delivery for Inside-Out Planet Formation” [1410.5819]. In a fiducial \(\dot m=10^{-9}\,M_\odot\,{\rm yr}^{-1}\), \(\alpha=10^{-3}\) disk, fixed-size pebbles starting at 10 AU require about \(10^5\) yr for \(0.1\) cm, \(10^4\) yr for \(1\) cm, and \(10^3\) yr for \(10\) cm sizes to reach the DZIB region. When growth during drift is included, delivery accelerates sharply: a pebble that begins at \(0.1\) cm at \(16.5\) AU reaches the DZIB in \(1.15\times10^4\) yr rather than \(2.88\times10^5\) yr, and a \(0.1\) cm pebble starting at 100 AU can reach the inner disk in only a few \(\times 10^4\) yr [1410.5819]. The 2017 global pebble-evolution models extended this conclusion, finding that pebbles typically grow to a few cm during inward drift from several tens of AU and that producing realistic STIPs within disk lifetimes requires \(\dot m\sim10^{-9}\,M_\odot\,{\rm yr}^{-1}\) and relatively low dead-zone viscosity, \(\alpha\sim10^{-4}\) [1709.10130].

The sequential character of IOPF follows from how growth terminates. Once the first planet becomes massive enough to perturb the gas disk strongly, the local pressure maximum is displaced, fresh pebbles are intercepted outside the planet, and the dead-zone inner boundary retreats outward. Pebbles drifting inward from the outer disk then collect at the new pressure maximum, where the next ring and next planet form [1411.2629]. This ring-to-planet-to-retreat cycle is the defining “inside-out” logic of the model.

## 4. Planet trapping, gap opening, and characteristic mass–radius scalings

For IOPF to be observationally meaningful, the planet must remain near its formation site. Hydrodynamic simulations of the DZIB transition confirm this requirement. In a 2D accretion-heated disk with \(\alpha\) rising inward from \(\alpha_{\rm DZIB}=0.001\) at 0.1 AU to \(\alpha_{\rm MRI}=0.01\) at 0.07 AU, the total disk torque on fixed planets crosses from positive to negative across the transition, so the zero-torque point near the original pressure maximum is a stable planet trap. This holds from \(0.1\,M_G\) to \(1.0\,M_G\), where \(M_G\) is the analytic viscous gap-opening scale [1508.02791]. In the canonical IOPF interpretation, Type I migration is therefore strongly suppressed at the trap, and later Type II migration remains limited because the gap-opening planet can exceed the gas mass remaining interior to it [1411.2629].

The mass scale that ends growth is set by gap opening. In the Vulcan-planet analysis the gap-opening mass is written
\[
M_p=\phi_G \frac{40\nu m_*}{r^2\Omega_K},
\]
which becomes, in the steady viscous disk model,
\[
M_p \rightarrow 5.67\, \phi_{G,0.3}\, \gamma_{1.4}^{4/5} \kappa_{10}^{1/5} \alpha_{-3}^{4/5} m_{*,1}^{3/10} \left(f_r\dot m_{-9}\right)^{2/5}\left(\frac{r}{0.1\,\mathrm{AU}}\right)^{1/10} M_\oplus
\]
[1411.2629]. Setting \(r=r_{1200\,\mathrm{K}}\) and eliminating \(f_r\dot m\) yields the central Vulcan prediction
\[
M_{p,1}=5.0\, \phi_{G,0.3}\,\phi_{{\rm DZIB},0.5}^{-9/10}\,\gamma_{1.4}\,\alpha_{-3}\left(\frac{r}{0.1\,\mathrm{AU}}\right) M_\oplus,
\]
so that for fiducial parameters
\[
M_{p,1}\simeq 5.0\left(\frac{r}{0.1\,\mathrm{AU}}\right)M_\oplus.
\]
A notable feature of this derivation is that the explicit \(\dot m\), \(\kappa\), and \(m_*\) dependence cancels, leaving the normalization controlled mainly by \(\phi_G\), \(\phi_{\rm DZIB}\), \(\gamma\), and especially the dead-zone viscosity \(\alpha\) [1411.2629].

Later work revised the gap criterion. Using the Duffell-based formulation and new hydrodynamic calibration, the 2017 paper found that a pressure-maximum displacement exceeding one Hill radius at 0.1 AU requires \(M_p=3.86\,M_\oplus\) for the fiducial \(\dot m=10^{-9}\,M_\odot\,{\rm yr}^{-1}\), \(\alpha_{\rm DZ}=10^{-4}\) case, and that the corresponding Vulcan relation becomes
\[
M_{p,1}=3.50\, \phi_{\rm G,D,1.44}\,\phi_{{\rm DZIB},0.5}^{-9/8}\,\gamma_{1.4}^{5/4}\,\alpha_{-4}^{1/2}\,m_{*,1}^{-1/4}\,r_{0.1{\rm AU}}^{5/4}\,M_\oplus
\]
[1709.10130]. The IOPF literature therefore contains two related but not identical Vulcan scalings: an earlier linear law derived from the viscous-thermal gap criterion and a later \(r^{5/4}\) law derived from the Duffell-based calibration.

Planet spacing is likewise tied to disk response. In fixed-\(\alpha\) simulations, the azimuthally averaged pressure maximum is not displaced substantially until \(M_p\gtrsim0.5\,M_G\), at which point the offset jumps to \(\sim 5\,R_H\); in evolving-\(\alpha\) models that treat DZIB retreat heuristically through X-ray penetration, the new pressure maximum appears much farther out, \(\sim 25\)–68 Hill radii beyond a \(0.5\,M_G=5.59\,M_\oplus\) planet at 0.1 AU [1508.02791]. This large first retreat is central to the prediction that the first planet pair in STIPs should often be spaced more widely, in Hill units, than later pairs.

## 5. Empirical tests with Kepler systems

The most direct population-level test in the series concerns Vulcan planets in Kepler multis. The 2014 study selected 629 systems from the NASA Exoplanet Archive containing at least two transiting planets and analyzed only the innermost transiting planet in each system [1411.2629]. Because Kepler measures radii rather than masses, the comparison was carried out through Monte Carlo forward modeling: theoretical masses from the IOPF relation were converted to radii by drawing from lognormal density PDFs calibrated on planets with measured \(M_p\) and \(R_p\), detectability was imposed using host-star \(K_p\), combined differential photometric precision, a \(3.5\)-yr baseline, and \(\mathrm{SNR}>7\), and the resulting synthetic detected sample was compared to the observed one [1411.2629].

In radius space, the observed innermost planets follow
\[
R_{p,1}/R_\oplus \propto r^{0.3\pm0.2},
\]
reported more explicitly as
\[
R_{p,1}/R_\oplus=(3.5\pm0.5)\,r^{\,0.3\pm0.2},
\]
with \(r\) in AU, and the synthetic IOPF population gives the same \(R_{p,1}\propto r^{0.3\pm0.2}\) scaling after density scatter and Kepler selection are applied [1411.2629]. In mass space, the intrinsic theoretical prediction is \(M_{p,1}\propto r^1\), but radius-to-mass inference and selection effects flatten the recovered trend. The synthetic population yields approximately \(M_{p,1}\propto r^{0.9\pm0.2}\), while the observed Vulcan-planet scaling is approximately \(M_{p,1}\propto r^{0.7\pm0.2}\); the three empirical \(M_p\)-\(R_p\) prescriptions used in the body of the paper give observed relations with slopes \(0.56\pm0.17\), \(0.49\pm0.17\), and \(0.72\pm0.17\) [1411.2629]. Within that forward-model comparison, the IOPF prediction is deemed consistent with the Kepler data.

Normalization provides a separate diagnostic. In the Vulcan relation the normalization scales roughly linearly with \(\alpha\), and the paper finds that \(\alpha=10^{-3}\) predicts masses too high by a factor of a few, whereas \(\alpha=2\times10^{-4}\) brings both scaling and normalization into much better agreement with the observed innermost Kepler planets [1411.2629]. The 2015 overview recast this as a preferred \(\alpha_{-3}=0.205\), that is \(\alpha\simeq2.05\times10^{-4}\), for the observed Vulcan normalization [1510.06703].

Spacing comparisons are qualitatively similar. The evolving-\(\alpha\) hydrodynamic models produce first-to-second trap separations of \(25.3\), \(56.3\), and \(67.6\) Hill radii, overlapping the broad observed peak of roughly \(20\)–\(60\) in Kepler multis [1508.02791]. The review article likewise summarizes that \(\phi_{\Delta r,1}\) peaks at \(\sim 20\)–40, whereas later pair spacings peak at \(\lesssim 20\), consistent with the expectation that the first gap-opening event produces the largest DZIB retreat [1510.06703]. The caveat, stressed repeatedly, is that these are statistical tests contingent on assumptions about intrinsic density distributions, the identification of the observed innermost transiting planet with the true physical innermost planet, and the smallness of post-formation migration [1411.2629].

## 6. Onset conditions, chemical consequences, and unresolved problems

Recent extensions have shifted part of the focus from the mature sequential mechanism to the boundary conditions for its onset. In the 2025 installment, the DZIB structure from the MRI-based inner-disk models was combined with a pebble trapping criterion to ask when the first trap becomes efficient. For a fiducial \(m_*=1\,M_\odot\), \(\alpha=10^{-4}\) disk, the minimum trapping size scales as
\[
a_{p,{\rm trap}}\propto \rho_p^{-1}k_P^{-1}\gamma^{-2/3}\kappa^{-1/3}\alpha^{1/3}f_r^{-2/3}\dot m^{1/3},
\]
and the numerical calculations show thresholds of \(\simeq0.10\) cm at \(\dot m=10^{-8}\), \(\simeq0.063\) cm at \(10^{-9}\), and \(\simeq0.035\) cm at \(10^{-10}\,M_\odot\,{\rm yr}^{-1}\) [2508.10755]. The proposed onset of IOPF occurs when \(\dot m\sim10^{-9}\,M_\odot\,{\rm yr}^{-1}\), the DZIB lies near \(\sim0.1\)–0.2 AU, and typical pebbles are of order \(0.5\) mm, so that trapping first becomes efficient at the epoch required to produce a few-\(M_\oplus\) Vulcan planet. The same paper further suggests that this trapping transition may coincide with the first emergence of the transition-disk phase [2508.10755].

Chemical modeling added a different dimension. The 2022 astrochemical study followed gas advection, pebble drift, and gas-grain chemistry from 300 AU to the DZIB and found that in outer cool disk regions carbon and up to 90% of oxygen nuclei begin locked in CO and \(\mathrm{O}_2\) ice, with \(\mathrm{CO}_2\) and \(\mathrm{H}_2\mathrm{O}\) lower by about an order of magnitude [2202.02483]. Pebble drift then enhances gas-phase volatile abundances by up to two orders of magnitude at ice lines; inside \(\lesssim1\) AU the models yield water-rich gas with \(\mathrm{C/O}\lesssim0.1\), while solids delivered to the hot inner trap are predicted to be volatile-poor. In the IOPF interpretation, close-in planets formed near the DZIB should therefore have volatile-poor interiors but may accrete primordial atmospheres from oxygen-rich, water-rich gas, with volatile mass fractions of order \(\sim10\%\) becoming plausible only when the formation region approaches the water ice line [2202.02483].

The main internal dynamical challenge to the standard one-ring/one-planet interpretation came from direct \(N\)-body evolution of a DZIB planetesimal ring. In that study, a \(1\,M_\oplus\) ring at \(\sim0.1\) AU almost never collapsed to a single dominant body. Instead, it underwent oligarchic growth and typically ended with 2 or 3 surviving oligarchs on nearly coplanar, circular orbits; 84.5% of the 360 runs ended with exactly two planets, 15.3% with three, and only one with four [2104.03128]. The primary usually acquired about 70% of the mass, but the secondary retained \(\sim30\)–65% of the primary’s mass, yielding period-ratio, Hill-spacing, and mass–radius trends inconsistent with observed innermost STIP pairs. The paper therefore did not reject IOPF as a whole, but it directly challenged the assumption that a pebble ring can pass through a planetesimal-oligarchic phase and still naturally deliver a single Vulcan planet without additional physics such as gas torques, gas drag, pebble filtering, or gradual seed formation [2104.03128].

A different challenge arose from non-axisymmetric burst dynamics at the dead-zone inner edge. High-resolution 2D multifluid radiation-hydrodynamic simulations of the inner 10 AU showed that accretion outbursts are highly unstable to the Rossby-wave instability, generating numerous vortices that merge and drive an effective burst-phase turbulent stress \(\alpha_{\rm vort}\gtrsim10^{-2}\) across the burst region [2606.11344]. Under these conditions, planetesimal formation by streaming-instability-style concentration is strongly suppressed during the burst itself: a 1D model that would convert \(0.12\)–\(0.19\,M_\oplus\) of solids during burst conditions yields only \(0.004\)–\(0.04\,M_\oplus\) in the fiducial 2D run and effectively zero in the highest-resolution run [2606.11344]. The same simulations, however, show that the disk returns after the burst to a low-turbulence, quasi-axisymmetric state with a surviving pressure bump near \(\sim0.75\)–0.8 AU, so efficient seed formation is delayed rather than eliminated. This suggests an episodic IOPF variant in which the dead-zone-edge trap remains viable over secular times, but planetesimal formation is favored in quiescent intervals rather than continuously.

Taken together, these developments have refined rather than abolished IOPF. The framework still offers a specific disk-structure-regulated route to compact inner planetary systems, with explicit predictions for mass scales, radial trends, and first-pair spacings. At the same time, the modern literature has made clear that its viability depends on more than the existence of a pressure bump. The detailed location of the trap depends on layered MRI physics and non-ideal MHD; the onset of trapping depends on pebble size and accretion rate; the ring-to-planet step may fail if oligarchic growth is not circumvented; and burst-driven non-axisymmetry can suppress seed formation during active phases [1712.07049], [2508.10755], [2104.03128], [2606.11344]. IOPF is therefore best understood not as a single closed analytic model, but as a continuing research program centered on one robust proposition: inner-disk pressure traps tied to dead-zone physics can regulate the sequential, inside-out assembly of close-in planetary systems.

Source: https://www.emergentmind.com/topics/inside-out-planet-formation-iopf