# Hamiltonian Dynamics and Adiabatic Invariants for Time-Dependent Superconducting Qubit-Oscillators and Resonators in Quantum Computing Systems

Authors: Jeong Ryeol Choi
Journal: Advances in Mathematical Physics
 (2015)

---

## Abstract

An adiabatic invariant, which is a conserved quantity, is useful for studying quantum and classical properties of dynamical systems. Adiabatic invariants for time-dependent superconducting qubit-oscillator systems and resonators are investigated using the Liouville-von Neumann equation. At first, we derive an invariant for a simple superconducting qubit-oscillator through the introduction of its reduced Hamiltonian. Afterwards, an adiabatic invariant for a nanomechanical resonator linearly interfaced with a superconducting circuit, via a coupling with a time-dependent strength, is evaluated using the technique of unitary transformation. The accuracy of conservation for such invariant quantities is represented in detail. Based on the results of our developments in this paper, perturbation theory is applicable to the research of quantum characteristics of more complicated qubit systems that are described by a time-dependent Hamiltonian involving nonlinear terms.


---

## Body

## 1. Introduction

Interactions of light wave with matter are ubiquitous and have many applications in modern science and technology. An interesting research topic in this line is cavity quantum electrodynamics (cavity QED) which is relevant to the interaction of an atom with light waves inside a cavity. Meanwhile, it has been shown that the role of the cavity can be played by a superconducting device or others such as a nanomechanical oscillator or a transmission line [1–3]. Superconducting devices exhibit quantum interference and the electrons in them have the ability to act like qubits that are necessary in quantum computing systems. We mainly concentrate on these devices in this research. It is demonstrated from corresponding experiments that entanglement and the Rabi oscillations appear in a strongly coupled atom-cavity system as a result of quantum coherence between light and matter [3]. These features of spatial correlations enable us to do specific tasks regarding quantum information techniques such as quantum computing [4].Recently, the study of superconducting qubit-oscillators in quantum computing systems with time-dependent parameters attracted considerable interest in the community of physics and information science [5–9]. In general, a high ability of computation with a tremendous speed can be achieved for particular mathematical problems from quantum computation on the basis of new paradigms known as superposition and entanglement. This is the most distinguishing feature of quantum computers, which cannot be mimicked by classical ones. The realization of producing universal quantum gates that fulfill quantum computations may require the use of geometric phase shifts of qubit systems along with the changes of suitable loops in the Hamiltonian under efficient quantum control. This can be performed through the use of superconducting nanocircuits [10].Meanwhile, obtaining conserved quantities, that is, adiabatic invariants, is of importance in the research of the dynamics of quantum computing systems, because not only are they helpful for deepening the understanding of the system but also they can be used for studying the quantum and classical properties of the system itself. If the time derivative of a quantity results in zero, this quantity is a kind of exact invariant. On the other hand, an adiabatic invariant is a quantity in which its time derivative gives nearly zero under a slow variation of the Hamiltonian with time. The first proposal of adiabatic invariants was sprung by a pioneering work of Burgers [11], which is for a particular class of Hamiltonian systems. Soon after, many other researches concerning adiabatic invariants were carried out via the introduction of the Hamiltonian of a given system [12–18]. Adiabatic invariants for a linear harmonic oscillator have been studied by Dykhne in connection with quantum transitions [15]. Quantum algorithms for adiabatic quantum computation, which are robust to decoherence and not affected by certain noises, were developed [19]. A hybrid quantum simulation model composed of adiabatically controlled simulation registers connected to a qubit was proposed as well [20].We will find adiabatic invariants for the Hamiltonians that describe the mechanics of qubit-oscillators and resonators in quantum computing systems that are characterized by time-varying parameters. Concerning this, the identification of adiabatic invariants for a time-dependent Hamiltonian system with complicated time-variable parameters may be a significant challenge. We will first address the invariant of a simple superconducting qubit-oscillator. Afterwards, an adiabatic invariant for a nanomechanical resonator coupled to a superconducting resonator, which is described by a complicated Hamiltonian involving a cross term, will be investigated. The Liouville-von Neumann equation will be used in order to derive conservative quantities for such dynamical systems.


## 2. A Simple Superconducting Qubit-Oscillator

Let us consider a simple superconducting qubit-oscillator of which the effective reduced Hamiltonian is given by Equation (2.2) of [21]. Typically, the two lowest energy levels of quantized states in the anharmonic LC resonator are regarded as a set of a qubit. In this case, the system is described by the Hamiltonian of the form [21] (1)
H
^
A
=
ħ
ω
a
^
†
a
^
+
1
2
+
ħ
δ
t
a
^
+
a
^
†
,where δ
t is the time-dependent coupling strength. Here, a
^ is the annihilation operator of the simple harmonic oscillator: (2)
a
^
=
L
ω
2
ħ
q
^
+
i
p
^
2
ħ
L
ω
,with p
^
=
-
i
ħ
∂
/
∂
q, where L is the inductance of the LC resonator. According to the usual convention, a
^
† is Hermitian adjoint of a
^, and this rule will also be applied to other ladder operators that will be introduced later. In many cases of previous researches, the factor 1
/
2 given in (1) is neglected for simplicity. Using (2) and its Hermitian adjoint, we can rewrite the Hamiltonian to be (3)
H
^
A
=
p
^
2
2
L
+
1
2
L
ω
2
q
^
2
+
2
L
ω
δ
t
q
^
.The classical equation of motion for this system is given by (4)
d
2
q
d
t
2
+
ω
2
q
=
-
2
ω
L
δ
t
.Although an (adiabatic) invariant is obtained under the assumption that its total change during a physical process with time-variable parameters is zero, it is represented in terms oft. A quantum (adiabatic) invariant is obtained from the following Liouville-von Neumann equation: (5)
d
I
^
A
d
t
=
∂
I
^
A
∂
t
+
I
^
A
,
H
^
A
i
ħ
=
0
.From a straightforward evaluation with this equation using (3), we have (6)
I
^
A
=
p
^
-
P
p
t
2
2
L
+
1
2
L
ω
2
q
^
-
Q
p
t
2
,where Q
p
(
t
) is a particular solution of (4) and P
p
(
t
)
=
L
d
Q
p
(
t
)
/
d
t. Most previous research findings for adiabatic invariants, such as the research based on the method of Noether’s theory [22], were performed under the supposition of the slow change of time-variables; however, the derivation of the invariant given in (6) does not require such a condition as you can see. This invariant in fact does not vary with time. It is an exact constant of motion for the given system, even for the case that the parameter δ
(
t
) undergoes fast change with time.


## 3. A Nanoresonator Coupled to a Superconducting Resonator

The theorem of adiabatic invariants originally developed for slowly time-varying systems in one-dimension may also be extended to more complicated cases. Now let us consider a nanomechanical resonator interfaced with a superconducting circuit in a linear way. In this case, the Hamiltonian is represented as (see Equation (2) of [23] and Equation (31) of [24]) (7)
H
^
B
=
ħ
ω
a
^
1
†
a
^
1
+
1
2
+
ħ
Ω
a
^
2
†
a
^
2
+
1
2
+
ħ
λ
t
a
^
1
+
a
^
1
†
a
^
2
+
a
^
2
†
,where a
^
1 and a
^
2 are annihilation operators for the nanoresonator and stripline, respectively, and λ
(
t
) is the time-dependent coupling strength between the two oscillators. We can represent the two annihilation operators to be (8)
a
^
1
=
L
n
ω
2
ħ
q
^
1
+
i
p
^
1
2
ħ
L
n
ω
,
a
^
2
=
L
s
Ω
2
ħ
q
^
2
+
i
p
^
2
2
ħ
L
s
Ω
,with p
^
j
=
-
i
ħ
∂
/
∂
q
j, where L
n and L
s are inductances for the nanomechanical resonator and the superconducting circuit, respectively. In terms of canonical variables, we can rewrite the Hamiltonian as (9)
H
^
B
=
p
^
1
2
2
L
n
+
p
^
2
2
2
L
s
+
1
2
L
n
ω
2
q
^
1
2
+
1
2
L
s
Ω
2
q
^
2
2
+
2
λ
t
L
n
L
s
ω
Ω
q
^
1
q
^
2
.If we introduce two linearly independent real classical solutions,ρ
j
,
I
(
t
) and ρ
j
,
I
I
(
t
), of the following equations [14] (10)
d
2
ρ
j
d
t
2
+
ω
~
j
2
t
ρ
j
=
0
,for each j, where j
=
1,2, and (11)
ω
~
1
2
t
=
ω
2
cos
2
⁡
φ
t
+
Ω
2
sin
2
⁡
φ
t
-
4
λ
t
ω
Ω
cos
⁡
φ
t
sin
⁡
φ
t
,
ω
~
2
2
t
=
ω
2
sin
2
⁡
φ
t
+
Ω
2
cos
2
⁡
φ
t
+
4
λ
t
ω
Ω
cos
⁡
φ
t
sin
⁡
φ
t
,
φ
t
=
1
2
tan
-
1
⁡
4
λ
t
ω
Ω
Ω
2
-
ω
2
,we can define the Wronskians in the form (12)
W
j
=
ρ
j
,
I
t
ρ
˙
j
,
I
I
t
-
ρ
˙
j
,
I
t
ρ
j
,
I
I
t
,which are time-constants. The quantities W
j defined in this way are necessary for evaluating the adiabatic invariant of the system. Hence, they will be used in the development of our theory for the adiabatic invariant later.Let us start from the Hamiltonian given in (9) in order to derive an adiabatic invariant of the system. This Hamiltonian is somewhat complicated because it is not only dependent on time, but also described in terms of a cross term q
^
1
q
^
2. Hence, to simplify the problem, we will transform the Hamiltonian to a simple form using the unitary transformation method. To do so, let us introduce the following unitary operator: (13)
U
^
B
=
exp
⁡
i
ħ
p
^
1
q
^
1
+
q
^
1
p
^
1
ln
⁡
L
n
L
s
1
/
8
exp
⁡
i
ħ
p
^
2
q
^
2
+
q
^
2
p
^
2
ln
⁡
L
s
L
n
1
/
8
exp
⁡
-
i
φ
t
ħ
p
^
1
q
^
2
-
p
^
2
q
^
1
.Then, the Hamiltonian can be transformed using the relation (14)
H
^
B
′
=
U
^
B
-
1
H
^
B
U
^
B
-
i
ħ
U
^
B
-
1
∂
U
^
B
∂
t
.Through a little evaluation using fundamental properties of quantum operators, we derive the transformed Hamiltonian in the form (15)
H
^
B
′
=
p
^
1
2
+
p
^
2
2
2
L
n
L
s
+
1
2
L
n
L
ω
~
1
2
t
q
^
1
2
+
1
2
L
n
L
ω
~
2
2
t
q
^
2
2
+
β
t
λ
˙
t
p
^
1
q
^
2
-
p
^
2
q
^
1
,where (16)
β
t
=
-
2
ω
Ω
Ω
2
-
ω
2
+
16
ω
Ω
λ
2
t
/
Ω
2
-
ω
2
.If we consider the adiabatic condition λ
˙
(
t
)
≃
0, this Hamiltonian reduces to (17)
H
^
B
′
≃
p
^
1
2
+
p
^
2
2
2
L
n
L
s
+
1
2
L
n
L
s
ω
~
1
2
t
q
^
1
2
+
1
2
L
n
L
s
ω
~
2
2
t
q
^
2
2
.This decoupled Hamiltonian is much more simple than the original one given in (9). We will use this Hamiltonian in the subsequent development of our quantum theory, instead of (15), for convenience.Similar to the previous case, the invariant operator in the transformed system should satisfy the relation(18)
d
I
^
B
′
d
t
=
∂
I
^
B
′
∂
t
+
I
^
B
′
,
H
^
B
′
i
ħ
=
0
.The direct evaluation of the invariant operator using this equation with (17) leads to (19)
I
^
B
′
=
L
n
L
s
2
∑
j
=
1,2
W
j
2
ρ
j
2
t
q
^
j
2
+
ρ
˙
j
t
q
^
j
-
ρ
j
t
L
n
L
s
p
^
j
2
,where (20)
ρ
j
t
=
ρ
j
,
I
2
t
+
ρ
j
,
I
I
2
t
1
/
2
.From the inverse transformation of (19) using I
^
B
=
U
^
B
I
^
B
′
U
^
B
-
1, we can obtain the adiabatic invariant in the original system as (21)
I
^
B
=
L
n
L
s
2
∑
j
=
1,2
W
j
2
ρ
j
2
t
Q
^
j
2
+
ρ
˙
j
t
Q
^
j
-
ρ
j
t
L
n
L
s
P
^
j
2
,where (22)
Q
^
1
=
L
n
L
s
1
/
4
cos
⁡
φ
t
q
^
1
-
L
s
L
n
1
/
4
sin
⁡
φ
t
q
^
2
,
Q
^
2
=
L
n
L
s
1
/
4
sin
⁡
φ
t
q
^
1
+
L
s
L
n
1
/
4
cos
⁡
φ
t
q
^
2
,
P
^
1
=
L
s
L
n
1
/
4
cos
⁡
φ
t
p
^
1
-
L
n
L
s
1
/
4
sin
⁡
φ
t
p
^
2
,
P
^
2
=
L
s
L
n
1
/
4
sin
⁡
φ
t
p
^
1
+
L
n
L
s
1
/
4
cos
⁡
φ
t
p
^
2
.Thus, the adiabatic invariant of the system is completely identified now. Because we have used the condition λ
˙
(
t
)
≃
0 in the derivation of the invariant quantity, (21) is valid only when λ
(
t
) varies sufficiently slowly. Hence, the accuracy of conservation for I
^
B is determined by the degree of time-variance of λ
(
t
). Such accuracy for a conserved quantity is of theoretical interest in the general physics literature [16–18].


## 4. Conclusion

Adiabatic invariants for superconducting qubit-oscillator systems and nanomechanical resonators have been investigated. An exact invariant for a simple superconducting qubit-oscillator was evaluated. An adiabatic invariant for a nanoresonator coupled to a superconducting resonator was also obtained as shown in (21). The invariant I
^
B for the latter case is more complicated than the former I
^
A. Notice that I
^
B was derived under the assumption that the time-dependent parameter λ
(
t
) varies sufficiently slowly, while the derivation of I
^
A required no such assumption. The concept of adiabatic invariants is important in both quantum and classical points of view [25, 26], because we can study the characteristics of the system by making use of them. In particular, these invariants enable us to derive Schrödinger solutions of the dynamical systems [13].For a more complicated system, the evaluation of adiabatic invariants may be not an easy task. We do not know how to derive adiabatic invariants for the case that the Hamiltonian involves higher order terms. For such example, let us consider the Hamiltonian:(23)
H
^
C
=
ħ
ω
a
^
1
†
a
^
1
+
1
2
+
ħ
Ω
a
^
2
†
a
^
2
+
1
2
+
ħ
3
/
2
g
t
a
^
2
†
a
^
2
a
^
1
+
a
^
1
†
.The last term in this equation represents a nonlinear interaction which typically appears in a mechanical resonator coupled on superconducting stripline resonator [23]. In this case, it is difficult to know the adiabatic invariant due to the nonlinear term, and hence, to our knowledge, the use of perturbation theory is the best policy for studying quantum mechanical properties of this system. For a more detailed illustration for this, let us rewrite (23) as (24)
H
^
C
=
H
^
A
′
+
g
t
L
n
ω
2
L
s
Ω
q
^
2
2
q
^
1
+
1
L
s
Ω
p
^
2
2
q
^
1
,where (25)
H
^
A
′
=
p
^
1
2
2
L
n
+
p
^
2
2
2
L
s
+
1
2
L
n
ω
2
q
^
1
2
+
1
2
L
s
Ω
2
q
^
2
2
-
ħ
g
t
L
n
ω
2
q
^
1
.By applying the method of previous development given for a simple superconducting resonator cavity, we can easily derive the adiabatic invariant of H
^
A
′. If it is possible to derive quantum solutions of the system described by H
^
A
′, we can use perturbation theory in order to derive the quantum solutions of H
^
C after regarding the last term in (24) as the perturbation term [27].The accuracy of conservation for invariant quantities is crucial for securing the validity of the associated analyses of classical and quantum mechanical characteristics of the system fulfilled on the basis of such invariants [16]. In this research, I
^
A given in (6) is an exact invariant while I
^
B given in (21) is an adiabatic invariant valid under the condition that λ
(
t
) is a slowly varying function. Both invariants are useful for investigating the classical and the quantum properties of each respective system. As a further task in the future, it may be a good subject to investigate quantum mechanical properties of superconducting qubit-oscillator systems employing the theory of invariants developed in this work.



---

*Source: 120573-2015-11-12.xml*