Search NASA⌕ Search

Engineering topics

Li, Fei

Publications and source records attributed to Li, Fei.

At least 37 records · Page 2

Accurately simulating nine-dimensional phase space of relativistic particles in strong fields

Next-generation high-power laser systems that can be focused to ultra-high intensities exceeding 10 23 W/cm 2 are enabling new physics regimes and applications. The physics of how these lasers interact with matter is highly nonlinear, relativistic, and can involve lowest-order quantum effects. The current tool of choice for modeling these interactions is the particle-in-cell (PIC) method. In the presence of strong electromagnetic fields, the motion of charged particles and their spin is affected by radiation reaction (either the semi-classical or the quantum limit). Standard (PIC) codes usually use Boris or similar operator-splitting methods to advance the particles in standard phase space. These methods have been shown to require very small time steps in the strong-field regime in order to obtain accurate results. In addition, some problems require tracking the spin of particles, which creates a nine-dimensional (9D) particle phase space, i.e., (x, u, s). Therefore, numerical algorithms that enable high-fidelity modeling of the 9D phase space in the strong-field regime (where both the spin and momentum evolution are affected by radiation reaction) are desired. Here we present a new particle pusher that works in 9D and 6D phase space (i.e., with and without spin) based on analytical rather than leapfrog solutions to the momentum and spin advance from the Lorentz force, together with the semi-classical form of radiation reaction in the Landau-Lifshitz equation and spin evolution given by the Bargmann-Michel-Telegdi equation. Analytical solutions for the position advance are also obtained, but these are not amenable to the staggering of space and time in standard PIC codes. These analytical solutions are obtained by assuming a locally uniform and constant electromagnetic field during a time step. The solutions provide the 9D phase space advance in terms of a particle's proper time, and a mapping is used to determine the proper time step duration for each particle as a function of the lab frame time step. Due to the analytical integration of particle trajectory and spin orbit, the constraint on the time step needed to resolve trajectories in ultra-high fields can be greatly reduced. The time step required in a PIC code for accurately advancing the fields may provide additional constraints. We present single-particle simulations to show that the proposed particle pusher can greatly improve the accuracy of particle trajectories in 6D or 9D phase space for given laser fields. We have implemented the new pusher into the PIC code Osiris. Example simulations show that the proposed pusher provides improvement for a given time step. A discussion on the numerical efficiency of the proposed pusher is also provided.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Direct observation of nanoscale dynamics of ferroelectric degradation

Failure of polarization reversal, i.e., ferroelectric degradation, induced by cyclic electric loadings in ferroelectric materials, has been a long-standing challenge that negatively impacts the application of ferroelectrics in devices where reliability is critical. It is generally believed that space charges or injected charges dominate the ferroelectric degradation. However, the physics behind the phenomenon remains unclear. Here, using in-situ biasing transmission electron microscopy, we discover change of charge distribution in thin ferroelectrics during cyclic electric loadings. Charge accumulation at domain walls is the main reason of the formation of c domains, which are less responsive to the applied electric field. The rapid growth of the frozen c domains leads to the ferroelectric degradation. This finding gives insights into the nature of ferroelectric degradation in nanodevices, and reveals the role of the injected charges in polarization reversal.

42 ENGINEERING↗

High-throughput injection-acceleration of electron bunches from a linear accelerator to a laser wakefield accelerator

Plasma-based accelerators (PBAs) driven by either intense lasers (laser wakefield accelerators, LWFAs) or particle beams (plasma wakefield accelerators, PWFAs), can accelerate charged particles at extremely high gradients compared to conventional radio-frequency (RF) accelerators. In the past two decades, great strides have been made in this field, making PBA a candidate for next-generation light sources and colliders. However, these challenging applications necessarily require beams with good stability, high quality, controllable polarization and excellent reproducibility. To date, such beams are generated only by conventional RF accelerators. As such, it is important to demonstrate the injection and acceleration of beams first produced using a conventional RF accelerator, by a PBA. In some recent studies on LWFA staging and external injection-acceleration in PWFA only a very small fraction (from below 0.1% to few percent) of the injected charge (the coupling efficiency) was accelerated. For future colliders where beam energy will need to be boosted using multiple stages, the coupling efficiency per stage must approach 100%. Here we report the first demonstration of external injection from a photocathode-RF-gun-based conventional linear accelerator (LINAC) into a LWFA and subsequent acceleration without any significant loss of charge or degradation of quality, which is achieved by properly shaping and matching the beam into the plasma structure. Furthermore, this is an important step towards realizing a high-throughput, multi-stage, high-energy, hybrid conventional-plasma accelerator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

In Situ Generation of High-Energy Spin-Polarized Electrons in a Beam-Driven Plasma Wakefield Accelerator

We report that in situ generation of a high-energy, high-current, spin-polarized electron beam is an outstanding scientific challenge to the development of plasma-based accelerators for high-energy colliders. In this Letter, we show how such a spin-polarized relativistic beam can be produced by ionization injection of electrons of certain atoms with a circularly polarized laser field into a beam-driven plasma wakefield accelerator, providing a much desired one-step solution to this challenge. Using time-dependent Schrödinger equation (TDSE) simulations, we show the propensity rule of spin-dependent ionization of xenon atoms can be reversed in the strong-field multiphoton regime compared with the non-adiabatic tunneling regime, leading to high total spin polarization. Furthermore, three-dimensional particle-in-cell simulations are incorporated with TDSE simulations, providing start-to-end simulations of spin-dependent strong-field ionization of xenon atoms and subsequent trapping, acceleration, and preservation of electron spin polarization in lithium plasma. We show the generation of a high-current (0.8 kA), ultralow-normalized-emittance (~ 37 nm ), and high-energy (2.7 GeV) electron beam within just 11 cm distance, with up to ~31% net spin polarization. Higher current, energy, and net spin-polarization beams are possible by optimizing this concept, thus solving a long-standing problem facing the development of plasma accelerators.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A quasi-static particle-in-cell algorithm based on an azimuthal Fourier decomposition for highly efficient simulations of plasma-based acceleration: QPAD

The three-dimensional (3D) quasi-static particle-in-cell (PIC) algorithm is a very efficient method for modeling short-pulse laser or relativistic charged particle beam–plasma interactions. In this algorithm, the plasma response, i.e., plasma wave wake, to a non-evolving laser or particle beam is calculated using a set of Maxwell’s equations based on the quasi-static approximate equations that exclude radiation. The plasma fields are then used to advance the laser or beam forward using a large time step. The algorithm is many orders of magnitude faster than a 3D fully explicit relativistic electromagnetic PIC algorithm. It has been shown to be capable to accurately model the evolution of lasers and particle beams in a variety of scenarios. Additionally, at the same time, an algorithm in which the fields, currents and Maxwell equations are decomposed into azimuthal harmonics has been shown to reduce the algorithmic complexity of a 3D explicit PIC algorithm to that of a 2D algorithm when the expansion is truncated while maintaining accuracy for problems with near azimuthal symmetry. This hybrid algorithm uses a PIC description in r–z and a gridless description in . We describe a novel method that combines the quasi-static and hybrid PIC methods. This algorithm expands the fields, charge and current density into azimuthal harmonics. A set of the quasi-static field equations is derived for each harmonic. The complex amplitudes of the fields are then solved using the finite difference method. The beam and plasma particles are advanced in Cartesian coordinates using the total fields. Details on how this algorithm was implemented using a similar workflow to an existing quasi-static code, QuickPIC, are presented. The new code is called QPAD for QuickPIC with Azimuthal Decomposition. Benchmarks and comparisons between a fully 3D explicit PIC code (OSIRIS), a full 3D quasi-static code (QuickPIC), and the new quasi-static PIC code with azimuthal decomposition (QPAD) are also presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An ultrapotent synthetic nanobody neutralizes SARS-CoV-2 by stabilizing inactive Spike

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus enters host cells via an interaction between its Spike protein and the host cell receptor angiotensin-converting enzyme 2 (ACE2). By screening a yeast surface-displayed library of synthetic nanobody sequences, we developed nanobodies that disrupt the interaction between Spike and ACE2. Cryo–electron microscopy (cryo-EM) revealed that one nanobody, Nb6, binds Spike in a fully inactive conformation with its receptor binding domains locked into their inaccessible down state, incapable of binding ACE2. Affinity maturation and structure-guided design of multivalency yielded a trivalent nanobody, mNb6-tri, with femtomolar affinity for Spike and picomolar neutralization of SARS-CoV-2 infection. mNb6-tri retains function after aerosolization, lyophilization, and heat treatment, which enables aerosol-mediated delivery of this potent neutralizer directly to the airway epithelia.

60 APPLIED LIFE SCIENCES↗

Perovskite relaxor-PBTI03 based ferroelectric ceramics with ultrahigh dielectric and piezoelectric properties through polar nanoregions engineering

Embodiments of the invention can be directed to controlling and/or engineering the size and/or volume of polar nanoregions (PNRs) of ferroelectric polycrystalline material systems. Some embodiments can achieved this via composition modifications to cause changes in the PNRs and/or local structure. Some embodiments can be used to control and/or engineer dielectric, piezoelectric, and/or electromechanical properties of polycrystalline materials. Controlling and/or engineering the PNRs may facilitate improvements to the dielectric, piezoelectric, and/or electromechanical properties of materials. Controlling and/or engineering the PNRs may further facilitate generating a piezoelectric material that may be useful for many different piezoelectric applications.

Li, Fei↗

A new field solver for modeling of relativistic particle-laser interactions using the particle-in-cell algorithm

A customized finite-difference field solver for the particle-in-cell (PIC) algorithm that provides higher fidelity for wave-particle interactions in intense electromagnetic waves is presented. In many problems of interest, particles with relativistic energies interact with intense electromagnetic fields that have phase velocities near the speed of light. Numerical errors can arise due to (1) dispersion errors in the phase velocity of the wave, (2) the staggering in time between the electric and magnetic fields and between particle velocity and position and (3) errors in the time derivative in the momentum advance. Errors of the first two kinds are analyzed in detail. It is shown that by using field solvers with different -space operators in Faraday’s and Ampere’s law, the dispersion errors and magnetic field time-staggering errors in the particle pusher can be simultaneously removed for electromagnetic waves moving primarily in a specific direction. Here, the new algorithm was implemented into Osiris by using customized higher-order finite-difference operators. Schemes using the proposed solver in combination with different particle pushers are compared through PIC simulation. It is shown that the use of the new algorithm, together with an analytic particle pusher (assuming constant fields over a time step), can lead to accurate modeling of the motion of a single electron in an intense laser field with normalized vector potentials, eA / mc 2 , exceeding for typical cell sizes and time steps.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A molecular tandem cell for efficient solar water splitting

Artificial photosynthesis provides a way to store solar energy in chemical bonds. Achieving water splitting without an applied external potential bias provides the key to artificial photosynthetic devices. We describe here a tandem photoelectrochemical cell design that combines a dye-sensitized photoelectrosynthesis cell (DSPEC) and an organic solar cell (OSC) in a photoanode for water oxidation. When combined with a Pt electrode for H 2 evolution, the electrode becomes part of a combined electrochemical cell for water splitting, 2H 2 O → O 2 + 2H 2 , by increasing the voltage of the photoanode sufficiently to drive bias-free reduction of H + to H 2 . The combined electrode gave a 1.5% solar conversion efficiency for water splitting with no external applied bias, providing a mimic for the tandem cell configuration of PSII in natural photosynthesis. The electrode provided sustained water splitting in the molecular photoelectrode with sustained photocurrent densities of 1.24 mA/cm 2 for 1 h under 1-sun illumination with no applied bias.

14 SOLAR ENERGY↗

On numerical errors to the fields surrounding a relativistically moving particle in PIC codes

The particle-in-cell (PIC) method is widely used to model the self-consistent interaction between discrete particles and electromagnetic fields. It has been successfully applied to problems across plasma physics including plasma based acceleration, inertial confinement fusion, magnetically confined fusion, space physics, astrophysics, high energy density plasmas. In many cases the physics involves how relativistic particles (those with high relativistic γ factors) are generated and interact with plasmas. However, when relativistic particles stream across the grid, both in “vacuum” and in plasma, many numerical issues may arise which can lead to unphysical results. We present a detailed analysis of how discretized Maxwell solvers used in PIC codes can lead to numerical errors to the fields that surround particles that move at relativistic speeds across the grid. Expressions for the axial electric field as integrals in k space are presented that reveal two types of errors. The first arises from errors to the numerator of the integrand and leads to unphysical fields that are antisymmetric about the particle. Furthermore, the second arises from errors to the denominator of the integrand and lead to Cherenkov like radiation in “vacuum”. These fields are not anti-symmetric, extend behind the particle, and cause the particle to accelerate or decelerate depending on the solver and parameters. The unphysical fields are studied in detail for two representative solvers - the Yee solver and the FFT based solver. Although the Cherenkov fields are absent, the space charge fields are still present in the fundamental Brillouin zone for the FFT based solvers. In addition, the Cherenkov fields are present in higher order zones for the FFT based solvers. Comparison between the analytical solutions and PIC simulation results are presented. A solution for eliminating these unphysical fields by modifying the k operator in the axial direction is also presented. Using a customized finite difference solver, this solution was successfully implemented into OSIRIS [1]. Results from the customized solver are also presented. Additionally, this solution will be useful for a beam of particles that all move in one direction with a small angular divergence.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Solution–Processable Thermally Crosslinked Organic Radical Polymer Battery Cathodes

Organic radical polymers are promising cathode materials for next-generation batteries because of their rapid charge transfer and high cycling stability. However, these organic polymer electrodes gradually dissolve in the electrolyte, resulting in capacity fade. Several crosslinking methods have been developed to improve the performance of these electrodes, but they are either not compatible with carbon additives or compromise the solution processability of the electrodes. A one-step post-synthetic, carbon-compatible crosslinking method was developed to effectively crosslink an organic polymer electrode and allow for easy solution processing. The highest electrode capacity of 104 mAh g –1 (vs. a theoretical capacity of 111 mAh g –1 ) is achieved by introducing 1 mol % of the crosslinker, whereas the highest capacity retention (99.6 %) is obtained with 3 mol % crosslinker. In addition, mass transfer was observed in situ by using electrochemical quartz crystal microbalance with dissipation monitoring. Furthermore, these results may guide future electrode design toward fast-charging and high-capacity organic electrodes.

25 ENERGY STORAGE↗

A stable dye-sensitized photoelectrosynthesis cell mediated by a NiO overlayer for water oxidation

In the development of photoelectrochemical cells for water splitting or CO 2 reduction, a major challenge is O 2 evolution at photoelectrodes that, in behavior, mimic photosystem II. At an appropriate semiconductor electrode, a water oxidation catalyst must be integrated with a visible light absorber in a stable half-cell configuration. Here, we describe an electrode consisting of a light absorber, an intermediate electron donor layer, and a water oxidation catalyst for sustained light driven water oxidation catalysis. In assembling the electrode on nanoparticle SnO 2 /TiO 2 electrodes, a Ru(II) polypyridyl complex was used as the light absorber, NiO was deposited as an overlayer, and a Ru(II) 2,2′-bipyridine-6,6′-dicarboxylate complex as the water oxidation catalyst. In the final electrode, addition of the NiO overlayer enhanced performance toward water oxidation with the final electrode operating with a 1.1 mA/cm 2 photocurrent density for 2 h without decomposition under one sun illumination in a pH 4.65 solution. We attribute the enhanced performance to the role of NiO as an electron transfer mediator between the light absorber and the catalyst.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mechanisms underpinning the ultrahigh piezoelectricity in Sm-doped 0.705Pb(Mg 1/3 Nb 2/3 )O 3 -0.295PbTiO 3 : Temperature-induced metastable local structure and field-induced polarization rotation

The solid solution of (100 – x )%Pb(Mg 1/3 Nb 2/3 )O 3 - x %PbTiO 3 (PMN- x PT) exhibits ultrahigh piezoelectric and dielectric properties near the morphotropic phase boundary compositions and, as such, has been extensively studied in recent years. Recently, 2.5 mol. % Sm-doped PMN-29PT polycrystalline ceramics were reported to possess the highest piezoelectric coefficients (~1500 pC/N) among all reported piezoceramics, but the atomic-scale mechanisms for such high piezoelectric properties are not yet clear. Here, in situ X-ray diffraction and X-ray total scattering measurements during the application of an electric field, together with in situ total scattering measurement at different temperatures, were conducted for 2.5 mol. % Sm-doped PMN-29.5PT (2.5Sm-PMN-29.5PT). Both the largest field-induced strain and the piezoelectric response were found in the crystallites oriented with their &#12296 100 &#12297 PC directions parallel to the applied field. The local Pb displacement was analyzed using the reverse Monte Carlo method based on the pair distribution functions at different temperatures, where a temperature-induced directional change of Pb displacement was observed. Based on the experimental observations, a field-induced polarization rotation is suggested to be the dominant mechanism for the ultrahigh piezoelectricity of the 2.5Sm-PMN-29.5PT ceramic, while the ease of polarization rotation is possibly attributed to the temperature-induced metastable local monoclinic symmetries with their polar axes close to &#12296 111 &#12297 PC .

36 MATERIALS SCIENCE↗

Effects of Wall Cooling on Supersonic Modes in High-Enthalpy Hypersonic Boundary Layers over a Cone

High-enthalpy hypersonic boundary-layer transition plays an important role in many entry/descent vehicles. Aerothermodynamic performance of these vehicles strongly depends on the transition location on the surface. However, detailed transition flow physics in these chemically reacting boundary layers are poorly understood and transition estimates during the design phase rely heavily on empirically-derived transition criteria such as 𝑹𝒆(𝜽)/𝑴(𝒆). One of the most intriguing characteristics in hypersonic boundary layers is the presence of unstable supersonic modes, first identified in the 1990s. Due to a recent surge in hypersonic applications, there has been renewed interest in studying the flow physics of supersonic modes using either theory or direct numerical simulations. This paper investigates the rise of supersonic modes in a high-enthalpy hypersonic flow over a 𝟓° half-angle cone at various wall temperatures using both quasiparallel linear stability theory and linear PSE. It was found that supersonic modes exist in all wall temperature conditions including the adiabatic wall case. Cooler wall temperature causes the second Mack mode to become an unstable supersonic mode naturally downstream of the upper-branch neutral location when the wall is sufficiently cooled. In terms of the integrated growth, the second mode is still the dominant mode. Nonetheless, supersonic modes can cause additional series of relatively weaker growth than that of the second mode beyond the peak amplitude location. According to the present linear PSE results, contrary to what was speculated in the literature, supersonic-mode pressure disturbance structure radiated into the freestream is nonacoustic in nature and the formation of unstable supersonic modes is mainly associated with the synchronization of phase speed between the instability and acoustic waves in a nonp ayer, not due to nonlinear modal interaction as suggested in the literature.

Chang, Chau-Lyan↗