This issue of Laser and Particle Beams includes papers presented at the “Workshop on Simulations of Ultra Intense Laser Beams Interaction with Matter” held in Bordeaux, France, from 14 to 16 April, 2003. This workshop was organized jointly by the Polytechnic University (Madrid, Spain), University of Bordeaux 1 (Talence, France), University of Milano-Bicocca (Milan, Italy), and Commissariat à l'Energie Atomique (CEA, France), sponsored by the European Science Foundation (ESF), University of Bordeaux 1, and CEA.

The transport of hot electrons in solid density plasma created by a high-intensity subpicosecond laser pulse and the resulting heating and ionization of the bulk of plasma are simulated with the electron kinetic code “FPI.” Both the hot and the thermal electrons are treated kinetically. An analysis of the results in terms of a two Maxwellians fit to the numerically obtained distribution function is made.

Nonlocal electron heat transport in plasmas heated by a high-intensity laser field is investigated. We show how the heat transport is strongly modified by the super-Gaussian character of the energy distribution caused by the strong collisional heating. The change in the collisional electron–ion energy exchange due to this modification of the shape of the electron distribution function is also studied.

In the context of the French Laser Mégajoule (LMJ) fusion research program, direct drive is an alternate to indirect drive to reach ignition and thermonuclear burn. We present recent progress in the direct-drive fusion studies for LMJ. Calculations have shown that the LMJ irradiation uniformity is characterized by long wavelength asymmetries compatible with direct drive requirements. Calculations of the irradiation uniformity in the context of indirect drive beam positioning have been done. We show that non-uniformity can be minimized by repointing the beams. Unfortunately, a time analysis shows that this nonuniformity increases strongly in time above levels usually considered inconsistent for direct drive. Finally, a recent baseline target design is presented and consists of a DT ice shell surrounded by a low-density CH foam wicked with cryogenic DT. This design can potentially reach a gain of 90 with a 1-MJ on-target laser driver. Hydrodynamic stability is increased at the ablation front and the laser–target coupling efficiency achieves 85%.

This paper emphasizes the energy dissipation through collective electromagnetic modes (mostly transverse to the incoming beam) of ultraintense relativistic electrons and nonrelativistic protons interacting with a supercompressed core of deuterium + tritium (DT) thermonuclear fuel. This pattern of beam–plasma interaction documents the fast ignition scenario for inertial confinement fusion.

The electronmagnetic Weibel instability is considered analytically in a linear approximation. Relevant growth rates parameters then highlight density ratios between target and particle beams, as well as transverse temperatures. Significant refinements include mode–mode couplings and collisions with target electrons. The former qualify the so-called quasi-linear (weakly turbulent) approach. Usually, it produces significantly lower growth rates than the linear ones. Collisions enhance them slightly for kc /ωp < 1, and dampen them strongly for kc /ωp ≥ 1. Those results simplify rather drastically for the laser-produced and nonrelativistic proton beams. In this case, those growth rates remain always negative through a wide range of beam–target parameters.

For the classic diffusion description of radiative transfer, the specific intensity can be represented by a small angular deviation of the local Planckian equilibrium. In a transparent media, the angular anisotropy becomes strong and one has to solve the general transfer equation. We propose a hierarchy of models that can describe the regime that lies between those two limits. Every member of this family is hyperbolic, flux-limited, and possesses a locally dissipated entropy. This hierarchy also formally recovers the diffusion limit. This study demonstrates that the two-polynomial model is already capable of capturing strong anisotropies.

Experiments of heating of solid targets by fast electrons have been analyzed by means of simulations with a recently developed hybrid code. Electron propagation, refluxing effects, relative importance of self-generated fields, and heating of targets are presented. We found a good agreement between simulations and experiments on the Kα yield.

Collective absorption of intense laser pulses by atomic clusters is studied by PIC simulations. The cluster is modeled in two-dimensional calculations as a cylindrical plasma column with a diameter of D = 6.4 nm and an initial electron density of ne0 = 1023 cm−3. The frequency and intensity dependence of absorption is discussed. It is found that nonresonant absorption by electron emission increases as a power law with the laser intensity. The absorbed energy per electron reaches a maximum of about Wmax = mωp2D2 (ωp: plasma frequency, m: electron mass) at the intensity where ionization saturates.

K-α emission is an intense short-pulse line source well suited for X-ray diagnostic techniques with subpicosecond and micrometer resolution. Numerical simulations are performed here in a search for laser–target interaction regimes where both high efficiency of laser energy transformation to X-ray emission and ultrashort X-ray pulses are achieved. We use the one-dimensional PIC code for the description of the laser interaction with the plasma layer at the target surface. Fast electron transport into the target is treated by our newly developed Monte Carlo code with temporal resolution that is described here in detail. Our simulations reveal extremely short ∼200 fs FWHM bright K-α X-ray pulses emitted from targets heated by 120-fs pulses of a table-top laser. Laser energy conversion efficiency to K-α line emission as high as 6 × 10−5 is noticed. Integration of the emitted energy over the focal spot is carried out to improve the simulation accord with published experimental data. Negligible impact of self-induced electric fields on K-α emission is found for conducting target materials at moderate laser intensities [lsim ]1017 W/cm2.

Hole boring and fast ignition seem to exclude each other: When there is hole boring, no ignition occurs, and vice versa. The laser beam pressure only causes a more or less deep cone-shaped critical surface that leads to better guidance of the beam and to improved laser–plasma coupling. At laser wavelengths of the order of 1 μm, successful fast ignition requires strong anomalous laser beam–pellet coupling.

Effective ion acceleration of picosecond-duration well-collimated bunches in the strong relativistic interaction of a short laser pulse with a thin solid target has been experimentally demonstrated. In this work, with reference to the sharp rear solid–vacuum interface, where ion energization takes place, the one-dimensional Poisson–Boltzmann equation is analytically solved on a finite spatial interval whose extension is determined by requiring electron energy conservation, resulting in the consistent spatial distributions of the hot electrons created by the laser and of the corresponding electrostatic potential. Then, the equation of motions for an ensemble of test ions, initially distributed in a thin layer of the rear target surface, with different initial conditions, is solved and the energy spectrum corresponding to a given initial ion distribution is determined.

Foams are an important component of various inertial confinement fusion target schemes. The propagation of shock waves in foams is also an important issue for many other laser experiments (e.g., laboratory astrophysics experiments). The usual approach is to assume that the foam can be considered as a homogeneous mixture. Taking into account the effects of foam heterogeneity leads to increased shock velocity and reduced compressibility. Evidence of this effect is obtained using two-dimensional adaptive mesh refinement Eulerian numerical simulations. Nonetheless, for the very low density foams filled with DT ice used in recent direct drive target designs, the homogeneous mixture model provides adequate shock timing, density, and pressure profiles.

We consider ion acceleration at the front surface of overdense plasma by a short laser pulse. In this configuration, the laser ponderomotive pressure pushes the background electrons, and a double layer is produced at the boundary of the overdense region. The ions are accelerated by the electrostatic field of the double layer. If the laser intensity is so large that the plasma becomes relativistically transparent, then ion trapping in the running double layer and acceleration to relativistic energies is possible. We study this physics using one-dimensional particle-in-cell simulations. Ion acceleration in one- and two-component plasmas is considered. It is shown that the proton acceleration is more effective when they represent only a small dope to the heavy background ions.

A one-dimensional model for a laser-ablated slab under acceleration g is developed. A characteristic value gc is found to separate two solutions: Lower g values allow sonic or subsonic flow at the critical surface; for higher g the sonic point approaches closer and closer to the slab surface. A significant reduction in the ablation pressure is found in comparison to the g = 0 case. A simple dependence law between the ablation pressure and the slab acceleration, from the initial value g0 to infinity, is identified. Results compared well with fully hydrodynamic computer simulations with Multi2D code. The model has also been found a key step to produce indefinitely steady numerical solutions to study Rayleigh–Taylor instabilities in heat ablation fronts, and validate other theoretical analysis of the problem.

The forthcoming laser installations related to inertial confinement fusion, Laser Mégajoule (LMJ) (France) and National Ignition Facility (NIF) (USA), require multidimensional numerical simulation tools for interpreting current experimental data and to perform predictive modeling for future experiments. Simulations of macroscopic plasma volumes of the order of 1 mm3 and laser exposure times of the order of hundreds of picoseconds are necessary. We present recent developments in the PARAX code towards this goal. The laser field is treated in a standard paraxial approximation in three dimensions. The plasma response is described by single-fluid, two-temperature, fully nonlinear hydrodynamical equations in the plane transverse to the laser propagation axis. The code also accounts for the dominant nonlocal transport terms in spectral form originating from a linearized solution to the Fokker–Planck equation. The simulations of interest are hohlraum plasmas in the case of indirect drive or the plasma corona for direct drive. Recent experimental results on plasma-induced smoothing of RPP laser beams are used to validate the code.