# Harmonic Fields

This continues from the optical fields and Maxwell's equations tutorial.

Optical fields are harmonic fields that vary sinusoidally with time. The field vectors defined in the preceding tutorials are all real quantities. For harmonic fields, it is always convenient to use *complex fields*. We define the space- and time-dependent complex electric field, **E**(**r**, *t*), through its relation to the real electric field, * E*(

**r**,

*t*):

\[\tag{39}\pmb{E}(\mathbf{r},t)=\mathbf{E}(\mathbf{r},t)+\mathbf{E}^*(\mathbf{r},t)=\mathbf{E}(\mathbf{r},t)+\text{c.c.}\]

where c.c. means the complex conjugate. In our convention, **E**(**r**, *t*) contains the complex field components that vary with time as exp(-i*ωt*) with positive values of *ω*, while **E***(**r**, *t*) contains those varying with time as exp(i*ωt*) with positive *ω*, or exp(-i*ωt*) with negative *ω*. The complex fields of other field quantities are similarly defined.

In some literature, the complex field is defined through a relation with the real field as \(\pmb{E}(\mathbf{r},t)=\frac{1}{2}[\mathbf{E}(\mathbf{r},t)+\mathbf{E}^*(\mathbf{r},t)]\), which differs from the relation in (39) by a factor of 1/2. The magnitude of the complex field defined through this alternative relation is twice that of the complex field defined through (39). As a result, expressions for many quantities may be different under the two different definitions. An example is that of the optical intensity. We have chosen to define the complex field through the relation in (39) without the factor of 1/2 primarily because this definition is more convenient and less confusing in expressing the nonlinear polarizations discussed in later tutorials.

With this definition for the complex fields, all of the linear field equations retain their forms. In particular, Maxwell's equations for the complex optical fields are

\[\tag{40}\boldsymbol{\nabla}\times\mathbf{E}=-\frac{\partial\mathbf{B}}{\partial t}\]\[\tag{41}\boldsymbol{\nabla}\times\mathbf{H}=\frac{\partial\mathbf{D}}{\partial t}\]\[\tag{42}\boldsymbol{\nabla}\cdot\mathbf{D}=0\]\[\tag{43}\boldsymbol{\nabla}\cdot\mathbf{B}=0\]

The wave equation in terms of the complex electric field is

\[\tag{44}\boldsymbol{\nabla}\times\boldsymbol{\nabla}\times\mathbf{E}+\frac{1}{c^2}\frac{\partial^2\mathbf{E}}{\partial t^2}=-\mu_0\frac{\partial^2\mathbf{P}}{\partial t^2}\]

while

\[\tag{45}\mathbf{P}(\mathbf{r},t)=\epsilon_0\displaystyle\int\limits_{-\infty}^{\infty}\text{d}\mathbf{r}'\int\limits_{-\infty}^t\text{d}t'\pmb{\chi}(\mathbf{r}-\mathbf{r}',t-t')\cdot\mathbf{E}(\mathbf{r},t')\]

and

\[\tag{46}\begin{align}\mathbf{D}(\mathbf{r},t)&=\epsilon_0\mathbf{E}(\mathbf{r},t)+\epsilon_0\displaystyle\int\limits_{-\infty}^{\infty}\text{d}\mathbf{r}'\int\limits_{-\infty}^t\text{d}t'\pmb{\chi}(\mathbf{r}-\mathbf{r}',t-t')\cdot\mathbf{E}(\mathbf{r}',t')\\&=\int\limits_{-\infty}^{\infty}\text{d}\mathbf{r}'\int\limits_{-\infty}^t\text{d}t'\pmb{\epsilon}(\mathbf{r}-\mathbf{r}',t-t')\cdot\mathbf{E}(\mathbf{r}',t')\end{align}\]

It is important to note that while **P**, **D**, and **E** are complex, **χ**(**r** - **r**', *t* - *t*') and **ε**(**r** - **r**', *t* - *t*') in (45) and (46) are always real.

For a harmonic optical field of wavevector **k** and angular frequency *ω*, its complex electric field can be further written as

\[\tag{47}\mathbf{E}(\mathbf{r},t)=\pmb{\mathcal{E}}(\mathbf{r},t)\exp(\text{i}\mathbf{k}\cdot\mathbf{r}-\text{i}\omega t)\]

where \(\pmb{\mathcal{E}}(\mathbf{r},t)\) is the space- and time-varying field envelope, such as that for a modulated field, a guided field, or an optical pulse. Other complex quantities, such as \(\mathbf{H}(\mathbf{r},t)\), can be similarly expressed. The phase factor in (47) indicates the direction of wave propagation:

\(\text{i}\mathbf{k}\cdot\mathbf{r}-\text{i}\omega t\), forward propagating in **k** direction;

\(-\text{i}\mathbf{k}\cdot\mathbf{r}-\text{i}\omega t\), backward propagating in -**k** direction.

The light intensity, or irradiance, is the power density of the harmonic optical field. It can be calculated by time averaging of the Poynting vector over one wave cycle:

\[\tag{48}\overline{\pmb{S}}=\frac{1}{T}\displaystyle\int\limits_0^T\pmb{E}\times\pmb{H}\text{d}t=2\text{ Re}(\mathbf{E}\times\mathbf{H}^*)\]

where Re() means taking the real part. We can define a *complex Poynting vector*:

\[\tag{49}\mathbf{S}=\mathbf{E}\times\mathbf{H}^*\]

so that

\[\tag{50}\overline{\pmb{S}}=\mathbf{S}+\mathbf{S}^*\]

which has the same form as the relation between the real and complex fields defined in (39) except that the real Poynting vector in this relation is time averaged. The light intensity, *I*, is simply the magnitude of the *real* time-averaged Poynting vector:

\[\tag{51}I = |\overline{\pmb{S}}|=|\mathbf{S}+\mathbf{S}^*|\]

where *I* is in watts per square meter.

For harmonic optical fields, it is often useful to consider the complex fields in the momentum space and frequency domain defined by the following Fourier-transform relations:

\[\tag{52}\mathbf{E}(\mathbf{k},\omega)=\displaystyle\int\limits_{-\infty}^{\infty}\text{d}\mathbf{r}\int\limits_{-\infty}^{\infty}\text{d}t\,\mathbf{E}(\mathbf{r},t)\exp(-\text{i}\mathbf{k}\cdot\mathbf{r}+\text{i}\omega t),\qquad\text{for }\omega\gt0\]\[\tag{53}\mathbf{E}(\mathbf{r},t)=\frac{1}{(2\pi)^4}\displaystyle\int\limits_{-\infty}^{\infty}\text{d}\mathbf{k}\int\limits_0^{\infty}\text{d}\omega\,\mathbf{E}(\mathbf{k},\omega)\exp(\text{i}\mathbf{k}\cdot\mathbf{r}-\text{i}\omega t)\]

Note that **E**(**k**, ω) in (52) is defined for ω > 0 only, and the integral for the time dependence of **E**(**r**, t) in (53) extends only over positive values of ω. This is in accordance with the convention we used to define the complex field **E**(**r**, t) in (39). All other space- and time-dependent quantities, including other field vectors and the permittivity and susceptibility tensors, are transformed in a similar manner. Through the Fourier transform, the convolution integrals in real space and time become simple products in the momentum space and frequency domain. Consequently, we have

\[\tag{54}\mathbf{P}(\mathbf{k},\omega)=\epsilon_0\pmb{\chi}(\mathbf{k},\omega)\cdot\mathbf{E}(\mathbf{k},\omega)\]\[\tag{55}\mathbf{D}(\mathbf{k},\omega)=\pmb{\epsilon}(\mathbf{k},\omega)\cdot\mathbf{E}(\mathbf{k},\omega)\]

The next part continues with the linear optical susceptibility tutorial.