---
title: Cluster D-TRILEX Framework
url: https://www.emergentmind.com/topics/cluster-extension-of-d-trilex
type: topic
---

# Cluster D-TRILEX Framework

Searching arXiv for the cited D-TRILEX cluster-extension paper and closely related identifiers.
The cluster extension of D-TRILEX is a diagrammatic framework for correlated electronic systems in which the dual TRILEX expansion is built around a cluster reference system rather than a single-site one. In the formulation introduced in "Cluster-diagrammatic D-TRILEX approach to non-local electronic correlations" [2507.06015], the method is designed to combine the exact treatment of short-range correlation effects within the cluster with an efficient diagrammatic description of long-range charge and spin collective fluctuations beyond the cluster. The construction is presented for the Hubbard model, benchmarked on the one-dimensional nano-ring Hubbard model, and analyzed in relation to periodization and translational-symmetry breaking.

## 1. Definition and conceptual setting

D-TRILEX is extended by replacing the single-site reference problem with a cluster impurity problem. The resulting scheme starts from cluster dynamical mean-field theory (CDMFT), but does not stop at the cluster solution: off-diagonal and longer-ranged correlations are generated afterwards by dual diagrammatics. In this sense, the method is a cluster-diagrammatic extension in which the cluster provides the nonperturbative short-distance input and the dual fermion–boson formalism supplies non-local corrections [2507.06015].

The central objective is to retain an exact numerical treatment of intra-cluster physics while restoring, at least partially, information lost in a purely cluster-local construction. The paper formulates this objective in terms of combining "the exact treatment of short-range correlation effects within the cluster" with "an efficient diagrammatic description of the long-range charge and spin collective fluctuations beyond the cluster" [2507.06015]. This suggests that the cluster extension is intended to interpolate between strongly local impurity solvers and explicitly momentum-sensitive many-body approximations.

A distinguishing feature of the construction is that it is organized around three-point vertices rather than four-point vertices. The data explicitly notes that the computational effort "remains far cheaper than handling four-point vertices" [2507.06015]. A plausible implication is that the method targets a compromise between non-local accuracy and computational tractability.

## 2. Cluster reference system

The starting point is the Hubbard model on the original lattice, written as an action
\[
{\cal S}[c^*,c]
=
-\sum_{k\sigma ll'}c_{k\sigma l}^{*}\bigl[(i\nu+\mu)\,\delta_{ll'}-\varepsilon^{ll'}_{k}\bigr]\,c_{k\sigma l'}
\;+\;\frac12\!\sum_{\substack{k,k',q\l_1\!-\!l_4,\sigma\sigma'}}
U^{\,l_1l_2l_3l_4}\,
c_{k\sigma l_1}^{*}c_{k+q,\sigma l_2}^{\phantom{*}}
c_{k'+q,\sigma' l_4}^{*}c_{k'\sigma' l_3}^{\phantom{*}}\,.
\]

Within CDMFT, the lattice is tiled into identical clusters of \(N_{\rm c}\) sites and represented by an impurity action
\[
{\cal S}_{\rm imp}[\bar c^*,\bar c]
= \sum_{\nu,\sigma,ll'}\bar c_{\nu\sigma l}^*\bigl[(i\nu+\mu)\,\delta_{ll'}-\Delta^{ll'}_{\nu}\bigr]\bar c_{\nu\sigma l'}
\;+\;
\tfrac12\!\sum_{\substack{\nu,\nu',\omega\l_1\!-\!l_4,\sigma\sigma'}}
U_{\,l_1l_2l_3l_4}\,
\bar c_{\nu\sigma l_1}^*
\bar c_{\nu+\omega,\sigma l_2}^{\phantom{*}}
\bar c_{\nu'+\omega,\sigma' l_4}^*
\bar c_{\nu'\sigma' l_3}^{\phantom{*}}\,.
\]

Here \(\Delta^{ll'}_{\nu}\) is the hybridization matrix, which "in general has off-diagonal cluster indices" [2507.06015]. The numerical setup then rotates to a basis \(\mathcal R\) that diagonalizes the local part of the cluster Hamiltonian,
\[
\overline\varepsilon_{K} \;=\;\mathcal R\,\varepsilon_K\,\mathcal R^\dagger,
\qquad
\overline\Delta_\nu
\;=\;\mathcal R\,\Delta_\nu\,\mathcal R^\dagger,
\]
and chooses as reference system the diagonal hybridization \(\Delta^{ll}_\nu\). The off-diagonal terms are not solved directly at the impurity level; instead, "off-diagonal terms are then generated later by the dual diagrammatics" [2507.06015].

This separation of roles is one of the defining structural choices of the method. The impurity problem supplies a cluster-resolved but diagonally hybridized starting point, while the subsequent expansion reconstructs the missing non-local structure.

## 3. Dual fermion–boson construction beyond the cluster

The extension beyond the cluster is built by subtracting and adding back the diagonal \(\Delta_\nu\), introducing dual fermions \(\tilde f,\tilde f^*\) and Hubbard–Stratonovich bosons \(\phi^{\varsigma}\) with \(\varsigma={\rm ch},x,y,z\), and decoupling charge and spin channels. The dual action takes the form
\[
{\cal S}_{\rm dual}
=
-\sum_{K\nu,ll'}\tilde f_{K\nu l}^*\bigl[\hat{\cal G}_{K\nu}^{-1}-\delta_{ll'}\,g_{\nu}^{-1}\bigr]_{ll'}\,\tilde f_{K\nu l'}
\;+\;
\tfrac12\sum_{Q\omega,\varsigma}
\phi_{Q\omega}^{\varsigma}\, [\hat W_{Q\omega}^{\varsigma}]^{-1}\,\phi_{-Q,-\omega}^{\varsigma}
\;+\;
{\cal V}[\tilde f,\phi]\,,
\]
where \(g_\nu\) and \(\hat W_{Q\omega}^{\varsigma}\) are the impurity Green’s function and renormalized interaction, and \({\cal V}\) is expressed through the three-point vertices \(\Lambda_{\nu\omega}^{\varsigma}\) [2507.06015].

From the converged cluster impurity problem one computes the three-point vertex
\[
\Lambda^{\varsigma\,l_1l_2;l_3l_4}_{\nu\omega}
\;=\;
\frac{\langle \bar c_{\nu+\omega,l_1}
\bar c^*_{\nu,l_2}
\phi^{\varsigma}_{\omega,l_3l_4}
\rangle_{\rm imp}}
{g_{\nu+\omega}^{l_1l_2}\;g_{\nu}^{l_2l_1}\;
\chi^{\varsigma}_{\omega,l_3l_4}},
\]
together with the cluster susceptibilities \(\chi^{\varsigma}_{\omega}\). These objects determine the dual propagators and dual self-energies.

At the one-particle level, the cluster self-energy is
\[
\Sigma_{\rm c}(K,i\nu)
\;=\;
\bigl[G_{0,\rm c}(K,i\nu)\bigr]^{-1}
\;-\;
\bigl[G_{\rm c}(K,i\nu)\bigr]^{-1},
\]
with \(G_{0,\rm c}^{-1}=(i\nu+\mu)\,\openone-\overline\varepsilon_{K}\). After summing dual diagrams, one obtains a non-local correction \(\Delta\Sigma_{\rm nonloc}(k,i\nu)\), and the full lattice self-energy is written as
\[
\Sigma(k,i\nu)
\;=\;
\Sigma_{\rm c}\bigl(K(k),i\nu\bigr)
\;+\;\Delta\Sigma_{\rm nonloc}(k,i\nu).
\]

The leading D-TRILEX correction has a \(GW\)-like form in each channel \(\varsigma\):
\[
\Delta\Sigma_{\rm nonloc}(k,i\nu)
=\;
-\,T\sum_{q,\omega,\varsigma}
\Lambda^{\varsigma}_{\nu\omega}(k,q)\;
\chi^{\varsigma}(q,i\omega)\;
G^0(k-q,i\nu-i\omega)\;
\Lambda^{\varsigma}_{\nu-\omega,-\omega}(k-q,-q)\,.
\]
The paper also states that the Bethe–Salpeter equation for the full lattice susceptibility is built from the cluster vertex and the lattice dual propagators [2507.06015].

Taken together, these formulas show that the cluster extension preserves the TRILEX-type coupling between fermionic propagation and bosonic collective modes, but embeds that coupling in a cluster-resolved dual formalism.

## 4. Self-consistent computational scheme

The computational workflow is organized into a cluster DMFT stage, a self-consistent dual loop, and a lattice reconstruction stage [2507.06015].

In the cluster DMFT step, one diagonalizes the local cluster Hamiltonian to obtain \(\mathcal R\), transforms the lattice problem to that basis, solves CDMFT with diagonal \(\Delta^{ll}_\nu\), and extracts the impurity quantities \(g_\nu,\;\Sigma^{\rm imp}_\nu,\;\chi^{\varsigma}_{\omega},\;\Lambda^{\varsigma}_{\nu\omega}\).

The self-consistent dual loop proceeds through the sequence listed in the paper:
1. Build bare dual propagators \(\tilde{\cal G}_{K\nu}=\hat G_{K\nu}-g_\nu\) and \(\tilde{\cal W}_{Q\omega}^{\varsigma}=\hat W_{Q\omega}^{\varsigma}-\tfrac12U^{\varsigma}\).
2. Compute the dual polarization \(\tilde\Pi_{Q\omega}^{\varsigma}\) via the two-\(\tilde G\)–two-\(\Lambda\) bubble.
3. Dress the dual interaction through \([\tilde W]^{-1}=[\tilde{\cal W}]^{-1}-\tilde\Pi\).
4. Compute the dual self-energy \(\tilde\Sigma_{K\nu}^{ll'}\) via the single-boson exchange.
5. Update the dual Green’s function according to \(\tilde G^{-1}=\tilde{\cal G}^{-1}-\tilde\Sigma\).
6. Iterate until \(\tilde\Sigma,\tilde\Pi\) converge.

The lattice reconstruction then uses
\[
\overline\Sigma_{K\nu}^{ll'}
=\delta_{ll'}\,\Sigma^{\rm imp}_{\nu,ll}
+\sum_{l_1}\tilde\Sigma_{K\nu}^{ll_1}
\bigl[\openone+g_\nu\tilde\Sigma_{K\nu}\bigr]^{-1}_{l_1l'},
\]
followed by
\[
\overline G_{K\nu}^{-1}=(i\nu+\mu)\openone-\overline\varepsilon_K-\overline\Sigma_{K\nu},
\]
and finally a rotation back to the original site-orbital basis,
\[
G_{K\nu}= \mathcal R^\dagger\,\overline G_{K\nu}\,\mathcal R,
\]
with an analogous transformation for \(\Sigma\).

Convergence is declared when changes in \(\tilde\Sigma\) and \(\Sigma\) between iterations fall below a chosen threshold, with \(10^{-5}\) given as an example. The cost per dual iteration scales as
\[
O(N_k\,N_q\,N_\nu\,N_\omega\,N_c^4),
\]
and the paper emphasizes that this remains much less demanding than schemes requiring four-point vertices [2507.06015].

## 5. Periodization and translational symmetry

A central issue in cluster methods is the reconstruction of lattice quantities from cluster-resolved objects. Once \(\Sigma^{ll'}_{K\nu}\) is known in the reduced Brillouin zone, translational invariance is imposed through
\[
O^{\rm latt}_{k\nu}
\;=\;
\frac1{N_{\rm c}}
\sum_{l,l'\in{\rm cluster}}
e^{-ik\,(r_{l}-r_{l'})}\;O^{ll'}_{K(k)\,\nu},
\]
where \(O\) denotes either \(\Sigma\) or \(G\) [2507.06015].

For a dimer this becomes
\[
O^{\rm latt}_{k\nu}
=\tfrac12\bigl(O^{11}+O^{22}\bigr)
+\Re\,O^{12}\cos(ka)\;+\;\Im\,O^{12}\sin(ka)\,.
\]

The paper contrasts pure CDMFT with cluster D-TRILEX on this point. In pure CDMFT, \(\Sigma^{12}_{\rm inter}=0\), and "the periodized \(\Sigma(k)\) is far from the intra-cluster value" [2507.06015]. In cluster D-TRILEX, the diagrammatic correction \(\Delta\Sigma_{\rm nonloc}\) generates a nonzero inter-cluster \(\Sigma^{12}_{\rm inter}\), and the difference between intra-, inter-, and periodized \(\Sigma\) is reported to be "drastically reduced compared to CDMFT alone" [2507.06015]. The paper interprets this as a partial restoration of translational symmetry by non-local diagrams.

The source of the remaining mismatch is identified explicitly: "the CDMFT impurity problem" is described as "the main source of the translational-symmetry breaking" [2507.06015]. A proposed remedy is an outer self-consistency loop in which the dual-corrected \(\Sigma\) or hybridization \(\Delta\) is fed back into the impurity solver. The stated aim is that this would enforce \(\Sigma^{\rm inter}=\Sigma^{\rm intra}\) exactly and thereby restore full translational invariance. Since the paper presents this as a proposal, it is most accurately described as an intended computational scheme rather than an implemented result.

## 6. Benchmarks and comparative performance

The numerical demonstration is carried out for the one-dimensional nanoring Hubbard model, with chains tiled as \(N_{\rm c}=4,6,8\) using two-site clusters, and all results reported at \(U=2\), \(\beta=10\) [2507.06015]. The principal observable is the imaginary part of the self-energy at the Fermi momentum \(k_F=\pi/2\).

The benchmark comparisons given in the paper can be organized as follows.

| Method | Reported behavior |
|---|---|
| Exact Hirsch–Fye QMC | Shows an insulating divergence \(\Im\Sigma\to +\infty\) as \(\nu\to0\) |
| Parquet DΓA | Fails to diverge at \(N=4\) and remains essentially at the DMFT level |
| Ladder DΓA | Yields a divergence but is quantitatively less accurate |
| Single-site D-TRILEX | Captures an insulating branch, but slightly overshoots at low \(\nu\) |
| Two-site cluster D-TRILEX | Lies almost on top of QMC for both \(N=4\) and \(N=8\) |

To quantify the error at the lowest Matsubara point \(\nu_0\), the paper defines
\[
\delta\Sigma
=\bigl|\Im\Sigma_{\rm method}(k_F,i\nu_0)
-\Im\Sigma_{\rm QMC}(k_F,i\nu_0)\bigr|\,.
\]
The reported values are \(\delta\Sigma\sim0.02\) for cluster D-TRILEX, \(\sim0.2\)–\(0.3\) for ladder DΓA, and \(>0.5\) for parquet DΓA at \(N=4\) [2507.06015].

These benchmarks are used to support two claims made explicitly in the paper: first, that the cluster extension of D-TRILEX "accurately reproduces the electronic self-energy at momenta corresponding to the Fermi energy, in good agreement with the numerically exact quantum Monte Carlo solution" [2507.06015]; second, that it "outperforms significantly more computationally demanding approach based on the parquet approximation" [2507.06015]. Within the scope of the presented data, the most prominent advantage is therefore momentum-resolved self-energy accuracy near the Fermi surface together with a marked reduction of periodization ambiguity.

## 7. Scope, interpretation, and open direction

The method is framed as a cluster-diagrammatic approach that extends CDMFT by non-local fermion–boson diagrams. Its short-range content comes from an exact cluster solution, while its long-range content comes from dual corrections built from three-point vertices and collective charge and spin fluctuations [2507.06015]. This division of labor is the main organizing principle of the formalism.

A frequent misconception in the cluster-method setting is that periodization ambiguities are only a post-processing artifact. The presentation here argues for a more specific diagnosis: the explicit translational-symmetry breaking originates primarily in the CDMFT impurity construction, where intra-cluster bonds are treated exactly and inter-cluster bonds only at the DMFT level, so that \(\Sigma^{\rm inter}=0\) in the reference problem [2507.06015]. The cluster D-TRILEX correction does not eliminate this source completely, but it does generate \(\Sigma^{\rm inter}\neq0\) and brings intra-, inter-, and periodized self-energies closer together.

The open direction identified in the paper is therefore an outer self-consistent reformulation in which the impurity reference itself is updated using the dual-corrected quantities. The paper states that implementing such an outer loop is "a natural next step," while noting the price of recalculating two-particle vertices at each update [2507.06015]. This suggests a research program aimed at a fully translationally invariant cluster-diagrammatic scheme that preserves the same basic fermion–boson architecture.

Source: https://www.emergentmind.com/topics/cluster-extension-of-d-trilex