We establish and analyse the performance and resource requirements of an end-to-end fault-tolerant quantum algorithm for computing the absorption spectrum and population dynamics of photoexcited pyrazine. The quantum circuit construction consists of initial state preparation using uniformly controlled rotations, the time-dependent Hamiltonian propagation based on the grid-based Split Operator Quantum Fourier Transform (SO-QFT) method, and cost-effective measurements including statistical and canonical phase estimation. We use classical emulations to validate the quantum resources required for the task, and propose generalised formulae for the qubit count and gate depth calculation. Simulating the vibronic dynamics of pyrazine in a low-dimensional abstraction requires 17-qubit circuits with a gate depth of $\mathcal{O}(10^5)$, whereas a full-dimensional simulation of pyrazine in 24 modes requires at least 97-qubit circuits with a gate depth of $\mathcal{O}(10^6)$. Our work provides a foundational framework for understanding high-dimensional wavepacket-based quantum simulations of photo-induced dynamics and vibronic spectra, anticipating future applications in the simulation of even larger molecular systems on fault-tolerant quantum computers.
We establish and analyse the performance and resource requirements of an end-to-end fault-tolerant quantum algorithm for computing the absorption spectrum and population dynamics of photoexcited pyrazine. The quantum circuit construction consists of initial state preparation using uniformly controlled rotations, the time-dependent Hamiltonian propagation based on the grid-based Split Operator Quantum Fourier Transform (SO-QFT) method, and cost-effective measurements including statistical and canonical phase estimation. We use classical emulations to validate the quantum resources required for the task, and propose generalised formulae for the qubit count and gate depth calculation. Simulating the vibronic dynamics of pyrazine in a low-dimensional abstraction requires 17-qubit circuits with a gate depth of O(105), whereas a full-dimensional simulation of pyrazine in 24 modes requires at least 97-qubit circuits with a gate depth of O(106). Our work provides a foundational framework for understanding high-dimensional wavepacket-based quantum simulations of photo-induced dynamics and vibronic spectra, anticipating future applications in the simulation of even larger molecular systems on fault-tolerant quantum computers.
This thesis investigates the grid-based Split Operator–Quantum Fourier Transform (SO-QFT) scheme as a general-purpose backbone for time-dependent quantum simulation. We develop an end-to-end algorithmic framework that unifies state preparation on discretised grids, SO-QFT time propagation, and the extraction of observables via autocorrelation functions, with a focus on efficient quantum circuit construction. Classical emulation is used to evaluate the performance and resource demands of this framework across diverse chemical contexts. For vibronic dynamics arising from conical intersections, we assemble a complete workflow that reproduces absorption spectra and population transfer for a well-characterised prototypical polyatomic system. On this basis, we estimate the required quantum resources by mapping grid size, time resolution and spectral precision to qubit budgets, circuit depths and measurement counts. The framework is further tailored to infrared spectroscopy through probabilistic encoding of non-unitary dipole operators and explicit treatment of higher-order Hamiltonian terms via multi-controlled operations. To minimise gate depth, we optimise time discretisation and employ approximation strategies for initial states and dipole truncations, while preserving the fidelity of computed spectra. Importantly, the resource scaling relations derived from the vibronic and infrared simulations generalise to comparable molecular systems, providing a transferable template for predicting resource estimates beyond the specific benchmark studied. Finally, the methodology is extended to Coulombic systems, whose singular potentials challenge grid-based representations. We improve the simulation fidelity by proposing two correction schemes: a Hamiltonian correction that accounts for grid-basis structure, and a modification of the initial wavefunction to better match the discretised Hamiltonian. Both corrections are compatible with compact quantum encodings, and their representation in truncated Fourier or Walsh series offers an additional reduction in circuit depth. Collectively, these studies advance the grid-based SO-QFT scheme from a promising construct to a systematically assessed framework, establishing methodological guidance, algorithmic variants and quantitative resource benchmarks for future fault-tolerant quantum applications.
We develop a time-dependent, grid-based framework for simulating infrared spectra that is specifically designed for quantum computers. The proposed circuit employs a probabilistic strategy for applying the non-unitary dipole operator and an Split Operator-Quantum Fourier Transform time evolution scheme. Using a vibrational model of the water molecule as a test system, our classical emulation results demonstrate accurate determination of fundamental and overtone band positions and intensities via Fourier-transformed dipole-dipole autocorrelation functions. We also identify the optimal time parameters that minimise gate depths while maintaining high fidelity. For further resource reduction, we validate the feasibility of utilising harmonic oscillator approximations in state preparation and dipole operator truncations. With its scalability to higher-dimensional normal mode spaces, this wavefunction-based approach establishes a robust foundation for studying IR spectra on future quantum hardware.
First-quantized, grid-based methods for chemistry modeling are a natural and elegant fit for quantum computers. However, it is infeasible to use today's quantum prototypes to explore the power of this approach because it requires a substantial number of near-perfect qubits. Here, we use exactly emulated quantum computers with up to 36 qubits to execute deep yet resource-frugal algorithms that model 2D and 3D atoms with single and paired particles. A range of tasks is explored, from ground state preparation and energy estimation to the dynamics of scattering and ionization; we evaluate various methods within the split-operator QFT (SO-QFT) Hamiltonian simulation paradigm, including protocols previously described in theoretical papers and our own techniques. While we identify certain restrictions and caveats, generally, the grid-based method is found to perform very well; our results are consistent with the view that first-quantized paradigms will be dominant from the early fault-tolerant quantum computing era onward.
Simulations of chemical dynamics are a powerful means for understanding chemistry. However, classical computers struggle to simulate many chemical processes, especially non-adiabatic ones, where the Born-Oppenheimer approximation breaks down. Quantum computers could simulate quantum-chemical dynamics more efficiently than classical computers, but there is currently no complete quantum algorithm for calculating dynamical observables to within a known error. Here, we develop an efficient, end-to-end quantum algorithm for simulating chemical dynamics that avoids all uncontrolled approximations (including the Born-Oppenheimer approximation) and whose error is bounded subject to mild assumptions. To do so, we treat the nuclei and the electrons on an equal footing and simulate the full molecular wavefunction on a momentum-space grid in first quantization, including all algorithmic steps: initial-state preparation, time evolution using qubitization, and measurement of chemical observables such as reaction yields and rates. Our work gives the first algorithm for quantum simulation of chemistry whose end-to-end complexity achieves sublinear scaling in the size of the grid. We achieve this by developing an exponentially faster method for initial-state-preparation. Photochemistry is a likely early application of our algorithm and we estimate resources required for end-to-end simulations of non-adiabatic dynamics of atmospherically important molecules. Classically intractable photochemical computations could be performed using resources comparable to those required for other chemical applications of quantum computing.
Quantum simulation of molecular electronic structure is one of the most promising applications of quantum computing. However, achieving chemically accurate predictions for strongly correlated systems requires quantum phase estimation (QPE) on fault-tolerant quantum computing (FTQC) devices. Existing resource estimates for typical FTQC architectures suggest that such calculations demand millions of physical qubits, thereby placing them beyond the reach of near-term devices. Here, we investigate the feasibility of performing QPE for chemically relevant molecular systems in an early-FTQC regime, characterized by partial fault tolerance, constrained qubit budgets, and limited circuit depth. Our framework is based on single-ancilla, Trotter-based QPE implementations combined with partially randomized time evolution. Within this framework, we develop a novel Hamiltonian optimization strategy, termed unitary weight concentration, that reduces algorithmic cost by reshaping linear-combination-of-unitaries representations. Applying this framework to active-space models of iron-sulfur clusters, cytochrome P450 active sites, and CO2-utilization catalysts, we perform end-to-end resource estimation using the latest version of the space-time efficient analog rotation (STAR) architecture. Our results indicate that ground-state energy estimation for active spaces of approximately 20 to 50 spatial orbitals, well beyond the reach of classical full configuration interaction, is achievable using ∼105 physical qubits, with runtimes on the order of days to weeks. These findings demonstrate that while full-fledged fault-tolerant quantum computers will be required for even larger molecular simulations, chemically meaningful quantum chemistry problems are already within reach in an experimentally relevant, early-FTQC regime.
Abstract We present a comprehensive end-to-end framework for simulating the real-time dynamics of chemical systems on a fault-tolerant quantum computer, incorporating both electronic and nuclear quantum degrees of freedom. An all-particle simulation is nominally efficient on a quantum computer, but practically infeasible. Hence, central to our approach is the construction of a first-quantized plane-wave algorithm making use of pseudoions. The latter consolidate chemically inactive electrons and the nucleus into a single effective dynamical ionic entity, extending the well-established concept of pseudopotentials in quantum chemistry to a two-body interaction. We explicitly describe efficient quantum circuits for initial state preparation across all degrees of freedom, as well as for block-encoding the Hamiltonian describing interacting pseudoions and chemically active electrons, by leveraging recent advances in quantum rejection sampling to optimize the implementations. To extract useful chemical information, we first design molecular fingerprints by combining density-functional calculations with machine learning techniques, and subsequently validate them through surrogate classical molecular dynamics simulations. These fingerprints are then coherently encoded on a quantum computer for efficient molecular identification via amplitude estimation. We provide an extensive analysis of the cost of running the algorithm on a fault-tolerant quantum computer for several chemically interesting systems. As an illustration, simulating the interaction between NH3 and BF3 (a 40-particle system) requires 808 logical qubits to encode the problem, and approximately 10^11 Toffoli gates per femtosecond of time evolution. Our results establish a foundation for further quantum algorithm development targeting chemical and material dynamics.