This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Differential algebra of polytopes and inversion formulas

V.M. Buchstaber Steklov Mathematical Institute and Moscow State University, Russia buchstab@mi-ras.ru  and  A.P. Veselov Department of Mathematical Sciences, Loughborough University, Loughborough LE11 3TU, UK A.P.Veselov@lboro.ac.uk
Abstract.

We use the differential algebra of polytopes to explain the known remarkable relation of the combinatorics of the associahedra and permutohedra with the universal compositional and multiplicative inversion formulas for the formal power series. This approach allows to single out the associahedra and permutohedra among all graph-associahedra and emphasizes the significance of the differential equations for special sequences of simple polytopes derived earlier by one of the authors. We discuss also the link with the geometry of Deligne-Mumford moduli spaces M¯0,n\bar{M}_{0,n} and permutohedral varieties, as well as the interpretation of the combinatorics of cyclohedra in relation with the classical Faà di Bruno’s formula.

1. Introduction

The following natural question goes back at least to Lagrange. For a given power series α(u)=u+n1anun+1\alpha(u)=u+\sum_{n\geq 1}a_{n}u^{n+1} find the inversion β(v)\beta(v) under the substitution: α(β(v))v,β(α(u))u.\alpha(\beta(v))\equiv v,\,\,\beta(\alpha(u))\equiv u. The coefficients of such series β(v)=v+n1bnvn+1\beta(v)=v+\sum_{n\geq 1}b_{n}v^{n+1} are known to be certain polynomials with integer coefficients:

b1=a1,b2=a2+2a12,b3=a3+5a1a25a13,b_{1}=-a_{1},\,b_{2}=-a_{2}+2a_{1}^{2},\,\,b_{3}=-a_{3}+5a_{1}a_{2}-5a_{1}^{3},
b4=a4+6a1a3+3a2221a12a2+14a14.b_{4}=-a_{4}+6a_{1}a_{3}+3a_{2}^{2}-21a_{1}^{2}a_{2}+14a_{1}^{4}.

Remarkably these coefficients can be interpreted in terms of the combinatorics of the associahedra (or, Stasheff polytopes), see [5, 41] for the details. For example, the formula for b4b_{4} says that 3D3D associahedron has 6 pentagonal and 3 quadrilateral faces, 21 edges and 14 vertices (see Fig. 1, where we show also its realisation as a truncated cube following [14]).

Refer to caption
Refer to caption
Figure 1. 3D associahedron and its realisation as a truncated cube.

It is difficult to trace who was the first to observe this remarkable fact. The first combinatorial description of the coefficients bnb_{n} seems to belong to Raney [53]. Stanley [56] provided 3 proofs in terms of tree combinatorics, known to be closely related to the associahedron. Loday [41] described this explicitly in terms of the associahedra, see also recent book by Aguiar and Ardila [5], who used the theory of Hopf monoids to derive this connection.

The aim of this paper is to provide one more derivation of this link using the differential equations for the generating function of the associahedra derived by one of the authors [11]. This allows us also to prove another known remarkable formula, expressing the multiplicative inversion of a formal power series in terms of combinatorics of permutohedra (see [5]) and to provide an interpretation of the combinatorics of cyclohedra in relation with the classical Faà di Bruno’s formula, which was first discovered by Aguiar and Bastidas [4]. In the case of the stellohedra we derive the relation with permutohedra, which seems to be new.

We discuss also the link with the rich theory of Deligne-Mumford moduli spaces M¯g,n\bar{M}_{g,n}, which are known to be related to associahedra [36], inversion formulas [44] and to the KdV hierarchy [60].

2. Differential algebra of simple polytopes

We start with a brief description of the differential algebra of simple polytopes following [11, 15]. For the general theory of convex polyhedra we refer to Ziegler [61].

Let 𝔓:=n0𝔓n\mathfrak{P}:=\sum_{n\geq 0}\mathfrak{P}^{n} be the free abelian group generated by all convex polytopes, considered up to combinatorial equivalence and naturally graded by the dimension of the polytopes. The product of polytopes turns 𝔓\mathfrak{P} into a graded commutative ring with the unit given by the point, considered as 0-dimensional polytope. Simple polytopes form a graded subring 𝔖𝔓.\mathfrak{S}\subset\mathfrak{P}. Note that the ring 𝔓\mathfrak{P} is different both from the polytope algebra [46] of convex polytopes in n\mathbb{R}^{n} with product given by the Minkowski sum and from the Grothendieck ring of Γ\Gamma-rational polytopes considered in [48].

A polytope is called indecomposable if it cannot be represented as a product of polytopes of positive dimension.

Theorem 2.1 ([13]).

The ring 𝔓\mathfrak{P} is a polynomial ring generated by indecomposable combinatorial polytopes, so any combinatorial polytope can be uniquely represented as a product of the indecomposable polytopes.

The proof can be found in [13] (see also Prop.1.7.2 in [15]). Note that if we introduce the second grading in 𝔓\mathfrak{P} by the number of the facets of polytopes, then the number of generators in any bi-graded component of 𝔓\mathfrak{P} will be finite [15].

Following [11] introduce the derivation dd in 𝔓\mathfrak{P} by defining the boundary dPdP of a given polytope PP simply as the sum of all facets of PP considered as elements in 𝔓\mathfrak{P}. It is easy to check that dd preserves the subring 𝔖𝔓\mathfrak{S}\subset\mathfrak{P} and satisfies the Leibnitz identity in 𝔓\mathfrak{P}

d(P1P2)=(dP1)P2+P1(dP2).d(P_{1}P_{2})=(dP_{1})P_{2}+P_{1}(dP_{2}).

This supplies 𝔓\mathfrak{P} with a structure of the differential ring, with 𝔖\mathfrak{S} being its differential subring. Note that in our case d20d^{2}\neq 0 in contrast with the differential introduced in [29].

We will be interested in the special class of simple polytopes known as graph-associahedra. They were introduced and studied in the work of Carr and Devadoss [20] and independently by Toledano Laredo [59] under the name De Concini–Procesi associahedra (see also Postnikov [51]). This class contains the classical series of permutohedra and associahedra.

Let Γ\Gamma be a connected simple (no loops or multiple edges) graph with the set of nn vertices, which we identify with [n]:={1,2,,n}.[n]:=\{1,2,\dots,n\}. The corresponding graph-associahedron PΓP_{\Gamma} is a particular case of the nestohedron [27] with the building set =(Γ)\mathcal{B}=\mathcal{B}(\Gamma) consisting of the connected induced subgraphs of Γ\Gamma, or equivalently as the subsets S[n]S\subset[n] such that the restriction Γ|S\Gamma|_{S} is connected. Explicitly PΓP_{\Gamma} can be realised as the following (n1)(n-1)-dimensional convex polytope

(1) PΓ={xn:i=1nxi=3n,iSxi3|S|,S,S[n]},P_{\Gamma}=\{x\in\mathbb{R}^{n}:\sum_{i=1}^{n}x_{i}=3^{n},\,\,\sum_{i\in S}x_{i}\geq 3^{|S|},\,S\in\mathcal{B},\,S\neq[n]\},

where |S||S| is the number of elements in SS (see [59]).

The boundary of the graph-associahedra can be given by the following formula (see e.g. [15], Prop. 1.7.16). For any connected induced subgraph SS of Γ\Gamma with vertex set [S][n][S]\subset[n] define the graph Γ/S\Gamma/S with the vertex set [n][S][n]\setminus[S] having an edge between two vertices ii and jj whenever they are path-connected in the restriction Γ|S{i,j}.\Gamma|_{S\cup\{i,j\}}. Then the boundary of the graph-associahedron PΓP_{\Gamma} can be given combinatorially as

(2) dPΓ=S(Γ),SΓPS×PΓ/S.dP_{\Gamma}=\sum_{S\in\mathcal{B}(\Gamma),\,S\neq\Gamma}P_{S}\times P_{\Gamma/S}.
Proposition 2.2.

For any connected graph Γ\Gamma the corresponding graph-assiciahedron PΓP_{\Gamma} is indecomposable.

Proof.

We can prove this by induction in the number kk of the vertices of Γ\Gamma. For k=1,2k=1,2 this is obvious. Let this be true for all k<nk<n and consider PΓP_{\Gamma} for Γ\Gamma with nn vertices. Assume that PΓ=P1×P2P_{\Gamma}=P_{1}\times P_{2} is decomposable, so that the boundary dPΓ=dP1×P2+P1×dP2.dP_{\Gamma}=dP_{1}\times P_{2}+P_{1}\times dP_{2}.

On the other hand the boundary can be given by formula (2), so

dP1×P2+P1×dP2=S(Γ),SΓPS×PΓ/S.dP_{1}\times P_{2}+P_{1}\times dP_{2}=\sum_{S\in\mathcal{B}(\Gamma),\,S\neq\Gamma}P_{S}\times P_{\Gamma/S}.

By induction in the right-hand side we have the products of the indecomposable polytopes, containing PS×PΓ/S=PΓ/SP_{S}\times P_{\Gamma/S}=P_{\Gamma/S} for any one-vertex SS. From the uniqueness part of Theorem 2.1 it follows that either P1P_{1} or P2P_{2} must be a segment, which easily leads to a contradiction. ∎

There are two famous particular cases of graph-associahedra: associahedra 𝔞n\mathfrak{a}_{n} (or Stasheff polytopes, traditionally denoted as Kn+1K_{n+1}) and permutohedra πn\pi_{n}, corresponding to path and complete graphs with nn vertices respectively.

For the associahedron 𝔞n\mathfrak{a}_{n} the corresponding graph Γ\Gamma is a path with edges (k,k+1), 1kn(k,k+1),\,1\leq k\leq n. The connected subgraphs SS are strings [i,j]=(i,i+1,,j1,j)[i,j]=(i,i+1,\dots,j-1,j), which can be naturally labelled by the brackets in the product x1x2xi1(xixj)xj+1xn+1x_{1}x_{2}\dots x_{i-1}(x_{i}\dots x_{j})x_{j+1}\dots x_{n+1}, linking this with the original Stasheff’s description (see [49, 61] for the history of this remarkable polytope and its role in combinatorics and algebra). The number of vertices of 𝔞n\mathfrak{a}_{n} is the Catalan number Cn=1n+1(2nn).C_{n}=\frac{1}{n+1}{2n\choose n}. One can also describe 𝔞n\mathfrak{a}_{n} also as the Newton polytope of the polynomial 1i<jn(xi+xi+1++xj1+xj).\prod_{1\leq i<j\leq n}(x_{i}+x_{i+1}+\dots+x_{j-1}+x_{j}). For n=3n=3 we have pentagon, for n=4n=4 3D associahedron, which is shown on Fig. 1 together with its realisation as 2-truncated cube found in [14].

Permutohedron πn\pi_{n} corresponds to the complete graph with nn vertices and can be described as the convex hull of the points σ(ρ)n,σSn\sigma(\rho)\in\mathbb{R}^{n},\,\sigma\in S_{n} with ρ=(1,2,,n),\rho=(1,2,\dots,n), or as the Newton polytope of the Vandermonde polynomial 1i<jn(xixj).\prod_{1\leq i<j\leq n}(x_{i}-x_{j}). For n=3n=3 we have hexagon, for n=4n=4 - the truncated octahedron shown on Fig. 2.

Refer to caption
Figure 2. 3D permutohedron π4\pi_{4}.

Our first observation is that these two sequences of combinatorial polytopes can be characterised in terms of the differential algebra of polytopes as follows.

Let Pk,k0P_{k},\,k\geq 0 be a sequence of simple polytopes with the dimension of PkP_{k} equal to k.k. Consider the following question: when the corresponding subring of 𝔖\mathfrak{S} generated by these polytopes is closed under differentiation dd (and thus is a differential subring of 𝔖\mathfrak{S})?

Assuming that Pk=PΓk+1P_{k}=P_{\Gamma_{k+1}} are graph-associahedra of the connected graphs Γk+1\Gamma_{k+1} with k+1k+1 vertices, we can give the following answer.

Theorem 2.3.

A sequence of graph-associahedra Pk,k0P_{k},\,k\geq 0 generates a subring of 𝔖\mathfrak{S} closed under differentiation dd only in two cases: associahedra Pk=𝔞k+1P_{k}=\mathfrak{a}_{k+1} and permutohedra Pk=πk+1.P_{k}=\pi_{k+1}.

Proof.

Let Γn,n1\Gamma_{n},\,n\geq 1 be the sequence of connected graphs, such that the corresponding graph-associahedra Pk=PΓk+1,k0P_{k}=P_{\Gamma_{k+1}},\,k\geq 0 (which are indecomposable by Proposition 2.2) freely generate a polynomial subring of 𝔖\mathfrak{S}, which is closed under differentiation d.d.

We need the following result from graph theory.

Lemma 2.4.

Let Γn,n\Gamma_{n},\,n\in\mathbb{N} be a sequence of connected graphs with nn vertices, such that any connected induced subgraph SS of Γn\Gamma_{n} and the corresponding graph Γn/S\Gamma_{n}/S are equivalent to some of Γj\Gamma_{j} with j<n.j<n. Then this must be the sequence of paths or complete graphs.

Proof.

We will prove this by induction in n.n. For n3n\leq 3 the claim is obvious since Γ3\Gamma_{3} is either a path or complete graph K3K_{3} with 3 vertices. Assume now that the claim is true for nkn\leq k and prove it for n=k+1.n=k+1.

Assume first that all induced subgraphs of Γ=Γn\Gamma=\Gamma_{n} are paths. Then Γ\Gamma must be either a tree, or a cycle. Indeed, if Γ\Gamma is not a tree, it must contain a cycle CΓn,C\subseteq\Gamma_{n}, so the corresponding induced subgraph C¯\bar{C} is not a path. The only exception is when C=Γ,C=\Gamma, so Γ\Gamma is a cycle.

Let Γ\Gamma be a tree. If Γ\Gamma has a vertex VV of index m3,m\geq 3, then Γ/V\Gamma/V contains a complete graph KmK_{m}, which is impossible. So all vertices of Γ\Gamma have index 1 or 2 and Γ\Gamma is a path. If Γ\Gamma is a cycle with n>3n>3 vertices, then for its subpath SS with n3n-3 vertices the graph Γ/S=K3\Gamma/S=K_{3}, which contradicts the assumption. Thus Γn\Gamma_{n} must be a path in this case.

Assume now that all induced subgraphs of Γ\Gamma are complete graphs. If Γ\Gamma has two vertices V1V_{1} and V2V_{2} which are not connected by an edge, then the shortest path connecting them is an induced subgraph, which is not complete. Thus Γn\Gamma_{n} must be a complete graph in this case. ∎

The theorem now follows from Lemma and the boundary formula (2) since the paths and complete graphs correspond to the associahedra and permutohedra respectively. ∎

Remark 2.5.

Note that there are other sequences of the simple polytopes PkP_{k} generating a subring in 𝔖\mathfrak{S} invariant under differentiation (e.g. simplices Δk\Delta_{k}, cubes k\Box_{k}), but the general problem of description of all such sequences seems to be open. We conjecture that in the indecomposable case the answer is given by the sequences of simplices, associahedra and permutohedra.

3. PDEs for graph-associahedra and inversion formulas

We start with the differential equations for the generating functions of some special polytopes, derived by one of the authors [11] (see also [10] and [15], section 1.8).

Let PP be a convex polytope of dimension nn and k(P)\mathcal{F}_{k}(P) be the set of all faces of PP of codimension kk considered as elements of 𝔓\mathfrak{P}. Consider the corresponding generating function FP(t)𝔓[t]F_{P}(t)\in\mathfrak{P}[t] defined by

(3) FP(t)=k=0ntkfk(P)f=f(P)ftcodimf.F_{P}(t)=\sum_{k=0}^{n}t^{k}\sum_{f\in\mathcal{F}_{k}(P)}f=\sum_{f\in\mathcal{F}(P)}f\,t^{\text{codim}f}.

One can check that the valuation map Ft0:𝔓𝔓F_{t_{0}}:\mathfrak{P}\to\mathfrak{P} sending PP to FP(t0)F_{P}(t_{0}) is a ring-homomorphism for all t0t_{0}.

Define now the motivic interior PP^{\circ} of a polytope PP as the alternating sum of all the faces of P:P:

(4) P:=FP(1)=f(P)(1)codimff𝔓.P^{\circ}:=F_{P}(-1)=\sum_{f\in\mathcal{F}(P)}(-1)^{\text{codim}f}f\in\mathfrak{P}.

At the level of the set-theoretical characteristic functions this formula agrees with a well-known result in the theory of finitely additive measures of convex polytopes [45, 52] with P=relint(P)P^{\circ}=relint(P) being the standard relative interior [61] of the polytope PP. Note that the polytope PP can be represented simply as

P=f(P)relint(f),P=\sqcup_{f\in\mathcal{F}(P)}relint(f),

where \sqcup denotes the disjoint union of sets (see [61]). More justification for the terminology is given by the theory of toric varieties, see Section 5 below.

For simple polytopes we have the following formula, which will be important for us.

Lemma 3.1.

For any simple convex polytope PP we have

(5) FP(t)=etdP:=k0tkk!dkP,F_{P}(t)=e^{td}P:=\sum_{k\geq 0}\frac{t^{k}}{k!}d^{k}P,

which can be used to characterise the simple polytopes among all polytopes.

The proof easily follows from the simplicity of the polytope. In particular, for simple polytope PP the motivic interior can be written as P=edP.P^{\circ}=e^{-d}P.

Following [11], consider now the generating functions of FP(t)F_{P}(t) for the families of associahedra and permutohedra

(6) 𝒜(t,x):=n=1F𝔞n(t)xn+1=𝔞1x2+(𝔞2+2𝔞1t)x3+(𝔞3+5𝔞2t+5𝔞1t2)x4+,\displaystyle\mathcal{A}(t,x):=\sum_{n=1}^{\infty}F_{\mathfrak{a}_{n}}(t)x^{n+1}=\mathfrak{a}_{1}x^{2}+(\mathfrak{a}_{2}+2\mathfrak{a}_{1}t)x^{3}+(\mathfrak{a}_{3}+5\mathfrak{a}_{2}t+5\mathfrak{a}_{1}t^{2})x^{4}+\dots,
𝒫(t,x):=n=1Fπn(t)xnn!=π1x+(π2+2π1t)x22!+(π3+6π2t+6π1t2)x33!+.\displaystyle\mathcal{P}(t,x):=\sum_{n=1}^{\infty}F_{\pi_{n}}(t)\frac{x^{n}}{n!}=\pi_{1}x+(\pi_{2}+2\pi_{1}t)\frac{x^{2}}{2!}+(\pi_{3}+6\pi_{2}t+6\pi_{1}t^{2})\frac{x^{3}}{3!}+\dots.

Here both 𝔞1\mathfrak{a}_{1} and π1\pi_{1} are points and thus are units in the polytopal ring, but for our purposes we will keep track of them considering them as independent variables. In particular, we define the corresponding weighted motivic interiors as

(7) 𝔞n:=F𝔞n(𝔞1)=f(𝔞n)(𝔞1)codimff,\displaystyle\mathfrak{a}_{n}^{\circ}:=F_{\mathfrak{a}_{n}}(-\mathfrak{a}_{1})=\sum_{f\in\mathcal{F}(\mathfrak{a}_{n})}(-\mathfrak{a}_{1})^{\text{codim}f}f,
πn:=Fπn(π1)=f(πn)(π1)codimff.\displaystyle\pi_{n}^{\circ}:=F_{\pi_{n}}(-\pi_{1})=\sum_{f\in\mathcal{F}(\pi_{n})}(-\pi_{1})^{\text{codim}f}f.

For example, for n=3n=3 we have 𝔞3=𝔞35𝔞2𝔞1+5𝔞13,π3=π36π2π1+6π13.\mathfrak{a}_{3}^{\circ}=\mathfrak{a}_{3}-5\mathfrak{a}_{2}\mathfrak{a}_{1}+5\mathfrak{a}_{1}^{3},\,\,\pi_{3}^{\circ}=\pi_{3}-6\pi_{2}\pi_{1}+6\pi_{1}^{3}. The corresponding motivic interiors (4) can be found from the weighted motivic interiors by setting 𝔞1=π1=1.\mathfrak{a}_{1}=\pi_{1}=1.

Denote by

ut(t,x):=tu(t,x),ux(t,x):=xu(t,x)u_{t}(t,x):=\frac{\partial}{\partial t}u(t,x),\,\,u_{x}(t,x):=\frac{\partial}{\partial x}u(t,x)

the corresponding partial derivatives of the function u(t,x).u(t,x).

Theorem 3.2 ([11, 15]).

The generating functions (6) satisfy the partial differential equations in 𝔓[t][[x]]\mathfrak{P}[t][[x]]

(8) 𝔞1𝒜t=𝒜𝒜x,𝒜(0,x)=A(x):=n=1𝔞nxn+1,\mathfrak{a}_{1}\mathcal{A}_{t}=\mathcal{A}\mathcal{A}_{x},\quad\mathcal{A}(0,x)=A(x):=\sum_{n=1}^{\infty}\mathfrak{a}_{n}x^{n+1},
(9) π1𝒫t=𝒫2,𝒫(0,x)=P(x):=n=1πnxnn!.\pi_{1}\mathcal{P}_{t}=\mathcal{P}^{2},\quad\mathcal{P}(0,x)=P(x):=\sum_{n=1}^{\infty}\pi_{n}\frac{x^{n}}{n!}.

Indeed, the boundary formula (2) in this particular case gives

(10) d𝔞n=i=1n1(i+1)𝔞i×𝔞ni,dπn=i=1n1(ni)πi×πni,d\mathfrak{a}_{n}=\sum_{i=1}^{n-1}(i+1)\mathfrak{a}_{i}\times\mathfrak{a}_{n-i},\quad d\pi_{n}=\sum_{i=1}^{n-1}{n\choose i}\pi_{i}\times\pi_{n-i},

which together with Lemma leads to the differential equations (8), (9).

The partial differential equation

(11) ut=uuxu_{t}=uu_{x}

is called Hopf equation (or inviscid Burgers’ equation) and used as the simplest model to describe the breaking of waves phenomenon.

It is well-known that its solution u(t,x)u(t,x) with initial value u(0,x)=f(x)u(0,x)=f(x) can be given implicitly by the formula

(12) u=f(x+ut),u=f(x+ut),

which can be easily checked directly. Applying this now for the equation (8), we have the identity

(13) 𝒜(𝔞1t,x)=A(x+𝒜(𝔞1t,x)t).\mathcal{A}(\mathfrak{a}_{1}t,x)=A(x+\mathcal{A}(\mathfrak{a}_{1}t,x)t).

Substituting here t=1t=-1 we see that the generating function

(14) A(x):=𝒜(𝔞1,x)=n=1𝔞nxn+1A^{\circ}(x):=\mathcal{A}(-\mathfrak{a}_{1},x)=\sum_{n=1}^{\infty}\mathfrak{a}^{\circ}_{n}x^{n+1}

satisfies the relation

(15) A(x)=A(xA(x)).A^{\circ}(x)=A(x-A^{\circ}(x)).

Consider the series ξ=xA(x)\xi=x-A^{\circ}(x), then

x=ξ+A(x)=ξ+A(ξ).x=\xi+A^{\circ}(x)=\xi+A(\xi).

Thus we have proved the following theorem, which is an interpretation of the result of Loday [41] (see also [5]).

Theorem 3.3.

The generating power series of associahedra and their weighted motivic interiors

(16) α(x)=x+n=1𝔞nxn+1,β(x)=xn=1𝔞nxn+1\alpha(x)=x+\sum_{n=1}^{\infty}\mathfrak{a}_{n}x^{n+1},\quad\beta(x)=x-\sum_{n=1}^{\infty}\mathfrak{a}^{\circ}_{n}x^{n+1}

are compositionally inverse to each other in the polytopal ring 𝔓[[x]]\mathfrak{P}[[x]]:

α(β(x))x,β(α(x))x.\alpha(\beta(x))\equiv x,\quad\beta(\alpha(x))\equiv x.

Note that since by Theorem 2.1 and Proposition 2.2 the associahedra are algebraically independent in 𝔓\mathfrak{P}, this implies the following universal inversion formulae: the coefficients of the formal series α(x)=x+n1anxn+1\alpha(x)=x+\sum_{n\geq 1}a_{n}x^{n+1} and its compositional inverse β(x)=x+n1bnxn+1\beta(x)=x+\sum_{n\geq 1}b_{n}x^{n+1} are related by the following formulae, expressing each other as certain polynomials with integer coefficients:

b1=a1,b2=a2+2a12,b3=a3+5a1a25a13,b_{1}=-a_{1},\quad b_{2}=-a_{2}+2a_{1}^{2},\quad b_{3}=-a_{3}+5a_{1}a_{2}-5a_{1}^{3},
b4=a4+6a1a3+3a2221a12a2+14a14,b_{4}=-a_{4}+6a_{1}a_{3}+3a_{2}^{2}-21a_{1}^{2}a_{2}+14a_{1}^{4},

in agreement with the formulas for the weighted motivic interiors of the associahedra

𝔞1=𝔞1,𝔞2=𝔞22𝔞12,𝔞3=𝔞35𝔞1𝔞2+5𝔞13,\mathfrak{a}^{\circ}_{1}=\mathfrak{a}_{1},\mathfrak{a}^{\circ}_{2}=\mathfrak{a}_{2}-2\mathfrak{a}_{1}^{2},\quad\mathfrak{a}^{\circ}_{3}=\mathfrak{a}_{3}-5\mathfrak{a}_{1}\mathfrak{a}_{2}+5\mathfrak{a}_{1}^{3},
𝔞4=𝔞46𝔞1𝔞33𝔞22+21𝔞12𝔞214𝔞14.\mathfrak{a}^{\circ}_{4}=\mathfrak{a}_{4}-6\mathfrak{a}_{1}\mathfrak{a}_{3}-3\mathfrak{a}_{2}^{2}+21\mathfrak{a}_{1}^{2}\mathfrak{a}_{2}-14\mathfrak{a}_{1}^{4}.

Note that the faces of 𝔞n\mathfrak{a}_{n} of codimension dd are in one-to-one correspondence with dissections of a based (n+3)(n+3)-gon by dd non-intersecting diagonals (see [30, 25]). Alternatively, the faces of 𝔞n\mathfrak{a}_{n} can be labelled by non-isomorphic planar rooted trees with n+1n+1 leaves [23], which can be used to rewrite the inversion formulas in these terms (see [56]).

The corresponding polynomials can also be expressed in terms of the partial Bell polynomials Bn,kB_{n,k} (see e.g. [21], Section 3.8) or, more explicitly, as

(17) bn=1(n+1)!k1+2k2+=n,ki0(1)k(n+k)!k1!k2!a1k1a2k2,b_{n}=\frac{1}{(n+1)!}\sum_{k_{1}+2k_{2}+\dots=n,\,k_{i}\geq 0}(-1)^{k}\frac{(n+k)!}{k_{1}!k_{2}!\dots}a_{1}^{k_{1}}a_{2}^{k_{2}}\dots,

where k=k1+k2+k=k_{1}+k_{2}+\dots (see e.g. [31]).

Similarly, the corresponding permutohedral differential equation (9) has the explicit solution

(18) 𝒫(π1t,x)=π(x)1tπ(x).\mathcal{P}(\pi_{1}t,x)=\frac{\pi(x)}{1-t\pi(x)}.

Substituting here t=1t=-1 we have the relation 𝒫(π1,x)=π(x)1+π(x).\mathcal{P}(-\pi_{1},x)=\frac{\pi(x)}{1+\pi(x)}. This implies that 1𝒫(π1,x)=11+π(x)1-\mathcal{P}(-\pi_{1},x)=\frac{1}{1+\pi(x)} and the following result, which is an interpretation of the known result from [5].

Theorem 3.4.

The exponential generating series of permutohedra and their weighted motivic interiors

(19) P(x)=1+n=1πnxnn!,P(x)=1n=1πnxnn!P(x)=1+\sum_{n=1}^{\infty}\pi_{n}\frac{x^{n}}{n!},\quad P^{\circ}(x)=1-\sum_{n=1}^{\infty}\pi_{n}^{\circ}\frac{x^{n}}{n!}

is multiplicatively inverse to each other in the polytopal ring 𝔓[[x]]\mathfrak{P}[[x]]:

P(x)P(x)1.P(x)P^{\circ}(x)\equiv 1.

Again since permotohedra are algebraically independent in 𝔓,\mathfrak{P}, this implies the universal multiplicative inversion formulae for the formal power series. Namely, the coefficients of the series p(x)=1+n=1pnxnn!p(x)=1+\sum_{n=1}^{\infty}p_{n}\frac{x^{n}}{n!} and its multiplicative inverse q(x)=1+n=1qnxnn!q(x)=1+\sum_{n=1}^{\infty}q_{n}\frac{x^{n}}{n!} are related by the polynomial formulas with integer coefficients: q1=p1,q2=p2+2p12,q_{1}=-p_{1},\,q_{2}=-p_{2}+2p_{1}^{2},

q3=p3+6p1p26p13,q4=p4+8p1p3+6p2236p12p2+24p14,q_{3}=-p_{3}+6p_{1}p_{2}-6p_{1}^{3},\,\,\,q_{4}=-p_{4}+8p_{1}p_{3}+6p_{2}^{2}-36p_{1}^{2}p_{2}+24p_{1}^{4},

in agreement with the fact that π3\pi_{3} is hexagon and π4\pi_{4} has 8 hexagonal and 6 quadrilateral faces, 36 edges and 24 vertices (see Fig. 2). This is also in a good agreement with the formulas for the weighted motivic interiors of the permutohedra: π1=p1,π2=p2+2p12,\pi^{\circ}_{1}=-p_{1},\,\pi^{\circ}_{2}=-p_{2}+2p_{1}^{2},

(20) π3=π36π1π2+6π13,π4=π48π1π36π22+36π12π224π14.\pi^{\circ}_{3}=\pi_{3}-6\pi_{1}\pi_{2}+6\pi_{1}^{3},\,\,\,\pi^{\circ}_{4}=\pi_{4}-8\pi_{1}\pi_{3}-6\pi_{2}^{2}+36\pi_{1}^{2}\pi_{2}-24\pi_{1}^{4}.

4. Inversion formula and the Deligne-Mumford moduli spaces

There is an interesting link with the moduli space M¯n=M¯0,n\bar{M}_{n}=\bar{M}_{0,n} of stable genus zero curves with nn ordered marked points. Its real version M¯n+1()\bar{M}_{n+1}(\mathbb{R}) is known to be tessellated by n!/2n!/2 copies of associahedron KnK_{n} (see [36, 23]). This explains the result of McMullen, who interpreted the compositional inversion formulas in terms of the corresponding real strata [44].

It is interesting that the geometry of the complex version M¯n+1=M¯n+1()\bar{M}_{n+1}=\bar{M}_{n+1}(\mathbb{C}) allows also to describe the compositional inversion of the exponential power series. More precisely, McMullen proved that the compositional inverse of the formal series f(x)=xn1anxn+1(n+1)!f(x)=x-\sum_{n\geq 1}a_{n}\frac{x^{n+1}}{(n+1)!} is given by

(21) g(x)=x+n1bnxn+1(n+1)!,bn=Nn1,,nsan1ans,g(x)=x+\sum_{n\geq 1}b_{n}\frac{x^{n+1}}{(n+1)!},\quad b_{n}=\sum N_{n_{1},\dots,n_{s}}a_{n_{1}}\dots a_{n_{s}},

where Nn1,,nsN_{n_{1},\dots,n_{s}} is the number of strata SM¯n()S\subset\bar{M}_{n}(\mathbb{C}) isomorphic to

Mn1()××Mns()M_{n_{1}}(\mathbb{C})\times\dots\times M_{n_{s}}(\mathbb{C})

(see Theorem 1 in [44]). The proof is based on the well-known bijection between the strata and the marked stable rooted trees, known also as modular graphs [40].

One can combine this with our result [16] claiming that the generating functions of the cobordism classes of the complex projective spaces Pn\mathbb{C}P^{n} and the theta-divisors Θn\Theta^{n}

α(z)=z+n=1[Pn]zn+1n+1,β(z)=z+n=1[Θn]zn+1(n+1)!,\alpha(z)=z+\sum_{n=1}^{\infty}[\mathbb{C}P^{n}]\frac{z^{n+1}}{n+1},\quad\beta(z)=z+\sum_{n=1}^{\infty}[\Theta^{n}]\frac{z^{n+1}}{(n+1)!},

are compositionally inverse to each other, to express the cobordism class [Θn1][\Theta^{n-1}] in terms of the cobordism classes [Pk],kn1[\mathbb{C}P^{k}],\,k\leq n-1 as follows

(22) [Θn]=Nn1,,nsan1ans,ak=k![Pk].[\Theta^{n}]=\sum N_{n_{1},\dots,n_{s}}a_{n_{1}}\dots a_{n_{s}},\quad a_{k}=-k![\mathbb{C}P^{k}].

One more remarkable fact here is due to Getzler [32]. Recall that M¯n\bar{M}_{n} is a special (Deligne-Mumford) compactification of the moduli space MnM_{n} of nn ordered distinct points on complex projective line considered up to projective equivalence. Using the framework of the operads, Getzler proved that the exponential generating functions of the corresponding (in case of MnM_{n}, motivic) Poincare polynomials

F(t,x):=xn2PMn+1(t)xnn!,G(t,y):=y+n2PM¯n+1(t)ynn!F(t,x):=x-\sum_{n\geq 2}P_{M_{n+1}}(t)\frac{x^{n}}{n!},\quad G(t,y):=y+\sum_{n\geq 2}P_{\bar{M}_{n+1}}(t)\frac{y^{n}}{n!}

are inverse to each other: F(G(t,y),t)y,G(F(t,x),t)x.F(G(t,y),t)\equiv y,\,\,G(F(t,x),t)\equiv x. Since the motivic Poincare polynomial

(23) PMn(t)=(t22)(t23)(t2n+2)P_{M_{n}}(t)=(t^{2}-2)(t^{2}-3)\dots(t^{2}-n+2)

is known explicitly (see e.g. [37]), this determines PM¯n(t).P_{\bar{M}_{n}}(t). Substituting here t=1t=-1 we have the same relation between exponential generating functions of the Euler characteristics

F(x):=xn2χ(Mn+1)xnn!,G(y):=y+n2χ(M¯n+1)ynn!.F(x):=x-\sum_{n\geq 2}\chi(M_{n+1})\frac{x^{n}}{n!},\quad G(y):=y+\sum_{n\geq 2}\chi(\bar{M}_{n+1})\frac{y^{n}}{n!}.

Since χ(Mn+1)=(1)n(n2)!\chi(M_{n+1})=(-1)^{n}(n-2)! we see that G(y)G(y) is the compositional inverse of the series F(x)=xn2(1)nn(n1)xn,F(x)=x-\sum_{n\geq 2}\frac{(-1)^{n}}{n(n-1)}x^{n}, which allows to compute χ(M¯n),\chi(\bar{M}_{n}), giving OEIS sequence A074059: 2, 7, 34, 213, 1630, 14747, 153946, 18214732,\,7,\,34,\,213,\,1630,\,14747,\,153946,\,1821473...

It is interesting that for the real version of Deligne-Mumford spaces M¯n()\bar{M}_{n}(\mathbb{R}) McMullen proved that the same relation

Bn=Nn1,,nsAn1AnsB_{n}=\sum N_{n_{1},\dots,n_{s}}A_{n_{1}}\dots A_{n_{s}}

describes the functional inversion G(x)=xn1Bnxn+1G(x)=x-\sum_{n\geq 1}B_{n}x^{n+1} of the usual power series F(x)=x+n1Anxn+1F(x)=x+\sum_{n\geq 1}A_{n}x^{n+1} (see Corollary 5 in [44]).

Over finite fields one can interpret these results within the approach of Weil and Deligne as follows.

Let 𝔽p\mathbb{F}_{p} with pp prime be the finite field with pp elements and M0,n(𝔽p),M0,n(𝔽p)¯M_{0,n}(\mathbb{F}_{p}),\overline{M_{0,n}(\mathbb{F}_{p})} be the corresponding varieties over 𝔽p.\mathbb{F}_{p}. Let |M0,n(𝔽p)||M_{0,n}(\mathbb{F}_{p})| and |M0,n(𝔽p)¯||\overline{M_{0,n}(\mathbb{F}_{p})}| be the number of points in these varieties respectively.

Theorem 4.1.

The generating power series of the number of points in M0,n(𝔽p)M_{0,n}(\mathbb{F}_{p}) and M0,n(𝔽p)¯\overline{M_{0,n}(\mathbb{F}_{p})}

(24) α(x)=x+n=2|M0,n+1(𝔽p)¯|xnn!,β(x)=xn=2|M0,n+1(𝔽p)|xnn!\alpha(x)=x+\sum_{n=2}^{\infty}|\overline{M_{0,n+1}(\mathbb{F}_{p})}|\frac{x^{n}}{n!},\quad\beta(x)=x-\sum_{n=2}^{\infty}|M_{0,n+1}(\mathbb{F}_{p})|\frac{x^{n}}{n!}

are compositionally inverse to each other for all primes p.p.

Proof.

The number of points in M0,n(𝔽p)M_{0,n}(\mathbb{F}_{p}) is easy to compute directly:

|M0,n(𝔽p)|=(p2)(p3)(pn+2),n4,|M_{0,n}(\mathbb{F}_{p})|=(p-2)(p-3)\dots(p-n+2),\quad n\geq 4,

with |M0,3(𝔽p)|=1.|M_{0,3}(\mathbb{F}_{p})|=1.

The number of points in M0,n(𝔽p)¯\overline{M_{0,n}(\mathbb{F}_{p})} was computed by Amburg, Kreines and Shabat in [3], who proved, in particular, that this number coincides with the value of the Poincare polynomial PM0,n()¯(t)P_{\overline{M_{0,n}(\mathbb{C})}}(t) when p=t2p=t^{2}:

|M0,n(𝔽p)¯|=PM0,n()¯(t),p=t2.|\overline{M_{0,n}(\mathbb{F}_{p})}|=P_{\overline{M_{0,n}(\mathbb{C})}}(t),\quad p=t^{2}.

Since motivic Poincare polynomial of M0,n()M_{0,n}(\mathbb{C}) is given by (23), we have the same relation for M0,n()M_{0,n}(\mathbb{C}): |M0,n(𝔽p)|=PM0,n()(t),p=t2.|M_{0,n}(\mathbb{F}_{p})|=P_{M_{0,n}(\mathbb{C})}(t),\quad p=t^{2}. Now the proof follows from Getzler’s result. ∎

5. Inversion formula and permutohedral varieties

The permutohedral variety XΠnX_{\Pi}^{n} is the toric variety [22], corresponding to the nn-dimensional permutohedron πn+1\pi_{n+1}. In particular, XΠ1=P1X_{\Pi}^{1}=\mathbb{C}P^{1}, XΠ2X_{\Pi}^{2} is the degree 6 del Pezzo surface.

Every toric variety XnX^{n} is a closure of the algebraic torus 𝕋n=()n\mathbb{T}_{\mathbb{C}}^{n}=(\mathbb{C}^{*})^{n}, which can be defined as the motivic interior of XnX^{n}.

To justify this consider the Grothendieck ring of complex quasi-projective varieties K0(𝒱)K_{0}(\mathcal{V}_{\mathbb{C}}) generated by the isomorphism classes of complex quasi-projective varieties modulo the relations

[X]=[Y]+[XY],[X]=[Y]+[X\setminus Y],

where YY is a Zariski locally closed subset of XX, with multiplication given by the formula [X]·[Y]=[X×Y][X]\textperiodcentered[Y]=[X\times Y] (see e.g. [33]).

Let XPX_{P} be the smooth toric variety corresponding to a Delzant polytope PP and YfXPY_{f}\subset X_{P} be the toric subvariety corresponding to the face f(P).f\in\mathcal{F}(P).

The following result provides more justification for our definition of the motivic interior. In analogy with (4) define the motivic interior of XPX_{P} as

(25) [XP]:=f(P)(1)codimf[Yf]K0(𝒱𝒞).[X_{P}^{\circ}]:=\sum_{f\in\mathcal{F}(P)}(-1)^{\text{codim}f}[Y_{f}]\in K_{0}(\mathcal{V}_{\mathcal{C}}).
Proposition 5.1.

In the Grothendieck ring K0(𝒱)K_{0}(\mathcal{V}_{\mathbb{C}}) we have the relation

(26) f(P)(1)codimf[Yf]=[𝕋n].\sum_{f\in\mathcal{F}(P)}(-1)^{\text{codim}f}[Y_{f}]=[\mathbb{T}_{\mathbb{C}}^{n}].

In other words, the motivic interior of toric variety XPX_{P} is [XP]=[𝕋n].[X_{P}^{\circ}]=[\mathbb{T}_{\mathbb{C}}^{n}].

The proof follows from the Proposition-definition 6 in [52] reformulated in terms of toric geometry [22] (see also Theorem 3.2.4 in [48]).

Note that similar relation holds for all convex polytopes if we extend it to the framework of the toric topology [15].

This implies that Theorem 3.5 can be reformulated as follows.

Theorem 5.2.

The exponential generating series of permutohedral varieties and algebraic tori

(27) π(x):=1+n=1[XΠn]xnn!,τ(x):=1n=1[𝕋n1]xnn!\pi(x):=1+\sum_{n=1}^{\infty}[X_{\Pi}^{n}]\frac{x^{n}}{n!},\quad\tau(x):=1-\sum_{n=1}^{\infty}[\mathbb{T}_{\mathbb{C}}^{n-1}]\frac{x^{n}}{n!}

are multiplicatively inverse to each other in the ring K0(𝒱)[[x]].K_{0}(\mathcal{V}_{\mathbb{C}})[[x]].

Motivic Poincare polynomial PZ(t),ZK0(𝒱)P_{Z}(t),\,Z\in K_{0}(\mathcal{V}_{\mathbb{C}}) defines the ring homomorphism K0(𝒱)[t],K_{0}(\mathcal{V}_{\mathbb{C}})\to\mathbb{Z}[t], where for the projective varieties X,YXX,Y\subseteq X by definition

P[XY](t)=PX(t)PY(t),P_{[X\setminus Y]}(t)=P_{X}(t)-P_{Y}(t),

where PX(t)P_{X}(t) is the usual Poincare polynomial. In particular, since \mathbb{C}^{*} is the complement of two points in P1\mathbb{C}P^{1}, the motivic Poincare polynomial P[](t)=PP1(t)2=(t2+1)2=t21,P_{[\mathbb{C}^{*}]}(t)=P_{\mathbb{C}P^{1}}(t)-2=(t^{2}+1)-2=t^{2}-1, and thus

P𝕋n(t)=(t21)n.P_{\mathbb{T}_{\mathbb{C}}^{n}}(t)=(t^{2}-1)^{n}.

As a corollary we have the following well-known fact (see e.g. [50]) relating permutohedral varieties with the classical Eulerian polynomials [19]. These polynomials were introduced by Euler in 1755 by the relation

k=1kntn=tAn(t)(1t)n+1.\sum_{k=1}^{\infty}k^{n}t^{n}=\frac{tA_{n}(t)}{(1-t)^{n+1}}.

They have the generating function

(28) n0An(s)xnn!=s1se(s1)x\sum_{n\geq 0}A_{n}(s)\frac{x^{n}}{n!}=\frac{s-1}{s-e^{(s-1)x}}

and can be computed recursively by

An+1(t)=[t(1t)ddt+nt+1]An(t),A0=A1=1:A_{n+1}(t)=[t(1-t)\frac{d}{dt}+nt+1]A_{n}(t),\quad A_{0}=A_{1}=1:
A1=1,A2=s+1,A3=s2+4s+1,A4=s3+11s2+11s+1,A_{1}=1,\,\,A_{2}=s+1,\,\,A_{3}=s^{2}+4s+1,\,\,A_{4}=s^{3}+11s^{2}+11s+1,
A5=s4+26s3+66s2+26s+1,A6=s5+57s4+302s3+302s2+57s+1.A_{5}=s^{4}+26s^{3}+66s^{2}+26s+1,\,\,A_{6}=s^{5}+57s^{4}+302s^{3}+302s^{2}+57s+1.

Their coefficients are known as Eulerian numbers and have natural combinatorial interpretations (see [56]).

Corollary 5.3.

The Poincare polynomial of permutohedral varieties is

(29) PXΠn(t)=An(t2),P_{X_{\Pi}^{n}}(t)=A_{n}(t^{2}),

where An(s)A_{n}(s) are the Eulerian polynomials.

Proof.

The generating function of the motivic Poincare polynomials of the algebraic tori can be computed explicitly as

n=0P𝕋n(t)xn+1(n+1)!=n=0(t21)nxn+1(n+1)!=e(t21)x1t21.\sum_{n=0}^{\infty}P_{\mathbb{T}_{\mathbb{C}}^{n}}(t)\frac{x^{n+1}}{(n+1)!}=\sum_{n=0}^{\infty}(t^{2}-1)^{n}\frac{x^{n+1}}{(n+1)!}=\frac{e^{(t^{2}-1)x}-1}{t^{2}-1}.

This means that

1n=0P𝕋n(t)xn+1(n+1)!==t2e(t21)xt21,1-\sum_{n=0}^{\infty}P_{\mathbb{T}_{\mathbb{C}}^{n}}(t)\frac{x^{n+1}}{(n+1)!}==\frac{t^{2}-e^{(t^{2}-1)x}}{t^{2}-1},

so from Theorem 5.1 we have

1+n=0PXΠn(t)xn+1(n+1)!=t21t2e(t21)x.1+\sum_{n=0}^{\infty}P_{X_{\Pi}^{n}}(t)\frac{x^{n+1}}{(n+1)!}=\frac{t^{2}-1}{t^{2}-e^{(t^{2}-1)x}}.

Comparing this with (28), we have the claim. ∎

Remark. It is interesting to mention that the permotuhedral variety XΠnX_{\Pi}^{n} is isomorphic to the Losev-Manin [42] compactification L¯0,n+3,2\bar{L}_{0,n+3,2} of the moduli space M0,n+3M_{0,n+3} (see more on this in the recent paper [18]).

6. Cyclohedra and Faà di Bruno’s formula

There is another remarkable family of the graph-polyhedra: cyclohedra (or Bott-Taubes polytopes) WnW_{n}, corresponding to the nn-cycle graph. In particular, W3W_{3} is hexagon and W4W_{4} is 3D3D polyhedron shown on Fig. 3.

These polytopes first appeared in Bott and Taubes [9] in connection with the link invariants, although implicitly they were already in the earlier work by Kontsevich [38] (see more detail in [24]).

Refer to caption
Figure 3. 3D cyclohedron W4W_{4}.

Its combinatorics was studied by Simion [55], who had shown, in particular, that the number fk(Wn+1)f_{k}(W_{n+1}) of faces of dimension kk of Wn+1W_{n+1} equals

fk(Wn+1)=(2nknk)(nk).f_{k}(W_{n+1})={2n-k\choose n-k}{n\choose k}.

In particular, Wn+1W_{n+1} has (2nn){2n\choose n} vertices and n(n+1)n(n+1) facets. The corresponding boundary formula dWn=ni=1n1𝔞i×WnidW_{n}=n\sum_{i=1}^{n-1}\mathfrak{a}_{i}\times W_{n-i} implies that the generating function 𝒞(t,x):=n=1FWn(t)xnn,\mathcal{C}(t,x):=\sum_{n=1}^{\infty}F_{W_{n}}(t)\frac{x^{n}}{n}, satisfies the following non-autonomous linear differential equation

𝒞t=𝒜(t,x)𝒞x,𝒞(0,x)=W(x):=n=1Wnxnn,\mathcal{C}_{t}=\mathcal{A}(t,x)\mathcal{C}_{x},\quad\mathcal{C}(0,x)=W(x):=\sum_{n=1}^{\infty}W_{n}\frac{x^{n}}{n},

where 𝒜(t,x)\mathcal{A}(t,x) are the corresponding function for the associahedra (6) and we set for simplicity 𝔞1=W1=1\mathfrak{a}_{1}=W_{1}=1. It is best to combine these two functions together as the solution of the following system of PDEs

(30) 𝒜t=𝒜𝒜x,𝒞t=𝒜𝒞x.\mathcal{A}_{t}=\mathcal{A}\mathcal{A}_{x},\quad\mathcal{C}_{t}=\mathcal{A}\mathcal{C}_{x}.

The solution of this system with the prescribed initial values is given by (15) and

𝒞(t,x)=W(x+t𝒜(t,x)).\mathcal{C}(t,x)=W(x+t\mathcal{A}(t,x)).

Substituting here t=1t=-1 we have the following result. Let WnW^{\circ}_{n} be the motivic interior (4) of the cyclohedron Wn.W_{n}.

Theorem 6.1.

The generating function of the motivic interiors of the cyclohedra

(31) W(x):=n=1Wnnxn=W(β(x))W^{\circ}(x):=\sum_{n=1}^{\infty}\frac{W^{\circ}_{n}}{n}x^{n}=W(\beta(x))

is the result of the substitution into W(x)W(x) the series β(x)\beta(x), which is the compositional inverse of the associahedral series α(x)\alpha(x) given by (16).

The result of the substitution of formal series can be given by the following combinatorial formula attributed to Faà di Bruno (see the history in [35]). It has several important interpretations within the theory of Lie and Hopf algebras and operads [28], as well as relations with integrable systems [54].

For the formal power series Faà di Bruno formula has the following form. Let f(x)=n1anxn,g(x)=n1bnxn,f(x)=\sum_{n\geq 1}a_{n}x^{n},\,\,g(x)=\sum_{n\geq 1}b_{n}x^{n}, then the coefficients of the composition f(g(x))=n1cnxnf(g(x))=\sum_{n\geq 1}c_{n}x^{n} can be written as

(32) cn=k!k1!k2!kn!akj=1nbjkj,c_{n}=\sum{\frac{k!}{k_{1}!\,k_{2}!\,\cdots\,k_{n}!}}a_{k}\prod_{j=1}^{n}b_{j}^{k_{j}},

where k=k1++knk=k_{1}+\dots+k_{n} and the summation is taken over all integer (k1,,kn),ki0(k_{1},\dots,k_{n}),k_{i}\geq 0 such that k1+2k2++nkn=n.k_{1}+2k_{2}+\dots+nk_{n}=n. In particular, we have

c1=a1b1,c2=a1b2+a2b1,c3=a1b3+2a2b1b3+a3b13,c_{1}=a_{1}b_{1},\,\,c_{2}=a_{1}b_{2}+a_{2}b_{1},\,\,c_{3}=a_{1}b_{3}+2a_{2}b_{1}b_{3}+a_{3}b_{1}^{3},
c4=a1b4+a2(b22+2b1b3)+3a3b12b2+a4b14.c_{4}=a_{1}b_{4}+a_{2}(b_{2}^{2}+2b_{1}b_{3})+3a_{3}b_{1}^{2}b_{2}+a_{4}b_{1}^{4}.

Combining this with the inversion formula (17), one can compute the coefficients dnd_{n} of the function

h(x)=f(g(1)(x))=n1dnxn,h(x)=f(g^{(-1)}(x))=\sum_{n\geq 1}d_{n}x^{n},

assuming for convenience that a1=b1=1a_{1}=b_{1}=1:

d1=1,d2=a2b2,d3=a32a2b2+2b22b3,d_{1}=1,\,\,d_{2}=a_{2}-b_{2},\,\,d_{3}=a_{3}-2a_{2}b_{2}+2b_{2}^{2}-b_{3},
d4=a43a3b3+a2(5b222b3)+5b23+5b2b3b4.d_{4}=a_{4}-3a_{3}b_{3}+a_{2}(5b_{2}^{2}-2b_{3})+-5b_{2}^{3}+5b_{2}b_{3}-b_{4}.

The combinatorial formulas for the corresponding modification of the Bell polynomials dn=Dn(a,b)d_{n}=D_{n}(a,b) in terms of the so-called pointed non-crossing partitions can be found in [4] (see Theorem 5.8).

As a corollary we have the following polytopal interpretation of the corresponding polynomials, which was first discovered by Aguiar and Bastidas in relation with the Faà di Bruno Hopf monoid [4] (see Theorem 5.6). Note that both composition and inversion of formal series are naturally embedded into the so-called Faà di Bruno Hopf algebra, playing important role in various problems of mathematics and physics, see [47] and references therein.

Corollary 6.2.

The motivic interior of the cyclohedron WnW_{n} can be written in the algebra of polytopes as

(33) Wn=nDn(a,b),aj=Wjj,bj=𝔞j1.W^{\circ}_{n}=nD_{n}(a,b),\quad a_{j}=\frac{W_{j}}{j},\,b_{j}=\mathfrak{a}_{j-1}.

In particular, we have the formulas

W1=W1,W2=W22𝔞1,W3=W33W2𝔞13𝔞2+6𝔞12,W^{\circ}_{1}=W_{1},\,\,W^{\circ}_{2}=W_{2}-2\mathfrak{a}_{1},\,\,W^{\circ}_{3}=W_{3}-3W_{2}\mathfrak{a}_{1}-3\mathfrak{a}_{2}+6\mathfrak{a}_{1}^{2},
W4=W44W3𝔞1+10W2𝔞124W2𝔞24W1𝔞3+20W1𝔞1𝔞220W1𝔞13.W^{\circ}_{4}=W_{4}-4W_{3}\mathfrak{a}_{1}+10W_{2}\mathfrak{a}_{1}^{2}-4W_{2}\mathfrak{a}_{2}-4W_{1}\mathfrak{a}_{3}+20W_{1}\mathfrak{a}_{1}\mathfrak{a}_{2}-20W_{1}\mathfrak{a}_{1}^{3}.

This is in agreement with the fact that W3W_{3} is a hexagon and W4W_{4} has 4 hexagonal, 4 pentagonal and 4 quadrilateral faces, 30 edges and 20 vertices (see Fig. 3). Again since cyclohedra are algebraically independent, these formulae are universal.

7. Equation and formulas for the stellohedra

There is another interesting graph-associahedron called stellohedron SnS_{n}, corresponding to the star-graph with nn vertices (see Fig.4). In particular, S3=𝔞3S_{3}=\mathfrak{a}_{3} is a pentagon and S4S_{4} is the polyhedron shown on Fig.4. For more details about combinatorics of stellohedra we refer to [11, 50].

Refer to caption
Refer to caption
Figure 4. Star-graph and 3D stellohedron S4S_{4}.

The corresponding boundary formula

dSn=(n1)Sn1+i=1n1(n1i)Si×πnidS_{n}=(n-1)S_{n-1}+\sum_{i=1}^{n-1}{n-1\choose i}S_{i}\times\pi_{n-i}

implies that the corresponding exponential generating functions of stellohedra and permutohedra

𝒮(t,x):=n=0FSn+1(t)xnn!,𝒫(t,x):=n=1Fπn(t)xnn!\mathcal{S}(t,x):=\sum_{n=0}^{\infty}F_{S_{n+1}}(t)\frac{x^{n}}{n!},\quad\mathcal{P}(t,x):=\sum_{n=1}^{\infty}F_{\pi_{n}}(t)\frac{x^{n}}{n!}

satisfy the following system of equations

(34) 𝒫t=𝒫2,𝒮t=(x+𝒫)𝒮.\mathcal{P}_{t}=\mathcal{P}^{2},\quad\mathcal{S}_{t}=(x+\mathcal{P})\mathcal{S}.

For given initial data 𝒫(0,x)=Π(x),𝒮(0,x)=S(x):=n=0Sn+1xnn!\mathcal{P}(0,x)=\Pi(x),\,\mathcal{S}(0,x)=S(x):=\sum_{n=0}^{\infty}S_{n+1}\frac{x^{n}}{n!} this system can be easily solved explicitly as

𝒫(t,x)=Π(x)1tΠ(x),𝒮(t,x)=etxS(x)1tΠ(x)\mathcal{P}(t,x)=\frac{\Pi(x)}{1-t\Pi(x)},\quad\mathcal{S}(t,x)=\frac{e^{tx}S(x)}{1-t\Pi(x)}

(see [11, 15]). Substituting here t=1t=-1 and using theorem 3.5 we have

Theorem 7.1.

The exponential generating function of the motivic interiors of the stellohedra can be expressed via permutohedral generating functions (19) by the formula

(35) S(x):=n=0Sn+1xnn!=exΠ(x)S(x)=exΠ1(x)S(x).S^{\circ}(x):=\sum_{n=0}^{\infty}S^{\circ}_{n+1}\frac{x^{n}}{n!}=e^{-x}\Pi^{\circ}(x)S(x)=e^{-x}\Pi^{-1}(x)S(x).

In particular, using formulas (20) we have

S4=S43S3(1+π1)+3S2(1+2π1+2π12π2)S^{\circ}_{4}=S_{4}-3S_{3}(1+\pi_{1})+3S_{2}(1+2\pi_{1}+2\pi_{1}^{2}-\pi_{2})
S1(1+3π1+6π12+6π133π26π1π2+π3)-S_{1}(1+3\pi_{1}+6\pi_{1}^{2}+6\pi_{1}^{3}-3\pi_{2}-6\pi_{1}\pi_{2}+\pi_{3})

in agreement with the fact that 3D stellohedron has 1 hexagonal, 6 pentagonal and 3 quadrilateral faces, 24 edges and 16 vertices (see Fig. 4).

Corollary 7.2.

The exponential generating function of the numbers of vertices VnV_{n} of the stellohedra SnS_{n} has the form

(36) n=0Vn+1xnn!=ex1x.\sum_{n=0}^{\infty}V_{n+1}\frac{x^{n}}{n!}=\frac{e^{x}}{1-x}.

In particular, this implies that the number of vertices is

Vn=k=0n1(n1)!k!.V_{n}=\sum_{k=0}^{n-1}\frac{(n-1)!}{k!}.

More general formulas, expressing the number of faces of stellohedra of any dimension, can be found in [11] (Section 9) and [50] (Section 10.4).

For completeness consider also the families of simplices with nn vertices Δn\Delta_{n} and nn-dimensional cubes n\Box_{n} and define

Δ(x)=n1Δnxnn!,(x)=n0nxnn!.\Delta(x)=\sum_{n\geq 1}\Delta_{n}\frac{x^{n}}{n!},\quad\Box(x)=\sum_{n\geq 0}\Box_{n}\frac{x^{n}}{n!}.

The corresponding exponential generating functions

Δ(t,x):=n1FΔn(t)xnn!,(t,x):=n0Fn(t)xnn!\Delta(t,x):=\sum_{n\geq 1}F_{\Delta_{n}}(t)\frac{x^{n}}{n!},\quad\Box(t,x):=\sum_{n\geq 0}F_{\Box_{n}}(t)\frac{x^{n}}{n!}

satisfy the differential equations [15]

Δt(t,x)=xΔ(t,x),t(t,x)=2x(t,x).\Delta_{t}(t,x)=x\Delta(t,x),\quad\Box_{t}(t,x)=2x\Box(t,x).

Solving them explicitly with initial data Δ(0,x)=Δ(x),(0,x)=(x)\Delta(0,x)=\Delta(x),\,\Box(0,x)=\Box(x):

Δ(t,x)=Δ(x)etx,(t,x)=(x)e2tx,\Delta(t,x)=\Delta(x)e^{tx},\,\,\Box(t,x)=\Box(x)e^{2tx},

we come to the binomial formulas for the interiors

Δn=k=0n1(1)k(nk)Δnk,n=k=0n(2)k(nk)nk.\Delta^{\circ}_{n}=\sum_{k=0}^{n-1}(-1)^{k}{n\choose k}\Delta_{n-k},\quad\Box^{\circ}_{n}=\sum_{k=0}^{n}(-2)^{k}{n\choose k}\Box_{n-k}.

It is interesting to note that the polynomials pn(t)=FΔn(t)p_{n}(t)=F_{\Delta_{n}}(t) form the universal Appell sequence [6], satisfying the characteristic relation pn(t)=npn1(t)p_{n}^{\prime}(t)=np_{n-1}(t) (see [7], Ch. 19.3).

8. Concluding remarks

We have seen that associahedra and permutohedra play a special role in the differential algebra of graph-associahedra, generating subalgebras of the polytopal algebra 𝔓\mathfrak{P}, which are invariant under the differentiation. This explains their relation with the differential equations and ultimately with the inversion formulas.

The families of cyclohedra and stellohedra generate the subalgebras of 𝔓\mathfrak{P}, which can be viewed as the modules over the differential subalgebras generated by associahedra and permutohedra respectively. In the language of operads this was pointed out in [43] (see also [24, 58]). The description of all such subalgebras and modules over them within the differential algebra of polytopes 𝔓\mathfrak{P} seems to be a very interesting open problem.

The equations which appeared here are easily solvable, which hints the link with the theory of integrable systems. At the combinatorial level this link was already discussed in the literature, see e.g. [26, 39, 47, 54] and references therein. However, a general picture of the relations between combinatorics and integrable systems seems to be quite rich and needs better understanding (see [1, 2] for some very interesting thoughts in this direction).

In this relation it is worthy to mention that the Hopf equation ut=uuxu_{t}=uu_{x}, which appeared in our paper in relation with associahedra, is the dispersionless limit ε0\varepsilon\to 0 of the celebrated KdV equation

ut=uux+ε2uxxx,u_{t}=uu_{x}+\varepsilon^{2}u_{xxx},

well-known in soliton theory and enumerative algebraic geometry. In particular, according to Witten [60] the KdV hierarchy determines a certain generating function of the characteristic numbers of M¯g,n\bar{M}_{g,n} (see [40] for the detail), which makes the appearance of the inversion formula in the theory of the moduli spaces M¯0,n\bar{M}_{0,n} (discussed in Section 4) a bit less mysterious.

We would like to comment on the choice of the coefficients λn\lambda_{n} in the generating functions FP(x)=λnPnxn.F_{P}(x)=\sum\lambda_{n}P_{n}x^{n}. Most common cases are λn=1\lambda_{n}=1 and λn=1/n!\lambda_{n}=1/n!, corresponding to the usual and exponential generating functions. However, we have seen that in the case of cyclohedra the most appropriate choice is λn=1/n\lambda_{n}=1/n, related to the usual choice by integration. The role of different choices of λn\lambda_{n} in Boas-Buck’s approach [8] to the generating functions was emphasized in [12]. Note that both composition and multiplication of power series play the central role in this approach, going back to Appell [6] and providing umbrella for many classical polynomials.

Since the graph-associahedra are indecomposable, it is natural to ask if the corresponding hh-polynomials are irreducible over \mathbb{Q}. For example, for the permutohedra the corresponding hh-polynomials are the Eulerian polynomials An(t)A_{n}(t). It is known that A2k(t)A_{2k}(t) is divisible by t+1t+1 (which is hh-polynomial of the segment), but there is a conjecture that A2k+1(t)A_{2k+1}(t) and A2k(t)/(t+1)A_{2k}(t)/(t+1) are irreducible over \mathbb{Q} (see [34]). It would be interesting to study similar question for other graph-assiciahedra from our paper.

9. Acknowledgements.

We are very grateful to Vsevolod Adler, Alexander Braverman, Nikolai Erokhovets, Alexander Gaifullin and Georgy Shabat for the useful discussions, to Jose Bastidas, who attracted our attention to the preprint [4] and to Jim Stasheff for the encouraging comments.

References

  • [1] V.E. Adler On the combinatorics of several integrable hierarchies. J. Phys. A 48 (26) (2015), 265203.
  • [2] V.E. Adler Set partitions and integrable hierarchies. Theor. Math.Physics 187:3 (2016), 842-870.
  • [3] N. Amburg, E. Kreines, G. Shabat Poincare polynomial for the moduli space M0,n()¯\overline{M_{0,n}(\mathbb{C})} and the number of points in M0,n(𝔽q)¯.\overline{M_{0,n}(\mathbb{F}_{q})}. Moscow Univ. Math. Bull., 72 (2017),154-160.
  • [4] M. Aguiar, J. Bastidas Associahedra, cyclohedra and inversion of power series. arXiv:2010.14283 (2020).
  • [5] M. Aguiar, F. Ardila Hopf Monoids and Generalized Permutahedra. Memoirs of the AMS. Vol. 289, N.1437. AMS, 2023. arXiv:1709.07504 (2017).
  • [6] P. Appell Sur une classe de polynômes. Annales Sci. École Norm. Supér. 2e Sèrie. 9 (1880), 119-144.
  • [7] H. Bateman and A. Erdélyi Higher Transcendental Functions. Vol.III. McGraw-Hill Book Company, 1955.
  • [8] R.P. Boas, R.C. Buck Polynomial expansion of analytic functions. Springer, 1958.
  • [9] R. Bott and C. Taubes On the self-linking of knots. Journal of Math. Physics 35(10) (1994), 5247-5287.
  • [10] V. M. Buchstaber, E. V. Koritskaya The quasi-linear Burgers-Hopf equation and the Stasheff polytopes. Funct. Anal. Appl. 41:3 (2007), 196-207.
  • [11] V.M. Buchstaber Ring of simple polyhedra and differential equations. Trudy MIAN 263 (2008), 18-43.
  • [12] V.M. Buchstaber, A.N. Kholodov Boas-Buck structures on sequences of polynomials. Funct. Anal. Appl. 23:4 (1989), 266-276.
  • [13] V.M. Buchstaber, N. Erokhovets. Polytopes, Fibonacci numbers, Hopf algebras, and quasi-symmetric functions. Russian Math. Surveys 66 (2011), no. 2, 271-367.
  • [14] V. M. Buchstaber, V. D. Volodin Combinatorial 2-truncated cubes and applications. Progress in Mathematics, 299 (2012), Birkhäuser, Basel, 161-186.
  • [15] V.M. Buchstaber, T.E. Panov Toric Topology. Mathematical Surveys and Monographs, 204, AMS, 2015.
  • [16] V.M. Buchstaber, A.P. Veselov Chern-Dold character in complex cobordisms and theta divisors. Advances in Math. 449 (2024), 109720, 1-35.
  • [17] V.M. Buchstaber, A.P. Veselov Theta divisors and permutohedra. arXiv:2211.16042.
  • [18] V.M. Buchstaber, S. Tersic̀ Moduli space of weighted pointed stable curves and toric topology of Grassmann manifolds. arXiv:2410.01059.
  • [19] L. Carlitz Eulerian numbers and polynomials. Mathematics Magazine 32:5 (1959), 247-260.
  • [20] M.P. Carr and S.L. Devadoss Coxeter complexes and graph-associahedra. Topology Appl. 153 (2006), no. 12, 2155-2168.
  • [21] L. Comtet Advanced Combinatorics. D. Reidel Publ. Comp., Dordrecht, 1974.
  • [22] D.A. Cox, J.B. Little, and H.K. Schenck Toric Varieties. Graduate Studies in Mathematics, vol. 124, American Mathematical Society, Providence, RI, 2011.
  • [23] S.L. Devadoss Tesselations of moduli spaces and the mosaic operad. Contemp. Math. 239 (1999), 91–114.
  • [24] S.L. Devadoss Space of cyclohedra. Discrete and Comput. Geometry 29 (2002), 61-75.
  • [25] S.L. Devadoss and R.C. Read Cellular structures determined by polygons and trees. Annals of Combinatorics, 5 (2001), 71–98.
  • [26] G. Falqui, C. Reina and A. Zampa Krichever maps, Faà di Bruno polynomials, and cohomology in KP theory. Lett. Math. Phys. 42 (1997), 349-361.
  • [27] E.M. Feichtner and B. Sturmfels Matroid polytopes, nested sets and Bergman fans. Port. Math. (N.S.) 62 (2005), no. 4, 437-468.
  • [28] A. Frabetti, D. Manchon Five interpretations of Faà di Bruno’s formula. arXiv:1402.5551 (2014).
  • [29] A.A. Gaifullin Local formulae for combinatorial Pontrjagin classes. Izvestiya: Mathematics 68:5 (2004), 861-910.
  • [30] I. M. Gelfand, M. M. Kapranov, and A. V. Zelevinsky Discriminants, Resultants, and Multidimensional Determinants. Birkhäuser Boston Inc., Boston, MA, 1994.
  • [31] I.M. Gessel Lagrange inversion. J. Comb. Theory, Series A 144 (2016), 212-249.
  • [32] E. Getzler Operads and moduli spaces of genus 0 Riemann surfaces. Progress in Math. 129 (1995), Birkhäuser Boston, Boston, MA,199-230.
  • [33] S. M. Gusein-Zade, I. Luengo and A. Melle-Hernández A power structure over the Grothendieck ring of varieties. Math. Res. Lett. 11 (2004), 49-57.
  • [34] A.J.J. Heidrich On the Factorization of Eulerian Polynomials. J. Number Theory 18 (1984), 157-168.
  • [35] W.P. Johnson The curious history of Faà di Bruno’s formula. Amer. Math. Monthly 109:3 (2002), 217-234.
  • [36] M.M. Kapranov The permutoassociahedron, Mac Lane’s coherence theorem and asymptotic zones for the KZ equation. J. Pure Appl. Algebra 85 (1993), no. 2, 119–142.
  • [37] M.E. Kazaryan, S.K. Lando and V.V. Prasolov Algebraic Curves: Towards Moduli Spaces. Springer, 2019.
  • [38] M. Kontsevich Feynman diagrams and low-dimensional topology. Progr. Math. 120 (1994), 97-121.
  • [39] F. Lambert, J. Springael Soliton equations and simple combinatorics. Acta Appl. Math. 102 (2008), 147-178.
  • [40] S.K. Lando, A.K. Zvonkin Graphs on Surfaces and Their Applications. Encyclopaedia of Mathematical Sciences, 141, Berlin, New York: Springer-Verlag, 2004.
  • [41] J.-L. Loday The multiple faces of the associahedron. Clay Mathematics Institute Publication, 2005.
  • [42] A. Losev and Y. Manin New moduli spaces of pointed curves and pencils of flat connec- tions. Michigan Math. J. 48 (2000), 443-472.
  • [43] M. Markl Simplex, associahedron, and cyclohedron. Contemp. Math. 227, Amer. Math. Soc., Providence, RI, 1999.
  • [44] C.T. McMullen Moduli spaces in genus zero and inversion of power series. L’Enseign. Math. (2) 60 (2014), 25–30.
  • [45] P. McMullen Valuations and Euler-type relations on certain classes of convex polytopes. Proc. London Math. Soc. 35 (1977), 113-130.
  • [46] P. McMullen The polytope algebra. Adv. Math. 78 (1989), 76-130.
  • [47] Y. Mvondo-She From Hurwitz numbers to Feynman diagrams: Counting rooted trees in log gravity. Nuclear Physics B 995 (2023), 116350.
  • [48] J. Nicaise Grothendieck rings of polytopes and non-archimedean semi-algebraic sets. arXiv:2402.11699.
  • [49] V. Pilaud, F. Santos, G.M. Ziegler Celebrating Loday’s associahedron. Arch. Math. 121 (2023), 559–601.
  • [50] A. Postnikov, V. Reiner, L. Williams Faces of generalized permutohedra. Doc. Math. 13 (2008), 207-273.
  • [51] A. Postnikov Permutohedra, associahedra, and beyond. IMRN, 2009 (6), 1026-1106.
  • [52] A.V. Pukhlikov, A.G. Khovanski Finitely additive measures of virtual polytopes. St. Petersburg Math. J. 4 (1993), no. 2, 337-356.
  • [53] G.N. Raney Functional composition patterns and power series reversion. Trans. Amer. Math. Soc. 94 (1960), 441-451.
  • [54] A.B. Shabat, M.Kh. Efendiev On applications of Faà-di-Bruno formula. Ufa Math. Journal 9:3 (2017), 131-136.
  • [55] R. Simion A type-B associahedron. Advances in Applied Math. 30 (2003), 2-25.
  • [56] R.P. Stanley Enumerative Combinatorics. Vol. 2. Cambridge University Press, 1999.
  • [57] J. Stasheff Homotopy associativity of HH-spaces. I, II. Trans. Amer. Math. Soc. 108 (1963), 275–312.
  • [58] J. D. Stasheff The pre-history of operads. Contemp. Math. 202 (1997), 53-81.
  • [59] V. Toledano Laredo Quasi-Coxeter algebras, Dynkin diagram cohomology, and quantum Weyl groups. Int. Math. Res. Pap. IMRP 2008, Article ID rpn009, 167 pp.
  • [60] E. Witten Two-dimensional gravity and intersection theory on moduli space. Surveys Diff.Geom. 1 (1991), 243-310.
  • [61] G.M. Ziegler Lectures on Polytopes. Springer-Verlag, New York, 1995.