Efficient quantum algorithm for computing n-time correlation functions

Simulating Non-Markovian Dynamics in Multidimensional Electronic Spectroscopy via Quantum Algorithm

Including the effect of the molecular environment in the numerical modeling of time-resolved electronic spectroscopy remains an important challenge in computational spectroscopy. In this contribution, we present a general approach for the simulation of the optical response of multichromophore systems in a structured environment and its implementation as a quantum algorithm. A key step of the procedure is the pseudomode embedding of the system-environment problem resulting in a finite set of quantum states evolving according to a Markovian quantum master equation. This formulation is then solved by a collision model integrated into a quantum algorithm designed to simulate linear and nonlinear response functions. The workflow is validated by simulating spectra for the prototypical excitonic dimer interacting with fast (memoryless) and finite-memory environments. The results demonstrate, on the one hand, the potential of the pseudomode embedding for simulating the dynamical features of nonlinear spectroscopy, including lineshape, spectral diffusion, and relaxations along delay times. On the other hand, the explicit synthesis of quantum circuits provides a fully quantum simulation protocol of nonlinear spectroscopy harnessing the efficient quantum simulation of many-body dynamics promised by the future generation of fault-tolerant quantum computers.

Hybrid Hamiltonian Simulation for Excitation Dynamics

Hamiltonian simulation is one of the most anticipated applications of quantum computing. Quantum circuit depth for implementing Hamiltonian simulation is commonly time dependent using Trotter-Suzuki product formulas so that long time quantum dynamic simulations (QDSs) become impratical for near-term quantum processors. Hamiltonian simulation based on Cartan decomposition (CD) provides an appealing scheme for QDSs with fixed-depth circuits, while it is limited to a time-independent Hamiltonian. In this work, we generalize this CD-based Hamiltonian simulation algorithm for studying time-dependent systems by combining it with variational quantum algorithms. The time-dependent and time-independent parts of the Hamiltonian are treated by using variational and CD-based Hamiltonian simulation algorithms, respectively. As such, this hybrid Hamiltonian simulation requires only fixed-depth quantum circuits to handle time-dependent cases while maintaining a high accuracy. We apply this new algorithm to study the response of spin and molecular systems to δ-kick electric fields and obtain accurate spectra for these excitation processes.

Snapshotting quantum dynamics at multiple time points

Measurement-induced state disturbance is a major challenge in obtaining quantum statistics at multiple time points. We propose a method to extract dynamic information from a quantum system at intermediate time points, namely snapshotting quantum dynamics. To this end, we apply classical post-processing after performing the ancilla-assisted measurements to cancel out the impact of the measurements at each time point. Based on this, we reconstruct a multi-time quasi-probability distribution (QPD) that correctly recovers the probability distributions at the respective time points. Our approach can also be applied to simultaneously extract exponentially many correlation functions with various time-orderings. We provide a proof-of-principle experimental demonstration of the proposed protocol using a dual-species trapped-ion system by employing 171Yb+ and 138Ba+ ions as the system and the ancilla, respectively. Multi-time measurements are performed by repeated initialization and detection of the ancilla state without directly measuring the system state. The two- and three-time QPDs and correlation functions are reconstructed reliably from the experiment, negativity and complex values in the QPDs clearly indicate a contribution of the quantum coherence throughout dynamics.

Adaptive Variational Quantum Computing Approaches for Green's Functions and Nonlinear Susceptibilities

We present and benchmark quantum computing approaches for calculating real-time single-particle Green's functions and nonlinear susceptibilities of Hamiltonian systems. The approaches leverage adaptive variational quantum algorithms for state preparation and propagation. Using automatically generated compact circuits, the dynamical evolution is performed over sufficiently long times to achieve adequate frequency resolution of the response functions. We showcase accurate Green's function calculations using a statevector simulator on classical hardware for Fermi-Hubbard chains of 4 and 6 sites, with maximal ansatz circuit depths of 65 and 424 layers, respectively, and for the molecule LiH with a maximal ansatz circuit depth of 81 layers. Additionally, we consider an antiferromagnetic quantum spin-1 model that incorporates the Dzyaloshinskii-Moriya interaction to illustrate calculations of the third-order nonlinear susceptibilities, which can be measured in two-dimensional coherent spectroscopy experiments. These results demonstrate that real-time approaches using adaptive parametrized circuits to evaluate linear and nonlinear response functions can be feasible with near-term quantum processors.

A Quantum Algorithm from Response Theory: Digital Quantum Simulation of Two-Dimensional Electronic Spectroscopy

Multidimensional optical spectroscopies are powerful techniques to investigate energy transfer pathways in natural and artificial systems. Because of the high information content of the spectra, numerical simulations of the optical response are of primary importance to assist the interpretation of spectral features. However, the increasing complexity of the investigated systems and their quantum dynamics call for the development of novel simulation strategies. In this work, we consider using digital quantum computers. By combining quantum dynamical simulation and nonlinear response theory, we present a quantum algorithm for computing the optical response of molecular systems. The quantum advantage stems from the efficient quantum simulation of the dynamics governed by the molecular Hamiltonian, and it is demonstrated by explicitly considering exciton-vibrational coupling. The protocol is tested on a near-term quantum device, providing the digital quantum simulation of the linear and nonlinear response of simple molecular models.

Predicting molecular vibronic spectra using time-domain analog quantum simulation

Spectroscopy is one of the most accurate probes of the molecular world. However, predicting molecular spectra accurately is computationally difficult because of the presence of entanglement between electronic and nuclear degrees of freedom. Although quantum computers promise to reduce this computational cost, existing quantum approaches rely on combining signals from individual eigenstates, an approach whose cost grows exponentially with molecule size. Here, we introduce a method for scalable analog quantum simulation of molecular spectroscopy: by performing simulations in the time domain, the number of required measurements depends on the desired spectral range and resolution, not molecular size. Our approach can treat more complicated molecular models than previous ones, requires fewer approximations, and can be extended to open quantum systems with minimal overhead. We present a direct mapping of the underlying problem of time-domain simulation of molecular spectra to the degrees of freedom and control fields available in a trapped-ion quantum simulator. We experimentally demonstrate our algorithm on a trapped-ion device, exploiting both intrinsic electronic and motional degrees of freedom, showing excellent quantitative agreement for a single-mode vibronic photoelectron spectrum of SO2.

Simulating Hydrodynamics on Noisy Intermediate-Scale Quantum Devices with Random Circuits

In a recent milestone experiment, Google's processor Sycamore heralded the era of "quantum supremacy" by sampling from the output of (pseudo-)random circuits. We show that such random circuits provide tailor-made building blocks for simulating quantum many-body systems on noisy intermediate-scale quantum (NISQ) devices. Specifically, we propose an algorithm consisting of a random circuit followed by a trotterized Hamiltonian time evolution to study hydrodynamics and to extract transport coefficients in the linear response regime. We numerically demonstrate the algorithm by simulating the buildup of spatiotemporal correlation functions in one- and two-dimensional quantum spin systems, where we particularly scrutinize the inevitable impact of errors present in any realistic implementation. Importantly, we find that the hydrodynamic scaling of the correlations is highly robust with respect to the size of the Trotter step, which opens the door to reach nontrivial time scales with a small number of gates. While errors within the random circuit are shown to be irrelevant, we furthermore unveil that meaningful results can be obtained for noisy time evolutions with error rates achievable on near-term hardware. Our work emphasizes the practical relevance of random circuits on NISQ devices beyond the abstract sampling task.

Simulation of molecular spectroscopy with circuit quantum electrodynamics

Spectroscopy is a crucial laboratory technique for understanding quantum systems through their interactions with the electromagnetic radiation. Particularly, spectroscopy is capable of revealing the physical structure of molecules, leading to the development of the maser-the forerunner of the laser. However, real-world applications of molecular spectroscopy are mostly confined to equilibrium states, due to computational and technological constraints; a potential breakthrough can be achieved by utilizing the emerging technology of quantum simulation. Here we experimentally demonstrate through a toy model, a superconducting quantum simulator capable of generating molecular spectra for both equilibrium and non-equilibrium states, reliably producing the vibronic structure of diatomic molecules. Furthermore, our quantum simulator is applicable not only to molecules with a wide range of electronic-vibronic coupling strength, characterized by the Huang-Rhys parameter, but also to molecular spectra not readily accessible under normal laboratory conditions. These results point to a new direction for predicting and understanding molecular spectroscopy, exploiting the power of quantum simulation.

Measurement of linear response functions in Nuclear Magnetic Resonance

We measure multi-time correlation functions of a set of Pauli operators on a two-level system, which can be used to retrieve its associated linear response functions. The two-level system is an effective spin constructed from the nuclear spins of ¹H atoms in a solution of ¹³C-labeled chloroform. Response functions characterize the linear response of the system to a family of perturbations, allowing us to compute physical quantities such as the magnetic susceptibility of the effective spin. We use techniques exported from quantum information to measure time correlations on the two-level system. This approach requires the use of an ancillary qubit encoded in the nuclear spins of the ¹³C atoms and a sequence of controlled operations. Moreover, we demonstrate the ability of such a quantum platform to compute time-correlation functions of arbitrary order, which relate to higher-order corrections of perturbative methods. Particularly, we show three-time correlation functions for arbitrary times, and we also measure time correlation functions at fixed times up to tenth order.

Digital Quantum Simulation of Minimal AdS/CFT

We propose the digital quantum simulation of a minimal AdS/CFT model in controllable quantum platforms. We consider the Sachdev-Ye-Kitaev model describing interacting Majorana fermions with randomly distributed all-to-all couplings, encoding nonlocal fermionic operators onto qubits to efficiently implement their dynamics via digital techniques. Moreover, we also give a method for probing nonequilibrium dynamics and the scrambling of information. Finally, our approach serves as a protocol for reproducing a simplified low-dimensional model of quantum gravity in advanced quantum platforms as trapped ions and superconducting circuits.

Quantum Simulation of Generic Many-Body Open System Dynamics Using Classical Noise

We introduce a scheme for the quantum simulation of many-body decoherence based on the unitary evolution of a stochastic Hamiltonian. Modulating the strength of the interactions with stochastic processes, we show that the noise-averaged density matrix simulates an effectively open dynamics governed by k-body Lindblad operators. Markovian dynamics can be accessed with white-noise fluctuations; non-Markovian dynamics requires colored noise. The time scale governing the fidelity decay under many-body decoherence is shown to scale as N^{-2k} with the system size N. Our proposal can be readily implemented in a variety of quantum platforms including optical lattices, superconducting circuits, and trapped ions.

Nonergodicity in the Anisotropic Dicke Model

We study the ergodic-nonergodic transition in a generalized Dicke model with independent corotating and counterrotating light-matter coupling terms. By studying level statistics, the average ratio of consecutive level spacings, and the quantum butterfly effect (out-of-time correlation) as a dynamical probe, we show that the ergodic-nonergodic transition in the Dicke model is a consequence of the proximity to the integrable limit of the model when one of the couplings is set to zero. This can be interpreted as a hint for the existence of a quantum analogue of the classical Kolmogorov-Arnold-Moser theorem. In addition, we show that there is no intrinsic relation between the ergodic-nonergodic transition and the precursors of the normal-superradiant quantum phase transition.

Quantum Simulation of Dissipative Processes without Reservoir Engineering

We present a quantum algorithm to simulate general finite dimensional Lindblad master equations without the requirement of engineering the system-environment interactions. The proposed method is able to simulate both Markovian and non-Markovian quantum dynamics. It consists in the quantum computation of the dissipative corrections to the unitary evolution of the system of interest, via the reconstruction of the response functions associated with the Lindblad operators. Our approach is equally applicable to dynamics generated by effectively non-Hermitian Hamiltonians. We confirm the quality of our method providing specific error bounds that quantify its accuracy.