🔢 mathematics

Exactly solved Schrödinger equations with time-dependent Hamiltonians

This paper presents exact, assumption-free analytical formulas for the evolution operators of four time-dependent Hamiltonians relevant to quantum spin batteries by employing a novel combination of $\star$ -algebra, path-sums, and Omega calculus, thereby enabling the recovery and extension of existing approximations and providing explicit exact formulas for Floquet Hamiltonians at all orders.

Original authors: Michael Warnock, Antônio Francisco Neto, Pierre-Louis Giscard, Omid Faizy, Christian Joachim

Published 2026-07-10

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Michael Warnock, Antônio Francisco Neto, Pierre-Louis Giscard, Omid Faizy, Christian Joachim

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine you have a tiny, two-story quantum house. On the ground floor, a particle lives in a cozy state called $|0\rangle$ . On the second floor, it lives in an excited state called $|1\rangle$ . Usually, to get the particle to jump up to the second floor, you have to shake the house with a specific rhythm. If you shake it too fast, too slow, or with the wrong pattern, the particle just stays put.

For decades, physicists have tried to predict exactly how this particle moves when the shaking gets complicated—like when the rhythm changes over time, or when the house is being jostled by random thermal noise (like a crowd of people bumping into the walls). The old way of doing this was like trying to guess the path of a drunk sailor by taking tiny steps and hoping you don't fall off the boat. Scientists would use "perturbation theory," which is basically a fancy way of saying, "Let's assume the shaking is weak and add up a bunch of small guesses." If the shaking got too strong, or the rhythm got too weird, those guesses would fall apart, and the math would break.

The Big Breakthrough
In this paper, the authors, led by Michael Warnock and his team, say: "No more guessing." They have found the exact, perfect map for how this particle moves in four different scenarios, including some where the house is being shaken by random noise.

They didn't just find an approximation; they found the exact solution. This means their formulas work no matter how hard you shake the house, how fast the rhythm changes, or how much noise is involved. There are no "ifs," "buts," or "assuming the shaking is small." The math works everywhere.

The Magic Toolkit
How did they do it? They didn't use the old, broken ladders. Instead, they combined three brand-new mathematical tools into a super-toolkit:

The $\star$ -Algebra: Think of this as a special translator. It takes the messy, time-changing problem of the shaking house and turns it into a clean, static algebra problem. It's like taking a movie of the particle's journey and turning it into a single, static blueprint that you can solve with a ruler and a calculator.
Path-Sums: Imagine the particle is a hiker trying to get from the ground floor to the second floor. It can take many different paths, walking up stairs, taking elevators, or even going back down and up again. The old methods tried to count every possible path forever, which is impossible. The "Path-Sum" method is like a smart GPS that realizes all those infinite paths can be grouped into a neat, finite list of loops. It turns an infinite mess into a clean, finished fraction.
Omega Calculus: This is the final chef. Once the GPS has the list of loops, Omega Calculus cooks them up. It takes the complicated ingredients (integrals and time-dependent functions) and turns them into simple, delicious sums of standard math functions. It's like turning a complex recipe into a simple list of ingredients you can actually buy at the store.

The Four Scenarios They Solved
The team tested their new map on four different types of "shaking":

The Rhythmic Shaker: The ground floor energy wiggles up and down like a sine wave ( $\cos(\omega_0 t)$ ). They found the exact formula for how the particle jumps, even if the frequency is weird or the shaking is strong.
The Double Shaker: Now, the connection between the floors also wiggles ( $\cos(\omega t)$ ). This is the famous "Rabi" and "Bloch-Siegert" scenario. Their exact formula shows that the old approximations (like the Rotating Wave Approximation) are just the first few terms of their much bigger, more accurate map. They even found a new, exact formula for the "effective Hamiltonian" (the rulebook for the particle's average behavior) that works for all orders, not just the first guess.
The Noisy Shaker: Instead of a smooth rhythm, the ground floor is hit by random white noise (like static on a radio). This models a quantum system sitting on a warm surface. They showed that even with this chaos, the exact solution exists. In the weak noise limit, the particle still jumps, but the "jumpiness" is slightly altered by the noise.
The Gaussian Shaker: The connection between floors is turned on and off like a smooth, bell-shaped pulse (a Gaussian), while the ground floor is still noisy. This is like trying to flip a switch quickly in a storm. They found the exact path the particle takes, showing how the noise can sometimes ruin the smooth flip, but the math still holds perfectly.

What They Ruled Out
The paper is very clear about what their method is not. They are not using approximations. They are not relying on the idea that the shaking must be weak. They explicitly state that their solutions work in regimes where the old "high-frequency" or "weak-coupling" approximations fail completely. If you've been using the old "guessing" methods for strong shaking, this paper says: "Stop. Use the exact map instead."

How Sure Are They?
The authors are as sure as math can be. They didn't just simulate this on a computer and hope it looked right. They proved the formulas are exact.

They derived the solutions using rigorous mathematical theorems (the $\star$ -algebra, path-sums, and Omega calculus).
They checked their math by running computer simulations. In the paper, you can see graphs where the "exact analytical solution" (the red dashed line) and the "fully numerical solver" (the solid blue line) are indistinguishable. They overlap perfectly, even in the most chaotic, noisy, and strong-shaking scenarios.
They proved that their series of numbers always converges (meaning the sum doesn't blow up to infinity) and works for any parameter values.

Why It Matters
This isn't just about two-story houses. The authors show that this "magic toolkit" can be used for any system that changes over time, not just quantum particles. Whether you are designing a quantum battery, controlling a superconducting qubit, or just trying to understand how energy moves in a noisy world, you now have a way to calculate the exact answer without guessing.

They even showed how to use this to design better quantum gates (the switches used in quantum computers). By knowing the exact "effective Hamiltonian," you can correct for errors caused by the "Bloch-Siegert shift" (a tiny wobble in the particle's rhythm) and make your quantum switches more accurate.

In short, the authors took a problem that was usually solved by making messy guesses and replaced it with a clean, exact, and universally valid map. They didn't just improve the map; they drew the whole territory from scratch.

Technical Summary: Exactly Solved Schrödinger Equations with Time-Dependent Hamiltonians

Problem Statement
This work addresses the challenge of finding exact, analytical, and explicit solutions to the non-autonomous Schrödinger equation for two-level quantum systems driven by time-dependent Hamiltonians. The authors focus on four specific physical instances relevant to quantum energy storage (quantum spin batteries), including two deterministic cases and two stochastic cases involving white noise. The Hamiltonians considered take the general form:
$H(t) = \begin{pmatrix} S_0 + \epsilon(t) & g f(t) \\ \bar{g} \bar{f}(t) & S_1 + i\Gamma \end{pmatrix}$
where $\epsilon(t)$ represents a time-dependent energy shift (either periodic or stochastic noise) and $f(t)$ modulates the coupling $g$ . The study aims to determine the evolution operator $U(t)$ , from which transition probabilities and stored energy can be derived, without relying on perturbative expansions, implicit coefficients, or unevaluated integrals. The solutions must remain valid across all parameter regimes, including non-periodic driving, comparable physical parameter magnitudes, and stochastic environments.

Methodology
The authors employ a novel combination of three mathematical techniques to transform the differential problem into a solvable algebraic one:

$\star$ -algebra: This framework converts linear differential systems into ordinary linear algebraic systems using the $\star$ -product (a Volterra composition defined on functions of two variables). The time-ordered exponential solution to the Schrödinger equation is expressed as a matricial $\star$ -resolvent, $G = (Id_\star - H(t)\Theta)^{\star -1}$ .
Method of Path-Sums: This technique interprets the $\star$ -resolvent as a sum of walks on a graph where vertices represent quantum states. By exploiting the combinatorial factorization of walks into simple paths and cycles, the infinite Dyson series is resummed into a finite-depth continued fraction involving only scalar $\star$ -products and resolvents.
Omega Calculus: To evaluate the $\star$ -products and inverses appearing in the path-sum solution, the authors utilize Omega calculus (MacMahon's Partition Analysis). This method maps time-dependent functions to "crude generating functions" in an Omega domain using a linear operator ( $\Omega$ ) that extracts constant terms from Laurent series. Crucially, this process avoids explicit integration, converting iterated integrals of exponentials into divided-difference exponential functions in the time domain.

Key Results
The paper presents exact, explicit formulas for the evolution operators of four specific cases:

Case 1 (Periodic Diagonal Modulation): $\epsilon(t) = \epsilon_0 \cos(\omega_0 t)$ and $f(t) = 1$ . The solution is expressed as an unconditionally convergent series of divided-difference exponentials involving products of Bessel functions. The authors demonstrate how to recover standard high-frequency approximations and derive explicit formulas for effective Hamiltonians at all orders.
Case 2 (Periodic Diagonal and Off-Diagonal Modulation): $\epsilon(t) = \epsilon_0 \cos(\omega_0 t)$ and $f(t) = \cos(\omega t)$ . This case includes the Bloch-Siegert Hamiltonian as a subcase. The authors derive an explicit, commutator-free formula for the effective Hamiltonian $H_{\text{eff}}$ valid to all orders in $\omega^{-1}$ . This formula recovers known Rabi (RWA) and Bloch-Siegert approximations while extending beyond them, providing exact expressions for resonance shifts and effective rotation axes on the Bloch sphere.
Case 3 (Stochastic Diagonal Modulation): $\epsilon(t)$ is white noise and $f(t)$ is either constant or $\cos(\omega t)$ . Using the Wong-Zakai theorem and a Karhunen-Loève expansion for the Brownian motion, the stochastic Schrödinger equation is mapped to a random ordinary differential equation. The solution involves Generalized Bessel Functions (GBF). The authors show that while resonances exist mathematically, the low-frequency nature of the noise-induced driving and the small amplitude of high-order GBFs often prevent significant population transfer within finite timeframes.
Case 4 (Stochastic Modulation with Gaussian Coupling): $\epsilon(t)$ is white noise and $f(t)$ is a Gaussian or modulated Gaussian pulse. The solution demonstrates that Gaussian envelopes can stabilize population in the excited state in the absence of noise. In the presence of noise, the spectral content spreads across GBFs, leading to destructive interference that degrades the smooth population inversion characteristic of the noiseless case.

Significance and Claims
The authors claim that their primary contribution is the provision of unconditionally convergent, analytical, and explicit solutions that do not rely on perturbative assumptions or parameter restrictions. Key aspects of their significance include:

Exactness vs. Approximation: Unlike standard approaches (e.g., Floquet-Magnus, van Vleck, or Brillouin-Wigner perturbation theory) which are asymptotic or valid only in specific regimes (e.g., weak coupling or high frequency), these solutions are valid across the entire parameter space.
Derivation of Approximations: The exact solutions serve as a rigorous foundation from which known approximations (such as the Rotating Wave Approximation and Bloch-Siegert shifts) can be recovered and systematically extended to higher orders.
Commutator-Free Effective Hamiltonians: The paper provides explicit, integral-free formulas for effective Hamiltonians that avoid the proliferation of nested commutators typical of high-frequency expansions.
Stochastic Rigor: The work offers a rigorous analytical treatment of stochastic quantum systems, moving beyond numerical simulations to provide exact series representations that capture noise-induced modifications to phenomena like the Bloch-Siegert shift.
General Applicability: The authors state that the mathematical machinery (the triptych of $\star$ -algebra, path-sums, and Omega calculus) is widely applicable to other non-autonomous differential systems beyond the specific quantum models presented.

The paper concludes that these exact representations allow for the construction of effective descriptions from first principles, free from uncontrolled approximations, and offer a powerful systematic approach for non-autonomous problems in physics and engineering.

Technical Summary: Exactly Solved Schrödinger Equations with Time-Dependent Hamiltonians

More like this