Robert Coquereaux
email: robert.coquereaux@gmail.com

Jean-Bernard Zuber
email: jean-bernard.zuber@upmc.fr

From orbital measures to Littlewood-Richardson coefficients and hive polytopes

Keywords: Horn problem. Honeycombs. Polytopes. SU(n) Littlewood-Richardson coefficients Mathematics Subject Classification 2010: 17B08, 17B10, 22E46, 43A75, 52Bxx

The volume of the hive polytope (or polytope of honeycombs) associated with a Littlewood-Richardson coefficient of SU(n), or with a given admissible triple of highest weights, is expressed, in the generic case, in terms of the Fourier transform of a convolution product of orbital measures. Several properties of this function -a function of three non-necessarily integral weights or of three multiplets of real eigenvalues for the associated Horn problem-are already known. In the integral case it can be thought of as a semi-classical approximation of Littlewood-Richardson coefficients. We prove that it may be expressed as a local average of a finite number of such coefficients. We also relate this function to the Littlewood-Richardson polynomials (stretching polynomials) i.e., to the Ehrhart polynomials of the relevant hive polytopes. Several SU(n) examples, for n " 2, 3, . . . , 6, are explicitly worked out.

Introduction

In a previous paper [START_REF] Zuber | Horn's problem and Harish-Chandra's integrals. Probability distribution functions[END_REF], the following classical Horn's problem was addressed. For two n by n Hermitian matrices A and B independently and uniformly distributed on their respective unitary coadjoint orbits O α and O β , labelled by their eigenvalues α and β, call ppγ|α, βq the probability distribution function (PDF) of the eigenvalues γ of their sum C " A`B. With no loss of generality, we assume throughout this paper that these eigenvalues are ordered,

α 1 ě α 2 ě ¨¨¨ě α n (1)
and likewise for β and γ. In plain (probabilistic) terms, p describes the conditional probability of γ, given α and β. The general expression of p was given in [START_REF] Zuber | Horn's problem and Harish-Chandra's integrals. Probability distribution functions[END_REF] in terms of orbital integrals and computed explicitly for low values of n.

The aim of the present paper is to study the relations between this function p, and the tensor product multiplicities for irreducible representations (irreps) of the Lie groups Upnq or SUpnq, encoded by the Littlewood-Richardson (LR) coefficients.

Our main results are the following. A central role is played by a function J n pα, β; γq proportional to p, times a ratio of Vandermonde determinants, see [START_REF] Coquereaux | On some properties of SU(3) Fusion Coefficients[END_REF]. This J n is identified with the volume of the hive polytope (also called polytope of honeycombs) associated with the triple pα, β; γq, see Proposition 4. It is thus known [START_REF] Heckman | Projections of Orbits and Asymptotic Behavior of Multiplicities for Compact Lie Groups[END_REF] to provide the asymptotic behavior of LR coefficients, for large weights. We find a relation between J n and a sum of LR coefficients over a local, finite, n-dependent, set of weights, which holds true irrespective of the asymptotic limit, see Theorem 1. In particular for SU(3), the sum is trivial and enables one to express the LR coefficient as a piecewise linear function of the weights, see Proposition 5 and Corollary 1. Implications on the stretching polynomial (sometimes called Littlewood-Richardson polynomial) and its coefficients are then investigated.

The content of this paper is as follows. In sec. 1, we recall some basic facts on the geometric setting and on tensor and hive polytopes. We also collect formulae and results obtained in [START_REF] Zuber | Horn's problem and Harish-Chandra's integrals. Probability distribution functions[END_REF] on the function J n . Section 2 is devoted to the connection between Harish-Chandra's orbital integrals and SUpnq character formulae, to its implication on the relation between J n and LR coefficients (Theorem 1), and to consequences of the latter. In sec. 3, we reexamine the interpretation of J n as the volume of the hive polytope in the generic case (Proposition 4), through the analysis of the asymptotic regime. In the last section (examples), we take n " 2, 3 . . . , 6, consider for each case the expression obtained for J n , give the local relation existing between the latter and LR coefficients (this involves two polynomials, that we call R n and p R n , expressed as characters of SUpnq), and study the corresponding stretching polynomials. Some of the features studied in the main body of this article are finally illustrated in the last subsection where we consider a few specific hive polytopes.

Convolution of orbital measures, density function and polytopes 1.Underlying geometrical picture

We consider a particular Gelfand pair pUpnq˙H n , Upnqq associated with the group action of the Lie group Upnq on the vector space of n by n Hermitian matrices. This geometrical setup allows one to develop a kind of harmonic analysis where "points" are replaced by coadjoint orbits of Upnq : the Dirac measure (delta function at the point a) is replaced by an orbital measure whose definition will be recalled below, and its Fourier transform, here an orbital transform, is given by the so-called Harish-Chandra orbital function. This theory of integral transforms can also be considered as a generalization of the usual Radon spherical transform (also called Funk transform). Contrarily to Dirac measures, orbital measures are not discrete, since their supports are orbits of the chosen Lie group. Such a measure is described by a probability density function (PDF), which is its Radon-Nikodym derivative with respect to the Lebesgue measure. In Fourier theory one may consider the measure formally defined as a convolution product of Dirac masses: ă δ a ‹ δ b , f ą" ş δ a`b pxqf pxqdx. Here we shall consider, instead, the convolution product of two orbital measures described by the orbital analog of δ a`b pcq, a probability density function labelled by three Upnq orbits of H n . These orbits and that function p may be considered as functions of three hermitian matrices (we shall write it ppC|A, Bq), and this answers a natural question in the context of the classical Horn problem, as mentioned above in the Introduction, see also sec. 1.1.4 below. This was spelled out in paper [START_REF] Zuber | Horn's problem and Harish-Chandra's integrals. Probability distribution functions[END_REF]. Our main concern, here, is the study of the relations that exist between this function p, and the tensor product multiplicities for irreducible representations (irreps) of the Lie groups Upnq or SUpnq, encoded by the Littlewood-Richardson (LR) coefficients N ν λµ . For small values of n the function p can be explicitly calculated; for integral values of its arguments, the related function J n can be considered as a semi-classical approximation of the LR coefficients.

Orbital measures

For F , a function on the space of orbits, and O A , the orbit going through A P H n , one could formally consider the "delta function" ă δ O A , F ą" F pO A q, but we shall use test functions defined on H n instead.

The orbital measure m A , that plays the role of δ O A , is therefore defined, for any continuous function f on H n , by

ă m A , f ą" ż Upnq f pu ‹ Auqdu
where the integral is taken with respect to the Haar mesure1 on Upnq, i.e., by averaging the function f on a Upnq coadjoint orbit.

Fourier transform of orbital measures

Despite the appearance of the Haar measure on the group Upnq entering the definition of m A , one should notice that this is a measure on the vector space H n , an abelian group. Being an analog of the Dirac measure, its orbital transform2 is a complex-valued function y m A pXq on H n defined by evaluating m A on the following exponential function: Y P H n Þ Ñ exppi tr pX Y qq P C. Hence we obtain :

y m A pXq " ż Upnq exppi tr pXu ‹ Auqqdu
As this quantity only depends on the respective eigenvalues of X and A, i.e., on the diagonal matrices x " px 1 , x 2 , . . . x n q, and α " pα 1 , α 2 , . . . , α n q, it is then standard to rename the previous Fourier transform and consider the following two-variable function, called the Harish-Chandra orbital function:

Hpα, i xq " ż Upnq exppi tr pxu ‹ αuqq du (2)

The HCIZ integral

The following explicit expression of H was found in [START_REF] Harish-Chandra | Differential Operators on a Semisimple Algebra[END_REF][START_REF] Itzykson | The planar approximation II[END_REF].

Hpα, i xq " sf pn ´1q pdet e i x i α j q 1ďi,jďn ∆pixq∆pαq

where ∆pxq " Π iăj px i ´xj q is the Vandermonde determinant of the x's.

Here and in the following we make use of the superfactorial sf pmq :"

m ź p"1 p! . (4
ă m A ‹ m B , f ą"ă m A b m B , pf q ą
where pf qpa, bq :" f pa `bq .

This orbital analog of δ a`b pcq has a non discrete support: for A, B P H n , the support of m A,B " m A ‹ m B is the set of uAu ‹ `vBv ‹ for u, v P Upnq. The probability density function p of m A,B is obtained by applying an inverse Fourier transformation to the product of Fourier transforms (calculated using y m A pXq) of the two measures:

ppγ|α, βq " 1 p2πq n ˆ∆pγq sf pnq ˙2 ż R n d n x ∆pxq 2 Hpα, i xqHpβ, i xqHpγ, i xq ‹ . (5)
Notice that p involves three copies of the HCIZ integral and that we wrote it as in integral on R n , whence the prefactor coming from the Jacobian of the change of variables. We shall see below (formulae extracted from [START_REF] Zuber | Horn's problem and Harish-Chandra's integrals. Probability distribution functions[END_REF]) how to obtain quite explicit formulae for this expression.

On polytopes

In the present context of orbit sums and representation theory, one encounters two kinds of polytopes, not to be confused with one another. On the one hand, given two multiplets α and β, ordered as in [START_REF] Beck | Coefficients and Roots of Ehrhart Polynomials[END_REF], we have what may be called the Horn polytope r H αβ , which is the convex hull of all possible ordered γ's that appear in the sum of the two orbits O α and O β . As proved by Knutson and Tao [20] that Horn polytope is identical to the convex set of real solutions to Horn's inequalities, including the inequalities (1), applied to γ. For SUpnq, this Horn polytope is pn ´1q-dimensional.

On the other hand, combinatorial models associate to such a triple pα, β; γq, with γ P r H αβ , a family of graphical objects that we call generically pictographs. This family depends on a number pn ´1qpn ´2q{2 of real parameters, subject to linear inequalities, thus defining a d-dimensional polytope r H γ α β , with d ď pn ´1qpn ´2q{2. These two types of polytopes are particularly useful in the discussion of highest weight representations of SUpnq and their tensor product decompositions.

Given two highest weight representations V λ and V µ of SUpnq, we look at the decomposition into irreps of V λ b V µ , or of λ b µ, in short, see below sec. 2. Consider a particular space of intertwiners (equivariant morphisms) associated with a certain "branching", i.e., a particular term ν in that decomposition, that we call an admissible triple pλ, µ; νq, see below Definition 3. Such ν's lie in the tensor polytope H λµ inside the weight space. The multiplicity N ν λµ of ν in the tensor product λ b µ is the dimension of the space of intertwiners determined by the admissible triple pλ, µ; νq.

As proved in [START_REF] Knutson | The honeycomb model of GL n pCq tensor products I: proof of the saturation conjecture[END_REF], is is also the number of pictographs with integral parameters. It is thus also the number of integral points in the second polytope that we now denote H ν λ µ . These integral points may be conveniently thought of as describing the different "couplings" of the three chosen irreducible representations.

Pictographs are of several kinds. All of them have three "sides" but one may distinguish two families: first we have those pictographs with sides labelled by integer partitions (KT-honeycombs [START_REF] Knutson | The honeycomb model of GL n pCq tensor products I: proof of the saturation conjecture[END_REF], KT-hives [START_REF] Knutson | The honeycomb model of GL n pCq tensor products II: Puzzles determine facets of the Littlewood-Richardson cone[END_REF]), then we have those pictographs with sides labelled by highest weight components of the chosen irreps (BZ-triangles [START_REF] Berenstein | Triple multiplicities for sl(r `1) and the spectrum of the external algebra in the adjoint representation[END_REF], O-blades [25], isometric honeycombs3 ). For convenience, we refer to H ν λ µ as the "hive polytope", or also "the polytope of honeycombs". As mentioned above, for SUpnq, and for an admissible triple pλ, µ; νq, the dimension of the hive polytope is pn ´1qpn ´2q{2: this may be taken as a definition of a "generic triple", but see below Lemma 1 for a more precise characterization. The cartesian equations for the boundary hyperplanes have integral coefficients, the hive polytope is therefore a rational polytope. All the hive polytopes that we consider in this article are "integral hive polytopes" in the terminology of [START_REF] King | Stretched Littlewood-Richardson and Kostka coefficients[END_REF], however the corners of all such polytopes (usually called "vertices") are not always integral points, therefore an "integral hive polytope" is not necessarily an integral polytope in the usual sense: the convex hull of its integral points is itself a polytope, but there are cases where the latter is strictly included in the former. We shall see an example of this situation in sec. 4.4.2.

We shall return later to these polytopes and to the counting functions of their integral points, in relation with stretched Littlewood-Richardson coefficients, see sec. 3.

1.3 Some formulae and results from paper [START_REF] Zuber | Horn's problem and Harish-Chandra's integrals. Probability distribution functions[END_REF] 1.3.1 Determination of the density p and of the kernel function J n Some general expressions for the three variable function p were obtained in [START_REF] Zuber | Horn's problem and Harish-Chandra's integrals. Probability distribution functions[END_REF]. For the convenience of the reader, we repeat them here.

The determinant entering the HCIZ integral is written as

det e i x i α j " e i 1 n ř n j"1 x j ř n k"1 α k det e i px i ´1 n ř x k qα j (6
)

" e i 1 n ř n j"1 x j ř n k"1 α k ÿ P PSn ε P n´1 ź j"1 e i px j ´xj`1 qp ř j k"1 α P pkq ´j n ř n k"1 α k q , (7)
where ε P is the signature of permutation P .

In the product of the three determinants entering (5), the prefactor e i ř n j"1 x j ř n k"1 pα k `βk ´γk q{n yields, upon integration over 1 n ř x j , 2π times a Dirac delta of ř k pα k `βk ´γk q, expressing the conservation of the trace in Horn's problem. One is left with an expression involving an integration over pn ´1q variables u j :" x j ´xj`1 .

ppγ|α, βq "

sf pn ´1q n! δp

ÿ k pα k `βk ´γk qq ∆pγq ∆pαq∆pβq J n pα, β; γq (8)
J n pα, β; γq " i ´npn´1q{2 2 n´1 n! π n´1 ÿ P,P 1 ,P 2 PSn

ε P ε P 1 ε P 2 ż R n´1 d n´1 u r ∆puq n´1 ź j"1
e i u j A j pP,P 1 ,P 2 q (9)

A j pP, P 1 , P 2 q " j ÿ k"1

pα P pkq `βP 1 pkq ´γP 2 pkq q ´j n n ÿ k"1 pα k `βk ´γk q, (10)
where the Vandermonde ∆pxq has been rewritten as r ∆puq :" ź 1ďiďjďn´1 pu i `ui`1 `¨¨¨u j q .

(11)

Discussion

Several properties of ppγ|α, βq and of J n are described in the paper [START_REF] Zuber | Horn's problem and Harish-Chandra's integrals. Probability distribution functions[END_REF]. We only summarize here the information that will be relevant for our discussion relating these functions to the Littlewood-Richardson multiplicity problem. Note that the above expression of A j is invariant under simultaneous translations of all γ's

@ i γ i Ñ γ i `c c P R .
In the original Horn problem, this reflects the fact that the PDF ppγ|α, βq of eigenvalues of C " A `B is the same as that of C `cI, with a shifted support. Therefore in the computation of J n pα, β; γq, one has a freedom in the choice of a "gauge" (a) either γ n " 0, (b) or γ such that ÿ

i γ i " ÿ i pα i `βi q , (12)
(c) or any other choice, provided one takes into account the second term in the rhs of [START_REF] Fulton | Representation Theory, A First Course[END_REF] (which vanishes in case (b)). Note also that enforcing [START_REF] Haase | Lecture Notes on Lattice Polytopes, Fall School on Polyhedral Combinatorics[END_REF] starting from an arbitrary γ implies to translate γ Ñ γ " γ `c, with c " 1 n p ř i α i `ři β i ´ři γi q. If the original γ has integral components, this is generally not the case for the final γ.

J n pα, β; γq has the following properties that will be used below: -(i) As apparent on [START_REF] Derksen | On the Littelwood-Richardson polynomials[END_REF], it is an antisymmetric function of α, β or γ under the action of the Weyl group of SUpnq (the symmetric group S n ). As already said, we choose throughout this paper the ordering (1) and likewise for β and γ. For pα, β; γq satisfying ( 12) -(ii) J n pα, β; γq is piecewise polynomial, homogeneous of degree 1 2 pn ´1qpn ´2q in α, β, γ in the generic case; -(iii) as a function of γ, it is of class C n´3 . This follows by the Riemann-Lebesgue theorem from the decay at large u of the integrand in [START_REF] Derksen | On the Littelwood-Richardson polynomials[END_REF], see [START_REF] Zuber | Horn's problem and Harish-Chandra's integrals. Probability distribution functions[END_REF]; -(iv) it is non negative inside the polytope r H αβ , cf sec. 1. By a slight abuse of language, when dealing with triples of highest weights pλ, µ; νq, we say that such an admissible triple is generic iff the associated triple pα, β; γq is, see below sec. 2.1. By another abuse of language, we also refer to a single highest weight λ as generic iff none of its Dynkin indices vanishes, i.e., iff λ does not lie on one of the walls of the dominant Weyl chamber, or if equivalently the associated α has no pair of equal components.

From its interpretation as a probability density (up to positive factors), it is clear that J n could vanish at most on subsets of measure zero inside the Horn (or tensor) polytope. Actually it does not vanish besides the cases mentioned in points (v-vii) of the previous list.

We want to construct the linear span of honeycombs r H γ αβ defined above in sect. 1.2. We first consider what may be called the "SUpnq case", where α n " β n " 0 and γ n is fixed by [START_REF] Haase | Lecture Notes on Lattice Polytopes, Fall School on Polyhedral Combinatorics[END_REF]. By relaxing the inequalities on the pn ´1qpn ´2q{2 parameters defining the usual honeycombs, one builds a vector space of dimension 1 2 pn ´1qpn `4q " 3pn ´1q `pn ´1qpn ´2q{2 whose elements are sometimes called real honeycombs. One may construct a basis of "fundamental honeycombs", see [START_REF] Fulton | Representation Theory, A First Course[END_REF], and consider arbitrary linear combinations, with real coefficients, of these basis vectors. The components of any admissible triplepα, β, γq, depend linearly of the components of the associated honeycombs along the chosen basis. In such a way, one obtains a surjective linear map, from the vector space of real honeycombs, to the vector space R 3pn´1q .

One sees immediately that its fibers are affine spaces of dimension d max " pn ´1qpn ´2q{2, and for fixed α, β they are indexed by γ, i.e., by points of R pn´1q . By taking into account the inequalities defining usual honeycombs, but still working with real coefficients, the fibers of this map restrict to compact polytopes whose affine dimension d is at most equal to d max (the dimension can be smaller, because of the inequalities that define bounding hyperplanes). For given α and β, if γ belongs to the Horn polytope r H αβ Ă R pn´1q , the corresponding restricted fiber is nothing else than the associated hive polytope r H γ αβ . We therefore obtain a map π whose target set is the Horn polytope, a convex set, and whose fibers are compact polytopes. We then make use of the following result5 : the dimension of the fibers of π is constant on the interiors of the faces of its target set. In particular, it is constant on the interior of its face of codimension 0, which is the interior of the Horn polytope r H αβ . In the present situation this tells us that the dimension of π ´1pγq " r H γ αβ which is the fiber above γ, is constant when γ belongs to the interior of the Horn polytope r H αβ . In particular, its d-dimensional volume, where d has its maximal value d " pn ´1qpn ´2q{2 for SUpnq, cannot vanish there. We shall see later (in section 3) that this volume is given by J n pα, β; γq.

In the case of GLpnq, (with α n , β n non fixed to 0), the argument is similar, so we have:

Lemma 1. For α and β with distinct components, the function J n pα, β; γq does not vanish for γ inside the polytope r H αβ .

2 From Horn to Littlewood-Richardson and from orbital transforms to characters

Young partitions and highest weights

An irreducible polynomial representation of GLpnq or an irrep of SUpnq, denoted V λ , is characterized by its highest weight λ. One may use alternative notations, describing this highest weight either by its Dynkin indices (components in a basis of fundamental weights) λ i , i " 1, ¨¨¨, n, and λ n " 0 in SUpnq ; or by its Young components, i.e., the lengths of rows of the corresponding Young diagram: α " pλq, i.e., i pλq "

n ÿ j"i λ j i " 1, ¨¨¨, n . (14)
Note that such an α " pλq satisfies the ordering condition (1).

In the decomposition into irreps of the tensor product of two such irreps V λ and V µ of GLpnq, we denote by N ν λµ the Littlewood-Richardson (LR) multiplicity of V ν . As recalled above, N ν λµ equals the number of honeycombs with integral labels and boundary conditions α " pλq, β " pµq, γ " pνq, i.e., the number of integral points in the polytope H ν λ µ [START_REF] Knutson | The honeycomb model of GL n pCq tensor products I: proof of the saturation conjecture[END_REF]. Given three U(n) (resp. SUpnq) weights λ, µ, ν, for instance described by their n (resp. n ´1) components along the basis of fundamental weights, invariance under the U(1) center of U(n) (resp. the Z n center of SUpnq), tells us that a necessary condition for the non-vanishing of N ν λµ is ř n j"1 jpλ j `µj ´νj q " 0 (resp. ř n´1 j"1 jpλ j `µj ´νj q " 0 mod n). Given three SUpnq weights λ, µ, ν obeying the above SUpnq condition, one can build three U(n) weights (still denoted λ, µ, ν) obeying the U(n) condition by setting λ n " µ n " 0 and ν n "

1 n ř n´1
j"1 jpλ j `µj ´νj q; in terms of partitions, with α " pλq, β " pµq and γ " pνq, the obtained triple pα, β; γq automatically obeys eq. ( 12). More generally we shall refer to a U(n) triple such that the equivalent U(n) conditions eq. ( 12), or eq. ( 15) below, hold true, as a U(n)-compatible triple, or a compatible triple, for short. Definition 2. A triple pλ, µ; νq of Upnq weights is said to be compatible iff

n ÿ k"1 kpλ k `µk ´νk q " 0 . (15)
For triples of SUpnq weights, we could use the same terminology, weakening the above condition (15) since it is then only assumed to hold modulo n, but in the following we shall always extend such SUpnq-compatible triples to U(n)-compatible triples, as was explained previously. We also recall another more traditional definition Definition 3. A triple pλ, µ; νq of Upnq or SUpnq weights is said to be admissible iff N ν λµ ‰ 0 . The reader should remember (at least in the context of this article !) the difference between compatibility and admissibility, the former being obviously a necessary condition for the latter.

For given λ and µ, or equivalently, given α and β, if N ν λµ ‰ 0 for some h.w. ν, the corresponding γ must lie inside or on the boundary of the Horn polytope r H αβ , by definition of the latter. Since for n ě 3 the function J n pα, β; γq is continuous and vanishes on the boundary of its support, evaluating it for α, β, γ does not provide a strong enough criterion to identify admissible triples pα, β; γq.

Relation between Weyl's character formula and the HCIZ integral

There is an obvious similarity between the general form (5) of the PDF ppγ|α, βq and the expression of the LR multiplicity N ν λµ as the integral of the product of characters χ λ χ µ χ ν over the unitary group SU(n) or over its Cartan torus

T n " Up1q n´1 N ν λµ " ż SUpnq du χ λ puqχ µ puqχ ν puq or N ν λµ " ż Tn dT χ λ pT qχ µ pT qχ ν pT q (16)
with the normalized Haar measure on T n ,

dT " 1 p2πq n´1 n! |∆pe i t q| 2 n´1 ź i"1 dt i , (17) for
T " diag pe i t j q j"1,¨¨¨,n with

n ÿ j"1 t j " 0 . (18)
This similarity finds its root in the Kirillov [START_REF] Kirillov | Lectures on the Orbit Method[END_REF] formula expressing χ λ as the orbital function H relative to O pλ`ρq , defined in [START_REF] Bégin | sup3q k fusion coefficients[END_REF], see below [START_REF] Knutson | The honeycomb model of GL n pCq tensor products II: Puzzles determine facets of the Littlewood-Richardson cone[END_REF][START_REF] Macdonald | The volume of a compact Lie group[END_REF]; note the shift of λ by the Weyl vector ρ, the half-sum of positive roots.

Recall Weyl's formula for the dimension of the vector space V λ of h.w. λ dim V λ " ∆pα 1 q sf pn ´1q with α 1 " pλ `ρq , and as defined in ( 14) .

From a geometrical point of view, this formula expresses dim V λ as the volume of a group orbit normalized by the volume of SUpnq, the latter being also equal to sf pn ´1q, once a natural Haar measure has been chosen, see [START_REF] Macdonald | The volume of a compact Lie group[END_REF].

From group characters to Harish-Chandra orbital functions

Kirillov's formula [START_REF] Kirillov | Lectures on the Orbit Method[END_REF] relates Weyl's SUpnq character formula with the orbital function of O α 1 . Here and below, the prime on α 1 refers to the value of α, for the shifted highest weights λ `ρ

α 1 " pλ `ρq , (20)
and likewise for β 1 , γ 1 . Indeed evaluated on an element T of the SUpnq Cartan torus as in ( 18), Weyl's character formula reads χ λ pT q :" tr V λ pT q " det e i t i α 1 j ∆pe i t q with ∆pe i t q "

ź 1ďiăjďn pe i t i ´ei t j q , (21)
or in terms of the orbital function H defined in (2) and made explicit in (3)

χ λ pT q " ∆pα 1 q sf pn ´1q ˜ź 1ďiăjďn i pt i ´tj q pe i t i ´ei t j q ¸Hpα 1 , i tq (22)
or, owing to the Weyl dimension formula ( 19)

χ λ pT q dim V λ " ∆pi tq ∆pe i t q Hpα 1 , i tq . (23)

The polynomial R n pT q

Consider the following (semi-convergent) integral

J " ż R du e i uA u A P R
a one-dimensional analogue of the integral encountered in [START_REF] Derksen | On the Littelwood-Richardson polynomials[END_REF]. If A is a half-integer, we may write according to a well-known identity. If A is an integer, the previous sum over n is understood as a principal value. Then

A half-integer J " ż π ´π du e i uA 8 ÿ n"´8 p´1q n u `np2πq " ż π ´π du e i uA
A integer J " ż π
´π du e i uA P.V.

8 ÿ n"´8 1 u `np2πq " ż π ´π du e i uA 1 2 tanpu{2q
We now repeat this simple calculation for the pn ´1q-dimensional integral appearing in [START_REF] Derksen | On the Littelwood-Richardson polynomials[END_REF], evaluated either for unshifted α, β, γ or for shifted α 1 , β 1 , γ 1 , associated as above with a compatible triple of highest weights pλ, µ; νq.

First we observe that the determinant det e i px i ´1 n ř x k qα 1 j that appears in the first line of ( 7) is nothing else than the numerator of Weyl's formula [START_REF] Knutson | Honeycombs and sums of Hermitian matrices[END_REF] for the SUpnq character χ λ pT q, evaluated for the unitary and unimodular matrix

T " diag `ei px i ´1 n ř x k q ˘. ( 24
)
Henceforth we take t i " px i ´1 n ř x k q, ř t i " 0. Consider now the product of three such determinants as they appear in the computation of J n pα 1 , β 1 ; γ 1 q, see [START_REF] Derksen | On the Littelwood-Richardson polynomials[END_REF]. Each factor e i ř j u j A j , under 2π-shifts of the variables u j :" t j ´tj`1 , u j Ñ u j `pj p2πq, is not necessarily periodic, because of the second term of A j in ( 10):

e i ř j u j A j Ñ e i ř j u j A j e ´2πi ř j jp j n ř k pα 1 k `β1 k ´γ1 k q .
Indeed, for α 1 " pλ `ρq, etc, we have

n ÿ k"1 pα 1 k `β1 k ´γ1 k q " n´1 ÿ k"1 kpλ k `µk ´νk q `npn ´1q 2 ,
the first term of which vanishes for a compatible triple pλ, µ; νq, see [START_REF] Ikenmeyer | Small Littlewood-Richardson coefficients[END_REF]. Thus we find that under the above shift, e i ř j u j A j Ñ e i ř j u j A j p´1q ř j jpn´1qp j . For n odd, like in SU(3), the numerator is 2πperiodic in each variable u j . For n even, however, we have a sign p´1q jp j . We may thus compactify the integration domain of the u-variables, bringing it from R n´1 back to p´π, πq n´1 by translations u j Ñ u j `p2πqp j , while taking the above sign into account. Thus for a compatible triple pλ, µ; νq and the A j 's standing for the expressions of (10) computed at shifted weights α 1 " pλ `ρq and likewise for β 1 and γ 1 , we have

ż R n´1 ś n´1 j"1 du j e i u j A j r ∆puq " ż p´π,πq n´1 n´1 ź j"1 du j e i u j A j D n where D n " 8 ÿ p 1 ,¨¨¨,p n´1 "´8 p´1q ř j jp j pn´1q ź 1ďiăi 1 ďn 1 u i `ui`1 `¨¨¨`u i 1 ´1 `pp i `¨¨¨`p i 1 ´1qp2πq , (25)
a sum that always converges. Now define n :"

ź 1ďiăi 1 ďn 2 sinp 1 2 pu i `ui`1 `¨¨¨`u i 1 ´1qq " i ´npn´1q{2 ∆pe i t i q (26) R n pT q :" D n n . (27)
R n , as defined by [START_REF] Rassart | A Polynomiality Property for Littlewood-Richardson Coefficients[END_REF], is a function of T with no singularity, since all the poles of the original expression r ∆puq ´1 have been embodied in the denominator ∆pe i t i q. It must be a polynomial in T and T ‹ , invariant under permutations and complex conjugation, hence a real symmetric polynomial of the e i t j . (Since det T " 1, T ‹ is itself a polynomial in T .) We conclude that R n pT q may be expanded on real characters χ κ pT q, κ P K, with K a finite n-dependent set of highest weights. Moreover R n pIq " 1, as may be seen by looking at the small t limit of [START_REF] Rassart | A Polynomiality Property for Littlewood-Richardson Coefficients[END_REF]. Thus Proposition 1. The integrals over R n´1 appearing in J n pα 1 , β 1 ; γ 1 q in (9), for α 1 " pλ `ρq, β 1 " pµ `ρq, γ 1 " pν `ρq, pλ, µ; νq a compatible triple, may be "compactified" in the form

ż R n´1 ś n´1 j"1 du j e i u j A j r ∆puq " i npn´1q{2 ż p´π,πq n´1 n´1 ź j"1
du j e i u j A j R n pT q ∆pe i t i q (28)

where the real polynomial R n pT q is defined through [START_REF] Rassart | A Polynomiality Property for Littlewood-Richardson Coefficients[END_REF]. There exists a finite, n-dependent set K of highest weights such that R n pT q may be written as a linear combination R n pT q " ř κPK r κ χ κ pT q of real characters. The coefficients r κ are rational and such that, when evaluated at the identity matrix, R n pIq " 1.

Consider now the similar computation, again for a compatible triple pλ, µ; νq but with the A j 's standing for the expressions of (10) computed at unshifted weights, i.e., with α " pλq and likewise for β and γ. If the triple pα, β; γq is non generic, J n pα, β; γq " 0. If it is generic, and n is odd, pα, β; γq may be thought of as associated with the shift of the compatible triple pλ ´ρ, µ ´ρ ; ν ´ρq. Thus for n odd, this new calculation yields the same result as above. For n even, however, the latter triple is no longer compatible and a separate calculation has to be carried out. It is easy to see that the same line of reasoning leads to a modification of the formula [START_REF] Rassart | A Polynomiality Property for Littlewood-Richardson Coefficients[END_REF] and to a new family of real symmetric polynomials p R n pT q, according to

ż R n´1 ś n´1 j"1 du j e i u j A j r ∆puq " ż p´π,πq n´1 n´1 ź j"1 du j e i u j A j p D n p D n :" 8 ÿ p 1 ,¨¨¨,p n´1 "´8 ź 1ďiăi 1 ďn 1 u i `ui`1 `¨¨¨`u i 1 ´1 `pp i `¨¨¨`p i 1 ´1qp2πq (29) p R n pT q :" p D n n , (30)
with the same n as in [START_REF] Racah | Group theoretical concepts and methods in elementary particle physics[END_REF]. Note that the sum in ( 29) is convergent for n ą 2. The case n " 2 requires a special treatment, see below in sec. 4.2.1.

Proposition 2. The integrals over R n´1 appearing in J n pα, β; γq in (9), for α " pλq, β " pµq, γ " pνq, pλ, µ; νq a compatible triple, may be compactified in the form

ż R n´1 ś n´1 j"1 du j e i u j A j r ∆puq " i npn´1q{2 ż p´π,πq n´1 n´1 ź j"1 du j e i u j A j p R n pT q ∆pe i t i q (31)
where the real polynomial p R n pT q is defined through [START_REF] Vergne | Poisson summation formula and box splines[END_REF]. There exists a finite n-dependent set p K of highest weights such that p R n pT q may be written as a linear combination p R n pT q " ř κP p K rκ χ κ pT q of real characters. The coefficients rκ are rational and such that, when evaluated at the identity matrix, p R n pIq " 1. For n odd, the following objects coincide with those of Proposition 1: p R n " R n , p K " K and r κ " rκ .

A method of calculation and explicit expressions for low values of n of the polynomials R n , p R n and of the sets K, p K will be given in sections 2.4 and 4.2, establishing the rationality of the coefficients r κ , rκ . We shall see that the polynomial R n is equal to 1 for n " 2 and n " 3, but non-trivial when n ě 4. In contrast, already for n " 2, p R 2 pT q " 1 2 χ 1 pT q. These expressions of R n and p R n for low n suggest the following conjecture Conjecture 1. The coefficients r κ and rκ are non negative.

As we shall see below in sec. 2.5 (v), this Conjecture 1 is related to Lemma 1.

Relation between J n and LR coefficients

We may now complete the computation of J n pα 1 , β 1 ; γ 1 q and J n pα, β; γq. We rewrite 1 ∆pe i t q " |∆pe i t q| 2 1 ∆pe i t q∆pe i t q∆pe i t q ˚, the first term |∆pe i t q| 2 is what is needed for writing the normalized Haar measure over the SUpnq Cartan torus T n , see [START_REF] King | Stretched Littlewood-Richardson and Kostka coefficients[END_REF], while the three Vandermonde determinants in the denominator provide the desired denominators of Weyl's character formula.

Putting everything together we find Theorem 1. 1. For a compatible triple pλ, µ; νq, the integral J n of (8-9), evaluated for the shifted weights λ `ρ etc, or for the corresponding α 1 " pλ `ρq, β 1 " pµ `ρq, γ 1 " pν `ρq, may be recast as

J n pα 1 , β 1 ; γ 1 q " ż Tn dT χ λ pT qχ µ pT qχ ν pT q R n pT q (32)
where the integration is carried out on the Cartan torus with its normalized Haar measure. Writing R n pT q " ř κPK r κ χ κ pT q as in Prop. 1, this may be rewritten as

J n pα 1 , β 1 ; γ 1 q " ÿ κPK ν 1 r κ N ν 1 λµ N ν 1 κν " ÿ ν 1 c pνq ν 1 N ν 1 λµ ( 33
)
where the sum runs over the finite set of irreps ν 1 obtained in the decomposition of ' κPK pν b κq, with rational coefficients c pνq ν 1 "

ř κPK N ν 1 κν r κ . 2.
For a compatible triple pλ, µ; νq of weights not on the boundary of the Weyl chamber, the integral J n of (8-9), evaluated for the unshifted weights λ, µ, ν, or for the corresponding α " pλq, β " pµq, γ " pνq, may be recast as J n pα, β; γq " ż Tn dT χ λ´ρ pT qχ µ´ρ pT qχ ν´ρ pT q p R n pT q (34)

where the integration is carried out on the Cartan torus with its normalized Haar measure. Writing p R n pT q " ř κP p K rκ χ κ pT q as in Prop. 2, this may be rewritten as

J n pα, β; γq " ÿ κP p K ν 1 rκ N ν 1 λ´ρ µ´ρ N ν 1 κ ν´ρ (35)
"

ÿ ν 1 ĉpνq ν 1 N ν 1 λ´ρ µ´ρ ( 36
)
where the sum runs over the finite set of irreps ν 1 obtained in the decomposition of

' κP p K `pν ´ρqbκ ˘,
with rational coefficients ĉpνq ν 1 "

ř κP p K N ν 1 κ ν´ρ rκ . Proof.
(32) and (34) result from the previous discussion. The product R n pT qχ ν pT q may then be decomposed on characters,

R n pT qχ ν pT q " ÿ κPK r κ χ κ pT qχ ν pT q " ÿ κPK ν 1 N ν 1 κν r κ χ ν 1 pT q " ÿ ν 1 c pνq ν 1 χ ν 1 pT q , with c pνq ν 1 " ř κPK N ν 1 κν r κ , which yields (33). Similarly, p R n χ ν´ρ " ř ν 1 ĉpνq ν 1 χ ν 1 with ĉpνq ν 1 " ř κP p K N ν 1
κ ν´ρ rκ , which gives (36). Recall that if either of λ, µ or ν lies on the boundary of the Weyl chamber, α, β or γ has at least two equal components and J n pα, β; γq " 0. Thus, in words, J n pα 1 , β 1 ; γ 1 q and J n pα, β; γq may be expressed as linear combinations of LR coefficients over "neighboring" weights ν 1 of ν. If Conjecture 1 is right, the coefficients c pνq ν 1 , ĉpνq ν 1 are also non negative. Remark. Note that even though the function J n pα, β; γq is defined for any triple pα, β; γq, compatible or not, integral or not, equations (33),(36) hold only for triples pα 1 , β 1 ; γ 1 q or pα, β; γq associated with compatible triples pλ, µ; νq. Recall also from the previous discussion that for n even, the triple pα 1 , β 1 ; γ 1 q is not integral and compatible if the triple pα, β; γq (or pλ, µ; νq) is. Comment. It would be interesting to invert relations (33,36) and to express the LR coefficients N ν λµ as linear combinations of the functions J n and their derivatives. In view of the considerations of [START_REF] Vergne | Poisson summation formula and box splines[END_REF], this doesn't seem inconceivable6 .

Expression of the R and p R polynomials

Here is the essence of the method used to compute R n and p R n , as defined through ( 27), [START_REF] Vergne | Poisson summation formula and box splines[END_REF]. We first introduce two families of functions, defined recursively

f pu, mq " ´1 m ´1 B Bv f pv, m ´1q| v"
f pu, 1q " 2u 8 ÿ m"1 1 u 2 ´p2πqm 2 `1 u " 1 2 tanpu{2q
and gpu, 1q " 2u

8 ÿ m"1 p´1q m u 2 ´p2πq 2 m 2 `1 u " 1 2 sinpu{2q
.

R n and p R n , defined in [START_REF] Rassart | A Polynomiality Property for Littlewood-Richardson Coefficients[END_REF][START_REF] Vergne | Poisson summation formula and box splines[END_REF], are obtained explicitly by an iterative procedure. We start from 1{ r ∆puq "

ź 1ďiďjďn´1 1 pu i `ui`1 `¨¨¨u j q
First we pick a variable in p r ∆puqq ´1, say u 1 , shift it by p 1 p2πq, perform a partial fraction expansion of the rational function ś 2ďjďn 1 u 1 `¨¨¨`u j´1 `p1 p2πq with respect to the variable u 1 and make use of the previous identities in the summation over p 1 . This produces a sum of trigonometric functions of u 1 , ¨¨¨, u n´1 which are p2πq periodic or anti-periodic in each of these variables, times rational functions of u 2 , ¨¨¨, u n´1 . Then iterate with the variable u 2 , say, shifting it by p 2 p2πq etc. (Of course the order of the variables is immaterial.) As explained in sec. 2.2.2, the final result has the general form

R n presp. p R n q ś 1ďiăi 1 ďn 2 sinp 1 2 pu i `ui`1 `¨¨¨`u i 1 ´1qq where R n , resp. p
R, is a (complicated) trigonometric function of the u variables, or alternatively a symmetric trigonometric function of the t variables. The latter is then recast as a sum of real characters of the matrix T . This procedure will be illustrated in sec. 4.2 on the first cases, for 2 ď n ď 6.

Remark. The reader may have noticed the parallel between this way of computing p R n and the computation of J n in [START_REF] Zuber | Horn's problem and Harish-Chandra's integrals. Probability distribution functions[END_REF]: both rely on an iterative partial fraction expansion, the connection between the two being the Poisson formula. As a consequence of this simple correspondence, J n pα, β; γq evaluated for a compatible triple and p R n have rational coefficients with the same least common denominator δ n , see below Prop. 3. Proposition 3. For any integral compatible triple pα, β; γq, J n pα, β; γq is an integral multiple of some rational number δ ´1 n . Proof. Call δ n the least common denominator of the coefficients ĉν 1 in (36). Then we see that J n pα, β; γq is an integral multiple of 1{δ n .

Unfortunately we have no general expression of δ n and rely on explicit calculations for low values of n: n 2 3 4 5 6 ¨¨δ n 1 1 6 360 9! (iv) Asymptotic behavior. The asymptotic regime is read off (32-36): heuristically, we expect that asymptotically, for rescaled weights, the t-integral in the computation of J n will be dominated by t « 0, hence T « I, for which R n " p R n " 1, whence the asymptotic equality, for λ, µ, ν large

J n pα 1 , β 1 ; γ 1 q « J n pα, β; γq « N ν λµ . (43)
More precisely, it is known [START_REF] Rassart | A Polynomiality Property for Littlewood-Richardson Coefficients[END_REF] that, as a function of ν 1 , N ν 1 λµ can be extended to a continuous piecewise polynomial function, thus for large ν, one approximates the rhs of (33) by

N ν λµ ř ν 1 c ν 1 « N ν
λµ since the coefficients sum up to 1, again as a consequence of R n pIq " 1:

ÿ ν 1 c ν 1 " ÿ κPK r κ ÿ ν 1 N ν 1 κν large ν « ÿ κPK r κ dim V κ " 1
as observed above in (37). We shall see below in sec. 3 that (32,33) enable us to go (a bit) beyond this leading asymptotic behavior.

(v) Compare Conjecture 1 and Lemma 1.

We just observe here that Conjecture 1 is consistent with Lemma 1. Indeed, if we apply (33) to an admissible (hence compatible) triple pλ, µ; νq, with the assumption that the sum over ν 1 includes ν with a non vanishing coefficient c ν , and using the non negativity of the other c ν 1 (as stated in Conj. 1), one obtains J n pα 1 , β 1 ; γ 1 q ě N ν λµ ą 0, in agreement with Lemma 1.

On polytopes and polynomials

The polytopes r H αβ and H ν λµ considered in this section have been introduced in sec. 1.2.

Ehrhart polynomials

Given some rational polytope P, call sP the s-fold dilation of P, i.e., the polytope obtained by scaling by a factor s the vertex coordinates (corners) of P in a basis of the underlying lattice.

The number of lattice points contained in the polytope sP is given by a quasi-polynomial called the Ehrhart quasi-polynomial of P, see for example [START_REF] Stanley | Enumerative Combinatorics[END_REF]. It is polynomial for integral polytopes but one can also find examples of rational non-integral polytopes, for which it is nevertheless a genuine polynomial. We remind the reader that the first two coefficients (of highest degree) of the Ehrhart polynomial of a polytope P of dimension d are given, up to simple normalizing constant factors, by the d-volume of P and by the pd ´1q-volume of the union of its facets; the coefficients of smaller degree are usually not simply related to the volumes of the faces of higher co-dimension. We finally mention the Ehrhart-Macdonald reciprocity theorem: the number of interior points of P, of dimension d, is given, up to the sign p´1q d , by the evaluation of the Ehrhart polynomial at the negative value s " ´1 of the scaling parameter.

Littlewood-Richardson polynomials

It is well known [START_REF] Heckman | Projections of Orbits and Asymptotic Behavior of Multiplicities for Compact Lie Groups[END_REF][START_REF] Guillemin | Symplectic Fibrations and Multiplicity Diagrams[END_REF] that multiplicities like the LR coefficients admit a semi-classical description for "large" representations. In the present context, there is an asymptotic equality of the LR multiplicity N ν λµ , when the weights λ, µ, ν are rescaled by a common large integer s, with the function J n . Here again we assume that the admissible triple pλ, µ; νq is generic, in the sense of Definition 1. Indeed, from (43), as s Ñ 8 N sν sλ sµ « J n p psλ `ρq, psµ `ρq; psν `ρqq « J n psα, sβ; sγq " s pn´1qpn´2q{2 J n pα, β; γq . (44)

The last equality just expresses the homogeneity of the function J n . These scaled or "stretched" LR coefficients have been proved to be polynomial ("Littlewood-Richardson polynomials") in the stretching parameter s [START_REF] Derksen | On the Littelwood-Richardson polynomials[END_REF][START_REF] Rassart | A Polynomiality Property for Littlewood-Richardson Coefficients[END_REF],

N sν sλ sµ " P ν λµ psq (45)
and it has been conjectured that the polynomial P ν λµ psq (of degree at most pn ´1qpn ´2q{2 by ( 44)), has non negative rational coefficients [START_REF] King | Stretched Littlewood-Richardson and Kostka coefficients[END_REF]. More properties of P ν λµ psq, namely their possible factorization and bounds on their degree have been discussed in [START_REF] King | The hive model and the polynomial nature of stretched Littlewood-Richardson coefficients[END_REF]. For a generic triple, our study leads to an explicit value (eq. ( 44)) for the coefficient of highest degree, namely the kernel function J n pα, β; γq, see eq. ( 9).

From the very definition of the hive polytope H ν λ µ associated with an admissible triple (each integral point of which is a honeycomb contributing to the multiplicity), with Littlewood-Richardson, or stretching, polynomial P ν λµ psq, and from the general definition of the Ehrhart polynomial, it is clear that both polynomials are equal. Notice that P ν λµ psq, defined as the Littlewood-Richardson polynomial of the triple pλ, µ; νq or as the Ehrhart polynomial of the polytope H ν λ µ , is polynomial even if the hive polytope happens not to be an integral polytope; on the other hand the Ehrhart polynomial of the polytope defined as the convex hull of the integral points of H ν λ µ will differ from P ν λµ psq if H ν λ µ is not integral, see two examples in sec. 4.4.2 and 4.4.3. From the volume interpretation of the first Ehrhart coefficient, which was recalled in sec. 3.1, we find: Proposition 4. For SUpnq, the normalized d-volume V of the hive polytope H ν λµ equals d! J n pα, β; γq, with d " pn ´1qpn ´2q{2, for a generic and admissible triple pλ, µ; νq, with α " pλq, β " pµq, γ " pνq, and with J n pα, β; γq given by eq. ( 9).

We use here the definition given by [START_REF] Haase | Lecture Notes on Lattice Polytopes, Fall School on Polyhedral Combinatorics[END_REF][START_REF] Bosma | The Magma algebra system. I. The user language[END_REF]: for a polytope of dimension d, the Euclidean volume v is related to the normalized volume V by v " V{d!. More generally the total normalized p-volume V p of the p-dimensional faces of a polytope is related to its total Euclidean p-volume v p by v p " V p {p!. This is consistent with the result [START_REF] Knutson | The honeycomb model of GL n pCq tensor products I: proof of the saturation conjecture[END_REF] that the LR coefficient is equal to the number of integral points in the hive polytope. In words, (44) says that the number of integral points of that polytope is asymptotically well approximated by its euclidean volume J n . The Blichfeldt inequality [START_REF] Blichfeldt | Notes on geometry of numbers, in the October meeting of the San Francisco section of the AMS[END_REF] valid for an integral polytope Q of dimension d, states that its number of integral points is smaller than V `d, where V is its normalized volume. This property, which a fortiori holds for a rational polytope H with integral part Q, together with Proposition 4, implies the following inequality for a generic hive polytope H ν λ µ of SUpnq:

d! J n pα, β; γq ě N ν λµ ´d (46)
with d " pn ´1qpn ´2q{2 and α " pλq, β " pµq, γ " pνq.

Polytopes versus symplectic quotients

Here is another argument relating the volume of the hive polytope with ppγ|α, βq, hence also with J n pα, β; γq, for α " pλq, β " pµq, γ " pνq, λ, µ, ν being dominant integral weights. It goes in two steps, as follows.

Step 1. N ν λµ is the number of integral points of the hive polytope. For large s, the coefficient N sν sλ sµ is approximated by s d times the volume of the same polytope.

Step 2. For large s, N sν sλ sµ is approximated 7 by the volume of a symplectic quotient of the product of three coadjoint orbits labelled by λ, µ, ν, where ν is the conjugate of ν. The same volume is given, up to known constants, by ppγ|α, βq, hence by J n pα, β; γq, see [START_REF] Knutson | Honeycombs and sums of Hermitian matrices[END_REF], Th4. Hence the result.

As already commented in [START_REF] Knutson | Honeycombs and sums of Hermitian matrices[END_REF], the equality between the two volumes is quite indirect and it would be nice to construct a measure preserving map between the hive polytope and the above symplectic quotient, or a variant thereof. To our knowledge, this is still an open problem.

The details of the first part of step 2 are worked out in [START_REF] Suzuki | Asymptotic dimension of invariant subspace in tensor product representation of compact Lie group[END_REF]. We should mention that this last reference also adresses the problem of calculating the function ppγ|α, βq, at least when the arguments are determined by dominant integral weights, and the authors present quite general formulae that are similar to ours. However, they do not use the explicit writing of the orbital measures using formula (3), which was a crucial ingredient of our approach and allowed us to obtain rather simple expressions for J n pα, β; γq.

Subleading term

From the asymptotic behavior (44), we have N sν sλ sµ " P ν λµ psq " s pn´1qpn´2q{2 J n p pλq, pµq; pνqqp1 `Ops ´1qq provided the leading coefficient J n p pλq, pµq; pνqq does not vanish. According to Lemma 1 the stretching polynomial P ν λµ psq is of degree pn ´1qpn ´2q{2 for ν inside the tensor polytope and for λ, µ R BC, but is of lower degree on the boundary of that polytope, or for λ or µ on BC.

Write (33) for stretched weights J n p psλ `ρq, psµ `ρq; psν `ρqq "

ÿ κPK r κ ÿ ν 1 N ν 1 sν κ N ν 1 sλ sµ .
For s large enough, all the weights ν 1 " sν `k, where k runs over the multiset tκu of weights (i.e., counted with their multiplicity) of the irrep with highest weight κ, are dominant and thus contribute to the multiplicity N ν 1 sν κ [START_REF] Racah | Group theoretical concepts and methods in elementary particle physics[END_REF]. Thus J n p psλ `ρq, psµ `ρq; psν `ρqq "

ÿ κPK r κ ÿ kPtκu N sν`k sλ sµ . ( 47
)
But as a function of λ, µ, ν, and in the case of SUpnq, the LR coefficient N ν λ µ is itself a piecewise polynomial [START_REF] Rassart | A Polynomiality Property for Littlewood-Richardson Coefficients[END_REF]: more precisely in the latter reference it is shown that, for the case of SUpnq, the quasi-polynomials giving the Littlewood-Richardson coefficients in the cones of the Kostant complex are indeed polynomials of total degree at most pn ´1qpn ´2q{2 in the three sets of variables defined as the components of the highest weights λ, µ, ν. 7 More precisely limsÑ8 1 s d N sν sλ sµ " ş ω d {d!, with d " pn ´1qpn ´2q{2, where ω is the symplectic 2-form on the symplectic and Kähler manifold of complex dimension d defined as pO λ ˆOµ ˆOν q{{SUpnq :" m ´1p0q{SUpnq, with m, the moment map m : pa1, a2, a3q P O λ ˆOµ ˆOν Þ Ñ a1 `a2 `a3 P LiepSUpnqq ˚.

Remark. The well known Kostant-Steinberg method for the evaluation of the LR coefficients (a method where one performs a Weyl group average over the Kostant function) is not used in our paper, or it is only used as a check. However we should stress that, even in the case of SU [START_REF] Berenstein | Triple multiplicities for sl(r `1) and the spectrum of the external algebra in the adjoint representation[END_REF] where the LR coefficients can be deduced from our kernel function J 3 , see below sec. 4.1.2, the expressions obtained for N ν λµ using the Kostant-Steinberg method differ from ours. If we assume that N ν λ µ may be extended to a function of the same class as J n , namely C n´3 , see above sec. 1.3.2, a Taylor expansion to second order of the rhs of (47) is possible for n ě 4. This leaves out the cases n " 2 and n " 3 which may be treated independently, see below sec. [START_REF] Betke | Lattice points in lattice polytopes[END_REF]

A case by case study for low values of n

We examine in turn the cases n " 2, ¨¨¨, 6.

Expression and properties of the J n function

The expressions of J 2 , J 3 and J 4 were already given in [START_REF] Zuber | Horn's problem and Harish-Chandra's integrals. Probability distribution functions[END_REF]. We repeat them below for the reader's convenience. Those of J 5 and J 6 , which are fairly cumbersome, are available on the web site http://www.lpthe.jussieu.fr/ ~zuber/Z_Unpub.html

The case of SU(2)

In the case of n " 2, the function J 2 reads J 2 pα, β; γq " p1 I pγ 12 q ´1´I pγ 12 qq (49)

where γ 12 :" γ 1 ´γ2 and 1 I is the characteristic function of the segment8 I " p|α 12 ´β12 |, α 12 `β12 q. Then, when evaluated for shifted weights, α 1 " α 12 `1 " λ 1 `1, β 1 " β 12 `1 " µ 1 `1, γ 1 " γ 12 `1 " ν 1 `1 ą 0, it takes the value 1 iff |α 12 ´β12 | ă γ 12 `1 ă α 12 `β12 `2, i.e., iff |α 12 ´β12 | ď γ 12 ď α 12 `β12 which is precisely the well known value of the LR coefficient,

N ν λµ " $ ' & ' % 1 if |α 12 ´β12 | " |λ 1 ´µ1 | ď γ 12 " ν 1 ď α 12 `β12 " λ 1 `µ1 and ν 1 ´|λ 1 ´µ1 | even 0 otherwise . Ê Ê Ê Ê Ê Ê Ê Ê Ê Ê Ê Ê Ê Ê Ê Ê Ê Ê Ê Ê Ê Ê Ê Ê Ê Ê ‡ ‡ ‡ ‡ ‡ ‡ ‡ ‡ ‡ ‡ ‡ ‡ ‡ ‡ ‡ ‡ ‡ ‡ ‡ ‡ ‡ ‡ Ï Ï Ï Ï Ï Ï Ï Ï Ï Ï Ï Ï Ï Ï Ï Ï Ï Ï Ú Ú Ú Ú Ú Ú Ú Ú Ú Ú Ú Ú Ú Ù Ù Ù Ù Ù Ù Ù Ù Á Á Á
Figure 1: The Horn-tensor polygon r H αβ " H λµ for the two SU(3) weights λ " p9, 5q µ " p6, 5q, hence α " p14, 5, 0q, β " p11, 5, 0q. The multiplicity increases from 1 to 6 inside the polygon, giving a matriochka pattern to the successive contours. Proposition 6.

1. For an admissible triple, the function J 3 pα, β; γq of eq. ( 54) takes only values that are integral and non negative; as just discussed, these values vanish by continuity along the edges of the polygon; the vertices of the boundary polygon are integral and give admissible γ's;

2. for α " pλq, β " pµq, γ " pνq, J 3 pα, β; γq " N ν λµ ´1; in particular, if some λ i or µ i vanishes, hence α or β are non generic, N ν λµ " 1, a well-known property of SU(3);

3. the points ν of value J 3 p pλq, pµq; pνqq " m, for 0 ď m ă m max form a "matriochka" pattern, see Fig.

4. Now evaluate J 3 at shifted weights λ 1 " λ `ρ, µ 1 " µ `ρ, ρ the Weyl vector p1, 1q, hence α 1 i " i pλq `3 ´i, β 1 i " i pµq `3 ´i and still α 1 3 " β 1 3 " 0. Then

J 3 pα 1 , β 1 ; γ 1 q " N ν λµ ( 56
)
with ν such that γ 1 i " i pνq `3 ´i, i " 1, 2, 3.

The sum ř

γP r H αβ XZ 2 ∆pγq ∆pαq∆pβq J 3 pα, β; γq equals 1 2 ; therefore replacing the sum by an integral over the domain γ 3 ď γ 2 ď γ 1 , see [START_REF] Harish-Chandra | Differential Operators on a Semisimple Algebra[END_REF], gives the same value (namely 1 2 ).

Proof. Point 1 follows from Proposition 3, with δ 3 " 1. Integrality of the vertices of the polygon is seen by inspection of Horn's inequalities. Point 4 follows from (32) together with the fact that for n " 3, the polynomial R 3 " 1, see below sec. 4.2. Points 2 follows from (56) and the observation made in [START_REF] Coquereaux | Conjugation properties of tensor product multiplicities[END_REF] that, for SU(3),

N ν`ρ λ`ρ µ`ρ " N ν λ µ `1 . (57)
The matriochka pattern of point 3 matches the similar pattern of points of multiplicity m `1 in the tensor product decomposition λ b µ (cf [START_REF] Coquereaux | Conjugation properties of tensor product multiplicities[END_REF], eq [START_REF] Knutson | The honeycomb model of GL n pCq tensor products II: Puzzles determine facets of the Littlewood-Richardson cone[END_REF]]). Point 5 has already been derived in sec. 2.5 and is here a direct consequence of

ř ν N ν λ µ dim V ν " dim V λ dim V µ .
We want to stress a remarkable consequence of the above eq. (54,55,56)

Corollary 1. The LR coefficients N ν λµ of SU(3) may be expressed as a piecewise linear function of the weights λ, µ, ν, sum of the four terms of (54).

To the best of our knowledge, this expression was never given before. Note that the lines of non differentiability of the expression (54) split the plane into at most 9 domains. In each domain, the function J 3 is linear. This is to be contrasted with the known expressions that follow from Kostant-Steinberg formula (see for example [START_REF] Fulton | Representation Theory, A First Course[END_REF], Prop. [25][START_REF] Racah | Group theoretical concepts and methods in elementary particle physics[END_REF][START_REF] Rassart | A Polynomiality Property for Littlewood-Richardson Coefficients[END_REF][START_REF] Stanley | Enumerative Combinatorics[END_REF][START_REF] Suzuki | Asymptotic dimension of invariant subspace in tensor product representation of compact Lie group[END_REF] and which involve a sum over two copies of the SUp3q Weyl group. We should also recall that there exist yet another formula for the multiplicity N ν λµ , stemming from its interpretation [START_REF] Knutson | The honeycomb model of GL n pCq tensor products I: proof of the saturation conjecture[END_REF] as the number of integral solutions to the inequalities on the honeycomb variable,

N ν λ µ " J 3 pα 1 , β 1 ; γ 1 q " minpα 1 1 , ´β1 3 `γ1 2 , α 1 1 `α1 2 `β1 1 ´γ1 1 q (58) ´maxpα 1 2 , γ 1 3 ´β1 3 , γ 1 2 ´β1 2 , α 1 1 `α1 3 `β1 1 ´γ1 1 , α 1 1 `α1 2 `β1 2 ´γ1 1 , α 1 1 ´γ1 1 `γ1 2 q " 1 `minpλ 1 `λ2 , ν 2 `σ, ν 2 ´µ2 `2σq ´maxpλ 2 , σ, ν 2 ´µ2 `σ, ν 2 ´λ2 ´µ2 `2σ, ν 2 ´µ1 ´µ2 `2σ, λ 1 `λ2 ´ν1 q ,
where σ :" 1 3 pλ 1 `2λ 2 `µ1 `2µ 2 ´ν1 ´2ν 2 q. See also [START_REF] Bégin | sup3q k fusion coefficients[END_REF][START_REF] Coquereaux | Conjugation properties of tensor product multiplicities[END_REF] for alternative and more symmetric formulae and [START_REF] Coquereaux | On some properties of SU(3) Fusion Coefficients[END_REF] for an expression in terms of a semi-magic square. Remark. The lines or half-lines of non-differentiability of J 3 , as they appear on expression (54), (see also Figures in [START_REF] Zuber | Horn's problem and Harish-Chandra's integrals. Probability distribution functions[END_REF]), are a subset of the lines along which two arguments of the min or of the max functions of (58) coincide.

The case of SU(4)

The case of SU( 4) is more complicated. Some known features of SU(3) are no longer true. In particular, it is generically not true that multiplicities N ν λµ are equal to 1 on the boundary of the polytope; there is no matriochka pattern, with multiplicities growing as one goes deeper inside the tensor polytope; and relation ( 57) is wrong and meaningless, since pλ `ρ, µ `ρ; ν `ρq cannot be compatible if pλ, µ; νq is.

We first recall the expression of J 4 pα, β; γq given in [START_REF] Zuber | Horn's problem and Harish-Chandra's integrals. Probability distribution functions[END_REF]. With A j standing for A j pP, P 1 , P 2 q in the notations of [START_REF] Fulton | Representation Theory, A First Course[END_REF],

J 4 pα, β; γq " 1 2 3 4! ÿ P,P 1 ,P 2 PS 4 ε P ε P 1 ε P 2 pA 1 q ˆ1 3! pA 2 ´A1 qp|A 3 ´A1 | 3 ´|A 3 ´A2 `A1 | 3 ´|A 3 ´A2 | 3 `|A 3 | 3 q ´1 3 pA 2 qp|A 3 | 3 ´|A 3 ´A2 | 3 q ´1 2 p|A 2 ´A1 | ´|A 2 |qp|A 3 ´A2 |pA 3 ´A2 q `|A 3 |A 3 q ˙. (59)
One can actually restrict the previous triple sum over the Weyl group to a double sum only while multiplying the obtained result by 4!, and this is quite useful for practical calculations. Then, we have, for an admissible triple pλ, µ; νq of h.w. of SUp4q (with λ, µ R BC, i.e., λ i , µ i ‰ 0), and α " pλq, β " pµq, γ " pνq, Proposition 7.

1. N ν λµ ě 4 inside the tensor polytope.

2. J 4 pα, β; γq vanishes when γ belongs to the faces of the polytope r H αβ ; conversely J 4 does not vanish inside the polytope.

3. At these interior points, 6J 4 pα, β; γq, which is the normalized 3-volume V of the hive polytope H ν λµ , is an integer.

4. That integer satisfies V " 6J 4 pα, β; γq ě N ν λµ ´3.

The sum ř

γP r H αβ XZ 3 ∆pγq ∆pαq∆pβq J 4 pα, β; γq equals 1 12 , which matches the normalization [START_REF] Harish-Chandra | Differential Operators on a Semisimple Algebra[END_REF].

Proof. Point 1 results from a general inequality in integral d-polytopes that asserts that their number of integral points is larger or equal to d `1, see [START_REF] Beck | Coefficients and Roots of Ehrhart Polynomials[END_REF], Theorem 3.5. Here for points ν inside the tensor polytope, the polytope H ν λ µ is integral and 3-dimensional, hence d " 3. The first part of point 2 has been already amply discussed, while the second one follows from Lemma 1. Points 3 and 5 have been established in sec. 2.5. Point 4 follows from Blichfeldt's inequality (46).

The consequences of Theorem 1 on the values of J 4 at shifted weights will be discussed in the next subsection. 1. Based on the study of numerous examples, it seems that for weights ν interior to the tensor polytope, we have the lower bound N ν λµ ě 8. Note that the afore mentioned inequality of Theorem 3.5 of [START_REF] Beck | Coefficients and Roots of Ehrhart Polynomials[END_REF] (which would give the weaker N ν λµ ě 7) is no longer applicable, since the hive polytope is not generally integral for n " 5, see a counter-example in sec. [START_REF] Betke | Lattice points in lattice polytopes[END_REF] For n " 2 and n " 3, the polynomial R n is equal to 1. Indeed:

8 ÿ p 1 "´8 p´1q p 1 u 1 `2πp 1 " 1 2 sinpu 1 {2q " i ∆pe i t j q (60) 8 ÿ p 1 ,p 2 "´8 1 pu 1 `2πp 1 qpu 2 `2πp 2 qpu 1 `u2 `2πpp 1 `p2 qq " 1 2 3 sinpu 1 {2q sinpu 2 {2q sinppu 1 `u2 q{2q " i 3 ∆pe i t j q . ( 61
)
On the other hand,

P.V. 8 ÿ p 1 "´8 1 u 1 `2πp 1 " 1 u 1 `8 ÿ p 1 "1 2u u 2 1 ´p2πp 1 q 2 " cospu 1 {2q 2 sinpu 1 {2q " 1 2 i tr T ∆pe i t j q , hence p R 2 pT q " 1 4.2.4 Case n " 6
We have found, after long and tedious calculations 2 7 9! R 6 " 31356 `cospx 1 `x2 `x3 ´x4 ´x5 ´x6 q `perm. : 10 terms in total cospx 1 `x2 `2x 3 ´x4 ´x5 ´2x 6 q `perm. : 90 terms in total

1923

`cospx 1 `x2 `x3 ´x4 ´2x 5 q `perm. : 120 terms in total 284238

`cospx 1 `x2 ´x3 ´x4 q `perm. : 45 terms in total 126 `cosp2x 1 `x2 ´2x 3 ´x4 q `perm. : 180 terms in total 18906

`cospx 1 `x2 ´2x 3 q `perm. : 60 terms in total 1362 `cosp2x 1 ´2x 2 q `perm. : 15 terms in total 1801128

`cospx 1 ´x2 q `perm. : 15 terms in total 4919130

. Alternatively 2 8 9! R 6 pT q " 1699488 `715852χ p1,0,0,0,0q pT qpc.c.q `860χ p2,0,0,0,0q pT qpc.c.q `12032 ´χp0,1,0,0,0q pT qχ p2,0,0,0,0q pT q `c.c. ¯`202683χ p0,1,0,0,0q pT qpc.c.q `124χ p1,1,0,0,0q pc.c.q ´5207 ´χp0,0,1,0,0q pT qχ p1,1,0,0,0q pT q `c.c. ¯`10414χ p0,0,1,0,0q pT qpc.c.q `χp1,0,1,0,0q pT qpc.c.q `6876 ´χp1,0,1,0,0q pT qχ p0,1,0,0,0q pT q `c.c.

"

2629422χ p0,0,0,0,0q pT q `1670 ´χp0,0,1,1,1q pT q `c.c. ¯`24167χ p0,0,2,0,0q pT q `13826 ´χp0,1,0,0,2q pT q `c.c. ¯`216561χ p0,1,0,1,0q pT q `957461χ p1,0,0,0,1q pT q `χp1,0,2,0,1q pT q `125χ p1,1,0,1,1q pT q `985χ p2,0,0,0,2q pT q .

where the last expression is a decomposition as a sum over real representations, with a total dimension 2 8 9!, as it should. We also found :

9! p R 6 pT q " 5422χ p0,0,1,0,0q pT q `χp0,1,1,1,0q pT q `13pχ p0,2,0,0,1q pT q `χp1,0,0,2,0q pT qq 186χ p1,0,1,0,1q pT q `982pχ p0,0,0,1,1q pT q `χp1,1,0,0,0q pT qq .

When evaluated at T " 1 we check that the dimension count is correct: p5422, 1, 13, 186, 982q.p20, 1960, 560 ˆ2, 540, 70 ˆ2q " 9!. We leave it to the reader to write the relations involving N ν λµ that follow from (33) and (36), see an example below in sec. 4.4.3.

Stretching polynomials

The case n " 2

This is a trivial case. Since for any admissible triple, N ν λµ " 1, we have, according to a general result [START_REF] King | Stretched Littlewood-Richardson and Kostka coefficients[END_REF], P ν λµ psq " 1.

4.3.2

The case n " 3

For n " 3, we have, from point 2. in sec. 4.1.2

N ν λµ ´1 " J 3 p pλq, pµq; pνqq and the latter is an homogeneous linear function of s, hence P ν λµ psq " N sν sλ sµ " J 3 p pλq, pµq; pνqq `1 " pN ν λµ ´1qs `1 .

This expression is also valid for weights λ and/or µ on the boundary of the Weyl chamber C, in which case, as is well known ("Pieri's rule"), all LR multiplicities equal 1, and then again by the same general result [START_REF] King | Stretched Littlewood-Richardson and Kostka coefficients[END_REF], P ν λµ psq " 1, while as noticed above, J 3 " 0. Likewise as noticed in sec. 2.2.2, if ν lies on the boundary of tensor polytope, (the outer matriochka), N ν λ µ " 1 and thus again, P ν λµ psq " 1. Remark. The property that P ν λµ psq " 1 `spN ν λµ ´1q had been proved in [START_REF] King | Stretched Littlewood-Richardson and Kostka coefficients[END_REF], then recovered in [START_REF] Rassart | A Polynomiality Property for Littlewood-Richardson Coefficients[END_REF] using vector partition functions.

The case n " 4

For n " 4, given weights λ, µ R BC, and weights ν interior to the polytope, J 4 p pλq pµq; pνqqq ‰ 0 (assuming that Lemma 1 holds true) and the stretching polynomial P ν λµ psq is of degree exactly 3. Now let us Taylor expand J 4 p psλ `ρq, psµ `ρq; psν `ρqq " s 3 J 4 p pλq, pµq; pνqq `1 2 s 2 a `Opsq , where the coefficient a, stemming here from the first order derivatives of J 4 , will receive shortly a geometric interpretation. The stretching polynomial P ν λµ psq must satisfy the three conditions 1. P ν λµ p1q " N ν λµ , by definition;

2. P ν λµ p0q " 1;

3. P ν λµ psq " J 4 p pλq, pµq; pνqqs 3 `1 2 s 2 a `Opsq, as discussed in (48).

Recall now the discussion of sec. 3.1 and 3.2 : J 4 p pλq, pµq; pνqq is 1 6 times the normalized volume V of the hive polytope, and a is half the total normalized area A. There is a unique polynomial satisfying these conditions, namely

P ν λµ psq " J 4 p pλq, pµq; pνqqs 3 `1 4 As 2 `´N ν λµ ´J4 p pλq, pµq; pνqq ´1 4 A ´1¯s `1 " 1 6 V s 3 `1 4 As 2 `pN ν λµ ´1 6 V ´1 4 A ´1qs `1 . (72)
Then the alleged non-negativity of the s coefficient [START_REF] King | Stretched Littlewood-Richardson and Kostka coefficients[END_REF] amounts to The SUp4q hive polytope H ν λµ associated with the branching rule: pp21, 13, 5q, p7, 10, 12q; p20, 11, 9qq. Each integral point (367 of them) stands for a pictograph describing an allowed coupling of this triple, for example the one given in fig. 3. be thought of as a particular point of the tensor polytope and stands itself for a hive polytope of dimension 3 (d " pn ´1qpn ´2q{2 " 3 for SUp4q). It is displayed in fig. 2, right. It has 367 integral points: 160 are interior points, in blue in the figure, and 207 are boundary points. Among the latter, 17 are vertices, in red in the figure, the other boundary points are in brown. The polytope is integral since its vertices are integral -it is always so for SU(4) (see [START_REF] Buch | The saturation conjecture (after A. Knutson and T. Tao)[END_REF], example 2). Every single one of the 367 points of the polytope displayed in fig. 2, right, stands for a pictograph contributing by 1 to the multiplicity of the chosen tensor product branching rule. For illustration, we display one of them on fig. 3; actually we give several versions of this pictograph: first, the isometric honeycomb version and its dual, the O-blade version, and then, the KT-honeycomb version and its corresponding hive. Notice that for the first two kinds of pictographs the external vertices are labelled by Dynkin components of the highest weights, whereas for the last two, they are labelled by Young partitions.

N ν λµ ? ě 1 6 V `1 4 A `1 , (73)
The hive polytope has 12 facets (eight quadrilaterals, three pentagons and one heptagon), 27 edges, and 17 vertices (and Euler's identity is satisfied: 12 ´27 `17 " 2). Its normalized volume and area are V " 1484 and A " 410. The number of pictographs with prescribed edges gives the following sequence of multiplicities N sν sλ sµ " t367, 2422, 7650, 17535, 33561, 57212, 89972, 133325, 188755, 257746, . . .u, for s " 1, 2 . . . Only the first three terms of this sequence are used to determine the LR polynomial if we impose that its constant term be equal to 1: P ν λ,µ psq " p5936s 3 `2460s 2 `388s`24q{4! From our discussion in sec. 3.1, P ν λ,µ psq should be equal to the Ehrhart polynomial Epsq of the hive polytope; using the computer algebra package Magma [START_REF] Bosma | The Magma algebra system. I. The user language[END_REF] we checked that it is indeed so.

The direct calculation of J 4 using (59) gives J 4 p pλq, pµq; pνqq " 742{3, and more generally J 4 p psλq, psµq; psνqq " 742 s 3 {3. Using the same eq. ( 59), we can also calculate J 4 for ρ-shifted arguments: J 4 p psλ `ρq, psµ `ρq; psν `ρqq " ment with our general discussion of sec. 3.4, the first two terms of P ν λµ psq and of J 4 p psλ ρq, psµ `ρq; psν `ρqq are identical, the leading term being also equal to J 4 p psλq, psµq; psνqq. One checks that the leading coefficient of Epsq, hence of P ν λ,µ psq, is equal to 1 3! of the normalized volume of the polytope and that the second coefficient is equal to 1 2 1 2! of the normalized 2-volume of its boundary. In accordance with Ehrhart-Macdonald reciprocity theorem, one also checks that ´P ν λµ p´1q " 160, the number of interior points in the polytope. Finally, on this example, one can test eq (66) which relates I 4 p pλ `ρq, pµ `ρq; pν `ρqq " 11592 to a sum of the Littlewood-Richardson coefficient N ν λµ and its twelve "neighbors" ν `α appearing in the tensor product ν b p1, 0, 1q. Likewise eq (67) relates 24J 4 p pλq, pµq; pνqq " 1484 to a sum over six weights ν 2 " p18, 10, 9q, p18, 11, 7q, p19, 9, 8q, p19, 11, 8q, p20, 9, 9q, p20, 10, 7q of the product N ν 2 λ´ρ µ´ρ N ν 2 ν´ρ p0,1,0q which takes the respective values 254, 235, 254, 243, 259, 239, the sum being indeed 1484.

An example in SU(5)

Consider the following tensor branching rule of SUp5q: pλ, µ; νq with λ " p1, 3, 2, 3q, µ " p2, 1, 4, 2q, ν " p3, 1, 4, 3q. The hive polytope H ν λµ has dimension d " 6. We shall see that it is not an integral polytope. We denote Q the convex hull of its integral points. H ν λµ has 66 vertices and 99 points, all of them being boundary points. Q has 64 vertices and 99 points (the latter being the same as for H ν λµ , by definition). Therefore we see that 2 vertices of H ν λµ are not (integral) points of H ν λµ . The normalized volume of H ν λµ is 2544 (it is 2538 for Q). The normalized volume of the boundary of H ν λµ is 3630 (it is 3618 for Q). The polynomial P ν λµ psq, i.e., the Ehrhart polynomial of H ν λµ , is 53s 6 {15 `121s 5 {8 `667s 4 {24 `679s 3 {24 `687s 2 {40 `73s{12 `1. In the case of H ν λµ , we check the first two coefficients related to the 6-volume of the polytope and to the 5-volume of the facets: 2544{6! " 53{15 and 1{2 ˆ3630{5! " 121{8. The Ehrhart polynomial of Q is 141s 6 {40 `603s 5 {40 `665s 4 {24 `679s 3 {24 `259s 2 {15 `92s{15 `1. In the case of Q, the same volume checks read: 2538{6! " 141{40 and 1{2 ˆ3618{5! " 603{40.

An independent calculation using the function J 5 gives J 5 p pλq, pµq; pνqq " 53{15, the leading coefficient of the stretching polynomial.

In the present example, where H ν λµ and Q differ, it is instructive to consider what happens under scaling. The two vertices of H ν λµ that are not integral points are actually half-integral points, so that they become integral by doubling. The polytope 2H ν λµ has again 66 vertices (by construction), it is integral, it has 1463 points, 18 being interior points and 1445 being boundary points. It could also be constructed as the hive polytope associated with the doubled branching rule p2λ, 2µ; 2νq, and its own Littlewood-Richardson (LR) polynomial, equal to its Ehrhart polynomial, can be obtained from the LR polynomial of H ν λµ by substituting s to 2s. The polytope 2Q has again 64 vertices (of course), it is integral, it has 1460 points, 18 being interior points ans 1442 being boundary points. Since Q Ă H ν λµ we have 2Q Ă 2H ν λµ , but now both polytopes are integral (and they are different). Q and H ν λµ have the same integral points, so, in a sense, they describe the same multiplicity for the chosen triple pλ, µ; νq, however, under stretching (here doubling) of the branching rule, we have to consider 2H ν λµ , not 2Q, otherwise we would miss three honeycombs (" 1463 ´1460) and find an erroneous multiplicity. These three honeycombs correspond to the two (integral) vertices of 2H ν λµ coming from the two (non integral) vertices of H ν λµ that became integral under doubling, plus one extra (integral) point, which is a convex combination of vertices. For illustration purposes we give below the three pictographs (in the O-blade version) that correspond to these three points. We consider the following tensor branching rule of SUp6q: pλ, µ; νq with λ " p1, 3, 1, 2, 1q, µ " p2, 1, 3, 2, 1q, ν " p4, 1, 6, 2, 1q. The multiplicity is 38.

For SU [START_REF] Buch | The saturation conjecture (after A. Knutson and T. Tao)[END_REF], the number of fundamental pictographs is 5 ˆ3 `2 ˆ10 but there are 10 syzygies (one for each inner hexagon in the honeycomb picture) so that a basis has 25 elements, the set of 38 (integral) honeycombs is then described as a 25 ˆ38 matrix. The convex hull of these 38 points is then calculated, one finds that it is a 10 dimensional polytope Q (in R 25 ). The obtained polytope -which has no interior point and 38 integral points, 36 of them being vertices-happens not to coincide with the hive polytope H (we are in a situation analogous to the one examined in the previous SU(5) example). A quick study of Q reveals that this polytope, and so H itself, has dimension 10, and that the chosen triple is therefore generic. The fact that H differs from Q can be seen in (at least) three different ways: 1) The Ehrhart polynomial of Q fails to recover the multiplicity of ps λ, s µ; s νq, already for s " 2 where the multiplicity is 511. 2) The leading coefficient (30{9!) of this polynomial, hence the normalized volume of Q, differs from J 6 p pλq, pµq; pνqq " 32{9! determined directly or from Theorem 1 (part 2), we shall come back to this below. 3) A direct determination of the polytope H obtained as an intersection of 45 half-spaces -interpreted for instance as the number of (positive) edges in the oblade picture-will show that H is not an integral polytope (its vertices, aka corners, are rational but not all integral) and its integral part is indeed Q. We leave this as an exercise to the reader. The LR-polynomial associated with the chosen triple, equivalently the Ehrhart polynomial of H, is equal to The coefficient of s 10 , equal to 1{11340 " 32{9! and interpreted as the normalized volume of H, can be obtained from a direct evaluation of the expression of J 6 , but it can also be obtained easily from Theorem 1 (part 2). This double sum (35) involves the seven weights κ together with the seven associated coefficients p r κ that appear in (70) and turns out to involve only the following weights ν 1 : p1, 2, 2, 2, 0q, p1, 2, 3, 0, 1q, p2, 1, 2, 1, 1q, p2, 1, 3, 0, 0q. Most terms are actually zero (because of the vanishing of many Littlewood-Richardson coefficients), and the result is p1`2`2`1`13`13q{9! " 32{9!.

s

4. 1 . 4 A

 14 few facts about SU[START_REF] Blichfeldt | Notes on geometry of numbers, in the October meeting of the San Francisco section of the AMS[END_REF]

2 A

 2 while the counting of interior points, through Ehrhart-Macdonald reciprocity theorem, gives us another lower bound on N ν λµ #pinterior pointsq " ´P ν λµ p´1q " N ν λµ ´p 1 `2q ě 0 .

Figure 2 :

 2 Figure2: Left: The SUp4q tensor polytope H λµ for λ " p21, 13, 5q, µ " p7, 10, 12q , and its 7092 integral points (distinct irreps). Each such point can itself be thought as a hive polytope, for example the one given on the right. Right:The SUp4q hive polytope H ν λµ associated with the branching rule: pp21, 13, 5q, p7, 10, 12q; p20, 11, 9qq. Each integral point (367 of them) stands for a pictograph describing an allowed coupling of this triple, for example the one given in fig.3.

742 3 s 3 7 O

 37 `205 2 s 2 `12s `5 12 . In agree--blade version: edges are non-negative integers, opposite angles (sum of adjacent edges) around the inner points are equal.

  Isometric honeycomb version: opposite angles (sum of adjacent edges) of hexagons are equal.

Figure 3 :

 3 Figure3: One of the 367 pictographs associated with pp21, 13, 5q, p7, 10, 12q; p20, 11, 9qq. For completeness we also give below the corresponding KT-honeycomb and its dual hive.

Figure 4 :

 4 Figure 4: The three SUp5q pictographs (O-blade version) associated with p2λ, 2µ; 2νq, with λ "p1, 3, 2, 3q, µ " p2, 1, 4, 2q, ν " p3, 1, 4, 3q that belong to the hive polytope of this doubled branching rule but that do not belong to the double of the integral part of the hive polytope of pλ, µ; νq.

29 4. 4 . 3

 2943 An example in SU[START_REF] Buch | The saturation conjecture (after A. Knutson and T. Tao)[END_REF]

  Take two orbits of the group Upnq acting on H n , labelled by Hermitian matrices A and B, and consider the corresponding orbital measures m A , m B . The convolution product of the latter is defined as usual: with f , a function on H n , one sets

	)
	1.1.4 Convolution product of orbital measures

  As mentioned above, it is natural to adopt the following definition Definition 1. A triple pα, β; γq is called generic if J n pα, β; γq is non vanishing.

	which equals 1, 1 2 , 1 12 , 1 288 , 1 34560 , ¨¨¨for n " 2, 3, 4, ¨¨¨.
				2;
	-(v) it vanishes for ordered γ outside r H αβ ;
	-(vi) by continuity (for n ě 3) it vanishes for γ at the boundary of r H αβ ;
	-(vii) it also vanishes whenever at least two components of α or of β coincide 4 : this follows from
	the antisymmetry mentionned above;
	-(viii) its normalization follows from that of the probability density p, (normalized of course by
	ş R n d n γ ppγ|α, βq " 1), hence
	ż H αβ r	d n´1 γ	∆pγq ∆pαq∆pβq	J n pα, β; γq "	1 sf pn ´1q	(13)

  Thus for generic points, the two polynomials J n p psλ `ρq, psµ `ρq; psν `ρqq and P ν λµ psq have the same two terms of highest degree d max " pn ´1qpn ´2q{2 and d max ´1. In the degenerate case where the term of degree d max vanishes and the next does not, the leading terms of degree d max ´1 are equal. If the degree is strictly lower than d max ´1, there is no obvious relation between the two polynomials, see examples at the end of sec. 4.3.3.

							.1.1 and
	4.1.2. We thus Taylor expand for large s
	J n p psλ `ρq, psµ `ρq; psν `ρqq "	ÿ	r κ	ÿ	P λ µ ν`k{s	psq
				κPK		kPtκu
			"	ÿ κPK	r κ ¨dim V κ P ν λ µ psq	`1 s	kPtκu ÿ	k∇ ν P ν λ µ psq `¨¨¨'
			" P ν λ µ psq ˆ1	`o´1 s	¯˙(48)
	since	ř	κPK r

κ dim V κ " 1 as noticed above in sec. 2.2, and ř kPtκu k " 0 in any irrep.

  .4.2. 2. J 5 pα, β; γq vanishes outside (and on the boundary) of the polytope, as already discussed. 3. For a compatible triple pα, β; γq and γ inside the polytope r H αβ , 360J 5 pα, β; γq is a positive integer (see sec. 2.5), provided α and β have only distinct components. It is non vanishing according to Lemma 1. Moreover N ν λµ ď 6! J 5 pα, β; γq `6 according to (46). As in section 2.2 the notation χ λ denotes the character of the Lie group SUpnq associated with the irrep of highest weight λ. Also recall that for n odd, p R n " R n .

	4.	ř γP r H αβ XZ 4 J 5 pα, β; γq ∆pγq ∆pαq∆pβq " 1 288 , see (13) again.
	4.2 The polynomials R n and p R n . Application of Theorem 1
	4.2.1 Cases n " 2 and n " 3

In practice we use the normalized Haar measure that makes the volume of Upnq equal to 1.

The context being specified, people often simply write "Fourier transform" or "Fourier orbital transform" rather than "spherical transform" or "orbital transform".

The reader may look at[START_REF] Coquereaux | Conjugation properties of tensor product multiplicities[END_REF] for an explicit descriptions and a few examples of O-blades and isometric honeycombs in the framework of the Lie group SU[START_REF] Berenstein | Triple multiplicities for sl(r `1) and the spectrum of the external algebra in the adjoint representation[END_REF]. See also our SUp4q example in sec.

4.4.1.

If α and β are Young partitions describing the highest weights λ, µ of two Upnq or SUpnq irreps, this occurs when some Dynkin label of λ or µ vanishes, i.e., when λ or µ belongs to a wall of the dominant Weyl chamber C.

We thank Allen Knutson for pointing this out to us.

Our thanks to Michèle Vergne for pointing to that possibility.

This result should be connected with the fact that the support of the convolution product of measures on concentric 2-spheres is an annulus.

χ 1 pT q, while p R

" R 3 " 1.

Acknowledgements

We acknowledge stimulating discussions with Olivier Babelon, Paul Zinn-Justin and especially Allen Knutson and Michèle Vergne.

Consequences of Theorem 1

(i) We start with a useful lemma Lemma 2. With the notations of Theorem 1, we have the relations ÿ κPK r κ dim V κ " 1 (37)

Proof. From the relation R n pT q " ř κPK r κ χ κ pT q evaluated at T " I, with R n pIq " 1, it follows that ř κPK r κ dim V κ " 1. Then

because of the reality of the irreps of h.w. κ, hence

The two relations (38) and (40) are proved in the same way.

(ii) Localization of the normalization integral of J n .

For two given integral (non negative) α and β, consider the sum of J n pα, β; γq∆pγq over the integral γ's inside the connected part r H αβ of the support of J n . If either α or β is non generic, (i.e., has two equal components), all J n pα, β; γq vanish.

Conversely if both α and β are generic, i.e., λ and µ are not on the boundary of the Weyl chamber, we make use of ( 19) and (36)

by Lemma 2. (The ν's on the boundary of the Weyl chamber, for which ν ´ρ is not dominant, do not contribute because of the vanishing of J n pα, β; γq.) Comparing with [START_REF] Harish-Chandra | Differential Operators on a Semisimple Algebra[END_REF], we find that

In others words, the normalization integral of J n over the sector γ n´1 ď ¨¨¨ď γ 1 localizes over the integral points of that sector.

(iii) Quantization of J n .

We conclude that

in agreement with the general formula (33), provided we assume that the indicator function vanishes at the end points of the interval I.

On the other hand, as we shall see below in sec. 4.2.1, p R 2 " 1 2 χ 1 pT q, so that (36) amounts to

which is consistent with (49) if we assume now that the indicator function takes the value 1 2 at the end points of the interval I. This rather peculiar situation is a consequence of the irregular, discontinuous, structure of J 2 .

The case of SU(3)

For n " 3, J 3 takes a simple form within the tensor polytope (here a polygon). In [START_REF] Zuber | Horn's problem and Harish-Chandra's integrals. Probability distribution functions[END_REF], the following was established.

The function

with A 1 and A 2 as in [START_REF] Fulton | Representation Theory, A First Course[END_REF], may be recast in a more compact form:

Proposition 5. Take α 1 ě α 2 ě α 3 , and likewise for β. For γ satisfying (12), Horn's inequalities and γ 1 ě γ 2 ě γ 3 ,

where

pγ 2 ´α3 ´β1 q ´pγ 1 ´α1 ´β2 q if γ 2 ´α3 ´β1 ě 0 and γ 1 ´α1 ´β2 ă 0 pγ 3 ´α2 ´β3 q ´pγ 2 ´α3 ´β1 q if γ 3 ´α2 ´β3 ě 0 and γ 2 ´α3 ´β1 ă 0 pγ 1 ´α1 ´β2 q ´pγ 3 ´α2 ´β3 q if γ 1 ´α1 ´β2 ě 0 and γ 3 ´α2 ´β3 ă 0

.

J 3 pα, β; γq takes non negative values inside the tensor polygon and vanishes by continuity along the edges of the polygon. It also vanishes whenever two components of α or β coincide (non generic orbits).

The non-negativity follows from the interpretation of J 3 as proportional with a positive coefficient to the PDF p.

Consider now an admissible triple pλ, µ; νq of highest weights of SUp3q. The associated triple pα, β; γq is defined as explained above,

3 pλ 1 `2λ 2 `µ1 `2µ 2 ´ν1 ´2ν 2 q, an integer, so that ř 3

i"1 pγ i ´αi ´βi q " 0. Then

In contrast, for n ě 4, one finds non trivial polynomials R n pT q and p R n pT q. For instance for n " 4, with the notations D 4 , p D 4 and 4 introduced in ( 26)

and likewise

Now, in SU(4), we can write χ p1,0,1q pT qχ ν pT q " χ ν pT q `ÿ ν 1 χ ν 1 pT q χ p0,1,0q pT qχ ν´ρ pT q " ÿ ν 2 χ ν 2 pT q with a sum over the h.w. ν 1 , resp. ν 2 , appearing in the decomposition of ν b p1, 0, 1q, resp. of pν ´ρq b p0, 1, 0q. Notice that p1, 0, 1q is the highest weight of the adjoint representation, hence one may write ν 1 " ν `α where α runs over the 12 non zero roots α for ν "deep enough" in the Weyl chamber, i.e., provided all ν `α are dominant weights, and over three times the weight 0 . Thus we may write J 4 p pλ `ρq, pµ `ρq; pν `ρqq "

where ∆N ν λ µ :"

´N ν λ µ q may be regarded as a second derivative term (a discretized Laplacian), while the "first derivative" term vanishes because of ř α " 0.

Example: Take λ " p1, 2, 2q, µ " p2, 2, 1q, ν " p1, 4, 1q, the ν 1 and their multiplicities read pν 1 , N ν 1 ν p1,0,1q q " tp0, 3, 2q, 1q, pp0, 4, 0q, 1q, pp0, 5, 2q, 1q, pp0, 6, 0q, 1q, pp1, 3, 3q, 1q, pp2, 2, 2q, 1q, pp2, 3, 0q, 1q, pp2, 4, 2q, 1q, pp2, 5, 0q, 1q, pp3, 3, 1q, 1q, pp1, 4, 1q, 3qu , J 4 p pλ `ρq, pµ `ρq; pν `ρqq " 97{24 while N ν λµ " 5, ř ν 1 N ν 1 ν p1,0,1q N ν 1 λ µ " 52, the rhs of (66) equals 97{24, and matches the lhs. Note that in that example, only 10 out of the 12 α contribute.

There is a second relation, which follows from (36) with the above expression of p R 4 J 4 p pλq, pµq; pνqq "

For the previous example λ " p1, 2, 2q, µ " p2, 2, 1q, ν " p1, 4, 1q, three weights ν 2 contribute N ν 2 ν´ρ p0,1,0q " 1, namely p0, 2, 0q, p1, 2, 1q, p0, 4, 0q, but only the first two give N ν 2 λ´ρ µ´ρ " 1, the third has N ν 2 λ´ρ µ´ρ " 0, and the rhs equals 1 3 , which is the value of J 4 p pλq, pµq; pνqq.

Case n " 5

For n " 5, likewise R 5 pT q " 1 180 " 45 `12 `cospx 1 ´x2 q `perm. : 10 terms in total cospx 1 `x2 ´x3 ´x4 q `perm. : 15 terms in total ˘ı " 7 72 `1 40 tr T tr T ‹ `1 1440 rptr T q 2 ´tr T 2 src.c.s " 7 72 `1 40 χ p1,0,0,0q pT qχ p1,0,0,0q pT ‹ q `1 360 χ p0,1,0,0q pT qχ p0,1,0,0q pT ‹ q " 1 360

`45 `10χ p1,0,0,1q pT q `χp0,1,1,0q pT q ˘.

Comment: note that at T " I, 45 `10 ˆ24 `75 " 360, R 5 pIq " 1, as it should. Then denoting the h.w. appearing in p1, 0, 0, 1q b ν, resp. p0, 1, 1, 0q b ν, by ν 1 , resp. ν 2 , R 5 pT qχ ν pT q " 1 360 ´45χ ν pT q `10 ÿ ν 1 χ ν 1 pT q `ÿ ν 2 χ ν 2 pT q ānd 360J 5 p pλ `ρq, pµ `ρq; pν `ρqq " 45N ν λ µ `10

Here again, for ν "deep enough" in C, we can make the formula more precise: ν 1 ´ν runs over the 24 weights (=roots) of the adjoint representation p1, 0, 0, 1q, including 4 copies of 0 and 20 non zero roots α; likewise ν 2 ´ν runs over the 75 weights of the p0, 1, 1, 0q representation, including 5 copies of 0, twice the 20 α and the 30 weights β of the form ˘pα ij ˘α kl q with 1 ď i ă j ă k ă l ď 5 or ˘pα ij `α kl q with 1 ď i ă k ă j ă l ď 5. Here we are making use of the notations αi , 1 ď i ď 4 for the simple roots, and αij " αi `¨¨¨`α j´1 with 1 ď i ă j ď 5 for the positive roots . Thus "deep enough" actually means: all ν `α and ν `β P C. Then (68) reads J 5 p pλ `ρq, pµ `ρq; pν `ρqq " N ν λ µ `1 30

(with 20{30 `30{360 " 3{4).

Example. λ " p2, 3, 3, 2q, µ " p3, 2, 3, 2q, ν " p5, 3, 2, 3q, N ν λµ " 211. We find in the lhs of (68) 360J 5 p pλ `ρq, pµ `ρq; pν `ρqq " 63213 while the three terms in the rhs equal respectively 9495, 42010, 11708 with a sum of 63213, qed.

In [START_REF] Betke | Lattice points in lattice polytopes[END_REF][START_REF] Beck | Coefficients and Roots of Ehrhart Polynomials[END_REF] inequalities were obtained between coefficients of the Ehrhart polynomial of an integral polytope. Recall that for n " 4, all hive polytopes are integral [START_REF] Buch | The saturation conjecture (after A. Knutson and T. Tao)[END_REF], and we may apply on (72) these inequalities which read

which is precisely the Blichfeldt inequality mentioned above at point 4 of sec. 4.1.3.

In contrast, for non generic triples pλ, µ; νq, J 4 p pλq, pµq; pνqq " 0, the stretching polynomial is of degree strictly less than 3, and reads in general

If the coefficient a is non vanishing, it has now to be interpreted as the normalized area of the 2-dimensional hive polytope (a polygon). If a " 0, either N ν λµ ě 2 and P ν λµ psq " pN ν λµ ´1qs `1, or N ν λµ " 1 and P ν λµ psq " 1, consistent with the result of sec. 3.4 and the two general results P " 1 if N ν λµ " 1 and P " s `1 if N ν λµ " 2. In the former case (dimension 2 polytope, degree 2 Ehrhart polynomial), Erhrart-Macdonald reciprocity theorem gives us an upper bound on N ν λµ ď a `2, while the alleged non-negativity of the s-coefficient gives a lower bound, N ν λµ ě 1 2 pa `2q. Thus one should have

Also denoting c :" #internal points " P p´1q " a ´N ν λµ `2, b " # boundary points, b `c :" #total of points " N ν λµ , hence a `2 " b `2c which is Pick's formula for the Euclidean area a{2 " b{2 `c ´1.

Examples: Here we denote for short J 1 4 " J 4 p psλ `ρq, psµ `ρq; psν `ρqq. Take λ " p2, 2, 1q, µ " p2, 1, 3q, for ν " p0, 1, 4q, N ν λµ " 3, P ν λµ psq " 1 2 ps `1qps `2q, J 1 4 " 1 12 p6s 2 `15s `7q while for ν " p2, 4, 0q, N ν λµ " 3, P ν λµ psq " 2s `1, J 1 4 " 1 2 p1 `4sq and for ν " p2, 0, 4q, N ν λµ " 4, P ν λµ psq " ps `1q 2 , J 1 4 " 1 4 p4s 2 `7s `2q. Take λ " p3, 0, 3q, µ " p2, 3, 1q, for ν " p3, 4, 0q, N ν λµ " 3, P ν λµ psq " 2s `1, J 1 4 " 1 8 p14s `5q, while for ν " p2, 3, 1q, N ν λµ " 6, P ν λµ psq " ps `1qp2s `1q, J 1 4 " 1 8 p16s 2 `18s `3q. Consider the irreps of highest weight λ " p21, 13, 5q and µ " p7, 10, 12q. Their tensor product contains 7092 distinct irreps ν with multiplicities ranging from 1 to 377. The tensor polytope H λµ is displayed in fig. 2, left. The total multiplicity (sum of multiplicities for the various ν's) is 537186.

Let us now consider a particular term in the decomposition of the tensor product into irreps: the admissible triple pλ, µ; νq, with ν " p20, 11, 9q, whose multiplicity is equal to 367. This term can