CyberSec Research

Recent studies have applied variational calculus, conformal mapping, and point transformations to generalize the one-dimensional (1D) space-charge limited current density (SCLCD) and electron emission mechanisms to nonplanar geometries; however, these assessments have focused on extending the Child-Langmuir law (CLL) for SCLCD in vacuum. Since the charge in the diode is independent of coordinate system (i.e., covariant), we apply bijective point transformations to extend the Mott-Gurney law (MGL) for the SCLCD in a collisional or semiconductor gap to nonplanar 1D geometries. This yields a modified MGL that replaces the Cartesian gap distance with a canonical gap distance that may be written generally in terms of geometric scale factors that are known for multiple geometries. We tabulate results for common geometries. Such an approach may be applied to any current density, including non-space-charge limited gaps and SCLCD that may fall between the CLL and MGL.

Magnetic helicity is a quantity that underpins many theories of magnetic relaxation in electrically conducting fluids, both laminar and turbulent. Although much theoretical effort has been expended on magnetic fields that are everywhere tangent to their domain boundaries, many applications, both in astrophysics and laboratories, actually involve magnetic fields that are line-tied to the boundary, i.e. with a non-trivial normal component on the boundary. This modification of the boundary condition requires a modification of magnetic helicity, whose suitable replacement is called relative magnetic helicity. In this work, we investigate rigorously the behaviour of relative magnetic helicity under turbulent relaxation. In particular, we specify the normal component of the magnetic field on the boundary and consider the \emph{ideal limit} of resistivity tending to zero in order to model the turbulent evolution in the sense of Onsager's theory of turbulence. We show that relative magnetic helicity is conserved in this distinguished limit and that, for constant viscosity, the magnetic field can relax asymptotically to a magnetohydrostatic equilibrium.

Whether for materials processing or medical applications, the use of atmospheric pressure plasma jets (APPJs) has emerged as a relevant alternative to conventional methods. Within the APPJs research field, the search for innovation aims not only to solve existing problems, but also to explore novel options for generating plasma jets and find new possible applications. In this work, the properties of $\rm{Ar-H_2}$ APPJs generated using two plasma sources, which differ basically in the generated voltage frequency, amplitude and waveform, were studied through electrical, thermal and optical characterization. Discharge and plasma parameters were analyzed as a function of the $\rm{H_2}$ content in the gas mixture, with this parameter varying from $0\%$ to $3.5\%$. In all cases, the discharge power, electron density as well as the rotational, vibrational and gas temperatures presented a trend of growing when the proportion of $\rm{H_2}$ in the gas composition was increased. Optical emission spectroscopy revealed that the same reactive species were produced for both plasma sources, except for nitric oxide (NO), which was observed only for the one operated at higher frequency (PS #1). Applications on polymer (polypropylene, PP) and water treatment were performed using PS #1 without $\rm{H_2}$ and with $3.5\%$ of $\rm{H_2}$ in the gas mixture. NH functional groups were detected on the PP surface in the presence of $\rm{H_2}$ in the gas composition. This indicates a possible way to increase the nitrogen content on polymer surfaces. The results of water treatment revealed that ammonia ($\rm{NH_3}$) is also produced when there is $\rm{H_2}$ in the working gas. This opens an alternative for the use of plasma treated water in agriculture.

A novel concept called Air-Breathing Electric Propulsion proposes to fly satellites at altitudes in the range 180-250 km, since this would have some advantages for the performance of radio communication and Earth observation equipment. The ABEP satellites compensate the atmospheric drag through a continuous thrust provided by collecting, ionizing and accelerating the residual atmospheric particles. It is clear that the feasibility of this concept will require a significant design and testing effort, performed first on ground and later in orbit. Plasma simulation tools play a fundamental role in the development of this technology, for two main reasons: (i) they can potentially increase dramatically the optimization and testing process of ABEP systems, since on-ground testing and in-orbit demonstrators are costly and time consuming, and (ii) the fidelity of on-ground testing is limited by the finite size and pumping speed of high-vacuum facilities, as well as the means through which the orbital flow is produced. In this paper, we demonstrate a one-way coupled, particle-based simulation strategy for a CubeSat sized ABEP system. The neutral flow in the full geometry of the ABEP system comprising the intake and the thruster is simulated first through Direct Simulation Monte Carlo. Then, the resulting neutral density is used as the input for a Particle-in-Cell simulation of the detailed thruster geometry. The simulations are performed in 3D and within the VKI in-house code Pantera, taking advantage of the fully-implicit energy-conserving scheme.

We propose a laser-plasma wakefield based schemes for in situ axion generation and detection through the Primakoff process. Strong electromagnetic fields ($\gtrsim 10^{9}\,$V/cm) in the wakefield enhance axion production rates by orders of magnitude compared to conventional light-shining-through-wall (LSW) experiments. By replacing the axion generation stage with laser-wakefield interaction, one can achieve the axion-photon coupling constraints to the level of $g_{a\gamma\gamma}\sim 10^{-12}\,\text{GeV}^{-1}$. Besides, the generated axions can convert back into photons in the background field, leading to axion-regenerated electromagnetic fields (AREM) with unique polarization, frequency, and transverse distribution properties. This allows for effective filtering of the AREM from the background field, enhancing signal-to-noise ratios. This approach establishes plasma wakefields as a promising platform for laboratory axion searches.

Plasma electrolytic oxidation (PEO) is a technique used to create oxide-ceramic coatings on lightweight metals, such as aluminium, magnesium, and titanium. PEO is known for producing coatings with high corrosion resistance and strong adhesion to the substrate. The process involves generating short-lived microdischarges on the material surface through anodic dielectric breakdown in a conductive aqueous solution. To investigate single microdischarges during PEO, a single microdischarge setup was developed, where the active anode surface is reduced to the tip of a wire with a diameter of 1 mm. In this work the focus is on the effect of electrolyte concentration, anode material, and electrical parameters on the microdischarges. The electrolyte is composed of distilled water with varying concentrations of potassium hydroxide (0.5 - 4 g/l). High-speed optical measurements are conducted to gain insights into the formation and temporal evolution of individual microdischarges and the induced gas bubble formation. Optical emission spectroscopy is used to estimate surface and electron temperatures by fitting Bremsstrahlung and Planck's law to the continuum spectrum of the microdischarges. To evaluate the impact of the microdischarges on coating morphology, the resulting oxide layers on the metal tips are analysed using scanning electron microscopy. The study demonstrates that microdischarge behaviour is significantly influenced by the substrate material, treatment time, and electrolyte concentration, all of which impact the coating morphology. Under the conditions studied in this work, aluminium exhibits longer microdischarge and bubble lifetimes, with fewer cracks on the top layer of the coating, whereas titanium showed faster, shorter-lived bubbles due to more rapid microdischarge events.

First-principle studies of radiative processes aimed at explaining the origin of type II and type III solar radio bursts raise questions on the implications of downshifted electron beam plasma excitations with frequency (slightly) below the plasma frequency ($\omega\lesssim\omega_{pe}$) in the generation of radio emissions. Unlike the beam-induced Langmuir waves ($\omega \gtrsim \omega_{pe}$) in the standard radio emission plasma model, the primary wave excitations of cooler and/or denser beams have predominantly downshifted frequencies. Broadbands of such downshifted excitations are also confirmed by in situ observations in association with terrestrial foreshock and electron beams (in contrast to narrowband Langmuir waves), but their involvement in radiative processes has not been examined so far. We revisit three radiative scenarios specific to downshifted primary excitations, and the results demonstrate their direct or indirect involvement in plasma radio emission. Downshifted excitations of an electron beam primarily play an indirect role, contributing to the relaxation to a plateau-on-tail still able to induce Langmuir beam waves that satisfy conditions for nonlinear wave-wave interactions leading to free radio waves. At longer time scales, the primary excitations can become predominantly downshifted, and then directly couple with the secondary (backscattered) Langmuir waves to generate the second harmonic of radio emissions. Two counterbeams are more efficient and lead to faster radiative mechanisms, involving counterpropagating downshifted excitations, which couple to each other and generate intense, broadband and isotropic radio spectra of downshifted second harmonics. Such a long-lasting (second) radio harmonic can thus be invoked to distinguish regimes with downshifted ($\omega \gtrsim \omega_{pe}$) primary excitations.

We have developed TorbeamNN: a machine learning surrogate model for the TORBEAM ray tracing code to predict electron cyclotron heating and current drive locations in tokamak plasmas. TorbeamNN provides more than a 100 times speed-up compared to the highly optimized and simplified real-time implementation of TORBEAM without any reduction in accuracy compared to the offline, full fidelity TORBEAM code. The model was trained using KSTAR electron cyclotron heating (ECH) mirror geometries and works for both O-mode and X-mode absorption. The TorbeamNN predictions have been validated both offline and real-time in experiment. TorbeamNN has been utilized to track an ECH absorption vertical position target in dynamic KSTAR plasmas as well as under varying toroidal mirror angles and with a minimal average tracking error of 0.5cm.

Mixtures of bare atomic nuclei on a nearly uniform degenerate electron background are a realistic model of matter in the interior of white dwarfs. Despite tremendous progress in understanding their phase diagrams achieved mainly via first-principle simulations, structural, thermodynamic, and kinetic properties of such mixtures are poorly understood. We develop a semi-analytic model of the crystal state of binary mixtures based on the concept of mutual short-range ordering of ions of different sorts. We derive analytic formulas for electrostatic energy of crystal mixtures, including the effect of static ion displacements from the lattice nodes, and estimate their residual entropy. Then we perform free energy minimization with respect to the order parameters for a C/O mixture at all relevant compositions and temperatures. The resulting C/O phase diagram is in a reasonable agreement with that obtained in the most recent first-principle study. The equilibrium microstructure of a crystallized mixture is shown to evolve with decrease of temperature which, in principle, can induce structural transitions. The latter will be accompanied by thermal energy release. The proposed theory opens up a path to analyze ordering and construct phase diagrams of ternary mixtures, which are of great practical interest in astrophysics, as well as to improve calcuations of electron-ion scattering rates and kinetic properties of dense crystallized matter.

The preheat and pre-magnetization of the fuel are essential steps in the design of Magnetized Liner Inertial Fusion (MagLIF) configurations. Typically, the energy of the preheat laser is deposited in a central region of the fuel and propagates outward generating magneto-hydrodynamic structures that impact the fuel mass distribution and magnetic flux compression during the subsequent implosion. We present a theoretical analysis of preheat propagation in a magnetized plasma under conditions typical for MagLIF. The analysis is based on the acoustic time scale for the propagation of pressure disturbances being much shorter than the conductive time scale for heat diffusion. In this regime, the preheat-driven expansion induces the stratification of fuel mass and magnetic field, which accumulate in a dense outer shelf bounded by the leading shock. We derive self-similar solutions of the mathematical model that describe the hydrodynamic profiles of the expansion, and evaluate the evolution of the magnetic field in this configuration. The model is supported by FLASH simulations of preheat propagation. Our analysis shows that the regions where the magnetization of the fuel is significant tend to become localized asymptotically in time at the interface separating the outer shelf from the inner hot core. We assess the implications of this stratification on the magnetic flux conservation and performance of fully integrated MagLIF FLASH simulations.

The solar wind is a medium characterized by strong turbulence and significant field fluctuations on various scales. Recent observations have revealed that magnetic turbulence exhibits a self-similar behavior. Similarly, high-resolution measurements of the proton density have shown comparable characteristics, prompting several studies into the multifractal properties of these density fluctuations. In this work, we show that low-resolution observations of the solar wind proton density over time, recorded by various spacecraft at Lagrange point L1, also exhibit non-linear and multifractal structures. The novelty of our study lies in the fact that this is the first systematic analysis of solar wind proton density using low-resolution (hourly) data collected by multiple spacecraft at the L1 Lagrange point over a span of 17 years. Furthermore, we interpret our results within the framework of non-extensive statistical mechanics, which appears to be consistent with the observed nonlinear behavior. Based on the data, we successfully validate the q-triplet predicted by non-extensive statistical theory. To the best of our knowledge, this represents the most rigorous and systematic validation to date of the q-triplet in the solar wind.

Mendeleev's periodic table successfully groups atomic elements according to their chemical and spectroscopic properties. However, it becomes less sufficient in describing the electronic properties of highly charged ions (HCIs) in which many of the outermost electrons are ionized. In this work, we put forward a periodic table particularly suitable for HCIs. It is constructed purely based on the successive electron occupation of relativistic orbitals. While providing a much-simplified description of the level structure of highly charged isoelectronic ions -- essential for laboratory and astrophysical plasma spectroscopies, such a periodic table predicts a large family of highly forbidden transitions suitable for the development of next-generation optical atomic clocks. Furthermore, we also identify universal linear $Z$ scaling laws ($Z$ is the nuclear charge) in the so-called ``Coulomb splittings'' between angular momentum multiplets along isoelectronic sequences, complementing the physics of electron-electron interactions in multielectron atomic systems.

The recent ITER re-baselining calls for new fusion-relevant research best carried out in a DT-capable tokamak device with similar characteristics. The present paper describes key issues that could be addressed in a Suitably Enhanced DT-capable Tokamak (SET), with tungsten plasma facing components, boronization systems, and 10 MW of ECRH, based on characteristics and knowledgebase of JET. We discuss hardware options, and show that fusion-relevant operational scenarios could be achieved. Notably, development, validation and testing of fusion and nuclear diagnostics, to be used in next generation devices, would require a D-T capable tokamak as described.

Global consumption of heat is vast and difficult to decarbonise, but it could present an opportunity for commercial fusion energy technology. The economics of supplying heat with fusion energy are explored in context of a future decarbonised energy system. A simple, generalised model is used to estimate the impact of selling heat on profitability, and compare it to selling electricity, for a variety of fusion proposed power plant permutations described in literature. Heat production has the potential to significantly improve the financial performance of fusion over selling electricity. Upon entering a highly electrified energy system, fusion should aim to operate as a grid-scale heat pump, avoiding both electrical conversion and recirculation costs whilst exploiting firm demand for high-value heat. This strategy is relatively high-risk, high-reward, but options are identified for hedging these risks. We also identify and discuss new avenues for competition in this domain, which would not exist if fusion supplies electricity only.

We briefly review the recent developments in magnetohydrodynamics, which in particular deal with the evolution of magnetic fields in turbulent plasmas. We especially emphasize (i) the necessity of renormalizing equations of motion in turbulence where velocity and magnetic fields become H\"older singular; (ii) the breakdown of Laplacian determinism (spontaneous stochasticity) for turbulent magnetic fields; and (iii) the possibility of eliminating the notion of magnetic field lines, using instead magnetic path lines as trajectories of Alfvenic wave-packets. These methodologies are then exemplified with their application to the problem of magnetic reconnection -- rapid change in magnetic field pattern that accelerates plasma -- a ubiquitous phenomenon in astrophysics and laboratory plasmas. The necessity of smoothing out rough velocity and magnetic fields on a finite scale L implies that magnetohydrodynamic equations should be regarded as effective field theories with running parameters depending upon the scale L.

FLASH is a widely available radiation magnetohydrodynamics code used for astrophysics, laboratory plasma science, high energy density physics, and inertial confinement fusion. Increasing interest in magnetically driven inertial confinement fusion (ICF), including Pacific Fusion's development of a 60 MA Demonstration System designed to achieve facility gain, motivates the improvement and validation of FLASH for modeling magnetically driven ICF concepts, such as MagLIF, at ignition scale. Here we present a collection of six validation benchmarks from experiments at the Z Pulsed Power Facility and theoretical and simulation studies of scaling MagLIF to high currents. The benchmarks range in complexity from focused experiments of linear hydrodynamic instabilities to fully integrated MagLIF fusion experiments. With the latest addition of physics capabilities, FLASH now obtains good agreement with the experimental data, theoretical results, and leading ICF target design simulation code results across all six benchmarks. These results establish confidence in FLASH as a useful tool for designing magnetically driven ICF targets on facilities like Z and Pacific Fusion's upcoming Demonstration System.

High-yield inertial fusion offers a transformative path to affordable clean firm power and advanced defense capabilities. Recent milestones at large facilities, particularly the National Ignition Facility (NIF), have demonstrated the feasibility of ignition but highlight the need for approaches that can deliver large amounts of energy to fusion targets at much higher efficiency and lower cost. We propose that pulser-driven inertial fusion energy (IFE), which uses high-current pulsed-power technology to compress targets to thermonuclear conditions, can achieve this goal. In this paper, we detail the physics basis for pulser IFE, focusing on magnetized liner inertial fusion (MagLIF), where cylindrical metal liners compress DT fuel under strong magnetic fields and pre-heat. We discuss how the low implosion velocities, direct-drive efficiency, and scalable pulser architecture can achieve ignition-level conditions at low capital cost. Our multi-dimensional simulations, benchmarked against experiments at the Z facility, show that scaling from 20 MA to 50-60 MA of current enables net facility gain. We then introduce our Demonstration System (DS), a pulsed-power driver designed to deliver more than 60 MA and store approximately 80 MJ of energy. The DS is designed to achieve a 1000x increase in effective performance compared to the NIF, delivering approximately 100x greater facility-level energy gain -- and importantly, achieving net facility gain, or Qf>1 -- at just 1/10 the capital cost. We also examine the engineering requirements for repetitive operation, target fabrication, and chamber maintenance, highlighting a practical roadmap to commercial power plants.

We derive analytic dispersion relations for cold, orbitally constrained systems governed by the Vlasov equation. For magnetized plasmas, we obtain the first explicit relation for two-dimensional anisotropic BGK modes with finite magnetic field, showing that only a finite number of angular modes can become unstable and identifying a magnetic-field threshold for stabilization. In the gravitational case, we establish a bound on the growth rate of core perturbations, set by the potential's curvature. These results clarify how orbital constraints shape the spectrum and growth of kinetic instabilities in cold, collisionless media.