Next Article in Journal
Longitudinal Variations in Equatorial Ionospheric TEC from GPS, Global Ionosphere Map and International Reference Ionosphere-2016 during the Descending and Minimum Phases of Solar Cycle 24
Previous Article in Journal
Rotational Energy Extraction from the Kerr Black Hole’s Mimickers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Experimentum crucis for Electromagnetic Response of Metals to Evanescent Waves and the Casimir Puzzle

by
Galina L. Klimchitskaya
1,2,
Vladimir M. Mostepanenko
1,2,3,* and
Vitaly B. Svetovoy
4
1
Central Astronomical Observatory at Pulkovo of the Russian Academy of Sciences, 196140 Saint Petersburg, Russia
2
Peter the Great Saint Petersburg Polytechnic University, 195251 Saint Petersburg, Russia
3
Kazan Federal University, 420008 Kazan, Russia
4
Frumkin Institute of Physical Chemistry and Electrochemistry, Russian Academy of Sciences, Leninsky Prospect 31 bld. 4, 119071 Moscow, Russia
*
Author to whom correspondence should be addressed.
Universe 2022, 8(11), 574; https://doi.org/10.3390/universe8110574
Submission received: 2 October 2022 / Revised: 26 October 2022 / Accepted: 27 October 2022 / Published: 31 October 2022
(This article belongs to the Section Field Theory)

Abstract

:
It is well known that the Casimir force calculated at large separations using the Lifshitz theory differs by a factor of 2 for metals described by the Drude or plasma models. We argue that this difference is entirely determined by the contribution of transverse electric (s) evanescent waves. Taking into account that there is a lack of experimental information on the electromagnetic response of metals to low-frequency evanescent waves, we propose an experiment on measuring the magnetic field of an oscillating magnetic dipole spaced in a vacuum above a thick metallic plate. According to our results, the lateral components of this field are governed by the transverse electric evanescent waves and may vary by orders of magnitude depending on the model describing the permittivity of the plates used in calculations and the oscillation frequency of the magnetic dipole. Measuring the lateral component of the magnetic field for typical parameters of the magnetic dipole designed in the form of a 1-mm coil, one could either validate or disprove the applicability of the Drude model as a response function of metal in the range of low-frequency evanescent waves. This will elucidate the roots of the Casimir puzzle lying in the fact that the theoretical predictions of the Lifshitz theory using the Drude model are in contradiction with the high-precision measurements of the Casimir force at separations exceeding 150 nm. Possible implications of the suggested experiment for a wide range of topics in optics and condensed matter physics dealing with evanescent waves are discussed.

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 k < ω / c where k 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 k > ω / c and, thus, are characterized by a pure imaginary k z . Allowing for only real k z , 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 k , 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 k only slightly exceeding ω / c . 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 R 1 / k 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 ( k ω / c ) , 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 10 4 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, ω l = i ξ l = 2 π i k B T l / with l = 0 , 1 , 2 , [4,5,6,7,8]
P ( a , T ) = k B T π l = 0 0 d k q l k α r α 2 ( i ξ l , k ) e 2 a q l 1 r α 2 ( i ξ l , k ) e 2 a q l ,
where the magnitude of the lateral wave vector component k = ( k x 2 + k y 2 ) 1 / 2 was already defined in Section 1, q l = ( k 2 + ξ l 2 / c 2 ) 1 / 2 , the prime on the summation sign in l divides the term with l = 0 by 2, k B 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:
r TM ( i ξ l , k ) = ε l q l p l ε l q l + p l , r TE ( i ξ l , k ) = q l p l q l + p l ,
where ε l ε ( i ξ l ) is the dielectric permittivity of the plates calculated at the Matsubara frequencies and p l = p l ( k ) = ( k 2 + ε l ξ l 2 / c 2 ) 1 / 2 . The signs of the coefficients (2) are chosen in such a way that in the limiting case of ideal metal planes r TM = 1 and r TE = 1 .
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,
P ( a , T ) = P prop ( a , T ) + P evan ( a , T ) ,
where
P prop ( a , T ) = 2 π 2 0 d ω coth ω 2 k B T 0 k 0 d k k α Im q r α 2 ( ω , k ) e 2 a q 1 r α 2 ( ω , k ) e 2 a q
and
P evan ( a , T ) = 2 π 2 0 d ω coth ω 2 k B T k 0 d k k q α Im r α 2 ( ω , k ) e 2 a q 1 r α 2 ( ω , k ) e 2 a q .
Here, q = q ( ω , k ) = ( k 2 k 0 2 ) 1 / 2 , k 0 = ω / c and the reflection coefficients (2) with p = p ( ω , k ) = [ k 2 ε ( ω ) k 0 2 ] 1 / 2 can be transformed into the familiar Fresnel coefficients [74]
r TM ( ω , k ) = ε ( ω ) q p ε ( ω ) q + p = ε ( ω ) q ˜ p ˜ ε ( ω ) q ˜ + p ˜ , r TE ( ω , k ) = q p q + p = q ˜ p ˜ q ˜ + p ˜ ,
where
q ˜ = q ˜ ( ω , k ) = ( k 0 2 k 2 ) 1 / 2 , p ˜ = p ˜ ( ω , k ) = [ ε ( ω ) k 0 2 k 2 ] 1 / 2 .
Note that q ˜ = i q and p ˜ = i p . Because of this, Equation (4), expressing the contribution of propagating waves in the Casimir pressure, contains the rapidly oscillating function exp ( 2 i a q ˜ ) 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 T c / ( 2 a k B ) . In fact, at room temperature T = 300 K this limit occurs at separations exceeding 6 μ 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 l = 0 along, whereas all terms with l 1 are exponentially small
P ( a , T ) = k B T 2 π 0 d k k 2 α r α 2 ( 0 , k ) e 2 a k 1 r α 2 ( 0 , k ) e 2 a k .
The condition of large separations can be identically rewritten in the form ω c 2 k B T , where ω c = c / ( 2 a ) 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 ω 2 k B T and replace coth ( ω / 2 k B T ) in (4) and (5) with the first expansion term in powers of the small parameter
coth ω 2 k B T 2 k B T ω .
Then, (4) and (5) take the form:
P prop ( a , T ) = k B T π 2 0 d ω ω 0 k 0 d k k q ˜ α Re r α 2 ( ω , k ) e 2 i a q ˜ 1 r α 2 ( ω , k ) e 2 i a q ˜ ,
and
P evan ( a , T ) = k B T π 2 0 d ω ω k 0 d k k q α Im r α 2 ( ω , k ) e 2 a q 1 r α 2 ( ω , k ) e 2 a q .
The quantity (8) can be easily calculated using the dielectric permittivities of the Drude and plasma models,
ε D ( ω ) = 1 ω p 2 ω ( ω + i γ ) , ε p l ( ω ) = 1 ω p 2 ω 2 ,
where ω p 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:
r D , TM ( 0 , k ) = 1 , r D , TE ( 0 , k ) = 0
and
r p l , TM ( 0 , k ) = 1 , r p l , TE ( 0 , k ) = c k c 2 k 2 + ω p 2 c k + c 2 k 2 + ω p 2 ,
respectively.
Then, substituting (13) and (14) in (8), one finds the following different results [4]:
P D ( a , T ) = k B T 8 π a 3 ζ ( 3 ) , P p l ( a , T ) = k B T 4 π a 3 ζ ( 3 ) ,
where ζ ( z ) 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
P D , TE prop ( a , T ) = k B T 8 π a 3 ζ ( 3 ) , P D , TE evan ( a , T ) = k B T 8 π a 3 ζ ( 3 ) .
Thus, the contribution of evanescent waves leads to the repulsive Casimir force and cancels the contribution of propagating waves with the result P D , TE = 0 . 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 P D , TE evan in (16). Then, using the first equality in (15) for the total pressure, one can conclude that:
P D , TM prop ( a , T ) = k B T 8 π a 3 ζ ( 3 ) , P D , TM evan ( a , T ) = 0 .
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 ε p l ( ω ) 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 ω c 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:
P p l , TE prop ( a , T ) = k B T 8 π a 3 ζ ( 3 ) , P p l , TE evan ( a , T ) = 0
and
P p l , TM prop ( a , T ) = k B T 8 π a 3 ζ ( 3 ) , P p l , TM evan ( a , T ) = 0 .
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 Z α ( ω , k ) . Thus, the TE reflection coefficient of our interest is given by [49,81]:
r TE ( ω , k ) = q Z TE ( ω , k ) + i k 0 q Z TE ( ω , k ) i k 0 ,
where in the approximation of specular reflection the TE impedance is expressed via the transverse dielectric permittivity ε Tr ( ω , k ) as [49,82]:
Z TE ( ω , k ) = i k 0 π d k z k 0 2 ε Tr ( ω , k ) k 2 .
Here, k z is the wave vector component perpendicular to the plate and k 2 = k 2 + k z 2 . 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— ε Tr ( ω , k ) —and longitudinal— ε L ( ω , k ) —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:
ε Tr ( ω , k ) = 1 ω p 2 ω ( ω + i γ ) 1 + i v Tr k ω ,
where k = | k | , the value of v Tr = 1.5 v F , and v F 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 ε L ( ω , k ) .
The distinctive feature of the permittivity (22) is that for the propagating waves satisfying a condition k < k 0 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:
v Tr k ω = v Tr c k k 0 = v Tr c 1 .
As to extremely evanescent waves for which k k 0 , 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 ε Tr is positive under the condition k < γ / v Tr . For good metals, such as Au or Cu, this condition is satisfied for a 1 / ( 2 k ) > v Tr / ( 2 γ ) , 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 ( x , y ) plane (see Figure 1):
m = ( 0 , 0 , m 0 e i ω d t ) .
The dipole with the oscillation frequency ω d 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 λ d = 2 π c / ω d 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 exp ( i ω d t ) . The monochromatic components of the magnetic field can be represented in the form:
H x ( ω d , r ) = m 0 x z r 2 k d 2 r + 3 i k d r 2 3 r 3 e i k d r , H y ( ω d , r ) = m 0 y z r 2 k d 2 r + 3 i k d r 2 3 r 3 e i k d r , H z ( ω d , r ) = m 0 k d 2 r + i k d r 2 1 r 3 z 2 r 2 k d 2 r + 3 i k d r 2 3 r 3 e i k d r ,
where the radius-vector r = ( x , y , z ) is directed from the magnetic dipole to the observation point, r = | r | , k d = ω d / c , and the Gaussian system of units is used. In a similar way, the monochromatic components of the electric field are given by
E ( ω d , r ) = i m 0 k d i k d r 2 1 r 3 e i k d r y x 0 .
Below, we consider the quasistationary case of low oscillation frequencies ω d and large λ d r . 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:
k d r = ω d r c = 2 π r λ d 1 .
Taking into account that in the configuration of suggested experiment k d r 10 9 (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 exp ( i k r ) where k = ( k x , k y ) and r = ( x , y ) ,
H ( ω d , r ) = 1 ( 2 π ) 2 d k e i k r H ( ω d , k , z ) .
It should be recalled that the third component of the wave vector, k z , is not independent but is expressed via k x and k y by the dispersion relation k z 2 q ˜ d 2 = k d 2 k 2 , see the first equality in (7), following from the Maxwell equations.
The explicit expressions for the components of the Fourier transform H ( ω d , k , z ) are (see the Appendix A where these expressions are proven):
H x ( y ) ( ω d , k , z ) = 2 π i m 0 k x ( y ) sign ( z ) e i q ˜ d | z | , H z ( ω d , k , z ) = 2 π i m 0 k 2 q ˜ d e i q ˜ d | z | .
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 r α ( ω d , k ) 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 α = TE . 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 x = y = 0 , z = h and the fictitious (image) dipole is in the point x = y = 0 , z = h . 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:
H x ( y ) ( p ) ( ω d , k , z ) = 2 π i m 0 k x ( y ) r TE ( ω d , k ) e i q ˜ d ( z + h ) + sign ( z h ) e i q ˜ d | z h | , H z ( p ) ( ω d , k , z ) = 2 π i m 0 k 2 q ˜ d r TE ( ω d , k ) e i q ˜ d ( z + h ) + e i q ˜ d | z h | ,
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:
H ( p ) ( ω d , r ) = H ( R ) ( ω d , r ) + H ( ω d , r ) .
Here, H ( ω d , r ) is the monochromatic field (25) of real magnetic dipole where z should be replaced with z h and H ( R ) ( ω d , r ) is the monochromatic field of the image dipole given by:
H x ( y ) ( R ) ( ω d , r ) = i 2 π m 0 d k r TE ( ω d , k ) k x ( y ) e i [ k r + q ˜ d ( z + h ) ] , H z ( R ) ( ω d , r ) = i 2 π m 0 d k r TE ( ω d , k ) k 2 q ˜ d e i [ k r + q ˜ d ( z + h ) ] .
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, H x ( R ) , is rewritten as:
H x ( R ) ( ω d , r ) = i 2 π m 0 0 d k k 2 r TE ( ω d , k ) e i q ˜ d ( z + h ) 0 2 π d ψ cos ψ e i k ρ cos ( ψ φ ) .
Substituting here (A5), one finds:
H x ( R ) ( ω d , r ) = m 0 x ρ 0 d k k 2 J 1 ( k ρ ) r TE ( ω d , k ) e i q ˜ d ( z + h ) .
Likewise, the y-component in (32), H y ( R ) , is expressed as:
H y ( R ) ( ω d , r ) = m 0 y ρ 0 d k k 2 J 1 ( k ρ ) r TE ( ω d , k ) e i q ˜ d ( z + h ) .
Using the polar coordinates, for the z-component in (32) we obtain:
H z ( R ) ( ω d , r ) = i m 0 2 π 0 d k k 3 q ˜ d r TE ( ω d , k ) e i q ˜ d ( z + h ) 0 2 π d ψ e i k ρ cos ( ψ φ ) .
Then, applying (A12), we rewrite (36) as:
H z ( R ) ( ω d , r ) = i m 0 0 d k k 3 J 0 ( k ρ ) q ˜ d r TE ( ω d , k ) e i q ˜ d ( z + h ) .
Needless to say that both the propagating ( k < k d ) and evanescent ( k > k d ) waves contribute to (34), (35), and (37). However, in the region r λ d 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 ( k d h ) 3 = ( 2 π h / λ d ) 3 . In the configuration of suggested experiment h r and ( k d h ) 3 10 27 (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 k d . 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:
H x ( p ) ( ω d , r ) = m 0 x ρ k d d k k 2 J 1 ( k ρ ) r TE ( ω d , k ) e k 2 k d 2 ( z + h ) m 0 x ( z h ) r 2 k d 2 r + 3 i k d r 2 3 r 3 e i k d r ,
H z ( p ) ( ω d , r ) = m 0 k d d k k 3 k 2 k d 2 J 0 ( k ρ ) r TE ( ω d , k ) e k 2 k d 2 ( z + h ) + m 0 k d 2 r + i k d r 2 1 r 3 ( z h ) 2 r 2 k d 2 r + 3 i k d r 2 3 r 3 e i k d r .
The expression for H y ( p ) ( ω d , r ) 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 ω d . Moreover, (10) and (11) contain the second power of the factor r α e a q 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 v = 2 a k 2 k 0 2 in (11) and v = ( z + h ) k 2 k d 2 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 v 1 . The lateral components of the field of magnetic dipole in the presence of metallic plate computed at z = h are determined by the TE-polarized evanescent waves alone. If we replace the lower integration limit k d 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 z = h . 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 ω p = 1.12 × 10 16 rad/s and γ p = 1.38 × 10 13 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,
w = h q d = h k 2 k d 2 ,
and rewrite the TE reflection coefficient defined in (6) as:
r TE ( ω d , k ) = w w 2 K ( ω d ) w + w 2 K ( ω d ) ,
where
K ( ω d ) = [ ε ( ω d ) 1 ] ω d 2 ω h 2
and ω h = c / h .
The major contribution to | K ( ω d ) | computed in the region of small frequencies under consideration here using the dielectric permittivity of the Drude model is given by its imaginary part,
| K ( ω d ) | = γ ω p 2 ω d ( ω d 2 + γ 2 ) ω h 2 ω p 2 ω d γ ω h 2 .
If | K ( ω d ) | 1 , the reflections on a metallic surface occur much as in the case when the plasma model is used. If | K ( ω d ) | 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
| K ( ω d ) | ω p 2 ω d γ ω h 2 1 , ω d Ω γ ω h 2 ω p 2 .
For the Drude parameters of Cu indicated above and the typical value of h = 1 cm (see Section 5), one obtains Ω Cu 100 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 | Re H x ( p ) | normalized to m 0 as a function of separation from the magnetic dipole along the x-axis. In this case, one should put in (38) r = ( x , 0 , h ) and ρ = x . We recall that all computations are performed at z = h = 1 cm. Thus, only the reflected field contributes to the result, i.e., H x ( p ) = H x ( R ) . The solid lines labeled 1, 2, 3, 4, and 5 are computed using the Drude model and the Fresnel reflection coefficient r TE 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 r TE 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 v F Cu = 1.6 × 10 6 m/s was determined under an assumption of the spherical Fermi surface. The top dashed line is also independent on the value of ω d 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 | Re H x ( p ) | 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 H x ( p ) = H x ( R ) normalized to m 0 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 r TE in (38) is real so that the obtained result is Im H x ( p ) / m 0 = 0 . If the spatially nonlocal dielectric permittivity defined in (22) is used in computations, the obtained values of Im H x ( p ) / m 0 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 Im H x ( p ) is much greater than | Re H x ( p ) | 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, | Re H x ( p ) | , is far in excess of Im H x ( p ) which is either equal to zero or negligibly small.
For completeness, we also consider the component H z ( p ) given in (39) which contains both contributions of the source dipole field and its part reflected from the plate. The computational results for | Re H z ( p ) | at z = h = 1 cm normalized to m 0 as the function of separation from the magnetic dipole along the x-axis are presented in Figure 4. The solid line is obtained for ω d = 100 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 | Re H z | / m 0 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 x = 5 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 ω d = 100 rad/s the largest deviation between different theoretical predictions for the lateral component occurs at x = 1 cm).
In the inset to Figure 4, we also plot the ratio of Re H z ( R ) to Re H z 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 Re H z ( R ) 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 Im H z are presented in Section 5.
It is also interesting to illustrate the dependence on separation of the total field magnitude:
H ˜ ( ω d , r ) = | H x ( p ) ( ω d , r ) | 2 + | H z ( p ) ( ω d , r ) | 2 .
In Figure 5, the quantity H ˜ ( ω , r ) normalized to m 0 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,
H ˜ H ( R ) ( ω d , r ) = | H x ( R ) ( ω d , r ) | 2 + | H z ( R ) ( ω d , r ) | 2 ,
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 H ˜ H ( R ) / H ˜ 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: ω p Si = 7.0 × 10 14 rad/s and γ p Si = 1.5 × 10 14 rad/s. Then, using (44), one obtains Ω Si = 2.75 × 10 5 rad/s. Because of this, by choosing the dipole oscillation frequencies in the same proportion to Ω Si as above to Ω Cu , i.e., Ω Si / 500 , Ω Si / 50 , Ω Si / 10 , Ω Si / 5 , and Ω Si , 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 5.5 × 10 2 , 5.5 × 10 3 , 2.75 × 10 4 , 5.5 × 10 4 , 2.75 × 10 5 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 I 0 exp ( i ω d t ) through a circular turn or a small coil. In this case,
m 0 = 1 c π N I 0 R 2 ,
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 z = h with a center at the point ( 0 , 0 , h ) . The surface of the metallic plate is in the plane z = 0 . As shown in Section 3, if the magnetic field is measured at the points with coordinates ( x , 0 , h ) , 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, r TE = 1 , 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 z = h . 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 z + h = 2 h ,
H x ( p ) ( x ) = 6 m 0 x h ( x 2 + 4 h 2 ) 5 / 2 = m 0 h 3 6 x ˜ ( 4 + x ˜ 2 ) 5 / 2 .
here, x ˜ = x / h .
The scale factor m 0 / h 3 in (48) determines the field magnitude. Thus, in order to have the greater field, one has to use larger m 0 , 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 I 0 = 3 × 10 9 statA (this is equal to 1 A in SI). The respective magnetic dipole moment computed by (47) is equal to:
m 0 = 3.14 × 10 2 e r g / O e = 3.14 × 10 5 A m 2 .
For estimation of the produced lateral field, we put x = h = 10 mm so that x ˜ = 1 . Then (48) leads to:
H x ( p ) = 3.36 × 10 3 O e = 3.36 × 10 7 T = 0.27 A / m .
By using the value of m 0 from (49), one can see that the obtained field is in a very good agreement with the top line in Figure 2 at x = 1 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 m 0 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:
| Re H x ( p ) | 3.36 × 10 4 O e = 3.36 × 10 8 T = 0.027 A / m ,
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 k 0 r defined in (27) is of the order of 10 9 . 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 10 13 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 H x ( p ) and H z ( p ) computed as the functions of dipole frequency for the above parameters of the experimental configuration, [i.e., for the value of m 0 defined in (49), z = h = 1 cm, y = 0 ] 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 | Re H x ( p ) | = | Re H x ( R ) | 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 | Im H x ( p ) | = | Im H x ( R ) | are presented with the same notations. If the plasma model is used in computations, one arrives to Im H x ( p ) = 0 . As to the spatially nonlocal dielectric permittivity (22), it gives Im H x ( p ) that is smaller by the factor of 10 8 comparing to those shown in Figure 6b.
As one can see from Figure 6a, different theoretical predictions for | Re H x ( p ) | can differ significantly. For x = 10 mm and ω d = 10 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 Re H x ( p ) increases with the frequency but stays constant for the plasma model. Although the difference between predictions of these models is reduced, for ω d = 100 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 Im H x ( p ) shown by any line computed using the Drude model are much outside the theoretical predictions based on the plasma and spatially nonlocal models where Im H x ( p ) vanishes or nearly vanishes.
In Figure 7a, using the same values of all experimental parameters, we plot the computational results for | Re H z ( p ) | 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 | Im H z ( p ) | obtained using the Drude model are shown in Figure 7b with the same notations. In common with above, the plasma model leads to Im H z ( p ) = 0 and the magnitudes of Im H z ( p ) 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 | Re H z ( p ) | obtained using the Drude and plasma models are not so much different as for | Re H x ( p ) | . 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 Im H z ( p ) , 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].

7. Conclusions

According to our results, the suggested experiment will be capable of clearly discriminating between the magnetic fields computed using the Drude and plasma models but not between the fields computed using the plasma model and the spatially nonlocal permittivity of [58]. The latter, in any case, does not claim a complete description of the electromagnetic response of metals but can be considered as an example that, through the spatial nonlocality, in the framework of the Lifshitz theory it is possible to reach an agreement with the requirements of thermodynamics and the measurement data of high-precision experiments having a regard for the relaxation properties of conduction electrons (see also [104], which considers other nonlocal permittivities with the same aim).
To conclude, the fundamental theoretical description of the response of metals to both the propagating and evanescent waves could be derived in the future using the methods of quantum field theory [105], as was already done for graphene. The first steps in this direction were made in [106]. In the meantime, the results of the proposed Experimentum crucis may not only shed light on the Casimir puzzle, but also find applications in a number of topics in optics and condensed matter physics which deal with evanescent waves including studies on surface plasmon polaritons, near-field optical microscopy, total internal reflection and many others.

Author Contributions

Conceptualization, G.L.K., V.M.M. and V.B.S.; Investigation, G.L.K., V.M.M. and V.B.S.; Writing—original draft, V.M.M.; Writing—review & editing, G.L.K. and V.B.S. All authors have read and agreed to the published version of the manuscript.

Funding

G.L.K. and V.M.M. were partially supported by the Peter the Great Saint Petersburg Polytechnic University in the framework of the Russian state assignment for basic research (Project No. FSEG-2020-0024). The work of V.M.M. was also supported by the Kazan Federal University Strategic Academic Leadership Program. V.B.S. was partially supported by the Russian Science Foundation, grant No. 20-19-00214.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Fourier Transform of the Field of Magnetic Dipole

Here, we present the proof of (29). For this purpose we substitute the first equality in (29) into the first component of the Fourier transform (28) and obtain
H x ( ω d , r ) = i 2 π sign ( z ) d k k x e i k r e i q ˜ d | z | .
Introducing the polar coordinates
k = ( k x , k y ) = ( k cos ψ , k sin ψ ) , r = ( x , y ) = ( ρ cos φ , ρ sin φ ) ,
one finds
k r = k x x + k y y = k ρ cos ( ψ φ ) .
Then, using (A3), one can rewrite (A1) in the form
H x ( ω d , r ) = i 2 π m 0 sign ( z ) 0 d k k 2 e i q ˜ d | z | 0 2 π d ψ cos ψ e i k ρ cos ( ψ φ ) .
The integral with respect to ψ in (A4) can be found using the integral representation for the Bessel functions J n ( x ) (see [107])
0 2 π d ψ cos ψ e i k ρ cos ( ψ φ ) = 2 π i J 1 ( k ρ ) cos φ .
Then (A4) takes the form
H x ( ω d , r ) = m 0 sign ( z ) cos φ 0 d k k 2 J 1 ( k ρ ) e i q ˜ d | z | .
The integral in (A6) is calculated by using the following result which is derived from [108] (Equation 2.12.10.10) by differentiation with respect to the parameter p:
0 d k k ν + 1 e | z | k 2 + τ 2 J ν ( k ρ ) = 2 π ρ ν | z | τ 3 2 + ν ( z 2 + ρ 2 ) ν 2 3 4 K ν + 3 2 ( τ z 2 + ρ 2 ) ,
where K ν ( x ) is the Bessel function of imaginary argument. Note that (A7) does not coincide with the expression 2.12.20.6 given in [108] just for this integral. The point is that the expression 2.12.20.6 is in contradiction with the correct expression 2.12.20.20 and contains two typos: the power of τ on the right-hand side, which is indicated as 3/2, must be replaced with 3 / 2 + ν and the power of ( z 2 + ρ 2 ) , which is shown as ν 3 / 4 , must be replaced with ν / 2 3 / 4 .
Our integral in (A6) is obtained from (A7) as a particular case ν = 1 and τ = i k 0
0 d k k 2 J 1 ( k ρ ) e | z | k 2 k d 2 = 2 π ρ | z | i k d 5 2 ( z 2 + ρ 2 ) 5 4 K 5 / 2 ( i k d z 2 + ρ 2 ) .
Notice that although (A7) was derived under a condition Re τ > 0 , it can be analytically continued to the point Re τ = 0 leading to (A8).
The Bessel function of imaginary argument in (A8) can be presented in terms of the elementary functions [108]
K 5 / 2 ( x ) = π 2 x 2 + 3 x + 3 x 5 / 2 e x .
Substituting (A9) taken for x = i k 0 r (we recall that z 2 + ρ 2 = r 2 ) in (A8), one obtains
0 d k k 2 J 1 ( k ρ ) e | z | k 2 k d 2 = ρ | z | 3 3 i k d r k d 2 r 2 r 5 e i k d r .
Finally, after the substitution of (A10) in (A6), we arrive to the sought for first equality in (25). The second expression in (29) for the Fourier component H y ( ω d , k , z ) is proven in perfect analogy to the above with a replacement of k x with k y .
Now we deal with the third expression in (29) for the Fourier component H z ( ω d , k , z ) . Substituting it in the z component of (28) and using the polar coordinates (A2) and (A3) one obtains
H z ( ω d , r ) = i m 0 2 π 0 d k k 3 q ˜ d e i q ˜ d | z | 0 2 π d ψ e i k ρ cos ( ψ φ ) .
Using the result 2.5.27.19. in [109]
0 2 π d ψ e i k ρ cos ( ψ φ ) = 2 π J 0 ( k ρ ) ,
Equation (A11) can be rewritten in the form
H z ( ω d , r ) = m 0 0 d k k 3 J 0 ( k ρ ) k 2 k d 2 e | z | k 2 k d 2 .
With the help of equality [108]
J 0 ( k ρ ) = 2 k ρ J 1 ( k ρ ) J 2 ( k ρ )
Equation (A13) can be transformed to
H z ( ω d , r ) = m 0 2 ρ 0 d k k 2 J 1 ( k ρ ) k 2 k d 2 e | z | k 2 k d 2 0 d k k 3 J 2 ( k ρ ) k 2 k d 2 e | z | k 2 k d 2 .
Both integrals in (A15) are calculated using the result 2.12.10.10 in [108] with ν = 1 and ν = 2 which is analytically continued to the point Re τ = Re ( i k d ) = 0 . Then (A15) is rewritten as
H z ( ω d , r ) = m 0 2 π i k d 3 / 2 ( z 2 + ρ 2 ) 3 / 4 2 K 3 / 2 ( i k d z 2 + ρ 2 ) i ρ 2 k d ( z 2 + ρ 2 ) 1 / 2 K 5 / 2 ( i k d z 2 + ρ 2 ) .
Taking into account (A9) and similar identity [108]
K 3 / 2 ( x ) = π 2 x + 1 x 3 / 2 e x
at x = i k d r = i k d z 2 + ρ 2 , one arrives from (A16) to the third expression in (25).
This concludes the proof of (29).

Appendix B. Field Derivation by the Green Function Method

Lifshitz and co-workers [6,7,8] have used the Green function method to deduce the famous formula for the Casimir force. In this Appendix, we use this method to find the field of a magnetic dipole above a conductive plane derived in the main text [see (30)] by the method of images. We do this in the general case of a spatially nonlocal response of the plane material. In accordance with the main text, we are interested only in TE (s) polarization.
In the Green function method, an elementary current is the source of the field. To reproduce the field of magnetic dipole, this elementary current can be chosen as (along with the main text, the Gaussian system of units is used)
J ( r ) = c m × δ ( r ) .
By applying the Biot-Savart law, it can be easily checked that this current reproduces the field of magnetic dipole in empty space [see (25) and (26)].
We are looking for the electromagnetic field above the conductive plane at z = 0 when the source is located in the point r 0 = ( 0 , 0 , h ) . If the magnetic moment m = m 0 e i ω d t oscillates with the frequency ω d , then all the fields oscillate with the same frequency and the corresponding monochromatic components obey the Maxwell equations
× E ( ω d , r ) = i k d H ( ω d , r ) ,
× H ( ω d , r ) = 4 π c J ( ω d , r ) i k d E ( ω d , r ) ,
where k d = ω d / c is the absolute value of the wave vector. In vacuum the magnetic fields B and H coincide in the Gaussian system. The electric field is related to the Green tensor G i j ( ω d , r , r ) as
E i ( ω d , r ) = d r G i j ( ω d , r , r ) J j ( ω d , r ) ,
where the Latin indices i , j = 1 , 2 , 3 enumerate the components of the fields and tensors and a summation is made over the repeated indices. Excluding H from (B2) and (B3), and expressing E via the Green tensor with (B4), one finds the equation for G i l :
i ( j G j l ) ( 2 + k d 2 ) G i l = i 4 π k d c δ i l δ ( r r ) ,
where the derivatives are taken on the argument r. Similar equation has been used to derive the Casimir force in [8].
Since in the ( x , y ) -plane the problem is homogeneous, it is useful to separate the lateral components and denote them by Greek indices α , β , = 1 , 2 , so that x 1 = x , x 2 = y , k 1 = k x , and k 2 = k y . The homogeneity also means that the Green tensor depends only on the difference x α x α and one can make the Fourier expansion in these coordinates
G i j ( ω d , r , r ) = d k ( 2 π ) 2 e i k α ( x α x α ) G i j ( ω d , k , z , z ) ,
where k = ( k 1 , k 2 ) . Applying this expansion to (B5), one finds a closed-form equation for the lateral components of the tensor
2 z 2 + q ˜ d 2 δ α γ + k α k γ q ˜ d 2 G γ β ( ω d , k , z , z ) = i 4 π k d c δ α β δ ( z z ) ,
where q ˜ d = k d 2 k 2 . The solution of this equation is
G α β = δ α β k α k β k d 2 G ( ω d , k , z , z ) ,
where the function G ( ω d , k , z , z ) is the solution of the equation
2 z 2 + q ˜ d 2 G ( ω d , k , z , z ) = i 4 π k d c δ ( z z )
with appropriate boundary conditions.
If the magnetic moment of the source is directed along the z-axis m = ( 0 , 0 , m 0 ) , then this source generates only TE (s) polarization. For this case it is sufficient to know only the lateral components of the Green tensor. Using (B1), (B4), and (B8), one finds the electric field components in the presence of metallic plate
E α ( p ) ( ω d , k , z ) = i c m 0 ϵ α β k β G ( ω d , k , z , h ) , E z = 0 ,
where ϵ α β is the 2D antisymmetric tensor (Levi-Civita symbol). Similarly, using (B2) the magnetic field is
H α ( p ) ( ω d , k , z ) = c m 0 k α k d z G ( ω d , k , z , h ) ,
H z ( p ) ( ω d , k , z ) = i c m 0 k 2 k d G ( ω d , k , z , h ) .
Thus, it is desirable to solve the boundary problem for (B9). The first boundary condition is G 0 for z . Furthermore, at the surface of the reflecting plane the lateral components of the electric and magnetic fields have to be continuous. However, such a boundary condition is very inconvenient for a spatially nonlocal material response. Since we are interested only in the field above the plane, it is more constructive to use the surface impedance boundary conditions, when there is no need to know the field in the material. For the TE (s) polarization the impedance is supposed to be a known function Z TE ( ω d , k ) that is defined as the ratio
Z TE ( ω d , k ) = E α ( ω d , k 0 ) ϵ α β H β ( ω d , k 0 ) .
Using here the expressions (B10) and (B11), one finds the second boundary condition
i k d G ( ω d , k , 0 , h ) + Z TE ( ω d , k ) z G ( ω d , k , z , h ) z = 0 = 0 .
The solution of the boundary problem for (B9) is the following:
G ( ω d , k , z , h ) = 2 π k d c q ˜ d r TE e i q ˜ d ( z + h ) e i q ˜ d | z h | ,
where r TE is the reflection coefficient expressed via the impedance of the material (B12) by (20).
Substituting the solution (B14) in (B11), we reproduce the expressions (30).

References

  1. Casimir, H.B.G. On the attraction between two perfectly conducting plates. Proc. Kon. Ned. Akad. Wet. B 1948, 51, 793–795. [Google Scholar]
  2. Mostepanenko, V.M.; Trunov, N.N. The Casimir Effect and Its Applications; Clarendon Press: Oxford, UK, 1997. [Google Scholar]
  3. Milton, K.A. The Casimir Effect: Physical Manifestations of Zero-Point Energy; World Scientific: Singapore, 2001. [Google Scholar]
  4. Bordag, M.; Klimchitskaya, G.L.; Mohideen, U.; Mostepanenko, V.M. Advances in the Casimir Effect; Oxford University Press: Oxford, UK, 2015. [Google Scholar]
  5. Sernelius, B.E. Fundamentals of van der Waals and Casimir Interactions; Springer: New York, NY, USA, 2018. [Google Scholar]
  6. Lifshitz, E.M. The theory of molecular attractive forces between solids. Zh. Eksp. Teor. Fiz. 1955, 29, 94–110, Translated: Sov. Phys. JETP 1956, 2, 73–83. [Google Scholar]
  7. Dzyaloshinskii, I.E.; Lifshitz, E.M.; Pitaevskii, L.P. The general theory of van der Waals forces. Usp. Fiz. Nauk 1961, 73, 381–422, Translated: Adv. Phys. 1961, 10, 165–209. [Google Scholar] [CrossRef]
  8. Lifshitz, E.M.; Pitaevskii, L.P. Statistical Physics, Part II; Pergamon: Oxford, UK, 1980. [Google Scholar]
  9. Emig, T.; Graham, N.; Jaffe, R.L.; Kardar, M. Casimir forces between compact objects: The scalar case. Phys. Rev. D 2008, 77, 025005. [Google Scholar] [CrossRef] [Green Version]
  10. Kenneth, O.; Klich, I. Casimir forces in a T-operator approach. Phys. Rev. B 2008, 78, 014103. [Google Scholar] [CrossRef] [Green Version]
  11. Rahi, S.J.; Emig, T.; Graham, N.; Jaffe, R.L.; Kardar, M. Scattering theory approach to electromagnetic Casimir forces. Phys. Rev. D 2009, 80, 085021. [Google Scholar] [CrossRef] [Green Version]
  12. Lamoreaux, S.K. Demonstration of the Casimir Force in the 0.6 to 6 μm Range. Phys. Rev. Lett. 1997, 78, 5–8. [Google Scholar] [CrossRef] [Green Version]
  13. Lamoreaux, S.K. Reanalysis of Casimir force measurements in the 0.6-to-6-μm range. Phys. Rev. A 2010, 82, 024102. [Google Scholar] [CrossRef]
  14. Decca, R.S.; Fischbach, E.; Klimchitskaya, G.L.; Krause, D.E.; López, D.; Mostepanenko, V.M. Improved tests of extra-dimensional physics and thermal quantum field theory from new Casimir force measurements. Phys. Rev. D 2003, 68, 116003. [Google Scholar] [CrossRef] [Green Version]
  15. Decca, R.S.; López, D.; Fischbach, E.; Klimchitskaya, G.L.; Krause, D.E.; Mostepanenko, V.M. Precise comparison of theory and new experiment for the Casimir force leads to stronger constraints on thermal quantum effects and long-range interactions. Ann. Phys. 2005, 318, 37–80. [Google Scholar] [CrossRef] [Green Version]
  16. Decca, R.S.; López, D.; Fischbach, E.; Klimchitskaya, G.L.; Krause, D.E.; Mostepanenko, V.M. Tests of new physics from precise measurements of the Casimir pressure between two gold-coated plates. Phys. Rev. D 2007, 75, 077101. [Google Scholar] [CrossRef] [Green Version]
  17. Decca, R.S.; López, D.; Fischbach, E.; Klimchitskaya, G.L.; Krause, D.E.; Mostepanenko, V.M. Novel constraints on light elementary particles and extra-dimensional physics from the Casimir effect. Eur. Phys. J. C 2007, 51, 963–975. [Google Scholar] [CrossRef]
  18. Bimonte, G.; López, D.; Decca, R.S. Isoelectronic determination of the thermal Casimir force. Phys. Rev. B 2016, 93, 184434. [Google Scholar] [CrossRef] [Green Version]
  19. Bimonte, G.; Spreng, B.; Maia Neto, P.A.; Ingold, G.-L.; Klimchitskaya, G.L.; Mostepanenko, V.M.; Decca, R.S. Measurement of the Casimir Force between 0.2 and 8 μm: Experimental Procedures and Comparison with Theory. Universe 2021, 7, 93. [Google Scholar] [CrossRef]
  20. Chang, C.-C.; Banishev, A.A.; Castillo-Garza, R.; Klimchitskaya, G.L.; Mostepanenko, V.M.; Mohideen, U. Gradient of the Casimir force between Au surfaces of a sphere and a plate measured using an atomic force microscope in a frequency-shift technique. Phys. Rev. B 2012, 85, 165443. [Google Scholar] [CrossRef] [Green Version]
  21. Banishev, A.A.; Chang, C.-C.; Klimchitskaya, G.L.; Mostepanenko, V.M.; Mohideen, U. Measurement of the gradient of the Casimir force between a nonmagnetic gold sphere and a magnetic nickel plate. Phys. Rev. B 2012, 85, 195422. [Google Scholar] [CrossRef] [Green Version]
  22. Banishev, A.A.; Klimchitskaya, G.L.; Mostepanenko, V.M.; Mohideen, U. Demonstration of the Casimir Force between Ferromagnetic Surfaces of a Ni-Coated Sphere and a Ni-Coated Plate. Phys. Rev. Lett. 2013, 110, 137401. [Google Scholar] [CrossRef] [Green Version]
  23. Banishev, A.A.; Klimchitskaya, G.L.; Mostepanenko, V.M.; Mohideen, U. Casimir interaction between two magnetic metals in comparison with nonmagnetic test bodies. Phys. Rev. B 2013, 88, 155410. [Google Scholar] [CrossRef] [Green Version]
  24. Xu, J.; Klimchitskaya, G.L.; Mostepanenko, V.M.; Mohideen, U. Reducing detrimental electrostatic effects in Casimir-force measurements and Casimir-force-based microdevices. Phys. Rev. A 2018, 97, 032501. [Google Scholar] [CrossRef] [Green Version]
  25. Liu, M.; Xu, J.; Klimchitskaya, G.L.; Mostepanenko, V.M.; Mohideen, U. Examining the Casimir puzzle with an upgraded AFM-based technique and advanced surface cleaning. Phys. Rev. B 2019, 100, 081406(R). [Google Scholar] [CrossRef] [Green Version]
  26. Liu, M.; Xu, J.; Klimchitskaya, G.L.; Mostepanenko, V.M.; Mohideen, U. Precision measurements of the gradient of the Casimir force between ultraclean metallic surfaces at larger separations. Phys. Rev. A 2019, 100, 052511. [Google Scholar] [CrossRef] [Green Version]
  27. Munday, J.N.; Capasso, F.; Parsegian, V.A.; Bezrukov, S.M. Measurements of the Casimir-Lifshitz force in fluids: The effect of electrostatic forces and Debye screening. Phys. Rev. A 2008, 78, 032109. [Google Scholar] [CrossRef]
  28. de Man, S.; Heeck, K.; Wijngaarden, R.J.; Iannuzzi, D. Halving the Casimir force with Conductive Oxides. Phys. Rev. Lett. 2009, 103, 040402. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Torricelli, G.; van Zwol, P.J.; Shpak, O.; Binns, C.; Palasantzas, G.; Kooi, B.J.; Svetovoy, V.B.; Wuttig, M. Switching Casimir forces with phase-change materials. Phys. Rev. A 2008, 82, 010101(R). [Google Scholar] [CrossRef] [Green Version]
  30. Sedighi, M.; Svetovoy, V.B.; Palasantzas, G. Casimir force measurements from silicon carbide surfaces. Phys. Rev. B 2016, 93, 085434. [Google Scholar] [CrossRef]
  31. Torgerson, J.R.; Lamoreaux, S.K. Low-frequency character of the Casimir force between metallic films. Phys. Rev. E 2004, 70, 047102. [Google Scholar] [CrossRef]
  32. Bimonte, G. Comment on “Low-frequency character of the Casimir force between metallic films”. Phys. Rev. E 2006, 73, 048101. [Google Scholar] [CrossRef] [Green Version]
  33. Svetovoy, V.B.; Esquivel, R. The Casimir free energy in high- and low-temperature limits. J. Phys. A Math. Gen. 2006, 39, 6777–6784. [Google Scholar] [CrossRef] [Green Version]
  34. Bezerra, V.B.; Klimchitskaya, G.L.; Mostepanenko, V.M.; Romero, C. Violation of the Nernst heat theorem in the theory of thermal Casimir force between Drude metals. Phys. Rev. A 2004, 69, 022119. [Google Scholar] [CrossRef] [Green Version]
  35. Bordag, M.; Pirozhenko, I. Casimir entropy for a ball in front of a plane. Phys. Rev. D 2010, 82, 125016. [Google Scholar] [CrossRef] [Green Version]
  36. Klimchitskaya, G.L.; Korikov, C.C. Analytic results for the Casimir free energy between ferromagnetic metals. Phys. Rev. A 2015, 91, 032119. [Google Scholar] [CrossRef] [Green Version]
  37. Klimchitskaya, G.L.; Mostepanenko, V.M. Low-temperature behavior of the Casimir free energy and entropy of metallic films. Phys. Rev. A 2017, 95, 012130. [Google Scholar] [CrossRef]
  38. Boström, S.; Sernelius, B.E. Entropy of the Casimir effect between real metal plates. Phys. A 2004, 339, 53–59. [Google Scholar] [CrossRef]
  39. Brevik, I.; Aarseth, J.B.; Høye, J.S.; Milton, K.A. Temperature dependence of the Casimir effect. Phys. Rev. E 2005, 71, 056101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Høye, J.S.; Brevik, I.; Ellingsen, S.A.; Aarseth, J.B. Analytical and numerical verification of the Nernst theorem for metals. Phys. Rev. E 2007, 75, 051127. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Bezerra, V.B.; Decca, R.S.; Fischbach, E.; Geyer, B.; Klimchitskaya, G.L.; Krause, D.E.; López, D.; Mostepanenko, V.M.; Romero, C. Comment on “Temperature dependence of the Casimir effect”. Phys. Rev. E 2006, 73, 028101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Svetovoy, V.B.; Lokhanin, M.V. Temperature correction to the Casimir force in cryogenic range and anomalous skin effect. Phys. Rev. A 2003, 67, 022113. [Google Scholar] [CrossRef] [Green Version]
  43. Klimchitskaya, G.L.; Mostepanenko, V.M. Experiment and theory in the Casimir effect. Contemp. Phys. 2006, 47, 131–144. [Google Scholar] [CrossRef] [Green Version]
  44. Bimonte, G.; Emig, T.; Kardar, M.; Krüger, M. Nonequilibrium Fluctuational Quantum Electrodynamics: Heat Radiation, Heat Transfer, and Force. Ann. Rev. Condens. Matter Phys. 2017, 8, 119–143. [Google Scholar] [CrossRef] [Green Version]
  45. Milton, K.A.; Li, Y.; Kalauni, P.; Parashar, P.; Guérout, P.; Ingold, G.-L.; Lambrecht, A.; Reynaud, S. Negative Entropies in Casimir and Casimir-Polder Interactions. Fortschr. Phys. 2017, 65, 1600047. [Google Scholar] [CrossRef] [Green Version]
  46. Klimchitskaya, G.L.; Mohideen, U.; Mostepanenko, V.M. The Casimir force between real materials: Experiment and theory. Rev. Mod. Phys. 2009, 81, 1827–1885. [Google Scholar] [CrossRef] [Green Version]
  47. Woods, L.M.; Dalvit, D.A.R.; Tkatchenko, A.; Rodriguez-Lopez, P.; Rodriguez, A.W.; Podgornik, R. Materials perspective on Casimir and van der Waals interactions. Rev. Mod. Phys. 2016, 88, 045003. [Google Scholar] [CrossRef]
  48. Mostepanenko, V.M. Casimir Puzzle and Conundrum: Discovery and Search for Resolution. Universe 2021, 7, 84. [Google Scholar] [CrossRef]
  49. Kliewer, K.L.; Fuchs, R. Anomalous Skin Effect for Specular Electron Scattering and Optical Experiments at Non-Normal Angles of Incidence. Phys. Rev. 1968, 172, 607–625. [Google Scholar] [CrossRef]
  50. Kats, E.I. Influence of non-locality effects on van der Waals interaction. Zh. Eksp. Teor. Fiz. 1977, 73, 212–220, Translated: Sov. Phys. JETP 1977, 46, 109–113. [Google Scholar]
  51. Esquivel, R.; Villarreal, C.; Mochán, W.L. Exact surface impedance formulation of the Casimir force: Application to spatially dispersive metals. Phys. Rev. A 2003, 68, 052103, Erratum: Phys. Rev. A 2005, 71, 029904.. [Google Scholar] [CrossRef] [Green Version]
  52. Esquivel, R.; Svetovoy, V.B. Correction to the Casimir force due to the anomalous skin effect. Phys. Rev. A 2004, 69, 062102. [Google Scholar] [CrossRef] [Green Version]
  53. Svetovoy, V.B.; Esquivel, R. Nonlocal impedances and the Casimir entropy at low temperatures. Phys. Rev. E 2005, 72, 036113. [Google Scholar] [CrossRef] [Green Version]
  54. Sernelius, B.E. Effects of spatial dispersion on electromagnetic surface modes and on modes associated with a gap between two half spaces. Phys. Rev. B 2005, 71, 235114. [Google Scholar] [CrossRef]
  55. Haakh, H.R.; Henkel, C. Magnetic near fields as a probe of charge transport in spatially dispersive conductors. Eur. Phys. J. B 2012, 85, 46. [Google Scholar] [CrossRef]
  56. Klimchitskaya, G.L.; Mostepanenko, V.M. An alternative response to the off-shell quantum fluctuations: A step forward in resolution of the Casimir puzzle. Eur. Phys. J. C 2020, 80, 900. [Google Scholar] [CrossRef]
  57. Klimchitskaya, G.L.; Mostepanenko, V.M. Casimir effect for magnetic media: Spatially non-local response to the off-shell quantum fluctuations. Phys. Rev. D 2021, 104, 085001. [Google Scholar] [CrossRef]
  58. Klimchitskaya, G.L.; Mostepanenko, V.M. Theory-experiment comparison for the Casimir force between metallic test bodies: A spatially non-local dielectric response. Phys. Rev. A 2022, 105, 012805. [Google Scholar] [CrossRef]
  59. Klimchitskaya, G.L.; Mostepanenko, V.M. Casimir entropy and non-local response functions to the off-shell quantum fluctuations. Phys. Rev. D 2021, 103, 096007. [Google Scholar] [CrossRef]
  60. Bordag, M.; Klimchitskaya, G.L.; Mostepanenko, V.M.; Petrov, V.M. Quantum field theoretical description for the reflectivity of graphene. Phys. Rev. D 2015, 91, 045037, Erratum: Phys. Rev. D 2016, 93, 089907. [Google Scholar] [CrossRef] [Green Version]
  61. Bordag, M.; Fialkovskiy, I.; Vassilevich, D. Enhanced Casimir effect for doped graphene. Phys. Rev. B 2016, 93, 075414, Erratum: Phys. Rev. B 2017, 95, 119905. [Google Scholar] [CrossRef] [Green Version]
  62. Liu, M.; Zhang, Y.; Klimchitskaya, G.L.; Mostepanenko, V.M.; Mohideen, U. Demonstration of an Unusual Thermal Effect in the Casimir Force from Graphene. Phys. Rev. Lett. 2021, 126, 206802. [Google Scholar] [CrossRef]
  63. Liu, M.; Zhang, Y.; Klimchitskaya, G.L.; Mostepanenko, V.M.; Mohideen, U. Experimental and theoretical investigation of the thermal effect in the Casimir interaction from graphene. Phys. Rev. B 2021, 104, 085436. [Google Scholar] [CrossRef]
  64. Klimchitskaya, G.L.; Mostepanenko, V.M. Quantum field theoretical description of the Casimir effect between two real graphene sheets and thermodynamics. Phys. Rev. D 2020, 102, 016006. [Google Scholar] [CrossRef]
  65. Klimchitskaya, G.L.; Mostepanenko, V.M. Casimir and Casimir-Polder Forces in Graphene Systems: Quantum Field Theoretical Description and Thermodynamics. Universe 2020, 6, 150. [Google Scholar] [CrossRef]
  66. Greffet, J.-J.; Carminati, R. Image formation in near-field optics. Prog. Surf. Sci. 1997, 56, 133–237. [Google Scholar] [CrossRef]
  67. Törmä, P.; Barnes, W.L. Strong coupling between surface plasmon polaritons and emitters: A review. Rep. Progr. Phys. 2015, 78, 013901. [Google Scholar] [CrossRef] [PubMed]
  68. Culshaw, W.; Jones, D.S. Effect of a Metal Plate on Total Reflection. Proc. Phys. Soc. B 1953, 66, 859–864. [Google Scholar] [CrossRef]
  69. Brady, J.J.; Brick, R.O.; Pearson, V.D. Penetration of Microwaves into the Rarer Medium in Total Reflection. J. Opt. Soc. Am. 1960, 50, 1080–1084. [Google Scholar] [CrossRef]
  70. Zhu, S.; Yu, A.W.; Hawley, D.; Roy, R. Frustrated total internal reflection: A demonstration and review. Am. J. Phys. 1986, 54, 601–606. [Google Scholar] [CrossRef]
  71. Hsu, J.W.P. Near-field scanning optical microscopy studies of electronic and photonic materials and devices. Mater. Sci. Engin R Rep. 2001, 33, 1–50. [Google Scholar] [CrossRef]
  72. Aigouy, L.; Lahrech, A.; Grésillon, S.; Cory, H.; Boccara, A.C.; Rivoal, J.C. Polarization effects in apertureless scanning near-field optical microscopy: An experimental study. Opt. Lett. 1999, 24, 187–189. [Google Scholar] [CrossRef]
  73. Klimchitskaya, G.L.; Mostepanenko, V.M.; Svetovoy, V.B. Probing the response of metals to low-frequency s-polarized evanescent fields. Europhys. Lett. 2022, 139, 66001. [Google Scholar] [CrossRef]
  74. Landau, L.D.; Lifshitz, E.M.; Pitaevskii, L.P. Electrodynamics of Continuous Media; Pergamon: Oxford, UK, 1984. [Google Scholar]
  75. Henkel, C.; Joulain, K.; Mulet, J.-P.; Greffet, J.J. Coupled surface polaritons and the Casimir force. Phys. Rev. A 2004, 69, 023808. [Google Scholar] [CrossRef] [Green Version]
  76. Intravaia, F.; Lambrecht, A. Surface Plasmon Modes and the Casimir Energy. Phys. Rev. Lett. 2005, 94, 110404. [Google Scholar] [CrossRef] [Green Version]
  77. Bordag, M. The Casimir effect for thin plasma sheets and the role of the surface plasmons. J. Phys. A Math. Gen. 2006, 39, 6173–6186. [Google Scholar] [CrossRef] [Green Version]
  78. Boström, M.; Sernelius, B.E. Thermal Effects on the Casimir Force in the 0.1–5 μm Range. Phys. Rev. Lett. 2000, 84, 4757–4760. [Google Scholar] [CrossRef] [PubMed]
  79. Bordag, M.; Geyer, B.; Klimchitskaya, G.L.; Mostepanenko, V.M. Casimir Force at Both Nonzero Temperature and Finite Conductivity. Phys. Rev. Lett. 2000, 85, 503–506. [Google Scholar] [CrossRef] [PubMed]
  80. Antezza, M.; Pitaevskii, L.P.; Stringari, S.; Svetovoy, V.B. Casimir-Lifshitz force out of thermal equilibrium. Phys. Rev. A 2008, 77, 022901. [Google Scholar] [CrossRef] [Green Version]
  81. Dressel, M.; Grüner, G. Electrodynamics of Solids: Optical Properties of Electrons in Metals; Cambridge University Press: Cambridge, UK, 2003. [Google Scholar]
  82. Silin, V.P.; Fetisov, E.P. Electromagnetic properties of a relativistic plasma, III. Zh. Eksp. Teor. Fiz. 1961, 41, 159–170, Translated: Sov. Phys. JETP 1962, 14, 115–122. [Google Scholar]
  83. Landau, L.D.; Lifshitz, E.M. The Classical Theory of Fields; Pergamon: Oxford, UK, 1971. [Google Scholar]
  84. Reitz, J.R. Forces on Moving Magnets due to Eddy Currents. J. Appl. Phys. 1970, 41, 2067–2071. [Google Scholar] [CrossRef]
  85. Saslow, W.M. Maxwell’s theory of eddy currents in thin conducting sheets, and applications to electromagnetic shielding and MAGLEV. Am. J. Phys. 1992, 60, 693–710. [Google Scholar] [CrossRef]
  86. Skaar, J. Fresnel’s equations in statics and quasistatics. Eur. J. Phys. 2019, 40, 045201. [Google Scholar] [CrossRef]
  87. Ordal, M.A.; Bell, R.J.; Alexander, R.W.; Long, L.L.; Querry, M.R. Optical properties of fourteen metals in the infrared and far infrared: Al, Co, Cu, Au, Fe, Pb, Mo, Ni, Pd, Pt, Ag, Ti, V, and W. Appl. Opt. 1985, 24, 4493–4499. [Google Scholar] [CrossRef]
  88. Chen, F.; Mohideen, U.; Klimchitskaya, G.L.; Mostepanenko, V.M. Experimental test for the conductivity properties from the Casimir force between metal and semiconductor. Phys. Rev. A 2006, 74, 022103. [Google Scholar] [CrossRef] [Green Version]
  89. Ulvr, M. Design of PCB search coils for AC magnetic flux density measurement. AIP Adv. 2018, 8, 047505. [Google Scholar] [CrossRef] [Green Version]
  90. Ramadan, Q.; Samper, V.; Poenar, D.; Yu, C. On-chip micro-electromagnets for magnetic-based bio-molecules separation. J. Magn. Magnet. Mater. 2004, 281, 150–172. [Google Scholar] [CrossRef]
  91. Wensink, H.; Benito-Lopez, F.; Hermes, D.C.; Verboom, W.; Gardeniers, H.J.G.E.; Reinhoudt, D.N.; van den Berg, A. Measuring reaction kinetics in a lab-on-a-chip by microcoil NMR. Lab Chip 2005, 5, 280–284. [Google Scholar] [CrossRef] [PubMed]
  92. Drndić, M.; Johnson, K.S.; Thywissen, J.H.; Prentiss, M.; Westervelt, R.M. Micro-electromagnets for atom manipulation. Appl. Phys. Lett. 1998, 72, 2906–2908. [Google Scholar] [CrossRef]
  93. Liu, Z.-X.; Wang, B.; Kong, C.; Si, L.-G.; Xiong, H.; Wu, Y. A proposed method to measure weak magnetic field based on a hybrid optomechanical system. Sci. Rep. 2017, 7, 12521. [Google Scholar] [CrossRef] [Green Version]
  94. Murzin, D.; Mapps, D.J.; Levada, K.; Belyaev, V.; Omelyanchik, A.; Panina, L.; Rodionova, V. Ultrasensitive Magnetic Field Sensors for Biomedical Applications. Sensors 2020, 20, 1569. [Google Scholar] [CrossRef] [Green Version]
  95. Huang, J.-H.; Duan, X.-Y.; Wang, G.-J.; Hu, X.-Y. Enhancing the precision of detecting weak magnetic fields based on weak-value amplification. J. Opt. Soc. Amer. B 2022, 39, 1289. [Google Scholar] [CrossRef]
  96. Korepanov, V.; Berkman, R.; Rakhlin, L.; Klymovych, Y.; Prystai, A.; Marussenkov, A.; Afanassenko, M. Advanced field magnetometers comparative study. Meas. J. Int. Meas. Confed. 2001, 29, 137–146. [Google Scholar] [CrossRef]
  97. Coillot, C.; Moutoussamy, J.; Lebourgeois, R.; Ruocco, S.; Chanteur, G. Principle and performance of a dual-band search coil magnetometer: A new instrument to investigate fluctuating magnetic fields in space. IEEE Sens. J. 2010, 10, 255–260. [Google Scholar] [CrossRef] [Green Version]
  98. Sudhakaran, K.; Komaragiri, R.; Suja, K.J. Fabrication, characterization, and modelling of a novel via-less single metal level magnetic microcoil sensor for biosensing applications. Sens. Actuators A Phys. 2019, 290, 190–197. [Google Scholar]
  99. Séran, H.C.; Fergeau, P. An optimized low-frequency three-axis search coil magnetometer for space research. Rev. Sci. Instr. 2005, 76, 044502. [Google Scholar] [CrossRef]
  100. Sternickel, K.; Braginski, A.I. Biomagnetism using SQUIDs: Status and perspectives. Supercond. Sci. Technol. 2006, 19, S160. [Google Scholar] [CrossRef]
  101. Yang, K.; Chen, H.; Lu, L.; Kong, X.; Yang, R.; Wang, J. SQUID Array with Optimal Compensating Configuration for Magnetocardiography Measurement in Different Environments. IEEE Trans. Appl. Supercond. 2019, 29, 1–7. [Google Scholar] [CrossRef]
  102. Buchner, M.; Höfler, K.; Henne, B.; Ney, V.; Ney, A. Tutorial: Basic principles, limits of detection, and pitfalls of highly sensitive SQUID magnetometry for nanomagnetism and spintronics. J. Appl. Phys. 2018, 124, 161101. [Google Scholar] [CrossRef] [Green Version]
  103. Asfour, A. (Ed.) Magnetic Sensors—Development Trends and Applications; IntechOpen: London, UK, 2017. [Google Scholar]
  104. Hannemann, M.; Wegner, G.; Henkel, C. No-Slip Boundary Conditions for Electron Hydrodynamics and the Thermal Casimir Pressure. Universe 2021, 7, 108. [Google Scholar] [CrossRef]
  105. Altland, A.; Simons, B.D. Condensed Matter Field Theory; Cambridge University Press: Cambridge, UK, 2010. [Google Scholar]
  106. Bordag, M.; Fialkovsky, I.V.; Khusnutdinov, N.; Vassilevich, D.V. Bulk contributions to Casimir interaction of Dirac materials. Phys. Rev. B 2021, 104, 195431. [Google Scholar] [CrossRef]
  107. Abramowitz, M.; Stegun, E.A. (Eds.) Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables; National Bureau of Standards: Washington, DC, USA, 1964.
  108. Prudnikov, A.P.; Brychkov, Y.A.; Marichev, O.I. Integrals and Series. Vol. 2: Special Functions; Gordon and Breach: New York, NY, USA, 1986. [Google Scholar]
  109. Prudnikov, A.P.; Brychkov, Y.A.; Marichev, O.I. Integrals and Series. Vol. 1: Elementary Functions; Gordon and Breach: New York, NY, USA, 1986. [Google Scholar]
Figure 1. Magnetic dipole spaced in a vacuum above thick metallic plate described by the dielectric permittivity ε .
Figure 1. Magnetic dipole spaced in a vacuum above thick metallic plate described by the dielectric permittivity ε .
Universe 08 00574 g001
Figure 2. The normalized magnitude of real part of the lateral component of monochromatic magnetic field reflected from the copper plate computed using the Drude model as the function of separation from the magnetic dipole at the dipole height of 1 cm above a plate is shown by the solid lines labeled 1, 2, 3, 4, and 5 for the oscillation frequency equal to 0.2, 2, 10, 20, and 100 rad/s, respectively. The top solid and overlapping with it dashed lines are computed using the plasma and spatially nonlocal models.
Figure 2. The normalized magnitude of real part of the lateral component of monochromatic magnetic field reflected from the copper plate computed using the Drude model as the function of separation from the magnetic dipole at the dipole height of 1 cm above a plate is shown by the solid lines labeled 1, 2, 3, 4, and 5 for the oscillation frequency equal to 0.2, 2, 10, 20, and 100 rad/s, respectively. The top solid and overlapping with it dashed lines are computed using the plasma and spatially nonlocal models.
Universe 08 00574 g002
Figure 3. The normalized imaginary part of the lateral component of monochromatic magnetic field reflected from the copper plate computed using the Drude model as the function of separation from the magnetic dipole at the dipole height of 1 cm above a plate is shown by the solid lines labeled 1, 2, 3, 4, and 5 for the oscillation frequency equal to 0.2, 2, 10, 20, and 100 rad/s, respectively.
Figure 3. The normalized imaginary part of the lateral component of monochromatic magnetic field reflected from the copper plate computed using the Drude model as the function of separation from the magnetic dipole at the dipole height of 1 cm above a plate is shown by the solid lines labeled 1, 2, 3, 4, and 5 for the oscillation frequency equal to 0.2, 2, 10, 20, and 100 rad/s, respectively.
Universe 08 00574 g003
Figure 4. The normalized magnitudes of real part of the z-component of monochromatic magnetic field above the copper plate computed using the Drude model and by the plasma and spatially nonlocal models at the dipole height of 1 cm for the oscillation frequency of 100 rad/s as the functions of separation from the magnetic dipole are shown by the solid line and by another solid line overlapping with the short-dashed one, respectively. The long-dashed line demonstrates the normalized magnitude of real part of the z-component of magnetic field of oscillating dipole in the absence of metallic plate. The inset shows the ratio of real part of the z-component of reflected field to the same quantity for the dipole field in the absence of metallic plate computed using the Drude model (the solid line) and the plasma or spatially nonlocal models (the solid line overlapping with the short-dashed one).
Figure 4. The normalized magnitudes of real part of the z-component of monochromatic magnetic field above the copper plate computed using the Drude model and by the plasma and spatially nonlocal models at the dipole height of 1 cm for the oscillation frequency of 100 rad/s as the functions of separation from the magnetic dipole are shown by the solid line and by another solid line overlapping with the short-dashed one, respectively. The long-dashed line demonstrates the normalized magnitude of real part of the z-component of magnetic field of oscillating dipole in the absence of metallic plate. The inset shows the ratio of real part of the z-component of reflected field to the same quantity for the dipole field in the absence of metallic plate computed using the Drude model (the solid line) and the plasma or spatially nonlocal models (the solid line overlapping with the short-dashed one).
Universe 08 00574 g004
Figure 5. The normalized magnitude of total monochromatic magnetic field above the copper plate computed using the Drude model at the dipole height of 1 cm for the oscillation frequency of 0.2 and 100 rad/s as the function of separation from the magnetic dipole is shown by the solid lines 1 and 2, respectively. Another solid line overlapping with the dashed one shows the same quantity computed using the plasma and spatially nonlocal models. The inset shows the ratio of magnitudes of the reflected and total field computed using the Drude model (lines 1 and 2 for the dipole oscillation frequencies of 0.2 and 100 rad/s, respectively) and the plasma and spatially nonlocal models (the solid line overlapping with the dashed one).
Figure 5. The normalized magnitude of total monochromatic magnetic field above the copper plate computed using the Drude model at the dipole height of 1 cm for the oscillation frequency of 0.2 and 100 rad/s as the function of separation from the magnetic dipole is shown by the solid lines 1 and 2, respectively. Another solid line overlapping with the dashed one shows the same quantity computed using the plasma and spatially nonlocal models. The inset shows the ratio of magnitudes of the reflected and total field computed using the Drude model (lines 1 and 2 for the dipole oscillation frequencies of 0.2 and 100 rad/s, respectively) and the plasma and spatially nonlocal models (the solid line overlapping with the dashed one).
Universe 08 00574 g005
Figure 6. The magnitude of the real part (a) and the imaginary part (b) of the lateral component of monochromatic magnetic field reflected from the copper plate computed using the Drude model at the dipole height of 1 cm above a plate is shown by the solid lines labeled 1, 2, and 3 as the function of oscillation frequency at separations from the magnetic dipole equal to 10, 20, and 30 mm, respectively. The three solid lines overlapping with the dashed ones counted from top to bottom (a) are computed using the plasma and spatially nonlocal models at separations of 10, 20, and 30 nm, respectively.
Figure 6. The magnitude of the real part (a) and the imaginary part (b) of the lateral component of monochromatic magnetic field reflected from the copper plate computed using the Drude model at the dipole height of 1 cm above a plate is shown by the solid lines labeled 1, 2, and 3 as the function of oscillation frequency at separations from the magnetic dipole equal to 10, 20, and 30 mm, respectively. The three solid lines overlapping with the dashed ones counted from top to bottom (a) are computed using the plasma and spatially nonlocal models at separations of 10, 20, and 30 nm, respectively.
Universe 08 00574 g006
Figure 7. The magnitudes of real part (a) and the imaginary part (b) of the normal component of monochromatic magnetic field in the domain above the copper plate computed using the Drude model at the dipole height of 1 cm are shown by the solid lines labeled 1, 2, and 3 as the functions of oscillation frequency at separations from the magnetic dipole equal to 30, 40, and 50 mm, respectively. The three solid lines labeled 1, 2, and 3 overlapping with the dashed ones (a) are computed using the plasma and spatially nonlocal models at separations of 30, 40, and 50 nm, respectively.
Figure 7. The magnitudes of real part (a) and the imaginary part (b) of the normal component of monochromatic magnetic field in the domain above the copper plate computed using the Drude model at the dipole height of 1 cm are shown by the solid lines labeled 1, 2, and 3 as the functions of oscillation frequency at separations from the magnetic dipole equal to 30, 40, and 50 mm, respectively. The three solid lines labeled 1, 2, and 3 overlapping with the dashed ones (a) are computed using the plasma and spatially nonlocal models at separations of 30, 40, and 50 nm, respectively.
Universe 08 00574 g007
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Klimchitskaya, G.L.; Mostepanenko, V.M.; Svetovoy, V.B. Experimentum crucis for Electromagnetic Response of Metals to Evanescent Waves and the Casimir Puzzle. Universe 2022, 8, 574. https://doi.org/10.3390/universe8110574

AMA Style

Klimchitskaya GL, Mostepanenko VM, Svetovoy VB. Experimentum crucis for Electromagnetic Response of Metals to Evanescent Waves and the Casimir Puzzle. Universe. 2022; 8(11):574. https://doi.org/10.3390/universe8110574

Chicago/Turabian Style

Klimchitskaya, Galina L., Vladimir M. Mostepanenko, and Vitaly B. Svetovoy. 2022. "Experimentum crucis for Electromagnetic Response of Metals to Evanescent Waves and the Casimir Puzzle" Universe 8, no. 11: 574. https://doi.org/10.3390/universe8110574

APA Style

Klimchitskaya, G. L., Mostepanenko, V. M., & Svetovoy, V. B. (2022). Experimentum crucis for Electromagnetic Response of Metals to Evanescent Waves and the Casimir Puzzle. Universe, 8(11), 574. https://doi.org/10.3390/universe8110574

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop