Green's vector identity

The second derivative of two vector functions is related to the divergence of the vector functions with first order operators. Namely,

Error al representar (SVG (MathML puede ser habilitado mediante un plugin de navegador): respuesta no válida («Math extension cannot connect to Restbase.») del servidor «https://en.wikipedia.org/api/rest_v1/»:): \mathbf{P}\cdot\nabla^{2}\mathbf{Q}-\mathbf{Q}\cdot\nabla^{2}\mathbf{P}=\nabla\cdot\left[\mathbf{P}\left(\nabla\cdot\mathbf{Q}\right)-\mathbf{Q}\left(\nabla\cdot\mathbf{P}\right)+\mathbf{P}\times\nabla\times\mathbf{Q}-\mathbf{Q}\times\nabla\times\mathbf{P}\right].

Introduction

Green’s second identity establishes a relationship between second and (the divergence of) first order derivatives of two scalar functions

$p_{m}\nabla ^{2}q_{m}-q_{m}\nabla ^{2}p_{m}=\nabla \cdot \left(p_{m}\nabla q_{m}-q_{m}\nabla p_{m}\right),$

where $p_{m}$ and $q_{m}$ are two arbitrary scalar fields. This identity is of great importance in physics because continuity equations can thus be established for scalar fields such as mass or energy .

Although the second Green’s identity is always presented in vector analysis, only a scalar version is found on textbooks. Even in the specialized literature, a vector version is not easily found. In vector diffraction theory, two versions of Green’s second identity are introduced. One variant invokes the divergence of a cross product and states a relationship in terms of the curl-curl of the field $\mathbf {P} \cdot \left(\nabla \times \nabla \times \mathbf {Q} \right)-\mathbf {Q} \cdot \left(\nabla \times \nabla \times \mathbf {P} \right)=\nabla \cdot \left(\mathbf {Q} \times \nabla \times \mathbf {P} -\mathbf {P} \times \nabla \times \mathbf {Q} \right)$. This equation can be written in terms of the Laplacians:

$\mathbf {P} \cdot \nabla ^{2}\mathbf {Q} -\mathbf {Q} \cdot \nabla ^{2}\mathbf {P} +\mathbf {Q} \cdot \left[\nabla \left(\nabla \cdot \mathbf {P} \right)\right]-\mathbf {P} \cdot \left[\nabla \left(\nabla \cdot \mathbf {Q} \right)\right]=\nabla \cdot \left(\mathbf {P} \times \nabla \times \mathbf {Q} -\mathbf {Q} \times \nabla \times \mathbf {P} \right).$

However, the terms Error al representar (SVG (MathML puede ser habilitado mediante un plugin de navegador): respuesta no válida («Math extension cannot connect to Restbase.») del servidor «https://en.wikipedia.org/api/rest_v1/»:): \mathbf{Q}\cdot\left[\nabla\left(\nabla\cdot\mathbf{P}\right)\right]-\mathbf{P}\cdot\left[\nabla\left(\nabla\cdot\mathbf{Q}\right)\right] , could not be readily written in terms of a divergence. The other approach introduces bi-vectors, this formulation requires a dyadic Green function .

Divergence of two vector fields

Consider that the scalar fields in Green's second identity are the Cartesian components of vector fields, i.e. $\mathbf {P} =\sum \limits _{m}p_{m}{\hat {\mathbf {e} }}_{m}$ and $\mathbf {Q} =\sum \limits _{m}q_{m}{\hat {\mathbf {e} }}_{m}$. Summing up the equations for each component, we obtain

$\sum \limits _{m}\left[p_{m}\nabla ^{2}q_{m}-q_{m}\nabla ^{2}p_{m}\right]=\sum \limits _{m}\left[\nabla \cdot \left(p_{m}\nabla q_{m}-q_{m}\nabla p_{m}\right)\right].$

The LHS according to the definition of the dot product may be written in vector form as

$\sum \limits _{m}\left[p_{m}\nabla ^{2}q_{m}-q_{m}\nabla ^{2}p_{m}\right]=\mathbf {P} \cdot \nabla ^{2}\mathbf {Q} -\mathbf {Q} \cdot \nabla ^{2}\mathbf {P} .$

The RHS is a bit more awkward to express in terms of vector operators. Due to the distributivity of the divergence operator over addition, the sum of the divergence is equal to the divergence of the sum, i.e. $\sum \limits _{m}\left[\nabla \cdot \left(p_{m}\nabla q_{m}-q_{m}\nabla p_{m}\right)\right]=\nabla \cdot \left(\sum \limits _{m}p_{m}\nabla q_{m}-\sum \limits _{m}q_{m}\nabla p_{m}\right)$. Recall the vector identity for the gradient of a dot product

$\nabla \left(\mathbf {P} \cdot \mathbf {Q} \right)=\left(\mathbf {P} \cdot \nabla \right)\mathbf {Q} +\left(\mathbf {Q} \cdot \nabla \right)\mathbf {P} +\mathbf {P} \times \nabla \times \mathbf {Q} +\mathbf {Q} \times \nabla \times \mathbf {P}$

which, written out in vector components is given by Error al representar (SVG (MathML puede ser habilitado mediante un plugin de navegador): respuesta no válida («Math extension cannot connect to Restbase.») del servidor «https://en.wikipedia.org/api/rest_v1/»:): \nabla\left(\mathbf{P}\cdot\mathbf{Q}\right)=\nabla\sum\limits _{m}p_{m}q_{m}=\sum\limits _{m}p_{m}\nabla q_{m}+\sum\limits _{m}q_{m}\nabla p_{m}. This result is similar to what we wish to evince in vector terms ’except’ for the minus sign. Since the differential operators in each term act either over one vector (say $p_{m}$’s) or the other (Error al representar (SVG (MathML puede ser habilitado mediante un plugin de navegador): respuesta no válida («Math extension cannot connect to Restbase.») del servidor «https://en.wikipedia.org/api/rest_v1/»:): q_{m} ’s) , the contribution to each term must be

$\sum \limits _{m}p_{m}\nabla q_{m}=\left(\mathbf {P} \cdot \nabla \right)\mathbf {Q} +\mathbf {P} \times \nabla \times \mathbf {Q} ,\sum \limits _{m}q_{m}\nabla p_{m}=\left(\mathbf {Q} \cdot \nabla \right)\mathbf {P} +\mathbf {Q} \times \nabla \times \mathbf {P} .$

These results can be rigorously proven to be correct in appendix [sec:Derivation-by-components] through evaluation of the vector components. Therefore, the RHS can be written in vector form as

$\sum \limits _{m}p_{m}\nabla q_{m}-\sum \limits _{m}q_{m}\nabla p_{m}=\left(\mathbf {P} \cdot \nabla \right)\mathbf {Q} +\mathbf {P} \times \nabla \times \mathbf {Q} -\left(\mathbf {Q} \cdot \nabla \right)\mathbf {P} -\mathbf {Q} \times \nabla \times \mathbf {P} .$

Putting together these two results, a theorem for vector fields analogous to Green’s theorem for scalar fields is obtained

$\mathbf {P} \cdot \nabla ^{2}\mathbf {Q} -\mathbf {Q} \cdot \nabla ^{2}\mathbf {P} =\nabla \cdot \left[\left(\mathbf {P} \cdot \nabla \right)\mathbf {Q} +\mathbf {P} \times \nabla \times \mathbf {Q} -\left(\mathbf {Q} \cdot \nabla \right)\mathbf {P} -\mathbf {Q} \times \nabla \times \mathbf {P} \right].$

Reassuringly, from the vector relationship ([eq:vec green]), we can go back to the scalar case as shown in appendix [sec:scalar-case]. The curl of a cross product can be written as Error al representar (SVG (MathML puede ser habilitado mediante un plugin de navegador): respuesta no válida («Math extension cannot connect to Restbase.») del servidor «https://en.wikipedia.org/api/rest_v1/»:): \nabla\times\left(\mathbf{P}\times\mathbf{Q}\right)=\left(\mathbf{Q}\cdot\nabla\right)\mathbf{P}-\left(\mathbf{P}\cdot\nabla\right)\mathbf{Q}+\mathbf{P}\left(\nabla\cdot\mathbf{Q}\right)-\mathbf{Q}\left(\nabla\cdot\mathbf{P}\right) ; Green’s vector identity can then be rewritten as

$\mathbf {P} \cdot \nabla ^{2}\mathbf {Q} -\mathbf {Q} \cdot \nabla ^{2}\mathbf {P} =\nabla \cdot \left[\mathbf {P} \left(\nabla \cdot \mathbf {Q} \right)-\mathbf {Q} \left(\nabla \cdot \mathbf {P} \right)-\nabla \times \left(\mathbf {P} \times \mathbf {Q} \right)+\mathbf {P} \times \nabla \times \mathbf {Q} -\mathbf {Q} \times \nabla \times \mathbf {P} \right].$

Since the divergence of a curl is zero, the third term vanishes and the identity can be written as

Error al representar (SVG (MathML puede ser habilitado mediante un plugin de navegador): respuesta no válida («Math extension cannot connect to Restbase.») del servidor «https://en.wikipedia.org/api/rest_v1/»:): \mathbf{P}\cdot\nabla^{2}\mathbf{Q}-\mathbf{Q}\cdot\nabla^{2}\mathbf{P}=\nabla\cdot\left[\mathbf{P}\left(\nabla\cdot\mathbf{Q}\right)-\mathbf{Q}\left(\nabla\cdot\mathbf{P}\right)+\mathbf{P}\times\nabla\times\mathbf{Q}-\mathbf{Q}\times\nabla\times\mathbf{P}\right].

This result should prove useful when the divergence and curl of the fields can be established in terms of other quantities. Another identity that may prove useful is obtained from the divergence of ([eq: grad prod]) $\nabla \cdot \nabla \left(\mathbf {P} \cdot \mathbf {Q} \right)=\nabla \cdot \left[\left(\mathbf {P} \cdot \nabla \right)\mathbf {Q} +\left(\mathbf {Q} \cdot \nabla \right)\mathbf {P} +\mathbf {P} \times \nabla \times \mathbf {Q} +\mathbf {Q} \times \nabla \times \mathbf {P} \right],$ invoking the Green’s vector identity ([eq:vec green]) derived above, the Laplacian of the dot product can be expressed in terms of the Laplacians of the factors

$\nabla ^{2}\left(\mathbf {P} \cdot \mathbf {Q} \right)=\mathbf {P} \cdot \nabla ^{2}\mathbf {Q} -\mathbf {Q} \cdot \nabla ^{2}\mathbf {P} +2\nabla \cdot \left[\left(\mathbf {Q} \cdot \nabla \right)\mathbf {P} +\mathbf {Q} \times \nabla \times \mathbf {P} \right].$

If the substitution of the vector identity ([eq:vec green]) is performed eliminating the terms $\left(\mathbf {Q} \cdot \nabla \right)\mathbf {P} +\mathbf {Q} \times \nabla \times \mathbf {P}$, the Laplacian of the dot product is

$\nabla ^{2}\left(\mathbf {P} \cdot \mathbf {Q} \right)=-\left(\mathbf {P} \cdot \nabla ^{2}\mathbf {Q} -\mathbf {Q} \cdot \nabla ^{2}\mathbf {P} \right)+2\nabla \cdot \left[\left(\mathbf {P} \cdot \nabla \right)\mathbf {Q} +\mathbf {P} \times \nabla \times \mathbf {Q} \right].$

Conclusions

Green’s second identity relating the Laplacians with the divergence has been derived for vector fields. No use of bi-vectors or dyadics has been made as in some previous approaches. In diffraction theory, the vector identity was stated before in terms of the curl. However, this earlier formulation had the drawback that the Laplacian could not be invoked without involving extra terms. As a corollary, the awkward terms in eq:vec difrac ident can now be written in terms of a divergence by comparison with eq:vec green 2

$\mathbf {P} \cdot \left[\nabla \left(\nabla \cdot \mathbf {Q} \right)\right]-\mathbf {Q} \cdot \left[\nabla \left(\nabla \cdot \mathbf {P} \right)\right]=\nabla \cdot \left[\mathbf {P} \left(\nabla \cdot \mathbf {Q} \right)-\mathbf {Q} \left(\nabla \cdot \mathbf {P} \right)\right].$

This result can be verified by expanding the divergence of a scalar times a vector on the RHS.

The condition imposed by Helmholtz equation $\nabla ^{2}\mathbf {P} =-k^{2}\mathbf {P}$ can be readily incorporated in the present formulation of Green’s second identity. This result is particularly useful if the vector fields satisfy the wave equation.

Derivation by components

In order to evaluate $\left(\mathbf {P} \cdot \nabla \right)\mathbf {Q} +\mathbf {P} \times \nabla \times \mathbf {Q}$ consider the first term in three dimensional Cartesian components Error al representar (SVG (MathML puede ser habilitado mediante un plugin de navegador): respuesta no válida («Math extension cannot connect to Restbase.») del servidor «https://en.wikipedia.org/api/rest_v1/»:): \left(\mathbf{P}\cdot\nabla\right)\mathbf{Q}=\left(p_{x}\frac{\partial}{\partial x}+p_{y}\frac{\partial}{\partial y}+p_{z}\frac{\partial}{\partial z}\right)\left(q_{x}\hat{\mathbf{e}}_{x}+q_{y}\hat{\mathbf{e}}_{y}+q_{z}\hat{\mathbf{e}}_{z}\right) that may be written as

Error al representar (SVG (MathML puede ser habilitado mediante un plugin de navegador): respuesta no válida («Math extension cannot connect to Restbase.») del servidor «https://en.wikipedia.org/api/rest_v1/»:): \left(\mathbf{P}\cdot\nabla\right)\mathbf{Q} = \left(p_{x}\frac{\partial q_{x}}{\partial x}+p_{y}\frac{\partial q_{x}}{\partial y}+p_{z}\frac{\partial q_{x}}{\partial z}\right)\hat{\mathbf{e}}_{x} + \left(p_{x}\frac{\partial q_{y}}{\partial x}+p_{y}\frac{\partial q_{y}}{\partial y}+p_{z}\frac{\partial q_{y}}{\partial z}\right)\hat{\mathbf{e}}_{y} + \left(p_{x}\frac{\partial q_{z}}{\partial x}+p_{y}\frac{\partial q_{z}}{\partial y}+p_{z}\frac{\partial q_{z}}{\partial z}\right)\hat{\mathbf{e}}_{z}.

The curl in the second term is $^{R}\mathbf {Q} =\nabla \times \mathbf {Q} =\left({\frac {\partial q_{z}}{\partial y}}-{\frac {\partial q_{y}}{\partial z}}\right){\hat {\mathbf {e} }}_{x}+\left({\frac {\partial q_{x}}{\partial z}}-{\frac {\partial q_{z}}{\partial x}}\right){\hat {\mathbf {e} }}_{y}+\left({\frac {\partial q_{y}}{\partial x}}-{\frac {\partial q_{x}}{\partial y}}\right){\hat {\mathbf {e} }}_{z}.$ The cross product is $\mathbf {P} \times {}^{R}\mathbf {Q} =\left(p_{y}{}^{R}q_{z}-p_{z}{}^{R}q_{y}\right){\hat {\mathbf {e} }}_{x}+\left(p_{z}{}^{R}q_{x}-p_{x}{}^{R}q_{z}\right){\hat {\mathbf {e} }}_{y}+\left(p_{x}{}^{R}q_{y}-p_{y}{}^{R}q_{x}\right){\hat {\mathbf {e} }}_{z};$ The second term is then

$\mathbf {P} \times \nabla \times \mathbf {Q} =\left(p_{y}\left({\frac {\partial q_{y}}{\partial x}}-{\frac {\partial q_{x}}{\partial y}}\right)-p_{z}\left({\frac {\partial q_{x}}{\partial z}}-{\frac {\partial q_{z}}{\partial x}}\right)\right){\hat {\mathbf {e} }}_{x}+\left(p_{z}\left({\frac {\partial q_{z}}{\partial y}}-{\frac {\partial q_{y}}{\partial z}}\right)-p_{x}\left({\frac {\partial q_{y}}{\partial x}}-{\frac {\partial q_{x}}{\partial y}}\right)\right){\hat {\mathbf {e} }}_{y}+\left(p_{x}\left({\frac {\partial q_{x}}{\partial z}}-{\frac {\partial q_{z}}{\partial x}}\right)-p_{y}\left({\frac {\partial q_{z}}{\partial y}}-{\frac {\partial q_{y}}{\partial z}}\right)\right){\hat {\mathbf {e} }}_{z}$

that expands to

$\mathbf {P} \times \nabla \times \mathbf {Q} =\left(p_{y}{\frac {\partial q_{y}}{\partial x}}-p_{y}{\frac {\partial q_{x}}{\partial y}}-p_{z}{\frac {\partial q_{x}}{\partial z}}+p_{z}{\frac {\partial q_{z}}{\partial x}}\right){\hat {\mathbf {e} }}_{x}+\left(p_{z}{\frac {\partial q_{z}}{\partial y}}-p_{z}{\frac {\partial q_{y}}{\partial z}}-p_{x}{\frac {\partial q_{y}}{\partial x}}+p_{x}{\frac {\partial q_{x}}{\partial y}}\right){\hat {\mathbf {e} }}_{y}+\left(p_{x}{\frac {\partial q_{x}}{\partial z}}-p_{x}{\frac {\partial q_{z}}{\partial x}}-p_{y}{\frac {\partial q_{z}}{\partial y}}+p_{y}{\frac {\partial q_{y}}{\partial z}}\right){\hat {\mathbf {e} }}_{z}.$

Evaluate Error al representar (SVG (MathML puede ser habilitado mediante un plugin de navegador): respuesta no válida («Math extension cannot connect to Restbase.») del servidor «https://en.wikipedia.org/api/rest_v1/»:): \left(\mathbf{P}\cdot\nabla\right)\mathbf{Q}+\mathbf{P}\times\nabla\times\mathbf{Q} in the x direction

$\left[\left(\mathbf {P} \cdot \nabla \right)\mathbf {Q} +\mathbf {P} \times \nabla \times \mathbf {Q} \right]_{x}=\left(p_{x}{\frac {\partial q_{x}}{\partial x}}+p_{y}{\frac {\partial q_{x}}{\partial y}}+p_{z}{\frac {\partial q_{x}}{\partial z}}+p_{y}{\frac {\partial q_{y}}{\partial x}}-p_{y}{\frac {\partial q_{x}}{\partial y}}-p_{z}{\frac {\partial q_{x}}{\partial z}}+p_{z}{\frac {\partial q_{z}}{\partial x}}\right){\hat {\mathbf {e} }}_{x},$

canceling out terms $\left[\left(\mathbf {P} \cdot \nabla \right)\mathbf {Q} +\mathbf {P} \times \nabla \times \mathbf {Q} \right]_{x}=\left(p_{x}{\frac {\partial q_{x}}{\partial x}}+p_{y}{\frac {\partial q_{y}}{\partial x}}+p_{z}{\frac {\partial q_{z}}{\partial x}}\right){\hat {\mathbf {e} }}_{x}.$

Analogous results are obtained in the other directions so that

$\left(\mathbf {P} \cdot \nabla \right)\mathbf {Q} +\mathbf {P} \times \nabla \times \mathbf {Q} =\left(p_{x}{\frac {\partial q_{x}}{\partial x}}+p_{y}{\frac {\partial q_{y}}{\partial x}}+p_{z}{\frac {\partial q_{z}}{\partial x}}\right){\hat {\mathbf {e} }}_{x}+\left(p_{x}{\frac {\partial q_{x}}{\partial y}}+p_{y}{\frac {\partial q_{y}}{\partial y}}+p_{z}{\frac {\partial q_{z}}{\partial y}}\right){\hat {\mathbf {e} }}_{y}+\left(p_{x}{\frac {\partial q_{x}}{\partial z}}+p_{y}{\frac {\partial q_{y}}{\partial z}}+p_{z}{\frac {\partial q_{z}}{\partial z}}\right){\hat {\mathbf {e} }}_{z}.$

that may be written out in vector form as $\left(\mathbf {P} \cdot \nabla \right)\mathbf {Q} +\mathbf {P} \times \nabla \times \mathbf {Q} =\mathbf {P} \cdot {\frac {\partial \mathbf {Q} }{\partial x}}+\mathbf {P} \cdot {\frac {\partial \mathbf {Q} }{\partial y}}+\mathbf {P} \cdot {\frac {\partial \mathbf {Q} }{\partial z}}.$ However, the terms can be rearranged as

$\left(\mathbf {P} \cdot \nabla \right)\mathbf {Q} +\mathbf {P} \times \nabla \times \mathbf {Q} =\left(p_{x}{\frac {\partial q_{x}}{\partial x}}{\hat {\mathbf {e} }}_{x}+p_{x}{\frac {\partial q_{x}}{\partial y}}{\hat {\mathbf {e} }}_{y}+p_{x}{\frac {\partial q_{x}}{\partial z}}{\hat {\mathbf {e} }}_{z}\right)+\left(p_{y}{\frac {\partial q_{y}}{\partial x}}{\hat {\mathbf {e} }}_{x}+p_{y}{\frac {\partial q_{y}}{\partial y}}{\hat {\mathbf {e} }}_{y}+p_{y}{\frac {\partial q_{y}}{\partial z}}{\hat {\mathbf {e} }}_{z}\right)+\left(p_{z}{\frac {\partial q_{z}}{\partial x}}{\hat {\mathbf {e} }}_{x}+p_{z}{\frac {\partial q_{z}}{\partial y}}{\hat {\mathbf {e} }}_{y}+p_{z}{\frac {\partial q_{z}}{\partial z}}{\hat {\mathbf {e} }}_{z}\right),$

and thus $\left(\mathbf {P} \cdot \nabla \right)\mathbf {Q} +\mathbf {P} \times \nabla \times \mathbf {Q} =\left(p_{x}\nabla q_{x}+p_{y}\nabla q_{y}+p_{z}\nabla q_{z}\right).$ An equivalent procedure for Error al representar (SVG (MathML puede ser habilitado mediante un plugin de navegador): respuesta no válida («Math extension cannot connect to Restbase.») del servidor «https://en.wikipedia.org/api/rest_v1/»:): \left(\mathbf{Q}\cdot\nabla\right)\mathbf{P}+\mathbf{Q}\times\nabla\times\mathbf{P} gives $\left(\mathbf {Q} \cdot \nabla \right)\mathbf {P} +\mathbf {Q} \times \nabla \times \mathbf {P} =q_{x}\nabla p_{x}+q_{y}\nabla p_{y}+q_{z}\nabla p_{z}.$

[sec:scalar-case]scalar case

If we take one component vectors, for example, $\mathbf {P} \rightarrow p_{x}{\hat {\mathbf {e} }}_{x},\mathbf {Q} \rightarrow q_{x}{\hat {\mathbf {e} }}_{x}$, the vector relationship ([eq:vec green]) becomes $p_{x}\nabla ^{2}q_{x}-q_{x}\nabla ^{2}p_{x}=\nabla \cdot \left[p_{x}{\frac {\partial q_{x}}{\partial x}}{\hat {\mathbf {e} }}_{x}-q_{x}{\frac {\partial p_{x}}{\partial x}}{\hat {\mathbf {e} }}_{x}+\mathbf {P} \times \nabla \times \mathbf {Q} -\mathbf {Q} \times \nabla \times \mathbf {P} \right].$ Since $\nabla \times \mathbf {Q} ={\frac {\partial q_{x}}{\partial z}}{\hat {\mathbf {e} }}_{y}-{\frac {\partial q_{x}}{\partial y}}{\hat {\mathbf {e} }}_{z}$,

Error al representar (SVG (MathML puede ser habilitado mediante un plugin de navegador): respuesta no válida («Math extension cannot connect to Restbase.») del servidor «https://en.wikipedia.org/api/rest_v1/»:): \mathbf{P}\times\nabla\times\mathbf{Q} = \left(p_{x}\frac{\partial q_{x}}{\partial y}\right)\hat{\mathbf{e}}_{y}+\left(p_{x}\frac{\partial q_{x}}{\partial z}\right)\hat{\mathbf{e}}_{z}

and $\mathbf {Q} \times \nabla \times \mathbf {P} =\left(q_{x}{\frac {\partial p_{x}}{\partial y}}\right){\hat {\mathbf {e} }}_{y}+\left(q_{x}{\frac {\partial p_{x}}{\partial z}}\right){\hat {\mathbf {e} }}_{z}$. Therefore Error al representar (SVG (MathML puede ser habilitado mediante un plugin de navegador): respuesta no válida («Math extension cannot connect to Restbase.») del servidor «https://en.wikipedia.org/api/rest_v1/»:): p_{x}\nabla^{2}q_{x}-q_{x}\nabla^{2}p_{x}=\nabla\cdot\left[p_{x}\nabla q_{x}-q_{x}\nabla p_{x}\right], we recover Green’s second identity for the functions $p_{x},\,q_{x}$.

Colegrave88 R. K. Colegrave and M. A. Mannan. Invariants for the time-dependent harmonic oscillator. J. Math. Phys., 29(7):1580–1587, 1988.

Fernandezguasti02 M. Fernández-Guasti and A. Gil-Villegas. Orthogonal functions invariant for the time dependent harmonic oscillator. Physics Letters A, 292:243–245., 2002.

Fernandezguasti03a M. Fernández-Guasti and H. Moya-Cessa. Amplitude and phase representation of quantum invariants for the time dependent harmonic oscillator. Physical Review A, 67(063803):063803–1–5, 2003.

Pedrosa04 I. A. Pedrosa and I. Guedes. Quantum states of a generalized time-dependent inverted harmonic oscillator. Int. J. of Mod. Phys. B, 18(9):1379–1385, 2004.

Lewis67 H. R. Lewis. Classical and quantum systems with time-dependent harmonic-oscillator-type Hamiltonians. Phys. Rev. Lett., 18(13):510–512, 1967.

love1901 A. E. H. Love. The Integration of the Equations of Propagation of Electric Waves. Philosophical Transactions of the Royal Society of London. Series A, Containing Papers of a Mathematical or Physical Character, 197:pp. 1–45, 1901.

stratton1939 J. A. Stratton and L. J. Chu. Diffraction Theory of Electromagnetic Waves. Phys. Rev., 56(1):99–107, Jul 1939.

bruce2010 N. C. Bruce. Double scatter vector-wave Kirchhoff scattering from perfectly conducting surfaces with infinite slopes. Journal of Optics, 12(8):085701, 2010.

franz1950 W Franz. On the Theory of Diffraction. Proceedings of the Physical Society. Section A, 63(9):925, 1950.

tai1972 Chen-To Tai. Kirchhoff theory: Scalar, vector, or dyadic? Antennas and Propagation, IEEE Transactions on, 20(1):114–115, jan 1972.

Arfken70 G. Arfken. Mathematical Methods for Physicists. Academic Press, NY, 1970.