We present trajectory simulation-based modeling to capture the interactions between ions and charged grains in dusty or complex plasmas. Our study is motivated by the need for a self-consistent and experimentally validated approach for accurately calculating the ion drag force and grain charge that determine grain collective behavior in plasmas. We implement Langevin dynamics in a computationally efficient predictor–corrector approach to capture multiscale ion and grain dynamics. Predictions of grain velocity, grain charge, and ion drag force are compared with prior measurements to assess our approach. The comparisons reveal excellent agreement to within ±20% between predicted and measured grain velocities [Yaroshenko et al., Phys. Plasmas 12, 093503 (2005) and Khrapak et al., Europhys. Lett. 97, 35001 (2012)] for 0.64, 1.25 μm grains at ∼20−500 Pa. Comparisons with the measured grain charge [Khrapak et al., Phys. Rev. E 72, 016406 (2005)] under similar conditions reveal agreement to within ∼20% as well. Measurements of the ion drag force [Hirt et al., Phys. Plasmas 11, 5690 (2004); IEEE Trans. Plasma Sci. 32, 582 (2004)] are used to assess the viability of the presented approach to calculate the ion drag force experienced by grains exposed to ion beams of well-defined energy. Excellent agreement between calculations and measurements is obtained for beam energies >10 eV, and the overprediction below 10 eV is attributed to the neglect of charge exchange collisions in our modeling. Along with critical assessments of our approach, suggestions for future experimental design to probe charging of and momentum transfer onto grains that capture the effect of space charge concentration and external fields are outlined.
more »
« less
Comment on “Ultracold plasma expansion in quadrupole magnetic field”
Bronin et al. [Phys. Rev. E 108, 045209 (2023)] recently reported molecular-dynamics simulations of ultracold neutral plasmas expanding in a quadrupole magnetic field. While the main results are in agreement with prior experimental measurements, we present data showing oscillations not captured in the simulations of Bronin et al. Plasmas formed using pulsed or continuous-wave ionization processes have similar confinement times.
more »
« less
- Award ID(s):
- 2009999
- PAR ID:
- 10588730
- Publisher / Repository:
- American Physical Society
- Date Published:
- Journal Name:
- Physical Review E
- Volume:
- 110
- Issue:
- 2
- ISSN:
- 2470-0045
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Previously, using an incompressible von Kármán–Howarth formalism, the behavior of cross-scale energy transfer in magnetic reconnection and turbulence was found to be essentially identical to each other, independent of an external magnetic (guide) field, in the inertial and energy-containing ranges [Adhikari et al., Phys. Plasmas 30, 082904 (2023)]. However, this description did not account for the energy transfer in the dissipation range for kinetic plasmas. In this Letter, we adopt a scale-filtering approach to investigate this previously unaccounted-for energy transfer channel in reconnection. Using kinetic particle-in-cell simulations of antiparallel and component reconnection, we show that the pressure–strain interaction becomes important at scales smaller than the ion inertial length, where the nonlinear energy transfer term drops off. Also, the presence of a guide field makes a significant difference in the morphology of the scale-filtered energy transfer. These results are consistent with kinetic turbulence simulations, suggesting that the pressure strain interaction is the dominant energy transfer channel between electron scales and ion scales.more » « less
-
We compare a variety of models used for the calculation of transport coefficients in dense plasmas, including average-atom models, models based on kinetic theory, structure matching effective potentials, and pair-potential molecular dynamics. In particular, we focus on the parameter space investigated in the second charged-particle transport coefficient code comparison workshop [Stanek et al., Phys. Plasmas 31, 052104 (2024)]. Each model is based on the self-consistent output of our average-atom calculations. Ionic transport properties are generated from implicit electron pair matched molecular dynamics simulations, bypassing the need for either dynamical electron simulations or on-the-fly electronic structure calculations. These matched pair potentials are generated in a nonlinear way using a classical mapping procedure, further avoiding an expensive force-matching procedure. We compare these results with the density functional theory data presented at the workshop, as well as a set of widely used parametric models, which we have modified to enhance accuracy, especially at the low- and high-temperature extremes of the parameter space. We also detail the non-trivial statistical aspect of converging ionic transport coefficients.more » « less
-
Ion wakes in two-dimensional complex plasma crystals play a crucial role in determining the structural and dynamical properties of these strongly coupled systems. Using the DRIAD molecular dynamics code [Matthews et al., Phys. Plasmas 27, 023703 (2020)], we investigate the formation and morphology of ion wakes beneath a hexagonal cluster of 19 microparticles under two limiting conditions: a collisionless ion flow and an electric-field-driven collisional flow. Our simulations reveal that, in both regimes, individual wakes form downstream of each microparticle, but their shape and intensity differ significantly. Collisionless flows produce elongated wakes with moderate ion-density enhancement, whereas collisional flows generate shorter, more concentrated wakes due to ion-neutral collisions. A striking feature observed in both cases is the emergence of a collective ion-density enhancement beneath the entire cluster, accompanied by a measurable increase in electric potential (up to ∼ 0.5–0.6 V).more » « less
-
This paper presents a summary of the element test simulations (calibration simulations) submitted by 11 numerical simulation (prediction) teams that participated in the LEAP-2017 prediction exercise. A significant number of monotonic and cyclic triaxial (Vasko, An investigation into the behavior of Ottawa sand through monotonic and cyclic shear tests. Masters Thesis, The George Washington University, 2015; Vasko et al., LEAP-GWU-2015 Laboratory Tests. DesignSafe-CI, Dataset, 2018; El Ghoraiby et al., LEAP 2017: Soil characterization and element tests for Ottawa F65 sand. The George Washington University, Washington, DC, 2017; El Ghoraiby et al., LEAP-2017 GWU Laboratory Tests. DesignSafe-CI, Dataset, 2018; El Ghoraiby et al., Physical and mechanical properties of Ottawa F65 Sand. In B. Kutter et al. (Eds.), Model tests and numerical simulations of liquefaction and lateral spreading: LEAP-UCD-2017. New York: Springer, 2019) and direct simple shear tests (Bastidas, Ottawa F-65 Sand Characterization. PhD Dissertation, University of California, Davis, 2016) are available for Ottawa F-65 sand. The focus of this element test simulation exercise is to assess the performance of the constitutive models used by participating team in simulating the results of undrained stress-controlled cyclic triaxial tests on Ottawa F-65 sand for three different void ratios (El Ghoraiby et al., LEAP 2017: Soil characterization and element tests for Ottawa F65 sand. The George Washington University, Washington, DC, 2017; El Ghoraiby et al., LEAP-2017 GWU Laboratory Tests. DesignSafe-CI, Dataset, 2018; El Ghoraiby et al., Physical and mechanical properties of Ottawa F65 Sand. In B. Kutter et al. (Eds.), Model tests and numerical simulations of liquefaction and lateral spreading: LEAP-UCD-2017. New York: Springer, 2019). The simulated stress paths, stress-strain responses, and liquefaction strength curves show that majority of the models used in this exercise are able to provide a reasonably good match to liquefaction strength curves for the highest void ratio (0.585) but the differences between the simulations and experiments become larger for the lower void ratios (0.542 and 0.515).more » « less
An official website of the United States government

