On the Schrödingerization method for linear non-unitary dynamics with optimal dependence on matrix queries

基本信息

• 论文ID: 2505.00370
• 标题: On the Schrödingerization method for linear non-unitary dynamics with optimal dependence on matrix queries
• 作者: Shi Jin, Nana Liu, Chuwen Ma, Yizhe Peng, Yue Yu
• 分类: math.NA cs.NA quant-ph
• 发表时间: October 14, 2025 (arXiv preprint)
• 论文链接: https://arxiv.org/abs/2505.00370

摘要

Schrödingerization方法通过所谓的扭曲相位变换将线性偏微分方程和常微分方程的非酉动力学转换为具有酉演化的Schrödinger型方程系统，该变换将原始方程映射到高一维的Schrödinger型方程。原始方法在辅助空间中使用特定的初始函数，未能实现精度的最优缩放。本文证明，通过在辅助空间中选择更光滑的初始函数，Schrödingerization方法实际上可以在矩阵查询中实现近最优甚至最优缩放。

研究背景与动机

问题背景

1. 非酉动力学的挑战：许多物理现象（如燃烧、大气和海洋环流、带物理边界的电磁波传播）表现出非酉动力学，传统的Hamiltonian模拟技术不适用。
2. 量子计算的应用需求：量子计算在处理大规模科学计算问题方面具有潜在的多项式甚至指数级计算优势，但需要酉演化算子。
3. 现有方法的局限性：
   • 原始Schrödingerization方法使用简单的初始函数ψ(p) = e^(-|p|)，由于缺乏正则性导致仅有一阶近似
   • 实现精度ε可能需要网格尺寸Δp = O(ε)，导致最大Fourier模式μ_max = O(1/ε)，这不是最优的

研究动机

通过采用更光滑的初始化函数来改善非最优的O(1/ε)缩放，从而实现对矩阵查询的最优依赖性。

核心贡献

1. 理论框架：建立了Schrödingerization方法复杂度分析的抽象框架（定理2.2）
2. 最优性条件：构造了初始辅助状态必须满足的三个必要条件(H1)-(H3)以实现最优性
3. 四种光滑初始化方法：
   • (a) 误差函数及相关函数（实现最优性）
   • (b) 截断函数（近最优）
   • (c) 高阶多项式插值（近最优）
   • (d) Fourier变换方法（近最优）
4. 最优复杂度：对于时间无关情况，实现了Õ(αₕT log(1/ε))的矩阵查询复杂度，达到最优依赖性

方法详解

任务定义

考虑线性动力学系统：

du/dt = A(t)u(t) + b(t), t ∈ (0,T)
u(0) = u₀

其中A一般不是反Hermitian矩阵，目标是在量子计算机上高效求解此系统。

Schrödingerization方法架构

1. 齐次化

通过引入辅助向量r(t)将非齐次系统转换为齐次系统：

d/dt u_f = A_f u_f, A_f = [A B; O O], u_f(0) = [u₀; r₀]

2. 扭曲相位变换

使用变换w(t,p) = e^(-p)u_f(t)对p ≥ 0，并对称扩展到p < 0：

∂w/∂t = -H₁∂_p w + iH₂w
w(0,p) = ψ(p)u_I

其中H₁ = (A_f + A_f†)/2, H₂ = (A_f - A_f†)/(2i)

3. 离散Fourier变换

在p方向应用离散Fourier变换得到：

d/dt W_h(t) = -i(P_μ ⊗ H₁)W_h + i(I ⊗ H₂)W_h

技术创新点

1. 光滑初始化函数设计

误差函数方法（最优）：

ψ(p) = φ(p)e^(-p), φ(p) = (erf(ap) + 1)/2

其中a = 2log^(1/2)(1/ε)，实现了‖ψ^(r)‖^(1/r)_(L²) ≤ Cr的最优界。

截断函数方法（近最优）： 使用mollifier与阶跃函数的卷积构造光滑扩展，但由于mollifier不是解析的，只能达到β = 1/2。

2. 复杂度分析框架

建立了抽象定理2.2，将查询复杂度与初始函数的正则性联系：

• 如果‖ψ^(r)‖^(1/r)_(L²) ≤ Cr^(1/β)，则μ_max ≲ (log(1/ε))^(1/β)
• 最优性要求β = 1

3. 最优性的三个必要条件

(H1) 指数衰减：ψ(p)在R上表现出指数衰减 (H2) 近似性：对p ∈ p*, R，|ψ(p) - e^(-p)| ≤ ε (H3) 正则性：‖ψ^(r)‖^(1/r)_(L²) ≤ Cr当r ≃ log(1/ε)

实验设置

理论分析方法

本文主要进行理论分析，通过数学证明验证：

1. 误差估计（定理2.1, 5.1-5.3）
2. 复杂度分析（定理2.2, 4.1）
3. 光滑初始化的构造和性质分析

对比方法

与现有量子ODE求解方法比较（表1）：

• 谱方法20
• 截断Dyson级数13
• 时间步进方法16
• 改进的LCHS方法27
• 最优LCHS方法43

实验结果

主要结果

复杂度比较（表1）：

• 本文方法（时间无关）：Õ(u_r α_A T log(1/ε))
• 改进LCHS（时间无关）：Õ(u_r α_A T (log(1/ε))^(1/β)), β < 1
• 最优LCHS：Õ(u_r α_A T log(1/ε))

本文方法在时间无关情况下达到与最优LCHS相同的复杂度。

理论发现

1. 正则性与复杂度的关系

通过Parseval恒等式：‖ψ^(r)‖(L²) = ‖w^r ψ̂‖(L²)

• C^∞但非解析函数：|ψ̂(w)| ≤ Ce^(-c|w|^β), β < 1 → 次最优
• 解析函数：|ψ̂(w)| ≤ Ce^(-c|w|) → 最优

2. 误差函数的优越性

误差函数erf(p)及其构造的ψ(p)满足Gevrey-1类条件，实现了所需的最优界：

‖ψ^(r)‖^(1/r)_(L²) ≲ ar^(1/2) ≃ log^(1/2)(1/ε) · r^(1/2) ≃ r

相关工作

主要研究方向

1. 量子线性系统算法：QLSA及其改进
2. Hamiltonian模拟：基于块编码的技术
3. 系统酉化：将非酉系统扩展为酉系统
4. LCHS方法：线性Hamiltonian模拟的线性组合

本文优势

• 相比原始Schrödingerization：从O(1/ε)改进到O(log(1/ε))
• 相比改进LCHS：从O((log(1/ε))^(1/β))改进到O(log(1/ε))
• 提供了实现最优性的具体构造方法

结论与讨论

主要结论

1. 最优性可达成：通过选择合适的光滑初始化函数，Schrödingerization方法可以实现对矩阵查询的最优依赖
2. 误差函数最优：基于误差函数的初始化实现了β = 1的最优缩放
3. 理论完备性：提供了完整的复杂度分析框架和误差估计

局限性

1. 解析性要求：实现最优性需要初始函数具有解析性，限制了函数选择
2. 时间依赖情况：时间依赖系统需要额外的对数因子O((log(1/ε))²)
3. 实际实现：理论结果需要在实际量子设备上验证

未来方向

1. 在NISQ设备上的实际实现和验证
2. 扩展到非线性微分方程
3. 优化时间依赖情况的复杂度

深度评价

优点

1. 理论严谨：提供了完整的数学证明和复杂度分析
2. 方法创新：首次实现Schrödingerization方法的最优复杂度
3. 实用价值：为量子PDE/ODE求解提供了理论最优的算法
4. 分析深入：深刻揭示了初始函数正则性与算法复杂度的关系

不足

1. 缺乏数值验证：主要是理论工作，缺乏数值实验验证
2. 构造复杂性：最优初始函数的构造相对复杂
3. 适用范围：仍限制于线性系统和特定的矩阵性质假设

影响力

1. 理论贡献：为量子微分方程求解领域提供了重要的理论突破
2. 方法指导：为设计高效量子算法提供了新的思路
3. 实用潜力：在量子计算成熟后具有重要应用价值

适用场景

• 大规模线性微分方程组求解
• 具有非酉动力学的物理系统模拟
• 需要高精度的量子科学计算应用

参考文献

本文引用了52篇相关文献，涵盖量子计算、数值分析、偏微分方程等多个领域的重要工作，为研究提供了坚实的理论基础。