Ivy-Mode Channels in Dielectric Breakdown

• Ivy-mode channels are nearly one-dimensional dielectric breakdown paths in electron-irradiated PMMA that display regular sinusoidal modulations due to a z-pinch entropy-mode instability.
• High-speed imaging, SEM, and Raman mapping reveal that these channels typically have a radius of about 35 µm and a modulation wavelength near 80 µm, capturing the plasma’s transient state.
• By excluding slower mechanical and thermal processes, the analysis confirms that the observed periodicity is generated during a nanosecond plasma discharge, linking channel morphology to specific plasma parameters.

Searching arXiv for the cited paper and closely related references. Ivy-mode channels are a newly identified type of dielectric breakdown channel observed in electron-irradiated polymethyl methacrylate (PMMA), a transparent plastic used in the study as a proxy for many spacecraft insulators. After intense charging with 5 MeV electron pulses from a linear accelerator and tip-triggered breakdown, trapped charge suddenly discharges through the bulk and forms conductive trees known as Lichtenberg figures. Among the observed morphologies, ivy-mode channels are highly directional, long, nearly straight channels with minimal branching; they represent the fastest solid-state discharge observed, with propagation speeds approaching $10^7~\text{m/s}$. Their defining additional feature is a regular periodic modulation of channel width, with alternating bulges and narrower segments over millimeter lengths. The reported interpretation is that these modulations are produced during the nanosecond plasma discharge by a z-pinch entropy-mode instability, rather than by post-discharge thermal or mechanical relaxation (Schwartz et al., 18 Aug 2025).

1. Experimental setting and defining morphology

The experiments concern dielectric breakdown in electron-irradiated PMMA under conditions intended to isolate the physics of ultra-fast solid-state discharge. Blocks of PMMA are intensely charged using 5 MeV electron pulses from a linear accelerator. After charging, a sharp metallic tip triggers dielectric breakdown, and the stored charge discharges through the bulk. High-speed imaging with nanosecond resolution resolves several distinct breakdown morphologies, among which ivy-mode channels form an advancing “breakdown front”: behind the front the material is carbonized, while ahead of it the material remains unbroken and charged (Schwartz et al., 18 Aug 2025).

Morphologically, ivy-mode channels differ from more familiar tree-like Lichtenberg figures. They are nearly one-dimensional, with few branches, and retain a single dominant direction over long distances. Many breakdown channels in dielectrics are roughly uniform in width; ivy-mode channels instead often show regular, quasi-sinusoidal width variations. In ivy-mode discharges, much of the trapped charge remains unperturbed after breakdown, indicating a front-like discharge rather than wholesale decharging. This combination of extreme speed, simple geometry, and regular modulation is the basis for treating ivy-mode channels as a comparatively clean system for identifying the instability that shapes the channel during the discharge phase.

The broader significance of this morphology lies in what it records. The channel does not merely mark where the dielectric failed; it preserves a spatial signature of the plasma state inside the solid during breakdown. The paper’s central claim is therefore not only taxonomic, but mechanistic: ivy-mode channels provide access to the physics of a dense, transient, current-carrying plasma confined within an insulating solid (Schwartz et al., 18 Aug 2025).

2. Geometric periodicity and characterization methods

Optical microscopy and scanning electron microscopy (SEM) show that a typical ivy-mode channel has mean radius

$R \approx 35~\mu\text{m},$

with characteristic wavelength

$\lambda \approx 80~\mu\text{m},$

so that the dimensionless wavenumber is

$kR = \frac{2\pi R}{\lambda} \approx 2.8.$

These modulations extend over many wavelengths, and the channel interior is rough and carbonized. A surrounding boundary layer appears as a slightly darker shell in SEM and back-lit optical images, and this boundary layer itself exhibits the same periodic modulation (Schwartz et al., 18 Aug 2025).

The reported workflow for characterization combines full-block back-lit photography, white- and polarized-light optical microscopy, SEM of fractured samples, and Raman mapping. For wavelength extraction, channel width is measured as a function of axial position. Because optical resolution is limited, the width profile is estimated from carefully hand-drawn bounds and periodicity is then inferred from the positions of peaks and troughs. Over the measurement window, this procedure gives

$\lambda_{\text{geom}} = 60.1 \pm 6.1~\mu\text{m},$

which the authors regard as consistent with the order-of-magnitude $80~\mu\text{m}$ value used in instability modeling (Schwartz et al., 18 Aug 2025).

A key geometric conclusion is that the periodicity is dominated by a single mode. The channel resembles a sinusoidal radius modulation with $kR \approx 2.8$, although other modes may also be present. The paper interprets this as the usual instability pattern-selection outcome in which the fastest-growing mode sets the primary morphology. This suggests that the observed structure is not a random irregularity of carbonization, but the macroscopic imprint of a specific linear instability with a finite preferred wavelength.

3. Temporal sequence and evidence for plasma-phase formation

High-speed camera frames with 3 ns exposure and 1 ns interframe spacing show the ivy-mode front moving at about $10^7~\text{m/s}$. Over a typical channel length of centimeters, the discharge therefore completes within tens of nanoseconds. The reported sequence is: discharge initiation at the metallic tip, plasma channel formation, high-current discharge through the channel, and plasma quench. Later evolution occurs on much longer timescales: evaporated gas flows through the channels on millisecond timescales, and the hot boundary layer cools and solidifies on second timescales (Schwartz et al., 18 Aug 2025).

The article’s principal interpretive alternative is that the periodic width modulation might arise after the discharge through slower mechanical or thermal processes in the solid or in a viscoelastic boundary layer. The paper rejects that interpretation and argues that the structure forms during the nanosecond plasma phase. Four lines of evidence are given. First, Raman spectroscopy shows that the channel interior is carbon rather than PMMA; carbon formation is treated as a plasma process requiring high-temperature conditions. Second, the intensities of the carbon D and G Raman modes peak where the channel is widest and trough where it is narrowest, implying a temperature pattern coupled to the width pattern. Third, the surrounding boundary layer is itself periodically modulated, which indicates that the instability was already active while radiation from the plasma was heating the walls. Fourth, the solid-state instability candidates examined in the paper operate on millisecond-to-second timescales and fail to reproduce the wavelength with realistic parameters (Schwartz et al., 18 Aug 2025).

The Raman measurements are spatially explicit. Mapping is performed over a $75 \times 148~\mu\text{m}^2$ area across a periodic channel section. The carbon D mode and carbon G mode have nearly identical spatial patterns, with intensity maxima at wider parts of the channel and minima at narrower parts. The PMMA $\nu$ stretching mode shows weaker but qualitatively aligned structure. Two one-sided t-tests (TOST) are used to test whether intensity periodicity and width periodicity are statistically equivalent within a chosen margin. Relative to a width periodicity of $R \approx 35~\mu\text{m},$0, the minimum equivalence margins are 32.4% for the D mode, 40.2% for the G mode, and 85.3% for the PMMA mode. The tighter correspondence of D and G with width is taken as evidence that carbon deposition is directly tied to the width modulation (Schwartz et al., 18 Aug 2025).

In this reading, the geometric pattern is a frozen-in record of hotter plasma in the wider segments and cooler plasma in the narrower ones. That conclusion is important because it shifts the problem of pattern formation from post-discharge materials mechanics to plasma instability theory.

4. Competing instability mechanisms

The paper evaluates three candidate mechanisms: the Asaro–Tiller–Grinfeld (ATG) instability, the Plateau–Rayleigh instability (PRI) in a viscoelastic boundary layer, and the z-pinch entropy mode. Their distinguishing features, as summarized in the study, are the dominant physics, timescale, and whether the instability is directional along the axis (Schwartz et al., 18 Aug 2025).

Mechanism	Competing physics	Timescale / directional character
Asaro–Tiller–Grinfeld	Elasticity vs surface tension	$R \approx 35~\mu\text{m},$1 seconds; not directional
Plateau–Rayleigh	Inertia (or elasticity) vs surface tension	$R \approx 35~\mu\text{m},$2 milliseconds; directional
z-pinch entropy mode	Density vs temperature exchange in a current-carrying plasma	$R \approx 35~\mu\text{m},$3 nanoseconds; directional

For ATG, the proposed context is a thin carbon film coating the channel walls and bonded to PMMA. During cooling, mismatch in thermal expansion coefficients generates stress, and a stressed interface could in principle undulate to reduce elastic energy. The paper adapts cylindrical ATG analysis to a cylindrical cavity and asks what parameters would be required to satisfy the observed $R \approx 35~\mu\text{m},$4 with $R \approx 35~\mu\text{m},$5. Using the known PMMA surface tension $R \approx 35~\mu\text{m},$6, the numerical analysis finds that matching the observation would require elastic moduli

$R \approx 35~\mu\text{m},$7

over realistic $R \approx 35~\mu\text{m},$8. PMMA, however, has

$R \approx 35~\mu\text{m},$9

so the required modulus is five to six orders of magnitude too small. On this basis, ATG is excluded. The paper further notes that ATG acts during slow cooling and does not account for the Raman-inferred temperature pattern (Schwartz et al., 18 Aug 2025).

For PRI, the proposed setting is the post-discharge channel containing hot gases and evaporated material, surrounded by a viscoelastic or semi-liquid PMMA boundary layer produced by UV radiation. The viscoelastic analysis follows Tamim & Bostwick (2021). In the purely fluid limit, the maximum unstable dimensionless wavenumber is

$\lambda \approx 80~\mu\text{m},$0

The measured value is

$\lambda \approx 80~\mu\text{m},$1

which lies in the stable regime even before viscoelasticity is added. Because viscoelasticity and elasticity stabilize the system further, the study concludes that PRI cannot generate the observed short-wavelength modulation (Schwartz et al., 18 Aug 2025).

By elimination and by timescale matching, the z-pinch entropy mode remains the only candidate that is both fast enough and able to select a finite wavelength consistent with observation. The result is not merely that the plasma explanation survives; it is that the solid-state alternatives fail both parametrically and temporally.

5. z-pinch entropy mode and inferred plasma state

The plasma model treats the breakdown channel as a z-pinch: a column of plasma carrying axial current $\lambda \approx 80~\mu\text{m},$2 along the channel. Ampère’s law gives the azimuthal magnetic field

$\lambda \approx 80~\mu\text{m},$3

and the magnetic pressure

$\lambda \approx 80~\mu\text{m},$4

compresses the plasma toward the axis while thermal pressure resists compression. The equilibrium condition is written

$\lambda \approx 80~\mu\text{m},$5

Classical ideal-MHD sausage modes favor very short wavelengths, which is inconsistent with a sharply defined finite wavelength. The paper therefore turns to the entropy mode, a kinetic, electrostatic, flute-like instability relevant when finite-gyroradius effects are important (Schwartz et al., 18 Aug 2025).

The operative regime is

$\lambda \approx 80~\mu\text{m},$6

where the ion gyroradius is

$\lambda \approx 80~\mu\text{m},$7

In this regime, ions cannot conserve magnetic moment adiabatically, and the entropy mode becomes the dominant $\lambda \approx 80~\mu\text{m},$8 mode. The physical picture advanced in the paper is that slightly wider channel segments host higher-temperature, lower-density plasma, while slightly narrower segments host lower-temperature, higher-density plasma. Misaligned density and temperature gradients then generate entropy fluctuations that grow into a periodic modulation of radius and temperature. This is presented as the plasma counterpart of the Raman observation that wider segments are more strongly carbonized (Schwartz et al., 18 Aug 2025).

To match the observed wavelength, the study numerically solves the entropy-mode dispersion relation cited from Angus et al. (2019), using $\lambda \approx 80~\mu\text{m},$9, observed $kR = \frac{2\pi R}{\lambda} \approx 2.8.$0, and an ion mixture based on PMMA chemistry: C$kR = \frac{2\pi R}{\lambda} \approx 2.8.$1, O$kR = \frac{2\pi R}{\lambda} \approx 2.8.$2, and H$kR = \frac{2\pi R}{\lambda} \approx 2.8.$3 in proportions matching C$kR = \frac{2\pi R}{\lambda} \approx 2.8.$4O$kR = \frac{2\pi R}{\lambda} \approx 2.8.$5H$kR = \frac{2\pi R}{\lambda} \approx 2.8.$6. Two scenarios are considered: fixed current with varying ion composition, and fixed PMMA stoichiometry with current varied from 100 to 1000 A. Most solutions that reproduce $kR = \frac{2\pi R}{\lambda} \approx 2.8.$7 lie in the regime $kR = \frac{2\pi R}{\lambda} \approx 2.8.$8. The inferred plasma conditions consistent with the observed wavelength are electron temperature

$kR = \frac{2\pi R}{\lambda} \approx 2.8.$9

electron density about $\lambda_{\text{geom}} = 60.1 \pm 6.1~\mu\text{m},$0 of solid density, and channel current

$\lambda_{\text{geom}} = 60.1 \pm 6.1~\mu\text{m},$1

per channel (Schwartz et al., 18 Aug 2025).

These values are checked against direct current measurements. Using highly charged PMMA samples at 1.5 $\lambda_{\text{geom}} = 60.1 \pm 6.1~\mu\text{m},$2C/m$\lambda_{\text{geom}} = 60.1 \pm 6.1~\mu\text{m},$3 electron fluence and experimental conditions that produce isolated one-dimensional channels, the study measures single-channel currents rather than the sum over many branches. For samples with breakdown velocity above the ivy-mode threshold,

$\lambda_{\text{geom}} = 60.1 \pm 6.1~\mu\text{m},$4

the inferred typical peak current is

$\lambda_{\text{geom}} = 60.1 \pm 6.1~\mu\text{m},$5

That measured value lies within the predicted $\lambda_{\text{geom}} = 60.1 \pm 6.1~\mu\text{m},$6 range, and the paper treats this agreement, together with the morphology and Raman results, as validation of the entropy-mode interpretation (Schwartz et al., 18 Aug 2025).

6. Predictive use, significance, and unresolved questions

From the combined modeling and measurements, the study infers a characteristic ivy-mode plasma state during the nanosecond discharge: channel radius $\lambda_{\text{geom}} = 60.1 \pm 6.1~\mu\text{m},$7, current per channel $\lambda_{\text{geom}} = 60.1 \pm 6.1~\mu\text{m},$8, electron temperature $\lambda_{\text{geom}} = 60.1 \pm 6.1~\mu\text{m},$9, electron density $80~\mu\text{m}$0, and an ion composition consistent with PMMA decomposition into C, O, and H species. These parameters place the discharge in a regime where the ion gyroradius is small but comparable to the perturbation scale and where $80~\mu\text{m}$1, so the entropy mode is expected to be the dominant $80~\mu\text{m}$2 instability (Schwartz et al., 18 Aug 2025).

On that basis, the paper proposes a framework for predicting discharge morphology. The procedure is: estimate or measure channel radius $80~\mu\text{m}$3, current per channel $80~\mu\text{m}$4, likely plasma density $80~\mu\text{m}$5 and temperature $80~\mu\text{m}$6, and ion composition from material chemistry; apply the entropy-mode dispersion relation

$80~\mu\text{m}$7

to determine the most unstable wavenumber $80~\mu\text{m}$8; then predict

$80~\mu\text{m}$9

while checking the rough onset condition

$kR \approx 2.8$0

If $kR \approx 2.8$1 lies in an observable range such as 1–10, the expectation is visible periodic modulation; if channel radius is too large, current per channel too small, or the plasma too cold or too dense, the drive may be weak and the channels may remain smooth (Schwartz et al., 18 Aug 2025).

The significance of this framework extends beyond morphology description. The periodic modulation is presented as a fingerprint of the entropy mode and therefore as a post-mortem diagnostic of plasma state during dielectric breakdown. Measuring $kR \approx 2.8$2 and $kR \approx 2.8$3 in breakdown channels may allow inference of discharge current, plasma temperature, and density through the instability model. For space electronics, the reported implication is that breakdown damage is shaped by plasma parameters, and that altering current density or channel-radius selection could in principle suppress or modify the instability and the resulting damage pattern (Schwartz et al., 18 Aug 2025).

The same paper also states several limitations. The work is confined to a single material system, PMMA. Plasma density and temperature are inferred indirectly rather than measured directly during the nanosecond discharge. The instability modeling imports drift-ideal MHD and gyro-kinetic theory formulated for classical z-pinches into a confined solid-state channel with complex boundary conditions. The mechanism that selects channel radius remains unresolved, and not all ivy-mode channels are periodic, indicating a threshold effect whose exact form has not yet been determined. The amorphous carbon is confirmed by Raman spectroscopy, but its detailed nanostructure and formation pathway are not fully characterized (Schwartz et al., 18 Aug 2025).

Taken together, these results define ivy-mode channels as ultra-fast, nearly one-dimensional dielectric breakdown paths in electron-irradiated insulators whose periodic width modulation is produced during the nanosecond discharge by a z-pinch entropy-mode plasma instability. A plausible implication is that channel morphology in irradiated insulators can serve not only as a failure signature but also as a quantitative record of the plasma regime that generated the breakdown.