Supercomputer simulation opens prospects for obtaining ultra-dense electron-positron plasmas
Physicists from the Institute of Applied Physics of the Russian Academy of Sciences, researchers from Chalmers University of Technology and computer scientists from Lobachevsky University have developed a new software tool called PICADOR for numerical modeling of laser plasmas on modern supercomputers.
The work on the PICADOR software system started in 2010. PICADOR is a parallel implementation of the particle-in-cell method that has been optimized for modern heterogeneous cluster systems. The project combined the competencies and efforts of experts from many fields, thus becoming the basis for the well-thought-out optimization and development of new computing approaches that take into account various physical processes. Eventually, this opened the way for a breakthrough in modeling capabilities in a number of research projects. The system's functional capabilities and performance make it possible to perform numerical simulations in a range of problems at the forefront of modern laser plasma physics.
In their article published in Scientific Reports, Nizhny Novgorod scientists formulated the conditions under which the avalanche-like generation of electrons and positrons in the focus of a high-power laser pulse yields an electron-positron plasma of record density. The study will make it possible to understand processes occurring in astrophysical objects and to study elementary particle production processes.
A well-known fact in quantum physics is the possibility of transforming certain particles into other particles. In particular, in a sufficiently strong electric or magnetic field, a gamma photon can decay into two particles, an electron and a positron. Until now, this effect was observed in laboratory conditions mainly when gamma radiation was transmitted through crystals in which sufficiently strong fields exist near atomic nuclei. Scientists seek a new tool for studying this phenomenon: lasers capable of generating short pulses with a power of more than 10 petawatts. This level of power is achieved by extreme focusing of radiation. For example, scientists suggest using a laser field configuration referred to as dipole focusing. In this case, the focus point is irradiated from all sides. It has been shown theoretically that electron-positron avalanches can be observed at the focus of such a laser facility. Particles created by the decay of a gamma photon will be accelerated by a laser field and will emit gamma photons, which in turn will give rise to new electrons and positrons. As a result, the number of particles in a short time should grow immensely giving rise to a superdense electron-positron plasma.
However, there are some limitations on the density of the plasma that can be obtained in this way. At some point, the laser radiation will not be able to penetrate the plasma that has become too dense, and the avalanche will subside. According to existing estimates, particle concentration in the laser focus will be just over 1024 particles per cubic centimeter. For comparison, approximately the same electron concentration is found in heavy metals, for example, in platinum or gold.
In their new paper, a team of authors headed by Professor A.M. Sergeev, Academician of the Russian Academy of Sciences, showed that under certain conditions, this number can be an order of magnitude higher.
Large-scale numerical simulation of the electron-positron avalanche development in a tightly focused laser field demonstrates a fundamentally new object of investigation, the quasistationary states of a dense electron-positron plasma. These states have a very interesting and unexpected structure. While the laser field in the form of a dipole wave has an axial symmetry, the distribution of electron-positron plasma resulting from the development of the current instability degenerates into two thin layers oriented at a random angle. The thickness of the layers and particle concentration in these layers is apparently limited only by the randomness of the radiation process, which leads to extreme plasma density values. With a total number of particles of the order of 1011, the density exceeds the value of 1026 particles per cubic centimeter, and in our case it was only limited by the resolution of numerical simulation.