Magnetic reconnection is an important process in astrophysical environments, as it re-configures magnetic field topology and converts magnetic energy into thermal and kinetic energy. In extreme astrophysical systems, such as black hole coronae and pulsar magnetospheres, radiative cooling modifies the energy partition by radiating away internal energy, which can lead to the radiative collapse of the layer. In this paper, we perform 2D & 3D simulations to model the MARZ (Magnetic Reconnection on Z) experiments, which are designed to access cooling rates in the laboratory necessary to investigate reconnection in a previously unexplored radiatively-cooled regime. These simulations are performed in GORGON, an Eulerian resistive magnetohydrodynamic code, which models the experimental geometry comprising two exploding wire arrays driven by 20 MA of current on the Z machine (Sandia National Laboratories). Radiative losses are implemented using non-local thermodynamic equilibrium tables computed using the atomic code Spk, and we probe the effects of radiation transport by implementing both a local radiation loss model and P$_{1/3}$ multi-group radiation transport. The load produces highly collisional, super-Alfv\'enic $(M_{A} \approx 1.5)$, supersonic $(M_S \approx 4-5)$ plasma flows which generate a reconnection layer ($L/{\delta} \approx 100, S_L \approx 400$). The reconnection layer undergoes radiative collapse when the radiative losses exceed Ohmic and compressional heating ${\tau}_{cool}^{-1}/{\tau}_A^{-1} \approx 100$; this generates a cold strongly compressed current sheet, leading to an accelerated reconnection rate, consistent with theoretical predictions. Finally, the current sheet is unstable to the plasmoid instability, but the magnetic islands are extinguished by strong radiative cooling before ejection from the layer.
We present evidence for strong radiative cooling in a pulsed-power-driven magnetic reconnection experiment. Two aluminum exploding wire arrays, driven by a 20 MA peak current, 300 ns rise time pulse from the Z machine (Sandia National Laboratories), generate strongly driven plasma flows (MA≈7) with anti-parallel magnetic fields, which form a reconnection layer (SL≈120) at the mid-plane. The net cooling rate far exceeds the Alfvénic transit rate (τcool−1/τA−1≫1), leading to strong cooling of the reconnection layer. We determine the advected magnetic field and flow velocity using inductive probes positioned in the inflow to the layer, and inflow ion density and temperature from analysis of visible emission spectroscopy. A sharp decrease in x-ray emission from the reconnection layer, measured using filtered diodes and time-gated x-ray imaging, provides evidence for strong cooling of the reconnection layer after its initial formation. X-ray images also show localized hotspots, regions of strong x-ray emission, with velocities comparable to the expected outflow velocity from the reconnection layer. These hotspots are consistent with plasmoids observed in 3D radiative resistive magnetohydrodynamic simulations of the experiment. X-ray spectroscopy further indicates that the hotspots have a temperature (170 eV) much higher than the bulk layer (≤75 eV) and inflow temperatures (about 2 eV) and that these hotspots generate the majority of the high-energy (>1 keV) emission.
An X-pinch load driven by an intense current pulse (>100 kA in ∼100 ns) can result in the formation of a small radius, runaway compressional micro-pinch. A micro-pinch is characterized by a hot (>1 keV), current-driven (>100 kA), high-density plasma column (near solid density) with a small neck diameter (1–10 µm), a short axial extent (<1 mm), and a short duration (≲1 ns). With material pressures often well into the multi-Mbar regime, a micro-pinch plasma often radiates an intense, sub-ns burst of sub-keV to multi-keV x rays. A low-density coronal plasma immediately surrounding the dense plasma neck could potentially shunt current away from the neck and thus reduce the magnetic drive pressure applied to the neck. To study the current distribution in the coronal plasma, a Faraday rotation imaging diagnostic (1064 nm) capable of producing simultaneous high-magnification polarimetric and interferometric images has been developed for the MAIZE facility at the University of Michigan. Designed with a variable magnification (1–10×), this diagnostic achieves a spatial resolution of ∼35 µm, which is useful for resolving the ∼100-μm-scale coronal plasma immediately surrounding the dense core. This system has now been used on a reduced-output MAIZE (100–200 kA, 150 ns) to assess the radial distribution of drive current immediately surrounding the dense micro-pinch neck. The total current enclosed was found to increase as a function of radius, r, from a value of ≈50±25 kA at r ≈ 140 µm (at the edge of the dense neck) to a maximal value of ≈150±75 kA for r ≥ 225 µm. This corresponds to a peak magnetic drive pressure of ≈75±50 kbar at r ≈ 225 µm. The limitations of these measurements are discussed in the paper.
We present the first experimental study of plasmoid formation in a magnetic reconnection layer undergoing rapid radiative cooling, a regime relevant to extreme astrophysical plasmas. Two exploding aluminum wire arrays, driven by the Z machine, generate a reconnection layer (S_{L}≈120) in which the cooling rate far exceeds the hydrodynamic transit rate (τ_{hydro}/τ_{cool}>100). The reconnection layer generates a transient burst of >1 keV x-ray emission, consistent with the formation and subsequent rapid cooling of the layer. Time-gated x-ray images show fast-moving (up to 50 km s^{-1}) hotspots in the layer, consistent with the presence of plasmoids in 3D resistive magnetohydrodynamic simulations. X-ray spectroscopy shows that these hotspots generate the majority of Al K-shell emission (around 1.6 keV) prior to the onset of cooling, and exhibit temperatures (170 eV) much greater than that of the plasma inflows and the rest of the reconnection layer, thus providing insight into the generation of high-energy radiation in radiatively cooled reconnection events.
We present a detailed characterization of the structure and evolution of differentially rotating plasmas driven on the MAGPIE pulsed-power generator (1.4 MA peak current, 240 ns rise time). The experiments were designed to simulate physics relevant to accretion discs and jets on laboratory scales. A cylindrical aluminium wire array Z pinch enclosed by return posts with an overall azimuthal off-set angle was driven to produce ablation plasma flows that propagate inwards in a slightly off-radial trajectory, injecting mass, angular momentum and confining ram pressure to a rotating plasma column on the axis. However, the plasma is free to expand axially, forming a collimated, differentially rotating axial jet that propagates at ${\approx }100\,{\rm km}\,{\rm s}<^>{-1}$. The density profile of the jet corresponds to a dense shell surrounding a low-density core, which is consistent with the centrifugal barrier effect being sustained along the jet's propagation. We show analytically that, as the rotating plasma accretes mass, conservation of mass and momentum implies plasma radial growth scaling as $r \propto t<^>{1/3}$. As the characteristic moment of inertia increases, the rotation velocity is predicted to decrease and settle on a characteristic value ${\approx }20\,{\rm km}\,{\rm s}<^>{-1}$. We find that both predictions are in agreement with Thomson scattering and optical self-emission imaging measurements.
Recent pulsed-power experiments have demonstrated the formation of astrophysically relevant, differentially rotating plasmas. Key features of the plasma flows are the discovery of a quasi-Keplerian rotation curve, the launching of highly collimated angular-momentum-transporting axial jets, and a hollow density structure sustained by the centrifugal barrier effect. In this communication, we discuss several features of the plasma structure in these experiments through order-of-magnitude models. First, we show that the observed rotation velocity would produce a centrifugal force strong enough to support the hollow density profile. Second, we show that the axial jet should diverge much faster than what was observed, were it not for a magnetized halo with 3 T which surrounds the jet and exerts pressure on the interface. Finally, we discuss the temperature structure in the axial jet and plasma halo. We show that a 3-T magnetic field would also suppress electron heat conduction, leading to the flat profile observed experimentally. We also find that the axial jet is efficiently radiatively cooled, whereas the halo is not, which would explain the thermal decoupling between the two regions.
In order to develop a better understanding of current losses in the magnetically insulated region of high-power pulsed power machines, it is crucial to be able to conduct experiments at scale in smaller facilities. Here, we present a novel experiment that has been tested on the MAGPIE driver at Imperial College. The targets are inexpensive and easy to customize for experiments ranging from power flow to warm dense matter. Simulations in COMSOL indicated electric fields of up to 600 MV/m and magnetic fields of up to 300 T could be produced on MAGPIE. In initial testing, Electric fields exceeding 100 MV/m and magnetic fields exceeding 50 T were generated, and both magnetically insulated transmission and plasma-shorted transmission were demonstrated.
We present results from pulsed-power driven differentially rotating plasma experiments designed to simulate physics relevant to astrophysical disks and jets. In these experiments, angular momentum is injected by the ram pressure of the ablation flows from a wire array Z pinch. In contrast to previous liquid metal and plasma experiments, rotation is not driven by boundary forces. Axial pressure gradients launch a rotating plasma jet upward, which is confined by a combination of ram, thermal, and magnetic pressure of a surrounding plasma halo. The jet has subsonic rotation, with a maximum rotation velocity 23±3 km/s. The rotational velocity profile is quasi-Keplerian with a positive Rayleigh discriminant κ^{2}∝r^{-2.8±0.8} rad^{2}/s^{2}. The plasma completes 0.5-2 full rotations in the experimental time frame (∼150 ns).
Shock structures formed in supersonic plasmas colliding against conducting obstacles can be significantly affected by radiative cooling and magnetic fields, leading to instabilities and turbulence. Under these conditions, experimental measurements of plasma parameters are very difficult to perform using conventional laser-based diagnostics, such as shadowgraphy, schlieren and interferometry.
We present a study of subcritical shocks in a highly collisional laboratory plasma with a dynamically significant magnetic field. Shocks were produced by placing cylindrical obstacles into the supermagnetosonic ( $M_{{\rm MS}} \sim 1.9$ ) outflow from an inverse wire array z-pinch at the MAGPIE pulsed power facility ( $n_{e} \sim 8.5 \times 10^{17}\,{\rm cm}^{-3}$ , $v \sim 45\,{\rm km}\,{\rm s}^{-1}$ ). We demonstrate the existence of subcritical shocks in this regime and find that secondary stagnation shocks form in the downstream which we infer from interferometry and optical Thomson scattering measurements are hydrodynamic in nature. The subcritical shock width is found to be approximately equal to the resistive diffusion length and we demonstrate the absence of a jump in hydrodynamic parameters. Temperature measurements by collective optical Thomson scattering showed little temperature change across the subcritical shock ( ${<}10\,\%$ of the ion kinetic energy) which is consistent with a balance between adiabatic and Ohmic heating and radiative cooling. We demonstrate the absence of subcritical shocks when the obstacle diameter is less than the resistive diffusion length due to decoupling of the magnetic field from the plasma. These findings are supported by magnetohydrodynamic simulations using the Gorgon and AstroBEAR codes and discrepancies between the simulations and experiment are discussed.
We study the structure of reverse shocks formed by the collision of supersonic, magnetized plasma flows driven by an inverse (or exploding) wire array with a planar conducting obstacle. We observe that the structure of these reverse shocks varies dramatically with wire material, despite the similar upstream flow velocities and mass densities. For aluminum wire arrays, the shock is sharp and well-defined, consistent with magneto-hydrodynamic theory. In contrast, we do not observe a well-defined shock using tungsten wires, and instead we see a broad region dominated by density fluctuations on a wide range of spatial scales. We diagnose these two very different interactions using interferometry, Thomson scattering, shadowgraphy, and a newly developed imaging refractometer that is sensitive to small deflections of the probing laser corresponding to small-scale density perturbations. We conclude that the differences in shock structure are most likely due to radiative cooling instabilities, which create small-scale density perturbations elongated along magnetic field lines in the tungsten plasma. These instabilities grow more slowly and are smoothed by thermal conduction in the aluminum plasma.
The physics of astrophysical jets can be divided into three regimes: (i) engine and launch (ii) propagation and collimation, (iii) dissipation and particle acceleration. Since astrophysical jets comprise a huge range of scales and phenomena, practicality dictates that most studies of jets intentionally or inadvertently focus on one of these regimes, and even therein, one body of work may be simply boundary condition for another. We first discuss long standing persistent mysteries that pertain the physics of each of these regimes, independent of the method used to study them. This discussion makes contact with frontiers of plasma astrophysics more generally. While observations theory, and simulations, and have long been the main tools of the trade, what about laboratory experiments? Jet related experiments have offered controlled studies of specific principles, physical processes, and benchmarks for numerical and theoretical calculations. We discuss what has been accomplished on these fronts. Although experiments have indeed helped us to understand certain processes, proof of principle concepts, and benchmarked codes, they have yet to solved an astrophysical jet mystery on their own. A challenge is that experimental tools used for jet-related experiments so far, are typically not machines originally designed for that purpose, or designed with specific astrophysical mysteries in mind. This presents an opportunity for a different way of thinking about the development of future platforms: start with the astrophysical mystery and build an experiment to address it.
We present a technique to measure the time-resolved velocity and ion sound speed in magnetized, supersonic high-energy-density plasmas. We place an inductive ("b-dot") probe in a supersonic pulsed-power-driven plasma flow and measure the magnetic field advected by the plasma. As the magnetic Reynolds number is large (RM > 10), the plasma flow advects a magnetic field proportional to the current at the load. This enables us to estimate the flow velocity as a function of time from the delay between the current at the load and the signal at the probe. The supersonic flow also generates a hydrodynamic bow shock around the probe, the structure of which depends on the upstream sonic Mach number. By imaging the shock around the probe with a Mach-Zehnder interferometer, we determine the upstream Mach number from the shock Mach angle, which we then use to determine the ion sound speed from the known upstream velocity. We use the sound speed to infer the value of Z̄Te, where Z̄ is the average ionization and Te is the electron temperature. We use this diagnostic to measure the time-resolved velocity and sound speed of a supersonic (MS ∼ 8), super-Alfvénic (MA ∼ 2) aluminum plasma generated during the ablation stage of an exploding wire array on the Magpie generator (1.4 MA, 250 ns). The velocity and Z̄Te measurements agree well with the optical Thompson scattering measurements reported in the literature and with 3D resistive magnetohydrodynamic simulations in GORGON.
We present a study of perpendicular subcritical shocks in a collisional laboratory plasma. Shocks are produced by placing obstacles into the supermagnetosonic outflow from an inverse wire array z pinch. We demonstrate the existence of subcritical shocks in this regime and find that secondary shocks form in the downstream. Detailed measurements of the subcritical shock structure confirm the absence of a hydrodynamic jump. We calculate the classical (Spitzer) resistive diffusion length and show that it is approximately equal to the shock width. We measure little heating across the shock (<10% of the ion kinetic energy) which is consistent with an absence of viscous dissipation.
Heavy ion beam probe (HIBP) diagnostics on TUMAN-3M tokamak was updated to provide measurements in the regime with neutral beam injection co-directed with plasma current (co-NBI). By means of HIBP plasma potential measurements in the center of plasma were carried out. Plasma potential evolution in the discharge with L-H transition (transition to improved confinement mode) is in good agreement with the concept of negative radial electric field generation during the formation of transport barrier.
V. Valenzuela-Villaseca, ∗ L. G. Suttle, F. Suzuki-Vidal, J. W. D. Halliday, S. Merlini, D. R. Russell, E. R. Tubman, † J. D. Hare, ‡ J. P. Chittenden, M. E. Koepke, E. G. Blackman, and S. V. Lebedev Blackett Laboratory, Imperial College, London SW7 2BW, UK. Department of Physics, West Virginia University, USA. Department of Physics and Astronomy, University of Rochester, USA. (Dated: January 26, 2022)