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 . 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,
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
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
Inserting both harmonic series into the wave equation gives
Next generalize (Aj and uj) to (Ajm and ujm) where m corresponds to the matching Bm.
Now specify j as corresponding to the order of the overtone of v,
where uo is the incoming E-wave direction. This gives
Now renumbering the j-1 term results in
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,
As an example, the j=0 term is
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.
 Kogelnik, Herwig. Coupled Wave Theory for Thick Hologram Gratings. The Bell System Technical Journal. Volume 48, number 9, page 2909-2947. November 1969.