The transmission-line shunted plasma oscillation qubit, or transmon in short, is a variant of superconducting charge qubit that is insensitive to charge noise. It consists of two small superconducting islands tunnel coupled through two Josephson junctions (forming a dc-SQUID) and shunted by a large capacitor (see
Figure 1
). The key feature that distinguishes the transmon from the so-called Cooper-pair box (Nakamura, Pashkin and Tsai, 1999) is to reduce the sensitivity to charge noise using the capacitive shunt while maintaining a sufficient anharmonicity for selective control of a pair of levels to form a quantum bit (Koch et al., 2007).
Figure 1. A schematic of transmon. It consists of two superconducting islands tunnel coupled through two Josephson junctions and shunted by a large capacitor.
The Transmon package is a part of QuantumPlaybook, but needs to be loaded separately.
In[]:=
PacletInstall["QuantumMob/QuantumPlaybook"]
Out[]=
PacletObject
Name: QuantumMob/QuantumPlaybook
Version: 4.0.13

In[]:=
Needs["QuantumMob`QuantumPlaybook`"]

Transmon Hamiltonian

A transmon is a small Josephson junction and is described by the Hamiltonian

H
4
E
C
2


n
-q
+
E
J
1-cos

ϕ

,
where
E
C
:
2
e
2C
is the charging energy associated with junction capacitance
C
,
E
J
is the Josephson energy, and
q:
C
g
V
g
2e
is the dimensionless gate charge for gate capacitance
C
g
and voltage
V
g
(see Figure 2). Two operators

n
and

ϕ
describe the number of Cooper pairs and the phase of the superconducting order parameter, respectively. The two operators

n
and

ϕ
satisfy the commutation relation


ϕ
,

n
i
.
Due to this commutation relation, the phase and number are subjected to quantum fluctuations satisfying the uncertainty relation
Δϕ Δn≥
1
2
.
In this sense, the above transmon Hamiltonian is also called the quantum phase model.
Figure 2. An equivalent circuit for transmon. A transmon is a small Josephson junction characterized by Josephson critical current
I
J
and shunted by a large capacitance
C
. The Josephson critical current
I
J
determines the Josephson energy
E
J
I/2e
while the capacitance
C
gives the charging energy scale
E
C

2
e
2C
. The transmon device is often tuned by an external gate voltage
V
g
via a gate capacitor of capacitance
C
g
, inducing the dimensionless gate charge
n
g

C
g
V
g
2e
.
The quantum fluctuations of the phase and number are governed by the two energy scales
E
C
and
E
J
in the transmon Hamiltonian. If
E
C
≫
E
J
, then the first term dominates, and the number is well defined while the phase strongly fluctuates. In the opposite limit,
E
C
≪
E
J
, the Josephson term dominates and the phase is well defined whereas the Cooper-pair number has a large uncertainty. It is thus convenient to characterize the quantum fluctuations with the ratio of the two energy scales, more precisely, with
α:
8
E
C
E
J
.
which is called the MacCumber parameter of the transmon device. Furthermore, we measure energies in units of the geometric mean of the two energy scales, that is,
8
E
C
E
J
8
E
C
αα
E
J
.
Then, the dimensionless transmon Hamiltonian reads as

H

1
2
α
2


n
-q
+
1
α
1-cos

ϕ

.
In the ϕ-space representation, the corresponding Schrödinger equation is given by

α
2
2
-i
∂
ϕ
-q
+
1
α
(1-cosϕ)ψ(ϕ)ϵψ(ϕ)
,
∂
ϕ
≡
∂
∂ϕ
.
Phases 0 and 2π are equivalent, so physically meaningful wave function must satisfy the periodic boundary condition
ψ(ϕ+2π)ψ(ϕ)
.

Bloch States

Here, we consider a slightly different model described by the Hamiltonian

H

1
2
α
2

n
+
1
α
1-cos

ϕ

and the open boundary condition (instead of the
periodic boundary condition
). In the ϕ-space representation, the Schrödinger equation reads as
-
1
2
α
2
∂
ϕ
+
1
α
(1-cosϕ)ψ(ϕ)ϵψ(ϕ)
.
We call the above Hamiltonian the Mathieu Hamiltonian and the model consisting of the Mathieu Hamiltonian (or equivalently the Mathieu-Schrödinger equation) with the open boundary condition the Mathieu model.
Notice the absence of the dimensionless gate charge
q
from this Hamiltonian (cf. the above specification of the
transmon Hamiltonian
or the corresponding
Schrödinger equation
). Without the periodic boundary condition, the dimensionless gate charge
q
can always be "gauged away". To see this, assume a solution of the form
ψ(ϕ)
iqϕ
e

ψ
(ϕ)
in the
Schrödinger equation
with finite
q
. Then, one can see that

ψ
(ϕ)
satisfies the above
Mathieu-Schrödinger equation
without
q
.

Bloch Theorem

Note that the Hamiltonian has a discrete translational symmetry. More specifically, let

T
(θ)exp(-iθ

n
)
be the translation operator by amount θ such that
†

T
(θ)

ϕ

T
(θ)

ϕ
+θ
.
The Hamiltonian commutes with any translation by integer multiples of 2π,


H
,

T
(2πm)0
for every integer
m
.
According to the Bloch (or Floquet) theorem, the discrete translational symmetry allows us to choose eigenstates of the Hamiltonian of the form
ψ
ν
(ϕ)
iν ϕ
e
f
ν
(ϕ)
,
where
f
ν
(ϕ)
is a periodic function of period 2π and real parameter ν is called the characteristic exponent (or the quasi-momentum) of the wave function. A state with wave function of the above form is called a Bloch (or Floquet) state. Note that the Bloch wave functions are defined only up to multiples of 2π; that is, if
ψ
ν
(ϕ)
is an eigenstate of the Bloch Hamiltonian

H
, then
2πmx
e
ψ
ν
(ϕ)
is also an eigenstate of

H
with the same energy as
ψ
ν+m
(ϕ)
for any integer
m
.
In the Transmon package, the wave function
ψ
ν
(ϕ)
for MacCumber parameter α is given by
BlochFunction[ν,α,ϕ]
, and the corresponding energy by
BlochEnergy[ν,α,ϕ]
.
In[]:=
Get["QuantumMob`Transmon`"]
Consider an eigenfunction with a specific characteristic exponent:
In[]:=
α=2.;​​ν=1.25;​​sol[x_]=BlochFunction[ν,α,x*Pi];
The wave function itself is not periodic due to the non-integer value of the characteristic exponent.
In[]:=
ReImPlot[sol[x],{x,-2,2},​​FrameLabel->{"ϕ / π","ψ"}]
Out[]=
However, with the plane-wave factor removed, the remaining part is periodic.
The energy depends on the characteristic exponent. The dispersion relation forms bands.

Mathieu Functions

Transmon Eigenstates

and wave function
As an example, consider the first excited state of the transmon Hamiltonian:
Clearly, the above wave function satisfies the Schrödinger equation:
Consider another example in the opposite limit.
As the MacCumber parameter α varies, the spectrum gradually changes from the Josephson to charging limit.
The degeneracy in the Josephson limit is lifted for finite dimensionless gate charge q.

Related Links

◼
  • Y. Nakamura, Y. A. Pashkin and J. S. Tsai (1999), Nature 398, 786 (1999), “Coherent control of macroscopic quantum states in a single-Cooper-pair box.”
  • ◼
  • J. Koch et al. (2007), Physical Review A 76, 042319 (2007), “Charge-insensitive qubit design derived from the Cooper pair box.”
  • ◼
  • G. Blanch (1972), “Mathieu Functions” in Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, edited by M. Abramowitz and I. A. Stegun (John Wiley & Sons, 1972).
  • CITE THIS NOTEBOOK

    Transmon: a system for superconducting quantum computers​
    by Mahn-Soo Choi​
    Wolfram Community, STAFF PICKS, August 7, 2025
    ​https://community.wolfram.com/groups/-/m/t/3526952