Arxiv:2005.06689V1 [Math.CO]

Home , Stirling number, Stirling polynomials

arXiv:2005.06689v1 [math.CO] 14 May 2020 ubro h eodkind second the of number rmas Partial grammars; re oitrrtthe interpret to order where ti elkon[6 that [16] known well is It tations: r set,1 ascents, are etenme facns(ep ltax ecns of descents) plateaux, (resp. ascents of number the be partial sa xr lmn mle hnalltesin letters all than smaller element extra an as ttsiso h ocle aoiSiln emttosintroduced permutations Jacobi-Stirling so-called the on statistics w nindex an polynomials i +1 e od n phrases. and words Key Let e sfis iea vriwo h aoiSiln emttos Th permutations. Jacobi-Stirling the of overview an give first us Let LTAXO EEAIE TRIGPRUAIN AND PERMUTATIONS STIRLING GENERALIZED ON PLATEAUX , w A S in ytesaitc fpaeu,dset n set r par are ascents and descents plateaux, of statistics the by tions Abstract. ino u eutt h aoiSiln emttoscnrsarec a confirms permutations Jacobi-Stirling the to result our of tion ojcuedet a e n h eodnmdato.Orpar Our author. named second the and Yeh Ma, to due conjecture aso,a ela obntra nepeainfrtecorresp the for interpretation combinatorial a as well as pansion, trigpruain.Bsds eas rvd natraieappro alternative ac an group provide a also and we grammars Besides, γ context-free of permutations. machine Stirling the via obtained γ i pstvt rmtesaiiyo eti utvraepolynomials. multivariate certain of stability the from -positivity ea lhbtwoeeeet r oal ree.Fraword a For ordered. totally are elements whose alphabet an be pstvt ojcueo a aadYh[6 ocrigtestud the concerning [26] Yeh and Ma Ma, of conjecture -positivity n w > i 0 , stesto l emttoso [ of permutations all of set the is A i ≤ , +1 n r ltaxad3 and plateaux are 4 ( i ,weew s h ovninthat convention the use we where ), t γ r endby defined are ) ≤ -positivity. epoeta h nmrtv oyoil fgnrlzdSiln per Stirling generalized of polynomials enumerative the that prove We n san is , eododrElra polynomials Eulerian second-order HCN I,JNM,ADPII .ZHANG B. PHILIP AND MA, JUN LIN, ZHICONG ltax trigpruain;Jcb-trigplnmas Contex polynomials; Jacobi-Stirling permutations; Stirling Plateaux; A n S ascent ( ( X t k PARTIAL ,k n, a eitrrtda h ecn oyoil vrpermu- over polynomials descent the as interpreted be can ) ≥ 0 S nmrtstestpriin f[ of partitions set the enumerates ) X k rs.a (resp. ( , ≥ A n 1. r ecnsof descents are 5 1 n + k ( Introduction t n ,k k, = ) t k γ plateau = ) -POSITIVITY π t X 1 ∈ k n (1 A S = := ] n o ntne if instance, For . A − t n (1 des( a , t ( ) t { − C n ) π descent w w +1 1 ) n , , t 0 e asc( Let . ( ) 2 . t 2 = n , . . . , w ) n C +1 hsppri oiae ya by motivated is paper This . w n ( of ) , n t tial +1 pern as appearing ) onding } n partial ent ealta the that Recall . w w γ = w pstv.Specializa- -positive. rs.plat( (resp. ) tial ino generalized on tion if n[19]. in .Here 0. ˆ 11211 = c otepartial the to ach γ w n cecet,are -coefficients, γ w i with ] e pstvt ex- -positivity w < = lsia Eulerian classical γ -positivity fteethree these of y w sconsidered is 0 ˆ i ∈ +1 1 k · · · P muta- w (resp. 5 lcs In blocks. hn0 then , ,des( ), w n Stirling A ∈ w t-free w i , = )) n 2 , 2 Z. LIN, J. MA, AND P.B. ZHANG

Gessel and Stanley [20] invented the Stirling permutations. A Stirling permutation of order n is a permutation of the multiset {1, 1, 2, 2,...,n,n} such that for each i ∈ [n], all entries between the two occurrences of i are larger than i. Gessel and Stanley [20] provided three diﬀerent proofs for the interpretation

des(π) Cn(t)= t , π∈Q Xn where Qn denotes the set of Stirling permutations of order n. Interestingly, the three statistics asc, plat and des are equidistributed over Qn, as was shown by Bóna in [5] where the statistic plat was first considered. The Jacobi-Stirling numbers JS(n, k; z), as a generalization of the Stirling number S(n, k), were introduced in the study of a problem involving the spectral theory of powers of the classical second-order Jacobi differential expression (see [1]). Write the Jacobi-Stirling n polynomial JS(n + k,k; z) as pn,0(k)+ pn,1(k)z + ··· + pn,n(k)z . Generalizing the above study for S(n+k,k), Gessel, Lin and Zeng [19] investigated the diagonal generating function A (t) p (k)tk = n,i n,i (1 − t)3n+1−i k X≥0 and showed that An,i(t) is the descent polynomial over the Jacobi-Stirling permutations J SP n,i defined in the same flavor as the Stirling permutations. Introduce the multiset

Mn := {1, 1, 1¯, 2, 2, 2¯, . . . , n, n, n¯}, where the elements are ordered by 1¯ < 1 < 2¯ < 2 < ··· < n¯ < n.

Let [¯n] := {1¯, 2¯,..., n¯}. For any subset S ⊆ [¯n], we set Mn,S = Mn \ S. A permutation of Mn,S is a Jacobi-Stirling permutation if for each i ∈ [n], all entries between the two occurrences of i are larger than i. We denote by J SP n,S the set of all Jacobi-Stirling permutations of Mn,S and set

J SP n,i = J SP n,S. S⊆[¯n] |[S|=i

Note that J SP n,n = Qn. Gamma-positive polynomials arise frequently in enumerative combinatorics and have recent motivation coming from geometry; see the survey of Athanasiadis [2]. A univariate polynomial f(x) is said to be γ-positive if it can be expanded as

n ⌊ 2 ⌋ k n−2k f(x)= γkx (1 + x) k X=0 with γk ≥ 0. A bivariate polynomial h(x, y) is said to be homogeneous γ-positive, if h(x, y) is homogeneous and h(x, 1) is γ-positive. This is equivalent to say that h(x, y) can be STIRLING PERMUTATIONS AND PARTIAL γ-POSITIVITY 3 expressed as n ⌊ 2 ⌋ k n−2k h(x, y)= γk(xy) (x + y) k X=0 with γk ≥ 0. One of the typical examples due to Foata and Sch¨uzenberger [16] is the bivariate Eulerian polynomial asc(π) des(π) An(x, y) = π∈Sn x y . A trivariate polynomial i p(x, y, z) = i si(x, y)z is said to be partial γ-positive if every si(x, y) is homogeneous γ-positive. Based on empirical evidence, Ma,P Yeh and the second named author [26] pro- posed the followingP conjecture. Conjecture 1.1 (Ma, Ma and Yeh [26]). Let

asc(π) des(π) plat(π) JSPn,i(x, y, z) := x y z π ∈JX SPn,i be a trivariate extension of An,i(y). For any n ≥ 1 and 0 ≤ i ≤ n, the polynomial JSPn,i(x, y, z) is partial γ-positive. This conjecture has already been confirmed by Ma et al. [26] in the special i = 0 and i = n cases. As a main contribution of this paper, we will prove a generalization of Conjecture 1.1 for the Stirling permutations of a fixed multiset. We need further definitions before we can state our main result. n For each vector m = (m1, m2,...,mn) ∈ P , denote by Mm the general multiset m1 m2 mn {1 , 2 ,...,n }, where i appears mi times. A permutation of Mm is a generalized Stirling permutation if all entries between any two occurrences of i are larger than i for each i ∈ [n]. Let Qm be the set of all generalized Stirling permutations of Mm. Note that Qm = Qn when m = (2, 2,..., 2). The generalized Stirling permutations and various statistics over them have been studied in [10,21,22]. In particular, Brenti [10] showed that n the descent polynomial over Qm has only real roots for each m ∈ P . We will consider more statistics on Qm. n Definition 1.2 (Statistics on Qm). Let π = π1π2 ··· πm ∈ Qm, where m = i=1 mi. As usual, we set π0 = πm+1 = 0. A letter k ∈ [n] is said to be multiple if mk > 1; and single, P otherwise. An index i ∈ [m] is called a multiple (resp. single) descent of π if πi > πi+1 and πi is multiple (resp. single). A double-ascent (resp. double-descent, peak, ascent-plateau, descent-plateau) of π is an index i ∈ [m] such that πi−1 < πi < πi+1 (resp. πi−1 > πi > πi+1, πi−1 < πi > πi+1, πi−1 < πi = πi+1, πi−1 > πi = πi+1). It is clear that if i is a peak, then πi must be single. We further distinguish a double-descent i to be single or multiple according to πi is single or multiple. A descent-plateau i is free if there does not exist an integer ℓ, 1 ≤ ℓ < i, such that πℓ = πi. A plateau i of π is said to be unmovable if i is neither a free descent-plateau nor an ascent-plateau. Let us introduce the statistics of π by

• dasc(π) = #{i ∈ [m]: πi−1 < πi < πi+1}, the number of double-ascents of π; • sddes(π) := #{i ∈ [m]: πi−1 > πi > πi+1 and πi is single}, the number of single double-descents of π; • fdesp(π) := #{i ∈ [m]: πi−1 > πi = πi+1 and πℓ =6 πi for all ℓ ∈ [i − 1]}, the number of free descent-plateaus of π; 4 Z. LIN, J. MA, AND P.B. ZHANG

• ascpp(π) := #{i ∈ [m]: πi−1 < πi = πi+1 or πi−1 < πi > πi+1}, the number of ascent-plateaus and peaks. • mdup(π) := #{i ∈ [m]: i is a multiple descent or a unmovable plateau}. For example, if π = 15565333124411, then dasc(π) = 2, sddes(π) = 0, fdesp(π) = 1, ascpp(π) = 4 and mdup(π) = 4. Now we are ready to state our main result.

n Theorem 1.3. For any m ∈ P with m1 + m2 + ··· mn = m, let

asc(π) des(π) plat(π) Sm(x, y, z)= x y z . π m X∈Q The polynomial Sm(x, y, z) is partial γ-positive and has the expansion

m+1−i m−1 ⌊ 2 ⌋ i j m+1−i−2j Sm(x, y, z)= z γm,i,j(xy) (x + y) , i=0 j=1 X X where

(1.1) γm,i,j = #{π ∈ Qm : sddes(π) = fdesp(π)=0, mdup(π)= i, ascpp(π)= j}. Remark 1.4. To see that Theorem 1.3 implies Conjecture 1.1, for any S ⊆ [n] with |S| = i, deﬁne m(S)=(m1,...,m2n−i) where

2, if ℓ = p + |{a ∈ [n] \ S : a ≤ p}| for some 1 ≤ p ≤ n; mℓ = (1, otherwise. For instance, if S = {1, 2, 5, 7} ⊆ [7], then m(S)=(2, 2, 1, 2, 1, 2, 2, 1, 2, 2). Let

asc(π) des(π) plat(π) JSPn,S(x, y, z) := x y z . π ∈JX SPn,S It is routine to check that JSPn,S(x, y, z) = Sm(x, y, z) with m = m(S). Therefore, by Theorem 1.3, JSPn,S(x, y, z) is partial γ-positive, namely,

3n−i+1−k 3n−i−1 ⌊ 2 ⌋ k j 3n−i+1−k−2j JSPn,S(x, y, z)= z γm,k,j(xy) (x + y) . k j=1 X=0 X Conjecture 1.1 then follows from JSPn,i(x, y, z)= S∈[n] JSPn,S(x, y, z). |S|=i The rest of this paper is organized as follows.P In Section 2, we provide a proof of Theorem 1.3 by using the context-free grammars and a generalization of the Foata–Strehl action on Stirling permutations. In Section 3, an analytic approach to the partial γ- positivity of Sm(x, y, z) is studied by using the stable theory of multivariate polynomials developed by Bocrea and Brändén. As a result, the descent polynomials over Stirling permutations with fixed number of plateaux are shown to be real-rooted. STIRLING PERMUTATIONS AND PARTIAL γ-POSITIVITY 5

2. Proof of Theorem 1.3 2.1. Context-free grammars and an equidistribution. For a set V = {x, y, z, . . .} of commutative variables, a context-free grammar G is a set of substitution rules that replace a variable in V by a Laurent polynomial of variables in V . The formal derivative D associated with a context-free grammar G (introduced by Chen in [11]) is defined by D(x)= G(x) for any x ∈ V and satisfies the following relations: D(u + v)= D(u)+ D(v), D(uv)= D(u)v + uD(v), where u and v are two Laurent polynomials of variables in V . The context-free grammars have been found useful in studying various combinatorial structures [11–15,26,27], including permutations, increasing trees, labeled rooted trees and set partitions. For example, if V = {x, y} and G = {x → xy,y → xy}, then D(x) = xy, D2(x) = xy(x + y) and D3(x)= D(xy)(x + y)+ D(x + y)xy = x3y +4x2y2 + xy3. This is the grammar introduced n by Dumont [15] to generate the bivariate Eulerian polynomials D (x)= An(x, y). n Let π = π1π2 ··· πm ∈ Qm, where m = i=1 mi. A plateau i of π is called a first plateau if πj =6 πi for all 1 ≤ j < i. Let fplat(π) be the number of first plateau of π. Denote by sdes(π) (resp. mdes(π)) the number of singleP (resp. multiple) descents of π. We will apply the context-free grammars to prove the following equidistribution.

n Lemma 2.1. For any m ∈ P with m1 + m2 + ··· + mn = m, the following two triplets of statistics are equidistributed on the Stirling permutations Qm: (des, plat, asc) and (fplat + sdes, mdup, asc). Proof. Note that i is a first plateau if and only if i is either a free descent-plateau or an ascent-plateau. In other words, a plateau is unmovable if and only if it is not a first plateau. Thus, plat(π) − fplat(π) = uplat(π), the number of unmovable plateau of π. We aim to show that the quintuplets (2.1) (sdes, mdes, fplat, uplat, asc) and (sdes, fplat, mdes, uplat, asc) are equidistributed on Qm, from which the lemma follows. For a Stirling permutation π ∈ Qm, we first introduce a grammatical labeling of π as follows:

(L1) If i is a single descent, then put a superscript label x right after πi; (L2) If i is a multiple descent, then put a superscript labelx ˜ right after πi; (L3) If i is a ﬁrst plateau, then put a superscript labely ˜ right after πi; (L4) If i is a unmovable plateau, then put a superscript label y right after πi; (L5) If i is an ascent, then put a superscript label z right after πi.

Recall that we always set π0 = πm+1 = 0. For example, if π = 15565333124411, then the labeling of π is 0z1z5y˜5z6x5x˜3y˜3y3x˜1z2z4y˜4x˜1y1x˜0. It is clear that the weight xsdes(π)x˜mdes(π)y˜fplat(π)yuplat(π)zasc(π) is the product of all the su- perscripts in the labelings of π. Let V = {x, x,˜ y, y,˜ z}. For integer k ≥ 2, introduce the 6 Z. LIN, J. MA, AND P.B. ZHANG context-free grammar k−2 k−2 k−2 k−2 k−2 Gk = {x → x˜yy˜ z, x˜ → x˜yy˜ z,y → x˜yy˜ z, y˜ → x˜yy˜ z, z → x˜yy˜ z}. Also, deﬁne the context-free grammar

G1 = {x → xz, x˜ → xz,y → xz, y˜ → xz, z → xz}.

For any k ≥ 1, let Dk be the formal derivative associated with the context-free grammar Gk. We claim that sdes(π) mdes(π) fplat(π) uplat(π) asc(π) Dmn Dmn−1 ··· Dm1 (z)= x x˜ y˜ y z . π∈Qm X The equidistribution (2.1) then follows from this claim and the fact that the context-free grammars Gk are symmetric inx ˜ andy ˜ for all k ≥ 1. It remains to show the claim. We proceed by induction on n. The statement is obvi- z ously true for the initial case n = 0, as Q0 = {0 0}. Suppose that we have all labeled ′ permutations in Qm′ , where m =(m1, m2,...,mn−1). Note that very permutation in Qm mn can be constructed from a permutation σ ∈ Qm′ by inserting n , mn copies of n, to the position between σi and σi+1 for some nonnegative integer i. The changes of labelings are illustrated as follows for any v ∈{x, x,˜ y, y,˜ z} (no matter what the label v is):

• If mn ≥ 2, then v z y˜ y y x˜ ··· σi σi+1 · · · 7→ · · · σi n n ··· n n σi+1 ··· , y where n appears (mn − 2) times. • Otherwise mn = 1, and v z x ··· σi σi+1 · · · 7→ · · · σi n σi+1 ··· .

mn In either case, the insertion of n corresponds to one substitution rule in Gmn . Since the action of Dmn on elements of Qm′ generates all elements of Qm, the claim holds. This completes the proof of the lemma. In order to ﬁnish the proof of Theorem 1.3, we need a generalization of the so-called Foata–Strehl group action on Stirling permutations that will be introduced below. 2.2. Generalized Foata–Strehl actions on Stirling permutations. The Foata–Strehl group action on permutations was invented by Foata and Strehl [17] (see also [28]) to prove combinatorially the homogeneous γ-positivity expansion

n+1 ⌊ 2 ⌋ k n+1−2k An(x, y)= γn,k(xy) (x + y) , k X=1 where γn,k enumerates the permutations in Sn with k descents and with no double descents. Since then some generalizations and analogues of this γ-positivity expansion, with or without combinatorial proofs, have been found [7, 18, 23–25]. In particular, although the γ-coefficients of the double Eulerian polynomials were known [24] to be nonnegative, to find a combinatorial interpretation of them is still widely open. Here we develop a generalization of the Foata–Strehl action for Stirling permutations to get the combinatorial interpretation of γm,i,j in (1.1), finalizing the proof of Theorem 1.3. STIRLING PERMUTATIONS AND PARTIAL γ-POSITIVITY 7

6s ❅ 5s 5s ❅s5 ✁ 4s 4s ✛✁ ❅❅3s s3 3s ✁ ❅ ❅ 2s ✲❅ ✁ ❅ 1s✁ ✲ ❅1s ❅❅1s s1 ✁ ❆ ✁ ❆ 0 0

Figure 1. Generalized Foata–Strehl actions on 15565333124411

For a Stirling permutation π ∈ Qm and a value x ∈ [n], suppose that πℓ is the leftmost occurrence of x in π. We call x a free descent-plateau value (resp. a single double descents value, a double ascents value) of π if ℓ is a free descent plateau (resp. a single double descents, a double ascents) of π. Then we introduce the Generalized Foata-Strehl action (GFS-action for short) ϕx as follows:

• If x is a free descent-plateau value or a single double descents value of π, then ϕx(π) is obtained from π by moving πℓ to the right of the letter πk, where k = max{a : 0 ≤ a ≤ ℓ − 2, πa < x}; • If x is a double ascents value of π, then ϕx(π) is obtained from π by moving πℓ to the left of the letter πk, where k = min{a : ℓ +2 ≤ a ≤ m +1, πa ≤ x}; • If x is not in the above two cases, then let ϕx(π)= π. The GFS-action has a nice visualization as depicted in Fig. 1. For example, if π = 15565333124411, then ϕ1(π) = 55653331124411, ϕ3(π)=13556533124411 and ϕ2(π) = 15565333144211. When m1 = m2 = ··· = mn = 1, the GFS-action becomes the version of Modiﬁed Foata-Strehl action introduced in [23]. We are now ready for the proof of Theorem 1.3. Proof of Theorem 1.3. By Lemma 2.1, we have asc(π) fplat(π)+sdes(π) mdup(π) Sm(x, y, z)= x y z . π∈Qm X Let us deﬁne the set Qm,i = {π ∈ Qm : mdup(π)= i}. Then Theorem 1.3 is equivalent to (2.2) xasc(π)yfplat(π)+sdes(π) = (xy)ascpp(π)(x + y)m+1−i−2×ascpp(π),

π∈Qm,i ^ X π∈XQm,i ] where Qm,i := {π ∈ Qm,i : sddes(π) = fdesp(π)=0}. Clearly, the GFS-actions ϕx’s are involutions and commute. Thus, for any S ⊆ [n] we can deﬁne the function ϕS : Qm → Qm by ϕS = x∈S ϕx, where the product is the functional compositions. For instance, continuing with the example in Fig. 1, we have Zn Q ϕ{1,3}(π) = 35565331124411. Hence the group 2 acts on Qm via the function ϕS. Since the statistic “mdup” is invariant under this group action, it divides Qm,i into some disjoint Zn orbits. For each π ∈ Qm,i, let Orb(π) = {g(π): g ∈ 2 } be the orbit of π under the GFS-action. Note that x is a free descent-plateau value or a single double descents value of π if and only if x is a double ascents value of ϕx(π). Thus, there exists a unique Stirling 8 Z. LIN, J. MA, AND P.B. ZHANG ] permutationπ ˜ in Orb(π) such that sddes(˜π) = fdesp(˜π) = 0, that is, Qm,i Orb(π)= {π˜}. Therefore, we have T xasc(σ)yfplat(σ)+sdes(σ) = xasc(˜π)−dasc(˜π)yfplat(˜π)+sdes(˜π)(x + y)dasc(˜π) σ∈XOrb(π) =(xy)ascpp(˜π)(x + y)m+1−i−2×ascpp(˜π), where the second equality follows from the following relationships (2.3) asc(˜π) − dasc(˜π) = fplat(˜π) + sdes(˜π) = ascpp(˜π) and (2.4) dasc(˜π)= m +1 − i − 2 × ascpp(˜π).

Summing over all orbits of Qm,i under the GFS-action then gives (2.2). It remains to show the above relationships. Sinceπ ˜ has neither free descent-plateau nor single double descents, every ﬁrst plateau must be an ascent-plateau and every single descent is a peak. Thus, we have fplat(˜π) + sdes(˜π) = ascpp(˜π). As each ascent is followed immediately by an ascent, a plateau or a descent, we have asc(˜π) = dasc(˜π) + ascpp(˜π), which proves (2.3). Clearly, we have asc(˜π) + plat(˜π) + des(˜π) = m + 1 and mdup(˜π)+ fplat(˜π) + sdes(˜π) = plat(˜π) + des(˜π). It then follows that mdup(˜π) + asc(˜π) + fplat(˜π)+ sdes(˜π)= m + 1 and so by (2.3), dasc(˜π) = asc(˜π) − ascpp(˜π)= m +1 − mdup(˜π) − 2 × ascpp(˜π), which is relationship (2.4). This completes the proof of the theorem.

3. Stability

In this section, we shall study the analytic property of the trivariate polynomial Sm(x, y, z). From the involution

π1π2 ··· πm 7→ πmπm−1 ··· π1, we know that Sm(x, y, z) is symmetric in x and y. Hence, the partial γ-positivity of Sm(x, y, z) is equivalent to the γ-positivity of the following reﬁned descent polynomials

des(π) Sm,i(x)= x , π X where the sum runs over all Stirling permutations π ∈ Qm with plat(π) = i. It is known (see [8, Remark 7.3.1]) that the real-rootedness of a polynomial with symmetric coefficients implies the γ-positivity of such polynomial. This motivates us to study an alternative approach to the partial γ-positivity of Sm(x, y, z). Let us recall the stability of multivariate polynomials, which has been developed enor- mously by Bocrea and Brändén [3,4]. A polynomial f ∈ R[x1,...,xn] is said to be stable if either f(x1,...,xn) =6 0 whenever Im(xi) > 0 for all i or f is identically zero. Note that a univariate real polynomial is stable if and only if it has only real roots. Several multivariate Eulerian polynomials have been shown to be stable, see [6, 9, 13, 29, 31]. Our main result of this section is stated as follows. STIRLING PERMUTATIONS AND PARTIAL γ-POSITIVITY 9

Theorem 3.1. Let p(x, y, z) be a trivariate homogeneous polynomial with nonnegative coef- ﬁcients. If p(x, y, z) is stable and symmetric in x and y, then p(x, y, z) is partial γ-positive. In order to prove Theorem 3.1, we shall introduce some basic results about stability- preserving linear operators. Lemma 3.2 (See [30, Lemma 2.4]). Given i, j ∈ [n], the following operations preserve stability of f ∈ R[x1,...,xn]:

(1) Diﬀerentiation: f 7→ ∂f/∂xi.

(2) Diagonalization: f 7→ f|xi=xj . R (3) Specialization: for a ∈ , f 7→ f|xi=a. Now it is time for us to prove Theorem 3.1.

d i Proof of Theorem 3.1. Suppose that p(x, y, z)= i=0 si(x, y)z . We ﬁrst prove that for each 0 ≤ k ≤ d, the polynomial sk(x, y) is stable. By taking the k-th order partial derivative with respect to z of the real stable polynomial p(x,P y, z) it follows from Lemma 3.2 that d i−k i=k(i)ksi(x, y)z is stable, where (i)k = i(i−1) ··· (i−k+1). Note that if z is in the upper half plane then so is −1/z. Hence, we obtain the stability of d (−1)i−k(i) s (x, y)zd−i. P i=k k i Similarly, by taking the (d − k)-th order the partial derivative with respect to z, we get the P stability of sk(x, y). We next prove that sk(x, y) is homogeneous γ-positive. By Lemma 3.2, we get that sk(x, 1) is real-rooted. Since sk(x, y) is symmetric in x and y, we know that the coeﬃcients of sk(x, 1) are symmetric. Hence, the polynomial sk(x, 1) is γ-positive. Since sk(x, y) is homogeneous in x and y, we obtain the homogeneous γ-positivity of sk(x, y). This completes the proof.

We proceed to use Theorem 3.1 to prove the partial γ-positivity of Sm(x, y, z). In order to do this, it suﬃces to show that the stability Sm(x, y, z). For an element π ∈ Qm, let D(π) and A(π) be the set of descents and ascents of π, respectively. If κ = maxi mi and 1 ≤ j < κ, deﬁne Pj(π) to be the set of indices i such that πi = πi+1 where π1 ··· πi−1 contains j − 1 instances of πi. Haglund and Visontai [21, Theorem 3.5] showed that the multivariate polynomial κ−1

xπi yπi zj,πi π∈Qm j=1 X i∈DY(π) i∈AY(π) Yi∈PYj (π) is stable. By diagonalizing the variables xπi , yπi and zj,πi to x, y and z respectively, it follows from Lemma 3.2 that the polynomial Sm(x, y, z) is stable. As an application of this result we get the real-rootedness of Sm,i(x). n Corollary 3.3. For any m ∈ P with m1 + m2 + ··· mn = m and any 0 ≤ i ≤ m − 1, the polynomial Sm,i(x) has only real roots.

Acknowledgments The ﬁrst author was supported by the National Science Foundation of China grant 11871247 and the project of Qilu Young Scholars of Shandong University. The second 10 Z. LIN, J. MA, AND P.B. ZHANG author was supported by the National Science Foundation of China grant 11571235. The third author was supported by the National Science Foundation of China grant 11701424.

References [1] G.E. Andrews, E. Egge, W. Wolfgang and L.L. Littlejohn, The Jacobi-Stirling numbers, J. Combin. Theory Ser. A, 120 (2013), 288–303. [2] C.A. Athanasiadis, Gamma-positivity in combinatorics and geometry, Sém. Lothar. Combin. 77 (2018), Article B77i, 64pp (electronic). [3] J. Borcea and P. Brändén, The Lee-Yang and Pólya-Schur programs. I. Linear operators preserving stability, Invent. Math., 177 (2009), 541–569. [4] J. Borcea and P. Brändén, The Lee-Yang and Plya-Schur programs. II. Theory of stable polynomials and applications, Comm. Pure Appl. Math., 62 (2009), 15951631. [5] M. Bona, Real zeros and normal distribution for statistics on Stirling permutations defined by Gessel and Stanley, SIAM J. Discrete Math., 23 (1) (2008/2009), 401–406. [6] P. Brändén, M. Leander, M. Visontai, Multivariate Eulerian polynomials and exclusion processes. Combin. Probab. Comput., 25 (2016), 486499. [7] P. Brändén, Actions on permutations and unimodality of descent polynomials, European J. Combin., 29 (2008), 514–531. [8] P. Brändén, Unimodality, log-concavity, real-rootedness and beyond, Handbook of Enumerative Com- binatorics, CRC Press Book. (arXiv:1410.6601) [9] P. Brändén, J. Haglund, M. Visontai, and D. G. Wagner, Proof of the monotone column per- manent conjecture, in Notions of positivity and the geometry of polynomials, Trends Math., Birkhäuser/Springer Basel AG, Basel, 2011, 63–78. [10] F. Brenti, Unimodal, Log-concave and Pólya frequency sequences in combinatorics, Mem. Amer. Math. Soc., 413 (1989). [11] W.Y.C. Chen, Context-free grammars, differential operators and formal power series, Theoret. Com- put. Sci., 117 (1993), 113–129. [12] W.Y.C. Chen and A.M. Fu, Context-free grammars for permutations and increasing trees, Adv. in Appl. Math., 39 (2017), 58–82. [13] W.Y.C. Chen, R.X.J. Hao, H.R.L. Yang, Context-free grammars and multivariate stable polynomials over Stirling permutations, arXiv:1208.1420. [14] W.Y.C. Chen and H.R.L. Yang, A context-free grammar for the Ramanujan-Shor polynomials, Adv. in Appl. Math., to appear (2019). [15] D. Dumont, Grammaires de William Chen et dérivations dans les arbres et arborescences, Sém. Lothar. Combin., 37 (1996), Art. B37a. [16] D. Foata and M.-P. Schüzenberger, Théorie Géométrique des Polynômes Eulériens, Lecture Notes in Mathematics, Vol. 138, Springer-Verlag, Berlin, 1970. [17] D. Foata and V. Strehl, Rearrangements of the symmetric group and enumerative properties of the tangent and secant numbers, Math. Z., 137 (1974), 257–264. [18] S. Fu, Z. Lin and J. Zeng, On two unimodal descent polynomials, Discrete Math., 341 (2018), 2616– 2626. [19] I. Gessel, Z. Lin and J. Zeng, Jacobi-Stirling polynomials and P -partitions, European J. Combin., 33 (2012), 1987–2000. [20] I. Gessel and R.P. Stanley, Stirling polynomials, J. Combin. Theory Ser. A, 24 (1978), 24–33. [21] J. Haglund, M. Visontai, Stable multivariate Eulerian polynomials and generalized Stirling permutations, European J. Combin., 33 (2012), 477–487. [22] S. Janson, M. Kuba and A. Panholzer, Generalized Stirling permutations, families of increasing trees and urn models, J. Combin. Theory Ser. A, 118 (2011), 94–114. [23] Z. Lin and D. Kim, A sextuple equidistribution arising in pattern avoidance, J. Combin. Theory Ser. A, 155 (2018), 267–286. [24] Z. Lin, Proof of Gessel’s γ-positivity conjecture, Electron. J. Combin., 23 (2016), #P3.15. STIRLING PERMUTATIONS AND PARTIAL γ-POSITIVITY 11

[25] Z. Lin and J. Zeng, The γ-positivity of basic Eulerian polynomials via group actions, J. Combin. Theory Ser. A, 135 (2015), 112–129. [26] S.-M. Ma, J. Ma and Y.-N. Yeh, γ-positivity and partial γ-positivity of descent-type polynomials, J. Combin. Theory Ser. A, 167 (2019), 257–293. [27] S.-M. Ma, J. Ma and Y.-N. Yeh and B.-X. Zhu, Context-free grammars for several polynomials associated with Eulerian polynomials, Electron. J. Combin., 25(1) (2018), #P1.31. [28] L. Shapiro, W.-J. Woan and S. Getu, Runs, slides and moments, SIAM J. Algebraic Discrete Methods, 4 (1983), 459–466. [29] M. Visontai and N. Williams, Stable multivariate W -Eulerian polynomials, J. Combin. Theory Ser. A, 120 (2013), 1929–1945. [30] D.G. Wagner, Multivariate stable polynomials: theory and applications, Bull. Amer. Math. Soc., 48 (2011), 53–84. [31] P.B. Zhang and X. Zhang, Multivariate stable Eulerian polynomials on segmented permutations, European J. Combin., 78 (2019), 155-162.

(Zhicong Lin) Research Center for Mathematics and Interdisciplinary Sciences, Shandong University, Qingdao 266237, P.R. China E-mail address: [email protected]

(Jun Ma) Department of mathematics, Shanghai Jiao Tong University, Shanghai 200240, P.R. China E-mail address: [email protected]

(Philip B. Zhang) College of Mathematical Science, Tianjin Normal University, Tianjin 300387, P.R. China E-mail address: [email protected]