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

Frank S. Tsung

Frank S. Tsung 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-ph

Catalog footprint

What is connected

11works
3topics
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

11 published item(s)

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
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

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
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

Exploiting multi-scale parallelism for large scale numerical modelling of laser wakefield accelerators

A new generation of laser wakefield accelerators, supported by the extreme accelerating fields generated in the interaction of PW-Class lasers and underdense targets, promises the production of high quality electron beams in short distances for multiple applications. Achieving this goal will rely heavily on numerical modeling for further understanding of the underlying physics and identification of optimal regimes, but large scale modeling of these scenarios is computationally heavy and requires efficient use of state-of-the-art Petascale supercomputing systems. We discuss the main difficulties involved in running these simulations and the new developments implemented in the OSIRIS framework to address these issues, ranging from multi-dimensional dynamic load balancing and hybrid distributed / shared memory parallelism to the vectorization of the PIC algorithm. We present the results of the OASCR Joule Metric program on the issue of large scale modeling of LWFA, demonstrating speedups of over 1 order of magnitude on the same hardware. Finally, scalability to over $\sim 10^6$ cores, and sustained performance over $\sim 2$ PFlops is demonstrated, opening the way for large scale modeling of laser wakefield accelerator scenarios.

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

Convective Raman Amplification of Light Pulses Causing Kinetic Inflation in Inertial Fusion Plasmas

We perform 1D particle-in-cell (PIC) simulations using OSIRIS, which model a short-duration (~500/ω0 FWHM) scattered light seed pulse in the presence of a constant counter-propagating pump laser with an intensity far below the absolute instability threshold. The seed undergoes linear convective Raman amplification and dominates over fluctuations due to particle discreteness. Our simulation results are in good agreement with results from a coupled mode solver when we take into account special relativity and the use of finite size PIC simulation particles. We present linear gain spectra including both effects. Extending the PIC simulations past when the seed exits the simulation domain reveals bursts of large-amplitude scattering in many cases, which does not occur in simulations without the seed pulse. These bursts can have amplitudes several times greater than the amplified seed pulse, and we demonstrate that this large-amplitude scattering is the result of kinetic inflation by examining trapped particle orbits. This large-amplitude scattering is caused by the seed modifying the distribution function earlier in the simulation. We perform some simulations with longer duration seeds, which lead to parts of the seeds undergoing kinetic inflation and reaching amplitudes several times more than the steady-state linear theory results. Simulations with continuous seeds demonstrate that the onset of inflation depends on seed wavelength and incident intensity, and we observe oscillations in the reflectivity at a frequency equal to the difference between the seed frequency and the frequency at which the inflationary SRS grows.

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