Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Chirped-quasi-periodic structure for quasi-phase-matching

Open Access Open Access

Abstract

We propose in this paper a chirped-quasi-periodic structure using the projection method. This type of new structure combines the advantages of chirped and quasi-periodic structures, and can be used for both multiple quasi-phase-matching and multiple bandwidths control. Numerical simulation of second-harmonic generation performance is in good agreement with the Fourier spectrum of the structure.

©2010 Optical Society of America

1. Introduction

One of the challenges in the field of nonlinear frequency conversion is to design a single poled grating structure that can simultaneously quasi-phase-match several different χ(2) processes [16]. A χ(2) process involves the nonlinear mixing of two waves to produce a third wave at the sum or difference frequency. The process proceeds efficiently if the quasi-phase-matching (QPM) condition is satisfied. For example, in a sum frequency generation processω1+ω2=ω3, the QPM condition requires the wave vector mismatch between three wavesΔk=k3k2k1be compensated by a reciprocal lattice vectorG, i.e. G=k3k2k1.

If a single χ(2) process is desired, it is well known that for obtaining the largest Fourier coefficient the grating need to be periodically poled with periodΛ=2π/G. If several χ(2) processes are desired to be cascaded into one grating, the structure of the grating needs to be fabricated by some new means, such as quasi-periodic structures [1,2], aperiodic structures [3,4], dual-periodic structures [5], numerically optimized phase modulation [6] and so forth.

However, in cascaded χ(2) processes, for effective generation of every process the matching temperatures of all these processes need to be the same. For example, in a third-harmonic generation (THG) process [1], the third-harmonic (TH) wave is achieved by two steps—a second harmonic generation (SHG) cascading a sum frequency generation (SFG). Theoretically, using Sellmeier equation of the nonlinear material we can design a structure to quasi-phase-match these two processes perfectly. Practically, however, the QPM condition is sensitive. Tiny inaccuracy of Sellmeier equation or slight error of domain period will cause great shifts of the matching temperature of these processes. Therefore, the two bandwidths will not overlap well for efficient THG [7]. Moreover, for other cascaded processes such as an optical parametric oscillation (OPO) cascading a SFG or SHG [8,9], bandwidth matching of the two cascading processes is a crucial factor for generating desired visible light effectively.

As is known, chirped structure [1012] can be used to broaden acceptance bandwidth in nonlinear optical interactions. To solve the problem mentioned above, here we propose a novel structure for both multiple QPM and bandwidths controlling. The key point of our idea is to introduce chirping into quasi-periodic structure to design a chirped-quasi-periodic (CQP) structure. This idea provides more flexibility for structure designing of QPM gratings.

The paper is organized as follows: In Section 2, the designing approach of CQP structure is derived using projection method. A concrete CQP structure is analyzed in Fourier space in Section 3. Numerical simulations using SHG to verify the validity of the structure are given in Section 4. And in Section 5 is our conclusion.

2. Structure design

A one-dimensional quasi-periodic (QP) structure can be obtained by projection from a square lattice [13], in which the ‘hidden symmetry’ of the one-dimensional QP structure is displayed. See Fig. 1 . The projection line ξ has an angle θ to x axis. We label the orthogonal line to ξ as η axis and let the width of the projection window equal to the sum of the projection from unit spacing of x and y axis to η axis, i.e. w = sinθ + cosθ, as shown in Fig. 1. As long as the projection angle θ is irrational, the arrangement of the structure will be quasi-periodic.

 figure: Fig. 1

Fig. 1 Schema of the projection method to obtain a quasi-periodic structure.

Download Full Size | PDF

From the projection method described above, we can obtain a two-component quasi-periodic structure. The two components—block A and block B—are projected by vertical and horizontal spacing, respectively. For dielectric optical superlattice such as periodically-poled lithium niobate (PPLN) or lithium tantalite (PPLT), we set both block A and B consisting of a positive domain and a negative domain, and the lengths of the positive domains in each block have the same value l. The reciprocal vectors of this structure are given in Ref [2]:

Gm,n=2πm+nττDA+DB,
where the arrangement parameter τ is equal to the proportion of block number of A, NA, to block number of B, NB , i.e. τ = NA / NB. DA, DB are the lengths of block A and B, respectively.

When the quasi-periodic arrangement is projected from a square lattice, the proportion DA / DB is fixed, and we have τ = tanθ. However, in a more general situation, DA / DB is an adjustable parameter [2]. Under such circumstance the quasi-periodic arrangement needs to be projected from a rectangular lattice instead, and τ still indicates the proportion of block numbers, but the relation τ = tanθ do not hold any more.

Considering the situation of a rectangular lattice, let dx and dy indicate the horizontal and the vertical spacing, respectively. Then we can obtain the expression of τ:

τ=NANB=dxdytanθ.
Equation (2) indicates the arrangement parameter τ varies with the spacing of the projection lattice.

We present here a new grating designing approach which introduces the chirping factor into QP grating. Specifically, in the projection method described above, if we chirp the spacing of the original projection lattice first, then the influence of chirping will be hidden in the projected quasi-periodic structure, so that we obtain a chirped-quasi-periodic structure. To obtain a CQP structure we could chirp its horizontal and/or vertical spacing.

For a linear chirping, we can define a chirp factor r to describe the chirp rate, see Ref [3]:

rx(y)=dx(y)(N)dx(y)(1)dx0(y0),
where di(1), di(N)and di0 represent the length of the first, the last and the nominal spacing of the lattice, i = x,y indicates horizontal and vertical, respectively. The nth spacing’s length is:
dx(y)(n)=dx(y)(1)+nNrx(y)dx0(y0).
However, chirping is not limited to linear [11]. Indeed, as long as the domain period varies with position, linearly or non-linearly, the grating can be considered as chirped.

Here we present the general situation. We set both spacing vary with position arbitrarily, i.e. dx = dx(ξ) and dy = dy(ξ), ξ is the coordinate along the direction of the QP grating. From Fig. 1 we can see that using variable ξ to describe the spacing variation is equivalent to using coordinate x and y since they are proportional. Moreover, here we use continuous functions to describe discrete spacing, which is the same as Eq. (4) in Ref.12. Actually, for a given starting point, using a recursion method the desired structure can be easily obtained.

Generally, introducing of chirping will bring two kinds of influences to the QP structure. First, chirping the projection lattice will apparently affect the length of the blocks. Since all blocks A are projected from vertical spacing while all blocks B are projected from horizontal spacing, the block lengths DA and DB will vary with position ξ as following forms:

DA(ξ)=dy(ξ)sinθ,
DB(ξ)=dx(ξ)cosθ.

Secondly, chirping will also affect the arrangement parameter τ. From Eq. (2) we know that for a CQP structure τ can be expressed as:

τ(ξ)=dx(ξ)dy(ξ)tanθ.

Noticing that τ is a localized parameter in a CQP structure, different positions ξ have different quasi-periodic arrangements, which is totally different from the traditional QP structure. Substituting Eq. (5-7) into the expression of reciprocal vectors Eq. (1), we obtain:

Gm,n(ξ)=2π(mcosθdx(ξ)+nsinθdy(ξ))=2π(mcos2θDB(ξ)+nsin2θDA(ξ)).

Equation (8) can be treated as a more general expression of reciprocal vectors Gm,n of QP structure, which includes arbitrary spacing nonuniform. In Eq. (8) we use two sets of parameters (dx, dy, θ) and (DA, DB, θ) to describe the CQP structure. The former one shows the origin of the structure and is convenient for analyzing in projection method, and the latter one depends directly on the real block lengths DA and DB, which is therefore easier to draw structure. In the rest of this paper, we use the latter form to describe the structure. The deriving using the former form will be similar.

From Eq. (8) the designing approach of the multiple-bandwidth-controlling-structure can be obtained. Specifically, if we want two arbitrary reciprocal vectors Gm 1, n 1 and Gm 2, n 2 to stretch with position ξ as arbitrary form Gm 1, n 1(ξ) and Gm 2, n 2(ξ), then resolving following equations:

{Gm1,n1(ξ)=2π(m1cos2θDB(ξ)+n1sin2θDA(ξ)),Gm2,n2(ξ)=2π(m2cos2θDB(ξ)+n2sin2θDA(ξ)),
we obtain:

{DA(ξ)=2πsin2θ(m2n1m1n2)m2Gm1,n1(ξ)m1Gm2,n2(ξ),DB(ξ)=2πcos2θ(m2n1m1n2)n1Gm2,n2(ξ)n2Gm1,n1(ξ).

Equation (10) give the form of block lengths varying with position ξ when two desired reciprocal vectors are given. As shown in Eq. (10), the two blocks vary with position as a quite complex form, which in most case is non-linear chirping.

3. Fourier transformation of CQP structure

In this section, we design a concrete CQP structure and analyze it in Fourier space. The quasi-periodic structure we choose is projected from a rectangular lattice with horizontal spacing dx0 = 19.1μm and vertical spacing dy0 = 17.5μm to a line with slope tanθ = 0.414, the initial block lengths are DA 0 = 6.69μm and DB 0 = 17.64μm and the lengths of positive domain in both blocks are l = 6μm. For a 10-mm-long grating, the Fourier transformation of this QP grating is:

Some strong reciprocal vectors are marked in Fig. 2 . From Eq. (10) we know that two bandwidths can be controlled simultaneously in a CQP structure. Here we select two arbitrary reciprocal vectors, G 11and G 21, and let G 11 holds its δ–function shape while G 21 linearly stretches to ± δG, the two reciprocal vectors will vary with position ξ as following forms:

 figure: Fig. 2

Fig. 2 Fourier transformation of a quasi-periodic grating.

Download Full Size | PDF

{G1,1(ξ)=2π(cos2θDB0+sin2θDA0),G2,1(ξ)=2π(2cos2θDB0+sin2θDA0)+ξL/2L/2δG.

Here L is the length of the grating, ξ is the position coordinate of the grating, varies from 0 to L. Substituting (11) into (10), we obtain the relation of two block lengths vary with position:

{DA(ξ)=1/(1DA0δG(ξL/2)πLsin2θ),DB(ξ)=1/(1DB0+δG(ξL/2)πLcos2θ).

As expected, DA and DB both non-linearly vary with position ξ. From Eq. (12) we could obtain the structure parameters of the CQP grating. Here we choose δG = 0.01μm −1. The Fourier transformation of this structure is shown in Fig. 3 and the detail of G 21 is shown in inset. From Fig. 3 we can see that G 11 holds its original shape and G 21 is stretched to ± 0.01μm −1.

 figure: Fig. 3

Fig. 3 Fourier transformation of the CQP grating. The inset shows the detail of G 21.

Download Full Size | PDF

4. Second-harmonic generation in CQP structure

From analysis above we know that in a CQP grating some reciprocal vectors could hold the δ-function shape while the others were stretched. Thus we could expect these two kinds of reciprocal vectors will have different bandwidths in second-harmonic generation (SHG) process due to the bandwidth-broadening effect of chirped gratings [12].

We choose two reciprocal vectors of the CQP structure described in Section 3, G 21, which was stretched to ± 0.01μm −1, and G 11, which was not stretched, for numerical SHG. And we choose lithium tantalate (LT) as the nonlinear crystal and set the working temperature at 180°C. Using Sellmeier equation. of LT, we can estimate that G 11 and G 21 can quasi-phase-match SHG process around 1334.4 and 1103.9 nm, respectively. Under the slow-varying amplitude approximation, a SHG process should satisfy the following couple-wave equations [14] during propagation:

{dE1/dx=iω1d33f(x)n1cE2E1*exp(idkx),dE2/dx=i2ω2d33f(x)n2cE12exp(idkx),
where Ei, ωi, ni, indicate the electric field, angular frequency and refractive index of waves, and the subscript i = 1,2 represents the fundamental and second-harmonic wave, respectively. dk = k 2 -2k 1 represents the wave vector mismatch between fundamental and SH wave. Here we choose x, instead of ξ, to describe the propagating direction of waves, and f(x) satisfies:

f(x)={1whenxisinthepositivedomain,1whenxisinthenegativedomain.

The numerical results of SH conversion efficiency spectrums are shown in Fig. 4 . From the figure we can see, the full width at half maximum (FWHM) of the δ-function-shaped reciprocal vector G 11 is around 0.3nm while the FWHM of the linearly-stretched reciprocal vector G 21 is about 9nm, which is 30 times that of G 11. The two SHG bandwidths are in good agreement with the Fourier spectrum in Fig. 3.

 figure: Fig. 4

Fig. 4 SHG conversion efficiency versus fundamental wavelength using reciprocal vectors G 11 (a) and G 21 (b) in a CQP structure (The initial intensity of fundamental wave is 30MW/cm 2).

Download Full Size | PDF

Since the temperature mismatch and the wavelength mismatch can both be considered as a kind of phase-mismatch [12], based on the wavelength bandwidth data above we can obtain the temperature bandwidth of two reciprocal vectors. For the unstretched G 11 the 0.3nm wavelength FWHM equals 2.4°C temperature FWHM, and for the stretched G 21 the 9nm wavelength FWHM equals 134.3°C-222.3°C, which is 88°C temperature FWHM.

The above example illustrates two bandwidths can be controlled simultaneously in a CQP grating. In cascaded processes such as THG, by designing a structure with unstretched SHG reciprocal vector and properly stretched SFG reciprocal vector (10°C for example), the two bandwidths can be overlapped with a high tolerance of temperature-shift, thus the bandwidth mismatching effect mentioned in Chapter 1 can be solved.

5. Conclusion

We have proposed a novel structure named chirped-quasi-periodic structure for both multi-QPM and bandwidth control. We have analyzed the structure in Fourier space, and performed numerical simulations of SHG in CQP. Numerical results are in good agreement with theory. For further development, the CQP structure can be used for wide-tunable visible laser source based on cascading χ(2) processes, ultrashort multi-wavelength pulse-compression, et. al.

Acknowledgements:

This work was supported by the National Natural Science Foundations of China (NNSFC) under contracts 10904066 and 10874082 and the State Key Program for Basic Research of China contracts 2006CB921804 and 2010CB630703.

References and links

1. S. N. Zhu, Y. Y. Zhu, and N. B. Ming, “Quasi-phase-matched third-harmonic generation in a quasi-periodic optical superlattice,” Science 278(5339), 843–846 (1997). [CrossRef]  

2. C. Zhang, H. Wei, Y. Y. Zhu, H. T. Wang, S. N. Zhu, and N. B. Ming, “Third-harmonic generation in a general two-component quasi-periodic optical superlattice,” Opt. Lett. 26(12), 899–901 (2001). [CrossRef]  

3. B. Y. Gu, B. Z. Dong, Y. Zhang, and G. Z. Yang, “Enhanced harmonic generation in aperiodic optical superlattices,” Appl. Phys. Lett. 75(15), 2175–2177 (1999). [CrossRef]  

4. A. H. Norton and C. M. de Sterke, “Aperiodic 1-dimensional structures for quasi-phase matching,” Opt. Express 12(5), 841–846 (2004). [CrossRef]   [PubMed]  

5. Z. W. Liu, Y. Du, J. Liao, S. N. Zhu, Y. Y. Zhu, Y. Q. Qin, H. T. Wang, J. L. He, C. Zhang, and N. B. Ming, “Engineering of a dual-periodic optical superlattice used in a coupled optical parametric interaction,” J. Opt. Soc. Am. B 19(7), 1676–1684 (2002). [CrossRef]  

6. M. Asobe, O. Tadanaga, H. Miyazawa, Y. Nishida, and H. Suzuki, “Multiple quasi-phase-matched LiNbO3 wavelength converter with a continuously phase-modulated domain structure,” Opt. Lett. 28(7), 558–560 (2003). [CrossRef]   [PubMed]  

7. X. P. Hu, G. Zhao, C. Zhang, Z. D. Xie, J. L. He, and S. N. Zhu, “High-power, blue-light generation in a dual-structure, periodically poled, stoichiometric LiTaO3 crystal,” Appl. Phys. B 87(1), 91–94 (2007). [CrossRef]  

8. Z. D. Gao, S. N. Zhu, S.-Y. Tu, and A. H. Kung, “Monolithic red-green-blue laser light source based on cascaded wavelength conversion in periodically poled stoichiometric lithium tantalite,” Appl. Phys. Lett. 89(18), 181101 (2006). [CrossRef]  

9. G. K. Samanta and M. Ebrahim-Zadeh, “Continuous-wave, single-frequency, solid-state blue source for the 425-489 nm spectral range,” Opt. Lett. 33(11), 1228–1230 (2008). [CrossRef]   [PubMed]  

10. M. A. Arbore, O. Marco, and M. M. Fejer, “Pulse compression during second-harmonic generation in aperiodic quasi-phase-matching gratings,” Opt. Lett. 22(12), 865–867 (1997). [CrossRef]   [PubMed]  

11. G. Imeshev, “Tailoring of Ultrafast Frequency Conversion with Quasi-Phase-Matching Gratings,” Ph.D. dissertation (Stanford University, 2000).

12. X. J. Lv, Z. Sui, Z. D. Gao, M. Z. Li, Q. H. Deng, and S. N. Zhu, “Bandwidth and stability enhancement of optical parametric amplification using chirped ferroelectric superlattice,” Opt. Laser Technol. 40(1), 21–29 (2008). [CrossRef]  

13. R. K. P. Zia and W. J. Dallas, “A simple derivation of quasi-crystalline spectra,” J. Phys. Math. Gen. 18(7), L341–L345 (1985). [CrossRef]  

14. A. Yariv, and P. Yeh, Optical Waves in Crystal (Wiley, 1984).

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (4)

Fig. 1
Fig. 1 Schema of the projection method to obtain a quasi-periodic structure.
Fig. 2
Fig. 2 Fourier transformation of a quasi-periodic grating.
Fig. 3
Fig. 3 Fourier transformation of the CQP grating. The inset shows the detail of G 21.
Fig. 4
Fig. 4 SHG conversion efficiency versus fundamental wavelength using reciprocal vectors G 11 (a) and G 21 (b) in a CQP structure (The initial intensity of fundamental wave is 30MW/cm 2).

Equations (14)

Equations on this page are rendered with MathJax. Learn more.

G m , n = 2 π m + n τ τ D A + D B ,
τ = N A N B = d x d y tan θ .
r x ( y ) = d x ( y ) ( N ) d x ( y ) ( 1 ) d x 0 ( y 0 ) ,
d x ( y ) ( n ) = d x ( y ) ( 1 ) + n N r x ( y ) d x 0 ( y 0 ) .
D A ( ξ ) = d y ( ξ ) sin θ ,
D B ( ξ ) = d x ( ξ ) cos θ .
τ ( ξ ) = d x ( ξ ) d y ( ξ ) tan θ .
G m , n ( ξ ) = 2 π ( m cos θ d x ( ξ ) + n sin θ d y ( ξ ) ) = 2 π ( m cos 2 θ D B ( ξ ) + n sin 2 θ D A ( ξ ) ) .
{ G m 1 , n 1 ( ξ ) = 2 π ( m 1 cos 2 θ D B ( ξ ) + n 1 sin 2 θ D A ( ξ ) ) , G m 2 , n 2 ( ξ ) = 2 π ( m 2 cos 2 θ D B ( ξ ) + n 2 sin 2 θ D A ( ξ ) ) ,
{ D A ( ξ ) = 2 π sin 2 θ ( m 2 n 1 m 1 n 2 ) m 2 G m 1 , n 1 ( ξ ) m 1 G m 2 , n 2 ( ξ ) , D B ( ξ ) = 2 π cos 2 θ ( m 2 n 1 m 1 n 2 ) n 1 G m 2 , n 2 ( ξ ) n 2 G m 1 , n 1 ( ξ ) .
{ G 1 , 1 ( ξ ) = 2 π ( cos 2 θ D B 0 + sin 2 θ D A 0 ) , G 2 , 1 ( ξ ) = 2 π ( 2 cos 2 θ D B 0 + sin 2 θ D A 0 ) + ξ L / 2 L / 2 δ G .
{ D A ( ξ ) = 1 / ( 1 D A 0 δ G ( ξ L / 2 ) π L sin 2 θ ) , D B ( ξ ) = 1 / ( 1 D B 0 + δ G ( ξ L / 2 ) π L cos 2 θ ) .
{ d E 1 / d x = i ω 1 d 33 f ( x ) n 1 c E 2 E 1 * exp ( i d k x ) , d E 2 / d x = i 2 ω 2 d 33 f ( x ) n 2 c E 1 2 exp ( i d k x ) ,
f ( x ) = { 1 w h e n x i s i n t h e p o s i t i v e d o m a i n , 1 w h e n x i s i n t h e n e g a t i v e d o m a i n .
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.