  The harmonic polytope and the bipermutahedron are two related polytopes which
arose in Ardila, Denham, and Huh's work on the Lagrangian geometry of matroids.
We study the bipermutahedron. We show that its faces are in bijection with the
vertex-labeled and edge-labeled multigraphs with no isolated vertices; the
generating function for its f-vector is a simple evaluation of the three
variable Rogers--Ramanujan function.
  We show that the h-polynomial of the bipermutahedral fan is the biEulerian
polynomial, which counts bipermutations according to their number of descents.
We construct a unimodular triangulation of the product of n triangles that is
combinatorially equivalent to (the triple cone over) the nth bipermutahedral
fan. Ehrhart theory then gives us a formula for the biEulerian polynomial,
which we use to show that this polynomial is real-rooted and that the h-vector
of the bipermutahedral fan is log-concave and unimodal.
  We describe all the deformations of the bipermutahedron; that is, the ample
cone of the bipermutahedral toric variety. We prove that among all polytopes in
this family, the bipermutahedron has the largest possible symmetry group.
Finally, we show that the Minkowski quotient of the bipermutahedron and the
harmonic polytope equals 2.
