Bulk Plasma Research Articles

Kinetic simulations of capacitively coupled plasmas driven by tailored voltage waveforms with multi-frequency matching

Impedance matching is crucial for optimizing plasma generation and reducing power reflection in capacitively coupled plasmas (CCP). Designing these matchings is challenging due to the varying and typically unknown impedance of the plasma, especially in the presence of multiple driving frequencies. Here, a computational design method for impedance matching networks (IMNs) for CCPs is proposed and applied to discharges driven by tailored voltage waveforms (TVW). This method is based on a self-consistent combination of particle in cell/Monte Carlo collision simulations of the plasma with Kirchhoff’s equations to describe the external electrical circuit. Two Foster second-form networks with the same structure are used to constitute an L-type matching network, and the matching capability is optimized by iteratively updating the values of variable capacitors inside the IMN. The results show that the plasma density and the power absorbed by the plasma continuously increase in the frame of this iterative process of adjusting the matching parameters until an excellent impedance matching capability is finally achieved. Impedance matching is found to affect the DC self-bias voltage, whose absolute value is maximized when the best matching is achieved. Additionally, a change in the quality of the impedance matching is found to cause an electron heating mode transition. Poor impedance matching results in a heating mode where electron power absorption in the plasma bulk by drift electric fields plays an important role, while good matching results in the classical α-mode operation, where electron power absorption by ambipolar electric fields at the sheath edges dominates. The method proposed in this work is expected to be of great significance in promoting TVW plasma sources from theory to industrial application, since it allows designing the required complex multi-frequency IMNs.

Radial and Vertical Structures of Plasma Disk in Jupiter's Middle Magnetosphere

AbstractThe Juno mission flew through the plasma disk near the equator in Jupiter's magnetosphere frequently. We identify 274 plasma disk crossings of Juno between 10 and 40 RJ from PJ5 to PJ44. Using a forward modeling method that combines the JADE‐I time‐of‐flight and SPECIES data sets, we perform a survey of ion properties in the plasma disk. Ions are heated from 1.5 to 6 keV between 15 and 30 RJ. Density and temperature are locally anti‐correlated. Assumed to be related to centrifugal instabilities, cold, dense plasma are commonly observed near midnight. Plasma corotates around Jupiter and the rigid corotation breaks down outside 15–20 RJ. The plasma bulk velocity increases from the post‐dusk sector to the pre‐dawn sector featuring injection flows in the pre‐dawn sector, which is consistent with the Vasyliunas cycle. Strong outflows (>100 km/s) are commonly observed outside 20 RJ and the average radial velocity increases with radial distance. The ion abundance changes between 10 and 18 RJ and that might indicate plasma sources and/or sinks near Europa and Ganymede. The vertical distribution of ions is controlled by the balance between centrifugal, pressure gradient, and ambipolar electric field forces. An example near the M‐shell of 13.5 shows that average plasma temperature increases by a factor of 10 from the disk center to edge, because cold ions are more confined near the equator. Lighter ions with higher charge states have more mobility along the field line and have larger scale heights. The observations are compared with multi‐species diffusive equilibrium model.

Plasma properties in the vicinity of the last closed flux surface in hydrogen and helium fusion plasma discharges

The origin of the bi-Maxwellian electron energy distribution function (EEDF) observed in the scrape-off layer (SOL) of tokamak plasmas by means of Langmuir probes is still under discussion. It has been assumed that the ionization of hydrogen and deuterium neutrals by thermal electrons penetrating the SOL from the bulk plasma is the main reason for the appearance of a second Maxwellian. To validate this assumption, radial measurements of the electron temperatures and densities, or the plasma properties in helium plasmas in the GOLEM tokamak and the TJ-II stellarator were performed. The radial profiles of the low-temperature electron group densities follow the trend of the calculated radial profiles of the electron sources arising from the ionization of neutrals in both deuterium and helium plasmas in TJ-II. The difference in the radial location where the bi-Maxwellian EEDF appears can be explained by the difference in the rate coefficients for ionization of deuterium and helium. The results of probe measurements in GOLEM and the WEST tokamak divertor, at one radial location in the SOL, are compatible with the hypothesis concerning the ionization of neutral atoms and the type of the EEDF.

Proper retention factor for Avanafil and Dapoxetine Hydrochloride Separation and Validation with RP-HPLC in Bulk, Dosage Form and Rat Plasma; Validation, Degradation, and Pharmacokinetic Application

Proper retention factor for Avanafil and Dapoxetine Hydrochloride Separation and Validation with RP-HPLC in Bulk, Dosage Form and Rat Plasma; Validation, Degradation, and Pharmacokinetic Application

Void closure in a pulsed complex plasma in microgravity

A new experimental method for creating void-free complex (dusty) plasmas under microgravity conditions is presented. The method is based on a pulsed operation mode of a four-channel radio frequency generator for plasma sustainment. A dust cloud of micrometer-sized particles can be immersed in the bulk of a low temperature plasma under microgravity conditions. It typically contains a central volume depleted of particles—the void—that prevents the generation of large, continuous clouds. Experiments performed at different neutral gas pressures and discharge volumes during the microgravity phase of a parabolic flight show that the central void is closed completely once the pulsed operation mode is applied. The particle cloud shape and the density distribution within the cloud are practically independent of the pulse period within the investigated parameter range and mainly depend on the overall discharge parameters neutral gas pressure and discharge volume. This indicates that the pulsed operation of the plasma source does not introduce new physical effects on the particles aside from the void closure. The proposed method has great potential for future application in experimental facilities dedicated to fundamental studies of large three-dimensional, homogeneous complex plasma systems in microgravity.

Well-controlled mucosal exudation of plasma proteins in airways with intact and regenerating epithelium.

Superficial, systemic microcirculations, distinct from the pulmonary circulation, supply the mucosae of human nasal and conducting airways. Non-injurious, inflammatory challenges of the airway mucosa cause extravasation without overt mucosal oedema. Instead, likely reflecting minimal increases in basolateral hydrostatic pressure, circulating proteins/peptides of all sizes are transmitted paracellularly across the juxtaposed epithelial barrier. Thus, small volumes of extravasated, unfiltered bulk plasma appear on the mucosal surface at nasal and bronchial sites of challenge. Importantly, the plasma-exuding mucosa maintains barrier integrity against penetrability of inhaled molecules. Thus, one-way epithelial penetrability, strict localization, and well-controlled magnitude and duration are basic characteristics of the plasma exudation response in human intact airways. Invivo experiments in human-like airways demonstrate that local plasma exudation is also induced by non-sanguineous removal of epithelium over an intact basement membrane. This humoral response results in a protective, repair-promoting barrier kept together by a fibrin-fibronectin net. Plasma exudation stops once the provisional barrier is substituted by a new cellular cover consisting of speedily migrating repair cells, which may emanate from all types of epithelial cells bordering the denuded patch. Exuded plasma on the surface of human airways reflects physiological microvascular-epithelial cooperation in first line mucosal defense at sites of intact and regenerating epithelium.

Parameters of dust particle chains levitated vertically in a gas discharge plasma

The structural parameters of the dust particle chains suspended vertically in the electric field of a gas discharge are studied here. The investigations are performed by a numerical multi-block model based on a mean field approximation. The model describes the movement of ions and dust particles under the action of an external electric field, the electric field (Coulomb) of each charged dust particle, and the field of bulk plasma charge (ions and electrons) that screens the charges of dust particles. The gravity and the ion drag forces acting on the dust particles are also taken into account. Self-consistent chain parameters are calculated and compared for different number (1, 3, 5, and 7) and two different diameters (3.4 and 8.94 mcm) of dust particles in the chains. It is shown that an “ion wake” is formed behind the dust structures, and it grows with the number of dust particles and their size. With an increase in the number of dust particles in the chain, the charge of the first (top) dust particle increases. In the case of small dust particles, with an addition of new particles into the chain, the centers of the chains almost remain at the position of a single particle in the vertical direction. The chains of big dust particles move as a whole down in the direction of gravity with an increase in the number of particles.

Enhancing simultaneous determination of some angiotensin II receptor antagonists and amlodipine in plasma using HPTLC with fluorescence densitometry: Independent fluorescence detection of the co-administrative drugs in the mixture across various pH conditions

A novel and highly sensitive high-performance thin-layer chromatographic (HPTLC) method was developed and validated to quantify a combination of five pharmaceutical mixtures spiked to human plasma. The compounds comprised Amlodipine (AML) along with five angiotensin II receptor antagonist drugs (AIIRAs), namely Olmesartan (OLM), Telmisartan (TLM), Candesartan (CAN), Losartan (LOS), and Irbesartan (IRB). HPTLC was performed on silica gel 60 F254 plates using a mobile phase of Toluene: ethyl acetate: methanol: acetone: acetic acid (6:1.5:1:0.5:1, v/v/v/v/v). In a pioneering move, a reflectance/fluorescence detection mode was employed to identify two concurrently administered drugs at different pH levels for the first time. This method utilized the same chromatographic system, incorporating a specific measurement for AML at a neutral medium to achieve its maximum fluorescence at a 360 nm excitation wavelength, and measuring emission using a 540 nm optical filter. The process involved obtaining a very low fluorescence response from AIIRA. Subsequently, to enhance AIIRA’s fluorescence, the plate was sprayed with perchloric acid to transition to a strong acidic medium, ultimately attaining the maximum fluorescence of AIIRA using various excitation wavelengths and a 400 nm emission filter. Through this strategic process, we could optimize the fluorescence signals of both drugs, thereby elevating the sensitivity of detection for this drug combination. AML demonstrated a linear range of 18–300 ng/band, while AIIRAs drugs exhibited a linear range of 6–150 ng/band. The method satisfied the International Conference on Harmonization (ICH) criteria for recovery, precision, repeatability, and robustness, showcasing exceptional sensitivity. The approach was successfully applied to quantify AML and AIIRAs drugs in both bulk drug and plasma samples, achieving high recovery percentages and minimal standard deviations.

Control of spontaneous charging of sliding water drops by plasma-surface treatment

Slide electrification is the spontaneous separation of electric charges at the rear of water drops sliding over solid surfaces. This study delves into how surfaces treated with a low-pressure plasma impact water slide electrification. Ar, O2, and N2 plasma treatment reduced the drop charge and contact angles on glass, quartz, and SU-8 coated with 1H,1H,2H,2H-perfluoroctyltrichlorosilane (PFOTS), and polystyrene. Conversely, 64% higher drop charge was achieved using electrode-facing treatment in plasma chamber. Based on the zeta potential, Kelvin potential, and XPS measurements, the plasma effects were attributed to alterations of the topmost layer's chemistry, such as oxidation and etching, and superficially charge deposition. The surface top layer charges were less negative after electrode-facing and more negative after bulk plasma treatment. As a result, the zeta potential was less negative after electrode-facing and more negative after bulk plasma treatment. Although the fluorinated layer was applied after plasma activation, we observed a discernible impact of plasma-glass treatment on drop charging. Plasma surface modification offers a means to adjust drop charges: electrode-facing treatment of the fluorinated layer leads to an enhanced drop charge, while plasma treatment on the substrate prior to fluorination diminishes drop charges, all without affecting contact angles or surface roughness.

Time-helicity de-resonation (T-H) diagram for energy-selective mixing of charged particles during sawtooth crashes in tokamaks

It was recently shown that there exists a narrow parameter window where benign sawtooth crashes cause only mixing of bulk plasma and slowed-down alpha particle ‘ash’, while leaving MeV-class fast alphas largely unperturbed (Bierwage et al 2022 Nat. Commun. 13 3941). Here, we revisit the underlying physical picture and reframe it in a manner that may be suitable for systematic analyses of this phenomenon in modeling, simulation and experimental studies. In particular, we propose a graph that we call ‘time-helicity de-resonation diagram’ (short: T-H diagram) that captures the physical essence of energy-selectivity of sawtooth-particle interactions and visualizes it in a compact, intuitive way. Moreover, the regimes of good confinement and strong mixing during a sawtooth crash can be discerned via a single figure of merit: the T-H radius. The concept is introduced here on the basis of simulation results and would eventually benefit from further validation when applied to suitable empirical data.

Dual-frequency sheath oscillations and consequences on the ion and electron transport in dielectric barrier discharges at atmospheric pressure

Current–voltage characteristics, space- and time-resolved optical emission spectroscopy, and 1D fluid modeling are used to examine the effect of dual-frequency sheath oscillations on the ion and electron transport in dielectric barrier discharges sustained by a combination of low frequency (LF, 50 kHz, 650 V) and radiofrequency (RF, 5.3 MHz, 195 V) voltages, exhibiting the α-to-γ mode transition. On one hand, when polarities of the LF and RF voltages are opposite, an electric field near the LF cathode (due to LF cathode sheath) drives the secondary electrons to the plasma bulk and an opposite electric field between the sheath edge and the LF anode attracts the electrons toward the LF cathode (to maintain quasi-neutrality in the plasma bulk). At the sheath edge, electrons become trapped and ions drift toward the cathode and the anode simultaneously according to their position in the gap. On the other hand, when the RF voltage has the same polarity as the LF voltage, the total applied voltage increases and this yields to enhanced production of electrons and ions in the sheath. To maintain quasi-neutrality in the bulk, the electric field along the gap exhibits the same polarity as the one in the sheath, allowing electrons created in the sheath to be evacuated toward the LF anode. The behavior of the LF cathode is, therefore, controlled by the LF sheath, and, thus, by the LF voltage amplitude, while the behavior in the bulk and at the anode alternates on the time scale of the RF voltage.

Characterization of hydrogen plasma bulk in 2.45 GHz ECR ion source

Electron cyclotron resonance ion source (ECRIS) as the front-end for the neutral beam injection (NBI) system has production potential of hydrogen ion beam in tokamaks. In this article, two methods are hired to evaluate the plasma ion fractions: H+, H2+ and H3+ with respect to the ECRIS operating parameters. The methods are matured by utilizing 28 important volume and surface collisions processes. We encrypted a global 0D numerical model (GM) based on solving the nonlinear equations of plasma species as the first method. Next, a 3D plasma simulation (PS) is applied through the finite element method (FEM) by COMSOL Multiphysics. Then, the results are compared to the experimental data.

An incoherent Thomson scattering system for measurements near plasma boundaries.

Laser Thomson scattering (LTS) is a minimally invasive measurement technique used for determining electron properties in plasma systems. Sheath model closure validation requires minimally invasive measurements of the electron properties that traverse the boundaries between the bulk plasma, the presheath, and the plasma sheath. Several studies have probed the radial properties along the surface of discharge electrodes with laser-based diagnostics and electrostatic probes. These measurements provide valuable insight into the electron properties in this dynamic region. However, sheath model calibration requires plasma property measurements perpendicular to plasma bounding surfaces, in this case, along the electrode normal vector between discharge electrodes. This work presents the development of a discharge plasma cell and laser Thomson scattering system with a measurement volume step of 1mm normal to plasma bounding surfaces. The laser Thomson scattering measurements are made between a set of discharge electrodes separated by ∼25mm that are used to generate a pulsed argon plasma. The spatial distribution of electron temperature and density is measured at several discharge voltages between 8 and 20kV at a pressure of 8 Torr-Ar. It is determined that the system is statistically stationary and resembles a classic DC discharge plasma. The results are some of the first laser diagnostic-based "between electrode" measurements made along the plasma bounding electrode normal vector. A one-dimensional sheath model is applied to determine the near cathode electron properties, and it is determined that the edge of the presheath is probed in the high-voltage cases. As the lengths of the presheath and sheath decrease with decreasing voltage, the region recedes below the closest probed point to the cathode. To improve the performance of the diagnostic, the step size of the interrogation volume should decrease by an order of magnitude from 1mm to less than 100 μm, and the data acquisition strategy should be revised to increase the signal-to-noise ratio.

Plasma density enhancement in radio-frequency hollow electrode discharge

The plasma density enhancement outside hollow electrodes in capacitively coupled radio-frequency (RF) discharges is investigated by a two-dimensional (2D) particle-in-cell/Monte-Carlo collision (PIC/MCC) model. Results show that plasma exists inside the cavity when the sheath inside the hollow electrode hole is fully collapsed, which is an essential condition for the plasma density enhancement outside hollow electrodes. In addition, the existence of the electron density peak at the orifice is generated via the hollow cathode effect (HCE), which plays an important role in the density enhancement. It is also found that the radial width of bulk plasma at the orifice affects the magnitude of the density enhancement, and narrow radial plasma bulk width at the orifice is not beneficial to obtain high-density plasma outside hollow electrodes. Higher electron density at the orifice, combined with larger radial plasma bulk width at the orifice, causes higher electron density outside hollow electrodes. The results also imply that the HCE strength inside the cavity cannot be determined by the magnitude of the electron density outside hollow electrodes.

Numerical study of the plasma meniscus shape and beam optics in RF negative ion sources

In order to extract intense ion beams with good beam optics from hydrogen negative ion sources, it is important to control the shape of the plasma meniscus (i.e. beam emission surface). Recently, it is suggested experimentally that the degradation of beam optics in the RF negative ion sources may be due to the fluctuation of the distance d eff between the meniscus and the extraction grid caused by the fluctuation of the plasma density np . The purpose of this study is to make clear the dependence of d eff on np in the presence of a large amount of surface produced H- ions in order to understand such fluctuation of beam optics in RF sources For the purpose, 3D electrostatic PIC simulation was conducted taking the bulk plasma density as a parameter, investigating the extraction region of a H- ion source. A large amount of the surface H- production on the PG has been taken into account under the assumption that the H- production rate is proportional to the bulk plasma density. The result shows that the effective distance d eff is proportional to np -1/2 even for a large amount of surface H- production. This dependence suggests that the bulk plasma density np is the key parameters to control d eff and the resultant beam optics extracted from the negative ion source.

Effects of low-frequency voltage on nonlinear standing wave excitation, plasma uniformity, and ion dynamics in dual-frequency asymmetric capacitive discharges

It is known that in very-high-frequency (VHF) capacitively coupled plasmas, the higher harmonics generated by nonlinear sheath motion can enhance the standing wave effect (SWE), which can lead to center-peaked plasma density profiles. In this work, an improved nonlinear electromagnetic model incorporating a transmission line model, an electron momentum balance model, a bulk plasma model, a collisionless nonlinear numerical sheath model, and an ion Monte-Carlo collision (MCC) model is developed to study the effects of low-frequency (LF) voltage VL on the nonlinear standing wave excitation, plasma uniformity, and ion energy and angular distribution functions (IEDFs and IADFs) in dual-frequency (DF) asymmetric capacitive argon discharges at relatively low pressure of 3 Pa. The plasma diffusion in the radial direction and ion dynamics within the LF oscillating sheath are self-consistently considered. The LF voltage VL at 2 MHz varies from 0 to 700 V while the HF voltage VH at 60 MHz is fixed at 100 V. Simulation results indicate that without the addition of an LF source (i.e. VL= 0 V), there are a considerable number of high-order harmonics with short wavelengths, leading to significant SWE and central peak in the radial plasma density profile. Nevertheless, the high-order harmonic excitations tend to be weakened and merely occur around the phase of the full LF sheath collapse due to a shorter characteristic damping time of the surface waves as VL increases. This, combined with increased surface wavelengths of both the driving frequency and the higher harmonics at a higher VL , leads to suppressed standing waves and improved plasma uniformity. Meanwhile, the simulations show that both the low and the high energy peaks of IEDF move towards higher energies, and the energy peak separation width ΔE becomes wider with the increase of VL . The IEDF at the radial center of the powered electrode exhibits a broader ΔE than that at the edge. For the IADF, an increased VL results in more ions incident on the electrode with a smaller deflection angle. Because of a thinner sheath and a higher sheath voltage at the electrode center, the peak value of IADF at the electrode center is greater than that at the edge.

Potentiometric detection of apomorphine in human plasma using a 3D printed sensor

Apomorphine is a dopamine agonist that is used for the management of Parkinson's disease and has been proven to effectively decrease the off-time duration, where the symptoms recur, in Parkinson's disease patients. This paper describes the design and fabrication of the first potentiometric sensor for the determination of apomorphine in bulk and human plasma samples. The fabrication protocol involves stereolithographic 3D printing, which is a unique tool for the rapid fabrication of low-cost sensors. The solid-contact apomorphine ion-selective electrode combines a carbon-mesh/thermoplastic composite as the ion-to-electron transducer and a 3D printed ion-selective membrane, doped with the ionophore calix[6]arene. The sensor selectively measures apomorphine in the presence of other biologically present cations – sodium, potassium, magnesium, and calcium – as well as the commonly prescribed Parkinson's pharmaceutical, levodopa (L-Dopa). The sensor demonstrated a linear, Nernstian response, with a slope of 58.8 mV/decade over the range of 5.0 mM–9.8 μM, which covers the biologically (and pharmaceutically) relevant ranges, with a limit of detection of 2.51 μM. Moreover, the apomorphine sensor exhibited good stability (minimal drift of just 188 μV/hour over 10 h) and a shelf-life of almost 4 weeks. Experiments performed in the presence of albumin, the main plasma protein to which apomorphine binds, demonstrate that the sensor responds selectively to free-apomorphine (i.e., not bound or complexed forms). The utility of the sensor was confirmed through the successful determination of apomorphine in spiked human plasma samples.

Multi‐Spacecraft Magnetic Field Reconstructions: A Cross‐Scale Comparison of Methods

AbstractSpace plasma studies frequently use in situ magnetic field measurements taken from many spacecraft simultaneously. A useful data product of these measurements is the reconstructed magnetic field in a volume near the spacecraft observatory. We compare a standard Linear method of computing the magnetic field at arbitrary spatial points to two novel approaches: a Radial Basis Function interpolation and a time‐dependent 2D inverse distance weighted interpolation scheme called Timesync. These three methods, which only require in situ measurements of the magnetic fields and bulk plasma velocities at a sparse set of spatial points, are implemented on synthetic data drawn from a time‐evolving numerical simulation of plasma turbulence. We compare both the topology of the reconstructed field to the ground truth of the simulation and the statistics of the fluctuations found in the reconstructed field to those from the simulated turbulence. We conclude that the Radial Basis Function and Timesync methods outperform the Linear method in both topological and statistical comparisons.

Commissioning and first results of the 174 GHz collective Thomson scattering diagnostic at Wendelstein 7-X

Collective Thomson Scattering (CTS) diagnostics measure the scattering spectrum of monochromatic incident radiation off collective fluctuations in the plasma. In this contribution, we present the first results from the upgraded CTS diagnostic at Wendelstein 7-X (W7-X) operating in the frequency range between 172 and 176 GHz. This frequency range allows for minimization of noise originating from the electron cyclotron emission in the plasma. Consequently, the good signal-to-noise ratio allows measurements of fast ions or bulk plasma parameters with higher temporal resolution compared with the previously used 140 GHz system.

Particle injection methods in 3D-PIC MCC simulations applied to plasma grid biasing

In negative ion sources for the ITER Neutral Beam Injection system, the co-extraction of electrons is one of the main limiting factors. The current of co-extracted electrons can be decreased by applying a positive bias voltage to the Plasma Grid (PG) with respect to its source walls. Simulations using three-dimensional Particle-in-Cell Monte Carlo Collision (3D-PIC MCC) model are a powerful tool for studying the extraction region of such ion sources. However, the inclusion of both PG and source walls in the simulation domain is difficult due to numerical constraints. This study uses the 3D-PIC MCC code ONIX to explore the effects of particle injection models on plasma characteristics, using a flux injection model to regulate particle influx for a flat transition in potential from the bulk plasma to the simulation domain. Biasing of the PG above floating potential is possible using the flux injection scheme and results in a notable reduction in co-extracted electrons, corroborating with established experimental observations.