-
A Novel Method for Solving the Linearized 1D Vlasov--Poisson Equation
Authors:
Frank M. Lee,
B. A. Shadwick
Abstract:
We present a novel method for solving the linearized Vlasov--Poisson equation, based on analyticity properties of the equilibrium and initial condition through Cauchy-type integrals, that produces algebraic expressions for the distribution and field, i.e., the solution is expressed without integrals. Standard extant approaches involve deformations of the Bromwich contour that give erroneous result…
▽ More
We present a novel method for solving the linearized Vlasov--Poisson equation, based on analyticity properties of the equilibrium and initial condition through Cauchy-type integrals, that produces algebraic expressions for the distribution and field, i.e., the solution is expressed without integrals. Standard extant approaches involve deformations of the Bromwich contour that give erroneous results for certain physically reasonable configurations or eigenfunction expansions that are misleading as to the temporal structure of the solution. Our method is more transparent, lacks these defects, and predicts previously unrecognized behavior.
△ Less
Submitted 5 May, 2023;
originally announced May 2023.
-
Four-field Hamiltonian fluid closures of the one-dimensional Vlasov-Poisson equation
Authors:
Cristel Chandre,
Bradley A Shadwick
Abstract:
We consider a reduced dynamics for the first four fluid moments of the onedimensional Vlasov-Poisson equation, namely, the fluid density, fluid velocity, pressure and heat flux. This dynamics depends on an equation of state to close the system. This equation of state (closure) connects the fifth order moment-related to the kurtosis in velocity of the Vlasov distribution-with the first four moments…
▽ More
We consider a reduced dynamics for the first four fluid moments of the onedimensional Vlasov-Poisson equation, namely, the fluid density, fluid velocity, pressure and heat flux. This dynamics depends on an equation of state to close the system. This equation of state (closure) connects the fifth order moment-related to the kurtosis in velocity of the Vlasov distribution-with the first four moments. By solving the Jacobi identity, we derive an equation of state which ensures that the resulting reduced fluid model is Hamiltonian. We show that this Hamiltonian closure allows symmetric homogeneous equilibria of the reduced fluid model to be stable.
△ Less
Submitted 14 June, 2022;
originally announced June 2022.
-
Snowmass Whitepaper AF6: Plasma-Based Particle Sources
Authors:
M. Fuchs,
B. A. Shadwick,
N. Vafaei-Najafabadi,
A. G. R. Thomas,
G. Andonian,
M. Büscher,
A. Lehrach,
O. Apsimon,
G. Xia,
D. Filippetto,
C. B. Schroeder,
M. C. Downer
Abstract:
High-brightness beams generated by particle sources based on advanced accelerator concepts have the potential to become an essential part of future accelerator technology. High-gradient accelerators can generate and rapidly accelerate particle beams to relativistic energies while minimizing irreversible detrimental effects to the beam brightness that occur at low beam energies. Due to the high acc…
▽ More
High-brightness beams generated by particle sources based on advanced accelerator concepts have the potential to become an essential part of future accelerator technology. High-gradient accelerators can generate and rapidly accelerate particle beams to relativistic energies while minimizing irreversible detrimental effects to the beam brightness that occur at low beam energies. Due to the high accelerating gradients, these novel accelerators are also significantly more compact than conventional technology. The beam parameters of these particle sources are largely determined by the injection and subsequent acceleration processes. While there has been significant progress crucial parameters that are required for a future collider or more near-term applications, including X-ray free-electron lasers (XFELs), such as a sufficiently small energy spread and small emittance for bunches with a high charge and at high pulse repetition rate. Major research and development efforts are required to realize these approaches for a front-end injector for a future collider in order to address these limitations. In particular, this includes methods to control and manipulate the phase-space and spin degrees-of-freedom of ultrashort LWFA electron bunches with high accuracy, methods that increase the laser-to-electron beam efficiency and increased repetition rate. This also includes the development of high-resolution diagnostics, such as full 6D phase-space measurements, beam polarimetry and high-fidelity simulation tools. A further increase in beam luminosity can be achieve through emittance damping. For future colliders, the damping rings might be replaced by a substantially more compact plasma-based approach. Here, plasma wigglers are used to achieve similar damping performance but over a two orders of magnitude reduced length.
△ Less
Submitted 17 March, 2022; v1 submitted 15 March, 2022;
originally announced March 2022.
-
Linear colliders based on laser-plasma accelerators
Authors:
C. Benedetti,
S. S. Bulanov,
E. Esarey,
C. G. R. Geddes,
A. J. Gonsalves,
A. Huebl,
R. Lehe,
K. Nakamura,
C. B. Schroeder,
D. Terzani,
J. van Tilborg,
M. Turner,
J. -L. Vay,
T. Zhou,
F. Albert,
J. Bromage,
E. M. Campbell,
D. H. Froula,
J. P. Palastro,
J. Zuegel,
D. Bruhwiler,
N. M. Cook,
B. Cros,
M. C. Downer,
M. Fuchs
, et al. (18 additional authors not shown)
Abstract:
White paper to the Proceedings of the U.S. Particle Physics Community Planning Exercise (Snowmass 2021): Linear colliders based on laser-plasma accelerators
White paper to the Proceedings of the U.S. Particle Physics Community Planning Exercise (Snowmass 2021): Linear colliders based on laser-plasma accelerators
△ Less
Submitted 4 July, 2022; v1 submitted 15 March, 2022;
originally announced March 2022.
-
Transverse Oscillating Bubble Enhanced Laser-driven Betatron X-ray Radiation Generation
Authors:
Rafal Rakowski,
Ping Zhang,
Kyle Jensen,
Brendan Kettle,
Tim Kawamoto,
Sudeep Banerjee,
Colton Fruhling,
Grigory Golovin,
Daniel Haden,
Matthew S. Robinson,
Donald Umstadter,
B. A. Shadwick,
Matthias Fuchs
Abstract:
Ultrafast high-brightness X-ray pulses have proven invaluable for a broad range of research. Such pulses are typically generated via synchrotron emission from relativistic electron bunches using large-scale facilities. Recently, significantly more compact X-ray sources based on laser-wakefield accelerated (LWFA) electron beams have been demonstrated. In particular, laser-driven sources, where the…
▽ More
Ultrafast high-brightness X-ray pulses have proven invaluable for a broad range of research. Such pulses are typically generated via synchrotron emission from relativistic electron bunches using large-scale facilities. Recently, significantly more compact X-ray sources based on laser-wakefield accelerated (LWFA) electron beams have been demonstrated. In particular, laser-driven sources, where the radiation is generated by transverse oscillations of electrons within the plasma accelerator structure (so-called betatron oscillations) can generate highly-brilliant ultrashort X-ray pulses using a comparably simple setup. Here, we experimentally demonstrate a method to markedly enhance and control the parameters of LWFA-driven betatron X-ray emission. With our novel Transverse Oscillating Bubble Enhanced Betatron Radiation (TOBER) scheme, we show a significant increase in the number of generated photons by specifically manipulating the amplitude of the betatron oscillations. We realize this through an orchestrated evolution of the temporal laser pulse shape and the accelerating plasma structure. This leads to controlled off-axis injection of electrons that perform large-amplitude collective transverse betatron oscillations, resulting in increased radiation emission. Our concept holds the promise for a method to optimize the X-ray parameters for specific applications, such as time-resolved investigations with spatial and temporal atomic resolution or advanced high-resolution imaging modalities, and the generation of X-ray beams with even higher peak and average brightness.
△ Less
Submitted 2 February, 2022;
originally announced February 2022.
-
Spin-dependent two-color Kapitza-Dirac effects
Authors:
Scot McGregor,
Wayne Cheng-Wei Huang,
Bradley A. Shadwick,
Herman Batelaan
Abstract:
In this paper we present an analysis of the spin behavior of electrons propagating through a laser field. We present an experimentally realizable scenario in which spin-dependent effects of the interaction between the laser and the electrons are dominant. The laser interaction strength and incident electron velocity are in the nonrelativistic domain. This analysis may thus lead to novel methods of…
▽ More
In this paper we present an analysis of the spin behavior of electrons propagating through a laser field. We present an experimentally realizable scenario in which spin-dependent effects of the interaction between the laser and the electrons are dominant. The laser interaction strength and incident electron velocity are in the nonrelativistic domain. This analysis may thus lead to novel methods of creating and characterizing spin-polarized nonrelativistic femtosecond electron pulses.
△ Less
Submitted 14 May, 2021;
originally announced May 2021.
-
Hamiltonian fluid reduction of the 1.5D Vlasov-Maxwell equations
Authors:
Cristel Chandre,
Bradley A. Shadwick
Abstract:
We consider the Vlasov-Maxwell equations with one spatial direction and two momenta, one in the longitudinal direction and one in the transverse direction. By solving the Jacobi identity, we derive reduced Hamiltonian fluid models for the density, the fluid momenta and the second order moments, related to the pressure tensor. We also provide the Casimir invariants of the reduced Poisson bracket.…
▽ More
We consider the Vlasov-Maxwell equations with one spatial direction and two momenta, one in the longitudinal direction and one in the transverse direction. By solving the Jacobi identity, we derive reduced Hamiltonian fluid models for the density, the fluid momenta and the second order moments, related to the pressure tensor. We also provide the Casimir invariants of the reduced Poisson bracket. We show that the linearization of the equations of motion around homogeneous equilibria reproduces some essential feature of the kinetic model, the Weibel instability.
△ Less
Submitted 7 May, 2021;
originally announced May 2021.
-
Noise and error analysis and optimization in particle-based kinetic plasma simulations
Authors:
E. G. Evstatiev,
J. M. Finn,
B. A. Shadwick,
N. Hengartner
Abstract:
We analyze the noise in macro-particle methods used in plasma physics and fluid dynamics, leading to approaches for minimizing the total error, focusing on electrostatic models in one dimension. We describe kernel density estimation for continuous values of the spatial variable, $x$, and compute the covariance matrix for uniform true density. The band width of the covariance matrix is related to t…
▽ More
We analyze the noise in macro-particle methods used in plasma physics and fluid dynamics, leading to approaches for minimizing the total error, focusing on electrostatic models in one dimension. We describe kernel density estimation for continuous values of the spatial variable, $x$, and compute the covariance matrix for uniform true density. The band width of the covariance matrix is related to the width of the kernel. We find the presence of constant negative terms in the elements of the covariance matrix both on and off-diagonal. These negative correlations are related to the fact that the total number of particles is fixed at each time step. The effect of these negative correlations on the electric field, computed by Gauss's law, is that the noise in the electric field is related to a process called the Ornstein-Uhlenbeck bridge.
For non-constant density still with continuous $x$, we analyze the total error in the density estimation and discuss it in terms of bias-variance optimization (BVO). For some characteristic length $l$ and kernel width $h$, having too few particles within $h$ leads to too much variance; for $h$ large relative to $l$, there is too much smoothing of the density. The optimum between these two limits is found by BVO.
We repeat the analysis for $x$ discretized on a grid, connecting particle shapes and estimation kernels. If the particle shape satisfies a sum rule, the charge deposited on the grid is conserved exactly. Particle shapes expressed as the convolution of a kernel with another kernel that satisfies the sum rule also obey the sum rule. This property holds for kernels of arbitrary width, including widths that are not integer multiples of the grid spacing.
We show good agreement with numerical computations.
△ Less
Submitted 3 February, 2021;
originally announced February 2021.
-
Multi-color, femtosecond $γ$-ray pulse trains driven by comb-like electron beams
Authors:
S. Y. Kalmykov,
X. Davoine,
I. Ghebregziabher,
B. A. Shadwick
Abstract:
Photon engineering can be exploited to control the nonlinear evolution of the drive pulse in a laser-plasma accelerator (LPA), offering new avenues to tailor electron beam phase space on a femtosecond time scale. One promising option is to drive an LPA with an incoherent stack of two sub-Joule, multi-TW pulses of different colors. Slow self-compression of the bi-color optical driver delays electro…
▽ More
Photon engineering can be exploited to control the nonlinear evolution of the drive pulse in a laser-plasma accelerator (LPA), offering new avenues to tailor electron beam phase space on a femtosecond time scale. One promising option is to drive an LPA with an incoherent stack of two sub-Joule, multi-TW pulses of different colors. Slow self-compression of the bi-color optical driver delays electron dephasing, boosting electron beam energy without accumulation of a massive low-energy tail. The modest energy of the stack affords kHz-scale repetition rate at manageable laser average power. Propagating the stack in a preformed plasma channel induces periodic self-focusing in the trailing pulse, causing oscillations in the size of accelerating bucket. The resulting periodic injection generates, over a mm-scale distance, a train of GeV-scale electron bunches with 5D brightness exceeding $10^{17}$ A/m$^2$. This unconventional comb-like beam, with femtosecond synchronization and controllable energy spacing of components, emits, via Thomson scattering, a train of highly collimated gigawatt $γ$-ray pulses. Each pulse, corresponding to a distinct energy band between 2.5 and 25 MeV, contains over $10^6$ photons.
△ Less
Submitted 21 December, 2017;
originally announced December 2017.
-
Variational Formulation of E & M Particle Simulation Algorithms in Cylindrical Geometry using an Angular Modal Decomposition
Authors:
A. B. Stamm,
B. A. Shadwick
Abstract:
Taking advantage of the flexibility of the variational method with coordinate transformations, we derive a self-consistent set of equations of motion from a discretized Lagrangian to study kinetic plasmas using a Fourier decomposed cylindrical coordinate system. The phase-space distribution function was reduced to a collection of finite-sized macro-particles of arbitrary shape moving on a virtual…
▽ More
Taking advantage of the flexibility of the variational method with coordinate transformations, we derive a self-consistent set of equations of motion from a discretized Lagrangian to study kinetic plasmas using a Fourier decomposed cylindrical coordinate system. The phase-space distribution function was reduced to a collection of finite-sized macro-particles of arbitrary shape moving on a virtual Cartesian grid. However, the discretization of field quantities was performed in cylindrical coordinates and decomposed into a truncated Fourier series in angle. A straightforward finite element interpolation scheme is used to transform between the two grids. The equations of motion were then obtained by demanding the action be stationary. The primary advantage of the variational approach is preservation of Lagrangian symmetries. In the present case, this leads to exact energy conservation, thus avoiding possible difficulties with grid heating.
△ Less
Submitted 3 November, 2014;
originally announced November 2014.
-
Variational Formulation of Macro-Particle Models for Electromagnetic Plasma Simulations
Authors:
A. B. Stamm,
B. A. Shadwick,
E. G. Evstatiev
Abstract:
A variational method is used to derive a self-consistent macro-particle model for relativistic electromagnetic kinetic plasma simulations. Extending earlier work [E. G. Evstatiev and B. A. Shadwick, J. Comput. Phys., vol. 245, pp. 376-398, 2013], the discretization of the electromagnetic Low Lagrangian is performed via a reduction of the phase-space distribution function onto a collection of finit…
▽ More
A variational method is used to derive a self-consistent macro-particle model for relativistic electromagnetic kinetic plasma simulations. Extending earlier work [E. G. Evstatiev and B. A. Shadwick, J. Comput. Phys., vol. 245, pp. 376-398, 2013], the discretization of the electromagnetic Low Lagrangian is performed via a reduction of the phase-space distribution function onto a collection of finite-sized macro-particles of arbitrary shape and discretization of field quantities onto a spatial grid. This approach may be used with both lab frame coordinates or moving window coordinates; the latter can greatly improve computational efficiency for studying some types of laser-plasma interactions. The primary advantage of the variational approach is the preservation of Lagrangian symmetries, which in our case leads to energy conservation and thus avoids difficulties with grid heating. Additionally, this approach decouples particle size from grid spacing and relaxes restrictions on particle shape, leading to low numerical noise. The variational approach also guarantees consistent approximations in the equations of motion and is amenable to higher order methods in both space and time. We restrict our attention to the 1-1/2 dimensional case (one coordinate and two momenta). Simulations are performed with the new models and demonstrate energy conservation and low noise.
△ Less
Submitted 21 April, 2014; v1 submitted 1 October, 2013;
originally announced October 2013.
-
Variational formulation of particle algorithms for kinetic plasma simulations
Authors:
E. G. Evstatiev,
B. A. Shadwick
Abstract:
Common time-explicit numerical methods for kinetic simulations of plasmas in the low-collisions limit fall into two classes of algorithms: momentum conserving and energy conserving. Each has certain drawbacks. The PIC algorithm does not conserve total energy, which may lead to spurious numerical heating (grid heating). Its overall accuracy is at most second due to the nature of the force interpola…
▽ More
Common time-explicit numerical methods for kinetic simulations of plasmas in the low-collisions limit fall into two classes of algorithms: momentum conserving and energy conserving. Each has certain drawbacks. The PIC algorithm does not conserve total energy, which may lead to spurious numerical heating (grid heating). Its overall accuracy is at most second due to the nature of the force interpolation between grid and particle position. Energy-conserving algorithms do not exhibit grid heating, but because their formulation uses potentials, computationally undesirable matrix inversions may be necessary. In addition, compared to PIC algorithms for the same accuracy, these algorithms have higher numerical noise due to the restricted choice of particle shapes. Here we formulate time-explicit, finite-size particle algorithms using particular reductions of the particle distribution function. These reductions are used in two variational principles, a Lagrangian-based and a Hamiltonian-based in conjunction with a non-canonical Poisson bracket. The Lagrangian formulations here generalize previous such formulations. The Hamiltonian formulation is presented here for the first time. Many drawbacks of the two classes of particle methods are mitigated. For example, restrictions on particle shapes are relaxed in energy conserving algorithms, which allows to decrease the numerical noise in these methods. The Hamiltonian formulation of particle algorithms is done in terms of fields instead of potentials, thus avoiding solving Poisson's equation. An algorithm that conserves both energy and momentum is presented. Other features of the algorithms include a natural way to perform coordinate transformations, the use of various time integrating methods, and the ability to increase the overall accuracy beyond second order, including all generalizations.
△ Less
Submitted 13 December, 2012; v1 submitted 13 October, 2012;
originally announced October 2012.
-
Computationally efficient methods for modelling laser wakefield acceleration in the blowout regime
Authors:
B. M. Cowan,
S. Y. Kalmykov,
A. Beck,
X. Davoine,
K. Bunkers,
A. F. Lifschitz,
E. Lefebvre,
D. L. Bruhwiler,
B. A. Shadwick,
D. P. Umstadter
Abstract:
Electron self-injection and acceleration until dephasing in the blowout regime is studied for a set of initial conditions typical of recent experiments with 100 terawatt-class lasers. Two different approaches to computationally efficient, fully explicit, three-dimensional particle-in-cell modelling are examined. First, the Cartesian code VORPAL using a perfect-dispersion electromagnetic solver pre…
▽ More
Electron self-injection and acceleration until dephasing in the blowout regime is studied for a set of initial conditions typical of recent experiments with 100 terawatt-class lasers. Two different approaches to computationally efficient, fully explicit, three-dimensional particle-in-cell modelling are examined. First, the Cartesian code VORPAL using a perfect-dispersion electromagnetic solver precisely describes the laser pulse and bubble dynamics, taking advantage of coarser resolution in the propagation direction, with a proportionally larger time step. Using third-order splines for macroparticles helps suppress the sampling noise while keeping the usage of computational resources modest. The second way to reduce the simulation load is using reduced-geometry codes. In our case, the quasi-cylindrical code CALDER-CIRC uses decomposition of fields and currents into a set of poloidal modes, while the macroparticles move in the Cartesian 3D space. Cylindrical symmetry of the interaction allows using just two modes, reducing the computational load to roughly that of a planar Cartesian simulation while preserving the 3D nature of the interaction. This significant economy of resources allows using fine resolution in the direction of propagation and a small time step, making numerical dispersion vanishingly small, together with a large number of particles per cell, enabling good particle statistics. Quantitative agreement of the two simulations indicates that they are free of numerical artefacts. Both approaches thus retrieve physically correct evolution of the plasma bubble, recovering the intrinsic connection of electron self-injection to the nonlinear optical evolution of the driver.
△ Less
Submitted 3 April, 2012;
originally announced April 2012.
-
Spectral reduction: a statistical description of turbulence
Authors:
John C. Bowman,
B. A. Shadwick,
P. J. Morrison
Abstract:
A method is described for predicting statistical properties of turbulence. Collections of Fourier amplitudes are represented by nonuniformly spaced modes with enhanced coupling coefficients. The statistics of the full dynamics can be recovered from the time-averaged predictions of the reduced model. A Liouville theorem leads to inviscid equipartition solutions. Excellent agreement is obtained wi…
▽ More
A method is described for predicting statistical properties of turbulence. Collections of Fourier amplitudes are represented by nonuniformly spaced modes with enhanced coupling coefficients. The statistics of the full dynamics can be recovered from the time-averaged predictions of the reduced model. A Liouville theorem leads to inviscid equipartition solutions. Excellent agreement is obtained with two-dimensional forced-dissipative pseudospectral simulations. For the two-dimensional enstrophy cascade, logarithmic corrections to the high-order structure functions are observed.
△ Less
Submitted 5 November, 1999; v1 submitted 5 March, 1999;
originally announced March 1999.
-
Exactly Conservative Integrators
Authors:
B. A. Shadwick,
John C. Bowman,
P. J. Morrison
Abstract:
Traditional numerical discretizations of conservative systems generically yield an artificial secular drift of any nonlinear invariants. In this work we present an explicit nontraditional algorithm that exactly conserves these invariants. We illustrate the general method by applying it to the three-wave truncation of the Euler equations, the Lotka--Volterra predator--prey model, and the Kepler pro…
▽ More
Traditional numerical discretizations of conservative systems generically yield an artificial secular drift of any nonlinear invariants. In this work we present an explicit nontraditional algorithm that exactly conserves these invariants. We illustrate the general method by applying it to the three-wave truncation of the Euler equations, the Lotka--Volterra predator--prey model, and the Kepler problem. This method is discussed in the context of symplectic (phase space conserving) integration methods as well as nonsymplectic conservative methods. We comment on the application of our method to general conservative systems.
△ Less
Submitted 2 August, 1995;
originally announced July 1995.