We present a numerical, non-perturbative framework for constructing effective Hamiltonians that describe the dynamics of low-energy degrees of freedom within a restricted Hilbert space in quantum many-body systems. The approach is based on block diagonalization guided by a minimal-deformation principle imposed within a selected target sector. The formulation is designed to remain compatible with the numerical linked-cluster expansion. For gapped systems, the relation between minimal deformation and cluster additivity requires careful treatment when excited eigenstates contain finite admixtures of the ground state. After establishing a cluster-additive basis that reproduces the Hörmann-Schmidt construction, the minimal-deformation criterion uniquely determines the effective Hamiltonian within each excitation sector. The same criterion also provides a practical numerical procedure for selecting relevant low-energy eigenstates, including regimes characterized by strong level mixing and avoided crossings. The framework is illustrated using two spin models: the one-dimensional transverse-field Ising model as a benchmark and the two-dimensional Shastry-Sutherland model with Dzyaloshinskii-Moriya interactions, relevant to SrCu _2 2 (BO _3 3 ) _2 2 . In both cases, the resulting effective Hamiltonians accurately capture the excitation dynamics and the associated band structures.
We present a numerical, non-perturbative framework for constructing effective Hamiltonians that describe the dynamics of low-energy degrees of freedom within a restricted Hilbert space in quantum many-body systems. The approach is based on block diagonalization guided by a minimal-deformation principle imposed within a selected target sector. The formulation is designed to remain compatible with the numerical linked-cluster expansion. For gapped systems, the relation between minimal deformation and cluster additivity requires careful treatment when excited eigenstates contain finite admixtures of the ground state. After establishing a cluster-additive basis that reproduces the Hörmann-Schmidt construction, the minimal-deformation criterion uniquely determines the effective Hamiltonian within each excitation sector. The same criterion also provides a practical numerical procedure for selecting relevant low-energy eigenstates, including regimes characterized by strong level mixing and avoided crossings. The framework is illustrated using two spin models: the one-dimensional transverse-field Ising model as a benchmark and the two-dimensional Shastry-Sutherland model with Dzyaloshinskii-Moriya interactions, relevant to SrCu _2 2 (BO _3 3 ) _2 2 . In both cases, the resulting effective Hamiltonians accurately capture the excitation dynamics and the associated band structures.
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.
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⌉ 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⌉ 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.
In this work ground and excited electronic states of Heisenberg cluster models, in the form of con- figuration interaction many-body wave functions, are characterized within the spin-adapted Graph- ical Unitary Group Approach framework, and relying on a novel combined unitary and symmetric group approach. Finite-size cluster models of well-defined point-group symmetry and of general local-spin Slocal ≥ 21 are presented, including J1–J2 triangular and tetrahedral clusters, which are often used to describe magnetic interactions in biological and bio-mimetic polynuclear transition metal clusters with unique catalytic activity, such as nitrogen fixation and photosynthesis. We show that a unique block-diagonal structure of the underlying Hamiltonian matrix in the spin-adapted basis emerges when an optimal lattice site ordering is chosen that reflects the internal symmetries of the model investigated. The block-diagonal structure is bound to the commutation relations between cumulative spin operators and the Hamiltonian operator, that in turn depend on the geometry of the cluster investigated. The many-body basis transformation, in the form of the orbital/site reorder- ing, exposes such commutation relations. These commutation relations represent a rigorous and formal demonstration of the block-diagonal structure in Hamiltonian matrices and the compression of the corresponding spin-adapted many-body wave functions. As a direct consequence of the block- diagonal structure of the Hamiltonian matrix it is possible to selectively optimize electronic excited states without the overhead of calculating the lower-energy states by simply relying on the initial ansatz for the targeted wave function. Additionally, more compact many-body wave functions are obtained. In extreme cases, electronic states are precisely described by a single configuration state function, despite the curse of dimensionality of the corresponding Hilbert space. These findings are crucial in the electronic structure theory framework, for they offer a conceptual route towards wave functions of reduced multi-reference character, that can be optimized more easily by approximated eigensolvers and are of more facile physical interpretation. They open the way to study larger ab initio and model Hamiltonians of increasingly larger number of correlated electrons, while keeping the computational costs at their lowest.
We demonstrate that a numerical linked cluster expansion method is a powerful tool to calculate quantum dynamics. We calculate the dynamics of the magnetization and spin correlations in the two-dimensional transverse field Ising and XXZ models evolved from a product state. Such dynamics are directly probed in ongoing experiments in ultracold atoms, molecules, and ions. We show that a numerical linked cluster expansion gives dramatically more accurate results at short-to-moderate times than exact diagonalization, and simultaneously requires fewer computational resources. More specifically, the cluster expansion frequently produces more accurate results than an exact diagonalization calculation that would require 105--1010 more computational operations and memory.
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 propose a framework to design concurrently a frustration-free quantum many-body Hamiltonian and its numerically exact ground states on a sufficiently large finite-size cluster in one and two dimensions using an elementary matrix product state (MPS) representation. Our approach strategically chooses a local cluster Hamiltonian, which is arranged to overlap with neighboring clusters on a designed lattice. The frustration-free Hamiltonian is given as the sum of the cluster Hamiltonians by ensuring that there exists a state that has its local submanifolds as the lowest-energy eigenstate of every cluster. The key to find such a solution is a systematic protocol, which projects out excited states on every cluster using MPS and effectively entangles the cluster states. This algorithm not only offers an exact solution to general frustration-free Hamiltonian, but even provide the excited eigenstates by setting a proper target. It has several advantages, including the ability to achieve the calculation of states at nearly equal cost in one and two dimensions including those belonging to gapless or long-range entangled classes, flexibility in designing Hamiltonians unbiasedly across various forms of models, and numerically feasible validation of exactness through energy calculations. It enables the exploration of exact phase boundaries and the analysis of even a spatially nonuniform random system, providing platforms for quantum simulations and benchmarks.