Background Plasma Research Articles

Sensitivity calculation for Monte Carlo particle simulations of neutrals in the plasma edge

AbstractAlgorithmic differentiation (AD) is a promising tool to accurately and efficiently compute sensitivities of several code outputs w.r.t. the code input parameters for the complex Monte Carlo code EIRENE. We illustrate the procedure by assessing the sensitivity of the total atom content w.r.t. the atomic and molecular cross‐sections and rate coefficients. AD does not suffer from decorrelation between primal and perturbed particle trajectories, which easily happens for finite difference (FD) sensitivity calculations. Consequently, the resulting statistical errors on AD sensitivities are up to five orders of magnitude smaller than the statistical errors on FD sensitivities for a JET low‐recycling background plasma. Although this statistical error reduction is promising for the application of the AD sensitivities in a future optimization or uncertainty quantification framework, we show that the AD sensitivities can blow up for a minority of long‐lived particles, especially occurring in detachment. Finally, we conclude that the Monte Carlo simulation and estimator type possibly need to be adapted to improve the accuracy of the sensitivities.

Nucleosome profiling at AR binding sites (ARBS) to provide a prognostic tool for metastatic castration-resistant prostate cancer (mCRPC).

42 Background: Prostate cancer is the second leading cause of cancer death for men in the U.S., and metastatic castration-resistant prostate cancer (mCRPC) is the most advanced form. In previous work, we have reported prognostic significance of PTEN-PI3K-AKT pathway and DNA damage response aberrations detected in cell-free DNA (cfDNA) of mCRPC. Here we utilized nucleosome profiling pattern at androgen receptor (AR) binding sites to further investigate AR signaling pathway dysregulation of mCRPC. Methods: Using the next-generation sequencing PredicineScore liquid biopsy assay (low-pass whole genome sequencing, LP-WGS), we profiled cfDNA fragment coverage centered around AR binding sites (ARBS). For 118 mCRPC patient samples, we quantified sample-level ARBS nucleosome profiling abnormality score (ARBS score) by comparing its centric fragment coverage with normal plasma background using standard Z-score. ARBS score was correlated with clinical outcomes including PSA response rate (PSA RR), progression-free survival and overall survival (OS). Results: By comparing ARBS nucleosome profiling and normal plasma background, we found that mCRPC plasma samples had significant ARBS nucleosome depletion, reflective of transcription factor binding and persistent AR signaling. Notably, ARBS nucleosome depletion was prostate cancer specific, and did not occur in other examined cancer types. ARBS score was highly correlated with tumor fraction (R=0.82), indicating that high tumor fraction mCRPC sample had higher AR binding. Importantly, we found that high ARBS score (>2) in pre-treatment samples (n=86) was associated with significantly shorter OS and lower PSA RR compared with low ARBS (≤2) group for patients receiving AR pathway inhibitor (ARPI) or taxane chemotherapy (13.4 vs. 19.4 months, p=7.3×10-3, 45.7% vs. 77.5%, p=5.2×10-3). By integrating genome-wide DNA methylation profiling, we further found that ARBS was hypo-methylated in mCRPC patients, and identified a strong negative correlation between DNA methylation beta value and ARBS score (R=-0.83). Conclusions: ARBS nucleosome profiling abnormality inferred from LP-WGS can serve as a key prognostic tool for mCRPC patients. Patients with high ARBS score had worse clinical outcomes, highlighting the potential clinical utility of cfDNA fragmentomics analysis for mCRPC patient stratification. Integration of mutation, copy number variation and fragmentomics profiling of ARBS may further enhance our understanding of molecular aberrations and treatment outcomes in mCRPC.

Coupling multi-fluid dynamics equipped with Landau closures to the particle-in-cell method

The particle-in-cell (PIC) method is successfully used to study magnetized plasmas. However, this requires large computational costs and limits simulations to short physical run times and often to set-ups of less than three spatial dimensions. Traditionally, this is circumvented either via hybrid-PIC methods (adopting massless electrons) or via magneto-hydrodynamic-PIC methods (modelling the background plasma as a single charge-neutral magneto-hydrodynamical fluid). Because both methods preclude modelling important plasma-kinetic effects, we introduce a new fluid-PIC code that couples a fully explicit and charge-conserving multi-fluid solver to the PIC code SHARP through a current-coupling scheme and solve the full set of Maxwell's equations. This avoids simplifications typically adopted for Ohm's law and enables us to fully resolve the electron temporal and spatial scales while retaining the versatility of initializing any number of ion, electron or neutral species with arbitrary velocity distributions. The fluid solver includes closures emulating Landau damping so that we can account for this important kinetic process in our fluid species. Our fluid-PIC code is second-order accurate in space and time. The code is successfully validated against several test problems, including the stability and accuracy of shocks and the dispersion relation and damping rates of waves in unmagnetized and magnetized plasmas. It also matches growth rates and saturation levels of the gyro-scale and intermediate-scale instabilities driven by drifting charged particles in magnetized thermal background plasmas in comparison with linear theory and PIC simulations. This new fluid-SHARP code is specially designed for studying high-energy cosmic rays interacting with thermal plasmas over macroscopic time scales.

Better Understanding of Hydrogen Pellet Ablation Cloud Spectra through the Occupation Probability Formalism in LHD

We have recently incorporated the occupation probability formalism (OPF) in the simulation model [C. Stehlé and S. Jacquemot, Astron. Astrophys. 271, 348 (1993)] to have a smooth transition from discrete lines to continuum spectrum in the wavelength range near the Balmer series limit. We have analyzed spectra measured for the hydrogen pellet ablation cloud in the Large Helical Device with the revised model, and have found that the electron density in the ablation cloud has a close correlation with the electron temperature of the background plasma. This type of correlation is first confirmed in the present analysis and should give a new insight in the simulation studies of pellet ablation for the magnetically confined fusion plasma.

Particle-in-cell simulations of electron–positron cyclotron maser forming pulsar radio zebras

Context. The microwave radio dynamic spectra of the Crab pulsar interpulse contain fine structures represented via narrowband quasiharmonic stripes. The pattern significantly constrains any potential emission mechanism. Similar to the zebra patterns observed, for example, in type IV solar radio bursts or decameter and kilometer Jupiter radio emission, the double plasma resonance (DPR) effect of the cyclotron maser instability may allow for interpretion of observations of pulsar radio zebras. Aims. We provide insight at kinetic microscales of the zebra structures in pulsar radio emissions originating close to or beyond the light cylinder. Methods. We present electromagnetic relativistic particle-in-cell (PIC) simulations of the electron–positron cyclotron maser for cyclotron frequency smaller than the plasma frequency. In four distinct simulation cycles, we focused on the effects of varying the plasma parameters on the instability growth rate and saturation energy. The physical parameters were the ratio between the plasma and cyclotron frequency, the density ratio of the “hot” loss-cone to the “cold” background plasma, the loss-cone characteristic velocity, and comparison with electron–proton plasma. Results. In contrast to the results obtained from electron–proton plasma simulations (for example, in solar system plasmas), we find that the pulsar electron–positron maser instability does not generate distinguishable X and Z modes. On the contrary, a singular electromagnetic XZ mode was generated in all studied configurations close to or above the plasma frequency. The highest instability growth rates were obtained for the simulations with integer plasma-to-cyclotron frequency ratios. The instability is most efficient for plasma with characteristic loss-cone velocity in the range vth = 0.2 − 0.3c. For low density ratios, the highest peak of the XZ mode is at double the frequency of the highest peak of the Bernstein modes, indicating that the radio emission is produced by a coalescence of two Bernstein modes with the same frequency and opposite wave numbers. Our estimate of the radiative flux generated from the simulation is up to ∼30 mJy from an area of 100 km2 for an observer at 1 kpc distance without the inclusion of relativistic beaming effects, which may account for multiple orders of magnitude.

Series of Small-scale Low Plasma β Magnetic Flux Ropes Originating from the Same Longitudinal Region: Parker Solar Probe Observations

In this study, we report on small-scale magnetic flux ropes (SMFRs) observed as a compact series in a narrow Carrington longitudinal range during three Parker Solar Probe (PSP) encounters. First, during ∼1.5 days of PSP's inbound part of Encounter 4, we identified a series of 11 SMFRs within 1.°4 in longitude over the radial distance of ∼8.4 R ☉ (from ∼44 to 35 R ☉). The identified SMFRs lasted from ∼0.5 to 1.8 hr, and adjacent events were separated mostly by a few hours and up to ∼6.5 hr at the longest, but some events were very closely spaced with intervals of a few ∼tens of minutes or less apart. Most of the identified SMFRs are successfully fitted to the force-free model. The SMFRs are clearly distinguished from the surroundings by a notable reduction in plasma β, which itself was comparably low (less than unity) in the background plasma. Furthermore, the magnetic field and plasma flow within the SMFRs fluctuated significantly less than the more turbulent surroundings. The fluctuations in the surrounding medium exhibited occasional Br polarity reversal (possibly switchbacks) and were Alfvénic to a large extent with far weaker compressional components. The majority of these key features with some differences have also been found in the series of SMFRs and their surroundings identified within 1.°3 or less in longitude during Encounters 1 and 5. We speculate that these SMFRs were repetitively generated by successive reconnection within a very narrow angular zone located close to the Sun but not necessarily at the same radial position.

Kinetic simulations comparing quasi-parallel and quasi-perpendicular piston-driven collisionless shock dynamics in magnetized laboratory plasmas

Magnetized collisionless shocks are common in astrophysical systems, and scaled versions can be created in laboratory experiments by utilizing laser-driven piston plasmas to create these shocks in a magnetized background plasma. A key parameter for these experiments is the angle θB between the shock propagation direction and the background magnetic field. We performed quasi-1D piston-driven shock simulations to explore shock formation, evolution, and key observables relevant to laboratory experiments for a range of shock angles between θB=90° to θB=30°. Our results show that the spatial and temporal scales of shock formation for all angles considered are similar when expressed in terms of the perpendicular component of the magnetic field. In a steady state, ion and electron temperatures become more isotropic, and the electron-to-ion temperature ratio is higher for smaller θB. At θB=30°, ion heating parallel to the magnetic field becomes dominant, associated with more ions being reflected at one discontinuity and subsequently trapped by the next discontinuity due to shock reformation.

Exploration on the Possible Bump-on-Tail Instability in VASIMR

The bump-on-tail (BOT) instability is generally caused by a beam of energetic particles existing in relatively cold background plasma. The employment of second-stage wave-driven module in the variable specific impulse magnetoplasma rocket (VASIMR) yields the production of energetic ions, which could drive the BOT instability. The present work explores this possibility for the first time via numerical simulations based on the experimental data on the VASIMR, i.e., referring to VX-50. It is found that the BOT instability does exist even in the plume region away from antenna. The results indicate that velocity space diffusion provides a stabilizing effect on the nonlinear evolution of waves, while dissipation in the bulk plasma essentially impedes it. To show the practical values implied by these computations, the influences of this BOT instability on the power coupling and thrust are investigated particularly. These findings are valuable for VASIMR, as well as other plasma thrusters that yield energetic particles inside relatively cold background plasma, to suppress BOT instability and thus increase the power coupling efficiency and thrust performance.

Fast Spacecraft Charging Mitigation With Plasma Contactors During High‐Power Electron Beam Emission in the GEO Environment

AbstractThe emission of artificial electron beams causes an active spacecraft charging effect that poses challenges for electron beam sounding experiments. Plasma contactors have been suggested as a potential solution to this issue based on theoretical and simulation studies, but the impact of these contactors on spacecraft potentials in complex space environments remains unexplored. This paper investigates the fast active charging effect of spacecraft under the action of plasma contactors in the geostationary Earth orbit environment by using a two‐dimensional Particle‐in‐Cell model. We analyze and compare the effects of the “ion emission” process of the plasma contactor and the “electron collection” process of the background plasma on the spacecraft potential at different beam and environment parameters, as well as the coupling effects when both processes occur simultaneously. Our findings indicate that the “ion emission” and “electron collection” processes are the primary factors in determining the spacecraft potential in relatively thin and dense plasma environments, respectively. When the two processes occur together, the contactor plasma cloud increases the “electron collection” return current by expanding the range of the positive potential area around the spacecraft, while the dense background plasma mediates and provides the return current to ensure that the plasma contactor can effectively suppress the active charging effect even after the ion sheath is present. These results suggest that a reasonable adjustment of the contactor parameters based on the relationship between the electron beam and the background environment can effectively extend the normal emission time of the electron beam.

Non-Gaussian Signal Statistics’ Impact on LIBS Analysis

A detailed study has been carried out to reveal signal statistics’ impact on analysis sensitivity in laser-induced breakdown spectroscopy (LIBS) measurements. For several signals measured simultaneously, it was demonstrated that space-, spectra- and time-integrated plasma emission followed a normal distribution while the spectra- and time-resolved LIBS signal (atomic line intensity, plasma background emissions) distribution functions were biased compared to a Gaussian distribution function. For the first time in LIBS, the impact of a non-Gaussian distribution function on the limit of detection (LOD)’s determination has been studied in detail for single-shot spectra as well as for averaged spectra. Here, we demonstrated that the non-symmetrical distribution of the LIBS signals influenced the estimated LODs, so knowledge of a LIBS signal’s distribution function provides more reliable results, and the analysis sensitivity can be wrongly estimated if Gaussian distribution is presumed.

Observations and simulations of large-scale traveling ionospheric disturbances during the January 14-15, 2022 geomagnetic storm

Using the total electron content (TEC) observations from GPS, and simulations from the Thermosphere Ionosphere Electrodynamic General Circulation Model (TIEGCM), this work investigates the large-scale traveling ionospheric disturbances (LSTIDs) and the possible involved drivers during the geomagnetic storm on January 14-15, 2022. Based on the term analysis of O+ continuity equation in TIEGCM, it is found that the traveling atmospheric disturbances in equatorward winds are responsible for the LSTIDs, with minor contributions from plasma drifts owing to the prompt penetration electric field. A strong interhemispheric asymmetry of the LSTIDs is observed, which might be attributed to both the equatorward wind disturbances and background plasma. The stronger wind (plasma) disturbances occurs in the winter hemisphere than that in the summer hemisphere. The maximum magnitude of LSTIDs in electron density disturbances occurs at ∼250 and ∼270 km in the northern and southern hemispheres, respectively, owing to both the thermospheric equatorward winds and background plasma. An interesting phenomenon that tail-like LSTIDs occur at the dip equator and low latitudes might be related to the eruption of the Tonga volcano, but it is not well reproduced in TIEGCM that deserves further exploration in a future study.

Unintended gas breakdowns in narrow gaps of advanced plasma sources for semiconductor fabrication industry

Occurrence of unintended gas breakdown in the narrow gaps of plasma processing chambers is one of the critical challenges in developing advanced plasma sources. We present a combined experimental and theoretical study of unintended discharges in the narrow gaps of plasma processing chambers and report significant drop of the gas breakdown voltage in the presence of a background plasma facing the gap. Experimentally measured breakdown voltages decrease in subsequent breakdown events due to wall erosion caused by the discharge. Therefore, preventing and mitigating the first discharge is of paramount importance. An analysis of kinetic simulation results indicates that the charged particle influx from the background plasma in the processing chamber into the gap is responsible for the onset of early breakdown: higher charged particle density within the gap modifies the electric field profile, allowing unintended breakdowns to occur at a significantly reduced threshold voltage.

Deciphering the physical basis of the intermediate-scale instability

We study the underlying physics of cosmic ray (CR)-driven instabilities that play a crucial role for CR transport across a wide range of scales, from interstellar to galaxy cluster environments. By examining the linear dispersion relation of CR-driven instabilities in a magnetised electron–ion background plasma, we establish that both the intermediate and gyroscale instabilities have a resonant origin, and show that these resonances can be understood via a simple graphical interpretation. These instabilities destabilise wave modes parallel to the large-scale background magnetic field at significantly distinct scales and with very different phase speeds. Furthermore, we show that approximating the electron–ion background plasma with either magnetohydrodynamics (MHD) or Hall-MHD fails to capture the fastest-growing instability in the linear regime, namely the intermediate-scale instability. This finding highlights the importance of accurately characterising the background plasma for resolving the most unstable wave modes. Finally, we discuss the implications of the different phase speeds of unstable modes on particle–wave scattering. Further work is needed to investigate the relative importance of these two instabilities in the nonlinear, saturated regime and to develop a physical understanding of the effective CR transport coefficients in large-scale CR hydrodynamics theories.

Preanalytical impact on the accuracy of measurements of glucagon, GLP-1 and GIP in clinical trials

Background Plasma concentrations of glucagon, GLP-1 and GIP are reported in numerous clinical trials as outcome measures but preanalytical guidelines are lacking. We addressed the impact of commonly used blood containers in metabolic research on measurements of glucagon, GLP-1 and GIP in humans. Methods Seventeen overweight individuals were subjected to an overnight fast followed by an intravenous infusion of amino acids to stimulate hormonal secretion. Blood was sampled into five containers: EDTA-coated tubes supplemented with DMSO (control), a neprilysin inhibitor, aprotinin (a kallikrein inhibitor) or a DPP-4 inhibitor, and P800 tubes. Plasma was kept on ice before and after centrifugation and stored at −80 Celsius until batch analysis using validated sandwich ELISAs or radioimmunoassays (RIA). Results Measures of fasting plasma glucagon did not depend on sampling containers, whether measured by ELISA or RIA. Amino acid-induced hyperglucagonemia was numerically higher when blood was collected into P800 tubes or tubes with aprotinin. The use of p800 tubes resulted in higher concentrations of GLP-1 by RIA compared to control tubes but not for measurements with sandwich ELISA. Plasma concentrations of GIP measured by ELISA were higher in control tubes and negatively affected by P800 and the addition of aprotinin. Conclusions The choice of blood containers impacts on measurements of plasma concentrations of glucagon, GLP-1 and GIP, and based on this study, we recommend using EDTA-coated tubes without protease inhibitors or P800 tubes for measurements of glucagon, GLP-1 and GIP in clinical trials.

Thermal quench induced by a composite pellet-produced plasmoid

Injecting shattered pellets is the critical concept of the envisaged ITER disruption mitigation system (DMS). Rapid deposition of large amounts of material should presumably result in controlled cooling of the entire plasma. A considerable transfer of thermal energy from the electrons of the background plasma to the ions accompanies a localized material injection due to the ambipolar expansion along the magnetic field line of the cold and dense plasmoid produced by the ablated pellet. Radiation initially plays the dominant role in the energy balance of a composite plasmoid containing high-Z impurities. A competition between the ambipolar expansion and the radiative losses defines the Thermal Quench scenario, including the amount of pre-quench thermal energy radiated on a short collisional timescale—possibly detrimental for the plasma-facing components. The present work quantifies plasmoid energy balance for disruption mitigation parameters. For pure hydrogen injection, up to 90% of the pre-pellet electron thermal energy may go to the newly injected ions. We also demonstrate that a moderate high-Z impurity content within the plasmoid can reduce highly localized radiation at the beginning of the expansion. The thermal energy will then dissipate on the much longer ion collisional timescale, which would be attractive for ITER DMS.

Dynamics of a Plasma Cloud Generated by a Compact Coaxial Gun upon Expansion into Vacuum and Large-Volume Background Plasma in an External Magnetic Field

Dynamics of a Plasma Cloud Generated by a Compact Coaxial Gun upon Expansion into Vacuum and Large-Volume Background Plasma in an External Magnetic Field

Dynamics of a Plasma Cloud Generated by a Compact Coaxial Gun upon Expansion into Vacuum and Large-Volume Background Plasma in an External Magnetic Field

Results of experiments on injection of dense plasma clouds created by a small-scale coaxial generator into vacuum and large-volume background plasma in an ambient magnetic field are presented. The regime of an “infinite” background medium that allows studying the plasma-cloud dynamics on the scale of about one meter in the directions perpendicular and parallel to a quasi-uniform magnetic field is realized on “Krot” plasma device. The dynamics of the diamagnetic cavity appearing upon magnetic-field expulsion by a plasma blob, the electromagnetic noise appearing in the cavity, along with the evolution of plasma-cloud structure during injection and at the stage of its decay, were studied. It is demonstrated that the key properties of the cloud dynamics that are typical of the active space and high-energy laboratory experiments, including complete expulsion of the magnetic field from the cloud and development of the flute instability at its boundary, are reproduced at low injection speed (below 30 km/s) and low plasma energy (on the order of 0.1 J).

Simulation study of the influence of E× B drift on tungsten impurity transport in the scrape-off layer

Tungsten (W) is used as the plasma-facing material in the divertor region of future fusion reactors, such as ITER; however, its concentration in the core plasma must be maintained at an extremely low level. W transport in the scrape-off layer (SOL), which is related to the source of core W contamination, has been extensively studied. In this study, the influence of E× B drift on the transport of W impurity in the SOL is studied via numerical simulations of a model case based on EAST upper single-null configuration with high recycling divertor plasma. W transport is simulated using DIVIMP on the background plasma obtained from scape-off layer plasma simulation-ITER simulation including drifts. The E× B drift of W ions is introduced based on the background electric field. Therefore, both the direct E× B drift effect of W ion and the indirect effect via background plasma on W transport in the SOL are studied. The influence on the flux of W impurities entering confined plasma across the last closed flux surface Γ enter is focused on, which is expected to be proportional to the core W concentration. Results reveal that Γ enter is mainly from the outer (inner) target under a favorable (unfavorable) toroidal field B T and can be increased by more than one order of magnitude compared with the case without drifts; this reflects the significant effect of E× B drift. The effects due to the background plasma and the poloidal and radial E× B drift of W ion, as well as the related mechanisms, are analyzed in detail for three stages of W transport in the SOL: effective sputtering from the target, leakage from the divertor, and entry into the confined plasma.

Association between atherogenic index of plasma and asymptomatic intracranial arterial stenosis in middle-aged and elderly women: A cross-sectional study in Shandong, China

Background and aimsThe atherogenic index of plasma (AIP) is associated with progression of atherosclerosis and used to describe how pro- or anti-atherogenic components are balanced. However, the association of AIP with asymptomatic intracranial arterial stenosis (aICAS) is uncertain. The purpose of this study is to investigate the association between AIP and aICAS in rural China. Methods and resultsA total of 1990 participants aged ≥40 years free of stroke or transient ischemic attack were enrolled in this study. The presence of aICAS was examined by Transcranial Doppler ultrasound and confirmed by magnetic resonance angiography. The adjusted AIP (aAIP) was calculated according to the ratio of TG and HDL-C and further separated into 4 quartiles. Multiple logistic regression was used to investigate the association between aAIP and aICAS, and the dose-response relationship was explored by restricted cubic spline. After adjusting for conventional confounders, aAIP was significantly higher in the aICAS group than that in the non-aICAS group. Furthermore, the common odds ratios for aICAS risk increased with increasing aAIP quartiles. Multivariate logistic regression revealed that aAIP was independently associated with aICAS in female or middle-aged and elderly (age ≥50 years), and superior to other lipid profiles. Multiple-adjusted spline regression showed the dose-response association between aAIP levels and aICAS prevalence. ConclusionsAIP might be independently and positively associated with the prevalence of aICAS in middle-aged and elderly women, which might be superior to traditional and nontraditional lipid profiles in rural China.

Observation of D2 molecule line emission after massive D2 injection into runaway electron plateaus in DIII-D

Molecular deuterium line emission is observed in both the visible and ultraviolet (UV) wavelength ranges after massive (> 100 Torr-L) injection of D2 gas into post-disruption runaway electron (RE) dominated plasmas in the DIII-D tokamak. D2 UV line emission is found to be the dominant source of radiated power, surpassing D Lyα. Interpretive modeling with a collisional-radiative model (CRM) indicates that D2 radiation surpasses D radiation because Lyα is strongly trapped, while D2 UV lines are mostly untrapped. The CRM also indicates that the D2 line emission is completely dominated by RE impact (rather than thermal electron impact), so the D2 line emission can serve as a good diagnostic for the spatial localization of REs. Analysis of D2 visible lines indicates that the D2 molecules in the plasma are thermally equilibrated with the background plasma, with vibrational, rotational, and kinetic temperatures all near 0.3 eV. D2 spectroscopy therefore serves as a convenient diagnostic of background plasma temperature. Measurement of D2 radiated power also appears to serve as a useful diagnostic for constraining neutral transport modeling.