This paper reports an efficient method of using a chirped grating to reduce the phase-matching sensitivity and increase the output power of continuous wave second harmonic generation in a MgO doped LiNbO3 waveguide. It was found that using a chirped grating with decreasing periods from 6.96 μm to 6.95 μm significantly improved the phase mismatch tolerance and increased the second harmonic generation output power from 0.64 W (using uniform grating) to 1.13 W. This design was shown to lead to a phase matching insensitive waveguide. This indicates that there is a possibility of obtaining higher output power with the MgO doped LiNbO3 waveguide.
© 2017 Optical Society of America
Second harmonic generation (SHG) has received much attention because it provides an efficient way for green and blue light generation [1–3], high-resolution imaging [4–6] and optical signal processing [7–9]. SHG occurs when a laser beam incidents upon a material with a nonzero second order susceptibility and has been demonstrated in various materials [10–12]. Among these materials, periodically poled 5 mol.% MgO doped LiNbO3 (PPMgLN) is the most promising one for realizing highly efficient and compact blue and green lasers with moderate output power and good beam quality. It is of great interest because of the combined advantage of a large nonlinear coefficient and a good compatibility with waveguide fabrication process [13–15]. Two-dimensional spatial confinement provided by the waveguide significantly reduces the power required for highly efficient SHG . However, it has been found that the SHG process using this material is inhibited for applications with relatively high power [14, 17].
The inhibition can be interpreted by the dephasing theory that absorptions induced inhomogeneous temperature rising and laser induced inhomogeneous refractive index changing result in considerable inhomogeneous phase mismatch in the irradiated zone. The phase mismatch will dramatically decrease the energy conversion from fundamental wave to SHG wave. For example, our past work reveals that inhibition of high power SHG in a PPMgLN waveguide is a result of significant phase mismatch induced by the photorefractive effect and the thermal effect along their radiated zone . Great efforts have been made to minimize the inhibition of high power second harmonic generation in PPMgLN [19–26].
One way to minimize the inhibition is to broaden the bandwidth of second harmonic generation by chirped gratings. Basedon nonlinear coupled mode equations and undepleted pump approximation, Suhara et al. found that phase-matching bandwidth can be effectively broadened by appropriate design of chirped gratings .Since realizing a broadband device with a linearly chirped grating needs to ensure the smallness of change(~100 picometer) in the grating periods, Tehranchi et al. proposed a step-chirped grating which is engineered to highly broaden the bandwidth with larger step changes in the chirped period . Although chirped gratings can increase the spectral or thermal bandwidth of the SHG, the conversion of fundamental wave into SHG wave becomes inefficient. Baranova et al. proposed an adiabatic SHG to resolve the bandwidth-efficiency tradeoff based on undepleted pump approximation . Based on an analytical theory of SHG in a generic nonuniform second order nonlinear medium, Yaakobi et al. showed that by varying the properties of the medium gradually enough, the system can enter an autoresonant state in which the phases of the fundamental wave and of the SHG wave are locked .Subsequently, comprehensive physical models of adiabatic three-wave mixing were developed without using the undepleted pump approximation [24, 25]. The removal of undepleted pump restriction enables the exploration of adiabatic process in the fully nonlinear dynamic regime of nonlinear optics. Recently, Leshem et al. experimentally studied fully nonlinear adiabatic process by implementation in type I and type II SHG in two KTP nonlinear crystals .In an adiabatic SHG, the initial phase mismatch should be large enough so that nonlinear coupling is negligible and the rate of variation of the phase mismatch is sufficiently slow .To realize full conversion in adiabatic SHG, the signal intensities should be several orders of magnitude higher than the damage limits signal intensities under continuous-wave (CW) illumination, which makes the adiabatic CW SHG inefficient. Besides, these studies focused on a SHG with phase mismatch much smaller than the bandwidth. Although the bandwidth of the SHG can be very large, the efficiency is very low for CW operation due to the low damage limits signal intensity, which is ~2 MW/cm2 at 532 nm . However, little attention has been paid to the selection of an appropriate linearly chirped grating for minimizing the inhibition of CW SHG to improve the conversion efficiency.
The present paper presents a theoretical investigation on phase mismatch insensitive high power CW SHG by a coupled model in which phase mismatch from thermal dephasing and photorefractive effect are included. On basis of this model, it then describes the use of chirped gratings with linearly decreasing grating periods to reduce the sensitivity of PPMgLN waveguides to phase-matching conditions and maintain the SHG conversion efficiency. We found that the chirped gratings significantly reduced the sensitivity of the waveguide to phase mismatch in which the full width at half maximum (FWHM) bandwidth of the temperature tuning curve increased from 2.4 °C to 11.7 °C and the SHG output power increased from 0.64 W to 1.13 W.
2. Materials and methods
The current investigation involved analyzing the use of linearly chirped gratings to reduce the sensitivity of a PPMgLN waveguide to phase-matching and meanwhile increase the SHG output power. The PPMgLN waveguide has been widely used for highly efficient SHG [20, 28]. Frequency doubling of a 1064 nm CW laser with a PPMgLN waveguide was modeled and then calculated with a three dimensional coupled model equation. This method has been used for analyzing the inhibition phenomena occurred in PPMgLN waveguides , which can simulate the effect of photorefraction and heating upon the SHG. In this study, we chose a similar waveguide as the one we employed in the previous work. Thus, we can compare the results obtained in this study with that obtained in the previous one. As shown in Fig. 1(a), the waveguide is a periodically poled Y-cut 5% mol. MgO doped LiNbO3 ridge waveguide [17, 18]. A polished PPMgLN crystal is bonded directly on the LN substrate. An airgap boundary, formed under the ridge structure, separates the PPMgLN waveguide from the LN substrate and performs strong confinement in the depth direction. The waveguide has a length l of 11.5 mm, a width of 5 μm and a height of 4 μm. It was chosen because notable inhibition of SHG has been observed for a fundamental power higher than 0.6 W . The refractive index of the LN substrate can be found in .The normalized conversion efficiency used in this paper is 199%/W cm2 which was determined through an experiment . During the fabrication of the periodic structures, errors may be introduced in it . Period errors that may affect a temperature acceptance bandwidth can be controlled very well with the rapid development of PPLN . Therefore, we assumed that period errors were well controlled in this study. Fabrication errors that only change the conversion efficiency were taken into account in the normalized conversion efficiency in our calculation.
A waveguide with a linearly chirped grating was designed and then calculated for the SHG. As is known, the increase in SHG power along the z direction will result in an increase refractive index difference between SHG wave and fundamental wave [17, 18]. According to quasi-phase matching condition, the wave vector mismatch caused by the material dispersion should be compensated by the grating wave vector . Hence, a larger grating wave vector is required to compensate the wave vector mismatch caused by material dispersion in the rear of the waveguide. The waveguide was designed to have linearly decreasing grating periods so that the grating has relatively larger wave vectors in the rear of the waveguide. This makes the grating could better compensate the wave vector mismatch of the whole waveguide. As shown in Fig. 1(a), the grating periods linearly vary from ΛF at the front facet to ΛE at the end facet. Since the wave vector mismatch is mainly caused by the second harmonic wave , the grating period ΛF at the front facet is set to 6.96 μm aiming to satisfy the quasi-phase matching condition at a temperature of 34 °C. The grating period ΛE at the end facet of the waveguide varied from 6.80 μm to 6.95 μm as shown in Fig. 1(b).The wavelength and temperature dependent Sellmeier equation for 5% MgO-doped CLN was taken from . It is valid for0.5–4 μm spectral range and 20–200 °C temperature range.
3. Results and discussions
Data obtained in previous studies using models with no phase mismatch approximation indicates that sufficient intense signal are required to achieve high efficiency in broadband SHG using linear chirped poling periods [20–26]. According to Leshem et al., a pump intensity as high as ~486 MW/cm2 was required to achieve a SHG conversion efficiency of 60% in a periodically poled KTP crystal with poling periods varied between 8.6 μm to 9.2 μm . In our study, phase mismatch insensitive high power SHG was investigated by using a 3D thermal-optical model in which phase mismatch from thermal dephasing and photorefractive effect are taken into account . This model has also been proven by comparing the simulation results with the experimental results in that paper. The calculated SHG powers were compared with that of the waveguide with uniform grating periods. Figure 1(b) shows the dependence of SHG power on the grating period ΛE at a fundamental power of 1.5 W. Each of the SHG power shown in Fig. 1(b) was maximized by optimizing the waveguide temperature.
As can be seen in Fig. 1(b), the SHG power first increased and then decreased with the increasing ΛE, and a peak SHG power which is around 1.13 W appeared at ΛE = 6.95 μm. The corresponding SHG intensity was below the damage threshold (~2 MW/cm2) of 5 mol.% MgO doped LiNbO3 at 532 nm . As the ΛE decreased to 6.94 μm, the SHG power declined to ~1.08 W. The SHG power decreased with the decreasing ΛE. At ΛE = 6.80μm, it declined to 0.22 W. To evaluate the effect of a chirped grating upon the SHG output power, we compared the SHG power versus fundamental power of the waveguide WL with ΛE = 6.95 μm with that of the waveguide WU with uniform grating period as shown in Fig. 2(a). Since the damage limits signal intensity is ~2 MW/cm2 , the maximum pump power was limited to 1.5 W. The curve of SHG power versus fundamental power for the waveguide WU below a fundamental power of 1 W was in good agreement with experimental results [17, 18].
In Fig. 2(a), it can be seen that the SHG power generated by the waveguide WL with the linearly chirped grating was significantly higher than that by waveguide WU using uniform grating for a fundamental power above 0.7 W. As the SHG process was inhibited by increasing phase mismatch in the WU for fundamental above 0.3 W, the conversion of the energy from the fundamental wave to SHG wave became less efficient . The SHG process of the WL was almost not affected by the phase mismatch and the SHG power increased steadily with increasing fundamental power as compared with that of WU. Therefore, the difference between the SHG power of WL and WU became greater with the increasing fundamental power. At a fundamental power of 1.5 W, the SHG power of WL increased to 1.13 W which was ~77% higher than that of WU.
The prior work done by Jedrzejczyk et al. indicates that the SHG using uniform grating is quite sensitive to temperature variation at high fundamental power . For example, it showed that the FWHM width of the temperature tuning curve of SHG at Pω = 0.747 W was ~45.8% narrower than that of SHG at Pω = 0.162 W . The narrowed FWHM width was a result of the photorefractive effect and the thermal effect induced imhomogeneous refractive index distribution . In our study, the use of chirped grating could broaden the FWHM width of the temperature tuning curve by improving the phase mismatch tolerance at high pump powers. Figure 2(b) shows the temperature tuning curves obtained using the waveguide WL.
As can be seen, the SHG at a fundamental power of 1.5 W had a much larger FWHM width which was ~4.5 times of that at a fundamental power of 0.2 W. Although the temperature tuning curve of the waveguide WL deviated from the sinc2 distribution of that in a waveguide with uniform grating, it had a nearly symmetric distribution. The FWHM width of temperature tuning curve of the waveguide WL at Pω = 0.2 W is ~2.6 °C which is higher than that (~2.4 °C) of the waveguide WU at Pω = 0.167 W . However, at high pump powers, a striking difference between the temperature tuning curves of the waveguide WL and WU was noted when the width and the shape of the tuning curves were considered. The waveguide WU has a narrowed FWHM width at a high pump power, whereas the waveguide WL has a significantly broadened FWHM width at a high pump power. As shown in Fig. 2(b), the waveguide WL has a temperature tuning curve with a FWHM width up to 11.7 °C which was ~4.5 times of that at a fundamental power of 0.2 W. As compared with that (~2.4 °C) of the waveguide WU, the FWHM width of temperature tuning curve of the waveguide WL is ~4.9 times. The result indicates that the chirped grating significantly improves the tolerance of the waveguide to phase mismatch. To find out how the chirped grating improved the SHG output power, we plotted the wave vector mismatch caused by material dispersion Δkn of the waveguide WL, the total wave vector mismatch ΔkT of the waveguide WL and WU and evolving of the SHG process along z axis of the waveguide WL in Fig. 3.
As can be seen in Fig. 3(a), the wave vector mismatch caused by material dispersion Δkn of the waveguide WLis inhomogeneous in the xy plane. At the plane of z = l/2 (5.75 mm), the peak wave vector mismatch which was ~9034 cm−1 appeared at x = 0.0 μm and y = 2.5 μm. In the xy plane, Δkn decreases as it was far away from the peak point. At the edge, it reached the minimum value of ~9028 cm−1. In a quasi-phase matched SHG, the efficiency is determined by total wave vector mismatch ΔkT which is equal to Δkn minus wave vector of a grating kG. According to the quasi-phase matching condition (ΔkT = 0), it is most desirable to achieve quasi-phase-matching by making the kG close to the Δkn. kG of the waveguide WL at z = l/2 was ~9030 cm−1 which is close to the intermediate value of the Δkn at that position. The peak wave vector mismatch of the waveguide WL increased to 9036 cm−1 at z = l as shown in Fig. 3(b). kG of the waveguide WL at z = l was ~9036 cm−1 which was equal to peak wave vector mismatch at that position. As shown in Fig. 3(c), the offset of ΔkT of WU and WL has been optimized to obtain maximum outputs by tuning the temperature. However, total wave vector mismatch imhomogeneities cannot be reduced in the waveguide WU. As can be seen, ΔkT of the waveguide WU increased ~6.0 cm−1 from z = 0.0 mm to z = 4.4 mm. The increased ΔkT inhibited the SHG process to a large extent. The ΔkT reach the peak at z = 4.4 mm and decreasing slightly from this point to end facet. Although, the variation in Δkn was small (0.05%), the variation in the ΔkT was significant. As with the tapered grating, kG also increased along the beam propagation direction. As can be seen, ΔkT of the waveguide WL increased only ~1.3 cm−1 from z = 0.0 mm to at z = 3.1 mm. Therefore, the total wave vector mismatch imhomogeneities is reduced to a large extent. The ΔkT reach the peak at z = 3.1 mm and decreasing slightly from this point to end facet. Although the SHG power of the waveguide WL was ~76% higher than that of the waveguide WU, the ΔkT of it was only slightly (< 27%) higher than that of the waveguide WU in the second half of the waveguide. Therefore, the fact that kG of the chirped grating grows in phase with Δkn makes the SHG process in the WL more efficient. There have ~1600 domains in the waveguide. The SHG output is a coherent superposition of the SHG wave generated at each domain. With reduced total wave vector mismatch inhomogeneities, the quasi-phase matching temperatures of the domains are close to each other at a high SHG power. Then, a wider temperature acceptance bandwidth can be obtained.
Data obtained in previous studies indicated that the strong inhibition of SHG began at the middle of the waveguide with uniform grating. According to our previous work, the SHG power reached the peak at the middle part of waveguide and declined in the rear of the waveguide . In this study, the inhibition was suppressed by using the chirped grating. The evolving of P2ω in the waveguide WL is shown in Fig. 3(d).
As can be seen, the SHG process was almost not inhibited. A wider temperature acceptance bandwidth indicates that the waveguide WL is insensitive to phase mismatch as compared with the WU. It makes the WL has a better performance than the WU despite that the WL has the higher phase mismatch. The SHG power grew steadily in the front part of the waveguide, which was consistent with that in the waveguide with uniform grating. However, in the rear of the waveguide, the evolving of the SHG power was not consistent with that with uniform grating. It can be observed that the SHG power kept growing in the rear of waveguide. Although the SHG power grew with a lower speed in middle part of the waveguide than that in the front part, the growing speed was recovered at the rear of the waveguide. Since the total wave vector mismatch did not increase exactly linearly along the waveguide, the chirped grating could not compensate the phase mismatching completely. There will be some residual total wave vector mismatch along the waveguide. Nevertheless, these results suggest that a chirped grating may provide a SHG waveguide that is insensitive to phase mismatch as compared with a waveguide with uniform grating, and thereby increasing the SHG output power.
Prior work has reported that the effectiveness of chirped grating in broadening the bandwidth of SHG. Leshem et al. , for example, demonstrated that thermal acceptance bandwidth can be increased to 100 °C with chirped poling periods. In the experiment, a pump intensity as high as 290 MW/cm2 was used to achieve a conversion efficiency of 60%.However, the damage limits the pump intensity to several orders of magnitude lower than the intensity required for efficient conversion for CW lasers [27, 29].
In this study, we investigated phase mismatch insensitive high power SHG by using a 3D thermal-optical model in which phase mismatch from thermal dephasing and photorefractive effect are taken into account. We found that a chirped grating with linearly decreasing periods can reduce the sensitivity of SHG to phase-matching and increase SHG output power. By a tapered grating with decreasing periods from 6.96 μm to 6.95 μm, the FWHM width of the temperature tuning curve was broadened to 11.7 °C (the waveguide WL) from 2.4 °C (the waveguide WU) and the SHG output power was increased to ~1.13 W from ~0.64 W at a fundamental power of 1.5 W. These findings extend those of Leshem et al. , confirming that reduction of the inhibition of SHG can be achieved by engineering the phase-matching condition at a moderate signal intensity. Our result provides compelling evidence for phase-matching insensitive SHG and suggests that this approach appears to be effective in increasing SHG power that is required in many applications.
National Natural Science Foundation of China (NSFC) (61605136 and 11204205); Natural Science Foundation of Taiyuan University of Technology (2013Z029).
References and links
1. T. Kobayashi, D. Akamatsu, Y. Nishida, T. Tanabe, M. Yasuda, F. Hong, and K. Hosaka, “Second harmonic generation at 399 nm resonant on the 1S0 - 1P1 transition of ytterbium using a periodically poled LiNbO3 waveguide,” Opt. Express 24(11), 12142–12150 (2016). [PubMed]
2. L. Chang, Y. Li, N. Volet, L. Wang, J. Peters, and J. E. Bowers, “Thin film wavelength converters for photonic integrated circuits,” Optica 3(5), 531–535 (2016).
3. D. Lehr, J. Reinhold, I. Thiele, H. Hartung, K. Dietrich, C. Menzel, T. Pertsch, E. B. Kley, and A. Tünnermann, “Enhancing second harmonic generation in gold nanoring resonators filled with lithium niobate,” Nano Lett. 15(2), 1025–1030 (2015). [PubMed]
4. Y. Nakayama, P. J. Pauzauskie, A. Radenovic, R. M. Onorato, R. J. Saykally, J. Liphardt, and P. Yang, “Tunable nanowire nonlinear optical probe,” Nature 447(7148), 1098–1101 (2007). [PubMed]
5. W. R. Zipfel, R. M. Williams, R. Christie, A. Y. Nikitin, B. T. Hyman, and W. W. Webb, “Live tissue intrinsic emission microscopy using multiphoton-excited native fluorescence and second harmonic generation,” Proc. Natl. Acad. Sci. U.S.A. 100(12), 7075–7080 (2003). [PubMed]
6. E. Brown, T. McKee, E. diTomaso, A. Pluen, B. Seed, Y. Boucher, and R. K. Jain, “Dynamic imaging of collagen and its modulation in tumors in vivo using second-harmonic generation,” Nat. Med. 9(6), 796–800 (2003). [PubMed]
7. K. L. Seyler, J. R. Schaibley, P. Gong, P. Rivera, A. M. Jones, S. Wu, J. Yan, D. G. Mandrus, W. Yao, and X. Xu, “Electrical control of second-harmonic generation in a WSe2 monolayer transistor,” Nat. Nanotechnol. 10(5), 407–411 (2015). [PubMed]
8. B. A. Ruzicka, L. K. Werake, G. Xu, J. B. Khurgin, E. Y. Sherman, J. Z. Wu, and H. Zhao, “Second-harmonic generation induced by electric currents in GaAs,” Phys. Rev. Lett. 108(7), 077403 (2012). [PubMed]
9. W. Ding, L. Zhou, and S. Y. Chou, “Enhancement and electric charge-assisted tuning of nonlinear light generation in bipolar plasmonics,” Nano Lett. 14(5), 2822–2830 (2014). [PubMed]
10. A. Rao, M. Malinowski, A. Honardoost, J. R. Talukder, P. Rabiei, P. Delfyett, and S. Fathpour, “Second-harmonic generation in periodically-poled thin film lithium niobate wafer-bonded on silicon,” Opt. Express 24(26), 29941–29947 (2016). [PubMed]
11. Y. Gao and I. V. Shadrivov, “Second harmonic generation in graphene-coated nanowires,” Opt. Lett. 41(15), 3623–3626 (2016). [PubMed]
12. Y. Sun, Z. Zheng, J. Cheng, G. Sun, G. Qiao, and G. Qiao, “Highly efficient second harmonic generation in hyperbolic metamaterial slot waveguides with large phase matching tolerance,” Opt. Express 23(5), 6370–6378 (2015). [PubMed]
13. V. Y. Shur, A. R. Akhmatkhanov, and I. S. Baturin, “Micro- and nano-domain engineering in lithium niobate,” Appl. Phys. Rev. 2, 040604 (2015).
14. J. Sun and C. Xu, “466 mW green light generation using annealed proton-exchanged periodically poled MgO: LiNbO3 ridge waveguides,” Opt. Lett. 37(11), 2028–2030 (2012). [PubMed]
15. P. Koch, F. Ruebel, J. Bartschke, and J. A. L’huillier, “5.7 W cw single-frequency laser at 671 nm by single-pass second harmonic generation of a 17.2 W injection-locked 1342 nm Nd : YVO4 ring laser using periodically poled MgO : LiNbO3.,” Appl. Opt. 54(33), 9954–9959 (2015). [PubMed]
16. K. R. Parameswaran, J. R. Kurz, R. V. Roussev, and M. M. Fejer, “Observation of 99% pump depletion in single-pass second-harmonic generation in a periodically poled lithium niobate waveguide,” Opt. Lett. 27(1), 43–45 (2002). [PubMed]
17. D. Jedrzejczyk, R. Güther, K. Paschke, G. Erbert, and G. Tränkle, “Diode laser frequency doubling in a ppMgO:LN ridge waveguide: influence of structural imperfection, optical absorption and heat generation,” Appl. Phys. B 109, 33–42 (2012).
18. G. Li, Y. Cui, and J. Wang, “Photorefractive inhibition of second harmonic generation in periodically poled MgO doped LiNbO3 waveguide,” Opt. Express 21(19), 21790–21799 (2013). [PubMed]
19. A. Sahm, M. Uebernickel, K. Paschke, G. Erbert, and G. Tränkle, “Thermal optimization of second harmonic generation at high pump powers,” Opt. Express 19(23), 23029–23035 (2011). [PubMed]
20. T. Suhara and H. Nishihara, “Theoretical analysis of waveguide second-harmonic generation phase matched with uniform and chirped gratings,” IEEE J. Quantum Electron. 26(7), 1265–1276 (1990).
21. A. Tehranchi and R. Kashyap, “Design of novel unapodized and apodized step-chirped quasi-phase matched gratings for broadband frequency converters based on second-harmonic generation,” J. Lightwave Technol. 26(3), 343–349 (2008).
22. N. B. Baranova, M. A. Bolshtyanskii, and B. Y. Zel’dovich, “Adiabatic energy transfer from a pump wave to its second harmonic,” Quantum Electron. 25(7), 638 (1995).
23. O. Yaakobi, M. Clerici, L. Caspani, F. Vidal, and R. Morandotti, “Complete pump depletion by autoresonant second harmonic generation in a nonuniform medium,” J. Opt. Soc. Am. B 30(6), 1637–1642 (2013).
24. H. Suchowski, G. Porat, and A. Arie, “Adiabatic processes in frequency conversion,” Laser Photonics Rev. 8(3), 333–367 (2014).
25. C. R. Phillips and M. M. Fejer, “Efficiency and phase of optical parametric amplification in chirped quasi-phase-matched gratings,” Opt. Lett. 35(18), 3093–3095 (2010). [PubMed]
26. A. Leshem, G. Meshulam, G. Porat, and A. Arie, “Adiabatic second-harmonic generation,” Opt. Lett. 41(6), 1229–1232 (2016). [PubMed]
27. D. N. Nikogosyan, Nonlinear Optical Crystals: A Complete Survey (Springer, 2005).
28. J. Lin, Y. Xu, Z. Fang, M. Wang, J. Song, N. Wang, L. Qiao, W. Fang, and Y. Cheng, “Fabrication of high-Q lithium niobate microresonators using femtosecond laser micromachining,” Sci. Rep. 5, 8072 (2015). [PubMed]
29. F. Charra, G. Gurzadyan, and D. F. Nelson, Springer Materials-The Landolt-Börnstein Database (Springer-Verlag, 2000).
30. M. M. Fejer, G. A. Magel, D. H. Jundt, and R. L. Byer, “Quasi-phase-matched second harmonic generation: tuning and tolerances,” IEEE J. Quantum Electron. 28, 2631–2654 (1992).
31. O. Gayer, Z. Sacks, E. Galun, and A. Arie, “Temperature and wavelength dependent refractive index equations for MgO-doped congruent and stoichiometric LiNbO3,” Appl. Phys. B 91, 343–348 (2008).