PPLN Λ = 6.82 µm · 1064 → 532 nm · 2 px/µm

3.6 Acousto-optics

L1 Phenomenon

Bragg diffraction by sound

A sound wave in a crystal is a moving grating of compressed and stretched regions, and the refractive index follows the strain (photoelastic effect). Light crossing it is diffracted and shifted in frequency by exactly the sound frequency. In LN the same crystal can launch the sound piezoelectrically, so a modulator, frequency shifter or deflector needs no glued-on transducer (acousto-optic effect).

Sweep the RF frequency: the diffracted beam swings like a scanner, the basis of an acousto-optic deflector, and the light frequency of the +1 order rises by exactly the RF frequency.

Bragg diffraction by a sound wave in LN
Sound

wavevectors to scale: kin + K = kout, |k| = 2πn/λ₀, K = 2π/Λ

k_ink_outK

K/k = 1.11e−1

n (e), Zelmon 1997
2.1376
v sound
6544.6 m/s
Λ = v/f
  6.54 µm
θB inside
 3.176°
θB in air
 6.801°
deflection 2θB, air
13.602°
shift of the +1 order
+1.000 GHz
1550 nm
1.000 GHz
Diffraction efficiency: not computed.
geometry derived velocities: Kushibiki et al. 1999 constants derived 128° Y-X SAW derived

Crystal block, 200 µm wide

soundtransducer+1 order, optical frequency + f0 order
light: spectral color below 760 nm, false color and dashed abovecompressed bands, spacing Λ to scalebeam angles: air values, not exaggerated

Bragg angle in air

 6.801°
acoustic wavelength Λ
  6.54 µm
light frequency shift
+1.000 GHz
inside the crystal
 3.176°
1.000 GHz
1550 nm
Sound
Polarization
geometry derivedn = 2.1376 (e), Zelmon et al. 1997 P-dbvelocities: Kushibiki et al. 1999 constants derived deriveddiffraction efficiency not computed

L2 Photoelastic tensor, figures of merit, devices

Weis & Gaylord 1985 Eqs. 18-20; 633 nm sets
Δ(1n2)ij=pijklSkl(1)

Photoelastic (strain-optic) relation (Weis & Gaylord 1985 Eq. 18 P). Nelson and Lax showed the independent variable is the displacement gradient, not only strain: rotations of volume elements contribute in birefringent crystals, so pijkl is not symmetric in kl for shear waves; in LN p44 = p55 has symmetric and antisymmetric parts, not resolved by the available data (Weis & Gaylord 1985 p. 200 P-comp).

3m form, eight independent coefficients p11, p12, p13, p14, p31, p33, p41, p44, with p22 = p11, p23 = p13, p24 = −p14, p32 = p31, p42 = −p41, p55 = p44, p56 = p41, p65 = p14, p66 = (p11 − p12)/2 (Weis & Gaylord 1985 Eq. 19 P). pIJ is not a symmetric matrix (p14 ≠ p41).

Setp11p12p13p14p31p33p41p44Status
Andrushchak 2009−0.021+0.060+0.172−0.052+0.141+0.118−0.109+0.121S
Avakyants 1976 (Set B)−0.026+0.090+0.133−0.075+0.179+0.071−0.151+0.146P-comp
Krupych 2014−0.023+0.076+0.147+0.057+0.157+0.141+0.051+0.126S
Coquin (Pinnow 1971)−0.020+0.080+0.130−0.080+0.170+0.070−0.150+0.120P-comp
TFLN fit, LN on sapphire (three components)–––––+0.038−0.038+0.085S

Avakyants and Andrushchak at 633 nm; the wavelengths of the Coquin and Krupych sets are not recorded in the sources held. No set at 1550 nm is known; TFLN devices at 1550 nm reuse 633 nm values. Uncertainties of Andrushchak 2009: ±0.007 to ±0.029 (via Sarabalis et al. 2020 Table I).

Indirect (piezo-electro-optic) contribution. Strain makes a field through e, the field changes the index through r. It depends on the acoustic propagation direction and matters for p13 and p33: p13(primary) = p13(measured) − 0.043, p33(primary) = p33(measured) − 0.154 (Weis & Gaylord 1985 p. 201 P-comp; propagation direction not stated). Applied to Set B this gives p33(primary) of about −0.08 derived, so which p33 a model uses matters. The offsets are shown, not applied.

Figures of merit

M1=n7p2ρv,M2=n6p2ρv3,M3=n7p2ρv2,M4=n8p2vρ(2)

LN, longitudinal wave along X with extraordinary light at 35° to y: n = 2.20, ρ = 4.7×10³ kg/m³, v = 6.57×10³ m/s (Dixon 1967 and Chang 1976 via Weis & Gaylord 1985 p. 201 P-comp).

QuantityPrintedFrom Eq. 2 with p implied by M2Status
M1 (m² s/kg)6.65×10−7 (Dixon 1967)6.64×10⁻⁷P-comp
M2 (s³/kg)6.99×10−15 (Dixon 1967)input; p = 0.287P-comp
M4 (m4/(s kg))62.9 (Chang 1976)63.0P-comp
M3printed “1.01×10−16 m s² kg”1.01×10⁻¹⁰ m s²/kgU

The four definitions with the stated n, ρ, v reproduce M1 and M4; the printed M3 exponent and unit are inconsistent with them, so the printed M3 is not used derived.

Bragg example. λ0 = 1550 nm, f = 1 GHz, X-longitudinal v = 6544.6 m/s: Λ = 6.54 µm and θB = 6.8° in air derived.

Integrated acousto-optics on LN. Collinear TE0-TE1 conversion by an SH0 wave at 440 MHz in a 1.25 µm wide suspended X-cut waveguide, up to 18% conversion (Sarabalis et al. 2021 reported); acousto-optic modulators in X-cut LN on sapphire (Sarabalis et al. 2020 reported); AO modulator VπL of 0.03 V cm in TFLN with chalcogenide, unsuspended (Wan et al. 2022 reported); a bulk-wafer longitudinal piezoelectric resonant photoelastic modulator (Atalar et al. 2022 S).

L3 Sign and frame pitfalls

p14, p41 vs c14

1. Pair signs only within one frame. The p14 and p41 signs flip with c14 (and e22, d22) under a 180° rotation about Z. Never combine p14, p41 from one paper with c14, e22 from another without checking the frame.

2. What is known. Weis & Gaylord 1985 list Avakyants' p14, p41 < 0 next to Smith and Welsh's c14 > 0 in a paper that states the IRE frame, without saying it checked Avakyants' own frame; Krupych 2014 pairs c14 > 0 with p14, p41 > 0 (via Mys et al. 2015 S); Mytsyk et al. 2021 use S14 > 0 (c14 < 0) S, a frame that is non-standard for LN if d22 > 0 is imposed (inference). The common TFLN pairing (Kushibiki et al. c14 > 0 with Set B) matches the Weis & Gaylord compilation but is not confirmed at the primary level.

3. What a mismatch does. For a film with in-plane unit vector b (light polarization and dominant strain both along b), the longitudinal coefficient is p′bbbb = pijklbibjbkbl. On an X-cut film with b = (0, cos φ, sin φ), φ from +Y toward +Z derived:

Set, X-cutp′ at φ = 0, 45, 90, 135°Maximum
Andrushchak 2009, A: as listed−0.021, +0.304, +0.118, +0.143+0.304 at 44.3°
Andrushchak 2009, B: p14, p41 negated−0.021, +0.143, +0.118, +0.304+0.304 at 135.7°
Avakyants 1976 (Set B), A: as listed−0.026, +0.348, +0.071, +0.122+0.352 at 41.7°
Avakyants 1976 (Set B), B: p14, p41 negated−0.026, +0.122, +0.071, +0.348+0.352 at 138.3°

On X-cut the flip mirrors the optimum orientation from about 45° to about 135°; the magnitude is unchanged. On Y-cut (b = (cos φ, 0, sin φ)) the flip changes nothing (Andrushchak: maximum +0.230 at 50.7 and 129.3°; Set B maximum +0.238 at 48.3 and 131.7°). On Z-cut p′bbbb = p11 for every in-plane direction. Which orientation is right on X-cut depends on the unresolved pairing in item 2.

The Andrushchak values at 45 and 135°, 0.304 (X-cut) and 0.223 (Y-cut), equal the “p22′” values in the author's public 2018 note on photoelastic tensor rotation (outside5sigma.com, 2018 S), where the X-cut angle was relabelled from 45 to 135°. A p14/p41 sign or frame inconsistency produces such a relabelling; a change of the θ = 0 reference, also recorded in that note, is an alternative explanation (inference).

Toggle the film cut: on X-cut the two wafers disagree about the bright direction; on Y-cut they agree.

Longitudinal photoelastic coefficient p′bbbb(φ) in the film plane, 633 nm

X-cut, φ from +Y toward +Z; radius |p′|

0.10.20.30.4+Y+Z

Y-cut, φ from +X toward +Z; radius |p′|

0.10.20.30.4+X+Z

solid A: p₁₄, p₄₁ as listed (< 0); dashed B: negated

X-cut φAB
0°−0.021−0.021
45°+0.304+0.143
90°+0.118+0.118
135°+0.143+0.304
Film
045.0°
p′ at φ, A / B
+0.3040 / +0.1430
maximum, A
+0.3041 at φ = 044.3°
maximum, B
+0.3041 at φ = 135.7°
Sign pairing not resolved. Whether p₁₄, p₄₁ as listed belong to the same frame as c₁₄ > 0 (IEEE 176-1987) is not confirmed at the primary level for any set; A and B are equally possible. On X-cut the choice moves the optimum between about 45° and about 135°; on Y-cut and Z-cut it changes nothing. Frame of this set: not stated in the sources (via Sarabalis et al. 2020 Table I and Mytsyk et al. 2021 Table 6).compiled in Weis & Gaylord 1985 (IRE 1949 axes) next to c14 > 0; the frame of the 1976 measurement is not stated.not stated; quoted with c14 > 0 (via Mys et al. 2015), with p14, p41 opposite in sign to the Weis & Gaylord 1985 compilation.IRE 1949 axes as compiled in Weis & Gaylord 1985.
Andrushchak 2009 SAvakyants 1976 (Set B) P-compKrupych 2014 SCoquin (Pinnow 1971) P-comp rotated coefficient derived frame IEEE-176-1987

Pairing A: p₁₄, p₄₁ as listed

+Y+Z
+0.304p′ along the beam

Pairing B: p₁₄, p₄₁ negated

+Y+Z
+0.143p′ along the beam
The sign pairing is open. p₁₄ and p₄₁ flip with c₁₄, e₂₂ and d₂₂ under a 180° rotation about Z. Which of A and B matches c₁₄ > 0 is not confirmed for any set; the size of the effect is the same, only the best direction on an X-cut chip moves, from about 45° to about 135°. Switch to Y-cut: both wafers agree. Frame of this set: not stated in the sources (via Sarabalis et al. 2020 Table I and Mytsyk et al. 2021 Table 6).compiled in Weis & Gaylord 1985 (IRE 1949 axes) next to c14 > 0; the frame of the 1976 measurement is not stated.not stated; quoted with c14 > 0 (via Mys et al. 2015), with p14, p41 opposite in sign to the Weis & Gaylord 1985 compilation.IRE 1949 axes as compiled in Weis & Gaylord 1985.
045.0°
Film
Coefficient set
Andrushchak 2009 SAvakyants 1976 (Set B) P-compKrupych 2014 SCoquin (Pinnow 1971) P-comprose: p′bbbb(φ) > 0, radius 0.1 per 30 px derivedwafer flat drawn for orientation only

4. Optical wavelength. The Avakyants and Andrushchak sets are at 633 nm (Weis & Gaylord 1985 P-comp; Sarabalis et al. 2020 Table I S); no 1550 nm set is known.

5. Rotation (Nelson-Lax) term. For shear waves in a birefringent crystal it adds Δε terms proportional to (εii − εjj) times the rotation; an estimate for one LN nanobeam geometry put it below 4% of the photoelastic coupling (author's public 2018 note, outside5sigma.com S; unreviewed). Optomechanical coupling calculations on TFLN include moving-boundary and photoelastic terms; the rotation term is usually omitted.

References

3.6
SourceReferenceLink
Weis & Gaylord 1985Weis & Gaylord, Applied Physics A Solids and Surfaces 37 191-203 (1985). Lithium niobate: Summary of physical properties and crystal structuredoi:10.1007/BF00614817
Avakyants et al. 1976Avakyants, Kiselev & Shchitkov, Soviet Physics - Solid State 18 899 (1976)none
Coquin 1971Coquin, cited in D. A. Pinnow, Handbook of Lasers (CRC) (1971)none
Andrushchak et al. 2009Andrushchak et al., Journal of Applied Physics 106 073510 (2009). Complete sets of elastic constants and photoelastic coefficients of pure and MgO-doped lithium niobate crystals at room temperaturedoi:10.1063/1.3238507
Krupych et al. 2014Krupych, Savaryn & Vlokh, Applied Optics 53 B1 (2014). Precise determination of full matrix of piezo-optic coefficients with a four-point bending technique: the example of lithium niobate crystalsdoi:10.1364/AO.53.0000B1
Mys et al. 2015Mys, Applied Optics 54, 8176 (2015); arXiv:1508.01075 (2015)none
Mytsyk et al. 2021Mytsyk et al., Crystals 11 1095 (2021). Photoelastic Properties of Trigonal Crystalsdoi:10.3390/cryst11091095
Sarabalis et al. 2020Sarabalis et al., APL Photonics 5 086104 (2020). Acousto-optic modulation in lithium niobate on sapphiredoi:10.1063/5.0012288
Nelson & Lax 1970Nelson & Lax, Physical Review Letters 24 379-380 (1970). New Symmetry for Acousto-Optic Scatteringdoi:10.1103/PhysRevLett.24.379
Nelson & Lax 1971Nelson & Lax, Physical Review B 3 2778-2794 (1971). Theory of the Photoelastic Interactiondoi:10.1103/PhysRevB.3.2778
Dixon 1967Dixon, Journal of Applied Physics 38 5149-5153 (1967). Photoelastic Properties of Selected Materials and Their Relevance for Applications to Acoustic Light Modulators and Scannersdoi:10.1063/1.1709293
Chang 1976Chang, IEEE Transactions on Sonics and Ultrasonics 23 2-21 (1976). I. Acoustooptic Devices and Applicationsdoi:10.1109/T-SU.1976.30835
Zelmon et al. 1997Zelmon, Small & Jundt, JOSA B 14 3319 (1997)doi:10.1364/JOSAB.14.003319
Sarabalis et al. 2021Sarabalis et al., Optica 8 477 (2021). Acousto-optic modulation of a wavelength-scale waveguidedoi:10.1364/OPTICA.413401
Wan et al. 2022Wan et al., Light: Science & Applications 11 145 (2022). Highly efficient acousto-optic modulation using nonsuspended thin-film lithium niobate-chalcogenide hybrid waveguidesdoi:10.1038/s41377-022-00840-6
Atalar et al. 2022Atalar et al., Nature Communications 13 1526 (2022). Longitudinal piezoelectric resonant photoelastic modulator for efficient intensity modulation at megahertz frequenciesdoi:10.1038/s41467-022-29204-9
Jiang et al. 2019Jiang et al., Optica 6 845 (2019). Lithium niobate piezo-optomechanical crystalsdoi:10.1364/OPTICA.6.000845
Jiang 2018Jiang (2018). Rotation of photoelastic tensor and piezoelectric tensor (outside5sigma.com)outside5sigma.com/tutorial/rotation-of-photoelastic-tensor-and-piezoelectric/