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:2511.06623 · Physical Review Research

Quantum algorithm for one-quasiparticle excitations in the thermodynamic limit via cluster-additive block diagonalization

Sumeet, M. Hörmann, K. P. Schmidt

quant-phcond-mat.str-elSubmitted Nov 10, 2025• Announced Nov 11, 2025• 0 citations• Revised Jul 13, 2026 (v2)
v1 · Nov 10, 2025v2 · Jul 13, 2026
arXivHTMLPDF

Abstract

We propose a quantum algorithm for computing one-quasiparticle excitation energies in the thermodynamic limit by combining numerical linked-cluster expansions (NLCEs) and the variational quantum eigensolver (VQE). Our approach uses VQE to block-diagonalize the cluster Hamiltonian through a single-unitary transformation. This unitary is then postprocessed using the projective cluster-additive transformation (PCAT) to ensure cluster additivity, a key requirement for NLCE convergence. We benchmark our method on the transverse-field Ising model (TFIM) in one and two dimensions, and with longitudinal field, computing one-quasiparticle dispersions in the high-field polarized phase. We compare two cost function classes, trace minimization and variance based, demonstrating their effectiveness with the Hamiltonian variational ansatz (HVA). For pure TFIM, $\lceil N/2 \rceil$ layers of HVA suffice: NLCE+VQE matches exact diagonalization. For TFIM with longitudinal field, where parity symmetry breaks and PCAT becomes essential, both $\lceil N/2 \rceil$ and $N$ layers of HVA converge with increasing cluster size, with $N$ layers providing improved accuracy. Our results establish PCAT as a cluster-additive framework that extends variational quantum algorithms to excited-state calculations in the thermodynamic limit via NLCE. While demonstrated with VQE, the PCAT postprocessing approach, which requires only low-energy eigenspace information, applies to any quantum eigenstate preparation method.

Journal reference
Phys. Rev. Research 8, 023358 (2026)
DOI
10.1103/sdpn-s42c
Comments
22 pages, 14 figures

Related papers

Thermodynamic-limit dispersion relations on trapped-ion quantum hardware

Lucas Marti, Sumeet, Stefan Wolf +2 more

We run a numerical linked-cluster expansion with a quantum algorithm (NLCE+QA), computing ground-state energies and one quasi-particle dispersions in the thermodynamic limit using a 20-qubit trapped-ion quantum processing unit (QPU). The NLCE+QA framework extracts thermodynamic-limit properties from small-cluster calculations, making it naturally suited for near-term quantum devices. Projector-based block-diagonalization schemes such as projective cluster-additive transformation (PCAT) are essential to NLCE+QA, and they involve matrix inversion and square root operations that amplify measurement noise. A central question is therefore whether current hardware can provide expectation values that are accurate enough to withstand non-linear classical post-processing. We explore this challenge for the transverse-field Ising model (TFIM) in one dimension, on a ladder geometry, as well as in a longitudinal field in one dimension. For the quantum algorithm, we consider adiabatic state preparation (ASP), as well as a variational quantum eigensolver (VQE) trained on a classical device. The final expectation values are obtained from the QPU, using a novel alternative to the Hadamard test that we name the CX-test. We explore the regimes currently attainable on quantum devices and comment on the improvements needed for quantum computers to achieve results beyond classical reach.

quant-pharXiv:2605.28599•May 28, 2026•0 citations

Hybrid quantum-classical algorithm for the transverse-field Ising model in the thermodynamic limit

Sumeet, Max Hörmann, Kai Phillip Schmidt

We describe a hybrid quantum-classical approach to treat quantum many-body systems in the thermodynamic limit. This is done by combining numerical linked-cluster expansions (NLCE) with the variational quantum eigensolver (VQE). Here, the VQE algorithm is used as a cluster solver within the NLCE. We test our hybrid quantum-classical algorithm ($\mathrm{NLCE}+\mathrm{VQE}$) for the ferromagnetic transverse-field Ising model on the one-dimensional chain and the two-dimensional square lattice. The calculation of ground-state energies on each open cluster demands a modified Hamiltonian variational ansatz for the VQE. One major finding is convergence of $\mathrm{NLCE}+\mathrm{VQE}$ to the conventional NLCE result in the thermodynamic limit when at least $N/2$ circuit repetitions (known as layers) are used in the VQE ansatz for each cluster with $N$ sites. Our approach demonstrates the fruitful connection of techniques known from correlated quantum many-body systems with hybrid algorithms explored on existing quantum-computing devices.

arXiv:2310.07600•Physical review. B./Physical review. B•Oct 15, 2024•9 citations

Quasiparticle pictures and graphs - from perturbative to non-perturbative linked-cluster expansions

Kris Cöster

A fundamental principle in condensed matter physics is the effective description in terms of quasiparticles. The high-energy part of the system is accounted for by renormalized properties of the quasiparticles, providing both, an accurate description and interpretation of the relevant low-energy physics. We consider two classes of methods suited to provide effective quasiparticle descriptions in the thermodynamic limit: (non-perturbative) linked-cluster expansions ((N)LCEs) and continuous unitary transformations (CUTs). We focus specifically on the combination of both methods, providing a different perspective and solutions to existing challenges. LCEs provide high-order series expansions in a perturbation parameter, by combining effective Hamiltonians determined on finite clusters. We introduce a white-graph expansion for the method of perturbative continuous unitary transformations when implemented as an LCE. The essential idea behind an expansion in white graphs is to perform an optimized bookkeeping during the calculation by exploiting the model-independent effective Hamiltonian in second quantization and the associated inherent cluster additivity. This approach is shown to be especially well suited for microscopic models with many coupling constants, since the total number of relevant graphs is drastically reduced. The white-graph expansion is exemplified for a two-dimensional quantum spin model, illustrating its efficiency. In NLCEs, the perturbative treatment of finite clusters is replaced by a numerical exact (block) diagonalization. While LCEs are restricted due to their perturbative nature, the non-perturbative variant is not. In graph-based continuous unitary transformations (gCUTs), the block diagonalization is achieved by a non-perturbative CUT performed on finite clusters. The central objective of this thesis is a modification of the gCUT scheme, allowing to treat two kinds of major issues, denoted by pseudo and genuine decay, occuring for NLCEs due to the non-perturbative treatment of clusters. Indeed, one finds surprising effects caused by the non-perturbative renormalization. In particular, we identify a fundamental challenge for any non-perturbative approach based on finite clusters resulting from the reduced symmetry on graphs, most importantly the breaking of translational symmetry when targeting the properties of excited states. This can be traced back to the appearance of intruder states in the low-energy spectrum, which represent a major obstacle in quasi-degenerate perturbation theory. Here, a generalized notion of cluster additivity is introduced, which is used to formulate an optimized scheme of gCUTs, allowing to solve and to physically understand this major issue. Most remarkably, our improved scheme demands to go beyond the paradigm of using the exact eigenvectors on graphs. We demonstrate that the modified scheme is correct in the non-perturbative regime. Even at quantum criticality, the scheme gives valid results. To determine the critical behavior, one must rely on extrapolation schemes. We introduce a generic approach to extract critical properties from sequences of numerical data which is directly relevant for NLCEs. The scheme is applied to the quantum phase transition between the dimerized and the isotropic spin 1/2 Heisenberg chain. Finally, we investigate a scenario where the gapped quasiparticle excitation is not stable for all momenta, i.e., the one-particle mode merges with the continuum for certain momenta and one observes quasiparticle decay. In this case, intruder states merging with the low-energy spectrum on finite clusters represent genuine physics of the system. Again, a proper renormalization on finite clusters satisfies the generalized cluster additivity. Our adjusted renormalization approach performed on finite clusters does not necessarily lead to a full decoupling of the quasiparticle subspaces and the remaining interactions are part of the effective description. The resulting effective Hamiltonian can be analyzed to describe quasiparticle decay in the thermodynamic limit. The modified scheme is applied to four-leg spin 1/2 Heisenberg ladders, providing insights into the spectral properties. Overall, this thesis presents several important developments for the derivation of effective quasiparticle pictures in quantum spin models via perturbative and non-perturbative LCEs, opening various opportunities for future investigations.

•Dec 1, 2015•2 citations

Numerical linked cluster expansions for quantum quenches in one-dimensional lattices

Krishnanand Mallayya, Marcos Rigol

We discuss the application of numerical linked cluster expansions (NLCEs) to study one dimensional lattice systems in thermal equilibrium and after quantum quenches from thermal equilibrium states. For the former, we calculate observables in the grand canonical ensemble, and for the latter we calculate observables in the diagonal ensemble. When converged, NLCEs provide results in the thermodynamic limit. We use two different NLCEs: a maximally connected expansion introduced in previous works and a site-based expansion. We compare the effectiveness of both NLCEs. The site-based NLCE is found to work best for systems in thermal equilibrium. However, in thermal equilibrium and after quantum quenches, the site-based NLCE can diverge when the maximally connected one converges. We relate this divergence to the exponentially large number of clusters in the site-based NLCE and the behavior of the weights of observables in those clusters. We discuss the effectiveness of resummations to cure the divergence. Our NLCE calculations are compared to exact diagonalization ones in lattices with periodic boundary conditions. NLCEs are found to outperform exact diagonalization in periodic systems for all quantities studied.

arXiv:1701.02758•Physical review. E•Mar 7, 2017•20 citations

Effective quasiparticle conserving Lindbladians in the thermodynamic limit

Lea Lenke, Kai Phillip Schmidt

Open quantum many-body systems are commonly described by Lindblad master equations, yet the treatment of Lindbladian operators in the thermodynamic limit remains a major challenge. We develop a framework for constructing effective quasiparticle-conserving Lindbladian operators directly in the thermodynamic limit. Our approach extends continuous similarity transformations to non-Hermitian open quantum systems and enables the systematic block diagonalization of Lindbladians with respect to the quasiparticle (qp) number. We formulate two complementary methods. The first, projective continuous similarity transformations (pcst++), generalizes perturbative continuous unitary transformations to Lindblad operators so that a linked-cluster expansion allows us to obtain high-order series expansions of the infinite system. The second, deepCST, extends directly evaluated enhanced perturbative continuous unitary transformations by combining the same qp-conserving generator with a perturbative truncation scheme that yields non-perturbative effective Lindbladian operators directly in thermodynamic limit. We apply pcst++ and deepCST to the dissipative transverse-field Ising chain with local dissipation. We focus on the low-Ising regime. We purify the Lindbladian by splitting each spin into two sites. A spin flip then corresponds to two qps. We derive and analyze the effective, qp-conserving Lindbladian in the sectors with zero to two qps. We show how the Ising interaction renormalizes the decay rates of elementary spin-flip excitations and provides the microscopic mechanism for a competition between coherent interactions and dissipation. Our work establishes perturbative and non-perturbative continuous similarity transformations as a versatile tool for deriving effective qp pictures of open quantum many-body systems in the thermodynamic limit.

quant-pharXiv:2609.02725•Sep 3, 2026•0 citations

Projective cluster-additive transformation for quantum lattice models

Max Hörmann, Kai Phillip Schmidt

We construct a projection-based cluster-additive transformation that block-diagonalizes wide classes of lattice Hamiltonians \mathcal{H}=\mathcal{H}_0 +V ℋ = ℋ 0 + V . Its cluster additivity is an essential ingredient to set up perturbative or non-perturbative linked-cluster expansions for degenerate excitation subspaces of \mathcal{H}_0 ℋ 0 . Our transformation generalizes the minimal transformation known amongst others under the names Takahashi’s transformation, Schrieffer-Wolff transformation, des Cloiseaux effective Hamiltonian, canonical van Vleck effective Hamiltonian or two-block orthogonalization method. The effective cluster-additive Hamiltonian and the transformation for a given subspace of \mathcal{H} ℋ , that is adiabatically connected to the eigenspace of \mathcal{H}_0 ℋ 0 with eigenvalue e_0^n e 0 n , solely depends on the eigenspaces of \mathcal{H} ℋ connected to e_0^m e 0 m with e_0^m\leq e_0^n e 0 m ≤ e 0 n . In contrast, other cluster-additive transformations like the multi-block orthogonalization method or perturbative continuous unitary transformations need a larger basis. This can be exploited to implement the transformation efficiently both perturbatively and non-perturbatively. As a benchmark, we perform perturbative and non-perturbative linked-cluster expansions in the low-field ordered phase of the transverse-field Ising model on the square lattice for single spin-flips and two spin-flip bound-states.

arXiv:2303.04774•SciPost Physics•Sep 15, 2023•8 citations
We propose a quantum algorithm for computing one-quasiparticle excitation energies in the thermodynamic limit by combining numerical linked-cluster expansions (NLCEs) and the variational quantum eigensolver (VQE). Our approach uses VQE to block-diagonalize the cluster Hamiltonian through a single-unitary transformation. This unitary is then postprocessed using the projective cluster-additive transformation (PCAT) to ensure cluster additivity, a key requirement for NLCE convergence. We benchmark our method on the transverse-field Ising model (TFIM) in one and two dimensions, and with longitudinal field, computing one-quasiparticle dispersions in the high-field polarized phase. We compare two cost function classes, trace minimization and variance based, demonstrating their effectiveness with the Hamiltonian variational ansatz (HVA). For pure TFIM, ⌈N/2⌉\lceil N/2 \rceil⌈N/2⌉ layers of HVA suffice: NLCE+VQE matches exact diagonalization. For TFIM with longitudinal field, where parity symmetry breaks and PCAT becomes essential, both ⌈N/2⌉\lceil N/2 \rceil⌈N/2⌉ and NNN layers of HVA converge with increasing cluster size, with NNN layers providing improved accuracy. Our results establish PCAT as a cluster-additive framework that extends variational quantum algorithms to excited-state calculations in the thermodynamic limit via NLCE. While demonstrated with VQE, the PCAT postprocessing approach, which requires only low-energy eigenspace information, applies to any quantum eigenstate preparation method.
We run a numerical linked-cluster expansion with a quantum algorithm (NLCE+QA), computing ground-state energies and one quasi-particle dispersions in the thermodynamic limit using a 20-qubit trapped-ion quantum processing unit (QPU). The NLCE+QA framework extracts thermodynamic-limit properties from small-cluster calculations, making it naturally suited for near-term quantum devices. Projector-based block-diagonalization schemes such as projective cluster-additive transformation (PCAT) are essential to NLCE+QA, and they involve matrix inversion and square root operations that amplify measurement noise. A central question is therefore whether current hardware can provide expectation values that are accurate enough to withstand non-linear classical post-processing. We explore this challenge for the transverse-field Ising model (TFIM) in one dimension, on a ladder geometry, as well as in a longitudinal field in one dimension. For the quantum algorithm, we consider adiabatic state preparation (ASP), as well as a variational quantum eigensolver (VQE) trained on a classical device. The final expectation values are obtained from the QPU, using a novel alternative to the Hadamard test that we name the CX-test. We explore the regimes currently attainable on quantum devices and comment on the improvements needed for quantum computers to achieve results beyond classical reach.
We describe a hybrid quantum-classical approach to treat quantum many-body systems in the thermodynamic limit. This is done by combining numerical linked-cluster expansions (NLCE) with the variational quantum eigensolver (VQE). Here, the VQE algorithm is used as a cluster solver within the NLCE. We test our hybrid quantum-classical algorithm (NLCE+VQE\mathrm{NLCE}+\mathrm{VQE}NLCE+VQE) for the ferromagnetic transverse-field Ising model on the one-dimensional chain and the two-dimensional square lattice. The calculation of ground-state energies on each open cluster demands a modified Hamiltonian variational ansatz for the VQE. One major finding is convergence of NLCE+VQE\mathrm{NLCE}+\mathrm{VQE}NLCE+VQE to the conventional NLCE result in the thermodynamic limit when at least N/2N/2N/2 circuit repetitions (known as layers) are used in the VQE ansatz for each cluster with NNN sites. Our approach demonstrates the fruitful connection of techniques known from correlated quantum many-body systems with hybrid algorithms explored on existing quantum-computing devices.
A fundamental principle in condensed matter physics is the effective description in terms of quasiparticles. The high-energy part of the system is accounted for by renormalized properties of the quasiparticles, providing both, an accurate description and interpretation of the relevant low-energy physics. We consider two classes of methods suited to provide effective quasiparticle descriptions in the thermodynamic limit: (non-perturbative) linked-cluster expansions ((N)LCEs) and continuous unitary transformations (CUTs). We focus specifically on the combination of both methods, providing a different perspective and solutions to existing challenges. LCEs provide high-order series expansions in a perturbation parameter, by combining effective Hamiltonians determined on finite clusters. We introduce a white-graph expansion for the method of perturbative continuous unitary transformations when implemented as an LCE. The essential idea behind an expansion in white graphs is to perform an optimized bookkeeping during the calculation by exploiting the model-independent effective Hamiltonian in second quantization and the associated inherent cluster additivity. This approach is shown to be especially well suited for microscopic models with many coupling constants, since the total number of relevant graphs is drastically reduced. The white-graph expansion is exemplified for a two-dimensional quantum spin model, illustrating its efficiency. In NLCEs, the perturbative treatment of finite clusters is replaced by a numerical exact (block) diagonalization. While LCEs are restricted due to their perturbative nature, the non-perturbative variant is not. In graph-based continuous unitary transformations (gCUTs), the block diagonalization is achieved by a non-perturbative CUT performed on finite clusters. The central objective of this thesis is a modification of the gCUT scheme, allowing to treat two kinds of major issues, denoted by pseudo and genuine decay, occuring for NLCEs due to the non-perturbative treatment of clusters. Indeed, one finds surprising effects caused by the non-perturbative renormalization. In particular, we identify a fundamental challenge for any non-perturbative approach based on finite clusters resulting from the reduced symmetry on graphs, most importantly the breaking of translational symmetry when targeting the properties of excited states. This can be traced back to the appearance of intruder states in the low-energy spectrum, which represent a major obstacle in quasi-degenerate perturbation theory. Here, a generalized notion of cluster additivity is introduced, which is used to formulate an optimized scheme of gCUTs, allowing to solve and to physically understand this major issue. Most remarkably, our improved scheme demands to go beyond the paradigm of using the exact eigenvectors on graphs. We demonstrate that the modified scheme is correct in the non-perturbative regime. Even at quantum criticality, the scheme gives valid results. To determine the critical behavior, one must rely on extrapolation schemes. We introduce a generic approach to extract critical properties from sequences of numerical data which is directly relevant for NLCEs. The scheme is applied to the quantum phase transition between the dimerized and the isotropic spin 1/2 Heisenberg chain. Finally, we investigate a scenario where the gapped quasiparticle excitation is not stable for all momenta, i.e., the one-particle mode merges with the continuum for certain momenta and one observes quasiparticle decay. In this case, intruder states merging with the low-energy spectrum on finite clusters represent genuine physics of the system. Again, a proper renormalization on finite clusters satisfies the generalized cluster additivity. Our adjusted renormalization approach performed on finite clusters does not necessarily lead to a full decoupling of the quasiparticle subspaces and the remaining interactions are part of the effective description. The resulting effective Hamiltonian can be analyzed to describe quasiparticle decay in the thermodynamic limit. The modified scheme is applied to four-leg spin 1/2 Heisenberg ladders, providing insights into the spectral properties. Overall, this thesis presents several important developments for the derivation of effective quasiparticle pictures in quantum spin models via perturbative and non-perturbative LCEs, opening various opportunities for future investigations.
We discuss the application of numerical linked cluster expansions (NLCEs) to study one dimensional lattice systems in thermal equilibrium and after quantum quenches from thermal equilibrium states. For the former, we calculate observables in the grand canonical ensemble, and for the latter we calculate observables in the diagonal ensemble. When converged, NLCEs provide results in the thermodynamic limit. We use two different NLCEs: a maximally connected expansion introduced in previous works and a site-based expansion. We compare the effectiveness of both NLCEs. The site-based NLCE is found to work best for systems in thermal equilibrium. However, in thermal equilibrium and after quantum quenches, the site-based NLCE can diverge when the maximally connected one converges. We relate this divergence to the exponentially large number of clusters in the site-based NLCE and the behavior of the weights of observables in those clusters. We discuss the effectiveness of resummations to cure the divergence. Our NLCE calculations are compared to exact diagonalization ones in lattices with periodic boundary conditions. NLCEs are found to outperform exact diagonalization in periodic systems for all quantities studied.
Open quantum many-body systems are commonly described by Lindblad master equations, yet the treatment of Lindbladian operators in the thermodynamic limit remains a major challenge. We develop a framework for constructing effective quasiparticle-conserving Lindbladian operators directly in the thermodynamic limit. Our approach extends continuous similarity transformations to non-Hermitian open quantum systems and enables the systematic block diagonalization of Lindbladians with respect to the quasiparticle (qp) number. We formulate two complementary methods. The first, projective continuous similarity transformations (pcst++), generalizes perturbative continuous unitary transformations to Lindblad operators so that a linked-cluster expansion allows us to obtain high-order series expansions of the infinite system. The second, deepCST, extends directly evaluated enhanced perturbative continuous unitary transformations by combining the same qp-conserving generator with a perturbative truncation scheme that yields non-perturbative effective Lindbladian operators directly in thermodynamic limit. We apply pcst++ and deepCST to the dissipative transverse-field Ising chain with local dissipation. We focus on the low-Ising regime. We purify the Lindbladian by splitting each spin into two sites. A spin flip then corresponds to two qps. We derive and analyze the effective, qp-conserving Lindbladian in the sectors with zero to two qps. We show how the Ising interaction renormalizes the decay rates of elementary spin-flip excitations and provides the microscopic mechanism for a competition between coherent interactions and dissipation. Our work establishes perturbative and non-perturbative continuous similarity transformations as a versatile tool for deriving effective qp pictures of open quantum many-body systems in the thermodynamic limit.
We construct a projection-based cluster-additive transformation that block-diagonalizes wide classes of lattice Hamiltonians \mathcal{H}=\mathcal{H}_0 +V ℋ = ℋ 0 + V . Its cluster additivity is an essential ingredient to set up perturbative or non-perturbative linked-cluster expansions for degenerate excitation subspaces of \mathcal{H}_0 ℋ 0 . Our transformation generalizes the minimal transformation known amongst others under the names Takahashi’s transformation, Schrieffer-Wolff transformation, des Cloiseaux effective Hamiltonian, canonical van Vleck effective Hamiltonian or two-block orthogonalization method. The effective cluster-additive Hamiltonian and the transformation for a given subspace of \mathcal{H} ℋ , that is adiabatically connected to the eigenspace of \mathcal{H}_0 ℋ 0 with eigenvalue e_0^n e 0 n , solely depends on the eigenspaces of \mathcal{H} ℋ connected to e_0^m e 0 m with e_0^m\leq e_0^n e 0 m ≤ e 0 n . In contrast, other cluster-additive transformations like the multi-block orthogonalization method or perturbative continuous unitary transformations need a larger basis. This can be exploited to implement the transformation efficiently both perturbatively and non-perturbatively. As a benchmark, we perform perturbative and non-perturbative linked-cluster expansions in the low-field ordered phase of the transverse-field Ising model on the square lattice for single spin-flips and two spin-flip bound-states.