In the following discussion, we assume that the quantum states are quasi-static. We solve the time-independent Schrödinger equation supposing that external potentials are independent of time . When we change them slowly enough, the corresponding quantum states (the effective wavefunction) changes so that they keep following the time-independent Schrödinger equation (with the given potentials) at each point of time.

First, we shows that the time-independent Schrödinger equation for the gauge-invariant wavefuntion has the following form:

--- (1)

**[Proof]**

We start from the time-dependent Schrödinger equation of the effective wafefunction:

--- (3.5)

When the potentials are static , the energy eigenstate is given by:

--- (2)

and (3.5) becomes:

--- (3)

From (2), the gauge-independent wavefunction has the following from:

So we define the time-independent part of as:

Equivalently, we have:

--- (4)

By substituting (4) into (3), we get (1). ■

Now we apply (1) to the one-dimensional potential barrier model of Josephson Junctions:

where .

For , the solution can be plane-wave functions:

The arbitrariness of the global phase is expressed by the constants . At this point, they can be taken independently, but later, we will show that they are related through the boundary condition of currents .

For , the solution can be a liner combination of and .

--- (5)

where .

From the boundary condition (continuity of ) at , we have:

By solving these linear equations, we have:

--- (3.37)

In general, the currents:

--- (3.23)

can be calculated from the gauge-invariant wavefunction as:

(You can confirm it with the direct calculation using the general form .)

In our particular case (5), we have:

Hence, we finally get:

--- (6)

where:

(6) suggests that the currents are constant and unique inside the insulator, and we enforce the boundary condition that it's the same as the currents outside the insulator:

--- (3.34)

From this condition, the phase difference can be decided.

In terms of the current intensity:

--- (3.32)

we have the similar relationship:

Now we assume that we slowly apply an external electric field over the circuit (outside the insulator) so that the voltage (the electric potential difference) between the contact points to the insulator becomes . By applying the discussion:

・Derivation of the gauge-independent relation between the phase and the electric field

the phase difference between two points follows the relationship:

--- (7)

(We're assuming that is very small and the change follows the quasi-static assumption described at the beginning of this note.)

So we conclude that we have the oscillating current inside the insulator:

#### Note

As in the discussion in the following link, I need an assumption that is constant and unique to use the relationship (7).

・Derivation of the gauge-independent relation between the phase and the electric field

So I applied the discussion to the wavefunction outside the insulator that (presumably) forms a closed loop. I'm not sure if this is really a common way to justify using (7) in this model.