Laser and Astrophysical Plasmas and Analogy between Similar Instabilities

: Multipulse laser–matter interactions initiate nonlinear and nonequilibrium plasma fluid flow dynamics and their instability creating microscale vortex filaments, loop-soliton chains, and helically paired structures, similar to those at the astrophysical mega scale. We show that the equation with the Hasimoto structure describes both, the creation of loop solitons by torsion of vortex filaments and the creation of solitons by helical winding of magnetic field lines in the Crab Nebula. Our experiments demonstrate that the breakup of the loop solitons creates vortex rings with (i) quasistatic toroidal Kelvin waves and (ii) parametric oscillatory modes—i.e., with the hierarchical instability order. For the first time, we show that the same hierarchical instability at the micro-and the megascale establishes the conceptual frame for their unique classification based on the hierarchical order of Bessel functions. Present findings reveal that conditions created in the laser-target regions of a high filament density lead to their collective behavior and formation of helically paired and filament-braided “complexes”. We also show, for the first time, that morphological and topological characteristics of the filament-bundle “complexes” with the loop solitons indicate the analogy between similar laser-induced plasma instabilities and those of the Crab and Double-Helix Nebulas—thus enabling conceptualization of fundamental characteristics. These results reveal that the same rotating metric accommodates the complexity of the instabilities of helical filaments, vortex rings, and filament jets in the plasmatic micro-and megascale astrophysical objects


Introduction
Many physical problems-like the plasma system dynamics and instability-are complex in nature and include various physical processes, taking place on a wide range of space and time scales.Their effects on plasma dynamics occur at the time scales from nanoseconds to several millennia and at the spatial scales from less than a micron to those measured in the light years.As mechanisms occur on different scales to the patterns they shape, this makes scale the key conceptual problem in the evolution of shock-induced flow in the plasmatic systems [1][2][3][4][5].Having established a relationship between the physical processes that govern plasma dynamics and the evolution of instabilities at the laboratory microscale, one can establish an analogy with the processes at a much larger scale of astrophysical objects, like Nebulas.Reflecting on experiences from studying laser-matter interactions (LMIs) and plasmatic systems at various scales, it is possible to conclude that there are some common features (a common "skeleton" in the multidisciplinary research-relating to multiphysics phenomena) [5], but that, each system will have its specific differences.
In this paper, we aim to give an overview of the main characteristics of plasma systems created by intense and ultraintense lasers, which cause plasma instabilities and structures generated by the shock-induced variable density flow and by the magnetohydrodynamicinduced instabilities.
The first part of this article summarizes features with regard to the (i) main characteristics of laser-generated plasmas by intense lasers (≥10 8 W/cm 2 ) and ultraintense lasers (≥10 15 W/cm 2 ) [2][3][4]6]; (ii) the physical processes that can be initiated at these power densities; and (iii) characteristics of laser-generated plasma instabilities, which are analyzed having in mind that similar phenomena occur in the astrophysical systems like stars and nebulas.Along this line of thinking, the conceptualization of fundamental characteristics of laser-induced plasma instabilities is extended and applied to those observed in the Crab Nebula and Double-Helix Nebula.

Short Outlines of Plasmas Generated by Intense and Ultraintense Lasers 2.1. Plasmas Created by Intense Lasers
Plasma formation in laser-matter interaction (LMI) is initiated when the laser power density reaches the optical breakdown threshold, which depends on the laser pulse (energy, wavelength, and pulse duration) and the target material (optical, thermal, and electrical properties and surface morphology).The early stages of plasma evolution created by intense power lasers of the pulse duration τ~10 −9 s and power density ~10 8 -10 11 W/cm 2 are characterized by the high electron density (~10 19 cm −3 ) and temperature ~10 4 K [7,8], which may reach ~10 5 K, and dominated by the Bremsstrahlung continuum emission [9].Both the electron density and temperature decrease with time due to plasma spheroid expansion and recombination phenomena.High collision frequency between the free electrons and heavy species (atoms and ions) leads to excitation/de-excitation and ionization/recombination equilibria [9].
Regarding this process, the excited ions and atoms in the laser-induced plasma soon de-excite by emitting radiation; atomic line emission forms the "laser-induced breakdown signal", which enables spectroscopic diagnostics.The line emission-superimposed on the continuum emission-can be observed after the plasma expansion associated with the fast decrease in temperature [9].Plasma spheroid evolution/expansion depends on the energy deposited into the target, as well as on the environmental conditions (vacuum or background gas).
In the vacuum environment, plasma spheroid evolution is adiabatic and depends on the laser pulse duration; it is therefore different for the ns, ps, and fs laser pulses.According to [10], laser pulses longer than about 5ps cause the laser-plasma interaction and plasma heating with the temperature increasing with pulse duration, which is well illustrated for the ns pulses."Due to the fact that the electron-lattice heating time is ~10 −12 s, which is ~1000 times shorter than the pulse duration-the ns-laser pulse causes domination of thermal effects over the ionization.The laser pulse melts and vaporizes target surface, while the temperature increase ionizes the atoms.The plasma spheroid becomes opaque for laser radiation between 10 −9 s and 10 −8 s, so that the tail of the laser pulse interacts with the plasma spheroid becoming absorbed.The absorbed nanosecond laser pulse causes plasma reheating, elongates its lifetime, increases line emission, but also increases the background emission" [10].
In the background gas environment, plasma spheroid evolution is non-adiabatic and depends on ablated mass, spot size, energy deposited into the target, and the pressure of background gas.Plasma creation and evolution in the background gas are shown in Figure 1.
Characteristics of plasma evolution also depend on the laser pulse duration and the wavelength.Considering the ns-LMIs and the wavelengths ranging from the far infrared CO 2 laser (λ = 10.8 µm, τ = 100 ns), to the infrared Nd-YAG laser (λ = 1.06 µm, τ = 40 ns), to the visible Ruby laser (τ = 30 ns, λ = 600 nm),and to the ultraviolet XeCl laser (λ = 308 nm, τ = 16-20 ns), one finds that in all cases plasma expansion depends on the interaction with the surrounding background gas [11,12].The expanding plasma becomes a mixture of atoms and ions of both, vaporized target material (In, Co, Fe, Ti, etc.), and ambient gas.The reaction with components of the background gas (usually air) like O 2 and CO 2 generates metal-oxide molecules, etc., as ionized species.Its lifetime-depending on density, the ambient gas, and the laser wavelength-ranges from ~300 ns to >40 µs [9].For dependence of the plasma lifetime on the wavelength, it was found that shortwavelength UV irradiation of the iron target creates plasma, which is initially dominated by the continuum background up to ~400 ns; then, the atomic line emission appears and extends to the µs time scale.The longer wavelength of an IR laser causes the continuum emission at longer times ranging to several microseconds [9].The evolution of the line emission intensity for two iron lines (Fe I 285.2 nm and Fe II 288.4 nm) shows characteristic behavior.During the first 3 µs, the intensity of the ionic line exceeds the neutral atoms by 50 times, whereas after ~10 µs, the neutral line intensity increases and becomes eight times that of the ionic [9].
The dependence of the plasma lifetime on the pressure of the background gas can be related to the expansion velocity-as observed from the ICCD (intensity field chargecoupled device) photographs of visible emission and the size of the aluminum plasma cloud.For the pressure of p = 100 Tor~1.33 × 10 4 Pa [13], it was found that the initial expansion velocity of ~3000-4000 m/s-(between 450 ns and 950 ns after the plasma formation)starts to decrease after ~1500 ns to ~900-1000 m/s.This continues to ~2000 ns (2 µs) when the light emission ends.However, for the atmospheric pressure p = 10 5 Pa, the expansion velocity is lower while the plasma lifetime is extended for a few orders of magnitude.For the giant laser pulses, the plasma lifetime extends to ~0.2 ms [14].In general, characteristics of laser-induced plasma change with time and with the spheroid expansion, between ns and µs time scale, when the plasma spheroid extends over the target surface [10].
Expansion of the plasma spheroid into the surrounding background gas forms the air shell [15] or the cocoon [16-19] (Figure 2a).With the expansion, its density decreases in the central region (CR) and increases at the periphery of the spot (in the NCR and NPR), resulting in the inversion of the initial density profile.Radial differences in the plasma Atoms 2024, 12, 23 4 of 41 density may become very large.The number density of the light-fluid particles in the NCR is higher for about two orders of magnitude with respect to the CR as shown in [15,20].Figure 2b shows the cross-section of the expanding plasma spheroid in the (z,x) plane.
.  Plasma-generated shock with "blast profile" or simply called blast wave "is a shock wave that consists of a leading shock front followed by expansion wave which decays its strength at a steady rate.This decay occurs across all the properties behind the leading shock.For a typical shock wave with 'step profile' or 'shock profile', the properties behind the shock remain constant" [21].The shock wave striking density interface of the vapor/plasma (low-density fluid, (ρL) and molten surface layer (high-density fluid, (ρH)) impulsively accelerates the light fluid into the heavy fluid causing the evolution of instabilities due to the transfer of mass, momentum, and energy to the density interface.
The pressure shock generated in the confined (or semiconfined) configuration may be estimated as follows [22].Assuming the laser power density of ~1 GW/cm 2 (which is higher than the plasma breakdown threshold of various metals of ~0.36 GW/cm 2 according to [23]), the plasma detonation generates a shock wave with the pressure peak p(kbar), given by the Raizer equation [22,24,25].
where Z is the reduced shock impedance of the two materials, I is the flux of the laser beam (GW/cm 2 ), and a is the corrective factor, a ≈ 0.20, and Z~390,625 (~4 × 10 5 ) [22].Estimation of the pressure generated in the confined configuration for the laser power density, I ~1 GW/cm 2 , and a and Z parameters similar to Devaux et al. [22], Equation ( 1) is as follows [17]: p~15.3 kbar~1.53 GPa. ( Detonation of plasma discoid in a spatially localized laser-heated region causes the overpressure within a background gas.However, the overpressure quickly diminishes, and the vapor/plasma blast wave is highly accelerated [26]. Impulsive fluid acceleration by the shock wave causes Richtmyer-Meshkov instability (RMI) as a vorticity-driven phenomenon, which starts with linear growth as small amplitude perturbations.It evolves with time when the incident shock wave propagates from a light to a heavy gas (∇ p .∇ρ < 0) or from a heavy to a light fluid (∇ p .∇ρ > 0), associated with phase reversal [27][28][29][30][31][32].The instability causes perturbation on the interface to grow, producing vortical structures followed by the generation of bubbles and spikes.According to Schiling and Jacobs [32], «after passage of a shock wave-the interface begins to decelerate.The pressure behind the shock front decreases monotonically with distance, so that the reversal of pressure and of density gradients occurs from the heavy into the light fluid».
Rayleigh-Taylor instability (RTI) is initiated with the constant acceleration by the shock wave of the heavy fluid into the lighter one, when the density gradient is mismatched with the direction of the acceleration, i.e., (∇p .∇ρ< 0) .The interaction of a pressure gradient (∇p, from shock) with a misaligned density gradient (∇ρ, from the interface) causes the evolution of the vorticity caused by the baroclinic torque at the interface of the fluid with the velocity υ.The deposition of baroclinic vorticity on the density interface can be described by the dynamics of the vorticity field (ω = ∇ × υ) [28,[31][32][33][34].The baroclinic vorticity deposition by the shock wave is governed by the 3D vorticity equation: where D/Dt = ∂/∂t + υ∇, and (∇ρ × ∇p)/ρ 2 is the baroclinic term, ω • ∇υ is vortex stretching and ω(∇υ) is vortex dilatation [31][32][33].The amount of vorticity deposition increases with the increase in density gradients (Atwood number), pressure gradients (shock strength or Mach number, Ma), and the mismatch angle α between the pressure and density gradients (as sinα) [27].
The perturbations of the density interface start to grow into spikes of heavy fluid into light fluid and bubbles of light fluid into heavy fluid, generating the Rayleigh-Taylor instability [33][34][35].The subsequent shear flow along the growing spikes causes Kelvin-Helmholtz (KH) instability and the formation of mushroom-shaped caps at the RTI spikes [17][18][19]36].The key parameter governing the instability is the shock Mach number, Ma, which determines the effects of compressibility in the flow.It means that for weak shocks, the fluid dynamics may be described as assuming incompressible flow [32].
In many studies, the assumption is that the shock-induced fluid acceleration is stationary in space-time; in that case, the energy transfer is Kolmogorov direct cascade characteristic for the isotropic homogenous turbulence.However, in a more general case, the acceleration may be variable (nonstationary), leading to different RMI and RTI dynamics and non-Kolmogorov inverse cascade [37][38][39][40].The variable acceleration is assumed to be a power law function of time, g = Gt a , t > 0, where a is the exponent of an acceleration power law, a ∈ (−∞, +∞), and G is the pre-factor, G > 0, and their dimensions are [G] = m/s 2+a and [a] = 1.The evolution of instability is considered on the fluid layer with the periodic flow in the plane (x,y) normal to the z direction of the acceleration g; |g| = g, and also the acceleration directed from the heavy to the light fluid, g = (0,0,−g) along the z-axis [39][40][41].The evolution of RMI/RTI starts from infinitesimal spatial perturbations at the density interface, governed by parameters like the initial perturbation (the height of the disturbance h and the wavelength, Λ of parallel scratches on modulated target surface), the Atwood number A, A = (ρ H − ρ L )/(ρ H + ρ L ), and the momentum transfer.
Initial perturbation conditions: The initial perturbation-which on metal targets is deposited by the set of parallel scratches-creates multimodal perturbation.Becoming the seed of interface instability, this perturbation initiates the growth of RMI/RTI at early times and depends on the amplitude and wavelength of the perturbation modes.The multimode perturbation is the basis of the interference concept of the RMI and RTI evolution.This concept-besides the amplitude and the wavelength-also takes into account the phase relations between the modes that lead to interference of the initial perturbation modes [42,43].
The laser-generated shock wave moving on the scratched target surface becomes dispersed into a fast oscillatory wave and a slow modulation wave [44].Thus, the shock and reshock waves in the fluid layer become oscillatory waves with variations in the phase and amplitude, wavelength and frequency, as well as the direction of propagation.When such oscillatory shock wave strikes the density interface, it causes fluid acceleration and shock-induced variable density flow with the formation of RTI as well as the shear layer roll-up into vortex filaments if the Re number is Re > Re crit .

Plasmas Created by Ultraintense Power Lasers
Ultraintense power lasers create plasmas with extremely high electron density and temperature as well as strong magnetic and electric fields.A high degree of ionization, high velocities, high density, and pressure in such plasmas are comparable to that in the laserinduced inertial confinement fusion and in the stars.This makes it possible to establish an inter-relationship between laboratory, magnetic fusion, inertial fusion plasma, space, and astrophysical experiments [2][3][4][5]45], as the high energy density phenomena.In this respect, it should be noticed that "a significant fraction of the visible Universe is composed by matter in extreme conditions of temperature, density and pressure.When the pressure in a physical system exceeds 1 Mbar, this is defined as a high energy density (HED) state, which corresponds to the pressure required to deform the water molecule or in other words the pressure at which water becomes compressible, corresponding to an energy density exceeding 10 11 J/m 3 . .." [46].
Creating such plasmas, lasers of the power density from 10 18 W/cm 2 [46] up to 10 24 W/cm 2 , with the pulse duration of τ < 10 −9 s (ns)-10 −12 s (ps) and to 10 −15 s (fs) make possible the study of matter under extreme conditions by monitoring the light emission, the emission of the X-rays and plasma debris.In the plasmas created by intensities of ~10 14 W/cm 2 , laser fields start to compete with intra-atomic fields, causing rapid ionization and complex dynamics of the electron wave function [6].An example is lithium plasma created by the ultrashort pulse of τ = 10 −12 s, of 2 × 10 15 W/cm 2 , with the electron temperature of 35 eV (4 × 10 5 K) and high electron density of ~10 21 cm −3 showing the soft X-ray emission between λ = 1 nm and 15 nm [47].A linear increase in the Li III Ly α -and Ly β -line intensity coincides with the plasma expansion length.For the higher intensity of ~10 16 W/cm 2 , the laser fields surpass the intra-atomic fields that band electrons and enable rapid ionization of various targets.At even higher intensities of ~10 23 W/cm 2 , plasma manifests the new phenomenon: individual incoherent emission (radiation damping) of electrons which start to affect the plasma dynamics suppressing the plasma instabilities.This causes an energy transformation into the gamma range and the emission of gamma radiation [6].
Dynamics of plasmas created by ultraintense lasers are intended to mimic plasma dynamics in cosmic explosions and planetary cores creating ultrahigh-pressure shocks [47] In this respect, irradiation of a small sphere by an ultraintense laser creates intense shock waves and ultrahigh pressures reaching ~10 12 atmospheres.Such gigantic shock fronts are similar to the thin shock regions at the boundary between a collapsed supernova and the surrounding material-creating a spheroid of super-hot plasma.Its hot steep front is associated with a turbulent magnetic field, which affects plasma behavior [48].
Magnetic fields in laser plasma are self-generated fields caused by direct laser acceleration [49] and driven by intense and ultraintense laser pulses.They vary in strength, topology, and timescale over many orders of magnitude [50].Laser-generated magnetized plasma can be scaled to astrophysical plasmas.A 0.5-2 kJ laser with a pulse duration of 1 ns creates a laser plasma and fast electrons that induce a current in the coil generating a magnetic field [50].Their strength ranges from 50-200 T [51,52], to 600-800 T [53,54], and even to the tens of kilotesla [46,55].
A kilotesla strength magnetic field was generated by fs ultraintense laser pulses of P = 65 TW and I ≈ 10 20 W/cm 2 [56].The strength of such a magnetic field can be estimated by considering the absorption rate of the laser beam: where c is the speed of light, and τ is the pulse duration; they used the model valid for low-order orbital angular momentum (OAM) modes, characterized by the parameters ℓ and σ z .Parameter ℓ relates to the order of azimuthal modes, while σ z is the spin number; σ z = [±1] denotes a circularly polarized beam, and σ z = 0 denotes a linearly polarized beam.For the fs pulse and multi-terawatt laser power and given electron plasma density n e , the peak axial magnetic field strength of the modes (σ z = 1, l = 0) and (σ z = 0, l = 1) is given by the relationship [56]: For the laser plasma parameters n e = 0.03 cr , beam wavelength λ = 1 µm, beam waist w 0 = 6 µm, a abs = 0.435 mm −1 , τ = 100 fs, and p = 65 TW, where n c is the critical plasma density, Equation (5) gives the peak of the axial magnetic field strength of B = 1 kT.This kilotesla, multipicosecond, axial magnetic field was found to extend over hundreds of microns in underdense plasma [56].
In astrophysically relevant laser experiments, such strong magnetic fields are measured by proton imaging, in which high-energy protons interact with the electromagnetic fields, creating spatial variations in the proton fluence on the diagnostic plane [48].The collisionless interaction of two counter-propagating, laser plasmas "results in the formation of small-scale magnetic filaments.This process created in the laboratory may explain the presence of magnetic fields throughout the intergalactic medium" [48].Proton images that indicate the amplitude of the field and its orientation relative to a proton's trajectory display features corresponding to the magnetized filaments and reveal randomized systems of filamentary magnetic fields [48].Magnetic fields in laser plasmas are found universally in all plasma regimes from astrophysical to fusion systems [48,50].While the magnetic fields govern the dynamics of plasma expansion and the evolution of plasma instabilities, very strong magnetic fields enhance the generation of ion beams; they affect the formation and propagation of plasma filament jets in nebulas.

Short Outlines of the Experiments
The experiments-shortly described below-have shown that various plasma flow instabilities can be created by laser-matter interaction, which on the metal targets stay permanently after pulse termination.Multipulse laser irradiation of metal surface modulated by parallel micro-scale scratches has been used in order to induce initial multimodal perturbation [12,[57][58][59] Interference of the initial perturbation modes results in the evolution of various plasma flow instabilities on pure metal targets and targets coated by layers of different metal types.The Co-coated steel plate targets of 1 × 1 × 0.05 cm were exposed (in the open configuration) to the beam of XeCl excimer laser (λ = 308 nm, of E = 250-300 mJ, Es ~5 J/cm 2 , I~2.4 × 10 8 W/cm 2 , τ ≤ 20 ns), of the rectangular cross-section with constant power density distribution (top hat profile), and irradiated by N = 12 pulses at the low repetition frequency of 10 Hz.Laser irradiation of the Co-coated steel plate from above in the air at the atmospheric pressure (p = 10 5 Pa) generated a high-temperature plasma spheroid, which expanded laterally along the "cold" target surface as described above, establishing the surface shear layer.Nonlinear and nonequilibrium processes at the power density scale of ≥10 8 W/cm 2 take place in the shear layer of the plasma/molten-target surface.Driven by the shock-induced variable density flow, these processes lead to the evolution of the Rayleigh-Taylor instability, vortex filaments, and filament loops.This is especially emphasized when the interaction occurs with multimodal surface perturbation when the interference of initial perturbation modes causes various flow instabilities and the formation of a number of fluid flow structures.Their relaxation time is longer than the liquid-layer solidification time, so that they solidify ultrafast after pulse termination and become "imprinted" in the surface morphology.The cooling interval of ~5 ns or longer is comparable to the heating time, while the cooling rate is ~10 9 -10 10 K/s [12].The structures-permanently solidified-make possible a posteriori analysis.Thus, the target surface represents a "diagnostic plane" for the interface dynamics.The structures "encoded" in the surface morphology represent time integrated picture of hydrodynamic instabilities.A posteriori optical, SEM, and AFM analysis enables the identification of characteristic RMI and RTI structures, the vortex filaments and jets, their wavelength, and amplitude, as well as their organization and symmetry.

Vortex Filament Formation in Laser Plasma
The shock wave acceleration of the shear layer (variable density flow) establishes the RT environment and causes the rollup of shear waves into a series of parallel micronscale vortex filaments and wave-vortex structures.Their formation on the target surface was promoted by the parallel micron-scale scratches (with the separation distance Λ, which causes transversal perturbation of the shock wave as Ψ = Aexp(k p x), where A is the perturbation amplitude equivalent to the height of the scratch wall, and k p is the perturbation wave vector; k p = 2π/Λ(µm −1 ).An oscillatory shock perturbs the shear layer transversally to the y-direction, causing the formation of waves and their roll-up into vortex filaments; their length reaches L ≳ 500 µm, while the core thickness is about σ ~5-7 µm [57-62] (Figure 3).Surface patterning and multipulse LMI create a series of shock waves of different wavevectors, phases, and intensities.Two kinds of vortex filament arrays are formed; the low-density one, in which the filament-filament distance, Λ is much larger than the core size σ , and a high-density array . The parameter Λ is a crucial parameter that determines whether the filaments in the array will manifest individual dynamics after N pulses, or they will interact and form secondary coherent structures low-density one, in which the filament-filament distance, Λ is much larger than the core size σ (Λ >> σ), and a high-density array (Λ ≤ σ).The parameter Λ is a crucial parameter that determines whether the filaments in the array will manifest individual dynamics after N pulses, or they will interact and form secondary coherent structures like thick-vortex filaments, ribbons, vortex filament-bundles, and braided and tangled structures [61].

Individual Dynamics of Vortex Filaments
Vortex filaments formed with a large separation distance, Λ ≳ 10 µm (a 1D vortexfilament lattice), extend in the y-direction on the target surface and manifest individual behavior.Their formation on the surface shear layer can be described by two dynamically important quantities; the square of the vorticity modulus [63] is as follows: and the square of the rate of strain tensor [64] is as follows: The latter quantity is linked to the local energy dissipation ε = νT 2 ρ, where ρ is the fluid density, and ν is the kinematic viscosity.Another quantity connected with the above two is the pressure field p, which can be written [63] as follows: The regions of the shear layer with large ω 2 and small T 2 become rolled up into vortex filaments.Namely, the horizontal shock-momentum component causes acceleration of the shear layer, while the vertical one causes Kelvin-Helmholtz instability and rollup into axisymmetric vortex filaments when the Reynolds number Re ≥ Re critical [63].
The velocity field of the vortex filament in the cylindrical coordinates, U = U(u r , u ν, u Θ ), has the components: u z = γz (γ = strain), and u Θ = Γ/2πr(1 − exp(−γr 2 /4ν)) where ν = kinematic viscosity and Γ is the circulation.The Gaussian vorticity distribution across the filament is [63] as follows: The fluid velocity U is comparable with the velocity of the expelled vapor of ~10 3 -10 4 m/s.Assuming U ~5 × 10 3 m/s, the kinematic viscosity of liquid metal, ν~10 −6 m 2 /s, and the thickness h of the molten layer is h ~(1-2) × 10 −6 m, one finds the Reynolds number Re = Uh/ν ≈ 5 × 10 3 to be 10 4 .At some places, a plasma cloud of ~10 4 -10 5 K is created above the liquefied-boiling surface (T B ~2-5 × 10 3 K), associated with high fluid velocity, and the critical Re number for the formation of vortex filaments, estimated to be Re crit ~10 4 .For the filament of the core size, σ (σ = 2r 0 ~5-7 µm), the circulation, Γ = 2πσU, is Γ ~0.3 m 2 /s.After formation, the vortex filament immersed in the background turbulent field starts to move on the surface.

Generation of Loop Solitons on Vortex Filaments
Motion of a vortex filament of the length L and the core diameter d immersed in an external flow U ext (background turbulent shearing field) at any point is affected by its shape along its entire length.However, it is reasonable to expect that the contribution to the vortex velocity from distant parts of the filament is small compared to the local contributions [64].The curvature of the vortex filament exposed to torsion due to interaction with a turbulent shearing field creates loops called the Hasimoto solitons [57,59,64] (Figure 4a,b).
If the diameter d of the filament does not change along the length L, the condition L/d = constant makes it possible to describe its motion by the localized induction approx-imation (LIA).The vortex filament experiences a self-induced velocity described by the equation of motion, i.e., by the LIE written as [57,59,64] follows: where k is the curvature of the filament curve, b the unit binormal vector, Γ the circulation of the vortex, and e x is the unit vector in the x-direction.
velocity, and the critical Re number for the formation of vortex filaments, estimated to be Recrit~10 4 .For the filament of the core size, σ (σ = 2r0~5-7 μm), the circulation, , is Γ~0.3 m 2 /s.After formation, the vortex filament immersed in the background turbulent field starts to move on the surface.

Generation of Loop Solitons on Vortex Filaments
Motion of a vortex filament of the length L and the core diameter d immersed in an external flow Uext (background turbulent shearing field) at any point is affected by its shape along its entire length.However, it is reasonable to expect that the contribution to the vortex velocity from distant parts of the filament is small compared to the local contributions [64].The curvature of the vortex filament exposed to torsion due to interaction with a turbulent shearing field creates loops called the Hasimoto solitons [57,59,64] (Figure 4a,b).If the diameter d of the filament does not change along the length L, the condition L/d = constant makes it possible to describe its motion by the localized induction approximation (LIA).The vortex filament experiences a self-induced velocity described by the equation of motion, i.e., by the LIE written as [57,59,64] follows: The loop solitons may travel without friction along the filament and appear in LMI and in other systems including those at the astrophysical scale.For the loop soliton of the radius R, with d < R < L, the dynamics on filament can be described on the basis of LIA [57][58][59].Neglecting the background flow, the nonlinear perturbations can be described by a coupled system of Betchov-Da Rios equations, written in terms of local curvature k and torsion T. Using the complex function, Ψ (the Hasimoto transformation) is as follows: These equations are transformed into the non-linear Schrödinger equation: where s is the arch length along the filament and A is an arbitrary function of time t.For a single-loop soliton, the curvature k may be written [57][58][59] as follows: where η is constant, the torsion T and the velocity of the curvature (soliton loop) are constant.
The motion of the soliton loop along the filament can be studied as the traveling wave solution.For the reason of simplicity, we look for the traveling wave solution in 2D.Using a modified Korteveg de Vries (mKdV) equation, one finds a 2D solution for a planar curve [65,66].The mKdV equation for the planar curvature k can be written as follows: with boundary conditions k → 0, s → ±∞ .The solution of the equation can be constructed by means of the symbolic programs based on the Hirota method [66,67] for the transformation of the mKdV equation into bilinear form, which is a complex procedure and out of the scope of this paper.Following the procedure of the Hirota method, one finds the curvature k of the loop in 2D [66] k = [2k exp(ωt + ks)]/[exp(2ωt) + exp(2ks)] = 2k sec Θ where Θ = ks − ωt, and dispersion ω = k 3 .After finding the radius of the loop, the solution of the equation yields the loop soliton as real and imaginary parts of the radius vector of the loop [66].
The loop soliton propagates with increasing time; its size decreases as the wave vector increases: Therefore, a smaller loop travels faster than a larger loop [66].

Loop Soliton Chain: Multisolitons
Multipulse laser-matter interactions with patterned target surfaces create a series of such loops-the loop-soliton chains on vortex filaments [57][58][59], as well as supercomplex 3D networks [61].The loop-soliton chain-such as the one generated on the Co-coated steel, may be compared with the chain of particles-a 1D periodic lattice (Figure 5).3D networks [61].The loop-soliton chain-such as the one generated on the Co-coated steel, may be compared with the chain of particles-a 1D periodic lattice (Figure 5).κ [66].The solution of the two loop solitons is obtained by the same Hirota transformation of the mKdV equation into bilinear form, as above, and the result is presented in Figure 6.This figure shows the propagation of two solitons with  κ are obtained by very complex procedures, detailed in [66].
When a large number of almost parallel vortex filaments is generated, the loop sol iton chains can be formed simultaneously on such an array.A quasi-linear array of par allel vortex filaments (which are slightly disturbed and make quasi-periodic 1D lattices is created by multipulse LMI with Co-coated steel.In such an array of parallel vortex For a simple illustration of the multisoliton dynamics, we present the motion of 2D loop solitons of different radii r 1 and r 2 , with the curvatures k 1 and k 2 [66].The solution of the two loop solitons is obtained by the same Hirota transformation of the mKdV equation into bilinear form, as above, and the result is presented in Figure 6.This figure shows the propagation of two solitons with k 1 = 3 and with k 2 = 1 in nondimensional time.Rather long expressions for k 1 and k 2 are obtained by very complex procedures, detailed in [66]. When a large number of almost parallel vortex filaments is generated, the loop soliton chains can be formed simultaneously on such an array.A quasi-linear array of parallel vortex filaments (which are slightly disturbed and make quasi-periodic 1D lattices) is created by multipulse LMI with Co-coated steel.In such an array of parallel vortex filaments, the multipulse LMI usually forms an inhomogeneous turbulent field, in which the torsion of vortex filaments T may be small or large and may vary with the series of pulses.As a consequence, the solitons of different radii and curvatures k 1 , k 2 , . . .will be formed on the filaments.In the oscillatory strain field associated with multipulse LMI, the situation may be created in which torsion on some filament causes (+) and (-) curvature and the loop solitons, which may cancel each other or may form complex loop-soliton chains.
Such a complex case of parallel loop-soliton chains on vortex filaments is shown in Figure 7.The radius of the loop solitons is R~25-30 μm, i.e., they are mostly of the same size meaning that they travel almost at the same velocity.However, a slight difference in th loop size and curvature causes the solitons of smaller radius to travel faster than th larger ones.As a result, they approach their neighbor solitons and-an initially periodi 1D loop soliton lattice-becomes quasiperiodic (Figure 7).
Loop solitons and multiple loop solitons generated by laser plasma at the microscal represent mimics of similar megascale structures appearing in astrophysical plasm systems like nebulas, as will be shown for the Crab Nebula.

Breakdown of the Loop Soliton Chain: Generation of Vortex Rings and Vortex Ring Instability
When the number of pulses reaches some critical level (N = Ncrit)-w show-consistent with our previous results [58] that LMI creates conditions that favo complex dynamics including a breakdown of the loop soliton lattice into individual ope loops and vortex rings.After collapse-and-reconnection, the loop solitons become unsta ble and yield a cascade of vortex rings.Then, simultaneous excitation of spectra on vorte rings occurs due to perturbation by a random inhomogeneous strain field [58] Such a scenario takes place for N ≳ 12, because pulsation causes the local strain fiel to become larger than the vorticity ω.Then, the strain field causes the collapse of th loops followed by reconnection, thus creating a cascade of vortex rings.They are furthe exposed to the straining fields and, in the cooling phase (at the end of every laser pulse to the viscosity effect; after formation, rings practically do not move from the former po sition in the string, as seen in Figure 8.The radius of the loop solitons is R ~25-30 µm, i.e., they are mostly of the same size, meaning that they travel almost at the same velocity.However, a slight difference in the loop size and curvature causes the solitons of smaller radius to travel faster than the larger ones.As a result, they approach their neighbor solitons and-an initially periodic 1D loop soliton lattice-becomes quasiperiodic (Figure 7).
Loop solitons and multiple loop solitons generated by laser plasma at the microscale represent mimics of similar megascale structures appearing in astrophysical plasma systems like nebulas, as will be shown for the Crab Nebula.

Breakdown of the Loop Soliton Chain: Generation of Vortex Rings and Vortex Ring Instability
When the number of pulses reaches some critical level (N = Ncrit)-we showconsistent with our previous results [58] that LMI creates conditions that favor complex dynamics including a breakdown of the loop soliton lattice into individual open loops and vortex rings.After collapse-and-reconnection, the loop solitons become unstable and yield a cascade of vortex rings.Then, simultaneous excitation of spectra on vortex rings occurs due to perturbation by a random inhomogeneous strain field [58].
Such a scenario takes place for N ≳ 12, because pulsation causes the local strain field to become larger than the vorticity ω.Then, the strain field causes the collapse of the loops followed by reconnection, thus creating a cascade of vortex rings.They are further exposed to the straining fields and, in the cooling phase (at the end of every laser pulse), to the viscosity effect; after formation, rings practically do not move from the former position in the string, as seen in Figure 8.
loops and vortex rings.After collapse-and-reconnection, the loop solitons become unsta-ble and yield a cascade of vortex rings.Then, simultaneous excitation of spectra on vortex rings occurs due to perturbation by a random inhomogeneous strain field [58] Such a scenario takes place for N ≳ 12, because pulsation causes the local strain field to become larger than the vorticity ω.Then, the strain field causes the collapse of the loops followed by reconnection, thus creating a cascade of vortex rings.They are further exposed to the straining fields and, in the cooling phase (at the end of every laser pulse) to the viscosity effect; after formation, rings practically do not move from the former position in the string, as seen in Figure 8. Regarding the source of the strain field, it may be said that "The source of the strain field may be classified into three categories: (a) internal or self-induced strain, (b) mutually induced strain field by companion rings, (c) external strain, varying from place to place, imposed by the pressure field of a laser pulse" [58] Subsequently, random local multipolar strains excite the unstable waves on vortex rings, which cause instability with hierarchical order.On some vortex rings, new types of instabilities appear, which have not been observed previously or have been predicted only theoretically [58].
Considering the effect of the strain field on the vortex ring instability, one can assume that in the absence of external shear, asymptotic expansions of the Navier-Stokes equations for a vortex ring may be carried out for a small parameter, ε = a/R, where a and R are the core and the ring radius, respectively.Viewed locally, the leading-order flow represents a columnar vortex with a Gaussian vorticity distribution.The velocity field V comprising radial, r, and azimuthal, Θ, velocity components can be expanded as [58] follows: where V 0 is the velocity component independent of Θ.The field V 1 component comprises terms proportional to cos Θ and sin Θ, while the quadrupole field V 2 comprises terms proportional to cos 2Θ and sin 2Θ [58] and causes deformation of the ring core into an elliptical shape; the associated instability, called elliptical instability, was discussed by Lugomer and Fukumoto [58], and references are cited there.In that case, the local strain fields-which have different symmetry properties-are arranged in increasing (higher) order, i.e., a dipole field, then a quadrupole field in second order and so on.Assuming that the "vortex ring is separated by a few diameters from the mother filament, and the ring radius R is smaller than the distance to the filament, than the selfinduced strain on a vortex ring outweighs the strain caused by neighbor filaments" [58].In that case, the frozen pictures of a distorted vortex ring can be characterized by the instability of an isolated vortex ring.A simplified stability analysis of a vortex ring, chosen as V 0 the Rankine vortex (a circular cylindrical vortex of uniform core), revealed that the instability is a parametric resonance between a pair of Kelvin waves.The analysis also revealed that the quadrupole field V 2 supports the resonance between two helical waves.In general, the resonance instability is possible at frequency and wavenumber k for which Kelvin waves with azimuthal wavenumbers m and m + 2 are simultaneously excited by disturbance υ 0 of the basic flow [58,68,69]: where the arclength parameter s is taken along the center circle of the torus [70].The major perturbation V 1 supports the resonance between pairs of Kelvin waves whose azimuthal wavenumbers differ by 1 [69,71].Curvature instability is a parametric instability of vortex rings as found by [72].Since the dipole field V 1 is tied with the vortex line curvature, this resonance is called curvature instability [68].
As shown in these laser experiments, even an isolated vortex ring can afford to accommodate a number of instability modes.Fukumoto made classification of the observed rings relies on the theories of parametric resonance and takes into account that longwavelength oscillations on isolated vortex rings-in the laser experiment-are excited by the additional external straining field [58].Namely, the multipulse LMI causes repetitive or time-periodic modulation of the multipolar strain field in the interaction space."This modulation is different in different regions due to in-homogeneity in the distribution of local strain-field symmetry.As a consequence, various types of instability are created in different regions of the quasi-linear array of rings" [58].
Vortex rings generated by the breakdown of loops under the action of cut-and-connect operation may be assumed-circular cylindrical vortices with a Gaussian vorticity (known as the Oseen vortices).Under local random multipolar strains in the laser spot, these vortices show the evolution of various types of instabilities.Comparison with theoretical results enabled hierarchical classification of these instabilities.It was shown that the velocity field of an unstable wave on Kelvin's vortex ring (that is, constant vorticity at leading order) can be expressed in terms of the Bessel and modified Bessel functions [69][70][71].The order of Bessel functions introduces azimuthal length scale, axial wavenumber, and axial length scale, as well as radial wavenumber for nodal structures of the core.They originate from parametric resonance instability caused by the effect of toroidal curvature [68][69][70][71][72][73] and are called Widnall-Bliss-Tsai (WBT) instability [74].
An important finding of this study-consistent with our previous results [58]-is that among many structures, four types of vortex ring instability modes can be identified in laser interactions as the basic structures; two of them (I)-left-and right-handed helical waves and (II)-curvature instability are parametric resonances, while III and IV are the long wavelength oscillations.
Type I of the vortex ring instability obtained by LMI s is based on "WBT instability", which originates from the elliptical core deformation and includes parametric resonance of the left-and right-handed helical waves (m = ±1).This type includes also helicalhelical wave resonance (m, m + 2) = (−1, 1), which consists of the left-and right-handed quasi-stationary helical waves.
This instability on the vortex rings of the diameter 2R~28-30 µm (comparable with the radius 2R~26 µm in ref. [58]) is given in Figure 9. Vortex rings generated by the breakdown of loops under the action of cut-and-connect operation may be assumed-circular cylindrical vortices with a Gaussian vorticity (known as the Oseen vortices).Under local random multipolar strains in the laser spot, these vortices show the evolution of various types of instabilities.Comparison with theoretical results enabled hierarchical classification of these instabilities.It was shown that the velocity field of an unstable wave on Kelvin's vortex ring (that is, constant vorticity at leading order) can be expressed in terms of the Bessel and modified Bessel functions [69][70][71].The order of Bessel functions introduces azimuthal length scale, axial wavenumber, and axial length scale, as well as radial wavenumber for nodal structures of the core.They originate from parametric resonance instability caused by the effect of toroidal curvature [68][69][70][71][72][73] and are called Widnall-Bliss-Tsai (WBT) instability [74].
An important finding of this study-consistent with our previous results [58]-is that among many structures, four types of vortex ring instability modes can be identified in laser interactions as the basic structures; two of them (I)-left-and right-handed helical waves and (II)-curvature instability are parametric resonances, while III and IV are the long wavelength oscillations.
Type I of the vortex ring instability obtained by LMI s is based on "WBT instability", which originates from the elliptical core deformation and includes parametric resonance of the left-and right-handed helical waves (m = ±1).This type includes also helical-helical wave resonance (m, m + 2) = (−1, 1), which consists of the left-and right-handed quasi-stationary helical waves.
This instability on the vortex rings of the diameter 2R~28-30 μm (comparable with the radius 2R~26 μm in ref. [58]) is given in Figure 9. Type II of the vortex ring instability is also based on WBT instability and identified as the "curvature instability", which includes parametric resonance between the helical (m = 1) and the bulge (m = 0) modes [58].Vortex ring instability of the helical-bulge wave resonance (m, m + 1) = (0,1) is shown in Figure 10.Type II of the vortex ring instability is also based on WBT instability and identified as the "curvature instability", which includes parametric resonance between the helical (m = 1) and the bulge (m = 0) modes [58].Vortex ring instability of the helical-bulge wave resonance (m, m + 1) = (0,1) is shown in Figure 10.Type III of vortex ring instability is identified as the Kelvin axisymmetric mode or the "bulge wave" (m = 0) [58].The "bulge wave" for which the S1-symmetry of the circular core is created by the combination of the axisymmetric (m = 0) and the bending (m = 1) waves (Figure 11). Figure 11a shows the initial stage of instability while Figure 11b shows the final stage after the breakup of the ring when the oscillation reaches the critical level.The schematic illustration in Figure 11c(i-iv) shows a standing deformation mode of the core cross-section shown in Figure 10 of Kop'ev and Chernyshev in Ref. [75].
Type IV of the vortex ring instability is identified as the "standing deformation mode", created by the long wavelength oscillations of multiple Kelvin waves combined Type III of vortex ring instability is identified as the Kelvin axisymmetric mode or the "bulge wave" (m = 0) [58].The "bulge wave" for which the S1-symmetry of the circular core is created by the combination of the axisymmetric (m = 0) and the bending (m = 1) waves (Figure 11).Type III of vortex ring instability is identified as the Kelvin axisymmetric mode or the "bulge wave" (m = 0) [58].The "bulge wave" for which the S1-symmetry of the circular core is created by the combination of the axisymmetric (m = 0) and the bending (m = 1) waves (Figure 11). Figure 11a shows the initial stage of instability while Figure 11b shows the final stage after the breakup of the ring when the oscillation reaches the critical level.The schematic illustration in Figure 11c(i-iv) shows a standing deformation mode of the core cross-section shown in Figure 10 of Kop'ev and Chernyshev in Ref. [75].
Type IV of the vortex ring instability is identified as the "standing deformation mode", created by the long wavelength oscillations of multiple Kelvin waves combined Figure 11a shows the initial stage of instability while Figure 11b shows the final stage after the breakup of the ring when the oscillation reaches the critical level.The schematic illustration in Figure 11c(i-iv) shows a standing deformation mode of the core cross-section shown in Figure 10 of Kop'ev and Chernyshev in Ref. [75].
Type IV of the vortex ring instability is identified as the "standing deformation mode", created by the long wavelength oscillations of multiple Kelvin waves combined of m = 0 and m = ±1 modes, shown in Figure 12a and schematic illustration, and Figure 12b of Kop'ev and Chernyshev [75].
of m = 0 and m = ±1 modes, shown in Figure 12a and schematic illustration, and Figure 12b of Kop'ev and Chernyshev [75].The findings of four types of vortex ring instabilities indicate that they are caused by nonselective, spatially inhomogeneous shock-induced excitation-in the regions of the low filament density (Λ > σ).Under conditions dominating in these regions, a vortex ring can accommodate a number of instability modes creating the structures related to a parametric resonance instability caused by dipolar and quadrupolar fields.This gives a clearly better understanding of the shock wave effects, which create vorticity, but also may alter the picture established by a theory of homogeneous media.

Vortex Filament "Complexes"
Besides the above findings, characteristic of the regions of a low filament density, our experiments also demonstrate that multipulse laser interaction (with a Co-coated steel target) creates very complex patterns-in the regions of a high vortex filament density (Λ~σ).Vortex filament dynamics leading to such structures are caused by inhomogeneous turbulence and various strain fields in the fluid layer.Under the action of strain fields, vortex filaments form coils, helices, bundles, and other structures, which we refer to as "complexes" [61] The strain fields can be (quasi-)static or oscillatory, creating torsional and twisting filamentary structures ranging from the nanometer to the millimeter scale.With increasing N, filaments become collected into axially merged bundles, merged into thick filaments, or transformed into ribbons and tubular ribbons.Their interaction with localized strain fields gives rise to various types of complex structures with complexity being controlled by the number of pulses and the density of the vortex filament array [61].

Formation of Vortex Filament "Complexes"
A series of laser pulses causes the motion of vortex filaments and collective behavior with the formation of vortex filament bundles, helically paired structures, and the formation of multiple loop solitons on bundles.While the first pulses cause the roll-up of the shear layer into vortex filaments (for Re > Recritical; Recritical~10 4 ), and the formation of a parallel array, the other pulses cause their motion.Strong transversal oscillation and motion of filaments cause their aggregation in the regions of the minimum strain [61].If lateral (parallel to the surface) and vertical (normal) strain field components are comparable, the filaments move closer and form a bundle.Immersed in the turbulent background fluid, vortex filament bundles also experience the oscillatory strain field of background fluid, which causes twisting and winding, helical pairing, formation of double helices of bundles, braiding, and tangling, and formation of loops on the bundles as well, giving rise to the above-mentioned complex structures, which we refer to as "complexes" [61].The findings of four types of vortex ring instabilities indicate that they are caused by nonselective, spatially inhomogeneous shock-induced excitation-in the regions of the low filament density (Λ > σ).Under conditions dominating in these regions, a vortex ring can accommodate a number of instability modes creating the structures related to a parametric resonance instability caused by dipolar and quadrupolar fields.This gives a clearly better understanding of the shock wave effects, which create vorticity, but also may alter the picture established by a theory of homogeneous media.

Vortex Filament "Complexes"
Besides the above findings, characteristic of the regions of a low filament density, our experiments also demonstrate that multipulse laser interaction (with a Co-coated steel target) creates very complex patterns-in the regions of a high vortex filament density (Λ ~σ).Vortex filament dynamics leading to such structures are caused by inhomogeneous turbulence and various strain fields in the fluid layer.Under the action of strain fields, vortex filaments form coils, helices, bundles, and other structures, which we refer to as "complexes" [61] The strain fields can be (quasi-)static or oscillatory, creating torsional and twisting filamentary structures ranging from the nanometer to the millimeter scale.With increasing N, filaments become collected into axially merged bundles, merged into thick filaments, or transformed into ribbons and tubular ribbons.Their interaction with localized strain fields gives rise to various types of complex structures with complexity being controlled by the number of pulses and the density of the vortex filament array [61].

Formation of Vortex Filament "Complexes"
A series of laser pulses causes the motion of vortex filaments and collective behavior with the formation of vortex filament bundles, helically paired structures, and the formation of multiple loop solitons on bundles.While the first pulses cause the roll-up of the shear layer into vortex filaments (for Re > Re critical ; Re critical ~10 4 ), and the formation of a parallel array, the other pulses cause their motion.Strong transversal oscillation and motion of filaments cause their aggregation in the regions of the minimum strain [61].If lateral (parallel to the surface) and vertical (normal) strain field components are comparable, the filaments move closer and form a bundle.Immersed in the turbulent background fluid, vortex filament bundles also experience the oscillatory strain field of background fluid, which causes twisting and winding, helical pairing, formation of double helices of bundles, braiding, and tangling, and formation of loops on the bundles as well, giving rise to the above-mentioned complex structures, which we refer to as "complexes" [61].
Such filament-bundle "complex", which forms a double helix shown in Figure 13a, (and in the refined micrograph in Figure 13b), is obtained for the first time.
Such filament-bundle "complex", which forms a double helix shown in Figure 13a, (and in the refined micrograph in Figure 13b), is obtained for the first time.Schematic reconstruction of vortex filaments organized into bundles, which form a double helix in Figure 13, is shown in three consecutive steps, starting with two bundles of parallel vortex filaments, which are exposed to torsion, is shown in Figure 14.Schematic reconstruction of the vortex filament "complex" in Figure 14.Formation of two bundles of vortex filaments (i).Turbulent background field causes dense packing of filaments; twisting of bundles causes helical winding of individual bundles but also their helical pairing (ii).The peripheral filament in the bundle experiences additional torsion, which starts to form a series of loop solitons (iii).Closely packed filaments under the action of the surrounding turbulent background field experience core spreading and axial merging.Schematic reconstruction of vortex filaments organized into bundles, which form a double helix in Figure 13, is shown in three consecutive steps, starting with two bundles of parallel vortex filaments, which are exposed to torsion, is shown in Figure 14.

Figure 14.
Schematic reconstruction of the vortex filament "complex" in Figure 14.Formation of two bundles of vortex filaments (i).Turbulent background field causes dense packing of filaments; twisting of bundles causes helical winding of individual bundles but also their helical pairing (ii).The peripheral filament in the bundle experiences additional torsion, which starts to form a series of loop solitons (iii).Closely packed filaments under the action of the surrounding turbulent background field experience core spreading and axial merging.
While many filaments are firmly incorporated into bundles, some of them (slightly separated from the mother bundle) show individual behavior with the formation of the loop solitons-thus giving rise to the filament-bundle "complex".
Extensive results carried out also show the formation of even more complex bundle structures in multipulse LMI. Figure 15 shows that vortex filament bundles exposed to torsion and twist-perturbation cause the successive formation of (+) and (−) soliton loops, as well as the filament coiling-which we obtained for the first time.
While many filaments are firmly incorporated into bundles, some of them (slightly separated from the mother bundle) show individual behavior with the formation of the loop solitons-thus giving rise to the filament-bundle "complex".
Extensive results carried out also show the formation of even more complex bundle structures in multipulse LMI. Figure 15 shows that vortex filament bundles exposed to torsion and twist-perturbation cause the successive formation of (+) and (−) soliton loops, as well as the filament coiling-which we obtained for the first time.This picture (left side) also shows two connected rings, which are formed after loop soliton breakup, meaning that they are formed by the "cut and connect" operation.Once formed, the rings stay in the vicinity of the mother filament, but the strong core diffusion makes the surrounding ambient in a turbulent background field, rather misty.An approximate reconstruction of this "complex" is given in Figure 16a.Schematic illustration of the interaction of (+) and (−) loop solitons (i) leads to their annihilation (ii) as shown in Figure 16b.With subsequent laser pulses, a bundle of vortex filaments is exposed to growing compression, which invites resistance to the filament motion and causes the axial to merge into a thick vortex filament [61].In that case, the formation of a vortex filament bundle consisting of thick and thin filaments (under an oscillating pressure field) is associated with the motion of the peripheral filaments, which This picture (left side) also shows two connected rings, which are formed after loop soliton breakup, meaning that they are formed by the "cut and connect" operation.Once formed, the rings stay in the vicinity of the mother filament, but the strong core diffusion makes the surrounding ambient in a turbulent background field, rather misty.An approximate reconstruction of this "complex" is given in Figure 16a.
Schematic illustration of the interaction of (+) and (−) loop solitons (i) leads to their annihilation (ii) as shown in Figure 16b.With subsequent laser pulses, a bundle of vortex filaments is exposed to growing compression, which invites resistance to the filament motion and causes the axial to merge into a thick vortex filament [61].In that case, the formation of a vortex filament bundle consisting of thick and thin filaments (under an oscillating pressure field) is associated with the motion of the peripheral filaments, which become slightly shifted from thick filaments.Under torsion, the helical motion of thick and thin filaments creates "complex" with (+) and (−) loop solitons (Figure 17).become slightly shifted from thick filaments.Under torsion, the helical motion of thick and thin filaments creates "complex" with (+) and (−) loop solitons (Figure 17).Under a series of pulses, thick and thin filaments may show individual translative motion and wiggling; however, a thin filament mostly follows the curvature of a thick one.In the inhomogeneous flow field of the laser spot, such vortex filaments experience different torsion, stretching, etc.Therefore, the loop solitons of different radii and curvature will be formed making a highly complex structure.In the oscillatory strain field of laser shocks, vortex filaments may also form a series of kinks and, if they come into close vicinity, may axially merge or form chaotic structures depending on the number N of laser pulses.When N reaches the critical level, N = Ncrit, the loop-soliton breakup will form open loops, which may become connected to vortex rings [61].
An important finding in the understanding of this phenomenon is that the process starts with one filament approaching its neighbor; the approaching filament longitudinally stretches while its cross-section narrows.As the angular momentum is conserved, the vortex filament spins faster enhancing the co-rotation of the same signed filamentary vortices, intensifying viscous diffusion.The crucial role in axial merging-and the formation of thick vortex filament-is played by the viscous diffusion term, a source term due to stretching, but also an advection term with variable advection speed [61].Intense laser pulses of ~10 8 W/cm 2 , create plasma-which together with molten target surface-establishes the layer of low-and high-density fluids, separated by the density interface.The shock-induced acceleration of density interface and variable density flow lead to the RM/RT instability and the formation of vortex-filament structures like bundles and their "complexes".Since very little is known about the filament-bundle "complexes"-the above results provide evidence on their formation and behavior in the oscillatory shock-induced environment.
The mechanisms responsible for their formation and behavior operate in and between many different scales.Along this line of thinking, the nonlinear and nonequilibrium processes that lead to the formation of the small-scale structures created by laser experiments on metal targets at medium power density are mimic of those characteristic for the large-scale near-peripheral region of the Crab Nebula.

Basic Aspects of Astrophysical Plasma Systems: Nebulas
Astrophysical plasma systems like Cosmic nebulas are generated by various mechanisms and appear in various shapes and sizes, which extend over many light years.Some nebulas are formed from interstellar gas (dust particles and molecular species), which under twisting magnetic lines creates helicoidal filamentary structures-such as Double-Helix Nebula.The other nebulas created by the collapse of stars represent the Under a series of pulses, thick and thin filaments may show individual translative motion and wiggling; however, a thin filament mostly follows the curvature of a thick one.In the inhomogeneous flow field of the laser spot, such vortex filaments experience different torsion, stretching, etc.Therefore, the loop solitons of different radii and curvature will be formed making a highly complex structure.In the oscillatory strain field of laser shocks, vortex filaments may also form a series of kinks and, if they come into close vicinity, may axially merge or form chaotic structures depending on the number N of laser pulses.When N reaches the critical level, N = N crit , the loop-soliton breakup will form open loops, which may become connected to vortex rings [61].
An important finding in the understanding of this phenomenon is that the process starts with one filament approaching its neighbor; the approaching filament longitudinally stretches while its cross-section narrows.As the angular momentum is conserved, the vortex filament spins faster enhancing the co-rotation of the same signed filamentary vortices, intensifying viscous diffusion.The crucial role in axial merging-and the formation of thick vortex filament-is played by the viscous diffusion term, a source term due to stretching, but also an advection term with variable advection speed [61].Intense laser pulses of ~10 8 W/cm 2 , create plasma-which together with molten target surface-establishes the layer of low-and high-density fluids, separated by the density interface.The shock-induced acceleration of density interface and variable density flow lead to the RM/RT instability and the formation of vortex-filament structures like bundles and their "complexes".Since very little is known about the filament-bundle "complexes"-the above results provide evidence on their formation and behavior in the oscillatory shock-induced environment.
The mechanisms responsible for their formation and behavior operate in and between many different scales.Along this line of thinking, the nonlinear and nonequilibrium processes that lead to the formation of the small-scale structures created by laser experiments on metal targets at medium power density are mimic of those characteristic for the large-scale near-peripheral region of the Crab Nebula.

Basic Aspects of Astrophysical Plasma Systems: Nebulas
Astrophysical plasma systems like Cosmic nebulas are generated by various mechanisms and appear in various shapes and sizes, which extend over many light years.Some nebulas are formed from interstellar gas (dust particles and molecular species), which under twisting magnetic lines creates helicoidal filamentary structures-such as Double-Helix Nebula.The other nebulas created by the collapse of stars represent the remnants of supernovae-actually, a thin plasmatic shock area at its boundary exposed to magnetic fields, which affect the behavior of expanding matter [1].A spherical blast wave propagates from the center of a star outwards through the star's layers of gradually less-dense gases, causing the shock-induced variable density flow thus generating the RT instability [1,76] and wave-vortex structures.This shock wave expands into the surrounding interstellar medium, sweeping up an expanding shell of gas and dust.The ejected debris mixes with the interstellar medium and forms a nebulas a supernova remnant [76,77], such as the Crab Nebula.
The picture of the Crab Nebula (Figure 18) reveals its composition and structure [78] and provides a starting point for the study of regular and irregular structures and cluster formations.The diffuse blue region of Crab Nebula is generated by the synchrotron radiation caused by electrons moving at relativistic speeds [79] along trajectories curved by the strong magnetic field, which originates from a spinning neutron star-a pulsar (the ultra-dense core of the exploded star)-at the center of the Nebula [80].Electrons circulating around the magnetic field lines create this blue light in the interior of the synchrotron nebula, which is bounded by the thermal ejecta.Both the synchrotron nebula and filament structures are embedded in the expanding supernova plasmatic remnant.According to Hester [81], these structures-called thermal filaments-are composed of ejecta from the explosion and located around the outer part of the synchrotron nebula.The ejecta evolving along the magnetic field lines take the form of vortex filaments, filament jets, knotted filaments, loops, etc., which are dominated by the line emission [79,82] They consist of ionized helium and hydrogen, as well as carbon, oxygen, nitrogen, iron, neon, and sulfur, which create colored structures [82] (Figure 18).Spectroscopic diagnostics reveal that orange filaments consist mostly of hydrogen, while green ones consist of singly ionized sulfur.In the outer part of the nebula, dominating blue filaments originate from neutral oxygen, while red ones originate from doubly ionized oxygen [78].This highly ionized red region makes a "skin" around the synchrotron emission.This red "skin" is driven into a larger, expanding remnant, which surrounds the crab.Within this expanding remnant, there are dense structures of "thermal filaments" with self-shielded cores that are highly ionized [59].Spectroscopy of multiply ionized Oxygen III reveals that temperatures of thermal filaments range from ~12,000 K to 20,000 K [59].Somewhat different values are obtained from other measurements, which give the filaments' temperatures about ~11,000-18,000 K, and their densities about ~1300 particles per cm 3 [63] Morphological properties of these particle structures in the Crab Nebula are described in [83].
Particles and the wind of waves streaming in a radial direction are emitted from the pulsar in the central part of the nebula, which loses its rotational energy emitting a relativistic wind of waves and particles.They interact with the surrounding ambient, generating pulsar wind nebulae (PWNe), observable from radio to γ-rays.PWNe often shows a torus-jet structure [84].Orientation of the torus-jet structures is referenced with respect to the line of sight to the Earth: For the view of the front side of the remnant, the north is up and the east is to the left.The other view 60 • E from the line of sight to the Earth, and 30 • S of the E-W plane gives a different picture.This is roughly equivalent to viewing the remnant from the southeast, sighting along the major optical axis.All jets and filament structures are arranged in a pattern that can best be described in relation to the plane of the E-W torus [85].In this respect, the most impressive magnetized jet is the "South-East (SE) jet" of the Crab Nebula, which was recently studied in scaled experiments by ultraintense laser [5].
Regarding the interaction of ejected plasma particles with the surrounding medium, the study of Kuranzet al. [86] reveals that various unstable processes in the nebula create convoluted structures at boundaries between regions, with significant interpenetration of the plasmatic matter from one region to the next.This is evident in the regions where spikes of high-density stellar ejecta penetrate into the PWN [86], initiating the Rayleigh-Taylor instability and mixing.Rayleigh-Taylor instability leads to intense mixing, the characteristics of which depend on the local velocity, scale of structure growth, and the scale of energy dissipation, which determine the width of the mixing zone-mostly different in different regions of the huge plasmatic system.In addition, there are suggestions that some other processes may also play a role in the material mixing scenario-besides RT and RM instabilities [87].
The velocity of the plasma jet at the Crab Nebula periphery is ~100 km/s with respect to the nebular expansion velocity [88].Based on the radio image, the equilibrium magnetic field of the "jet" appears well ordered, while the density of the thermal gas inside the plasma jet is estimated to be ~0.7/cm 3 , and the Alfven velocity is more than ~700 km/s.It is suggested that "The "jet" is probably confined most likely by a toroidal magnetic field which must be of the same order of magnitude, i.e., 3 × 10 −4 G, resulting in a helical field" [88] For comparison, in solar flares, the magnetic fields are of the order of hundreds of Gauss with a maximum of ~10-50 G [89].Also, the plasma density in solar flares is of the order of 10 9 to 10 10 cm −3 [89], while in the central part of the nebula (pulsar wind nebulas), the ambient number density is only ~0.1 cm −3 [90,91].
The surface magnetic dipole of radio pulsars has a strength of B~10 11−13 G [92],which for the Crab Nebula is estimated to be ~7.6 × 10 12 G [93] , while in its peripheral region the strength is ~300-500 µG [94].The strength of the magnetic field is the parameter, which-in combination with velocity and plasma particle density-dictates the dynamics, morphological, and topological characteristics of ejecta.The ejecta take the form of long magnetized jets in the near central region of the strong magnetic field, and the form of filaments, filament loops, etc., in the near peripheral region of the weak magnetic field, in the Crab Nebula.These plasmatic structures in the nebula can be analyzed relaying on the analogy with similar structures evolving from nonlinear and nonequilibrium processes at the microscale created by intense and ultraintense lasers.

Crab Nebula: Filament "Complexes" in the near Periphery Region
Crab Nebula (Figure 18) consists of (i) an external supernova remnant (SNR) as a diffuse, expanding nebula plasma shell, and (ii) an inner pulsar wind nebula (PWN) inside the shell of a supernova remnant (SNR), driven by winds from the central pulsar.

Magnetized Vortex Filaments and Loop Solitons in the Peripheral Region of the Crab Nebula: Weak Magnetic Field
A diffuse supernova remnant of expanding nebula shell at the external side of the PWN/SNR boundary, where the magnetic field is weak, represents the peripheral and the near peripheral region of the Crab Nebula.A large magnification (besides filaments) also shows filament loop solitons, vortex rings, and broken rings (Figure 19).
These loop solitons (created by winding of particle trajectories along the magnetic field lines) are similar to those in LMI (created by winding of vortex filaments by mechanical torsional strain).A number of random loop solitons in the right peripheral crab region can also be seen in Figure 19.These loop solitons appear in different colors in two different parts of this region, characteristic of different plasma compositions of ionized gases.They are mostly closed loops, indicating equal (constant) plasma acceleration (equal magnetic field strength and winding), meaning that the instability of these filaments did not reach the critical level for their breakup.Crab Nebula (Figure 18) consists of (i) an external supernova remnant (SNR) as a diffuse, expanding nebula plasma shell, and (ii) an inner pulsar wind nebula (PWN) inside the shell of a supernova remnant (SNR), driven by winds from the central pulsar.

Magnetized Vortex Filaments and Loop Solitons in the Peripheral Region of the Crab Nebula: Weak Magnetic Field
A diffuse supernova remnant of expanding nebula shell at the external side of the PWN/SNR boundary, where the magnetic field is weak, represents the peripheral and the near peripheral region of the Crab Nebula.A large magnification (besides filaments) also shows filament loop solitons, vortex rings, and broken rings (Figure 19).These loop solitons (created by winding of particle trajectories along the magnetic field lines) are similar to those in LMI (created by winding of vortex filaments by mechanical torsional strain).A number of random loop solitons in the right peripheral crab region can also be seen in Figure 19.These loop solitons appear in different colors in two different parts of this region, characteristic of different plasma compositions of ionized gases.They are mostly closed loops, indicating equal (constant) plasma acceleration (equal magnetic field strength and winding), meaning that the instability of these filaments did not reach the critical level for their breakup.
Very similar loop soliton structures can be seen in other peripheral regions of the Crab Nebula.An example is the upper peripheral left region.The loop solitons also appear in two different colors, but in contrast to the right region, the loops are mostly broken (in the green part), while some of the loops are closed or even appear as rings (in the brown part) of this region.This indicates spatially variable plasma acceleration (or variable magnetic field strength), meaning that the instability of filaments reached the critical Very similar loop soliton structures can be seen in other peripheral regions of the Crab Nebula.An example is the upper peripheral left region.The loop solitons also appear in two different colors, but in contrast to the right region, the loops are mostly broken (in the green part), while some of the loops are closed or even appear as rings (in the brown part) of this region.This indicates spatially variable plasma acceleration (or variable magnetic field strength), meaning that the instability of filaments reached the critical level for their breakup in the green part, and did not reach it in the brown part of the region.This may indicate the inhomogeneity of the plasmatic and magnetic field(Figure 20).Conjecture may be established with the multipulse LMI, where the loop solitons on vortex filaments reach the instability threshold after a critical number of pulses and spatially variable acceleration-becoming broken into open loops and later on into vortex rings.Common to both cases is the underlying helical motion of the vortex filament (in LMI) and of the magnetized filament (in nebula), which is characteristic of many astrophysical systems.In this respect, the formation of the loop solitons on magnetized vortex filaments in the Crab Nebula-as the nonlinear and nonequilibrium phenomenon-is analogous to the formation of the knot solitons [95] to the visible hot spot soliton-like points on magnetized jets and to the solitons caused by torsion of the magnetic field lines (filament jets) in "magnetars"-the magnetic stars [96][97][98].Considering magnetized knotted jets in active galactic nuclei (radio galaxy 3C303), Lapenta and Kronberg [95] defined knots as "visual manifestation of an underlying helical form of a jet, which can also be explained as shocks along the jets that heat the plasma particles, creating visible hot spot".The other interpretations are based on "the Kelvin-Helmholtz instability, and on current-driven modes, but also on the model which considers the jets with knotted solitary structure as actual plasma features, which are ejected from the accretion disc and propagate along the jet axis" [95].Based on this approach, Lapenta [95] used a mathematical analogy between the stationary magnetohydrodynamic (MHD) equations (Grad-Shafranov equation) and the nonlinear Schrodinger equation for soliton propagation.Taking the magnetized jet that propagates along the z-direction perpendicular to the (x,y) plane, and neglecting the y-direction, the problem is reduced to 2D solutions in the (x,z) plane, with the system size Lx and Lz [95].The solution is constructed as a mathematical analogy to the soliton equation for nonlinear wave propagation (cubic Schrodinger equation).The magnetic field is obtained from a flux function as described by Lapenta and Kromberg [95] ikz p e L x h Ψ ) / ( sec = Ψ (20) where ᴪp is the on-plane magnetic field.This equation for the magnetized jet with knotty structure is just the same equation of Hasimoto transformation (11) for the loop soliton on Conjecture may be established with the multipulse LMI, where the loop solitons on vortex filaments reach the instability threshold after a critical number of pulses and spatially variable acceleration-becoming broken into open loops and later on into vortex rings.Common to both cases is the underlying helical motion of the vortex filament (in LMI) and of the magnetized filament (in nebula), which is characteristic of many astrophysical systems.In this respect, the formation of the loop solitons on magnetized vortex filaments in the Crab Nebula-as the nonlinear and nonequilibrium phenomenon-is analogous to the formation of the knot solitons [95] to the visible hot spot soliton-like points on magnetized jets and to the solitons caused by torsion of the magnetic field lines (filament jets) in "magnetars"-the magnetic stars [96][97][98].Considering magnetized knotted jets in active galactic nuclei (radio galaxy 3C303), Lapenta and Kronberg [95] defined knots as "visual manifestation of an underlying helical form of a jet, which can also be explained as shocks along the jets that heat the plasma particles, creating visible hot spot".The other interpretations are based on "the Kelvin-Helmholtz instability, and on current-driven modes, but also on the model which considers the jets with knotted solitary structure as actual plasma features, which are ejected from the accretion disc and propagate along the jet axis" [95].Based on this approach, Lapenta [95] used a mathematical analogy between the stationary magnetohydrodynamic (MHD) equations (Grad-Shafranov equation) and the nonlinear Schrodinger equation for soliton propagation.Taking the magnetized jet that propagates along the z-direction perpendicular to the (x,y) plane, and neglecting the y-direction, the problem is reduced to 2D solutions in the (x,z) plane, with the system size Lx and Lz [95].The solution is constructed as a mathematical analogy to the soliton equation for nonlinear wave propagation (cubic Schrodinger equation).The magnetic field is obtained from a flux function as described by Lapenta and Kromberg [95].
where Ψ p is the on-plane magnetic field.This equation for the magnetized jet with knotty structure is just the same equation of Hasimoto transformation (11) for the loop soliton on vortex filament in LMI (where Ψ p corresponds to torsion) [95], which makes an analogy between similar instabilities at the micro-and megascale.
The system evolution has two phases: a long period of stability followed by sudden instability.The stable phase corresponds to the long collimation phase of the magnetized filament jet, appearing as a straight regular structure.The knots characterize the periodic "bubbly" structure of the solitons.Such soliton-like solutions reveal that a magnetic field amplitude grows through periodic alternating minima and maxima.This feature is common to the well-known coherence of solitons in many areas of physics, in which nonlinear equations admit soliton solutions [95].
The origin of loop-soliton-like structures was also studied as the result of the torsion of magnetic vortex line curves [96].According to Ricca [97], torsion in twisted vortex filaments affects significantly the motion of helical vortex filaments, in the fluid, where the binormal component is responsible for the displacement of the vortex filament.Based on this approach, Andarde [96] solved the equilibrium equations of magnetic stars for vortex filaments and wrote the magnetar equations in the Frenet-Serret frame for constant torsion.It was shown that the magnetic field oscillates in a helical form, and the behavior of magnetized vortex filament follows the Hasimoto transformation of the Schrodinger equation for the constant torsion creating loop solitary structures.Long vortex filaments and underlying helical motion due to magnetic torsional strain field(s) lead to the formation of the series of loop solitons, which are quite similar to the string of the loop solitons on vortex filaments in LMI.The creation of multisolitons of different radii and curvature, which travel at different velocities, and their nonlinear interaction may cause their complex organization in space.A complex case occurs when two solitons of the same polarity of the wavevector (either + or −) collide from behind, and a more complex case occurs when they have different polarity [66].
Such a loop-soliton chain on magnetized filamentary jets in the peripheral crab region is shown in Figure 21.vortex filament in LMI (where ᴪp corresponds to torsion) [95], which makes an analogy between similar instabilities at the micro-and megascale.
The system evolution has two phases: a long period of stability followed by sudden instability.The stable phase corresponds to the long collimation phase of the magnetized filament jet, appearing as a straight regular structure.The knots characterize the periodic "bubbly" structure of the solitons.Such soliton-like solutions reveal that a magnetic field amplitude grows through periodic alternating minima and maxima.This feature is common to the well-known coherence of solitons in many areas of physics, in which nonlinear equations admit soliton solutions [95] The origin of loop-soliton-like structures was also studied as the result of the torsion of magnetic vortex line curves [96].According to Ricca [97], torsion in twisted vortex filaments affects significantly the motion of helical vortex filaments, in the fluid, where the binormal component is responsible for the displacement of the vortex filament.Based on this approach, Andarde [96] solved the equilibrium equations of magnetic stars for vortex filaments and wrote the magnetar equations in the Frenet-Serret frame for constant torsion.It was shown that the magnetic field oscillates in a helical form, and the behavior of magnetized vortex filament follows the Hasimoto transformation of the Schrodinger equation for the constant torsion creating loop solitary structures.Long vortex filaments and underlying helical motion due to magnetic torsional strain field(s) lead to the formation of the series of loop solitons, which are quite similar to the string of the loop solitons on vortex filaments in LMI.The creation of multisolitons of different radii and curvature, which travel at different velocities, and their nonlinear interaction may cause their complex organization in space.A complex case occurs when two solitons of the same polarity of the wavevector (either + or −) collide from behind, and a more complex case occurs when they have different polarity [66].
Such a loop-soliton chain on magnetized filamentary jets in the peripheral crab region is shown in Figure 21.

Vortex Ring Instabilities in the near Peripheral Region of the Crab Nebula
Peripheral regions of the Crab Nebula under magnification also reveal deformed and broken rings as well as the open loops in the vicinity of the string of the loop solitons on filaments (Figure 22).Morphological and topological characteristics of broken loops, rings, and deformed vortex rings of the Crab Nebula indicate that these ring instabilities may eventually be attributed to the same classification as those generated in the multipulse LMI, shown in Figures 10-13.
Based on this classification, characteristics of the shape, symmetry, and oscillation modes of broken rings and open loops of the crab in Figure 22 indicate the ring axisymmetric instability of type III, called the bulge wave.This instability still has circular symmetry and corresponds to the case of the vortex axis tilted [58].However, broken rings indicate that this instability has reached a critical level for the integrity of the rings.

Evolution of Rayleigh-Taylor Instability at the PWN/SNR Boundary of the Crab Nebula
A diffuse supernova remnant of an expanding nebula shell at the PWN/SNR boundary experiences the pulsar wind from the central region, which expands radially driving a hot plasmatic bubble (PWN-bubble) into a colder external SNR shell.The high pressure of the hot PWN bubble drives a shock wave into the high-density cold ejecta causing the formation of a density interface [99,100].The interface perturbation evolves into Rayleigh-Taylor (RT) instability with the growth of spike fingers oriented downwards and bubbles growing between the fingers (RTI spikes).According to [101,102], the RTI is the origin of the "thermal filaments" observed in the near peripheral regions of the Crab Nebula-the outer side of the PWN/SNR boundary.
At the contact discontinuity between PWN and SNR regions of the Crab Nebula, these RTI structures dominate.Using the self-similar model of PWN inflating the ejecta with density ρ ∝r −α by a pulsar wind, it was found that the shock speed that causes growth of RTI can be expressed [103,104] as follows: where α is the index of the jet ejecta density distribution; for uniform ejecta α = 0. Starting with the linear regime, the RTI growth switches into the non-linear regime when the Morphological and topological characteristics of broken loops, rings, and deformed vortex rings of the Crab Nebula indicate that these ring instabilities may eventually be attributed to the same classification as those generated in the multipulse LMI, shown in Figures 10-13.
Based on this classification, characteristics of the shape, symmetry, and oscillation modes of broken rings and open loops of the crab in Figure 22 indicate the ring axisymmetric instability of type III, called the bulge wave.This instability still has circular symmetry and corresponds to the case of the vortex axis tilted [58].However, broken rings indicate that this instability has reached a critical level for the integrity of the rings.

Evolution of Rayleigh-Taylor Instability at the PWN/SNR Boundary of the Crab Nebula
A diffuse supernova remnant of an expanding nebula shell at the PWN/SNR boundary experiences the pulsar wind from the central region, which expands radially driving a hot plasmatic bubble (PWN-bubble) into a colder external SNR shell.The high pressure of the hot PWN bubble drives a shock wave into the high-density cold ejecta causing the formation of a density interface [99,100].The interface perturbation evolves into Rayleigh-Taylor (RT) instability with the growth of spike fingers oriented downwards and bubbles growing between the fingers (RTI spikes).According to [101,102], the RTI is the origin of the "thermal filaments" observed in the near peripheral regions of the Crab Nebula-the outer side of the PWN/SNR boundary.
At the contact discontinuity between PWN and SNR regions of the Crab Nebula, these RTI structures dominate.Using the self-similar model of PWN inflating the ejecta with density ρ ∝r −α by a pulsar wind, it was found that the shock speed that causes growth of RTI can be expressed [103,104] as follows: where α is the index of the jet ejecta density distribution; for uniform ejecta α = 0. Starting with the linear regime, the RTI growth switches into the non-linear regime when the amplitude of the interface distortion becomes comparable to the wavelength λ.At the onset of this regime, the light fluid creates bubbles/columns of diameter ~λ, which steadily rise with the speed: while the heavy fluid forms thin fingers in the Crab Nebula approaching the free-fall regime, as observed by a number of authors between 1950 and 2012 (listed by Porth et al. [85] and later shown by more complex 3D simulations. Regarding the location (regions) of their appearance, there are statements like "Most individual filaments are small-scale structures but some are much longer and appear to cross almost the entire nebula", but also "the filaments with low line-of-sight speed avoid the central region of the nebula image" [85].This indicates that the crab filaments do not penetrate the whole volume of the nebula, as otherwise low line-of-sight speed emission would be seen there.Instead, "the filaments reside near its outer edge, where they occupy a thick shell of thickness about one-third of the nebula radius" [85,104], which is of the size ~1.83 ly.In addition, some authors when describing tiny filamentary structures refer to them as "threads"-a characteristic of the near central region (the PWN dominant region)-and not for the shell region.
Numerical simulations reveal that the evolution of RTI at the contact discontinuity between PWN and SNR regions creates coherent filamentary structures.These simulated coherent RT finger-spike structures or filaments resemble some of the real crab filaments; the longest ones reach the length of ~1/4 of the nebula radius (r = 5.5 ly).Regarding the formation of such long filaments, the picture is not quite clear.Namely, as mentioned by Porth et al. [105], the thickness of the mixing layer occupied by the RTI fingers is much smaller, only approximately 1/15 of the PWN radius-which is about five times below the thickness of the Crab's filamentary shell.
However, in spite of some uncertainness, the point to be emphasized is that filaments are identified as the "Rayleigh-Taylor filaments" (actually tiny RTI spikes), resulting from the nonlinear evolution in the presence of the magnetic field.Hester et al. [94] applied the theory of magnetic RTI, assuming that the smallest structures of the Crab's filamentary network resemble the RT bubbles and fingers with the wavelength of λ = 2λ c , for ρ 2 >> ρ 1 .
Simulations have shown different effects on the RT modes with respect to the direction of the magnetic field.
For the RT modes normal to the magnetic field, the growth rate of the non-magnetic case is recovered-the magnetic field does not suppress modes normal to the magnetic field.
For RT modes parallel to the field, the magnetic tension suppresses the perturbations with wavelengths below the critical one, λ c (full suppression).
The presence of a magnetic field is expected to have a significant effect because pulsar winds inject highly magnetized plasma into the PWN bubble.Such a strong magnetic field should suppress the evolution of RTI and of the filamentary structures.However, "the 3D simulations have shown that because of the strong magnetic dissipation and field randomizationthe magnetic tension at the contact discontinuity between PWN and SNR regions is not strong enough to suppress the growth of RT filaments", so that prominent filaments appear in the Crab Nebula [103].

The Crab Nebula near Central Region: Formation of Filament "Complexes" 7.1. Bundles of Filaments Helically Paired: Very Strong Magnetic Field
Inside the shell of a supernova remnant (SNR)-in the near central region of the Crab Nebula-different phenomena dominate driven by winds from a central pulsar (PWN), where the magnetic field is very strong.A strong magnetic field overwhelms the plasma shear flow instability, suppresses the RTI evolution, and has a dominant role in the formation of magnetized filament jets [105].Magnification of this region of the Crab Nebula reveals long filament jets organized into bundles, which form double helices-a kind of "complexes" (Figures 23 and 24).
Nebula reveals long filament jets organized into bundles, which form double helices-a kind of "complexes" (Figures 23 and 24).Double helices-a kind of "complexes on bundles"-are caused by twisting of magnetized bundles of filament jets.Magnification of two segments of the double helix in Figure 24a,b reveals vortex rings and loops on those magnetized filaments, which are slightly separated from bundles and manifest individual behavior.These loop solitons, as well as vortex rings and broken rings in the near central region of the Crab Nebula, are quite similar to those observed on vortex filament bundles created by multipulse LMI.
Another type of magnetized filament "bundles" that form a double helix (selected from the upper central part of the Crab Nebula in Figure 18) is shown in Figure 25.Double helices-a kind of "complexes on bundles"-are caused by twisting of magnetized bundles of filament jets.Magnification of two segments of the double helix in Figure 24a,b reveals vortex rings and loops on those magnetized filaments, which are slightly separated from bundles and manifest individual behavior.These loop solitons, as well as vortex rings and broken rings in the near central region of the Crab Nebula, are quite similar to those observed on vortex filament bundles created by multipulse LMI.Nebula reveals long filament jets organized into bundles, which form double helices-a kind of "complexes" (Figures 23 and 24).Double helices-a kind of "complexes on bundles"-are caused by twisting of magnetized bundles of filament jets.Magnification of two segments of the double helix in Figure 24a,b reveals vortex rings and loops on those magnetized filaments, which are slightly separated from bundles and manifest individual behavior.These loop solitons, as well as vortex rings and broken rings in the near central region of the Crab Nebula, are quite similar to those observed on vortex filament bundles created by multipulse LMI.
Another type of magnetized filament "bundles" that form a double helix (selected from the upper central part of the Crab Nebula in Figure 18) is shown in Figure 25.Another type of magnetized filament "bundles" that form a double helix (selected from the upper central part of the Crab Nebula in Figure 18) is shown in Figure 25.
The double helix consists of fragmented magnetized filaments, which connects the blue central part and the brown-yellow peripheral part in the Crab Nebula.In contrast to bundles of the long filament jets in Figures 23 and 24, these "bundles" are composed of fine filaments (like threads) much shorter than the bundle itself; they look like agglomerates or clusters rather than regular bundles of filaments.Many fine filaments are about 60 or more times shorter than this "bundle-like cluster", and in addition, many of them are arc-bended, also showing a series of kinks.Magnification reveals reconnection between the filaments and their segments, axial merging, formation of thick filaments, core diffusion, as well as vortex rings.The double helix consists of fragmented magnetized filaments, which connects the blue central part and the brown-yellow peripheral part in the Crab Nebula.In contrast to bundles of the long filament jets in Figures 23 and 24, these "bundles" are composed of fine filaments (like threads) much shorter than the bundle itself; they look like erates or clusters rather than regular bundles of filaments.Many fine filaments are about 60 or more times shorter than this "bundle-like cluster", and in addition, many of them are arc-bended, also showing a series of kinks.Magnification reveals reconnection between the filaments and their segments, axial merging, formation of thick filaments, core diffusion, as well as vortex rings.

Vortex Ring Instabilities in the near Central Region of the Crab Nebula
Besides magnetized filaments, filament bundles, and clusters, the arc-bended segments, broken rings, and open-loop structures appear in the near central crab region.Various types of ring instabilities similar to those observed in multipulse LMI, including parametric oscillations and deformation Kelvin-like modes, can be identified.Many of these structures showing pinching are similar to the vortex ring instabilities in the near peripheral region of the very weak magnetic field described above and to the ring-like structures generated by multipulse LMI.This similarity is especially seen in the rings and broken rings (Figure 26a-d).Various types of ring instabilities similar to those observed in multipulse LMI, including parametric oscillations and deformation Kelvin-like modes, can be identified.Many of these structures showing pinching are similar to the vortex ring instabilities in the near peripheral region of the very weak magnetic field described above and to the ring-like structures generated by multipulse LMI.This similarity is especially seen in the rings and broken rings (Figure 26a-d).The double helix consists of fragmented magnetized filaments, which connects the blue central part and the brown-yellow peripheral part in the Crab Nebula.In contrast to bundles of the long filament jets in Figures 23 and 24, these "bundles" are composed of fine filaments (like threads) much shorter than the bundle itself; they look like agglomerates or clusters rather than regular bundles of filaments.Many fine filaments are about 60 or more times shorter than this "bundle-like cluster", and in addition, many of them are arc-bended, also showing a series of kinks.Magnification reveals reconnection between the filaments and their segments, axial merging, formation of thick filaments, core diffusion, as well as vortex rings.

Vortex Ring Instabilities in the near Central Region of the Crab Nebula
Besides magnetized filaments, filament bundles, and clusters, the arc-bended segments, broken rings, and open-loop structures appear in the near central crab region.Various types of ring instabilities similar to those observed in multipulse LMI, including parametric oscillations and deformation Kelvin-like modes, can be identified.Many of these structures showing pinching are similar to the vortex ring instabilities in the near peripheral region of the very weak magnetic field described above and to the ring-like structures generated by multipulse LMI.This similarity is especially seen in the rings and broken rings (Figure 26a-d).Keeping in mind the morphological and topological analogy with the vortex ring instabilities in LMI, one can say that the onset of the oscillations of the vortex ring and ring-like structures in the near central crab region corresponds to the onset of instabilities with mode m = 0 and m = 1, which deform the ring or cause its breakup.However, the onset of the core thinning (radial contraction) of the magnetized filament and its deformation may lead to other types of instability, including the formation of kinks, arc-shaped kinks, and open loops.The evolution of kinks with arch shape is observed on the magnetized above-mentioned "South-East (SE) jet" of the Crab Nebula.Such instabilities cause the SE jet to change the direction of propagation [65].The characteristic feature that repeats along the jet represents the pinching intervals corresponding to the generation of magnetized shock waves.The evolution of this 3D-magnetized SE jet from the Crab pulsar and the formation of kink were studied by three-dimensional relativistic magnetohydrodynamic (MHD) numerical simulations [65].
High-energy-density plasmas with strong magnetic fields and electric current and filamentary-jet structures similar to that of the Crab Nebula can be created by ultraintense lasers.A special design of the scaled laser experiments can help quantify RT/RM dynamics of the matter ejecta due to the magnetized plasma fluid acceleration and the jet instabilities.The laser-driven 3D-magnetized plasma jet modeled with three-dimensional numerical simulations can show the evolution of the kink instability.Such a study of the kink instability evolution on the SE jet was performed by the simulation-scaled experiment driven by an ultraintense laser of 10 14 W/cm 2 [5] It has created two beams that produced two plasma plumes; their collision generated a high Mach number plasma jet of velocity ~10 6 m/s and of very high temperature ~300 eV ~3.5 × 10 6 K.With the jet propagation, the onset of instabilities causes variation in the propagation direction and the kink behavior of magnetized jet [5].
In the simulated laser experiment, the Reynolds number was Re ~2 × 10 3 , and the magnetic Reynolds number was Re mag ~3 × 10 3 ; at the scale of the Crab Nebula Re ~2 × 10 17 and Re mag ~1 × 10 22 [5].Laser-created plasma jets had the specific configuration of self-generated spontaneous magnetic field (B = 10 6 G = 100 T), which consisted of poloidal (B p ) and toroidal (B ϕ ) fields [5].Magnetic fields and current-driven MHD instabilities in the jet were visualized and measured by proton radiography."With the jet propagation, plasma instabilities cause multiple deflections of the propagation direction, resembling the kink instability of the Crab jet" [5].
Various kinds of magnetized jet instability-including kink instability-are caused by the interplay of radial, azimuthal, and toroidal components of the magnetic field of the jet.The kinks evolve on the plasma jet at the places where the toroidal B φ component of the magnetic field becomes weak.In this scenario, the jet is confined by a toroidal magnetic field and accelerated outwards by the magnetic fields [5].The kink instabilities are stabilized by high jet velocity, meaning that instabilities alter the jet orientation, but do not disrupt the overall structure of the jet [5].We present the formation of the kink on the magnetized jet, which keeps its direction of propagation and takes the serpentine shape in the toroidal magnetic field (Figure 27).
Magnetic field (B) of the jet consists of poloidal (B p ) and toroidal (B ϕ ) fields: B = B p + B ϕ , where B p = B R + B z , i.e., the poloidal field consists of radial and azimuthal components as shown in Figure 28a.
Proton radio graphics revealed that the jet "is collimated throughout its propagation but has a sequence of clumps and changes of direction along its length.These features reflect perturbations in the magnetic field structure around the jet"; they grow locally and expand at each axial position, where the jet is unstable.The shape of the jet becomes serpentine due to the kink instability, creating a chain of kinks [5], and thus similar to that of the magnetized "South-East (SE) jet" of the Crab Nebula.As concluded in [5], the jet filament shows two low-order fast-growing MHD currentdriven instabilities with an increase in the B ϕ toroidal component.Mode m = 0 is initiated by an increase in B ϕ tension due to the radial contraction caused by vortex filament core thinning and the axis pinching when the criterion B φ > 2|B p is fulfilled (Figure 28b)."Mode m = 1, kink, is initiated when the strength of the pressure B ϕ increases at the inside of the kinks and decreases outside" [5] (Figure 28c).In that case, the product of toroidal and poloidal components becomes larger of some critical value α.The exact expression can be written as B φ B p −1 λ(2πr j ) −1 > α , where α is the criterion for the kink instability, λ is the jet modulation wavelength, and r j is the radius [5].
Regarding the astrophysical plasmas where the expansion of plasmatic matter "can tend to stabilize the jet, resulting in α larger than of order unity.The plasma jet is stabilized when the magnetic field is overwhelmed by the parallel components as the toroidal components around the jet are too weak to excite the MHD instability.When the field is sufficiently large and has nonuniform toroidal components B φ , current driven MHD kinks are excited" [5].
The fields are embedded in ("frozen-in") and advected with the fast-moving magnetized plasma flow and give an insight into the kink instability evolution on the South-East Crab jet.According to Li et al. [5], "the plasmatic vortex filament-jets are caused by current-curing magnetic filaments in which spontaneously formed currents, magnetic fields, and plasma flow vectors are all collinear".It is important to emphasize that in the laser-generated plasma, the jet kink instabilities are stabilized by high jet velocity, meaning that instabilities alter the jet orientation, but do not disrupt its overall structure-in agreement with the jet characteristics of the South-East Crab jet [5].
Therefore, the kink mode on the magnetized jet-which is similar to the bump-is initiated when the strength and pressure of B φ increase at the inside of the kinks and decrease outside.

Double-Helix Nebula: Formation of Filament "Complex"
Different filamentary organization-a double-helix "complex"-at the astrophysical scale can be seen in the image of a Double-Helix Nebula (DHN), which is an infrared (IR) emission nebula created by filaments of ionized molecular gasses and dust, which under torsional magnetic field form a double-helix spiral similar to DNA molecule (Figure 29).This feature is understood as a torsional Alfvén wave propagating from the Galactic disk and driven by rotation of the magnetized circumnuclear gas disk.It is proposed that the DHN is a hydrodynamic torsional Alfvén wave propagating vertically out of the Galactic plane, along magnetic field lines, from the near vicinity of the Galactic center This feature is understood as a torsional Alfvén wave propagating from the Galactic disk and driven by rotation of the magnetized circumnuclear gas disk.It is proposed that the DHN is a hydrodynamic torsional Alfvén wave propagating vertically out of the Galactic plane, along magnetic field lines, from the near vicinity of the Galactic center [106] (Figure 29).
The IR image of the nebula is colored to be observed in visible light.The nebula is about ~80 ly long and at 300 ly from the huge black hole thought to be at the galaxy's center.The formation of a double-helix shape of the nebula is caused by twisting magnetic field lines originating at the center of the Milky Way Galaxy, where the magnetic fields are extremely strong.Magnetic lines of this field are oriented perpendicular to the plane of the galaxy.It is speculated that these strong magnetic fields are created by the massive disc of gas circulating around the central supermassive black hole.This huge disk of gas (which orbits the black hole once every 10,000 years) is known as the circumnuclear disk (CND).It is strongly magnetized, and its predominant shear-induced azimuthal field is believed to merge with the ambient vertical galactic center field [107][108][109] The magnetic field lines are anchored in the CND and assumed to be the origin of the twisting of the magnetic field lines.These lines that are twisted at their base send a torsional wave in space.The magnetic field lines affect the orbits of molecular clouds and direct wind of hot gas and dust particles away from the central region.This twisting torsional wave carrying small dust particles and propagating very fast (about 1000 km/s) creates filaments that can absorb and emit infrared radiation giving rise to the Double-Helix Nebula shape.According to Morris et al. [106] the emission from the double helix most probably arises from thermal dust emission and molecular gas composition of filaments.The best-fit color temperature is ~630 K, while the thermal emission from dust gives ~410 K [106].
Large magnification of the image in Figure 29 reveals loops and kinks on the gaseous filament strands that make a double helix.Magnification reveals that filament strands are not long single filaments but eventually a group of filaments-"filament clusters".Even more, some of the tiny filaments in the clusters are not firmly and closely incorporated into the common structure, but although slightly shifted-they follow the same organization principle.Thus, the "filaments" may eventually be taken as "filament clusters", so that Double-Helix Nebula may be assumed a double-helix structure of filament clusters.
The double helix of vortex filament bundles (DHV)-which resembles the Double-Helix Nebula (DHN) of filament clusters-can be created by multipulse interaction of intense power laser with a Co-coated steel target.While the filament cluster in DHN is formed of "cold" dust particles at T = 400-600 K along helical magnetic lines, the DHV is formed by helical twisting of very hot fluid at T~10 4 K.Although these double helix structures are created by quite different physical processes, they form plasmatic structures and their self-organization with the same (or similar) morphological and topological characteristics.

Formation of Complex Filamentary Structures in LMI
The double helix of vortex filament bundles can be created by the multipulse interaction of an intense power laser with a Co-coated steel target on a parallel scratched target surface.Such surface creates initial conditions with multimodal perturbation of the shear layer, which behaves, in a certain stage of the nanosecond LMI, like viscoelastic fluid.Normal stress difference in shearing flows is a fundamental property of viscoelastic fluids [105] The long filaments and ribbons with strongly anisotropic order parameter texture within the vortex sheet are characteristic of the anisotropic fluid mobility and are consistent with driving oscillatory strain field.This strain field changes in intensity and topology with the number of pulses N and causes a cascade of vortex filament and ribbon transformations [60,105].Oscillatory strain fields cause motion and self-organization of filaments into helically paired, braided, and tangled patterns with topological complexity depending on the filament-filament separation distance and on the number of pulses.Interaction of vortex filaments with these strain fields generates various complex structures-"complexes"-in different regions of the target.Generation of localized coherent filamentary structures associated with turbulent mixing is a general phenomenon in accelerated flow environments from the atomic to nano-and micro-scales up to astrophysical scale as a result of nonlinear and nonequilibrium dynamics.consistent with driving oscillatory strain field.This strain field changes in intensity and topology with the number of pulses N and causes a cascade of vortex filament and ribbon transformations [60,105].Oscillatory strain fields cause motion and self-organization of filaments into helically paired, braided, and tangled patterns with topological complexity depending on the filament-filament separation distance and on the number of pulses.Interaction of vortex filaments with these strain fields generates various complex structures-"complexes"-in different regions of the target.Generation of localized coherent filamentary structures associated with turbulent mixing is a general phenomenon in accelerated flow environments from the atomic to nano-and micro-scales up to astrophysical scale as a result of nonlinear and nonequilibrium dynamics.

Helically Paired, Braided, and Tangled Structures of Vortex Filaments in LMI
Vortex filament bundles in Double-Helix Nebula organized along the helical magnetic field lines are very similar to double-helix vortex filament structures created by the power laser on the Co-coated steel target in the areas of a high density of vortex filaments.Under twisting, an array of closely packed parallel vortex filaments creates helices, double helices, and braided and tangled "complexes" such as the one shown in the center of Figure 30.Helically paired and braided vortex filaments generated by multiple LMI are more complex self-organized structures than those generated by the single-pulse interaction.In fact, the first pulse in a series generates vortex filaments in a parallel array or the unbraided braid in the terminology of Tufillaro et al. [110] The subsequent pulses cause motion of vortex filaments, which results in braiding, the complexity of which varies over the interaction space.Vortex filaments generated in some surface region either loop through the surface to the position of another filament or simply step to it and continue in the same direction to another surface element.As these filaments are forced to move about each other under a series of laser pulses, the loops become entangled with the same topological structure.In many cases, a set of vortices intervene with each other while staying more or less parallel.At some places, fast reconnection and axial merging of nearby filaments disturb the topologic structure.Thus, filament dynamic is inhomogeneous in the interaction space during multipulse LMI, creating the complex double-helix structure of many vortex filament bundles, which actually reveals all the above-mentioned interactions as can be seen in Figure 31a.
In the above figure, one can identify the 3D vortex reconnection process seen at various stages as well as the series of vortex solitons.To analyze these structures, taking into account all involved processes would be an enormous task.Very idealized theoretical work has to be performed to relate these structures to vortex filaments and their multistage evolutions.This approach ignores (for simplicity) many physical processes, e.g., shock wave baroclinic vortex deposition, viscous dissipative processes, thermal transport, and change of state.However, a simplified approach enables insights into the structure of this "complex" and enables its topological analysis with transformation, which represents it-in the mathematical sense-by the irreducible braid.
Vortex filaments in bundles are very close and tend to axial merging with diffuse core spreading.However, some filaments that are slightly separated from the bundles make a series of random loops, kinks, and other deformations, which-due to diffuse core spreading-appear misty.Approximate reconstruction of this "complex"-assuming that the core size of individual filaments is σ ~5-6 µm-reveals 10-11 vortex filaments organized into 5 bundles (Figure 31b).fact, the first pulse in a series generates vortex filaments in a parallel array or the unbraided braid in the terminology of Tufillaro et al. [110] The subsequent pulses cause motion o vortex filaments, which results in braiding, the complexity of which varies over the in teraction space.Vortex filaments generated in some surface region either loop through the surface to the position of another filament or simply step to it and continue in the same direction to another surface element.As these filaments are forced to move abou each other under a series of laser pulses, the loops become entangled with the same top ological structure.In many cases, a set of vortices intervene with each other while staying more or less parallel.At some places, fast reconnection and axial merging of nearby fil aments disturb the topologic structure.Thus, filament dynamic is inhomogeneous in the interaction space during multipulse LMI, creating the complex double-helix structure o many vortex filament bundles, which actually reveals all the above-mentioned interac tions as can be seen in Figure 31a.In the above figure, one can identify the 3D vortex reconnection process seen a various stages as well as the series of vortex solitons.To analyze these structures, taking into account all involved processes would be an enormous task.Very idealized theoreti The reconstruction reveals that some filaments that form bundles at a larger distance become split and then joined to another bundle.An example is the bundle composed of light-blue and dark-blue filaments, which are very close to the left side; approaching the right side of the diagram, this bundle becomes separated into two independent bundles that have different pairing dynamics.Similar behavior can be seen on the black bundle, which after the formation of two loops also splits into two different bundles (Figure 31b).

Topological Transformations of Braided Vortex Filaments
To analyze the characteristics of a double helix in the reconstructed braid of vortexfilament strings in Figure 31b, the topological operations should be applied to reach the "irreducible structure of the braid".Successive operations start with untwist-called the Reidemeister move type I-which gives the new braid.The next operations, pulling apart and sliding through the middle-Reidemeister moves II and III, respectively-result in the final braid diagram.These different planar diagrams of the same knot represent topologically equivalent knots [110,111].Reidemeister moves of types I, II, and III simplify the study of knot equivalence by reducing it to a two-dimensional problem.The type I move untwists a section of a string, the type II move pulls apart two strands, and the type III move acts on three strings sliding the middle strand between the outer strands [110,111].
According to the braid theory, any oriented link of the strings can be represented by a closed braid.The identification between links, braids, and the braid group allows us to pass back and forth between the geometric braids and the algebraic braid group.A geometric braid is constructed between two level lines with n base points chosen on these level lines.If the i + 1st strand passes over the i-th strand, there is a positive crossing between the two strands.An overcrossing (or right-crossing) between the i + 1st and i-th string is denoted by the symbol σ i .The inverse σ i −1 represents an undercross (or left-cross), i.e., the i + 1st strand goes under the i-th strand.
A general n-braid can be built up from successive applications of the operators σ i and σ i −1 .The first crossing between the second and third strands is represented by the operator σ 2 .The next crossing between the third and fourth strands is negative, σ 3 −1 , while the last crossing is represented by the operator σ.Each geometric diagram for a braid is equivalent to an algebraic braid word constructed from the overcrossing and undercrossing operators.The construction of a braid word is based on the following conventions.In each operation (σ i ,σ i −1 ), only the i-th and i + 1st strands are involved-not any other; it is not the string, but the base point which is numbered [110,111].By the application of these topological operations, the reconstructed braid in Figure 31b is first-by the untwist operation-transformed into the new equivalent geometric braid (Figure 32a).The successive operations, pulling apart and sliding through the middle, transform this braid into another new equivalent braid (Figure 32b).Finally, the identification of overcrossings σ i and undercrossings σ i −1 of the strings in this diagram makes it possible to construct the algebraic braid word.Figure 32c shows the final reconstruction of the braid "complex" with 35 over-and undercrossings of five strings represented by the algebraic braid word.
Since braiding is not homogeneous, the strings (4,1) and (4,2) form (irregularly) alternating segments of dense and rare helical pairing, while the string ( 5) is not paired except with the string (1).
Having identified braid words and transferring them into matrix form (matrix representation of braid) enable one to define the set of differential equations, which describe filament dynamics.This would be a serious task for the braid (23) because its structure does not offer the possibility for the formation of smaller simple block matrices.Namely, helical pairing alternates with variable phase and wavelength along the braid axis; a dense and rare pairing alternate in opposite phases somewhat irregularly, so that the braid (Figure 32b,c) is not symmetric, and its theoretical consideration is out of the scope of this paper.
However, we can mention some important topological characteristics of such braids.The topological structure of braided vortex filaments is created by their twisting caused by vertical motions.Because twist commutes with other structures, the opposite senses of twist on a tube cancel [112] The crucial role is played by the commutativity of twist and the relaxation of braids into minimal patterns-irreducible topologic structures [111,112]."Uniform twists commute with all other braid structures.However, more complex braid structures generally do not commute.In principle, the structures do have same freedom to relax to minimal patterns -irreducible structures, corresponding to minimum energy states" [112].
The braids generated by multipulse LMI are mostly random braids frozen (solidified) in a turbulent state, i.e., not in the state of the minimum energy [111,112].It happens that the flux within one element can break, loose, and wander across the surface, but not far from its original location.If the elements break up and reform before the sequence of vorticity changes sign, then the twist will be trapped within a more complicated braid structure, preventing cancelation [111,112].
With regard to the micrograph shown in Figure 31a, it can be analyzed and described at two levels: (i) at the level of filaments of decomposed braids and represented by the irreducible graph in Figure 32b,c and (ii) at the level of braids which are helically paired giving rise to the double-helix structure similar to that of the Crab Nebula.

Conclusions
This study of plasma instabilities brings together the results of LMI at the microscale and plasmatic structures in astrophysical systems at the megascale.Demonstrating a variety of plasma fluid flow instabilities in LMI, their morphologic and topologic characteristics are used as a reference for the classification of similar megascale instabilities in the Crab Nebula and Double-Helix Nebula.
Using multipulse laser interaction, we created plasma-fluid-layer instability, which resulted in the formation of RTI, vortex filaments, loop solitons, and the loop-soliton chains on vortex filaments.They are similar to those in the near peripheral region of the Crab Nebula, where the loop-soliton chains evolve due to helicoidal magnetic field lines.
We found that the same Hasimoto-type function that describes the formation of the loop solitons (caused by torsion of vortex filaments) also describes the formation of the loop solitons on magnetized plasma filaments (caused by helical winding of magnetic field lines) in the NPR of the Crab Nebula.
Our findings reveal that multipulse laser irradiation causes breakup of the loop solitons, creates vortex rings, and creates four types of vortex ring instabilities.We demonstrate, for the first time, that morphological and topological characteristics of these ring instabilities at the microscale are similar to those at the megascale and establish the conceptual frame for their unique classification based on the hierarchical order of Bessel functions.This classification reveals that two types of vortex ring instabilities belong to the parametric resonance of helicoidal modes, while the other two are a combination of quasistatic deformation modes of Kelvin waves on the vortex ring.We also demonstrate that although all four types of instabilities are present, only the bulge wave is dominant in the peripheral region of the Crab Nebula.
Our results also demonstrate that multipulse LMI may create a quite different type of instability, critically depending on the filament-filament distance.Governed by the characteristics of the plasma-fluid layer, where the vortex filament distance is small (about one core radius), the instability is initiated leading to the formation of braids and helical pairing.In such regions, the created conditions favor the motion of vortex filaments close to each other, leading to the onset of collective behavior and the formation of bundles.Subjected to twisting and torsional strain fields, bundles become helically paired while some-almost free vortex filaments in the bundles-form loop solitons creating "complexes".The formation and behavior of filament-bundle "complexes" in the shock-induced environment are described for the first time.We also show that these complexes are very similar to the magnetized jet filaments, filament bundles, and helically paired bundles-actually clusters-in the near central region of the Crab Nebula.Strong MHD radial and azimuthal instabilities on the filament jets in the bundles cause the formation of kinks, open loops, rings, and ring instabilities in the Crab Nebula.
Our overall findings show that the evolution of complexes may create a variety of filamentary architectures.Some of the complexes of closely packed vortex filaments tend to form a double-helix structure (DHS), which may be compared with a Double-Helix Nebula (DHN) bundle structure.Topological transformations of (symmetric) double-helix bundles lead to the irreducible braid graph, which reveals that DHN is actually the braided filamentary structure of much lower symmetry.

Figure 1 .
Figure 1.Figure 1. Characteristics of laser-generated plasma on a solid target.Schematic diagram of expanding laser-produced plasma in a background gas.The plasma plume, first, expands vertically to the target; colored zones illustrate variations in temperature and density in the expanding plasma plume.

Figure 1 .
Figure 1.Figure 1. Characteristics of laser-generated plasma on a solid target.Schematic diagram of expanding laser-produced plasma in a background gas.The plasma plume, first, expands vertically to the target; colored zones illustrate variations in temperature and density in the expanding plasma plume.

Figure 2 .
Figure 2. Plasma spheroid expansion along the target surface.Plasma density distribution after radial expansion with the formation of characteristic regions below Gaussian laser power density distribution.Schematically, (a) a cross-section of the plasma spheroid in the (z,x) plane; lateral compression of background gas (air), and formation of the air shell (cocoon) around the plasma spheroid.(b) Formation of central minimum in the plasma density distribution due to expansion and compression of the background gas (ρ = plasma density; CR = central region; NCR = near central region; Y-direction is parallel to the target surface and perpendicular to the plane of the paper).

Atoms 2024 , 43 Figure 3 .
Figure 3. Optical micrograph showing an array of parallel vortex filaments generated by N = 12 pulses, on a Co-coated steel target in the air, at the atmospheric pressure by XeCl laser with rectangular beam and the top hat profile; E = 0.25 J; τ = 16 ns (filament thickness ~5-7 μm; loop diameter 2R~28-30 μm).

Figure 3 .
Figure 3. Optical micrograph showing an array of parallel vortex filaments generated by N = 12 pulses, on a Co-coated steel target in the air, at the atmospheric pressure by XeCl laser with rectangular beam and the top hat profile; E = 0.25 J; τ = 16 ns (filament thickness ~5-7 µm; loop diameter 2R~28-30 µm).Surface patterning and multipulse LMI create a series of shock waves of different wavevectors, phases, and intensities.Two kinds of vortex filament arrays are formed; the

Figure 4 .
Figure 4. Loop soliton on vortex filament.Optical micrograph of the loop soliton on Cr-coated steel generated in the air (a); SEM micrograph of the loop soliton on Titanium target generated in vacuum (b).

Figure 4 .
Figure 4. Loop soliton on vortex filament.Optical micrograph of the loop soliton on Cr-coated steel generated in the air (a); SEM micrograph of the loop soliton on Titanium target generated in vacuum (b).

Figure 5 .
Figure 5. Optical micrograph of the chain of vortex filament loops (Hasimoto solitons) generated on Co-coated steel target.For a simple illustration of the multisoliton dynamics, we present the motion of 2D loop solitons of different radii r1 and r2, with the curvatures

Figure 6 .
Figure 6.Two-soliton solution for k1 = 3, k2 = 1, t = −0.30(a); t = −0.045(b); t = 0.30 (c).Interaction o two loop solitons.Small loop soliton moving at high velocity approaches to the large slower one (a).Reaching the larger soliton, it becomes incorporated into the larger soliton (b); overtaking the large soliton, it continues to move at a higher velocity (c).
in nondimensional time.Rather long expressions for 1 κ and

Figure 5 .
Figure 5. Optical micrograph of the chain of vortex filament loops (Hasimoto solitons) generated on Co-coated steel target.

Figure 6 .
Figure 6.Figure 6. Two-soliton solution for k 1 = 3, k 2 = 1, t = −0.30(a); t = −0.045(b); t = 0.30 (c).Interaction of two loop solitons.Small loop soliton moving at high velocity approaches to the large slower one (a).Reaching the larger soliton, it becomes incorporated into the larger soliton (b); overtaking the large soliton, it continues to move at a higher velocity (c).

Figure 6 .
Figure 6.Figure 6. Two-soliton solution for k 1 = 3, k 2 = 1, t = −0.30(a); t = −0.045(b); t = 0.30 (c).Interaction of two loop solitons.Small loop soliton moving at high velocity approaches to the large slower one (a).Reaching the larger soliton, it becomes incorporated into the larger soliton (b); overtaking the large soliton, it continues to move at a higher velocity (c).

Atoms 2024 , 4 Figure 7 .
Figure 7. Array of closely packed parallel quasi-periodic loop soliton chains-a quasiperiodic la tice.The strings are disturbed and indicate a tendency to catastrophic dynamics of the loop soliton after N = 12 to 13 pulses.

Figure 7 .
Figure 7. Array of closely packed parallel quasi-periodic loop soliton chains-a quasiperiodic lattice.The strings are disturbed and indicate a tendency to catastrophic dynamics of the loop solitons after N = 12 to 13 pulses.

Figure 8 .
Figure 8. Breakdown of 1D loop chain into a set of open loops and vortex rings on Co-coated steel target.Vortex rings after breakdown do not move much from the former positions in the string.

Atoms 2024 ,
12, x FOR PEER REVIEW 15 of 43 space."This modulation is different in different regions due to in-homogeneity in the distribution of local strain-field symmetry.As a consequence, various types of instability are created in different regions of the quasi-linear array of rings" [58].

Figure 9 .
Figure 9. Type I: Vortex ring instability identified as the "WBT instability", parametric resonance between the left-and right-handed helical waves (m = ±1).The evolution of type I instability.The radius a of the core cross-section grows as (a) a~2.5 μm, and (b) a schematic illustration.

Figure 9 .
Figure 9. Type I: Vortex ring instability identified as the "WBT instability", parametric resonance between the left-and right-handed helical waves (m = ±1).The evolution of type I instability.The radius a of the core cross-section grows as (a) a~2.5 µm, and (b) a schematic illustration.

Figure 10 .
Figure 10.Type II: Vortex ring instability identified as the "curvature instability", parametric resonance between the helical (m = 1) and the bulge (m = 0) modes.This instability may be represented as a combination of the axisymmetric (m = 0) and the bending (m = 1) modes (a).A schematic illustration of this type of instability is shown in (b).

Figure 11 .
Figure 11.Type III: Vortex ring oscillation identified as the Kelvin axisymmetric mode (m = 0).Schematic illustration (a-c) shows a standing deformation mode of the core as provided by Fig-ure10 of Kop'ev and Chernyshev[75].Illustration (i-iv) shows the stages of the deformation mode.

Figure 10 .
Figure 10.Type II: Vortex ring instability identified as the "curvature instability", parametric resonance between the helical (m = 1) and the bulge (m = 0) modes.This instability may be represented as a combination of the axisymmetric (m = 0) and the bending (m = 1) modes (a).A schematic illustration of this type of instability is shown in (b).

Figure 10 .
Figure 10.Type II: Vortex ring instability identified as the "curvature instability", parametric resonance between the helical (m = 1) and the bulge (m = 0) modes.This instability may be represented as a combination of the axisymmetric (m = 0) and the bending (m = 1) modes (a).A schematic illustration of this type of instability is shown in (b).

Figure 11 .
Figure 11.Type III: Vortex ring oscillation identified as the Kelvin axisymmetric mode (m = 0).Schematic illustration (a-c) shows a standing deformation mode of the core as provided by Fig-ure10 of Kop'ev and Chernyshev [75].Illustration (i-iv) shows the stages of the deformation mode.

Figure 11 .
Figure 11.Type III: Vortex ring oscillation identified as the Kelvin axisymmetric mode (m = 0).Schematic illustration(a-c) shows a standing deformation mode of the core as provided by Figure10of Kop'ev and Chernyshev[75].Illustration (i-iv) shows the stages of the deformation mode.

Figure 13 .
Figure 13.Optical micrograph showing a bundle of vortex filaments and a string of solitons created by laser on Co-coated steel (a).Rendering of (a) reveals the vortex filament "complex" showing multiloop solitons on bundles (b).While some vortex filaments are incorporated into bundles, those on the outer side of the bundles show individual motion.

Figure 14 .
Figure14.Schematic reconstruction of the vortex filament "complex" in Figure14.Formation of two bundles of vortex filaments (i).Turbulent background field causes dense packing of filaments; twisting of bundles causes helical winding of individual bundles but also their helical pairing (ii).The peripheral filament in the bundle experiences additional torsion, which starts to form a series of loop solitons (iii).Closely packed filaments under the action of the surrounding turbulent background field experience core spreading and axial merging.

Figure 13 .
Figure 13.Optical micrograph showing a bundle of vortex filaments and a string of solitons created by laser on Co-coated steel (a).Rendering of (a) reveals the vortex filament "complex" showing multiloop solitons on bundles (b).While some vortex filaments are incorporated into bundles, those on the outer side of the bundles show individual motion.

Figure 16 .
Figure 16.Reconstruction of the vortex filament "complex" in Figure 16 showing (+) and (−) soliton loops and filament coiling.Two connected rings represent the Hopf link created after loop soliton breakup and the formation of vortex rings.(a) The ring instability causes a breakup and anew connection with the nearest neighbor ring.(b) Schematic illustration of the interaction of (+) and (−) loop solitons (i) and their annihilation (ii).

Figure 16 .
Figure 16.Reconstruction of the vortex filament "complex" in Figure 16 showing (+) and (−) soliton loops and filament coiling.Two connected rings represent the Hopf link created after loop soliton breakup and the formation of vortex rings.(a) The ring instability causes a breakup and anew connection with the nearest neighbor ring.(b) Schematic illustration of the interaction of (+) and (−) loop solitons (i) and their annihilation (ii).

Figure 17 .
Figure 17.Optical micrograph of vortex filament "complex" showing the series of disarrayed Hasimoto solitons on two tangled vortex filaments.Small loop solitons are formed at the thin vortex filament separated from the thicker one.This indicates that thin and thick filaments under torsion have different dynamics although both form loop solitons.The structure is generated on the Co-coated-steel-plate in the air.

Figure 17 .
Figure 17.Optical micrograph of vortex filament "complex" showing the series of disarrayed Hasimoto solitons on two tangled vortex filaments.Small loop solitons are formed at the thin vortex filament separated from the thicker one.This indicates that thin and thick filaments under torsion have different dynamics although both form loop solitons.The structure is generated on the Co-coated-steel-plate in the air.

Figure 18 .
Figure 18.Crab Nebulain the Taurus constellation(as the large mosaic of the Crab Nebula assembled from 24 individual exposures captured by Hubble Space telescope over three months).Image courtesy of NASA/ESA and Allison Loll/Jeff Hester (Arizona State University).Acknowledgment: Davide De Martin (ESA/Hubble).

Figure 19 .
Figure 19.The upperright peripheral region of the Crab Nebula showing random loop solitons on magnetized vortex filaments.(Enlarged segment of the Crab Nebula in Figure 18.Image courtesy of NASA/ESA and Allison Loll/Jeff Hester (Arizona State University).Acknowledgment: Davide De Martin (ESA/Hubble).)

Figure 19 .
Figure 19.The upperright peripheral region of the Crab Nebula showing random loop solitons on magnetized vortex filaments.(Enlarged segment of the Crab Nebula in Figure 18.Image courtesy of NASA/ESA and Allison Loll/Jeff Hester (Arizona State University).Acknowledgment: Davide De Martin (ESA/Hubble).)

Figure 20 .
Figure 20.The upperleft peripheral region of the Crab Nebula showing loop solitons in the chaotic organization of vortex filaments and diffusively spreaded cores of the vortex filaments.(Enlarged segment of the Crab Nebula in Figure 18.Image courtesy of NASA/ESA and Allison Loll/Jeff Hester (Arizona State University).Acknowledgment: Davide De Martin (ESA/Hubble).)

Figure 20 .
Figure 20.The upperleft peripheral region of the Crab Nebula showing loop solitons in the chaotic organization of vortex filaments and diffusively spreaded cores of the vortex filaments.(Enlarged segment of the Crab Nebula in Figure 18.Image courtesy of NASA/ESA and Allison Loll/Jeff Hester (Arizona State University).Acknowledgment: Davide De Martin (ESA/Hubble).)

Figure 21 .
Figure 21.Series of almost periodic loops on vortex filament (chain of the loop solitons) on the segment of the Crab Nebula.Vortex rings separated from the parent filament are seen in the vicinity.The picture is taken from the segment of the Crab Nebula.(Enlarged segment of the Crab Nebula in Figure 18.Image courtesy of NASA/ESA and Allison Loll/Jeff Hester (Arizona State University).Acknowledgment: Davide De Martin (ESA/Hubble).)

Figure 21 .
Figure 21.Series of almost periodic loops on vortex filament (chain of the loop solitons) on the segment of the Crab Nebula.Vortex rings separated from the parent filament are seen in the vicinity.The picture is taken from the segment of the Crab Nebula.(Enlarged segment of the Crab Nebula in Figure 18.Image courtesy of NASA/ESA and Allison Loll/Jeff Hester (Arizona State University).Acknowledgment: Davide De Martin (ESA/Hubble).)6.2.Vortex Ring Instabilities in the Near Peripheral Region of the Crab NebulaPeripheral regions of the Crab Nebula under magnification also reveal deformed and broken rings as well as the open loops in the vicinity of the string of the loop solitons on filaments (Figure22).

Figure 22 .
Figure 22.The right peripheral region of the Crab Nebula.Broken filament loops and rings indicate the ring instability of type III (a,b); formation of clusters of vortex rings (c);and cluster of broken rings (d).(Enlarged segment of the Crab Nebula in Figure 18.Image courtesy of NASA/ESA and Allison Loll/Jeff Hester (Arizona State University).Acknowledgment: Davide De Martin (ESA/Hubble).)

Figure 22 .
Figure 22.The right peripheral region of the Crab Nebula.Broken filament loops and rings indicate the ring instability of type III (a,b); formation of clusters of vortex rings (c);and cluster of broken rings (d).(Enlarged segment of the Crab Nebula in Figure 18.Image courtesy of NASA/ESA and Allison Loll/Jeff Hester (Arizona State University).Acknowledgment: Davide De Martin (ESA/Hubble).)

Figure 23 .Figure 24 .
Figure 23.The double helix of gaseous magnetized bundles of filament jets in the near central region of the Crab Nebula.Some of the thin filament jets slightly separated from the bundles show the formation of loops.Filament jets are helically paired, and closed loops (similar to the loop solitons) are formed on a few filaments.(Enlarged segment of the Crab Nebula in Figure 18.Image courtesy of NASA/ESA and Allison Loll/Jeff Hester (Arizona State University).Acknowledgment: Davide De Martin (ESA/Hubble).)

Figure 23 .
Figure 23.The double helix of gaseous magnetized bundles of filament jets in the near central region of the Crab Nebula.Some of the thin filament jets slightly separated from bundles show the formation of loops.Filament jets are helically paired, and closed loops (similar to the loop solitons) are formed on a few filaments.(Enlarged segment of the Crab Nebula in Figure 18.Image courtesy of NASA/ESA and Allison Loll/Jeff Hester (Arizona State University).Acknowledgment: Davide De Martin (ESA/Hubble).)

Figure 23 .Figure 24 .
Figure 23.The double helix of gaseous magnetized bundles of filament jets in the near central region of the Crab Nebula.Some of the thin filament jets slightly separated from the bundles show the formation of loops.Filament jets are helically paired, and closed loops (similar to the loop solitons) are formed on a few filaments.(Enlarged segment of the Crab Nebula in Figure 18.Image courtesy of NASA/ESA and Allison Loll/Jeff Hester (Arizona State University).Acknowledgment: Davide De Martin (ESA/Hubble).)

Figure 24 .
Figure 24.Enlarged details of Figure 23 showing vortex filament loops and rings on the filaments not firmly incorporated into the bundle, which manifest an individual behavior.(Enlarged segment of Crab Nebula in Figure 18.Image courtesy of NASA/ESA and Allison Loll/Jeff Hester (Arizona State University).Acknowledgment: Davide De Martin (ESA/Hubble).)

Figure 25 .
Figure 25.The double helix of magnetized tiny filaments, which connects the blue central part and the brown-yellow peripheral part in the Crab Nebula.Large magnification reveals small loops.(Enlarged segment of the Crab Nebula in Figure 18.Image courtesy of NASA/ESA and Allison Loll/Jeff Hester (Arizona State University).Acknowledgment: Davide De Martin (ESA/Hubble).)

Figure 25 .
Figure 25.The double helix of magnetized tiny filaments, which connects the blue central part and the brown-yellow peripheral part in the Crab Nebula.Large magnification small loops.(Enlarged segment of the Crab Nebula in Figure 18.Image courtesy of NASA/ESA and Allison Loll/Jeff Hester (Arizona State University).Acknowledgment: Davide De Martin (ESA/Hubble).)7.2.Vortex Ring Instabilities in the Near Central Region of the Crab NebulaBesides magnetized filaments, filament bundles, and clusters, the arc-bended segments, broken rings, and open-loop structures appear in the near central crab region.Various types of ring instabilities similar to those observed in multipulse LMI, including parametric oscillations and deformation Kelvin-like modes, can be identified.Many of these structures showing pinching are similar to the vortex ring instabilities in the near peripheral region of the very weak magnetic field described above and to the ring-like structures generated by multipulse LMI.This similarity is especially seen in the rings and broken rings (Figure26a-d).

Figure 25 .
Figure 25.The double helix of magnetized tiny filaments, which connects the blue central part and the brown-yellow peripheral part in the Crab Nebula.Large magnification reveals small loops.(Enlarged segment of the Crab Nebula in Figure 18.Image courtesy of NASA/ESA and Allison Loll/Jeff Hester (Arizona State University).Acknowledgment: Davide De Martin (ESA/Hubble).)

Figure 26 .
Figure 26.Near the central region of the Crab Nebula.Open loops, broken rings, and closed rings appearing with bundles of magnetic filament jets in the region of the very magnetic field.The ring instability (a) is similar to the type I parametric instability; the ring instabilities in (b,c) are similar to the type III, i.e., to the axisymmetric "bulge wave" (m = 0) that has reached the catastrophic limit and caused the ring breakdown.The instability (d) is similar to the vortex ring instability type IV of the long wavelength ring oscillations.The instabilities (e,f) do not correspond to the hierarchical ring instabilities but rather indicate the open loops.(Enlarged segment of the Crab Nebula in Figure 18.Image courtesy of NASA/ESA and Allison Loll/Jeff Hester (Arizona State University).Acknowledgment: Davide De Martin (ESA/Hubble).)

7. 3 .
Simulation of the Kink Instability in Magnetized Jet of the Crab Nebula: Plasma Created by Ultraintense Lasers

Figure 27 .
Figure27.Schematic of a plasma jet with the resulting toroidal magnetic field configuration.Plasma jet propagates in this magnetic field configuration with the formation of kinks that cause its serpentine shape.The kinks evolve at the places-in the planes-where the toroidal B ϕ component of the magnetic field is weak.

Figure 29 .
Figure 29.Double-Helix Nebula is the IR emission nebula in Sagittarius.The image was taken by the Multiband Imaging Photometer for Spitzer (MIPS).(Image courtesy of NASA/JPL-Caltech/M.Morris (UCLA).Acknowledgment: Steward Wolpert, UCLA.)

9. 1 .
Helically Paired, Braided, and Tangled Structures of Vortex Filaments in LMI Vortex filament bundles in Double-Helix Nebula organized along the helical magnetic field lines are very similar to double-helix vortex filament structures created by the power laser on the Co-coated steel target in the areas of a high density of vortex filaments.Under twisting, an array of closely packed parallel vortex filaments creates helices, double helices, and braided and tangled "complexes" such as the one shown in the center of Figure 30.

Figure 30 .
Figure 30.Optical micrograph of a segment of vortex filaments that follow parallel surface scratches.Interaction between vortex filaments causes the formation of tangled filament bundles-the "complexes" (central region)-generated by XeCl laser on Co-coated steel by N = 13 pulses.Transversal oscillation causes motion and perturbation of filaments.

Figure 31 .
Figure 31.The double helix of vortex filament bundles generated by XeCl laser on Co-coated steel E = 0.25 J, N = 13, τ = 16-20 ns.Many filaments in the bundles are very close and show axial merging and the core spreading.Vortex filaments show strong diffuse core spreading in turbulent fields (a) Approximate reconstruction reveals ~11 vortex filaments organized into 5 tangled bundles, many of them with the loop solitons.A few bundles with distance from the left side show the splitting o some filaments, which then become part of another bundle (b).

Figure 31 .
Figure 31.The double helix of vortex filament bundles generated by XeCl laser on Co-coated steel: E = 0.25 J, N = 13, τ = 16-20 ns.Many filaments in the bundles are very close and show axial merging and the core spreading.Vortex filaments show strong diffuse core spreading in turbulent fields (a).Approximate reconstruction reveals ~11 vortex filaments organized into 5 tangled bundles, many of them with the loop solitons.A few bundles with distance from the left side show the splitting of some filaments, which then become part of another bundle (b).

Figure 32 .
Figure 32.Figure 32.Topological reconstruction of vortex filament bundles organized into a braid and identification of the algebraic braid word.Application of topological operations (Reidemeister move type I) creates the new braid (a).Subsequent operations pull apart and the middle slide (Reidemeister moves II and III) give the final braid shown in (b) with identified over-and undercrossings (c).

Figure 32 .
Figure 32.Figure 32.Topological reconstruction of vortex filament bundles organized into a braid and identification of the algebraic braid word.Application of topological operations (Reidemeister move type I) creates the new braid (a).Subsequent operations pull apart and the middle slide (Reidemeister moves II and III) give the final braid shown in (b) with identified over-and undercrossings (c).

Funding:
This work is supported by the Croatian Science Foundation under the project: IP-2014-09-7046.