Wave Equation for Non-Linear Holographic Materials


Posted

A few years ago I needed to create a tool to compute the incoming angles of the laser beams required to create a specified hologram. The tool was a 3-dimensional generalization of the equations in Kogelnik [1]. At the time I wondered about the effect that material non-linearity might have on the results. Kogelnik assumes that the material has a single sine-wave perturbation of the index of refraction in the material. The higher spatial harmonics will change the relative power of the outgoing diffraction orders, much like the blaze of a grating. Here I work out a generalization of the wave equations, insofar as to include the harmonics.

Light interaction with a hologram is described by the scalar wave equation,

∇ 2 E + k 2 E = 0 , nabla^2 E + k^2 E = 0 ","

(1)

where E (the electric field) and k (related to the dielectric constant and conductivity) are complex functions over space. Since the material making up a hologram is periodic, k can be expanded as a harmonic series

k 2 = ∑ m B m ( r ) e i v m ⋅ r . k^2 = sum from m B sub m ( vec r ) e^ { i vec v sub m cdot vec r } "."

(2)

Bm is a slowly changing function of space. Given that the material is periodic, it is somewhat intuitive that the electric field will be periodic as well, so E can also be expanded as the harmonic series. At this point we shall also assume that the hologram is thick with respect to the wavelength, so that the electric field intensity is changes slowly over space. This gives

E = ∑ j A j ( r ) e i u j ⋅ r . E = sum from j A sub j ( vec r ) e^{ i vec u sub j cdot vec r } "."

(3)

Inserting both harmonic series into the wave equation gives

∑ j [ − u j 2 A j + 2 i u j ⋅ ∇ A j + ∇ 2 A j + A j ∑ m B m e i v m ⋅ r ] e i u j ⋅ r = 0 . sum from j left [ - u sub j sup 2 A sub j + 2 i vec u sub j cdot vec nabla A sub j + nabla^2 A sub j + A sub j sum from m B sub m e^ { i vec v sub m cdot vec r } right ] e^{ i vec u sub j cdot vec r } = 0 "."

(4)

Next generalize (Aj and uj) to (Ajm and ujm) where m corresponds to the matching Bm.

∑ j , m [ − u jm 2 A jm + 2 i u jm ⋅ ∇ A jm + ∇ 2 A jm ] e i u jm ⋅ r + A jm B m e i ( u jm + v m ) ⋅ r = 0 . sum from {j,m} left [ - u sub jm sup 2 A sub jm + 2 i { vec u } sub jm cdot vec nabla A sub jm + nabla^2 A sub jm right ] e^{ i { vec u } sub jm cdot vec r } + A sub jm B sub m e^ { i ( { vec u } sub jm + { vec v } sub m ) cdot vec r } = 0 "."

(5)

Now specify j as corresponding to the order of the overtone of v,

u jm = u o + j v m . { vec u } sub jm = { vec u } sub o + j { vec v } sub m "."

(6)

where uo is the incoming E-wave direction. This gives

∑ j , m [ − ( u o 2 + 2 j u o ⋅ v m + j 2 v m 2 ) A jm + 2 i ( u o + j v m ) ⋅ ∇ A jm + ∇ 2 A jm ] e i ( u o + j v m ) ⋅ r + A jm B m e i ( u o + ( j + 1 ) v m ) ⋅ r = 0 . sum from {j,m} left [ -( u sub o sup 2 + 2 j {vec u} sub o cdot {vec v} sub m + j^2 v sup 2 sub m ) A sub jm + 2 i ( {vec u} sub o + j {vec v} sub m ) cdot vec nabla A sub jm + nabla^2 A sub jm right ] e^{ i ( {vec u} sub o + j {vec v} sub m ) cdot vec r } newline + A sub jm B sub m e^ { i ( {vec u} sub o + (j+1) {vec v} sub m ) cdot vec r } = 0 "."

(7)

Now renumbering the j-1 term results in

∑ j , m [ − ( u o 2 + 2 j u o ⋅ v m + j 2 v m 2 ) A jm + 2 i ( u o + j v m ) ⋅ ∇ A jm + ∇ 2 A jm + B m A j − 1 , m ] e i ( u o + j v m ) ⋅ r = 0 . sum from {j,m} left [ -( u sub o sup 2 + 2 j {vec u} sub o cdot {vec v} sub m + j^2 v sup 2 sub m ) A sub jm + 2 i ( {vec u} sub o + j {vec v} sub m ) cdot vec nabla A sub jm + nabla^2 A sub jm + B sub m A sub {j-1,m} right ] e^{ i ( {vec u} sub o + j {vec v} sub m ) cdot vec r } = 0 "."

(8)

Since the various vm terms are most likely harmonically related, sums across equal jvm's are zero. Therefore the expression is zero everywhere if and only if, for every unique C,

∑ j v m = C − ( u o 2 + 2 j u o ⋅ v m + j 2 v m 2 ) A jm + 2 i ( u o + j v m ) ⋅ ∇ A jm + ∇ 2 A jm + B m A j − 1 , m = 0 . sum from { j vec v sub m = vec C } { -( u sub o sup 2 + 2 j {vec u} sub o cdot {vec v} sub m + j^2 v sup 2 sub m ) A sub jm + 2 i ( {vec u} sub o + j {vec v} sub m ) cdot vec nabla A sub jm + nabla^2 A sub jm + B sub m A sub {j-1,m} } = 0 "."

(9)

As an example, the j=0 term is

− u o 2 A o + 2 i u o ⋅ ∇ A o + ∇ 2 A o + ∑ m B m A − 1 , m = 0 . - u sub o sup 2 A sub o + 2 i {vec u} sub o cdot vec nabla A sub o + nabla^2 A sub o + sum from m B sub m A sub {-1,m} = 0 "."

(10)

The A's are solved for a given set of B's by an eigenvalue problem represented by a matrix which is mostly zeros. Only the diagonal and one off-diagonal are populated.

For a thick hologram, the A's are constant across the material and only change from one side to the other. This results in a second-order differential equation in the material depth. The lowest B-term (after the constant) is the wave programmed in by the original interference setup, while higher terms are due to non-linearity in the material.

Reference

[1] Kogelnik, Herwig. Coupled Wave Theory for Thick Hologram Gratings. The Bell System Technical Journal. Volume 48, number 9, page 2909-2947. November 1969.


Categories Holographics, Optics

← Older Newer →