Effective locality and chiral symmetry breaking in QCD
Thierry Grandou, Peter H. Tsang

TL;DR
This paper explores how effective locality, a non-perturbative property of QCD's fermionic Green's functions, relates to dynamical chiral symmetry breaking under certain approximations.
Contribution
It investigates the connection between effective locality and chiral symmetry breaking in QCD using eikonal and quenching approximations.
Findings
Effective locality is linked to chiral symmetry breaking in QCD.
The study provides insights into non-perturbative aspects of QCD.
Results are obtained within specific approximation schemes.
Abstract
A few years ago the use of standard functional manipulations was demonstrated to imply an unexpected property satisfied by the fermionic Green's functions of QCD, and called effective locality. This feature of QCD is non-perturbative as it results from a full integration of the gluonic degrees of freedom. In this paper, at eikonal and quenching approximation at least, the relation of effective locality to dynamical chiral symmetry breaking is examined.
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.
Effective locality and chiral symmetry breaking in
T. Grandou
Université Côte d’Azur,
INPHYNI UMR CNRS 7010; 1361 routes des Lucioles, 06560 Valbonne, France
P.H. Tsang
Physics Department, Brown University, Providence, RI 02912, USA
Abstract
A few years ago the use of standard functional manipulations was demonstrated to imply an unexpected property satisfied by the fermionic Green’s functions of QCD, and called effective locality. This feature of QCD is non-perturbative as it results from a full integration of the gluonic degrees of freedom. In this paper, at eikonal and quenching approximation at least, the relation of effective locality to dynamical chiral symmetry breaking is examined
Non-perturbative QCD, functional methods, random matrices, chiral symmetry breaking.
pacs:
12.38.Cy
I Introduction
The two major issues of QCD in its non-perturbative regime are confinement and chiral symmetry breaking. Some years ago a new non-perturbative property of QCD was disclosed and shown to be exact QCD1 ; QCD-II . It is a non-perturbative property in the sense that gluonic degrees of freedom are integrated out from the onset. Therefore, this new property, called effective locality, bears on fermionic Green’s functions and it can be summarised as follows.
For any fermionic -point Green’s functions and related amplitudes, the full gauge-fixed sum of cubic and quartic gluonic interactions, fermionic loops included, results in a local contact-type interaction. This local interaction is mediated by a tensorial field which is antisymmetric both in Lorentz and color indices. Moreover, the resulting sum appears to be fully gauge-fixing independent, that is, gauge-invariant.
Years ago, a similar property was discovered in an euclidean formulation of the pure Yang-Mills theory where it complied with a form of duality, established at least at first non-trivial orders of a semi-classical expansion Hug .
Recently, it has been suggested that effective locality could represent the very way non-abelian gauge invariance is realised in the non-perturbative regime of , avoiding the intractable and never ending issue of Gribov’s copies Thess . If this property is really relevant to , it should also shed some light on the two fundamental issues just alluded to above, confinement and dynamical chiral symmetry breaking. It is the latter which is the purpose of the current letter.
Calculations being quite involved, the chiral condensate will be estimated using the quenching approximation and some light form of eikonal approximation, admittedly well suited to high energy scattering processes. Present calculations will be carried on for the two-point Greens function, even though it cannot be conclusive in itself, due to a trivial algebraic statement. This obstruction which, besides its algebraic triviality, presents also meaningful aspects, is circumvented by a more involved four-point Green’s function calculation whose full details will be given elsewhere tgpt .
The paper is organised as follows. Calculations are set up in the next section, revisiting somewhat the effective locality context where they are carried out. The chiral condensate is defined as the limit of the Green’s function , and is evaluated in the case of one quark flavour. It is shown that the result is trivialised algebraically. As in the case of massive two-dimensional , but for different reasons (contrarily to the case of where it can be invoked FGHF , the cluster decomposition property doesn’t hold in ), this difficulty is circumvented by a four-point Green’s function calculation to which an element will be borrowed.
The analytically continued Random Matrix treatment of the matter involves a Vandermonde determinant comprising as much as monomials, alternate in sign, to be integrated upon. Even to most powerful computers, this seems to be an impossible task to carry out.
Now, relying on the properties of random matrice’s spectra, it is shown how to use Wigner’s semicircle law so as to frame the overall contribution of so many monomials and to offer an estimate of the chiral condensate.
These results are commented in a conclusion while some extra technical details are deferred to an Appendix in order to alleviate an already technical enough paper.
II Calculations
II.1 Two terms, one of them vanishing in the chiral limit
At quenching, the relevant expression to be coped with reads,
[TABLE]
where stands for a normalisation factor, taken at zero coupling , and where,
[TABLE]
in which is the usual Feynman propagator, while any other choice would yield the same results tg2 . The functional can be written with the help of a standard Schwinger/Fradkin representation Herb ,
[TABLE]
where and the subscript indicates Schwinger’s proper-time ordering. Ignoring the term involving the part in view of the eikonal approximation, allows one to deal with the simplified expression of,
[TABLE]
Inserting this expression into (1), it is necessary to take the -dependences outside of the ordered exponential so as to permit the functional operator of (1) to act upon linear and quadratic -field dependences. This can be achieved by introducing two extra functional integrations,
[TABLE]
Even within the approximations being used, the net expression for the order parameter is rather cumbersome. Written in extenso it reads as,
[TABLE]
followed, in the last line, by the prescription of . The prescription of trace bears on spinorial and color indices. Out of (II.1), two contributions are obtained,
[TABLE]
and,
[TABLE]
where in both expressions (II.1) and (II.1) the following definitions apply,
[TABLE]
with,
[TABLE]
For the first contribution (II.1) one gets for the last line,
[TABLE]
This is the known effective locality phenomenon of the resulting effective interaction QCD1 ; QCD-II where between currents (z) a local interaction is mediated by the term . In the strong coupling limit, , relevant to the non-perturbative regime of , the -piece of the current can be neglected in front of , and in view of (9) and (10), the effective interaction term reduces to,
[TABLE]
Theorems in Wiener functional space Lapidus , where Fradkin’s variables are taken to live, can be used to evaluate the product of distributions appearing in (12) and one gets QCD6 ,
[TABLE]
with a function to be identified to , out of a -point fermionic Green’s function calculation carried out at same eikonal approximation QCD6 ; tgpt . In this case, one has and thus , in some reference frame. The interaction is peaked at the value of , QCD1 ; QCD6 , and the remaining integrations can be carried through with the help of Random Matrices, as will be shown below for the second term (II.1). Being regular, these integrations lead to a result on the order of which vanishes in the chiral limit, . Details are given in the appendix. This is in agreement with a long known result established in QED, that scalar electrons/positrons do not produce chiral symmetry breaking Schwinger . The same applies here, where contrarily to (II.1), in (II.1) nothing refers to any spinorial structure and thus to any chiral symmetry possible breaking.
In other words, contribution (II.1) remains on the order of the explicit chiral symmetry breaking corresponding to non zero quark masses, and vanishes in the chiral limit where quarks masses are taken to zero.
For the second contribution (II.1), the last line can be evaluated in two different ways leading to the same result,
[TABLE]
whose evaluation yields,
[TABLE]
where the first term is to be part of the normalization constant . As before, the strong coupling limit is invoked to simplify (II.1) to,
[TABLE]
This is of the form,
[TABLE]
and in view of given in (10), involves the product of delta distributions,
[TABLE]
From Wiener space results again, the evaluation of (18) can be inferred from the -point fermionic Green’s function’s case QCD6 to read,
[TABLE]
where stands for the mass squared associated to the property of effective locality QCD1 ; QCD-II ; QCD5 ; QCD6 . Relying on the functional identities (57) and (59), it is straight forward to obtain for (II.1) the expression,
[TABLE]
The interaction term being peaked at , only the contribution at of will play a role, all of the other points giving an infinite continuous product of factors due, at each point, to a normalisation taken at , QCD1 ; QCD-II ; QCD5 ; QCD6 . The small increment is introduced so as to keep the exponent dimensionless, and will play no role otherwise (see comment after (II.2)). For the last line of (II.1) this gives,
[TABLE]
where , of dimension , will disappear from the final result (56). This dimensionful parameter can be used to turn the original -field into the dimensionless field, . As the normalisation factor entails both and , the way in which this elementary space-time cell extension comes into play is read out of the last line of (II.1),
[TABLE]
where the same symbol of is eventually written instead of .
To proceed with a(n analytically continued QCD6 ) Random Matrix integration of (22) it is convenient to define where the are the eight Lie algebra generators of taken in the adjoint representation. Matrices turn out to be real symmetric traceless matrices of format , that is , at space-time dimensions and colours, as . Matrices are also skew-symmetric.
Likewise one can define the 32-component vector and the matrices , where the are the eight Lie algebra generators of taken, this time, in the fundamental representation, i.e., the Gell-Mann matrices. Also we define and . As recalled in Appendix C, the linear spaces which are constructed in this way are canonically endowed with a metric structure. A one-to-one relation exists between new and original indices. With and , one has by construction,
[TABLE]
where the notation stands for -modulo-n. In this way (II.1) can be re-written in a unified covariant way as,
[TABLE]
where, as will be commented later, the dimensional quantity is identified out of a -point calculation QCD6 . The expression (II.1) represents the contributions of (II.1) and (II.1) to the chiral limit of at strong coupling, quenching and a mild eikonal approximations. The and of the integration measures and respectively, simplify to unity with the normalisation and leave an order parameter of dimension as it should.
Now, the matter is that one has,
[TABLE]
where the trace-prescription bears now on -indices, , a mix of colour (), and spin indices () in view of the definition of matrices . This algebraic statement completely trivialises the order parameter of the chiral symmetry, and it is worth observing that the situation is not improved by relaxing this part of the eikonal approximation. Including the -term of (II.1), in effect, amounts to redefine vectors s according to,
[TABLE]
whose components, with , are now matrix-valuated. This additional structure, however, does not improve on the trivial result (25) since one obtains again,
[TABLE]
with, again, the factor a -point calculation element QCD6 . Of course, this does not mean that the order parameter be trivial, but imposes to circumvent the trivial algebraic vanishing of (25).
As in massive at space-time dimensions FGHF , a trick allows to circumvent this issue and consists in calculating the -point fermionic Green’s function in the limits of , . In this way in effect, the trivial left hand side of (25) gets basically replaced by expressions like whose non-triviality can be checked at , and and/or .
Though straight forward, calculations are quite long and cumbersome and will be presented elsewhere in extenso as they appear to testify to non-perturbative peculiar features tgpt . As discussed in Section III, the necessity of proceeding through a -point fermionic Green’s function calculation in order to obtain the chiral condensate, will be put forth as a possibly meaningful aspect of the non-perturbative regime of .
Like in massive two-dimensional , but for different reasons FGHF ; tgpt , the -point fermionic Green’s function calculation displays the squared value of the chiral condensate. In this letter, accordingly, we will content ourselves of a calculation based on (II.1), and ignore that it is trivialised by the trace prescription. The difference will come out as some missing accurate dependences, to be corrected later tgpt , and which fall beyond the scope of our present focus, i.e., the relation of the enigmatic effective locality mass scale, , to the dynamical breaking of chiral symmetry.
Moreover, in the context of these calculations and beyond, that is for the sake of further experimental predictions PeterThesis , this will allow one to introduce a manipulation of Wigner’s semicircle law Mehta .
II.2 Using Wigner’s semicircle distribution
In the second line of (II.1), one has for the measure of integration Mehta ,
[TABLE]
where the s are the eigenvalues of and where . In term of this parametrisation of , (II.1) can be expressed as,
[TABLE]
where , a shorthand for the last factors of (28), is a Haar measure of integration on the orthogonal group , tg . In (II.2) the dimensions of s and of s have been taken out from the numerator and denominator , by defining the dimensionless -vectors and , while a subsidiary integration on is introduced both in (II.2) and in the normalisation factor , in order to implement the traceless property of -matrices.
Differentiating with respect to and integrating on the orthogonal group tg , one obtains,
[TABLE]
retaining the first non trivial contribution to the expansion of the interaction term in (II.2). Integrations on and then on are easily carried out and give (up to the dots),
[TABLE]
As shown in Appendix C, the presence of the metric tensor element is crucial to allow the replacement of by the constant ,
and to write for (II.2),
[TABLE]
This result however is vanishing as can be seen in various ways tgpt , the simplest being to rely on the structure of , the spectrum of the matrices defined in (27), displaying pairs of equal and opposite real eigenvalues tg . The expression (II.2) can accordingly be re-written as tg ,
[TABLE]
whose vanishing is manifest, but makes sense only if the same measure of integration on does not lead to a vanishing result for the normalisation constant either. This latter condition allowing one to introduce some basic elements of a Random Matrix treatment to be used at next order, , with defined in Equation (46) below, is analysed in some details now.
That the normalisation be non-zero can be seen as follows,
[TABLE]
where the prefactor stands for , the volume of the orthogonal group vol . After rescaling and analytical continuation QCD6 , , with , one obtains,
[TABLE]
The exponential suppression of any -contribution beyond is used now. First, one has clearly,
[TABLE]
because . Given the symmetry of the integrand one may write eventually, with the exponent ,
[TABLE]
Introducing now the standard Random Matrix definitions,
[TABLE]
[TABLE]
and,
[TABLE]
in accord to another standard definition of Random Matrix theory Mehta , one can write (37),
[TABLE]
where the last and crucial inequality is obvious. Just before it, the approximate equality means ‘up to sub-leading terms in the large -limit’. In (II.2) in effect, use has been made of Wigner’s semicircle law, valid at large -values Mehta ,
[TABLE]
and in this way the intractable sum of monomials generated by expanding the Vandermonde determinant of (38) is circumvented, while corrections to the asymptotics of (42) can be calculated in a systematic fashion A8 . This is the more appropriate as the universality of Wigner’s semi-circle law is now recognised to extend far beyond the realm of its original derivation Krajewski .
At next order (, and omitting the bars, i.e., , one needs the following averages tg ,
[TABLE]
which result from (II.2) as the derivation has been performed. Under integration on the right hand side of (43) reduces to tgpt ,
[TABLE]
and one has therefore to evaluate,
[TABLE]
having rescaled the eigenvalues according to as in (II.2). The dimensionless ratio has been introduced,
[TABLE]
which, even at large enough coupling constant, , remains a small number in the range of the -momenta which are relevant to this non-perturbative case (see Appendix E in tg ). The second line of (II.2) yields a factor of
[TABLE]
and in this paper, one will calculate the contribution due to the first term in the big parenthesis, of order , because the second one, of order , represents only a much smaller contribution. Then, by analytically continuing according to with as was done to pass from (34) to (35), one gets,
[TABLE]
This estimate of order satisfies,
[TABLE]
where, in (II.2), a change of integration variables has been performed, explicitly, , and where . Using now Random Matrix definitions (38), (40) and (42),
[TABLE]
In (II.2), the two last integrals can be decomposed into a sum , with,
[TABLE]
[TABLE]
and while as proven in the Appendix, one gets for the finite result,
[TABLE]
where is the probability integral GR . This is obtained by using formulae and of Ref.GR . The normalisation is evaluated with formula of GR and the same inequality as used for (II.2). Up to the same overall multiplicative factor of as the one in (II.2) with which it cancels out, one obtains,
[TABLE]
leading for to an estimation,
[TABLE]
Out of a -point calculation relevant to a 2-by-2 scattering quark process and in the same (mild) eikonal approximation, one can identify QCD6 , where and are the quark’s energies and momenta in their center of mass system, while the elementary space-time cell probed by effective locality has volume , QCD-II ; QCD6 ; tgpt ; tg . Now, by the same token, this means also that in the limit , cannot approach beyond this elementary space-time cell extension, that is, one must have essentially , and the chiral condensate reads eventually as,
[TABLE]
The absolute value of (56) represents a lower boundary of . As can be proven easily tgpt , an upper and finite boundary value can be found also for and matters to make sure that one deals with a definite quantity. This estimate of the chiral condensate comes out as a product of three groups of terms which can be given the following comments.
(i) First and most importantly, effective locality involves dynamical chiral symmetry breaking (through non normal-ordered Green’s function calculations jamin ; sjbetal ).
(ii) The chiral condensate goes like the third power of the effective locality mass scale , suggesting a close connection between effective locality and chiral symmetry breaking.
(iii) Then one observes something one could dub a partonic depleting function, here expressed in the center of mass system of a two quark scattering processs, i.e., identified out of a -point Green’s function QCD6 ; tgpt . One can see that the chiral condensate magnitude is modulated by a function of energy and momentum variables and attached to the quarks, taken as particles (partons). Would the latter variables be in a range of magnitude corresponding to the perturbative regime of QCD, i.e., say, then the particle-like character of quarks would allow to be replaced by , the squared quark mass.
What (56) suggests anyway is that as the energy and momentum of the scattering quarks increase, the chiral condensate gets substantially depleted in agreement with the fact that chiral symmetry breaking does not result from a perturbative mechanism.
(iv) Eventually comes an involved combination of numerical factors depending on and (). Now, approximations have been used as well as, in the end, Wigner’s semi-circle law which applies to very large values of . What is worth pointing out, however, is that using Wigner’s law allows one to frame the contribution attached to the intractable sum over the too large number of monomials coming from a Vandermonde determinant expansion, and that corrections to this law can be calculated in a systematic way A8 . Numerical estimates could be given elsewhere tgpt , as the refining of the integration measure on - matrices will allow better estimates of both (II.2) and (54).
III Conclusion
Remarks are in order concerning the result (56) and the four comments given above.
-
Effective locality, an exact non-perturbative property of QCD-II , involves dynamical chiral symmetry breaking. In this calculation this is established in the case of one flavour of quark, at large coupling, quenching and a partial eikonal approximation. While the latter is recognised relevant to the matter of high energy scattering, the quenching approximation shouldn’t be taken as a proviso to the result. In effect, as supported by models hftg , numerical simulations and other non-perturbative approaches such as the Dyakonov instanton calculus Dyakonov , whenever chiral symmetry breaking can be established, it is an outcome of both quenched and non-quenched calculations.
-
Since it was discovered QCD1 , effective locality comes about with a mass scale, , necessarily associated to it in order to avoid triviality QCD11 . The relation of non-perturbative phenomena to this mass scale is of utmost interest, and motivated the current analysis. The fermionic condensate, a measure of dynamical chiral symmetry breaking, comes out to be on the order of the third power of it, a reasonable answer in view of a pure dimensional argument. However, estimation (56) tells more.
-
The result’s magnitude is modulated by a decreasing function of and , the energy and momentum of scattering quarks. The more increase, that is, the more ‘partonic’ the situation becomes, the more the chiral condensate is depleted, down to zero, recovering the perturbative result. This is also in full agreement with experimental facts according to which non-perturbative aspects are washed out beyond an energy scale of about Kos .
In the formal context of effective locality, it is remarkable that there is no way to identify such a depleting function out of a -point Green’s function calculation, such as , and this, irrespectively of the ‘algebraic trivialisation’ of the result which is inherent to the -point case and is circumvented by a -point calculation (in the context of effective locality, in effect, a -point calculation only involves a self-interaction whose renormalisation is achieved perturbatively along the lines described in QCD5' and in agreement with a prediction of Light Front SJB ).
One can venture the idea that this fact, not innocuous, may be meaningful.
Effective locality shows up once gluonic degrees of freedom have been integrated out, and the resulting effective interaction between quarks is then given by -by- body local interactions tgD . This is clearly reminiscent of the reduction of low-energy to Quark Flavour Dynamics HugoR ; Schaden , itself reminiscent of an extended Nambu-Jona Lasinio model NJL . Thus, if this is the relevant quark dynamical structure in the non-perturbative regime of , it should not come as a surprise that dynamical chiral symmetry breaking be revealed by a -point fermionic Green’s function calculation.
In some sense, it is as if effective locality did not favour much the long held point of view relating dynamical chiral symmetry breaking to the quark dressed propagator, in particular to its massive pole sjbetal (in this respect it is worth observing also that in a most recent axiomatic analysis, the definite existence of such massive poles has been recognised a questionable, gauge-dependent issue PL ). Quite on the contrary, this very aspect of dynamical chiral symmetry breaking would appear to be consistent with the manner confinement seems to show up in this same effective locality context, as could be argued elsewhere QCD11 .
- So far, out of effective locality, some qualitative features could be derived from the general form of fermionic -point Green’s functions QCD-II ; QCD6 ; tg ; tg2 , whereas phenomenological applications, however promising, required further simplifications QCD5 ; QCD5' . A major issue has been the very big number of monomials comprised in the Vandermonde determinant (28), alternate in sign, whose global behaviour when summed up is difficult to appreciate. The present calculation, though a bit technical, avoids this issue by relying on Wigner’s semi-circle law. In principle, the latter applies to the cases of very large -values while the current value of is only . It would be interesting to know how fast, as a function of , does the semi-circle law ‘saturates’, as this law would matter to improve on possible phenomenological applications. Though a bit lengthy, it is to be noted that corrections to the asymptotics of Wigner’s law can also be dealt with analytically. This and other related improvements will be coped with in a future publication, when some approximations will be relaxed tgpt .
For our present concern, relying on Wigner’s law allows one to prove the important fact that dynamical chiral symmetry breaking is an outcome of effective locality, at least within a set of approximations known to preserve chiral symmetry breaking if any; and accordingly, to continue proposing effective locality as a possibly interesting access to the non-perturbative regime of .
IV appendix
IV.1 On the chiral limit of Equation (II.1)
In Section II, it is stated that in the chiral limit of , no contribution is generated by the term (II.1). This can be seen along the following steps. One can rely on a functional identity,
[TABLE]
where , and a representation of the distribution,
[TABLE]
This allows to take advantage of another functional identity,
[TABLE]
and to write the first contribution (II.1) as,
[TABLE]
where the Fradkin’s field variables are not treated on the same footing as their derivatives, , which is consistent with the assumption of taking them to live in a Wiener functional space. Now as argued in the main text, the interaction is peaked at and (IV.1) is reduced to,
[TABLE]
As proven in the main text the integration on the -fields is regular. In the limit , the first line is divergent and is taken care of by renormalisation. Therefore, as the quark mass parameter is finally taken to zero, the contribution of (IV.1) and thus of (II.1) vanishes.
IV.2 On the vanishing of
Recalling ,
[TABLE]
integration on can be carried out with the help of formula in GR . One gets,
[TABLE]
whose subsequent integration on is vanishing.
IV.3 On Equation (II.2)
In the construction outlined right after (22), which aims also at dealing with the required calculations in some unified and covariant way, the tensorial product space of vectors has dimension , and from the metrics on the composing spaces, is endowed with a metric tensor g_{ii}\equiv\bigr{[}(g_{\mu\nu})\otimes(\delta_{ab})\bigl{]}_{ii} which inherits the Minkowskian character of spacetime, . That is, it is a diagonal metric tensor with components and components . This is at the origin of a point which would not come about with a pure Euclidean metric tensor but requires some care in the current pseudo-euclidean case.
In (II.2), it is easy to check that the integration,
[TABLE]
does depend on the sign of . This dependence, erased within the subsequent integration on , pops out again in the end, to wit,
[TABLE]
which of course is zero in view of (25). In the -point calculation designed to circumvent this algebraic triviality, what comes about, as announced after (27), is the trace whose only cases of non-triviality require . Now, in this way one recovers independence on the sign of since for all . As a result, once summed upon , this allows to replace the vanishing by , as is claimed after (II.2).
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1(1) H.M. Fried, Y. Gabellini, T. Grandou, and Y.-M. Sheu, Eur. Phys. J. C 65 , 395 (2010).
- 2(2) H.M. Fried, T. Grandou, and Y.-M. Sheu., Ann. Phys. 327 , 2666 (2012).
- 3(3) H. Reinhardt, K. Langfeld and L. v. Smekal, Phys. Lett. B, 300 , 11 (1993); H. Reinhardt, "Dual description of QCD", ar Xiv:hep-th/9608191 v 1 (1996).
- 4(4) See A. Maas talk in XI Ith Quark Confinement and the Hadron Spectrum , Thessaloniki, Greece, 28 August-02 September 2016.
- 5(5) T. Grandou and P.H. Tsang, work in completion.
- 6(6) F. Guerin and H.M. Fried, Phys.Rev.D 33 (1986) 3039.
- 7(7) H.M. Fried, T.Grandou and R. Hofmann, Mod. Phys. Lett. A, Vol.32(2017) 1730030 (16 pages) (2017); ar Xiv:1706.02264 v 1 [hep-th].
- 8(8) H.M. Fried, Functional Methods and Models in Quantum Field Theory , (the MIT Press, Cambridge, Massachusetts, 1972). Basics of Functional Methods and Eikonal Models (Editions Frontières, Gif-sur-Yvette, France,1990).
