Computing parametric weighted Ehrhart polynomials of smooth polytopes
Daniel Hwang, Juliet Whidden, Josephine Yu
TL;DR
This paper shows that when integral polytopes $P_A(\mathbf{b})$ deform with fixed facet normals, the coefficients of weighted Ehrhart and $h^*$-polynomials are piecewise polynomial in $\mathbf{b}$. It develops and implements an algorithm based on the Khovanskii–Pukhlikov weighted Euler–Maclaurin formula to compute these polynomials for smooth polytopes, including type $A$ alcoved polytopes, and provides a catalog of results. It also investigates the signs of weighted $h^*$-coefficients, proving that positivity regions are not universally convex and analyzing the asymptotic behavior under large dilates, where roots converge to those of Eulerian polynomials. The work offers practical parametric tools for computing weighted Ehrhart polynomials and deepens understanding of positivity and root structure in weighted Ehrhart theory, with implications for combinatorial and polyhedral geometry.
Abstract
We show that when integral polytopes are deformed while keeping the same facet normal vectors, the coefficients of weighted Ehrhart and $h^*$-polynomials are piecewise polynomial functions in the ``right hand sides'' of the linear inequalities defining the polytopes. We give an algorithm and an implementation in SageMath for computing these polynomials for smooth polytopes, such as type $A$ alcoved polytopes, using a weighted Euler-Maclaurin type formula by Khovanskiǐ and Pukhlikov. We discuss some natural questions concerning signs of the coefficients of the weighted $h^*$-polynomials.
