2 Methodology

The fully relativistic particle-in-cell (PIC) code EPOCH^[ Reference Arber, Bennett, Brady, Lawrence-Douglas, Ramsay, Sircombe, Gillies, Evans, Schmitz, Bell and Ridgers ^{78 ]} was used in two and three dimensions to model the gamma ray generation in the interaction. Output data were generated individually for synchrotron and bremsstrahlung photons above 100 keV. The techniques used for calculating synchrotron emission in the EPOCH code are described in Ref. [Reference Ridgers, Kirk, Duclous, Blackburn, Brady, Bennett, Arber and Bell79], and those used for calculating bremsstrahlung emission are described in Refs. [Reference Vyskočil, Klimo and Weber56,Reference Wu, He, Yu and Fritzsche80]. For the 2D simulations, the spatial grid had dimensions of 30 μm × 20 μm ( $x\times y$ ) and cell size 8 nm × 12 nm, with free-space boundaries. To simulate an experimentally practical low Z target material, solid density plastic (CH) was chosen and modelled as a uniform, fully ionized plasma with electron density ${n}_{\mathrm{e}} = 3.5\times {10}^{29}\kern0.1em$ m⁻³, neutralized by an equal ratio of C⁶⁺ and H⁺ ions. The initial electron and ion temperatures were 3 keV and 100 eV, respectively, with 50 electron macroparticles per cell and 10 ion macroparticles per cell per species. The target parameters were the thickness, $l$ , and the angle between the normal to the target surface and the $x$ -axis, ${\theta}_{\mathrm{i}}$ . The laser pulse was incident upon the target from the minimum $x$ boundary, and had wavelength ${\lambda}_{\mathrm{L}} = 800\kern0.1em$ nm and a Gaussian temporal-intensity profile with full width at half maximum (FWHM) ${\tau}_{\mathrm{L}}$ . The laser pulse was focused at $x = {x}_{\mathrm{f}}$ , where $x = 0$ corresponds to the target irradiated (front) surface (the laser pulse focuses behind the target surface for ${x}_{\mathrm{f}}>0$ and in front of the target for ${x}_{\mathrm{f}}<0$ ) and ${x}_{\mathrm{f}}$ is the defocus parameter. At focus, the laser pulse has a Gaussian spatial-intensity profile with diameter ${\phi}_{\mathrm{L}}$ (FWHM). The values of $l$ , ${\theta}_{\mathrm{i}}$ , ${\tau}_{\mathrm{L}}$ , ${x}_{\mathrm{f}}$ and ${\phi}_{\mathrm{L}}$ used are stated in each section.

The 3D simulations used a spatial grid with dimensions 20 μm × 15 μm × 15 μm ( $x\times y\times z$ ) and cell size 10 nm $\times$ 30 nm $\times$ 30 nm. The target density was reduced to an electron density of ${n}_{\mathrm{e}} = 1.74\times {10}^{29}\kern0.1em$ m⁻³ due to the increased cell size, whilst maintaining the areal density of the quoted target thicknesses in Section 7 by modelling an increased thickness. There were 12 electron macroparticles per cell and 6 ion macroparticles per cell per species.

The BISHOP code was used in conjunction with EPOCH to automate the 2D PIC simulations in 2D grid scans of various parameters, and for the Bayesian optimization of various objectives using a GP regression algorithm^[ Reference Pedregosa, Varoquaux, Gramfort, Michel, Thirion, Grisel, Blondel, Prettenhofer, Weiss, Dubourg, Vanderplas, Passos, Cournapeau, Brucher, Perrot, Duchesnay and Mach ^{81 ]}, in the same manner as the code used in Ref. [Reference Dolier, King, Wilson, Gray and McKenna67]. For all optimization scans, 30 initial simulations are performed with randomly generated input parameters (10 more than used previously^[ Reference Dolier, King, Wilson, Gray and McKenna ^{67 ]} due to one more input parameter), the objective function is evaluated from the simulation data and the algorithm produces a probability distribution of all potential functions that could fit the results to create a model of the objective function. An acquisition function calculated from the model then determines the next set of input parameters to simulate. Acquisition functions corresponding to the upper confidence bound, expected improvement and probability of improvement methods are calculated and one of these is chosen at each iteration in a process known as hedging, which outperforms the use of individual acquisition functions in identifying the optimum in the minimum number of iterations^[ Reference Hoffman, Brochu and de Freitas ^{82 ]}. With each iteration the model is updated and the acquisition function guides the parameters towards the optimum of the objective function. Up to 200 simulations were run for each objective function used to identify their respective optimum. The optimization parameters were  ${\theta}_{\mathrm{i}}$ , log ${}_{10}$ $l$ , ${\tau}_{\mathrm{L}}$ , ${\phi}_{\mathrm{L}}$ and ${x}_{\mathrm{f}}$ , illustrated in Figure 1(a). Several of the objective functions used (defined in Sections 4 and 5), which correspond to important parameters of the synchrotron emission, are also shown: the synchrotron yield, peak angle-resolved yield and number of photons in the high-energy spectral tail for different cut-off values. The optimum parameters found in Sections 4 and 5 generally match what is expected from the 2D parameter space scans in Section 3. An example demonstration of Bayesian optimization of an arbitrary 1D function susceptible to noise is shown in Figure 1(b) for different numbers of iterations.

Figure 1 (a) The Bayesian optimization loop and schematic of the simulation setup. The synchrotron photon energy spectrum ( ${\mathrm{d}N}_{\mathrm{sy}}/ \mathrm{d}\varepsilon$ ) and angle-resolved yield ( $\mathrm{d}\sum {\varepsilon}_{\mathrm{sy}}/ \mathrm{d}\theta$ ) generated in each simulation are depicted to illustrate several of the objective functions. (b) An example of Bayesian optimization of a noisy 1D function showing the true function (black), the model (red) and the acquisition function (blue) for different numbers of iterations (n).

Previously, the same approach was applied for the optimization of laser-driven ion acceleration^[ Reference Dolier, King, Wilson, Gray and McKenna ^{67 ]} using a single physical parameter in the objective function (maximum ion energy). Here this method is applied both for single and multiple physical parameters in the objective function, enabling the exploration of different objective functions to influence the trade-off in one required beam parameter against another.

3 2D parameter space scans of gamma ray emission

Before optimizing the interaction, the influence of a number of input parameters, including each of the chosen optimization parameters, was first explored in 2D parameter space scans. The target thickness was one of the varied parameters in each scan, enabling separation of the effect of RSIT, which can be controlled with target thickness^[ Reference Goodman, King, Wilson, Gray and McKenna ^{83 ]}, from the effect of varying each of the other parameters. Initially, in Figures 2(a)–2(c) a pulse with ${\tau}_{\mathrm{L}} = 30\kern0.1em$ fs, ${\phi}_{\mathrm{L}} = 3\kern0.1em \unicode{x3bc}$ m and ${x}_{\mathrm{f}} = 0$ incident at target normal ( ${\theta}_{\mathrm{i}} = 0^{\circ}$ ) was considered, and simulated for peak laser intensities ${I}_{\mathrm{L}} = 3.16\times {10}^{21}$ – $3.16\times {10}^{23}\kern0.1em$ W cm⁻², covering a range for which synchrotron radiation may become measurable, up to where it is expected to dominate. Target thicknesses between 50 nm and 5 μm were simulated, encompassing relativistically transparent to opaque targets across the intensity range, demonstrated by the percentage laser energy transmission in Figure 2(a) and also shown by the white contours of these values in Figures 2(a)–2(c).

Figure 2 (a) Percentage transmission of the laser pulse, (b) total electron energy in front of the plasma critical surface and in the laser skin depth averaged over the period of synchrotron emission and (c) laser-to-synchrotron photon energy conversion efficiency, all for varying target thickness and laser intensity. (d)–(f) Laser-to-synchrotron photon energy conversion efficiency for varying pulse duration, focal spot size and defocus, respectively, with target thickness.

To produce synchrotron radiation in these interactions, highly relativistic electrons are required in an intense electromagnetic field, and as changes to the target thickness cause a transition between transparent and opaque targets, the coupling of the laser energy to the electrons also changes^[ Reference Williamson, Gray, King, Wilson, Dance, Armstrong, Rusby, Brabetz, Wagner, Zielbauer, Bagnoud, Neely and McKenna ^{84 ]}. Therefore, the total electron energy sampled along the $x$ -axis ( $\mid y\mid <6\kern0.1em$ nm) in the transparent plasma in front of the relativistic critical density surface ( ${x}_{\mathrm{c}}$ , where ${n}_{\mathrm{e}} = {\gamma}_{\mathrm{e}}{n}_{\mathrm{c}}$ ) and in the laser skin depth ( ${\delta}_{\mathrm{s}} = c/{\omega}_{\mathrm{pe}}$ ), integrated over the period of synchrotron emission (estimated as $-{\tau}_{\mathrm{L}}/2<t<{\tau}_{\mathrm{L}}$ , where the laser peak intensity reaches $x = 0$ at $t = 0$ ), was calculated for each of the simulations and is shown in Figure 2(b). The highest total electron energy in the plasma accessible to the laser field for each ${I}_{\mathrm{L}}$ simulated is on average found when the laser energy transmission is $13\%$ , where the points either side of the maximum provide transmission bounds of 2.6% and 29%, respectively, for the maximum. This agrees well with the laser-to-synchrotron radiation energy conversion efficiency, ${\eta}_{\mathrm{sy}}$ , values in Figure 2(c), which are maximized for 11% transmission, with the adjacent points providing a range between 1.9% and 26%, indicating that the improved coupling of energy to electrons in the laser field for thin targets that undergo RSIT enhances the synchrotron conversion efficiency.

3.4 Scaling of the synchrotron conversion efficiency

The 2D parameter space scans provide target thickness-dependent scalings of ${\eta}_{\mathrm{sy}}$ with each of the other parameters. In Figure 3, these are presented for peak laser intensity, pulse duration and focal spot size, with power law fit values given for the value of $l$ that maximizes ${\eta}_{\mathrm{sy}}$ (black line) and for an indicative opaque target (red line; $l = 5\kern0.22em \unicode{x3bc}$ m for Figure 3(a) and $l = 3\kern0.22em \unicode{x3bc}$ m for Figures 3(b) and 3(c)). In Figure 3(a), ${\eta}_{\mathrm{sy}}\propto {I}_{\mathrm{L}}^{1.8}$ for $l = 5\kern0.22em \unicode{x3bc}$ m and ${\eta}_{\mathrm{sy}}\propto {I}_{\mathrm{L}}^{1.4}$ for the optimum thicknesses, in agreement with the ${\eta}_{\mathrm{sy}}\propto {I}_{\mathrm{L}}^{1.5}$ scalings reported in Refs. [Reference Ji, Pukhov, Nerush, Kostyukov, Shen and Akli27,Reference Vyskočil, Gelfer and Klimo57]. In both cases, the rate of increase slows as ${I}_{\mathrm{L}}$ approaches ${10}^{23}\kern0.1em$ W cm⁻². The faster scaling for $l = 5\kern0.22em \unicode{x3bc}$ m and the convergence of the two lines of fit in Figure 3(a) result from the reduced role of target expansion in the optimization of the synchrotron emission with increasing laser intensity, due to the increasing importance of the radiation pressure and the increasing relativistically corrected critical density.

Figure 3 Scaling of the laser-to-synchrotron energy conversion efficiency with (a) peak laser intensity, (b) pulse duration and (c) focal spot FWHM, for varying target thickness. Power law fits are shown for the optimum target thicknesses (black) and for the thickest targets used (red; $l = 5\ \unicode{x3bc}\mathrm{m} $ for (a) and $l = 3\ \unicode{x3bc}\mathrm{m} $ for (b) and (c)).

The use of ultrathin targets can increase ${\eta}_{\mathrm{sy}}$ by almost an order of magnitude compared with $l = 5\kern0.22em \unicode{x3bc}$ m for the lowest intensities considered here, and therefore such targets provide the most accessible route to generating synchrotron radiation in laser–solid interactions. Synchrotron emission can be increased for the case of thick targets if the laser pulse interacts with a significant pre-plasma at the front surface, which can be produced by the laser light preceding the main pulse if the temporal-intensity contrast is low enough. This was not modelled in our simulations, which only included the main Gaussian peak in the temporal profile of the laser pulse. In experiments, increasing the peak laser intensity by orders of magnitude may require increasing the temporal-intensity contrast to prevent pre-expansion of the target or the creation of a significant pre-plasma at the front surface. Otherwise, the optimum conditions for synchrotron generation are expected to change.

For constant laser pulse energy, the pulse duration is found to play a weaker role, with ${\eta}_{\mathrm{sy}}\propto {\tau}_{\mathrm{L}}^{-0.8}$ for both $l = 3\kern0.22em \unicode{x3bc}$ m and the transparent optimum thicknesses in Figure 3(b). In contrast, the conversion efficiency depends strongly on ${\phi}_{\mathrm{L}}$ , with ${\eta}_{\mathrm{sy}}\propto {\phi}_{\mathrm{L}}^{-1.9}$ for $l = 3\kern0.22em \unicode{x3bc}$ m and ${\eta}_{\mathrm{sy}}\propto {\phi}_{\mathrm{L}}^{-1.3}$ for the optimum thicknesses in Figure 3(c). In a similar manner to the intensity scan, the transition from RESE-dominated synchrotron emission for ${\phi}_{\mathrm{L}} = 10\kern0.22em \unicode{x3bc}$ m to laser-injection-dominated emission for ${\phi}_{\mathrm{L}} = 1\kern0.22em \unicode{x3bc}$ m produces these different scalings. The faster ${I}_{\mathrm{L}}\propto {\phi}_{\mathrm{L}}^{-2}$ dependence of the pulse intensity on the focal spot size in 3D for constant pulse energy may provide a different ${\eta}_{\mathrm{sy}}$ scaling with ${\phi}_{\mathrm{L}}$ .

3.5 Angle-of-incidence

The results for the final optimization parameter, the angle-of-incidence of the laser pulse on the target, are shown in Figure 4. In Figure 4(a), the synchrotron conversion efficiency is maximized for ${\theta}_{\mathrm{i}} = 45{}^{\circ}$ and $l = 216\kern0.1em$ nm, where it is 60% higher than the maximum value reached for normal incidence. Similar to the other parameters, ${\eta}_{\mathrm{sy}}$ is highest for targets that on average transmit $14\%$ of the laser light. To explore why the angle-of-incidence improves ${\eta}_{\mathrm{sy}}$ , the electron energy spectra and ${\chi}_{\mathrm{e}}$ spectra (both averaged over the period of synchrotron emission) are plotted in Figures 4(b) and 4(c), respectively, for additional simulations using $l = 200\kern0.1em$ nm and ${\theta}_{\mathrm{i}} = 0{}^{\circ}$ , 45°. The fast electron population contains 36% more energy for ${\theta}_{\mathrm{i}} = 45{}^{\circ}$ compared with 0°, because the p-polarized light improves energy coupling to electrons, and a much larger number of electrons with energy greater than 100 MeV are produced. This contributes to the enhanced spectrum of ${\chi}_{\mathrm{e}}$ values, increasing the generation of synchrotron radiation.

Figure 4 (a) Laser-to-synchrotron photon energy conversion efficiency for varying angle-of-incidence and target thickness. (b) Electron spectra, sampled over the whole simulation space, averaged over the period of synchrotron emission for a 200 nm foil at normal and 45° incidence, and (c) the corresponding time-averaged ${\chi}_{\mathrm{e}}$ spectra.

3.6 Bremsstrahlung emission for varied laser intensity and target thickness

Until now, the discussion has focused on the generation of synchrotron radiation and how this depends on key laser and plasma parameters. Gamma radiation will, however, also be produced via bremsstrahlung emission in these interactions. Distinguishing between these two photon sources is important for the design of experiments that aim to investigate either mechanism.

In this investigation, thin foils of a low Z material have been selected to minimize the production of bremsstrahlung radiation. In Figure 5(a), the conversion efficiency to bremsstrahlung radiation, ${\eta}_{\mathrm{br}}$ , is shown for different values of the target thickness and peak laser intensity. Generally, the conversion efficiency is observed to increase with each of these parameters, except where the laser energy starts to be transmitted. In such targets the bremsstrahlung emission is expected to decrease due to the increased loss of fast electrons and the reduced absorption of the pulse into electrons as the transmission becomes large.

Figure 5 (a) Laser-to-bremsstrahlung radiation energy conversion efficiency for varying laser intensity and target thickness. (b) Energy spectra for bremsstrahlung photons (solid) and synchrotron photons (dotted) for different target thicknesses. (c) The rate of energy conversion to bremsstrahlung radiation.

In Figure 5(b), the energy spectra of the bremsstrahlung and synchrotron radiation are shown for various target thicknesses, where ${I}_{\mathrm{L}} = {10}^{22}\kern0.1em$ W cm⁻². In all cases, the synchrotron radiation dominates by many orders of magnitude. The conversion efficiency to synchrotron radiation scales very quickly with peak laser intensity, and therefore the difference is expected to become even greater for higher intensities. The bremsstrahlung emission can take place for many picoseconds in laser–solid interactions^[ Reference Morris, Robinson and Ridgers ^{58 ]}, due to the presence of a hot electron population trapped inside the target^[ Reference Nilson, Davies, Theobald, Jaanimagi, Mileham, Jungquist, Stoeckl, Begishev, Solodov, Myatt, Zuegel, Sangster, Betti and Meyerhofer ^{87 ]}. However, in Figure 5(c), the rate of energy conversion to bremsstrahlung radiation, $\mathrm{d}\sum {\varepsilon}_{\mathrm{br}}/ \mathrm{d}t$ , where ${\varepsilon}_{\mathrm{br}}$ is the bremsstrahlung photon energy, for different target thicknesses is shown to quickly reduce after the peak of the laser temporal profile reaches the target at $t = 0$ . Only for the thickest target tested, $l = 3\kern0.22em \unicode{x3bc}$ m does the emission rate remain noticeable at the end of the simulation, and extrapolating this for a further 6 ps (the longest duration of ${K}_{\alpha }$ emission measured in Ref. [Reference Nilson, Davies, Theobald, Jaanimagi, Mileham, Jungquist, Stoeckl, Begishev, Solodov, Myatt, Zuegel, Sangster, Betti and Meyerhofer87]) increases ${\eta}_{\mathrm{br}}$ from $1.0\times {10}^{-3}\%$ to $8.6\times {10}^{-3}\%$ . This is still dominated by the corresponding value of ${\eta}_{\mathrm{sy}} = 0.27\%$ . Our modelling captures the brightest period of bremsstrahlung emission, enabling us to find parameters that minimize its production whilst maximizing properties of the synchrotron emission, informing the design of experiments to generate the purest source and clearest signature of synchrotron radiation.

4 Application of Bayesian optimization

The scans discussed thus far vary only two input parameters at a time, over a $15\times 15$ grid corresponding to 225 simulations. To perform a systematic grid scan of this type with 15 values for each of the five input parameters under consideration would require ${15}^5 = \mathrm{759,375}$ simulations, which was not feasible with the computational resources available. Instead, GP regression was used to identify the values of the input parameters that maximize a given objective function within 200 simulations.

The parameter search space was defined as $0{}^{\circ}\le {\theta}_{\mathrm{i}}\le 70{}^{\circ}$ , $50\kern0.1em$ nm $\le l\le 10\kern0.22em \unicode{x3bc}$ m, $30\kern0.1em$ fs $\le {\tau}_{\mathrm{L}}\le 100\kern0.1em$ fs, $1\kern0.22em \unicode{x3bc}$ m $\le {\phi}_{\mathrm{L}}\le 6\kern0.22em \unicode{x3bc}$ m and $-50\kern0.22em \unicode{x3bc}$ m $\le {x}_{\mathrm{f}}\le 50\kern0.22em \unicode{x3bc}$ m. The laser energy was kept constant in the same way as in the previous section, by setting ${I}_{\mathrm{L}} = 3\times {10}^{22}(30\kern0.1em$ fs $/{\tau}_{\mathrm{L}})(1\kern0.22em \unicode{x3bc}$ m $/{\phi}_{\mathrm{L}})\kern0.1em$ W cm⁻². The objective functions used are shown in Table 1, with the values of the input parameters for their optima. The objective functions that include only one physical property of the photon emission are as follows: the total energy of synchrotron emission ( ${f}_{{\mathrm{O}}1}$ ), where ${\varepsilon}_{\mathrm{sy}}$ is the synchrotron photon energy; the maximum of the angle-resolved energy of synchrotron emission, $\mathrm{d}\sum {\varepsilon}_{\mathrm{sy}}/ \mathrm{d}\theta$ ( ${f}_{{\mathrm{O}}2}$ ), where $\theta$ is the angle in the $xy$ -plane from the positive $x$ -axis; and the total number of synchrotron photons, ${N}_{\mathrm{sy}}$ , with energy exceeding 10 MeV ( ${f}_{{\mathrm{O}}3}$ ).

Table 1 The objective functions maximized with Bayesian optimization and the parameters of the found optimum for each.

To optimize multiple properties of the photon emission, a series of different objective functions were considered. These include the ratio of the peak angle-resolved energy of synchrotron emission to the total bremsstrahlung radiation energy ( ${f}_{{\mathrm{M}}1}$ ), and the same function multiplied by the acceptance function $A\left({f}_{{\mathrm{O}}2}\right)$   $\left({f}_{{\mathrm{M}}2}\right)$ . The acceptance function is defined as $A(f) = 1/\left(1+\exp \left(-\left(20/{f}_{\mathrm{max}}\right)\left(f-0.5{f}_{{\mathrm{max}}}\right)\right)\right)$ , where $f$ is a given objective function and ${f}_{{\mathrm{max}}}$ is the maximum value of the objective function found separately in a previous optimization scan. The acceptance function strongly reduces the value of the overall function for $f<{f}_{{\mathrm{max}}}/2$ to guide the optimization towards results where $f>{f}_{{\mathrm{max}}}/2$ , and it is plotted in the Supplementary Materials. In addition, the product of the peak angle-resolved energy of synchrotron emission and the number of synchrotron photons above 10 MeV ( ${f}_{{\mathrm{M}}3}$ ), and this function divided by the total bremsstrahlung radiation energy ( ${f}_{{\mathrm{M}}4}$ ), were optimized.

4.1 Optimization of individual synchrotron emission properties

The resulting synchrotron and bremsstrahlung photon energy spectra and the angle-resolved energy of synchrotron radiation are shown in Figures 6(a)–6(c), respectively, for the optimum of each objective function. In maximizing the total energy of synchrotron radiation with ${f}_{{\mathrm{O}}1}$ , the synchrotron spectrum with the highest number of photons less than 2 MeV is produced; such photons require less electron energy to be generated and thus can be generated in such high numbers to dominate the total synchrotron energy, accounting for $54\%$ for ${f}_{{\mathrm{O}}1}$ . This also corresponds to the only optimum where the pulse duration is greater than the minimum of 30 fs, at 42.8 fs. In optimizing the maximum of $\mathrm{d}\sum {\varepsilon}_{\mathrm{sy}}/ \mathrm{d}\theta$ and the number of synchrotron photons with energy greater than 10 MeV, similar peaks in $\mathrm{d}\sum {\varepsilon}_{\mathrm{sy}}/ \mathrm{d}\theta$ at $\theta \approx -50{}^{\circ}$ are found in Figure 6(c). The angular profiles differ mainly by the energy emitted in the backward direction ( $\mid \theta \mid >90{}^{\circ}$ ) due to the dependency of the RESE mechanism on target thickness, enhanced for ${f}_{{\mathrm{O}}3}$ . The optimum for ${f}_{{\mathrm{O}}1}$ only reaches approximately half of the peak value of $\mathrm{d}\sum {\varepsilon}_{\mathrm{sy}}/ \mathrm{d}\theta$ compared with ${f}_{{\mathrm{O}}2}$ and ${f}_{{\mathrm{O}}3}$ . The only sizeable difference in the parameters here, besides the small increase in ${\tau}_{\mathrm{L}}$ , is the change from ${\theta}_{\mathrm{i}} = 41.1{}^{\circ}$ to ${\theta}_{\mathrm{i}}\approx 25{}^{\circ}$ . The production of more directional synchrotron emission with changes to the angle-of-incidence has been reported in Refs. [Reference Serebryakov and Nerush31,Reference Zhao, Liu, Li and Xia46,Reference Nerush, Kostyukov, Ji and Pukhov88], in which a single lobe structure is also reported.

Figure 6 Synchrotron and bremsstrahlung radiation for the objective function optima in Table 1, for which ${I}_{\mathrm{L}} = 3\times {10}^{22}(30\kern0.1em$ fs $/{\tau}_{\mathrm{L}})(1\kern0.22em \unicode{x3bc}$ m $/{\phi}_{\mathrm{L}})\kern0.1em$ W cm⁻². (a) Synchrotron photon energy spectra, (b) bremsstrahlung photon energy spectra and (c) angular profiles of total emitted synchrotron photon energy.

5 Optimization in the synchrotron dominated regime

The optimization results discussed thus far are for conditions in which the synchrotron emission is relatively weak. The laser energy was increased by a factor of 10, corresponding to ${I}_{\mathrm{L}} = 3\times {10}^{23}(30\kern0.1em$ fs $/{\tau}_{\mathrm{L}}\left)\right(1\kern0.22em \unicode{x3bc}$ m $/{\phi}_{\mathrm{L}})\kern0.1em$ W cm⁻² for which the synchrotron emission is expected to become extremely strong, and further optimization scans were performed to identify the impact on the optimum input parameters. The objective functions ${f}_{{\mathrm{O}}1}$ , ${f}_{{\mathrm{O}}2}$ , ${f}_{{\mathrm{M}}1}$ and ${f}_{{\mathrm{M}}2}$ were used again. Instead of ${f}_{{\mathrm{O}}3}$ , the total number of synchrotron photons with energy above 50 MeV ( ${f}_{{\mathrm{O}}4}$ ) is optimized, due to the higher energy photon spectra produced under these conditions. Similar optimum parameters were obtained, as listed in Table 2, with the highest intensity combination of ${\tau}_{\mathrm{L}}$ , ${\phi}_{\mathrm{L}}$ and ${x}_{\mathrm{f}}$ values for most objective functions. An optimum defocus of close to 2 μm is found for objective functions ${f}_{{\mathrm{O}}1}$ , ${f}_{{\mathrm{O}}2}$ and ${f}_{{\mathrm{O}}4}$ , several times greater than the lower intensity cases in Table 1. This is due to the increased hole boring velocity, induced by the higher laser intensity, leading to greater recession of the critical surface at the target front side.

Table 2 The objective functions used for optimization with laser intensity of $3\times {10}^{23}(30\kern0.1em$ fs $/{\tau}_{\mathrm{L}}\left)\right(1\kern0.22em \unicode{x3bc}$ m $/{\phi}_{\mathrm{L}})\kern0.1em$ W cm⁻², and the parameters of the found optima.

5.1 Angle-resolved synchrotron emission

The optimum of ${f}_{{\mathrm{O}}2}$ corresponds to the maximum pulse duration of $100\kern0.1em$ fs. Similar peak values of $\mathrm{d}\sum {\varepsilon}_{\mathrm{sy}}/ \mathrm{d}\theta$ were obtained across the permitted range of values for the pulse duration. However, the direction of peak emission moves from $\theta \approx -50{}^{\circ}$ for ${\tau}_{\mathrm{L}}\approx 30\kern0.1em$ fs to $\theta \approx {0}^{\circ }$ for ${\tau}_{\mathrm{L}}\approx 100\kern0.1em$ fs. Emission in the same direction as the laser pulse propagation ( $\theta = 0{}^{\circ}$ ) has previously been associated with skin depth emission from the electrons near the front of the laser pulse in the reflected light^[ Reference Brady, Ridgers, Arber and Bell ^{26 ,} Reference Brady, Ridgers, Arber and Bell ^{62 ]}. However, the peak in the angle-resolved emission for the optimum of ${f}_{{\mathrm{O}}2}$ is a result of the stronger deformation of the target for the longer pulse duration. A deeper cavity in the target is formed, and electrons injected into the laser pulse from the edge of the cavity on the $y>0$ side of the focal spot by the positive half-cycle of the laser electric field are caught up by the next half-cycle of the laser pulse as they propagate. The change in the direction of the laser electric field experienced by the injected electrons forces them to propagate close to parallel with the laser pulse, emitting in this direction. In comparison with the lower laser intensity case in Figure 6(c), with the same normalization constant, a $100\times $ enhancement in the peak value is obtained for only a $10 \times $ higher energy laser pulse. Although the maximum synchrotron photon energies obtained for this optimum set of parameters are similar to the other spectra in Figure 7(a), the spectrum is shifted to lower photon energies, containing the most radiation below $2\kern0.1em$ MeV.

Figure 7 Synchrotron and bremsstrahlung radiation for the objective function optima in Table 2, for which ${I}_{\mathrm{L}} = 3\times {10}^{23}(30\kern0.1em$ fs $/{\tau}_{\mathrm{L}})(1\kern0.22em \mu$ m $/{\phi}_{\mathrm{L}})\kern0.1em$ W cm⁻². (a) Synchrotron photon energy spectra, (b) bremsstrahlung photon energy spectra and (c) angular profiles of total emitted synchrotron photon energy.

Optimization of total synchrotron emission with ${f}_{{\mathrm{O}}1}$ produces a peak in the angle-resolved synchrotron emission of similar magnitude to the optimum for ${f}_{{\mathrm{O}}2}$ in Figure 7(c), inadvertently maximizing both objectives. This optimum exhibits a pronounced double-peaked structure in the angular emission characteristic of synchrotron radiation from solid targets^[ Reference Brady, Ridgers, Arber and Bell ^{26 –} Reference Hadjisolomou, Jeong and Bulanov ^{28 ,} Reference Stark, Toncian and Arefiev ^{32 ,} Reference Vyskočil, Gelfer and Klimo ^{57 ,} Reference Brady, Ridgers, Arber and Bell ^{62 ]}. However, the 38° angle-of-incidence has resulted in the suppression of the lobe in the direction closer to the perpendicular of the rotated target surface and the enhancement of the lobe closer to the parallel. This behaviour is also present in Figure 6(c), although the weakened lobe is less discernible due to the more comparable backward-directed radiation. The effect of varying the angle-of-incidence on the synchrotron angular profile is discussed in more detail in the following sections. The optimum for ${f}_{{\mathrm{O}}1}$ also produces a near-identical synchrotron spectrum to that for ${f}_{{\mathrm{O}}4}$ , where the high-energy spectral tail is optimized with a thicker target and normal incidence.

6 Angle-of-incidence dependence of the forward synchrotron emission

Through applying Bayesian optimization, the angle-of-incidence has been identified as a critical parameter influencing the spatial profile of the synchrotron emission. This was investigated in more detail using laser parameters ${I}_{\mathrm{L}} = 3\times {10}^{22}\kern0.1em$ W cm⁻², ${\tau}_{\mathrm{L}} = 30\kern0.1em$ fs, ${\phi}_{\mathrm{L}} = 1\kern0.22em \unicode{x3bc}$ m and ${x}_{\mathrm{f}} = 0$ (close to most of the optima) and $l = 3\kern0.22em \unicode{x3bc}$ m (an arbitrary thickness opaque target).

In Figure 8(a), the maximum value of $\mathrm{d}\sum {\varepsilon}_{\mathrm{sy}}/ \mathrm{d}\theta$ is shown for different directions of the emitted synchrotron photons, ${\theta}_{\mathrm{sy}}$ , as a function of the angle-of-incidence. The photons are grouped into those propagating with direction $90{}^{\circ}>{\theta}_{\mathrm{sy}}>0{}^{\circ}$ and $0{}^{\circ}>{\theta}_{\mathrm{sy}}>-90{}^{\circ}$ , corresponding to the directions in which the two lobes are usually produced due to laser-injected emission; these directions are hereafter referred to as ${\theta}_{90,0}$ and ${\theta}_{0,-90}$ for brevity. The optimization results from Table 1 and Figure 6(c) are added as data points. The observed lobes are of similar magnitude for normal incidence. As the target is rotated to give the normal to the front surface direction ${\theta}_{\mathrm{n}} = {\theta}_{\mathrm{i}}-180{}^{\circ}$ , the magnitude of the lobe propagating in ${\theta}_{0,-90}$ (closer to parallel with the target surface) is enhanced by more than a factor of two when it peaks for ${\theta}_{\mathrm{i}} = 22.5{}^{\circ}$ . Many of the Bayesian optimization results are for similar values of ${\theta}_{\mathrm{i}}$ , albeit with greater directional synchrotron emission achieved due to changes to some of the parameters and many iterations. At the same time the magnitude of the ${\theta}_{90,0}$ -directed emission is reduced, and this begins to reverse for ${\theta}_{\mathrm{i}} = 37.5{}^{\circ}$ –52.5° until the emission in both directions shown reduces in magnitude and begins to converge for increasing ${\theta}_{\mathrm{i}}$ .

Figure 8 (a) Maximum value of $\mathrm{d}\sum {\varepsilon}_{\mathrm{sy}}/ \mathrm{d}\theta$ as a function of the angle-of-incidence for synchrotron photons emitted in angular ranges ${\theta}_{90,0}$ (black) and ${\theta}_{0,-90}$ (blue), where $l = 3\kern0.22em \unicode{x3bc}$ m, ${I}_{\mathrm{L}} = 3\times {10}^{22}\kern0.1em$ W cm⁻², ${\phi}_{\mathrm{L}} = 1\kern0.22em \unicode{x3bc}$ m, ${\tau}_{\mathrm{L}} = 30\kern0.1em$ fs and ${x}_{\mathrm{f}} = 0$ . The optima in Figure 6 are also shown (diamonds). (b) Total energy in electrons more than 10 MeV in a local intensity more than 10²¹ W cm⁻² propagating with angle ${\theta}_{\mathrm{e}}$ in the ranges ${\theta}_{90,0}$ (dashed) and ${\theta}_{0,-90}$ (solid) averaged over the period of synchrotron emission. (c) Energy-weighted mean angle between the electron trajectory and the propagation direction of the local electromagnetic field (left-hand axis) and mean electron quantum parameter (right-hand axis) for each group of electrons in (b). (d)–(f) The electron density for ${\theta}_{\mathrm{i}} = 0{}^{\circ}$ , 22.5° and 60°, respectively, where the total momentum of fast electrons (arrows) and the $I = {10}^{21}$ W cm⁻² contour (red) are also shown.

The synchrotron radiation is caused by highly relativistic electrons within the laser field, and the direction of the emitted radiation is predominantly in the electron direction of motion, ${\theta}_{\mathrm{e}}$ . The effect of varying the angle-of-incidence is examined by considering the electrons with energy greater than 10 MeV in local electromagnetic fields with intensity greater than ${10}^{21}\kern0.1em$ W cm⁻². The total energy of such electrons propagating in the same direction as each of the synchrotron lobes is averaged over the period of synchrotron emission and is shown in Figure 8(b). For normal incidence, the total electron energy in each direction is equal, and as the target is rotated, the electron population that propagates in the ${\theta}_{0,-90}$ direction rapidly acquires several times more total energy and continues to gain energy until ${\theta}_{\mathrm{i}} = 60{}^{\circ}$ . By contrast, the total energy of the electron population that propagates in the ${\theta}_{90,0}$ direction over the interaction changes very little with ${\theta}_{\mathrm{i}}$ , with a minimum for ${\theta}_{\mathrm{i}}\approx 30{}^{\circ}$ and a second peak at ${\theta}_{\mathrm{i}}\approx 52.5{}^{\circ}$ . Although these changes match some of those for the peak angle-resolved synchrotron emission in Figure 8(a), they do not explain why the ${\theta}_{0,-90}$ -directed emission peaks for ${\theta}_{\mathrm{i}} = 22.5{}^{\circ}$ and why the peak emission in both directions reduces for ${\theta}_{\mathrm{i}}$ approaching 70°.

The equation for ${\chi}_{\mathrm{e}}$ indicates that the synchrotron emission increases for increasing ${\gamma}_{\mathrm{e}}$ and for increasing ${\mathbf{E}}_{\perp }+{\mathbf{v}}_{\mathrm{e}}\times \mathbf{B}$ . The orientation of the electron motion to the fields, $\Delta \theta$ , is therefore important to consider. Here, ${\chi}_{\mathrm{e}}$ is maximized for antiparallel propagation $\Delta \theta = {180}^{\circ }$ and minimized for parallel propagation $\Delta \theta = 0^{\circ }$ . Figure 8(c) shows the energy-weighted average absolute angle between the direction of motion of the electrons and the local electromagnetic field, calculated from the Poynting vector $\mathbf{S} = \mathbf{E}\times \mathbf{B}/{\mu}_0$ . For the ${\theta}_{0,-90}$ -propagating electron population, the average angle between the electrons and the field exhibits an approximately linear decrease from ${\theta}_{\mathrm{i}} = 0^{\circ}$ to 70°. For the same field magnitude and ${\gamma}_{\mathrm{e}}$ , this would reduce the synchrotron radiation. The energy-weighted average value of ${\chi}_{\mathrm{e}}$ for the electrons propagating in each direction is also shown in Figure 8(c), where it can be seen that ${\chi}_{{\mathrm{e}}, {\mathrm{avg}}}$ changes very little for ${\theta}_{\mathrm{i}} = 0^{\circ}$ –30°, but is strongly reduced for greater ${\theta}_{\mathrm{i}}$ . The optimization of the magnitude of the angle-resolved synchrotron emission in ${\theta}_{0,-90}$ is a result of the balance of the greater energy in the population of electrons propagating in the same direction with the reduced average emission parameter. For the ${\theta}_{90,0}$ -directed electron population, the angle to the field $\Delta \theta$ and the average emission parameter ${\chi}_{{\mathrm{e}}, {\mathrm{avg}}}$ only begin to strongly reduce for ${\theta}_{\mathrm{i}}>52.5{}^{\circ}$ , which reflects the behaviour of the associated synchrotron emission in Figure 8(a).

Example snapshots of the electron density on the simulation grid are shown in Figures 8(d)–8(f) close to the time of peak synchrotron emission for ${\theta}_{\mathrm{i}} = 0^{\circ}$ , 22.5° and 60°, respectively. The red contours correspond to $I = {10}^{21}\kern0.1em$ W cm⁻², and the arrows show the total momentum of fast electrons in each 150 nm × 150 nm region. For normal incidence in Figure 8(d), in every half laser cycle the electric field pulls electrons from one side of the focal spot towards the centre, with the electrons periodically coming from the different sides of the focal spot due to the oscillating direction of the laser electric field. This results in periodic synchrotron radiation at the laser frequency with the lobes separated by half a laser cycle, as observed in Ref. [Reference Hadjisolomou, Jeong and Bulanov28]. The shallow depth of the front surface plasma cavity means that only electrons propagating at a steep angle with the incoming laser pulse move into and experience the highest fields.

When the target is rotated to ${\theta}_{\mathrm{i}} = 22.5{}^{\circ}$ , as in Figure 8(b), the laser fields can pull electrons from a greater area and along the surface of the target further forwards with the pulse. This enhances the coupling of laser energy to ${\theta}_{0,-90}$ -propagating electrons, and as a result enhances the synchrotron emission in this direction. When ${\theta}_{\mathrm{i}}$ is increased further to 60°, as shown in Figure 8(c), the laser pulse is no longer reflected back in the $\mid \theta \mid >90{}^{\circ}$ direction, but along the surface of the target, reducing enhancement of the synchrotron radiation from counter-propagation of the reflected light with the laser-injected electrons propagating along the target surface. These electrons move with the laser pulse, which reduces their values of ${\chi}_{\mathrm{e}}$ , as shown in Figure 8(c), and causes them to produce less synchrotron radiation.

The Bayesian optimization results in Table 2 and Figure 7(c) indicate a laser intensity dependence of the optimum angle-of-incidence for producing the highest peak angle-resolved synchrotron emission, with the best results for ${{\theta}_{\mathrm{i}}\approx 42{}^{\circ}}$ in comparison with ${\theta}_{\mathrm{i}}\approx 25{}^{\circ}$ for the lower laser intensity case. The hole boring velocity increases for increasing laser intensity, causing the formation of a deeper cavity in the target and changing the evolution of the geometry of the interaction, and hence the optimum value of ${\theta}_{\mathrm{i}}$ . The hole boring velocity also reduces for higher densities. The optimum value of ${\theta}_{\mathrm{i}}$ is therefore also expected to be dependent on the target density.

7 Spatial control of synchrotron emission in 3D simulations

The investigation of the angle-of-incidence, which is highly influential on the spatial profile of synchrotron emission, identified through the optimization scans and explored in 2D in the previous section, was extended with 3D simulations. The synchrotron emission is dominated by electrons accelerated and injected from the edges further into the laser spatial profile due to interaction with the laser electric field, and therefore the spatial profile of synchrotron emission should change as the polarization of the laser light changes. This motivates simulations in which we separately test the effect of p-, s-, left-hand c- and right-hand c- polarized laser light on the spatial profile of synchrotron emission for varying angle-of-incidence.

A lower plasma density of $100\kern0.1em {n}_{\mathrm{c}}$ was simulated due to the high computational resources required to model higher plasma densities accurately in 3D. This may produce differences in the variation of the laser-injected synchrotron emission with the angle-of-incidence, as discussed at the end of the previous section. A single target thickness of 1 μm was chosen, due to the increased resources required to simulate thicker targets. Therefore, the synchrotron emission was not optimized with target thickness (or transparency) for each angle-of-incidence tested. Transparency occurs for each polarization state at normal incidence, eventually becoming opaque as the angle-of-incidence is increased, and is polarization dependent. The laser parameters were ${{I}_{\mathrm{L}} = 3\times {10}^{22}}$  W cm⁻², ${\tau}_{\mathrm{L}} = 30\kern0.1em$ fs, ${\phi}_{\mathrm{L}} = 1\kern0.22em \unicode{x3bc}$ m and ${x}_{\mathrm{f}} = 0$ , corresponding to a pulse energy of 10.9 J.

The spatial profiles of synchrotron emission are shown as a function of $\theta$ , the azimuthal angle from the positive $x$ -axis in the $xy$ plane where $-180{}^{\circ}\le \theta \le 180{}^{\circ}$ , and $\phi$ , the polar angle from the positive $z$ -axis where $0^{\circ}\le \phi \le 180{}^{\circ}$ . The angle-of-incidence corresponds to rotation in the $xy$ plane as in the 2D simulations.

Figure 9(a) shows the maximum angle-resolved energy of synchrotron emission of the two azimuthally defined regions corresponding to ${\theta}_{90,0}$ and ${\theta}_{0,-90}$ , as in the 2D results, as a function of the angle-of-incidence for p-polarization. In a similar manner to the 2D results, the magnitude of the emission in each direction for near-normal incidence is approximately equal. For increasing values of ${\theta}_{\mathrm{i}}$ , the magnitude of the ${\theta}_{90,0}$ -directed emission drops very rapidly, and the magnitude of the ${\theta}_{0,-90}$ -directed emission achieves a peak at ${\theta}_{\mathrm{i}} = 52.5{}^{\circ}$ that is only approximately 0.5 J sr⁻¹ higher than the result for normal incidence, less than the factor of two improvement observed in 2D. The magnitude of the angle-resolved emission in each direction appears to converge around ${\theta}_{\mathrm{i}} = 60{}^{\circ}$ , before the magnitude of the ${\theta}_{0,-90}$ lobe unexpectedly peaks again for ${\theta}_{\mathrm{i}} = 67.5{}^{\circ}$ ; such behaviour was not observed in the 2D simulations shown in Figure 8(a). The conversion efficiency for p-polarization is shown in Figure 9(d), where the overall synchrotron conversion efficiency peaks for ${\theta}_{\mathrm{i}} = 45{}^{\circ}$ while the conversion efficiency to ${\theta}_{0,-90}$ -directed emission peaks at approximately twice its normal incidence value.

For an s- (along $z$ ) polarized laser pulse, the synchrotron lobes are oriented perpendicular to those for p-polarization because they are generated along the axis of polarization of the laser light. Therefore, each lobe is defined by their polar angle, $\phi$ . In Figure 9(b), the lobes begin at normal incidence with approximately equal magnitude and have similar values to those for p-polarization. However, for increasing ${\theta}_{\mathrm{i}}$ the magnitude of the emission in each polar direction tends to reduce. The overall conversion efficiency to synchrotron radiation also reduces in a similar way in Figure 9(e) and remains equal for each polar direction. Rotating the target causes the laser pulse to be spread over a larger area, and to be reflected obliquely. As a result, the radiation pressure that is required for forming a cavity in the target and producing strong laser-injected emission is lowered. Furthermore, because the laser pulse is polarized parallel with the target surface, its electric field cannot directly pull electrons away from the surface and further into its spatial profile without deformation of the surface.

Figure 9 3D simulation results for synchrotron photon emission for different laser light polarization states. Peak angle-resolved synchrotron energy emitted in each direction for (a) p-polarization, (b) s-polarization and (c) left-hand and right-hand circular polarization. (d)–(f) Conversion efficiency to synchrotron radiation for p-, s- and both left-hand and right-hand circular polarization, respectively.

Finally, in Figures 9(c) and 9(f), results for both left-hand circular polarization (LHCP) and right-hand circular polarization (RHCP) are shown. For c-polarization, two lobes are not produced for normal incidence, but an annular structure instead^[ Reference Ji, Pukhov, Nerush, Kostyukov, Shen and Akli ^{27 ]}, due to the rotating electric field pulling electrons into the focal spot from all around the sides of the hole bored plasma cavity by the time it has completed one cycle. The magnitude of the angle-resolved energy of synchrotron emission in the full angular range is shown in Figure 9(c). For normal incidence it is approximately 0.5 J sr⁻¹, which is approximately 1 J sr⁻¹ lower than for linear polarization. However, this more than doubles in rotating the target to ${\theta}_{\mathrm{i}} = 45{}^{\circ}$ . The magnitude of the emission reduces for ${\theta}_{\mathrm{i}} = 60{}^{\circ}$ before peaking again for ${67.5}^{\circ }$ in a similar way to the p-polarization results. The conversion efficiency in Figure 9(f) remains approximately constant for ${\theta}_{\mathrm{i}} = 0^{\circ}$ –45°, before it quickly reduces for larger ${\theta}_{\mathrm{i}}$ . The results here for LHCP and RHCP are almost identical.

The results in Figure 9 indicate that p-polarization and ${\theta}_{\mathrm{i}}\approx 45{}^{\circ}$ –52.5° produce the optimal combination of peak angle-resolved energy of synchrotron emission and conversion efficiency for the target simulated. However, they do not capture all of the changes to the spatial structure of the synchrotron emission. In Figure 10 the angle-resolved synchrotron emission in the forward direction ( $\mid \theta \mid <90{}^{\circ}$ ) is shown for each polarization and a number of different values of ${\theta}_{\mathrm{i}}$ .

Figure 10 Angular profiles of the total energy of synchrotron emission in the forward direction ( $\mid \theta \mid <90{}^{\circ}$ ) in 3D simulations for different laser light polarization states and angles-of-incidence.

For p-polarization a double lobe structure is shown for $\theta = 0^{\circ}$ –45°, in which the ${\theta}_{0,-90}$ lobe gains energy up to 45°, whilst the ${\theta}_{90,0}$ lobe disappears. When ${\theta}_{\mathrm{i}} = 60{}^{\circ}$ the angular structure changes, a beam of synchrotron radiation narrow in $\phi$ ( $\mid \phi -90{}^{\circ}\mid <5{}^{\circ}$ ) is produced for $\theta = -25{}^{\circ}$ –0°, becoming wider like the original lobe structure for lower $\theta$ , and the magnitude is strongly reduced. Lastly, for ${\theta}_{\mathrm{i}} = 67.5{}^{\circ}$ , the full synchrotron emission profile becomes narrower in $\phi$ and strongest for $\theta$ between $-20{}^{\circ}$ and $-10{}^{\circ}$ . This variation with ${\theta}_{\mathrm{i}}$ is not replicated for s-polarization, where only the double lobe structure is present and reduces in magnitude for increasing ${\theta}_{\mathrm{i}}$ , although it curves towards the $-\theta$ direction.

For ${\theta}_{\mathrm{i}} = 67.5{}^{\circ}$ the angle between the laser propagation direction and target surface becomes 22.5°, approaching the divergence half-angle of the laser pulse, ${\theta}_{{\mathrm{div}}} \,{=}\, \sqrt{2\ln \kern0.1em 2}{\lambda}_{\mathrm{L}}/{\pi \phi}_{\mathrm{L}} = 17.2{}^{\circ}$ . When a component of the laser polarization is in the plane of target rotation, as is the case for p- and c-polarization, the laser electric field can inject electrons from close to the target surface further into the laser pulse as it focuses. The electrons are accelerated along the target surface, over a distance much greater than ${\phi}_{\mathrm{L}}$ , towards the focal point where they strongly emit synchrotron radiation in the intense fields. Because the accelerated electron bunches that reach the centre of the laser focus have a narrow divergence, the beam of synchrotron radiation also has a narrow divergence. For the larger ${\theta}_{\mathrm{i}}$ value of 75°, this mechanism is degraded, as shown in Figure 9(a). The expansion of electrons from the target surface affects the laser beam propagation, and the reduced difference in the propagation direction of the laser pulse and electrons moving parallel with the surface will act to reduce ${\chi}_{\mathrm{e}}$ and the synchrotron emission.

For LHCP and RHCP in Figure 10, synchrotron radiation is generated with an annular spatial profile for normal incidence. Rotation of the target, however, produces an asymmetry in the angular distribution of radiation. One section of the annular structure becomes thicker and brighter, the hole of the annulus moves away from $\theta = 0^{\circ}$ and $\phi = 90{}^{\circ}$ and the opposite side of the annulus becomes thinner and dimmer. The spatial profile of the synchrotron radiation now shows its dependence upon the direction of rotation of circular polarization, and the profile is flipped in $\phi$ from LHCP to RHCP. For increasing target rotation, the annular structure disappears and a broad arc of synchrotron radiation is formed with a brighter spot in a single direction. Rotation of the target reduces the number of electrons on one side of the focal spot, thus reducing the number of electrons drawn into the spot by the laser electric field on this side. As the electric field rotates it pulls most electrons from the rotated ( ${\theta}_{\mathrm{i}}>0^{\circ}$ ) target into the focal spot after it points in the $+y$ direction, which results in the production of angular distributions of electrons in the field that are either directed in the $+z$ ( $\phi <90{}^{\circ}$ ) or $-z$ ( $\phi >90{}^{\circ}$ ) directions, depending on the direction of rotation of the electric field. For LHCP, the electric field next rotates towards the $-z$ direction, accelerating electrons in the $+z$ direction and generating synchrotron radiation in the same direction as those electrons. The situation is reversed for RHCP, directing electrons in the field in the $-z$ direction and the concomitant synchrotron radiation.

For ${\theta}_{\mathrm{i}}\ge 60{}^{\circ}$ , the narrow enhanced beam of synchrotron radiation obtained from the acceleration of electrons along the target surface with the edge of the focusing laser beam appears, similar to the changes to the angular synchrotron profile for p-polarization. However, due to the rotation of the electric field vector, the synchrotron radiation is emitted in a small angle above the $xy$ plane for LHCP and below the $xy$ plane for RHCP.

A limited number of 3D simulations were performed for ${I}_{\mathrm{L}} = 3\times {10}^{23}\kern0.1em$ W cm⁻², corresponding to a laser energy of 109 J, to identify if the spatial structures in the synchrotron emission change when the target is strongly deformed by the radiation pressure. The angle-resolved total energy of synchrotron emission for these simulations is shown in the Supplementary Materials. An $l = 1\ \unicode{x3bc}\mathrm{m} $ and 3 μm target was simulated due to the high laser transmission of the former target for normal incidence at this intensity. The values of ${\theta}_{\mathrm{i}}$ were limited to 0°, 45° and 67.5°.

The strongest angle-resolved synchrotron emission of 160 J sr⁻¹ was generated for p-polarization, ${\theta}_{\mathrm{i}} = 45{}^{\circ}$ and $l = 3\kern0.22em \unicode{x3bc}$ m. However, this corresponds to a narrow beam centred at $\theta \approx 0^{\circ}$ and $\phi \approx 90{}^{\circ}$ , instead of an oblique single lobe as in Figure 10. The emission here is driven by electrons from the $y>0$ side of the target injected into the laser pulse and co-propagating with it as it moves further into the target. This is the same mechanism that drives similar emission for the optimum of ${f}_{{\mathrm{O}}2}$ in Figure 7(c), enabled here by the reduced density modelled in 3D. The double lobes for p- and s-polarization are generated closer to the laser propagation direction, and the synchrotron emission for s-polarization is less strongly suppressed for $l = 3\kern0.22em \unicode{x3bc}$ m with increasing ${\theta}_{\mathrm{i}}$ due to the increased target deformation. For c-polarization, the emission is also more collimated with the laser pulse, and the annular structure for normal incidence only appears for $l = 1\kern0.22em \unicode{x3bc}$ m, with the $l = 3\kern0.22em \unicode{x3bc}$ m target producing a spot centred close to the laser propagation direction. For all of these polarization states and $l = 1\kern0.22em \unicode{x3bc}$ m, the peak angle-resolved emission is enhanced for the larger values of ${\theta}_{\mathrm{i}}$ due to the increased target material within the path of the laser pulse.

8 Summary

In summary, Bayesian optimization has been applied to the generation of synchrotron radiation in ultrahigh-intensity laser interactions with CH foils in 2D simulations. Optimization of individual properties is shown, and control of simultaneous optimization of the conflicting objectives of maximizing synchrotron production and minimizing bremsstrahlung emission is demonstrated with changes to the objective function, including the use of an acceptance function.

Through the use of machine learning-based optimization, the angle-of-incidence is identified as a critical parameter in achieving the greatest angle-resolved synchrotron emission. Further 2D and 3D simulations of a $3\times {10}^{22}\kern0.1em$ W cm⁻² intensity, 30 fs pulse duration, $1\kern0.22em \unicode{x3bc}$ m spot size and 10.9 J energy laser pulse support the optimization of the synchrotron emission spatial profile into a single forward-directed lobe. The polarization of the laser pulse is shown to control the angular distribution of synchrotron radiation, due to the synchrotron emission being dominated by electrons pulled into the laser pulse by the laser electric field. Changing the direction of rotation of the electric field vector for circular polarization is shown to flip the spatial profile of the synchrotron emission in the $z$ (polar) direction. Furthermore, rotation of the target to bring the surface close to the divergence half-angle of the tightly focused laser pulse is observed to produce a narrow beam of synchrotron radiation from the electrons accelerated along the surface, provided the laser light is not s-polarized.

The changes induced in the synchrotron spatial profile with laser polarization and angle-of-incidence may enable such radiation to be more easily distinguishable from bremsstrahlung emission, and enhance studies of the QED-plasma physics in these interactions. The demonstrated control of the synchrotron emission is also useful for the application of this intense source of high-energy photons. The work presented here may be further extended by considering the generation of electron–positron pairs from the high-energy photons in the intense laser fields, and the influence of parameters not explored, such as the front surface density scale length, the spatial-intensity contrast^[ Reference Wilson, King, Butler, Carroll, Frazer, Duff, Higginson, Dance, Jarrett, Davidson, Armstrong, Liu, Hawkes, Clarke, Neely, Gray and McKenna ^{89 ]} and different target structures.