Pseudospectral Maxwell solvers for an accurate modeling of Doppler harmonic generation on plasma mirrors with particle-in-cell codes

Blaclard, G.; Vincenti, H.; Lehe, R.; Vay, J. L.

doi:10.1103/physreve.96.033305

Title: Pseudospectral Maxwell solvers for an accurate modeling of Doppler harmonic generation on plasma mirrors with particle-in-cell codes

Abstract

With the advent of petawatt class lasers, the very large laser intensities attainable on target should enable the production of intense high-order Doppler harmonics from relativistic laser-plasma mirror 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 to hundreds of harmonic orders on a broad range of angles. In particular, we show here that due to the numerical dispersion of waves they induce in vacuum, standard finite difference time domain (FDTD) Maxwell solvers employed in most PIC codes can induce a spurious angular deviation of harmonic beams potentially degrading simulation results. This effect was extensively studied and a simple toy model based on the Snell-Descartes law was developed that allows us to finely predict the angular deviation of harmonics depending on the spatiotemporal 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 spatiotemporal resolution preventing realistic three-dimensional (3D) simulations even on the largest computers available at the time of writing. Finally, we finally show that nondispersive pseudospectral analytical time domain solvers can considerably reduce the spatiotemporal resolution required to mitigatemore »« less

Authors:

Blaclard, G. ^[1]; Vincenti, H. ^[2]; Lehe, R. ^[1]; Vay, J. L. ^[1]

Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States)
Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States); CEA Saclay, Gif-sur-Yvette (France)

Publication Date:: Mon Sep 11 00:00:00 EDT 2017

Research Org.:: Lawrence Berkeley National Laboratory (LBNL), Berkeley, CA (United States)

Sponsoring Org.:: USDOE Office of Science (SC)

OSTI Identifier:: 1582612

Alternate Identifier(s):: OSTI ID: 1389153

Grant/Contract Number:: AC02-05CH11231; AC02-06CH11357

Resource Type:: Accepted Manuscript

Journal Name:: Physical Review. E

Additional Journal Information:: Journal Volume: 96; Journal Issue: 3; Journal ID: ISSN 2470-0045

Publisher:: American Physical Society (APS)

Country of Publication:: United States

Language:: English

Subject:: 70 PLASMA PHYSICS AND FUSION TECHNOLOGY; high intensity laser-plasma interactions; high-harmonic generation; particle-in-cell methods

Citation Formats


                    Blaclard, G., Vincenti, H., Lehe, R., and Vay, J. L. Pseudospectral Maxwell solvers for an accurate modeling of Doppler harmonic generation on plasma mirrors with particle-in-cell codes.  United States: N. p., 2017. 
Web.  doi:10.1103/physreve.96.033305.

Copy to clipboard


                    Blaclard, G., Vincenti, H., Lehe, R., & Vay, J. L. Pseudospectral Maxwell solvers for an accurate modeling of Doppler harmonic generation on plasma mirrors with particle-in-cell codes.  United States.  https://doi.org/10.1103/physreve.96.033305

Copy to clipboard


                    Blaclard, G., Vincenti, H., Lehe, R., and Vay, J. L. Mon .  
"Pseudospectral Maxwell solvers for an accurate modeling of Doppler harmonic generation on plasma mirrors with particle-in-cell codes".  United States.  https://doi.org/10.1103/physreve.96.033305.  https://www.osti.gov/servlets/purl/1582612.

Copy to clipboard


                    
@article{osti_1582612,

  title        = {Pseudospectral Maxwell solvers for an accurate modeling of Doppler harmonic generation on plasma mirrors with particle-in-cell codes},

  author       = {Blaclard, G. and Vincenti, H. and Lehe, R. and Vay, J. L.},

  abstractNote = {With the advent of petawatt class lasers, the very large laser intensities attainable on target should enable the production of intense high-order Doppler harmonics from relativistic laser-plasma mirror 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 to hundreds of harmonic orders on a broad range of angles. In particular, we show here that due to the numerical dispersion of waves they induce in vacuum, standard finite difference time domain (FDTD) Maxwell solvers employed in most PIC codes can induce a spurious angular deviation of harmonic beams potentially degrading simulation results. This effect was extensively studied and a simple toy model based on the Snell-Descartes law was developed that allows us to finely predict the angular deviation of harmonics depending on the spatiotemporal 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 spatiotemporal resolution preventing realistic three-dimensional (3D) simulations even on the largest computers available at the time of writing. Finally, we finally show that nondispersive pseudospectral analytical time domain solvers can considerably reduce the spatiotemporal resolution required to mitigate this spurious deviation and should enable in the near future 3D accurate modeling on supercomputers in a realistic time to solution.},

  doi          = {10.1103/physreve.96.033305},

  journal      = {Physical Review. E},

  number       = 3,

  volume       = 96,

  place        = {United States},

  year         = {Mon Sep 11 00:00:00 EDT 2017},

  month        = {Mon Sep 11 00:00:00 EDT 2017}

}

Copy to clipboard

Journal Article:

Free Publicly Available Full Text

Accepted Manuscript (DOE)

Publisher's Version of Record

https://doi.org/10.1103/physreve.96.033305

Other availability

Search WorldCat to find libraries that may hold this journal

Citation Metrics:

Cited by: 17 works

Citation information provided by
Web of Science

Save / Share:

Export Metadata

Save to My Library

Works referenced in this record:

Attosecond lighthouses from plasma mirrors
journal, November 2012

Wheeler, Jonathan A.; Borot, Antonin; Monchocé, Sylvain
Nature Photonics, Vol. 6, Issue 12
DOI: 10.1038/nphoton.2012.284

A domain decomposition method for pseudo-spectral electromagnetic simulations of plasmas
journal, June 2013

Vay, Jean-Luc; Haber, Irving; Godfrey, Brendan B.
Journal of Computational Physics, Vol. 243
DOI: 10.1016/j.jcp.2013.03.010

Publisher’s Note: Relativistic Generation of Isolated Attosecond Pulses in a $λ^{3}$ Focal Volume [Phys. Rev. Lett. 92 , 063902 (2004)]
journal, February 2004

Naumova, N. M.; Nees, J. A.; Sokolov, I. V.
Physical Review Letters, Vol. 92, Issue 8
DOI: 10.1103/PhysRevLett.92.089901

The PSTD algorithm: A time-domain method requiring only two cells per wavelength
journal, June 1997

Liu, Q. H.
Microwave and Optical Technology Letters, Vol. 15, Issue 3
DOI: 10.1002/(SICI)1098-2760(19970620)15:3%3C158::AID-MOP11%3E3.0.CO;2-3

Attosecond Lighthouses: How To Use Spatiotemporally Coupled Light Fields To Generate Isolated Attosecond Pulses
journal, March 2012

Vincenti, H.; Quéré, F.
Physical Review Letters, Vol. 108, Issue 11
DOI: 10.1103/PhysRevLett.108.113904

Novel methods in the Particle-In-Cell accelerator Code-Framework Warp
journal, January 2012

Vay, J-L; Grote, D. P.; Cohen, R. H.
Computational Science & Discovery, Vol. 5, Issue 1
DOI: 10.1088/1749-4699/5/1/014019

Optics in the relativistic regime
journal, April 2006

Mourou, Gerard A.; Tajima, Toshiki; Bulanov, Sergei V.
Reviews of Modern Physics, Vol. 78, Issue 2
DOI: 10.1103/RevModPhys.78.309

Characteristics of an envelope model for laser–plasma accelerator simulation
journal, January 2011

Cowan, Benjamin M.; Bruhwiler, David L.; Cormier-Michel, Estelle
Journal of Computational Physics, Vol. 230, Issue 1
DOI: 10.1016/j.jcp.2010.09.009

Optical properties of relativistic plasma mirrors
journal, March 2014

Vincenti, H.; Monchocé, S.; Kahaly, S.
Nature Communications, Vol. 5, Issue 1
DOI: 10.1038/ncomms4403

Ultrarelativistic nanoplasmonics as a route towards extreme-intensity attosecond pulses
journal, October 2011

Gonoskov, A. A.; Korzhimanov, A. V.; Kim, A. V.
Physical Review E, Vol. 84, Issue 4
DOI: 10.1103/PhysRevE.84.046403

Compression of amplified chirped optical pulses
journal, December 1985

Strickland, Donna; Mourou, Gerard
Optics Communications, Vol. 56, Issue 3, p. 219-221
DOI: 10.1016/0030-4018(85)90120-8

Theory of high-order harmonic generation in relativistic laser interaction with overdense plasma
journal, October 2006

Baeva, T.; Gordienko, S.; Pukhov, A.
Physical Review E, Vol. 74, Issue 4
DOI: 10.1103/PhysRevE.74.046404

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

Vincenti, H.; Lobet, M.; Lehe, R.
Computer Physics Communications, Vol. 210
DOI: 10.1016/j.cpc.2016.08.023

Route to intense single attosecond pulses
journal, January 2006

Tsakiris, George D.; Eidmann, Klaus; Meyer-ter-Vehn, Jürgen
New Journal of Physics, Vol. 8
DOI: 10.1088/1367-2630/8/1/019

Bright Multi-keV Harmonic Generation from Relativistically Oscillating Plasma Surfaces
journal, August 2007

Dromey, B.; Kar, S.; Bellei, C.
Physical Review Letters, Vol. 99, Issue 8
DOI: 10.1103/PhysRevLett.99.085001

Numerical methods for instability mitigation in the modeling of laser wakefield accelerators in a Lorentz-boosted frame
journal, July 2011

Vay, J. -L.; Geddes, C. G. R.; Cormier-Michel, E.
Journal of Computational Physics, Vol. 230, Issue 15
DOI: 10.1016/j.jcp.2011.04.003

High-order harmonic and attosecond pulse generation on plasma mirrors: basic mechanisms
journal, October 2010

Thaury, C.; Quéré, F.
Journal of Physics B: Atomic, Molecular and Optical Physics, Vol. 43, Issue 21
DOI: 10.1088/0953-4075/43/21/213001

Detailed analysis of the effects of stencil spatial variations with arbitrary high-order finite-difference Maxwell solver
journal, March 2016

Vincenti, H.; Vay, J. -L.
Computer Physics Communications, Vol. 200
DOI: 10.1016/j.cpc.2015.11.009

High-Order Finite Differences and the Pseudospectral Method on Staggered Grids
journal, August 1990

Fornberg, Bengt
SIAM Journal on Numerical Analysis, Vol. 27, Issue 4
DOI: 10.1137/0727052

High harmonic generation in the relativistic limit
journal, July 2006

Dromey, B.; Zepf, M.; Gopal, A.
Nature Physics, Vol. 2, Issue 7
DOI: 10.1038/nphys338

Numerical solution of initial boundary value problems involving maxwell's equations in isotropic media
journal, May 1966

Kane Yee,
IEEE Transactions on Antennas and Propagation, Vol. 14, Issue 3
DOI: 10.1109/TAP.1966.1138693

Plasma mirrors for ultrahigh-intensity optics
journal, April 2007

Thaury, C.; Quéré, F.; Geindre, J. -P.
Nature Physics, Vol. 3, Issue 6
DOI: 10.1038/nphys595

A high-accuracy realization of the Yee algorithm using non-standard finite differences
journal, June 1997

Cole, J. B.
IEEE Transactions on Microwave Theory and Techniques, Vol. 45, Issue 6
DOI: 10.1109/22.588615

High-accuracy Yee algorithm based on nonstandard finite differences: new developments and verifications
journal, September 2002

Cole, J. B.
IEEE Transactions on Antennas and Propagation, Vol. 50, Issue 9
DOI: 10.1109/TAP.2002.801268

Complete characterization of a plasma mirror for the production of high-contrast ultraintense laser pulses
journal, February 2004

Doumy, G.; Quéré, F.; Gobert, O.
Physical Review E, Vol. 69, Issue 2
DOI: 10.1103/PhysRevE.69.026402

Prepulse energy suppression for high-energy ultrashort pulses using self-induced plasma shuttering
journal, January 1991

Kapteyn, Henry C.; Szoke, Abraham; Falcone, Roger W.
Optics Letters, Vol. 16, Issue 7
DOI: 10.1364/OL.16.000490

Short‐pulse laser harmonics from oscillating plasma surfaces driven at relativistic intensity
journal, September 1996

Lichters, R.; Meyer‐ter‐Vehn, J.; Pukhov, A.
Physics of Plasmas, Vol. 3, Issue 9
DOI: 10.1063/1.871619

Bubble regime of wake field acceleration: similarity theory and optimal scalings
journal, January 2006

Pukhov, A.; Gordienko, S.
Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, Vol. 364, Issue 1840
DOI: 10.1098/rsta.2005.1727

Works referencing / citing this record:

Identification of Coupling Mechanisms between Ultraintense Laser Light and Dense Plasmas
journal, March 2019

Chopineau, L.; Leblanc, A.; Blaclard, G.
Physical Review X, Vol. 9, Issue 1
DOI: 10.1103/physrevx.9.011050

Identification of coupling mechanisms between ultraintense laser light and dense plasmas
text, January 2018

Chopineau, L.; Leblanc, A.; Blaclard, G.
arXiv
DOI: 10.48550/arxiv.1809.03903

Probing Strong-Field QED with Doppler-Boosted Petawatt-Class Lasers
journal, September 2021

Fedeli, L.; Sainte-Marie, A.; Zaim, N.
Physical Review Letters, Vol. 127, Issue 11
DOI: 10.1103/physrevlett.127.114801

A generalized massively parallel ultra-high order FFT-based Maxwell solver
text, January 2018

Kallala, Haithem; Vay, Jean-Luc; Vincenti, Henri
arXiv
DOI: 10.48550/arxiv.1812.07357

Similar Records in DOE PAGES and OSTI.GOV collections:

Pseudospectral Maxwell solvers for an accurate modeling of Doppler harmonic generation on plasma mirrors with particle-in-cell codes

Journal Article Blaclard, G. ; Vincenti, H. ; Lehe, R. ; ... - Physical Review. E

With the advent of petawatt class lasers, the very large laser intensities attainable on target should enable the production of intense high-order Doppler harmonics from relativistic laser-plasma mirror 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 to hundreds of harmonic orders on a broad range of angles. In particular, we show here that due to the numerical dispersion of waves they induce in vacuum, standard finite difference time domain (FDTD) Maxwell solvers employed in most PIC codes can induce a spurious angular deviation of harmonicmore »« less
https://doi.org/10.1103/PhysRevE.96.033305
Controlling the numerical Cerenkov instability in PIC simulations using a customized finite difference Maxwell solver and a local FFT based current correction

Journal Article Li, Fei ; Yu, Peicheng ; Xu, Xinlu ; ... - Computer Physics Communications

In this study we present a customized finite-difference-time-domain (FDTD) Maxwell solver for the particle-in-cell (PIC) algorithm. The solver is customized to effectively eliminate the numerical Cerenkov instability (NCI) which arises when a plasma (neutral or non-neutral) relativistically drifts on a grid when using the PIC algorithm. We control the EM dispersion curve in the direction of the plasma drift of a FDTD Maxwell solver by using a customized higher order finite difference operator for the spatial derivative along the direction of the drift (1ˆ direction). We show that this eliminates the main NCI modes with moderate |k₁|, while keeps additionalmore »« less
Cited by 24
https://doi.org/10.1016/j.cpc.2017.01.001

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

Journal Article Kallala, Haithem ; Vay, Jean-Luc ; Vincenti, Henri - Computer Physics Communications

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 boundary. 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 regionsmore »« less
Cited by 8
https://doi.org/10.1016/j.cpc.2019.07.009

Full Text Available
Detailed analysis of the effects of stencil spatial variations with arbitrary high-order finite-difference Maxwell solver

Journal Article Vincenti, H. ; Vay, J. -L. - Computer Physics Communications

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 themore »« less
Cited by 29
https://doi.org/10.1016/j.cpc.2015.11.009

Full Text Available
Quasi-Helmholtz decomposition, Gauss' laws and charge conservation for finite element particle-in-cell

Journal Article O'Connor, Scott ; Crawford, Zane D. ; Ramachandran, O. H. ; ... - Computer Physics Communications

Development of particle-in-cell (PIC) methods using finite element based methods (FEMs) have been a topic of renewed interest; this has largely been driven by (a) the ability of finite element methods to better model geometry, (b) better understanding of function spaces that are necessary to represent all Maxwell quantities, and (c) more recently, the fundamental rubrics that should be obeyed in space and time so as to satisfy Gauss' laws and the equation of continuity. In that vein, methods have been developed recently that satisfy these equations and are agnostic to time stepping methods. While this development is indeed amore »« less
https://doi.org/10.1016/j.cpc.2022.108345

Full Text Available

Similar Records

Title: Pseudospectral Maxwell solvers for an accurate modeling of Doppler harmonic generation on plasma mirrors with particle-in-cell codes

Abstract

Citation Formats

Attosecond lighthouses from plasma mirrors journal, November 2012

A domain decomposition method for pseudo-spectral electromagnetic simulations of plasmas journal, June 2013

Publisher’s Note: Relativistic Generation of Isolated Attosecond Pulses in a λ 3 Focal Volume [Phys. Rev. Lett. 92 , 063902 (2004)] journal, February 2004

The PSTD algorithm: A time-domain method requiring only two cells per wavelength journal, June 1997

Attosecond Lighthouses: How To Use Spatiotemporally Coupled Light Fields To Generate Isolated Attosecond Pulses journal, March 2012

Novel methods in the Particle-In-Cell accelerator Code-Framework Warp journal, January 2012

Optics in the relativistic regime journal, April 2006

Characteristics of an envelope model for laser–plasma accelerator simulation journal, January 2011

Optical properties of relativistic plasma mirrors journal, March 2014

Ultrarelativistic nanoplasmonics as a route towards extreme-intensity attosecond pulses journal, October 2011

Compression of amplified chirped optical pulses journal, December 1985

Theory of high-order harmonic generation in relativistic laser interaction with overdense plasma journal, October 2006

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

Route to intense single attosecond pulses journal, January 2006

Bright Multi-keV Harmonic Generation from Relativistically Oscillating Plasma Surfaces journal, August 2007

Numerical methods for instability mitigation in the modeling of laser wakefield accelerators in a Lorentz-boosted frame journal, July 2011

High-order harmonic and attosecond pulse generation on plasma mirrors: basic mechanisms journal, October 2010

Detailed analysis of the effects of stencil spatial variations with arbitrary high-order finite-difference Maxwell solver journal, March 2016

High-Order Finite Differences and the Pseudospectral Method on Staggered Grids journal, August 1990

High harmonic generation in the relativistic limit journal, July 2006

Numerical solution of initial boundary value problems involving maxwell's equations in isotropic media journal, May 1966

Plasma mirrors for ultrahigh-intensity optics journal, April 2007

A high-accuracy realization of the Yee algorithm using non-standard finite differences journal, June 1997

High-accuracy Yee algorithm based on nonstandard finite differences: new developments and verifications journal, September 2002

Complete characterization of a plasma mirror for the production of high-contrast ultraintense laser pulses journal, February 2004

Prepulse energy suppression for high-energy ultrashort pulses using self-induced plasma shuttering journal, January 1991

Short‐pulse laser harmonics from oscillating plasma surfaces driven at relativistic intensity journal, September 1996

Bubble regime of wake field acceleration: similarity theory and optimal scalings journal, January 2006

Identification of Coupling Mechanisms between Ultraintense Laser Light and Dense Plasmas journal, March 2019

Identification of coupling mechanisms between ultraintense laser light and dense plasmas text, January 2018

Probing Strong-Field QED with Doppler-Boosted Petawatt-Class Lasers journal, September 2021

A generalized massively parallel ultra-high order FFT-based Maxwell solver text, January 2018

Attosecond lighthouses from plasma mirrors
journal, November 2012

A domain decomposition method for pseudo-spectral electromagnetic simulations of plasmas
journal, June 2013

Publisher’s Note: Relativistic Generation of Isolated Attosecond Pulses in a $λ^{3}$ Focal Volume [Phys. Rev. Lett. 92 , 063902 (2004)]
journal, February 2004

The PSTD algorithm: A time-domain method requiring only two cells per wavelength
journal, June 1997

Attosecond Lighthouses: How To Use Spatiotemporally Coupled Light Fields To Generate Isolated Attosecond Pulses
journal, March 2012

Novel methods in the Particle-In-Cell accelerator Code-Framework Warp
journal, January 2012

Optics in the relativistic regime
journal, April 2006

Characteristics of an envelope model for laser–plasma accelerator simulation
journal, January 2011

Optical properties of relativistic plasma mirrors
journal, March 2014

Ultrarelativistic nanoplasmonics as a route towards extreme-intensity attosecond pulses
journal, October 2011

Compression of amplified chirped optical pulses
journal, December 1985

Theory of high-order harmonic generation in relativistic laser interaction with overdense plasma
journal, October 2006

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

Route to intense single attosecond pulses
journal, January 2006

Bright Multi-keV Harmonic Generation from Relativistically Oscillating Plasma Surfaces
journal, August 2007

Numerical methods for instability mitigation in the modeling of laser wakefield accelerators in a Lorentz-boosted frame
journal, July 2011

High-order harmonic and attosecond pulse generation on plasma mirrors: basic mechanisms
journal, October 2010

Detailed analysis of the effects of stencil spatial variations with arbitrary high-order finite-difference Maxwell solver
journal, March 2016

High-Order Finite Differences and the Pseudospectral Method on Staggered Grids
journal, August 1990

High harmonic generation in the relativistic limit
journal, July 2006

Numerical solution of initial boundary value problems involving maxwell's equations in isotropic media
journal, May 1966

Plasma mirrors for ultrahigh-intensity optics
journal, April 2007

A high-accuracy realization of the Yee algorithm using non-standard finite differences
journal, June 1997

High-accuracy Yee algorithm based on nonstandard finite differences: new developments and verifications
journal, September 2002

Complete characterization of a plasma mirror for the production of high-contrast ultraintense laser pulses
journal, February 2004

Prepulse energy suppression for high-energy ultrashort pulses using self-induced plasma shuttering
journal, January 1991

Short‐pulse laser harmonics from oscillating plasma surfaces driven at relativistic intensity
journal, September 1996

Bubble regime of wake field acceleration: similarity theory and optimal scalings
journal, January 2006

Identification of Coupling Mechanisms between Ultraintense Laser Light and Dense Plasmas
journal, March 2019

Identification of coupling mechanisms between ultraintense laser light and dense plasmas
text, January 2018

Probing Strong-Field QED with Doppler-Boosted Petawatt-Class Lasers
journal, September 2021

A generalized massively parallel ultra-high order FFT-based Maxwell solver
text, January 2018