---
title: Multichannel Dyson Equation Overview
url: https://www.emergentmind.com/topics/multichannel-dyson-equation
type: topic
---

# Multichannel Dyson Equation Overview

Searching arXiv for recent papers on the multichannel Dyson equation and closely related Dyson-equation generalizations.
The multichannel Dyson equation is a matrix generalization of the Dyson equation in which several many-body Green’s functions, or several channel sectors of higher-order Green’s functions, are coupled and solved within a single enlarged propagator. In the formulation introduced for photoemission, the one-body Green’s function is coupled to the \(2h1e\) and \(2e1h\) sectors of the three-body Green’s function, so that quasiparticle poles and satellite structures arise from the same effective problem rather than from a single one-body self-energy alone [2310.01204]. Subsequent developments extended this construction to a broader hierarchy of even- and odd-order Green’s functions, to screened formulations for solids, and to neutral- and double-ionization spectroscopies [2406.12485, 2501.04176, 2603.27329, 2604.01085].

## 1. Definition and conceptual scope

In standard many-body perturbation theory, a single-channel Dyson equation relates one interacting Green’s function to a reference propagator and a self-energy,
\[
G(\omega)=G_0(\omega)+G_0(\omega)\,\Sigma(\omega)\,G(\omega),
\]
or equivalently \(G^{-1}=G_0^{-1}-\Sigma\). In that setting, one usually works with the one-body Green’s function, and all many-body structure beyond the independent-particle poles is encoded indirectly in the frequency dependence of a one-body self-energy.

The multichannel Dyson equation replaces that closure on a single propagator by a coupled problem in an enlarged channel space. A “channel” denotes a sector with definite particle-hole content that shares the same final particle number and therefore the same excitation energies. For final states with \(M=N\pm1\), the one-body Green’s function and the \(2e1h/2h1e\) sectors of the three-body Green’s function have the same poles; for neutral excitations with \(M=N\), the electron-hole channel of the two-body Green’s function and the \(2e2h\) sector of the four-body Green’s function share the same poles [2501.04176].

This common pole structure is the central rationale of the formalism. Instead of integrating higher configurations out into a dynamical kernel, the multichannel Dyson equation keeps them explicitly. In the original photoemission application, this means that quasiparticles and satellites are treated on equal footing because both are eigenmodes of the same coupled problem [2310.01204]. In later neutral-excitation and double-ionization formulations, the same logic is used to represent double excitations, biexcitons, and shake-up structures as explicit multichannel configurations rather than as residual structures hidden in a frequency-dependent kernel [2501.04176, 2604.01085].

A useful classification introduced in the generalized theory is the \((n,s)\)-MCDE notation, where \(n\) denotes the highest-order Green’s function retained and \(s=|M-N|\) the change in particle number. Thus \((3,1)\)-MCDE targets photoemission and inverse photoemission, \((4,0)\)-MCDE targets neutral excitations, and \((4,2)\)-MCDE targets double addition or double ionization [2501.04176, 2604.01085].

## 2. Block-matrix formulation

The multichannel Dyson equation is written for a block Green’s function,
\[
\mathbf G(\omega)=\mathbf G_0(\omega)+\mathbf G_0(\omega)\,\mathbf \Sigma(\omega)\,\mathbf G(\omega),
\]
where \(\mathbf G_0\) is block diagonal in channel space and \(\mathbf \Sigma\) contains both intra-channel and inter-channel blocks.

For the original \((3,1)\)-MCDE, the noninteracting propagator is
\[
G_{03}(\omega)=
\begin{pmatrix}
G_{01}(\omega) & 0\\
0 & G_{03}^{3\mathrm p}(\omega)
\end{pmatrix},
\]
where \(G_{01}\) is the one-body Green’s function and \(G_{03}^{3\mathrm p}\) is the three-body Green’s function restricted to the \(2h1e/2e1h\) sector. The corresponding multichannel self-energy has the block form
\[
\Sigma_3=
\begin{pmatrix}
\Sigma^{1\mathrm p} & \Sigma^{\mathrm c}\\
\tilde\Sigma^{\mathrm c} & \Sigma^{3\mathrm p}
\end{pmatrix},
\]
with \(\Sigma^{\mathrm c}\) and \(\tilde\Sigma^{\mathrm c}\) coupling one-particle and three-particle sectors [2310.01204, 2406.12485].

For neutral excitations, the \((4,0)\)-MCDE couples the electron-hole channel of the two-body Green’s function to the \(2e2h\) channel of the four-body Green’s function. For double ionization, the \((4,2)\)-MCDE couples the particle-particle channel of the two-body Green’s function to the \(3h1e/3e1h\) channel of the four-body Green’s function [2501.04176, 2604.01085].

| Variant | Coupled channels | Principal target |
|---|---|---|
| \((3,1)\)-MCDE | \(1\mathrm p\) and \(2h1e/2e1h\) sectors of \(G_3\) | Photoemission and inverse photoemission |
| \((4,0)\)-MCDE | \(1e1h\) sector of \(G_2\) and \(2e2h\) sector of \(G_4\) | Neutral excitations, biexcitons, double excitations |
| \((4,2)\)-MCDE | Particle-particle \(G_2\) and \(3h1e/3e1h\) sector of \(G_4\) | Double ionization and double addition |

The physical interpretation of the off-diagonal blocks is direct. In the \((3,1)\) case, they convert a one-particle excitation into a three-particle one and back, i.e. a quasiparticle into a shake-up configuration and back. In the \((4,0)\) case, they convert a single electron-hole pair into a double electron-hole configuration and back. In the \((4,2)\) case, they convert a double-ionization quasiparticle into a configuration with an extra electron-hole pair. This explicit conversion mechanism is what makes satellite and multi-excitation physics appear as ordinary poles of the coupled Green’s function rather than as emergent anomalies of a reduced self-energy [2406.12485, 2501.04176, 2604.01085].

## 3. Multichannel self-energy and generated correlations

A distinctive feature of the MCDE literature is the use of a static multichannel self-energy constructed at first order in the bare Coulomb interaction. In the \((3,1)\)-MCDE, the proposed approximation sets
\[
\Sigma^{1\mathrm p}_{i;m}=0,\qquad
\Sigma^{\mathrm c}_{i;mok}=\bar v_{ikom},\qquad
\tilde\Sigma^{\mathrm c}_{ijl;m}=\bar v_{ijlm},
\]
with \(\bar v\) the antisymmetrized Coulomb matrix element, while the three-particle block \(\Sigma^{3\mathrm p}\) is a static first-order expression built from antisymmetrized Coulomb integrals and occupation-number prefactors that separate the \(2e1h\) and \(2h1e\) sectors [2310.01204].

This static construction has two immediate consequences. First, no explicit frequency convolutions appear at the multichannel level. Second, no self-consistency loop is required in the original approximation, because the multichannel self-energy is not taken as a functional of the interacting multichannel Green’s function [2406.12485]. Nevertheless, when the higher channel is downfolded, the effective lower-channel self-energy becomes dynamical. In the screened and unscreened photoemission formulations, the effective one-body correlation self-energy takes the form
\[
\Sigma_c(\omega)=\Sigma_{1p/3p}\,G_{3P}(\omega)\,\Sigma_{3p/1p},
\]
so its frequency dependence is entirely inherited from propagation in the higher channel [2603.27329].

Diagrammatically, the resulting content is substantially richer than the static input might suggest. For the \((3,1)\)-MCDE, the method is complete up to second order in the one-body self-energy, yields all ten third-order proper skeleton diagrams, includes all \(GW\) diagrams to infinite order, and additionally generates ladder and vertex-correction diagrams beyond \(GW\) [2310.01204]. The detailed derivation traces these structures to repeated insertions of the coupling blocks and the three-particle block, with the head of the multichannel Dyson series reproducing the ordinary one-body self-energy expansion [2406.12485].

The same first-order-in-\(v\) strategy is extended systematically to higher-order MCDEs. In the \((4,0)\)-MCDE, the head block is the standard RPAx Bethe–Salpeter kernel in the electron-hole sector, while the four-particle block and coupling blocks add all first-order direct and exchange interactions within and between \(1e1h\) and \(2e2h\) configurations [2501.04176]. In the \((4,2)\)-MCDE, the two-particle head is the pp-RPAx kernel, while the four-particle block and couplings add first-order pairwise Coulomb interactions in the \(3h1e/3e1h\) sector and between that sector and the pp channel [2604.01085].

A screened variant modifies this strategy for extended systems. The screened multichannel Dyson equation keeps the block structure of the \((3,1)\)-MCDE but replaces selected bare Coulomb matrix elements in the three-particle block by statically screened matrix elements \(W_{ijkl}\), with \(W(\mathbf r_1,\mathbf r_2)\) evaluated at \(\omega=0\) from an RPA dielectric function. The motivation is that the unscreened MCDE implicitly screens some, but not all, Coulomb lines, whereas the screened version aims to ensure that all interactions are either explicitly or effectively screened [2603.27329].

## 4. Effective-Hamiltonian formulation

Because the multichannel self-energy is static in these constructions, the MCDE can be rewritten as a resolvent of an effective Hamiltonian. For the \((3,1)\)-MCDE,
\[
G_3(\omega)=\left[\omega I-H^{\mathrm{eff}}\right]^{-1},
\]
with
\[
H^{\mathrm{eff}}=
\begin{pmatrix}
H^{1\mathrm p} & \Sigma^{\mathrm c}\\
\tilde\Sigma^{\mathrm c} & H^{3\mathrm p}
\end{pmatrix},
\]
where the one-particle block is diagonal in the reference orbital energies and the three-particle block is the noninteracting \(2h1e/2e1h\) Hamiltonian plus \(\Sigma^{3\mathrm p}\) [2310.01204]. The screened MCDE has an analogous Hermitian effective Hamiltonian, with the screened three-particle block \(H_{3p}^{\mathrm{scr}}\) replacing the unscreened one [2603.27329].

This representation is structurally important. It makes all poles of the coupled Green’s function eigenvalues of a single static operator. Quasiparticles correspond to eigenvectors with large projection on the lowest channel, while satellites or multiple-excitation states correspond to eigenvectors with large projection on higher channels. The physical spectral function is obtained by projecting the resolvent onto the relevant head block; for photoemission, that is the one-particle component, and for neutral excitations it is the two-particle electron-hole component [2310.01204, 2501.04176].

The effective-Hamiltonian mapping also changes the numerical character of the problem. In the original \((3,1)\)-MCDE, direct diagonalization of the effective Hamiltonian scales as \(N^9\), while iterative methods such as Haydock–Lanczos reduce the scaling to \(N^6\) [2310.01204]. In the screened silicon application, the effective Hamiltonian is solved with the Lanczos–Haydock method to obtain spectral information without full diagonalization [2603.27329]. In the \((4,2)\)-MCDE for double ionization, the dominant cost is the construction of the effective Hamiltonian, which scales as \(N^6\) in the number of electrons or orbitals, and iterative solvers are again proposed as the practical route for large systems [2604.01085].

This suggests a general tradeoff that recurs throughout the literature: the MCDE replaces a nonlinear or highly dynamical reduced problem by a larger but static eigenvalue problem. The computational burden moves from frequency-dependent kernels and self-consistency to enlarged configuration spaces and sparse linear algebra.

## 5. Applications and physical regimes

The first benchmark application was the Hubbard dimer. In the original \((3,1)\)-MCDE, coupling the one-body and three-body channels yields the exact spectral function of the Hubbard dimer at both \(1/4\) and \(1/2\) filling, while \(GW\) fails to reproduce the satellites correctly and yields unphysical extra peaks [2310.01204]. The same work also states that the approximate multichannel self-energy recovers the exact three-body Green’s function \(G_3\) for the dimer at those fillings.

A later application to the extended Hubbard dimer used the MCDE to compute the potential energy surface, spectral functions, and HOMO-LUMO gaps. The reported outcome is that the MCDE gives overall very good results for all properties considered and outperforms both \(GW\) and second Born; in particular, it yields the correct ground-state energy and HOMO-LUMO gap in the dissociation limit, contrary to \(GW\) [2507.12437]. This case is significant because dissociation is precisely the regime in which single-channel perturbative self-energies are known to become qualitatively unreliable.

For realistic solids, the screened MCDE was introduced to describe direct and inverse photoemission in bulk silicon. In that case, \(GW\), the unscreened MCDE, and the second-Born approximation all strongly overestimate the binding energy of the plasmon satellite, placing it around \(25\) eV, whereas the experimental valence plasmon satellite is about \(17\) eV below the valence-band maximum. The screened MCDE captures the correct satellite position and also reproduces the main features of the direct and inverse photoemission spectra [2603.27329]. The same study interprets this as evidence that properly screening all particle-particle and electron-hole interactions is essential in extended systems.

For neutral excitations, the \((4,0)\)-MCDE was applied to a two-level Helium-like model. Standard static BSE@GW does not produce a genuine double-excitation pole, whereas the MCDE yields both single and double excitations and gives excitation energies in good agreement with exact results [2501.04176]. The analysis also shows that the quantitative position of the double excitation is sensitive to the quasiparticle gap used in the higher-channel block, which indicates that the quality of the underlying one-particle input remains important even when the multichannel coupling is correct.

For double-ionization spectroscopies, the \((4,2)\)-MCDE provides the formal extension beyond pp-RPA by coupling the particle-particle two-body Green’s function to the \(3h1e/3e1h\) channels of the four-body Green’s function. In this formulation, double-ionization quasiparticles and satellites are described within the same multichannel framework, and the downfolded two-particle self-energy becomes dynamical through the higher-channel propagator, despite the static first-order input kernel [2604.01085].

Across these applications, a common pattern emerges. Whenever the relevant spectroscopy is strongly influenced by shake-up configurations, multiple electron-hole pairs, or dissociation-induced multiconfigurational character, the explicit retention of higher channels improves qualitative fidelity relative to single-channel approximations. This suggests that the principal strength of the MCDE is not merely numerical but representational: it enlarges the propagating space so that physically important configurations appear as primary degrees of freedom.

## 6. Related meanings, limitations, and extensions

The term “multichannel Dyson equation” is used in at least two related but distinct senses. In the MCDE literature discussed above, it denotes a coupled Dyson equation for different many-body Green’s functions or channel sectors. In another usage, “multichannel” refers to ordinary Dyson equations whose matrix indices already span several coupled physical channels such as momentum, orbital, atom, and spin. The fully self-consistent finite-temperature \(GW\) study of NdNiO\(_2\) is an example: there the Dyson equation is solved for \(G_{pq}(i\omega_n)\) in the space of momentum-dependent spin-orbitals, and the authors explicitly interpret this as a multichannel problem in \(k\)-point, orbital, and spin space. Their homotopy-continuation method is used to track multiple self-consistent Dyson solutions and branch switching in this high-dimensional channel space [2507.00290]. This is mathematically related, but it is not the same formalism as the block-coupled MCDE.

A further related usage appears in non-equilibrium quantum transport with several leads. In a rigorous partition-free formulation, the retarded Green operator acts on the full one-particle Hilbert space of an interacting sample coupled to several non-interacting reservoirs. The resulting Dyson equation is multichannel in the sense that the Green function has sample–lead and lead–lead blocks, while the interaction self-energy is supported only on the sample subspace [1901.07863]. Again, this is a multichannel Dyson equation in a transport-channel basis, not in the explicit many-body-channel sense of the MCDE.

Within the explicit MCDE program, several limitations are repeatedly emphasized. The first is scaling: enlarging the propagator space to \(3p\), \(4p\), or higher channels quickly becomes expensive [2310.01204, 2604.01085]. The second is truncation: a \((3,1)\)-MCDE omits \(5\)-body channels, a \((4,0)\)-MCDE omits \(6p\) channels, and so on, so the method is controlled by channel selection rather than by a small parameter alone [2501.04176]. The third is kernel approximation: most current implementations use static first-order kernels, which are exact only through low perturbative order and may require screened interactions or higher-order corrections for quantitative accuracy in general systems [2603.27329, 2507.12437].

The extension strategy is correspondingly explicit. The general formalism allows odd- and even-order channel towers, dressing of selected Coulomb lines by screened interactions, replacement of Hartree–Fock higher-channel propagators by second Born, \(GW\), or \(T\)-matrix quasiparticle propagators, and higher-channel constructions such as \((5,1)\)-MCDE for double satellites or \((4,2)\)-MCDE for Auger and related spectroscopies [2310.01204, 2501.04176, 2603.27329, 2604.01085]. A plausible implication is that the MCDE should be viewed less as a single approximation than as a general architecture for embedding lower-order observables into a systematically enlarged Green’s-function hierarchy.

Source: https://www.emergentmind.com/topics/multichannel-dyson-equation