Determining the N-representability of a reduced density matrix via unitary evolution and stochastic sampling

Abstract

The N-representability problem consists in determining whether, for a given p-body matrix, there exists at least one N-body density matrix from which the p-body matrix can be obtained by contraction, that is, if the given matrix is a p-body reduced density matrix (p-RDM). The knowledge of all necessary and sufficient conditions for a p-body matrix to be N-representable allows the constrained minimization of a many-body Hamiltonian expectation value with respect to the p-body density matrix and, thus, the determination of its exact ground state. However, the number of constraints that complete the N-representability conditions grows exponentially with system size, and hence the procedure quickly becomes intractable in practice. This work introduces a hybrid quantum-stochastic algorithm to effectively replace the N-representability conditions. The algorithm consists of applying to an initial N-body density matrix a sequence of unitary evolution operators constructed from a stochastic process that successively approaches the reduced state of the density matrix on a p-body subsystem, represented by a p-RDM, to a target p-body matrix, potentially a p-RDM. The generators of the evolution operators follow the adaptive derivative-assembled pseudo-Trotter method (ADAPT), while the stochastic component is implemented using a simulated annealing process. The resulting algorithm is independent of any underlying Hamiltonian, and it can be used to decide if a given p-body matrix is N-representable, establishing a criterion to determine its quality and correcting it. We apply this hybrid ADAPT algorithm to alleged reduced density matrices from a quantum chemistry electronic Hamiltonian, the reduced BCS model with constant pairing, and the Heisenberg XXZ spin model. In all cases, the proposed method behaves as expected for 1-RDMs and 2-RDMs, evolving the initial matrices towards different targets.

References in corpus (33)

• Barren plateaus in quantum neural network training landscapes
• Noisy intermediate-scale quantum (NISQ) algorithms
• Recent developments in the PySCF program package
• The Variational Quantum Eigensolver: a review of methods and best practices
• An adaptive variational algorithm for exact molecular simulations on a quantum computer
• Scalable Quantum Simulation of Molecular Energies
• Pairing in nuclear systems: from neutron stars to finite nuclei
• Generalized Unitary Coupled Cluster Wavefunctions for Quantum Computation
• qubit-ADAPT-VQE: An adaptive algorithm for constructing hardware-efficient ansatze on a quantum processor
• Quantum algorithms for electronic structure calculations: particle/hole Hamiltonian and optimized wavefunction expansions
• Structure of Fermionic Density Matrices: Complete N-representability Conditions
• Exact Parameterization of Fermionic Wave Functions via Unitary Coupled Cluster Theory
• N-representability is QMA-complete
• Parity-Affected Superconductivity in Ultrasmall Metallic Grains
• The Pauli principle revisited
• Qubit-excitation-based adaptive variational quantum eigensolver
• Quantum Solver of Contracted Eigenvalue Equations for Scalable Molecular Simulations on Quantum Computing Devices
• On the difference between variational and unitary coupled cluster theories
• Pinning of Fermionic Occupation Numbers
• Exact Study of the Effect of Level Statistics in Ultrasmall Superconducting Grains
• Prospects of Quantum Computing for Molecular Sciences
• Quantum-classical hybrid algorithm using an error-mitigating $N$-representability condition to compute the Mott metal-insulator transition
• Relating the pure and ensemble density matrix functional
• Generalized Pauli conditions on the spectra of one-electron reduced density matrices of atoms and molecules
• Generalized Pauli constraints in reduced density matrix functional theory
• Resolving Correlated States of Benzyne on a Quantum Computer with an Error-Mitigated Quantum Contracted Eigenvalue Solver
• Solving condensed-matter ground-state problems by semidefinite relaxations
• Quantum Many-body Theory from a Solution of the $N$-representability Problem
• Experimental Data from a Quantum Computer Verifies the Generalized Pauli Exclusion Principle
• Exact Ansatz of Fermion-Boson Systems for a Quantum Device
• Which options exist for NISQ-friendly linear response formulations?
• Benchmarking the variational reduced density matrix theory in the doubly-occupied configuration interaction space with integrable pairing models
• Physically motivated improvements of Variational Quantum Eigensolvers

Cited by in corpus (2)

• Determining the ensemble N-representability of Reduced Density Matrices
• Is the matrix completion of reduced density matrices unique?