Reach Us
+44-1522-440391

**Huo G ^{1*}, Chen H^{1} and Zhang M^{2}**

^{1}College of Control Engineering, Xijing University, Xi’an, China

^{2}School of Electronics Engineering, Xi’an University of Posts and Telecommunications, Xi’an, China

- *Corresponding Author:
- Huo G

College of Control Engineering

Xijing University, Xi’an 710123, China

**Tel:**(8629) 85628198

**E-mail:**[email protected]

**Received Date**: November 11, 2016; **Accepted Date:** December 13, 2016; **Published Date**: December 30, 2016

**Citation: **Huo G, Chen H, Zhang M (2016) A Generalized Method for Calculating Phase Matching Conditions in Biaxial Crystals. J Laser Opt Photonics 3: 143. doi: 10.4172/2469-410X.1000143

**Copyright:** © 2016 Huo G, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Visit for more related articles at** Journal of Lasers, Optics & Photonics

We present a generalized method for calculating phase matching conditions in biaxial crystals, especially in nonlinear monoclinic and triclinic crystals. Exploiting the angle definition introduced by Japanese mathematician Kodaira Kunihiko, we deduce the angular relations in geometry and obtain the expressions of refractive indices depending on angular orientation of wave vector k and optical axis angle. Then we calculate the phase matching parameters with bismuth borate BiB3O6 (BIBO) crystal in spontaneous parametric down conversion (SPDC) process for the type I and type II. On its basis, we discuss the angular gradients of the pump and emission wave refractive index near the exact phase matching direction, and compare the SPDC with double frequency process in geometrical relations of the refractive index ellipsoids. It indicates the anisotropic structures of nonlinear crystals making a capability of filtering effect, and the nonzero linear mismatch described by refractive index angular gradients can be used to estimate the beam width in crystal. This method is convenient to calculate the phase matching parameters in orthorhombic crystals without solving the quadratic Fresnel equations.

Nonlinear crystals; Phase matching; Angle-dependent refractive index; Spatial and time walk-off

Since the second-harmonic generation is first found in 1961 [1], other parameter processes were studied following it, such as sum frequency, difference frequency, optical parametric amplifiers, and optical parametric oscillator. Now, nonlinear frequency conversion has made rapid progress in laser sources owing to its simplicity and wavelength flexibility [2-4]. In 1995, the high efficiency polarization entanglement generation is first demonstrated with BBO crystal [5]. Spontaneous parametric down conversion (SPDC) in nonlinear crystals becomes the mainstay for entanglement production in quantum optics [6,7]. In the process of frequency conversion, phase matching is a key element in improving the conversion efficiency. The imperfect phase matching originates from wave vector direction not coinciding with that of pointing vector inside nonlinear crystal, depressing the conversion efficiency. It results from the detuning of central wavelength and modulation of pump pulse. It is also affected by the energy of pump and emission beams propagating along any phase matching direction with slightly deviation. We need take into account the dielectric dispersion and temperature to calculate accurately the phase matching. Comparing with uniaxial crystal, biaxial crystal, such as Mica, Sapphire, Olivine, Sulfur, has two axes along which the two linearly polarized optical waves have same propagation velocity (or refractive index) [8-11]. Under common conditions in biaxial crystals, two extraordinary linearly polarized beams exist for an incident light, and the refractive index in space distribution are two different ellipsoids [12]. Biaxial nonlinear crystal has excellent optical properties, especially for relatively large effective nonlinear coefficient. Taking the BiB_{3}O_{6} (BIBO) as an example, it has wide range of optical transparency, high damage threshold, inertness to moisture and the effective nonlinear coefficient attaining 2.5 pm/V [13,14]. These properties make biaxial crystal versatile in phase matching [15-25]. Phase matching in biaxial crystals is first calculated in 1984 [26], and the improved methods are worked out in the following years [27-30]. All the work was done by solving quadratic Fresnel equations. As the complex optical structure of crystals, a convenient approach for calculating phase matching conditions is worth exploring.

Here, we give out a convenient generalized method for calculating the phase matching parameters of nonlinear frequency conversion in nonlinear crystals, especially for biaxial crystal. In virtue of exploiting the angle definition introduced by Japanese mathematician Kodaira Kunihiko, we obtain the refractive indices contain orientation angles (θ, φ). It can be directly applied in phase matching conditions. Taking SPDC with BIBO crystal as an example, we calculate the effective nonlinear coefficients for type I and type II phase matching. On its basis, we discuss the angular gradients of the pump and emission wave refractive index near the exact phase matching direction, and compare the SPDC with double frequency process in geometrical relations of the refractive index ellipsoids. It indicates the anisotropic structures of nonlinear crystals making a capability of filtering effect, and the nonzero linear mismatch described by refractive index angular gradients can be used to estimate the beam width in crystal. This method is convenient to calculate the phase matching parameters in orthorhombic crystals without solving the quadratic Fresnel equations.

In biaxial crystals, we can calculate refractive index by wave normal Fresnel equations in principal coordinates system (xyz), while the crystal parameters in fact are measured in crystallographic system (abc). These axes in two different coordinates system are not parallel to each other. In order to find the angle dependence of refractive index, we rotate the principal coordinates system (xyz) φ around z axis, and acquire coordinates x′y′z′; then we rotate the coordinates (x′y′z′) θ around y′axis, obtianing coordinates x"y"z", where θ is polar angle relative to the optical z-axis in plane x′z′, and φ is azimuthal angle measured from x axis in xy plane. The schematic in **Figure 1** shows the coordinate transformation of dielectric axis (xyz), crystal principal axis (abc) and lab frame (x"y"z" )The two optical axes (C_{1}, C_{2}) distribute symmetrically about z axis, and the angles between z and C_{1}, C_{2} are Ω, −Ω respectively. Moreover θ_{1},θ_{2} are angles between wave vector k and C_{1}, C_{2}, and, θ ′,θ " are angles between z axis and the projections of k in yz and xz planes.

Proved by wave normal Fresnel equations [31], the refractive indexes of two extraordinary polarized beamsin n_{1,2} biaxial crystal can be calculated by

(1)

Where nx ny are the indexes in x and z directions. We can get n_{1} (n_{2}) by taking the upper (lower) sign “+” (“-”) in right side. The angles θ_{1} and θ_{2} are key parameters for calculating the refractive index. What’s more, the refractive index exerts an influence on the phase matching conditions in nonlinear frequency conversion, and its derivative decides the birefringence anisotropy in quantum entangled wave packet distribution.

For nonlinear crystals, it is convenient to calculate phase matching parameters in lab frame (x"y"z") during frequency conversion process, while the refractive index is solved easier in crystal’s dielectric axes (xyz). Then, we can seek angular conversion relations between the two coordinates systems in order to utilize the crystal parameters at will. With the method of angle projecting and the angle definition introduced by Japanese mathematician Kodaira Kunihiko [32], the θ_{1}, θ_{2} can be expressed below [19]:

(2a)

(2b)

In Eq (2), θ "=*arctan(tanθ. cosφ)*, and θ ,φ meet the conditions of 0 ≤θ < 90, 0 ≤φ ≤ 90. The optical axis angle can be expressed by where n_{x}, n_{y} and n_{z} are principal refractive indices at a given wavelength [20]. Nevertheless, the Eqs (2) remain correct in the range of and 90 < θ ≤180 and 0° ≤φ ≤ 90° , replacing with θ "=−arctan(tanθ ⋅cosφ ).

From Eqs (1 and 2), the refractive indices are function of Ω ,θ , φ ultimately in monolithic and orthorhombic crystal. The optical axis angle Ω is influenced by wavelength, temperature and so on. The refractive indices here contain two space angles (θ_{1}, θ_{2}). We can take advantage of it to study the characters of phase matching conditions and refractive index angular distribution. The above results are based on monolithic crystal, yet this method can be extended to triclinic crystal through a more angular projecting process, shown in **Figure 2**. The three crystal axes a, b, c are not vertical to each other in triclinic crystal. We can make b axis projected similarly, decomposed it into two segments which are vertical and parallel to the plane constituted of b and c, separately. According to the cosine relationship cosβ=cosβ ⋅cosβ the angles β_{1}, β_{2} can be expressed with angles Ω ,θ and φ. We can further obtain the refractive indices in triclinic crystal system through Eq (1 and 2).

In nonlinear optics, phase matching conditions are crucial to threewave frequency mixing process. It is a complex problem especially in biaxial crystal. At present, the calculation of phase matching in biaxial crystal needs to solve quadratic Fresnel equations. Here we calculate the phase matching conditions with a direct method of angular dependence of refractive index.

In three wave mixing process, it obeys the conservation laws of energy and momentum under phase matching conditions, which meets the relations below

(3a)

(3b)

Where ω_{p} ω_{s} ω_{i} the frequency of pump, signal, and idler beams, respectively. The k_{p}, k_{s}, k_{i} are corresponding wave vectors of the three waves. Generally, there are two beams in biaxial crystal. One is called the ‘‘fast’’ with a smaller refractive index, and the other is named as the ‘‘slow’’ with a bigger one. The momentum conservation in Eq (3) has following forms:

(4a)

(4b)

Eq (4a) is the type I phase matching configuration, and Eq (4b) is the type II.

The wave vectors can be expressed as (l=p, s, i) for convenience, where n_{1} is the refractive index in phase matching direction, k_{1} is the unit vector, and θ_{l}, φ_{l} are polar and azimuthal angles. The Eq (4) become

(5a)

(5b)

In Eq (5), the refractive index n_{l, m}(*l=p, s, i; m=slow, fast*) relies on the phase matching angles (θ_{l},φ_{l}) **Figure 3** describes the schematic of three wave mixing nonlinear process in biaxial crystal.

For the collinear degenerate SPDC process, the three unit vectors are in the same direction, and one incident photon is split into a pair of lower energy correlated photons. The two down-converted photons have equal frequency, satisfying the relation ω_{i}=ω_{s}=ω_{p}/2. Then the Eqs (5) are simplified in refractive indices forms [16].

(6a)

(6b)

Where n_{l,m} represent the refractive indices of fundamental and SPDC waves.

The refractive index in Eq (6) at any phase matching direction can be calculated by Eq (1). It is written by:

(7a)

(7b)

where l=p, s, i, and θ_{1}, θ_{2} are deduced in Eq (2) containing the phase matching angles (θ ,φ ). Comparing the two refractive indices in Eqs. (7), the smaller (bigger) one is the “fast” (“slow”) indicated as n_{lf} (n_{lf}).

BIBO is a new type biaxial crystal with good optical properties. Using the parameters from Hellwig’s data at room temperature [14], the dispersion curves are plotted in **Figure 4**. From the dispersion relation in **Figure 4**, we can get the Sellmeier coefficients at arbitrary wavelength.

Combining Eq (6) and Eq (7) together, we can get the relations of phase matching angles (θ ,φ ) in SPDC process in the range of θ = 0° − 90° and θ = 0° − 90°. **Figure 5a** shows the type I phase matching angles at fundamental wavelength of 400nm, while **Figure 5b** describes the type II at 800nm. It is fundamental to calculate the phase matching characters, such as the effective nonlinear coefficients, the spatial walk-off angle, acceptance angle, spectral acceptance, and group velocity dispersion.

From the discussion above, we solve the linear relations of refractive indices in phase matching conditions, avoiding complex quadratic Fresnel equations. It is a convenient method for calculating the phase matching parameters, especially for biaxial crystal. We can see that the spatial walk-off angles are bigger for type II phase matching, depressing the efficiency of SPDC. Thus, the angle acceptance and spectral acceptance are larger than that of type I, which is good for improving the efficiency. These results indicate that the type II phase matching in SPDC process with BIBO crystal is a good choice. Owing to afore mentioned monoclinic symmetry, the distributions of the effective nonlinear coefficient are symmetrical around the axes φ=0, φ=90 and φ=−90 [33]. We can obtain the optimal phase matching conditions easily in other principle planes.

The specified phase matching parameters here are acquired by taking advantage of angular-dependent refractive indices of biaxial crystals in Eq (7). In the near axis and wide crystal approximation, the function n(θ ,φ) can be approximated by the zero and first orders of its expansion. The linear detuning is

(8)

Where

(9)

is the gradient of refractive index at (θ_{0},φ_{0}) and the linear detuning Δ does not equal zero.

The angular gradient of pump wave refractive index n′_{p} can be used to explain three wave interactions both in and out of principal planes. It is easy to ascertain the directions of pump and emission waves by two refractive index ellipsoids. More importantly, the linear detuning offered a method for describing the beam width in biaxial crystals, affording a quantification parameter to estimate the space distribution of pump and emission waves in quantum optics.

In conclusion, we have presented an effective method for calculating phase matching with biaxial crystals in three-wave frequency mixing. The refractive indices depending on angular orientation (θ,φ) can be directly applied in common phase matching conditions without solving quadratic Fresnel equations. Taking collinear SPDC in BIBO for an example, we calculated the phase matching angles for the type I and type II based on angle-dependent refractive index. It is fundamental for calculating all the other phase matching parameters, such as the effective nonlinear coefficients, the spatial and time walkoff. Moreover, the angular gradient of refractive index near the phase matching direction can be used to decide the optimum phase matching conditions and the space distribution of pump and emission waves. This method is convenient to calculate the phase matching parameters in orthorhombic crystals, and it can be used to estimate the beam width in crystal.

The authors thank the Xijing University Research Foundation for Talented Scholars (No. XJ15B02), and the Scientific Research Program Funded by Shaanxi Provincial Education Department (No. 16JK2247).

- Franken PA, Hill AE, Peters CW, Weinreich G (1961) Generation of optical harmonics. Phys Rev Lett 7: 118-119.
- Herault E, Balembois F, Georges P, Georges T (2008) 1064 nm Nd: YVO4 laser intracavity pumped at 912 nm and sum-frequency mixing for an emission at 491 nm. Opt Lett 33: 1632-1634.
- Zou S, Li P, Wang L, Chen M, Li G (2009) 980 nm Yb-doped single-mode fiber laser and its frequency-doubling with BIBO. Appl Phys B 95: 685-690.
- Petrov V, Ghotbi M, Kokabee O, Esteban-Martin A, Noack F, et al. (2010) Femtosecond nonlinear frequency conversion based on BiB3O6. Laser & Photon Rev 4: 53-98.
- Kwiat PG, Mattle K, Weinfurter H, Zeilinger A, Sergienko AV, et al. (1995) New High-Intensity Source of Polarization-Entangled Photon Pairs Phys Rev Lett 75: 4337-4341.
- Rangarajan R, Goggin M, Kwiat P (2009) Optimizing type-I polarization-entangled photons. Opt Express 17: 18920-18933.
- Halevy A, Megidish E, Dovrat L, Eisenberg HS, Becker P, et al. (2011) The biaxial nonlinear crystal BiB3O6 as a polarization entangled photon source using non-collinear type-II parametric down-conversion. Opt Express 19: 20420-20434.
- Wyncke B, Brehat F (1989) Calculation of the effective second-order non-linear coeffcients along the phase matching directions in acentric orthorhombic biaxial crystals. J Phys B At Mol Opt Phys 22: 363.
- Roberts DA (1992) Simplified characterization of uniaxial and biaxial nonlinear optical crystals: a plea for standardization of nomenclature and conventions. IEEE J Quantum Electron 28: 2057-2074.
- Ito H, Naito H, Inaba H (1975) Generalized study on angular dependence of induced second-order nonlinear optical polarizations and phase matching in biaxial crystals. J Appl Phys 46: 3992-3998.
- Grechin SG, Grechin SS, and Dmitriev VG (2000) Complete classification of interaction types for the second-harmonic generation in biaxial nonlinear crystals. Quantum Electronics 30: 377-386.
- Dmitriev VG, Gurzadyan GG, Nikogosyan DN, Handbook of Nonlinear Optical Crystals, 3rd revised.
- Hellwig H, Liebertz J, Bohaty L (1999) Exceptional large nonlinear optical coefficients in the monoclinic bismuth borate BiB3O6 (BIBO). Solid State Commun 109: 249-251.
- Hellwig H, Liebertz J, Bohaty L (2000) Linear optical properties of the monoclinic bismuth borate BiB3O6. J Appl Phys 88: 240-244.
- Huo G, Zhang T, Wan R, Cheng G, Zhao W (2013) The Factors of Spontaneous Parametric Down Conversion Process in BiB3O6 Crystal with a Broadband Pump. Optik 124: 6627-6630.
- Huo G, Wan RG, Zhang MZ, Cheng GH, Zhao W (2012) Optimizing the Phase-Matching Parameters of a BIBO Crystal for Ultrafast Spontaneous Parametric Down Conversion. Proc of SPIE 8333: 833313.
- Huo G, Zhang T, Cheng G, Zhao W (2013) Calculation of effective nonlinear coeffective in BIBO for spontaneous parametric down conversion. J Nonlinear Opt Phys 22: 1350010.
- Zhang M, Huo G (2014) Calculation of phase matching conditions based on angle-dependent refractive index Opt Eng 53: 086105.
- Huo G, Wang Y, Zhang M (2015) An effective method for calculating phase matching conditions in biaxial crystals. Appl Phys B 120: 239-246.
- Du C, Wang Z, Xu X, Teng B, Fu K, et al. (2001) Efficient intracavity second harmonic generation at 1.06 μm in a BiB3O6 (BIBO) crystal. Appl Phys B 73: 215-217.
- Wang Z, Teng B, Fu K, Xu X, Song R (2002) Efficient second harmonic generation of pulsed laser radiation in BiB3O6 (BIBO) crystal with different phase matching directions. Opt Commun 202: 217-220.
- Kityk IV, Imiolek W, Majchrowski A, Michalski E (2003) Photo-induced second harmonic generation in partially crystallized BiB3O6 glass. Opt Commun 219: 421-426.
- Du C, Teng B, Wang Z, Liu J, Xu X, et al. (2002) Actively Q-switched intra-cavity second-harmonic generation of 1.06 μm in BiB3O6 crystal. Optics & Laser Technology 34: 343-346.
- Czeranowsky C, Heumann E, Huber G (2003) All-solid-state continuous-wave frequency-doubled Nd:YAG-BiBO laser with 2.8 W output power at 473 nm. Opt Lett 28: 432-434.
- Ghotbi M, Ebrahim-Zadeh M (2004) Optical second harmonic generation properties of BiB3O6. Opt Express12: 6002-6019.
- Yao JQ, Fahlen TS (1984) Calculations of optimum phase match parameters for the biaxial crystal KTiOPO4. J Appl Phys 55: 65-68.
- Kaschke M, Koch C (1989) Calculation of nonlinear optical polarization and phase matching in biaxial crystal. Appl Phys B 49: 419-423.
- Dreger MA, Erkkila JH (1992) Improved method for calculating phase-matching criteria in biaxial nonlinear materials. Opt Lett 17: 787-788.
- Yao JQ, Sheng WD, Shi WQ (1992) Accurate calculation of the optimum phase-matching parameters in three-wave interactions with biaxial nonlinear-optical crystals. J Opt Soc Am B 9: 891-902.
- Zadorozhnii VI (2000) Improved analytical method for calculating the parameters of phase-matching nonlinear-optical interactions in biaxial crystal. Opt Commun 176: 489-501.
- Born M, Wolf E (1970) Principles of Optics Cambridge University Press, Cambridge, UK.
- Kodaira K (2003) An Introduction to Calculus (Originally published in Japanese by Iwanami Shoten, Publishers, Tokyo.
- Huo G, Zhang M (2015) Monolithic symmetry effect on optimizing the phase matching parameters of a biaxial crystal for spontaneous parametric down conversion. J Optoelectrion Adv M 17: 116-121

Select your language of interest to view the total content in your interested language

- Capacitor
- Circuits
- Coal Consumption
- Contingency Analysis
- Current
- Directional Microphone
- Dynamic Stabilizer
- Electrical Engineering
- Electricity
- Electromagnetism
- Femtosecond Laser
- Fibre Optical Sensors
- Infrared
- Integrated Circuit
- Laser Diode
- Laser Spectrum
- Microphone Array
- Nano Photonics
- Nonlinear Optics
- Omnidirectional Microphone
- Optical Fiber
- Optical Switch
- Optics
- Oscillations
- Oscillator
- Photon
- Photonic Band Gap
- Photonic Crystal
- Photonic Crystal Fiber
- Photonics
- Power Analysis
- Power System Stabilizers
- Quantum Dot
- Quantum Well
- Quaternary Logics
- Refractive Index
- Robot Audition
- Semiconductor
- Semiconductor Laser
- Sensor
- Signal Crystal
- Single Mode
- Spectrum
- Two Photonabsorption
- Wave Theory

- Total views:
**2308** - [From(publication date):

December-2016 - Dec 11, 2019] - Breakdown by view type
- HTML page views :
**2136** - PDF downloads :
**172**

**Make the best use of Scientific Research and information from our 700 + peer reviewed, Open Access Journals**

International Conferences 2019-20