Ion Density Profiles Research Articles

Variability of Earth’s ionospheric outflow in response to the dynamic terrestrial exosphere

The most abundant neutral constituent in the exospheric region (i.e., beyond ≈ 500 km altitude) is the atomic hydrogen (H); however, its density distributions predicted by physics-based models have been challenged by satellite-based observations of its far ultraviolet emissions. This discrepancy may impact magnetospheric ions’ densities and velocities since numerous chemistry and ion-neutral coupling interactions rely sensitively on the underlying neutral hydrogen population. The Polar Wind Outflow Model a first-principled model for relevant ion species in the high-latitude ionosphere, is employed to investigate the role of neutral H on the ionospheric outflow. Specifically, variability in the outflow of ionospheric H+, He+, N+, and O+ as a response to systematic enhancement and depletion of H number densities were simulated. The altitude-dependent ion density and energy partition profiles vary with neutral H density, solar activities, and ion species. These findings suggest that the exosphere plays a crucial role in controlling the production and loss of ions through ionospheric chemistry, as well as the energy contributions by altering ion-neutral-electron collisions and the ambipolar electric field to the high-latitude ionospheric outflow. As a result, the escape rates of the ionospheric outflow are directly associated with exospheric distributions. This work potentially helps understand the dominant mechanisms of atmospheric escape, particularly during a hydrogen-rich early Earth’s and exoplanet’s atmosphere, which is known to play a significant role in understanding the evolution of Earth’s atmosphere.

Kinetic effects on the interaction of counter-propagating plasma shocks inside an ICF hohlraum

The interaction and interpenetration of two counter-propagating plasma shocks are investigated via hybrid fluid-PIC (particle-in-cell) simulations. This study seeks to probe the kinetic effects and ion collisions on the structure of colliding plasma shocks in complex multi-ion-species plasma, in particular, the presence of the expansion of high-Z plasma bubbles against the low-Z filled gas inside an ICF hohlraum. The superposition of shock wave results in a wave-like electric field in the downstream region. The electric field can further reduce the kinetic energy of the incoming particles, and modulate the ion density profile. It finally generates a new downstream platform of high temperature and high density. However, on the hundred-ps time scale, cumulative ion collisions can still significantly alter the structure of the shock wave and the reflection of ions by the shock front. This study will help to improve the predictions of hohlraum plasma states and the understanding of the shock wave interactions.

Magnetized multi-component plasmas sheath characteristics with three isothermal ion species

The dynamics of magnetized plasma sheath with three isothermal ions is studied numerically using the fluid model. The three ions ( Ar+,Kr+, and Xe+ ) are assumed to be singly charged and have the same temperature. In the beginning, the concentrations for Ar+ and Xe+ were taken in the same proportions, and then the Kr+ ion density is introduced slowly in the plasma. It is found that the presence of a third ion has changed the dynamics of ion densities inside the sheath and simultaneously their respective velocities. Ion density and electron density profiles were estimated inside the sheath region in the presence of a weak and strong magnetic field as well at different Kr+ ion concentrations. Further, the dependence of magnetic field obliqueness on the ion density profile and its velocities inside the sheath is also studied, and observed ion density bunching at different conditions. Based on the results, we have tried to introduce a concept of the derivation of system Bohm velocity by considering the ion mass as equivalent mass for the three interacting ion systems in the plasma. It is found that the system Bohm velocity with the equivalent mass effect is larger than their individual thermal velocities and close to Ar+ Bohm velocity, and higher than the system Bohm velocity with reference to the density ratio and their respective individual Bohm velocities in an isothermal multi-component plasma.

Core density profile control by energetic ion anisotropy in LHD

Electron and impurity ion density profiles have been controlled by using tangential and perpendicular neutral beams for plasma heating in a stellarator/heliotron for the first time. Reduced anisotropy of stored energy for energetic ion En⊥/Enǁ has resulted in an inward electron and impurity transport, forming a core electron density peaking. Increased anisotropy leads to a flat or hollowed electron density profile with an impurity exhaust in a core region [Yoshinuma et al., Nucl. Fusion 49, 062002 (2009)]. A high confinement state of particles in LHD has yet to be achieved, except for a temporal state of an electron density peaking created by a pellet injection. As a pioneering and crucial research result, the operation of energetic ion anisotropy by neutral beams has newly demonstrated that the direction of the radial transport of bulk and impurity ions can be controlled. At the same time, the overall plasma performance rises in neutron flux and stored energy. On the other hand, the increase in the anisotropy flattens the density profile. This new finding holds promise for a control knob of nuclear fusion reactors to enhance fusion power output.

One-dimensional, multi-fluid model of the plasma-wall transition. II. Negative ions

The plasma-wall transition is investigated by a one-dimensional steady-state multifluid model, which was presented in detail in Part I [T. Gyergyek et al., AIP Adv. 14, 045201 (2024)]. In this work, the plasma-wall transition is analyzed for the case where the plasma consists of singly charged positive ions, electrons, and singly charged negative ions. When the temperature and initial density of the negative ions are varied, a transition between two types of solutions of the model is observed. We call them the low and high solution, with respect to the absolute value of the potential drop. When the density and temperature of the negative ions are above a critical value, the low solution is observed. As the mass of the positive ions increases, these critical values also increase, but only until the ion mass is below about 1000 electron masses. With larger ion masses, the critical density of the negative ions and the temperature no longer change. In the low solution, the potential drop in front of the sheath is determined by the negative ions and is smaller in absolute terms than in the case of the high solution, where the potential drop in front of the sheath is determined by the electrons. If the problem is analyzed on the pre-sheath scale, the transition between the low and high solution is very sharp. However, when the neutrality condition is replaced by the Poisson equation, this transition becomes blurred and the solutions of the model equations exhibit oscillations. The role of the smallness parameter is highlighted. It is shown how the initial electric field is determined. Deviation of the negative ion density profile from the Boltzmann relation is discussed.

Impact of Hydrostatic Pressure on Molecular Structure and Dynamics of the Sodium and Chloride Ions in Portlandite Nanopores.

In order to address the issues of energy depletion, more resources are being searched for in the deep sea. Therefore, research into how the deep-sea environment affects cement-based materials for underwater infrastructure is required. This paper examines the impact of ocean depth (0, 500, 1000, and 1500 m) on the ion interaction processes in concrete nanopores using molecular dynamics simulations. At the portlandite interface, the local structural and kinetic characteristics of ions and water molecules are examined. The findings show that the portlandite surface hydrophilicity is unaffected by increasing depth. The density profile and coordination number of ions alter as depth increases, and the diffusion speed noticeably decreases. The main cause of the ions' reduced diffusion velocity is expected to be the low temperature. This work offers a thorough understanding of the cement hydration products' microstructure in deep sea, which may help explain why cement-based underwater infrastructure deteriorates over time.

Axisymmetric hybrid Vlasov equilibria with applications to tokamak plasmas

We derive axisymmetric equilibrium equations in the context of the hybrid Vlasov model with kinetic ions and massless fluid electrons, assuming isothermal electrons and deformed Maxwellian distribution functions for the kinetic ions. The equilibrium system comprises a Grad–Shafranov partial differential equation and an integral equation. These equations can be utilized to calculate the equilibrium magnetic field and ion distribution function, respectively, for given particle density or given ion and electron toroidal current density profiles. The resulting solutions describe states characterized by toroidal plasma rotation and toroidal electric current density. Additionally, due to the presence of fluid electrons, these equilibria also exhibit a poloidal current density component. This is in contrast to the fully kinetic Vlasov model, where axisymmetric Jeans equilibria can only accommodate toroidal currents and flows, given the absence of a third integral of the microscopic motion.

Characterization of N+ Abundances in the Terrestrial Polar Wind Using the Multiscale Atmosphere‐Geospace Environment

AbstractThe High‐latitude Ionosphere Dynamics for Research Applications (HIDRA) model is part of the Multiscale Atmosphere‐Geospace Environment model under development by the Center for Geospace Storms NASA DRIVE Science Center. This study employs HIDRA to simulate upflows of H+, He+, O+, and N+ ions, with a particular focus on the relative N+ concentrations, production and loss mechanisms, and thermal upflow drivers as functions of season, solar activity, and magnetospheric convection. The simulation results demonstrate that N+ densities typically exceed He+ densities, N+ densities are typically ∼10% O+ densities, and N+ concentrations at quiet‐time are approximately 50%–100% of N+ concentrations during storm‐time. Furthermore, the N+ and O+ upflow fluxes show similar trends with variations in magnetospheric driving. The inclusion of ion‐neutral chemical reactions involving metastable atoms is shown to have significant effects on N+ production rates. With this metastable chemistry included, the simulated ion density profiles compare favorably with satellite measurements from Atmosphere Explorer C and Orbiting Geophysical Observatory 6.

Molecular dynamics simulations of the effect of starch on transport of water and ions through graphene nanopores.

We use molecular dynamics simulations to unravel the molecular level mechanisms underlying the structure and dynamics of water and ions flowing through nanoporous starch-graphene membranes. Our findings indicate that there is a significant tendency for the formation of short-range order in close proximity to the graphene membrane surface. This leads to a greater concentration of water and ions, suggesting strong interactions between the membrane and the saltwater solution. Furthermore, we found that the starch-graphene membrane was most efficient in sieving out ions when the starch loading is 15 wt.%, and the pore diameter is 14Å. At these conditions, the starch-graphene membrane showed a high water transport rate and maintained a high level of ion rejection. We investigated the effect of loading of starch and the pore diameter on the pressure-induced transport, structure, and dynamics of Na+, Cl-, and water using the GROMACS 2021.4 package. We further analyze the density profiles of water and ions in the context of ion-polymer and water-polymer interactions and provide mechanistic insights into the piston-induced flow of saltwater through the starch-graphene membranes using Visual Molecular Dynamics (VMD) software.

Laboratory observation of electron energy distribution near three-dimensional magnetic nulls

The acceleration of electrons near three-dimensional (3D) magnetic nulls is crucial to the energy conversion mechanism in the 3D magnetic reconnection process. To explore electron acceleration in a 3D magnetic null topology, we constructed a pair of 3D magnetic nulls in the PKU Plasma Test (PPT) device and observed acceleration of electrons near magnetic nulls. This study measured the plasma floating potential and ion density profiles around the 3D magnetic null. The potential wells near nulls may be related to the energy variations of electrons, so we measured the electron distribution functions (EDFs) at different spatial positions. The axial variation of EDF shows that the electrons deviate from the Maxwell distribution near magnetic nulls. With scanning probes that can directionally measure and theoretically analyze based on curve fitting, the variations of EDFs are linked to the changes of plasma potential under 3D magnetic null topology. The kinetic energy of electrons accelerated by the electric field is 6 eV ( ) and the scale of the region where accelerating electrons exist is in the order of serval electron skin depths.

Variation of Surface Propensity of Halides with Droplet Size and Temperature: The Planar Interface Limit.

The radial number density profiles of halide and alkali ions in aqueous clusters with equimolar radius ≲1.4 nm, which correspond to ≲255 H2O molecules, have been extensively studied by computations. However, the surface abundance of Cl-, Br-, and I- relative to the bulk interior in these smaller clusters may not be representative of the larger systems. Indeed, here we show that the larger the cluster is, the lower the relative surface abundance of chaotropic halides is. In droplets with an equimolar radius of ≈2.45 nm, which corresponds to ≈2000 H2O molecules, the polarizable halides show a clear number density maximum in the droplet's bulk-like interior. A similar pattern is observed in simulations of the aqueous planar interface with halide salts at room temperature. At elevated temperature the surface propensity of Cl- decreases gradually, while that of I- is partially preserved. The change in the chaotropic halide location at higher temperatures than the room temperature may considerably affect photochemical reactivity in atmospheric aerosols, vapor-liquid nucleation and growth mechanisms, and salt crystallization via solvent evaporation. We argue that the commonly used approach of nullifying parameters in a force field in order to find the factors that determine the ion location does not provide transferable insight into other force fields.

Surface tension between liquids containing antagonistic and regular salts.

We investigate theoretically the surface tension between two immiscible liquids containing antagonistic and regular salts. The mean-field model includes the electrostatic energy and the Gibbs transfer energy of the ions. We find the Donnan potential difference between the liquids and solve the Poisson-Boltzmann equation to calculate the potential and ion density profiles. When the liquids contain only a regular salt, the surface tension increases when compared to the no-salt case; in contrast, the surface tension is reduced when they contain only an antagonistic salt. When both salts are present, the surface tension can either further decrease or increase depending on the preferential solvation of the ions, in agreement with published experimental observations.

Turbulence stabilization in tokamak plasmas with high population of fast ions

This letter provides a new physical insight into the fast ion effects on turbulence in plasmas having a high fast ion fraction and peaked fast ion density profile. We elucidate turbulence stabilization mechanisms by fast ions that result in internal transport barrier formation in the fast ion regulated enhancement mode plasma. Both linear and nonlinear gyrokinetic simulations show that the dominant turbulence suppression mechanisms are the dilution effects. In particular, we find that turbulence can be sufficiently suppressed solely by an inverted main ion density gradient due to fast ions, for the first time. New physical findings reported here improve our understanding of fast ion effects on turbulence, essential for fusion energy production where . Moreover, they will open up a new methodology to control plasma turbulence applicable to a wide range of plasma confinement regimes.

Are the oscillations found in magnetized collisional electronegative warm plasma artifact during plasma-surface interaction?

The formation of the sheath and its characteristics are investigated in an electronegative plasma under the impact of an oblique and constant magnetic field. Specifically, the role of the magnetic field and initial conditions of positive ion velocity on the sheath characteristics are uncovered by adopting the Sagdeev potential approach for the determination of the modified Bohm’s criterion. At the sheath edge, the upper and lower components of the positive ion velocity perpendicular to the probe/wall are determined and examined in greater detail. The sheath structure is found to be modified considerably by the initial conditions. If the magnetic field is strong near the sheath edge, then the profile of the positive ion density shows a pulse-like structure for the case when we neglect the x-component (parallel to the wall) of positive ion velocity at the sheath edge. This pulse-like structure disappears when all the components of the positive ion velocity at the sheath edge are included, i.e. the ions are considered to move in any direction but towards the probe. The occurrence or disappearance of pulse-like structures or oscillations completely depend upon the initial conditions. Since in real situations the ions are expected to move in any direction (carrying all their velocity components) and are not aligned only in the direction perpendicular to the wall, these oscillations are artifact in nature. The sheath thickness is observed to be reduced with increasing magnetic field strength and its inclination. In order to better understand the phenomenon, isothermal and adiabatic cases are also discussed.

Predictive simulations for plasma scenarios in the SMART tokamak

The SMall Aspect Ratio Tokamak (SMART) is a new spherical machine that is currently being constructed at the University of Seville (Mancini et al., 2021; Agredano-Torres et al., 2021). The operation of SMART will cover three different phases reaching an inductive plasma current (IP) of more than 500 kA, a toroidal magnetic field (BT) of 1 T and a pulse length of 500 ms (Mancini et al., 2021; Agredano-Torres et al., 2021). The main goal of the SMART tokamak is to study high plasma confinement regimes in a broad triangularity range (-0.5≤δ≤0.5) (Doyle et al., 2021; Doyle et al., 2021). While in phase 1 the ohmic heating alone is expected to provide enough power to access the H-mode, in phase 2 and phase 3 the access to the H-mode will be ensured by applying Neutral Beam Injection (NBI) as external heating system. The NBI will consist of one injector at 25 keV and 1 MW of power. The overall design of the NBI, including injection geometry, energy and power have been optimized using the ASCOT5 code (Hirvijoki et al., 2021). The SMART scenarios have been developed with the help of the free boundary equilibrium solver code FIESTA (Cunningham, 2013) coupled to the linear time independent, rigid plasma model RZIP (Lazarus et al., 1990) to calculate the target equilibria for all the different operational phases. To assess the feasibility of those scenarios, predictive modelling needs to be included to evaluate properly the evolution of the temperatures, density profiles for both electrons and ions. To this extent, the 1.5D transport code ASTRA (Pereverzev and Yushmanov, 2002) has been used including models for the ohmic current, bootstrap current and current driven by NBI. This contribution discusses the electron and ion density and temperature profiles obtained for various scenarios for phase 1 and 2 and presents the design study of the NBI.

Molecular insights into nuclear-magnetic-resonance properties of NaCl solution confined within calcite nanopores

Nuclear magnetic resonance (NMR) relaxation time is an important petrophysical property probing fluid dynamics within pores and providing valuable information to understand solid–fluid interactions and ion-related effects. In this study, we investigate NMR and diffusion properties of confined sodium chloride (NaCl) solutions within calcite nanopores via molecular dynamics (MD) simulations. Density profiles, self-diffusion coefficients, and NMR relaxation times were analyzed as a function of pore size, pore geometry, and NaCl concentration. MD results indicate that the presence of NaCl hinders water movement, resulting in reduced self-diffusion coefficients and relaxation times of water. Nevertheless, the influence of NaCl on water mobility is less than the effects caused by confinement and geometry. Strength of oscillations in water and ion density profiles increases when the slit pore size is reduced to 1 nm, suggesting that water-calcite interactions near pore surfaces are stronger in smaller pore sizes. Analysis of dipole–dipole correlation functions and probability distribution of correlation times shows that ions in aqueous solution have an impact on both intra- and intermolecular interactions. A comprehensive comparison of NMR relaxation times of water in both slit and spherical nanopores reveals that spherical models provide a more accurate representation of the confinement effect induced by pore networks within tight rocks. These findings provide a novel quantitative understanding of the impact of confinement, pore geometry, and ion concentration on the translational and rotational dynamics of confined fluids in calcite nanopores.

The Effects of the Upper Atmosphere and Corona on the Solar Wind Interaction With Venus

AbstractSince Venus has no substantial planetary magnetic field, the fast‐flowing solar wind plasma interacts directly with its ionosphere, upper atmosphere, and corona. The thermal atmosphere and hot oxygen corona of Venus are expected to play essential roles in the interaction process. To quantify their effects, we combine three major state‐of‐the‐art three‐dimensional numerical models, including (a) The Venus thermosphere general circulation model (VTGCM), (b) The adaptive mesh particle simulator (AMPS), and (c) A multi‐species Magnetohydrodynamics (MHD) model of Venus based on Block Adaptive Tree Solar‐wind Roe Upwind Scheme. The global multi‐species MHD model is used to simulate the planet's interaction with the ambient solar wind, where the background neutral densities are provided by VTGCM and AMPS. Using the above‐mentioned numerical models, the effect of the upper atmosphere and corona on solar wind driven ion escape rates are examined for three typical solar cycle conditions, including solar maximum, solar median, and solar minimum conditions. The model results suggest that even though the hot oxygen corona only has a minor effect on the ion density profiles in the ionosphere, it can enhance ion loss rate up to ∼25% for high solar activity.

Preliminary observation of tungsten-impurity suppression using on-axis ECRH by X-ray crystal spectroscopy in EAST

Impurities degrade tokamak plasma confinement by causing energy loss, diluting fuel concentration, and even terminating discharge in some extreme cases. Previously, the suppression effects of impurity accumulation due to on-axis electron cyclotron resonance heating (ECRH) have been studied on Experimental and Advanced Superconducting Tokamak (EAST) using extreme ultraviolet (EUV) spectroscopy. However, it is difficult to quantify changes in the tungsten (W) impurity profile since W-line emission in the EUV range cannot be easily resolved. X-ray crystal spectroscopy (XCS) is widely used to measure the ion temperature and rotation velocity of plasmas by using line emission in the soft X-ray range. In addition, the XCS can also be used to study the behavior of impurities. An in situ absolute intensity calibration of tangential XCS was conducted by analyzing calculations and measurements of bremsstrahlung radiation. After obtaining the calibration coefficient, the W44+-ion-density profiles were evaluated using Abel inversion operations and the spectral line of W XLV (W44+, 3.9095 Å). Thus, a direct observation of the W44+-impurity concentration suppressed by ECRH was accomplished. Such W44+-density profiles can be used in the future to analyze W transport in combination with impurity transport codes.

Inference of main ion particle transport coefficients with experimentally constrained neutral ionization during edge localized mode recovery on DIII-D

The plasma and neutral density dynamics after an edge localized mode are investigated and utilized to infer the plasma transport coefficients for the density pedestal. The Lyman-Alpha Measurement Apparatus (LLAMA) diagnostic provides sub-millisecond profile measurements of the ionization and neutral density and shows significant poloidal asymmetries in both. Exploiting the absolute calibration of the LLAMA diagnostic allows quantitative comparison to the electron and main ion density profiles determined by charge-exchange recombination, Thomson scattering and interferometry. Separation of diffusion and convection contributions to the density pedestal transport are investigated through flux gradient methods and time-dependent forward modeling with Bayesian inference by adaptation of the Aurora transport code and IMPRAD framework to main ion particle transport. Both methods suggest time-dependent transport coefficients and are consistent with an inward particle pinch on the order of 1 m s−1 and diffusion coefficient of 0.05 m2 s−1 in the steep density gradient region of the pedestal. While it is possible to recreate the experimentally observed phenomena with no pinch in the pedestal, low diffusion in the core and high outward convection in the near scrape-off layer are required without an inward pedestal pinch.

Decoupling of dust cloud and embedding plasma for high electron depletion in nanodusty plasmas

In recent years nanoparticles have become key technological products, e.g., as coatings with tunable optical gap in third generation solar cells, as nanocrystals for photonic applications, and as pharmaceutical nanocarriers. In particle sources, that use reactive, nanodusty plasmas, a high dust density changes the properties of the dusty plasma compared to a dust free plasma considerably, as electron depletion leads to a reduced number of free electrons. This is called the Havnes effect and was central for the understanding of the famous spokes in Saturns rings. We see here, that it is also important for technological applications. Using self excited dust density waves (DDW) as a diagnostic tool, we completely characterize an argon discharge with embedded amorphous hydrocarbon nanoparticles of different size and density. The results show, that electron depletion governs the charge of dust grains, while the size of the particles has only a weak influence. The ion density and electric potential profile are almost independent of both, dust size as well as dust density. This suggests, that the rf generated plasma and the dust cloud seem to coexist and coupling of both is weak.