Magnetized liner inertial fusion experiments on the Z accelerator suffer from magneto-Rayleigh–Taylor instabilities (MRTI) that compromise integrity of the imploding cylindrical liner, limiting achievable fusion fuel conditions and ultimately reducing magneto-inertial fusion target performance. Dynamic screw pinches (DSP) provide a method to reduce MRTI in-flight via application of magnetic field line tension to the imploding liner outer surface. In contrast with z-pinches that drive implosions with an azimuthal magnetic field, dynamic screw pinches enforce an additional axial drive magnetic field component, making the overall drive magnetic field helical. As the liner implodes, cumulative MRTI development is reduced by dynamically shifting the orientation of the fastest growing instability modes. Three-dimensional magnetohydrodynamic simulations show that the DSP mechanism effectively stabilizes initially solid cylindrical liner implosions driven by Z-scale current pulses, indicating that MRTI mitigation increases with the ratio of axial to azimuthal drive magnetic field components (i.e., the drive field ratio). We also performed a spectral analysis of the simulated imploding density distributions, extracting wavelength and pitch angle of the simulated MRTI structures to study their dynamics during the implosion. Simulations of liners initially perturbed with drive-field-aligned sinusoidal structures indicate that MRTI mitigation in DSP implosions decreases with perturbation wavelength, once again suggestive of magnetic field line tension effects.

## I. INTRODUCTION

Improving implosion stability in magnetized liner inertial fusion (MagLIF)^{1,2} is predicted to improve fusion fuel compression and confinement. Magneto-Rayleigh–Taylor instabilities (MRTI) limit the compressive work that the imploding liner can do on the fuel, reducing achievable fuel temperatures and pressures and therefore thermonuclear yields. Additionally, instability of the imploding liner results in degradation of fuel confinement,^{3} reducing the time during which thermonuclear reactions can take place prior to disassembly.

Dynamic screw pinches (DSP)^{4–8} are designed to stabilize magnetically driven cylindrical implosions by employing a dynamic helical magnetic drive field. As the target implodes, the polarization of the magnetic field driving the implosion rotates, resulting in the dynamic reorientation of magnetic field line tension over the outer surface of the imploding liner, continuously shifting the fastest growing instability modes and reducing cumulative MRTI growth. Linear theory^{4} indicates how the DSP mechanism is predicted to mitigate MRTI and experiments studying DSP-driven thin foil implosions^{6} demonstrate stabilization compared to thin foil z-pinches. The numerical study of liner implosions more directly relevant to MagLIF experiments on the Z accelerator^{9,10} provides assessment of the DSP mechanism in the presence of physical effects not captured in the linear theory such as non-linear instability growth and mode merging,^{11,12} magnetic field diffusion into finite electrical conductivity materials, and the effects of shocks transmitting through liner material and material interfaces.

In this manuscript, we utilize three-dimensional magnetohydrodynamic (3D MHD) simulations to explore how effectively the DSP mechanism stabilizes magnetically driven cylindrical implosions. In Sec. II, we compare the stability of simulated implosions on the Z accelerator with respect to the initial ratio of the axial and azimuthal drive field components evaluated at the initial outer liner surface; a larger ratio of axial to azimuthal drive magnetic field components is predicted to result in more effective stabilization. We then compare the spectral characteristics of the imploding liner perturbations, focusing particularly on the evolution of the dominant wavelength and average pitch angle of instability structures and how these correspond to the effectiveness of the DSP mechanism. Section III briefly describes results from simulations exploring a special case: liner implosions with initial sinusoidal perturbations on the outer liner surface. By prescribing a dominant wavelength via a single mode sinusoidal perturbation, the dependence of the DSP mechanism on wavenumber can be explored without relying on instability structures growing from a white noise random MRTI seed with maximum wavenumber limited by the simulation mesh resolution.

## II. SIMULATED DYNAMIC SCREW PINCH STABILIZATION IN INITIALLY-SOLID METALLIC LINER IMPLOSIONS

### A. Cumulative MRTI growth vs magnetic drive field ratio

To assess the stability of magnetically driven cylindrical implosions, we primarily used three-dimensional (3D) ALEGRA.^{13,14} ALEGRA is an arbitrary Lagrangian–Eulerian (ALE) hydrodynamics code that utilizes resistive MHD with thermal conduction. Simulations employed SESAME equation of state (EOS) and coupled thermal/electrical conductivity tables. The “baseline” simulation geometry used for this study (Fig. 1) was a 1-mm tall simulation region with periodic top and bottom boundary conditions and a base mesh resolution of 20 $\mu $m. Simulations employed a “radial trisection” mesh, with radial mesh cells (i.e., mesh cell edges oriented along the radial, azimuthal, and axial directions) converted to Cartesian geometry near the axis. Note that radial mesh cells correspond to increasing mesh cell resolution in the azimuthal direction as radius decreases (until the conversion to Cartesian geometry). The liner was modeled as beryllium with an initial outer radius $ R out 0$ = 2.616 mm and an initial inner radius $ R in 0$ = 2.325 mm, corresponding to an aspect ratio (AR) of 9, where AR $= R in 0 / R out 0 \u2212 R in 0$. Simulations were driven with a prescribed 20-MA peak current, 100 ns rise time sine squared pulse: $I t= 20 \u2009 MA* sin 2 \pi t / \tau $, where $\tau =2 t r$ = 200 ns and $ t r$ is the rise time. A “white noise” surface roughness was applied to the initial outer surface of the liner to seed MRTI growth by randomly perturbing mesh nodes conformal to $ R out 0$ radially over a ±1.25 $\mu $m range (uniform distribution).

While this MRTI seeding mechanism is consistent across the main ALEGRA simulations presented in this study (we explored several variations to a limited degree), we note that it is not rigorously benchmarked against prior implosion data from unmagnetized^{11} and magnetized^{15,16} initially solid liner implosion experiments on Z. This is partially because the mechanisms that give rise to the observed MRTI in prior experiments are not fully understood in terms of providing a self-consistent seeding protocol to reproduce the high levels of azimuthal correlation observed in unmagnetized experiments^{11} or the helical structures^{15–17} observed in magnetized implosions. We note that several hypotheses have been provided to explain the discrepancies such as electrothermal instabilities (ETI) exacerbating the effective surface roughness^{18,19} and local plasma induced flux compression^{20–22} in the target region giving rise to helical structures. Future work will focus on studying the DSP mechanism in the context of these MRTI seeding hypotheses, but we restrict our study here to a seeding protocol implemented in a resistive MHD numerical framework kept consistent across the main simulations presented in this study.

We focus our study on “empty” beryllium liners (in contrast with premagnetized, preheated fusion-fuel-filled liners associated with “integrated” MagLIF experiments^{1,2}) and restrict our assessment of MRTI up to levels of convergence lower than those associated with peak compression in MagLIF. Specifically, we report on MRTI development up until the liner has imploded to CR_{in} < 10 in this study vs CR_{in} > 30 associated with peak compression in MagLIF experiments, where CR_{in} is the convergence ratio of the inner liner surface, $ R in 0/ R in t$. This is primarily due to the emphasis of this study on understanding acceleration driven MRTI development at the outer liner surface instead of instability feedthrough to the inner surface and its effects on fusion fuel conditions. Such development of inner surface perturbations is expected to be much less in our simulated system compared to MagLIF experiments due to the lower convergences studied and the lack of hydrodynamic pressure from the preheated/premagnetized fusion fuel that would induce deceleration of the inner liner surface at these higher convergences. The assessment of the development of inner surface perturbations resulting from feedthrough of MRTI from the outer liner surface to the inner liner surface during a fusion-fuel-filled liner implosion and the effects on fuel compression and confinement will be the subject of future study.

. | $ \Xi o=0$ . | $ \Xi o=0.25$ . | $ \Xi o=0.5$ . | $ \Xi o=1.0$ . |
---|---|---|---|---|

$ c 1$ | 1.0 | 1.0 | 1.025 | 1.075 |

$ c 2$ | 0 | 0.243 | 0.447 | 0.707 |

$ c 3$ | 1.0 | 0.970 | 0.894 | 0.707 |

. | $ \Xi o=0$ . | $ \Xi o=0.25$ . | $ \Xi o=0.5$ . | $ \Xi o=1.0$ . |
---|---|---|---|---|

$ c 1$ | 1.0 | 1.0 | 1.025 | 1.075 |

$ c 2$ | 0 | 0.243 | 0.447 | 0.707 |

$ c 3$ | 1.0 | 0.970 | 0.894 | 0.707 |

^{3}density contour. A 2D map of the areal density $\rho r$ of the imploding liner was produced by radially integrating the mass density at each ( $\theta $, $z$) coordinate in the simulation, where $\theta $ is the azimuthal coordinate and $z$ is the axial coordinate. $\sigma (\rho r)$ was calculated for each 2D $\rho r$ map at select simulation times.

. | $ G \sigma , 0.25$ . | $ G \sigma , 0.5$ . | $ G \sigma , 1.0$ . |
---|---|---|---|

ALEGRA | 3.0 | 6.0 | 16.0 |

KRAKEN | 2.3 | 4.2 | 7.4 |

. | $ G \sigma , 0.25$ . | $ G \sigma , 0.5$ . | $ G \sigma , 1.0$ . |
---|---|---|---|

ALEGRA | 3.0 | 6.0 | 16.0 |

KRAKEN | 2.3 | 4.2 | 7.4 |

Simulated DSP implosions with $ \Xi o$ = 0.25, 0.5, and 1.0 demonstrate 3X, 6X, and 16X lower cumulative MRTI, respectively, once the liner has imploded to $\Delta R$ = 0.8 according to the $\sigma (\rho r)$ analysis (Fig. 5, Table II). While the 16X reduction according to $\sigma (\rho r$) analysis from $ \Xi o$ = 1 simulations does not match the 170X theoretical prediction^{4} for a $\lambda $ = 200 $\mu $m instability mode, we emphasize that the theory is focused on linear mode amplitude for a single wavelength mode without accounting for effects like nonlinear instability saturation, mode coupling, material compression, and magnetic diffusion. The theoretical predictions^{4} for cumulative MRTI reduction of $\lambda =400\u2009\mu $m and $\lambda =1000\u2009\mu $m modes (based on the assessment of the reduction in MRTI amplitude growth), 18X and 3X, respectively, can be more closely compared to $\sigma (\rho r$) analysis from simulation results in this study, potentially highlighting how the dynamic multi-mode instability spectrum in simulations influences the DSP mechanism. However, spectral analysis described in the next subsection indicates that the dominant instability wavelength reached in DSP implosion simulations is in the 80–200 $\mu $m range, so comparison with theoretical predictions of single mode MRTI reduction for $\lambda >200\u2009\mu $m remains dubious. The lack of quantitative agreement of cumulative MRTI reduction from theoretical predictions with simulation predictions does not diminish the fact that simulations agree with the qualitative prediction from the theory: the effectiveness of the DSP mechanism to reduce MRTI growth increases as the initial drive field ration $ \Xi o$ is greater.

To assess the robustness of the trends observed in the “baseline” simulations, we completed >80 additional 3D ALEGRA simulations focused on varying initial conditions and numerical modeling choices. Initial conditions and physics models we varied included the mesh resolution (12.5, 16, 20 $\mu $m), the surface roughness perturbation random seed, the uniform distribution surface roughness perturbation magnitude (±0.625, ±1.25, ±2.5 $\mu $m), the simulation region height (1 mm vs 1.5 mm), the vacuum treatment (i.e., the density, electrical conductivity, and Joule heating “floors” that serve to turn off physical mechanisms in mesh cells when density falls below a specified value), the type of mesh (i.e., Cartesian mesh vs radial trisection mesh), the type of MRTI seed (i.e., random temperature perturbation vs random density perturbation vs surface roughness perturbation), magnetic field effects on thermal conduction (i.e., anisotropic vs isotropic thermal conductivity), radiation transport (i.e., none vs gray diffusion), EOS tables, strength/yield models, and two-temperature physics. Throughout all the variations studied, the cumulative MRTI growth still decreased with $ \Xi o$ as we observed in the baseline simulations and qualitatively reproduce the trends observed in baseline simulations. Future exploration will focus on determining which simulation settings, MRTI seeding mechanism, mesh type, etc., best reproduce experimental data from Z implosions.

We also executed simulations including a background axial magnetic field similar in magnitude to that used in historical MagLIF experiments ( $ B z$ = 15 T) to compare to the unmagnetized baseline simulations ( $ B z$ = 0) of this study. Such background magnetic fields have been shown to influence implosion dynamics^{15,16,20–22} in z-pinch implosions. However, in the resistive MHD numerical system, we study here and without introduction of conductive plasmas or higher order transport terms (e.g., Hall), we observe no significant differences in implosion dynamics and MRTI development in magnetized simulations compared to unmagnetized simulations. Assessment of the DSP mechanism using extended MHD and/or particle-in-cell numerical tools should be the subject of future study, especially given the experimental evidence for the effects of background axial magnetic fields in solid liner z-pinch implosions.

To further explore how the DSP mechanism stabilizes implosions, we simulated the same target geometry driven with the same current pulse using KRAKEN, a 3D Eulerian resistive MHD code that implements the same system of electromagnetic equations from GORGON^{23} with finite volume hydrodynamics. KRAKEN's numerical algorithms are fundamentally different from ALEGRA's, providing an opportunity to test whether the trends extracted from ALEGRA DSP simulations are independent of the choice of simulation tool. KRAKEN simulations employed 20 $\mu $m Cartesian mesh cells and seeded MRTI via an initial ±1% random temperature perturbation applied to all beryllium cells in the simulation. In addition, KRAKEN's Cartesian mesh results in “stair stepping” of the mesh cell boundaries near the cylindrical outer liner initial surface, which imparts an additional, artificial, or numerical surface roughness component absent in ALEGRA simulations. KRAKEN simulations did not employ radiation transport or thermal conduction. The KRAKEN simulation region was 12 mm tall (*without* periodicity assigned to the top and bottom boundaries). As in ALEGRA simulations, $ \Xi o= 0 , \u2009 0.25 , \u2009 0.5 , \u2009 1.0$ cases were studied, again with scale factors applied to the drive magnetic field boundary conditions [Eqs. (1) and (2)] to approximately match implosion time and acceleration time history. $\sigma (\rho r)$ was extracted in the same manner as in ALEGRA simulations to assess MRTI cumulative growth throughout the simulated implosions.

KRAKEN simulations indicate levels of MRTI stabilization similar to ALEGRA for each DSP case compared to the z-pinch case (Figs. 5 and 6, Table II). The z-pinch KRAKEN simulation demonstrates lower cumulative MRTI development compared to ALEGRA, which is thought to be a result of the choice of MRTI seeding protocol. Although an exhaustive study of MRTI seeding protocols particular to ALEGRA and to KRAKEN has not been completed as a part of this study, future study could focus on improving agreement of the z-pinch cumulative MRTI results through fine-tuning of the seed (e.g., adjusting the amplitude and spectral content of the node perturbation in ALEGRA or adjusting the magnitude and spectral content of the temperature perturbation in KRAKEN). Despite the differences in numerical algorithms, MRTI seeding methods, and mesh geometry, the results from 3D ALEGRA and 3D KRAKEN display overall agreement. The largest discrepancy between ALEGRA and KRAKEN is a 27% difference for $\sigma (\rho r)$ between the $ \Xi o$ = 0 z-pinch cases at $\Delta R=0.8$. This difference reduces to 18% for the $\Xi $ = 0.25 case, 9% for the $ \Xi o$ = 0.5 case, and 20% for the $ \Xi o$ = 1.0 case. More importantly, the prediction of cumulative MRTI reduction $ G \sigma $ from 3D ALEGRA is comparable to 3D KRAKEN (Table II). ALEGRA predicts more DSP stabilization compared to KRAKEN for all of the DSP implosions compared to the z-pinch case (Table II). However, note that the largest source for the discrepancy in $ G \sigma $ between the codes comes from the observed difference in $\sigma (\rho r)$ for the $ \Xi o$ = 0 z-pinch case, which functions as the baseline for calculating $ G \sigma $.

_{in}$\u2272$ 10) compared to peak compression in MagLIF experiments (CR

_{in}> 30) and our simulated liners lack preheated/premagnetized fusion fuel to provide internal hydrodynamic pressure and deceleration-induced growth of inner surface perturbations. We assess the inner liner surface perturbation as a function of convergence by assigning $ R in t$ to an 1800 kg/m

^{3}density contour nearest to the central axis. We then compute an average perturbation amplitude of the inner surface,

### B. Spectral analysis of MRTI structures

A 2D fast Fourier transform (FFT) was performed on the 2D maps of $\rho r$ extracted at different points of convergence throughout each 3D ALEGRA simulation (Fig. 8). The 2D FFT provides the axial and azimuthal wavenumber information from which one can compute the dominant pitch angle of the instability structures and the wavelength of the dominant instability mode. The azimuthal angular coordinates for each $\rho r$ map are converted to distance along the circumference $C t$ by using the liner inner radius $ R in t$ at each time step via $C t=2\pi R in t$.

^{24}Specifically, the wavelength component in the azimuthal direction (i.e., along the circumference) $ \lambda c=2\pi / k c$ varies with the liner radius: $ \lambda c t\u2009\u221d R l t$. Axially oriented “flute” modes have wavelengths that vary as $ R l t$, while azimuthally oriented “sausage” modes have no direct dependence on $ R l t$. Therefore, the peak in wavelength occurs earlier (at lower convergence) for simulations with instabilities that have larger $ \lambda c$ produced by larger $ \Xi o$. As to why the instability structures initially increase in wavelength, this is due primarily to mode merging

^{11,12}of the initial continuum spectrum of unstable Rayleigh–Taylor modes.

The z-pinch simulations indeed show increasing dominant instability wavelength as the liner converges (Fig. 9, top) as expected.^{11} In contrast, while mode merging occurs in the DSP simulations and results in the *initial* increase in wavelength during implosion, the $ R l t$ dependence of the azimuthal wavelength component later dominates and results in decreasing wavelength at a higher liner convergence. Furthermore, feedthrough of MRTI is approximately proportional to $ e \u2212 k \Delta r t$, where $\Delta r t$ is the in-flight liner thickness and $k$ is the magnitude of the perturbation wavevector, so larger primary $k$ for a given in-flight liner thickness should generate much lower variations in $\rho r$. Importantly, the expected increase in dominant wavelength over time is not captured in the linear theory,^{4} so simulations of DSP implosions provide crucial information about the predicted effectiveness of the DSP mechanism in experiments.

^{25}for MRTI: $ \gamma 2=kg\u2212 2 k \u21c0 \xb7 B \u21c0 2 \mu o \rho $, where $\gamma $ is the exponential growth rate of the MRTI, $g$ is the acceleration of the liner, $\rho $ is the liner material density, and $ \psi k B$ is the angle between $ k \u21c0$ and $ B \u21c0$. The second term in this relation is the principal motivation behind the DSP mechanism and corresponds to the stabilizing effects of magnetic field line tension; dynamically maximizing this term, particularly as $ k \u21c0$ shifts throughout the implosion, corresponds to maximizing MRTI mitigation via the DSP mechanism.

## III. FIELD-ALIGNED SINGLE MODE PERTURBED IMPLOSION SIMULATIONS

To assess development of the most deleterious MRTI modes in z-pinches compared to dynamic screw pinches, namely, those aligned with the drive magnetic field lines, we used 3D ALEGRA to simulate liner implosions with $ \Xi o= 0 , \u2009 0.25 , \u2009 0.5$ in which sinusoidal perturbations aligned with the initial magnetic drive field lines at $ R out t$ were initially applied to the outer liner surface (Fig. 10). A 6 $\mu m$ amplitude was used for the sinusoidal perturbation for all simulations, and four different wavelengths were simulated: $\lambda = 100 , \u2009 150 , \u2009 200 , 250\u2009\mu m$. While field-aligned MRTI modes tend to grow fastest in magnetically driven implosions, we expected lower cumulative MRTI growth in DSP simulations compared to z-pinch simulations. This is because the wavevector of the sinusoidal perturbations applied to the outer surface in the z-pinch case ( $ \Xi o=0$) remains always perpendicular to the azimuthal drive magnetic field throughout the implosion. In contrast, the wavevector of the helically oriented sinusoidal perturbation in the DSP calculations is only initially perpendicular with the direction of the drive magnetic field, and magnetic field line tension effects grow as the implosion progresses. As the DSP implosions proceed, the perturbation pitch angle $\varphi t$ shifts toward the vertical and the drive magnetic field polarization $ \theta B t$ shifts toward the horizontal, increasing stabilizing field-line tension effects and reducing MRTI cumulative growth. The dominant wavenumber (and thus wavelength) remains approximately constant throughout the simulation for all $ \varphi 0$= 0 single-mode simulations, while the dominant wavelength in $ \varphi 0\u2260$ 0 single mode simulations decreases with $ R l t$ as expected due to convergence effects.

The linear theory^{4} suggests that the DSP mechanism is more effective at mitigating growth of shorter wavelength instability structures. This essentially corresponds to stabilizing force imparted by magnetic field line tension scaling inversely with length scale: $ F \u21c0 T= \mu o \u2212 1( B \u21c0 \u2009 \xb7 \u2009 \u2207) B \u21c0\u221d L \u2212 1$, where $ F \u21c0 T$ is the force exerted perpendicular to $ B \u21c0$ magnetic field lines (i.e., due to magnetic field line tension) and *L* is the characteristic length of the system (e.g., the dominant instability wavelength). Accordingly, the DSP mechanism more effectively stabilizes the simulated DSP implosions with shorter wavelengths (Fig. 11 and Table III), except for when comparing the 100 $\mu m$ wavelength perturbation case to the 150 $\mu m$ case. We postulate that this departure from the expected trend of more effective stabilization at shorter wavelengths arises from the increasing ratio of the perturbation amplitude to the wavelength ( $ 6 \u2009 \mu m / 150 \u2009 \mu m=0.04$ compared to $ 6 \u2009 \mu m / 100 \u2009 \mu m=0.06$), which causes the 100 $\mu m$ case to enter the nonlinear growth phase more quickly, diminishing the screw pinch mechanism's cumulative effectiveness. Future study will focus on normalizing for the perturbation amplitude to wavelength ratio to observe whether the trend of more effective DSP stabilization with shorter wavelength is preserved. More generally, this trend should also correspond to a changing effectiveness of the DSP mechanism as the dominant wavelength of instability structures changes during implosion. Simulated implosions initialized with a white noise surface roughness and implosions in reality that have a physical surface roughness are expected to demonstrate an increasing dominant wavelength over time, primarily due to processes like mode merging.^{11,12}

. | $ G \sigma $ at $\Delta R=0.4$ $\lambda =100\u2009\mu $m . | $ G \sigma $ at $\Delta R=0.4$ $\lambda =150\u2009\mu $m . | $ G \sigma $ at $\Delta R=0.4$ $\lambda =200\u2009\mu $m . | $ G \sigma $ at $\Delta R=0.4$ $\lambda =250\u2009\mu $m . |
---|---|---|---|---|

$ \Xi o=0.25$ | 1.29 | 1.30 | 1.23 | 1.09 |

$ \Xi o=0.5$ | 1.50 | 1.65 | 1.47 | 1.38 |

. | $ G \sigma $ at $\Delta R=0.4$ $\lambda =100\u2009\mu $m . | $ G \sigma $ at $\Delta R=0.4$ $\lambda =150\u2009\mu $m . | $ G \sigma $ at $\Delta R=0.4$ $\lambda =200\u2009\mu $m . | $ G \sigma $ at $\Delta R=0.4$ $\lambda =250\u2009\mu $m . |
---|---|---|---|---|

$ \Xi o=0.25$ | 1.29 | 1.30 | 1.23 | 1.09 |

$ \Xi o=0.5$ | 1.50 | 1.65 | 1.47 | 1.38 |

## IV. CONCLUDING REMARKS

3D simulations of dynamic screw pinch (DSP) implosions executed the ALEGRA and KRAKEN MHD codes demonstrate that DSP stabilization increases with the ratio of axial to azimuthal drive magnetic field components (i.e., the drive field ratio). Our results are in qualitative agreement with the linear theory.^{4} The observed stabilization by the DSP mechanism is encouraging, especially since simulations capture physical effects beyond those captured in the linear theory. Spectral analysis of the implosion instability structures indicates that the time-dependent dominant instability wavelength and average instability pitch angle contribute to the operation of the DSP mechanism, specifically by enabling magnetic field line tension to more effectively stabilize shorter wavelength structures compared to longer wavelength structures.

This study suggests that DSP stabilization is a powerful mechanism that can be leveraged to reduce cumulative MRTI development. Experiments to further determine whether and to what extent DSP implosions stabilize MRTI compared to z-pinch implosions will be necessary to improve models and determine best practices for establishing initial conditions such as MRTI seeding mechanisms. Potential improvements for modeling practices could focus on moving away from *ad hoc* MRTI seeding mechanisms such as artificial surface roughness initialization or volumetric temperature or density perturbations. Such improvements could involve finer mesh resolutions, perhaps enabled by local or adaptive mesh refinement to resolve known seeding processes such as the electrothermal instability that operate on the <1 $\mu $m scale.^{18} Additionally, the local plasma environment^{26} is not captured in the simulations shown in this study, but future examination of the effects of plasmas in the target region, such as those sourced by the pulsed power transmission lines, will provide information on whether and how the DSP mechanism is affected.

In experiments, DSP implosions will utilize helical return current geometries^{4,5} to produce the axial drive magnetic field component and will, therefore, introduce additional inductance into the pulsed power circuit. To enable a clean comparison, the implosion time was maintained for all simulations presented in this manuscript. However, a more complete evaluation of the effects of added inductance combined with the added magnetic drive pressure from the additional axial drive component should be pursued. A trade-off exists between added inductance from the axial drive magnetic field component, which will reduce circuit energy coupling from the generator to the imploding z-pinch liner, and the additional magnetic pressure from the axial magnetic drive component. Future circuit-coupled, 3D radiation-MHD simulations will focus on modeling the entirety of the MagLIF system with premagnetized, laser-preheated fusion fuel to evaluate thermonuclear performance in DSP implosions on Z.

## ACKNOWLEDGMENTS

The authors gratefully acknowledge helpful discussions with A. Porwitzky, W. Lewis, L. Stanek, J. Niederhaus, A.C. Robinson, E.P. Yu, M.R. Weis, D.C. Lamppa, T.J. Awe, A.J. Harvey-Thompson, and R.D. McBride. The authors also gratefully acknowledge programmatic support and encouragement from K. Beckwith, L. Shulenburger, B. Jones, D. Ampleford, K. Peterson, G. Rochau, and D. Sinars.

One of the authors (Shipley) was partially funded by the National Nuclear Security Administration through the Krell Institute via the Stewardship Science Graduate Fellowship (Grant No. DE-NA0003864). One of the authors (Shipley) was also funded in part by Sandia's Laboratory Directed Research and Development program via the Truman Fellowship, Project No. 226067. Sandia National Laboratories is a multi-mission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC (NTESS), a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy's National Nuclear Security Administration (DOE/NNSA) under contract DE-NA0003525. This written work is authored by an employee of NTESS. The employee, not NTESS, owns the right, title, and interest in and to the written work and is responsible for its contents. Any subjective views or opinions that might be expressed in the written work do not necessarily represent the views of the U.S. Government. The publisher acknowledges that the U.S. Government retains a non-exclusive, paid-up, irrevocable, world-wide license to publish or reproduce the published form of this written work or allow others to do so, for U.S. Government purpose. The DOE will provide public access to results of federally sponsored research in accordance with the DOE Public Access Plan.

## AUTHOR DECLARATIONS

### Conflict of Interest

The authors have no conflicts to disclose.

### Author Contributions

**Gabriel Alan Shipley:** Conceptualization (equal); Data curation (lead); Formal analysis (lead); Investigation (lead); Methodology (lead); Writing – original draft (lead); Writing – review & editing (lead). **Daniel Edward Ruiz:** Conceptualization (supporting); Formal analysis (supporting); Methodology (supporting); Writing – review & editing (supporting). **Christopher A. Jennings:** Formal analysis (supporting); Methodology (supporting); Writing – review & editing (supporting). **David Yager-Elorriaga:** Conceptualization (supporting); Methodology (supporting); Writing – review & editing (supporting). **Paul F. Schmit:** Conceptualization (equal); Methodology (supporting); Writing – review & editing (supporting).

## DATA AVAILABILITY

The data that support the findings of this study are available within the article.