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

Basic Information

• Paper ID: 2505.00370
• Title: On the Schrödingerization method for linear non-unitary dynamics with optimal dependence on matrix queries
• Authors: Shi Jin, Nana Liu, Chuwen Ma, Yizhe Peng, Yue Yu
• Classification: math.NA cs.NA quant-ph
• Publication Date: October 14, 2025 (arXiv preprint)
• Paper Link: https://arxiv.org/abs/2505.00370

Abstract

The Schrödingerization method converts non-unitary dynamics of linear partial and ordinary differential equations into Schrödinger-type equation systems with unitary evolution through a so-called twisted phase transformation, which maps the original equations to higher-dimensional Schrödinger-type equations. The original method employs a specific initial function in the auxiliary space, failing to achieve optimal scaling of accuracy. This paper demonstrates that by selecting smoother initial functions in the auxiliary space, the Schrödingerization method can indeed achieve near-optimal or even optimal scaling in matrix queries.

Research Background and Motivation

Problem Background

1. Challenges of Non-unitary Dynamics: Many physical phenomena (such as combustion, atmospheric and oceanic circulation, electromagnetic wave propagation with physical boundaries) exhibit non-unitary dynamics, for which traditional Hamiltonian simulation techniques are inapplicable.
2. Application Requirements in Quantum Computing: Quantum computing possesses potential polynomial or even exponential computational advantages for large-scale scientific computing problems, but requires unitary evolution operators.
3. Limitations of Existing Methods:
   • The original Schrödingerization method uses a simple initial function ψ(p) = e^(-|p|), achieving only first-order approximation due to lack of regularity
   • Achieving precision ε may require grid size Δp = O(ε), leading to maximum Fourier mode μ_max = O(1/ε), which is suboptimal

Research Motivation

To improve the suboptimal O(1/ε) scaling by adopting smoother initialization functions, thereby achieving optimal dependence on matrix queries.

Core Contributions

1. Theoretical Framework: Establishes an abstract framework for complexity analysis of the Schrödingerization method (Theorem 2.2)
2. Optimality Conditions: Constructs three necessary conditions (H1)-(H3) that the initial auxiliary state must satisfy to achieve optimality
3. Four Smooth Initialization Methods:
   • (a) Error function and related functions (achieves optimality)
   • (b) Truncation functions (near-optimal)
   • (c) High-order polynomial interpolation (near-optimal)
   • (d) Fourier transform method (near-optimal)
4. Optimal Complexity: For time-independent cases, achieves Õ(α_H T log(1/ε)) matrix query complexity, attaining optimal dependence

Methodology Details

Task Definition

Consider the linear dynamical system:

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

where A is generally not an anti-Hermitian matrix, with the goal of efficiently solving this system on a quantum computer.

Schrödingerization Method Architecture

1. Homogenization

Convert the non-homogeneous system to a homogeneous system by introducing auxiliary vector r(t):

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

2. Twisted Phase Transformation

Use the transformation w(t,p) = e^(-p)u_f(t) for p ≥ 0, with symmetric extension to p < 0:

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

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

3. Discrete Fourier Transform

Apply discrete Fourier transform in the p direction to obtain:

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

Technical Innovations

1. Smooth Initialization Function Design

Error Function Method (Optimal):

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

where a = 2log^(1/2)(1/ε), achieving optimal bounds ‖ψ^(r)‖^(1/r)_(L²) ≤ Cr.

Truncation Function Method (Near-optimal): Construct smooth extensions using convolution of mollifiers with step functions, but achieving only β = 1/2 due to mollifiers not being analytic.

2. Complexity Analysis Framework

Establishes abstract Theorem 2.2, linking query complexity to regularity of the initial function:

• If ‖ψ^(r)‖^(1/r)_(L²) ≤ Cr^(1/β), then μ_max ≲ (log(1/ε))^(1/β)
• Optimality requires β = 1

3. Three Necessary Conditions for Optimality

(H1) Exponential Decay: ψ(p) exhibits exponential decay on ℝ (H2) Approximation: For p ∈ p*, R, |ψ(p) - e^(-p)| ≤ ε (H3) Regularity: ‖ψ^(r)‖^(1/r)_(L²) ≤ Cr when r ≃ log(1/ε)

Experimental Setup

Theoretical Analysis Method

This paper primarily conducts theoretical analysis through mathematical proofs to verify:

1. Error estimates (Theorems 2.1, 5.1-5.3)
2. Complexity analysis (Theorems 2.2, 4.1)
3. Construction and property analysis of smooth initializations

Comparison Methods

Comparison with existing quantum ODE solvers (Table 1):

• Spectral methods 20
• Truncated Dyson series 13
• Time-stepping methods 16
• Improved LCHS method 27
• Optimal LCHS method 43

Experimental Results

Main Results

Complexity Comparison (Table 1):

• This paper's method (time-independent): Õ(u_r α_A T log(1/ε))
• Improved LCHS (time-independent): Õ(u_r α_A T (log(1/ε))^(1/β)), β < 1
• Optimal LCHS: Õ(u_r α_A T log(1/ε))

This paper's method achieves the same complexity as optimal LCHS in time-independent cases.

Theoretical Findings

1. Relationship Between Regularity and Complexity

Via Parseval's identity: ‖ψ^(r)‖(L²) = ‖w^r ψ̂‖(L²)

• C^∞ but non-analytic functions: |ψ̂(w)| ≤ Ce^(-c|w|^β), β < 1 → suboptimal
• Analytic functions: |ψ̂(w)| ≤ Ce^(-c|w|) → optimal

2. Superiority of Error Function

The error function erf(p) and its constructed ψ(p) satisfy Gevrey-1 class conditions, achieving the required optimal bounds:

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

Related Work

Main Research Directions

1. Quantum Linear System Algorithms: QLSA and its improvements
2. Hamiltonian Simulation: Block-encoding-based techniques
3. System Unitarization: Extending non-unitary systems to unitary systems
4. LCHS Method: Linear combinations of Hamiltonian simulations

Advantages of This Paper

• Compared to original Schrödingerization: improvement from O(1/ε) to O(log(1/ε))
• Compared to improved LCHS: improvement from O((log(1/ε))^(1/β)) to O(log(1/ε))
• Provides concrete construction methods for achieving optimality

Conclusions and Discussion

Main Conclusions

1. Optimality is Achievable: By selecting appropriate smooth initialization functions, the Schrödingerization method can achieve optimal dependence on matrix queries
2. Error Function is Optimal: Initialization based on error functions achieves optimal scaling with β = 1
3. Theoretical Completeness: Provides a complete complexity analysis framework and error estimates

Limitations

1. Analyticity Requirement: Achieving optimality requires the initial function to be analytic, restricting function choices
2. Time-Dependent Cases: Time-dependent systems require an additional logarithmic factor O((log(1/ε))²)
3. Practical Implementation: Theoretical results require verification on actual quantum devices

Future Directions

1. Practical implementation and verification on NISQ devices
2. Extension to nonlinear differential equations
3. Optimization of complexity for time-dependent cases

In-Depth Evaluation

Strengths

1. Theoretical Rigor: Provides complete mathematical proofs and complexity analysis
2. Methodological Innovation: First to achieve optimal complexity for the Schrödingerization method
3. Practical Value: Provides theoretically optimal algorithms for quantum PDE/ODE solving
4. Deep Analysis: Profoundly reveals the relationship between initial function regularity and algorithm complexity

Weaknesses

1. Lack of Numerical Verification: Primarily theoretical work lacking numerical experimental validation
2. Construction Complexity: Construction of optimal initial functions is relatively complex
3. Limited Scope: Still restricted to linear systems and specific matrix property assumptions

Impact

1. Theoretical Contribution: Provides important theoretical breakthrough in quantum differential equation solving
2. Methodological Guidance: Offers new insights for designing efficient quantum algorithms
3. Practical Potential: Possesses significant application value once quantum computing matures

Applicable Scenarios

• Solving large-scale linear differential equation systems
• Simulating physical systems with non-unitary dynamics
• Quantum scientific computing applications requiring high precision

References

This paper cites 52 relevant references, covering important works in quantum computing, numerical analysis, partial differential equations, and other fields, providing a solid theoretical foundation for the research.