Tearing of charged current layers

Using PIC simulations, this study reveals that electrically charged current layers in astrophysical plasmas exhibit an intricate interplay between electrostatic Bernstein waves and the tearing instability, where charge redistribution modifies the initial tearing stage and can significantly enhance plasmoid growth rates depending on the layer configuration and plasma temperature.

The Gist

Imagine the universe as a giant, cosmic kitchen where invisible "current sheets" (thin layers of electrically charged gas) act like the dividers between different flavors of soup. In places like pulsar winds (the super-fast streams of particles shooting out of dead stars), these dividers can sometimes get a little "charged up," meaning they have an extra bit of electric imbalance, even though the soup itself is usually a perfect mix of positive and negative ingredients.

This paper is a scientific investigation into what happens when these cosmic dividers are charged. The researchers used powerful computer simulations to watch how these layers behave, specifically looking at how they "tear" apart and form clumps of plasma (called plasmoids). They compared two main types of cosmic dividers: the Harris sheet (a straight, flat divider) and the Rotational sheet (a swirling, twisting divider).

Here is what they found, explained through simple analogies:

1. The "Static Shock" in the Straight Divider (The Harris Sheet)

Imagine a straight line of people holding hands, but the people on the left are holding positive balloons and the people on the right are holding negative ones. In the middle, the magnetic force holding them together is zero.

• The Problem: Because the positive people are stuck at the "top" of an energy hill and the negative people are stuck at the "bottom," the system is unstable. It's like balancing a pencil on its tip.
• The Reaction: Almost immediately, the system creates a "static shock." The paper calls these Bernstein waves. Think of these as rapid, vibrating ripples of charge that bounce back and forth inside the layer, like a guitar string plucked and trapped inside a box.
• The Result: These vibrations act like a fast-acting reset button. They quickly shuffle the charges around until the layer becomes electrically neutral again.
• The Tearing: Once the static shock settles down, the layer tears apart just like it would if it had never been charged in the first place. The "clumps" (plasmoids) that form are only slightly charged.
• The Temperature Factor: If the "soup" is hotter (the particles are moving faster), these static shocks happen even faster, like a hot metal cooling down quicker than a cold one.

2. The "Swirling Storm" in the Twisting Divider (The Rotational Sheet)

Now, imagine a divider that isn't straight but is swirling like a tornado.

• The Surprise: Even if you start with a perfectly balanced, neutral swirl, the act of it tearing apart naturally creates huge, temporary surges of electric charge. It's like a calm river suddenly developing massive, chaotic whirlpools of static electricity just because the water is moving fast.
• The Speed Boost: Here is the big discovery: If you start with a charged, swirling layer, it tears apart much faster than a neutral one. It's as if adding a little bit of extra static electricity to a swirling storm makes the storm explode into pieces much more violently and quickly.
• The Heat Factor: Just like the straight divider, if the swirling layer is hotter, the charge fluctuations happen faster.

3. What This Means for Pulsars

The paper connects this to pulsars (rapidly spinning neutron stars). The famous "Michel solution" is a mathematical model describing how the current sheet around a pulsar should look.

• The Reality Check: The researchers found that this mathematical model is actually unstable. It's like drawing a perfect circle on a piece of paper that is actually made of jelly; the jelly would immediately wobble and change shape.
• The Conclusion: A perfectly charged current sheet, as described in the old math, probably never actually forms in nature. Instead, the moment it tries to form, those "static shocks" (Bernstein waves) kick in, scramble the charges, and prevent the sheet from ever reaching that perfect, charged state. The universe seems to prefer a slightly messy, neutral state over a perfectly charged, unstable one.

Summary

In short, the paper shows that when these cosmic current layers get electrically charged:

1. Straight layers quickly shake off the charge with fast vibrations before tearing apart normally.
2. Swirling layers generate massive charge surges on their own, and if they start charged, they tear apart much faster.
3. Nature likely prevents the formation of the perfectly charged models we see in textbooks because these layers are too unstable to hold that charge for long.

Technical Summary

Technical Summary: Tearing of Charged Current Layers

Problem Statement
The paper investigates the dynamics and plasmoid formation (tearing instability) in astrophysical current layers that are electrically charged. While standard tearing mode theory typically assumes charge-symmetric plasmas (e.g., $e^\pm$), the author notes that in relativistic astrophysical environments, such as pulsar winds (specifically the Michel solution), current layers may possess a net electric charge. The study addresses a gap in understanding how this net charge, and the associated electrostatic forces, interacts with the magnetic tearing mode. The author distinguishes these cases from trivial Lorentz-transformed charges, focusing on configurations where spatially distributed currents and charges cannot be locally transformed away.

Methodology
The study employs Particle-in-Cell (PIC) simulations to model two distinct types of charged current layers:

1. Charged Harris Current Sheet (CHCS): Adapted from the classic Harris configuration to include non-zero charge density while maintaining global force balance.
2. Charged Rotational Current Sheet (CRCS): Adapted from force-free, sheared rotational configurations (Lyutikov 2003).

The simulations utilize a two-fluid approach for pair plasmas ($e^\pm$). The setup involves defining a global equilibrium where electromagnetic fields, total charges, and currents are balanced, followed by a local prescription for species velocity, density, and pressure. Both single and double current layer configurations are simulated to manage boundary conditions and potential differences. The study includes parameter scans varying plasma temperature ($\Theta$), magnetic field strength ($b_0$), and initial charge levels ($\gamma_0$). Both 2D and 3D simulations are performed, alongside quasi-1D tests to isolate electrostatic relaxation from plasmoid instability.

Key Results

• Charged Harris Current Sheets (CHCS):

  • Electrostatic Instability: Initially force-balanced charged Harris sheets are found to be electrostatically unstable. The configuration places opposite charges at electric potential extrema (maxima/minima) where the magnetic field vanishes.
  • Bernstein Waves (BWs): This instability rapidly generates fast electrostatic oscillations identified as Bernstein waves. These waves are trapped within the current layer, reflecting at the upper hybrid resonance.
  • Charge Relaxation: The BWs quickly redistribute charge, causing the system to relax to a state where the net charge density is effectively zero ($\rho_e \sim 0$). This relaxation occurs on a timescale faster than the tearing instability.
  • Tearing Dynamics: Once the charge is redistributed, the subsequent tearing instability and plasmoid growth resemble the standard uncharged case. The resulting plasmoids are only mildly charged. The initial charge does not significantly alter the long-term tearing rate or overall plasmoid growth.

• Charged Rotational Current Sheets (CRCS):

  • Charge Fluctuations: Even in initially charge-neutral configurations, large fluctuations in charge density develop during the instability.
  • Enhanced Tearing Rate: In contrast to the Harris sheet, charged rotational layers exhibit a significantly increased tearing rate compared to their uncharged counterparts.
  • Temperature Dependence: The dynamics are strongly dependent on initial plasma temperature. In warmer plasmas, charge relaxation via BWs occurs much faster due to higher phase speeds. In cold cases, relaxation is slow, and BW oscillations may not complete a full cycle before other dynamics dominate.
  • 3D Consistency: 3D simulations confirm the 2D results, showing large charge density fluctuations appearing at similar timescales ($L/c \sim 20-30$).

Significance and Claims
The paper claims that the presence of net electric charge fundamentally alters the initial evolution of current layers through the generation of Bernstein waves, though the long-term impact differs by geometry:

• For Harris-type sheets (relevant to pulsar equatorial current sheets), the electrostatic instability acts as a rapid relaxation mechanism that neutralizes the layer before the tearing mode fully develops. This suggests that a "proper" charged Michel current sheet may never fully assemble in pulsar winds; instead, efficient electrostatic dissipation via BWs at the upper hybrid resonance would occur during assembly.
• For Rotational sheets, the charged state does not merely relax but actively accelerates the tearing instability, leading to faster plasmoid formation.

The author concludes that while the initial charge state triggers complex electrostatic dynamics (Bernstein waves), the ultimate plasmoid formation in Harris sheets proceeds similarly to the uncharged case, whereas rotational sheets experience a greatly enhanced tearing rate in the charged regime. The study highlights the importance of plasma temperature and the interplay between electrostatic waves and magnetic tearing in relativistic, charged astrophysical plasmas.