Determining the ensemble N-representability of Reduced Density Matrices

Abstract

The N-representability problem for reduced density matrices remains a fundamental challenge in electronic structure theory. Following our previous work that employs a unitary-evolution algorithm based on an adaptive derivative-assembled pseudo-Trotter variational quantum algorithm to probe pure-state N-representability of reduced density matrices [J. Chem. Theory Comput. 2024, 20, 9968], in this work we propose a practical framework for determining the ensemble N-representability of a p-body matrix. This is accomplished using a purification strategy consisting of embedding an ensemble state into a pure state defined on an extended Hilbert space, such that the reduced density matrices of the purified state reproduce those of the original ensemble. By iteratively applying variational unitaries to an initial purified state, the proposed algorithm minimizes the Hilbert-Schmidt distance between its p-body reduced density matrix and a specified target p-body matrix, which serves as a measure of the N-representability of the target. This methodology facilitates both error correction of defective ensemble reduced density matrices, and quantum-state reconstruction on a quantum computer, offering a route for density-matrix refinement. We validate the algorithm with numerical simulations on systems of two, three, and four electrons in both, simple models as well as molecular systems at finite temperature, demonstrating its robustness.

References in corpus (32)

• 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
• 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
• Quantum marginal problem and N-representability
• Exact Parameterization of Fermionic Wave Functions via Unitary Coupled Cluster Theory
• N-representability is QMA-complete
• 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
• Pinning of Fermionic Occupation Numbers
• Prospects of Quantum Computing for Molecular Sciences
• Quantum state tomography via reduced density matrices
• 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
• Quasi-pinning and entanglement in the lithium isoelectronic series
• Thermofield Theory for Finite-Temperature Quantum Chemistry
• Quantum Many-body Theory from a Solution of the $N$-representability Problem
• Experimental Data from a Quantum Computer Verifies the Generalized Pauli Exclusion Principle
• Dual-Cone Variational Calculation of the 2-Electron Reduced Density Matrix
• Wave function methods for canonical ensemble thermal averages in correlated many-fermion systems
• Exact Ansatz of Fermion-Boson Systems for a Quantum Device
• Which options exist for NISQ-friendly linear response formulations?
• Physically motivated improvements of Variational Quantum Eigensolvers
• Thermofield theory for finite-temperature electronic structure
• EQC : Ensembled Quantum Computing for Variational Quantum Algorithms
• Determining the N-representability of a reduced density matrix via unitary evolution and stochastic sampling