Magnetoplasmadynamic Thruster

A Magnetoplasmadynamic Thruster (MPD thruster or MPDT, also historically termed a Lorentz Force Accelerator) is a high-power electromagnetic spacecraft propulsion system that ionizes a gaseous or liquid-metal propellant and accelerates the resulting quasi-neutral plasma to exhaust velocities of 15 to 80 km s⁻¹ via the electromagnetic Lorentz body force f = J × B. Unlike electrostatic ion thrusters, whose thrust density is limited by space-charge saturation (the Child–Langmuir law) between extractor grids, or Hall-effect thrusters, which operate at moderate power densities, an MPD thruster drives a high-current radial electric arc (10² to 10⁵ A) between a central refractory cathode and a coaxial annular anode. Operating in either self-field mode (where the discharge current itself induces a purely azimuthal magnetic field B_θ) or applied-field mode (augmented by an external solenoid), MPD thrusters combine high specific impulse (I_sp ≈ 1,500–8,000 s) with thrust densities orders of magnitude greater than electrostatic engines, positioning them as primary candidates for megawatt-class nuclear-electric and solar-electric deep-space cargo and crewed missions.

An MPD thruster consists of a coaxial discharge chamber featuring a central cylindrical rod or hollow cathode (typically thoriated or lanthanated tungsten, or barium-calcium-aluminate impregnated tungsten) surrounded by a concentric cylindrical or flared anode nozzle. Propellant—such as argon, xenon, hydrogen (H₂), nitrogen (N₂), hydrazine decomposition products, or vaporized lithium (Li)—is injected at the upstream backplate and ionized into a highly conducting plasma by a high-current radial discharge I_d flowing inward from the anode to the central cathode.

Self-Field Lorentz Acceleration and the Maecker Formula

In a self-field MPD thruster (SF-MPDT), no external magnet is required. By Ampère's circuital law, the axial and radial current density J = (J_r, 0, J_z) flowing onto the central cathode of radius r_c generates an axisymmetric azimuthal magnetic field B_θ(r, z) = −μ₀I_enc(r, z) / (2*πr*). Crossing the current density J with the self-induced azimuthal field B = (0, B_θ, 0) yields two complementary Lorentz body-force components (f = J × B):

1. Axial Blowing Force (f_z = −J_r B_θ): The radial current J_r crossing the azimuthal field B_θ pushes the plasma directly downstream toward the exhaust plane (analogous to the magnetic piston in a coaxial railgun). Integrating f_z over an annular discharge volume between cathode radius r_c and anode radius r_a yields Maecker's (1955) electromagnetic thrust formula:

F_blow = (μ₀ I_d² / 4π) ln(r_a / r_c)

1. Radial Pinch / Pumping Force (f_r = J_z B_θ): Where current lines curve axially along the cathode (J_z ≠ 0), the radial inward pinch force compresses the plasma against the cathode surface and tip, creating a high-pressure stagnation spike on the cathode face that converts into axial kinetic momentum via pressure balances:

F_EM = (μ₀ I_d² / 4π) [ ln(r_a / r_c) + k_tip ]

where k_tip ∈ [0, 3/4] is a dimensionless cathode-tip geometry constant (equal to 3/4 when current attaches uniformly across a flat hemispherical cathode face). Notice that the electromagnetic thrust F_EM = b I_d² (with geometric inductance gradient b = (μ₀ / 4π)[ln(r_a / r_c) + k_tip] ≈ 10⁻⁷ N A⁻²) depends purely on the square of the discharge current I_d² and the electrode radius ratio r_a / r_c, independent of propellant atomic mass or mass flow rate ṁ.

Applied-Field MPD Thrusters (AF-MPDT)

Because self-field operation requires kiloampere currents (I_d ≳ 5–20 kA, corresponding to megawatt instantaneous powers) for B_θ to exceed a few tenths of a tesla, steady-state thrusters operating below 200 kW incorporate an external electromagnetic solenoid or permanent magnet surrounding the anode to impose an axial and radial applied magnetic field B_app = (B_r,app, 0, B_z,app) of 0.05 to 1.0 T. In an Applied-Field MPD Thruster, four acceleration mechanisms act in concert:

• Swirl (Azimuthal) Acceleration: The radial discharge current J_r crossed with the applied axial field B_z,app exerts an azimuthal torque f_θ = −J_rB_z,app, spinning the plasma column to high rotational speeds; subsequent expansion through a diverging magnetic nozzle converts azimuthal rotational energy and thermal enthalpy into axial thrust.
• Hall Acceleration: Electron azimuthal Hall currents J_θ = σ_eω_ceτ_e(E_r + u_zB_r) crossing the diverging radial magnetic field B_r,app generate direct axial thrust f_z,Hall = J_θB_r,app.
• Self-Field Blowing/Pumping and Electrothermal Gasdynamic Expansion (Joule heating recovery in the magnetic nozzle).

For a given mass flow rate ṁ, increasing the discharge current I_d increases the specific impulse I_sp = F / (ṁg₀) ≈ b I_d² / (ṁg₀). In 1964, Rudolf Hügel and collaborators discovered that every MPD thruster exhibits a fundamental stability boundary governed by the dimensional ratio J_d² / ṁ (or equivalently by the ratio of exhaust velocity u_ex to Hannes Alfvén's Critical Ionization Velocity, u_civ = √(2*eε*_i / m_i), where ε_i is the first ionization potential and m_i is the ion mass).

When I_d² / ṁ exceeds a threshold (I_d² / ṁ)_onset, the thruster enters the "Onset" regime, characterized by:

1. Anode Starvation: The inward radial Lorentz pumping force (J_zB_θ) pinches plasma toward the central axis so strongly that the region directly adjacent to the anode surface becomes depleted of charge carriers (electrons and ions).
2. Anode Fall Voltage Spike and Lower-Hybrid Drift Instabilities: To sustain the imposed current across the carrier-depleted anode sheath, the anode fall voltage V_anode spikes from ~10–20 V to >50–100 V, exciting high-frequency (MHz) electroacoustic and modified two-stream plasma micro-instabilities (Choueiri 1998).
3. Accelerated Anode and Cathode Erosion: High-energy electron bombardment and arc-spot filamentation melt and vaporize the anode lips and erode the tungsten cathode at rates incompatible with multi-thousand-hour deep-space lifetimes.

To suppress the onset instability and reduce frozen-flow ionization losses (the fraction of electrical input power wasted ionizing the propellant rather than accelerating it), researchers in Russia (Moscow Aviation Institute and RKK Energiya) and the United States (Princeton University and NASA JPL) developed the Lithium Lorentz Force Accelerator (Li-LFA) utilizing multi-channel hollow tungsten rods fed with liquid lithium.

Lithium offers three fundamental thermophysical advantages over noble gases or molecular propellants:

• Low First Ionization Potential vs. High Second Ionization Potential: Lithium's first ionization energy is only ε_i,1 = 5.39 eV (compared to 15.76 eV for argon and 13.60 eV for hydrogen), ensuring nearly 100% single ionization with minimal frozen-flow loss, while its second ionization potential is exceptionally high (ε_i,2 = 75.64 eV), preventing energy waste into doubly ionized Li²⁺ states.
• Cathode Work-Function Lowering: Adsorbed monolayer films of barium or lithium on multi-channel tungsten cathodes lower the thermionic work function from 4.5 eV to ~2.1 eV, enabling diffuse thermionic emission at lower surface temperatures (~1,500–1,800 K) and reducing cathode evaporation by two orders of magnitude.
• Condensable Exhaust in Ground Vacuum Chambers: Unlike hydrogen or argon, which require immense cryopumping speeds at megawatt flow rates, lithium condenses as a solid on cooled vacuum-chamber walls.

On orbit, pulsed quasi-steady self-field MPD thrusters have been demonstrated in space by Japan's ISAS/JAXA aboard the MS-T4 (Tansei-4) satellite in 1980, the Space Shuttle Spacelab-1 SEPAC experiment in 1983, and the Space Flyer Unit (SFU) Electric Propulsion Experiment (EPEX) in 1995 (which completed over 40,000 hydrazine MPD pulses in low Earth orbit). In 2021, the University of Stuttgart's 100-kW-class steady-state applied-field SX3 thruster demonstrated 62% thrust efficiency at I_sp = 4,665 s and 2.75 N thrust on argon, complementing alternative high-power architectures such as the Variable Specific Impulse Magnetoplasma Rocket and pulsed inductive thrusters.

Why can an MPD thruster produce much higher thrust per unit area than an electrostatic ion thruster?

In an electrostatic ion thruster, electrons are stripped out of the acceleration zone, leaving a net positive ion cloud between the grids. Electrostatic repulsion among those positive ions (the Child–Langmuir space-charge limit) restricts maximum thrust density to roughly 1 to 10 N m⁻². In an MPD thruster, the plasma remains quasi-neutral (n_e ≈ n_i) throughout the electromagnetic acceleration channel, eliminating space-charge saturation and permitting thrust densities exceeding 10⁴ to 10⁵ N m⁻².

Why are pulsed quasi-steady MPD thrusters used on kilowatt-class solar-powered spacecraft?

Self-field MPD thrusters only reach high electromagnetic efficiency when the discharge current exceeds several kiloamperes (corresponding to instantaneous electrical powers of 0.5 to 5 MW), whereas solar arrays on standard satellites provide only 1 to 20 kW continuously. By charging a capacitor-inductor pulse-forming network over ~0.5–1 second at low average power (1–5 kW) and discharging it through the thruster in a 1-millisecond, 10-kA rectangular pulse, a quasi-steady pulsed MPD thruster achieves megawatt-level Lorentz acceleration during each pulse while matching available spacecraft solar power.

What is the difference between a pulsed plasma thruster (PPT) and an MPD thruster?

Both accelerate plasma using the J × B Lorentz force, but a classic Pulsed Plasma Thruster (PPT) is a micro-propulsion device (1–100 W average power, μN to mN thrust) that ablates solid polytetrafluoroethylene (Teflon) across parallel rail electrodes in microsecond bursts for satellite attitude control, whereas an MPD thruster is a coaxial, continuous or quasi-steady gas/liquid-metal engine designed for primary orbit transfer and interplanetary propulsion from 10 kW to 10 MW.

1. Jahn, R. G. (1968). Physics of Electric Propulsion. McGraw-Hill Book Company (Dover Publications reprint, 2006). https://ui.adsabs.harvard.edu/abs/1968pep..book.....J
2. Maecker, H. (1955). Plasmaströmungen in Lichtbögen infolge eigenmagnetischer Kompression. Zeitschrift für Physik, 141(1–2), 198–216. https://doi.org/10.1007/BF01327300
3. Choueiri, E. Y. (1998). Scaling of thrust in self-field magnetoplasmadynamic thrusters. Journal of Propulsion and Power, 14(5), 744–753. https://doi.org/10.2514/2.5337
4. Kodys, A. D., & Choueiri, E. Y. (2005). A critical review of the state-of-the-art in the performance of applied-field magnetoplasmadynamic thrusters. 41st AIAA/ASME/SAE/ASEE Joint Propulsion Conference, AIAA 2005-4247. https://doi.org/10.2514/6.2005-4247
5. Boxberger, A., Behnke, A., & Herdrich, G. (2021). Current advances in optimization of operational envelope of steady-state applied-field MPD thruster SX3. International Electric Propulsion Conference (IEPC), IEPC-2022-304. https://doi.org/10.13009/EUCASS2022-7325