## Towards an in situ , full-power gauge of the focal-volume intensity of petawatt-class lasers

C. Z. HE, 1,2 A. LONGMAN, 3 J. A. PÉREZ-HERNÁNDEZ, 4 M. DE MARCO, 4 , C. SALGADO, 4 G. ZERAOULI, 4 G. GATTI, 4 L. ROSO, 4 R. FEDOSEJEVS, 3 AND W. T. HILL, III 1,2,5,*

1 Joint Quantum Institute, University of Maryland, College Park, MD 20742, USA

- 2 Institute for Physical Science and Technology, University of Maryland, College Park, MD 20742, USA
- 3 Electrical and Computer Engineering, University of Alberta, Edmonton, Alberta T6G 2V4, Canada
- 4 Centro de Láseres Pulsados (CLPU), 37185 Villamayor, Salamanca, Spain
- 5 Department of Physics, University of Maryland, College Park, MD 20742, USA
- * wth@umd.edu

Abstract: About 50 years ago, Sarachick and Schappert [Phys. Rev. D. 1 , 2738-2752 (1970)] showed that relativistic Thomson scattering leads to wavelength shifts that are proportional to the laser intensity. About 28 years later Chen et al. [Nature 396 , 653-655 (1998)] used these shifts to estimate their laser intensity near 10 18 W/cm 2 . More recently there have been several theoretical studies aimed at exploiting nonlinear Thomson scattering as a tool for direct measurement of intensities well into the relativistic regime. We present the first quantitative study of this approach for intensities between 10 18 and 10 19 W/cm 2 . We show that the spectral shifts are in reasonable agreement with estimates of the peak intensity extracted from images of the focal area obtained at reduced power. Finally, we discuss the viability of the approach, its range of usefulness and how it might be extended to gauge intensities well in excess of 10 19 W/cm 2 .

© 2019 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

## 1. Introduction

Since the invention of the laser, physicists have been dreaming of ways to exploit super-intense laser fields to explore new physics in the relativistic regime. Arguably, the most fundamental idea today centers on the nature of the quantum vacuum, which could hold key information necessary for unraveling the mystery of dark matter [1]. Even with current advances, we are still orders of magnitude away from reaching the threshold intensity for creating electron-positron pairs directly via photon-photon collisions, the non-perturbative production of which is estimated to be GLYPH&lt;24&gt; 2 GLYPH&lt;2&gt; 10 29 W/cm 2 . Nevertheless, petawatt-class lasers of short duration have placed us on the cusp of being able to examine experimentally, nonlinear aspects of electrodynamics [1-4], photon-photon interaction and precursor processes to pair production [5-10], which comprise tests of quantum electrodynamics (QED) in ways it has never been tested (see, for example [4,11] and reference therein). Precision measurements of QED require accurate knowledge of the intensity. Measuring the intensity in the focus at relativistic strengths is challenging. We define the threshold for relativistic strength as I GLYPH&lt;24&gt; 10 18 W/cm 2 , the point at which the motion of free electrons becomes relativistic when exposed to a laser field of this strength. A measurement technique must not only be able to determine the true intensity correctly, it also should be minimally intrusive and cause no optical damage to instrumentation. In addition, the approach would ideally be single-shot capable, allowing the intensity to be determined with each laser pulse. Developing suitable techniques is a long-standing desire of the laser community at large. A host of other studies from generating secondary sources [12] to high-energy density physics [13-15] to medical applications [16-21] also would benefit from better knowledge of the intensity.

A number of techniques and proposals for intensity assessment have emerged over the years. Most, however, rely on approximating the peak intensity with methods requiring the laser to be run at significantly reduced powers. A one-to-one correspondence between what is extracted in these cases and the actual intensity, has never been established experimentally. Since the early observation of Ref. [22], there have been three proposals for direct-measurement that can be found in the literature, which, in principle, meet the criteria above. Two are based on relativistic Thomson scattering (RTS) [23-29] and the third on the appearance intensity for inner-shell tunnel ionization [30]. To our knowledge, there have been no reports of experimental verification of any of these approaches for I &gt; 10 18 W/cm 2 . This communication is devoted to addressing that void. Specifically, we present the first peak intensity assessment using the RTS spectrum significantly above 10 18 W/cm 2 . We compare the RTS peak intensity ( I RTS pk ) with estimates extracted from low-power images of the spatial profile of the focal spot ( I Im pk ). Our analysis shows that the spectral and temporal phase must be monitored in the target area to ensure that the Strehl ratio and pulse widths, typically measured in the laser bay, are close to their optimal values after beam transport and pulse compression to maximize the energy in the central Airy spot.

We begin our discussion by briefly outlining the principle of the RTS intensity gauge, which exploits free electrons that are always produced via ionization when I &amp; 10 14 W/cm 2 . Electrons undergoing RTS produce a rich spectrum of intensity-dependent, Doppler-shifted laser wavelengths and their harmonics [24]. In a frame moving in the direction of the laser fi k vector (the average drift frame), the electron approximately executes a figure-8 (see Fig. 1 in [23]) while in the lab frame the electron oscillates in the fi fi Ek -plane while being propelled along fi k , punctuated by cusps (see Fig. 2 in [28]) where large acceleration exists. The benefit of exploiting RTS radiation is threefold:

- · it can be generated from low-density, background gas in the experimental chamber;
- · it can be captured with off-the-shelf optics near, but outside, the focal volume; and
- · the Doppler shifts are straightforward to calculate, based on earlier theoretical work [23-26], leading to a simple expression for the intensity as a function of the observed n th -harmonic wavelength, GLYPH&lt;21&gt; ' n ' ,

$$\L, \, \lambda ^ { \dots \nu }, \, & I ( \lambda ^ { ( n ) }, \theta ; n ) = \frac { 2 \pi m _ { e } c ^ { 3 } } { r _ { 0 } \lambda _ { 0 } ^ { 2 } ( 1 - \cos \theta ) } \left ( \frac { n \lambda ^ { ( n ) } } { \lambda _ { 0 } } - 1 \right ), & & ( 1 )$$

where GLYPH&lt;21&gt; 0 , GLYPH&lt;18&gt; , m e , r 0 and c are the laser wavelength, the polar angle relative to fi k , the electron mass, the classical electron radius and the vacuum light speed respectively.

At relativistic intensities, when the normalized vector potential is of order unity or greater, GLYPH&lt;21&gt; ' n ' differs measurably from GLYPH&lt;21&gt; 0 GLYPH&lt;157&gt; n . The normalized vector potential is defined as

$$a _ { 0 } = e E \lambda _ { 0 } / 2 \pi m _ { e } c ^ { 2 } \approx 0. 8 5 5 \sqrt { I [ \times 1 0 ^ { 1 8 } \, \mathbf W / \text{cm} ^ { 2 } ] } \lambda _ { 0 } [ \mu m ],$$

where e is the elementary charge. In addition, for intensities in excess of 10 25 W/cm 2 such that a 0 ' m e GLYPH&lt;157&gt; m p ' GLYPH&lt;24&gt; 1, where m p is the proton mass, the technique might be extendable to an effective proton RTS. These features make RTS a potentially useful diagnostic tool for a variety of experiments at relativistic intensities.

## 2. Experimental method

Our measurements were done at the Centro de Láseres Pulsados (CLPU) [31] with a 200 TW Ti:sapphire laser (VEGA 2) running at 10 Hz with GLYPH&lt;21&gt; 0 GLYPH&lt;24&gt; 790 nm and compressed to nominally 30 fs. For our experiments, linearly polarized pulses were used with uncompressed energies between 2 and 6 J. About 47% of this energy, the energy on target ( U ), was delivered to the focal volume as compressed pulses after being reflected by an f GLYPH&lt;157&gt; 13, focal length (f.l.) = 130 cm,

Fig. 1. (a) Cross sectional image of the focal volume and an expanded view of the central Airy disk and portions of three secondary rings. The solid (dashed) white circle indicates the extent of the pulse (central Airy disk). The signal in each image was smoothed and filtered to remove hot spots and reduce the noise, and clipped at the background count GLYPH&lt;24&gt; 12. (b) Schematic of the focal region showing the orientation of both microscope objectives ( E ? -obj and E k -obj), gas nozzle, fi E and fi k vectors. Note, the ˆ and x ˆ y axes correspond to ( GLYPH&lt;18&gt;; GLYPH&lt;30&gt; ) = ( GLYPH&lt;25&gt; GLYPH&lt;157&gt; 2 0) and ( ; GLYPH&lt;25&gt; GLYPH&lt;157&gt; 2 ; GLYPH&lt;25&gt; GLYPH&lt;157&gt; 2) respectively, with GLYPH&lt;18&gt; and GLYPH&lt;30&gt; being the polar and azimuthal angles respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00110v3-figure-a5f1a1d8-f941-42f8-a589-f83840cf3a91.jpg)

off-axis parabola (OAP). The diameter (FWHM, i.e., 2 w 0 ) of the central Airy peak, measured at low power, was GLYPH&lt;24&gt; 22 GLYPH&lt;22&gt; m( GLYPH&lt;24&gt; 2 GLYPH&lt;2&gt; the diffraction limit) and contained GLYPH&lt;24&gt; 33% of U , i.e., GLYPH&lt;24&gt; 16% of the uncompressed energy. Figure 1 (a) shows an image of the focal spot taken with a microscope objective at low intensity.

The RTS radiation was captured from two orthogonal directions: ( GLYPH&lt;18&gt;; GLYPH&lt;30&gt; ) = ( GLYPH&lt;25&gt; GLYPH&lt;157&gt; 2 0) and ( ; GLYPH&lt;25&gt; GLYPH&lt;157&gt; 2 ; GLYPH&lt;25&gt; GLYPH&lt;157&gt; 2) respectively, see Fig. 1 (b). Infinity-corrected microscope objectives (Mitutoyo MLWD-10 GLYPH&lt;2&gt; and Motic ELWD-10 GLYPH&lt;2&gt; ) were used, each having an effective f : l : = 20 mm and numerical aperture 0 28, corresponding to a half angle for light collection of : GLYPH&lt;24&gt; 16 . The objectives were adjusted GLYPH&lt;14&gt; with mirror mounts placed on stepper-motor controlled xyz stages. Preliminary alignment of the objectives was enabled by scatting HeNe and diode beams, collinear with the VEGA-2 beam, off of a needle placed at the intended OAP focus. Next, a plasma column was generated by running VEGA 2 at low energy, with the focal volume filled with a few tens of millibars of N 2 . The emission column was then focused onto the detectors by adjusting both the OAP and the objectives position. Final adjustments were done by optimizing the RTS signal at full power. Light captured by the objectives was directed along light-shielded paths (inside and outside the chamber) through a viewport to be analyzed. A second component of the scattered-light

Electrons were generated from N 2 (which yielded up to ten electrons per molecule), delivered to the focal volume by a custom-machined conical nozzle (0.8 mm throat orifice diameter) placed GLYPH&lt;24&gt; 1 8 cm below the focal point, a distance chosen to minimize scattered light from the nozzle : while maximizing gas delivered to the interaction region. The nozzle backing pressure was 10 GLYPH&lt;6&gt; 3 mbar. The orifice was open for 4 ms and triggered GLYPH&lt;24&gt; 2 5 ms before each pulse. : The pressure in the focal volume was estimated to be GLYPH&lt;24&gt; 10 GLYPH&lt;0&gt; 3 mbar, yielding a density of GLYPH&lt;24&gt; 2 4 : GLYPH&lt;2&gt; 10 13 cm GLYPH&lt;0&gt; 3 , assuming ideal gas behavior. The focal volume was taken to be a cylindrical column with a diameter equal to the FWHM diameter of the central Airy peak ( GLYPH&lt;24&gt; 22 GLYPH&lt;22&gt; m). The effective length of the column was defined by the field of view of the collection optics, GLYPH&lt;24&gt; 85 GLYPH&lt;22&gt; m(see text below and Fig. 2), leading to a focal volume of GLYPH&lt;24&gt; 3 GLYPH&lt;2&gt; 10 GLYPH&lt;0&gt; 8 cm 3 . When the intensity reaches GLYPH&lt;24&gt; 10 16 W/cm 3 , a combination of mechanisms including enhanced ionization [32], Coulomb explosion (once three or more electrons are removed) and barrier suppression/tunnel (so-called ADK) ionization [33, 34] induce N 2 to lose all six valence ( p ) electrons leading to N 3 + ions. Additional (inner-shell) electrons can be removed via post-dissociation ADK ionization as the intensity rises through 10 17 and 10 19 W/cm 2 .

Fig. 2. Composite images (500 shots each) obtained with E ? -obj (a) and E k -obj (b), for U = 2 67 J, which includes both RTS and recombination radiation, where : L , h and b are the FWHM length k to fi k , multi-shot width ? to fi k and single-shot width ? to fi k (see text) respectively with values L = 526 85 GLYPH&lt;157&gt; GLYPH&lt;22&gt; m, h = 66 41 GLYPH&lt;157&gt; GLYPH&lt;22&gt; m and b = 15 GLYPH&lt;22&gt; m N.A. for GLYPH&lt;157&gt; image (a)/image (b) respectively. Composite images of the dispersed light and associated integrated profiles (c) for the prompt signal (A), obtained within a 5 ns temporal window that starts GLYPH&lt;24&gt; 2 5 ns before the laser pulse, and the delayed signal (B), obtained within a 50 ns : temporal window starting GLYPH&lt;24&gt; 2 5 ns after the laser pulse. :

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00110v3-figure-8c967436-bfd0-496d-b9dc-2b14b613da03.jpg)

remediation was a beam dump housed in a separate chamber GLYPH&lt;24&gt; 1 m after the focal volume, consisting of flat-black aluminum foil oriented at 45 GLYPH&lt;14&gt; with respect to fi k and enclosed by two coaxial cylindrical walls (flat-black foil) with an entrance hole just large enough for the laser to pass unobstructed.

Data was collected in 10 s bursts (100 shots/burst); the bursts were repeated multiple times (typically five) to create a shot batch. Each shot batch, obtained with the nozzle open, was followed by a second series with the nozzle closed - a background batch. Composite images were generated by combining the two batches: (shot batch) - (background batch). The experimental data set consisted of composite images obtained under different conditions. Specifically, we varied the pulse energy and chirp, the optical filters (i.e., no filter, with BG39 filter or with RG850 filter), the temporal collection window for the Andor iStar CCD camera and the target pressure. Examples of composite images of the focal region, 500 shots without spectral dispersion, in the 2 nd harmonic ( n = 2) spectral region, are shown in Fig. 2 (a) and (b). Examples of dispersed composite images for the same spectral region are shown in Fig. 2 (c).

We detected E k -obj light either directly with a cooled, low-noise CMOS camera (Zyla 4.2MPPLUS) or by a gateable (5 ns resolution) intensified CCD camera (Andor iStar DH 720-25U-03 ICCD) after being dispersed by a 0 25-m spectrometer (Oriel MS260i). The desired : detector was selected with a flip mirror. We detected E ? -obj light either with a cooled, back-side thinned CCD camera (Andor DV440-BU2-2492) or by the same spectrometer-CCD combination just mentioned; selection was again mediated by a flip mirror. Lenses with f : l : = 300 mm were used to create the images on the cameras/spectrometer, giving an approximate effective magnification of 15. A Schott BG39 filter [35], with a transmittance &gt; 0 1 for 320 : &lt; GLYPH&lt;21&gt; &lt; 660 nm and &lt; 10 GLYPH&lt;0&gt; 4 for GLYPH&lt;21&gt; &lt; 320 nm and &gt; 780 nm (see Fig. 3 (a) bottom trace), was placed in the path before the detectors to minimize the laser light contribution to the images.

To facilitate a comparison between the direct RTS determination of the intensity and an indirect approach, we will analyze the focal spot image (Fig. 2 (a)). We captured this low-intensity image by reflecting a small portion of the beam, with a wedge beamsplitter after the OAP, into a third infinity-corrected microscope objective placed near the focus of the reflected beam. We

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00110v3-figure-37fd5fb8-7e84-4c29-9bfb-613d7cf78926.jpg)

Fig. 3. (a) Respectively from the top, RTS-spectral profiles (solid curves) captured by E k -obj for three values of U (2 67 : GLYPH&lt;6&gt; 0 08 : ; 1 73 : GLYPH&lt;6&gt; 0 05 : ; and 0 99 : GLYPH&lt;6&gt; 0 03 J) and the plasma : recombination lines (solid curve) along with the BG39 filter transmission (dashed red curve). The red arrows in the insets indicate GLYPH&lt;21&gt; ' 2 ' r as determined when the spectrum first reaches Sig bk while the dot-dash line indicates GLYPH&lt;21&gt; ' 2 ' r from the numerical fit to the shaded regions (see text). (b) Comparison between I RTS pk for the U values (black solid circles) and I Im pk for GLYPH&lt;1&gt; t = 30 fs (blue solid line) and 60 fs (red dashed line). The shaded areas are the uncertainties in I Im pk ( GLYPH&lt;24&gt; 16.1%); the dot-dashed lines are the uncertainties in I Im pk if A pix were known exactly ( GLYPH&lt;24&gt; 3.6%, see text). The open green triangle I RTS pk point was taken at U = 2 67 J but after the : low energy points (see text).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00110v3-figure-bd6c92c9-b3b4-4701-8507-65e52ab76cf0.jpg)

detected the magnified image of the focal spot with a Blackfly PGE 23S6C CMOS camera (having square pixels with 5.86 GLYPH&lt;22&gt; msides) after reducing the intensity further with neutral density filters. We imaged the spot size at a series of locations around the focus by adjusting the position of the microscope objective with a translation stage. We calibrated the Blackfly magnification by measuring the experimental focal spot directly (without the beamsplitter) with a Chamelon CCD camera (having square pixels with 3.75 GLYPH&lt;22&gt; msides) leading to 0 51 : GLYPH&lt;6&gt; 0 04 : GLYPH&lt;22&gt; m/pixel for the Blackfly camera.

## 3. Results and analysis

Before we discuss the analysis and how we can extract the peak intensity from composite images and Eq. 1, there are a few salient observations we can make. First, Fig. 2 (c) makes clear the need to gate the signal in order to isolate RTS and plasma recombination signals. The RTS radiation is prompt, only occurring during the laser pulse, while the plasma radiation begins a few tens of nanoseconds later and can be substantially stronger. Gating not only helps to minimize the recombination signal, it serves to block any laser light scattered from the chamber walls. This emphasizes the need for employing a large chamber so that the walls are far from the collection optics. It is clear when comparing the two signals in Fig. 2 (a) and (b), that only the central portion of the radiation captured by E k -obj makes it to the detector; the outer portion was clipped along the transfer path. We note that the weak signal above the primary signal in Fig. 2 (b) is due to residual, off-axis scattered light. To ensure E k -obj was viewing the focal spot, this objective was adjusted to maximize the on-axis signal. The signal from E ? -obj makes clear h is composed of multiple thinner lines of widths GLYPH&lt;24&gt; 15 GLYPH&lt;22&gt; m, closer to the anticipated focal width. The multiple lines are most likely caused by pointing instability on the order of 1 to 2 beam diameters; this instability appears to occur in both horizontal and vertical directions.

To that end, we display the prompt spectra for three U values (0 99 : GLYPH&lt;6&gt; 0 03 : ; 1 73 : GLYPH&lt;6&gt; 0 05 and 2 67 : : GLYPH&lt;6&gt; 0 08 J) in Fig. 3 (a). : We take GLYPH&lt;21&gt; ' 2 ' r to be the point where the spectrum falls to the background, Sigbk . For U = 2 67 J (Fig. 3 (a) top), : GLYPH&lt;21&gt; ' 2 ' r GLYPH&lt;24&gt; 700 nm. The values for GLYPH&lt;21&gt; ' 2 ' r are less clear for the other two energies, because of interference from the plasma recombination radiation. For these energies, we estimated GLYPH&lt;21&gt; ' 2 ' r two ways: (i) from the location where the curve falls to Sig bk and (ii) from a numerical fit of the tail of the distribution, identified by the shaded section, to a rational polynomial, f ' GLYPH&lt;21&gt; ' = a + b GLYPH&lt;157&gt;' GLYPH&lt;21&gt; + c ' with a , b and c being fitting parameters. In the latter we solve f ' GLYPH&lt;21&gt; ' = Sigbk to determine GLYPH&lt;21&gt; ' 2 ' r . This is also complicated by the presence of the plasma recombination lines distorting the shape, which is likely to cause the fit curve to decay to Sig bk prematurely. Because of the plasma contribution, the first method may select a GLYPH&lt;21&gt; ' 2 ' r that is too large at low intensities. Our estimates of GLYPH&lt;21&gt; ' 2 ' r for these intensities are averages of the two values. Even though GLYPH&lt;21&gt; ' 2 ' r is clear in the top trace of Fig. 3 (a), f ' GLYPH&lt;21&gt; ' GLYPH&lt;25&gt; Sigbk at a wavelength where the BG39 transmission is GLYPH&lt;24&gt; 1% and falling sharply. Thus, GLYPH&lt;21&gt; ' 2 ' r GLYPH&lt;24&gt; 700 nm may be a lower limit of the red shift. (We note that a numerical fit to the tail of the U = 2 67 J plot also gives : GLYPH&lt;21&gt; ' 2 ' r GLYPH&lt;24&gt; 700 nm.) Using GLYPH&lt;21&gt; 0 = 790 nm, obtained from a Wizzler [37] pulse characterization, we estimate I RTS pk ' GLYPH&lt;21&gt; ' 2 ' r ' = ' 1 71 : GLYPH&lt;6&gt; 0 40, 1 94 : : GLYPH&lt;6&gt; 0 41, 6 77 : : GLYPH&lt;6&gt; 0 72 : ' GLYPH&lt;2&gt; 10 18 W/cm 2 from lowest to highest U . The uncertainty in I RTS pk is based on half of the wavelength spread ( GLYPH&lt;24&gt; 39 2 nm GLYPH&lt;157&gt; = 19 5 nm) of the : laser, while the uncertainty in U is based on the RMS energy fluctuations of GLYPH&lt;24&gt; 2.65%.

Turning our attention to estimating the RTS intensity, we focus more closely on the prompt spectrum displayed in Fig. 2 (c). This spectrum contains information about the peak intensity and spatial intensity profile of the pulse, both of which are extractable, in principle, from Eq. 1. Determining the intensity profile, however, requires detailed knowledge of when in time and where in space each electron contributing to the measured RTS radiation is born (i.e., knowledge of the appearance intensities such as discussed in Ref. [30]), its trajectory induced by the field, and a reliable numerical RTS model providing quantitative information about its contribution to the captured light in various wavelength bands. With such, a functional fit of the data to Eq. 1 would be possible. Currently, a suitable model does not exist; however, this team is actively working on this problem. As a first step in developing a diagnostic tool, we will focus on the RTS peak intensity, I RTS pk , which corresponds to the red spectral edge, GLYPH&lt;21&gt; ' n ' r , in the prompt spectrum. We will confine our attention to n = 2. Because the n = 2 signal captured by E k -obj is larger than that captured by E ? -obj, as expected (see Fig. 5 in [36]), we will focus on composite images from E k -obj in the remainder of this manuscript.

We now compare the three I RTS pk values with I Im pk values obtained from Fig. 1 (a). There are three basic requirements for a good intensity estimate from such an image. First, to prevent image saturation so that the intensity profile can be measured over many decades, the camera must have a large dynamic range. Second, because the optical path contains lenses, windows and filters causing unknown magnification, the physical size of the camera pixels must be calibrated as mentioned earlier. Finally, assuming the spatial and temporal parts of the intensity profile are separable, an integrable temporal profile function is needed. When these conditions are met it is straightforward to show that the intensity associated with the energy U on target is given by

$$I _ { p k } ^ { \text{Im} } = \frac { \ U C _ { p k } } { T _ { \text{eff} } C _ { \text{sum} } A _ { \text{pix} } },$$

where C pk , T eff , C sum and A pix are, respectively, the count associated with the peak intensity, an effective pulse width, the total count summed of the pixels comprising the entire pulse and the area of a camera pixel in physical units. We derive this equation in the appendix. The value attached to C pk is the average value of nine pixels in a 3 GLYPH&lt;2&gt; 3 block of pixels near the the peak of the pulse on the image in Fig. 1 (a). To estimate T eff , we first integrated the Wizzler pulse profile, which was approximately Gaussian; however, the Wizzler measurement was made in the laser bay before beam transport to the target area. Because the pulse width estimated by the Wizzler ( GLYPH&lt;1&gt; t W) was smaller than that estimated by an autocorrelator ( GLYPH&lt;1&gt; t AC) in the target area by GLYPH&lt;24&gt; 20% (e.g., GLYPH&lt;1&gt; t W = 25 2 : GLYPH&lt;6&gt; 0 4 fs and : GLYPH&lt;1&gt; t AC = 31 0 : GLYPH&lt;6&gt; 0 4 fs), we rescaled : T eff by the ratio GLYPH&lt;1&gt; t GLYPH&lt;157&gt; GLYPH&lt;1&gt; t W, where GLYPH&lt;1&gt; t is a nominal pulse width of 30 fs or 60 fs in this case. The value of C sum is obtained by summing the pixels on the image that comprise the pulse (i.e., inside the white circle in Fig. 2 (a)). From the magnification calibration, A pix = 0 26 : GLYPH&lt;6&gt; 0 04 : GLYPH&lt;22&gt; m 2 . Our best estimates for I Im pk for GLYPH&lt;1&gt; t = 30 fs are ' 2 2 : GLYPH&lt;6&gt; 0 4 : ; 3 9 : GLYPH&lt;6&gt; 0 6 and 6 0 : : GLYPH&lt;6&gt; 1 0 : ' GLYPH&lt;2&gt; 10 18 W/cm 2 from lowest to highest U .

## 4. Discussion

We compare I RTS pk and I Im pk graphically in Fig. 3 (b), from which we can draw several conclusions. First, given the uncertainties in I Im pk and I RTS pk , the two are in general agreement. Second, the uncertainty in I Im pk is quite large ( GLYPH&lt;24&gt; 16 1%), with the largest contribution coming from the : calibration of the camera ( GLYPH&lt;24&gt; 15.7%). Consequently, a higher precision comparison is possible if an improved pixel area calibration is made. It is interesting to note that the largest and smallest I RTS pk datum points are in good statistical agreement with I Im pk (30 fs), even ignoring the A pix uncertainty; I RTS pk ( U = 2 67 J : ' &gt; I Im pk (30 fs) while I RTS pk ( U = 0 99 J : ' &lt; I Im pk (30 fs). At the same time, I RTS pk ( U = 1 73 J) is lower than : I Im pk (30 fs) by more than 3 GLYPH&lt;27&gt; , and consistent with I Im pk (60 fs). There are three contributions to this observation. First, as mentioned earlier, the spectra are contaminated with plasma recombination lines, making the determination of GLYPH&lt;21&gt; ' 2 ' r difficult; the observed onset wavelengths might be shifted to lower values. Second, the I RTS pk measurements were taken over several hours in the following order: U = 2 67, 0 99, and 1 73 J. The pulse width of this laser, : : : which will affect the peak intensity, is known to drift with time; we had no capacity to monitor the pulse width in the target area in real time so the laser was not re-optimized after the initial optimization at the beginning of the day. Third, changing the pulse energy typically requires re-optimization for peak performance. Consequently, it is not surprising that the intensity was lower later in the day and after adjusting the energy. To check for possible issues of this kind, we took a final I RTS pk measurement at U = 2 67 J after the : U = 1 73 J measurement, the green : triangle in Fig. 3 (b). The fact that the green U = 2 67 J point is lower than the black point is : supporting evidence that the lower intensities after energy adjustment is likely caused by the laser not being optimized. Evidently I RTS pk is sensitive to changes in the beam properties leading to changes in the intensity, whereas I Im pk is not. While our results suggest the indirect image and RTS methods are in reasonable agreement, they also point out the need for real-time monitoring of the laser characteristics if the former is employed for determining the peak intensity.

Our study clearly shows that the RTS approach has promise as a diagnostic tool. Nevertheless, there are issues that demand further investigation. The most critical of these is a lack of knowledge of the conditions necessary to ensure that the GLYPH&lt;21&gt; ' n ' r observed is associated with Thomson scattering of electrons experiencing the peak intensity of the pulse. Electrons born on axis near the peak of the pulse would seem to have the best chance of contributing to the RTS signal associated with the peak intensity. An electron born off axis with near zero velocity will not experience the peak intensity because (1) the ponderomotive energy ( U p ) it picks up from the field will at best propel it no closer to the axis than its birth position and (2) the ponderomotive force tends to push it further and further away from the axis with each oscillation [26]. Many, electrons fall into this category. It is unclear if electrons born on the axis, but early in time before the pulse reaches its peak intensity, will remain close to the axis long enough to experience the peak intensity or be scattered away from the axis by the increasing ponderomotive force. More subtle issues associated with the carrier phase and the time of birth of the electron also may play a role. For example, electrons born either stationary at a peak of an optical cycle or with an initial velocity transverse to the laser axis at the right time during a cycle so that it is turned around by the field, both have a chance to remain near the laser axis as they undergo Thomson scattering. (Only the initial transverse velocity of the electron is important because the initial longitudinal velocity does not affect whether the electron is scattered away from the axis.) The dependence on phase has been discussed theoretically; see, for example, [25] and references therein. To summarize, many scattered electrons either (1) leave the beam axis before the peak of the pulse envelope, never experiencing the peak intensity or leave the axis before the pulse is complete. Either way many electrons may contribute Doppler-shifted components to the RTS spectra that are smaller than expected from a plane-wave analysis. Consequently, I RTS pk might only be a lower limit to the true intensity. In order to understand these issues more fully, and to determine which RTS electrons properly reflect the peak intensity, both a numerical model and more experimental investigations are required.

## 5. Conclusions

The measurements reported in this manuscript were performed with N2. In addition to N 2 yielding ten electrons/molecule, the electrons are liberated at several different intensities, a few near the expected peak intensity. This gives us confidence that at least a few electrons are present when the intensity is at its peak. As argued above, electrons liberated near the peak intensity may have a better chance of experiencing the peak intensity. Measuring the RTS signal for different gases with different ionization potentials and number of electrons liberated (e.g., He, Ar, N2, etc.) could provide a better understanding of the role the ponderomotive force plays in the RTS, particularly for steep field gradients. It might also be possible to perform a pump/probe investigation to determine the time history of the electrons.

Capturing the RTS signal at ' GLYPH&lt;18&gt;; GLYPH&lt;30&gt; ' = ' GLYPH&lt;25&gt; GLYPH&lt;157&gt; 2 0 ; ' is most useful between 10 18 and 10 19 W/cm 2 . At some point, which we estimate to be GLYPH&lt;24&gt; 10 19 W/cm 2 , the harmonics will start to overlap, making it difficult to extract GLYPH&lt;21&gt; ' n ' r as we have done. It is clear from inverting Eq. 1 that the overlap between harmonics can be slowed by observing the radiation at smaller GLYPH&lt;18&gt; . While additional information about the RTS is contained in the angular distribution of the radiation, as suggested by Harvey [29], this approach will run into issues when GLYPH&lt;18&gt; gets too small, making it difficult to prevent the laser light from being captured with the RTS radiation. We estimate that this will occur before 10 21 W/cm 2 . A second possibility would be to pick off a calibrated portion of the energy and always measure the intensity between 10 18 and 10 19 W/cm 2 . A combination of picking off GLYPH&lt;24&gt; 1 % and capturing the radiation at smaller angles should allow intensities to be measured up to about 10 23 W/cm 2 . At this stage, quantum effects and radiation reaction will need to be considered. A third approach is to monitor the infrared ( n = 1) directly. As the electrons become more relativistic, GLYPH&lt;18&gt; again becomes smaller. Between 10 23 and 10 25 W/cm 2 , it is not clear how to implement the red shifts with initially cold electrons to monitor the intensities. The more complicated blue-shifts from higher-order Thomson scattering may need to be employed [38]. As the intensity approaches 10 25 W/cm 2 , a new paradigm for characterizing the intensity based on proton RTS may open. When an effective normalized vector potential for protons approaches unity (i.e., a p = eE GLYPH&lt;21&gt; GLYPH&lt;157&gt; 2 GLYPH&lt;25&gt; m c p 2 GLYPH&lt;24&gt; 1), the motion of free protons will become relativistic in the laser field and could emit radiation similar to that emitted by the electron, but for intensities about ' m p GLYPH&lt;157&gt; m e ' 2 times larger than the electron case. To our knowledge, this has never been discussed in the literature.

## Appendix

In this appendix, we derive of Eq. 2. We start by defining the intensity profile, which is a function of time and space in the focal plane,

$$I ( t, x, y ) = I _ { 0 } g ( t ) f ( x, y ), & & ( \text{A1} )$$

where I 0 is the peak intensity ( I 0 GLYPH&lt;17&gt; I ' t 0 ; x 0 ; y 0 ' , with t 0 , x 0 , and y 0 being the time and position of peak intensity respectively); g ' ' t and f ' x ; y ' are respectively the temporal and spatial shape functions with g ' t 0 ' = f ' x 0 ; y 0 ' = 1. The total pulse energy on target is

$$U \equiv \iiint I ( t, x, y ) d t d x d y = I _ { 0 } \int g ( t ) d t \iint f ( x, y ) d x d y.$$

We define the time integral of g ' ' t as an effective pulse width,

$$T _ { \text{eff} } & = \int g ( t ) d t, & ( \text{A3} )$$

where T eff / GLYPH&lt;1&gt; t ' FWHM . ' For a Gaussian shape ( g ' ' t = e GLYPH&lt;0&gt; 4 ln 2 ' t GLYPH&lt;157&gt; GLYPH&lt;1&gt; t ' 2 ), for example, T eff = p GLYPH&lt;25&gt; GLYPH&lt;157&gt; ln 2 GLYPH&lt;1&gt; t GLYPH&lt;157&gt; 2 GLYPH&lt;25&gt; 1 064 : GLYPH&lt;1&gt; t . We next convert the spatial integral over f ' x ; y ' to a discrete sum,

$$\iint f ( x, y ) d x d y & \approx \sum _ { i, j } f ( x _ { i }, y _ { j } ) \Delta x \Delta y, & ( A 4 )$$

where GLYPH&lt;1&gt; x GLYPH&lt;1&gt; y corresponds to A pix , the area in the focal plane captured by a pixel. We identify f ' x i ; y j ' as the count at pixel ' x i ; y j ' , divided by the peak count C pk . We determined C pk from the focal spot image by taking the set of nine pixels (3 GLYPH&lt;2&gt; 3 region) near the vicinity of the peak with the largest total count, and letting C pk be equal to the average of the nine pixels. Summing over all pixel comprising the pulse, we therefore obtain

$$\sum _ { i, j } f ( x _ { i }, y _ { i } ) = \frac { C _ { \text{sum} } } { C _ { \text{pk} } },$$

where C sum is the sum of the counts of the image pixels comprising the pulse. Thus, we have

$$U = I _ { 0 } T _ { \text{eff} } \frac { C _ { \text{sum} } } { C _ { \text{pk} } } A _ { \text{pix} }.$$

Solving for I 0 = I Im pk , we arrive at Eq. 3.

## Funding

National Science Foundation (PHY1806584); Natural Sciences and Engineering Research Council of Canada (NSERC) (RGPIN-2014-05736); Spanish Ministerio de Ciencia, Innovacion y Universidades (PALMA Grant FIS2016-81056-R); LaserLab Europe IV (654148); Junta de Castilla y León Consolidated Research, Junta de Castilla y Leon (167).

## Acknowledgements

We thank Enam Chowdhury for helpful discussion regarding proton RTS; Mauricio Rico, Diego de Luis, Juan Hernandez, Diego Arana and Jon Imanol Apiñaniz Aginako for technical help before and during the experiment; Cruz Mendez and the VEGA laser team for running the laser during the experiment; Emma McMillian, Max Szostak and Jacob Smith for help with the data analysis.

## References

- 1. W. T. Hill, III and L. Roso, 'Probing the quantum vacuum with petawatt lasers,' J. Phys.: Conf. Ser. 869 , 012015 (2017).
- 3. G. V. Dunne, 'New strong-field QED effects at extreme light infrastructure,' Eur. Phys. J. D 55 , 327-340 (2009).
- 2. D. Tommasini, A. Ferrando, H. Michinel, and M. Seco, 'Precision tests of QED and non-standard models by searching photon-photon scattering in vacuum with high power lasers,' J. High Energy Phys. 2009 , 043-043 (2009).
- 4. National Academies of Sciences, Engineering, and Medicine, Opportunities in Intense Ultrafast Lasers: Reaching for the Brightest Light (The National Academies Press, Washington, DC, 2018).
- 6. H. Euler and B. Kockel, 'Über die Streuung von Licht an Licht nach der Diracschen Theorie (The Scattering of Light by Light in the Dirac Theory),' Naturwissenschaften 23 , 246-247 (1935).
- 5. O. Halpern, 'Scattering Processes Produced by Electrons in Negative Energy States,' Phys. Rev. 44 , 855-856 (1933).
- 7. H. Euler, 'Über die Streuung von Licht an Licht nach der Diracschen Theorie (On the Scattering of Light by Light in Dirac's Theory),' Ann. Phys. 418 , 398-448 (1936).
- 9. T. Heinzl, B. Liesfeld, K.-U. Amthor, H. Schwoerer, R. Sauerbrey, and A. Wipf, 'On the observation of vacuum birefringence,' Opt. Commun. 267 , 318-321 (2006).
- 8. Y. Liang and A. Czarnecki, 'Photon-photon scattering: a tutorial,' Can. J. Phys. 90 , 11-16 (2012).
- 10. D. Tommasini, D. Novoa, and L. Roso, 'Quantum vacuum polarization searches with high power lasers below the pair production regime,' in Progress in Ultrafast Intense Laser Science X, vol. 106 of Springer Series in Chemical Physics K. Yamanouchi, G. Paulus, and D. Mathur, eds. (Springer Science, 2014), pp. 137-153.
- 12. M. Roth, D. Jung, K. Falk, N. Guler, O. Deppert, M. Devlin, A. Favalli, J. Fernandez, D. Gautier, M. Geissel, R. Haight, C. E. Hamiltonand, B. M. Hegelich, R. P. Johnson, F. Merrill, G. Schaumann, K. Schoenberg, M. Schollmeier, T. Shimada, T. Taddeucci, J. L. Tybo, F. Wagner, S. A. Wender, C. H. Wilde, and G. A. Wurden, 'Bright Laser-Driven Neutron Source Based on the Relativistic Transparency of Solids,' Phys. Rev. Lett. 110 , 044802 (2013).
- 11. A. Blinne, H. Gies, F. Karbstein, C. Kohlfürst, and M. Zepf, 'All-optical signatures of quantum vacuum nonlinearities in generic laser fields,' Phys. Rev. D 99 , 016006 (2019).
- 13. G. M. Dyer, A. C. Bernstein, B. I. Cho, J. Osterholz, W. Grigsby, A. Dalton, R. Shepherd, Y. Ping, H. Chen, K. Widmann, and T. Ditmire, 'Equation of State Measurement of Dense Plasmas Heated with Fast Protons,' Phys. Rev. Lett. 101 , 015002 (2008).
- 15. S. Feldman, G. Dyer, D. Kuk, and T. Ditmire, 'Measurement of the equation of state of solid-density copper heated with laser-accelerated protons,' Phys. Rev. E 95 , 031201(R) (2017).
- 14. M. Roth, I. Alber, V. Bagnoud, C. R. D. Brown, R. Clarke, H. Daido, J. Fernandez, K. Flippo, S. Gaillard, C. Gauthier, M. Geissel, S. Glenzer, G. Gregori, M. Günther, K. Harres, R. Heathcote, A. Kritcher, N. Kugland, S. LePape, B. Li, M. Makita, J. Mithen, C. Niemann, F. Nürnberg, D. Offermann, A. Otten, A. Pelka, D. Riley, G. Schaumann, M. Schollmeier, J. Schütrumpf, M. Tampo, A. Tauschwitz, and A. Tauschwitz, 'Proton acceleration experiments and warm dense matter research using high power lasers,' Plasma Phys. Control. Fusion 51 , 124039 (2009).
- 16. M. Roth, A. Blazevic, M. Geissel, T. Schlegel, T. E. Cowan, M. Allen, J.-C. Gauthier, P. Audebert, J. Fuchs, J. Meyer-ter-Vehn, M. Hegelich, S. Karsch, and A. Pukhov, 'Energetic ions generated by laser pulses: A detailed study on target properties,' Phys. Rev. ST Accel. Beams 5 , 061301 (2002).
- 18. E. Lefebvre, 'Numerical simulation of isotope production for positron emission tomography with laser-accelerated ions,' J. Appl. Phys. 100 , 113308 (2006).
- 17. B. Jones, 'Radiotherapy for the future: Protons and ions hold much promise,' BMJ 330 , 979-980 (2005).

19. K. W. D. Ledingham, W. Galster, and R. Sauerbrey, 'Laser-driven proton oncology - a unique new cancer therapy?'

80

Br. J. Radiol.

, 855-858 (2007).

- 21. F. Tommasino and M. Durante, 'Proton Radiobiology,' Cancers 7 , 353-381 (2015).
- 20. T.-C. Liu, X. Shao, C.-S. Liu, B. Eliasson, W. T. Hill, III, J. Wang, and S.-H. Chen, 'Laser acceleration of protons using multi-ion plasma gaseous targets,' New J. Phys. 17 , 023018 (2015).
- 22. S. Chen, A. Maksimchuk, and D. Umstadter, 'Experimental observation of relativistic nonlinear Thomson scattering,' Nature 396 , 653-655 (1998).
- 24. C. I. Castillo-Herrera and T. W. Johnston, 'Incoherent harmonic emission from strong electromagnetic waves in plasmas,' IEEE Trans. Plasma Sci. 21 , 125-135 (1993).
- 23. E. S. Sarachick and G. T. Schappert, 'Classical Theory of the Scattering of Intense Laser Radiation by Free Electrons,' Phys. Rev. D 1 , 2738-2752 (1970).
- 25. J. Gao, 'Thomson Scattering from Ultrashort and Ultraintense Laser Pulses,' Phys. Rev. Lett. 93 , 243001 (2004).
- 27. O. Har-Shemesh and A. Piazza, 'Peak intensity measurement of relativistic lasers via nonlinear Thomson scattering,' Opt. Lett. 37 , 1352-1354 (2012).
- 26. J. Gao, 'Laser intensity measurement by Thomson scattering,' Appl. Phys. Lett. 88 , 091105 (2006).
- 28. G. Tarbox, E. Cunningham, R. Sandberg, J. Peatross, and M. Ware, 'Radiation from free electrons in a laser focus at 10 18 w/cm 2 : modeling of photon yields and required focal conditions,' J. Opt. Soc. Am. B 32 , 743-750 (2015).
- 30. M. F. Ciappina, S. V. Popruzhenko, S. V. Bulanov, T. Ditmire, G. Korn, and S. Weber, 'Progress toward atomic diagnostics of ultrahigh laser intensities,' Phys. Rev. A 99 , 043405 (2019).
- 29. C. N. Harvey, ' In situ characterization of ultraintense laser pulses,' Phys. Rev. Accel. Beams 21 , 114001 (2018).
- 31. Centro de Láseres Pulsados (CLPU), https://www.clpu.es/en, (2019).
- 32. T. Zuo and A. D. Bandrauk, 'Charge-resonance-enhanced ionization of diatomic molecular ions by intense lasers,'

Phys. Rev. A 52 , 2511-2514 (1995).

- 34. X. M. Tong, Z. X. Zhao, and C. D. Lin, 'Theory of molecular tunneling ionization,' Phys. Rev. A 66 , 033402 (2002).
- 33. M. Ammosov, N. Delone, and V. Krainov, 'Tunnel ionization of complex atoms and of atomic ions in an alternating electromagnetic field,' J. Exp. Theor. Phys. 64 , 1191-1194 (1986).
- 35. Schott, 'BG39 Data Sheet,' https://www.sydor.com/wp-content/uploads/ SCHOTT-BG39-Bandpass-Filter.pdf, (2019).
- 37. FASTLITE, 'Wizzler,' http://www.fastlite.com/en/ar824651-823532-Wizzler.html, (2019).
- 36. A. Paredes, D. Novoa, and D. Tommasini, 'Measuring extreme vacuum pressure with ultraintense lasers,' Phys. Rev. Lett. 109 , 253903 (2012).
- 38. W. Yan, C. Fruhling, G. Golovin, D. Haden, J. Luo, P. Zhang, B. Zhao, J. Zhang, C. Liu, M. Chen, S. Chen, S. Banerjee, and D. Umstadter, 'High-order multiphoton Thomson scattering,' Nat. Photonics 11 , 514-520 (2017).
C. Z. He A. Longman J. A. Pérez-Hernández M. De Marco C. Salgado G. Zeraouli G. Gatti L. Roso R. Fedosejevs W. T. Hill III
Towards an in situ, full-power gauge of the focal-volume intensity of petawatt-class lasers
About 50 years ago, Sarachick and Schappert [Phys. Rev. D. 1, 2738 (1970)]
showed that relativistic Thomson scattering leads to wavelength shifts that are
proportional to the laser intensity. About 28 years later Chen et al. [Nature
396, 653 (1998)] used these shifts to estimate their laser intensity near 10^18
W/cm^2. More recently there have been several theoretical studies aimed at
exploiting nonlinear Thomson scattering as a tool for direct measurement of
intensities well into the relativistic regime. We present the first
quantitative study of this approach for intensities between 10^18 and 10^19
W/cm^2. We show that the spectral shifts are in reasonable agreement with
estimates of the peak intensity extracted from images of the focal area
obtained at reduced power. Finally, we discuss the viability of the approach,
its range of usefulness and how it might be extended to gauge intensities well
in excess of 10^19 W/cm^2.