PHYS598500 · Week 3 · T2
Course Duration: 120 minutes
Core Question: How can we understand fluxonium through noise operators, matrix elements, and symmetry, then generalize those principles into Hamiltonian engineering?
Difficulty: Very high
Prerequisites: W2 circuit quantization, Transmon energy scales, external flux, matrix element
Original Source Material: T2 modules 5, 6
5 Questions
Under what conditions can an external perturbation cause a quantum transition?
What quantum-state structures can suppress noise transitions while preserving control transitions?
What does adding a superinductor contribute to the Hamiltonian of a Josephson circuit?
How does fluxonium use \(E_C,E_J,E_L\) and external flux to shape its wavefunctions and realize protection conditions?
Hamiltonian engineering: how can we design the Hamiltonian of a new qubit starting from the error operator we want to suppress?
Last week, the transmon showed how choosing \(E_J/E_C\) balances charge sensitivity against anharmonicity. But reduced parameter sensitivity does not mean that every error transition disappears. This week asks a deeper question: can we directly design the quantum states, coupling operators, and Hamiltonian terms so that the environment has difficulty causing relaxation? And if the environment cannot readily drive a transition, will that transition also become harder for us to control?
Driving question for this week: By combining a large inductance with a Josephson junction, fluxonium creates a spectrum unlike the transmon's. Which transitions are protected, and which control or readout operations become more difficult? How do we avoid mistaking a single coherence record for an overall architectural advantage?
| Time | Question | Outcome |
|---|---|---|
| 0–10 min | W2 review and prerequisite diagnostic | Transmon trade-off → error-target design |
| 10–30 min | Question 1: How noise causes a transition | operator × matrix element × spectrum |
| 30–50 min | Question 2: Quantum states that are protected yet controllable | sweet spot/selection rule/suppression |
| 50–70 min | Question 3: Adding a superinductor | circuit modification → Hamiltonian term |
| 70–95 min | Question 4: Energy-scale and wavefunction engineering | \(E_C,E_J,E_L,\Phi_{\rm ext}\) → spectrum/matrix element |
| 95–118 min | Question 5: Hamiltonian engineering | Four criteria, a reverse-design method, and three advanced directions |
| 118–120 min | Exit ticket | Leave with an entry point for continued self-study |
First ask students, without consulting references, to draw on one page the flow of information from the abstract model to the real hardware for this week's topic. Do not reveal the answer immediately; instead, collect which layers students mistakenly equate. At the end of the class, have them redraw the diagram and use the differences between the two versions to assess whether they have developed genuine cross-layer understanding.
Do not begin by memorizing a transition-rate formula. Begin with the system, the environment, and their coupling. The most general weak-coupling form is
\[ \hat H_{\rm tot} =\hat H_{\rm S}(\lambda_0)+\hat H_{\rm B}+\hat H_{\rm SB}, \qquad \hat H_{\rm S}(\lambda_0)\ket n=E_n\ket n. \]
\(\hat H_{\rm S}\) determines the energy levels and eigenstates, while \(\hat H_{\rm B}\) describes the environment. What actually determines whether the environment can change the qubit state is \(\hat H_{\rm SB}\). If a circuit parameter experiences a small perturbation \(\lambda(t)=\lambda_0+\delta\lambda(t)\), expand the Hamiltonian directly:
\[ \begin{aligned} \hat H_{\rm S}(\lambda_0+\delta\lambda) &\simeq \hat H_0 +\delta\lambda(t) \left.\frac{\partial\hat H_{\rm S}}{\partial\lambda}\right|_{\lambda_0} +O(\delta\lambda^2)\\ &\equiv \hat H_0+F(t)\hat O. \end{aligned} \]
Thus \(F(t)\) is not merely an abstract “noise”; it is the actual fluctuation, while \(\hat O=\partial_\lambda\hat H_{\rm S}\) is determined by the circuit Hamiltonian. If the environment is also quantized, write \(\hat H_{\rm SB}=\sum_\alpha\hat O_\alpha\otimes\hat B_\alpha\); in a semiclassical noise model, \(\hat B_\alpha\) is replaced by a stochastic function \(F_\alpha(t)\). Both descriptions retain the same system-side coupling operator.
\[ \hat O=\sum_{m,n}O_{mn}\ket m\!\bra n, \qquad O_{mn}=\langle m|\hat O|n\rangle. \]
This immediately separates the same perturbation into two physical effects:
\[ \hat V(t)=F(t) \left[ \underbrace{\sum_n O_{nn}\ket n\!\bra n}_{\text{diagonal: changes phase and transition frequency}} + \underbrace{\sum_{m\ne n}O_{mn}\ket m\!\bra n}_{\text{off-diagonal: changes population}} \right]. \]
Define the Bohr frequency \(\omega_{mn}=(E_m-E_n)/\hbar\). In the interaction picture of \(\hat H_0\),
\[ \hat V_I(t)=F(t)\sum_{m,n}O_{mn}e^{i\omega_{mn}t}\ket m\!\bra n. \]
The interaction-picture Schrödinger equation and the first-order term of the Dyson expansion are
\[ i\hbar\frac{d}{dt}\ket{\psi_I(t)}=\hat V_I(t)\ket{\psi_I(t)}, \qquad \hat U_I(t,0)\simeq 1-\frac{i}{\hbar}\int_0^t dt'\,\hat V_I(t'). \]
If the system begins in \(\ket i\), projection onto \(\ket f\) gives the first-order \(i\to f\) amplitude:
\[ c_f^{(1)}(t) =-\frac{i}{\hbar}O_{fi} \int_0^t dt'\,F(t')e^{i\omega_{fi}t'}. \]
The transition probability is therefore
\[ P_{i\to f}(t) =\frac{|O_{fi}|^2}{\hbar^2} \left| \int_0^t dt'\,F(t')e^{i\omega_{fi}t'} \right|^2. \]
This line already gives two necessary conditions. First, symmetry must not force \(O_{fi}\) to zero. Second, \(F(t)\) must have a Fourier component near \(|\omega_{fi}|\). For a finite observation time, the resonance width is approximately \(1/t\); strict energy conservation emerges only gradually in the long-time limit.
For an intentional drive, let \(F(t)=F_d\cos\omega_dt\). When \(\omega_d\simeq|\omega_{fi}|\) and the rotating-wave approximation holds,
\[ \Omega_{fi}=\frac{F_d|O_{fi}|}{\hbar}, \qquad P_{i\to f}(t)= \frac{\Omega_{fi}^2}{\Omega_{fi}^2+\Delta^2} \sin^2\!\left(\frac{t}{2}\sqrt{\Omega_{fi}^2+\Delta^2}\right), \quad \Delta=\omega_d-|\omega_{fi}|. \]
This shows that the strength of a control transition is set by the drive-operator matrix element. For stationary random noise, define the correlation function and two-sided spectral density:
\[ C_F(\tau)=\langle F(\tau)F(0)\rangle, \qquad S_F(\omega)=\int_{-\infty}^{\infty}d\tau\, e^{i\omega\tau}C_F(\tau). \]
Ensemble-averaging the amplitude gives
\[ \begin{aligned} \big\langle P_{i\to f}(t)\big\rangle &=\frac{|O_{fi}|^2}{\hbar^2} \int_0^t dt_1\int_0^t dt_2\, e^{i\omega_{fi}(t_1-t_2)}C_F(t_2-t_1)\\ &\xrightarrow[t\gg\tau_c]{} t\,\frac{|O_{fi}|^2}{\hbar^2}S_F(\omega_{if}), \qquad \omega_{if}=\frac{E_i-E_f}{\hbar}. \end{aligned} \]
The Fermi-golden-rule transition rate is therefore
\[ \boxed{ \Gamma_{i\to f} =\frac{1}{\hbar^2} \left|\langle f|\hat O|i\rangle\right|^2 S_F(\omega_{if})}. \]
This result assumes weak coupling, a stationary bath, \(t\) longer than the bath correlation time, and negligible higher-order back-action. If a text absorbs \(2\pi\), \(\hbar\), or the density of states into the definition of spectral density, the prefactor may instead appear in forms such as \(2\pi/\hbar\). Always compare the definition of \(S_F\) before comparing the visual form of equations.
Let \(\omega_{01}=(E_1-E_0)/\hbar>0\). The same noise channel gives
\[ \Gamma_\downarrow =\frac{|O_{01}|^2}{\hbar^2}S_F(+\omega_{01}), \qquad \Gamma_\uparrow =\frac{|O_{01}|^2}{\hbar^2}S_F(-\omega_{01}), \qquad \frac{1}{T_1}=\Gamma_\downarrow+\Gamma_\uparrow. \]
For a thermal bath, detailed balance gives \(S_F(-\omega)=e^{-\hbar\omega/k_BT}S_F(+\omega)\). At low temperature, usually \(\Gamma_\uparrow\ll\Gamma_\downarrow\), but upward transitions and thermal population cannot be neglected for a low-frequency qubit or a high effective environmental temperature. For multiple uncorrelated channels, the total rate is the sum of the channel rates; correlated channels require the cross-spectral density \(S_{\alpha\beta}(\omega)\). These upward and downward rates use the unsymmetrized quantum-noise spectrum. A symmetrized spectrum loses the direction of energy flow and cannot distinguish the two on its own.
The diagonal part does not directly exchange population; it changes the energy difference between the two levels:
\[ \delta\omega_{01}(t) =\frac{F(t)}{\hbar}\bigl(O_{11}-O_{00}\bigr), \qquad \delta\phi(t)=\int_0^t dt'\,\delta\omega_{01}(t'). \]
If \(F(t)=\delta\lambda(t)\) and \(\hat O=\partial_\lambda\hat H_0\), the Hellmann–Feynman theorem directly connects the diagonal matrix element to frequency sensitivity:
\[ O_{nn}=\left\langle n\left|\frac{\partial\hat H_0}{\partial\lambda}\right|n\right\rangle =\frac{\partial E_n}{\partial\lambda}, \qquad \frac{O_{11}-O_{00}}{\hbar} =\frac{\partial\omega_{01}}{\partial\lambda}. \]
Therefore,
\[ \delta\omega_{01}(t) \simeq \frac{\partial\omega_{01}}{\partial\lambda}\delta\lambda(t) +\frac12\frac{\partial^2\omega_{01}}{\partial\lambda^2}\delta\lambda^2(t)+\cdots. \]
A sweet spot \(\partial_\lambda\omega_{01}=0\) eliminates only first-order frequency fluctuations. Second-order sensitivity, other noise parameters, and finite-frequency relaxation remain. Dephasing measured by Ramsey, echo, or dynamical-decoupling sequences is the low-frequency \(S_\lambda(\omega)\) weighted by the corresponding filter function.
Consider a Hamiltonian containing charging, Josephson, and inductive energies:
\[ \hat H_0 =4E_C(\hat n-n_g)^2 -E_J\cos\hat\varphi +\frac12E_L(\hat\varphi-\varphi_{\rm ext})^2. \]
A small parameter fluctuation does not require guessing an extra operator; differentiate directly:
\[ \delta\hat H(t) \simeq \delta n_g(t)\underbrace{\big[-8E_C(\hat n-n_g)\big]}_{\hat O_{n_g}} +\delta\varphi_{\rm ext}(t) \underbrace{\big[-E_L(\hat\varphi-\varphi_{\rm ext})\big]}_{\hat O_{\Phi}} +\delta E_J(t)\underbrace{\big[-\cos\hat\varphi\big]}_{\hat O_{E_J}}. \]
If the physical flux \(\Phi_{\rm ext}\) is used, multiply additionally by \(\partial\varphi_{\rm ext}/\partial\Phi_{\rm ext}=2\pi/\Phi_0\). Likewise, a capacitive microwave drive usually couples through \(\hat n\). The diagonal sensitivity and off-diagonal matrix element must therefore be computed separately for different operators using the same set of eigenstates.
Charge noise can arise from oxide defects, two-level systems (TLSs), or gate-voltage fluctuations and couples to the Cooper-pair number:
\[ \hat V_{\rm charge}(t) =-8E_C\,\delta n_g(t)(\hat n-n_g) \simeq-8E_C\,\delta n_g(t)\hat n, \qquad \delta n_g=\frac{C_g\delta V_g}{2e}. \]
Low-frequency charge noise often has \(1/f\) behavior. After increasing \(E_J/E_C\) in a transmon, what is exponentially suppressed is the charge dispersion—the dependence of energy on offset charge \(n_g\). This is a flattening of frequency sensitivity and must not be equated with exponential suppression of \(\langle0|\hat n|1\rangle\).
Flux noise can arise from a surface-spin bath. Its low-frequency component primarily causes dephasing through frequency sensitivity:
\[ \text{dephasing weight} \propto \left(\frac{\partial\omega_{01}}{\partial\Phi_{\rm ext}}\right)^2 S_\Phi(\omega\approx0). \]
Noise in the same flux channel near \(\omega_{01}\) instead causes transitions through the off-diagonal matrix element:
\[ \Gamma_{1\to0}^{(\Phi)} =\frac{1}{\hbar^2} \left| \left\langle0\left| \frac{\partial\hat H_0}{\partial\Phi_{\rm ext}} \right|1\right\rangle \right|^2 S_\Phi(+\omega_{01}). \]
A sweet spot and relaxation protection are therefore not the same statement. The former controls diagonal sensitivity; the latter also requires examining the finite-frequency spectrum and off-diagonal matrix element.
Under the weak-coupling, Markovian exponential-decay approximation, the full derivation can be compressed to:
\[ \boxed{ \text{circuit and bias} \Rightarrow \hat H_0(\lambda) \Rightarrow \{E_n,\ket n\} \Rightarrow \hat O_\lambda=\partial_\lambda\hat H_0 \Rightarrow \begin{cases} O_{nn}\ \text{and}\ S_\lambda(0) &\Rightarrow T_\phi,\\ O_{mn}\ \text{and}\ S_\lambda(\omega_{nm}) &\Rightarrow T_1\ \text{and excitation}, \end{cases} \Rightarrow T_2^{-1}=(2T_1)^{-1}+T_\phi^{-1}.} \]
For an external perturbation to cause a particular transition, three conditions must hold simultaneously: the environment must couple through a definite \(\hat O\); \(\langle f|\hat O|i\rangle\neq0\); and \(S_F\) must have spectral weight at the required Bohr frequency. “Allowed” means only that symmetry does not forbid the transition, not that it is necessarily strong. “Forbidden” means exactly zero; “suppressed” means nonzero but small. Coherence comparisons must also go beyond a single \(T_1\): matrix-element suppression and reduced environmental spectral density have different implications for device scaling.
Conclusion to Question 1: Decoherence cannot be attributed merely to “a noisy environment.” A complete answer must begin from \(\hat H_0\) and identify the fluctuating parameter; the coupling operator derived from \(\partial_\lambda\hat H_0\); the eigenstates; the diagonal and off-diagonal matrix elements; and the relevant frequency in \(S_F(0)\) or \(S_F(\omega_{01})\).
Question 1 identified two multiplicative factors in a transition: the environmental spectrum and the system matrix element. In Question 2, assume that the environment cannot be eliminated entirely and ask instead whether the wavefunctions, symmetry, and logical encoding of \(\ket0,\ket1\) can be designed so that the dominant noise operator barely sees the logical transition while another controllable operator remains effective.
Environment-induced relaxation and intentional control have the same mathematical skeleton:
\[ \hat V_{\rm noise}(t)=\hat O_{\rm noise}F_{\rm noise}(t), \qquad \Gamma_\downarrow= \frac{1}{\hbar^2} \left|\langle0|\hat O_{\rm noise}|1\rangle\right|^2 S_{\rm noise}(+\omega_{01}), \]
\[ \hat V_{\rm drive}(t)=\hat O_{\rm drive}A\cos\omega_dt, \qquad \Omega_R=\frac{A}{\hbar} \left|\langle0|\hat O_{\rm drive}|1\rangle\right|. \]
The difference between these expressions is not that “noise is bad and drive is good,” but whether they couple through the same physical port and operator. If \(\hat O_{\rm noise}=\hat O_{\rm drive}\), reducing the qubit matrix element simultaneously lengthens \(T_1\) and lowers the Rabi rate. Only when the two have different symmetry, spatial profiles, or frequency dependence can coupling be weak to the environment yet strong to control.
Define the logical projector
\[ \hat P_{\rm L}=\ket0\!\bra0+\ket1\!\bra1. \]
Any Hermitian operator can be decomposed in the logical subspace as
\[ \hat P_{\rm L}\hat O\hat P_{\rm L} =\bar O\,\hat I +\frac{O_{11}-O_{00}}{2}\hat\sigma_z +\operatorname{Re}O_{01}\,\hat\sigma_x -\operatorname{Im}O_{01}\,\hat\sigma_y, \qquad \bar O=\frac{O_{00}+O_{11}}{2}. \]
The physical requirements for state protection can therefore be read directly as:
The ideal case can be summarized as
\[ \boxed{ \hat P_{\rm L}\hat O_{\rm noise}\hat P_{\rm L} \simeq c\hat I, \qquad \left|\langle0|\hat O_{\rm drive}|1\rangle\right| \neq0.} \]
This condition must be applied separately to every dominant noise channel. Holding for charge noise does not imply that it also holds for flux, critical-current, dielectric-loss, or quasiparticle channels.
Suppose a parity-type unitary symmetry \(\hat{\mathcal P}\) exists and satisfies \(\hat{\mathcal P}^2=\hat I\) and \([\hat H_0,\hat{\mathcal P}]=0\). The energy eigenstates can be chosen simultaneously as symmetry eigenstates:
\[ \hat{\mathcal P}\ket n=p_n\ket n, \qquad p_n=\pm1. \]
If the coupling operator has parity \(\eta_O=\pm1\) under this symmetry,
\[ \hat{\mathcal P}\hat O\hat{\mathcal P}^{\dagger}=\eta_O\hat O, \]
then the matrix element must satisfy
\[ \begin{aligned} O_{mn} &=\langle m|\hat O|n\rangle\\ &=\langle m|\hat{\mathcal P}^{\dagger} (\hat{\mathcal P}\hat O\hat{\mathcal P}^{\dagger}) \hat{\mathcal P}|n\rangle\\ &=p_mp_n\eta_O\,O_{mn}. \end{aligned} \]
Therefore,
\[ \boxed{p_mp_n\eta_O=-1\quad\Longrightarrow\quad O_{mn}=0.} \]
This is the first way to separate noise from control: give \(\hat O_{\rm noise}\) a forbidden parity and \(\hat O_{\rm drive}\) an allowed parity. In the symmetric half-flux fluxonium potential, states have definite parity. Note carefully, however, that \(\hat n\) and the phase operator relative to the symmetry center are both odd operators, so the \(0\leftrightarrow1\) transition between states of opposite parity is parity-allowed. If that matrix element is small, the physical origin is wavefunction or tunneling suppression; it must not be mislabeled parity-forbidden.
If the real device contains a small symmetry-breaking term,
\[ \hat H=\hat H_0+\epsilon\hat W, \qquad [\hat W,\hat{\mathcal P}]\neq0, \]
first-order perturbation theory gives
\[ \ket{\widetilde n} \simeq\ket n+\epsilon\sum_{k\neq n} \frac{\langle k|\hat W|n\rangle}{E_n-E_k}\ket k. \]
A \(\langle0|\hat O|1\rangle\) that was exactly zero becomes \(O(\epsilon)\) through opposite-symmetry admixture, and the corresponding transition rate generally begins at \(O(\epsilon^2)\) because it depends on the squared matrix element. Thus, a claim of “symmetry protection” must also specify how flux detuning, junction mismatch, stray capacitance, offset charge, or control-line asymmetry breaks the symmetry; showing only the ideal Hamiltonian is insufficient.
In the representation of generalized coordinate \(x\),
\[ O_{01}=\langle0|\hat O|1\rangle =\int dx\,\psi_0^*(x)\,\hat O(x,-i\partial_x)\,\psi_1(x). \]
If logical wavefunctions are localized in different wells or regions, or have nearly disjoint support in a multidimensional configuration space, the integral for certain local or smooth operators can be extremely small. This suppression may scale exponentially with the barrier action, effective mass, or an energy-scale ratio, but it is usually not exactly zero: finite tunneling, disorder, mode mixing, and the derivative structure of the operator can all leave a residual matrix element.
The two statements must therefore be kept separate:selection-rule forbidden A zero obtained from symmetry algebra; overlap suppressed A small value obtained from an explicit eigenfunction integral. The former is governed primarily by symmetry-breaking perturbations; the latter by parameter scaling and fabrication variation.
A multi-degree-of-freedom or encoded-state structure can also make a dominant fluctuation act in common mode, so that \(\hat P_{\rm L}(\sum_j c_j\hat O_j)\hat P_{\rm L}\simeq c\hat I\). This protects only the specified collective channel; local disorder or differential noise that projects onto \(\sigma_z\) or \(\sigma_{x,y}\) still causes decoherence.
For parameter noise \(\lambda\), the Hellmann–Feynman relation from Question 1 gives
\[ \frac{\partial\omega_{01}}{\partial\lambda} =\frac{1}{\hbar} \left( \left\langle1\left|\frac{\partial\hat H_0}{\partial\lambda}\right|1\right\rangle - \left\langle0\left|\frac{\partial\hat H_0}{\partial\lambda}\right|0\right\rangle \right). \]
A sweet spot \(\partial_\lambda\omega_{01}=0\) means that the two logical states have the same first-order diagonal response to that fluctuation, making them harder for the environment to distinguish through a frequency shift. It does not require \(\langle0|\partial_\lambda\hat H_0|1\rangle=0\) and therefore does not automatically imply a longer \(T_1\). Conversely, a small off-diagonal matrix element does not guarantee that \(O_{11}-O_{00}\) is small.
Three forms of protection must therefore be distinguished:
Preserving controllability does not require making the protected \(0\leftrightarrow1\) transition strong. Common methods include:
For a detuned auxiliary state \(\ket a\), eliminating \(\ket a\) gives an effective Raman coupling of the form
\[ \Omega_{\rm eff}\sim \frac{\Omega_{0a}\Omega_{a1}}{2\Delta_a}, \]
The costs include auxiliary-state population, AC Stark shifts, additional decay paths, and calibration complexity. Temporarily breaking symmetry is not free either: a control pulse may also open a previously forbidden noise channel.
Dispersive readout can likewise use virtual transitions through higher states without requiring a strong direct qubit transition. Under the rotating-wave and dispersive approximations, it can be represented schematically as
\[
\chi_j\simeq
\sum_{k
This route still requires checking the critical photon number, Purcell decay, near-resonant higher transitions, and measurement-induced leakage.
| Noise channel | Primary transmon strategy | Strategies available to fluxonium |
|---|---|---|
| Charge noise | Large \(E_J/E_C\) suppresses charge dispersion of the energy with respect to \(n_g\) | A multiwell wavefunction structure can modify the \(\hat n\) matrix element |
| Flux noise | A fixed-frequency design does not directly use flux bias; a tunable design requires a separate sweet-spot analysis | The half-flux sweet spot makes \(\partial\omega_{01}/\partial\Phi=0\) |
| Even-operator perturbation | Determined by the specific symmetry | Half-flux parity can strictly forbid certain transitions |
| Protection nature | Parameter-regime and energy-dispersion engineering | Symmetry, sweet spots, and wavefunction engineering |
| Control operator | Usually controlled directly through a capacitive \(\hat n\) drive on \(0\leftrightarrow1\) | Choose charge drive, flux drive, or a higher transition according to the matrix-element hierarchy |
| Leakage | Weaker anharmonicity constrains drive bandwidth | Stronger anharmonicity helps, but every transition in the rich multilevel spectrum must still be checked |
If \(O_{01}^{\rm drive}\) is very small, increasing amplitude \(A\) can indeed restore \(\Omega_R\), but the same pulse also couples to \(\ket{k\ge2}\). Define the closest leakage detuning from the drive:
\[ \Delta_{\rm leakage} =\min_{i\in\{0,1\},\,k\ge2} \left|\omega_{ki}-\omega_d\right|. \]
For an off-resonant leakage channel, the lowest-order estimate is
\[ P_{\rm leakage} \sim \sum_{i\in\{0,1\},\,k\ge2} \left| \frac{A\langle k|\hat O_{\rm drive}|i\rangle} {\hbar(\omega_{ki}-\omega_d)} \right|^2. \]
A useful state structure must therefore do more than make \(O_{01}^{\rm noise}\) small: it must create a sufficiently large hierarchy between the desired drive matrix element and undesired leakage matrix elements. Direct resonant control generally requires \(\Omega_R\ll\Delta_{\rm leakage}\); an actual gate must also account for the pulse spectrum, AC Stark shift, counter-rotating terms, and calibration error.
Conclusion to Question 2: An ideal state structure makes the dominant noise operator approximately the identity in the logical subspace while retaining a transverse component in a selectable drive operator. An implementation must jointly verify the origin of symmetry or overlap, sensitivity to symmetry breaking, Rabi rate, auxiliary paths, and leakage. Protection is operator-specific state engineering, not the elimination of all coupling.
Fluxonium can be viewed as a rearrangement of charge- and flux-qubit design. It does not merely substitute a different set of numbers into the transmon Hamiltonian; it adds a very large inductance, creating new global confinement, flux dependence, and a multiwell potential landscape. Compared with the transmon, fluxonium commonly has a lower computational frequency, greater anharmonicity, and a longer relaxation time at selected operating points.
A superinductor must combine very large inductance, low loss, and low parasitic capacitance. It may be implemented with a Josephson-junction array, a high-kinetic-inductance material, or another high-impedance structure. These implementations also introduce array modes, process variation, and parasitic coupling.
Choose the branch flux across the small junction as the generalized coordinate:
\[ \Phi(t)=\int^t V(t')dt', \qquad \varphi=\frac{2\pi\Phi}{\Phi_0}, \qquad \varphi_{\rm ext}=\frac{2\pi\Phi_{\rm ext}}{\Phi_0}. \]
The capacitor, Josephson junction, and linear inductor respectively contribute
\[ T_C=\frac12C_\Sigma\dot\Phi^2, \qquad U_J=-E_J\cos\varphi, \qquad U_L=\frac{(\Phi-\Phi_{\rm ext})^2}{2L}. \]
The location of external flux in this expression depends on the gauge and coordinate choice. The physical content is that the fluxoid constraint makes the small-junction phase, superinductor phase, and external flux mutually dependent.
The Lagrangian is kinetic energy minus potential energy:
\[ \mathcal L =\frac12C_\Sigma\dot\Phi^2 +E_J\cos\!\left(\frac{2\pi\Phi}{\Phi_0}\right) -\frac{(\Phi-\Phi_{\rm ext})^2}{2L}. \]
The canonical charge conjugate to \(\Phi\) is
\[ Q=\frac{\partial\mathcal L}{\partial\dot\Phi} =C_\Sigma\dot\Phi, \qquad n=\frac{Q}{2e}, \qquad [\hat\Phi,\hat Q]=i\hbar \Longleftrightarrow [\hat\varphi,\hat n]=i. \]
Perform the Legendre transform:
\[ \begin{aligned} \hat H &=\hat Q\dot{\hat\Phi}-\mathcal L\\ &=\frac{\hat Q^2}{2C_\Sigma} -E_J\cos\hat\varphi +\frac{(\hat\Phi-\Phi_{\rm ext})^2}{2L}\\ &=4E_C\hat n^2 -E_J\cos\hat\varphi +\frac12E_L(\hat\varphi-\varphi_{\rm ext})^2, \end{aligned} \]
where
\[ E_C=\frac{e^2}{2C_\Sigma}, \qquad E_L=\frac{1}{L}\left(\frac{\Phi_0}{2\pi}\right)^2. \]
The most direct answer to Question 3 is therefore:A superinductor adds a quadratic inductive-energy term \(\tfrac12E_L(\hat\varphi-\varphi_{\rm ext})^2\). More importantly and counterintuitively, larger \(L\) means smaller \(E_L\). A superinductor supplies global confinement with low curvature that extends across many cosine wells.
Define the new coordinate as \(\hat\theta=\hat\varphi-\varphi_{\rm ext}\); the Hamiltonian can be rewritten as
\[ \hat H =4E_C\hat n^2 +\frac12E_L\hat\theta^2 -E_J\cos(\hat\theta+\varphi_{\rm ext}). \]
Thus the two forms commonly seen in the literature
\[ -E_J\cos\hat\varphi +\frac12E_L(\hat\varphi-\varphi_{\rm ext})^2 \qquad\text{and}\qquad -E_J\cos(\hat\theta+\varphi_{\rm ext}) +\frac12E_L\hat\theta^2 \]
differ only by a coordinate shift. A derivation or numerical calculation must use one gauge consistently from beginning to end; the external-flux shift must not be added to both the cosine and the parabola.
A transmon-like circuit can be written as
\[ \hat H_{\rm T}=4E_C(\hat n-n_g)^2-E_J\cos\hat\varphi. \]
Without an inductive term, the potential is \(2\pi\)-periodic, phase is a compact coordinate, and \(n_g\) enters the wavefunction boundary condition. Adding the inductor makes \(\varphi\) an extended coordinate on the full real line: \(\varphi\to\varphi+2\pi\) moves to another cosine well and is no longer the same configuration point.
If offset charge is retained temporarily,
\[ \hat H=4E_C(\hat n-n_g)^2+U(\hat\varphi), \qquad \hat U_g=e^{in_g\hat\varphi}, \qquad \hat U_g^\dagger(\hat n-n_g)\hat U_g=\hat n. \]
Because \(\varphi\in\mathbb R\), this unitary transformation does not introduce the twisted boundary condition of a compact coordinate, and static \(n_g\) can be removed from the low-energy Hamiltonian. This is the structural origin of fluxonium being “free of charge offsets.” It does not mean that every capacitive-loss channel or finite-frequency noise coupled through \(\hat n\) disappears.
In practice, \(N\) larger Josephson junctions are often connected in series. Let the phase drop across each array junction be \(\theta_j\), with \(E_{J,a}\gg E_{C,a}\) so that quantum phase slips in each junction are rare. Expanding for small \(\theta_j\):
\[ -E_{J,a}\cos\theta_j \simeq-E_{J,a}+\frac12E_{J,a}\theta_j^2. \]
If the collective phase drop is \(\Theta=\sum_{j=1}^{N}\theta_j\), the lowest-energy distribution is approximately \(\theta_j=\Theta/N\), so
\[ U_{\rm array} \simeq\sum_{j=1}^{N}\frac12E_{J,a}\left(\frac{\Theta}{N}\right)^2 =\frac12\frac{E_{J,a}}{N}\Theta^2 \equiv\frac12E_L\Theta^2. \]
Therefore, in the simplest uniform-array approximation,
\[ E_L\simeq\frac{E_{J,a}}{N}, \qquad L_{\rm array}\simeq N\frac{(\Phi_0/2\pi)^2}{E_{J,a}}. \]
Adding series junctions increases total inductance and lowers effective \(E_L\). But \(N\) cannot be increased arbitrarily: junction capacitance, island-to-ground capacitance, parameter disorder, and array length together create a family of internal plasma or array modes.
The effective potential is
\[ U(\varphi) =-E_J\cos\varphi +\frac12E_L(\varphi-\varphi_{\rm ext})^2. \]
Local minima \(\varphi_m\) satisfy
\[ E_J\sin\varphi_m +E_L(\varphi_m-\varphi_{\rm ext})=0, \]
The curvature near each minimum and the small-oscillation energy scale are approximately
\[ U''(\varphi_m)=E_J\cos\varphi_m+E_L, \qquad \hbar\omega_{p,m}\simeq \sqrt{8E_C\big(E_J\cos\varphi_m+E_L\big)}. \]
When \(E_J/E_L\) is large, the cosine produces many local wells while the slowly varying parabola sets their global energy ordering. Ignoring local zero-point corrections,
\[ U_m\approx-E_J+\frac12E_L(2\pi m-\varphi_{\rm ext})^2. \]
External flux can therefore reorder different fluxoid-like wells. At \(\varphi_{\rm ext}=\pi\), neighboring wells can become a symmetric degenerate pair. \(E_C\) and the barrier determine interwell tunneling—the wavefunction engineering pursued in Question 4.
A real array is not an ideal coil without dynamics. When all array-island fluxes are retained, the Hamiltonian contains a capacitance matrix and multiple Josephson cosines:
\[ \hat H_{\rm array} =\frac12\hat{\mathbf Q}^{\mathsf T}\mathbf C^{-1}\hat{\mathbf Q} -\sum_{j=1}^{N}E_{J,a}^{(j)}\cos\hat\theta_j. \]
Reducing it to one quadratic term requires at least:
Once these conditions fail, “adding a superinductor” no longer means adding only \(\tfrac12E_L\varphi^2\). It also adds parasitic modes, avoided crossings, additional decay and leakage paths, and fabrication sensitivity.
The goal of superinductance is zero DC resistance, low dissipation, and an impedance above the Cooper-pair resistance quantum \(R_Q=h/(2e)^2\approx6.45\,\mathrm{k}\Omega\) across the operating band, while pushing self-resonance outside that band. A large resistor may also have high impedance, but it introduces real admittance, Johnson–Nyquist noise, and energy dissipation, and cannot provide the same coherent Hamiltonian term.
Core statement: Fluxonium is not merely a different parameter set; it uses a superinductor to add global quadratic confinement to the Hamiltonian. The ideal model adds \(\tfrac12E_L(\hat\varphi-\varphi_{\rm ext})^2\). A real device must additionally verify that phase slips, self-resonance, internal modes, and loss in the JJ array can be safely eliminated.
| Layout/device structure | Primary parameter | Further impact |
|---|---|---|
| Capacitor pad geometry | Total capacitance and \(E_C\) | Phase fluctuation, anharmonicity, coupling capacitance |
| Small-junction area/oxide | \(I_c\) and \(E_J\) | Cosine-well depth, spectrum, fabrication spread |
| JJ array or high-impedance structure | \(L\) and \(E_L\) | Global confinement, array parasitic modes, loss |
| Array junction count/area | \(E_L\), phase-slip amplitude | Inductance, linearity, array length, and accumulated process error |
| Array island ground capacitance | Internal-mode frequencies | Self-resonance, mode hybridization, readout/drive leakage |
| Flux loop/flux line | \(\Phi_{\rm ext}\) | Sweet-spot position, bias precision, bandwidth, flux noise |
Chip design does not merely draw a circuit; it directly sculpts the Hamiltonian's parameter space. At the same time, process variation becomes \(E_J,E_C,E_L\) spread, frequency offsets, and yield problems.
Classify the symbols in this question into four groups: selectable design parameters (\(E_J,E_L,E_C,\Phi_{\rm ext}\)); fabrication/device-generated parameters (junction area, oxide thickness); parameters that must be calibrated (actual \(\omega_{01}\), sweet-spot location); and states that can only be estimated by measurement (thermal population, residual flux offset).
Question 3 derived the fluxonium Hamiltonian from the circuit. Question 4 goes beyond asking “what term was added” and tracks how each energy scale changes the eigenvalue problem:
\[ \hat H =4E_C\hat n^2 -E_J\cos\hat\varphi +\frac12E_L(\hat\varphi-\varphi_{\rm ext})^2, \qquad [\hat\varphi,\hat n]=i. \]
In the phase basis, \(\hat n=-i\partial_\varphi\), so
\[ \left[ -4E_C\frac{\partial^2}{\partial\varphi^2} +U(\varphi) \right]\psi_k(\varphi) =E_k\psi_k(\varphi), \]
where
\[ U(\varphi)=-E_J\cos\varphi +\frac12E_L(\varphi-\varphi_{\rm ext})^2. \]
Comparing with an ordinary one-dimensional Schrödinger equation, the effective mass of the phase particle is
\[ m_\varphi=\frac{\hbar^2}{8E_C}. \]
Thus \(E_C\) is more than a “charging-energy number”: it sets the quantum kinetic energy and effective mass of the phase particle, while \(E_J,E_L,\varphi_{\rm ext}\) jointly determine the potential landscape it sees.
| Energy scale | Physical origin | Role in fluxonium |
|---|---|---|
| \(E_C\) | Capacitance | Sets quantum fluctuations, spread, and interwell tunneling of the phase wavefunction |
| \(E_J\) | Josephson junction | Forms periodic cosine wells and controls well depth and nonlinearity |
| \(E_L\) | Superinductor | Provides slowly varying parabolic confinement and sets the relative energies of different wells |
| \(\Phi_{\rm ext}\) | External flux | Shifts the parabolic term and controls potential symmetry and the operating point |
Fluxonium's physics comes not from one parameter but from competition among \(E_C,E_J,E_L\). Tuning them simultaneously changes wavefunction localization, tunneling, transition frequencies, anharmonicity, operator matrix elements, and the leakage spectrum.
Let \(\varphi_m\) denote a local minimum and define
\[ K_m=U''(\varphi_m)=E_J\cos\varphi_m+E_L. \]
A harmonic expansion around \(\varphi_m\) gives
\[ \hbar\omega_{p,m}\simeq\sqrt{8E_CK_m}, \qquad \varphi_{{\rm zpf},m}=\left(\frac{2E_C}{K_m}\right)^{1/4}, \qquad n_{{\rm zpf},m}=\left(\frac{K_m}{32E_C}\right)^{1/4}. \]
These three quantities provide direct parameter intuition:
In the deep-well regime, the tunnel splitting can be written schematically as \(\Delta\sim\hbar\omega_p e^{-S}\), where the semiclassical action \(S\) increases with barrier height and effective mass. Larger \(E_J/E_C\) therefore generally localizes the states more strongly and reduces \(\Delta\); the exact value still requires solving the full potential.
When \(E_J/E_L\gg1\), local minima lie near \(\varphi_m\simeq2\pi m\) and their localized-state energies can first be estimated as
\[ E_m^{\rm loc} \approx -E_J +\frac12E_L(2\pi m-\varphi_{\rm ext})^2 +\frac12\hbar\omega_{p,m}. \]
The slope with respect to external flux corresponds to the persistent current:
\[ I_m=-\frac{\partial E_m^{\rm loc}}{\partial\Phi_{\rm ext}} \approx \frac{2\pi E_L}{\Phi_0}(2\pi m-\varphi_{\rm ext}). \]
When the energies of two localized wells approach each other, tunneling turns the crossing into an avoided crossing. External flux controls well detuning; \(E_J/E_C\) primarily controls the tunnel gap; and \(E_L\) controls both the detuning scale and the persistent-current slope.
At \(\varphi_{\rm ext}=\pi\), let \(x=\varphi-\pi\). The Hamiltonian becomes
\[ \hat H_{1/2} =-4E_C\frac{\partial^2}{\partial x^2} +E_J\cos x+\frac12E_Lx^2, \qquad U(x)=U(-x). \]
The potential has exact reflection symmetry, so eigenstates can have definite parity: the ground state \(\ket0\) is even and the first excited state \(\ket1\) is odd. The bonding/antibonding picture provides intuition, but a rigorous selection rule follows from symmetry and operator parity.
Retaining only the left and right localized states \(\ket L,\ket R\), the low-energy Hamiltonian can be written as
\[ \hat H_{\rm 2w} =-\frac12\left(\Delta\hat\sigma_x+\epsilon\hat\sigma_z\right), \qquad \epsilon\simeq2I_p\,\delta\Phi, \qquad \delta\Phi=\Phi_{\rm ext}-\frac{\Phi_0}{2}. \]
Its qubit splitting and half-flux eigenstates are
\[ \hbar\omega_{01}=\sqrt{\Delta^2+\epsilon^2}, \qquad \ket0\simeq\frac{\ket L+\ket R}{\sqrt2}, \qquad \ket1\simeq\frac{\ket L-\ket R}{\sqrt2}. \]
Away from half flux, \(|\epsilon|\) increases and the eigenstates gradually localize; at half flux they form even and odd superpositions. This two-well model supplies intuition, but higher plasmon states, other wells, and actual matrix elements must be calculated from the complete Hamiltonian.
Exact parity follows from \(U(x)=U(-x)\) at half flux. Conditions such as \(E_J\gg E_L\) and \(E_J\gg E_C\) are useful in the deep-well/bonding picture but are not required for parity symmetry. Away from half flux, symmetry is broken and first-order flux sensitivity returns, with a magnitude that depends on device parameters.
The external-flux coupling operator is
\[ \frac{\partial\hat H}{\partial\Phi_{\rm ext}} =-\frac{2\pi E_L}{\Phi_0} (\hat\varphi-\varphi_{\rm ext}). \]
At half flux, \(x=\varphi-\pi\) is an odd operator, and every parity eigenstate satisfies \(\langle n|x|n\rangle=0\). By the Hellmann–Feynman theorem,
\[ \left.\frac{\partial E_n}{\partial\Phi_{\rm ext}}\right|_{\Phi_0/2}=0, \qquad \boxed{ \left.\frac{\partial\omega_{01}}{\partial\Phi_{\rm ext}} \right|_{\Phi_0/2}=0.} \]
The two-well model gives the same conclusion:
\[ \frac{\partial\omega_{01}}{\partial\Phi_{\rm ext}} =\frac{1}{\hbar} \frac{\epsilon}{\sqrt{\Delta^2+\epsilon^2}} \frac{\partial\epsilon}{\partial\Phi_{\rm ext}} \xrightarrow[\epsilon\to0]{}0. \]
At this operating point, first-order flux sensitivity vanishes and dephasing is dominated by higher-order terms. But a sweet spot does not automatically ensure relaxation protection; relaxation still requires the relevant operator matrix element and \(S_F(\omega_{01})\) to be checked.
Moreover, in the two-well approximation,
\[ \left.\frac{\partial^2\omega_{01}}{\partial\Phi_{\rm ext}^2} \right|_{\Phi_0/2} \simeq\frac{(2I_p)^2}{\hbar\Delta}. \]
Making the tunnel gap \(\Delta\) extremely small can reduce \(\omega_{01}\) and some transition matrix elements, but it increases the second-order curvature at the sweet spot and aggravates thermal-initialization and slow-control problems. Protection remains a multiobjective trade-off.
| Fluctuation | System operator(half flux) | Parity | Direct conclusion |
|---|---|---|---|
| External-flux noise | \(x=\varphi-\pi\) | Odd | The diagonal element vanishes, providing a first-order dephasing sweet spot; \(0\leftrightarrow1\) relaxation is allowed |
| Capacitive/electric loss | \(\hat n\) | Odd | \(0\leftrightarrow1\) is allowed, but its matrix element may be suppressed for a low-frequency interwell transition |
| Critical-current/\(E_J\) noise | \(-\cos\varphi\) | Even | Transitions between opposite parities are forbidden; diagonal sensitivity must still be calculated |
| Inductive-energy/\(E_L\) noise | \(\tfrac12x^2\) | Even | Transitions between opposite parities are forbidden; frequency fluctuations may remain |
Fluxonium's low-energy states can occupy different phase wells. If two wavefunctions have little overlap, transition matrix elements of some operators also shrink, reducing the probability of \(\ket1\rightarrow\ket0\) through that noise channel. This is wavefunction engineering and must be named separately from an exact parity selection rule.
For Hamiltonian \(\hat H=4E_C\hat n^2+U(\hat\varphi)\), the commutator provides a very useful exact relation:
\[ [\hat H,\hat\varphi]=-8iE_C\hat n, \qquad \left|\langle0|\hat n|1\rangle\right| =\frac{\hbar\omega_{01}}{8E_C} \left|\langle0|\hat\varphi|1\rangle\right|. \]
This shows that, even when the half-flux \(0\leftrightarrow1\) transition has a sizable phase dipole, low \(\omega_{01}\) still suppresses the charge matrix element. Capacitive/dielectric relaxation can therefore be suppressed, whereas relaxation through a flux-like operator does not automatically receive the same suppression; each must still be multiplied by its own \(S_{\rm noise}(\omega_{01})\).
For capacitive and flux-like drives, the direct Rabi rates respectively take the forms
\[ \Omega_R^{(n)}\propto A_n\left|\langle0|\hat n|1\rangle\right|, \qquad \Omega_R^{(\varphi)}\propto A_\varphi\left|\langle0|x|1\rangle\right|. \]
A small charge matrix element slows a direct capacitive gate. Higher transitions in fluxonium, such as \(1\leftrightarrow2\), may have much larger charge matrix elements, so control and readout can use the operator hierarchy, an auxiliary level, or a different coupling port. The cost is more complicated pulse sequences, frequency allocation, and leakage analysis.
At minimum, compute simultaneously
\[ \alpha=\omega_{12}-\omega_{01}, \qquad \Delta_{\rm leakage} =\min_{i\in\{0,1\},\,k\ge2} |\omega_{ki}-\omega_d|, \qquad \Omega_R\ll\Delta_{\rm leakage}. \]
Fluxonium often has greater anharmonicity but also a rich multilevel spectrum. “Large anharmonicity” cannot replace checking every drive matrix element, near resonance, array mode, and multiphoton leakage channel.
Fluxonium's protected regime often comes with a lower computational frequency. For example, with \(f_{01}=500\,\mathrm{MHz}\), \(hf/k_B\approx24\,\mathrm{mK}\): at \(T=50\,\mathrm{mK}\), \(p_1/p_0\approx e^{-24/50}\approx0.62\) and \(p_1\approx0.38\); at \(T=20\,\mathrm{mK}\), \(p_1/p_0\approx0.30\) and \(p_1\approx0.23\).
Another useful scale is \(f_{01}=800\,\mathrm{MHz}\), for which \(hf/k_B\approx38\,\mathrm{mK}\); at \(20\,\mathrm{mK}\), \(p_1\approx0.13\) (only a thermal-order estimate; the actual value also depends on nonequilibrium photons and effective temperature). A low-frequency transition turns active reset, thermalization, and additional readout transitions into design problems.
Even if the mixing chamber reads 10 mK, a 500 MHz fluxonium cannot be assumed to initialize naturally in the ground state if its effective temperature is 50 mK. The same device may have both a low-frequency computational transition and higher-frequency transitions useful for readout and control. One must therefore choose which transition to use, prevent spectator excitation, and verify that the electronics and resonator bandwidth support the plan.
The first four questions used fluxonium to build the complete chain from noise source and operator to Hamiltonian, eigenstates, and matrix elements. Question 5 treats fluxonium as the first complete case and reverses the direction: If the dominant error operator is specified first, should we change parameters, add degrees of freedom, create a new symmetry, or modify the Josephson element itself?
\[ \boxed{ \begin{aligned} \text{Dephasing protection:}\quad &\frac{\partial\omega_{01}}{\partial\lambda}\approx0,\\[4pt] \text{Relaxation protection:}\quad &\langle0|\hat O_{\rm noise}|1\rangle\approx0,\\[4pt] \text{Controllability:}\quad &\langle0|\hat O_{\rm drive}|1\rangle\neq0,\\[4pt] \text{Leakage suppression:}\quad &\Delta_{\rm leakage}\ \text{sufficiently large}. \end{aligned}} \]
These are not four optional advantages; they form a design specification that every protected-qubit proposal must answer in full:
Reverse-design chain: error operator → desired eigenstate structure → Hamiltonian term/symmetry → spectrum and matrix-element test → control/leakage test → circuit/material implementation.
The teaching objective here is explicitly to create entry points into the terminology and concepts. Students need not complete a full derivation of every architecture this week, but they must see which layer of the Hamiltonian each direction modifies. When they later encounter a new qubit proposal, they should be able to analyze it independently with the same four criteria.
Fluxonium can still be viewed primarily as an effective single-degree-of-freedom problem \(H=H(\varphi,n)\). A \(0-\pi\) circuit promotes the problem to a multimode Hamiltonian, for example \(H=H(\phi,\theta)\). The physical circuit uses two Josephson junctions, two large inductors, and two large capacitors. Its goal is not to find another ordinary sweet spot, but to give the logical states nearly disjoint support in generalized-coordinate space while bringing their logical energies close to degeneracy:
\[ \text{disjoint support}\Rightarrow \text{small transition matrix element}, \qquad \text{near degeneracy}\Rightarrow \text{small frequency susceptibility}. \]
Self-study question: Why can adding one degree of freedom enable protection that is difficult to achieve simultaneously with a single phase coordinate?
An ordinary Josephson junction provides \(-E_J\cos\varphi\). More generally, the Fourier components of the current–phase relation can be written as:
\[ U_J(\varphi)=-\sum_kE_k\cos(k\varphi), \qquad U(\varphi)=-E_2\cos(2\varphi). \]
\(\cos\varphi\) corresponds to single-Cooper-pair tunneling; an effective \(\cos(2\varphi)\) may correspond to a two-Cooper-pair process, creating a different periodicity and Cooper-pair parity sectors so that symmetry suppresses some charge-induced transitions. This elevates Hamiltonian engineering from “adding \(L\) and \(C\)” to “directly selecting the Fourier content of the Josephson potential.”
Self-study question: Must we accept only the \(-E_J\cos\varphi\) naturally supplied by a junction, or can its Fourier components be engineered?
\[ E_J\longrightarrow E_J(M) \]
The ferromagnetic Josephson-junction direction does not simply add a linear circuit element. It makes the Josephson energy, or more generally the current–phase relation, depend on a magnetic state. A hybrid ferro-transmon can exploit hysteretic behavior so that magnetic pulses provide approximately digital tuning of qubit frequency. But magnetic degrees of freedom simultaneously introduce new noise sources, including magnetic fluctuations and domain dynamics.
Core reminder: Hamiltonian engineering is not necessarily protection engineering; a new control knob may simultaneously introduce a new noise operator.
| Direction | Primary layer modified | Protection intuition | First cost to investigate |
|---|---|---|---|
| Fluxonium | Add an inductive-energy term | Sweet spot, multi-well wavefunction, matrix-element engineering | Low frequency, thermal population, weak drive, and superinductor modes |
| \(0-\pi\) | Add degrees of freedom and hardware symmetry | Disjoint support + near degeneracy | Multimode leakage, disorder, large \(L/C\) |
| \(\cos2\varphi\) | Modify the periodicity of the Josephson potential | Cooper-pair parity sectors | Residual \(\cos\varphi\), control-induced symmetry breaking |
| Ferro-transmon | Modify a material-dependent Josephson element | A new hysteretic/nonvolatile control knob | New magnetic noise; not automatically a protected qubit |
Evaluating fluxonium with the four criteria requires listing the operating frequency, anharmonicity, sweet spot, dominant noise operator, gate primitive, readout contrast, reset, fabrication complexity, coupler compatibility, and frequency allocation. Comparing only the best \(T_1\) ignores the control stack and yield.
| Dimension | Transmon | Fluxonium |
|---|---|---|
| Hamiltonian | \(4E_C\hat n^2-E_J\cos\hat\varphi\) | \(4E_C\hat n^2+\frac12E_L(\hat\varphi-\varphi_{\rm ext})^2-E_J\cos\hat\varphi\) |
| Primary design strategy | Large \(E_J/E_C\) to reduce charge dispersion | Add a superinductor to shape the phase potential and wavefunction |
| Sweet spot | Evaluate charge and flux sensitivity separately for fixed and tunable designs | \(\Phi_{\rm ext}\approx\Phi_0/2\) flattens first-order flux sensitivity |
| Typical frequency | Usually several GHz | May be as low as 0.5–2 GHz, depending on the design |
| Anharmonicity | Moderate, set approximately by \(E_C\) | Can usually provide greater anharmonicity |
| Thermal population | Usually lower at high frequency | Higher at low frequency; active reset may be required |
| Gate speed | Mature microwave control, still limited by leakage | The drive matrix element may be smaller and requires device-specific design |
| Readout | Mature dispersive-readout techniques | Using \(\ket1\leftrightarrow\ket2\) or another higher transition may require an additional control line, frequency band, and transition selection |
| Primary noise/loss | Charge, dielectric, Purcell, flux(tunable), quasiparticle | Flux, dielectric, inductor/array modes, quasiparticle |
| Fabrication | Single JJ or SQUID; mature fabrication | The superinductor/JJ array adds complexity and parameter spread |
| Scalability | A complete ecosystem for control, readout, and couplers, while frequency collisions, crosstalk, leakage, and packaging constraints remain | Coupler compatibility, frequency allocation, array modes, parameter stability, and control-channel count must be handled together |
Each group proposes one real advantage of fluxonium, one cost arising from the same feature, and one experiment that distinguishes a protection mechanism from an accidentally low-noise sample. For example, the half-flux sweet spot reduces first-order flux sensitivity, but a low-frequency transition increases the burdens of thermal initialization and frequency-band integration.
Validation should not cite only the best single-point \(T_1\). It should compare coherence, gate fidelity, thermal population, and yield across flux biases, times, and multiple devices. Grading asks whether students can close the loop from hypothesis to observable, discriminating experiment, and decision.
Group A uses \(\hat O_{\rm noise}=\hat n\) and Group B uses \(\hat O_{\rm noise}=\hat\varphi\). Each group selects one Hamiltonian-engineering direction and answers:
Under what conditions should Fluxonium be chosen over Transmon?
This question is retained but does not replace this week's fifth question on protected Hamiltonians. W3 records only a provisional answer. After students study control in W4, coupling in W5, readout in W6, thermalization in W11, and fabrication/yield in W14–W15, they revisit it in W16 and form a formal architecture decision.
Conclusion for this week: Fluxonium is not the endpoint of the protected-qubit problem; it is the first complete case of Hamiltonian engineering. Students need not leave class able to design every advanced qubit. They should, however, be able to use the four criteria to identify what a new architecture attempts to protect, how it retains control, where leakage appears, and which new degree of freedom or fabrication step hides the cost.
Cross-layer assessment: Redraw the “abstract model → real hardware” information-flow diagram from the start of the class. The difference between the two diagrams is the evidence of this week's learning. Completion requires a checkable derivation, relation diagram, comparison table, or provisional decision, together with an explanation of its physical assumptions and engineering costs.
Next week: W4 will begin with a driven two-level system and derive the rotating frame and RWA. Students should carry this week's intuition for the qubit Hamiltonian, matrix elements, and leakage gap into W4 to understand why pulse control is constrained by the device spectrum.