Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2509.07084 · PRX Quantum

Fourier Neural Operators for Time-Periodic Quantum Systems: Learning Floquet Hamiltonians, Observable Dynamics, and Operator Growth

Zihao Qi, Yang Peng, Christopher Earls

Published Apr 13, 2026• 2 citations
arXivHTMLPDF

Abstract

Time-periodic quantum systems exhibit a rich variety of far-from-equilibrium phenomena and serve as ideal platforms for quantum engineering and control. However, simulating their dynamics with conventional numerical methods remains challenging due to the exponential growth of Hilbert space dimension and rapid spreading of entanglement. In this work, we introduce Fourier neural operators (FNOs) as an efficient, accurate, and scalable framework for nonequilibrium quantum dynamics. Parameterized in Fourier space, FNO naturally captures temporal correlations and remains minimally dependent on discretization of time. We demonstrate the versatility of FNO through three complementary learning paradigms: reconstructing effective Floquet Hamiltonians, predicting expectation values of local observables, and learning quantum information spreading. For each learning task, FNO achieves remarkable accuracy, while attaining a significant speedup, compared to exact numerical methods. Moreover, FNO possesses capabilities beyond that of conventional methods, such as predicting all local observables from a subset of measurements without information about the Hamiltonian, as well as extrapolating beyond the time window provided by training data, enabling access to observables and operator-spreading dynamics that might be beyond the coherence time. By employing a spatially local basis, we argue that the computational cost of FNOs scales only polynomially with the system size. Our results establish FNO as a versatile and scalable computation framework that integrates numerical simulations and experimental data seamlessly, with direct implications for extracting meaningful physics from measurements by near-term quantum computers.

DOI
10.1103/2tt3-yndr

Related papers

Fourier neural operators for learning dynamics in quantum spin systems

Freya Shah, Taylor L. Patti, Julius Berner +3 more

Abstract Fourier Neural Operators (FNOs) excel on functional data tasks, such as those originating from partial differential equations. This renders them effective at simulating the time evolution of quantum wavefunctions, a computationally challenging, yet coveted task. In this manuscript, we use FNOs to model the evolution of random quantum spin systems, chosen due to their representative dynamics. We find that traditional neural networks, such as U-Net, exhibit limited extrapolation beyond the training time interval, whereas FNOs capture the underlying time-evolution operator, generalizing to unseen times. Additionally, we apply FNOs to a compact set of Hamiltonian observables (~ poly( n )) instead of 2 n -component wavefunctions, greatly reducing the size of our inputs, outputs, and FNO models. This Hamiltonian observable-based method demonstrates that FNOs can distill information from high to low-dimensional spaces. We perform numerical experiments on a 20-qubit system, extrapolating Hamiltonian observables to twice the training time with a relative error of 5.8%. Notably, relative to numerical time-evolution methods, FNO achieves an inference speedup of approximately 10 4 × for 20-qubit systems, underscoring its computational efficiency. This extrapolation of observables past training times stands to fundamentally increase the simulatability of quantum systems beyond the limitations of quantum device coherence and tensor network circuit-depth.

arXiv:2409.03302•Communications Physics•Apr 29, 2026•3 citations

Inverse Design of Quantum Control Sequences with Fourier Neural Operators

Anastasia Pipi, Valentin Duruisseaux, Emily Been +4 more

Quantum optimal control is a key tool for steering quantum dynamics, but its computational cost grows rapidly with the Hilbert space dimension. Here, we introduce a Fourier Neural Operator (FNO)-based framework for learning high dimensional molecular quantum dynamics and accelerating the inverse design of control protocols. Given an initial molecular population distribution, laser frequency, and polarization, the FNO predicts molecular-motional population dynamics up to $10^7$ times faster than GPU-accelerated numerical propagation with CUDA-Q Dynamics. Using this fast and differentiable surrogate, we develop the FNO stochastic pulse-measurement planner (FNO-SPMP), which constructs pulse sequences to purify an initially mixed Boltzmann distribution. We demonstrate the protocol in an 888-dimensional subspace of the hydronium molecule at 20 K, achieving a target-state population of 0.98 with a sequence success rate of up to 86.2%. In a shared discrete control space, FNO-SPMP achieves nearly twice the success rate of a reinforcement-learning baseline while using roughly half as many quantum control pulses and reducing pulse-sequence generation time from approximately 10 hours to 10-20 minutes. These results show that operator-learning surrogates can enable inverse design in quantum systems whose Hilbert spaces are too large for conventional direct optimization.

quant-pharXiv:2608.03702•Aug 5, 2026•0 citations

Floquet-informed Learning of Periodically Driven Hamiltonians

KaiMeng Li

Characterizing time-periodic Hamiltonians is pivotal for validating and controlling driven quantum platforms, yet prevailing and unadjusted reconstruction methods demand dense time-domain sampling and heavy post-processing. We introduce a scalable Floquet-informed learning algorithm that represents the Hamiltonian as a truncated Fourier series and recasts parameter estimation as a compact linear inverse problem in the Floquet band picture. The algorithm is well-suited to problems satisfying mild smoothness/band-limiting. In this regime, its sample and runtime complexities scale polynomially with the Fourier cutoff, time resolution, and the number of unknown coefficients. For local Hamiltonian models, the coefficients grows polynomially with system size, yielding at most polylogarithmic dependence on the Hilbert-space dimension. Furthermore, numerical experiments on one- and two-dimensional Ising and Heisenberg lattices show fast convergence to time resolution and robustness to higher-order perturbations. An adaptive rule learns the cutoff of Fourier series, removing the need to set known truncation \textit{a priori}. These features enable practical certification and benchmarking of periodically driven platforms with rapidly decaying higher-order content, and extend naturally to near-periodic drives.

arXiv:2509.02331•Sep 2, 2025•0 citations

Learning Quantum Operator Dynamics from Short-Time Data

Jinyang Li, Satoshi Iso, Shunji Matsuura +2 more

Real-time dynamics of quantum observables provide direct access to excitation spectra and correlation functions in quantum many-body systems, but currently available quantum devices are limited to short evolution times due to decoherence. We propose a neural ordinary differential equation (Neural ODE) framework with a physics-driven neural network to reconstruct long-time operator dynamics from short-time quantum data. By expanding observables in the Pauli basis and exploiting locality and symmetry constraints, the operator evolution is reduced to a tractable set of coefficients whose dynamics are learned from the measurement data. Applied to the Ising spin chain with mixed fields, the operator learning method accurately extrapolates its long-time correlation function and resolves its excitation spectrum. We validate the framework by collecting short-time data from IBM's superconducting quantum processor and successfully reconstruct the Fourier spectrum of an $N=8$ Ising spin chain with 5 times enhanced precision. Our results demonstrate a scalable and data-efficient strategy for extracting dynamical and spectral information from practical quantum hardware.

quant-pharXiv:2603.14699•Mar 17, 2026•0 citations

Toward the Thermodynamic Limit: Neural Operators for Non-equilibrium Dynamics of Mott Insulators

Miles Waugh, Chuwei Wang, Radu Andrei +4 more

Mott insulators exhibit complex photoexcitation dynamics under intense optical driving, with potential implications for carrier multiplication beyond the Shockley-Queisser limit. Probing these nonequilibrium processes requires access to the thermodynamic limit, where the number of lattice sites becomes arbitrarily large, but conventional solvers are constrained to small systems due to the exponential growth of the Hilbert space. Fourier Neural Operators (FNOs), originally developed for solving partial differential equations, naturally accommodate inputs of varying resolution and are capable of capturing nonlocal effects. Here, we employ FNOs to learn the mapping from noise-perturbed ground-state momentum distributions to their post-pulse counterparts across a range of interaction strengths and driving parameters. Trained only on small lattices, the model generalizes zero-shot to much larger systems, producing physically reasonable momentum distributions well beyond the reach of numerical solvers. Specifically, the model can predict momentum distribution for a 1024x1024 system within a few seconds that matches the theoretical behavior of key observables, whereas direct numerical simulations have so far been restricted to edge sizes of ~30. These results demonstrate the potential of neural operators to directly access large-scale nonequilibrium dynamics, providing a new pathway toward the thermodynamic limit in strongly correlated materials.

arXiv:2602.19484•Feb 23, 2026•0 citations

Universal Neural Propagator: Learning Time Evolution in Many-Body Quantum Systems

Zihao Qi, Christopher Earls, Yang Peng

Conventional approaches to simulating quantum many-body dynamics produce a single trajectory: if the Hamiltonian or the initial state is changed, the computation must be re-performed. Recent efforts toward foundation models have begun to address this limitation, yet existing methods transfer across either Hamiltonians or initial states, but not both. In this work, we introduce the Universal Neural Propagator (UNP), a single, unified model that learns the functional mapping from driving protocols to time-evolution propagators. Trained in an entirely self-supervised way, a single UNP model predicts dynamics across a function space of driving protocols and an exponentially large Hilbert space of initial states simultaneously. We benchmark on a two-dimensional driven Ising model and demonstrate the UNP's accuracy and transferability across product and entangled initial states, as well as for both in- and out-of-distribution driving protocols. The UNP remains accurate at system sizes beyond exact diagonalization, and can be efficiently fine-tuned across all initial states using observable data. By shifting the object of learning from quantum states to operators, this work opens a route toward transferable simulation of driven quantum matter.

quant-pharXiv:2605.05299•May 8, 2026•0 citations
Time-periodic quantum systems exhibit a rich variety of far-from-equilibrium phenomena and serve as ideal platforms for quantum engineering and control. However, simulating their dynamics with conventional numerical methods remains challenging due to the exponential growth of Hilbert space dimension and rapid spreading of entanglement. In this work, we introduce Fourier neural operators (FNOs) as an efficient, accurate, and scalable framework for nonequilibrium quantum dynamics. Parameterized in Fourier space, FNO naturally captures temporal correlations and remains minimally dependent on discretization of time. We demonstrate the versatility of FNO through three complementary learning paradigms: reconstructing effective Floquet Hamiltonians, predicting expectation values of local observables, and learning quantum information spreading. For each learning task, FNO achieves remarkable accuracy, while attaining a significant speedup, compared to exact numerical methods. Moreover, FNO possesses capabilities beyond that of conventional methods, such as predicting all local observables from a subset of measurements without information about the Hamiltonian, as well as extrapolating beyond the time window provided by training data, enabling access to observables and operator-spreading dynamics that might be beyond the coherence time. By employing a spatially local basis, we argue that the computational cost of FNOs scales only polynomially with the system size. Our results establish FNO as a versatile and scalable computation framework that integrates numerical simulations and experimental data seamlessly, with direct implications for extracting meaningful physics from measurements by near-term quantum computers.
Abstract Fourier Neural Operators (FNOs) excel on functional data tasks, such as those originating from partial differential equations. This renders them effective at simulating the time evolution of quantum wavefunctions, a computationally challenging, yet coveted task. In this manuscript, we use FNOs to model the evolution of random quantum spin systems, chosen due to their representative dynamics. We find that traditional neural networks, such as U-Net, exhibit limited extrapolation beyond the training time interval, whereas FNOs capture the underlying time-evolution operator, generalizing to unseen times. Additionally, we apply FNOs to a compact set of Hamiltonian observables (~ poly( n )) instead of 2 n -component wavefunctions, greatly reducing the size of our inputs, outputs, and FNO models. This Hamiltonian observable-based method demonstrates that FNOs can distill information from high to low-dimensional spaces. We perform numerical experiments on a 20-qubit system, extrapolating Hamiltonian observables to twice the training time with a relative error of 5.8%. Notably, relative to numerical time-evolution methods, FNO achieves an inference speedup of approximately 10 4 × for 20-qubit systems, underscoring its computational efficiency. This extrapolation of observables past training times stands to fundamentally increase the simulatability of quantum systems beyond the limitations of quantum device coherence and tensor network circuit-depth.
Quantum optimal control is a key tool for steering quantum dynamics, but its computational cost grows rapidly with the Hilbert space dimension. Here, we introduce a Fourier Neural Operator (FNO)-based framework for learning high dimensional molecular quantum dynamics and accelerating the inverse design of control protocols. Given an initial molecular population distribution, laser frequency, and polarization, the FNO predicts molecular-motional population dynamics up to 10710^7107 times faster than GPU-accelerated numerical propagation with CUDA-Q Dynamics. Using this fast and differentiable surrogate, we develop the FNO stochastic pulse-measurement planner (FNO-SPMP), which constructs pulse sequences to purify an initially mixed Boltzmann distribution. We demonstrate the protocol in an 888-dimensional subspace of the hydronium molecule at 20 K, achieving a target-state population of 0.98 with a sequence success rate of up to 86.2%. In a shared discrete control space, FNO-SPMP achieves nearly twice the success rate of a reinforcement-learning baseline while using roughly half as many quantum control pulses and reducing pulse-sequence generation time from approximately 10 hours to 10-20 minutes. These results show that operator-learning surrogates can enable inverse design in quantum systems whose Hilbert spaces are too large for conventional direct optimization.
Characterizing time-periodic Hamiltonians is pivotal for validating and controlling driven quantum platforms, yet prevailing and unadjusted reconstruction methods demand dense time-domain sampling and heavy post-processing. We introduce a scalable Floquet-informed learning algorithm that represents the Hamiltonian as a truncated Fourier series and recasts parameter estimation as a compact linear inverse problem in the Floquet band picture. The algorithm is well-suited to problems satisfying mild smoothness/band-limiting. In this regime, its sample and runtime complexities scale polynomially with the Fourier cutoff, time resolution, and the number of unknown coefficients. For local Hamiltonian models, the coefficients grows polynomially with system size, yielding at most polylogarithmic dependence on the Hilbert-space dimension. Furthermore, numerical experiments on one- and two-dimensional Ising and Heisenberg lattices show fast convergence to time resolution and robustness to higher-order perturbations. An adaptive rule learns the cutoff of Fourier series, removing the need to set known truncation \textit{a priori}. These features enable practical certification and benchmarking of periodically driven platforms with rapidly decaying higher-order content, and extend naturally to near-periodic drives.
Real-time dynamics of quantum observables provide direct access to excitation spectra and correlation functions in quantum many-body systems, but currently available quantum devices are limited to short evolution times due to decoherence. We propose a neural ordinary differential equation (Neural ODE) framework with a physics-driven neural network to reconstruct long-time operator dynamics from short-time quantum data. By expanding observables in the Pauli basis and exploiting locality and symmetry constraints, the operator evolution is reduced to a tractable set of coefficients whose dynamics are learned from the measurement data. Applied to the Ising spin chain with mixed fields, the operator learning method accurately extrapolates its long-time correlation function and resolves its excitation spectrum. We validate the framework by collecting short-time data from IBM's superconducting quantum processor and successfully reconstruct the Fourier spectrum of an N=8N=8N=8 Ising spin chain with 5 times enhanced precision. Our results demonstrate a scalable and data-efficient strategy for extracting dynamical and spectral information from practical quantum hardware.
Mott insulators exhibit complex photoexcitation dynamics under intense optical driving, with potential implications for carrier multiplication beyond the Shockley-Queisser limit. Probing these nonequilibrium processes requires access to the thermodynamic limit, where the number of lattice sites becomes arbitrarily large, but conventional solvers are constrained to small systems due to the exponential growth of the Hilbert space. Fourier Neural Operators (FNOs), originally developed for solving partial differential equations, naturally accommodate inputs of varying resolution and are capable of capturing nonlocal effects. Here, we employ FNOs to learn the mapping from noise-perturbed ground-state momentum distributions to their post-pulse counterparts across a range of interaction strengths and driving parameters. Trained only on small lattices, the model generalizes zero-shot to much larger systems, producing physically reasonable momentum distributions well beyond the reach of numerical solvers. Specifically, the model can predict momentum distribution for a 1024x1024 system within a few seconds that matches the theoretical behavior of key observables, whereas direct numerical simulations have so far been restricted to edge sizes of ~30. These results demonstrate the potential of neural operators to directly access large-scale nonequilibrium dynamics, providing a new pathway toward the thermodynamic limit in strongly correlated materials.
Conventional approaches to simulating quantum many-body dynamics produce a single trajectory: if the Hamiltonian or the initial state is changed, the computation must be re-performed. Recent efforts toward foundation models have begun to address this limitation, yet existing methods transfer across either Hamiltonians or initial states, but not both. In this work, we introduce the Universal Neural Propagator (UNP), a single, unified model that learns the functional mapping from driving protocols to time-evolution propagators. Trained in an entirely self-supervised way, a single UNP model predicts dynamics across a function space of driving protocols and an exponentially large Hilbert space of initial states simultaneously. We benchmark on a two-dimensional driven Ising model and demonstrate the UNP's accuracy and transferability across product and entangled initial states, as well as for both in- and out-of-distribution driving protocols. The UNP remains accurate at system sizes beyond exact diagonalization, and can be efficiently fine-tuned across all initial states using observable data. By shifting the object of learning from quantum states to operators, this work opens a route toward transferable simulation of driven quantum matter.