1. Introduction
Casimir’s discovery [
1] that two parallel uncharged ideal-metal planes attract each other marked the beginning of active research and led to new concepts and far-reaching implications in both fundamental and applied physics (see, e.g., the monographs [
2,
3,
4,
5]). The Casimir force results from zero-point and thermal (if the planes are kept at nonzero temperature) fluctuations of the electromagnetic field. On ideal-metal planes the tangential component of electric field and the normal component of magnetic induction vanish leading to the discrete oscillation frequencies. A finite difference between the infinite free energies in the presence and in the absence of the planes determines the finite Casimir free energy per unit area. The derivative of this free energy with respect to the separation between the planes with a minus sign is the Casimir force.
The Lifshitz theory generalized the Casimir force for two thick material plates (semispaces) in thermal equilibrium with the environment [
6,
7,
8]. In this case, the electromagnetic response of the plate materials is given by the standard continuity boundary conditions formulated in terms of the frequency-dependent dielectric permittivity and (for magnetic plates) magnetic permeability. The Lifshitz theory treats the van der Waals force as a particular case of the Casimir force at separations of a few nanometers where one can neglect the effects of relativistic retardation. During the last two decades the Lifshitz theory was generalized for material bodies of arbitrary shape [
9,
10,
11].
The interest in measuring the Casimir force was rekindled by the experiment of [
12] (an unaccounted systematic error which was present in this experiment was corrected by its author in [
13]). The beginning of the 21st century was marked by performing several high-precision experiments on measuring the Casimir force between metallic test bodies. These experiments were made by means of a micromechanical torsional oscillator [
14,
15,
16,
17,
18,
19] and an atomic force microscope [
20,
21,
22,
23,
24,
25,
26] at the separation distances exceeding 150 nm. A comparison between the measurement data and theoretical predictions of the Lifshitz theory showed that these predictions are experimentally excluded if the electromagnetic response of metals (Au, Ni) at low frequencies is described by the well-established Drude model which takes the proper account of the relaxation properties of conduction electrons [
14,
15,
16,
17,
18,
19,
20,
21,
22,
23,
24,
25,
26]. If, however, the response of metals at low frequencies is described by the dissipationless plasma model—a very good agreement between the measurement data and theoretical predictions was revealed [
14,
15,
16,
17,
18,
19,
20,
21,
22,
23,
24,
25,
26].
For distances below 100 nm, a comparison between the Lifshitz theory and the measurement data demonstrated good agreement (see, e.g., [
27,
28,
29,
30]). Analysis made it apparent that the problem is related to the low-frequency contribution to the Casimir force, which is determined by conduction electrons [
31,
32,
33]. At separations below 100 nm the major contribution to the Casimir force is given by the core electrons whereas with increasing separation the role of conduction electrons progressively increases. Thus, it is very surprising that the Drude model does not work at large separations because it provides a well-established description of collisions of conduction electrons which are responsible for the resistance of metals (Ohm’s law).
No less surprise is the fact that, for metals with perfect crystal lattices described by the Drude model, the Casimir entropy goes with vanishing temperature to the negative quantity which depends on the parameters of a system, such as the separation between the plates and the plasma frequency, i.e., the Nernst heat theorem is violated [
34,
35,
36,
37]. If the experimentally consistent plasma model is used, the Nernst heat theorem is satisfied [
34,
35,
36,
37] in spite of the fact that this model disregards the dissipation of conduction electrons, which plays an important role at low frequencies. To satisfy the Nernst heat theorem with the Drude model, it was suggested to take into account a nonzero residual relaxation at zero temperature, which exists for crystal lattices with defects [
38,
39,
40]. It was noted, however, that the perfect crystal lattice is an equilibrium system with the nondegenerate ground state so that the Nernst heat theorem must be satisfied in this case [
41]. Moreover, in equilibrium there must be no defects [
42]. The unresolved problems of disagreement between the theoretical predictions of the Lifshitz theory and measurement data, on the one hand, and the requirements of thermodynamics, on the other hand, were named the Casimir puzzle [
43,
44,
45].
During the last 20 years, many attempts to resolve the Casimir puzzle have been undertaken based on the assumptions that there might be some unaccounted systematic effects in the performed experiments or inaccuracies in the comparison between the measurement data and theoretical predictions (see [
4,
46,
47,
48] for review). These efforts, although they did not resolve the Casimir puzzle, have had a profound impact on future investigations of the problem. Specifically, it was underlined that, in the frequency region of the anomalous skin effect, the concept of the dielectric permittivity depending only on frequency used in the standard Lifshitz theory loses its meaning and one has to take into account the effects of spatial dispersion [
42,
49,
50,
51,
52,
53,
54,
55]. In doing so, the reflection coefficients entering the Lifshitz formula for the Casimir pressure are expressed via the nonlocal surface impedances [
49,
52]. This approach allowed some progress in resolution of thermodynamic part of the Casimir puzzle [
53,
54] but did not remove the contradiction between experiment and theory because for metallic test bodies corrections to the Casimir pressure due to spatial nonlocality in the region of the anomalous skin effect turned out to be too small [
51,
52].
To find a solution for both components of the Casimir puzzle, recently the phenomenological spatially nonlocal dielectric permittivities have been proposed [
56,
57,
58], which almost coincide with the permittivity of the standard Drude model for the propagating waves satisfying the condition
where
is the magnitude of the wave vector projection on the plane of Casimir plates and
is the frequency. The suggested permittivities, however, can deviate from the Drude model significantly for the evanescent waves which obey the inequality
and, thus, are characterized by a pure imaginary
. Allowing for only real
, one can say that the evanescent waves are off the mass shell in the free space. It was shown that the Lifshitz theory using the permittivities of this kind is in agreement with the Nernst heat theorem [
59] and with the measurement data of all high-precision experiments on measuring the Casimir force [
56,
57,
58]. Thus, we can say that the nonlocal permittivities, proposed for metals in [
56,
57,
58], are to some extent analogous to the permittivities of graphene, which are also spatially nonlocal but were derived on the basis of first principles of quantum electrodynamics at nonzero temperature using the formalism of polarization tensor [
60,
61]. It was shown that, for graphene described by these permittivities, there is no Casimir puzzle, i.e., the Lifshitz theory is consistent with both the experimental results [
62,
63] and the requirements of thermodynamics [
64,
65].
The proposed phenomenological permittivities [
56,
57,
58] take into account the dissipation of conduction electrons in metals and leave almost unchanged the transverse magnetic (
p) reflection coefficient as compared to the standard Drude model. As to the transverse electric (
s) reflection coefficient, its value at zero frequency becomes nonzero for both propagating and evanescent waves as it holds if the plasma model is used (recall that for evanescent waves the transverse electric reflection coefficient calculated using the Drude model vanishes at zero frequency). Keeping in mind that the suggested nonlocal permittivities are of entirely phenomenological character, it is desirable to determine their possible role in some physical phenomena other than the Casimir effect.
The direct measurement of the
s-type reflection coefficient in the range of evanescent waves cannot be performed because all conventional methods (for example, ellipsometry) are designed for propagating waves. Nevertheless, the evanescent fields are actively discussed in the literature in relation to the possibility to overcome the diffraction limit in optics [
66]. Although the physics of surface plasmon polaritons [
67] provides a considerable body of data related to the region of rather large
, it is restricted to only the transverse magnetic (
p) reflection coefficient. The widely used techniques of total internal reflection and frustrated total internal reflection give the possibility to test the reflection properties in the region of
only slightly exceeding
. This is connected to the fact that there are no transparent media with sufficiently large index of refraction (see the long-performed experiments [
68,
69,
70]). The method of nano frustrated total internal reflection exploits the illumination of a nanoparticle of size
placed near the material surface under investigation. This approach is employed in the near-field optical microscopy, which allows exceeding of the standard resolution limit using the evanescent waves [
66,
71]. The near-field optical microscopes, however, are mostly employed in various technological applications [
71] and are more sensitive to the
p-polarized electromagnetic field [
72].
In this paper, we suggest the
Experimentum crucis, which can reliably select between different response functions of metals to the low-frequency electromagnetic field in the area of evanescent waves. To do so, we consider an oscillating magnetic dipole spaced at a distance of a few millimeters from a metallic plate. The electric field of this dipole is negligibly small as compared to the magnetic one. We obtain explicit expressions for the latter via the transverse electric (
s) reflection coefficient of the plate without specifying its form. According to our results, the lateral field components are fully determined by the contribution of extremely evanescent waves
, whereas the contribution of propagating waves is negligibly small. The components of the magnetic field are calculated as functions of lateral separation from the dipole for different oscillation frequencies using the dielectric permittivities of metal described by the Drude model, plasma model, and the spatially nonlocal model of [
56,
57,
58,
59]. It is shown that the lateral field components, which are fully determined by reflections on the plate, strongly depend on the model of used dielectric permittivity. For example, the lateral field component computed using the Drude model is smaller by up to a factor of
than that computed using the plasma model or the spatially nonlocal model of [
56,
57,
58,
59]. Thus, by measuring the field of oscillating magnetic dipole spaced near a metallic plate, one can probe the response function of metal in the range of evanescent waves.
The above results are used to propose the
Experimentum crucis for the Casimir puzzle. For this purpose, we consider the Casimir pressure at large separations between the plates. It is demonstrated that the differences in theoretical predictions obtained using the Drude model, plasma model and the spatially nonlocal model of [
56,
57,
58,
59] are fully determined by different contributions from the evanescent waves to the transverse electric (
s) reflection coefficient. Possible parameters of the experiment, which could determine the response function of metals to the extremely evanescent waves, are suggested. The role of this experiment not only for the Casimir physics, but for a wide range of topics in optics and condensed matter physics is discussed.
The brief presentation of the part of obtained results related to the lateral component of the field of magnetic dipole is published in the form of a letter in [
73].
The paper is organized as follows. In
Section 2, it is shown that the conflict between experiment and theory in the Casimir physics is determined by a contribution from the evanescent waves.
Section 3 contains the derivation of the field of oscillating magnetic dipole spaced near a metallic plate. In
Section 4, we compute the dipole field for different frequencies and find how the contribution from the evanescent waves depends on the response function of a metal.
Section 5 suggests the
Experimentum crucis for the Casimir puzzle. In
Section 6 and
Section 7, the reader will find a discussion and our conclusions.
Appendix A and
Appendix B contain some necessary details about the mathematical derivations.
2. Role of Evanescent Waves in the Conflict between Experiment and Theory in Casimir Physics
It has long been known that the Casimir pressure between two thick material plates (semispaces) maintained at separation
a at temperature
T in thermal equilibrium with the environment is given by the Lifshitz formula which can be written in two mathematically equivalent representations. The representation that has enjoyed the widest application is in terms of the pure imaginary discrete Matsubara frequencies,
with
[
4,
5,
6,
7,
8]
where the magnitude of the lateral wave vector component
was already defined in
Section 1,
, the prime on the summation sign in
l divides the term with
by 2,
is the Boltzmann constant, and the sum in
is over the transverse magnetic (TM or
p) and transverse electric (TE or
s) polarizations of the electromagnetic field. The respective reflection coefficients are given by:
where
is the dielectric permittivity of the plates calculated at the Matsubara frequencies and
. The signs of the coefficients (
2) are chosen in such a way that in the limiting case of ideal metal planes
and
.
Another representation of the Lifshitz formula is in terms or real frequencies [
4,
5,
6,
7]. It expresses the Casimir pressure as the sum of contributions from the propagating and evanescent waves,
where
and
Here,
,
and the reflection coefficients (
2) with
can be transformed into the familiar Fresnel coefficients [
74]
where
Note that
and
. Because of this, Equation (
4), expressing the contribution of propagating waves in the Casimir pressure, contains the rapidly oscillating function
and is not convenient for numerical computations.The contribution of evanescent waves (
5) does not have this disadvantage (different aspects of evanescent waves with respect to the Casimir force are considered in [
31,
32,
75,
76,
77]).
The problems arising in the Casimir physics become more pronounced in the limit of large separations (high temperatures) satisfying the condition
. In fact, at room temperature
K this limit occurs at separations exceeding
m. In the case of large separations, (
1) and (
4), (
5) become much simplified. Thus, in (
1) the total result is given by the zero-frequency term with
along, whereas all terms with
are exponentially small
The condition of large separations can be identically rewritten in the form
, where
is the characteristic frequency giving the major contribution to the Casimir pressure. Because of this, in the limiting case of large separations one can use the inequality
and replace
in (
4) and (
5) with the first expansion term in powers of the small parameter
Then, (
4) and (
5) take the form:
and
The quantity (
8) can be easily calculated using the dielectric permittivities of the Drude and plasma models,
where
is the plasma frequency and
is the relaxation parameter. Substituting (
12) in (
2), one obtains the values of the reflection coefficients at zero Matsubara frequency for the Drude and plasma models:
and
respectively.
Then, substituting (
13) and (
14) in (
8), one finds the following different results [
4]:
where
is the Riemann zeta function. Note that differences between the theoretical predictions for the thermal Casimir force using the Drude and plasma models were discussed long ago in [
78] and [
79], respectively.
Now we consider what the relative role of the propagating and evanescent waves in the Casimir pressures is (
15). Let us start with the pressures calculated using the Drude model. For the TE polarization of the electromagnetic field, this question was solved in [
33,
80] using (
10) and (
11). It was shown that
Thus, the contribution of evanescent waves leads to the repulsive Casimir force and cancels the contribution of propagating waves with the result
. As for the TM mode, numerical computations performed for the evanescent waves using (
11) result in the attractive Casimir pressure whose magnitude is negligibly small as compared to
in (
16). Then, using the first equality in (
15) for the total pressure, one can conclude that:
We are coming now to the Casimir pressure calculated using the plasma model. It is given by the second equality in (
15). The contributions of evanescent waves (
11) to this result for both the TM and TE modes are equal to zero because the dielectric permittivity
in (
12) and the reflection coefficients in (
6) are the real-valued functions. In regard to the propagating waves, the TM and TE modes contribute to (
10) equally because at large separations the characteristic frequency
goes to zero and the TM and TE reflection coefficients (
16) go to 1 and −1, respectively, if the plasma model is used in calculations.
As a result, taking into account the second equality in (
15), one obtains:
and
From (
16)–(
19), we conclude that, at large separations between metallic plates, both the TM and TE modes make equal contributions to the Casimir pressures determined by the propagating waves irrespective of whether the Drude or the plasma model is used in calculations. As is seen from (
17) and (
19), the TM mode also makes equal (zeroth) contributions to the Casimir pressures computed using either the Drude or the plasma model.
Thus, according to (
16) and (
18), the heart of the problem is in the TE contribution to the Casimir pressure determined by the low-frequency evanescent waves. This contribution essentially depends on the model of dielectric permittivity used. However, as discussed in
Section 1, the required experimental data concerning the dielectric response of metal to the evanescent waves are few.
Before proceeding to the proposed experiment which could supply us with these data, we briefly dwell on the phenomenological spatially nonlocal dielectric permittivity suggested in [
56,
57,
58,
59]. This permittivity combines the merits of the Drude model by taking into account the dissipation of conduction electrons and of the plasma model by bringing the Lifshitz theory in agreement with experiments on measuring the Casimir force and with the requirements of thermodynamics.
In the spatially nonlocal case it is convenient to express the reflection coefficients in terms of the surface impedances
. Thus, the TE reflection coefficient of our interest is given by [
49,
81]:
where in the approximation of specular reflection the TE impedance is expressed via the transverse dielectric permittivity
as [
49,
82]:
Here,
is the wave vector component perpendicular to the plate and
. Recall that in the presence of spatial dispersion the response of metals to perpendicular and parallel to
k electric fields is described by the transverse—
—and longitudinal—
—dielectric permittivities [
74,
81]. In so doing, the TM impedance and respective TM reflection coefficient are expressed via both of them.
In [
58], the following transverse dielectric permittivity was suggested to describe the response of metals to the low-frequency electromagnetic field:
where
, the value of
, and
is the Fermi velocity. This permittivity brings theoretical predictions of the Lifshitz theory in good agreement with the experimental data of all high-precision measurements of the Casimir force. In so doing, the theoretical predictions are scarcely affected by the form of
.
The distinctive feature of the permittivity (
22) is that for the propagating waves satisfying a condition
it differs little from the standard Drude model in (
12). This is because for the propagating waves the magnitude of an addition to unity in the parentheses of (
22) becomes negligibly small:
As to extremely evanescent waves for which
, the permittivity (
22) may significantly deviate from the Drude model.
It should be noted that the permittivity (
22) is of entirely phenomenological character and does not lay claim that it is generally applicable. Specifically, the imaginary part of
is positive under the condition
. For good metals, such as Au or Cu, this condition is satisfied for
, i.e., at separations exceeding a few tens of nanometers (we recall that the high-precision experiments of [
14,
15,
16,
17,
18,
19,
20,
21,
22,
23,
24,
25,
26] were performed at much larger separations).
Below, we consider the experimental test for the response of metals in the area of evanescent waves unrelated to the Casimir effect which may give additional information regarding the validity of dielectric functions (
12) and (
22) in this area.
3. Field of Oscillating Magnetic Dipole near a Metallic Plate
We consider the field of an oscillating magnetic dipole situated above a thick metallic plate and show that it bears a close analogy to the Casimir pressure at large separation between two metallic plates discussed in the previous section.
Let the magnetic moment of our dipole be directed along the
z axis which is perpendicular to the surface of metallic plate situated in the
plane (see
Figure 1):
The dipole with the oscillation frequency is spaced in a vacuum at height h above the plate. It is assumed that h is much larger than the dipole size.
Let us consider first the electromagnetic field of the magnetic dipole (
24) in free space, i.e., in the absence of a conducting plate. Under a condition that the wavelength
is much larger than the dipole size the solution of this problem is contained in [
83] (§ 72, problem 1). All fields considered below depend on
t as
. The monochromatic components of the magnetic field can be represented in the form:
where the radius-vector
is directed from the magnetic dipole to the observation point,
,
, and the Gaussian system of units is used. In a similar way, the monochromatic components of the electric field are given by
Below, we consider the quasistationary case of low oscillation frequencies
and large
. In this case, the magnitudes of components of the electric field (
26) are smaller than those of the magnetic field (
25) by the factor of:
Taking into account that in the configuration of suggested experiment
(see
Section 5), one can safely neglect by the electric field as compared to the magnetic one.
Below we also use the Fourier expansion of the field of magnetic dipole (
25) in plane waves
where
and
,
It should be recalled that the third component of the wave vector,
, is not independent but is expressed via
and
by the dispersion relation
, see the first equality in (
7), following from the Maxwell equations.
The explicit expressions for the components of the Fourier transform
are (see the
Appendix A where these expressions are proven):
Now we find the field of the magnetic dipole in the presence of a thick metallic plate (see
Figure 1). We do not specify the explicit form of the reflection coefficients
to take into account both cases of spatially local and nonlocal response of the plate to the electromagnetic field.
In fact in our case the reflected wave is fully determined by only the TE reflection coefficient . As was shown above, the electric field of our magnetic dipole is negligibly small in comparison with the magnetic one. At the same time, the magnetic field of the TM-polarized wave is perpendicular to the plane of incidence, i.e., parallel to the surface of metallic plate. Thus, it is not reflected from this plate.
The magnetic field of magnetic dipole in a vacuum in the domain above the metallic plate can be calculated most simply by the method of images [
84,
85,
86]. As shown in the
Appendix B, the same result can be obtained using the method of the Green tensor. To apply the method of images, we use the coordinate system illustrated in
Figure 1 where the magnetic dipole is positioned at the point
,
and the fictitious (image) dipole is in the point
,
. The magnetic moment of the image dipole is in opposition to the real one and its magnitude depends on the reflectivity properties of the metallic plate. In the domain above the plate the resulting magnetic field is found as a superposition of the fields of real and image dipoles. Using (
29), we find the Fourier transform of the resulting magnetic field in the presence of a plate:
where the first terms on the right-hand side are related to the image dipole and the second to the real dipole.
The presence of only TE reflection coefficient in (
30) is explained by the fact that the TE polarized field has a nonzero component of magnetic field perpendicular to the surface. Note that for the oscillating electric dipole the magnetic field is negligibly small, as compared to the electric one, and, vice versa, the TM polarized field has a nonzero component of the electric field perpendicular to the surface.
Now we substitute (
30) into (
28) and represent the total magnetic field above the plate in the form:
Here,
is the monochromatic field (
25) of real magnetic dipole where
z should be replaced with
and
is the monochromatic field of the image dipole given by:
The integrals in (
32) can be calculated as is described in
Appendix A for the case of magnetic dipole in the absence of a plate. Thus, using the polar coordinates defined in (
A2) and (
A3), the
x-component,
, is rewritten as:
Substituting here (
A5), one finds:
Likewise, the
y-component in (
32),
, is expressed as:
Using the polar coordinates, for the
z-component in (
32) we obtain:
Then, applying (
A12), we rewrite (
36) as:
Needless to say that both the propagating (
) and evanescent (
) waves contribute to (
34), (
35), and (
37). However, in the region
under consideration here one can neglect by the phase in the exponential factors under the integrals. Then, the contribution of propagating waves in the reflected field, as compared to the contribution of evanescent ones, is suppressed by the factor
. In the configuration of suggested experiment
and
(see
Section 5). This means that the contribution of propagating waves to (
34), (
35), and (
37) does not play any role and one can replace the lower integration limit in these equations with
. The results of numerical computations confirm this conclusion.
Finally, we write out the full components of the monochromatic field (
31) of the magnetic dipole in the presence of metallic plate in a form convenient for computations:
The expression for
is similar to (
38) where
x is replaced with
y.
In fact, (
34), (
35) and (
37) have much in common with (
10) and (
11) presenting the contributions of propagating and evanescent waves to the Casimir pressure between two parallel metallic plates at large separations. At first glance (
10) and (
11) are quite different from the components of magnetic field (
38) and (
39) because the former are presented as the integrals over real frequencies whereas the latter are measured as the function of dipole oscillation frequency
. Moreover, (
10) and (
11) contain the second power of the factor
and take into account multiple reflections on two surfaces whereas (
38) and (
39) describe the single reflection on one surface.
What is in common, however, is the role of the TE-polarized evanescent waves in both cases. The contribution of evanescent waves to the Casimir force at large separations is given by (
11). In the case of magnetic field of an oscillating magnetic dipole located above metallic plate, the contribution of evanescent waves is given by the first terms in (
38) and (
39). By using the dimensionless variables
in (
11) and
in (
38) and (
39), one can see that the values of the Casimir force and of the dipole magnetic field are determined by one and the same integration domain
. The lateral components of the field of magnetic dipole in the presence of metallic plate computed at
are determined by the TE-polarized evanescent waves alone. If we replace the lower integration limit
in (
38) and (
39) with 0, this does not influence the obtained results, i.e., the propagating waves do not contribute to them. According to
Section 2, the contribution of propagating waves (
10) to the Casimir pressure is common for both polarizations and for any model of dielectric permittivity used. The TM-polarized evanescent waves do not contribute to the Casimir pressure irrespective of the used dielectric function. Thus, again, only the contribution of TE-polarized evanescent waves eventually determines the value of the total pressure which is different for different models of the dielectric permittivity. This means that by measuring magnetic field of the magnetic dipole with fixed oscillation frequency spaced near a metallic plate one can test the response of metals to the TE-polarized evanescent waves in the wide region of real frequencies contributing to the magnetic field and simultaneously to the Casimir pressure. The obtained results can be used to resolve the Casimir puzzle.
4. Influence of the Response Function in the Evanescent Domain on the Field of Magnetic Dipole near Metallic Plate
Here, we perform numerical computations of the dipole field (
38) and (
39) for different oscillator frequencies and different models of the dielectric permittivity. Taking into account that an impact of metal on the dipole field is performed through the reflected field in (
31), it is advantageous to make all computations at the dipole height
. At this height, the lateral component (
38) of the dipole field coincides with that of the reflected field whereas the relative role of reflected field in
z-component (
39) is maximized.
As a metal of the plate, we consider Cu with the Drude parameters
rad/s and
rad/s [
87]. The oscillation frequencies of the magnetic dipole should be chosen so as to gain sufficient reflectance and simultaneously large enough deviation between theoretical predictions obtained using the Drude and plasma models.
To estimate the corresponding parameters, one can use the spatially local reflection coefficients (
6). We introduce a new variable,
and rewrite the TE reflection coefficient defined in (
6) as:
where
and
.
The major contribution to
computed in the region of small frequencies under consideration here using the dielectric permittivity of the Drude model is given by its imaginary part,
If
, the reflections on a metallic surface occur much as in the case when the plasma model is used. If
is much less than unity, only a small fraction of the magnetic field is reflected from the surface. Thus, for our purposes, it would be best to choose
For the Drude parameters of Cu indicated above and the typical value of
cm (see
Section 5), one obtains
rad/s. Based on this result, the computations below are performed for the oscillation frequencies of magnetic dipole equal to 0.2, 2, 10, 20, and 100 rad/s.
In
Figure 2, we present the computational results for
normalized to
as a function of separation from the magnetic dipole along the
x-axis. In this case, one should put in (
38)
and
. We recall that all computations are performed at
cm. Thus, only the reflected field contributes to the result, i.e.,
=
. The solid lines labeled 1, 2, 3, 4, and 5 are computed using the Drude model and the Fresnel reflection coefficient
defined in (
6) for the dipole oscillation frequencies equal to 0.2, 2, 10, 20, and 100 rad/s, respectively. The top solid line is computed using the plasma model and the coefficient
defined in (
6). Under the condition (
27) it does not depend on the value of
. The top dashed line overlapping with the solid one is computed using the spatially nonlocal dielectric permittivity (
22) and the impedance reflection coefficient (
20). The Fermi velocity
m/s was determined under an assumption of the spherical Fermi surface. The top dashed line is also independent on the value of
as well as in the case of the top solid line plotted for the plasma model.
From
Figure 2, one can see that, as expected, the difference between the values of
computed using the Drude and plasma models increases with decreasing oscillation frequency of the magnetic dipole. This difference reaches several orders of magnitude and can be observed experimentally (see
Section 5).
Similar results for the imaginary part of
=
normalized to
are shown in
Figure 3. Here, again, the lines labeled 1, 2, 3, 4, and 5 are computed by (
38) using the Drude model for the dipole oscillation frequencies of 0.2, 2, 10, 20, and 100 rad/s, respectively. For the plasma model, the reflection coefficient
in (
38) is real so that the obtained result is
. If the spatially nonlocal dielectric permittivity defined in (
22) is used in computations, the obtained values of
are by ten orders of magnitudes smaller than those shown by the line labeled 1 in
Figure 3.
By comparing
Figure 3 with
Figure 2, it is seen that
is much greater than
for each value of oscillation frequency if the Drude model is used in computations. If, however, the plasma or spatially nonlocal models are used, the magnitude of the real part,
, is far in excess of
which is either equal to zero or negligibly small.
For completeness, we also consider the component
given in (
39) which contains both contributions of the source dipole field and its part reflected from the plate. The computational results for
at
cm normalized to
as the function of separation from the magnetic dipole along the
x-axis are presented in
Figure 4. The solid line is obtained for
rad/s using the Drude model and another solid line overlapping with the short-dashed one, which does not depend on frequency, is found by means of the plasma model and the spatially nonlocal model (
22). For comparison purposes, the long-dashed line shows the quantity
computed in the absence of metallic plate for a magnetic dipole in free space. It is also frequency-independent for the frequency region under consideration.
From
Figure 4, it is seen that, at separations exceeding 2 cm, there is a considerable deviation between the theoretical predictions of the Drude model, on the one hand, and the plasma and our spatially nonlocal models, on the other hand. This deviation takes the largest value at
cm. It is seen also that at large separations the impact of metallic plate on theoretical predictions obtained using the Drude model is far less than can be obtained when the plasma or spatially non-local models are used. Because of this, the measurement of the
z-component of magnetic field can be considered as complementary to the measurement of the lateral component (according to
Figure 2 for
rad/s the largest deviation between different theoretical predictions for the lateral component occurs at
cm).
In the inset to
Figure 4, we also plot the ratio of
to
as the function of separation where the first quantity characterizes the reflected field and the second one is for the dipole field in free space in the absence of metallic plate. The solid line is computed using the Drude model and another solid line overlapping with the dashed one is found by means of the plasma and spatially nonlocal models. According to the inset, the quantity
determined by the reflected field changes its sign with separation if it is computed using the plasma model and the spatially nonlocal model (
22) but preserves its sign if the Drude model is used in computations. Some computational results for
are presented in
Section 5.
It is also interesting to illustrate the dependence on separation of the total field magnitude:
In
Figure 5, the quantity
normalized to
is shown as the function of separation from the magnetic dipole by the solid lines labeled 1 and 2, which are computed using the Drude model for the dipole oscillation frequencies of 0.2 and 100 rad/s, respectively. Another solid line overlapping with the dashed one is computed employing the plasma and the spatially nonlocal models. These lines do not depend on the value of oscillation frequency.
From
Figure 5, it is seen that at separations of 2.67 and 4.85 cm the theoretical predictions of all three models coincide for the dipole oscillation frequency equal to 100 and 0.2 rad/s, respectively. At these separations, there are, however, significant deviations between the theoretical predictions of the Drude mode, on the one hand, and the plasma and spatially nonlocal models, on the other hand, if the oscillation frequencies of 0.2 and 100 rad/s are used, respectively.
In the inset to
Figure 5, we plot as a function of separation the ratio of the reflected,
to the total (
45) magnitudes of the field. The solid lines 1 and 2 are computed using the Drude model for the dipole oscillation frequencies of 0.2 and 100 rad/s, respectively, whereas another solid line overlapping with the dashed one using the plasma and the spatially nonlocal models. As is seen from the inset, for the oscillation frequency of 0.2 rad/s the quantity
is separation-independent. The largest difference between the theoretical predictions is reached at 5 cm for the dipole oscillation frequency of 100 rad/s.
The above computations are for a plate made of copper. The obtained results are, however, applicable to any metallic material with appropriately chosen frequencies of the magnetic dipole. Thus, for the plate made of B-doped Si in the metallic state used in the experiment of [
88], it holds:
rad/s and
rad/s. Then, using (
44), one obtains
rad/s. Because of this, by choosing the dipole oscillation frequencies in the same proportion to
as above to
, i.e.,
,
,
,
, and
, one arrives to the same computational results as were obtained for a copper plate. This means that
Figure 2,
Figure 3,
Figure 4 and
Figure 5 plotted for the case of Cu plate are also valid for a B-doped Si plate if the oscillation frequencies of the magnetic dipole 0.2, 2, 10, 20, and 100 rad/s are replaced with
,
,
,
,
rad/s, respectively.
5. Experimentum crucis for the Casimir Puzzle
The oscillating magnetic dipole with a magnetic moment (
24) is generated by an alternating current
through a circular turn or a small coil. In this case,
where
R is the turn radius and
N is the number of turns. In accordance to
Figure 1, the loop should be located in the plane
with a center at the point
. The surface of the metallic plate is in the plane
. As shown in
Section 3, if the magnetic field is measured at the points with coordinates
, i.e., along the line parallel to the
x-axis at height
h, its lateral component is completely determined by reflections from the surface of a metal [see (
38)].
For a qualitative estimation of the coil parameters, we can consider the case of total reflection,
, which is close to metal plate described by the plasma model. Then the Fourier transform of the
x-component of the magnetic field is obtained from (
30) where the second term does not contribute because
. As a result, under the condition (
27), the
x-component of the field is given by the last term on the right-hand side of (
25) where one should replace
z with
,
here,
.
The scale factor
in (
48) determines the field magnitude. Thus, in order to have the greater field, one has to use larger
, i.e., the larger current, turn radius, and the number of turns, but smaller distance
h between the magnetic dipole and the metallic plate. These conditions are somewhat contradictory. The calculations made in
Section 3 assume that both
h and the distance
x between the source and the observation point are large in comparison with the dipole size. This means that, in order to decrease
h, one should decrease the dipole size accordingly. For too-small coils, however, it would be difficult to ensure large currents and the number of turns.
In the literature, one can find examples of small coils manufactured by the methods of mini- [
89] and micro- [
90,
91] technologies (micro-electromagnets for atom manipulation have long been created [
92]). The coils of such kind contain of about 10 turns with a radius of 1 mm and have height of about 1 mm. These coils are able to support the current up to
statA (this is equal to 1 A in SI). The respective magnetic dipole moment computed by (
47) is equal to:
For estimation of the produced lateral field, we put
mm so that
. Then (
48) leads to:
By using the value of
from (
49), one can see that the obtained field is in a very good agreement with the top line in
Figure 2 at
cm which was computed using the plasma model and the spatially nonlocal model of [
58].
Thus, the numerical data in
Figure 2,
Figure 3,
Figure 4 and
Figure 5 in combination with the value of
in (
49) can be used for the calculation of the magnetic field produced by the oscillating magnetic dipole at different separations and oscillation frequencies depending on the used model of electromagnetic response of metal to the low-frequency evanescent waves. As an example, if the Drude model is used and the dipole oscillation frequency is equal to 100 rad/s (see the line 5 in
Figure 2), one has:
i.e., by more than an order of magnitude smaller field than for the plasma model in the separation region from 1 to 2 cm. In this case, the parameter
defined in (
27) is of the order of
. The Drude model predicts all the more damping of the lateral magnetic field for smaller oscillation frequencies (see the lines 1–4 in
Figure 2).
The obtained values of the magnetic field are quite measurable. Using different methods and laboratory techniques, the presently available limit for the magnetic field resolution is down to
T [
93,
94,
95]. This resolution is already reached [
96,
97] by means of traditional search coil magnetometers consisting of several micro induction coils [
98,
99]. The same sensitivity is provided by the SQUID magnetometers [
100,
101,
102]. Some novel types of magnetometers using the magnetoresistance and magnetoelectric effects, as well as the optically pumped atomic magnetometers, possess an even higher resolution when measuring alternating magnetic fields (see [
94,
103] for a review). In the proposed experiment one should measure the magnetic fields of the order of tens of nT which are by several orders of magnitude larger than the resolution limit of currently available magnetic sensors. Therefore, the systematic effects, i.e., from the industrial or thermal noise, although deserve careful attention, should not be a crucial problem. Based on this, measuring the magnetic field of an oscillating magnetic dipole can be considered the
Experimentum crucis, which can provide valuable information concerning the electromagnetic response of metals to the evanescent waves of low frequency. Depending on the magnitude of measured magnetic field, one will have the possibility to clearly discriminate between the description of metal in the range of evanescent waves by means of the Drude model, as the first option, or by means of the plasma model or spatially nonlocal model as the second option.
Direct demonstration of an inadequacy of the Drude model in the range of low-frequency evanescent waves by means of an experiment performed in the area of classical electrodynamic phenomena would be of great help to resolution of the Casimir puzzle discussed in
Section 1 and
Section 2. The Casimir effect is the quantum phenomenon dealing with the concept of a quantum vacuum. As indicated above, the high precision experiments of [
14,
15,
16,
17,
18,
19,
20,
21,
22,
23,
24,
25,
26], measuring the Casimir force, have already demonstrated that the theoretical predictions of the Lifshitz theory obtained using the Drude model at separations exceeding 150 nm are excluded by the measurement data. An independent confirmation of a conclusion that the Drude model does not work in the area of evanescent waves would also be of prime importance for several divisions of optics and condensed matter physics dealing with evanescent waves.
Keeping in mind subsequent experiments, in addition to the computational results presented in
Section 4, here we report the real and imaginary parts of
and
computed as the functions of dipole frequency for the above parameters of the experimental configuration, [i.e., for the value of
defined in (
49),
cm,
] with different response functions of a metal and for several values of the separation
x between the dipole and the point of observation. The computations are performed using (
38) and (
39) for a copper plate.
In
Figure 6a, we plot the computational results for
=
as a function of dipole oscillation frequency. The solid lines labeled 1, 2, and 3 are computed by using the Drude model at separations from the magnetic dipole equal to
x = 10, 20, and 30 mm, respectively. The computational results obtained using the plasma model and the spatially nonlocal model (
22) at 10, 20, and 30 mm are shown by the three solid lines and overlapping with them dashed lines counted from the top to bottom, respectively. In
Figure 6b, similar computational results obtained using the Drude model for
=
are presented with the same notations. If the plasma model is used in computations, one arrives to
. As to the spatially nonlocal dielectric permittivity (
22), it gives
that is smaller by the factor of
comparing to those shown in
Figure 6b.
As one can see from
Figure 6a, different theoretical predictions for
can differ significantly. For
mm and
rad/s the field computed using the plasma model is 280 times larger than that for the Drude model. For the Drude model, the magnitude of
increases with the frequency but stays constant for the plasma model. Although the difference between predictions of these models is reduced, for
rad/s the predictions still differ by a factor of 14 which is sufficient for the experimental discrimination. As to
Figure 6b, the values of
shown by any line computed using the Drude model are much outside the theoretical predictions based on the plasma and spatially nonlocal models where
vanishes or nearly vanishes.
In
Figure 7a, using the same values of all experimental parameters, we plot the computational results for
as a function of the oscillation frequency of the magnetic dipole. The results obtained using the Drude and plasma (or the spatially nonlocal) models are shown by the three pairs of solid and solid-dashed lines, respectively. These pairs are labeled as 1, 2, and 3 for separations from the magnetic dipole equal to 30, 40, and 50 mm, respectively. The choice of separations, which is different from that in
Figure 6, is caused by the fact that at larger separations one obtains greater differences between the
z-components of magnetic field computed using different models of the dielectric permittivity (see
Figure 4). The computational results for
obtained using the Drude model are shown in
Figure 7b with the same notations. In common with above, the plasma model leads to
and the magnitudes of
computed using the spatially nonlocal model (
22) are much below those shown in
Figure 7b.
As is seen in
Figure 7a, the theoretical predictions for
obtained using the Drude and plasma models are not so much different as for
. This is caused by a contribution to the
z-component of the dipole field in the absence of metallic plate. Thus, measuring the lateral component of the magnetic field is preferable for a clear discrimination between different models of the dielectric permittivity. It should be noted, however, that the imaginary part of the field can be determined by simultaneous measurements of the field magnitude and the phase shift. In doing so, if some nonzero phase shift relative to the magnetic moment is registered, this would be in favor of the Drude model. From
Figure 7b, one can also see that
, computed using the Drude model, possesses an interesting property by changing its sign from positive to negative at the dipole oscillation frequencies equal to 29.4, 58.9, and 162.8 rad/s when the measurements are performed at separations from the dipole of 50, 40, and 30 mm, respectively.
6. Discussion
In the foregoing, we have demonstrated that the longstanding problem of a disagreement between the theoretical predictions of fundamental Lifshitz theory and numerous high-precision experiments on measuring the Casimir force at separations exceeding 150 nm is determined by the role of TE-polarized low-frequency evanescent waves. The analysis of literature sources shows that the electromagnetic response of metals to the evanescent waves of this type has not yet been sufficiently explored experimentally. Taking into account that the predictions of the Lifshitz theory come into conflict with the measurement data when the low-frequency electromagnetic response of metals is described by the well-established Drude model taking into account the relaxation properties of conduction electrons, it is very desirable to independently verify whether or not this model is applicable in the range of low-frequency evanescent waves.
In this paper, we have derived an expression for the electromagnetic field of an oscillating magnetic dipole spaced above a thick metallic plate. According to our results obtained under the reasonable experimental conditions, the lateral components of magnetic field of this dipole are entirely determined by the TE-polarized evanescent waves, whereas the electric field turns out to be negligibly small. The derivation was performed both by the method of images and using the Green tensor of the boundary problem in the general case of a spatially nonlocal electromagnetic response of the plate metal with the coinciding results.
The real and imaginary parts of the components of magnetic field above the copper plate were computed by describing the electromagnetic response of copper using the dissipative Drude model, dissipationless plasma model and the spatially nonlocal phenomenological dielectric permittivity of [
58] which takes into account the relaxation properties of conduction electrons and simultaneously brings the theoretical predictions of the Lifshitz theory in agreement with all high-precision experiments on measuring the Casimir force. It was shown that for the typical parameters of the magnetic dipole designed in the form of 1-mm coil the lateral components of magnetic field computed using the Drude model, as the first option, and the plasma model or the spatially nonlocal model of [
58], as the second option, differ by up to several orders of magnitude depending on the oscillation frequency of the dipole. In doing so, the magnitudes of the field component remain quite measurable using the available laboratory techniques. Based on this, we have proposed an experiment which can independently either validate or disprove an applicability of the Drude model in the range of low-frequency TE-polarized evanescent waves.
The suggested experiment would elucidate the roots of the Casimir puzzle by making clear whether a disagreement between the theoretical predictions and the measurement data is caused by an inadequacy of the Drude model in the range of evanescent waves or by some other problems. Because of this, it could be called the
Experimentum crucis. Based on the studies of the Casimir force in graphene systems, one could argue in favor of the inadequacy of the Drude model. The point is that the spatially nonlocal response of graphene was derived on the basis of first principles of quantum electrodynamics via the formalism of the polarization tensor [
60,
61]. As mentioned in
Section 1, using the respective nonlocal dielectric functions, the theoretical predictions of the Lifshitz theory were shown to be in perfect agreement with the requirements of thermodynamics [
64,
65] and with experiments on measuring the Casimir force in graphene systems [
62,
63].