BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Xinlu Xu

Xinlu Xu appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.plasm-phphysics.comp-phphysics.acc-phhep-phphysics.optics

Catalog footprint

What is connected

15works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2022arXiv

Efficient Generation of Tunable Magnetic and Optical Vortices Using Plasmas

Plasma is an attractive medium for generating strong microscopic magnetic structures and tunable electromagnetic radiation with predictable topologies due to its extraordinary ability to sustain and manipulate high currents and strong fields. Here, using theory and simulations, we show efficient generation of multi-megagauss magnetic and tunable optical vortices when a sharp relativistic ionization front (IF) passes through a relatively long-wavelength Laguerre-Gaussian (LG) laser pulse with orbital angular momentum (OAM). The optical vortex is frequency upshifted within a wide spectral range simply by changing the plasma density and compressed in duration. The topological charges of both vortices can be manipulated by controlling the OAM mode of the incident LG laser and/or by controlling the topology and density of the IF. For relatively high (low) plasma densities, most energy of the incident LG laser pulse is converted to the magnetic (optical) vortex.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Ultra-Bright Electron Bunch Injection in a Plasma Wakefield Driven by a Superluminal Flying Focus Electron Beam

We propose a new method for self-injection of high-quality electron bunches in the plasma wakefield structure in the blowout regime utilizing a "flying focus" produced by a drive beam with an energy chirp. In a flying focus the speed of the density centroid of the drive bunch can be superluminal or subluminal by utilizing the chromatic dependence of the focusing optics. We first derive the focal velocity and the characteristic length of the focal spot in terms of the focal length and an energy chirp. We then demonstrate using multidimensional particle-in-cell simulations that a wake driven by a superluminally propagating flying focus of an electron beam can generate GeV-level electron bunches with ultralow normalized slice emittance ($\sim$30 nm rad), high current ($\sim$ 17 kA), low slice energy-spread ($\sim$0.1%) and therefore high normalized brightness ($>10^{19}$ A/rad$^2$/m$^2$) in a plasma of density $\sim10^{19}$ cm$^{-3}$. The injection process is highly controllable and tunable by changing the focal velocity and shaping the drive beam current. Near-term experiments at FACET II where the capabilities to generate tens of kA, <10 fs drivers are planned, could potentially produce beams with brightness near $10^{20}$ A/rad$^2$/m$^2$.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Extremely Dense Gamma-Ray Pulses in Electron Beam-Multifoil Collisions

Sources of high-energy photons have important applications in almost all areas of research. However, the photon flux and intensity of existing sources is strongly limited for photon energies above a few hundred keV. Here we show that a high-current ultrarelativistic electron beam interacting with multiple submicrometer-thick conducting foils can undergo strong self-focusing accompanied by efficient emission of gamma-ray synchrotron photons. Physically, self-focusing and high-energy photon emission originate from the beam interaction with the near-field transition radiation accompanying the beam-foil collision. This near field radiation is of amplitude comparable with the beam self-field, and can be strong enough that a single emitted photon can carry away a significant fraction of the emitting electron energy. After beam collision with multiple foils, femtosecond collimated electron and photon beams with number density exceeding that of a solid are obtained. The relative simplicity, unique properties, and high efficiency of this gamma-ray source open up new opportunities for both applied and fundamental research including laserless investigations of strong-field QED processes with a single electron beam.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Generation of topologically complex three-dimensional electron beams in a plasma photocathode

Laser-triggered ionization injection is a promising way of generating controllable high-quality electrons in plasma-based acceleration. We show that ionization injection of electrons into a fully nonlinear plasma wave wake using a laser pulse comprising of one or more Laguerre-Gaussian modes with combinations of spin and orbital angular momentum can generate exotic three-dimensional (3D) spatial distributions of high-quality relativistic electrons. The phase dependent residual momenta and initial positions of the ionized electrons are encoded into their final phase space distributions, leading to complex spatiotemporal structures. The structures are formed as a result of the transverse (betatron) and longitudinal (phase slippage and energy gain) dynamics of the electrons in the wake immediately after the electrons are injected. Theoretical analysis and 3D simulations verify this mapping process leads to the generation of these complex topological beams. These beams may trigger novel beam-plasma interactions as well as produce coherent radiation with orbital angular momentum when sent through a resonant undulator.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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 $\mathbf{k}$-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. 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 $10^4$ for typical cell sizes and time steps.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

The 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 to a non-evolving laser or particle beam is calculated using 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. 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 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 are 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, a full 3D quasi-static code, and the new quasi-static PIC code with azimuthal decomposition are also presented.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Generation and acceleration of high brightness electrons beams bunched at X-ray wavelengths using plasma-based acceleration

We show using particle-in-cell (PIC) simulations and theoretical analysis that a high-quality electron beam whose density is modulated at angstrom scales can be generated directly using density downramp injection in a periodically modulated density in nonlinear plasma wave wakefields. The density modulation turns on and off the injection of electrons at the period of the modulation. Due to the unique longitudinal mapping between the electrons' initial positions and their final trapped positions inside the wake, this results in an electron beam with density modulation at a wavelength orders of magnitude shorter than the plasma density modulation. The ponderomotive force of two counter propagating lasers of the same frequency can generate a density modulation at half the laser wavelength. Assuming a laser wavelength of $0.8\micro\meter$, fully self-consistent OSIRIS PIC simulations show that this scheme can generate high quality beams modulated at wavelengths between 10s and 100 angstroms. Such beams could produce fully coherent, stable, hundreds of GW X-rays by going through a resonant undulator.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Generation of terawatt, attosecond pulses from relativistic transition radiation

When a fs duration and hundreds of kA peak current electron beam traverses the vacuum and high-density plasma interface a new process, that we call relativistic transition radiation (R-TR) generates an intense $\sim100$ as pulse containing $\sim$ TW power of coherent VUV radiation accompanied by several smaller fs duration satellite pulses. This pulse inherits the radial polarization of the incident beam field and has a ring intensity distribution. This R-TR is emitted when the beam density is comparable to the plasma density and the spot size much larger than the plasma skin depth. Physically, it arises from the return current or backward relativistic motion of electrons starting just inside the plasma that Doppler up-shifts the emitted photons. The number of R-TR pulses is determined by the number of groups of plasma electrons that originate at different depths within the first plasma wake period and emit coherently before phase mixing.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Generating high quality ultra-relativistic electron beams using an evolving electron beam driver

A new method of controllable injection to generate high quality electron bunches in the nonlinear blowout regime driven by electron beams is proposed and demonstrated using particle-in-cell simulations. Injection is facilitated by decreasing the wake phase velocity through varying the spot size of the drive beam and can be tuned through the Courant-Snyder (CS) parameters. Two regimes are examined. In the first, the spot size is focused according to the vacuum CS beta function, while in the second, it is focused by the plasma ion column. The effects of the driver intensity and vacuum CS parameters on the wake velocity and injected beam parameters are examined via theory and simulations. For plasma densities of $\sim 10^{19} ~\text{cm}^{-3}$, particle-in-cell (PIC) simulations demonstrate that peak normalized brightnesses $\gtrsim 10^{20}~\text{A}/\text{m}^2/\text{rad}^2$ can be obtained with projected energy spreads of $\lesssim 1\%$ within the middle section of the injected beam, and with normalized slice emittances as low as $\sim 10 ~\text{nm}$.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

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 are generated and interact with plasmas. However, when relativistic particles stream across the grid both in vacuum and in plasma there are many numerical issues that may arise which can lead to incorrect physics. 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. Two types of errors to these expressions are identified. The first arises from errors to the numerator of the integrand and leads to unphysical fields that are antisymmetric about the particle. The second arises from errors to the denominator of the integrand and lead to Cerenkov 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. 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. Results from the customized solver are also presented. This solution will be useful for a beam of particles that all move in one direction with a small angular divergence.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Elimination of the numerical Cerenkov instability for spectral EM-PIC codes

When using an electromagnetic particle-in-cell (EM-PIC) code to simulate a relativistically drifting plasma, a violent numerical instability known as the numerical Cerenkov instability (NCI) occurs. The NCI is due to the unphysical coupling of electromagnetic waves on a grid to wave-particle resonances, including aliased resonances, i.e., $ω+ 2πμ/Δt=(k_1+ 2πν_1/Δx_1)v_0$, where $μ$ and $ν_1$ refer to the time and space aliases and the plasma is drifting relativistically at velocity $v_0$ in the $\hat{1}$-direction. Recent studies have shown that an EM-PIC code which uses a spectral field solver and a low pass filter can eliminate the fastest growing modes of the NCI. Based on these studies a new spectral PIC code for studying laser wakefield acceleration (LWFA) in the Lorentz boosted frame was developed. However, we show that for parameters of relevance for LWFA simulations in the boosted frame, a relativistically drifting plasma is susceptible to a host of additional unstable modes with lower growth rates, and that these modes appear when the fastest growing unstable modes are filtered out. We show that these modes are most easily identified as the coupling between modes which are purely transverse (EM) and purely longitudinal (Langmuir) in the rest frame of the plasma for specific time and space aliases. We rewrite the dispersion relation of the drifting plasma for a general field solver and obtain analytic expressions for the location and growth rate for each unstable mode, i.e, for each time and space aliased resonances. We show for the spectral solver that when the fastest growing mode is eliminated a new mode at the fundamental resonance ($μ=ν_1=0$) can be seen. (Please check the whole abstract in the paper).

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Lorentz boosted frame simulation of Laser wakefield acceleration in quasi-3D geometry

When modeling laser wakefield acceleration (LWFA) using the particle-in-cell (PIC) algorithm in a Lorentz boosted frame, the plasma is drifting relativistically at $β_b c$ towards the laser, which can lead to a computational speedup of $\sim γ_b^2=(1-β_b^2)^{-1}$. Meanwhile, when LWFA is modeled in the quasi-3D geometry in which the electromagnetic fields and current are decomposed into a limited number of azimuthal harmonics, speedups are achieved by modeling three dimensional problems with the computation load on the order of two dimensional $r-z$ simulations. Here, we describe how to combine the speed ups from the Lorentz boosted frame and quasi-3D algorithms. The key to the combination is the use of a hybrid Yee-FFT solver in the quasi-3D geometry that can be used to effectively eliminate the Numerical Cerenkov Instability (NCI) that inevitably arises in a Lorentz boosted frame due to the unphysical coupling of Langmuir modes and EM modes of the relativistically drifting plasma in these simulations. In addition, based on the space-time distribution of the LWFA data in the lab and boosted frame, we propose to use a moving window to follow the drifting plasma to further reduce the computational load. We describe the details of how the NCI is eliminated for the quasi-3D geometry, the setups for simulations which combine the Lorentz boosted frame and quasi-3D geometry, the use of a moving window, and compare the results from these simulations against their corresponding lab frame cases. Good agreement is obtained, particularly when there is no self-trapping, which demonstrates it is possible to combine the Lorentz boosted frame and the quasi-3D algorithms when modeling LWFA to achieve unprecedented speedups.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Mitigation of numerical Cerenkov radiation and instability using a hybrid finite difference-FFT Maxwell solver and a local charge conserving current deposit

A hybrid Maxwell solver for fully relativistic and electromagnetic (EM) particle-in-cell (PIC) codes is described. In this solver, the EM fields are solved in $k$ space by performing an FFT in one direction, while using finite difference operators in the other direction(s). This solver eliminates the numerical Cerenkov radiation for particles moving in the preferred direction. Moreover, the numerical Cerenkov instability (NCI) induced by the relativistically drifting plasma and beam can be eliminated using this hybrid solver by applying strategies that are similar to those recently developed for pure FFT solvers. A current correction is applied for the charge conserving current deposit to correctly account for the EM calculation in hybrid Yee-FFT solver. A theoretical analysis of the dispersion properties in vacuum and in a drifting plasma for the hybrid solver is presented, and compared with PIC simulations with good agreement obtained. This hybrid solver is applied to both 2D and 3D Cartesian and quasi-3D (in which the fields and current are decomposed into azimuthal harmonics) geometries. Illustrative results for laser wakefield accelerator simulation in a Lorentz boosted frame using the hybrid solver in the 2D Cartesian geometry are presented, and compared against results from 2D UPIC-EMMA simulation which uses a pure spectral Maxwell solver, and from OSIRIS 2D lab frame simulation using the standard Yee solver. Very good agreement is obtained which demonstrates the feasibility of using the hybrid solver for high fidelity simulation of relativistically drifting plasma with no evidence of the numerical Cerenkov instability.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Modeling of Laser wakefield acceleration in Lorentz boosted frame using EM-PIC code with spectral solver

Simulating laser wakefield acceleration (LWFA) in a Lorentz boosted frame in which the plasma drifts towards the laser with $v_b$ can speedup the simulation by factors of $γ^2_b=(1-v^2_b/c^2)^{-1}$. In these simulations the relativistic drifting plasma inevitably induces a high frequency numerical instability that contaminates the interested physics. Various approaches have been proposed to mitigate this instability. One approach is to solve Maxwell equations in Fourier space (a spectral solver) as this has been shown to suppress the fastest growing modes of this instability in simple test problems using a simple low pass, ring (in two dimensions), or shell (in three dimensions) filter in Fourier space. We describe the development of a fully parallelized, multi-dimensional, particle-in-cell code that uses a spectral solver to solve Maxwell's equations and that includes the ability to launch a laser using a moving antenna. This new EM-PIC code is called UPIC-EMMA and it is based on the components of the UCLA PIC framework (UPIC). We show that by using UPIC-EMMA, LWFA simulations in the boosted frames with arbitrary $γ_b$ can be conducted without the presence of the numerical instability. We also compare the results of a few LWFA cases for several values of $γ_b$, including lab frame simulations using OSIRIS, a EM-PIC code with a finite difference time domain (FDTD) Maxwell solver. These comparisons include cases in both linear, and nonlinear regimes. We also investigate some issues associated with numerical dispersion in lab and boosted frame simulations and between FDTD and spectral solvers.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Numerical instability due to relativistic plasma drift in EM-PIC simulations

The numerical instability observed in the Electromagnetic-Particle-in-cell (EM-PIC) simulations with a plasma drifting with relativistic velocities is studied using both theory and computer simulations. We derive the numerical dispersion relation for a cold plasma drifting with a relativistic velocity and find an instability attributed to the coupling between the beam modes of the drifting plasma and the electromagnetic modes in the system. The characteristic pattern of the instability in Fourier space for various simulation setups and Maxwell Equation solvers are explored by solving the corresponding numerical dispersion relations. Furthermore, based upon these characteristic patterns we derive an asymptotic expression for the instability growth rate. The asymptotic expression greatly speeds up the calculation of instability growth rate and makes the parameter scan for minimal growth rate feasible even for full three dimensions. The results are compared against simulation results and good agreement is found. These results can be used as a guide to develop possible approaches to mitigate the instability. We examine the use of a spectral solver and show that such a solver when combined with a low pass filter with a cutoff value of $|\vec{k}|$ essentially eliminates the instability while not modifying modes of physical interest. The use of spectral solver also provides minimal errors to electromagnetic modes in the lowest Brillouin zones.

0 citations0 reviewsNo code linkedOpen access