## Abstract

We report for the first time measurements and modelization of the angular distributions of absorption and fluorescence in a monoclinic crystal. Studies on Nd:YCOB revealed specific topologies with ombilics. These new data upgrade the knowledge on low symmetry crystal optics.

©2008 Optical Society of America

## 1. Introduction

The potentiality of new laser crystals relies on a complete characterization of their absorption and fluorescence properties. These two effects are governed by the imaginary part of the complex permittivity *$\widehat{\epsilon}$* that is described by a polar second rank tensor, the real part of *$\widehat{\epsilon}$* being related to the refractive index. When laser crystals belong to cubic crystal classes, like Nd:YAG for example [1], the three main values of *$\widehat{\epsilon}$* are equal, leading to optical properties that are independent to the crystal orientation. But, nowadays, there are promising laser materials that are anisotropic, like for example the uniaxial crystals Nd:YLiF_{4} and Nd:YVO_{4} [1], or the biaxial ones Yb:GdCOB [2] and Nd:YCOB [3]. For these crystals the tensors of the real and imaginary parts of *$\widehat{\epsilon}$* can have up to three different main values each, so that the absorption and fluorescence angular distributions are expected to be anisotropic, as it is the case for the refractive index [4]. Here we theoretically and experimentally considered the case of a monoclinic crystal class, and Nd:YCOB is investigated. We report, for the first time to the best of our knowledge, that the angular distributions of the two polarization components of the absorption coefficient and of the fluorescence power exhibit singularities of the same kind than those of the biaxial index surface. Furthermore, we show that the main values of absorption and fluorescence are not along the principal axes of the dielectric frame, but respectively along the principal axes of two new frames tilted from the dielectric frame.

## 2. State of the Art : absorption of anisotropic high symmetry crystals

A complete description of absorbing anisotropic crystals had been reported by introducing a complex permittivity *$\widehat{\epsilon}$*=*ε*+*jε*′ in the equations of the propagation of light, where *ε* is related to the refractive index and *ε*′ governs absorption [4]. The discussion had been confined to anisotropic crystals belonging to high symmetry crystal classes, *i.e.* hexagonal, tetragonal, trigonal or orthorhombic classes. In these cases, the real and imaginary parts of *$\widehat{\epsilon}$* are described by polar second rank tensors that are both diagonal in the dielectric frame [4]. The analysis was considerably simplified by assuming *ε*′≪*ε*. Then it relies on solutions of the complex Fresnel’s equation where real and imaginary parts can be analytically expressed. Here we focus on these solutions as a function of the spherical angle θ in the X-Z principal plane of the dielectric frame (X, Y, Z) of orthorhombic crystals, which are biaxial crystals. We consider the main values of the real and imaginary parts of the complex index *n̂* that are defined from the main values of the complex permittivity tensor by: ${\widehat{n}}^{2}=\widehat{\epsilon}{\epsilon}_{0}^{-1}$, where *ε*
_{0} is the free-space permittivity, leading to *n*
^{2}=*εε*
^{-1}
_{0} and 2*nn*′=*ε*′*ε*
^{-1}
_{0} when *n*′≪*n* is assumed, n and n′ being the real and imaginary parts respectively of *n*̂. They write, according to the perpendicular (⊥) and parallel (//) polarization components with respect to the X-Z plane [4]:

where *n _{1}*,

*n*and

_{2}*n*are the three main values of the real part of

_{3}*n̂*,

*i.e.*the three main refractive indices, and

*n*′

_{1},

*n*′

_{2}and

*n*′

_{3}are the three main values of the imaginary part of

*n*̂. The (⊥) polarization is defined as orthogonal to the X-Z plane: it is thus collinear with the Y-axis. The (//) component is located in the X-Z plane. Equations (1) describe the (⊥) and (//) polarization components of the two sections of the well-known index surface of a biaxial crystal in the X-Z plane of the dielectric frame. They are given in Fig. 1(a) and show the two Optical Axes (OA), which correspond by pair to the four ombilics of the index surface. Equations (2) describe the (⊥) and (//) polarization components of the angular distributions of the imaginary part of

*n̂*in the X-Z plane. They are given in Fig. 1(b) and exhibit two directions (OA′) where the (⊥) and (//) polarization components are equal, as in the case of the index surface. Figures 1(a) and 1(b) indicate that the (⊥) polarized angular distribution is a circle for both the real and imaginary parts of

*n̂*. But there is a difference between the (//) components of polarization, the angular distribution being elliptical for the real part, and “bilobar” for the imaginary one.

According to Eq. (1) and (2), a complete picture of the optical properties of absorbing anisotropic high symmetry crystals can be then obtained from measurements along the principal axes of the dielectric frame, which is completely linked up to the crystallographic frame for the high symmetry crystal classes mentioned above. Classically such measurements are performed with a polarized light beam properly focused in a polished slab oriented along the principal axes of the dielectric frame. It provides the determination of the main values of the refractive index and of the absorption, as well as the recording of their spectra. They can be discriminated by adjusting successively the orientation of the linear polarization of light along the (⊥) and (//) components associated with the considered direction of propagation. Note that fluorescence is also governed by the imaginary part of the complex index, so that the angular distributions of its (⊥) and (//) components are described by Eq. 2, in the case of high symmetry crystals.

## 3. Direct measurement of absorption and fluorescence in monoclinic crystals: method and experimental setup

Our goal was to perform direct measurements of the angular distributions of absorption and fluorescence of low symmetry crystals belonging to a monoclinic crystal class. A full access to the space is possible if the studied crystal is cut as a polished and oriented sphere. We report here such study for the monoclinic Nd:YCOB crystal cut as a sphere with a diameter of 7.44 mm. Nd:YCOB belongs to the *Cm* space group. Then the only symmetry element of this crystal is a mirror plane, *m*, which is perpendicular to the **b**-axis of the crystallographic frame (**a**, **b**, **c**). The angles between the crystallographic axes are: *β _{ab}*=90°,

*β*=90° and

_{bc}*β*=120.26° [5]. The Y-axis of the dielectric frame (X, Y, Z) is parallel to the

_{ac}**b**-axis, while the X-and Z-axes are not collinear with the

**a**- and

**b**-axes, exhibiting the following angles:

*γ*=24.7° and

_{aX}*γ*=13.4°, these angles having no observable dependence with wavelength [5]. We studied in the Nd:YCOB sphere the absorption at 0.812 µm on the one hand, and the fluorescence at 1.061 µm in on the other hand when pumping the sphere with a laser source at 0.812 µm. These two wavelengths are relative to the Nd

_{cZ}^{3+}ions doping the YCOB matrix.

We implemented the same experimental setup for both measurements, which top view is depicted in Fig. 2. It was coupled with a 10Hz-repetition rate CONTINUUM Panther Optical Parametric Oscillator (OPO) emitting a beam with a FWHM pulse duration of 5 ns which power is recorded on photodiode 1. The Nd:YCOB sphere was stuck along the Y-axis and placed on a rotation stage inserted behind the lens L_{1}. It was then possible to propagate the OPO beam parallel inside the sphere, successively over θ=360° in the X-Z plane that is also the symmetry mirror plane *m* of the crystal. The linear polarization of the OPO beam can be adjusted by using the half-wave plate L_{H} from the (⊥) component of polarization, along the Y-axis, to the (//) one oriented in the X-Z plane.

For the absorption measurements, we considered successively these two polarization components, and their transmission through the sphere was recorded on photodiode 2 coupled with lens L_{2}. For the fluorescence studies, the OPO beam was kept (⊥) polarized for any considered direction of propagation *θ* in the X-Z plane, in order to avoid any spatial walk-off or anisotropic effects. By this way, the power at 0.812 µm absorbed in the sphere remains the same for any direction of propagation. The pump beam was cut after the sample by a chromatic filter F (HR at 0.812 µm and HT at 1.061 µm) so that fluorescence power only is recorded through lens L_{2} on photodiode 2. The (//) and (⊥) polarization components of absorption and fluorescence were discriminated by using a Glan-Taylor polarizer P.

The (//) and (⊥) components of the absorption coefficient at 0.812 µm, *α*
_{//} and *α*
_{⊥} respectively, were determined from measurements of the total transmission coefficients, *T*
_{//,⊥}, and of Fresnel transmission relative to the entrance and exit surfaces of the sphere, *T*
^{F}
_{//,⊥}, according to
${\alpha}_{//,\perp}=-\frac{1}{L}\mathrm{ln}\left(\frac{{T}_{//,\perp}}{{T}_{//,\perp}^{F}}\right)$
where *L* is the sphere diameter. The (//) and (⊥) components of the power of fluorescence at 1.061 µm, *P*
_{//} and *P*
_{⊥} respectively, were directly recorded and corrected by Fresnel transmission relative to the exit surface of the sphere. Absorption and fluorescence measurements were performed as a function of the propagation angle *θ*, successively over 360°, in the X-Z plane.

## 4. Experimental results and their modelization

Figure 3(a) gives the angular distribution of the (⊥) and (//) components of polarization of the absorption coefficient, *α*
_{⊥}(*θ*) and *α*
_{//}(*θ*), at 0.812 µm of Nd:YCOB. It is the first time, to the best of our knowledge, that this kind of measurement is performed. It shows that *α*
_{⊥}(_{θ}) remains constant as a function of the direction of propagation *θ*, and *α*
_{//}(*θ*) exhibits a bi-lobar angular distribution, which is in a very good agreement with the behaviour of the calculated patterns of Fig. 1(b) corresponding to orthorhombic crystals. But the symmetry axes of the extraordinary angular distribution of Fig. 3(a), labelled as X’_{abs} and Z’_{abs}, do not correspond to the axes X and Z of the dielectric frame that are the symmetry axes of the extraordinary pattern of Fig. 1(b). So for monoclinic crystals, which is the topics of our study, we have to define a new frame (X′_{abs}, Y′_{abs}, Z′_{abs}), the Y′_{abs}-axis being collinear with the Y-axis of the dielectric frame. We called the “Absorption Axes” (AA), the directions along which the (//) and (⊥) polarized absorption coefficients are equal, the angle between these two axes being written 2V_{abs}. We found the same type of angular distribution for the fluorescence at 1.061 µm, as shown in Fig. 3(b) where the powers are normalized to that of the (⊥) component emitted along the Z-axis, *i.e.*
*P*
_{⊥}(*θ*=0°). The angular distribution of the (⊥) polarized component of fluorescence defined as *p*
_{⊥}(*θ*)=*P*
_{⊥}(*θ*)/*P*
_{⊥}(*θ*=0°) is described by a circle, while that of the (//) one, *p*
_{//}(*θ*)=*P*
_{//}(*θ*)/*P*
_{⊥}(*θ*=0°), is a bi-lobar pattern, which in accordance with the fact that fluorescence is also governed by the imaginary part of the complex index, as for absorption. Fig. 3(b) also shows that the symmetry axes of the (//) fluorescence pattern, *i.e.* X′_{fluo} and Z′_{fluo}, are different than X′_{abs} and Z′_{abs} on the one hand, and than X and Z on the other hand, while Y′_{fluo} is collinear with Y′_{abs} and Y. The two axes for which the two components of fluorescence are equal are called the “Fluorescence Axes” (FA), the angle 2V_{fluo} between these two axes being different than the angle 2V_{abs} between the two “Absorption Axes” (AA).

Then our experimental results suggest the introduction of two specific principal frames for absorption and fluorescence of a monoclinic crystal, *i.e.* (X′_{abs}, Y′_{abs}, Z′_{abs}) and (X′_{fluo}, Y′_{fluo}, Z′_{fluo}) respectively, these frames being different and also not in coincidence with the dielectric frame (X, Y, Z). As a consequence, the imaginary part of the *$\widehat{\epsilon}$* tensor is diagonal only in (X′_{abs}, Y′_{abs}, Z′_{abs}) for absorption, and in (X′_{fluo}, Y′_{fluo}, Z′_{fluo}) for fluorescence. Then the model developed in reference [4] for high symmetry crystals, describing the imaginary part of *$\widehat{\epsilon}$* as diagonal in the dielectric frame, does not apply to monoclinic crystals. We thus wrote a new model able to describe our measurements, considering the imaginary part of *$\widehat{\epsilon}$* to be not diagonal in the dielectric frame:

*ε _{xx}*,

*ε*and

_{yy}*ε*are the three main values of the real part of the permittivity, and

_{zz}*ε*′

_{xx},

*ε*′

_{yy},

*ε*′

_{zz},

*ε*′

_{xz}=

*ε*′

_{zx}are the coefficients of the imaginary part. According to Eq. (3), the projection of the light propagation equation along the three principal axes of the dielectric frame (X, Y, Z) leads to the following single equation that relates the complex index

*n̂*to the different coefficients of the complex permittivity tensor, and that writes for any direction of propagation in the X-Z plane,

*i.e.*for 0°<

*θ*<360°

*ϕ*=0°:

The angle *θ* is between the direction of propagation and the Z-axis. The resolution of Eq. (4) provides two complex solutions. Assuming *n*′≪n, the (⊥) and (//) components of the real part of the two solutions are given by Eq. (1), the same ones than for high symmetry crystals. But the equations relative to the (⊥) and (//) components of the imaginary part are specific to monoclinic crystals, and write in the X-Z plane:

The two Eq. (5) confirm circle and bi-lobar patterns for the angular distributions of the (⊥) and (//) components of the imaginary part of *n̂* respectively. They were used to interpolate our experimental data, by assuming *α*
_{//,⊥}(*θ*)=4*πn*′_{//,⊥}(*θ*)*λ*
^{-1} for the absorption at 0.812 µm, and *p*
_{//,⊥}=*p*
_{//,⊥}(*θ*)/*P*⊥(*θ*=0°)=*n*′_{//,⊥}(*θ*)/*n*′_{⊥}(*θ*=0°) for the fluorescence at 1.061 µm, the fitting parameters being: *ε*′_{xx}, *ε*′_{yy}, *ε*′_{zz} and *ε*′_{xz}=*ε*′_{zx.}. The interpolations are shown in Fig. 3(a) and Fig. 3(b), and the agreement between the model and the experiment is very satisfying. The interpolation of absorption data allowed us to determine the tilt angle between the absorption frame (X′_{abs}, Y′_{abs}, Z′_{abs}) and the dielectric frame (X, Y, Z), *i.e.*
*δθ _{abs}*=31.1±0.7°. We got

*δθ*=-6.4±0.9° for the tilt angle between the fluorescence frame (X′

_{fluo}_{fluo}, Y′

_{fluo}, Z′

_{fluo}) and the dielectric frame (X, Y, Z). The interpolation of our data with Eq. (5) also lead to 2V

_{abs}=72°±4° and 2V

_{fluo}=64°±3° as the angles between the two “Absorption Axes” (AA), and between the two “Fluorescence Axes” (FA), respectively. Note that they are also different from the angle between the two “Optical Axes” (OA),

*i.e.*2V

_{z}=120.6° [5]. From the knowledge of the orientation of the absorption and fluorescence frames, we can determine the main values of the absorption and fluorescence in their specific frames. We got for the main values of the absorption coefficient at 0.812 µm:

*α*

_{1}=2.1±0.8 cm

^{-1},

*α*

_{2}=3.5±0.4 cm

^{-1}and

*α*

_{3}=6.2±0.8 cm

^{-1}, where the indices 1, 2 and 3 refer to X′

_{abs}, Y′

_{abs}and Z′

_{abs}respectively. For the fluorescence at 1.061 µm we got for the main values: p

_{1}=0.65±0.04, p

_{2}=1.00±0.02 and p

_{3}=2.01±0.04, where the indices 1,2 and 3 correspond to X′

_{fluo}, Y′

_{fluo}and Z′

_{fluo}respectively. We see that these values are different than the diagonal values that are determined from measurements relative to the dielectric frame,

*i.e.*:

*α*=3.2±0.4 cm

_{1}^{-1},

*α*=3.5±0.4 cm

_{2}^{-1}and

*α*=5.1±0.6 cm

_{3}^{-1}for absorption, and p

_{1}=0.67±0.02 cm

^{-1}, p

_{2}=1.00±0.02 cm

^{-1}and p

_{3}=1.99±0.02 cm

^{-1}, where the indices 1, 2 and 3 correspond in that case to X, Y and Z respectively. It corroborates that the main values of absorption and fluorescence can be obtained only along the principal axes of the absorption and fluorescence frames respectively, but not along those of the dielectric frame, as classically assumed prior to the present work.

## 5. Conclusion

We report, for the first time to the best of our knowledge, measurements and calculations relative to the angular distributions of the (⊥) and (//) components of the absorption and fluorescence in a monoclinic crystal. We investigated Nd:YCOB and showed that it is necessary to introduce new specific frames leading to the main values of absorption and fluorescence. The practical consequence to fully characterize absorption and fluorescence properties is that it is necessary to measure their main values in these frames instead of doing it in the dielectric frame. This feature is also of prime importance for devices purpose, for which it is necessary to optimize the absorption of the pump as well as the fluorescence. Then the choice of the direction in which the monoclinic laser material has to be cut must take into account the difference of symmetry axes of the two involved properties. This choice is even more complicated for self-doubling, which also occurs in Nd:YCOB, where a third property has also to be optimized, *i.e.* the phase-matching of second harmonic generation of the fluorescence, which principal frame is the dielectric frame. Finally, the part of our work devoted to absorption is of general interest for any monoclinic crystal, since absorption has to be considered and minimized in a lot of situations, like for example: second and third harmonic generations in un-doped YCOB [6], parametric up- and down-conversions in BiB_{3}0_{5} (BiBO) [7], or photorefractivity in Sn_{2}P_{2}S_{6} [8, 9].

## Aknowledgment

This work was supported by the collaborative National French framework “Bulk Crystals and Optical Devices” of Centre National de la Recherche Scientifique.

## References and links

**1. **A. I. Zagumennyi, V. A. Mikhailov, and I. A. Shcherbakov, “Rare earth ion lasers - Nd^{3+},” Handbook of Laser Technology and Applications, Vol. II, Edited by C.E. Webb and J.D.C. Jones, IOP Publishing, 353 (2004).

**2. **F. Mougel, K. Dardenne, G. Aka, A. Kahn-Harari, and D. Vivien, “Ytterbium-doped Ca_{4}GdO(BO_{3})_{3}: an efficient infrared laser and self-frequency doubling crystal,” J. Opt. Soc. Am. B **14**, 164 (1999). [CrossRef]

**3. **J. M. Eichenholz, D. A. Hammons, L. Shah, Q. Ye, R. E Peale, M. Ridcharson, and B. H. T. Chai, “Diodepumped self-frequency doubling in a Nd^{3+}:YCa_{4}O(BO_{3})_{3} laser,” Appl. Phys. Lett. **74**, 1954 (1999). [CrossRef]

**4. **Wolf and Born
, *Principles of Optics* (Oxford Pergamon press, 1965).

**5. **P. Segonds, B. Boulanger, B. Ménaert, J. Zaccaro, J. P. Salvestrini, M. D. Fontana, R. Moncorgé, F. Porée, G. Gadret, J. Mangin, A. Brenier, G. Boulon, G. Aka, and D. Pelenc, “Optical characterization of Ca_{4}YO(BO_{3})_{3} and Nd:Ca_{4}YO(BO_{3})_{3} crystals,” Opt. Mater. **29**, 975 (2007). [CrossRef]

**6. **C. A. Ebbers, A. J. Bayramian, R. Campbell, R. Cross, B. L. Freitas, Z. Liao, K. I. Schaffers, S. Sutton, J. A. Caird, C. P. J. Barty, Y. Fei, and B. H. T. Chai, “High Average Power Frequency Conversion with Large Aperture YCOB,” in *Proceedings of Advanced Solid State Photonisc (ASSP) 2008*, WD3 (2008).

**7. **H. Hellwig, J. Liebertz, and L. Bohaty, “Linear optical properties of the monoclinic bismuth BiB_{3}O_{6},” J. Appl. Phys. **88**, 240 (2000). [CrossRef]

**8. **A. M. Weber, G. Von Bally, A. Shumelyuk, and S. Odoulov, “Reflection-type photorefractive gratings in the hypotiodiphosphate,” Appl. Phys. B Lasers Opt. **74**, 29 (2002). [CrossRef]

**9. **A. Shumelyuk, A. Volkov, A. Selinger, M. Imlau, and S. Odoulov, “Frequency-degenerate nonlinear light scattering in low-symmetry crystals,” Opt. Lett. **33**, 150 (2008). [CrossRef] [PubMed]