Time-Dependent Effective Hamiltonians for Light–Matter Interactions

In this paper, we present a systematic approach to building useful time-dependent effective Hamiltonians in molecular quantum electrodynamics. The method is based on considering part of the system as an open quantum system and choosing a convenient unitary transformation based on the evolution operator. We illustrate our formalism by obtaining four Hamiltonians, each suitable to a different class of applications. We show that we may treat several effects of molecular quantum electrodynamics with a direct first-order perturbation theory. In addition, our effective Hamiltonians shed light on interesting physical aspects that are not explicit when employing more standard approaches. As applications, we discuss three examples: two-photon spontaneous emission, resonance energy transfer, and dispersion interactions.


Introduction
In molecular quantum electrodynamics, atoms and molecules are treated within nonrelativistic quantum mechanics, and the electromagnetic field mediating the interactions is quantized.This approach, born with Dirac's seminal treatment of spontaneous emission [1], is still an ongoing and intense research field, especially with the unprecedented control of light-matter interactions at the atomic scale reached in the last decades.
All phenomena in this topic can be fully understood by starting with the classical minimal coupling Hamiltonian and quantizing it.For molecular quantum electrodynamics, the most convenient approach is to work in the Coulomb gauge.Throughout this work, we shall deal with neutral molecules, in which case we can make a unitary transformation on the minimal coupling Hamiltonian and work with the equivalent multipolar Hamiltonian [2][3][4][5][6].Nonetheless, the generality of these Hamiltonians is also their main weakness, since we must perform extensive calculations to obtain the quantities describing most of the effects.For instance, the interaction between two nonpolar molecules in their ground state results from a tedious fourth-order perturbative calculation.
Here comes the convenience of working with effective Hamiltonians, which are tailored for each specific application, bringing several physical insights and shortening the technical calculation to a much simpler and lower perturbative order analysis.An insightful example is the dynamical polarizability (DP) Hamiltonian, obtained by R. Passante and collaborators [7,8], which is built directly on the molecular dynamical polarizability instead of its electric dipole operator, capturing better the physics governing the interaction.Indeed, nonpolar molecules do not possess permanent electric dipole moments, and their interaction is possible only due to virtual internal transitions that are automatically taken into account by the dynamical polarizability.This is the main message of effective Hamiltonians: building the relevant physical mechanism into the Hamiltonian lowers the perturbative order required in calculations.With the DP Hamiltonian, intermolecular interactions are determined by means of a second-order calculation.Effective Hamiltonians and actions are also useful in describing non-stationary systems and have been employed to develop a multipolar approach to the dynamical Casimir effect [9] and to understand its microscopic origin [10][11][12][13][14].
Effective Hamiltonians are easily constructed from unitary transformations [15], but there is no general recipe to generate useful ones.In this paper, we fill this gap by presenting a systematic method to obtain convenient effective Hamiltonians and by discussing their physical implications.To illustrate our approach, we derive four general effective Hamiltonians allowing us to extend the scope of applications to important phenomena in molecular quantum electrodynamics.
Our method is based on choosing a unitary transformation inspired by (but not equal to) the Hermitian conjugate of the evolution operator for the system.A key concept for our formalism is that the linear susceptibility χ of a quantum system is built from the unequal time commutator of an appropriate operator O describing the system.For the context explored in this paper, we shall take O as being (i) the molecular dipole operator, in which case χ is related to the molecular polarizability, and (ii) the electric field operator, with χ representing the electric field generated by a point dipole (Green function), as discussed in Appendix A.
The importance of the unequal time commutators of the electric field operators in connection to the measurability of the fields was stressed in the literature [16,17].Here, we show that the unequal time commutators can be taken as the basis to generate convenient effective Hamiltonians by allowing a given degree of freedom to effectively dress another one.Physically, this is equivalent to considering part of our system (a molecule or the electromagnetic field) as an open quantum system that is effectively dressed by an appropriate unitary evolution operator, thus yielding an effective time-dependent Hamiltonian for the system.
In Section 2, we employ our formalism to set up the Hamiltonian H eff M , in which the molecular degrees of freedom are dressed by the field.This Hamiltonian generalizes the DP one in two aspects: (i) it naturally accounts for internal dissipation in the molecules, and (ii) it does not require the molecules to remain in the same internal state (usually the ground state) during the process to be described.The latter aspect is a key element and enables us to evaluate the two-photon spontaneous emission (TPSE) in Section 3 within first-order perturbation theory-a much simpler route than the one commonly followed in the literature.
From Section 4 on, we work with situations involving more than one molecule.We employ our method to build the second effective Hamiltonian H eff F , where one molecule, say molecule B, dresses the field.With this dressing, the field acting on the other molecules is given by the superposition of the vacuum electric field with the electric dipole field generated by B. In Section 5, we demonstrate the convenience of H eff F by directly computing the resonance energy transfer (RET) between two quantum emitters in first order.
Then, in Section 6, we analyze dipole-dipole correlation effects and show that different effective Hamiltonians are convenient depending on the distance separating the molecules.In the asymptotic long-distance limit, we demonstrate the Hamiltonian H eff MF , in which the field dresses the molecules, and, in turn, one of the dressed molecules dresses back the field.This Hamiltonian is similar to H eff F , but now the electric field generated by molecule B does not depend on the molecular dipole operator.Instead, it is produced by the dipole induced by the vacuum itself.In the particular case where we assume the molecules to be in the ground state, in the long-distance limit, we recover the Hamiltonian originally proposed by P.W. Milonni [18].
Finally, for the short-distance limit (the non-retarded regime), we follow a complementary route: first, one molecule dresses the field, and then the dressed field dresses back the molecules.This leads us to a new effective Hamiltonian H eff FM .We demonstrate its convenience by applying it in Section 7 to obtain the London interaction energy in first-order perturbation theory and show that this route provides some relevant physical insights.Indeed, our approach can quantify the contributions to the interaction energy coming from the dipole fluctuations of each molecule.The above examples surely do not exhaust the list of useful effective Hamiltonians, and our method should be valuable for several additional applications.

The Field Dresses the Molecules
We begin with a single neutral molecule at position R in the presence of the quantized electromagnetic field.In the dipole approximation, the Hamiltonian describing the system is where d is the molecular dipole operator, and E(R) is the quantized electric field evaluated at the molecule's center-of-mass position R.Note that in the dipole approximation, the electric field can be taken as uniform over the scale of the molecule.H 0 stands for the free Hamiltonians and is given by In this expression a kσ (a † kσ ) stands for the annihilation (creation) operator for a photon with wavector k and polarization σ, whose electromagnetic field oscillates with a frequency ω k = |k|c.H 0m is the free molecular Hamiltonian whose eigenstates and eigenenergies are assumed to be known.The coupling between the molecule and field, given by the second term on the right-hand side of Equation ( 1), can be treated as a perturbation.Therefore, it is convenient to work in the interaction picture with the interaction Hamiltonian The time dependence is obtained by evolving the operators with the free Hamiltonian H 0 .A nonpolar molecule is characterized by not having a permanent electric dipole in its ground state |g⟩, i.e., ⟨g|d|g⟩ = 0. Thus, any process during which the molecule does not excite, such as the Stark shift, must be obtained at least through second-order perturbation theory.If |ψ(t)⟩ symbolizes the state of the molecule-field system in the interaction picture, then its evolution can be written as An equivalent description is generated once we apply a unitary transformation to the state.We choose it as Note that transformation (5) implements the Heisenberg picture to first order in H, thus canceling, at this order, the time evolution of |ψ(t)⟩, which is precisely our goal.If [H(t), H(t ′ )] = 0, then the transformation given in Equation ( 5) would implement the Heisenberg picture exactly.Therefore, in the representation defined by U M , the time evolution of |ψ(t)⟩ results from the non-vanishing value of the commutator [H(t), H(t ′ )], which is consistent with the discussion on linear susceptibilities outlined in Section 1.As will become clear below, this transformation effectively dresses the molecular degree of freedom indicated by the subscript M. Next, we derive the equation satisfied by |ψ M (t)⟩ = U M |ψ(t)⟩.
From Equation (4), we find with Here enters the fact that we are not looking for an equivalent Hamiltonian but rather an effective one.We desire an equivalent Hamiltonian only up to quadratic order in the dipole operator, and thus, we are allowed to expand U M and collect results up to the second order in H, obtaining Note that the expansion in Equation (8) does not correspond to a Dyson series, since the unitary transformation given in Equation ( 5) is not an evolution operator (it lacks a timeordering operator).These expansions differ in the second-order term, and we stress that the third term on the right side of Equation ( 8) is proportional to the square of the second term, which is crucial for the results we will obtain.We define the effective Hamiltonian H eff M (t) as the second-order approximation of H M , which is obtained by substituting the previous relations into Equation (7): where we used the identity 2H(t)H(t Notice that the linear term in the dipole vanished.From Equation (3) and since electric field operators at the same spatial point commute at all times (see Appendix A), we are left with where we employed Einstein notation and denoted by j, l = 1, 2, 3 the Cartesian components of the operators.The great convenience of this Hamiltonian is that it is quadratic in the operators, thus halving the required perturbation order in comparison to the Hamiltonian given by Equation (3).We point out that our demonstration remains the same whether the electric field is quantized or not.As an example, with this effective Hamiltonian, the Stark effect can be obtained from a first-order perturbative calculation.We emphasize that this Hamiltonian is valid only within first-order perturbation theory, but improvements can be made if one keeps extra terms in Equations ( 8) and (9).H eff M mixes both the materials' and fields' degrees of freedom.When the atom is assumed to remain in the ground state, we may take the expectation value of H eff M in the molecular's subspace defined by the ground state through the evaluation of ⟨g|H eff M (t)|g⟩.We stress here that we are not acting on the field subspace, and thus, this average is still an operator in the field variables, which we denote by with being the molecular dynamical polarizability tensor for the ground state describing its linear response to an applied electric field-see Appendix A for details.For practical applications and some physical interpretations, it is generally more suitable to work with the dynamical polarizability in the Fourier representation rather than in the time domain.
To do so, we write the free electromagnetic field in the usual form: where ω k = c|k|, σ is the polarization degree of freedom, and the superscript + (−) refers to positive (negative) frequencies of the field.Substituting Equation ( 14) into (12), we arrive at where stands for the vacuum-induced dipole operator, and Physically, this is the quantum counterpart of the interaction energy of a polarizable system without permanent electric dipole moments in the presence of an external electric field.
In the case without dissipation, the dynamical polarizability in the Fourier space is given by (see Appendix A) where r denotes the excited internal molecular states, and d gr = ⟨g|d|r⟩ = d rg * is the transition dipole moment between states g and r, while ω rg is the corresponding transition frequency.
There are some subtleties concerning the unitary transformation employed in this section.One could argue that once the integration present in Equation ( 5) starts from −∞, our truncation in Equation ( 9) is not rigorous.The point is that the molecule has a finite memory, characterized by a time scale τ.This means that α jl (t − t ′ ) vanishes for t − t ′ ≫ τ in Equation (13), enabling the lower limit in (5) to be replaced by t − τ.The validity of this truncation is then tantamount to the validity of the perturbative method in the moleculefield interaction, justifying our approach.This argument also underlies the convenience of working in the Fourier space.Convergence of the time integration in Equation ( 12) requires that we account for dissipation in the polarizability.Nevertheless, in many cases of interest, the most relevant Fourier modes are far from molecular resonances, and we may neglect dissipation when using Equations ( 15) and (16).
Another important aspect is that the effective Hamiltonian ( 11) is convenient only when first-order perturbation theory in Hamiltonian (3) vanishes, even though regularization techniques may render it applicable if this is not the case [8].We may separate the main applications of the effective Hamiltonian H eff M into two groups: (i) the molecule remains in the same internal state during the entire process, and the expectation value of the electric dipole operator in this state is zero; (ii) the molecule undergoes a transition between two internal states, but the electric dipole operator is unable to connect these two states.Examples involving (i) have already been discussed in the literature [7,8], in contrast with case (ii).One fascinating example of this second group is the two-photon spontaneous emission, with selection rules forbidding the one-photon transition.In the next section, we explore this example from the perspective of the effective Hamiltonian derived in this section.

Application to the Two-Photon Spontaneous Emission
An excited molecule may decay to its ground state through the emission of two photons in the so-called two-photon spontaneous emission (TPSE).This phenomenon is particularly interesting when the one-photon transition is forbidden.The TPSE makes the vacuum unstable and is responsible for the initial buildup of the intracavity field in two-photon micromasers [19][20][21].More recently, it was shown that the simultaneously emitted photons can be indistinguishable and entangled in time and frequency [22][23][24], renewing the interest in this phenomenon [25][26][27][28][29].This section aims to obtain the TPSE rate directly from first-order perturbation theory in Hamiltonian (11).Let us consider that a molecule in an internal state |e⟩ decays in vacuum to its ground state |g⟩ through the emission of two photons with wavevectors k and k ′ and polarizations σ and σ ′ .To this end, it suffices to analyze the matrix element of H eff M connecting the initial and the final states, given by ⟨g; where we chose the origin of our coordinate system at the position of the molecule.Here, |0⟩ denotes the vacuum state of the electromagnetic field.We also define which involves only the molecular degrees of freedom and quantifies the linear response of the molecule to an applied field connecting internal states |e⟩ and |g⟩.Note that, when taking |e⟩ = |g⟩ in Equation ( 20), the tensor ← → D yields as a particular case the polarizability of the molecule, which is given by Equation (13).The TPSE rate is immediately obtained in the long-time limit by substituting Equation ( 19) into Fermi's golden rule.In general, it is more convenient to represent ← → D in Fourier space.We begin by writing and, at the end, we take t 0 → −∞.In this expression, H 0m denotes the free molecular Hamiltonian, with eigenstates satisfying H 0m |r⟩ = hω r |r⟩, so that by inserting a closure relation I = ∑ r |r⟩⟨r| into Equation (20), we obtain The instant t 0 plays a role only in an unimportant global phase, which was discarded.
In Fourier space, Equation (23) becomes The desired matrix element given in Equation ( 19) can be obtained from Equation (14).For emission processes, only positive frequency modes contribute.Using Equation ( 22), we also obtain ⟨g; where ϵ j kσ is the jth Cartesian component of the polarization unit vector for the mode with wavevector k and polarization σ, and V is the volume of the quantization box.In the long-time limit, we are interested in photon pairs satisfying the condition , where we used Equation (24) in the last equality.Therefore, we may simplify Equation (25) to ⟨g; The probability rate of emitting one photon in a solid angle dΩ around k and with a frequency in the interval (ω, ω + dω) and another in a solid angle dΩ ′ around k′ , as well as with a frequency in the interval (ω ′ , ω ′ + dω ′ ), is given by (from now on we denote Employing Fermi's golden rule, we arrive at Integration over the solid angles may be readily evaluated from the identity We also integrate over ω ′ to find the photon emission rate: which is equivalent to the result of Ref. [2].When performing second-order perturbation theory, the usual notation is to describe the molecular response in terms not of α ge -which is a function of a single frequency variable, but rather a function of two frequency variables, which are obtained from Equation ( 24) by replacing ω re + ω by ω rg − ω ′ in the second term.The calculation from the new effective Hamiltonian (11) not only yields the two-photon spontaneous emission rate with a much shorter first-order calculation but also describes the results in terms of a single frequency variable function that sheds an interesting light on the physical mechanism involved in the phenomenon.In order to unveil the physical significance of α ge jl , let us project the Hamiltonian H eff M into the field's Hilbert space, thus generalizing Equation ( 12) for situations where the molecule undergoes an internal transition.This is done by defining the new effective Hamiltonian from Equation ( 11): which involves only electric field operators.Note that D ge jl (t, t ′ ), given by Equation ( 20), is not a real number, and, therefore, H eff(ge) M is non-Hermitian.This is due to the fact that the field degrees of freedom alone constitute an open quantum system, extracting energy from the drive provided by the molecular internal transitions encapsulated in D ge jl (t, t ′ ).The non-hermiticity of Hamiltonian (31) also reflects the break of time inversion symmetry imposed by the two-photon decay.
Following the same steps that led us from Equations ( 12) to (15), we obtain where the induced dipole for the transition |e⟩ −→ |g⟩ is given by For this reason, we denominate α ge jl as the transition polarizability tensor.This is a useful concept whenever the transition dipole element of a given internal molecular transition vanishes but can be induced by an external electric field.It generalizes the concept of the polarizability tensor, which stands for the dipole induced for a fixed internal molecular state.This induced transition dipole acts as an external source oscillating with frequency ω eg and driving the field appearing in the effective Hamiltonian (31).Here, ω eg > 0 indicates that energy-conserving processes must be accompanied by photon creation, as can be verified in Equation (14).In this case, only the last term in Equation (33) contributes to the process.The other term is relevant for two-photon absorption, and the calculation presented in the section applies with minor modifications to this case.

The Molecules Dress the Field
In the previous section, we investigated the convenience of employing effective Hamiltonians in which the electric field dresses the molecules.Now, we shall analyze the opposite case and present an effective Hamiltonian that describes a molecule dressing the electric field operator.Consider two nonpolar molecules A and B. The electric dipole Hamiltonian describing this system in the interaction picture is where and d ζ is the electric dipole operator of molecule ζ = A, B, whose center of mass is at position R ζ .We again represent the system's state with |ψ(t)⟩, satisfying Equation ( 4), but implicitly including a tensor product of both the molecules' and fields' states.We follow the same reasoning as in the previous section.In this case, however, we want the molecule B to dress the electric field operator.Hence, we choose as the unitary transformation the inverse of the evolution operator for the coupling between molecule B and the field: where T is the anti-time ordering operator (earlier-time operators on the left).Its presence implies a crucial difference in comparison with Equation (5), and its purpose is to eliminate H B so that the entire role played by molecule B will be through the field it produces.
If only molecule B were present, the unitary transformation U F would take the interaction picture into the Heisenberg picture.Nonetheless, in the presence of atom A, this unitary transformation yields a new effective Hamiltonian, to which we now turn our attention.Following steps analogous to those in Section 2, the equivalent Hamiltonian is given by Due to the anti-time ordering operator, we have ih∂ t U F = −U F H B , and, thus, canceling H B as anticipated.Expanding up to the linear term in d B , we obtain (see Appendix A) where with r = R A − R B and t r = t − r/c being the retarded time.This expression corresponds to the electric field generated by dipole B at the position of molecule A. This result is readily extended to any number of molecules by exchanging 35), thus obtaining While Equation ( 37) is exact and constitutes an equivalent Hamiltonian, Equation ( 40) is effective and valid only up to linear order in d B , d C , etc.It is worth mentioning that Equation ( 40) can be generalized to other situations.For instance, if atom A is in the presence of a magnetically polarizable atom [30,31], we have to add the electric field produced by the magnetic dipole of atom B in Equation (38).In the next section, we demonstrate the convenience of the new effective Hamiltonian H eff F by obtaining the RET rate in a first-order calculation.

Application to the Resonance Energy Transfer
In a resonance energy transfer (RET) process, an excited molecule decays through nonradiative channels, transferring its energy to a molecule in the ground state [32][33][34][35][36][37][38][39][40][41].This phenomenon is of notable importance to many areas of science due to its broad range of applications across fields such as chemistry [42], medicine [43], and biology [44].Throughout this section, we discuss the probability that an excited molecule A decays, exciting an identical molecule B that was initially in its ground state, placed at a distance r from A, with both in vacuum.
Up to the second order in perturbation theory, the probability amplitude of interest can be calculated as [15] In this expression, |ψ i ⟩ = |e A , g B , 0 kσ ⟩ (with energy E i ) and |ψ f ⟩ = |g A , e B , 0 kσ ⟩ describe, respectively, the system's initial and final states, ĤInt is the interaction Hamiltonian, and |ψ r ⟩ are the intermediate states with energy E r .In the standard approach, the interaction Hamiltonian is taken as the dipolar Hamiltonian given by Equation ( 34): H int = H (2) .With this choice, however, the first term in Equation ( 41) vanishes, and the RET rate is obtained from second-order perturbation theory.Here, we offer an alternative and simpler approach by letting atom B dress the field and taking H int = H eff F , as in Equation (40).In this case, it suffices to calculate the first-order matrix element Following Equation ( 21), the first matrix element on the right-hand side of the previous equation becomes (more precisely, we should consider the evolution beginning at time t 0 ; however, as explained in Section 3, this would only contribute as an irrelevant global phase) and, by using Equation ( 39), the terms contained in the second matrix element give the contributions ⟨e where k = ω eg /c.Replacing these results in Equation ( 42), we arrive at By applying Fermi's golden rule, where ρ(ω f ) is the density of final states with energy E f = hω f , and we directly recover the well-known result for the RET rate [45,46].

The Dressed Molecules Dress the Field-And the Reverse
The first-order perturbation in the Hamiltonian (40) vanishes whenever the electric dipole operator of one of the molecules cannot connect the involved molecular states.An example is the force between molecules in their ground state to be analyzed in the next section.Here, instead, we focus on a general discussion without specifying the molecular internal state.The physical mechanism that limits the dipole-dipole correlation depends strongly on the distance R separating the molecules.Indeed, two characteristic time scales are key to understanding the two different regimes: the time it takes light to travel between the molecules, t γ = r/c, and the characteristic time for dipole fluctuations, t d = 1/ω 0 , where ω 0 is a typical transition frequency for the molecules.In the asymptotic longdistance regime, t γ ≫ t d , it is the electrodynamical retardation that limits the dipole-dipole correlation, and we may neglect dispersion in the atomic response.In the opposite shortdistance regime, electrodynamical retardation is negligible, and it is now the delay in the molecular response that limits the dipole-dipole correlation.Now, the molecular dispersion is crucial, but we can take the electrostatic limit for the electric field produced by each electric dipole.To go deeper into the physical particularities of these two complementary regimes, we shall develop a different effective Hamiltonian appropriate to each case.

The Dressed Molecules Dress the Field: Retarded Long-Distance Regime
In the long-distance regime, we may neglect dispersion in the molecules, which is tantamount to considering an instantaneous molecular response.This means that the time-scale variation for the electric field is much slower than the molecular response, enabling us to approximate when evaluating the effective Hamiltonian H eff M (t) given by Equation (11).We have defined the molecular operator Π jl as From Equation ( 13), we see that the expectation value ⟨g|Π jl |g⟩ is the static polarizability of the molecule in its ground state, given by setting ω = 0 in Equation ( 18).The tensor operator Π jl generalizes this concept by enabling us to capture the static response even for processes involving changes in the molecular internal state.Substituting ( 48) into (11) leads to which corresponds to the static response limit of H eff M (t).In the case of two molecules, dipole-dipole correlations arise in second-order perturbation theory in the Hamiltonian where Π ζ,jl is the operator (49) for molecule ζ = A, B.
An equivalent Hamiltonian where the molecules couple directly with each other will be able to capture the dipole-dipole correlation in first-order perturbation theory.This can be done by employing the unitary transformation which mimics the one employed in Section 4, with the difference that it is the dressed molecule (through operator Π jl ), instead of the naked molecule, that dresses the field.Following the same steps that led us from Equation (34) for molecule B into Equation (37), we get H eff,s MF is very suitable for handling effects related to the interaction between atoms A and B because the expansion of U MF in Equation ( 53) contains terms combining the product Π A,jl Π B,mn .Such terms also appear through a fourth-order perturbation theory in Hamiltonian (34) but already appear in first order here.To obtain H eff,s MF , it is enough to implement the transformation rule for the electric field at R A , since Π A,jl commutes with Π B,mn .Substituting Equation ( 50) into (52), we obtain the following up to linear order in Π B (see Appendix A): where E ind dip,B (R A , t) is given by Equation ( 39) with the substitution d B,j (t r ) −→ Π B,jk E k (R B , t r ).This result is mathematically similar to Equation ( 38) but with a remarkable physical difference.Here, the field is dressed not by a naked molecule but by a dressed one.This means that the source of the electric dipole field is not the molecular dipole operator but rather a vacuum-induced dipole-as indicated by the superscript "ind" in Equation (54).From Equations ( 53) and (54), we arrive at where we have kept only the terms capturing the dipole-dipole correlation and neglected the higher-order term Π A,jk E ind dip,B,k (R A , t)E ind dip,B,j (R A , t).This expression is manifestly symmetric upon the exchange A ←→ B, as can be verified by substituting Equation ( 39) into (55).
As an application of the effective Hamiltonian (55), we may consider that both molecules are at their ground states throughout the entire process.In this case, we may take the expectation value of H eff,s MF (t) in the ground states of the molecules, which is equivalent to substituting the operators Π jl with the static polarizability tensor of the corresponding molecule in Equation ( 55).In the particular case of isotropic molecules, this reproduces the asymptotic long-distance limit of the effective Hamiltonian originally employed by P.W. Milonni, which readily yields the known Casimir-Polder result in first order [18].For comparison, when taking the average of H eff,s M (2) (see Equation ( 51)) over the molecular ground state, the resulting effective Hamiltonian [47] yields the Casimir-Polder energy only to the second order of perturbation theory.
In the opposite short-distance limit, molecular dispersion is essential, and we cannot make the approximation given by Equation (48).In this case, it is more convenient to work with a different effective Hamiltonian, which we shall present in the next subsection.

The Dressed Field Dresses the Molecules: Non-Retarded Regime
We now turn to the opposite regime, in which the intermolecular distance is so small that we may neglect the electromagnetic retardation in comparison to the molecular response time.Unlike the other examples in this paper, this case does not require quantization of the electromagnetic field.On the other hand, the molecular dispersion is crucial in this regime.The dipole-dipole correlation is usually obtained from a second-order perturbation theory in the dipole-dipole Hamiltonian [48] where r = rr is the position of molecule B with respect to molecule A. Notice that this Hamiltonian is a particular case of Equation ( 40) without the vacuum electric field operator and with the dipole electric field taken in the electrostatic approximation.In the nonretarded regime, the field and the molecular operators switch roles with respect to what happens in the retarded asymptotic regime.While, in the latter, we could begin with Hamiltonian (48), which approximates the molecular response with its dressed static response, Equation ( 56) is precisely the opposite: now it is the field that is dressed by its static response.Indeed, the electrostatic dipole field leading to (56) corresponds to the zero-frequency limit of the Green function of the wave equation, which plays the role of the field susceptibility, as discussed in Appendix B.
Mirroring the procedure of the previous subsection, we let the already-dressed field dress the molecules, thus leading to a new effective Hamiltonian.In the long-distance regime, we employed a unitary transformation extending the formalism of Section 4. Here, on the other hand, we want to extend the formalism developed in Section 2. As a first step, we write Hamiltonian (56) in the interaction picture and then take as the unitary transformation the operator which should be compared to Equation (5).Following steps analogous to the ones leading to Equation ( 10) yields Since operators involving different molecules commute, we have By substituting Equation ( 56) and the previous identity into Equation (58), we obtain the effective Hamiltonian The new effective Hamiltonian (60) has two terms that capture the physics involved in the dipole-dipole correlation.The product [d A,j (t), d A,m (t ′ )]{d B,k (t), d B,n (t ′ )} measures how the dipole fluctuations of molecule B induces a dipole in molecule A, while the other term is its reciprocal.This decomposition is possible because, differently from the standard approach based on second-order perturbation theory with the time-independent dipoledipole Hamiltonian (56) where the two molecules are considered as an isolated system, here, we take the complementary approach of considering each molecule separately as an open quantum system.This perspective offers two main novelties: (i) H eff FM brings to light the dynamical character of the dispersion interaction by making an explicit connection with dipole fluctuations, and (ii) H eff FM enables us to assess the contribution from the fluctuations of each molecule separately.In the next section, we analyze an example that illustrates these advantages.

Application to the London Interaction Energy
In this section, we consider the dispersion interaction between two ground-state nonpolar molecules A and B in vacuum, which interact due to correlations between their fluctuating electric dipoles.As discussed in the previous section, the physical mechanism limiting the dipole-dipole correlation strongly depends on comparing the distance separating the molecules and their internal transition wavelengths.For ground-state molecules, the resulting intermolecular interaction energy exhibits a different power-law dependence with distance in each of the two regimes discussed in Section 6.
As originally demonstrated by London [49], the non-retarded interaction energy can be obtained without quantizing the electromagnetic field and scales with 1/r 6 .The asymptotic long-distance limit was first obtained in the seminal paper by Casimir and Polder [50], where they showed that retardation imposes the necessity of quantizing the electromagnetic field and demonstrated that the interaction energy scales asymptotically with 1/r 7 .
In the previous section, we discussed how the Casimir-Polder result for the asymptotic long-distance limit can be derived from the effective Hamiltonian H eff MF (55).In this section, we obtain the London result for the short-distance non-retarded limit from the new effective Hamiltonian H eff FM (60), which provides new physical insights into the dipole-dipole correlation present in the non-retarded regime.By taking the average of Hamiltonian (60) in the ground state of the molecules and employing the result (13) for the molecular polarizability tensors α A jl , α B jl , we obtain where we defined the symmetrical dipole correlation function with ζ = A, B. To work in the Fourier space, we can apply Parseval's theorem, so Equation ( 61) becomes where we used that η A,B (ω) are real functions, as they are Fourier transforms of real even functions.This result is the analog for two molecules of the decomposition obtained for an atom coupled to the vacuum electric field [115][116][117].In the latter, the field susceptibility captures the field radiation reaction.More recently, an analogous decomposition was also obtained for atoms interacting with a scalar quantum field [118].A decomposition similar to (63) was employed to derive a nonlocal phase for a moving atom interacting with a planar surface [119] and a Sagnac-like atomic phase induced by a rotating nanosphere [120].
In the isotropic case, the polarizability tensors simplify to α A(B) rs (τ) = δ rs α A(B) (τ), and the symmetric correlation functions simplify to η A(B) rs (τ) = δ rs η A(B) (τ).Then, Equation (63) leads to In Appendix B, we employ the analytical properties of the correlation functions to demonstrate that our results are equivalent to the standard way of expressing the London interaction energy for any molecular model of the polarizabilities.Here, we show the convenience of Equation ( 64) by considering the simple case of two-level atoms, for which (ζ = A, B) [18] Let us analyze each contribution to the London interaction energy in Equation ( 64) separately.We define as the contribution arising from the dipole induced at atom B by the dipole fluctuations of atom A. From Equations ( 65) and ( 66), we obtain The interaction energy is E London = E B→A London + E A→B London , with E B→A London being obtained by interchanging the roles of A and B in Equation (68).
Let us consider ω 0A > ω 0B .In that case, from Equation (68), we see that E A→B London > 0, indicating that the dipole induced at a slower atom by a faster one generates a repulsive contribution to the dispersion force.This is due to the fact that the polarizability given by Equation (65) becomes negative for frequencies higher than the atomic transition frequency.Indeed, the induced dipole at the slower atom B cannot follow the fast oscillation of the fluctuating dipole of atom A. The induced dipole at B lags behind the field of atom A and points opposite to its direction at a given time.As a consequence, the induced dipole at B repels the fluctuating dipole at A. However, the opposite holds for the complementary term E B→A London : the dipole induced in the faster atom A can follow the dipole fluctuations of atom B in phase, leading to an attractive contribution.Attraction overcomes repulsion by a factor ω 0A /ω 0B , since the slower atom couples less effectively to the field than the faster one.If ω 0A = ω 0B , each contribution diverges due to a resonant response.This divergence would be avoided if dissipation were taken into account.Nevertheless, it is remarkable that the divergence cancels once we sum E B→A  London and E A→B London , leaving us with the well-behaved total interaction energy which agrees with the result [2] calculated from second-order perturbation theory based on the dipole-dipole Hamiltonian (56).Notice that varying ω 0B while keeping the other parameters fixed shows that the attraction is maximal when ω 0B → ∞.The previous decomposition clearly illustrates the physical mechanism involved.From Equation (68), we see that in this limit, the repulsive contribution E B→A  London vanishes, indicating that atom A is effectively transparent, decoupling from the rapid oscillating field produced by B. The attractive term in Equation (68), on the other hand, takes its maximal absolute value in this limit, since the response of atom B is so fast that it perfectly mirrors the fluctuations of the other atom.In this sense, we may conclude that atom B in the limit ω 0B → ∞ is the atomic analog of a perfect conductor.
As was true with the other effective Hamiltonians discussed in this paper, we see that the convenience of employing H eff FM is twofold: (i) it lowers the perturbation order required to obtain the London dispersion energy from second to first order, and (ii) it offers physical insights into the mechanisms involved in the phenomenon.The results in this section can be readily extended for multilevel atoms.To this end, it suffices to substitute Equations ( 65) and ( 66) with a summation over all internal transition frequencies.

Final Remarks and Conclusions
All phenomena in molecular quantum electrodynamics can be obtained from the multipolar Hamiltonian.In this paper, we restricted our attention to phenomena that can be treated perturbatively (which includes the vast majority of cases in this field).In most situations, the dominant effect is obtained from a high-order perturbation theory, requiring intermediate states to connect the initial and the final states.A clear example is the interatomic interaction.While in classical electrodynamics, we may always take the field at each charged particle as the superposition of the field generated by all other particles, in standard quantum electrodynamics, each particle couples only to the free electromagnetic field.Consequently, we must go up to the fourth order to obtain the dominant contribution when considering molecules without permanent electric dipole moments.An alternative is to build effective Hamiltonians.They are customized for each specific application and lose meaning and validity after some point in the perturbative expansion.Not only do they greatly simplify the technical difficulties involved in calculations using the multipolar Hamiltonian, but, equally importantly, the effective Hamiltonians cast the phenomena in a new light, offering insightful physical interpretations.
Several effective Hamiltonians have successfully been employed by many authors in the last decades.In this paper, we have developed a systematic approach to constructing effective Hamiltonians, which allowed us to derive a number of new ones, choosing as a unitary transformation the Hermitian conjugate of the evolution operator for part of the system.This transfers part of the time evolution from the vector state to the operators, dressing them and providing a Hamiltonian that requires a lower order in perturbation theory to account for the process of interest.This method can always be used when the firstorder perturbation theory vanishes.Our approach yields time-dependent Hamiltonians that enable us to follow the energy exchange between matter and the field, with each subsystem constituting an open quantum system.We emphasize here that our system of interest is the entire molecule-field system, and it is not interesting to trace over any of the subsystems, as is common in an open quantum system approach [121][122][123].
As a first application, we have derived the Hamiltonian H eff M , where the field dresses the molecule and the dipole operators are replaced by their commutator at different times.If we project the commutator in an internal molecular state, it yields the dynamical polarizability of the molecule for the corresponding state.Its nondiagonal elements, on the other hand, allow the dressing to leave the molecule in a final state that is different from the initial one.We have demonstrated that this new time-dependent Hamiltonian provides a simpler treatment of the two-photon spontaneous emission, as the dominant contribution is obtained in first-order perturbation theory.In addition, our formalism introduces the concept of an induced dipole transition, which generalizes the notion of an induced dipole for a given internal state.
Then, we discussed applications involving two molecules A and B. We constructed the new effective Hamiltonian H eff F through a unitary transformation that transfers all of the effects related to molecule B to the electric field it generates.In this way, molecule A feels an effectively dressed electric field given by the superposition of the free vacuum electric field, and the one generated by the dipole operator of molecule B. H eff F allows for the description of the resonance energy transfer rate in first-order perturbation theory.
Lastly, we derived two additional Hamiltonians that merge aspects of the previous two, where each one is appropriate for a different intermolecular distance regime.In the asymptotic long-distance regime, molecular dispersion is negligible, enabling us to derive an effective Hamiltonian H eff MF which is formally similar to H eff M .In this new case, however, the field acting on molecule A is given by the superposition of the free electric field and the one produced by the vacuum-induced dipole generated on molecule B. When we average H eff MF over the molecular ground state, we re-obtain the asymptotic limit of the Hamiltonian employed by P. Milonni [18].
Finally, for the short-distance non-retarded limit, we derived our fourth and last effective Hamiltonian H eff FM based on the fact that, in this limit, we do not need to quantize the electric field.This effective Hamiltonian enables us to clearly identify the different physical mechanisms involved in the correlations responsible for the interaction, separating one term where the dipole fluctuations of molecular A induce a dipole on molecule B and another term where the roles are exchanged.
As a final application, we employed H eff FM to obtain the London dispersion interaction energy in first-order perturbation theory.We showed that, for two-level atoms, the dipole fluctuations of the atom with the higher transition frequency give rise to a repulsive term, since its fast fluctuations cannot be followed by the slower atom.Nonetheless, the force between two isotropic atoms is always attractive, since the fluctuations of the slower atom are strongly correlated and easily followed by the faster atom, overcoming the repulsive contribution.The possibility of quantitatively and separately analyzing the contributions arising from each mechanism correlating fluctuating systems is an advantage of H eff FM .The Hamiltonians presented in this paper can be employed in a great variety of situations.For instance, one may treat the effects of boundaries in the two-photon spontaneous emission or resonance energy transfer by simply introducing the appropriate field modes.As in the examples discussed in this paper, these effective Hamiltonians allow for a direct first-order calculation within perturbation theory.More notably, the methodology introduced here can be applied to generate other effective Hamiltonians that may optimize calculations and provide physical intuition.
[7]iω k t is the l-th Cartesian component of E Notice that d ind (t) acts on the field's Hilbert space.Due to the reality of α(t − t ′ ), α jl (−ω k ) = α * jl (ω k ), and, thus, d ind (t) is an Hermitian operator.If dissipation is negligible, we re-obtain as a particular case the DP Hamiltonian[7]