Particle-in-cell (PIC) methods have a long history in the study of laser-plasma interactions. Early electromagnetic codes used the Yee staggered grid for field variables combined with a leapfrog EM-field update and the Boris algorithm for particle pushing. The general properties of such schemes are well documented. Modern PIC codes tend to add to these high-order shape functions for particles, Poisson preserving field updates, collisions, ionisation, a hybrid scheme for solid density and high-field QED effects. In addition to these physics packages, the increase in computing power now allows simulations with real mass ratios, full 3D dynamics and multi-speckle interaction. This paper presents a review of the core algorithms used in current laser-plasma specific PIC codes. Also reported are estimates of self-heating rates, convergence of collisional routines and test of ionisation models which are not readily available elsewhere. Having reviewed the status of PIC algorithms we present a summary of recent applications of such codes in laser-plasma physics, concentrating on SRS, short-pulse laser-solid interactions, fast-electron transport, and QED effects.

/pdf/contemporary-particle-in-cell-approach-to-laser-plasma-1hb69lir4e.pdf

Contemporary particle-in-cell approach to laser-plasma modelling

Magnetic reconnection is important to the dynamics of many astrophysical and fusion plasmas but our understanding of it is incomplete. Petaflop-scale simulations of the evolution of turbulent magnetic reconnection in a three-dimensional plasma indicate that it proceeds in a way that is dramatically different from classical theory.

/pdf/role-of-electron-physics-in-the-development-of-turbulent-4dvxxy5olc.pdf

Role of electron physics in the development of turbulent magnetic reconnection in collisionless plasmas

The direct-drive, laser-based approach to inertial confinement fusion (ICF) is reviewed from its inception following the demonstration of the first laser to its implementation on the present generation of high-power lasers. The review focuses on the evolution of scientific understanding gained from target-physics experiments in many areas, identifying problems that were demonstrated and the solutions implemented. The review starts with the basic understanding of laser–plasma interactions that was obtained before the declassification of laser-induced compression in the early 1970s and continues with the compression experiments using infrared lasers in the late 1970s that produced thermonuclear neutrons. The problem of suprathermal electrons and the target preheat that they caused, associated with the infrared laser wavelength, led to lasers being built after 1980 to operate at shorter wavelengths, especially 0.35 μm—the third harmonic of the Nd:glass laser—and 0.248 μm (the KrF gas laser). The main physics areas relevant to direct drive are reviewed. The primary absorption mechanism at short wavelengths is classical inverse bremsstrahlung. Nonuniformities imprinted on the target by laser irradiation have been addressed by the development of a number of beam-smoothing techniques and imprint-mitigation strategies. The effects of hydrodynamic instabilities are mitigated by a combination of imprint reduction and target designs that minimize the instability growth rates. Several coronal plasma physics processes are reviewed. The two-plasmon–decay instability, stimulated Brillouin scattering (together with cross-beam energy transfer), and (possibly) stimulated Raman scattering are identified as potential concerns, placing constraints on the laser intensities used in target designs, while other processes (self-focusing and filamentation, the parametric decay instability, and magnetic fields), once considered important, are now of lesser concern for mainline direct-drive target concepts. Filamentation is largely suppressed by beam smoothing. Thermal transport modeling, important to the interpretation of experiments and to target design, has been found to be nonlocal in nature. Advances in shock timing and equation-of-state measurements relevant to direct-drive ICF are reported. Room-temperature implosions have provided an increased understanding of the importance of stability and uniformity. The evolution of cryogenic implosion capabilities, leading to an extensive series carried out on the 60-beam OMEGA laser [Boehly et al., Opt. Commun. 133, 495 (1997)], is reviewed together with major advances in cryogenic target formation. A polar-drive concept has been developed that will enable direct-drive–ignition experiments to be performed on the National Ignition Facility [Haynam et al., Appl. Opt. 46(16), 3276 (2007)]. The advantages offered by the alternative approaches of fast ignition and shock ignition and the issues associated with these concepts are described. The lessons learned from target-physics and implosion experiments are taken into account in ignition and high-gain target designs for laser wavelengths of 1/3 μm and 1/4 μm. Substantial advances in direct-drive inertial fusion reactor concepts are reviewed. Overall, the progress in scientific understanding over the past five decades has been enormous, to the point that inertial fusion energy using direct drive shows significant promise as a future environmentally attractive energy source.

Direct-drive inertial confinement fusion: A review

The plasma membrane of eukaryotic cells contains several types of lipids displaying high biochemical variability in both their apolar moiety (e.g. the acyl chain of glycerolipids) and their polar head (e.g. the sugar structure of glycosphingolipids). Among these lipids, cholesterol is unique because its biochemical variability is almost exclusively restricted to the oxidation of its polar -OH group. Although generally considered the most rigid membrane lipid, cholesterol can adopt a broad range of conformations due to the flexibility of its isooctyl chain linked to the polycyclic sterane backbone. Moreover, cholesterol is an asymmetric molecule displaying a planar face and a rough  face. Overall, these structural features open up a number of possible interactions between cholesterol and membrane lipids and proteins, consistent with the prominent regulatory functions that this unique lipid exerts on membrane components. The aim of this review is to describe how cholesterol interacts with membrane lipids and proteins at the molecular/atomic scale, with special emphasis on transmembrane domains of proteins containing either the consensus cholesterol-binding motifs CRAC and CARC or a tilted peptide. Despite their broad structural diversity, all these domains bind cholesterol through common molecular mechanisms, leading to the identification of a subset of amino acid residues that are overrepresented in both linear and three-dimensional membrane cholesterol-binding sites.

/pdf/how-cholesterol-interacts-with-membrane-proteins-an-4fbhlpho06.pdf

How cholesterol interacts with membrane proteins: an exploration of cholesterol-binding sites including CRAC, CARC, and tilted domains

An unsolved problem in plasma turbulence is how energy is dissipated at small scales. Particle collisions are too infrequent in hot plasmas to provide the necessary dissipation. Simulations either treat the fluid scales and impose an ad hoc form of dissipation (e.g., resistivity) or consider dissipation arising from resonant damping of small amplitude disturbances where damping rates are found to be comparable to that predicted from linear theory. Here, we report kinetic simulations that span the macroscopic fluid scales down to the motion of electrons. We find that turbulent cascade leads to generation of coherent structures in the form of current sheets that steepen to electron scales, triggering strong localized heating of the plasma. The dominant heating mechanism is due to parallel electric fields associated with the current sheets, leading to anisotropic electron and ion distributions which can be measured with NASA's upcoming Magnetospheric Multiscale mission. The motion of coherent structures also generates waves that are emitted into the ambient plasma in form of highly oblique compressional and shear Alfven modes. In 3D, modes propagating at other angles can also be generated. This indicates that intermittent plasma turbulence will in general consist of both coherent structures and waves. However, the current sheet heating is found to be locally several orders of magnitude more efficient than wave damping and is sufficient to explain the observed heating rates in the solar wind.

Coherent structures, intermittent turbulence, and dissipation in high-temperature plasmas

The algorithms, implementation details, and applications of VPIC, a state-of-the-art first principles 3D electromagnetic relativistic kinetic particle-in-cell code, are discussed. Unlike most codes, VPIC is designed to minimize data motion, as, due to physical limitations (including the speed of light!), moving data between and even within modern microprocessors is more time consuming than performing computations. As a result, VPIC has achieved unprecedented levels of performance. For example, VPIC can perform ∼0.17 billion cold particles pushed and charge conserving accumulated per second per processor on IBM’s Cell microprocessor—equivalent to sustaining Los Alamos’s planned Roadrunner supercomputer at ∼0.56 petaflop (quadrillion floating point operations per second). VPIC has enabled previously intractable simulations in numerous areas of plasma physics, including magnetic reconnection and laser plasma interactions; next generation supercomputers like Roadrunner will enable further advances.

Ultrahigh performance three-dimensional electromagnetic relativistic kinetic plasma simulationa)

VPIC [1], a first-principles 3d electromagnetic charge-conserving relativistic kinetic particle-in-cell code, was recently adapted to run on Los Alamos's Roadrunner [2], the first supercomputer to break a petaflop (1015 floating point operations per second) in the TOP500 supercomputer performance rankings. [3] We summarize VPIC's modeling capabilities, VPIC's optimization techniques and Roadrunner's computational characteristics. We then discuss three applications enabled by VPIC's unprecedented performance on Roadrunner: modeling laser plasma interaction in upcoming inertial confinement fusion experiments at the National Ignition Facility, modeling short-pulse laser GeV ion acceleration and modeling reconnection in space and laboratory plasmas.

https://iopscience.iop.org/article/10.1088/1742-6596/180/1/012055/pdf

Advances in petascale kinetic plasma simulation with VPIC and Roadrunner

We demonstrate the outstanding performance and scalability of the VPIC kinetic plasma modeling code on the heterogeneous IBM Roadrunner supercomputer at Los Alamos National Laboratory. VPIC is a three-dimensional, relativistic, electromagnetic, particle-in-cell (PIC) code that self-consistently evolves a kinetic plasma. VPIC simulations of laser plasma interaction were conducted at unprecedented fidelity and scale---up to 1.0 x 1012 particles on as many as 136 x 106 voxels---to model accurately the particle trapping physics occurring within a laser-driven hohlraum in an inertial confinement fusion experiment. During a parameter study of laser reflectivity as a function of laser intensity under experimentally realizable hohlraum conditions [1], we measured sustained performance exceeding 0.374 Pflop/s (s.p.) with the inner loop itself achieving 0.488 Pflop/s (s.p.). Given the increasing importance of data motion limitations, it is notable that this was measured in a PIC calculation---a technique that typically requires more data motion per computation than other techniques (such as dense matrix calculations, molecular dynamics N-body calculations and Monte-Carlo calculations) often used to demonstrate supercomputer performance. This capability opens up the exciting possibility of using VPIC to model, from first-principles, an issue critical to the success of the multi-billion dollar DOE/NNSA National Ignition Facility.

0.374 Pflop/s trillion-particle kinetic modeling of laser plasma interaction on Roadrunner

A suite of three-dimensional (3D) VPIC [K. J. Bowers et al., Phys. Plasmas 15, 055703 (2008)] particle-in-cell simulations of backward stimulated Raman scattering (SRS) in inertial confinement fusion hohlraum plasma has been performed on the heterogeneous multicore supercomputer, Roadrunner, presently the world’s most powerful supercomputer. These calculations reveal the complex nonlinear behavior of SRS and point to a new era of “at scale” 3D modeling of SRS in solitary and multiple laser speckles. The physics governing nonlinear saturation of SRS in a laser speckle in 3D is consistent with that of prior two-dimensional (2D) studies [L. Yin et al., Phys. Rev. Lett. 99, 265004 (2007)], but with important differences arising from enhanced diffraction and side loss in 3D compared with 2D. In addition to wave front bowing of electron plasma waves (EPWs) due to trapped electron nonlinear frequency shift and amplitude-dependent damping, we find for the first time that EPW self-focusing, which evolved from trap...

B. Bergen

Papers

Role of electron physics in the development of turbulent magnetic reconnection in collisionless plasmas

Ultrahigh performance three-dimensional electromagnetic relativistic kinetic plasma simulationa)

Advances in petascale kinetic plasma simulation with VPIC and Roadrunner

0.374 Pflop/s trillion-particle kinetic modeling of laser plasma interaction on Roadrunner

Onset and saturation of backward stimulated Raman scattering of laser in trapping regime in three spatial dimensions