skip to main content


Title: Particle-in-Cell Simulation of Quasi-Neutral Plasma Trapping by RF Multipole Electric Fields
Radio-frequency (RF) charged particle traps, such as the Paul trap or higher order RF multipole traps, may be used to trap quasi-neutral plasma. The presence of positive and negative plasma species mitigates the ejection of particles that occurs due to space charge repulsion. For symmetric species, such as a pair plasma, the trapped particle distribution is essentially equal for both species. For plasma with species of disparate charge-to-mass ratio, the RF parameters are chosen to directly trap the lighter species, leading to loss of the heavier species until sufficient net space charge develops to prevent further loss. Two-dimensional (2D) electrostatic particle-in-cell simulations are performed of cases with mass ratio m+/m− = 10, and also with ion–electron plasma. Multipole cases including order N = 2 (quadrupole) and higher order N = 8 (hexadecapole) are considered. The light ion-heavy ion N = 8 case exhibits particles losses less than 5% over 2500 RF periods, but the N = 8 ion–electron case exhibits a higher loss rate, likely due to non-adiabaticity of electron trajectories at the boundary, but still with low total electron loss current on the order of 10 μA. The N = 2 ion-electron case is adiabatic and stable, but is subject to a smaller trapping volume and greater initial perturbation of the bulk plasma by the trapping field.  more » « less
Award ID(s):
1806113 1619615
NSF-PAR ID:
10170175
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Physics
Volume:
1
Issue:
3
ISSN:
2624-8174
Page Range / eLocation ID:
392 to 401
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    A large aspect ratio (AR) leads to higher ion capacity in miniaturized ion trap mass spectrometers. The AR of an ion trap represents the ratio between an extended trapping dimension and the characteristic trapping dimension. In contrast to linear and rectilinear traps, changing the AR of a toroidal ion trap(TorIT) results in changes to the degree of curvature and shape of the trapping potential, and hence, on performance as a mass analyzer. SIMION simulations show that higher-order terms in the trapping potential vary strongly for small and moderate AR values (below ~10), with the effects asymptotically flattening for larger AR values. Because of the asymmetry in electrode geometry, the trapping center does not coincide with the geometric center of the trap, and this displacement also varies with AR. For instance, in the asymmetric TorIT, the saddle point in the trapping potential and the geometric trap center differ from þ0.6 to 0.4 mm depending on AR. Ion secular frequencies also change with the AR. Whereas ions in the simplified TorIT have stable trajectories for any value of AR, ions in the asymmetric TorIT become unstable at large AR values. Variations in high-order terms, the trapping center, and secular frequencies with AR are a unique feature of toroidal traps, and require significant changes in trap design and operation as the AR is changed. 
    more » « less
  2. All-solid-state batteries (ASSBs) have garnered increasing attention due to the enhanced safety, featuring nonflammable solid electrolytes as well as the potential to achieve high energy density. 1 The advancement of the ASSBs is expected to provide, arguably, the most straightforward path towards practical, high-energy, and rechargeable batteries based on metallic anodes. 1 However, the sluggish ion transmission at the cathode-electrolyte (solid/solid) interface would result in the high resistant at the contact and limit the practical implementation of these all solid-state materials in real world batteries. 2 Several methods were suggested to enhance the kinetic condition of the ion migration between the cathode and the solid electrolyte (SE). 3 A composite strategy that mixes active materials and SEs for the cathode is a general way to decrease the ion transmission barrier at the cathode-electrolyte interface. 3 The active material concentration in the cathode is reduced as much as the SE portion increases by which the energy density of the ASSB is restricted. In addition, the mixing approach generally accompanies lattice mismatches between the cathode active materials and the SE, thus providing only limited improvements, which is imputed by random contacts between the cathode active materials and the SE during the mixing process. Implementing high-pressure for the electrode and electrolyte of ASSB in the assembling process has been verified is a but effective way to boost the ion transmission ability between the cathode active materials and the SE by decreasing the grain boundary impedance. Whereas the short-circuit of the battery would be induced by the mechanical deformation of the electrolyte under high pressure. 4 Herein, we demonstrate a novel way to address the ion transmission problem at the cathode-electrolyte interface in ASSBs. Starting from the cathode configuration, the finite element method (FEM) was employed to evaluate the current concentration and the distribution of the space charge layer at the cathode-electrolyte interface. Hierarchical three-dimensional (HTD) structures are found to have a higher Li + transfer number (t Li+ ), fewer free anions, and the weaker space-charge layer at the cathode-electrolyte interface in the resulting FEM simulation. To take advantage of the HTD structure, stereolithography is adopted as a manufacturing technique and single-crystalline Ni-rich (SCN) materials are selected as the active materials. Next, the manufactured HTD cathode is sintered at 600 °C in an N 2 atmosphere for the carbonization of the resin, which induces sufficient electronic conductivity for the cathode. Then, the gel-like Li 1.4 Al 0.4 Ti 1.6 (PO 4 ) 3 (LATP) precursor is synthesized and filled into the voids of the HTD structure cathode sufficiently. And the filled HTD structure cathodes are sintered at 900 °C to achieve the crystallization of the LATP gel. Scanning transmission electron microscopy (STEM) is used to unveil the morphology of the cathode-electrolyte interface between the sintered HTD cathode and the in-situ generated electrolyte (LATP). A transient phase has been found generated at the interface and matched with both lattices of the SCN and the SE, accelerating the transmission of the Li-ions, which is further verified by density functional theory calculations. In addition, Electron Energy Loss Spectroscopy demonstrates the preserved interface between HTD cathode and SEs. Atomic force microscopy is employed to measure the potential image of the cross-sectional interface by the peak force tapping mode. The average potential of modified samples is lower than the sample that mix SCN and SEs simply in the 2D planar structure, which confirms a weakened space charge layer by the enhanced contact capability as well as the ion transmission ability. To see if the demonstrated method is universally applicable, LiNi 0.8 Co 0.1 Mn 0.1 O 2 (NCM811) is selected as the cathode active material and manufactured in the same way as the SCN. The HTD cathode based on NCM811 exhibits higher electrochemical performance compared with the reference sample based on the 2D planar mixing-type cathode. We believe such a demonstrated universal strategy provides a new guideline to engineer the cathode/electrolyte interface by revolutionizing electrode structures that can be applicable to all-solid-state batteries. Figure 1. Schematic of comparing of traditional 2D planar cathode and HTD cathode in ASSB Tikekar, M. D. , et al. , Nature Energy (2016) 1 (9), 16114 Banerjee, A. , et al. , Chem Rev (2020) 120 (14), 6878 Chen, R. , et al. , Chem Rev (2020) 120 (14), 6820 Cheng, X. , et al. , Advanced Energy Materials (2018) 8 (7) Figure 1 
    more » « less
  3. Metal-mediated cross-coupling reactions offer organic chemists a wide array of stereo- and chemically-selective reactions with broad applications in fine chemical and pharmaceutical synthesis.1 Current batch-based synthesis methods are beginning to be replaced with flow chemistry strategies to take advantage of the improved consistency and process control methods offered by continuous flow systems.2,3 Most cross-coupling chemistries still encounter several issues in flow using homogeneous catalysis, including expensive catalyst recovery and air sensitivity due to the chemical nature of the catalyst ligands.1 To mitigate some of these issues, a ligand-free heterogeneous catalysis reaction was developed using palladium (Pd) loaded into a polymeric network of a silicone elastomer, poly(hydromethylsiloxane) (PHMS), that is not air sensitive and can be used with mild reaction solvents (ethanol and water).4 In this work we present a novel method of producing soft catalytic microparticles using a multiphase flow-focusing microreactor and demonstrate their application for continuous Suzuki-Miyaura cross-coupling reactions. The catalytic microparticles are produced in a coaxial glass capillary-based 3D flow-focusing microreactor. The microreactor consists of two precursors, a cross-linking catalyst in toluene and a mixture of the PHMS polymer and a divinyl cross-linker. The dispersed phase containing the polymer, cross-linker, and cross-linking catalyst is continuously mixed and then formed into microdroplets by the continuous phase of water and surfactant (sodium dodecyl sulfate) introduced in a counter-flow configuration. Elastomeric microdroplets with a diameter ranging between 50 to 300 micron are produced at 25 to 250 Hz with a size polydispersity less than 3% in single stream production. The physicochemical properties of the elastomeric microparticles such as particle swelling/softness can be tuned using the ratio of cross-linker to polymer as well as the ratio of polymer mixture to solvent during the particle formation. Swelling in toluene can be tuned up to 400% of the initial particle volume by reducing the concentration of cross-linker in the mixture and increasing the ratio of polymer to solvent during production.5 After the particles are produced and collected, they are transferred into toluene containing palladium acetate, allowing the particles to incorporate the palladium into the polymer network and then reduce the palladium to Pd0 with the Si-H functionality present on the PHMS backbones. After the reduction, the Pd-loaded particles can be washed and dried for storage or switched into an ethanol/water solution for loading into a micro-packed bed reactor (µ-PBR) for continuous organic synthesis. The in-situ reduction of Pd within the PHMS microparticles was confirmed using energy dispersive X-ray spectroscopy (EDS), X-ray photoelectron spectroscopy (XPS) and focused ion beam-SEM, and TEM techniques. In the next step, we used the developed µ-PBR to conduct continuous organic synthesis of 4-phenyltoluene by Suzuki-Miyaura cross-coupling of 4-iodotoluene and phenylboronic acid using potassium carbonate as the base. Catalyst leaching was determined to only occur at sub ppm concentrations even at high solvent flow rates after 24 h of continuous run using inductively coupled plasma mass spectrometry (ICP-MS). The developed µ-PBR using the elastomeric microparticles is an important initial step towards the development of highly-efficient and green continuous manufacturing technologies in the pharma industry. In addition, the developed elastomeric microparticle synthesis technique can be utilized for the development of a library of other chemically cross-linkable polymer/cross-linker pairs for applications in organic synthesis, targeted drug delivery, cell encapsulation, or biomedical imaging. References 1. Ruiz-Castillo P, Buchwald SL. Applications of Palladium-Catalyzed C-N Cross-Coupling Reactions. Chem Rev. 2016;116(19):12564-12649. 2. Adamo A, Beingessner RL, Behnam M, et al. On-demand continuous-flow production of pharmaceuticals in a compact, reconfigurable system. Science. 2016;352(6281):61 LP-67. 3. Jensen KF. Flow Chemistry — Microreaction Technology Comes of Age. 2017;63(3). 4. Stibingerova I, Voltrova S, Kocova S, Lindale M, Srogl J. Modular Approach to Heterogenous Catalysis. Manipulation of Cross-Coupling Catalyst Activity. Org Lett. 2016;18(2):312-315. 5. Bennett JA, Kristof AJ, Vasudevan V, Genzer J, Srogl J, Abolhasani M. Microfluidic synthesis of elastomeric microparticles: A case study in catalysis of palladium-mediated cross-coupling. AIChE J. 2018;0(0):1-10. 
    more » « less
  4. null (Ed.)
    Particle charging in the afterglows of non-thermal plasmas typically take place in a non-neutral space charge environment. We model the same by incorporating particle-ion collision rate constant models, developed in prior work by analyzing particle-ion trajectories calculated using Langevin Dynamics simulations, into species transport equations for ions, electrons and charged particles in the afterglow. A scaling analysis of particle charging and additional Langevin Dynamics calculations of the particle-ion collision rate constant are presented to extend the range of applicability to ion electrostatic to thermal energy ratios of 300 and diffusive Knudsen number (that scales inversely with gas pressure) up to 2000. The developed collision rate constant models are first validated by comparing predictions of particle charge against measured values in a stationary, non-thermal DC plasma from past PK-4 campaigns published in Phys. Rev. Lett. 93(8): 085001 and Phys. Rev. E 72(1): 016406). The comparisons reveal excellent agreement within ±35% for particles of radius 0.6,1.0,1.3 μm in the gas pressure range of ~20-150 Pa. The experiments to probe particle charge distributions by Sharma et al. (J. Physics D: Appl. Phys. 53(24): 245204) are modeled using the validated particle-ion collision rate constant models and the calculated charge fractions are compared with measurements. The comparisons reveal that the ion/electron concentration and gas temperature in the afterglow critically influence the particle charge and the predictions are generally in qualitative agreement with the measurements. Along with critical assessment of the modeling assumptions, several recommendations are presented for future experimental design to probe charging in afterglows. 
    more » « less
  5. Abstract In the presence of a strong electric field perpendicular to the magnetic field, the electron cross-field (E × B) flow relative to the unmagnetized ions can cause the so-called electron cyclotron drift instability (ECDI) due to resonances of the ion acoustic mode and the electron cyclotron harmonics. This occurs in, for example, collisionless shock ramps in space, and in E × B discharge devices such as Hall thrusters. A prominent feature of ECDI is its capability to induce an electron flow parallel to the background E field at a speed greatly exceeding predictions by classical collision theory. Such anomalous transport is important due to its role in particle thermalization at space shocks, and in causing plasma flows towards the walls of E × B devices, leading to unfavorable erosion and performance degradation, etc. The development of ECDI and anomalous transport is often considered requiring a fully kinetic treatment. In this work, however, we demonstrate that a reduced variant of this instability, and more importantly, the associated anomalous transport, can be treated self-consistently in a collisionless two-fluid framework without any adjustable collision parameter. By treating both electron and ion species on an equal footing, the free energy due to the inter-species velocity shear allows the growth of an anomalous electron flow parallel to the background E field. We will first present linear analyses of the instability in the two-fluid five- and ten-moment models, and compare them against the fully-kinetic theory. At low temperatures, the two-fluid models predict the fastest-growing mode in good agreement with the kinetic result. Also, by including more ( > = 10 ) moments, secondary (and possibly higher) unstable branches can be recovered. The dependence of the instability on ion-to-electron mass ratio, plasma temperature, and background B field strength is also thoroughly explored. We then carry out direct numerical simulations of the cross-field setup using the five-moment model. The development of the instability, as well as the anomalous transport, is confirmed and in excellent agreement with theoretical predictions. The force balance properties are also studied using the five-moment simulation data. This work casts new insights into the nature of ECDI and the associated anomalous transport and demonstrates the potential of the two-fluid moment model in efficient modeling of E × B plasmas. 
    more » « less