3. 超导量子电路哈密顿量：LC振荡器、非线性电感、量子化与能级

好，咱们进入正题。超导量子电路的核心，说白了就是怎么用电路元件搭出一个“人工原子”。这个原子得有分立的能级，我们才能用它做量子比特。那第一步，得从最基础的LC振荡器说起。

3.1 经典LC振荡器：谐振子的电路版

一个LC电路，就是一个电感L和一个电容C并联或串联。你想想看，它和力学里的弹簧振子是一模一样的。电荷q在电容上积累，相当于弹簧的位移；电流i流过电感，相当于振子的速度。

经典力学里，谐振子的哈密顿量是：

H = p²/(2m) + (1/2) m ω² x²

换成电路语言，我们把电荷q当作坐标，磁通φ当作动量。嗯，这里有个习惯问题——我个人喜欢用磁通φ作为广义坐标，这样电感上的能量项看起来更顺眼。那么LC振荡器的经典哈密顿量就是：

H = q²/(2C) + φ²/(2L)

你看，第一项是电容上的电能，第二项是电感上的磁能。两者此消彼长，能量在电场和磁场之间来回振荡，频率就是ω = 1/√(LC)。

核心要点：LC振荡器是线性系统，它的能级是等间距的。这意味着你没法单独操控某个能级——这恰恰是我们不想要的。量子比特需要非线性，需要不等间距的能级。

3.2 非线性电感：约瑟夫森结登场

那怎么引入非线性呢？答案就是约瑟夫森结。这东西是超导量子电路的心脏。

约瑟夫森结本质上是一个“非线性电感”。它的电感值不是常数，而是随电流变化的。我当年第一次看到它的本构方程时，觉得这玩意儿真巧妙：

I = I₀ sin(δ)
V = (Φ₀/2π) dδ/dt

这里δ是结两端超导波函数的相位差，I₀是临界电流，Φ₀ = h/2e是磁通量子。从这两个方程，你可以推导出结的势能：

U_J = -E_J cos(δ)

其中E_J = I₀Φ₀/2π是约瑟夫森能量。这个余弦势能，就是非线性的来源。

我的经验：在项目中选约瑟夫森结参数时，E_J和E_C（充电能量）的比值至关重要。我一般先根据目标频率估算E_J，再通过调整结面积来微调。千万别搞反了顺序，否则后面调参数会非常痛苦。

3.3 量子化：从经典到量子

好，现在我们有非线性元件了。怎么把它量子化呢？

标准做法是：把经典哈密顿量中的变量换成算符。对于电荷q和磁通φ，它们满足对易关系：

[φ, q] = iħ

这就像位置和动量的对易关系一样。于是，一个包含约瑟夫森结的电路，其量子哈密顿量可以写成：

H = 4E_C n² - E_J cos(δ)

这里n = q/2e是库珀对数算符，E_C = e²/(2C)是充电能量。注意，我用了约化单位，让事情看起来简洁一些。

你可能会问：为什么是4E_C？因为每个库珀对携带2e的电荷，所以n²前面的系数是(2e)²/(2C) = 4e²/(2C) = 4E_C。嗯，这个细节容易搞混，我刚开始学的时候也绕了半天。

3.4 能级结构：从等间距到非等间距

现在到了最精彩的部分。我们来看看这个量子化后的系统能级长什么样。

对于纯LC振荡器（E_J = 0的情况），能级是等间距的：

E_n = ħω (n + 1/2)

但加上约瑟夫森结后，余弦势能破坏了这种对称性。能级变得不等间距了。具体来说，当E_J >> E_C时，系统处于“相量子比特”区域，能级在底部近似等间距，但越往上间距越小。

我画了一张图，展示这个能级结构的变化：

<svg viewBox="0 0 600 400" xmlns="http://www.w3.org/2000/svg">
  <!-- 背景 -->
  <rect width="600" height="400" fill="#f8f9fa" rx="10"/>
  
  <!-- 标题 -->
  <text x="300" y="30" text-anchor="middle" font-size="16" font-weight="bold" fill="#333">超导量子电路能级结构对比</text>
  
  <!-- 左侧：LC振荡器（等间距） -->
  <text x="150" y="60" text-anchor="middle" font-size="14" fill="#555">LC振荡器（线性）</text>
  
  <!-- 能级线 -->
  <line x1="80" y1="100" x2="220" y2="100" stroke="#4a90d9" stroke-width="3"/>
  <text x="230" y="105" font-size="12" fill="#4a90d9">n=2</text>
  
  <line x1="80" y1="150" x2="220" y2="150" stroke="#4a90d9" stroke-width="3"/>
  <text x="230" y="155" font-size="12" fill="#4a90d9">n=1</text>
  
  <line x1="80" y1="200" x2="220" y2="200" stroke="#4a90d9" stroke-width="3"/>
  <text x="230" y="205" font-size="12" fill="#4a90d9">n=0</text>
  
  <!-- 等间距标注 -->
  <line x1="60" y1="100" x2="60" y2="200" stroke="#999" stroke-dasharray="4,4"/>
  <text x="55" y="155" text-anchor="end" font-size="11" fill="#999">ΔE=ħω</text>
  
  <!-- 右侧：约瑟夫森结电路（非等间距） -->
  <text x="450" y="60" text-anchor="middle" font-size="14" fill="#555">约瑟夫森结电路（非线性）</text>
  
  <!-- 能级线 -->
  <line x1="350" y1="90" x2="490" y2="90" stroke="#e74c3c" stroke-width="3"/>
  <text x="500" y="95" font-size="12" fill="#e74c3c">|2⟩</text>
  
  <line x1="350" y1="140" x2="490" y2="140" stroke="#e74c3c" stroke-width="3"/>
  <text x="500" y="145" font-size="12" fill="#e74c3c">|1⟩</text>
  
  <line x1="350" y1="200" x2="490" y2="200" stroke="#e74c3c" stroke-width="3"/>
  <text x="500" y="205" font-size="12" fill="#e74c3c">|0⟩</text>
  
  <!-- 非等间距标注 -->
  <line x1="340" y1="90" x2="340" y2="140" stroke="#999" stroke-dasharray="4,4"/>
  <text x="335" y="118" text-anchor="end" font-size="11" fill="#999">ΔE₁₂</text>
  
  <line x1="340" y1="140" x2="340" y2="200" stroke="#999" stroke-dasharray="4,4"/>
  <text x="335" y="173" text-anchor="end" font-size="11" fill="#999">ΔE₀₁</text>
  
  <!-- 不等间距说明 -->
  <text x="450" y="260" text-anchor="middle" font-size="12" fill="#e74c3c">ΔE₀₁ ≠ ΔE₁₂</text>
  
  <!-- 箭头 -->
  <defs>
    <marker id="arrow" markerWidth="10" markerHeight="7" refX="10" refY="3.5" orient="auto">
      <polygon points="0 0, 10 3.5, 0 7" fill="#666"/>
    </marker>
  </defs>
  
  <line x1="280" y1="150" x2="320" y2="150" stroke="#666" stroke-width="2" marker-end="url(#arrow)"/>
  
  <!-- 底部说明 -->
  <text x="300" y="310" text-anchor="middle" font-size="12" fill="#666">非线性使能级不再等间距，实现可寻址的量子比特</text>
  
  <!-- 参数区域 -->
  <rect x="150" y="330" width="300" height="50" rx="8" fill="#e8f4f8" stroke="#b8d4e8" stroke-width="1"/>
  <text x="300" y="352" text-anchor="middle" font-size="11" fill="#333">关键参数：E_J/E_C 比值决定非线性强度</text>
  <text x="300" y="370" text-anchor="middle" font-size="11" fill="#666">典型值：E_J/E_C ≈ 50~100（Transmon）</text>
</svg>

你看，左边是LC振荡器的等间距能级，右边是加了约瑟夫森结后的非等间距能级。正是这个非等间距，让我们可以用微波脉冲单独操控|0⟩→|1⟩的跃迁，而不影响|1⟩→|2⟩。

注意：非等间距不是越大越好。如果非线性太强（E_J/E_C太小），能级间距差异过大，会导致退相干时间变短。我见过有人为了追求大非线性，把E_J/E_C做到10以下，结果T1时间掉到几十纳秒，根本没法用。这是个trade-off。

3.5 实际电路中的哈密顿量

在实际的超导量子电路中，我们通常用Transmon（传输子量子比特）架构。它的哈密顿量是：

H = 4E_C n² - E_J cos(δ) + 其他耦合项

在E_J >> E_C的极限下，我们可以把余弦势能在底部做泰勒展开：

H ≈ 4E_C n² + (1/2)E_J δ² - (1/24)E_J δ⁴ + ...

前两项就是谐振子，第三项是四阶非线性修正。这个非线性项就是让能级不等间距的罪魁祸首。

我习惯把Transmon的能级近似写成：

E_m ≈ ħω_p (m + 1/2) - (E_C/12)(6m² + 6m + 3)

其中ω_p = √(8E_CE_J)/ħ是等离子体频率。你看，第二项就是非线性的贡献，它让能级向下弯曲。

参数 符号 典型值 作用
约瑟夫森能量 E_J 10~20 GHz 决定比特频率
充电能量 E_C 0.1~0.3 GHz 决定非线性强度
非谐性 α = -E_C -0.1~-0.3 GHz 能级间距差异

避坑指南：我曾经在设计一个5量子比特芯片时，忽略了相邻比特之间的串扰耦合。结果发现，当两个比特频率接近时，它们的能级会发生avoided crossing，导致操控出错。后来我学乖了，在设计阶段就用哈密顿量对角化算一遍，确保所有比特的频率间隔大于200 MHz。

好了，这一章的内容就到这里。我们从经典的LC振荡器出发，引入了约瑟夫森结这个非线性元件，然后做了量子化，最后分析了能级结构。这些是理解超导量子电路的基础，也是后续做仿真和优化的起点。

本章核心公式：H = 4E_C n² - E_J cos(δ) —— 这是超导量子电路的“标准模型”，几乎所有量子比特设计都从这里出发。

公众号:蓝海资料掘金营，微信deep3321

← 上一章 📖 返回目录 下一章 →