Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2508.07296

Reduced density matrix and cumulant approximations of quantum linear response

Theo Juncker von Buchwald, Erik Rosendahl Kjellgren, Jacob Kongsted, Stephan P. A. Sauer, Sonia Coriani, Karl Michael Ziems

Published Aug 10, 2025• 0 citations
arXivHTMLPDF

Abstract

Linear response (LR) is an important tool in the computational chemist's toolbox. It is therefore no surprise that the emergence of quantum computers has led to a quantum version, quantum LR (qLR). However, the current quantum era of near-term intermediary scale quantum (NISQ) computers is dominated by noise, short decoherence times, and slow measurement speed. It is therefore of interest to find approximations that greatly reduce the quantum workload while only slightly impacting the quality of a method. In an effort to achieve this, we approximate the naive qLR with singles and doubles (qLRSD) method by either directly approximating the reduced density matrices (RDMs) or indirectly through their respective reduced density cumulants (RDCs). We present an analysis of the measurement costs behind qLR with RDMs, and report qLR results for model Hydrogen ladder systems; for varying active space sizes of OCS, SeH$_2$, and H$_2$S; and for symmetrically stretched H$_2$O and BeH$_2$. Discouragingly, while approximations to the 4-body RDMs and RDCs seem to produce good results for systems at the equilibrium geometry and for some types of core excitations, they both tend to fail when the system exhibits strong correlation. All approximations to the 3-body RDMs and/or RDCs severely affect the results and cannot be applied.

DOI
10.48550/arxiv.2508.07296

Related papers

Reduced Density Matrix Formulation of Quantum Linear Response

Theo Juncker von Buchwald, Karl Michael Ziems, Erik Rosendahl Kjellgren +3 more

The prediction of spectral properties via linear response (LR) theory is an important tool in quantum chemistry for understanding photoinduced processes in molecular systems. With the advances of quantum computing, we recently adapted this method for near-term quantum hardware using a truncated active space approximation with orbital rotation, named quantum linear response (qLR). In an effort to reduce the classic cost of this hybrid approach, we here derive and implement a reduced density matrix (RDM) driven approach of qLR. This allows for the calculation of spectral properties of moderately sized molecules with much larger basis sets than so far possible. We report qLR results for benzene and R -methyloxirane with a cc-pVTZ basis set and study the effect of shot noise on the valence and oxygen K-edge absorption spectra of H 2 O in the cc-pVTZ basis.

arXiv:2404.16586•Journal of Chemical Theory and Computation•Aug 6, 2024•9 citations

Which Options Exist for NISQ-Friendly Linear Response Formulations?

Karl Michael Ziems, Erik Rosendahl Kjellgren, Peter Reinholdt +4 more

Linear response (LR) theory is a powerful tool in classic quantum chemistry crucial to understanding photoinduced processes in chemistry and biology. However, performing simulations for large systems and in the case of strong electron correlation remains challenging. Quantum computers are poised to facilitate the simulation of such systems, and recently, a quantum linear response formulation (qLR) was introduced [Kumar et al., J. Chem. Theory Comput. 2023, 19, 9136–9150]. To apply qLR to near-term quantum computers beyond a minimal basis set, we here introduce a resource-efficient qLR theory, using a truncated active-space version of the multiconfigurational self-consistent field LR ansatz. Therein, we investigate eight different near-term qLR formalisms that utilize novel operator transformations that allow the qLR equations to be performed on near-term hardware. Simulating excited state potential energy curves and absorption spectra for various test cases, we identify two promising candidates, dubbed “proj LRSD” and “all-proj LRSD”.

arXiv:2312.13937•Journal of Chemical Theory and Computation•Apr 25, 2024•17 citations

Quantum Simulation of Molecular Response Properties in the NISQ Era

Ashutosh Kumar, Ayush Asthana, Vibin Abraham +6 more

Accurate modeling of the response of molecular systems to an external electromagnetic field is challenging on classical computers, especially in the regime of strong electronic correlation. In this article, we develop a quantum linear response (qLR) theory to calculate molecular response properties on near-term quantum computers. Inspired by the recently developed variants of the quantum counterpart of equation of motion (qEOM) theory, the qLR formalism employs "killer condition" satisfying excitation operator manifolds that offer a number of theoretical advantages along with reduced quantum resource requirements. We also used the qEOM framework in this work to calculate the state-specific response properties. Further, through noiseless quantum simulations, we show that response properties calculated using the qLR approach are more accurate than the ones obtained from the classical coupled-cluster-based linear response models due to the improved quality of the ground-state wave function obtained using the ADAPT-VQE algorithm.

Journal of Chemical Theory and Computation•Dec 6, 2023•32 citations

Quantum simulation of molecular response properties

Ashutosh Kumar, Ayush Asthana, Vibin Abraham +6 more

Accurate modeling of the response of molecular systems to an external electromagnetic field is challenging on classical computers, especially in the regime of strong electronic correlation. In this paper, we develop a quantum linear response (qLR) theory to calculate molecular response properties on near-term quantum computers. Inspired by the recently developed variants of the quantum counterpart of equation of motion (qEOM) theory, the qLR formalism employs "killer condition" satisfying excitation operator manifolds that offers a number of theoretical advantages along with reduced quantum resource requirements. We also used the qEOM framework in this work to calculate state-specific response properties. Further, through noise-less quantum simulations, we show that response properties calculated using the qLR approach are more accurate than the ones obtained from the classical coupled-cluster based linear response models due to the improved quality of the ground-state wavefunction obtained using the ADAPT-VQE algorithm.

arXiv:2301.06260•Jan 16, 2023•2 citations

Understanding and mitigating noise in molecular quantum linear response for spectroscopic properties on quantum computers

Karl Michael Ziems, Erik Rosendahl Kjellgren, Stephan P. A. Sauer +2 more

The promise of quantum computing to circumvent the exponential scaling of quantum chemistry has sparked a race to develop chemistry algorithms for quantum architecture. However, most works neglect the quantum-inherent shot noise, let alone the effect of current noisy devices. Here, we present a comprehensive study of quantum linear response (qLR) theory obtaining spectroscopic properties on simulated fault-tolerant quantum computers and present-day near-term quantum hardware. This work introduces novel metrics to analyze and predict the origins of noise in the quantum algorithm, proposes an Ansatz-based error mitigation technique, and reveals the significant impact of Pauli saving in reducing measurement costs and noise in subspace methods. Our hardware results using up to cc-pVTZ basis set serve as proof of principle for obtaining absorption spectra on quantum hardware in a general approach with the accuracy of classical multi-configurational methods. Importantly, our results exemplify that substantial improvements in hardware error rates and measurement speed are necessary to lift quantum computational chemistry from proof of concept to an actual impact in the field.

arXiv:2408.09308•Chemical Science•Jan 1, 2025•9 citations

Cumulant approach to the direct calculation of reduced density matrices: A critical analysis

Marcel Nooijen, Mark Wladyslawski, Anirban Hazra

In this paper we address a number of topics regarding the contracted Schrödinger equation (CSE): (i) The cumulant expansion of reduced density matrices (RDMs), and the reduction of the CSE to a connected form. (ii) The comparison of the resulting formalism to coupled cluster theory. (iii) The direct calculation of the two-particle RDM in terms of the cumulants, without explicitly imposing any N-representability conditions. We explore a number of approximations, the most elaborate of which is the 3,5-CSE in which the cumulant expansion is carried out to the three-particle level, and the only approximation made is to equate the four- and five-particle cumulants to zero. Overall the numerical results we present for selected small model systems are found to be disappointing and we provide a critical analysis of the perceived shortcomings of our cumulant approach to the CSE.

The Journal of Chemical Physics•Feb 28, 2003•42 citations
Linear response (LR) is an important tool in the computational chemist's toolbox. It is therefore no surprise that the emergence of quantum computers has led to a quantum version, quantum LR (qLR). However, the current quantum era of near-term intermediary scale quantum (NISQ) computers is dominated by noise, short decoherence times, and slow measurement speed. It is therefore of interest to find approximations that greatly reduce the quantum workload while only slightly impacting the quality of a method. In an effort to achieve this, we approximate the naive qLR with singles and doubles (qLRSD) method by either directly approximating the reduced density matrices (RDMs) or indirectly through their respective reduced density cumulants (RDCs). We present an analysis of the measurement costs behind qLR with RDMs, and report qLR results for model Hydrogen ladder systems; for varying active space sizes of OCS, SeH2_22​, and H2_22​S; and for symmetrically stretched H2_22​O and BeH2_22​. Discouragingly, while approximations to the 4-body RDMs and RDCs seem to produce good results for systems at the equilibrium geometry and for some types of core excitations, they both tend to fail when the system exhibits strong correlation. All approximations to the 3-body RDMs and/or RDCs severely affect the results and cannot be applied.
The prediction of spectral properties via linear response (LR) theory is an important tool in quantum chemistry for understanding photoinduced processes in molecular systems. With the advances of quantum computing, we recently adapted this method for near-term quantum hardware using a truncated active space approximation with orbital rotation, named quantum linear response (qLR). In an effort to reduce the classic cost of this hybrid approach, we here derive and implement a reduced density matrix (RDM) driven approach of qLR. This allows for the calculation of spectral properties of moderately sized molecules with much larger basis sets than so far possible. We report qLR results for benzene and R -methyloxirane with a cc-pVTZ basis set and study the effect of shot noise on the valence and oxygen K-edge absorption spectra of H 2 O in the cc-pVTZ basis.
Linear response (LR) theory is a powerful tool in classic quantum chemistry crucial to understanding photoinduced processes in chemistry and biology. However, performing simulations for large systems and in the case of strong electron correlation remains challenging. Quantum computers are poised to facilitate the simulation of such systems, and recently, a quantum linear response formulation (qLR) was introduced [Kumar et al., J. Chem. Theory Comput. 2023, 19, 9136–9150]. To apply qLR to near-term quantum computers beyond a minimal basis set, we here introduce a resource-efficient qLR theory, using a truncated active-space version of the multiconfigurational self-consistent field LR ansatz. Therein, we investigate eight different near-term qLR formalisms that utilize novel operator transformations that allow the qLR equations to be performed on near-term hardware. Simulating excited state potential energy curves and absorption spectra for various test cases, we identify two promising candidates, dubbed “proj LRSD” and “all-proj LRSD”.
Accurate modeling of the response of molecular systems to an external electromagnetic field is challenging on classical computers, especially in the regime of strong electronic correlation. In this article, we develop a quantum linear response (qLR) theory to calculate molecular response properties on near-term quantum computers. Inspired by the recently developed variants of the quantum counterpart of equation of motion (qEOM) theory, the qLR formalism employs "killer condition" satisfying excitation operator manifolds that offer a number of theoretical advantages along with reduced quantum resource requirements. We also used the qEOM framework in this work to calculate the state-specific response properties. Further, through noiseless quantum simulations, we show that response properties calculated using the qLR approach are more accurate than the ones obtained from the classical coupled-cluster-based linear response models due to the improved quality of the ground-state wave function obtained using the ADAPT-VQE algorithm.
Accurate modeling of the response of molecular systems to an external electromagnetic field is challenging on classical computers, especially in the regime of strong electronic correlation. In this paper, we develop a quantum linear response (qLR) theory to calculate molecular response properties on near-term quantum computers. Inspired by the recently developed variants of the quantum counterpart of equation of motion (qEOM) theory, the qLR formalism employs "killer condition" satisfying excitation operator manifolds that offers a number of theoretical advantages along with reduced quantum resource requirements. We also used the qEOM framework in this work to calculate state-specific response properties. Further, through noise-less quantum simulations, we show that response properties calculated using the qLR approach are more accurate than the ones obtained from the classical coupled-cluster based linear response models due to the improved quality of the ground-state wavefunction obtained using the ADAPT-VQE algorithm.
The promise of quantum computing to circumvent the exponential scaling of quantum chemistry has sparked a race to develop chemistry algorithms for quantum architecture. However, most works neglect the quantum-inherent shot noise, let alone the effect of current noisy devices. Here, we present a comprehensive study of quantum linear response (qLR) theory obtaining spectroscopic properties on simulated fault-tolerant quantum computers and present-day near-term quantum hardware. This work introduces novel metrics to analyze and predict the origins of noise in the quantum algorithm, proposes an Ansatz-based error mitigation technique, and reveals the significant impact of Pauli saving in reducing measurement costs and noise in subspace methods. Our hardware results using up to cc-pVTZ basis set serve as proof of principle for obtaining absorption spectra on quantum hardware in a general approach with the accuracy of classical multi-configurational methods. Importantly, our results exemplify that substantial improvements in hardware error rates and measurement speed are necessary to lift quantum computational chemistry from proof of concept to an actual impact in the field.
In this paper we address a number of topics regarding the contracted Schrödinger equation (CSE): (i) The cumulant expansion of reduced density matrices (RDMs), and the reduction of the CSE to a connected form. (ii) The comparison of the resulting formalism to coupled cluster theory. (iii) The direct calculation of the two-particle RDM in terms of the cumulants, without explicitly imposing any N-representability conditions. We explore a number of approximations, the most elaborate of which is the 3,5-CSE in which the cumulant expansion is carried out to the three-particle level, and the only approximation made is to equate the four- and five-particle cumulants to zero. Overall the numerical results we present for selected small model systems are found to be disappointing and we provide a critical analysis of the perceived shortcomings of our cumulant approach to the CSE.