2020/10/31 by Elena Panova, Valentin Volokitin, Evgeny Efimenko +9 · 6 citations
Physics and Astronomy · #Computation #Computational electromagnetics #Electromagnetic field #Field (mathematics) #High resolution #Laser-Matter Interactions and Applications #Laser-Plasma Interactions and Diagnostics #Node (physics) #Orbital Angular Momentum in Optics #Python (programming language) #Stability (learning theory) #physics.comp-ph
paper · pdf · doi:10.3390/app11030956
published in Applied Sciences 11(3), 956 (Multidisciplinary Digital Publishing Institute) · 27 pages, 8 figures
openalex created_date 2020/10/08 · openalex publication_date 2021/01/21 · arxiv created 2021/01/22 · arxiv updated 2021/01/25 · openalex updated_date 2026/08/05
When a pulsed, few-cycle electromagnetic wave is focused by optics with f-number smaller than two, the frequency components it contains are focused to different regions of space, building up a complex electromagnetic field structure. Accurate numerical computation of this structure is essential for many applications such as the analysis, diagnostics, and control of high-intensity laser-matter interactions. However, straightforward use of finite-difference methods can impose unacceptably high demands on computational resources, owing to the necessity of resolving far-field and near-field zones at sufficiently high resolution to overcome numerical dispersion effects. Here, we present a procedure for fast computation of tight focusing by mapping a spherically curved far-field region to periodic space, where the field can be advanced by a dispersion-free spectral solver. In many cases of interest, the mapping reduces both run time and memory requirements by a factor of order 10, making it possible to carry out simulations on a desktop machine or a single node of a supercomputer. We provide an open-source C++ implementation with Python bindings and demonstrate its use for a desktop machine, where the routine provides the opportunity to use the resolution sufficient for handling the pulses with spectra spanning over several octaves. The described approach can facilitate the stability analysis of theoretical proposals, the studies based on statistical inferences, as well as the overall development and analysis of experiments with tightly-focused short laser pulses.