Next Article in Journal / Special Issue
Axion Electrodynamics and the Casimir Effect
Previous Article in Journal
Electric Octupole-Dependent Contributions to Optical Binding Energy
Previous Article in Special Issue
Casimir Physics beyond the Proximity Force Approximation: The Derivative Expansion
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Casimir Forces with Periodic Structures: Abrikosov Flux Lattices

by
Shunashi Guadalupe Castillo-López
,
Raúl Esquivel-Sirvent
,
Giuseppe Pirruccio
and
Carlos Villarreal
*
Instituto de Física, Universidad Nacional Autónoma de México, Apdo. Postal 20-364, Mexico City 01000, Mexico
*
Author to whom correspondence should be addressed.
Physics 2024, 6(1), 394-406; https://doi.org/10.3390/physics6010026
Submission received: 4 December 2023 / Revised: 27 January 2024 / Accepted: 30 January 2024 / Published: 12 March 2024
(This article belongs to the Special Issue 75 Years of the Casimir Effect: Advances and Prospects)

Abstract

:
We investigate the influence of the Abrikosov vortex lattice on the Casimir force in a setup constituted by high-temperature superconductors subject to an external magnetic field. The Abrikosov lattice is a property of type II superconductors in which normal and superconducting carriers coexist and the latter define a periodic pattern with square symmetry. We find that the optical properties determined by spatial redistribution of the superconducting order parameter induce Casimir forces with a periodic structure whose minimal strengths coincide with the vortex cores.

1. Introduction

Seventy-five years have passed since, motivated by a suggestion by Niels Bohr during a walk, Hendrik Casimir proposed that vacuum fluctuations could induce an attractive force per unit area between two perfectly conducting parallel plates, a distance d apart, given by F = c π 2 / 240 d 4 [1] with = h / 2 π the reduced Planck constant and c the speed of light. A more realistic theory was proposed by Evgeny Lifshitz in 1956 by considering fluctuating electrodynamics, based on the fluctuation–dissipation mechanism. Lifshitz theory allows to determine the Casimir force in terms of the dispersive and dissipative properties of the materials [2], as described by its optical properties. Grounded on Lifshitz formulation, numerous experiments [3,4,5,6,7,8,9,10,11,12,13] have been performed on measuring the Casimir forces involving a diversity of experimental arrangements and materials [14,15]. The influence and taming of these forces in the design and construction of micro- and nanodevices is a current field of research. Extensive treatments of the Casimir effect and its applications, are presented in Refs. [16,17,18,19,20].
In spite of the fruitful advances in the investigations of the Casimir effect there exist yet pending fundamental problems on the basic theory, concerning the role of dissipative mechanisms on the strength of the force between metallic bodies. These may involve electron scattering by impurities, other electrons, phonons, etc., yielding a total scattering rate τ 1 = τ el imp 1 + τ el el 1 ( T ) + τ el ph 1 ( T ) + . Actually, the inclusion of this kind of contributions in the theoretical characterization of the optical response of materials involved in a given setup should be necessary to achieve congruence with the fluctuation–dissipation theorem underlying the Lifshitz theory. However, earlier measurements of Casimir forces in metals at room temperature at body separations d 50 –600 nm showed consistence with theoretical predictions if dissipative effects are neglected. This hypothesis is also consistent with more recent experimental studies of the gradient of the Casimir force between metallic surfaces performed at larger separations, d > 700 nm, such that k B T c / d [21,22,23], with k B the Boltzmann constant. On the contrary, experiments carried out at micrometer distances to measure the total magnitude of the Casimir force display a better agreement with predictions including electronic relaxation [19]. The main difference between the results at small and large separations, is that in the latter thermal effects become important.
It has been proposed that the study of the Casimir effect in superconducting (SC) materials may constitute an excellent scenario to asses the influence of relaxation phenomena on the strength of the Casimir force between metallic bodies [24,25,26,27,28]. This is motivated by the finding that charge carriers in these materials exhibit a transition from dissipative transport to a dissipationless coherent behavior at a critical temperature T = T c . However, measurements of the influence of the SC transition on the Casimir force in setups involving conventional BCS (Bardeen–Cooper–Schrieffer) superconductors turn out to be extremely difficult, since for typical values T c 1 K, and k B T c c / d for sub- μ m body separations. Therefore, indirect approaches have been proposed based on observations of the Casimir-induced shift of the critical magnetic field H c of a thin superconducting film, or differential measurements of the Casimir force [29,30].
This suggests that the use of high-temperature superconductors (HTSCs), with T c 100 K, could constitute a suitable alternative to perform a direct analysis of the effect of the SC transition on the Casimir effect. In previous studies, we investigated the Casimir forces between objects made of optimally doped YBa 2 Cu 3 O 7 δ (YBCO), with T c = 93 K, either in thermal [31], or out of thermal equilibrium [32]. In the first case, we found that the Casimir force displays an abrupt increment as as T approaches T c from above, T T c . On the other hand, for T < T c , the (force × distance) approaches a constant in the limit T T c . In the second case, each slab was in local equilibrium with a thermal reservoir at respective temperatures, T 1 = 300 K and T 2 , where 300 T 2 0 K. In contrast with the thermal equilibrium situation, the Casimir force displays an abrupt decrement in the transit from normal metal to the SC state as T 2 T c . The low-temperature asymptotic behavior of the force is similar to that displayed in the equilibrium situation.
To get further insight on the influence of superconductivity-related effects on the Casimir effect, in this paper, we study the effect of the Abrikosov lattice (AL) [33] on the local properties of Casimir forces associated to high-temperature superconductors (HTSCs). The AL is a manifestation of the Meissner effect, in which the presence of an external magnetic field induces surface screening supercurrents, which expel out the magnetic field lines from the material’s interior within a London penetration length λ L ( T ) n s ( T ) 1 / 2 . Here, n s ( T ) is the number density of Cooper pairs (CPs) at a temperature T. In the case of type-II superconductors, such as YBCO, the Meissner effect involves the existence of a mixed phase of coexistence of normal and SC charge carriers determined by two critical magnetic fields, H c 1 < H c 2 . For values of the applied field higher than H c 1 , magnetic flux lines penetrate the sample in the form of quantum vortexes, Φ 0 = h / 2 e , with e the elemantary charge, thus inducing local screening currents to overcome the applied field [34]. Upon increasing magnitude of the field, the vortex density increases and saturates at the upper critical field H c 2 , where superconductivity disappears. Remarkably, as shown by Alexei Abrikosov [33], for intensities of the applied field just below H c 2 the vortexes align in a compact square lattice with period L x = L y = 2 π ξ ( T ) , where ξ ( T ) is the CP coherence length. In the case of YBCO, ξ ( 0 ) ξ 0 1.65 nm, and λ L ( 0 ) λ 0 156 nm, while the temperature-independent ratio κ = λ L ( T ) / ξ ( T ) 95 [34].
It can be shown that the AL vortexes strongly repel each other, giving rise to highly correlated configurations which are stable when thermal fluctuations and disorder are both negligible. A measure of the magnitude of the energy associated to thermal fluctuations with respect to the magnetic condensation energy is provided by the Ginzburg number [35], G i = 2 γ 2 κ 4 k B T c / H c 2 ( 0 ) 2 ξ 0 3 2 , where γ is a measure of the anisotropic conductivity (discussed below). In the case of conventional BCS superconductors, G i = 10 7 . In comparison, in the case of HTSCs, T c 10 2 K and κ 10 2 , implying that G i 10 2 . This relatively large value of G i joined with the feature that these materials display a layered anisotropic structure at the atomic level, leads to the manifestation of a manifold of phenomena generally termed as vortex matter, encompassing a complex phase diagram under different environmental conditions and material compositions [36]. Thus, thermal fluctuations may significantly alter the properties of the AL, generally leading to melting towards a liquid phase displaying vortex deformation, entanglement or migration. Superposed with repulsive interactions and thermal fluctuations, disorder due to material imperfections induce vortex pinning, which may conduce to the formation of glassy configurations [35,36,37]. Vortex matter has been investigated by recurring to techniques such as scanning tunneling microscopy [38] or muon spin rotation [39]. In addition, the possibility of using cold atoms to asses vortex-noise has been considered in the literature [40].
In order to examine the influence of the AL on the Casimir force, in this study we consider a setup, depicted in Figure 1, constituted by a spherical Au nanoparticle located at a minimum distance d from a planar YBCO substrate, in presence of an applied magnetic field directed along the z-axis. We show that the force acquires a spatial structure congruent with the AL due to the modulation imprinted by the vortexes on the dielectric permittivity.
In what follows, we present the formalism aimed to evaluate the Casimir force between a planar substrate and a nanosphere, which relies upon the frequency-dependent optical properties of the involved materials. To investigate the optical response of YBCO under the action of an external magnetic field, we then discuss a generalization of the Ginzburg–Landau (GL) theory of superconductivity, which allows the consideration of the anisotropic properties HTSCs, allowing the characterization of the number density of SC pairs, n s ( r ) , with a spatial dependence induced by the action of the applied magnetic field, as provided by the Abrikosov solution of the GL equation. We then discuss the thermal behavior of n s ( r ) by taking into account that SC charge carriers in HTSCs may be described as a 2D (2-dimensional) gas of weakly-interacting CPs able to form a Bose–Einstein condensate (BEC). In Section 2, we model the optical response of YBCO by combining the derived expression for n s ( r , T ) with experimental data for the YBCO dielectric function obtained in the normal ( T = 100 K) and SC regime ( T = 2 K). We then integrate the former antecedents to evaluate the Casimir force associated to the AL, and discuss the derived results.

2. Theory and Definitions

2.1. Casimir Force between a Nanosphere and a Planar Substrate

The theory of the Casimir effect between a sphere and a planar surface beyond the Proximity Force Approximation has been investigated within alternative perspectives, including some developed by authors of the present study [18,41,42,43,44,45,46,47]. In this Section, we extend the formalism previously presented in Refs. [41,42,43] to calculate the finite-temperature Casimir force for the nanosphere-substrate setup displayed in Figure 1, with the corresponding dielectric function of the nanosphere, ε ns ( ω ) , and substrate, ε sub ( r , ω ) . We first evaluate the zero-temperature interaction energy as a sum over proper frequencies, ω k ( r ; d ) , of the considered configuration:
E ( r ; d ) = 1 2 k ω k ( r ; d ) ω k ( r ; d ) ,
where r ( x , y ) . Straightforward use of the argument principle method lets us express the sum over proper mode frequencies as a sum over the zeros of a spectral function G ( ω ; r ; d ) (discussed below). This is determined by the solutions of Maxwell equations with boundary conditions satisfied by the plate-sphere setup:
E ( r ; d ) = 1 2 π i C ω 2 ω log G ( ω ; r , d ) log G ( ω ; r , d ) d ω .
Here, the (counter clock-wise) contour C is defined along the imaginary axis of the complex ω -plane and a semicircle in the right hand of this plane with its radius tending to infinity. The integral along the semicircle yields a null contribution, and the integral in Equation (3) may be evaluated by considering a contour between i and i . An integration by parts leads to
E ( r ; d ) = 4 π log G ( ω ; r , d ) log G ( i ζ ; r , ) d ξ .
The Casimir force is then obtained by calculating the derivative F ( r ; d ) = E ( r ; d ) / d . The eigenfrequency set of the sphere-substrate setup { ω k ( r ; d ) } is obtained by assuming that the vacuum fluctuations induce a charge distribution on the sphere, described at lowest level, as a point dipole moment located at its center:
p ns 0 ( ω ) = α ( ω ) E vac ( ω ) ,
where α ( ω ) = 4 π R 3 ( ε ns 1 ) / ( ε ns + 2 ) denotes the nanosphere polarizability with R the nanosphere radius and E vac denoting the electric field in vacuum. This dipole moment introduces in turn a charge distribution in the YBCO half-space. By using the images method it follows that the total induced dipole moment on the sphere is
p ns ( ω ; r ) = α ( ω ) E vac ( ω ) + T · p sub ( r ) .
Here, T is the dipole-dipole interaction tensor T = ( 3 r 0 r 0 r 0 2 I ) / r 0 5 , with r 0 a vector joining the center of the sphere and its image dipole below the substrate surface. In turn, the relation between the dipole moment on the sphere and the dipole moment induced on the YBCO substrate is p sub ( ω ; r ) = f c ( ω ; r ) M · p ns ( r ) , where M = diag ( 1 , 1 , 1 ) is a diagonal matrix in the chosen coordinate frame, and the contrast factor f c ( ω ; r ) ( 1 ϵ sub ( ω ; r ) ) / ( 1 + ϵ sub ( ω ; r ) ) . By substituting p sub ( ω ; r ) ) into Equation (5), one obtains:
α 1 ( ω ) I + f c ( ω ; r ) M · T · p ns ( ω ; r ) ) = E vac ( ω ) ,
which, by introducing the function u ( ω ) = [ 1 ϵ ns ( ω ) ) ] 1 and explicitly substituting α ( ω ) , may be re-expressed as the secular equation:
u ( ω ) I + H ( r ) · p ns ( ω ; r ) = E ˜ vac ( ω ) ,
where H ( r ) = ( 1 / 3 ) [ I + R 3 f c ( ω ; r ) M · T ] , and E ˜ vac = ( 1 / 3 ) R 3 E vac . By performing the change of variable ω i ζ , it follows that the matrix H in Equation (7) is real. This allows us to introduce the spectral function such that
G ( i ζ ; r , d ) l [ u ( i ζ ) + η l ( i ζ ; r , d ) ] = 0 ,
which in the present case implies three eigenvalues,
η 1 , 2 ( i ζ ; r , d ) = 1 3 1 + f c ( i ζ ; r ) [ 2 ( 1 + d / R ) ] 3 ; η 3 ( i ζ ; r , d ) = 1 3 1 + 2 f c ( i ζ ; r ) [ 2 ( 1 + d / R ) ] 3 ,
whose structure reflects the dipole–dipole interaction described by the tensor T , and η 1 , 2 correspond to dipoles parallel to the surface, while η 3 is perpendicular. Substitution of G ( i ζ ; r , d ) into Equation (2) leads to a final expression for Casimir force at null temperature:
F ( r , d ) = 4 π d l log [ u ( i ζ ) + η l ( i ζ ; r , d ) ] d ζ = 16 π R 1 ( 1 + d / R ) 4 f c ( i ζ ; r ) u ( i ζ ) + η 1 ( i ζ ; r , d ) + f c ( i ζ ; r ) u ( i ζ ) + η 3 ( i ζ ; r , d ) d ζ ,
where the feature that η 1 ( i ζ ; r , d ) = η 2 ( i ζ ; r , d ) has been considered.
This result may be generalized to the finite-temperature regime by use of the Matsubara formalism. In this approach, the frequency integration is replaced by a summation over discrete frequencies ζ n = 2 π k B T n / , with n an integer number. In that case, the final expression for temperature-dependent Casimir force is:
F ( r , d ; T ) = k B T 4 R 1 ( 1 + d / R ) 4 n = 0 f c ( i ζ n ; r ) u ( i ζ n ) + η 1 ( i ζ n ; r ) + f c ( i ζ n ; r ) u ( i ζ n ) + η 3 ( i ζ n ; r ) ,
where the prime implies that the n = 0 term should be multiplied by a factor 1 / 2 .
The influence of higher-order multipoles on the Casimir force for the present configuration has been studied in Refs. [42,43]. It follows that for sphere–substrate multipolar interaction of order l and l , respectively, the interaction energy E ( r ; d ) f ( r ) / ( 2 [ R + d ] ) l + l + 1 , and therefore, F ( r ; d ) f ( r ) / ( 2 [ R + d ] ) l + l + 2 . However, specific calculations [42] indicate that the contributions for l , l > 1 , only become relevant for distances d / R < 2 , consistently with the dipolar approximation considered in this study.

2.2. Ginzburg–Landau Theory and the Optical Response of the YBCO Substrate

In the expressions derived in Section 2.1, the functions u ( ω ) and f c ( ω ) are respectively determined by the dielectric response of the nanosphere, ε ns , and the planar substrate, ε sub . The dielectric properties of the gold nanosphere may be straightforwardly represented by a Drude function ε ns ( ω ) = 1 ω A u 2 / ( ω 2 + i γ A u ω ) , where the gold plasma frequency, ω A u = 9.1 eV, and the inverse scattering rate, γ A u = 0.02 eV. Therefore, u ( i ζ ) = ( ξ 2 + γ A u ζ ) / ω A u 2 .
On the other hand, to characterize the dielectric response of the YBCO subject to the action of an applied magnetic field, H , we put forth a straightforward variation of the GL theory of superconductivity that takes into account specific SC features of HTSCs, such as YBCO. In the GL theory, the transit to the SC state is described as a second-order phase transition determined by a complex order parameter, Ψ ( r ) , null in the normal phase, but finite in the SC phase, characterizing a long-range order specified by the number density of SC pairs, n s ( r ) = | Ψ ( r ) | 2 . We incorporate in the GL formalism the following findings [48]: (i) cuprate superconductors exhibit a layered crystallographic structure in which superfluid transport of CPs occurs mainly along CuO 2 planes (the a b -planes), whereas CP transport along the perpendicular c-axis occurs due to Josephson tunneling, with respective effective masses, m a b and m c , reflecting an anisotropic charge transport measured by the ratio γ 2 = m c / m a b 50 , (ii) an extremely short CP coherence length, such that the GL parameter κ 1 . The finding that γ 2 1 implies that YBCO is an uniaxial material with a concomitant anisotropic dielectric response characterized by the diagonal tensor ε = diag ( ε a b , ε a b , ε c ) in the coordinate system depicted in Figure 1. However, a related study on the radiative heat transfer between nearby YBCO surfaces [49] revealed that electromagnetic mode contributions involving ε c are relevant only for thin film configurations. Therefore, in the following, we assume that ε sub = ε a b . On the other hand, the joint conditions (i) and (ii) are indicative of a strongly binding pair interaction, leading the formation of a 2D gas of weakly-interacting particles, able to form a BEC [50].
In this framework, the anisotropic free energy density in the SC state in presence of a magnetic field B ( r ) = × A ( r ) , with A the electromagnetic potential, can be expressed as
f s ( T ) = f n ( T ) + 1 2 ( π ^ Ψ ) · 1 m · ( π ^ Ψ ) + a T | Ψ | 2 + b 2 | Ψ | 4 + 1 2 μ 0 B 2 ( r ) ,
where f n ( T ) is the normal state contribution, Ψ is the superconducting order parameter, the symbol   denotes hermitian conjugate, μ 0 is magnetic permeability of free space, the canonical momentum π ^ = i + e A , with e = 2 e the CP charge, and the reciprocal mass tensor is also diagonal in the chosen coordinate system: ( 1 / m ) = diag ( 1 / m a b , 1 / m a b , 1 / m c ) . In the original GL approach, a T = a 0 ( T T c ) and b > 0 ; however, as shown in Section 2.3, a more accurate description of the physical properties of SC materials can be achieved by considering alternative temperature parameterizations.
In the SC state, the total energy is obtained by spatial integration of Equation (12), F s = f s ( T ) d 3 r . It follows that the functional differentiation, δ F s ( T ) / δ Ψ ( r ) , then leads to the anisotropic GL equation
2 2 i e A · 1 m · i e A Ψ + a T + b | Ψ | 2 Ψ = 0 ,
while the differentiation δ F s ( T ) / δ A ( r ) yields the current density
J s = i e 2 1 m · Ψ Ψ Ψ Ψ 1 m · e 2 A | Ψ | 2 .
It can be shown that the mass anisotropy induces in turn anisotropic coherence lengths, ξ j ( T ) , satisfying the relations ξ a b 2 ( T ) = 2 / 2 m a b | a T | , and ξ c 2 ( T ) = 2 / 2 m c | a T | [48]. In absence of external fields or material boundaries, the former equations predict a second-order phase transition, with an order parameter magnitude, | Ψ | 2 = 0 , for a T > 0 , and n s = | Ψ | 2 = a T / b , for a T < 0 , and a null current density. In this case, the condensation energy of the superconductor (per unit volume) is given by [48]
f s ( T ) f n ( T ) = a T 2 / 2 b = μ 0 H c 2 2 ( T ) / 4 κ 2 ,
with H c 2 ( T ) = Φ 0 / 2 π ξ a b 2 ( T ) . In the general case, Equations (13) and (14) represent a coupled equation system; however, for type II materials this system uncouples by taking into account that a thermodynamic phase transition also occurs for magnetic fields H = H c 2 . Since the phase transition is of second order, Ψ is small and so is the magnetization M . Therefore, B = μ 0 ( H + M ) μ 0 H , where H is the applied field. This implies that, near H c 2 , the spatial variation in B ( r ) can be neglected and consider that it has the form B = B e z , where e z is directed along the crystallographic c-axis. Then, in the Landau gauge, the vector potential A = B x e y . A further consequence of the smallness of Ψ is that the GL equation may be linearized by dropping cubic term in Equation (13). Then, the GL equation becomes
2 2 m a b i e B x e y 2 Ψ 2 2 m c 2 z 2 Ψ = | a T | Ψ .
As shown by Abrikosov [33], Equation (16) admits Landau-level solutions of the form Ψ ( r ) = Φ ( r ) e i k z z , where the ground state is given by k z = 0 , and
Φ ( r ) = n = C n e i ( 2 π n y ) / L y exp [ ( x + n Φ 0 / B L y ) 2 / ξ a b 2 ( T ) ] ,
which is a periodic function provided C n = C n + ν , for some integer ν . The case ν = 1 corresponds to a square lattice as displayed by HTSCs, while the case ν = 2 yields a triangular lattice, characteristic of conventional superconductors [34]. The stability of these solutions can be analyzed by incorporating the effects of nonlinear terms in the description. This is determined by the Abrikosov parameter, β A = | Ψ | 4 / | Ψ | 2 . In the homogeneous case, β A = 1 , whereas β A > 1 implies less favorable energies. In particular, the energy of the triangular configuration is smaller than that of the square one by less than 1 % [34]. In the square lattice configuration, the summation may be re-expressed in terms of a Jacobi theta function, θ 3 ( z ; q ) = n = q n 2 e 2 i n z , so that
Φ ( x ˜ , y ˜ ) = C e 1 2 κ 2 x ˜ 2 θ 3 1 ; 2 π κ i ( x ˜ + i y ˜ ) .
In this case, κ = λ a b / ξ a b where λ a b is the a b -plane penetration length, while x ˜ = x / λ a b , y ˜ = y / λ a b . Notice that λ a b describes the magnetic field screening by supercurrents flowing along the a b -plane. Figure 1 depicts the resulting contours of constant probability defined by | Φ ( r ˜ ) | 2 . We observe a lattice structure with square elementary cells with dimensionless periods L x = L y = 2 π / κ . By writing Φ ( r ˜ ) = | Φ ( r ˜ ) | e i χ ( r ˜ ) it follows that the GL current density is given by J s = ( e / 2 m a b ) | Φ | 2 ( χ ( e / ) A ) , indicating that super-current lines coincide with the equi-probability contours, being the vortex cores located at the darkest zones of the figure. Notice that the vortex size can be tuned by the substrate temperature T, which modulates the lattice parameter. In normal units, it follows that L x ( T ) = L y ( T ) = 2 π ξ a b ( T ) , so that L x ( 2 K ) 4 nm, whereas L x ( 90 K ) 16 nm.

2.3. Thermal Properties of the Order Parameter

We describe the thermal properties of the order parameter by recurring to the London two-fluid model of superconductivity, inspired in turn by the two-fluid model of the superfluid He 4 . This model assumes that the charge density, n, at a fixed position, is split as a sum of normal and SC contributions, n = n n ( T ) + n s ( T ) . The temperature behavior of n s ( T ) is frequently described according to the Casimir–Gorter model, n s ( T ) / n = 1 ( T / T c ) 4 [51]. However, this empirical relation is in clear disagreement with experimental observations [48].
We consider instead, as mentioned above, that CPs define a 2D gas of weakly-interacting particles with BE statistical properties. As a consequence, the energy excitation spectrum is given by Bogoliubov expression E k = ( c s k ) 2 + ( 2 k 2 / 2 m ) 2 1 / 2 which, in the low-momentum limit, leads to a phonon dispersion relation, E k c s k , with c s the sound’s speed, consistent with Landau’s criterion for superfluid particle transport [52]. In that case, we assume the pair occupancy density at any given position can be expressed, for T < T BEC , in the form [53]
n 2 D = n 0 2 D ( T ) + k 0 1 exp c s k / k B T 1 .
Here, T BEC = 2 π 2 c s 2 n 2 D / k B 2 ζ ( 2 ) 1 / 2 , and ζ ( 2 ) is a Riemann’s zeta function. The summation in Equation (19) can be straightforwardly evaluated by integrating over a 2D density of states. It follows that n 0 2 D ( T ) / n 2 D = 1 ( T / T c ) 2 [53,54]. Taking into account that in the dilute regime the condensate and the density of superfluid charge carriers coincide, we identify T BEC = T c . This latter result implies, in particular, that the penetration length satisfies the universal relation λ a b 2 ( 0 ) / λ a b 2 ( T ) = 1 ( T / T c ) 2 , which has shown to represent an accurate description of experimental data for YBCO systems under an ample range of doping levels in the whole temperature interval 0 < T < T c [54,55,56]. By assuming now that the former results hold at every position over the SC substrate, n 0 2 D ( r , T ) n s 2 D ( r , T ) , and then
n s 2 D ( r , T ) = | Φ ( r ) | 2 1 ( T / T c ) 2 .

2.4. YBCO Dielectric Response

The optical properties of HTSCs have been experimentally investigated for different compounds at several temperatures and frequencies using reflectivity and impedance-type measurements [57,58,59]. In particular, the dielectric function, ε ( ω ) , of YBCO samples has been measured in the normal and SC states at T = 100 K and T = 2 K, respectively. Notably, the estimated values of the optical parameters remain practically unaltered in these temperature extremes. In particular, the plasma frequency, ω p 2 = e 2 n / ε 0 m , satisfies ω p ( 2 K ) = ω p ( 100 K ) = 0.75 eV, suggesting that London’s two-fluid model of superconductivity can be employed to derive an interpolation, ε ( ω ; T ) , valid in the entire temperature range 2 < T < 100 K. For T > T c , an accurate representation of the dielectric response includes a constant term, ε , in the optical spectrum, as well as Drude, mid-infrared (MIR), and N ph phonon contributions:
ε a b ( n ) ( ω ) = ε ω p 2 ω 2 + i γ ( T ) ω Ω MIR 2 ω 2 ω MIR 2 + i Γ MIR ω r = 1 N ph Ω ph 2 ω 2 ω r 2 + i γ r ω .
Here, ε = 3.8 , the inverse scattering rate γ ( T ) = 0.037 + γ 1 T eV, with γ 1 = 8 × 10 15 eV/K, and the MIR parameters Ω MIR = 2.6 eV, ω MIR = 0.26 eV, Γ MIR = 1 eV, whereas the phonon parameters are given in Refs. [57,58]. In the SC regime, dissipative scattering does not occur, so that γ 0 . In that limit, ω ± i γ 1 P 1 / ω i π δ ( ω ) , with δ ( ω ) the Dirac delta function, and the dielectric function becomes:
ε a b ( s ) ( ω ; r ) = ε i π ω p 2 2 ω δ ( ω ) + ω p 2 ω 2 1 ( T / T c ) 2 | Ψ ( r ) | 2 ω p 2 ( T / T c ) 2 ω 2 + i γ ( T ) ω Ω MIR 2 ω 2 ω MIR 2 + i Γ MIR ω r = 1 N ph Ω ph 2 ω 2 ω r 2 + i γ r ω .
As explained in Section 2.1, the dielectric function has to be calculated in the rotated frequency space i ζ where the details can be found in Ref. [60]. In this scheme, the nanosphere and substrate permittivities are respectively given by ε ns ( i ζ ) = ε A u ( i ζ ) , ε sub ( i ζ , r , T ) = ε a b ( n ) ( i ζ , r , T > T c ) and ε sub ( i ζ , r , T ) = ε a b ( s ) ( i ζ , r , T < T c ) .

3. Results

We show in Figure 2 the structure of Casimir force at T = 2 K, as a function of the position of the Au nanosphere over the AL. Here, the nanosphere radius R = 4 nm, and d = 2 R . Although the dipole approximation requires, actually, d R , as mentioned above, in a previous study, it was found that higher multipolar contributions to the force become negligible at this separation [42]. This figure reveals that the Casimir force displays a periodic structure congruent with the spatial charge distribution induced by the AL. It can be observed that the modulation amplitude Δ F = | F max F min | is maximized at regions corresponding to the vortex cores, consistently with the finding that the material reflectivity is strongly reduced at these zones. In order to compare how these results are altered with increasing temperature, we present in Figure 3 a cross-section of the Casimir force surface at a fixed value of L x = 0.5 , for three different temperatures: T = 2 K, T = 40 K, and T = 90 K, with corresponding lattice size: L x ( 2 K ) = 4.1 nm, L x ( 40 K ) = 4.6 nm, and L x ( 90 K ) = 16.4 nm. We observe that in the low-temperature regime, 2 T 40 K, quite similar periodic patterns arise, essentially independent of the temperature, with a relatively small modulation amplitude Δ F 0.04 pN. On the other hand, for T T c the vortex cell size increases, but the force modulation is drastically reduced.

4. Discussion and Conclusions

The former results have been derived within a mean-field approach that neglects thermal fluctuations of the order parameter and pinning disorder. However, in the weak-interacting limit of Cooper pairs, thermal effects can be taken into account [36]. The effect of disorder in the vortex array can be introduced by adding white noise to the coefficients of the GL free energy and performing the Z functional integration, or by performing vortex matter simulations based on the numerical analysis of the time-dependent GL formulation [35,36].
Although the present study was focused on the action of magnetic fields just below the upper critical field, H H c 2 , the effect of fields near the lower critical value, H H c 1 can be straightforwardly discussed within the clean-limit of the London theory [34]. In that case, the order parameter is given by | ϕ ( r ) | 2 ( 1 + 2 ξ a b 2 / r 2 ) 1 , while the local magnetic induction B ( r ) = Φ 0 K 0 r 2 + 2 ξ a b 2 / λ a b / 2 π λ a b 2 , where K 0 ( x ) is a modified Bessel function. Then, the total order parameter can be built as the product i ϕ ( | r r i | ) , where r i denotes the localization of the different vortexes, whereas the total magnetic induction B ( r ) = i B ( | r r i | ) [61].
The periodic structure of the system could, actually, induce diffractive effects in the predicted spatial-pattern observed in the Casimir force. Indeed it has been shown how diffraction gratings modify the force [62]. We did not consider this kind of effects because they would occur for wavelengths of the order of the period of the Abrikosov lattice, about 4–16 nm in our case. This corresponds to frequencies higher than 300 eV. For these high values of energy, the dielectric function is constant (tends to the high-frequency value, ε ), and no position dependence will be observed.
In conclusion, we presented a general methodology aimed to evaluate the Casimir force in setups constituted by SC materials under the action of an external magnetic field. We have shown that the Abrikosov vortex lattice displayed by a type II superconductor induces Casimir forces with a periodic structure that mirrors the local charge redistribution due to superconducting currents conducing to magnetic fluxon confinement within the vortex cores. This approach may be applied to SC systems under different conditions of temperature, oxygen doping, and magnetic field configurations, allowing the analysis of alternative orderings competing with vortex matter, such as charge density waves [63], or the investigation of normal matter inside the vortexes subject to multiple Andreev reflections [64].

Author Contributions

Conceptualization, S.G.C.-L., R.E.-S., G.P. and C.V.; formal analysis, S.G.C.-L., R.E.-S., G.P. and C.V.; methodology, S.G.C.-L., R.E.-S., G.P. and C.V.; software, S.G.C.-L. and C.V.; writing—original draft, S.G.C.-L., R.E.-S., G.P. and C.V.; writing—review and editing, S.G.C.-L., R.E.-S., G.P. and C.V. All authors have read and agreed to the published version of the manuscript.

Funding

S.G.C-L acknowledges financial support from UNAM DGAPA PAPIIT Grant No. TA100724. G.P. acknowledges financial support from UNAM DGAPA PAPIIT Grant No. IN104522, UNAM DGAPA PAPIME PE101223 and CONACyT Projects No. 1564464 and No. 1098652. R.E.-S. acknowledges partial support from CONACyT Project No. A1-S-10537.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Casimir, H.B.G. On the attraction between two perfectly conducting plates. Proc. Kon. Ned. Akad. Wet. B 1948, 51, 793–795. Available online: https://dwc.knaw.nl/DL/publications/PU00018547.pdf (accessed on 29 January 2024).
  2. Lifshitz, E.M. The theory of molecular attractive forces between solids. Sov. Phys. JETP 1956, 2, 73–83. Available online: http://jetp.ras.ru/cgi-bin/e/index/e/2/1/p73?a=list (accessed on 29 January 2024).
  3. 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]
  4. Harris, B.W.; Chen, F.; Mohideen, U. Precision measurement of the Casimir force using gold surfaces. Phys. Rev. A 2000, 62, 052109. [Google Scholar] [CrossRef]
  5. Bressi, G.; Carugno, G.; Onofrio, R.; Ruoso, G. Measurement of the Casimir force between parallel metallic surfaces. Phys. Rev. Lett. 2002, 88, 041804. [Google Scholar] [CrossRef] [PubMed]
  6. Decca, R.S.; López, D.; Fischbach, E.; Krause, D.E. Measurement of the Casimir force between dissimilar metals. Phys. Rev. Lett. 2003, 91, 050402. [Google Scholar] [CrossRef] [PubMed]
  7. Chen, F.; Klimchitskaya, G.L.; Mohideen, U.; Mostepanenko, V.M. Theory confronts experiment in the Casimir force measurements: Quantification of errors and precision. Phys. Rev. A 2004, 69, 022117. [Google Scholar] [CrossRef]
  8. 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]
  9. Krause, D.E.; Decca, R.S.; López, D.; Fischbach, E. Experimental investigation of the Casimir force beyond the proximity-force approximation. Phys. Rev. Lett. 2007, 98, 050403. [Google Scholar] [CrossRef]
  10. Jourdan, G.; Lambrecht, A.; Comin, F.; Chevrier, J. Quantitative non-contact dynamic Casimir force measurements. EPL (Europhys. Lett.) 2009, 85, 31001. [Google Scholar] [CrossRef]
  11. 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]
  12. 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] [PubMed]
  13. Castillo-Garza, R.; Xu, J.; Klimchitskaya, G.L.; Mostepanenko, V.M.; Mohideen, U. Casimir interaction at liquid nitrogen temperature: Comparison between experiment and theory. Phys. Rev. B 2013, 88, 075402. [Google Scholar] [CrossRef]
  14. 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]
  15. Rodriguez, A.W.; Capasso, F.; Johnson, S.G. The Casimir effect in microstructured geometries. Nat. Phot. 2011, 5, 211–221. [Google Scholar] [CrossRef]
  16. Milonni, P.W. The Quantum Vacuum: An Introduction to Quantum Electrodynamics; Academic Press, Inc.: San Diego, CA, USA, 2013. [Google Scholar] [CrossRef]
  17. Milton, K.A. The Casimir Effect: Physical Manifestations of Zero-Point Energy; World Scientific Co., Ltd.: Singapore, 2001. [Google Scholar] [CrossRef]
  18. Bordag, M.; Klimchitskaya, G.L.; Mohideen, U.; Mostepanenko, V.M. Advances in the Casimir Effect; Oxford University Press Inc.: New York, NY, USA, 2009. [Google Scholar] [CrossRef]
  19. Sushkov, A.; Kim, W.; Dalvit, D.; Lamoreaux, S. Observation of the thermal Casimir force. Nat. Phys. 2011, 7, 230–233. [Google Scholar] [CrossRef]
  20. Simpson, W.M.R.; Leonhardt, U. Forces of the Quantum Vacuum: An Introduction to Casimir Physics; World Scientific: Singapore, 2015. [Google Scholar] [CrossRef]
  21. Bimonte, G.; López, D.; Decca, R.S. Isoelectronic determination of the thermal Casimir force. Phys. Rev. B 2016, 93, 184434. [Google Scholar] [CrossRef]
  22. 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]
  23. 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 mm: Experimental procedures and comparison with theory. Universe 2021, 7, 93. [Google Scholar] [CrossRef]
  24. Bimonte, G.; Calloni, E.; Esposito, G.; Milano, L.; Rosa, L. Towards measuring variations of Casimir energy by a superconducting cavity. Phys. Rev. Lett. 2005, 94, 180402. [Google Scholar] [CrossRef]
  25. Bimonte, G.; Born, D.; Calloni, E.; Esposito, G.; Il’ichev, E.; Rosa, L.; Scaldaferri, O.; Tafuri, F.; Vaglio, R.; Hübner, U. The Aladin2 experiment: Status and perspectives. J. Phys. A Math. Gen. 2006, 39, 6153–6159. [Google Scholar] [CrossRef]
  26. Bimonte, G.; Calloni, E.; Esposito, G.; Rosa, L. Casimir energy and the superconducting phase transition. J. Phys. A Math. Gen. 2006, 39, 6161–6171. [Google Scholar] [CrossRef]
  27. Bimonte, G.; Born, D.; Calloni, E.; Esposito, G.; Huebner, U.; Il’ichev, E.; Rosa, L.; Tafuri, F.; Vaglio, R. Low noise cryogenic system for the measurement of the Casimir energy in rigid cavities. J. Phys. A Math. Theor. 2008, 41, 164023. [Google Scholar] [CrossRef]
  28. Bimonte, G. Casimir effect in a superconducting cavity and the thermal controversy. Phys. Rev. A 2008, 78, 062101. [Google Scholar] [CrossRef]
  29. Norte, R.A.; Forsch, M.; Wallucks, A.; Marinkovi´c, I.; Gröblacher, S. Platform for measurements of the Casimir force between two superconductors. Phys. Rev. Lett. 2018, 121, 030405. [Google Scholar] [CrossRef]
  30. Bimonte, G. Casimir effect between superconductors. Phys. Rev. A 2019, 99, 052507. [Google Scholar] [CrossRef]
  31. Villarreal, C.; Caballero-Benitez, S.F. Casimir forces and high-Tc superconductors. Phys. Rev. A 2019, 100, 042504. [Google Scholar] [CrossRef]
  32. Castillo-López, S.G.; Esquivel-Sirvent, R.; Pirruccio, G.; Villarreal, C. Casimir forces out of thermal equilibrium near a superconducting transition. Sci. Rep. 2022, 12, 2905. [Google Scholar] [CrossRef]
  33. Abrikosov, A.A. Nobel Lecture: Type-II superconductors and the vortex lattice. Rev. Mod. Phys. 2004, 76, 975–979. [Google Scholar] [CrossRef]
  34. Annett, J.F. Superconductivity, Superfluids and Condensates; Oxford University Press Inc.: New York, NY, USA, 2004; Available online: https://archive.org/details/superconductivit0000anne/ (accessed on 29 January 2024).
  35. Blatter, G.; Feigel’man, M.V.; Geshkenbein, V.B.; Larkin, A.I.; Vinokur, V.M. Vortices in hightemperature superconductors. Rev. Mod. Phys. 1994, 66, 1125–1388. [Google Scholar] [CrossRef]
  36. Rosenstein, B.; Li, D. Ginzburg-Landau theory of type II superconductors in magnetic field. Rev. Mod. Phys. 2010, 82, 109–168. [Google Scholar] [CrossRef]
  37. Kwok, W.K.; Welp, U.; Glatz, A.; Koshelev, A.E.; Kihlstrom, K.J.; Crabtree, G.W. Vortices in high-performance high-temperature superconductors. Rep. Prog. Phys. 2016, 79, 116501. [Google Scholar] [CrossRef] [PubMed]
  38. Suderow, H.; Guillamón, I.; Rodrigo, J.G.; Vieira, S. Imaging superconducting vortex cores and lattices with a scanning tunneling microscope. SuST 2014, 27, 063001. [Google Scholar] [CrossRef]
  39. Sonier, J.E.; Brewer, J.H.; Kiefl, R.F. μSR studies of the vortex state in type-II superconductors. Rev. Mod. Phys. 2000, 72, 769–811. [Google Scholar] [CrossRef]
  40. Scheel, S.; Fermani, R.; Hinds, E. Feasibility of studying vortex noise in two-dimensional superconductors with cold atoms. Phys. Rev. A 2007, 75, 064901. [Google Scholar] [CrossRef]
  41. Román-Velázquez, C.E.; Noguez, C.; Villarreal, C.; Esquivel-Sirvent, R. Spectral representation of the nonretarded dispersive force between a sphere and a substrate. Phys. Rev. A 2004, 69, 042109. [Google Scholar] [CrossRef]
  42. Noguez, C.; Román-Velázquez, C.E.; Esquivel-Sirvent, R.; Villarreal, C. High-multipolar effects on the Casimir force: The non-retarded limit. EPL (Europhys. Lett.) 2004, 67, 191–197. [Google Scholar] [CrossRef]
  43. Noguez, C.; Román-Velázquez, C.E. Dispersive force between dissimilar materials: Geometrical effects. Phys. Rev. B 2004, 70, 195412. [Google Scholar] [CrossRef]
  44. Neto, P.A.M.; Lambrecht, A.; Reynaud, S. Casimir energy between a plane and a sphere in electromagnetic vacuum. Phys. Rev. A 2008, 78, 012115. [Google Scholar] [CrossRef]
  45. Canaguier-Durand, A.; Neto, P.A.M.; Cavero-Pelaez, I.; Lambrecht, A.; Reynaud, S. Casimir interaction between plane and spherical metallic surfaces. Phys. Rev. Lett. 2009, 102, 230404. [Google Scholar] [CrossRef]
  46. Canaguier-Durand, A.; Neto, P.A.M.; Lambrecht, A.; Reynaud, S. Thermal Casimir effect in the plane-sphere geometry. Phys. Rev. Lett. 2010, 104, 040403. [Google Scholar] [CrossRef]
  47. Bimonte, G. Going beyond PFA: A precise formula for the sphere-plate Casimir force. EPL (Europhys. Lett.) 2017, 118, 20002. [Google Scholar] [CrossRef]
  48. Tinkham, M. Introduction to Superconductivity; Dover Publications, Inc.: Mineola, NY, USA, 2004. [Google Scholar]
  49. Castillo-López, S.; Villarreal, C.; Esquivel-Sirvent, R.; Pirruccio, G. Enhancing near-field radiative heat transfer by means of superconducting thin films. Int. J. Heat Mass Transf. 2022, 182, 121922. [Google Scholar] [CrossRef]
  50. Nozières, P.; Schmitt-Rink, S. Bose condensation in an attractive fermion gas: From weak to strong coupling superconductivity. J. Low Temp. Phys. 1985, 59, 195–211. [Google Scholar] [CrossRef]
  51. Pal, S.; Ganguly, K.; Basu, A.; Sharma, U.D. Gorter-Casimir two fluid model revisited and possible applications to superconductivity. Int. J. Innov. Res. Phys. 2019, 1, 17–26. [Google Scholar] [CrossRef]
  52. Landau, L.D.; Lifshitz, E.M. Statistical Physics; Elsevier Butterworth-Heinemann: Oxford, UK, 2013. [Google Scholar] [CrossRef]
  53. Fujita, S.; Godoy, S. Quantum Statistical Theory of Superconductivity; Springer Science + Business Media: New York, NY, USA, 1996. [Google Scholar] [CrossRef]
  54. Lomnitz, M.; Villarreal, C.; de Llano, M. BEC model of high-Tc superconductivity in layered cuprates. Int. J. Mod. Phys. B 2013, 27, 1347001. [Google Scholar] [CrossRef]
  55. Zuev, Y.; Kim, M.S.; Lemberger, T.R. Correlation between superfluid density and Tc of underdoped YBa2 Cu3 O6+x near the superconductor-insulator transition. Phys. Rev. Lett. 2005, 95, 137002. [Google Scholar] [CrossRef] [PubMed]
  56. Chen, Q.; Stajic, J.; Tan, S.; Levin, K. BCS–BEC crossover: From high temperature superconductors to ultracold superfluids. Phys. Rep. 2005, 412, 1–88. [Google Scholar] [CrossRef]
  57. Bonn, D.A.; O’Reilly, A.H.; Greedan, J.E.; Stager, C.V.; Timusk, T.; Kamarás, K.; Tanner, D.B. Far-infrared properties of ab-plane oriented YBa2Cu3O7-δ. Phys. Rev. B 1988, 37, 1574–1579. [Google Scholar] [CrossRef] [PubMed]
  58. Timusk, T.; Herr, S.L.; Kamarás, K.; Porter, C.D.; Tanner, D.B.; Bonn, D.A.; Garrett, J.D.; Stager, C.V.; Greedan, J.E.; Reedyk, M. Infrared studies of ab-plane oriented oxide superconductors. Phys. Rev. B 1988, 38, 6683–6688. [Google Scholar] [CrossRef] [PubMed]
  59. Basov, D.N.; Timusk, T. Electrodynamics of high-Tc superconductors. Rev. Mod. Phys. 2005, 77, 721–779. [Google Scholar] [CrossRef]
  60. Klimchitskaya, G.L.; Mohideen, U.; Mostepanenko, V.M. Kramers–Kronig relations for plasma-like permittivities and the Casimir force. J. Phys. A Math. Theor. 2007, 40, F339–F346. [Google Scholar] [CrossRef]
  61. Brandt, E.H. The flux-line lattice in superconductors. Rep. Prog. Phys. 1995, 58, 1465–1594. [Google Scholar] [CrossRef]
  62. Intravaia, F.; Koev, S.; Jung, I.W.; Talin, A.A.; Davids, P.S.; Decca, R.S.; Aksyuk, V.A.; Dalvit, D.A.; López, D. Strong Casimir force reduction through metallic surface nanostructuring. Nat. Commun. 2013, 4, 2515. [Google Scholar] [CrossRef]
  63. Chang, J.; Blackburn, E.; Holmes, A.T.; Christensen, N.B.; Larsen, J.; Mesot, J.; Liang, R.; Bonn, D.A.; Hardy, W.N.; Watenphul, A.; et al. Direct observation of competition between superconductivity and charge density wave order in YBa2Cu3O6.67. Nat. Phys. 2012, 8, 871–876. [Google Scholar] [CrossRef]
  64. Huebener, R.P. The Abrikosov vortex lattice: Its discovery and impact. J. Supercond. Nov. Magn. 2019, 32, 475–481. [Google Scholar] [CrossRef]
Figure 1. Setup consisting of a spherical Au nanoparticle of radius R located at a minimum distance d from a planar substrate made of optimally-doped YBCO ( YBa 2 Cu 3 O 7 δ ), both at the same temperature T T ns and T sub , respectively,—in presence of an applied magnetic field directed along the z axis. We show the resulting lattice structure with elementary cells of nondimensional spatial periods a = L x = L y . Induced supercurrents coincide with equiprobability contours defined by the number density of Cooper pairs, n s = | Ψ ( r ˜ ) | 2 , where Ψ ( r ˜ ) is the superconducting order parameter. The vortex cores of radius ξ correspond to the darkest inner zones. B denotes the magnetic field, and the tilde stands for adimensional variables, x ˜ = x / λ a b and y ˜ = y / λ a b , where λ a b repersents the a b -plane penetration length.
Figure 1. Setup consisting of a spherical Au nanoparticle of radius R located at a minimum distance d from a planar substrate made of optimally-doped YBCO ( YBa 2 Cu 3 O 7 δ ), both at the same temperature T T ns and T sub , respectively,—in presence of an applied magnetic field directed along the z axis. We show the resulting lattice structure with elementary cells of nondimensional spatial periods a = L x = L y . Induced supercurrents coincide with equiprobability contours defined by the number density of Cooper pairs, n s = | Ψ ( r ˜ ) | 2 , where Ψ ( r ˜ ) is the superconducting order parameter. The vortex cores of radius ξ correspond to the darkest inner zones. B denotes the magnetic field, and the tilde stands for adimensional variables, x ˜ = x / λ a b and y ˜ = y / λ a b , where λ a b repersents the a b -plane penetration length.
Physics 06 00026 g001
Figure 2. Periodic structure of the Casimir force as a function of the location of the Au nanosphere over the Abrikosov lattice at T = 2 K, for a fixed distance d = 2 R . Here, R = L x ( 2 K ) 4 nm. It can be observed that the minimal strength of the Casimir force corresponds to the vortex cores.
Figure 2. Periodic structure of the Casimir force as a function of the location of the Au nanosphere over the Abrikosov lattice at T = 2 K, for a fixed distance d = 2 R . Here, R = L x ( 2 K ) 4 nm. It can be observed that the minimal strength of the Casimir force corresponds to the vortex cores.
Physics 06 00026 g002
Figure 3. Comparison of the Casimir force profiles as a function of the normalized coordinate y ˜ / L y along a line passing right above the vortex core ( x ˜ / L x = 0.5 ), and d = 2 R = 8 nm, at three different temperatures: T = 2 K, T = 40 K and T = 90 K as indicated. Depending on the temperature, each unit period L y corresponds to spatial periods: L y ( 2 K ) = 4.1 nm, L y ( 40 K ) = 4.6 nm, and L y ( 90 K ) = 16.4 nm. We observe that in the low-temperature regime the force magnitude shows almost coincident values up to T = 40 K, consistent with expectations that vacuum fluctuations ( c / d ) overwhelm thermal fluctuations ( k B T ) at nanometer separations. Here, , c, and k B denote the reduced Planck constant, the speed of light, and the Boltzmann constant, respectively.
Figure 3. Comparison of the Casimir force profiles as a function of the normalized coordinate y ˜ / L y along a line passing right above the vortex core ( x ˜ / L x = 0.5 ), and d = 2 R = 8 nm, at three different temperatures: T = 2 K, T = 40 K and T = 90 K as indicated. Depending on the temperature, each unit period L y corresponds to spatial periods: L y ( 2 K ) = 4.1 nm, L y ( 40 K ) = 4.6 nm, and L y ( 90 K ) = 16.4 nm. We observe that in the low-temperature regime the force magnitude shows almost coincident values up to T = 40 K, consistent with expectations that vacuum fluctuations ( c / d ) overwhelm thermal fluctuations ( k B T ) at nanometer separations. Here, , c, and k B denote the reduced Planck constant, the speed of light, and the Boltzmann constant, respectively.
Physics 06 00026 g003
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Castillo-López, S.G.; Esquivel-Sirvent, R.; Pirruccio, G.; Villarreal, C. Casimir Forces with Periodic Structures: Abrikosov Flux Lattices. Physics 2024, 6, 394-406. https://doi.org/10.3390/physics6010026

AMA Style

Castillo-López SG, Esquivel-Sirvent R, Pirruccio G, Villarreal C. Casimir Forces with Periodic Structures: Abrikosov Flux Lattices. Physics. 2024; 6(1):394-406. https://doi.org/10.3390/physics6010026

Chicago/Turabian Style

Castillo-López, Shunashi Guadalupe, Raúl Esquivel-Sirvent, Giuseppe Pirruccio, and Carlos Villarreal. 2024. "Casimir Forces with Periodic Structures: Abrikosov Flux Lattices" Physics 6, no. 1: 394-406. https://doi.org/10.3390/physics6010026

APA Style

Castillo-López, S. G., Esquivel-Sirvent, R., Pirruccio, G., & Villarreal, C. (2024). Casimir Forces with Periodic Structures: Abrikosov Flux Lattices. Physics, 6(1), 394-406. https://doi.org/10.3390/physics6010026

Article Metrics

Back to TopTop