BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

David A. Mazziotti

David A. Mazziotti appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.chem-phquant-phphysics.comp-phcond-mat.mes-hallcond-mat.mtrl-scicond-mat.str-elmath-phmath.MPBiological Physicscond-mat.supr-conphysics.atom-phphysics.optics

Catalog footprint

What is connected

42works
12topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

42 published item(s)

preprint2022arXiv

A general quantum algorithm for open quantum dynamics demonstrated with the Fenna-Matthews-Olson complex

Using quantum algorithms to simulate complex physical processes and correlations in quantum matter has been a major direction of quantum computing research, towards the promise of a quantum advantage over classical approaches. In this work we develop a generalized quantum algorithm to simulate any dynamical process represented by either the operator sum representation or the Lindblad master equation. We then demonstrate the quantum algorithm by simulating the dynamics of the Fenna-Matthews-Olson (FMO) complex on the IBM QASM quantum simulator. This work represents a first demonstration of a quantum algorithm for open quantum dynamics with a moderately sophisticated dynamical process involving a realistic biological structure. We discuss the complexity of the quantum algorithm relative to the classical method for the same purpose, presenting a decisive query complexity advantage of the quantum approach based on the unique property of quantum measurement.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Accelerated Convergence of Contracted Quantum Eigensolvers through a Quasi-Second-Order, Locally Parameterized Optimization

A contracted quantum eigensolver (CQE) finds a solution to the many-electron Schrödinger equation by solving its integration (or contraction) to the 2-electron space -- a contracted Schrödinger equation (CSE) -- on a quantum computer. When applied to the anti-Hermitian part of the CSE (ACSE), the CQE iterations optimize the wave function with respect to a general product ansatz of two-body exponential unitary transformations that can exactly solve the Schrödinger equation. In this work, we accelerate the convergence of the CQE and its wavefunction ansatz via tools from classical optimization theory. By treating the CQE algorithm as an optimization in a local parameter space, we can apply quasi-second-order optimization techniques, such as quasi-Newton approaches or non-linear conjugate gradient approaches. Practically these algorithms result in superlinear convergence of the wavefunction to a solution of the ACSE. Convergence acceleration is important because it can both minimize the accumulation of noise on near-term intermediate-scale quantum (NISQ) computers and achieve highly accurate solutions on future fault-tolerant quantum devices. We demonstrate the algorithm, as well as some heuristic implementations relevant for cost-reduction considerations, comparisons with other common methods such as variational quantum eigensolvers, and a fermionic-encoding-free form of the CQE.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Beginnings of Exciton Condensation in Coronene Analog of Graphene Double Layer

Exciton condensation, a Bose-Einstein condensation of excitons into a single quantum state, has recently been achieved in low-dimensional materials including twin layers of graphene and van der Waals heterostructures. Here we examine computationally the beginnings of exciton condensation in a double layer comprised of coronene, a seven-benzene-ring patch of graphene. As a function of interlayer separation, we compute the exciton population in a single coherent quantum state, showing that the population peaks around 1.8 at distances near 2 Å. Visualization reveals interlayer excitons at the separation distance of the condensate. We determine the exciton population as a function of the twist angle between the two coronene layers to reveal the magic angles at which the condensation peaks. As with previous recent calculations showing some exciton condensation in hexacene double layers and benzene stacks, the present two-electron reduced-density-matrix calculations with coronene provide computational evidence for the ability to realize exciton condensation in molecular-scale analogs of extended systems like the graphene double layer.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Comparison of Density-Matrix Corrections to Density Functional Theory

Density functional theory (DFT), one of the most widely utilized methods available to computational chemistry, fails to describe systems with statically correlated electrons. To address this shortcoming, in previous work we transformed DFT into a one-electron reduced density matrix theory (1-RDMFT) via the inclusion of a quadratic one-electron reduced density matrix (1-RDM) correction. Here, we combine our 1-RDMFT approach with different DFT functionals as well as Hartree-Fock to elucidate the method's dependence on the underlying functional selection. Furthermore, we generalize the information density matrix functional theory (iDMFT), recently developed as a correction to the Hartree-Fock method, by incorporating density functionals in place of the Hartree-Fock functional. We relate iDMFT mathematically to our approach and benchmark the two with a common set of functionals and systems.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Cooper-Pair Condensates with Non-Classical Long-Range Order on Quantum Devices

An important problem in quantum information is the practical demonstration of non-classical long-range order on quantum computers. One of the best known examples of a quantum system with non-classical long-range order is a superconductor. Here we achieve Cooper pairing of qubits on a quantum computer to represent superconducting or superfluid states. We rigorously confirm the quantum long-range order by measuring the large $O(N)$ eigenvalue of the two-electron reduced density matrix. The demonstration of maximal quantum long-range order is an important step towards more complex modeling of superconductivity and superfluidity as well as other phenomena with significant quantum long-range order on quantum computers.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Density Functional Theory Transformed into a One-electron Reduced Density Matrix Functional Theory for the Capture of Static Correlation

Density functional theory (DFT), the most widely adopted method in modern computational chemistry, fails to describe accurately the electronic structure of strongly correlated systems. Here we show that DFT can be formally and practically transformed into a one-electron reduced-density-matrix (1-RDM) functional theory, which can address the limitations of DFT while retaining favorable computational scaling compared to wavefunction-based approaches. In addition to relaxing the idempotency restriction on the 1-RDM in the kinetic energy term, we add a quadratic 1-RDM-based term to DFT's density-based exchange-correlation functional. Our approach, which we implement by quadratic semidefinite programming at DFT's computational scaling of $O(r^{3})$, yields substantial improvements over traditional DFT in the description of static correlation in chemical structures and processes such as singlet biradicals and bond dissociations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Elucidating the molecular orbital dependence of the total electronic energy in multireference problems

The accurate resolution of the chemical properties of strongly correlated systems, such as biradicals, requires the use of electronic structure theories that account for both multi-reference as well as dynamic correlation effects. A variety of methods exist that aim to resolve the dynamic correlation in multi-reference problems, commonly relying on an exponentially scaling complete-active-space self-consistent-field (CASSCF) calculation to generate reference molecular orbitals (MOs). However, while CASSCF orbitals provide the optimal solution for a selected set of correlated (active) orbitals, their suitability in the quest for the resolution of the total correlation energy has not been thoroughly investigated. Recent research has shown the ability of Kohn-Shan density functional theory (KS-DFT) to provide improved orbitals for coupled cluster (CC) and Møller-Plesset perturbation theory (MP) calculations. Here we extend the search for optimal and more cost effective MOs to post-configuration-interaction (post-CI) methods, surveying the ability of the MOs obtained with various DFT functionals, as well as Hartree-Fock, and CC and MP calculations to accurately capture the total electronic correlation energy. Applying the anti-Hermitian contracted Schrödinger equation (ACSE) to the dissociation of N$_2$, the calculation of biradical singlet-triplet gaps and the transition states of the bicylobutane isomerization, we demonstrate DFT provides a cost-effective alternative to CASSCF in providing reference orbitals for post-CI dynamic correlation calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Entangled phase of simultaneous fermion and exciton condensations realized

Fermion-exciton condensates (FECs) -- computationally and theoretically predicted states that simultaneously exhibit the character of superconducting states and exciton condensates -- are novel quantum states whose properties may involve a hybridization of superconductivity and the dissipationless flow of energy. Here, we exploit prior investigations of superconducting states and exciton condensates on quantum devices to identify a tuneable quantum state preparation entangling the wave functions of the individual condensate states. Utilizing this state preparation, we prepare a variety of FEC states on quantum computers -- realizing strongly correlated FEC states on current, noisy intermediate-scale quantum devices -- and verify the presence of the dual condensate via postmeasurement analysis. This confirmation of the previously predicted condensate state on quantum devices as well as the form of its wave function motivates further theoretical and experimental exploration of the properties, applications, and stability of FECs.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Interplay of Electronic and Geometric Structure Tunes Organic Biradical Character in Bimetallic Tetrathiafulvalene Tetrathiolate Complexes

The synthesis and design of organic biradicals with tunable singlet-triplet gaps has become the subject of significant research interest, owing to their possible photochemical applications and use in the development of molecular switches and conductors. Recently, tetrathiafulvalene tetrathiolate (TTFtt) has been demonstrated to exhibit such organic biradical character in doubly ionized bimetallic complexes. In this article we use high-level {\em ab initio} calculations to interrogate the electronic structure of a series of TTFtt-bridged metal complexes, resolving the factors governing their biradical character and singlet-triplet gaps. We show that the degree of biradical character correlates with a readily measured experimental predictor, the central TTFtt C-C bond length, and that it may be described by a one-parameter model, providing valuable insight for the future rational design of TTFtt based biradical compounds and materials.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Large Cumulant Eigenvalue as a Signature of Exciton Condensation

The Bose-Einstein condensation of excitons into a single quantum state is known as exciton condensation. Exciton condensation, which potentially supports the frictionless flow of energy, has recently been realized in graphene bilayers and van der Waals heterostructures. Here we show that exciton condensates can be predicted from a combination of reduced density matrix theory and cumulant theory. We show that exciton condensation occurs if and only if there exists a large eigenvalue in the cumulant part of the particle-hole reduced density matrix. In the thermodynamic limit we show that the large eigenvalue is bounded from above by the number of excitons. In contrast to the eigenvalues of the particle-hole matrix, the large eigenvalue of the cumulant matrix has the advantage of providing a size-extensive measure of the extent of condensation. Here we apply this signature to predict exciton condensation in both the Lipkin model and molecular stacks of benzene. The computational signature has applications to the prediction of exciton condensation in both molecules and materials.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Many-Fermion Simulation from the Contracted Quantum Eigensolver without Fermionic Encoding of the Wave Function

Quantum computers potentially have an exponential advantage over classical computers for the quantum simulation of many-fermion quantum systems. Nonetheless, fermions are more expensive to simulate than bosons due to the fermionic encoding -- a mapping by which the qubits are encoded with fermion statistics. Here we generalize the contracted quantum eigensolver (CQE) to avoid fermionic encoding of the wave function. In contrast to the variational quantum eigensolver, the CQE solves for a many-fermion stationary state by minimizing the contraction (projection) of the Schrödinger equation onto two fermions. We avoid fermionic encoding of the wave function by contracting the Schrödinger equation onto an unencoded pair of particles. Solution of the resulting contracted equation by a series of unencoded two-body exponential transformations generates an unencoded wave function from which the energy and two-fermion reduced density matrix (2-RDM) can be computed. We apply the unencoded and the encoded CQE algorithms to the hydrogen fluoride molecule, the dissociation of oxygen O$_{2}$, and a series of hydrogen chains. Both algorithms show comparable convergence towards the exact ground-state energies and 2-RDMs, but the unencoded algorithm has computational advantages in terms of state preparation and tomography.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Quantum Phase Transitions in a Model Hamiltonian Exhibiting Entangled Simultaneous Fermion-Pair and Exciton Condensations

Quantum states of a novel Bose-Einstein condensate, in which both fermion-pair and exciton condensations are simultaneously present, have recently been realized theoretically in a model Hamiltonian system. Here we identify quantum phase transitions in that model between fermion-pair and exciton condensations based on a geometric analysis of the convex set of ground-state 2-particle reduced density matrices (2-RDMs). The 2-RDM set provides a finite representation of the infinite parameter space of Hamiltonians that readily reveals a fermion-pair condensate phase and two distinct exciton condensate phases, as well as the emergence of first- and second-order phase transitions as the particle number of the system is increased. The set, furthermore, shows that the fermion-exciton condensate (FEC) lies along the second-order phase transition between the exciton and fermion-pair condensate phases. The detailed information about the exciton and fermion-pair phases, the forces behind these phase, as well as their associated transitions provides additional insight into the formation of the FEC condensate, which we anticipate will prove useful in its experimental realization.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Quantum Simulation of Open Quantum Systems Using Density-Matrix Purification

Electronic structure and transport in realistically-sized systems often require an open quantum system (OQS) treatment, where the system is defined in the context of an environment. As OQS evolution is non-unitary, implementation on quantum computers -- limited to unitary operations -- is challenging. We present a general framework for OQSs where the system's $d \times d$ density matrix is recast as a $d^{2}$ wavefunction which can be evolved by unitary transformations. This theory has two significant advantages over conventional approaches: (i) the wavefunction requires only an $n$-qubit, compared to $2n$-qubit, bath for an $n$-qubit system and (ii) the purification includes dynamics of any pure-state universe. We demonstrate this method on a two-level system in a zero temperature amplitude damping channel and a two-site quantum Ising model. Quantum simulation and experimental-device results agree with classical calculations, showing promise in simulating non-unitary operations on NISQ quantum devices.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Quantum Simulation of Quantum Phase Transitions Using the Convex Geometry of Reduced Density Matrices

Transitions of many-particle quantum systems between distinct phases at absolute-zero temperature, known as quantum phase transitions, require an exacting treatment of particle correlations. In this work, we present a general quantum-computing approach to quantum phase transitions that exploits the geometric structure of reduced density matrices. While typical approaches to quantum phase transitions examine discontinuities in the order parameters, the origin of phase transitions -- their order parameters and symmetry breaking -- can be understood geometrically in terms of the set of two-particle reduced density matrices (2-RDMs). The convex set of 2-RDMs provides a comprehensive map of the quantum system including its distinct phases as well as the transitions connecting these phases. Because 2-RDMs can potentially be computed on quantum computers at non-exponential cost, even when the quantum system is strongly correlated, they are ideally suited for a quantum-computing approach to quantum phase transitions. We compute the convex set of 2-RDMs for a Lipkin-Meshkov-Glick spin model on IBM superconducting-qubit quantum processors. Even though computations are limited to few-particle models due to device noise, comparisons with a classically solvable 1000-particle model reveal that the finite-particle quantum solutions capture the key features of the phase transitions including the strong correlation and the symmetry breaking.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Quantum State Preparation and Non-Unitary Evolution with Diagonal Operators

Realizing non-unitary transformations on unitary-gate based quantum devices is critically important for simulating a variety of physical problems including open quantum systems and subnormalized quantum states. We present a dilation based algorithm to simulate non-unitary operations using probabilistic quantum computing with only one ancilla qubit. We utilize the singular-value decomposition (SVD) to decompose any general quantum operator into a product of two unitary operators and a diagonal non-unitary operator, which we show can be implemented by a diagonal unitary operator in a 1-qubit dilated space. While dilation techniques increase the number of qubits in the calculation, and thus the gate complexity, our algorithm limits the operations required in the dilated space to a diagonal unitary operator, which has known circuit decompositions. We use this algorithm to prepare random sub-normalized two-level states on a quantum device with high fidelity. Furthermore, we present the accurate non-unitary dynamics of two-level open quantum systems in a dephasing channel and an amplitude damping channel computed on a quantum device. The algorithm presented will be most useful for implementing general non-unitary operations when the SVD can be readily computed, which is the case with most operators in the noisy intermediate-scale quantum computing era.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Reducing the Quantum Many-electron Problem to Two Electrons with Machine Learning

An outstanding challenge in chemical computation is the many-electron problem where computational methodologies scale prohibitively with system size. The energy of any molecule can be expressed as a weighted sum of the energies of two-electron wave functions that are computable from only a two-electron calculation. Despite the physical elegance of this extended ``aufbau'' principle, the determination of the distribution of weights -- geminal occupations -- for general molecular systems has remained elusive. Here we introduce a new paradigm for electronic structure where approximate geminal-occupation distributions are ``learned'' via a convolutional neural network. We show that the neural network learns the $N$-representability conditions, constraints on the distribution for it to represent an $N$-electron system. By training on hydrocarbon isomers with only 2-7 carbon atoms, we are able to predict the energies for isomers of octane as well as hydrocarbons with 8-15 carbons. The present work demonstrates that machine learning can be used to reduce the many-electron problem to an effective two-electron problem, opening new opportunities for accurately predicting electronic structure.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Simultaneous fermion and exciton condensations from a model Hamiltonian

Fermion-exciton condensation in which both fermion-pair (i.e., superconductivity) and exciton condensations occur simultaneously in a single coherent quantum state has recently been conjectured to exist. Here, we capture the fermion-exciton condensation through a model Hamiltonian that can recreate the physics of this new class of highly correlated condensation phenomena. We demonstrate that the Hamiltonian generates the large-eigenvalue signatures of fermion-pair and exciton condensations for a series of states with increasing particle numbers. The results confirm that the dual-condensate wave function arises from the entanglement of fermion-pair and exciton wave functions, which we previously predicted in the thermodynamic limit. This model Hamiltonian -- generalizing well-known model Hamiltonians for either superconductivity or exciton condensation -- can explore a wide variety of condensation behavior. It provides significant insights into the required forces for generating a fermion-exciton condensate, which will likely be invaluable for realizing such condensations in realistic materials with applications from superconductors to excitonic materials.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Correlation-Driven Phenomena in Periodic Molecular Systems from Variational Two-electron Reduced Density Matrix Theory

Correlation-driven phenomena in molecular periodic systems are challenging to predict computationally not only because such systems are periodically infinite but also because they are typically strongly correlated. Here we generalize the variational two-electron reduced density matrix (2-RDM) theory to compute the energies and properties of strongly correlated periodic systems. The 2-RDM of the unit cell is directly computed subject to necessary $N$-representability conditions such that the unit-cell 2-RDM represents at least one $N$-electron density matrix. Two canonical but non-trivial systems, periodic metallic hydrogen chains and periodic acenes, are treated to demonstrate the methodology. We show that, while single-reference correlation theories do not capture the strong (static) correlation effects in either of these molecular systems, the periodic variational 2-RDM theory predicts the Mott metal-to-insulator transition in the hydrogen chains and the length-dependent polyradical formation in acenes. For both hydrogen chains and acenes the periodic calculations are compared with previous non-periodic calculations with the results showing a significant change in energies and increase in the electron correlation from the periodic boundary conditions. The 2-RDM theory, which allows for much larger active spaces than are traditionally possible, is applicable to studying correlation-driven phenomena in general periodic molecular solids and materials.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Quantum Simulation of Molecules without Fermionic Encoding of the Wave Function

Molecular simulations generally require fermionic encoding in which fermion statistics are encoded into the qubit representation of the wave function. Recent calculations suggest that fermionic encoding of the wave function can be bypassed, leading to more efficient quantum computations. Here we show that the energy can be expressed as a functional of the two-electron reduced density matrix (2-RDM) where the 2-RDM is a unique functional of the unencoded $N$-qubit-particle wave function. Contrasts are made with current hardware-efficient methods. An application to computing the ground-state energy and 2-RDM of H$_{4}$ is presented.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Quantum Simulation of Open Quantum Systems Using a Unitary Decomposition of Operators

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Quantum-Classical Hybrid Algorithm for the Simulation of All-Electron Correlation

While the treatment of chemically relevant systems containing hundreds or even thousands of electrons remains beyond the reach of quantum devices, the development of quantum-classical hybrid algorithms to resolve electronic correlation presents a promising pathway toward a quantum advantage in the computation of molecular electronic structure. Such hybrid algorithms treat the exponentially scaling part of the calculation -- the static (multireference) correlation -- on the quantum computer and the non-exponentially scaling part -- the dynamic correlation -- on the classical computer. While a variety of such algorithms have been proposed, due to the dependence on the wave function of most classical methods for dynamic correlation, the development of easy-to-use classical post-processing implementations has been limited. Here we present a novel hybrid-classical algorithm that computes a molecule's all-electron energy and properties on the classical computer from a critically important simulation of the static correlation on the quantum computer. Significantly, for the all-electron calculations we circumvent the wave function by using density-matrix methods that only require input of the statically correlated two-electron reduced density matrix (2-RDM), which can be efficiently measured in the quantum simulation. Although the algorithm is completely general, we test it with two classical 2-RDM methods, the anti-Hermitian contracted Schrödinger equation (ACSE) theory and multiconfiguration pair-density functional theory (MC-PDFT), using the recently developed quantum ACSE method for the simulation of the statically correlated 2-RDM. We obtain experimental accuracy for the relative energies of all three benzyne isomers and thereby, demonstrate the ability of the quantum-classical hybrid algorithms to achieve chemically relevant results and accuracy on currently available quantum computers.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Relaxation of a Stationary State on a Quantum Computer Yields Unique Spectroscopic Fingerprint of the Computer's Noise

Quantum computing has the potential to revolutionize computing for certain classes of problems with exponential scaling, and yet this potential is accompanied by significant sensitivity to noise, requiring sophisticated error correction and mitigation strategies. Here we simulate the relaxations of stationary states at different frequencies on several quantum computers to obtain unique spectroscopic fingerprints of their noise. Response functions generated from the data reveal a clear signature of non-Markovian dynamics, demonstrating that each of the quantum computers acts as a non-Markovian bath with a unique colored noise profile. The study suggest that noisy intermediate-scale quantum computers (NISQ) provide a built-in noisy bath that can be analyzed from their simulation of closed quantum systems with the results potentially being harnessed for error mitigation or open-system simulation.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Resolving Correlated States of Benzyne on a Quantum Computer with an Error-Mitigated Quantum Contracted Eigenvalue Solver

The simulation of strongly correlated many-electron systems is one of the most promising applications for near-term quantum devices. Here we use a class of eigenvalue solvers (presented in Phys. Rev. Lett. 126, 070504 (2021)) in which a contraction of the Schrödinger equation is solved for the two-electron reduced density matrix (2-RDM) to resolve the energy splittings of ortho-, meta-, and para-isomers of benzyne ${\textrm C_6} {\textrm H_4}$. In contrast to the traditional variational quantum eigensolver, the contracted quantum eigensolver solves an integration (or contraction) of the many-electron Schrödinger equation onto the two-electron space. The quantum solution of the anti-Hermitian part of the contracted Schrödinger equation (qACSE) provides a scalable approach with variational parameters that has its foundations in 2-RDM theory. Experimentally, a variety of error mitigation strategies enable the calculation, including a linear shift in the 2-RDM targeting the iterative nature of the algorithm as well as a projection of the 2-RDM onto the convex set of approximately $N$-representable 2-RDMs defined by the 2-positive (DQG) $N$-representability conditions. The relative energies exhibit single-digit millihartree errors, capturing a large part of the electron correlation energy, and the computed natural orbital occupations reflect the significant differences in the electron correlation of the isomers.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Efficient Two-Electron Ansatz for Benchmarking Quantum Chemistry on a Quantum Computer

Quantum chemistry provides key applications for near-term quantum computing, but these are greatly complicated by the presence of noise. In this work we present an efficient ansatz for the computation of two-electron atoms and molecules within a hybrid quantum-classical algorithm. The ansatz exploits the fundamental structure of the two-electron system, and treating the nonlocal and local degrees of freedom on the quantum and classical computers, respectively. Here the nonlocal degrees of freedom scale linearly with respect to basis-set size, giving a linear ansatz with only $\mathcal{O}(1)$ circuit preparations required for reduced state tomography. We implement this benchmark with error mitigation on two publicly available quantum computers, calculating accurate dissociation curves for 4- and 6- qubit calculations of ${\rm H}_\textrm{2}^{}$ and ${\rm H}_\textrm{3}^+$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Experimental Data from a Quantum Computer Verifies the Generalized Pauli Exclusion Principle

"What are the consequences ... that Fermi particles cannot get into the same state ... " R. P. Feynman wrote of the Pauli exclusion principle, "In fact, almost all the peculiarities of the material world hinge on this wonderful fact." In 1972 Borland and Dennis showed that there exist powerful constraints beyond the Pauli exclusion principle on the orbital occupations of Fermi particles, providing important restrictions on quantum correlation and entanglement. Here we use computations on quantum computers to experimentally verify the existence of these additional constraints. Quantum many-fermion states are randomly prepared on the quantum computer and tested for constraint violations. Measurements show no violation and confirm the generalized Pauli exclusion principle with an error of one part in one quintillion.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Potential Coexistence of Exciton and Fermion Pair Condensations

Extensive theoretical and experimental investigation has been conducted on fermion pair condensation and exciton condensation as distinct classes of Bose-Einstein-like condensation. In this work, the existence of a fermion-exciton condensate---a single quantum state in which character of both fermion pair and exciton condensates coexist---is established computationally in the low-particle-number ($N$) limit and theoretically in the large-$N$ thermodynamic limit. The trade-off between the fermion pair and excitonic character of the fermion-exciton condensate is shown to be elliptic in nature. The possibility that the properties of fermion-exciton condensates could be a hybrid of the properties of fermion pair condensates and exciton condensates is discussed, and future experimental and computational exploration of this new class of condensate, which may potentially be realizable in a bilayer of superconductors, is anticipated.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Quantum-classical hybrid algorithm using an error-mitigating $N$-representability condition to compute the Mott metal-insulator transition

Quantum algorithms for molecular electronic structure have been developed with lower computational scaling than their classical counterparts, but emerging quantum hardware is far from being capable of the coherence,connectivity and gate errors required for their experimental realization. Here we propose a class of quantum-classical hybrid algorithms that compute the energy from a two-electron reduced density matrix (2-RDM). The 2-RDM is constrained by $N$-representability conditions, conditions for representing an $N$-electron wavefunction, that mitigates noise from the quantum circuit. We compute the strongly correlated dissociation of doublet H$_{3}$ into three hydrogen atoms. The hybrid quantum-classical computer matches the energies from full configuration interaction to 0.1 kcal/mol, one-tenth of "chemical accuracy," even in the strongly correlated limit of dissociation. Furthermore, the spatial locality of the computed one-electron RDM reveals that the quantum computer accurately predicts the Mott metal-insulator transition.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Accurate Non-adiabatic Quantum Dynamics from Pseudospectral Sampling of Time-dependent Gaussian Basis Sets

Quantum molecular dynamics requires an accurate representation of the molecular potential energy surface from a minimal number of electronic structure calculations, particularly for nonadiabatic dynamics where excited states are required. In this paper, we employ pseudospectral sampling of time-dependent Gaussian basis functions for the simulation of non-adiabatic dynamics. Unlike other methods, the pseudospectral Gaussian molecular dynamics tests the Schrödinger equation with $N$ Dirac delta functions located at the centers of the Gaussian functions reducing the scaling of potential energy evaluations from $\mathcal{O}(N^2)$ to $\mathcal{O}(N)$. By projecting the Gaussian basis onto discrete points in space, the method is capable of efficiently and quantitatively describing nonadiabatic population transfer and intra-surface quantum coherence. We investigate three model systems; the photodissociation of three coupled Morse oscillators, the bound state dynamics of two coupled Morse oscillators, and a two-dimensional model for collinear triatomic vibrational dynamics. In all cases, the pseudospectral Gaussian method is in quantitative agreement with numerically exact calculations. The results are promising for nonadiabatic molecular dynamics in molecular systems where strongly correlated ground or excited states require expensive electronic structure calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Necessary N-representability Constraints from Time-reversal Symmetry for Periodic Systems

The variational calculation of the two-electron reduced density matrix (2-RDM) is extended to periodic molecular systems. If the 2-RDM theory is extended to the periodic case without consideration of time-reversal symmetry, however, it can yields energies that are significantly lower than the correct energies. We derive and implement linear constraints that enforce time-reversal symmetry on the 2-RDM without destroying its computationally favorable block-diagonal structure from translational invariance. Time-reversal symmetry is distinct from space-group or spin (SU(2)) symmetries which can be expressed by unitary transformations. The time-reversal symmetry constraints are demonstrated through calculations of the metallic hydrogen chain and the one-dimensional lithium hydride crystal.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Comparison of one-dimensional and quasi-one-dimensional Hubbard models from the variational two-electron reduced-density-matrix method

Minimizing the energy of an $N$-electron system as a functional of a two-electron reduced density matrix (2-RDM), constrained by necessary $N$-representability conditions (conditions for the 2-RDM to represent an ensemble $N$-electron quantum system), yields a rigorous lower bound to the ground-state energy in contrast to variational wavefunction methods. We characterize the performance of two sets of approximate constraints, (2,2)-positivity (DQG) and approximate (2,3)-positivity (DQGT) conditions, at capturing correlation in one-dimensional and quasi-one-dimensional (ladder) Hubbard models. We find that, while both the DQG and DQGT conditions capture both the weak and strong correlation limits, the more stringent DQGT conditions improve the ground-state energies, the natural occupation numbers, the pair correlation function, the effective hopping, and the connected (cumulant) part of the 2-RDM. We observe that the DQGT conditions are effective at capturing strong electron correlation effects in both one- and quasi-one-dimensional lattices for both half filling and less-than-half filling.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Generalized Pauli conditions on the spectra of one-electron reduced density matrices of atoms and molecules

The Pauli exclusion principle requires the spectrum of the occupation numbers of the one-electron reduced density matrix (1-RDM) to be bounded by one and zero. However, for a 1-RDM from a wave function, there exist additional conditions on the spectrum of occupation numbers, known as pure N-representability conditions or generalized Pauli conditions. For atoms and molecules, we measure through a Euclidean-distance metric the proximity of the 1-RDM spectrum to the facets of the convex set (polytope) generated by the generalized Pauli conditions. For the ground state of any spin symmetry, as long as time-reversal symmetry is considered in the definition of the polytope, we find that the 1-RDM's spectrum is pinned to the boundary of the polytope. In contrast, for excited states, we find that the 1-RDM spectrum is not pinned. Proximity of the 1-RDM to the boundary of the polytope provides a measurement and classification of electron correlation and entanglement within the quantum system. For comparison, this distance to the boundary of the generalized Pauli conditions is also compared to the distance to the polytope of the traditional Pauli conditions, and the distance to the nearest 1-RDM spectrum from a Slater determinant. We explain the difference in pinning in the ground- and excited-state 1-RDMs through a connection to the N-representability conditions of the two-electron reduced density matrix.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Global Solutions of Hartree-Fock Theory and their Consequences for Strongly Correlated Quantum Systems

We present a density matrix approach for computing global solutions of Hartree-Fock theory, based on semidefinite programming (SDP), that gives upper and lower bounds on the Hartree-Fock energy of quantum systems. Equality of the upper- and lower-bound energies guarantees that the computed solution is the globally optimal solution of Hartree-Fock theory. For strongly correlated systems the SDP approach requires us to reassess the accuracy of the Hartree-Fock energies and densities from standard software packages for electronic structure theory. Calculations of the H$_{4}$ dimer and N$_{2}$ molecule show that the energies from SDP Hartree-Fock are lower than those from standard Hartree-Fock methods by 100-200 kcal/mol in the dissociation region. The present findings have important consequences for the computation and interpretation of electron correlation, which is typically defined relative to the Hartree-Fock energy and density.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Global solutions of restricted open-shell Hartree-Fock theory from semidefinite programming with applications to strongly correlated quantum systems

We present a density matrix approach for computing global solutions of restricted open-shell Hartree-Fock theory, based on semidefinite programming (SDP), that gives upper and lower bounds on the Hartree-Fock energy of quantum systems. While wave function approaches to Hartree-Fock theory yield an upper bound to the Hartree-Fock energy, we derive a semidefinite relaxation of Hartree-Fock theory that yields a rigorous lower bound on the Hartree-Fock energy. We also develop an upper-bound algorithm in which Hartree-Fock theory is cast as a SDP with a nonconvex constraint on the rank of the matrix variable. Equality of the upper- and lower-bound energies guarantees that the computed solution is the globally optimal solution of Hartree-Fock theory. The work extends a previously presented method for closed-shell systems [S. Veeraraghavan and D. A. Mazziotti, Phys. Rev. A89, 010502R (2014)]. For strongly correlated systems the SDP approach provides an alternative to the locally optimized Hartree-Fock energies and densities with a certificate of global optimality. Applications are made to the potential energy curves of C2, CN, Cr2, and NO2.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Molecule-optimized Basis Sets and Hamiltonians for Accelerated Electronic Structure Calculations of Atoms and Molecules

Molecule-optimized basis sets, based on approximate natural orbitals, are developed for accelerating the convergence of quantum calculations with strongly correlated (multi-referenced) electrons. We use a low-cost approximate solution of the anti-Hermitian contracted Schr{ö}dinger equation (ACSE) for the one- and two-electron reduced density matrices (RDMs) to generate an approximate set of natural orbitals for strongly correlated quantum systems. The natural-orbital basis set is truncated to generate a molecule-optimized basis set whose rank matches that of a standard correlation-consistent basis set optimized for the atoms. We show that basis-set truncation by approximate natural orbitals can be viewed as a one-electron unitary transformation of the Hamiltonian operator and suggest an extension of approximate natural-orbital truncations through two-electron unitary transformations of the Hamiltonian operator, such as those employed in the solution of the ACSE. The molecule-optimized basis set from the ACSE improves the accuracy of the equivalent standard atom-optimized basis set at little additional computational cost. We illustrate the method with the potential energy curves of hydrogen fluoride and diatomic nitrogen. Relative to the hydrogen fluoride potential energy curve from the ACSE in a polarized triple-zeta basis set, the ACSE curve in a molecule-optimized basis set, equivalent in size to a polarized double-zeta basis, has a nonparallelity error of 0.0154 a.u. which is significantly better than the nonparallelity error of 0.0252 a.u. from the polarized double-zeta basis set.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Effect of Strong Electron Correlation on the Efficiency of Photosynthetic Light Harvesting

Research into the efficiency of photosynthetic light harvesting has focused on two factors: (1) entanglement of chromophores, and (2) environmental noise. While chromophores are conjugated $π$-bonding molecules with strongly correlated electrons, previous models have treated this correlation implicitly without a mathematical variable to gauge correlation-enhanced efficiency. Here we generalize the single-electron/exciton models to a multi-electron/exciton model that explicitly shows the effects of enhanced electron correlation within chromophores on the efficiency of energy transfer. The model provides more detailed insight into the interplay of electron correlation within chromophores and electron entanglement between chromophores. Exploiting this interplay is assisting in the design of new energy-efficient materials, which are just beginning to emerge.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Significant Conditions on the Two-electron Reduced Density Matrix from the Constructive Solution of N-representability

We recently presented a constructive solution to the N-representability problem of the two-electron reduced density matrix (2-RDM)---a systematic approach to constructing complete conditions to ensure that the 2-RDM represents a realistic N-electron quantum system [D. A. Mazziotti, Phys. Rev. Lett. 108, 263002 (2012)]. In this paper we provide additional details and derive further N-representability conditions on the 2-RDM that follow from the constructive solution. The resulting conditions can be classified into a hierarchy of constraints, known as the (2,q)-positivity conditions where the q indicates their derivation from the nonnegativity of q-body operators. In addition to the known T1 and T2 conditions, we derive a new class of (2,3)-positivity conditions. We also derive 3 classes of (2,4)-positivity conditions, 6 classes of (2,5)-positivity conditions, and 24 classes of (2,6)-positivity conditions. The constraints obtained can be divided into two general types: (i) lifting conditions, that is conditions which arise from lifting lower (2,q)-positivity conditions to higher (2,q+1)-positivity conditions and (ii) pure conditions, that is conditions which cannot be derived from a simple lifting of the lower conditions. All of the lifting conditions and the pure (2,q)-positivity conditions for q>3 require tensor decompositions of the coefficients in the model Hamiltonians. Subsets of the new N-representability conditions can be employed with the previously known conditions to achieve polynomially scaling calculations of ground-state energies and 2-RDMs of many-electron quantum systems even in the presence of strong electron correlation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Structure of Fermionic Density Matrices: Complete N-representability Conditions

We present a constructive solution to the N-representability problem---a full characterization of the conditions for constraining the two-electron reduced density matrix (2-RDM) to represent an N-electron density matrix. Previously known conditions, while rigorous, were incomplete. Here we derive a hierarchy of constraints built upon (i) the bipolar theorem and (ii) tensor decompositions of model Hamiltonians. Existing conditions D, Q, G, T1, and T2, known classical conditions, and new conditions appear naturally. Subsets of the conditions are amenable to polynomial-time computations of strongly correlated systems.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Impact of Multichannel and Multipole Effects on the Cooper Minimum in the High-Harmonics Spectrum of Argon

We investigate the relevance of multiple-orbital and multipole effects during high-harmonic generation (HHG). The time-dependent configuration-interaction singles (TDCIS) approach is used to study the impact of the detailed description of the residual electron-ion interaction on the HHG spectrum. We find that the shape and position of the Cooper minimum in the HHG spectrum of argon changes significantly whether or not interchannel interactions are taken into account. The HHG yield can be underestimated by up to 2 orders of magnitude in the energy regio of 30-50 eV. We show that the argument of low ionization probability is not sufficient to justify ignoring multiple-orbital contributions. Additionally, we find the HHG yield is sensitive to the nonspherical multipole character of the electron-ion interaction.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Decoherence in Attosecond Photoionization

The creation of superpositions of hole states via single-photon ionization using attosecond extreme-ultraviolet pulses is studied with the time-dependent configuration interaction singles (TDCIS) method. Specifically, the degree of coherence between hole states in atomic xenon is investigated. We find that interchannel coupling not only affects the hole populations, it also enhances the entanglement between the photoelectron and the remaining ion, thereby reducing the coherence within the ion. As a consequence, even if the spectral bandwidth of the ionizing pulse exceeds the energy splittings among the hole states involved, perfectly coherent hole wave packets cannot be formed. For sufficiently large spectral bandwidth, the coherence can only be increased by increasing the mean photon energy.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Functional Subsystems and Quantum Redundancy in Photosynthetic Light Harvesting

The Fenna-Matthews-Olson (FMO) antennae complex, responsible for light harvesting in green sulfur bacteria, consists of three monomers, each with seven chromophores. Here we show that multiple subsystems of the seven chromophores can transfer energy from either chromophore 1 or 6 to the reaction center with an efficiency matching or in many cases exceeding that of the full seven chromophore system. In the FMO complex these functional subsystems support multiple quantum pathways for efficient energy transfer that provide a built-in quantum redundancy. There are many instances of redundancy in nature, providing reliability and protection, and in photosynthetic light harvesting this quantum redundancy provides protection against the temporary or permanent loss of one or more chromophores. The complete characterization of functional subsystems within the FMO complex offers a detailed map of the energy flow within the FMO complex, which has potential applications to the design of more efficient photovoltaic devices.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Rank Restriction for the Variational Calculation of Two-electron Reduced Density Matrices of Many-electron Atoms and Molecules

Variational minimization of the ground-state energy as a function of the two-electron reduced density matrix (2-RDM), constrained by necessary $N$-representability conditions, provides a polynomial-scaling approach to studying strongly correlated molecules without computing the many-electron wavefunction. Here we introduce a new route to enhancing necessary conditions for $N$-representability through rank restriction of the 2-RDM. Rather than adding computationally more expensive $N$-representability conditions, we directly enhance the accuracy of two-particle (2-positivity) conditions through rank restriction, which removes degrees of freedom in the 2-RDM that are not sufficiently constrained. We select the rank of the particle-hole 2-RDM by deriving the ranks associated with model wavefunctions including both mean-field and antisymmetrized geminal power (AGP) wave functions. Because the 2-positivity conditions are exact for quantum systems with AGP ground states, the rank of the particle-hole 2-RDM from the AGP ansatz provides a minimum for its value in variational 2-RDM calculations of general quantum systems. To implement the rank-restricted conditions, we extend a first-order algorithm for large-scale semidefinite programming. The rank-restricted conditions significantly improve the accuracy of the energies; for example, the percentages of correlation energies recovered for HF, CO, and N$_{2}$ improve from 115.2%, 121.7%, and 121.5% without rank restriction to 97.8%, 101.1%, and 100.0% with rank restriction. Similar results are found at both equilibrium and non-equilibrium geometries. While more accurate, the rank-restricted $N$-representability conditions are less expensive computationally than the full-rank conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

On the relationship between parametric two-electron reduced-density-matrix methods and the coupled electron pair approximation

Parametric two-electron reduced-density-matrix (p-2RDM) methods have enjoyed much success in recent years; the methods have been shown to exhibit accuracies greater than coupled cluster with single and double substitutions (CCSD) for both closed- and open-shell ground-state energies, properties, geometric parameters, and harmonic frequencies. The class of methods is herein discussed within the context of the coupled electron pair approximation (CEPA), and several CEPA-like topological factors are presented for use within the p-2RDM framework. The resulting p-2RDM/n methods can be viewed as a density-based generalization of CEPA/n family that are numerically very similar to traditional CEPA methodologies. We cite the important distinction that the obtained energies represent stationary points, facilitating the efficient evaluation of properties and geometric derivatives. The p-2RDM/n formalism is generalized for an equal treatment of exclusion-principle-violating (EPV) diagrams that occur in the occupied and virtual spaces. One of these general topological factors is shown to be identical to that proposed by Kollmar [C. Kollmar, J. Chem. Phys. 125, 084108 (2006)], derived in an effort to approximately enforce the D, Q, and G conditions for N-representability in his size-extensive density matrix functional.

0 citations0 reviewsNo code linkedOpen access