Next Article in Journal
Mean Remaining Strength Estimation of Multi-State System Based on Nonparametric Bayesian Method
Next Article in Special Issue
Dependence Relations and Grade Fuzzy Set
Previous Article in Journal
Homotopy Coherent Representations
Previous Article in Special Issue
An Optimization Model of Integrated AGVs Scheduling and Container Storage Problems for Automated Container Terminal Considering Uncertainty
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Superconductors and Gravity

by
Antonio Gallerati
1,2,* and
Giovanni Alberto Ummarino
1,3
1
Dipartimento di Scienza Applicata e Tecnologia, Politecnico di Torino, Corso Duca degli Abruzzi 24, 10129 Torino, Italy
2
Istituto Nazionale di Fisica Nucleare, Sezione di Torino, Via Pietro Giuria 1, 10125 Torino, Italy
3
Institute for Physics and Engineering, National Research Nuclear University MEPhI, Kashirskoe hwy 31, 115409 Moscow, Russia
*
Author to whom correspondence should be addressed.
Symmetry 2022, 14(3), 554; https://doi.org/10.3390/sym14030554
Submission received: 20 January 2022 / Revised: 18 February 2022 / Accepted: 8 March 2022 / Published: 10 March 2022

Abstract

:
We review and discuss some recent developments on the unconventional interaction between superconducting systems and the local gravitational field. While it is known that gravitational perturbations (such as gravitational waves) can affect supercondensates and supercurrents dynamics, we want to focus here on the more subtle superfluid back-reaction acting on the surrounding gravitational field, analysing some specific favourable situations. To this end, we will consider suitable quantum macrosystems in a coherent state, immersed in the static weak Earth’s gravitational field, investigating possible slight local alterations of the latter not explained in terms of classical physics.

1. Introduction

The gravitational force has the distinctive feature of universal interaction with all forms of matter and energy. It dominates at large-scales where it is well described by general relativity. In the latter theory, gravity is not interpreted as a standard force acting on different masses, but as a direct affection of the geometry of the spacetime: masses generate curvature, which in turn dictates the motion of the masses. The spacetime then plays a dynamical role and it is not a rigid background structure.
While classical general relativity gives a consistent description of the large-scale dynamics dominated by gravity, we know that quantum field theory is the fundamental formulation to describe physics at the microscopic scale, where the effects gravity are in general negligible. In the last few decades, different quantum gravity formulations have been proposed to consistently describe the physics of the particles when the gravitational field is so intense as to affect the motion of elementary particles (presumably in the vicinity of a black hole or a neutron star, as well as in the early stages of the evolution of our Universe). This is clearly an ambitious target, since it will imply a fundamental knowledge about the functioning of the laws of nature. However, experimental verification of this kind of theories is really hard to realize, since this would in general imply very high ranges of energy. Then, direct observation of quantum gravity effects, involving gravitons’ dynamical interactions with other quantum fields at the microscopic level, is a very difficult task.
A different approach could originate from the study of unconventional, macroscopic states of matter. In this regard, one should consider quantum macrosystems existing in nature, like superconductors and superfluids. The latter can be thought as large systems featuring a macroscopical coherent phase, suitably described by order parameters. It could be then possible to formulate (and observe) a possible interplay between the extended, coherent system and the surrounding gravitational field [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21]. In this regard, the coupling with the current flow without resistance in superconductors was exploited to use the latter as a sensitive detection systems, in particular for gravitational waves [22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38]. Another remarkable phenomenon, showing quantum effects originating from the interaction of quantum particles with a weak-field gravitational background, is the gravity-induced quantum interference [39,40,41,42,43,44,45,46,47,48]. This effect takes place in the presence of a gravitational potential, to be considered in the Schrodinger equation [49,50] and giving rise to a phase shift for elementary particles (The experimental effect can be measured splitting a nearly monoenergetic beam of thermal particles and considering the produced interference paths: a gravity-induced quantum mechanical phase shift is observed, due to the presence of the Earth’s gravitational field [39,40]).
Inspired by the above results, we then also want to consider the back-reaction of superfluids and supercurrents on the local gravitational field in some specific, favourable situation. The first step to achieve the goal will be to formulate an appropriate theoretical model justifying this anomalous coupling. In the following subsections, we will briefly discuss the most convincing theoretical basis and experimental evidence in favour of the existence of this unconventional interaction.

1.1. Theoretical Foundations

We now want to characterize a possible interplay between superfluids and the local gravitational field in the framework of a quantum gravity theory or, at least, in a suitable approximation of the latter for weak fields.
Let us first consider the classical picture. Clearly, the absence of (gravitational) charges of opposite sign excludes the possibility of counteracting the field inside the medium by a local redistribution, ruling out dielectric-type effects. If we then take the medium to be a standard quantum mechanical system, the smallness of the gravitational coupling strongly suppresses the possibility of a (graviton) excitation for a medium particle and any subsequent affection of the local field. We are then led to consider the interaction of the gravitational field with an anomalous external source, that is, an unconventional state of matter exhibiting quantization on a macroscopic scale, like a Bose condensate or a more generic superfluid.
Let us then consider a quantum gravity framework and write the Lagrangian for this coherent macrosystem; immersed in the Earth’s gravitational field, (we work in the “mostly plus” convention, where the Minkowski metric is η μ v = diag ( 1 , + 1 , + 1 , + 1 ) , and set c = = 1 ) [8,9]:
L = L E H + L ϕ = 1 8 π G R 2 Λ 1 2 g μ v μ ϕ v ϕ + 1 2 m 2 ϕ ϕ .
The first term is the Einstein–Hilbert contribution, R being the Ricci scalar and Λ the cosmological constant. The other terms describe the dynamics of the medium supercondensate (for example, Cooper pairs of mass m) that we can characterize as a bosonic field ϕ with non-vanishing vacuum expectation value (vev) ϕ 0 = 0 | ϕ | 0 . We assume this vev to be forced from the outside to a certain value, as it happens, for instance, in a superconductor subjected to external electromagnetic fields. In the weak gravity limit, the metric g μ v can be written as
g μ v ( x ) = η μ v + h μ v ( x ) ,
sum of the flat Minkowski background η μ v plus perturbations given by the h μ v ( x ) contribution. Now, we expand the bosonic field as
ϕ ( x ) = ϕ 0 ( x ) + ϕ ¯ ( x ) ,
where the ϕ 0 vev depends on the medium characteristics and can be seen as an external source, while the ϕ ¯ contribution is included in the integration variables. The scalar field ϕ then suitably describes a superfluid with ground state density ϕ 0 fixed by external conditions (for example, Cooper pairs density in a supercondensate in the presence of external EM fields). In the weak gravity limit, the L ϕ term then reads
L ϕ = L ϕ + L h + L 0 .
In the above expansion, the first term is related to the ϕ ¯ contributions, involved in the negligible excitation processes related to the graviton emission–absorption mechanism, and several vertices of interaction that turn out to be irrelevant due to the smallness of the gravitational coupling. The second term takes into account the coupling of the condensate with the h μ v metric fluctuations and is written as
L h h μ v μ ϕ 0 v ϕ 0 ,
which determines corrections to the gravitational propagator, which is again a negligible contribution. Finally, the last term determines a local supercondensate contribution to the total effective cosmological term of the form
L 0 = 1 2 μ ϕ 0 μ ϕ 0 + 1 2 m 2 ϕ 0 2 ,
connected to the coherent vacuum energy density and depending on the fixed external source ϕ 0 . The above coupling has the correct structure to produce possible, localized instabilities in superfluid regions featuring larger condensate density [8,9]: this could determine detectable effects in spite of the smallness of the gravitational coupling. We should also note that, in the latter instable regions, some physical cutoff or regularizing process should come into play, preventing local contribution of arbitrary intensity (This can be considered a gravitational analog of the Casimir effect, where observable evidence originates from inhomogeneities in the vacuum fluctuations. In the latter case, the metallic conductors impose a cutoff on the electromagnetic vacuum fluctuations, while the same role is played here by the coherent superfluid). The field then tends to be pinned, assuming fixed extremal values which are independent from those in the neighbouring regions. One could expect, as a physical effect, some kind of slight partial shielding (“absorption”) locally affecting fields propagation and potentials. As we already pointed out, the introduced superfluid density ϕ 0 ( x ) is related to the microscopic structure of the involved sample, as well as to the presence of currents, vortex lattices and electromagnetic fields in the supercondensate (It has also been conjectured that high-frequency electromagnetic fields could provide the required energy to enhance the described gravitational field affection [8,15]).
We have then described a theoretical quantum gravity model with an unconventional coupling between the local gravitational field and the superfluid. The existence of strong variations of the supercondensate components density (for example, Cooper pairs) produce small regions with higher density, where a criticality condition could take place, giving rise to localized instabilities. This gives us a possible way to elude the weakness of the standard coupling and produce a related affection of the local gravitational field. The key ingredient is the macroscopic quantum coherence of the condensate that is taken into account when computing the anomalous interplay, at a fundamental level, between the superfluid and the external gravitational field.

1.2. Experimental Evidence

The discussed formulation laid the foundations of a theoretical approach to an unconventional coupling between superfluids and gravity, in the framework of a quantum model. This, however, involves a formalism that makes it almost impossible to extract quantitative predictions, to be tested in a laboratory experiment. For this reason, one is then led to also consider many phenomenological research reports and evidence, to better understand the proposed interplay and obtain an effective theory leading to more explicit experimental predictions.
One the first attempts to formulate an effective quantum model describing the interaction between conductors and the local gravitational field was given in [51], where a quantum-mechanical formalism is developed to calculate an (additional) electric field component, generated in the vicinity of a conductor by the presence of the Earth’s gravity. The main consequence of this formulation is the definition of generalized electric-type fields and potentials, existing near the surface of a conductor and featuring a gravitationally-induced component. We can schematically express this generalized field as
E = E e + E ind V = V e + V ind
where E e is the standard electric field, while E ind is the gravitationally-induced component. In [52,53], the induced E ind and V ind were experimentally detected as a direct affection to the free fall of electrons in the presence of conductors. This evidence of existence of generalized fields and potentials was then theoretically analysed and experimentally verified in subsequent works [54,55,56,57,58,59,60,61].
Analogous concepts were subsequently extended to superconductors, obtaining similar results [11,33,62,63,64,65,66,67,68,69,70,71,72,73]: generalized gravitoelectric and gravitomagnetic fields can be induced by the presence of a local gravitational field coupled to the supercondensate.
In the following section, we discuss a formal derivation of a consistent form for these generalized fields and potentials, exploiting a weak field expansion for the local gravitational field. This approach will lead us to the definition of a generalized form for Maxwell equations.

2. Linearized Gravity: Gravito–Maxwell Fields

It is well known that gravity is in general mediated by a symmetric g μ v tensor field, featuring 10 independent components (potentials). However, under certain approximations and suitable gauge choice, the gravitational field behaviour can be described in an electromagnetic-like fashion, by means of vector-like field equations instead of the corresponding tensorial expressions. In particular, linearized gravity can be considered as a consistent weak-field limit of the complete tensorial theory, in the regime where nonlinear effects can be ignored (In the linear-order assumption, the gravitational field does not transfer energy to the gravitational sources; this also cancels matter–gravity coupling from the domain of linear approximations).
Let us the consider a nearly-flat spacetime, characterized by the presence of a weak and static gravitational field. This means we can consider small perturbation of the Minkowski metric η μ v and express the spacetime metric g μ v as
g μ v η μ v + h μ v ,
where the symmetric tensor h μ v is a small perturbation of the constant flat η μ v in the mostly plus convention, η μ v = diag ( 1 , + 1 , + 1 , + 1 ) . The inverse metric, in linear approximation, is given by
g μ v η μ v h μ v .
while the metric determinant can be expanded as
g = det [ g μ v ] = ε μ v ρ σ g 1 μ g 2 v g 3 ρ g 4 σ 1 h g 1 + 1 2 h ,
where h = h σ σ .
We are now going to exploit the above weak-limit expansion of the metric to obtain a linearized form for the Einstein and London equations. We will then take advantage of the obtained results to define suitable backgrounds to test the discussed gravity/superfluid interplay.

2.1. Generalizing Maxwell Equations

Let us put ourselves in an inertial coordinate system. To first order in h μ v , the connection is expanded as
Γ λ μ v 1 2 η λ ρ μ h v ρ + v h ρ μ ρ h μ v .
The Riemann tensor is defined as:
R σ μ λ v = λ Γ σ μ v v Γ σ μ λ + Γ σ ρ λ Γ ρ v μ Γ σ ρ v Γ ρ λ μ ,
while the Ricci tensor is obtained from the contraction
R μ v = R σ μ σ v ,
and, to linear order in h μ v , it is expressed as
R μ v σ Γ σ μ v + μ Γ σ σ v + Γ Γ Γ Γ = 1 2 μ ρ h v ρ + v ρ h μ ρ 1 2 ρ ρ h μ v 1 2 μ v h = ρ ( μ h v ) ρ 1 2 2 h μ v 1 2 μ v h ,
having used Equation (11).
The Einstein equations are written as:
R μ v 1 2 g μ v R = 8 π G T μ v ,
where R = g μ v R μ v is the Ricci scalar. In linear-order approximation, we have
1 2 g μ v R 1 2 η μ v η ρ σ R ρ σ = 1 2 η μ v ρ σ h ρ σ 2 h ,
having used Equation (14). The l.h.s. of (15) then reads
R μ v 1 2 g μ v R ρ ( μ h v ) ρ 1 2 2 h μ v 1 2 μ v h 1 2 η μ v ρ σ h ρ σ 2 h .
Let us now introduce the symmetric traceless tensor
h ¯ μ v = h μ v 1 2 η μ v h ,
so that (17) is rewritten as
R μ v 1 2 g μ v R 1 2 ρ μ h ¯ v ρ + ρ v h ¯ μ ρ ρ ρ h ¯ μ v η μ v ρ σ h ¯ ρ σ = ρ [ v h ¯ ρ ] μ + ρ σ η μ [ σ h ¯ v ] ρ = ρ [ v h ¯ ρ ] μ + σ η μ [ ρ h ¯ v ] σ .
We also define the tensor
G μ v ρ [ v h ¯ ρ ] μ + σ η μ [ ρ h ¯ v ] σ ,
in terms of which the Einstein equations take the compact form:
ρ G μ v ρ = 8 π G T μ v .

2.1.1. Gauge Fixing

We now consider the harmonic coordinate condition, expressed by the relation [74,75]:
μ g g μ v = 0 x μ = 0 ,
that in turn can be rewritten in the form
g μ v Γ λ μ v = 0 ,
also known as De Donder gauge. The requirement of the above coordinate condition (22) then plays the role of gauge fixing. In particular, in harmonic coordinates, the metric satisfies a manifestly Lorenz-covariant condition, so that the De Donder gauge becomes a natural choice. Moreover, if one considers the weak-field expansion of the Einstein–Hilbert action in De Donder gauge, the action itself (as well as the graviton propagator) takes a particularly simple form.
Using Equations (8) and (11) together with the above gauge fixing (23), in first-order approximation, we find:
0 1 2 η μ v η λ ρ μ h v ρ + v h ρ μ ρ h μ v = μ h μ λ 1 2 λ h ,
that in turn implies the condition
μ h μ v 1 2 v h μ h μ v 1 2 v h ,
We can also write
μ h μ v = μ h ¯ μ v + 1 2 η μ v h = μ h ¯ μ v + 1 2 v h ,
so that, using Equation (25), we obtain the Lorentz gauge condition:
μ h ¯ μ v 0 .
This condition further simplifies Equation (20) for G μ v ρ , which takes the simple form
G μ v ρ [ v h ¯ ρ ] μ ,
and satisfies the relation
[ λ G 0 μ v ] = 0 S 0 μ v μ A v v A μ ,
The above expression then suggests the existence of a potential. In the following paragraphs, we are going to analyse in detail suitable expressions for fields and potentials defining the desired formalism.

2.1.2. Gravito–Maxwell Equations

Let us now define the following fields [19,63,68,76,77] (for the sake of simplicity, we initially set the physical charge e = m = 1 )
E g E i = 1 2 G 00 i = 1 2 [ 0 h ¯ i ] 0 ,
A g A i = 1 4 h ¯ 0 i ,
B g B i = 1 4 ε i j k G 0 j k ,
with i = 1 , 2 , 3 and
G 0 i j = [ i h ¯ j ] 0 = 1 2 i h ¯ j 0 j h ¯ i 0 = 4 [ i A j ] .
From the above definitions, it follows that
B g = 1 4 ε i j k 4 [ j A k ] = ε i j k j A k = × A g ,
that also implies
· B g = 0 .
We then also find
· E g = i E i = i G 00 i 2 = 8 π G T 00 2 = 4 π G ρ g ,
having used Equation (21) and defined the mass density as ρ g T 00 .
We then consider the curl of E g :
× E g = ε i j k j E k = ε i j k j G 00 k 2 = 1 2 ε i j k j [ 0 h ¯ k ] 0 = 1 4 4 0 ε i j k j A k = 0 B i = B g t .
Finally, for the curl of B g , we find
× B g = ε i j k j B k = 1 4 ε i j k ε k m j G 0 m = 1 4 δ i δ j m δ i m δ j j G 0 m = 1 2 j G 0 i j = 1 2 μ G 0 i μ + 0 G 0 i 0 = 1 2 μ G 0 i μ 0 G 00 i = 1 2 8 π G T 0 i 0 G 00 i = 4 π G j i + E i t = 4 π G j g + E g t ,
having used again Equation (21) and defined the mass gravito-current density vector as j g j i T 0 i .
In summary, the fields (30) are defined, one can write the field equations [11,19,28,32,63,67,68,76,78,79,80,81,82,83,84,85,86,87,88,89,90,91]:
· E g = 4 π G ρ g , · B g = 0 , × E g = B g t , × B g = 4 π G c 2 j g + 1 c 2 E g t ,
having restored physical units. The above expressions are formally equivalent to Maxwell equations, with E g and B g being the gravitoelectric and gravitomagnetic field, respectively (For instance, on the Earth’s surface, E g corresponds to the Newtonian gravitational acceleration, while B g is related to angular momentum interactions [63,68,78,81]). The mass current density vector j g can also be written as:
j g = ρ g v ,
in terms of the mass density and velocity v .

2.1.3. Generalized Maxwell Equations

Inspired by the discussion of Section 1.2, it is now straightforward to extend the above results and introduce generalized electric/magnetic fields, scalar and vector potentials. The latter feature both electromagnetic and gravitational contributions and can be written as:
E = E e + m e E g B = B e + m e B g V = V e + m e V g A = A e + m e A g
where m and e are the electron mass and charge, respectively [51].
The generalized Maxwell equations then become:
· E = 1 ε g + 1 ε 0 ρ , · B = 0 , × E = B t , × B = μ g + μ 0 j + 1 c 2 E t ,
where ε 0 and μ 0 are the standard electric permittivity and magnetic permeability in the vacuum, and where we have set
ρ g = m e ρ , j g = m e j ,
ρ and j being the electric charge density and electric current density, respectively. The introduced vacuum gravitational permittivity ε g and vacuum gravitational permeability μ g are defined as
ε g = 1 4 π G e 2 m 2 , μ g = 4 π G c 2 m 2 e 2 .
The obtained generalized Maxwell equations turn out to be a consistent approximation to the complete tensorial theory, valid in the limit of the weak gravitational field (like the static, weak Earth’s gravity). It is then possible to take advantage of the obtained results and consider suitable situations and parameters regime where the gravitoelectric field plays a fundamental role and/or where gravitomagnetic effects are not negligible.

2.2. Generalizing London Equations

The London equations for a superconductor in stationary state characterize an analogous Ohm’s law (zero resistivity) and Meissner effect (expulsion of the magnetic field from the interior sample) for the superfluid. They can be explicitly written as [92,93,94]:
E e = m n s e 2 j t ;
B e = m n s e 2 × j .
where j = n s e v s is the supercurrent and n s is the superelectron density.
The Ampère’s law for a superconductor in stationary state (no displacement current) has the form
× B e = μ 0 j ,
so that taking the curl gives
× × B e = · B e 2 B e = μ 0 × j = μ 0 n s e 2 m B e ,
that is,
2 B e = 1 λ e 2 B e ,
having introduced the penetration depth
λ e = m μ 0 n s e 2 .
The above quantity gives an estimate of the mean distance the magnetic field B e can penetrate the sample. Since the values of the λ e parameter vary from 2 nm (low T c superconductors) to 200 nm (high T c superconductors), the above Equations (47) and (48) quantitatively characterize the Meissner effect.
The two London Equations (43) and (44) can be now rewritten in terms of the vector potential A e in the (not gauge-invariant) form:
j = 1 μ 0 λ e 2 A e
with B e = × A e and expressing the electric field as E e = A e t .

Generalized London Equations

Let us now take into account gravitational contributions, and consider for the fields and potentials the generalized form (39). In particular, we consider the generalized potential A minimally coupled to the wave function
ψ = ψ 0 exp ( i φ ) ψ 0 2 | ψ | 2 = n s .
The second London equation can be derived from the quantum mechanical current density
j = i 2 m ψ ˜ ψ ψ ˜ ψ ,
where ˜ is the covariant derivative for the minimal coupling:
˜ = i g ˜ A ,
with unknown coupling constant g ˜ . We then find for the current
j = i 2 m ψ ψ ψ ψ g ˜ m A | ψ | 2 = 1 m | ψ | 2 φ g ˜ A .
Now, taking the curl of the previous expression gives
B = m g ˜ | ψ | 2 × j = 1 ζ × j ,
which is the generalized form of the second London Equation (44) [79].
We now want to fix the values of the ζ parameter and coupling constant g ˜ . To this end, let us restrict to the case B g = 0 :
B = B e + m e B g = 1 ζ × j ,
so that, using (44), (48) and (50), we find
g ˜ = e 2 1 ζ = μ 0 λ e 2 .
In order to define an analogue gravitational penetration depth, we now consider the case B e = 0 :
B = B e + m e B g = μ 0 λ e 2 × j = μ 0 λ e 2 m e × j g ,
the gravito-Ampère’s law (37) in stationary state reading
× B g = μ g j g .
Taking the curl of the above equation, we have
× × B g = 2 B g = μ g × j g = μ g 1 μ 0 λ e 2 B g = 1 λ g 2 B g ,
having introduced the gravitational penetration depth
λ g = μ 0 λ e 2 μ g = c 2 4 π G m n s .
Writing now the stationary generalized Ampère’s law (40) and using Equation (60), we obtain
× B = μ 0 + μ g j = μ 0 1 + λ e 2 λ g 2 j ,
and, taking the curl, we find the general form
2 B = μ 0 1 + λ e 2 λ g 2 × j = μ 0 1 μ 0 λ e 2 1 + λ e 2 λ g 2 B = = 1 λ e 2 + 1 λ g 2 B = 1 λ 2 B ,
where we have introduced the generalized penetration depth  λ :
λ = λ g λ e λ g 2 + λ e 2 λ e , with λ g λ e 10 21 .
Finally, we can recast Equation (49) in the form
j = ζ A ,
with B = × A . Moreover, since charge-conservation requires the condition · j = 0 , we obtain for the vector potential
· A = 0 ,
that is, the so-called Coulomb gauge (or London gauge).
In the following sections, we are going to consider suitable frameworks where the proposed gravity/superfluid interplay can in principle be detected, precisely characterizing the physical system and optimizing the range of parameters in order to maximize the effect. We will also exploit the described formalism and introduce generalized fields.

3. A Simple Application: Josephson Effect

The Josephson effect consists of the transmission of supercurrents through thin insulating barriers by means of quantum-mechanical tunnelling [95]. The phenomenon can be seen as a general property of coupled superconducting systems and could take place in suitable tunnel junctions, where quantum interference appears. In particular, when the states of two superconductors are assumed to be coherent (that is, coherent superpositions of states with different numbers of particle pairs), there exists a phase-dependent coupling energy between the two. The latter then implies the possibility of a supercurrent flowing across the junction [96].

3.1. Josephson Junction

If two superconductors are put in contact and the critical current in the contact region is well below those of the individual constituents, the configuration is defined as weak link. Once the weak link is formed, coherence is established across the barrier, with a phase difference Δ φ causing interference between the previously independent wavefunctions, so that the system can be described with a single wavefunction as a whole.
When two superconducting samples are connected through a weak Josephson link, the response of the supercondensate (through the corresponding coupling energy) keeps the macroscopic internal coherence of the system, allowing direct observation of coherence–interference phenomena. In particular, a simple manifestation of the Josephson effect can be observed in a circuit closed on a superconductor–insulator–superconductor (SIS) junction, to which a constant potential difference Δ V is applied. The voltage, in turn, produces a sinusoidal superconductive current across the junction with pulsation [95,96,97]
ω = 2 e Δ V .
Let us briefly discuss the phenomenon.

Josephson AC Current

Let us consider a weak link between two superconductors. The latter, when taken separately, are described by wavefunctions of phases φ 1 , φ 2 and amplitudes | ψ 1 | , | ψ 2 | . We can explicitly write
ψ 1 = | ψ 1 | exp i φ 1 = ρ 1 exp i φ 1 , ψ 2 = | ψ 2 | exp i φ 2 = ρ 2 exp i φ 2 ,
where ρ 1 , ρ 2 are the probability amplitudes of Cooper pair densities.
Once the weak link is formed, coherence is established across the barrier, and the phase difference
Δ φ = φ 2 φ 1 = γ
determines interference between the (previously independent wavefunctions), so that the system can be described by means of a single wavefunction as a whole.
An SIS tunnel Josephson junction is a typical weak link consisting of two superconductors (that we take equal for simplicity) of thickness L and surface A, separated by a thin oxide layer of thickness L , see Figure 1. The time dependent Schrodinger equation can be used to characterize the system evolution as
i ψ t = E ψ .
As we already discussed, in the weak link, a coherent overlap takes place between the two wavefunctions, and an additional term must be added to take into account the interaction. In particular, the rate of change of ψ 1 is proportional to its coupling to ψ 2 , the same happening for ψ 2 on the other side. It is then possible to write the relations [98]
i ψ 1 t = E 1 ψ 1 + K ψ 2 ,
i ψ 2 t = E 2 ψ 2 + K ψ 1 .
For the sake of simplicity, we will consider a superconductor of the same kind, so that the probability amplitudes of Cooper pair densities are equal, ρ 1 = ρ 2 = ρ . The quantities E 1 and E 2 are the ground state energies of the unperturbed system (i.e. when K = 0 ) and we choose the zero of energy to be halfway between E 1 and E 2 , the evolution of the system then depending on the difference Δ E = E 2 E 1 .
Now, we put expressions (66) in the evolution relations Equations (69) and (70), assuming each wavefunction having a well-defined Cooper pair density and space-independent macroscopic phase. Separating the real and imaginary part, we find:
γ t = Δ E = 0 ,
ρ t = 2 K ρ sin ( γ ) ,
with γ = φ 2 φ 1 , the expression being valid in the absence of applied voltage of any kind (electric or gravitational-like) (Equation (72) is written in the standard Josephson formalism [99] with a little abuse of notation: the ρ -density involved in the time-derivative on the l.h.s. refers to the superconducting current density across the interface, while the density on the r.h.s. refers to the global density of Cooper pairs in the system; the latter is a conserved (constant) quantity, being the system in the superconductive state). The supercurrent across the contact then is written
J s = 2 e ρ t = 4 e K ρ sin ( γ ) = J 0 sin ( γ ) ,
flowing through the thin layer separating the superconductors and depending on the phase difference across the barrier.
If we apply a constant voltage Δ V across the junction, an oscillatory variation of phase difference takes place. A corresponding AC supercurrent then appears in the weak link, due to the existing finite potential difference in the junction. The phenomenon is a manifestation of the Josephson–Gor’kov principle [96,100,101], which simply states that the oscillation frequency of the coherent matter field is driven by the existing chemical potential (corresponding to pairs of electrons in the case of superconductivity). Since the supercurrent is a periodic function of Δ φ , AC supercurrents must be associated with any applied voltage difference (general principle of gauge invariance dictates that all physical properties must be periodic functions of the phase with period 2 π ). Equation (71) is then modified in order to take into account the applied voltage, resulting in a time-dependent relation of the form [95,96,97,101,102]:
γ t = 2 e Δ V ,
relating the phase difference variation on opposite sides to the existing potential difference across the junction. After integration, the above (74) gives
γ ( t ) = γ 0 + 2 e Δ V t ,
γ 0 being an integration constant. The supercurrent density turns out to be
J s = J 0 sin γ 0 + 2 e Δ V t .
and, below the critical temperature, the amplitude of the corresponding tunnelling supercurrent I 0 = J 0 A is temperature-dependent and is expressed by the Ambegaokar–Baratoff formula [103,104]:
I 0 = π Δ S ( T ) 2 e R N tanh Δ S ( T ) 2 k B T ,
where R N is the junction resistance in the normal state and Δ S ( T ) is the superconductive gap. The described AC signal, coming from the applied DC voltage, may be understood as the result of the energy conversion of electron pairs into photons [105].
In the following, we will examine the possibility of a Josephson-like effect induced by the weak-static Earth’s gravitational field, also analysing suitable experimental settings and parameter optimization.

3.2. Josephson Effect Induced by Gravity

We have discussed and motivated in the previous sections the introduction of generalized electric field and potential of the form
E = E e + m e E g V = V e + m e V g .
If we now restrict to a simple situation in which it is present only the Earth’s static gravitational field ( E e = 0 ), we have that
E = m e E g = m e g ,
while the corresponding potential difference reads
Δ V = m e Δ V g = 0 d z m e g = m e g ,
having chosen the z-axis along the direction of the gravitational field, see Figure 1. The resulting induced Josephson current [77] then has the form
I s ( t ) = I 0 sin 2 e Δ V t + φ = I 0 sin ( ω t + φ ) .
We also expect the induced effect to disappear when the junction is rotated in a position where the normal vector to the surface is perpendicular to the gravitational field direction.

Experimental Settings

Let us first consider a junction involving high- T c superconductors (HTSC). The latter have a coherence length ξ of the order of 10 9 m that fixes the thickness of the insulating layer to be ξ . Then, if we consider the pulsation
ω = 2 e Δ V = 2 m g ,
a junction with an insulating layer of thickness 1 nm would result in ω 1.7 × 10 4 s 1 , determining a corresponding period for the Josephson current T = 2 π / ω 3.7 × 10 4 s . This implies that the distinctive oscillatory behaviour can be observed only in very stable junctions, since a reasonable duration of the experiment turns out to be longer than one day, see Figure 2.
To reduce the time duration of the experiment, it is necessary to increase the voltage V g . Clearly, it is impossible to vary the intensity of the local gravitational field, so that the only strategy left is to have larger . This means that we need tunnelling junctions working in the presence of a thicker insulating layer: this is possible using low- T c superconductors (LTSC), and the latter can feature a larger coherent length, of the order of 10 3 nm . In this case, we can take an insulating layer of thickness 300 nm and obtain for the voltage V g 1.67 × 10 17 Volt . The pulsation turns out to be ω 0.05 s 1 and the corresponding period T 123 s , strongly reducing the experiment duration, see Figure 3.
From a practical point of view, it would be preferable to work with experimental setups stable enough to allow accurate oscillations measurements, but that, at the same time, give rise to a Josephson current of detectable intensity. If we increase the junction thickness using low- T c superconductors, the time duration for the experiment decreases and a stable setting is then possible, but the associated Josephson current becomes very weak and difficult to measure. Currently, the best choice to observe experimental evidence is to realize the most stable setup possible with HTCS, and then make long-time measurements of stronger Josephson currents.
We have seen how the proposed theoretical model provides the possibility to investigate the discussed interplay between gravitation and a superconductive condensate. For the simple case of the Josephson junction, the difficulties lie in the experimental setup that has to be stable in time to allow for careful observations of the oscillatory behaviour, and very sensitive to the induced voltage. In the following sections, we are going to analyse a more detailed microscopic description of the superfluid, exploiting a mean–field theory formulation for the system thermodynamics, including the effects of thermal fluctuations. In particular, we will analyse how the local gravitational field can be affected by the presence of a supercondensate exploiting the time-dependent Ginzburg–Landau equations in the regime of fluctuations.

4. Affecting the Field Just Outside the Sample: Ginzburg–Landau Formulation

We now want to better characterize the interaction between the superfluid and the local gravitational field. To this end, we need a microscopic quantum model describing the supercondensate behaviour. However, the formalism that characterizes the material superconductive state is in general very complicated, so that extracting quantitative predictions (or even just qualitative descriptions) for the interplay turns out to be an almost impossible task.
A simpler framework for analysing the interaction mechanism is given by a superconducting sample in the vicinity of its critical temperature T c . In particular, for T near T c , the system can be described by the Ginzburg–Landau equations, for which analytic solutions could be found.

4.1. Thermodynamic Fluctuations vs. Mean–Field Theory

The physics of low-temperature condensed matter systems is based on two fundamental notions: the low-energy long-living excitations (quasiparticles) and the mean field approximation. For instance, the BCS theory of superconductivity [106] is a paradigmatic example of the exploitation of both approaches mentioned.
Physical situations which cannot be consistently described in terms of the quasiparticle method or the mean field approximation are called fluctuations. The regime in which the fluctuations come into play is, in general, a very narrow temperature range around the critical temperature (On the contrary, for high temperature cuprate superconductors, organic superconductors, iron pnictides, low dimensional and amorphous superconducting systems, the situation changes radically due to the very small value of the coherence length, so that the temperature range of fluctuation is considerably larger). In particular, many effects on the superconducting phase occur while the system is still in the normal phase (just above the critical temperature) and originate from the appearance of the superconducting fluctuations themselves. In this regard, diamagnetic susceptibility, conductivity, heat capacity and other physical quantities may increase considerably near the transition temperature.
If we consider a range of temperature sufficiently far from the critical T c , the fluctuation regime ceases and the physics of the system is described in terms of a mean field formulation. The latter approach approximates the physics by averaging over the degrees of freedom of the system, that is, by approximating all the interactions acting on a single component with a single averaged effect. The technique allows for map a multi-body problem onto a one-body problem. In particular, the thermodynamic properties of the system are obtained by treating the order parameter as spatially constant, the spatial fluctuations being negligible. Many predictions can therefore be obtained by exploiting a much simpler mathematical formulation: this is a great advantage when dealing with new and unconventional systems, for which a complete description is not known.

GL Equations

As we have briefly discussed, the analysis of condensed matter systems in general involves the study of complicated, many electron states. A certain number of phenomenological approaches, based on classical field theory, were then developed to address the problem. A possibility is to consider a slowly-varying density of fields, carrying sufficient quantum information to write down an energy function for the system to be minimized: this corresponds to the celebrated Ginzburg–Landau (GL) formulation [107,108], based on a mean field approach. Its most notable use is in the theory of superconductors, where a complex scalar field ψ is used to characterize the density of the superconducting paired electrons. Even if the GL approach is, in general, superseded by the more fundamental BCS theory, it is a powerful tool in the vicinity of the critical temperature, where a more fundamental theory is lacking or the formulation is too complicated.

4.2. Ginzburg–Landau Formulation

Let us consider a superconductive sample near its critical temperature. At the microscopic level, thermodynamic fluctuations of the order parameter ψ ( x , t ) describing superconducting electrons occur, giving rise to localized regions of accelerated charge carriers [109,110,111,112,113]. From a physical point of view, ψ can be thought as the pseudowavefunction characterizing the motion of the center of mass of the Cooper pairs. The average size of these regions is much greater than the mean free path for a certain range of temperature above T c , while it decreases for larger temperature [114]. Moreover, we are going to consider sufficiently dirty materials, so that the effects of the fluctuations can be observed over a sizable range of temperature [111] (In order to have a sufficiently large temperature interval, the electronic mean free path characterizing the material in the normal state should be less than 10 Å ).
We now want to characterize in more detail the behaviour of the superconductive sample, also analysing its possible interaction with the surrounding gravitational field. If the sample is put at a temperature T slightly greater than T c but sufficiently far from the transition point (mean field regime), the system can be described in terms of linearized time-dependent Ginzburg–Landau equations (The order parameter being very small in the thermodynamic fluctuations regime, a linear order formulation can be exploited). The latter can be expressed in the gauge-invariant form as [115,116,117]:
Γ t 2 i e ϕ ψ = 1 2 m ϕ 2 i e A 2 ψ + α ψ ( T > T c ) .
where ψ ( x , t ) is the order parameter, ϕ ( x , t ) the electric potential and A ( x , t ) is the vector potential. We also introduce the quantities:
Γ = α ϵ ( T ) π 8 k B T c ϵ ( T ) = T T c T c ξ ( T ) = ξ 0 ϵ ( T ) α = 2 2 m ξ ( T ) 2
ξ 0 being the BCS intrinsic coherence length, roughly characterizing the smallest size of a wave packet formed by superconducting charge carriers. It plays a role analogous to the mean free path in the nonlocal electrodynamics of normal metals and is in general larger in metal superconductors (In spite of the fact that the two electrons in a Cooper pair can be far apart from each other, other electrons belonging to different Cooper pairs are usually closer). The temperature-dependent Ginzburg–Landau coherence length ξ ( T ) provides a measure of the distance over which the order parameter can vary without undue energy increase, for a given temperature T. Alternatively, it can be thought of as a characterization of the distance from the surface over which the order parameter is close to its bulk value.
We now consider the following ansatz for the solution:
ψ ( x , t ) = f ( x , t ) exp i g ( x , t ) ,
and obtain from (83) the relations
Γ f t = α f 1 2 m v s 2 f + 2 2 m Δ f ,
Γ f g t = 2 e Γ ϕ f 2 2 m f Δ g 2 v s · f
The superfluid speed v s has the form
v s = 1 m g + 2 e c A ,
and the associated supercurrent density j s reads
j s = 2 e m | ψ | 2 h g + 2 e c A = 2 e f 2 v s .

4.2.1. Thermodynamic Fluctuations

The presence of a thermal energy of the order of k B T implies that the system could fluctuate in different low-lying states with a non-zero probability. Let us then use f k to define the value of f for a fluctuation of the wave vector k . The above (86) can be recast in the form
Γ f k t = α f k 2 2 m k 2 f k 1 2 m v s 2 f k ,
v s t = 2 e m E
having used Equations (87), (86b) and
ϕ = E A t .
Equation (89b) can be easily integrated and the resulting expression for the superfluid speed can be used in (89a) giving
Γ f k t = α 2 2 m k 2 2 e 2 m E 2 t 2 f k .
We then find for f k
f k ( t ) = f k ( 0 ) exp α 2 2 m k 2 t 2 3 e 2 m E 2 t 3 Γ ,
with
f k 2 ( 0 ) = k B T 2 | α | + 2 2 m k 2 ,
and the associated current density j s k ( t ) can be written as
j s k ( t ) = 4 e 2 m E t f k 2 ( 0 ) exp 2 α 2 2 m k 2 t 2 3 e 2 m E 2 t 3 Γ .
Finally, the explicit expression for the physical supercurrent density j s [76] can be found integrating over k :
j s ( t ) = 1 8 π 3 0 + d k 4 π k 2 j s k ( k , t ) ,
where we have considered a three-dimensional sample of dirty material, whose dimensions are larger than the correlation length.

4.2.2. Generalized EM Fields

The above expression for the supercurrent density allows for extracting the explicit form of the generalized electromagnetic fields and potentials characterizing the physical evolution of the system. First of all, the vector potential A ( x , y , z , t ) is obtained from
A ( x , y , z , t ) = μ 0 4 π j s ( t ) d x d y d z ( x x ) 2 + ( y y ) 2 + ( z z ) 2 ,
where t is the retarded time
t = t ( x x ) 2 + ( y y ) 2 + ( z z ) 2 c .
The generalized electric field E ( x , y , z , t ) (39) is obtained from:
E ( x , y , z , t ) = A ( x , y , z , t ) t + m e g .
As we can appreciate, the generalized gravito–Maxwell E ( x , y , z , t ) features two contributions. In particular, the second term is the standard, constant weak Earth’s gravity contribution. On the other hand, the unconventional first term originates from the presence of the (non-constant) supercurrent density and can determine a local, additional contribution to the constant gravitational field g . The final result clearly depends on the superconducting sample shape and dimensions, as well as on the space point (outside the sample) where the gravitational fluctuation is measured.

4.3. Expected Effects

Let us now study in detail a suitable experimental setting to evaluate the proposed interplay. Here, we consider a superconductive disk at a temperature higher but very close to T c . The sample is kept in the normal state by a weak magnetic field that is then turned off at the time t = 0 , where the superconductive transition occurs. The axis of the disk is aligned with the direction of the gravitational field, the bases being parallel to the ground.
The chosen temperature regime ( T T c ) corresponds to the thermodynamic fluctuations regime we discussed in the previous section, so that we can exploit the corresponding results for the supercurrent and generalized EM fields expressions, see Equations (95), (96) and (98). We are interested in the gravitational correction along the axis of the disk, just above the upper base of cylindrical sample.
First, we consider the local alteration of the gravitational as a function of time. In Figure 4, we show the computed effect for an In sample. The latter is a low- T c metallic superconductor, then featuring a large intrinsic coherence length ξ 0 . The same analysis is then performed in Figure 5 for a Ba 0.4 K 0.6 Fe 2 As 2 sample, an high- T c superconductor with small ξ 0 . In both cases, the variation is measured along the disk axis, at a fixed distance d above the base surface. We can note that the local gravitational field is initially reduced with respect to the unperturbed value; then, it increases up to a maximum g + Δ for t = τ 0 and it finally relaxes to the standard unperturbed value g (We can also note that, for a very short time interval, the local field seems to change sign: this can be prevented by means of appropriate physical cutoffs, excluding the arbitrary growth of instabilities which would give rise to negative values [8]).
We then focus on the local alteration as a function of the distance from the sample for fixed time. In particular, we choose to maximize the effect putting ourselves at t = τ 0 . In Figure 6 and Figure 7, it is shown the variation, measured along the axis of the disk above the base surface, for the same In and Ba 0.4 K 0.6 Fe 2 As 2 samples. In both cases, the effect is stronger in the vicinity of the sample, as it seems reasonable.
From a preliminary qualitative analysis, it is possible to show that the maximum perturbation value Δ of the local field is proportional to inverse of the coherence length,
Δ ξ ( T ) 1 .
This suggests that a stronger affection can be obtained by using high– T c superconductors (the latter featuring smaller coherence length) and can be appreciated comparing the strength of the perturbation for low and high– T c superconducting samples in the presented Figures.
On the other hand, it is easily demonstrated that the maximal effect occurs after a time interval
τ 0 ( T T c ) 1 .
This means that the time range in which the perturbation takes place can be extended keeping the sample at a temperature close to the transition temperature. From this point of view, if we want to be very close to the effective critical T c , it could be easier to consider a low– T c sample, the temperature transition range being very narrow for the latter. However, this in turns results in a reduced alteration of the local field, since, close to T c , the Ginzburg–Landau coherence length ξ ( T ) diverges, see Equation (84).
In light of the above discussion, an optimized experimental settings should involve a large high– T c superconducting sample at a temperature very close to T c . The latter condition could help in extending the time range in which the effect takes place, while choosing a high– T c superconductor would determine an enhanced local alteration due to the short intrinsic coherence length. Finally, large dimensions for the sample give a larger integration range and a resulting stronger contribution.
The above considerations show how a careful arrangement of the experimental setup is very important, since the material parameters and the sample geometry, dimensions and temperature directly affect the magnitude of the interaction and the related time scales. In this regard, the very short time intervals in which the effect occurs complicate direct measurements.
In the following chapter, we will consider the possible affection of the local gravitational field in the sample interior exploiting again the effective framework of the gravito–Maxwell formulation combined with the Ginzburg–Landau formalism. The analysis will suggest that, in the superfluid region, a slight affection of the local field could take place, as we have discussed in Section 1 considering the formal quantum gravity point of view. A possibility to enhance the effect comes from the presence of suitable electric and magnetic fields, determining the formation of moving vortices and giving rise to a further interaction with the local gravitational field.

5. Affecting the Field Inside the Sample: Vortex Lattice

Now, we want to consider the possible alteration of the local static gravitational field in the region inside the superfluid. To this end, we will exploit the time-dependent Ginzburg–Landau equations for the supercondensate order parameter, looking for analytic solutions in the weak field condition. First, we will restrict to the simpler case of an isolated isotropic superconductor immersed in the Earth’s gravity in the absence of external EM fields. Then, we will analyse a more complicated setup, switching on suitable electric and magnetic fields: this will give rise to the formation of a vortex lattice inside the superfluid, possibly determining stronger effects for the proposed interplay.

5.1. Time-Dependent Ginzburg–Landau Formulation

Let us consider the case of a superconducting sample on the Earth surface. We already pointed out that the situation leads to the appearance of effective, generalized Maxwell fields. In particular, the local static weak gravitational field is treated as the gravitational component of the generalized gravitoelectric field, exploiting the formal analogy discussed in the previous Section 2.
The chosen physical system can be characterized in terms of time-dependent Ginzburg–Landau equations (TDGL). The latter are derived minimizing the total Gibbs free energy of the system [92,93,94], and can be written in a general explicit form as [120,121,122,123,124,125,126]:
2 2 m i + 2 e A 2 ψ a ψ + b | ψ | 2 ψ = 2 2 m D t + 2 i e ϕ ψ ,
× × A × B = μ 0 j n + j s
where j n and j s are expressed as
j n = σ A t + ϕ , j s = i e m ψ ψ ψ ψ 4 e 2 m | ψ | 2 A .
and correspond to the contributions of the normal current and supercurrent densities, respectively. In the above expressions, σ is the conductivity in the normal phase, D is the diffusion coefficient, B is the applied field, and the vector potential A is minimally coupled to ψ . The coefficients a and b in (101a) can be written as:
a = a ( T ) = a 0 ( T T c ) b = b ( T c )
where is T c the critical temperature of the superconductor, while a 0 and b are positive constant quantities. We can write consistent boundary and initial conditions for the system as
i ψ + 2 e A ψ · n = 0 × A · n = B · n A · n = 0 on Ω × ( 0 , t ) , ψ ( x , 0 ) = ψ 0 ( x ) A ( x , 0 ) = A 0 ( x ) on Ω ,
where Ω is the boundary of a smooth and simply connected domain in R N .

Dimensionless TDGL

The above Equations (101) can be recast in a useful dimensionless form. To this end, we define the following quantities:
Ψ 2 ( T ) = | a ( T ) | b ξ ( T ) = 2 m | a ( T ) | λ ( T ) = b m 4 μ 0 | a ( T ) | e 2 κ = λ ( T ) ξ ( T ) τ ( T ) = λ 2 ( T ) D η = μ 0 σ D B C ( T ) = μ 0 | a ( T ) | 2 b = 2 2 e λ ( T ) ξ ( T )
where λ ( T ) , ξ ( T ) and B C ( T ) are the penetration depth, coherence length and thermodynamic critical field, respectively. We also introduce the dimensionless quantities
t = t τ x = x λ y = y λ ψ = ψ Ψ
and the new dimensionless fields and currents
A = A κ 2 B C λ ϕ = ϕ κ 2 B C D E = E λ κ 2 B C D B = B κ 2 B C j = j μ 0 λ κ 2 B C .
We then insert the above Equations (106) and (107) in Equations (101) (we also drop the primes for the sake of notational simplicity) and get the dimensionless TDGL equations in a bounded, smooth and simply connected domain in R N [121,123]:
ψ t + i ϕ ψ + κ 2 | ψ | 2 1 ψ + ( i + A ) 2 ψ = 0 ,
× × A × B = j n + j s = η A t + ϕ i 2 ψ ψ ψ ψ | ψ | 2 A
while the boundary and initial conditions (104) in the dimensionless form read
i ψ + A ψ · n = 0 × A · n = B · n A · n = 0 on Ω × ( 0 , t ) ; ψ ( x , 0 ) = ψ 0 ( x ) A ( x , 0 ) = A 0 ( x ) on Ω .

5.2. Isolated Superconductor in the Weak Gravitational Field

Let us now now try to solve the above equations for a superconductor immersed in the Earth’s static gravity in the absence of external electromagnetic fields (39):
E e = 0 B e = 0 E = m e E g B = 0
having also set to zero the B g contribution that is negligible in the Solar system [28,127].

5.2.1. Solving TDGL Equations

A convenient gauge choice for subsequent calculations turns out to be ϕ = 0 , i.e., the vanishing of the scalar potential (Clearly, any alternative gauge shall not influence any physical results, the equations being gauge-invariant). From a physical point of view, this choice also reflects the absence of localized charges inside the superfluid, while contributions to the total gravitational field originating from the sample mass are clearly totally irrelevant. The dimensionless TDGL then explicitly read [128]:
ψ t = ( i + A ) 2 ψ κ 2 | ψ | 2 1 ψ ,
η A t = × × A + × B | ψ | 2 ( A θ )
where ψ ψ ( x , t ) is a complex function that we can write as
ψ = | ψ | exp ( i θ ) = Re ψ + i Im ψ = ψ 1 + i ψ 2 ,
so that (111a) splits into two distinct equations for the real and imaginary parts ψ 1 and ψ 2 .
Let us now restrict to a one-dimensional field configuration, so that one has
/ x A A x A .
In this simplified framework, (111) reads:
ψ 1 t = 2 ψ 1 x 2 + A ψ 2 x + ψ 2 A x ψ 1 A 2 κ 2 | ψ 1 | 2 + | ψ 2 | 2 1 ψ 1 , ψ 2 t = 2 ψ 2 x 2 A ψ 1 x ψ 1 A x ψ 2 A 2 κ 2 | ψ 1 | 2 + | ψ 2 | 2 1 ψ 2 , η A t = ψ 2 ψ 1 x ψ 1 ψ 2 x ψ 1 2 + ψ 2 2 A ,
since, in one dimension, 2 A = x · A and then
× × A = · A 2 A = 1 d 0 .
Then, let us consider an ideal, half-infinite superconductive region, see Figure 8. The u x direction is orthogonal to the superconducting separation surface, corresponding to the y z plane and parallel to the ground, so that, for x > 0 , we find an empty space, while the superfluid region is located at x 0 . The whole setting is immersed in the Earth’s uniform and static gravitational field that is captured by the gravitoelectric component
E g EXT = g u x ,
g being the standard gravity acceleration.
The dimensional form of the gravitoelectric field inside the superfluid region
E g = A g ( t ) t ,
while (116) suggests for the external (outside) gravitational vector potential the form
A g EXT ( t ) = g C + t u x ,
C being a constant.
In the 1D setup, the generalized external potential in the dimensionless form reads
A EXT = m e A g EXT κ 2 B C λ = g c 1 + t ,
where we have dropped the primes for notational simplicity. Using (105), we can also explicitly write
c 1 = C τ , g = m κ λ ( T ) g 2 e D B C ( T ) 1 .
Next, we express ψ 1 , ψ 2 and A as:
ψ 1 ( x , t ) = ψ 10 ( x ) + g γ 1 ( x , t ) ,
ψ 2 ( x , t ) = ψ 20 ( x ) + g γ 2 ( x , t ) ,
A ( x , t ) = g β ( x , t )
where ψ 10 and ψ 20 characterize the unperturbed system and satisfy
0 = 1 κ 2 2 ψ 10 x 2 + ψ 10 ψ 10 ψ 10 2 + ψ 20 2 ,
0 = 1 κ 2 2 ψ 20 x 2 + ψ 20 ψ 20 ψ 10 2 + ψ 20 2 .
the ψ 10 and ψ 20 behaviour therefore being described by equations of the same type.
We now choose to set
ψ 20 = 0 ψ 0 = ψ 10 + i ψ 20 = ψ 10 R ,
so that (122a) reads
0 = 1 κ 2 2 ψ 10 x 2 + ψ 10 ψ 10 3 ,
and is solved by [94]
ψ 10 = tanh κ x 2 .
We are therefore left with the following set of equations:
γ 1 t = 2 γ 1 x 2 + κ 2 1 3 ψ 10 2 γ 1 ,
γ 2 t = 2 γ 2 x 2 + κ 2 1 3 ψ 10 2 γ 2 β ψ 10 x ψ 10 β x ,
η β t = γ 2 ψ 10 x ψ 10 γ 2 x ψ 10 2 β ,
the last (126c) implying that β ( x , t ) does not depend on γ 1 ( x , t ) .
If we now decide to put ourselves away from borders, we can set ψ 10 1 in (126), resulting in
γ 1 t 2 γ 1 x 2 2 κ 2 γ 1 ,
γ 2 t 2 γ 2 x 2 2 κ 2 γ 2 β x ,
η β t γ 2 x β .
We then find for β the solution
β ( x , t ) = e t η b 1 ( x ) + 1 η 0 t d t e t η γ 2 ( x , t ) x .
where b 1 ( x ) = c 1 , as it is implied by Equation (121c) for t 0 .
Let us imagine that the sample transition to the superconducting state occurs at t = 0 . We also make the natural assumption that, before the transition, no alteration of the gravitational field takes place (material in the normal state), the gravitational field assuming the same value inside and outside the sample region for t < 0 . This results in the following boundary and initial conditions:
ψ ( 0 , t ) = 0 ψ ( x , 0 ) = ψ 10 ( x ) ψ 1 x ( x , 0 ) = 0 γ 1 ( 0 , t ) = 0 γ 1 ( x , 0 ) = 0 γ 1 x ( x , 0 ) = 0 γ 2 ( 0 , t ) = 0 γ 2 ( x , 0 ) = 0 γ 2 x ( x , 0 ) = 0
together with the condition for β
lim t 0 g β t ( x , t ) = g .
implying that the interplay occurs only in the presence of a superconducting phase.
In order to fix the dimensionless constant c 1 , we use Equations (117), (121c) and (126c) to write the relation between E g and β as
E g g = β t = 1 η γ 2 ψ 10 x ψ 10 γ 2 x + ψ 10 2 η β .
To satisfy the hypothesis that any affection of the gravitational field occurs only after the appearance of a superconducting phase ( t > 0 ), we assume
lim t 0 E g g = 1 ,
while, from the initial conditions in (129), we also have
lim t 0 γ 2 ( x , t ) = 0 lim t 0 γ 2 x ( x , t ) = 0 .
We then obtain
t 0 : 1 = ψ 10 2 η β ( x , 0 ) = ψ 10 2 η A EXT ( 0 ) g = ψ 10 2 η c 1 c 1 = η ψ 10 2 .
This c 1 constant is ineffective in empty space, while it is responsible for the desired, unconventional effects in the presence of the superconductor.
Finally, we can write the final form for β ( x , t ) away from borders ( ψ 10 1 , c 1 η ) :
β ( x , t ) = e t η η + 1 η 0 t d t e t η γ 2 ( x , t ) x ,
from which we obtain the ratio
E g g = β ( x , t ) t = 1 η e t η η + 1 η 0 t d t e t η γ 2 ( x , t ) x 1 η γ 2 ( x , t ) x .
The discussed formulation characterizes more explicitly the proposed interplay between gravity and supercondensates in the presented, simplified setup. First, we see that the external gravitational vector potential seems to play a role in the superconducting transition: in particular, the external constant c 1 tends to assume fixed values depending on the specific properties of the sample undergoing the superconducting transition. On the other hand, we expect the back-reaction on the local gravitational to take place only after the transition itself, when the vector potential begins to “perceive” the presence of a superfluid phase.

5.2.2. Expected Effects

The above (136) for the ratio E g / g can be used to estimate the value of gravitational field inside the superconductor just after the superconducting phase transition:
t 0 + : E g g 1 t η 1 η γ 2 ( x , 0 + ) x .
In the superconducting state, the alteration of the local field depends on physical characteristic of the involved sample. In particular, (137) shows that the relevant quantities are η , and the spatial derivative of γ 2 .
In order to enhance the interaction, we should maximize the variation γ 2 x , an effect than can be achieved by introducing suitable disorder in the material sample (This can be obtained, for instance, by means of chemical doping or proton irradiation). A maximized effect would also require small values for η . The latter is proportional to the product of the diffusion coefficient D times the conductivity just above T c , see (105). This would suggest to consider materials that are bad conductors in the normal state and have low Fermi energies (for example, cuprates).
Finally, we have to take into account the (usually very small) time scales in which the effect occurs, expressed by the τ coefficient
τ ( T ) = λ 2 ( T ) D with λ ( T ) λ 0 T c T T c .
The latter can be maximized with a reduced diffusion coefficient and large penetration length, as occurs in superconducting cuprates with internal disorder.
Performing measurements at a temperature close to T c would give rise to enhanced effects: for example, in the case of Bi 2 Sr 2 CaCu 2 O 8 ( T c 109 K , λ 0 500 nm , σ 1 3.6 × 10 6 Ω m , D 10 3 m 2 / s , ξ 0 1.4 nm [118]) for T 105 K , we find
T 105 K : τ 6.8 × 10 9 η 3.5 × 10 4 .
This would determine a reduction of the local gravitational field of the order of 2 × 10 5 , see Equation (137) neglecting the last term (Non-irradiated high– T c superconductors (like BSCCO) usually feature low disorder, resulting in reduced values for the spatial derivative of γ 2 ).
The above analysis shows how a perceptible affection of the local field inside the sample is possible even in a simplified setup (zero EM fields). Experimental difficulties may still arise from the short time intervals in which the effect manifests itself (see the previous Section 4.3). In addition, in this case, an appropriate choice of the material parameters is essential, in order to enhance the interaction and extend the time ranges to workable scales.
In the following section, we will analyse a more complicated setup involving external electric and magnetic fields, which in turn determine the presence of moving vortices. The new configuration will not only result in an additional affection of the local gravitational field, but also in the appearance of a new component of the generalized electric field inside the sample, parallel to the superconductor surface.

5.3. Switching on EM Fields: Vortex Lattice

We now consider a superconducting sample with finite thickness L and very large dimensions along u z and u y directions. The sample is immersed in an external magnetic field B 0 and has a square lattice of vortices, whose axes are directed along B 0 . We choose the latter as
B 0 = B 0 u z ,
together with a vector potential A of the form
A = B 0 x u y .
Here, we decide to work in the Coulomb gauge · A = 0 , where
2 A = μ 0 j n + j s .
We also allow for the presence of a constant external (standard) electric field E 0 ( e ) along the u x direction. Given the simultaneous presence of the Earth’s static gravity, the situation gives rise to a generalized static field E 0 of the form
E 0 = E 0 ( e ) + E 0 ( g ) = E 0 ( e ) E 0 ( g ) u x = E 0 ( e ) m e g u x = E 0 u x ,
and a related scalar potential
ϕ 0 = E 0 x .
As in the previous case, the transition takes place at t = 0 . In particular, for t < 0 , we also have T < T c and B > B C 2 , while, at t = 0 , we still have T < T c but B B C 2 . The external vector potential (outside the superfluid) is denoted by A 0 , and coincides with the inside value for t < 0 (sample in the normal state and very weakly diamagnetic material).

5.3.1. Linearized TDGL

In the new setup with non-zero external EM fields, it is possible to write an analytic approximate solution of the TDGL (101) for the order parameter as
ψ ( x , y , t ) = n = c n exp i q n y + E 0 B 0 t exp 1 2 ξ ( T ) x q n 2 e B 0 2 + i e E 0 ξ 2 ( T ) D x q n 2 e B 0 .
The expression is valid for an external magnetic field B 0 B C 2 and is then a solution of linearized TDGL equations [129,130] describing the behaviour of an ordered vortex lattice, moving under the influence of the external E 0 .
The above solution does not necessarily hold for different values of the magnetic field (for instance, B 0 B C 1 ), where the order parameter values are bigger and the linearized approximation does not hold. Moreover, close to B C 2 , the vortices are densely packed and the distance between them can be estimated to be of the order the coherence length ξ ( T ) . This can be then used to precisely characterize the vortex lattice, while this is not possible for generic values of B 0 [94]. (The presence of the external electric fields causes vortices motion and determines dissipative phenomena even in the superconducting state; it is possible to prevent it and anchor the vortices (vortex pinning) by introducing defects in the sample, thus reducing or eliminating energy dissipation [129]).
From an experimental point of view, in high- T c superconductors, the formation of a square lattice seems to be energetically favourable, and, in the following, we will restrict to this possibility (This is not the case for low- T c superconductors, where a triangular lattice formation usually occurs). We denote by q the distance between adjacent vortices that, for a square lattice, reads [131]
q 2 π ξ ( T ) ,
and the general c n coefficients could be replaced by the correspondent c expression for the square lattice:
c n c = 2 2 π ξ 2 ( T ) ,
the c coefficients being then independent of n.

5.3.2. Dimensionless Framework

Let us now consider the useful introduced dimensionless formulation. Working in the dimensionless version of the chosen Coulomb gauge · A = 0 it is possible to write a first-order expression for the dimensionless order parameter satisfying a linearized form for adimensional TDGL Equation (108) as [132]:
ψ ( x , y , t ) = n = | c n | exp i q n y + E 0 B 0 t exp κ 2 2 x n x 0 2 + i E 0 κ x n x 0 ,
with
| ψ | 2 = n = | c n | 2 exp κ 2 x n x 0 2 .
The equations for the vector potential components read
2 A x ( x , t ) x 2 = η A x ( x , t ) t E 0 + A x ( x , t ) E 0 κ n = | c n | 2 exp κ 2 ( x n x 0 ) 2 , 2 A y ( x , t ) x 2 = η A y ( x , t ) t + n = A y ( x , t ) 2 π κ n | c n | 2 exp κ 2 ( x n x 0 ) 2 , 2 A z ( x , t ) x 2 = η A z ( x , t ) t + n = | c n | 2 exp κ 2 ( x n x 0 ) 2 .
Let us now consider an expansion to linear order in E 0 . In order to obtain a more explicit solution for the order parameter (148), we have to estimate the summations
n = | c n | 2 exp κ 2 x n x 0 2 n = n | c n | 2 exp κ 2 x n x 0 2 .
Since we are interested in high- T c superconductors featuring a square vortex lattice, we replace the general coefficients c n with the correspondent c that, in the considered framework, reads [131]
c 2 = 2 2 π κ 2 .
being then a constant function of κ = λ / ξ .
For high- T c superconductors, the κ parameter is usually large, κ 2 10 4 : this in turn implies for the above (151):
n = | c n | 2 exp κ 2 x n x 0 2 = c 2 e κ 2 x 2 n = e κ 2 n 2 x 0 2 e 2 x x 0 n κ 2 c 2 e κ 2 x 2 , n = n | c n | 2 exp κ 2 x n x 0 2 0 ,
where the summation on the first line receives a non-negligible contribution only from the n = 0 term.
The Equation (150) for the vector potential can be now recast as
2 A x ( x , t ) x 2 = η A x ( x , t ) t E 0 + A x ( x , t ) E 0 κ c 2 e κ 2 x 2 , 2 A y ( x , t ) x 2 = η A y ( x , t ) t + A y ( x , t ) c 2 e κ 2 x 2 , 2 A z ( x , t ) x 2 = η A z ( x , t ) t + c 2 e κ 2 x 2 .
Since we are considering high- T c superconductors ( κ 2 10 4 ), it is also possible to approximate
e κ 2 x 2 π κ δ ( x ) ,
so that the above expressions read
A x ( x , t ) t 1 η 2 A x ( x , t ) x 2 A x ( x , t ) E 0 κ c 2 π η κ δ ( x ) + E 0
A y ( x , t ) t 1 η 2 A y ( x , t ) x 2 A y ( x , t ) c 2 π η κ δ ( x )
A z ( x , t ) t 1 η 2 A z ( x , t ) x 2 c 2 π η κ δ ( x )
The initial conditions for the vector potential components are:
A x x , 0 = 0 A y x , 0 = B 0 x A z x , 0 = 0
and the generalized electric field E inside the superfluid is given by
E = A t ϕ .

5.3.3. Averaged Solutions

We now consider the spatial averaged effects, determined by the presence of generalized field and poptentials, inside the supercondensate region. This can be obtained integrating the vector potential components (156) over the x-variable [133].
First, we integrate Equation (156c) over x in the interval x [ L / 2 , L / 2 ] , obtaining
A ¯ z ( t ) t = c 2 π η κ L ,
having introduced the averaged component
A ¯ z ( t ) = 1 L L / 2 L / 2 d x A z ( x , t ) ,
and taking advantage of symmetric conditions for the first derivatives with respect to x. Let us also keep in mind that we are dealing with the dimensionless quantities, having dropped the primes for the sake of notational simplicity (In particular, the x coordinate corresponds to the dimensionless x of (106), while one would explicitly have for the dimensionless thickness L = L / λ , L being the physical thickness and λ the penetration depth). The above (160) is solved by
A ¯ z ( t ) = c 2 π η κ L t + A ¯ z ( 0 ) = c 2 π η κ L t ,
where initial conditions (157) implies A ¯ z ( 0 ) = 0 . The averaged, generalized electric field E ¯ z component is then given by
E ¯ z = c 2 π η κ L = 2 2 π κ η L .
having used (158) and (144).
The averaged differential equation for the A ¯ y ( t ) component, defined in the same way as (160), is obtained from (156b) and reads
A ¯ y ( t ) t = A ¯ y ( t ) c 2 π η κ L ,
having used the approximation A y ( 0 , t ) A ¯ y ( t ) . The resulting averaged component reads
A ¯ y ( t ) = A ¯ y ( 0 ) exp c 2 π η κ L t = 0 ,
having again used initial condition (157). This also implies that the electric field E ¯ y ( t ) component is vanishing,
E ¯ y ( t ) = 0 .
The equation for the vertical component comes from the (156a) expression and reads
A ¯ x ( t ) t = A ¯ x ( t ) L E 0 κ c 2 π η κ + E 0 ,
Using again the approximation A x ( 0 , t ) A ¯ x ( t ) and the initial conditions (157), we find for the A ¯ x ( t ) solution
A ¯ x ( t ) = A ¯ x ( 0 ) exp c 2 π η κ L t + E 0 L κ + η κ L c 2 π 1 exp c 2 π η κ L t = = E 0 L κ + η κ L c 2 π 1 exp c 2 π η κ L t
Finally, the averaged E x ( t ) component along the vertical direction for the generalized electric field comes from Formulas (158) and (144) and reads
E ¯ x ( t ) = E 0 E 0 L κ + η κ L c 2 π c 2 π η κ L exp c 2 π η κ L t = = E 0 E 0 2 2 π η + 1 exp 2 2 π κ η L t .

5.3.4. Expected Effects

The analysis of the averaged effect inside the supercondensate region shows some interesting predictions.
The first effect is the emergence of a new component of the (generalized) electric field, parallel to the superconductor surface and directed along the external applied magnetic field. The value of this new contribution is found using the (dimensionless) result (162) together with Formula (107), and in dimensional units reads
E z = 4 π B C ( T ) D η L .
If we consider a Bi 2 Sr 2 Ca 3 Cu 3 O 10 sample ( T c 107 K , λ 0 2.4 × 10 7 m , ξ 0 1 nm , σ 1 3.6 × 10 6 Ω m , D 10 3 m 2 / s [134,135]) of thickness L = 15 cm at a temperature T = 102 K ; this would correspond to a resulting field
T 102 K : E z = 4 π B C ( T ) D η L = 4 π B C ( T ) μ 0 σ L 77 V m ,
with B c ( T ) 0.32 Tesla .
The second expected effect is affection of the local gravitational field along the x direction in the supercondensate region. The averaged effect is expressed by Equation (168), from which it is possible to appreciate the predicted, temporary alteration of the local field.
In Figure 9, we plot the field variation inside the superfluid region for two samples of different dimensions. Analogous with the results of Section 4 about the local alteration outside the material, we can see that, for very short time scales, the gravitational field has a non-negligible reduction. Clearly, sample dimensions and chemical composition play a key role in maximizing the effect.
First, we can appreciate that larger samples (i.e., larger values of L) would determine an increase of the time scales in which the effect manifests itself. In the same way, (168) suggests that large values of the η parameter, sample characteristics, determine an analogous increase of time ranges. The analysis then shows that L and η parameters determine similar effects: choosing a sample of reduced dimension (small L) of disordered material (small η , bad conductors in the normal state) would result in very short time scales, with a slight enhancement of the effect. Again, appropriate physical cutoffs should come into play, preventing non-physical growth of instabilities within the supercondensate, which would in turn lead to local field alterations of arbitrary intensity.
Since experimental issues would reside in the very short observation times, it is useful to take advantage of effects determined by the internal disorder. The effects of the latter can be easily understand, since material disorder causes an increase of the λ penetration depth. This, in turn, dictates an extension of the typical time scale τ of duration being τ λ 2 , see definitions (105) and (138).
Finally, if the system is put at temperatures very close to T c ; this is again an increase of the λ parameter and related larger time scales. In the latter case, however, the effects of thermal fluctuations should also be taken into account [136]. The described effects occur analogous with what we found in Section 4 for the affection of the local field just outside the sample.

6. Conclusions

A deeper intertwining of different scientific areas has always proved to be a powerful tool for improving our understanding of many fascinating physical aspects of our world, see e.g., [137,138,139,140,141,142,143,144,145,146,147,148,149,150,151,152,153,154,155,156,157,158,159]. The intriguing existence of an interplay between gravity and superconductivity has been investigated by many researchers in the last decades, due to the enormous conceptual implications and many possible applications. In particular, the interaction has been theoretically predicted by numerous authors, with very different approaches and techniques. The phenomenon was then successfully tested in relation to the effects of gravitational perturbation on supercurrents and supercondensates, having used the latter as “gravitational antennas” for the detection of gravitational waves.
In this review, we mainly focused on the possible back-reaction exerted by the superfluid on the surrounding gravitational field, trying to provide qualitative and quantitative predictions about the extent of the proposed effect. Inspired by theoretical and experimental studies on gravity-induced generalized fields in superconductors, we studied the possible alterations exploiting a gravito–Maxwell formalism, integrated with the Ginzburg–Landau theory of phase transitions for superconducting systems. The latter formalism is a phenomenological theory, the superconducting materials being characterizing by parameters which, in principle, can be optimized to enhance specific effects.
Clearly, there is still a lot of work to be done in order to better define the ranges and magnitude of the effect, as well as to determine optimal situations from an experimental point of view. In this regard, a crucial role would be played by suitable samples geometry, external electromagnetic fields of adequate frequency and appropriate characteristics of the material. In the future, 2D materials with variable number of layers should also be taken into consideration, in order to exploit their peculiar properties [160,161].

Author Contributions

Both authors contributed equally to conceptualization, methodology, formal analysis, original draft preparation and writing. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

G.A. Ummarino acknowledges partial support from the MEPhI. We also thank Fondazione CRT that partially supported this work for A. Gallerati.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. DeWitt, B.S. Superconductors and gravitational drag. Phys. Rev. Lett. 1966, 16, 1092–1093. [Google Scholar] [CrossRef]
  2. Papini, G. London moment of rotating superconductors and Lense-Thirring fields of general relativity. Il Nuovo Cimento B 1966, 45, 66–68. [Google Scholar] [CrossRef]
  3. Papini, G. Detection of inertial effects with superconducting interferometers. Phys. Lett. A 1967, 24, 32–33. [Google Scholar] [CrossRef]
  4. Hirakawa, H. Superconductors in gravitational field. Phys. Lett. A 1975, 53, 395–396. [Google Scholar] [CrossRef]
  5. Ciubotariu, C. Absorption of gravitational waves. Phys. Lett. A 1991, 158, 27–30. [Google Scholar] [CrossRef]
  6. Anandan, J. Relativistic gravitation and superconductors. Class. Quant. Grav. 1994, 11, 23. [Google Scholar] [CrossRef]
  7. Podkletnov, E.; Nieminen, R. A possibility of gravitational force shielding by bulk YBa2Cu3O7-X superconductor. Phys. C Supercond. 1992, 203, 441–444. [Google Scholar] [CrossRef]
  8. Modanese, G. Theoretical analysis of a reported weak gravitational shielding effect. Europhys. Lett. 1996, 35, 413–418. [Google Scholar] [CrossRef] [Green Version]
  9. Modanese, G. Role of a ‘local’ cosmological constant in Euclidean quantum gravity. Phys. Rev. D 1996, 54, 5002–5009. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Agop, M.; Buzea, C.; Griga, V.; Ciubotariu, C.; Stan, C.; Jatomir, D. Gravitational paramagnetism, diamagnetism and gravitational superconductivity. Aust. J. Phys. 1996, 49, 1063–1074. [Google Scholar] [CrossRef] [Green Version]
  11. Li, N.; Torr, D. Effects of a gravitomagnetic field on pure superconductors. Phys. Rev. D 1991, 43, 457. [Google Scholar] [CrossRef] [PubMed]
  12. Ahmedov, B. General relativistic thermoelectric effects in superconductors. Gen. Relativ. Gravit. 1999, 31, 357–369. [Google Scholar] [CrossRef]
  13. Agop, M.; Ioannou, P.; Diaconu, F. Some implications of gravitational superconductivity. Prog. Theor. Phys. 2000, 104, 733–742. [Google Scholar] [CrossRef] [Green Version]
  14. Modanese, G. Local contribution of a quantum condensate to the vacuum energy density. Mod. Phys. Lett. A 2003, 18, 683–690. [Google Scholar] [CrossRef]
  15. Wu, N. Gravitational shielding effects in gauge theory of gravity. Commun. Theor. Phys. 2004, 41, 567–572. [Google Scholar] [CrossRef] [Green Version]
  16. Hathaway, G.; Cleveland, B.; Bao, Y. Gravity modification experiment using a rotating superconducting disk and radio frequency fields. Phys. C Supercond. 2003, 385, 488. [Google Scholar] [CrossRef]
  17. Kiefer, C.; Weber, C. On the interaction of mesoscopic quantum systems with gravity. Ann. Phys. 2005, 14, 253–278. [Google Scholar] [CrossRef] [Green Version]
  18. Quach, J.Q. Gravitational Casimir effect. Phys. Rev. Lett. 2015, 114, 081104, Erratum: Phys. Rev. Lett. 2017, 118, 139901. [Google Scholar] [CrossRef]
  19. Ummarino, G.A.; Gallerati, A. Superconductor in a weak static gravitational field. Eur. Phys. J. C 2017, 77, 549. [Google Scholar] [CrossRef]
  20. Atanasov, V. The geometric field (gravity) as an electro-chemical potential in a Ginzburg-Landau theory of superconductivity. Phys. B Condens. Matter 2017, 517, 53–58. [Google Scholar] [CrossRef] [Green Version]
  21. Atanasov, V. Gravitation at the Josephson junction. Adv. Cond. Matt. Phys. 2018, 2018, 1618252. [Google Scholar] [CrossRef] [Green Version]
  22. Papini, G. Superconducting and normal metals as detectors of gravitational waves. Lett. Nuovo Cim. 1970, 4S1, 1027–1032. [Google Scholar] [CrossRef]
  23. Adler, R.J. Long conductors as antennae for gravitational radiation. Nature 1976, 259, 296–297. [Google Scholar] [CrossRef]
  24. Anandan, J. Relativistic thermoelectromagnetic gravitational effects in normal conductors and superconductors. Phys. Lett. A 1984, 105, 280–284. [Google Scholar] [CrossRef]
  25. Anandan, J. Detection of gravitational radiation using superconducting circuits. Phys. Lett. A 1985, 110, 446–450. [Google Scholar] [CrossRef]
  26. Carelli, P.; Castellano, M.; Cosmelli, C.; Foglietti, V.; Modena, I. Coupling of a high-sensitivity superconducting amplifier to a gravitational-wave antenna. Phys. Rev. A 1985, 32, 3258. [Google Scholar] [CrossRef] [PubMed]
  27. Chan, H.; Paik, H. Superconducting gravity gradiometer for sensitive gravity measurements. I. Theory. Phys. Rev. D 1987, 35, 3551. [Google Scholar] [CrossRef] [PubMed]
  28. Mashhoon, B.; Paik, H.J.; Will, C.M. Detection of the gravitomagnetic field using an orbiting superconducting gravity gradiometer. Theoretical principles. Phys. Rev. D 1989, 39, 2825. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Preparata, G. ‘Superradiance’ Effects in a Gravitational Antenna. Mod. Phys. Lett. A 1990, 5, 1. [Google Scholar] [CrossRef]
  30. Peng, H. The effects of gravitational waves on a superconducting antenna and its sensitivity. Gen. Rel. Grav. 1990, 22, 33–43. [Google Scholar] [CrossRef]
  31. Peng, H.; Torr, D. The Electric field induced by a gravitational wave in a superconductor: A Principle for a new gravitational wave antenna. Gen. Rel. Grav. 1990, 22, 53–59. [Google Scholar] [CrossRef]
  32. Peng, H.; Chin, Y.; Lind, G. Interaction between gravity and moving superconductors. Gen. Rel. Grav. 1991, 23, 1231–1250. [Google Scholar] [CrossRef]
  33. Peng, H.; Torr, D.; Hu, E.; Peng, B. Electrodynamics of moving superconductors and superconductors under the influence of external forces. Phys. Rev. B 1991, 43, 2700. [Google Scholar] [CrossRef] [PubMed]
  34. Li, F.; Baker, R., Jr. Detection of high-frequency gravitational waves by superconductors. Int. J. Mod. Phys. B 2007, 21, 3274–3278. [Google Scholar] [CrossRef]
  35. Minter, S.J.; Wegter-McNelly, K.; Chiao, R.Y. Do Mirrors for Gravitational Waves Exist? Phys. E 2010, 42, 234. [Google Scholar] [CrossRef] [Green Version]
  36. Inan, N.; Thompson, J.; Chiao, R. Interaction of gravitational waves with superconductors. Fortschr. Phys. 2017, 65, 1600066. [Google Scholar] [CrossRef]
  37. Inan, N. A new approach to detecting gravitational waves via the coupling of gravity to the zero-point energy of the phonon modes of a superconductor. Int. J. Mod. Phys. D 2017, 26, 1743031. [Google Scholar] [CrossRef]
  38. Hammad, F.; Landry, A. A simple superconductor quantum interference device for testing gravity. Mod. Phys. Lett. A 2020, 35, 2050171. [Google Scholar] [CrossRef]
  39. Overhauser, A.; Colella, R. Experimental test of gravitationally induced quantum interference. Phys. Rev. Lett. 1974, 33, 1237. [Google Scholar] [CrossRef]
  40. Colella, R.; Overhauser, A.; Werner, S. Observation of gravitationally induced quantum interference. Phys. Rev. Lett. 1975, 34, 1472–1474. [Google Scholar] [CrossRef]
  41. Anandan, J. Gravitational and Rotational Effects in Quantum Interference. Phys. Rev. D 1977, 15, 1448–1457. [Google Scholar] [CrossRef]
  42. Anandan, J. Interference, Gravity and Gauge Fields. Nuovo Cim. A 1979, 53, 221. [Google Scholar] [CrossRef]
  43. Anandan, J. Gravitationally Coupled Electromagnetic Systems and Quantum Interference. Class. Quant. Grav. 1984, 1, L51. [Google Scholar] [CrossRef]
  44. Cai, Y.; Papini, G. Particle Interferometry in Weak Gravitational Fields. Class. Quant. Grav. 1989, 6, 407. [Google Scholar] [CrossRef]
  45. Ahluwalia, D.; Burgard, C. Gravitationally induced quantum mechanical phases and neutrino oscillations in astrophysical environments. Gen. Rel. Grav. 1996, 28, 1161–1170. [Google Scholar] [CrossRef]
  46. Bhattacharya, T.; Habib, S.; Mottola, E. Gravitationally induced neutrino oscillation phases in static space-times. Phys. Rev. D 1999, 59, 067301. [Google Scholar] [CrossRef] [Green Version]
  47. Müntinga, H.; Ahlers, H.; Krutzik, M.; Wenzlawski, A.; Arnold, S.; Becker, D.; Bongs, K.; Dittus, H.; Duncker, H.; Gaaloul, N.; et al. Interferometry with Bose-Einstein Condensates in Microgravity. Phys. Rev. Lett. 2013, 110, 093602. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Asenbaum, P.; Overstreet, C.; Kovachy, T.; Brown, D.D.; Hogan, J.M.; Kasevich, M.A. Phase Shift in an Atom Interferometer due to Spacetime Curvature across its Wave Function. Phys. Rev. Lett. 2017, 118, 183602. [Google Scholar] [CrossRef] [Green Version]
  49. Kiefer, C.; Singh, T.P. Quantum gravitational corrections to the functional Schrodinger equation. Phys. Rev. D 1991, 44, 1067–1076. [Google Scholar] [CrossRef] [PubMed]
  50. Sakurai, J.J.; Napolitano, J. Modern Quantum Mechanics; Cambridge University Press: Cambridge, UK, 2017. [Google Scholar]
  51. Schiff, L.; Barnhill, M. Gravitation-induced electric field near a metal. Phys. Rev. 1966, 151, 1067. [Google Scholar] [CrossRef]
  52. Witteborn, F.; Fairbank, W. Experimental comparison of the gravitational force on freely falling electrons and metallic electrons. Phys. Rev. Lett. 1967, 19, 1049. [Google Scholar] [CrossRef]
  53. Witteborn, F.; Fairbank, W. Experiments to determine the force of gravity on single electrons and positrons. Nature 1968, 220, 436–440. [Google Scholar] [CrossRef]
  54. Beams, J. Potentials on rotor surfaces. Phys. Rev. Lett. 1968, 21, 1093. [Google Scholar] [CrossRef]
  55. Herring, C. Gravitationally induced electric field near a conductor, and its relation to the surface-stress concept. Phys. Rev. 1968, 171, 1361. [Google Scholar] [CrossRef]
  56. Peshkin, M. Gravity-induced electric field near a conductor. Ann. Phys. 1968, 46, 1–11. [Google Scholar] [CrossRef]
  57. Peshkin, M. Gravity-induced electric field near a conductor. Phys. Lett. A 1969, 29, 181–182. [Google Scholar] [CrossRef]
  58. Craig, P.P. Direct observation of stress-induced shifts in contact potentials. Phys. Rev. Lett. 1969, 22, 700. [Google Scholar] [CrossRef]
  59. Rieger, T. Gravitationally Induced Electric Field in Metals. Phys. Rev. B 1970, 2, 825. [Google Scholar] [CrossRef]
  60. Leung, M. Electric fields induced by gravitational fields in metals. Il Nuovo Cimento B 1972, 7, 220–224. [Google Scholar] [CrossRef]
  61. Lockhart, J.; Witteborn, F.; Fairbank, W. Evidence for a temperature-dependent surface shielding effect in Cu. Phys. Rev. Lett. 1977, 38, 1220. [Google Scholar] [CrossRef]
  62. Anandan, J. New relativistic gravitational effects using charged-particle interferometry. Gen. Rel. Grav. 1984, 16, 33–41. [Google Scholar] [CrossRef]
  63. Peng, H. On calculation of magnetic-type gravitation and experiments. Gen. Relativ. Gravit. 1983, 15, 725–735. [Google Scholar] [CrossRef]
  64. Jain, A.; Lukens, J.; Tsai, J. Test for relativistic gravitational effects on charged particles. Phys. Rev. Lett. 1987, 58, 1165–1168. [Google Scholar] [CrossRef] [PubMed]
  65. Li, N.; Torr, D. Gravitational effects on the magnetic attenuation of superconductors. Phys. Rev. B 1992, 46, 5489. [Google Scholar] [CrossRef] [PubMed]
  66. Harris, E. Analogy between general relativity and electromagnetism for slowly moving particles in weak gravitational fields. Am. J. Phys. 1991, 59, 421–425. [Google Scholar] [CrossRef]
  67. Torr, D.; Li, N. Gravitoelectric-electric coupling via superconductivity. Found. Phys. Lett. 1993, 6, 371–383. [Google Scholar] [CrossRef]
  68. Agop, M.; Buzea, C.; Nica, P. Local gravitoelectromagnetic effects on a superconductor. Phys. C Supercond. 2000, 339, 120–128. [Google Scholar] [CrossRef]
  69. Tajmar, M.; De Matos, C. Gravitomagnetic field of a rotating superconductor and of a rotating superfluid. Phys. C 2003, 385, 551–554. [Google Scholar] [CrossRef] [Green Version]
  70. Tajmar, M.; de Matos, C. Extended analysis of gravitomagnetic fields in rotating superconductors and superfluids. Phys. C 2005, 420, 56. [Google Scholar] [CrossRef] [Green Version]
  71. Ahmedov, B.; Kagramanova, V. Electromagnetic effects in superconductors in stationary gravitational field. Int. J. Mod. Phys. D 2005, 14, 837–847. [Google Scholar] [CrossRef] [Green Version]
  72. De Matos, C.J. Gravitational force between two electrons in superconductors. Phys. C Supercond. 2008, 468, 229–232. [Google Scholar] [CrossRef] [Green Version]
  73. Tajmar, M. Electrodynamics in superconductors explained by Proca equations. Phys. Lett. A 2008, 372, 3289–3291. [Google Scholar] [CrossRef] [Green Version]
  74. Misner, C.W.; Thorne, K.; Wheeler, J. Gravitation; W. H. Freeman: San Francisco, CA, USA, 1973. [Google Scholar]
  75. Wald, R.M. General Relativity; Chicago University Press: Chicago, IL, USA, 1984. [Google Scholar] [CrossRef]
  76. Ummarino, G.A.; Gallerati, A. Exploiting weak field gravity-Maxwell symmetry in superconductive fluctuations regime. Symmetry 2019, 11, 1341. [Google Scholar] [CrossRef] [Green Version]
  77. Ummarino, G.A.; Gallerati, A. Josephson AC effect induced by weak gravitational field. Class. Quant. Grav. 2020, 37, 217001. [Google Scholar] [CrossRef]
  78. Braginsky, V.B.; Caves, C.M.; Thorne, K.S. Laboratory Experiments to Test Relativistic Gravity. Phys. Rev. D 1977, 15, 2047. [Google Scholar] [CrossRef] [Green Version]
  79. Ross, D. The London equations for superconductors in a gravitational field. J. Phys. A 1983, 16, 1331. [Google Scholar] [CrossRef]
  80. Thorne, K. Gravitomagnetism, jets in quasars, and the stanford gyroscope experiment. In Near Zero: New Frontiers of Physics; W.H. Freeman & Co.: New York, NY, USA, 1988; pp. 573–586. [Google Scholar]
  81. Peng, H. A new approach to studying local gravitomagnetic effects on a superconductor. Gen. Relativ. Gravit. 1990, 22, 609–617. [Google Scholar] [CrossRef]
  82. Ruggiero, M.L.; Tartaglia, A. Gravitomagnetic effects. Nuovo Cim. B 2002, 117, 743–768. [Google Scholar]
  83. Tartaglia, A.; Ruggiero, M.L. Gravitoelectromagnetism versus electromagnetism. Eur. J. Phys. 2004, 25, 203–210. [Google Scholar] [CrossRef]
  84. Vieira, R.; Brentan, H. Covariant theory of gravitation in the framework of special relativity. Eur. Phys. J. Plus 2018, 133, 165. [Google Scholar] [CrossRef] [Green Version]
  85. Behera, H. Comments on gravitoelectromagnetism of Ummarino and Gallerati in “Superconductor in a weak static gravitational field” vs other versions. Eur. Phys. J. C 2017, 77, 822. [Google Scholar] [CrossRef] [Green Version]
  86. Giardino, S. A novel covariant approach to gravito-electromagnetism. Braz. J. Phys. 2020, 50, 372–378. [Google Scholar] [CrossRef] [Green Version]
  87. Sbitnev, V.I. Quaternion algebra on 4D superfluid quantum space-time. Gravitomagnetism. Found. Phys. 2019, 49, 107–143. [Google Scholar] [CrossRef] [Green Version]
  88. Gallerati, A. Interaction between superconductors and weak gravitational field. J. Phys. Conf. Ser. 2020, 1690, 012141. [Google Scholar] [CrossRef]
  89. Williams, L.L.; Inan, N. Maxwellian mirages in general relativity. New J. Phys. 2021, 23, 053019. [Google Scholar] [CrossRef]
  90. Gallerati, A. Local affection of weak gravitational field from supercondensates. Phys. Scripta 2021, 96, 064001. [Google Scholar] [CrossRef]
  91. Toth, G.Z. Energy-momentum tensor and duality symmetry of linearized gravity in a Maxwellian formalism. arXiv 2021, arXiv:2108.02124. [Google Scholar]
  92. De Gennes, P.G. Superconductivity of Metals and Alloys; Taylor & Francis Ltd.: London, UK, 2018. [Google Scholar] [CrossRef]
  93. Tinkham, M. Introduction to Superconductivity; Dover Publications Inc.: New York, NY, USA, 2004. [Google Scholar]
  94. Ketterson, J.; Song, S. Superconductivity; Cambridge University Press: Cambridge, UK, 1999. [Google Scholar] [CrossRef]
  95. Josephson, B. Possible new effects in superconductive tunnelling. Phys. Lett. 1962, 1, 251–253. [Google Scholar] [CrossRef]
  96. Anderson, P. The Josephson Effect and Quantum Coherence Measurements in Superconductors and Superfluids. In Progress in Low Temperature Physics; Elsevier: Amsterdam, NL, USA, 1967; Volume 5, pp. 1–43. [Google Scholar]
  97. Barone, A.; Paternò, G. Physics and Applications of the Josephson Effect; John Wiley & Sons: New York, NY, USA, 1982. [Google Scholar] [CrossRef]
  98. Feynman, R.; Leighton, R.; Sands, M. The Josephson junction. In The Feynman Lectures on Physics; Addison-Wesley Publ. Comp.: New York, NY, USA, 1965; Chapter 21; Volume III. [Google Scholar]
  99. Fossheim, K.; Sudbø, A. Superconductivity: Physics and Applications; John Wiley & Sons Ltd.: Hoboken, NJ, USA, 2004. [Google Scholar] [CrossRef]
  100. Gor’kov, L. Microscopic derivation of the Ginzburg-Landau equations in the theory of superconductivity. Sov. Phys. JETP 1959, 9, 1364–1367. [Google Scholar]
  101. Josephson, B. Supercurrents through barriers. Adv. Phys. 1965, 14, 419–451. [Google Scholar] [CrossRef]
  102. Josephson, B. Coupled superconductors. Rev. Mod. Phys. 1964, 36, 216. [Google Scholar] [CrossRef]
  103. Ambegaokar, V.; Baratoff, A. tunnelling Between Superconductors. Phys. Rev. Lett. 1963, 10, 486. [Google Scholar] [CrossRef]
  104. Ambegaokar, V.; Baratoff, A. tunnelling Between Superconductors (Errata). Phys. Rev. Lett. 1963, 11, 104. [Google Scholar] [CrossRef]
  105. Saxena, A.K. The Proximity and Josephson Effects. In High-Temperature Superconductors; Springer: Berlin/Heidelberg, Germany, 2009; pp. 147–198. [Google Scholar] [CrossRef]
  106. Bardeen, J.; Cooper, L.; Schrieffer, J. Theory of superconductivity. Phys. Rev. 1957, 108, 1175–1204. [Google Scholar] [CrossRef] [Green Version]
  107. Ginzburg, V.; Landau, L. On the Theory of superconductivity. Zh. Eksp. Teor. Fiz. 1950, 20, 1064–1082. [Google Scholar]
  108. Ginzburg, V.; Landau, L. On the Theory of Superconductivity. In On Superconductivity and Superfluidity; Springer: Berlin/Heidelberg, Germany, 2009; pp. 113–137. [Google Scholar] [CrossRef]
  109. Ginzburg, V. Some remarks on phase transitions of the second kind and the microscopic theory of ferroelectric materials. Soviet Phys. Solid State 1961, 2, 1824–1834. [Google Scholar]
  110. Thouless, D.J. Perturbation theory in statistical mechanics and the theory of superconductivity. Ann. Phys. 1960, 10, 553–588. [Google Scholar] [CrossRef] [Green Version]
  111. Shier, J.; Ginsberg, D. Superconducting transitions of amorphous Bismuth alloys. Phys. Rev. 1966, 147, 384. [Google Scholar] [CrossRef]
  112. Glover, R. Ideal resistive transition of a superconductor. Phys. Lett. A 1967, 25, 542–544. [Google Scholar] [CrossRef]
  113. Strongin, M.; Kammerer, O.; Crow, J.; Thompson, R.; Fine, H. ‘Curie-Weiss’ behavior and fluctuation phenomena in the resistive transitions of dirty superconductors. Phys. Rev. Lett. 1968, 20, 922. [Google Scholar] [CrossRef]
  114. Ferrell, R.; Schmidt, H. Predicted critical behavior near the superconducting phase transition. Phys. Lett. A 1967, 25, 544–545. [Google Scholar] [CrossRef]
  115. Schmid, A. A time dependent Ginzburg-Landau equation and its application to the problem of resistivity in the mixed state. Phys. Der Kondens. Mater. 1966, 5, 302–317. [Google Scholar] [CrossRef]
  116. Hurault, J. Nonlinear Effects on the Conductivity of a Superconductor above Its Transition Temperature. Phys. Rev. 1969, 179, 494. [Google Scholar] [CrossRef]
  117. Schmid, A. Diamagnetic susceptibility at the transition to the superconducting state. Phys. Rev. 1969, 180, 527. [Google Scholar] [CrossRef]
  118. Poole, C.K.; Farach, H.A.; Creswick, R.J. Handbook of Superconductivity; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  119. Welp, U.; Xie, R.; Koshelev, A.; Kwok, W.; Luo, H.; Wang, Z.; Mu, G.; Wen, H.H. Anisotropic phase diagram and strong coupling effects in Ba1-xKxFe2As2 from specific-heat measurements. Phys. Rev. B 2009, 79, 094505. [Google Scholar] [CrossRef]
  120. Ullah, S.; Dorsey, A. Effect of fluctuations on the transport properties of type-II superconductors in a magnetic field. Phys. Rev. B 1991, 44, 262. [Google Scholar] [CrossRef] [PubMed]
  121. Tang, Q.; Wang, S. Time dependent Ginzburg-Landau equations of superconductivity. Phys. D Nonlinear Phenom. 1995, 88, 139–166. [Google Scholar] [CrossRef]
  122. Du, Q.; Gray, P. High-kappa limits of the time-dependent Ginzburg-Landau model. SIAM J. Appl. Math. 1996, 56, 1060–1093. [Google Scholar] [CrossRef]
  123. Lin, F.H.; Du, Q. Ginzburg-Landau vortices: Dynamics, pinning, and hysteresis. SIAM J. Math. Anal. 1997, 28, 1265–1293. [Google Scholar] [CrossRef] [Green Version]
  124. Fleckinger-Pellé, J.; Kaper, H.G.; Takáč, P. Dynamics of the Ginzburg-Landau equations of superconductivity. Nonlinear Anal. Theory Methods Appl. 1998, 32, 647–665. [Google Scholar] [CrossRef] [Green Version]
  125. Kopnin, N.; Thuneberg, E. Time-dependent Ginzburg-Landau analysis of inhomogeneous normal-superfluid transitions. Phys. Rev. Lett. 1999, 83, 116. [Google Scholar] [CrossRef]
  126. Ghinovker, M.; Shapiro, I.; Shapiro, B.Y. Explosive nucleation of superconductivity in a magnetic field. Phys. Rev. B 1999, 59, 9514. [Google Scholar] [CrossRef]
  127. Ljubičić, A.; Logan, B. A proposed test of the general validity of Mach’s principle. Phys. Lett. A 1992, 172, 3–5. [Google Scholar] [CrossRef]
  128. Ummarino, G.A.; Gallerati, A. Possible alterations of local gravitational field inside a superconductor. Entropy 2021, 23, 193. [Google Scholar] [CrossRef] [PubMed]
  129. Kopnin, N.; Ivlev, B.; Kalatsky, V. The flux-flow Hall effect in type II superconductors. An explanation of the sign reversal. J. Low Temp. Phys. 1993, 90, 1–13. [Google Scholar] [CrossRef]
  130. Kopnin, N. Theory of Nonequilibrium Superconductivity; Oxford University Press: Oxford, UK, 2001. [Google Scholar] [CrossRef] [Green Version]
  131. Hoffmann, K.; Tang, Q. Ginzburg-Landau Phase Transition Theory and Superconductivity; Springer: Basel, Switzerland, 2012. [Google Scholar] [CrossRef]
  132. Ummarino, G.A.; Gallerati, A. Superconductor in static gravitational, electric and magnetic fields with vortex lattice. Results Phys. 2021, 30, 104838. [Google Scholar] [CrossRef]
  133. Sanders, J.A.; Verhulst, F.; Murdock, J. Averaging Methods in Nonlinear Dynamical Systems; Springer: New York, NY, USA, 2007. [Google Scholar] [CrossRef]
  134. Weigand, M.; Eisterer, M.; Giannini, E.; Weber, H. Mixed state properties of Bi2Sr2Ca2Cu3O10+δ single crystals before and after neutron irradiation. Phys. Rev. B 2010, 81, 014516. [Google Scholar] [CrossRef]
  135. Piriou, A.; Fasano, Y.; Giannini, E.; Fischer, Ø. Effect of oxygen-doping on Bi2Sr2Ca2Cu3O10+δ vortex matter: Crossover from electromagnetic to Josephson interlayer coupling. Phys. Rev. B 2008, 77, 184508. [Google Scholar] [CrossRef] [Green Version]
  136. Larkin, A.; Varlamov, A. Theory of fluctuations in superconductors; Oxford University Press: Oxford, UK, 2005. [Google Scholar] [CrossRef]
  137. Zurek, W. Cosmological experiments in condensed matter systems. Phys. Rept. 1996, 276, 177–221. [Google Scholar] [CrossRef] [Green Version]
  138. Volovik, G. Superfluid 3He-B and gravity. Phys. B Condens. Matter 1990, 162, 222–230. [Google Scholar] [CrossRef]
  139. Volovik, G. Superfluid analogies of cosmological phenomena. Phys. Rept. 2001, 351, 195–348. [Google Scholar] [CrossRef] [Green Version]
  140. Baeuerle, C.; Bunkov, Y.M.; Fisher, S.N.; Godfrin, H.; Pickett, G.R. Laboratory simulation of cosmic string formation in the early Universe using superfluid He-3. Nature 1996, 382, 332–334. [Google Scholar] [CrossRef]
  141. Ruutu, V.; Eltsov, V.; Gill, A.; Kibble, T.; Krusius, M.; Makhlin, Y.; Placais, B.; Volovik, G.; Xu, W. Big bang simulation in superfluid He-3-b: Vortex nucleation in neutron irradiated superflow. Nature 1996, 382, 334. [Google Scholar] [CrossRef] [Green Version]
  142. Garay, L.J.; Anglin, J.R.; Cirac, J.I.; Zoller, P. Black holes in Bose-Einstein condensates. Phys. Rev. Lett. 2000, 85, 4643–4647. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  143. Jacobson, T.A.; Volovik, G.E. Event horizons and ergoregions in He-3. Phys. Rev. D 1998, 58, 064021. [Google Scholar] [CrossRef] [Green Version]
  144. Barcelo, C.; Liberati, S.; Visser, M. Analog gravity from Bose-Einstein condensates. Class. Quant. Grav. 2001, 18, 1137. [Google Scholar] [CrossRef] [Green Version]
  145. Novello, M.; Visser, M.; Volovik, G.E. Artificial Black Holes; World Scientific: Singapore, 2002. [Google Scholar] [CrossRef]
  146. Barcelo, C.; Liberati, S.; Visser, M. Analogue gravity. Living Rev. Rel. 2005, 8, 12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  147. Carusotto, I.; Fagnocchi, S.; Recati, A.; Balbinot, R.; Fabbri, A. Numerical observation of Hawking radiation from acoustic black holes in atomic Bose-Einstein condensates. New J. Phys. 2008, 10, 103001. [Google Scholar] [CrossRef]
  148. Mannarelli, M.; Manuel, C. Transport theory for cold relativistic superfluids from an analogue model of gravity. Phys. Rev. D 2008, 77, 103014. [Google Scholar] [CrossRef] [Green Version]
  149. Boada, O.; Celi, A.; Latorre, J.I.; Lewenstein, M. Dirac Equation For Cold Atoms In Artificial Curved Spacetimes. New J. Phys. 2011, 13, 035002. [Google Scholar] [CrossRef] [Green Version]
  150. Gallerati, A. Graphene properties from curved space Dirac equation. Eur. Phys. J. Plus 2019, 134, 202. [Google Scholar] [CrossRef]
  151. Capozziello, S.; Pincak, R.; Saridakis, E.N. Constructing superconductors by graphene Chern-Simons wormholes. Annals Phys. 2018, 390, 303–333. [Google Scholar] [CrossRef]
  152. Andrianopoli, L.; Cerchiai, B.L.; D’Auria, R.; Gallerati, A.; Noris, R.; Trigiante, M.; Zanelli, J. N-extended D=4 supergravity, unconventional SUSY and graphene. JHEP 2020, 1, 84. [Google Scholar] [CrossRef] [Green Version]
  153. Gallerati, A. Supersymmetric theories and graphene. PoS 2021, 390, 662. [Google Scholar] [CrossRef]
  154. Zaanen, J.; Liu, Y.; Sun, Y.W.; Schalm, K. Holographic Duality in Condensed Matter Physics; Cambridge University Press: Cambridge, UK, 2015. [Google Scholar] [CrossRef]
  155. Franz, M.; Rozali, M. Mimicking black hole event horizons in atomic and solid-state systems. Nat. Rev. Mater. 2018, 3, 491–501. [Google Scholar] [CrossRef]
  156. Kolobov, V.I.; Golubkov, K.; Muñoz de Nova, J.R.; Steinhauer, J. Observation of stationary spontaneous Hawking radiation and the time evolution of an analogue black hole. Nat. Phys. 2021, 17, 362–367. [Google Scholar] [CrossRef]
  157. Sbitnev, V.I. Quaternion Algebra on 4D Superfluid Quantum Space-Time. Dirac’s Ghost Fermion Fields. Found. Phys. 2022, 52, 19. [Google Scholar] [CrossRef]
  158. Gallerati, A. Negative-curvature spacetime solutions for graphene. J. Phys. Condens. Matter 2021, 33, 135501. [Google Scholar] [CrossRef] [PubMed]
  159. Lambiase, G.; Papini, G. The Interaction of Spin with Gravity in Particle Physics; Springer Nature: Cham, Switzerland, 2021. [Google Scholar] [CrossRef]
  160. Clark, J.; Khodel, V.; Zverev, M.; Yakovenko, V. Unconventional superconductivity in two-dimensional electron systems with long-range correlations. Phys. Rep. 2004, 391, 123–156. [Google Scholar] [CrossRef]
  161. Uchihashi, T. Two-dimensional superconductors with atomic-scale thickness. Supercond. Sci. Technol. 2016, 30, 013002. [Google Scholar] [CrossRef]
Figure 1. SIS junction with an axis directed along the z-direction, parallel to the local Earth’s gravitational field.
Figure 1. SIS junction with an axis directed along the z-direction, parallel to the local Earth’s gravitational field.
Symmetry 14 00554 g001
Figure 2. Time dependence of the Josephson current for an insulating layer of thickness = 1 nm .
Figure 2. Time dependence of the Josephson current for an insulating layer of thickness = 1 nm .
Symmetry 14 00554 g002
Figure 3. Time dependence of the Josephson current for an insulating layer of thickness = 300 nm .
Figure 3. Time dependence of the Josephson current for an insulating layer of thickness = 300 nm .
Symmetry 14 00554 g003
Figure 4. Local gravitational field variation as a function of time for a In sample ( ξ 0 = 360 nm , T c = 3.410 K , Δ T = 10 3 K [118]) measured along the axis of a superconductive disk at fixed distance d = 0.1 cm above the base surface. The disk radius is R = 15 cm and the disk thickness is h = 3 cm .
Figure 4. Local gravitational field variation as a function of time for a In sample ( ξ 0 = 360 nm , T c = 3.410 K , Δ T = 10 3 K [118]) measured along the axis of a superconductive disk at fixed distance d = 0.1 cm above the base surface. The disk radius is R = 15 cm and the disk thickness is h = 3 cm .
Symmetry 14 00554 g004
Figure 5. Local gravitational field variation as a function of time for a Ba 0.4 K 0.6 Fe 2 As 2 sample ( ξ 0 = 1.20 nm , T c = 37.0 K , Δ T = 0.1 K [119]) measured along the axis of a superconductive disk at fixed distance d = 0.1 cm above the base surface. The disk radius is R = 15 cm and the disk thickness is h = 3 cm .
Figure 5. Local gravitational field variation as a function of time for a Ba 0.4 K 0.6 Fe 2 As 2 sample ( ξ 0 = 1.20 nm , T c = 37.0 K , Δ T = 0.1 K [119]) measured along the axis of a superconductive disk at fixed distance d = 0.1 cm above the base surface. The disk radius is R = 15 cm and the disk thickness is h = 3 cm .
Symmetry 14 00554 g005
Figure 6. Local gravitational field variation as a function of distance for the same In sample, measured along the disk axis above the base surface, at fixed time t = τ 0 = 1.64 ns . The disk radius is R = 15 cm and the disk thickness is h = 3 cm .
Figure 6. Local gravitational field variation as a function of distance for the same In sample, measured along the disk axis above the base surface, at fixed time t = τ 0 = 1.64 ns . The disk radius is R = 15 cm and the disk thickness is h = 3 cm .
Symmetry 14 00554 g006
Figure 7. Local gravitational field variation as a function of distance for the same Ba 0.4 K 0.6 Fe 2 As 2 sample, measured along the disk axis above the base surface at fixed time t = τ 0 = 7.50 × 10 3 ns . The disk radius is R = 15 cm and the disk thickness is h = 3 cm .
Figure 7. Local gravitational field variation as a function of distance for the same Ba 0.4 K 0.6 Fe 2 As 2 sample, measured along the disk axis above the base surface at fixed time t = τ 0 = 7.50 × 10 3 ns . The disk radius is R = 15 cm and the disk thickness is h = 3 cm .
Symmetry 14 00554 g007
Figure 8. Half-infinite superconductor approximation. The Earth’s gravitational field is parallel to the u x direction.
Figure 8. Half-infinite superconductor approximation. The Earth’s gravitational field is parallel to the u x direction.
Symmetry 14 00554 g008
Figure 9. Local field variation as a function of time for a Bi 2 Sr 2 Ca 3 Cu 3 O 10 sample ( T c 107 K , λ 0 2.4 × 10 7 m , σ 1 3.6 × 10 6 Ω m , D 10 3 m 2 / s , ξ 0 1 nm [134,135]) at a temperature T = 102 K . The red solid line refers to a sample of thickness L = 15 cm , while the blue dashed line shows the result for L = 5 cm .
Figure 9. Local field variation as a function of time for a Bi 2 Sr 2 Ca 3 Cu 3 O 10 sample ( T c 107 K , λ 0 2.4 × 10 7 m , σ 1 3.6 × 10 6 Ω m , D 10 3 m 2 / s , ξ 0 1 nm [134,135]) at a temperature T = 102 K . The red solid line refers to a sample of thickness L = 15 cm , while the blue dashed line shows the result for L = 5 cm .
Symmetry 14 00554 g009
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gallerati, A.; Ummarino, G.A. Superconductors and Gravity. Symmetry 2022, 14, 554. https://doi.org/10.3390/sym14030554

AMA Style

Gallerati A, Ummarino GA. Superconductors and Gravity. Symmetry. 2022; 14(3):554. https://doi.org/10.3390/sym14030554

Chicago/Turabian Style

Gallerati, Antonio, and Giovanni Alberto Ummarino. 2022. "Superconductors and Gravity" Symmetry 14, no. 3: 554. https://doi.org/10.3390/sym14030554

APA Style

Gallerati, A., & Ummarino, G. A. (2022). Superconductors and Gravity. Symmetry, 14(3), 554. https://doi.org/10.3390/sym14030554

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