📚 课程目录第2章双量子比特耦合的物理基础

2. 双量子比特耦合的物理基础

要分析双量子比特的耦合动力学，第一步不是急着写演化算符，而是把“两个独立的二能级系统如何变成一个四维整体”这件事在物理上说清楚。本节从单比特回顾出发，依次建立张量积空间、耦合机制与哈密顿量框架，为后续的能级劈裂、交换振荡和门操作分析打下基础。

一、从单比特到双比特：张量积空间

单个量子比特是一个二维希尔伯特空间，取计算基 |0⟩、|1⟩，任意态写作 α|0⟩+β|1⟩，满足 |α|²+|β|²=1。物理上它可以是超导transmon的两个最低能级、自旋-1/2 的取向、或囚禁离子的两个内态。单比特的哈密顿量在无驱动时通常写成 H = −(ħω/2)σ_z（以能隙 ħω 为标度），本征态就是 |0⟩ 与 |1⟩。

两个比特放在一起时，复合系统的态空间不是二维加二维，而是张量积：ℋ = ℋ_A ⊗ ℋ_B，维数为 2×2 = 4。基矢取 |00⟩、|01⟩、|10⟩、|11⟩，其中 |01⟩ 表示 A 处于 |0⟩、B 处于 |1⟩。这个约定必须在全文与图示中保持一致：左位对应比特 A，右位对应比特 B。张量积的关键性质是“局域算符只作用于自己那一侧”：σ_z⊗I 只翻转 A 的相位，I⊗σ_x 只翻转 B 的比特。正是这种“各管一边”的结构，使得无耦合时系统可分解，而耦合项一旦出现，就会把四个基矢重新混合，产生纠缠与交换。

二、耦合机制：三种常见物理来源

耦合项的本质是让 A 的状态影响 B 的演化能量。工程上最常见的有三类：

交换相互作用（exchange / XY 型）：形式为 H_ex = (ħJ/2)(σ_x⊗σ_x + σ_y⊗σ_y)。它把 |01⟩ 与 |10⟩ 耦合起来，物理图像是一个激发在两个比特间来回传递，是超导比特通过电容耦合、以及自旋比特通过交换作用实现 iSWAP 类门的核心项。
偶极相互作用（dipole / ZZ 型）：形式为 H_zz = (ħζ/2)(σ_z⊗σ_z)。它不交换激发，只让 |11⟩ 相对 |01⟩、|10⟩ 产生额外能量偏移，表现为条件相位，是构造受控相位门与读取“耦合强度”的直接观测量。
纵向/横向驱动耦合：当用微波驱动时，还会出现 σ_x⊗I 一类的局域项，用于单比特旋转；若驱动频率接近比特间失谐，则可能引入非共振的寄生耦合，需要在设计时抑制。

需要强调：真实器件中交换项与 ZZ 项往往同时存在，且耦合强度 J、ζ 随比特间距、耦合器频率、失谐量变化。因此“耦合机制”不是抽象名词，而是决定后续动力学时间尺度的可调参数。

三、耦合双比特哈密顿量框架

把局域项与耦合项合起来，得到本节的核心表达式（取 ħ=1）：

H = −(ω_A/2) σ_z⊗I − (ω_B/2) I⊗σ_z
    + (J/2)(σ_x⊗σ_x + σ_y⊗σ_y)
    + (ζ/2)(σ_z⊗σ_z)

其中 ω_A、ω_B 为两比特的本征频率，J 为交换耦合强度，ζ 为 ZZ 耦合强度。这个框架的价值在于：它把“器件参数”与“可观测动力学”直接连起来。例如当 J=0、ζ=0 时，四个基矢就是本征态，能级为 ±ω_A/2 ± ω_B/2；一旦 J≠0，|01⟩ 与 |10⟩ 形成简并子空间并被混合，本征态变为 (|01⟩±|10⟩)/√2，能级劈裂大小约为 2J。这正是后续章节分析交换振荡频率的出发点。

约定提醒：本文统一采用 σ_z 本征值 +1 对应 |0⟩、−1 对应 |1⟩ 的写法，因此 H 中的负号使 |0⟩ 能量高于 |1⟩；若换用相反约定，所有频率符号需同步翻转，但物理结论不变。

四、一个可手算的示例：从参数到能级

设某超导双比特示例参数为 ω_A/2π = 5.00 GHz、ω_B/2π = 5.20 GHz、J/2π = 20 MHz、ζ/2π = 2 MHz（均为合理假设的示例值，非实测）。在 |01⟩、|10⟩ 子空间内，忽略 ζ 的一阶影响，有效哈密顿量矩阵为：

H_sub = [ [−ω_B/2,      J/2   ],
          [   J/2,    −ω_A/2  ] ]   （单位：2π·GHz）

代入数值，对角元分别为 −2.60 与 −2.50 GHz，非对角元 0.010 GHz。本征值近似为 −2.55 ± 0.05 GHz，即两本征态能量间隔约 100 MHz，对应交换振荡角频率约 2π×100 MHz。这个手算过程说明：失谐 Δ=ω_A−ω_B 决定子空间的对角劈裂，J 决定混合强度；当 |Δ|≫J 时混合被抑制，激发基本锁在单个比特上，交换振荡幅度变小、频率升高。工程上正是通过调节失谐来“开关”耦合。

五、框架的使用要点与常见误区

拿到上述哈密顿量后，实际分析通常分三步走。第一步，确定子空间：若只关心单激发动力学，可把 |00⟩ 与 |11⟩ 投影掉，得到 2×2 有效矩阵；若关心条件相位，则必须保留 |11⟩。第二步，对角化或做旋转波近似：当 |Δ|≫J 时，用 Schrieffer-Wolff 变换可把交换项等效为 ZZ 项，得到 J²/Δ 量级的残余耦合，这正是“耦合关断”时仍存在寄生 ZZ 的来源。第三步，核对量纲与符号：所有频率以 2π·GHz 或 MHz 给出时，代入前先统一，避免把角频率与普通频率混用。

常见误区有三：其一，把张量积顺序写反，导致 |01⟩ 与 |10⟩ 的物理含义颠倒，后续振荡相位全错；其二，忽略失谐，直接令 ω_A=ω_B，得到理想 iSWAP，却与器件实际不符；其三，把 ζ 当作可随意忽略的小量，但在多比特链中残余 ZZ 会累积成显著的门误差。建议在建模时先写出完整四维矩阵，再按需投影，而不是一开始就凭直觉删项。

至此，双比特耦合的物理基础已经搭好：二维张量积给出四维空间，交换与偶极两类耦合给出混合与相位两种效应，哈密顿量框架把器件参数映射为可计算的能级与振荡频率。后续的动力学分析，都将在这个框架上展开。