Quantum computation is an attractive front for many problems that are intractable for computers today. One such problem is nonadiabatic quantum molecular dynamics, where quantized internal states coupling to parameterized modes result in a Hamiltonian resistant to oracle-based models and spectral decomposition. This dissertation applies diabatic Hamiltonian operators directly to the computational basis as first-quantized split-operator propagators, validated with dynamic observables including absorption and recurrence spectra, scattering cross-sections, population dynamics, and quantum scars. Circuits are derived and specified, with focused circuit optimization in multi-mode and multi-channel extensions, including multivariate potential energy terms and graph theoretic optimization from molecular symmetry. Resource estimation shows circuit depth advantage against QROM-loading architectures on a fault-tolerant scale, and a quantitative comparison against quantum signal processing variants confirms that a Trotter-based architecture retains a scalable T-gate advantage. Expanding beyond electronic states demonstrates that duality between finite basis and discrete variable representations permits congruent structural decompositions into quantum circuits, expanding the use of multi-channel dynamics far beyond chemistry.
Quantum computation is an attractive front for many problems that are intractable for computers today. One such problem is nonadiabatic quantum molecular dynamics, where quantized internal states coupling to parameterized modes result in a Hamiltonian resistant to oracle-based models and spectral decomposition. This dissertation applies diabatic Hamiltonian operators directly to the computational basis as first-quantized split-operator propagators, validated with dynamic observables including absorption and recurrence spectra, scattering cross-sections, population dynamics, and quantum scars. Circuits are derived and specified, with focused circuit optimization in multi-mode and multi-channel extensions, including multivariate potential energy terms and graph theoretic optimization from molecular symmetry. Resource estimation shows circuit depth advantage against QROM-loading architectures on a fault-tolerant scale, and a quantitative comparison against quantum signal processing variants confirms that a Trotter-based architecture retains a scalable T-gate advantage. Expanding beyond electronic states demonstrates that duality between finite basis and discrete variable representations permits congruent structural decompositions into quantum circuits, expanding the use of multi-channel dynamics far beyond chemistry.
Quantum advantage claims in chemical dynamics are difficult to evaluate when critical assumptions---model reductions, discretizations, Hamiltonian access models, numerical precision, and validation practices---are implicit.Focusing on discretized diabatic nonadiabatic Hamiltonians Hdia(Q)=TQ⊗IN+Vdia(Q), we propose a unified theoretical and methodological framework that standardizes:(i) qubit register encodings and physical-subspace embeddings including a leakage diagnostic when dimensions are not powers of two;(ii) Hamiltonian access models, i.e., explicit decompositions versus sparse-oracle access with fixed-point value semantics;(iii) instance-parameterized resource envelopes for Hamiltonian simulation and observable estimation;(iv) an additive error-budget ledger separating modeling/truncation error from oracle precision, unitary simulation error, estimation error, leakage, and hardware noise; and(v) a benchmark problem family Pnonadchem together with prototype instance ladders and a reporting checklist for matched classical--quantum comparisons.To support reproducibility without requiring software dependencies, we further specify a lightweight intermediate representation (IR) as a proposed data standard for nonadiabatic benchmark instances, including units, conventions, discretizations, embeddings, access declarations, and precision metadata.This paper is intentionally free of numerical results and implementation claims; it is a proposed baseline intended to make future algorithmic and experimental studies comparable and reviewable.
With quantum computers on the horizon, it is essential to understand what problems they can be used to solve. One potential application is simulating quantum systems, in particular, systems with many interacting electrons. However, there is a disparity between the computational resources required to study these problems and the capabilities of modern experimental quantum devices. In this thesis, we make progress toward closing this gap by reducing the cost of quantum simulation of quantum chemistry on a number of fronts. In the first half, we develop three improved quantum simulation algorithms by approximating time evolution under the Hamiltonian using a truncated Taylor series rather than the Trotter-Suzuki decomposition. We do this first for the molecular electronic structure Hamiltonian in two representations, obtaining asymptotic improvements in terms of both the number of orbitals and the required precision. We also apply the truncated Taylor series paradigm to quantum simulation in real space, where we consider simulation of general wave functions discretized on a uniform grid. In the second half, by working with Hamiltonians tailored specifically to quantum algorithms, we demonstrate that Trotter steps of the electronic structure Hamiltonian can be simulated in exactly N depth and with only linear connectivity, where N is the number of spin-orbitals. This resolves the problem of non-locality in fermion-to-qubit mappings, and has wide-ranging applications for variational and phase estimation-based quantum simulation of quantum chemistry. We further study this algorithm on an error-corrected quantum computer using the surface code. Benchmarking it alongside other product formula approaches, we show that quantum simulation of scientifically important instances of the electronic structure problem that are classically intractable can be performed with only a few hundred thousand physical qubits and fewer than one billion costly fault-tolerant gates. Finally, we study the effect of using randomized Hamiltonians in iterative phase estimation, and how this allows us to reduce the costs of simulation as well as to exploit knowledge of the ground state in how we design our simulation algorithm.
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.
The theoretical investigation of nonadiabatic processes is hampered by the complexity of the coupled electron-nuclear dynamics beyond the Born-Oppenheimer approximation. Classically, the simulation of such reactions is limited by the unfavorable scaling of the computational resources as a function of the system size. While quantum computing exhibits proven quantum advantage for the simulation of real-time dynamics, the study of quantum algorithms for the description of nonadiabatic phenomena is still unexplored. In this Letter, we propose a quantum algorithm for the simulation of fast nonadiabatic chemical processes together with an initialization scheme for quantum hardware calculations. In particular, we introduce a first-quantization method for the time evolution of a wave packet on two coupled harmonic potential energy surfaces (Marcus model). In our approach, the computational resources scale polynomially in the system dimensions, opening up new avenues for the study of photophysical processes that are classically intractable.
In quantum chemistry, the "electronic structure problem" refers to the process of numerically solving Schrodinger's equation for a given molecular system. These calculations rank among the most CPU-intensive computations that are performed nowadays, yet they are also crucial to the design of new drugs and materials. The idea behind quantum simulation of quantum chemistry is that a quantum computer could perform such calculations in a manner that is exponentially more efficient than what a classical computer could accomplish. In fact, quantum simulation is quite possibly the application most naturally suited to a quantum computer. In this work we introduce novel algorithms for the simulation of quantum chemistry, algorithms that are exponentially more precise than previous algorithms in the literature. Previous algorithms are based on the Trotter-Suzuki decomposition and scale like the inverse of the desired precision, while our algorithms represent the first practical application of the recently developed truncated Taylor series approach, and they achieve scaling that is logarithmic in the inverse of the precision. We explore algorithms for both second-quantized and first-quantized encodings of the chemistry Hamiltonian, where the second-quantized encoding requires mathcalO(N) qubits for an N spin-orbital system, and the highly compressed first-quantized encoding requires mathcalO(eta) qubits, where eta<<N is the number of electrons in the system. We first introduce a second-quantized algorithm that uses pre-computed molecular integrals, resulting in a gate count that scales like widetildecalO(N8t). Next we introduce algorithms that compute molecular integrals on the fly by carefully discretizing these integrals. Indeed, figuring out how to evaluate the integrals in an efficient manner is one of the main challenges of applying the truncated Taylor series approach to the chemistry problem: the complexity of a general quantum algorithm is traditionally formulated in terms of the number of queries made to an oracle whose assumption we just assume, but in this case, we actually construct the oracle for our specific problem. The evaluation of the oracle itself requires evaluating molecular integrals, which adds additional complexity to an algorithm. We find that the gate count of the second-quantized on-the-fly algorithm scales like widetildecalO(N5t). We also find that the first-quantized on-the-fly algorithm, which uses the configuration interaction (CI) matrix representation of the Hamiltonian, requires at most widetildecalO(eta2N3t) gates. In general, the discovery of exponentially more precise algorithms will allow for the simulation of larger molecular systems than was previously possible under Trotter-based methods. In addition, the results here can be generalized to the simulation of other fermionic systems beyond the electronic structure problem.
ConspectusSimulating molecular dynamics (MD) within a comprehensive quantum framework has been a long-standing challenge in computational chemistry. An exponential scaling of computational cost renders solving the time dependent Schrödinger equation (TDSE) of a molecular Hamiltonian, including both electronic and nuclear degrees of freedom (DOFs), as well as their couplings, infeasible for more than a few DOFs. In the Born-Oppenheimer (BO), or adiabatic, picture, electronic and nuclear parts of the wave function are decoupled and treated separately. Within this framework, the nuclear wave function evolves along potential energy surfaces (PESs) computed as solutions to the electronic Schrödinger equation parametrized in the nuclear DOFs. This approximation, together with increasingly elaborate numerical approaches to solve the nuclear time dependent Schrödinger equation (TDSE), enabled the treatment of up to a few dozens of degrees of freedom (DOFs). However, for particular applications, such as photochemistry, the BO approximation breaks down. In this regime of non-adiabatic dynamics, solving the full molecular problem including electron-nuclear couplings becomes essential, further increasing the complexity of the numerical solution. Although valuable methods such as multiconfigurational time-dependent Hartree (MCTDH) have been proposed for the solution of the coupled electron-nuclear dynamics, they remain hampered by an exponential scaling in the number of nuclear DOFs and by the difficulty of finding universal variational forms.In this Account, we present a perspective on novel quantum computational algorithms, aiming to alleviate the exponential scaling inherent to the simulation of many-body quantum dynamics. In particular, we focus on the derivation and application of quantum algorithms for adiabatic and non-adiabatic quantum dynamics, which include efficient approaches for the calculation of the BO potential energy surfaces (PESs). Thereafter, we study the time-evolution of a model system consisting of two coupled PESs in first and second quantization. In a first application, we discuss a recently introduced quantum algorithm for the evolution of a wavepacket in first quantization and exploit the potential quantum advantage of mapping its spatial grid representation to logarithmically many qubits. For the second demonstration, we move to the second quantization framework and review the scaling properties of two alternative time-evolution algorithms, namely, a variational quantum algorithm (VQA) (based on the McLachlan variational principle) and conventional Trotter-type evolution (based on a Lie-Trotter-Suzuki formula). Both methods clearly demonstrate the potential of quantum algorithms and their favorable scaling compared to the available classical approaches. However, a clear demonstration of quantum advantage in the context of molecular quantum dynamics may require the implementation of these algorithms in fault-tolerant quantum computers, while their application in near-term, noisy quantum devices is still unclear and deserves further investigation.