Hostname: page-component-7c8c6479df-24hb2 Total loading time: 0 Render date: 2024-03-17T06:38:25.879Z Has data issue: false hasContentIssue false

Reflection of intense laser light from microstructured targets as a potential diagnostic of laser focus and plasma temperature

Published online by Cambridge University Press:  27 December 2018

J. Jarrett
Affiliation:
Department of Physics SUPA, University of Strathclyde, G4 0NG, UK
M. King
Affiliation:
Department of Physics SUPA, University of Strathclyde, G4 0NG, UK
R. J. Gray
Affiliation:
Department of Physics SUPA, University of Strathclyde, G4 0NG, UK
N. Neumann
Affiliation:
Institut für Kernphysik, Technische Universität Darmstadt, 64289 Darmstadt, Germany
L. Döhl
Affiliation:
York Plasma Institute, University of York, YO10 5DD, UK
C. D. Baird
Affiliation:
York Plasma Institute, University of York, YO10 5DD, UK
T. Ebert
Affiliation:
Institut für Kernphysik, Technische Universität Darmstadt, 64289 Darmstadt, Germany
M. Hesse
Affiliation:
Institut für Kernphysik, Technische Universität Darmstadt, 64289 Darmstadt, Germany
A. Tebartz
Affiliation:
Institut für Kernphysik, Technische Universität Darmstadt, 64289 Darmstadt, Germany
D. R. Rusby
Affiliation:
Department of Physics SUPA, University of Strathclyde, G4 0NG, UK STFC-Rutherford Appleton Laboratory, Didcot, OX11 0QX, UK
N. C. Woolsey
Affiliation:
York Plasma Institute, University of York, YO10 5DD, UK
D. Neely
Affiliation:
Department of Physics SUPA, University of Strathclyde, G4 0NG, UK STFC-Rutherford Appleton Laboratory, Didcot, OX11 0QX, UK
M. Roth
Affiliation:
Institut für Kernphysik, Technische Universität Darmstadt, 64289 Darmstadt, Germany
P. McKenna*
Affiliation:
Department of Physics SUPA, University of Strathclyde, G4 0NG, UK
*
Correspondence to:  P. McKenna, Department of Physics SUPA, University of Strathclyde, G4 0NG, UK. Email: paul.mckenna@strath.ac.uk

Abstract

The spatial-intensity profile of light reflected during the interaction of an intense laser pulse with a microstructured target is investigated experimentally and the potential to apply this as a diagnostic of the interaction physics is explored numerically. Diffraction and speckle patterns are measured in the specularly reflected light in the cases of targets with regular groove and needle-like structures, respectively, highlighting the potential to use this as a diagnostic of the evolving plasma surface. It is shown, via ray-tracing and numerical modelling, that for a laser focal spot diameter smaller than the periodicity of the target structure, the reflected light patterns can potentially be used to diagnose the degree of plasma expansion, and by extension the local plasma temperature, at the focus of the intense laser light. The reflected patterns could also be used to diagnose the size of the laser focal spot during a high-intensity interaction when using a regular structure with known spacing.

Type
Research Article
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution, and reproduction in any medium, provided the original work is properly cited.
Copyright
© The Author(s) 2019

1 Introduction

The interaction of high-intensity laser pulses with solid targets results in the production of compact high energy ion sources[Reference Daido, Nishiuchi and Pirozhkov1Reference Higginson, Gray, King, Dance, Williamson, Butler, Wilson, Capdessus, Armstrong, Green, Hawkes, Martin, Wei, Mirfayzi, Yuan, Kar, Borghesi, Clarke, Neely and McKenna3], as well as X-ray[Reference Murnane, Kapteyn and Falcone4], THz[Reference Singh and Sharma5], EMP[Reference Mead, Neely, Gauoin, Heathcote and Patel6, Reference Bradford, Woolsey, Scott, Liao, Liu, Zhang, Zhu, Armstrong, Astbury, Brenner, Brummitt, Consoli, East, Gray, Haddock, Huggard, Jones, Montgomery, Musgrave, Oliveira, Rusby, Spindloe, Summers, Zemaityte, Zhang, Li, McKenna and Neely7] and high harmonic generation[Reference Dromey, Zepf, Gopal, Lancaster, Wei, Krushelnick, Tatarakis, Vakakis, Moustai‘is, Kodama, Tampo, Stoeckl, Clarke, Habara, Neely, Karsch and Norreys8, Reference Ganeev9]. These novel sources are being applied for probing transient laser–plasma interaction physics on picosecond timescales[Reference Borghesi, Campbell, Schiavi, Haines, Willi, MacKinnon, Patel, Gizzi, Galimberti, Clarke, Pegoraro, Ruhl and Bulanov10], micron-scale resolution radiographic density diagnosis[Reference MacKinnon, Patel, Borghesi, Clarke, Freeman, Habara, Hatchett, Hey, Hicks, Kar, Key, King, Lancaster, Neely, Nikkro, Norreys, Notley, Phillips, Romagnani, Snavely, Stephens and Town11] and the isochoric heating of matter[Reference Patel, MacKinnon, Key, Cowan, Foord, Allen, Price, Ruhl, Springer and Stephens12]. In addition, a wide range of potential societal applications, including medical oncology[Reference Bulanov and Khoroshkov13] and fast-ignition inertial confinement fusion[Reference Tabak, Hammer, Glinsky, Kruer, Wilks, Woodworth, Campbell and Perry14, Reference Roth, Cowan, Key, Hatchett, Brown, Fountain, Johnson, Pennington, Snavely, Wilks, Yasuike, Ruhl, Pegoraro, Bulanov, Campbell, Perry and Powell15], have been proposed, all of which motivate investigation of the fundamental physics of intense laser–solid interactions.

The temperature of beams of electrons, ions or photons generated in intense laser–solid interactions is typically diagnosed by spectral measurements using spectrometers[Reference Carroll, Brummitt, Neely, Lindau, Lundh, Wahlström and McKenna16], stacked dosimetry film[Reference Rusby, Wilson, Gray, Dance, Butler, MacLellan, Scott, Bagnoud, Zielbauer, McKenna and Neely17] or nuclear activation[Reference Zepf, Clark, Krushelnick, Beg, Escoda, Dangor, Santala, Tatarakis, Watts, Norreys, Clarke, Davies, Sinclair, Edwards, Goldsack, Spencer and Ledingham18]. Diagnosing the temperature of the dense plasma heated by the intense laser light is more difficult and often involves measuring the ratio of intensities of discrete X-ray emission lines[Reference Zhao, Deng, Xiong, Hu, Wei, Zhu, Shang, Li, Yang, Zhang and Yang19]. This approach involves a spatially averaged measurement over the region of the target defined by the collection angle of the X-ray spectrometer. It remains difficult to diagnose the temperature of the plasma specifically within the spatially localized region of the laser focus, near the critical density ( $n_{c}$ ), and at a time limited to the laser pulse interaction time. We note here that the critical density surface is the density contour at which the plasma becomes opaque to the laser light and at which significant energy coupling to fast electrons occurs, via mechanisms such as resonance absorption[Reference Freidberg, Mitchell, Morse and Rudsinski20], vacuum heating[Reference Brunel21] and J $\times$ B heating[Reference Kruer and Estabrook22]. The critical density is given by $n_{c}=\unicode[STIX]{x1D714}_{L}^{2}\unicode[STIX]{x1D6FE}m_{e}\unicode[STIX]{x1D716}_{0}/e^{2}$ , where $\unicode[STIX]{x1D714}_{L}$ is the laser angular frequency, $m_{e}$ is the rest mass of the electron, $\unicode[STIX]{x1D716}_{0}$ is the permittivity of free space, $-e$ is the electron charge and $\unicode[STIX]{x1D6FE}=1/\sqrt{1-(v/c)^{2}}$ , where $v$ is the quiver velocity of electrons in the laser field.

In addition to difficulties in diagnosing the properties of dense plasma within the spatially localized region of the laser focus at high intensities, there are also challenges in characterizing the on-shot laser focal spot size. This is particularly true in the case of tight focusing geometries, for which the highest laser intensities are achieved, due to divergence of the beam after focus. The focal spot is typically characterized at low powers before and after a high power laser shot. However, thermal-induced lensing or aberrations in the laser amplifying crystals, or self-focusing in preformed plasma at the front of the target can change the focal spot size on a full power laser shot. An ability to confirm the laser focal spot size at high intensity would provide additional confidence in calculated on-shot intensities. This is particularly important when physical processes which scale nonlinearly with intensity occur. For example, ultrathin foils undergoing relativistic self-induced transparency have been found to be highly sensitive to the laser intensity and focal spot distribution[Reference Gonzalez-Izquierdo, Gray, King, Dance, Wilson, McCreadie, Butler, Capdessus, Hawkes, Green, Borghesi, Neely and McKenna23, Reference Gonzalez-Izquierdo, King, Gray, Wilson, Dance, Powell, MacLellan, McCreadie, Butler, Hawkes, Green, Murphy, Stockhausen, Carroll, Booth, Scott, Borghesi, Neely and McKenna24].

In this article, we measure the spatial-intensity profile of light reflected during the interaction of intense laser pulses with microstructured targets (which are being investigated as a route to increase laser energy absorption[Reference Jiang, Krygier, Schumacher, Akli and Freeman25, Reference Margarone, Klimo, Kim, Prokupek, Limpouch, Jeong, Mocek, Psikal, Kim, Proska, Nam, Stolcova, Choi, Lee, Sung, Yu and Korn26]). Patterns in both the reflected laser light at the fundamental harmonic and second harmonic light generated at the critical density surface[Reference Streeter, Foster, Cameron, Borghesi, Brenner, Carroll, Divall, Dover, Dromey, Gallegos, Green, Hawkes, Hooker, Kar, McKenna, Nagel, Najmudin, Palmer, Prasad, Quinn, Rajeev, Robinson, Romagnani, Schreiber, Spindloe, Ter-Avetisyan, Tresca, Zepf and Neely27] are characterized experimentally. We explore the possibility to use these patterns to diagnose the temperature of the plasma, localized spatially to the region of the focal spot and temporally to the peak of the interaction. Through the use of ray-tracing and numerical modelling, we explore how the spatial-intensity profile evolves as the microstructures change due to plasma expansion and relate that to the plasma temperature within the high-intensity region of the focused laser light. The influence of the size of the laser focus relative to the target structures on the reflected light patterns is also investigated, highlighting the potential to use this approach to diagnose the laser focus on high-intensity laser shots.

2 Experiment

The Vulcan laser at the Rutherford Appleton Laboratory, UK, was used in this study. The light was p-polarized at a central wavelength of 1053 nm. The pulses, with duration equal to ( $0.9\pm 0.1$ ) ps (full width at half maximum, FWHM), were focused by an $f/3$ off-axis parabolic mirror and reflected from two consecutive plasma mirrors in order to produce a very high temporal intensity contrast[Reference Dromey, Kar, Zepf and Foster28, Reference Scott, Bagnoud, Brabetz, Clarke, Green, Heathcote, Powell, Zielbauer, Arber, McKenna and Neely29]. An overview of the experimental arrangement is shown in Figures 1(a) and 1(b). This configuration resulted in laser pulses with energy $(160\pm 30)~\text{J}$ incident onto the target at an angle of $20^{\circ }$ to the target normal, focused to a spot diameter of ${\sim}7~\unicode[STIX]{x03BC}\text{m}$ FWHM (as measured at low power). The resulting laser intensity is calculated to be $(5\pm 0.5)\times 10^{19}~\text{W}\cdot \text{cm}^{-2}$ .

Figure 1. (a) Plan view of the experiment arrangement inside the vacuum chamber. The incoming laser beam is shown in red and light reflected out of chamber to the CCDs is shown in blue. (b) Schematic showing the path of the incoming laser beam (solid red line), from the double plasma mirror onto target and finally onto the scatter screen. The imaging line is shown by the dashed red line. (c) Schematic illustrating the four types of targets employed; from left to right: flat foil, grooves, pillars and needles.

Figure 2. Measurements of the spatial-intensity distribution of the laser light reflected from the plasma critical density surface, at fundamental and second harmonic frequencies, as captured on a scatter screen, with dashed red line denoting the expected specular direction. Images (a) and (b) correspond to the $\unicode[STIX]{x1D714}_{L}$ and $2\unicode[STIX]{x1D714}_{L}$ signals for the flat foil target, respectively. (c) and (d) are the same for the groove target, (e) and (f) are for the pillar target, and (g) and (h) are obtained with the needle target. The scale presented in (a) and (b) is the same for all $\unicode[STIX]{x1D714}_{L}$ and $2\unicode[STIX]{x1D714}_{L}$ images, respectively.

Silicon foil targets with four distinct surface profiles, as illustrated in Figure 1(c), were used in the investigation: (1) a flat planar surface; (2) etched parallel linear grooves; (3) a uniform grid of vertical pillars; and, (4) a uniform grid of needles or cone-shaped columns. All targets were comprised of $25~\unicode[STIX]{x03BC}\text{m}$ -thick silicon wafers. The structured targets consisted of a $10~\unicode[STIX]{x03BC}\text{m}$ -thick base, with a $15~\unicode[STIX]{x03BC}\text{m}$ -high structure. The groove structures had a width of $2~\unicode[STIX]{x03BC}\text{m}$ , with $5~\unicode[STIX]{x03BC}\text{m}$ centre-to-centre spacing. The pillar structures had a diameter of $3~\unicode[STIX]{x03BC}\text{m}$ , arranged in a grid layout, again with $5~\unicode[STIX]{x03BC}\text{m}$ centre-to-centre spacing, and the needle structures had a similar grid layout but with a base diameter of $5~\unicode[STIX]{x03BC}\text{m}$ and $5~\unicode[STIX]{x03BC}\text{m}$ centre-to-centre spacing (the base of each needle in contact with adjacent needles). The targets were produced at the Detector and Target Laboratory, at the Technical University of Darmstadt, using repetitive femtosecond laser pulse irradiation in a sulphur hexafluoride environment[Reference Her, Finlay, Wu, Deliwala and Mazur30]. The fabrication method is described in Ref. [Reference Ebert, Neumann, Abel, Schaumann and Roth31].

CCD cameras mounted outside the target chamber were used to image the light reflected onto a ground-glass scatter screen, which had dimensions of $45~\text{cm}\times 45~\text{cm}$ , positioned 50 cm from the target. The reflected light at the fundamental frequency ( $\unicode[STIX]{x1D714}_{L}$ ) and second harmonic light ( $2\unicode[STIX]{x1D714}_{L}$ ) produced at the critical density surface in the most intense region of the focal spot were both measured. This was achieved by using two CCD cameras with separate band-pass filters, with transmission centred at the two chosen wavelengths, respectively.

The patterns measured in the reflected $\unicode[STIX]{x1D714}_{L}$ and generated $2\unicode[STIX]{x1D714}_{L}$ light are shown in Figure 2. A circular disc of light is imaged at both wavelengths for the planar (unstructured) foils. In the case of the periodic groove target, the light forms a structured pattern similar to the diffraction patterns formed by light reflected from a grating structure. The spatial frequency of the structures in the $2\unicode[STIX]{x1D714}_{L}$ light is double that of the $\unicode[STIX]{x1D714}_{L}$ , which is consistent with the patterns being produced by diffraction. In the case of the pillar targets, the periodic pattern is observed in two dimensions, due to the periodicity of the structure along the orthogonal axes. However, in the case of the needle targets a speckle-like pattern is imaged and arises because the light is reflected by the sides of the needles and undergoes multiple reflections.

The line-outs shown in Figure 3, for the $\unicode[STIX]{x1D714}_{L}$ and $2\unicode[STIX]{x1D714}_{L}$ light produced by the groove target, confirm the changes to the periodicity of the pattern expected from diffraction. As the second harmonic light is predominately produced at the critical surface in the intense regions of the plasma, the presence of the periodic pattern indicates that the plasma surface retains some degree of the initial target structure during the interaction with the intense laser light. These reflected light patterns may thus provide valuable insight into the evolution of the critical plasma surface in the localized region of the laser focal spot during the interactions. As such, the following sections detail numerical modelling of the laser–solid interaction to explore the critical surface evolution, along with simplified ray-tracing modelling to investigate the potential of using the specularly reflected light as a diagnostic tool.

Figure 3. Normalized line-outs from the measured (a) $\unicode[STIX]{x1D714}_{L}$ and (b) $2\unicode[STIX]{x1D714}_{L}$ reflected light patterns produced by the groove target. The dashed lines correspond to the expected positions of light maxima from diffraction theory.

3 Modelling

To investigate this potential, first 2D particle-in-cell (PIC) simulations were performed to investigate how the target critical density surface profile changes with time as it is irradiated by a high-intensity laser pulse. The Extendible Particle-in-cell Open Collaboration (EPOCH[Reference Arber, Bennett, Brady, Lawrence-Douglas, Ramsay, Sircombe, Gillies, Evans, Schmitz, Bell and Ridgers32]) PIC code was used to model the interaction between a groove-structured target and a laser pulse defined as a Gaussian beam focused to a spot diameter of $7~\unicode[STIX]{x03BC}\text{m}$ (FWHM) and with a Gaussian temporal profile of 500 fs FWHM. The time $t=0$ fs is defined as the time at which the peak of the pulse reaches the target surface. The peak intensity of the laser pulse was $5\times 10^{19}~\text{W}\cdot \text{cm}^{-2}$ .

The simulations were run with a $40~\unicode[STIX]{x03BC}\text{m}$ $\times 28.8~\unicode[STIX]{x03BC}\text{m}$ simulation box with a mesh resolution of $2~\text{nm}\times 2~\text{nm}$ , for a total simulation time equal to 1.2 ps. The plasma comprised of $\text{Si}^{12+}$ ions, with the electron and ion densities selected to ensure overall neutrality. The initial electron temperature was set at 30 keV and was varied between 30 keV and 80 keV for the case with no laser present. The initial electron density was set to $200n_{c}$ . The groove target structure was initialized with the same periodicity and other dimensions as the targets used in the experiment.

Figure 4. PIC simulation results showing electron density (and thus the groove expansion) for a laser intensity of $5\times 10^{19}~\text{W}\cdot \text{cm}^{-2}$ , pulse duration of 500 fs (FWHM) and $7~\unicode[STIX]{x03BC}\text{m}$ focal spot (FWHM), for (a) $t=-500~\text{fs}$ and (b) $t=0$ . Overdense plasma is shown in red while underdense plasma is shown in blue. The green line in (b) shows the critical density trace as used in the ray-tracing model.

Figure 4(a) shows the initialized cold target electron density from the PIC simulation at $t=-500~\text{fs}$ , and Figure 4(b) shows the corresponding profile at $t=0$ (i.e. the peak of the laser pulse interaction). It is clear that as the laser pulse interacts with the target, the increase in temperature and resulting plasma expansion changes the groove structure. The gaps between the grooves fill with plasma such that the effective groove depth decreases with time.

The evolution of the plasma critical density surface due to thermal expansion was also modelled in the absence of the laser pulse. A total of $10^{8}$ particles were distributed such that they match the geometry of the groove target. The particles were defined with an initial temperature and a normal distribution of random momenta calculated using that temperature, and were then allowed to propagate ballistically with time. The model traces the outline of the critical surface as a function of time for a given starting temperature. Example plots of the critical surface at different time steps are shown in Figure 5. In order to quantify the extent of the plasma expansion at any given time, the depth of the groove structure is determined by calculating the vertical separation between the highest and lowest points (peaks and troughs) of the structure.

Figure 5. Contour plot showing the evolution of the target profile (the plasma critical density surface) as determined from modelling the plasma thermal expansion.

Next, a simplified ray-tracing model was developed to investigate how the profile of the reflected light (at a detector screen) evolves as the plasma expands and the groove depth decreases. This package, based on scripts available from the MathWorks File Exchange[33], is capable of tracing a large number of light rays as they reflect from a user-defined mirror surface. A surface closely matching that of the expanded groove profile in the PIC simulations at the peak of the interaction ( $t=0$ ) was created. The profile is defined as $F(y)=A\sin [(\unicode[STIX]{x1D70B}/S)y+\unicode[STIX]{x1D719}]^{C}$ , where $A$ , $S$ and $C$ are parameters that vary the groove depth, spacing and shape, respectively, and $\unicode[STIX]{x1D719}$ is the phase of the structure. The resultant accumulation of rays reflected from the defined surface produces a series of maxima (at a distance corresponding to the experimental scatter screen distance) that varies depending upon $A$ , $S$ and $C$ .

Figure 6(a) shows a magnified view of the top of the grooves, illustrating how the local curvature, and thus the degree of plasma expansion, change the direction of the reflected light rays. This example is shown for $A=3.7~\unicode[STIX]{x03BC}\text{m}$ , $S=5~\unicode[STIX]{x03BC}\text{m}$ , $C=6$ and $\unicode[STIX]{x1D719}=\unicode[STIX]{x1D70B}$ , such that the separation between the maxima (at the scatter screen distance) is similar to that measured in the experiment. This example follows the surface profile (green curve) in the simulation results in Figure 4(b). Keeping the shape parameter constant, the separation in the maxima varies with groove depth as shown in Figure 6(b).

Figure 6. (a) Magnified view of the top of three groove structures showing reflected light rays, for light incident vertically downwards. (b) Separation of light maxima at the distance of the scatter screen as a function of the groove depth, as determined from the ray-tracing model.

4 Diagnosing plasma temperature and focal spot size

The numerical modelling of the thermal expansion, along with PIC simulations run without the laser pulse present for a variety of different initial temperatures, was used in order to obtain an estimate of the groove depth as a function of plasma temperature after a period of the FWHM of the laser pulse. The results of this model are plotted in Figure 7(a).

The combination of the functional dependence of the reflected light profile on the groove depth and the relationship between plasma temperature and groove depth enables the structure in the spatial-intensity distribution of the reflected light to be correlated to the plasma temperature, as shown by the model calculations in Figure 7(b) along with the corresponding PIC simulation data. The two models are in good agreement and highlight how the profile of the reflected laser light may be developed as a diagnostic of the localized plasma electron temperature.

Figure 7. (a) Groove depth as a function of electron temperature, as determined from the numerical thermal expansion model. (b) Plot of results from numerical modelling, showing expected separation between maxima in reflected light (at the distance of the scatter screen) as a function of plasma electron temperature. The red line represents the numerical model and black dots are data points from the PIC simulations.

In addition to this simple ray-tracing approach, the influence of the wave behaviour of the reflected light was also investigated by utilizing a model based on the Huygens–Fresnel principle. By assuming each point on the surface emits as a spherical wave with the intensity of each point defined by a spatial Gaussian distribution (to model an effective focal spot), it is possible to numerically integrate the resultant intensity distribution on a plane at a distance from the surface. This was set at 1 mm to ensure the plane is at a sufficient distance from the surface and yet limit the computational requirements. Figure 8(a) illustrates the behaviour of the intensity distribution for a focal spot size of $2~\unicode[STIX]{x03BC}\text{m}$ as the groove depth, $A$ , is varied in the same manner as the ray-tracing model. The spacing between the maxima is observed to decrease as the groove depth is increased, tending to the separation expected from simple diffraction theory. This behaviour is similar to that of the ray-tracing model, but also shows the expected destructive interference produced as $A$ varies over one wavelength, resulting in the diffraction pattern appearing to oscillate with $A$ . Experimentally, this behaviour may lead to ambiguous results, however, as $A$ evolves over the course of the interaction with the laser pulse. When the focal spot diameter is increased to $7~\unicode[STIX]{x03BC}\text{m}$ , there is no change in behaviour with $A$ , as shown in Figure 8(b). This indicates that in order to observe the changes in the diffraction pattern seen in the ray-tracing model, the focal spot size must be reduced to the order of the groove spacing or below.

Although this highlights a limitation in using this approach to diagnose plasma expansion and thus temperature, the dependence of the pattern on laser focal spot size, or more correctly the ratio of the focal spot diameter to the groove separation, suggests that this approach may be used to diagnose the size of the focal spot on a full power laser shot.

To explore this aspect further, Figure 9(a) shows the effect of increasing the groove spacing, $S$ , on the reflected intensity profile, for a fixed focal spot FWHM of $7~\unicode[STIX]{x03BC}\text{m}$ and $A=4~\unicode[STIX]{x03BC}\text{m}$ . With increasing $S$ the diffraction pattern becomes less pronounced, until it disappears for values significantly above the focal spot size. Thus the profile of the reflected light depends on the ratio of the laser focus to the groove spacing and it follows that by engineering the target to have specific values of $S$ the size of the laser focus can be determined from the measured reflected light profile.

In all of the Huygens–Fresnel analysis above, it is assumed that the pulse is centred on the middle of the groove spacing. The precise position of the focus with respect to the groove profile in an experiment is subject to the laser pointing stability. To investigate what effect this has on this diagnostic approach, we next vary the relative position of the laser focus and the target structure in our model. For convenience in the application of our model, we simply vary the phase term in the structure, keeping the centre of the laser focus fixed at $y=0$ . Figure 9(b) shows the change in reflected intensity distribution as the defined structure is moved in phase, $\unicode[STIX]{x1D719}$ , through $\unicode[STIX]{x1D70B}$ (as the structure is defined as a $\sin ^{2}$ function) for a focal spot of $7~\unicode[STIX]{x03BC}\text{m}$ and $S=35~\unicode[STIX]{x03BC}\text{m}$ . As the spot size is much smaller than $S$ , when the centre of the focal spot moves away from the middle of the groove structure, the reflection from the structure slope dominates, resulting in a displacement of the central maxima in the reflected pattern. Additional maxima are observed, skewed in one direction. As the relative displacement is increased, the slope of the structure changes direction resulting in the central maxima displacement inverting. This indicates that although there is a variation depending upon where on the structure the focal spot is centred, the variation is predictable. When $S$ is reduced to $20~\unicode[STIX]{x03BC}\text{m}$ , as shown in Figure 9(c), the central maximum is displaced to a similar extent, but the diffraction effects become stronger. Finally, as shown in Figure 9(d), reducing $S$ to $5~\unicode[STIX]{x03BC}\text{m}$ , a value smaller than the focal spot size, the relative position of the focal spot has no effect on the diffraction pattern, as expected.

Figure 8. Intensity distribution determined from a Huygens–Fresnel model at a plane 1 mm from an evolved groove structure as a function of $A$ , with $S=5~\unicode[STIX]{x03BC}\text{m}$ and laser focal spot FWHM equal to (a) $2~\unicode[STIX]{x03BC}\text{m}$ and (b) $7~\unicode[STIX]{x03BC}\text{m}$ . The solid lines indicate the expected first order diffraction position and the dashed lines correspond to the results determined from the model.

Figure 9. Intensity distribution determined from a Huygens–Fresnel model at a plane 1 mm from an evolved groove structure as a function of $S$ , with $A=4~\unicode[STIX]{x03BC}\text{m}$ and laser focal spot FWHM equal to $7~\unicode[STIX]{x03BC}\text{m}$ . (b)–(d) The intensity profile as the phase of the structure is varied for (b) $S=35~\unicode[STIX]{x03BC}\text{m}$ , (c) $S=20~\unicode[STIX]{x03BC}\text{m}$ and (d) $S=5~\unicode[STIX]{x03BC}\text{m}$ .

5 Summary

This study shows that by analysing the pattern of the laser light reflected from a microstructured target it may be possible to diagnose the local plasma temperature in the intense region of the laser focus. This approach can also potentially be applied to determine the size of the laser focal spot at full power. It has been verified experimentally that the target structure remains intact long enough during the interaction with the intense laser pulse for periodic structures to be produced in the reflected light at both $\unicode[STIX]{x1D714}_{L}$ and $2\unicode[STIX]{x1D714}_{L}$ . The plasma surface evolves over the interaction, but the temporal window over which the expansion is probed can be reduced by using second (or higher order) harmonic light, which is only produced above a threshold intensity. The approach is limited, however, to laser intensities below the threshold for which the laser radiation pressure exceeds the thermal plasma pressure. Above this, plasma expansion is replaced by laser hole-boring into the target, resulting in additional deformation of the critical density surface[Reference Gonzalez-Izquierdo, Capdessus, King, Gray, Wilson, Dance, McCreadie, Butler, Hawkes, Green, Booth, Borghesi, Neely and McKenna34]. By employing multiple microstructured targets with different groove spacings, it may also be possible to determine the actual focal spot size due to the presence or absence of the expected diffraction structures in the specularly reflected light. The approach could also potentially be extended by engineering variations into the groove pattern to enable the degree of asymmetry in the spatial-intensity distribution of the laser focus to be determined.

Acknowledgements

The authors acknowledge the expertise of the STFC Central Laser Facility staff and the Detector and Target Laboratory staff at the Institut fur Kernphysik, Technische Universitat Darmstadt. Simulations were performed using the EPOCH PIC code (developed under EPSRC grant number EP/G054940/1), running on the ARCHER high performance computer, with access provided via the EPSRC-funded Plasma Physics HEC Consortia (EP/L000237/1). This work is financially supported by EPSRC (grant numbers EP/R006202/1 and EP/K022415/1) and the European Union Horizon 2020 research and innovation programme under grant agreement number 654148 Laserlab-Europe. The data presented in this paper is available at http://dx.doi.org/10.15129/6805696a-13c1-4c56-a2db-d5eb5a471d33.

References

Daido, H. Nishiuchi, M. and Pirozhkov, A. Rep. Prog. Phys. 75, 056401 (2012).Google Scholar
Macchi, A. Borghesi, M. and Passoni, M. Rev. Mod. Phys. 85, 751 (2013).Google Scholar
Higginson, A. Gray, R. J. King, M. Dance, R. J. Williamson, S. D. R. Butler, N. M. H. Wilson, R. Capdessus, R. Armstrong, C. Green, J. S. Hawkes, S. J. Martin, P. Wei, W. Q. Mirfayzi, S. R. Yuan, X. H. Kar, S. Borghesi, M. Clarke, R. J. Neely, D. and McKenna, P. Nature Commun. 9, 724 (2018).Google Scholar
Murnane, M. M. Kapteyn, H. C. and Falcone, R. W. Phys. Rev. Lett. 62, 155 (1989).Google Scholar
Singh, M. and Sharma, R. P. Contrib. Plasma Phys. 53, 540 (2013).Google Scholar
Mead, M. J. Neely, D. Gauoin, J. Heathcote, R. and Patel, P. Rev. Sci. Instrum. 75, 4225 (2004).Google Scholar
Bradford, P. Woolsey, N. C. Scott, G. G. Liao, G. Liu, H. Zhang, Y. Zhu, B. Armstrong, C. Astbury, S. Brenner, C. Brummitt, P. Consoli, F. East, I. Gray, R. Haddock, D. Huggard, P. Jones, P. J. R. Montgomery, E. Musgrave, I. Oliveira, P. Rusby, D. R. Spindloe, C. Summers, B. Zemaityte, E. Zhang, Z. Li, Y. McKenna, P. and Neely, D. High Power Laser Sci. Eng. 6, e21 (2018).Google Scholar
Dromey, B. Zepf, M. Gopal, A. Lancaster, K. Wei, M. S. Krushelnick, K. Tatarakis, M. Vakakis, N. Moustai‘is, S. Kodama, R. Tampo, M. Stoeckl, C. Clarke, R. Habara, H. Neely, D. Karsch, S. and Norreys, P. Nature Phys. 2, 456 (2006).Google Scholar
Ganeev, R. A. J. Phys. B 40, R213 (2007).Google Scholar
Borghesi, M. Campbell, D. H. Schiavi, A. Haines, M. G. Willi, O. MacKinnon, A. J. Patel, P. Gizzi, L. A. Galimberti, M. Clarke, R. J. Pegoraro, F. Ruhl, H. and Bulanov, S. Phys. Plasmas 9, 2214 (2002).Google Scholar
MacKinnon, A. J. Patel, P. K. Borghesi, M. Clarke, R. C. Freeman, R. R. Habara, H. Hatchett, S. P. Hey, D. Hicks, D. G. Kar, S. Key, M. H. King, J. A. Lancaster, K. Neely, D. Nikkro, A. Norreys, P. A. Notley, M. M. Phillips, T. W. Romagnani, L. Snavely, R. A. Stephens, R. B. and Town, R. P. J. Phys. Rev. Lett. 97, 045001 (2006).Google Scholar
Patel, P. K. MacKinnon, A. J. Key, M. H. Cowan, T. E. Foord, M. E. Allen, M. Price, D. F. Ruhl, H. Springer, P. T. and Stephens, R. Phys. Rev. Lett. 91, 125004 (2003).Google Scholar
Bulanov, S. and Khoroshkov, V. S. Plasma Phys. Rep. 28, 453 (2002).Google Scholar
Tabak, M. Hammer, J. Glinsky, M. E. Kruer, W. L. Wilks, S. C. Woodworth, J. Campbell, E. M. and Perry, M. D. Phys. Plasmas 1, 1626 (1994).Google Scholar
Roth, M. Cowan, T. E. Key, M. H. Hatchett, S. P. Brown, C. Fountain, W. Johnson, J. Pennington, D. M. Snavely, R. A. Wilks, S. C. Yasuike, K. Ruhl, H. Pegoraro, F. Bulanov, S. V. Campbell, E. M. Perry, M. D. and Powell, H. Phys. Rev. Lett. 86, 436 (2001).Google Scholar
Carroll, D. C. Brummitt, P. Neely, D. Lindau, F. Lundh, O. Wahlström, C. G. and McKenna, P. Nucl. Instrum. Methods Phys. A 620, 23 (2010).Google Scholar
Rusby, D. Wilson, L. A. Gray, R. J. Dance, R. J. Butler, N. M. H. MacLellan, D. A. Scott, G. G. Bagnoud, V. Zielbauer, B. McKenna, P. and Neely, D. J. Plasma Phys 81, 475810505 (2015).Google Scholar
Zepf, M. Clark, E. L. Krushelnick, K. Beg, F. N. Escoda, C. Dangor, A. E. Santala, M. I. K. Tatarakis, M. Watts, I. F. Norreys, P. A. Clarke, R. J. Davies, J. R. Sinclair, M. A. Edwards, R. D. Goldsack, T. J. Spencer, I. and Ledingham, K. W. D. Appl. Phys. Lett. 8, 2323 (2001).Google Scholar
Zhao, Y. Deng, B. Xiong, G. Hu, Z. M. Wei, M. X. Zhu, T. Shang, W. L. Li, J. Yang, G. H. Zhang, J. Y. and Yang, J. M. Phys. Lett. 28, 065201 (2011).Google Scholar
Freidberg, J. P. Mitchell, R. W. Morse, R. L. and Rudsinski, L. I. Phys. Rev. Lett. 28, 795 (1972).Google Scholar
Brunel, F. Phys. Rev. Lett. 59, 52 (1987).Google Scholar
Kruer, L. and Estabrook, K. Phys. Fluids. 28, 430 (1985).Google Scholar
Gonzalez-Izquierdo, B. Gray, R. J. King, M. Dance, R. J. Wilson, R. McCreadie, J. Butler, N. M. H. Capdessus, R. Hawkes, S. Green, J. S. Borghesi, M. Neely, D. and McKenna, P. Nature Phys. 12, 505 (2016).Google Scholar
Gonzalez-Izquierdo, B. King, M. Gray, R. J. Wilson, R. Dance, R. J. Powell, H. MacLellan, D. A. McCreadie, J. Butler, N. M. H. Hawkes, S. Green, J. S. Murphy, C. D. Stockhausen, L. C. Carroll, D. C. Booth, N. Scott, G. G. Borghesi, M. Neely, D. and McKenna, P. Nature Commun. 7, 12891 (2016).Google Scholar
Jiang, S. Krygier, A. G. Schumacher, D. W. Akli, K. U. and Freeman, R. R. Phys. Rev. E. 89, 013106 (2014).Google Scholar
Margarone, D. Klimo, O. Kim, I. J. Prokupek, J. Limpouch, J. Jeong, T. M. Mocek, T. Psikal, J. Kim, H. T. Proska, J. Nam, K. H. Stolcova, L. Choi, I. W. Lee, S. K. Sung, J. H. Yu, T. J. and Korn, G. Phys. Rev. Lett. 109, 234801 (2012).Google Scholar
Streeter, M. J. V. Foster, P. S. Cameron, F. H. Borghesi, M. Brenner, C. Carroll, D. C. Divall, E. Dover, N. P. Dromey, B. Gallegos, P. Green, J. S. Hawkes, S. Hooker, C. J. Kar, S. McKenna, P. Nagel, S. R. Najmudin, Z. Palmer, C. A. J. Prasad, R. Quinn, K. E. Rajeev, P. P. Robinson, A. P. L. Romagnani, L. Schreiber, J. Spindloe, C. Ter-Avetisyan, S. Tresca, O. Zepf, M. and Neely, D. New J. Physics. 13, 023041 (2011).Google Scholar
Dromey, B. Kar, S. Zepf, M. and Foster, P. Rev. Sci. Instrum. 75, 645 (2004).Google Scholar
Scott, G. G. Bagnoud, V. Brabetz, C. Clarke, R. J. Green, J. S. Heathcote, R. I. Powell, H. W. Zielbauer, B. Arber, T. D. McKenna, P. and Neely, D New J. Phys. 17, 033027 (2015).Google Scholar
Her, T.-H. Finlay, R. J. Wu, C. Deliwala, S. and Mazur, E. Appl. Phys. Lett. 73, 1673 (1998).Google Scholar
Ebert, T. Neumann, N. W. Abel, T. Schaumann, G. and Roth, M. High Power Laser Sci. Eng. 13, e18 (2017).Google Scholar
Arber, T. D. Bennett, K. Brady, C. S. Lawrence-Douglas, A. Ramsay, M. G. Sircombe, N. J. Gillies, P. Evans, R. G. Schmitz, H. Bell, A. R. and Ridgers, C. P. Plasma Phys. Control. Fusion 57, 113001 (2015).Google Scholar
Gonzalez-Izquierdo, B. Capdessus, R. King, M. Gray, R. J. Wilson, R. Dance, R. J. McCreadie, J. Butler, N. M. H. Hawkes, St. J. Green, J. S. Booth, N. Borghesi, M. Neely, D. and McKenna, P. Appl. Sci. 8, 336 (2018).Google Scholar
Figure 0

Figure 1. (a) Plan view of the experiment arrangement inside the vacuum chamber. The incoming laser beam is shown in red and light reflected out of chamber to the CCDs is shown in blue. (b) Schematic showing the path of the incoming laser beam (solid red line), from the double plasma mirror onto target and finally onto the scatter screen. The imaging line is shown by the dashed red line. (c) Schematic illustrating the four types of targets employed; from left to right: flat foil, grooves, pillars and needles.

Figure 1

Figure 2. Measurements of the spatial-intensity distribution of the laser light reflected from the plasma critical density surface, at fundamental and second harmonic frequencies, as captured on a scatter screen, with dashed red line denoting the expected specular direction. Images (a) and (b) correspond to the $\unicode[STIX]{x1D714}_{L}$ and $2\unicode[STIX]{x1D714}_{L}$ signals for the flat foil target, respectively. (c) and (d) are the same for the groove target, (e) and (f) are for the pillar target, and (g) and (h) are obtained with the needle target. The scale presented in (a) and (b) is the same for all $\unicode[STIX]{x1D714}_{L}$ and $2\unicode[STIX]{x1D714}_{L}$ images, respectively.

Figure 2

Figure 3. Normalized line-outs from the measured (a) $\unicode[STIX]{x1D714}_{L}$ and (b) $2\unicode[STIX]{x1D714}_{L}$ reflected light patterns produced by the groove target. The dashed lines correspond to the expected positions of light maxima from diffraction theory.

Figure 3

Figure 4. PIC simulation results showing electron density (and thus the groove expansion) for a laser intensity of $5\times 10^{19}~\text{W}\cdot \text{cm}^{-2}$, pulse duration of 500 fs (FWHM) and $7~\unicode[STIX]{x03BC}\text{m}$ focal spot (FWHM), for (a) $t=-500~\text{fs}$ and (b) $t=0$. Overdense plasma is shown in red while underdense plasma is shown in blue. The green line in (b) shows the critical density trace as used in the ray-tracing model.

Figure 4

Figure 5. Contour plot showing the evolution of the target profile (the plasma critical density surface) as determined from modelling the plasma thermal expansion.

Figure 5

Figure 6. (a) Magnified view of the top of three groove structures showing reflected light rays, for light incident vertically downwards. (b) Separation of light maxima at the distance of the scatter screen as a function of the groove depth, as determined from the ray-tracing model.

Figure 6

Figure 7. (a) Groove depth as a function of electron temperature, as determined from the numerical thermal expansion model. (b) Plot of results from numerical modelling, showing expected separation between maxima in reflected light (at the distance of the scatter screen) as a function of plasma electron temperature. The red line represents the numerical model and black dots are data points from the PIC simulations.

Figure 7

Figure 8. Intensity distribution determined from a Huygens–Fresnel model at a plane 1 mm from an evolved groove structure as a function of $A$, with $S=5~\unicode[STIX]{x03BC}\text{m}$ and laser focal spot FWHM equal to (a) $2~\unicode[STIX]{x03BC}\text{m}$ and (b) $7~\unicode[STIX]{x03BC}\text{m}$. The solid lines indicate the expected first order diffraction position and the dashed lines correspond to the results determined from the model.

Figure 8

Figure 9. Intensity distribution determined from a Huygens–Fresnel model at a plane 1 mm from an evolved groove structure as a function of $S$, with $A=4~\unicode[STIX]{x03BC}\text{m}$ and laser focal spot FWHM equal to $7~\unicode[STIX]{x03BC}\text{m}$. (b)–(d) The intensity profile as the phase of the structure is varied for (b) $S=35~\unicode[STIX]{x03BC}\text{m}$, (c) $S=20~\unicode[STIX]{x03BC}\text{m}$ and (d) $S=5~\unicode[STIX]{x03BC}\text{m}$.