Low-Recycling Regime: Plasma and Planetary Dynamics

• Low-Recycling Regime is defined by suppressed exchange with boundaries, resulting in altered transport properties and unique operational conditions.
• In tokamak SOLs, low recycling produces high edge temperatures and low densities as neutrals from targets are intercepted by low-recycling side walls.
• In protoplanetary systems, insufficient gas exchange leads to radiative cooling and quasi-1D contraction of the inner envelope.

Searching arXiv for the cited papers and related usage of “low-recycling regime.” arXiv_search("(Moldenhauer et al., 2021)")

Low-recycling regime denotes a class of operating conditions in which exchange with a boundary or external reservoir is sufficiently weak that the usual replenishment cycle is suppressed. The term has distinct technical meanings in different subfields. In tokamak edge and scrape-off-layer (SOL) physics, it refers to the case where plasma particles hitting material surfaces are not efficiently returned as neutrals, producing a hot, rarefied, long-mean-free-path SOL (Shukla et al., 12 Nov 2025). In planetary-atmosphere formation, the phrase is not used explicitly in the source paper, but the framework implies a converse to efficient atmosphere–disk recycling: a state in which gas exchange with the disk is too slow or too shallow to offset radiative cooling, so parts of the envelope contract quasi-spherically rather than remaining in a recycling-supported steady state (Moldenhauer et al., 2021).

1. Terminology and domain-specific meanings

The technical content of “recycling” depends on the open system under consideration. In tokamak SOL studies, recycling concerns the return of neutrals from plasma-facing components; in protoplanetary-envelope studies, recycling concerns hydrodynamic replacement of cooled atmospheric gas by high-entropy disk gas.

Domain	Recycling process	Low-recycling condition
Tokamak SOL	Ions hit surfaces, become neutrals, and re-enter plasma	$R \ll 1$, weak neutral return, hot low-density SOL
Protoplanetary envelope	Low-entropy envelope gas is exchanged with high-entropy disk gas	Inferred case where exchange cannot offset cooling

In the tokamak usage, the recycling coefficient $R$ parameterizes surface return. The cited STEP study gives $R_{\text{target}}=0.99$ for the divertor plate and $R_{\text{wall}}=0.5$ for lithium-coated side walls, and associates low recycling with high edge/SOL temperatures, low edge/SOL densities, and a long-mean-free-path regime in which fluid closures become questionable (Shukla et al., 12 Nov 2025). By contrast, the planetary paper studies an efficient, deep-reaching recycling regime for a $1\,M_\oplus$ core at $0.1\,\mathrm{au}$; a “low-recycling regime” there is an inferred extrapolation in which the exchange becomes too slow or too shallow to compensate radiative losses, so the inner envelope behaves more like a 1D cooling atmosphere (Moldenhauer et al., 2021).

2. Low-recycling regime in tokamak scrape-off-layer physics

In the STEP geometry, low recycling is operationally realized not by imposing a low-recycling divertor target, but by combining a high-recycling target with low-recycling side walls. The relevant magnetic configuration is a spherical tokamak with a long outboard divertor leg, strong flux expansion, and a narrow plasma channel near the divertor plate. The SOL width in the divertor leg is $\sim 15\ \mathrm{mm}$, while the ionization mean free path of neutral deuterium is tens of meters, so most neutral deuterium ejected from the plate and walls passes through the plasma without undergoing a reaction (Shukla et al., 12 Nov 2025).

This geometry changes the global particle balance. Neutral deuterium can traverse the divertor leg many times and strike the side walls repeatedly before reacting with the plasma. With lithium-coated walls at $R_{\text{wall}}=0.5$, many neutrals are absorbed rather than recycled repeatedly. The result is low neutral density near the plates, weak neutral cooling, and low SOL density, even though the target itself remains high recycling with $R_{\text{target}}=0.99$ (Shukla et al., 12 Nov 2025).

The simulated $10\ \mathrm{keV}$ STEP case exhibits the canonical low-recycling signatures emphasized in the paper. Near the separatrix, the electron density is $R$0 and the electron temperature is $R$1. Near the divertor plate, the electron density is $R$2, the electron temperature is $R$3, and the ion temperature at the peak heat-flux channel exceeds $R$4. These profiles indicate a hot, attached, weakly collisional SOL with no significant cooling in the divertor leg (Shukla et al., 12 Nov 2025).

A frequent misconception is that low recycling requires a low-recycling target material. The STEP study argues the opposite for this geometry: low-recycling SOL conditions can be realized by controlling wall recycling, provided that long neutral mean free paths and a narrow divertor-leg plasma channel allow target-born neutrals to be intercepted by low-$R$5 side walls before they repopulate the SOL (Shukla et al., 12 Nov 2025).

3. Kinetic modeling, fluid closures, and the long-mean-free-path regime

The low-recycling SOL studied for STEP lies in a parameter regime where kinetic effects become decisive. The paper uses Gkeyll’s axisymmetric gyrokinetic solver coupled to EIRENE, with a Maxwellian source at the inner radial boundary of $R$6 and input power $R$7. Cross-separatrix transport is modeled with $R$8 and $R$9, and drifts are turned off in the main run used to design a $R_{\text{target}}=0.99$0 heat-flux width (Shukla et al., 12 Nov 2025).

The gyrokinetic formulation evolves the gyrocenter distribution function $R_{\text{target}}=0.99$1 with collisions, sources, electrostatic fields, and magnetic-mirror effects. In schematic form,

$R_{\text{target}}=0.99$2

Within the STEP application, the crucial kinetic ingredients are $R_{\text{target}}=0.99$3, grad-$R_{\text{target}}=0.99$4, curvature and parallel motion, the mirror force, and the parallel electric field $R_{\text{target}}=0.99$5. These are essential because the SOL is weakly collisional, neutral mean free paths are extremely long, and non-Maxwellian distributions, trapped-passing dynamics, and finite-orbit-width effects materially affect heat transport and impurity motion (Shukla et al., 12 Nov 2025).

The comparison model is SOLPS, coupling B2.5 to EIRENE. Its fluid equations include continuity, parallel momentum, and energy evolution with Braginskii-type parallel heat-flux closures such as $R_{\text{target}}=0.99$6. The paper emphasizes that this description neglects the mirror force in the parallel momentum equation, cannot represent trapped populations in velocity space, and does not reproduce collisionless banana-orbit spreading of heat flux. The assumption $R_{\text{target}}=0.99$7 for impurities also fails in the low-collisionality regime studied (Shukla et al., 12 Nov 2025).

These modeling differences translate directly into distinct predictions. In a related $R_{\text{target}}=0.99$8 SOL-only case, Gkeyll gives an upstream ion temperature of $R_{\text{target}}=0.99$9, whereas SOLPS gives $R_{\text{wall}}=0.5$0, because mirror trapping reduces heat conduction in the kinetic calculation. Along the outboard separatrix, the upstream charged-argon density in Gkeyll is 100 times lower than in SOLPS, indicating much stronger kinetic confinement of impurities to the divertor region. For heat-flux widths, both models give $R_{\text{wall}}=0.5$1 without drifts, but with drifts SOLPS remains at $R_{\text{wall}}=0.5$2 while Gkeyll broadens to $R_{\text{wall}}=0.5$3, equal to the ion banana width, with a factor-of-2 reduction in peak heat flux (Shukla et al., 12 Nov 2025).

4. Heat exhaust, impurity confinement, and material constraints

The low-recycling SOL is attractive for core confinement because it combines high edge temperature with low edge density, but the same regime intensifies heat-exhaust constraints. In the $R_{\text{wall}}=0.5$4 STEP case without drifts, the upstream heat-flux width is $R_{\text{wall}}=0.5$5, and the peak heat flux at the outboard plate is $R_{\text{wall}}=0.5$6 (Shukla et al., 12 Nov 2025).

The paper argues that kinetic physics mitigates two standard objections to such operation. First, in a spherical tokamak with large trapped fraction and substantial mirror ratio, hot low-collisionality plasmas develop strong parallel potential drops and banana-orbit spreading. The authors note that for the $R_{\text{wall}}=0.5$7 case the ion banana width is $R_{\text{wall}}=0.5$8, so a full kinetic calculation including drifts and mirror force would likely broaden $R_{\text{wall}}=0.5$9 far beyond $1\,M_\oplus$0–$1\,M_\oplus$1, thereby reducing peak heat loads (Shukla et al., 12 Nov 2025). Second, impurity confinement strengthens rather than weakens in this regime. The combination of a large parallel potential drop and impurity temperatures lower than main-ion temperatures means that only a small fraction of impurities can climb back upstream; in the cited kinetic simulations the upstream impurity density is two orders of magnitude below the fluid prediction (Shukla et al., 12 Nov 2025).

Material selection then becomes a controlling issue. Lithium is a classic low-recycling material, but the paper states that it tends to evaporate at high heat fluxes and can only go up to $1\,M_\oplus$2 without substantial evaporation. Using the simulated target conditions, a tungsten plate with a $1\,M_\oplus$3 lithium coating would reach $1\,M_\oplus$4. The estimated lithium evaporation rates are $1\,M_\oplus$5 at $1\,M_\oplus$6, $1\,M_\oplus$7 at $1\,M_\oplus$8, and $1\,M_\oplus$9 at $0.1\,\mathrm{au}$0, compared with a deuterium ion flux of $0.1\,\mathrm{au}$1 to one outboard plate (Shukla et al., 12 Nov 2025). At $0.1\,\mathrm{au}$2, the lithium evaporation is therefore two orders of magnitude larger than the deuterium ion flux. The paper further notes that even if only $0.1\,\mathrm{au}$3 of the lithium is ionized before leaving the plasma, the ionized-lithium source becomes comparable to the entire deuterium influx, strongly cooling and densifying the SOL and destroying the desired regime (Shukla et al., 12 Nov 2025).

The resulting design logic is specific: lithium is retained only on side walls, where the heat flux is low enough to avoid strong evaporation, while the target is assigned high-heat-flux-compatible materials such as tungsten or alternative liquid metals. In the paper’s interpretation, low-recycling SOL conditions are then maintained globally by the low-$0.1\,\mathrm{au}$4 walls rather than by a low-$0.1\,\mathrm{au}$5 target (Shukla et al., 12 Nov 2025).

5. Contrastive usage in protoplanetary atmosphere–disk systems

In planetary-envelope studies, “recycling” refers to continual exchange between a forming proto-atmosphere and the surrounding protoplanetary disk. The cited three-dimensional radiation-hydrodynamic simulations solve the equations of mass continuity, momentum, and gas energy with radiative transfer via flux-limited diffusion in a local frame centered on a $0.1\,\mathrm{au}$6 core at $0.1\,\mathrm{au}$7, embedded in an MMEN-type disk with midplane density $0.1\,\mathrm{au}$8, temperature $0.1\,\mathrm{au}$9, and constant opacity $\sim 15\ \mathrm{mm}$0 (Moldenhauer et al., 2021).

The simulated flow topology is open rather than closed. After time averaging over small oscillations, the velocity field shows polar inflow and midplane outflow, with no permanently bound “closed bubble” enclosing the core. Recycling time $\sim 15\ \mathrm{mm}$1 is defined by backward integration of dense tracer populations along steady-state streamlines until each tracer exits the Hill sphere. In the fiducial close-in case, recycling is shortest in the outer Hill sphere and near the poles, longer in the dense inner region, but still finite even down to the core surface; only a narrow ring at the interface between polar inflow and midplane outflow fails to recycle within $\sim 15\ \mathrm{mm}$2, and that region contains negligible mass (Moldenhauer et al., 2021).

Quantitatively, the bulk of the Hill-sphere mass is recycled on $\sim 15\ \mathrm{mm}$3, inner atmospheric regions on $\sim 15\ \mathrm{mm}$4, and more than $\sim 15\ \mathrm{mm}$5 of the mass inside the Hill sphere is exchanged with disk gas within $\sim 15\ \mathrm{mm}$6. The spherically averaged entropy initially drops by radiative cooling, but after $\sim 15\ \mathrm{mm}$7 the entropy, temperature, and density profiles become time independent up to tiny oscillations; fluctuations at fixed radii remain $\sim 15\ \mathrm{mm}$8 without secular drift. The total gas mass inside the Hill sphere converges to

$\sim 15\ \mathrm{mm}$9

In this steady state, radiative cooling is exactly compensated by inward advective transport of enthalpy, summarized as

$R_{\text{wall}}=0.5$0

Under the same core mass, opacity, and disk conditions, the 1D model instead continues cooling and contracting, reaching $R_{\text{wall}}=0.5$1 at $R_{\text{wall}}=0.5$2 and $R_{\text{wall}}=0.5$3 at $R_{\text{wall}}=0.5$4, indicating runaway gas accretion within $R_{\text{wall}}=0.5$5–$R_{\text{wall}}=0.5$6 orbital periods (Moldenhauer et al., 2021).

Within that framework, a low-recycling regime is an inferred contrast rather than a directly simulated state. The source material suggests that such a regime would arise when exchange with the disk no longer penetrates deeply enough or rapidly enough to balance radiative losses. The implied layer-by-layer criterion is

$R_{\text{wall}}=0.5$7

versus

$R_{\text{wall}}=0.5$8

This suggests a thermally isolated inner region that cools and contracts quasi-1D, while only outer layers remain dynamically connected to the disk. The discussion in the paper further indicates that recycling should weaken at larger orbital radii because the disk impact rate

$R_{\text{wall}}=0.5$9

decreases outward, while gravitational binding at the planetary surface steepens relative to the disk environment (Moldenhauer et al., 2021). A plausible implication is that the efficient close-in regime naturally favors mini-Neptune-like envelopes, whereas low recycling at larger radii or higher core masses would favor bound inner envelopes and eventual giant-planet growth.

6. Conceptual synthesis, misconceptions, and open problems

Across these two domains, low recycling is not merely “little exchange” in a generic sense; it denotes a specific timescale hierarchy in an open system. In the STEP SOL, low recycling means that neutrals generated at material boundaries do not efficiently return to the plasma, yielding high $R_{\text{target}}=0.99$0, low $R_{\text{target}}=0.99$1, and a kinetic, long-mean-free-path edge (Shukla et al., 12 Nov 2025). In the protoplanetary-envelope context, the inferred low-recycling limit is the opposite hierarchy: gas parcels remain in the envelope long enough to cool, so radiative contraction is no longer arrested by exchange with high-entropy surroundings (Moldenhauer et al., 2021).

The two literatures also correct different misconceptions. In SOL physics, low recycling does not require a low-recycling target; the STEP study is explicitly organized around a high-recycling target and low-recycling walls, with geometry and neutral mean free paths determining the global outcome (Shukla et al., 12 Nov 2025). In planetary atmospheres, the converse point is that an apparently “bound” inner envelope need not be dynamically isolated: in the simulated close-in case, even gas down to the core surface recycles on finite timescales, and essentially the entire Hill-sphere mass participates in the steady state (Moldenhauer et al., 2021).

Open problems remain domain specific. For STEP-like low-recycling SOLs, the paper identifies helium exhaust, full drift-included kinetic modeling at $R_{\text{target}}=0.99$2, more realistic wall and liquid-metal surface physics, and 3D or turbulent effects as unresolved directions (Shukla et al., 12 Nov 2025). For protoplanetary envelopes, the outstanding issue is parameter mapping: the paper explicitly points to future work on planet mass, orbital location, and disk headwind as the route to determining where the efficient recycling steady state gives way to the inferred low-recycling regime (Moldenhauer et al., 2021).

In both cases, the low-recycling regime is best understood not as an isolated material property, but as a systems-level state determined by transport topology, boundary interaction, and the competition between exchange times and internal relaxation times.