Casimir Forces with Periodic Structures: Abrikosov Flux Lattices

: 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.


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 = −hcπ 2 /240d 4  [1] with h = 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 = τ −1 el−imp + τ −1 el−el (T) + τ −1 el−ph (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 ∼ hc/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 ≪ hc/d for sub-µm body separations.Therefore, indirect approaches have been proposed based on observations of the Casimirinduced 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 c1 < H c2 .For values of the applied field higher than H c1 , magnetic flux lines penetrate the sample in the form of quantum vortexes, Φ 0 = h/2e, 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 c2 , where superconductivity disappears.Remarkably, as shown by Alexei Abrikosov [33], for intensities of the applied field just below H c2 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], , 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.

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 finitetemperature 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: 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: 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 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: 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 Here, T is the dipole-dipole interaction tensor T = (3r 0 r 0 − r 2 0 I)/r 5 0 , 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: which, by introducing the function u(ω) = [1 − ϵ ns (ω))] −1 and explicitly substituting α(ω), may be re-expressed as the secular equation: where and 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 which in the present case implies three eigenvalues, (9)   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: 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 Tn/h, with n an integer number.In that case, the final expression for temperature-dependent Casimir force is: 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 . 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.

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 Au /(ω 2 + iγ Au ω), where the gold plasma frequency, ω Au = 9.1 eV, and the inverse scattering rate, γ Au = 0.02 eV.Therefore, 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 ab-planes), whereas CP transport along the perpendicular c-axis occurs due to Josephson tunneling, with respective effective masses, m * ab and m * c , reflecting an anisotropic charge transport measured by the ratio γ 2 = m * c /m * ab ≈ 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(ε ab , ε ab , ε 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 = ε ab .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 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 π = −ih∇ + e * A, with e * = 2e the CP charge, and the reciprocal mass tensor is also diagonal in the chosen coordinate system: (1/m * ) = diag(1/m * ab , 1/m * ab , 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 while the differentiation δF s (T)/δA(r) yields the current density It can be shown that the mass anisotropy induces in turn anisotropic coherence lengths, ξ j (T), satisfying the relations ξ 2 ab (T) = h2 /2m * ab |a T |, and ξ 2 c (T) = h2 /2m * c |a T | [48].In absence of external fields or material boundaries, the former equations predict a secondorder 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) with H c2 (T) = Φ 0 /2πξ 2 ab (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 c2 .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 c2 , 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 = Bx 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 As shown by Abrikosov [33], Equation ( 16) admits Landau-level solutions of the form Ψ(r) = Φ(r ⊥ )e ik z z , where the ground state is given by k z = 0, and 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 2inz , so that In this case, κ = λ ab /ξ ab where λ ab is the ab-plane penetration length, while x = x/λ ab , ỹ = y/λ ab .Notice that λ ab describes the magnetic field screening by supercurrents flowing along the ab-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 = (he * /2m * ab )|Φ| 2 (∇χ − (e * /h)A), indicating that supercurrent 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

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 weaklyinteracting particles with BE statistical properties.As a consequence, the energy excitation spectrum is given by Bogoliubov expression which, in the low-momentum limit, leads to a phonon dispersion relation, E k ≈ hc 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] Here, , 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 2D 0 (T)/n 2D = 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 λ 2 ab (0)/λ 2 ab (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 2D 0 (r, T) ≈ n 2D s (r, T), and then

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, ω 2 p = e * 2 n/ε 0 m * , satisfies ω p (2K) = ω 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: 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: 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ζ) = ε Au (iζ), ab (iζ, r ⊥ , T < T c ).

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 = 2R.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.

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 weakinteracting 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 c2 , the effect of fields near the lower critical value, H ≥ H c1 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ξ 2 ab /r 2 ) −1 , while the local magnetic induction B(r) = Φ 0 K 0 r 2 + 2ξ 2 ab /λ ab /2πλ 2 ab , 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].

Figure 2 .
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 = 2R.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 3 .
Figure 3.Comparison of the Casimir force profiles as a function of the normalized coordinate ỹ/L y along a line passing right above the vortex core ( x/L x = 0.5), and d = 2R = 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 (∼ hc/d) overwhelm thermal fluctuations (∼ k B T) at nanometer separations.Here, h, c, and k B denote the reduced Planck constant, the speed of light, and the Boltzmann constant, respectively.