Abstract
We performed the first study of phase-matched second harmonic generation and difference frequency generation in the acentric nonlinear crystal NaI3O8. By fitting all the recorded phase-matching tuning curves, we determined the Sellmeier equations describing the wavelength dispersion of the ordinary and extraordinary principal refractive indices. We also measured the absolute value of the quadratic nonlinear coefficients and we calculated the conditions of supercontinuum generation in NaI3O8.
© 2017 Optical Society of America
1. Introduction
In 2007, we synthesized and studied the crystal structure of the new acentric crystal NaI3O8 [1]. Three years later, another group showed that this material is a negative uniaxial crystal with potential nonlinear optical properties. They performed powder second harmonic generation (SHG) measurements and calculations based on the density functional theory (DFT) [2]. Recent advances in crystal growth allowed us to get single crystals of very high quality suited for the first exhaustive study of linear and nonlinear optical properties described in the present letter. We recorded the transmission spectra in polarized light and determined the damage threshold in NaI3O8 slabs polished to optical quality. We directly measured phase-matching tuning curves of SHG and difference frequency generation (DFG) by using a method where NaI3O8 was cut as a cylinder polished on its curved surface [3]. From these experimental data we determined accurate Sellmeier equations describing the dispersion of the ordinary and extraordinary principal refractive indices, no and ne (no > ne), as a function of the wavelength. We also determined the absolute values of the two independent nonlinear coefficients d14 and d15 using a phase-matching technique [4].
2. Crystal growth, transmission spectra and damage threshold
NaI3O8 single crystals were grown by slow evaporation of Nitric Acid aqueous solutions (7M) at fixed temperature ranging between 60 and 70 °C. We obtained samples of several millimeter dimensions as the one shown in the insert of Fig. 1(a). While NaI3O8 displays a significant solubility in acidic solutions of lower molarity, such solutions were proved to be unstable toward the reduction of the iodate species into iodine. The later gets incorporated in the growing crystals giving them a yellow tint. Since NaI3O8 belongs to the tetragonal point group, its crystallographic frame (a, b, c) is orthonormal and fully coincides with the dielectric frame (x, y, z) [1, 2].
The transmission spectra depicted in Fig. 1(a) were recorded in polarized light through a 3-mm-thick slab. The two faces were cut perpendicularly to the x-axis, polished but uncoated. The linear polarization of light oriented successively perpendicularly and collinear with the z-axis, led to the ordinary and extraordinary transmission coefficients respectively. We used a Perkin-Elmer Lambda 900 spectrometer to record spectra as a function of the wavelength between 0.175 and 3.300 µm, and a Bruker FT-IR above 3.3 µm. Figure 1(a) shows that NaI3O8 is transparent between 0.32 and 6 µm despite strong absorption bands above 4 µm.
We determined the surface damage threshold of NaI3O8 compared with that of KTiOPO4 (KTP). The two slabs were illuminated by a 1.064-µm laser with 5-ns FWHM, 10-Hz repetition rate and 30-µm beam waist radius. NaI3O8 was damaged at an input energy of 140 J, corresponding to a peak power density of 0.79 GW/cm2. It was only 5 times lower than the damage threshold in KTP, which has been observed at 760 µJ i.e. 4.29 GW/cm2.
3. Phase-matching properties, Sellmeier equations and nonlinear coefficients
Under Kleinman symmetry, the second-order electric susceptibility tensor of the NaI3O8 crystal has two independent nonzero coefficients, i.e. dxyz = dxzy = dyzx = dyxz = dzxy = dzyx ( = d14); dxzx = dxxz = dzxx = - dyzy = - dyyz = - dzyy ( = d15), where d14 and d15 stand for the contracted notation [5]. The five following types of SHG and DFG processes can be phase-matched with a nonzero effective coefficient identically in the xz- and yz-plane: type I SHG, type II SHG , type I DFG type II DFG and type III DFG . Superscripts o and e stand for the ordinary and extraordinary waves, respectively. and (with k being o or e) are the fundamental and second harmonic wavelengths, respectively. and are the pump and signal input wavelengths, and is the idler wavelength generated by DFG, with the following relation of order: . The corresponding tuning curves were directly measured by using a method where the NaI3O8 crystal was cut as a cylinder with a diameter D = 4.13 mm and an acylindricity ΔD/D below 1%. The curved surface was polished to optical quality as shown in the insert of Fig. 1(b). We previously used this sample shape to study many nonlinear crystals, as RTP for example [3].
Using an X-ray backscattered Laue method, the cylinder was cut with its rotation axis oriented along the y-axis with a precision better than 0.5°. Then it was stuck on a goniometric head along this axis. When placed at the center of an Euler circle, the cylinder can rotate on it-self to access any direction of the xz-plane. Note that this plane is equivalent to yz-plane, since NaI3O8 belongs to the tetragonal point group. Using a focusing 100-mm-focal length lens properly placed, an incoming beam remains in normal incidence and propagates parallel to the diameter of the rotating cylinder as shown in Fig. 1(a) [3].
One incoming tunable beam was necessary for SHG. It was emitted by a 5-ns (FWHM) and 10-Hz repetition rate optical parametric oscillator (OPO) from Continuum Company. The OPO was pumped at 1.064 µm and tunable between 0.4 µm and 2.4 µm. The OPO beam and part of the 1.064-µm beam were collinearly combined for achieving DFG. They were put in spatial and temporal coincidence inside the cylinder using mirrors and the delay line shown in Fig. 1(b). All beams being linearly polarized, achromatic half-wave plates (HWP) provided their polarization rotation, in order to achieve the different types of phase-matching, i.e. type I, II and III. Input energies were measured at the entrance of the cylinder using a J4-09 Molectron pyroelectric joulemeter coupled to a beam splitter and a 50-mm focal length CaF2 lens. They were removed at the output of the cylinder by using a polarizer and a filter. Thus only beams generated by SHG and DFG were measured on the J3-05 Molectron joulemeter combined with a PEM531 amplifier. Phase-matching angles corresponding to maximum values of the conversion efficiency were directly read on the Euler circle with an accuracy of ± 0.5°. The corresponding phase-matching wavelengths were recorded with a precision of ± 1nm using a HR 4000 Ocean Optics spectrometer. Figures 2 show the measured types I and II SHG tuning curves, types II and III DFG being given in Figs. 3. For SHG, it was not possible to record phase-matching angles for fundamental wavelengths lower than 0.7 µm because the corresponding generated wavelengths were below the ultra-violet cut-off wavelength. Similarly, we did not study type I DFG tuning curve since the incoming tunable pump wavelength being below 0.5 µm, it could damage the surface of the cylinder. The tuning curves of Figs. 2 and 3 indicate that NaI3O8 allows phase-matching conditions over its entire transparency range, which is also the spectral range of reliability of its two principal refractive indices. These data can be used directly to cut the crystal at a phase-matching angle corresponding to the targeted parametric process and phase-matching wavelength. But the combination of the simultaneous fit of all the phase-matching curves with the magnitude of the ordinary refractive index no at a given wavelength, can lead to the determination of the Sellmeier equations [3]. We used no(λ = 0.671 µm) = 1.6 that we determined from the direct measurement of a Brewster angle of 58° in the xy-plane of NaI3O8 cut as a slab. For that purpose, we used the Levenberg-Marquardt algorithm encoded with Matlab. We tried different dispersion equations, the best result being obtained with the following dual oscillator model:
λ is in µm, and j stands for o or e. The Sellmeier coefficients Aj, Bj, Cj, Dj and Ej are given in Table 1. Since the precision of the measured phase-matching angles was of about ± 0.5°, the relative accuracy of the determined refractive indices is better than 10−4. It is corroborated in Figs. 2 and 3 by the very good agreement between our data and calculations. On the contrary, the comparison of our refractive indices at 1.064 µm with calculations using the model of ref. [2] that is shown in Table 2 highlights a strong discrepancy. This feature is not so surprising since the modelling of the first order electric susceptibility is very complicated.The following step was to determine the magnitude of the two independent nonlinear coefficients of NaI3O8, and. For that purpose, we measured conversion efficiencies in NaI3O8 slabs cut at critical phase-matching angles for types I-SHG and type II-SHG, and for type II-SHG in KTP that was taken as a reference [4]. We chose all fundamental wavelengths as near as possible, to get rid of the experimental setup spectral response. With ρe as the spatial walk-off angle, the corresponding effective coefficients are in NaI3O8 [4,5]:
Using Eq. (1) and Tab. 2, we calculated the corresponding phase-matching and walk-off angles in NaI3O8 at the phase-matching wavelength . We found for type I SHG that and ρe( = 1.38°; ρe( = 1.52° and ρe( = 1.51° in the case of type II SHG. We chose type II SHG in the xy-plane of KTP as a reference. The absolute value of the associated effective coefficient is = 2.43 pm/V, the corresponding phase-matching and walk-off angles being φPM = 23.1°, ρe( = 0.17° and ρe( = 0.29° [6]. We cut two NaI3O8 slabs 800 µm, 680 µm) and one KTP slab 800 µm) at the phase-matching angles given above. We took such a small interacting length in order to avoid any spatial walk-off attenuation. Since the incoming fundamental beam was focused through a 100-mm focal length CaF2 lens, the corresponding beam waist radius was around wo = 60 µm. It leads to a Rayleigh length 2xZR = 22 mm, which is much longer than the thickness of the slabs and ensures a parallel beam propagation. Types I and II SHG conversion efficiencies of NaI3O8, recorded relatively to KTP as a function of the fundamental wavelength are shown in Fig. 4: and correspond to types I and II SHG conversion efficiencies of NaI3O8 respectively, and stands for type II SHG conversion efficiency of KTP.
The absolute value of the effective nonlinear coefficient relatively to of KTP can be determined from the maximum of each curve shown in Fig. 4 and by using:
With
and stand for the ordinary refractive indices of NaI3O8 and KTP respectively. and are the extraordinary refractive indices. They are involved at the fundamental (a = 1) or second harmonic (a = 2) wavelengths, and or phase-matching angles. To and Te are the associated Fresnel coefficients. The spatial walk-off attenuations are negligible since we found: for type I SHG in NaI3O8, for type II SHG, and = 0.999 for type II SHG in KTP. The phase-matching wavelengths of NaI3O8 are = 1.080 µm for type I SHG and = 1.058 µm for type II SHG as shown in Figs. 4. They are slightly different from 1.064 µm, showing the difficulty of cutting very thin oriented slabs. Note that they are just slightly different from the calculated ones using our Sellmeier Eq. (1) and Tab. 1, which give the very good order of magnitude reached with the equations from the present work. From the maximum value of the curves shown in Fig. 4 and using Eqs. (4)-(7), we found that = 0.22 ± 0.02 pm/V and = 0.53 ± 0.05 pm/V. Then it comes from Eqs. (2) and (3): = 0.54 ± 0.08 pm/V and = 0.85 ± 0.13 pm/V. The discrepancy is strong when compared with calculations using the model of ref. [2] as shown in Tab. 2 since the modelling of the second-order electric susceptibility remains an open question. The corresponding Miller indices determined using Eq. (1) are [7]: , , and [6].Reliable NaI3O8-OPG tuning curves can be calculated using our Sellmeier equations. Thus, we found that a supercontinuum can be generated from type II OPG in the dielectric xz-plane of NaI3O8. The supercontinuum is the broadest one, ranging between 2.55 µm and 4.65 µm, when the pump wavelength is equal to 1.16 µm and the crystal oriented at θ = 18.5° from the z-axis as shown in Fig. 5. A supercontinuum can also be generated when the crystal is pumped by the Ti:Sapphire at 0.8 µm and the Nd:YAG at 1.064 µm, also shown in Fig. 5.
4. Conclusion
We measured for the first time to the best of our knowledge the transmission spectra in polarized light, the damage threshold, and the phase-matching conditions of SHG and DFG in the new acentric uniaxial NaI3O8 crystal. It allowed us to determine the Sellmeier equations and the absolute value of the two nonzero second-order nonlinear coefficients. Our data show that NaI3O8 is very attractive since using our Sellmeier’s equations we showed that a NaI3O8-based OPG can generate a supercontinuum in the mid-IR when pumped in the near-IR. We also bring a new input to refine the modelling of the second-order electric susceptibility in iodates compounds.
Funding
The authors thank the China Scholarship Council (CSC) for the financial support of Feng Guo.
References and links
1. D. Phanon and I. Gautier-Luneau, “Promising material for infrared nonlinear optics: NaI3O8 salt containing an octaoxotriiodate(V) anion formed from condensation of [IO3-] ions,” Angew. Chem. Int. Ed. Engl. 46(44), 8488–8491 (2007). [CrossRef] [PubMed]
2. X. Xu, C. L. Hu, B. X. Li, B. P. Yang, and J. G. Mao, “α-AgI3O8 and β-AgI3O8 with large SHG responses: polymerization of IO3 groups into the I3O8 polyiodate anion,” Chem. Mater. 26(10), 3219–3230 (2014). [CrossRef]
3. Y. Guillien, B. Ménaert, J. P. Fève, P. Segonds, J. Douady, B. Boulanger, and O. Pacaud, “Crystal growth and refined Sellmeier equations over the complete transparency range of RbTiOPO4,” Opt. Mater. 22(2), 155–162 (2003). [CrossRef]
4. I. Shoji, T. Kondo, A. Kitamoto, M. Shirane, and R. Ito, “Absolute scale of second-order nonlinear-optical coefficients,” J. Opt. Soc. Am. B 14(9), 2268–2294 (1997). [CrossRef]
5. B. Boulanger and J. Zyss, in International Tables for Crystallography, A. Authier, ed., Vol. D of Physical Properties of Crystals (Kluwer Academic, 2004), pp. 178–219.
6. B. Boulanger, J. P. Fève, G. Marnier, B. Ménaert, X. Cabirol, P. Villeval, and C. Bonnin, “Relative sign and absolute magnitude of d(2) nonlinear coefficients of KTP from second-harmonic-generation measurements,” J. Opt. Soc. Am. B 11(5), 750–757 (1994). [CrossRef]
7. R. C. Miller, “Optical second harmonic generation in piezoelectric crystals,” Appl. Phys. Lett. 5(1), 17–19 (1964). [CrossRef]