research on quantum statistical properties of light has been a subject of intense investigations and discussions . in view of the very wide - spread potential applications ,
the influence of material bodies on the quantum features of light must be thoroughly investigated .
typical quantum effects that are closely related to the change of the quantum vacuum due to the presence of material bodies are the casimir effect @xcite ( for a review , see @xcite ) and the purcell effect @xcite ( see also @xcite ) . though there has been a large body of work on the problem of quantization of the electromagnetic field in dielectric media ( for a review , see @xcite ) , most work has been concentrated on nonabsorbing media .
roughly speaking , there have been two routes to treat radiation in absorbing media . in the first
, attention has been restricted to equilibrium field correlation functions that are calculated by employing the dissipation - fluctuation theorem @xcite . in the second
, explicit field quantization has been performed ( for a review , see @xcite ) .
a consistent quantization scheme for the electromagnetic field in bulk material has been given by huttner and barnett @xcite .
using the mesoscopic hopfield model @xcite for a dielectric , they have diagonalized the bilinear hamiltonian of the system that consists of the electromagnetic field , a harmonic - oscillator polarization field , and an infinite set of harmonic - oscillator reservoir fields responsible for absorption . by starting directly with the phenomenological maxwell equations for the macroscopic electromagnetic field , the scheme can naturally be extended to arbitrary inhomogeneous media @xcite characterized by a space- and frequency - dependent complex permittivity that satisfies the kramers - kronig relations .
the theory is based on a source - quantity representation of the electromagnetic field , where the field operators are expressed in terms of a continuous set of fundamental bose fields via the green tensor of the classical problem ( for details , see also @xcite ) .
the green - tensor formalism is well suited to study the behavior of the quantized electromagnetic field in the presence of dispersing and absorbing bodies .
in particular , it has successfully been applied to the study of input - output relations @xcite , the spontaneous emission @xcite , the resonant energy exchange @xcite , the resonant dipole - dipole interaction @xcite , and the casimir effect @xcite , and also various geometries ( including multilayer structures ) have been considered .
it is worth noting that the theory not only takes into account material absorption in a consistent way but it also includes automatically evanescent - field contributions generated , e.g. , by radiating atoms that are very close to the bodies under consideration .    for the particular case of a dispersing and absorbing dielectric structure in vacuum , without active light sources at any finite distance from the structure ( and without active light sources inside the structure )
, a ratherinvolved hybrid quantization scheme has been proposed @xcite and used to study the problem of three - dimensional input - output relations @xcite . in the scheme
, the whole space is divided into two regions , namely a dielectric scattering region and a vacuum region that surrounds the scattering region .
the quantization in the vacuum region is performed on the basis of a mode decomposition of the electromagnetic - field operators in a similar way as in vacuum quantum electrodynamics , whereas the contribution of the dielectric scattering region is taken into account by applying the green tensor formalism @xcite .
the argument is @xcite that the green tensor formalism alone is not complete , because it does not explicitly treat the boundary of the medium with free space and thus does not describe explicitly the scattering and emission processes which are usually experimentally investigated .
this is of course not the case as the above mentioned applications @xcite and the following study of the three - dimensional input - output relations at multilayer plates clearly show .
moreover , when there are active light sources at _ finite distances _ from the scattering region , then the incoming field contains both propagating and evanescent components .
the latter ones , which are typically observed for small distances , are not included in the hybrid quantization scheme @xcite and the input - output relations derived from it @xcite .
however , they are automatically included in the green tensor formalism @xcite .    in the following we consider the problem of the three - dimensional input - output relations for the electromagnetic - field operators at the boundaries of a dielectric planar multilayer structure in more detail , by applying the green tensor formalism and extending our previous one - dimensional analysis @xcite to three dimensions . to be quite general ,
we allow ( i ) for the presence of active light sources at arbitrary positions inside and/or outside the multilayer plate and thus ( ii ) for both propagating - field components and evanescent - field components .
in particular , in cavity qed the active sources are typically inside the plate , and the outgoing fields are not only determined with the incoming fields and the noise fields ( unavoidably associated with material absorption ) but also with the fields generated by the sources inside the plate . on the contrary , in optical near - field microscopy the active ( probe )
sources are typically outside the plate , but near its surface .    employing the well - known three - dimensional green tensor for a multilayer plate @xcite
, we introduce amplitude operators for the input and output fields as well as for the fields inside the plate and derive input - output relations both in the two - dimensional fourier space and in the coordinate space .
finally , the problem of introduction of bosonic input and output operators in the two - dimensional fourier space is considered in detail .
the paper is organized as follows . in sec .
[ sec : ba ] the quantization scheme is briefly summarized and the green tensor for a multilayer plate is introduced . in sec .
[ sec : gr ] the generally valid input - output relations are derived , and the problem of formulating input - output relations for bosonic field operators is studied .
a summary and some concluding remarks are given in sec .
[ sec : co ] followed by an appendix , in which relevant commutation relations are derived .
to study optical fields interacting with active sources in the presence of dispersing and absorbing ( linear ) dielectric bodies , we first note that on a length scale that is large compared with interatomic distances in the bodies , the effect of the bodies can be described within the frame of macroscopic maxwell equations in terms of a spatially varying permittivity which is a complex function of frequency .
this concept , which is widely used in classical optics also applies in quantum optics .
let us assume that the active light sources are neutral atoms and consider an arbitrarily inhomogeneous medium characterized by a permittivity @xmath1 where the real part @xmath2 and the imaginary part @xmath3 are necessarily related to each other through the kramers - kronig relations , due to the causality principle .
the motion of the atoms and the medium - assisted electromagnetic field is then governed by the multipolar - coupling hamiltonian @xcite    @xmath4        \!\!\bigg\}^2 \nonumber \\[1ex ] & &        + \,\frac{1}{2\epsilon_0 }        \sum_{a , a ' } \int { \mbox{\rm d}}^3{r } \ ,        \hat{\bf p}_{a}({\bf r})\cdot \hat{\bf p}_{a'}({\bf r } )        -\sum_a \int { \mbox{\rm d}}^3{r}\ , \hat{{\bf p}}_{a}({\bf r})\cdot        \hat{\bf e}({\bf r } ) ,        \end{aligned}\ ] ]    where @xmath5 numbers the atoms , and @xmath6 numbers the charged particles inside the @xmath5-th atom .
further , @xmath7 and @xmath8 are respectively the operators of coordinates and canonical momenta of the particles , and @xmath9\ ] ] is the operator of the polarization of the @xmath5-th atom at position @xmath10 .
the @xmath11 [ and @xmath12 are bosonic field operators that play the role of fundamental variables of the electromagnetic field and the medium , including a reservoir necessarily associated with material absorption .
the commutation relations for the operators @xmath13 and @xmath14 are @xmath15        = \delta _ { \mu \mu'}\delta ( \omega - \omega   ' )        \delta ^{(3)}({\bf r } - { \bf r } ' ) , \ ] ] @xmath16 = 0,\ ] ] where the greek letters label the cartesian coordinates @xmath17 . in eq .
( [ 2 ] ) , the operators @xmath18 and @xmath19 of the medium - assisted electromagnetic field are expressed in terms of the fundamental variables as follows : @xmath20 @xmath21 @xmath22 @xmath23 where the integration should be performed over all space , and @xmath24 is the classical green tensor , which can be found from the equation @xmath25 together with appropriate boundary conditions at infinity . equations ( [ 6 ] - [ 9 ] ) can be considered as generalization of the familiar mode decomposition that would apply if dispersion and absorption could be disregarded . instead of dealing with equations of motion for mode operators , equations of motion for the fields @xmath11 must be treated .
it should be pointed out that since the real part @xmath2 and the imaginary part @xmath3 of the permittivity are related to each other through the kramers - kronig relations , @xmath3 can not vanish identically for really existing media .
clearly , @xmath26 can be very small , so that @xmath27 in eq .  ( [ 9 ] ) is very small in certain areas of space ( @xmath28 ) .
however , this statement says nothing about the magnitude of the total integral over the coordinate @xmath28 , because of the following integral relation for the classical green tensor @xmath29 ,        \end{aligned}\ ] ] where we have adopted the convention of summation over repeated vector - component indices .
in particular , this relation enables one to include also vacuum - like areas in the consideration .
thus , all the calculations are to be performed by assuming a permittivity close to unity with a small but finite imaginary part in those areas , and at the end the permittivity may be set equal to unity . in practice ,
experimental realization of ( macroscopic ) vacuum areas is of course fictional .
scheme of the multilayer dielectric plate .
the hatched regions indicate the presence of active light sources .
, width=326 ]      as mentioned , the quantization scheme is valid for an arbitrary space dependence of the permittivity . here
, we consider a multilayered planar structure ( fig .  [ fig ] ) , whose permittivity is defined in a stepwise fashion ( the @xmath30-direction is chosen to be perpendicular to the layers ) : @xmath31 where @xmath32                  0 , \ , & \mbox{otherwise } ,          \end{array } \right.\ ] ] and @xmath33 is the ( complex ) permittivity of the @xmath34-th layer . in the above
, the index @xmath34 labels the region on the left of the plate , the region on the right of the plate , and the layers of the plate . for simplicity , we express the @xmath30-coordinate dependence in shifted coordinate systems , introduced in each layer separately , so that the range of the @xmath30-coordinate is taken to be for the region on the left of the plate , for the region on the right of the plate , and for the @xmath34-th layer of the plate with thickness @xmath35 .    exploiting the translational symmetry in the lane , we may represent the green tensor as a two - dimensional fourier integral @xmath36 where , and is the wave vector parallel to the layers .
the notations @xmath37 and @xmath38 indicate that @xmath30 varies in the @xmath34-th layer and @xmath39 in the @xmath40-th layer . inserting eq .
( [ 14 ] ) into eq .
( [ 9 ] ) , we may write the electric - field operator @xmath41 as a two - fold fourier transform , @xmath42 where @xmath43        & & \times \sum _ { j'=0}^n \int _ { [ j ' ] } { \mbox{\rm d}}z'\ ,        { \mbox{\boldmath $ g$}}^{(jj')}(z , z ' , { \bf k   } , \omega)\cdot        \hat { \bf j}^{(j ' ) } ( z ' , { \bf k   } , \omega )        \end{aligned}\ ] ] ( @xmath44 $ ] indicates integration over the @xmath40-th region ) . here , the fourier transformed current operators @xmath45 have been introduced , where the operators @xmath46 are the fourier transforms of the original bosonic field operators @xmath47 , @xmath48    the green tensor @xmath49 taken from ref .
@xcite may be written as @xmath50 & &        = - \,{\bf e}_z \,\frac { \delta_{jj ' } } { k_j ^2 } \,{\bf e } _
z\ ,        \delta ( z - z ' ) +    { \mbox{\boldmath $ g$}}^{(jj')}(z , z ' , { \bf k   } ,   \omega ) , \end{aligned}\ ] ] where @xmath51 & &        =        \frac { i } { 2 } \sum _ { q = p , s}\sigma _ q \left [        \bm { \mathcal{e } } ^{j > } _ { q }   ( z , { \bf k   } ,   \omega ) \ ,
\xi^{jj ' } _ { q }   \bm{\mathcal{e } } ^{j ' < } _ { q }   ( z ' , -{\bf k   } ,   \omega )        \theta ( j\!-\!j ' )        \right .
\nonumber \\[1ex ] & &         \left .
\,+\ ,        \bm{\mathcal{e}}^{j < } _ { q }   ( z , { \bf k   } ,
\omega ) \ ,        \xi^{j'j } _ { q } \bm{\mathcal{e } } ^{j ' > } _ { q } ( z ' , -{\bf k   } ,   \omega )        \theta ( j'-j )        \right ]         \end{aligned}\ ] ] ( , ) .
note that for one should write instead of [ and instead of ] .
( [ 20 ] ) , the functions @xmath52 and @xmath53 denote waves of unit strength traveling , respectively , rightward and leftward in the @xmath34-th layer and being reflected at the boundary , @xmath54 @xmath55 and @xmath56 where @xmath57 . here , @xmath58
( @xmath59 @xmath60 @xmath61 ) , where @xmath62 and and @xmath63 are , respectively , the transmission and reflection coefficients between the layers @xmath40 and @xmath34 . finally , the unit vectors @xmath64 in eqs .
( [ 21 ] ) and ( [ 22 ] ) are the polarization unit vectors for te and tm waves , @xmath65        \label{28 }        { \bf e}^{(j)}_{p\pm}({\bf k } )        & = &        \frac { 1 } { k_j } \left(\mp \beta _ j\frac{{\bf k}}{k }        + k { \bf e } _ z \right ) .
\end{aligned}\ ] ] the ` propagation constant ' @xmath66 determines the propagation behavior in @xmath30-direction of the waves in the @xmath34-th region .
note that in case of vacuum the waves are propagating only for ( i.e. , ) .
they are evanescent for ( i.e. , ) .
the equations of motion which follow from the hamiltonian ( [ 2 ] ) in sec .
[ sec : ba.1 ] together with the green tensor as given in sec .
[ sec : ba.2 ] can be used to study various effects of the atom - light interaction in the presence of dispersing and absorbing multilayered planar structures ( see the examples mentioned in sec .
[ sec1 ] ) . here , we restrict our attention to the general relation between the fields outside and inside a multilayer plate regardless of the generation of the fields by the active light sources .
the theory thus includes the general case , where sources are present in the both regions , outside and inside the plate .
if the incoming fields , incident on the two boundary planes of the plate ( i.e. , @xmath30 @xmath60 @xmath67 for and @xmath30 @xmath60 @xmath68 for ; cf . fig .  [ fig ] ) , are known as well as the fields generated inside the plate , one can calculate the fields outgoing from the two boundary planes by means of input - output relations .
note that any two planes @xmath30
@xmath60 @xmath69 @xmath70 @xmath71 and @xmath30 @xmath60 @xmath72 @xmath73 @xmath74 for and , respectively , also can be used in principle . as we shall see , these input - output relations are valid for both passive and active devices , for arbitrary layer materials and arbitrary media surrounding the plate and for both propagating - field components and evanescent - field components .      here and in the following
we restrict our attention to the electric field , noting that the corresponding expressions for the magnetic field can readily be obtained by applying eq .
( [ 8 ] ) . following the line suggested in ref .
@xcite and substituting the green tensor ( [ 19 ] ) for and into eq .
( [ 16 ] ) , we may decompose the field operators @xmath76 and @xmath77 in the form of @xmath78        & & \hspace{15ex }        \left .
+ \ , { \bf e}^{(0)}_{q- } ( { \bf k   } )        \hat{\mathrm{e}}_{q\,{\rm out } } ^{(0 ) } ( z,{\bf k   } , \omega)\right],\end{aligned}\ ] ] @xmath79        & & \hspace{15ex }        \left .        + \ , { \bf e}^{(n)}_{q+ } ( { \bf k   } )        \hat{\mathrm{e}}_{q\,{\rm out}}^{(n ) } ( z,{\bf k } , \omega)\right ] .        \end{aligned}\ ] ] here , the operators @xmath80        & & \times\ ;        \int _ { -\infty } ^{z}\ ! \ ! { \mbox{\rm d}}z'\ , e^{- i \beta _ 0 z ' }        \,\hat { \bf j } ^{(0 ) } ( z ' , { \bf k } , \omega )        \cdot { \bf e}^{(0)}_{q+}({\bf k } ) ,    \\[1ex ]        \label{32 } \lefteqn {        \hat{\mathrm{e}}_{q\,{\rm in}}^{(n ) } ( z , { \bf k   } , \omega )        = - \frac { \mu _ 0 \omega } { 2 \beta _ n }
e^{- i \beta _ n z } } \nonumber \\[1ex ]        & & \times\ ;        \int _ { z } ^{\infty }   { \mbox{\rm d}}z'\ ,        e^ { i \beta _ n z ' }        \,\hat { \bf j } ^{(n)}(z ' , { \bf k   } , \omega )        \cdot   { \bf e}^{(n)}_{q-}({\bf
k } )         \end{aligned}\ ] ] and @xmath81        & &        + \,e^ { - i \beta _ 0 z }        \int _ { z } ^{0 }   { \mbox{\rm d}}z'\,e^ { i \beta _ 0 z ' }        \,\hat { \bf j } ^{(0 ) } ( z ' , { \bf k } , \omega ) \cdot        { \bf e}^{(0)}_{q-}({\bf k } ) , \quad \\[1ex ]        \label{34 } \lefteqn {        \hat{\mathrm{e}}_{q\,{\rm out}}^{(n ) } ( z,{\bf k } , \omega )        = e^ { i \beta _ n z }        \hat{\mathrm{e}}_{q\,{\rm out}}^{(n ) } ( { \bf k } , \omega ) } \nonumber\\[1ex ]        & &        + \,e^ { i \beta _ n z }        \int _ { 0 } ^{z }   { \mbox{\rm d}}z'\ ,        e^{- i \beta _ n z ' }        \,\hat { \bf j } ^{(n ) } ( z ' , { \bf k } , \omega)\cdot        { \bf e}^{(n)}_{q+}({\bf k } ) , \qquad        \end{aligned}\ ] ] respectively , play the role of input and output amplitude operators , where the input - output relations    @xmath82        \hat{\mathrm{e}}_{q\,{\rm out}}^{(n)}({\bf k   } , \omega )        \end{array }
\right )        = \left (        \begin{array}{cc }        r^q_{0/n}({\bf k},\omega ) & t_{n/0 } ^q({\bf k},\omega)\\[1ex ]        t_{0/n } ^q({\bf k},\omega ) &   r ^q _ { n/0}({\bf k},\omega )        \end{array } \right )        \left (        \begin{array}{c }        \hat{\mathrm{e}}_{q\,{\rm in } } ^{(0)}({\bf k   } ,   \omega ) \\[1ex ]        \hat{\mathrm{e}}_{q\,{\rm in}}^{(n)}({\bf k   } ,   \omega )        \end{array } \right )        + \sum _ { j = 1 } ^{n-1 } \left (        \begin{array}{cc }        \phi _ { q\ , 0 + } ^{(j)}({\bf k},\omega )        & \phi _ { q\ , 0- } ^{(j)}({\bf k},\omega ) \\[1ex ]        \phi _ { q\ , n+ } ^{(j)}({\bf k},\omega )        & \phi _ { q\ ,
n- } ^{(j)}({\bf k},\omega )        \end{array } \right )        \left (        \begin{array}{c }        \hat{\mathrm{e}}_{q+}^{(j)}({\bf k   } , \omega)\\[1ex ]        \hat{\mathrm{e}}_{q-}^{(j)}({\bf k   } , \omega )        \end{array } \right)\!,\quad        \end{aligned}\ ] ]    are valid ( for the commutation relations , see the appendix ) .
they relate the output amplitude operators at the boundary planes of the multilayer plate to the input amplitude operators at the boundary planes , @xmath83        \label{37 }        \hat{\mathrm{e}}_{q\,{\rm in , out}}^{(n ) } ( { \bf k   } , \omega )        & = &        \left .
\hat{\mathrm{e}}_{q\,{\rm in , out}}^{(n ) } ( z , { \bf k   } , \omega )        \right|_{z=0^+ } ,        \end{aligned}\ ] ] and the intraplate amplitude operators @xmath84        & & \times\ ;        \int _ { 0 } ^{d_j }   { \mbox{\rm d}}z ' \ ,        e^{\mp i \beta _ j z ' }        \,\hat { \bf j}^{(j ) } ( z ' , { \bf k   } , \omega)\cdot        { \bf e}^{(j ) } _ { q\pm}({\bf k})\end{aligned}\ ] ] ( @xmath34 @xmath60 @xmath85 ) , which are associated with the excitation inside the layers of the plate . in eq .
( [ 35 ] ) , the @xmath86-coefficients read @xmath87 @xmath88 it should be pointed out that the first term on the right - hand side in eq .
( [ 19 ] ) , which gives rise to a local contribution to the electric field , has been omitted in eqs .
( [ 29 ] ) and ( [ 30 ] ) .
though this contribution is irrelevant for the incoming and outgoing fields , it must be included in the overall field operator in general , even if there are effectively no sources at the points of observations .
it is not difficult to prove that ( similar to the one - dimensional case @xcite ) the @xmath30-dependent amplitude operators ( [ 31 ] )  ( [ 34 ] ) obey quantum langevin equations , @xmath89        & & \hspace{10ex }        -\,\frac { \mu _ 0 \omega } { 2 \beta _ 0 }        \,\hat { \bf j } ^{(0 ) } ( z , { \bf k } , \omega)\cdot        { \bf e}^{(0)}_{q+}({\bf k } ) ,         \end{aligned}\ ] ] @xmath90        & & \hspace{10ex }        + \,\frac { \mu _ 0 \omega } { 2 \beta _ 0 }        \,\hat { \bf j } ^{(0 ) } ( z , { \bf k } , \omega)\cdot        { \bf e}^{(0)}_{q-}({\bf k } ) ,         \end{aligned}\ ] ] and similar equations are valid for @xmath91 and @xmath92 . these equations together with eq .
( [ 35 ] ) render it possible to easily calculate the input and output fields at any position outside the plate .
needless to say that other than the boundary planes and of the plate can be chosen as reference planes for formulating the input - output relations .
equation ( [ 35 ] ) represents the basic input - output relations for the amplitude operators in the fourier space .
the coefficients therein are determined only by the ( complex ) permittivities and thicknesses of the layers of the plate and the permittivities of the surrounding media . in particular , for ( i.e. , ) eq .
( [ 35 ] ) immediately yields the input - output relations for the special case of the plate being surrounded by vacuum .
the input - output relations in the form of eq .
( [ 35 ] ) are generally valid , independent of the mechanism of creation of the incoming fields and the fields inside the layers of the plate .
it is worth noting that they take into account both propagating waves and evanescent waves .
in particular in the case when the plate is surrounded by vacuum , then the input and output amplitude operators are associated with propagating waves or evanescent waves in @xmath30-direction .    the temporal evolution ( in the heisenberg picture ) of the amplitude operators
is determined by the time dependence of the basic variables @xmath11 [ or equivalently , @xmath93 , which is governed by the hamiltonian ( [ 2 ] ) . in the special case
when the plate is free of active atomic sources , then the fields inside the layers of the plate represent langevin noise sources associated with material absorption . if absorption is disregarded and the plate is surrounded by vacuum that may contain active sources , then eq
( [ 35 ] ) solves a scattering problem of the type considered in ref .
@xcite , namely scattering of fields containing evanescent components .
in general , absorption can not be disregarded and both atomic sources and langevin noise sources contribute to a field ( e.g. , in cavity qed ) , which can contain both propagating and evanescent components ( e.g. , in near - field scanning probe microscopy ) .    the input - output relations ( [ 35 ] ) enable one to calculate correlation functions of the output field amplitudes in terms of those of the input field amplitudes and the amplitudes of the fields inside the plate .
the simplest case is the calculation of the expectation values of the field amplitudes .
for example , in a typical scattering arrangement the active light sources are located outside the plate , so that the field inside the plate is the absorption - assisted random field whose thermal - equilibrium expectation value vanishes .
application of eq .
( [ 35 ] ) thus leads to the expectation - value relations @xmath94 @xmath95 which exactly correspond to standard results in classical optics . on the other hand , when the active sources are located ( in a cavity - like system ) inside the plate , we find that @xmath96 , \end{split}\ ] ] @xmath97 . \qquad
\ ;     \end{split}\ ] ] needless to say that when there are no active sources inside the plate , then is valid and thus . clearly , higher - order correlation functions of the outgoing amplitude operators do not vanish in this case in general .
for example , the spectral intensity ( in the @xmath0-space ) of the radiation outgoing from a plane of a plate in thermal equilibrium at temperature @xmath98 , for chosen polarization , is proportional to @xmath99 , where @xmath100 & & = w_{q\,{\rm out}}^{(0,n)}({\bf k},\omega)\,\delta(\omega-\omega ' ) \delta({\bf k}-{\bf k}').\end{aligned}\ ] ] applying eq .
( [ 35 ] ) together with eqs .
( [ 38 ] ) and ( [ 17 ] ) and making use of @xmath101 & & = n(\omega , t)\ , \delta_{\mu\mu ' } \delta(\omega-\omega')\delta({\bf k}-{\bf k}')\end{aligned}\ ] ] ( @xmath34 @xmath102 @xmath103 ) , where @xmath104^{-1}\ ] ] ( @xmath105 , boltzmann constant ) is the well - known bose - einstein distribution function , we derive @xmath106 & & \left . + \ , \bigl|r_{j/0}^q\bigr|^2 c_{q--}^{(j)}({\bf k},\omega ) + \left[r_{j/0}^q c_{q-+}^{(j)}({\bf k},\omega ) + { \rm c.c.}\right]\right\}\qquad\end{aligned}\ ] ] ( and @xmath107 accordingly ) , with the coefficients @xmath108 being given in eqs .
( [ a15 ] ) and ( [ a16 ] ) in the appendix .      the input - output relations ( [ 35 ] ) in the @xmath110-space can be transformed , according to eq .
( [ 15 ] ) , into the @xmath109-space in a straightforward manner .
writing the electric - field operator at the boundaries of the plate as @xmath111        & = &        \hat{\bf{e}}_{q\,{\rm in } } ^{(0 ) } ( { \bm \rho   } , \omega )        +        \hat{\bf{e}}_{q\,{\rm out } } ^{(0 ) } ( { \bm \rho   } , \omega ) ,       \\[1ex ] \nonumber\\              \label{48 }        \hat{\bf e }    ^{(n)}({\bm \rho   } , \omega )        & = & \left.\hat{\bf e }   ^{(n)}(z , { \bm \rho   } , \omega ) \right|_{z=0^+ }        \nonumber\\*[1ex ]        & = & \
\hat{\bf{e}}_{q\,{\rm in}}^{(n)}({\bm \rho   } , \omega )        +       \hat{\bf{e}}_{q\,{\rm out}}^{(n ) } ( { \bm \rho } , \omega ) ,        \end{aligned}\ ] ] we derive    @xmath112&&\hspace{10ex } + \sum _ { j=1 } ^{n-1 } \int   { \mbox{\rm d}}^2 { \rho}'\ ,      \left\lbrace { \mbox{\boldmath $ \phi$ } } _ { 0 + } ^{(j ) }          ( { \bm \rho } , { \bm \rho ' } , \omega)\cdot     \hat{\bf e } _ { + }   ^{(j)}({\bm \rho ' } , \omega ) +   { \mbox{\boldmath $ \phi$ } } _ { 0- } ^{(j ) } ( { \bm \rho } , { \bm \rho ' } , \omega)\cdot     \hat{\bf e } _ { - }   ^{(j)}({\bm \rho ' } , \omega ) \right\rbrace , \\[1ex ] \label{50 } \lefteqn {     \hat{\bf e } _ { { \rm out } }   ^{(n)}({\bm \rho   } , \omega )     = \int { \mbox{\rm d}}^2 { \rho}'\ ,     { \mbox{\boldmath $ r$}}_{n/0 } ( { \bm \rho } , { \bm \rho ' } , \omega)\cdot     \hat{\bf e } _ { { \rm in } }   ^{(n)}({\bm \rho ' } , \omega )     + \int { \mbox{\rm d}}^2 { \rho}'\ ,     { \mbox{\boldmath $ t$}}_{0/n } ( { \bm \rho } , { \bm \rho ' } , \omega)\cdot     \hat{\bf e } _ { { \rm in } }   ^{(0)}({\bm \rho ' } , \omega ) } \nonumber\\[1ex]&&\hspace{10ex } + \sum _ { j=1 } ^{n-1 } \int   { \mbox{\rm d}}^2 { \rho}'\ ,      \left\lbrace { \mbox{\boldmath $ \phi$ } } _ { n+ } ^{(j ) }          ( { \bm \rho } , { \bm \rho ' } , \omega)\cdot     \hat{\bf e } _ { + }   ^{(j)}({\bm \rho ' } , \omega ) +   { \mbox{\boldmath $ \phi$ } } _ { n- } ^{(j ) } ( { \bm \rho } , { \bm \rho ' } , \omega)\cdot     \hat{\bf e } _ { - }   ^{(j)}({\bm \rho ' } , \omega ) \right\rbrace,\end{aligned}\ ] ]    where the electric - field operators @xmath113 and @xmath114 are respectively related to the amplitude operators @xmath115 and @xmath116 as @xmath117 \label{52 }     \hat{\bf e } _ { { \rm in } }   ^{(n)}({\bm \rho   } , \omega )     \!=\ ! \frac { 1 } { ( 2 \pi ) ^2}\!\int \!{\mbox{\rm d}}^2 { k}\sum _ { q = p , s}\ !     { \bf e}^{(n)}_{q- } ( { \bf k } )     \hat{\mathrm{e}}_{q\,{\rm in}}^{(n ) } ( { \bf k   } , \omega )     e^{i{\bf k } { \bm { \rho } } } \ ! , \quad \\[1ex ] \label{53 }     \hat{\bf e } _ { { \rm \pm } }   ^{(j)}({\bm \rho   } , \omega )     \!=\ ! \frac { 1 } { ( 2 \pi ) ^2}\!\int \!{\mbox{\rm d}}^2 { k}\sum _ { q = p , s}\ !     { \bf e}^{(j)}_{q\pm } ( { \bf k } )     \hat{\mathrm{e}}_{q\,{\rm \pm}}^{(j ) } ( { \bf k   } , \omega )     e^{i{\bf k } { \bm { \rho } } } \!. \quad \;\end{aligned}\ ] ] in the above , the tensor - valued reflection and transmission integral kernels , respectively , read as @xmath118 \label{55 } \lefteqn {        { \bf r}_{n/0 } ( { \bm \rho } , { \bm \rho ' } , \omega ) } \nonumber \\ & &        = \frac { 1 } { ( 2 \pi)^2}\int { \mbox{\rm d}}^2 { k}\sum _ { q = p , s }        { \bf{e}}^{(n)}_{q+ } ( { \bf k   } ) r^q _ { n/0}{\bf{e}}^{(n)}_{q- }        ( { \bf k } )        e^{i{\bf k } ( { \bm \rho}-{\bm \rho'})}\qquad       \end{aligned}\ ] ] and @xmath119 \label{57 } \lefteqn {        { \bf t}_{n/0 } ( { \bm \rho } , { \bm \rho ' } , \omega ) } \nonumber \\ & &        = \frac { 1 } { ( 2 \pi)^2}\int { \mbox{\rm d}}^2 { k } \sum _ { q = p , s }        { \bf{e}}^{(0)}_{q- } ( { \bf k   } ) t^q _ { n/0}{\bf{e}}^{(n)}_{q-}({\bf k } )        e^{i{\bf k } ( { \bm \rho}-{\bm \rho'})},\qquad\end{aligned}\ ] ] and the integral kernels that relate the outgoing fields to the intraplate fields are given by @xmath120 \lefteqn { \label{59 }        { \bf \phi}_{n\pm } ( { \bm \rho } , { \bm \rho ' } , \omega )   } \nonumber \\ & &        = \frac { 1 } { ( 2 \pi)^2}\int { \mbox{\rm d}}^2 { k } \sum _ { q = p , s }        { \bf e}^{(j)}_{q+ } ( { \bf k   } ) \phi^{(j ) } _ { q\ , n \pm }        { \bf e}^{(j)}_{q\pm}({\bf k } )        e^{i{\bf k } ( { \bm \rho}-{\bm \rho'})}.\qquad\end{aligned}\ ] ] note that the integral kernels depend only on the difference , because of the translational symmetry in the @xmath121- and @xmath122-directions . the input - output relations in the form of eqs .  ( [ 49 ] ) , ( [ 50 ] ) can be considered as a special realization of the general relations given for a dielectric scattering region of arbitrary shape in ref .
@xcite .        in many applications it may be advantageously to express the incoming and the outgoing field operators in terms of appropriately chosen bosonic operators .
recalling the basic commutation relations ( [ 4 ] ) and ( [ 5 ] ) and the orthogonality of polarization unit vectors as given by eqs .
( [ 27 ] ) and ( [ 28 ] ) , from eqs .
( [ 31 ] ) and ( [ 32 ] ) together with eq .
( [ 17 ] ) and ( [ 18 ] ) and eqs .
( [ 36 ] ) and ( [ 37 ] ) we find that the input amplitude operators satisfy the commutation relations @xmath123 } \nonumber \\[1ex ] & &        \,=\ , c_{q\,{\rm in}}^{(0,n)}({\bf k},\omega)\ ,
\delta_{qq'}\delta(\omega-\omega')\delta({\bf k}-{\bf k}'),\end{aligned}\ ] ] where @xmath124        \label{62 }        c_{q\,{\rm in}}^{(n ) }         = \frac{\pi \hbar}{\epsilon_0 }       \left(\frac{\omega}{c}\right)^2        \frac{\beta_{n}'}{|\beta_{n}|^2 }        \,{\bf e}^{(n)}_{q- } ( { \bf k   } ) \cdot{\bf e}^{(n)*}_{q- }        ( { \bf   k } ) .        \end{aligned}\ ] ] similarly , the output amplitude operators and the intraplate amplitude operators can be shown to satisfy commutation relations of the same type , i.e. , @xmath125 } \nonumber \\[1ex ] & &        = \ , c_{q\,{\rm out}}^{(0,n)}({\bf k},\omega)\ ,        \delta_{qq'}\delta(\omega-\omega')\delta({\bf k}-{\bf k } ' ) , \\[1ex ]        \label{64 } \lefteqn {        \left[\hat{\mathrm{e}}_{q\lambda}^{(j ) }        ( { \bf k   } , \omega ) ,        \hat{\mathrm{e}}_{q'\lambda'}^{(j)\dagger }        ( { \bf k   } ' , \omega')\right ] } \nonumber \\[1ex ] & &         = \,c_{q\lambda\lambda'}^{(j)}({\bf k},\omega)\ ,         \delta_{qq'}\delta(\omega-\omega')\delta({\bf k}-{\bf k}')\end{aligned}\ ] ] [ for the coefficients @xmath126 and @xmath108 , respectively , see eqs .
( [ a12 ] ) , ( [ a13 ] ) and eqs .
( [ a15 ] ) , ( [ a16 ] ) in the appendix ] . needless to say that input amplitude operators , that refer to different sides of the plate , commute , input amplitude operators and intraplate amplitude operators commute , and intraplate amplitude operators that refer to different layers also commute .
note that output amplitude operators that refer to different sides of the plate do not commute in general [ see eqs .
( [ a9 ] ) and ( [ a10 ] ) in the appendix ] , which is similar as in the one - dimensional situation ( cf .
@xcite ) .    from eqs .
( [ 60 ] ) and ( [ 63 ] ) it is seen that bosonic input and output operators can be introduced according to @xmath127 thus @xmath128 } \nonumber \\[1ex]&&\hspace{6ex }        = \
,         \delta_{qq'}\delta(\omega-\omega')\delta({\bf k}-{\bf k } ' ) .
\end{aligned}\ ] ] since for in eq .
( [ 64 ] ) is valid in general , it is useful to introduce intraplate bosonic operators according to the superposition @xmath129 and choosing the coefficients @xmath130 [ eqs .
( [ a20 ] ) and ( [ a21 ] ) in the appendix ] in such a way that @xmath131 } \nonumber\\&&\hspace{6ex }        = \,\delta_{\lambda\lambda'}\delta_{qq ' }        \delta ( \omega -\omega ' ) \delta({\bf k } - { \bf k ' } ) .
\end{aligned}\ ] ]    substituting eqs .
( [ 65 ] ) and ( [ 67 ] ) into eq .
( [ 35 ] ) , we may express the input - output relations in the @xmath0-space in terms of bosonic operators ,    @xmath132        \hat{a}_{q\,{\rm out}}^{(n ) } ( { \bf k   } , \omega )        \end{array } \right )        = \left (        \begin{array}{cc }        \tilde{r}^q_{0/n}({\bf k},\omega )        & \tilde{t}_{n/0 } ^q({\bf k},\omega)\\[1ex ]        \tilde{t}_{0/n } ^q({\bf k},\omega )        &   \tilde{r } ^q _ { n/0}({\bf k},\omega )        \end{array } \right )        \left (        \begin{array}{c }        \hat{a}_{q\,{\rm in}}^{(0 ) } ( { \bf k   } ,   \omega ) \\[1ex ]        \hat{a}_{q\,{\rm in}}^{(n ) } ( { \bf k   } ,   \omega )         \end{array } \right )        + \sum _ { j = 1 } ^{n-1 } \left (        \begin{array}{cc }        \tilde{\phi } _ { q\ , 0 + } ^{(j)}({\bf k},\omega )        & \tilde{\phi } _ { q\ , 0- } ^{(j)}({\bf k},\omega ) \\[1ex ]        \tilde{\phi } _ { q\ , n+ } ^{(j)}({\bf k},\omega )        & \tilde{\phi } _ { q\ , n- } ^{(j)}({\bf k},\omega )        \end{array } \right )        \left (        \begin{array}{c }        \hat { a}_{q+}^{(j ) } ( { \bf k   } , \omega)\\[1ex ]        \hat { a}_{q-}^{(j ) } ( { \bf k   } , \omega )        \end{array } \right)\!,\quad        \end{aligned}\ ] ]    where the coefficients [ modified in comparison with eq .  ( [ 35 ] ) ] read    @xmath133        \label{71 }        \tilde{t } ^q _ { 0/n } =           \sqrt { \frac{c_{q\,{\rm in}}^{(0 ) } }        { c_{q\,{\rm out}}^{(n ) } } }        \,t ^q _ { n/0 } \ ,
,        \qquad        \tilde{r}_{n/0}^q =         \sqrt { \frac{c_{q\,{\rm in}}^{(n ) } }        { c_{q\,{\rm out}}^{(n ) } } }        \,r_{n/0 } ^q \,,\end{aligned}\ ] ]    and @xmath134 note that for ( i.e. , perpendicular incidence of light on the plate ) , the input and output relations ( [ 69 ] ) exactly agree with previous results @xcite obtained within a one - dimensional treatment .
let us turn to the limiting case that the space outside the plate  except for possible active atomic sources  may be regarded as being vacuum , i.e. , @xmath135 @xmath136 \sqrt{\omega^2/c^2 - k^2 } & \omega / c > k. \end{array } \right.\ ] ]    for the propagating - field components observed for , the coefficients @xmath137 and @xmath138 , respectively , read @xmath139 and @xmath140 as it can be seen from eqs .
( [ 61 ] ) and ( [ 62 ] ) and eqs .
( [ a23 ] ) and ( [ a24 ] ) in the appendix .
hence , bosonic input and output operators can be defined according to eq .
( [ 65 ] ) . moreover , the output amplitude operators that refer to different sides of the plate commute [ eq .
( [ a22 ] ) in the appendix ] , which implies that the associated bosonic operators commute as well . making use of eqs .
( [ 75 ] ) and ( [ 76 ] ) , we see that the modified reflection and transmission coefficients ( [ 70 ] ) and ( [ 71 ] ) , respectively , become identical with the original ones . for , eq .
( [ 69 ] ) thus reduces to    @xmath141        \hat{a}_{q\,{\rm out}}^{(n ) } ( { \bf k   } , \omega )        \end{array } \right )        = \biggl (        \begin{array}{cc }        r^q_{0/n}({\bf k   } , \omega ) & t_{n/0 } ^q({\bf k   } , \omega)\\[1ex ]        t_{0/n } ^q({\bf k   } , \omega ) &   r ^q _ { n/0}({\bf k   } , \omega )        \end{array } \biggr )        \left (        \begin{array}{c }        \hat{a}_{q\,{\rm in}}^{(0 ) } ( { \bf k   } ,   \omega ) \\[1ex ]        \hat{a}_{q\,{\rm in}}^{(n ) } ( { \bf k   } ,   \omega )         \end{array } \right )        + \sum _ { j = 1 } ^{n-1 } \left (        \begin{array}{cc }        \tilde{\phi } _ { q\ , 0 + } ^{(j)}({\bf k   } , \omega )        & \tilde{\phi } _ { q\ , 0- } ^{(j)}({\bf k   } , \omega ) \\[1ex ]        \tilde{\phi } _
{ q\ , n+ } ^{(j)}({\bf k   } , \omega )        & \tilde{\phi } _ { q\ , n- } ^{(j)}({\bf k   } , \omega )        \end{array } \right )        \left (        \begin{array}{c }        \hat { a}_{q+}^{(j ) } ( { \bf k   } , \omega)\\[1ex ]        \hat { a}_{q-}^{(j ) } ( { \bf k   } , \omega )        \end{array } \right),\quad        \end{aligned}\ ] ]    where the transformation matrix connecting the bosonic output operators with the bosonic input operators is exactly the same as that for the corresponding amplitude operators in eq .  ( [ 35 ] ) , and the matrix equation    @xmath142        t_{0/n } ^q &   r ^q _ { n/0 }        \end{array } \right )        \left (        \begin{array}{cc }        r^{q\ast}_{0/n } & t_{0/n } ^{q\ast}\\[1ex ]        t_{n/0 } ^{q\ast } &   r ^{q\ast } _ { n/0 }        \end{array } \right )        + \sum _ { j = 1 } ^{n-1 } \left (        \begin{array}{cc }        \tilde{\phi } _ { q\ , 0 + } ^{(j ) }        & \tilde{\phi } _ { q\ , 0- } ^{(j ) } \\[1ex ]        \tilde{\phi } _ { q\ , n+ } ^{(j ) }        & \tilde{\phi } _ { q\ , n- } ^{(j ) }        \end{array } \right )        \left (        \begin{array}{cc }        \tilde{\phi } _ { q\ , 0 + } ^{(j)\ast }        & \tilde{\phi } _ { q\ , n+ } ^{(j)\ast } \\[1ex ]        \tilde{\phi } _ { q\ , 0- } ^{(j)\ast }        & \tilde{\phi } _ { q\ , n- } ^{(j)\ast }        \end{array } \right )        = \left (        \begin{array}{cc }        1 & 0 \\[1ex ]        0 & 1         \end{array } \right )        \end{aligned}\ ] ]    is valid . in this case , the field operators @xmath143 and @xmath144 at the boundaries of the multilayer plate [ see eqs .
( [ 29 ] ) , ( [ 30 ] ) , ( [ 36 ] ) , and ( [ 37 ] ) ] can be represented as    @xmath145 ,   \label{79}\\[1ex ]   { \bf \hat{e } }    ^{(n)}(0^+,{\bf k   } , \omega ) =   -\frac { \omega } { c } \sqrt{\frac { \pi \hbar    } { \beta _ n\epsilon _ 0 } }   \sum _ { q = p , s}\ ,   \left [ { \bf{e}}_{q-}^{(n ) } ( { \bf k } ) \hat{a}_{q\,{\rm in}}^{(n ) }   ( { \bf k   } , \omega )   +   { \bf{e}}_{q+}^{(n ) } ( { \bf k   } ) \hat{a}_{q\,{\rm out}}^{(n ) }   ( { \bf k   } , \omega)\right ]    \label{80}\end{aligned}\ ] ]    ( @xmath146 @xmath60 @xmath147 @xmath60 @xmath148 , @xmath149 @xmath150 @xmath151 ) .
note that in eqs .
( [ 79 ] ) and ( [ 80 ] ) is nothing said about the location of the active light sources .
clearly , if there are no active light sources outside the plate , then @xmath143 and @xmath144 evolve effectively freely there [ cf .
( [ 33 ] ) and ( [ 34 ] ) ] .    for the evanescent - field components observed for , the coefficients @xmath152 [ eq .  (
[ 61 ] ) ] and @xmath153 [ eq .  (
[ 62 ] ) ] identically vanish , because of @xmath154 @xmath60 @xmath155 @xmath60 @xmath156 .
recalling eq .
( [ 65 ] ) , we see that bosonic input operators can not be introduced for the evanescent - field components . hence , it is impossible to extend the validity of eq .
( [ 77 ] ) to the evanescent - field components . to treat them
, one has therefore to go back to the generally valid input - output relations ( [ 35 ] ) for the amplitude operators .
clearly , when there are no active light sources at any finite distance from the plate , then the input field can be regarded as being effectively a propagating field [ see eqs .
( [ 31 ] ) and ( [ 32 ] ) ] , that is to say , a field that does not contain evanescent input components
. however , the output field may contain evanescent output components resulting from the intraplate field .
applying the green - tensor formalism of the quantization of the electromagnetic field in the presence of inhomogeneous , causal , linear dielectric bodies , we have derived the most general three - dimensional quantum input - output relations for the electromagnetic - field operators at an arbitrary multilayer dielectric plate , and we have given the commutation relations needed .
taking into account ( i ) dispersion and absorption of the plate and the surrounding medium , ( ii ) possible active light sources inside and/or outside the plate , and ( iii ) both propagating - field components and evanescent - field components , the input - output relations connect the output field operators at the boundaries of the plate with the input field operators at the boundaries of the plate and the fields inside the layers of the plate .
the outgoing fields outside the plate can then be obtained from the ones at the boundaries of the plate via quantum langevin equations , with the langevin noise sources being determined by the active and passive light sources in the respective half - space outside the plate .
we have given the input - output relations in terms of both algebraic relations for amplitude operators in the @xmath0-space and integral relations for field operators in the coordinate space .
further , we have shown that in the @xmath0-space it is possible , provided that the incoming fields contain effectively propagating components , to rewrite the input - output relations in order to obtain them in terms of relations between bosonic operators instead of amplitude operators .
therefore the input - output relations for the amplitude operators are more general than those for the bosonic operators .
finally , we have studied the limiting case of the plate being surrounded by vacuum . in this case
, the input - output relations in terms of bosonic operators only apply to the propagating - field components , for which is valid .
an application to the evanescent - field components would run into contradictions . to treat them
, one has therefore to go back to the amplitude operators .
the derived input - output relations generalize previous results obtained from a one - dimensional calculation @xcite ) . in the special case of a plate that is embedded in vacuum , without active light sources at any finite distance from the plate ( and without active light inside the plate ) , they reduce to those given in ref .
@xcite , if evanescent - field components are disregarded .    in conclusion ,
the derived input - output relations are suited for studying the quantum statistical properties of electromagnetic fields in the presence of planar multilayered structures ( such as cavity - like systems or photonic crystals ) , with special emphasis on absorption - assisted quantum decoherence , including evanescent - field effects .
we would like to thank ho trung dung and c. raabe for useful discussions .
let us begin with the calculation of the commutators of the output amplitude operators @xmath157 . from eq .
( [ 35 ] ) it follows that @xmath158&&\hspace{6ex }        + \,t_{n/0}^q   \hat{\mathrm{e}}_{q\,{\rm in}}^{(n ) }        ( { \bf k   } ,   \omega ) +        \hat{\mathrm{f}}_{q}^{(0)}({\bf k   } ,   \omega )        \end{aligned}\ ] ] and @xmath159&&\hspace{6ex }        + \ ,        r_{n/0}^q   \hat{\mathrm{e}}_{q\,{\rm in}}^{(n ) }        ( { \bf k   } ,   \omega ) +        \hat{\mathrm{f}}_{q}^{(n)}({\bf k   } ,   \omega ) ,        \end{aligned}\ ] ] where @xmath160 \\      & = i \omega \mu _ 0 { \bf e}^{(0)}_{q-}({\bf
k } )          \cdot\sum _ { j=1}^{n-1 } \int _ { [ j ] } { \mbox{\rm d}}z\ ,          { \mbox{\boldmath $ g$}}^{(0j)}(0 , z , { \bf k   } , \omega)\cdot          \hat { \bf j}^{(j ) } ( z , { \bf k   } ,
\omega ) , \end{split}\ ] ] @xmath161\\      & = i \omega \mu _ 0 \,{\bf e}^{(n)}_{q+}({\bf k } )      \cdot\sum _ { j=1}^{n-1 } \int _ { [ j ] } { \mbox{\rm d}}z\ ,          { \mbox{\boldmath $ g$}}^{(nj)}(0 , z , { \bf k   } , \omega)\cdot          \hat { \bf j}^{(j ) } ( z , { \bf k   } , \omega ) . \end{split}\ ] ] hence , the ( relevant ) commutators of the output amplitude operators that refer to different sides of the plate can be given by    @xmath162 } \nonumber\\[1ex ] & &         = \,r^q_{0/n } t_{0/n } ^{q\ , * }          \left [         \hat{\mathrm{e}}_{q\,{\rm in } } ^{(0 ) } ( { \bf k   } , \omega ) ,         \hat{\mathrm{e}}_{q'\,{\rm in } } ^{(0)\dagger } ( { \bf k   } ' , \omega ' )         \right ]          + t_{n/0 } ^q r^{q\,*}_{n/0 }          \left [         \hat{\mathrm{e}}_{q\,{\rm in } } ^{(n ) } ( { \bf k   } , \omega ) ,         \hat{\mathrm{e}}_{q'\,{\rm in } } ^{(n)\dagger } ( { \bf k   } ' , \omega ' )         \right ]          +         \left [         \hat{\mathrm{f}}_{q } ^{(0 ) } ( { \bf k   } , \omega ) ,         \hat{\mathrm{f}}_{q ' } ^{(n)\dagger } ( { \bf k   } ' , \omega ' )         \right ]   \nonumber\\[1ex ] & &         = \,\delta_{qq'}\ , \delta(\omega-\omega')\,\delta({\bf k}-{\bf k}')\ ,         \left ( c_{q\,{\rm in}}^{(0 ) } r^q_{0/n } t_{0/n } ^{q\ , * }       +      c_{q\,{\rm in}}^{(n ) }      t_{n/0 } ^q r^{q\,*}_{n/0 } \right )      +      \left [          \hat{\mathrm{f}}_{q } ^{(0 ) } ( { \bf k   } , \omega ) ,          \hat{\mathrm{f}}_{q ' } ^{(n)\dagger } ( { \bf k   } ' , \omega ' )          \right ] ,          \end{aligned}\ ] ]    where eqs .  ( [ 60 ] )  ( [ 62 ] ) have been used . making use of eqs .
( [ 17 ] ) and ( [ 18 ] ) and recalling the basic commutation relations ( [ 4 ] ) and ( [ 5 ] ) , we derive @xmath163 } \nonumber\\*[1ex ] & &       = \,\delta_{qq'}\,\delta(\omega-\omega')\ ,           \delta({\bf k}-{\bf   k}')\ ,           \frac{4\pi \hbar}{\epsilon _ 0}\left(\frac { \omega}{c}\right ) ^2           \ , e^{(0)}_{q-\ , \mu } ( { \bf k   } ) \ ,           e^{(n)\ , * } _ { q+\ , \mu ' } ( { \bf k   } )           \sum_{j=1}^{n-1}\,\int_{[j ] } { \mbox{\rm d}}z\ ,           g^{(0j)}_{\mu\nu}(0 , z , { \bf k } , \omega ) \ ,           \frac { \omega ^2 } { c^2 } \varepsilon_{j}''\ ,           g^{(nj)*}_{\mu'\nu}(0 , z , { \bf k } , \omega ) \nonumber\\[1ex ] & &       = \ , \delta_{qq'}\,\delta(\omega-\omega')\ ,           \delta({\bf k}-{\bf   k}')\ ,           \frac{4\pi \hbar}{\epsilon _ 0}\left(\frac { \omega}{c}\right ) ^2           \ , e^{(0)}_{q-\ , \mu } ( { \bf k   } ) \ ,             e^{(n)\ , * } _ { q+\ , \mu ' } ( { \bf k   } ) \ ,       \biggl\lbrace \sum_{j=0}^{n}\,\int_{[j ] } { \mbox{\rm d}}z\ ,           g^{(0j)}_{\mu\nu}(0 , z , { \bf k } , \omega ) \ ,           \frac { \omega ^2 } { c^2 } \varepsilon_{j}''\ ,           g^{(nj)*}_{\mu'\nu}(0 , z , { \bf k } ,
\omega )        \biggr .
\nonumber\\[1ex]&&\qquad           \biggl
.       - \int_{-\infty}^{0 } { \mbox{\rm d}}z\ ,           g^{(00)}_{\mu\nu}(0 , z , { \bf k } , \omega ) \ ,           \frac { \omega ^2 } { c^2 } \varepsilon_{0}''\ ,           g^{(n0)*}_{\mu'\nu}(0 , z , { \bf k } , \omega )       - \int^{\infty}_{0 } { \mbox{\rm d}}z\ ,           g^{(0n)}_{\mu\nu}(0 , z , { \bf k } , \omega ) \ ,           \frac { \omega ^2 } { c^2 } \varepsilon_{n}''\ ,           g^{(nn)*}_{\mu'\nu}(0 , z , { \bf k } , \omega )           \biggr\rbrace.\end{aligned}\ ] ]    it is not difficult to calculate the last two integrals in eq .
( [ a6 ] ) , by using the explicit expression ( [ 20 ] ) for @xmath164 . to calculate the sum of integrals @xmath165 }
{ \mbox{\rm d}}z\,\ldots$ ] , we employ the integral relation ( [ 11 ] ) for the classical green tensor @xmath166 , rewritten in terms of @xmath167 ,    @xmath168 } { \mbox{\rm d}}z''\ ,        g^{(jj'')}_{\mu\nu}(z , z '' , { \bf k } , \omega ) \ ,        \frac { \omega ^2 } { c^2 } \varepsilon_{j''}''\ ,        g^{(j'j'')*}_{\mu'\nu}(z ' , z '' , { \bf k } , \omega )        =   \frac{1}{2i }        g^{(jj')}_{\mu\mu'}(z ,
z ' , { \bf k } , \omega ) } \nonumber\\[1ex ] & &        -\,\frac{1}{2i }        g^{(j'j)*}_{\mu'\mu}(z ' , z , { \bf k } , \omega )        + g^{(jj')}_{\mu\nu}(z , z ' , { \bf k } , \omega ) e_{z\,\nu }        \,\frac{\varepsilon_{j'}''}{\varepsilon_{j'}^*}\,e_{z\ , \mu ' }        + e_{z\,\mu } \ ,
\frac{\varepsilon_{j}''}{\varepsilon_{j}}\ ,        g^{(j'j)*}_{\mu'\nu}(z ' , z , { \bf k } , \omega ) e_{z\,\nu }        \ , .
\end{aligned}\ ] ]    after lengthy , but straightforward calculations we then derive , on using eqs .
( [ 61 ] ) and ( [ 62 ] ) , @xmath169 } \nonumber\\[1ex ] & &        = \left [   \hat{\mathrm{e}}_{p\,{\rm out } } ^{(0 ) } ( { \bf k   } ,
\omega ) ,     \hat{\mathrm{e}}_{s\,{\rm out } } ^{(n)\dagger } ( { \bf k   } ' , \omega ' )     \right]=0,\end{aligned}\ ] ] @xmath170     = \delta(\omega-\omega')\,\delta({\bf k}-{\bf k } ' ) } \qquad \nonumber\\[1ex]&&\times\ ;     \frac{\pi \hbar}{\epsilon _ 0}\left(\frac { \omega}{c}\right)^2     \left ( \frac{i \beta_0 ' ' } { { |\beta _ 0 |}^2 } t^{s\,*}_{0/n}-     \frac { i \beta _ n''}{{|\beta _ n |}^2 } t_{n/0}^{s}\right ) ,     \qquad\end{aligned}\ ] ]    @xmath171=     \frac{\pi \hbar}{\epsilon _ 0}\left(\frac { \omega}{c}\right ) ^2     \left\ {     \frac { t_{n/0 } ^p}{\beta _ n { |k_n| } ^2 }     \left(k^2\frac{k_n ^2}{k_n^{*\,2 } } - { |\beta _ n |}^2 \right )      \right . }   \nonumber\\[1ex ] & &     \left .
+ \,\frac { t_{0/n } ^{p\ , * } } { \beta _ 0 ^ * { |k_0| } ^2 }     \left(k^2\frac{k_0^{*\,2}}{k_0 ^2 } - { |\beta _ 0 |}^2 \right )      - \frac{\beta_n'}{{|\beta_n|}^2}\,t^{p}_{n/0 }     \left[{\bf e}^{(n)}_{p- } ( { \bf k   } ) \cdot{\bf e}^{(n)*}_{p- } ( { \bf k   } ) \right ]     \left[{\bf e}^{(n)*}_{p- } ( { \bf k})\cdot\,{\bf e}^{(n)*}_{p+ } ( { \bf k   } ) \right ]     \right . \nonumber\\[1ex ] & &     -\ , \frac{\beta_0'}{{|\beta_0|}^2}\,t^{q*}_{0/n}\ ,     \left[{\bf e}^{(0)*}_{q+ } ( { \bf k   } ) \cdot{\bf e}^{(0)}_{q+ } ( { \bf k   } ) \right ]     \left[{\bf e}^{(0)}_{q+ } ( { \bf k   } ) \cdot{\bf e}^{(0)}_{q- } ( { \bf k   } ) \right ]     \biggr\}\ ,     \delta(\omega-\omega')\,\delta({\bf k}-{\bf k}')\ , .   \end{aligned}\ ] ]    next , let us consider the commutators of the output amplitude operators that refer to the same sides of the plate . performing the same steps as before , we now arrive at @xmath172 } \nonumber\\[1ex ] & &     = \left [   \hat{\mathrm{e}}_{p\,{\rm out } } ^{(0 ) } ( { \bf k   } , \omega ) ,     \hat{\mathrm{e}}_{s\,{\rm out } } ^{(0)\dagger } ( { \bf k   } ' , \omega ' )     \right ] = 0,\end{aligned}\ ] ] @xmath173     = \delta(\omega-\omega')\,\delta({\bf k}-{\bf k } ' ) } \qquad \nonumber\\[1ex]&&\times\ ;     \frac{\pi \hbar}{\epsilon _ 0}\left(\frac { \omega}{c}\right ) ^2     \left ( \frac{\beta_0 ' } { { |\beta _ 0 |}^2 } +     \frac { 2 \beta _ 0''}{{|\beta _ 0 |}^2 } r_{0/n}^{s\,\prime \prime}\right ) ,     \qquad\end{aligned}\ ] ]    @xmath174= \frac{\pi \hbar}{\epsilon _ 0}\left(\frac{\omega}{c}\right ) ^2     \left\ {     \frac { r_{0/n } ^p}{\beta _ 0 { |k_0| } ^2 }     \left(k^2\frac{k_0 ^2}{k_0^{*\,2 } } - { |\beta _ 0 |}^2 \right )     \right .
} \nonumber\\[1ex ] & &      \left .
+ \,\frac { r_{0/n } ^{p\ , * } } { \beta _ 0 ^ * { |k_0| } ^2 }     \left(k^2\frac{k_0^{*\,2}}{k_0 ^2 } - { |\beta _ 0 |}^2 \right )     - \frac { k ^2 } { { |k_0| } ^2}\left(\frac{\beta _ 0}{k_0^{*\,2}}+     \frac{\beta _ 0 ^*}{k_0^{2 } }   \right )     + \frac { 2 \beta _ 0 ' } { { |k_0| } ^2}\left ( \frac { k ^4 }     { { |\beta_0| } ^2 { |k_0| } ^2 } + 1 \right )      \right .
\nonumber\\[1ex ] & &      + \,\frac{\beta_0 ' } { { |\beta _ 0 |}^2 }     \left [ { \bf e}^{(0)}_{p+ } ( { \bf k   } )     \cdot{\bf e}^{(0)*}_{p+ } ( { \bf k   } ) \right ]     \left [ \left| r_{0/n}^q\right|^2     - \left|{\bf e}^{(0)}_{p- } ( { \bf k } ) \cdot     { \bf e}^{(0)}_{p+ } ( { \bf k   } ) + r_{0/n } ^p\right| ^2 \right ]      \biggr\}\ ,     \delta(\omega-\omega')\,\delta({\bf k}-{\bf k}').\end{aligned}\ ] ]    the commutators @xmath175 $ ] are obtained from eqs .
( [ a11 ] )  ( [ a13 ] ) , by making the replacements @xmath146 @xmath176 @xmath147 , @xmath177 @xmath176 @xmath178 , @xmath179 @xmath176 @xmath180 , and @xmath181 @xmath176 @xmath182 .    finally , it can easily be proved that the intraplate amplitude operators ( [ 38 ] ) satisfy the commutation relations (; ) @xmath183 } \nonumber\\[1ex ] & &     = c_{q\lambda\lambda'}^{(j)}({\bf k},\omega)\ ,     \delta_{qq'}\,\delta(\omega-\omega')\,\delta({\bf k}-{\bf k}'),\end{aligned}\ ] ] where @xmath184 & & \times\ ;        \left ( e^{\pm 2\beta _ j '' d_j } - 1 \right )        { \bf e}^{(j)}_{q \pm } ( { \bf k   } )        \cdot{\bf e}^{(j)*}_{q\pm } ( { \bf k } ) , \\[1ex ]        \label{a16 } \lefteqn {       c_{q\,{\rm \pm \mp } } ^{(j ) } = \pm i\frac{\pi \hbar}{\epsilon_0 }       \left(\frac{\omega}{c}\right)^2       \frac{\beta_{j}''}{|\beta_{j}|^2 } } \nonumber\\[1ex ] & & \times\ ;     \left ( e^{\mp 2i \beta _ j ' d_j } - 1 \right )     { \bf e}^{(j)}_{q \pm } ( { \bf k   } )     \cdot{\bf e}^{(j)*}_{q\mp } ( { \bf    k}).\end{aligned}\ ] ] using eqs .  ( [ a14 ] )  ( [ a16 ] ) , we find that the operators @xmath185&&\times\ ;       \left [ e^{i \beta_{j } d_{j } }     \hat{\mathrm{e}}_{q+}^{(j ) } ( { \bf k   } , \omega )   \pm     \hat{\mathrm{e}}_{q-}^{(j ) } ( { \bf k   } , \omega ) \right ] ,      \label{a17}\end{aligned}\ ] ] where @xmath186&&\,\times\ ;     \left\ {     \beta _ j ' \sinh(\beta_j ' ' d_j)\ ,     { \bf{e}}^{(j)}_{q+}({\bf k } )     \cdot { \bf{e}}^{(j)\,*}_{q+}({\bf k } )     \right . \nonumber\\[1ex]&&\hspace{2ex }      \left .
\pm\ , \beta _ j ''     \sin(\beta_j ' d_j)\ ,     { \bf{e}}^{(j)}_{q+}({\bf k } )     \cdot { \bf{e}}^{(j)\,*}_{q-}({\bf k } )     \right\ } ^{1/2}\!\qquad\end{aligned}\ ] ] satisfy bosonic commutation relations . from eq .
( [ a17 ] ) it then follows that @xmath187 where @xmath188 \label{a21 } & & \tau_{q - \pm}^{(j ) } = \pm { \textstyle\frac { 1 } { 2 } } \xi^{(j)}_{q \pm}({\bf k   } , \omega).\end{aligned}\ ] ]    in the limiting case ( when the plate is surrounded by vacuum ) one has to distinguish between ( propagating - field components ) and ( evanescent - field components ) . in the first case we have , so that eqs .
( [ a8 ] )  ( [ a10 ] ) reduce to @xmath189        = 0,\ ] ] and eqs .
( [ a11 ] )  ( [ a13 ] ) simplify to @xmath190 } \nonumber\\[1ex ] & &        = \,\frac{\pi \hbar}{\epsilon_0}\left(\frac{\omega}{c}\right)^2        \frac{1}{\beta_0}\,\delta_{qq'}\,\delta(\omega-\omega')\ ,        \delta({\bf k}-{\bf k}').\qquad        \end{aligned}\ ] ] @xmath191 } \nonumber\\[1ex ] & &        = \,\frac{\pi \hbar}{\epsilon_0}\left(\frac{\omega}{c}\right)^2        \frac{1}{\beta_n}\,\delta_{qq'}\ ,        \delta(\omega-\omega')\,\delta({\bf k}-{\bf   k}').\qquad        \end{aligned}\ ] ] in the second case we have . equations ( [ a8 ] )  ( [ a10 ] ) then lead to @xmath192 } \nonumber\\[1ex ] & &        = \,\frac{\pi \hbar}{\epsilon_0 }        \left(\frac{\omega}{c}\right)^2        \frac{2 t^{q\ , \prime \prime}_{0/n}}{{|\beta_0| } }        \,\delta_{qq'}\ , \delta(\omega-\omega')\,\delta({\bf k}-{\bf   k } ' ) ,        \qquad\end{aligned}\ ] ] and from eqs .
( [ a11 ] )  ( [ a13 ] ) we find @xmath193 } \nonumber\\[1ex ] & &        = \,\frac{\pi \hbar}{\epsilon_0 }        \left(\frac{\omega}{c}\right)^2        \frac{2 r_{0/n}^{q\,\prime \prime}}{{|\beta_0|}}\ ,        \delta_{qq'}\,\delta(\omega-\omega')\,\delta({\bf k}-{\bf   k } ' ) ,        \qquad\end{aligned}\ ] ] and similarly @xmath194 } \nonumber\\[1ex ] & &        = \,\frac{\pi \hbar}{\epsilon_0 }        \left(\frac{\omega}{c}\right)^2        \frac{2 r_{n/0}^{q\,\prime \prime}}{{|\beta_n|}}\ ,        \delta_{qq'}\,\delta(\omega-\omega')\,\delta({\bf k}-{\bf   k } ' ) .
\qquad\end{aligned}\ ] ]    note that the commutation relations ( 42 ) in ref .
@xcite for the bosonic output operators are only valid for propagating waves , not for evanescent ones .
it should also be noted that the free - space light modes @xmath195 used in ref .
@xcite are introduced by explicitly requiring that , where @xmath196 corresponds to @xmath0 in the present paper .
hence , these modes solely represent propagating waves .