The Schl\"afli Fan
Michael Joswig, Marta Panizzut, and Bernd Sturmfels

TL;DR
The paper introduces the Schl"afli fan, a detailed combinatorial structure that refines the classification of tropical cubic surfaces, revealing all line patterns and serving as a foundational tool for tropical geometry data analysis.
Contribution
It develops the theory of the Schl"afli fan and provides a framework for analyzing large-scale data in tropical geometry.
Findings
Refined classification of tropical cubic surfaces via the Schl"afli fan
Identification of all possible line patterns on tropical cubic surfaces
A blueprint for big data analysis in tropical geometry
Abstract
Smooth tropical cubic surfaces are parametrized by maximal cones in the unimodular secondary fan of the triple tetrahedron. There are such cones, organized into a database of symmetry classes. The Schl\"afli fan gives a further refinement of these cones. It reveals all possible patterns of lines on tropical cubic surfaces, thus serving as a combinatorial base space for the universal Fano variety. This article develops the relevant theory and offers a blueprint for the analysis of big data in tropical geometry.
| Marked Lines | Associated Motifs | Necessary Conditions |
| Isolated Lines | ||
| Exits: and . | ||
|---|---|---|
| Exits: , , and . | ||
| Exits: , . | ||
| Exits: , . | ||
| Exits: . | ||
| Exits: . | ||
| Exits: has exits also in and , . | ||
| Exits: has exits also in and , , . | ||
| Families of Lines | ||
| Exits: has exits also in and . | ||
| Exits: , | ||
| 3 | 4 | 6 | 8 | 12 | 24 |
| 3 | 15 | 25 | 82 | 10124 | 14363396 |
| Identifier | Canonical Hash | Altshuler Determinant |
| Index | Points | Exits |
| Motifs 3B | ||
| 0 | 9, 15, 7, 1, 18, 19 | 0, 1, 2, 3 |
| Motifs 3D | ||
| 1 | 9, 15, 2, 11, 1, 9, 15 | 1, 0, 2, 3 |
| 2 | 3, 14, 2, 11, 1, 15, 18 | 1, 0, 2, 3 |
| 3 | 9, 15, 2, 11, 1, 15, 18 | 1, 0, 2, 3 |
| 4 | 14, 15, 2, 11, 1, 15, 18 | 1, 0, 2, 3 |
| 5 | 3, 14, 2, 11, 1, 18, 19 | 1, 0, 2, 3 |
| 6 | 9, 15, 2, 11, 1, 18, 19 | 1, 0, 2, 3 |
| 7 | 14, 15, 2, 11, 1, 18, 19 | 1, 0, 2, 3 |
| 8 | 9, 15, 1, 11, 2, 3, 14 | 1, 3, 2, 0 |
| 9 | 9, 15, 1, 11, 2, 14, 15 | 1, 3, 2, 0 |
| 10 | 9, 15, 1, 11, 2, 9, 15 | 1, 3, 2, 0 |
| 11 | 2, 3, 14, 11, 17, 15, 18 | 0, 1, 2, 3 |
| 12 | 2, 3, 14, 11, 17, 18, 19 | 0, 1, 2, 3 |
| Motifs 3F | ||
| 13 | 15, 18, 11, 17, 14, 15, 2, 9 | 2, 3, 0, 1 |
| 14 | 18, 19, 11, 17, 14, 15, 2, 9 | 2, 3, 0, 1 |
| Motifs 3G | ||
| 15 | 9, 15, 2, 11, 3, 14 | 1, 3, 2, 0 |
| 16 | 9, 15, 2, 11, 14, 15 | 1, 3, 2, 0 |
| 17 | 9, 15, 1, 11, 15, 18 | 0, 1, 2, 3 |
| 18 | 9, 15, 1, 11, 18, 19 | 0, 1, 2, 3 |
| Motifs 3H | ||
| 19 | 7, 15, 1, 18, 19 | 0, 1, 2, 3 |
| 20 | 9, 15, 2, 14, 3 | 1, 3, 2, 0 |
| Motifs 3I | ||
| 21 | 1, 11, 9, 15 | 1, 2, 0, 3 |
| 22 | 2, 11, 9, 15 | 1, 2, 0, 3 |
| Motifs 3J | ||
| 23 | 11, 9, 15, 1, 2 | 0, 3, 1, 2 |
| Index | Points | Exits | Schläfli walls | ||
| Motifs 3A | |||||
| 0 | 18, 17, 15, 11, 2, 9 | 0, 2, 3, 1 |
|
||
| 1 | 18, 19, 15, 11, 2, 9 | 3, 2, 0, 1 | |||
| 2 | 18, 19, 15, 11, 2, 9 | 0, 2, 3, 1 | |||
| 3 | 18, 17, 15, 11, 1, 9 | 0, 2, 3, 1 |
|
||
| 4 | 18, 19, 15, 11, 1, 9 | 3, 2, 0, 1 | |||
| 5 | 18, 19, 15, 11, 1, 9 | 0, 2, 3, 1 | |||
| Motifs 3B | |||||
| 6 | 17, 18, 11, 1, 15, 7 | 0, 2, 1, 3 |
|
||
| 7 | 17, 18, 11, 1, 15, 9 | 0, 2, 1, 3 |
|
||
| 8 | 19, 18, 11, 1, 15, 7 | 0, 2, 1, 3 |
|
||
| 9 | 19, 18, 11, 1, 15, 9 | 0, 2, 1, 3 |
|
||
| Motifs 3D | |||||
| 10 | 1, 9, 15, 11, 17, 15, 18 | 0, 1, 2, 3 | |||
| 11 | 2, 9, 15, 11, 17, 15, 18 | 0, 1, 2, 3 | |||
| 12 | 1, 9, 15, 11, 17, 18, 19 | 0, 1, 2, 3 | |||
| 13 | 2, 9, 15, 11, 17, 18, 19 | 0, 1, 2, 3 | |||
| Motifs 3H | |||||
| 14 | 18, 19, 1, 13, 4 | 0, 2, 1, 3 | |||
| 15 | 11, 18, 1, 15, 9 | 0, 2, 1, 3 | |||
| 16 | 18, 19, 1, 13, 7 | 0, 2, 1, 3 | |||
| 17 | 11, 18, 1, 15, 7 | 0, 2, 1, 3 | |||
| Index | Points | Exits | Schläfli walls | |||||||||
| Motifs 3D | ||||||||||||
| 0 | 3, 14, 2, 11, 1, 9, 15 | 1, 0, 2, 3 |
|
|||||||||
| 1 | 14, 15, 2, 11, 1, 9, 15 | 1, 0, 2, 3 |
|
|||||||||
| 2 | 15, 18, 1, 11, 2, 3, 14 | 1, 3, 2, 0 |
|
|||||||||
| 3 | 18, 19, 1, 11, 2, 3, 14 | 1, 3, 2, 0 |
|
|||||||||
| 4 | 15, 18, 1, 11, 2, 9, 15 | 1, 3, 2, 0 |
|
|||||||||
| 5 | 18, 19, 1, 11, 2, 9, 15 | 1, 3, 2, 0 |
|
|||||||||
| 6 | 15, 18, 1, 11, 2, 14, 15 | 1, 3, 2, 0 |
|
|||||||||
| 7 | 18, 19, 1, 11, 2, 14, 15 | 1, 3, 2, 0 |
|
|||||||||
| Motifs 3H | ||||||||||||
| 8 | 9, 11, 1, 15, 7 | 0, 2, 1, 3 |
|
|||||||||
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
The Schläfli Fan
Michael Joswig
TU Berlin
Germany
,
Marta Panizzut
TU Berlin
Germany
and
Bernd Sturmfels
MPI Leipzig
Germany and UC Berkeley
USA
Abstract.
Smooth tropical cubic surfaces are parametrized by maximal cones in the unimodular secondary fan of the triple tetrahedron. There are such cones, organized into a database of symmetry classes. The Schläfli fan gives a further refinement of these cones. It reveals all possible patterns of lines on tropical cubic surfaces, thus serving as a combinatorial base space for the universal Fano variety. This article develops the relevant theory and offers a blueprint for the analysis of big data in tropical geometry.
1. Introduction
A cubic surface in projective -space is the zero set of a cubic polynomial
[TABLE]
Here are homogeneous coordinates on . George Salmon and Arthur Cayley discovered in the 1840s that every smooth cubic surface contains lines. Ludwig Schläfli studied the combinatorics of the lines in his 1858 article [23]. The name of that Swiss mathematician appears in our title.
This article is dedicated to the memory of Branko Grünbaum. Grünbaum is famous for his work on polytopes and arrangements, especially those that admit a high degree of symmetry. In the literature on these geometric figures, one sees a direct line connecting Ludwig Schläfli to Branko Grünbaum. This is highlighted by the use of the Schläfli symbol for symmetries of polyhedra.
The combinatorial strand of algebraic geometry underwent a major shift during the past two decades, thanks to the advent of tropical geometry [18]. The following question emerged early on during the tropical revolution: What are all shapes of smooth cubic surfaces in tropical -space, and which arrangements of tropical lines occur on such surfaces? A first guess is that there are lines, just like in the classical case. But this is false. Vigeland [24] showed that the number of lines can be infinite. A textbook reference is [18, Theorem 4.5.8].
The aim of this article is to give a comprehensive answer to the questions above. We will do so via a computational study of all smooth tropical cubic surfaces. These surfaces are dual to unimodular regular triangulations of the triple tetrahedron , which is the Newton polytope of the cubic polynomial seen in (1). The relevant definitions will be reviewed in Section 2.
Our point of departure is the article [22], which classifies the ten motifs that describe the potential positions of a tropical line on a cubic surface. These motifs are denoted 3A, 3B, , 3J. They are shown in Table 1. The advance we report in this paper is a large-scale computation that identifies the motifs of all lines that actually occur on the many tropical smooth cubic surfaces.
Our contribution rests on earlier work by Jordan et al. [16] who developed highly efficient tools for enumerating triangulations. Their count for in [16, Theorem 19] shows that there are combinatorial types of smooth tropical cubic surfaces. Here, the types are the orbits of the symmetric group by permuting in the terms of (1). Adding up the sizes of all -orbits, we obtain the total number of smooth tropical cubics.
This article is organized as follows. In Section 2 we fix notation, we discuss unimodular triangulations of the tetrahedron , and we review basics on lines and surfaces in tropical projective space . We also recall the classification of motifs in [22]. Section 3 furnishes our classification of smooth tropical cubic surfaces. This is presented in Theorem 3.1, and it is followed by a detailed explanation of the methodology that underlies our work and its results.
Section 4 studies occurrences of motifs in the unimodular triangulations of . Our main result is Theorem 4.1. We present an algorithm for computing occurrences. This rests on several lemmas that describe geometric constraints. The algorithm is applied to all triangulations in Theorem 3.1. As a consequence, we get a complete list of occurrences of motifs for each of the types.
In Section 5 we zoom in on particular secondary cones. For each cubic surface of one type, an occurrence of a motif may be visible or not. Being visible means that there exists a line for that motif. Hence, for any specific surface, only a subset of the motifs occurring in the triangulation is visible. The regions on which that subset is constant are convex polyhedral cones. These form the Schläfli fan. Thus, each of the secondary cones is divided into its Schläfli cones. We present and discuss the result of that computation.
Our combinatorial and computational study in this paper lays the foundation for future work on the nonarchimedean geometry of classical cubic surfaces over a valued field. In Section 6 we take a step into that direction. We discuss the universal Fano variety and the universal Brill variety, and we examine the tropical discriminants of these universal families. The first version of this article had a Section 7 which proposed a normal form for cubic surfaces, called the eight-point model. This was deleted in this final version because an even better such model was found in the subsequent project [21] with Emre Sertöz.
The methods from computer algebra and polyhedral geometry which led to our results are at the forefront of what is currently possible in terms of hardware, algorithms and software. For instance, to determine and analyze the regular unimodular triangulations of took more than 200 CPU days on an Intel Xeon E5-2630 v2 cluster. Yet the most difficult question we had to answer was how to make the results of such a large computation available to others. For this we set up a polymake extension TropicalCubics [17] and a database within the polyDB framework [20]. They can be accessed via polymake [8]. The database can also be used via an independent API. We believe that this approach can serve as a model for sharing “big data” in mathematical research.
2. Triangulations, Cubic Surfaces and Tropical Lines
In this section we review the basics and known results on which our study rests. For conventions on tropical geometry we follow the textbook by Maclagan and Sturmfels [18]. Our tropical semiring is the min-plus algebra . We use upper case letters to denote tropical variables and coefficients. Our orderings of variables and monomials are consistent with the conventions used by polymake [8]. For instance, here is a homogeneous tropical cubic polynomial:
[TABLE]
The expression (2) is evaluated in classical arithmetic as follows:
[TABLE]
The surface defined by (2) is the set of all points for which this minimum is attained at least twice. That tropical cubic surface lives in the tropical projective torus , but it also has a natural compactification in the tropical projective space . The latter is described in [18, Chapter 6].
A standard reference for the material that follows next is the textbook by De Loera, Rambau and Santos [4]. Reading the coefficients of the tropical polynomial as a height function defines a regular polyhedral subdivision of the lattice points in . If the coefficients are generic enough then the dual subdivision is a triangulation. For now the latter property may be taken as a definition for generic; it is a main point of later sections to refine this. If each of its tetrahedra has unit normalized volume, then the triangulation is unimodular and the tropical cubic surface is smooth. Every unimodular triangulation of the configuration has the same f-vector . Its boundary has the f-vector . From this we conclude that every smooth tropical cubic surface has vertices, edges, rays, bounded 2-cells, and unbounded 2-cells. This is the case in [18, Theorem 4.5.2]. Specifically, the interior edges of correspond to the bounded polygons in the surface. These polygons form the bounded complex of the tropical surface. This is also known as the tight span. For cubics, it is contractible. We define the B-vector of the triangulation to be , where denotes the number of -gons in the tight span. The GKZ-vector is , where is the number of tetrahedra containing point .
Example 2.1**.**
The tropical cubic polynomial in (2) is identified with its coefficient vector . This defines a unimodular triangulation of . Its tetrahedra are given by their labels:
[TABLE]
The GKZ-vector equals . The last entry means that the label occurs five times in (3). The B-vector of (3) is . To see this, we list the ten interior edges and their links:
[TABLE]
The link of an edge in is the graph of all edges in whose union with is a tetrahedron in . If is an interior edge of , then this graph is a cycle. For instance, the link of is the -cycle , , , . The corresponding bounded -cell in the tropical cubic surface is a quadrilateral. The triangulation (3) lies in the same -orbit as the one featured in [14, §6.2].
Each of the bounded edges of the surface determines a linear inequality among the coefficients , expressing that the edge has positive length. The secondary cone is the set of solutions to these inequalities. This is a full-dimensional cone in with -dimensional lineality space. The number of facets of is between and . The secondary cone of the triangulation (3) has facets. It contains the coefficient vector of (2).
The symmetric group acts naturally on the points in . This induces an action on the set of all triangulations. Note that also acts on the set of GKZ-vectors. The -orbit of the triangulation from (3) has size 24. Equivalently, the stabilizer of is trivial. The census of unimodular triangulations and associated cubic surfaces is presented in Theorem 3.1.
We now come to tropical lines in 3-space. Vigeland started the classification of how such lines can lie on generic smooth tropical cubic surfaces. Based on a massive random search with polymake, Simon Hampe realized that the classification was not complete. The triangulation (3) occurred in the joint article [14] as the first explicit counter-example to Vigeland’s list. The final characterization is due to Panizzut and Vigeland [22]. Their list of ten motifs is reproduced in Table 1. This table forms the foundation for our present study.
We identify with by setting , , and . A tropical line in is a balanced polyhedral complex given by two -valent adjacent vertices, joined by one bounded edge, and four rays with directions , , and . If the bounded edge has length zero, the tropical line is degenerate. Non-degenerate lines come in three labeled types, given by the direction of the bounded edge. This direction is either or or . We denote these three types by , and . This is shown in Figure 1.
Each tropical line in is encoded (up to tropical scaling) by its tropical Plücker vector . The six are the tropical minors of a -matrix. A vector is the tropical Plücker vector of a line if and only if it lies on the tropical hypersurface given by
[TABLE]
This means that the minimum in (4) is attained at least twice. Equivalently, is a height function on the six vertices of the regular octahedron which induces a split into two Egyptian pyramids [18, Figure 4.4.1]. The tropical hypersurface defined by (4) is the tropical Grassmannian Trop(G).
A tropical line is recovered from its Plücker vector as follows. We start by identifying the pair of terms in (4) which attains the minimum. Suppose , i.e., the labeled type is . Then, by [18, Example 4.3.19], consists of the segment joining the two points
[TABLE]
and the four rays , , , . The formulas for the other two labeled types, and , are analogous.
In summary, the vertices and of a tropical line are computed from the Plücker coordinates in (5). Conversely, the Plücker vector is obtained by taking the tropical minors of the -matrix with rows and .
The article [22] describes the various ways in which a tropical line can lie on a smooth cubic surface in -space. Here we require to be generic in the precise sense of Section 5. On the line we mark the points where intersects edges or vertices of the surface . These are the bars and dots indicated on the tropical lines in the left column of Table 1. Each bar is dual to a triangle in , and each dot is dual to a tetrahedron in . Formally, a motif of a tropical cubic surface is one of the ten abstract simplicial complexes 3A, 3B, , 3J which are listed in the middle column of Table 1. Each is equipped with a labeling of its vertices by and a marking of precisely four edges by . That this list of ten motifs is complete is the main result of [22].
The number of vertices of the ten motifs range between four and eight; the marked edges are the exits of the motif. The names of the motifs all start with the digit 3 to indicate the degree of the tropical surface; there are more motifs for other degrees [22, Table 2]. The article [22] distinguishes between “primal motifs” and “dual motifs”. We use the term motif for what is called “dual motif” in [22]. Our Table 1 uses for the homogeneous coordinates of the lattice points in , and it uses the notation for the facets of . The third column of Table 1 explicates additional conditions to be satisfied by some edges in order for the motif to occur in . These are derived in [22, Proposition 23]. They will become important in Section 4.
3. Data, Software, and Lines on Cubics
A primary goal of the present work is to present a database for smooth tropical cubic surfaces. We now explain our database and the underlying methodology. We start with the classification of combinatorial types. The proof of this result is the computation reported in [16, Theorem 19], plus an analysis of the orbits.
Theorem 3.1**.**
The triple tetrahedron has precisely regular unimodular triangulations. These are grouped into orbits with respect to the natural action of . The distribution of orbit sizes is shown in Table 2.
Remark 3.2**.**
Each smooth tropical cubic surface in has four elliptic curves in its boundary in . These are the tropical plane cubics which are dual to the induced triangulations of the ten lattice points in the triple triangle . That configuration has precisely unimodular triangulations, all of which are regular. They are grouped into orbits with respect to the natural action of . Hence, we encounter at most triangulations of the boundary . This means that, on the average, more than eight regular unimodular triangulations of induce the same boundary triangulation.
Before we enter the technical details, we briefly pause to reflect on the nature of a result like Theorem 3.1, how it can be useful, and to what extent it can be trusted. Theorem 3.1 is a highly condensed statement which was derived from massive computations, partially on large clusters, and the total time spent exceeds several months. Most readers will not have access to these types of hardware and technical resources and therefore will be unable to repeat these computations on their own. As we see it, the bulk of the data is the actual theorem. Theorem 3.1 is a mere corollary which follows from something which is too large to write down in any article. That data and more is made publically available at
[TABLE]
to allow everyone to derive their own corollaries. We stress that all the software that was used in the process is open source. Therefore, the entire proof of Theorem 3.1, which consists of software and data (in addition to this text), is available for scrutiny. Ideally, such a computer proof would be formalized, but currently this seems to be out of scope for a project of this size. Turning this into a formal proof would be a large project on its own, probably much larger than flyspeck [7], if feasible at all. This leaves the question of correctness.
As we see it, making data available and documenting this in an article is a necessary first step. Everyone is invited to probe the data for its correctness; we prepared various tools, explained below, to help with the probing. Any errors found in the future will be corrected in the data base. It would be desirable to have a general mechanism for this, accepted by the mathematical community. Finally, we would like to point out that it was a massive polymake experiment run by Simon Hampe which lead to the triangulation (3), which exhibited a flaw in a first version of [22]. That may be seen as a predecessor to this project.
High-level view on the data computed. For each of the triangulations in our database, the following annotations are reported: the GKZ-vector, the B-vector, the orbit size with respect to the -action, and a unique identifier. The identifier is an integer between and , which can be used to retrieve the triangulation and data derived. Frequently we will use the symbol ‘#’ for marking identifiers. The triangulation (3) has the identifier .
The facets of each triangulation are listed in lexicographic order. The representative for a combinatorial type is chosen such that the GKZ-vector is lexicographically minimal. Another important item in our database is a vector of minimum coordinate sum in the interior of each secondary cone. In order to find this vector, we had to solve an integer linear programming problem. We did this using the software SCIP [11]. The coefficients of the tropical polynomial (2) were derived from the triangulation (3) in this way. Note that, by construction, is always generic in the sense that the regular subdivision induced is a triangulation. However, it is not generic as defined in Section 5.
Exploring the database. We now describe how to access the data we produced. We offer a collection SchlaefliFan within the database Tropical of polyDB [20]. The simplest possible access is by directing a standard web browser to (6). However, for best results, we recommend the concurrent use of a recent version of polymake [8]. The new polymake extension TropicalCubics [17] is the software companion to this paper. It is available from and further explained at https://polymake.org/doku.php/extensions/tropicalcubics. Future additions will deal with other aspects of tropical cubic surfaces.
One pertinent question is how to find a given triangulation in the database. The user is unlikely to know the search key, and may be given by its list of facets as in (3). One way is to compute the GKZ-vector and to then generate the lexicographically minimal representative within its -orbit. This is the preferred method since it identifies the regular triangulation uniquely. Thus, in practice, the lex-minimal GKZ-vector works as another search key. An alternative method is to find a canonical form of as a simplicial complex. This means identifying the isomorphism type of the incidence graph of the vertices and the tetrahedra. The software nauty [19] is a standard tool for this task. It computes a canonical hash value, which is a 64-Bit integer that encodes the isomorphism type. This hash value is also stored in our database. It can be used as an index to retrieve a triangulation instantly; cf. Table 3.
The canonical hash value is a combinatorial invariant, but it is not unique. Table 3 shows two triangulations with the same hash value. Nonetheless, they are not isomorphic as abstract simplicial complexes, as can be seen as follows. Let and be an ordering of the vertices and the facets, respectively, of a simplicial complex . The incidence matrix is the -matrix with if vertex lies on the facet and otherwise. We define the Altshuler determinant of to be \max\bigl{(}\bigl{|}\det(J{J}^{\top})\bigr{|},\bigl{|}\det({J}^{\top}J)\bigr{|}\bigr{)}. This number does not depend on the orderings [1, Theorem 3]. It is a combinatorial invariant of . This distinguishes the third and fourth triangulations in Table 3. Our database can be queried for Altshuler determinants directly.
It also happens that two abstractly isomorphic triangulations lie in different -orbits. A pair of examples is given at the end of Table 3. Altogether there are hash values (i.e., about ) that correspond to two or more -orbits of triangulations. The maximal multiplicity of any hash value is four. So, with high probability, nauty identifies the triangulation uniquely.
Lines in surfaces. We now shift gears, with a discussion of the following basic problem. Given a non-degenerate tropical line and a tropical cubic surface , decide whether contains . We present an algorithm that solves this.
Let be an ordered list of linear polynomials . An interval in is covered by if the minimum value in the list is attained at least twice for all . This can only happen if some appear multiple times in . We introduce the coincidence partition
[TABLE]
where ( and ) implies ( if and only if ). We write for the linear function with . The tropical polynomial function defines a partition into smaller intervals,
[TABLE]
with the following property: on each precisely one function attains the minimum among our linear functions. Then covers if and only if
[TABLE]
Our discussion translates into an algorithm called the Covering Subroutine. Its input is an interval in and a list of linear polynomials, and its output is a yes-no decision whether is covered by . In the no-case, the Covering Subroutine also outputs a rational number such that the minimum in is attained only once. In the yes-case, the Covering Subroutine outputs the list of index sets , along with the corresponding tropical roots of . We call this list the covering certificate.
We next present an algorithm that decides whether a given non-degenerate tropical line lies on a given tropical cubic surface. It makes five calls to the Covering Subroutine. An illustration of Algorithm 1 is given in Example 3.3.
Example 3.3**.**
Fix the line with and the cubic with
[TABLE]
*This vector induces the honeycomb triangulation from [22, §6]: *
[TABLE]
The tropical line is non-degenerate and of labeled type because . Using (5) we find and . In all five iterations through steps 4–11, the answer is yes. The covering certificates are:
[TABLE]
There are two special points where is attained four times. At the point , the minimum is attained thrice. The relevant index sets are cells in the triangulation: two tetrahedra and , and the triangle . These data identify an occurrence of the motif 3D in Table 1.
Remark 3.4**.**
Algorithm 1 can be turned into a method for identifying all non-degenerate tropical lines in a given tropical surface in . Here is an alternative method for the same task. Let be the tropical polynomial defining . First we compute the dome . This is an unbounded polyhedron in which represents . We obtain a description of the surface as a polyhedral complex by projecting the codimension 2 skeleton of the dome. The maximal cells of are obtained by a convex hull computation [14, §3]. From this we enumerate the poset of all cells of ; cf. [15, Algorithm 1]. Each pair of cells is a candidate for possible locations of the two vertices and . These points are described as convex combinations of the cells’ vertices with unknown coefficients. Whether or not they form the two vertices of a tropical line in can be decided by checking the feasibility of a linear program.
Simon Hampe implemented a similar approach for tropical cubic surfaces. This is the function lines_in_cubic in the polymake extension a-tint [13], which is slightly different from our Algorithm 1. First, lines_in_cubic also computes degerate lines; second, that function is tailored to the cubic case.
4. Motifs and their Occurrences
We now turn to the ten motifs in Table 1. We are interested in their occurrences in the unimodular regular triangulations of . As before, our goal is the complete classification of all possibilities. We begin by stating our main result. The proof is given by exhaustive computations using Algorithm 2.
Theorem 4.1**.**
The number of occurrences of all motifs in the unimodular regular triangulations of varies between and , as shown in Figure 2. There are no triangulations with precisely , , or occurrences.
We now define the notion of occurrence. Fix a regular unimodular triangulation of . Let be a motif, viewed as a labeled simplicial complex. An occurrence of in is a simplicial map from to that satisfies the conditions in the third column of Table 1. These conditions include a bijection between the set of exits and the four facets of . Such a simplicial map sends vertices of to vertices of , while faces are mapped to faces. Often occurrences are embeddings, but it can happen that two vertices of are mapped to the same vertex of . We shall see this in Example 4.4.
An occurrence of a motif in is a map of simplicial complexes. The definition above is subtle. One might think that such a map is determined by the image of the set of vertices of . This is not true! The same subcomplex of may support several occurrences of a motif. We now present an example.
Example 4.2**.**
The line in Example 3.3 gives an occurrence of the motif 3D in the honeycomb triangulation. The corresponding simplicial map is given by
[TABLE]
*This uses our fixed ordering of the lattice points in , so the vertices are *
[TABLE]
The left diagram in Figure 3 helps in verifying the conditions from Table 1:
[TABLE]
The motif (12) is made visible in Example 3.3 by the line in the surface . The motif occurrence is seen in the covering certificates (11) given by Algorithm 1.
The above occurrence is special in that the exit edge lies in the edge of . We can relabel the points and the exits as follows:
[TABLE]
This is another occurrence of a 3D motif in , shown on the right in Figure 3.
In conclusion, the same subcomplex of the honeycomb triangulation supports two distinct occurrences of the motif 3D. However, it is impossible for both to be visible in the same cubic surface. To ascertain whether an occurrence of a motif is visible in a specific cubic surface is our problem in Section 5.
We now show all motif occurrences in a given triangulation. As it stands, Algorithm 2 is too naïve to be useful. The number of vertices of a motif varies between four (type 3I) and eight (type 3F). For the 3F motif alone we would need to enumerate and check potential simplicial maps into .
In practice, it is essential to exploit symmetries and other simplifications. A symmetry of a motif is a simplicial bijection from the labeled simplicial complex to itself such that the conditions in the third column of Table 1 are preserved. Two symmetric occurrences of a motif yield the same line in a given tropical surface (or none). The symmetries of a motif form a group. The following lemma is derived by direct inspection from the data in Table 1.
Lemma 4.3**.**
*The ten motifs of tropical cubic surfaces have the following symmetry groups. In each case, generators and a description are given: *
- (3A)
. Cyclic group of order .
- (3B)
, . Dihedral group of order .
- (3C)
, , . Elementary abelian group of order .
- (3D)
, . Elementary abelian group of order .
- (3E)
*, , . Nonabelian group of order : direct product of an order group and a dihedral group of order . *
- (3F)
*, , , , . Nonabelian group of order . Here span an abelian subgroup of order . *
- (3G)
, , . Elementary abelian group of order .
- (3H)
, . Elementary abelian group of order .
- (3I)
, . Elementary abelian group of order .
- (3J)
, . Elementary abelian group of order .
We next show that occurrences of motifs are generally not embeddings.
Example 4.4**.**
*The motif 3F occurs in the triangulation (3) via the labeling *
[TABLE]
In this occurrence, and are mapped to the same point, labeled by .
To obtain Theorem 4.1, we developed a highly efficient version of Algorithm 2, we implemented it in polymake, and we applied it to millions of triangulations. This required substantial speed-ups, based on structural constraints that control the combinatorial explosion. In the rest of this section, we present a sample of such constraints, and we discuss how they are used.
Lemma 4.5**.**
Vertex is distinct from and in any occurrence of motif 3A.
Proof.
If coincides with or then has coordinates , and equal to zero. Moreover, the condition implies that the coordinate is equal to one. This is impossible, since the four coordinates sum to three. ∎
Our strategy for enumerating motif occurrences is to find the possible ways in which the simplices of a motif are mapped into the given triangulation. This leads to more book-keeping in Algorithm 2, to be used for shortcuts. We exploit the following features in the various motifs. A tetrahedron is called sided if it has one edge on a facet of and the opposite edge lies on the plane . The associated tropical line contains the vertex of the surface dual to in the interior of the ray in direction . We call split if it has two opposite edges with prescribed exits. Here there are two possibilities. The line has two adjacent rays in directions given by the exits, and one ray contains in its interior the vertex dual to the split tetrahedron. Or the bounded edge contains the vertex dual to the split tetrahedron in its interior, and the rays in directions given by the exits are not adjacent. We say that is centered if the constraints in Table 1 induce a bijection between its vertices and the facets of . Its dual vertex lies in the interior of the bounded edge of the tropical line. Finally, a triangle in a motif is dangling if it has two edges with required exits. The tropical line has a vertex in the interior of an edge of the surface. The two rays adjacent to that vertex have direction given by the exits of the dangling triangle.
The features we defined above occur in the ten motifs as follows:
- •
The following tetrahedra are sided: in motif 3A, in 3C, and in 3E, and in 3F, in 3G, and in 3J.
- •
Tetrahedron in 3D is split; so are in 3F and in 3G.
- •
Tetrahedron in motif 3B is centered.
- •
Triangle in motif 3A is dangling, likewise and in 3B, in 3C, in 3D, in 3E, and in 3H.
Our strategy for Theorem 4.1 is to first enumerate the features of a triangulation, i.e. its sided, split and centered tetrahedra, and its dangling triangles. This is combined with searching for occurrences of a motif by local extensions.
We illustrate this for the 3A motif with a heuristic estimate for the number of subcases arising. Let be the triangulation in (3) and in Example 5.3 below. We start out by finding the candidates for the sided tetrahedron , with the exit on facet . Considering all labelings, there are choices for this in . Next we need to find the candidates for . Here it suffices to consider those which are in the link of the edge . For instance, the link for has six vertices. The number six appears to be typical and we use this number for our estimate. By Lemma 4.5, must be distinct from and , reducing the number of candidates to four. We further exclude any where does not lie in the boundary of . For the remaining ones we try the three directions other than , which is already fixed. The only item missing is the vertex . Assuming, e.g., we need to check three candidates in the link of (four minus one for , because ) and two remaining exits. This leads to cases, including all possible labelings. In fact, the enumeration is even faster, as many of these cases can be ruled out early while the various conditions in Table 1 are being checked. Summing up, the number of subcases considered by this approach is much smaller than 3.3 million subcases for one 3A motif one sees in a naïve backtracking search.
5. Schläfli Cones
In Section 4 we studied the occurrences of motifs in the types of regular unimodular triangulations of . Their number per type ranges between and . In this section we focus on individual smooth tropical cubic surfaces from a fixed secondary cone . Every tropical line on a generic surface gives a motif that occurs in . But the converse is not true. An occurrence of a motif need not contribute a tropical line to a given surface.
Let be a regular unimodular triangulation of . Each point in the open secondary cone specifies a smooth tropical cubic surface which is dual to the triangulation . Given an occurrence of a motif in , we say that is visible in if there is a tropical line in that has the dual complex . We write for the set of all motifs that are visible in .
We regard two vectors and in as equivalent if . Each equivalence class is a finite union of relatively open convex polyhedral cones in . The full-dimensional cones among these are the Schläfli cones. Each facet of a Schläfli cone is defined by a linear form in . This linear form is unique up to scaling. We identify this linear form with the hyperplane it defines, and we call it a Schläfli wall for the type . The collection of all Schläfli walls defines a hyperplane arrangement in .
The Schläfli fan of the combinatorial type is the subdivision of induced by the Schläfli walls of type . Every maximal cone of the Schläfli fan is fully contained in a Schläfli cone. Hence, the set is constant for all surfaces in a fixed maximal cone of the Schläfli fan. A tropical cubic surface is generic if its coefficient vector is in the interior of a Schläfli cone. Notice that, in general, a Schläfli cones need not be a cone of the Schläfli fan.
There are distinct Schläfli fans. Algorithm 3 finds their Schläfli walls. We coded this in Macaulay2 [12]. Here is one of the results we found:
Theorem 5.1**.**
For each of the types in Theorem 4.1 with exactly motifs, the secondary cone remains undivided in the Schläfli fan. Among these, types feature isolated tropical lines only. The remaining have precisely one occurrence of motif 3I; in particular, motif 3J does not occur at all.
The situation is different for many triangulations with more than motif occurrences. Here, the Schläfli fan is nontrivial; it does divide into smaller cones, according to which tropical lines lie on the various cubic surfaces. Each Schläfli wall arises (non-uniquely) from some motif that occurs in . If one crosses from one Schläfli cone to a neighboring one through a shared facet, then the set of visible motifs changes. If a motif is no longer visible then the Schläfli wall gives a linear inequality that is necessary for to be visible. We write for the set of Schläfli walls arising from the motif .
Lemma 5.2**.**
Let be an occurrence of a motif 3F, 3G or 3I in a type . Then . In other words, is visible in every tropical cubic surface of type .
Proof.
This was shown for the motifs 3G and 3I in [22, Proposition 23]. Now consider the motif 3F. Suppose that is an occurrence of 3F. The three tetrahedra , and are dual to three vertices of . The necessary conditions on the edges in Table 1 allow trespassing segments respectively in the directions , and . Thanks to the exits of the three tetrahedra, these segments can always be completed to a tropical line, irrespective of the specific values of the parameters . ∎
We now discuss how the set of walls can be computed for the other motifs. The basic idea is this. Given , we compute a tropical line that matches the combinatorics in . The line is uniquely determined by its two vertices. Their coordinates are linear forms in . In this section we use lowercase letters instead of uppercase letters for the coordinates of the tropical coefficient vector , so as to make our tables more readable.
If the take on values in then the tropical line may or may not be contained in . We require that it lies in as prescribed by . Each vertex must lie on a cell of that is specified by the tropical cubic polynomial. These linear forms must be equal and bounded above by the other ones. We consider these linear inequalities together with those that define the secondary cone. They define the visibility cone of in . The irredundant linear inequalities for this cone give us the linear forms in . A Schläfli cone is the intersection of the full dimensional visibility cones of the visible motifs.
If all linear inequalities we found are redundant, then the visibility cone equals the secondary cone. In that case, the motif is visible in each surface with , and the motif is globally visible. This holds in Theorem 5.1.
If Algorithm 3 finds irredundant linear forms, then we distinguish two cases, according to the dimension of the visibility cone. If the visibility cone is full dimensional, then the motif is partially visible. Finally, a visibility cone might not be full dimensional. This means that it is contained in a linear space of positive codimension. A motif with visibility cone of lower dimension is not visible in generic surfaces. We therefore call it hardly visible.
We now illustrate these concepts for the tropical cubic surface from (3).
Example 5.3**.**
The triangulation has occurrences of the motifs 3A, 3B, , 3J. Their frequencies are . Lemma 5.2 says that the motifs 3F, 3G and 3I are globally visible. In Tables 4, 5 and 6 we list all motifs, together with their sets of Schläfli walls . We describe how the Schläfli walls are computed for the motifs of type 3H. The motif consists of a tetrahedron and a dangling triangle . One of the vertices of the tropical line defined by is dual to the tetrahedron. In order for the line to be contained in the surface, the other vertex must lie on the edge dual to the dangling triangle, i.e., the minimum in the tropical polynomial must be achieved at the monomials corresponding to , and . These linear inequalities define the visibility cone. Note that the occurrence 8 of motif 3H in Table 6 is hardly visible, since its visibility cone is not full dimensional.
The list of partially visible motifs in Table 5 shows that the Schläfli walls generate a hyperplane arrangement defined by the seven linear forms:
[TABLE]
We write and for the two halfspaces defined by these linear forms.
Let us look at the Schläfli walls from partially visible motifs of type 3B. The hyperplanes for the Schläfli walls of these motifs are and . They divide the secondary cone into four cells , , , . These four cells correspond in Table 5 to the occurrences 8, 6, 7 and 9, in this order. Each motif occurrence is visible in precisely that cell.
For the motifs of type 3D, we also have two hyperplanes and . These give the Schläfli walls that divide the secondary cone into four cells. In the cells and the motifs 11, 13 and 10, 12 are visible, respectively. In the cell none of the partially visible motif is visible. Finally, on the cell all the partially visible motifs are visible. Moreover, when we pass through the Schläfli wall from to , the motifs 0 and 1 of type 3A are no longer visible, while the motifs 11 and 13 of type 3D become visible.
Remark 5.4**.**
In this section we have considered the problem whether a motif is visible on a certain tropical surface. If the motif is visible, the next natural step is to ask whether the tropical line defined by the motif is realizabile on a cubic surface defined over a field with valuation. More precisely, given a tropical line on a tropical surface , we ask whether there exits a line on a cubic surface such that and . This realizability problem has been studied in [2] and [3]. The authors show that non-degenerate lines in families of type 3I are not realizable on surfaces over a valued field of characteristic zero. Moreover, in the recent article [9] the result is extended to valued field with residue field of characteristic different from two. The paper also provides an example of a line of type 3J which is realizable on a cubic surface defined over the field of -adic numbers.
6. The Universal Fano Variety and its Tropical Discriminant
We now relate our combinatorial results to classical algebraic geometry. The natural parameter space for our problem is the universal Fano variety. Its points are pairs consisting of a line and a cubic surface that contains it. The map onto the second factor is a -to- cover of . The fiber over a smooth cubic surface, regarded as a point in , is the Fano variety on that surface, i.e. the set of its lines. The branch locus of the -to- map is its discriminant, a hypersurface in . We shall see that the codimension one skeleton of the Schläfli fan plays the role of the tropical discriminant for this map.
We follow the approach to tropical geometry in the textbook [18]. One starts with a classical variety, defined by an ideal in a (Laurent) polynomial ring over a field with valuation. The tropical variety is the set of all weight vectors whose initial ideal contains no monomials. We would like to apply this to the universal Fano variety for lines on cubic surfaces, represented by an ideal in the polynomial ring in the unknowns and . This is the homogeneous coordinate ring of . The first factor contains the Grassmannian of lines in as a quadratic hypersurface in .
The quadric defining is the Pfaffian of the skew-symmetric matrix
[TABLE]
We have . The line with Plücker coordinates is the image in of the column span of the associated rank matrix .
The second factor parametrizes cubic forms . Its coordinates are the coefficients . Fix a row vector of unknowns and form the vector-matrix product . We write for the polynomial obtained by replacing with . Thus, is a homogeneous cubic in . Its coefficients are bihomogeneous forms of degree , like
[TABLE]
We write for the ideal in that is generated by these polynomials together with the Plücker quadric .
The zero set of in is the universal Fano variety of lines on cubic surfaces. We verified by computations on affine charts that the ideal defines the correct scheme. We consider the tropical universal Fano variety
[TABLE]
By the Structure Theorem [18, Theorem 3.3.5], is a pure -dimensional balanced fan. For simplicity, we disregard boundary phenomena, and we replace each tropical projective space with its dense tropical torus . The former is compact while the latter is not. For a detailed discussion see [18, §6.2]. The points in are the pairs consisting of a line in and a cubic surface that contains the line. The tropical line is represented by its Plücker vector . The cubic is represented by its coefficient vector . Unlike in previous sections, this tropical cubic not be tropically smooth. A pair lies in if and only if contains no monomial. We take this initial ideal in the Laurent polynomial ring.
Example 6.1**.**
The line given by lies on the surface given by . This pair corresponds to the motif of type 3D in Example 3.3; see the diagram on the left-hand side of Figure 3. We verify the containment algebraically by checking that
[TABLE]
contains no monomial. This initial ideal lives in the Laurent polynomial ring. For instance, the ten terms in (13) have weights , in this order, and the resulting initial form equals .
The point lies in the relative interior of a maximal cell of . The inequality description of this cell is read off from a Gröbner basis of . For instance, the polynomial (13) contributes the equation and eight inequalities, namely, is bounded above by
[TABLE]
Such constraints, derived from polynomials in , define the cells of .
The maximal cones of represent occurrences of motifs in . In particular, if we could compute this fan, then this would be an ab initio derivation of the motifs 3A, 3B, , 3J. These were found geometrically in [22].
Remark 6.2**.**
Motifs and their occurrences can be identified from initial ideals . For instance, the indices of the unknowns in (14) form the list we saw in Example 4.2.
Unfortunately, it is very difficult to compute with the ideal . Even finding a single Gröbner basis is hard. For instance, the computation of (14) only terminated after we imposed some degree constraints in Macaulay2. One open problem naturally arising here is to find a tropical basis of . The Schläfli fan fits into a broader theory, yet to be developed, for discriminants of morphisms in tropical algebraic geometry. We propose the following approach. Let be a tropical variety in and the projection from onto the second factor . We assume that is onto, so . Let be the subcomplex of consisting of all cells of dimension at most . If this is the ramification locus then plays the role of the branch locus.
Example 6.3**.**
Tropical discriminants [5] are a special case of this construction. Let be a subset of elements in . Consider hypersurfaces in -space defined by Laurent polynomials with these terms. We write for the vector of coefficients, and for a point in . The universal tropical hypersurface is the tropical variety defined by
[TABLE]
The map is surjective. The fiber is the hypersurface in whose tropical polynomial has coefficients . The tropical variety has dimension . It is a fan with maximal cones, one for each pair of terms in (15). The subfan consists of cones of dimension . On each cone, the minimum among the terms in (15) is attained by a fixed set of terms. Hence the regular subdivision of defined by is not a triangulation. The image consists of the cones of codimension in the secondary fan of . In particular, the tropical discriminant defined above contains that of [5]. The difference arises from the distinction between the -discriminant and the principal -determinant; see [6] and [10].
Remark 6.4**.**
The number of cones in is much smaller than that of its image under the projection . This phenomenon is familiar from computer algebra (cf. elimination theory) and optimization (cf. extended formulations). In our context, take in Example 6.3. The universal cubic surface has only maximal cones, whereas its discriminant forms the walls between many more than cones.
We now come to the main theoretical result in this section. The role of points will be played by lines. An analogous result holds for Fano varieties of arbitrary hypersurfaces (15). We focus on the case of cubic surfaces in .
Proposition 6.5**.**
Let be the tropical universal Fano variety in and the map onto the second factor (space of tropical cubics). The tropical discriminant of is contained in the union of the codimension cones in the Schläfli fan. The latter is a subset of the union of all Schläfli walls.
Proof.
All cubics in the interior of one fixed Schläfli cone have the same visible motifs. The Plücker vectors of the lines are linear functions in the entries of , as long as staying within one Schläfli cone. Hence the set of cells in that are intersected by the fiber remains constant throughout that Schläfli cone. These cells all have the full dimension . In particular, is disjoint from for in the interior of a Schläfli cone. This shows that this interior is disjoint from the tropical discriminant of . ∎
We conclude with a brief discussion of a related universal family. It lives in , where now parametrizes planes in the ambient -space. Each plane intersects a cubic surface in a plane cubic curve. The plane is a tritangent plane if the plane cubic decomposes into three lines. The universal Brill variety is the -dimensional irreducible variety consisting of all pairs , where is a tritangent plane to the cubic surface . The map from this variety onto is a -to- covering, since a general cubic surface has tritangent planes.
We introduce an ideal that defines the universal Brill variety. It lives in the ring , where the last ten unknowns are the coefficients of a ternary cubic. In these unknowns, we consider the prime ideal of codimension and degree that defines the factorizable cubics. Its variety is an instance of a Chow variety, and the equations are known as Brill equations [10, §I.4.H]. This prime ideal is generated by quartics in the unknowns.
We now derive generators of . Set in , and clear denominators to get a ternary cubic with coefficients . We substitute these cubics into the Brill equations and we remove factors of . The resulting polynomials of bidegree in generate our ideal .
We are interested in the resulting tropical universal Brill variety
[TABLE]
Its points are pairs consisting of a tropical cubic and a tritangent plane. The maximal cones of represent occurrences of triple motifs in . It would be desirable to compute these. We note that the tritangent planes correspond to the triangles in the Schläfli graph. This is the -regular graph whose vertices are the lines, and whose edges are incident pairs of lines. The motifs and the triple motifs that occur in a triangulation can be seen as a tropical structure for annotating and extending the Schläfli graph.
Acknowledgements
We are very grateful to Lars Kastner, Benjamin Lorenz and Andreas Paffenholz for their help with the computations for this project. We thank Sara Lamboglia, Yue Ren and Emre Sertöz for their comments on a manuscript version of this article. We are also grateful to the two anonymous referees whose comments helped us improving the article. Michael Joswig was supported by Deutsche Forschungsgemeinschaft (EXC 2046: "MATH+", SFB-TRR 195: "Symbolic Tools in Mathematics and their Application", and GRK 2434: "Facets of Complexity").
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] Altshuler, A., Steinberg, L.: Neighborly 4 4 4 -polytopes with 9 9 9 vertices. J. Combinatorial Theory Ser. A 15 , 270–287 (1973)
- 2[2] Bogart, T., Katz, E.: Obstructions to lifting tropical curves in surfaces in 3-space. SIAM J. Discrete Math. 26 (3), 1050–1067 (2012)
- 3[3] Brugallé, E., Shaw, K.: Obstructions to approximating tropical curves in surfaces via intersection theory. Canad. J. Math. 67 (3), 527–572 (2015)
- 4[4] De Loera, J.A., Rambau, J., Santos, F.: Triangulations, Algorithms and Computation in Mathematics , vol. 25. Springer-Verlag, Berlin (2010)
- 5[5] Dickenstein, A., Feichtner, E.M., Sturmfels, B.: Tropical discriminants. J. Amer. Math. Soc. 20 (4), 1111–1133 (2007)
- 6[6] Dickenstein, A., Tabera, L.F.: Singular tropical hypersurfaces. Discrete Comput. Geom. 47 (2), 430–453 (2012)
- 7[7] https://github.com/flyspeck/flyspeck
- 8[8] Gawrilow, E., Joswig, M.: polymake : a framework for analyzing convex polytopes. In: Polytopes – combinatorics and computation (Oberwolfach, 1997), DMV Sem. , vol. 29, pp. 43–73. Birkhäuser (2000)
