Determining the N-representability of a reduced density matrix via unitary evolution and stochastic sampling
- URL: http://arxiv.org/abs/2503.17303v1
- Date: Fri, 21 Mar 2025 16:52:22 GMT
- Title: Determining the N-representability of a reduced density matrix via unitary evolution and stochastic sampling
- Authors: Gustavo. E. Massaccesi, Ofelia. B. Oña, Pablo Capuzzi, Juan I. Melo, Luis Lain, Alicia Torre, Juan E. Peralta, Diego R. Alcoba, Gustavo E. Scuseria,
- Abstract summary: This work introduces a hybrid quantum-stochastic algorithm to effectively replace the N-representability conditions.<n>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.
- Score: 0.0
- License: http://arxiv.org/licenses/nonexclusive-distrib/1.0/
- 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.
Related papers
- Determining the ensemble N-representability of Reduced Density Matrices [0.0]
We propose a framework for determining the ensemble N-representability of a p-body matrix.<n>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.
arXiv Detail & Related papers (2026-02-05T20:11:03Z) - Quantum Hermitian conjugate and encoding unnormalized matrices [49.494595696663524]
We develop the family of matrix-manipulation algorithms based on the encoding the matrix elements into the probability amplitudes of the pure superposition state of a certain quantum system.
We introduce two extensions to these algorithms which allow (i) to perform Hermitian conjugation of matrices under consideration and (ii) to weaken the restriction to the absolute values of matrix elements unavoidably imposed by the normalization condition for a pure quantum state.
arXiv Detail & Related papers (2025-03-27T08:49:59Z) - Matrix encoding method in variational quantum singular value decomposition [49.494595696663524]
Conditional measurement is involved to avoid small success probability in ancilla measurement.
The objective function for the algorithm can be obtained probabilistically via measurement of the state of a one-qubit subsystem.
arXiv Detail & Related papers (2025-03-19T07:01:38Z) - Projective purification of correlated reduced density matrices [0.0]
We present an algorithm capable of performing all of the following tasks in the least invasive manner.<n>We demonstrate the superiority of the present purification algorithm over previous ones in the context of the Fermi-Hubbard model.
arXiv Detail & Related papers (2024-12-18T07:33:51Z) - Simulating NMR Spectra with a Quantum Computer [49.1574468325115]
This paper provides a formalization of the complete procedure of the simulation of a spin system's NMR spectrum.
We also explain how to diagonalize the Hamiltonian matrix with a quantum computer, thus enhancing the overall process's performance.
arXiv Detail & Related papers (2024-10-28T08:43:40Z) - Semi-supervised Symmetric Non-negative Matrix Factorization with Low-Rank Tensor Representation [27.14442336413482]
Semi-supervised symmetric non-negative matrix factorization (SNMF)
We propose a novel SNMF model by seeking low-rank representation for the tensor synthesized by the pairwise constraint matrix.
We then propose an enhanced SNMF model, making the embedding matrix tailored to the above tensor low-rank representation.
arXiv Detail & Related papers (2024-05-04T14:58:47Z) - Vectorization of the density matrix and quantum simulation of the von
Neumann equation of time-dependent Hamiltonians [65.268245109828]
We develop a general framework to linearize the von-Neumann equation rendering it in a suitable form for quantum simulations.
We show that one of these linearizations of the von-Neumann equation corresponds to the standard case in which the state vector becomes the column stacked elements of the density matrix.
A quantum algorithm to simulate the dynamics of the density matrix is proposed.
arXiv Detail & Related papers (2023-06-14T23:08:51Z) - Quantum Many-body Theory from a Solution of the $N$-representability
Problem [0.0]
We derive an equation that re-expresses physical constraints on higher-order RDMs to generate direct constraints on the 2-RDM.
We illustrate by computing the ground-state electronic energy and properties of the H$_8$ ring.
arXiv Detail & Related papers (2023-04-17T19:19:31Z) - A hybrid quantum-classical algorithm for multichannel quantum scattering
of atoms and molecules [62.997667081978825]
We propose a hybrid quantum-classical algorithm for solving the Schr"odinger equation for atomic and molecular collisions.
The algorithm is based on the $S$-matrix version of the Kohn variational principle, which computes the fundamental scattering $S$-matrix.
We show how the algorithm could be scaled up to simulate collisions of large polyatomic molecules.
arXiv Detail & Related papers (2023-04-12T18:10:47Z) - A Novel Maximum-Entropy-Driven Technique for Low-Rank Orthogonal
Nonnegative Matrix Factorization with $\ell_0$-Norm sparsity Constraint [0.0]
In data-driven control and machine learning, a common requirement involves breaking down large matrices into smaller, low-rank factors.
This paper introduces an innovative solution to the orthogonal nonnegative matrix factorization (ONMF) problem.
The proposed method achieves comparable or improved reconstruction errors in line with the literature.
arXiv Detail & Related papers (2022-10-06T04:30:59Z) - Direct solution of multiple excitations in a matrix product state with
block Lanczos [62.997667081978825]
We introduce the multi-targeted density matrix renormalization group method that acts on a bundled matrix product state, holding many excitations.
A large number of excitations can be obtained at a small bond dimension with highly reliable local observables throughout the chain.
arXiv Detail & Related papers (2021-09-16T18:36:36Z) - Understanding Implicit Regularization in Over-Parameterized Single Index
Model [55.41685740015095]
We design regularization-free algorithms for the high-dimensional single index model.
We provide theoretical guarantees for the induced implicit regularization phenomenon.
arXiv Detail & Related papers (2020-07-16T13:27:47Z) - Variational optimization of continuous matrix product states [0.0]
We show how to optimize continuous matrix product states for systems with inhomogeneous external potentials.
We show how both the energy and its backwards derivative can be calculated exactly and at a cost that scales as the cube of the bond dimension.
We illustrate this by finding ground states of interacting bosons in external potentials, and by calculating boundary or Casimir energy corrections of continuous many-body systems.
arXiv Detail & Related papers (2020-06-02T17:35:41Z) - One-body entanglement as a quantum resource in fermionic systems [0.0]
We show that one-body entanglement can be considered as a quantum resource.
It is consistent with a model of fermionic quantum computation which requires correlations beyond antisymmetrization.
arXiv Detail & Related papers (2020-01-10T17:27:09Z)
This list is automatically generated from the titles and abstracts of the papers in this site.
This site does not guarantee the quality of this site (including all information) and is not responsible for any consequences.