-
Dephasing of ion beams as Magnetic Vortex Acceleration regime transitions into a bubble-like field structure
Authors:
Sahel Hakimi,
Stepan S. Bulanov,
Axel Huebl,
Lieselotte Obst-Huebl,
Kei Nakamura,
Anthony Gonsalves,
Thomas Schenkel,
Jeroen van Tilborg,
Jean-Luc Vay,
Carl B. Schroeder,
Eric Esarey,
Cameron R. Geddes
Abstract:
The interaction of an ultra-intense laser pulse with a near critical density target results in the formation of a plasma channel, a strong azimuthal magnetic field and moving vortices. An application of this is the generation of energetic and collimated ion beams via Magnetic Vortex Acceleration. The optimized regime of Magnetic Vortex Acceleration is becoming experimentally accessible with new hi…
▽ More
The interaction of an ultra-intense laser pulse with a near critical density target results in the formation of a plasma channel, a strong azimuthal magnetic field and moving vortices. An application of this is the generation of energetic and collimated ion beams via Magnetic Vortex Acceleration. The optimized regime of Magnetic Vortex Acceleration is becoming experimentally accessible with new high intensity laser beamlines coming online and advances made in near critical density target fabrication. The robustness of the acceleration mechanism with realistic experimental conditions is examined with three-dimensional simulations. Of particular interest is the acceleration performance with different laser temporal contrast conditions, in some cases leading to pre-expanded target profiles prior to the arrival of the main pulse. Preplasma effects on the structure of the accelerating fields is explored, including a detailed analysis of the ion beam properties and the efficiency of the process. Improved scaling laws for the MVA mechanism, including the laser focal spot size effects, are presented.
△ Less
Submitted 13 September, 2024;
originally announced September 2024.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
High Energy Physics Forum for Computational Excellence: Working Group Reports (I. Applications Software II. Software Libraries and Tools III. Systems)
Authors:
Salman Habib,
Robert Roser,
Tom LeCompte,
Zach Marshall,
Anders Borgland,
Brett Viren,
Peter Nugent,
Makoto Asai,
Lothar Bauerdick,
Hal Finkel,
Steve Gottlieb,
Stefan Hoeche,
Paul Sheldon,
Jean-Luc Vay,
Peter Elmer,
Michael Kirby,
Simon Patton,
Maxim Potekhin,
Brian Yanny,
Paolo Calafiura,
Eli Dart,
Oliver Gutsche,
Taku Izubuchi,
Adam Lyon,
Don Petravick
Abstract:
Computing plays an essential role in all aspects of high energy physics. As computational technology evolves rapidly in new directions, and data throughput and volume continue to follow a steep trend-line, it is important for the HEP community to develop an effective response to a series of expected challenges. In order to help shape the desired response, the HEP Forum for Computational Excellence…
▽ More
Computing plays an essential role in all aspects of high energy physics. As computational technology evolves rapidly in new directions, and data throughput and volume continue to follow a steep trend-line, it is important for the HEP community to develop an effective response to a series of expected challenges. In order to help shape the desired response, the HEP Forum for Computational Excellence (HEP-FCE) initiated a roadmap planning activity with two key overlapping drivers -- 1) software effectiveness, and 2) infrastructure and expertise advancement. The HEP-FCE formed three working groups, 1) Applications Software, 2) Software Libraries and Tools, and 3) Systems (including systems software), to provide an overview of the current status of HEP computing and to present findings and opportunities for the desired HEP computational roadmap. The final versions of the reports are combined in this document, and are presented along with introductory material.
△ Less
Submitted 28 October, 2015;
originally announced October 2015.
-
Detailed analysis of the effects of stencil spatial variations with arbitrary high-order finite-difference Maxwell solver
Authors:
H. Vincenti,
J-L. Vay
Abstract:
Due to discretization effects and truncation to finite domains, many electromagnetic simulations present non-physical modifications of Maxwell's equations in space that may generate spurious signals affecting the overall accuracy of the result. Such modifications for instance occur when Perfectly Matched Layers (PMLs) are used at simulation domain boundaries to simulate open media. Another example…
▽ More
Due to discretization effects and truncation to finite domains, many electromagnetic simulations present non-physical modifications of Maxwell's equations in space that may generate spurious signals affecting the overall accuracy of the result. Such modifications for instance occur when Perfectly Matched Layers (PMLs) are used at simulation domain boundaries to simulate open media. Another example is the use of arbitrary order Maxwell solver with domain decomposition technique that may under some condition involve stencil truncations at subdomain boundaries, resulting in small spurious errors that do eventually build up. In each case, a careful evaluation of the characteristics and magnitude of the errors resulting from these approximations, and their impact at any frequency and angle, requires detailed analytical and numerical studies. To this end, we present a general analytical approach that enables the evaluation of numerical discretization errors of fully three-dimensional arbitrary order finite-difference Maxwell solver, with arbitrary modification of the local stencil in the simulation domain. The analytical model is validated against simulations of domain decomposition technique and PMLs, when these are used with very high-order Maxwell solver, as well as in the infinite order limit of pseudo-spectral solvers. Results confirm that the new analytical approach enables exact predictions in each case. It also confirms that the domain decomposition technique can be used with very high-order Maxwell solver and a reasonably low number of guard cells with negligible effects on the whole accuracy of the simulation.
△ Less
Submitted 20 July, 2015;
originally announced July 2015.