arxiv:quant-ph/ v1 11 Nov 1995

Transcription

1 Quantization of electromagnetic fields in a circular cylindrical cavity arxiv:quant-ph/950v Nov 995 K. Kakazu Department of Physics, University of the Ryukyus, Okinawa 903-0, Japan Y. S. Kim Department of Physics, University of Maryland, College Park, Maryland 074 We present a quantization procedure for the electromagnetic field in a circular cylindrical cavity with perfectly conducting walls, which is based on the decomposition of the field. A new decomposition procedure is proposed; all vector mode functions satisfying the boundary conditions are obtained with the help of this decomposition. After expanding the quantized field in terms of the vector mode functions, it is possible to derive the Hamiltonian for this quantized system. PACS:4.50.-p, k

2 Introduction The behavior of an atom depends on the mode structure of the electromagnetic field surrounding it []. The mode structure in a cavity is determined by its boundary conditions which are quite different from the periodic boundary conditions used in Lorentz-covariant electrodynamics. There are many interesting effects on the interaction of the atom with the field with various boundary conditions. It is possible to calculate them and compare them with experimental data []. Thanks to many new high-tech experimental instruments, cavity electrodynamics has become one of the most active branches of physics in recent years. In quantum electrodynamics, we use the plane waves satisfying the periodic boundary condition whose quantization volume eventually becomes infinite. As a result, the orientation of the field is not essential. On the other hand, the field in a cavity must satisfy the boundary conditions and the orientation of the field with respect to cavity walls becomes an important factor. For this reason, we have to use vector mode functions to deal with this problem. These functions satisfy the boundary conditions and the boundary effect in turn is contained in the functions. Thus, for a given cavity, the most important step in the quantization process is to construct the vector mode functions. Indeed, several authors have investigated the quantization procedure for the field in cavities [3]. They have studied the field quantization, the emission rate, and the atomic level shifts in the cavities with one or two infinite plane mirrors. Cavity electrodynamics with rectangular coordinate system does not present further significant mathematical problems. However, the case is quite different for cavities with curved boundary conditions. Even for circular cylindrical or spherical cavities, the quantization procedure has not yet been firmly established. The purpose of the present paper is to discuss in detail the mathematical problems in constructing the vector mode functions with a circular cylindrical boundary condition. We hope to discuss the spherical cavity in a later paper. Recently, we have carried out the field quantization for several different rectangular boundary conditions including a rectangular tube using the vector mode functions [4]. These functions have been derived with the help of an orthogonal matrix. However, the procedure developed there is not applicable to other cavities in a straightforward manner. In this paper, we present a more general method of obtaining the vector mode functions and apply it to the field quantization in a circular cylindrical cavity with perfectly conducting walls. This will be summarized in two theorems. With this in mind, we first decompose the field into three parts in the circular cylindrical coordinates in the following section. Then, in Sec. 3, we apply it to the circular cylindrical cavity. After obtaining all vector mode functions by using the decomposition, we arrive at the quantized field and the quantized electromagnetic Hamiltonian in Sec. 4. Decomposition of Electromagnetic fields Let us now derive the decomposition formula for the electromagnetic field from Maxwell s equations. Using the cylindrical coordinate system, we shall decompose the solutions of Maxwell s equations into three terms, in preparation for the field quantization within the

3 circular cylindrical cavity in Sec. 3. We emphasize here again that this decomposition process is the key to the quantization. Maxwell s equations for the electric field E and the magnetic field B in free space are given by E = 0, (.) B = 0, (.) E + t B = 0, (.3) B c te = 0, (.4) where c is the velocity of light in free space and t = / t. As is well known, it follows from Maxwell s equations that the fields E and B satisfy the wave equations: ( ) ( c t E = 0, ) c t B = 0, (.5) where is the Laplacian operator. The electromagnetic field E and B can be written in the cylindrical coordinates (r, ϕ, z) as E = E T + E z, B = B T + B z, (.6) where E T = e r E r + e ϕ E ϕ is the transverse component of the field and E z = e z E z. Here e r, e ϕ, and e z are the unit vectors in the r, ϕ, and z directions, respectively. Similarly, E, E, and φ (φ a scaler function) can also be divided into two parts. For example, E = T E + z E, where the first term contains the derivatives with respect to r and ϕ, while the second term contains the derivative with respect to z. For simplicity, the derivatives with respect to r, ϕ, and z are described as / r = r, / ϕ = ϕ and / z = z, respectively. Since z E z = 0, Eq. (.3) gives Similarly, from Eq. (.4), we have t B z = T E T, t B T = T E z z E T. (.7) c te z = T B T, Equations (.7) and (.8) give which leads to c te T = T B z + z B T. (.8) c t E T = T t B z z ( T E z + z E T ), c t B T = c T t E z z ( T B z + z B T ), (.9) c t E T + z z E T T T E z = t B z E z, c t B T + z z B T T T B z = c te z B z,(.0) 3

4 where we have used T E z = E z. Let us note that z z E T = z E T, T T E z = T E z, c t E T = ( E) T = E T, (.) where ( E) T is the transverse part of E; the last equation holds in the present coordinates, although it is not correct in general. (For example, it is not correct in the spherical coordinates.) Then, from Eq. (.0), we arrive at T E = E z + t B z, T B = c te z B z. (.) To rewrite Eq. (.), we must obtain an expression for the components E z and B z. Suppose that the field is in a finite region and expand E z and B z in terms of a certain complete system of functions with mode s: E z (r, t) = s B z (r, t) = s (E zs (r, t) + c.c.), (B zs (r, t) + c.c.), (.3) where E zs (r, t) = Ẽzs(r)e iω st, B zs (r, t) = B zs (r)e iω st. (.4) Here ω sσ (ω sσ 0, σ =, ) is determined by using given boundary conditions. Since E z and B z satisfy the wave equation, their components satisfy the Helmholtz equations: where We assume that the components satisfy E zs = k s E zs, B zs = k s B zs, (.5) k sσ = ω sσ c. (.6) z E zs = h s E zs, z B zs = h s B zs, (.7) where h sσ is also determined by the boundary conditions. Then we have two dimensional Helmholtz equations: where T E zs = g se zs, T B zs = g sb zs, (.8) g sσ = k sσ h sσ. (.9) Now we define two functions F σ from E zs and B zs with gsσ 0 as F σ (r, t) = [F sσ (r, t) + c.c.] = [ F sσ (r) e iωsσt + c.c.], (.0) g sσ 0 4 g sσ 0

5 where F s = E zs /g s, F s = B zs /g s. (.) The functions F σ and their components F sσ satisfy the same equations as the z components of the field: F σ = c t F σ, T F sσ = g sσf sσ (.) From Eqs. (.) and (.) the component F sσ is a solution of the Poisson equation and then the functions T F σ satisfy T F s = E zs, T F s = B zs (.3) T F = (gs F s + c.c.) = (E zs + c.c.), g s 0 g s 0 T F = g s 0 (g sf s + c.c.) = g s 0 (B zs + c.c.). (.4) On the other hand, if there is a component E zs or B zs with gsσ = 0, Eq. (.8) reduces to the two dimensional Laplace equation: T E zs = 0 or T B zs = 0. Now define E 0z and B 0z as E 0z = (E zs + c.c.), B 0z = (B zs + c.c.), (.5) gs =0 gs =0 which satisfy T E 0z = 0, T B 0z = 0. (.6) Then we have a useful expression for E z and B z ; for any E z and B z there exist functions F σ given by Eq. (.0) and functions E 0z and B 0z are given by Eq. (.5) such that E z = T F + E 0z, B z = T F + B 0z. (.7) It is worth emphasizing again that the functions F σ are constructed with the components E zs and B zs with g sσ 0, while E 0z and B 0z consist of the components with g sσ = 0. Using Eq. (.7) and defining F σ by we can rewrite Eq. (.) as F σ = e z F σ, (.8) T ( E F + t F ) = E0z + t B 0z, (.9) T ( B c tf F ) = c te 0z B 0z, (.30) where E 0z = e z E 0z and B 0z = e z B 0z. Let us define E 0 and B 0 as the quantities in the parentheses at the left hand side in Eqs. (.9) and (.30), respectively. The results of this section are then summarized in the following theorem. 5

6 Theorem : The field can be decomposed into three components as follows: E = F t F + E 0, B = c tf + F + B 0, (.3) where E 0 and B 0 satisfy T E 0 = E 0z + t B 0z, T B 0 = c te 0z B 0z. (.3) Theorem plays a central role in performing field quantization in this paper. Take the z component of Eq. (.3), we find that the z components of E 0 and B 0 are given by E 0z and B 0z in Eq. (.7), respectively. Furthermore, Eq. (.3) is consistent with Eq. (.). Let us consider the physical meaning of the above decomposition formula. First it is worth emphasizing that each term in Eq. (.3) is a solution to Maxwell s equations. That is, each term containing F σ (F σ term) satisfies Eqs. (.) - (.4), which is easily shown by using the fact that F σ satisfies the wave equation (.). As a result, the third term E 0,B 0 is also a solution to Maxwell s equations, so that it satisfies Eq. (.), as mentioned above. Next take a particular example of the boundary condition satisfying E 0z, B 0z = 0. Then the field E 0,B 0 describes so-called the TEM (transverse electromagnetic) mode [5]. Since the F term has no z component of the magnetic field, it becomes TM (transverse magnetic), while the F term becomes TE (transverse electric), because it does not contain the z component of the electric field. In the above particular case where E 0z, B 0z = 0, Maxwell s equations for the field E 0,B 0 reduce to T E 0 = 0, T E 0 = 0, T B 0 = 0, T B 0 = 0, z E 0 + t B 0 = 0, z B 0 c te 0 = 0. (.33) 3 Determination of Functions F σ We are considering here the cavity enclosed by a circular cylindrical wall with radius a and height L: r < a, 0 < z < L. We assume that the cavity has perfectly conducting walls at z = 0, L and at r = a. The tangential component of the electric field E tan and the normal component of the magnetic field B norm must accordingly vanish at the boundary of the cavity. The above boundary condition reduces to that for the z components E z = 0, r B z = 0, (r = a), (3.) B z = 0, z E z = 0, (z = 0, L). (3.) 6

7 It is easy to get the second condition in Eq. (3.) if we take the ϕ component of Eq. (.4). Equation (.) leads to the second condition in Eq. (3.). Before obtaining the components E zs and B zs, we show that the folowing lemma. Lemma : E 0 = 0, B 0 = 0. Proof: First we show that E 0z, B 0z = 0. Consider the component E 0z, which satisfies the Laplace equation T E 0z = 0 [see Eq. (.6)] with the boundary condition E 0z = 0 at the boundary (r = a). Then we have E 0z = 0; this is a well known property of the Laplace equation. Next consider the component B 0z satisfying T B 0z = 0 with the boundary condition r B 0z = 0 (r = a). Then B 0z must be a constant, i.e., independent of the variables r and ϕ: B 0z = f(z, t). Applying the two-dimensional Gauss theorem to B 0 = T B 0T + z B 0z = 0, we get z B 0z = 0, because B 0r = 0 (r = a) and, as a result, B 0z is independent of the variables r, ϕ, and z. On the other hand, B 0z must be zero at z = 0, L, which leads to B 0z = 0. As mentioned in the preceding section, the field E 0,B 0 becomes TEM satisfying Eq. (.33), because E 0z, B 0z = 0. There exists no TEM in the present cavity. That is, from Eq. (.33) we have T B 0 = 0 and T B 0 = 0. As a result, there exists a function ϕ such that B 0 = T ϕ, which gives B 0 = 0 by the boundary condition at r = a. Similarly, we get a vector function A such that E 0 = A. The boundary condition leads us to E 0 = 0. Q.E.D. Since E 0,B 0 = 0, Theorem shows that the field is constructed from the functions F σ. Consequently, from Eq. (.0), the field only contains the components with gsσ 0. Moreover, it is easy to prove that gsσ > 0 by using Gauss theorem and the boundary condition. Next we solve the Helmholtz equation (.5) for the components E zs and B zs under the boundary conditions (3.) and (3.). Since gsσ > 0, the solution is given by E zs (r, t) = C s (t)j m (χ s r/a) e imϕ cos(nπz/l), B zs (r, t) = C s (t)j m (χ s r/a) e imϕ sin(nπz/l), (3.3) where the mode index is s = (m, µ, n) (m = 0, ±, ±, ; µ =,, 3, ; n = 0,,, ), J m is the Bessel function of the first kind, χ s χ mµ is the µth zero point of J m, and χ s χ mµ is the µth zero point of J m, the derivative of J m. It follows from Eq. (.4) that C sσ (t) exp( iω sσ t). From the solution (3.3) we have which gives T E zs = (k s (nπ/l) )E zs = (χ s /a) E zs, T B zs = (k s (nπ/l) )B zs = (χ s /a) B zs, (3.4) g sσ = (χ sσ/a), k sσ = (χ sσ/a) + (nπ/l). (3.5) Let us next obtain the functions F σ, which are defined in Eqs. (.0) and (.); they are given by F s (r, t) = a E zs (r, t) hωs i a χ s (t)ψ s (r), s ε 0 7

8 F s (r, t) = a B zs (r, t) hωs i a χ s (t)ψ s (r), (3.6) s ε 0 where we have introduced a sσ and ψ sσ : a sσ (t) = a sσ (0)e iωsσt, (3.7) ψ s (r, t) = c s J m (χ s r/a) e imϕ cos(nπz/l), (3.8) ψ s (r, t) = c s J m (χ s r/a) e imϕ sin(nπz/l), (3.9) where c sσ are normalization constants. The functions ψ sσ have the orthonormality property, which is used in quantization in the next section. Lemma : c where c dr = cavity rdrdϕdz, V is the cavity volume, and dr ψ sσ (r)ψ s σ(r) = c sσ V α sσ δ ss, (3.0) α s α mµ = J m+(χ mµ ), α s α mµ = J m(χ mµ ) J m+(χ mµ ). (3.) Here the quantity cos(nπz/l) in Eq. (3.8) must be changed to / when n = 0. Proof: The orthonormality property (3.0) is easily derived from the following equality for the Bessel functions: a dr rj m (χ mµσ r/a)j m (χ mµ σr/a) = a α sσ δ µµ. (3.) 0 To prove Eq. (3.), let us first observe the following two integrals [6] [ dxxj m (χ i x)j m (χ j x) = χi J 0 χ i χ m+ (χ i )J m (χ j ) χ j J m+ (χ j )J m (χ i ) ] j [ = χi J χ i χ m(χ i )J m (χ j ) + χ j J m (χ i )J m(χ j ) ] (3.3) j 0 dxxj m (χ ix) = [ J m (χ i ) J m (χ i )J m+ (χ i ) ] = [ J m (χ i ) J m+(χ i ) J m(χ i )J m+ (χ i ) ], (3.4) where χ i (i =,, ) is a positive real number and χ i χ j. Considering that χ mµ is the µth zero point of J m (x) and that J m (χ mµ )+J m+ (χ mµ ) = 0, from Eqs. (3.3) and (3.4), we have dxxj m (χ mµ x)j m (χ mµ x) = α mµ δ µµ, (3.5) 0 where α mµ is given by Eq. (3.). If we change the variable x to r/a, Eq. (3.5) gives Eq. (3.) for σ =. Similarly, Eqs. (3.3) and (3.4) lead to dxxj m (χ mµ x)j m (χ mµ x) = α mµ δ µµ, (3.6) 0 where α mµ is given by Eq. (3.). Equation (3.6) gives us Eq. (3.) for σ =. Q.E.D. 8

9 4 Vector Mode Functions and Field Quantization After defining the vector mode functions, we quantize the field and obtain the quantized Hamiltonian. The mode functions are constructed according to the decomposition formula given in Eq. (.3) In the cavity with the perfectly conducting walls, the decomposition (.3) in Theorem reduces to E = F t F, B = c tf + F, (4.) because of E 0,B 0 = 0 [see Lemma ]. Substituting the functions F σ in Eq. (3.6) into E in Eq. (4.), we find E(r, t) = i sσ hωsσ ε 0 [ asσ (t)u sσ (r) a sσ (t)u sσ (r)], (4.) where the vector mode functions u sσ are given by u s = e z ψ s, u s = iω s e z ψ s. (4.3) Similarly, substituting F σ into B and considering the equality u s = k s e zψ s, (4.4) we obtain B(r, t) = h [ asσ (t) u sσ (r) + a sσ ε 0 ω sσ(t) u sσ(r) ]. (4.5) sσ Before obtaining the quantized field and Hamiltonian, let us present some properties of the mode functions. Lemma 3: u sσ = 0, (4.6) ( + k sσ )u sσ = 0, (4.7) u sσ tan = 0, u sσ norm = 0, (on walls), (4.8) dru sσ(r) u s σ (r) = δ ss δ σσ, (4.9) c Proof: We give a proof of the orthonormality property (4.9), because it is obvious to prove the other equations. From Eq. (4.3), the mode function u s is rewritten as which gives u s = ( e z ψ s ) e z ψ s = k se z ψ s + ( z ψ s ), (4.0) u s u s = [ k s k s k s h s ] ψ s ψ s +k s z[ ( z ψ s )ψ s ] + [ ( z ψ s ) zψ s ]. (4.) 9

10 Taking into account Gauss theorem and the boundary conditions, we find that the last two terms in Eq. (4.) have no effect on the volume integration of u s u s. Equation (4.) gives the orthonormality relation for u s : dru s(r) u s (r) = δ ss, (4.) c where we have used Eq. (3.0) and set c a c s =. (4.3) V α s χ sωs In the case of u s, we make use of the equality ( e z ψs ) ( e zψ s ) = e z ψs ( e zψ s ) + [ e z ψs ( e zψ s ) ] = χ s a ψ s ψ s + (ψs Tψ s ). (4.4) The orthonormality property of u s is then obtained from Eq. (3.0) if we set the normalization constant: a c s =. (4.5) V α s χ sωs Finally we show that the mode functions u s and u s are orthogonal to each other. Since u s has no z component, we see that u s u s = iω s ( k s e z ψ s + zψ s) ( ez ψ s ) = iω s ( z ψ s ) ( e zψ s ) = iω s [( z ψ s)( e z ψ s ) ], (4.6) which has no effect on the volume integration of u s u s. Q.E.D. We are now ready to carry out the field quantization. To get the quantized field the functions a sσ (t) are regarded as annihilation operators, which annihilate photons with indices sσ and satisfy the commutation relation: [a sσ (t), a s σ (t)] = δ ss δ σσ, (4.7) where a s σ (t) are creation operators. Then the field given by Eqs. (4.) and (4.5) becomes operators. The quantized electromagnetic Hamiltonian H R is obtained by using the equality dr ( u sσ (r)) ( u s σ (r)) = k s σ dru sσ (r) u s σ (r). (4.8) c These are summarized in the following theorem. Theorem : The quantized field and the Hamiltonian are given by E(r, t) = i hωsσ [ asσ (t)u sσ (r) a sσ ε sσ(t)u sσ(r) ], (4.9) 0 B(r, t) = h [ asσ (t) u sσ (r) + a sσ ε 0 ω sσ(t) u sσ(r) ]. (4.0) sσ H R = hω ( sσ a sσ a sσ + a sσ a ( sσ) = hω sσ a sσ a sσ + ). (4.) sσ sσ 0 c

11 Theorem will play a fundamental role when we investigate the electromagnetic property of a physical system in the circular cylindrical cavity. There are also important issues in the present cavity including, for example, spontaneous (stimulated) emission, absorption, the atomic energy shifts, and the Casimir effect; they will be investigated with the help of Theorem. It is necessary to expand the field in terms of the vector mode functions as in Theorem. According to Fourier analysis, the field can be expanded in terms of, for instance, the periodic functions as in ordinary quantum electrodynamics. However, in this case, the operators a sσ and a sσ cannot describe the physical photons. Consider the situation where the field consists of only one mode corresponding to the photon with a sσ u sσ, where u sσ satisfies the periodic boundary condition. This must be impossible, because this field does not satisfy the boundary conditions; the photon is then fictitious. 5 Conclusions The quantization procedure of the field in the cavity without E 0,B 0 is as follows: (a) obtain the decomposition formula (.3) in the circular cylindrical coordinates; (b) solve the Helmholtz equations (.5) for the components E zs and B zs under the boundary conditions (3.) and (3.); (c) determine the functions F σ, which are the solutions of the Poisson equation (.3), from E zs, B zs ; (d) substitute F σ into the decomposition formula (.3) and obtain the vector mode functions satisfying the orthonormality property (4.9); (e) then we arrive at the quantized field and Hamiltonian in Theorem. As has been shown explicitly in Theorem in Sec. 4, the energy hω sσ of the photon depends on the polarization index as well as the wave number index. This is one of the characteristics of the circular cylindrical cavity. The photon energy is independent of the polarization index in a rectangular cavity and in free space. In the entire process of quantization, the decomposition formula in Theorem plays the essential role. We must find such a formula in an appropriate coordinate system if we study other types of cavities. Also, if a cavity allows the third term E 0,B 0 in (.3), we have to construct a methodology to deal with the problem. We hope to discuss the quantization of electromagnetic fields inside a spherical cavity in a future publication. Acknowledgments We would like to thank Prof. Y. Ezawa, Dr. R. Ray, Dr. G. Gat, Dr. S. Kobayashi, Dr. Kurennoy, and Mr. K. Ohshiro for numerous valuable discussions. We are also grateful to Prof. M. S. Kim for bringing some of the references to our attention. References [] E. M. Purcell, Phys. Rev. 69, 68 (946); D. Kleppner, Phys. Rev. Lett. 47, 33 (98).

12 [] G. Barton and N. S. J. Fawcett, Phys. Reports 70, (988); P. Meystre, ibid 9, 43 (99); H. Walther, ibid 9, 63 (99); Cavity Quantum Electrodynamics, edited by P. R. Berman (Academic Press, New York, 994); P. W. Milonni, The Quantum Vacuum: An Introduction to Quantum Electrodynamics (Academic Press, New York, 994). [3] C. K. Carniglia and L. Mandel, Phys. Rev. D 3, 80 (97); P. W. Milonni, Phys. Rev. A 5, 35 (98); E. A. Power and T. Thirunamachandran, Phys. Rev. A 5, 473 (98). [4] K. Kakazu and Y. S. Kim, Phys. Rev. A 50, 830 (994). [5] J. D. Jackson, Classical Electrodynamics (John Wiley and Sons, New York, 975). [6] G. N. Watson, A Treatice on the Theory of Bessel Functions (Cambridge University Press, Cambridge, 96).