1. Introduction
For quite a few decades, Franco Persico admirably led the theoretical physics division at the University of Palermo in Sicily. Together with his students and colleagues, he made notable advances in the areas of quantum optics, strong field–matter interactions, many-body forces and dressed states [
1,
2,
3,
4]. An often-employed theory that has yielded solutions to many fundamental problems, including those involving atoms and molecules, is quantum electrodynamics (QED) [
5,
6,
7,
8]. Electromagnetic radiation and matter are treated equally in this framework, with both being subject to the principles of quantum mechanics.
A particularly interesting phenomenon that has been successfully explained by QED theory is the radiation field-induced shift in the dispersion energy between two atoms or molecules, colloquially termed the optical binding energy [
9,
10,
11]. Thirunamachandran first predicted this in 1980 [
12]. He used QED to calculate the potential energy of a pair of interacting particles in the ground state when subject to an intense laser field. This process provides an example featuring both real and virtual photons. The former originate from the applied radiation field and are scattered by the two objects, which are themselves coupled through the exchange of a single virtual photon. Experimentally, this manifests as the mutual attraction of neutral polarisable particles located in the focal area of a beam of light, in the process generating optomechanical forces and torques via optical spanner or wrench techniques that lead to the manipulation of pairs of particles [
13], extending the earlier trapping, control and movement of individual entities via optical tweezer methods which culminated in the award of the 2018 Nobel Prize in Physics to Ashkin [
14].
The original calculation of the change in the energy shift was carried out using time-dependent perturbation theory and involved the summation over 48 diagrams of a so-called “dynamic” mechanism, where off-resonant laser light is absorbed and emitted at different centres, and the later evaluation of an additional 48 diagrams required when the real photons are scattered at the same site in a second so-called “static” mechanism [
15].
The QED theory of this effect has also been extended to describe the optical binding energy between a pair of optically active molecules [
16,
17,
18]. To account for the discriminatory change in the potential of such a system required dropping the electric dipole approximation and including both the electric and magnetic dipole coupling terms. The interference term between these two moments yields a contribution to the energy shift that is dependent on the chirality of the molecule that changes sign on replacing one enantiomer with its mirror-image structure. A diagrammatic time-dependent perturbation theory calculation of the change in the energy shift is necessarily prohibitive for the case of two chiral molecules due to the large number (192) of diagrams that need to be drawn, individually evaluated and finally summed. The result was instead obtained by calculating the coupling of the dipoles induced at each site by the time-varying electromagnetic fields. Similar to London’s view of the origin of the dispersion force, this method was adopted by Power and Thirunamachandran to evaluate the Casimir–Polder potential [
19] and applied later to the pair radiation-induced interaction energy [
20,
21], and the optical binding energy in three- and
N-body systems in the electric dipole approximation [
22].
An often-overlooked contribution to many processes is the electric octupole coupling term. While of a higher order relative to the magnetic dipole and electric quadrupole interaction terms, when the octupole moment,
, is decomposed into its irreducible components of weight-1 and weight-3,
and
, respectively, the former has the properties of a vector and may be viewed as a higher-order correction to the electric dipole interaction term. Its effect on two- [
23,
24] and three- [
25,
26,
27] body dispersion potentials and the rate of resonance energy transfer [
28,
29] has been shown previously.
This partitioning may be written as
where
and
with
e denoting the elementary charge,
being the
i-th Cartesian component of the electron coordinate, and Einstein’s summation convention is employed for Latin indices in the space-fixed frame of reference or Greek subscripts in the molecular frame that repeat.
δij is the Kronecker delta.
We extend the coupling of fluctuating moments at each centre to calculate the radiation field-induced energy shift between a pair of molecules involving electric dipole and octupole coupling terms and compare results with purely electric dipole potentials obtained previously [
6,
11,
12] to see if higher-order corrections occur and what explicit form they take. A complete polarisation analysis is also carried out for linearly and circularly polarised light propagating in directions parallel and perpendicular to the inter-nuclear separation distance vector to aid experimental investigation. Results are also obtained for fully isotropic situations applicable to species in the fluid phase.
2. Optical Binding Energy between an Electric Dipole-Polarisable Molecule and a Mixed Dipole–Octupole-Polarisable Molecule
The picture to be adopted in the calculation of field-induced inter-particle energy shifts involving molecules dependent upon the electric octupole coupling term is the one in which an applied electric field induces an electric multipole moment at a specific centre. This depends on the polarisability characteristics of the reacting species. The induced moments interact via their mutual multipole–multipole coupling tensor. The energy shift is then obtained on taking the expectation value for both particles in the ground electronic state and the radiation field in a state corresponding to a laser containing
N photons. Results for the Casimir–van der Waals dispersion potential are obtained immediately on letting
N → 0, which were previously calculated using diagrammatic perturbation theory techniques [
30], serving as a useful check of the work.
The first octupole-dependent case to be considered is the radiation field-induced interaction energy between an electric dipole-polarisable molecule,
A, positioned at
, and a mixed electric dipole–octupole-polarisable molecule,
B, located at
. In the presence of a transverse electric displacement field,
, the electric dipole moment induced at
A,
, is
(with
denoting the permittivity of the free space), where the dynamic pure electric dipole polarisability tensor of
A, at the frequency
, where
k is the magnitude of the wave vector, is given by
In Equation (5),
is the transition electric dipole moment matrix element between the ground state |0> and the virtual state |
a> of
A,
, with similar expressions for other transition multipole moments, with |
b> denoting the virtual electronic states of
B. In the denominator of Equation (5), the energy difference between states |0> and |
a> is
. Due to the polarisability characteristics of
B, two different multipole moments are induced there, namely the electric dipole moment
and the electric octupole moment
where the frequency-dependent mixed electric dipole–octupole polarisability tensor is defined as
The multipole moments induced at each site couple via the relevant retarded multipole–multipole interaction tensor. For the coupling of an electric multipole moment of order
m at
A and order
n at
B, its generalised form is [
29]
where the tensor
is
with
Thus the two induced electric dipoles couple through the retarded electric dipole–dipole interaction tensor,
, and the induced dipole at
A and the induced octupole at
B interact via the retarded electric dipole–octupole coupling tensor,
giving rise to the following formula for the energy shift:
with
and
Substituting Equations (4), (6) and (7) into Equation (11) provides an expression of the energy shift explicitly in terms of the molecular polarisabilities (5) and (8) of
A and
B,
Next, the expectation value of Formula (14) is taken over the state
corresponding to both molecules in the electronic ground state and the field containing
N photons of mode
, where
is the direction of propagation of the applied radiation field and
λ is its index of polarisation. For the matter sub-system, this yields the ground state polarisabilities (5) and (8). For the radiation field, use is made of the microscopic second quantised transverse electric displacement field operator of QED [
31,
32], whose form for a specific mode is
where
is the complex unit electric polarisation vector,
is the photon annihilation operator,
is the creation operator and
V is the volume of the quantisation box. Employing Equation (16) leads straightforwardly to the following expectation value for the product of the
m-th- and
n-th-order gradients of the electric displacement field evaluated at two spatially distinct points for a state containing
photons:
Concentrating on the first term of the energy shift (14), the expectation value of the field-dependent factor is readily obtained from Equation (17),
which, after carrying out the polarisation sum using the identity [
5]
yields
2.1. Orientationally Averaged Energy Shift
To obtain the isotropic interaction potential, a pair orientational average of the radiation field vectors relative to the inter-particle displacement vector must be performed. The contributing averages, denoted by angular brackets, may be extracted from the generalised relation
so that the average featuring in Equation (20) is
where from Equation (10)
Thus, contribution (20) becomes
on defining the intensity of the incident field as
On substituting Equation (12), Equation (24) becomes
The rotational averages for the two particles are given by the formulae [
33]
and
Note that a factor of 2/15 has been absorbed into the isotropic mixed dipole–octupole polarisability,
In these and subsequent results, Greek subscripts refer to Cartesian tensor components in the body-fixed frame of reference. Interestingly, there is no contribution to the orientationally averaged mixed electric dipole–octupole polarizability (8), due to the octupole weight-3 term (3). After contracting tensors (25) becomes
Evaluating
and
,
and
noting that
Equation (28) becomes
Returning to Equation (14) and calculating the second term, the field part from Equation (18) is
Substituting Equation (32) and employing Equation (19) produces
From Equation (21), the pair orientational average is
and is given by Equation (30). With
given by Equation (13), Equation (33) becomes
which, after deploying Equations (26) and (27), becomes
and a term identical to Equation (31) results. Hence, the isotropic energy shift induced by an external electric field between an electric dipole-polarisable molecule and a mixed electric dipole–octupole-polarisable one is twice expression (31),
It is interesting to compare expression (37) with the analogous result between a pair of electrically-polarisable molecules [
12]:
It can be seen that the coefficients of the terms within square brackets match each other and agree with those occurring in the Casimir–Polder dispersion potential, allowing Equation (37) to be viewed as a higher-order correction to the leading pure electric dipole contribution to the radiation-induced dispersion energy. Asymptotic limits of the result Equation (37) readily follow. In the near zone,
kR << 1, yielding an inverse separation distance dependence
while at the opposite extreme of the far zone,
kR >> 1, resulting in a modulated inverse square behaviour on relative separation
2.2. Polarisation Analysis
It is possible to obtain expressions for the change in energy shift that illustrate the dependence on different polarisation states and configurations of the incoming radiation field. This offers enhanced possibilities for experimental control and detection. Going back to the first term of Equation (14) and inserting relation (18) produces
For large occupation numbers of the photon field, it is sufficient to approximate
N + 1 by
N so that Equation (41) becomes, after molecular orientational averaging,
One possibility is for the propagation direction of the laser to be parallel to the inter-molecular join, resulting in the polarisation vector being perpendicular to Then, A second possibility is for the laser to be oriented orthogonal to in which case the polarisation vector is parallel to the inter-particle separation distance vector. This results in In addition, the laser may be linearly or circularly polarised. We examine each of these four possible combinations for each of the examples considered.
- (a)
Linearly polarised light.
- (i)
In this scenario, Equation (42) becomes
giving rise to near- and far-zone limiting forms
and
displaying
R−3 and (cos
2kR)/
R dependences, respectively.
- (ii)
Expression (42) then becomes
from which inverse cubic short-range and modulated inverse square long-range dependences on
R are found.
- (b)
Circularly polarised light.
To proceed, use is made of the relation
where the upper and lower signs refer to left (L)- and right (R)-handed circular polarisations, and
is the Levi-Civita tensor. Note that only the
i,
j-symmetric part of (47) persists, since the geometric tensor appearing in Equation (42) is
i,
j-symmetric.
- (i)
For this case, an identical contribution to that arising from this configuration with linearly polarised light is found, namely Equation (43). Identical asymptotic behaviour follows.
- (ii)
In this scenario, Equation (42) results in
yielding modulated inverse
R behaviour in the far zone and
R−3 dependence in the near zone.
For the second term of Equation (14), substituting Equation (32), taking
N + 1 ~
N and carrying out the molecular rotational averages with relations (26) and (27), produces
Equation (49) is found to be identical to Equation (42). Hence, identical polarisation analysis results and limiting forms follow, with expressions (43)–(46) and (48) all doubling up.
4. Field-Induced Potential between an Electric Dipole-Polarisable Molecule and an Electric Octupole-Polarisable Molecule
The final radiation-induced energy shift to be calculated that depends on electric octupole coupling and involves the scattering of a real photon at different centres is that between an electric dipole-polarisable molecule,
A, and its pure octupole analogue,
B. This potential is of the same order of magnitude as that considered in the previous section, and therefore should be included for consistency. The electric dipole induced at
A is given by Equation (4), with pure electric dipole polarisability (5). By dint of the polarisability characteristic of
B, only an electric octupole moment is induced there,
with the electric octupole polarisability defined as
The two moments couple via the tensor
written in Equation (13), giving rise to an interaction energy of the form
Substituting relation (18), carrying out the polarisation sum using Equation (19) and performing the pair average via the result (22), and rewriting the coupling tensor using Equation (13), the energy shift expression (72) becomes
To obtain the isotropic potential energy, use is made of the result (26) for the orientational average of the electric dipole polarisability, while for pure octupole polarisability, we employ the result [
23,
33]
Noting that [
23]
and
the energy shift expression (73) becomes
Evaluating the product of geometric tensors finally results in
Inspection of expression (78) shows that this contribution to the optical binding energy depends on both octupole weight-1 and weight-3 terms. As in the previous examples considered, the weight-1-dependent term may be viewed as a higher-order correction to the leading pure electric dipole contribution, the coefficients preceding each term in square brackets being identical. The second term of Equation (78) is proportional to
, which is a third-rank tensor with 27 independent components. The coefficients within square brackets of this term are identical to that found in the dispersion potential between an electric dipole-polarisable molecule and an octupole-polarisable one [
23].
Polarisation Analysis
Substituting relation (18) into Equation (72) and performing the molecular orientational averages yields
- (a)
Linearly polarised beam.
- (i)
: Inserting the geometric tensors produces
with asymptotic limits
and
with the last of these limits independent of
- (ii)
: Now,
with limiting forms
and
- (b)
Circularly polarised radiation.
- (i)
: For this arrangement, given by Equation (80), and with identical asymptotic limits (81) and (82).
- (ii)
: This results in
with far- and near-zone asymptotes
and
The last expression is again independent of the octupole weight-1 term.
5. Octupole-Dependent Contributions to Optical Binding Arising from the Static Mechanism
The three examples presented thus far involve the scattering of the real photon at different molecular sites. This is frequently referred to as the dynamic mechanism. We now consider the situation in which the incident laser light is emitted and absorbed by the same particle, either
A or
B. This is often termed the static mechanism. As well as the static electric dipole moment, higher-order response tensors participate in the induction of multipole moments. In the case of the induced electric dipole moment, this is given by
where
is the static electric dipole moment,
is the pure electric dipole polarisability and given earlier by Equation (5),
is the mixed electric dipole–octupole polarisability tensor defined by Equation (8),
is the pure electric dipole first hyperpolarisability tensor and
is the electric dipole–dipole–octupole first hyperpolarisability tensor of molecule
ξ. Similarly, for the induced octupole moment,
where
is the pure electric octupole polarisability (71) and
is the pure electric octupole first hyperpolarisability tensor.
Let
A be purely electric dipole-hyperpolarisable and
B be described by a permanent electric octupole moment. From Equations (89) and (90), the interaction energy coupling the electric dipole induced at
A and the octupole induced at
B can be written as
On account of
B being polar, with its ground state permanent octupole moment appearing to leading order, the static coupling tensor
features in Equation (91). A corresponding term arises on interchanging
A and
B, so that
with the factor of ½ due to index symmetry, allowing for
A and
B to be identical or distinct.
Taking the expectation value of Equation (92) over the ground state of both particles and the radiation field with state
photons as before, noting that
on taking
N >> 1, the first term of expression (92) is
The molecular factors are averaged using the results [
33] for a third-rank Cartesian tensor
and
yielding the isotropic potential
This vanishes because the laser polarisation factor is symmetric in
j,
k while the Levi-Civita tensor is antisymmetric in this pair. The second term of Equation (92) similarly vanishes. Hence,
A second situation occurs when
A is electric dipole–dipole–octupole-first-hyperpolarisable and
B is characterised by a permanent electric dipole moment. From Equations (89) and (90), coupling of an induced electric octupole at
A with an induced electric dipole at
B, and between induced electric dipoles at each centre, leads to
Examining the first term of Equation (99), substituting the relation (93) produces
Now,
on taking
N + 1 ~
N, so that the second term of Equation (99) is
Both Equations (100) and (102) vanish on orientational averaging because [
34] (p. 66)
To this order of approximation, no term dependent on a single octupole moment contributes to the static optical binding mechanism.
6. Conclusions
While the dominant contribution to the change in energy shift between two interacting neutral polarisable molecules induced by an intense radiation field is that proportional to the electric dipole polarisablility of each species, interesting effects arise when higher multipole moment terms are accounted for. When magnetic dipole coupling terms are included, for instance, the contribution to the potential that is dependent on both transition electric and magnetic dipole moments of either or both species is discriminatory, changing sign when one enantiomer is switched to its isomer of opposite handedness. In this work, new terms contributing to the optical binding energy that are dependent upon one or two electric octupole moments at each centre were calculated. These include the field-induced energy shift between an electric dipole-polarisable molecule and a second that is either mixed electric dipole–octupole-polarisable or purely electric octupole-polarisable, and between two mixed electric dipole–octupole-polarisable molecules.
Results were obtained for both particles randomly oriented relative to each other, and for all directions of propagation of the incident beam relative to the inter-nuclear displacement vector. The first and third examples mentioned above were found to depend only on the octupole weight-1 moment, which has three independent components and the properties of a vector. They may be viewed as higher-order corrections to the leading electric dipole–dipole energy shift. The second system’s potential depended on both octupole weight-1 and weight-3 terms, with the former also interpreted as a higher-order correction. All three of these examples fell under the dynamic mechanism classification of the process as the laser beam is absorbed at one site and emitted at the other. Also examined were two examples that belong to the static mechanism categorisation. In this case, the beam of light is absorbed and emitted at the same centre. This necessitates one molecule—the one emitting or absorbing only a virtual photon—being polar, and possessing a static multipole moment. In one case, this was a permanent electric dipole moment so that the other particle was characterised by its mixed electric dipole–dipole–octupole first hyperpolarisability. In the second case, the polar molecule was described by its ground state permanent octupole moment while the other molecule possessed purely electric dipole first hyperpolarisability. Both energy shifts vanished after orientational averaging. For the surviving isotropic potentials, a complete polarisation analysis was also performed. Linearly and circularly polarised incident radiation propagating parallel and perpendicular to the inter-nuclear axis were four scenarios considered that could stimulate possible experimental detection.
Instead of perturbation theory, which is prohibitive when considering higher-order multipole moments due to the excessive number of diagrams that need to be drawn and summed to yield the matrix element or energy shift for a process, especially one like optical binding that is intrinsically higher-order in its number of light–matter couplings, an alternative computational method was employed. This involved first calculating the multipole moments induced at each centre by the applied electromagnetic field based on the polarisability characteristics of the responder. These moments then couple via their respective retarded multipole–multipole interaction tensor. An expectation value taken over the ground electronic state of each molecule with the radiation field in the state containing N photons yielded the field-induced energy shift. The respective dispersion potential readily follows on letting N = 0.