Quantum Interferometry for Different Energy Landscapes in a Tuneable Josephson Junction Circuit

Pernel Nguenang^{1}^{*},
Michael Nana Jipdi^{1,2},
Patrick Louodop^{1,3},
Martin Tchoffo^{1,4},
Lukong Cornelius Fai^{1},
Hilda A. Cerdeira^{3}

Show more

1. Introduction

The Josephson effect, predicted in 1962 by B.D. Josephson [1], is a phenomenon that occurs when two superconductors are separated by a very thin layer of non-superconductive material. A superconducting current passes through the barrier and the electrical properties of this system are very precisely defined [2]. This current depends on the phase difference between the superconductors and it has been observed experimentally by P.L. Anderson and J.W Rowel [3] and called Josephson’s current [4]. Furthermore, if we apply a constant voltage to this junction, an alternating current is created and it was experimentally demonstrated by S. Saphiro in 1963 [5].

It has been demonstrated that the well-known Josephson effect is a result of the quantum tunnelling [2] [6] [7] [8] [9]. In fact, Josephson’s effect is a macroscopic observation of the tunnelling effect of condensed electron pairs between the Fermi surfaces of two metals through a thin barrier. From another perspective, quantum tunnelling is a well-known phenomenon where a particle crosses a potential barrier higher than its own energy [10]. There exist applications in many fields of physics, chemistry and biology [11] [12]. The quantum tunnelling phenomenon is the key point of quantum computing and has also been used in characterizing dynamical behaviour across Josephson junctions [13] [14] [15] [16]. The phenomenon gained a lot of interest in the two last decades [17] in the investigation of properties of superconducting devices [18] - [27]. One of the prominent applications of the Josephson junction is found in quantum computing where the resulting physical system when modelled as a qubit is helpful in implementing logic gates, as well as quantum transistors [28] [29]. This is achieved by modelling Josephson circuits as two or multi-state systems [30] [31]; some examples include traditional two-state systems [32], flux qubits [33], charge qubits [34], Cooper-pair box [35] [36] and three state systems [37]. In these set-ups, due to external sources (current or voltage), the tunnelling occurring is of the Landau Zener type.

In the early 1930’s, Landau and Zener developed a possible transition between two similar levels because of a control-value scanned at the point of minimum energy splitting [38] [39] [35]. The Landau-Zener tunneling is at the basis of several quantum mechanical processes, and it was recently discovered in periodic structures, with applications in driven superlattices [40], current driven Josephson junction [41] and bosonic systems [42] [43] [44]. The Landau-Zener tunneling occurs at the crossing of two energy levels which move away due to a weak interaction [35], and it has largely been investigated for particular configurations such as external field influence, periodic modulation [45] [46], non-linear models [42] [43] [47] [48], multi-state models [37] [49] and multi-crossing model in its multiple configurations (real or fictitious crossing) [50] [51] [52] and it still gains a lot of attention due to its prominent uses.

From another perspective, level crossing systems often show avoided energy level (anti-crossing) which can be handled using an external control-value [53]. When the energy level in a given region has a double crossing, the accumulation phase between transitions generates either a destructive or constructive interference in the time domain, called Stückelberg oscillations [54] [55]. When such a level crossing system is subjected to a periodic driving force, the physical observables of the system show a periodic dependence that mimics the Landau-Zener-Stückelberg interferometry [51]. Landau-Zener-Stückelberg interferometry is therefore the realization of an interferometer with an energy spectrum having at least two bands or branches, separated by a gap [56]. The Landau-Zener Stückelberg interferometry has been demonstrated for superconducting qubits [53] [57], a Cooper-pair box [35] [36], nitrogen deficiency centers [58], and quantum dots [59] [60]. The Landau-Zener-Stückelberg interferometry is of major importance in quantum computation since it helps in the design of controllable modulus and also confers a better read out for quantum logic gates [50] [61].

In the present work, we shall show that a simple electrical circuit with a single Josephson junction assisted by controllable parameters (inductance and capacitance) and under a periodic force will present the abovementioned phenomena. The main objective of this work is to show how the resulting modified Josephson junction can model either a controllable two or multi-state system where quantum tunnelling is induced by the Landau-Zener scenario. The resulting quantum interferometer is defined and the condition for constructive and destructive interferences is investigated.

The rest of this paper is organized as follows: Section 2 presents the electronic circuit and describes the elements. In Section 3, we focus on the construction of the Hamiltonian of the system based on the range of values of controllable parameters. Section 4 is devoted to quantum tunnelling and the evaluation of state occupation in different set-ups (two-state and multi-state); interference is also investigated. Section 5 presents the summary of the work.

2. Circuit and Model: Description of Circuit Elements

The electrical circuit that we are analyzing is given by the following graph (See Figure 1) as presented in [62].

Summarizing, a single-component Josephson junction (JJ) (blue part of the circuit [63] - [68] ) connected in parallel to a coil of inductance L, which is part of the C_{1}-L-C_{2} tank circuit series with a periodic force, while the Josephson junction,

Figure 1. Josephson junction based circuit [62].

the nonlinear element, plays the feedback loop [62]. Applying Kirchhoff’s law on the circuit on Figure 1, we obtain the following expressions:

$\frac{{\Phi}_{0}}{2\pi}\stackrel{\dot{}}{\varphi}\left(t\right)={V}_{j}\left(t\right)$ (1a)

$\begin{array}{c}{\stackrel{\dot{}}{V}}_{j}\left(t\right)=\frac{2C}{C+2{C}_{j}}\left(\frac{\stackrel{\dot{}}{e}\left(t\right)}{2}-\frac{{i}_{s}}{C}\mathrm{sin}\left(\varphi \right)-\frac{{V}_{j}\left(t\right)}{{R}_{j}C}\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-\frac{2C}{C+2{C}_{j}}\left(\frac{{\Phi}_{0}}{2\pi}\frac{1}{LC}\varphi \left(t\right)-p\right)\mathrm{,}\end{array}$ (1b)

where $\varphi $ and ${V}_{j}$ are respectively the phase and the voltage of the Josephson junction and the constant p is expressed as: $p=\frac{1}{LC}\varphi \left(0\right)-\frac{1}{C}i\left(0\right)$. The capacitors ${C}_{1}$ and ${C}_{2}$ are assumed identical and equal to C. $\frac{{\Phi}_{0}}{2\pi}$ is the quantum flux.

To facilitate discussions, we define new variables: ${x}_{1}=\varphi \left(t\right)$, ${P}_{1}\left(t\right)=\frac{{\Phi}_{0}}{2\pi}{C}_{j}{V}_{j}\left(t\right)$, ${e}_{0}$ is the amplitude of $\frac{\stackrel{\dot{}}{e}\left(t\right)}{2}$, $\tau =\Omega t$, $\beta =\frac{2C}{C+2{C}_{j}}\frac{{\Phi}_{0}}{2\pi}$, ${a}_{1}=\frac{{i}_{s}}{C}$, $\mu ={\left(\frac{{\Phi}_{0}}{2\pi}\right)}^{2}{C}_{j}$, ${a}_{3}=\frac{{\Phi}_{0}}{2\pi}\frac{1}{LC}$, such that the circuit represented in Figure 1 is well represented by the dimensionless generalized coordinates:

$\begin{array}{l}\mu {\stackrel{\dot{}}{x}}_{1}\left(\tau \right)={P}_{1}\left(\tau \right)\\ {\stackrel{\dot{}}{P}}_{1}\left(\tau \right)=\beta \left({e}_{0}\mathrm{sin}\left(\Omega \tau \right)-{a}_{1}\mathrm{sin}\left({x}_{1}\left(\tau \right)\right)-{a}_{3}{x}_{1}\left(\tau \right)-p\right).\end{array}$ (2)

Here, ${C}_{j}$ is the capacitance of the capacitor parallel to the Josephson junction component. It is the first controllable parameter that conditions the value of

the charging energy $\frac{{\hslash}^{2}}{2\mu}$ as well as the energy stored in the junction $\beta {a}_{1}$ through parameters $\beta $ and $\mu $.

L is the inductance of the coil parallel to the Josephson junction component; it is the second controllable parameter that conditions the shape of the potential energy and therefore, the level configuration as well as the number of states through the ${a}_{3}$ parameter.

The parameter p carries the initial conditions in the inductive branch. It is helpful in the design of constructive and destructive interferences. ${e}_{0}\mathrm{sin}\left(\tau \right)$ is the periodic force (external source) that induces the time dependence of energy levels and conditions the avoided crossing on the energy-time graph.

3. Hamiltonian Modelling

To study the microscopic phenomenon that conditions the behaviour of the circuit, we focus on Hermitian contributions and ignore dissipative terms (that means ${R}_{j}$ is very large). Representing Equation (3) as the Hamilton’s equation with generalized coordinates $\left({x}_{1}\mathrm{,}{P}_{1}\right)$, the Hamiltonian of the system is given by:

$H\left({x}_{1},{P}_{1}\right)=\frac{1}{2\mu}{P}_{1}{}^{2}+\beta \left(p-f\left(t\right)\right){x}_{1}+\beta {a}_{1}\mathrm{cos}\left({x}_{1}\right)+\beta {a}_{3}\frac{{x}_{1}^{2}}{2},$ (3)

and the Lagrangian

$L\left({x}_{1},{P}_{1}\right)=\frac{1}{2\mu}{P}_{1}{}^{2}-\beta \left(p-f\left(t\right)\right){x}_{1}-\beta {a}_{1}\mathrm{cos}\left({x}_{1}\right)-\beta {a}_{3}\frac{{x}_{1}^{2}}{2}.$ (4)

Here, ${x}_{1}$ plays the role of coordinates while ${P}_{1}$ is the particle momentun and from canonical quantization (first quantization) $\left[{x}_{1},{P}_{1}\right]=i\hslash $ and ${P}_{1}=-i\hslash \frac{\partial}{\partial {x}_{1}}$. From the Lagrangian of the system, we deduce in terms of ${x}_{1}$ the dimensionless confining static potential:

$U\left({x}_{1}\right)=\frac{2\mu}{{\hslash}^{2}}\left(\beta p-\beta f\left(t\right)\right){x}_{1}+\frac{2\mu}{{\hslash}^{2}}\beta {a}_{1}\mathrm{cos}\left({x}_{1}\right)+\frac{2\mu}{{\hslash}^{2}}\beta {a}_{3}\frac{{x}_{1}^{2}}{2}.$ (5)

The shape of the confining potential is conditioned by the tunable parameters (inductance L, capacitance C and ${C}_{j}$ ) through parameters $\beta $, ${a}_{1}$, $\mu $ and ${a}_{3}$ ; two particular cases are observed:

Firstly, the effective charging energy of the system is greater than the effective energy stored in the Josephson junction ( $\frac{2\mu}{{\hslash}^{2}}\beta {a}_{1}\ll 1$ and ${a}_{1}\gg {a}_{3}$ ). Here, the last two terms of the Hamiltonian act as a perturbation to the self consistent part.

In the second case, the effective charging energy of the system is smaller than the effective energy stored in the Josephson junction ( $\frac{2\mu}{{\hslash}^{2}}\beta {a}_{1}\gg 1$ and ${a}_{1}\gg {a}_{3}$ ).

Here, the last two terms of the Hamiltonian can no more be treated as a perturbation but as a part of the self consistent Hamiltonian. It is instructive to mention that the condition ${a}_{1}>{a}_{3}$ is very crucial for the appearance of multiple

potential wells; if ${a}_{1}\le {a}_{3}$ and no matter the value of $\frac{2\mu}{{\hslash}^{2}}\beta {a}_{1}$, we obtain a

modified single harmonic potential and therefore no tunnelling barrier. These situations are depicted in Figure 2 and Figure 3.

In Figure 2, we present the case where the effective charging energy of the system is greater than the effective energy stored in the Josephson junction

( $\frac{2\mu}{{\hslash}^{2}}\beta {a}_{1}\ll 1$ ). The diagram presents several minima with different energies;

each minimum represents a possible phase state and the system can be modelled as a multi-state one.

In Figure 3 ( $\frac{2\mu}{{\hslash}^{2}}\beta {a}_{1}\gg 1$ ), we find only two relevant energy minima; the

remaining curvature do not present any dept and cannot be considered as a confining well. The system can be modeled as a two-state system with different ground state energy [69].

It is instructive to mention that, the number of states in the system is conditioned

Figure 2. Schematic representation of the potential with respect to ${x}_{1}$. $\frac{2\mu}{{\hslash}^{2}}\beta =0.25$ ; ${a}_{1}=2$ ; $p=1$ ; ${a}_{3}=0.01$.

Figure 3. Schematic representation of the potential with respect to ${x}_{1}$. $\frac{2\mu}{{\hslash}^{2}}\beta =9$ ; ${a}_{1}=2$ ; $p=1$ ; ${a}_{3}=0.5$.

by the inductance value. From a deeper investigation, if $\frac{{a}_{3}}{{a}_{1}}>1$, the system mimics a simple harmonic potential, the system is devoid of potential barrier. If $0.1284374<\frac{{a}_{3}}{{a}_{1}}<1$ ( ${a}_{3}<{a}_{1}$ ), the system presents two energy minima; the supplementary curvature do not present any deepness and the confining potential is

a double well; this is helpful in designing a two-state system. If

$0.0709134<\frac{{a}_{3}}{{a}_{1}}<0.1284374$ ( ${a}_{3}<{a}_{1}$ ), the system has four energy minima and can therefore be modelled as a four-states system; $0.0490296<\frac{{a}_{3}}{{a}_{1}}<0.0709134$, we have six minima so a six-states system; $0.0374745<\frac{{a}_{3}}{{a}_{1}}<0.0490296$, the system presents eight-states; we can easily realize that the smaller the factor $\frac{{a}_{3}}{{a}_{1}}$,

the larger the number of energy minima and consequently the number of states. If that factor is exceedingly small, we have an infinite number of states describing the situation in Figure 2.

These particular behaviours are obtained from a single Josephson junction circuit only by playing on tuneable parameters that are the inductance L and the

capacitance C and ${C}_{j}$ ). Therefore, when $\frac{2\mu}{{\hslash}^{2}}\beta {a}_{1}<1$ we resume the condition as $\frac{4{C}_{j}}{C+2{C}_{j}}<\frac{{\hslash}^{2}}{{i}_{s}}{\left(\frac{2\pi}{{\Phi}_{0}}\right)}^{3}$ ; for ${a}_{3}<{a}_{1}$, we design as $\frac{{\varphi}_{0}}{2\pi}<L{i}_{S}$ ; here, the effective

charging energy of the system is larger than the effective energy stored in the junction. This condition is achieved by reducing the frequency $\Omega $ of generator for the process to be adiabatic and reducing the capacitance ${C}_{j}$ and enhancing

the inductance L. When the capacitance ${C}_{j}$ is enhanced, $\frac{4{C}_{j}}{C+2{C}_{j}}$ is large;

the effective energy stored in the junction becomes large and can no more be considered as a perturbation; the number of states of system is reduced and a

two-state system is achieved for $0.1284374<\frac{{a}_{3}}{{a}_{1}}<1$ with the confining potential represented by a double well (see Figure 3).

3.1. Large Effective Charging Energy

When the effective charging energy of the system is greater than the effective energy stored in the Josephson junction then $\frac{2\mu}{{\hslash}^{2}}\beta {a}_{1}\ll 1$ and ${a}_{1}\gg {a}_{3}$.

In this situation, the shape of the potential is given in Figure 2 and the Hamiltonian is given by

$H=\frac{1}{2\mu}{P}_{1}{}^{2}+\beta {a}_{1}\mathrm{cos}\left({x}_{1}\right)+\beta {a}_{3}\frac{{x}_{1}^{2}}{2}+\beta \left(p-f\left(t\right)\right){x}_{1}.$ (6)

Since $\frac{2\mu}{{\hslash}^{2}}\beta {a}_{1}\ll 1$ and ${a}_{1}\gg {a}_{3}$ then, the two last terms of the Hamiltonian are considered as a perturbation and the unperturbed part can be labeled as ${H}_{0}$ therefore,

$H={H}_{0}+{V}_{1},$

where

${H}_{0}=\frac{1}{2\mu}{P}_{1}{}^{2}+\beta {a}_{1}\mathrm{cos}\left({x}_{1}\right),$

${V}_{1}=\beta {a}_{3}\frac{{x}_{1}^{2}}{2}+\beta \left(p-f\left(t\right)\right){x}_{1}.$

To quantize this Hamiltonian, we use the displacement operator and write down the wave function as follows:

$\psi \left({x}_{1},t\right)=\phi \left({x}_{1},t\right)\mathrm{exp}\left[-\frac{i}{\hslash}{\displaystyle \int \beta \left(p-f\left(t\right)\right){x}_{1}\text{d}t}\right].$ (7)

The effective Hamiltonian of the perturbed system is reduced to:

$H=\frac{1}{2\mu}{\left({P}_{1}-F\left(t\right)\right)}^{2}+\beta {a}_{1}\mathrm{cos}\left({x}_{1}\right)+\beta {a}_{3}\frac{{x}_{1}^{2}}{2},$ (8)

with

$F\left(t\right)={\displaystyle \int \beta}\left(p-{e}_{0}\mathrm{sin}\left(\Omega t\right)\right)\text{d}t\mathrm{.}$ (9)

The Schrödinger equation of the unperturbed part is

$\frac{{\partial}^{2}{\psi}_{m}}{\partial {x}^{2}}+\left(\frac{2\mu {E}_{m}}{{\hslash}^{2}}-\frac{2\mu \beta {a}_{1}\mathrm{cos}\left(x\right)}{{\hslash}^{2}}\right){\psi}_{m}=0.$ (10)

Because of intrinsic parity of the potential, the solution is characterized by an even function. The solution of this equation is the well known periodic Mathieu function [70] given by

${\psi}_{m}=AC{e}_{{i}_{m}}\left(\frac{2\mu \beta {a}_{1}}{{\hslash}^{2}},x\right),$ (11)

or in terms of series

$C{e}_{{i}_{m}}={\displaystyle \underset{n=0}{\overset{\infty}{\sum}}}\text{\hspace{0.05em}}{A}_{n}^{\left(m\right)}\left(\frac{2\mu \beta {a}_{1}}{{\hslash}^{2}}\right)\mathrm{cos}\left(mx\right).$

The unperturbed wave function is then a linear combination of cosine functions; since $\frac{2\mu \beta {a}_{1}}{{\hslash}^{2}}\ll 1$, the energy of the unperturbed part is [71] [72].

$\frac{2\mu {E}_{m}}{{\hslash}^{2}}\sim {m}^{2}+\frac{{\left(2\mu \beta {a}_{1}\right)}^{2}}{2{\hslash}^{4}\left({m}^{2}-1\right)}$ (12)

or

${E}_{m}\sim \frac{{\hslash}^{2}}{2\mu}{m}^{2}+\frac{\mu {\left(\beta {a}_{1}\right)}^{2}}{{\hslash}^{2}\left({m}^{2}-1\right)}$.

The second part of the energy is just the perturbation (very small) and the system mimics a high energy quantum rotor. From the wave function in Equation (11), we can build the solution of the system inducing the perturbation ( ${V}_{1}$ ) as a linear combination of cosine functions but with variable coefficients. Therefore,

${\psi}_{m}={\displaystyle \underset{m}{\sum}}\frac{{A}_{m}\left(t\right)}{\sqrt{\pi}}\mathrm{cos}mx.$ (13)

This wave function satisfies the time dependent Schrödinger equation; which when multiplied by the factor $\frac{1}{\sqrt{\pi}}\mathrm{cos}nk$ and the integral over $0\prec x\prec 2\pi $ one has

$\begin{array}{c}i\hslash \frac{\partial}{\partial t}{\stackrel{\dot{}}{A}}_{n}\left(t\right)=\frac{\beta {a}_{1}}{2}\left({A}_{n-1}+{A}_{n+1}\right)+{\displaystyle \underset{m}{\sum}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\beta {a}_{3}\frac{1}{{\left(n-m\right)}^{2}}{A}_{n-m}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+{\displaystyle \underset{n}{\sum}}\frac{{\left(\hslash m-F\left(t\right)\right)}^{2}}{2\mu}{A}_{n}.\end{array}$ (14)

We observe that the additional coupling factor( $\beta {a}_{3}\frac{1}{{\left(n-m\right)}^{2}}$ ) decays quickly

when the difference between quantum numbers ( $m-n$ ) is large; since ${a}_{3}$ is chosen to be very small, only the nearest neighbour are relevant and enables us to write the effective Hamiltonian

$\mathcal{H}={\displaystyle \underset{m}{\sum}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{\u03f5}_{m}{A}_{m}^{\u2020}{A}_{m}+\Delta \left({A}_{m}^{\u2020}{A}_{m-1}+{A}_{m}^{\u2020}{A}_{m+1}\right),$ (15)

or in terms of states:

$\mathcal{H}={\displaystyle \underset{m}{\sum}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{\u03f5}_{m}|m\rangle \langle m|+\Delta \left(|m\rangle \langle m-1|+|m\rangle \langle m+1|\right)$ (16)

with

${\epsilon}_{m}=\frac{4{\pi}^{2}}{3}\beta {a}_{3}+\frac{1}{2\mu}{\left(\hslash m-F\left(t\right)\right)}^{2},$

$\Delta =\frac{\beta {a}_{1}}{2}+\frac{\beta {a}_{3}}{2}.$

Focusing only on states $m=0$, $m=1$ and $m=-1$, the electronic circuit mimics a three-state system:

$\mathcal{H}=\left(\begin{array}{ccc}{\u03f5}_{-1}& \Delta & 0\\ \Delta & {\u03f5}_{0}& \Delta \\ 0& \Delta & {\u03f5}_{1}\end{array}\right).$ (17)

${\u03f5}_{-1}$, ${\u03f5}_{0}$ and ${\u03f5}_{1}$ are the energies of the three relevant states; $\Delta $ stands for the energy gap.

If all the states are considered as relevant, the system mimics a multi-state system with an infinite number of states represented by the Hamiltonian as follow:

$\mathcal{H}={\displaystyle \underset{m}{\sum}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{\u03f5}_{m}|m\rangle \langle m\mathrm{|}+\Delta \left(|m\rangle \langle m-\mathrm{1|}+|m\rangle \langle m+1|\right)$

3.2. Low Effective Charging Energy

In the case where the effective charging energy of the system is less than the effective energy stored in the Josephson junction, the factor $\frac{2\mu}{{\hslash}^{2}}\beta {a}_{3}$ is large and

the shape of the confining potential changes drastically (double well potential Figure 3). This is achieved by reducing the value of the capacitance C and the inductance L. The Hamiltonian of the system is given as in Equation (8)

To analyse the form of the wave function, we take

${H}_{0}=\frac{1}{2\mu}{P}_{1}{}^{2}+\beta {a}_{1}\mathrm{cos}\left({x}_{1}\right)\mathrm{,}$ (18)

the corresponding eigenfunctions of this part are periodic Mathieu wave function:

$C{e}_{i}{}_{m}={\displaystyle \underset{m}{\overset{\infty}{\sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{A}_{n}^{m}\left(q\right)\mathrm{cos}\left(mx\right)\mathrm{,}$

but in this configuration $\frac{2\mu}{{\hslash}^{2}}\beta {a}_{1}\gg 1$ therefore, the energy spectrum of the system [71] [72]:

$\frac{2\mu {\u03f5}_{m}}{{\hslash}^{2}}\simeq -\left(\frac{4\mu \beta {a}_{1}}{{\hslash}^{2}}\right)+\left(2n+1\right){\left(\frac{2\mu \beta {a}_{1}}{{\hslash}^{2}}\right)}^{\frac{1}{2}}\mathrm{,}$

so

${\u03f5}_{m}\simeq 2\beta {a}_{1}+\hslash {\Omega}_{\beta}\left(n+\frac{1}{2}\right)\mathrm{,}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\Omega}_{\beta}={\left(\frac{\beta {a}_{1}}{8\mu}\right)}^{\frac{1}{2}}\mathrm{.}$ (19)

The latter corresponds approximatively to the solution of a particle in a harmonic potential. The eigenfunction of the whole system can then be described by the Hermite polynomial

${\varphi}_{n}\left(x\right)\sim {\left(\frac{\beta}{{2}^{n}n\mathrm{!}\sqrt{\pi}}\right)}^{\frac{1}{2}}{H}_{n}\left(\beta n\right)\mathrm{exp}\left(\frac{-{\beta}^{2}{x}^{2}}{2}\right)\mathrm{.}$ (20)

Since the junction operates at very low temperature, only the ground state ( $n=0$ ) is relevant. Observing Figure 3, we have two states corresponding to the ground state of each potential well and described by the left and right wave functions

${\varphi}_{R}\left(y\right)={\left(\frac{k}{\sqrt{\pi}}\right)}^{\frac{1}{2}}\mathrm{exp}\left(-\frac{{k}^{2}}{2}{y}^{2}\right),$ (21)

${\varphi}_{L}\left(z\right)={\left(\frac{k}{\sqrt{\pi}}\right)}^{\frac{1}{2}}\mathrm{exp}\left(-\frac{{k}^{2}}{2}{z}^{2}\right),$ (22)

with $k={\left(\frac{\mu {\Omega}_{\alpha}}{\hslash}\right)}^{\frac{1}{2}}$ and ${\Omega}_{\alpha}^{2}=\frac{\beta}{\mu}\left({a}_{3}-\sqrt{{a}_{1}^{2}-{a}_{3}^{2}{\alpha}^{2}}\right)$. $y={x}_{1}+\alpha $ and $z={x}_{1}-\alpha $ are two new variables that help to center the wave function with the bottoms of the double well. The bottoms are located at ${x}_{1}=\pm \alpha $.

The global wave function of our system is now written as a superposition of all possible states:

$\psi ={C}_{R}|{\psi}_{R}\left(y\right)\rangle +{C}_{L}|{\psi}_{L}\left(z\right)\rangle \mathrm{.}$ (23)

The wave function (24) can now be used and the time dependent Schrödinger equation for the perturbed system gives:

$\begin{array}{l}i\hslash \frac{\text{d}{C}_{L}}{\text{d}t}={\u03f5}_{L}\left(t\right){C}_{L}+{\Delta}_{0}{C}_{R}\mathrm{,}\\ i\hslash \frac{\text{d}{C}_{R}}{\text{d}t}={\u03f5}_{R}\left(t\right){C}_{R}+{\Delta}_{0}{C}_{L}\mathrm{,}\end{array}$ (24)

with

${\u03f5}_{L}={\u03f5}_{0}-\alpha \beta \left(p-f\left(t\right)\right)\mathrm{,}$

${\u03f5}_{R}={\u03f5}_{0}+\alpha \beta \left(p-f\left(t\right)\right)\mathrm{,}$

${\Delta}_{0}=\mathrm{exp}\left(-k{\alpha}^{2}\right)\left(\frac{\beta {a}_{3}}{2k}+\beta {a}_{1}\mathrm{exp}\left(-\frac{1}{4k}\right)+\frac{{\hslash}^{2}}{4\mu}\left(1-2k{\alpha}^{2}\right)\right),$

where $\alpha $ is the solution of the nonlinear equation ${a}_{1}\mathrm{sin}\left(\alpha \right)-{a}_{3}\alpha =0$ ; the term $2\alpha $ is the width of the barrier between the two wells. ${\u03f5}_{0}$ is a constant energy shift with no relevance as far as probability is concerned.

The Hamiltonian of the system is then written as follows:

${\mathcal{H}}_{eff}={\u03f5}_{L}\left(t\right){C}_{L}^{\u2020}{C}_{L}+{\u03f5}_{R}\left(t\right){C}_{R}^{\u2020}{C}_{R}+{\Delta}_{0}\left(\alpha \right){C}_{R}^{\u2020}{C}_{L}+{\Delta}_{0}\left(\alpha \right){C}_{L}^{\u2020}{C}_{R}\mathrm{.}$ (25)

This is the Hamiltonian of a two-state system that can also be rewritten in terms of Pauli matrix

${\mathcal{H}}_{eff}=\alpha \beta \left(p-f\left(t\right)\right){\sigma}_{z}+{\Delta}_{0}{\sigma}_{x}+{\u03f5}_{0}\mathrm{.}$ (26)

The system mimics the flux qubit or Cooper-pair box with two relevant states [33] [35] [36].

4. Quantum Tunneling: Survival and Transition Probabilities

In order to investigate the quantum tunneling occurring in the system, we focus on the two models mentioned above. For each of these configurations, the system is isolated and therefore the study related to decoherence are omitted; the transition probabilities will be investigated for each cases.

4.1. Two-State

The two-state system is obtained when $\frac{2\mu}{{\hslash}^{2}}\beta {a}_{1}\gg 1$ and ${a}_{1}\gg {a}_{3}$. To achieve this, we reduce the capacitance ${C}_{j}$ as well as the inductance L so that $\frac{4{C}_{j}}{C+2{C}_{j}}\frac{{i}_{s}}{{\hslash}^{2}}{\left(\frac{{\Phi}_{0}}{2\pi}\right)}^{2}\gg 1$. The Hamiltonian of the system is given by

$\begin{array}{c}{\mathcal{H}}_{eff}=\left({\u03f5}_{0}-\alpha \beta \left(f\left(t\right)-p\right)\right){C}_{L}^{\u2020}{C}_{L}+{\Delta}_{0}{C}_{L}^{\u2020}{C}_{R}+{\Delta}_{0}{C}_{R}^{\u2020}{C}_{L}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\left({\u03f5}_{0}+\alpha \beta \left(f\left(t\right)-p\right)\right){C}_{R}^{\u2020}{C}_{R}\mathrm{.}\end{array}$

To solve the problem, we make use of the dynamic matrix approach [73]. We introduce the field vector of the system as $\psi =\left(\begin{array}{c}{C}_{L}|L\rangle \\ {C}_{R}|R\rangle \end{array}\right)$ with ${\psi}_{0}=\left(\begin{array}{c}|L\rangle \\ |R\rangle \end{array}\right)$

being the static field function; we construct the Hamiltonian of the system in matrix representation as

$H={\psi}^{+}\epsilon \left(t\right)\psi +{\psi}^{+}G\psi \mathrm{,}$ (27)

where, $\epsilon \left(t\right)=\left(\begin{array}{cc}{\u03f5}_{0}+\alpha \beta \left(p-f\left(t\right)\right)& 0\\ 0& {\u03f5}_{0}-\alpha \beta \left(p-f\left(t\right)\right)\end{array}\right)$ and $G=\left(\begin{array}{cc}0& {\Delta}_{0}\left(\alpha \right)\\ {\Delta}_{0}\left(\alpha \right)& 0\end{array}\right)$ are respectively the energy and tunneling matrix. The equation of motion is given as follows:

$i\hslash {\psi}_{\mathrm{,0}}=\epsilon \left(t\right)\psi +G\psi \mathrm{.}$ (28)

Considering the transformation

$\psi \to \mathrm{exp}\left(-\frac{i}{\hslash}{\displaystyle \underset{0}{\overset{t}{\int}}\epsilon \left(\tau \right)\text{d}\tau}\right)\psi \mathrm{,}$ (29)

and with the help of the Campbell-Baker-Hausdorff Formula [74], the equation of motion expressed through the dynamic matrix is as follows

${\psi}_{\mathrm{,0}}=\frac{1}{i\hslash}\mathrm{exp}\left(\frac{i}{\hslash}{\displaystyle \underset{0}{\overset{t}{\int}}\epsilon \left(\tau \right)\text{d}\tau}\right)G\mathrm{exp}\left(-\frac{i}{\hslash}{\displaystyle \underset{0}{\overset{t}{\int}}\epsilon \left(\tau \right)\text{d}\tau}\right)\psi \mathrm{,}$ (30)

where

${\psi}_{\mathrm{,0}}=\frac{\text{d}\xi}{\text{d}t}\psi \mathrm{,}$ (31)

with the dynamic matrix being

$\xi =\left(\begin{array}{cc}0& a\\ -{a}^{\mathrm{*}}& 0\end{array}\right)\mathrm{,}$ (32)

and

$a={\displaystyle \underset{-\infty}{\overset{t}{\int}}}\frac{{\Delta}_{0}\left(\alpha \right)}{i\hslash}\mathrm{exp}\left(\frac{2i\alpha \beta}{\hslash}{\displaystyle \int}\left(p-f\left(\tau \right)\right)\text{d}\tau \right)\text{d}\tau ,$ (33)

We find the properties and exact solution of the transition matrix ${T}^{M}$ of this problem. The analytical expression of ${T}^{M}$ is obtained from the solution of:

$\psi ={T}^{M}{\psi}_{0}=\stackrel{^}{T}\text{\hspace{0.05em}}\mathrm{exp}\left(\xi \right){\psi}_{0}.$ (34)

where
$\stackrel{^}{T}$ is the time ordering operator. We remark that the Hamiltonian describes a periodic system. Therefore, if the transition matrix
${T}^{M}$ describes the evolution of the system after a single period, the evolution after n periods is obtained as the n^{th} power of the transition matrix. Therefore, if T is the period and n an arbitrary integer,

$\psi \left(T\right)={T}^{M}{\psi}_{0},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\psi \left(nT\right)={\left({T}^{M}\right)}^{n}{\psi}_{0},$ (35)

$\psi \left(nT+\tau \right)={\left({T}^{M}\right)}^{n}\psi \left(\tau \right),$ (36)

Therefore, the main difficulty is to find out ${T}^{M}$ for a single period. To analyse that term, we plot the diabatic and adiabatic energies of the system in Figure 4.

It is seen that we have two crossings, the first playing the role of a beam splitter while the second recombine the energy lines giving rise to an interferometer. Such a system has been largely studied in the literature [53] [57].

The time interval two crossings is known as relaxation time and conditions the appearance of interferences. For small relaxation time, the system undergoes a diabatic change resulting in a sudden death of the survival probability at the effective crossing point. The system mimics a single passage problem that resembles the so-called slow driving regime [51]. On the other hand, for large relaxation time, the system undergoes an adiabatic change and saturates prior to the second crossing thereby guaranteeing a survival probability. Here, the dynamic matrix of the system can be decomposed and reordered in term of crossing time; the later mimics the so-called fast driving regime see ref. [51].

Therefore, the solution (34) can be reordered in terms of the two crossing times. The field vector of the system is given by

Figure 4. Diabatic and adiabatic energies versus dimensionless time for a two-state configuration. The dotted black line is diabatic energy and the solid blue line is adiabatic energies.

$\psi =\stackrel{^}{T}\mathrm{exp}\left({\xi}_{1}\left(t\right)\right)\mathrm{exp}\left({\xi}_{2}\left(t\right)\right){\psi}_{0},$ (37)

and the transition matrix

${T}^{M}=\stackrel{^}{T}\mathrm{exp}\left({\xi}_{1}\left(t\right)\right)\mathrm{exp}\left({\xi}_{2}\left(t\right)\right),$ (38)

with

${\xi}_{1}\left(t\right)=\left(\begin{array}{cc}0& a\\ {a}^{*}& 0\end{array}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\xi}_{2}\left(t\right)=\left(\begin{array}{cc}0& b\\ {b}^{*}& 0\end{array}\right),$

$a={\displaystyle \underset{-\infty}{\overset{t}{\int}}}\frac{{\Delta}_{0}(\alpha )}{i\hslash}\mathrm{exp}\left(\frac{2i\alpha \beta}{\hslash}{\displaystyle \int}\left(p-f(\tau )\right)\text{d}\tau \right)\text{d}\tau ,$

$b={\displaystyle \underset{-\infty}{\overset{t}{\int}}}\frac{{\Delta}_{0}\left(\alpha \right)}{i\hslash}\mathrm{exp}\left(\frac{2i\alpha \beta p\pi}{\hslash \Omega}\right)\mathrm{exp}\left(\frac{2i\alpha \beta}{\hslash}{\displaystyle \int}\left(p-f(\tau )\right)\text{d}\tau \right)\text{d}\tau .$

The action of the time ordering operator help in the derivation of the transition matrix of the system [50] [75]; it is a $2\times 2$ matrix with elements:

${\left({T}^{M}\right)}_{11}={p}_{1}{p}_{2}-{q}_{1}{q}_{2}\mathrm{exp}\left(-2i{\Phi}_{\text{ac}}\right)\mathrm{,}$ (39)

${\left({T}^{M}\right)}_{12}=i\left({p}_{1}{q}_{2}\mathrm{exp}\left(i{\Phi}_{\text{ac}}\right)-{q}_{1}{p}_{2}\mathrm{exp}\left(-i{\Phi}_{\text{ac}}\right)\right)\mathrm{exp}\left(-i{\Phi}_{0}\right)\mathrm{,}$ (40)

${\left({T}^{M}\right)}_{21}=-{\left({\left({T}^{M}\right)}_{12}\right)}^{*},{\left({T}^{M}\right)}_{11}={\left({\left({T}^{M}\right)}_{22}\right)}^{*},$ (41)

with

${q}_{k}=i\sqrt{1-\mathrm{exp}\left(-2\pi {\lambda}_{\text{\hspace{0.05em}}v}\Lambda \left({\left(-1\right)}^{k}v\mathrm{,}\frac{2\alpha \beta {e}_{0}}{\hslash \Omega}\right)\right)}\text{\hspace{0.05em}}\mathrm{exp}\left(i{\varphi}_{\text{ks}}\right)\mathrm{,}$

${p}_{k}=\mathrm{exp}\left(-\pi {\lambda}_{\text{\hspace{0.05em}}v}\Lambda \left({\left(-1\right)}^{k}v\mathrm{,}\frac{2\alpha \beta {e}_{0}}{\hslash \Omega}\right)\right)\mathrm{,}k=\mathrm{1,2,}$

$v=\frac{2\alpha \beta p}{\hslash \Omega}\mathrm{,}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.17em}}2{\Phi}_{0}={\varphi}_{\text{\hspace{0.05em}}\text{2s}}+{\varphi}_{\text{\hspace{0.05em}}\text{1s}}+\pi v\mathrm{,}$

here,
$2{\Phi}_{\text{ac}}={\varphi}_{\text{\hspace{0.05em}}\text{2s}}-{\varphi}_{\text{\hspace{0.05em}}\text{1s}}+\pi v$ is the phase accumulated after the two passages,
${\varphi}_{\text{\hspace{0.05em}}\text{ks}}=\lambda \mathrm{ln}\lambda +\lambda +\Gamma \left(i\lambda \right)$ is the Stock’s phase of the k^{th} passage

${\lambda}_{v}=\frac{{\Delta}_{0}^{2}\left(\alpha \right)}{\hslash \Omega \text{\hspace{0.05em}}\alpha \beta {e}_{0}},\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\lambda ={\lambda}_{v}\Lambda \left(v,\frac{2\alpha \beta {e}_{0}}{\hslash \Omega}\right).$ (42)

are respectively the Landau-Zener and the renormalized Landau-Zener parameters of the system,

$\Lambda \left(v,\frac{2\alpha \beta {e}_{0}}{\hslash \Omega}\right)=\frac{\pi \alpha \beta {e}_{0}}{\hslash \Omega}{A}_{J}^{2}\left(v,\frac{2\alpha \beta {e}_{0}}{\hslash \Omega}\right)+\frac{\pi \alpha \beta {e}_{0}}{\hslash \Omega}{W}_{E}^{2}\left(v,\frac{2\alpha \beta {e}_{0}}{\hslash \Omega}\right)$ (43)

is the adiabatic power of the system that renormalizes the Landau-Zener parameter and appears as a consequence of the nonlinearity of energy level.

Under these considerations, the survival and transition probabilities of the system are respectively given by

${P}_{s}={\left({p}_{1}{p}_{2}-{q}_{1}{q}_{2}\mathrm{cos}\left(2{\Phi}_{\text{ac}}\right)\right)}^{2}+{q}_{2}^{2}{q}_{1}^{2}{\mathrm{sin}}^{2}\left(2{\Phi}_{\text{ac}}\right)\mathrm{,}$ (44)

${P}_{\text{tr}}={q}_{1}^{2}{p}_{2}^{2}+{p}_{1}^{2}{q}_{2}^{2}+2{p}_{1}{p}_{2}{q}_{1}{q}_{2}\mathrm{cos}\left(2{\Phi}_{\text{ac}}\right)\mathrm{.}$ (45)

It is instructive to mention that the probability as well as the phase accumulated is different from one passage to another. In the particular case where the system is design in such a way that ${p}_{1}={p}_{2}=p$ then the survival and transition probabilities become

${P}_{s}={\left(2{p}^{2}-1\right)}^{2}+4{p}^{2}\left(1-{p}^{2}\right){\mathrm{cos}}^{2}\left({\Phi}_{\text{st}}\right)\mathrm{,}$ (46)

${P}_{\text{tr}}=4{p}^{2}\left(1-{p}^{2}\right){\mathrm{sin}}^{2}\left({\Phi}_{\text{st}}\right)\mathrm{,}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.17em}}{\Phi}_{\text{st}}={\Phi}_{\text{ac}}+\frac{\pi}{2}\mathrm{.}$ (47)

The system under this condition mimics a Landau-Zener-Stückelberg interferometer with the accumulated phase being the Stückelberg phase [51]. For the

critical value $p=\frac{\sqrt{2}}{2}$, any anti-crossing behaves like a 50 - 50 beam splitter and the probabilities are given by

${P}_{s}={\mathrm{cos}}^{2}\left({\Phi}_{\text{st}}\right)\mathrm{,}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.17em}}{P}_{\text{tr}}={\mathrm{sin}}^{2}\left({\Phi}_{\text{st}}\right)\mathrm{.}$ (48)

the system oscillates between the two basis states and the interference in this case is constructive. In the particular case where $v=2n$ (n being an integer), the two passages are identical and therefore for small values of the gap

${\varphi}_{\text{\hspace{0.05em}}1\text{s}}-{\varphi}_{\text{\hspace{0.05em}}2\text{s}}=0$ (49)

and the probabilities are given by

${P}_{\text{tr}}=4\mathrm{exp}\left(-2\pi {\lambda}_{2n}\Lambda \left(2n\mathrm{,}\frac{2\alpha \beta {e}_{0}}{\hslash \Omega}\right)\right)\left(1-\mathrm{exp}\left(-2\pi {\lambda}_{2n}\Lambda \left(2n\mathrm{,}\frac{2\alpha \beta {e}_{0}}{\hslash \Omega}\right)\right)\right)$ (50)

${P}_{s}={\left(2\mathrm{exp}\left(-2\pi {\lambda}_{2n}\Lambda \left(2n\mathrm{,}\frac{2\alpha \beta {e}_{0}}{\hslash \Omega}\right)\right)-1\right)}^{2}\mathrm{.}$

The latter mimics the result derived by Zener in the double passage problem [76]. In a general frame, constructive interference will be observed when the imaginary part of the diagonal elements of the transition matrix tends to zero. Therefore the condition of appearance of constructive interference is resumed by

${\varphi}_{\text{2s}}-{\varphi}_{\text{1s}}+\pi v=\pi n\mathrm{,}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.17em}}n=\mathrm{0,1,2,3,4},\cdots $ (51)

Furthermore, the parameter v helps to control the style of interference. If v is an integer, the transitions are identical and the conditions above are almost satisfied. The plotting of survival and transition probabilities for constructive and destructive interferences is depicted on Figure 5 and Figure 6.

Figure 5 and Figure 6 show the plotting of survival and transition probabilities versus time for several periods and different values of parameters. From Figure 5, it is seen that after any crossing, the probability add up and enable the oscillation of probabilities between the two basis states (Stückelberg oscillations) [54] [55]. The observation in Figure 5 show that, after any half-period even if two consecutive crossings can add up there will be a cancellation that destroys the construction of probabilities; the system cannot oscillate between the two basis states and the interference is known as destructive. It can be seen that if v is an integer, the transitions under the same period are identical ( ${p}_{1}={p}_{2}$ ) and the phase accumulated become small. We register an oscillation of probabilities (see panel 5(a)); for some particular values of v, the second crossing behaves as a reflector but even in those cases the overall probability oscillates and the interference is still constructive. From another perspective, the types of interferences are conditioned by the parameter v; if v is an even integer, the interferences are purely constructive; the probabilities add up and we register an oscillation of probabilities. If v is an odd integer, the phase accumulated has a $\left(2n+1\right)\pi $ dephasing after a single passage. As a result, we register a certain sending back or cancelation after any construction, however, the overall probability oscillates around the two basis states with interferences being constructive. In the case where v is not an integer, the interferences are almost destructive unless condition (51) is satisfied. So the appearances of constructive and destructive interferences are controlled with the turnable parameters (p, L and ${C}_{j}$ ).

From another perspective, if the parameter v is larger than the oscillation amplitude ( $p>{e}_{0}$ ) there is no crossing in the energy spectrum and therefore, no tunneling. A similar observation is achieved when we reduce the value of the inductance (L). Since the width of the barrier ( $2\alpha $ ) is an inverse proportion of the inductance, any enhancement in the inductance reduces the barrier width. The plotting of the probability limit for a single transition against the width of the barrier is depicted in Figure 7.

Figure 7 shows the variation of the probability limits for a single passage versus the width of the barrier; the observation of this figure shows that for small values of $\alpha $ the transition probability is unity while the survival probability is zero. This only highlights that by reducing the width of the barrier (enhancing

Figure 5. Survival Probability (solid blue line) and transition probability (dotted-black line) versus dimensionless time ( $\tau $ ) for a two-state configuration (constructive interference) for parameters: (a) $v=0$, $\frac{2\alpha \beta {e}_{0}}{\hslash \Omega}=\frac{180}{\pi}$ ; (b) $v=3$, $\frac{2\alpha \beta {e}_{0}}{\hslash \Omega}=\frac{180}{\pi}$.

Figure 6. Survival probability (solid blue line) and transition probability (dotted-black line) versus dimensionless time ( $\tau $ ) for a two-state configuration (destructive interference) for parameters: (a) $v=1.25$, $\frac{2\alpha \beta {e}_{0}}{\hslash \Omega}=\frac{180}{\pi}$ ; (b) $v=3.125$, $\frac{2\alpha \beta {e}_{0}}{\hslash \Omega}=\frac{180}{\pi}$.

L), we easily achieve a perfect transfer of the population. However, for large values of $\alpha $ the transition probability vanishes while the survival probability is unity. This materializes the complete blockage of the system achieved for large values of $\alpha $. It is instructive to mention that, in this model, the total population

Figure 7. Survival probability (solid blue line) and transition probability (dotted-black line) versus width of the barrier ( $\alpha $ ) in a two-state configuration.

transfer as well as the complete blockage are obtained in permissible range of parameters only by changing the values of the inductance.

4.2. Three-State System

This configuration is obtained under the condition $\frac{2\mu}{{\hslash}^{2}}\beta {a}_{1}\ll 1$ but with ${a}_{1}\gg {a}_{3}$ where only the two relevant neighbouring states are considered. The Hamiltonian of the model is given by a $3\times 3$ matrix:

$\mathcal{H}=\left(\begin{array}{ccc}{\u03f5}_{-1}& {\Delta}_{0}& 0\\ {\Delta}_{0}& {\u03f5}_{0}& {\Delta}_{0}\\ 0& {\Delta}_{0}& {\u03f5}_{1}\end{array}\right).$

with

$\begin{array}{l}{\u03f5}_{m}\sim \frac{\hslash m}{\mu}\left(\alpha \beta pt-\frac{\alpha \beta {e}_{0}}{\Omega}\mathrm{cos}\left(\Omega t\right)\right)+\frac{\hslash}{2\mu}{\left(\alpha \beta pt-\frac{\alpha \beta {e}_{0}}{\Omega}\mathrm{cos}\left(\Omega t\right)\right)}^{2}\mathrm{.}\hfill \end{array}$

It can be seen that the energies of the system have both linear and a sinusoidal time dependence. The graphical representation of the diabatic and adiabatic energies is given in Figure 8.

This figure shows the plots of diabatic energies (dot black curves) as well as adiabatic (solid blue curves) as a function of the dimensionless time; the parameter

that makes time dimensionless is: $\tau =\sqrt{\frac{{\beta}_{0}}{\hslash}}t$ ( ${\beta}_{0}=\frac{\hslash \alpha \beta p}{\mu}$ ). The system

Figure 8. Diabatic and adiabatic energies versus dimensionless time. The dotted black line is diabatic energy and the solid blue line is adiabatic energies.

presents a single crossing and the probability of tunnelling through avoided crossing is the quantity of primary interests. This figure also shows the energy to be non-linear. The non-linearity may change the transition probability limit. When the system is driven slowly, it follows an adiabatic scenario; the resulting level crossing system is resumed by a three-state Landau-Zener problem. The investigation of the tunneling probability is done via dynamic matrix approach [52] [75]. In order to find the transition probability, we use the Hamiltonian (17).

Considering ${\psi}_{0}=\left(\begin{array}{c}|-1\rangle \\ |0\rangle \\ |1\rangle \end{array}\right)$ being the static field function then, the energy and tunneling matrix are respectively given by

$\epsilon \left(t\right)=\left(\begin{array}{ccc}{\u03f5}_{-1}& 0& 0\\ 0& {\u03f5}_{0}& 0\\ 0& 0& {\u03f5}_{1}\end{array}\right)$ and $G=\left(\begin{array}{ccc}0& {\Delta}_{0}& 0\\ {\Delta}_{0}& 0& {\Delta}_{0}\\ 0& {\Delta}_{0}& 0\end{array}\right)$

the master equation becomes

${\psi}_{\mathrm{,0}}=\frac{\text{d}\xi}{\text{d}t}\psi \mathrm{,}$ (52)

with the dynamic matrix being

$\xi =\left(\begin{array}{ccc}0& a& 0\\ -{a}^{\mathrm{*}}& 0& b\\ 0& -{b}^{\mathrm{*}}& 0\end{array}\right)\mathrm{,}$ (53)

where

$a={\displaystyle \underset{-\infty}{\overset{t}{\int}}}\frac{\Delta}{i\hslash}{\displaystyle \underset{n=-\infty}{\overset{n=+\infty}{\sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{\text{J}}_{n}\left(z\right)\mathrm{exp}\left(\frac{i{\beta}_{0}{\tau}^{2}}{2\hslash}+{v}_{1n}\tau \right)\text{d}\tau \mathrm{,}$

$b={\displaystyle \underset{-\infty}{\overset{t}{\int}}}\frac{\Delta}{i\hslash}{\displaystyle \underset{n=-\infty}{\overset{n=+\infty}{\sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{\text{J}}_{n}\left(z\right)\mathrm{exp}\left(\frac{i{\beta}_{0}{\tau}^{2}}{2\hslash}+{v}_{2n}\tau \right)\text{d}\tau \mathrm{,}$

${\beta}_{0}=\frac{\hslash \alpha \beta p}{\mu}\mathrm{,}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{v}_{kn}={\left(-1\right)}^{k}\frac{{\hslash}^{2}}{2\mu}+n\hslash \Omega \mathrm{,}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}z=\frac{\alpha \beta {e}_{0}}{\mu \text{\hspace{0.05em}}{\Omega}^{2}}\mathrm{.}$

The analytical expression of ${T}^{M}$ is obtained from the solution of:

$\psi ={T}^{M}{\psi}_{0}=\stackrel{^}{T}\text{\hspace{0.05em}}\mathrm{exp}\left(\xi \right){\psi}_{0}\mathrm{,}$ (54)

where $\stackrel{^}{T}$ is the time ordering operator. With the help of Wick’s theorem and Dyson series [77] [78] [79] [80], we compute the latter equation and obtain

$\psi =\left[\begin{array}{ccc}{T}_{11}^{M}& {T}_{12}^{M}& {T}_{13}^{M}\\ {T}_{21}^{M}& {T}_{22}^{M}& {T}_{23}^{M}\\ {T}_{31}^{M}& {T}_{32}^{M}& {T}_{33}^{M}\end{array}\right]{\psi}_{0}\mathrm{,}$ (55)

where

${T}_{11}^{M}=1-\frac{{\lambda}_{00}}{{\lambda}_{00}+{\lambda}_{1}}\left(1-\mathrm{exp}\left(-\pi \left({\lambda}_{00}+{\lambda}_{01}\right)\right)\right)\mathrm{,}$

${T}_{12}^{M}={T}_{21}^{M\ast}=-i{\left(\frac{{\lambda}_{00}}{{\lambda}_{00}+{\lambda}_{01}}\right)}^{\frac{1}{2}}\sqrt{1-\mathrm{exp}\left(-2\pi \left({\lambda}_{00}+{\lambda}_{01}\right)\right)}\mathrm{,}$

${T}_{13}^{M}={T}_{31}^{M\ast}=\frac{\sqrt{{\lambda}_{00}{\lambda}_{01}}}{{\lambda}_{00}+{\lambda}_{01}}\left(1-\mathrm{exp}\left(-\pi \left({\lambda}_{00}+{\lambda}_{01}\right)\right)\right)$

${T}_{22}^{M}=\mathrm{exp}\left(-\pi \left({\lambda}_{00}+{\lambda}_{01}\right)\right)\mathrm{,}$

${T}_{23}^{M}={T}_{32}^{M\ast}=-i{\left(\frac{{\lambda}_{1}}{{\lambda}_{00}+{\lambda}_{01}}\right)}^{\frac{1}{2}}\sqrt{1-\mathrm{exp}\left(-2\pi \left({\lambda}_{00}+{\lambda}_{01}\right)\right)},$

${T}_{33}^{M}=1-\frac{{\lambda}_{01}}{{\lambda}_{00}+{\lambda}_{01}}\left(1-\mathrm{exp}\left(-\pi \left({\lambda}_{00}+{\lambda}_{01}\right)\right)\right)\mathrm{,}$

$\pi {\lambda}_{00}={\left(a{a}^{*}\right)}_{ord}\text{\hspace{0.05em}},\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.17em}}\pi {\lambda}_{01}={\left(b{b}^{*}\right)}_{ord},$

and the symbol ord representing the mean time order.

${\lambda}_{00}={\left|{\displaystyle \underset{-\infty}{\overset{+\infty}{\sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{\text{J}}_{n}\left(z\right)\mathrm{exp}\left(-\frac{i{v}_{1n}^{2}}{2\hslash {\beta}_{0}}\right)\right|}^{2}\mathrm{,}$ (56)

${\lambda}_{11}={\left|{\displaystyle \underset{-\infty}{\overset{+\infty}{\sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{\text{J}}_{n}\left(z\right)\mathrm{exp}\left(-\frac{i{v}_{2n}^{2}}{2\hslash {\beta}_{0}}\right)\right|}^{2}\mathrm{,}$ (57)

The field function of the system can now be written as

$\psi =\left[\begin{array}{ccc}{T}_{11}^{M}& {T}_{12}^{M}& {T}_{13}^{M}\\ {T}_{21}^{M}& {T}_{22}^{M}& {T}_{23}^{M}\\ {T}_{31}^{M}& {T}_{32}^{M}& {T}_{33}^{M}\end{array}\right]\left[\begin{array}{c}|-1\rangle \\ |0\rangle \\ |1\rangle \end{array}\right]\mathrm{,}$ (58)

From the field vector (or transition matrix) above it is seen that the probable occupation varies with the initial occupied state; therefore, if the initial occupation is $|-1\rangle $, then the survival probability is given by the quantities

${P}_{s\left(-1\mapsto -1\right)}={\left|{T}_{11}^{M}\right|}^{2}={\left(1-\frac{{\lambda}_{00}}{{\lambda}_{00}+{\lambda}_{01}}\left(1-\mathrm{exp}\left(-\pi \left({\lambda}_{00}+{\lambda}_{01}\right)\right)\right)\right)}^{2}\mathrm{,}$ (59)

while the transition probabilities of the system are

$\begin{array}{l}{P}_{-1\mapsto 0}={\left|{T}_{12}^{M}\right|}^{2}=\frac{{\lambda}_{00}^{2}}{{\left({\lambda}_{00}+{\lambda}_{01}\right)}^{2}}\left(1-\mathrm{exp}\left(-2\pi \left({\lambda}_{00}+{\lambda}_{01}\right)\right)\right)\mathrm{,}\hfill \end{array}$ (60)

$\begin{array}{l}{P}_{-1\mapsto 1}={\left|{T}_{13}^{M}\right|}^{2}=\frac{{\lambda}_{00}{\lambda}_{01}}{{\left({\lambda}_{00}+{\lambda}_{01}\right)}^{2}}{\left(1-\mathrm{exp}\left(-\pi \left({\lambda}_{00}+{\lambda}_{01}\right)\right)\right)}^{2}\mathrm{,}\hfill \end{array}$ (61)

When the system is prepared in such a way that the initial occupation is $|0\rangle $, the survival probability is given by the quantities

${P}_{s\left(0\mapsto 0\right)}={\left|{T}_{22}^{M}\right|}^{2}=\mathrm{exp}\left(-2\pi \left({\lambda}_{00}+{\lambda}_{01}\right)\right),$ (62)

while the transition probabilities of the system are

${P}_{0\mapsto -1}={\left|{T}_{21}^{M}\right|}^{2}=\frac{{\lambda}_{00}}{{\lambda}_{00}+{\lambda}_{01}}\left(1-\mathrm{exp}\left(-2\pi \left({\lambda}_{00}+{\lambda}_{01}\right)\right)\right),$ (63)

${P}_{0\mapsto 1}={\left|{T}_{23}^{M}\right|}^{2}=\frac{{\lambda}_{1}}{{\lambda}_{00}+{\lambda}_{01}}\left(1-\mathrm{exp}\left(-2\pi \left({\lambda}_{00}+{\lambda}_{01}\right)\right)\right),$ (64)

The last case is the most relevant since for any central lobe, we have two relevant lobes. The multi-crossing behaviour is obtained when the relaxation time v, which depends on the external capacitor, is large. In fact, the system carries intrinsically a multiple crossing behaviour. If the relaxation time v is small, the induced energy gap ${J}_{n}\left(z\right)$ operates in a really small time interval; therefore fictitious crossings are not expected. However, if v is large, the induced energy gaps around $n=0$ mean that $\cdots ,{J}_{-2}\left(z\right)\mathrm{,}{J}_{-1}\left(z\right)\mathrm{,}{J}_{1}\left(z\right)\mathrm{,}{J}_{2}\left(z\right),\cdots $ operate separately. Therefore, the system behaves at each of these points as if they where an anti-crossing (avoided crossing); the fictitious crossing are expected and the probability changes drastically [52] [75].

The survival probability as function of dimensionless time for different values of v is illustrated in Figure 9.

Figure 9 shows the variation of the survival probability as a function of time for different values of v. It is observed that as v increases, we register a drastic change in the values of the probabilities. This phenomenon can easily be understood when the notion of fictitious crossings is introduced. In fact, for small relaxation time, the system does not have time to adapt to change before it faces the next crossing. Therefore each of the induced gap adds in complimentary manner and it results in a deep fall in survival probability at the effective crossing point (see Figure 9(a), dotted black curve).

In the case of large relaxation time, the system after facing the first crossing adapts to the change imposed, and saturates before it comes to the next crossing;

Figure 9. Survival probability (solid blue line) and transition probability (dotted-black line) versus dimensionless time ( $\tau $ ) for a two-state configuration (constructive interference) for parameters: (a) $v=0$, $v=1.25$ $\frac{2\alpha \beta {e}_{0}}{\hslash \Omega}=\frac{180}{\pi}$ ; (b) $v=3$, $v=5$, $\frac{2\alpha \beta {e}_{0}}{\hslash \Omega}=\frac{180}{\pi}$.

each of the induced gaps adds separately and because of the phase accumulated after each of these crossing, the transition probabilities can be constructed or destroyed and the overall is the probability limit. From the energy time diagram, the system presents only three crossings. For small relaxation those three crossings occur at around $t=0$ ; this is justified by the single crossing observed on the diagram Figure 9(a) (dotted black curve). For large relaxation, additionally to the two crossings observed on energy curves, we have many other ones. Since they are not part of the energy diagram, they are known as fictitious crossings. Their appearance is conditioned by the relaxation time v and the effective gap $\Delta {J}_{n}\left(z\right)$. Observing the solid blue line of Figure 9, we can count the number of crossings in the system; this is due to the fact that the gap $\Delta {J}_{n}\left(z\right)$ decays as n increases. Therefore, when n is large, the gap vanishes and the probability saturates.

4.3. Multi-State System

Here, we consider the same situation as previously (section 3.B) but, we take into consideration all relevant states. So $\beta {a}_{1}\frac{2\mu}{{\hslash}^{2}}\ll 1$ and the Hamiltonian of the system is labeled as in Equation (16)

$H={\displaystyle \underset{m}{\sum}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{\u03f5}_{m}|m\rangle \langle m|+\Delta |m+1\rangle \langle m|+\Delta |m-1\rangle \langle m|$.

The field vector of the system is given by:

$\psi =\left(\begin{array}{c}{C}_{-l}\\ {C}_{-l+1}\\ \vdots \\ {C}_{l}\end{array}\right)\mathrm{,}$ (65)

with $l\to +\infty $ ; we then have an infinite number of relevant states. The evaluation of the dynamic matrix in this case is via:

$i\hslash {\psi}_{\mathrm{,0}}=H{\psi}_{0}\mathrm{.}$ (66)

The transition amplitude ${C}_{m}$ can be generalized as a single differential equation

$i\hslash {\stackrel{\dot{}}{C}}_{m}={\u03f5}_{m}{C}_{m}+\Delta \left({C}_{m+1}+{C}_{m-1}\right)\mathrm{.}$ (67)

with $m=-l,\cdots ,l$, $l\to +\infty $. Because of the infinite number of values taken by m, we can rewrite that equation in terms of

$i\hslash {\psi}_{\mathrm{,0}}={\u03f5}_{m}\psi +\Delta \left({J}_{+}+{J}_{-}\right)\psi \mathrm{.}$ (68)

where

${\u03f5}_{m}\sim \frac{m}{2}\u03f5\left(t\right)+\frac{1}{2}{\u03f5}_{0}\left(t\right)\mathrm{,}$

is the elements of the energy matrix $\epsilon $ and so ${\epsilon}_{m\mathrm{,}n}={\u03f5}_{m}{\delta}_{m\mathrm{,}n}$ and

$\epsilon \left(t\right)=\left({\beta}_{0}t-z\Omega \mathrm{cos}\left(\Omega t\right)\right)\mathrm{,}$

${\epsilon}_{0}\left(t\right)={\left({\beta}_{0}t-z\Omega \mathrm{cos}\left(\Omega t\right)\right)}^{2}\mathrm{,}$

The term $\Delta \left({J}_{+}+{J}_{-}\right)$ represents the tunneling matrix. Here,

${J}_{+}=\left(\begin{array}{ccccccc}0& 1& 0& 0& 0& 0& \cdots \\ 0& 0& 1& 0& 0& 0& \cdots \\ 0& 0& 0& 1& 0& 0& \cdots \\ 0& 0& 0& 0& 1& 0& \cdots \\ 0& 0& 0& 0& 0& 1& \cdots \\ 0& 0& 0& 0& 0& 0& \cdots \\ \vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \ddots \end{array}\right)\mathrm{,}$ (69)

is an upper triangular matrix while ${J}_{-}={\left({J}_{+}\right)}^{t}$ is a lower triangular matrix with constant elements 1. Moreover, ${J}_{+}$ and ${J}_{-}$ are shifting operators acting on the k space. Taking that

$\psi =\mathrm{exp}\left(i{\displaystyle \int \epsilon \left(t\right)\text{d}t}\right)M$

Then Equation (69) is reduced as

$\stackrel{\dot{}}{M}=\frac{\text{d}}{\text{d}t}\left(\alpha \left(t\right){J}_{+}-\beta \left(t\right){J}_{-}\right)M\mathrm{.}$ (70)

with

$\Xi \left(t\right)={\displaystyle {\int}_{-\infty}^{t}}\frac{\Delta}{i\hslash}\mathrm{exp}\left(\frac{i}{\hslash}{\displaystyle {\int}_{-\infty}^{{t}^{\prime}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{\u03f5}_{m}\left(\tau \right)-{\u03f5}_{m-1}\left({\tau}^{\prime}\right)\text{d}\tau \right)\text{d}{t}^{\prime}\mathrm{.}$

and

$\Upsilon \left(t\right)={\displaystyle {\int}_{-\infty}^{t}}\frac{\Delta}{i\hslash}\mathrm{exp}\left(\frac{i}{\hslash}{\displaystyle {\int}_{-\infty}^{{t}^{\prime}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{\u03f5}_{m}\left(\tau \right)-{\u03f5}_{m+1}\left({\tau}^{\prime}\right)\text{d}\tau \right)\text{d}{t}^{\prime}\mathrm{.}$

So the dynamics matrix is expressed by

$\xi =\Xi \left(t\right){J}_{+}-\Upsilon \left(t\right){J}_{-}\mathrm{.}$ (71)

and the solution is obtained as:

$M=\stackrel{^}{T}\mathrm{exp}\left(\xi \right){M}_{-\infty}\mathrm{.}$ (72)

or

$M=\stackrel{^}{T}\mathrm{exp}\left(\Xi \left(t\right){J}_{+}-\Upsilon \left(t\right){J}_{-}\right){M}_{-\infty}\mathrm{.}$ (73)

from the above solution, the transition matrix is given by:

${T}^{M}=\stackrel{^}{T}\mathrm{exp}\left(\Xi \left(t\right){J}_{+}-\Upsilon \left(t\right){J}_{-}\right){M}_{-\infty}\mathrm{.}$ (74)

The solution is therefore being rewritten in terms of the state vector $\psi $

$\psi =\mathrm{exp}\left(\frac{i}{\hslash}{\displaystyle \int}\epsilon \left(t\right)\text{d}t\right)\cdot \stackrel{^}{T}\mathrm{exp}\left(\alpha \left(t\right){J}_{+}-\beta \left(t\right){J}_{-}\right){M}_{-\infty}\mathrm{.}$ (75)

Taking into consideration that the operators ${J}_{+}$ and ${J}_{-}$ and $\stackrel{^}{m}$ obey to the following commutation relation:

$\left[{J}_{\pm}\mathrm{,}\stackrel{^}{m}\right]=\pm {J}_{\pm}\mathrm{.}$ (76)

we use Dyson ordering or Magnus expansion to evaluate the transition matrix [50] [73] [81]; taking ${M}_{-\infty}={\displaystyle \underset{k}{\sum}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{C}_{m}{\delta}_{km}$, then the elements of the field vector are:

${C}_{m}={\displaystyle \underset{k}{\sum}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{C}_{k}\mathrm{exp}\left(\frac{ik}{\hslash}{\displaystyle {\int}_{-\infty}^{t}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{\u03f5}_{0}\left(\tau \right)+{\u03f5}_{0}\left(\tau \right)\text{d}\tau \right){\left(\frac{\Xi \left(t\right)}{\Upsilon \left(t\right)}\right)}^{\frac{m-k}{2}}{J}_{m-k}\left[2\sqrt{\alpha \left(t\right)\beta \left(t\right)}\right]\mathrm{.}$ (77)

Taking that ${C}_{m}\left(-\infty \right)={\delta}_{m0}$ then

${C}_{m}=\mathrm{exp}\left(\frac{i}{\hslash}{\displaystyle {\int}_{-\infty}^{t}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\u03f5\left(\tau \right)+{\u03f5}_{0}\left(\tau \right)\text{d}\tau \right){\left(\frac{\Xi \left(t\right)}{\Upsilon \left(t\right)}\right)}^{\frac{l}{2}}{J}_{l}\left[2\sqrt{\alpha \left(t\right)\beta \left(t\right)}\right]\mathrm{.}$ (78)

this gives us the generalized transition amplitude for an arbitrary state $\text{\hspace{0.05em}}|m\rangle $ where the probabilities are found as the square of modulus.

5. Conclusions

In this paper, we have studied a modified Josephson junction based circuit and the resulting Landau-Zener scenario. We show that the Josephson junction based circuit can be modeled as a multi-state system with crossing energy levels. The charging energy, the energy stored in the junction as well as the inductance is controllable parameters that condition the shape of the potential energy and therefore the state configuration of the system. The initial condition of the inductive branch helps in designing constructive and destructive interferences. Through the relaxation parameter, constructive interference is obtained when $\mu $ is an integer while destructive interference is obtained when $\mu $ is not an integer.

Therefore, playing on the initial condition of the coil, we can set the type of interference while, the value of the inductance controlled the energy gap as well as the frequency of probability oscillation. Moreover, the nonlinearity of energy levels modifies the transition probability and the derived adiabatic parameter helps to redefine The Landau-Zener probability. Therefore, the total population transfer as well as the complete blocking of the system is obtained in a permissible range of parameters only by changing the values of the inductance of the coil. Thus, the system models a controllable level-crossing system where the additional branches (inductive and capacitive) help to design the number of states, the type of interferometry as well as the control of state’s occupation.

Data Availability Statement

The data that supports the findings of this study are available within this article.

Acknowledgements

HAC thanks ICTP-SAIFR and FAPESP grant 2016/01343-7 for partial support. P.L. acknowledges support by the FAPESP Grant No. 2014/13272-1. HAC and PL acknowledge hospitality from the Abdus Salam ICTP.

Sincere thanks to the members of JAMP for their professional performance, and special thanks to managing editor *Hellen** XU* for a rare attitude of high quality.

References

[1] Josephson, B.D. (1962) Possible New Effects in Superconductive Tunnelling. Physics Letters, 1, 251-253.

https://doi.org/10.1016/0031-9163(62)91369-0

[2] Asai, H., Ota, Y., Kawabata, S., Machida, M. and Nori, F. (2014) Theory of Macroscopic Quantum Tunneling with Josephson-Leggett Collective Excitations in Multiband Superconducting Josephson Junctions. Physical Review B, 89, 224507.

https://doi.org/10.1103/PhysRevB.89.224507

[3] Barone, A. and Paternò, G. (1982) Physics and Applications of the Josephson Effect. John Wiley & Sons, Inc., Hoboken.

https://doi.org/10.1002/352760278X

[4] Anderson, P.W. and Rowell, J.M. (1963) Probable Observation of the Josephson Superconducting Tunneling Effect. Physical Review Letters, 10, 230-232.

https://doi.org/10.1103/PhysRevLett.10.230

[5] Shapiro, S. (1963) Josephson Currents in Superconducting Tunneling: The Effect of Microwaves and Other Observations. Physical Review Letters, 11, 80-82.

https://doi.org/10.1103/PhysRevLett.11.80

[6] Jäck, B., Eltschka, M., Assig, M., Hardock, A., Etzkorn, M., Ast, C.R. and Kern, K. (2015) A Nanoscale Gigahertz Source Realized with Josephson Scanning Tunneling Microscopy. Applied Physics Letters, 106, 013109.

https://doi.org/10.1063/1.4905322

[7] van den Brink, A.M., Odintsov, A.A., Bobbert, P.A. and Schön, G. (1991) Coherent Cooper Pair Tunneling in Systems of Josephson Junctions: Effects of Quasiparticle Tunneling and of the Electromagnetic Environment. Zeitschrift für Physik B Condensed Matter, 85, 459-467.

https://doi.org/10.1007/BF01307644

[8] Vávra, O., Soni, R., Petraru, A., Himmel, N., Vávra, I., Fabian, J., Kohlstedt, H. and Strunk, C. (2017) Coexistence of Tunneling Magnetoresistance and Josephson Effects in SFIFS Junctions. AIP Advances, 7, 025008.

https://doi.org/10.1063/1.4976822

[9] Paternò, G., Cucolo, A.M. and Modestino, G. (1985) Resonant Modes in NB Baselayer Interferometers with Two Josephson Junctions. Journal of Applied Physics, 57, 1680-1685.

https://doi.org/10.1063/1.334437

[10] Gurney, R.W. and Condon, E.U. (1929) Quantum Mechanics and Radioactive Disintegration. Physical Review, 33, 127-140.

https://doi.org/10.1103/PhysRev.33.127

[11] Taylor, J. (2004) Modern Physics for Scientists and Engineers. Pearson Prentice Hall, Upper Saddle River.

[12] Rommel, J.B., Liu, Y., Werner, H.-J. and Kästner, J. (2012) Role of Tunneling in the Enzyme Glutamate Mutase. The Journal of Physical Chemistry B, 116, 13682-13689.

https://doi.org/10.1021/jp308526t

[13] Leggett, A.J. (1984) Quantum Tunneling in the Presence of an Arbitrary Linear Dissipation Mechanism. Physical Review B, 30, 1208-1218.

https://doi.org/10.1103/PhysRevB.30.1208

[14] Simánek, E. (1994) Inhomogeneous Superconductors. Granular and Quantum Effects: by Eugen Simanek. Oxford University Press, New York, 1994, 383 Pages, Hardback ISBN 0-19-507828-4, Price $75. Materials Research Bulletin, 32, 1137-1138.

https://doi.org/10.1016/S0025-5408(97)00068-8

[15] Takagi, S. (2002) Macroscopic Quantum Tunneling. Cambridge University Press, Cambridge.

[16] Weiss, U. (2012) Quantum Dissipative Systems. 4th Edition, World Scientific.

https://doi.org/10.1142/8334

[17] Tafuri, F. (2019) Fundamentals and Frontiers of the Josephson Effect. Springer-Verlag GmbH.

https://www.ebook.de/de/product/37902246/fundamentals_and_frontiers_of_the_josephson_effect

.html

https://doi.org/10.1007/978-3-030-20726-7

[18] Kato, T. and Imada, M. (1996) Macroscopic Quantum Tunneling of a Fluxon in a Long Josephson Junction. Journal of the Physical Society of Japan, 65, 2963-2975.

https://doi.org/10.1143/JPSJ.65.2963

[19] Kawabata, S., Kashiwaya, S., Asano, Y. and Tanaka, Y. (2005) Effect of Zero-Energy Bound States on Macroscopic Quantum Tunneling in High-Tc Superconductor Junctions. Physical Review B, 72, 052506.

https://doi.org/10.1103/PhysRevB.72.052506

[20] Machida, M. and Koyama, T. (2007) Theory for Collective Macroscopic Tunneling in High-Tc Intrinsic Josephson Junctions. Physica C: Superconductivity and Its Applications, 463-465, 84-88.

https://doi.org/10.1016/j.physc.2007.05.008

[21] Savel’ev, S., Rakhmanov, A.L. and Nori, F. (2007) Quantum Terahertz Electrodynamics and Macroscopic Quantum Tunneling in Layered Superconductors. Physical Review Letters, 98, 077002.

https://doi.org/10.1103/PhysRevLett.98.077002

[22] Sboychakov, A.O., Savelev, S., Rakhmanov, A.L. and Nori, F. (2007) Why Macroscopic Quantum Tunnelling in Josephson Junctions Differs from Tunnelling of a Quantum Particle. Europhysics Letters (EPL), 80, 17009.

https://doi.org/10.1209/0295-5075/80/17009

[23] Sboychakov, A.O., Savel’ev, S. and Nori, F. (2008) Quantum Electrodynamics and Photon-Assisted Tunneling in Long Josephson Junctions. Physical Review B, 78, 134518.

https://doi.org/10.1103/PhysRevB.78.134518

[24] Savel’ev, S., Sboychakov, A.O., Rakhmanov, A.L. and Nori, F. (2008) Nonlocal Macroscopic Quantum Tunneling and Quantum Terahertz Electrodynamics in Layered Superconductors: Theory and Simulations. Physical Review B, 77, 014509.

https://doi.org/10.1103/PhysRevB.77.014509

[25] Kawabata, S., Bauch, T. and Kato, T. (2009) Theory of Two-Dimensional Macroscopic Quantum Tunneling YBa2Cu3O7-δ Josephson Junctions Coupled to an LC Circuit. Physical Review B, 80, 174513.

https://doi.org/10.1103/PhysRevB.80.174513

[26] Ota, Y., Machida, M. and Koyama, T. (2011) Macroscopic Quantum Tunneling in Multigap Superconducting Josephson Junctions: Enhancement of Escape Rate via Quantum Fluctuations of the Josephson-Leggett Mode. Physical Review B, 83, 060503(R).

https://doi.org/10.1103/PhysRevB.83.060503

[27] Leoni, R., Carelli, P. and Foglietti, V. (1988) Stray Capacitance Effect in Superconducting Quantum Interferometers. Journal of Applied Physics, 64, 2527-2532.

https://doi.org/10.1063/1.341636

[28] Zalba, M. F. G., Haigh, J., Yang, T.-Y. and Ibberson, L. (2019) Quantum Computing Using Silicon Transistors. Hitachi Review, 68, 548-549.

[29] Vrijen, R., Yablonovitch, E., Wang, K., Jiang, H.W., Balandin, A., Roychowdhury, V., Mor, T. and DiVincenzo, D. (2000) Electron-Spin-Resonance Transistors for Quantum Computing in Silicon-Germanium Heterostructures. Physical Review A, 62, 012306.

https://doi.org/10.1103/PhysRevA.62.012306

[30] Devoret, M.H. and Schoelkopf, R.J. (2013) Superconducting Circuits for Quantum Information: An Outlook. Science, 339, 1169-1174.

https://doi.org/10.1126/science.1231930

[31] Wendin, G. and Shumeiko, V. (2007) Quantum Bits with Josephson Junctions (Review Article). Low Temperature Physics, 33, 724-744.

https://doi.org/10.1063/1.2780165

[32] Makhlin, Y., Schön, G. and Shnirman, A. (2000) Quantum-State Engineering with Josephson-Junction Devices. Reviews of Modern Physics, 73, 357-400.

https://doi.org/10.1103/RevModPhys.73.357

[33] Mooij, J., Orlando, T., Levitov, L., Tian, L., Van der Wal, C.H. and Lloyd, S. (1999) Josephson Persistent-Current Qubit. Science, 285, 1036-1039.

https://doi.org/10.1126/science.285.5430.1036

[34] Nakamura, Y., Pashkin, Y.A. and Tsai, J.S. (1999) Coherent Control of Macroscopic Quantum States in a Single-Cooper-Pair Box. Nature, 398, 786-788.

https://doi.org/10.1038/19718

[35] Sillanpää, M., Lehtinen, T., Paila, A., Makhlin, Y. and Hakonen, P. (2006) Continuous-Time Monitoring of Landau-Zener Interference in a Cooper-Pair Box. Physical Review Letters, 96, 187002.

https://doi.org/10.1103/PhysRevLett.96.187002

[36] Bladh, K., Duty, T., Gunnarsson, D. and Delsing, P. (2005) The Single Cooper-Pair Box as a Charge Qubit. New Journal of Physics, 7, 180-180.

https://doi.org/10.1088/1367-2630/7/1/180

[37] Steffen, M., Martinis, J.M. and Chuang, I.L. (2003) Accurate Control of Josephson Phase Qubits. Physical Review B, 68, 224518.

https://doi.org/10.1103/PhysRevB.68.224518

[38] Landau, L.D. (1932) Zur theorie der energieubertragung ii. Zeitschrift fur Physik (Sowjetunion), 2, 46.

[39] Zener, C. (1932) Non-Adiabatic Crossing of Energy Levels, Proceedings of the Royal Society of London. Series A, Containing Papers of a Mathematical and Physical Character 137 (833) 696-702.

https://doi.org/10.1098/rspa.1932.0165

[40] Sibille, A., Palmier, J.F. and Laruelle, F. (1998) Zener Interminiband Resonant Breakdown in Superlattices. Physical Review Letters, 80, 4506-4509.

https://doi.org/10.1103/PhysRevLett.80.4506

[41] Mullen, K., Ben-Jacob, E. and Schuss, Z. (1988) Combined Effect of Zener and Quasiparticle Transitions on the Dynamics of Mesoscopic Josephson Junctions. Physical Review Letters, 60, 1097-1100.

https://doi.org/10.1103/PhysRevLett.60.1097

[42] Wu, B. and Niu, Q. (2000) Nonlinear Landau-Zener Tunneling. Physical Review A, 61, 023402.

https://doi.org/10.1103/PhysRevA.61.023402

[43] Liu, J., Fu, L., Ou, B.-Y., Chen, S.-G., Choi, D.-I., Wu, B. and Niu, Q. (2002) Theory of Nonlinear Landau-Zener Tunneling. Physical Review A, 66, 023404.

https://doi.org/10.1103/PhysRevA.66.023404

[44] Chen, Y.-A., Huber, S. D., Trotzky, S., Bloch, I. and Altman, E. (2010) Many-Body Landau-Zener Dynamics in Coupled One-Dimensional Bose Liquids. Nature Physics, 7, 61-67.

https://doi.org/10.1038/nphys1801

[45] Suqing, D., Fu, L.-B., Liu, J. and Zhao, X.-G. (2005) Effects of Periodic Modulation on the Landau-Zener Transition. Physics Letters A, 346, 315-320.

https://doi.org/10.1016/j.physleta.2005.07.086

[46] Wang, G.-F., Fu, L.-B. and Liu, J. (2006) Periodic Modulation Effect on Self-Trapping of Two Weakly Coupled Bose-Einstein Condensates. Physical Review A, 73, 013619.

https://doi.org/10.1103/PhysRevA.73.013619

[47] Raghavan, S., Smerzi, A., Fantoni, S. and Shenoy, S.R. (1999) Coherent Oscillations between Two Weakly Coupled Bose-Einstein Condensates: Josephson Effects, π Oscillations, and Macroscopic Quantum Self-Trapping. Physical Review A, 59, 620-633.

https://doi.org/10.1103/PhysRevA.59.620

[48] Zobay, O. and Garraway, B.M. (2000) Time-Dependent Tunneling of Bose-Einstein Condensates. Physical Review A, 61, 033603.

https://doi.org/10.1103/PhysRevA.61.033603

[49] Wang, G.-F., Ye, D.-F., Fu, L.-B., Chen, X.-Z. and Liu, J. (2006) Landau-Zener Tunneling in a Nonlinear Three-Level System. Physical Review A, 74, 033414.

https://doi.org/10.1103/PhysRevA.74.033414

[50] Jipdi, M., Tchoffo, M. and Fai, L. (2018) Adiabatically Modeling Quantum Gates with Two-Site Heisenberg Spins Chain: Noise vs Interferometry. Physica E: Low-Dimensional Systems and Nanostructures, 96, 36-45.

https://doi.org/10.1016/j.physe.2017.09.023

[51] Shevchenko, S., Ashhab, S. and Nori, F. (2010) Landau-Zener-Stückelberg Interferometry. Physics Reports, 492, 1-30.

https://doi.org/10.1016/j.physrep.2010.03.002

[52] Jipdi, M., Fai, L. and Tchoffo, M. (2016) Quantum System Driven by Incoherent a.c Fields: Multi-Crossing Landau Zener Dynamics. Physics Letters A, 380, 3665-3671.

https://doi.org/10.1016/j.physleta.2016.08.044

[53] Zhang, J.-N., Sun, C.-P., Yi, S. and Nori, F. (2011) Spatial Landau-Zener-Stückelberg Interference in Spinor Bose-Einstein Condensates. Physical Review A, 83, 033614.

https://doi.org/10.1103/PhysRevA.83.033614

[54] Stückelberg, E. (1932) Theory of Inelastic Collisions between Atoms (Theory of Inelastic Collisions between Atoms, Using Two Simultaneous Differential Equations. Helvetica Physica Acta (Basel), 5, 369-422.

[55] Majorana, E. (1932) Atomi Orientati in Campo Magnetico Variabile. Il Nuovo Cimento, 9, 43-50.

https://doi.org/10.1007/BF02960953

[56] Lim, L.-K., Fuchs, J.-N. and Montambaux, G. (2015) Geometric Phase in Stückelberg Interferometry. Physical Review A, 91, 042119.

https://doi.org/10.1103/PhysRevA.91.042119

[57] Tan, X., Zhang, D.-W., Zhang, Z., Yu, Y., Han, S. and Zhu, S.-L. (2014) Demonstration of Geometric Landau-Zener Interferometry in a Superconducting Qubit. Physical Review Letters, 112, 027001.

https://doi.org/10.1103/PhysRevLett.112.027001

[58] Fuchs, G.D., Dobrovitski, V.V., Toyli, D.M., Heremans, F.J. and Awschalom, D.D. (2009) Gigahertz Dynamics of a Strongly Driven Single Quantum Spin. Science, 326, 1520-1522.

https://doi.org/10.1126/science.1181193

[59] Ribeiro, H. and Burkard, G. (2009) Nuclear State Preparation via Landau-Zener-Stückelberg Transitions in Double Quantum Dots. Physical Review Letters, 102, 216802.

https://doi.org/10.1103/PhysRevLett.102.216802

[60] Petta, J.R., Lu, H. and Gossard, A.C. (2010) A Coherent Beam Splitter for Electronic Spin States. Science, 327, 669-672.

https://doi.org/10.1126/science.1183628

[61] Wang, L., Tu, T., Gong, B., Zhou, C. and Guo, G.-C. (2016) Experimental Realization of Non-Adiabatic Universal Quantum Gates Using Geometric Landau-Zener-Stückelberg Interferometry. Scientific Reports, 6, Article No. 19048.

https://doi.org/10.1038/srep19048

[62] Louodop, P., Tchitnga, R., Fagundes, F.F., Kountchou, M., Tamba, V.K., Pando L., C.L. and Cerdeira, H.A. (2019) Extreme Multistability in a Josephson-Junction-Based Circuit. Physical Review E, 99, 042208.

https://doi.org/10.1103/PhysRevE.99.042208

[63] Sun, H., Scott, S.K. and Showalter, K. (1999) Uncertain Destination Dynamics. Physical Review E, 60, 3876-3880.

https://doi.org/10.1103/PhysRevE.60.3876

[64] Hongray, T., Balakrishnan, J. and Dana, S.K. (2015) Bursting Behaviour in Coupled Josephson Junctions. Chaos, 25, 123104.

https://doi.org/10.1063/1.4936675

[65] Mishra, A., Saha, S., Hens, C., Roy, P.K., Bose, M., Louodop, P., Cerdeira, H.A. and Dana, S.K. (2017) Coherent Libration to Coherent Rotational Dynamics via Chimeralike States and Clustering in a Josephson Junction Array. Physical Review E, 95, 010201(R).

https://doi.org/10.1103/PhysRevE.95.010201

[66] Domnguez, D. and Cerdeira, H.A. (1995) Spatiotemporal Chaos in RF-Driven Josephson Junction Series Arrays. Physical Review B, 52, 513-526.

https://doi.org/10.1103/PhysRevB.52.513

[67] Domnguez, D. and Cerdeira, H.A. (1993) Order and Turbulence in RF-Driven Josephson Junction Series Arrays. Physical Review Letters, 71, 3359-3362.

https://doi.org/10.1103/PhysRevLett.71.3359

[68] Domnguez, D. and Cerdeira, H.A. (1995) Turbulence, Chaos and Noise in Globally Coupled Josephson Junction Arrays. Physics Letters A, 200, 43-50.

https://doi.org/10.1016/0375-9601(95)00172-Y

[69] Caldeira, A.O. (2009) An Introduction to Macroscopic Quantum Phenomena and Quantum Dissipation. Cambridge University Press, Cambridge.

https://doi.org/10.1017/CBO9781139035439

[70] Hochstadt, H. (1961) Special Functions of Mathematical Physics Hardcover. Holt, Rinehart and Winston, New York.

[71] Ayub, M., Naseer, K., Ali, M. and Saif, F. (2009) Atom Optics Quantum Pendulum. Journal of Russian Laser Research, 30, 205-223.

https://doi.org/10.1007/s10946-009-9078-x

[72] Doncheski, M.A. and Robinett, R.W. (2003) Wave Packet Revivals and the Energy Eigenvalue Spectrum of the Quantum Pendulum. Annals of Physics, 308, 578-598.

https://doi.org/10.1016/S0003-4916(03)00171-4

[73] Dattoli, G., Gallardo, J.C. and Torre, A. (1988) An Algebraic View to the Operatorial Ordering and Its Applications to Optics. La Rivista Del Nuovo Cimento (1978-1999), 11, 1-79.

https://doi.org/10.1007/BF02724503

[74] Loday, J.-L., Frabetti, A., Chapoton, F. and Goichot, F. (2001) Dialgebras and Related Operads. Part of Lecture Notes in Mathematics (Volume 1763), Springer-Verlag, Berlin.

https://doi.org/10.1007/b80864

[75] Fai, L.C., Tchoffo, M. and Jipdi, M.N. (2015) Multi-Particle and Multi-State Landau-Zener Model: Dynamic Matrix Approach. The European Physical Journal Plus, 130, Article No. 71.

https://doi.org/10.1140/epjp/i2015-15071-y

[76] Ateuafack, M., Diffo, J. and Fai, L. (2016) Dynamics of a Landau-Zener Transitions in a Two-Level System Driven by a Dissipative Environment. Physica B: Condensed Matter, 483, 37-43.

https://doi.org/10.1016/j.physb.2015.12.012

[77] Dokas, I. (2009) Zinbiel Algebras and Commutative Algebras with Divided Powers. Glasgow Mathematical Journal, 52, 303-313.

https://doi.org/10.1017/S0017089509990358

[78] Novelli, J.-C. and Thibon, J.-Y. (2008) Noncommutative Symmetric Functions and Lagrange Inversion. Advances in Applied Mathematics, 40, 8-35.

https://doi.org/10.1016/j.aam.2007.05.005

[79] Wick, G.C. (1950) The Evaluation of the Collision Matrix. Physical Review, 80, 268-272.

https://doi.org/10.1103/PhysRev.80.268

[80] Bauer, M., Chetrite, R., Ebrahimi-Fard, K. and Patras, F. (2012) Time-Ordering and a Generalized Magnus Expansion. Letters in Mathematical Physics, 103, 331-350.

https://doi.org/10.1007/s11005-012-0596-z

[81] Dattoli, G. and Torre, A. (1987) SU(2) and SU(1,1) Time-Ordering Theorems and Bloch-Type Equations. Journal of Mathematical Physics, 28, 618-621.

https://doi.org/10.1063/1.527648