2026-08-19 | | Total: 9
Magnetic reconnection governs the explosive release of magnetic energy in systems from the solar corona to fusion plasmas, yet controlling it in the laboratory has remained out of reach. Here we demonstrate active control of reconnection in high-power laser-driven plasmas using a third, relativistic-intensity laser pulse that injects filaments of electron current into the reconnecting system. Two moderate-intensity lasers drive colliding magnetized plumes that reconnect, forming plasmoids in the current sheet as seen in proton deflectometry. The relativistic laser generates magnetic fields matching the polarity on either side of the layer, and, depending on its arrival time, either accelerates the breakup of the current sheet or suppresses reconnection. Arriving early, before the plumes strongly interact, it builds a pocket of magnetic pressure that repels them via flux pileup; arriving after the current sheet forms, it accelerates electrons that extend current filamentation instabilities into the upstream, causing rapid dissipation of the reconnecting magnetic field. This approach opens a route to steering magnetic energy flow in fusion plasmas and broadens the range of systems accessible to laboratory astrophysics.
Direct laser acceleration (DLA) in near-critical-density (NCD) plasmas can efficiently generate high-charge relativistic electron beams, yet beam parameters depend critically on precise plasma state manipulation. Solid-ablation NCD plasmas evolve rapidly, posing severe controllability challenges. We produce NCD plasma via indirectly heating foam targets with ns laser driven hohlraum soft X-ray. Electrons are generated through irradiating the plasma with another picosecond laser. Tuning the laser pulse delay $τ$ enables control of plasma profiles and beam parameters. Experiments show that when the foam is heated ($τ$ = 6 ns, 9 ns), the beam exhibits $T \sim 13$ MeV effective temperature, $E_k \sim 80$ MeV cutoff energy, and hundreds of nC/sr charge for $E_k > 7.5$ MeV. These values are significantly higher than those from solid-foil ($T$ $\sim$ 2.7 MeV, $E_k$ $\sim$ 20 MeV, $Q$ $\sim$ 9 nC/sr) and cold-foam ($T$ $\sim$ 12 MeV, $E_k$ $\sim$ 50 MeV, $Q$ $\sim$ 5 nC/sr) interactions. At a longer delay of $τ$ = 15 ns, the charge increases further while the temperature decreases, and at a shorter delay of $τ$ = 3 ns, both temperature and charge are lower. 3D PIC simulations link these observations to the interplay between the microstructure of the cold foam and the evolving plasma density profile at different delay times, which together determine the beam charge, effective temperature, and divergence. The finding provides a routine to generate and tailor the relativistic electron beams, which is essential for designing laser-driven electron sources for high energy density physics and photonuclear reaction applications.
In this paper, the authors delve into the collisionless kinetic theory of planar vacuum diodes to investigate how a pulse-periodic (ectonic) cathode plasma emission mechanism influences the non-thermal expansion of the cathode plasma and "anomalous" ion acceleration. Another aspect explored in this paper is how pulse-periodic emission affects the comparative dynamics of motion of multi-charged ions from the cathode material. Using kinetic theory, we conclude that the ectonic emission mechanism affects both the macroscopic picture of plasma expansion and the intricacies of electric potential redistribution over one complete emission period.
Pump-probe experiments using dual ultrashort X-ray pulses provide unique opportunities for resolving non-equilibrium dynamics initiated by intense X-ray excitation. Betatron radiation from laser wakefield accelerators offers femtosecond duration, micrometer-scale source size, and intrinsic synchronization with the driving laser, making it a promising candidate for compact ultrafast X-ray sources. Here, we experimentally demonstrate a high-flux, dual-pulse betatron X-ray source based on a density-tailored gas-mixture target. Two electron bunches are generated within a single plasma wakefield through ionization-induced and shock-front-triggered injection, subsequently producing twin X-ray pulses. The measured electron spectra and dual-component X-ray angular profiles, together with particle-in-cell simulations, identify the contributions of the two electron populations to the radiation. The total X-ray photon yield reaches the level of 10^{10} photons per shot with a 40-TW laser system. These results establish a compact, single-stage route toward high-flux dual-pulse betatron sources for laboratory-scale ultrafast X-ray spectroscopy.
Fast analysis of microscopic drift-wave instabilities based on linear gyrokinetic simulation is desirable for modeling anomalous transport in fusion device. In this work, we present an orbit-invariant decomposition method for solving collisionless gyrokinetic eigenvalue problems. By discretizing velocity space along orbit invariants using particle energy and magnetic moment, the full eigenvalue matrix is separated into independent orbit blocks that couple with each other through the field equation, greatly reducing both matrix dimension and computational cost without sacrificing physics. Based on this method, we extend the MGK code [Phys.\ Plasmas 24, 072106 (2017)] with both CPU and GPU implementations, supporting collisionless electrostatic linear simulations in $s$--$α$ and Miller equilibrium model with kinetic. For kinetic ion temperature gradient (ITG) and trapped electron mode (TEM) eigenvalue problems, the solver reduces single-solution times to the 0.01--0.1s range---more than three orders of magnitude faster than CGYRO on the same hardware---enabling efficient large-scale parameter scans. The eigenfrequencies and mode structures are verified by comparing with CGYRO results. The method is generally adapt to to all collisionless gyrokinetic eigenvalue formulations and can be extended to fully electromagnetic simulations.
The Facility for Laboratory Reconnection Experiments (FLARE) has been constructed to study magnetic reconnection in multiple X-line regimes relevant to space, astrophysical, and fusion plasmas. Building upon the successful design of the Magnetic Reconnection Experiment (MRX), FLARE features a larger physical volume, stronger magnetic fields, and an independent ohmic heating drive to significantly extend the accessible parameter space, targeting Lundquist numbers up to S ~ 10^5 and normalized system sizes up to λ~ 10^3. This paper details the facility's core engineering components, including the primary vacuum vessel, internal flux cores, highly segmented external coil systems, modular capacitor banks, and the safety interlock and control architecture. An initial diagnostic suite is presented, comprising high-resolution 2D magnetic probe arrays, triple Langmuir probes, a fully fiber-coupled interferometer, ion Doppler spectroscopy, and fast camera imaging. Initial operations demonstrate the device's experimental flexibility and reliability, successfully executing symmetric push-pull reconnection, spheromak merging, and asymmetric downstream configurations. Currently operating within "Stage 2.5" with S ~ 2,500 and λ~ 60 for anti-parallel reconnection, FLARE provides immediate access to the multiple X-line regimes. Planned hardware upgrades, advanced diagnostic additions, and integration with fully kinetic simulations will further expand its capabilities as it transitions into a collaborative user facility for the broader plasma science community.
Stellarator design explores a vast space of three-dimensional plasma boundaries, only a small fraction of which yields usable equilibria. Data-driven models can narrow this search by learning from existing optimized configurations. Building on the ConStellaration database, we extend conditional boundary generation to four-field-period QI configurations, focusing on the sparsely sampled low-aspect-ratio regime. The approach learns the common geometric structure of known QI equilibria and then adapts it using a small high-fidelity compact dataset, allowing target magnetic properties to guide generation beyond the original data distribution. This adaptation reduces the compact-domain test loss by approximately 87% and yields converged ultra-compact candidates consistent with the prescribed conditions. Several candidates show favorable confinement indicators, and one provides a useful seed for further QI optimization and finite-beta assessment. The method therefore serves as a data-informed front end to high-fidelity physics and optimization.
Successive innovations in particle accelerators have continually expanded the frontiers of scientific discovery. Laser wakefield accelerators promise to transform science, medicine, and industry, yet moving them from laboratory demonstrations to reliable real-world operation has remained a central, long-standing challenge. Here we report a field-deployable system that produced 100-MeV-class electron beams with 1%-level energy stability during 72 hours of continuous operation and supported routine full-power use throughout a seven-month field trial in an industrial setting. Applied to in situ micro-nondestructive testing, the system generated tens-of-MeV bremsstrahlung X-rays that enabled three-dimensional microtomography of dense materials at sub-50-μm spatial resolution and revealed 100-μm-scale internal defects in large composite structures, extending the capabilities beyond those of existing high-energy X-ray sources. These results mark a transition of laser wakefield acceleration from laboratory proof of concept toward practical deployment in scientific and industrial applications.
This work applies Bayesian optimization to a loop composed of PIC simulations of laser electron acceleration and Monte Carlo (MC) simulations of bremsstrahlung-induced nuclear reactions, to maximize the production of molybdenum-99, the precursor of the most used radiopharmaceutical in nuclear medicine, metastable technetium-99. PIC and MC simulations are computationally intensive, and besides reducing the time spent, the Bayesian optimization coupling both simulations resulted in an improvement of an order of magnitude in the $^\text{99}$Mo yield over a previous work, in which the output of an optimization loop based solely on PIC simulations was used a posteriori to estimate $^{99}\mathrm{Mo}$ production through a MC simulation.