The unimodality conjecture for the Boţ-Nguyen coefficients

Less than 1 year old · traced to

Let calpha>2calpha>2, and for each nn consider the coefficient sequence (cn,k)k∈{1,2,…,n}(c_{n,k})_{k\in\{1,2,\ldots,n\}}. Unimodality conjecture. For every choice of calpha>2calpha>2, the sequence

(cn,k)k∈{1,2,…,n}(c_{n,k})_{k\in\{1,2,\ldots,n\}}

is unimodal. Numerical experiments motivate this claim, while the preceding discussion shows that including k=0k=0 would fail to give unimodality for small calphacalpha; no resolution is supplied.

References

Primary source

Heinz H. Bauschke and Yuan Gao, “Boţ-Nguyen Acceleration, Weighted Mean Ergodic Iteration, and the Beta-Binomial Distribution”, arXiv:2604.17084 (2026).

Progress summary

Refreshed
Claimed solved

Bauschke and Gao’s conjecture is settled in one special case, while a reader-written argument now claims a complete proof for all parameters but has not been independently verified.

Bauschke and Gao (2026) conjectured that, for every α>2\alpha>2 and every nn, the coefficient row from k=1k=1 through k=nk=n is unimodal. Their paper supplied numerical motivation and no general proof.

Known results

  • For α=4\alpha=4, Proposition 5.3 proves unimodality of the full coefficient row, including the small cases.
  • The source explains why k=0k=0 is excluded: for some α>2\alpha>2, the constant coefficient disrupts unimodality; more generally, cn,1−cn,0<0c_{n,1}-c_{n,0}<0 when α≤1+5\alpha\leq1+\sqrt{5}.

Posted attempt

A reader-written argument claims a complete proof for every real α>2\alpha>2 and n≥1n\geq1, with mode at max⁡{1,⌊n/2⌋}\max\{1,\lfloor n/2\rfloor\}. It uses a normalized recurrence and a path-matching expansion, but the argument has not been independently verified.

Current status (as of August 2026): The α=4\alpha=4 case is proved, while the general conjecture has an unverified complete-proof claim and therefore remains unresolved.

Sources

Solutions 1

ProofThis solution needs a summarySee full solutionHide full solution

Unimodality of the Boţ–Nguyen coefficient rows

Bauschke and Gao's Conjecture 6.5 in arXiv:2604.17084v2 concerns the coefficients cn,kc_{n,k} defined by their Lemma 4.1, equations (44) and (45b). The conjecture excludes cn,0c_{n,0}. That exclusion matters: the source explains that including the constant coefficient can destroy unimodality. The case α=4\alpha=4 is already proved in the source.

Theorem. For every real α>2\alpha>2 and every integer n≥1n\geq1, the sequence

cn,1,cn,2,…,cn,nc_{n,1},c_{n,2},\ldots,c_{n,n}

is unimodal. More precisely, it is nondecreasing up to

μn=max⁡{1,⌊n/2⌋}\mu_n=\max\{1,\lfloor n/2\rfloor\}

and nonincreasing thereafter. The inequalities asserted here are weak; uniqueness of the mode is not needed.

Main idea. Normalize each row by its positive constant coefficient. The resulting recurrence has a nonnegative expansion indexed by matchings of a path. Grouping terms with all but the first marked edge fixed lets the decreasing edge weights produce sums of symmetric unimodal polynomials whose middle modes agree. Only the constant coefficient needs to be adjusted, so every coefficient in the conjectured range is unchanged.

1. A normalized recurrence

Put β=α/2>1\beta=\alpha/2>1, and let

pn(t)=∑k=0ncn,ktk.p_n(t)=\sum_{k=0}^n c_{n,k}t^k.

The source recurrence is

p0(t)=1,p1(t)=β+1+βt2β+1,\begin{aligned} p_0(t)&=1,\\ p_1(t)&=\frac{\beta+1+\beta t}{2\beta+1}, \end{aligned}

and, for n≥2n\geq2,

pn(t)=n+βn+2β(1+t)pn−1(t)−nn+2βtpn−2(t).(1)\begin{aligned} p_n(t)&=\frac{n+\beta}{n+2\beta}(1+t)p_{n-1}(t)\\ &\quad-\frac{n}{n+2\beta}t p_{n-2}(t). \end{aligned} \tag{1}

For the rising factorial (x)n=x(x+1)⋯(x+n−1)(x)_n=x(x+1)\cdots(x+n-1), with (x)0=1(x)_0=1, the constant term of (1) gives

cn,0=(β+1)n(2β+1)n>0.c_{n,0}=\frac{(\beta+1)_n}{(2\beta+1)_n}>0.

Consequently, the normalized polynomials Hn=pn/cn,0H_n=p_n/c_{n,0} satisfy

H0=1,H1=1+at,a=ββ+1∈(1/2,1),\begin{gathered} H_0=1,\qquad H_1=1+a t,\\ a=\frac{\beta}{\beta+1}\in(1/2,1), \end{gathered} Hn=(1+t)Hn−1−(1−rn)tHn−2,rn=β(β−1)(n+β)(n+β−1).(2)\begin{aligned} H_n&=(1+t)H_{n-1}\\ &\quad-(1-r_n)tH_{n-2},\\[3pt] r_n&=\frac{\beta(\beta-1)}{(n+\beta)(n+\beta-1)}. \end{aligned} \tag{2}

Indeed,

1−rn=n(n+2β−1)(n+β)(n+β−1).1-r_n=\frac{n(n+2\beta-1)}{(n+\beta)(n+\beta-1)}.

The quantities rnr_n are positive and strictly decreasing in nn. We will prove a recurrence lemma that only requires

12≤a≤1,r2≥r3≥⋯≥rn≥0.(3)\begin{gathered} \frac12\leq a\leq1,\\ r_2\geq r_3\geq\cdots\geq r_n\geq0. \end{gathered} \tag{3}

No upper bound on the rir_i is required for that lemma.

2. Symmetric unimodal polynomials

Write

SL(t)=1+t+⋯+tL(L≥0).S_L(t)=1+t+\cdots+t^L\qquad(L\geq0).

A polynomial is called symmetric of order NN if its coefficient array, padded by zeros to indices 0,…,N0,\ldots,N, is symmetric under k↦N−kk\mapsto N-k. Its actual degree may be less than NN.

Nonnegative symmetric unimodal polynomials of order NN are precisely the nonnegative linear combinations of

tjSN−2j(t),0≤j≤⌊N/2⌋.(4)\begin{gathered} t^jS_{N-2j}(t),\\ 0\leq j\leq\lfloor N/2\rfloor. \end{gathered} \tag{4}

To see this, the coefficient multiplying the term with index jj is the difference between the coefficients at jj and j−1j-1, with the coefficient at −1-1 set to zero. The converse follows directly from the interval of ones in each summand.

Products of such polynomials are again nonnegative, symmetric and unimodal, with their orders added. One elementary proof uses

Sb(t)Sc(t)=∑j=0min⁡(b,c)tjSb+c−2j(t),(5)\begin{aligned} &S_b(t)S_c(t)\\ &\quad=\sum_{j=0}^{\min(b,c)}t^jS_{b+c-2j}(t), \end{aligned} \tag{5}

whose coefficients on both sides count the pairs of integers in [0,b]×[0,c][0,b]\times[0,c] with prescribed sum, and then applies (4).

In particular, a nonnegative symmetric unimodal polynomial of order NN, or of order N−1N-1, has a mode at ⌊N/2⌋\lfloor N/2\rfloor after padding to 0,…,N0,\ldots,N. For even NN, an order-N−1N-1 polynomial has equal middle coefficients at N/2−1,N/2N/2-1,N/2. For odd NN, an order-NN polynomial has equal middle coefficients at (N−1)/2,(N+1)/2(N-1)/2,(N+1)/2. Thus nonnegative sums of these two types are unimodal with the same specified mode.

3. A prefix-sum identity

For integers L≥0L\geq0 and 0≤h≤L0\leq h\leq L, set

UL,h(t)=∑ℓ=0h(Sℓ(t)+1)SL−ℓ(t).\begin{aligned} &U_{L,h}(t)\\ &\quad=\sum_{\ell=0}^h\bigl(S_\ell(t)+1\bigr)S_{L-\ell}(t). \end{aligned}

The following identity will supply the needed positivity:

UL,h(t)=(h+2)SL(t)+∑j=1min⁡(h,L−h)tjSL−2j(t)+∑ℓ=1hSℓ−1(t)SL−ℓ(t).(6)\begin{aligned} U_{L,h}(t)&=(h+2)S_L(t)\\ &\quad+\sum_{j=1}^{\min(h,L-h)}t^jS_{L-2j}(t)\\ &\quad+\sum_{\ell=1}^h S_{\ell-1}(t)S_{L-\ell}(t). \end{aligned} \tag{6}

Empty sums are zero. The first two terms on the right together are symmetric unimodal of order LL; the last sum is symmetric unimodal of order L−1L-1. When L=0L=0, the latter is zero.

Here is an algebraic verification valid for every L,hL,h. Subtract the last sum in (6) from the definition of UL,hU_{L,h}, and use Sℓ−Sℓ−1=tℓS_\ell-S_{\ell-1}=t^\ell. The result is

2SL+∑ℓ=1h(1+tℓ)SL−ℓ.(7)2S_L+\sum_{\ell=1}^h(1+t^\ell)S_{L-\ell}. \tag{7}

For 1≤ℓ≤L1\leq\ell\leq L, define

Eℓ=(1+tℓ)SL−ℓ−SL.E_\ell=(1+t^\ell)S_{L-\ell}-S_L.

If 2ℓ≤L2\ell\leq L, then Eℓ=tℓSL−2ℓE_\ell=t^\ell S_{L-2\ell}. If 2ℓ=L+12\ell=L+1, then Eℓ=0E_\ell=0. The remaining terms obey

Eℓ=−EL+1−ℓ.E_\ell=-E_{L+1-\ell}.

Pairing these terms in the prefix ℓ=1,…,h\ell=1,\ldots,h leaves exactly

∑ℓ=1hEℓ=∑j=1min⁡(h,L−h)tjSL−2j.\begin{aligned} &\sum_{\ell=1}^h E_\ell\\ &\quad=\sum_{j=1}^{\min(h,L-h)}t^jS_{L-2j}. \end{aligned}

Substitution into (7) proves (6), including h=0h=0, h=Lh=L, and both parities of LL.

If w0≥⋯≥wL≥0w_0\geq\cdots\geq w_L\geq0, ordinary finite summation by parts gives

∑ℓ=0Lwℓ(Sℓ+1)SL−ℓ=wLUL,L+∑h=0L−1(wh−wh+1)UL,h.(8)\begin{aligned} &\sum_{\ell=0}^L w_\ell(S_\ell+1)S_{L-\ell}\\ &\quad=w_L U_{L,L}\\ &\qquad+\sum_{h=0}^{L-1}(w_h-w_{h+1})U_{L,h}. \end{aligned} \tag{8}

All weights on the right are nonnegative, including the final boundary weight wLw_L. Equations (6) and (8) therefore decompose the left side into nonnegative symmetric unimodal parts of orders LL and L−1L-1.

4. The path-matching expansion

Fix n≥1n\geq1, and consider a path with vertices 1,…,n1,\ldots,n. Label its edge {i−1,i}\{i-1,i\} by ii, for 2≤i≤n2\leq i\leq n. A matching is a set M⊆{2,…,n}M\subseteq\{2,\ldots,n\} with no consecutive indices. Write

s=∣M∣,rM=∏i∈Mri,s=|M|,\qquad r_M=\prod_{i\in M}r_i,

where an empty product is one. After deleting the endpoints of its marked edges, let L0,L1,…,LsL_0,L_1,\ldots,L_s be the lengths of the remaining consecutive vertex segments, in order; zero lengths are allowed. Thus

L0+⋯+Ls=n−2s.L_0+\cdots+L_s=n-2s.

The recurrence (2), for arbitrary aa, has the expansion

Hn(t)=∑MrMts⋅(aSL0(t)+1−a)⋅∏j=1sSLj(t).(9)\begin{aligned} H_n(t)&=\sum_M r_M t^s\\ &\quad\cdot\bigl(aS_{L_0}(t)+1-a\bigr)\\ &\quad\cdot\prod_{j=1}^s S_{L_j}(t). \end{aligned} \tag{9}

For completeness, (9) can be derived from a tridiagonal determinant. Take diagonal entries 1+at,1+t,…,1+t1+at,1+t,\ldots,1+t, superdiagonal entries one, and subdiagonal entry (1−ri)t(1-r_i)t at edge ii. Expansion along the last row gives exactly (2), with the stated initial conditions. In the determinant's permutation expansion, every nontrivial cycle is an adjacent transposition, so the edge transpositions form a matching. Expand each factor −(1−ri)t-(1-r_i)t into −t+rit-t+r_i t, and fix the edges where the second term was chosen. These edges form the marked matching MM. The remaining determinant factors over the undeleted path segments, with all remaining edge parameters equal to zero. A standard segment of length LL has determinant SLS_L, since it satisfies SL=(1+t)SL−1−tSL−2S_L=(1+t)S_{L-1}-tS_{L-2}; the initial segment has determinant aSL+1−aaS_L+1-a. This proves (9).

Let QnQ_n denote the same recurrence with a=1a=1. Formula (9) shows that QnQ_n is a nonnegative sum of products

ts∏j=0sSLj,t^s\prod_{j=0}^sS_{L_j},

each symmetric unimodal of order nn. Hence QnQ_n has that property.

It remains to treat a=1/2a=1/2; denote the corresponding polynomial by KnK_n. In (9), the empty matching contributes (Sn+1)/2(S_n+1)/2. For a nonempty matching, multiplication by two changes the initial factor into SL0+1S_{L_0}+1.

Group the nonempty matchings by keeping all marked edges except the first fixed. There are two cases.

  • If there are no other marked edges, the first edge ranges over i=2,…,ni=2,\ldots,n. Put L=n−2L=n-2, ℓ=i−2\ell=i-2, and C(t)=1C(t)=1.
  • Otherwise, let jj be the first of the fixed later marked edges. The first edge ranges over i=2,…,j−2i=2,\ldots,j-2. Put L=j−4L=j-4, ℓ=i−2\ell=i-2, and let C(t)C(t) be the product of the SS-polynomials for the fixed segments after edge jj.

Groups with no admissible first edge are omitted. In either case, 0≤ℓ≤L0\leq\ell\leq L, the two initial segment lengths are ℓ,L−ℓ\ell,L-\ell, and the group's contribution to 2Kn2K_n is

rM′tsC(t)⋅∑ℓ=0Lrℓ+2(Sℓ+1)SL−ℓ,(10)\begin{aligned} &r_{M'}t^s C(t)\\ &\quad\cdot\sum_{\ell=0}^L r_{\ell+2}(S_\ell+1)S_{L-\ell}, \end{aligned} \tag{10}

where M′M' is the fixed set of later marked edges and s=∣M′∣+1s=|M'|+1. If DD is the degree of CC, then

L+D=n−2s.(11)L+D=n-2s. \tag{11}

The weights rℓ+2r_{\ell+2} in (10) are nonnegative and nonincreasing by (3). Apply (8) and then (6). The factor CC is a product of SS-polynomials, and multiplying also by tst^s adds 2s2s to the symmetry order. By (11), every term obtained is symmetric unimodal of order nn or n−1n-1, with a nonnegative scalar coefficient.

Adding the empty matching therefore proves that

2Kn(t)−1=Sn(t)+An(t)+Bn(t),(12)\begin{aligned} 2K_n(t)-1&=S_n(t)\\ &\quad+A_n(t)+B_n(t), \end{aligned} \tag{12}

where AnA_n and BnB_n are nonnegative symmetric unimodal polynomials of orders nn and n−1n-1, respectively. Either may be zero. Thus 2Kn−12K_n-1 is unimodal with a mode at ⌊n/2⌋\lfloor n/2\rfloor.

For n=1n=1, there are no nonempty matchings and (12) simply reads 2K1−1=S12K_1-1=S_1. For n=2n=2, the sole nonempty matching has L=0L=0, and (10) is 2r2t2r_2t, a nonnegative order-two term. These cases require no negative-index SS-polynomial.

Finally, the recurrence is linear in its initial conditions. For every a∈[1/2,1]a\in[1/2,1],

Hn=(2a−1)Qn+2(1−a)Kn.H_n=(2a-1)Q_n+2(1-a)K_n.

Consequently,

Hn−(1−a)=(2a−1)Qn+(1−a)(2Kn−1).(13)\begin{aligned} H_n-(1-a)&=(2a-1)Q_n\\ &\quad+(1-a)(2K_n-1). \end{aligned} \tag{13}

Both scalar weights are nonnegative, and both polynomials on the right are unimodal with a mode at ⌊n/2⌋\lfloor n/2\rfloor. This proves that the full coefficient sequence of Hn−(1−a)H_n-(1-a) has that mode.

5. Conclusion for the original coefficients

Equation (13) changes only the coefficient at index zero. Hence, for n≥2n\geq2, the coefficients of HnH_n at indices 1,…,n1,\ldots,n are nondecreasing up to ⌊n/2⌋\lfloor n/2\rfloor and nonincreasing afterward. For n=1n=1, the sequence has one term and is unimodal. Multiplication by the positive constant cn,0c_{n,0} preserves every inequality. Applying (3) to the parameters in (2) proves the theorem for all real α>2\alpha>2 and all n≥1n\geq1.