1. 课程导入与数值分析全景

超导量子比特的能级结构，是理解其量子行为的起点。无论是设计一个Transmon比特、评估其相干性能，还是标定其工作频率与门操作条件，最终都绕不开对能级间距、能级分布以及能级对外部参数响应的定量把握。实验上我们通过微波谱测到的是频率，而理论上真正决定量子比特“好不好用”的，是能级本身如何随电路参数变化。数值分析正是连接这两端的桥梁：它把哈密顿量、电路参数与可观测频率之间的关系，变成可以计算、可以扫描、可以验证的数值结果。

本节作为全课程的开篇，不急于推导具体方程，而是先建立一张“全景图”：我们要分析哪些物理量、这些量之间如何关联、数值任务按什么顺序推进、每一步的输出如何被下一步使用。只有把任务路线理清楚，后续的建模、离散化、求解与验证才不会变成零散公式的堆砌。

为什么超导量子比特的能级必须靠数值方法

超导量子比特本质上是非线性量子谐振电路。以最常用的Transmon为例，其核心是一个约瑟夫森结并联一个大电容，哈密顿量可写成充电能与约瑟夫森能量竞争的形式。若把相位算符按电荷基展开，得到的是一组耦合的差分方程，而不是解析可解的简谐振子。即便在Transmon近似下，能级间距也带有非谐性修正，高阶能级、能级随磁通或栅压的调谐、以及多能级泄漏，都需要在截断的希尔伯特空间中数值求解本征值问题。

更实际的是，真实器件存在参数离散性：结电阻、结电容、并联电容都有制造误差，环境磁通也有漂移。我们关心的往往不是某一组参数下的一个能级值，而是一段参数区间内能级如何移动、非谐性是否足够大、工作点是否落在对磁通噪声不敏感的位置。这类“参数—能级”关系只能通过批量数值扫描获得。因此，数值分析在超导量子比特研究中不是可选项，而是从器件设计到实验标定的常规工具。

本节涉及的典型物理量

在进入具体计算前，需要明确我们要盯住哪些量。它们大致可以分为三类。第一类是电路与材料参数：约瑟夫森能量 E_J、充电能量 E_C、并联电容 C、结电阻 R_n、临界电流 I_c，以及外部磁通偏置 Φ。第二类是由哈密顿量导出的能级量：基态与第一激发态能量 E_0、E_1，跃迁频率 f_01=(E_1-E_0)/h，非谐性 α=(E_2-E_1)-(E_1-E_0)，以及更高能级 E_2、E_3 用于评估泄漏。第三类是派生指标：能级对磁通的敏感度、电荷色散、以及由能级差换算出的门操作频率。

这些量之间存在明确的依赖链：电路参数决定哈密顿量矩阵元，矩阵元决定本征值，本征值之差决定频率与非谐性，频率与非谐性又反过来约束器件设计是否满足“非谐性远大于线宽、工作频率落在微波频段”的工程要求。数值任务的组织，本质上就是沿着这条依赖链安排计算与验证。

数值任务路线图

把上面的依赖链落成可执行的任务，大致分为五步，但每一步都有明确的输入输出，而不是笼统的“建模—仿真”。

第一步，参数整理与单位统一。把设计值或实验反推值整理成同一单位制。常见做法是以频率为单位，把 E_J、E_C 换算成 GHz：E_J/h、E_C/h。例如一组教学示例参数取 E_J/h = 15 GHz、E_C/h = 0.3 GHz，则 E_J/E_C = 50，处于Transmon的典型区间。这一步的输出是一张参数表，后续所有计算都引用它，避免中途换单位导致量级错误。

第二步，构建截断哈密顿量。在电荷基 |n⟩ 下，充电能项是对角的，约瑟夫森项耦合相邻电荷态。取电荷数范围 n = -N … N，矩阵维数为 2N+1。教学示例中先取 N=15，即31维。矩阵元由 E_C 与 E_J 直接给出，磁通偏置通过约瑟夫森项前的系数进入。输出是一个稀疏或稠密矩阵，规模不大时稠密矩阵即可。

第三步，本征值求解与收敛检查。对矩阵做对角化，取最低若干本征值。关键不是算一次，而是把 N 从15加到25、35，观察 f_01 是否稳定。若两次截断间 f_01 变化小于约1 MHz，可认为截断足够；若变化明显，说明波函数尾部被截掉，需要扩大维数。这一步的输出是收敛后的 E_0…E_3 及频率、非谐性。

第四步，参数扫描。固定 E_C，让 E_J/E_C 在一段区间内变化，或固定 E_J 扫描磁通偏置，重复第二、三步，得到 f_01 与非谐性随参数的变化曲线。扫描时要注意保持截断策略一致，否则曲线上的跳变可能来自数值设置而非物理。

第五步，结果判读与实验对照。判读至少看三件事：非谐性是否足够大，使 f_12 与 f_01 可分辨；工作点是否落在磁通敏感度低的位置；高能级泄漏在门操作时间尺度上是否可接受。随后把计算得到的 f_01 与谱线测量值比对，若偏差在几个百分点内，可进一步用测得的 f_01 和 α 反推 E_J、E_C，形成“计算—测量—修正参数”的闭环。

从本节走向后续章节

本节建立的全景图，为后续内容提供了坐标：哈密顿量的具体形式、矩阵元的推导、对角化算法的选择、扫描策略的设计、以及误差来源的分解，都会沿着这张图逐层展开。读者在后续每一节中，都可以回到这张路线图，确认当前计算处在哪一环、输入来自哪里、输出将交给谁使用。能级数值分析的价值，最终体现在它能否支撑一个可制造、可标定、可复现的超导量子比特设计，而这张任务路线正是通向该目标的施工图。