Models and experiments were developed to study femtosecond laser ablation of silicon using 800 nm, 40 fs pulses with fluences ranging from 0.5 to 35 J/cm2. At low fluences, ablation was found to occur due to bubble formation and splashing within the melt layer. At higher fluences, it was found that the ablation depth exceeded the melt layer depth due to shockwave ablation. The variation in ion flux and ion velocity was also studied both experimentally and theoretically. It was found that the variation in ion flux is mainly dependent on the variation in the average charge state, with only a small variation in the total number of ions above $\sim \!\!1.5\; \,{\rm J/c}{{\rm m}^2}$. Comparisons between the theoretical and experimental ion flux showed that higher charge state ions received greater portion of the laser energy compared with lower charge state ions.

We analyze the transverse stability for a configuration of multiple Gaussian bunches subject to the self-generated plasma wakefield. Through a semi-analytical approach we first study the equilibrium configuration for the modulated beam and then we investigate the evolution of the equilibrium configuration due to the emittance-driven expansion of the beam front that results in a rigid backward shift. The rear-directed shift brings the modulated beam out of the equilibrium, with the possibility for some of the bunch particles to be lost with a consequent deterioration of the driver. We look therefore for the proper position of the single bunches that maximize the stability without severely affecting the accelerating field behind the driver. We then compare the results with three-dimensional particle in cell simulations.

We use an Eulerian Vlasov code, which solves the one-dimensional relativistic Vlasov–Maxwell equations for both electrons and ions, to follow in details the evolution of the distribution functions and the mechanism of the formation and evolution of double layers during ion acceleration driven by a high-intensity circularly polarized short laser pulse (12 ω−1 where ω is the laser angular frequency) normally incident on a thin dense foil. We compare three cases with a high-density deuterium plasma target of total thickness 1.767 cω−1 and constant n/n cr = 100, where n cr is the critical density, and where the laser intensity is varied from a situation where the target is opaque to the laser pulse (normalized vector potential or quiver momentum a 0 = 80), to a situation where, above a critical laser intensity, a very small fraction of the laser pulse is transmitted through the target (a 0 = 90), and finally to a situation where a more important fraction is transmitted through the target (a 0 = 100). The dynamics of ion and electron acceleration are quite different in the three cases, and are followed in detail by the Eulerian Vlasov code, which allows an accurate representation of the distribution function. In the intermediate case, the Vlasov code has revealed a remarkably well-developed spiral structure in the phase space of the electron distribution function, which is associated with large sawtooth modulations in the electron density profiles.

Development of a detonation wave due to α heating following short pulse laser irradiation in pre-compressed deuterium–tritium (DT) plasma is considered. The laser parameters required for development of a detonation wave are calculated. We find that a laser irradiance and energy of IL = 1.75 × 1023 W/cm2 and 12.8 kJ accordingly during 1.0 ps in a pre-compressed target at 900 g/cm3 creates an α heating fusion detonation wave. In this case, the nuclear fusion ignition conditions for the pre-compressed DT plasma are achieved along the detonation wave orbit.

Terahertz (THz) radiation generation by nonlinear mixing of two laser beams, obliquely incident on an overdense plasma is investigated. In an overdense plasma, the laser beams penetrate to only thin layer of a plasma surface and reflected. At this thin layer, the laser beams exert a ponderomotive force on the electrons of plasma and impart them oscillatory velocity at the different frequency of lasers. THz waves appear in the reflected component from the plasma surface. The amplitude of THz waves can be augmented by applying the magnetic field perpendicular to the direction of propagation of lasers. It is found that the field strength of the emitted THz radiations is sensitive to the angle of incident of the laser beams, beat frequency, and magnetic field strength. In this scheme, the magnetic field strength plays an important role for strong THz wave generation.

Usually, the high-power microwave (HPM) devices suffer from impedance collapse and cathode material degradation or even failure. When the intense electron beam bombards the anode (or named as collector in HPM device), an anode plasma could appear under certain conditions. In this case, the impedance collapse is caused by the expansions of the cathode and anode plasmas and diode current overshot caused by the bipolar flow. In this paper, characterization of a short-pulse high-power diode operated with anode effects with a dielectric fiber (velvet) cathode is discussed. The bipolar flow (or anode plasma) is indeed evident at beam power densities ~11 MW/cm2 and the pulse durations of ~50 ns. The analysis results of the deposit dose and thermal regime of the anode show that the electron stimulated desorption played an important role in the generation of anode plasma in this case. With the effect of anode plasma, the appearance of local cathode plasma flares (or nonuniform electron emission) is particularly detrimental for the diode closure. Micro-structure and elemental surface compositions of cathode are changed by the anode splashing, which is very harmful to the performance of cathode.

The paper presents the results of calculation and optimization of a structure of a radial insulating magnetic field in an acceleration gap of a high-power ion diode. A diode configuration with an induction plasma source and an anode configuration with azimuthally symmetrical slots and a pair of cathode coils of a magnetic diode system have been studied. When the size of the slots is ≤5 mm and codirectional magnetic fields of a diode and shock induction coil, the perturbation of the B-field does not exceed ≤20% and is located in the region near the anode. In this condition, the topology of the magnetic field В = f(1/r) is maintained in the acceleration gap. It was shown that the required radial distribution of the B-field can be optimized by varying the anode profile in the region opposite to two cathode coils of the diode magnetic system.

In inertial fusion, one of scientific issues is to reduce an implosion non-uniformity of a spherical fuel target. The implosion non-uniformity is caused by several factors, including the driver beam illumination non-uniformity, the Rayleigh–Taylor instability (RTI) growth, etc. In this paper, we propose a new control method to reduce the implosion non-uniformity; the oscillating implosion acceleration δg(t) is created by pulsating and dephasing heavy-ion beams (HIBs) in heavy-ion inertial fusion (HIF). The δg(t) would reduce the RTI growth effectively. The original concept of the non-uniformity control in inertial fusion was proposed in [Laser Part. Beams (1993) 11, 757–768]. In this paper, it was found that the pulsating and dephasing HIBs illumination provide successfully the controlled δg(t) and that δg(t) induced by the pulsating HIBs reduces well the implosion non-uniformity. Consequently the pulsating HIBs improve a pellet gain remarkably in HIF.

In this paper, the robustness of the dynamic instability mitigation mechanism is first examined, and then the instability mitigation phenomenon is demonstrated in a deuterium–tritium (DT) fuel target implosion by wobbling heavy-ion beams (HIBs). The results presented here show that the mechanism of the dynamic instability mitigation is rather robust against changes in the phase, the amplitude and the wavelength of the wobbling perturbation applied. In general instability would emerge from the perturbation of the physical quantity. Normally the perturbation phase is unknown, so that the instability growth rate is discussed. However, if the perturbation phase is known, the instability growth can be controlled by a superposition of perturbations imposed actively: if the perturbation is induced by, for example, a driving beam axis oscillation or wobbling, the perturbation phase could be controlled and the instability growth is mitigated by the superposition of the growing perturbations. In this paper, we realize the superposition of the perturbation by the wobbling HIBs’ illumination onto a DT fuel target in heavy-ion inertial fusion (HIF). Our numerical fluid implosion simulations present that the implosion non-uniformity is mitigated successfully by the wobbling HIBs illumination in HIF.

This paper aims at a test of a plasma focus (PF) device as a metallic plasma jet generator. The experiment was carried out at the DPF-1000U device in which the inner electrode face was conically shaped. Deuterium (D2) with the initial pressure of 0.9 Torr was used as a filling gas. In the experiment, a metallic plasma was ensured due to erosion of the inner electrode during a PF discharge. To create the metallic plasma jet, the eroded copper (Cu) plasma, swept by the deuterium plasma sheath, was accelerated axially and compressed to very small radius (about 1–2 mm). The Cu plasma jet achieved a velocity of 3 × 107 cm/s. To study processes of the plasma jet creation and propagation a 16-frame laser interferometer and a four-frame X-ray pinhole camera were used. Recorded images prove a successful adaptation of the PF device to the metallic plasma jet generator.

The emission of X-rays from solid tin targets irradiated by low-energy (few mJ) femtosecond laser pulses propagated through air plasma sparks is investigated. The aim is that to better understand the X-ray emission mechanism and to show the possibility to produce proper radiation for spectroscopic and imaging applications with a table-top laser system. The utilization of a controlled ultrashort prepulse is found necessary to optimize the conversion efficiency of laser energy into Lα radiation. The optimum contrast between the main pulse and the controlled prepulse is found about 102. A correlation between the laser contrast value and the laser near-infrared spectra at the exit of the plasma spark is observed.

We developed a comprehensive model for simulating laser/target interaction in the presence of external axial magnetic fields. The model was integrated into the framework of the HEIGHTS-LPP computer simulation package and benchmarked with recent experimental results. The package was then used to study the angular distribution of extreme ultraviolet (EUV) photon output in plasmas produced in tin planar targets by a Nd:YAG laser. A moderate (0.5 T) permanent magnetic field does not affect EUV source evolution and output. More effective control of plasma plume expansion should be based on magnetic field gradients, that is, on the temporary varying magnetic fields as a magnetic pinch. Analysis of angular EUV output showed strong anisotropy of photon emissions. We found that the correct monitoring angle (i.e., at which the measured EUV flux corresponds to the averaged value after the correctly integrated angular distribution) does not depend on laser irradiance in the studied range and is equivalent to ~60°. We recommend arranging the EUV sensors accordingly in experiments with planar tin targets.

As the energy spread of intense pulsed electron beams (IPEB) strongly influences the irradiation effects, it has been of great importance to characterize the IPEB energy spectrum. With the combination of Child–Langmuir law and Monte Carlo simulation, the IPEB energy spectrum has been obtained in this work by transformation from the accelerating voltage applied to the diode. To verify the accuracy of this simple algorithm, a magnetic spectrometer with an imaging plate was designed to test the IPEB energy spectrum. The measurement was completed with IPEB generated by explosive emission electron diode, the pulse duration, maximum electron energy, total beam current being 80 ns, 450 keV, and 1 kA, respectively. The results verified the reliability of the above analysis method for energy spectrum, which can avoid intercepting the beam, and at the same time significantly improved the energy resolution. Some calculation and experimental details are discussed in this paper.

Hot-spot path in the thermodynamic space $({\rm \rho} R,T_{\rm i} )_{{\rm hs}} $ is investigated for direct-drive scaled-target family covering a huge interval of kinetic energy on both sides of kinetic threshold for ignition. Different peak implosion velocities and two initial aspect ratios have been considered. It is shown that hot spot follows almost the same path during deceleration up to stagnation whatever the target is. As attended, after stagnation, a clear distinction is done between non-, marginally-, or fully igniting targets. For the last, ionic temperature can reach very high values when the thermonuclear energy becomes very high.

Collisionless shocks are shocks in which the mean-free path is much larger than the shock front. They are ubiquitous in astrophysics and the object of much current attention as they are known to be excellent particle accelerators that could be the key to the cosmic rays enigma. While the scenario leading to the formation of a fluid shock is well known, less is known about the formation of a collisionless shock. We present theoretical and numerical results on the formation of such shocks when two relativistic and symmetric plasma shells (pair or electron/proton) collide. As the two shells start to interpenetrate, the overlapping region turns Weibel unstable. A key concept is the one of trapping time τp, which is the time when the turbulence in the central region has grown enough to trap the incoming flow. For the pair case, this time is simply the saturation time of the Weibel instability. For the electron/proton case, the filaments resulting from the growth of the electronic and protonic Weibel instabilities, need to grow further for the trapping time to be reached. In either case, the shock formation time is 2τp in two-dimensional (2D), and 3τp in 3D. Our results are successfully checked by particle-in-cell simulations and may help designing experiments aiming at producing such shocks in the laboratory.

A model is presented for the resonant third harmonic generation of short pulse lasers in cluster plasma in the presence of density ripple. Because of ripple in cluster density and plasma electron density outside the cluster, the phase-matching condition for the third harmonic process is satisfied, leading to resonant enhancement of harmonic generation. We explore the impact of laser intensity, cluster size, and collisional frequency of electrons on the efficiency of third harmonic generation. Moreover, since the group velocity of the third harmonic wave is greater than that of the fundamental wave, it causes the slippage of the generated harmonic pulse out of the fundamental laser pulse and its amplitude increases with time.

The dispersion relation of electrostatic waves with phase velocities smaller than the electron thermal velocity is investigated in relativistic temperature plasmas. The model equations are the electron relativistic collisionless hydrodynamic equations and the ion non-relativistic Vlasov equation, coupled to the Poisson equation. The complex frequency of electrostatic modes are calculated numerically as a function of the relevant parameters kλDe and ZTe/Ti where k is the wavenumber, λDe, the electron Debye length, Te and Ti the electron and ion temperature, and Z, the ion charge number. Useful analytic expressions of the real and imaginary parts of frequency are also proposed. The non-relativistic results established in the literature from the kinetic theory are recovered and the role of the relativistic effects on the dispersion and the damping rate of electrostatic modes is discussed. In particular, it is shown that in highly relativistic regime the electrostatic waves are strongly damped.

The paper presents the results of the analysis of the influence of curvature of the electrons trajectory in the anode–cathode gap of an ion diode on the diode impedance and evaluation of applicability of the one-dimensional (1D) Child–Langmuir (CL) ratio for calculation of the electron current. Investigations of an ion diode with a graphite anode in self-magnetic insulation mode were carried out. Experiments were performed on the TEMP-4M ion accelerator set in a double pulse mode, with the first negative pulse (150–200 kV, 400–600 ns) followed by the second positive pulse (250–300 kV, 150 ns). The result of this study is that we have determined the boundary conditions for the applicability of 1D CL ratio for calculation of the electron current in the ion diode. It was found that the deviation of the diode current–voltage characteristics from CL ratio will be observed only with a significant change in the acceleration voltage during electron drift or when the electron drift time exceeds the transit time of ions.

The effect of transverse magnetic fields on surface high-harmonic generation in intense laser–solid interactions is investigated. It is shown that the longitudinal motion of electrons can be coupled with the transverse motion via the magnetic fields, which lead to even-order harmonics under normal laser incidence. The dependence of the coupling efficiency and hence even harmonic generation with preplasma scale length and magnetic field strength are presented based upon particle-in-cell simulations. When the magnetic field is parallel to the laser electric field, the spectral intensity of the second harmonic is proportional to the magnetic field strength in a wide range up to 160 MG, while the situation with the magnetic field perpendicular to the laser electric field is more complicated. The second harmonic generation due to the magnetic field also tends to increase with the plasma density scale lengths, which is different from the high-harmonic generation by the oscillating mirror mechanism. With the increase of the laser spot size from a laser wavelength λL, both the magnetic field-induced harmonics and oscillating mirror high harmonics tend to increase first and then become saturated after 3λL. The magnetic field-induced second harmonic may be used to evaluate large self-generated magnetic fields developed near the critical density region and the preplasma conditions.

Results of experimental studies of the amplitude–temporal characteristics of a runaway electron (RE) beam, as well as breakdown voltage and discharge current with a picosecond time resolution are presented. The maximum pressure, at which a RE beam is detectable, decreases with increasing the voltage rise time. The waveforms of the discharge and RE beam currents are synchronized with those of the voltage pulses. It is shown that the amplitude–temporal characteristics of the RE beam depend on the designs of the gas-filled diode and cathode, as well as the gap length. The mechanism for the generation of REs in atmospheric-pressure gases is analyzed on the basis of the obtained experimental data.