Search

Most quantum algorithms designed to generate or probe properties of the ground state of a quantum many-body system require as input an initial state with a large overlap with the desired ground state. One approach for preparing such a ground state is Imaginary Time Evolution (ITE). Recent work by [Motta, M., Sun, C., Tan, A.T.K. et al. (2020)] introduced an algorithm -- which we will refer to as Quantum Imaginary Time Evolution (QITE) -- that shows how ITE can be approximated by a sequence of unitary operators, making QITE potentially implementable on early fault-tolerant quantum computers. In this work, we provide a heuristic study of the capabilities of the QITE algorithm in approximating the ITE of lattice and molecular electronic structure Hamiltonians. We numerically study the performance of the QITE algorithm when provided with a good classical initial state for a large class of systems, some of which are of interest to industrial applications, and check if QITE is able to qualitatively replicate the ITE behavior and improve over a classical mean-field solution. The systems we consider in this work range from one- and two-dimensional lattice systems of various lattice geometries displaying short- and long-range interactions, to active spaces of molecular electronic structure Hamiltonians. In addition to the comparison of QITE and ITE, we explicitly show how imaginary time evolved fermionic Gaussian states can serve as initial states which can be efficiently computed on classical computers and efficiently implemented on quantum computers for generic spin Hamiltonians in arbitrary lattice geometries and dimensions, which can be of independent interest.

Quantum chemical calculations are among the most promising applications for quantum computing. Implementations of dedicated quantum algorithms on available quantum hardware were so far, however, mostly limited to comparatively simple systems without strong correlations. As such, they can also be addressed by classically efficient single-reference methods. In this work, we calculate the lowest energy eigenvalue of active space Hamiltonians of industrially relevant and strongly correlated metal chelates on trapped ion quantum hardware, and integrate the results into a typical industrial quantum chemical workflow to arrive at chemically meaningful properties. We are able to achieve chemical accuracy by training a variational quantum algorithm on quantum hardware, followed by a classical diagonalization in the subspace of states measured as outputs of the quantum circuit. This approach is particularly measurement-efficient, requiring 600 single-shot measurements per cost function evaluation on a ten qubit system, and allows for efficient post-processing to handle erroneous runs., Comment: 12 pages, 5 figures

Charge-carrier transport at the semiconductor-gate dielectric interface in organic field-effect transistors is critically dependent on the degree of disorder in the typically semi-crystalline semiconductor layer. Charge trapping can occur at the interface as well as in the current-carrying semiconductor layer itself. A detailed and systematic understanding of the role of grain boundaries between crystallites and how to avoid their potentially detrimental effects is still an important focus of research in the organic electronics community. A typical macroscopic measurement technique to extract information about the energetics of the grain boundaries is an activation energy measurement. Here, we compare detailed experiments on the energetic properties of monolayer thin films implemented in organic field-effect transistors, having controlled numbers of grain boundaries within the channel region to kinetic Monte-Carlo simulations of charge-carrier transport to elucidate the influence of grain boundaries on the extracted activation energies. Two important findings are: 1) whereas the energy at the grain boundary does not change with the number of grain boundaries in a thin film, both the measured and simulated activation energy increases with the number of grain boundaries. 2) In simulations where both energy barriers and valleys are present at the grain boundaries there is no systematic relation between the number of grain boundaries and extracted activation energies. We conclude, that a macroscopic measurement of the activation energy can serve as general quality indicator of the thin film, but does not allow microscopic conclusions about the energy landscape of the thin film.

The simulation of chemistry is among the most promising applications of quantum computing. However, most prior work exploring algorithms for block-encoding, time-evolving, and sampling in the eigenbasis of electronic structure Hamiltonians has either focused on modeling finite-sized systems, or has required a large number of plane wave basis functions. In this work, we extend methods for quantum simulation with Bloch orbitals constructed from symmetry-adapted atom-centered orbitals so that one can model periodic \textit{ab initio} Hamiltonians using only a modest number of basis functions. We focus on adapting existing algorithms based on combining qubitization with tensor factorizations of the Coulomb operator. Significant modifications of those algorithms are required to obtain an asymptotic speedup leveraging translational (or, more broadly, Abelian) symmetries. We implement block encodings using known tensor factorizations and a new Bloch orbital form of tensor hypercontraction. Finally, we estimate the resources required to deploy our algorithms to classically challenging model materials relevant to the chemistry of Lithium Nickel Oxide battery cathodes within the surface code.

In this work, we test a recently developed method to enhance classical auxiliary-field quantum Monte Carlo (AFQMC) calculations with quantum computers against examples from chemistry and material science, representatives of classes of industry-relevant systems. As molecular test cases, we calculate the energy curve of H4 and relative energies of ozone and singlet molecular oxygen with respect to triplet molecular oxygen, which are industrially relevant in organic oxidation reactions. We find that trial wave functions beyond single Slater determinants improve the performance of AFQMC and allow to generate energies close to chemical accuracy compared to full configuration interaction (FCI) or experimental results. As a representative for material science we study a quasi-1D Fermi-Hubbard model derived from CuBr2, a compound displaying electronic structure properties analogous to cuprates. We find that trial wave functions with both, significantly larger fidelities and lower energies over a Hartree-Fock solution, do not necessarily lead to better AFQMC results.

Spontaneous emission of quantum emitters can be modified by engineering their optical environment. This allows a resonant nanoantenna to significantly modify the radiative properties of a quantum emitter. In this article, we go beyond the common electric dipole approximation for the molecular electronic transition and take light-matter coupling through higher order multipoles into account. We investigate, by means of theory and numerical simulations, a strong enhancement of the magnetic dipole and electric quadrupole emission channels of a molecule adjacent to a plasmonic patch nanoantenna. While this on its own had been considered, the assumption in prior work usually has been that each molecular transition is dominated only by one of those multipolar emission channels. This leads naturally to the notion of discussing the modified emission in terms of a modified local density of states defined for each specific multipolar transition. In reality, this restricts the applicability of the approach, since specific molecular transitions occur via multiple multipolar pathways that have to be considered all at once. Here, we introduce a framework to study interference effects between higher order transitions in molecules by (a) a rigorous quantum-chemical calculation of their multipolar moments and (b) by a consecutive investigation of the transition rate upon coupling to an arbitrarily shaped nanoantenna. Based on that formalism we predict interference effects between these transition channels. This allows for a strong suppression of radiation by exploiting destructive interference. Our work suggests that placing a suitably chosen molecule at a well defined position and at a well defined orientation relative to a nanoantenna can fully suppress the transition probability., Comment: 30 pages, 8 figures

The study and prediction of chemical reactivity is one of the most important application areas of molecular quantum chemistry. Large-scale, fully error-tolerant quantum computers could provide exact or near-exact solutions to the underlying electronic structure problem with exponentially less effort than a classical computer thus enabling highly accurate predictions for comparably large molecular systems. In the nearer future, however, only "noisy" devices with a limited number of qubits that are subject to decoherence will be available. For such near-term quantum computers the hybrid quantum-classical variational quantum eigensolver algorithm in combination with the unitary coupled-cluster ansatz (UCCSD-VQE) has become an intensively discussed approach that could provide accurate results before the dawn of error-tolerant quantum computing. In this work we present an implementation of UCCSD-VQE that allows for the first time to treat both open- and closed-shell molecules. We study the accuracy of the obtained energies for nine small molecular systems as well as for four exemplary chemical reactions by comparing to well-established electronic structure methods like (non-unitary) coupled-cluster and density functional theory. Finally, we roughly estimate the required quantum hardware resources to obtain "useful" results for practical purposes., Comment: 26 pages, 15 figures

In this article we consider a special type of degenerate elliptic partial differential equations of second order in convex domains that satisfy the interior sphere condition. We show that any positive viscosity solution $u$ of $-|\nabla u|^\alpha \Delta_p^N u = 1$ has the property that $u^\frac{\alpha + 1}{\alpha + 2}$ is a concave function. Secondly we consider positive solutions of the eigenvalue problem $-|\nabla u|^\alpha \Delta_p^N u = \lambda |u|^{\alpha} u$, in which case $\log u$ turns out to be concave. The methods provided include a weak comparison principle and a Hopf-type Lemma.