Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ion number density”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

Atom probe tomography characterization of ion and neutron irradiated Alloy 800H

Alloy 800H is a high temperature and creep resistant alloy, and is considered as a candidate alloy for use in Generation IV nuclear reactor systems. To clarify the alloy's behavior under irradiation, samples were subjected to single ion irradiations up to 20 dpa at 440°C, a dual ion/He irradiation to 17 dpa at 460°C, and a neutron irradiation to 17 dpa at 385°C. The irradiated microstructures were characterized using atom probe tomography to complement previously published transmission electron microscopy studies. Additionally, after single ion irradiation, sparse fine Al and Ti clusters were observed after 1dpa, while high number densities of nanoscale Ni-Al-Ti clusters, Cr-Ti rich carbides, and Si-decorated dislocation loops developed after 10 and 20 dpa. The microstructure formed during single ion irradiation exhibited a strong depth dependence. No significant differences were observed between the single and dual ion irradiations. While Ni-Al-Ti clusters, Cr-Ti rich carbides, and Si-decorated dislocation loops were also observed after neutron irradiation, the neutron-irradiated microstructure differed from those found in the ion-irradiated samples.

36 MATERIALS SCIENCE↗

Nanoscale Ion Transport Enhances Conductivity in Solid Polymer-Ceramic Lithium Electrolytes

The predictive design of flexible and solvent-free polymer electrolytes for solid-state batteries requires an understanding of the fundamental principles governing the ion transport. In this work, we establish a correlation among the composite structures, polymer segmental dynamics, and lithium ion (Li + ) transport in a ceramic-polymer composite. Elucidating this structure–property relationship will allow tailoring of the Li + conductivity by optimizing the macroscopic electrochemical stability of the electrolyte. The ion dissociation from the slow polymer segmental dynamics was found to be enhanced by controlling the morphology and functionality of the polymer/ceramic interface. The chemical structure of the Li + salt in the composite electrolyte was correlated with the size of the ionic cluster domains, the conductivity mechanism, and the electrochemical stability of the electrolyte. Polyethylene oxide (PEO) filled with lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) or lithium bis(fluorosulfonyl) imide (LiFSI) salts was used as a matrix. A garnet electrolyte, aluminum substituted lithium lanthanum zirconium oxide (Al-LLZO) with a planar geometry, was used for the ceramic nanoparticle moieties. Further, the dynamics of the strongly bound and highly mobile Li + were investigated using dielectric relaxation spectroscopy. The incorporation of the Al-LLZO platelets increased the number density of more mobile Li + . The structure of the nanoscale ion-agglomeration was investigated by small-angle X-ray scattering, while molecular dynamics (MD) simulation studies were conducted to obtain the fundamental mechanism of the decorrelation of the Li + in the LiTFSI and LiFSI salts from the long PEO chain.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Constitutional under-potential plating (CUP) – New insights for predicting the morphological stability of deposited lithium anodes in lithium metal batteries

Lithium (Li) cycling in lithium metal batteries (LMBs) inevitably witnesses a non-planar Li/electrolyte interface (cellular/dendritic microstructure) affecting the performance, safety and reliability of LMBs. Herein, the “constitutional under-potential plating” (CUP) concept is presented for the first time to understand the evolution, formation and predict the morphological stability/instability of the Li/electrolyte interface. Here, the CUP condition, a measure of the driving force for interface perturbation, occurs when the true potential gradient at the interface is smaller than the equilibrium electrodeposition potential gradient at the Li/electrolyte interface resulting in the electrolyte ahead of the Li-electrolyte interface to be in an under-potential state. Detailed CUP criterion calculations conducted herein establishes the relationship of battery operating conditions (e.g. current density, potential gradient) with the electrolyte properties (e.g. Li-ion transference number and diffusion coefficient) predicting the morphological stability condition at the interface. In addition to the driving force for interface perturbation, the capillarity effect or the Gibbs Thomson effect, a measure of the resistance to perturbation of the perturbed interface is discussed for understanding the influence of interfacial energy/adhesion energy between the current collector and the deposited Li on the overall morphological stability of the Li/electrolyte interface, a critical criterion for the unwanted dendritic morphology.

25 ENERGY STORAGE↗

Spectroscopic Insight on Neodymium Solvation in Lithium Borohydride-Supported Electrolyte

Borohydride-based electrolytes have recently emerged as promising media for the electrodeposition of electropositive metals, including rare earth (RE) elements. While the presence of supporting alkali metal cations and RE counteranions provides essential electrochemical conductivity for achieving fast metal electrodeposition, interactions between the host ligand and solvated neodymium (Nd) complexes remain unclear. This study provides insights into the coordination structure of a concentrated and directly solvated Nd salt in a lithium borohydride-supported electrolyte. Our spectroscopic results indicate that the RE coordination environment is significantly influenced by the solvation mechanism, which can vary between metathesis and complexation pathways, primarily dictated by stoichiometric factors. Under dilute conditions, nearly complete metathesis of anions leads to a high coordination number for the host ligand (borohydride), consistent with the previously reported solvated Nd speciation in chlorine-free electrolytes. In contrast, concentrated dissolution of the Nd salt in the supported electrolyte is dominated by a complexation pathway featuring a Li-ion-paired complex with a low coordination number of the host ligand. Density functional theory (DFT) calculations indicated that the observed blue shift in the borohydride vibration was the result of an increase in electron density drawn into the terminal B–H interbond region from the hydride as the coordination changed from Li to Nd. In conjunction with DFT results, vibrational analyses allowed correlation of the experimental shifts associated with changes in Nd ligation and coordination spheres, further consolidating the prevalence of highly chloride-coordinated species under concentrated conditions. In conclusion, the outcomes of this work illuminate the distinctive and heterogeneous coordination structures that the electroactive RE species can adopt at high concentrations in lithium borohydride-supported electrolytes, as a key step to comprehend the reported metal electrodeposition performance in these media.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Comprehensive Characterization of Water Group Ion Composition and Distributions in Saturn's Magnetosphere With Cassini Plasma Spectrometer Data

Saturn's magnetosphere is continuously supplied with neutrals from the Enceladus plume and the icy rings, which undergo ionization and charge-exchange to form a complex water-group plasma environment. While the Cassini Plasma Spectrometer (CAPS) instrument has provided extensive compositional information, detailed separation of individual water-group ion species in time-of-flight (TOF) data has not previously been achieved. In this study, we perform forward modeling of CAPS-IMS energy-per-charge (E/Q) and TOF spectra obtained between 2004 and 2012 to resolve O + , OH + , H 2 O + , and H 3 O + and to characterize their plasma properties, including number density, temperature, and thermodynamic κ. Our results demonstrate that O + is the dominant thermal ion species throughout Saturn's magnetosphere, comprising up to ∼70% of the total ion population beyond ∼5 Saturn radii (R S ). In contrast, molecular ions such as OH + , H 2 O + , and H 3 O + dominate closer to Enceladus but rapidly dissociate into atomic ions between ∼5 and 10 R S . This radial region is also characterized by the steepest increase in plasma flow speed, which rises from ∼40% to ∼80% of rigid corotation. Simultaneously, ion velocity distributions approach Maxwell–Boltzmann equilibrium, as indicated by high kappa values. These findings provide new constraints on the ion–neutral chemistry that regulates the balance between molecular and atomic ions in Saturn's magnetosphere. They also emphasize the critical role of the 5–10 R S region as a transition zone for both plasma composition and dynamics. Our results refine previous CAPS-based studies and underscore the need to incorporate seasonal variability and ionospheric coupling into future global models of Saturn's plasma environment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Yttrium Contenting Compositionally Complex Medium-Entropy Li-Garnet Electrolyte with Improved Ionic Conductivity

The compositionally complex medium/high-entropy design concept can greatly expand the categories and affect the properties of the materials. With such a designing concept, a medium-entropy Li-garnet electrolyte with appropriate yttrium content (formula Li 6.6 La 3 ZrNb 0.3 Ta 0.3 Hf 0.3 Y 0.1 O 12 ) shows a record-high ionic conductivity of ∼5.7 × 10 –4 S/cm, the highest reported for any single-site substituted high/medium-entropy Li-garnet. The assembled Li metal symmetric cells also show stable long-term cycling (0.1 mA/cm 2 for over 200 h). Neutron powder diffraction and Rietveld refinement results indicate that a competing conduction mechanism between (1) occupancy on high mobility of 96h sites and (2) the associated site vacancies and the bond length requires an appropriate content of Y for enhanced ionic conductivity. Li-ion hopping through the bottleneck can also contribute to the conductivity. Finally, density functional theory and Born–Oppenheimer molecular dynamics simulations also indicate the high mobility and number of hopping transitions of Li ions, contributing to the high ionic conductivity.

Li-garnet↗

Cumulants: it’s more than you think

Cumulants of baryon number are given considerable attention in analyses of heavy-ion collision experiments as possible signatures of the QCD critical point. In this work, we show that the values of the lowest three cumulants can also be utilized to recover information about the isothermal speed of sound and its logarithmic derivative with respect to the baryon number density. This result provides a new method for obtaining information about fundamental properties of nuclear matter studied in heavy-ion collisions, with consequences for both the search for the QCD CP and neutron star studies. While the approximations and the model comparison we considered apply to experiments at low energies, the approach itself can be used at any collision energy provided that measurements of cumulants of baryon number distribution aswell as their temperature dependence are available.

Sorensen, Agnieszka↗

Dislocation loop evolution in Kr-irradiated ThO 2

Here, the early stage of microstructural evolution of ThO 2 , under krypton irradiation at 600, 800, and 1000°C, was investigated using in situ transmission electron microscopy (TEM). Dislocation loops grew faster, whereas their number density decreased with increasing irradiation temperature. Loop density was found to decrease with ion dose. Interstitial dislocation loops, including Frank loops with Burgers vector of $\textit{a}/3\langle111\rangle$ and perfect loops with Burgers vector of $\textit{a}/2\langle110\rangle$, were determined by traditional TEM and atomic resolution–scanning TEM techniques. Atomistic and mesoscale level modeling are performed to interpret experimental observations. The migration energy barriers of defects in ThO 2 were calculated using density-functional theory. The energetics of different dislocation loop types were studied using molecular dynamics simulations. Loop density and diameter were analyzed using a kinetic rate theory model that considers stoichiometric loop evolution. This analysis reveals that loop growth is governed by the mobility of cation interstitials, whereas loop nucleation is determined by the mobility of anion defects. Lastly, a rate theory model was used to extract the diffusion coefficients of thorium interstitials, oxygen interstitials, and vacancies.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Fluctuations and correlations of baryonic chiral partners

The exploration of critical phenomena in phase transitions of strongly interacting matter governed by quantum chromodynamics (QCD) is one of the goals of present ultrarelativistic heavy-ion collision experiments at BNL and CERN. The key research direction is to locate the putative critical point on the phase diagram of QCD linked to the chiral symmetry restoration at finite temperature and/or density. One of the main theoretical tools used for this purpose is the fluctuations of conserved charges, such as the net-baryon number. However, due to experimental limitations, analyses of heavy-ion collision data suffer from a very doubtful basing of the net-proton number being a proxy for the total net-baryon number fluctuations. In this work, we use the parity doublet model to investigate the fluctuations of the net-baryon number density in hot and dense hadronic matter. The model accounts for chiral criticality within the mean-field approximation. We focus on the qualitative properties and systematics of the first- and second-order susceptibility of the net-baryon number density, and their ratios for nucleons of positive and negative parity, as well as their correlator. We show that the fluctuations of the positive-parity nucleon do not necessarily reflect the fluctuations of the total net-baryon number density at the phase boundary of the chiral phase transition. We also investigate the nontrivial structure of the correlator. Furthermore, we discuss and quantify the differences between the fluctuations of the net-baryon number density in the vicinity of the chiral and liquid-gas phase transition in nuclear matter. We indicate a possible relevance of our results with the interpretation of the experimental data on net-proton number fluctuations in heavy-ion collisions. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Accelerating kinetic plasma simulations with machine-learning-generated initial conditions

Computational models of plasma technologies often solve for the system operating conditions by time-stepping an initial value problem to a quasi-steady solution. However, the strongly nonlinear and multi-timescale nature of plasma dynamics often necessitate millions, or even hundreds of millions, of steps to reach convergence, reducing the effectiveness of these simulations for computer-aided engineering. We consider acceleration of kinetic plasma simulations via data-driven machine-learning-generated initial conditions, which initialize the simulations close to their final quasi-steady-state, thereby reducing the number of steps to reach convergence. Three machine-learning models are developed to predict the density and ion kinetic profiles of capacitively coupled plasma discharges relevant to the microelectronics industry. The models are trained on kinetic simulations over a range of device operating frequencies and pressures. Best performance was observed when simulations were initialized with ion kinetic profiles generated by a convolutional neural network, reducing the mean number of steps to reach convergence by 17.1× when compared to initialization with a zero-dimensional global model. We also outline a workflow for continuous data-driven model improvement and simulation speedup, with the aim of generating sufficient data for full device digital twins.

Artificial neural networks↗

Higher order cumulants of net baryon-number distributions at non-zero μ

Using recent results on higher order cumulants of conserved charge fluctuations from lattice QCD, we construct mean, variance, skewness, kurtosis, hyper-skewness and hyper-kurtosis of net-baryon number distributions for small baryon chemical potentials μ B . For the strangeness neutral case ( μ S = 0) at fixed ratio of electric charge to baryon number density ( n Q n B = 0.4 ), which is appropriate for a comparison with heavy ion collisions, we present results for κ B σ 2 B , S B σ 3 B / M B , κ H B σ 4 B and S H B σ 5 B / M B on the crossover line for the chiral transition, T pc ( μ B ). Furthermore, continuum extrapolations for this pseudo-critical transition line have recently been reported by HotQCD up to baryon chemical potentials μ B ≃ 300 MeV [A. Bazavov, et al., Chiral crossover in QCD at zero and non-zero chemical potentials, Phys. Lett. B 795 (2019) 15–21. arXiv:1812.08235 .]. These cumulant ratios are of direct relevance for comparisons with corresponding ratios measured by STAR in the BES-I and II runs at beam energies s NN ≥ 20 GeV . In particular, we point out that recent high statistics results on skewness and kurtosis of net-baryon number distributions obtained by STAR at s NN = 54.4 GeV put strong constraints on freeze-out parameters and are consistent with predictions from thermal QCD.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Low-Energy Ion Scattering Intensities from Supported Nanoparticles: The Spherical Cap Model

Supported nanoparticles are of great importance to many technologies like fuel processing and chemical synthesis using catalysts and electrocatalysts, energy storage and generation using fuel cells and batteries, electrochemistry, magnetic information storage, and more. Low-energy ion scattering spectroscopy (LEIS) with noble gas ions like He+ is a powerful tool for the characterization of nanoparticles dispersed across flat support surfaces due to its ability to probe the elemental composition in the topmost atomic layer of a surface, providing quantitative information regarding the size and number density of nanoparticles. In this work, we present a derivation of the LEIS intensities expected from nanoparticles and the support material as a function of the average particle size, their number per unit area, and their contact angle with the support when modeled as spherical caps of the nanoparticle material dispersed over the surface of a flat support. The model assumes that the ion intensities are determined only by the physical blocking of linear ion trajectories and independent of the tilt angle of the local surface relative to the incident and scattered ion directions, an assumption we support by quantitative modeling of published data which tested tilt-angle effects. The model is a generalization to arbitrary contact angles of the hemispherical cap model, which assumes 90° contact angle and has been widely used to model spectroscopic signals in LEIS (and also in Auger and photoelectron spectroscopies) during nanoparticle growth. Here, this new model quantitatively reveals how LEIS signals are sensitive not only to the diameter and number density of the nanoparticle but also to their contact angle (or height/diameter ratio). With the use of additional data (e.g., from microscopy or adsorption microcalorimetry), the model presented here will enable more accurate determination of the average size, shape, and number density of supported nanoparticles based on LEIS intensity measurements.

99 GENERAL AND MISCELLANEOUS↗

An experimental characterization of core turbulence regimes in Wendelstein 7-X

First results from the optimized helias Wendelstein 7-X stellarator (W7-X) have shown that core transport is no longer mostly neoclassical, as is the case in previous kinds of stellarators. Instead, power balance analysis has shown that turbulent transport poses a serious limitation to the global performance of the machine. Furthermore, several studies have found this particularly relevant for ion transport, with core ion temperatures becoming clamped at relatively low values of T i ≃ 1.7 keV, except in the few scenarios in which turbulence can be suppressed. In order to understand the precise turbulent mechanisms at play and thus design improved performance scenarios, it is important to have a clear understanding of the parametric dependencies of turbulent fluctuations, and the relation between them and turbulent transport. As a first step in this direction, in this work we use Doppler reflectometry measurements carried out during a number of relevant operational scenarios to provide a systematic characterization of ion-scale ( k ⊥ ρ i ≃ 1) density fluctuations in the core of W7-X. Then, we study the relation between fluctuation amplitude and plasma profiles and show how distinct regimes can be defined for the former, depending on normalized gradients and . Furthermore, we discuss the importance of other potentially relevant parameters such as T e / T i , E r or collisionality. Comparing the different regimes, we find that turbulence amplitude depends generally on the gradient ratio η i = , as would be expected for ITG modes, with the exception of a range of discharges, for which turbulence suppression may be better explained by an ITG to TEM transition triggered by a drop in collisionality. Finally, we show a number of scenarios under which T i,core > 1.7 keV is achieved and how core fluctuations are suppressed in all of them, thus providing experimental evidence of microturbulence being the main responsible for the limited ion confinement in W7-X.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Solid-State, Single-Ion Conducting, Polymer Blend Electrolytes with Enhanced Li + Conduction, Electrochemical Stability, and Limiting Current Density

The development of solid-state polymer electrolytes with high lithium conductivity is crucial to improve lithium-ion battery performance and ameliorate the safety challenges associated with current solvent-based electrolytes. Unfortunately, sluggish polymer segmental dynamics are known to constrain conductivity enhancements in solid-state polymer electrolyte systems, limiting overall performance. In this work, a glassy single-ion-conducting polymer, poly[lithium sulfonyl(trifluoromethane sulfonyl)imide methacrylate] (PLiMTFSI), was blended with a flexible polymer, poly(oligo-oxyethylene methyl ether methacrylate) (POEM), and the impact of PLiMTFSI molecular weight and ion concentration on the thermal and ion-conducting behavior of blend electrolytes was investigated. High ionic conductivities approaching 1 × 10 -2 S/cm at 150 °C were realized in this polymer blend electrolyte system as a result of decoupling Li + transport from polymer segmental dynamics. The decoupled ion transport was attributed to the packing frustration of the glassy PLiMTFSI – sufficient percolating free volume was generated to produce effective ion diffusion pathways. This decoupling was tunable as the ion transport could be altered from being closely coupled to the polymer segmental dynamics (Vogel–Tammann–Fulcher-like) to hopping (Arrhenius-like) by increasing the PLiMTFSI molecular weight and ion concentration. Moreover, the immobilized TFSI anion resulted in high Li + selectivity (Li + transference number = 0.9), high electrochemical stability (up to 4.7 V against Li + / Li), and limiting current density of 1.8 mA/cm 2 (electrolyte thickness = 0.05 cm). These features suggest that this single-ion-conducting, polymer blend electrolyte might be a promising alternative to a benchmark system – salt-doped polyethylene oxide. Moreover, the above characteristics can support the battery operation at higher voltages using energy-dense Li metal anodes, with faster charging rates and enhanced energy/power densities. Altogether, the results suggest that polymer chain packing frustration can be exploited to overcome the constraints of slow polymer segmental relaxations to achieve rapid and highly selective ion transport and enhanced performance in solid-state polymer electrolytes.

36 MATERIALS SCIENCE↗

Single-Ion Conducting Polymer Electrolyte Enabled via Aza-Michael Addition

Single-ion conducting polymer electrolytes present a possible route for achieving high energy density next-generation batteries due to their flexible nature and high cation transport number. The trifluoromethanesulfonimide (TFSI) functional group represents one of the most efficient tethered anions for Li-ion conductivity. However, the covalent attachment of TFSI groups into polymer electrolytes has been challenging and costly due to its synthetic difficulty and limited commercial availability of building blocks. Here, we present a new polymer electrolyte synthesized by the Michael addition reaction between poly(allylamine) (PAA) and a vinyl TFSI anion under mild reaction conditions. The resulting PAATFSI exhibits a lower glass transition temperature and several orders of magnitude higher Li-ion conductivity than similar TFSI-based single-ion conducting homopolymers in the dry state. Moreover, PAATFSI exhibits excellent Li-ion conductivity (2.7 × 10 –4 S/cm at 30 °C) with the addition of a plasticizer (60 wt % of carbonate solvent), which enables stable lithium–metal battery performance. In conclusion, this study demonstrates a facile route for the synthesis of new single-ion conducting polymer electrolytes that opens a library of possibilities to enable the realization of polymer-based batteries.

25 ENERGY STORAGE↗

Behavior of helium cavities in ion-irradiated W-Ni-Fe ductile-phase toughened tungsten

This study reports on the distribution of helium (He) cavities in a hot-rolled W-Ni-Fe ductile-phase toughened tungsten (DPT W) composite irradiated to a dose and a helium concentration that are comparable to those in the material after 5-year irradiation in a conceptual fusion power plant. Here, the DPT W sample consists of W particles embedded in a ductile-phase NiFeW matrix with a nominal composition of 90W-7Ni-3Fe by weight. It was hot-rolled to a thickness reduction by 87% (87R DPT W). Sequential irradiation of the material with 1.2 MeV Ni + ions to a fluence of 2.15×10 16 Ni + /cm 2 and 90 keV He + ions to 6.5×10 15 He + /cm 2 was performed at 973 K. Larger He cavities with a lower number density are observed in NiFeW than W. Helium cavities are aggregated preferentially along the NiFeW/W interphase boundary. This behavior is not observed along the W/W grain boundary under the same irradiation conditions. A similar distribution of He cavities with comparable sizes is observed in 87R DPT W irradiated with He + ions only, suggesting that the He cavity distribution is not significantly affected by pre-irradiation of Ni + ions at 973 K. The data from this study does not provide any evidence for formation of visible voids or Ni precipitates in W phase irradiated with Ni + ions only. Diffusion and clustering of He atoms and mobile vacancies during ion irradiation at 973 K are responsible for the formation and growth of the He cavities.

36 MATERIALS SCIENCE↗

Probing the high-density nuclear symmetry energy with the Ξ – / Ξ 0 ratio in heavy-ion collisions at s N N ≈ 3 GeV

Recent beam energy scan experiments at the BNL Relativistic Heavy Ion Collider by the STAR Collaboration found that hadronic interactions dominate the collective flow and the proton cumulant ratios are driven by baryon number conservation in a region of high baryon density in √ s NN = 3 GeV Au + Au reactions, indicating that the dense medium formed in such collisions is likely hadronic matter. Within an updated a relativistic transport model with momentum dependent isoscalar and isovector single-nucleon mean-field potentials corresponding to different symmetry energies at suprasaturation densities, the n/p, π – /π + , $K$ $^{0}_{s}$/K + , Σ – /Σ + , and Ξ – /Ξ 0 ratios are studied for central Au + Au collisions at √ s NN = 3 GeV, where the maximum central density reaches about (3.6–4.0) ρ 0 . The doubly strange Ξ – /Ξ 0 ratio is found to have the strongest sensitivity to the variation of high-density nuclear symmetry energy. Furthermore, the Ξ – /Ξ 0 ratio in relativistic heavy-ion reactions at √ s NN~3 GeV may help probe sensitively the poorly known symmetry energy of dense neutron-rich matter critically important for understanding various properties of neutron stars.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Wave steepening and shock formation in ultracold neutral plasmas

We present observations of wave steepening and signatures of shock formation during expansion of ultracold neutral plasmas formed with an initial density distribution that is centrally peaked and decays exponentially with distance. The plasma acceleration and velocity decrease at large distance from the plasma center, leading to central ions overtaking ions in the outer regions and the development of a steepening front that is narrow compared to the size of the plasma. The density and velocity change dramatically across the front, and significant heating of the ions is observed in the region of steepest gradients. For a reasonable estimate of electron temperature, the relative velocity of ions on either side of the front modestly exceeds the local sound speed (Mach number M≳1). This indicates that by sculpting steep density gradients, it is possible to create the conditions for shock formation, or very close to it, opening a new avenue of research for ultracold neutral plasmas.

Computational fluid dynamics↗