Equipartition of a Segment: Every Continuous Function on Subsegments of [0, 1] Vanishing on Degenerate Segments Takes Equal Values on Some Partition of [0, 1] into m SegmentsResearch Paper
Motivation
In fair division, a value may be assigned to each interval that a cut produces. If the value is continuous in the interval's endpoints, an equal-value partition is a natural target. Some valuations are nonnegative and additive, but a difference of competing measures or a signed utility need have neither property. Avvakumov and Karasev establish equal-value partitions even when the value changes sign and is not additive, provided it vanishes on degenerate intervals. Their result is used by Bogomolnaia and Moulin in work on guarantees for fair division. Avvakumov and Karasev, pp. 1–2; Bogomolnaia and Moulin.
The paper identifies earlier special cases: a nonnegative valuation follows from the Knaster–Kuratowski–Mazurkiewicz theorem; prime-power numbers of parts with a valuation of varying sign were known through degree-theoretic results; additive valuations admit an elementary argument. The remaining case permits a sign-changing, nonadditive valuation and an arbitrary positive number of parts. The authors distinguish this from an envy-free division problem with separate valuations for several players, which they state remains open outside prime-power cases. Avvakumov and Karasev, p. 2; Avvakumov and Karasev, Envy-free division using mapping degree; Meunier and Zerbib.
Setting
Write I for the set of all subsegments I=[a,b] of [0,1] with 0≤a≤b≤1. An interval may be degenerate, meaning a=b. A real-valued function f on I is continuous when its value varies continuously with the two endpoints. The required boundary condition is f([a,a])=0 for every a∈[0,1]. No sign, monotonicity, or additivity condition is imposed. A partition into m parts is a sequence 0=x0≤x1≤⋯≤xm=1, so its parts [xj−1,xj] may also be degenerate. Avvakumov and Karasev, Theorem 1.1, pp. 1–2.
The intermediate setting is a nice multivalued function. Its graph is a compact subset Z of I×(−1,1). The graph separates the bottom I×{−1} from the top I×{1} inside the cylinder C=I×[−1,1]: those two sides lie in different connected components of C∖Z. A representation of a nice graph is a continuous real function φ on C that is negative on the bottom, positive on the top, and whose zero set is the graph. The paper passes between graphs and representations in Lemma 2.2. Avvakumov and Karasev, Definition 2.1 and Lemma 2.2, p. 2.
For a prime p, a segment I=[a,b], and a value y, a zero partition is a sequence a=x0≤⋯≤xp=b such that φ([xj−1,xj],y)=0 for each j=1,…,p. The set Zp(φ) contains exactly those pairs (I,y) that admit a zero partition. This is the projection to the cylinder of the paper's simultaneous-solution set S for labeled partitions. Avvakumov and Karasev, equation (4.1) and Claim 4.4, pp. 7–9.
Formalization targets
The goal is Theorem 1.1: for every positive integer m and every continuous f satisfying the degenerate-interval condition, there are possibly repeated cut points with a common value:
0=x0≤x1≤⋯≤xm=1,f([x0,x1])=⋯=f([xm−1,xm]).
The theorem is already posed as the Open platform statement FairGuarantees.DivideChoose.theorem_of_4; this mission references that exact statement as its goal. It supplies the paper's supporting statements, including Lemma 2.2 on graph representations, Claim 4.4 on separation by projected zero partitions, and the prime-step Lemma 2.3. The latter asserts that, when φ([a,a],y)=y, there is a represented nice graph with the same degenerate identity whose zeros each yield a zero partition for φ. The milestone list includes the fixed-segment and continuous-family statements used to characterize Zp(φ), and equation (4.2) describing its degenerate boundary. Avvakumov and Karasev, Lemmas 2.2–2.3, pp. 2–3; Claims 4.3–4.4 and equation (4.2), pp. 8–9.
Significance
The theorem guarantees equal values for any positive number of consecutive parts under only endpoint continuity and the zero-on-degenerates condition. It covers valuations for which nonnegative or additive special-case results give no conclusion. It also provides the theorem invoked in the fair-guarantee application cited by the authors. Avvakumov and Karasev, pp. 1–2; Bogomolnaia and Moulin.
The mathematical result was proved in the paper; the platform's theorem statement remains Open and lacks a machine-checked proof there. Formalizing the paper's intermediate claims makes its topological assumptions and its treatment of degenerate parts explicit. The definitions of a separating graph, its representation, and zero partitions can also be reused for related partition problems. The present mission states the paper's prime-step targets, leaving the proofs and the underlying equivariant topology as work for solvers.
Difficulty
For a fixed signed valuation, continuity alone does not select one common value that works across every piece of an arbitrary partition. The nonnegative case has a simplex-covering argument, but that argument does not directly control values of both signs. In the paper's general case, the difficult claim is that the projected simultaneous-zero set separates an entire cylinder, even as the segment varies. Its source proof requires a mod-p statement for a space of labeled partitions and control of the resulting solution set. A further boundary issue arises on degenerate segments: a representation of the projected graph must satisfy ψ([a,a],y)=y, not merely have the right zeros. Avvakumov and Karasev, pp. 2, 7–10.
Formalization scope
Lean represents I=[a,b] by the endpoint pair (a,b) with 0≤a≤b≤1. The cylinder uses the closed value range [−1,1]; a nice graph lies in its open interior (−1,1). Functions are total Lean functions, but only their values on this cylinder enter the definitions and results. Continuity is required on the cylinder. Separation uses connected components of its complement, as in Definition 2.1, and explicitly requires the bottom and top to avoid the graph. This prevents a graph containing a boundary side from satisfying separation vacuously.
Every prime-step claim assumes p prime, the representation's strict boundary signs, and the degenerate identity. Zero partitions have exactly p consecutive parts whose endpoints are pinned to the original segment; repeated cuts are allowed. The goal reference likewise permits degenerate parts and assumes m≥1. The projected set is defined through existence of an ordinary ordered partition, independent of how the paper labels partitions in its configuration space. The needed existing Lean notions include ContinuousOn, IsCompact, connectedComponentIn, finite index types, and Nat.Prime. A complete proof development also needs results about the paper's labeled-partition complex, equivariant cycles and transverse maps. Contributions establishing that reusable topology, as well as proofs of the stated milestones, are within scope. Avvakumov and Karasev, §§2–4.
Selected references
Sergey Avvakumov and Roman Karasev, Equipartition of a Segment, Mathematics of Operations Research, 2023. arXiv:2009.09862v4.
Anna Bogomolnaia and Hervé Moulin, Guarantees in Fair Division: General or Monotone Preferences, 2019. arXiv:1911.10009v3.
Sergey Avvakumov and Roman Karasev, Envy-Free Division Using Mapping Degree, Mathematika 67, 2021. arXiv:1907.11183.
Frédéric Meunier and Shira Zerbib, Envy-Free Cake Division Without Assuming the Players Prefer Nonempty Pieces, Israel Journal of Mathematics 234, 2019. arXiv:1804.00449.
A Dynamic Contagion Risk Model with Recovery Features: Fluid Limit for the Final Fraction of Failed Nodes in Threshold Contagion with Growth on the Configuration ModelResearch Paper
Motivation
Financial networks, supply chains and contact networks share a failure mode: a node fails when enough of its counterparties have failed, and its failure in turn weakens its own counterparties. Threshold contagion on random graphs is the standard mathematical model of this. Amini, Cont and Minca (Resilience to contagion in financial networks, Mathematical Finance 2016) proved that on a configuration-model network the final fraction of failed nodes concentrates, as the network grows, around an explicit value determined by the smallest fixed point of a one-dimensional map.
In that model thresholds are fixed. In practice institutions recover while a crisis unfolds: capital is raised, earnings accrue, and a node that survives long enough can absorb more losses than at the start. Amini, Minca and Sulem (A Dynamic Contagion Risk Model with Recovery Features, Mathematics of Operations Research 2022; the mission follows the authors' preprint of February 2021) add this recovery feature: every surviving node's threshold grows linearly in calendar time, at a rate proportional to its connectivity. The main result, Theorem 5, identifies the scaling limit of the cascade with growth. It is the basis of the paper's analysis of equilibrium connectivity choice (§3), where agents trade off the rewards of links against contagion risk.
Setting
A network has n nodes. Node i has a connectivityλi∈N (in-degree and out-degree both equal λi) and an integer initial thresholdθi∈N. The empirical type distribution is μλ,θ(n)=#{i:θi=θ,λi=λ}/n and the mean connectivity is λˉ(n)=∑λ,θλμλ,θ(n). Two regularity conditions are assumed: (C1) μλ,θ(n)→μλ,θ for a probability distribution μ on N2; (C2) λˉ(n)→λˉ=∑λμλ,θ∈(0,∞).
Each node owns λi out-going and λi in-coming half-edges. The directed configuration model matches out-half-edges to in-half-edges by a uniformly random bijection; self-loops and multiple links are allowed.
The cascade: at time 0 the nodes with θi=0 fail. Interactions happen at times Tk(n) whose gaps are i.i.d. exponential with rate n. At each interaction one unrevealed out-link of a failed node, chosen uniformly, is revealed; its counterparty j gains one failed link and fails if its count of failed links reaches its current threshold θj+αλjTk(n)/λˉ(n), where 0<α<1 is the growth rate. Only surviving nodes grow. The cascade stops at kstop(n), when no unrevealed failed link remains; ∣Df(n)∣ is the number of failed nodes and Sλ,θ(n)(kstop(n)) the number of survivors of type (λ,θ).
The limit objects: Pλ,θ,ℓ(π) is the probability that the order statistics of ℓ i.i.d. uniform variables on [0,π] satisfy U(u)>(u−θ)/(αλ) for u=θ+1,…,ℓ; Bλ,0α=1 and, for θ≥1,
Jα(π)=∑λ,θλμλ,θBλ,θα(π)/λˉ; and the relaxed fixed pointπ∗=min{π∈[0,1]:Jα(π)≤π}.
Formalization targets
Goal: Theorem 5 (corrected)
With L=∑λ,θμλ,θBλ,θα(π∗): for every ε>0, with high probability ∣Df(n)∣/n≥L−ε. (i) If π∗=1 and μ0,θ=0 for θ≥1, then ∣Df(n)∣/n→1 in probability. (ii) If π∗<1 is a stable fixed point (Jα(π)<π just to the right of π∗), then
Lemma 16: a recursion in ℓ for Pλ,θ,ℓ(t), for 0<α<1 and θ<ℓ≤λ, valid for t≥t^λ,θ,ℓ=(ℓ−θ)λˉ/(αλ).
Lemma 15: for 0<α<1 and 0<θ≤λ, the piecewise linear system (DE) for the rescaled numbers of surviving nodes with ℓ failed links has a unique continuous solution on [0,λˉ), in closed form.
(24)–(25): the limiting failed fraction of a class is μλ,θBλ,θα(t/λˉ), and the limiting number of unrevealed failed links is λˉ(Jα(t/λˉ)−t/λˉ).
(27)–(28): for k=γλˉ(n)n+o(n), Tk(n)→γλˉ in probability.
Remark 6: the minimum defining π∗ is attained.
Significance
Theorem 5 gives the final size of a cascade with recovery as an explicit function of the type distribution and the growth rate. The paper uses it to define each agent's expected reward as a function of its connectivity and to study the equilibrium of the resulting connectivity game (§3). Without the theorem, the equilibrium analysis would rest on the heuristic of §2.3 alone.
On the formal side, the platform has no random graph model of configuration type, no network contagion process and no fluid limit of a random cascade. This mission introduces all three. The result is proved on paper (by Wormald's differential-equation method); none of it is machine-checked. The work is to formalize the proof, and the milestones are self-contained pieces of it: an order-statistic identity, an explicit ODE solution, an algebraic identification of the limit, a concentration bound for sums of exponentials, and a fixed-point existence result.
Difficulty
The obvious approach is to view the cascade as a Markov chain in the counts of surviving nodes by type and number of failed links, and to apply Wormald's theorem directly. Two things break. First, with growth a node's survival depends on the whole history of its failed-link count, not on its current value, so the state space grows with time: new states become relevant at each growth time t^λ,θ,ℓ, and the fluid limit has to be established inductively on successive time intervals. Second, the cascade stops at a random time, near the first zero of the limiting number of unrevealed failed links. Showing that the stopping time concentrates requires the stability condition on π∗ and a truncation of the type distribution in the degree, since μ may have infinite support. The fixed point need not be a crossing point of a smooth map, so the identification of the stopping time is delicate even when the drift is known.
Formalization scope
Network. Half-edges of node i are pairs (i,a) with a<λi; the matching is a uniformly random permutation of the half-edges. Each n has its own probability space: the matching, m=nλˉ(n) i.i.d. Exp(n) gaps, and m i.i.d. uniform variables on [0,1) that choose the revealed link (uniform among the available ones given the past). The probability of an event averages over permutations the product law of gaps and choices. "With high probability" and "in probability" are stated with explicit ε and limits of these probabilities.
Cascade. Computed step by step from the randomness and run for m steps; ∣Df(n)∣ and Sλ,θ(n)(kstop) are read off the final state. A node fails when its failed links reach (not exceed) its threshold, as in every formula of the paper. Initially failed nodes are those with θi=0.
Limit objects. Order statistics are encoded by counting (U(u)>s iff fewer than u variables are ≤s), and probabilities as Lebesgue measures of subsets of [0,π]ℓ. Sums over types are tsums with summability hypotheses; μ is nonnegative with total mass 1.
Corrections, all labelled. (a) Bλ,0α:=1; the printed formula gives a value below 1 for π>0, contradicting the paper's p. 9, and makes the theorem false. (b) Part (i) assumes no isolated nodes with a positive threshold. (c) Lemma 16 assumes t≥t^λ,θ,ℓ. (d) Lemma 15 is stated for 0<θ≤λ, with uniqueness among continuous solutions, and reads the printed λˉ−u as λˉ−t. (e) (24) ends in μλ,θBλ,θα(t/λˉ), not μλ,θ(1−Bλ,θα(t/λˉ)). "Stable fixed point" means Jα(π)<π on some interval (π∗,π∗+η).
No trivialization.∣Df(n)∣ and Sλ,θ(n) are computed from the random matching, the gaps and the revelation order, never from Bα or Jα; π∗ is defined from μ and α, not taken as a parameter.
Not included. Theorem 3 (the case α=0, cited from earlier work), the equilibrium results of §3, Wormald's theorem (Lemma 17, an external result), the in/out-degree extension (Theorem 18) and the displays (32)–(41) of the inductive argument.
Infrastructure. A complete proof needs Wormald's differential-equation method for Markov chains with bounded increments, concentration of the configuration model, and the elementary facts in the milestones. A general Wormald theorem and the configuration model are reusable well beyond this mission; contributions of either, as separate theorems, are welcome.
Selected references
H. Amini, A. Minca, A. Sulem, A dynamic contagion risk model with recovery features, Mathematics of Operations Research, 2022 (preprint February 2021). https://doi.org/10.1287/moor.2021.1174
H. Amini, R. Cont, A. Minca, Resilience to contagion in financial networks, Mathematical Finance, 2016. https://doi.org/10.1111/mafi.12051
B. Bollobás, A probabilistic proof of an asymptotic formula for the number of labelled regular graphs, European Journal of Combinatorics, 1980. https://doi.org/10.1016/S0195-6698(80)80030-8
Online Learning over a Finite Action Set with Limited Switching: An Algorithm for Prediction From Experts Whose Regret and Number of Switches Are Both O(√(T log n log(1/δ))) with Probability 1 − δResearch Paper
Motivation
In prediction from experts (PFE), an online learner chooses one of n actions in each of T rounds, and an adversary assigns each action a loss in [0,1]. Performance is measured by regret, the learner's total loss minus that of the best single action in hindsight. In many applications changing the action is itself costly: re-routing traffic, re-balancing a portfolio, re-deploying a model. The number of switches, rounds in which the action changes, is then a second resource to keep small.
It has long been known that both quantities can be kept at O(Tlogn)in expectation. Kalai and Vempala's multiplicative Follow the Perturbed Leader (Kalai–Vempala 2005) and Devroye, Lugosi and Neu's random-walk perturbation (Devroye et al. 2015) both achieve this. Devroye et al. asked for guarantees that hold with high probability. A bound in expectation says little about a single run: the switch count of these algorithms can have heavy upper tails. Altschuler and Talwar (arXiv:1803.01548, COLT 2018; Mathematics of Operations Research 2021) gave the first algorithms whose regret and switches are both O(Tlognlog(1/δ)) with probability 1−δ. They used this to characterise the optimal regret under a hard switching budget, both for PFE and for multi-armed bandits.
Setting
A loss sequence is (ℓt)t≥1 with ℓt:[n]→[0,1]. The adversary is oblivious: it fixes the whole sequence before the game, independently of the learner's randomness. A randomized online algorithmA draws internal randomness ω from a probability space (Ω,μ). At round t it plays an action it∈[n] that is a function of ω and of the past losses ℓ1,…,ℓt−1 only. This is full-information feedback: after round t the whole vector ℓt is revealed. Then
Follow the Perturbed Leader (FPL) with perturbations Pt(i) plays at round t the minimiser of ∑s=0t−1(ℓs(i)+Ps+1(i)), with ℓ0:=0. FPLε∗ perturbs only once, with P1(i)=−R(i)/ε, where R(1),…,R(n) are i.i.d. standard exponential.
Framework 1 turns an algorithm A into a batched algorithm B. B runs a copy of A with fresh randomness. The copy sees only the losses since its epoch began. The epoch ends just before the copy would exceed S′ switches, and a new copy starts there. E is the number of epochs used in T rounds. BFPLδ∗ is Framework 1 applied to FPLε∗ with
N=⌈log(2/δ)⌉,ε=21Nlogn/T,S′=135Tlogn/N.
Formalization targets
Goal: Theorem 1 (p. 8)
There is an absolute constant C>0 such that for all T,n≥1 and δ∈(0,21) some randomized online algorithm A satisfies, for every oblivious loss sequence,
The algorithm may depend on T, n and δ, as on the page, but not on the losses. The constant is left unspecified, so a sharper constant would not change the goal.
Lemma 3 (p. 10): if E[Switchesτ0(A)]≤cAτ0 with τ0=min(T,⌈8T/N⌉), then Framework 1 with S′=23cAT/N has P(E>N)≤e−N and P(SwitchesT>N(S′+1))≤e−N.
Corollary 6 for BFPL∗ (p. 14): RegretT≤SwitchesT+E+∑e<EmaxiRe(i)/ε.
Lemma 18 (p. 28): E[etX]≤nt/(1−t) for the maximum X of n standard exponentials.
Lemma 9 (p. 14): P(∑e=1NmaxiRe(i)>6Nlogn)≤e−N for N,n≥2.
Theorem 8 (p. 13): BFPLδ∗ itself satisfies the bound of Theorem 1 for n≥2.
Significance
The result separates two notions of "optimal". Expected-value bounds on switching do not control a single run, and the paper shows (Appendix B) that naive restarting does not concentrate. Theorem 1 is the black box behind the paper's later results: the Θ~ characterisation of minimax regret under a switching budget S, with a phase transition at S≈Tlogn, and the switching-cost corollary obtained by mini-batching. Framework 1 is generic: any algorithm with an expected switching bound of order τ gains a high-probability one.
The paper result is proved, but none of it has been machine-checked. A formalization yields a verified online-learning protocol with oblivious adversaries, a restart framework with fresh randomness, and verified tail bounds for maxima of exponentials. Formalizing it also repairs four printed slips, listed under Formalization scope.
Difficulty
The obvious approach runs A about log(1/δ) times in fixed-length blocks and applies Markov's inequality in each. It fails: a single block can use far too many switches with probability that is not small enough. Framework 1 therefore uses epochs of random length, ended by the switch budget. The epoch lengths are then dependent, and the final epoch is truncated at T. The concentration argument has to handle both. On the regret side, Lemma 5 is pointwise, but summing it over random epochs and then controlling the sum of perturbation maxima needs a sub-exponential tail bound for maxima of exponentials.
Formalization scope
Rounds are 1-based, and losses are ℓ:N→Finn→R with values in [0,1]. An online algorithm has type Ω→(t:N)→(Fin(t−1)→Finn→R)→Finn: at round t it sees only ω and ℓ1,…,ℓt−1. The adversary is a deterministic loss sequence, quantified after the algorithm, which is the page's own reduction (p. 5). High-probability statements bound the outer measure of the failure event by δ. Expectations are lower Lebesgue integrals of natural-valued switch counts. The O(⋅) of Theorems 1 and 8 is an absolute constant C>0 chosen before T, n and δ. Argmins break ties towards the lowest index. Framework 1 takes randomness ω:N→Ω under the infinite product measure, one coordinate per epoch.
Three formalizations of the goal would be trivial, and all are excluded. An algorithm that could see ℓt or the future would simply play the best action in hindsight. A constant C chosen after T could be C=T. A loss sequence chosen before the algorithm would let the algorithm depend on it.
The printed slips and how they are handled:
FPL* sign. The paper perturbs by +R/ε. With one-sided exponentials that variant violates Lemma 7, so the perturbation is −R/ε.
Lemma 5 and Corollary 6. These need M(SwitchesT+1) and +E; the printed forms fail at T=1.
Lemma 3. It takes an integer N≥2 with δ/2=e−N, it assumes the switching hypothesis only at the horizon the proof uses, and it counts one boundary switch per epoch.
BFPLδ∗. It uses N=⌈log(2/δ)⌉.
Measurability. Lemma 3 assumes that A's actions are measurable in ω.
The development needs:
product and infinite-product measures of exponential laws (in Mathlib);
Markov's inequality and Azuma–Hoeffding for submartingales with bounded differences;
the memoryless property of the exponential law;
Chernoff bounds.
The online-learning layer (protocol, regret, switches, FPL, the restart framework) is reusable for other switching-cost results in the paper and elsewhere. Proofs of any milestone are welcome, as are alternative routes to Theorem 1, such as the random-walk perturbation of §2.5.
Selected references
J. M. Altschuler and K. Talwar, Online learning over a finite action set with limited switching, COLT 2018; Mathematics of Operations Research 46(1), 2021. arXiv:1803.01548v2
A. Kalai and S. Vempala, Efficient algorithms for online decision problems, Journal of Computer and System Sciences 71(3), 2005. doi:10.1016/j.jcss.2004.10.016
L. Devroye, G. Lugosi and G. Neu, Random-walk perturbations for online combinatorial optimization, IEEE Transactions on Information Theory 61(7):4099–4106, 2015. doi:10.1109/TIT.2015.2428253
N. Cesa-Bianchi and G. Lugosi, Prediction, Learning, and Games, Cambridge University Press, 2006. doi:10.1017/CBO9780511546921
Learning Markov Models Via Low-Rank Optimization: KL Error Bound for the Rank-Constrained Maximum Likelihood Estimator of a Low-Rank Markov ChainResearch Paper
Motivation
Many systems studied in operations research and reinforcement learning are modelled as finite Markov chains with a large number of states p: the stops of a taxi in a city, a random walk on a network, the conformations of a molecule. Estimating the p×p transition matrix from one observed trajectory needs, in general, a trajectory long compared with the p2 unknown entries. In many such systems the transition matrix has, at least approximately, a small rankr≪p: the observed chain is driven by a low-dimensional latent process, as in state aggregation, where P(st+1∣st)=∑zP(zt=z∣st)P(st+1∣zt=z) (Zhu, Li, Wang & Zhang, arXiv:1907.00113, §1.1). Zhu, Li, Wang and Zhang show that a low-rank transition matrix can be estimated from a trajectory whose length is proportional to p (up to logarithmic factors), and they give two maximum-likelihood estimators with explicit rates. This mission formalizes the rate of the second, rank-constrained, estimator in Kullback–Leibler divergence.
Setting
Let X0,X1,…,Xn be a Markov chain on the states [p]={s1,…,sp} with transition matrixP (Pij≥0, ∑jPij=1, Pij=P(X1=sj∣X0=si)) and stationary distributionπ (π⊤P=π⊤). The chain is primitive (some power of P has all entries positive), and X0 is drawn from π. Write πmin=minjπj and πmax=maxjπj.
The right spectral gap is 1−ρ+, where, with ⟨h1,h2⟩π=∑jh1jh2jπj and the adjoint P∗=Diag(π)−1P⊤Diag(π),
ρ+=⟨h,1⟩π=0,⟨h,h⟩π=1sup21⟨(P+P∗)h,h⟩π.
The KL divergence of a transition matrix Q from P is DKL(P,Q)=∑i,jπiPijlog(Pij/Qij)1{Pij=0}. From the trajectory, nij counts the transitions si→sj, and the averaged negative log-likelihood is ℓn(Q)=−n1∑i,jnijlogQij. The empirical KL divergence is DKL(P,Q)=ℓn(Q)−ℓn(P).
Assumption 1 of the paper: rank(P)=r, and Pjk∈{0}∪[α/p,β/p] for constants α,β>0. The (α/p)-minoration (Remark 2) removes the 0: α/p≤Pjk≤β/p for all j,k. The rank-constrained MLE (5) is any
In Lean these objects are IsStochastic, IsStationary, IsPrimitive, piMin, piMax, trajProb, rhoPlus, klDiv, negLogLik, empKL, Assumption1, Minorated, RankFeasible and IsRankMLE in the namespace LowRankMarkov.RankMLE.
Formalization targets
Goal: Theorem 2, KL display (p. 10)
There are universal constants C1,C2>0 such that, for p≥2, under minoration with rankP=r and nπmax(1−ρ+)>max(20logp,logn), every rank-constrained MLE satisfies, for all ξ>0,
The constants are not fixed numerically, so a sharper analysis does not invalidate the goal.
Milestones
Lemma EC.1, both forms: DKL(u,v)≥2β2pα∥u−v∥22 for distributions with uj∈{0}∪[α/p,β/p], vj∈[α/p,β/p], and its matrix form ∥P−Q∥F2≤απminp2β2DKL(P,Q).
The rank remark (p. 14): rank(Pr−P)≤2r and ∥Pr−P∥∗≤(2r)1/2∥Pr−P∥F.
(EC.14): on the shell Q∈C(η), l≥1 and 2l−1η≤DKL(P,Q)≤2lη, ∥Q−P∥∗≤2β(2lηr/(pαπmin))1/2.
(EC.15), (EC.17): Bernstein-type tails p−(ξ−1) for Wn(1)=maxj∑i≤nPXi−1,j and Wn(2)=maxj∑i≤n1{Xi−1=sj}.
Lemma 3: with probability at least 1−C2exp(−ηπmaxrplogp/πmin), ∣DKL(P,Q)−DKL(P,Q)∣≤21DKL(P,Q)+C1πmaxβ2rplogp/(πminα3n) for every Q∈C with DKL(P,Q)≥η.
(EC.3): DKL(P,Pr)≤0. Together with DKL(P,Pr)=L(Pr)−L(P) it gives (EC.4), DKL(P,Pr)≤DKL(P,Pr)−DKL(P,Pr).
Significance
The bound says that, with high probability, the KL error of Pr is of order rplogp/n up to the factors πmax/πmin and β2/α3, plus a term independent of n controlled by ξ. A trajectory of length of order rplogp, rather than p2, therefore suffices: the low rank is what makes the estimation problem sample-efficient. Lemma EC.1 converts the KL bound into a Frobenius-norm bound, and Lemma 3 is a uniform law for the empirical KL process of a dependent sequence over the paper’s cone class.
The result is proved in the paper; nothing in it has been machine-checked. The mission asks for a formal proof of the KL display. Two points need care beyond transcription: the printed hypotheses do not suffice (see Formalization scope), and the paper's proof of Lemma 3 symmetrizes the dependent summands as if they were independent (EC.8), so a formal proof must supply a valid argument for that step.
Difficulty
Comparing ℓn(Pr) with ℓn(P) is immediate; the difficulty is to control DKL−DKL uniformly over a non-convex set of matrices, from one dependent trajectory. A union bound over a net of the feasible set loses polynomial factors in p, and i.i.d. empirical-process tools do not apply directly to a Markov chain. The printed Lemma 3 covers the cone class C, while the appendix proof restricts its shells to rank-r matrices. This leaves a gap for matrices in C of higher rank and for the claimed application to Pr when r>16. The paper combines a tail inequality for suprema of empirical processes of Markov chains (Adamczak 2008), a contraction step, the matrix Freedman inequality (Tropp 2011) for a matrix martingale, and a Bernstein inequality for Markov chains with a spectral gap (Jiang, Fan & Sun 2018) for the quadratic variations Wn(1),Wn(2). None of these is in Mathlib.
Formalization scope
States are Fin p, matrices Matrix (Fin p) (Fin p) ℝ, logarithms natural. A trajectory x=(X0,…,Xn) makes n transitions, and its probability is π(x0)∏k<nPxkxk+1; probabilities of events are finite sums, so no measure theory is needed. ρ+ is a real sSup; "ρ+<1" is not built into its definition. The nuclear and Frobenius norms come from the published definition CaiCandesShen.Convergence.Basic.
Hypotheses fixed relative to the page:
Minoration instead of Assumption 1(ii), in the goal and (EC.3). The printed Theorem 2 allows Pjk=0, and it is false then. For even p, let rows 1,…,p/2 of P equal 2/p on columns p/2+1,…,p and 0 elsewhere, and rows p/2+1,…,p equal 1/p. Then rankP=2, α=1, β=2, ρ+=0, π puts mass 1/3 on the first half and πmax/πmin=2. Every feasible Q has Qij≥1/p, so on each of the first p/2 rows Q puts mass at most 1/2 on the support of P, and DKL(P,Q)≥31log2. With ξ=0.9plogp the threshold is below 31log2 for all large n, so the probability is 1, while C2e−ξ<1 for large p. The broken step is (EC.3), which needs P feasible for (5); minoration with rankP=r makes it feasible. Lemma EC.1, (EC.14), (EC.15), (EC.17) and Lemma 3 keep the printed Assumption 1.
X0∼π, stated in Theorem 1 and used by the proof of Theorem 2 through Lemma 3.
Ergodic is read as primitive; πi>0 is stated (it is implied); p≥2 keeps logp>0 in denominators.
Lemma 3's set is the printed C(η), whose cone is ∥Q−P∥∗≤4×21/2∥Q−P∥F and which has no rank constraint. The appendix proof uses rank-restricted shells; for r>16 the rank remark alone does not place Pr in C. This is a proof gap, while the formal statement retains the page's quantification.
Lemma EC.1's matrix form and (EC.14) include the printed cone condition. (EC.14) is stated pointwise for members of Cl.
The statements rule out trivializing readings: C1,C2 are quantified before all data; the goal holds for every selection of minimizers, not for one fixed estimator; p≥2 and π>0 exclude the junk values x/0=0 and log1=0 that would collapse the thresholds; a sorry-free check shows the goal's hypotheses are jointly satisfiable (p=2, uniform P, n=28).
Not posed: Theorem 1 (nuclear-norm MLE), the Frobenius display of Theorem 2 (Lemma EC.1 turns the KL threshold into twice the printed Frobenius threshold, so the paper's one-line derivation gives C2e−ξ/2), Theorem 3 (lower bound), and the DC-programming results of §4. Useful reusable contributions: a Bernstein inequality for finite Markov chains with a spectral gap, the matrix Freedman inequality, and the reverse-Pinsker bound of Lemma EC.1.
Selected references
Z. Zhu, X. Li, M. Wang, A. R. Zhang, Learning Markov models via low-rank optimization, Operations Research, 2022; preprint arXiv:1907.00113v2, 2020. https://arxiv.org/abs/1907.00113
B. Jiang, J. Fan, Q. Sun, Bernstein's inequality for general Markov chains, arXiv:1805.10721, 2018. https://arxiv.org/abs/1805.10721
J. A. Tropp, Freedman's inequality for matrix martingales, Electronic Communications in Probability 16:262–270, 2011. https://arxiv.org/abs/1101.3039
R. Adamczak, A tail inequality for suprema of unbounded empirical processes with applications to Markov chains, Electronic Journal of Probability 13(34):1000–1034, 2008. https://arxiv.org/abs/0709.3110
J. A. Fill, Eigenvalue bounds on convergence to stationarity for nonreversible Markov chains, with an application to the exclusion process, The Annals of Applied Probability 1(1):62–87, 1991. https://doi.org/10.1214/aoap/1177005981
Skill Acquisition and Data Sales: Under Uncertain Skill Acquisition and Accurate Data, Higher Skill-Acquisition Costs Raise Data Clarity, Data Price and Price InformativenessResearch Paper
Motivation
Asset managers increasingly buy alternative data: satellite images, credit-card panels, web traffic. Raw data of this kind is worthless until someone has the skills to process it, and those skills are costly and uncertain to acquire. A fund can hire data scientists and still find that its models extract little. Huang, Xiong and Yang (Management Science, 2022) build an equilibrium model of a data seller who sells to investors who must first acquire such skills. They ask how the cost of acquiring skills feeds through the data market into the financial market: into how much the price reveals, the cost of capital and return volatility.
The model joins two lines of work. One is the classical noisy rational-expectations market of Grossman and Stiglitz (1980), in which prices partly reveal private information. The other is the literature on selling information started by Admati and Pfleiderer (1986, 1990), in which a monopolist seller degrades the information it sells to limit what the price leaks. The new feature is the "gold rush" of skill acquisition: an investor who pays the cost learns only afterwards whether her skill is high or low.
Setting
A risky asset pays θ~∼N(0,τθ−1) and there is noise demand u~∼N(0,τu−1). A data seller observes θ~+η~ with η~∼N(0,τη−1); the precision τη is the data accuracy. The seller chooses a clarityx∈(0,∞] and passes the data with personalized noise, so that a buyer of skill z receives a signal of precision xz about θ~+η~. The value x=∞ means the data are sold "as is".
A continuum of CARA investors with risk aversion γ is indexed by i∈[0,1]. Investor i can acquire data-processing skills at cost C(c,i), where c>0 is a cost parameter. The cost is continuous and strictly increasing in i, with C(c,0)=0 and C(c,i)→∞ as i→1, and it is strictly increasing in c. A skilled investor turns out high-skilled (ZH) or low-skilled (ZL) with probability 21 each, with 0<ZL<ZH. A mass λ of investors becomes skilled.
In the asset market, the price reveals θ~+η~ with precision α2τu, and α is pinned down by the market-clearing equation (14). Skilled investors buy the data at a Nash-bargained price q=βG, where G(x,α) is the certainty-equivalent trading gain (15) and β∈(0,1) is the seller's bargaining power. Investors acquire skills until the marginal investor is indifferent, (1−β)G(x,α)=C(c,λ) (19). The seller chooses x to maximize the profit λq. An equilibrium is a clarity x∗ with its (λ∗,α∗) such that no other feasible clarity gives more profit.
The market outcomes are:
the data priceq∗;
price informativeness1/Var(θ~∣p~);
the cost of capitalE(θ~−p~);
return volatilityσ(θ~−p~);
skilled investors' performance, the average certainty-equivalent advantage (21) of a skilled over an unskilled investor.
Formalization targets
Goal: Proposition 5
Fix 0<c1<c2. There is a threshold T such that, for every data accuracy τη≥T and every pair of equilibria at c1 and c2, both clarities are finite and
The paper writes these as signs of derivatives in c. The goal states them as comparisons between two cost levels, which is what the prose of the proposition asserts.
Milestones
Proposition 1: for given (λ,x), equation (14) has a unique positive root α.
Proof of Proposition 2: for each finite clarity x, a unique (λ,α) solves (14) and (19).
Proposition 2 (1): an equilibrium exists.
(C5): on feasible tuples the profit equals βλC(c,λ)/(1−β).
Proposition 3: in equilibrium, x∗<∞ if λ∗τη>α∗γ and x∗=∞ otherwise.
The clarity limit: as τη→∞, every solution x of (14) tends to 2αγ/(λ(ZH+ZL)).
The α limit: the equilibrium α∗ tends to τθ/τu.
(C18): the equilibrium λ∗ tends to a root of (C18).
Significance
The result. Proposition 5 is the paper's headline. When skills are uncertain, a rise in skill-acquisition costs makes prices more informative and lowers the cost of capital and return volatility. When skills are certain, these outcomes stay unchanged (Proposition 4). The mechanism is that the seller cannot tailor one clarity to buyers of unknown skill. To keep buyers interested as costs rise, she provides clearer data, and more information leaks into the price. The result separates the data market from a pure cost shock to information production. It gives empirical predictions for markets where skill acquisition is costlier, such as private or international equities.
Formalizing it. The proposition is proved in the paper by expansion in w=1/τη around w=0. A machine-checked version must make "sufficiently high τη" precise, with a threshold that holds uniformly over all equilibria. It must also control the exact finite-τη equilibrium rather than its limit. The mission also contributes the first noisy rational-expectations model with a strategic data seller on the platform: posterior variances, the equilibrium equation (14), the certainty-equivalent gain, and the seller's choice between finite clarity and selling "as is". No part of the paper has been formalized before.
Difficulty
Several comparisons follow at leading order in 1/τη. As τη→∞ the equilibrium mass converges to the root of (C18), which is decreasing in c, and the clarity, the data price and performance are, to first order, functions of that mass. The price-informativeness, cost-of-capital and volatility comparisons do not follow at leading order. At τη=∞ the ratio α∗ tends to τθ/τu whatever the cost, so these three outcomes are flat in c to leading order. Their signs come from the O(1/τη) correction to α∗, and the correction must be controlled jointly with the seller's optimal clarity. A newcomer's first idea is to set τη=∞ in the equilibrium conditions and differentiate. That loses exactly the term that carries the result. A further difficulty is that the seller's optimum is defined by comparison with every feasible clarity, including selling "as is", not by a first-order condition.
Formalization scope
All quantities are real numbers. The parameters satisfy τθ,τu,γ>0, 0<β<1, 0<ZL<ZH and τη>0. Two of these go beyond the paper's §2 domains, and both are needed:
β>0: with β=0 the profit is identically zero and every clarity is optimal;
ZL>0: this is the paper's own condition on p. 23.
The cost family is any admissible C(c,i): continuous and strictly increasing in i on [0,1), C(c,0)=0, C(c,i)→∞ as i→1−, and strictly increasing in c on (0,1). "Increasing" in the paper is read as strictly increasing. Clarity is either a finite x>0 or "as is". At "as is", (14) becomes γα=λτη and both conditional variances given a signal equal 1/(τθ+τη). These are the x→∞ limits of the finite-clarity formulas. The value x=0 is excluded from the choice set because it gives no positive root of (14).
The posterior variances (8), (10), (12), the gain (15), the price coefficients (C1)–(C4), the cost of capital (C12), volatility (C13) and performance (21) are taken as closed-form definitions. Their CARA-normal derivation (Bayes' rule, market clearing, certainty equivalents) is not formalized. Equilibria are defined relationally: a feasible tuple whose profit is at least that of every feasible tuple, with no choice functions.
The goal concludes, rather than assumes, that the equilibrium clarities are finite. It quantifies over every equilibrium at each cost level. It does not define the equilibrium through the τη=∞ economy or through the paper's expansions, which would replace the theorem by its heuristic.
A complete development needs real analysis only: the intermediate value theorem, monotonicity, compactness on (0,∞], and asymptotics in 1/τη. Contributions are welcome on any milestone, on an independent proof of Proposition 4's Δ=0 counterpart, and on a formal derivation of the closed forms from Gaussian conditioning and CARA utility, which would be reusable for other rational-expectations models.
Selected references
S. Huang, Y. Xiong, L. Yang, Skill Acquisition and Data Sales, Management Science 68(8), 2022; working paper version of July 2019. https://doi.org/10.1287/mnsc.2021.4117
S. J. Grossman, J. E. Stiglitz, On the Impossibility of Informationally Efficient Markets, American Economic Review 70(3), 1980. https://www.jstor.org/stable/1805228
A. R. Admati, P. Pfleiderer, Direct and Indirect Sale of Information, Econometrica 58(4), 1990. https://doi.org/10.2307/2938330
N. Gârleanu, L. H. Pedersen, Efficiently Inefficient Markets for Assets and Asset Management, Journal of Finance 73(4), 2018. https://doi.org/10.1111/jofi.12696
Managing Self-Replicating Innovative Goods: Characterization of the Optimal Production and Sales Policies under Strong ReplicabilityResearch Paper
Why self-replicating goods
Some goods can only be made by copies of themselves. RepRap and Prusa 3D printers print the parts of new printers, seed stock is grown from seed, and livestock herds breed. For such a good the inventory is also the production capacity: a unit sold today is a unit that will not replicate tomorrow. A firm launching such a good has to decide, at every instant, how much of its stock to use for replication and how much demand to serve. Demand for an innovative good is not exogenous either. It diffuses through the market by word-of-mouth, so sales held back today also slow future demand.
B. Hu and Z. Sun (Management Science) model this as a deterministic optimal control problem. Their demand side is the Bass diffusion model (Bass, 1969) in the supply-constrained form of Ho, Savin and Terwiesch (2002) and Kumar and Swaminathan (2003), where cumulative sales rather than cumulative demand drive imitation. Hu and Sun add diminishing word-of-mouth. Their main structural result, Theorem 1, says that under a condition they call Strong Replicability the optimal policy switches only once and can be found by a one-dimensional search.
The model
Time runs over [0,T], where T is the life cycle. The firm controls a production ratep(t) and a sales rates(t). The state consists of the cumulative demandD(t), the adjusted cumulative salesS(t) and the inventoryI(t). The demand rate is
d(t)=[α+mβS(t)][m−D(t)],
where m is the market size, α the coefficient of innovation and β the coefficient of imitation. Problem (1) of the paper is
Here π is the price, c the unit replication cost, h the holding cost rate, r the replication rate, ρ≥0 the discount rate and w≥0 the rate at which word-of-mouth decays. Unmet demand is lost. The paper assumes π>c and r>ρ. Production is on when p=rI and off when p=0. Sales are on when s=d and off when s=0.
Strong Replicability is condition (2),
(r−ρ)π>cr+h.
It says that letting a unit replicate for a short time and selling the offspring later is worth more than selling the unit now.
Two levels of the launch inventory organise the answer. I is the level at which keeping production and sales always on depletes the inventory exactly at T. I is the level at which keeping production off and sales on depletes it exactly at T.
Formalization targets
Goal: Theorem 1
Assume Strong Replicability. Then:
(a) if IL<I, the policy that keeps production always on and switches sales on once, at the time that depletes the inventory exactly at T, is optimal;
(b) if I<IL<I, the policy that keeps sales always on and switches production on once, at the time that depletes the inventory exactly at T, is optimal;
(c) if IL>I, the policy that keeps sales always on and production off is optimal, and it leaves inventory at T.
In each case, optimal means that the profit is at least
∫0T[πs~−cp~−hI~]e−ρtdt
for every admissible trajectory (D~,S~,I~,p~,s~), not only for policies of the stated shape. The switching times are not given in closed form, and the theorem asserts only that they exist.
Milestones
The thresholds are ordered, I<I (preamble of Theorem 1).
Proposition 1. Once production is switched on, an optimal policy keeps it on until T or inventory depletion.
Proposition 2, first half. Under Strong Replicability, every optimal trajectory has I(t)>0 for all t<T, whatever T is.
Significance
Theorem 1 turns an infinite-dimensional control problem into a choice of one switching time. Its managerial content is that a firm with a replicating good should, when stock is scarce, run replication at full capacity and hold back sales, losing some early customers, as long as the stock built up will eventually be sold. Proposition 2 says that under Strong Replicability the firm never sells out before the end of the life cycle. The paper's numerical study (§4) and its extensions (§5: backlogged demand, liquidity, stochastic diffusion) are built on this characterization.
The result has a published informal proof (pp. 36–37 of the manuscript), which applies necessary conditions from Pontryagin's maximum principle and then a perturbation and exchange argument. None of it is machine-checked. A complete formal proof would establish global optimality, against all integrable controls, for a nonlinear controlled ODE with mixed state–control constraints. That is a sufficiency statement, which the necessary conditions alone do not give.
Difficulty
The maximum principle gives necessary conditions only, and Problem (1) is not concave: the dynamics of D are bilinear in S and D, and the constraint p≤rI couples control and state. So the obvious route, solving the Pontryagin system and reading off the optimum, does not prove optimality. The paper's own argument is an exchange argument. It postpones sales or shifts production forward and claims each change improves profit, and it applies the change recursively ("by making the same argument recursively"). Making it rigorous requires comparison results for the Bass system under perturbed sales, which are delicate. Postponing sales lowers the adjusted sales S, and with it all future demand. The existence of the switching times also needs a continuity argument on the depletion time as a function of the switch.
Formalization scope
The Lean development lives in the namespace SelfReplic.Strong and uses real numbers throughout.
Parameters form a structure Params with fields m, price, c, h, r, alpha, beta, rho, w, T, IL. The predicate Params.Standing collects the page's assumptions IL>0, π>c, r>ρ≥0, w≥0. It adds the sign conditions m>0, T>0, α>0, β≥0, c≥0, h≥0, which the paper uses but never states. Params.StrongRepl is condition (2).
Admissible trajectories (IsAdmissible) have integrable controls p,s and continuous states D,S,I on [0,T]. The three ODEs hold in integral form at every t∈[0,T], the constraints 0≤p≤rI, 0≤s≤d hold at every t∈[0,T], and I(T)≥0. The integral form is used because bang-bang controls are discontinuous. Continuity and integrability prevent any integral from silently taking Lean's default value 0. Nonnegativity of I follows from 0≤p≤rI, as on p. 8.
Optimality (IsOptimal) compares the discounted profit with every admissible trajectory of the same parameters.
Policy shapes are pointwise statements on [0,T], and a switch instant counts as "on". The thresholds I,I enter as real numbers with their defining property (IsLowLevel, IsHighLevel), evaluated at the same parameters with IL replaced.
Proposition 1 is read up to null sets ("production has run before t0", conclusion almost everywhere). It carries the added hypothesis ρc+h>0, which the paper's proof uses. Without it (ρ=h=0), pausing production is costless and the statement fails.
A formalization that states Theorem 1 inside its own policy class (the best policy of shape (a) is the best of shape (a)) would be empty. Here the competitor in every case ranges over all admissible trajectories.
Not formalized: Theorem 2 (Weak Replicability), whose closed form (3) for the production interval is inexact when ρ>0; Corollaries 1–2; the second half of Proposition 2, which fails when IL≥m; the unnumbered bang-bang claim of p. 9; Propositions 3–4 (§5).
Contributions that would be reusable beyond this mission include existence and uniqueness for the controlled Bass system in integral form, a comparison lemma for D and S under pointwise-ordered sales, and a sufficiency (verification) theorem for bang-bang controls of problems that are linear in the control.
Selected references
B. Hu, Z. Sun, Managing self-replicating innovative goods, Management Science (doi:10.1287/mnsc.2020.3936); authors' manuscript of 2020-10-18. https://doi.org/10.1287/mnsc.2020.3936
T.-H. Ho, S. Savin, C. Terwiesch, Managing demand and sales dynamics in new product diffusion under supply constraint, Management Science 48(2), 187–206, 2002. https://doi.org/10.1287/mnsc.48.2.187.255
Supplier Encroachment in a Non-Exclusive Reselling Channel: The Encroaching Supplier May Sell Only Directly at a Positive Direct Cost, or Through Both Channels at a Negative OneResearch Paper
Why encroachment through a shared reseller matters
Manufacturers that sell through retailers increasingly also sell directly to consumers, through their own web stores or outlets. The operations literature calls this supplier encroachment. The classical analysis (Arya, Mittendorf and Sappington, 2007) assumes an exclusive reseller: the reseller carries only the encroaching supplier's product. In that setting a supplier whose direct channel costs more per unit than the reseller's still sells through the reseller, and a supplier whose direct channel is cheaper never does.
Real resellers are rarely exclusive. A retailer that carries one brand usually carries competing brands too, and its order from one supplier depends on its order from the others. Hotkar and Gilbert (Management Science, 2021) study what changes when the reseller is non-exclusive. Their headline finding, Proposition 5, is that both conclusions of the exclusive model can reverse.
Setting
Two suppliers, i∈{1,2}, sell partially substitutable products through one reseller. Supplier 1 can also sell directly. The reseller orders q1 and q2; supplier 1 sells y units directly. The substitutability between the two suppliers' products is θ12∈(0,1), written θ below. Supplier 1's direct product is a perfect substitute for her product at the reseller. Prices clear the market linearly:
Production and reselling cost nothing. A direct sale costs supplier 1 an extra c per unit, and c may be negative: the direct channel can be more efficient than the reseller. With wholesale pricesw1,w2≥0, the reseller earns q1(p1−w1)+q2(p2−w2), supplier 1 earns y(pd−c)+q1w1, and supplier 2 earns q2w2.
The game has three stages:
the suppliers set w1 and w2 simultaneously, and neither observes the other's price;
the reseller, who sees both prices, orders q1 and q2;
supplier 1 observes q1, but not q2 or w2, and chooses y.
Supplier 1 holds a beliefw2′ about supplier 2's price. In stage III her quantity y and her conjecture q2′ of the reseller's order from supplier 2 are mutual best responses computed with w2′ (Lemma 1). In an equilibrium (the paper's (4.7)–(4.16)), beliefs equal the actual prices. The reseller best-responds at every price pair, anticipating supplier 1's reaction. Neither supplier gains from a unilateral price change. After a deviation by supplier 2, supplier 1's belief stays where it was. This is not a subgame-perfect equilibrium, because supplier 1 cannot observe w2.
The paper restricts the direct cost to c0(θ)<c<c3(θ) (p. 9). It identifies thresholds c0≤c1≤c2≤c3 (Table 2) that separate three regimes:
Direct Only (c≤c1) and Direct Only Boundary (c1≤c≤c2), in which the reseller orders nothing from supplier 1;
Reseller & Direct (c>c2), in which supplier 1 uses both channels.
The sign of c2 is that of
Num1(θ)=32−θ(80+θ(32+θ(14+θ(8−(4−θ)θ)))),
which has a unique root θ^≈0.3443 in (0,1).
Formalization targets
Goal: Proposition 5 (p. 22)
∃θ^∈(0,1),Num1(θ^)=0,such that
for 0<θ<θ^: c2(θ)>0, and for every c∈(0,c2(θ)] an equilibrium exists and every equilibrium has q1∗=0<y∗. Supplier 1 sells only directly, although her direct channel is less efficient;
for θ^<θ<1: c2(θ)<0. For θ^<θ≤9/10 and every c∈(c2(θ),0) there is an equilibrium with q1∗>0 and y∗>0. Supplier 1 sells through both channels, although her direct channel is more efficient.
Milestones
Lemma 1: the stage-III pair (yIII,q2′), which is unique.
Lemma 2: the reseller's conditionally optimal q2III.
Corollary 1: q2III=q2′ when w2=w2′.
Lemma 4: the threshold order c0≤c1≤c2≤c3.
Lemma 4: for c≤c2, the Direct Only and Direct Only Boundary prices are equilibria, and every equilibrium has q1∗=0<y∗.
Lemma 4: the Reseller & Direct equilibrium for c2<c<0 (stated for θ≤9/10).
The sign facts about Num1 used in the proof of Proposition 5.
Significance
With an exclusive reseller, efficiency alone decides channel use: a costlier direct channel is used alongside the reseller, and a cheaper one replaces it. Proposition 5 shows that this rule fails once the reseller also carries a competitor.
When the products are weakly substitutable, supplier 1 may abandon the reseller although it is the cheaper channel. Supplier 2 then prices just low enough to keep her out.
When they are strongly substitutable, supplier 1 keeps selling through the reseller although her own channel is cheaper. The reseller's order shields her from supplier 2.
The result explains why observed channel structures need not reflect relative channel costs. It also shows which assumption of the exclusive model drives its conclusions.
As far as the mission's searches found, no part of this model has been machine-checked. The mission formalizes a game with unobservable prices and belief-based stage-III conjectures. Such games occur throughout the supply-chain contracting literature, and the equilibrium predicate here can be reused for them. The paper's proofs are algebraic, but they verify deviation incentives across a piecewise-defined reseller response, with polynomials of degree eight in θ. A formal proof has to check every branch.
Difficulty
The reseller's best response is piecewise in q1: supplier 1 sells directly only while q1 is below a threshold, and her conjecture about q2 can hit zero. The paper analyzes the reseller's choice under a restriction (Lemmas 1–3) that rules out the conjecture q2′=0 for small q1, and argues that the excluded region cannot occur in equilibrium. Deviation prices, however, can push the reseller into that region. There the reseller's objective is not concave in q1, and supplier 1 may find a price cut that triggers a jump in the reseller's order. A proof of Lemma 4 cannot simply cite first-order conditions inside the restricted region: it has to rule out such jumps for every deviation price. Numerically this succeeds for θ≤9/10, but it fails for θ close to 1.
The uniqueness claim in part (a) is a statement about every strategy profile, including the reseller's off-path rules. It cannot be obtained by exhibiting one profile. The equilibrium prices themselves are not unique for c≤c2, so this claim concerns the channel use q1∗=0<y∗ only.
Formalization scope
Everything is real-valued and lives in NonExclEncroach.Channels. The conventions are:
the paper's case θ1d=1 (§4), hence θ2d=θ; §5 (θ1d=0.8) is out of scope;
the standing range 0<θ<1 and c0<c<c3 (p. 9), with c allowed to be negative;
prices and quantities in [0,∞);
best responses as argmaxes over [0,∞), with no tie-breaking rule;
a strategy profile made of the two prices, supplier 1's direct-quantity and conjecture rules as functions of the observed q1, and the reseller's order rule at every price pair. Supplier 1's belief is the equilibrium price of supplier 2 and does not move after supplier 2 deviates.
The goal quantifies over this equilibrium predicate. Stating Proposition 5 with the quantities read off Table 1, or as the sign of c2 alone, would reduce it to the sign of Num1, and is ruled out.
Three statements depart from the page, and each departure is disclosed in its Formalization Note:
Part (b) of the goal and the Direct Only/Boundary part of Lemma 4 are stated for θ≤9/10. The Reseller & Direct part of Lemma 4 is stated for θ≤9/10 and c2<c<0. A numerical solution of the game shows that the paper's equilibrium prices fail outside these ranges: for θ near 1, for θ≥0.85 with c>0.3, and just below c3 (where supplier 2 gains by raising her price until the reseller shuts down the direct channel). The paper also announces that the Direct Only and Direct Only Boundary prices are the unique equilibrium prices for c≤c2. That is false: a lower w2, with w1 at the reseller's threshold, can also be an equilibrium. Those prices are therefore stated as existence.
"There exist positive (negative) values of c" is strengthened to every c in (0,c2] (in (c2,0)).
Lemma 2 assumes that the upper branch of (4.19) is non-negative.
Not formalized: Lemma 3, whose threshold 2c−1 is correct only when w2=w2′; the appendix lemma eliminating equilibria with q2′=0; and Propositions 1–4 and 6–8.
Contributions welcome: proofs of the milestones, a general existence theorem for the reseller's order rule, and an analysis of the game for θ near 1.
A. Arya, B. Mittendorf and D. E. M. Sappington, The Bright Side of Supplier Encroachment, Marketing Science 26(5), 2007. https://doi.org/10.1287/mksc.1070.0280
D. G. Pearce, Rationalizable Strategic Behavior and the Problem of Perfection, Econometrica 52(4), 1984. https://doi.org/10.2307/1911197
Design and Pricing of Discretionary Service Lines: Under SPT Without Customer Heterogeneity, the Optimal Service Line Has Variety, No Cutoff and a Concave Price Falling at the TopResearch Paper
Motivation
Many services are discretionary: the provider chooses how long to spend with a customer, and a longer service is worth more to the customer. Medical consultations, legal and tax advice, repair and diagnostic work, and personal services such as hair styling all have this shape. A firm offering such services designs a service line, a menu of service durations with a price for each, and runs it on one shared resource. Longer services are more valuable, but they also slow the queue for everyone else, and customers dislike waiting.
L. Debo and C. Li, in Design and Pricing of Discretionary Service Lines (SSRN 3054741, 2017; later published in Management Science, doi:10.1287/mnsc.2020.3670), study this problem by combining a Mussa–Rosen style product-line model with a single-server queue. A central question is whether congestion alone, with all customers identical, can make variety optimal. Under first-in-first-out (FIFO) service the answer is no: with homogeneous customers the optimal line collapses to a single service. Under the shortest-processing-time (SPT) priority rule the paper shows that the optimal line does offer variety. This mission formalizes that result.
Setting
Customers arrive as a Poisson process of rate Λ>0. A customer who receives a service of duration t gets value V(t,θˉ)=θˉt−bt2 with b>0, and pays a waiting cost c per unit of sojourn time, waiting plus service. All customers have the same quality sensitivity θˉ: the paper's heterogeneity parameter Δ is 0. The standing assumption is 0<c<θˉ.
A service line is a shortest duration t, a longest duration tˉ, and a demand measureμ on [t,tˉ]⊆[0,∞). When demand has a density λ(t), customers choosing durations in [t,t+dt] arrive at rate λ(t)dt. Its moments are
Λk=∫τkμ(dτ),k=0,1,2,
so Λ0 is the throughput and Λ1 the load on the server. A line is feasible if Λ0≤Λ (customers may be turned away) and Λ1<1 (the queue is stable).
The server uses non-preemptive SPT: whenever it frees up, it starts the shortest waiting service. By Phipps and Van Voorhis (1956), a customer with duration t waits in queue for
W(t)=(1−∫ttτλ(τ)dτ)2Λ2/2
on average, which is eq. (7) of the paper. The expected queue length is Q=∫Wdμ. With Δ=0 the paper's mechanism-design problem (its Proposition 1) reduces to the firm's problem
maxΠ=(θˉΛ1−bΛ2)−c(Λ1+Q)over feasible service lines,(9)
and the price of duration t is
P(t)=θˉt−bt2−c(t+W(t)),
which leaves every customer zero surplus.
Formalization targets
Goal: Proposition 5
For every Λ>0, problem (9) has an optimal service line, and every optimal line has
service variety, t<tˉ;
no customer cutoff, Λ0=Λ;
a price P that is concave on [t,tˉ], with
dtdP(tˉ)<0.
There are also thresholds Λ,Λˉ>0, depending on (θˉ,b,c), such that
dtdP(t)<0 for Λ∈(0,Λ),dtdP(t)>0 for Λ∈(Λˉ,∞).
The thresholds are existential, so the goal does not commit to their values.
Milestones
Corollary 1. For every feasible line, Q≥21{ln(1−Λ1)}2. The bound does not depend on Λ0, and it is attained by λ(t)∝1/t3 on a line with t<tˉ.
Root of (A–35). The equation F(ω)=0 of the paper's Lemma A7 has a root ω∈(0,1).
Lemma A7. Every optimal line is given in closed form by a root ω of (A–35): Λ0=Λ, Λ1=1−ω, and explicit Λ2, t, tˉ.
Quadratic waiting time. At an optimum, W(t)=−α0+α2t2 on [t,tˉ] with α2>0 and 2α2tˉ=θˉ/c−1.
Proposition 5 split into its parts: (i)–(ii), (iii) at the top end, and (iii) at the bottom end.
Significance
The result separates two sources of variety. Classical product-line theory obtains variety from heterogeneous tastes. Proposition 5 obtains it from operations alone: under SPT, short services jump the queue, so a spread of durations lowers total waiting below what any single duration achieves (Corollary 1 makes this exact). This outweighs the loss in service value from moving away from a single duration. The price result is the operational counterpart. Customers who choose the longest services wait the longest, so the price falls at the top of the line, and in small markets it falls along the whole line. This contrasts with FIFO, where the price always increases in duration.
Formalizing these results adds the following. The paper's proofs are long first-order-condition calculations over a constrained calculus-of-variations problem. They take existence of an optimum for granted, and they argue the bottom-end sign through Taylor approximations in the limits 2b/(Λc)→0 and →∞. A machine-checked proof would establish existence, rigorously derive the optimality conditions over all demand measures, and turn the asymptotic sign argument into a proof. As far as is known, none of these results, nor the SPT waiting-time formula, has a machine-checked proof.
Difficulty
The obvious approach is to optimize pointwise in λ(t). This fails because W at duration t depends on the whole demand below t through the cumulative load, and on Λ2, so the objective is a nonlocal, nonconvex functional of the demand schedule. A second obvious idea is to compare with the best single service. That comparison shows only that some non-degenerate line beats one particular collapsed line. It neither produces an optimum nor rules out an optimum with an atom. Existence also needs an argument: the feasible set of measures is not compact in any naive sense, because the support can drift and the load constraint is strict. Finally, the sign of P′(t) depends on where the root of the transcendental equation (A–35) lies, and the paper controls this only asymptotically.
Formalization scope
Everything is real-valued. A service line is a Lean structure with endpoints lo, hi and a measure μ : Measure ℝ. Admissibility requires:
μ is finite and lives on [t,tˉ] with 0≤t≤tˉ;
μ charges every right-neighbourhood of t and every left-neighbourhood of tˉ, so the endpoints are the shortest and longest services actually offered;
Λ1<1.
Feasibility adds Λ0≤Λ. The waiting time is
W(t)=(1−L−(t))(1−L+(t))Λ2/2,
with L− and L+ the loads of durations <t and ≤t. This is eq. (7) for densities. At atoms it is Cobham's non-preemptive priority wait, with ties served FIFO, and for a single service it is the M/D/1 value Λ2/(2(1−Λ1)).
Point masses are allowed on purpose. If demand were restricted to densities, a collapsed line would be infeasible and t<tˉ would hold by construction. Letting the endpoints float freely, beyond the support, would trivialize it in the same way. Both encodings are therefore excluded.
The following conventions are committed:
the reduced problem (9) is used instead of the mechanism of Proposition 1;
Δ=0 is built in, with no parameter for θ;
V is uncapped;
the price derivatives at t and tˉ are one-sided, taken within [t,tˉ], because P has a kink at tˉ;
existence of an optimum is stated explicitly;
uniqueness of the optimum or of the root of (A–35) is not claimed;
the hypotheses b>0 and Λ>0 are added, since they are implicit in the paper.
A complete development needs: integration against finite measures on R; existence of optimal measures under moment and load constraints; first-order conditions for the nonlocal functional Q; and real analysis of F. The SPT waiting-time functional and Corollary 1's queue-length bound are reusable for other priority-queue design problems. Contributions are welcome on any milestone, in particular Corollary 1 and the existence part of Lemma A7.
T. E. Phipps Jr. and W. R. Van Voorhis, Machine repair as a priority waiting-line problem, Operations Research 4(1), 76–86, 1956. https://doi.org/10.1287/opre.4.1.76
A. Cobham, Priority assignment in waiting line problems, Journal of the Operations Research Society of America 2(1), 70–76, 1954. https://doi.org/10.1287/opre.2.1.70
Nonzero-Sum Stochastic Games and Mean-Field Games with Impulse Controls: A Regular Solution of the N-Player Quasi-Variational Inequalities Gives a Nash EquilibriumResearch Paper
Why this game matters
An impulse control changes a system's state at a chosen instant, rather than continuously adjusting its drift. When several decision makers can intervene in the same state process, an intervention changes both the acting player's future cost and every opponent's. This makes the equilibrium question different from separate single-player control problems. Basei, Cao and Guo formulate a multidimensional, N-player stochastic impulse game and give a verification theorem for Nash equilibria in their preprint. This mission formalizes that theorem, Theorem 2.12, and the path and QVI statements on which its proof depends. It is the cost-minimizing, N-player generalization of the two-player impulse game of Aïd et al., posed independently because the conventions and equilibrium inequalities differ; that two-player verification theorem is itself an open target on the platform.
Game and costs
There are N≥1 players and a state X in a nonempty open domain S⊆Rd. Between interventions the state follows an Itô diffusion with drift b, volatility σ, and an M-dimensional Brownian driver. Player i chooses a strategy(Ai,ξi): a closed action region Ai and a continuous impulse rule ξi:S→Zi. When the state enters Ai, that player shifts it from y to Γi(y,ξi(y)). The process stops at its exit time τS from S. The action regions of an admissible profile are pairwise disjoint, although consecutive interventions at one instant can still occur when an impulse lands in another player's region. These rules are given in Definitions 2.1, 2.2 and 2.7 of the source.
Every player minimizes cost. Player i pays a discounted running cost fi, an own-intervention cost ϕi, a cross cost ψi,j when player j=i intervenes, and terminal cost hi if the game exits. With τi,k and ξi,k denoting player i's intervention times and impulses, the expected cost is
An admissible profile satisfies the source's integrability, moment and nonaccumulation conditions. A Nash equilibrium is an admissible profile for which no player can lower their own cost through an admissible unilateral change Definitions 2.6–2.8.
Formalization targets
For candidate values Vi and minimizing impulses ξi, define the own intervention operatorMiVi(y)=Vi(Γi(y,ξi(y)))+ϕi(y,ξi(y)) and the cross intervention operatorHi,jVi(y)=Vi(Γj(y,ξj(y)))+ψi,j(y,ξj(y)). Put Di={y∈S:MiVi(y)−Vi(y)>0}, Ai∗={y∈S:MiVi(y)−Vi(y)=0}, and Wi=S∩⋂j=iDj. The QVI system requires a boundary limit Vi→hi, nonnegative own intervention gaps, zero cross gaps when an opponent acts, and the diffusion condition
The goal, Theorem 2.12, says that if every ξi uniquely minimizes the own intervention expression, the Vi solve the full QVI system and have the stated smoothness and boundary regularity, and φ∗=((Ai∗,ξi))i is admissible at x, then (with one added hypothesis, explained under Formalization scope)
φ∗ is a Nash equilibrium,Vi(x)=Ji(x;φ∗)for every i.
The milestone list follows the proof's path: the controlled-process representation and jumps in Lemma 2.3; the candidate opponents' waiting and intervention properties in (2.15); the smooth approximations Vi,m of (a)–(c), which the paper takes from Øksendal, proof of Thm. 10.4.1 and App. D; the three pointwise comparisons (2.19)–(2.21), a thin bridge from the QVI to the proof; the deviation bound in Step 1; and the value identity in Step 2. The localized Itô identity (2.18) is not a separate milestone: as printed, its expectation is not shown to be finite near the boundary of Wi, and with simultaneous interventions its jump terms must be read stage by stage, so a faithful statement needs choices the page does not make. It remains the natural intermediate target between (a)–(c) and Step 1.
What the result provides
The theorem turns a regular solution of a coupled, nonlocal QVI system into a certificate for an equilibrium and its expected costs. It supplies the bridge from analytic candidate value functions to the original stochastic game. The result is proved in the preprint; the Lean declarations here are open statements, so the remaining work is a machine-checked development of the path construction, Itô formula across interventions, approximation near action boundaries and limiting arguments. The stochastic basis and diffusion generator reused here are published definitions from the two-player development, while the N-player game and its minimizing QVI are new reusable interfaces.
Main obstacles
Applying a smooth Itô formula directly to Vi(Xt) fails because X jumps at interventions and Vi is only C1 across ∂Di. The source therefore uses smooth approximations and localizes before taking limits Theorem 2.12, proof. The analytic challenge is to control the generator of the approximants near a Lipschitz surface; the probabilistic challenge is to justify the expectation and the limits from the exact admissibility conditions. Treating the QVI as an inequality everywhere in S would impose a generator value where the second derivative need not exist.
Formalization scope
Lean represents the state as a Euclidean space, players as Fin N with N≥1, costs and values as real numbers, and times as [0,∞]. The published SolvesFrom predicate encodes the restarted SDE through a Brownian Itô integral; the filtration satisfies the usual conditions and supports an F-Brownian motion. The source names the Brownian natural filtration, so this is a disclosed generalization. An admissible realization is a family of diffusion pieces solving that same SDE. Membership in Φx asks for at least one realization, while cost identities and Nash comparisons cover every admissible realization. The expected cost includes an integrability requirement, and admissibility asks that the running-cost integrand be integrable in time along almost every path, which is how "provided that the right-hand side exists and is finite" (Definition 2.6) is read. This prevents Lean's integral default on a nonintegrable payoff from making the theorem empty.
S is open: the page calls it "a non-empty subset", but p. 9 says that Vi "is not defined in ∂S", which is openness, and the C2 condition (ii) and the generator in (2.13d) need an open domain. One hypothesis is added to Theorem 2.12, Step 1 and Step 2: along every admissible deviation, the time the state spends on ∂Di has Lebesgue measure zero almost surely. This is the occupation-time hypothesis of Øksendal's Theorem 10.4.1, the result whose proof the paper cites for (a)–(c); the passage from (2.17) to (2.18) uses it, because (c) controls LVi,m only off ∂Di. The page omits it, and the theorem fails without it when σ is degenerate: a path can run along ∂Di, where (2.13d) imposes nothing on fi. For a nondegenerate diffusion the hypothesis holds automatically. In (a)–(c), Wi is assumed open, as the page's C2(∩j=iDj) presupposes; in the proof it is open because S is open and the Aj∗ are closed.
The path construction uses entry into closed action regions, choosing the least index when multiple entry times agree; exit from S uses the next stage's exit time. Discount at infinite time is zero. The path norm in the moment condition is taken through τS, because a process continued beyond game exit could have an infinite global supremum and make Φx empty. Cross costs use player j's intervention time. Boundary data in (2.13a) are limits from within S; the generator equation and derivative bounds are required where the derivatives exist, off ∂Di. The common waiting region includes S, so it remains S for one player. Action regions are built from Vi and ξi, not free variables. The source's printed formulas also call scalar costs vector valued and state some jump identities without accounting for consecutive interventions; the definitions and Lemma 2.3 follow the operational readings described in the item notes. In Lemma 2.3 the avoidance of the action regions is stated at positive times, since the initial state may lie in an action region without an intervention at time 0. A statement that drops the QVI clauses, fixes one realization for the conclusion, restricts the deviations, or reads (2.8) as a supremum over all times (which makes Φx empty for bounded S) would trivialize the goal and is ruled out by the definitions.
The definitions of the N-player path, admissibility, cost and QVI are intended for reuse. Contributions that establish the localized Itô identity, smooth approximation and dominated-convergence steps are especially useful. Every proposed proof must preserve the full QVI system, all admissible unilateral deviations, and the source's standing assumptions.
Selected references
M. Basei, H. Cao and X. Guo, Nonzero-sum stochastic games and mean-field games with impulse controls, arXiv:1901.08085v4, 2020. https://arxiv.org/abs/1901.08085v4
R. Aïd et al., Nonzero-Sum Stochastic Differential Games with Impulse Controls: A Verification Theorem with Applications, Mathematics of Operations Research 45(1), 2020. https://doi.org/10.1287/moor.2019.0982
A General Analysis of Sequential Social Learning: Asymptotic Learning Holds for Every Utility Function if and only if Private Signals Are Totally UnboundedResearch Paper
Motivation
In sequential social learning, agents act one after another. Each agent observes the actions of its predecessors and a private signal about an unknown state of the world. The model explains herd behaviour, information cascades and the adoption of products, technologies and opinions when people imitate one another. The basic question is whether the population eventually acts as if it knew the state, or whether public information can lock agents into a wrong action.
The answer depends on the setting:
Bikhchandani, Hirshleifer and Welch (1992) and Banerjee (1992) showed that with finitely many signals agents eventually ignore their private information and may herd on a wrong action.
Smith and Sorensen (2000) characterized learning for binary states and actions: asymptotic learning holds if and only if private beliefs are unbounded, meaning a single signal can make an agent arbitrarily sure of either state.
Arieli and Mueller-Frank (preprint SSRN 2466903, 2019; Mathematics of Operations Research, 2021) treat general state spaces, compact action spaces and all bounded utilities that are continuous in the action. They characterize when learning holds for every such utility.
Setting
A social learning model(Ω,μ,A,S,F) consists of the following:
a separable complete metric state space Ω with a prior μ;
a compact metric action space A with at least two actions;
a standard Borel signal space S;
a Markov kernel F:Ω→Δ(S), ω↦Fω. Any two Fω,Fω′ are mutually absolutely continuous, and distinct states have distinct signal laws.
A state ω∼μ is drawn once. Agents t=0,1,2,… act in sequence. Agent t observes the history ht=(a0,…,at−1) and a private signal st∼Fω, drawn conditionally i.i.d. given ω. It then chooses an action by a mixed strategy σt(ht,st)∈Δ(A). All agents share a utility u from the class U of bounded measurable functions A×Ω→R that are continuous in the action. The profile σ, the prior and the kernel determine the play measurePσ on states, signals and actions.
A profile is a perfect Bayesian equilibrium (PBE) if no agent can raise its expected utility Eσ[u(at,ω)] by changing its own strategy. Write O(ω)=maxau(a,ω) and yt=u(at,ω).
Asymptotic learning holds for u if yt→O(ω)Pσ-almost surely in every PBE.
Adequate learning holds if asymptotic learning holds for every u∈U.
Let μF(E∣s) be the posterior probability of an event E⊆Ω after one signal s. The information structure has totally unbounded signals if the support of the law of μF(E∣s) contains both 0 and 1 for every event E with 0<μ(E)<1.
Formalization targets
Goal: Theorem 1
adequate learning holds⟺(F,μ) features totally unbounded signals.
The statement is an equivalence between two properties of the model. It has no constants.
Milestones
The milestones are the results the paper's proof uses, in attack order:
Lemma 3. A utility in U is jointly uniformly continuous on a compact set of measure at least 1−ε.
Lemma 4. A single agent observing two signals has an ε-optimal strategy that takes finitely many actions and is constant on the cells of a finite product partition.
Lemma 1. The vanishing value of private information: dt=Eσ[yt−y^t∣ht]→0 almost surely in every PBE, where y^t=maxaEσ[u(a,ω)∣ht].
Lemma 2. The belief at infinity μ∞(⋅∣h) is almost surely a degenerate prior: one more signal has zero value against it.
Lemma 5 (first part). Total unboundedness passes from μ to every ν≪μ with dν/dμ≤C.
Lemma 7. Conditioning P∈Δ(Δ(Ω)) on a set K gives a barycenter whose density with respect to P~ is at most 1/P(K).
Lemma 8. If learning fails with probability above r′, there is a set B with μ(B)>0 on which the failure persists, and every belief concentrated near B acts almost optimally on B.
Necessity. Without totally unbounded signals, some u∈U has no asymptotic learning.
Significance
Theorem 1 separates two questions. Whether learning succeeds for a particular utility depends on the utility, and for finite actions the paper gives weaker sufficient conditions (its Theorems 2 and 3). Whether learning succeeds for every utility depends on the information structure alone, and total unboundedness is exactly the condition. The theorem extends Smith and Sorensen's binary characterization to arbitrary separable complete metric state spaces. It also shows that their convex-hull condition does not carry over to more states. Lemma 1, the vanishing value of private information, holds in every sequential social learning game and is a reusable tool beyond this paper.
The results are proved on paper. To our knowledge none of them has a machine-checked proof. Formalizing them requires:
the equilibrium play measure, built by the Ionescu-Tulcea theorem;
martingale convergence of public beliefs;
a regular-conditional-probability treatment of beliefs over general Polish state spaces.
The paper's proofs pass over some of this measure theory, for example the measurability of best replies and of the maximal posterior payoff. A formal proof has to supply it.
Difficulty
Necessity is a direct construction. The difficulty is sufficiency. The natural first idea is to copy the binary argument: the public belief is a martingale, so it converges, and at a limit belief in the interior of the simplex unbounded signals would give private information positive value. With a general state space and a general utility this breaks in two places:
the limit belief μ∞ is a random measure on Ω, not a number in (0,1);
"positive value of a signal" has to hold uniformly over agents who face different histories.
The paper handles this with three ingredients:
Lemma 2, whose proof needs an imitation argument through Lemma 4's finite discretization;
Lemma 8, which localizes the failure of learning to a small set B;
Lemmas 5 and 7, which transfer total unboundedness from the prior to an averaged limit belief.
Each step uses facts about conditional distributions that hold only almost everywhere.
Formalization scope
Everything is in the namespace SeqSocialLearning.Adequate on top of Mathlib. The conventions are:
Agents are indexed from 0: Lean's agent t is the paper's agent t+1 and sees Fin t → A.
Strategies are Markov kernels Kernel ((Fin t → A) × S) A, so they are mixed, as on p. 4. Restricting to pure strategies would change the equilibrium set.
The play measurePσ is constructed with Kernel.trajMeasure on trajectories of triples (ω,st,at). The state is copied into every coordinate.
PBE is the paper's ex-ante definition: no profitable unilateral measurable deviation. No off-path optimality is imposed.
Posteriors use Mathlib's posterior and condDistrib. "The support contains 0 and 1" is stated in the equivalent ε-form.
Degeneracy of a prior is the equality of the one-signal decision value with the no-signal value.
Lemma 8'sminω′∈B is an infimum.
Lemma 4's product partition is a grid given by two finite measurable maps.
Nonempty Ω appears as a hypothesis; it follows from μ being a probability measure.
Two formalizations would make the goal trivial, and both are excluded:
quantifying total unboundedness over all events, including those of prior probability 0 or 1, would make it always false;
positing the play measure by a predicate instead of constructing it would let it be arbitrary.
Contributions are welcome on two kinds of work:
reusable infrastructure: measurable selection of best replies on compact action spaces, martingale convergence of condDistrib along a filtration, and the identification of the play measure's marginals;
the proofs of the eight milestones and of the goal.
Theorems 2 and 3 (finite actions), Lemma 6, the second part of Lemma 5 and Example 1 are not part of this mission.
S. Bikhchandani, D. Hirshleifer and I. Welch, A Theory of Fads, Fashion, Custom, and Cultural Change as Informational Cascades, Journal of Political Economy 100(5), 1992. https://doi.org/10.1086/261849
A. V. Banerjee, A Simple Model of Herd Behavior, Quarterly Journal of Economics 107(3), 1992. https://doi.org/10.2307/2118364
D. Rosenberg, E. Solan and N. Vieille, Informational Externalities and Emergence of Consensus, Games and Economic Behavior 66(2), 2009. https://doi.org/10.1016/j.geb.2008.07.007
Finite State Mean Field Games with Wright–Fisher Common Noise as Limits of N-Player Weighted Games: Normalized N-Player Equilibrium Values Converge to the Master Equation at a Polynomial RateResearch Paper
Motivation
Mean field games describe Nash equilibria of games with very many symmetric players by a single limiting object. Whether, and how fast, the equilibria of the N-player game converge to that limit is the convergence problem. For games with common noise, a source of randomness shared by all players, the standard route goes through the master equation, a PDE on the space of probability measures whose solution U gives the equilibrium value of a representative player as a function of the population's distribution (Cardaliaguet, Delarue, Lasry and Lions, 2019, arXiv:1509.02505).
For finite state spaces the distribution is a point of the simplex. Bayraktar, Cecchin, Cohen and Delarue (J. Math. Pures Appl. 147, 2021) introduced a common noise of Wright–Fisher type, the noise of population genetics, which keeps the distribution in the simplex and makes the master equation uniquely solvable. Its N-player counterpart is not obvious, because the Wright–Fisher diffusion is the limit of a resampling of the population, not of independent player noises. The paper behind this mission (Bayraktar, Cecchin, Cohen and Delarue, Mathematics of Operations Research, 2022, arXiv:2012.04845) builds such a game, in which players carry weights that are resampled at random times, and proves that its normalized equilibrium values converge to U at a polynomial rate.
Setting
There are d≥1 states [[d]]={1,…,d} and N≥1 players. Player l is at a state xl and carries a weight yl≥0; the weights are rational and sum to N, and the set of such weight tuples is Y. The players interact through the weighted empirical measure
μx,yN[i]=N1l∑yl1{xl=i},
a point of the simplex Sd−1={p∈Rd:pi≥0,∑ipi=1}.
Each player controls the rates at which it jumps between states, at quadratic cost. A drift φ(μ[j]) pushes players towards nearly empty states, where φ is non-increasing and Lipschitz, equal to κ on [0,δ] and to 0 beyond 2δ. At rate εN the common noise draws S∼Multinomial(N,μx,yN) and replaces every weight yl by ylS[xl]/(Nμx,yN[xl]). Players pay ylf(xl,μ) per unit time and ylg(xl,μT) at time T.
The equilibrium values divided by the weights, wN,l(t,x,y), solve the normalized Nash system: a countable family of ordinary differential equations in t, indexed by (l,x,y), with the Hamiltonian H(i,u)=−21∑j=i(ui−uj)+2 and the terminal condition wN,l(T,x,y)=g(xl,μx,yN).
The limit object is U=(U1,…,Ud) on [0,T]×Sd−1, solving the master equation
with Ui(T,⋅)=gi. Here di are the intrinsic derivatives on the simplex, dih(p)=d1∑m∂ei−emh(p), and the last term is the Kimura operator of Wright–Fisher diffusions. The mean-field prediction for player l is zN,l(t,x,y)=Uxl(t,μx,yN).
Formalization targets
Goal: Theorem 2.4
For κ large enough, there are χ>0 and C such that for every N, l, t, x and y∈Y,
The exponent χ is existential, so an improvement of the rate does not invalidate the target.
Milestones
Proposition 3.1: the normalized Nash system has a unique bounded solution.
Proposition 3.2: that solution is bounded by T∥f∥∞+∥g∥∞, uniformly in N.
Lemma 4.5: moment and deviation bounds for (Sμ[i]/(Nμ[i]))ℓ−1.
Lemmas 4.8–4.11: the common-noise, discrete-gradient, drift and Hamiltonian terms of the Nash system along zN,l against those of the master equation.
Proposition 4.6 with Remark 4.7: zN,l solves the Nash system up to a rest O(N−η) on the set TN where miniμx,yN[i]≥N−ϵ and maxlyl≤21N1−ϵ.
Significance
Theorem 2.4 gives a convergence rate for a finite-state mean field game whose common noise keeps the population on the simplex. It shows that the master equation of the earlier paper is the right limit of a concrete finite game, and with a diffusion-approximation argument it yields the paper's second main result, the weak convergence of the equilibrium empirical measure to the Wright–Fisher-type SDE (Theorem 2.6). The bound is not uniform in (x,y): it degenerates near the boundary of the simplex and for spread-out weights, which reflects the singular coefficient 1/μ[xl] in the resampling.
None of these results has a machine-checked proof. The mission produces a formal statement of the convergence problem for a game with weighted players, a formal Nash system on a countable state space, intrinsic derivatives on the simplex, and multinomial moment estimates.
Difficulty
The usual proof of convergence plugs zN,l into the Nash system, controls the rest, and closes with a Gronwall argument along the N-player dynamics. Here the rest is small only on TN, away from the boundary of the simplex and for moderate weights; the second derivatives of U may blow up at the boundary, and the resampling ratio S[i]/(Nμ[i]) is unbounded as μ[i]→0. A uniform estimate of supx,y∣zN,l−wN,l∣ is therefore not available, and the comparison must be weighted by functions of (x,y) whose expectations along the dynamics can be controlled. The weights are unbounded and enter through their moments.
Formalization scope
States are Fin d and players Fin N, with d≥1 and N≥1 (no N≥2). Weights are rational, and every equation and conclusion is required only for y∈Y and t∈[0,T]. The simplex is Mathlib's stdSimplex. Intrinsic derivatives use lineDeriv; multinomial expectations are finite sums over Finset.Nat.antidiagonalTuple. The resampling ratio is 1 when μ[xl]=0, as in the paper. Solutions of the Nash system and of the master equation are classical in t, with one-sided derivatives at 0 and T. The Nash system is solved by any bounded classical solution; Proposition 3.1 says there is exactly one.
The formalization commits to the following readings of the source.
ε for ε2. The master equation carries ε and ε/2 in front of the common-noise terms where (2.3) prints ε2 and ε2/2. The N-player noise rings at rate εN, and the limit SDE (2.4) has intensity ε; with ε2, Proposition 4.6 and Theorem 2.4 would be false.
Regularity class. The Wright–Fisher Hölder space of Proposition 2.1, which the paper only describes, is replaced by an explicit class: U continuous on [0,T]×Sd−1, differentiable along the simplex in the interior together with dU, with U and dU bounded, and dU and pjpkdj,k2U bounded and γ′/2-Hölder. These are the properties §7 says the proofs use. The Hölder assumptions on f,g are replaced by boundedness.
Constants.χ is chosen after d and γ′, the threshold κˉ0 after d,ε,T and the bounds of f,g, and C after all model data and U, before N and the point.
Remark 4.7. Lemmas 4.9–4.11 and Proposition 4.6 are stated for zN,l, so no division by a vanishing weight occurs.
Two printed slips. Lemma 4.9's (4.28) uses dUj, not dUxl, and Lemma 4.5's (4.19) assumes Nμ[i]≥1, without which it fails.
The bound is not trivial. Since ∣z−w∣≤2(T∥f∥∞+∥g∥∞) and the right side is at least of order CN−χ, a constant C chosen after N would make the theorem empty; C precedes N. The functions w and U enter only through the Nash system and the master equation, and the class of U is not empty: for d=1 the function U1(t)=g1+(T−t)f1 is in it.
The mission takes U as given. Its existence for d≥2 is Theorems 3.2 and 3.4 of Bayraktar, Cecchin, Cohen and Delarue (2021), quoted as Proposition 2.1, which is not a target. Also outside the mission are the N-player game itself (Propositions 2.2–2.3), Theorem 2.6, and the process estimates of §4.1 (Theorem 4.1, Remark 4.2, Propositions 4.3–4.4), which need a continuous-time Markov chain with a martingale-problem characterization. Contributions are welcome on every milestone; the multinomial estimates and the simplex calculus are reusable elsewhere.
Selected references
E. Bayraktar, A. Cecchin, A. Cohen, F. Delarue, Finite state mean field games with Wright–Fisher common noise as limits of N-player weighted games, Mathematics of Operations Research, 2022. arXiv:2012.04845v2
E. Bayraktar, A. Cecchin, A. Cohen, F. Delarue, Finite state mean field games with Wright–Fisher common noise, Journal de Mathématiques Pures et Appliquées 147:98–162, 2021. search
P. Cardaliaguet, F. Delarue, J.-M. Lasry, P.-L. Lions, The Master Equation and the Convergence Problem in Mean Field Games, Annals of Mathematics Studies 201, Princeton University Press, 2019. arXiv:1509.02505
C. L. Epstein, R. Mazzeo, Degenerate Diffusion Operators Arising in Population Biology, Annals of Mathematics Studies 185, Princeton University Press, 2013. search
S. N. Ethier, T. G. Kurtz, Markov Processes: Characterization and Convergence, Wiley, 1986. doi:10.1002/9780470316658
A Dynkin Game on Assets with Incomplete Information on the Return: The Entry Times of the Two Stopping Regions Form a Saddle Point When k > 0Research Paper
Motivation
A game option (or Israeli option, introduced by Kifer, Game options, 2000) is an American option whose seller may also cancel the contract at any time, paying the exercise payoff plus a penalty. Pricing it is a Dynkin game: a zero-sum optimal stopping game in which the buyer chooses an exercise time to maximise the discounted payoff and the seller chooses a cancellation time to minimise it. Such contracts model, for instance, a public concession that the granting authority may revoke against compensation.
In the classical theory the drift of the underlying asset is known. De Angelis, Gensbittel and Villeneuve (arXiv:1705.07352v4, 2019; Mathematics of Operations Research, 2021) study the perpetual game call option when neither player knows whether the asset pays a dividend: the return is r−δ0D with an unobservable D∈{0,1}. Both players learn about D only by watching the price. The paper's headline results are that the game has a value, that the value is globally C1, and that under explicit conditions on the parameters the game has a Nash equilibrium in pure strategies given by entry times into the two players' stopping regions. This mission formalizes the equilibrium result for the case k>0.
Timeline. Kifer (2000) introduced game options. Ekström and Villeneuve (AAP 2006) proved existence of the value for one-dimensional diffusions without uniform integrability and analysed the perpetual call with y=0. Ekström and Peskir (SICON 2008) proved existence of the value and of Nash equilibria for Markov processes under uniform integrability. Yam, Yung and Zhou (Game call options revisited, Mathematical Finance, 2012) treated the perpetual call with a known dividend (y=1). De Angelis, Gensbittel and Villeneuve (2017–2019) treated the case with an unknown dividend, y∈(0,1).
Setting
Fix a discount rate r>0, a dividend jump δ0>0, a volatility σ>0, a strike K>0 and a penalty ε0>0, and set G1(x)=(x−K)+ and G2(x)=(x−K)++ε0. On a probability space carrying a Brownian motion B and an independent random variable D∈{0,1} with P(D=1)=y, the asset is
Stx=xexp(σBt+(r−δ0D)t−2tσ2).
The players observe only S: their stopping times are stopping times of FS, the filtration generated by S and augmented by the null sets, with values in [0,∞]. For such τ (buyer) and γ (seller) the expected payoff is
with payoff 0 on {τ∧γ=∞}. The upper value is V=infγsupτM and the lower valueV=supτinfγM. A saddle point is a pair (τ∗,γ∗) with M(τ,γ∗)≤M(τ∗,γ∗)≤M(τ∗,γ) for all τ,γ.
The filterYt=E[D∣FtS] is the players' posterior belief that a dividend is paid; it has the closed form Yt=yℓt/(yℓt+1−y) with ℓt the likelihood ratio of the observed path. The pair (X,Y)=(S,Y) is Markov. With V=V the value, the stopping regions are S1={V=G1} (buyer) and S2={V=G2} (seller), and the entry times are
Finally k:=r−σ2/2−δ0/2. In the coordinates Zt=lnXt+δ0σ2ln1−YtYt one has Zt=Z0+kt, so the sign of k decides whether this time-like variable increases or decreases.
Formalization targets
Goal: Proposition 7.3 (p. 26)
Under Assumption 7.1, σ2/δ0>1, and k>0: for every x>0 and y∈(0,1),
The goal asserts the shape of the equilibrium (entry times into the model's own stopping regions), with no constant to be improved.
Milestones
In attack order: the filter identity (§2, p. 5); the regularity of V,V (Lemma 3.1); the value, the seller's optimality, and the super- and submartingale properties (Theorem 3.2); positivity of V (Remarks 3.3); the location of the stopping regions (Proposition 4.1, Lemma 4.2 (i), (iii)); the boundaries b1,b2 (Corollary 4.5) and finiteness of b1 (Lemma 4.6); the parabolic change of coordinates (5.1)–(5.4); connectedness of the regions in parabolic coordinates (Lemmas 5.2, 5.4); the boundary c1 of S1′ (Corollary 5.5, (5.20)); and positivity of c1 (Lemma 7.2).
Significance
The result. Existence of the value does not give an equilibrium: without uniform integrability the buyer's natural candidate τ∗ may fail to be optimal, and in the complete-information case y=0 the pair (τ∗,γ∗) is indeed not a saddle point (Proposition C.1). Proposition 7.3 shows that partial information about the dividend restores a pure-strategy equilibrium when k>0, and identifies it: each player stops on entering their own stopping region, which is described by monotone free boundaries. This makes the option's hedging and cancellation policies explicit.
Formalizing it. The result is proved on paper; nothing of it is machine-checked. The mission provides a Lean model of a continuous-time game with an unobservable parameter, built from Mathlib's Brownian motion, stopping times and martingales, and a decomposition of the proof into the paper's own lemmas. The theory the proofs use (optimal stopping games for Markov processes, Itô's formula, Girsanov's theorem) is largely absent from Mathlib.
Difficulty
The obvious argument fails at the buyer's side. The seller's optimality comes with the value (Theorem 3.2); for the buyer one needs V(x,y)≤M(τ∗,γ) for every γ, i.e. that the submartingale e−r(t∧τ∗)V can be sampled at τ∗∧γ. The payoffs are not uniformly integrable (Esupte−rtGi(St)=+∞), so optional sampling cannot be applied directly and the term E[e−rtV(Xt,Yt)1{t≤τ∗∧γ}] must be shown to vanish. That requires a pathwise bound on the price before τ∗, which in turn needs the geometry of S1 in both coordinate systems, monotonicity of Z (k>0), and strict positivity of the boundary c1 (Lemma 7.2, which needs σ2/δ0>1). Theorem 3.2 itself rests on the Ekström–Peskir theory of Markovian Dynkin games applied to a degenerate two-dimensional diffusion.
Formalization scope
Model. The paper's working process (X,Y) solves an SDE ((2.3)); Mathlib has no stochastic integral. The mission therefore states everything on the paper's original space (§2): X:=S as the explicit geometric Brownian motion (2.1), and Y the closed-form Bayes filter, which has the same law as the paper's pair (p. 6). The item filter_condexp ties Y to E[D∣FtS] and records its observation measurability and continuous paths, while eq_5_3 ties the encoding to the paper's coordinates (Z,Y).
Probability space. A Setting carries a measure P, a Brownian motion B with all paths continuous (a continuous modification), and U uniform on [0,1] independent of B; then Dy=1{U<y} realises every prior y on one space, so V(x′,y′) is defined at every state.
Information. The ambient probability measure is complete. Thus FtS is generated by (Ss)s≤t and all its null sets; it depends on (x,y). Stopping times take values in WithTop ℝ≥0, so infinite stopping times are allowed.
Payoff and values.M is a Lebesgue integral (lintegral) in [0,∞], with payoff 0 on {τ∧γ=∞}; V,V∈[0,∞], and V:=V as a real number (finite since V≤G2(x)).
Entry times are Mathlib's hittingAfter … 0 into {G2≤V} and {V≤G1}, equal to +∞ if never reached.
Boundaries.b1,b2 take values in [0,∞] with inf∅=+∞, sup∅=0. c1(z)=sup({0}∪{y∈(0,1)∣(z,y)∈S1′}) is an explicit choice of the boundary whose existence Corollary 5.5 asserts.
Domains and readings. The goal is posed for x>0, y∈(0,1), which the proof needs; it is false at y=0. Remarks 3.3 is posed for x>0 (V(0,y)=0). For the same reason the stopping regions and the boundaries b1,b2 live on (0,∞)×[0,1]: at the degenerate price x=0 one has V(0,y)=0=G1(0), and the page's literal x∈[0,∞) in (4.2) and (4.11) would give b1≡0. Lemma 3.1 (iii)'s "a given constant C" is read as an existential constant. The standing assumption k=0 of §5 is a hypothesis of Lemmas 5.2, 5.4, Corollary 5.5 and Lemma 7.2. Only clause (5.20) of Corollary 5.5 is formalized. Proposition 4.1 uses the explicit forms of A1,A2 printed on p. 9.
No trivialization. The filtration is that of the price, never one containing D (with D observed the game is the complete-information one, where the result fails); Y is the posterior, not the prior; V in the entry times is the game's own value, never a free function; stopping times include infinite ones; and the payoff is not a Bochner integral, which would silently be 0 for non-integrable payoffs.
Not formalized. Proposition 6.4 (C1 regularity), Proposition 7.6 (the case k<0), Proposition 5.6 and Lemmas 4.3, 4.4, 4.7 are not used by the goal's proof.
Infrastructure needed. Optional sampling for continuous-time super/submartingales with unbounded stopping times; hitting times of closed sets by continuous processes as stopping times of an augmented filtration; filtering for Brownian motion with a two-point drift; Itô's formula and Girsanov's theorem for the lemmas of §4–§7; existence of the value for Markovian Dynkin games. Each is reusable well beyond this mission, and contributions of any of them are welcome.
Selected references
T. De Angelis, F. Gensbittel, S. Villeneuve, A Dynkin Game on Assets with Incomplete Information on the Return, arXiv:1705.07352v4, 2019; Mathematics of Operations Research, 2021. https://arxiv.org/abs/1705.07352v4
E. Ekström, G. Peskir, Optimal stopping games for Markov processes, SIAM Journal on Control and Optimization 47, 684–702, 2008. https://doi.org/10.1137/060673916
S. C. P. Yam, S. P. Yung, W. Zhou, Game call options revisited, Mathematical Finance, 2012.
R. S. Liptser, A. N. Shiryaev, Statistics of Random Processes I: General Theory, 2nd ed., Springer, 2001.
On Uniform Exponential Ergodicity of Markovian Multiclass Many-Server Queues in the Halfin–Whitt Regime: Positive Spare Capacity Makes the Limiting Diffusion Uniformly Exponentially ErgodicResearch Paper
Motivation
Call centres and hospital wards are run in the Halfin–Whitt regime: the number of servers n exceeds the offered load by a margin of order n, so that delays are short and servers are nearly always busy. With several customer classes, the operator must also choose a scheduling policy: which waiting class is served next. In the diffusion limit, every such policy becomes a control of an m-dimensional piecewise-linear diffusion (Atar, Mandelbaum and Reiman 2004; Harrison and Zeevi 2004). Ergodic control of this diffusion, and the justification that its optimal policies are asymptotically optimal for the queues, both need the diffusion to be stable under every admissible policy, with bounds that do not depend on the policy.
A. Arapostathis, H. Hmedi and G. Pang (arXiv:1812.03528v3) prove that positive spare capacity alone gives this: the limiting diffusion is exponentially ergodic, at a rate and with constants that are uniform over all stationary Markov controls.
Timeline.
2004: Atar–Mandelbaum–Reiman and Harrison–Zeevi derive the controlled diffusion limit of the multiclass M/M/n queue in the Halfin–Whitt regime.
2013: Dieker and Gao prove positive recurrence and, with abandonment, exponential ergodicity of piecewise Ornstein–Uhlenbeck processes, for one specific static priority control.
2019: Arapostathis, Pang and Sandrić show positive recurrence under every constant control when the spare capacity ϱ is positive, and transience under Markov controls when ϱ<0 (without abandonment).
2019: this paper extends exponential ergodicity to all stationary Markov controls, with uniform constants.
Setting
Fix m≥1 classes, service rates μi>0, abandonment rates γi≥0, and write M=diag(μi), Γ=diag(γi), e=(1,…,1), βi=γi/μi. Control values lie in the simplex Δ={u∈Rm:u≥0,⟨e,u⟩=1}; ui is the fraction of the queue made up of class i. The spare capacity is ϱ∈R. The controlled drift is
b(x,u)=−mϱMe−M(x−⟨e,x⟩+u)−⟨e,x⟩+Γu,
and the diffusion is dXt=b(Xt,Ut)dt+σdWt with a constant positive diagonal σ=diag(ri) and a standard Brownian motion W. The generator under the control value u is Luf=21trace(σσT∇2f)+⟨b(⋅,u),∇f⟩.
A stationary Markov control is a Borel map v:Rm→Δ, used as Ut=v(Xt); the class is Usm. Its transition probability is Ptv(x,⋅), the law of Xt started at x. An invariant probability measureπv satisfies ∫Ptv(x,⋅)πv(dx)=πv.
The Lyapunov function is built from the convex C2cut-offψ (ψ=−21 on (−∞,−1], ψ(t)=t on [0,∞)): for ε,θ>0,
The rate γ and the constant C do not depend on the control: that uniformity is the content of the theorem.
Milestones, in the order the proof uses them
p. 5: under every v∈Usm the SDE has a unique strong solution and a strong Feller transition function.
p. 7: the properties of ψ, including ψε′′≤23ε2.
(2.13a)–(2.13b), p. 8: the identities for ⟨∇Ψε,b⟩ and ⟨∇Ψ(−⋅),b⟩.
Lemma 2.1, p. 8: the drift inequalities (2.17) for the truncated drift bc, c∈[1,∞].
(2.31), p. 10: 21trace(a∇2V)≤ε(23(ε+θ)+εCˉ)V.
(2.28), p. 10: the Foster–Lyapunov equation above (part (a) of the goal).
(2.10), p. 7: Ψε,θ∗ is comparable to ∥x∥1.
(2.34), p. 11: ∫Ptv(x,dy)V(y)≤κ0/κ1+e−κ1tV(x) whenever LuV≤κ0−κ1V.
p. 11: the uniform minorization P1v(x,B)≥δ0 on a ball, from the parabolic Harnack inequality.
Significance
The result. Uniform exponential ergodicity makes the ergodic control problem for the multiclass Halfin–Whitt diffusion well posed over all stationary Markov policies: invariant measures form a tight family with uniformly integrable exponential moments, and long-run averages converge at a policy-independent geometric rate. In the paper it is the base for the prelimit results of Section 3, where the drift inequalities of Lemma 2.1 are compared with those of the n-server queues. With Γ=0 the condition ϱ>0 is sharp, since ϱ<0 gives transience under every Markov control (Arapostathis, Pang and Sandrić 2019).
Formalizing it. The theorem is proved in the paper; none of it is formalized. The mission produces a machine-checked version with two printed slips corrected (see Formalization scope). Its analytic half, Lemma 2.1, (2.13), (2.10), (2.31) and the Foster–Lyapunov equation, is finite-dimensional calculus and is reachable now. The probabilistic half needs Itô's formula for SDEs with additive noise, existence and uniqueness for SDEs with Borel drifts, and a Harnack-type minorization, none of which Mathlib has.
Difficulty
The control enters the drift only through ⟨e,x⟩+u, but a Markov control may switch arbitrarily, so no single linear or quadratic Lyapunov function works for all policies; Dieker and Gao's common quadratic Lyapunov functions cover one static priority only. The function V mixes an ε-scaled workload term with a θ-weighted idleness term, and the drift inequality must hold on both half-spaces {⟨e,x⟩≤0} and {⟨e,x⟩≥0} with explicit constants, which forces case splits over the classes with βi≤1 and βi>1.
For part (b), a Foster–Lyapunov inequality alone gives ergodicity for each control but not a control-independent rate: the constants of the standard Meyn–Tweedie argument depend on a minorization on a small set, and that minorization must be shown uniformly over all Borel controls, for which there is no regularity in v to exploit.
Formalization scope
Representation. States are EuclideanSpace ℝ (Fin m); classes are Fin m; control values are vectors in stdSimplex ℝ (Fin m). The model data is the published ManyServerQED.Scheduling.DiffusionData with its drift ℓ+(μ−γ)⟨e,x⟩+u−μx; the SDE equation, invariant measures, generator and weighted total-variation norm use published GPhNGIDiffusion.EMError definitions. The local IsStrongSolution and IsMarkovTransition specify measurable, adapted solutions and their transition laws. The cut-off ψ is the published cutoff.
Committed conventions.
The drift is centred, ℓ=−mϱMe (Remark 2.3, "without loss of generality").
σ=diag(ri) is constant (the proof uses σ=diag(2λ~i)1/2), and ∑iri2/(2μi)=1, i.e. ∑iλ~i/μi=1, the identity the proof invokes in (2.31). Without it the Foster–Lyapunov equation fails numerically.
θ(βmax−1)+≤1 is written as θ(βi−1)+≤1 for all i.
A transition function is the law of every measurable strong solution on every probability space; "positive recurrent" is read as the existence of an invariant probability measure. The published transition predicate quantifies over arbitrary pathwise solutions, including nonmeasurable ones for which Measure.map returns a junk zero measure, so a local predicate is used here.
The V-norm uses the weight 1+V, since V can be below 1 while the paper's f-norm needs f≥1; the norm is +∞ unless both measures integrate 1+V.
The truncation level c∈[1,∞] of Lemma 2.1 is an extended real.
Printed slips, corrected.
(2.28) is quoted as printed in its milestone but posed with 2θ in place of θ on ∥x−∥1: as printed it fails along x=−Re, R→∞, where the left side exceeds the right by ε(θϱ+2mϱ+εθ2Cˉ)V→∞. The rest of the theorem is unaffected.
The lower bound of (2.10) is quoted as printed but posed with the constant m(1∨εθ)/(2μmin) in place of 2m, which fails when μmin<1.
Ruling out trivial readings. The transition-function hypothesis is not vacuous: milestone 1 asserts existence of solutions and of a transition function under every Borel control. The constants γ,C are quantified before the control, the Foster–Lyapunov inequality ranges over u∈Δ only, the norm is never a junk value, and the moment bound on πv uses the same κ0 as the Foster–Lyapunov equation.
Constructed ball. The minorization uses an existential ball whose complement satisfies the inequality chosen after (2.34), with its lower bound uniform over controls.
Needed infrastructure and welcome contributions. Calculus for ψ and V (reusable for every Lyapunov argument built on this cut-off); Itô's formula with localization for additive-noise SDEs; strong existence and pathwise uniqueness for Borel drifts of linear growth; a Harnack or Girsanov lower bound for transition probabilities on balls; the Meyn–Tweedie computable geometric bound. Each is reusable beyond this mission, and proofs of any milestone are welcome independently.
Selected references
A. Arapostathis, H. Hmedi, G. Pang, On uniform exponential ergodicity of Markovian multiclass many-server queues in the Halfin–Whitt regime, arXiv:1812.03528v3, 2019 (Mathematics of Operations Research, 2021). https://arxiv.org/abs/1812.03528
A. Arapostathis, G. Pang, N. Sandrić, Ergodicity of a Lévy-driven SDE arising from multiclass many-server queues, Ann. Appl. Probab. 29(2), 2019. https://doi.org/10.1214/18-AAP1430
R. Atar, A. Mandelbaum, M. I. Reiman, Scheduling a multi class queue with many exponential servers: asymptotic optimality in heavy traffic, Ann. Appl. Probab. 14(3), 2004. https://doi.org/10.1214/105051604000000233
A. B. Dieker, X. Gao, Positive recurrence of piecewise Ornstein–Uhlenbeck processes and common quadratic Lyapunov functions, Ann. Appl. Probab. 23(4), 2013. https://doi.org/10.1214/12-aap870
M. Gruber, Harnack inequalities for solutions of general second order parabolic equations and estimates of their Hölder constants, Math. Z. 185(1), 1984. https://doi.org/10.1007/BF01214972
I. Gyöngy, N. Krylov, Existence of strong solutions for Itô's stochastic equations via approximations, Probab. Theory Related Fields 105(2), 1996. https://doi.org/10.1007/BF01203833
J. M. Harrison, A. Zeevi, Dynamic scheduling of a multiclass queue in the Halfin–Whitt heavy traffic regime, Oper. Res. 52(2), 2004. https://doi.org/10.1287/opre.1030.0084
S. P. Meyn, R. L. Tweedie, Computable bounds for geometric convergence rates of Markov chains, Ann. Appl. Probab. 4(4), 1994. https://doi.org/10.1214/aoap/1177004900
Quantum Approximate Counting, Simplified: A Grover-Only Algorithm Estimates K to Within Relative Error ε with Probability 1 − δ Using O(√(N/K) ε⁻¹ log(1/δ)) QueriesResearch Paper
Motivation
Approximate counting asks for the number K of marked items in a list of N items, to within a multiplicative factor 1±ε, when the marked set can only be accessed through a membership oracle. Brassard, Høyer, Mosca and Tapp (arXiv:quant-ph/0005055) gave a quantum algorithm that uses O(N/Kε−1) queries, which is optimal. Their algorithm combines Grover iterations with phase estimation, and so it needs the quantum Fourier transform (QFT). Counting is a basic subroutine in quantum algorithms for estimation, numerical integration and Monte Carlo speed-ups, so the circuit it requires matters.
Aaronson and Rall (arXiv:1908.10846, SOSA 2020) give an algorithm with the same query complexity that uses only Grover iterations and classical post-processing of measurement outcomes. Its analysis is short and elementary, and every constant in it is explicit. This mission formalizes their Theorem 1, the rotation lemma it rests on, and the numbered estimates of its proof (arXiv:1908.10846v6, pp. 4–10).
Timeline (as recounted on pp. 1–3 of the paper):
1996: Grover introduced the iteration whose powers rotate the uniform superposition towards the marked items.
1998: Brassard, Høyer, Mosca and Tapp gave counting with O(N/Kε−1) queries by phase estimation of the Grover operator. The journal version appeared in 2002.
1997–1999: Grover outlined a broadly similar Grover-only approach. Abrams and Williams discussed it further, for estimating integrals to additive error with O(ε−1) queries. Neither proved the approach correct.
2019: Suzuki et al. gave an O(ε−1N/K)-query QFT-free algorithm that uses maximum-likelihood post-processing, supported by numerical experiments rather than a proof. Wie sketched a QFT-free algorithm based on Hadamard tests, also without proof.
2019–2021: Aaronson and Rall proved a Grover-only algorithm correct with the optimal O(N/Kε−1logδ−1) bound. The constants used here are those of v6.
Setting
Let S⊆[N]={0,…,N−1} be a nonempty marked set with K=∣S∣. The algorithm first pads the list to N′=(106+1)N items, the added items being unmarked, and then works on CN′. The membership oracle is the diagonal unitary U∣x⟩=(−1)[x∈S]∣x⟩. The uniform superposition is ∣ψ⟩=N′−1/2∑x∣x⟩, and the Grover operator is
G=(I−2∣ψ⟩⟨ψ∣)U.
Preparing G(r−1)/2∣ψ⟩ for an odd r costs (r−1)/2 queries. Measuring it in the computational basis yields a marked item with probability p given by the Born rule. With θ=arcsinK/N′≤0.001 this probability is sin2(rθ), so m repetitions give a Bin(m,sin2(rθ)) count of marked outcomes.
The algorithm has two loops. Step 1 tries rk, the largest odd integer at most 1.05k, for k=0,1,…, with ⌈5000ln(5/δ)⌉ measurements per round. It stops at the first k=:kend at which at least 95% of the outcomes are marked. Step 2 keeps a bracket [θmin,θmax], initially θmin=0.9⋅1.05−kend and θmax=1.65θmin. In round t it computes rt from the bracket by the rotation lemma (Lemma 2):
k=round(2(θmax−θmin)θmin),r=closest odd integer to θminπk.
It then measures ⌈250lnδt−1⌉ times, with δt=(δε/65)0.9−t, and shrinks γ=θmax/θmin−1 by the factor 0.9: it raises θmin if at least 12% of the outcomes are marked and lowers θmax otherwise. It stops once θmax≤(1+ε/5)θmin and returns K^=N′sin2θmax.
Formalization targets
Goal: Theorem 1
For nonempty S⊆[N], K=∣S∣:
for 0<ε,δ<1 there is a fuel bound F such that
Pr[the algorithm outputs K^ with K(1−ε)<K^<K(1+ε)]≥1−δ;
there is an absolute constant C such that for 0<ε,δ≤1/2, every run that outputs an accurate K^ has made fewer than CN/Kε−1lnδ−1 queries.
Milestones
In the order of the proof:
eq. (1), the Grover state and its marked probability sin2(rθ);
eqs. (2)–(6), the ideal bracket after step 1;
eqs. (7)–(12) and (15)–(28), no late and no early termination of step 1, with the Kullback–Leibler form (13)–(14) of the Chernoff bound;
eqs. (52), (57), (58)–(66) and (68)–(74), the estimates behind the rotation lemma;
Lemma 2 itself;
eqs. (30)–(32), the union bound over step 2;
footnote 2, the conversion from angle to count;
eqs. (37)–(44), the arithmetico-geometric sum in the query bound.
Significance
Theorem 1 shows that the optimal query complexity of approximate counting is achieved without the quantum Fourier transform. The circuits involved are powers of the Grover operator followed by computational-basis measurements, and the classical control is a pair of loops with explicit thresholds. The same argument gives QFT-free amplitude estimation (Theorem 3 of the paper), and its rotation lemma is a reusable tool for any algorithm that brackets a rotation angle by biased coin tosses.
The result is proved on paper but not formalized. The platform has Grover's iteration for a single marked item; it has no Grover state for a marked set, no counting algorithm and no query-complexity statement of this kind. A complete development would provide a machine-checked quantum query algorithm whose output law is derived from the Born rule and the actual Grover operator. It would also provide a Chernoff bound in Kullback–Leibler form for the binomial distribution, which Mathlib does not have.
Difficulty
The deterministic estimates are elementary but numerous, and each constant (0.96, 7.5%, 16.5%, 1.102, 0.27, 0.42) must be checked against interval arithmetic for sin and π. Two printed steps are numerically loose, (8) and (26), and must be replaced by correct intermediate values.
The probabilistic part has two difficulties. First, the success event of step 1 is not a single tail bound. Termination too early is controlled by a union bound over rounds k≤k0, and the additive Hoeffding bound is too weak there: it needs the KL form (13). Second, step 2's output law is defined by a recursion whose bracket is random, and Lemma 2 applies to a round only while the bracket still contains θ. The union bound over rounds has to be carried through this recursion. Third, the paper bounds the queries only on its "ideal" event. Part 2 of the goal requires the bound on every run that ends with an accurate estimate, including runs in which step 1 stopped late.
Formalization scope
States are vectors Fin N' → ℂ, operators are Matrix (Fin N') (Fin N') ℂ, and ∣ψ⟩⟨ψ∣ is vecMulVec ψ (star ψ). The padding embeds Fin N as the first N indices of Fin N'.
The algorithm approxCount N S ε δ F : PMF (Option (ℝ × ℕ)) is defined, not quantified. Its outputs are some (K̂, queries), or none when either loop runs out of its fuel F. Fuel exhaustion counts as failure. Measurement counts are Mathlib's binomial law with the Born probability ∑x∈S∣⟨x∣Gj∣ψ⟩∣2 computed from U, ψ and G.
The algorithm sees S only through these probabilities. rk depends on k only, rt on the bracket only, and the repetition counts and thresholds on ε,δ,t only. A formalization in which the algorithm reads K or θ directly, or in which "there exists an algorithm" is quantified, is trivial (output K with zero queries) and is ruled out.
"At least m times" means exactly ⌈m⌉ times. Mathlib's round breaks ties upward, a choice the paper permits. The step-2 threshold is "at least 12%", as in the algorithm box; Lemma 2 is stated for both "more than" and "at least".
Eq. (1) is stated in corrected form: the printed identity misses the global sign (−1)(r−1)/2 coming from G=−(2∣ψ⟩⟨ψ∣−I)U.
The O(⋅) of Theorem 1 is one absolute constant chosen before N,S,ε,δ,F, in the regime ε,δ≤1/2. The rotation lemma's O(θ−1γ−1) is pinned to r≤5/(γθ).
Not formalized: the O(logN)-qubit space bound (there is no circuit model), Theorem 3 (amplitude estimation), and the unquantified remark on K=0.
Contributions of reusable infrastructure are welcome, in particular the KL-form Chernoff bound for ProbabilityTheory.binomial and lemmas on powers of a reflection product acting on a two-dimensional invariant subspace.
Selected references
S. Aaronson, P. Rall, Quantum Approximate Counting, Simplified, Symposium on Simplicity in Algorithms (SOSA) 2020; arXiv:1908.10846v6, 2021. https://arxiv.org/abs/1908.10846
G. Brassard, P. Høyer, M. Mosca, A. Tapp, Quantum Amplitude Amplification and Estimation, Contemporary Mathematics 305, 2002. https://arxiv.org/abs/quant-ph/0005055
Shipment Policies for Products with Fixed Shelf Lives: Impact on Profits and Waste: Shipping Oldest First Minimizes the Long-Run Average Warehouse Waste and the Total WasteResearch Paper
Motivation
Manufacturers of products with fixed shelf lives must decide both how much to send from a warehouse to a store and which age cohorts to send. A shipment can reduce expiry at the warehouse while changing the age mix seen by shoppers at the store. In the two-echelon system studied by Akkaş and Honhon, the store does not control which items each shopper takes: some shoppers select older items first and the remainder select newer items first. The resulting waste at the warehouse and store can move in different directions when the manufacturer changes a shipment. The paper uses this setting to ask whether ship oldest first (SOF), a simple age-selection rule, minimizes waste over an infinite horizon. Akkaş and Honhon, SSRN manuscript.
The paper reports an inventory-age audit of a consumer packaged goods supplier and then studies the shipment decision analytically. Its Theorem 1 separates warehouse waste from the sum of warehouse and store waste. The former has a specified SOF policy that ships in full; the latter has an SOF age-selection rule but allows the shipped quantity to vary with the state. This mission formalizes those two claims, with a correction to one printed clause described below. Akkaş and Honhon, §3 and Theorem 1.
Setting
A unit remains saleable for N periods, where N≥1. At the start of a period the warehouse holds xi units of age i, and the store holds yi units of that age, for i=1,…,N−1. Both age-zero entries are zero before replenishment. A plant then delivers fresh units so the warehouse's total post-replenishment stock reaches the exogenous order-up-to levelSW. Thus its age-zero stock is x0post=SW−∑i=1N−1xi. The manufacturer sends qi units of age i from the warehouse to the store, subject to warehouse availability and the store's shelf capacityL. In particular, qi≤xipost and ∑iqi+∑iyi≤L. Akkaş and Honhon, §4.1, pp. 10–13, and §4.2, p. 15.
The period's shopper population is λ. A fraction γ∈[0,1] selects inventory first in, first out (FIFO), meaning older units first; the remaining shoppers select it last in, first out (LIFO), meaning newer units first. The discrete FIFO population is ⌈γλ⌉. If zi=yi+qi is the store stock after delivery, the paper's age-specific demand is
Di(z,λ;γ)=(⌈γλ⌉−j>i∑zj)++(λ−⌈γλ⌉−j<i∑zj)+.
Here u+=max(u,0). Only units of age N−1 expire at the end of the period. Warehouse waste is WW=xN−1post−qN−1, store waste is WR=(zN−1−DN−1)+, and total waste is W=WW+WR. The remaining units advance one age before the next period. An SOF shipment uses a younger unit only when every older available unit has been shipped. An SOF shipment in full fills the remaining store space or ships the entire available warehouse stock. Akkaş and Honhon, §4.1, pp. 11–12, and §4.2.1, p. 15.
Formalization targets
Warehouse waste
Let WW,+ and WW,− denote the upper and lower limits of expected Cesàro averages of warehouse waste. The first sentence of Theorem 1 is the mission's sole milestone:
WW,+(SOF-full;s0)≤WW,−(σ;s0)for every feasible initial state s0 and every admissible policy σ.
Total waste and mission goal
The full goal also asserts that one stationary, feasible SOF policy f minimizes total waste from every feasible initial state:
∃fstationary and SOF,W+(f;s0)≤W−(σ;s0)for all feasible s0 and admissible σ.
The policy f chooses its shipping quantity from the current state; the theorem does not prescribe shipping in full for this second objective. These are waste objectives, independent of the paper's separate revenue, labor, and expiration-cost parameters. Akkaş and Honhon, Theorem 1, p. 15.
Significance
The warehouse result identifies a fixed, implementable shipment rule for minimizing waste at that location across the paper's whole FIFO/LIFO range. The total-waste result says that optimal age selection can remain SOF even when the quantity decision must respond to the store inventory and shopper behavior. Together they distinguish the two decisions in a setting where moving an expiring unit out of the warehouse can merely move the expiry to the store. The paper later uses its waste comparison in the analysis of profit policies under additional cost assumptions. Akkaş and Honhon, Appendix A.2–A.3.
The theorem is a published paper result, but this mission starts with open Lean statements rather than a machine-checked proof. A complete development would provide a reusable formal model of age-structured two-echelon inventory, stochastic shopper depletion, feasible shipment policies, and long-run waste criteria. The general controlled-process and Cesàro-average definitions can also support other finite inventory models.
Difficulty
An age-selection decision changes the stock carried into future periods, and a quantity decision changes the next warehouse replenishment. A comparison of immediate waste therefore does not by itself compare long-run averages. In particular, simply reordering each shipment to send older units first while keeping its quantity fixed can increase later warehouse waste for some policies; it cannot serve as an unrestricted policy-level argument. The store's FIFO/LIFO mixture introduces another intertemporal effect because a shipped younger unit can be sold ahead of an older one. Akkaş and Honhon, Appendix A.1, Appendix p. 1.
Formalization scope
Lean uses finite age vectors indexed by Fin N, natural-number stocks and shipments, a real γ constrained to [0,1], and an arbitrary probability mass function for shopper population. The standing condition N≥1 makes the final age cohort defined; SW=0 and L=0 remain allowed. Demand samples are independent and identically distributed across periods, as in the stationary transition model of §4.2. Policies are time-dependent deterministic Markov rules: at period t they choose a shipment from the current inventory state. This is the finite-state class used by the paper's dynamic program. The comparison includes every feasible rule in that class, every feasible initial state, and a common stationary SOF witness for total waste; it does not restrict competitors to SOF or full shipments. Standard finite-MDP results explain why stationary deterministic optimal policies suffice among broader history-dependent randomized classes. Puterman, Markov Decision Processes, Theorem 9.1.8.
Costs and expectations live in [0,∞], and long-run criteria use upper and lower limits of expected Cesàro averages. The paper prints a plain limit without specifying what happens when a policy's average fails to converge. The upper-versus-lower comparison states the stronger optimality claim without relying on that convergence. Lean time starts at zero. Natural-number subtraction implements the paper's positive-part operation, and the post-replenishment stock formula is used only for feasible states. The model follows the printed balance equations when discarding the final-age store stock, despite a tension with the narrative timing of expired-unit pickup.
Three transcription corrections are explicit. The demand formula's sums use zj, the running index, where the page prints zi; the action sum runs through age N−1; and the waste objectives are minimized despite a printed “max” in §4.2.1. More consequentially, Theorem 1's last clause says shipping is in full under pure FIFO. This is false and is omitted: for γ=1, N=2, SW=1, L=2, and a Bernoulli(1/2) shopper population, SOF in full has average total waste 1/2 while an SOF policy that holds stock back attains 1/4. The formal goal neither assumes this clause nor asserts its negation. Its feasible-policy and feasible-initial-state quantifiers rule out vacuous or artificially easy readings.
A complete proof may need finite-state inventory invariants, admissibility of the SOF-in-full rule, comparisons of age profiles over time, and results about average costs of controlled chains. Contributions establishing those reusable facts, or the stated warehouse and total-waste targets, are within scope. The separate profit theorems and numerical studies of the paper are outside this mission.
Selected references
Akkaş, A., and D. Honhon. Shipment Policies for Products with Fixed Shelf Lives: Impact on Profits and Waste. Authors' SSRN manuscript 3247290, 2021; published in Manufacturing & Service Operations Management. SSRN manuscript; journal DOI.
Puterman, M. L. Markov Decision Processes: Discrete Stochastic Dynamic Programming. Wiley, 1994. DOI.
Game of Variable Contributions to the Common Good Under Uncertainty: The Symmetric Game Has a Regular Control Markov Perfect EquilibriumResearch Paper
Motivation
Many shared assets deteriorate unless someone pays to maintain them: a public infrastructure stock, a jointly used ecosystem, a joint venture's capital, a common research base. When several parties benefit from the same stock and each pays for its own contributions, every party is tempted to wait for the others to pay. Kwon (arXiv:1904.00500, published in Operations Research, 2022) models this as a two-player stochastic differential game in continuous time, in which each player may raise the stock at any rate, and studies the Markov perfect equilibria of the game.
The single-player version is a classical monotone follower problem of singular stochastic control (Alvarez 2001; Øksendal and Sulem 2005; Lon and Zervos 2011): the optimal policy reflects the stock upward at a threshold. In the game, a natural first guess is that free riding lowers that threshold, or that the players take turns reflecting. The paper shows instead that, in the symmetric game, there is an equilibrium in which both players contribute gradually, at a finite rate, below the same threshold as the single decision maker. The equilibrium is a continuous-state analogue of the mixed-strategy delay equilibrium of the war of attrition.
Setting
The stock Zt lives in an interval I=(a,b)⊆R (a=−∞ or b=∞ allowed) and satisfies
dZt=μ(Zt)dt+σ(Zt)dWt+dξ1t+dξ2t,
on a filtered probability space satisfying the usual conditions, where W is a Wiener process, μ<0 (the stock deteriorates) and σ>0. Player i's controlξi is any adapted, nondecreasing, right-continuous process: its cumulative contribution. A control may jump, including at time 0. Player i's payoff from the initial stock x is
Vi(x;ξ)=Ex[∫0∞e−rtπ(Zt)dt−∫0∞e−rtkdξit],
where π is an increasing, bounded-above flow profit, k>0 the unit cost and r>0 the discount rate. In the symmetric game both players have the same π and k.
The analysis uses:
the uncontrolled processXx, solving dX=μ(X)dt+σ(X)dW, X0=x;
the resolvent(Rg)(x)=E[∫0∞e−rtg(Xtx)dt];
the generatorA=21σ(x)2∂x2+μ(x)∂x−r;
the increasing and decreasing fundamental solutionsψ,φ of Aψ=Aφ=0;
the function q(x)=π(x)+[μ(x)−rx]k.
Assumptions 1–3 of the paper require Lipschitz coefficients with linear growth, a uniform integrability condition on e−rτ(Xτ)−, integrability of e−rtπ(Xt), and boundary behaviour of Rπ and q at a and b; in particular q increases below a point x∗ and decreases above it.
The single decision maker's value is
V(x)={k(x−θ)+V(θ)(Rπ)(x)+Aφ(x)x<θ,x≥θ,(4)
for a threshold θ∈I and a coefficient A that make V twice continuously differentiable (Lemma 1). The candidate equilibrium rate is
u(x)={−k1(AV(x)+π(x))0x<θ,x≥θ.(6)
Formalization targets
Goal: Theorem 2 (pp. 15–16)
If u of (6) is Lipschitz on I with linear growth, then the profile in which both players contribute at the regular rate dξit=u(Zt)dt is a Markov perfect equilibrium, and
V1(x;ξ)=V2(x;ξ)=V(x)for all x∈I.
The equilibrium is against all adapted deviations, from every initial state, with the opponent keeping its feedback rule u.
Milestones
Proof of Lemma 1, part (i) (pp. 33–34): for V of (4) in C2(I), AV+π=0 on [θ,b), AV+π<0 on (a,θ), V′=k on (a,θ] and V′<k on (θ,b).
Lemma 1 (p. 10): such θ and A exist, and V is the single decision maker's optimal value, attained by reflection at θ.
Theorem 1 (pp. 14–15), in the instance where both the opponent and the candidate contribute at regular feedback rates: a C2, nondecreasing, bounded-above U satisfying conditions (iii)–(iv) is the candidate's payoff and dominates every deviation in the class Σi.
The optimality condition of §2.2 (p. 10), for V of the form (4) and the operator Aξj=21σ2∂x2+[μ+uj]∂x−r of an opponent playing a regular rate uj (as applied on p. 38): conditions (i)–(ii) make V the optimal value over all controls.
Significance
Theorem 2 says that the free-rider problem in this game does not shift the threshold below which the stock is maintained; it changes how the stock is maintained, from instantaneous reflection by one party to gradual, rate-limited contribution by both. The equilibrium value equals the single decision maker's value, so the players lose nothing relative to a single owner who pays for everything, although each pays only part. This contrasts with the asymmetric case, where the paper's Theorem 3 rules out a regular-control equilibrium with smooth payoffs. The result is one of the few explicit equilibria of a nonzero-sum singular-control game with linear costs.
The paper's proofs are on paper only. A formalization requires machine-checked versions of a verification theorem for one-dimensional singular control (the optimality condition), of the analysis of the free-boundary point θ through the fundamental solutions, and of the passage from best response within a restricted class Σi to best response among all adapted controls. The verification steps for a controlled diffusion with jumps, Itô's formula for semimartingales with a nondecreasing part, and the comparison Z≥X with the uncontrolled process are reusable for other singular-control and monotone-follower problems on the platform.
Difficulty
The obvious route, Itô's formula applied to e−rtV(Zt) along an arbitrary deviation, requires the transversality limit E[e−rtV(Zt)]→0 and uniform integrability of the localised terms. V is unbounded below when a=−∞, and the deviation is arbitrary, so this needs the comparison of every controlled state with the uncontrolled process and Assumption 1(ii). The proof also passes through a best-response theorem for the restricted class Σi (Theorem 1), whose conclusion is weaker than a Nash equilibrium; the upgrade to all of Ξi is a separate argument. Finally, the identification of θ uses the boundary behaviour of ψ, φ and the scale density at the ends of I, which Mathlib does not provide for general diffusions.
Formalization scope
Time is [0,∞). The probabilistic data (measure, filtration, Wiener process, the uncontrolled family Xx, ψ, φ) form one structure, and the standing assumptions of §2.1 together with Assumptions 1–3 form one predicate, a hypothesis of every theorem. The Itô integral of the state equation is the published relation EthierKurtz.HasBrownianItoIntegral; Mathlib supplies Brownian motion, stopping times, uniform integrability and Stieltjes measures.
Conventions the formalization commits to:
The initial value x is the pre-jump value Z0−; the discounted cost is a Lebesgue–Stieltjes integral that counts the jump at time 0.
Statements hold for every solution of the relevant state equation; existence and uniqueness of strong solutions (the paper's appeal to SDE theory) is not posed.
Payoffs are compared only for admissible controls: finite expected discounted cost and integrable discounted running profit. Since π is bounded above, every other control has payoff −∞, so no deviation is lost.
Markov perfection is the best-response property from every initial state against the opponent's feedback rule, among all adapted controls.
The fundamental solutions are pinned by the Laplace transforms of hitting times, E[e−rτy]=ψ(x)/ψ(y) for x≤y and φ(x)/φ(y) for y≤x (Borodin and Salminen).
The scale density of Assumption 3(iii) uses a base point x0∈I in place of the printed 0, and the limits of Assumption 3 that may be infinite are limits in the extended reals.
Added, and disclosed: e−rtq(Xtx) is integrable, so that (Rq) of (3) exists, as the paper presupposes.
The optimality condition is stated for V of the form (4) only: for an arbitrary C2 function it is false (an affine V=kx+c with large c satisfies its conditions on I=R but is unbounded).
The class U in Theorem 2 is its state-only version (Lipschitz on I with linear growth).
The threshold θ and coefficient A are not chosen by the formalization: the goal holds for every pair making V twice continuously differentiable on I, and smoothness at θ determines the pair. The rate u is computed from V through the generator, not supplied as a free function, and its nonnegativity is part of the conclusion. A payoff is never a Bochner integral of a non-integrable function, so no junk value can enter an equilibrium comparison.
Contributions welcome: proofs of the milestones, and reusable infrastructure for Itô's formula with nondecreasing càdlàg integrators, comparison theorems for one-dimensional SDEs, and the boundary classification of one-dimensional diffusions.
Selected references
H. Dharma Kwon, Game of Variable Contributions to the Common Good under Uncertainty, arXiv preprint, 2019 (published in Operations Research, 2022). https://arxiv.org/abs/1904.00500v1
L. H. R. Alvarez, Singular stochastic control, linear diffusions, and optimal stopping: A class of solvable problems, SIAM Journal on Control and Optimization 39(6), 1697–1710, 2001. https://doi.org/10.1137/S0363012900367825
L. H. R. Alvarez, On the properties of r-excessive mappings for a class of diffusions, Annals of Applied Probability 13, 1517–1533, 2003. https://doi.org/10.1214/aoap/1069786505
P. C. Lon and M. Zervos, A model for optimally advertising and launching a product, Mathematics of Operations Research 36(2), 363–376, 2011. https://doi.org/10.1287/moor.1110.0493
Minimax-Optimal Policy Learning Under Unobserved Confounding: Minimizing the Worst-Case Empirical Regret Attains the Population Minimax Regret up to an Explicit n^(-1/2) TermResearch Paper
Motivation
Personalized decisions in medicine, marketing and public policy are increasingly learned from observational data: records of who received which treatment and how it turned out. Standard policy-learning methods assume unconfoundedness, that treatment assignment depends only on observed covariates. When an unobserved factor influences both the treatment and the outcome, these methods can return a policy that is worse than the status quo, and nothing in the data reveals it.
Kallus and Zhou (arXiv:1805.08593, published in Management Science 2021 as Minimax-Optimal Policy Learning Under Unobserved Confounding) take the robust route. They fix a bound Γ≥1 on how far the unobserved confounding may distort the propensity scores, in the marginal sensitivity model of Tan (2006), and learn the policy that minimizes the worst-case regret against a baseline policy π0 over every distortion the model allows. The paper proves two guarantees for this policy: it never does worse than the baseline up to vanishing terms (Theorem 2), and it is asymptotically minimax-optimal (Corollary 7). This mission is about the second.
Setting
Data are iid triples (X,T,Y): covariates X∈X, one of m≥2 treatments T∈{0,…,m−1}, and a real outcome Y, read as a loss. The nominal propensity is e~t(x)=P(T=t∣X=x). A policyπ assigns to each x a probability vector π(⋅∣x) over treatments; Π is a class of policies and π0 a baseline.
Inverse propensity weights W(t,x,y) are unknown under confounding. The marginal sensitivity model with parameter Γ confines them to the weight box
The learned policy πˉ^ minimizes Rˉ^π0(⋅;WnΓ) over Π; the population minimax regret is infπ∈ΠRˉπ0(π;WΓ).
The complexity of Π is measured by its VC-major dimension (Definition 1): the largest v such that some v points (tj,xj) have all 2v subsets cut out by superlevel sets {(t,x):π(t∣x)>θ}, π∈Π, θ∈R. The standing assumptions are bounded outcomes ∣Y∣≤B (Assumption 1), overlap with constant ν>0 (Assumption 2), and VC-major dimension of Π at most v (Assumption 3).
Formalization targets
Goal: Corollary 7 (p. 15)
There is a constant K depending only on v such that, with probability at least 1−δ,
Proposition 3 and Corollary 4: each treatment's worst case over the box is attained among monotone weights atΓ(x)+u(y(π(t∣x)−π0(t∣x)))(btΓ(x)−atΓ(x)), u nondecreasing into [0,1], in the population and in every sample.
Lemma 5:suph−supg≤sup(h−g).
Lemma 2: uniform deviation of (π(T∣X)−1/m)Y over Π, at rate 29CK(v+2)log(5/δ)/n.
Lemma 3 and Corollary 15: for a fixed policy, uniform deviation over monotone u, at rate 92πlog(1/p)/n, overall and per treatment.
Lemma 4: the same jointly over Π and u, given the VC-major dimension of the composition class.
Theorem 6 (as proved): with probability at least 1−δ, supπ∈Π∣Rˉ^π0(π;WnΓ)−Rˉπ0(π;WΓ)∣ is at most the goal's rate with 18 in place of 36.
Significance
The result. Corollary 7 says that the confounding-robust policy loses at most O(n−1/2) against the best possible worst-case guarantee in its class. The worst-case regret Rˉπ0(π;WΓ) cannot be identified from data. It is the regret the decision maker must budget for when the confounding is only known to lie within Γ. With Theorem 2, it makes the procedure both safe and asymptotically optimal among safe procedures, for randomized and multi-treatment policy classes of finite VC-major dimension.
Formalizing it. The paper's proof of the dimension-only dependence of K is incomplete. Its Proposition 5, which states that the monotone weight class WˉΓ=⋃π∈ΠWˉΓ(π) has finite VC-major dimension under Assumption 3, is false as printed. A class of monotone policies on [0,1] with VC-major dimension at most 2 gives a WˉΓ that shatters arbitrarily large grids. The printed proof of Theorem 6 runs through Proposition 5 and Lemma 4, so Theorem 6 and Corollary 7 are believed true but are not established by the paper. A formal proof would settle the claim. None of the statements has a machine-checked proof.
Difficulty
The worst case over WΓ is a supremum over an essentially arbitrary bounded function class, which is far too large for a uniform law of large numbers. Proposition 3 reduces it to monotone transformations of the score y(π(t∣x)−π0(t∣x)). This is a linear-fractional program whose population version lives in a function space. The natural next step is to apply a uniform law over the pairs (π,u). That step needs a complexity bound on the composition u∘scoreπ in terms of the dimension of Π alone, and this is exactly where the printed argument fails. A different route is needed. One candidate exploits the Lipschitz dependence of the worst-case Hájek value on the scores. The Hájek ratios also need denominators bounded away from zero uniformly. The paper's decomposition (p. 42) uses E[WI[T=t]]=1, which holds for the true weights only.
Formalization scope
All statements live in the namespace ConfoundRobust.Minimax, with one definitions file Setting.
Data and sampling. An observation is a point of 𝒳 × Fin m × ℝ, with treatments 0-indexed as on the page. The population is a probability measure μ, and a sample is ω : Fin n → Obs 𝒳 m under Measure.pi. En is (1/n)∑i.
Policies and propensity. A policy is π : 𝒳 → Fin m → ℝ, with π x t=π(t∣x), measurable and simplex-valued. The nominal propensity is a measurable version of P(T=t∣X=x), given by its defining integral identity, with ν≤e~≤1.
Assumptions. Assumption 1 is ∣Y∣≤B almost surely. Assumption 2 enters only through its observable consequence e~t(x)≥ν; this is a weaker hypothesis, so the statements are at least as strong. Assumption 3 is the predicate "VC-major dimension at most v" on the functions (t,x)↦π(t∣x).
Weights and suprema. Weight functions in WΓ are measurable. Suprema and infima are real sSup/sInf; the weight sets are nonempty (they contain 1/e~) and the regrets are bounded, so these are true extrema. An empty treatment group gives a Hájek ratio 0/0=0.
Probability statements. "With probability at least 1−δ" means a measurable event of that product measure on which the bound holds for every policy. πˉ^ is any sample minimizer.
Constants. In the goal, Theorem 6, Lemma 2 and Lemma 4 the constant K is quantified as ∀ v, ∃ K, ∀ (everything else). Letting K depend on Π, μ, n or δ would trivialize the bound and is ruled out.
Disclosed readings.
Lemmas 3, 4 and Corollary 15 are posed for p≤1/2: as printed (every p) they are false for p near 1.
The clipped bounds of Corollary 15 and Lemma 4 (PDF p. 41) are reconstructed.
Lemma 2 drops "for n large enough", which its proof does not use.
Lemma 2, Theorem 6 and the goal assume Π is pointwise separable (a countable subclass approximates every policy pointwise), and Lemma 4 assumes the same of its composition class. This is the measurability condition of the empirical-process tools the proofs use. Without it the uniform statements fail (initial segments of ω1 under the countable–cocountable measure).
Theorem 6 uses 18, the constant its proof ends with. The goal's 36 follows from it, but not from the printed Theorem 6.
The printed proof of Theorem 6 passes through Proposition 5, which is false as printed. Do not expect that route to work. Lemma 4 is kept as a statement about the composition class's own dimension.
A complete development needs: linear-fractional optimization over function boxes; Dudley-type entropy bounds for VC-major and VC-hull classes; the Dvoretzky–Kiefer–Wolfowitz inequality with Massart's constant (posed separately on the platform as MassartDKW.Tight.corollary_1); and a uniform bound for Hájek ratios. The empirical-process pieces (VC-major classes, monotone classes) are reusable well beyond this mission. Contributions to any milestone, and alternative routes to Theorem 6 that avoid Proposition 5, are welcome.
Selected references
N. Kallus and A. Zhou, Confounding-Robust Policy Improvement, arXiv:1805.08593v3, 2019; published as Minimax-Optimal Policy Learning Under Unobserved Confounding, Management Science, 2021. https://arxiv.org/abs/1805.08593v3
Z. Tan, A distributional approach for causal inference using propensity scores, Journal of the American Statistical Association 101(476), 2006. https://doi.org/10.1198/016214506000000023
Predicting Many Properties of a Quantum System from Very Few Measurements: Median of Means over Classical Shadows Predicts M Linear Properties to Error ε with Probability 1 − δResearch Paper
Motivation
Characterizing the state ρ of an n-qubit device is a basic task in quantum computing and simulation. Full quantum state tomography reconstructs ρ completely, but it needs a number of measurements exponential in n (Haah et al., 2017; O'Donnell–Wright, 2016). In practice one usually wants only a list of properties, such as the expectation values tr(O1ρ),…,tr(OMρ) of observables (fidelities with target states, entanglement witnesses, local energies). Aaronson's shadow tomography (Aaronson, 2018) showed that M such properties can be learned from a number of copies polylogarithmic in M, but with collective measurements across many copies that are far beyond present hardware.
Huang, Kueng and Preskill (arXiv:2002.08953, Nature Physics 16, 2020) gave a procedure that uses only single-copy randomized measurements: rotate ρ by a random unitary from a fixed ensemble, measure in the computational basis, and store the outcome classically. From these records, called classical shadows, a median-of-means estimator predicts all M properties at once. Their Theorem 1 bounds the number of measurements by O(logM) times a state-independent norm of the observables. The method has since become a standard tool for benchmarking and for estimating many observables in experiments.
Setting
Fix n qubits. Computational-basis labels are bit strings b∈{0,1}n, and operators are 2n×2n complex matrices. A state is a positive semidefinite matrix of trace one.
A measurement primitiveU (Definition 1) is a finite ensemble of unitaries Ui drawn with probabilities wi. One snapshot of ρ draws i with probability wi, applies ρ↦UiρUi†, and measures, obtaining b with probability ⟨b∣UiρUi†∣b⟩ (S1). The primitive is tomographically complete if any two distinct states σ=ρ give different outcome probabilities ⟨b∣UiσUi†∣b⟩=⟨b∣UiρUi†∣b⟩ for some b and some i of positive weight.
The measurement channel (S2) is the linear map
M(X)=i∑wib∑⟨b∣UiXUi†∣b⟩Ui†∣b⟩⟨b∣Ui,
and the classical shadow of the outcome (i,b) is ρ^=M−1(Ui†∣b⟩⟨b∣Ui) (S3). For an observable O the single-shadow prediction is o^=tr(Oρ^). The shadow norm (S7) is
The median-of-means prediction (S12) takes NK independent snapshots, splits them into K batches of N, averages tr(Oρ^j) within each batch, and outputs the median o^(N,K) of the K batch means.
Formalization targets
Goal: Theorem 1 (p. 13)
For Hermitian O1,…,OM, ϵ,δ∈(0,1] and integers
K≥2logδ2M,N≥ϵ2341≤i≤MmaxOi,0shadow2,N≥1,
the median-of-means predictions from NK independent classical shadows satisfy
Pr[∣o^i(N,K)−tr(Oiρ)∣≤ϵ for all 1≤i≤M]≥1−δ.
Milestones
Tomographic completeness makes M invertible (p. 11).
tr(ρ^)=1 and E[ρ^]=ρ (p. 11).
E[o^]=tr(Oρ), (S5).
M−1 is self-adjoint for the trace pairing (p. 12).
Lemma 1: Var[o^]≤∥O0∥shadow2, (S6).
The median-of-means guarantee: for N≥34σ2/ϵ2, Pr[∣μ^(N,K)−EX∣≥ϵ]≤2e−K/2 (p. 13).
The per-observable bound Pr[∣o^i(N,K)−tr(Oiρ)∣≥ϵ]≤δ/M (p. 13).
The constants 34, 2 and 2M/δ are the paper's own. Remark 1 of the paper says they are chosen for simplicity, not tightness.
Significance
Theorem 1 gives a sample complexity NK=O(log(M)maxi∥Oi,0∥shadow2/ϵ2). It grows only logarithmically in the number of properties and does not depend on the dimension 2n except through the shadow norm. For random global Clifford measurements the shadow norm is at most 3tr(O02). For random single-qubit Pauli measurements it is at most 4k∥O∥∞2 for an observable acting on k qubits (Propositions 1–3 of the paper). These instantiations are what make the protocol useful for fidelity estimation and for local observables. The theorem also separates two independent ingredients, a state-independent variance bound and a generic median-of-means concentration bound, and each is reusable.
The result is proved in the paper. As far as the platform's corpus shows, it has not been machine-checked. This mission produces a formal statement of the protocol and its guarantee, and its definitions (measurement channel, classical shadow, shadow norm) can be reused for the Clifford and Pauli instantiations and for the paper's quadratic-function estimator (Theorem 4).
Difficulty
The linear-algebra layer is elementary but has to be set up carefully. The channel must be shown invertible from completeness on states alone. Its inverse must be shown trace-preserving and self-adjoint, and the variance has to be rewritten as an expectation of squared diagonal entries.
The probabilistic layer is where the naive argument fails. Chebyshev's inequality on a single empirical mean gives a sample size proportional to 1/δ, not log(1/δ). Applying Hoeffding's inequality to the number of failing batches gives exponent 2K(1/2−1/34)2≈0.44K. That is weaker than K/2, so it does not reach the stated constant 2e−K/2. The bound needs the relative-entropy form of the Chernoff bound for a binomial tail, which Mathlib does not provide in this form. The median also has to be related to a count of deviating batch means without sign or tie errors.
Formalization scope
Model. Operators are Matrix (Fin n → Fin 2) (Fin n → Fin 2) ℂ. States are the published WildeQIT.IsDensityOperator.
Ensemble. The ensemble is finite: a published WildeQIT.FinDist of weights and a family in Matrix.unitaryGroup. The paper's ensemble is general, but both of its examples are finite with uniform weights. "U∈U" in completeness means positive weight.
Channel and inverse.M is a genuine C-linear endomorphism and M−1 is Ring.inverse of it. Every statement using M−1 assumes tomographic completeness.
Real parts. Predictions and targets are real parts of traces, which are real for Hermitian observables.
Shadow norm. The maximum in (S7) is a real supremum over the nonempty, compact set of states, of a function bounded above.
Probability. The NK snapshots are modelled exactly. Probabilities are finite sums over all K×N outcome arrays of the product of single-snapshot probabilities wi⟨b∣UiρUi†∣b⟩, and batch k is row k of the array. The median-of-means milestone is stated for general real random variables in Mathlib's measure-theoretic language (independence, identical distribution, L2).
Median. The median is the lower median.
Parameters.N and K are natural numbers at least the printed real values. N≥1, ϵ>0, δ>0 and M≥1 are added, because at ϵ=0 or δ=0 the printed values are undefined.
Theorem 1 must not be restated about an abstract estimator with a variance bound; that is milestone 6. The goal is about the median of means of classical shadows of a concrete primitive, with the shadow norm of (S7). Tomographic completeness is satisfiable non-trivially: the three single-qubit Pauli bases form a complete primitive, so the goal's hypotheses are not vacuous.
Needed infrastructure: elementary facts about the trace pairing and diagonal entries of UXU†; finite product measures, or direct finite-sum manipulations; Chebyshev's inequality; a binomial Chernoff bound in relative-entropy form; and order-statistic facts about the median. The binomial tail bound and the median-of-means lemma are independent of quantum information and reusable elsewhere. Contributions to any milestone are welcome. The Clifford and Pauli shadow-norm computations are natural follow-up work on top of these definitions.
Selected references
H.-Y. Huang, R. Kueng, J. Preskill, Predicting many properties of a quantum system from very few measurements, Nature Physics 16, 1050–1057 (2020); arXiv:2002.08953v2. https://arxiv.org/abs/2002.08953
J. Haah, A. W. Harrow, Z. Ji, X. Wu, N. Yu, Sample-optimal tomography of quantum states, IEEE Trans. Inf. Theory 63 (2017). https://arxiv.org/abs/1508.01797
M. R. Jerrum, L. G. Valiant, V. V. Vazirani, Random generation of combinatorial structures from a uniform distribution, Theoret. Comput. Sci. 43 (1986), 169–188. https://doi.org/10.1016/0304-3975(86)90174-X
A. S. Nemirovsky, D. B. Yudin, Problem Complexity and Method Efficiency in Optimization, Wiley (1983). (Jerrum–Valiant–Vazirani and Nemirovsky–Yudin are the paper's references [41, 55] for median of means.)
NLTS Hamiltonians from good quantum codesResearch Paper
Motivation
NLTS Hamiltonians from good quantum codes establishes the selected mathematical guarantees below. The mission turns those original paper results into precise Lean proof targets.
Setting
Fix positive clustering constants c1,c2,δ0, a bound s on the number of nonzero entries in each row and column of the binary check matrices, and a positive constant B bounding the total number of checks by Bn. There exist a>0 and an integer n0≥2, uniform over the code, state, energy density and preparing circuit, with the following property.
For n≥n0, let Hx∈F2mx×n and Hz∈F2mz×n be check matrices with mx,mz>0, sparsity at most s, and mx+mz≤Bn. Put Cx=kerHx and Cz=kerHz, and assume the CSS condition Cx⊥⊆Cz. Write k for the encoded dimension, and assume k≥2. Assume Property 1: for every 0≤δ<δ0 and every y∈F2n, ∣Hzy∣≤δmz implies ∣y∣Cx⊥≤c1δn or ∣y∣Cx⊥≥c2n, and the same implication holds with x,z interchanged. Here ∣y∣S=minv∈S∣y+v∣, and ∣⋅∣ is Hamming weight. Define the unnormalized check Hamiltonian H=∑w row of Hx(I−Xw)/2+∑w row of Hz(I−Zw)/2, counting all supplied rows.
For every density matrix ρ and every 0≤ϵ<Eeff(C),
tr(Hρ)≤ϵn⟹depth(U)≥alogn
for every preparing circuit U. A preparing circuit consists of layers of disjoint arbitrary one- and two-qubit unitary gates on n+w qubits, for any w≥0, starts in ∣0⟩⊗(n+w), and has reduced output ρ on the first n qubits after discarding the other qubits. Depth is its number of layers; log is the natural logarithm.
Moreover, for every fixed r,q>0, there exist e∗>0 and n1≥2, uniform over all code instances, such that whenever n≥n1, k≥rn and min{mx,mz}≥qn, one has
k≥2,Eeff(C)≥e∗.
Consequently the allowed energy density can be chosen as a positive constant for the constant-rate, linear-distance qLDPC families satisfying Property 1 considered in the paper. The distance parameter does not enter the displayed cutoff; the paper's footnote 3 notes that Property 1 at δ=0 already implies linear distance.
This gives a depth obstruction for every low-energy density matrix, including mixed states, rather than only ground states.
Source correction. The paper prints only ϵ<E(C). Its proof uses Property 1 at 2ϵ1 with ϵ1=200nϵ/min(mx,mz), so we also require 400nϵ<δ0min(mx,mz). This is exactly the missing domain condition; it is automatic if c1≥1. No normalization c1≥1 is imposed, and the frozen printed cutoff is retained as the first entry of the minimum.
Correction to the printed source.
Printed: Theorem 5 permits eps < E(C), where E(C) is the printed energyCutoff, for arbitrary positive c1.
Corrected: Use eps < min(E(C), delta0min(mx,mz)/(400n)); give the eventual positive floor for this same minimum. No other hypothesis, constant, or conclusion is changed.
Reason: The source proof sets epsilon1=200neps/min(mx,mz) on p.3 and applies Property 1 at 2epsilon1 on p.4. Its domain is delta<delta0, but the printed cutoff only supplies 2c1epsilon1<delta0. The additional strict upper bound is exactly 2epsilon1<delta0. It is redundant for c1>=1 and retains every c1<1 instance satisfying the actual proof-domain condition.
Target
Theorem 5 (Formal statement of the NLTS theorem), §2
Fix positive clustering constants c1,c2,δ0, a bound s on the number of nonzero entries in each row and column of the binary check matrices, and a positive constant B bounding the total number of checks by Bn. There exist a>0 and an integer n0≥2, uniform over the code, state, energy density and preparing circuit, with the following property.
For n≥n0, let Hx∈F2mx×n and Hz∈F2mz×n be check matrices with mx,mz>0, sparsity at most s, and mx+mz≤Bn. Put Cx=kerHx and Cz=kerHz, and assume the CSS condition Cx⊥⊆Cz. Write k for the encoded dimension, and assume k≥2. Assume Property 1: for every 0≤δ<δ0 and every y∈F2n, ∣Hzy∣≤δmz implies ∣y∣Cx⊥≤c1δn or ∣y∣Cx⊥≥c2n, and the same implication holds with x,z interchanged. Here ∣y∣S=minv∈S∣y+v∣, and ∣⋅∣ is Hamming weight. Define the unnormalized check Hamiltonian H=∑w row of Hx(I−Xw)/2+∑w row of Hz(I−Zw)/2, counting all supplied rows.
For every density matrix ρ and every 0≤ϵ<Eeff(C),
tr(Hρ)≤ϵn⟹depth(U)≥alogn
for every preparing circuit U. A preparing circuit consists of layers of disjoint arbitrary one- and two-qubit unitary gates on n+w qubits, for any w≥0, starts in ∣0⟩⊗(n+w), and has reduced output ρ on the first n qubits after discarding the other qubits. Depth is its number of layers; log is the natural logarithm.
Moreover, for every fixed r,q>0, there exist e∗>0 and n1≥2, uniform over all code instances, such that whenever n≥n1, k≥rn and min{mx,mz}≥qn, one has
k≥2,Eeff(C)≥e∗.
Consequently the allowed energy density can be chosen as a positive constant for the constant-rate, linear-distance qLDPC families satisfying Property 1 considered in the paper. The distance parameter does not enter the displayed cutoff; the paper's footnote 3 notes that Property 1 at δ=0 already implies linear distance.
This gives a depth obstruction for every low-energy density matrix, including mixed states, rather than only ground states.
Source correction. The paper prints only ϵ<E(C). Its proof uses Property 1 at 2ϵ1 with ϵ1=200nϵ/min(mx,mz), so we also require 400nϵ<δ0min(mx,mz). This is exactly the missing domain condition; it is automatic if c1≥1. No normalization c1≥1 is imposed, and the frozen printed cutoff is retained as the first entry of the minimum.
Correction to the printed source.
Printed: Theorem 5 permits eps < E(C), where E(C) is the printed energyCutoff, for arbitrary positive c1.
Corrected: Use eps < min(E(C), delta0min(mx,mz)/(400n)); give the eventual positive floor for this same minimum. No other hypothesis, constant, or conclusion is changed.
Reason: The source proof sets epsilon1=200neps/min(mx,mz) on p.3 and applies Property 1 at 2epsilon1 on p.4. Its domain is delta<delta0, but the printed cutoff only supplies 2c1epsilon1<delta0. The additional strict upper bound is exactly 2epsilon1<delta0. It is redundant for c1>=1 and retains every c1<1 instance satisfying the actual proof-domain condition.
Significance
These are known results from the cited paper. The current contribution consists of source-reviewed, compiler-checked theorem statements and mathematical definitions. Formal proofs remain open. Only the listed selected results are in scope; no whole-paper completeness is claimed.
Difficulty
A complete solution must establish every quantified conclusion under the exact stated hypotheses, including error, probability and resource conditions wherever present. The source assumptions and conventions are recorded in the attached statements and independent readbacks.
Formalization scope
Definitions preserve the existing reviewed Lean interfaces. Every target proof is currently sorry, so each target is Open. Printed-source corrections, where necessary, are labeled with the printed claim, the change and its reason.
There is a deterministically explicit family of binary CSS check descriptions on n qubits, and constants ϵ,c,a,b>0, a locality bound s, and a threshold n0≥1, all independent of n. Explicit means that a deterministic polynomial-time function on unary n outputs the actual row counts and both binary check matrices. For every n≥1, both check matrices have row and column weights at most s, their total row count mn satisfies an≤mn≤bn, and their associated projector terms hi=(I−Xwi)/2 or (I−Zwi)/2 act on at most s qubits, commute, and annihilate a common unit vector. Let Hn=∑ihi be this unnormalized Hamiltonian. Every circuit P of parallel one- and two-qubit unitaries on disjoint wires, starting from zero qubits and allowing arbitrary clean ancillas and discarded qubits, has an output density operator ρP satisfying, for n≥n0,
Tr(HnρP)<ϵn⟹depth(P)≥clnn.
This is the existence of an NLTS Hamiltonian family. Quantifying over every preparation circuit gives the lower bound for mixed as well as pure low-energy states.
Formalization Note The deterministic interpretation of explicit and the CSS witness implementation are recorded for independent moderation; the theorem does not assume an NLTS guarantee on its family.
Section 2: support of both measurement distributions
Let HX,HZ be a binary CSS check presentation with n,mx,mz>0, and let ρ be any density operator with Tr(Hρ)≤ϵn, where ϵ≥0 and H is the sum of the actual check projectors. Write Dz,Dx for its standard- and product-Hadamard-basis distributions and Gbδ for the words violating at most a δ fraction of the actual checks on side b. For every ϵ1>0,
Dz(Gzϵ1)≥1−ϵ1mzϵn,Dx(Gxϵ1)≥1−ϵ1mxϵn.
If ϵ=0, both Gz0 and Gx0 have probability one. If ϵ>0, choosing ϵ1=200nϵ/min(mx,mz) gives both masses at least 199/200.
These bounds place low-energy measurement outcomes in the approximate classical codes. The energy identity and probability bounds are conclusions of the source argument, not hypotheses of the formal statement.
Section 2: equivalence classes and separated clusters
Let CX=kerHX and CZ=kerHZ be a binary CSS pair satisfying Property 1 with positive constants c1,c2,δ0. Assume ϵ1≥0, 2ϵ1<δ0, and 4c1ϵ1<c2. On Gzϵ1 define x∼zy when ∣x+y∣CX⊥≤2c1ϵ1n, with the corresponding definition on Gxϵ1 modulo CZ⊥. Here ∣v∣S=mins∈S∣v+s∣. For either basis this is an equivalence relation. Its distinct computed classes are nonempty, pairwise disjoint, and cover the corresponding approximate code. Distinct classes satisfy
x∈B,y∈B′,B=B′⟹∣x−y∣≥c2n.
This constructs the separated clusters needed for the depth argument.
Formalization Note The strict domain follows Property 1 on page 2; the page-4 prose writes a non-strict endpoint. The classes are computed from the relation rather than supplied with a separation guarantee.
Correction to the printed source.
Printed: Page 4 says 'assuming 2 epsilon1 <= delta0', but Property 1 on page 2 is stated only for delta < delta0.
Corrected: Use 0 <= epsilon1 and the strict inequality 2 epsilon1 < delta0 in the clustering regime.
Reason: The proof invokes Property 1 at delta = 2 epsilon1; the strict domain is necessary. This is a standing-domain clarification pending independent moderation.
Lemma 2: one basis has no highly concentrated cluster
Consider a binary CSS presentation with n,mx,mz>0 and logical dimension k≥1, satisfying Property 1 with c1,c2,δ0>0. Let ρ be any density operator. Assume ϵ1≥0, 2ϵ1<δ0, 4c1ϵ1<c2, and
2c1ϵ1≤(4nk−1)2.
For the computed clusters in the preceding result, either every standard-basis cluster has Dz-mass strictly less than 99/100, or every Hadamard-basis cluster has Dx-mass strictly less than 99/100.
This uncertainty consequence selects a basis in which no single cluster captures almost all probability. It requires no low-energy hypothesis once the clusters are defined.
Formalization Note The subtraction k−1 takes place in the reals. The positive-k reading comes from the section's positive-rate setting and is flagged for independent moderation.
Correction to the printed source.
Printed: The displayed finite lemma does not explicitly say k>=1, while Section 2 is introduced for constant-rate code families.
Corrected: State the finite lemma for k>=1 and n,mx,mz>0, keeping ((k-1)/(4n))^2 with real subtraction.
Reason: The proof uses a positive logical dimension; at k=0 the zero-radius clusters can concentrate in both bases. The finite positive-rate reading requires independent moderation.
Fact 3: uncertainty for subsets in conjugate bases
Let ρ be any n-qubit density operator, and let Dz,Dx be its standard- and product-Hadamard-basis measurement distributions. For all sets S,T⊆{0,1}n,
Dx(T)≤21−Dz(S)+2n∣S∣∣T∣.
This bounds concentration in conjugate bases and applies to mixed states. No positive-mass hypothesis on S is required; empty sets and n=0 are allowed.
Section 2: positive-mass separated unions
Let a binary CSS presentation on n>0 qubits satisfy Property 1 with c1,c2,δ0>0, and assume ϵ1≥0, 2ϵ1<δ0, and 4c1ϵ1<c2. Fix either basis and any density operator ρ. Suppose its measurement distribution D gives the approximate-code set mass at least 199/200, while each computed class has mass less than 99/100. Then there are disjoint subcollections M,M′ of the computed classes whose unions S,S′ are disjoint and satisfy
D(S)≥4001,D(S′)≥4001,∣x−y∣≥c2n(x∈S,y∈S′).
This extracts two positive-mass separated sets in whichever basis Lemma 2 selects. Both probability constants and the cross-union separation are retained.
Fact 4: depth lower bound for well-spread distributions
Let P be a circuit of depth t consisting of layers of disjoint one- and two-qubit unitaries, starting in the all-zero product state. Allow any number of clean ancillas and measure only its n>0 output qubits in the standard basis, giving distribution D. Let 0<μ<1, and suppose D(S),D(T)≥μ. For any u≥0 such that every x∈S,y∈T has Hamming distance at least u,
u2≤400nln(1/μ)23t.
For positive set distance this is equivalently t≥31log2(dist(S,T)2/(400nln(1/μ))).
This is the quantitative circuit-depth obstruction for well-spread output distributions. The exponential formulation also covers overlapping sets without taking the logarithm of zero; the original circuit's workspace is unrestricted.
Lemma 8: expansion implies clustering
Let A∈F2m×n and V⊆F2n. Let δ>0 and 0<μ<ν. Suppose every x∈V with ∣x∣≤νn satisfies ∣Ax∣≥(δ/μ)∣x∣. Then
x∈V,∣Ax∣≤δn⟹∣x∣≤μnor∣x∣≥νn.
This transfers small-set expansion to the clustering dichotomy. Both the size cutoff and syndrome cutoff use the number of columns n, even for rectangular matrices.
Formalization Note Positivity is taken from the proof and the paragraph immediately following Lemma 8; the unsupported zero-δ extension is excluded.
Correction to the printed source.
Printed: The displayed Lemma 8 omits its parameter domain; its proof uses delta>0 and mu>0 and the following paragraph takes mu<nu.
Corrected: State Lemma 8 on the meaningful positive regime delta>0 and 0<mu<nu, as pinned by the plan and proof.
Reason: For delta=0, an expansion hypothesis with beta=0 supplies no information and the claimed nontrivial dichotomy can fail. Independent moderation must confirm the intended domain.
Let F be a family of binary CSS check presentations. Fix positive αx,βx,αz,βz,δˉ,a,b and n0≥1. For all n≥n0, suppose an≤mx(n),mz(n)≤bn. For every 0≤δ<δˉ, define Vzδ as the words in Gzδ of minimum Hamming weight within their coset modulo CX⊥, and define Vxδ analogously modulo CZ⊥. Assume HX is (αx,βx)-expanding on Vxδ and HZ is (αz,βz)-expanding on Vzδ. Then there are constants c1,c2,δ0>0, independent of n, such that for all n≥n0 and 0≤δ<δ0,
y∈Gzδ⟹∣y∣CX⊥≤c1δnor∣y∣CX⊥≥c2n,
and the same dichotomy holds on the X side modulo CZ⊥.
This derives Property 1 from uniform expansion of minimal representatives, without assuming the desired dichotomy. The check-density factors converting δmside to a threshold relative to n are retained.
These additional statements have current no-edit source review. Their formal proofs remain open; deferred results are excluded.
Wigner quasiprobability: nonnegativity and Hudson's theoremTextbook
Motivation
Classical statistical mechanics describes a particle by a probability density on phase space; quantum mechanics replaces this by a wavefunction, and the uncertainty principle forbids any honest probability interpretation of simultaneous position and momentum. In 1932, Eugene Wigner introduced the transform that has become the standard way of placing quantum states on classical phase space anyway: a quasiprobability distribution W(x,p), defined by an explicit integral formula, whose marginals reproduce the Born position and momentum distributions of the wavefunction (Wigner 1932). Unlike a probability density, W generally takes negative values, and this negativity is the standard marker of genuinely quantum behavior — of interference — in quantum optics, quantum information, and the study of the quantum-to-classical transition.
The central question about the probabilistic character of W is: which quantum states have an everywhere nonnegative Wigner function? Hudson's theorem (1974) gives the complete answer in one dimension: only Gaussian wave packets (Hudson 1974). This mission formalizes the definition, the basic properties, and the positivity characterization of the Wigner function, following the source text, the Wikipedia article Wigner quasiprobability distribution (oldid 1376057304).
Timeline of the key results:
1932 — E. Wigner introduces the distribution to study quantum corrections to thermodynamic equilibrium (Physical Review 40, 749).
1946 — H. J. Groenewold computes the Wigner function of the m-th Fock (number) state, this mission's key worked example (Physica 12, 405).
1949 — J. E. Moyal recognizes W as the moment-generating functional of quantum mechanics in phase space (Proc. Camb. Phil. Soc. 45, 99).
1974 — R. L. Hudson proves that Gaussian pure states are the only pure states with nonnegative Wigner function in one dimension (Rep. Math. Phys. 6, 249).
1983 — F. Soto and P. Claverie extend the characterization to arbitrary dimension via the Segal–Bargmann transform (J. Math. Phys. 24, 97).
Setting
A wavefunction is a function ψ:R→C, and ℏ>0 denotes the reduced Planck constant. The Wigner function of ψ at the phase-space point (x,p)∈R2 is
W(x,p)=πℏ1∫−∞∞ψ(x+y)ψ(x−y)e2ipy/ℏdy,
a Lebesgue (Bochner) integral over y. The Lean development factors this into three definitions mirroring the formula: the integrand (Wigner.wignerIntegrand), the integral (Wigner.wignerIntegral), and the Wigner function itself (Wigner.wigner), defined as the real part of the integral divided by πℏ; that the discarded imaginary part vanishes is itself milestone 1. The Fock states of the harmonic oscillator enter through the physicists' Hermite polynomials Hm and the ordinary Laguerre polynomials Lm, defined by their standard three-term recurrences, and the normalized m-th Hermite function
um(x)=(2mm!)−1/2π−1/4Hm(x)e−x2/2.
Formalization targets
Goal — Hudson's theorem
∫R∣ψ∣2=1 and Wψ≥0 on R2⟹∃a,b,c∈C,Rea>0,ψ=e−ax2+bx+ca.e.
This is exactly the source's statement of Hudson's theorem; nothing is left unfixed. The mission is complete when the goal is proved.
Milestones
The milestone list is the source's property list, ordered as the attack path: (1) reality of the defining integral; (2) the position marginal ∫Wdp=∣ψ(x)∣2; (3) the total integral ∬W=∫∣ψ∣2; (4) Groenewold's closed formula for the Fock states; (5) a concrete negativity witness, the first excited Fock state at the origin; (6) the pointwise bound ∣W∣≤2/h; (7) the converse of Hudson's theorem, that Gaussian pure states have nonnegative Wigner function.
Significance
The result itself. Hudson's theorem draws the exact boundary between the classical and the quantum in phase space: Gaussian states — coherent and squeezed states — are precisely the pure states that can be described by an honest phase-space probability density. This is used in quantum optics as a nonclassicality criterion, in tomography and state-reconstruction theory, and in the study of the quantum-to-classical transition: negativity of W is shielded by the uncertainty principle (negative regions cannot extend beyond a few ℏ) and disappears in the classical limit, as the source emphasizes.
What formalizing it adds. The Wigner transform and this circle of ideas are absent from current Mathlib, so the mission builds the foundations itself: the integral transform and its basic identities (reality, marginals, normalization, boundedness), the Hermite/Laguerre recurrences and Groenewold's Fock-state formula, and finally Hudson's theorem with its converse. All of these are reusable well beyond the mission — for quantum optics, for the Wigner–Ville time–frequency distribution used in signal analysis, and for the wider phase-space formulation of quantum mechanics (Moyal brackets, star products). These are known results: this is formalization of known mathematics, not new mathematics, and no machine-checked proof of Hudson's theorem exists in Mathlib today.
Difficulty
The naive approach to Hudson's theorem is doomed from the start, and solvers should know it: since the Husimi Q function (a Gaussian smoothing of W) is positive for every quantum state, no smoothing or positivity-of-convolution argument can extract any information from W≥0. The known route (Soto–Claverie) is: W≥0 everywhere forces the Husimi function strictly positive everywhere; the Husimi function is the squared modulus of the Segal–Bargmann transform F of ψ, so F is a zero-free entire function; hence F=eg for some entire g, and square-integrability of F against a Gaussian measure forces g to have at most quadratic growth, hence to be a quadratic polynomial; inverting the Segal–Bargmann transform then exhibits ψ as a Gaussian. The hard Lean steps are the growth-to-polynomial (Liouville-type) step, the inversion of the transform, and, for the milestone properties, the Fubini and Fourier-inversion machinery. Mathlib contains the raw harmonic analysis (Fourier inversion for Schwartz functions) but none of the assembled phase-space layer.
Formalization scope
Integrals are Bochner integrals with respect to Lebesgue measure on R; statements carry the integrability or regularity hypotheses under which the source's displayed integrals exist.
ℏ is an explicit real parameter, and every statement assumes 0<ℏ (the standing convention of the source). The Fock-state section uses the source's dimensionless units, i.e. ℏ=1.
The Wigner function is defined as the real part of the defining integral divided by πℏ; milestone 1 proves the imaginary part vanishes, so this is the same function as the source's.
Pure states are encoded by the normalization ∫∣ψ∣2=1.
The conclusion of Hudson's theorem holds only up to modification of ψ on a Lebesgue-null set (W is insensitive to null-set changes of ψ), so the goal states almost-everywhere equality; this is the faithful reading of "the wavefunction must have the form".
The Fock states are normalized: the source's displayed um=π−1/4Hme−x2/2 omits the factor (2mm!)−1/2 that its Fock-state formula requires; the definition includes it, since Groenewold's formula is stated for the normalized number states ∣m⟩=(a†)m/m!∣0⟩.
There is no trivializing encoding: Hudson's hypothesis ranges over all normalized wavefunctions, and the conclusion is not weakened beyond the almost-everywhere reading.
Infrastructure contributed and reusable beyond this mission: hermiteP, laguerreL, fockState, wignerIntegrand, wignerIntegral, wigner. Welcome additional contributions: the momentum marginal, the state-overlap formula (Property 6), Galilei covariance (Property 4), and harmonic-oscillator time evolution, as further theorems in the same namespace.
J. E. Moyal, Quantum mechanics as a statistical theory, Mathematical Proceedings of the Cambridge Philosophical Society 45(1), 99–124, 1949. https://doi.org/10.1017/S0305004100000487
Let Φ:M2(C)→M2(C) be complex-linear, completely positive and trace preserving. Let ∣β⟩=(∣00⟩+∣11⟩)/2 and ΓΦ=(id⊗Φ)(∣β⟩⟨β∣). Then
Φ is entanglement breaking⟺ΓΦ≤21I4.
The inequality is the positive-semidefinite operator order. This provides an operator test on the normalized Choi state.
Theorem 3
Let t,λ∈R3, and let Φt,λ be the trace-preserving canonical Bloch map v↦t+diag(λ)v. If it is entanglement breaking, then
j=1∑3∣λj∣≤1.
There is no unitality assumption. This is a necessary condition for the general translated channel.
Theorem 4
Let Φ be a completely positive trace-preserving qubit map with Φ(I)=I. Suppose independent unitary changes of input and output bases put it in canonical form with parameters t,λ∈R3. Then
Φ is entanglement breaking⟺j=1∑3∣λj∣≤1.
This makes the necessary condition sufficient in the unital case.
Amplitude damping counterexample
For 0≤α≤1, let Φα have translation (0,0,1−α2) and diagonal Bloch parameters (α,α,α2). These maps are completely positive and trace preserving, and
Φα is entanglement breaking⟺α=0.
Moreover, for every η>0 there is α∈(0,1] with 2α+α2<η and Φα not entanglement breaking. Thus small diagonal-parameter sum does not suffice for translated channels.
Linear channels and equation 9
For t∈R3 and a∈R, let Φ=Φt,(0,0,a). Then
Φ is CPT⟺t12+t22+(∣t3∣+∣a∣)2≤1.
Whenever this holds, Φ is CQ and equals X↦∑b=01Tr(∣b⟩⟨b∣X)Rb, where R0=(I+(t+ae3)⋅σ)/2 and R1=(I+(t−ae3)⋅σ)/2. This includes a=0, when the image is a point.
QC channels are unitarily equivalent to CQ channels
Let Φ be a QC qubit channel: it measures a POVM and prepares one of two orthogonal pure states. There are qubit unitaries U,V and a CQ map Ψ such that
Φ=AdV∘Ψ∘AdU,AdU(X)=UXU†.
A CQ map measures rank-one orthogonal input projectors and prepares density matrices. The assertion relates these two channel subclasses without identifying their definitions.
Extreme POVM need not give an extreme channel
Set w1=(1,0,0), w2=(−1/2,0,3/2), w3=(−1/2,0,−3/2), and Ek=(I+wk⋅σ)/3. Set Rk=(I+σk)/2, with σ1=X,σ2=Y,σ3=Z, and define
Φ(X)=k=1∑3Tr(EkX)Rk.
The effects form an extreme three-outcome POVM, each Rk is a density matrix, and Φ is EBT, neither CQ nor a point channel, and not extreme in the EBT set. Thus extreme measurement data need not produce an extreme channel.
Proposition 7
Let Φ=Φt,λ be a CPT qubit map in the fixed canonical basis, and let T(X)=XT. Then T∘Φ∘T is completely positive, and changing any individual translation coordinate tj to −tj preserves complete positivity. In addition,
CP(Φ∘T)⟺[∀j,CP(Φt,λ[j↦−λj])]⟺CP(T∘Φ).
Here CP denotes complete positivity. These are all three numbered assertions, together with the any-coordinate translation invariance explained immediately before the proposition.
Planar channels are entanglement breaking
Let Φt,λ be a CPT qubit map in canonical form. Then
(∃j∈{1,2,3}:λj=0)⟹Φt,λ is entanglement breaking.
The planar regime includes degenerate lines and points. This is a sufficient condition only.
Pure point channels are extreme
Let t∈R3 with ∥t∥2≤1, and set R=(I+t⋅σ)/2. The canonical point map is CPT and CQ, R is a density matrix, and on every operator X it satisfies
Φt,0(X)=Tr(X)R.
If ∥t∥2=1, there is a unit vector ψ with R=∣ψ⟩⟨ψ∣, and the map is extreme in both the CPT and EBT sets. This identifies the pure-output endpoint of the point-channel family.
Truly planar channels do not touch the Bloch sphere
Let Φt,λ be a CPT canonical qubit map. If exactly one coordinate of λ is zero, then
∥t+diag(λ)v∥2<1(∥v∥2≤1).
Also, when λ3=0, (∣λ1∣+∣λ2∣)2+∥t∥22≤1. If both λ1,λ2 are nonzero and equality holds, then t1=t2=0. Thus a truly planar image stays strictly inside the Bloch sphere.
Theorem 8
In a fixed diagonal basis, let U={λ∈R3:Φ0,λ is CPT and EBT}, and V={±e1,±e2,±e3}. Then
U=convR(V),Ext(U)=V.
The theorem determines both the entire parameter octahedron and exactly its six extreme points.
Unital EBT channels are the convex hull of unital CQ maps
In the real vector space of complex-linear qubit maps, let Eu be the unital trace-preserving entanglement-breaking maps, and Cu the unital CQ maps. Then
Eu=convR(Cu).
The channel sets allow arbitrary input and output bases. This is a statement in map space, not only in one diagonal parameter slice.
Theorem 9
Fix any translation t∈R3, and define Λt={λ:Φt,λ is completely positive}. Then
Λt∩(−Λt)={λ:Φt,λ is EBT}.
Inversion negates λ and leaves t fixed. Infeasible translations are permitted; both sides can be empty.
Theorem 10
Let E be the trace-preserving entanglement-breaking complex-linear qubit maps, and C the CQ maps. Then
Ext(E)⊆C,E=convR(C).
These are both conclusions of the main theorem. It does not assert that every CQ map is extreme. The proof gap concerning intersecting planes on p.15 is recorded separately, without becoming a theorem hypothesis.
Remark a planar same-basis convex decomposition
Suppose Φt,λ is EBT and λj=0 for a chosen coordinate j. It belongs to the real convex hull of canonical CQ maps Φt′,λ′ satisfying λj′=0 and having at most one nonzero coordinate of λ′. Moreover, if Ξt3,0 denotes the CPT parameter slice fixing only t3 and λ3=0, then
Ξt3,0=convR{(t,λ)∈Ξt3,0:Φt,λ is CQ}.
This preserves the canonical coordinate plane, including point-channel endpoints; the other translation coordinates may vary.
Equation 29 axial-translation EBT criterion
For a CPT canonical qubit map with t=(0,0,t3),
Φt,λ is EBT⟺(∣λ1∣+∣λ2∣)2≤(1−∣λ3∣)2−t32.
For fixed t3,l3 with ∣l3∣≤1 and q=(1−∣l3∣)2−t32≥0, its EBT slice is exactly {(v1,v2):∣v1∣+∣v2∣≤q}. This includes the filled square and its collapsed endpoints; a negative radicand is not totalized into an extra radius-zero slice.
Remark c parallelogram of same-basis decompositions
Fix t∈R3 and l3∈R. Let P consist of canonical EBT maps with these fixed parameters and either λ1=0 or λ2=0. Let Cdiag be all CQ maps canonical in the same basis. Then
convR(P)⊆convR(Cdiag).
This gives the same-basis decomposition on the illustrated parallelogram. It neither covers the entire slice nor requires CQ summands to retain l3.
Remark d equal singular magnitudes
Let Φ=Φt,λ be EBT, with ∣λj∣=a>0 for all three coordinates. There are input/output unitaries U,V and s∈{a,−a} such that, writing B(Ψ)=AdV∘Ψ∘AdU,
B(Φ)=Φ(∥t∥2,0,0),(s,s,s)∈convR(A).
Here A consists of canonical CQ maps for which tj=λj=0 for some j∈{2,3}. Also Φ∈convR(B−1(A)). Thus summands are canonical in the rotated basis and need not be canonical in the original one.
Formalization scope
Only these reviewed original-paper statements are in scope. Definitions have reviewed bodies; target proofs remain open (sorry). Reviewed deferrals are excluded. No theorem proof or whole-paper formal proof is claimed.
Product Decomposition of Periodic Functions in Quantum Signal Processing: Every SU(2)-Valued Laurent Polynomial of Degree n in E_n Is Uniquely a Product of n Primitive Matrices and a UnitaryResearch Paper
Motivation
Quantum signal processing transforms the eigenphases of a unitary operator by a sequence of controlled operations and single qubit rotations. To know which transformations such a sequence can realize, one must characterize the matrix functions produced by the sequence. Haah's product decomposition theorem gives such a characterization without the restriction on the projectors used in some earlier quantum signal processing constructions. The result concerns exact functions; the numerical synthesis algorithm and its precision analysis are separate parts of the paper.
Setting
Let t be a nonzero complex number. A rank one orthogonal projectorP on C2 is a Hermitian matrix with P2=P and trace one. Its complementary projector is I−P. The associated primitive matrix is
EP(t)=tP+t−1(I−P).
For ∣t∣=1, this matrix is unitary with determinant one. Given n such projectors, indexed in their multiplication order, and a constant matrix E0∈SU(2), the product of interest is
F(t)=E0EP1(t)EP2(t)⋯EPn(t).
The leftmost primitive factor is EP1, and the product has exactly n primitive factors. At n=0 the product of primitive factors is the identity and F is the constant matrix E0. Although the physical interpretation uses points on the unit circle, the formula defines F at every nonzero complex t. These conventions follow §2, equations (2) and (3).
A matrix Laurent polynomial of width n is written
F(t)=j=−n∑nCjtj,
where each Cj is a complex 2×2 matrix and all coefficients outside [−n,n] vanish. The set Pn consists of these polynomials for which F(η)∈SU(2) whenever ∣η∣=1. A member has degree n when Cn=0 or C−n=0. The subset En consists of members whose nonzero exponents have the parity of n, so C−n+k=0 for every odd integer k. The notation F†(t) conjugate transposes each coefficient while leaving its exponent unchanged: F†(t)=∑jCj†tj. These are the conventions of Definition 1 and equations (4) and (5).
Formalization targets
The goal is all of Theorem 2. First, every product of n primitive matrices belongs to En. Conversely, every F∈En of degree n has a unique ordered tuple (E0,P1,…,Pn) satisfying the product formula, with E0∈SU(2) and each Pj a rank one orthogonal projector. For n>0, the final projector is determined by the leading coefficient:
Pn=Tr(Cn†Cn)Cn†Cn.
The formula includes the assertion that its denominator is nonzero under the theorem's hypotheses; that is a consequence to establish, rather than an extra assumption. The n=0 case retains the product characterization and uniqueness, while having no final projector. The milestones state the source's primitive matrix observation, polynomial identities, leading coefficient conditions, degree reduction, and the relation between the two extreme coefficients.
Significance
The theorem identifies the Laurent polynomials that admit the specified quantum signal processing product and says that the factors are uniquely recoverable from the matrix function when the number of factors is fixed. This supports the exact synthesis question motivating the paper: whether a desired SU(2) valued function on the unit circle is achievable by the primitive operations. The formula for Pn also identifies the last factor directly from a coefficient of the target function. Haah's §4 algorithm uses the decomposition, but its numerical guarantees are outside this mission.
The paper proves Theorem 2. The formalization target is a machine checked version of that known theorem and its source stated intermediate claims. The definitions of matrix Laurent polynomial evaluation, parity restricted coefficient families, and primitive projector factors can also support later work on polynomial matrix functions. The mission does not claim that these statements have already been proved in Lean.
Difficulty
Membership in SU(2) is given only on the unit circle, while coefficient identities needed for the decomposition concern Laurent polynomials throughout C∖{0}. A pointwise unitary identity on the circle must therefore be promoted to an identity of Laurent polynomials before its extreme coefficients can be compared. The main structural difficulty is that a matrix product may lose its highest and lowest powers even when its factors are nonzero. Identifying a projector that lowers the degree must control both extremes at once and preserve the parity and SU(2) conditions. Uniqueness also depends on fixing exactly n factors: extra complementary factors can cancel.
Formalization scope
Lean represents a Laurent polynomial by a coefficient family C:Z→Mat2(C) with explicit width support, and evaluates it by a finite sum of integer powers. Equalities of Laurent polynomials are stated for every nonzero complex t. The positive integer n is represented as a natural number, with integer casts for exponents. The special unitary group is Mathlib's SU(2); the paper's rank one orthogonal projectors are represented by Hermitian idempotence together with trace one. This trace condition is equivalent to rank one for these 2×2 projectors. An ordered function on Fin(n) supplies precisely n factors. The last factor has index n−1 in Lean.
The set En retains both the width and SU(2) conditions of Pn and adds the parity condition. It is defined from coefficients, not from the existence of a product decomposition. Allowing rank zero or rank two projectors, dropping unitarity, or allowing a variable factor count would change the theorem. The source's remark that any orthogonal projector yields a member of E1 is read with the rank one condition of equation (2): the literal unrestricted remark fails for 0 and I. The singular leading coefficient milestones require n>0, and the last projector formula has the same domain.
The development needs facts about 2×2 determinants, conjugate transpose, projectors, finite matrix products, and uniqueness of finite Laurent expansions. The Laurent polynomial infrastructure and rank one matrix facts are reusable beyond this theorem. Contributions proving the source stated milestones, their faithful supporting algebraic lemmas, and the main product characterization are in scope. The parity refinement of §3.1, Lemma 4, and the §4 synthesis algorithm with its running time claim are outside this mission.
Selected references
J. Haah, Product Decomposition of Periodic Functions in Quantum Signal Processing, Quantum 3, 190 (2019); formalization source: arXiv:1806.10236v4, 2020, pp. 3–5.
Quantum Singular Value Transformation and Beyond: Exponential Improvements for Quantum Matrix Arithmetics: Singular Value Transformation by Real Polynomials via Alternating Phase ModulationResearch Paper
Motivation
Many quantum algorithms act on a matrix A that is not itself unitary: a Hamiltonian to be simulated, a linear system to be inverted, a projector to be amplified. The standard way to hand such a matrix to a quantum computer is to embed it as a block of a unitary U, so that A=ΠUΠ for two orthogonal projectors Π,Π. The question is then which functions f(A) can be produced from a few uses of U and U†.
Gilyén, Su, Low and Wiebe (arXiv:1806.01838, STOC 2019) answered it for polynomials: interleaving U and U† with phase rotations about Π and Π applies a polynomial to every singular value of A. This quantum singular value transformation (QSVT) recovers Hamiltonian simulation, amplitude amplification, fixed-point search, quantum walks and matrix inversion as instances of one identity, and is now the common language of these algorithms (see the survey of Martyn, Rossi, Tan and Chuang, PRX Quantum 2021).
Timeline.
2016: Low, Yoder and Chuang characterize the single-qubit sequences eiφ0σz∏jW(x)eiφjσz, which they call quantum signal processing (QSP) (PRX 2016).
2016–2017: Low and Chuang lift QSP to Hermitian blocks by qubitization (Quantum 2019; PRL 2017).
2018: this paper replaces eigenvalues by singular values, so that arbitrary (even rectangular) blocks are allowed, and removes the earlier constraint Pℜ(0)=0 for even degree.
Setting
Let HU be a finite-dimensional complex Hilbert space, U a unitary on it, and Π,Π orthogonal projectors. The projected unitary is A=ΠUΠ. Put d=rankΠ, d~=rankΠ, dmin=min(d,d~). A singular value decomposition of A consists of orthonormal bases (ψi)i≤d of imgΠ and (ψ~i)i≤d~ of imgΠ and numbers ς1≥⋯≥ςdmin≥0 with
A=i=1∑dminςi∣ψ~i⟩⟨ψi∣.
For f:R→C, the singular value transformation is f(SV)(A)=∑i≤dminf(ςi)∣ψ~i⟩⟨ψi∣ when f is odd and f(SV)(A)=∑i≤df(ςi)∣ψi⟩⟨ψi∣ when f is even, with ςi:=0 for i>dmin.
The phase operator of a projector is eiφ(2Π−I)=eiφΠ+e−iφ(I−Π). For Φ=(φ1,…,φn)∈Rn, the alternating phase modulation sequenceUΦ is the product, from j=1 on the left to j=n on the right, of eiφj(2Π−I)U when n−j is even and of eiφj(2Π−I)U† when n−j is odd. On the single-qubit side, R(x)=(x1−x21−x2−x) for x∈[−1,1], and Φ is said to realize P∈C[x] if P(x) is the top-left entry of ∏j=1neiφjσzR(x) for every x∈[−1,1].
In Lean these are ProjUnitarySVD U Proj Projt (the decomposition data), svtOdd/svtEven, phaseRefl, altPhaseSeq, reflR and qspRefl, in the namespace QSVT19.Main.
Formalization targets
Goal: Corollary 18 (singular value transformation by real polynomials)
Let n≥1 and let Pℜ∈R[x] have degree n, parity nmod2, and ∣Pℜ(x)∣≤1 on [−1,1]. Then there is Φ∈Rn such that for every singular value decomposition of ΠUΠ
with Π′=Π for odd n and Π′=Π for even n. The only constraints on Pℜ are those that are necessary; no degree-dependent constant appears.
Milestones
Lemma 6, Theorem 5, Theorem 4 (first claim): sum-of-squares factorization and the characterization of which polynomials, and which real parts, arise from QSP.
Corollary 8: every P meeting four explicit conditions is realized by some Φ∈Rd in the reflection form; Corollary 10: every Pℜ as in the goal is ℜ[P] for such a P.
Definition 11 (existence of the decomposition), Table 2 (orthonormality of the spanning vectors of Definition 12), Lemma 14 (the block decomposition of U into one- and two-dimensional pieces).
Theorem 17: if Φ realizes P, then P(SV)(ΠUΠ)=ΠUΦΠ for odd n and ΠUΦΠ for even n.
Lemma 9 (optional): explicit phases realizing the Chebyshev polynomial Td.
Significance
The result. Corollary 18 says that every real polynomial of degree n, of definite parity and bounded by 1 on [−1,1], can be applied to the singular values of ΠUΠ with n uses of U and U† and one ancilla qubit. Every application in §3.4–§5 of the paper (fixed-point amplitude amplification, singular vector transformation, the pseudoinverse, Hamiltonian simulation, Gibbs sampling) is this corollary applied to a suitable approximating polynomial, so a formal proof makes all of them conditional only on polynomial approximation bounds.
Formalizing it. The result is proved in the paper; none of the statements of this mission has a machine-checked proof on the platform. A formal development would provide reusable infrastructure absent from Mathlib: singular value decompositions over sub-spaces, the two-dimensional invariant-subspace decomposition of a projected unitary (a quantitative form of Jordan's lemma), and the polynomial characterization of single-qubit QSP. Theorem 3 of the paper, the W(x) form of the QSP characterization, is already posed on the platform as QAlgorithms.Childs.qsp_characterization; Corollary 8 depends on it.
Difficulty
The obvious argument for Theorem 17 is "restrict UΦ to each two-dimensional subspace and read off the QSP product". Making that rigorous means proving that the subspaces of Definition 12 are pairwise orthogonal, invariant under every factor of UΦ, and that U acts on each pair of bases by exactly the reflection R(ςi), including the degenerate blocks with ςi∈{0,1}, the one-dimensional blocks HiR,HiL where UΦ contributes a phase e±iφ0, and the complement H⊥. The polynomial side is harder still: Lemma 6 factors a nonnegative even polynomial through its multiset of complex roots, grouping roots in (0,1), in [1,∞), on the imaginary axis and in the open quadrant, while tracking degree and parity. Neither step reduces to a library lemma, and an argument that works only for Hermitian A, or only when Π=Π, proves a different theorem.
Formalization scope
HU is ι → ℂ for a Fintype ι with the inner product star v ⬝ᵥ w; operators are Matrix ι ι ℂ; "unitary" is U ∈ Matrix.unitaryGroup ι ℂ; "orthogonal projector" is IsStarProjection. The projectors are named Proj, Projt because Π is Lean notation. Indices are 0-based.
A singular value decomposition is the structure ProjUnitarySVD U Proj Projt of data: orthonormal families with Π=∑i∣ψi⟩⟨ψi∣, Π=∑i∣ψ~i⟩⟨ψ~i∣, nonnegative non-increasing ς, and (17). The goal is posed for every such decomposition, with Φ chosen first; the Definition 11 milestone shows that one exists. ςi≤1 is a consequence and is not assumed.
f(SV) is the sum over the decomposition (Definition 16), never the right-hand side of (32) or (33); a formalization that defines it by an operator expression trivializes the goal.
eiφ(2Π−I) is written eiφΠ+e−iφ(I−Π), equal to the matrix exponential for an idempotent Π. The ancilla of (33) is the first tensor factor.
eiφσz and the parity predicate are the published definitions of QAlgorithms_Defs_Hamiltonian. "deg ≤ k − 1" is degree < k and "parity k−1" is parity k+1, avoiding natural-number subtraction.
Printed slips that affect statements:
Corollary 10 omits ℜ[P]=Pℜ; it is restored, as used in the proof of Corollary 18.
Corollary 8 and Corollary 18 are false for degree 0 (a constant eiα, respectively a constant c=1); both are posed with d≥1, n≥1.
Theorem 17's "P and Φ as in Corollary 8" is read as (14), the only property its proof uses.
Lemma 14's (25)–(28) index the HL blocks by [d]∖[r]; this is read as [d~]∖[r], as in (24).
Contributions welcome: proofs of any milestone, general lemmas on SVDs of compressions ΠUΠ, and a proof of the platform's open QSP characterization.
Selected references
A. Gilyén, Y. Su, G. H. Low, N. Wiebe, Quantum singular value transformation and beyond: exponential improvements for quantum matrix arithmetics, arXiv:1806.01838v1, 2018; STOC 2019. https://arxiv.org/abs/1806.01838
G. H. Low, T. J. Yoder, I. L. Chuang, Methodology of resonant equiangular composite quantum gates, Physical Review X 6, 041067, 2016. https://doi.org/10.1103/PhysRevX.6.041067
Persuading Customers to Buy Early: The Value of Personalized Information Provisioning: The Personalized Buy/Wait/No Mechanism Is Optimal and Strictly Beats Public SignalingResearch Paper
Motivation
A seller who sells over two periods, at a full price now and a discount later, loses revenue to strategic customers who wait for the discount. Waiting is costly to the customer only through availability risk: if inventory is scarce, the product may be gone by the second period. Inventory is usually known to the firm and not to its customers, so the firm can try to shape how scarce the product looks. Drakopoulos, Jain and Randhawa (SSRN 3191629, version of April 17, 2020) study this as a problem of Bayesian persuasion (Kamenica and Gentzkow, AER 2011) with many receivers: the firm commits to prices and to a signaling policy before it learns its inventory, and then sends each customer a recommendation.
Their two main results separate public from private information. Theorem 1 (p. 15): with a public signal, seen by all customers, the firm earns no more than with prices alone. Theorem 2 (p. 17): with personalized signals, which may depend on each customer's valuation and be randomized independently across customers, a simple Buy/Wait/No mechanism is optimal and strictly beats every public mechanism when the low inventory level is small. Personalized recommendations are the kind of information online retailers and travel platforms send (low-stock alerts, buy-now prompts).
Setting
There is a unit mass of customers with valuations distributed on V=[0,b] with density f and complementary distribution Fˉ(v)=∫vbf; the hazard rate f/Fˉ is non-decreasing. The firm's inventory Q is high, qH=1 (enough for everyone), with probability P(H)∈(0,1), and low, qL∈(0,1], with probability P(L)=1−P(H). The monopoly price is pm=argmaxppFˉ(p), and Fˉ−1(q) is the valuation at which a mass q of customers lies above.
A straightforward mechanism fixes prices p1,p2≥0 and, after seeing Q, draws a latent index k from a finite set with a state-dependent law; given (Q,k), each v-customer is independently told Buy (in period 1), Wait (buy in period 2) or No. A customer who requests the product in a period with s units left and demand of mass D from the others gets it with probability avail(s,D): 0 if s≤0, 1 if D≤s, s/D otherwise (random rationing). The mechanism is obedient if every customer (for almost every v) weakly prefers to follow her recommendation, given the availabilities she can expect conditional on it: (Buy-1)–(No), p. 40. The firm's revenue is the expected value of rationed sales p1min{D1,Q}+p2min{D2,(Q−D1)+}. A mechanism is public when the recommendations depend only on (k,v) and obedience holds for each realized k.
The personalized mechanism (PS1)–(PS2), p. 17, sets p1=Fˉ−1(qL). In state L it tells customers with v≥p1 to buy and the rest No; in state H it tells v≥p1 to buy with probability
σH(v)=min{P(H)P(L)p1−p2v−p1,1}(13)
and to wait otherwise, tells p2≤v<p1 to wait, and v<p2No. The price p2 solves Fˉ(p2)−p2f(p2)=Fˉ(vˉ) (14) with vˉ=sup{v:σH(v)<1} (15).
Formalization targets
Goal: Theorem 2
If qL≤Fˉ(pm), then (14) has a root p2∈(0,p1), and for every such root the mechanism (PS2) is obedient, its revenue Rprivate is at least that of every obedient mechanism at any prices, and
Lemma 2 (p. 35): the L-state profit ϕ(t,p) of a public threshold t satisfies ϕ(t,p)≤max{(Fˉ−1(qL)−p)qL,0}, and for qL<Fˉ(p) its unique maximizer is Fˉ−1(qL).
Theorem 1 (c) (p. 15): Rϕ∗ is the largest revenue of a public mechanism.
Lemma 3 (p. 44): at fixed prices p1>p2, the (t∗,p1,p2)-mechanism built from an optimal threshold t∗ of the one-parameter program (A.16) is obedient, earns its value Dˉ(p1,p2), and is optimal among obedient mechanisms with those prices.
Lemma 4 (p. 44, corrected reading): if Fˉ(p2)>qL, then Dˉ(p1,p2)≤Dˉ(Fˉ−1(qL),p2).
The result says that the value of information in this market lies entirely in its personalization: public information is worthless (Theorem 1), and private recommendations strictly raise revenue whenever low inventory is scarce relative to the monopoly market. The optimal mechanism is explicit and has a simple form, with the low-type firm selling out exactly in period 1 and the high-type firm mimicking it as far as incentive compatibility allows; in the paper's uniform example (qL=P(H)=1/2) revenue rises from 1/4 to 13/48. The proof technique, a relaxation to a linear program over signal profiles followed by a dual reduction to a single threshold, is a template for many-receiver persuasion problems with congestion.
The theorems are proved in the paper; none of them has a machine-checked proof. The mission produces a formal model of strategic customers under rationing and personalized signaling, and the formal versions of the five steps. Two printed steps need correction to be true and are stated in corrected form (Lemma 4's tie case, the location of the root of (14)); formal proofs settle that the corrections suffice.
Difficulty
The obvious argument is a pointwise one: for each valuation, the Buy incentive constraint caps the high-state Buy probability given the low-state one, so set every cap tight. It fails because the availabilities in the constraints depend on the aggregate masses that the recommendations themselves create, so the constraints for different valuations are coupled; the paper resorts to an infinite-dimensional linear program and asserts strong duality without proof (p. 41). Lemma 4 is the algebraically heavy step (pp. 46–50): an envelope argument over the threshold combined with the hazard-rate monotonicity. Theorem 2 then needs a one-dimensional optimization in p2, the existence of a root of (14) (a fixed-point equation, since vˉ depends on p2), and the case p2>p1, which the paper never treats.
Formalization scope
The paper's game (Definitions 1–2: a continuum of customers, Bayesian Nash equilibrium, sender-preferred selection) is re-expressed as the class its proofs reduce to by the revelation principle (p. 39, footnote 5): obedient straightforward mechanisms with aggregate masses. Revenue optimality over this class is what Theorem 2 asserts after that reduction. The latent index is finite (Fin n, for every n); it carries the correlation between customers' signals that public signals need. Obedience is required for almost every valuation, as in Definition 1.
Committed conventions: V=[0,b], b>0, f continuous on [0,b] and positive on [0,b); 0<P(H)<1, 0<qL≤1; qH=1 and customer mass 1. Availability is avail above; at zero demand it is 1 with stock left and 0 without. (14) is written Fˉ(p2)−p2f(p2)=Fˉ(vˉ). In (PS2)(a) "v>p1" is taken as v≥p1 (a null set). (B) is read with min{⋅,1}. The program (A.16) keeps the no-Buy mechanism (Fˉ(t)=0), which the printed version drops. "The optimal mechanism" is stated as "an optimal mechanism" (ties exist). Lemmas 3–4 take p2<p1, implicit in the paper's divisions; the goal ranges over all prices p1,p2≥0, including p2>p1 (p. 8). Milestone (16) is stated at the prices of (PS1): p1=Fˉ−1(qL) and p2∈[0,p1) a root of (14).
Trivializing formalizations are ruled out: revenue counts rationed sales, not p1D1+p2D2; obedience is the printed conditional constraint multiplied by the probability of the recommendation, never divided by a possibly zero probability; availability never uses Lean's x/0=0; the goal carries the existence of a root of (14), so "for every root" is not vacuous; the mechanism class contains public and correlated mechanisms; and Rprivate is the revenue of the mechanism (PS2), tied to the formula (16) by a milestone.
A complete development needs: the density measure on [0,b] and its tail function Fˉ, monotonicity and continuity of Fˉ and its inverse, the quasi-concavity of pFˉ(p) under a non-decreasing hazard rate, and bounds for the rationing function. The hazard-rate and Fˉ−1 lemmas are reusable for other pricing models. Contributions of helper lemmas, of a primal proof of Lemma 3 avoiding the duality argument, and of the case p2>p1 are welcome.
Selected references
K. Drakopoulos, S. Jain, R. S. Randhawa, Persuading Customers to Buy Early: The Value of Personalized Information Provisioning, SSRN working paper 3191629, version of April 17, 2020. https://ssrn.com/abstract=3191629
D. Bergemann, S. Morris, Information Design: A Unified Perspective, Journal of Economic Literature 57(1):44–95, 2019. https://doi.org/10.1257/jel.20181489
Addictive Platforms 1: Monopoly versus Duopoly — Competition Raises Addictiveness and Hurts the Consumer When Attention Is Scarce, and Lowers Addictiveness When It Is AbundantResearch Paper
Motivation
Digital services are paid for with attention: an ad-financed platform earns more the longer its users stay. Policy reports on digital markets (Crémer, de Montjoye and Schweitzer, 2019; the U.K. Digital Competition Expert Panel, 2019; Scott Morton et al., 2019) worry that such firms degrade their services with features that capture attention at the user's expense, and they usually recommend more competition. Ichihashi and Kim, Addictive Platforms (Bank of Canada Staff Working Paper 2022-16; Management Science, 2023), ask whether competition actually helps. Their answer depends on how scarce attention is. When attention is scarce, platforms compete by making their services more addictive rather than better, and a monopoly serves the consumer better than a duopoly. This mission formalizes that comparison, the paper's headline result, together with the two equilibrium characterizations it rests on.
Setting
There are K platforms and a single consumer. Each platform k chooses an addictivenessdk≥0. The consumer observes d=(d1,…,dK), joins a set J of platforms and allocates attention ak≥0 to each k∈J. She receives 0 if J=∅, and otherwise
k∈J∑u(ak,dk)−C(k∈J∑ak),
subject to the attention cap∑k∈Jak≤A, where A∈R+∪{∞}. Platform k earns the attention ak if the consumer joins it and 0 otherwise.
The service utilityu(a,d) is twice differentiable, strictly increasing and concave in attention a, and strictly decreasing in d. The marginal utility u1(a,d)=∂u/∂a is strictly increasing in d, so a more addictive service is worse to join but harder to stop using. There is some d at which no amount of attention makes the service worth joining, and u(0,0)≥0. The attention costC:R+→R+ is strictly increasing, convex and twice differentiable, with lima→0[u1(a,0)−C′(a)]>0. These are Assumptions 1 and 2 of the paper.
Platforms choose addictiveness simultaneously, and then the consumer chooses J and a. The solution concept is pure-strategy subgame perfect equilibrium: the consumer chooses optimally at every profile d, and no platform gains by changing its own dk. A monopolist (K=1) breaks ties in the consumer's favour. Consumer surplus is the consumer's payoff in the equilibrium outcome.
The monopoly benchmark uses three quantities. The first is A(d)=argmaxA≥0u(A,d)−C(A), the consumer's attention on one platform when there is no cap. The second is dA(A)=min{d∈[0,∞]:A(d)≥A}, the lowest addictiveness at which she exhausts her attention. The third is dP(A), the addictiveness at which her best capped payoff from one platform falls to zero. Under competition, AK(d)=argmaxA∈[0,A]Ku(A/K,d)−C(A) is the total attention she spends on K platforms of common addictiveness d.
Formalization targets
Goal: Proposition 3 (monopoly versus duopoly)
Let AM be the least cap from which on monopoly consumer surplus is zero. For every monopoly equilibrium with addictiveness dM and every duopoly equilibrium with (d1,d2):
The proposition is silent about intermediate caps, and so is the goal.
Milestones
Proposition 1 (monopoly). The monopolist sets min{dA(A),dP(A)}, which is increasing in A. The consumer's payoff rises on [0,A(0)], falls on [A(0),AM] and is zero from some AM>A(0) on.
Lemmas B.1–B.4. These are comparative statics of the consumer's problem: the tangent intercept u(x)−xu′(x) is increasing, the value of spreading attention over y platforms is concave in y, and the incremental gain from one more platform falls with the addictiveness and the number of the platforms already joined.
Step 1 of the proof of Proposition 2. For K≥2 there is a unique d∗>0 solving the indifference condition (3) between joining K and K−1 platforms.
Proposition 2 (competition). For K≥2, every platform chooses d∗ in the unique equilibrium.
Significance
Proposition 3 separates two effects of competition. Business stealing pushes duopolists to raise addictiveness, because a more addictive service captures a larger share of a fixed attention budget. Competition for participation restrains them, because a consumer can drop a platform and spend her time on its rival. When attention is scarce the first effect dominates and the consumer is strictly worse off under duopoly. When it is abundant the monopolist extracts the whole surplus and duopoly cannot be worse. The result is the basis of the paper's policy discussion: digital curfews and the comparison between the equilibrium and the joint-profit outcome (Propositions 4 and 6). Propositions 1 and 2 are used throughout the paper.
The paper proves these results on paper. No machine-checked proof exists. A formal proof would check several delicate points. One is the role of the monopolist's tie-break, without which every d∈[0,dP] is a monopoly equilibrium when attention is scarce. Another is that the uniqueness argument for the duopoly needs a consumer who best-responds off the equilibrium path. The appendix also asserts a few claims without argument or slightly too strongly, recorded below. Verifying them is part of the work.
Difficulty
The obvious argument for Proposition 3 compares equilibrium payoffs through the closed forms of Figure 3, and closed forms do not exist in general. Every object is defined by an optimization problem whose maximizer has no formula. The value functions are concave but need not be differentiable, so the envelope-theorem steps of Appendix B cannot be applied as written; they have to be replaced by monotone comparative statics or one-sided derivatives. Proposition 2's uniqueness is the hardest part. It rules out asymmetric equilibria by a deviation argument that controls which platforms the consumer joins after a small change in one platform's addictiveness (Lemma B.4). It also needs the consumer's whole strategy, not only her choice on the path.
Formalization scope
Platforms are Fin K; the paper's platform 1 is the index 0. Addictiveness and attention are real numbers restricted to [0,∞), and only the values of u and C there matter. The cap A is an extended nonnegative real (ℝ≥0∞), so A=∞ is a genuine case and "A≥AM" includes it. Derivatives on [0,∞) are one-sided at 0.
The definition file states the paper's Assumptions 1–2 as printed and adds disclosed pins:
u(⋅,d) is strictly concave. The page says "concave", but Lemma B.1 needs strict concavity and Lemma B.2 applies Lemma B.1 to u(⋅,d).
C(0)=0. Without it dP need not exist.
The uncapped consumer problem has an optimal attention profile for every profile and every nonempty J. The page (p. 6) assumes this for A=∞, and A(d) presupposes it.
A>0 where the result needs it: Proposition 1's characterization, Step 1, Proposition 2 and the first half of Proposition 3. At A=0 every profile is an equilibrium.
"Increasing" and "decreasing" are weak unless the page says "strictly". The one exception is "higher addictiveness" in the first half of Proposition 3, which is stated strictly because the proof gives dM=0<d∗.
Three statements differ from the printed text. Lemma B.2's "V∗ is strictly concave" is false on the paper's own Figure 3 instance and is stated as concavity. The strict clause of Lemma B.4 fails when C is linear and is not posed. The derivative claims of Lemmas B.2–B.4 are stated as decreasing differences of value functions.
Several trivializing formalizations are excluded:
An equilibrium requires the consumer to best-respond at every profile, not only at the equilibrium one.
The tie-break is part of the monopoly equilibrium.
Equilibria exist (Propositions 1 and 2 assert existence), and a sanity file checks the standing hypotheses on u(a,d)=1−e−(a−d), C(a)=a/2.
No value function is a real supremum; maxima enter through maximizers.
Not posed: Proposition 4 and later results, and the side models of Remarks 2–3. A complete development needs maximizers of concave problems on intervals and simplices, monotone comparative statics (Topkis, Milgrom–Shannon), and envelope arguments for concave value functions. These pieces are reusable well beyond this paper, and contributions of them as separate lemmas are welcome.
J. Crémer, Y.-A. de Montjoye and H. Schweitzer, Competition Policy for the Digital Era, report for the European Commission, 2019. https://doi.org/10.2763/407537