1. 单光束拉曼跃迁的物理原理与实现在量子计算中单比特门操作通常通过拉曼跃迁实现。拉曼跃迁是一种双光子过程通过虚拟中间态耦合两个量子态。对于173Yb原子的1S0-3P1跃迁单光束拉曼跃迁技术的核心在于利用激光诱导的光频移Light Shift来实现量子态的相干操控。1.1 光频移哈密顿量的数学表达光频移哈密顿量可以表示为ˆHLS 3πc²Γ/(2ω₀³) × Σ[⟨Fmb|e∗L·ˆD(q)FF′|F′me⟩⟨F′me|eL·ˆD†(q)F′F|Fma⟩/(ΔL−ΔHFS(F′))] × I(q)L × |Fmb⟩⟨Fma|其中关键参数包括c光速ω₀1S0-3P1跃迁的共振频率Γ3P1态的自然线宽ΔL控制激光频率失谐ΔHFS(F′)超精细能级F′的共振频率这个哈密顿量可以通过三个参数进行调控激光强度IL激光频率失谐ΔL激光偏振eL1.2 轨道湮灭算符的展开轨道湮灭算符ˆD(q)FF′可以进一步展开为ˆD(q)FF′ eq(-1)F′J1Ip√[(2F′1)(2J1)] × Σ[⟨Fmg|F′me;1−q⟩{...}|Fmg⟩⟨F′me|其中{...}表示Wigner 6j符号⟨Fmg|F′me;1−q⟩是Clebsch-Gordan系数可以用Wigner 3j符号表示。2. 自旋动力学与量子主方程2.1 量子主方程的构建单光束拉曼跃迁诱导的多自旋动力学可以用量子主方程描述∂ˆρ(t)/∂t -i/ℏ[ˆHrot, ˆρ(t)] Σ[2ˆCnˆρ(t)ˆC†n - ˆρ(t)ˆC†nˆCn - ˆC†nˆCnˆρ(t)]/2其中旋转框架下的哈密顿量为ˆHrot ˆd(F)(π/2) ˆHLS ˆd(F)†(π/2)2.2 坍塌算符的表达式来自控制激光光子散射的坍塌算符定义为ˆCq √Γ × Σ[Ω/2/(ΔL−ΔHFS(F′)iΓ/2)] × (e∗q·ˆD(q)FF′)(eL·ˆD†(q)F′F)在实际模拟中我们使用固定的激光失谐协变SU(2)旋转11.217 GHz相对于1S0-3P1(F′7/2)跃迁非线性旋转-5.005 GHz 激光采用σ圆偏振并利用实验获得的初始自旋布居作为ˆρ(0)。3. 最优微分光频移条件3.1 微分光频移的必要条件为了实现ˆRx(π)和ˆR(cat)x(π/2)门操作微分光频移(DLS)需要满足特定比例关系对于ˆRx(π)门Δ1:Δ2:Δ3:Δ4:Δ5 (2n11):(2n21):(2n31):(2n41):(2n51)对于ˆR(cat)x(π/2)门Δ1:Δ2:Δ3:Δ4:Δ5 (4n1±1):(4n2∓1):(4n3±1):(4n4∓1):(4n5±1)这些条件保证了在特定时间t幺正时间演化算符ˆU(F)rot(t)能够精确实现目标门操作。3.2 充分性证明通过Wigner小d矩阵的显式表达式我们可以证明当DLS满足上述比例关系时对于ˆRx(π)门在t(2n11)π/Δ1时刻实现对于ˆR(cat)x(π/2)门在t(4n1∓1)π/(2Δ1)时刻实现虽然这些最优DLS条件是在F5/2系统中推导的但它们同样适用于其他自旋系统如87Sr的1S0基态F9/2。4. 最优激光失谐的数值确定4.1 数值模拟方法我们通过量子主方程模拟计算了173Yb原子1S0-3P1跃迁中不同激光失谐下的门保真度。门保真度定义为F(ˆρ1,ˆρ2) : |Tr[√√ˆρ1ˆρ2√ˆρ1]|²4.2 关键实验结果图10(a,b)和表I、II展示了在特定失谐下门错误率显著降低的区域。这些结果验证了推导的最优DLS条件在大自旋系统中实现高保真度单比特门的有效性。特别值得注意的是对于ˆRx(π)门在11.217 GHz失谐处获得最低错误率3.59×10⁻⁵对于ˆR(cat)x(π/2)门在-5.005 GHz失谐处获得错误率5.15×10⁻⁴5. 对DLS比例偏差的鲁棒性5.1 偏差表征定义DLS比例偏差的标准差 对于ˆRx(π)门σX √[1/(2F) Σ(Δk/[Δ1/(2n11)] - (2nk1))²]对于ˆR(cat)x(π/2)门σH √[1/(2F) Σ(Δk/[Δ1/(4n1±1)] - (4nk∓(-1)k))²]5.2 鲁棒性分析如图10(c,d)所示当标准偏差σX和σH小于约0.01时门错误率可以保持在10⁻³以下。这意味着即使激光失谐偏离最优值几百MHz仍能实现高保真度门操作。值得注意的是光子散射效应也会显著影响门错误率。例如在-11.9 GHz和11.2 GHz失谐处虽然σX相当但由于散射事件数量的差异分别为6.4×10⁻²和1.6×10⁻³门错误率相差约一个数量级。6. 单比特门的秩保持条件6.1 容错量子计算的要求在容错量子计算中单组件故障最多只能在每个编码块中产生一个错误。对于自旋猫编码某些理想的单比特门可能会增加跳跃错误数量将可纠正错误转变为不可纠正错误。6.2 秩保持条件秩保持条件要求对于每个操作ˆE∈EK满足ˆUˆEˆU† ⊂ EK其中EK是可纠正错误的集合。自旋F的SU(2)旋转协变SU(2)旋转满足这一秩保持条件。其旋转矩阵D(k)q′q(α,β,γ)就是著名的Wigner D矩阵。6.3 逻辑门实现Pauli门可以通过自旋F的SU(2)旋转以秩保持方式实现ˆX ˆD(π,π,0) ˆY ˆD(0,π,0) ˆZ ˆD(π,0,0)然而Hadamard门等不能仅用SU(2)旋转表示的操作是否满足秩保持条件仍需进一步研究。图11的模拟结果表明只有当F1/2时Hadamard门才能用单一旋转算符表示。7. 噪声偏置二面体随机基准测试7.1 DRB协议步骤制备ˆZ或ˆX算符的本征态通常为|0⟩或|⟩应用从单比特Pauli群中随机采样的酉操作ˆP应用从D8二面体群中随机采样的m个门序列应用将量子态返回初始态的反向操作在z或x基下测量最终态对不同长度m2重复上述步骤对结果取平均得到每个长度m2的存活概率7.2 数据分析模型z基和x基DRB的平均返回概率遵循指数衰减模型P1(m) A1λ1^(m2) B1 P2(m) A2λ2^(m4) B2通过拟合可以得到退相干错误概率pD和非退相干错误概率pND进而计算偏置参数η pD/pND。8. 误差预算分析我们考虑了多种误差源对单比特门的影响包括激光强度波动脉冲面积波动的标准差σ为平均脉冲面积的0.2%QB1和0.56%QB2激光偏振波动方位角φ和椭圆率χ的标准差分别为3.1°(QB1)/0.8°(QB2)和0.13°(QB1)/0.02°(QB2)激光频率波动线宽标准差σ小于1MHz正交性缺陷QB1和QB2的交叉角与90°的偏差为2°有限塞曼分裂CRB测量使用1.01mT磁场DRB测量使用1.35mT退相干使用实验测量的退相干时间T*2,k 251(21)ms/(F-k) (k0,1,2)光子散射通过坍塌算符ˆCq建模通过技术改进如稳定激光强度和偏振、改善正交性对准、实现更快的门操作以减少塞曼分裂影响我们预期可以实现超过0.999的门保真度。