Search tips
Search criteria 


Logo of scirepAboutEditorial BoardFor AuthorsScientific Reports
Sci Rep. 2017; 7: 6579.
Published online 2017 July 26. doi:  10.1038/s41598-017-06771-2
PMCID: PMC5529533

Analytical investigation of nonreciprocal response in 1D nonlinear photonic crystals


The nonreciprocal response of the SHG process in 1D periodical nonlinear photonic crystals with a defect embedded has been theoretically studied by solving the nonlinear coupled equations. The nonreciprocal response has been deduced analytically with the solution of non-reciprocity parameters obtained. The result shows that as the non-reciprocity approaches 100%, the crystal length and the input power needed increase at a logarithmic rate. Any target nonreciprocal response can be reached in this structure by adjusting the structure parameters.


Recent years have witnessed a growing tendency to study nonreciprocal structures and devices in many different areas13. Inspired by the development of photoelectric diode4, 5, researchers have successively invented the acoustic diode6, 7 and the thermal diode8, 9. According to Lorentz reciprocity theorem10, the light propagation process is usually reversible. Therefore, it is difficult to achieve non-reciprocity in optical field. However, several effective ways have been figured out to overcome this problem in recent years. For example, a metal-silicon waveguide system was proposed in 2011, which enabled light to propagate unidirectionally on silicon chips by adjusting light potential energy11. In nonlinear optics, an all-optical diode made of a defect embedded in lithium niobate channel waveguide was numerically demonstrated with a spatially nonreciprocal response12, 13. External electric field can be employed to further improve the isolating performance for this structure1416, and multiple quasi-phase-matching (QPM) technique was also introduced to realize all-optical isolating action17.

In these researches, numerical computation is the main approach used to study the nonreciprocal process, which may not directly explain why it shows a nonreciprocal response as well as how to get a bigger response. In this paper, we have studied the above 1D defective structure theoretically and obtained a series of analytical solutions not only for the second-harmonic output but also for the nonreciprocal response. Further derivations have been implemented and some new results have been revealed from the solutions. The role of the non-reciprocity parameter P as well as its corresponding realization conditions has been discussed in detail.

SHG output solutions and simulations for forward and backward processes

It is well known that the field distributions of the SHG process in homogenous nonlinear crystals can be solved exactly with the help of Jacobi elliptic functions18. In this paper we show that the nonreciprocal SHG process in 1D defective nonlinear photonics crystal (NPC) structures can be solved in a similar manner. The schematic of the defective NPC structure we studied is shown in Fig. 1, where the up and down arrows represent the domain poling directions. Supposing the whole length of the NPC is L, a defect with a width of δL is embedded at the position x = L 1. Theoretically, the whole structure can be treated as three sections, where the first and third sections have the same period (denoted as Λ) but different initial conditions and the second section is the embedded defect inducing phase shift between the fundamental wave (FW) and second-harmonic wave (SHW). The nonreciprocal responses of the structure can be revealed by studying the forward and backward SHG processes separately, which can be treated analytically with the nonlinear coupled equations.

Figure 1
A schematic of the forward and backward SHG processes in a 1D nonlinear photonic crystal with a defect embedded.

Firstly, we focus on the forward SHG process. The nonlinear coupled equations for this process are usually expressed as19:


where A 1, A 2 refer to the field distributions of FW and SHW, respectively. K is the nonlinear coupling coefficient and f(x) represents the structure function of the NPC. Δk = k 2  2k 1 describes the phase mismatching in the SHG process, where k 1 and k 2 are respectively the wave vectors of FW and SHW.

For the first section of the NPC (0 ~ L 1), when the QPM condition is satisfied2022, Eq. (1) can be simplified as:


where K 1 = K · g m, and g m is the Fourier coefficient of the reciprocal vector used in the QPM process. Eq. (2) can be solved easily with the initial conditions of A 1(0) = A 0, A 2(0) = 0, and the following solutions for A 1 and A 2 at the position x = L 1 can be obtained17:


For the second section (L 1 ~ L 1 + δL), a phase shift (denoted as Δ[var phi]) is induced between the FW and SHW, which can be written as Δ[var phi] = (2δL  1)π.

For the third section (L 1 + δL ~ L), the simplified coupled equations under the QPM condition are similar with Eq. (2) except additional terms exp(iΔ[var phi]) and exp(−iΔ[var phi]) resulting from the second section, which can be expressed as:


In this situation, the initial conditions can be determined by Eq. (3), where not only the FW but also the SHW has a non-zero initial value.

If we concentrate on the variation of amplitudes and phases respectively, Eq. (4) can be further reduced to the following ones:



where y 1, y 2 represent the amplitudes of the FW and SHW respectively while [var phi] 1, [var phi] 2 represent the corresponding phases and θ = 2[var phi] 1  [var phi] 2  Δ[var phi].

From the above equations and QPM conditions, two first integrals can be drawn. They are


where Γ is a constant and its value can be obtained by substituting Eq. (3) to the above equation:


On the basis of Eqs (6) and (5a) can be further solved as:


Supposing that the term under the square root sign can be written as (A022y22)2y22Γ2=4(y22A) (y22B)(y22C), the parameters A, B, C can be expressed as23, 24:


where ζ=arccos(1+27Γ2/A06)/3 and A  B  C. Here A, B, C are all real numbers since 027Γ2/A062 is always satisfied and thus obviously ζ [set membership] [0, π/3]. Eq. (8) can be further solved as:

y22(L) = a2sn2[bK1(L − L1) ± F(γφ), γ] + A

where F(γ,φ)=F(a/b,arcsin(y22(L1)A/a)) is the Legendre’s incomplete elliptic integral of the first kind, and a, b are two real constants which are defined by a=BA and b=CA. sn(z,γ) represents Jacobi elliptic function with modulus γ = a/b. In Eq. (10), we should choose a “+” sign if y 2 is to increase with increasing x and a “−” sign if y 2 is to decrease with increasing x.

Secondly, for the backward transmission, the results are similar to the forward ones, thus the amplitude for SHW can be expressed as:

y22(0) = a′2sn2[bK1L1 ± F(γ′, φ′), γ′] + A

where F(γ,φ)=F(a/b,arcsin(y22(LL1)A/a)). Parameters a′, b′, A′ can be obtained by repeating the above deduction processes with the defect position L  L 1 instead of L 1.

Numerical calculations have been conducted together with the above analytical solutions to demonstrate the nonreciprocal response of the defective NPC structure. Assuming that the defect position is L 1 = 0.2L and δL = 5Λ/6, that is the phase shift Δ[var phi] equals to 2π/3, the normalized amplitudes of FW and SHW in the NPC are shown in Fig. 2. Figure 2(a) and (b) show the forward and backward transmission processes calculated on the basis of the original coupled equations Eq. (1). Figure 2(c) and (d) are corresponding processes obtained by numerical calculation of the simplified equations Eqs (2) and (4). Whereas Fig. 2(e) and (f) show the corresponding results of the analytical solutions Eqs (10) and (11). We can see that they match well with each other.

Figure 2
Normalized amplitudes of FW and SHW for the forward and backward processes in the defective NPC. (af) are the corresponding forward and backward propagations calculated on the basis of the original coupled equations, the simplified coupled equations ...

Analytical and numerical derivations of the non-reciprocity parameter

According to [12–15], the non-reciprocity parameter P in this SHG process can be defined as:


Here, superscripts “+” and “−” represent the forward and backward processes, respectively. In the situation that L  L 1 [dbl greater-than sign] L 1, Eq. (12) can be simplified as:


where (η+)2=2|y2+|2/A02 describes the SHG conversion efficiency of the forward process.

Substituting Eq. (10) into Eq. (13), it is found that when the forward SHG output reaches its minimum, the non-reciprocity achieves the maximum, and vice versa. According to the properties of Jacobi elliptic function, when the forward output achieves the minimum or maximum, the following expressions for the total length of the crystal L can be obtained respectively:


where T(γ) = F(γ, π/2) describes the period of the first kind of complete elliptic integral and n represents positive integers. When n is an even number, the forward SHG output reaches the minimum value A and otherwise achieves the maximum value B. The corresponding extremes of P can be solved as


It can be analytically solved from Eq. (15) that P max equals to P min when the condition 27Γ2/A06=2 is satisfied. That means with specific values of A 0 K 1 L 1 and Δ[var phi], no matter what the length of the total crystal, the value of P keeps the same.

Surface plot of the extremes of is calculated and shown in Fig. 3. It can be seen that high nonreciprocal responses can be obtained when Δ[var phi] is close to 0 or π, which corresponds to the situation that Γ is small. Furthermore, if Γ  0, ideal non-reciprocity can be obtained. The expression of P max can be approximately simplified as follows:


Figure 3
Surface plot of the extremes of the non-reciprocity parameter with different structure parameters.

In this situation, γ can be expressed as γ=132Γ/A03. Substituting Eq. (16) into Eq. (14), it can be concluded that the total NPC length needed increases logarithmically with P max since the corresponding period T(γ) can be approximately expressed as


Figure 4 shows the dependence of P max and Δ[var phi], where Δ[var phi] varies in the range of 0 and π, and A 0 K 1 L 1 is set to be 0.3π. Figure 4(b) is the partial magnification of Fig. 4(a). It can be seen that the non-reciprocity parameter P is periodically distributed for a given Δ[var phi], which is resulted from the periodicity of the Jacobi elliptic function. Besides, when Δ[var phi]  0 or Δ[var phi]  π, the corresponding period will increase logarithmically, which matches well with the analytical discussions.

Figure 4
Distribution of non-reciprocity parameter with different values of Δ[var phi] under the condition that A 0 K 1 L 1 = 0.3π. Different color represents different nonreciprocal responses.

Simulation results show that the contrast ratio of the all-optical diode can maintain close to 1. Thus, a near complete nonreciprocal all-optical diode based on a NPC can be realized. As is shown in the Fig. 4(b), the same color describes the same nonreciprocal response. It can be seen that for any target non-reciprocity, a NPC structure with a defect embedded can be designed with specific values of Δ[var phi] and A 0 K 1 L 1, which are realized by adjusting the defect width, the crystal length and the FW input. As Δ[var phi] moves away from zero, the same nonreciprocal response is kept as long as the total length of the crystal decreases. It means that the target non-reciprocity could be achieved with a shorter crystal length and a lower FW input by adjusting the defect width. For practical applications in a defective NPC structured on the lithium niobate (LN), the nonlinear coefficient d 33 = 27 pm/V. If a 70% nonreciprocal contrast is required, that is P max = 0.7, A 0 K 1 L  6.191 and Δ[var phi]  0.532 can be carried out according to Eqs (7) and (14), sketched by point ‘A’ in Fig. 4(b). The data show that if we set the FW (1064 nm) intensity to be 60 MW/cm 2, the needed total length is 1.3 cm and the required defect width is 3.8 μm with the QPM period being 6.6 μm. Similarly, if an 80% nonreciprocal contrast is required, the corresponding parameters can be obtained by the same method, shown by point ‘B’ in the figure. Discussions above have explained a universal design of nonlinear crystals. In comparison with previously reported NPC diodes, our results reveal that the design of diodes can be more flexible and easier to control.


In all, we have theoretically studied the non-reciprocity properties of SHG process realized by a defective 1D NPC. Exact solutions for the SHG output in both forward and backward transmission processes have been derived with elliptic function forms. According to the property of the elliptic function, it can be concluded that for specific values of Δ[var phi] and A 0 K 1 L 1, the non-reciprocity parameter has a maximum value as well as a minimum value. It is found that the complete non-reciprocity cannot be realized in this structure but can be approached infinitely. In this case, the crystal length needed increases with the nonreciprocal parameter at a logarithmic rate. Any target non-reciprocity can be obtained by adjusting the defect width, the crystal length and the FW input. These results may lead to a variety of relevant applications, including all-optical diode, optical isolator, amplifiers and so on.


It is obvious that there is no non-reciprocity when the defect is at the middle position, which means L 1 = L/2. Asymmetric structures should be adopted for significant responses. Hereafter we present a detailed analysis under the condition that L  L 1 [dbl greater-than sign] L 1. In this situation, it is found that the SHG output in backward transmission process will tend to 2A0/2 as long as the intensity of the input FW is sufficiently high. Thus Eq. (12) can be simplified as:


Ordering (η+)2=2|y2+|2/A02 to describe the SHG conversion efficiency of the forward process, Eq. (18) can be written as P=1η+1+η+. As a result, complete reciprocity and non-reciprocity are corresponding to P = 0 and P = 1, respectively.

Figures for the discussion of the non-reciprocity parameter are all calculated by MATLAB on the basis of the analytical solutions. Figure 3 is obtained with Eq. (15) under the approximate condition that L  L 1 [dbl greater-than sign] L 1 while Fig. 4 is obtained by accurately substituting the Eqs (10) and (11) into Eq. (12). The digital trends for the corresponding points match well between the two figures, which demonstrated that the approximation above is reasonable.


This work is supported by the National Natural Science Foundation of China (Grant Nos. 11574146, 11374150 and 11504166), the Jiangsu Science Foundation (BK20150563), and Priority Academic Program Development of Jiangsu Higher Education Institutions of China (PAPD).

Author Contributions

Author Contributions

C.Z., Y.Q.Q. and Y.Y.Z. conceived the theory. J.C.J. helped to complete the theory and derived the equations. R.E.L. did numerical simulations and wrote the main manuscript text. X.F., R.Z.Z. and X.H.H. helped to improve the numerical calculations and prepare the figures. All authors reviewed the manuscript and gave final approval for publication.


Competing Interests

The authors declare that they have no competing interests.


Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Contributor Information

Chao Zhang, nc.ude.ujn@hchz.

Yiqiang Qin, nc.ude.ujn@niqqy.


1. Kulishov M, Laniel J, Bélanger N. Nonreciprocal waveguide Bragg gratings. Opt. Express. 2005;13:3068–3078. doi: 10.1364/OPEX.13.003068. [PubMed] [Cross Ref]
2. Khanikaev AB, Steel MJ. Low-symmetry magnetic photonic crystals for nonreciprocal and unidirectional devices. Opt. Express. 2009;17:5265–5272. doi: 10.1364/OE.17.005265. [PubMed] [Cross Ref]
3. Peng B, et al. Nonreciprocal light transmission in parity-time-symmetric whispering-gallery microcavities. Nature Phys. 2013;10:394–398. doi: 10.1038/nphys2927. [Cross Ref]
4. Pačebutas V, Krotkus A, Šimkienė I, Viselga R. Electric and photoelectric properties of diode structures in porous silicon. J. Appl. Phys. 1995;77:2501–2507. doi: 10.1063/1.358779. [Cross Ref]
5. Xue J, Zhao W, Su X. Research on laser damage of avalanche photoelectric diode. Appl. Laser. 2000;20:127–128.
6. Liang B, Yuan B, Cheng JC. Acoustic diode: rectification of acoustic energy flux in one-dimensional systems. Phys. Rev. Lett. 2009;103:104301. doi: 10.1103/PhysRevLett.103.104301. [PubMed] [Cross Ref]
7. Liang B, Guo XS, Tu J, Zhang D, Cheng JC. An acoustic rectifier. Nature Mater. 2010;9:989–992. doi: 10.1038/nmat2881. [PubMed] [Cross Ref]
8. Li B, Wang L, Casati G. Thermal diode: rectification of heat flux. Phys. Rev. Lett. 2004;93:10445–10458. [PubMed]
9. Chen Z, et al. A photon thermal diode. Nat. Commun. 2014;5:5446–5446. doi: 10.1038/ncomms6446. [PMC free article] [PubMed] [Cross Ref]
10. Kerr, D. E. Appendix A: Application of the Lorentz reciprocity theorem to scattering (IET Digital Library, 1987).
11. Feng L, et al. Nonreciprocal light propagation in a silicon photonic circuit. Science. 2011;333:729–733. doi: 10.1126/science.1206038. [PubMed] [Cross Ref]
12. Gallo K, Assanto G, Parameswaran KR, Fejer MM. All-optical diode in a periodically poled lithium niobate waveguide. Appl. Phys. Lett. 2001;79:314–316. doi: 10.1063/1.1386407. [Cross Ref]
13. Gallo K, Assanto G. All-optical diode based on second-harmonic generation in an asymmetric waveguide. J. Opt. Soc. Am. B. 1999;16:267–269. doi: 10.1364/JOSAB.16.000267. [Cross Ref]
14. Qian X, et al. Electro-optic tunable optical isolator in periodically poled LiNbO3. J. Appl. Phys. 2011;109:053111. doi: 10.1063/1.3549185. [Cross Ref]
15. Yu ZY, et al. Tunable broadband isolator based on electro-optically induced linear gratings in a nonlinear photonic crystal. Opt. Lett. 2010;35:3327–3329. doi: 10.1364/OL.35.003327. [PubMed] [Cross Ref]
16. Jalas D, et al. What is - and what is not - an optical isolator. Nature Photon. 2013;7:579–582. doi: 10.1038/nphoton.2013.185. [Cross Ref]
17. Yuan L, Shi J, Chen X. All-optical isolator under arbitrary linearly polarized fundamental wave in an optical superlattice. Appl. Opt. 2011;50:6352–6357. doi: 10.1364/AO.50.006352. [PubMed] [Cross Ref]
18. Armstrong JA, Bloembergen N, Ducuing J, Pershan PS. Interactions between light waves in a nonlinear dielectric. Phys. Rev. 1962;127:1918–1939. doi: 10.1103/PhysRev.127.1918. [Cross Ref]
19. Liu X, Zhang H, Zhang M. Exact analytical solutions and their applications for interacting waves in quadratic nonlinear medium. Opt. Express. 2002;10:83–97. doi: 10.1364/OE.10.000083. [PubMed] [Cross Ref]
20. Zhu SN, Zhu YY, Ming NB. Quasi-phase-matched third-harmonic generation in a quasi-periodic optical superlattice. Science. 1997;278:843–846. doi: 10.1126/science.278.5339.843. [Cross Ref]
21. Chen X, et al. Multiple quasi phase-matching in a nonperiodic domain-inverted optical superlattice. Phys. Rev. A. 2004;69:013818. doi: 10.1103/PhysRevA.69.013818. [Cross Ref]
22. Ellenbogen T, Arie A, Saltiel SM. Noncollinear double quasi phase matching in one-dimensional poled crystals. Opt. Lett. 2007;32:262–264. doi: 10.1364/OL.32.000262. [PubMed] [Cross Ref]
23. Rustagi KC, Mehendale SC, Meenakshi S. Optical frequency conversion in quasi-phase-matched stacks of nonlinear crystals. IEEE J. Quantum Electron. 1982;18:1029–1041. doi: 10.1109/JQE.1982.1071650. [Cross Ref]
24. Milton, A., Irene, A. S. & Robert, H. R. Handbook of mathematical functions with formulas, graphs, and mathematical tables (U.S. Govt. Print. Off., 1964).

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group