1. 量子计算赋能动态平均场理论的技术突破在凝聚态物理和计算化学领域强关联电子体系的理论描述一直是个巨大挑战。动态平均场理论(DMFT)作为处理这类问题的有效方法通过将复杂晶格模型映射到更简单的杂质模型来预测系统行为。传统DMFT计算的核心瓶颈在于杂质格林函数的求解——量子蒙特卡洛(QMC)面临符号问题精确对角化(ED)遭遇希尔伯特空间指数增长即使张量网络等现代方法也只能部分缓解这些限制。我们的研究团队开发了一套基于量子计算的创新解决方案其技术内核包含三个关键突破采用高斯子空间(SGS)表示基态将基态准备复杂度从指数级降为多项式级设计部分线路压缩算法大幅减少时间演化所需的量子门数量开发正定矩阵降噪技术有效提升含噪量子硬件的信号质量这套混合量子-经典框架已在IBM量子处理器上成功实现8量子比特(含1辅助比特)的杂质格林函数计算为材料科学中的关联体系模拟开辟了新途径。关键技术指标对于单杂质耦合3个浴轨道(共8量子比特)的系统我们的方法将基态表示所需资源从传统ED的4900维希尔伯特空间压缩至仅需12个高斯态同时通过线路压缩将双量子比特门数量减少60%以上。2. 高斯子空间表示的理论基础与实现2.1 费米高斯态(FGS)的数学特性费米高斯态是描述非相互作用费米子系统的精确基态其独特性质使其成为理想的基函数协方差矩阵表示每个N轨道系统的FGS可由2N×2N协方差矩阵完全描述计算高效性非正交FGS间的期望值和重叠可通过计算最多4×4子矩阵的Pfaffian获得量子线路友好任何FGS都可用多项式深度的自由费米子演化线路精确制备我们构建的SGS方法利用这些特性通过线性组合多个FGS来近似相互作用系统的基态。具体实现时每个FGS通过调节杂质问题的非相互作用参数如跃迁强度和在位能生成而库仑排斥参数U设为零。2.2 子空间选择算法子空间构建过程采用增量式优化策略候选池生成预生成包含1000个FGS的候选池M通过均匀变化非相互作用参数获得基态逼近首个FGS |ϕ₁⟩选择使⟨ϕ₁|Ĥimp|ϕ₁⟩最小的状态正交性优化后续状态选择与现有子空间保持最大正交性的FGS收敛判断当添加新FGS引起的能量变化可忽略或候选池耗尽时停止这一过程的关键创新在于利用特征向量延拓技术(Eigenvector Continuation)——当杂质哈密顿量参数缓慢变化时基态会保持在低能子空间内使得同一SGS可重复用于不同问题实例。# 伪代码子空间选择算法 def build_subspace(H_imp, candidate_pool): subspace [] energy_history [] while not converged(energy_history): next_state select_most_orthogonal(candidate_pool, subspace) subspace.append(next_state) # 子空间对角化 H_sub construct_subspace_matrix(H_imp, subspace) E_gs, alpha diagonalize(H_sub) energy_history.append(E_gs) if len(subspace) max_size or is_ill_conditioned(H_sub): break return subspace, alpha2.3 精度验证与性能表现我们通过两个指标验证SGS的准确性基态能量相对误差EGS |Ẽ(0)-E(0)|/|E(0)|杂质格林函数特征再现能力测试结果显示对于NI1、NB7的系统希尔伯特空间维度4900仅需12个FGS即可将EGS控制在10⁻²以内。更令人振奋的是即使当NB7时希尔伯特空间维度远超SGS秩GRimp(ω)仍能准确再现ED结果的所有特征峰包括那些强度较小的高频峰图3b。表1展示了不同系统规模下的SGS性能数据(NI,NB)希尔伯特空间维度平均SGS秩占比(1,3)3611.130.8%(1,7)49006.70.14%(2,3)49009.40.19%数据揭示了一个重要规律随着系统规模增大所需SGS秩与希尔伯特空间维度的比值迅速下降展现出优异的可扩展性。3. 量子线路设计与压缩技术3.1 杂质格林函数的量子测量方案计算杂质格林函数的核心是评估以下关联函数GRimp(t) Σᵢⱼ(αᵢ(0))*αⱼ(0)Cᵢⱼ(t)其中Cᵢⱼ(t) ⟨ϕᵢ|{d̂(t),d̂†}|ϕⱼ⟩。由于湮灭和产生算子非酉我们引入马约拉纳算子γ̂₊d̂†d̂和γ̂₋i(d̂-d̂†)进行转换Cᵢⱼ(t) ¼Σₐₑsₐₑ⟨ϕᵢ|{γ̂ₐ(t),γ̂ₑ}|ϕⱼ⟩每个项可通过Hadamard测试测量其量子线路如图5所示。传统方法需要深度随系统规模快速增长的线路我们的创新在于开发了部分压缩技术。3.2 线路压缩的四阶段优化控制结构简化利用Hadamard测试的测量特性移除冗余控制门将状态准备移至线路起始端消除一个时间演化算子图5bTrotter步压缩将每个Trotter步压缩为紧凑的匹配门结构图5c状态准备融合将两个受控自由费米演化酉算子Uᵢ和Uⱼ合并为单个短深度酉算子图5dTrotter吸收将初始Trotter步的一半吸收到状态准备中图5e对于NI1、NB3的系统8量子比特1辅助比特经过压缩后的线路如图6所示。实测显示这种策略可减少约60%的双量子比特门数量大幅降低噪声影响。操作提示在实现匹配门时建议使用原生量子门集进行分解。例如IBM量子处理器上的CNOT单量子比特门组合需注意不同量子芯片的拓扑结构对门分解效率的影响。3.3 硬件实现与误差缓解我们在IBM Sherbrooke处理器上执行了绝缘相(U5.337)的关联函数测量采用多层误差缓解技术泡利随机编译进行10轮泡利扭曲(Pauli Twirling)平均化系统误差动态解耦采用XY脉冲序列抑制退相干效应基于门数的重标定根据门数量和错误率调整测量结果对于χ4的SGS基和nₜ20个时间点共需要χ²nₜ320次电路执行。图7展示了原始硬件数据、误差缓解后数据以及经正定降噪处理后的结果对比。4. 数据处理与DMFT集成4.1 正定矩阵降噪技术量子硬件输出的GRimp(t)数据通常包含噪声导致对应的Gram矩阵失去正定性。我们采用以下处理流程构造Gram矩阵将离散时间点的GRimp(t)排列成Toeplitz矩阵正定投影找到与噪声Gram矩阵最接近的正定矩阵正定延拓利用正定函数的性质将有限时间数据延拓到更长时域图7c展示了处理效果——延拓后的数据不仅能锐化主要特征峰还能揭示大峰肩部(ω≈3.5)的次级峰结构尽管在ω5区域仍存在噪声引起的伪峰。4.2 Matsubara频率转换DMFT通常在虚频轴上进行自洽计算我们从时域数据提取Matsubara格林函数的步骤傅里叶分析对GRimp(t)做离散傅里叶变换提取频率fᵣ和振幅Aᵣ解析延拓利用正定延拓获得更长时域动态提高频率分辨率极点展开按Gimp(iωₙ)ΣᵣAᵣ/(iωₙ-fᵣ)构造Matsubara格林函数这种方法避免了直接测量虚时关联函数的困难充分利用了量子计算机在实时间演化中的优势。4.3 DMFT自洽验证为验证SGS在完整DMFT流程中的可靠性我们在无噪声环境下测试了Bethe晶格上Hubbard模型的DMFT收敛初始化构建包含24个FGS的SGS(U4.449时)固定子空间在整个DMFT循环中保持SGS不变比较基准与ED结果对比晶格态密度(DOS)和准粒子权重Z图4显示SGS与ED获得的晶格DOS高度一致金属-绝缘体相变的序参数Z也呈现相同行为证实了方法的可靠性。5. 技术局限与发展前景5.1 当前方法的局限性尽管取得了显著进展我们的框架仍存在一些限制候选池效率生成大量FGS可能导致冗余和病态问题需要正则化处理初始参数敏感性浴轨道参数的初始猜测质量显著影响子空间有效性硬件噪声缺乏量子纠错下大系统模拟仍具挑战性5.2 未来改进方向自适应子空间优化结合代理模型优化初始浴参数减少人工干预团簇扩展将方法推广到团簇DMFT处理更复杂的多杂质系统含时哈密顿量探索显式时间依赖的相互作用浴构建方法错误抑制技术集成最近发展的错误缓解方案如零噪声外推在实际操作中我们建议研究人员先在小规模系统上验证SGS的适用性再逐步扩大规模。对于NI≥2的系统要特别注意子空间秩随杂质数量的增长关系。这项工作的真正价值在于为量子计算在凝聚态物理中的应用提供了可行路径。通过精心设计的经典-量子分工我们既利用了经典计算机在基态准备方面的效率又发挥了量子处理器在时间演化中的独特优势。随着量子硬件的持续进步这种方法有望在描述真实材料系统方面展现出量子优势。