No Arabic abstract
We report the first electronic structure calculation performed on a quantum computer without exponentially costly precompilation. We use a programmable array of superconducting qubits to compute the energy surface of molecular hydrogen using two distinct quantum algorithms. First, we experimentally execute the unitary coupled cluster method using the variational quantum eigensolver. Our efficient implementation predicts the correct dissociation energy to within chemical accuracy of the numerically exact result. Second, we experimentally demonstrate the canonical quantum algorithm for chemistry, which consists of Trotterization and quantum phase estimation. We compare the experimental performance of these approaches to show clear evidence that the variational quantum eigensolver is robust to certain errors. This error tolerance inspires hope that variational quantum simulations of classically intractable molecules may be viable in the near future.
In this work, we present a linear optical implementation for analog quantum simulation of molecular vibronic spectra, incorporating the non-Condon scattering operation with a quadratically small truncation error. Thus far, analog and digital quantum algorithms for achieving quantum speedup have been suggested only in the Condon regime, which refers to a transition dipole moment that is independent of nuclear coordinates. For analog quantum optical simulation beyond the Condon regime (i.e., non-Condon transitions) the resulting non-unitary scattering operations must be handled appropriately in a linear optical network. In this paper, we consider the first and second-order Herzberg-Teller expansions of the transition dipole moment operator for the non-Condon effect, for implementation on linear optical quantum hardware. We believe the method opens a new way to approximate arbitrary non-unitary operations in analog and digital quantum simulations. We report in-silico simulations of the vibronic spectra for naphthalene, phenanthrene, and benzene to support our findings.
Ultrafast chemical reactions are difficult to simulate because they involve entangled, many-body wavefunctions whose computational complexity grows rapidly with molecular size. In photochemistry, the breakdown of the Born-Oppenheimer approximation further complicates the problem by entangling nuclear and electronic degrees of freedom. Here, we show that analog quantum simulators can efficiently simulate molecular dynamics using commonly available bosonic modes to represent molecular vibrations. Our approach can be implemented in any device with a qudit controllably coupled to bosonic oscillators and with quantum hardware resources that scale linearly with molecular size, and offers significant resource savings compared to digital quantum simulation algorithms. Advantages of our approach include a time resolution orders of magnitude better than ultrafast spectroscopy, the ability to simulate large molecules with limited hardware using a Suzuki-Trotter expansion, and the ability to implement realistic system-bath interactions with only one additional interaction per mode. Our approach can be implemented with current technology; e.g., the conical intersection in pyrazine can be simulated using a single trapped ion. Therefore, we expect our method will enable classically intractable chemical dynamics simulations in the near term.
We propose the use of 2-dimensional Penning trap arrays as a scalable platform for quantum simulation and quantum computing with trapped atomic ions. This approach involves placing arrays of micro-structured electrodes defining static electric quadrupole sites in a magnetic field, with single ions trapped at each site and coupled to neighbors via the Coulomb interaction. We solve for the normal modes of ion motion in such arrays, and derive a generalized multi-ion invariance theorem for stable motion even in the presence of trap imperfections. We use these techniques to investigate the feasibility of quantum simulation and quantum computation in fixed ion lattices. In homogeneous arrays, we show that sufficiently dense arrays are achievable, with axial, magnetron and cyclotron motions exhibiting inter-ion dipolar coupling with rates significantly higher than expected decoherence. With the addition of laser fields these can realize tunable-range interacting spin Hamiltonians. We also show how local control of potentials allows isolation of small numbers of ions in a fixed array and can be used to implement high fidelity gates. The use of static trapping fields means that our approach is not limited by power requirements as system size increases, removing a major challenge for scaling which is present in standard radio-frequency traps. Thus the architecture and methods provided here appear to open a path for trapped-ion quantum computing to reach fault-tolerant scale devices.
Methods for electronic structure based on Gaussian and molecular orbital discretizations offer a well established, compact representation that forms much of the foundation of correlated quantum chemistry calculations on both classical and quantum computers. Despite their ability to describe essential physics with relatively few basis functions, these representations can suffer from a quartic growth of the number of integrals. Recent results have shown that, for some quantum and classical algorithms, moving to representations with diagonal two-body operators can result in dramatically lower asymptotic costs, even if the number of functions required increases significantly. We introduce a way to interpolate between the two regimes in a systematic and controllable manner, such that the number of functions is minimized while maintaining a block diagonal structure of the two-body operator and desirable properties of an original, primitive basis. Techniques are analyzed for leveraging the structure of this new representation on quantum computers. Empirical results for hydrogen chains suggest a scaling improvement from $O(N^{4.5})$ in molecular orbital representations to $O(N^{2.6})$ in our representation for quantum evolution in a fault-tolerant setting, and exhibit a constant factor crossover at 15 to 20 atoms. Moreover, we test these methods using modern density matrix renormalization group methods classically, and achieve excellent accuracy with respect to the complete basis set limit with a speedup of 1-2 orders of magnitude with respect to using the primitive or Gaussian basis sets alone. These results suggest our representation provides significant cost reductions while maintaining accuracy relative to molecular orbital or strictly diagonal approaches for modest-sized systems in both classical and quantum computation for correlated systems.
Electron transport in realistic physical and chemical systems often involves the non-trivial exchange of energy with a large environment, requiring the definition and treatment of open quantum systems. Because the time evolution of an open quantum system employs a non-unitary operator, the simulation of open quantum systems presents a challenge for universal quantum computers constructed from only unitary operators or gates. Here we present a general algorithm for implementing the action of any non-unitary operator on an arbitrary state on a quantum device. We show that any quantum operator can be exactly decomposed as a linear combination of at most four unitary operators. We demonstrate this method on a two-level system in both zero and finite temperature amplitude damping channels. The results are in agreement with classical calculations, showing promise in simulating non-unitary operations on intermediate-term and future quantum devices.