We consider the propagation of powerful laser radiation in an anisotropic medium with natural birefringence and cubic nonlinearity. By the example of a quarter-wave plate, we show theoretically and experimentally that, under the simultaneous influence of linear birefringence and nonlinearity, the accuracy of polarization transformation decreases in proportion to squared В-integral.
© 2013 OSA
Powerful laser systems generating petawatt power femtosecond pulses having intensities of 1020 W/cm2 and higher open a way to the rapidly advancing trend in laser physics – creation of compact accelerators of high-energy particles, such as electrons, protons (ions) and neutrons. The principle of operation of these accelerators is based on the interaction of a powerful laser pulse with a thin solid target (e.g., aluminum or gold foil). During several femtoseconds the target is ionized so that free electrons may be accelerated under the action of ponderomotive force of the laser pulse, giving rise to strong charge separation into electrons and heavier, hence, immobile ions. The resulting ambipolar electric field leads to acceleration of protons in the target plasma in different directions.
Of major interest are mechanisms of proton acceleration as proton beams with energies from several tens to hundreds of MeV may be used for proton magnetic resonance imaging of small-scale objects , hadron therapy in oncology , and fast ignition of a thermonuclear target . According to the works by [4–6] laser radiation polarization is an important parameter affecting the process of proton acceleration. Investigation of the intense laser pulse interaction with a thin foil demonstrated that proton acceleration efficiency with the use of circularly polarized radiation is higher than with linearly polarized radiation . In the case of circular polarization, ambipolar electrostatic field has positive projection onto the direction of laser pulse propagation; hence, the protons are accelerated only forward. In the case of linear polarization, ponderomotive force has an oscillatory longitudinal component providing both forward and backward proton acceleration. Moreover, it was found that the use of circularly polarized radiation allows more efficient proton beam collimation, as compared to linearly polarized radiation .
Thus, it is necessary to produce circularly polarized radiation with the intensity of about 1020 W/cm2. Consider the most broadly used method – the transformation of linear polarization to a circular one by means of a quarter-wave plate. In this case, it is important that polarization transformation should be accomplished in the region of intense radiation. Indeed, the CPA (Chirp Pulse Amplification) as well as OPCPA (Optical Parametric Chirped-Pulse Amplification) technology used for laser radiation amplification requires linear polarization, as diffraction gratings that are inherent components of the system operate with maximum efficiency with polarization of this type. Consequently, it is interesting to investigate the simultaneous influence of the properties of the wave plate material (natural birefringence and nonlinearity) on laser radiation polarization.
We will be interested in wave plates cut from crystalline quartz, possessing both quadratic and cubic nonlinearity. In this work we consider plates for polarization rather than frequency transformation. The wave plates position relative to the direction of radiation propagation is such that the phase synchronism condition responsible for the effects related to quadratic nonlinearity is not fulfilled a fortiori. Therefore, we will neglect nonlinearity of this type.
Effect of polarization ellipse rotation caused by cubic nonlinearity is well-known [8–12]. In paper  authors investigated the isotropic cubic nonlinearity influence on accuracy of polarization transformation. Reduction of polarization transformation accuracy for the В-integral more than unity was shown. In this paper we considered the anisotropic cubic nonlinearity because crystalline quartz is positive uniaxial crystal possessing trigonal crystalline structure. Therefore, we had to take into account not only diagonal elements of fourth-order nonlinear susceptibility tensor χijkl but also the nondiagonal ones. This correction was very important because we proofed experimentally that the effect of depolarization growing is stronger than we calculated using model of isotropic nonlinearity. In other words, according to calculation within model of isotropic nonlinearity the 1% depolarization degree achieved for higher intensity (I ~4 TW/cm2) then in case of the anisotropic nonlinearity (I ~700 GW/cm2).
In the second and third sections of the paper we will consider nonlinear optical properties of crystalline quartz and propose a system of differential equations describing laser radiation propagation in an anisotropic medium with cubic nonlinearity. Results of theoretical and experimental investigation of the influence of cubic nonlinearity on accuracy of polarization transformation by means of a quarter-wave plate will be presented in Sec. 4, where it will be shown that the degree of depolarization (energy fraction in the radiation with polarization orthogonal to the polarization in the absence of nonlinearity) and the ellipticity of output radiation polarization at the output of the quarter-wave plate is proportional to squared intensity (В-integral).
2. Nonlinear optical properties of crystalline quartz
Crystalline quartz belongs to the group of positive uniaxial crystals. At room temperature, quartz consists of α-modification crystals possessing trigonal crystalline structure . The permittivity and gyration tensors for a structure of this type are written in the formEq. (2) permits characterizing the principal properties of quartz. For example, if a quartz plate is cut parallel to the optical axis, then the birefringent properties of quartz are shown due to the difference between ε11 and ε22. If the plate is cut normally to the optical axis, then quartz possesses gyrotropic properties, which is used in manufacturing 90° polarization rotators. In the present work we will be interested in plates cut parallel to the optical axis of the crystal – Oz.
To reduce the influence of refractive index dispersion on radiation polarization, zero-order plates are used. In λ/4, λ/2, λ zero-order plates the phase difference between the extraordinary and ordinary waves is π/2, π, and 2π, respectively. However, the thickness of the λ plate would then be L = 90 μm (at λ = 795 nm), which is hard to achieve in practice. Hence, two multiple-order wave plates are used, the difference in thicknesses of which produces path difference of about 1/4 or 1/2 of wavelength, depending on the plate type. Optical axes of these multiple-order plates make an angle of 90° with each other (Fig. 1 ).This is the type of wave plate of interest to us.
The influence of cubic nonlinearity in a wave plate of crystalline quartz is specified by nonzero components of the fourth-order nonlinear susceptibility tensor χijkl:Eq. (3) it follows that nonlinear properties of the crystal also depend on the way a wave plate is cut relative to the Ох and Оу axes. Note that there is no symmetry in tensor χ with respect to subscripts x and y. For instance, the tensor component with subscripts yyyz has no pair component with subscript xxxz. Consequently, if the plate is cut along the Оу axis, components with indices y and z play an important role in the system of differential equations. Quite the contrary, if the plate is cut along the Ox axis, the components with indices х and z are of primary importance. It is worth mentioning that the choice of crystal cutting along the Оу and Ox axes does not affect the birefringent properties of the plate, which complicates identification of the direction cutting.
As a measure of nonlinearity researchers usually use В-integral – phase incursion related to the nonlinear part of refractive index n. For example, in an isotropic medium with cubic nonlinearity, n linearly depends on intensity n(I) = n0 + γNLI, where n0 is a linear refractive index and γNL is the characteristic of the nonlinear medium. In this case, for a medium with length L, the magnitude of В-integral is defined by
The magnitude of В for an anisotropic medium depends on polarization and direction of radiation propagation. Hence, in the present work we will determine В for radiation with linear polarization coinciding with the polarization of an ordinary or extraordinary wave of crystalline quartz. The magnitudes of γNL corresponding to eigen modes of the crystal, according to , differ less than by 4%; so we set the largest γNL = 3.2⋅10−7cm2/GW.
3. System of differential equations
With the crystal cut parallel to the Ох axis, the Оу axis will be taken to be the direction of radiation propagation and Ох and Оz will be directed along the optical axes in the first and second parts of a composite wave plate (Fig. 1). Then, the impact of cubic nonlinearity will be described by the χijkl components, where i, j, k, l are х and z. Whereas with the crystal cut parallel to the Oy axis, formal substitution of x for y should be done in the equations taking into account the new terms containing components of permittivity tensor Eq. (3) with indices i, j, k, l are y and z.
The change of the parameters of radiation is described through vector projections onto the principal directions of the plate for the two cases of plate cutting:16,17].
Our estimates show that the length of pulse broadening due to the difference in group velocities of the “blue” and “red” wavelengths that limit the radiation spectrum used in our experiment is much larger than the thickness of the considered wave plate. Consequently, neglecting dispersion and diffraction, a steady-state system of differential equations for electric field vector in the case of the radiation propagating in a nonlinear medium with a trigonal crystal lattice will be written asEq. (6а) and δ(x) = k(nz – ny)x in Eq. (6b) are the phase differences of the extraordinary and ordinary waves. System Eq. (6а) corresponds to the case of the plate cut parallel to the Ох axis, and Eq. (6b) parallel to the Оу axis. From  it follows that χzzzz = 1.89⋅10−14 esu, χyyyy = χxxxx = 1.82⋅10−14 esu. The magnitudes of χxxzz, χyyzz and χyyyz are unknown to the best of our knowledge.
We remind the readers that in the present work we consider the transformation of linear polarization to the circular one using a zero-order quarter-wave plate. Assume that the vector of the electric field incident on the plate makes an angle of 45° with the optical axes of the plate (Fig. 1). In this case projections of vector onto the plate axes are equal. Since χzzzz ≅□χyyyy = χxxxx, from previous remark and system Eq. (6) it follows that it is the terms containing nondiagonal elements of tensor χ that will determine the change of polarization.
The system Eq. (6а) or Eq. (6b) allows calculating energy redistribution between the Ех and Еу components as well as Еz at the output of the plate with preset δ and В. Such a change will be characterized by the degree of local (at each point of the cross-section) depolarization – the fraction in the output radiation Eout(B ≠ 0) of the intensity with polarization orthogonal to the polarization in the absence of nonlinearity Eout⊥(B = 0) (i.e., circular polarization but with a different direction of rotation):Eq. (6) demonstrated that the impact of cubic nonlinearity reduces accuracy of polarization transformation by means of wave plates. In other words, if the radiation polarization at the input of a quarter-wave plate is such that clockwise circular polarization must be obtained at the output, then the depolarization (energy fraction contained in the counterclockwise circular polarization) at the plate output is proportional to squared В-integral. Experimental verification of this result is presented in Section 4.
Calculation of Γ allows us finding out the ellipticity Σ of the polarization at the output of the wave plate
4. Experimental investigation of the influence of cubic nonlinearity on accuracy of transformation polarization by means of a quarter-wave plate
We measured radiation polarization ellipticity of a terawatt-power titanium-sapphire laser with central wavelength λ = 795 nm . In other words, we determined depolarization of the laser system specified by accuracy of obtaining linear polarization. The depolarization degree was calculated as a fraction of the energy contained in the radiation component with polarization corresponding to the small semi-axis of the polarization ellipse. After that, we experimentally investigated the dependence of the energy fraction of this component in the output radiation on the magnitude of В-integral.
4.1 Schematic of the experiment
The schematic of the experimental setup is shown in Fig. 2 (а). On reflection from mirror 1 the main portion of the radiation passed through Galilean telescope 2-3 with magnification 1:2.5. This telescope permits increasing intensity up to the calculated level (~700 GW/cm2) at which the influence of cubic nonlinearity on polarization transformation in quarter-wave plate 4 was calculated to be high enough. One more plate 6 placed in lower-intensity radiation reflected from glass wedge 5 was used to increase accuracy of depolarization measurements. Hence, cubic nonlinearity could not affect polarization transformation in plate 6. After that the principal radiation component with polarization corresponding to a large semi-axis of the polarization ellipse and the depolarized component with orthogonal polarization were separated in space by calcite wedge 7. The calcite wedge was adjusted so that its principal axes coincided with the principal axes of the polarization ellipse at В = 0. The image was transferred from the plane of wave plate 4 to the matrix of CCD-camera 9 by lens 8. Thanks to the additional wave plate 6 radiation polarization at the output of the experimental setup at В = 0 coincides with the initial polarization. In this case, two-dimensional intensity distributions of the principal (bright spot) and depolarized (spot at the level of CCD-camera noise) components of radiation were recorded on the CCD-camera (Fig. 2(b), 2(c)). The brightness of the spot corresponding to the depolarized component increased at В > 1 (Fig. 2(e)).
The transverse intensity distribution at the input had a Gaussian profile ~exp(−r/R0)2, where R0 = 0.24 cm. The radiation time distribution is described well by the function 1/ch2(t/t0), where t0 = 28 fs, which corresponds to the pulse duration of 75 fs FWHM. Pulse energy was measured within the 1-7 mJ range.
4.2 Experimental results
Results of the experiments are presented in Figs. 3(a) , 4(a) (the same experimental points are used in these plots). The local depolarization degree Γ and ellipticity Σ calculated by Eq. (9) during processing of two-dimensional intensity distributions depicted in Figs. 2(с), 2(е) were calculated in the central region of the beam. The fraction of energy contained in the depolarized radiation component as a function of radiation intensity without wave plates is plotted by blue dots. In this fashion we determined ellipticity of the initial polarization, i.e. its difference from the linear polarization. One can readily see that the level of cold/initial depolarization is 0.2%. Red triangles were used to plot analogous dependences after wave plates 4 and 6 were added to the scheme. In this case the dependence of depolarization on intensity I (В-integral) is nonlinear. The value of Γ(B = 0) equal to 0.1% is less than the analogous value in the absence of wave plates 4 and 6, which may be attributed to reduced polarization ellipticity due to the plates. The maximum intensity magnitude of 700 GW/cm2 achieved in experiment corresponds to Г = 1.1% and Σ = −0.87.
For comparison of the obtained experimental and theoretical data it is necessary to determine the region of possible values of components of nonlinear susceptibility tensor χ unknown for crystalline quartz. Let us take I = 700 GW/cm2 and plot a one-dimensional Γ(χxxzz) and a two-dimensional Γ(χyyzz, χyyyz) functions for each type of wave plate cutting, respectively (Figs. 3(b), 4(b)). The resulting dependences will be used to determine values of the sought components of tensor χ, for which Г(I = 700 GW/cm2) = 1.1%. Further, we will search among the selected values of χ components the ones at which the theoretical and experimental curves for Г(I) coincide at all points (Fig. 3(a), 4(a)). Note, we assume that the both parts (multiple-order plates) of the quarter-wave plates used in experiments have the identical direction cutting.
Two possible values of χ: χxxzz = 0.76 χxxxx and χxxzz = −0.28 χxxxx were found for the plate cut parallel to the Ох axis (Fig. 3(b)). The curves in Fig. 3(b) are plotted for different values of radiation polarization ellipticity at the input, i.e., for the depolarization Г(В = 0) = 0%, 0.01% and 0.1%. Thus, we have found that the uncertainty of determining Г(В = 0) in the given range of values leads to 8% error in calculations of χxxzz. Note that the obtained values of χxxzz demonstrate a good agreement of experimental and theoretical data for the clockwise polarized radiation (black solid and dashed curves in Fig. 3(а)). For the case of the counterclockwise polarized radiation, good agreement between experiment and theory is attained only at one point I = 700 GW/cm2 (Fig. 3(а), grey solid line).
For a plate cut parallel to the Оу axis, similarly to the considered case of clock- and counterclockwise polarized radiation, there exist several sets of pairs χyyzz, χyyyz, at which experiment and theory agree well. Values of tensor χ components corresponding to the clockwise polarized radiation are shown by the solid curve in Fig. 4(b), where the dashed curve mark pairs χyyzz, χyyyz for the counterclockwise polarization.
Note that the principal goal of the current work is experimental verification of the supposition that accuracy of polarization transformation by means of wave plates reduces at В-integral larger than unity, rather than determination of unknown components of tensor χ. We have found the region of possible values of parameters χxxzz, χyyzz and χyyyz at which the experimental data agree with the theoretical calculations.
To conclude this section, we will make some estimates demonstrating why in our theoretical analysis we neglected spectrum broadening caused by cubic nonlinearity as the pulse was passing through the quartz zero-order quarter-wave plate. Let a laser pulse at the input to a plate having thickness L = 0.18 cm have the following shape А0(t) = 1/ch2(t/t0), where t0 = 28 fs. For the radiation with linear polarization coinciding with one of the eigen polarizations of the plate, a phase additive appears at the plate output defined by the nonlinear part of the refractive index: А(t) = 1/ch2(t/t0)⋅exp(– ikγNLI0(А0(t))2L), where I0 = 700 GW/cm2 is the maximum value obtained in the experiment. The spectra corresponding to the pulses А0(t) and А(t) are plotted in Fig. 5 .
We also plotted depolarization degree Γ0 versus λ specified by δ(λ) = k(nz(λ) – ny(λ))L. This function crosses zero at the point λ = 795 nm as the plate is fabricated to for radiation at this central wavelength. Integral of depolarization Γ0(λ) over the spectrum Аλ(λ) is much less than the Γ measured in the experiment (0.1-1%). Hence, the radiation spectrum broadening due to cubic nonlinearity may be neglected even at I0 = 700 GW/cm2.
The method of obtaining the circular polarization by means of a quarter-wave plate for a beam intensity at the plate 1 TW/cm2 and higher was investigated. A specific feature in this case is that, despite the small thickness of the wave plate, we have to consider the influence of ellipse polarization rotation caused by cubic nonlinearity on the polarization transformation.
Within the framework of the formulated problem we proposed a system of steady-state differential equations describing the propagation of laser radiation with arbitrary polarization in an anisotropic medium with cubic nonlinearity. Using this system of equations we calculated local depolarization degree Γ (fraction of energy in the polarization orthogonal to the polarization in the absence of nonlinearity) and polarization ellipticity at the output of a quarter-wave plate of crystalline quartz in the presence of cubic nonlinearity. Experiment on measuring depolarization at the output of a zero-order λ/4 plate 0.18 cm thick was conducted. It was found that the depolarization degree Γ is proportional to squared В-integral (nonlinear phase incursion). The maximum measured value of Γ was 1.1% at the laser radiation intensity of 700 GW/cm2. Good agreement between the experimental and the theoretical data was obtained.
Most promising for practical applications are wave plates made of KDP (or DKDP) crystals, especially for powerful lasers. According to [15,19], values of the components of nonlinear susceptibility tensor in these crystals are smaller than in crystalline quartz. Moreover, in contrast to quartz these crystals possess no gyrotropy, hence, whole λ/4 plates of arbitrary thickness may be cut from them. Study of the specific features of using such wave plates will be the next step in the work on this problem.
References and links
1. M. Borghesi, A. Schiavi, D. H. Campbell, M. G. Haines, O. Willi, A. J. Mackinnon, P. Patel, M. Galimberti, and L. A. Gizzi, “Proton imaging detection of transient electromagnetic fields in laser-plasma interactions,” Rev. Sci. Instrum. 74(3), 1688–1694 (2003). [CrossRef]
2. S. V. Bulanov and V. S. Khoroshkov, “Feasibility of using laser ion accelerators in proton therapy,” Plasma Phys. Rep. 28(5), 453–456 (2002). [CrossRef]
3. M. Roth, T. E. Cowan, M. H. Key, S. P. Hatchett, C. Brown, W. Fountain, J. Johnson, D. M. Pennington, R. A. Snavely, S. C. Wilks, K. Yasuike, H. Ruhl, F. Pegoraro, S. V. Bulanov, E. M. Campbell, M. D. Perry, and H. Powell, “Fast ignition by intense laser-accelerated proton beams,” Phys. Rev. Lett. 86(3), 436–439 (2001). [CrossRef]
5. A. Macchi, “A femtosecond neutron source,” Appl. Phys. B 82(3), 337–340 (2006). [CrossRef]
6. X. Q. Yan, C. Lin, Z. M. Sheng, Z. Y. Guo, B. C. Liu, Y. R. Lu, J. X. Fang, and J. E. Chen, “Generating high-current monoenergetic proton beams by a circularly polarized laser pulse in the phase-stable acceleration regime,” Phys. Rev. Lett. 100(13), 135003 (2008). [CrossRef]
7. D.-H. Kwon, K. Lee, S. H. Park, and Y. U. Jeong, “The effect of the transparency of an overdense plasma on proton beam generation by an intense ultra-short laser pulse,” J. Korean Phys. Soc. 49, 347–353 (2006).
8. S. N. Vlasov, V. I. Kryzhanovskiĭ, and V. E. Yashin, “Use of circularly polarized optical beams to suppress selffocusing instability in a nonlinear cubic medium with repeaters,” Sov. J. Quantum Electron. 12(1), 7–10 (1982). [CrossRef]
9. Y. B. Zel'dovich and Y. P. Raizer, “Self-focusing of light. Role of Kerr effect and striction,” JETP Lett. 3, 86–89 (1966).
10. P. D. Maker, R. W. Terhune, and C. M. Savage, “Intensity-dependent changes in the refractive index of liquids,” Phys. Rev. Lett. 12(18), 507–509 (1964). [CrossRef]
11. A. L. Berkhoer and V. E. Zakharov, “Self excitation of waves with different polarizations in nonlinear media,” Sov. Phys. JETP 31, 903–911 (1970).
12. D. Auric and A. Labadens, “On the use of circulary polarized beam to reduce the self-focusing effect in a glass rod amplifier,” Opt. Commun. 21(2), 241–242 (1977). [CrossRef]
13. M. S. Kuzmina, M. A. Martyanov, A. K. Poteomkin, E. A. Khazanov, and A. A. Shaykin, “Theoretical and experimental study of laser radiation propagating in a medium with thermally induced birefringence and cubic nonlinearity,” Opt. Express 19(22), 21977–21988 (2011). [CrossRef]
14. R. J. Pressley, Handbook of Lasers With selected Data on Optical Technology (Chemical Rubber Co, Cleveland, 1979).
15. R. L. Suthrland, Handbook of nonlinear optics, Marcel Dekker (New York, 2003).
16. H. R. Phillip, Handbook of Optical Constants of Solids (Academic Press, 1985).
17. S. S. Ballard, K. A. McCarthy, and W. L. Wolfe, Optical Materials for Infrared Instrumentation (IRIA-University of Michigan, Report #2389–11-S, 1959).
18. A. A. Babin, A. M. Kiselev, A. M. Sergeev, and A. N. Stepanov, “Terawatt femtosecond Ti:sapphire laser system,” Quantum Electron. 31(7), 623–626 (2001). [CrossRef]
19. N. G. Bondarenko, I. V. Eremina, and A. I. Makarov, “Measurement of nonlinear susceptibility tensor χijkl for DKDP crystals,” in Proceedings of the Five All-Union Conference on Laser Optics (Gosudarstvenny Optichesky Institut, Leningrad, 1987), pp. 37–40.