A mathematical derivation of zero-temperature 2D superconductivity from microscopic Bardeen-Cooper-Schrieffer model
J. Magnen, J. Unterberger

TL;DR
This paper rigorously derives the zero-temperature 2D superconducting state from the microscopic BCS model, demonstrating the transition from a Fermi liquid to a Cooper pair condensate and analyzing the associated Goldstone boson.
Contribution
It provides a mathematically rigorous derivation of 2D superconductivity from the BCS model at zero temperature, including symmetry breaking and Goldstone boson analysis.
Findings
Transition energy gap $ o ext{mass gap}$ formula derived
Goldstone boson characterized by non-linear sigma-model
Ladder diagrams dominate infra-red behavior
Abstract
Starting from H. Fr\"ohlich's second-quantized Hamiltonian for a -dimensional electron gas in interaction with lattice phonons describing the quantum vibrations of a metal, we present a rigorous mathematical derivation of the superconducting state, following the principles laid out originally in 1957 by J. Bardeen, L. Cooper and J. Schrieffer. As in the series of papers written on the subject in the 90es, of which the present paper is a continuation, the representation of ions as a uniform charge background allows for a -dimensional fermionic quantum-field theoretic reformulation of the model at equilibrium. For simplicity, we restrict in this article to dimensions and zero temperature, and disregard effects due to electromagnetic interactions. Under these assumptions, we prove transition from a Fermi liquid state to a superconducting state made up of Cooper pairs of…
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsPhysics of Superconductivity and Magnetism · Quantum and electron transport phenomena · Cold Atom Physics and Bose-Einstein Condensates
A mathematical derivation of zero-temperature 2D superconductivity from
microscopic Bardeen-Cooper-Schrieffer model
Jacques Magnena and Jérémie Unterbergerb
aCentre de Physique Théorique,111Laboratoire associé au CNRS UMR 7644 Ecole Polytechnique,
91128 Palaiseau Cedex, France
bInstitut Elie Cartan,222Laboratoire associé au CNRS UMR 7502. J. Unterberger acknowledges the support of the ANR, via the ANR project ANR-16-CE40-0020-01. Université de Lorraine,
B.P. 239, F – 54506 Vandœuvre-lès-Nancy Cedex, France
Starting from H. Fröhlich’s second-quantized Hamiltonian for a -dimensional electron gas in interaction with lattice phonons describing the quantum vibrations of a metal, we present a rigorous mathematical derivation of the superconducting state, following the principles laid out originally in 1957 by J. Bardeen, L. Cooper and J. Schrieffer. As in the series of papers written on the subject in the 90es [22], [25], [24], [23], [56], [17], [18], of which the present paper is a continuation, the representation of ions as a uniform charge background allows for a -dimensional fermionic quantum-field theoretic reformulation of the model at equilibrium. For simplicity, we restrict in this article to dimensions and zero temperature, and disregard effects due to electromagnetic interactions. Under these assumptions, we prove transition from a Fermi liquid state to a superconducting state made up of Cooper pairs of electrons at an energy level equal to the mass gap, expressed in terms of the Debye frequency , electron mass and coupling constant . The dynamical -symmetry breaking produces at energies lower than the energy gap a Goldstone boson, a non-massive particle described by an effective -dimensional non-linear sigma-model, whose parameters and correlations are computed. The proof relies on a mixture of general concepts and tools (multi-scale cluster expansions, Ward identities), adapted to this quantum many-body problem with its extended infra-red singularity located on the Fermi circle, and a specific -expansion giving the leading diagrams at intermediate energies. Ladder diagrams are proved to provide the leading behavior in the infra-red limit, in agreement with mean-field theory predictions.
Some insights about expected extensions of our method to a rigorous study of real-world, low-temperature superconductivity are provided.
Keywords: BCS theory, constructive field theory, renormalization, cluster expansions, non-linear sigma-model, Goldstone boson, low-temperature superconductivity, Cooper pairs, Ward identities, Bethe-Salpeter kernel, Fermi liquids.
Mathematics Subject Classification (2010): 81T08, 81V70, 82D55.
Contents
-
0.1 The Bardeen-Cooper-Schrieffer model of electrons and phonons
-
1.4 Bubble diagrams and fermion four-point function in Coooper channel
-
2.3.4 Momentum-decoupling expansion and displacement of external legs
-
3.4 Bosonic fixed-point procedure and low-energy multi-scale bounds
-
3.8 Proof of Theorem 5: local transverse behavior of the Goldstone boson
Index of notations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Index of notations
0 Introduction
0.1 The Bardeen-Cooper-Schrieffer model of electrons and phonons
The Bardeen-Cooper-Schrieffer (BCS) theory has proved extremely successful at predicting the main characteristic features of conventional, low-temperature superconductivity of metals such as tin, lead, aluminium… see e.g. textbooks [65], [12], [67], [68]. The original objective of the Nobel prizes in their ground-breaking paper [3] was to account for the formation of a condensate of Cooper pairs, bound states made of pairs of electrons which behave like bosons, and are in particular responsible for the Meissner effect 333to which is traditionally attributed the spectacular magnetic levitation. However, disorder-induced vortex trapping supercurrents are more directly responsible for this effect in real materials [5].. The explanation for this bound-state is the existence of an effective attractive force between pairs of electrons due to the quantum oscillations of the ionic lattice. This remarkable effect can be easily understood starting from H. Fröhlich’s [30, 31] second-quantized Hamiltonian , depending on a constant ,
[TABLE]
where Einstein’s implicit summation convention is used for spin states . This Hamiltonian involves:
• an electron field describing a spin , mass , fermionic particle with dispersion relation , where is a chemical potential regulating the density;
• a phonon field, namely, a bosonic quasi-particle describing the oscillations of the lattice,
[TABLE]
with dispersion relation , where is a smoothened Heaviside function. The superconducting material is assumed to take up a roughly cubic volume () with , implying a discrete sum over momenta rougly multiples of . The sum is really a finite sum over momenta such that or equivalently , where is roughly equal to the mean spacing between ions. The frequency is called Debye frequency. We shall eventually restrict to space dimension (possible extension to is briefly discussed in section 4), but let us keep arbitrary till then.
In practice, only electrons with momentum close to the Fermi circle
[TABLE]
defined by the vanishing of the energy participate in the interactions. We assume here that
[TABLE]
a condition satisfied in usual materials (see discussion in [26], §37), and consider only couplings of the phonon field to electrons with momentum such that
[TABLE]
i.e. to electrons with momenta in an annulus of radius around the Fermi circle. The reader may read with profit the classical book by A. Fetter and J. Walecka [26] for: the second quantization formalism applied to interacting many-particle systems (Chapter 1); phonons interacting with electrons (Chapter 12), see in particular §45 where a value for is derived in terms of the bulk compressibility of the ionic background seen as an elastic medium; fundamental properties of superconductors, and a semi-rigorous explanation of these using the above model (0.1), see Chapter 13, following arguments due to Bardeen, Cooper and Schrieffer and detailed computations done by L. P. Gor’kov using finite temperature Green’s functions [39]. For a more mathematical, modern presentation of 3D BCS theory, see the series of papers by R. Frank, C. Hainzl, E. Hamza, B. Schlein, R. Seiringer and J. P. Solovej [44], [28], [42], [43] where the Ginzburg-Landau theory is derived in an -small vicinity of the critical temperature by minimizing in the semi-classical limit a functional introduced by A. J. Leggett [51] (called BCS functional) built out of quasi-free trial states.
At zero temperature, the lowest energy state of the system is obtained by filling the energy levels up to , i.e. summing over all such that . The self-consistent quadratic approximation of introduced by Bogoliubov [7, 8] involves elementary excitations known as quasi-particles, with a minimum energy known as the energy gap, see e.g. [2], §6.4. A phenomenologically important parameter is the density of states at the Fermi surface (compare to [26], p. 333), at temperature , N(0)\equiv N(T=0)=\frac{2}{2\pi}k_{F}\frac{dk}{d\epsilon}\Big{|}_{\epsilon=\mu}=\frac{m}{\pi\hbar^{2}}. Fetter and Walecka find the following value for the energy gap in three dimensions:
[TABLE]
This energy gap may also be interpreted as a critical temperature , where is the Boltzmann constant. In 2D, however, the critical temperature is known to be zero because of the celebrated Mermin-Wagner theorem, which forbids continuous symmetry-breaking for any . At , the theory is effectively -dimensional, as we shall presently see, so the Mermin-Wagner theorem does not apply.
In this article, we use the well-known equivalence (shown by going over to interaction picture and using Wick’s theorem) of the quantum model with a functional integral representation in terms of Grassmann fields , living in a -dimensional space. The supplementary coordinate plays formally the rôle of an imaginary time, since dynamics are retrieved (at least formally) by letting . Space-time points are denoted by . We consider here only vacuum expectation values at zero temperature. Let us consider the infinite volume limit , so that becomes a continuous momentum variable; and integrate out the phonon field (see [26], §46). Then the ground-state expectation \Big{\langle}\Big{(}\prod_{i=1}^{n}\bar{\psi}_{\sigma_{i}}(\boldsymbol{x}_{i})\Big{)}\Big{(}\prod_{i^{\prime}=1}^{n}\psi_{\sigma_{i^{\prime}}}(\boldsymbol{x}^{\prime}_{i^{\prime}})\Big{\rangle} becomes the Grassmann integral
[TABLE]
which may be expressed in terms of the normalized Grassmann measure
[TABLE]
where:
• are now fields living in -dimensional space-time with coordinate ;
• is a Grassmann Gaussian measure with covariance kernel given by the inverse of the quadratic form ,
[TABLE]
for fields with Fourier support satisfying the cut-off condition
[TABLE]
see (0.5). A precise definition of the model will be given only at the end of section 1 (see Definition 1.1). It is enough to say here that the measure involves only Fourier scales , i.e. momenta such that , \Big{|}|\boldsymbol{p}|-p_{F}\Big{|}=O(2^{-j_{D}}p_{F}), where and
[TABLE]
by (0.4);
• and the interaction , formally defined in Fourier coordinates as
[TABLE]
for a general, spin-neutral two-body potential written in second-quantized form using the positive electron density operator , may be chosen in this context in the form
[TABLE]
where is (the Fourier transform of) a static, rotation-invariant, spin-neutral two-body potential, say, with high enough infra-red cut-off (see discussion in section 4), and (mean interparticle spacing, divided by the Bohr radius) is assumed to be , corresponding to a high-density (also called: degenerate) regime of the electron gas, see [26], §3. Up to inessential issues regarding cut-offs (with but not quite), (0.1) reduces exactly to this model if one sets and lets
[TABLE]
The well-known and essential observation is that the above kernel is , bounded above by 1, and in average, under our cut-off conditions. In practice, we shall simply replace by , implying an attractive -interaction between electrons. We still denote by the corresponding interaction,
[TABLE]
and argue in section 4 that a small enough, short-ranged two-body potential can be added to without altering the general conclusions, provided the overall effect near the Fermi sphere remains that of an attractive potential.
Let us mention the following result by W. Kohn and J. M. Luttinger [47, 52]: in dimensions, they proved that for essentially any arbitrary, even purely repulsive, rotation-invariant interaction, the scattering amplitude for pairs of quasi-particles of opposite momenta had poles on the Fermi surface, implying the possibility of creation of a superconducting bound state in large enough odd angular momentum sectors. This is shown by analyzing the sign of non-zero angular momentum second-order contributions such as those of Fig. 1.7.2, which however vanish for . On the other hand, in dimensions, J. Feldman, H. Knörrer, R. Sinclair and E. Trubowitz [20] proved that third-order contributions, namely, triangle diagrams
p_{1}$$\xi_{1}$$p_{2}$$\xi_{2}$$\xi_{3}\!=\!\xi_{4}$$p_{3}$$p_{4}
Fig. 0.1. Triangle diagram.
created superconducting bound states in the angular momentum sector . Hence one can expect that our conclusions extend to very general, not purely attractive interactions.
0.2 Constructive approaches for superconductivity
Constructive methods (actually, multi-scale cluster expansions), see e.g. [1, 19, 33, 34, 21, 35, 36, 54, 55, 60, 61, 70, 71], consist in implementing rigorously Wilson’s renormalization group ideas developed to study small perturbations of Gaussian models. In the case at hand, correlation functions are computed by averaging w.r. to a Grassmann Gaussian measure perturbed by a quartic interaction , with small enough. The bare covariance is infra-red singular on the set . In principle, bare parameters of the model, , become running coupling constants through the renormalization procedure. In our case, is not renormalized (see below why), and the model is directly rewritten in terms of its renormalized parameters and
[TABLE]
The first step consists in splitting the covariance of the measure according to the distance to the singularity; namely, one rewrites as sums of independent fields, , (), whose covariance is supported on the set |p^{0}|+\frac{p^{*}_{F}}{\mu^{*}}\Big{|}|\boldsymbol{p}|-p^{*}_{F}\Big{|}\approx 2^{-j}\mu. Splitting each field of each vertex into its components, and splitting accordingly perturbative graphs by letting scale indices grow from the top to the bottom (see e.g. Fig. 2.3.1), one gets the following picture: high-momentum diagrams of lowest scale (i.e. with all covariance indices ) are quasi-local w.r. to low-momentum diagrams of highest scale . For small enough, this allows the sum of all perturbative graphs of lowest scale to be resummed into effective corrections to low-momentum vertices, amounting to a scale-by-scale renormalization of parameters. Graphs of a given scale are resummed using a cluster expansion; as well-known, for fermionic theories, the exponential of the interaction may be expanded to infinity, see e.g. [55], [61]. The above scheme works provided effective corrections can be shown to remain small; in particular, if the running coupling constant remains . In the present case, however, the running coupling constant becomes large around some transition scale (logarithm of the inverse of the energy gap) called , implying that the perturbation around the Grassmann Gaussian measure is not pertinent any more. This is interpreted as the formation of a bound state made up of Cooper pairs. Instead of merely relying on perturbation theory, the idea is therefore to resum explicitly a class of four-point diagrams, forming the Bethe-Salpeter kernel, which contribute to this non-perturbative effect. The sum – mathematically, a kernel denoted , one of the components of a two-by-two matrix-valued kernel –, may be interpreted as the two-point function of the bound state. The leading behavior of this kernel may be captured by looking simply at the geometric series of ladder (bubble) diagrams, which can be explicitly computed. These two-point functions turn out to be the only non-massive (i.e. long-range) ones, hence they give the main contribution to the theory at large scale, an effective bosonic theory in the same class as the * non-linear sigma model*. The underlying non-massive boson may be called Goldstone boson, by reference to the general theory of continuous symmetry-breaking (see discussion of Regime II below), and the kernel Goldstone boson propagator, the other non-vanishing component, , being massive. A full understanding of the model at all scales can be obtained by considering ”mixed” diagrams featuring both fermionic propagators and vertices, and some bound-state two-point functions. Contrary to fermionic vertices, though, -kernels must not be systematically expanded, but rather (as for most non-massive bosonic models, see e.g. [21]) through some careful cluster expansion.
Let us emphasize one specific aspect of the present model, which ultimately explains why non-perturbative effects can be dealt with at an analytical level. The small parameter here is the expected mean value of the interaction, restricted to some scale and integrated over a scaled box with sides of length , up to spurious logarithmic corrections, see (2.23), (2.26), independently of the choice of the scale defined for typical momenta as or , which reflects the fact that the theory near the Fermi sphere is just renormalizable (independently of space dimension). Note that this very fact is actually not straightforward, and more easily established in dimensions than for (see brief discussion in section 4). Therefore we restrict to in the sequel. -point functions within a given scale will be reexpressed in terms of a geometric-like entire series in the non-dimensional parameter
[TABLE]
hence converge provided
[TABLE]
(2.26) is based on the ”sector-counting proposition”, see Proposition 2.5, which may be rephrased as follows (see [23]): The theory near the Fermi surface (in this context, Fermi circle) may be with remarkable accuracy reformulated after a scale-dependent rescaling as a large vector model with action
[TABLE]
, where the fields , () are fermion fields restricted to a given momentum angular sector indexed by and rescaled in such a way that . The momentum-scale dependent number increases exponentially as one gets nearer to the Fermi circle defined by . Expanding the interaction and using Wick’s theorem yields Feynman diagrams made up of vertices and fermion loops, each in a given angular sector, which may be thought of as a ”color”. Performing the sum over colors produces for generic vertices a factor per fermion loop, and per vertex. Alternatively, following a fermion loop, one can prove that there is at most one sum over sectors per vertex in a given diagram. Then dominant diagrams in an -expansion are chains of bubbles (see §1.4), as confirmed perturbatively by a Feynman diagram expansion; the simplicity of the theory is due to the fact that chains of bubbles make up a geometric series which can be resummed explicitly. The idea of the ”-expansion” is old and has been used in many different contexts; we refer the reader to [48] and [46] respectively for a rigorous analysis of the -component Gross-Neveu and non-linear sigma model. Instead of -component vectors, one sometimes also considers large -matrices, either in the context of random matrices or Schrödinger operators [6], [58] – in connection to two-dimensional gravity [14], since leading terms are then planar diagrams – or as an approximation to gauge theories, following a seminal paper by ’t Hooft, see [66] or [29], chap. 7 for a review of two-dimensional quantum-field theory models in this limit. R. Gurau, V. Rivasseau et al. have also applied these ideas to tensor field models, see [40], [41].
Under Hypothesis (0.18), the following scenario – in accordance with the original idea by Bardeen, Cooper and Schrieffer, but going beyond mean-field regime predictions, which are valid only for an infinitesimal interaction, physically, in an ideally degenerate regime, see below (0.5) – was explored in the 90’es by various theoretical physicists, including one of the authors of the present work, see articles by M. Disertori, J. Feldman, J. Magnen, V. Rivasseau, E. Trubowitz [22], [25], [24], [23], [56], [16],[17], [18], [15]. Notations are as follows: energy and momentum scales are labeled by an integer index ranging from to ; typical energies , resp. transverse momenta of scale are , resp. . The highest scale is defined in agreement with the above cut-off hypotheses (0.4), (0.5). Physical parameters, in particular, the electron mass, the coupling constant and the Fermi radius (or, equivalently, the chemical potential) are renormalized à la Wilson, defining scale parameters . Three energy regimes were singled out:
- (i)
(Regime I, high-energy regime) At high enough energy, i.e. for small enough, no bound states can form, and electrons are still in their normal, Fermi liquid phase, where they behave essentially like free fermions, see [16]. The scale-by-scale renormalization of the model à la Wilson yields the following flow for the coupling constant,
[TABLE]
where is the amputated Cooper pair bubble diagram of §1.4 (see Fig. 1.4.3), and is the evaluation at zero external momentum of computed for internal momenta of scale , see §1.4. The difference is approximately scale independent (which is essentially tantamount to saying that the theory is just renormalizable) and . As long as remains , the solution of the flow equation is well approximated by
[TABLE]
or equivalently,
[TABLE]
Just a few scales above a transition scale defined by , the renormalized coupling constant becomes large. Thus the above perturbative regime analysis breaks down. Since (as shown in §1.4) , the above condition holds for
[TABLE]
This defines an energy level for the transition,
[TABLE]
Let us add here some necessary precisions. Because one renormalizes in the vicinity of the Fermi circle, which is an extended singularity, it is easy to see that leading corrections to the vertex are ”pinched” diagrams of the form
\psi_{3}$$\psi_{4}$$p_{3}$$p_{1}$$\psi_{1}$$\psi_{2}
Fig. 0.2. Bethe-Salpeter kernel.
(see Fig. 1.7.5 for details), with external momenta near the Fermi circle, depending only on the relative angle . Two-particle irreducible diagrams of this type with Cooper pair external structure, i.e. with , form the so-called Bethe-Salpeter kernel, which may be summed into a geometric series. Now, the bare theory has an ultra-local, angle-independent vertex , and Cooper pair bubbles lead to angle-independent, -wave corrections. However, more complicated diagrams with extra vertices lead to angle-dependent corrections, that may be seen as effective -wave, -wave, etc. effective vertices. To leading order these do not interfere, because they are orthogonal Fourier modes. Hence we get instead of (0.20) an infinite series of flows,
[TABLE]
compare with (0.20), and
[TABLE]
The difference between (0.25) and (0.26) is that whereas , Thus coupling constants other than increase much more slowly, with leading term at scale bounded by the sum , whereas . This is not sufficient to conclude (see discussion in §1.7, and around eq. (3.37)), but strongly hints at an -wave superconducting behavior.
- (ii)
(Regime II, near transition scale)
In regimes II and III, that is, around and below the gap energy, the phonon field , initially centered around 0, is seen by semi-perturbative arguments to choose a more favorable position on a circle of radius , which breaks the fermionic charge symmetry \left(\begin{array}[]{c}\psi(\xi)\\ \bar{\psi}(\xi)\end{array}\right)\mapsto\left(\begin{array}[]{c}e^{{\rm i}\alpha}\psi(\xi)\\ \bar{\psi}(\xi)e^{-{\rm i}\alpha}\end{array}\right). At a semi-rigorous level, this can be seen using a Hubbard-Stratonovich transform,
[TABLE]
{\mathbb{\Gamma}}(\xi)\equiv\left(\begin{array}[]{cc}0&\Gamma^{*}(\xi)\\ \Gamma(\xi)&0\end{array}\right), where is the probability measure on complex-valued fields defined by
[TABLE]
corresponding to a quadratic action . This formula is very much related to the original BCS interaction in , see (0.1), but now the random potential has become an off-diagonal Hermitian matrix , coupling Cooper pairs or . Fermionic charge symmetries translate into rotations of the -field. The fact that the values of the field concentrate statistically on a small neighborhood of a circle is by itself a good argument to try to split locally into the sum of a tangential component (parallel to the circle) and a transversal (orthogonal) component ; from simpler models with the same symmetries, see e.g. [49], §13.3.1 for an example inspired by the Glashow-Salam-Weinberg model of electroweak interactions, it is understood by elementary computations that the tangential, resp. transversal components should be massless, resp. massive. The tangential component is traditionally called Goldstone boson by reference to the well-known Goldstone’s theorem [38] stating that for every spontaneously broken continuous symmetry, a given theory must contain a massless particle. Here we find as effective theory a non-linear sigma model, see e.g. [57], Chapter 11 and §13.3. Applying Gaussian integrations by parts on the measure (0.27), one easily sees that -point functions of the Goldstone boson are directly related to -point functions of electron pairs in the Cooper pair channel, or . Integrating instead w.r. to the fermions yields now a bosonic theory in terms of , with total (i.e., quadratic part + interaction) action
[TABLE]
up to cut-offs, where is the covariance kernel of the fermions. Using the identity \log\det({\mathrm{Id}}-A)={\mathrm{Tr}}\log({\mathrm{Id}}-A)=-{\mathrm{Tr}}\Big{(}A+\frac{A^{2}}{2}+\frac{A^{3}}{3}+\ldots\Big{)}, noting that odd powers do not contribute to the trace, assuming to be constant, and dividing by the volume, one obtains the so-called effective potential , , for which an explicit formula was obtained in [22], section 4,
[TABLE]
with an ultra-violet cut-off in at . (The exact correspondence is in the notations of [22], p. 48, with ). By construction, is due to the ladder diagrams obtained by resumming Cooper pair bubble diagrams, see §1.4. Searching for the minimum argmin of gives an implicit equation, known as the gap equation, which coincides with the one found in mean-field theory , i.e. starting from the Bogoliubov-De Gennes Hamiltonian, see e.g. [2], §6.4, and the one given in our §1.5. The solution of the gap equation is
[TABLE]
where is as in (0.24) and (0.6).
As was already well understood from the -asymptotic expansion in previous work on the subject, dominant contributions due to ladder diagrams are actually the only divergent ones, thus explaining in particular the essential fact that only electron pairs in the Cooper channel contribute to the superconductive phase.
Allowing to fluctuate, one obtains an action for tangential fluctuations which is roughly in Fourier,
[TABLE]
where and
[TABLE]
(Fermi velocity) is a velocity, see e.g. [23].
Although , the zero-temperature theory is effectively a -dimensional theory. Hence Mermin-Wagner’s argument does not apply, and infinitesimal bulk or boundary terms are in principle enough to imply a symmetry-breaking in some direction , by which we mean that
[TABLE]
- (iii)
(Regime III, low-energy regime) Around the transition scale, the Goldstone boson favors a position on a circle of radius , thus conferring the electron an effective squared mass which dominates the kinetic energy term. Thus electrons are ”hooked up” to a fixed scale : Regime III may be thought of from the point of view of unpaired fermions as a single scale extending from to . Dominant low-momentum contributions to Feynman diagrams come therefore exclusively from the Goldstone boson. Now, Ward identities associated to the (broken) number symmetry prove (more or less as in the case of diagrams with low-momentum photons in quantum electrodynamics) that such diagrams vanish up to error terms coming from cut-off effects, which decrease exponentially as one lowers the energy scale. Thus, below an energy level
[TABLE]
the effective coupling constant of diagrams is once again small enough to sum series produced by the cluster expansion. Contrary to individual fermions, the Goldstone boson (or equivalently, electron pairs in the Cooper channel) remains non-massive in the Euclidean infra-red defined not by the vicinity to the Fermi circle, but by the vicinity to [math].
Unfortunately, the writing of the above program stopped at the end of the 90’es somewhere between Regime I and Regime II for lack of a sound mathematical proof, due to technical difficulties coming mainly from the necessity of integrating simultaneously the fermion fields and the Goldstone boson field . At each scale , one had to distinguish ”small-field”, perturbative regions, from ”large-field” regions, which could be handled only through large deviation arguments. In theory, the idea was to split the -field into the sum of a ”fast” perturbation field with scale index , and of a ”slow” background field with scale index , giving a local orientation and making it possible to define a ”tangential” and a ”transversal” component. In the small field region, is small, and minimizes to a good approximation the effective potential computed with the contributions of the scales . The main difficulty was to prove that a given region was ”small-field” with high probability.
Also, both the implementation of Ward identities and the use of the -expansion are awkward in this representation, requiring partial integrations of the bosonic field in order to go back to a purely fermionic representation. Finally, integrating out fermions in order to deal with the infra-red behavior of the boson field naturally leads to functional determinants such as (0.29) which are delicate to define properly and handle in combination with cluster decompositions (see the above cited paper [46] by C. Kopper for a successful result in this direction for a scalar intermediate boson field).
Our strategy here is different; let us describe it briefly, highlighting differences with previous attempts. We introduce by hand a fixed symmetry-breaking term (”adding and subtracting” it) with a given, arbitrary but fixed, orientation (see (0.34)), and a module computed by a fixed-point argument; this spares us the trouble of having to define a local orientation around every point. Thus the transversal, resp. tangential, direction is parallel , resp. perpendicular , to . We carefully avoid introducing the Goldstone field in the first place, producing instead through an explicit expansion an effective, non-local interaction kernel in the perpendicular direction, called Goldstone boson propagator, which is obtained by resumming a geometric series made up of alternating Cooper pair bubble diagrams and Bethe-Salpeter kernels, themselves sums of two-particle irreducible diagrams as in Fig. 0.2; the -kernel may be interpreted as the propagator of a bosonic particle which is never introduced. Decomposing fields into angular sectors, we are able to prove that the two-point function of Cooper pairs diverge only in the -wave. Then – taking the -wave projection – we write down in terms of an exact, non-perturbative version of the gap equation mentioned in the discussion of the neighborhood of the transition scale (see Regime II above). The actual value of the gap is shown to be close to the mean-field one (0.31) by a fixed-point argument. Cooper pair two-point functions are shown to be proportional to the -kernel, itself roughly inverse of the conjectural action associated to tangential fluctuations of the -field (again, see Regime II).
Proceeding this way, we fill this important gap in the literature, thus hopefully laying the foundational stones of a rigorous mathematical analysis of 2D and 3D low-temperature superconductivity from first principles, going beyond Ginzburg-Landau theory. We also conjecture that our methods will help understand related models featuring a dynamical symmetry-breaking transition in a large regime, like 3D Anderson’s model in its delocalized phase, or the integer quantum Hall effect.
0.3 Our results
Instead of singling out three regimes, we proceed as follows.
A. (fermionic regime, or high-momentum theory) We choose a value for such that (0.23) holds, and consider the fermionic theory directly with effective parameters , , , namely, we modify the Grassmann Gaussian measure of eq. (0.9) by multiplying it with the exponential of a quadratic weight,
[TABLE]
where , and {\mathbb{\Gamma}}(\theta)=\left(\begin{array}[]{cc}0&\Gamma^{*}\\ \Gamma&0\end{array}\right)\equiv\left(\begin{array}[]{cc}0&\Gamma_{\phi}e^{-{\rm i}\theta}\\ \Gamma_{\phi}e^{{\rm i}\theta}&0\end{array}\right), is a constant off-diagonal, Hermitian matrix; see Definition 1.1 for a precise definition. Up to cut-offs, the covariance kernel of is the inverse of the quadratic form
[TABLE]
compare with (0.9), where is the effective (or renormalized) dispersion relation. This extra term is compensated by subtracting it to the interaction, which takes the form of a bare, scale counterterm added to the interaction defined in (0.15), namely,
[TABLE]
and
[TABLE]
Because of the ultra-violet, scale cut-off, the resulting measure
[TABLE]
depends on , on , on and on . In particular – as will be shown –, it induces an infinite-volume symmetry-breaking precisely in the direction . Then we study the renormalization flow from scale to scale . Fermions are shown to become massive at energy scales , making it possible to integrate out fermions with energy scales as if they were of the same scale, provided they do not form Cooper pairs with low transfer momentum.
The considerable difference with the above discussion of Regime I is that, should one start from , then the number symmetry would imply that is not renormalized: the fact that in average is a non-perturbative effect. Proceeding that way, the symmetry-breaking would be read only indirectly through the apparition of a large, renormalized coupling constant depending on the momentum angular sectors of its four legs, favoring in Cooper channels. These, in turn, might be subtracted by changing the value of . Our approach looks much simpler.
Let us remark at this point that another (maybe physically more natural) procedure would have led to an equivalent effective, large-scale behavior. Namely, we could have chosen a lattice ultra-violet cut-off with lattice constant on a large volume span, with inverse covariance kernel {\cal A}(\psi,\bar{\psi}):=\sum_{i}\Big{(}{1\over 2}(\bar{\psi}_{\sigma}(\xi_{i})\psi_{\sigma}(\xi_{i}+a\boldsymbol{e}_{0})-\bar{\psi}_{\sigma}(\xi_{i}+a\boldsymbol{e}_{0})\psi_{\sigma}(\xi_{i}))-\frac{1}{4m}\sum_{j=1,2}(\bar{\psi}_{\sigma}(\xi_{i}+a\boldsymbol{e}_{j})-\bar{\psi}_{\sigma}(\xi_{i}))({\psi}_{\sigma}(\xi_{i}+a\boldsymbol{e}_{j})-{\psi}_{\sigma}(\xi_{i}))-\mu\bar{\psi}_{\sigma}(\xi_{i})\psi_{\sigma}(\xi_{i})\Big{)}. Adding \sum_{i}\left(\begin{array}[]{c}\bar{\psi}_{\uparrow}\\ \psi_{\downarrow}\end{array}\right)^{t}(\xi_{i}){\mathbb{\Gamma}}(\theta)\left(\begin{array}[]{c}\psi_{\uparrow}\\ \bar{\psi}_{\downarrow}\end{array}\right)(\xi_{i}) to , and subtracting the same term to , would have led back formally to the initial theory – up to boundary terms on . Then it would have been natural to add to the interaction an infinitesimal bulk term, {\cal V}+\delta{\cal V}_{\theta}\to{\cal V}+\delta{\cal V}_{\theta}+\sum_{i}\left(\begin{array}[]{c}\bar{\psi}_{\uparrow}\\ \psi_{\downarrow}\end{array}\right)^{t}(\xi_{i})\,\left(\begin{array}[]{cc}0&\varepsilon e^{-{\rm i}\theta}\\ \varepsilon e^{{\rm i}\theta}&0\end{array}\right)\,\left(\begin{array}[]{c}\psi_{\uparrow}\\ \bar{\psi}_{\downarrow}\end{array}\right)(\xi_{i}) and let
[TABLE]
where is the normalized measure obtained by modifying as indicated. Our (unproven but unlikely to be difficult) claim is that models (0.39) and (0.40) have the same infra-red behavior.
B. (bosonic regime, or low-energy theory) Diagrams involving a small number () of low-momentum Cooper pairs are a priori divergent, because the Goldstone boson propagator is not integrable (see Theorem 2 below). Ward identities, however, imply that these are actually convergent. On the other hand, it would be a very bad idea to reorganize the series of perturbations in such a way as to produce systematically all Goldstone bosons, since this would lead to a local accumulation of low-momentum bosonic fields which cannot be controlled perturbatively. Thus what is needed here is a complementary, multi-scale Cooper pair expansion procedure producing only a limited number of Goldstone bosons. This procedure makes it possible to give at last the large scale behavior of -point functions of the theory.
Let us now present the main results of our article. Our starting point, namely, the symmetry-broken measure , is equal to the measure defined in (0.39), with a precise choice of cut-offs. First, we have:
Theorem 1 (construction of the model). Assume is small enough. Consider the Grassmann measure of Definition 1.1, re-expressed in terms of modified parameters . Denote by the expectation of a fermionic functional with respect to . Let
[TABLE]
and
[TABLE]
Then, for an adequate choice of parameters such that
[TABLE]
*the renormalization flow is well-defined down to scale , leading to vanishing renormalized counterterms , and a small value of . The coefficient of in the effective covariance kernel, see (0.29), is , with . The parameters , together with the pre-gap energy – defined as solution of an approximate gap equation, see Definition 1.2 – are fixed simultaneously by a fixed-point argument. *
Theorem 1 covers exactly the fermionic theory of A., and spans the entire section 2. Because fermions become massive in the infra-red region, the ”true” effective mass (defined in terms of decay of fermionic two-point functions, see Theorem 3) is not directly accessible (it is given in terms of the distance to the real axis of the first pole of the Fourier-transformed two-point function). Nevertheless, may be thought as effective parameters giving the essentially exact location of the renormalized Fermi circle. Going slightly beyond scale , downto scale – which is required for technical reasons pertaining to the bosonic regime – can only improve the precision. As for the scale counterterm coefficient , thanks to its supplementary prefactor, it is sufficiently close to 0 to enter only as a small, single-scale correction to the value of , which is chosen so as to satisfy the gap equation (see Theorem 2). The gap energy may equivalently be expressed in terms of the renormalized mass, since .
Cooper pairs are defined as linear combinations of electron pairs in Cooper channel, and , which may be represented in matrix form as \left(\begin{array}[]{c}\bar{\psi}_{\uparrow}\\ \psi_{\downarrow}\end{array}\right)^{t}{\mathbb{\Gamma}}\left(\begin{array}[]{c}\psi_{\uparrow}\\ \bar{\psi}_{\downarrow}\end{array}\right) where {\mathbb{\Gamma}}=\left(\begin{array}[]{cc}0&\Gamma^{*}\\ \Gamma&0\end{array}\right) is an arbitrary off-diagonal Hermitian matrix. Cooper pairs in the channel, i.e. with associated matrix {\mathbb{\Gamma}}^{\perp}=\Gamma_{\phi}\left(\begin{array}[]{cc}0&e^{-{\rm i}(\theta+\frac{\pi}{2})}\\ e^{{\rm i}(\theta+\frac{\pi}{2})}&0\end{array}\right), are dominant only in the low-momentum regime (Regime II and III, or B. in our terminology), so relevant observables are averages of these composite fields in a large volume, or correlations at distances .
Theorem 2 (gap equation, Cooper pair correlations at large scale).
- (i)
The gap equation for (see Remark below Theorem 1), defined by fixed point as the value of for which the local part of the fermionic two-point function vanishes, see Definition 3.2, has a unique solution
[TABLE]
- (ii)
Fix as in (i), and consider the associated theory with measure . Let (). If is a two-by-two matrix, we define for (small)
[TABLE]
Then -point functions of Cooper pairs have the following asymptotic large-scale behavior:
[TABLE]
where {\bf{\Gamma}}^{\perp}:=\Gamma_{\phi}\left(\begin{array}[]{c}e^{-{\rm i}(\theta+\frac{\pi}{2})}\\ e^{{\rm i}(\theta+\frac{\pi}{2})}\end{array}\right) is perpendicular to {\bf{\Gamma}}^{//}:=\Gamma_{\phi}\left(\begin{array}[]{c}e^{-{\rm i}\theta}\\ e^{{\rm i}\theta}\end{array}\right); the sum ranges over all pairings \Big{\{}\{i_{1},i_{2}\},\cdots,\{i_{2n-1},i_{2n}\}\Big{\}} of indices, and the two-by-two matrix kernel , whose entry in the -channel is interpreted as propagator of the Goldstone boson, is defined and studied in §3.2 and 3.4. In particular, see (3.13) and (3.112), if ,
[TABLE]
where
[TABLE]
is roughly equal to the (renormalized) Fermi velocity ; and
[TABLE]
The kernel is not integrable, in the sense that the integral diverges.
The gap parameter , as well as the Goldstone boson parameters , are obtained as limits of scale-dependent converging series, , , see Definition 3.4 and §3.4.
Remark 1. The Fourier transform of , as follows by straightforward computations from (0.66), is of the form
[TABLE]
and diverges for . Quite remarkably, the fact that has a pôle precisely at , implying a non-massive behavior, is shown by a Ward identity (see §3.3) to be equivalent to the gap equation. We actually use most of the time the latter characterization in terms of as a definition of the gap equation, a convenient choice for estimates.
Remark 2. Error terms in Theorem 2 (ii) go to 0 as . Details of the proof, see in particular (2.123) in section 2 and (3.113) in section 3, give error terms denoted by ” which are bounded by or equivalently ; a more detailed computation of lowest-order terms would allow to replace by . The leading order of in (i), on the other hand, depends both on the ultra-violet cut-off and on subleading diagrams.
Remark 3. The above Theorem only involves the -component of . Intermediate computations produce a two-by-two matrix propagator , diagonal in the basis . However, the other diagonal coefficient is massive, which translates into a quasi-exponential decay at large distances, similar to that of the fermions. This case is studied alongside the case of isolated fermions in Theorem 3 below.
We can now tentatively recast our results in terms of the Hubbard-Stratonovich transformed theory (0.27) involving the complex boson field . Namely, by Gaussian integration by parts (see Appendix), -point functions of the -field are in one-to-one correspondence with the Cooper pair -point functions of ((ii)), from which we can conclude (leaving aside mathematical rigor) to the following. The Goldstone boson behaves like a non-massive particle described by a non-linear sigma model with coupling constant having the dimension of an energy, namely, the statistical weight for a function with values on the circle is , where \nabla:=\left(\begin{array}[]{c}\partial_{\tau}\\ v_{\phi}\boldsymbol{\nabla}\end{array}\right) . In rescaled, non-dimensional units, , with , the weight for becomes with
[TABLE]
This is the signature of a weakly coupled, super-critical -dimensional non-linear sigma model, known (at least perturbatively) to be a free theory at large scales (see e.g. [57], §13.3), from which the qualitative behavior ((ii)) can be expected. Note, however, that the interaction between Cooper pairs is mediated by massive individual electrons.
Theorem 3 (isolated electrons). -point functions of electrons which are not in the Cooper pair channel decay quasi-exponentially at large distance with a decay rate , namely, for every , there exists such that, if ,
[TABLE]
and similarly,
[TABLE]
and, letting {\mathbb{\Gamma}}^{//}:=\Gamma_{\phi}\left(\begin{array}[]{cc}0&e^{-{\rm i}\theta}\\ e^{{\rm i}\theta}&0\end{array}\right),
[TABLE]
Remark. We fall short of proving exponential decay for two-point functions of electrons. This is mainly due to our choice of cut-offs for the Gaussian covariance kernel, which yields quasi-exponentially decaying scaled kernels. However, non-Gaussian cut-offs of the form , similar to those used in [15], would produce exponentially decaying kernels, and most likely allow us to prove an exponential decay in , resp. for two-point, resp. four-point functions, instead of the quasi-exponential decay factors \Big{(}1+\Gamma_{\phi}\eta^{-1}(|\tau-\tau^{\prime}|+\frac{m^{*}}{p_{F}^{*}}|\boldsymbol{x}-\boldsymbol{x}^{\prime}|)\Big{)}^{-n} in (0.71), resp. (0.72), implying that the effective mass of electrons is .
Let us add a few comments. Cooper pairs are the main observables at low energy. On the other hand, isolated electrons have become massive. The large-scale behavior of isolated fermions and Cooper pairs written down in Theorem 3 points out to a superconductive phase, qualitatively different from the normal phase, which is that of the free electron gas obtained for , or more generally of Fermi liquids. Namely, for (i.e. in the normal phase), specializing to vanishing -coordinates, the one-point density kernel is proved by a standard computation (see (5.32) and (5.33)),
[TABLE]
and the density-density correlation kernel (connected two-point function)
[TABLE]
These correlation functions decrease polynomially at large distances, to be compared with the quasi-exponential decrease exhibited in Theorem 3. Two-point functions of Cooper pairs in arbitrary channels behave similarly with a cubic inverse distance decrease,
[TABLE]
where {\mathbb{\Gamma}}(\theta_{i}):=\Gamma_{\phi}\left(\begin{array}[]{cc}0&e^{-{\rm i}\theta_{i}}\\ e^{{\rm i}\theta_{i}}&0\end{array}\right), , as opposed to the linear inverse distance decrease shown in Theorem 2 for the superconductive phase.
Theorem 4 (phase transition). Let and . Then
[TABLE]
*with . *
This is in agreement with the mean-field prediction for the average of Cooper pairs in a superconducting ground-state oriented in direction , namely, e.g.
[TABLE]
if . Thus an infinitesimally small ”Cooper magnetic field” (i.e. parameter conjugate to Cooper pair field, not to be confused with a true magnetic field pairing to individual electrons) in a given direction is enough to orient Cooper pairs in that direction.
In the following theorem, we consider the integral of a Cooper pair field over a space-time domain ”comparable” to a scale box. Polynomial and exponential moments of , i.e. expectations , may be computed as the integral w.r. to a probability measure denoted by , e.g. . Then upper bounds for quantities , , may be interpreted as large deviation estimates.
**Theorem 5 (local transverse behavior of Goldstone bosons).
**Let be a domain such that \{\xi=(\tau,\boldsymbol{x})\ |\ \Gamma_{\phi}\Big{(}|\tau|+\frac{m^{*}}{p^{*}_{F}}|\boldsymbol{x}|\Big{)}\leq c\}\subset\Omega\subset\{\xi=(\tau,\boldsymbol{x})\ |\ \Gamma_{\phi}\Big{(}|\tau|+\frac{m^{*}}{p^{*}_{F}}|\boldsymbol{x}|\Big{)}\leq C\}, where are constants, and
[TABLE]
where \sigma(\theta):=\cos(\theta)\sigma^{1}+\sin(\theta)\sigma^{2}=\left(\begin{array}[]{cc}0&e^{-{\rm i}\theta}\\ e^{{\rm i}\theta}&0\end{array}\right) is a linear combination of off-diagonal Pauli matrices, see (1.19). Let
[TABLE]
Then, for small enough,
[TABLE]
where
[TABLE]
This is in general agreement with mean-field estimates for the curvature of the effective potential (see (0.30)) near its minimum,
[TABLE]
in a neighborhood of , see [22]. Namely, the field may be considered as roughly constant over . Through the correspondence between the Cooper pair field and its conjugate field , one may assimilate the event \big{\{}|X_{\Omega}-\bar{X}_{\Omega}|>\eta\frac{\Gamma_{\phi}}{\lambda}\,{\mathrm{Vol}}(\Omega)|\big{\}} with the event \big{\{}\big{|}\,|\Gamma(\xi)|-\Gamma_{\phi}\big{|}>\eta\Gamma_{\phi},\ \xi\in\Omega\big{\}}, which has probability . For larger domains , the large deviation rate -\log\Big{(}{\mathbb{P}}\Big{[}\int_{\Omega}d\xi\,\big{|}\,|\Gamma(\xi)|-\Gamma_{\phi}\big{|}>\eta\Gamma_{\phi}{\mathrm{Vol}}(\Omega)\Big{]}\Big{)} is expected from (0.83) to increase linearly in .
0.4 Outline and notations
The paper is organized as follows.
Section 1 is a long, introductory section, where we gradually introduce notations, recall standard facts about the BCS model, reformulate the initial model as a multi-scale model, and evaluate bubble diagrams. Main points are Definition 1.1, where our model is precisely defined, and §1.7, where a two-by-two matrix-valued kernel – whose -component is interpreted as propagator of the Goldstone boson – and a gap equation are introduced in perturbative terms (Definition 1.9). Section §1.7 goes well beyond mean-field theory, but remains perturbative; hence the way in which the gap energy and the -kernel are defined differs in details from the non-perturbative definition of section 3. The precise construction in section 3 is heavily inspired from, but does not rely upon, the one in §1.7, so details in §1.7 are not really important, and the whole subsection – apart from the displacement procedure, which is accurately described there – may be skipped from a mathematical point of view. Yet we hope that §1.7 will help the reader find his way through the much more technical section 3.
Section 2 is dedicated to the high-momentum or fermionic theory (part A. of our scheme). It relies mainly on expansions. Formally, these apply to the dressed model, see Definition 2.6, which interpolates between the coupled model of section 1 and a fully-decoupled model. The horizontal (cluster) and momentum-decoupling (vertical) expansions implement the decoupling between degrees of freedom required to renormalize the diverging contributions, and show convergence of the series obtained by formally expanding . A simple bubble resummation procedure provides the ”-argument” required to control the series in the neighborhood of the transition scale.
Section 3 is concerned with the low-momentum or bosonic theory (part B.) As in QED, Ward identities are proved, which show that – despite the fact that Goldstone bosons is non-massive – all diagrams are infra-red convergent. By a careful complementary multi-scale expansion of Cooper pairs, -point functions are rewritten in terms of a sum over multi-scale trees with vertices connected by -kernels. The end of the section is dedicated to the proof of the Theorems stated in §0.3.
Perspectives are presented in section 4.
Finally, some technical lemmas are collected in section 5.
Important notations. The reader will find many more notations in the body of the article. We hope that this selection will help the reader find his way through the article and prevent any confusion. A more comprehensive index of notations with reference pages is given at the end of the text.
Following the physicists’ convention, the star (∗) denotes complex conjugation. Alternatively, real parameters with a star denote effective parameters; Gaussian measures and covariance kernels with a star are computed with the effective parameters.
The gap parameter is a positive real number.
– Space-time points are typically denoted by , , .
– Fermions are denoted either by \psi=(\psi_{\sigma})_{\sigma=\uparrow,\downarrow}=\left(\begin{array}[]{c}\psi_{\uparrow}\\ \psi_{\downarrow}\end{array}\right), \bar{\psi}=(\bar{\psi}_{\sigma})_{\sigma=\uparrow,\downarrow}=\left(\begin{array}[]{cc}\bar{\psi}_{\uparrow}&\bar{\psi}_{\downarrow}\end{array}\right) or \Psi=(\Psi_{\sigma})_{\sigma=\uparrow,\downarrow}=\left(\begin{array}[]{c}\psi_{\uparrow}\\ \bar{\psi}_{\downarrow}\end{array}\right), \bar{\Psi}=(\bar{\Psi}_{\sigma})_{\sigma=\uparrow,\downarrow}=\left(\begin{array}[]{cc}\bar{\psi}_{\uparrow}&{\psi}_{\downarrow}\end{array}\right). The latter are Nambu spinors, which mix and -components. Thus Cooper pairings are (depending on the context) written as or , or (in matrix notation) \bar{\Psi}\left(\begin{array}[]{cc}0&0\\ 1&0\end{array}\right)\Psi, \bar{\Psi}\left(\begin{array}[]{cc}0&1\\ 0&0\end{array}\right)\Psi, combining into (in terms of Pauli matrices, see (1.19)), or more generally, using an off-diagonal hermitian matrix, , where {\mathbb{\Gamma}}:=\left(\begin{array}[]{cc}0&\Gamma^{*}\\ \Gamma&0\end{array}\right), with . If , resp. , then we note , , resp. , . We also consider the complex-valued vector \boldsymbol{\Gamma}=\left(\begin{array}[]{c}\Gamma^{*}\\ \Gamma\end{array}\right), in particular, and . Seeing either as a complex number, a vector or a matrix will turn out to be useful. In particular, is the scalar product of with , or indifferently with , or with . All quantities are equivariant w.r. to rotations by an angle , acting as follows,
[TABLE]
and
– Part A. is dedicated in a large part to a multi-scale analysis, in which power-counting is directly dependent on the vicinity to the Fermi circle. Thus the norm on momenta,
[TABLE]
with , measures the distance to the (extended to -dimensional description) Fermi circle defined by . Dualizing, we let
[TABLE]
– On the other hand, the analysis in part B. focuses on Cooper pairs with small total momentum , entering into fermionic diagrams of A. as a small transfer momentum. Therefore, the Euclidean norm adapted to this context,
[TABLE]
vanishes when . Dualizing, we let
[TABLE]
Note that and are as different as can be. In order to further avoid confusions, we systematically denote electron momenta by the letter and Cooper pair momenta by the letter . On the other hand, may be used interchangeably most of the time.
– A fermion momentum has scale if . A Cooper pair momentum has Euclidean scale if . Momentum cut-offs , resp. select momenta of scale , resp. Euclidean scale , namely , resp. vanish except if , resp. . The sum , resp. are infra-red, resp. ultra-violet cut-offs; similarly for , resp. . Infra-red, resp. ultra-violet cut-off fermion fields are similarly denoted by , resp. .
– The amplitude of an amputated Feynman diagram is denoted by . Letting , resp. be the location of the internal, resp. external vertices of , is equal to the integral \big{(}\prod_{i=2}^{N_{ext}}\int d\xi_{ext,i})\big{(}\prod_{i=1}^{N_{int}}\int d\xi_{int,i}) of the product of the internal vertices by the internal propagators. Because of translation invariance, the location of one of the external vertices has been kept fixed. Restricting to propagators of momentum scales , one obtains instead a scale infra-red cut-off amplitude denoted by .
– A Cooper pair with momentum (or equivalently, a pair of electrons in Cooper channel) is an integrated composite field or with total momentum such that . Such fields are the main subject of part B.. However, they also contribute in an unessential way to the analysis of part A., since they can enter even high-energy diagrams as external legs of fermion loops, though with a very small relative volume. By abuse of language, we shall sometimes also speak of ”Cooper pairs” even when the total momentum is such that ; such pairs are simply unphysical quantities entering as virtual particles into diagrams.
– Sign conventions regarding Grassmann integrals are recalled in Appendix, §5.1.
– Finally, is a constant. Bounds for multi-scale polymers with , resp. external legs proceed differently due to the complicated angular dependence of external legs.
Acknowledgements. One of us (J. Magnen) would like to express his warmest and most sincere thanks to J. Feldman, V. Rivasseau and E. Trubowitz for extended discussions on the topics of this article and related subjects.
1 The gap equation
This long section starts with a series of reformulations of the initial model defined by the Grassmann measure
[TABLE]
recall that is Grassmann Gaussian measure cut-off at some scale , with covariance kernel given (up to cut-off issues) by the operator \Big{(}{\rm i}p^{0}-e(\boldsymbol{p})\Big{)}^{-1}, see (0.9), and
[TABLE]
see (0.15), again up to cut-off issues. The issue is to motivate and define precisely the symmetry-broken Grassmann measure (see Theorems in the Introduction) defining the model.
Let us describe informally the various stages:
• first (see §1.1), the -number (charge) symmetry is broken at energies comparable to the energy gap (i.e. in Regime II), yielding a non-zero value in the ground state for products of spin-neutral annihilation operators, for instance
[TABLE]
This is a direct consequence of the fact that Cooper pairs, i.e. neutral pairs of electrons form bound states. On the other hand, non spin-neutral bound states do not form. This motivates the introduction of Nambu fields . The well-known Bogolioubov-De Gennes-Gor’kov mean-field, quadratic Hamiltonian [7, 8] (which we do not discuss here) is diagonalized in a rotated basis expressed in terms of , and the mean-field ground-state is similarly obtained from the normal state ground-state by applying products of ’s, see e.g. [2], §6.4.
• second (see §1.2), the Cooper pair contribution to effective action is written in the Nambu basis in terms of an off-diagonal matrix \int d\xi\,\bar{\Psi}(\xi)\left(\begin{array}[]{cc}0&\Gamma^{*}\\ \Gamma&0\end{array}\right)\Psi(\xi), where the energy gap (coinciding with the above ground state expectation value) satisfies a consistency equation called gap equation (discussed later on). Pauli matrices help encode the separation of into its real/imaginary parts. This leads to a modified Grassmann Gaussian measure with covariance kernel in the form of a rational function, whose denominator is bounded near the Fermi circle. Its inverse Fourier transform therefore decreases exponentially at large distances like , with . Thus the energy gap also plays the rôle of an inverse correlation length for the fermions.
• further (see §1.3), the energy gap (as follows from the analysis in [24]) may be defined (neglecting small corrections due to scales ) as minus the inverse of the off-diagonal interacting Green function evaluated at such that and is located on the Fermi sphere. However, the radius of the Fermi sphere – defined as the singular locus of the theory – must be renormalized. Since renormalization is implemented à la Wilson as a flow of the parameters of the theory, we actually rewrite the model in terms of its effective (i.e. far infra-red) parameters and . This implies introducing quadratic counterterms proportional to at bare scale . If have been chosen correctly, then counterterms , for parameters vanish. Now, the flow is computed using a multi-scale decomposition of the fields,
[TABLE]
Scale counter-terms are then determined by computing two-point functions of the scale infra-red cut-off theory with external momenta located on the renormalized Fermi circle. The above arguments explain why the precise form of the UV cut-off (0.4) (in itself largely irrelevant) is not fixed from the beginning, but rather as part of (1.4).
The model used as basis for discussion in the present article is that introduced in §1.3, though (for pedagogical reasons) some of its features are specified only in §2.1. As for the rest of the section, it is dedicated to the study of the gap equation.
Starting with perturbation theory, we discuss in §1.4 and §1.5 ladder diagrams (see [26], §11), which we call here for obvious diagrammatic reasons bubble diagrams. The evaluation of these diagrams with transfer momentum close to zero yields the main contribution to the energy gap . As a concession to the mean-field approximation scheme, we resum ladder diagrams (see §1.6) in an intermediate step, and deduce the approximate value for known in textbooks. In our rigorous approach, the sum of ladder diagrams plays a quintessential rôle – as we shall see, ladder diagrams in their Cooper pairing channel are alone responsible for the superconducting transition, other diagrams yielding only finite corrections – but does not give the exact value of .
The last subsection (§1.7) – which lies half-way betweeen mean-field theory and rigorous arguments – shows how to re-sum perturbatively all diagrams forming the Goldstone boson propagator (the task is done beyond perturbation theory in §3.2). The series is formed by alternating bubble diagrams with two-point irreducible diagrams forming the Bethe-Salpeter kernel. Definition 1.9 defines as the value of the parameter such that diverges when , a requirement that – as mentioned in a Remark in the introduction – is equivalent to the exact gap equation. This means that Goldstone bosons, as expected, are not massive. The diagram series naturally forms a two-by-two, diagonal matrix kernel \Sigma=\left(\begin{array}[]{cc}\Sigma_{//,//}&0\\ 0&\Sigma_{\perp,\perp}\end{array}\right), whose component is shown to be massive, contrary to .
1.1 Green functions and Nambu formalism
The main quantity controlling the electron pairing is the Green function,
[TABLE]
In terms of the so-called Nambu fields,
[TABLE]
the above matrix is simply the tensor quantity . The particle number symmetry
[TABLE]
is broken by Cooper pairing. However, the remaining symmetries, namely, * spin symmetry*
[TABLE]
invariance under * involution*
[TABLE]
and spatial rotational invariance imply (see [24], Lemma II.1) that G=\left(\begin{array}[]{cc}G_{1,1}&G_{1,2}\\ G_{2,1}&G_{2,2}\end{array}\right) is rotationally invariant, and (in Fourier coordinates)
[TABLE]
[TABLE]
In particular, * is a Hermitian matrix of the form \left(\begin{array}[]{cc}a&b\\ b^{*}&-a\end{array}\right) with .*
In this language, the quadratic part of the action, see (0.9), is (up to the ultra-violet cut-off)
[TABLE]
So the Grassmann Gaussian measure obtained by letting is characterized (again, up to UV cut-off) by
[TABLE]
Then the interaction , see (1.2), rewrites as , with
[TABLE]
where or indifferently (since the square depends only on the determinant of the matrix inserted between and ). Choosing or more generally emphasizes the rôle of Cooper pairings in a given direction .
1.2 Symmetry-broken Grassman Gaussian measure
The BCS phase transition and the above symmetry considerations suggest that we:
- (i)
substitute to the bare dispersion relation the effective dispersion relation,
[TABLE]
in terms of the effective parameters where:
[TABLE]
is the Fermi sphere radius renormalization; and
[TABLE]
is the mass renormalization defined in such a way that ;
- (ii)
and rewrite the functional integral in terms of the unknown renormalized parameters and unknown gap as
[TABLE]
where:
•
[TABLE]
are the usual Pauli matrices, and \left(\begin{array}[]{c}\Gamma_{1}(\theta)\\ \Gamma_{2}(\theta)\end{array}\right)=\left(\begin{array}[]{cc}\cos\theta&\ -\sin\theta\\ \sin\theta&\cos\theta\end{array}\right)\left(\begin{array}[]{c}\Gamma_{\phi}\\ 0\end{array}\right)=\Gamma_{\phi}\left(\begin{array}[]{c}\cos\theta\\ \sin\theta\end{array}\right) is a vector in of norm pointing in the direction ;
• * is the free Grassmann measure characterized by its covariance kernel,*
[TABLE]
and
[TABLE]
In (1.20), {\mathbb{\Gamma}}(\theta):=\Gamma_{1}(\theta)\sigma^{1}+\Gamma_{2}(\theta)\sigma^{2}=\left(\begin{array}[]{cc}0&\bar{\Gamma}(\theta)\\ \Gamma(\theta)&0\end{array}\right) is an off-diagonal Hermitian matrix, in conformity with (1.11).
Let us emphasize that the measure \exp\Big{(}+\int dp\,\bar{\Psi}(-p)\Big{(}\Gamma_{1}(\theta)\sigma^{1}+\Gamma_{2}(\theta)\sigma^{2}+\frac{\delta m}{m^{2}}\frac{|\boldsymbol{p}|^{2}}{2}+\delta\mu\Big{)}\Psi(p)\Big{)}d\mu^{*}_{\theta}(\Psi,\bar{\Psi}) is by construction independent of and , and equal to the original measure of (0.8), up to error terms due to the ultra-violet cut-off.
1.3 Multi-scale analysis
Let be the radius of the effective Fermi circle
[TABLE]
defined by for . Near the Fermi circle, one has
[TABLE]
where is a signed measure of the distance of to the Fermi circle. Let be some momentum such that , see Assumption (0.5). Given the form of the denominator of , see (1.20), it is natural to impose the UV cut-off conditions ; namely, it can easily be seen that the region plays practically no rôle in the theory. Assumption (0.4) implies
[TABLE]
A generic point in will in general be denoted by , or , , etc. By convention is the Minkowski scalar product, while
[TABLE]
is a homogeneous norm chosen in such a way that if is close to the Fermi circle, see (1.23).
Remark. Note that has nothing to do with the Euclidean norm used in the low-momentum regime (see section 3).
Eq. (1.24) yields some upper momentum scale for and ,
[TABLE]
Finally, we note that for . Thus
[TABLE]
fixes an energy scale associated to the phase transition. On the other hand, if , then is essentially constant; thus, when (i.e. for the free theory), all momenta such that may be treated as being of the same scale. It turns out not to be the case for the interacting theory. In order to make the bridge between the high- and the low-momentum regimes, we fix a lower energy scale ,
[TABLE]
In this section, we only discuss momenta such that
[TABLE]
We now define a multi-scale version of our model according to the above principles. Fix some smooth function supported on such that on , where . Assume is further decomposed as a finite sum of smooth functions such that ; the will be specified only in §2.1 (with indexing momentum angular sectors), since it does not matter at this stage. Let
[TABLE]
so that is a scale infra-red cut-off (coupled with a scale UV cut-off), resp. is a scale UV cut-off. Define as in [25] the matrix-valued kernel
[TABLE]
[TABLE]
and similarly (assuming )
[TABLE]
[TABLE]
and let be independent free spin fermion fields, either on or on a finite volume , with covariance kernel
[TABLE]
If is a finite volume, then is simply extended outside of by [math]. Let
[TABLE]
Then
[TABLE]
has covariance kernel
[TABLE]
Its Fourier transform is equal to for .
Summarizing: for , is given by an integral over momenta such that |p|\approx\max\Big{(}|p^{0}|,|e^{*}(\boldsymbol{p})|\Big{)}\approx 2^{-j}\mu, or equivalently: .
Definition 1.1** (Grassmann measure)**
- (i)
Let
[TABLE]
where is the Grassmann Gaussian measure associated to the fields as in (1.31), (1.37), and
[TABLE]
[TABLE]
- (ii)
(same with scale infra-red cut-off) For , we let
[TABLE]
where is the Grassmann Gaussian measure associated to the fields as in (1.31), (1.37), and
[TABLE]
[TABLE]
The whole point in introducing the ’s at this very early stage is that the quadratic part in the above Definition is diagonal in that basis.
Up to change of normalization constant , is simply obtained from by letting for . In practice, this -dependent measure will result from the momentum-decoupling expansion and renormalization by replacing with , and letting . It is the right place to mention that these infra-red cut-off theories enjoy the same invariance properties as the original theory, see §1.1, which is crucial for the renormalization step (see §2.3).
1.4 Bubble diagrams and fermion four-point function in Coooper channel
We primarily consider in this subsection the simplest Feynman diagrams of the theory – the bubble diagrams. Feynman rules can be found in Chapter 46 of the book by Fetter and Walecka [26]. The amplitude of an amputated bubble diagram , computed using the free Grassmann measure , is denoted in whole generality by , where is the external momentum flowing into the diagram, or if one wishes to emphasize the dependence on the parameter . Replacing with , one obtains the scale infra-red cut-off evaluation .
The symmetry-broken theory with gap parameter involves
• diagonal or normal propagators,
\bar{\psi}_{\sigma_{1}}$$\psi_{\sigma_{2}}
evaluated as . Graphically, we orient such propagators using arrows from to , following the convention of e.g. quantum electrodynamics, in order to follow fermion loops;
• off-diagonal or symmetry-broken propagators,
\psi_{\sigma_{1}}$$\psi_{\sigma_{2}},\bar{\psi}_{\sigma_{1}}$$\bar{\psi}_{\sigma_{2}}
evaluated as , resp. . Such propagators are non-oriented. Because we introduced them by hand in the Gaussian measure, they exist even in Regime I, i.e. above the gap energy, but they are very small for scales due to the coefficient in the numerator.
The names diagonal and off-diagonal refer to the location of these quantities inside the two-by-two free Green function G\Big{|}_{\lambda=0}, see (1.5).
• -counterterms
-$$\Gamma
coming from the term in ((ii)).
In subsequent computations, we rewrite diagonal propagators as follows,
[TABLE]
with
[TABLE]
The integration measure for angle-independent quantities is of the two-dimensional form . Thus the theory is effectively -dimensional in the neighborhood of the Fermi sphere; this remarkable fact actually holds true for any dimension .
Vacuum polarization bubbles. In an intermediate boson picture, where density fields , interact via a delta-potential represented by a wavy line, such bubbles may be represented as
yielding a renormalization of the intermediate boson.
Bubble diagrams are evaluated without taking into account the usual -factor per loop, characteristic of fermionic theories. One obtains three possibilities:
• The normal -bubble diagram
\bar{\psi}_{\downarrow}$$\psi_{\downarrow}$$\psi_{\downarrow}$$\bar{\psi}_{\downarrow}$$\bar{\psi}_{\uparrow}$$\psi_{\uparrow}$$\psi_{\uparrow}$$\bar{\psi}_{\uparrow}$$q$$q
Fig. 1.4.1. Normal bubble diagram .
is evaluated up to error terms (see (1.45) and below) as
[TABLE]
A naive power-counting argument in terms of the integration measure, also valid for the Cooper-pair bubble below, yields
[TABLE]
a logarithmically divergent integral. However, is UV convergent, as can be seen from a contour deformation into the upper/lower half-plane of the integral in . Namely, \int\frac{dp^{0}}{({\rm i}p^{0}+e^{*}_{|\Gamma|}(\boldsymbol{p}))^{2}}={\rm i}\int dp^{0}\,\frac{d}{dp^{0}}\Big{(}\frac{1}{{\rm i}p^{0}+e^{*}_{|\Gamma|}(\boldsymbol{p}))}\Big{)}=0, so (by integration by parts in and parity in )
[TABLE]
Now, vanishes except for momenta of scales or , which thus suppresses the diverging sum over scales. Generalizing, we let be a small transfer momentum, , and subtract , yielding
[TABLE]
with . Restricting the integration domain to momenta of scale and moving the contour in to a distance from the real axis yields , with , summing to .
The same bound may also be obtained by noting that , , and
[TABLE]
which vanishes by translation invariance upon integration in up to cut-off issues; we skip details since we do not need a second argument, but refer the reader to Ward identities in quantum electrodynamics which can be proved in a similar way, see e.g. [57].
Error terms involve the subtacted terms, or or both. However, the small factor appearing in (1.46) transforms the logarithmically divergent estimates (1.48) into a convergent integral; the same argument holds for below.
• The off diagonal or symmetry-broken -bubble diagrams
\psi_{\downarrow}$$\psi_{\uparrow}$$\psi_{\uparrow}$$\psi_{\downarrow}$$\bar{\psi}_{\downarrow}$$\bar{\psi}_{\downarrow}$$\bar{\psi}_{\uparrow}$$\bar{\psi}_{\uparrow}$$q$$\bar{\psi}_{\uparrow}$$\bar{\psi}_{\downarrow}$$\psi_{\downarrow}$$\psi_{\uparrow}$$\bar{\psi}_{\downarrow}$$\bar{\psi}_{\uparrow}$${\psi}_{\uparrow}$${\psi}_{\downarrow}$$q
Fig. 1.4.2. Symmetry broken bubble diagrams , resp. (from left to right).
are evaluated at zero external momentum as
[TABLE]
which are equal by simple scaling to a constant times , resp. times . Letting and summing over scales (in this computation one may simply set ), one gets
[TABLE]
a finite contribution.
The diagram has a Cooper pair external structure; it connects a Cooper pair to another one. Later on (see §1.7), following the Nambu convention, see (1.5), we shall denote with a lower index ”diag” diagrams connecting a Cooper pair to the adjoint pair or conversely, and ”off” diagrams connecting two Cooper pairs of the same kind. This accounts for the notation . On the other hand, is a companion diagram to , connecting two conventional pairs . Rewriting the product of two vertices
[TABLE]
in the form (\bar{\psi}_{\uparrow}\bar{\psi}_{\downarrow})(\xi)\ \big{\{}(\psi_{\uparrow}(\xi)\psi_{\downarrow}(\xi^{\prime}))\,\cdot\,(\psi_{\downarrow}(\xi)\psi_{\uparrow}(\xi^{\prime}))\big{\}}\ (\bar{\psi}_{\uparrow}\bar{\psi}_{\downarrow})(\xi^{\prime}), one sees that is a contribution to the off-diagonal channel , equal to times a quantity which is negative, equal to , when .
• The main diagram in BCS theory is the Cooper pair bubble ,
\psi_{\downarrow}$$\bar{\psi}_{\downarrow}$$\psi_{\uparrow}$$\bar{\psi}_{\uparrow}$$\bar{\psi}_{\downarrow}$$\psi_{\downarrow}$$\bar{\psi}_{\uparrow}$$\psi_{\uparrow}$$-p$$p-q$$q
Fig. 1.4.3. Cooper pair bubble.
Remark. Ladder diagrams. In the intermediate boson picture as above, such bubbles take the form
where wiggling lines
\xi$$\xi^{\prime}
represent the -interaction ,
and may be concatenated in the form of ladders of arbitrary length,
Fig. 1.4.4. Ladder diagram.
conventionally called ladder diagrams in connection to the Bethe-Salpeter kernel, see [26], or equivalently in the following dual form,
Fig. 1.4.5. Bubble chain.
The above Cooper pair bubble diagram is evaluated up to convergent error terms as
[TABLE]
It is UV divergent: at zero external momenta,
[TABLE]
Rewriting (1.54) in the form (\bar{\psi}_{\uparrow}\bar{\psi}_{\downarrow})(\xi)\ \big{\{}(\psi_{\uparrow}(\xi)\bar{\psi}_{\uparrow}(\xi^{\prime}))\,\cdot\,(\psi_{\downarrow}(\xi)\bar{\psi}_{\downarrow}(\xi^{\prime}))\big{\}}\ (\psi_{\downarrow}\psi_{\uparrow})(\xi^{\prime}), one sees that is a contribution to the diagonal channel , which is positive when .
This quantity diverges logarithmically in the neighbourhood of the Fermi circle when . The estimates (1.56) also holds for when and . On the other hand, if , then the above integral has an infra-red cut-off of scale , so that
[TABLE]
is essentially independent of .
Both the transfer momentum and the infra-red cut-off at fermionic scale play the rôle of an infra-red cut-off; the resulting effective IR cut-off fermionic scale is k\simeq\min\Big{(}j,\lfloor\log(\mu/|q|_{+})\rfloor\Big{)}. In particular,
[TABLE]
if the transfer momentum is small, i.e. .
On the other hand, as proved in Appendix, see §5.4, if ,
[TABLE]
so is roughly equal to the evaluation of a bubble with zero transfer momentum and infra-red cut-off ; the reason is that (say, for ), the integrand in (a generalization of (1.56)) vanishes when and are not on the same side of the Fermi sphere.
It is a remarkable fact per se that, due to symmetry-breaking, the above integral (1.56) has an effective cut-off, making the diagram finite.
1.5 The gap equation at lowest order
We consider here only the Cooper pair bubble diagram and interpret the gap equation as a self-consistent identity ensuring the vanishing of the one-point function in the direction .
We first remark that, as an immediate consequence of (1.20),
[TABLE]
with error terms coming from , see (1.45), graphically,
•\times$$\simeq\ \int d\xi^{\prime}\,$$\xi$$\lambda$$\xi$$\lambda$$\xi^{\prime}••
Fig. 1.5.1. Off-diagonal covariance is equal to Cooper pair bubble diagram.
The diagram on the left-hand side may be interpreted as one of the two main contributions to the one-point function of the Goldstone boson, the other coming from the counterterm insertion at , -$$\Gamma
. Then these two main contributions cancel provided
[TABLE]
or equivalently, taking (1.57) into account,
[TABLE]
This is the crudest approximation to the gap equation. Letting implies that the bubble series,
\bar{\psi}_{\downarrow}$$\psi_{\downarrow}$$\bar{\psi}_{\uparrow}$${\psi}_{\uparrow}$$+\lambda$$\bar{\psi}_{\downarrow}$$\psi_{\downarrow}$$\bar{\psi}_{\uparrow}$${\psi}_{\uparrow}$$+\lambda^{2}$$\bar{\psi}_{\downarrow}$$\psi_{\downarrow}$$\bar{\psi}_{\uparrow}$${\psi}_{\uparrow}$$+\cdots
Fig. 1.5.2. Cooper pair bubble geometric series.
which is the leading order approximation of the Goldstone boson propagator, diverges at zero transfer momentum. These two properties – vanishing of the one-point function in the direction , singularity of the Goldstone boson propagator – remain simultaneously true beyond mean-field theory and even non-perturbatively, as we shall see later on, see §1.7 and section 3. Whatever the exact value of , and the cut-off scale , one has as soon as
[TABLE]
The infra-red behavior of chains of bubbles, and later on, of Goldstone boson propagators, are best understood after a change of basis, choosing for external structures the linear combinations instead of the diagonal channels and the off-diagonal channels , where . (Explicit change-of-basis formulae are provided in §3.1 in the case , see in particular (3.10) and below). Letting , , one obtains in the -basis a diagonal matrix =diag, with
[TABLE]
and , .
It is useful at this point to introduce the following
Definition 1.2** (pre-gap equation)**
Let be the solution of the pre-gap equation,
[TABLE]
Because , may also be understood as the solution with largest module of the eigenvalue equation,
[TABLE]
The matrix has eigenvectors \left(\begin{array}[]{c}\Gamma^{//}\\ (\Gamma^{//})^{*}\end{array}\right),\left(\begin{array}[]{c}\Gamma^{\perp}\\ (\Gamma^{\perp})^{*}\end{array}\right), with respective eigenvalues and [math]. As emphasized in §1.6, this means that the geometric series of Cooper pair bubble diagrams is massive in the parallel direction (), i.e. along , whereas it is non-massive in the perpendicular direction (), i.e. along , corresponding to the longitudinal direction along the circle Cst, in line with Goldstone’s insight. Compare with Definition 1.9. As shown in §1.7 for a corrected version of the gap equation including the Bethe-Salpeter kernel , this equation has a solution , as expected. Computations in 1.7 reduce to the results of §1.6 if one sets . The pre-gap equation is corrected in the bosonic regime, see Definition 3.3, leading to a sequence , converging when to the correct value of .
1.6 The pre-Goldstone boson propagator
The pre-Goldstone boson propagator is directly related to the second derivatives of and w.r. to around . A general argument, relying on the fact that
[TABLE]
and
[TABLE]
implies that \nabla^{2}{\cal A}_{q}^{j\to}(\Gamma,\Upsilon_{3,diag})\Big{|}_{q=0}\leq 0, resp. \nabla^{2}(-{\cal A}_{q}^{j\to}(\Gamma,\Upsilon_{3,off}))\Big{|}_{q=0}\leq 0 . By symmetry, \partial_{q^{0}}{\cal A}_{q}(\Gamma,\Upsilon_{3,\varepsilon})\Big{|}_{q=0}=0, \boldsymbol{\nabla}{\cal A}_{q}(\Gamma,\Upsilon_{3,\varepsilon})\Big{|}_{q=0}=0, \partial_{q^{0}}\boldsymbol{\nabla}{\cal A}_{q}(\Gamma,\Upsilon_{3,\varepsilon})\Big{|}_{q=0}=0 and \boldsymbol{\nabla}^{2}{\cal A}_{q}(\Gamma,\Upsilon_{3,\varepsilon})\Big{|}_{q=0} is a scalar matrix, where =diag,off.
Lemma 1.3
[TABLE]
[TABLE]
where
[TABLE]
have the dimension of an energy, and
[TABLE]
are velocities.
The proof is given in Appendix.
Definition 1.4** (Goldstone pre-boson coupling constant and velocity)**
Let
[TABLE]
The Goldstone pre-boson coupling constant and velocity are obtained by inverting F(q):=\Big{(}{\cal A}_{q}(\Gamma^{(j^{\prime}_{\phi}-1)\to},\Upsilon_{3,diag})-{\cal A}_{0}(\Gamma^{(j^{\prime}_{\phi}-1)\to},\Upsilon_{3,diag})\Big{)}-\Big{(}{\cal A}_{q}(\Gamma^{(j^{\prime}_{\phi}-1)\to},\Upsilon_{3,off})-{\cal A}_{0}(\Gamma^{(j^{\prime}_{\phi}-1)\to},\Upsilon_{3,diag})\Big{)} for , namely,
[TABLE]
where
[TABLE]
similarly to (0.87).
The meaning of Definition 1.4 is explained in the next subsection.
Summing the geometric series of Cooper pair bubble diagrams, one gets a matrix kernel, called pre-Goldstone boson propagator or simply pre-kernel,
[TABLE]
with (choosing a basis)
[TABLE]
This kernel is in sandwich between adjoint fermion pairs in Cooper pairings, or . Note that the amplitude is computed not for the exact value of , but for in such a way that diverges when .
As seen in the previous subsection, the above kernel is divergent in the infra-red limit, i.e. when , in the -channel. Let us now consider but small. The behavior of Pre is obtained (see (1.64)) as the inverse of the quadratic form -{1\over 2}\sum_{i,j=0}^{3}\partial_{q_{i}}\partial_{q_{j}}\Big{\{}({\cal A}_{q}(\Upsilon_{3,diag})+{\cal A}_{q}(\Upsilon_{3,off})\Big{\}}\Big{|}_{q=0}q^{i}q^{j}. Lemma 1.3 yields:
[TABLE]
with the notations of Definition 1.4. On the other hand,
[TABLE]
As we shall see in the next subsection, the kernel component is a good approximation of the (scalar) Goldstone boson propagator for large transfer momenta . For this reason, and also because it is explicitly produced through the bubble resummations in the fermionic regime, we may consider it as the covariance of fictitious bosonic particles which we call pre-Goldstone bosons, whence the name ; in Feynman diagrams it is represented as a dashed wavy line,
. This kernel is dressed by more complicated diagrams, which play an important rôle only in the bosonic regime, but do not modify the general picture described in the previous paragraphs; perturbatively (as shown in (1.141)),
[TABLE]
where , , see (1.135), and , see (0.87). The rigorous definition of in section 3 – which relies instead on a somewhat technical multi-scale construction – will be slightly different, but (1.80) will be shown to hold for some modified values of and up to error terms, for which we only prove a crude bound in (though one should be able to prove that they are by looking more precisely at lowest-order terms with a small number of vertices).
Let us now describe the (non-rigorous) connection to the effective non-linear sigma model theory described in the Introduction (see Regime II). Let . Bubbles in the channel may be resummed into the following counterterm,
[TABLE]
implying an effective action, see (0.27), of the form
[TABLE]
where by definition
[TABLE]
Removing the upper [math] indices, we get a similar resummation of the Goldstone boson. So plays the rôle of a coupling constant for the non-linear sigma model describing the infra-red behavior of the Goldstone boson.
Let us investigate corrections to the leading-order behavior of when (as we shall see in §3, the same estimates hold for , up to the replacement of by their counterparts ). Corrections involve fourth-order derivatives of ,
[TABLE]
Let
[TABLE]
(see Notations in the introduction), . An inverse Fourier transform yields
[TABLE]
Error terms can be bounded, yielding
[TABLE]
as proved in the Appendix.
Finally, let us bound the infra-red cut-off matrix kernel kernel with . For that, we need further estimates on the Fourier transformed kernel . Consider first the case when . We deduce from the computations of §1.4 that . Then, as proved in the Appendix, if ,
[TABLE]
where . The other case, namely when , is easier. As already observed in §1.4 , and similarly,
[TABLE]
so that , and finally, by an elementary computation,
[TABLE]
Introduce a reduced partition of unity for transfer momenta, fix , and let , . Then
[TABLE]
for every , as shown by successive integrations by parts. Whence, in particular,
[TABLE]
When , a good approximation for is simply its leading order, which is the -vertex . The bound (1.92) is coherent with the value of the integrated kernel , since .
On the other hand, the kernel , which is the Fourier transform of the massive kernel
[TABLE]
has a quasi-exponential decay at distances ,
[TABLE]
see (2.17), as can be shown in a standard way by repeated integrations by parts.
1.7 Energy gap and Goldstone boson
The above Cooper pair bubble diagrams , diag,off, may be composed an arbitrary number of times, yielding times the geometric series of bubbles, see Fig. 1.5.1 with transfer momentum , evaluated as . A quick but not rigorous computation, taking into account only these diagrams, yields an approximate value for , based on the requirement that the denominator of the series of bubbles vanishes for , see Gap equation in Definition 1.9 below, namely (considering only the -diagram, and summarizing briefly the findings of the two previous subsections)
[TABLE]
from which
[TABLE]
and
[TABLE]
As explained in the Introduction to this section, what we really need in section 3 is the effective contribution to the Cooper pair bound state (Goldstone boson) propagator of fermionic four-point functions. We need some preliminary explanations before we give a closed formula for this quantity, but let us mention already at this stage that the approximation (1.96) gives a correct order of magnitude for – as shall be checked self-consistently –, and that the value obtained for the bubble series (1.97) for small is a good approximation for – except that the latter is a two-by-two matrix.
First, a Feynman diagram with Cooper pair external structure
[TABLE]
(non-mixing case), or
[TABLE]
(mixing case), is called two-particle irreducible if the compound made up of vertices is connected to the compound made up of vertices . Note that two-particle reducible diagrams may be generated from the partition function by a double Legendre transform, see e.g. [59], eq. (10.132). We could, but have chosen not to, use this nice trick. As detailed in §2.1, external legs possibly have indices – scale and momentum angular indices for sectors with , scale , momentum angular indices and an extra integer index for micro-sectors – , numbered as the corresponding vertex locations, e.g. (\bar{\psi}^{j_{1},\alpha_{1}}_{\uparrow}\bar{\psi}^{j_{2},\alpha_{2}}_{\uparrow})\otimes(\psi^{j_{3},\alpha_{3}}_{\downarrow}\psi^{j_{4},\alpha_{4}}_{\uparrow})=\Big{(}(\bar{\psi}^{j_{1},\alpha_{1}}_{\uparrow}(\xi_{1})\bar{\psi}^{j_{2},\alpha_{2}}_{\uparrow}(\xi_{2})\Big{)}\Big{(}(\psi^{j_{3},\alpha_{3}}_{\downarrow}(\xi_{3})\psi^{j_{4},\alpha_{4}}_{\uparrow}(\xi_{4})\Big{)}. As explained later on in §2.3, only -kernels with external scales need be considered; their external momenta are very close to, but not on, the Fermi circle. Triples index regions of the momentum space of dimensions scaling like , so that, denoting by the center of , for lying inside . None of this is really required to read the present subsection, which remains at a descriptive level based on finite number of graphs produced by perturbation.
We let be the sum of all two-particle irreducible, four-point Feynman diagrams with fixed external structure; it may be represented as a two-by-two matrix, with non-mixing coefficients on the diagonal, and mixing coefficients off-diagonal,
[TABLE]
compare with (1.77). Lowest-order terms are
\lambda$$\bullet$$\alpha_{1}$$\alpha_{3}$$\alpha_{2}$$\alpha_{4}$$+\ \lambda^{2}$$\bullet$$\xi_{1}=\xi_{3}$$\bullet$$\xi_{2}=\xi_{4}$$\alpha_{1}$$\alpha_{3}$$\alpha_{2}$$\alpha_{4}$$+$$\bullet$$\xi_{1}=\xi_{3}$$\bullet$$\xi_{2}=\xi_{4}$$\times$$\times$$\alpha_{1}$$\alpha_{3}$$\alpha_{2}$$\alpha_{4}$$+\ \lambda^{3}$$\xi_{1}=\xi_{3}$$\xi_{2}$$\xi_{4}$$\bullet$$\bullet$$\bullet$$+$$\xi_{1}=\xi_{3}$$\xi_{2}$$\xi_{4}$$\bullet$$\bullet$$\bullet$$\times$$\times
+\ \lambda^{4}$$\bullet$$\bullet$$\bullet$$\bullet$$\xi_{1}$$\xi_{3}$$\xi_{2}$$\xi_{4}$$\alpha_{1}$$\alpha_{2}$$\alpha_{3}$$\alpha_{4}$$+\,\cdots
Fig. 1.7.1. Lowest-order terms of the Bethe-Salpeter kernel.
Note that the triangle diagram, Fig. 0.1, is a two-particle reducible diagram
obtained
by composingand
On the other hand, the diagrams on the first line of Fig. 1.7.1 (obtained by tilting the triangle diagram) are two-particle irreducible – but note that by cutting the two long slanted lines, vertices are cut off from the other two, .
The first diagram is the vertex, simply evaluated as the coupling constant . The second and third diagrams look superficially like bubble diagrams tilted to the side, but the impression is wrong because the momentum circulating inside the ”tilded bubble” diagrams , named with a tilde by reference to the original bubble diagrams of the previous subsection,
p_{1}$$p_{3}$$\psi$$\bar{\psi}$$\psi$$\bar{\psi}$$\bullet$$\bullet$$q\!-\!p_{1}$$-\!q-\!p_{3}$$p_{1}\!-\!p_{3}+r$$r$$p_{1}$$p_{3}$$\psi$$\bar{\psi}$$\psi$$\bar{\psi}$$\bullet$$\bullet$$q\!-\!p_{1}$$-\!q-\!p_{3}$$\times$$\times$$p_{1}\!-\!p_{3}+r$$r$$p_{1}$$p_{3}$$\bar{\psi}$$\bar{\psi}$$\bar{\psi}$$\bar{\psi}$$\bullet$$\bullet$$q\!-\!p_{1}$$-\!q-\!p_{3}$$p_{1}\!-\!p_{3}+r$$r$$\times$$\times
Fig. 1.7.2. ”Tilted bubble” diagrams (from left to right).
isn’t as in the -diagram, see Fig. 1.4.1, but , therefore, generically a large momentum. Note incidentally that the third diagram of Fig. 1.7.1 decomposes as the sum of two diagrams, and .
The first three diagrams, namely, the vertex, and , are of non-mixing type. The fourth one, , is mixing, proportional (when considered together with the conjugate diagram as entries of a matrix) to the off-diagonal matrix \left(\begin{array}[]{cc}0&\Gamma^{2}\\ (\Gamma^{*})^{2}&0\end{array}\right).
Finally, the last ones (the terms in and in Fig. 1.7.1), can be either of non-mixing or mixing type, see the figure below for the mixing version of the -diagram.
\times$$\Gamma^{*}$$\times$$\Gamma^{*}
Fig. 1.7.3. Mixing diagram .
Consider a general (mixing or non-mixing) diagram with an arbitrary number of - or -insertions. Cutting the symmetry-broken propagators in their middle, thus leaving only neutral vertices, i.e. four-leg vertices with the same number of entering and exiting lines, one finds: entering lines= exiting lines, whence (considering the whole diagram): , resp. [math], resp. depending on whether the external structure is of the type , resp. or , resp. . If , then the corresponding two-by-two symmetrized matrix {1\over 2}\Big{\{}\tilde{\Pi}(\xi_{1},\xi_{2};\xi_{3},\xi_{4})+\tilde{\Pi}(\xi_{2},\xi_{1};\xi_{4},\xi_{3})\Big{\}} is therefore of the form \left(\begin{array}[]{cc}c_{1,1}&c_{1,2}e^{-2{\rm i}\theta}\\ c_{1,2}e^{2{\rm i}\theta}&c_{2,2}\end{array}\right), where are real-valued functions.
Define now to be times the sum of all connected but not necessarily two-particle irreducible Feynman diagrams with external stucture as in (1.98,1.99), but with , :
\tilde{\Sigma}(\xi_{1},\xi_{3})=\delta(\xi_{1}-\xi_{3})+\lambda\Big{\{}$$\bullet$$\bullet$$\xi_{1}$$\xi_{3}$$+$$+$$\bullet$$\bullet$$\xi_{1}$$\xi_{3}$$\times$$\times$$\bullet$$\bullet$$\bullet$$\xi_{1}$$\xi^{\prime}_{1}$$\xi^{\prime}_{2}$$\tilde{\Pi}(\xi^{\prime}_{1},\xi^{\prime}_{2};\xi^{\prime}_{3},\xi^{\prime}_{4})$$\bullet$$\bullet$$\bullet$$\xi^{\prime}_{3}$$\xi^{\prime}_{4}$$\xi_{3}
+\$$\bullet$$\bullet$$\bullet$$\xi_{1}$$\xi^{\prime}_{1}$$\xi^{\prime}_{2}$$\tilde{\Pi}(\xi^{\prime}_{1},\xi^{\prime}_{2};\xi^{\prime}_{3},\xi^{\prime}_{4})$$\bullet$$\bullet$$\bullet$$\bullet$$\xi^{\prime}_{3}$$\xi^{\prime\prime}_{1}$$\xi^{\prime}_{4}$$\xi^{\prime\prime}_{2}$$+$$\bullet$$\bullet$$\bullet$$\bullet$$\xi^{\prime}_{3}$$\xi^{\prime\prime}_{1}$$\xi^{\prime}_{4}$$\xi^{\prime\prime}_{2}$$\times$$\times$$\tilde{\Pi}(\xi^{\prime\prime}_{1},\xi^{\prime\prime}_{2};\xi^{\prime\prime}_{3},\xi^{\prime\prime}_{4})$$\bullet$$\bullet$$\bullet$$\xi^{\prime\prime}_{3}$$\xi^{\prime\prime}_{4}$$\xi_{3}
+\cdots\qquad\Big{\}}
Fig. 1.7.4. Formal expansion of the two-particle propagator in terms of two-particle irreducible diagrams.
and let be its Fourier transform, with transfer momentum entering the chain of diagrams at the left end , and exiting at the right end . Note that connecting pairs of propagators preserve momenta, so that corresponding pairs of external momenta for neighboring pairs of -kernels are equal, e.g. .
Let us now introduce the following projection procedure, a simplified procedure giving the order of magnitude of the main contributions to the Bethe-Salpeter kernel near the infra-red singularity, i.e. for a transfer momentum (of which a corrected version, called averaging procedure, is presented in §3.3 D. (ii)).
Definition 1.5** (Fermi surface projection of 2 P.I. diagrams)**
For any two-particle irreducible, four-point diagram with Cooper pair external structure, we denote by the averaged evaluation of on the Fermi circle, namely,
[TABLE]
Lemma 1.6
[TABLE]
[TABLE]
[TABLE]
[TABLE]
where , see §2.1.
Proof. The bound on derives by simple averaging from the bound in showed for . Consider now the -momentum scale contribution to or : one finds or times
[TABLE]
Briefly said, both propagators have a scale infra-red cut-off; main contribution comes a priori from the region
[TABLE]
in which the product of the two propagators by the integration volume yields . Taking into account the prefactor or , the missing prefactor comes from the constraint (1.107) on . For a more precise computation, one notes that the sector of is fixed by the sectors of . Hence
[TABLE]
The bounds for , =diag,off, follow from similar bounds for the quantities , themselves consequences of the general bounds of section 2, in particular of the subsection (§3.2) on the complementary -expansion. As in the case of or , in the case of scale inner momenta, one gets compared to the naive power-counting estimates in two supplementary small prefactors, (i) one in , ; (ii) and one in for each symmetry-broken propagator since . Multi-scale diagrams with inner momenta of scales varying between and , enjoy a supplementary ”spring factor” , allowing a multi-scale generalization of the previous argument.
Finally, the bound for , say, may be obtained along the same line by noting that \Gamma_{\phi}|\frac{d}{d\Gamma_{\phi}}C^{*}_{diag}(p)|\approx\Gamma_{\phi}\Big{|}\frac{d}{d\Gamma_{\phi}}\Big{(}\frac{1}{-{\rm i}p^{0}+e^{*}_{|\Gamma|}(\boldsymbol{p})}\Big{)}\Big{|}\lesssim\frac{\Gamma_{\phi}}{|p|^{2}+\Gamma^{2}} is bounded by a symmetry-broken propagator. Thus, using similar arguments as for , one obtains a bound in . The bound for and follows more directly from the inequality by noting that .
In order to get the leading behavior of in the infra-red limit , we first want to displace the lower external legs, namely, those on the second line, of Fig. 1.7.4, to the corresponding locations on the first line , so that intermediate ladders form bubbles inserted between -kernels, see Fig. 1.7.5 below. To this end, we need to specify both the scales and the angular sector indices , , of the lower momenta entering and exiting the -kernels, see (1.33); the Bethe-Salpeter kernel is then obtained by summing over indices of all lower momenta. Next, we rewrite lower external propagators of, say,
[TABLE]
as convolutional squares,
[TABLE]
Then we displace lower external legs to the corresponding upper locations , while taking into account the oscillations. The simplest displacement procedure, independent – contrary to the one used for renormalization in §2.2, see e.g. (2.55) – of a choice of angular sector momenta , is done in Fourier coordinates. Namely, fix the momenta of the upper external propagators , and of the lower external propagators,
[TABLE]
We may restrict to , (otherwise by construction, see Definition 2.1). Then we change integration variables in (1.111), , , so that , are now attached to the -kernel at the locations of the upper vertices . Next, we replace the resulting extra oscillations , by , and rewrite , resp. , as , resp. . The remaining oscillations may now be attached to the partially Fourier transformed kernel
[TABLE]
By momentum conservation, the above kernel depends only on , so its integral
[TABLE]
is independent of . We call it the displaced kernel. Error terms due to extra oscillations , involve quantities
[TABLE]
or conjugate, with , , , resp. , , . The derivative acting on a sector covariance or a micro-sector covariance generates (see (2.11), (2.12) or (2.15)) a small prefactor proportional to the diameter of the (micro-)sector, where or (in the case of a micro-sector) , while the displacement distance , associated with the decay factors of the two-particle irreducible diagram, produces a factor , where is the lowest internal scale of the diagram, all together an extra small prefactor
[TABLE]
playing the rôle of a spring factor in section 3, see in particular 3.2 D.
Note that the -kernels are now attached to the upper external vertices of . Contracting them with their counterparts yields pairs of -kernels connecting upper vertices,
\xi^{\prime}_{1}$$\xi^{\prime}_{3}$$C^{\alpha^{\prime}_{3}}$$C^{\alpha^{\prime}_{4}}$$\xi^{\prime\prime}_{1}$$\xi^{\prime\prime}_{3}$$\xi^{\prime}_{2}$$\xi^{\prime}_{4}$$\xi^{\prime\prime}_{2}$$\xi^{\prime\prime}_{4}
Fig. 1.7.5. Chains contributing to . Gray squares stand for -kernels.
This procedure leads to a kernel made up of kernels , alternating with pairs of propagators , which may be arrayed into the Hermitian matrix
{\cal A}_{q}(p_{3}):=\left(\begin{array}[]{cc}{\cal A}_{q}(\Upsilon_{3,diag};p_{3})&{\cal A}_{q}(\Upsilon_{3,off};p_{3})\\ {\cal A}_{q}^{*}(\Upsilon_{3,off};p_{3})&{\cal A}_{q}(\Upsilon_{3,diag};p_{3})\end{array}\right), Note that -momenta cannot be integrated over in since the -kernels also depend on them, except for the leading-order contribution to , which is proportional to the Hermitian bubble matrix of eq. (1.77).
Definition 1.7
(-kernel, preliminary version)*
Let*
[TABLE]
As stated below (1.113), the above displacement procedure neglects error terms smaller by a prefactor that goes to [math] with . This requires however angular (micro-)sectors whose thickness in all directions is comparable to ; this is explained in details in section 3.
The Fermi and -wave projections. The above description of the makes it plain that the two-particle irreducible contributions are not factorized from the bubbles, making the above expressions intractable and not ready for use in the gap equation. In order to obtain a factorization, while keeping the leading-order contribution when , we take two further steps.
- (i)
(Fermi projection) We replace external momenta by their orthogonal projection onto the Fermi circle . This defines the Fermi projected kernel ,
[TABLE]
As mentioned before, the correct procedure, described in section 3, is rather an averaging procedure in the neighborhood of the Fermi circle, but it does not matter at this stage.
- (ii)
(-wave projection) If , momenta are obtained from by a rotation of angle ; the kernel depends only on and , not separately on the angular directions of and . (For but close to [math], this is only approximately true.) Let be the angles of w.r. to a fixed direction , and (The definition makes sense also for .) Assume first that . Bubble diagrams have by construction . Hence – taking a Fourier transform – the convolution of kernels turns into a geometric series which may be resummed for each individual Fourier mode , different Fourier modes being prevented from interacting by orthogonality. From the arguments in the Introduction, it may be conjectured that only the mode (or -wave) diverges in the infra-red limit, which is itself a consequence that the interaction vertex itself is in the -wave. Thus the kernel reduces in the end to the averaged matrix . When , we introduce similarly the averaged matrix
[TABLE]
Now, bubbles may be factorized, yielding
Definition 1.8
(-kernel, or Goldstone boson propagator)* Let*
[TABLE]
where is the Hermitian bubble matrix of (1.77), which we represent for simplicity as a bubble
: graphically (compare to Fig. 1.7.5)
{\Sigma}(q)\,=\,\lambda{\mathbb{1}}\,+\,\lambda^{2}\ \Big{\{}••••••+\ \cdots\ \Big{\}}(q)
Let us consider specifically the case of zero transfer momentum . Then, replacing by in (1.115), the -wave projection is exact, and defines
[TABLE]
where and
[TABLE]
First-order terms are
\lambda{\mathbb{1}}\,+\,\lambda^{2}$$\bullet$$\bullet$$+\,\lambda^{3}\,$$\bullet$$\bullet$$\bullet$$+\,\lambda^{4}\Big{(}$$\bullet$$\bullet$$\bullet$$\bullet$$\bullet$$\bullet$$+\,$$\bullet$$\bullet$$\bullet$$\bullet$$\Big{)}
+\,\lambda^{5}$$\Big{(}$$\bullet$$\bullet$$\bullet$$\bullet$$\bullet$$\bullet$$\bullet$$+\,$$\bullet$$\bullet$$\bullet$$\bullet$$\bullet$$+\,$$\bullet$$\bullet$$\bullet$$\bullet$$\bullet$$+$$\bullet$$\bullet$$\bullet$$\bullet$$\bullet$$\Big{)}
where propagators are either diagonal or off-diagonal.
The one-point counterterm coefficient is chosen to be the solution with largest module of the
Definition 1.9** (gap equation)**
[TABLE]
i.e. \left(\begin{array}[]{c}\Gamma^{\perp}\\ (\Gamma^{\perp})^{*}\end{array}\right) is a null eigenvector of the matrix . Recall . Graphically,
-$$\Gamma$$\Sigma_{\perp,\perp}$$q=0$$+$$q=0$$\Sigma_{\perp,\perp}••••
Fig. 1.7.6. Gap equation.
By the arguments below Fig. 1.7.3, the matrix is of the form \left(\begin{array}[]{cc}a&-b\omega\\ -b\omega^{-1}&a\end{array}\right) with if . Hence the matrix can be made real by the conjugation
[TABLE]
a condition which fixes mod . Then eigenvectors \left(\begin{array}[]{c}|\Gamma|\,e^{{\rm i}(\theta+\frac{\pi}{2})}\\ |\Gamma|e^{-{\rm i}(\theta+\frac{\pi}{2})}\end{array}\right), resp. \left(\begin{array}[]{c}|\Gamma|\,e^{{\rm i}\theta}\\ |\Gamma|e^{-{\rm i}\theta}\end{array}\right) are eigenvalues for two different values, and , of obtained by solving the equations , or equivalently (conjugating as in (1.122) to reduce to the case )
[TABLE]
for , yielding two solutions, . Using previous computations, see (1.56), (1.53), one obtains leading-order terms as functions of the non-dimensional variables :
[TABLE]
with
[TABLE]
[TABLE]
see Lemma 1.6, with . The coefficients that we introduced – and the coefficients that we introduce below – are . Finally, . The dependence on of and is secondary since vanishes to all orders at , implying in particular that is essentially diagonal. The change of variable allows one to rewrite the function as
[TABLE]
whence
[TABLE]
[TABLE]
with provided is small enough. Thus, by the implicit function theorem,
[TABLE]
with
[TABLE]
In particular,
[TABLE]
We note for further use that, choosing , one gets modulo the conjugation (1.122), or for ,
[TABLE]
with , implying a null eigenvector proportional to \left(\begin{array}[]{c}-{\rm i}\\ {\rm i}\end{array}\right).
Thus the net effect of diagrams other than the Cooper pair bubble is to change by a finite amount the pre-factor in front of and produce a finite multiplicative splitting between two solutions.
Estimates for near . We first need a Taylor expansion of near . Lemma 1.3 and power-counting estimates for imply that
[TABLE]
=diag,off, where \left(\begin{array}[]{c}Q_{diag}\\ Q_{off}\end{array}\right) is a couple of positive-definite quadratic forms, , Q_{off}(q)=-\big{(}\frac{1}{g_{\phi,off}}(q^{0})^{2}+\frac{v^{2}_{\phi,off}}{g_{\phi,off}}|\boldsymbol{q}|^{2}\big{)}, see Lemma 1.3, with
[TABLE]
Remaining terms involving have an extra prefactor. Using (1.133), one gets for
[TABLE]
where
[TABLE]
and by Definition 1.4. At this point it is natural to shift to the -basis (compare with §1.5), in which =diag, with
[TABLE]
[TABLE]
Multiplying by yields (see (1.118), (1.124), (1.132))
[TABLE]
and
[TABLE]
Connection to pre-kernel. Far enough from , it is possible to expand further the denominator: letting , and representing by
the pre-kernel as in the previous subsection, we can write down the geometric series
\frac{\lambda}{1-\Pi(q){\cal A}(q)}\ =\$$+\,$$\frac{1}{\lambda}{\cal A}(q)R(q)$$+\ \cdots
However, the above expansion is faulty for (i.e. in our Regime B.), since there (replacing by the kernel , of which it is assumed to be an approximation)
2 Fermionic theory
The section is organized as follows. Recall . We concentrate exclusively on momenta with scales .
First (§2.1), we refine the multi-scale analysis of §1.3 by decomposing scale propagators according to angular sectors (see Definitions 2.1 and 2.2).
Then (§2.2), we introduce a dressed action , see Definition 2.6, which is the basis of all subsequent computations. This action is dressed by scale dependent -parameters which enact the scale decoupling necessary to separate local parts of two-point functions and renormalize. We take the opportunity to give qualitative statements about the renormalization flow of parameters.
With at hand, it is possible to present the general cluster expansion scheme (§2.3), involving scale-by-scale horizontal expansions, a bubble resummation, a vertical momentum-decoupling expansion, the separation of local parts of divergent diagrams, and a Mayer expansion. On the way, we present a preliminary, single-scale version of the bounds for polymers.
Finally, general bounds for multi-scale fermionic polymers, based on the preliminary bounds proved in §2.3, and confirming at the same time the predictions of §2.2 about the renormalization flow, are proved in §2.4.
2.1 Angular decompositions
We continue here the analysis of §1.3 by decomposing the -th momentum shell into spatial sectors of aperture angle , where
[TABLE]
Technically, one uses a smooth partition of unity
[TABLE]
where the support of intersects the Fermi circle roughly on the circular arc
, in such a way that the distance dist() , with . For definiteness, we choose (defined below, see Definition 2.2) in such a way that , where
Definition 2.1** (angular sectors)**
- (i)
[TABLE]
[TABLE]
[TABLE]
where are the polar coordinates of , namely, \frac{\boldsymbol{p}}{|\boldsymbol{p}|}=\left(\begin{array}[]{c}\cos\theta(\boldsymbol{p})\\ \sin\theta(\boldsymbol{p})\end{array}\right).
- (ii)
* and are the spatial projections of .*
The number of angular sectors is chosen in such a way that sectors are essentially isotropic: they are small deformations of 3-dimensional cubes of side scaling like . Note that the momentum scale and angular sector of a given spatial momentum is defined up to , namely, if , , then
Remark. As discussed in ([25], Lemma 3) in details, see also Proposition 2.5 below, this choice of isotropic sectors is not optimal, and produces spurious (but not really disturbing) prefactors when resumming over angular sectors. However, the choice of curvelet-like anisotropic sectors (see e.g. [10, 11, 69]) with angular width instead of – taking into account the curvature of the Fermi circle – would have led to a tedious wavelet-like phase-space analysis which we want to spare to the reader.
Definition 2.2
Let
[TABLE]
[TABLE]
Let us rewrite (1.33),(1.34) for the sake of the reader,
[TABLE]
[TABLE]
For every , one chooses a vector lying inside the support of , and lets . As explained in §1.3, see discussion in the paragraph between (1.27) and (1.29), all fermions of scale may actually be considered as (non-optimal) decompositions of scale fermions.
Remark. In section 3, we shall need to further decompose scale sectors , , into smaller angular sectors of angle aperture , with ; we write
[TABLE]
and define accordingly angular sectors roughly equal to the support of the cut-off functions , and covariance kernels , as in (2.9). Note that such angular sectors are strongly anisotropic if . If necessary, they may be further ”chopped” into isotropic ”microsectors” , with size in the direction transverse to the Fermi circle. We then let be some momentum lying inside , and define covariance kernels with characteristic functions instead of .
Proposition 2.3
For every , there exists a constant such that, for every scale , sector and multi-index ,
- (i)
if ,
[TABLE]
- (ii)
if ,
[TABLE]
Eq. (2.11) is proved in ([25], Lemma 4). If and , then the factor is equal to the scaling factor of , times the scaling of the volume of the sector, . The supplementary prefactor comes from the Fourier multiplier associated to the operator , since is bounded over the support of by the diameter of the sector, which is proportional to . Assume then that , so that , \Big{|}\frac{1}{{\rm i}p^{0}-e^{*}_{|\Gamma|}(\boldsymbol{p})\sigma^{3}-{\mathbb{\Gamma}}(\theta)}\Big{|}\lesssim\Gamma_{\phi}^{-1}\approx(2^{-j_{\phi}}\mu)^{-1}. Hence . The summable spring factor shows clearly that this further decomposition would be useless for the free theory. The decay rate is only instead of the expected fermion mass because of the scale cut-offs, which are thus clearly not optimal.
In other words: the average magnitude of – defined as the square-root of for nearby – is
[TABLE]
A similar computation – with a supplementary relative factor coming from integration in Fourier coordinates – holds for the sector-symmetric covariance kernel , yielding
[TABLE]
Remark. The above Proposition generalizes to microsector covariances (see previous Remark),
[TABLE]
if () and Namely, compared to (2.12), the volume of the microsector scales now like , and its size along any direction like .
Boxes. For every , we decompose into a disjoint union of scale boxes defined as follows.
Definition 2.4** (boxes)**
A scale box is a ”cube” , . We denote by the set of scale boxes, so that (up to a subset of measure zero) defines a partition of . Let finally .
In order not to spoil completely momentum preservation at vertices, we further introduce for each a smoothened partition of unity, , where each is a smooth, compactly supported function and :
• (scaling property) if , ;
• (translation invariance) if \boldsymbol{v}=2^{j}\left(\begin{array}[]{c}i_{0}\mu^{-1}\\ i_{1}(p^{*}_{F})^{-1}\\ i_{2}(p^{*}_{F})^{-1}\end{array}\right) is a scaled integer translation, then ;
• (support) except if belongs to or its (direct of diagonal) neighbors;
• in some neighbourhood of the center of .
Then (by the scaling and translation invariance properties), for every multi-index ,
[TABLE]
The volume of a scale box is roughly inverse to that of a sector,
[TABLE]
The bound (2.16) implies a bound for the Fourier transform with a quasi-exponential decay,
[TABLE]
for every .
Sector counting. A characteristic feature of dimension 2 is that (generically, as we shall see) four spatial momenta in the immediate neighborhood of the Fermi circle and such that are essentially two-by-two equal, namely, or , by which we mean (letting be the angular sector of , ) that, to leading order,
[TABLE]
as illustrated by the following picture
\Sigma^{*}_{F}$$\boldsymbol{p}_{1}$$\boldsymbol{p}_{2}$$\boldsymbol{p}_{4}$$\boldsymbol{p}_{3}$$\theta
Fig. 2.1.1. Close to the Fermi circle, there are only two free sectors per vertex in dimension 2.
Thus, splitting each of the four fields of a vertex into its -components leads for ”quasi-scale-diagonal” terms in a near-Cooper pairing,
[TABLE]
with , to essentially only two sums over angular sectors. More precisely:
Proposition 2.5
(see [25], Lemma 3).
- (i)
Let , be the spatial projection of four angular sectors of scales and . Then the set
[TABLE]
is empty unless there exists such that
[TABLE]
- (ii)
(generalization) Let , be the spatial projection of angular sectors of scales . Then the set has cardinal .
Sketch of proof for (i) (case of a vertex). The condition ensures that the aperture angle is smaller than the relative pair angle . Since , we deduce that . In other words, is the sector closest to , and the pair is essentially obtained from through a rotation of angle . (Permuting with or yields two equivalent possibilities). Using the approximate equality, , and projecting onto the axis perpendicular to yields , whence
The most important case is , for the above Proposition gives then the sector counting factor of a vertex: accepting an error on the momentum preservation condition, one has possibilities for the choices of the sector indices ; namely, for fixed , there are choices for , possibilities for , and finally, possibilities for given . Summing over yields the result. We call generic a vertex for which the relative pair angle is large, which implies that the aperture angle is small, of order ; diagrams with only generic vertices have the following very important property: vertices may be cut into two parts in such a way that the transfer momentum between the half-vertices is very small. In (2.20), half-vertices are and . Sectors are very close. Splitting all vertices in this way yields fermion loops on which all sectors are very close, so that there is essentially only one sum over sectors per loop. This gives birth to the ”1/N” factors. This approach is pursued in a systematic way in §2.3.2. In this ideal case, the -factor is . Taking into consideration the opposite degenerate case when is closer to than to , the small factor is really (at worst) .
Alternatively, assuming to be fixed, this leaves only possibilities. Later on, we shall see that (at least perturbatively) sectors are in average shared by two vertices, so fixing one of the sectors is the correct way of counting; thus, for a vertex , there is – up to the logarithmic correction – only one sum over sectors. Also, as mentioned in the Remark just above Definition 2.2, choosing anisotropic sectors with angular width instead, the logarithmic correction disappears, so there is exactly one sum over sectors.
An immediate corollary is the following. Let
[TABLE]
, with fixed, and as above. Then, taking a partial Fourier transform w.r. to spatial coordinates, and considering a family of smooth, scaled Fourier cut-offs such that and except if , so that (in particular) for all (compare with the bounds of Proposition 2.3), we get
[TABLE]
where has cardinal for ; in practice, we restrict the sum in (2.24) to and replace without further mention by . Then
[TABLE]
Taking into account the scaling of each field, one obtains finally .
In other words: The vertex, i.e. the quartic term in the interaction, integrated over a scale box , fixing one of the sectors,
[TABLE]
is small – up to a logarithmic correction –, as follows from our Assumption on the coupling constant .
We shall see in §2.4 how to get rid of this logarithmic correction in diagram bounds. Note that there is no logarithmic correction as soon as .
2.2 Dressed action
The principle of multi-scale expansions is to rewrite the bare action of Definition 1.1 as a sum over of effective actions with highest scale , and renormalized, scale-dependent coefficients .
Field dressing. Recall . Given a function , we let ,
[TABLE]
with
[TABLE]
and
[TABLE]
The sum in (2.28) contains only at most terms by the support condition of . This definition implies
[TABLE]
The above expression vanishes except in some neighborhood of . The adjoint field is dressed similarly.
Let
[TABLE]
Definition 2.6** (dressed action)**
Let
[TABLE]
where
[TABLE]
[TABLE]
[TABLE]
and, for ,
[TABLE]
where , .
Compare with Definition 1.1. The principle of the dressing is the following:
• If , then , i.e. one retrieves the bare theory (see (2.41) below);
• If for all , then is the tensor product of two independent measures, a measure for the high-momentum fields , and a measure for the low-momentum fields . On the other hand, derivatives , acting on produce by construction either (i) mixed terms involving both high-momentum fields and low-momentum fields ; or (ii) subtracted terms or (T\bar{\Psi})^{\to(j+1)}\Big{\{}(\delta Z^{j}-\delta Z^{j+1})\partial_{\tau}-\Big{(}(\delta m^{j}-\delta m^{j+1})\frac{|\boldsymbol{\nabla}|^{2}}{2m^{2}}+(\delta\mu^{j}-\delta\mu^{j+1})\Big{)}\sigma^{3}+(\delta{\mathbb{\Gamma}}^{j}(\theta)-\delta{\mathbb{\Gamma}}^{j+1}(\theta))\Big{\}}(T\Psi)^{\to(j+1)}. Subtracted terms (ii) associated with parameter renormalization are there precisely to compensate mixed terms of type (i).
Looking at the theory at scale , and combining the covariance with the quadratic term in , we obtain an effective fermion covariance
[TABLE]
where
[TABLE]
[TABLE]
[TABLE]
In particular, when , one retrieves the bare parameters,
[TABLE]
*We determine , , and simultaneously , by a fixed point argument in such a way that *
[TABLE]
and the pre-gap equation (see Definition 1.2) holds. As emphasized in the Introduction, this procedure implies that the effective radius of the Fermi circle is .
As we shall see (see §2.3.2), the order of divergence of a lowest scale amputated diagram with external fermion legs is above symmetry-breaking momentum scale , which means that scales at most like . Thus diagrams with external legs must be renormalized. This power-counting estimate is not necessarily true for the amplitude of a polymer, which suffers from non-overlapping constraints; we shall therefore be led to renormalize diagrams with up to external legs, whose local parts may be freed from non-overlapping constraints by the Mayerization procedure described in §2.3.5. In this theory, one renormalizes with external legs on the effective Fermi circle where the IR singularity lies. By symmetry, non-vanishing -point functions must have external -fields and external -fields, with even. So we need only discuss two-, four- and six-point functions. Counterterms are associated to two-point functions only; they require a precise analysis, which will lead to the renormalization flow.
(two-point functions) The two external legs – , or – are in the vicinity of the Fermi sphere. Let be a unit vector, and its image through the rotation of angle . Expanding with respect to an external momentum in space direction , with , , and using , one gets
[TABLE]
and similarly (summing over all diagrams with internal legs of momentum scales ), letting be as in §1.1 the Fourier transform of the two-point function ,
[TABLE]
The difference selects the contribution of diagrams of lowest internal scale precisely equal to . Second-order derivatives need not be renormalized because the effective order of divergence is . By rotation invariance, \frac{d}{dp_{//}}\delta G^{j}(p)\Big{|}_{p=(0,p^{*}_{F}\boldsymbol{i})}=0.
Let us discuss to begin with diagonal components , . First, and \frac{d}{dp_{\perp}}\delta G^{j}_{i,i}(p)\Big{|}_{p=(0,p_{F}^{*}\boldsymbol{i})} yield spin-independent counterterms for the chemical potential and for the mass. Finally, \frac{d}{dp^{0}}\delta G^{j}_{i,i}(p)\Big{|}_{p=(0,p_{F}^{*}\boldsymbol{i})} yields a -counterterm .
The off-diagonal part \delta G^{j}_{off}=\left(\begin{array}[]{cc}0&\delta G^{j}_{1,2}\\ \delta G^{j}_{2,1}&0\end{array}\right) of involves an odd number of off-diagonal propagators with matrix coefficient or , , whence (as shown by induction on ) one extra prefactor , thus lowering the effective degree of divergence of diagrams by . Consequently, we need only Taylor expand to order [math],
[TABLE]
It follows crucially from the fact that the scale infra-red cut-off theory, see Definition 1.1 (ii), retains the symmetries of the original theory, so that (by (1.10) and (1.11))
[TABLE]
[TABLE]
[TABLE]
[TABLE]
Hence: , and are all real parameters, while is a Hermitian matrix of the form \delta{\mathbb{\Gamma}}^{j+1}-\delta{\mathbb{\Gamma}}^{j}\equiv\left(\begin{array}[]{cc}0&\bar{\delta\Gamma}^{j+1}-\bar{\delta\Gamma^{j}}\\ \delta\Gamma^{j+1}-\delta\Gamma^{j}&0\end{array}\right), with .
Remark. Note that diagrams contributing to -counterterms are not necessarily in the -direction () because of mixing diagrams. However, the leading-order contribution is a tadpole, see Fig. 2.2.4 (iv), which is in the -direction. Furthermore, a Ward identity (see §3.3) ensures that the local part of the fermionic two-point function vanishes in the direction perpendicular () to .
Summarizing:
Definition 2.7** (counterterms)**
Let be a unit vector. The following definitions do not depend on the orientation of :
- (i)
(chemical potential counterterm) let
[TABLE]
- (ii)
(mass counterterm)
[TABLE]
- (iii)
(energy gap counterterm) let
[TABLE]
- (iv)
(* counterterm) let*
[TABLE]
Subtracting counterterms is equivalent to subtracting local parts of diagrams, which is done as follows. Let
[TABLE]
be the amplitude of in direct space. The external leg at is paired to a field or , with , giving rise to an external covariance kernel . The direct space renormalization algorithm proceeds by applying Taylor’s formula to the above kernel, so as to displace its attachment point from to . For diagrams contributing to , we need to expand to order 1 so as to lower the order of divergence by two. For diagrams contributing to , a simple subtraction is enough:
- (i)
(off-diagonal diagrams) Setting apart only the leading-order term, we get
[TABLE]
a sum of a local part as in (2.45), and of a renormalized amplitude involving the scalar product is manifest. Since the lowest internal scale of is , distances between vertices contribute a scaling factor to the expected power-counting, as can be shown using part of the decay factors in the covariance (2.11); this is a standard argument, see e.g. [54]. On the other hand, see again (2.11), the gradient operator produces a small factor . All together, we have obtained a small spring factor , having the effect of lowering the degree of divergence of the diagram from [math] to .
- (ii)
(diagonal diagrams) For such diagrams, we need to set apart the two leading-order terms,
[TABLE]
where subtracts the local part of the diagram to order 1. Compared to (i), the renormalized amplitude involves now the square of the scalar product , hence a spring factor , having the effect of lowering the degree of divergence of the diagram from to . 2. 2.
(-point functions with ) A simple subtraction (as for off-diagonal two-point diagrams) is required. Fixing one of the external vertices , we rewrite \Big{(}\prod_{i=2}^{N_{ext}}\int d\xi_{ext,i}\Big{)}{\cal A}^{j\to}(\Upsilon;\boldsymbol{\xi}_{ext})\ \prod_{i=2}^{N_{ext}}C^{k_{i},\alpha_{i}}(\xi_{ext,i},\xi^{\prime\prime}_{i}) as a local part,
\Big{\{}\Big{(}\prod_{i=2}^{N_{ext}}\int d\xi_{ext,i}\,e^{{\rm i}(p^{k_{i},\alpha_{i}},\xi_{ext,i}-\xi_{ext,1})}\Big{)}{\cal A}^{j\to}(\Upsilon;\boldsymbol{\xi}_{ext})\Big{\}}\ \prod_{i=2}^{N_{ext}}C^{k_{i},\alpha_{i}}(\xi_{ext,1},\xi^{\prime\prime}_{i}), plus an error term which may be written in integral form as in (i) above.
Four-point diagrams are in principle logarithmically divergent. However, as we shall prove in §3.2, the only logarithmically divergent contributions come from Goldstone boson insertions, which are explicitly resummed into the -kernel. So no counterterms need to be introduced for four-point functions: the coupling constant keeps its bare value throughout.
Leading-order diagrams in or are obtained as follows,
(i)
-$$\bullet$$\bar{\psi}$$\psi$$\bar{\psi}$$\psi
Fig. 2.2.1. Leading contribution to chemical potential renormalization .
(ii)
-\frac{(m^{*})^{2}}{p_{F}^{*}}\frac{d}{dp_{\perp}}$$\psi$$\bar{\psi}$$\bullet$$\xi$$\bullet$$\xi^{\prime}$$\beta$$\alpha_{1}$$\alpha_{2}$$\alpha_{3}
Fig. 2.2.2. Leading contribution to mass renormalization .
(iii)
\psi$$\bar{\psi}$$-{\rm i}\frac{d}{dp^{0}}$$\bullet$$\bullet
Fig. 2.2.3. Leading contribution to renormalization .
(iv)
\times$$\Gamma^{j}$$\bar{\psi}$$\bar{\psi}$$\psi$$\psi$$\bullet
Fig. 2.2.4. Leading contribution to gap renormalization .
We denote by and their respective values. We can now restate Theorem 1 of the Introduction.
Lemma 2.8
(see Theorem 1) Assume . Then:
For , the renormalization flow is to leading order
[TABLE]
[TABLE]
[TABLE]
[TABLE]
For , one has instead
[TABLE]
[TABLE]
In particular (summing over scales )
[TABLE] 2. 2.
Furthermore, for an adequate choice of , such that
[TABLE]
the scale mass and chemical potential counterterms vanish, more precisely,
[TABLE]
and the pre-gap equation, see Definition 1.2, holds. Also,
[TABLE]
Proof.
We shall be content here with computing the values of ; the arguments of §1.3 show that more complicated diagrams give contributions which are smaller by a factor . The evaluation of the diagrams of Fig. 2.2.1,2.2.2,2.2.3,2.2.4 (see [25] for details) is based on the decomposition (1.45) of diagonal propagators, implying
C^{*}_{\theta}\approx\left(\begin{array}[]{cc}\frac{1}{-{\rm i}p^{0}+e^{*}(\boldsymbol{p})}&\\ &\frac{1}{{\rm i}p^{0}+e^{*}(\boldsymbol{p})}\end{array}\right)\approx\left(\begin{array}[]{cc}\frac{1}{-{\rm i}p^{0}+p_{\perp}}&\\ &\frac{1}{{\rm i}p^{0}+p_{\perp}}\end{array}\right) for momenta of scale . Recall that the integration measure averaged over the angular coordinates (i.e. summed over sectors) is . We do the computations for ; for larger , one need just remark that the covariance kernel (2.12) has a supplementary prefactor.
- (i)
(chemical potential) Evaluate the amputated diagram over momenta . By symmetry, the result is real. Using , one gets In particular, .
- (ii)
(mass) As in §5.4, the best way to get a correct order of magnitude of is to choose anisotropic sectors . If , then the main contribution to the integrand comes from sectors , roughly parallel to . Thus: the sum over sectors produces a factor ; the product , a factor ; the derivative , a factor ; the integral over , a factor , i.e. the volume of an anisotropic box. Whence all together, a contribution . The integral over yields a contribution smaller by a factor , since: each propagator is a instead of , see (5.28); the integration volume is instead of , while (for fixed) the sum over sectors , , reduces essentially in fact to a single sum, yields a supplementary factor only (see Proposition 2.5, with logarithmic factor removed for anisotropic sectors, following the Remark just above Definition 2.2).
- (iii)
(-coefficient) Similar to (ii).
- (iv)
(-counterterm renormalization) Main scale diagram is a tadpole involving the symmetry-broken propagator,
\times$$\Gamma$$\bullet$$\bullet$$\xi$$j$$j^{\prime}>j
As already noted in §1.5, see Fig. 1.5.1, this diagram involves to leading order the same integral as the Cooper pair bubble diagram , and is evaluated (see (1.56)) as . 2. 2.
(fixed point) Summing over scales , one gets
[TABLE]
[TABLE]
[TABLE]
(see proof of 1.(iv) for the last equality), where are non-dimensional constants, , and , , for any such that . Eq. (2.69), combined with (1.62), imply
[TABLE]
for such values of . Then the implicit function theorem yields for small enough a unique solution to the system of equations complemented with the pre-gap equation (1.65). Furthermore, they satisfy (2.64) and (2.66).
2.3 Cluster expansions
The art of cluster expansions for a general renormalizable theory consists in resumming separately the leading contributions of diverging -point functions, and absorbing them into a renormalization of the measure. By ”renormalization of the measure”, we intend a change of normalization by a constant prefactor (which disappears when one evaluates connected correlations), and a shift of its parameters, here and . This is done in an organized way by evaluating -point functions in terms of a sum over clusters or polymers, which themselves may be expanded perturbatively into an infinite series of diagrams involving only fields located inside the image of the polymer. Clusters are arrays of scaled boxes connected by a tree of links – the linking consisting of propagators between two boxes of . Polymers are synonymous for clusters, or multi-scale versions of them, involving both horizontal links connecting boxes of the same scale, and vertical links capturing the interaction between scales through multi-scale vertices. Here is a graphical representation of a multi-scale vertex , with the four half-propagators connecting it to other vertices; the scale of the vertex is by definition the scale of the half-propagator with highest momentum, here .
\Delta^{j}_{\xi}$$\Delta^{k}_{\xi}$$\Delta^{k^{\prime}}_{\xi}$$j$$k$$k^{\prime}•\xi$$\xi•
**Fig. 2.3.1. Two multi-scale vertices of scale , from left to right:
, and . By assumption .**
Horizontal cluster expansions (one per scale) produce horizontal links. Vertical cluster expansions produce vertical links, and single out in particular divergent polymers, which are polymers with a small (here ) number of external legs. The renormalization step resums the leading term – called local part – of divergent polymers into a scale-dependent redefinition of the parameters of the interaction. As a minor correction, a further cluster-like expansion (called Mayer expansion) is required to get rid of the non-overlapping conditions between boxes belonging to two different polymers, so as to regain translation invariance for divergent polymers. So far for fundamentals.
We start from scale and apply inductively the following sequence of fermionic expansions, down to a scale ,
(Horizontal cluster expansion of scale ) (Resummation of (lowest) scale chains of bubbles) (Momentum-decoupling expansion of scale ) (Displacement of external legs and Mayer expansion of scale )
(Horizontal cluster expansion of scale ) (Resummation of lowest scale chains of bubbles) (Momentum-decoupling expansion of scale ) (Displacement of external legs and Mayer expansion of scale )
(Horizontal cluster expansion of scale ) (Resummation of lowest scale chains of bubbles) (Momentum-decoupling expansion of scale ) (Displacement of external legs and Mayer expansion of scale )
Once this program has been completed, we shall be led to apply in the next section (see §3) a sequence of complementary horizontal or bosonic expansions of bosonic scales ,
(Horizontal and Complementary horizontal expansions of scale ) (Complementary horizontal expansion of bosonic scale ) (Complementary horizontal expansion of bosonic scale )
with scales ranging from to . Scale is very particular, since purely fermionic connectivity issues (see introduction to this subsection) mingle with the necessity of factorizing Goldstone boson propagators to show infra-red summability of the expansion. Therefore the scale horizontal cluster expansion is postponed to section 3.
Apart from the intertwining sequence of bubble resummations, the sequence of fermionic expansions is fairly standard, see e.g. [70] for a review. Through the horizontal and momentum-decoupling expansions, -point functions are rewritten as sums of products of polymer-dependent expressions of the type . A polymer is by definition a union of boxes in with the following properties:
(i) two boxes may or may not be connected by a horizontal link produced by the scale horizontal cluster expansion. The set of all boxes in , connected by horizontal links, is a cluster forest, i.e. a disjoint union of a finite number of cluster trees;
(ii) two boxes such that , i.e. lies just below , may or may not be connected by a vertical link produced by the scale momentum-decoupling (or vertical) expansion. The set of all boxes in , connected by cluster links, i.e. horizontal and vertical links, is connected.
The set of all such polymers with lowest scale will be denoted by . Furthermore, if , then
(iii) (external structure of a polymer) there exists a subset of scale boxes vertically connected to the boxes in lying just below them.
The function F(\mathbb{P})=F(\mathbb{P};\Psi^{\to j_{max}}\Big{|}_{{\bf\Delta}^{j_{max}}},\bar{\Psi}^{\to j_{max}}\Big{|}_{{\bf\Delta}^{j_{max}}}) is obtained by Gaussian integrations by parts, which produce a number of vertices, each lying in one of the boxes contained in ; say, vertices in . If are connected, then contains a propagator , with . Furthermore, contains at least one external field or per external box , namely, \Big{|}\frac{\delta F(\mathbb{P};\Psi^{\to j_{max}}|_{{\bf\Delta}^{j_{max}}},\bar{\Psi}^{\to j_{max}}|_{{\bf\Delta}^{j_{max}}})}{\delta\Psi^{\to j_{max}}(\xi)}\Big{|}+\Big{|}\frac{\delta F(\mathbb{P};\Psi^{\to j_{max}}|_{{\bf\Delta}^{j_{max}}},\bar{\Psi}^{\to j_{max}}|_{{\bf\Delta}^{j_{max}}})}{\delta\bar{\Psi}^{\to j_{max}}(\xi)}\Big{|}\,\not=0 if . The number of external fields of is denoted by .
In principle, the sum over all polymers can be bounded using perturbative arguments as follows. First we fix the location of one of the lowest-scale boxes, say, , in order to avoid introducing diverging volume factors proportional to , and let be the set of polymers in such that . In practice, may be seen as the approximate location of one of the external legs of the -point function. Then the sum is controlled by two types of factors,
– local factors, for each vertex, in particular, for each box;
– geometric factors due to the fast-decaying propagators .
Geometric factors make it possible to sum over all possible geometries and vertex number assignations for a given number of vertices, say, . Local factors ensure that the sum over is bounded by a geometric series of the type .
Things are however not so simple. We have overlooked in the above quick overview extra volume factors due to the necessity, for fixed, of summing over all lying above ; this necessity comes from the fact that only the locations of the external legs of the polymer (or, more or less equivalently, of ) are fixed. If there is a vertex in , then this can be remedied by simply replacing by a constant times . However, it may well be that successive vertical links connect to through boxes void of any vertex. The sum over such that becomes an arbitrarily large volume factor as the scale difference increases, eventually swallowing up the small factor . However, multi-scale diagrams extending over several scales can be shown by simple power-counting arguments to be smaller than expected by some similar, exponential factor , where depends on the number of external legs. This is typically a situation which calls for renormalization; if there are enough links between the part of lying above and , then the vertex power-counting beats the volume factor because ; in this theory links are a priori enough, but for technical reasons (due to the complexity of the sum over angular sectors of external legs, see §2.4 below), we shall follow the same procedure whenever . When , the renormalization procedure described in §2.2 should be applied. Subtracting the local part of a polymer is equivalent to displacing all its external legs to the location of one of them. Error terms (polymer) - (local part) have the same relative power-counting as if they had external legs, whence they are indeed considered as if they did have external legs. On the other hand, local parts are resummed by defining scale-dependent constants , at the price of some easy partial resummation of the series of polymers, called Mayer expansion, which takes into account the non-overlapping constraints between boxes, and whose outcome is a sum over Mayer polymers. At this stage, it is enough to know that the Mayerization step involves a ”multicolor” version of the theory, with one color per polymer. Then previously developed links connect boxes of the same color, while Mayer polymers are multicolor polymers connected by overlap links between different colored versions of the same two boxes. Painting all boxes gray, Mayer polymers are then polymers in the usual sense, with some extra overlap links connecting a box to itself. Summing over Mayer polymers in , one obtains e.g. for the theory restricted to scale
[TABLE]
in the sense that , where (fixing some box in an arbitrary way) is the volume of , defined as the number of scale boxes lyind inside , and is interpreted as the scale free energy density. Going further down the scales, the general task is to prove that , where – similarly interpreted as a scale free energy density – is obtained by summing over Mayer polymers in .
In general, the coupling constant should also be made scale-dependent, following the above local part resummation pattern. However, the specificity of this model (as explained in the Introduction) is that the ”local part” of the main diagrams contributing to the renormalization of , making up the Bethe-Salpeter kernel, must be resummed nonperturbatively as a series. This strongly model-dependent, nonperturbative step is performed separately, by expanding successively Cooper pairs with transfer momentum of scales downto the lowest infra-red scales. It will be described in the section about the bosonic regime (section 3).
Downto scale , however, namely, in the present section, much simpler arguments suffice. Namely, one may systematically expand chains of bubbles. Looking more specifically e.g. at the first resummation step of the sequence of fermionic expansions, chains of bubbles with scale fermionic propagators are first considered as contracted external legs of scale fermionic diagrams, and resummed systematically before expanding external legs of scale ; in some sense, this procedure may be seen as a expansion in itself, of scale intermediate between and . Similarly, the resummation of lowest scale chains of bubbles may be considered as an expansion step of scale .
The reader will have to wait until §2.4 to have a comprehensive image of how polymers of lowest scale are bounded. It is however not straightforward to understand how all arguments fit together. Therefore, we shall first present a single-scale bound for the outcome of the horizontal cluster expansion in §2.3.1. The bounds of §2.4 cover the whole fermionic regime down to scale and take into account bubble resummations.
2.3.1 Horizontal cluster expansion
Let be a family of functions , with
[TABLE]
extended trivially to the diagonal by letting . (Later on, we shall use slightly generalized functions, defined on pairs of clusters, where clusters are unions of boxes of the same scale). As we shall see, the set of coefficients for fixed, called scale weakening coefficients, defines an interpolation between the initial covariance and the box-diagonal covariances for which fields lying in different boxes have been made independent, namely, if , .
Index notations: it is often useful or simply convenient to use different indices for external - and -fields; thus we shall often index the -fields and their coordinates by an index ranging in some abstract set , and the -fields and their coordinates by an index ranging in some abstract set , where is a copy of . On the other hand, on other occasions it will be more convenient to use ”neutral” notations which do not distinguish between ’s and ’s; for instance, if , , , then , .
Horizontal cluster expansion of scale . In order to compute the connected -point function
[TABLE]
one applies a Brydges-Kennedy-type expansion, see Appendix, §5.2, to the numerator
[TABLE]
where is the free Grassmann measure on with weakened two-by-two matrix covariance kernel
[TABLE]
[TABLE]
where are sets of weakening coefficients which (as their name indicates) weaken the covariance kernels when lie in different boxes. The function is implicitly trivially extended to the diagonal by letting . The unnormalized -point functions interpolate between the original -point function and the unnormalized -point function computed in a totally decoupled theory. The BKAR formula (§5.2, Proposition 5.2) yields
[TABLE]
Induction on . Let , resp. be the set of Mayer polymers with lowest scale and , resp. external fields. We assume by induction on that has been rewritten as
[TABLE]
where
[TABLE]
and is expressed in terms of a sum over Mayer polymers extending down to scale ,
[TABLE]
The expression involves a sum over all possible assignments of the external field indices to the various polymers, so that ; the function is obtained by contracting in a certain way the external fields with the interaction or between themselves. The function depends (polynomially) on the fields through the external structure of . By assumption only Mayer polymers with external fields are considered in the sum; the reason is that
(i) vacuum polymers (i.e. polymers with ) have been resummed into the scale free energies;
(ii) two-point polymers (i.e. polymers with ) have been absorbed into the renormalization of the parameters ;
We shall not decouple boxes vertically connected to two boxes of scale belonging to the same polymer , since they are connected from above. Identifying all such sets of boxes defines a set of objects , which are unions of boxes in . Thus it is enough to introduce a weakening function , with . Replace in (2.79) the integral with
[TABLE]
Horizontal cluster expansion of scale . Start from the expression (2.81), and proceed as we did at scale :
[TABLE]
The above scale horizontal cluster expansion has produced step polymers in , described as follows:
– is a union of the polymers in obtained at the end of the scale perturbative expansions;
– is a union of clusters belonging to .
More precisely: two step polymers are connected after step if and only if there exist two horizontally connected boxes (i.e. boxes belonging to the same step polymer) and two boxes such that , resp. is vertically connected to , resp. ; in other words, if are connected from below.
Polymers with external legs have been ”mayerized” at scales , but not at scale , which means that they are non--overlapping: , are free to overlap, but if .
The final outcome of the set of scale expansions for a given scale can be found in §2.3.4, see (2.107).
2.3.2 Single-scale bounds
In this preliminary version of our final bounds, we present a bound for a single-scale amputated diagram. It is correct, but not fully conclusive, because it fails to take into account bubbles in an effective way, which will be done only in §2.3.3; as we shall see, the core principles are left unchanged by bubble resummations, but power-counting estimates are modified in an essential way. Fix , and assume we want to bound some connected, amputated -point function
[TABLE]
a sum of connected diagrams. In this subsection, we do not want to distinguish between ’s and ’s, and shift consequently to a neutral notation, . Scale cluster expansions do not produce directly , but a sum over scale polymers with external fields located in , . Polymers are in principle evaluated by expanding as a sum over Feynman diagrams; however, diagrams differing by merely exchanging vertices belonging to the same box should be evaluated together: their sum is much smaller than expected due to sign compensations (see arguments in §5.3 and later explanations). In coherence with the logic of multi-scale cluster expansions, we are really interested in power-counting estimates for the integrated quantities
[TABLE]
The dimensional factors , are chosen in such a way that is a non-dimensional scaling factor, so that has the same dimensionality as that of an amputated -point function. The fixed vertex location will be integrated later on at some lower scale (equal to the scale of one of the low-momentum external fields attached at ).
Single-scale bounds for a given diagram with fixed sectors. Let us first bound the contribution to of a single diagram , with angular sectors of all internal lines fixed. Because of the sign issue mentioned above (easily solved by summing over permutations of vertices located in the same box), but mainly because the present power-counting principle overlooks the main problem – which consists in summing over the different choices of angular sectors –, this paragraph can only be of pedagogical value, therefore it can be skipped. Here are two examples of polymers with external legs. Black lines are cluster links. We actually selected contributions to the polymers containing only one vertex per box (or cluster),
••\Psi^{2}_{ext}$$\Psi^{2}_{ext}$$\Delta_{1}$$\Delta_{2}$$\Delta_{1}$$\Delta_{2}$$\Delta_{3}$$\Delta_{4}•••••••\Psi_{ext}$$\Psi_{ext}$$\Psi_{ext}$$\Psi_{ext}
Fig. 2.3.2. Examples of cluster trees for a diagram with four external legs.
with , . Locations of external fields are indicated by the symbol (in green). Dashed lines (in blue) indicate unpaired fields. In general there may be more than one vertex per box, hence many more unpaired lines than cluster links. Pairing unpaired fields using Wick’s theorem, we obtain Feynman diagrams, e.g.
••\Psi^{2}_{ext}$$\Psi^{2}_{ext}$$\Delta_{1}$$\Delta_{2}$$\Delta_{1}$$\Delta_{2}$$\Delta_{3}$$\Delta_{4}•••••••\Psi_{ext}$$\Psi_{ext}$$\Psi_{ext}$$\Psi_{ext}
Fig. 2.3.3. For the bubble diagram (left), (in green), , (in black), (in blue), . Similarly, for the other diagram, , , , , .
General topological properties of a Feynman graph with four-valent vertices, loops, internal lines ( internal fields) and external legs imply that
[TABLE]
These relations give a quick power-counting argument for the evaluation of an amputated diagram with fixed angular sectors. Choose some spanning tree (which may be e.g. the cluster tree in the above examples, but in general, has much more vertices and edges than ; in fact one may take ). Internal lines are split into internal lines on the spanning tree (in black) – in this very particular example, cluster links – and lines generating the loops (in blue). All together, .
Now, fixing and integrating the other ’s in their respectives boxes, one obtains an overall scaling factor
[TABLE]
which takes into account the integration of vertices over their respectives boxes, and the pre-scaling factor of each propagator , see Proposition 2.3.
Counting one sum over angular sector assignments per loop yields a total scaling factor , characteristic of a just renormalizable theory with four-valent vertices. However, our aim is to prove that diagrams are smaller than that by a factor of the type for some constant , with , providing the basis for the -expansion.
Single-scale bounds for a given polymer. We shall actually use a different power-counting argument, not based on general topological properties of a Feynman graph, but on the concept of total weight per vertex (see below). Previous considerations – see Proposition 2.5 and the bound just below for the integrated vertex – imply that the leading behavior is obtained by assuming that sectors on a given vertex satisfy up to some permutation
[TABLE]
for some . We then decide that are on the same fermion loop, and similarly for . Generically, , , so it makes sense (at least graphically) to assume that the angular sector along a fermion loop is constant. Following loops, one obtains e.g.
••\alpha_{1}$$\alpha_{1}$$\beta_{1}$$\beta_{1}$$\beta_{2}$$\beta_{2}$$\Upsilon_{3}••\alpha_{1}$$\alpha_{2}$$\alpha_{1}$$\alpha_{2}$$\alpha_{1}$$\alpha_{2}$$\tilde{\Upsilon}_{3}
for the bubble diagram, and
•••••••\Psi_{ext}$$\Psi_{ext}$$\Psi_{ext}$$\Psi_{ext}$$\alpha_{3}$$\alpha_{3}$$\alpha_{3}$$\alpha_{4}$$\alpha_{4}$$\alpha_{4}$$\alpha_{1}$$\alpha_{1}$$\alpha_{2}$$\alpha_{2}$$\alpha_{2}$$\alpha_{2}$$\alpha_{1}$$\alpha_{1}$$\Upsilon_{6}
Fig. 2.3.4. Loop decompositions. Loop lengths are (), (), , , , .
for the other diagram. The diagram is the usual Cooper bubble diagram. Define the length of a loop to be the number of internal propagators belonging to it, see Figures above. Considering the generic case, we see that there are e.g. only independent sector assignments for . In , the external sectors keep outside the diagram; in all other cases, there are broken loops, finishing with two external legs, for instance, the loop with sector in .
Note that, by splitting vertices into half vertices, one has rewritten the original bubble diagram into the sum of two contributions: a bubble, , and another diagram, , which is of a different nature. Bubbles, in this sense, are characterized that they make up a loop of length 2. Loops of length 1, such as in or in , are called incomplete bubbles because external propagators may contract at some lower scale to form a full, multi-scale bubble, implying the possibility of obtaining chains of multi-scale bubbles of arbitrary lengths,
••••••
Fig. 2.3.5. A chain of 3 incomplete upper scale bubbles, possibly completed at lower scales into a chain of multi-scale bubbles.
Let us start with the argument. First of all, single scale polymers must be prepared as in [45], Proposition 5.3, by displacing all fields lying in a scale box or sub-box to the same location; see Appendix, §5.3 for a description of the procedure. Because of the fermionic nature of the fields , there remains only contributions with at most one field or per box, and terms featuring gradient fields , integrated over sub-boxes. The net outcome is a sum over Feynman diagrams with supplementary gradients, inducing an overall scaling factor bounded by
[TABLE]
for every , where , resp. is the number of , resp. -fields in the scale box .
Choose now a set of (not necessarily distinct) boxes such that , and some loop structure , namely, an equivalence set of Feynman diagrams with fixed topological structure – allowing the permutation of vertices located in the same box – and sector assignment compatible with the above set of loops, i.e. satisfying (2.87). We shall compute to begin with an upper bound to the power-counting of the sum of all Feynman diagrams in the equivalence set , including the integration of in its box, the power-counting of external legs, and the sum over external sectors. The main task consists in summing over angular sector attributions of each vertex. In order not to lose the benefit of (2.88), we fix a set of integers , and restrict to diagrams such that each list is some permutation of the list . The power-counting proceeds as follows:
(i) First select a ”loop spanning tree”, i.e. a tree connecting the loops. This may be done inductively by choosing for a root; then connecting successively to all loops such that there is a vertex connecting two successive moments along to two successive moments along , ; then connecting successively to all loops other than those previously chosen such that there is a vertex connecting two successive moments along to two successive moments along , ; and so on, until all loops have been exhausted. For instance, choosing to be the loop with sector in , the loop spanning tree associated to is •1$$4$$3•••
(ii) Reorder loops by using depth-first search algorithm on . In the case of , the algorithm visits successively vertices and , so loops are reordered by permuting indices and .
(iii) Order the vertices , along in such a way that vertices connecting to its descendant loops on the loop spanning tree (e.g. to for , with its loop reindexing defined in (ii)) are . Then choose some propagator along , denote by its sector, go around the loop in some arbitrary direction starting from , and index the successive sectors by . Write if is a vertex along connecting the momentum in the sector to the momentum in the sector , and \lfloor\log_{2}\Big{(}1+|\alpha_{i,p}-\alpha_{i,p+1}|\Big{)}\rfloor=j-k, . This means the following: there are possibilities for for fixed.
(iv) Let . The power-counting associated to a vertex with its 4 legs (see (2.13)) is (hence roughly a factor per vertex or per loop, compare with (2.86)), times , all together , that is, a -type factor, generically but not always .
(v) (overall rotation factor ) Consider one of the connecting vertices along , see (iii). Denote by , its two halves, one on , the other on . Summing over for fixed (see (2.87)) yields a factor compensating the factor in (iv). Now local rotation invariance has almost been fixed: choose some sector along ( possibilities), then the factors associated to connecting vertices fix one sector on every other loop.
(vi) (total weight associated to a connecting vertex) Let , be a connecting vertex. Given (iv), (v), there remains only a small factor , and a sum over all possible values of in , all together a factor if . If then is small, but this is actually very deceptive because of bubble resummations (see §2.3.3). The power-counting for connecting vertices associated to bubbles is better than that and computed in (viii) below.
(vii) (total weight associated to a non-connecting vertex) Let , be a non-connecting vertex. Given (iv), its total weight is , featuring a -prefactor.
(viii) (subcase of (vi): bubbles) A loop is called a bubble when , . As explained above, it can be complete or incomplete, depending on whether or . Then the outgoing transfer momentum is equal to the ingoing transfer momentum, so that in (vi) no sum over is required. Hence the total weight associated to the connecting vertex is instead of . Call strongly connecting vertex a connecting vertex which is not originated from a bubble as just described. Note that a vertex which is not a strongly connecting vertex (see cases (vii) and (viii)) has a small weight attached to it. On the other hand, strongly connecting vertices have a weight .
Assume for simplicity that . Let be the number of strongly connecting vertices, and . The total weight obtained by multiplying all the above factors is . Let be the number of bubbles (in the sense defined in (viii)). Using the equality (number of loops) – which is modified into in the case of diagrams containing only bubbles, e.g. – and the bound (expressing the fact that each loop different from a bubble has at least three half-vertices attached to it), one obtains
[TABLE]
or in the case of diagrams containing only bubbles. For instance, and , resp. for , resp. ; and and (depending on how the bubble tree is generated) or for .
The lowest possible value of for small diagrams with which are not only made up of bubbles is . Namely, such diagrams have at least two loops and two external vertices. One (connecting) vertex is needed to connect and , while another, non-connecting one is needed to assure that these do not make a chain of two bubbles. At least two other non-connecting vertices are needed to connect either to other loops or to external legs. Therefore (disregarding diagrams containing only bubbles), one may replace the exponent by , so that, improving on (2.89),
[TABLE]
Individual pre-factors , which can be interpreted as the weight of bare vertices, compare with (2.23), are modified in multi-scale bounds because effective vertices are scale-dependent, but the exponent remains the same, and is essential for polymer bounds.
One also gets a product of ””-factors, per non-connecting vertex, see (vii), which however is simply for diagrams consisting only of bubbles (chains of bubbles); this feature will however be exploited only later on, once chains of bubbles will have been resummed, see §2.3.3. The overall scaling factor must however be corrected: first, the location of is fixed, leading to a correcting factor . Second, the above weight mistakingly includes the power-counting of external fields together with their sectors, see (2.14), which leads to a correcting factor . Finally, assuming the underlying cluster tree is , the above may be multiplied by the decay factor C_{p}^{n}\prod_{\ell\in L({\mathbb{T}})}\Big{(}1+2^{-j}\frac{|\xi_{\ell}-\xi^{\prime}_{\ell}|}{\mu}\Big{)}^{-p} for every .
All together, we have proved the following bound for a single-scale diagram not made up only of bubbles, with underlying cluster tree and angular sector distribution :
[TABLE]
where are the centers of the boxes . Diagrams made up only of bubbles feature a small factor instead.
Sum over trees. This is a standard argument in statistical mechanics and constructive field theory; see e.g. [54], Corollary 5.3. Recall that there is no sum over permutations of vertices. So all trees may be generated by the following algorithm. (All boxes in the argument are in ). Fix a box . Sum at step 1 over all possible boxes (including ), and add a link between and , i.e. a pairing or between fields located at . The corresponding multiplicative spatial decrease factor is C:=\sum_{\Delta^{\prime}_{1}}\Big{(}1+d^{j}(\Delta_{1},\Delta^{\prime}_{1})\Big{)}^{-p}, where
[TABLE]
is a scaled distance between two scale boxes. Ordering boxes by their distance to , one obtains
[TABLE]
provided . Continue by picking a second link ) between and a box , and so on, until all pairings between fields in and fields either in or in any other box have been exhausted. There are local factorials involved, since fields located inside may be permuted, but inverse local factorials beat them for ; a precise argument may be found in [54], Corollary (5.3) 2., where inverse local factorials are deduced from the quasi-exponential decay of propagators, and not from the fermionic nature of the theory (see discussion in §2.4). Then, at step 2, one looks for all possible pairings between a field located in defined as , where , and a field located in , and so on. Once all vertices of the tree have been explored, the procedure stops. In the end, one finds for the sum over all scale polymers with fixed number of vertices:
[TABLE]
which is the general term of a converging series for small enough. Hence, finally,
[TABLE]
2.3.3 Bubble resummations
Let . We shall now detail the resummation of chains of bubbles with lowest fermionic scale , also called ””-scale expansion. The notation ” for a fermionic scale seemingly contradicts the principles laid out in the ”Important notations” section of the Introduction. However, as emphasized in §1.4, see (1.58), the kernel has an effective infra-red cut-off for . Therefore it makes sense in this context to identify the lowest fermionic scale with a bosonic cut-off scale.
This step proceeds ”by inspection”, namely, it requires no supplementary expansion, rather an explicit resummation by hand of structures found in the sum of perturbative diagrams.
The structures we want to single out are chains of bubbles made up of pairs of propagators with lowest scale , see Fig. 1.4.5. The lowest-order term, , connecting the two halves of any vertex,
\bar{\psi}^{j_{+}}$$\psi^{{j_{+}}\to}$$\bar{\psi}^{{j_{+}}\to}$$\psi^{{j_{+}}\to}$$\xi$$\lambda$$\equiv$$\bar{\psi}^{j_{+}}$$\psi^{{j_{+}}\to}$$\bar{\psi}^{{j_{+}}\to}$$\psi^{{j_{+}}\to}$$\xi$$\xi^{\prime}$$\lambda\delta(\xi-\xi^{\prime})
Fig. 2.3.6. Two equivalent ways a representing a local vertex.
is complemented by a sum of terms involving bubbles forming a ladder diagram,
\xi$$\bar{\psi}^{j_{+}\to}$$\psi^{j_{+}\to}$$\bar{\psi}^{j_{+}\to}$$\psi^{j_{+}\to}$$\cdots$$\cdots$$\psi^{j_{+}\to}$$\psi^{j_{+}\to}$$\bar{\psi}^{j_{+}\to}$$\bar{\psi}^{j_{+}\to}$$\psi^{j_{+}\to}$$\psi^{j_{+}\to}$$\bar{\psi}^{j_{+}\to}$$\bar{\psi}^{j_{+}\to}$$\xi^{\prime}
Fig. 2.3.7. Bubble resummation. At least one of the fields or is assumed to be of scale .
Resumming the above series yields – up to the discrepancy between the current value and the asymptotic infra-red value of the energy gap – the kernel Pre-Pre. The error term due to the difference will be bounded only in section 3, following the arguments of §3.4 C., 2. and 3. Now internal vertices of the chain may in turn be replaced by chains of bubbles with lowest momentum of scale , yielding all together Pre.
This has the following consequence; see the preliminary bounds of §2.3.2 for notations. All loops have at least vertices along them (), so case (viii) is absent. On the other hand, half-vertices are now possibly connected by Pre-kernels with , which is materialized by the substitution of delocalized vertices to the original, local vertices,
\xi$$\xi^{\prime}$$\lambda\delta(\xi-\xi^{\prime})Pre\Sigma^{j_{+}\to}(\xi-\xi^{\prime})$$\xi$$\xi^{\prime}
Fig. 2.3.8. Local vertices (on the left), substituted by delocalized vertices (on the right).
Let be a scale box, and its volume. *Assume first that . * Then the integrated vertex is , see (1.92), a bound substituting itself to the bound obtained by considering only the lowest-order term in in the Goldstone boson propagator, . For the highest scales, say , , the integrated vertex is since . However, for larger scale indices , that is, close to the transition scale, the bound is only . For , the bounds of §1.6 yield instead a large bound . Recapitulating, one may define an equivalent scale coupling constant replacing the bare constant :
Definition 2.9** (scale coupling constant )**
Let .
Previous computations show that
[TABLE]
Thus the weight of a vertex other than strongly connecting including its four half-propagators (called effective vertex in the literature), compare with §2.3.2 (viii), is (), where one has set
Definition 2.10** (scale effective vertex)**
Let
[TABLE]
and
[TABLE]
Strongly connecting vertices, see §2.3.2 (vi) and (viii), have an extra logarithmic prefactor .
Above the transition scale, i.e. if , grows logarithmically in terms of from the UV cut-off scale to the transition scale . Below transition scale, i.e. if , is the inverse of a spring factor between scales and . This shows clearly why other arguments (namely, Ward identities, developed in section 3) are required to bound polymers with . However, keeping , in line with (0.35), ensures that large factors , , bounded by , may be absorbed by putting aside some power of , say,
[TABLE]
for all provided . For definiteness, we choose
Definition 2.11** (choice of )**
[TABLE]
Hence the total weight of for a diagram with vertices is . For single-scale diagrams, or multi-scale diagrams with no external legs in Cooper pair configuration, later to form a pre-Goldstone boson, only small weights show up. The problem of potentially large weights is therefore raised up again only in §2.4.
Summarizing, and referring to the details of the computations in §2.3.2, the small factor in (2.91) undergoes the following transformations:
(vi)(vi)’ (total weight associated to a connecting vertex) , featuring now a (probably spurious but not very disturbing) logarithmic factor ;
(vii)(vii)’ (total weight associated to a non-connecting vertex)
[TABLE]
featuring a ”-type factor .
Most importantly, we note here that – after bubble chains have been resummed – connecting or strongly connecting vertices are now the same. Hence the number of factors in (vi) is n^{\prime}=n-n^{\prime\prime}\leq 2n/3\leq 2\Big{(}2+\max(1,\frac{n}{3}-1\Big{)} (since ). Combining these with the factors coming from (vii), we finally get instead of (2.91)
[TABLE]
By the same argument as in §2.3.2, one may sum over diagrams. This may be done in two stages; first one sums over all diagrams spanning a given polymer of scale , yielding
[TABLE]
where is a cluster tree spanning , and is the number of boxes of . Then the sum over polymers may be bounded using space-time decay, so that
[TABLE]
and finally, summing over all polymers of scale with one external leg located at
[TABLE]
since (considering the worst case, for which external momenta are grouped 3 by 3, see Figure below) .
Fig. 2.3.9. One example of external momentum configuration for a polymer with . The minimal number of vertices is 2.
2.3.4 Momentum-decoupling expansion and displacement of external legs
One must now test how step polymers are connected from below to boxes in . To this end, we apply to each polymer in the following expansion operator
[TABLE]
Expanding the product yields a sum of terms with , , ranging in , coding for the integral remainder term. Let , and denote by the boxes in containing the external fields. If , then all are , so is equal to the number of external fields. If then (i) either all are and the same conclusion holds; or (ii) some is equal to , which means that : when this happens, we decide somewhat arbitrarily that the polymer has exactly vertical links connecting to the box below it in . As we shall see, polymers with will undergo a further treatment.
The final outcome of the horizontal and momentum-decoupling expansion may be written as
[TABLE]
where:
– ranges over the set of scale cluster forests;
– are the non-overlapping polymers introduced at the end of the last subsection;
– and is as in (2.82).
Instead of summing over cluster forests, it is equivalent, but more convenient, to sum over non--overlapping polymers, so that
[TABLE]
in (2.107), where , are cluster trees compatible with the choice of non--overlapping polymers , namely: are the connected components of , and . It is in this form that we shall be using (2.107). Thus:
[TABLE]
Displacement of external legs. Let us discuss the smallest values for :
(i) if , then has no external field, so is a vacuum polymer;
(ii) if , then is a two-point polymer, with corresponding contribution
[TABLE]
where
[TABLE]
see (2.56), is the local part of the polymer, obtained by displacing symmetrically one of the external fields to the location of the other. The error term , once its external legs are contracted (at some later stage, namely, at scales and ), may be rewritten as a Taylor remainder as in (2.56), and comes with an extra spring factor . Since two-point functions are linearly divergent, subtracting local part to order 1 makes the remainder convergent. Referring to the logarithmically divergent, four-point bubble diagram, and to the power-counting of a contracted pair , may be thought of as having external legs; however, for pure book-keeping reasons, we shall consider it as a polymer with external legs – which simply means that it is convergent enough for our purposes. For future use we indeed replace with a polymer which is identical to , except that it has external legs in the box where the external fields have been displaced to, , resp. , and none in the other box, , resp. .
(iii) if , then we proceed as in (2.111), Taylor expanding to order 2 symmetrically in the neighborhood of each external vertex location . All such diagrams – recall Cooper pair bubble diagrams have already been resummed by hand – are already convergent before displacing external legs. However (see §2.4), the above operation will help us sum over external angular sectors in the multi-scale bounds. Local contributions are translation invariant polymers with external legs; the error term , on the other hand, is considered as as a polymer with external legs.
2.3.5 Mayer expansion
We shall now apply the restricted cluster expansion, see Proposition 5.4, to the result of our expansion. Cluster expansions have allowed us to rewrite Green functions as sums over non--overlapping polymers denoted by . The objects are now polymers in ; a link is a pair of polymers , . Objects of type 2 are polymers with external legs, whose non-overlap conditions we shall not remove. Then objects of type are polymers with external legs. Due to the displacement of external legs operated in §2.3.4, all external legs are located in the same scale box.
Implicit in the outcome of the cluster expansions is the non-overlapping condition,
[TABLE]
stating that a scale box belonging to and a scale box belonging to are necessarily distinct. Similarly to what we did during the horizontal cluster expansion, we choose some polymer with external legs, say , and weaken the non-overlap condition between and all the other polymers by introducing a parameter ,
[TABLE]
where is the subset of boxes with external legs - i.e. that have been differentiated with respect to -, and Taylor expand in to order ; each factor
[TABLE]
produced by differentiation is a Mayer link between and some , or more precisely, some box and some box , implying an explicit overlap between and , and adding a link to the forest . Iterating the procedure and applying Proposition 5.4 to the weakened non-overlap condition
[TABLE]
The outcome is a sum
[TABLE]
Links are obtained as links between polymers, however the corresponding differentiation is immediately rewritten as a sum over pairs of boxes . Thus we see Mayer links as links between boxes. As such they add up to the set of links produced by the horizontal cluster expansion, producing a forest with same vertices as but larger set of links , where (in bijection with ) is the set of Mayer links. Since a forest is characterized by its set of links, we rewrite in practice (2.117) as
[TABLE]
The number of external legs of a set of polymers connected by Mayer links is the sum of the number of external legs of each of the polymers. In particular,
- (i)
new vacuum polymers without any non-overlap conditions have been produced; they are resummed into the scale free energy function ;
- (ii)
two-point polymers have been dressed by a cloud of vacuum polymers; they can now be resummed into a renormalization of the two-point function;
- (iii)
links from polymers of type 1 to polymers of type 2 produce new polymers with external legs, whose overlap conditions have not been removed.
Other possibilities include e.g. polymers with external legs without any non-overlap condition, produced by Mayer-linking two-point polymers.
Let us now give some necessary precisions. Since the Mayer expansion is really applied to the non-overlap function NonOverlap and not to the outcome of the expansion, one must still extend the outcome of the expansion to the case when the , have some overlap. The natural way to do this is to assume that the fields (\Psi^{j,\alpha}\big{|}_{\mathbb{P}_{i}},\bar{\Psi}^{j,\alpha}\big{|}_{\mathbb{P}_{i}})_{i=1,\ldots,n} remain independent even when they overlap. This may be understood in the following way. Choose a different color for each polymer , and paint with that color all boxes . If , then its external vertical links to are left in black. The previous discussion implies that boxes with different colors may superpose; on the other hand, external inclusion links may not, so that low-momentum fields do not superpose and may be left in black.
Hence one must see as living on a two-dimensional set, , so that copies of with different colors are independent of each other. This defines new, extended fields , and Mayer-extended polymers. By abuse of notation, we shall skip the tilde in the sequel, and always implicitly extend the fields and the measures by taking into account colors.
As a general principle (see [9] for a single-scale version, or [54] for the general, multi-scale version), sums over Mayer-extended polymers are bounded exactly as sums over non-Mayer-extended polymers, , see (2.103), (2.123), with the constant of non-Mayer-extended polymers replaced by .
2.4 Multi-scale bounds and fermionic fixed-point
Our starting point is the single-scale bound
[TABLE]
see (2.103), complemented by multi-scale bounds for the effective vertex, (2.104) and (2.105), and Definition 2.10, which imply a total supplementary weight in a multi-scale setting. Recall . We shall now prove by scale induction the following multi-scale bound. Fix and consider a polymer with external legs , , . Let be the highest external momentum scale, considered as the external scale of . The procedure considered in §2.3.4 has produced three types of polymers:
- (i)
local parts of polymers with external legs;
- (ii)
polymers with external legs; and
- (ii’)
error terms associated with polymers with external legs, which have a spring factor , and have been considered as polymers with external legs in the above book-keeping.
As we shall see, the power-counting associated to external leg angular sectors is more favorable for translation-invariant polymers (case (i)), for which Proposition 2.5 holds. On the other hand, polymer contributions of type (ii), (ii’) enjoy (though for different reasons) a spring factor , making it possible to control them in a similar way.
Note that, for every ,
[TABLE]
and (by (2.100)) , hence
[TABLE]
For multi-scale polymers of lowest internal scale and highest external scale , the factor in (2.120) is multiplied by the spring factor . For , the sum
[TABLE]
is only , thus one expects the decay exponent to be lowered to at most. With a decay however, the sum (2.121) yields but with a logarithmic correction. Hence it is safer to say that a lower decay exponent holds, say, only, compare with the exponent in (2.120).
Our claims are therefore the following. (2.103), (2.104), (2.105) and (2.120) ensure that they hold for a single-scale polymer. Let and
[TABLE]
General fermionic multi-scale bound. Let not made up only of bubbles, then
[TABLE]
Furthermore, if ranges in the set of polymers in with vertices,
[TABLE]
from which, finally, summing over all polymers in with one external vertex located in ,
[TABLE]
compare with the discussion above Fig. 2.3.9. The single remaining factor , coming originally from the -expansion, is sufficient to perform the sum over scales, as we shall see.
We now consider polymers with vertices, forming at scale a polymer through pairings of (some of) their external legs between themselves or with scale propagators expanded at scale .
j$$k_{1}$$k_{2}$$k_{3}$$\mathbb{P}^{k_{1}\to}$$\mathbb{P}^{k_{2}\to}$$\mathbb{P}^{k_{3}\to}
Fig. 2.4.1. Multi-scale polymer with lowest scale and , obtained by contracting , , with , resp. , external legs. The blob represents a polymer expanded at scale and connecting the scale external legs of and .
Assuming (2.123) holds for the connected components of , we want to show that it also holds for . For that, we must take into account several factors related to the scale integrations. The following discussion is relative to a polymer with external scale , which is one of the ’s. By assumption . For simplicity, we remove all dimensional constants () from our estimates.
Volume factors. External legs of which are contracted at scale are integrated in a box of scale . This implies an absolute volume factor .
Field scaling factors. Each external field , of which is contracted at scale comes with a supplementary absolute prefactor .
Taking into account volume and field scaling factors yields the general bound
[TABLE]
where is a bound on the sum over all possible sector assignments of external legs of polymers contracted at scale . The prefactor comes from an over-counting: the number of scale contracted external fields is not , but \Big{\{}\sum_{i}N_{ext,i}\Big{\}}-N_{ext}, which changes the overall field scaling factor; and the overall volume factor is corrected by a factor , since amputated diagrams are evaluated by definition with one vertex location fixed. Polymers with 4 external legs may be considered as four-valent vertex insertions in a scale diagram. From the diagrammatic analysis of §2.3.2 and §2.3.3, it follows in particular that a compound polymer obtained by contracting only polymers with 4 external legs has at least one extra non-connecting vertex, which produces an extra factor , see (2.120), w.r. to the bound (LABEL:eq:2.106).
Sector assignment factors. We rewrite the overall sector assignment factor as , where , and
[TABLE]
and
[TABLE]
Factors , combine nicely with the factors between brackets \Big{\{}\ \cdot\ \Big{\}} in (LABEL:eq:2.106): while \Big{\{}\prod_{i=1}^{p}2^{-j\frac{N_{ext,i}-2}{2}}\Big{\}}\times{\cal N}^{\prime}_{sec}=\prod_{i=1}^{p}\overline{\cal N}_{sec,i}. Thus
[TABLE]
The value of the correcting factors and will be presently computed. The induction is successful provided we manage to prove that
[TABLE]
if , or
[TABLE]
in the specific case when for all . Sums take care of all possible choices for scales .
Scale momentum conservation for local parts of polymers with . Local parts may be seen as sums of Fourier diagrams with external momenta , all coming from the same external vertex. Consider a local part whose external momenta are all contracted at scale . Let be their angular sectors. Since (by momentum conservation) , from which \Big{|}\sum_{i=1}^{N_{ext}}p^{j,\alpha_{i}}\Big{|}\lesssim 2^{-j}p^{*}_{F}, and sectors of contracted legs are shared between the two contracted fields, Proposition 2.5 gives an external sector assignment . Thus one sets for such polymers. Then
[TABLE]
Polymers with external legs. For polymers of type (ii) or (ii’), momentum conservation does not hold in general, so that (disregarding the sharing issue) the number of scale external sector assignments is . The factor in the term between brackets \big{\{}\ \cdot\ \big{\}} in (2.130) can however be put to good use, enhancing the factor to , with the same exponent as for local parts of polymers.
Next, each possible scale sector assignment for each of the external legs must be redivided into scale subsectors. Assume once again that all external momenta are contracted at scale . Compared to the external sector assignment factor for local parts (see previous paragraph), the outcome is that polymers of type (ii) have an extra (shared) relative factor . For polymers of type (ii’), this factor is compensated by the spring factor , so one may set as for polymers of type (i). So let us concentrate on the case of polymers of type (ii), for which is precisely this extra relative factor. The factor is then compensated by the factor found in the terms between brackets \big{\{}\ \cdot\ \big{\}} in (LABEL:eq:2.106) or (2.130): namely,
[TABLE]
if . The relative factor makes it possible to sum over scales . The extra remaining factor, per polymer with external legs, may be used to produce the required overall factor in (2.123), except, of course, if all polymers are local parts.
Scale external legs. The above arguments double-count the external sector assignment for legs which are not contracted at scale , and therefore remain external legs of the polymer . Removing this double-counting produces an extra sector assignment factor when scale momentum conservation does not hold globally, i.e. when there is at least one polymer with external legs. On the other hand, if all polymers are local parts, then momentum conservation for each of them implies momentum conservation for , hence the extra sector assignment factor is . Thus in the first case, in the second case.
Case when all polymers are local parts. In that case, all correcting factors are equal to 1, so (2.131) holds when all . Otherwise, assuming e.g. that , one has , whence (2.130) holds.
Combinatorial factors. Compared to single-scale trees, multi-scale trees involve supplementary combinatorial factors. This is a well-known problem in constructive field theory, which has been e.g. exposed in details and some generality for bosonic theories in [54]. The fact that the present model is fermionic, hence obeys Pauli’s principle, leads to substantial simplifications.
Accumulation of low-momentum fields (see [54], Remark below Definition 2.9). Let , , and , , be the boxes of given scale included in . The scale -derivatives may have produced up to low-momentum fields or in each box . Decomposing into scales, this means that there are for each sector up to fermion fields located in a single box . The inverse local factorial factor of (2.105) makes it possible to sum over all pairing possibilities.
3 Low-energy theory for Cooper pairs
The picture emerging after the expansions of section 2 is that of fermionic polymers connected by chains of bubbles evaluated as Pre-kernels. The infra-red cut-off at scale keeps away infra-red singularities, and the expansions ensures that the series of perturbations converges. This is still clearly one step behind the conjectured, divergent infra-red behavior in of the Goldstone boson propagator (or kernel), found in principle by means of the gap equation as in §1.7.
We need in the first place a proper definition of the Goldstone boson propagator that goes beyond perturbation theory. Instead of developing blindly an infinite sum of individual four-point diagrams, and selecting those which are two-particle irreducible by checking all possible ways of cutting it into two pieces, we implement a decision rule relying only on the number of links explicitly constructed by hand between pairs of polymers seen as small, almost pointwise vertex insertions. This is the object of the complementary cluster expansions of §3.2. Those are horizontal cluster expansions, between scale fermions only, obtained by expanding links to third order between boxes of pairs of polymers. The idea is that we exhaust explicitly in this way all possible configurations of the type
\Delta_{1}$$\Delta_{2}$$\Delta_{3}$$\Delta_{4}polymerpolymer
where two polymers are connected by exactly two links; displacing one of the blue fermion links so as to form a bubble as in Fig. 1.7.5, one obtains eventually a rung connecting Bethe-Salpeter kernels. Such configurations contribute to the infra-red divergence of the Goldstone boson propagator, therefore the two polymers are not sufficiently connected to be able to sum over the location of one of them w.r. to that of the other. On the other hand, three links suffice to ensure strong connectivity between these. Building upon this concept of strong connectivity, we develop multi-scale trees whose vertices are polymers, and edges – three-fold links. Because polymers with external legs are not diverging, emphasis is laid on four-point polymers. Those are considered to have a bosonic scale given by the scale of their transfer momentum , which may be determined to some precision by refining the sector decomposition of the fermions.
The outcome of these complementary expansions is a expression of the kernel for in terms of a sum over polymer configurations, -point functions are then rewritten in terms of sums over simplified trees called bosonic trees, whose vertices are still polymers, but edges are now single-scale Goldstone boson propagators (not yet Goldstone boson propagators, see below) of given scales, see Fig. 3.2.14.
However, even at this stage, we are not able to prove that -point functions are infra-red summable, because the theory is a priori plagued with infra-red divergences coming from Goldstone boson propagators with low transfer momentum, connecting arbitrarily far-away fermionic polymers. This problem is solved in §3.3 using Ward identities similar to those used to solve the problem of ”soft photons” (infra-red photons) in QED. First of all, multi-scale two-point functions of Cooper pairs, inductively extracted from bosonic trees (see Fig. 3.2.14), have no associated Ward identity, hence must be resummed by hand. This is a simple geometric series, who sum gives at last the Goldstone boson propagator we were looking for. As a consequence of these partial resummations, previous bosonic trees are replaced by very similar bosonic trees with the same type of vertices, but whose edges are Goldstone boson propagators, and that furthermore have no internal two-point functions. We may now solve the gap equation by an easy fixed-point argument, and verify the predictions of §1.7 concerning the infra-red behavior of the geometric series .
In the large-scale limit, i.e. at space-time scales much larger than , polymers of bosonic trees appear in our multi-scale picture as tiny islands giving subleading contributions to -point functions. This makes plausible the general structure of the formula of Theorem 2,
[TABLE]
As apparent in (3.1), the effective theory is roughly that of a free bosonic particle – the Goldstone boson – with covariance kernel . This Theorem is proved in §3.5.
Let us now turn to another question. The effective potential approach described in the Introduction (Regime II) suggested that the field conjugate to Cooper pairs should behave like a non-linear sigma-model valued in the immediate neighborhood of the circle of radius . Thus, it is in principle possible to distinguish transverse (i.e. orthogonal to the circle) fluctuations – ”slow” degrees of freedom – from tangential fluctuations – ”fast” degrees of freedom.
Here, by comparison, a preferential direction has been fixed by the symmetry-breaking term. Thus all there remains to be done is to check that with high probability, the Cooper-pair field is to a good approximation aligned with this direction, and has a modulus fixed by the gap equation. This is the content of Theorems 4 and 5 in the Introduction, proved resp. in §3.7 and §3.8.
The introductory, heuristic subsection §3.1, recapitulates the main findings of mean-field theory concerning the infra-red behavior of the theory. Classical results presented there without proof are only half-quantitative since the Goldstone boson propagator has been replaced by its leading order obtained by resumming Cooper pair bubbles, but notations introduced there are important for the subsection on Ward identities.
3.1 A heuristic introduction
We present in this subsection a heuristic derivation of the large-scale behavior of the theory in the Cooper channel, based on elementary considerations, and the fundamental assumption that asymptotics may be obtained by considering only chains of Cooper pair bubble diagrams. Notations introduced here are fundamental for §3.3 where Ward identities will be proved, justifying the present heuristics beyond its pedagogical value.
Let us first write down the various quadratic fields built out of one and one :
• Cooper pairs:
[TABLE]
dividing by , we get a different possible normalization, \bar{\Psi}\sigma(\theta)\Psi=\bar{\Psi}\left(\begin{array}[]{cc}0&e^{-{\rm i}\theta}\\ e^{{\rm i}\theta}&0\end{array}\right)\Psi, where
[TABLE]
• normal pairs:
[TABLE]
• identity pairings:
[TABLE]
Let us now translate expressions (1.78,1.79) for in a new base: denoting by the direction parallel to , and the direction rotated by an angle , so that
[TABLE]
we consider the Pre-kernel in the basis:
[TABLE]
When , one gets , and
[TABLE]
and similarly
[TABLE]
[TABLE]
Using rotation equivariance, it is easy to prove that
[TABLE]
[TABLE]
Hermitian symmetry implies that is diagonal in the -basis. Now, (1.79,1.78) implies that
[TABLE]
The kernel has a characteristic massive decay for , with mass . Hence the main term in the large-scale behavior of -point functions of Cooper pairs is obtained by considering only bound states propagating through space-time through chains of Cooper space bubble diagrams
\sum_{k\geq 0}$$\sigma(\theta+\frac{\pi}{2})$$\sigma(\theta+\frac{\pi}{2})$$\sigma(\theta+\frac{\pi}{2})$$\sigma(\theta+\frac{\pi}{2})$$\sigma(\theta+\frac{\pi}{2})$$..........$$\sigma(\theta+\frac{\pi}{2})$$\sigma(\theta+\frac{\pi}{2})$$1$$2$$k$$\sigma(\theta+\frac{\pi}{2})
3.1**.1. Bound state propagation.**
which is simply the Pre-kernel
in the perpendicular direction.
If , the leading-order term of the function
\Big{\langle}\ :\,(\bar{\Psi}{\mathbb{\Gamma}}(\theta^{\prime})\Psi)(\xi)\,:\ \ :\,(\bar{\Psi}{\mathbb{\Gamma}}(\theta^{\prime\prime})\Psi)(\xi^{\prime})\,:\ \Big{\rangle}_{\theta;\lambda} is (letting )
\Gamma(\theta^{\prime})$$\Gamma(\theta^{\prime\prime})$$+\sum_{k\geq 0}$$\Gamma(\theta^{\prime})$$\sigma^{\perp}$$\sigma^{\perp}$$\sigma^{\perp}$$\sigma^{\perp}$$\sigma^{\perp}$$..........$$\sigma^{\perp}$$\sigma^{\perp}$$1$$2$$k$$\sigma^{\perp}$$\Gamma(\theta^{\prime\prime})
Fig. 3.1.2. Two-point function.
evaluated as the inverse Fourier transform of
[TABLE]
where one has neglected the difference in the term without kernel. Since are not singular in the infra-red limit, this gives to leading order in if
[TABLE]
Graphically,
\Big{\langle}\ :\,(\bar{\Psi}{\mathbb{\Gamma}}(\theta^{\prime})\Psi)(\xi)\,:\ \ :\,(\bar{\Psi}{\mathbb{\Gamma}}(\theta^{\prime\prime})\Psi)(\xi^{\prime})\,:\ \Big{\rangle}_{\theta;\lambda}\approx••\xi$$\xi^{\prime}$$\frac{\Gamma}{\lambda}$$\frac{\Gamma}{\lambda}
Generalizing, one obtains a semi-perturbative version of Theorem 2: the leading term of \Big{\langle}\prod_{i=1}^{2n}\,:\,(\bar{\Psi}{\mathbb{\Gamma}}^{\perp}\Psi)(\xi_{i})\,:\ \ \Big{\rangle}_{\theta;\lambda} for widely separated points is graphically (taking e.g. )
\frac{\Gamma}{\lambda}$$\frac{\Gamma}{\lambda}$$\frac{\Gamma}{\lambda}$$\frac{\Gamma}{\lambda}$$\frac{\Gamma}{\lambda}$$\frac{\Gamma}{\lambda}••••••\xi_{1}$$\xi_{4}$$\xi_{2}$$\xi_{5}$$\xi_{3}$$\xi_{6}+ perm.
Fig. 3.1.2. -point function of Cooper pairs.
where ”perm.” indicates the sum over all possible pairings.
Let us now consider the one-point function in the parallel direction . Main term is
\langle(\bar{\Psi}\sigma(\theta)\Psi)(\xi)\rangle_{\theta;\lambda}=$${\mathrm{Tr}}((\sigma(\theta))^{2})•
where the coefficient 2 comes from the trace . The gap equation, see Fig. 1.7.6, implies that the correction due to the Goldstone boson
•\xi$$\sigma(\theta)$$\sigma(\theta)$$\sigma(\theta)$$\sigma(\theta)$$\sigma(\theta)$$..........$$\sigma(\theta)$$\sigma(\theta)$$\sigma(\theta)$$-$$\Gamma$$+••
vanishes identically.
Let us finally mention another identity,
[TABLE]
obvious by the symmetry exchanging - and -states.
3.2 Complementary expansion
After these heuristics and general power-counting considerations, let us now come back to where we left the story. At the end of section 2, the situation was as follows: -point functions have been rewritten as a sum over polymers with lowest scale . Because two-point functions have been renormalized at scale , and vacuum polymers have been divided out using Mayer’s expansion, there remain only polymers with external legs, and polymers with (not fully produced) external legs. Since there is only one scale left, all external legs are of scale . The measure has retained an exponential weight , where ranges over the set of scale boxes – thereafter called isolated boxes – that are not connected by a vertical link to boxes of higher scale , and
[TABLE]
compare with (2.36) and (2.66), since by hypothesis . Recall that has an extra pre-factor. As apparent from §1.7, a particular rôle will be played by what are called polymers of type 2 in the following
Definition 3.1** (polymers of type 1 and type 2)**
Let be a four-point polymer with scale external legs in angular sectors . Then is said to be of type 2 if all the following conditions are satisfied:
- (i)
;
- (ii)
.
Its external structure is either non-mixing, , or mixing, or . Furthermore, if the polymer is mixing, we require that
- (iii)
.
All other polymers in , and also isolated scale boxes, are said to be of type 1.
Polymers of type 2 are exactly those which make part of the Bethe-Salpeter kernel, accounting for the large-scale behavior of the theory. Roughly speaking: is a polymer of type 2 if it is a four-point polymer with transfer momentum such that . Note that conditions (ii), (iii) have some degree of arbitrariness. For instance, because angular sectors do not have sharp cut-offs, one may have but , or vice-versa. Also, in the case of a mixing diagram, say with external structure , one may exchange and . Condition (iii) ensures that our criterion is unambiguous: there can be only one numbering of external legs such that is of type 2. As we shall see later on, such subtleties are irrelevant in the infra-red limit.
A. General introduction. The general philosophy in this subsection is to perform a partial, carefully devised separation of fermionic degrees of freedom from bosonic degrees of freedom. In theory, this means the following. Imagine the interaction has been wholly expanded into series, after which Wick’s formula has turned an -point function into an (infinite) sum of diagrams. Then sequences
\bullet$$\bullet$$\bullet$$\bullet
see Fig. 1.7.4, of two-particle irreducible, four-point diagrams with Cooper pair external structure connected by two fermionic lines, should be resummed into a Goldstone boson propagator, plus error terms due to displacements and projections, as in §1.7. Then diagrams would be bound as in §2.3 by following fermionic loops.
Things are however not so simple, because of the well-known problem of accumulation of low-momentum fields characteristic of non-massive bosonic theories (see e.g. [54]). Let us give a sketchy description of this problem, assuming for simplicity that . Suppose that one has produced one Cooper pair with transfer momentum () per scale box included in a large scale box ; in total, fields Cooper pairs. Then the (supposedly asymptotically Gaussian) contribution of all possible -pairings of Cooper pairs through Wick’s formula would lead to a combinatorial factor , in other words, per Cooper pair. This is more than can be compensated by the weight due to the products of the -kernels, which is (leaving out dimensional constants) roughly . Taking further into account small factors due to Ward identities, see §3.3, would contribute in the best case per Cooper pair: a global counting, not sufficient though to sum over all scales .
Strategies often developed to deal with this problem in equilibrium statistical physics, relying mainly on large-deviation estimates (see [54]), are difficult to implement here because the Goldstone boson has not been introduced in the first place. Instead we choose to expand some and not all Goldstone boson propagators, in such a way as to form a multi-scale forest, where scales refer to the bosonic scales of transfer momenta. As briefly mentioned in the introduction to §2.3, we develop links between objects, which are polymers of lowest fermionic scale , with external legs of scale . Fix some bosonic scale . At that scale, we must still deal with potentially diverging four-point functions with transfer momentum of scale . In order to read easily the transfer momentum of a Cooper pair or , we decompose into angular sectors , following (2.10), , . Then the total momentum of a pair or is if and only if . Apart if (in which case we also want to deal with the remaining interaction ), we mainly want to understand how polymers of type 2 are connected at scale . Compared to the fermionic cluster expansion of §2.3, things are more involved, because expansion rules are strongly dependent on the scales of transfer momenta. Furthermore, we shall also need to resum (by inspection) two-particle irreducible four-point diagrams with Cooper pair external structure, so as to produce the Goldstone boson propagator .
Extending Definition 3.1, we may say that a four-point polymer with external structure
[TABLE]
() is a scale polymer of type 2 provided , and (if the polymer is mixing) , so that, in particular, the transfer momentum has scale . However, as we shall see, this notion is essentially redundant as soon as , because four-point polymers not satisfying these conditions have already been taken care of at some higher scale .
As a general principle, polymers will be of two types:
– scale four-point polymers (), generically called low transfer momentum four-point polymers, whose transfer momentum is ;
– and ”scale-neutral” polymers – that is, all other types of polymers –, whose external leg transfer momenta we do not need to determine, because they cannot by themselves produce divergences in the infra-red limit.
Particularizing, we have irreducible polymers,
scale-neutral4$$\geq 6$$2$$j^{\prime}_{\phi}$$\vdots$$j_{+}
Fig. 3.2.1. Irreducible polymers.
Polymers are represented by blobs.
Figures at the center of scale-neutral polymers indicate the number of external legs. Polymers with 2 scale external legs contribute to an inessential renormalization of the two-point function (whereas two-point polymers with higher momentum external legs have already been resummed in section 2). Scale-neutral four-point polymers are four-point polymers with a large transfer momentum . Six-point polymers are convergent, whence one does not need to investigate the scale of transfer momenta of external leg pairs. The same holds for lonely Goldstone pre-boson propagators coming from ”half-integer” fermionic scales, which have not become part of a fermionic polymer: they may couple to scale four-point polymers because their transfer momentum may be arbitrarily small; however, due to their scale infra-red cut-off, such propagators have quasi-exponential decay at distances larger than , hence they may be considered at lower energies as irreducible.
The dotted line at the center of a scale four-point polymer separates external legs with index from external legs with index , in such a way that the transfer momentum flowing from either side is .
Assembling irreducible polymers yields irreducible and reducible polymers. There are no scale-neutral reducible polymers; contracting scale-neutral polymers simply gives a new scale-neutral polymer, whose number of external legs is at least equal to the sum of external legs of each constituent polymer, minus twice the number of pairings (”at least” because pairings of not fully expanded polymers may yield new vertices instead of contracting previously developed external legs),
4$$\times$$\times$$\geq 6$$2$$\longrightarrow$$\geq 6$$4$$2
Fig. 3.2.2. Assembly rules of scale-neutral polymers: an example.
Scale reducible polymers are made up of scale four-point polymers and possibly scale-neutral polymers, interspersed with bubbles: pure scale polymers generalizing chains of bubbles,
+$$+\cdots
Fig. 3.2.3. Generalized chains of bubbles.
and chains including one or more contribution from scale-neutral polymers with external legs,
scale-neutral
Fig. 3.2.4. Insertion of scale-neutral polymers into chains of bubbles.
Double lines are Cooper bubbles with dressed propagators obtained by resumming scale-neutral two-point polymers. Fig. 3.2.3 may be seen as the leading term of Fig. 3.2.4, the scale-neutral polymer being just a four-point vertex with a coefficient. After a Mayer expansion and an -wave projection (see D. below for the detailed procedure), all these chains are resummed to form so-called single-scale Goldstone boson propagators with transfer momentum , denoted by a dotted wiggling line
(not to be confused with the previous Pre- kernel dashed wiggling lines). A final resummation procedure produces at last our Goldstone boson propagator , which contain internal single-scale Goldstone boson propagators with transfer momentum with higher scales .
There is a slight inaccuracy in Fig. 3.2.1 and Fig. 3.2.2, in that the scales of the low transfer momentum polymers need not be exactly the same; because momentum cut-offs are not sharp, one has instead . Note also that a Mayer resummation is necessary to have momentum conservation.
If the inserted scale-neutral polymer has exactly 4 external legs, then this configuration forms a chain. If, however, it has external legs, then its supplementary legs may also contract to an arbitrary number of low transfer momentum four-point polymers according to the following pattern (where one has chosen ),
scale-neutralj_{+,1}$$j_{+,2}$$j_{+,3}
Fig. 3.2.5 (same as Fig. 3.2.6).
Scale-neutral polymers can also appear at either end of a chain.
Since external legs of a scale-neutral polymer with external legs – which is virtually connected to an infinite number of polymers, since the scale interaction has not been expanded inside – can be paired with an arbitrary number of low transfer momentum four-point polymers with arbitrary scales, it is only at the very end of the expansion – i.e., once all low transfer momentum four-point polymers have been connected – that the general connectivity structure of emerges. Then the scale-neutral polymer of Fig. 3.2.5 will be interpreted as a vertex of the multi-scale bosonic cluster tree,
j_{+,1}$$j_{+,2}$$j_{+,3}
Fig. 3.2.6. A vertex of the multi-scale bosonic cluster tree (same as Fig. 3.2.5).
Finally (see next Figure) a low transfer momentum four-point polymer may contract one of its two Cooper pairs to two different (scale neutral or low transfer momentum) polymers (left part of the figure), in which case the three polymers may be considered as a single, compound scale-neutral polymer (blue rectangle on the right part of the figure).
p$$-p+q$$q$$j_{+}
Fig. 3.2.7. One scale-neutral Cooper pair contraction. Blobs with ascending/descending arrows are located above (), resp. below () the double line, and represent scale-neutral, resp. low transfer momentum four-point polymers.
The remaining Cooper pair may contract to a single scale four-point polymer, yielding the first bricks of a scale Goldstone boson attached to a scale-neutral polymer; otherwise it does not show up any more as a Cooper pair, and becomes part of a scale-neutral polymer.
Let us start with a precise typology of irreducible diagrams, depending on the scale:
– at scale , there are polymers in with external legs, four-point polymers of type 1, and four-point polymers of type 2. Four-point polymers of type 1 are not divergent in the infra-red limit Four-point polymers of type 2 are further split (see below) into four-point polymers of type 2 with transfer momentum of scale – whose connections are closely examined through the scale cluster expansion – and four-point polymers of type 2 with low transfer momentum of scale . The latter ones are left for further investigation at scale . As for four-point polymers of type 2 with transfer momentum of scale , they are assembled together in a line into a single-scale Goldstone boson propagtor with transfer momentum of scale . The two ends of each scale Goldstone boson propagator are either connected to one (or several) of the external legs ; or connected to diagrams with external legs, whose links with four-point polymers of type 2 with low transfer momentum have not been clarified at this stage. As mentioned earlier, scale Goldstone boson propagators may also be connected to Goldstone boson propagators of lower scale because momentum cut-offs are not sharp;
– the final diagrammatical outcome of the scale cluster expansion is:
(i) a new set of polymers in with external legs;
(ii) the left-over set of four-point polymers of type 2 with low transfer momentum;
(iii) newly developed vacuum polymers (which must be divided out by a Mayer expansion), and two-point polymers, which contribute to a small, finite dressing of fermion propagators (therefore requiring no mass or chemical potential renormalization), resulting in the replacement of by the dressed propagator ;
– at scale , as discussed in the two previous paragraphs, all scale four-point polymers of type 2 with transfer momentum have been resummed and integrated into larger polymers. Therefore there remain only four-point polymers with transfer momentum scale : all four-point polymers (save for inessential cut-off effects) are already of type 2 for the scale . After performing (as at scale ) a splitting according to whether the transfer momentum is of scale or higher, and leaving the latter case for further investigation at scale , we are left with an assembly composed of scale four-point polymers, and scale-neutral two-point polymers and polymers with external legs.
By convention, scale-neutral polymers (i.e. two-point polymers, polymers with external legs, and and scale four-point polymers of type 1) are drawn on multi-scale figures in a separate upper line (above the double blue line).
– the above description of the final diagrammatical outcome of the scale cluster expansion, and of the scale -cluster expansion, extends to lower scales without modification.
Let us first describe diagrammatically the various expansions that shall be needed. We add some short comments and refer to B. and C. below for detailed explanations.
Preliminary second-order expansion in isolated boxes (step ) at scale .
=\ 1\ +$$+•••\cdots$$\cdots$$=\ 1\ +$$+$$\geq 6•
Fig. 3.2.8. Preliminary second-order expansion in isolated boxes.
Expanding the interaction to order 2 in each isolated box yields , plus a term with one vertex – interpreted as a scale polymer with external legs, on equal footing with four-point polymers in –, plus a term with vertices – interpreted as a scale polymer with external legs –. This is done only once – as a preliminary step to scale cluster expansion.
Scale splitting phase for four-point polymers of type .
••••\xi_{1}$$\xi_{3}$$\xi_{2}$$\xi_{4}$$\simeq••\xi_{1}$$\xi_{3}$$=••\xi_{1}$$\xi_{3}$$+$$\chi_{+}^{j_{+}}••\xi_{1}$$\xi_{3}$$\chi_{+}^{\to(j_{+}+1)}
Fig. 3.2.9. Scale splitting phase for four-point polymers.
Consider a four-point polymer of type with external legs located at which have not been contracted at scales . Displace external legs as in §1.7, and . Then insert characteristic functions at singling out the case when the transfer momentum is of scale . Since the polymer is of type , the remainder (last diagram on the right) involves only low transfer momenta of scale . Low transfer momentum four-point diagrams are left for future investigation at scale .
A single-scale diagram with trivial strongly connected components. Specifically at scale , we develop links between scale isolated boxes and fermionic polymers in , except four-point polymers with low transfer momentum, . The algorithm is as follows. Choose some arbitrary ordering of the set of polymers (four-point polymers with low transfer momentum excluded), . Start from and try to develop links between and . Say that and are strongly connected if there exist such links. In the example below – where all polymers are four-point polymers with scale transfer momenta, represented by numbered blobs – no pair of polymers is strongly connected. The general algorithm is simple: expand all links between and , (first figure). Then expand all links between and (second figure), and so on, until one has obtained a new polymer which is isolated at scale . The figures below present the successive stages of the expansion. Arrows point successively from , uncovering links between and new polymers .
1$$2$$3$$1$$2$$j_{+}
1$$2$$3$$1$$2$$4$$5$$3$$4$$5$$j_{+}
1$$2$$3$$1$$2$$4$$5$$3$$4$$5$$6$$j_{+}
1$$2$$3$$1$$2$$4$$5$$3$$4$$5$$8$$6$$7$$j_{+}
Fig. 3.2.10. Dressed propagator, dressed bubble. Wiggling lines at the two ends stand for low transfer momentum Cooper pairs.
This diagram is actually a ”dressed” bubble, involving one bare propagator (the one with index ), and a dressed propagator from blob number to blob number , involving blobs and .
The above algorithm also produces non-trivial strongly connected clusters of polymers, discussed in the next two examples.
A single-scale diagram with a single, non-trivial strongly connected component.
1$$2$$3$$j_{+}
Fig. 3.2.11. One contribution to the -kernel.
Blob has links to blob . Hence blobs 1 and 2 are assembled into a single cluster . Then cluster has links to blob 3, thereby forming a new cluster . This four-point irreducible polymer contributes to the -kernel.
Other contributions to the dressed propagators. This dressed bubble diagram involves strongly connected two-point polymers.
1$$4$$2$$3$$7$$5$$6$$j_{+}
Fig. 3.2.12. Dressed propagators.
At scales , we keep the idea of a third-order cluster expansion, but explore only links starting from four-point polymers with low transfer momentum of scale , and connecting them to arbitrary polymers, except four-point polymers with transfer momenta of higher scale .
B. Scale complementary expansion.
The scale horizontal cluster expansion.
We resume from expression (2.107) for .
Step . Preliminary second-order expansion in isolated boxes (see Fig. 3.2.8). Let . Expand , see (3.18), into
[TABLE]
The first non-trivial term is seen as a four-point polymer. Neglected terms involve at least 8 fields in , hence are considered as polymers with external legs.
Step . Splitting phase for four-point polymers of type (see Fig. 3.2.9).
Let be a four-point polymer of type , with external structure as in (3.19), with external sectors . By definition, . The argument below can be repeated without any modification at scale with , so we state the general version. Consider first the left part of . Take the Fourier transform,
[TABLE]
or conjugate. Insert inside the last expression the kernel . Because , there are two cases: (i) either the transfer momentum has bosonic scale , and ; or (ii) , and . Hence
[TABLE]
with a convolution kernel essentially supported on .
As a result, has been rewritten as the sum of two terms interpreted as follows:
(i) a polymer with upper transfer momentum scale (term with in factor);
(ii) a polymer with lower transfer momentum scale (term with in factor).
Step . External leg displacement for four-point polymers of type 2 with upper transfer momentum scale (see Fig. 3.2.9).
Consider four-point polymers obtained in (i) just above. Displace the fields , resp. on the lower line to the location of the fields , resp. on the upper line, see Fig. 1.7.5 and eq. (1.112), namely, replace by and similarly for .
Step . Scale expansion.
One step of expansion. Fermionic expansions have produced scale clusters, as described in §2.3.1. Consider two scale clusters, , and , where are arbitrary polymers, except four-point polymers with low transfer momentum scale . Introduce an interpolation coefficient between and , namely, multiply by (); and differentiate in (2.107) w.r. to . The connecting operator associated to the cluster link as in Proposition 5.2,
[TABLE]
either (i) transforms a pair of dangling fields into a cluster link, or produces (ii) one or (iii) two extra vertices, one in , one in . The procedure may be iterated an arbitrary number of terms by differentiating w.r. to We shall actually do so at most times per pair of boxes (see below).
We say that * and are strongly connected at scale * if there exist at least 3 cluster or complementary cluster links between and .
We now introduce a linking number , which is initially zero, and is incremented as links are added. If there is one link between and , we set . If there are two links, set .
We repeat the above procedure until LN()=, so that and are strongly connected. By choosing one of the terms in the cluster expansion in which an interpolation coefficient has been set to [math], one also gets factorized expressions of the form , times a quantity averaged w.r. to a measure , with , for which fields in and are independent one from the other. The quantity is either 1 (in which case and are totally factorized), a single propagator, or a product of two propagators. In the latter case, we consider that we have one external Cooper pair propagator connecting to .
After completing the cluster expansion between and , we define a new set of objects as follows. First of all, we are not really interested in the linking number between two boxes belonging to two given polymers , but rather in the total number of links connecting and (see Figure in the introduction to section 3), which is defined as
[TABLE]
Then we stop the cluster expansion between and as soon as . This mutualized stopping rule is more economical, sufficient for our purposes and consistent with the general organization of the expansion in terms of polymers.
Let us come back to the outcome of the first step of expansion. As in §2.3.1, we describe (super)-clusters, called objects, produced by connecting clusters. Scale clusters are objects. If LN(,, then clusters are also objects; otherwise a new object has been formed. The expansion goes on by considering links between the object containing and clusters . Explore in this way all clusters . In the end, one has an object containing , and a set of clusters which are not strongly connected to . Then one starts from a new box and explores in the same manner its links to clusters , and so on. Here is one example.
\Delta_{1}$$\Delta_{2}$$\Delta_{3}123
45312
Fig. 3.2.13. Third order cluster graph.
Blue lines are cluster links. Because we expand to third order, we obtain in general a graph, not a tree. However, it is easy to extract a tree out of this graph (in red on the figure). Decide that every third link between two objects is a tree link (red link number 1). It may happen that, upon introducing a third link (here between and ), and forming the corresponding cluster (here ), previously non strongly-connected objects become strongly connected, here and . One then adds a tree link, here between and ; for the sake of clarity, it is convenient to represent it as a link between and one of the boxes of , here (but it is an arbitrary choice) . Thus a new cluster is formed, and so on. These rules generate a tree.
C. Complementary expansions of scales . The series of steps Step Step Step may be repeated at scale until all Cooper pairs of all scales have been connected. Note that the preliminary second-order expansion in isolated boxes (Step ) is performed only once, at the beginning of the scale complementary expansion. As discussed above, Step (splitting phase for low transfer momentum for four-point polymers) is applied to Cooper pairs belonging to four-point polymers of type 2 (in the sense of Definition 3.1) which already have in factor, therefore there is no need to add a further type specification at scales . Step is somewhat simplified w.r. to the case , since (as explained above) there is no need to test pairings between scale-neutral polymers any more; so the exploration process starts exclusively from scale four-point polymers, whose sectors range not in , but in . Generally speaking, the scale complementary expansion inserts scale Cooper pairs and the four-point polymers they belong to inside a multi-scale tree (see Fig. 3.2.6). Also, newly developed Goldstone boson propagators of scales produce inductively scale contributions to the two-point function, e.g.
j^{\prime}_{\phi}$$j_{+}
Fig. 3.2.14. Scale dressing of the fermion propagator by internal Goldstone boson propagators.
yielding in the end of the story the dressed fermion propagator .
D. Resummation of Goldstone boson propagators. Once all Cooper pairs of all scales have been connected, the remains at each scale a sum of structures as in Fig. 3.2.3 or 3.2.4. However, it isn’t yet in the form of a geometric series; rather of a non-factorized expression of the form
[TABLE]
as in Definition 1.7 (and §1.7 for notations), except that is now computed as a product of two dressed propagators, . The displacement procedure of §1.7, see eq. (1.112), is defined in terms of isotropic ”micro-sectors” obtained by chopping the angular sectors as in the Remark following Definition 2.2. Error terms have by (1.114) a supplementary small pre-factor corresponding to the difference of scales between the external legs, and the internal fermions, which are of scale . One must first
- (i)
resum over the locations of all intermediate scale-neutral polymers , implying in particular the possibility of overlaps between polymers. This is a standard problem, solved by means of a Mayer expansion as in §2.3.5.
Once this is settled, one must still take into account the dependence of in and . In order to extract factorized contributions, we mimic the procedure sketched in a perturbative context in §1.7 under the name of Fermi and s-wave projections. As suggested in Definition 1.5, it is rather an averaging procedure, a double one in fact. we must get rid of the dependence on the norms , and (if ) on the angles . (If then depends separately on and ). Let .
- (ii)
(norm averaging procedure) The kernel has a support centered on , resp. . We absorb the intermediate -functions by rewriting products
[TABLE]
as
[TABLE]
where
[TABLE]
Define , to be the projection on the function w.r. to the -norm in momentum space,
[TABLE]
Let be the volume of for any sector . Then a short-hand for (3.26) is simply
[TABLE]
Define
[TABLE]
;
[TABLE]
and a similar formula by symmetry for K^{\alpha_{i,1},\alpha_{i,3}}_{0|1}(p_{1},p_{3}):=\Big{(}P^{\alpha_{i,1}}\otimes(1-P^{\alpha_{i,3}})\Big{)}(K)(p_{1},p_{3}); and finally,
[TABLE]
By -orthogonality, (3.30) reduces to a sum over chains of intermediate states,
[TABLE]
Also – and this is the crucial point –, the covariance kernels are nearly constant,
[TABLE]
hence all kernels except have a small prefactor, namely, for , and for . We may therefore resum all contributions corresponding to chains of length , i.e. including at least one intermediate state , and reconsider them as scale-neutral four-point polymers with a small factor
[TABLE]
Thus we are left (up to a redefinition of including these chains, for which we choose not to introduce a new notation) with purely chains: there is no norm dependence any more.
- (iii)
(angle averaging procedure) One must still a priori sum over intermediate sector variables . Assume first that . Then , depends only (by global rotation invariance) on angle differences. Using a discrete Fourier transform, , so that
[TABLE]
we immediately see by standard orthogonality properties of Fourier series that only diagonal terms contribute. The terms are called resp. -, -, -wave contributions, following the standard chemists’ notation for angular momentum. Dominant contributions here are due to bubbles, which are in the -wave because the interaction itself is. Hence the sum of all contributions other than that in the -wave has a supplementary factor , and may be reconsidered as in (ii) as scale-neutral four-point polymers with a small factor . Thus there remains only the -wave contribution.
One then obtains a geometric series resummed as the single-scale Goldstone boson propagator.
An easy supplementary step consists in
- (iv)
resumming single-scale Goldstone boson propagators connecting a polymer to itself. Namely, these tadpole like bosonic insertions are convergent.
- (v)
The last step consists in resumming by scale induction two-point functions of single-scale Goldstone boson propagators. For instance, the following scale two-point function may be built out of the bosonic tree of Fig. 3.2.6 ”glued” with its mirror image,
j_{+,1}$$j_{+,2}$$j_{+,3}
Fig. 3.2.15. Two-point function of the single-scale Goldstone boson propagator.
Resumming the geometric series of such diagrams with lowest transfer momentum scale yields the Goldstone boson propagator, denoted by a wiggling line
The final outcome of this series of expansion is a multi-scale tree , with scale-neutral polymers as vertices, and Goldstone boson propagators as edges. Vertices have an arbitrary coordination number . The case is excluded by the gap equation; the case is excluded since polymers in sandwich between two Goldstone boson propagators have already been included in the geometric series. The only cases when or are allowed is when external legs – either isolated fermions or Cooper pairs – are connected to the polymer. See Figure below.
\times$$\Psi_{ext}$$+$$\Psi_{ext}$$\times$$\times$$(\bar{\Psi}{\mathbb{\Gamma}}\Psi)_{ext}
Fig. 3.2.16. Final bosonic multi-scale tree.
E. General principle of the bounds.
The next lines are concerned with the power-counting of a polymer with Cooper pair external structure
[TABLE]
with , , see splitting phase, step #1 (3.22). As in the fermionic scales, see §2.3, one is led to separate the absolute power-counting, obtained by assuming that , from the relative power-counting, i.e. the power of the prefactor in , where is the highest external scale.
It is instructive at this point to redo the power-counting of §2.3.2 for a connected fermion diagram with external Cooper pair legs located at , and with scales , i.e. (considering the Fourier transformed graph) such that transfer momenta . The worst case is when Cooper pairs are connected by Goldstone boson propagators , scaling like \big{(}\Gamma_{\phi}2^{-j_{+}/2}\big{)}^{2}. Hence the power-counting of a Cooper pair , see (3.38), summed over all possible sectors scales, like ; note that this contrasts with the naive scaling (2.13) for the product of two fermions, . If , then the -integration in a box costs a volume factor . Thus, the relative power-counting for a polymer in with highest external scale is
[TABLE]
=relative degree of divergence, which implies that only diagrams with external Cooper pair legs are potentially divergent when with fixed.
Let us now bound instead the absolute power-counting of non-amputated polymers (i.e. including external bosonic legs) for is fixed; the total power-counting of the full diagram for arbitrary is by construction the product of the relative power-counting (3.39) by the absolute power-counting which we shall presently bound.
We restrict here to the case of a fermionic polymer in with external bosonic legs, hence external fermionic legs. The corresponding amputated power-counting (as shown in (2.123)) is . Polymers in , with have been resummed into the fermionic two-point function, hence we may assume that either or ; thus, , and the worst bound is obtained when , which we now assume. The power-counting of the corresponding full diagram is obtained by multiplying by the weight of the external bosonic legs, and by the volume of a scale box, all together (taking into account that ) at most:
[TABLE]
For (excluding the case of a Goldstone boson one-point function, or of a Goldstone boson propagator, for which ) this is small.
Ward identities – proved in §3.3 – will show that Goldstone boson one-point function vanish, and that diagrams with are actually . They will also allow to compute the absolute power-counting of non-amputated polymers in including scale fermions.
3.3 Ward identities
*The general purpose of this subsection is to obtain supplementary spring factors for external Cooper pair legs of -point functions or bosonic polymers with lowest internal scale and highest external scale .
- As is apparent from the previous discussion in §3.2 E., only a small (actually two, as we shall see), finite number of spring factors is really required to show that polymers converge. In this subsection, we shall be uniquely concerned with showing that the naive relative degree of divergence , of (3.39) can be enhanced to , with (called: relative gain) large enough so that .
In guise of introduction: gauge symmetries. The Grassmann measure is invariant under number symmetry,
[TABLE]
or equivalently
[TABLE]
if is a constant. Furthermore, the interaction is invariant under gauge transformations
[TABLE]
where now is allowed to depend on space-time.
Multiplying by and letting , we obtain an infinitesimal transformation such that
[TABLE]
[TABLE]
where {\boldsymbol{\Gamma}}^{\perp}(\theta):=\left(\begin{array}[]{c}-\Gamma_{2}(\theta)\\ \Gamma_{1}(\theta)\end{array}\right) is the image of the vector by a rotation of angle , and the associated off-diagonal Hermitian matrix. The inserted kernel is essentially the inverse of the fermion covariance kernel. This is true up to the ultra-violet cut-off , however. If is constant, then the multiplication by commutes with Fourier cut-offs; however, this is not true for non-constant gauge transformations. This generates error terms denoted by ”err.” in the expressions below, which will be shown in §5.5 to be of the same order or smaller than the terms in (3.45). Exponentiating the infinitesimal transformation , one obtains a rotated order parameter, ; equivalently (see (0.84)), gauge transformations locally rotate -vectors, .
Applying gauge transformations specifically to products of Cooper pair fields
yields up to error terms, denoted by ”err.”,
[TABLE]
where acts as on , and as on .
We are however interested in infra-red cut-off quantities obtained by integrating Goldstone boson propagators with scales , and leaving out -kernels with scales . This introduces further error terms due to the infra-red cut-off on the -kernel this time. Otherwise (3.46) still holds true with replaced by .
We choose to apply Ward identities to *amputated bosonic -point functions
*. By definition, compare with (2.83) and (3.22),
[TABLE]
where is the generating function of connected diagrams excluding four-point diagrams with transfer momentum scale , and additive sources, . The multi-scale cluster expansions of §3.2 downto scale typically produce – after displacing the -legs of each external Cooper pair – polymers with low-momentum external -kernels, summed as
[TABLE]
However, there is the possibility that some of the Cooper pairs are in the -channel, or that the two fermions composing the Cooper pair connect to two different polymers. These ”massive” cases are treated in C. below. Using linear combinations, one may restrict to or even to .
A. One-point functions and gap equation. Let us for a start consider the amputated one-point function , or . The latter function is conveniently rewritten as
-$$\times$$\sigma(\theta^{\prime})$$+$$\sigma(\theta^{\prime})•+$$-$$\times$$\sigma(\theta^{\prime})$$+$$\sigma(\theta^{\prime})•
Fig. 3.3.1. Amputated one-point function.
where gray blobs denote strongly connected (see 3.2), bosonic one-point functions, excluding four-point
subdiagrams with transfer momentum scale ; and
are two-point vertex insertions coming from off-diagonal propagators computed in the -channel and with the -coefficient computed at scale , therefore evaluated as when , resp. [math] when .
Consider to begin with the one-point function in the parallel direction . Setting to zero the first line of Figure 3.3.1
-\Gamma^{j_{+}\to}$$+$$\sigma(\theta)•
Fig. 3.3.2 Gap equation, one-point version.
yields an implicit equation for ,
Definition 3.2** (gap equation, one-point version)**
[TABLE]
being the sum of strongly connected (s.c.), bosonic one-point function diagrams (drawn above as a gray blob), which is – as we shall presently prove, using a Ward identity – exactly equivalent in the infra-red limit to its two-point version (a refined version of Definition 1.9) stating that the Goldstone boson propagator has a pole in the transverse direction. Let us note that, provided the above gap equation is verified, *all lines
- in Fig 3.3.1 vanish up to error terms due to external leg displacements which are at most of order , hence vanish in the infra-red limit (), e.g.
-$$\times$$\sigma(\theta^{\prime})$$+$$\sigma(\theta^{\prime})•
Next, we apply to the one-point function in the parallel direction, see Fig. 3.3.2, the global (constant) Ward identity associated to . We immediately obtain
•,
in other words,
[TABLE]
the one-point function in the perpendicular direction vanishes.
Finally, we consider the sum of strongly connected (s.c.) contributions to the one-point function in the perpendicular direction . Applying to it likewise the global (constant) Ward identity associated to , we get
-$$\sigma(\theta)$$+\int d\xi$$\sigma(\theta+\frac{\pi}{2})•••\xi_{3}$$\xi_{4}$${\mathbb{\Gamma}}^{\perp}•\xi$$=0
Fig. 3.3.3. Ward identity connecting one- and two-point functions of the Goldstone boson.
The left term is immediately seen to be , which is equal to by Definition 3.2. On the other hand, the gray blob in the right term is now a strongly connected four-point diagram. Displacing its external leg to following the displacement procedure of the previous subsection yields a product of a -kernel by a bubble in the -channel. The error term due to the displacement procedure has by symmetry an extra spring factor . One has found:
[TABLE]
As asserted above, this is – up to error terms in – equivalent to the gap equation of §1.7. For convenience, we shall define in §2.4 – see Definition 3.3 – in such a way that the infra-red cut-off Goldstone boson propagator in the -channel, , has a pole precisely at . This implies in turn that the one-point function version of the gap equation, see Definition 3.2, is satisfied up to error terms,
[TABLE]
B. Local Ward identities. Let . We shall now apply to
[TABLE]
a general gauge transformation along . Note that we do not restrict to strongly connected diagrams any more. At this point it is useful to introduce the following general notation,
[TABLE]
generalizing (3.47), where are general Hermitian (not even necessarily off-diagonal) two-by-two matrices. Gauge transformations rotate the matrix to The lower pair of indices points to the fact that all angles , are in the perpendicular direction.
Then – generalizing Fig. 3.3.3 – we obtain the following order 1 local Ward identities,
[TABLE]
Take such that , . Then the first term on the r.-h.s. of (3.67) vanishes, so that
[TABLE]
Let us now take the functional derivative of (3.74) w.r. to :
[TABLE]
Finally, integrating w.r. to the external legs except , one finds in particular by integrations by parts,
[TABLE]
The second-order operator contributes a factor of order , a squared spring factor; on the other hand, the first-order operator contributes a single spring factor , which is not enough for our purposes. Fortunately, the whole procedure may be iterated twice, using successively two derivatives , with and , and taking the functional derivative . Considering only time-derivatives, one has schematically
[TABLE]
[TABLE]
Thus , has been rewritten as a sum of four terms involving I^{j_{+}\to}_{n}\left(\begin{array}[]{ccccc}{\mathbb{\Gamma}}^{\perp}&\cdots&{\mathbb{\Gamma}}^{\perp}&*&*\end{array}\right) with , and product differential operators of the type \big{(}\frac{1}{\Gamma_{\phi}}\frac{\partial}{\partial\tau_{n-1}},\frac{1}{2m^{*}\Gamma_{\phi}}\boldsymbol{\nabla}_{\xi_{n-1}}^{2}\big{)}\times\big{(}\frac{1}{\Gamma_{\phi}}\frac{\partial}{\partial\tau_{n-1}},\frac{1}{2m^{*}\Gamma_{\phi}}\boldsymbol{\nabla}_{\xi_{n-1}}^{2}\big{)}, therefore of order . Thus one has decreased the relative degree of divergence of the original -point function by at least two units. This is sufficient to show the convergence of amputated -point functions provided , since then the corrected degree of divergence is .
C. Massive configurations. We have not considered in previous computations the following two cases:
- (i)
-point functions of the type
[TABLE]
with , ;
- (ii)
-point functions such that at least one of the external Cooper pairs connects to two different polymers, as in Fig. 3.2.7.
Both cases are massive, featuring at least one external fermionic propagator with quasi-exponential decay at distances , implying that the corresponding polymers are convergent, and furthermore, their contribution decreases exponentially as . More precisely:
- (i)
external Cooper pairs along the parallel () direction form a -kernel which is bounded as in (1.94).
- (ii)
Referring to the notations on Fig. 3.3.7, we have the following scaling factors for the compound polymer connected by two fermions (assuming fermion angular sectors are fixed):
– for the integral over ;
– for the integral over constrained by the fact that ;
– for the two fermionic propagators;
– a scaled quasi-exponential decay factor of scale for the propagator with momentum , resp. for the propagator with momentum ;
– and finally the squared spring factor coming from the Ward identities. The total factor is , multiplied by the two quasi-exponential decay factors. The first factor is the expected weight for four fermions . The bad scale decay factor for the second propagator yields an extra volume factor compared to the expected volume implied by a scale decay factor. Decomposing into , one sees that the extra volume factor is compensated, and that one has a further spring factor which makes it possible to sum over the scale .
3.4 Bosonic fixed-point procedure and low-energy multi-scale bounds
Consider a bosonic multi-scale tree as in Fig. 3.2.15. Bounds for fermionic polymers with low-momentum Cooper pair legs in 3.2 E. was based on the naive power-counting for a Goldstone boson propagator of scale ,
[TABLE]
therefore, per Cooper pair. Ward identities (see 3.3), as proved just above (see below (3.88)), imply however that one can get one or two supplementary spring factors per polymer, thus modifying the scaling of the corresponding Cooper pairs,
[TABLE]
Thus the sum over bosonic multi-scale trees may be bounded as follows. We first need some notations. Let be the scale-neutral polymers. By assumption, we fix the location of one of the vertices of , so that their power counting is , where . The scale Goldstone boson propagators merge these polymers into bosonic polymers with lowest scale . The power-counting of each of these is computed with the same convention, namely, by fixing the location of one of the vertices. Repeating this procedure, we get a sequence of polymers , which becomes stationary when all Goldstone boson propagators have been integrated into polymers.
We proceed as in §1.3 (compare with (2.123), (2.124), (2.125)).
General bosonic multi-scale bound. The actual energy gaps , resp. non-linear sigma model coupling constant , resp. velocity are obtained inductively as the limits of the converging series , resp. , resp. resp. .
We start in A. by stating the difference estimates. Such estimates make it possible to solve inductively a sequence of approximate gap equations (see Definition 3.3) yielding a Cauchy sequence . As an intermediate step, we first prove in B. polymer bounds by scale induction. Difference estimates are finally proved in C..
A. Statement of difference estimates and gap equations. For each , we search for the solution of largest module of the
Definition 3.3** (scale gap equation)**
[TABLE]
compare with Definition 1.9 and (1.121). The parameters of the Goldstone boson propagator are defined as in Lemma 1.3 and Definition 1.4 by
Definition 3.4** (scale Goldstone boson propagator
parameters)**
Let
[TABLE]
[TABLE]
[TABLE]
[TABLE]
and
[TABLE]
The scale gap equations are solved by induction on using difference bounds of two types:
- (i)
(scale differences) let such that , then
[TABLE]
[TABLE]
with ;
- (ii)
( dependence) let such that , and , then
[TABLE]
[TABLE]
If difference bounds hold, then subtracting (3.99) from (3.99) and applying the implicit function theorem yields after an elementary computations
[TABLE]
where . In particular, .
Also, by elementary computation, , , and
[TABLE]
[TABLE]
[TABLE]
as stated in Theorem 2.
B. Proof of polymer bounds.
Let , then we want to prove, in analogy with the bounds (2.123),(2.123),(2.125) of section 2 for fermionic polymers in , , but now with , , and bosonic polymers obtained by merging polymers also including scale fermions, following the procedure of §3.2:
General bosonic multi-scale bound.
For every belonging to the set of polymers with lowest scale and external legs,
[TABLE]
The exponent of has been chosen in such a way that the factor , times the weight of the external structure of a polymer with external Cooper pair legs (see below), is scale independent. Then, if ranges in the set of polymers in with vertices,
[TABLE]
from which, finally, summing over all polymers in with one external leg located at ,
[TABLE]
The proof of bosonic polymer bounds (3.113),(3.114),(3.115) is by induction on . Note first that (3.113),(3.114),(3.115) follow from the general fermionic multi-scale bound (2.123),(2.123),(2.125) if one sets , which means that all fermionic polymers (i.e. all polymers having all internal propagators of scales ) are assimilated to polymers with lowest scale ; this is a convenient choice for the induction, since field scaling factors (see below) have a scale-independent prefactor proportional to . Namely, identifying the bosonic legs of a fermionic polymer with
[TABLE]
fermionic legs, we remark that
[TABLE]
which shows that (3.113) holds with prefactor instead of for a fermionic polymer in . Since , the fermionic bound increases with , hence (3.113) is all the more true for a fermionic polymer in if . Finally, for , a correcting factor is required to obtain the fermionic bounds , which is compensated by since .
We now suppose as in §2.4 that is obtained by merging at scale polymers with lowest bosonic scales and external Cooper pair legs.
Ward spring factor. As discussed in the introduction to this subsection, we have an extra spring factor for each polymer .
Volume factors. External legs of which are contracted at scale are integrated in a box of scale . This yields a volume factor .
Field scaling factors. Each external Goldstone boson half-propagator of which is contracted at scale comes with a supplementary prefactor .
The product of the volume factors and of the field scaling factors gives the above mentioned weight . Taking these into account for the polymers and (in the denominator, in order to avoid double-counting) for , yields the general multi-scale bound (dismissing dimensional factors)
[TABLE]
where
[TABLE]
is the scale-dependent part of (3.113). We show that this is , where . For this we consider separately scale-independent factors, scale factors, and rescaling factors:
- (i)
(scale-independent factors) Each polymer is connected at least to one other polymer, hence at least one of the external legs of each polymer does not belong to the external legs of . Thus . The total exponent of the scale-independent factor in (3.118) is
[TABLE]
- (ii)
(rescaling factors) rewrite as times
. Since by hypothesis , the exponent is . The remaining spring factor makes it possible to sum over all scales .
- (iii)
(scale factors) For each , the left-over factor in (ii) cancels with the scale-dependent part of the weight of . Therefore, there remains only the required overall prefactor .
C. Proof of difference estimates.
- Difference estimates for . The difference , see (3.105), involves terms of three types which we present by decreasing order:
- (i)
pure scale irreducible four-point polymers, like the first term of Fig. §3.2.3. These contain no inner Goldstone boson, hence enjoy by Ward identities a spring factor , multiplied by the scaling (2.123) of an amputated, four-point, fermionic polymer of lowest scale , ; all together, one gets the expected factor .
The two remaining types are multi-scale bosonic trees obtained in the last two resummation steps, Step (iv), resp. Step (v) of p. 3.1, involving at least one Goldstone boson of scale :
- (ii)
(internal Goldstone boson propagator)
\mathbb{P}_{1}^{(j_{+}-1)\to}$$j_{+}\!-\!1
Fig. 3.4.1. Case of a scale internal Goldstone boson propagator.
involving one internal Goldstone boson propagator connecting a polymer to itself. Leading terms are as in the above Figure, leading to a contribution larger than in case (i),
[TABLE]
that is, (Amplitude of a polymer with 4 external Cooper pair legs \mathbb{P}_{1}^{(j_{+}-1)\to}\big{)}\times\big{(}Scale Goldstone boson propagator). See (3.113).
- (iii)
(Goldstone boson connecting two different polymers)
\mathbb{P}_{1}^{(j_{+}-1)\to}$$\mathbb{P}_{2}^{(j_{+}-1)\to}$$j_{+}\!-\!1$$j^{\prime}_{+}••
Fig. 3.4.2. Case of a scale Goldstone boson connecting two different polymers.
involving two Goldstone boson propagators connecting the same two polymers , one of scale , the other of scale . Spring factors make it possible to sum over for fixed, hence we can assume that .
Other more complicated diagrams with as in B., connecting , with and internal Goldstone boson propagators, are possible. The power-counting of such diagrams is as in the General bosonic multi-scale bound, eq. (3.113). Eq. (3.120) and point (i) around it show that the largest power-counting is obtained for , that is, for the situation depicted in the above Figure, for which one obtains a contribution to of order at most
[TABLE]
that is, (Integration volume of w.r. to fixed )\times\big{(}Product of the amplitudes of the polymers with 3 external Cooper pair legs \mathbb{P}_{1}^{(j_{+}-1)\to},\mathbb{P}_{2}^{(j_{+}-1)\to}\big{)}\times\big{(}Squared scale Goldstone boson propagator).
- Difference estimates for . These are obtained in a straightforward way using (3.100), (3.101), (3.102), (3.103) from \nabla_{q}^{2}\Big{\{}\big{(}\bar{\Pi}_{q}^{j_{+}\to}(\Gamma)-\bar{\Pi}_{q}^{(j_{+}-1)\to}(\Gamma)\big{)}{\cal A}_{q}(\Gamma^{j_{+}\to},\Upsilon_{3})\Big{\}}\Big{|}_{q=0}. Because the bubble diagram is made up of fermion propagators with a scale effective infra-red cut-off due to the -term in the denominator of (1.20), operators change its scaling by a factor . So let us study instead the effect of the action of on polymers of the type studied in 1. As a general rule, gradients w.r. to the transfer momentum of a diagram with two external Cooper pairs can be evaluated by choosing a line of propagators connecting the two external vertices ,
q+\delta q$$q+\delta q•••••••p_{2}+\delta q$$p_{1}+\delta q
Fig. 3.4.3. Changing the transfer momentum flowing inside a diagram.
see Fig. 2.3.3. Translating by inside the diagram is equivalent to substituting Taking two successive derivatives w.r. to yields sums of terms featuring the product \big{\{}\prod_{k\not=i,i^{\prime}}C^{*}_{\theta}(p_{k})\big{\}}\ \nabla C^{*}_{\theta}(p_{i})\,\nabla C^{*}_{\theta}(p_{i^{\prime}}). In principle, the inverse Fourier transform of a gradient fermion propagator , , is , which is bounded as in (2.11), but with a further prefactor scaling at most like , induced by the scale quasi-exponential decay. This implies directly for scale irreducible diagrams of type (i) as above a relative scaling of w.r. to of order at most. Consider now a diagram or type (ii) or (iii), involving a scale Goldstone boson propagator, assume that the line contains a Goldstone boson propagator of scale . Differentiating w.r. to yields (by the resolvent identity)
[TABLE]
Proceeding by induction, we see now that involves a supplementary scaling factor , exactly as .
Concluding, and letting , and similarly : (see (3.100), (3.101), (3.102), (3.103)) scale at most like , namely, (scaling factor due to )scaling factor of , whence the estimates (3.106).
- -dependence. Exactly as for , and for the same reasons (see arguments in 2.), the -derivative operator costs a scaling factor or (reintroducting dimensional constants) . Thus
[TABLE]
Summing over scales yields (3.107). The other estimate (3.108) follows as in 2. from the scale effective cut-off scale in the denominator of (1.20).
Error terms \big{|}{\mathrm{Pre}}\Sigma(\Gamma^{j_{+}\to};\xi-\xi^{\prime})-{\mathrm{Pre}}\Sigma(\Gamma^{(j_{+}-1)\to};\xi-\xi^{\prime})\big{|}, (see §2.3.3), modifying the value of the pre-Goldstone boson propagators inserted inside fermionic polymers, are similarly bounded.
Next subsections are dedicated to the proof of Theorems 2–5 giving estimates for various -point functions. The general scheme of proof is as follows. Define infra-red cut-off -point functions
[TABLE]
for . When , is the sum of entirely fermionic diagrams obtained by the fermionic expansion of \Big{\langle}\prod_{i=1}^{n}\Psi(\xi_{i})\prod_{i=1}^{\bar{N}}\bar{\Psi}(\bar{\xi}_{i})\Big{\rangle}_{\theta;\lambda}, i.e. of all fermionic polymer contributions including no Goldstone boson propagator. When , the sum is rather obtained from the fermionic polymer expansion of \Big{\langle}\prod_{i=1}^{n}\Psi(\xi_{i})\prod_{i=1}^{\bar{N}}\bar{\Psi}(\bar{\xi}_{i})\Big{\rangle}_{\theta;\lambda} by substituting to all Goldstone boson propagators .
To go from the purely fermionic contribution to the -point function \Big{\langle}\prod_{i=1}^{n}\Psi(\xi_{i})\prod_{i=1}^{\bar{N}}\bar{\Psi}(\bar{\xi}_{i})\Big{\rangle}_{\theta;\lambda}=\lim_{j_{+}\to+\infty}f^{j_{+}\to}(\boldsymbol{\xi},\boldsymbol{\bar{\xi}}), one proceeds by induction, extending downwards the polymers to the bosonic scales. Namely, let be the fermionic polymers of a given contribution to the -point function, and decide that and are connected at scale (which we denote by ) if there exists a scale propagator connecting to . Then we define
[TABLE]
[TABLE]
External legs of a polymer are then made up of external -propagators with scales . It is important to realize that no supplementary cluster expansion is needed here: going down the scales inductively, polymers simply coalesce, until there is only one polymer left at some scale .
3.5 Proof of Theorem 2: Cooper pair -point functions
In this section we prove Theorem 2. Fix and . The main contribution to
[TABLE]
is as in Fig. 3.1.2, except that -kernels are replaced by -kernels,
••••••••••••\xi_{1}$$\xi_{4}$$\xi_{2}$$\xi_{5}$$\xi_{3}$$\xi_{6}$${\mathbb{\Gamma}}^{\perp}$${\mathbb{\Gamma}}^{\perp}$${\mathbb{\Gamma}}^{\perp}$${\mathbb{\Gamma}}^{\perp}$${\mathbb{\Gamma}}^{\perp}$${\mathbb{\Gamma}}^{\perp}+ perm.
Fig. 3.5.2. -point function of Cooper pairs.
where ”perm.” indicates the sum over all possible pairings. Each bubble diagram is equal to by the gap equation. This yields the leading term in Theorem 2. Corrections are due to more complicated fermionic diagrams than bubble diagrams, possibly containing other Goldstone bosons. They all contain at least one more vertex, implying a further prefactor in .
3.6 Proof of Theorem 3: fermion quasi-exponential decay
Main diagrams contributing to (0.71), resp. (0.72) in Theorem 3 are
••\xi$$\xi^{\prime}
Fig. 3.6.1. Two-point functions for isolated fermions.
••\xi_{1}$$\xi_{2}••\xi_{3}$$\xi_{4}
Fig. 3.6.2. Four-point functions for fermions either in non-Cooper pairing, or in parallel Cooper pairing.
in the limit when and , resp. .
This, and the quasi-exponential decay of the -kernels (see Proposition 2.3) and of the -kernel (see (1.94) and discussion in §3.3 C.), accounts for the quasi-exponential decay of such correlation functions.
3.7 Proof of Theorem 4: phase transition
In this section we prove Theorem 4. By global rotation invariance, we may assume that . We want to prove
[TABLE]
On the other hand, , as follows from the vanishing of the bosonic one-point function in the perpendicular direction (see (3.50)), proved using a Ward identity. Note that, for any cut-off scale , the gray blob in the figure above (3.50) cannot be connected to the external Cooper pair by a Goldstone boson by momentum conservation, hence the limit can be taken; also, error terms in vanish in that limit. All together, this implies Theorem 4. Clearly (by translation invariance) is a constant.
Main contribution to is the bubble diagram
•\xi$$\times$$\Gamma
which is indeed equal to by the gap equation. By number conservation (see e.g. discussion before Definition 1.5), all diagrams contributing to contain at least one off-diagonal fermionic propagator, hence they all have at least one in factor. Since the one-point function of the Goldstone boson vanishes by construction, the Cooper pair may not be directly connected to a Goldstone boson propagator. On the other hand, more complicated fermionic diagrams than the bubble diagram, potentially containing Goldstone bosons are possible. However, they all contain at least one more vertex, implying a further prefactor in .
3.8 Proof of Theorem 5: local transverse behavior of the Goldstone boson
We are interested in this subsection in large deviations for the random variable
[TABLE]
– interpreted as the integral over of the projection along the direction of the Goldstone boson –, where , and is some volume comparable to a box of transition scale , say, is a roughly cubic connected union of boxes with dimensions , . By rotation invariance, we may and shall assume that .
We let be the ratio . The validity of the argument below extends in part to , so we simply assume that . Note that .
*Standard large-deviation theory *(see e.g. Cramer’s theorem in [13], §2.2) implies the following large deviation bound. Let be the average of in the symmetry-broken measure along ; as proved in §3.7, . Let
[TABLE]
be the log-cumulant of , and
[TABLE]
its Legendre transform. Then, for ,
[TABLE]
Similarly,
[TABLE]
We concentrate on the consequences of (3.133) in the sequel; similar bounds, obtained by substituting , hold for .
Note that and (by Jensen’s inequality) for all . As a consequence, , so that , and , where
[TABLE]
is the unique realizing the maximum of the function ; in particular, . Thus
[TABLE]
where is the connected expectation w.r. to the perturbed measure .
Now, define
[TABLE]
By (3.135),
[TABLE]
equivalently, .
By construction, is equal to the sum of truncated expectation values . If , the main contributions to are in the form of bubble diagrams connected through scale Goldstone bosons to a fermionic polymer of scale , represented as the central gray blob in the Figure below.
•••••
Fig. 3.8.1. Main contribution to for .
Similarly, the truncated expectation may be expanded into .
Let
[TABLE]
be the Cooper pair two-point function (see (3.16) for the intuition, or Theorem 2 for a more precise statement) integrated over the volume . Since , .
A diagram as in the above Figure contributes at most times : (i) per Goldstone boson propagator with one fixed extremity, terminated in the other direction by a bubble, all together: ; (ii) (power-counting (2.123) for an amputated fermionic diagram with lowest scale and external legs), times a non-dimensional localization factor breaking translation invariance, , all together : . Hence, assuming (which will be checked later on by self-consistency as soon as )
[TABLE]
where if .
On the other hand, the terms of order are resp. and . Concluding:
[TABLE]
Comparing with (3.138), we get, neglecting the error term in (3.141):
[TABLE]
Taking this as the value for , we check that the error term in (3.141) is , which is \ll\Big{|}\langle X_{\Omega}\rangle_{\lambda,\Omega,A_{\eta}}-\bar{X}_{\Omega}\Big{|}=\eta\frac{\Gamma_{\phi}}{\lambda}{\mathrm{Vol}}(\Omega)\approx\eta\frac{\Gamma_{\phi}}{\lambda}\ell^{3}\frac{v_{\phi}^{2}}{\Gamma_{\phi}^{3}}\approx\ell^{3}m^{*}\left(\frac{v_{\phi}}{\Gamma_{\phi}}\right)^{2}\frac{1}{g}\eta provided .
Now the covariance in (3.136) is similarly computed as
[TABLE]
where . The error term, computed using (3.142), is now , provided once again.
Concluding: assume , then
[TABLE]
For one has proved Theorem 5. For larger one observes that the large deviation rate (3.144) increases like , and not as expected from the discussion in the Introduction based on the effective potential approach. Note however that this volume cubic root prefactor is not reliable, since it is valid only for small deviations . Another method would be needed to study large transverse deviations of the Goldstone boson over a large volume.
4 Generalizations and perspectives
This work is meant to be the first of a series of articles investigating the superconducting phase transition and quantum fluctuations of superconducting materials in various situations. The spectacular experimental and theoretical advances of the last 30 years or so, see e.g. [27, 64], is a strong incentive to draw detailed predictions from our microscopic approach in a variety of situations.
A first generalization. As mentioned in the Introduction, our approach extends in a straightforward way to more general two-body potentials as in (0.12). Substituting the kernel to the delta-interaction only changes numerical constants, since the phonon interaction is a short-range one. The same conclusion holds if one adds a small two-body potential with high enough infra-red cut-off.
Three-dimensional theory. The phase space analysis of the 3D model is more complicated because the sector counting of Proposition 2.5 is not available, due to the fact that fixing two sectors does not fix the remaining two sectors of a vertex any more. Thus a precise sector-counting argument is needed already in Regime I. However, it can be proved perturbatively that there is just at most one sum over sectors per vertex in any dimension, as proved here in §2.3.2 by following a loop; we plan to extend this argument non-perturbatively (see [56] for a first step in that direction). After resumming explicitly the leading contribution in the form of chains of Cooper pair bubbles as in dimension 2, remaining terms should be small by virtue of a 1/N argument. Consequently, we expect the present proof to extend more or less straightforwardly to the 3d case.
BCS theory at positive temperature. The finite-temperature formalism using Matsubara frequences, as presented e.g. in Chapter 7 of the book by Fetter and Walecka, makes it possible to compute Green’s function of the system in thermal equilibrium in terms of a modified functional integral with discretized energies now ranging in the set , where is the inverse temperature. At first sight, the same scheme should produce similar conclusions when the lowest Matsubara frequencies are smaller than the energy gap , i.e. when . In the inverse regime (), Cooper pairing should be shown to be statistically unfavorable. Interesting things should occur in 3D around a critical temperature , where one should be able to investigate more precisely the quantum second-order phase transition, in connection to ideas due to M. Salmhofer [61, 62].
Real time. The Hamitonian dynamics may be discussed in a -dimensional quantum-field theoretic framework with real-time variable using advanced and retarded propagators at temperature zero, or more generally the Keldysh formalism [59, 50] at . The major difficulty is now that the infra-red singularity is carried by a paraboloid in energy-momentum variables, , a two-dimensional manifold instead of a circle.
5 Appendix
We collect in this section some technical formulas and bounds required for the proof. An essential bound is Proposition 5.3 below, ensuring the convergence of the series of perturbations.
5.1 Grassmann integrals
We briefly recall here usual conventions regarding Grassmann integrals, see e.g. [57], §9.5. By definition, if are Grassmann anticommuting variables, and is a polynomial function, then is the coefficient of in . In particular, if is a couple of conjugate Grassmann variables, i.e. , then
[TABLE]
and
[TABLE]
In case there are several couples , the above rules generalized easily using some linear algebra,
[TABLE]
The measure \frac{1}{\det(B)}\Big{(}\prod_{i}\int d\bar{\psi}_{i}\,d\psi_{i}\Big{)}\,e^{-\sum_{i,j}\bar{\psi}_{i}B_{ij}\psi_{j}} is called the (normalized) Grassmann Gaussian measure with covariance kernel , where .
5.2 Integration by parts formulas
Proposition 5.1
Let be a (fermionic or bosonic) Gaussian field with translation-invariant covariance kernel and functional measure . Then, letting be a functional of ,
[TABLE] 2. 2.
Consider a smooth variation of the covariance kernel, and the associated functional measures , with . Then
[TABLE]
The proof of these elementary identities in a quantum field theoretic context can be found e.g. in [70]. They imply in particular the following elegant formula for horizontal cluster expansions, called Brydges-Kennedy-Abdesselam-Rivasseau (BKAR) formula. For that we first need some definitions. Introduce a partition of into a disjoint union of ”objects” , which are clusters, i.e. connected unions of scale boxes; let denote the set of objects. (The simplest example is when objects are elementary scale boxes, and ). A weakening function is a function assigning a weight to each link between two different objects; by hypothesis, and . Alternatively, may be considered as a weakening function which is a constant on a box: letting being the cluster containing , . Then is the Hermitian kernel defined by
[TABLE]
A weakening function may be chosen independently for each scale . The associated Grassmann Gaussian measure is denoted by , where . We now present the outcome of a single application of the BKAR formula at scale , so that if . The result is in terms of a sum over scale forests, whose set is denoted by . Forests are disjoint unions of trees. We write if and are in the same connected component of .
Proposition 5.2** (BKAR formula)**
(see [54], Proposition 2.6) Let be a functional of the scale fermion fields. Then
[TABLE]
where
[TABLE]
and , being the infimum of the for running over the unique path from to in if , and [math] else.
5.3 Local bounds for determinants
We prove here the statement (2.88) made in §2.3.2 concerning bounds for diagrams with a large number of vertices per box. Before we do so however, let us discuss this problem briefly in general terms, and cite a result (Proposition 5.3 below) proved in a different context in [45]. The proposition is not directly applicable in our context, so we shall give an independent proof of (2.88), which however relies on the same arguments.
The expansion produces multilinear expressions of the type
[TABLE]
. Applying Wick’s formula and taking signs into account, one gets
[TABLE]
where
[TABLE]
Expanding naively the determinant, one gets a sum over terms indexed by permutations of the set of -indices, making the sum over naively divergent. Looking at this factorial into details, however, one sees that this apparent divergence comes from accumulations of large numbers of fields inside boxes. Namely, assuming e.g. that the ’s and ’s are all located in different boxes, Proposition 2.3 implies
[TABLE]
for some constant , since the sum between \Big{\{}\ \Big{\}} converges for . There remains to discard the possibility of local factorials in the numerator, i.e. combinatorial factors of the form to a certain power, where is the number of fields contained in a box . This does happen in the case of bosonic theories, which requires more elaborated, truncated expansions of the exponentiated action, in order to produce converging series. However, here, due to the fermionic character of the theory, Pauli’s principle stronly suppresses local accumulations of fields. As a matter of fact, one has the following
Proposition 5.3** (Local bounds for determinants)**
[45]** Consider an complex-valued matrix , of the same form as (5.12). Assume that the kernel is smooth, and that it has a scale quasi-exponential decay , i.e. that, for every integer and multi-index ,
[TABLE]
for some (compare with Proposition 2.3). Pave by cubes of side and define , resp. to be the number of ’s, resp. ’s inside . Then, for every couple of integers , there exists a constant , such that
[TABLE]
Proposition 5.3 applies in our context to the evaluation of determinants, with , and is the covariance kernel of the fields , accompanied by their sector oscillation, so that . The above Proposition has been proven only for fields with a single component (in our context, lying in a single angular sector), however the extension to fields with arbitrary angular sectors is straightforward (as can be checked from the arguments below).
As mentioned before, we cannot use this result directly in the context of §2.3.2, because we want to bound sums of diagrams with a given loop structure, sector assignment and choice of scale boxes for vertex locations. So let us rewrite briefly the main arguments of the proof of Proposition 5.3 and prove that they do apply to such sums of diagrams. Fix some integer (later on to be identified with some multiple of in (5.15)). The idea is, for each sector , to split each scale box into equal sub-boxes of side length scaling like , centered in , containing each fields, with . Then each -sector propagator line connecting a -field located in is rewritten using Taylor expansion as
[TABLE]
This leads to a rewriting of an individual Feynman diagram as a sum of Feynman diagrams with displaced vertices and differentiated -fields. Because of Pauli’s constraint, \Big{(}\nabla^{\kappa}(e^{-{\rm i}(p^{j,\alpha},\xi_{\delta})}\Psi_{\sigma}^{j,\alpha}(\xi_{\delta}))\Big{)}^{2}=0. This means that Feynman diagrams with ( being roughly the number of terms on the first line of (5.16)) involve a large number of -order derivative fields inside , enjoying by (2.11) a supplementary prefactor , multiplied by , all together a small factor . Multiplying all these factors yields an overall factor
[TABLE]
where is the total number of -fields, and . Summing all exponents corresponding to sub-boxes included in a given box yields for some constant , or equivalently, inverse factorials to a power .
5.4 Cooper pair bubble and -kernel estimates
Proof of (1.59). We need to show that, assuming ,
[TABLE]
Assume (the easy generalization to is left to the reader). Copying the derivation of (1.56), we have , with
[TABLE]
If both signs are equal, then f_{\boldsymbol{q}}(\boldsymbol{p})=\Big{|}\frac{1}{e^{*}_{|\Gamma|}(\boldsymbol{p}+\boldsymbol{q})+e^{*}_{|\Gamma|}(\boldsymbol{p})}\Big{|}. Split into . In the first regime, both signs are always equal, and ; more precisely, letting be the angle , and ,
[TABLE]
In the opposite regime (), implies either (i) , which amounts to , ; (ii) or , e^{*}_{|\Gamma|}(\boldsymbol{p}+\boldsymbol{q})\gtrsim\frac{1}{m^{*}}\Big{(}2p^{*}_{F}|\boldsymbol{q}|\cos(\theta)+|\boldsymbol{q}|^{2}\Big{)}\gtrsim e^{*}_{|\Gamma|}(\boldsymbol{p}) (briefly said, this is a sub-regime where dominates w.r. to ).
The small integration volume in in regime (i) makes the integral restricted to convergent,
[TABLE]
In sub-regime (ii), letting ,
[TABLE]
too.
Proof of Lemma 1.3.
We must evaluate the second derivatives at zero momentum of the Cooper pair bubble diagram of §1.4, \partial_{q^{0}}^{2}{\cal A}_{q}^{j\to}(\Upsilon_{3})\Big{|}_{q=0} and \boldsymbol{\nabla}^{2}{\cal A}_{q}^{j\to}(\Upsilon_{3})\Big{|}_{q=0}.
This is a tedious task in Fourier coordinates. It is easier to use an argument in direct space, where the positivity argument (1.67) can be used. First
[TABLE]
[TABLE]
We now need an estimate of . The covariance function is decomposed in [25] into anisotropic angular sectors of size along the Fermi circle. This makes it easy to understand where the major contribution to comes. Decompose into , where the are anisotropic as indicated, but otherwise similar to the of section 2. Take , where is an orthonormal basis, with following the direction of some sector , and . Let be the index of some angular sector, and the projection onto the Fermi sphere of some momentum contained in the given angular sector. Then (as proven in [25]), for every , there exists s.t.
[TABLE]
(compare with Proposition 2.3) where is the projection of along , and its projection along the orthogonal direction. So decreases quasi-exponentially outside a box of dimensions and volume . For , this means that the main contribution comes from sectors with and . Then
[TABLE]
where and (to leading order) has been replaced with . If , then the first-order Taylor expansion is a good approximation. By symmetry, the term of order [math] vanishes, so the maximum order of magnitude of is obtained not in a neighborhood of [math], but for . Replacing by , one obtains for such
[TABLE]
as expected from (5.25), whereas
[TABLE]
Hence (integrating over a box and summing over the angular sectors ) one may conjecture that
[TABLE]
Letting be near the transition scale, this yields -\boldsymbol{\nabla}^{2}{\cal A}_{q}(\Upsilon_{3})\Big{|}_{q=0}\approx(\frac{p^{*}_{F}}{m^{*}})^{2}\frac{m^{*}}{\Gamma_{\phi}^{2}}\,{\mathrm{Id}}, which is what one wanted to prove. Considering instead \partial^{2}_{q^{0}}{\cal A}_{q}(\Upsilon_{3})\big{|}_{q=0}, one expands similarly , replaces by and obtains a bound of the same magnitude but without the dimensionful prefactor .
This is however not a rigorous proof, since there is some overlap between sectors, namely, e.g. in general for neighboring sectors , with . It is simpler – though less instructive – to integrate over the angular coordinate along the Fermi sphere; making computations for an arbitrary value of and will allow us to produce at the same time the asymptotics (0.74) for the one-point density kernel, and the desired order of magnitude for the second derivatives of . Let . We use the standard formula
[TABLE]
in terms of Bessel functions of the first kind, and neglect the ultra-violet cut-off at scale . Then the theorem of residues yields
[TABLE]
where . Replacing by , and letting (as a vestige of the smoothing due to ultra-violet cut-off) C^{j\to}_{1,1}(0,\boldsymbol{x}):={1\over 2}\lim_{\tau\to 0^{+}}\Big{(}C_{1,1}^{j\to}(\tau,\boldsymbol{x})+C_{1,1}^{j\to}(-\tau,\boldsymbol{x})\Big{)}, so that the coefficient is replaced by the even function in the above integrals, one obtains (0.74) in the limit by remarking that
[TABLE]
Assume on the other hand that . Hence one may (up to an exponentially small error) replace in the above expressions by , and by , with . Hence
[TABLE]
with , where is a modified Bessel function (exponentially decreasing at infinity for ), see [4] (7) p. 11. All together, one has found that is equal for to , in conformity with (5.27), times a fast oscillating phase function.
The above argument is easily adapted to prove that for . Namely, the change of variable () reduces the problem to the evaluation of the Laplace transform of the function , which is also equal to the function , see [4] (29) p. 136.
We further prove eq. (1.87) and (1.88). The infra-red cut-off bound (1.88) actually plays a minor rôle for the proof of (1.87), so let us prove it first. We concentrate on spatial gradients , with , Pre or , but the bounds we write down also hold for homogeneized gradients . Let . Then as shown by following the arguments in the proof of (1.59); ; and , from which by an easy induction (1.88) holds.
We now return to corrections to leading-order behavior of ,
[TABLE]
Error terms involve both the correction to leading-order term in (1.84) for , and a subleading contribution obtained by integrating over large transfer momenta such that . In either case, the idea is to apply repeated integrations by parts on or ,
[TABLE]
where is arbitrarily large and . The contribution of large transfer momenta is easily bounded using (1.88): if ,
[TABLE]
So assume that ; we leave it to the reader to check that
[TABLE]
for . Hence
[TABLE]
from which, finally, (1.87) holds.
5.5 Error terms for Ward identities
We bound here the error terms ”err.” due to the propagator cut-offs neglected in §3.3. We consider to begin with the variation of the fermionic covariance kernel , under an infinitesimal gauge transformation . The commutator is a sum of two terms, . Since acts multiplicatively on fields, its Fourier transform acts by convolution,
[TABLE]
Write the contribution of this term to (3.67), and proceed as done in the proof of Ward identities, see §3.3 B., namely, take the functional derivative , multiply by the low-momentum kernel , and integrate w.r. to . After a Fourier transform, this is equivalent to the expression
[TABLE]
The above integral is restricted to . By symmetry, the difference in (5.40) vanishes to order 1 when . Hence the kernel in (5.40) may be replaced by a second-order derivative bounded like . Thus the error term exhibits a squared spring term as in §3.3, and a further absolute small factor .
Let us illustrate this on an example.
••••q_{1}$$q_{4}$$q_{3}$$q_{2}$$p_{13}$$j_{+}$$k_{+}
Fig. 5.5.1. An example of multi-scale diagram.
Internal momenta are , connecting with external transfer momentum to with external transfer momentum . Let . The amplitude of the amputated diagram is given by an integral in ,
[TABLE]
Let Then
[TABLE]
A similar analysis may be done for the variation under of the Goldstone boson propagator . This time, in (5.42) may be replaced by a second-order derivative bounded like , yielding precisely the desired squared spring factor.
5.6 Mayer expansion
We introduce here a variant of the BKAR formula already described in §5.2, to be used for the Mayer expansion (see §2.3.5). The idea is to test overlaps between polymers in with a small number of external legs, so as to extract leading-order translation-invariant quantities.
Let be the set of polymers in . Among these polymers, there are polymers with external legs, making up a subset . The complementary set is made up of polymers with external legs, which need no particular further treatment. The following variant of BKAR’s formula, found originally in [1], is stated in the present form in [54]. We now denote by a pair of polymers connected by a link .
Proposition 5.4** (restricted 2-type cluster or BKAR2 formula)**
Assume . Choose as initial object an object of type 1, and stop the Brydges-Kennedy-Abdesselam-Rivasseau expansion as soon as a link to an object of type 2 has appeared. Then choose a new object of type , and so on. This leads to a restricted expansion, for which only the link variables , with , have been weakened. The following closed formula holds. Let be a link weakening of , and a smooth function. Let be the set of forests on , each component of which is (i) either a tree of objects of type , called unrooted tree; (ii)or a rooted tree such that only the root is of type . Then
[TABLE]
where is either [math] or the minimum of the -variables running along the unique path in from to , and is the forest obtained from by merging all roots of into a single vertex.
The above Proposition is applied to the non-overlap function =NonOverlap in §2.3.5.
Index of notations
[TABLE]
[TABLE]
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] A. Abdesselam, V. Rivasseau. Trees, forests and jungles: A botanical garden for cluster expansions , Lecture Notes in Physics 446 , Springer Berlin/Heidelberg (1995).
- 2[2] A. Altland, B. Simons. Condensed matter field theory , Cambridge University Press (2010).
- 3[3] J. Bardeen, L. N. Cooper, J. R. Schrieffer. Theory of superconductivity , Phys. Rev. 108 , 1175– (1957).
- 4[4] H. Bateman. Tables of Integral Transforms (Vol. I), Mc Graw-Hill Book Company (1954).
- 5[5] K. van der Beek, F. Lévy-Bertrand. Les vortex dans les supraconducteurs. Une matière molle au sein de la matière dure , Reflets de la Physique, Société française de physique, January 2012.
- 6[6] J. Bellissard, J. Magnen, V. Rivasseau. Supersymmetric analysis of a simplified two-dimensional Anderson model at small disorder , Markov Process. Related Fields 9 , 261–278 (2003).
- 7[7] N. N. Bogoliubov, On the theory of superfluidity , J. Phys. 11 , 23– (1947).
- 8[8] N. N. Bogoliubov. A new method in the theory of superconductivity I , Soviet Phys. JETP 34 , 41–46 (1958).
