Showing 1–50 of 67 results for author: Vay, J

Application of mesh refinement to relativistic magnetic reconnection

Authors: Revathi Jambunathan, Henry Jones, Lizzette Corrales, Hannah Klion, Michael Rowan, Andrew Myers, Weiqun Zhang, Jean-Luc Vay

Abstract: During relativistic magnetic reconnection, antiparallel magnetic fields undergo a rapid change in topology, releasing a large amount of energy in the form of non-thermal particle acceleration. This work explores the application of mesh refinement to 2D reconnection simulations to efficiently model the ineherent disparity in length-scales. We have systematically investigated the effects of mesh ref… ▽ More During relativistic magnetic reconnection, antiparallel magnetic fields undergo a rapid change in topology, releasing a large amount of energy in the form of non-thermal particle acceleration. This work explores the application of mesh refinement to 2D reconnection simulations to efficiently model the ineherent disparity in length-scales. We have systematically investigated the effects of mesh refinement and determined necessary modifications to the algorithm required to mitigate non-physical artifacts at the coarse-fine interface. We have used the ultrahigh-order Pseudo-Spectral Analytical Time-Domain (PSATD) Maxwell solver to analyze how its use can mitigate the numerical dispersion that occurs with the finite-difference time-domain (FDTD) (or ``Yee'') method. Absorbing layers are introduced at the coarse-fine interface to eliminate spurious effects that occur with mesh refinement. We also study how damping the electromagnetic fields and current density in the absorbing layer can help prevent the non-physical accumulation of charge and current density at the coarse-fine interface. Using a mesh refinement ratio of 8 for two-dimensional magnetic reconnection simulations, we obtained good agreement with the high resolution baseline simulation, using only 36% of the macroparticles and 71% of the node-hours needed for the baseline. The methods presented here are especially applicable to 3D systems where higher memory savings are expected than in 2D, enabling comprehensive, computationally efficient 3D reconnection studies in the future. △ Less

Submitted 16 August, 2024; originally announced August 2024.

Report on the Advanced Linear Collider Study Group (ALEGRO) Workshop 2024

Authors: J. Vieira, B. Cros, P. Muggli, I. A. Andriyash, O. Apsimon, M. Backhouse, C. Benedetti, S. S. Bulanov, A. Caldwell, Min Chen, V. Cilento, S. Corde, R. D'Arcy, S. Diederichs, E. Ericson, E. Esarey, J. Farmer, L. Fedeli, A. Formenti, B. Foster, M. Garten, C. G. R. Geddes, T. Grismayer, M. J. Hogan, S. Hooker , et al. (19 additional authors not shown)

Abstract: The workshop focused on the application of ANAs to particle physics keeping in mind the ultimate goal of a collider at the energy frontier (10\,TeV, e$^+$/e$^-$, e$^-$/e$^-$, or $γγ$). The development of ANAs is conducted at universities and national laboratories worldwide. The community is thematically broad and diverse, in particular since lasers suitable for ANA research (multi-hundred-terawatt… ▽ More The workshop focused on the application of ANAs to particle physics keeping in mind the ultimate goal of a collider at the energy frontier (10\,TeV, e$^+$/e$^-$, e$^-$/e$^-$, or $γγ$). The development of ANAs is conducted at universities and national laboratories worldwide. The community is thematically broad and diverse, in particular since lasers suitable for ANA research (multi-hundred-terawatt peak power, a few tens of femtosecond-long pulses) and acceleration of electrons to hundreds of mega electron volts to multi giga electron volts became commercially available. The community spans several continents (Europe, America, Asia), including more than 62 laboratories in more than 20 countries. It is among the missions of the ICFA-ANA panel to feature the amazing progress made with ANAs, to provide international coordination and to foster international collaborations towards a future HEP collider. The scope of this edition of the workshop was to discuss the recent progress and necessary steps towards realizing a linear collider for particle physics based on novel-accelerator technologies (laser or beam driven in plasma or structures). Updates on the relevant aspects of the European Strategy for Particle Physics (ESPP) Roadmap Process as well as of the P5 (in the US) were presented, and ample time was dedicated to discussions. The major outcome of the workshop is the decision for ALEGRO to coordinate efforts in Europe, in the US, and in Asia towards a pre-CDR for an ANA-based, 10\,TeV CM collider. This goal of this coordination is to lead to a funding proposal to be submitted to both EU and EU/US funding agencies. This document presents a summary of the workshop, as seen by the co-chairs, as well as short 'one-pagers' written by the presenters at the workshop. △ Less

Submitted 15 August, 2024; v1 submitted 6 August, 2024; originally announced August 2024.

Comments: 72 pages

Comparison of WarpX and GUINEA-PIG for electron positron collisions

Authors: Bao Nguyen, Arianna Formenti, Remi Lehe, Jean-Luc Vay, Spencer Gessner, Luca Fedeli

Abstract: As part of the Snowmass'21 planning exercise, the Advanced Accelerator Concepts community proposed developing multi-TeV linear colliders and considered beam-beam effects for these machines. Such colliders operate under a high disruption regime with an enormous number of electron-positron pairs produced from QED effects. Thus, it requires a self-consistent treatment of the fields produced by the pa… ▽ More As part of the Snowmass'21 planning exercise, the Advanced Accelerator Concepts community proposed developing multi-TeV linear colliders and considered beam-beam effects for these machines. Such colliders operate under a high disruption regime with an enormous number of electron-positron pairs produced from QED effects. Thus, it requires a self-consistent treatment of the fields produced by the pairs, which is not implemented in state-of-the-art beam-beam codes such as GUINEA-PIG. WarpX is a parallel, open-source, and portable particle-in-cell code with an active developer community that models QED processes with photon and pair generation in relativistic laser-beam interactions. However, its application to beam-beam collisions has yet to be fully explored. In this work, we benchmark the luminosity spectra, photon spectra, and coherent production process from WarpX against GUINEA-PIG in the ILC and ultra-tight collision scenarios. Our performance comparison demonstrates a significant speed-up advantage of WarpX, ensuring a more robust and efficient modeling of electron-positron collisions at multi-TeV energies. △ Less

Submitted 14 May, 2024; originally announced May 2024.

Comments: 3 pages conference proceeding. 15th International Particle Accelerator Conference (IPAC'24). Paper ID is WEPC84

Report number: doi:10.18429/JACoW-IPAC2024-WEPC84

Journal ref: International Particle Accelerator Conference 2024

LASY: LAser manipulations made eaSY

Authors: Maxence Thévenet, Igor A. Andriyash, Luca Fedeli, Ángel Ferran Pousa, Axel Huebl, Sören Jalas, Manuel Kirchen, Remi Lehe, Rob J. Shalloo, Alexander Sinn, Jean-Luc Vay

Abstract: Using realistic laser profiles for simulations of laser-plasma interaction is critical to reproduce experimental measurements, but the interface between experiments and simulations can be challenging. Similarly, start-to-end simulations with different codes may require error-prone manipulations to convert between different representations of a laser pulse. In this work, we propose LASY, an open-so… ▽ More Using realistic laser profiles for simulations of laser-plasma interaction is critical to reproduce experimental measurements, but the interface between experiments and simulations can be challenging. Similarly, start-to-end simulations with different codes may require error-prone manipulations to convert between different representations of a laser pulse. In this work, we propose LASY, an open-source Python library to simplify these workflows. Developed through an international collaboration between experimental, theoretical and computational physicists, LASY can be used to initialize a laser profile from an experimental measurement, from a simulation, or from analytics, manipulate it, and write it into a file in compliance with the openPMD standard. This profile can then be used as an input of a simulation code. △ Less

Submitted 18 March, 2024; originally announced March 2024.

Comments: 6 pages, 4 figures, EAAC2023 conference proceedings

Synthesizing Particle-in-Cell Simulations Through Learning and GPU Computing for Hybrid Particle Accelerator Beamlines

Authors: Ryan T. Sandberg, Remi Lehe, Chad E. Mitchell, Marco Garten, Andrew Myers, Ji Qiang, Jean-Luc Vay, Axel Huebl

Abstract: Particle accelerator modeling is an important field of research and development, essential to investigating, designing and operating some of the most complex scientific devices ever built. Kinetic simulations of relativistic, charged particle beams and advanced plasma accelerator elements are often performed with high-fidelity particle-in-cell simulations, some of which fill the largest GPU superc… ▽ More Particle accelerator modeling is an important field of research and development, essential to investigating, designing and operating some of the most complex scientific devices ever built. Kinetic simulations of relativistic, charged particle beams and advanced plasma accelerator elements are often performed with high-fidelity particle-in-cell simulations, some of which fill the largest GPU supercomputers. Start-to-end modeling of a particle accelerator includes many elements and it is desirable to integrate and model advanced accelerator elements fast, in effective models. Traditionally, analytical and reduced-physics models fill this role. The vast data from high-fidelity simulations and power of GPU-accelerated computation open a new opportunity to complement traditional modeling without approximations: surrogate modeling through machine learning. In this paper, we implement, present and benchmark such a data-driven workflow, synthesising a fully GPU-accelerated, conventional-surrogate simulation for hybrid particle accelerator beamlines. △ Less

Submitted 30 April, 2024; v1 submitted 27 February, 2024; originally announced February 2024.

Comments: 11 pages, 10 figures. Accepted to PASC24 proceedings

Journal ref: Proceedings of the Platform for Advanced Scientific Computing Conference PASC24 (2024)

PIC-JR$_{ho}$m: a pseudo-spectral Particle-In-Cell formulation with arbitrary charge and current densities time dependencies for the modeling of relativistic plasmas

Authors: Olga Shapoval, Edoardo Zoni, Remi Lehe, Maxence Thevenet, Jean-Luc Vay

Abstract: This paper introduces a novel formulation of the Particle-In-Cell (PIC) method for the modeling of relativistic plasmas, which leverages the ability of the Pseudo-Spectral Analytical Time-Domain solver (PSATD) to handle arbitrary time dependencies of the charge and current densities during one PIC cycle. The new formulation is applied to a modified set of Maxwell's equations that was proposed earl… ▽ More This paper introduces a novel formulation of the Particle-In-Cell (PIC) method for the modeling of relativistic plasmas, which leverages the ability of the Pseudo-Spectral Analytical Time-Domain solver (PSATD) to handle arbitrary time dependencies of the charge and current densities during one PIC cycle. The new formulation is applied to a modified set of Maxwell's equations that was proposed earlier in the context of divergence cleaning, and to recently proposed extensions of the PSATD-PIC algorithm. Detailed analysis and testings revealed that, under some condition, the new formulation can expand the range of numerical parameters under which PIC simulations are stable and accurate when modeling relativistic plasmas such as, e.g., plasma-based particle accelerators. △ Less

Submitted 26 December, 2023; originally announced December 2023.

State of In Situ Visualization in Simulations: We are fast. But are we inspiring?

Authors: Axel Huebl, Arianna Formenti, Marco Garten, Jean-Luc Vay

Abstract: Visualization of dynamic processes in scientific high-performance computing is an immensely data intensive endeavor. Application codes have recently demonstrated scaling to full-size Exascale machines, and generating high-quality data for visualization is consequently on the machine-scale, easily spanning 100s of TBytes of input to generate a single video frame. In situ visualization, the techniqu… ▽ More Visualization of dynamic processes in scientific high-performance computing is an immensely data intensive endeavor. Application codes have recently demonstrated scaling to full-size Exascale machines, and generating high-quality data for visualization is consequently on the machine-scale, easily spanning 100s of TBytes of input to generate a single video frame. In situ visualization, the technique to consume the many-node decomposed data in-memory, as exposed by applications, is the dominant workflow. Although in situ visualization has achieved tremendous progress in the last decade, scaling to system-size together with the application codes that produce its data, there is one important question that we cannot skip: is what we produce insightful and inspiring? △ Less

Submitted 30 September, 2023; originally announced October 2023.

Comments: 2 pages + references, 1 figure, accepted lightning talk abstract for ISAV23 (in conjunction with SC23)

Laser-Plasma Ion Beam Booster Based on Hollow-Channel Magnetic Vortex Acceleration

Authors: Marco Garten, Stepan S. Bulanov, Sahel Hakimi, Lieselotte Obst-Huebl, Chad E. Mitchell, Carl Schroeder, Eric Esarey, Cameron G. R. Geddes, Jean-Luc Vay, Axel Huebl

Abstract: Laser-driven ion acceleration provides ultra-short, high-charge, low-emittance beams, which are desirable for a wide range of high-impact applications. Yet after decades of research, a significant increase in maximum ion energy is still needed. This work introduces a quality-preserving staging concept for ultra-intense ion bunches that is seamlessly applicable from the non-relativistic plasma sour… ▽ More Laser-driven ion acceleration provides ultra-short, high-charge, low-emittance beams, which are desirable for a wide range of high-impact applications. Yet after decades of research, a significant increase in maximum ion energy is still needed. This work introduces a quality-preserving staging concept for ultra-intense ion bunches that is seamlessly applicable from the non-relativistic plasma source to the relativistic regime. Full 3D particle-in-cell simulations prove robustness and capture of a high-charge proton bunch, suitable for readily available and near-term laser facilities. △ Less

Submitted 5 May, 2024; v1 submitted 9 August, 2023; originally announced August 2023.

Comments: 4 pages, 4 figures, under peer review

Journal ref: Physical Review Research 6, 033148 (2024)

Beam Delivery and Beamstrahlung Considerations for Ultra-High Energy Linear Colliders

Authors: Tim Barklow, Spencer Gessner, Mark Hogan, Cho-Kuen Ng, Michael Peskin, Tor Raubenheimer, Glen White, Erik Adli, Gevy Jiawei Cao, Carl A. Lindstrom, Kyrre Sjobak, Sam Barber, Cameron Geddes, Arianna Formenti, Remi Lehe, Carl Schroeder, Davide Terzani, Jeroen van Tilborg, Jean-Luc Vay, Edoardo Zoni, Chris Doss, Michael Litos, Ihar Lobach, John Power, Maximilian Swiatlowski , et al. (5 additional authors not shown)

Abstract: As part of the Snowmass'21 community planning excercise, the Advanced Accelerator Concepts (AAC) community proposed future linear colliders with center-of-mass energies up to 15 TeV and luminosities up to 50$\times10^{34}$ cm$^{-2}$s$^{-1}$ in a compact footprint. In addition to being compact, these machines must also be energy efficient. We identify two challenges that must be addressed in the de… ▽ More As part of the Snowmass'21 community planning excercise, the Advanced Accelerator Concepts (AAC) community proposed future linear colliders with center-of-mass energies up to 15 TeV and luminosities up to 50$\times10^{34}$ cm$^{-2}$s$^{-1}$ in a compact footprint. In addition to being compact, these machines must also be energy efficient. We identify two challenges that must be addressed in the design of these machines. First, the Beam Delivery System (BDS) must not add significant length to the accelerator complex. Second, beam parameters must be chosen to mitigate beamstrahlung effects and maximize the luminosity-per-power of the machine. In this paper, we review advances in plasma lens technology that will help to reduce the length of the BDS system and we detail new Particle-in-Cell simulation studies that will provide insight into beamstrahlung mitigation techniques. We apply our analysis to both $e^+e^-$ and $γγ$ colliders. △ Less

Submitted 31 May, 2023; v1 submitted 30 April, 2023; originally announced May 2023.

Particle-in-Cell Simulations of Relativistic Magnetic Reconnection with Advanced Maxwell Solver Algorithms

Authors: Hannah Klion, Revathi Jambunathan, Michael E. Rowan, Eloise Yang, Donald Willcox, Jean-Luc Vay, Remi Lehe, Andrew Myers, Axel Huebl, Weiqun Zhang

Abstract: Relativistic magnetic reconnection is a non-ideal plasma process that is a source of non-thermal particle acceleration in many high-energy astrophysical systems. Particle-in-cell (PIC) methods are commonly used for simulating reconnection from first principles. While much progress has been made in understanding the physics of reconnection, especially in 2D, the adoption of advanced algorithms and… ▽ More Relativistic magnetic reconnection is a non-ideal plasma process that is a source of non-thermal particle acceleration in many high-energy astrophysical systems. Particle-in-cell (PIC) methods are commonly used for simulating reconnection from first principles. While much progress has been made in understanding the physics of reconnection, especially in 2D, the adoption of advanced algorithms and numerical techniques for efficiently modeling such systems has been limited. With the GPU-accelerated PIC code WarpX, we explore the accuracy and potential performance benefits of two advanced Maxwell solver algorithms: a non-standard finite difference scheme (CKC) and an ultrahigh-order pseudo-spectral method (PSATD). We find that for the relativistic reconnection problem, CKC and PSATD qualitatively and quantitatively match the standard Yee-grid finite-difference method. CKC and PSATD both admit a time step that is 40% longer than Yee, resulting in a ~40% faster time to solution for CKC, but no performance benefit for PSATD when using a current deposition scheme that satisfies Gauss's law. Relaxing this constraint maintains accuracy and yields a 30% speedup. Unlike Yee and CKC, PSATD is numerically stable at any time step, allowing for a larger time step than with the finite-difference methods. We found that increasing the time step 2.4-3 times over the standard Yee step still yields accurate results, but only translates to modest performance improvements over CKC due to the current deposition scheme used with PSATD. Further optimization of this scheme will likely improve the effective performance of PSATD. △ Less

Submitted 20 April, 2023; originally announced April 2023.

Comments: 19 pages, 10 figures. Submitted to ApJ

Light-Matter Interaction Near the Schwinger Limit Using Tightly Focused Doppler-Boosted Lasers

Authors: Neïl Zaïm, Antonin Sainte-Marie, Luca Fedeli, Pierre Bartoli, Axel Huebl, Jean-Luc Vay, Henri Vincenti

Abstract: The Schwinger limit could be approached by focusing to its diffraction limit the light reflected by a plasma mirror irradiated by a multi-petawatt laser. We explore numerically the interaction between such intense light and matter. We find that the interaction with a relativistic counterpropagative electron beam would enable the exploration of the fully nonperturbative regime of strong-field quant… ▽ More The Schwinger limit could be approached by focusing to its diffraction limit the light reflected by a plasma mirror irradiated by a multi-petawatt laser. We explore numerically the interaction between such intense light and matter. We find that the interaction with a relativistic counterpropagative electron beam would enable the exploration of the fully nonperturbative regime of strong-field quantum electrodynamics (SF-QED), while the interaction with an initially solid target leads to a profusion of SF-QED effects that retroact on the laser-plasma interaction. We observe in both scenarios the formation of relativistic attosecond electron-positron jets with very high densities. △ Less

Submitted 30 March, 2023; originally announced March 2023.

Comments: 7 pages, 3 figures

From Compact Plasma Particle Sources to Advanced Accelerators with Modeling at Exascale

Authors: Axel Huebl, Remi Lehe, Edoardo Zoni, Olga Shapoval, Ryan T. Sandberg, Marco Garten, Arianna Formenti, Revathi Jambunathan, Prabhat Kumar, Kevin Gott, Andrew Myers, Weiqun Zhang, Ann Almgren, Chad E. Mitchell, Ji Qiang, David Grote, Alexander Sinn, Severin Diederichs, Maxence Thevenet, Luca Fedeli, Thomas Clark, Neil Zaim, Henri Vincenti, Jean-Luc Vay

Abstract: Developing complex, reliable advanced accelerators requires a coordinated, extensible, and comprehensive approach in modeling, from source to the end of beam lifetime. We present highlights in Exascale Computing to scale accelerator modeling software to the requirements set for contemporary science drivers. In particular, we present the first laser-plasma modeling on an exaflop supercomputer using… ▽ More Developing complex, reliable advanced accelerators requires a coordinated, extensible, and comprehensive approach in modeling, from source to the end of beam lifetime. We present highlights in Exascale Computing to scale accelerator modeling software to the requirements set for contemporary science drivers. In particular, we present the first laser-plasma modeling on an exaflop supercomputer using the US DOE Exascale Computing Project WarpX. Leveraging developments for Exascale, the new DOE SCIDAC-5 Consortium for Advanced Modeling of Particle Accelerators (CAMPA) will advance numerical algorithms and accelerate community modeling codes in a cohesive manner: from beam source, over energy boost, transport, injection, storage, to application or interaction. Such start-to-end modeling will enable the exploration of hybrid accelerators, with conventional and advanced elements, as the next step for advanced accelerator modeling. Following open community standards, we seed an open ecosystem of codes that can be readily combined with each other and machine learning frameworks. These will cover ultrafast to ultraprecise modeling for future hybrid accelerator design, even enabling virtual test stands and twins of accelerators that can be used in operations. △ Less

Submitted 18 April, 2023; v1 submitted 22 March, 2023; originally announced March 2023.

Comments: 4 pages, 3 figures, presented at the 20th Advanced Accelerator Concepts Workshop (AAC22)

Bayesian optimization of laser-plasma accelerators assisted by reduced physical models

Authors: A. Ferran Pousa, S. Jalas, M. Kirchen, A. Martinez de la Ossa, M. Thévenet, S. Hudson, J. Larson, A. Huebl, J. -L. Vay, R. Lehe

Abstract: Particle-in-cell simulations are among the most essential tools for the modeling and optimization of laser-plasma accelerators, since they reproduce the physics from first principles. However, the high computational cost associated with them can severely limit the scope of parameter and design optimization studies. Here, we show that a multitask Bayesian optimization algorithm can be used to mitig… ▽ More Particle-in-cell simulations are among the most essential tools for the modeling and optimization of laser-plasma accelerators, since they reproduce the physics from first principles. However, the high computational cost associated with them can severely limit the scope of parameter and design optimization studies. Here, we show that a multitask Bayesian optimization algorithm can be used to mitigate the need for such high-fidelity simulations by incorporating information from inexpensive evaluations of reduced physical models. In a proof-of-principle study, where a high-fidelity optimization with FBPIC is assisted by reduced-model simulations with Wake-T, the algorithm demonstrates an order-of-magnitude speedup. This opens a path for the cost-effective optimization of laser-plasma accelerators in large parameter spaces, an important step towards fulfilling the high beam quality requirements of future applications. △ Less

Submitted 23 December, 2022; originally announced December 2022.

RF Accelerator Technology R&D: Report of AF7-rf Topical Group to Snowmass 2021

Authors: Sergey Belomestnykh, Emilio A. Nanni, Hans Weise, Sergey V. Baryshev, Pashupati Dhakal, Rongli Geng, Bianca Giaccone, Chunguang Jing, Matthias Liepe, Xueying Lu, Tianhuan Luo, Ganapati Myneni, Alireza Nassiri, David Neuffer, Cho-Kuen Ng, Sam Posen, Sami Tantawi, Anne-Marie Valente-Feliciano, Jean-Luc Vay, Brandon Weatherford, Akira Yamamoto

Abstract: Accelerator radio frequency (RF) technology has been and remains critical for modern high energy physics (HEP) experiments based on particle accelerators. Tremendous progress in advancing this technology has been achieved over the past decade in several areas highlighted in this report. These achievements and new results expected from continued R&D efforts could pave the way for upgrades of existi… ▽ More Accelerator radio frequency (RF) technology has been and remains critical for modern high energy physics (HEP) experiments based on particle accelerators. Tremendous progress in advancing this technology has been achieved over the past decade in several areas highlighted in this report. These achievements and new results expected from continued R&D efforts could pave the way for upgrades of existing facilities, improvements to accelerators already under construction (e.g., PIP-II), well-developed proposals (e.g., ILC, CLIC), and/or enable concepts under development, such as FCC-ee, CEPC, C3, HELEN, multi-MW Fermilab Proton Intensity Upgrade, future Muon Colloder, etc. Advances in RF technology have impact beyond HEP on accelerators built for nuclear physics, basic energy sciences, and other areas. Recent examples of such accelerators are European XFEL, LCLS-II and LCLS-II-HE, SHINE, SNS, ESS, FRIB, and EIC. To support and enable new accelerator-based applications and even make some of them feasible, we must continue addressing their challenges via a comprehensive RF R&D program that would advance the existing RF technologies and explore the nascent ones. △ Less

Submitted 25 August, 2022; originally announced August 2022.

Comments: contribution to Snowmass 2021

Next Generation Computational Tools for the Modeling and Design of Particle Accelerators at Exascale

Authors: Axel Huebl, Remi Lehe, Chad E. Mitchell, Ji Qiang, Robert D. Ryne, Ryan T. Sandberg, Jean-Luc Vay

Abstract: Particle accelerators are among the largest, most complex devices. To meet the challenges of increasing energy, intensity, accuracy, compactness, complexity and efficiency, increasingly sophisticated computational tools are required for their design and optimization. It is key that contemporary software take advantage of the latest advances in computer hardware and scientific software engineering… ▽ More Particle accelerators are among the largest, most complex devices. To meet the challenges of increasing energy, intensity, accuracy, compactness, complexity and efficiency, increasingly sophisticated computational tools are required for their design and optimization. It is key that contemporary software take advantage of the latest advances in computer hardware and scientific software engineering practices, delivering speed, reproducibility and feature composability for the aforementioned challenges. A new open source software stack is being developed at the heart of the Beam pLasma Accelerator Simulation Toolkit (BLAST) by LBNL and collaborators, providing new particle-in-cell modeling codes capable of exploiting the power of GPUs on Exascale supercomputers. Combined with advanced numerical techniques, such as mesh-refinement, and intrinsic support for machine learning, these codes are primed to provide ultrafast to ultraprecise modeling for future accelerator design and operations. △ Less

Submitted 9 August, 2022; v1 submitted 3 August, 2022; originally announced August 2022.

Comments: 4 pages, 8 figures; NAPAC22, Invited Oral, TUYE2

MSC Class: 78-10 ACM Class: I.6.0; D.2.12; D.2.13

Journal ref: NAPAC22, 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.

Comments: Contribution to Snowmass 2021, Accelerator Frontier

Snowmass21 Accelerator Modeling Community White Paper

Authors: S. Biedron, L. Brouwer, D. L. Bruhwiler, N. M. Cook, A. L. Edelen, D. Filippetto, C. -K. Huang, A. Huebl, T. Katsouleas, N. Kuklev, R. Lehe, S. Lund, C. Messe, W. Mori, C. -K. Ng, D. Perez, P. Piot, J. Qiang, R. Roussel, D. Sagan, A. Sahai, A. Scheinker, M. Thévenet, F. Tsung, J. -L. Vay , et al. (2 additional authors not shown)

Abstract: After a summary of relevant comments and recommendations from various reports over the last ten years, this paper examines the modeling needs in accelerator physics, from the modeling of single beams and individual accelerator elements, to the realization of virtual twins that replicate all the complexity to model a particle accelerator complex as accurately as possible. We then discuss cutting-ed… ▽ More After a summary of relevant comments and recommendations from various reports over the last ten years, this paper examines the modeling needs in accelerator physics, from the modeling of single beams and individual accelerator elements, to the realization of virtual twins that replicate all the complexity to model a particle accelerator complex as accurately as possible. We then discuss cutting-edge and emerging computing opportunities, such as advanced algorithms, AI/ML and quantum computing, computational needs in hardware, software performance, portability and scalability, and needs for scalable I/O and in-situ analysis. Considerations of reliability, long-term sustainability, user support and training are considered next, before discussing the benefits of ecosystems with integrated workflows based on standardized input and output, and with integrated frameworks and data repositories developed as a community. Last, we highlight how the community can work more collaboratively and efficiently through the development of consortia and centers, and via collaboration with industry. △ Less

Submitted 22 September, 2022; v1 submitted 15 March, 2022; originally announced March 2022.

Comments: contribution to Snowmass 2021

Accelerator and Beam Physics: Grand Challenges and Research Opportunities

Authors: S. Nagaitsev, V. Shiltsev, A. Valishev, T. Zolkin, J. -L. Vay, M. Bai, Y. Cai, M. J. Hogan, Z. Huang, J. Seeman, B. Dunham, X. Huang, T. Roser, M. Minty, J. Rosenzweig, P. Piot, J. Power, J. M. Byrd, A. Seryi, S. Lund, J. R. Patterson

Abstract: Accelerator and beam physics (ABP) is the science of the motion, generation, acceleration, manipulation, prediction, observation and use of charged particle beams. It focuses on fundamental long-term accelerator and beam physics research and development. Accelerator and beam physics research has resulted in important advances in accelerator science, yet support for this research is declining. NSF… ▽ More Accelerator and beam physics (ABP) is the science of the motion, generation, acceleration, manipulation, prediction, observation and use of charged particle beams. It focuses on fundamental long-term accelerator and beam physics research and development. Accelerator and beam physics research has resulted in important advances in accelerator science, yet support for this research is declining. NSF has terminated its program in accelerator Science and funding by DOE through GARD and Accelerator Stewardship has been steady or declining. The declining support for accelerator research will slow advances and threaten student training and work-force development in accelerator science. We propose a robust and scientifically challenging program in accelerator and beam physics, which will position the field of US High Energy Physics to be productive and competitive for decades to come. △ Less

Submitted 13 March, 2022; originally announced March 2022.

Comments: Contribution to Snowmass 2021. arXiv admin note: substantial text overlap with arXiv:2101.04107

Collective Effects and Intense Beam-Plasma Interactions in Ion-Beam-Driven High Energy Density Matter and Inertial Fusion Energy

Authors: Igor D. Kaganovich, Edward A. Startsev, Hong Qin, Erik Gilson, Thomas Schenkel, Jean-Luc Vay, Ed P. Lee, William Waldron, Roger Bangerter, Arun Persaud, Peter Seidl, Qing Ji, Alex Friedman, Dave P. Grote, John Barnard

Abstract: For the successful generation of ion-beam-driven high energy density matter and heavy ion fusion energy, intense ion beams must be transported and focused onto a target with small spot size. One of the successful approaches to achieve this goal is to accelerate and transport intense ion charge bunches in an accelerator and then focus the charge bunches ballistically in a section of the accelerator… ▽ More For the successful generation of ion-beam-driven high energy density matter and heavy ion fusion energy, intense ion beams must be transported and focused onto a target with small spot size. One of the successful approaches to achieve this goal is to accelerate and transport intense ion charge bunches in an accelerator and then focus the charge bunches ballistically in a section of the accelerator that contains a neutralizing background plasma. This requires the ability to control space-charge effects during un-neutralized (non-neutral) beam transport in the accelerator and transport sections, and the ability to effectively neutralize the space charge and current by propagating the beam through background plasma. As the beam intensity and energy are increased in future heavy ion fusion (HIF) drivers and Fast Ignition (FI) approaches, it is expected that nonlinear processes and collective effects will become much more pronounced than in previous experiments. Making use of 3D electromagnetic particle-in-cell simulation (PIC) codes (BEST, WARP-X, and LTP-PIC, etc.), the theory and modelling studies will be validated by comparing with experimental data on the 100kV Princeton Advanced Test Stand, and future experiments at the FAIR facility. The theoretical predictions that are developed will be scaled to the beam and plasma parameters relevant to heavy ion fusion drivers and Fast Ignition scenarios. Therefore, the theoretical results will also contribute significantly toward the long-term goal of fusion energy production by ion-beam-driven inertial confinement fusion. △ Less

Submitted 31 January, 2022; originally announced January 2022.

Comments: White paper submitted to the IFE Science & Technology Community Strategic Planning Workshop https://lasers.llnl.gov/nif-workshops/ife-workshop-2022/white-papers (2022)

Absorption of charged particles in Perfectly-Matched-Layers by optimal damping of the deposited current

Authors: Remi Lehe, Aurore Blelly, Lorenzo Giacomel, Revathi Jambunathan, Jean-Luc Vay

Abstract: Perfectly-Matched Layers (PML) are widely used in Particle-In-Cell simulations, in order to absorb electromagnetic waves that propagate out of the simulation domain. However, when charged particles cross the interface between the simulation domain and the PMLs, a number of numerical artifacts can arise. In order to mitigate these artifacts, we introduce a new PML algorithm whereby the current depo… ▽ More Perfectly-Matched Layers (PML) are widely used in Particle-In-Cell simulations, in order to absorb electromagnetic waves that propagate out of the simulation domain. However, when charged particles cross the interface between the simulation domain and the PMLs, a number of numerical artifacts can arise. In order to mitigate these artifacts, we introduce a new PML algorithm whereby the current deposited by the macroparticles in the PML is damped by an analytically-derived, optimal coefficient. The benefits of this new algorithm is illustrated in practical simulations. △ Less

Submitted 31 August, 2022; v1 submitted 22 January, 2022; originally announced January 2022.

Development of a Virtual EM Detector for the Advanced Particle Accelerator Modeling Code WarpX

Authors: Elisa Rheaume, Lorenzo Giacomel, Jean-Luc Vay, Axel Huebl

Abstract: In physics research particle accelerators are highly valued, and extraordinarily expensive, technical instruments. The high cost of particle accelerators results from the immense lengths required to accelerate particles to high energies, using radio frequency cavities. A current promising field of research, laser-driven particle acceleration has the potential to reduce the size, cost, and energy c… ▽ More In physics research particle accelerators are highly valued, and extraordinarily expensive, technical instruments. The high cost of particle accelerators results from the immense lengths required to accelerate particles to high energies, using radio frequency cavities. A current promising field of research, laser-driven particle acceleration has the potential to reduce the size, cost, and energy consumption of particle accelerators by orders of magnitude. To understand and control particle acceleration in plasmas using ultra-small spatial configurations, researchers have been developing computational models to simulate the acceleration environment. Within these models, computational scientists have introduced virtual diagnostics to serve as the digital parallel to experimental detectors. Using WarpX\cite{Vay:2021}, an advanced Particle-in-Cell code that simulates laser-driven particle acceleration, we have developed a virtual diagnostic to measure electromagnetic radiation. Such radiation can for instance be produced from scattered and transmitted laser beams. This \textit{FieldProbe} diagnostic measures individual field components for electric and magnetic fields at a point or at all points along a specified geometry (line or plane). This diagnostic is GPU-accelerated and parallelized using the Message Passing Interface (MPI) and can thus run on High Performance Computing Centers such as NERSC. △ Less

Submitted 24 October, 2022; v1 submitted 20 January, 2022; originally announced January 2022.

Comments: Poster based on this paper can be found at https://zenodo.org/record/5750895 (DOI:10.5281/zenodo.5750895) Data set used in the creation of this paper can be found at https://zenodo.org/record/5847120 (DOI:10.5281/zenodo.5847120)

PICSAR-QED: a Monte Carlo module to simulate Strong-Field Quantum Electrodynamics in Particle-In-Cell codes for exascale architectures

Authors: Luca Fedeli, Neïl Zaïm, Antonin Sainte-Marie, Maxence Thévenet, Axel Huebl, Andrew Myers, Jean-Luc Vay, Henri Vincenti

Abstract: Physical scenarios where the electromagnetic fields are so strong that Quantum ElectroDynamics (QED) plays a substantial role are one of the frontiers of contemporary plasma physics research. Investigating those scenarios requires state-of-the-art Particle-In-Cell (PIC) codes able to run on top high-performance computing machines and, at the same time, able to simulate strong-field QED processes.… ▽ More Physical scenarios where the electromagnetic fields are so strong that Quantum ElectroDynamics (QED) plays a substantial role are one of the frontiers of contemporary plasma physics research. Investigating those scenarios requires state-of-the-art Particle-In-Cell (PIC) codes able to run on top high-performance computing machines and, at the same time, able to simulate strong-field QED processes. This work presents the PICSAR-QED library, an open-source, portable implementation of a Monte Carlo module designed to provide modern PIC codes with the capability to simulate such processes, and optimized for high-performance computing. Detailed tests and benchmarks are carried out to validate the physical models in PICSAR-QED, to study how numerical parameters affect such models, and to demonstrate its capability to run on different architectures (CPUs and GPUs). Its integration with WarpX, a state-of-the-art PIC code designed to deliver scalable performance on upcoming exascale supercomputers, is also discussed and validated against results from the existing literature. △ Less

Submitted 1 October, 2021; originally announced October 2021.

Comments: 33 pages, 7 figures, 4 tables

HiPACE++: a portable, 3D quasi-static Particle-in-Cell code

Authors: Severin Diederichs, Carlo Benedetti, Axel Huebl, Rémi Lehe, Andrew Myers, Alexander Sinn, Jean-Luc Vay, Weiqun Zhang, Maxence Thévenet

Abstract: Modeling plasma accelerators is a computationally challenging task and the quasi-static particle-in-cell algorithm is a method of choice in a wide range of situations. In this work, we present the first performance-portable, quasi-static, three-dimensional particle-in-cell code HiPACE++. By decomposing all the computation of a 3D domain in successive 2D transverse operations and choosing appropria… ▽ More Modeling plasma accelerators is a computationally challenging task and the quasi-static particle-in-cell algorithm is a method of choice in a wide range of situations. In this work, we present the first performance-portable, quasi-static, three-dimensional particle-in-cell code HiPACE++. By decomposing all the computation of a 3D domain in successive 2D transverse operations and choosing appropriate memory management, HiPACE++ demonstrates orders-of-magnitude speedups on modern scientific GPUs over CPU-only implementations. The 2D transverse operations are performed on a single GPU, avoiding time-consuming communications. The longitudinal parallelization is done through temporal domain decomposition, enabling near-optimal strong scaling from 1 to 512 GPUs. HiPACE++ is a modular, open-source code enabling efficient modeling of plasma accelerators from laptops to state-of-the-art supercomputers. △ Less

Submitted 21 September, 2021; originally announced September 2021.

Modeling of Advanced Accelerator Concepts

Authors: J. -L. Vay, A. Huebl, R. Lehe, N. M. Cook, R. J. England, U. Niedermayer, P. Piot, F. Tsung, D. Winklehner

Abstract: Computer modeling is essential to research on Advanced Accelerator Concepts (AAC), as well as to their design and operation. This paper summarizes the current status and future needs of AAC systems and reports on several key aspects of (i) high-performance computing (including performance, portability, scalability, advanced algorithms, scalable I/Os and In-Situ analysis), (ii) the benefits of ecos… ▽ More Computer modeling is essential to research on Advanced Accelerator Concepts (AAC), as well as to their design and operation. This paper summarizes the current status and future needs of AAC systems and reports on several key aspects of (i) high-performance computing (including performance, portability, scalability, advanced algorithms, scalable I/Os and In-Situ analysis), (ii) the benefits of ecosystems with integrated workflows based on standardized input and output and with integrated frameworks developed as a community, and (iii) sustainability and reliability (including code robustness and usability). △ Less

Submitted 14 September, 2021; v1 submitted 10 September, 2021; originally announced September 2021.

Simulations of Future Particle Accelerators: Issues and Mitigations

Authors: D. Sagan, M. Berz, N. M. Cook, Y. Hao, G. Hoffstaetter, A. Huebl, C. -K. Huang, M. H. Langston, C. E. Mayes, C. E. Mitchell, C. -K. Ng, J. Qiang, R. D. Ryne, A. Scheinker, E. Stern, J. -L. Vay, D. Winklehner, H. Zhang

Abstract: The ever increasing demands placed upon machine performance have resulted in the need for more comprehensive particle accelerator modeling. Computer simulations are key to the success of particle accelerators. Many aspects of particle accelerators rely on computer modeling at some point, sometimes requiring complex simulation tools and massively parallel supercomputing. Examples include the modeli… ▽ More The ever increasing demands placed upon machine performance have resulted in the need for more comprehensive particle accelerator modeling. Computer simulations are key to the success of particle accelerators. Many aspects of particle accelerators rely on computer modeling at some point, sometimes requiring complex simulation tools and massively parallel supercomputing. Examples include the modeling of beams at extreme intensities and densities (toward the quantum degeneracy limit), and with ultra-fine control (down to the level of individual particles). In the future, adaptively tuned models might also be relied upon to provide beam measurements beyond the resolution of existing diagnostics. Much time and effort has been put into creating accelerator software tools, some of which are highly successful. However, there are also shortcomings such as the general inability of existing software to be easily modified to meet changing simulation needs. In this paper possible mitigating strategies are discussed for issues faced by the accelerator community as it endeavors to produce better and more comprehensive modeling tools. This includes lack of coordination between code developers, lack of standards to make codes portable and/or reusable, lack of documentation, among others. △ Less

Submitted 24 August, 2021; originally announced August 2021.

Comments: 21 pages, 1 figure. To be published in JINST

Improving I/O Performance for Exascale Applications through Online Data Layout Reorganization

Authors: Lipeng Wan, Axel Huebl, Junmin Gu, Franz Poeschel, Ana Gainaru, Ruonan Wang, Jieyang Chen, Xin Liang, Dmitry Ganyushin, Todd Munson, Ian Foster, Jean-Luc Vay, Norbert Podhorszki, Kesheng Wu, Scott Klasky

Abstract: The applications being developed within the U.S. Exascale Computing Project (ECP) to run on imminent Exascale computers will generate scientific results with unprecedented fidelity and record turn-around time. Many of these codes are based on particle-mesh methods and use advanced algorithms, especially dynamic load-balancing and mesh-refinement, to achieve high performance on Exascale machines. Y… ▽ More The applications being developed within the U.S. Exascale Computing Project (ECP) to run on imminent Exascale computers will generate scientific results with unprecedented fidelity and record turn-around time. Many of these codes are based on particle-mesh methods and use advanced algorithms, especially dynamic load-balancing and mesh-refinement, to achieve high performance on Exascale machines. Yet, as such algorithms improve parallel application efficiency, they raise new challenges for I/O logic due to their irregular and dynamic data distributions. Thus, while the enormous data rates of Exascale simulations already challenge existing file system write strategies, the need for efficient read and processing of generated data introduces additional constraints on the data layout strategies that can be used when writing data to secondary storage. We review these I/O challenges and introduce two online data layout reorganization approaches for achieving good tradeoffs between read and write performance. We demonstrate the benefits of using these two approaches for the ECP particle-in-cell simulation WarpX, which serves as a motif for a large class of important Exascale applications. We show that by understanding application I/O patterns and carefully designing data layouts we can increase read performance by more than 80%. △ Less

Submitted 15 July, 2021; originally announced July 2021.

Comments: 12 pages, 15 figures, accepted by IEEE Transactions on Parallel and Distributed Systems

Journal ref: IEEE Transactions on Parallel and Distributed Systems, 2021

A Hybrid Nodal-Staggered Pseudo-Spectral Electromagnetic Particle-In-Cell Method with Finite-Order Centering

Authors: Edoardo Zoni, Remi Lehe, Olga Shapoval, Daniel Belkin, Neil Zaïm, Luca Fedeli, Henri Vincenti, Jean-Luc Vay

Abstract: Electromagnetic particle-in-cell (PIC) codes are widely used to perform computer simulations of a variety of physical systems, including fusion plasmas, astrophysical plasmas, plasma wakefield particle accelerators, and secondary photon sources driven by ultra-intense lasers. In a PIC code, Maxwell's equations are solved on a grid with a numerical method of choice. This article focuses on pseudo-s… ▽ More Electromagnetic particle-in-cell (PIC) codes are widely used to perform computer simulations of a variety of physical systems, including fusion plasmas, astrophysical plasmas, plasma wakefield particle accelerators, and secondary photon sources driven by ultra-intense lasers. In a PIC code, Maxwell's equations are solved on a grid with a numerical method of choice. This article focuses on pseudo-spectral analytical time-domain (PSATD) algorithms and presents a novel hybrid PSATD PIC scheme that combines the respective advantages of standard nodal and staggered methods. The novelty of the hybrid scheme consists in using finite-order centering of grid quantities between nodal and staggered grids, in order to combine the solution of Maxwell's equations on a staggered grid with the deposition of charges and currents and the gathering of electromagnetic forces on a nodal grid. The correctness and performance of the novel hybrid scheme are assessed by means of numerical tests that employ different classes of PSATD equations in a variety of physical scenarios, ranging from the modeling of electron-positron pair creation in vacuum to the simulation of laser-driven and particle beam-driven plasma wakefield acceleration. It is shown that the novel hybrid scheme offers significant numerical and computational advantages, compared to purely nodal or staggered methods, for all the test cases presented. △ Less

Submitted 23 June, 2021; originally announced June 2021.

Comments: 39 pages, 15 figures, submitted to Computer Physics Communications

Overcoming timestep limitations in boosted-frame Particle-In-Cell simulations of plasma-based acceleration

Authors: Olga Shapoval, Remi Lehe, Maxence Thévenet, Edoardo Zoni, Yinjian Zhao, Jean-Luc Vay

Abstract: Explicit electromagnetic Particle-In-Cell (PIC) codes are typically limited by the Courant- Friedrichs-Lewy (CFL) condition, which implies that the timestep multiplied by the speed of light must be smaller than the smallest cell size. In the case of boosted-frame PIC simulations of plasma-based acceleration, this limitation can be a major hinderance as the cells are often very elongated along the… ▽ More Explicit electromagnetic Particle-In-Cell (PIC) codes are typically limited by the Courant- Friedrichs-Lewy (CFL) condition, which implies that the timestep multiplied by the speed of light must be smaller than the smallest cell size. In the case of boosted-frame PIC simulations of plasma-based acceleration, this limitation can be a major hinderance as the cells are often very elongated along the longitudinal direction and the timestep is thus limited by the small, transverse cell size. This entails many small-timestep PIC iterations, and can limit the potential speed-up of the boosted-frame technique. Here, by using a CFL-free analytical spectral solver, and by mitigating additional numerical instabilities that arise at large timestep, we show that it is possible to overcome traditional limitations on the timestep and thereby realize the full potential of the boosted-frame technique over a much wider range of parameters. △ Less

Submitted 28 April, 2021; originally announced April 2021.

In-Situ Assessment of Device-Side Compute Work for Dynamic Load Balancing in a GPU-Accelerated PIC Code

Authors: Michael E. Rowan, Axel Huebl, Kevin N. Gott, Jack Deslippe, Maxence Thévenet, Remi Lehe, Jean-Luc Vay

Abstract: Maintaining computational load balance is important to the performant behavior of codes which operate under a distributed computing model. This is especially true for GPU architectures, which can suffer from memory oversubscription if improperly load balanced. We present enhancements to traditional load balancing approaches and explicitly target GPU architectures, exploring the resulting performan… ▽ More Maintaining computational load balance is important to the performant behavior of codes which operate under a distributed computing model. This is especially true for GPU architectures, which can suffer from memory oversubscription if improperly load balanced. We present enhancements to traditional load balancing approaches and explicitly target GPU architectures, exploring the resulting performance. A key component of our enhancements is the introduction of several GPU-amenable strategies for assessing compute work. These strategies are implemented and benchmarked to find the most optimal data collection methodology for in-situ assessment of GPU compute work. For the fully kinetic particle-in-cell code WarpX, which supports MPI+CUDA parallelism, we investigate the performance of the improved dynamic load balancing via a strong scaling-based performance model and show that, for a laser-ion acceleration test problem run with up to 6144 GPUs on Summit, the enhanced dynamic load balancing achieves from 62%--74% (88% when running on 6 GPUs) of the theoretically predicted maximum speedup; for the 96-GPU case, we find that dynamic load balancing improves performance relative to baselines without load balancing (3.8x speedup) and with static load balancing (1.2x speedup). Our results provide important insights into dynamic load balancing and performance assessment, and are particularly relevant in the context of distributed memory applications ran on GPUs. △ Less

Submitted 22 April, 2021; originally announced April 2021.

Comments: 11 pages, 8 figures. Paper accepted in the Platform for Advanced Scientific Computing Conference (PASC '21), July 5 to 9, 2021, Geneva, Switzerland

Journal ref: PASC 2021: Proceedings of the Platform for Advanced Scientific Computing Conference

Porting WarpX to GPU-accelerated platforms

Authors: A. Myers, A. Almgren, L. D. Amorim, J. Bell, L. Fedeli, L. Ge, K. Gott, D. P. Grote, M. Hogan, A. Huebl, R. Jambunathan, R. Lehe, C. Ng, M. Rowan, O. Shapoval, M. Thévenet, J. -L. Vay, H. Vincenti, E. Yang, N. Zaïm, W. Zhang, Y. Zhao, E. Zoni

Abstract: WarpX is a general purpose electromagnetic particle-in-cell code that was originally designed to run on many-core CPU architectures. We describe the strategy followed to allow WarpX to use the GPU-accelerated nodes on OLCF's Summit supercomputer, a strategy we believe will extend to the upcoming machines Frontier and Aurora. We summarize the challenges encountered, lessons learned, and give curren… ▽ More WarpX is a general purpose electromagnetic particle-in-cell code that was originally designed to run on many-core CPU architectures. We describe the strategy followed to allow WarpX to use the GPU-accelerated nodes on OLCF's Summit supercomputer, a strategy we believe will extend to the upcoming machines Frontier and Aurora. We summarize the challenges encountered, lessons learned, and give current performance results on a series of relevant benchmark problems. △ Less

Submitted 2 September, 2021; v1 submitted 28 January, 2021; originally announced January 2021.

Comments: 11 pages, 5 figures, accepted by Parallel Computing. Minor revisions, results unchanged

Journal ref: Parallel Computing, Volume 108, 2021, 102833

Accelerator and Beam Physics Research Goals and Opportunities

Authors: S. Nagaitsev, Z. Huang, J. Power, J. -L. Vay, P. Piot, L. Spentzouris, J. Rosenzweig, Y. Cai, S. Cousineau, M. Conde, M. Hogan, A. Valishev, M. Minty, T. Zolkin, X. Huang, V. Shiltsev, J. Seeman, J. Byrd, Y. Hao, B. Dunham, B. Carlsten, A. Seryi, R. Patterson

Abstract: This report is a summary of two preparatory workshops, documenting the community vision for the national accelerator and beam physics research program. It identifies the Grand Challenges of accelerator and beam physics (ABP) field and documents research opportunities to address these Grand Challenges. This report will be used to develop a strategic research roadmap for the field of accelerator sci… ▽ More This report is a summary of two preparatory workshops, documenting the community vision for the national accelerator and beam physics research program. It identifies the Grand Challenges of accelerator and beam physics (ABP) field and documents research opportunities to address these Grand Challenges. This report will be used to develop a strategic research roadmap for the field of accelerator science. △ Less

Submitted 25 January, 2021; v1 submitted 11 January, 2021; originally announced January 2021.

Report number: FERMILAB-FN-1123

Probing strong-field QED with Doppler-boosted PetaWatt-class lasers

Authors: Luca Fedeli, Antonin Sainte-Marie, Neil Zaïm, Maxence Thévenet, Jean-Luc Vay, Andrew Myers, Fabien Quéré, Henri Vincenti

Abstract: We propose a scheme to explore regimes of strong-field Quantum Electrodynamics (SF-QED) otherwise unattainable with the currently available laser technology. The scheme relies on relativistic plasma mirrors curved by radiation pressure to boost the intensity of PetaWatt-class laser pulses by Doppler effect and focus them to extreme field intensities. We show that very clear SF-QED signatures could… ▽ More We propose a scheme to explore regimes of strong-field Quantum Electrodynamics (SF-QED) otherwise unattainable with the currently available laser technology. The scheme relies on relativistic plasma mirrors curved by radiation pressure to boost the intensity of PetaWatt-class laser pulses by Doppler effect and focus them to extreme field intensities. We show that very clear SF-QED signatures could be observed by placing a secondary target where the boosted beam is focused. △ Less

Submitted 14 December, 2020; originally announced December 2020.

Comments: 5 figures

Journal ref: Phys. Rev. Lett. 127, 114801 (2021)

Simulation of plasma accelerators with the Particle-In-Cell method

Authors: J. L. Vay

Abstract: We present the standard electromagnetic Particle-in-Cell method, starting from the discrete approximation of derivatives on a uniform grid. The application to second-order, centered, finite-difference discretization of the equations of motion and of Maxwells equations is then described in one dimension, followed by two and three dimensions. Various algorithms are presented, for which we discuss th… ▽ More We present the standard electromagnetic Particle-in-Cell method, starting from the discrete approximation of derivatives on a uniform grid. The application to second-order, centered, finite-difference discretization of the equations of motion and of Maxwells equations is then described in one dimension, followed by two and three dimensions. Various algorithms are presented, for which we discuss the stability and accuracy, introducing and elucidating concepts like numerical stochastic heating, CFL limit and numerical dispersion. The coupling of the particles and field quantities via interpolation at various orders is detailed, together with its implication on energy and momentum conserving. Special topics of relevance to the modeling of plasma accelerators are discussed, such as moving window, optimal Lorentz boosted frame, the numerical Cherenkov instability and its mitigation. Examples of simulations of laser-driven and particle beam-driven accelerators are given, including with mesh refinement. We conclude with a discussion on high-performance computing and a brief outlook. △ Less

Submitted 17 August, 2020; originally announced August 2020.

Comments: 39 pages, 24 figures

Report number: Proceedings of the 2019 CERN--Accelerator--School course on High Gradient Wakefield Accelerators, Sesimbra, (Portugal)

Intra-bunch feedback system developments at DAFNE

Authors: Alessandro Drago, INFN-LNF, Frascati, Tor Vergata University, Rome, Italy D. Alesini, S. Caschera, A. Gallo, INFN-LNF, Frascati, Italy J. D. Fox, Stanford University, Stanford, USA J. Cesaratto, J. Dusatko, J. Olsen, C. Rivetta, O. Turgut, SLAC, Menlo Park, USA W. Hofle, G. Iadarola, K. Li, E. Metral, E. Montesinos , et al. (12 additional authors not shown)

Abstract: This paper presents history and evolution of the intra-bunch feedback system for circular accelerators. This pro-ject has been presented by John D. Fox (SLAC/Stanford Un.) at the IPAC2010 held in Kyoto. The idea of the pro-posal is to build a flexible and powerful instrument to mit-igate the parasitic e-cloud effects on the proton (and poten-tially positron) beams in storage rings. Being a new and… ▽ More This paper presents history and evolution of the intra-bunch feedback system for circular accelerators. This pro-ject has been presented by John D. Fox (SLAC/Stanford Un.) at the IPAC2010 held in Kyoto. The idea of the pro-posal is to build a flexible and powerful instrument to mit-igate the parasitic e-cloud effects on the proton (and poten-tially positron) beams in storage rings. Being a new and ambitious project, the financial issues have been quite im-portant. US LHC Accelerator Research Program (LARP) and other institution funding sources have assured the de-velopment of the design for implementing the feedback in the SPS ring at CERN. Here the intra-bunch feedback sys-tem has been installed and tested in the frame of the LIU (LHC Injector Upgrade) program. After the end of the LARP funding, a possible new inter-esting chance to continue the R&D activity, could be by implementing the system in a lepton storage ring affected by e-cloud effects. For achieving this goal, a possible ex-periment could be carried out in the positron ring of DAFNE at Frascati, Italy. The feasibility of the proposal is evaluated in the following sections. In case of approval of the experiment, indeed the project could be inserted in the DAFNE-TF (DAFNE Test Facility) program that is fore-seen after the 2020 for the following 3-5 years. △ Less

Submitted 24 April, 2020; originally announced April 2020.

Comments: 6 pages, 5 figures, Sixth e-cloud workshop, ECLOUD'18 , 3-7 June, 2018, La Biodola (Isola d'Elba ), Italy

Investigation of light ion fusion reactions with plasma discharges

Authors: T. Schenkel, A. Persaud, H. Wang, P. A. Seidl, R. MacFadyen, C. Nelson, W. L. Waldron, J. -L. Vay, G. Deblonde, B. Wen, Y. -M. Chiang, B. P. MacLeod, Q. Ji

Abstract: The scaling of reaction yields in light ion fusion to low reaction energies is important for our understanding of stellar fuel chains and the development of future energy technologies. Experiments become progressively more challenging at lower reaction energies due to the exponential drop of fusion cross sections below the Coulomb barrier. We report on experiments where deuterium-deuterium (D-D) f… ▽ More The scaling of reaction yields in light ion fusion to low reaction energies is important for our understanding of stellar fuel chains and the development of future energy technologies. Experiments become progressively more challenging at lower reaction energies due to the exponential drop of fusion cross sections below the Coulomb barrier. We report on experiments where deuterium-deuterium (D-D) fusion reactions are studied in a pulsed plasma in the glow discharge regime using a benchtop apparatus. We model plasma conditions using particle-in-cell codes. Advantages of this approach are relatively high peak ion currents and current densities (0.1 to several A/cm^2) that can be applied to metal wire cathodes for several days. We detect neutrons from D-D reactions with scintillator-based detectors. For palladium targets, we find neutron yields as a function of cathode voltage that are over 100 times higher than yields expected for bare nuclei fusion at ion energies below 2 keV (center of mass frame). A possible explanation is a correction to the ion energy due to an electron screening potential of 1000+/-250 eV, which increases the probability for tunneling through the repulsive Coulomb barrier. Our compact, robust setup enables parametric studies of this effect at relatively low reaction energies. △ Less

Submitted 24 May, 2019; v1 submitted 8 May, 2019; originally announced May 2019.

Journal ref: Journal of Applied Physics 126, 203302 (2019)

Ion Acceleration in Laser Generated Mega Tesla Magnetic Vortex

Authors: Jaehong Park, Stepan S. Bulanov, Jianhui Bin, Qing Ji, Sven Steinke, Jean-Luc Vay, Cameron G. R. Geddes, Carl B. Schroeder, Wim P. Leemans, Thomas Schenkel, Eric Esarey

Abstract: Magnetic Vortex Acceleration (MVA) from near critical density targets is one of the promising schemes of laser-driven ion acceleration. 3D particle-in-cell simulations are used to explore a more extensive laser-target parameter space than previously reported on in the literature as well as to study the laser pulse coupling to the target, the structure of the fields, and the properties of the accel… ▽ More Magnetic Vortex Acceleration (MVA) from near critical density targets is one of the promising schemes of laser-driven ion acceleration. 3D particle-in-cell simulations are used to explore a more extensive laser-target parameter space than previously reported on in the literature as well as to study the laser pulse coupling to the target, the structure of the fields, and the properties of the accelerated ion beam in the MVA scheme. The efficiency of acceleration depends on the coupling of the laser energy to the self-generated channel in the target. The accelerated proton beams demonstrate high level of collimation with achromatic angular divergence, and carry a significant amount of charge. For PW-class lasers, this acceleration regime provides favorable scaling of maximum ion energy with laser power for optimized interaction parameters. The mega Tesla-level magnetic fields generated by the laser-driven co-axial plasma structure in the target are prerequisite for accelerating protons to the energy of several hundred MeV. △ Less

Submitted 5 April, 2019; originally announced April 2019.

Comments: submitted to Physics of Plasmas

A generalized massively parallel ultra-high order FFT-based Maxwell solver

Authors: Haithem Kallala, Jean-Luc Vay, Henri Vincenti

Abstract: Dispersion-free ultra-high order FFT-based Maxwell solvers have recently proven to be paramount to a large range of applications, including the high-fidelity modeling of high-intensity laser-matter interactions with Particle-In-Cell (PIC) codes. To enable a massively parallel scaling of these solvers, a novel parallelization technique was recently proposed, which consists in splitting the simulati… ▽ More Dispersion-free ultra-high order FFT-based Maxwell solvers have recently proven to be paramount to a large range of applications, including the high-fidelity modeling of high-intensity laser-matter interactions with Particle-In-Cell (PIC) codes. To enable a massively parallel scaling of these solvers, a novel parallelization technique was recently proposed, which consists in splitting the simulation domain into several processor sub-domains, with guard regions appended at each sub-domain boundaries. Maxwell's equations are advanced independently on each sub-domain using local shared-memory FFTs (instead of a single distributed global FFT). This implies small truncation errors at sub-domain boundaries, the amplitude of which depends on guard regions sizes and order of the Maxwell solver. For a moderate number of guard cells (<10) , this 'local' technique proved to be highly scalable on up to a million cores. Yet, depending on the targeted applications, the number of guard cells required to mitigate truncations errors might be very large, which would severely limit the parallel efficiency of this technique due to the large volume of guard cells to be exchanged between sub-domains. In this context, we propose a novel 'hybrid' parallelization technique that ensures very good scaling of FFT-based solvers with an arbitrarily high number of guard cells. It consists in performing distributed FFTs on local groups of processors with guard regions now appended to boundaries of each group of processors. A dual domain decomposition method is used for the Maxwell solver and other parts of the PIC cycle to keep the simulation load-balanced. This 'hybrid' technique was implemented in the open source exascale library PICSAR. Benchmarks show that for a large number of guard cells ($>16$), the 'hybrid' technique offers a $\times 3$ speed-up and $\times 8$ memory savings compared to the 'local' one. △ Less

Submitted 18 December, 2018; originally announced December 2018.

Comments: 11 pages, 9 figures

Identification of coupling mechanisms between ultraintense laser light and dense plasmas

Authors: L. Chopineau, A. Leblanc, G. Blaclard, A. Denoeud, M. Thévenet, J-L. Vay, G. Bonnaud, Ph. Martin, H. Vincenti, F. Quéré

Abstract: The interaction of intense laser beams with plasmas created on solid targets involves a rich non-linear physics. Because such dense plasmas are reflective for laser light, the coupling with the incident beam occurs within a thin layer at the interface between plasma and vacuum. One of the main paradigms used to understand this coupling, known as Brunel mechanism, is expected to be valid only for v… ▽ More The interaction of intense laser beams with plasmas created on solid targets involves a rich non-linear physics. Because such dense plasmas are reflective for laser light, the coupling with the incident beam occurs within a thin layer at the interface between plasma and vacuum. One of the main paradigms used to understand this coupling, known as Brunel mechanism, is expected to be valid only for very steep plasma surfaces. Despite innumerable studies, its validity range remains uncertain, and the physics involved for smoother plasma-vacuum interfaces is unclear, especially for ultrahigh laser intensities. We report the first comprehensive experimental and numerical study of the laser-plasma coupling mechanisms as a function of the plasma interface steepness, in the relativistic interaction regime. Our results reveal a clear transition from the temporally-periodic Brunel mechanism to a chaotic dynamic associated to stochastic heating. By revealing the key signatures of these two distinct regimes on experimental observables, we provide an important landmark for the interpretation of future experiments. △ Less

Submitted 11 September, 2018; originally announced September 2018.

Comments: 16 pages, 12 figures

Journal ref: Phys. Rev. X 9, 011050 (2019)

Convergence in nonlinear laser wakefield accelerators modeling in a Lorentz-boosted frame

Authors: Patrick Lee, Jean-Luc Vay

Abstract: Laser wakefield acceleration modeling using the Lorentz-boosted frame technique in the particle-in-cell code has demonstrated orders of magnitude speedups. A convergence study was previously conducted in cases with external injection in the linear regime and without injection in the nonlinear regime, and the obtained results have shown a convergence within the percentage level. In this article, a… ▽ More Laser wakefield acceleration modeling using the Lorentz-boosted frame technique in the particle-in-cell code has demonstrated orders of magnitude speedups. A convergence study was previously conducted in cases with external injection in the linear regime and without injection in the nonlinear regime, and the obtained results have shown a convergence within the percentage level. In this article, a convergence study is carried out to model electron self-injection in the 2-1/2D configuration. It is observed that the Lorentz-boosted frame technique is capable of modeling complex particle dynamics with a significant speedup. This result is crucial to curtail the computational time of the modeling of future chains of $10\,\mathrm{GeV}$ laser wakefield accelerator stages with high accuracy. △ Less

Submitted 5 March, 2018; originally announced March 2018.

Comments: 11 pages, 14 figures

Investigation of the dynamics of ionization induced injected electrons under the influence of beam loading effects

Authors: P. Lee, T. L. Audet, R. Lehe, J. -L. Vay, G. Maynard, B. Cros

Abstract: In laser-driven wakefield, ionization induced injection is an efficient way to inject electrons in the plasma wave. A detailed study on the beam dynamics under the influence of beam loading effects, which can be controlled by the concentration of nitrogen impurity introduced in the hydrogen gas was conducted. For a specific value of this percentage, the final energy of the high-energy electron bun… ▽ More In laser-driven wakefield, ionization induced injection is an efficient way to inject electrons in the plasma wave. A detailed study on the beam dynamics under the influence of beam loading effects, which can be controlled by the concentration of nitrogen impurity introduced in the hydrogen gas was conducted. For a specific value of this percentage, the final energy of the high-energy electron bunch becomes nearly independent of the trapped positions, thus leading to a small energy dispersion. We also show that the final beam emittance is mainly determined by the injection process. △ Less

Submitted 1 February, 2018; originally announced February 2018.

Warp-X: a new exascale computing platform for beam-plasma simulations

Authors: J. -L. Vay, A. Almgren, J. Bell, L. Ge, D. P. Grote, M. Hogan, O. Kononenko, R. Lehe, A. Myers, C. Ng, J. Park, R. Ryne, O. Shapoval, M. Thevenet, W. Zhang

Abstract: Turning the current experimental plasma accelerator state-of-the-art from a promising technology into mainstream scientific tools depends critically on high-performance, high-fidelity modeling of complex processes that develop over a wide range of space and time scales. As part of the U.S. Department of Energy's Exascale Computing Project, a team from Lawrence Berkeley National Laboratory, in coll… ▽ More Turning the current experimental plasma accelerator state-of-the-art from a promising technology into mainstream scientific tools depends critically on high-performance, high-fidelity modeling of complex processes that develop over a wide range of space and time scales. As part of the U.S. Department of Energy's Exascale Computing Project, a team from Lawrence Berkeley National Laboratory, in collaboration with teams from SLAC National Accelerator Laboratory and Lawrence Livermore National Laboratory, is developing a new plasma accelerator simulation tool that will harness the power of future exascale supercomputers for high-performance modeling of plasma accelerators. We present the various components of the codes such as the new Particle-In-Cell Scalable Application Resource (PICSAR) and the redesigned adaptive mesh refinement library AMReX, which are combined with redesigned elements of the Warp code, in the new WarpX software. The code structure, status, early examples of applications and plans are discussed. △ Less

Submitted 8 January, 2018; originally announced January 2018.

Saturation of the hosing instability in quasi-linear plasma accelerators

Authors: Remi Lehe, Carl B. Schroeder, Jean-Luc Vay, Eric Esarey, Wim P. Leemans

Abstract: The beam hosing instability is analyzed theoretically for a witness beam in the quasi-linear regime of plasma accelerators. In this regime, the hosing instability saturates, even for a monoenergetic bunch, at a level much less than standard scalings predict. Analytic expressions are derived for the saturation distance and amplitude and are in agreement with numerical results. Saturation is due to… ▽ More The beam hosing instability is analyzed theoretically for a witness beam in the quasi-linear regime of plasma accelerators. In this regime, the hosing instability saturates, even for a monoenergetic bunch, at a level much less than standard scalings predict. Analytic expressions are derived for the saturation distance and amplitude and are in agreement with numerical results. Saturation is due to the natural head-to-tail variations in the focusing force, including the self-consistent transverse beam loading. △ Less

Submitted 3 January, 2018; originally announced January 2018.

Journal ref: 119, PRL 119, 244801 (2017)

Optimization of laser-plasma injector via beam loading effects using ionization-induced injection

Authors: P. Lee, G. Maynard, T. L. Audet, R. Lehe, J. L. Vay, B. Cros

Abstract: Simulations of ionization induced injection in a laser driven plasma wakefield show that high-quality electron injectors in the 50-200 MeV range can be achieved in a gas cell with a tailored density profile. Using the PIC code Warp with parameters close to existing experimental conditions, we show that the concentration of $\mathrm{N_2}$ in a hydrogen plasma with a tailored density profile is an e… ▽ More Simulations of ionization induced injection in a laser driven plasma wakefield show that high-quality electron injectors in the 50-200 MeV range can be achieved in a gas cell with a tailored density profile. Using the PIC code Warp with parameters close to existing experimental conditions, we show that the concentration of $\mathrm{N_2}$ in a hydrogen plasma with a tailored density profile is an efficient parameter to tune electron beam properties through the control of the interplay between beam loading effects and varying accelerating field in the density profile. For a given laser plasma configuration, with moderate normalized laser amplitude, $a_0=1.6$ and maximum electron plasma density, $n_{e0}=4\times 10^{18}\,\mathrm{cm^{-3}}$, the optimum concentration results in a robust configuration to generate electrons at 150~MeV with a rms energy spread of 4\% and a spectral charge density of 1.8~pC/MeV. △ Less

Submitted 5 November, 2017; originally announced November 2017.

Comments: 13 pages, 10 figures

Ultrahigh-order Maxwell solver with extreme scalability for electromagnetic PIC simulations of plasmas

Authors: Henri Vincenti, Jean-Luc Vay

Abstract: The advent of massively parallel supercomputers, with their distributed-memory technology using many processing units, has favored the development of highly-scalable local low-order solvers at the expense of harder-to-scale global very high-order spectral methods. Indeed, FFT-based methods, which were very popular on shared memory computers, have been largely replaced by finite-difference (FD) met… ▽ More The advent of massively parallel supercomputers, with their distributed-memory technology using many processing units, has favored the development of highly-scalable local low-order solvers at the expense of harder-to-scale global very high-order spectral methods. Indeed, FFT-based methods, which were very popular on shared memory computers, have been largely replaced by finite-difference (FD) methods for the solution of many problems, including plasmas simulations with electromagnetic Particle-In-Cell methods. For some problems, such as the modeling of so-called "plasma mirrors" for the generation of high-energy particles and ultra-short radiations, we have shown that the inaccuracies of standard FD-based PIC methods prevent the modeling on present supercomputers at sufficient accuracy. We demonstrate here that a new method, based on the use of local FFTs, enables ultrahigh-order accuracy with unprecedented scalability, and thus for the first time the accurate modeling of plasma mirrors in 3D. △ Less

Submitted 26 July, 2017; originally announced July 2017.

Comments: 9 pages, 6 figures

Accurate modeling of plasma acceleration with arbitrary order pseudo-spectral particle-in-cell methods

Authors: Sören Jalas, Irene Dornmair, Rémi Lehe, Henri Vincenti, Jean-Luc Vay, Manuel Kirchen, Andreas R. Maier

Abstract: Particle in Cell (PIC) simulations are a widely used tool for the investigation of both laser- and beam-driven plasma acceleration. It is a known issue that the beam quality can be artificially degraded by numerical Cherenkov radiation (NCR) resulting primarily from an incorrectly modeled dispersion relation. Pseudo-spectral solvers featuring infinite order stencils can strongly reduce NCR, or eve… ▽ More Particle in Cell (PIC) simulations are a widely used tool for the investigation of both laser- and beam-driven plasma acceleration. It is a known issue that the beam quality can be artificially degraded by numerical Cherenkov radiation (NCR) resulting primarily from an incorrectly modeled dispersion relation. Pseudo-spectral solvers featuring infinite order stencils can strongly reduce NCR, or even suppress it, and are therefore well suited to correctly model the beam properties. For efficient parallelization of the PIC algorithm, however, localized solvers are inevitable. Arbitrary order pseudo-spectral methods provide this needed locality. Yet, these methods can again be prone to NCR. Here, we show that acceptably low solver orders are sufficient to correctly model the physics of interest, while allowing for parallel computation by domain decomposition. △ Less

Submitted 29 March, 2017; v1 submitted 17 November, 2016; originally announced November 2016.

Journal ref: Physics of Plasmas 24, 033115 (2017)

Pseudo-spectral Maxwell solvers for an accurate modeling of Doppler harmonic generation on plasma mirrors with Particle-In-Cell codes

Authors: G. Blaclard, H. Vincenti, R. Lehe, J. L. Vay

Abstract: With the advent of PW class lasers, the very large laser intensities attainable on-target should enable the production of intense high order Doppler harmonics from relativistic laser-plasma mirrors interactions. At present, the modeling of these harmonics with Particle-In-Cell (PIC) codes is extremely challenging as it implies an accurate description of tens of harmonic orders on a a broad range o… ▽ More With the advent of PW class lasers, the very large laser intensities attainable on-target should enable the production of intense high order Doppler harmonics from relativistic laser-plasma mirrors interactions. At present, the modeling of these harmonics with Particle-In-Cell (PIC) codes is extremely challenging as it implies an accurate description of tens of harmonic orders on a a broad range of angles. In particular, we show here that standard Finite Difference Time Domain (FDTD) Maxwell solvers used in most PIC codes partly fail to model Doppler harmonic generation because they induce numerical dispersion of electromagnetic waves in vacuum which is responsible for a spurious angular deviation of harmonic beams. This effect was extensively studied and a simple toy-model based on Snell-Descartes law was developed that allows us to finely predict the angular deviation of harmonics depending on the spatio-temporal resolution and the Maxwell solver used in the simulations. Our model demonstrates that the mitigation of this numerical artifact with FDTD solvers mandates very high spatio-temporal resolution preventing doing realistic 3D simulations. We finally show that non dispersive pseudo-spectral analytical time domain solvers can considerably reduce the spatio-temporal resolution required to mitigate this spurious deviation and should enable in the the near future 3D accurate modeling on supercomputers in a realistic time-to-solution. △ Less

Submitted 1 August, 2017; v1 submitted 19 August, 2016; originally announced August 2016.

Comments: 12 pages, 7 figures, 2 tables

Elimination of Numerical Cherenkov Instability in flowing-plasma Particle-In-Cell simulations by using Galilean coordinates

Authors: Remi Lehe, Manuel Kirchen, Brendan B. Godfrey, Andreas R. Maier, Jean-Luc Vay

Abstract: Particle-In-Cell (PIC) simulations of relativistic flowing plasmas are of key interest to several fields of physics (including e.g. laser-wakefield acceleration, when viewed in a Lorentz-boosted frame), but remain sometimes infeasible due to the well-known numerical Cherenkov instability (NCI). In this article, we show that, for a plasma drifting at a uniform relativistic velocity, the NCI can be… ▽ More Particle-In-Cell (PIC) simulations of relativistic flowing plasmas are of key interest to several fields of physics (including e.g. laser-wakefield acceleration, when viewed in a Lorentz-boosted frame), but remain sometimes infeasible due to the well-known numerical Cherenkov instability (NCI). In this article, we show that, for a plasma drifting at a uniform relativistic velocity, the NCI can be eliminated by simply integrating the PIC equations in Galilean coordinates that follow the plasma (also sometimes known as comoving coordinates) within a spectral analytical framework. The elimination of the NCI is verified empirically and confirmed by a theoretical analysis of the instability. Moreover, it is shown that this method is applicable both to Cartesian geometry and to cylindrical geometry with azimuthal Fourier decomposition. △ Less

Submitted 31 July, 2016; originally announced August 2016.

Comments: 18 pages, 6 figures

Journal ref: Phys. Rev. E 94, 053305 (2016)

Stable discrete representation of relativistically drifting plasmas

Authors: Manuel Kirchen, Remi Lehe, Brendan B. Godfrey, Irene Dornmair, Soeren Jalas, Kevin Peters, Jean-Luc Vay, Andreas R. Maier

Abstract: Representing the electrodynamics of relativistically drifting particle ensembles in discrete, co-propagating Galilean coordinates enables the derivation of a Particle-in-Cell algorithm that is intrinsically free of the Numerical Cherenkov Instability, for plasmas flowing at a uniform velocity. Application of the method is shown by modeling plasma accelerators in a Lorentz-transformed optimal frame… ▽ More Representing the electrodynamics of relativistically drifting particle ensembles in discrete, co-propagating Galilean coordinates enables the derivation of a Particle-in-Cell algorithm that is intrinsically free of the Numerical Cherenkov Instability, for plasmas flowing at a uniform velocity. Application of the method is shown by modeling plasma accelerators in a Lorentz-transformed optimal frame of reference. △ Less

Submitted 31 July, 2016; originally announced August 2016.

ASCR/HEP Exascale Requirements Review Report

Authors: Salman Habib, Robert Roser, Richard Gerber, Katie Antypas, Katherine Riley, Tim Williams, Jack Wells, Tjerk Straatsma, A. Almgren, J. Amundson, S. Bailey, D. Bard, K. Bloom, B. Bockelman, A. Borgland, J. Borrill, R. Boughezal, R. Brower, B. Cowan, H. Finkel, N. Frontiere, S. Fuess, L. Ge, N. Gnedin, S. Gottlieb , et al. (29 additional authors not shown)

Abstract: This draft report summarizes and details the findings, results, and recommendations derived from the ASCR/HEP Exascale Requirements Review meeting held in June, 2015. The main conclusions are as follows. 1) Larger, more capable computing and data facilities are needed to support HEP science goals in all three frontiers: Energy, Intensity, and Cosmic. The expected scale of the demand at the 2025 ti… ▽ More This draft report summarizes and details the findings, results, and recommendations derived from the ASCR/HEP Exascale Requirements Review meeting held in June, 2015. The main conclusions are as follows. 1) Larger, more capable computing and data facilities are needed to support HEP science goals in all three frontiers: Energy, Intensity, and Cosmic. The expected scale of the demand at the 2025 timescale is at least two orders of magnitude -- and in some cases greater -- than that available currently. 2) The growth rate of data produced by simulations is overwhelming the current ability, of both facilities and researchers, to store and analyze it. Additional resources and new techniques for data analysis are urgently needed. 3) Data rates and volumes from HEP experimental facilities are also straining the ability to store and analyze large and complex data volumes. Appropriately configured leadership-class facilities can play a transformational role in enabling scientific discovery from these datasets. 4) A close integration of HPC simulation and data analysis will aid greatly in interpreting results from HEP experiments. Such an integration will minimize data movement and facilitate interdependent workflows. 5) Long-range planning between HEP and ASCR will be required to meet HEP's research needs. To best use ASCR HPC resources the experimental HEP program needs a) an established long-term plan for access to ASCR computational and data resources, b) an ability to map workflows onto HPC resources, c) the ability for ASCR facilities to accommodate workflows run by collaborations that can have thousands of individual members, d) to transition codes to the next-generation HPC platforms that will be available at ASCR facilities, e) to build up and train a workforce capable of developing and using simulations and analysis to support HEP scientific research on next-generation systems. △ Less

Submitted 31 March, 2016; v1 submitted 30 March, 2016; originally announced March 2016.

Comments: 77 pages, 13 Figures; draft report, subject to further revision

An efficient and portable SIMD algorithm for charge/current deposition in Particle-In-Cell codes

Authors: H. Vincenti, R. Lehe, R. Sasanka, J-L. Vay

Abstract: In current computer architectures, data movement (from die to network) is by far the most energy consuming part of an algorithm (10pJ/word on-die to 10,000pJ/word on the network). To increase memory locality at the hardware level and reduce energy consumption related to data movement, future exascale computers tend to use more and more cores on each compute nodes ("fat nodes") that will have a red… ▽ More In current computer architectures, data movement (from die to network) is by far the most energy consuming part of an algorithm (10pJ/word on-die to 10,000pJ/word on the network). To increase memory locality at the hardware level and reduce energy consumption related to data movement, future exascale computers tend to use more and more cores on each compute nodes ("fat nodes") that will have a reduced clock speed to allow for efficient cooling. To compensate for frequency decrease, machine vendors are making use of long SIMD instruction registers that are able to process multiple data with one arithmetic operator in one clock cycle. SIMD register length is expected to double every four years. As a consequence, Particle-In-Cell (PIC) codes will have to achieve good vectorization to fully take advantage of these upcoming architectures. In this paper, we present a new algorithm that allows for efficient and portable SIMD vectorization of current/charge deposition routines that are, along with the field gathering routines, among the most time consuming parts of the PIC algorithm. Our new algorithm uses a particular data structure that takes into account memory alignement constraints and avoids gather/scatter instructions that can significantly affect vectorization performances on current CPUs. The new algorithm was successfully implemented in the 3D skeleton PIC code PICSAR and tested on Haswell Xeon processors (AVX2-256 bits wide data registers). Results show a factor of $\times 2$ to $\times 2.5$ speed-up in double precision for particle shape factor of order $1$ to $3$. The new algorithm can be applied as is on future KNL (Knights Landing) architectures that will include AVX-512 instruction sets with 512 bits register lengths (8 doubles/16 singles). △ Less

Submitted 8 January, 2016; originally announced January 2016.

Comments: 36 pages, 5 figures