arxiv: v2 [math.ac] 5 Jan 2018

Transcription

1 Random Monomial Ideals Jesús A. De Loera, Sonja Petrović, Lily Silverstein, Desina Stasi, Dane Wilburne arxiv:70.070v [math.ac] Jan 8 Abstract: Insired by the study of random grahs and simlicial comlexes, and motivated by the need to understand average behavior of ideals, we roose and study robabilistic models of random monomial ideals. We rove theorems about the robability distributions, exectations and thresholds for events involving monomial ideals with given Hilbert function, Krull dimension, first graded Betti numbers, and resent several exerimentally-backed conjectures about ularity, rojective dimension, strong genericity, and Cohen-Macaulayness of random monomial ideals.. Introduction Randomness has long been used to study olynomials. A commutative algebraist s interest in randomness stems from the desire to understand average or tyical behavior of ideals and rings. A natural aroach to the roblem is to define a robability distribution on a set of ideals or varieties, which, in turn, induces a distribution on algebraic invariants or roerties of interest. In such a formal setu, questions of exected (tyical) or unlikely (rare, non-generic) behavior can be stated formally using robability. Let us consider some examles. Already in the 90 s, Littlewood and Offord [] studied the exected number of real roots of a random algebraic equation defined by random coefficients. The investigations on random varieties, defined by random coefficients on a fixed Newton olytoe suort, have generated a lot of work now considered classical (see, e.g., [, 9, 9] and the references therein). A robabilistic analysis of algorithms called smooth analysis has been used in algebraic geometry, see [, 9]. Roughly seaking, smooth analysis measures the exected erformance of an algorithm under slight random erturbations of worst-case inuts. Another more recent examle of a notion of algebraic randomness aears in [, ], where they consider the distribution of Betti numbers generated randomly using concets from the Boij-Söderberg theory [8]. Still, there are many other examles of algebraic toics in which robabilistic analysis lays a useful role (see e.g., [] and the references therein). Our aer introduces a new robabilistic model on monomial ideals inside a olynomial ring. Why work on a robabilistic model for monomial ideals? There are at least three good reasons: First, monomial ideals are the simlest of ideals and lay a fundamental role in commutative algebra, most notably as Gröbner degenerations of general ideals, caturing the range of values shown in all ideals for several invariants (see [, 7]). Second, monomial ideals rovide a strong link to algebraic combinatorics (see [,, 8]). Third, monomial ideals naturally generalize grahs, hyergrahs, and simlicial comlexes; and, since the seminal aer of Erdös and Rényi [9], robabilistic methods have been successfully alied to study those objects (see, e.g., [, 7, 7] and the references therein). Thus, our work is an extension of both classical robabilistic combinatorics and the new trends in stochastic toology. Our goal is to rovide a formal robabilistic setu for constructing and understanding distributions of monomial ideals and the induced distributions on their invariants (degree, dimension, Hilbert This work is suorted by the NSF collaborative grants DMS- to Illinois Institute of Technology and DMS-8 to the Univ. of California, Davis.

2 De Loera, Petrović, Silverstein, Stasi, Wilburne/Random monomial ideals function, ularity, etc.). To this end, we define below a simle robability distribution on the set of monomial ideals. Drawing a random monomial ideal from this distribution allows us to study the average behavior of algebraic invariants. While we are interested in more general robability distributions on ideals, some of which we describe in Section, we begin with the most basic model one can consider to generate random monomial ideals; insired by classical work, we call this family the Erdős-Rényi-tye model, or the ER-tye model for random monomial ideals. The ER-tye model for random monomial ideals Let k be a field and let S = k[x,..., x n ] be the olynomial ring in n indeterminates. We construct a random monomial ideal in S by roducing a random set of generators as follows: Given an integer D and a arameter = (n, D), 0, we include indeendently with robability each non-constant monomial of total degree at most D in n variables in a generating set of the ideal. In other words, starting with B =, each monomial in S of degree at most D is added to the set B indeendently with equal robability. The resulting random monomial ideal is then simly I = (B); if B =, then we let I = (0). Henceforth, we will denote by B(n, D, ) the resulting Erdős-Rényi-tye distribution on the sets of monomials. Since sets of monomials are now random variables, a random set of monomials B drawn from the distribution B(n, D, ) will be denoted by the standard notation for distributions B B(n, D, ). Note that if B S is any fixed set of monomials of degree at most D each and B B(n, D, ), then P (B = B) = B ( ) (D+n D ) B. In turn, the distribution B(n, D, ) induces a distribution I(n, D, ) on ideals. We will use the notation I I(n, D, ) to indicate that B B(n, D, ) and I = (B) is a random monomial ideal generated by the ER-tye model. Before we state our results, let us establish some necessary robabilistic notation and background. Readers already familiar with random structures and robability background may ski this art, but those interested in further details are directed to many excellent texts on random grahs and the robabilistic method including [, 7, ]. Given a random variable X, we denote its exected value by E [X], its variance by Var [X], and its covariance with random variable Y as Cov [X, Y ]. We will often use in our roofs four well-known facts: the linearity roerty of exectation, i.e., E [X + Y ] = E [X] + E [Y ], the first moment method (a secial case of Markov s inequality), which states that P (X > 0) E [X], and the second moment method (a secial case of Chebyshev s inequality), which states that P (X = 0) Var [X]/(E [X]) for X a non-negative integer-valued random variable. Finally, an indicator random variable A for an event A is a random variable such that A =, if event A occurs, and A = 0 otherwise. Indicator random variables behave articularly nicely with ard to taking exectations: for any event A, E [ A ] = P (A) and Var [ A ] = P (A)( P (A)). Following convention, we abbreviate by saying a roerty holds a.a.s, to mean that a roerty holds asymtotically almost surely if, over a sequence of sets, the robability of having the roerty converges to. An imortant oint is that I(n, D, ) is a arametric robability distributions on the set of monomial ideals, because it clearly deends on the value of the robability arameter. A key concern of our aer is to investigate how invariants evolve as changes. To exress some of our results, we will need to use asymtotic analysis of robabilistic events: For functions f, g : N R, we write f(n) = o (g(n)), and g(n) = ω (f(n)) if lim n f(n)/g(n) = 0. We also write f(n) g(n) if lim n f(n)/g(n) =, and f(n) g(n) if there exist ositive constants n 0, c, c such that c g(n) f(n) c g(n) when n n 0. When a sequence of robabilistic events is given by f(n) for n N, we say that f haens asymtotically almost surely, abbreviated a.a.s., if lim n P (f(n)) =.

3 De Loera, Petrović, Silverstein, Stasi, Wilburne/Random monomial ideals In analogy to grah-theoretic roerties, we define an (monomial) ideal-theoretic roerty Q to be the set of all (monomial) ideals that have roerty Q. A roerty Q is monotone increasing if for any monomial ideals I and J such that I J, I Q imlies J Q as well. We will see that several algebraic invariants on monomial ideals are monotone. Let I I(n, D, ). A function f : N R is a threshold function for a monotone roerty Q if for any : N [0, ]: lim P (I Q) = D { 0, if = o (f(d)),, if = ω (f(d)). A similar definition holds for the case when f is a function of n and n. A threshold function gives a zero/one law that secifies when a certain behavior aears. Our results Our results, summarized in (A)-(E) below, describe the kind of monomial ideals generated by the ER-tye model, in the sense that we can get a handle on both the average and the extreme behavior of these random ideals and how various ranges of the robability arameter control those roerties. The following can also serve as an outline of this aer. (A) Hilbert functions and the distribution of monomial ideals In Section, we begin our investigation of random monomial ideals by showing that I(n, D, ) is not a uniform distribution on all monomial ideals, but instead, the robability of choosing a articular monomial ideal under the ER-tye model is comletely determined by its Hilbert function and the first total Betti number of its quotient ring. For an ideal I, denote by β i,j = β i,j (S/I) the (i, j)-th graded Betti number of S/I, that is, the number of syzygies of degree j at ste i of the minimal free resolution. We will collect the first graded Betti numbers, counting the minimal generators of I, in a vector ˆβ = (β,, β,,..., β,0 ), and denote by β = β, + β, + + β,0 the first total Betti number of S/I. We will denote by h I ( ), or simly h I, the Hilbert function of the ideal I. We can now state the following foundational result: Theorem.. Let I S be a fixed monomial ideal generated in degree at most D. The robability that the random monomial ideal I I(n, D, ) equals I is P (I = I) = β (S/I) ( ) + D d= h I(d). In this way, two monomial ideals with the same Hilbert function and the same number of minimal generators have the same robability of occurring. Then, the following natural question arises: How many monomial ideals in n variables with generators of degree less than or equal to D and with a given Hilbert function are there? It is well-known that one can comute the Hilbert function from the graded Betti numbers (see e.g., [8]). Thus it is no surrise that we can state a combinatorial lemma to count monomial ideals in terms of their graded Betti numbers. If we denote by NMon(n, D, h, ˆβ ) the number of ossible monomial ideals in n variables, generated in degree no more than D, Hilbert function h, and first graded Betti numbers ˆβ = (β,, β,,..., β,0 ), then NMon(n, D, h, ˆβ ) is equal to the number of 0 vertices of a certain convex olytoe (Lemma.). Based on Theorem. and Lemma., we can rovide a closed formula for the induced distribution on Hilbert functions under the ER-tye model, that is, the robability that the ER-tye model laces on any given Hilbert function. This is resented in Theorem..

4 De Loera, Petrović, Silverstein, Stasi, Wilburne/Random monomial ideals (B) The Krull dimension and random monomial ideals Section is dedicated to a first fundamental ring invariant: the Krull dimension dim S/I. By encoding the Krull dimension as the transversal number of a certain hyergrah, we can show (see Theorem. for the recise statement) that the robability that dim(s/i) is equal to t for 0 t n, for a random monomial ideal I I(n, D, ), is given by a olynomial in of degree t+ ( D )( n i= i i). This formula is exonentially large but we make it exlicit in some interesting values of t (Theorem.); which in articular give a comlete descrition of the case for two- and three-variable olynomials. Turning to asymtotic behavior, we rove that the Krull dimension of (S/I) can be controlled by bounding the asymtotic growth of the robability arameter = (D) as D. The evolution of the Krull dimension in terms of is illustrated in Figure. The result is obtained by combining the family of threshold results from Theorem. and is stated in the following corollary: Corollary.. Let I I(n, D, ), n be fixed, and 0 t < n. If the arameter = (D) is such that = ω ( D (t+)) and = o ( D t) as D, then dim(s/i) = t asymtotically almost surely. It is very useful to consider the evolution as the robability increases from 0 to of the random monomial ideal from the ER-tye model and its random generating set B. For very small values of, B is all but guaranteed to be emty, and the random monomial ideal is asymtotically almost surely the zero ideal. As increases, the random monomial ideal evolves into a more comlex ideal generated in increasingly smaller degrees and suort. Simultaneously, as the density of B continues to increase with, smaller-degree generators aear; these divide increasingly larger numbers of monomials and the random ideal starts to become less comlex as its minimal generators begin to have smaller and smaller suort, causing the Krull dimension to dro. Finally, the random ideal becomes 0-dimensional and continues to evolve towards the maximal ideal. n (zero ideal) n n 0 = 0 = D n = D (n ) = D (n ) = D = D = Fig. Evolution of the Krull dimension of S/I, where I is the random monomial ideal I I(n, D, ), as the robability arameter = (D) changes. (C) First Betti numbers and random monomial ideals The ER-tye model induces a distribution on the Betti numbers of the coordinate ring of random monomial ideals. To study this distribution, we first ask: what are the first Betti numbers generated under the model? A natural way to understand a distribution is to comute its exected value, i.e., the average first Betti number. In Theorem. we establish the asymtotic behavior for the exected number of minimal generators in a random monomial ideal. We also establish and quantify threshold behavior of first graded Betti numbers of random monomial ideals I I(n, D, ) in two imes, as n or D go to infinity, in Theorem.. As in the Krull dimension case, we combine the threshold results to obtain the following corollary: Corollary.. Let I I(n, D, ). a) Let D be fixed, and d be a constant such that < d D. If the arameter = (n) is such that (n) = ω ( n d) and (n) = o ( n d+) then initdeg(i) = d asymtotically almost surely. b) Let n be fixed. Suose that d i = d i (D), d i D and lim D d i (D) = for i {, }. If the arameter = (D) is such that (D) = ω ( d n ) and (D) = o ( d n ) as D, then

5 De Loera, Petrović, Silverstein, Stasi, Wilburne/Random monomial ideals d initdeg(i) d asymtotically almost surely. Note that this is attainable when d (D) = ω (d (D)). While Corollary. exlains how the arameter controls the smallest non-zero first Betti number, it gives no information arding the largest such number; the latter, of course, leads the comlexity of the first Betti numbers. To that end, we study the degree comlexity 0 (I), introduced by Bayer and Mumford in [] as the maximal degree of any reduced Gröbner basis of I; for monomial ideals, this is simly the highest degree of a minimal generator. In Theorem.7 we show that can be secified so that, even though B contains many large degree monomials, none are minimal. Intuitively, one can think of the above results in following terms: when the growth of is bounded above by that of /D, Corollary. rovides the minimum degree of a minimal generator of I. When the growth of is bounded below by that of /D, Theorem.7 rovides the maximum degree of a minimal generator of I. The evolution of the minimum degree of a minimal generator for the case when D is fixed and n tends to infinity is illustrated in Figure : zero ideal D D D d = 0 = n D = n (D ) = n (D ) = n (D ) = n d = n (d ) = Fig. Evolution of the initial degree of S/I, where I is the random monomial ideal I I(n, D, ), as the model arameter = (n) changes, for fixed D. The initial degree, that is, the smallest integer k such that β,j = 0 for j < k and β,k 0, is the value dislayed in each box. (D) Other robabilistic models for generating monomial ideals In Section we define more general models for random monomial ideals. This is useful to study distributions on ideals that lace varying robabilities on monomials deending on a roerty of interest. For examle, one may wish to change the likelihood of monomials according to their total degree. To that end, we define a general multiarameter model for random monomials, show how the ER-tye model is a secial case of it, and introduce the graded model. In Theorem. we show there is a choice of robability arameters in the general model for random monomial ideals that recovers, as a secial case, the multiarameter model for random simlicial comlexes of [], which itself generalizes various models for random grahs and clique comlexes resented in [7, ]. (E) Exeriments & conjectures Section contains a summary of comuter simulations of additional algebraic roerties of random monomial ideals using the ER-tye model. For each trile of selected values of model arameters, n, and D, we generate 000 monomial ideals I I(n, D, ) and comute their various roerties. Our simulation results are summarized in Figures,,,, 7, and 8. The exeriments suggest several trends and conjectures, which we state exlicitly.. Hilbert functions and the Erdős-Rényi distribution on monomial ideals The Hilbert function, and the Hilbert olynomial that it determines, are imortant tools in the classification of ideals via the Hilbert scheme [7]. Of course, monomial ideals are key to this enterrise: given an arbitrary homogeneous ideal I and a monomial order on k[x,..., x n ], the monomial ideal generated by the leading terms of the olynomials in I has the same Hilbert function as I. Therefore the study of all Hilbert functions reduces to the study of Hilbert functions of monomial ideals.

6 De Loera, Petrović, Silverstein, Stasi, Wilburne/Random monomial ideals Moreover, the Hilbert function h I ( ) of a monomial ideal I S has a useful combinatorial meaning: the value of h I (d), d D, counts the number of standard monomials - monomials of degree exactly d that are not contained in I. In this section we exlore the Erdős-Rényi distribution of monomial ideals and the induced distribution on Hilbert functions. It turns out they are intimately related... Distribution of monomial ideals Our first theorem recisely describes the distribution I(n, D, ) on ideals secified by the ER-tye model. Recall that Theorem. states that the Hilbert function and first total Betti number determine the robability of any given ideal in n variables generated in degree at most D. Proof of Theorem.. Fix I S generated in degree at most D and let M be the unique minimal set of generators of I. Then, I = I if and only if B M and no monomial x α such that x α I is in B. Let A denote the event that each of the β (S/I) elements of M is in B and let A denote the event that no monomial x α such that x α I is in B. Then, the event I = I is equivalent to the event A A. Since the events A and A are indeendent, P (I = I) = P (A A ) = P (A )P (A ). Observe that P (A ) = β(s/i), since each of the β (S/I) elements of M is chosen to be in B indeendently with robability and P (A ) = ( ) + D d= hi(d), since there are exactly + D d= h I(d) monomials of degree at most D not contained in I and each of them is excluded from B indeendently with robability. Consider a secial case of this theorem in the following examle. Examle. (Princial random monomial ideals in variables). Fix n = and D > 0. In this examle, we calculate the robability of observing the rincial random monomial ideal (B) = (x α y γ ) k[x, y], α + γ D. We observe the ideal (x α y γ ) exactly when x α y γ B and x α y γ B for every monomial x α y γ with α < α D or γ < γ D. Thus, we must count the lattice oints below the staircase whose only corner is at (α, γ) Z 0. A simle counting argument shows that there are exactly γ (D γ + ) + α (D α + ) αγ lattice oints. Hence, P (I = (x α y γ )) = ( ) γ (D γ+)+ α (D α+) αγ. Note that the right hand side of this exression is maximized at α =, γ = 0 or α = 0, γ =, so that among the rincial ideals, the most likely to aear under the ER-tye model are (x) and (y). A corollary of Theorem. follows directly from the well-known formulas for Stanley-Reisner rings: Corollary.. If I is a square-free monomial ideal, then the robability of I under the ER-tye model is determined by the number of minimal non-faces and faces of the associated simlicial comlex... Distribution of Hilbert functions Since the Erdős-Rényi-tye model secifies a distribution on monomial ideals, it rovides a formal robabilistic rocedure to generate random Hilbert functions, in turn roviding a very natural question: what is the induced robability of observing a articular Hilbert function? In other words, what is the most likely Hilbert function under the ER-tye model? If we randomly generate a monomial ideal, are all Hilbert functions equally likely? Not all ossible non-negative vectors are first graded Betti vectors, but Lemma. below gives an algorithmic way to determine, using olyhedral geometry, whether a articular otential set of first

7 De Loera, Petrović, Silverstein, Stasi, Wilburne/Random monomial ideals 7 graded Betti numbers can occur for a monomial ideal in n variables. We note that Onn and Sturmfels introduced other olytoes useful in the study of initial (monomial) ideals of a zero-dimensional ideal for n generic oints in affine d-dimensional sace: the staircase olytoe and the simler corner cut olyhedron []. Lemma.. Denote by NMon(n, D, h) the number of ossible monomial ideals in n variables, with generating monomials of degree no more than D and given Hilbert function h. Then NMon(n, D, h) is equal to the number of vertices of the 0 convex olytoe Q(n, D, h) defined by α=l ( ) n + L x α = h(l), L x α x γ, L =,... D α γ, α + = γ, where α, γ denote exonent vectors of monomials with n variables and total degree no more than D, thus the system has ( ) n+d D variables. Proof of Lemma.. The variables x α are indicator variables recording when a monomial x α, for 0 < α D, is chosen to be in the ideal (x α = ) or is not chosen to be in the ideal (x α = 0). The first tye of equation forces the values of the Hilbert function at degree L to be satisfied. The inequalities x α x γ, for all α γ, ensure the set of chosen monomials is closed under divisibility and thus forms an ideal. A non-redundant subset of these satisfying α + = γ suffices to cut out the olyhedron. Theorem.. Let h be a Hilbert function for an ideal of the olynomial ring in n variables. As before, Q(n, D, h) is the olytoe whose vertices are in one-to-one corresondence with the distinct monomial ideals in n variables, with generators of degree less than or equal to D, and Hilbert function h. Then the robability of generating a random monomial ideal I from the ER-tye model with Hilbert function h is exressed by the following formula P (h I = h) = ( ) + D d= h(d) v V (Q(n,D,h)) β (v), (.) where β (v) is the total first Betti number of the monomial ideal associated to the vertex v. Proof of Theorem.. To derive the above formula for the robability of randomly generating a monomial ideal with a articular Hilbert function h using the ER model, we decomose the random event into robability events that are combinatorially easy to count. The robability of generating random monomial ideals with a given Hilbert function h is the sum of the robabilities of disjoint events enumerated by the vertices v in the olytoe Q(n, D, h) (see Lemma.). For each choice of v, the robability of the corresonding ideal being randomly generated by our model is given in Theorem. and is comletely determined by its Hilbert function and first total Betti number. These ideals share a Hilbert function, h, and so ( ) + D d= h(d) is a common factor for all terms. The choice of v determines the ideal and its first Betti number. Some final remarks about Theorem.. First of all, one can find examles that show that the Hilbert function of an ideal does not uniquely determine the value of the number of minimal generators. Second, the values for the Betti numbers have exlicit tight bounds from the Bigatti-Hulett- Pardue theorem (see [],[], []). Finally, note that, from the oint of view of comuter science, it is useful to know that the Hilbert function h can be secified with one olynomial (the Hilbert olynomial) accomanied by a list of finitely many values that do not fit the olynomial.

8 . Krull dimension De Loera, Petrović, Silverstein, Stasi, Wilburne/Random monomial ideals 8 In this section, we study the Krull dimension dim S/I of a random monomial ideal under the model I(n, D, ). First, in Section., we give an exlicit formula to calculate the robability of I having a given Krull dimension. Second, in Section., we establish threshold results characterizing the asymtotic behavior of the Krull dimension and we rovides ranges of values of that, asymtotically, control the Krull dimension of the random monomial ideal... Hyergrah Transversals and Krull dimension Let us recall some necessary definitions. A hyergrah H is a air H = (V, E) such that V is a finite vertex set and E is a collection of non-emty subsets of V, called edges; we allow edges of size one. A clutter is a hyergrah where no edge is contained in any other (see, e.g., [0, Chater ]). A transversal of H, also called a hitting set or a vertex cover, is a set T V such that no edge in E has a nonemty intersection with T. The transversal number of H, also known as its vertex cover number and denoted c(h), is the minimum cardinality of any transversal of H. Given a collection B of monomials in k[x,..., x n ], the suort hyergrah of B is the hyergrah G(B) with vertex set V = {x,..., x n }, and edge set E = {su(x α ) x α B}, where su(x α ) denotes the suort of the monomial x α, i.e., x i su(x α ), if and only if x i x α. What we require here is that for I = (B), the Krull dimension of S/I is determined by the transversal number of G(B) through the formula dim S/I + c (G(B)) = n. (.) This formula is known to hold for square-free monomial ideals minimally generated by B (see e.g., []). Since the suort hyergrah of B is equal to that of its maximal square-free subset and c (G(B)) is invariant under the choice of generating set B of I, the equation, obviously, also holds for general ideals... Probability distribution of Krull dimension The relationshi between Krull dimension and hyergrah transversals leads to a comlete characterization of the robability of roducing a monomial ideal with a articular fixed Krull dimension in the ER-tye model (Theorem.). Exlicitly comuting the values given by Theorem. requires an exhaustive combinatorial enumeration, which we demonstrate for several secial cases in Theorem.. Nevertheless, where the size of the roblem makes enumeration rohibitive, Theorem. gives that P (dim S/I = t) is always a olynomial in of a secific degree, and thus can be aroximated by numerically evaluating P (dim S/I = t) for a sufficient number of values, and interolating. Theorem.. Let I I(n, D, ). For any integer t, 0 t n, the robability that S/I has Krull dimension t is given by a olynomial in of degree t+ i= ( D i )( n i). More recisely, P (dim S/I = t) = C C n t σ E(C) ( ) ( D σ) σ {x,...,x n} σ σ σ E(C) ( ) ( D σ ), where C n t is the set of all clutters on {x,..., x n } with transversal number n t. Proof. Let B be the random generating set that gave rise to I. By Equation (.), the robability that S/I has Krull dimension t is equal to the robability that c(g(b)) = n t. Let G min (B) be

9 De Loera, Petrović, Silverstein, Stasi, Wilburne/Random monomial ideals 9 the hyergrah obtained from G(B) by deleting all edges of G(B) that are strict suersets of other edges in G(B): in other words, G min (B) contains all edges σ of G(B) for which there is no σ in G(B) such that σ σ. Also let C n t denote the set of clutters on {x,..., x n } with transversal number n t. Then c(g(b)) = n t if and only if G min (B) = C for some C C n t. As these are disjoint events, we have P (dim S/I = t) = Let σ {x,..., x n }. There are ( D σ) monomials suorted on σ, so C C n t P (G min (B) = C). (.) P (σ E(G(B))) = ( ) ( D σ), and P (σ E(G(B))) = ( ) ( D σ). Note that both exressions are olynomials in of degree exactly ( D σ). Now Gmin (B) = C if only if:. every edge σ of C is an edge of G(B), and. every σ {x,..., x n } such that σ does not contain any σ E(C) is not an edge of G(B). Note that condition is the contraositive of: any edge of G(B) which is not an edge of C is not minimal. Since edges are included in G(B) indeendently, the above are also indeendent events and the robability of their intersection is P (G min (B) = C) = P (σ E(G(B))) P ( σ E(G(B)) ) = σ E(C) σ E(C) s.t. σ t+ P (σ E(G(B))) σ {x,...,x n} σ σ σ E(C) σ {x,...,x n} s.t. σ t and σ σ σ E(C) P ( σ E(G(B)) ). (.) The last equality follows because for C C n t, if σ is an edge of C then σ t+, and also if σ is a subset of {x,..., x n } satisfying σ σ for all σ E(C), then σ t. To show the first statement, suose σ E(C) with σ > t +. Since C is a clutter, no roer subset of σ is an edge of C, so for every σ E(C), σ σ, σ contains at least one vertex not in σ. Hence the set T = {x,..., x n }\σ intersects every edge of C excet σ. By taking the union of T with any one vertex in σ, we create a transversal of C of cardinality T + = n σ + < n t, contradicting c(c) = n t. For the second statement, suose σ > t. By assumtion no edge of C is a subset of σ, so every edge of C contains at least one vertex in the set T = {x,..., x n }\σ. Hence T is a transversal of C with T = n σ < n t, again a contradiction. No subset of {x,..., x n } can aear in both index sets of (.), and each subset that does aear contributes a olynomial in of degree ( ) D i. It follows that (.) is a olynomial in of degree no greater than t+ ( D )( n i= i i), and hence so is P (dim S/I = t) as (.) is a sum of such olynomials. To rove this bound is in fact the recise degree of the olynomial, we show there is a articular clutter for which (.) has degree exactly t+ ( D )( n i= i i), and furthermore, that for every other clutter the exression in (.) is of strictly lower degree. Hence the sum over clutters in (.) has no cancellation of this leading term.

10 Consider the hyergrah Kn t+ C n t and Then K t+ n De Loera, Petrović, Silverstein, Stasi, Wilburne/Random monomial ideals 0 P ( G min (B) = Kn t+ ) = σ {x,...,x n} σ=t+ that contains all ( n t+) edges of cardinality t+ and no other edges. P (σ E(G(B))) σ {x,...,x n} σ t P (σ E(G(B))) = ( ( ) ( D t+) ) ( n t+) ( ) Dn+( D )( n )+ +( D t )( n t) which has the correct degree. On the other hand, if C C n t and C Kn t+, then at least one edge σ of C is not an edge of Kn t+ ; hence σ t. All subsets roerly containing σ are neither edges of C, nor do they satisfy condition above, hence these subsets are not indexed by either roduct in.. In articular there are ositively many subsets of cardinality t + which do not contribute factors to P (G min (B) = C). To use the main formula of Theorem., one must be able to enumerate all clutters on n vertices with transversal number n t. When t is very small or very close to n this is tractable, as we see in the following theorem. Theorem.. Let I I(n, D, ). Then, a) P (dim S/I = 0) = ( ( ) D) n. b) P (dim S/I = ) = n ( n j=0 j) ( ( ) D ) j ( ) ( D(n j) ( ) (D ) ) ( n j ). c) P (dim S/I = n ) = ( ) (n+d n ) + n j= ( )j ( ) n j ( ) ( n+d n ) ( n+d j n ). d) P (dim S/I = n) = ( ) (n+d n ). For n, the listed formulas gives the comlete robability distribution of Krull dimension induced by I(n, D, ), for any integer D. Proof. Proof of Part (a) For B B(n, D, ), I = (B), S/I is zero dimensional if and only if G min (B) C n. There is a single clutter on n vertices with transversal number n: the one with edge set {{x }, {x },..., {x n }}. Hence by Theorem., P (dim S/I = 0) = n i= ( ) ( D {x i }) = ( ( ) D ) n. Proof of Part (b): S/I is one-dimensional if and only if G min (B) C n. We wish to describe C n. Suose C is a clutter on n vertices and exactly j of the vertices are contained in a -edge. Then j n else C would be the clutter from art a, so let 0 j n, and denote by V the set of these j vertices. Then V is a subset of any transversal of C. Let W = {x,..., x n }\V, then it can be shown that c(c) = n if and only E(C) = {{x i } x i V } {{x i, x k } x i, x k W, x i x k }. Hence P (G min (B) = C) equals P ({x i } E(G(B))) P ({x i, x k } E(G(B))) P ({x i } E(G(B))) x i V x i,x k W x i W = ( ( ) D ) j ( ( ) (D ) ) ( n j ) ( ) D(n j). The exression for P (dim S/I = ) is obtained by summing over all ( n j) ways of selecting the j -edges, for each 0 j n.

11 De Loera, Petrović, Silverstein, Stasi, Wilburne/Random monomial ideals Proof of Part (c): For the case of (n ) dimensionality, Theorem. requires us to consider clutters with transversal number : clutters where some vertex aears in all the edges. However, for this case we can give a simler argument by looking at the monomials in B directly. Now the condition equivalent to (n ) dimensionality is that there is some x i that divides every monomial in B. Fix an i, then there are ( ) ( n+d D n+d ) D monomials that xi does not divide. If F i is the event that x i divides every monomial in B, then P (F i ) = ( ) (n+d D ) ( n+d D ). To get an exression for (n )-dimensionality, we need to take the union over all F i, which we can do using an inclusionexclusion formula considering the events that two variables divide every monomial, three variables divide every monomial, etc. Finally, we subtract the robability that B is emty. Proof of Part (d): Since only the zero ideal has Krull dimension n, this occurs if and only if B is emty, which has robability ( ) (n+d D )... Threshold and asymtotic results Next we study the asymtotic behavior of dim(s/i). In articular, we derive ranges of values of that control the Krull dimension as illustrated in Figure. As a secial case we obtain a threshold result for the random monomial ideal being zero-dimensional. Lemma.. Let n, t be fixed integers, with 0 t n. If B B(n, D, ) and (D) = o ( D t), then a.a.s. B will contain no monomials of suort size t or less as D tends to infinity. Proof. By the first moment method, the robability that B contains some monomial of suort at most t is bounded above by the exected number of such monomials. As the number of monomials in n variables with suort of size at most t is strictly less than ( )( n D+t ) t t, the exectation is bounded above by the quantity ( )( n D+t ) ( t t. This quantity tends to zero when (D) = o D t ) and n and t are constants, thus establishing the lemma. Theorem.. Let n, t be integers and = (D) with t n and 0. Suose that B B(n, D, ) and I I(n, D, ). Then D t is a threshold function for the roerty that dim(s/i) t. In other words, lim P (dim(s/i) t ) = D {0, if = o ( D t), if = ω ( D t). Proof. Let = o ( D t). Consider the suort hyergrah of B, G(B). If every monomial in B has suort of size t + or more, then any (n t)-vertex set of G(B) is a transversal, giving c(g(b)) n t. Thus Equation. and Lemma. give that dim(s/i) t holds with robability tending to 0. Now, let = ω ( D t). Consider σ {x,..., x n } with σ = t and let X σ be the random variable ) that records the number of monomials in B with suort exactly the set σ. There are exactly ( D t such monomials of degree at most D: each monomial can be divided by all variables in σ and the remainder is a monomial in (some or all of) the variables of σ of degree no greater than D t. Thus, E [X σ ] = ( D t ) tends to infinity with D, when = ω ( D t ). Further, Var [X σ ] E [X σ ], as X σ is a sum of indeendent indicator random variables, one for each monomial with suort exactly σ. Hence we can aly the second moment method to conclude that as D tends to infinity, P (X σ = 0) Var [X σ] E [X σ ] ( D t ) 0.

12 De Loera, Petrović, Silverstein, Stasi, Wilburne/Random monomial ideals In other words, a.a.s. B will contain a monomial with suort σ. Since σ was an arbitrary t-subset of the variables {x,..., x n }, it follows that the same holds for every choice of σ. The robability that any of the events X σ = 0 occurs is bounded above by the sum of their robabilities and tends to zero, as there are a finite number of sets σ. Then, with robability tending to, G(B) contains all size-t edges, which imlies that, for every vertex set of size n t or less, G(B) contains an edge of size t that is disjoint from it. Thus the transversal number of the hyergrah is at least n t +, and, by Formula., dim(s/i) = n c(g(b)) t a.a.s. The smallest choice of t-values in Theorem., t =, gives a threshold function for the event that I I(n, D, ) has zero Krull dimension: Corollary.. Let I I(n, D, ), n be fixed and = (D). Then the function /D is a threshold function for the 0-dimensionality of S/I. Now we move to obtain Corollary. announced in the Introduction. Theorem. establishes a threshold result for each choice of constant t. But, if both = ω ( D (t+)) and = o ( D t) hold, then the theorem gives that events dim(s/i) > t and dim(s/i) < t each hold with robability tending to 0. Therefore the robability of their union, in other words the robability that dim(s/i) is either strictly larger than t or strictly smaller than t, also tends to 0. By combining the threshold results alied to two consecutive constants t and t+ in this way, we have established the Corollary.. Minimal generators of random monomial ideals from the ER-tye model A way to measure the comlexity of a monomial ideal is to study the number of minimal generators. As before, we will denote by β := β (S/I) the first total Betti number of the quotient ring of the ideal I, that is, the total number of minimal generators of I. Subsection. addresses this invariant and rovides a formula for the asymtotic average value of β. A more refined invariant is, of course, the set of first graded Betti numbers β,d (S/I), recording the number of minimal generators of I of degree d. Subsection. rovides threshold results that exlain how the asymtotic growth of the arameter value influences the aearance of minimal generators of given degree. Subsection. rovides an analogous result for the degree comlexity 0 (I) of a monomial ideal I, that is, the maximum degree of a minimal generator... The first total Betti number Consider the random variable β (S/I) for I I(n, D, ) in the ime D. Before stating the results of this section, we introduce a concet from multilicative number theory. The order r divisor function, denoted by τ r (k), is a function that records the number of ordered factorizations of an integer k into exactly r arts (e.g., τ () = ). Theorem.. Let I I(n, D, ) and let β (S/I) be the random variable denoting the number of minimal generators of the random monomial ideal I. Then, lim E [β (S/I)] = D τ n (k + )( ) k. () Proof. For each α = (α,..., α n ) Z n 0, let α be the indicator random variable for the event that x α is a minimal generator of I. Excluding the constant monomial and x α itself, there are k=0

13 De Loera, Petrović, Silverstein, Stasi, Wilburne/Random monomial ideals n i= (α i + ) monomials in S that divide x α. Therefore, E [ α ] = ( ) n i= (αi+) and letting D shows that lim E [β (S/I)] = D 0 α Z n 0 E [ α ] = 0 α Z n 0 ( ) n i= (α i+). Since n i= (α i +) is also the number of lattice oints contained in the box n i= [0, α i], excluding the origin and α, it follows that every integer k 0 aears as an exonent of q in the infinite sum above exactly as many times as k + can be factored into a roduct of n ordered factors, so that the above sum may be rewritten as k=0 τ n(k + )( ) k. When n =, τ (k) is simly the usual divisor function which counts the number of divisors of an integer k. The ordinary generating function of τ (k) has the form of a Lambert series (see []), which leads to the following result in the two-variable case. Corollary.. Let I I(, D, ). Then, lim E [β (S/I)] = D ( ) [ k= Proof. This follows from the Lambert series identity k= ] ( ) k ( ) k. ( ) k ( ) k = k= τ (k)( ) k in (). Corollary. rovides an exression for the asymtotic exected number of minimal generators of a random monomial ideal in two variables that can be easily evaluated for any fixed value of. For examle, if I I(, D, 0 ), then asymtotically (i.e., for very large values of D) one exects I to have about generators on average. On the other hand, when n, the right-hand side of () is difficult to comute exactly, but one can obtain bounds from the fact that n τ n (k) k n, which are valid for every k. Hence n lim D E [β (S/I)] ( ) Li n( ), () where Li s (x) denotes the olylogarithm function of order s defined by Li s (x) = k= xk k. A olylogarithm of negative integer order can be exressed as a rational function of its argument: s n ( ) j+ Li n ( ) = j!s(n +, j + ), j=0 where S(n +, k + ) denotes a Stirling number of the second kind (see [0]). The uer bound in () is thus of order O ( n )... The first graded Betti numbers Here we study the behavior of the random variables β,k (S/I) for I I(n, D, ). Secifically, we show how the choice of in the ER-tye model controls the minimum degree of a generator of I, by controlling the degrees of the monomials in the random set B. Lemma.. Let B B(n, D, ), I = B and B d := {x α B : α d}. a) If d, D is fixed, 0 < d D and = (n), then f(n) = n d is a threshold function for the roerty that B d = 0.

14 De Loera, Petrović, Silverstein, Stasi, Wilburne/Random monomial ideals b) If n is fixed, d = d(d) D is such that lim D d(d) and = (D), then f(d) = d n is a threshold function for the roerty that B d = 0. Proof. For art a), we must show that when (n) = o ( n d), then B d = 0 with robability tending to, and when (n) = ω ( n D), then B d > 0 with robability tending to. For d fixed, the number of monomials of degree at most d in n variables is ( ) n+d d n d (by considering the usual asymtotic bounds for binomial coefficients). For any α Z n 0, with 0 < α d, let α be the indicator random variable for the event that x α B, so that B d = α α, and therefore E [B d ] = ( (n+d α E [ ) ) α] =. By the first moment method: d (( ) ) n + d P (B d > 0) E [B d ] = e d n d. d If (n) = o ( n d), then e d n d 0 as n and hence lim n P (B d > 0) = 0. Now suose (n) = ω ( n d). In this case E [B d ], since: (( ) ) n + d E [B d ] = d ( ) d n d. d The variance of B d is calculated as follows: Var [B d ] = Var [ α ] + (( ) ) n + d Cov [ α, α ] = ( ) = E [B d ] ( ), d α α α where we have again used that Var [ α ] = E [ α] E [α ] = and for α α, α and α are indeendent and thus have 0 covariance. Then by the second moment method and the fact that E [B d ] when = ω ( n d) : P (B d = 0) Var [B d] (E [B d ]) = E [B d ] 0 as n. E [B d ] Part b) follows in a similar fashion, once one notices that, for fixed n, ( ) ( n+d d = n+d ) n d n as d tends to infinity. From Lemma. we immediately obtain a threshold result for the first graded Betti numbers: Theorem.. Let B B(n, D, ) and I = (B). a) If d D is fixed and = (n), then the function f(n) = n d is a threshold function for the roerty that initdeg(i) > d. b) If n is fixed, d = d(d) D is such that lim D d(d) = and = (D), then f(d) = d n is a threshold function for the roerty that initdeg(i) > d. Proof. By definition, initdeg(i) > d if and only if d k= β,k(s/i) = 0. Equivalently, initdeg(i) > d if and only if B and B contains no monomials of degree k for each k d and therefore ( d ) P β,k (S/I) = 0 = P (B d = 0). k= Both a) and b) now follow directly from Lemma..

15 De Loera, Petrović, Silverstein, Stasi, Wilburne/Random monomial ideals The case d = D of Theorem. deserves secial mention, as it rovides a threshold function for the roerty that the random ideal generated by the ER-tye model is the zero ideal. Corollary.. Let B B(n, D, ) and I = B. a) Let D be fixed and = (n). Then n D is a threshold function for the roerty that I = (0). b) Let n be fixed and = (D). Then D n is a threshold function for the roerty that I = (0). On the other hand, in the critical ion of this threshold function, both the exected number of monomials in B and the variance of the number of monomials in B is constant. Proosition.. Let B B(n, D, ). ( (n+d ) ) a) For any n, D, if = / then E [B] = and Var [B] =. D b) For D fixed, if = (n) n D then k lim n E [B] k and c lim n Var [B] c, for some constants k, k, c, c. c) For n fixed, if = (D) D n then k lim D E [B] k and c lim D Var [B] c, for some constants k, k, c, c. Proof. Since B = B D, each of these claims follows directly from the calculated values of E [B D ] and Var [B D ] from the roof of Lemma.. By combining threshold functions from Theorem. we can obtain Corollary. announced in the Introduction: the result secifies ranges of for which initdeg(i) = d asymtotically almost surely. The argument is identical to the one used in Section. to obtain Corollary. from Theorem.... Degree comlexity In the revious subsection, we saw how the choice of influences the smallest degree of a minimal generator of a random monomial ideal. In this subsection, we ask the comlementary question: what can one say about the largest degree of a minimal generator of a random monomial ideal? The following theorem establishes an asymtotic bound for the degree comlexity 0 (I) of I I(n, D, ) for certain choices of robability arameter. Theorem.7. Let n be fixed, I I(n, D, ), = (D) and let r = r(d) be a function tending to infinity as D. If = ω ( r ), then 0 (I) nr a.a.s. Proof. By the roof of Theorem., for each variable x i, a.a.s. B will contain a monomial of the form x j i where j r. If j,..., j n are integers such that j i r for each i n, then the ideal (x j,..., x jn n ) contains all monomials of degree rn in S, since if deg(x α ) r, by the igeonhole rincile at least one α i r n. Hence, with robability tending to, I will contain every monomial in S of degree nr and thus every minimal generator has degree at most nr. To illustrate this result, suose that = ω (/log D) and let I I(n, D, ). Then, by Theorem.7, as D tends to infinity, one should exect 0 (I) to be at most n log D. One should kee in mind that I never contains generators of degree larger than D, and so the bound given in this theorem is not otimal for certain choices of r.

16 . Other robabilistic models De Loera, Petrović, Silverstein, Stasi, Wilburne/Random monomial ideals As mentioned in the Introduction, the study of tyical ideals from a family of interest may require not only the ER-tye model, but also more general models of random monomial ideals. To that end, we define the most general robabilistic model on sets of monomials which, a riori, rovides no structure, but it does rovide a framework within which other models can be recovered. A general model. Fix a degree bound D > 0. To each monomial x α k[x,..., x n ] with 0 < deg(x α ) D, the general model for random monomials assigns an arbitrary robability 0 α of selecting the monomial x α : P (x α ) = α. (.) Hence the general model has many arameters, namely { α : α N n \ {0}, α D}. It is clear that that, for a fixed n and fixed degree bound D, the ER-tye model is a secial case of the general model, where α = (n, D), for all α N n such that 0 < α D, does not deend on α. Of course, there are many other interesting models one can consider. We define another natural extension of ER-tye model. The graded model. Fix a degree bound D > 0. The graded model for random monomials laces the following robability on each monomial x α k[x,..., x n ] with 0 < deg(x α ) D: α = f(n, D, α), (.) where 0 α. Note that, given deg(x α ) = α, the robability of each monomial is otherwise constant in α. The graded model is a D-arameter family of robability distributions on random sets of monomials that induces a distribution on monomial ideals in the same natural way as the ER-tye model. Namely, the robability model selects random monomials and they are added to a generating set with robability according to the model, thus roducing a random monomial ideal in k[x,..., x n ] whose generators are of degree at most D. More recisely, let (n, D) = ( (n, D),..., D (n, D)) be such that 0 i (n, D) < for all i D. Initialize B = and for each monomial x α of degree i, i D, add x α to B with robability i (n, D). The resulting random monomial ideal is I = (B) (again with the convention that if B =, then we set I = (0)). Denote by G(n, D, ) the resulting induced distribution on monomial ideals. As before, since ideals are now random variables, we will write I G(n, D, ) for a random monomial ideal generated using the graded model. The analogue of Theorem. for G(n, D, ), which has an almost identical roof, is as follows: Theorem.. Fix n, D, and the graded model arameters (n, D) = (,..., D ). For any fixed monomial ideal I S, random monomial ideals from the graded model distribution I G(n, D, ) satisfy the following: P (I = I) = D d= β,d(s/i) d ( d ) h I(d), where β,d (S/I) is the number of degree-d minimal generators of I (that is, the first graded Betti number of S/I), and h I (d) is its Hilbert function... Random monomial ideals generalize random simlicial comlexes An abstract simlicial comlex on the vertex set [n] is a collection of subsets (called faces) of [n] closed under the oeration of taking subsets. The dimension dim(σ) of a face σ [n] is σ. If Y is an

17 De Loera, Petrović, Silverstein, Stasi, Wilburne/Random monomial ideals 7 abstract simlicial comlex on [n], then the dimension dim(y ) of Y is max{dim(σ) : σ is a face of Y } and the Stanley-Reisner ideal I Y in the olynomial ring k[x,..., x n ] is the square-free monomial ideal generated by the monomials corresonding to the non-faces τ of Y : I Y := τ : τ [n], τ Y. For examle, if n = and Y = {, {}, {}}, then I Y = x x k[x, x ]. A fundamental result in combinatorial commutative algebra is that the Stanley-Reisner corresondence constitutes a bijection between abstract simlicial comlexes on [n] and the square-free monomial ideals in k[x,..., x n ]. See, for examle, [, Theorem.7]. In [], Costa and Farber introduced a model for generating abstract simlicial comlexes at random. Their model, which we call the Costa-Farber model, is a robability distribution on (r) n, the set of all abstract simlicial comlexes on [n] of dimension at most r. In it, one selects a vector = ( 0,..., r, 0,..., 0) of robabilities such that 0 l for all 0 l r. Then, one retains each of the n vertices with robability 0, and for each air ij of remaining vertices, one forms the edge ij with robability, and for each trile of vertices ijk that are airwise joined by an edge, one forms the dimensional face ijk with robability, and so on. One can study the combinatorial and toological roerties of the random simlicial comlexes generated by the Costa-Farber model. In view of the Stanley-Reisner corresondence, this model can be seen as a model for generating random square-free monomial ideals. Thus, the general model for generating random monomial ideals can be viewed as a generalization of the Costa-Farber model, which in turn has been roved to generalize many other models of random combinatorial objects, for examle, the Erdős-Rényi model for random grahs and Kahle s model for random clique comlexes [7]. See [, Section.] for details. The relationshi between our models for random monomial ideals and the Costa-Farber model can be made recise in the following way: there exists a choice of arameters α in (.) such that the resulting distribution on square-free monomial ideals in S = k[x,..., x n ] is recisely the distribution on the abstract simlicial comlexes on [n] under the Costa-Farber model. Theorem.. Let = ( 0,,..., r, 0,..., 0) denote the n-vector of robabilities in the Costa- Farber model for random simlicial comlexes. Let Y (r) n be a simlicial comlex on [n] of dimension at most r and let I Y be the Stanley-Reisner ideal corresonding to Y. Fix D = r + and secify the following robabilities α, where α N n and 0 < α r +, for the general monomial model (.): { α = deg(x α ), if 0 α {0, } n, (.) 0, otherwise. Then, P CF (Y ) = P (I Y ), where the former is robability under the Costa-Farber model and the latter under the distribution on random monomial ideals induced by the general model (.). In other words, the secification of robabilities in Theorem. recovers the Costa-Farber model on random simlicial comlexes as a sub-model of the model (.) on random monomial ideals. Note that this secification of robabilities can be considered as an instance of the graded model with suort restricted to the vertices of the unit hyercube. Proof. From [, Equation ()], the following robability holds under the Costa-Farber model: P CF (Y ) = r i=0 f i(y ) i ( i ) e i(y ), where f i (Y ) denotes the number of i-dimensional faces of Y and e i (Y ) denotes the number of i- dimensional minimal non-faces of Y (i.e., the number of i-dimensional non-faces of Y that do not