1. Introduction
Currently, significant efforts are devoted to the creation and development of laser systems with both high pulse energy and a high pulse repetition rate. In many cases, such systems include amplifiers based on media doped with Yb
3+ ions. This is due to the fact that Yb
3+ ions have a low quantum defect of about 9%. In addition, many widely used active media, such as Yb:YAG, Yb:CaF
2, Yb:YLF [
1,
2], and other materials doped with Yb
3+ ions, have good thermal conductivity, which increases significantly when the active element is cooled to cryogenic temperatures. Nevertheless, it is not possible to avoid the negative influence of thermal effects if the average pump radiation power is high enough. There are methods for a further decrease in a thermal load on active elements, for example, the movement of a pumped region in an active medium [
3] or the application of several spatially separated slabs with different doping levels [
4]. Another approach is the application of active elements with an ion doping concentration gradient in the direction of pump propagation. Gradient doping makes it possible to optimize the temperature distribution in the active element, thus reducing a thermally induced wavefront distortion of the amplified radiation [
5].
In laser systems with high pump radiation power, changes in the wavefront of the amplified radiation in the active elements lead to the appearance of pump-induced lenses [
6,
7] as well as to the depolarization of radiation [
8,
9]. These effects are especially detrimental in the case of multielement and multipass systems. The wavefront distortion of amplified radiation limits the output energy and the maximum reachable diameter of the focused radiation, which reduces the peak power and intensity [
10,
11]. Beam wavefront has a great influence on the performance of coherent beam combining [
12]. The compensation of distortion caused by pumping the active element in laser systems is very important when scaling to multi-kilowatt average power levels. Taking this into account, wavefront distortion is necessary when developing the element base of a laser system [
5] and tailoring pump and injection radiation parameters [
13].
Typically, a Shack–Hartmann sensor [
14,
15,
16] is used to measure wavefront distortion. This is a well-known, simple, and commercially available device. There are also wavefront measurement sensors based on lateral shearing interferometry [
17]. Research using interferometers [
18,
19] is also common. For these approaches, the stability of the probe laser beam wavefront and intensity is most important. To separate the processes of wavefront distortion and the processes of changes in the spatial distribution of radiation intensity due to amplification or absorption, it is necessary to use probe laser radiation that is not absorbed or amplified in the active medium. The wavelength of the scanning radiation must also be in accordance within the spectral range for which the dielectric coatings of the optical elements are designed. Since the parameters of the medium and introduced distortion are wavelength-dependent, it is preferable to carry out the study at the amplification wavelength. In the case of measurements using interferometric methods, the disadvantage is also the complex setup and sensitivity to a drift of various mechanical units of the system [
20]. Mechanical issues are especially significant in high-power laser systems due to the operation of the cooling system of the active elements.
If high stability of radiation parameters is required, it is also necessary to separate the contributions of various mechanisms in order to compensate them. However, in most cases, the methods discussed above do not allow decomposing the change in the refractive index caused by the thermal effects and polarizability difference in the excited and ground states of the active ions.
Wavefront distortion in pumped active elements can be calculated from the experimentally determined distribution of temperature and population density of the upper laser level in the active element [
18,
21]. On the one hand, with such a calculation, the results’ accuracy depends on the accuracy of the laser and thermophysical characteristics of the medium. On the other hand, such a study makes it possible to distinguish between the thermal contribution to the change in the refractive index and the electronic one. In addition, such a calculation essentially complements the experimental studies of the pump radiation effect on the active medium. Such experimental studies are necessary to verify numerical models of temperature distribution and radiation amplification in active elements and estimate the quality of the thermal contact between the active element and the heat sink.
Thus, progress in the development of modern laser systems with high radiation power requires the development of new methods for studying the interaction of pump radiation, active medium, and amplified radiation. Currently, one of the directions of research is the development of novel techniques applicable in the case of cryogenic cooling, for the experimental study of temperature distribution in active elements and distortion of the wavefront of amplified radiation [
21,
22].
In this work, a laser scanning method was proposed for a simultaneous time-resolved study of the temperature distribution and the upper level population density in the active element of a high-power laser amplifier. The accuracy of the method was analyzed. The method was applied to study the wavefront distortion in a crystalline Yb:YAG active mirror of a cryogenically cooled multi-disk laser amplifier with a diode pump per element up to 200 W [
13]. The thermal and electronic contributions to optical path difference (OPD) profiles introduced by the pumped active elements to the wavefront of amplified radiation were experimentally determined and analyzed.
2. Laser Scanning Method
There are two contributions to the distortion of radiation wavefront in active elements. Thermal contribution (Δ
lr,θ(th)) is associated with a temperature gradient in an active element, and electronic contribution (Δ
l(el)) occurs because of a change in the polarizability of excited dopant ions. Thermal contribution to the change in the optical path length for radiation components with radial and tangential polarizations are different and are indicated by
r and
θ, respectively. OPD can be calculated with respect to any region of an active element. So, the center of a pumped area in an active element (
x0,
y0) is chosen. For the calculation of wavefront distortion caused by the temperature gradient in an active element, a common approach is the application of the plane strain approximation (for thin rods) and the plane stress approximation (for thin disks). Thus, the plane stress approximation is valid for the case of active elements in the form of active mirrors, which are often used in laser amplifiers with both high average and high peak power levels. The electronic contribution to OPD and thermal contribution calculated in the plane stress approximation can be described by the following formulae [
23,
24]:
where
x,
y are the Cartesian coordinates in the plane perpendicular to the
z-axis;
L is the optical path length in an active element, m;
is the partial derivative of refractive index with respect to temperature taken at a zero strain, K
−1;
n0 is the refractive index of the medium before pumping is turned on;
αT(
T) is the thermal expansion coefficient, K
−1;
Cr,θ are photoelastic constants for radiation components with radial and tangential polarizations, respectively;
ν is Poisson’s ratio;
is the temperature difference averaged along the beam propagation axis (
z-axis), K; Δ
p is the polarizability difference in dopant ions in the excited and ground states, m
3; and
N(
x,
y) is the population of the excited laser state, ions/m
3.
In Equation (1), different terms account for different mechanisms involved into the formation of wavefront distortion of amplified radiation. The first term is associated with the refractive index dependence on temperature; the term with χbg takes into account the end faces bulging caused by the thermal expansion of the active medium. The term with χpe is related to thermal stresses in the medium resulting in the induced birefringence due to the photoelastic effect.
Since the partial derivative of the refractive index taken at zero stress
, was measured directly in the published data, the following relation should be used [
25]:
where p
11 and p
12 are the coefficients related to the components of a fourth-rank elasto-optical tensor. The population of the exited laser state can be easily obtained from the experimentally investigated gain coefficient
G (
x,
y) [
26]:
where
σL(T) is the luminescence cross-section at the gain wavelength, m
2;
Ntot is the dopant ion concentration, ions/m
3; and f
11(T) and
f03(T) are Boltzmann occupation factors.
In formulae 1–4, the dependence of laser and thermophysical parameters of a gain medium is taken into account. Generally, elasto-optical coefficients
pij and, calculated from them, photoelastic constants
Cr,θ are temperature-dependent [
27]. But, according to our knowledge, there are no measurements of temperature dependence or even values of elasto-optical coefficients at cryogenic temperatures for Yb:YAG crystals or ceramics. So, in our calculations, we used elasto-optical coefficients and photoelastic constants determined at room temperature. There was a polarizability difference Δ
p at room temperatures [
28] and cryogenic temperatures [
29] in earlier published data, and the values differed by approximately 2 times. Thus, the polarizability difference is strongly dependent on the temperature value. To the best of our knowledge, there are no published data concerning the dependence of the polarizability difference on temperature. So, it is considered as a constant, and the value 88 × 10
−34 m
2 reported in [
29] was used.
The distributions of temperature and population density of the excited laser state were determined experimentally in the pump region in accordance with the laser scanning technique developed by the authors. The method is a development of the previously reported method for studying the temperature in the active elements of laser amplifiers [
22]. The method allows one to carry out a simultaneous experimental investigation of temperature and population density based on measuring the transmittance and gain in the active element. The method can be used in the case of an active medium in which the absorption cross-section at the gain wavelength is temperature-dependent. However, if it is necessary to measure only the temperature distribution, it is sufficient if a temperature dependence of the absorption cross-section at the wavelength of the scanning laser is present.
The experimental setup for scanning an active element is shown in
Figure 1.
Such an experimental setup can be used for different geometries of active elements, including active mirrors and slabs. In the experiment, the laser radiation at the gain wavelength was focused into an active element and translated in parallel by movable mirrors M
2 and M
3 to scan a certain region in the active element. The radiation passed through the lens (L
1) before being reflected by the mirrors for scanning. The first mirror (M
2) displaced the beam in the vertical direction, and the second one (M
3), in the horizontal direction. The horizontal plane was parallel to the plane of the optical table. To accelerate the scanning process, the motion of translation stages, on which mirrors were mounted, was automated. The setup also contained a CCD camera to monitor the stability of the beam position during scanning and a photodetector to monitor the scanning radiation power. A polarizing beamsplitter (PBS), photodetector I, and photodetector II were used to measure the change in the polarization of scanning radiation in an active element (described in the section titled Results). For the measurement of temperature and gain coefficient, the experimental setup without the PBS and photodetector II was used. Such a measurement can be carried out using a CCD camera [
30] and a large-diameter probe beam without scanning. In this case, the wavefront distortion introduced by an active element and entire optic system affected the measured absorption and gain profiles. So, laser scanning allows one to measure the temperature and gain coefficient profiles without the influence of wavefront distortions. Also, a CCD camera typically has a much lower temporal resolution compared to photodiode-based detectors.
The investigation algorithm was as follows. The transmittance coefficient dependence on the active element temperature was investigated in the absence of pump radiation , where PI and Pc are the radiation power measured by the photodetector I and control photodetector, respectively. The dependence was then calculated from it.
So, since not only the active element was scanned by the laser beam but also the mirrors, lens, and the photosensitive active regions of the photodetector, it was necessary to take into account temperature-independent losses in the optic system. For this purpose, the transmittance coefficient distribution was measured in a cooled active element without pumping.
Consider the situation when the active element is pumped. In the case of a high pulse repetition rate or large duty cycle, the upper laser level does not have time to completely decay between pump pulses. If so, the pump radiation should be interrupted for a period several times longer than the lifetime of the upper-level state. This allows one to measure the absorption coefficient not affected by the pumping radiation. Since the average radiation power decreases when the pump radiation is interrupted, the pump interruption frequency should be chosen so as to have a minimal effect on the steady-state temperature field in the active element. In our studies, the pumping process was interrupted 2 times per second for 10 ms, and it corresponded to a decrease in the average pump radiation power by 2%. The transmittance distribution
and gain distribution
in the active element were investigated. A typical signal in the vicinity of a pump interruption is shown in
Figure 2. In this case, pump pulses had a repetition rate of 1 kHz and pulse duration 0.6 ms. The scanning radiation was continuous wave. To avoid numerical errors, the signal was fitted with a decaying exponent.
The investigated distributions were divided to the transmittance distribution in the cooled active element (from step 2)
and
, where
. The equivalent temperature was determined from the measured transmittance using the dependence
(from step 1). The equivalent temperature was the temperature of the uniformly heated active element with the same transmittance as the measured one [
22]. The temperature dependence of the absorption coefficient in the element is specified as follows:
Then, the equivalent temperature is:
where
A−1 is the inverse function of
A.
For an arbitrary temperature distribution and absorption coefficient temperature dependence, the equivalent temperature
Tc and average temperature
Ta are different. If the temperature is constant along the
z-axis in the active element, then
Tc and
Ta are the same. Moreover, if absorption depends linearly on temperature, the quantities are also the same regardless of the temperature distribution. Let the absorption coefficient per unit length to be written as
, where
a,
b are arbitrary real-valued coefficients. Then, the absorption in the element over the path
L is
This relation is held for an arbitrary temperature distribution.
In case the absorption coefficient does not comply with the linearity requirement, then the equivalent temperature, as well as the average temperature, lies in the range between minimum and maximum temperatures. In the Yb:YAG crystal, there was an almost linear region of the absorption coefficient (
Figure 3).
As can be seen, for the temperature values in the range of 180–280 K, the dependence was almost linear and the maximal relative error of fit was less than 7%. In this region, the temperature measured by the laser scanning method was, therefore, the average temperature, which was conventionally used to calculate the wavefront distortion.
The uncertainty was calculated by the means of error propagation formula:
where
σk is the transmittance standard deviation,
σC is the standard deviation of the transmittance measured in the absence of pump radiation,
σg is the standard deviation of the gain coefficient, and
ts is the coefficient calculated using the Student distribution with a 97.5% confidence interval. To determine the number of repeated measurements at each point on the map, the dependence of the error in determining the temperature on the number of transmittance measurements in different regions of the active element was estimated in different map points. A typical dependence is presented in
Figure 4.
The error in temperature determining decreased as 1/N, where N is the number of measurements. For the measurement presented in the graph, the minimum value of the temperature determination error based on the curve fitting results was 1.37 K. The non-zero limit value was due to the influence of the instrumental error of the recording system. In our study, the number of repeated measurements was chosen to be 15. Therefore, the degree of freedom used for the calculation was 14 for the temperature measurements. The number of peak gain coefficient measurements in a single map point was more than 10,000.
3. Results
To verify the method, the study of the amplified radiation wavefront distortion was carried out. The distortion was measured in one of the eight active elements of a 16-pass laser amplifier with cryogenic cooling and an average pump radiation power of up to 200 W [
13]. The active elements were diffusion-bonded, disk-shaped YAG-Yb:YAG crystals (Yb
3+ ion concentration 9.8 at.%), which were attached in pairs on the opposite sides to cryogenically cooled copper holders. The disk diameter was 25.4 mm, the thickness of the doped part was 3.75 mm, and the thickness of the undoped part was 2 mm. During the laser scanning, absorption and amplification occurred only in the doped part. Accordingly, when calculating the OPD, the optical radiation path length was taken into account only in the doped part. The undoped part did not affect the electronic component of the OPD due to the lack of an upper-level population distribution. According to the simulated data, the temperature distribution in the undoped part had significantly smaller gradients in comparison with the doped part, and, therefore, it induced a negligible effect on the result of the OPD calculation.
A continuous wave laser with a central wavelength of 1030 nm and a spectral width of less than 0.3 nm was used for scanning. The beam diameter in the active element was 220 ± 14 µm. The angle between the direction of propagation of the scanning radiation and the normal to the surface was 9.5° in the horizontal plane. The pump beam had a hyper-Gaussian intensity profile, the pump beam diameter was 2 mm (1/e
2 level), pump pulses had a square shape in the time domain, the pulse repetition rate was 1000 Hz, the peak power was 180 W, the average absorbed power of the pump radiation was 110 W, and the duty cycle was 60%. The experimental scheme is shown in
Figure 1.
The field with a diameter of 4.4 mm was scanned. The maximum measured temperature was 192 ± 0.6 K, and the temperature at the edges of the scanned area was ~165 K. The error in determining the temperature did not exceed 2.5 K. The maximum gain was 2.9 ± 0.07, and the error in determining the gain did not exceed 0.09. The thermal contribution determination error took into account the temperature determination error. The electronic contribution was calculated using the temperature error as well as the gain coefficient error. In
Figure 5A,B are the thermal and electronic contributions to the OPD profile introduced by the active element. The corresponding horizontal sections of the OPD are represented in
Figure 5C and 5D, respectively.
The electronic contribution to the OPD was less than the temperature contribution by 2 orders of magnitude, which was in agreement with the data available for laser systems with a high average power [
21,
29]. For systems with low pulse repetition rates and high peak power, the electronic contribution to the wavefront distortion was comparable to the thermal contribution and may even have exceeded it [
18].
The OPD profiles were analyzed using Zernike polynomials. The thermal profile of the OPD was defined mostly by the following terms: vertical tilt, 51.1 nm; defocus, 498.2 nm; primary spherical aberration, 171.4 nm; and second spherical aberration, 65.5 nm. The electronic profile was almost not tilted; so, the next terms were significant in the electronic profile: defocus, 4.9 nm; primary spherical aberration, 4.2 nm; and second spherical aberration, 2.2 nm. Since the OPD contained significant high-order coefficients, namely, the second spherical aberration, the equivalent focal length estimated using Zernike coefficients could have been inaccurate. To determine the focal length of the thermal and electron lenses, the corresponding contributions to the OPD were considered as the effect of a thin spherical lens and approximated by an elliptic paraboloid in the region limited by the pump beam diameter of 2 mm. Thus, the focal lengths of the thermal lens were
fx = 0.40 ± 0.03 m and
fy = 0.47 ± 0.05 m, and the focal lengths of the electronic length were estimated to be
fx = 21.3 ± 1.9 m and
fy = 20.9 ± 1.7 m. For the thermal OPD profile, the equivalent focal length calculated according to [
31] was two times larger than one calculated by means of a parabolic approximation.
It is known that the laser and thermal characteristics of Yb:YAG are strongly dependent on temperature. Therefore, in the case of a large temperature gradient, the accounting for this dependence becomes necessary, especially for the electronic contribution to the wavefront distortion. The influence of temperature-dependent parameters’ distribution on the determination of the OPD profile was analyzed. The thermal and electronic contributions introduced by the active element were calculated for the case of distributed medium parameters and for the case of the parameters taken at the minimal measured temperature of 165 K. The thermal contributions calculated in these ways were slightly different. The maximal difference between the calculated profiles was 57.5 nm, which was 2 orders lower than the amplitude of thermal contribution to the OPD profile. For the electronic contribution, the neglect of temperature dependence resulted in a more significant error. The difference in the calculated profiles increased from the center to the edges to 5.6 nm, which corresponded to a decrease in amplitude by 20% without taking into account the temperature dependence.
The difference between the thermal contribution to the OPD profile for radial and tangential polarizations, calculated from formula (1)
, increased from the center of the active element to the edge, from 0 nm to 150 nm, which was an order of magnitude smaller than the amplitude of the thermal contribution. The maximal error in determining the difference between the OPD for different polarizations was 10.5 nm. In the Yb:YAG case, for the radial polarization, the photoelastic effect contributed to the OPD with the same sign as other mechanisms, forming a positive lens; for the tangential component, it had a different sign. The difference in the OPD for the radial and tangential polarizations led to the radiation depolarization effect [
8,
9]. A thermally induced stress led to the birefringence due to the photoelastic effect. Radiation, passing through various regions of the active element, was decomposed into radial and tangential polarizations. The directions of the vectors, along which the radiation was decomposed, the phase difference and, consequently, the polarization of the radiation that passed through the active element, had spatial distributions. The depolarization effect is especially detrimental in cases where polarization-sensitive elements are used in laser systems. To calculate the radiation depolarization, the active element was considered as a set of wave plates, the action of which can be described by Jones matrices [
32] (p. 447).
The depolarization of radiation was calculated from the difference in the OPD for the radial and tangential components and measured experimentally by the setup shown in
Figure 1. The fraction of radiation power with a polarization state orthogonal to the initial state was determined (hereinafter, the fraction of depolarized radiation). For the measurement, incident radiation was horizontally polarized, and a polarizing beamsplitter was introduced into the radiation passed through the active element. The beamsplitter with a transmitted beam extinction ratio of 1000:1 and reflected beam extinction ratio of 100:1 was used. The power of the beams after the beamsplitter was measured. The fraction of depolarized radiation was determined as
, where
is the power of the horizontal and vertical polarizations. The depolarization of radiation caused by constant stresses was measured in the cooled active element without pumping radiation. The depolarization in the absence of pumping was no more than a half percent. In
Figure 6 are the calculated and measured maps of the fraction of depolarized radiation.
The maximal value of the fraction of depolarized radiation calculated from the OPD profiles was 18.3%. For the experimentally measured radiation depolarization profile, the maximal value was 17.9%. Thus, the results were in a good agreement. A slight difference in the profile shapes may have been caused by the fact that the calculation did not take into account the angle between the direction of propagation of the scanning radiation and the normal to the surface of the active element. In the region limited by the pump beam diameter of 2 mm, the maximum value of the experimentally measured fraction of depolarized radiation was 8.5%. This confirmed that it is necessary to take into account the depolarization of radiation and compensate for this effect in a 16-pass laser amplifier.