---
title: Atmosphere-Breathing Electric Propulsion (ABEP)
url: https://www.emergentmind.com/topics/atmosphere-breathing-electric-propulsion-abep
type: topic
---

# Atmosphere-Breathing Electric Propulsion (ABEP)

Atmosphere-Breathing Electric Propulsion (ABEP) is an electric-propulsion architecture in which a spacecraft intake collects residual atmospheric particles encountered along the ram direction and delivers them to a thruster, so that the ambient atmosphere itself serves as propellant for continuous drag compensation in very low orbital regimes. In the literature surveyed here, ABEP is primarily formulated for Very Low Earth Orbit (VLEO, below about \(450\) km) and, by extension, for very low orbits around other atmospheric bodies. The system is usually decomposed into an intake, a compression or flow-conditioning stage, an ionization stage, and an acceleration stage, with feasibility governed by the coupled balance among rarefied-flow capture, thruster efficiency, and spacecraft drag [2106.15912], [2211.09493], [2607.02635].

## 1. Concept and orbital regime

ABEP is motivated by the fact that at very low altitudes the residual atmosphere generates substantial drag, while conventional electric propulsion remains constrained by finite onboard propellant. In the VLEO regime, the flow is free molecular or highly rarefied, with \(Kn > 10\) emphasized in recent gas–surface interaction work, so the dominant physics is not continuum compression but wall-mediated transport, scattering, and backflow. Atmospheric particles strike spacecraft surfaces at hypersonic orbital speeds of \(6{,}000\)–\(10{,}000\) m/s; representative values cited across the ABEP literature are \(u_\infty \approx 7.8\) km/s, \(v \approx 7.5\)–\(7.8\) km/s, and \(T_\infty\) in the several-hundred-Kelvin range [2606.31928], [2504.12829].

The ambient composition is altitude dependent. For Earth VLEO, the cited models and simulations include atomic oxygen \(O\), \(N_2\), \(O_2\), He, Ar, H, and N, with \(O\) and \(N_2\) usually dominating mass flow; recent intake and scattering studies also treat \(O_2\) formation through surface recombination and explicitly note the erosive role of atomic oxygen in materials selection [2106.15912], [2504.12829]. The same ABEP logic is extended in mission studies to Mars, where very low Mars orbit is analyzed with \(CO_2\)-dominated atmospheres [2103.02328], [2211.09493].

At system level, the defining balance is between aerodynamic drag and propulsive thrust. The recurring relations are

\[
D = \frac{1}{2}\,\rho(h)\,v(h)^2\,C_D\,A_f
\]

and

\[
T \approx \dot{m}\,v_e,
\]

with stationkeeping requiring \(T \ge D\) [2103.02328], [2211.09493], [2607.02635]. Because \(D\) scales with density and frontal area, and \(\dot{m}\) scales with intake capture, ABEP feasibility is intrinsically coupled to intake geometry, gas–surface interaction, and available power rather than to thruster performance alone. Earth-focused mission analyses place a practical lower bound near \(120\) km because of heating and re-entry onset, while upper ABEP design envelopes such as \(150\)–\(250\) km are chosen because at higher altitude collectible mass flow becomes marginal for typical spacecraft [2103.02328], [2211.10079].

## 2. Intake physics, capture, and gas–surface interaction

The intake is the defining ABEP subsystem because it must collect a highly collimated, hyperthermal beam and deliver it to the thruster under free-molecular conditions. The standard intake performance measure is the collection efficiency

\[
\eta_c(h) = \frac{\sum_{i=1}^{N_s}\dot{N}_{\mathrm{out},i}(h)}{\sum_{i=1}^{N_s}\dot{N}_{\mathrm{in},i}(h)},
\]

with the delivered mass flow

\[
\dot{m}_{thr}(h) = \sum_{i=1}^{N_s} m_{p,i}\,\dot{N}_{\mathrm{out},i}(h).
\]

A closely related free-stream estimate is \(\dot{m}_{\mathrm{col}} = \rho_\infty V_\infty A_{\mathrm{in}}\eta_c\) [2106.15912], [2103.02328], [2211.10079].

The reviewed intake literature repeatedly treats gas–surface interaction through idealized limiting models. In the PICLas intake-design study, Maxwell boundary conditions are implemented with `MomentumACC = 0.0` for specular reflection and `MomentumACC = 1.0` for diffuse reflection [2106.15912]. In a separate intake-thruster design paper, energy exchange is described through

\[
\alpha = \frac{E - E'}{E - E_W},
\]

with the paper stating that specular reflection corresponds to \(\alpha = 1\) and diffuse reflection to \(\alpha = 0\) [2211.10079]. The reviewed papers therefore do not use a single accommodation-coefficient convention, but they agree on the physical distinction: specular walls preserve directionality and energy, whereas diffuse walls thermalize and randomize the reflected flux.

That distinction leads directly to the observed geometry classes. Diffuse intakes in the DISCOVERER studies use hexagonal honeycomb front ducts followed by a conical chamber; the ducts act as molecular traps, larger central ducts increase direct transmission, and narrower outer ducts suppress perimeter escape. Under these assumptions, the D-MDAR concept at \(h=150\) km yields \(A_{\mathrm{intake}}=0.008\) m\(^2\), \(A_{in}=0.004\) m\(^2\), \(\dot{N}_{in}=1.34\times 10^{18}\) s\(^{-1}\), \(\dot{N}_{out}=6.14\times 10^{17}\) s\(^{-1}\), \(\eta_c=0.458\), \(\dot{m}_{thr}=0.0240\) mg/s, and \(p_{ch}\approx 0.27\) Pa, with \(\eta_c\) decreasing to \(0.416\) by \(250\) km [2106.15912]. Specular intakes instead use parabolic scoops that optically focus incident trajectories into the discharge channel; the S-FT design at \(150\) km reaches \(A_{in}=0.019\) m\(^2\), \(\dot{N}_{in}=6.33\times 10^{18}\) s\(^{-1}\), \(\dot{N}_{out}=5.97\times 10^{18}\) s\(^{-1}\), \(\eta_c=0.943\), \(\dot{m}_{thr}=0.232\) mg/s, and \(p_{ch}\approx 0.3\) Pa [2106.15912].

Misalignment sensitivity is also largely a gas–surface problem. For the diffuse D-MDAR intake at \(150\) km, \(\eta_c\) falls from \(0.456\) at \(\alpha=0^\circ\) to \(0.378\), \(0.270\), \(0.194\), and \(0.150\) at \(5^\circ\), \(10^\circ\), \(15^\circ\), and \(20^\circ\), respectively [2106.15912]. The specular S-FT intake maintains nearly constant \(\eta_c\) up to about \(10^\circ\), then drops sharply; by about \(20^\circ\), \(\eta_c \approx 0.13\) [2106.15912]. The literature explicitly notes that real surfaces deviate from ideal specular and diffuse limits because roughness, contamination, and atomic-oxygen erosion drive reflection behavior over time [2106.15912].

A further refinement appears in recent molecular-dynamics-based work on gas–surface scattering for ABEP intakes. There the central object is the scattering kernel \(K\), defined in the local frame by

\[
f_{out}(v') = \int K(v \rightarrow v')\,f_{in}(v)\,dv,
\]

with normalization, positivity, and detailed-balance constraints imposed so that multi-bounce deceleration can be modeled from thermal to hypersonic velocities [2606.31928]. This is particularly important because ABEP intakes rely on successive wall collisions to decelerate particles from orbital speed to near-thermal velocity before they reach the plasma stage.

## 3. Thruster architectures and intake–thruster coupling

The thruster side of ABEP is dominated by electrodeless RF concepts because chemically aggressive species, especially atomic oxygen, erode electrodes, grids, and discharge channels in conventional electric propulsion. The RF Helicon-based Inductive Plasma Thruster (IPT) developed at IRS within the DISCOVERER project is explicitly designed as a contactless, neutralizer-free device: an intake feeds neutrals into a discharge channel, an RF antenna and axial magnetic field sustain a helicon discharge, and the exhaust is a quasi-neutral plasma plume [2007.06397], [2211.10079], [2607.02635].

The laboratory IPT reported by IRS uses a birdcage antenna resonant at \(40.68\) MHz, an external solenoid capable of \(B_0 < 70\) mT, a vacuum chamber with base pressure \(p < 0.01\) Pa, and operation on Ar, \(N_2\), and \(O_2\). The design is centered on resonance and low reflected power; in the 2026 dissertation, the finalized birdcage has measured \(S_{11} \approx -24.8\) dB at \(40.68\) MHz with \(Z_{IPT} \approx 44.56 + j0.17\ \Omega\), implying \(>99\%\) coupling in vacuum, and post-test measurements improved to \(S_{11} = -43.4\) dB with \(Z_{IPT} \approx 49.84 - j0.66\ \Omega\) [2607.02635]. Experimentally, successful ignition and operation across relevant mass flows were demonstrated with low power consumption, reported as \(P \sim 60\) W or \(P_f < 60\) W depending on the study [2211.10079], [2607.02635].

A second IRS line of work uses the inductively heated IPG6-S as an ABEP-relevant test-bed for an Inductive Plasma Thruster. In that configuration the RF source operates at about \(4\) MHz with active RF power up to about \(3.5\) kW, using \(O_2\) and air as working gases. The measured electric-to-thermal coupling efficiency is reported up to about \(30\%\) for \(O_2\) and about \(25\%\) for air, with a maximum specific enthalpy \(h_{cal} \approx 7.5\) MJ/kg at \(\dot{m}_{gas} \approx 60\) mg/s. Under the paper’s upper-limit assumption that the measured enthalpy is fully converted to directed kinetic energy,

\[
v_e \approx \sqrt{2\,h_{cal}},
\]

which gives \(v_e \approx 3.9\) km/s and an upper-limit thrust of about \(230\)–\(250\) mN at low altitudes [2103.02328].

A distinct particle-based simulation of a CubeSat-scale ABEP system resolves the intake and thruster jointly with DSMC for neutrals and PIC-MCC for the plasma. In that model, the RF coil is driven at \(13.56\) MHz, the screen grid is RF biased with \(V_{pp}=1500\) V and \(V_{DC}=730\) V, the absorbed RF power stabilizes near \(20\) W, and the steady coil current is about \(1.83\) A. The simulated ion exhaust speed for \(N_2^+\) under \(\Delta V \approx 730\)–\(1480\) V is \(71\)–\(100\) km/s, corresponding to \(I_{sp} \approx 7200\)–\(10{,}200\) s [2504.12829]. That study separates the architecture into intake, ionization stage, acceleration grids, and neutral exhaust, while explicitly noting that one-way coupling neglects neutral depletion and back-pressure feedbacks.

A related but distinct electrodeless RF concept is the SUb-atmospheric Radio-frequency Engine (SURE), designed for near-space pressures of \(32\)–\(5332\) Pa rather than VLEO free-molecular conditions. SURE uses two axially separated solenoid coils at \(13.56\) MHz and up to \(1\) kW, combining inductive coupling under each leg with capacitive coupling between the legs. It reports plasma density up to \(2.23\times 10^{18}\) m\(^{-3}\) and electron temperature up to \(2.79\) eV without plasma-facing electrodes or external magnetic field [2006.09842]. Although its target regime differs from orbital ABEP, it belongs to the broader class of air-breathing electrodeless RF propulsion schemes.

## 4. Modeling hierarchy: DSMC, PIC, and learned scattering kernels

ABEP has been developed through a layered modeling hierarchy because no single solver spans free-molecular intake transport, gas–surface scattering, plasma production, and mission-level drag compensation. On the neutral side, DSMC is the dominant intake tool. PICLas is used in the DISCOVERER intake-design work for free-molecular or highly rarefied intake transport with Maxwell wall models [2106.15912]. Pantera is used for a one-way coupled DSMC\(\rightarrow\)PIC-MCC workflow in a full 3D CubeSat-sized ABEP geometry, where DSMC delivers the neutral density to the plasma solver on the same unstructured mesh [2504.12829]. A hybrid PIC–DSMC framework has also been used to analyze ground testing of ABEP intakes under electric-propulsion plumes, including ion-surface neutralization and sputtering [2509.17547].

Recent work extends this hierarchy by replacing idealized wall models with a data-driven scattering kernel learned from molecular dynamics. The 2026 scattering-kernel paper parameterizes the conditional density \(p(v'|v)\) with a conditional Real-valued Non-Volume Preserving flow, cRealNVP, and enforces positivity and normalization by construction. The kernel is trained on molecular-dynamics data for atomic oxygen on amorphous \(Al_2O_3\), with incident velocities ranging from thermal speeds to \(10{,}000\) m/s and nine polar angles \(\theta_i \in [0^\circ,80^\circ]\) [2606.31928].

The thermodynamic constraint central to that work is detailed balance. In its probability form,

\[
p_{eq}(v)\,K(v \rightarrow v') = p_{eq}(v')\,K(v' \rightarrow v),
\]

and in the flux-weighted form,

\[
|v \cdot n|\,p_{eq}(v)\,K(v \rightarrow v')
=
|v' \cdot n|\,p_{eq}(v')\,K(-v' \rightarrow -v),
\]

with

\[
p_{eq}(v)=\left(\frac{m}{2\pi k_B T}\right)^{3/2}\exp\!\left(-\frac{m\|v\|^2}{2k_B T}\right).
\]

The total loss is

\[
L_{total}=L_{NLL}+\lambda_{DB}L_{DB},
\]

where \(L_{DB}\) penalizes deviations from detailed balance on equilibrium samples [2606.31928]. This matters operationally because ABEP intake particles undergo multiple collisions and eventually enter near-thermal regimes where equilibrium consistency is no longer optional.

Validation results are explicit. When a Maxwell flux at \(T_w=300\) K is scattered with cRealNVP, reservoir simulation starting from \(T_0=700\) K converges to a mean of about \(300.29\) K over the last \(3000\) steps, corresponding to about \(0.1\%\) error. The prior CVAE-based kernel stabilizes at about \(454.19\) K, which the paper presents as evidence of equilibrium inconsistency [2606.31928]. The cRealNVP model also matches molecular-dynamics reflected distributions across incident magnitudes \(3000\), \(7261.3\), and \(8585.9\) m/s and extrapolates successfully to \(11{,}000\) m/s, while flat-plate DSMC tests show that both cRealNVP and CVAE reproduce MD-derived lift, drag, and \(L/D\) trends better than a fixed-accommodation Maxwell model [2606.31928].

Because the learned \(K\) is explicit over the full speed range, it can be embedded directly in multi-bounce intake simulation. The reported workflow initializes particle state at the intake entrance, propagates ballistically to the next wall, samples outgoing velocity from \(K(v\rightarrow v';c)\), repeats over successive wall interactions, and terminates when the particle is captured through the throat, lost to free stream, sticks if that is modeled, or falls below a thermal threshold [2606.31928]. This formulation connects molecular-scale scattering directly to capture efficiency, compression, and the phase-space distribution entering the ionization stage.

## 5. Performance envelopes, active intakes, and mission applications

ABEP performance is strongly assumption dependent because altitude thresholds move with \(A_f\), \(A_{in}\), \(\eta_c\), \(\eta_T\), and available electrical power. In the passive-intake DISCOVERER study, the mass flow and required exhaust velocity already show a large specular–diffuse separation. At \(150\) km, the S-FT specular intake provides \(\dot{m}_{thr} \approx 2.32\times 10^{-7}\) kg/s, while D-MDAR gives \(\dot{m}_{thr} \approx 2.40\times 10^{-8}\) kg/s. For full drag compensation with the assumed frontal areas, the required exhaust velocity is about \(9.1\)–\(9.2\) km/s for S-FT over \(150\)–\(250\) km, but about \(35.8\)–\(39.8\) km/s for D-MDAR [2106.15912]. This is why the same paper concludes that the specular parabolic scoop with focus inside the discharge channel is the most promising intake in its design set.

Mission-level Earth analyses confirm that conclusion under broader spacecraft assumptions. In the Earth circular-orbit study with \(A_f=A_{in}=1\) m\(^2\), \(C_D=3.7\), and \(\eta_T=0.20\), a \(200\) km orbit requires \(v_e=33.49\) km/s and peak power \(P_{req,peak}\approx 2.87\) kW for \(\eta_c=0.43\), whereas \(\eta_c=0.70\) reduces these to \(20.57\) km/s and about \(1.75\) kW [2211.09493]. The same study finds that full drag compensation at \(250\) km is feasible with \(P_{GOCE}=1.6\) kW for \(\eta_t=0.20\) and \(\eta_c\) in the range \(0.43\)–\(0.70\), while at \(180\) km feasibility at that power level requires reducing \(A_f\) to about \(0.465\) m\(^2\) [2211.09493].

A complementary GOCE case study in the 2026 dissertation uses actual intake areas rather than \(A_f=A_{in}\). With \(A_f=1.1\) m\(^2\), \(C_D \approx 3.7\), \(\eta_T=0.2\), and clustered intakes, it reports \(D \approx 6\)–\(200\) mN and \(\dot{m}_{thr} \approx 0.2\)–\(20\) mg/s across \(150\)–\(250\) km, with required exhaust velocity \(c_e \approx 15\)–\(32\) km/s. Under those assumptions, continuous operation within \(1.6\) kW is feasible above about \(215\) km for diffuse or EFD intakes and above about \(190\) km for the specular intake, yielding \(T/P \approx 10\)–\(24\) mN/kW and \(I_{sp} \approx 1500\)–\(3000\) s [2607.02635].

Not all ABEP studies are positive under compact-power assumptions. In the 3U CubeSat particle-based simulation, the intake compresses neutrals in the ionization chamber to \(n_{N_2} \approx 1.0\times 10^{20}\) m\(^{-3}\), \(n_{O_2} \approx 1.9\times 10^{19}\) m\(^{-3}\), and \(n_O \approx 2.0\times 10^{16}\) m\(^{-3}\), but a power-limited estimate for a \(\sim 20\) W thruster gives \(T \approx 0.15\)–\(0.25\) mN, compared with a drag estimate \(F_d \approx 9.3\) mN for a \(10\times 10\) cm frontal area [2504.12829]. That study therefore concludes that the modeled CubeSat-scale system cannot station-keep near \(200\) km unless the frontal area is reduced, the intake capture is increased, the power is raised to the \(100\)–\(200\) W class, or a lower-\(v_e\) operating point is adopted [2504.12829]. The literature thus separates low-power plasma operation from full system drag compensation.

A qualitatively different path is the cryocondensation-regeneration active intake device, CRAID. Instead of relying only on passive multi-bounce transport, CRAID stores propellant as condensed film on a cryopanel at \(20\) K, then sublimates it at \(54.5\) K and meters it to the thrusters through narrow injection tubes [2503.02021]. At \(200\) km, the conceptual prototype model gives \(t_{cond}^\ast = 42.4\) min, \(t_{regen}^\ast = 31.7\) min, \(t_{cool}^\ast = 1.2\) min, and \(t_{total}^\ast = 75.3\) min, with effective cycle-averaged \(\eta_c^{eff}=26.1\%\) and \(CR^{eff} \approx 3.0\times 10^7\) [2503.02021]. The paper states that this compression is at least \(1000\) times higher than that of prevalent passive intake devices and reports complete drag compensation for \(h \ge 190\) km under normal solar activity, with no upper altitude bound identified in the \(150\)–\(300\) km range analyzed [2503.02021].

Mission applications consequently extend beyond pure stationkeeping. Earth and Mars analyses include circular and elliptical orbit maintenance, orbit raising, constant-rate de-orbit, and atmospheric-propellant collection for refueling or space-tug missions [2211.09493]. On Mars, the same study reports full drag compensation near \(133\) km with \(\eta_c=0.43\) and about \(1\) kW, or near \(129\) km with \(\eta_c=0.70\) and about \(1\) kW [2211.09493].

## 6. Ground testing, uncertainties, and open directions

Ground validation is unusually difficult for ABEP because the target environment combines hyperthermal directed flow, extreme rarefaction, and reactive species. The Rarefied Orbital Aerodynamics Research (ROAR) facility is designed around neutral atomic oxygen produced by electron-stimulated desorption, with a baseline DSMC inflow of \(V_\infty = 7340\) m/s, \(T_\infty = 523\) K, \(\Gamma_\infty = 1.028\times 10^{18}\) m\(^{-2}\) s\(^{-1}\), and \(n_\infty = 1.4\times 10^{14}\) m\(^{-3}\) [2406.06299]. Two intake-test methodologies are proposed there: a pressure-difference method using the rise of \(p_{aux}(t)\) in an auxiliary chamber, and a gas-sensor method using a quartz crystal microbalance at the outlet [2406.06299].

For the pressure-based method, the paper derives

\[
p_{aux}(t)=
\frac{\sqrt{2\pi m k_B T_{aux}}}{\eta_{rpc}A_{aux}}\,N_{in}\left(1-e^{-t/\tau}\right),
\]

with

\[
\tau=
\frac{V_{aux}}{\eta_{rpc}A_{aux}}
\sqrt{\frac{2\pi m}{k_B T_{aux}}},
\]

so that fitting the transient directly yields \(\eta_c\) once the reverse-performance coefficient \(\eta_{rpc}\) is known [2406.06299]. In the same study, a perfectly specular parabolic inlet gives \(\eta_c > 90\%\) at \(\alpha \approx 0^\circ\), while among the tested subscale geometries the Open channel with \(L_c=5\) cm and \(D_c=10\) mm reaches \(\eta_c=76.6\%\) and \(\eta_{rpc}=6.8\%\) [2406.06299].

A separate numerical analysis studies intake testing with an electric-propulsion plasma plume rather than a neutral AO beam. That hybrid PIC–DSMC model shows that the outlet flow can become neutral-dominated because ions are neutralized on the intake surface, and that sputtering becomes non-negligible at high beam energies. The central experimental conclusion is that simultaneous ion and neutral diagnostics are required for reliable capture-efficiency evaluation when EP plumes are used as a flow generator [2509.17547]. This is an important correction to any testing approach that would infer intake performance from ion current alone.

The principal uncertainties identified across the literature are not merely engineering details. In the learned-kernel work, the molecular-dynamics database is trained primarily for atomic oxygen on amorphous \(Al_2O_3\); chemistry, adsorption, contamination, and anisotropic roughness are omitted, and temperature dependence at high velocity is approximated through scaling [2606.31928]. In passive-intake studies, real surfaces are acknowledged to deviate from ideal specular or diffuse behavior under AO exposure, contamination, and microroughness growth [2106.15912]. In coupled neutral–plasma simulations, one-way coupling neglects neutral depletion by ionization and back-pressure from the thruster [2504.12829]. In active-intake work, CRAID introduces cryocoolers, valves, and mass-flow controllers, so valve reliability, AO-mediated surface chemistry, and long-duration cryosurface contamination become part of the propulsion problem [2503.02021].

Future directions are correspondingly multi-scale. The gas–surface community names multi-species reactive kernels for \(AO\), \(N_2\), and \(O_2\), explicit conditioning on wall temperature, material state, and anisotropic roughness, uncertainty quantification, and coupling to plasma solvers and optimization loops [2606.31928]. System simulations call for two-way neutral–plasma coupling, more complete multi-species chemistry, refined wall models, and integrated power-system modeling [2504.12829]. Experimental programs prioritize AO wind-tunnel campaigns in ROAR, calibration against ROAR data and SOAR in-orbit measurements, and end-to-end intake–thruster validation under representative rarefied conditions [2106.15912], [2406.06299].

Taken together, the current ABEP literature defines a field that is no longer limited to the original passive-intake concept. It now includes high-\(\eta_c\) specular parabolic intakes, diffuse molecular-trap geometries, active cryogenic intakes with cycle-averaged compression ratios near \(3.0\times 10^7\), electrodeless RF helicon thrusters with laboratory operation near \(60\) W, system analyses for Earth and Mars, and thermodynamically constrained learned scattering kernels that close a longstanding gap between molecular dynamics and intake-level rarefied-flow prediction [2106.15912], [2503.02021], [2607.02635], [2606.31928].

Source: https://www.emergentmind.com/topics/atmosphere-breathing-electric-propulsion-abep