Skip to content

A many-body diagram is a compact representation of one algebraic term or one class of terms in an expansion of a correlation function, energy, free energy, response, or effective interaction. Lines encode specified two-point functions or interactions. Vertices encode specified couplings. Topology organizes contractions, signs, conservation laws, and repeated structures.

A diagram is not:

  • a literal trajectory followed by a microscopic particle;
  • a complete physical process independent of its mathematical rules;
  • an approximation merely because it has been drawn;
  • uniquely defined without an ordering convention, reference state, and observable;
  • interchangeable across real time, imaginary time, stationary perturbation theory, and QFT.

The basic chain is

ordered operator expansion⟶Wick contractions⟶algebraic terms⟶diagram topologies.\text{ordered operator expansion} \longrightarrow \text{Wick contractions} \longrightarrow \text{algebraic terms} \longrightarrow \text{diagram topologies}.

The value lies in the last step: many long sums differ only by labels while sharing the same topology. Once the rules are fixed, diagrams expose which terms are connected, reducible, repeated, dominant, or required by symmetry.

Helpful background. Wick contractions generate the graphs, and many-body perturbation theory explains the reference and volume counting.

A line acquires a meaning from the Hamiltonian, reference and ordering convention. In the example below a solid directed line is a fermion propagator and a dashed line is a bare instantaneous interaction. Other texts may use the same styles differently.

Use ℏ=kB=1\hbar=k_{\mathrm B}=1 and a normal-state quadratic fermion reference at T>0T>0, with β=1/T\beta=1/T. The reference expectation is the full grand-canonical Fock trace of e−βK0e^{-\beta K_0}, K0=H0−μNK_0=H_0-\mu N. A finite fixed-number projection is generally non-Gaussian and needs modified contraction rules. Define

G0(1,2)=−⟨Tτψ(1)ψ†(2)⟩0.G_0(1,2) = -\left\langle \mathcal T_\tau \psi(1) \psi^\dagger(2) \right\rangle_0.

For a translation-invariant diagonal reference,

G0(k,iωn)=1iωn−ξk,G_0(\mathbf k,i\omega_n) = \frac{1}{ i\omega_n-\xi_{\mathbf k} },

where

ωn=(2n+1)πT\omega_n = (2n+1)\pi T

is a fermionic Matsubara frequency and

ξk=ϵk−μ.\xi_{\mathbf k} = \epsilon_{\mathbf k}-\mu.

This formula defines the solid directed line used below. It is a reference propagator, not yet a statement about the exact interacting spectrum.

Take a two-body interaction

V=12V∑k,k′,q∑σ,σ′vq×ck+q,σ†ck′−q,σ′†ck′,σ′ck,σ.V = \frac{1}{2\mathcal V} \sum_{\mathbf k,\mathbf k',\mathbf q} \sum_{\sigma,\sigma'} v_{\mathbf q} \times c_{\mathbf k+\mathbf q,\sigma}^\dagger c_{\mathbf k'-\mathbf q,\sigma'}^\dagger c_{\mathbf k',\sigma'} c_{\mathbf k,\sigma}.

One may represent vqv_{\mathbf q} by a dashed interaction line joining two density vertices. Alternatively, one may use one antisymmetrized four-fermion vertex. These are different bookkeeping conventions; exchange must not be built into the vertex and then added again as a separate exchange graph.

For an imaginary-time observable OO,

⟨O⟩=⟨TτOIexp⁡ ⁣[−λ∫0βdτ VI(τ)]⟩0⟨Tτexp⁡ ⁣[−λ∫0βdτ VI(τ)]⟩0.\langle O\rangle = \frac{ \left\langle \mathcal T_\tau O_I \exp\!\left[ -\lambda \int_0^\beta d\tau\, V_I(\tau) \right] \right\rangle_0 }{ \left\langle \mathcal T_\tau \exp\!\left[ -\lambda \int_0^\beta d\tau\, V_I(\tau) \right] \right\rangle_0 }.

Expanding the numerator gives

⟨O⟩num=∑n=0∞(−λ)nn!×∫0βdτ1⋯dτn⟨TτOIVI(τ1)⋯VI(τn)⟩0.\langle O\rangle_{\mathrm{num}} = \sum_{n=0}^{\infty} \frac{(-\lambda)^n}{n!} \times \int_0^\beta d\tau_1\cdots d\tau_n \left\langle \mathcal T_\tau O_I V_I(\tau_1)\cdots V_I(\tau_n) \right\rangle_0.

If the reference is Gaussian, Wick’s theorem rewrites each ordered expectation value as a sum over complete contraction patterns. Every contraction contributes a reference two-point function. Every insertion of VV contributes a vertex. A diagram records the connectivity of one contraction pattern after equivalent labelings are grouped.

The denominator contains diagrams with no connection to the external operator OO. These vacuum pieces multiply every numerator contribution. Dividing removes them.

Equivalently, the logarithm of a partition or projection amplitude retains connected vacuum diagrams. Vacuum normalization and connectedness are related but distinct operations:

  • the denominator removes vacuum factors from normalized correlators;
  • selecting connected diagrams removes pieces disconnected from the external insertions;
  • taking a logarithm selects connected vacuum contributions to free energies.

The factor 1/n!1/n! compensates permutations of nn identical interaction insertions. Many permutations generate the same labeled contraction topology and cancel part or all of this factorial. A residual symmetry factor remains when a diagram has automorphisms that leave it unchanged.

One cannot assign a symmetry factor from the picture alone unless the vertex normalization, line types, and whether labels are distinguishable have been specified.

For the working convention:

Graphical object Algebraic meaning
directed solid internal line G0(k,iωn)G_0(\mathbf k,i\omega_n) or a declared dressed GG
dashed interaction line vqv_{\mathbf q}, including its normalization and sign convention
interaction vertex coupling factor plus momentum, frequency, spin, and flavor conservation
external solid leg field insertion defining a one-particle correlator
external source leg density, spin, current, or other operator insertion
closed fermion loop internal momentum-frequency sum and one extra factor −1-1
repeated identical topology symmetry or combinatorial factor

Each internal fermion line contributes its propagator factor. Each independent loop variable, after vertex conservation laws are imposed, contributes a momentum-frequency sum or integral. For one such variable the measure and a line factor look like

T∑ωn∫kG0(k,iωn),T \sum_{\omega_n} \int_{\mathbf k} G_0(\mathbf k,i\omega_n),

only as part of the complete product dictated by the graph. Here

∫k≡∫BZddk(2π)d\int_{\mathbf k} \equiv \int_{\mathrm{BZ}} \frac{d^d k}{(2\pi)^d}

for a lattice Brillouin zone, or a regulated continuum integral.

For a translation-invariant instantaneous two-body vertex,

k1+k2=k3+k4\mathbf k_1+\mathbf k_2 = \mathbf k_3+\mathbf k_4

and

ω1+ω2=ω3+ω4.\omega_1+\omega_2 = \omega_3+\omega_4.

On a crystal, momentum conservation is modulo a reciprocal-lattice vector, allowing Umklapp processes. Spin or flavor is conserved or transformed according to the interaction tensor. Drawing an unlabeled vertex does not remove these Kronecker or Dirac delta functions; it suppresses them visually.

For a number-conserving fermion propagator, an arrow records the orientation of fermionic contraction and conserved charge flow. In a hole interpretation, the same line segment can be read relative to a filled reference. In Nambu notation, line and arrow conventions change because normal and anomalous propagators are combined.

An arrow is not a classical path through space.

A particle-hole bubble with two density insertions and internal lines labeled k plus q and k

Two density insertions join the fermion propagators G0(k+q)G_0(k+q) and G0(k)G_0(k). Here k=(k,iωn)k=(\mathbf k,i\omega_n) and q=(q,iνm)q=(\mathbf q,i\nu_m). There is one independent internal momentum-frequency sum, even though the graph has two internal lines. The solid arrows follow the fermion contractions; the external arrows show transferred momentum and frequency.

If every interaction insertion carries λ\lambda, a diagram with VdV_d interaction vertices is nominally order

λVd.\lambda^{V_d}.

That statement is only the first layer. The actual size also depends on:

  • powers of volume from loop sums and line normalization;
  • spin, flavor, or component sums;
  • energy denominators or low-frequency propagators;
  • small external momentum or frequency;
  • long-range interaction kernels;
  • occupation factors and Fermi-surface geometry;
  • symmetry factors and cancellations;
  • logarithms generated by broad scale intervals.

A high-order diagram can dominate a lower-order one in an infrared-enhanced channel. Conversely, a nominally allowed diagram can vanish by symmetry, kinematics, or antisymmetry.

If a closed matter loop sums over NN flavors, it can contribute a factor NN. A vertex may carry 1/N1/N. Counting loops and vertices then organizes diagrams by powers of 1/N1/N.

Large-N and Saddle-Point Methods Preview explains the model scaling and fluctuation hierarchy. The diagram by itself does not reveal which NN is large.

Each independent internal momentum sum contributes a power of volume in finite normalization; each translation-invariant interaction matrix element carries compensating inverse powers. Connected vacuum or energy diagrams should combine to an extensive result when infrared and ultraviolet behavior is regular.

Disconnected products can carry higher powers of volume. Their cancellation belongs to the linked-cluster structure developed in Perturbation Theory in Many-Body Systems.

Fermionic signs come from anticommuting operators into the order needed to form contractions.

Each closed fermion loop contributes an additional factor

−1.-1.

For density nq=∑kσckσ†ck+q,σn_{\mathbf q}=\sum_{\mathbf k\sigma}c_{\mathbf k\sigma}^\dagger c_{\mathbf k+\mathbf q,\sigma}, define the connected Matsubara susceptibility per volume by

χ0(q,τ)=1V⟨Tτδnq(τ)δn−q(0)⟩0.\chi_0(\mathbf q,\tau)=\frac1{\mathcal V} \langle\mathcal T_\tau\delta n_{\mathbf q}(\tau)\delta n_{-\mathbf q}(0)\rangle_0.

Its regular retarded continuation is the physical susceptibility for a source K′=−fqn−qK'=-f_{\mathbf q}n_{-\mathbf q}, with χR=+iθ(t)⟨[nq(t),n−q]⟩/V\chi^R=+i\theta(t)\langle[n_{\mathbf q}(t),n_{-\mathbf q}]\rangle/\mathcal V. Wick contraction gives

χ0(q,iνm)=−gT∑ωn∫kG0(k+q,iωn+iνm)×G0(k,iωn),\chi_0( \mathbf q, i\nu_m ) = -g T \sum_{\omega_n} \int_{\mathbf k} G_0( \mathbf k+\mathbf q, i\omega_n+i\nu_m ) \times G_0( \mathbf k, i\omega_n ),

where gg is a spin or flavor degeneracy when all components have the same propagator.

The loop minus sign gives a positive static susceptibility for the stable ideal gas. The polarization for a potential-energy source +Un+Un is Π0=−χ0\Pi_0=-\chi_0, as in RPA. This response-definition sign is separate from the fermion-loop sign.

Interchanging two external or internal fermionic labels introduces an antisymmetry sign. One may handle exchange by:

  • drawing direct and exchange contractions separately with their signs; or
  • using antisymmetrized vertices that combine them.

Doing both double counts exchange.

The full sign can also depend on:

  • the ordering of external fermionic operators;
  • anomalous Nambu conventions;
  • how an antisymmetrized matrix element is defined;
  • permutation of identical external legs;
  • contour-branch indices in nonequilibrium theory.

“Minus one per loop” is a rule within a complete convention, not a substitute for deriving that convention.

Connected, Reducible, and Skeleton Diagrams

Section titled “Connected, Reducible, and Skeleton Diagrams”

Several topological classifications answer different questions.

Class Definition Typical role
disconnected graph has more than one connected component cancels from normalized connected correlators or exponentiates in ZZ
connected every part is linked to every other part contributes to connected correlators
one-particle reducible cutting one internal single-particle line disconnects the graph generated by repeated self-energy insertion
one-particle irreducible cannot be disconnected by cutting one internal single-particle line belongs to the proper self-energy or 1PI vertex
amputated external propagator factors have been removed yields a vertex function, which need not be irreducible
skeleton built from dressed lines with no explicit self-energy insertion subgraph used in self-consistent functionals

Connected does not imply one-particle irreducible. A chain of two proper self-energy insertions on one propagator is connected but one-particle reducible.

Two-particle reducibility is channel dependent

Section titled “Two-particle reducibility is channel dependent”

A four-point graph can be reducible by cutting two lines in:

  • the direct particle–hole channel;
  • the crossed particle–hole channel;
  • the particle–particle channel.

The same graph can be irreducible in one channel and reducible in another. Bethe–Salpeter and parquet methods organize this channel dependence. Their full equations lie beyond this preview.

The proper self-energy Σ\Sigma is the sum of one-particle-irreducible two-point insertions in a declared expansion.

The Dyson equation is

G=G0+G0ΣG.G = G_0 + G_0\Sigma G.

In a translation-invariant scalar notation,

G−1(k,iωn)=G0−1(k,iωn)−Σ(k,iωn).G^{-1}( \mathbf k, i\omega_n ) = G_0^{-1}( \mathbf k, i\omega_n ) -\Sigma( \mathbf k, i\omega_n ).

Iterating once gives

G=G0+G0ΣG0+G0ΣG0ΣG0+⋯ .\begin{aligned} G ={}& G_0 + G_0\Sigma G_0 + G_0\Sigma G_0\Sigma G_0 + \cdots. \end{aligned}

Thus Σ\Sigma is not merely an energy shift. It is the irreducible kernel whose repeated insertion dresses propagation.

For a two-body fermion interaction, first-order self-energy terms commonly include:

  • a Hartree direct term determined by an average density;
  • a Fock exchange term produced by antisymmetry.

For an instantaneous interaction and a stationary normal reference, these terms can be frequency independent. Hartree–Fock Approximation derives these fields variationally.

Higher-order self-energies can depend on momentum and frequency. After a justified analytic continuation, their real and imaginary parts affect dispersion and damping. Green Functions in Many-Body QM develops the pole and Lehmann interpretation. Spectral Functions explains how the self-energy is converted into quasiparticle positions, residues, linewidths, and line shapes.

If G0G_0 already contains a Hartree–Fock potential, the perturbation must include a counterterm that subtracts it. Otherwise the first-order self-energy is counted once in G0G_0 and again as an insertion.

This is the diagrammatic form of the repartitioning rule:

H=(Hfree+Uref)+(Hint−Uref).H = \left( H_{\mathrm{free}}+U_{\mathrm{ref}} \right) + \left( H_{\mathrm{int}}-U_{\mathrm{ref}} \right).

A bubble is a closed loop made from two single-particle propagators joined by two operator or interaction vertices. In a density channel, it describes the independent propagation of a particle–hole pair between two density insertions.

The Matsubara sum gives

χ0(q,iνm)=g∫knF(ξk+q)−nF(ξk)iνm−ξk+q+ξk.\begin{aligned} \chi_0( \mathbf q, i\nu_m ) &= g \int_{\mathbf k} \frac{ n_{\mathrm F}( \xi_{\mathbf k+\mathbf q} ) -n_{\mathrm F}( \xi_{\mathbf k} ) }{ i\nu_m -\xi_{\mathbf k+\mathbf q} +\xi_{\mathbf k} }. \end{aligned}

For coincident energies and νm=0\nu_m=0, take the divided-difference limit: the integrand is −nF′(ξ)>0-n_{\mathrm F}'(\xi)>0. For exactly coincident energies and νm≠0\nu_m\ne0, it is zero. In particular, the conserved total number has a static Matsubara fluctuation term at m=0m=0 but no nonzero-frequency retarded commutator response. Such exactly zero-energy equilibrium terms must be separated before continuing the regular dynamic part; they cannot be recovered by blindly substituting into a zero denominator.

One bubble contains:

  • Pauli occupation factors;
  • the independent particle–hole continuum;
  • momentum and energy transfer from the external source;
  • spin or flavor multiplicity;
  • the reference dispersion.

It does not by itself contain:

  • screening to all orders;
  • self-energy broadening;
  • vertex corrections;
  • multiple particle–hole pairs;
  • a collective pole outside the independent continuum.

For repulsive vv, the induced Hartree potential opposes the physical chemical-potential source:

δn=χ0(f−vδn),χRPA=χ01+vχ0.\delta n=\chi_0(f-v\delta n),\qquad \chi_{\mathrm{RPA}}=\frac{\chi_0}{1+v\chi_0}.

Equivalently, ΠRPA=Π0/(1−vΠ0)\Pi_{\mathrm{RPA}}=\Pi_0/(1-v\Pi_0) with Π0=−χ0\Pi_0=-\chi_0. RPA derives this feedback, screening and collective poles. A bubble is the two-propagator building block; the self-consistent chain is the additional approximation.

Not every two-line loop is the same channel

Section titled “Not every two-line loop is the same channel”

Pairing diagrams contain two fermion lines but organize particle–particle propagation rather than particle–hole propagation. Nambu bubbles can include anomalous propagators. A visual loop is not enough to identify the channel; arrows, external labels, and momentum routing matter.

A vertex correction changes the coupling between an external source or interaction and the propagating particles. It can encode:

  • exchange between a particle and a hole;
  • repeated scattering in one channel;
  • screening or local-field effects;
  • coupling to collective modes;
  • constraints required by a conservation law.

Let Γμ\Gamma^\mu be the full vertex for a conserved U(1)U(1) current. A Ward identity has the structural form

qμΓμ(k+q,k)=G−1(k+q)−G−1(k),q_\mu \Gamma^\mu(k+q,k) = G^{-1}(k+q) -G^{-1}(k),

with components, factors of charge, ii, and metric conventions fixed by the real- or imaginary-time formulation.

This relation says that self-energy and vertex corrections are not independent decorations. Changing GG changes the vertex structure required by charge conservation.

Replacing G0G_0 by a dressed GG inside a bubble while keeping the bare density vertex can violate a Ward identity. It may be a useful approximation in a specified regime, but it is not automatically conserving.

The required correction depends on how the self-energy was generated. There is no universal vertex patch that repairs every dressed propagator.

A bare expansion uses:

  • fixed G0G_0 lines;
  • fixed bare interactions;
  • all diagrams through a declared order.

Its order counting is transparent, but a low-order truncation may poorly describe shifted poles or collective screening.

One may solve a Dyson equation and use GG in internal lines. This resums infinitely many bare self-energy insertions. If the original insertion diagrams are then retained explicitly, they are double counted.

A skeleton expansion uses dressed propagators but retains only diagrams without explicit self-energy insertion subgraphs. The self-energy is obtained self-consistently from those skeletons.

This can improve consistency, but:

  • self-consistency does not guarantee accuracy;
  • multiple solutions can occur;
  • vertex corrections must still be matched;
  • nominal bare perturbative order is no longer the only ordering variable.

In a conserving construction, one builds a functional Φ[G]\Phi[G] from selected closed skeleton diagrams and defines

Σ=δΦδG\Sigma = \frac{\delta\Phi}{\delta G}

up to convention-dependent signs and factors. Further derivatives generate compatible vertex kernels.

This Baym–Kadanoff or Luttinger–Ward logic explains why a coherent family of diagrams can preserve conservation laws better than an arbitrary list of individually plausible graphs. Its existence, uniqueness, and practical validity are nontrivial and belong to advanced many-body field theory.

Goldstone, Matsubara, and Real-Time Diagrams

Section titled “Goldstone, Matsubara, and Real-Time Diagrams”

Stationary many-body perturbation theory often uses Goldstone diagrams:

  • time flows in one direction;
  • particle and hole lines are defined relative to a reference determinant;
  • intermediate states carry explicit energy denominators;
  • each time ordering can be displayed separately.

They are natural for ground-state energies and wavefunction amplitudes.

Equilibrium finite-temperature diagrams use imaginary time:

  • fermionic frequencies are (2n+1)πT(2n+1)\pi T;
  • bosonic frequencies are 2mπT2m\pi T;
  • internal frequencies are summed;
  • time orderings are combined into propagators;
  • analytic continuation may later connect to real-frequency functions.

One Matsubara or Feynman graph can combine several Goldstone time orderings. Their numerical factors and denominators should not be transferred by visual analogy.

Zero-temperature time-ordered propagators use real frequencies with causal prescriptions. Retarded, advanced, lesser, greater, and time-ordered functions are distinct. A time-ordered self-energy cannot be inserted blindly into a retarded Dyson equation without the correct continuation or contour relation.

Driven or transient systems often require a forward-and-backward time contour. Every line and vertex then carries a contour index or Keldysh component. Real-Time Thermal Dynamics Preview explains the initial-state and contour motivation; full Schwinger–Keldysh rules require a dedicated nonequilibrium field-theory treatment.

A Matsubara function is known initially at discrete imaginary frequencies. For a function with the required analyticity and spectral behavior, the retarded quantity is obtained by continuation

iνm⟶ω+i0+.i\nu_m \longrightarrow \omega+i0^+.

Several cautions matter:

  • perform internal Matsubara sums and derive the analytic function before substituting;
  • preserve branch cuts and thresholds, not only poles;
  • continue the complete conserving combination when cancellations matter;
  • numerical continuation from noisy data is ill conditioned;
  • a finite broadening η\eta is not automatically a physical lifetime.

Retarded and Advanced Green Functions explains causal boundary values in the one-particle setting. Green Functions in Many-Body QM develops the normal many-body retarded and Matsubara spectral relation, including the ill-conditioned numerical-continuation warning.

Diagrams expose divergences; they do not remove them.

Continuum loops can diverge at large internal momentum or frequency. A calculation must state:

  • regulator;
  • bare parameters;
  • matching or renormalization prescription;
  • counterterms;
  • which observable is cutoff independent.

Normal ordering can remove selected vacuum contractions but does not renormalize every local product.

Gapless lines, Fermi surfaces, Goldstone modes, or long-range interactions can enhance low-energy regions. Look for:

ln⁡ ⁣(ΛE),E−α,ln⁡L,T−γ.\ln\!\left( \frac{\Lambda}{E} \right), \qquad E^{-\alpha}, \qquad \ln L, \qquad T^{-\gamma}.

If a family of diagrams has the same enhanced scaling at all orders, fixed-order truncation is not controlled. One must identify and resum the relevant channel.

For response functions,

lim⁡ω→0lim⁡q→0\lim_{\omega\to0} \lim_{\mathbf q\to0}

can differ from

lim⁡q→0lim⁡ω→0.\lim_{\mathbf q\to0} \lim_{\omega\to0}.

The diagram formula does not choose the physical limit. The measurement or thermodynamic question does.

  1. Evaluate the Matsubara bubble sum. Starting from

    χ0(q,iνm)=−gT∑ωn∫k1iωn+iνm−ξk+q1iωn−ξk,\chi_0(\mathbf q,i\nu_m) = -gT \sum_{\omega_n} \int_{\mathbf k} \frac{1}{ i\omega_n+i\nu_m-\xi_{\mathbf k+\mathbf q} } \frac{1}{ i\omega_n-\xi_{\mathbf k} },

    use the standard fermionic sum to obtain the occupation-number form.

Solution

Use

T∑ωn1(iωn−a)(iωn+iνm−b)=nF(a)−nF(b)iνm+a−b.T \sum_{\omega_n} \frac{1}{ (i\omega_n-a) (i\omega_n+i\nu_m-b) } = \frac{ n_{\mathrm F}(a)-n_{\mathrm F}(b) }{ i\nu_m+a-b }.

Set

a=ξk,b=ξk+q.a=\xi_{\mathbf k}, \qquad b=\xi_{\mathbf k+\mathbf q}.

Including the closed-loop minus sign gives

χ0(q,iνm)=g∫knF(ξk+q)−nF(ξk)iνm−ξk+q+ξk.\chi_0( \mathbf q, i\nu_m ) = g \int_{\mathbf k} \frac{ n_{\mathrm F}( \xi_{\mathbf k+\mathbf q} ) -n_{\mathrm F}( \xi_{\mathbf k} ) }{ i\nu_m -\xi_{\mathbf k+\mathbf q} +\xi_{\mathbf k} }.

This is the susceptibility for −fn-fn. The potential-energy polarization is its negative. When a=ba=b and νm=0\nu_m=0, the sum is nF′(a)n_{\mathrm F}^\prime(a), so the loop minus sign gives −nF′(a)>0-n_{\mathrm F}^\prime(a)>0.

  1. Classify topology. A two-point graph consists of an external propagator, then a proper self-energy block, another propagator, a second proper self-energy block, and a final external propagator. Is the whole graph connected? Is it one-particle irreducible? Should it be included in Σ\Sigma?
Solution

The graph is connected because every piece belongs to one component between the two external insertions.

It is one-particle reducible: cutting the propagator between the two self-energy blocks disconnects the graph. Therefore the full chain does not belong to the proper self-energy. It is generated automatically by iterating Dyson’s equation with each proper block included in Σ\Sigma.

Putting the whole chain into Σ\Sigma and also solving Dyson’s equation would double count it.

  1. Bare versus dressed double counting. Suppose

    G=G0+G0Σ(1)G0+O(λ2).G = G_0+G_0\Sigma^{(1)}G_0+O(\lambda^2).

    A calculation uses GG on every internal line and also adds every explicit first-order self-energy insertion graph. Show why an O(λ)O(\lambda) term is counted twice.

Solution

Replacing an internal bare line by GG already gives

G=G0+G0Σ(1)G0+O(λ2).G = G_0 + G_0\Sigma^{(1)}G_0 + O(\lambda^2).

Thus the dressed-line diagram contains the bare topology plus all versions in which one internal line carries the first-order insertion. Adding those insertion graphs explicitly contributes the same

G0Σ(1)G0G_0\Sigma^{(1)}G_0

terms again. A skeleton scheme avoids this by using dressed lines and omitting explicit self-energy insertion subgraphs.

  1. Ward-identity limit. Use the structural Ward identity

    qμΓμ(k+q,k)=G−1(k+q)−G−1(k)q_\mu\Gamma^\mu(k+q,k) = G^{-1}(k+q)-G^{-1}(k)

    to explain qualitatively why a frequency-dependent self-energy generally requires a corrected density vertex in the small-frequency limit.

Solution

After the appropriate analytic continuation, take the dynamic small-frequency limit at zero spatial momentum. This is distinct from the strictly static thermal derivative. Then

G−1(iω+iν)−G−1(iω)≃iν∂G−1∂(iω).G^{-1}(i\omega+i\nu)-G^{-1}(i\omega) \simeq i\nu \frac{\partial G^{-1}}{ \partial(i\omega) }.

More explicitly,

G−1(iω+iν)−G−1(iω)iν⟶∂G−1∂(iω).\frac{ G^{-1}(i\omega+i\nu)-G^{-1}(i\omega) }{ i\nu } \longrightarrow \frac{\partial G^{-1}}{ \partial(i\omega) }.

Since

G−1=G0−1−Σ,G^{-1} = G_0^{-1}-\Sigma,

the derivative contains

−∂Σ∂(iω).-\frac{\partial\Sigma}{ \partial(i\omega) }.

The density vertex must contain the corresponding correction. Using a frequency-dependent dressed GG with an unchanged bare vertex generally fails this relation.