Effective Hamiltonians in Nonrelativistic Quantum Electrodynamics

In this paper, we consider some second-order effective Hamiltonians describing the interaction of the quantum electromagnetic field with atoms or molecules in the nonrelativistic limit. Our procedure is valid only for off-energy-shell processes, specifically virtual processes such as those relevant for ground-state energy shifts and dispersion van der Waals and Casimir-Polder interactions, while on-energy-shell processes are excluded. These effective Hamiltonians allow for a considerable simplification of the calculation of radiative energy shifts, dispersion, and Casimir-Polder interactions, including in the presence of boundary conditions. They can also provide clear physical insights into the processes involved. We clarify that the form of the effective Hamiltonian depends on the field states considered, and consequently different expressions can be obtained, each of them with a well-defined range of validity and possible applications. We also apply our results to some specific cases, mainly the Lamb shift, the Casimir-Polder atom-surface interaction, and the dispersion interactions between atoms, molecules, or, in general, polarizable bodies.


I. INTRODUCTION
In molecular quantum electrodynamics, that is, the quantum theory of atoms and molecules interacting with the electromagnetic field in the nonrelativistic limit, several processes of great interest are of a high order in atomfield coupling [1][2][3]. For example, interatomic dispersion interactions such as van der Waals and Casimir-Polder interactions between two atoms or molecules are fourth-order processes [4][5][6], and, in the case of three or more atoms, many-body effects start from the sixth order [1,7,8]. Additionally, the resonance energy transfer between molecules may involve high-order perturbative calculations in the atom-field coupling [9]. In such cases, the number of relevant Feynman diagrams rapidly grows with the perturbative order, with consequent increasing complexity of the calculations. For this reason, the possibility of finding approximated effective Hamiltonians allowing for the simplification of calculations is very important, and, hopefully, this will also yield a transparent interpretation and physical insights into the relevant physical processes involved [6,[10][11][12][13][14]. Effective Hamiltonians can be also used in dynamical (time-dependent) nonequilibrium situations [15]. All these possibilities have fostered the investigation of effective Hamiltonians containing an interaction term that is at least quadratic in the atom-field coupling, where the response of the atom is included in quantities such as, for example, its polarizability, thus allowing for a considerable reduction of the perturbative order required for calculating specific processes and of the number of relevant Feynman diagrams [6]. Very recently, resummation techniques, in which the polarizability is summed to any order, have been developed, and they could also be of great importance for the evaluation of radiative processes such as the Lamb shift and van der Waals interactions for nanostructured materials [16][17][18].
In this paper, we obtain and review different forms of effective Hamiltonians, at the second order in the atom-field coupling, used in nonrelativistic quantum electrodynamics, stressing the range and limit of validity of each of them, as well as their physical interpretation. Application to retarded and nonretarded dispersion interactions and to the atom-surface Casimir-Polder interaction is also discussed. A physical interpretation in terms of the response of the atoms or molecules to vacuum field fluctuations or to real photons is also outlined, as well as the possible role of dissipation in the response function. More specifically, we here present in a clear and organic way the different forms of second-order effective Hamiltonians that can be obtained, depending on the specific field states considered. We also show and stress that the response of the atom is different in the case of vacuum fluctuations or when photons are present; in fact, we find that in the first case, the atom responds through a real function of the frequency, while in the second case through its dynamical polarizability, that has poles in the real frequency axis. This paper is organized as follows. In Section II, the general expression of the effective Hamiltonian and of its matrix elements between atom+field states is obtained. In Section III, specific forms of the effective Hamiltonian in relevant cases are given, and their limits of validity are stressed, as well as the response function of the atoms in the various cases considered, and the possible role of dissipation. Application to interatomic dispersion interactions and to the Casimir-Polder atom-surface interaction, both in the nonretarded and retarded cases, is outlined. Finally, Section IV is devoted to our conclusive remarks.

II. THE GENERAL EXPRESSION OF THE EFFECTIVE HAMILTONIAN IN MOLECULAR QUANTUM ELECTRODYNAMICS
We first start from the Hamiltonian of one atom, placed at R, interacting with the quantum electromagnetic field, in the multipolar coupling scheme and within dipole approximation [5]. Extension to the case of two or more atoms is straightforward, at the order considered. The Hamiltonian of the system is where states |φ ℓ are a complete set of atomic states with energies E ℓ , and a † kλ and a kλ are, respectively, creation and annihilation operators relative to the field mode (kλ), with λ = 1, 2 being the polarization index and satisfying the usual bosonic commutation relations, and ω k = ck; µ = eq is the electric dipole moment operator of the atom, with q being the electron coordinate. For simplicity, we assume only one atomic electron is taking part in the radiative process considered, and q ℓm = φ ℓ |q|φ m is its matrix element between atomic states. Although this simplification, which is strictly valid for hydrogen-like systems, is not essential in the present calculation, it allows one to include one-electron dipole moment matrix elements only and simplify the discussion of the results; it can, however, be simply overtaken by introducing a full dipole moment operator.
Finally, E(r) is the electric field operator where E(kλ; r) is a Fourier component of the electric field, with f kλ (r) the field mode functions taking into account the boundary conditions present. In the free space, we have whereê kλ (λ = 1, 2) are polarization unit vectors, such thatê kλ ·ê kλ ′ = δ λλ ′ ,ê kλ ·k = 0, and V is the quantization volume. The multipolar-coupling Hamiltonian (1) contains also a term equal to 2π d 3 r(P ⊥ (r)) 2 , where P ⊥ (r) is the transverse part of the polarization field P(r) = i erδ(r − r i ), with r i being the position of an atomic electron [1,19,20]. This is a second-order term in the electric charge that, although it is important for the Lamb shift, does not contribute to the dispersion interactions [1,6,20,21], and for this reason we do not include it, except whenever necessary.
We now follow and briefly review the standard general procedure used in nonrelativistic quantum electrodynamics to obtain a second-order effective Hamiltonian, which consists of applying the following transformation to the Hamiltonian (1) where z is an Hermitian operator (with some limitations on the state space where it is defined, as specified later on), chosen in such a way to eliminate, in the transformed Hamiltonian, the first-order terms in the electron charge e [10,13,14]. Up to the second order in e, we have The operator z is chosen in such a way to make vanishing the first-order terms in the transformed Hamiltonian (5), that is so that the transformed Hamiltonian becomes with the effective Hamiltonian, expressed in terms of the Hermitian operator z, given bỹ This Hamiltonian is at the second order in the atom-field coupling. A second-order energy shift of the system due to the atom-field coupling is then given by the average value ofH ef f on the state of the system at hand. Fourth-order corrections, as in the evaluation of interatomic dispersion interactions, can be obtained by a second-order approach if the effective Hamiltonian is used, as explicitly discussed in the next section.
Equation (6) is an implicit relation defining the operator z of the transformation. From this relation, we can obtain its matrix elements between atom+field states of the generic form |{p}, φ P , {p}, and φ P , denoting, respectively, generic field and atom states ({p} denotes a general set of number states of the field for all the modes allowed by the boundary conditions), with energy E p = i ω i p i and E P , respectively. We obtain Assuming E n − E m + E N − E M = 0, that is, that the matrix element is taken between unperturbed states with a different energy, we get the matrix elements of the operator z Equation (10) defines all matrix elements of z, except those on the energy shell, that are undefined; thus, the transformation (4) with the matrix elements (10) is not defined on the energy shell. Therefore, the effective Hamiltonian so obtained is not valid when energy-conserving processes are relevant, unless some regularization procedure of the vanishing energy denominators, taking into account dissipation, is used [22][23][24]. This regularization can be done by adding an imaginary part to the energy denominators shifting the poles to the lower complex frequency half-plane, in agreement with the causality requirement; in our case, this can be done only phenomenologically because we are working within a second-order perturbative approach for a closed system. Because we will obtain effective Hamiltonians mainly to calculate atomic energy shifts and (many-body) dispersion interactions for ground-state systems, where only virtual processes are involved, or when resonant processes are suppressed (in a cavity, for example), in general, this is not a strong limitation for our purposes. However, we wish to point out that our results cannot be directly extended to the case of excited atomic states when resonant processes are present.
For the reasons mentioned above, substitution of the matrix elements (10) of z into (5) allows us to obtain only the off-energy-shell matrix elements of the effective Hamiltonian and not its full expression as an operator in the Hilbert space of the system. For this reason, we can obtain different expressions in terms of field operators, according to the type of field states at hand. The matrix elements of the effective Hamiltonian are where the index s indicates photonic modes with frequency ω s = ck s , and the sum over s runs over all field modes. Using the mode expansion (2) for the electric field operator, the matrix elements appearing in the RHS of (11) can be cast in the form From Equation (12), it is easy to see that, in order that the quantity above be nonvanishing, it is necessary that the total number of photons in {m} and {n} (that is, the photon number summed over all field modes) must be equal or differ by two. There are thus three possibilities: (i) {n} = {m}, that is elements diagonal in the photon space state; (ii) the difference in {n} and {m} is all in a single mode (kλ), and thus mkλ = nkλ ± 2, with all other modes containing the same number of photons; (iii) the difference is by one photon in each of the two modes (kλ) and (kλ), with (kλ) = nkλ ± 1 and (kλ) = nkλ ± 1 (upper or lower sign for both modes), while all other modes in {m} and {n} have the same photon number.
In the next section, we will explicitly find the form of the effective Hamiltonian for some specific cases, which is relevant in molecular quantum electrodynamics, in particular for atom-surface and (many-body) atom-atom interactions for ground-state atoms or molecules. On the basis of the results obtained in this section, we will now find the relative effective Hamiltonians in the various cases, according to the relevant photon states involved, and point out the range and limit of application of the specific forms obtained.

III. DIAGONAL AND OFF-DIAGONAL MATRIX ELEMENTS OF THE EFFECTIVE HAMILTONIAN
We now evaluate the matrix elements of the effective Hamiltonian introduced in the previous section. As mentioned in the previous section, the expression of the effective Hamiltonian operator that we obtain, which acts only in the field space, and of its matrix elements, can differ according to the subspace of photon states considered. We work in the Schrödinger representation, and the energy shifts are to be evaluated by time-independent perturbation theory. Related effective Hamiltonians for evaluating energy shifts and dispersion interactions have been also used with a different approach, specifically in the Heisenberg representation in terms of time-dependent field operators, and by separating free and source (scattering) fields [25][26][27][28][29].
We now separately consider different cases; we obtain the relative explicit expressions of the effective Hamiltonian operator (off the energy shell), and apply them to some relevant physical systems, mainly van der Waals and Casimir-Polder dispersion interactions. The possible role of dissipation will be also discussed. The results we obtain are relative to a microscopic description of the system; however, we argue that some aspects of our results could be highly relevant also for macroscopic polarizable bodies and Casimir interactions between macroscopic dielectric bodies in the vacuum.

A. Diagonal Elements of the Effective Hamiltonian and Ground-State Systems (Zero Photons)
Let us first consider the unperturbed ground state of the system, that is, |{0 kλ }, g , where |{0 kλ } denotes the photon vacuum state and |g is the atom's ground state. Due to the atom-field coupling, this state is not an eigenstate of the interacting Hamiltonian, and, as it is well known, the interaction leads to a second-order energy shift that, after mass renormalization, yields the nonrelativistic Lamb shift of the ground-state atom [5,19,30,31]. This energy shift is a second-order effect in the electron charge.
Using our effective Hamiltonian (11), this energy shift can be obtained just through a simple first-order calculation, namely, as the expectation value of the effective Hamiltonian on the unperturbed state [14]. In fact, from (11), and using (2), we immediately obtain (Einstein's convention of summation over the repeated indices i, j is assumed) where |1 kλ are one-photon states, |φ L are intermediate atomic states, and E Lg = E L − E g ; also, E Lg > 0, and therefore the quantity in the round brackets in (13) has no poles. It involves off-energy-shell processes only. Thus, in (13) there is no need to introduce dissipation, and, very importantly, the real quantity appears, not the atomic dynamic polarizability (see also the following discussion). All relevant properties of the atom are embedded in this function. Equation (13) is the average value on the vacuum state of the electromagnetic field {0 kλ } of the following field operator that we can take as our effective Hamiltonian within the zero-photon subspace, as specified by the apex 0. By taking the average value of (15) on the field vacuum state |{0 kλ } , we immediately obtain the second-order energy shift due to the atom-field interaction where f kλ (R) are the appropriate mode functions for the system at hand, as introduced in Equation (2), evaluated at the position R of the atom. For example, in the case of an atom in the unbounded space, Equation (16) yields its ground-state Lamb shift [14] (as previously mentioned, the second-order term proportional to P 2 ⊥ (r) should be added in this case), while for an atom in front of a reflecting plate it yields the atom-surface Casimir-Polder interaction energy at zero temperature [32].
The effective Hamiltonian (15) has a clear physical interpretation. Vacuum electric field fluctuations with wavevector k induce a dipole moment on the atom given by and, in turn, this fluctuating dipole moment interacts with the vacuum field fluctuations yielding (we assume that a vacuum expectation value is taken) finally yielding (15) and (16). In this case, it is as if a Fourier component of the induced dipole moment with a given k interacts only with the Fourier component of the electric field having the same k. As we will show in the next subsection, this does not occur when off-diagonal elements of the effective Hamiltonian are relevant, as in the case of two-and many-body dispersion interactions between atoms or molecules. A very important point is that the response of the atom to the vacuum field fluctuations is not through its dynamical polarizability but through the function β g ij (ω k ), as defined by (14), which is a real quantity at any frequency without poles and related dissipative properties. This does not contrast with the fluctuation-dissipation theorem and the Kramers-Krönig dispersion relations, which indeed consider the response of a system to an external applied field (linear response theory), while in our case only vacuum field fluctuations act, which cannot induce real transitions. This result could be important also for macroscopic polarizable bodies and their role as boundary conditions in the Casimir effect, as well as for the long-lasting dispute in the literature about the most appropriate dielectric model (plasma or Drude model, for example) to be used in the Casimir effect for dielectrics [33][34][35]. We will address this point in a forthcoming publication.
From Equations (11) and (12), we can also obtain diagonal matrix elements of the effective Hamiltonian in the case of one mode populated with n kλ photons, with all other field modes being in their vacuum state, obtaining where is the ground-state atomic dynamic polarizability, with poles at ω k = E LN and obeying the dispersion relations, and is a function analogous to (14) for the generic atomic state N . Equation (21) has no poles if N is the ground state, as already discussed, while it can have poles in the case of other atomic states. We wish to point out and stress that Equation (19) clearly shows that the atom responds to real photons through its dynamical polarizability α ij (ω k ) and to the vacuum fluctuating field through the β N ij (ω k ) function. For atoms or molecules with a random orientation, being µ gL i µ Lg j =| µ gL | 2 δ ij /3, Equation (15) becomes with If the relevant frequencies of the field are such that ω k ≪ E Lg / , then β g av (ω k ) ≃ β g av (0) = α g , where α g is the isotropic static polarizability of the ground-state atom. Then, Equation (22) reduces to H (0) ef f ≃ − 1 2 αE 2 (R); we wish to stress that only in this limiting case, ω k → 0, our β function coincides with the (static) polarizability of the atom. This form of the effective Hamiltonian has been used, for example, for evaluating the far-zone (retarded) dispersion interaction of an atom placed near a conducting wall [27].

B. Off-Diagonal Elements of the Effective Hamiltonian
We now consider the off-diagonal elements of the effective Hamiltonian between atom-field states, which we can obtain from (11) and (12); as mentioned at the end of the previous section, in this case the total number of photons in the two states must differ by two. We find where {p} indicates all modes different from kλ and k ′ λ ′ , kλ = k ′ λ ′ . We can specialize Equation (24) to the case n kλ = m k ′ λ ′ = 0, {p} = {0} (vacuum state for all other field modes), φ M = φ N , relevant for the calculation of two-or three-body dispersion interactions between atoms or molecules in the vacuum space at zero temperature, allowing one to reduce, respectively, a fourth or sixth-order calculation to a second-or third-order calculation [10], where the symmetry between primed and non primed (kλ)s has been exploited, and the atomic matrix elements µ LN can be taken real without loss of generality. This quantity is the matrix element on the field states considered, |{0}, 0 kλ , 0 k ′ λ ′ , φ N and |{0}, 1 kλ , 1 k ′ λ ′ , φ N of the following operator that we can take as the effective Hamiltonian in the field subspace spanned by the states considered. The properties of the atoms or molecules relevant for the processes involved are contained in their dynamical polarizability, α ij (ω k ), and the effective Hamiltonian acts on the photon subspace only. The effective Hamiltonian (26) is different from (15), and its physical interpretation is different. The response of the atom in the present case is not through the β ij (ω k ) function but through its dynamical polarizability, and the interaction of the k Fourier component of the induced dipole moment is with the total electric field and not with only its k component, as in the previous diagonal case (18).
If the frequency of the (virtual) photons involved is much smaller than the relevant atomic transition frequencies, we can approximate the dynamical polarizabiblity with the static one, α ij (ω k ) ≃ α ij (0) = α ij , and (26) reduce it to which can be used to obtain retarded far-zone dispersion interactions, where only the contribution of low-frequency virtual photons is relevant [13]. A straightworward application of the effective Hamiltonian (26) is the evaluation of dispersion (van der Waals and Casimir-Polder) interactions between two ground-state atoms in the vacuum (zero temperature), even when boundary conditions are present. In this case, we have two atoms, labeled as A and B, respectively, located at R A and R b . The standard quantum electrodynamical calculation of their dispersion interaction involves a fourth-order perturbative calculation [4,21], while the use of the effective Hamiltonian (26) allows us to reduce it to a much simpler second-order calculation. The Hamiltonian of the system is We wish to point out that, although our effective Hamiltonian is a second-order one in the atomic dipole moments (see Equation (5)), in the present case of two atoms we can take the two dipole moments µ A and µ B as independent expansion parameters, and thus (28) is a second-order quantity in µ A and in µ B . Fourth-order terms in µ A and in µ B , which would yield the fourth-order correction to the Lamb shift (not depending from the interatomic distance), would require a fourth-order effective Hamiltonian and are here of course neglected, coherently with our approximations.
Taking into account the form (26) of the effective Hamiltonian, and using (2), we get where α A (ω k ) and α B (ω k ) are, respectively, the dynamic polarizability of atoms A and B. For atoms in the free space and using the mode functions (3), we obtain where R = R A − R B is the distance between the two atoms or molecules (we have assumed real polarization unit vector). Explicit evaluation of (30), using λ (ê kλ ) i (ê kλ ) j = δ ij −k ikj for the polarization sums and in the continuum limit k → (V /(2π) 3 ) d 3 k, yields the well-known dispersion interaction energy between two isotropic atoms as obtained through standard fourth-order perturbation theory, by summing the contributions of the twelve relevant time-ordered Feynman diagrams [4,6,21]. It scales as R −6 in the near zone (nonretarded regime: R ≪ c/ω 0 , ω 0 being a main transition angular frequency of the atom) and as R −7 in the far zone (retarded Casimir-Polder regime, R ≫ c/ω 0 ); thus, the (attractive) force between the two atoms scales as R −7 and R −8 in the near and far zone, respectively. Typical values of the dispersion force between two hydrogen atoms are: ∼10 −30 N for R = 10 −7 m (near zone) and ∼10 −37 N for R = 10 −6 m (far zone). If the appropriate mode functions are used, Equation (29) can be also exploited to evaluate the dispersion interaction in the presence of generic metallic boundaries, for example, a reflecting mirror [36,37]. These results could in principle be extended to generic electrically polarizable bodies. For systems of three or more macroscopic bodies, the non-additivity of dispersion interactions should, however, be taken into account [27,38], and effective Hamiltonians can be very useful in dealing with nonadditive interactions.

IV. CONCLUSIONS
In this paper, we have considered second-order effective Hamiltonians in nonrelativistic quantum electrodynamics, which are very useful to obtain several radiative processes through simpler calculations, compared to the usual evaluations based on the multipolar-coupling Hamiltonian. Examples are energy shifts and Casimir-Polder atomsurface interactions, as well as two-and many-body dispersion interactions between atoms or molecules, in both the nonretarded (near zone) and the retarded (far zone) regime, even in the presence of reflecting boundary conditions. We have reviewed and analyzed different expressions of the effective Hamiltonian, according to the specific system considered and the field states, stressing their limit of validity. The presence or absence of poles in the response function of the atoms has been carefully analyzed, showing that in the case of the photon vacuum the response function has no poles and it is a real function, while in the other cases it is the atomic dynamical polarizability, obeying the Kramers-Krönig dispersion relations. In a forthcoming publication, we will address possible extensions of these results to macroscopic polarizable bodies and Casimir interactions between macroscopic dielectric bodies, stressing in particular their relevance for the dielectric model to be used for the evaluation of such interactions between macroscopic bodies.