4. Transmon哈密顿量建模

前三节我们建立了超导量子比特的物理直觉：约瑟夫森结提供非线性电感，电容提供动能项，二者并联构成一个非线性LC振荡器。本节要把这个物理图像变成可计算的哈密顿量，并最终得到能级与波函数。核心路线是：写出经典拉格朗日量 → 识别广义坐标与共轭动量 → 正则量子化 → 截断到有限维希尔伯特空间 → 数值对角化。这条路线对Transmon、Fluxonium、0-π等所有超导电路都通用，区别只在于电路拓扑和参数区间。

一、从电路到拉格朗日量：节点磁通作为广义坐标

量子化电路方法的第一步是选广义坐标。对超导电路，最自然的选择是节点磁通Φ_n(t)，定义为该节点相对于地节点的磁通，满足Φ_n = ∫V_ndt。节点电压V_n = dΦ_n/dt，因此节点磁通扮演“位置”的角色，而节点电荷Q_n扮演“动量”的角色。这个约定必须全文一致：我们取Φ为位置、Q为动量，符号不随章节切换。

对最简单的Transmon电路——一个约瑟夫森结（临界电流I_c，结电容C_J）并联一个分流电容C_s——只有一个非地节点。总电容C_Σ = C_J + C_s。电路的拉格朗日量为：

L = (1/2) C_Σ Φ̇² + E_J cos(2πΦ/Φ₀)

其中 Φ₀ = h/(2e) ≈ 2.0678×10⁻¹⁵ Wb 为磁通量子，
E_J = I_c Φ₀/(2π) 为约瑟夫森能量。

注意符号约定：我们取势能项为 +E_Jcos(2πΦ/Φ₀)，这样在Φ=0附近势能取极大值，对应势阱底部在Φ=Φ₀/2处。若采用相反符号，只需在后续定义中同步翻转，物理结果不变。本文统一采用上述约定。

二、正则量子化：共轭动量与哈密顿量

广义动量定义为p = ∂L/∂Φ̇ = C_ΣΦ̇，这正是节点电荷Q。做勒让德变换H = pΦ̇ − L，得到经典哈密顿量：

H = Q²/(2C_Σ) − E_J cos(2πΦ/Φ₀)

正则量子化把Φ和Q提升为算符，并施加对易关系[Φ, Q] = iℏ。为方便，引入无量纲变量：节点磁通相位φ = 2πΦ/Φ₀，节点电荷数n = Q/(2e)。对易关系变为[φ, n] = i。哈密顿量写成：

H = 4E_C n² − E_J cos φ

其中 E_C = e²/(2C_Σ) 为充电能。

这里E_C和E_J是Transmon的两个核心参数。Transmon的设计区间是E_J/E_C ≈ 50～100，此时电荷色散极小，对1/f电荷噪声不敏感，但非线性仍然足够强（约200～300 MHz），可以单独寻址。若E_J/E_C过小，回到电荷量子比特，对电荷噪声敏感；若过大，非线性减弱，泄漏到高能级的问题变严重。

三、截断：从无限维到有限维

φ和n都是连续或无限维算符，数值计算必须截断。两种常用基：

• 电荷基：以n的本征态|n⟩为基，n取整数。φ算符在此基下是升降算符的组合：e^±iφ|n⟩ = |n±1⟩。哈密顿量矩阵为三对角形式，对角元4E_Cn²，次对角元−E_J/2。截断到n ∈ [−N, N]，矩阵维度2N+1。
• 相位基：以φ的离散采样点为基，n算符为−i∂/∂φ。适合势阱形状复杂、需要看波函数空间分布的场景。

截断维度N的选择有明确判据：先取N=15计算前5个能级，再取N=25重算，若最低5个能级变化小于1 MHz，则N足够。对E_J/E_C=50的典型Transmon，N=20～30即可收敛。截断过小会导致高能级被“压扁”，影响门保真度估计。

四、数值对角化：能级、波函数与参数扫描

在电荷基下构造哈密顿量矩阵后，用标准本征值求解器（如LAPACK的zheev或Python的numpy.linalg.eigh）对角化。以下是一个可直接运行的最小示例，参数取E_C/h=0.2 GHz、E_J/h=10 GHz（即E_J/E_C=50），截断N=25：

import numpy as np

E_C = 0.2   # GHz
E_J = 10.0  # GHz
N = 25
dim = 2*N + 1

# 电荷基：n = -N, ..., N
n_vals = np.arange(-N, N+1)
H = np.diag(4*E_C * n_vals**2)

# 次对角元：-E_J/2
for i in range(dim-1):
    H[i, i+1] = -E_J/2
    H[i+1, i] = -E_J/2

evals, evecs = np.linalg.eigh(H)
# 前5个能级（GHz）
print(evals[:5])
# 典型输出（示例，非实测）：
# [0.        , 6.72, 13.18, 19.38, 25.32]

从示例输出可以读出Transmon的关键特征：基态能量设为0，第一激发态约6.72 GHz，第二激发态约13.18 GHz。跃迁频率ω₀₁≈6.72 GHz，ω₁₂≈6.46 GHz，非线性α = ω₁₂ − ω₀₁ ≈ −0.26 GHz，即约−260 MHz。这个负号表示能级间距随激发数增加而减小，是Transmon作为量子比特可寻址的物理基础。若把E_J/E_C降到20，α会增大到约−500 MHz，但电荷噪声敏感度上升；升到100，α缩小到约−130 MHz，泄漏风险增加。设计时需要在噪声敏感度和非线性之间折中。

波函数也可以从本征矢中提取。在电荷基下，第k个本征矢的分量ψ_k(n)给出该能级在电荷数上的分布。对Transmon，基态波函数在n=0附近集中，随E_J/E_C增大而变窄，这正是电荷噪声不敏感的直接体现。若需要看相位空间的波函数，可对ψ_k(n)做离散傅里叶变换得到ψ_k(φ)。

五、参数扫描与设计判断

实际设计中，我们不会只算一组参数。典型做法是固定E_C（由电容决定，通常取0.15～0.3 GHz），扫描E_J/E_C从20到120，记录ω₀₁、α和电荷色散。判断标准有三条：ω₀₁落在微波频段（4～8 GHz）以便与读出腔和微波线匹配；|α|在150～350 MHz之间，既能区分0-1和1-2跃迁，又不会因太小而泄漏；电荷色散∂ω₀₁/∂n_g在n_g=0.5处小于1 MHz，保证对电荷噪声不敏感。这三条同时满足的E_J/E_C区间通常就是50～80。

若需要更精确的模型，可以在哈密顿量中加入读出腔耦合项g(a†+a)n、驱动项Ω(t)n，以及结电容的非线性修正。但作为设计起点，上述单节点模型已经能给出足够准确的ω₀₁和α，后续用微扰或全电路量子化修正即可。数值对角化的优势在于：一旦矩阵构造正确，改变参数只需重新赋值E_C、E_J并重算，几分钟内就能完成一轮参数扫描，比解析近似更可靠。