pith. sign in

arxiv: 2604.14079 · v1 · submitted 2026-04-15 · 🪐 quant-ph

Hybrid quantum-classical algorithms for complex nonlinear partial differential equations with Ginzburg-Landau potential and vortex motion laws

Shi Jin , Nana Liu , Chuwen Ma This is my paper

Pith reviewed 2026-05-10 13:49 UTC · model grok-4.3

classification 🪐 quant-ph
keywords hybrid quantum-classical algorithmsnonlinear Schrödinger equationGinzburg-Landau potentialvortex dynamicsquantum algorithms for PDEsexponential speedupSchrödingerization
0
0 comments X

The pith

Hybrid quantum-classical algorithms solve nonlinear PDEs by classically advancing vortices and quantumly solving coupled elliptic equations, with exponential spatial scaling gains in 2D.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops hybrid algorithms for complex nonlinear partial differential equations such as the nonlinear Schrödinger equation in strongly nonlinear regimes dominated by vortex cores and phase singularities. It reduces the problem asymptotically to linear elliptic equations paired with low-dimensional vortex dynamics, advancing the vortices with classical methods while handling the elliptic boundary-value problems via quantum techniques including BPX preconditioning and Schrödingerization. For the two-dimensional case this produces an exponential improvement in scaling with spatial problem size while accuracy dependence stays nearly linear up to polylog factors. Readers may care because these equations describe superconductors, superfluids, and nonlinear optics where large-scale direct simulation remains difficult.

Core claim

In the strongly nonlinear regime the solutions to these PDEs are asymptotically governed in leading order by linear elliptic equations coupled with low-dimensional vortex dynamics. The hybrid algorithm advances the vortex motion classically and solves the linear elliptic problems with quantum algorithms. For the two-dimensional nonlinear Schrödinger equation, combining quantum BPX preconditioning with Schrödingerization to estimate observables in the small-output regime yields an exponential improvement in the dependence on spatial problem size, while the dependence on target accuracy remains essentially linear up to polylogarithmic factors.

What carries the argument

Asymptotic reduction of the nonlinear PDE to a linear elliptic boundary-value problem coupled with low-dimensional vortex dynamics, enabling classical vortex advancement and quantum solution of the elliptic system.

If this is right

  • The same hybrid principle extends to dissipative Ginzburg-Landau vortex dynamics.
  • The approach applies to vortex filaments in three-dimensional superconductivity.
  • Physically relevant observables can be estimated efficiently in the small-output regime with the stated scaling.
  • Numerical results confirm both the validity of the PDE reduction and the effectiveness of the hybrid method.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Identifying low-dimensional classical structures such as vortices could serve as a general route to quantum advantage for other nonlinear PDEs in physics and engineering.
  • Further combination with improved quantum linear solvers might reduce the remaining linear dependence on target accuracy.
  • The method could enable simulation of larger-scale vortex systems in superconductivity models than feasible with purely classical techniques.

Load-bearing premise

The solutions to these complex nonlinear PDEs in the strongly nonlinear regime are asymptotically governed in leading order by linear elliptic equations coupled with low-dimensional vortex dynamics.

What would settle it

A numerical or theoretical demonstration that the hybrid algorithm's runtime for the two-dimensional nonlinear Schrödinger equation scales only polynomially with spatial problem size rather than exponentially would disprove the claimed improvement.

Figures

Figures reproduced from arXiv: 2604.14079 by Chuwen Ma, Nana Liu, Shi Jin.

Figure 1
Figure 1. Figure 1: Vortex trajectories generated by M1 and M2 for ε = 0.025 [PITH_FULL_IMAGE:figures/full_fig_p025_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: shows the common masked region and the phase comparison at the final time. In the left panel, the background colors represent the modulus |u ε |: bright colors correspond to values close to 1, while dark blue regions indicate small values of |u ε |. These dark blue regions correspond to the vortex cores and their nearby neighborhoods, since |u ε | ≈ 0 near a vortex center. The dashed circles mark the exclu… view at source ↗
Figure 3
Figure 3. Figure 3: Final-time masked errors Eε M1 (T) and Eε M2 (T) versus ε. Finally, [PITH_FULL_IMAGE:figures/full_fig_p026_3.png] view at source ↗
read the original abstract

We propose quantum algorithms for complex-valued nonlinear partial differential equations in the strongly nonlinear regime, where the dynamics is governed by vortex cores, phase singularities, and nonlinear vortex interactions. Examples include the complex-valued nonlinear Schr\"odinger equation, as well as nonlinear heat and wave equations with Ginzburg--Landau-type nonlinearity. In the strongly nonlinear regime, the solutions to these equations are asymptotically governed by, in leading order, linear elliptic equations, coupled with low-dimensional vortex dynamics, where the vortex cores correspond to topological defects in superconductors. Our hybrid quantum-classical algorithms utilize this asymptotic property, in which the vortex dynamic is advanced classically while the boundary-value problem of linear elliptic equation is handled by quantum algorithms. For the two-dimensional nonlinear Schr\"odinger equation, we also combine quantum BPX preconditioning with Schr\"odingerization to estimate physically relevant observables in the small-output regime. This yields, already in two dimensions, an {\it exponential} improvement in the dependence on the spatial problem size, while the dependence on the target accuracy remains essentially linear up to polylogarithmic factors. We further show that the same principle extends to dissipative Ginzburg--Landau vortex dynamics and to vortex filaments in three-dimensional superconductivity. Numerical results support the validity of this PDE reduction and the effectiveness of the proposed approach.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper proposes hybrid quantum-classical algorithms for complex nonlinear PDEs (nonlinear Schrödinger, heat, and wave equations with Ginzburg-Landau nonlinearity) in the strongly nonlinear regime. It exploits an asymptotic reduction in which solutions are governed at leading order by linear elliptic equations coupled to low-dimensional vortex dynamics; vortex motion is advanced classically while the elliptic boundary-value problems are solved with quantum linear-system algorithms. For the 2D nonlinear Schrödinger equation the authors combine quantum BPX preconditioning with Schrödingerization to estimate observables, claiming an exponential improvement in the dependence on spatial problem size while the dependence on target accuracy remains essentially linear (up to polylog factors). The same principle is asserted to extend to dissipative Ginzburg-Landau dynamics and 3D vortex filaments; numerical results are presented in support.

Significance. If the asymptotic reduction is accompanied by explicit, ε-dependent error bounds that guarantee the reduction error remains o(ε) uniformly in the regime of interest, the work would constitute a meaningful advance in quantum algorithms for nonlinear PDEs arising in superconductivity and fluid dynamics. The hybrid separation of low-dimensional classical vortex tracking from quantum-solvable linear elliptic subproblems is conceptually clean, and the incorporation of BPX preconditioning together with Schrödingerization for small-output observables is a technical contribution that could be reusable.

major comments (2)
  1. [§2 (Asymptotic Reduction)] §2 (Asymptotic Reduction) and abstract: The central claim that solutions are 'asymptotically governed, in leading order, by linear elliptic equations, coupled with low-dimensional vortex dynamics' is asserted without a quantitative error estimate. No lemma or theorem supplies a concrete rate (e.g., O(δ^k) with δ the core size or inter-vortex distance) that is guaranteed to be smaller than the target accuracy ε. Because the quantum linear solver is applied to the reduced model, the absence of such a bound means the stated complexity for the original nonlinear PDE is not yet justified.
  2. [§4 (Complexity Analysis)] §4 (Complexity Analysis): The exponential improvement in spatial-size dependence for the 2D NLS is derived for the hybrid reduced system. Without an explicit relation showing that the reduction error is o(ε) uniformly, the overall runtime bound claimed for the original nonlinear problem does not follow.
minor comments (2)
  1. [Abstract] Abstract: the phrase 'essentially linear up to polylogarithmic factors' should be replaced by the precise scaling (including the polylog degree) once the analysis is complete.
  2. [Numerical Results] Numerical section: more quantitative tables comparing the reduction error to the target ε, together with details on how the quantum linear solver is simulated (gate counts, condition-number estimates), would make the supporting evidence easier to assess.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive feedback. The points raised highlight the need for explicit error control in the asymptotic reduction to rigorously support the complexity claims. We will revise the manuscript to address these concerns directly.

read point-by-point responses

  1. Referee: §2 (Asymptotic Reduction) and abstract: The central claim that solutions are 'asymptotically governed, in leading order, by linear elliptic equations, coupled with low-dimensional vortex dynamics' is asserted without a quantitative error estimate. No lemma or theorem supplies a concrete rate (e.g., O(δ^k) with δ the core size or inter-vortex distance) that is guaranteed to be smaller than the target accuracy ε. Because the quantum linear solver is applied to the reduced model, the absence of such a bound means the stated complexity for the original nonlinear PDE is not yet justified.

    Authors: We agree that a self-contained quantitative error bound is required. In the revised manuscript we will insert a new lemma in §2 that derives an explicit O(δ) reduction error (δ = vortex core radius) under the standard assumptions of the Ginzburg-Landau regime, together with a short argument showing that this error is o(ε) uniformly for the target accuracy ε of interest. The lemma will be supported by classical vortex-asymptotics references and a brief proof outline. This will make the subsequent complexity statements for the original nonlinear PDE fully justified. revision: yes

  2. Referee: §4 (Complexity Analysis): The exponential improvement in spatial-size dependence for the 2D NLS is derived for the hybrid reduced system. Without an explicit relation showing that the reduction error is o(ε) uniformly, the overall runtime bound claimed for the original nonlinear problem does not follow.

    Authors: We concur. The runtime analysis in §4 currently applies to the reduced hybrid system. We will add a paragraph immediately after the new lemma that composes the O(δ) reduction error with the quantum-solver error, verifies that the total error remains ≤ ε, and therefore transfers the exponential spatial-size improvement to the original nonlinear PDE. Corresponding clarifications will be made in the abstract and the complexity theorem statement. revision: yes

Circularity Check

0 steps flagged

No circularity: hybrid speedup follows from standard quantum linear solver applied to imported asymptotic reduction

full rationale

The paper's central claim rests on utilizing a known asymptotic property of the Ginzburg-Landau-type nonlinear PDEs (solutions governed in leading order by linear elliptic equations plus low-dimensional vortex dynamics). Vortex motion is advanced classically and the elliptic boundary-value problem is solved via quantum algorithms (including BPX preconditioning and Schrödingerization). The exponential improvement in spatial-size dependence is the standard complexity of quantum linear-system solvers applied to the reduced linear problem; it is not obtained by fitting parameters to the algorithm's own outputs, redefining the target quantity, or chaining self-citations that themselves reduce to the present result. No self-definitional, fitted-prediction, or ansatz-smuggling steps appear in the derivation chain. The approach is therefore self-contained against external benchmarks for the linear solver and the cited asymptotic regime.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the domain assumption that the strongly nonlinear regime permits an asymptotic split into linear elliptic problems plus low-dimensional vortex dynamics; no free parameters or invented entities are introduced in the abstract.

axioms (1)
  • domain assumption In the strongly nonlinear regime, solutions to the complex nonlinear Schrödinger equation and related Ginzburg-Landau-type equations are asymptotically governed by linear elliptic equations coupled with low-dimensional vortex dynamics.
    Invoked in the abstract as the foundation for advancing vortex motion classically and solving the elliptic part quantum-mechanically.

pith-pipeline@v0.9.0 · 5539 in / 1440 out tokens · 26294 ms · 2026-05-10T13:49:45.261549+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

41 extracted references · 41 canonical work pages

  1. [1]

    D. An, A. W. Childs, and L. Lin. Quantum algorithm for linear non-unitary dynamics with near-optimal dependence on all parameters. 2023

    work page 2023

  2. [2]

    D. An, J. Liu, and L. Lin. Linear combination of Hamiltonian simulation for non-unitary dynamics with optimal state preparation cost.Phys. Rev. Lett., 131(15):150603, 2023

    work page 2023

  3. [3]

    Bao and Q

    W. Bao and Q. Tang. Numerical study of quantized vortex interactions in the nonlinear schr¨ odinger equation on bounded domains.Multiscale Model. Simul., 12(2):411–439, 2014

    work page 2014

  4. [4]

    D. W. Berry. High-order quantum algorithm for solving linear differential equations.J. Phys. A: Math. Theor., 47(10):105301, 2014

    work page 2014

  5. [5]

    D. W. Berry, A. M. Childs, A. Ostrander, and G. Wang. Quantum algorithm for linear differential equations with exponentially improved dependence on precision.Commun. Math. Phys., 356(3):1057–1081, 2017

    work page 2017

  6. [6]

    J. H. Bramble, J. E. Pasciak, and J. Xu. Parallel multilevel preconditioners.Math. Comp., 55(191):1–22, 1990. 27

    work page 1990

  7. [7]

    Bravyi, R

    S. Bravyi, R. Manson-Sawko, M. Zayats, and S. Zhuk. Quantum simulation of a noisy classical nonlinear dynamics. 2025

    work page 2025

  8. [8]

    A. M. Childs, R. Kothari, and R. D. Somma. Quantum algorithm for systems of linear equa- tions with exponentially improved dependence on precision.SIAM J. Comput., 46(6):1920– 1950, 2017

    work page 1920

  9. [9]

    A. W. Childs and J. Liu. Quantum spectral methods for differential equations.Commun. Math. Phys., 375(2):1427–1457, 2020

    work page 2020

  10. [10]

    B. D. Clader, B. C. Jacobs, and C. R. Sprouse. Preconditioned quantum linear system algo- rithm.Phys. Rev. Lett., 110(25):250504, 2013

    work page 2013

  11. [11]

    D. Dong, Y. Li, and J. Xue. A quantum algorithm for linear autonomous differential equations via Pad´ e approximation.Quantum, 9:1770, 2025

    work page 2025

  12. [12]

    W. E. Dynamics of vortices in Ginzburg–Landau theories with applications to superconduc- tivity.Physica D, 77:383–404, 1994

    work page 1994

  13. [13]

    Gily´ en, Y

    A. Gily´ en, Y. Su, G. H. Low, and N. Wiebe. Quantum singular value transformation and beyond: exponential improvements for quantum matrix arithmetics. InProceedings of the 51st Annual ACM SIGACT Symposium on Theory of Computing, pages 193–204, 2019

    work page 2019

  14. [14]

    A. W. Harrow, A. Hassidim, and S. Lloyd. Quantum algorithm for linear systems of equations. Phys. Rev. Lett., 103(15):150502, 2009

    work page 2009

  15. [15]

    S. Jin, X. Li, N. Liu, and Y. Yu. Quantum simulation for quantum dynamics with artificial boundary conditions.SIAM J. Sci. Comput., 46:B403–B421, 2024

    work page 2024

  16. [16]

    Jin and N

    S. Jin and N. Liu. Quantum algorithms for computing observables of nonlinear partial differ- ential equations. 2022

    work page 2022

  17. [17]

    Jin and N

    S. Jin and N. Liu. Analog quantum simulation of partial differential equations.Quantum Sci. Technol., 2023

    work page 2023

  18. [18]

    Jin and N

    S. Jin and N. Liu. Analog quantum simulation of parabolic partial differential equations using Jaynes–Cummings-like models. 2024

    work page 2024

  19. [19]

    Jin and N

    S. Jin and N. Liu. Quantum algorithms for nonlinear partial differential equations.Bull. Sci. Math., 194:103457, 2024

    work page 2024

  20. [20]

    Jin and N

    S. Jin and N. Liu. Quantum simulation of discrete linear dynamical systems and simple iter- ative methods in linear algebra via Schr¨ odingerisation.Proc. R. Soc. A, 480(2291):20230370, 2024

    work page 2024

  21. [21]

    Jin and N

    S. Jin and N. Liu. Quantum algorithms for viscosity solutions to nonlinear Hamilton–Jacobi equations based on an entropy penalisation method. 2025. preprint. 28

    work page 2025

  22. [22]

    S. Jin, N. Liu, and C. Ma. Quantum simulation of Maxwell’s equations via Schr¨ odingerization. ESAIM Math. Model. Numer. Anal., 58:1853–1879, 2024

    work page 2024

  23. [23]

    S. Jin, N. Liu, and C. Ma. On Schr¨ odingerization based quantum algorithms for linear dy- namical systems with inhomogeneous terms.SIAM J. Numer. Anal., 63(4):1861–1885, 2025

    work page 2025

  24. [24]

    S. Jin, N. Liu, and C. Ma. Schr¨ odingerization based computationally stable algorithms for ill-posed problems in partial differential equations.SIAM J. Sci. Comput., 47(4):B976–B1000, 2025

    work page 2025

  25. [25]

    S. Jin, N. Liu, C. Ma, Y. Peng, and Y. Yu. On the Schr¨ odingerization method for linear non-unitary dynamics with optimal dependence on matrix queries.Commun. Math. Sci., 24(2):565–592, 2026

    work page 2026

  26. [26]

    S. Jin, N. Liu, C. Ma, and Y. Yu. Quantum preconditioning method for linear systems problems via Schr¨ odingerization, 2025. arXiv:2505.06866

    work page arXiv 2025

  27. [27]

    S. Jin, N. Liu, and Y. Yu. Time complexity analysis of quantum difference methods for linear high dimensional and multiscale partial differential equations.J. Comput. Phys., 471:111641, 2022

    work page 2022

  28. [28]

    S. Jin, N. Liu, and Y. Yu. Quantum simulation of partial differential equations: Applications and detailed analysis.Phys. Rev. A, 108:032603, 2023

    work page 2023

  29. [29]

    S. Jin, N. Liu, and Y. Yu. Quantum simulation of partial differential equations via Schr¨ odingerization.Phys. Rev. Lett., 133(23):230602, 2024

    work page 2024

  30. [30]

    Jin and S

    S. Jin and S. Osher. A level set method for the computation of multivalued solutions to quasi- linear hyperbolic PDEs and Hamilton–Jacobi equations.Commun. Math. Sci., 1(3):575–591, 2003

    work page 2003

  31. [31]

    I. Joseph. Koopman–von neumann approach to quantum simulation of nonlinear classical dynamics.Phys. Rev. Res., 2:043102, 2020

    work page 2020

  32. [32]

    F. H. Lin and J. X. Xin. A unified approach to vortex motion laws of complex scalar field equations.Math. Res. Lett., 5(4):455–460, 1998

    work page 1998

  33. [33]

    F. H. Lin and J. X. Xin. On the incompressible fluid limit and the vortex motion law of the nonlinear Schr¨ odinger equation.Commun. Math. Phys., 200(2):249–274, 1999

    work page 1999

  34. [34]

    Lin and Y

    L. Lin and Y. Tong. Optimal polynomial based quantum eigenstate filtering with application to solving quantum linear systems.Quantum, 4:361, 2020

    work page 2020

  35. [35]

    J.-P. Liu, D. An, D. Fang, J. Wang, G. H. Low, and S. Jordan. Efficient quantum algo- rithm for nonlinear reaction-diffusion equations and energy estimation.Commun. Math. Phys., 404(2):963–1020, 2023. 29

    work page 2023

  36. [36]

    J.-P. Liu, H. O. Kolden, H. K. Krovi, N. F. Loureiro, K. Trivisa, and A. M. Childs. Efficient quantum algorithm for dissipative nonlinear differential equations.Proc. Natl. Acad. Sci. U.S.A., 118(35):e2026805118, 2021

    work page 2021

  37. [37]

    C. Ma, S. Jin, N. Liu, K. Wang, and L. Zhang. Schr¨ odingerization based quantum circuits for Maxwell’s equation with time-dependent source terms, 2024

    work page 2024

  38. [38]

    Meng and Y

    Z. Meng and Y. Yang. Quantum computing of fluid dynamics using the hydrodynamic Schr¨ odinger equation.Phys. Rev. Res., 5:033182, 2023

    work page 2023

  39. [39]

    J. C. Neu. Vortices in complex scalar fields.Physica D, 43:385–406, 1990

    work page 1990

  40. [40]

    Shao and H

    C. Shao and H. Xiang. Quantum circulant preconditioner for a linear system of equations. Phys. Rev. A, 98:062321, 2018

    work page 2018

  41. [41]

    Wossnig, Z

    L. Wossnig, Z. Zhao, and A. Prakash. Quantum linear system algorithm for dense matrices. Phys. Rev. Lett., 120(5):050502, 2018. 30

    work page 2018