REVIEW 1 major objections 55 references
Exactly solved Schr\"odinger equations with time-dependent Hamiltonians
T0 review · 1 major / 0 minor · reviewed 2026-07-10 · grok-4.5
Pith's one-line read Four time-dependent 2×2 Hamiltonians, including noisy ones, now have exact evolution operators written as convergent series of ordinary functions.
desk verdict Exact series for four driven/noisy two-level systems plus an all-orders commutator-free Floquet formula; the stochastic half leans on a cited Wong–Zakai limit without a self-contained remainder check. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The star-product algebra converts the time-ordered exponential into a matrix star-resolvent; path-sums rewrite every entry of that resolvent as a finite continued fraction of scalar star-products; Omega calculus evaluates the fractions into series of divided-difference exponentials.
What would settle it
Direct numerical integration of any of the four Hamiltonians for a generic set of parameters must reproduce, to machine precision, the truncated series given in the corresponding equation of the paper; a clear mismatch at moderate truncation order would falsify the claimed exactness.
Extended reading notes
Core claim
The evolution operators of the four Hamiltonians (constant and cosine couplings, each with either a deterministic or a stochastic diagonal drive) are given exactly by the series of divided-difference exponentials written in Eqs. (32), (42), (65), (70), (81) and (88). From those series an explicit all-orders formula for the Floquet effective Hamiltonian follows (Eq. 48).
Load-bearing premise
The noisy solutions rest on the claim that a finite Karhunen–Loève truncation of Brownian motion, once solved exactly, converges almost surely to the original Stratonovich equation when the truncation rank goes to infinity.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript derives exact, explicit, unconditionally convergent series expressions for the evolution operators of four families of 2 imes2 time-dependent Hamiltonians relevant to quantum spin batteries (Eqs. 24, 38, 61/67, 76/85). The solutions (Eqs. 32, 42, 65, 70, 81, 88) are obtained by converting the Schrödinger equation into a ★-resolvent, evaluating the resolvent via path-sums, and reducing the resulting ★-products to series of divided-difference exponentials with Omega calculus. From the exact series the authors recover known RWA/Bloch–Siegert approximations, extract resonance conditions and effective Rabi frequencies, and supply an explicit all-orders, commutator-free formula for the Floquet Hamiltonian (Eq. 48). Two of the models are stochastic; they are treated by Karhunen–Loève truncation of Brownian motion followed by an appeal to the Wong–Zakai theorem.
Significance. If the derivations hold, the work supplies the first closed-form, non-perturbative evolution operators for a set of physically relevant driven and noisy two-level systems, valid across the entire parameter space (including non-periodic and strongly driven regimes). The explicit Floquet formula (Eq. 48) and the systematic extraction of multi-photon resonances and noise-renormalized Bloch–Siegert shifts are concrete advances over existing high-frequency and Magnus expansions. The combination of ★-algebra, path-sums and Omega calculus is presented as a general toolkit for non-autonomous linear systems; the appendices contain the necessary induction proofs and reductions to known limits, and the numerical comparisons (Figs. 1–8) provide independent verification of the truncated series.
major comments (1)
- Section V and Appendix I invoke the Wong–Zakai theorem to pass from the finite-rank Karhunen–Loève random ODEs to the original Stratonovich SDEs, yet supply no remainder estimate or continuity argument showing that the specific ★-resolvents and divided-difference series remain continuous in the truncation topology. The deterministic solutions and the Floquet formula (Eq. 48) are unaffected, but the claim of exact, assumption-free evolution operators for the two stochastic Hamiltonians (Eqs. 61/67, 76/85) is therefore only formal until this gap is closed or the claim is appropriately qualified.
Circularity Check
No load-bearing circularity; exact series follow from Schrödinger equation via independently established ★/path-sum/Omega tools, with only ordinary self-citation of the foundational method papers.
full rationale
The derivation chain begins from the non-autonomous Schrödinger equation (1)–(2), rewrites the evolution operator as a ★-resolvent (11), evaluates the resolvent entries by the path-sum theorem on the two-vertex graph (Appendix A, Eqs. (A3)–(A15)), and converts the resulting ★-products/resolvents into unconditionally convergent series of divided-difference exponentials by Omega calculus (Appendix C, Identities 1–5 and induction (C11)). None of these steps assumes the final series (32), (42), (65), (70), (81) or (88); the series are obtained by direct evaluation. The Floquet Hamiltonian (48) is likewise extracted from the already-derived exact U_rot(T) by the elementary limit T o0 after setting t=T (Appendix G). Self-citations appear only for the three tools themselves (★-algebra [35], path-sums [17,18,20], Omega calculus [13–15]); those prior works rest on combinatorial and analytic arguments independent of the present physical models and are re-derived in the appendices for the 2 imes2 case. The stochastic half invokes the external Wong–Zakai theorem after a Karhunen–Loève truncation; that is a completeness issue, not a circular reduction of the claimed series to their own inputs. No parameters are fitted to data and re-labeled as predictions, no uniqueness theorem is imported to forbid alternatives, and no known empirical pattern is merely renamed. The paper is therefore free of the six circularity patterns.
Assumptions & free parameters
assumptions (4)
- standard math The ★-product turns the time-ordered exponential into a matricial ★-resolvent (Eq. 11).
- standard math Path-sums express every entry of a ★-resolvent as a finite continued fraction over simple cycles of the interaction graph.
- standard math Omega calculus evaluates ★-products and ★-resolvents by converting them into rational functions whose constant terms are divided-difference exponentials.
- domain assumption Wong–Zakai theorem: smooth approximations of Brownian motion converge to the Stratonovich SDE.
Cite this review
Pith. "Pith review of Exactly solved Schr\"odinger equations with time-dependent Hamiltonians." pith.science (2026). https://pith.science/paper/53TNQ34D
@misc{pith2026260708450,
author = {Pith},
title = {Pith review of: Exactly solved Schr\"odinger equations with time-dependent Hamiltonians},
year = {2026},
howpublished = {\url{https://pith.science/paper/53TNQ34D}},
note = {Machine review of arXiv:2607.08450}
}
abstract
We present the analytical, exact, explicit, and assumption free formulas for the evolution operators corresponding to four instances of time-dependent Hamiltonians relevant to quantum spin batteries including two stochastic cases. We demonstrate how to recover and go beyond existing expansions and approximations directly from the exact solutions giving, for example, an explicit exact formula for Floquet Hamiltonians at all orders. The exact solutions are obtained through a completely novel combination of three mathematical techniques, the $\star$-algebra, path-sums and Omega calculus, which we briefly overview. These are widely applicable to other non-autonomous differential systems.
Figures
Figures from the paper (5 more)
Reference graph
Works this paper leans on
-
[1]
General casee 0 ‰0. Just as in the previous case and for the same reasons, the quantity max t ∆Eptqundergoes resonances as bothω 0 andωare tuned. As indicated by the exact solution in Eq. (42b), the resonances occur in the presence of repeated 15 n=-1 n=-2 n=0 n=1 n=2 n=3 n=4 n=5 0.0 0.5 1.0 1.5 2.0 2.5 3.00.0 0.5 1.0 1.5 2.0 2.5 3.0 ω/ε 0 ω0/ε0 0.2 0.5 1...
-
[2]
Physical Interpretation Similarly to Case 1, the term of orderkinU 11 involves a product of 2kGeneralized Bessel functions and divided- difference exponentials with 2k`1 arguments. ForU 12, the term of orderkcomprises a product of 2k`1 Generalized Bessel Functions (GBF) and divided-difference exponentials with 2k`2 arguments. The series given here areunco...
-
[3]
The difference with Case 1 here is that the derived amplitude and frequencies are noise-dependent and now vary slightly over different realizations of the noisy parameters,Z i. Even for moderate increases of the noise, more GBF terms must be considered, invalidating the above approximation and requiring novel re-summation strategies or a full retreat to t...
-
[4]
Physical Interpretations: resonances and the role of the noise For zero noise (γ“0), we recover the Bloch-Siegert Hamiltonian and its unitary evolution, reproducing for small gthe Rabi resonance frequency Ω eff “ |g{2|. In other regimes—specifically when analyzing the population transfer from the ground to the excited state and the resonances of the Bloch...
-
[5]
Physical interpretation Forγ“0, the exact solutions as described by Eqs. (81) simplifies to U11ptq “ ÿ kě0 ˆc π 2 1 σ ˙2k ÿ mmmk,nnnk cnnnk,mmmk eirA1,ε0`B1,A2,ε0`B2,...,Ak,ε0`Bk,0st,(82a) U12ptq “ ÿ kě0 ˆc π 2 1 σ ˙2k`1 ÿ n,mmmk,nnnk eipε0`nωηqt ˆc n,nnnk,mmmk e´irA1,ε0`B1,A2,ε0`B2,...,Ak,ε0`Bk,ε0`nωη,0st.(82b) In these equations, we define cnnnk,mmmk :“...
-
[6]
8ÿ k“´8 Jn´ℓkpxqJkpyq,(B4) whereℓmay be any integer and we used the notationJ npx, yq
Physical Interpretation As considered previously, this model undergoes numerous resonances whose mathematical signature is the presence of repeated arguments in the divided-difference exponentials of the path-sum kernel Eqs. (87). Given thatα“ ˘1 24 andnPZwe obtain ε0 ˘ω ωη “nPZzt0u,orn“0 that isω“ε 0.(89) The mathematical analysis of the resonance is sim...
-
[7]
λn p1´a 1{λq p1´a 2{λq ¨ ¨ ¨ p1´a m{λq “
Preliminary results concerning Omega calculus We here collect four identities pertaining to Omega calculus that are useful in the proofs of the results of the main text. 28 •Identity 1 (Multiplication invariance under the Omega operator) Multiplying Omega variables by auxiliary vari- ables does not change the result after the elimination of those variable...
-
[8]
General treatment ofU ij The proof of the form for the solution relies on the‹-Neumann series expansion for the‹-resolvents appearing in any of theU ij. First of all, we observe that in all cases the diagonal termsU 11 “Θ‹ p1 ‹ ´KΘq ‹´1 andU 22 are always produced by path-sum kernelsKof the form Kpt, sq “ Nÿ j“1 cj eaj terbj ,0spt´sq.(C8) It follows that ...
Show all 55 references
-
[9]
Abramowitz and I.A
M. Abramowitz and I.A. Stegun.Handbook of Mathematical Functions: With Formulas, Graphs, and Mathematical Tables. Applied mathematics series. Dover Publications, 1965. 40
1965
-
[10]
G. E. Andrews, P. Paule, and A. Riese. MacMahon’s Partition Analysis: The Omega Package.European Journal of Combinatorics, 22(7):887–904, 2001
2001
-
[11]
Barchielli, C
A. Barchielli, C. Pellegrini, and F. Petruccione. Stochastic schr¨ odinger equations with coloured noise.EPL (Europhysics Letters), 91(2):24001, 2010
2010
-
[12]
B. J. C. Baxter and R. Brummelhuis. Functionals of exponential brownian motion and divided differences.Journal of Computational and Applied Mathematics, 236(4):424–433, 2011. International Workshop on Multivariate Approximation and Interpolation with Applications (MAIA 2010)
2011
-
[13]
Bialynicki-Birula, B
I. Bialynicki-Birula, B. Mielnik, and J. Pleba´ nski. Explicit solution of the continuous baker-campbell-hausdorff problem and a new expression for the phase operator.Annals of Physics, 51(1):187–200, 1969
1969
-
[14]
Biggs.Algebraic Graph Theory
N. Biggs.Algebraic Graph Theory. Cambridge University Press, Cambridge, 2nd edition, 1993
1993
-
[15]
Bloch and A
F. Bloch and A. Siegert. Magnetic resonance for nonrotating fields.Physical Review, 57:522–527, 1940
1940
-
[16]
D. P. Burum. Magnus expansion generator.Phys. Rev. B, 24:3684–3692, Oct 1981
1981
-
[17]
Dattoli and A
G. Dattoli and A. Torre.Theory and applications of generalized Bessel functions. Aracne, 1996
1996
-
[18]
A. Dey, D. Lonigro, K. Yuasa, and D. Burgarth. Error bounds for the Floquet-Magnus expansion and their application to the semiclassical quantum Rabi model.Phys. Rev. A, 112:053723, Nov 2025
2025
-
[19]
Eckardt and E
A. Eckardt and E. Anisimovas. High-frequency approximation for periodically driven quantum systems from a Floquet- space perspective.New Journal of Physics, 17(9):093039, sep 2015
2015
-
[20]
Flajolet and R
P. Flajolet and R. Sedgewick.Analytic Combinatorics. Cambridge University Press, Cambridge, 1st edition, 2009
2009
-
[21]
Francisco Neto
A. Francisco Neto. Matrix Analysis and Omega Calculus.SIAM Review, 62(1):264–280, 2020
2020
-
[22]
Francisco Neto
A. Francisco Neto. A basis- and integral-free representation of time-dependent perturbation theory via the Omega matrix calculus.Annales de l’Institut Henri Poincare D, 11(2):383–407, July 2023
2023
-
[23]
Francisco Neto and B
A. Francisco Neto and B. M. Villegas-Mart´ ınez. A new basis- and integral-free approach to perturbation theory: The Schr¨ odinger dynamics of N trapped ions in the high-intensity regime.Annals of Physics, 474:169917, 2025
2025
-
[24]
Gemme, M
G. Gemme, M. Grossi, S. Vallecorsa, M. Sassetti, and D. Ferraro. Qutrit quantum battery: Comparing different charging protocols.Phys. Rev. Res., 6:023091, Apr 2024
2024
-
[25]
Giscard and C
P.-L. Giscard and C. Bonhomme. Dynamics of quantum systems driven by time-varying Hamiltonians: Solution for the Bloch-Siegert Hamiltonian and applications to NMR.Physical Review Research, 2:023081, 2020
2020
-
[26]
Giscard, K
P.-L. Giscard, K. Lui, S. J. Thwaite, and D. Jaksch. An exact formulation of the time-ordered exponential using path-sums. Journal of Mathematical Physics, 56(5):053503, 2015
2015
-
[27]
Giscard and S
P.-L. Giscard and S. Pozza. Lanczos-like algorithm for the time-ordered exponential: The‹-inverse problem.Applications of Mathematics, 65:807–827, 2020
2020
-
[28]
Giscard, S
P.-L. Giscard, S. J. Thwaite, and D. Jaksch. Walk-sums, continued fractions and unique factorisation on digraphs. arXiv:1202.5523 [cs.DM], 2012
2012 arXiv
-
[29]
Haeberlen and J
U. Haeberlen and J. S. Waugh. Coherent averaging effects in Magnetic Resonance.Phys. Rev., 175:453–467, Nov 1968
1968
-
[30]
G.-N. Han. A General Algorithm for the MacMahon Omega Operator.Annals of Combinatorics, 7(4):467–480, Dec 2003
2003
-
[31]
Stochastic normalizing flows.arXiv preprint arXiv:2002.09547, 2020
Liam Hodgkinson, Chris van der Heide, Fred Roosta, and Michael W Mahoney. Stochastic normalizing flows.arXiv preprint arXiv:2002.09547, 2020
2002 arXiv
-
[32]
Kalev and I
A. Kalev and I. Hen. An integral-free representation of the Dyson series using divided differences.New Journal of Physics, 23(10):103035, oct 2021
2021
-
[33]
Korenev.Bessel Functions and Their Applications
B.G. Korenev.Bessel Functions and Their Applications. CRC Press, 2019
2019
-
[34]
H. J. Korsch, A. Klumpp, and D. Witthaut. On two-dimensional bessel functions.Journal of Physics A: Mathematical and General, 39(48):14947, nov 2006
2006
-
[35]
M. Lesch. Divided differences in noncommutative geometry: Rearrangement lemma, functional calculus and expansional formula.J. Noncommut. Geom., 23(11):193–223, 2017
2017
-
[36]
Y. Li, D. Sauzin, and S. Sun. The baker–campbell–hausdorff formula via mould calculus.Letters in Mathematical Physics, 109(3):725–746, Mar 2019
2019
-
[37]
I. G. MacDonald.Symmetric functions and Hall polynomials.Oxford: Clarendon Press, 2nd ed. edition, 1998
1998
-
[38]
P. A. MacMahon.Combinatory Analysis, Volumes I and II. AMS Chelsea Publishing. AMS Chelsea Pub., 2001
2001
-
[39]
Mikami, S
T. Mikami, S. Kitamura, K. Yasuda, N. Tsuji, T. Oka, and H. Aoki. Brillouin-Wigner theory for high-frequency expansion in periodically driven systems: Application to Floquet topological insulators.Phys. Rev. B, 93:144307, Apr 2016
2016
-
[40]
Milne-Thomson.The Calculus of Finite Differences
L.M. Milne-Thomson.The Calculus of Finite Differences. AMS Chelsea Publishing Series. AMS Chelsea Pub., 2000
2000
-
[41]
Pozza and N
S. Pozza and N. Van Buggenhout. A new matrix equation expression for the solution of non-autonomous linear systems of ODEs.PAMM, 22:e202200117, 2023
2023
-
[42]
Pozza and N
S. Pozza and N. Van Buggenhout. A new legendre polynomial-based approach for non-autonomous lin- ear odes.ETNA - Electronic Transactions on Numerical Analysis, pages 292–326, 2024. Online available: https://epub.oeaw.ac.at/?arp=0x003f234e - Last access:16.6.2026
2024
-
[43]
Ryckebusch, A
M. Ryckebusch, A. Bouhamidi, and P.-L. Giscard. A Fr´ echet Lie group on distributions.J. Math. Anal. Appl., 546(1):129195, 2025
2025
-
[44]
Stochastic schr¨ odinger equations for markovian and non-markovian cases.Open Systems & Information Dynamics, 21(01n02):1440008, 2014
I Semina, Vitalii Semin, Francesco Petruccione, and Alberto Barchielli. Stochastic schr¨ odinger equations for markovian and non-markovian cases.Open Systems & Information Dynamics, 21(01n02):1440008, 2014
2014
-
[45]
J. H. Shirley. Solution of the Schr¨ odinger equation with a Hamiltonian periodic in time.Physical Review, 138:B979–B987, 1965
1965
-
[46]
Tuorila, M
J. Tuorila, M. Silveri, M. Sillanp¨ a¨ a, E. Thuneberg, Y. Makhlin, and P. Hakonen. Stark effect and generalized Bloch–Siegert shift in a strongly driven two-level system.Phys. Rev. Lett., 105:257003, 2010. 41
2010
-
[47]
Wong-zakai approximations for stochastic differential equations.Acta Applicandae Mathematica, 43(3):317–359, 1996
Krystyna Twardowska. Wong-zakai approximations for stochastic differential equations.Acta Applicandae Mathematica, 43(3):317–359, 1996
1996
-
[48]
M. Vogl, P. Laurell, A. D. Barr, and G. A. Fiete. Flow equation approach to periodically driven quantum systems.Physical Review X, 9:021037, 2019
2019
-
[49]
Volterra and J
V. Volterra and J. P´ er` es.Le¸ cons sur la composition et les fonctions permutables. Gauthier-Villars, Paris, 1924. Recent edition by J. Gabay Ed., 2008, isbn 9782876473232
1924
-
[50]
Warnock, D
M. Warnock, D. A. Hague, and V. F. Mitrovic. Multi-directional periodic driving of a two-level system beyond floquet formalism, 2025
2025
-
[51]
Polynomial chaos functions and stochastic differential equations.Annals of Nuclear Energy, 33(9):774–785, 2006
MMR Williams. Polynomial chaos functions and stochastic differential equations.Annals of Nuclear Energy, 33(9):774–785, 2006
2006
-
[52]
´Etude th´ eorique et exp´ erimentale des transitions ` a plusieurs quanta entre les sous-niveaux Zeeman d’un atome.Ann
Winter, J.-M. ´Etude th´ eorique et exp´ erimentale des transitions ` a plusieurs quanta entre les sous-niveaux Zeeman d’un atome.Ann. Phys., 13(4):745–811, 1959
1959
-
[53]
Y. Yan, Z. L¨ u, and H. Zheng. Bloch-Siegert shift of the rabi model.Physical Review A, 91:053834, 2015
2015
-
[54]
Q. Zeng, N. Ezzell, A. Babakhani, I. Hen, and L. Barash. Inequalities, identities, and bounds for divided differences of the exponential function, 2025
2025
-
[55]
Zeuch, F
D. Zeuch, F. Hassler, J. J. Slim, and D. P. DiVincenzo. Exact rotating wave approximation.Annals of Physics, 423:168327, 2020
2020
Reviewed July 10, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.