2.3 Second-order nonlinear optics
L1 Phenomenon
The polarization of LN contains a term quadratic in the optical field (χ(2)), so two photons merge into one at the sum frequency (SFG; SHG when both are equal) or one photon splits into two (DFG; parametric amplification). The largest coefficient, d33, couples fields all polarized along c.
Dispersion makes 775 nm light slower than 1550 nm light in bulk CLN: the second harmonic slips out of phase after the coherence length Lc = π/Δk = 9.495 µm derived and converts back. Reversing the ferroelectric domain every Lc (period Λ = 18.99 µm, PPLN) restarts the growth in every domain, and the harmonic power grows as L2. In a TFLN waveguide etched from a sub-micrometre film the period is near 4 µm, set by waveguide dispersion reported; measured normalized efficiencies reach thousands of percent per watt per square centimetre (L2).
d(z) |±±±±±±±±±±±±±±±±±±±±±±±±±±±±±±±±±±±±±±±±±±±±±±±±±±±±±±±±±±±±|
0 z = L = 10.0 mm | P2ω(L), RK4 | 2.208e−03 W | derived |
| P2ω(L) / Pω(0) | 2.208e+01 % | derived |
| eq. 2.13, undepleted | 2.600e−03 W | derived |
| Pω(L) / Pω(0) | 77.922 % | derived |
| ηnorm (input) | 2600 %/W/cm² | reported |
| ηnorm sin²(πD) L², QPM | 2600.0 %/W | derived |
| Lc = Λ/2 (unpoled) | 9.495 µm | derived |
| |Pω + P2ω − P0| / P0, lossless | 9.5e−15 | derived |
Drawntop: the guide from z = 0 at 2 px per µm (at most the first 450 µm); middle: the whole guide, 775 nm glow scaled to its maximum in each view; dashed: 1550 nm pump, width = remaining power. Poling off: the glow pulses every coherence length. Periodic poling: it grows monotonically; at watts the curve bends from L² to tanh² as the pump thins. Curves: RK4 of the coupled-mode equations (solid), undepleted eq. 2.13 (dashed), P0 tanh² (dotted, QPM). QPM uses the domain-averaged coupling κ² = ηnorm sin²(πD); without poling the full d33 acts with the bulk or TFLN mismatch. Loss and duty cycle are inputs template.
- Second harmonic at L, RK4
- 2.208e−03 W
- Converted P2ω/Pω(0)
- 2.208e+01 %
- Undepleted, eq. 2.13
- 2.600e−03 W
- Pump remaining
- 77.922 %
- Measured ηnorm (input)
- 2600 %/W/cm² reported
- Low-power efficiency, QPM
- 2600.0 %/W
- Coherence length Lc (unpoled)
- 9.495 µm
- Energy check, lossless
- 9.5e−15
L2 Quantitative model
d tensor
Contracted index l: 1 = xx, 2 = yy, 3 = zz, 4 = yz, 5 = xz, 6 = xy. Point group 3m with Kleinman symmetry (d15 = d24 = d31; Kleinman 1962):
Complex-amplitude convention of Boyd 2008; χ(2) = 2d. For two distinct input frequencies (SFG) the prefactor is 4ε0.
Values (pm/V)
No source held tabulates a post-1997 absolute d set: every value below is a legacy-scale compilation, a secondary table without per-cell references, or unverified. No set is used as default.
| Source | Fundamental | d33 | d31 | d22 | Scale | Status |
|---|---|---|---|---|---|---|
| Volk & Wohlecke Table 1.2 (CLN), citing Nikogosyan 2005 | 1.06 µm | 34 | 6 | 3 | legacy | P-comp |
| same, near-SLN | 1.06 µm | 42 | 5 | 2.5 | legacy | P-comp |
| Zhu et al. review Table 1 (no per-cell reference) | 1064 nm | −27.0 | −4.3 | 2.10 | not stated | S |
| Shoji et al. 1997, absolute scale | 1.064 µm | 25.2 claimed | about 4.6 claimed | – | post-1997 | U |
| Kitamoto et al. 1995, parametric fluorescence | 488 / 532 nm pump | – | 4.8 ± 0.5 / 4.3 ± 0.5 claimed | – | post-1997 method | U |
Weis & Gaylord 1985 has no SHG d table. The legacy d33 = 34 pm/V is 35% above the claimed 25.2 pm/V U, which attributes the older values to multiple-reflection errors (claim, abstract only). Miller's rule scales χ(2) with the product of the linear susceptibilities at the three frequencies; Shoji et al. report that it holds poorly for LN (claim) U.
Field mapping (eq. 2.10, no further assumption)
| Cut, mode | Fields | Active coefficient | Phase matching |
|---|---|---|---|
| x-cut, TE, propagation y | E along z at ω and 2ω | d33 | QPM, type 0 (e + e → e) |
| z-cut, TM | E mostly along z | d33 | QPM, type 0 |
| z-cut, TE, in-plane field at α from x | E in xy | TE SH: −2ε0d22E² sin 3α, largest for E along y (propagation x), zero for E along x (propagation y); TM SH via d31 | modal or QPM |
| bulk, type I (o + o → e) | ordinary inputs, no z component | d31, d22 only | birefringent, angle-tuned |
Birefringent phase matching never uses d33: ordinary inputs have no z component, so only the xx, yy, xy columns act.
| P/(2ε0E²) | × d33 | × d31 | × d22 | pm/V |
|---|---|---|---|---|
| x | 0.000 | 0.000 | 0.000 | – |
| y | 0.000 | 0.000 | 0.000 | – |
| z | +1.000 | 0.000 | 0.000 | – |
| onto TE out | +1.000 | 0.000 | 0.000 | – |
| onto TM out | 0.000 | 0.000 | 0.000 | – |
| 1 xx | 2 yy | 3 zz | 4 yz | 5 xz | 6 xy | |
|---|---|---|---|---|---|---|
| x | 0 | 0 | 0 | 0 | d31 | −d22 |
| y | −d22 | d22 | 0 | d31 | 0 | 0 |
| z | d31 | d31 | d33 | 0 | 0 | 0 |
z
|
| o: E along x (TM)
| +: origin
············+----------ky k at 0° from y
|
|
|
x-cut, normal x z-cut TE puts E in the xy plane, the field of an ordinary wave in bulk: rotating it through 180° never drives d33, only d31 and d22, and the TE response repeats every 120°. x-cut TE is driven by d33 alone. Propagation angle from the first in-plane axis (x-cut: y; y-cut and z-cut: x) toward the second; TE: in-plane field perpendicular to k; TM: field along the wafer normal. Numbers: 1.06 µm values of the selected set.
- Propagation k
- ( 0.00, +1.00, 0.00)
- Field E
- ( 0.00, 0.00, +1.00)
- P/(2ε0E²)
- d33 ( 0.00, 0.00, +1.00)
- SH field along TE
- d33
- SH field along TM
- 0
- Coefficients driven
- d33
Phase matching
Wave-vector mismatch with an m-th order grating of period Λ (phase matching at Δk = 0):
Birefringent, type I SHG, bulk CLN, 21 °C. 1550 nm: θpm = 49.55° from c, SH walk-off 2.034° derived; 2000 nm: 45.01° derived. 1064 nm: no angle exists, since ne(532 nm) − no(1064 nm) = +1.46×10−3 > 0 derived; 1064 nm type-I SHG in CLN needs temperature tuning (noncritical phase matching).
Quasi-phase matching. Order m, duty cycle D (Fourier coefficient of the square-wave sign of d; form as in Fejer et al. 1992), so (2/π) d33 at D = 0.5, m = 1:
| Process, bulk CLN, e + e → e, 21 °C | Λ (µm) | Lc = Λ/2 (µm) | Status |
|---|---|---|---|
| SHG 1064 → 532 nm | 6.818 | 3.409 | derived |
| SHG 1550 → 775 nm (also degenerate OPO pumped at 775 nm) | 18.989 | 9.495 | derived |
| SHG 2090 → 1045 nm | 31.07 | 15.54 | derived |
| DFG / OPO: pump 1064, signal 1550, idler 3393.4 nm | 29.91 | 14.96 | derived |
Jundt ne at 24.5 °C gives 18.976 µm for 1550 nm SHG derived. TFLN waveguides shift the period to a few µm because neff disperses differently from bulk ne; x-cut films have been poled at 600 nm reported and 747 nm periods reported. Early thin-film PPLN converters: Chang et al. 2016.
Undepleted and depleted SHG
ηnorm computed from d33 inherits the d spread. Eq. 2.13 at 1550 nm with bulk indices, Aeff = 1 µm² template, deff = (2/π) d33:
| d33 source | d33 (pm/V) | deff (pm/V) | ηnorm (%/W/cm²) | Status |
|---|---|---|---|---|
| Volk & Wohlecke, CLN, legacy scale P-comp | 34.0 | 21.6 | 5828 | derived |
| same, near-SLN P-comp | 42.0 | 26.7 | 8893 | derived |
| Zhu et al. review, |d33| S | 27.0 | 17.2 | 3675 | derived |
| Shoji et al. 1997, claimed U | 25.2 claimed | 16.0 | 3201 claimed | U |
Depleted, at Δk = 0 (Armstrong et al. 1962):
Bulk focusing cannot reach the L2 law: a Gaussian beam diffracts, the optimum confocal parameter is b = L/2.84 reported, and P2ω then grows only linearly in L. Waveguides hold Aeff fixed over centimetres, hence the L2 law and the 1/Aeff gain.
Measured normalized efficiencies
| Device | Value | Normalization | Date | Status |
|---|---|---|---|---|
| LiNbO3 waveguide, field-poled at room temperature, first-order QPM, blue SHG | 600 %/W/cm²; 20.7 mW blue | waveguide | 1993-02-01 | reported |
| TFLN PPLN, 1.5 µm pump, 4 µm period | 2600 %/W/cm² | waveguide | 2018-11-07 | reported |
| TFLN, shallow etch, MgO, 5 mm | 3757 %/W/cm² (939 %/W) | waveguide | 2020-06-19 | reported |
| TFLN, actively monitored poling | 4600 %/W/cm² | waveguide | 2019-08-28 | reported |
| double-layer, polarization-reversed TFLN | 9600 %/W/cm² measured; >10,000 theory; 85% absolute at high power | waveguide | 2024-11-07 | reported |
| TFLN, blue SHG | 33,000 %/W/cm² | waveguide | 2022-05-19 | reported |
| TFLN on sapphire, DFG to 3.66 µm | up to 200 %/W/cm² | waveguide (DFG) | 2021-06-16 | reported |
| TFLN, modal phase matching (no poling), 8 mm | 4.7 %/W | waveguide, %/W | 2018-08-16 | reported |
| z-cut PPLN microring | 250,000 %/W | ring | 2019-11-21 | reported |
| PPLN microring; g/2π about 1.2 MHz | 5,000,000 %/W | ring | 2020-11-23 | reported |
Waveguides: %/W/cm² = P2ω/(Pω²L²); rings: %/W. The two units are not convertible without a model (L3). Dates: issue dates.
Acceptance bandwidths
sinc²(x/2) falls to one half at x = ±2.783, so the full width at half maximum in any parameter q is
| SHG pump, bulk CLN, first-order QPM | FWHM in λω × L, Zelmon 21 °C (nm cm) | FWHM in T × L, Jundt 25 °C (K cm) | dλω/dT, Jundt 25 °C (nm/K) |
|---|---|---|---|
| 1064 nm | 0.205 derived | 3.21 derived | 0.064 derived |
| 1550 nm | 1.176 derived | 10.12 derived | 0.115 derived |
Thermal expansion of the period not included. The temperature acceptance shifts with the reference temperature: 1550 nm SHG 10.14 K cm at 24.5 °C, 10.26 K cm at 21 °C derived.
Group-velocity mismatch, bulk, e-polarized 1550/775 nm: Δng = 0.0905, 302 fs/mm derived. Comparison with a TFLN device: Zhao et al. report a 1 nm linewidth for 5 mm and 0.1 nm/°C tuning reported; bulk CLN scaled to 5 mm gives 2.35 nm derived and 0.115 nm/K derived. Waveguide dispersion lets nanophotonic PPLN cancel group-velocity mismatch and run femtosecond SHG with < 100 fJ reported.
| λSH | 0.775 µm | derived |
| Λ, order m = 1 | 18.989 µm | derived |
| Lc = π/Δk | 9.495 µm | derived |
| deff/d33 = 2 sin(mπD)/(mπ) | 0.6366 | derived |
| θpm type I (o + o → e) | 49.55° | derived |
| walk-off ρ of the e wave | 2.034° | derived |
| FWHM in λω at L (SHG) | 1.176 nm | derived |
| FWHM in T at L (SHG, Jundt) | 9.64 K | derived |
| dλω/dT at fixed Λ (SHG, Jundt) | 0.1205 nm/K | derived |
| GVM 1/vg,2ω − 1/vg,ω (SHG) | 301.9 fs/mm | derived |
| grating Λg = mΛ(40 °C) + ΔΛ (SHG, Jundt) | 18.9264 µm | derived |
| Δk at λω, T, Λg (eq. 2.12) | 0.000 rad/mm | derived |
| ΔkL/2 | 0.000 rad | derived |
| sinc²(ΔkL/2) | 1.0000 | derived |
| SH power / matched 1st order, D = 0.5 | 1.0000 | derived |
| FWHM in Λ at L | 31.73 nm | derived |
bulk CLN, order 1, Λ 18.989 µm, D 0.50 |+++++++++----------+++++++++----------+++++++++---| TFLN, Λ 4 µm reported (Wang et al. 2018, 1.5 µm pump) |++--++--++--++--++--++--++--++--++--++--++--++--++| 0 10 20 30 40 50 µm SH 775 nm |########################################| 1.000
Drawndomains at 2 px per µm: the computed bulk period (1550 nm SHG: 38 px) above the reported TFLN period (8 px). Below, a fixed grating Λg = mΛ1(40 °C, Jundt) + ΔΛ: the generated beam, in the color of its wavelength, is as bright as sinc²(ΔkL/2) allows; heat the crystal or offset the period and it fades, then flickers through the side lobes. The curve shows the same sinc² versus the selected parameter, area inside the FWHM filled, diamond = operating point. First wavelength: fundamental (SHG), first input (SFG) or pump (DFG); second: second input (SFG) or signal (DFG, idler from energy conservation). Curve and beam: Jundt ne, defined from 20 to 250 °C, so a curve ends where T leaves that window; brightness relative to a matched first-order grating with D = 0.5. Period expansion with T not included. Type-I angle: Zelmon 21 °C.
- Output wavelength
- 0.775 µm
- Period Λ, order 1
- 18.989 µm
- Coherence length
- 9.495 µm
- deff / d33
- 0.6366
- Type-I angle from c
- 49.55°
- FWHM in wavelength
- 1.176 nm
- FWHM in temperature
- 9.64 K
- Tuning dλ/dT
- 0.1205 nm/K
- Group-velocity mismatch
- 301.9 fs/mm
- Grating Λg
- 18.9264 µm
- Δk
- 0.000 rad/mm
- ΔkL/2
- 0.000 rad
- sinc²(ΔkL/2)
- 1.0000
- SH power / matched 1st order
- 1.0000
- FWHM in Λ
- 31.73 nm
L3 Normalizations and conventions
Why the efficiency numbers do not compare directly
| Normalization | Definition | Valid when | Breaks when |
|---|---|---|---|
| %/W/cm² (ηnorm) | P2ω/(Pω²L²), low-power slope, phase matched | uniform QPM, no loss, undepleted; on-chip powers | nonuniform film or poling (effective length < L, so ηnorm falls with L); loss; quoted with fiber-coupled powers |
| %/W | P2ω/Pω² for one device | rings and resonators; also whole waveguides | compares devices of different length or Q; a ring value scales with resonant build-up |
| absolute % | P2ω/Pω at a stated pump | high-power, depleted regime | says nothing about low-power efficiency |
| pulse energy (fJ, pJ) | energy at which a stated conversion or gain is reached | femtosecond, dispersion-engineered devices | not convertible to CW numbers without the pulse shape |
Ring build-up from coupled-mode theory, H/ħ = g(a²b† + h.c.), energy decay rates κa (fundamental), κb (SH), external rates κa,ex, κb,ex, zero detuning derived:
At fixed g a ring's %/W therefore scales as Q3. Published definitions of g differ by numerical factors, so g values from different papers are compared only after conversion to one Hamiltonian. ηnorm scales as 1/(λ²Aeff) at fixed deff (eq. 2.13); Aeff itself scales roughly with λ² in a geometrically scaled guide, and deff disperses (Miller). Blue SHG at 33,000 %/W/cm² reported and 1550 nm SHG at 2600 to 9600 %/W/cm² are therefore not ranked by one number.
Overlap
For dominant field components normalized to unit power, the integral in the numerator runs where d is nonzero (LN, with sign flips in poled regions). Some authors fold neff into this definition.
Coupled-mode equations
|A|² = power, fields written as A exp(i(ωt − kz)), the complex conjugate of the convention in Boyd, so phases of A carry the opposite sign (Armstrong et al. 1962). The equations conserve Pω + P2ω (Manley-Rowe) and reproduce eq. 2.13 and the tanh² law; the SHG figure in L1 integrates them by RK4.
Error sources in measured ηnorm
Facet-loss calibration at two wavelengths (the SH often exits multimode); poling duty-cycle errors and missing domains (deff falls as sin πD; Fejer et al. 1992); film-thickness variation that moves the phase-matching wavelength along the device, which adapted poling compensates (Chen et al. 2024); photorefraction and green-induced infrared absorption at high SH power (chapter 04). The d scale carries the largest systematic uncertainty: 25 vs 34 pm/V changes ηnorm by 1.85× derived.
Reviews: TFLN platform (Zhu et al. 2021); LN photonics from UV to microwave (Boes et al. 2023); LNOI phase-matching schemes (BPM, modal, QPM) and tuning effects (You et al. 2024).
| ηnorm L² at the set L | 650.0 %/W | derived |
| ring, critical coupling, Qa = Qb = Q | 1.373e+05 %/W | derived |
| κa/2π = νa/Q (1550 nm) | 1934.1 MHz | derived |
TFLN PPLN, 1.5 µm pump, 4 µm period reported
Drawnupper panel per watt and square centimetre, lower panel per watt for one device; a point never appears on a panel of another unit. Drag Q: the ring line climbs as Q³ past the ring records while the upper panel stays still. Marker fill: color of the generated light where the source gives the pump wavelength; IR output dashed. Values: the measured-efficiency table in L2, each as reported. Ring model: 1550 nm fundamental, critical coupling, Qa = Qb = Q; g/2π is a display parameter template.
- Record
- Wang et al. 2018
- Issue date
- 2018-11-07
- Value
- 2600 %/W/cm²
- Normalization
- waveguide, normalized
- ηnorm L² at the set L
- 650.0 %/W
- Ring model %/W (critical, Qa = Qb)
- 1.373e+05
- κa/2π = νa/Q (1550 nm)
- 1934.1 MHz
References
| Armstrong et al. 1962 | Interactions between Light Waves in a Nonlinear Dielectric. Physical Review 127 1918-1939 (1962). | doi:10.1103/PhysRev.127.1918 |
| Boes et al. 2023 | Lithium niobate photonics: Unlocking the electromagnetic spectrum. Science 379 eabj4396 (2023). | doi:10.1126/science.abj4396 |
| Boyd & Kleinman 1968 | Parametric Interaction of Focused Gaussian Light Beams. Journal of Applied Physics 39 3597-3639 (1968). | doi:10.1063/1.1656831 |
| Boyd 2008 | Nonlinear Optics, 3rd ed.. Academic Press (2008). | doi:10.1016/b978-0-12-369470-6.00014-9 |
| Chang et al. 2016 | Thin film wavelength converters for photonic integrated circuits. Optica 3 531 (2016). | doi:10.1364/OPTICA.3.000531 |
| Chen et al. 2024 | Adapted poling to break the nonlinear efficiency limit in nanophotonic lithium niobate waveguides. Nature Nanotechnology 19 44-50 (2024). | doi:10.1038/s41565-023-01525-w |
| Fejer et al. 1992 | Quasi-phase-matched second harmonic generation: tuning and tolerances. IEEE Journal of Quantum Electronics 28 2631-2654 (1992). | doi:10.1109/3.161322 |
| Jankowski et al. 2020 | Ultrabroadband nonlinear optics in nanophotonic periodically poled lithium niobate waveguides. Optica 7 40 (2020). | doi:10.1364/OPTICA.7.000040 |
| Kitamoto et al. 1995 | Optical Review 2 280 (1995). | no DOI |
| Kleinman 1962 | Nonlinear Dielectric Polarization in Optical Media. Physical Review 126 1977-1979 (1962). | doi:10.1103/PhysRev.126.1977 |
| Li et al. 2024 | Highly efficient second-harmonic generation in a double-layer thin-film lithium niobate waveguide. Advanced Photonics Nexus 3 1-7 (2024). | doi:10.1117/1.APN.3.6.066009 |
| Lu et al. 2019 | Periodically poled thin-film lithium niobate microring resonators with a second-harmonic generation efficiency of 250,000%/W. Optica 6 1455 (2019). | doi:10.1364/OPTICA.6.001455 |
| Lu et al. 2020 | Toward 1% single-photon anharmonicity with periodically poled lithium niobate microring resonators. Optica 7 1654 (2020). | doi:10.1364/OPTICA.403931 |
| Luo et al. 2018 | Highly tunable efficient second-harmonic generation in a lithium niobate nanophotonic waveguide. Optica 5 1006 (2018). | doi:10.1364/OPTICA.5.001006 |
| Mishra et al. 2021 | Mid-infrared nonlinear optics in thin-film lithium niobate on sapphire. Optica 8 921 (2021). | doi:10.1364/OPTICA.427428 |
| Nagy & Reano 2020 | Submicrometer periodic poling of lithium niobate thin films with bipolar preconditioning pulses. Optical Materials Express 10 1911 (2020). | doi:10.1364/OME.394724 |
| Park et al. 2022 | High-efficiency second harmonic generation of blue light on thin-film lithium niobate. Optics Letters 47 2706 (2022). | doi:10.1364/OL.455046 |
| Rao et al. 2019 | Actively-monitored periodic-poling in thin-film lithium niobate photonic waveguides with ultrahigh nonlinear conversion efficiency of 4600 %W−1cm−2. Optics Express 27 25920 (2019). | doi:10.1364/OE.27.025920 |
| Shoji et al. 1997 | Absolute scale of second-order nonlinear-optical coefficients. Journal of the Optical Society of America B 14 2268 (1997). | doi:10.1364/JOSAB.14.002268 |
| Volk & Wöhlecke 2009 | Lithium Niobate. Springer Series in Materials Science (2009). | doi:10.1007/978-3-540-70766-0 |
| Wang et al. 2018 | Ultrahigh-efficiency wavelength conversion in nanophotonic periodically poled lithium niobate waveguides. Optica 5 1438 (2018). | doi:10.1364/OPTICA.5.001438 |
| Yamada et al. 1993 | First-order quasi-phase matched LiNbO3 waveguide periodically poled by applying an external field for efficient blue second-harmonic generation. Applied Physics Letters 62 435-436 (1993). | doi:10.1063/1.108925 |
| You et al. 2024 | Lithium niobate on insulator – fundamental opto‐electronic properties and photonic device prospects. Nanophotonics 13 3037-3057 (2024). | doi:10.1515/nanoph-2024-0132 |
| Zhao et al. 2020 | Shallow-etched thin-film lithium niobate waveguides for highly-efficient second-harmonic generation. Optics Express 28 19669 (2020). | doi:10.1364/OE.395545 |
| Zhao et al. 2020 | Poling thin-film x-cut lithium niobate for quasi-phase matching with sub-micrometer periodicity. Journal of Applied Physics 127 193104 (2020). | doi:10.1063/1.5143266 |
| Zhu et al. 2021 | Integrated photonics on thin-film lithium niobate. Advances in Optics and Photonics 13 242 (2021). | doi:10.1364/AOP.411024 |