In large-aperture, high-energy laser systems, such as those used in inertial confinement fusion (ICF)-class laser systems, potassium dihydrogen phosphate (KDP) and its deuterated analog (DKDP) are commonly used owing to their ability to grow in large-size crystal boules and because of their high nonlinear coefficients^[1^]. As a result, KDP and DKDP crystal plates are utilized for second- and third-harmonic frequency generation (THG) and polarization control. As high-fluence laser light propagates through large-aperture (greater than 300 mm) crystal plates, spontaneously scattered light from the dominant A₁ Raman mode, associated with the totally symmetric ‘breathing mode’ of the PO₄ group, can experience high gain across the transverse dimensions of the plates. This stimulated Raman scattering process involves all photon paths, or propagation angles, confined within the plates via internal reflections from their surfaces, which, in the presence of the large-aperture laser beam propagating throughout the optic, introduces long amplification paths and, thus, high gain^[2^]. For a specific plate size and shape, the gain of this transverse stimulated Raman scattering (TSRS) is a function of the crystal cut orientation, which governs the 3D distribution of the spontaneous Raman cross-section and the laser parameters, including the wavelength and polarization state. As the spontaneous Raman cross-section rapidly increases with decreasing laser wavelength, TSRS in ICF-class systems becomes of concern in the ultraviolet section of the laser system, starting with the crystal plate converting the laser energy to the third harmonic^[3^]. Due to the splitting of the A₁ mode with the introduction of deuterium in the crystal, the peak Raman cross-section is lower in DKDP compared to that in KDP. As a result, DKDP is used for THG in the more powerful laser systems, such as the National Ignition Facility (NIF) and Laser Mégajoule (LMJ)^[4^,5^]. Furthermore, limited experiments exploring the use of DKDP for polarization control at the third harmonic showed prohibitively high TSRS gain and damaged optics and surrounding holders. On the other hand, KDP is successfully employed for harmonic generation and polarization control in lower-output-power laser systems, such as the OMEGA Laser System^[6^].

TSRS amplification is a coherent process that grows exponentially and is distributed nonuniformly in the crystal and at the crystal surfaces. Due to this exponential growth, the fluence level of the TSRS signal reaching the crystal perimeter can vary from harmless levels to ones that cause damage to the crystal or ablation of mechanical components in the vicinity. Methods to manage the TSRS effect have been developed for the THG plate by beveling the edges of the crystal to avoid back-reflections, and thus additional amplification, during the laser pulse^[3^]. More recently, innovative ideas for TSRS management have emerged and include the introduction of a laser-induced damage array composed of numerous pinpoints inside the plate to form a barrier to prevent signal amplification^[7^]. A second approach introduces a polarization element at the fundamental frequency, which divides each beam into several sub-beams whose polarization directions are orthogonal as they propagate through the two subsequent frequency triplers^[8^]. However, these proposed solutions can introduce significant unintended consequences, such as beam-intensity modulations.

Detailed studies to understand such TSRS and its distribution for a specific crystal application configuration have been very limited^[9^]. To understand the performance limits of current-generation systems and in anticipation of the next-generation ICF lasers for basic research and possible inertial fusion energy production, the understanding and the management of the growth and spatial distribution of TSRS energy are of fundamental importance. In this work, a modeling approach is discussed and used to simulate TSRS generation for several practical configurations, such as harmonic generation and polarization control. The results of this study provide information on the angular dependence of TSRS as a function of the crystal cut orientation and pump polarization. Based on these results, we explore approaches to minimize the TSRS gain and growth in current-generation designs and configurations.

Stimulated Raman scattering is a process that involves amplification of the spontaneous Raman signal of a particular vibration mode, typically the mode that produces the highest intensity. The spatial distribution of the spontaneous Raman energy is a function of the incident and emitted electric-field polarizations and the directions of incident and emitted waves and quantified by the scattering cross-section, σ_mol, of one crystal molecule in the following form:

$$\begin{align}{\sigma}_{\mathrm{mol}}={A}_{\mathrm{mol}}\cdot {\left|\boldsymbol{E}_{\mathrm{P}}\cdot \boldsymbol{R}\cdot \boldsymbol{E}_{{R}}\right|}^2,\end{align}$$  (1)

Raman scattering in KDP and DKDP has been extensively studied, but the Raman scattering tensor of the dominant A₁ vibration mode was only recently accurately measured following more than four decades of effort. The complexity of the measurement was due to the birefringence in KDP that causes signal artifacts, as detailed in previous work^[10^]. Efforts also focused on the measurement of the TSRS gain in KDP and KDP, including a quantitative measurement of the peak value of the Raman cross-section at different wavelengths^[11^,12^]. The normalized Raman tensor for this mode in the crystal frame is as follows:

$$\begin{align}\boldsymbol{R}\left(\mathrm{A}_1\right)=\left(\begin{array}{ccc}1& 0& 0\\ {}0& 1& 0\\ {}0& 0& B\end{array}\right),\end{align}$$  (2)

KDP and DKDP are birefringent materials in which the polarizations of orthogonal modes are along the ordinary (o) and extraordinary (e) directions. The o-polarization electric-field vector is normal to both the k vector and the crystal optic axis, and the e-polarization electric-field vector is normal to both the k vector and the o-polarization vector (Figure 1). The index of refraction for the e-polarization depends on the direction of propagation, so the o- and e-components travel at different velocities. During the pump laser pulse, Raman photons are generated into these two polarization modes and amplified by separate gain coefficients (to be discussed later). Because of the velocity difference, the two pump polarization components quickly lose their phase relationship as they travel from the point at which the pump ray enters the crystal to the point in the crystal that generates the spontaneous Raman scattering. Thus, the Raman scattering is effectively separated into the two incoherent (o and e) components as determined by the pump beam orientation (see Section 6).

The pump o- and e-polarization unit vectors are defined as the following vector products:

$$\begin{align}{\boldsymbol{o}}_{\mathrm{P}}={\boldsymbol{k}}_{\mathrm{P}}\times \overrightarrow{\mathrm{OA}},\quad{\boldsymbol{e}}_{\mathrm{P}}={\boldsymbol{k}}_{\mathrm{P}}\times {\boldsymbol{o}}_{\mathrm{P}},\end{align}$$  (3)

In a similar manner, the Raman scattering (or TSRS) generated and propagating in a specific direction is also broken into o- and e-components. As with the pump light, the Raman o- and e-components generated have no fixed phase relationship as they propagate from the source point to the point at which they exit the crystal:

$$\begin{align}{\boldsymbol{o}}_{\mathrm{R}}={\boldsymbol{k}}_{\mathrm{R}}\times \overrightarrow{\mathrm{OA}},\quad{\boldsymbol{e}}_{\mathrm{R}}={\boldsymbol{k}}_{\mathrm{R}}\times {\boldsymbol{o}}_{\mathrm{R}}.\end{align}$$  (4)

The generated o- and e-components of the spontaneous Raman (or TSRS) signal at any point receive contributions from both the o- and e-components of the pump excitation. These values are determined by the Raman cross-section in the direction of the ray for each Raman and pump polarization. However, as discussed above, both the pump polarization and Raman polarization undergo polarization evolution due to the birefringence of the material. As a result, the o- and e-components of the integrated (total) Raman cross-section generated within a unit volume can be expressed as two separate functions.

For the o-Raman polarization:

$$\begin{align}{\sigma}_{\mathrm{o}}\left({\boldsymbol{k}}_{\mathrm{R}}\right)&={A}_{\mathrm{vol}}\cdot \left[{\left|\left({\boldsymbol{E}}_{\mathrm{P}}\cdot {\boldsymbol{o}}_{\mathrm{P}}\right)\cdot \left({\boldsymbol{o}}_{\mathrm{P}}\cdot \boldsymbol{R}\cdot {\boldsymbol{o}}_{\mathrm{R}}\right)\right|}^2\right.\nonumber\\ &\quad\left.+{\left|\left(\boldsymbol{E}_{\mathrm{P}}\cdot {\boldsymbol{e}}_{\mathrm{P}}\right)\cdot \left({\boldsymbol{e}}_{\mathrm{P}}\cdot \boldsymbol{R}\cdot {\boldsymbol{o}}_{\mathrm{R}}\right)\right|}^2\right]={A}_{\mathrm{vol}}\cdot {Q}_{\mathrm{o}}\left({\boldsymbol{k}}_{\mathrm{R}}\right).\end{align}$$  (5)

For the e-Raman polarization:

$$\begin{align}{\sigma}_{\mathrm{e}}\left({\boldsymbol{k}}_{\mathrm{R}}\right)&={A}_{\mathrm{vol}}\cdot \left[{\left|\left(\boldsymbol{E}_{\mathrm{P}}\cdot {\boldsymbol{o}}_{\mathrm{P}}\right)\cdot \left({\boldsymbol{o}}_{\mathrm{P}}\cdot \boldsymbol{R}\cdot {\boldsymbol{e}}_{\mathrm{R}}\right)\right|}^2\right.\nonumber\\ &\left.\quad+{\left|\left(\boldsymbol{E}_{\mathrm{P}}\cdot {\boldsymbol{e}}_{\mathrm{P}}\right)\cdot \left({\boldsymbol{e}}_{\mathrm{P}}\cdot \boldsymbol{R}\cdot {\boldsymbol{e}}_{\mathrm{R}}\right)\right|}^2 \right]={A}_{\mathrm{vol}}\cdot {Q}_{\mathrm{e}}\left({\boldsymbol{k}}_{\mathrm{R}}\right).\end{align}$$  (6)

To verify these formulae, cross-section values for the in-plane directions were experimentally measured with spherical crystal samples^[13^]. This measurement was carried out in the x–z plane (which means the Raman k vector is in this plane). Separate measurements were made to quantify the vertical (z) and the horizontal (in-plane) polarization components. Agreement between experiments and theory provides the basis for quantitative evaluation of the Raman scattering cross-section and TSRS gain in three dimensions.

The azimuthal angle ( $\phi$ ) and the internal angle (α) are used to define a Raman propagation direction over 4π solid angles, as depicted in Figure 2. In the latter part of this paper, we are only concerned with rays that undergo total internal reflections (TIRs) inside the plate (thus for α < 42°) because rays outside this range will lose energy after experiencing multiple reflections between the top and bottom surfaces of the crystal plate. Ray paths that do not support TIR are not considered as contributors to TSRS.

Figures 3 and 4 show the 3D normalized cross-section functions for the full 4π solid angle. The strong dependence for both polarizations provides the basis for the later discussion of optimization by choosing the OA tilt angle to reduce TSRS.

Let us consider a Raman ray that has been generated in a volume element dV during a time interval dt. This ray travels in the direction of k_R as a cone for a distance l inside the crystal, and has a solid angle element of dΩ. The initial spontaneous Raman emission ( $\mathrm{d}{E}_{\mathrm{TSRS}\;\mathrm{o},\mathrm{e}}^{\mathrm{o}}$ ) has two components, o and e, that are related to the normalized cross-sections for o and e as discussed in the previous section:

$$\begin{align}\mathrm{d}{E}_{\mathrm{TSRS}\;\mathrm{o},\mathrm{e}}^{\mathrm{o}}\left({\boldsymbol{k}}_{\mathrm{R}}\right)={A}_{\mathrm{vol}}\cdot {Q}_{\mathrm{o},\mathrm{e}}\left({\boldsymbol{k}}_{\mathrm{R}}\right)\cdot \mathrm{d}t\cdot {I}_{\mathrm{pump}}\cdot \mathrm{d}V\cdot \mathrm{d}\Omega,\end{align}$$  (7)

After this ray propagates by a length l, its Raman radiation energies ( $\mathrm{d}{E}_{\mathrm{TSRS}\;\mathrm{o},\mathrm{e}}^{\mathrm{o}}$ ) are the corresponding initial emitted energies (Equation (7)) multiplied by the exponential gain factor for each individual polarization, which is governed by the propagation length, the laser intensity I_pump and the Raman gain coefficient:

$$\begin{align}\mathrm{d}{E}_{\mathrm{TSRS}\;\mathrm{o},\mathrm{e}}^l\left({\boldsymbol{k}}_{\mathrm{R}}\right)&=\mathrm{d}{E}_{\mathrm{TSRS}\;\mathrm{o},\mathrm{e}}^{\mathrm{o}}\left({\boldsymbol{k}}_{\mathrm{R}}\right)\nonumber\\[2pt] &\quad \cdot \exp \left[\frac{8\pi cM}{h{\omega}_{\mathrm{R}}^3{n}_{\mathrm{o},\mathrm{e}}^2\Delta \nu }{I}_{\mathrm{pump}}{A}_{\mathrm{mol}}{Q}_{\mathrm{o},\mathrm{e}}\left({\boldsymbol{k}}_{\mathrm{R}}\right)l\right]\nonumber\\[2pt] &=\mathrm{d}{E}_{\mathrm{TSRS}\;\mathrm{o},\mathrm{e}}^{\mathrm{o}}\left({\boldsymbol{k}}_{\mathrm{R}}\right)\cdot \exp \left[{I}_{\mathrm{pump}}{gQ}_{\mathrm{o},\mathrm{e}}\left({\boldsymbol{k}}_{\mathrm{R}}\right)l\right].\end{align}$$  (8)

Here, g is the generally accepted notation for the gain coefficient with the dimension of cm/GW:

$$\begin{align}g=\frac{8\pi cM}{h{\omega}_{\mathrm{R}}^3{n}_{\mathrm{o},\mathrm{e}}^2\Delta \nu }{A}_{\mathrm{mol}},\end{align}$$  (9)

We now combine the contributions of the normalized cross-section functions and the propagation length to define G_o,e = Q_o,el as the gain factor. This factor contains all the dependence of vectors in Figure 1 and includes the contribution of the path length.

In a ray-tracing calculation, we consider N rays generated in a crystal unit volume dV and during a time interval dτ that is a fraction of the pump pulse duration τ. The generated rays are propagating in a 4π solid angle with the initial Raman energy (generated within ΔV and Δτ) in the o- or e-polarization, expressed as follows:

$$\begin{align}{E}_{\mathrm{TSRS}\;\mathrm{o},\mathrm{e}}^{\mathrm{o}}\left({\boldsymbol{k}}_{\mathrm{R}}\right)=4\pi \Delta V\cdot \Delta \tau \cdot {I}_{\mathrm{pump}}\cdot {A}_{\mathrm{vol}}\cdot {Q}_{\mathrm{o},\mathrm{e}}.\end{align}$$  (10)

For conservation of energy, we assume that the energy of each ray propagating along a direction k_R is ${E}_{\mathrm{TSRS}\;\mathrm{o},\mathrm{e}}^{\mathrm{o}}\left({\boldsymbol{k}}_{\mathrm{R}}\right)/N.$ The total Raman energy generated by a pulse of length τ is the sum (integral) over all volume elements, solid angles and initial time when the Raman seed signal is generated for both o- and e-ray spontaneous Raman energy (Equation (7)) multiplied by the gain factor exp(I_pump ⋅ g ⋅ G_o,e). This sum should converge to be independent of N when the sample ray density is sufficiently large. In the numerical calculations in Section 5 over 30,000 rays were used for each data point and doubling the number of rays changed the results by less than 1%.

As an illustration to help understand the distribution of generated Raman scattering at the edge of the crystal, let us consider a source point in the middle of a square plate. The normalized gain factor is calculated for all TIR directions (both α and $\phi$ ) from this source point to the edges of the plate. Figure 5 shows this normalized gain factor for the case of a square plate with the OA on the y–z plane with tilt angle of 60° and pump polarization along the diagonal direction of the plate ( $\phi$ = 45°). The location of the maximum value of G_o is at $\phi$ = 90° and $\phi$ = 270° (indicated with arrows), while G_e presents a much lower value. We notice that the range of angles with a significant gain factor is restricted to a rather narrow area (solid angle) near the limit of the TIR condition. This pattern depends on the location of the source point. For each different configuration (i.e., crystal geometry, OA cut and pump laser polarization), the distribution patterns of the gain factors can be determined and used to calculate the accumulated TSRS energy around the edges of the crystal.

Additional considerations for modeling the gain factor assuming square optics are as follows:

1. (i)the path length of a specific ray from one side to the other side is D/cos(α)cos ( $\phi$ ) (Figure 5);
2. (ii)the gain of rays confined by TIRs is the averaged values between the two principal directions (toward and away from the input surface of the plate);
3. (iii)combining these two effects, the maximum path length of all rays is approximately 2D.

This is shown in Figures 6(a) and 6(b).

The ray-tracing model utilizes a volume grid (element) for the bulk of the crystal and a surface grid for the side surfaces of the crystal. The amount of Raman signal (energy) originated from each ray (starting from a volume element during a time interval to eventually arrive at the surface element) is proportional to the pump intensity, the volume element size, the time interval, the cross-section for this ray direction and the solid angle element associated with the surface element extending from the source point. For a ray that is confined via TIR between the top and bottom surfaces, this solid angle is the surface element area projected to the normal direction of the ray, divided by the total path length. Rays from different locations in the crystal plate arrive at the same sampling point at a side surface. This process is depicted in Figure 6(b) where each arrow represents many rays that have different TIR internal angles, as shown in Figure 6(a).

For each source point and a given internal angle, a pair of rays (o and e) is generated at one time interval during the pump pulse. Each is amplified by its corresponding gain until they arrive at the next TIR point on the top or bottom surface. At this point the o- and e-rays are projected to the s and p polarization directions (related to the beam incidence on the surface). Following a TIR, the reflected rays become s^′ and p^′, where s^′ remains the same as s but p^′ changes direction from p by 2 $\alpha$ . From s^′ and p^′ we reconstitute o^′ and e^′ intensities and continue the ray propagation using the gains for the new direction. The total TSRS reaching a sampling area element is the sum of all rays reaching this location from all source points and angles and includes consideration of the time at which the initial spontaneous Raman photon is generated during the pulse length. Rays initiated during the pulse are amplified for the part of pulse between its starting moment until the end of the pulse or until it hits the sampling area, whichever is earlier.

Within this modeling, side surfaces are considered but reflections at the edge surfaces are not. When a ray reaches the sampling point on the edge, we add both o- and e-intensities. We assume a 40 cm × 40 cm × 1 cm plate and use a 3-ns pump pulse, as used by Dixit et al.^[9^] . This interval is the time required for light to travel across the plate parallel to the surface. There are, however, rays that take longer than 3 ns to reach an edge. For example, a ray that zigzags from one corner to the opposite corner takes approximately 4 ns to propagate. When the 3-ns pump pulse is over, those rays that have not reached any of the four edge surfaces continue until they reach one of those edges, but there is no further amplification. As we showed in the previous section, the rays with a maximum internal angle, and consequently a longer path length, contribute to TSRS most significantly. We first use this model to simulate the most detailed results available to date on TSRS presented by Dixit et al.^[9^]. Although the exact locations of the measurements are not included in this manuscript, we use the gains and cross-sections that would fit their data, as follows.

KDP: cross-section A_vol = 3.47 × 10^–7 cm^–1 sr^–1, gain g = 0.347 cm/GW

DKDP (70%): cross-section A_vol = 2.9 × 10^–7 cm^–1 sr^–1, gain g = 0.203 cm/GW

These values are then used in our modeling to generate detailed data for various situations.

Let us first consider the third-harmonic crystal configuration used on the NIF laser, which has the 3ω polarization oriented in the y direction (see Figure 2). The incident 1ω and 2ω rays at the front surface are converted into 3ω as they propagate through the plate^[14^] and the 3ω light only reaches its full power at the exiting surface of the crystal plate. We incorporate this 3ω distribution in our modeling to calculate the TSRS distribution in the side surfaces of the crystal. In the plots shown in Figure 7, the TSRS values are given in arbitrary units since the focus of this study is to determine the comparative TSRS values between the existing THG and proposed polarization rotating devices. Also, because the details of the NIF measurements presented in Ref. [9], such as how and where the measurement was obtained (there is a great variation of intensity distribution at the edges, as shown by the results), are not available and the model results are based on certain assumptions discussed earlier, we think it is more prudent not to use absolute values.

Figure 7 shows the TSRS distributions for the THG configuration for different pump intensity levels. The inset shows a schematic depiction of the crystal configuration, where the pump polarization is along the y-axis and the OA orientation is at an angle of θ = 59° along the z–y plane. The results show that the TSRS intensity varies within the crystal side surfaces and reaches maximum intensity in the middle regions of the x side along the y-axis. It is thus expected that significant TSRS gain will be first generated along this direction.

Next, we consider the TSRS for a wave-plate configuration suitable for 90° rotation of polarization at 3ω. In this configuration (Figure 8), the crystal OA is in the plane that cuts through the diagonal line of the plate. With the OA at θ = 10° (Figure 8(a)), two corners receive several orders of magnitude higher TSRS compared to that observed for the THG configuration. With the OA at θ = 90° (Figure 8(b)), the variation of TSRS along each side is not significant, but the average value is even higher. It is appreciated that the spatial distribution along the edge surfaces parallel to both the x-axis and the y-axis is symmetric and that the maximum intensity is observed along the diagonal direction.

Since the TSRS distribution depends on the direction of the OA, we explored the change of its peak value as a function of the OA tilt angle, as shown in Figure 9. These results show that by using a θ = 90° (OA tilt angle), the TSRS can be reduced by more than a factor of 10. However, the TSRS level is still several orders of magnitude higher than the tripler configuration (see Figure 7(b)).

An alternative configuration for the wave plate can be considered by assuming the beam polarization is at 45° and the OA is in the y–z plane, as depicted by the inset of Figure 10, where the OA tilt angle is used as a parameter for optimization. As can be seen from the results shown in Figure 10, this configuration makes it possible to have a wave plate that produces maximum TSRS at 2 GW/cm², similar to that of the tripler.

These results can be viewed in two different ways. Firstly, we can examine the maximum TSRS fluence for a fixed pump level (of 2 GW/cm²) and consider that of the tripler crystal as the reference value, since the tripler is inevitably going to be part of the laser system (see Table 1). We can also examine the maximum possible pump level for each configuration that will facilitate the same maximum TSRS, which in this case is considered to be that of the tripler for a pump intensity of 2 GW/cm². Arguably, these results provide a better understanding on the TSRS-induced limitations of current systems and guide potential paths forward.

The TSRS can be treated as an incoherent superposition of the two polarization components, both temporarily and spatially. In reality, Raman emission and amplification have finite spatial and temporal coherences, which have been discussed extensively by Raymer et al.^[15^]. The degree of spatial coherence and intensity modulation depend on the geometry of the medium where the Raman scattering is generated and amplified. A square meter-scale crystalline plate, used in large aperture solid-state laser systems, generally has a low degree of Raman spatial coherence at the plate’s edge. On the other hand, a long cylindrical medium, such as laser propagation over long air paths, tends to have a higher degree of Raman spatial coherence at its end. Based on Raymer’s theory, we conclude that the TSRS in the ICF plates has weak coherent effects, both spatially and temporally. The time-integrated speckle intensity contrast is calculated to be in the range of a few percent or lower. Thus, this additional factor does not significantly affect our conclusions.

The aim of this work is to demonstrate a method to quantitatively evaluate TSRS in a configuration relevant to large-aperture laser systems. The specific case examples presented herein were focused on simulating the most detailed results available to date on TSRS generation in large-aperture laser systems presented by Dixit et al.^[9^]. For this reason, we assumed a 40 cm × 40 cm × 1 cm plate and used a 3-ns flat in time pump pulse with a uniform spatial beam profile. We are also seeking the TSRS conditions at the onset of strong nonlinear amplification. Thus, we do not consider the depletion of the pump or filamentation of the pump. The model can be readily modified to include temporal and spatial variations of the pump. However, depletion of the pump is a more difficult problem because it breaks the independent evolution of Raman rays and makes parallel computation invalid.

In this work we are also not considering amplification of the reflected TSRS back to the crystal at the edges. This can be a major technical challenge and can be mitigated using various methods that have already been utilized, such as edge cladding, beveling, antireflection coatings or special cuts of the crystal to direct reflection of the ‘worst rays’ in directions that do not support additional gain. Another solution to be explored is to apply an index-matched absorptive coating on plate edges to further reduce TSRS propagation. Such index-adjustable absorptive paint has been developed and is available^[16^,17^]. Furthermore, controlling the pulse shape can be an additional tool. Ultimately, there are a number of ways to design wave plates that suppress TSRS gain that may also require control of the incoming polarization state with respect to the geometric orientation of the plate (such as the rotation of pump polarization by 45° discussed above). The general method presented in this work can be applied in a variety of designs and specific excitation parameters.

Finally, it must be noted that the laser-induced damage on the material is based on the total local fluence with contributions from all wavelengths. Given the small shift in energy between the pump beam and the generated TSRS (associated with a shift in wavelength of about 12 nm), the damage probability should be considered as the sum of TSRS and pump fluences. In addition, unconverted second-harmonic light will also contribute, to a lesser extent, to the damage initiation process^[18^]. Once damage is initiated and plasma is formed, the remaining unconverted fundamental and second-harmonic beams will strongly contribute to the energy deposition process due to absorption by the formed plasma, and thus the final size of the damage site. One should also consider that the laser-ablation threshold of metals that may be used in optic mounts is very low (of the order of 250 mJ/cm²). As a result, even a small amount of TSRS, which does not pose a direct concern for the optic that escapes from the side surfaces of the plate, can initiate contamination of the optic by particles generated from the ablation of metal holders. These particles can subsequently initiate damage (and downstream beam modulation) that can also have a detrimental effect on the lifetime of the optic and the laser system in general. Thus, it is imperative that the spatial distribution of the TSRS for a particular system geometry is considered in the design of the mounts and the selection of materials used for mounting.

TSRS in laser crystal optics grows exponentially with the pump fluence and can cause damage to laser components. An incoherent propagation model for TSRS that allows the determination of the Raman cross-section and gain in any direction in three dimensions is constructed. Ray-tracing calculations utilizing these cross-sections and gain values lead to the determination of TSRS generated in a crystal during a pulse as well as the distribution of TSRS fluence at the edges of the crystal plate.

By varying crystal configuration parameters, we can minimize the TSRS in distributed polarization ratios to be used in current high-power laser systems.