## 1. Introduction

With the development of chirped-pulse amplification (CPA) and optical parametric CPA (OPCPA), extremely intense short pulse laser systems are increasingly applied to experimental investigations in ‘fast ignition’ inertial confinement fusion (ICF), laser–plasma interactions and strong field physics^{[Reference Kiriyama, Shimomura, Mori, Nakai, Tanoue, Kondo, Kanazawa, Pirozhkov, Esirkepov, Hayashi, Ogura, Kotaki, Suzuki, Daito, Okada, Kosuge, Fukuda, Nishiuchi, Kando, Bulanov, Nagashima, Yamagiwa, Kondo, Sugiyama, Bolton, Matsuoka and Kan1]}. The European Extreme Light Infrastructure (ELI) has designed a 25 PW, $10^{25}~\text{W cm}^{-2}$ single laser beamline (ILE)^{[Reference Chambaret, Canova, Martens, Cheriaux, Mourou, Cotel, Le Blanc, Druon, Georges, Forget, Ple and Pittman2]}, in which the high temporal contrast required by pre-plasma-free interactions^{[Reference Kaluza, Schreiber, Santala, Tsakiris, Eidmann, Meyer-ter-Vehn and Witte3]} will be one of the major challenges. Various factors resulting in temporal contrast degradation have been analyzed and verified, such as gain saturation^{[Reference Li, Zhang and Xu4]}, asymmetric spectral clipping^{[Reference Yang, Guo, Xie, Zhang, Sun, Gao, Li and Lin5]}, self-phase modulation ^{[Reference Kalashnikov, Andreev and Schönnagel6]}, amplified spontaneous emission (ASE) in the CPA system^{[Reference Kalashnikov, Risse, Schönnagel and Sandner7]} and pump distortion-induced noise (PDN)^{[Reference Forget, Cotel, Brambrink, Audebert, Le Blanc, Jullien and Chériaux8]}, surface-reflection-initiated pre-pulses^{[Reference Wang, Yuan, Ma, Wang, Xie and Qian9]} and parametric fluorescence (PF)^{[Reference Manzoni, Moses, Kärtner and Cerullo10]} in the OPCPA system, etc.

High energy laser systems usually enlarge the beam diameter to reduce the energy density on optical elements to mitigate the limited damage threshold problem, which inevitably leads to worse near-field quality, such as near-field intensity modulation and wavefront deviation. Means have existed for some time to analyze the influence on spatial output characteristics by wavefront deviation^{[Reference Guo, Zhu and Yang11]}. However, the discussed temporal characteristics in the wavefront deviation situation only cover the pulse width^{[Reference Liu, Cen and Li12]}, leaving the temporal contrast degradation by intensity modulation and wavefront deviation unclear.

In large aperture high power short pulse laser systems, the output is focused by an off-axis parabolic mirror (OAP), and this focusing propagation can be well simulated by a two-step coordinate transformed focusing method, not only taking advantage of using coordinate transformation but also avoiding the contradiction between the space scale and the diffraction limit^{[Reference Yang, Zhong and Ren13]}. Using this method, this paper simulates the output temporal contrast of a large aperture high power short pulse laser system. The temporal contrast degradation by intensity modulation and wavefront deviation is analyzed, the influence of the diameter on the temporal contrast degradation is revealed and the simulation and assumption results based on the parameters of the Shen Guang-II (SG-II) laser system are compared with the online experimental temporal contrast data. The conclusion can be drawn that the near-field quality is one of the possible factors hindering higher temporal contrast in large aperture high power short pulse laser systems.

## 2. Theory and simulation

### 2.1. Model

In practical short pulse applications, the amplified pulse by CPA or OPCPA will be focused onto the target using an OAP after being compressed by a grating compressor. The analogy can be made between the OAP focusing subsystem and an ideal lens, because it linearly reflects propagation in free space without dispersion, chromatic aberration or B-integral nonlinear effects, the spatial intensity distribution of which can be calculated by the Fresnel diffraction integral formula:

where $x$ and $y$ are the coordinates in the observation plane $(z=L)$, $x_{0}$ and $y_{0}$ are the coordinates in the incident plane, ${\it\lambda}=2{\rm\pi}c/{\it\omega}$ is the wavelength of the incident pulse, $k=2{\rm\pi}/{\it\lambda}$ is the wavevector and $f$ is the focal length of the OAP.

For higher spatial resolution and a more accurate result, on applying the coordinate transforms $x^{\prime }=x/(1-L/f)$, $y^{\prime }=y/(1-L/f)$ and $z^{\prime }=z/(1-L/f)$, Equation (1) becomes

Here, $f_{X}=x/{\it\lambda}z$ and $f_{Y}=y/{\it\lambda}z$ are the spatial frequencies in the $x$ and $y$ directions, respectively, and the spatial Fourier transform and its inverse form are introduced to reduce the computing time by the fast Fourier transform (FFT) algorithm:

When Equation (2) is applied to focusing propagation ($L=f$), a problem of infinite $x^{\prime }$, $y^{\prime }$ and $z^{\prime }$ is encountered, which can be ingeniously resolved by a two-step simulation, first calculating from $z=0$ to $z=L$ with coordinate transforms and finally calculating from $z=L$ to $z=f$ without coordinate transforms, in which angular spectrum diffraction theory is applied. The value of $L$ is related to the coordinate magnification, which determines whether there are sufficient effective sampling points in the focal spot area. At the same time, the value of $L$ determines whether the situation meets the requirement for the Fresnel approximation in step one with coordinate transforms. With this method in two spatial dimensions, the calculation of the focusing characteristics of a short wavelength laser beam through a short focus system has been successfully made^{[Reference Yang, Zhong and Ren13]}.

Thus, we can use an operator $P_{\text{space}}$ to represent the above focusing calculation:

The output field as a function of time is given by the inverse time Fourier transform:

Then, the temporal contrast can be derived by the temporal distribution of the normalized pulse intensity. In this paper, simulation of one dimension in space is used to demonstrate the influence of spatial near-field quality on temporal contrast. The idea of the two-step coordinate transformed space–time propagation simulation is shown in Figure 1.

### 2.2. Near-field intensity modulation

In large aperture high power laser systems, light beams are always shaped into flattened Gaussian beams (FGBs) for higher energy extraction efficiency and slower nonlinear growth. A FGB of order $N$ and spot size ${\it\omega}_{0}$ has the following normalized field distribution^{[Reference Bagini, Borghi, Gori, Pacileo, Santarsiero, Ambrosini and Spagnolo14]}:

where $r$ is the radial coordinate.

In practical applications, the real near-field intensity is always modulated rather than an ideal FGB due to multiple causes, such as diffraction and nonlinear growth. The fill factor (FF) of the laser beam is usually used to describe the uniformity of the modulated near-field intensity, which is defined as the ratio of the average intensity and the peak intensity in the flattened area:

### 2.3. Wavefront deviation

The ideal incident wavefront is plane; however, it can be easily deviated by phase shifts from the optical system and optical elements. Phase shifts from the optical system, such as thermal deformation and inhomogeneity of the medium, mainly contribute to dynamic wavefront deviation on a large spatial scale, which can be detected online by a Shack–Hartmann wavefront sensor. According to research carried out at Lawrence Livermore National Laboratory^{[Reference Wolfe and Lawson15, Reference Lawson, Auerbach and English16]}, the phase shifts from optical elements, as the main component of static wavefront deviation, can be divided into three spatial frequency regions: figure ($0{-}0.03~\text{mm}^{-1}$), waviness ($0.03{-}8~\text{mm}^{-1}$) and roughness (higher than $8~\text{mm}^{-1}$), which can be expressed by the improved noise model^{[Reference Chen, Zheng, Chen, He, Su, Chen and Yuan17]}:

where $C_{i}$ is the noise magnitude, $R_{i}(-1,1)$ is a random function distributed between $-1$ and 1, $\ast$ represents a convolution operation, $s_{ix}$ and $s_{iy}$ are the spatial scales in the $x$ and $y$ directions, respectively, determining the noise spatial length. The total wavefront deviation is the sum of noise of every spatial length. Common evaluation indices of wavefront deviation are the peak to valley (PV) value, root mean square (RMS), RMS gradient (GRMS) and power spectral density (PSD) curve.

## 3. Results and discussion

The parameters of the analyzed large aperture high power short pulse model are those of the SG-II laser system^{[Reference Xu, Wang, Li, Dai, Lin, Gu and Zhu18]}. A short pulse of 290 mm near-field diameter (ninth-order FGB) and 5 ps pulse width (Gaussian) at a center wavelength of 1053 nm is focused by an OAP of 800 mm focal length and 320 mm diameter. The two-step coordinate transformed focusing method is used, while the computer samples 2048 points in space from $-160$ to 160 mm and 2048 points in time from $-100$ to 100 ps. The value of $L$ is the key parameter of the simulation due to the two reasons mentioned above. Here, we use $L=0.9997f{-}0.9999f$ to reach the requirements of both high resolution and the Fresnel approximation, then we obtain the convergence results.

### 3.1. Temporal contrast degradation by intensity modulation

First, the influence of intensity modulation on temporal contrast is analyzed. The model applied here is a ninth-order FGB with a diameter of 290 mm, and its near-field intensity is modulated randomly in three situations for which the FFs are 0.4, 0.6 and 0.8, respectively, as shown in Figure 2(a). The temporal profile of the incident pulse is Gaussian type with a pulse width (FWHM) of 5 ps. No wavefront deviation is introduced during the focusing propagation.

Figure 2(b) shows the temporal contrast degradation by intensity modulation with different FFs. If there is no intensity modulation, the temporal contrast is better than $10^{9}$ (blue dashed line). A pre-pulse can be seen 50 ps before the main pulse when the near-field intensity is modulated, which leads to a poor temporal contrast. Moreover, the front edge of the main pulse becomes moderate 10–20 ps ahead of the peak point, which means that the pulse width of the main pulse will be broadened at the same time. When the intensity is modulated by a FF of 0.6, the temporal contrast will be degraded by more than one order of magnitude to worse than $10^{8}$ in the time window of 50 ps before the main pulse. The tendency can be seen from this figure that the smaller the FF of the modulated intensity is the worse the output temporal contrast is.

### 3.2. Temporal contrast degradation by wavefront deviation

Second, using the same model as described above, the effect of wavefront deviation on temporal contrast is discussed. At this time the incident intensity has no modulation, but various types of random wavefront deviations are introduced in the focusing propagation, which have different magnitudes and different spatial frequencies. As shown in Figure 3(a), these four wavefront deviations are illustrated in the form of PSD curves.

Figure 3(b) shows the temporal contrast degradation by different types of wavefront deviations. If there is no wavefront deviation, the temporal contrast will be better than $10^{9}$. Comparing the blue line $(\text{PV}=0.5{\it\lambda})$ with the green line $(\text{PV}=0.1{\it\lambda})$ in the same spatial frequency range (the same minimum spatial period $T$ is 40 mm, which means that the same maximum spatial frequency is $0.025~\text{mm}^{-1}$), the wavefront deviation with larger PV value will cause more severe temporal contrast degradation and in this spatial frequency range the temporal contrast will be two orders of magnitude worse than that of the ideal situation.

When the wavefront deviations have the same magnitude $(\text{PV}=0.1{\it\lambda})$ but different spatial frequency ranges, as shown in Figure 3(b), there are large differences between the temporal contrast degradations. The wavefront deviation with the highest maximum spatial frequency of $6.4~\text{mm}^{-1}$ (cyan line; minimum spatial period of 0.15625 mm) causes huge degradation of more than five orders of magnitude from $10^{9}$ to $10^{4}$, while the wavefront deviation with the lowest maximum spatial frequency of $0.025~\text{mm}^{-1}$ (green line; minimum spatial period of 40 mm) induces smaller degradation of one order. Interestingly, temporal contrast degradation by the wavefront deviation with the middle maximum spatial frequency of $1.6~\text{mm}^{-1}$ (red line; minimum spatial period of 0.625 mm) is quite notable within 40 ps before the main pulse; however, it becomes less significant beyond 40 ps before the main pulse. Therefore, higher spatial frequency contributes to temporal contrast degradation in a larger time window and for the same magnitude the wavefront deviation with low maximum spatial frequency has a slight influence on temporal contrast.

### 3.3. Influence of diameter on temporal contrast degradation

Third, the importance of discussing temporal contrast in a system with large diameter is demonstrated by analyzing the influence of diameter on temporal contrast degradation. In this case, introducing the same magnitude and same spatial frequency wavefront deviation as shown in Figures 4(a) and (b), we consider three models with different diameters: 320 (blue solid line), 160 (green solid line) and 80 mm (red solid line). The focal length is fixed at 800 mm.

The result can be seen in Figure 4(c). The temporal contrast for 80 mm diameter is three orders of magnitude better than that for 320 mm diameter and the model with larger diameter suffers more severe temporal contrast degradation. The temporal contrast is around $10^{10}$ in the model with 320 mm diameter without any near-field quality deterioration (cyan dashed line). However, for an ideal thin beam that does not need to be focused, the temporal contrast will be promoted to the computer simulation resolution limit, about $10^{30}$ (purple dotted line). Therefore, temporal contrast degradation deserves more attention in larger aperture high power short pulse laser systems.

### 3.4. Simulation of temporal contrast degradation compared with experimental measurement results

Finally, the simulation model is combined with the near-field quality measured in the established SG-II large aperture high power short pulse laser system. Figure 5(a) shows the experimental near-field intensity detected in the incident plane and Figure 5(b) shows the experimental intensity modulation (green solid line) and ideal spatial profile (blue dashed line) in one dimension. The calculated FF for the experimental near-field intensity modulation is 0.73. Figure 5(c) shows the experimental dynamic wavefront deviation detected online by a Shack–Hartmann wavefront sensor, which is unable to show waviness or roughness in the intermediate and higher frequency regions. Thus, in this case, no waviness and roughness are introduced into the simulation model. The experimental dynamic wavefront deviation in one dimension and its PSD curve are illustrated as the violet lines in Figures 5(d) and (e). After the propagation through the OAP focusing subsystem, the output temporal contrast can be obtained. In Figure 5(f), the cyan dashed line is the ideal output temporal contrast; the red solid line is the online experimental measurement result and the blue solid line is the calculation result. The online experimental measurement result^{[Reference Wang, Ouyang, Ma, Yuan, Xu and Qian19]} shows a degraded temporal contrast of $10^{4}$, while the calculation result shows a temporal contrast degradation of less than two orders of magnitude from $10^{9}$. There are two reasons why the two results differ. (1) In this model beam quality is considered as the only influencing factor; however, in real experiments there are many other characteristics, such as ASE, PF and spectral modulation. (2) Only near-field intensity and low spatial frequency wavefront deviation are introduced, without considering waviness and roughness.

To introduce the waviness and roughness into the simulation model, an assumption is made based on the SG-II short pulse laser system configuration. In consideration of the cutoff frequency of the spatial filter, a wavefront deviation of waviness and roughness whose maximum spatial frequency is $0.6~\text{mm}^{-1}$ is added, as shown in Figure 5(d) by the green line. The PV of total wavefront deviation is $2.5{\it\lambda}$, the RMS is $0.5{\it\lambda}$, the GRMS is $3.7{\it\lambda}~\text{cm}^{-1}$ and the PSD curve is shown in Figure 5(e) by the green line. These values are gradually adjusted since there are no online measurement data for the system. After considering the waviness and roughness, there is serious temporal contrast degradation of more than three orders of magnitude from the temporal contrast degradation level due to only intensity modulation and low spatial frequency wavefront deviation. In addition, the assumption result has a similar trend to the online experimental measurement result, in spite of the discrepancy due to some other factors. From this point of view, high spatial frequency wavefront deviations, such as waviness and roughness, influence the temporal contrast more seriously; therefore, waviness and roughness deserve careful attention, and methods of measurement and control of the wavefront deviation within the abovementioned spatial frequency region should be investigated further.

## 4. Conclusion

In this paper, a two-step coordinate transformed focusing method is used to simulate the output focused by an OAP in a large aperture high power short pulse laser system, and thus the output temporal contrast is calculated. The temporal contrast degradation by intensity modulation and wavefront deviation is analyzed. It is found that intensity modulation with a smaller FF and larger wavefront deviation causes more severe temporal contrast degradation, and this degradation is closely related to the spatial frequency of the total wavefront deviation. Then, the importance of discussing temporal contrast in a system with large diameter is demonstrated by analyzing the influence of diameter on temporal contrast degradation. Finally, with experimental intensity modulation, experimental low spatial frequency wavefront deviation and assumed higher spatial frequency wavefront deviation using the parameters of the SG-II laser system, the output temporal contrast is derived, which shows that waviness and roughness influence the temporal contrast more seriously. Therefore, the conclusion can be drawn that near-field quality deterioration might lead to temporal contrast degradation, hindering higher temporal contrast in large aperture high power short pulse laser systems.

Due to the lack of effective tools to measure intermediate and higher spatial frequency wavefront deviation, the need for online monitoring of the wavefront deviation in these frequency regions in large aperture high power short pulse laser systems is put forward. Improvement of the laser near-field quality might be a new method towards high temporal contrast in large aperture high power short pulse laser systems.

## Acknowledgement

The presented work is supported by the National High-Tech Committee of China (Grant no. 2013AA8044010).