MULTIVARIATE VOLUME, EHRHART, AND h∗-POLYNOMIALS OF POLYTROPES MARIE-CHARLOTTE BRANDENBURG, SOPHIA ELIA, LEON ZHANG ∗ ABSTRACT. The univariate Ehrhart and h -polynomials of lattice polytopes have been widely stud- ied. We describe methods from toric geometry for computing multivariate versions of volume, Ehrhart and h∗-polynomials of lattice polytropes, which are both tropically and classically convex. These algorithms are applied to all polytropes of dimensions 2, 3 and 4, yielding a large class of integer polynomials. We give a complete combinatorial description of the coefficients of volume polynomials of 3-dimensional polytropes in terms of regular central subdivisions of the fundamental polytope. Finally, we provide a partial characterization of the analogous coefficients in dimension 4. 1. INTRODUCTION Tropical polytopes are analogues of ordinary polytopes arising naturally both in tropical geometry and in a range of algorithmic applications to other fields, including optimization [JS19], phyloge- netics [YZZ19], mechanism design [CT18], and building theory [JSY07]. Several recent articles have focused on the question of how to define and compute the volume of a tropical polytopes. One approach is simply to consider their discrete volume, which is relevant to a number of appli- cations: for instance, the irreducible components of a Mustafin variety correspond to lattice points of a certain computable tropical polytope [Car+11; Zha20]. Polytropes, which are both tropically and classically convex [JK08], form a fundamental class of tropical polytopes. They also appear in the literature as alcoved polytopes of type A. Polytropes are closely related to the notion of Kleene stars and the problem of finding shortest paths in weighted graphs [Tra17; JS19]. Both the volume and the nonnegativity of the h-vector of poly- tropes have been studied in the language of affine Coxeter arrangements [LP07]. A certain family of polytropes called Lipschitz polytopes corresponding to finite metric spaces has also been heavily investigated (e.g. [Ver15; C¸el+20]). In this article we examine the multivariate Ehrhart theory of polytropes, and describe methods for computing the multivariate volume, Ehrhart, and h∗-polynomials for all polytropes. We begin by arXiv:2006.01920v1 [math.CO] 2 Jun 2020 describing the Ehrhart theory, tropical geometry, and Grobner¨ theory necessary for our methods in Section 2. In Section 3 we describe methods for computing multivariate Ehrhart polynomials and volume polynomials of polytropes. In Section 4 we apply these methods to compute the volume, Ehrhart, and h∗-polynomials for all polytropes of dimension at most 4, and we study the combinatorics deriving from the coefficients of the volume polynomials. In recent work, Loho and Schymura [LS19] developed a notion of volume for tropical polytopes driven by a tropical version of dilation, which yields an Ehrhart theory for a new class of tropical lattices. This notion of volume is intrinsically tropical and exhibits many natural properties of a volume measure, such as being monotonic and rotation-invariant. However, their notion dis- cards the classical lattice Zd in favor of tropical lattices and thus is less relevant for applications involving ordinary lattices. 1 De la Puente and Claver´ıa proved the Mahler conjecture for centrally symmetric 3-dimensional polytropes by computing their volume polynomials in [PC18]. Their approach to volume polyno- mials differs from ours, relying instead on a description of 3-dimensional polytropes as perturbed 3-dimensional boxes. By contrast, our approach works in arbitrary dimension; it produces six dif- ferent volume polynomials in 12 variables to cover all combinatorial types of polytropes, whereas de la Puente and Claver´ıa use a single volume polynomial in 21 variables. Acknowledgements. The authors thank Christian Haase, Michael Joswig, Benjamin Schroter,¨ Rainer Sinn, and Bernd Sturmfels for many helpful conversations. They also thank Michael Joswig, Lars Kastner, and Benjamin Schroter¨ for providing a dataset for 4-dimensional polytropes, and the Max Planck Institute for Mathematics in the Sciences for its hospitality while working on this project. Leon Zhang was partially supported by an NSF Graduate Research Fellowship. Sophia Elia was supported by the Deutsche Forschungsgemeinschaft (DFG) Graduiertenkolleg “Facets of Complexity”. 2. BACKGROUND In this section we give a brief overview of the background material we use for our results. Note that throughout this article, we assume P is a lattice polytope unless stated otherwise. 2.1. Ehrhart theory. In this subsection, we recall some essentials of Ehrhart theory. For more details, we refer the reader to [BR15]. A lattice polytope P ⊆ Rn is the convex hull of a finite set of points in Zn. The Ehrhart counting function of P, written ehrP(k), gives the number of lattice points in the k-th dilate of P for k 2 Z≥1: n ehrP(k) = jkP \ Z j. Ehrhart’s theorem [Ehr62] says that ehrP(k) is given by a polynomial in k of degree equal to the dimension of P. Furthermore, the constant term of this polynomial is equal to 1 and the coefficient of the leading term is equal to the Euclidean volume of P within its affine span. The interpretation of other coefficients of the Ehrhart polynomial is an active direction of research. Generating functions play a central role in Ehrhart theory. The Ehrhart series of a polytope P, written EhrP(t), is the formal power series given by: k EhrP(t) = 1 + ∑ ehrP(k)t . k≥1 For a d-dimensional lattice polytope, the Ehrhart series has the following rational expression: h∗ (t) Ehr (t) = 1 + ehr (k)tk = P , P ∑ P ( − )d+1 k≥1 1 t ∗ d ∗ i ∗ where hP(t) = ∑i hi t is a polynomial in t of degree at most d. Furthermore, each hi is a non- negative integer [Sta80]. The normalized volume of P is defined as Vol(P) = dim(P)! vol(P), where vol(P) is the Euclidean volume of P within its affine span. It is equal to the sum of the coefficients of the h∗-polynomial. The Ehrhart polynomial may be recovered from the h∗-vector through the following transforma- tion: d ∗ k + d − i ehrP(k) = ∑ hi . i=0 d 2 Polytopes also admit an inequality description as bounded intersections of finitely many closed half-spaces: P = fx 2 Rn : Ax ≤ ag for some A 2 Rm×n and a 2 Rm. The univariate Ehrhart polynomial is a function in a single variable k, which counts the number of lattice points in kP = fx 2 Rn : Ax ≤ kag, in which all facet-defining hyperplanes of P are moved uniformly. A multivariate Ehrhart polynomial is polynomial in n variables a = (a1,..., an), which allows for the independent movement of facets. 2.2. Tropical convexity. In this subsection we review some basics of tropical arithmetic and trop- ical convexity. We refer readers to [DS04] or [Jos] for a more detailed exposition. Over the min-plus tropical semiring T = R [ f¥g, we define for a, b 2 T the operations of addition a ⊕ b and multiplication a b by a ⊕ b = min(a, b), a b = a + b. We can similarly define vector addition and scalar multiplication: for any scalars a, b 2 T and for n any vectors v = (v1,..., vn), w = (w1,..., wn) 2 T , we define a v = (a + v1, a + v2,..., a + vn) a v ⊕ b w = (min(a + v1, b + w1), . , min(a + vn, b + wn)). n Let V = fv1,..., vrg ⊆ R be a finite set of points. The tropical convex hull of V is given by the set of all tropical linear combinations tconv(V) = fa1 v1 ⊕ · · · ⊕ ar vrj a1,..., ar 2 Rg. A tropically convex set in Rn is closed under tropical scalar multiplication. As a consequence, we can identify a tropically convex set P contained in Rn with its image in the tropical projective torus TPn−1 = Rn/R(1, . , 1).A tropical polytope is the tropical convex hull of a finite set V in TPn−1. A tropical lattice polytope is a tropical polytope whose spanning points are all contained in Zn. Let P = tconv(V) ⊆ TPn−1 be a tropical polytope. The (tropical) type of a point x in TPn−1 with respect to V is the collection of sets S = (S1,..., Sn), where an index i is contained in Sj if (vi)j − xj = min((vi)1 − x1,..., (vi)n − xn). i 2 S2 i 2 S2 \ S3 vi i 2 S1 i 2 S3 2 FIGURE 1. The max-tropical hyperplane Hi ⊆ TP in the chart where the third coordinate is 0, with faces labeled for type identification. n−1 Geometrically, we can view the type of x as follows: a max-tropical hyperplane Ha ⊆ TP with n−1 n−1 apex at a 2 TP is the set of points y 2 TP such that the maximum of fai + yi : i 2 [n]g is attained at least twice. The max-tropical hyperplane H0 induces a complete polyhedral fan 3 n−1 F0 in TP . For each i 2 [r], let Hi be a max-tropical hyperplane with apex vi. Each of these n−1 hyperplanes determines a translate FHi of F0. Two points x, y 2 TP lie in the same face of FHi if and only if vi − x and vi − y achieve their minima in the same coordinates. For a point n−1 x 2 TP with type S = (S1,..., Sn), the set Sj records for which hyperplanes Hi the point x lies in a face such that vi − x is minimal in coordinate j. Figure 1 shows when i is contained in Sj based on the position of x in TP2. The tropical polytope P consists of all points x whose type S = (S1,..., Sn) has all Si nonempty.