December 03, 2023
Let \(\mathcal{M}\) denote the class of randomised monotone functions on \(\mathbb{R}\) with values in \([0,1]\), and let \(U_{\mathcal{M}}\colon \mathbb{R}_+\to \mathbb{R}_+\) be the minimal function for which \[\mathbb{P}\left\{ \sqrt{\eta_f}\, \sup_{t\in\mathbb{R}} \left| f_Z(t) - \mathbb{E}\left[f_Z(t)\right] \right| \ge \varepsilon\sqrt{U_{\mathcal{M}}(\eta_f)} \right\} \le 2\mathrm{e}^{-2\varepsilon^2}\] holds for every member \(f_Z\) of \(\mathcal{M}\) of finite effective sample size \(\eta_f\) and every positive \(\varepsilon\). We prove that for every \(x> 1\), \[\left| \sqrt{U_{\mathcal{M}}(x)} - \sqrt{\log_4 x} \right| \le 2 \min\!\left\{ 1,\, \frac{2 \ln(\mathrm{e} + \ln x)}{\sqrt{\ln x}} \right\}\,.\] The optimal scale \(\sqrt{U_{\mathcal{M}}(x)}\) is sharply tied, uniformly at finite sample sizes, to \(\frac{1}{\sqrt{2\ln 2}}\sqrt{\ln x}\).
Massart [1] establishes that for every empirical distribution function \(F_n \colon \mathbb{R} \to [0,1]\) based on \(n\) independent and identically distributed observations, \[\mathbb{P}\left\{ \sqrt{n} \sup_{t\in\mathbb{R}} \left|F_n(t) - \mathbb{E}\left[F_n(t)\right]\right| \ge \varepsilon \right\} \le 2\mathrm{e}^{-2\varepsilon^2} \,, \qquad \text{for all } \varepsilon > 0 \,.\] The normalisation \(\sqrt{n}\) is the same as the optimal distribution-free normalisation implied by Hoeffding’s inequality [2] for a single threshold: \[\mathbb{P}\left\{ \sqrt{n} \left|F_n(t) - \mathbb{E}\left[F_n(t)\right]\right| \ge \varepsilon \right\} \le 2\mathrm{e}^{-2\varepsilon^2} \,, \qquad \text{for all } \varepsilon > 0 \,.\]
The Massart bound is not robust to mild departures from the structure of the randomised function: \[F_n(t) = \frac{1}{n} \sum_{i=1}^n \mathbf{1}_{\{Z_i \le t\}} \,, \qquad \mathbb{P}\{Z_i \le t\}= \mathbb{P}\{Z_j \le t\}\,,\qquad \text{for each } t\in \mathbb{R}\,.\] In this paper we determine the optimal sub-Gaussian normalisation for the class of all randomised \([0,1]\)-valued monotone functions.
Let \[f \colon \mathcal{Z}_1 \times \cdots \times \mathcal{Z}_n \times \Theta \to \mathbb{R}\] be a bounded function on the product \(\mathcal{Z} = \prod_{i=1}^n \mathcal{Z}_i\) of the measurable spaces \(\mathcal{Z}_i\) and an index set \(\Theta\). For each \(z \in \mathcal{Z}\) and \(\theta \in \Theta\), define the sections \(f_\theta(z) := f(z, \theta)\) and \(f_z(\theta) := f(z, \theta)\). We suppose throughout that, for every \(\theta \in \Theta\), the section \(f_\theta \colon \mathcal{Z} \to \mathbb{R}\) is measurable.
For a random vector \(Z = (Z_1, \dots, Z_n)\) with independent components \(Z_i \in \mathcal{Z}_i\), the mapping \(f_Z \colon \Theta \to \mathbb{R}\) is a randomised function, and the pair \((f, Z)\) is a randomisation.
Define \(\eta_f \in (0, \infty]\) by \[\eta_f := \left( \sum_{i=1}^n c_i(f)^2 \right)^{-1} \,, \qquad \text{where} \qquad c_i(f) := \sup_{(\theta, z_i^\ast, z) \in \Theta \times \mathcal{Z}_i \times \mathcal{Z}} \left| f_\theta(z_i^\ast, z_{-i}) - f_\theta(z) \right| \,.\] If \(\eta_f\) is finite, then the McDiarmid–Hoeffding inequality [3] gives, for a single index, \[\label{eq:macd} \mathbb{P}\left\{ \sqrt{\eta_f} \left| f_Z(\theta) - \mathbb{E}\left[f_Z(\theta)\right] \right| \ge \varepsilon \right\} \le 2\mathrm{e}^{-2\varepsilon^2} \,, \qquad \text{for all } \varepsilon > 0\,.\tag{1}\] For the empirical distribution \(F_n\) we have \(\eta_{F_n} = n\). We call \(\eta_f\) the effective sample size and take \(\sqrt{\eta_f}\) as the baseline distribution-free normalisation. When \(\eta_f=\infty\), the function \(f_Z\) is deterministic for any \(Z\).
Let \(\mathsf{M}\) be the space of all deterministic non-decreasing functions on \(\mathbb{R}\) with values in \([0,1]\). Write \(\mathcal{M}\) for the class of all randomisations \((f, Z)\), taken over all finite product measurable spaces, that satisfy \[f_z \in \mathsf{M}\,, \qquad \text{for each } z \in \mathcal{Z}\,.\] We also, for convenience, write \(f_Z \in \mathcal{M}\) and call \(f_Z \colon \mathbb{R} \to [0,1]\) a randomised monotone function.
Let \(U_{\mathcal{M}} \colon \mathbb{R}_+ \to \mathbb{R}_+\) be the infimum of all functions \(U \colon \mathbb{R}_+ \to \mathbb{R}_+\) that satisfy \[\label{eq:admis} \mathbb{P}\left\{ \sqrt{\eta_f} \sup_{t \in \mathbb{R}} \left| f_Z(t) - \mathbb{E}\left[f_Z(t)\right] \right| \ge \varepsilon \sqrt{U(\eta_f)} \right\} \le 2\mathrm{e}^{-2\varepsilon^2} \,, \qquad \text{for all } \varepsilon > 0\,,\tag{2}\] for every \(f_Z\in \mathcal{M}\) satisfying \(\eta_f<\infty\). The function \(U_{\mathcal{M}}\) is well-defined, monotone, non-zero except at zero, and itself satisfies 2 (Proposition 1). Thus, \(U_\mathcal{M}\) provides the optimal distribution-free normalisation for the class: \[\sqrt{\eta_f / U_{\mathcal{M}}(\eta_f)}, \qquad \text{0<\eta_f<\infty}\,.\] The main concern of this paper is to uniformly and sharply bound the optimal adjustment function \(\sqrt{U_\mathcal{M}}\). We prove the following result.
Theorem 1 (Main theorem). For every \(x > 1\), \[\left| \sqrt{U_{\mathcal{M}}(x)} - \sqrt{\log_4 x} \right| \le 2 \min\!\left\{ 1, \, \frac{2 \ln(\mathrm{e} + \ln x)}{\sqrt{\ln x}} \right\} \,.\]
The normalisation may be taken to be \[\frac{\sqrt{\eta_f}}{\sqrt{\log_4 \eta_f} + 2} \,, \qquad \text{for all } 1 < \eta_f < \infty \,,\] and rearranging the bound of Theorem 1, we get a tight bound on logarithmic growth \[\left| \sqrt{\frac{U_{\mathcal{M}}(x)}{\ln x}} - \frac{1}{\sqrt{2\ln 2}} \right| \le 2 \min\!\left\{ \frac{1}{\sqrt{\ln x }} , \, \frac{2 \ln(\mathrm{e} + \ln x )}{\ln x} \right\} , \qquad \text{for all } x > 1 \,.\] The exact constant \((2 \ln 2)^{-1/2}\) is familiar from asymptotics of binomial moderate deviations [4].
The lower bound is obtained by an embedding into the Hamming cube along monotone paths. Systems of integer equations are tuned using binomial coefficient estimates. For each of these we exhibit \(f_Z \in \mathcal{M}\) whose deterministic separation from its expectation gives the lower bound. The upper bound proceeds from the elementary propagation property of the deterministic monotone functions in \(\mathsf{M}\) (Figure 1). Parameter-side uniformisation isolates this property by means of an envelope of gauge functionals; optimising the envelope at each \(x\) gives the stated upper bound. Coincidence of the upper and lower bounds gives the constant \(\sqrt{2\ln2}\) and the iterated logarithmic term.
The class \(\mathcal{M}\) exhibits complex dependence across \(t\). Standard entropy and chaining methods from the concentration-of-measure literature [5]–[9] do not, on their own, identify the sharp leading constant for the full class \(\mathcal{M}\). Directly comparable are the standard bracketing bounds for the subclass of averages of monotone summands [10] in \(\mathcal{M}\). Those bounds are exponentially loose relative to the bound of Theorem 1 for the subclass. The reason is that those bounds normalise by \(\sqrt{n}\), the number of independent summands, rather than by \(\sqrt{\eta_f}\). For the empirical distribution the two normalisations coincide, because \(\eta_{F_n} = n\). However, the effective sample size may far outstrip the number of independent summands: under the strict stratification of Section 6.1, where each of the \(n\) independent blocks contributes \(m_n\) dependent ordered observations, it is \(\eta_f = n m_n^2\), and the bracketing normalisation is loose by a factor that grows with the block size (Remark 30). This same construction attains the sharp bound of Theorem 1: the configuration of average summands on which bracketing is loose is the one on which our logarithmic bound is sharp.
The smooth Gaussian process of Diebolt and Posse [11] is another comparison. Their suprema satisfy an \(\ell_2\) arc-length tail bound, derived from the differential geometry of the level sets. The bounded distribution-free analogue, with smooth non-monotone paths, is Corollary 5, obtained from our main theorem. The two bounds have the same form, the arc length now measured in \(\ell_1\) rather than \(\ell_2\): their bound rests on the curvature of a curve on the \(\ell_2\)-sphere, ours on the variation of a curve on the \(\ell_1\)-sphere. They differ only in the constant, which for us is the sharp logarithmic penalty of the main theorem (Remark 31). Monotonicity buys Gaussian-like geometry across all distributions, at the cost of boundedness and the logarithmic penalty.
The remainder of the paper is organised as follows. Section 2 establishes the basic properties of \(U_{\mathcal{M}}\), including the equivalence framework and the measurability of the centred supremum. Section 3 proves the lower bound using Hamming-cube embeddings, tuned admissible constellations, and an interpolation estimate for the resulting discrete sequence. Section 4 develops a parameter-side uniformisation, the classical data-side approach being unavailable here. Section 5 uses it to prove the upper bound, through a propagation property of the monotone functions in \(\mathsf{M}\). Section 6 gives applications. We first present an empirical distribution function for stratified clustered data that asymptotically saturates the constant of the main theorem. We then extend the bounds to non-monotone processes, bounding the suprema of smooth bounded processes by their \(\ell_1\) arc length. The Appendix collects auxiliary results and deferred proofs.
This section establishes the following preliminary result.
Proposition 1. The function \(U_{\mathcal{M}}\) is well defined and monotone on \(\mathbb{R}_+\); \(U_{\mathcal{M}}(0)=0\), \(U_{\mathcal{M}}(x)>0\) for all \(x>0\), and \[U_{\mathcal{M}}(x) \ge \frac{1}{2 \ln 2}, \qquad \text{for every } x \ge 1\,.\] Moreover, \[\mathbb{P}\left\{ \sqrt{\eta_f}\, \sup_{t \in \mathbb{R}} \left| f_Z(t) - \mathbb{E}\left[f_Z(t)\right] \right| \ge \varepsilon \sqrt{U_{\mathcal{M}}(\eta_f)} \right\} \le 2 \mathrm{e}^{-2\varepsilon^2} \,, \qquad \text{for all } \varepsilon > 0 \,,\] for every \(f_Z \in \mathcal{M}\) with \(\eta_f < \infty\).
The proof draws on several lemmas that will be used throughout the paper.
The following notion of equivalence is purely parameter-side.
Definition 1. Two randomised functions \(f_Z \colon \Theta \to \mathbb{R}\) and \(g_Z \colon T \to \mathbb{R}\) over a common product measurable space \(\mathcal{Z}\) and random vector \(Z\) are equivalent if \[\sqrt{\eta_f}\, \sup_{\theta \in \Theta} \left| f(z, \theta) - \mathbb{E}\left[f_Z(\theta)\right] \right| = \sqrt{\eta_g}\, \sup_{t \in T} \left| g(z, t) - \mathbb{E}\left[g_{Z}(t)\right] \right|, \qquad \text{for each } z \in \mathcal{Z}\,.\]
The defining equation is not defined when \(\eta_f = \infty\) or \(\eta_g = \infty\); we adopt the convention that the functions with \(\eta = \infty\) form a single equivalence class. Such functions are deterministic for every choice of \(Z\), so their centred supremum \[\sup_{\theta \in \Theta} \left| f(z, \theta) - \mathbb{E}\left[f_Z(\theta)\right] \right|\] is zero for each \(z \in \mathcal{Z}\).
We give a construction that preserves equivalence and will be used often in the sequel; its proof is immediate.
Lemma 2. Let \(T\) be a set, let \(\psi \colon T \to \Theta\) be onto, let \(h \colon T \to \mathbb{R}\) be a bounded deterministic function, and let \(a \neq 0\). If \(f_Z \colon \Theta \to \mathbb{R}\) is a randomised function, then \[g(z, t) := a\, f(z, \psi(t)) + h(t)\] defines a randomised function \(g_Z \colon T \to \mathbb{R}\) equivalent to \(f_Z\); moreover, \(a^2 \eta_g = \eta_f\).
We next define the optimal admissible constant \(u_{f_Z}\) attached to a single randomised function.
Definition 2. For an arbitrary randomised function \(f_Z \colon \Theta \to \mathbb{R}\) with finite \(\eta_f\), denote by \(u_{f_Z}\) the infimum of all \(u \ge 0\) satisfying \[\mathbb{P}\left\{ \sqrt{\eta_f} \sup_{\theta \in \Theta} \left| f_Z(\theta) - \mathbb{E}\left[f_Z(\theta)\right] \right| \ge \varepsilon \sqrt{u} \right\} \le 2\mathrm{e}^{-2\varepsilon^2}, \qquad \text{for all } \varepsilon > 0\,.\] Any \(u\) satisfying the above condition is called admissible for \(f_Z\). A function \(U \colon \mathbb{R}_+ \to \mathbb{R}_+\) is pointwise admissible if \(U(\eta_f)\) is admissible for every \(f_Z \in \mathcal{M}\) with \(\eta_f < \infty\).
We record the optimality property, which is immediate.
Lemma 3. Let \(X\) be a bounded non-negative random variable, and let \(A\) be the set of \(u \ge 0\) satisfying \[\mathbb{P}\{ X \ge \varepsilon \sqrt{u} \} \le 2 \mathrm{e}^{-2\varepsilon^2}, \qquad \text{for all } \varepsilon > 0 \,.\] If \(\mathbb{P}\{X > 0\} > 0\), then \(A = [u^*, \infty)\), where \[u^* = \sup_{\left\{t > 0 \colon \mathbb{P}\{X \ge t\} > 0\right\}} \frac{2t^2}{\ln\left(\frac{2}{\mathbb{P}\{X \ge t\}}\right)} > 0 \,,\] and otherwise \(A = (0,\infty)\).
The next result is a consequence of the previous characterisation of the admissible set.
Lemma 4. If \(f_Z \colon \Theta \to [0,1]\) is a randomised function with finite \(\eta_f\) and measurable centred supremum \[\xi = \sup_{\theta \in \Theta} |f_Z(\theta) - \mathbb{E}\left[f_Z(\theta)\right]|\,,\] then \(u_{f_Z} \le 2\eta_f / \ln 2\). Further, \(\mathbb{P}\{\xi > 0\} > 0\) if and only if \(u_{f_Z} > 0\), in which case \(u_{f_Z}\) is admissible.
The class \(\mathcal{M}\) contains scaled indicators of coordinate maxima, and these realise every effective sample size.
Lemma 5. For every \(x > 0\) there exists \(f_Z \in \mathcal{M}\) with \(\eta_f = x\) and \(u_{f_Z}>0\).
Proof. Choose \(n \in \mathbb{N}\) with \(nx \ge 1\), let \(\mathcal{Z}_i = \{0, 1\}\), and let \(Z_1, \ldots, Z_n\) be independent Bernoulli\((1/2)\) \(\{0,1\}\)-valued random variables. Define \[f_Z(t) := \frac{1}{\sqrt{nx}} \prod_{i=1}^n \mathbf{1}_{\{Z_i \le t\}}\,.\] Because \(nx \ge 1\), the function \(f_Z\) takes values in \([0,1]\) and is non-decreasing, so \(f_Z \in \mathcal{M}\). Flipping a single coordinate \(Z_i\) alters the product of indicators by at most \(1\), with equality attained, so \(c_i(f) = 1/\sqrt{nx}\). Summing the squares gives \(\eta_f = \bigl( n \cdot (1/\sqrt{nx})^2 \bigr)^{-1} = x\).
For every realisation \(z\), the function \(f_z\) is not equal to its expectation function, so the centred supremum is positive. Thus, \(u_{f_Z}>0\) by Lemma 4. ◻
We show the measurability of the centred supremum of randomised monotone functions.
Lemma 6. For every \(f_Z \in \mathcal{M}\), the map \[z \mapsto \sup_{t \in \mathbb{R}} \left| f(z, t) - \mathbb{E}\left[f_Z(t)\right] \right|\] is measurable. In particular, \(u_{f_Z}\) is well defined for every \(f_Z \in \mathcal{M}\) with finite \(\eta_f\).
Proof. We use the following elementary fact, whose proof is a standard exercise in monotonicity: if \(g, h \in \mathsf{M}\) and \(Q_g := \mathbb{Q} \cup D_g\), where \(D_g\) is the countable set of discontinuities of \(g\), then \[\sup_{t \in \mathbb{R}} \left| h(t) - g(t) \right| = \sup_{t \in Q_g} \left| h(t) - g(t) \right|\,.\] Let \(f_Z \in \mathcal{M}\), and write \(\mathbf{F}_Z(t) := \mathbb{E}\left[f_Z(t)\right]\) for the expectation of the randomised monotone function. Because \(f_z \in \mathsf{M}\) for each \(z\), the function \(\mathbf{F}_Z\) also belongs to \(\mathsf{M}\), as a pointwise expectation of members of \(\mathsf{M}\). Hence, for each \(z \in \mathcal{Z}\), \[\sup_{t \in \mathbb{R}} \left| f_z(t) - \mathbf{F}_Z(t) \right| = \sup_{t \in Q_{\mathbf{F}_Z}} \left| f_z(t) - \mathbf{F}_Z(t) \right|\,.\] The set \(Q_{\mathbf{F}_Z}\) is countable and does not depend on \(z\). For each \(t \in Q_{\mathbf{F}_Z}\), the map \[z \mapsto |f(z, t) - \mathbf{F}_Z(t)|\] is measurable, since \(z \mapsto f(z, t)\) is measurable and \(\mathbf{F}_Z(t)\) does not depend on \(z\). Therefore \[z \mapsto \sup_{t \in Q_{\mathbf{F}_Z}} \left| f(z, t) - \mathbf{F}_Z(t) \right|\] is measurable, being the supremum of a countable family of measurable functions, and the first displayed equality completes the proof. ◻
We restate previous results for the monotone class, in light of the established measurability.
Lemma 7. For any \(f_Z \in \mathcal{M}\) with finite \(\eta_f\), every \(u \ge u_{f_Z}\) with \(u > 0\) is admissible for \(f_Z\).
Proof. By Lemma 6 the centred supremum is measurable, so Lemma 4 applies. If \(u_{f_Z} = 0\), the centred supremum is almost surely zero, and every \(u > 0\) is admissible. If \(u_{f_Z} > 0\), then \(u_{f_Z}\) is admissible by Lemma 4 and so is \(u\). ◻
Proof of Proposition 1. Define \(U \colon \mathbb{R}_+ \to \mathbb{R}_+\) by \(U(0) = 0\) and, for \(x > 0\), \[U(x) := \sup\{u_{f_Z} : f_Z \in \mathcal{M},\;\eta_f = x\}\,.\] By Lemma 5 the set is non-empty, and by Lemma 4 it is bounded above by \(2x/\ln 2\), so \(U\) is defined.
Admissibility of \(U\). The value at zero is irrelevant for admissibility, because no \(f_Z\in\mathcal{M}\) has \(\eta_f=0\). For \(x>0\), \(U(x) \ge u_{f_Z}\) for every \(f_Z \in \mathcal{M}\) with \(\eta_f=x\). Because \(x>0\), there exists by Lemma 5 \(f_Z\in \mathcal{M}\) with \(\eta_f=x\) and \(u_{f_Z}>0\). Thus, \(U(x)>0\) and is admissible by Lemma 7.
Identification \(U_\mathcal{M} = U\). For any pointwise admissible \(V \colon \mathbb{R}_+ \to \mathbb{R}_+\) and any \(f_Z \in \mathcal{M}\) with \(\eta_f = x\), the bound at \(V(x)\) holds, so \(V(x) \ge u_{f_Z}\). Taking the supremum over \(f_Z\) gives \(V(x) \ge U(x)\), and taking the infimum over pointwise admissible \(V\) gives \(U_\mathcal{M}(x) \ge U(x)\). Conversely, \(U\) is itself pointwise admissible by the previous paragraph, so \(U_\mathcal{M}(x) \le U(x)\). Hence \(U_\mathcal{M} = U\), and the displayed inequality of the proposition is the admissibility of \(U_\mathcal{M}\) at \(\eta_f\).
Monotonicity. Let \(y \ge x > 0\) and let \(f_Z \in \mathcal{M}\) with \(\eta_f = x\). Define \(g_Z := \sqrt{x/y}\, f_Z\). Because \(0<\sqrt{x/y}\le 1\), the function \(g_Z\) still takes values in \([0,1]\) and is non-decreasing. Thus, \(g_Z \in \mathcal{M}\). By Lemma 2 with \(a = \sqrt{x/y}\), \(\eta_g = y\) and \(g_Z\) is equivalent to \(f_Z\), hence \(u_{g_Z} = u_{f_Z}\). Taking the supremum over \(f_Z\) with \(\eta_f = x\) gives \(U_\mathcal{M}(y) \ge U_\mathcal{M}(x)\).
Lower bound. Specialise the construction of Lemma 5 to \(n = 1\), so \(f_Z(t) = \mathbf{1}_{\{Z_1 \le t\}}\) with \(\eta_f = 1\), and let \(Z_1\) be uniform on \(\{0, 1\}\). For each \(t \in [0, 1)\), \(f_Z(t) = \mathbf{1}_{\{Z_1 = 0\}}\) takes the values \(0\) and \(1\) with equal probability, so \(|f_Z(t) - \mathbb{E}\left[f_Z(t)\right]| = \tfrac{1}{2}\). Hence \(\sqrt{\eta_f}\, \sup_t |f_Z(t) - \mathbb{E}\left[f_Z(t)\right]| = \tfrac{1}{2}\) surely.
Because \(\sqrt{\eta_f}\,\sup_{t \in \mathbb{R}}|f_Z(t) - \mathbb{E}\left[f_Z(t)\right]| = 1/2\) surely, Lemma 3 gives \(u_{f_Z} = \frac{2(1/2)^2}{\ln 2} = 1/(2\ln 2)\). Hence \[U_\mathcal{M}(1) \ge 1/(2 \ln 2)\,,\] and monotonicity extends this to every \(x \ge 1\). ◻
Remark 8 ().
The class \(\mathcal{M}\) includes randomised monotone functions \(f_Z\) that need not be jointly measurable in \((z, t)\). Were the class restricted to
jointly measurable members, all results would hold, and no improvement would be obtained in Theorem 1: the lower bound constructions used in its proof are jointly measurable, the
randomisations having finite support; and neither Proposition 1 nor the upper bound argument requires joint measurability. Monotonicity does the work usually done by joint measurability.
\(\diamond\)
Remark 9 ().
The lower bound \(1/(2\ln 2)\) for \(U_{\mathcal{M}}(1)\) can be improved by applying the optimal adjustment for the one-sample empirical distribution function.
Take \(Z_1\) uniform on \([0,1]\) and let \[F_1(t) := \mathbf{1}_{\{Z_1 \le t\}}\,, \qquad t \in \mathbb{R}\,,\] the one-sample empirical distribution function. Thus, \(F_1 \in \mathcal{M}\) and \(\eta_{F_1} = 1\). Because \(\mathbb{E}\left[F_1(t)\right] = \mathbb{P}\{Z_1 \le t\} = t\) for \(t \in [0,1]\), and the difference vanishes outside \([0,1]\), for every realisation \(Z_1=z\), \[\sup_{t \in \mathbb{R}} \left| F_1(t) - \mathbb{E}\left[F_1(t)\right] \right| = \max\{z,\, 1 - z\} \,.\] Thus, the centred supremum is \[\xi := \max\{Z_1,\, 1 - Z_1\},\] the two-sided one-sample Kolmogorov–Smirnov statistic. For \(s \in [\tfrac{1}{2}, 1]\), \[\mathbb{P}\{\xi \le s\} = \mathbb{P}\{1 - s \le Z_1 \le s\} = 2s - 1\,,\] and therefore \(\mathbb{P}\{\xi \ge s\} = 2(1 - s)\) on the same interval, while \(\mathbb{P}\{\xi \ge s\} = 1\) for \(s \in (0, \tfrac{1}{2}]\).
Because \(\eta_{F_1} = 1\), Lemma 3 applied to \(X = \xi\) gives \[u_{F_1} = \sup_{\{s > 0\,\colon\, \mathbb{P}\{\xi \ge s\} > 0\}} \frac{2 s^2}{\ln\!\big(2/\mathbb{P}\{\xi \ge s\}\big)}\,.\] For \(s \in (0, \tfrac{1}{2}]\), \(\mathbb{P}\{\xi \ge s\} = 1\) and the ratio is \(2s^2/\ln 2\), taking value \(1/(2 \ln 2)\) at \(s = \tfrac{1}{2}\). For \(s \in (\tfrac{1}{2}, 1)\), \[\frac{2 s^2}{\ln\!\big(2/2(1-s)\big)} = \frac{2 s^2}{\ln\!\big(1/(1-s)\big)}\,.\] Writing \[h(s) = \frac{2s^2}{\ln(1/(1-s))}, \qquad s\in(\tfrac{1}{2},1)\,,\] the equation \(h'(s)=0\) reduces to \[2(1-s)\ln\!\big(1/(1-s)\big) = s\,.\] Equivalently, with \(r=1-s\in(0,\tfrac{1}{2})\), \[2r\ln(1/r)+r-1=0\,.\] The left-hand side is strictly increasing in \(r\) on \((0,\tfrac{1}{2})\), has limit \(-1\) as \(r\downarrow0\), and is positive at \(r=\tfrac{1}{2}\); hence there is a unique root, equivalently a unique \(s^* \in (\tfrac{1}{2},1)\). At this root, \[h(s^*) = \frac{2(s^*)^2}{\ln(1/(1-s^*))} = \frac{2(s^*)^2 \cdot 2(1-s^*)}{s^*} = 4\,s^*(1 - s^*)\,.\] Numerically \(s^* \approx 0.7153\), so \[h(s^*) \approx 0.8145 > \frac{1}{2\ln 2}\,.\] Because \(h(s)\to0\) as \(s\to1^-\) and \(h(\tfrac{1}{2}) = 1/(2\ln2)\), the supremum is attained at \(s^*\). Hence \[u_{F_1} = 4\,s^*(1 - s^*) \approx 0.8145 \,.\] This is the optimal adjustment for the one-sample empirical-distribution case. Therefore \[U_{\mathcal{M}}(1) \ge u_{F_1} = 4\,s^*(1-s^*) \approx 0.8145 > \frac{1}{2\ln 2}\,.\] This refinement does not affect the sharp bound of Theorem 1; its interest is pedagogical. \(\diamond\)
We prove the following theorem.
Theorem 2. For all \(x > 1\), \[\left( \sqrt{U_{\mathcal{M}}(x)} - \sqrt{\log_4 x} \right)^- \le 2 \min\left\{1, \frac{2\ln(\mathrm{e}+ \ln x)}{\sqrt{\ln x}}\right\} \,.\]
The proof uses Lemma 10. The lemma applies to every nonempty subset of the Hamming cube and converts its covering radius and shortest-path traversal length into a lower bound on \(\sqrt{U_{\mathcal{M}}}\).
For \(n \in \mathbb{N}\), let \(\mathcal{Q}_n := \{0,1\}^n\) denote the discrete cube equipped with the Hamming distance \(d_H(q,q') = \#\{i \colon q_i \neq q'_i\}\). For any nonempty \(S \subseteq \mathcal{Q}_n\), define the covering radius \[\rho(S) = \max_{q \in \mathcal{Q}_n} \min_{s \in S} d_H(q, s) \,,\] and the shortest path traversal length \[\mathcal{L}(S) = \min \left\{ \sum_{j=2}^{|S|} d_H(s_j, s_{j-1}) \colon (s_1, s_2, \dots, s_{|S|}) \text{ is an ordering of } S \right\} \,.\]
The following result ties Hamming cube geometry to randomised monotone functions.
Lemma 10. For every nonempty \(S \subseteq \mathcal{Q}_n\), \[\sqrt{U_{\mathcal{M}}\left( \frac{(\mathcal{L}(S) + n)^2}{n} \right)} \ge \frac{n - 2\rho(S)}{\sqrt{2n \ln 2}} \,.\]
We record a consequence.
Corollary 1. For all \(x \ge 4\), \[\sqrt{U_{\mathcal{M}}(x)} \ge \frac{1}{\sqrt{2\ln 2}} \sqrt{\log_4 x - 1} \,.\]
Proof. Let \(x \ge 4\) and define \(n_x := \lfloor \log_4 x \rfloor\). Set \(S = \mathcal{Q}_{n_x}\). We have \(\rho(S) = 0\). The hypercube \(\mathcal{Q}_{n_x}\) admits a Hamiltonian path with \(\mathcal{L}(S) = 2^{n_x} - 1\).
We have \[(2^{n_x} - 1 + n_x)^2 \le n_x 4^{n_x}\,.\] Indeed, \[2^{n_x}+n_x-1 \le \sqrt{n_x}\,2^{n_x}\,.\] For \(n_x=1\) this is equality. For \(n_x\ge2\), it follows because \[n_x-1=(\sqrt{n_x}-1)(\sqrt{n_x}+1) \le (\sqrt{n_x}-1)2^{n_x}\,.\] Thus, because \(\mathcal{L}(S)=2^{n_x}-1\) and \(4^{n_x}\le x\), \[\frac{(\mathcal{L}(S)+n_x)^2}{n_x} \le 4^{n_x} \le x \,.\]
Lemma 10 and the monotonicity of \(U_{\mathcal{M}}\) imply \[\sqrt{U_{\mathcal{M}}(x)} \ge \frac{n_x}{\sqrt{2n_x\ln 2}} = \frac{1}{\sqrt{2\ln 2}}\sqrt{n_x} \,.\] Because \(n_x \ge \log_4 x - 1\), the claim follows. ◻
Proof of Lemma 10. If \(n - 2\rho(S) \le 0\), the lower bound is non-positive and holds because \(U_{\mathcal{M}}\) is non-negative. Assume \(n - 2\rho(S) > 0\).
Identify \(\{0,1\}^n\) with the oriented cube \(\mathcal{Z} = \{-1,1\}^n\) via the standard bijection. Identify \(S\) with a subset \(\Theta \subseteq \mathcal{Z}\), preserving \(\rho(\Theta) = \rho(S)\) and \(\mathcal{L}(\Theta) = \mathcal{L}(S)\). Define \[\label{eq:f-extreme} f(z,\theta) = \langle z, \theta \rangle \,.\tag{3}\]
Because \(\langle z, \theta \rangle = n - 2d_H(z,\theta)\), we have \[\min_{z\in \mathcal{Z}} \max_{\theta\in \Theta} \langle z, \theta \rangle = n - 2\rho(\Theta) \,.\] This establishes \(\sup_{\theta\in \Theta} f(z,\theta) \ge n - 2\rho(\Theta)\) for all \(z\in \mathcal{Z}\).
We have \(c_i(f) = 2\) for all \(i\), because changing \(z_i\) to \(-z_i\) alters \(\langle z, \theta\rangle\) by \(\pm 2\theta_i\) with \(\theta_i \in \{-1, 1\}\). This gives \(\eta_f = 1/(4n)\). Thus, \[\sqrt{\eta_f} \sup_{\theta\in \Theta} f(z,\theta) \ge \frac{n - 2\rho(\Theta)}{2\sqrt{n}}, \qquad \text{for every } z \in \mathcal{Z} \,.\]
Let \(Z\) be uniformly distributed over \(\mathcal{Z}\). Because \(\mathbb{E}\left[f_Z(\theta)\right]=0\) and \(n - 2\rho(\Theta) > 0\), we have \[\sqrt{\eta_f} \sup_{\theta\in\Theta} |f(z,\theta) - \mathbb{E}\left[f_Z(\theta)\right]| \ge \sqrt{\eta_f} \max_{\theta\in\Theta} f(z,\theta) \ge \frac{n - 2\rho(\Theta)}{2\sqrt{n}}\] for every \(z\in\mathcal{Z}\). In particular, \[\mathbb{P}\left\{ \sqrt{\eta_f} \sup_{\theta\in \Theta} \left| f_Z(\theta) - \mathbb{E}\left[f_Z(\theta)\right] \right| \ge \frac{n - 2\rho(\Theta)}{2\sqrt{n}} \right\} = 1 \,.\] Therefore, \[\sqrt{u_{f_Z}} \ge \frac{n - 2\rho(\Theta)}{\sqrt{2n \ln 2}} \,.\] To see this note that for arbitrary \(0<u < (n - 2\rho(\Theta))^2/(2n\ln 2)\), taking \(\varepsilon = (n - 2\rho(\Theta))/(2\sqrt{n}\sqrt{u})\) gives \(\varepsilon\sqrt{u} = (n - 2\rho(\Theta))/(2\sqrt{n})\) and \(2\mathrm{e}^{-2\varepsilon^2} < 1\). Thus, \(u_{f_Z}\ge u\).
Let \(m = |\Theta|\) and order \(\Theta = \{\theta_1, \theta_2, \ldots, \theta_m\}\) to match the shortest path \(\mathcal{L}(\Theta)\). Define the surjective mapping \(\sigma \colon \mathbb{R} \to \Theta\) by \(\sigma(t) = \theta_{j(t)}\), where \[j(t) = \max\{1, \min\{m, \lfloor t \rfloor\}\} \,.\] Define \(g \colon \mathcal{Z} \times \mathbb{R} \to \mathbb{R}\) by \[g(z,t) = \frac{1}{2(\mathcal{L}(\Theta) + n)}\left( n + f(z, \sigma(t)) + 2\sum_{k=2}^{j(t)} d_H(\theta_k, \theta_{k-1}) \right) \,,\] where the sum evaluates to zero for \(j(t)=1\).
Because \(f(z, \theta) \in [-n, n]\) and the sum of successive distances is at most \(\mathcal{L}(\Theta)\), the mapping \(g_Z\) takes values in \([0,1]\). If \(j\ge2\), then \[f(z,\theta_j)-f(z,\theta_{j-1}) \ge -2d_H(\theta_j,\theta_{j-1})\,.\] The added term therefore compensates for any fall in the inner product, so \(g_Z(t)\) is non-decreasing in \(t\). Thus, \(g_Z \in \mathcal{M}\).
By Lemma 2, \(g_Z\) is equivalent to \(f_Z\) with \[\eta_{g} = 4(\mathcal{L}(\Theta) + n)^2 \eta_f = \frac{(\mathcal{L}(\Theta)+n)^2}{n} \,.\] Evaluating \(U_{\mathcal{M}}\) gives \[\sqrt{U_{\mathcal{M}}\left(\frac{(\mathcal{L}(\Theta)+n)^2}{n}\right)} = \sqrt{U_{\mathcal{M}}(\eta_g)} \ge \sqrt{u_{g_Z}} = \sqrt{u_{f_Z}} \ge \frac{n - 2\rho(\Theta)}{\sqrt{2n \ln 2}} \,,\] completing the proof. ◻
We prove a consequence of Proposition 1 and Corollary 1.
Lemma 11. For all \(1 < x \le 4^{111}\), the bound in Theorem 2 is satisfied.
Proof. For \(1 < x \le 4^8\), the bound holds by the lower estimate in Proposition 1. Specifically, Proposition 1 establishes \(U_{\mathcal{M}}(x) \ge (2\ln 2)^{-1}\) for all \(x\ge 1\). Note that \[\sqrt{\log_4(4^8)} - (2\ln 2)^{-1/2} = \sqrt{8} - (2\ln 2)^{-1/2} < 2\,,\] and \[2 \min\left\{1, \frac{2\ln(\mathrm{e}+ \ln x)}{\sqrt{\ln x}}\right\} = 2, \qquad \text{for } 1< x\le 4^{55}\,.\] Together these establish the bound for \(1 < x \le 4^8\).
By Corollary 1, for \(x \ge 4^8\), the residual gap is bounded by \[\sqrt{\log_4 x} - \sqrt{U_{\mathcal{M}}(x)} \le \sqrt{\log_4 x} - \frac{1}{\sqrt{2\ln 2}}\sqrt{\log_4 x - 1} \,.\] Letting \(y = \log_4 x\) and \(a = 2\ln 2\), direct calculation shows that \[\sqrt{y} - \frac{\sqrt{y-1}}{\sqrt{a}} \le 2, \qquad 1 \le y \le 170\,,\] and \[\sqrt{a}y -\sqrt{y^2-y} \le 4\ln(\mathrm{e}+ay), \qquad 1 \le y \le 111\,.\] These together give \[\sqrt{y} - \frac{\sqrt{y-1}}{\sqrt{a}} \le 2\min\left\{1, \frac{2\ln(\mathrm{e}+ay)}{\sqrt{ay}}\right\} \,,\] for the remaining range \(8 \le y \le 111\). ◻
Admissible constellations are the integer parameter systems used to extend the lower bound beyond the deterministic range.
Definition 3. A tuple \((k,t_k,n_k,r_k,m_k)\in\mathbb{N}^5\) is a constellation if \[\begin{gather} n_k = k\,t_k^2, \qquad r_k = \frac{k\,t_k(t_k-1)}{2}, \\[2mm] t_k \ge 2\,. \end{gather}\] It is an admissible constellation if \(m_k\le 2^{n_k}\).
Here, \(k\) is the sequence index, \(t_k\) is the tuning parameter, \(n_k\) is the dimension of the Hamming cube, \(r_k\) is the target covering radius, and \(m_k\) is the cardinality of the target subset within the \(n_k\)-dimensional cube. The associated volume is \(n_k m_k^2\).
The next lemma records that admissible constellations are generated by the three integer parameters \(k\), \(t_k\), and \(m_k\).
Lemma 12. Let \(k,t_k\in\mathbb{N}\) with \(t_k\ge2\). Define \[n_k := k\,t_k^2, \qquad r_k := \frac{k\,t_k(t_k-1)}{2}\,.\] Hence \(n_k,r_k\in\mathbb{N}\). Consequently, for every \(m_k\in\mathbb{N}\) satisfying \(m_k\le 2^{n_k}\), the tuple \[(k,t_k,n_k,r_k,m_k)\] is an admissible constellation.
Proof. That \(n_k\) is an integer is immediate. The number \(r_k\) is an integer because \(t_k(t_k-1)\) is even. The remaining defining conditions are \(t_k\ge2\) and \(m_k\le2^{n_k}\). ◻
We now apply Lemma 10 to admissible constellations.
Lemma 13. For an admissible constellation \((k, t_k, n_k, r_k, m_k)\), if there exists a nonempty \(S \subseteq \mathcal{Q}_{n_k}\) satisfying \(|S| \le m_k\) and \(\rho(S) \le r_k\), then \[\sqrt{U_{\mathcal{M}}(n_k m_k^2 )} \ge \frac{1}{\sqrt{2\ln 2}} \sqrt{k} \,.\]
Proof. The maximal path length is bounded by \[\mathcal{L}(S) \le n_k(|S|-1) \,,\] which implies \((\mathcal{L}(S) + n_k)^2/n_k \le n_k|S|^2\). Thus, the monotonicity of \(U_{\mathcal{M}}\) and Lemma 10 imply \[\begin{align} \sqrt{U_{\mathcal{M}}(n_k m_k^2 )} & \ge \sqrt{U_{\mathcal{M}} (n_k |S|^2 )} \\ & \ge \sqrt{U_{\mathcal{M}} \left( \frac{(\mathcal{L}(S) + n_k)^2}{n_k} \right)} \\ & \ge \frac{1}{\sqrt{2\ln 2}} \frac{n_k - 2\rho(S)}{\sqrt{n_k}} \\ & \ge \frac{1}{\sqrt{2\ln 2}} \frac{n_k - 2r_k}{\sqrt{n_k}} \,. \end{align}\] Because any constellation satisfies \[\frac{n_k - 2r_k}{\sqrt{n_k}} = \sqrt{k} \,,\] the stated bound is established. ◻
A standard probabilistic argument establishes sufficient conditions for the existence of the requisite subset of Lemma 13.
Lemma 14. For an admissible constellation \((k, t_k, n_k, r_k, m_k)\), if \[\sum_{j=0}^{r_k} \binom{n_k}{j} \ge \frac{2^{n_k}}{m_k} n_k \ln 2 \,,\] then \[\sqrt{U_{\mathcal{M}} (n_k m_k^2 )} \ge \frac{1}{\sqrt{2\ln 2}} \sqrt{k} \,.\]
Proof. Let \(s_1, \dots, s_{m_k}\) be independent random variables uniformly distributed over \(\mathcal{Q}_{n_k}\). Applying the union bound over all \(q \in \mathcal{Q}_{n_k}\) we obtain \[\begin{align} \mathbb{P}\left\{ \max_{q\in \mathcal{Q}_{n_k}} \min_{1 \le j \le m_k} d_H(q,s_j) > r_k \right\} & \le 2^{n_k}\left(1 - \frac{1}{2^{n_k}} \sum_{j=0}^{r_k}\binom{n_k}{j}\right)^{m_k} \\ & < 2^{n_k} \exp\left[- \frac{m_k}{2^{n_k}} \sum_{j=0}^{r_k}\binom{n_k}{j}\right] \,. \end{align}\] Under the stated hypothesis, this probability is strictly less than \(1\). This establishes the existence of a set \(S \subseteq \mathcal{Q}_{n_k}\) satisfying \(|S| \le m_k\) and \(\rho(S) \le r_k\). The conclusion follows directly from Lemma 13. ◻
Retaining only one binomial term, at the cost of slack, gives a closed-form sufficient condition for the preceding lemma.
Lemma 15. For an admissible constellation \((k, t_k, n_k, r_k, m_k)\), if \[\ln(m_k^2) \ge k + 3\ln k + \left[6\ln t_k + \frac{k}{6(t_k^2-1)}\right] + \left[\ln 2 + 2\ln(\ln 2)\right] \,,\] then \[\sqrt{U_{\mathcal{M}}(n_k m_k^2)} \ge \frac{1}{\sqrt{2\ln 2}} \sqrt{k} \,.\]
Proof. By Lemma 14, the conclusion holds if \[\frac{1}{2^{n_k}} \sum_{j=0}^{r_k} \binom{n_k}{j} \ge \frac{1}{m_k} n_k \ln 2 \,.\] Lemma 32 in Appendix 7.1 provides a lower bound on the relevant binomial shell: \[\frac{1}{2^{n_k}}\binom{n_k}{r_k} \ge \frac{1}{\sqrt{2k}\,t_k} \exp\left(-\frac{k}{12(t_k^2-1)}\right) \mathrm{e}^{-k/2} \,.\] Bounding the summation from below by this single term, the required condition is satisfied if \[\frac{1}{\sqrt{2k}\,t_k} \exp\left(-\frac{k}{12(t_k^2-1)}\right) \mathrm{e}^{-k/2} \ge \frac{1}{m_k} n_k \ln 2 \,.\] Substituting \(n_k = k t_k^2\), taking logarithms, and isolating \(\ln m_k\) we obtain \[\begin{align} \ln m_k & \ge \ln(k t_k^2 \ln 2) + \frac{1}{2}\ln(2k) + \ln t_k + \frac{k}{12(t_k^2-1)} + \frac{k}{2} \\ & = \frac{k}{2} + \frac{3}{2}\ln k + 3\ln t_k + \frac{k}{12(t_k^2-1)} + \frac{1}{2}\ln 2 + \ln(\ln 2) \,. \end{align}\] Multiplying by \(2\) establishes the stated bound. ◻
The sufficient condition in Lemma 15 is equivalently a lower bound on \(m_k\). We define the real parameter \(\tilde{m}_k\) by taking equality in its logarithmic form: \[2\ln(\tilde{m}_k) = k + 3\ln k + \left[6\ln t_k + \frac{k}{6(t_k^2-1)}\right] + \left[\ln 2 + 2\ln(\ln 2)\right] \,.\] Applying the identity \(n_k = k t_k^2\) we obtain \[\begin{align} \ln(n_k \tilde{m}_k^2) & = k + 4\ln k + 8\ln t_k + \frac{k}{6(t_k^2-1)} + \left[\ln 2 + 2\ln(\ln 2)\right] \\ & < k + 4\ln k + \left[8\ln t_k + \frac{k}{6(t_k^2-1)}\right] \,, \end{align}\] where the strict inequality follows from \(\ln 2 + 2\ln(\ln 2) < 0\).
Minimising the leading dominant terms of the tuning penalty \[8\ln t_k + \frac{k}{6(t_k^2-1)} \,,\] gives \[t_k \sim \frac{\sqrt{k}}{\sqrt{24}} \,.\] We set the calculation-friendly real tuning parameter \[\tilde{t}_k = \frac{\sqrt{k}}{5} \,.\]
Definition 4. A constellation \((k, t_k, n_k, r_k, m_k)\) is tuned if \[t_k = \left\lceil \tilde{t}_k \right\rceil \quad\text{and}\quad m_k = \lceil \tilde{m}_k \rceil \,.\] If it is also admissible, then it is an admissible tuned constellation.
We now establish the existence of admissible tuned constellations.
Lemma 16. For all \(k \ge 26\), there exists a unique admissible tuned constellation \[(k, t_k, n_k, r_k, m_k) \,.\] Furthermore, this constellation satisfies \[\sqrt{U_{\mathcal{M}}(n_k m_k^2 )} \ge \frac{1}{\sqrt{2\ln 2}} \sqrt{k} \,.\]
Proof. Fix \(k \ge 26\). The formulas \[t_k=\left\lceil\tilde{t}_k\right\rceil, \qquad m_k=\left\lceil\tilde{m}_k\right\rceil\] uniquely determine \(n_k\) and \(r_k\). It remains to verify admissibility.
Because \(\tilde{t}_k=\sqrt{k}/5\) and \(k\ge26\), we have \[2\le t_k\le\sqrt{k}\,.\] Note also that \(\tilde{m}_k>1\). By its defining equation and \(t_k\ge2\) we have \[2\ln\tilde{m}_k \ge k+3\ln k+\ln2+2\ln(\ln2)>0\] for \(k\ge26\).
We have \[\begin{align} \ln m_k & = \ln \lceil \tilde{m}_k \rceil \le \ln(\tilde{m}_k + 1) \le \ln(2\tilde{m}_k) = \ln \tilde{m}_k + \ln 2 \\ & < \frac{k}{2} + \frac{3}{2}\ln k + 3\ln t_k + \frac{k}{12(t_k^2-1)} + \ln 2 \\ & \le \frac{k}{2} + \frac{3}{2}\ln k + 3\ln t_k + \frac{k}{36} + \ln 2 \\ & \le \frac{k}{2} + \frac{3}{2}\ln k + 3\ln \sqrt{k} + \frac{k}{36} + \ln 2 \\ & = \frac{19}{36}k + 3\ln k + \ln 2 \\ & \le \frac{19}{36}k + k \\ & = \frac{55}{36}k < 4k\ln 2 \le k t_k^2 \ln 2 = n_k\ln 2\,. \end{align}\] This establishes \(m_k \le 2^{n_k}\), and hence the tuned constellation is admissible.
Also, \(m_k=\lceil\tilde{m}_k\rceil\ge\tilde{m}_k\), so \(\ln(m_k^2)\ge\ln(\tilde{m}_k^2)\). By the definition of \(\tilde{m}_k\), the hypothesis of Lemma 15 is satisfied, and the final statement follows from that lemma. ◻
For each \(k \ge 26\), let \[(k, t_k, n_k, r_k, m_k)\] be the unique admissible tuned constellation, the existence of which is guaranteed by Lemma 16.
The next result locates each \(x>4^{111}\) between two consecutive tuned volumes.
Lemma 17. For all \(x > 4^{111}\), there exists an integer \(k_x \ge 120\) satisfying \[n_{k_x} m_{k_x}^2 \le x < n_{k_x+1} m_{k_x+1}^2 \,.\]
Proof. Fix \(x > 4^{111}\) and define the index set \[\mathcal{I}_x = \{k \ge 26 \colon n_k m_k^2 \le x \} \,.\] Because the sequence \((n_k m_k^2)\) diverges to infinity, the set \(\mathcal{I}_x\) has finite cardinality. We now show that it is not empty.
At the index \(k = 120\), we have \(t_{120} = 3\), \(n_{120} = 1080\), and expanding the logarithm gives \[\begin{align} \ln(n_{120} m_{120}^2) & \le \ln(n_{120}(\tilde{m}_{120}+1)^2) \\ & = \ln(n_{120}\tilde{m}_{120}^2) + 2\ln\left(1 + \frac{1}{\tilde{m}_{120}}\right) \\ & < 150.4 + 0.6 < 151 < 222\ln 2 < \ln x \,. \end{align}\] Thus, \(120 \in \mathcal{I}_x\), so the set is non-empty.
Define \(k_x = \max \mathcal{I}_x\). The maximality of \(k_x \ge 120\) ensures \(k_x+1 \notin \mathcal{I}_x\), establishing the bounds \[n_{k_x} m_{k_x}^2 \le x < n_{k_x+1} m_{k_x+1}^2 \,,\] and completing the proof. ◻
Combining Lemmas 15, 16, and 17, for each \(x > 4^{111}\) we have \[\begin{align} \sqrt{\log_4 x} - \sqrt{U_{\mathcal{M}}(x)} & \le \frac{1}{\sqrt{2\ln 2}} \left( \sqrt{\ln x} - \sqrt{k_x} \right) \\ & = \frac{1}{\sqrt{2\ln 2}} \left( \frac{\ln x - k_x}{\sqrt{\ln x} + \sqrt{\ln x - (\ln x - k_x)}} \right) \,, \end{align}\] where \(k_x \ge 120\) is the index specified in Lemma 17. The rearrangement is valid because \(k_x\le \ln x\). Indeed, \(m_{k_x}\ge\tilde{m}_{k_x}\), and the definition of \(\tilde{m}_{k_x}\) gives \[\ln(n_{k_x}m_{k_x}^2) \ge \ln(n_{k_x}\tilde{m}_{k_x}^2) \ge k_x\,.\] Because \(n_{k_x}m_{k_x}^2\le x\), it follows that \(k_x\le\ln x\).
We now estimate the gap \(\ln x - k_x\) using an interpolation based on Lemma 17.
Lemma 18. For \(x > 4^{111}\), we have \[\ln x - k_x < 8\ln(\ln x) \,.\]
Proof. Let \(y = \ln x\). By Lemma 17, we have \[\begin{align} y & \ge \ln(n_{k_x} m_{k_x}^2) \ge \ln(n_{k_x} \tilde{m}_{k_x}^2) \\ & = k_x + 4 \ln k_x + 8 \ln t_{k_x} + \frac{k_x}{6(t_{k_x}^2-1)} + \ln 2 + 2\ln(\ln 2) \\ & > k_x + 4 \ln k_x + 8 \ln t_{k_x} + \ln 2 + 2\ln(\ln 2) \,. \end{align}\] Substituting the lower bounds \(k_x \ge 120\) and \(t_{k_x} \ge 3\) establishes \(y > k_x + 27\). This verifies the strict inequality \[\label{eq:1kx} y > k_x + 1 \,.\tag{4}\]
Define \(K_x = k_x + 1\). By Lemma 17, we have \[\begin{align} y & \le \ln(n_{K_x} m_{K_x}^2) \le \ln(n_{K_x} \tilde{m}_{K_x}^2) + \frac{2}{\tilde{m}_{K_x}} \\ & < K_x + 4 \ln K_x + 8 \ln t_{K_x} + \frac{K_x}{6(t_{K_x}^2-1)} + \frac{2}{\tilde{m}_{K_x}} \,. \end{align}\] We bound the two tuning penalties. By definition, \[t_{K_x} \le \frac{\sqrt{K_x}}{5} + 1 = \frac{\sqrt{K_x}}{5}\left(1 + \frac{5}{\sqrt{K_x}}\right) \,.\] Combining this with the elementary inequality \(\ln(1+z) \le z\) for \(z > -1\) gives \[8 \ln t_{K_x} \le 4 \ln K_x - 8 \ln 5 + \frac{40}{\sqrt{K_x}} \,.\] Furthermore, the lower bound \(t_{K_x} \ge \frac{\sqrt{K_x}}{5}\) ensures \[\frac{K_x}{6(t_{K_x}^2-1)} \le \frac{25 K_x}{6(K_x-25)} = \frac{25}{6}\left(1 + \frac{25}{K_x-25}\right) \,.\]
For \(K_x\ge121\), the first four terms in the bracket are bounded below by their values at \(K_x=121\). Also \(\tilde{m}_{K_x}>1\), so \(2/\tilde{m}_{K_x}<2\). Direct calculation gives \[8\ln5-1-\frac{40}{11} -\frac{25}{6}\left(1+\frac{25}{96}\right) -2 > 0 \,.\] Substituting \(K_x = k_x + 1\) and evaluating the constants for \(K_x \ge 121\) gives \[\begin{align} y - k_x & \le 1 + 8 \ln K_x + \frac{40}{\sqrt{K_x}} - 8\ln 5 + \frac{25}{6}\left(1 + \frac{25}{K_x-25}\right) + \frac{2}{\tilde{m}_{K_x}} \\ & \le 8 \ln(k_x+1) - \left[ 8\ln 5 - 1 - \frac{40}{\sqrt{K_x}} - \frac{25}{6}\left(1 + \frac{25}{K_x-25}\right) - \frac{2}{\tilde{m}_{K_x}} \right] \\ & < 8\ln(k_x+1) \,. \end{align}\] Applying the logarithmic substitution \(\ln(k_x+1) < \ln y = \ln(\ln x)\) derived from 4 establishes the stated result. ◻
Proof of Theorem 2. For \(1 < x \le 4^{111}\) the bound holds by Lemma 11. We show the bound for \(x > 4^{111}\).
Fix \(x > 4^{111}\) and let \(y = \ln x\). By Lemmas 16 and 17, \[\sqrt{\log_4 x} - \sqrt{U_{\mathcal{M}}(x)} \le \frac{1}{\sqrt{2\ln 2}} \left( \frac{y - k_x}{\sqrt{y} + \sqrt{y - (y - k_x)}} \right) \,.\] The map \(u \mapsto \frac{u}{\sqrt{y} + \sqrt{y - u}}\) is strictly increasing for \(0 \le u < y\). By Lemma 18 and the condition \(y \ge 222\ln 2\), \[y - k_x < 8\ln y < y \,,\] so applying the monotonicity above, \[\begin{align} \frac{1}{\sqrt{2\ln 2}} \frac{y - k_x}{\sqrt{y} + \sqrt{y - (y - k_x)}} & < \frac{1}{\sqrt{2\ln 2}} \frac{8\ln y}{\sqrt{y} + \sqrt{y - 8\ln y}} \\ & = \frac{8\ln y}{\sqrt{y}} \left[ \sqrt{2\ln 2} \left( 1 + \sqrt{1 - \frac{8\ln y}{y}} \right) \right]^{-1} \,. \end{align}\] The function \(y \mapsto \frac{\ln y}{y}\) is strictly decreasing for \(y \ge \mathrm{e}\). For \(y > 222\ln 2\), the bracketed factor in the inverted expression exceeds \(2\). Therefore, \[\frac{1}{\sqrt{2\ln 2}} \frac{y - k_x}{\sqrt{y} + \sqrt{y - (y - k_x)}} < \frac{4\ln y}{\sqrt{y}} \,.\] Because \(\ln y < \ln(\mathrm{e} + y) = \ln(\mathrm{e} + \ln x)\), this establishes the bound stated in Theorem 2 for \(x > 4^{111}\). ◻
Remark 19 ().
The same extremal geometry in the proof of Lemma 10 is realised, up to the equivalence developed in Section 6.1, by
empirical distribution functions arising from \(n\) independent blocks of \(m_n\) strictly stratified observations. \(\diamond\)
Remark 20 ().
The deterministic range of Lemma 11 uses the full cube \(\mathcal{Q}_n\) for \(1 \le n
\le 111\); here \(n = 111\) is the largest dimension for which Corollary 1 verifies Theorem 2. The path-indexed formulation 3 extends to proper subsets \(S \subset \mathcal{Q}_n\) at
any \(n\). An admissible constellation asks for a subset of \(\mathcal{Q}_{n_k}\) of cardinality at most \(m_k\) and covering radius at most \(r_k\); this is a covering-code existence problem. Such subsets may be obtained from deterministic covering codes [12], [13] or exhaustive search; we do not pursue these. The non-constructive argument of Lemma 14 lets us tune the constellations. \(\diamond\)
Massart’s [1] upper bound argument proceeds by uniformising the data. That approach is not available for a general randomised monotone function. The construction in this section is reminiscent of Debreu’s gap lemma [14]: the expectation map is made strictly increasing, and its gaps are filled by interpolation. Equivalence is preserved throughout.
Let \(\mathcal{M}_{\mathsf{id}}([0,1])\) be the set of randomised monotone functions \(f_Z\colon [0,1]\to [0,1]\) satisfying the diagonal expectation condition: \[\mathbb{E}\left[f_Z(t)\right] = t \,,\quad \text{for all } t\in [0,1] \,.\] Let \(U_{\mathcal{M}_{\mathsf{id}}}(0)=0\), and for each finite \(x>0\) define \[U_{\mathcal{M}_{\mathsf{id}}}(x) := \sup\left\{ u_{f_Z} \colon f_Z \in \mathcal{M}_{\mathsf{id}}([0,1]), \;\eta_f = x \right\} \,.\]
The properties of \(U_{\mathcal{M}_{\mathsf{id}}}\) are established in the theorem of the present section.
Theorem 3. The function \(U_{\mathcal{M}_{\mathsf{id}}}\colon \mathbb{R}_+\to \mathbb{R}_+\) is well defined and monotone. Further, the right-continuous versions of \(U_{\mathcal{M}_{\mathsf{id}}}\) and \(U_{\mathcal{M}}\) are equal everywhere.
The proof requires an interpolation argument.
Let \(\Lambda\) be an arbitrary set and \(f_Z \colon \Lambda \times \Theta \to [0,1]\) be such that \(\theta \mapsto f_Z(\lambda, \theta)\) is a randomised function for each fixed \(\lambda \in \Lambda\). Assume that there is a maximal index \(\lambda^\ast \in \Lambda\) satisfying \[c_i(f(\lambda^\ast, \cdot)) \ge c_i(f(\lambda, \cdot))\,,\] for all \(i \in \{1, \dots, n\}\) and all \(\lambda \in \Lambda\), and \[\sup_{\theta\in \Theta} \left|f_Z(\lambda^\ast, \theta) - \mathbb{E}\left[f_Z(\lambda^\ast,\theta)\right]\right| \ge \sup_{\theta\in \Theta}\left|f_Z(\lambda, \theta) - \mathbb{E}\left[f_Z(\lambda, \theta)\right]\right| \,,\] for all \(\lambda \in \Lambda\).
Let \(T\) be a set satisfying \(\Theta \subseteq T\), and let \(h \colon T \to [0,1]^{\Lambda\times \Theta}\) be a partition of unity satisfying:
\(h_{(\lambda, \theta)}(t) = 0\), for all but finitely many \((\lambda,\theta)\), and \(\sum_{(\lambda,\theta)\in \Lambda\times \Theta} h_{(\lambda,\theta)} (t) = 1\) for all \(t \in T\).
For every \(\theta \in \Theta\), we have \(h_{(\lambda^\ast, \theta)}(\theta) = 1\).
Define the function \(g_Z \colon T \to [0,1]\) by \[g_Z(t) := \sum_{(\lambda,\theta) \in \Lambda\times \Theta} f_Z(\lambda, \theta) h_{(\lambda, \theta)}(t) \,.\]
Lemma 21. The function \(g_Z\) is a randomised function equivalent to \(f_Z(\lambda^\ast, \cdot)\), and \[c_i(g) = c_i(f(\lambda^\ast, \cdot)) \,,\] for all \(i \in \{1, \dots, n\}\).
Proof. Because \(h(t)\) is a partition of unity, \(g_Z(t)\) is a finite convex combination of the variables \(f_Z(\lambda, \theta) \in [0,1]\). For any coordinate \(i\), the triangle inequality ensures the variation of the sum is bounded by the convex combination of the maximal variations, giving \[c_i(g) \le c_i(f(\lambda^\ast, \cdot))\,.\] Similarly, centring the function and applying the triangle inequality gives \[\sup_{t \in T} \left|g_Z(t) - \mathbb{E}\left[g_Z(t)\right]\right| \le \sup_{\theta \in \Theta} \left|f_Z(\lambda^\ast, \theta) - \mathbb{E}\left[f_Z(\lambda^\ast, \theta)\right]\right|\,.\] Equality in both inequalities follows from the second property of the partition of unity. For any \(\theta \in \Theta\), evaluating at \(t = \theta\) gives \(h_{(\lambda^\ast, \theta)}(\theta) = 1\) and \(h_{(\lambda', \theta')}(\theta) = 0\) for all other pairs \((\lambda', \theta')\). Thus, \(g_Z(\theta) = f_Z(\lambda^\ast, \theta)\). Because the supremum over \(T\) is bounded below by the supremum over the subset \(\Theta \subseteq T\), exact equality holds. ◻
The proof of the theorem is an immediate consequence of the following constructive lemma.
Lemma 22. For any randomised function \(f_Z\) in \(\mathcal{M}\) and any \(a>0\), there is an equivalent randomised function \(g_Z\) in \(\mathcal{M}_{\mathsf{id}}([0,1])\) satisfying \[(1+ a)^{-2}\eta_g = \eta_f \,.\]
Proof. Restricting the domain. Fix \(f_Z \in \mathcal{M}\) and let \(a > 0\). Let \(\psi \colon (1/4,3/4) \to \mathbb{R}\) be a strictly increasing surjective function. Define \[f^{(1)}(z,t) := \begin{cases} 0, & t \le 1/4 \,, \\[4pt] f\left(z,\psi(t)\right), & t\in(1/4,3/4)\,, \\[4pt] 1, & t \ge 3/4 \,. \end{cases}\] By Lemma 2, and because outside \((1/4,3/4)\) the function is constant, the randomised function \(f^{(1)}_Z\colon [0,1]\to [0,1]\) is equivalent to \(f_Z\). Padding the domain ensures that subsequent steps do not require boundary conditions at the endpoints.
Strictly increasing expectations. Define the affine scaling, \[f^{(2)}(z,t) := \frac{1}{1+a} f^{(1)}(z,t) + \frac{a}{1+a} t \,,\qquad t\in[0,1] \,.\] By Lemma 2, the function \(f^{(2)}_Z\) is equivalent to \(f^{(1)}_Z\), and \[(1+a)^{-2} \eta_{f^{(2)}} = \eta_{f^{(1)}} \,.\] Hence, \(\eta_{f^{(2)}} = (1+a)^2 \eta_f\). Because \(a>0\), the added linear drift ensures the mapping \(t \mapsto f^{(2)}(z,t)\) is strictly increasing for every \(z\). Consequently, the expected function \(\mathbb{E}\left[f^{(2)}_Z(t)\right]\) is strictly increasing. Strict monotonicity ensures that the generalised inverse below is regular.
Generalised inverse. For every \(\vartheta\in [0,1]\), define the inverse mapping \[t_\vartheta = \inf \left\{ t \in [0,1] \colon \vartheta \le \mathbb{E}\left[f^{(2)}_Z(t)\right] \right\} = \sup \left\{ t \in [0,1] \colon \mathbb{E}\left[f^{(2)}_Z(t)\right] \le \vartheta \right\} \,.\] Let \(\Theta \subseteq [0,1]\) be the set of values \(\theta\) satisfying \(\theta = \mathbb{E}\left[f^{(2)}_Z(t_\theta)\right]\). Notice that \(\theta\mapsto t_\theta\) maps \(\Theta\) onto \([0,1]\), thus in particular \(f^{(3)}_Z(\theta) := f^{(2)}_Z(t_\theta)\) is equivalent to \(f^{(2)}_Z\) with the same coefficient. Further, \(\mathbb{E}\left[f^{(3)}_Z(\theta)\right] = \theta\) for all \(\theta\in \Theta\).
Interpolation. If the expected path is discontinuous, then \(\Theta\) has gaps. We use interpolation to fill these gaps and extend the domain from \(\Theta\) to \([0,1]\).
Define the left- and right-continuous limits \(f_Z^{(2)\mathrm{L}}\) and \(f_Z^{(2)\mathrm{R}}\), respectively, of \(f_Z^{(2)}\), given by \[f_Z^{(2)\mathrm{L}}(z,x) := \lim_{y\uparrow x} f^{(2)}(z,y) \quad \text{and}\quad f_Z^{(2)\mathrm{R}}(z,x) := \lim_{y\downarrow x} f^{(2)}(z,y), \qquad \forall z\in \mathcal{Z},\;x\in[0,1].\] Because \(f^{(2)}_Z\) is deterministic on \(\{0,1\}\), \(f_Z^{(2)\mathrm{L}}\) and \(f_Z^{(2)\mathrm{R}}\) are equivalent with \(c_i(f^{(2)\mathrm{L}}) = c_i(f^{(2)\mathrm{R}})\) for all \(i\). Furthermore, because pointwise limits preserve non-strict inequalities, we trivially have for all \(i\): \[\begin{align} c_i(f^{(2)\mathrm{R}}) & \le c_i(f^{(2)}), \quad \text{and} \\ \sup_{x\in [0,1]}\left| f_Z^{(2)\mathrm{R}}(z,x)- \mathbb{E}\left[f^{(2)\mathrm{R}}_Z(x)\right]\right| & \le \sup_{x\in [0,1]}\left| f_Z^{(2)}(z,x)-\mathbb{E}\left[f_Z^{(2)}(x)\right]\right|, \quad \forall z\in \mathcal{Z}. \end{align}\] Identical bounds hold for \(f_Z^{(2)\mathrm{L}}\).
Because the expected limits bound \(\vartheta\), exactly one of the following mutually exclusive cases must hold for any \(\vartheta \in [0,1]\): \[\begin{align} \mathbb{E}\left[f^{(2)}_Z(t_\vartheta)\right] & < \vartheta \le \mathbb{E}\left[f^{(2)\mathrm{R}}_Z(t_\vartheta)\right],\\ \mathbb{E}\left[f^{(2)\mathrm{L}}_Z(t_\vartheta)\right] & \le \vartheta < \mathbb{E}\left[f^{(2)}_Z(t_\vartheta)\right],\\ \mathbb{E}\left[f^{(2)}_Z(t_\vartheta)\right] & = \vartheta . \end{align}\] Thus, there exists a unique weight vector \((\alpha_{\vartheta}, \alpha_{\mathrm{R} \vartheta}, \alpha_{\mathrm{L} \vartheta}) \in \mathbb{R}^3_+\) satisfying the simplex conditions: \[\alpha_{\vartheta}+ \alpha_{\mathrm{R} \vartheta}+ \alpha_{\mathrm{L} \vartheta}=1, \qquad \alpha_{\mathrm{L} \vartheta} \alpha_{\mathrm{R} \vartheta}=0,\] that interpolates the target expectation: \[\vartheta = \alpha_{\vartheta} \mathbb{E}\left[f^{(2)}_Z(t_\vartheta)\right] + \alpha_{\mathrm{R} \vartheta} \mathbb{E}\left[f^{(2)\mathrm{R}}_Z(t_\vartheta)\right] + \alpha_{\mathrm{L} \vartheta} \mathbb{E}\left[f^{(2)\mathrm{L}}_Z(t_\vartheta)\right] .\] Clearly \(\alpha_\vartheta=1\) if and only if \(\vartheta = \mathbb{E}\left[f^{(2)}_Z(t_\vartheta)\right]\). Also, in case (R), \(\alpha_{\mathrm{L} \vartheta}=0\); in case (L) \(\alpha_{\mathrm{R} \vartheta}=0\).
Define the target randomised function \(g_Z\) in \(\mathcal{M}_{\mathsf{id}}([0,1])\) via \[g(z,\vartheta) := \alpha_{\vartheta} f^{(2)}(z,t_\vartheta) + \alpha_{\mathrm{R} \vartheta} f^{(2)\mathrm{R}}(z,t_\vartheta) + \alpha_{\mathrm{L} \vartheta} f^{(2)\mathrm{L}}(z,t_\vartheta).\] By Lemma 33 of the Appendix, the map \(\vartheta \mapsto g(z, \vartheta)\) is non-decreasing for each \(z \in \mathcal{Z}\). Linearity of expectation implies \(\mathbb{E}\left[g_Z(\vartheta)\right] = \vartheta\) for all \(\vartheta \in [0,1]\), establishing \(g_Z \in \mathcal{M}_{\mathsf{id}}([0,1])\).
Applying Lemma 21, \(g_Z\) is equivalent to \(f^{(3)}_Z\). Because \(f^{(3)}_Z\) is equivalent to \(f^{(2)}_Z\), \[\eta_g = \eta_{f^{(2)}} = (1+a)^2 \eta_f\,.\] Equivalence is transitive; \(g_Z\) is equivalent to \(f_Z\). ◻
Proof of Theorem 3. By Lemma 22, \(U_{\mathcal{M}_{\mathsf{id}}}\) is well defined on the whole of \(\mathbb{R}_+\), because \(U_\mathcal{M}\) is well defined. Further, \[U_{\mathcal{M}_{\mathsf{id}}}(x) \le U_{\mathcal{M}}(x) \le U_{\mathcal{M}_{\mathsf{id}}}\big((1+a)^2 x\big)\,,\] for any \(x\ge 0\) and any \(a>0\). Because \(a\) is arbitrary, \(U_{\mathcal{M}_{\mathsf{id}}}\) is monotone and its right-continuous version equals the right-continuous version of \(U_{\mathcal{M}}\). ◻
This section establishes the following theorem.
Theorem 4. For every \(x> 1\), \[\left( \sqrt{U_{\mathcal{M}}(x)} - \sqrt{\log_4 x} \right)^+ \le 2\min\left\{1,\frac{\ln(\mathrm{e}+\ln x)}{\sqrt{\ln x}}\right\} \,.\]
Because \(U_{\mathcal{M}}\) is monotone, its right-continuous version dominates it, and by Theorem 3 it, in turn, equals the right-continuous version of \(U_{\mathcal{M}_{\mathsf{id}}}\). As the envelope on the right of the theorem is continuous, it suffices to establish the bound for \(U_{\mathcal{M}_{\mathsf{id}}}\). We begin with a geometric property of deterministic non-decreasing functions.
The deviation of a non-decreasing function \(f\colon[0,1]\to[0,1]\) from the identity satisfies the following property, whose straightforward proof is deferred to the Appendix.
Lemma 23 (Ramp property). Let \(f\colon[0,1]\to[0,1]\) be a non-decreasing function.
If \(\sup_{t\in[0,1]} (f(t)-t)^+ \ge \varepsilon\), then there exists \(t^{*}\in[0,1]\) such that for all \(0<\alpha\le \varepsilon\), \[(f(t^{*}+\alpha)-(t^{*}+\alpha))^+ \ge \varepsilon-\alpha \,.\]
If \(\sup_{t\in[0,1]} (f(t)-t)^- \ge \varepsilon\), then there exists \(t^{*}\in[0,1]\) such that for all \(0<\alpha\le \varepsilon\), \[(f(t^{*}-\alpha)-(t^{*}-\alpha))^- \ge \varepsilon-\alpha \,.\]
A pointwise deviation induces a forward-propagating ramp of deviations, bounding a triangular region adjacent to the identity (Figure 1).
Let \(\mathcal{H}\) be the cone of \(C^1\), strictly convex, strictly increasing Young functions \(H\colon\mathbb{R}_+\to\mathbb{R}_+\) satisfying \(H(0)=0\). For \(H\in\mathcal{H}\), let \(h:=H'\) denote its strictly increasing gauge.
Corollary 2. If \(f\colon[0,1]\to[0,1]\) is a deterministic non-decreasing function and \(r>0\), then \[r \sup_{t\in[0,1]} |f(t)-t| \;\le\; \inf_{H\in\mathcal{H}} H^{-1}\!\left(r\int_0^1 h\!\left(r\,|f(t)-t|\right)\,\mathrm{d}t\right) \,.\]
Proof. Assume \(m := r \sup_{t\in[0,1]} |f(t)-t| > 0\); the case \(m=0\) is trivial. By definition either \(\sup_t (f(t)-t)^+ = m/r\) or \(\sup_t (f(t)-t)^- = m/r\); consider the former, the latter being symmetric. By Lemma 23 there is \(t^\ast\) with \(t^\ast+\alpha\in[0,1]\) for \(0\le\alpha\le m/r\) and \[r\bigl(f(t)-t\bigr)\;\ge\; m-r\,(t-t^\ast), \qquad t\in[t^\ast,\,t^\ast+m/r] \,.\] Fix \(H\in\mathcal{H}\) with gauge \(h=H'\ge 0\) strictly increasing. Restricting the integral to this interval and using monotonicity of \(h\), \[\begin{align} r\int_0^1 h\!\left(r|f(t)-t|\right)\mathrm{d}t & \ge r\int_{t^\ast}^{t^\ast+m/r} h\!\bigl(m-r(t-t^\ast)\bigr)\mathrm{d}t \\ & = \int_0^m h(s)\,\mathrm{d}s = H(m) \,, \end{align}\] by the substitution \(s=m-r(t-t^\ast)\). Because \(H^{-1}\) is increasing, \[H^{-1}\!\left(r\int_0^1 h\!\left(r|f(t)-t|\right)\mathrm{d}t\right)\ge m \,.\] As \(H\in\mathcal{H}\) was arbitrary, the infimum is \(\ge m\). ◻
We record the following corollary, whose proof is now immediate.
Corollary 3. If \(f_Z\in\mathcal{M}_{\mathsf{id}}([0,1])\) and \(0<\eta_f<\infty\), then for every realisation \(z\) of \(Z\), \[\sqrt{\eta_f}\,\sup_{t\in[0,1]} |f_z(t)-t| \;\le\; \inf_{H\in\mathcal{H}} H^{-1}\!\left(\sqrt{\eta_f}\int_0^1 h\!\left(\sqrt{\eta_f}\,|f_z(t)-t|\right)\mathrm{d}t\right) \,.\]
Proof. Fix a realisation \(z\). We have \(f_z\colon[0,1]\to[0,1]\) is deterministic and non-decreasing, so Corollary 2 applies with \(r=\sqrt{\eta_f}\) and \(f=f_z\). ◻
For any \(f_Z \in \mathcal{M}_{\mathsf{id}}\), \(\mathbb{E}\left[f_Z(t)\right] = t\). By McDiarmid’s inequality, the pointwise absolute deviation satisfies the sub-Gaussian tail: \[\sup_{t\in[0,1]} \mathbb{P}\left\{ \sqrt{\eta_f}|f_Z(t)-t| \ge \varepsilon \right\} \le 2\mathrm{e}^{-2\varepsilon^2} \,, \qquad \varepsilon > 0 \,.\] Let \(\mathcal{H}_0 \subset \mathcal{H}\) denote the subclass of gauges satisfying the growth condition \[\lim_{\varepsilon \to \infty} h(\varepsilon)\mathrm{e}^{-2\varepsilon^2} = 0 \,.\] Corollary 3 provides a bound for the expected absolute supremum.
Lemma 24. For any \(f_Z \in \mathcal{M}_{\mathsf{id}}\), we have \[\mathbb{E}\left[ \sqrt{\eta_f} \sup_{t\in[0,1]} |f_Z(t)-t| \right] \le \inf_{H\in\mathcal{H}_0} H^{-1}\left( \sqrt{\eta_f} \int_0^\infty 8\varepsilon\, h(\varepsilon) \mathrm{e}^{-2\varepsilon^2} \,\mathrm{d}\varepsilon \right) \,.\]
Proof. Set \(r := \sqrt{\eta_f}\) and define the random variable \(X_t := r|f_Z(t)-t|\) for \(t\in[0,1]\), with realisation \(X_t(z) = r|f_z(t)-t|\). By McDiarmid’s inequality, \(\mathbb{P}\{X_t \ge \varepsilon\} \le 2\mathrm{e}^{-2\varepsilon^2}\) for all \(\varepsilon>0\). By Corollary 3, for each realisation \(z\) and any \(H \in \mathcal{H}_0\), \[\sup_{t\in[0,1]}X_t(z) \le H^{-1}\left( r\int_0^1 h(X_t(z))\,\mathrm{d}t \right) \,.\] Taking expectations (the left side being independent of \(H\)), then applying Jensen’s inequality to the strictly concave \(H^{-1}\) and Tonelli’s theorem to the inner integral, we have \[\begin{align} \mathbb{E}\left[\sup_{t\in[0,1]}X_t\right] & \le \inf_{H\in\mathcal{H}_0} \mathbb{E}\left[ H^{-1}\left( r\int_0^1 h(X_t)\,\mathrm{d}t \right)\right] \\ & \le \inf_{H\in\mathcal{H}_0} H^{-1}\left( r\int_0^1 \mathbb{E}\left[h(X_t)\right]\,\mathrm{d}t \right) \,. \end{align}\] For fixed \(t\), the layer-cake representation with respect to the Lebesgue–Stieltjes measure \(\mathrm{d}h\) provides \[\mathbb{E}\left[h(X_t)\right] = h(0)+\int_0^\infty \mathbb{P}\{X_t\ge \varepsilon\}\,\mathrm{d} h(\varepsilon) \le h(0)+\int_0^\infty 2\mathrm{e}^{-2\varepsilon^2}\,\mathrm{d} h(\varepsilon) \,.\] By integration by parts, \[\begin{align} \int_0^\infty 2\mathrm{e}^{-2\varepsilon^2}\,\mathrm{d}h(\varepsilon) & = \left[2h(\varepsilon)\mathrm{e}^{-2\varepsilon^2}\right]_{0}^{\infty} - \int_0^\infty h(\varepsilon)\,\mathrm{d}\left(2\mathrm{e}^{-2\varepsilon^2}\right) \\ & = -2h(0) + \int_0^\infty 8\varepsilon\,h(\varepsilon)\mathrm{e}^{-2\varepsilon^2}\,\mathrm{d}\varepsilon \,. \end{align}\] Because \(H \in \mathcal{H}_0\), the upper boundary term is zero. Combining the residual terms, \(h(0) - 2h(0) = -h(0) \le 0\). Upper-bounding this non-positive term by zero, we have \[\mathbb{E}\left[h(X_t)\right] \le \int_0^\infty 8\varepsilon\,h(\varepsilon)\mathrm{e}^{-2\varepsilon^2}\,\mathrm{d}\varepsilon \,,\] uniformly in \(t\). Substituting this bound back into the integral inside \(H^{-1}\) completes the proof. ◻
For each \(x>0\), let \[\beta(x) := \inf_{H\in\mathcal{H}_0} H^{-1}\left( \sqrt{x} \int_0^\infty 8\varepsilon\, h(\varepsilon) \mathrm{e}^{-2\varepsilon^2} \,\mathrm{d}\varepsilon \right)\,.\] By Lemma 24, any \(f_Z \in \mathcal{M}_{\mathsf{id}}\) with finite \(\eta_f>0\) satisfies \[0 \le \sqrt{\eta_f}\, \mathbb{E}\left[ \sup_{t\in[0,1]} |f_Z(t)-t| \right] \le \beta(\eta_f)\,.\] Define the envelope \[\mathcal{E}(x) := \sqrt{\frac{2}{\ln 2}\beta(x)^2 +1} \,,\] for \(x>0\). We obtain the following bound.
Lemma 25. For every \(x>0\), \[\sqrt{U_{\mathcal{M}_{\mathsf{id}}}(x)} \le \mathcal{E}(x)\,.\]
Proof. Fix \(f_Z\in\mathcal{M}_{\mathsf{id}}\) and set \(x=\eta_f\). We establish that for all \(\varepsilon>0\), \[\mathbb{P}\left\{ \sqrt{\eta_f}\sup_{t\in[0,1]}|f_Z(t)-t| \ge \varepsilon\mathcal{E}(\eta_f) \right\} \le 2\mathrm{e}^{-2\varepsilon^2}.\] If \(\varepsilon\le \sqrt{(\ln 2)/2}\), then \(2\mathrm{e}^{-2\varepsilon^2}\ge 1\), and the inequality holds trivially. Assume \(\varepsilon>\sqrt{(\ln 2)/2}\) and set \(a := \sqrt{(\ln 2)/2}\). The \(\ell_2\)-norm identities are \[\mathcal{E}(\eta_f) = \left\| \left(1, \frac{\beta(\eta_f)}{a}\right) \right\|_2 \,, \qquad \varepsilon = \left\| \left(\sqrt{\varepsilon^2-a^2}, a\right) \right\|_2 \,.\] By the Cauchy–Schwarz inequality, \[\begin{align} \varepsilon\mathcal{E}(\eta_f) & \ge \left\langle \left(\sqrt{\varepsilon^2-a^2}, a\right), \left(1, \frac{\beta(\eta_f)}{a}\right) \right\rangle \\ & = \sqrt{\varepsilon^2-a^2} + \beta(\eta_f) \,. \end{align}\] Define the maximal deviation \[S(Z) := \sup_{t\in[0,1]}|f_Z(t)-t| \,.\] Applying Corollary 6 with the deterministic shift \(g(t) = t\) provides the centred upper tail concentration \[\mathbb{P}\left\{ \sqrt{\eta_f}\left( S(Z) - \mathbb{E}\left[S(Z)\right] \right) \ge v \right\} \le \mathrm{e}^{-2v^2}, \qquad v>0.\] Substituting the expectation bound \(\sqrt{\eta_f}\mathbb{E}\left[S(Z)\right] \le \beta(\eta_f)\) of Lemma 24, we have \[\mathbb{P}\left\{ \sqrt{\eta_f}S(Z) \ge \beta(\eta_f)+v \right\} \le \mathrm{e}^{-2v^2}, \qquad v>0.\] Evaluating at \(v = \sqrt{\varepsilon^2-a^2}\), \[\begin{align} \mathbb{P}\left\{ \sqrt{\eta_f}S(Z) \ge \varepsilon\mathcal{E}(\eta_f) \right\} & \le \mathbb{P}\left\{ \sqrt{\eta_f}S(Z) \ge \beta(\eta_f)+\sqrt{\varepsilon^2-a^2} \right\} \\ & \le \mathrm{e}^{-2(\varepsilon^2-a^2)} = \mathrm{e}^{-2\varepsilon^2} \mathrm{e}^{2a^2} = 2\mathrm{e}^{-2\varepsilon^2}, \end{align}\] where \(\mathrm{e}^{2a^2}=\mathrm{e}^{\ln2}=2\). Hence \(\mathcal{E}(x)^2\) is admissible, in the sense defined in Section 2, for every \(f_Z\in\mathcal{M}_{\mathsf{id}}\) with \(\eta_f=x\), so \(\mathcal{E}(x)^2\ge u_{f_Z}\); taking the supremum, \(\sqrt{U_{\mathcal{M}_{\mathsf{id}}}(x)}\le\mathcal{E}(x)\). ◻
We extract an explicit analytical envelope from \(\beta(x)\). For \(p>0\), define the Young function \(H_p(y) = \mathrm{e}^{py} - 1\), with gauge \(h_p(\varepsilon) = p \mathrm{e}^{p\varepsilon}\) and inverse \(H_p^{-1}(z) = \frac{1}{p}\ln(1+z)\). Because \(p\mathrm{e}^{p\varepsilon} = o(\mathrm{e}^{2\varepsilon^2})\), we have \(H_p \in \mathcal{H}_0\). We establish the following envelope result.
Proposition 26. For any \(p>0\) and \(x>0\), the absolute expected supremum envelope satisfies: \[\beta(x) \le \psi(x,p) := \frac{1}{p} \ln\left( 1 + 2p\sqrt{x} + p^2 \sqrt{2\pi x} \, \mathrm{e}^{p^2/8} \right) .\] In particular, for all \(x, p>0\), \[\sqrt{U_{\mathcal{M}}(x)} \le \mathcal{E}(x,p) := \sqrt{\frac{2}{\ln 2} \psi(x,p)^2 + 1}\,.\]
Proof. We upper-bound the infimum in the definition of \(\beta(x)\) by evaluating the exact integral for the specific gauge \(h_p(\varepsilon) = p \mathrm{e}^{p\varepsilon}\): \[I_p := \int_0^\infty 8\varepsilon\, h_p(\varepsilon) \mathrm{e}^{-2\varepsilon^2} \,\mathrm{d}\varepsilon = \int_0^\infty 8\varepsilon p \mathrm{e}^{p\varepsilon - 2\varepsilon^2} \,\mathrm{d}\varepsilon \,.\] To evaluate this, we apply integration by parts. Let \(u = 2p \mathrm{e}^{p\varepsilon}\) and \(\mathrm{d} v = 4\varepsilon \mathrm{e}^{-2\varepsilon^2}\,\mathrm{d}\varepsilon\), which integrates to \(v = -\mathrm{e}^{-2\varepsilon^2}\). This implies: \[\begin{align} I_p & = \left[ -2p \mathrm{e}^{p\varepsilon - 2\varepsilon^2} \right]_0^\infty - \int_0^\infty \left(2p^2 \mathrm{e}^{p\varepsilon}\right) \left(-\mathrm{e}^{-2\varepsilon^2}\right) \,\mathrm{d}\varepsilon \\ & = 2p + 2p^2 \int_0^\infty \mathrm{e}^{p\varepsilon - 2\varepsilon^2} \,\mathrm{d}\varepsilon \,. \end{align}\] We isolate the Gaussian core by completing the square in the exponent: \(-2\varepsilon^2 + p\varepsilon = -2\left(\varepsilon - \frac{p}{4}\right)^2 + \frac{p^2}{8}\). Substituting this back into the integral and factoring out the exponential, we have \[I_p = 2p + 2p^2 \mathrm{e}^{p^2/8} \int_0^\infty \mathrm{e}^{-2(\varepsilon - p/4)^2} \,\mathrm{d}\varepsilon \,.\] We upper-bound the remaining shifted Gaussian integral by extending its lower limit to \(-\infty\), which evaluates exactly to \(\sqrt{\pi/2}\). Substituting this, we have \[I_p \le 2p + 2p^2 \sqrt{\frac{\pi}{2}} \mathrm{e}^{p^2/8} = 2p + p^2\sqrt{2\pi} \, \mathrm{e}^{p^2/8} \,.\] Because \(\beta(x)\) is defined as an infimum, it is bounded by the evaluation at \(H_p\), so that \(\beta(x) \le H_p^{-1}(\sqrt{x} I_p) = \frac{1}{p} \ln(1 + \sqrt{x} I_p)\).
The “in particular” claim follows because \(\sqrt{U_{\mathcal{M}_{\mathsf{id}}}(x)}\le\mathcal{E}(x)\le\mathcal{E}(x,p)\) by Lemma 25, and the uniformisation Theorem 3 transports the bound to \(U_{\mathcal{M}}\), the envelope \(\mathcal{E}(\cdot,p)\) being continuous. ◻
Tuning \(p(x)=2\sqrt{\ln x}\) in Proposition 26, and writing \(y=\ln x\), the choice \(p=2\sqrt{y}\) has \(p^2/8=y/2\), so that \(\psi(e^y,2\sqrt{y})=\tfrac{\sqrt{y}}{2}\bigl(1+\Lambda(y)\bigr)\) and hence \[\label{eq:lemasymp} \sqrt{U_{\mathcal{M}}(e^y)}\le \mathcal{E}(e^y,2\sqrt{y}) = \sqrt{\frac{y}{2\ln 2}}\, \sqrt{(1+\Lambda(y))^2 + \frac{2\ln 2}{y}}\,,\tag{5}\] where \[\Lambda(y) := \frac{1}{y} \ln\left( 4\sqrt{2\pi}\,y + 4\sqrt{y}\,\mathrm{e}^{-y/2} + \mathrm{e}^{-y} \right) \,.\]
The proof of the following computational lemma is deferred to the appendix.
Lemma 27. With \(\Lambda\) as in 5 , the bound \[\sqrt{\frac{y}{2\ln 2}} \left[ \sqrt{ (1+\Lambda(y))^2 + \frac{2\ln 2}{y} } - 1 \right] \le \frac{2\ln(\mathrm{e}+y)}{\sqrt{y}}\] holds for all \(y \ge y_0\), where \(y_0\) is the unique solution of \(\ln(\mathrm{e}+y)/\sqrt{y} = 1\), \(y \ge 1\).
We obtain the needed bound after threshold, for \(y\ge y_0\).
Corollary 4. For every \(x\) with \(\ln x\ge y_0\), \[\sqrt{U_{\mathcal{M}}(x)}-\sqrt{\log_4 x} \le \frac{2\ln(\mathrm{e}+\ln x)}{\sqrt{\ln x}}\,.\]
Proof. Write \(y=\ln x\), so \(\sqrt{\log_4 x}=\sqrt{y/(2\ln2)}\). Subtracting this leading term from 5 , we have \[\sqrt{U_{\mathcal{M}}(e^y)}-\sqrt{\tfrac{y}{2\ln2}} \le \sqrt{\tfrac{y}{2\ln2}} \left[\sqrt{(1+\Lambda(y))^2+\tfrac{2\ln2}{y}}-1\right],\] which by Lemma 27 is at most \(2\ln(\mathrm{e}+y)/\sqrt{y}\) for \(y\ge y_0\). ◻
For small \(x\), let \(p=5\) in Proposition 26, chosen heuristically. The substitution provides \(\psi(x,5)=\tfrac15\ln(1+K\sqrt{x})\) with \(K=10+25\sqrt{2\pi}\,\mathrm e^{25/8}\), hence \[\label{eq:lemsmall} \sqrt{U_{\mathcal{M}}(x)}\le \mathcal{E}(x,5) = \sqrt{\frac{2}{\ln 2}\,\psi(x,5)^2 + 1}\,.\tag{6}\]
The following is a direct calculation over the range \([2.5, e^{y_0})\).
Lemma 28. For every \(x\) with \(2.5\le x\) and \(\ln x<y_0\), \[\sqrt{U_{\mathcal{M}}(x)}-\sqrt{\log_4 x} < 2 \,.\]
Proof. By 6 , \[\sqrt{U_{\mathcal{M}}(x)}-\sqrt{\log_4 x} \le \mathcal{E}(x,5)-\sqrt{\log_4 x}=:S(x) \,.\] Because \(\mathcal{E}(x,5)\ge\sqrt{\tfrac{2}{\ln2}}\,\psi(x,5)\), on \(2.5\le x<\mathrm{e}^{y_0}\) \[\frac{\mathrm d}{\mathrm dx}\,\mathcal{E}(x,5) = \frac{(2/\ln2)\,\psi(x,5)\,\psi'(x)}{\mathcal{E}(x,5)} \le \sqrt{\tfrac{2}{\ln2}}\,\psi'(x) < \frac{\mathrm d}{\mathrm dx}\sqrt{\log_4 x}\,,\] so \(S\) is decreasing and attains its maximum at \(x=2.5\): \[S(2.5)=\mathcal{E}(2.5,5)-\sqrt{\log_4 2.5}=1.9966\ldots<2 \,.\] This establishes the result. ◻
We also know the trivial bound of Lemma 4, which holds for all classes of \([0,1]\)-valued randomised functions: \[\sqrt{U_\mathcal{M}(x)} \le \sqrt{\frac{2}{\ln2}\,x}, \qquad \text{for all } x\ge 0\,.\] It is the sharpest of the three on the micro-\(x\) range \([1,2.5)\).
Lemma 29. For every \(x\) with \(1\le x < 2.5\), \[\sqrt{U_{\mathcal{M}}(x)}-\sqrt{\log_4 x} < 2 \,.\]
Proof. By the trivial bound of Lemma 4, \(\sqrt{U_{\mathcal{M}}(x)}-\sqrt{\log_4 x}\le T(x):= \sqrt{\tfrac{2}{\ln2}\,x}-\sqrt{\log_4 x}\). A direct verification provides \(T(x)\le T(2.5)=1.8728\ldots<2\) on \([1,2.5]\). ◻
Proof of Theorem 4. Put \(y=\ln x\), so that \(\sqrt{\log_4 x}=\sqrt{y/(2\ln2)}\), and let \(x_0=\mathrm e^{y_0}\). Because the right-hand side of the asserted bound is nonnegative, it suffices to show \[\sqrt{U_{\mathcal{M}}(x)}-\sqrt{\log_4 x} \le 2\min\!\left\{1,\frac{\ln(\mathrm{e}+y)}{\sqrt{y}}\right\}\] for all \(x>1\); the positive part then follows at once. The function \(g(y):=\ln(\mathrm{e}+y)/\sqrt{y}\) is strictly decreasing on \((0,\infty)\), because \(g'(y)<0\) is equivalent to the elementary inequality \(2y/(\mathrm{e}+y)<\ln(\mathrm{e}+y)\); and \(g(y_0)=1\), so \(g\le1\) for \(y\ge y_0\) and \(g\ge1\) for \(0<y\le y_0\). Hence \[2\min\!\left\{1,\frac{\ln(\mathrm{e}+y)}{\sqrt{y}}\right\} =\begin{cases} \dfrac{2\ln(\mathrm{e}+y)}{\sqrt{y}}, & y\ge y_0 \quad(x\ge x_0), \\[8pt] 2, & 0<y< y_0 \quad(1<x<x_0). \end{cases}\] Case \(x\ge x_0\). By Corollary 4, \[\sqrt{U_{\mathcal{M}}(x)}-\sqrt{\log_4 x} \le \frac{2\ln(\mathrm{e}+\ln x)}{\sqrt{\ln x}} = 2\min\!\left\{1,\frac{\ln(\mathrm{e}+y)}{\sqrt{y}}\right\}.\] Case \(1<x<x_0\). Here the right-hand side is \(2\). The ranges \([1,2.5)\) and \([2.5,x_0)\) partition \((1,x_0)\): on \([1,2.5)\), Lemma 29 provides \(\sqrt{U_{\mathcal{M}}(x)}-\sqrt{\log_4 x}<2\), and on \([2.5,x_0)\), Lemma 28 provides the same bound. Thus \(\sqrt{U_{\mathcal{M}}(x)}-\sqrt{\log_4 x}<2\) throughout \((1,x_0)\). Combining the two cases proves the bound for every \(x>1\). ◻
We open with the empirical distribution function that motivated the present study: a stratified block construction that asymptotically saturates the Main Theorem and so exhibits the lower bound (Theorem 2) as sharp. We then give two applications of the upper bound (Theorem 4) to non-monotone bounded processes: an analogue of the smooth Gaussian setting of Diebolt and Posse [11] with the \(\ell_1\) arc length in place of the \(\ell_2\) arc length used there.
The methods developed in the proof, which transform randomised functions, extend the scope of the main theorem beyond manifestly monotone processes. We do not pursue this here beyond the smooth processes example.
We present an empirical distribution function that asymptotically saturates the Main Theorem. It arose in the analysis of a randomised field experiment with clustered structure, the cluster sizes growing faster than the number of clusters, and motivated the present study. The classical uniform law of large numbers does not in general apply to such data (cf. [15]); the stratification below restores it.
Let \[X_{i,j}, \quad 1 \le i \le n, ~1 \le j \le m_n\,,\] be real-valued variables, the blocks \((X_{i,1}, \dots, X_{i,m_n})\) being independent across \(i\), and write \[\label{eq:F} F_{nm}(t) = \frac{1}{n m_n} \sum_{i=1}^n \sum_{j=1}^{m_n} \mathbf{1}_{\{X_{i,j} \le t\}}.\tag{7}\] Suppose the supports are disjoint and ordered within each block: \[\label{eq:stratification} \operatorname*{ess\,sup} X_{i,j} \le \operatorname*{ess\,inf} X_{i,j+1}.\tag{8}\] Any threshold \(t\) lies in at most one of the intervals \([\operatorname*{ess\,inf} X_{i,j}, \operatorname*{ess\,sup} X_{i,j}]\), and resampling the \(i\)-th block alters the inner sum \(\sum_{j} \mathbf{1}_{\{X_{i,j} \le t\}}\) by at most \(1\). Hence \(\eta_{F_{nm}} = n m_n^2\).
An instance of \(F_{nm}\) realises the lower-bound construction of Section 3. Take \(\mathcal{Z} = \{-1, +1\}^n\) and a nonempty \(\Theta \subseteq \mathcal{Z}\) with \(|\Theta| = m_n\). Arrange \(\Theta = [\theta_{ij}]_{i,j=1}^{n, m_n}\) as a matrix, with columns \(\theta_j\). Draw \(Z\) uniformly on \(\mathcal{Z}\) and set \[X_{i,j} = Z_i \theta_{ij} + 4j.\] The supports \(\{4j-1, 4j+1\}\) are disjoint and ordered, so 8 holds. For \(t \in [4j-1, 4j+1)\), \[F_{nm}(t) - \mathbb{E}\left[F_{nm}(t)\right] = -\frac{1}{2 n m_n} \langle Z, \theta_j \rangle,\] and the centred difference vanishes elsewhere. The right-hand side is the inner product \(f\) of 3 , scaled by \(-(2 n m_n)^{-1}\); by Lemma 2, \(F_{nm}\) is equivalent to \(f_Z\) on \(\Theta\), with \(\eta_{F_{nm}} = n m_n^2\) as before.
Consider the regime \[\frac{\ln n}{\ln m_n} \to 0,\] in which \(\ln n / \ln \eta_{F_{nm}} \to 0\) and \(\ln m_n / \ln \eta_{F_{nm}} \to 1/2\). The tuned constellation construction of Section 3.5, applied with \(\ln m_n \sim \tfrac{5}{2}\sqrt{n}\), gives \[\lim_{n \to \infty} \frac{u_{F_{nm}}}{\ln \eta_{F_{nm}}} = \frac{1}{2 \ln 2},\] the constant of the Main Theorem; with this growth we realise finite-sample lower bounds in the proof of Theorem 2.
The cluster growth regime \[\frac{\ln n}{\ln m_n} \to 0,\] attaining the optimal normalisation \(\sqrt{\eta_{F_{nm}}/\log_4 \eta_{F_{nm}}}\) of the main theorem, is natural in cluster sampling: for the Ewens sampling formula, for example, the \(k\)th largest cluster grows linearly while the number of clusters grows logarithmically.
The intermediate regime, \[\frac{\ln m_n}{\ln n} \to c \in (0,\infty),\] for which the bracketing normalisation \(\sqrt{n}\) is informative, is incompatible with cluster sampling models admitting a Kingman paintbox representation [16], [17] and is realised by nonexchangeable microclustering sampling [18].
The reverse regime, in which \(n\) outpaces \(m_n\), \[\frac{\ln m_n}{\ln n} \to 0,\] for which the bracketing normalisation \(\sqrt{n}\) is tight (Remark 30), is unattainable with microclustering sampling, and for exchangeable sampling is realised only in the trivial case \[m_n = 1,\] recovering Massart’s framework of scalar i.i.d.empirical distribution functions.
We turn to applications of the upper bound (Theorem 4), and begin with a class of smooth bounded processes whose suprema are controlled by the \(\ell_1\) arc length of the coordinate functions. This is the distribution-free counterpart of a classical Gaussian extreme-value setting, discussed in Remark 31.
Consider a randomised function \[X_Z(t) = \sum_{j=1}^m g_j(Z) \psi_j(t) \,, \qquad t \in [0, T] \,,\] where \(g_j \colon \mathcal{Z} \to [-1, 1]\) are measurable, and \(\psi_j \colon [0, T] \to \mathbb{R}\) are continuously differentiable with \(\sum_{j=1}^m |\psi_j(t)| = 1\) for all \(t \in [0, T]\).
Assume \(\mathbb{E}\left[X_Z(t)\right] = 0\) for all \(t \in [0, T]\). Define the \(\ell_1\) arc length \[\mathrm{Arc}_{\ell_1}(t) := \int_0^t \sum_{j=1}^m |\psi'_j(s)| \, \mathrm{d}s \,, \qquad \mathrm{Arc}_{\ell_1} := \mathrm{Arc}_{\ell_1}(T) \,.\]
Corollary 5. If \(1< \eta_X<\infty\), then for all \(\varepsilon > 0\), \[\mathbb{P}\!\left\{ \frac{\sqrt{\eta_X}}{R\bigl((2 + \mathrm{Arc}_{\ell_1})^2 \eta_X\bigr)} \sup_{t \in [0, T]} |X_Z(t)| \ge \varepsilon \right\} \le 2 \mathrm{e}^{-2\varepsilon^2}\,,\] where \[R(x) := \sqrt{\log_4 x} + 2\min\!\left\{1, \frac{\ln(\mathrm{e}+\ln x)}{\sqrt{\ln x}}\right\} \,, \qquad \text{for all x> 1} \,.\]
Proof. The bound \(|X_Z(t)| \le 1\) holds by the normalisation \(\sum_{j=1}^m |\psi_j(t)| = 1\), and the local-derivative bound \(\kappa(t) := \sum_{j=1}^m |\psi'_j(t)|\) satisfies \(|X'_Z(t)| \le \kappa(t)\) for every realisation. The drifted process \[Y_Z(t) := X_Z(t) + \mathrm{Arc}_{\ell_1}(t) + 1\] has derivative \(Y'_Z(t) = X'_Z(t) + \kappa(t) \ge 0\), so for every realisation \(Z = z\) the path \(Y_z\) is non-decreasing in \(t\), and \(Y_Z(t) \in [0,\, 2 + \mathrm{Arc}_{\ell_1}]\). Normalise to define \[f_Z(t) := \frac{Y_Z(t)}{2 + \mathrm{Arc}_{\ell_1}} \,,\] giving \(f_Z \in \mathcal{M}\). Applying Lemma 2 twice — first with \(a = 1\) and deterministic shift \(h(t) = \mathrm{Arc}_{\ell_1}(t) + 1\), then with \(a = (2 + \mathrm{Arc}_{\ell_1})^{-1}\) and \(h \equiv 0\) — shows that \(f_Z\) is equivalent to \(X_Z\) with \(\eta_{f_Z} = (2 + \mathrm{Arc}_{\ell_1})^2 \eta_X>1\), so Theorem 4 applied to \(f_Z\) gives the result. ◻
Remark 30 ().
The cluster empirical distribution \(F_{nm}\) is an average of \(n\) independent monotone summands. For the identically distributed case, the bracketing estimates [10] normalise by \(\sqrt{n}\), the number of independent summands, ensuring \[\mathbb{E}\left[\sup_t \left|F_{nm}(t) - \mathbb{E}\left[F_{nm}(t)\right]\right|\right] \le C/\sqrt{n}\,,\] for a universal constant \(C\). The optimal normalisation of Theorem 1 instead uses the effective sample size \(\eta_{F_{nm}} = n m_n^2\), normalising by \(\sqrt{\eta_{F_{nm}}/U_{\mathcal{M}}(\eta_{F_{nm}})} \sim \sqrt{\eta_{F_{nm}}/\log_4 \eta_{F_{nm}}}\). The bracketing normalisation is therefore loose by the factor \[\frac{\sqrt{\eta_{F_{nm}}/\log_4 \eta_{F_{nm}}}}{\sqrt{n}} = \frac{m_n}{\sqrt{\log_4 \eta_{F_{nm}}}} \sim \frac{m_n}{\sqrt{2 \log_4 m_n}}\] in the regime \(\ln n / \ln m_n \to 0\). In a regime realising the lower-bound construction (Section 3), \(\ln m_n \sim \tfrac{5}{2}\sqrt{n}\), so bracketing’s \(\sqrt{n}\) normalisation is exponentially loose at the scale of the effective sample size \(\sqrt{\eta_{F_{nm}}} = \sqrt{n}\, m_n\). \(\diamond\)
Remark 31 ().
Corollary 5 is the distribution-free analogue of the smooth Gaussian setting of Diebolt and Posse [11], who consider the process \[X_\xi(t) = \tau_X^{-1}(t) \sum_{j=1}^n \xi_j\, g_j(t) \,, \qquad t \in [0, T] \,,\] with
independent \(\xi_j \sim \mathcal{N}(0, 1)\), unit-sphere basis \(\sum_{j=1}^n g_j^2(t) \equiv 1\), and inverse standard deviation \(\tau_X(t) \ge 1\), so
that \(\operatorname{Var}\big(X_\xi(t)\big) = \tau_X^{-2}(t) \le 1\). Both \(\tau_X\) and the \(g_j\) are taken of class \(C^3\). The mapping \[t \mapsto g(t) := (g_1(t), \dots, g_n(t)) \in \mathbb{S}^{n-1}_2\] traces a curve \(g([0,T])\) in general position on the unit sphere \(\mathbb{S}^{n-1}_2\) of \(\mathbb{R}^n\), with \(\ell_2\) arc length \[\mathrm{Arc}_{\ell_2} := \int_0^T \sqrt{\sum_{j=1}^n
g'_j(t)^2}\, \mathrm{d}t \,.\] Diebolt and Posse obtain matching non-asymptotic upper and lower bounds on the density of \(\sup_{t \in [0, T]} X_\xi(t)\) via the differential geometry of the level manifolds of
the supremum functional; in the unit-variance stationary case \(\tau_X \equiv 1\) these integrate to the classical tail \[\mathbb{P}\!\left\{ \sup_{t \in [0, T]} X_\xi(t) > \varepsilon \right\}
\sim \frac{\mathrm{Arc}_{\ell_2}}{2\pi}\, \mathrm{e}^{-\varepsilon^2/2} \,, \qquad \varepsilon \to \infty \,.\] The maximum variance along the path, \(\sup_{t} \tau_X^{-2}(t)\), plays the role for \(X_\xi\) that the variance proxy \(\eta_X^{-1}\) plays for the bounded process \(X_Z\). Setting the tail probability to order one locates the typical maximum,
which scales as \[\sqrt{\eta_X^{-1}\ln \mathrm{Arc}_{\ell_1}} \quad \text{for } X_Z, \qquad\text{and}\qquad \sqrt{\big(\sup_{t} \tau_X^{-2}(t)\big)\,\ln \mathrm{Arc}_{\ell_2}} \quad \text{for } X_\xi,\] sharing the leading
\(\sqrt{\ln \mathrm{Arc}}\) factor. The correspondence is term-by-term: \[\begin{array}{ccc} X_Z \text{ (distribution-free)} & \longleftrightarrow & X_\xi \text{ (Gaussian)} \\[4pt]
\mathrm{Arc}_{\ell_1} & \longleftrightarrow & \mathrm{Arc}_{\ell_2} \\[2pt] \eta_X^{-1} & \longleftrightarrow & \sup_{t} \tau_X^{-2}(t) \\[2pt] \log_4 & \longleftrightarrow & 2\pi \end{array}\] the arc length, the worst-case
variance, and the tail constant respectively; only the last, the sharp logarithmic penalty in our main theorem, reflects the loss of Gaussianity. We do not pursue the exact connection between our non-asymptotic bounds and the density bounds of Diebolt and
Posse [11]. \(\diamond\)
To prove the required binomial lower bound for a constellation, we use the bounds \[\label{eq:binbound} \frac{2^{n\mathcal{H}_2(\lambda)}}{\sqrt{8n\lambda(1-\lambda)}} \le \binom{n}{\lambda n} \le \frac{2^{n\mathcal{H}_2(\lambda)}}{\sqrt{2\pi n\lambda(1-\lambda)}}, \qquad 0<\lambda<1,~ n\lambda\in\mathbb{N},~ n\in\mathbb{N} \,,\tag{9}\] where \(\mathcal{H}_2 \colon [0,1] \to [0,1]\), \[\mathcal{H}_2(x) := -x\log_2 x - (1-x)\log_2(1-x)\,,\] is the binary entropy function [13]. This function is symmetric about its maximiser \(x=\tfrac12\) and the exact Taylor expansion about \(\tfrac{1}{2}\) has no odd terms, \[\label{eq:H295taylor} 1-\mathcal{H}_2\left(\frac{1}{2}-\delta\right) = \frac{1}{\ln 2} \sum_{j=1}^{\infty} \frac{(2\delta)^{2j}}{2j(2j-1)}, \qquad \text{for } |\delta| < 1/2 \,.\tag{10}\]
We are ready to prove the required bound.
Lemma 32. For a constellation \((k, t_k, n_k, r_k, m_k)\), we have \[\frac{1}{2^{n_k}}\binom{n_k}{r_k} \ge \frac{1}{\sqrt{2k}\,t_k} \exp\left(-\frac{k}{12(t_k^2-1)}\right) \mathrm{e}^{-k/2} \,.\]
Proof. For a constellation, by definition, we have \(k,t_k\in\mathbb{N}\) with \(t_k\ge2\), and \[n_k := k\,t_k^2, \qquad r_k := \frac{k\,t_k(t_k-1)}{2}\,.\] Because \(t_k\) and \(t_k-1\) are consecutive integers, their product \(t_k(t_k-1)\) is an even integer, ensuring that \(r_k\ge 1\) is an integer. We also clearly have that \(n_k\ge 4\) is an integer. Set \[\lambda_k := \frac{r_k}{n_k} = \frac{k\,t_k(t_k-1)}{2k\,t_k^2}= \frac{t_k-1}{2t_k} = \frac{1}{2} - \frac{1}{2t_k}\,.\] Note that \(1/4\le\lambda_k <1/2\), because \(t_k\ge 2\). Thus, we apply the lower bound of 9 : \[\binom{n_k}{r_k}= \binom{n_k}{\lambda_k n_k} \ge \frac{1}{\sqrt{8n_k\lambda_k(1-\lambda_k)}}\, 2^{n_k \mathcal{H}_2(\lambda_k)} \,.\]
Dividing by \(2^{n_k}\) gives \[\label{eq:brackterm} \frac{1}{2^{n_k}}\binom{n_k}{r_k} \ge \frac{1}{\sqrt{8n_k\lambda_k(1-\lambda_k)}} \, \left[2^{-n_k(1-\mathcal{H}_2(\lambda_k))}\right] \,.\tag{11}\]
Now \[\lambda_k(1-\lambda_k) = \left(\frac{1}{2}-\frac{1}{2t_k}\right) \left(\frac{1}{2}+\frac{1}{2t_k}\right) = \frac{1}{4}-\frac{1}{4t_k^2} \le \frac{1}{4} \,,\] and therefore \[\label{eq:firstest} \frac{1}{\sqrt{8n_k\lambda_k(1-\lambda_k)}} \ge \frac{1}{\sqrt{2n_k}} = \frac{1}{\sqrt{2k}\,t_k} \,.\tag{12}\]
It remains to bound the bracketed term in 11 . Writing \(\delta:=\frac{1}{2t_k}\) and noting that \(0<\delta \le 1/4\), we have \[\lambda_k = \frac{1}{2} - \delta\,.\] By the exact Taylor expansion in 10 , we have \[1-\mathcal{H}_2\left(\frac{1}{2}-\delta\right) = \frac{1}{\ln 2} \sum_{j=1}^{\infty} \frac{(2\delta)^{2j}}{2j(2j-1)}\,.\] Because \(2\delta=1/t_k\), this becomes \[1-\mathcal{H}_2(\lambda_k) = \frac{1}{\ln 2} \sum_{j=1}^{\infty} \frac{t_k^{-2j}}{2j(2j-1)} \,.\] Separating the first term and bounding the remainder by a geometric series, \[\begin{align} 1-\mathcal{H}_2(\lambda_k) & \le \frac{1}{\ln 2} \left[ \frac{1}{2t_k^2} + \frac{1}{12}\sum_{j=2}^{\infty}t_k^{-2j} \right] \\ & = \frac{1}{\ln 2} \left[ \frac{1}{2t_k^2} + \frac{1}{12t_k^4}\sum_{m=0}^{\infty}t_k^{-2m} \right] \,. \end{align}\] Hence \[1-\mathcal{H}_2(\lambda_k) \le \frac{1}{\ln 2} \left[ \frac{1}{2t_k^2} + \frac{1}{12t_k^2(t_k^2-1)} \right] \,,\] and multiplying by \(n_k\ln 2=k\,t_k^2\ln 2\) gives \[n_k(1-\mathcal{H}_2(\lambda_k))\ln 2 \le k\,t_k^2\left(\frac{1}{2t_k^2} + \frac{1}{12t_k^2(t_k^2-1)}\right) = \frac{k}{2} + \frac{k}{12(t_k^2-1)} \,.\] Thus \[\begin{align} 2^{-n_k(1-\mathcal{H}_2(\lambda_k))} & = \exp\left(-n_k(1-\mathcal{H}_2(\lambda_k))\ln 2\right) \\ & \ge \exp\left(-\frac{k}{2}-\frac{k}{12(t_k^2-1)}\right) = \exp\left(-\frac{k}{12(t_k^2-1)}\right)\mathrm{e}^{-k/2} \,. \end{align}\]
Substituting this and 12 into 11 , we obtain \[\frac{1}{2^{n_k}}\binom{n_k}{r_k} \ge \frac{1}{\sqrt{2k}\,t_k} \exp\left(-\frac{k}{12(t_k^2-1)}\right) \mathrm{e}^{-k/2} \,.\] This proves the claim. ◻
Lemma 33. Let \(g_Z(\vartheta)\) be the interpolated function constructed in the proof of Lemma 22. For each fixed \(z \in \mathcal{Z}\), the mapping \(\vartheta \mapsto g(z, \vartheta)\) is non-decreasing.
Proof. Let \(\vartheta_1 < \vartheta_2\). Because the expectation \(\mathbb{E}\left[f^{(2)}_Z(t)\right]\) is strictly increasing, the generalised inverse gives \(t_{\vartheta_1} \le t_{\vartheta_2}\). Two cases arise:
If \(t_{\vartheta_1} < t_{\vartheta_2}\), the monotonicity of \(t \mapsto f^{(2)}(z,t)\) bounds the maximum possible interpolated value at \(t_{\vartheta_1}\) by the minimum possible interpolated value at \(t_{\vartheta_2}\). Thus, \[g(z, \vartheta_1) \le f^{(2)\mathrm{R}}(z, t_{\vartheta_1}) \le f^{(2)\mathrm{L}}(z, t_{\vartheta_2}) \le g(z, \vartheta_2) \,.\]
If \(t_{\vartheta_1} = t_{\vartheta_2} = t^\ast\), the interpolation occurs across a single vertical jump. Because \(\vartheta_1 < \vartheta_2\), the unique interpolation weights for \(\vartheta_2\) shift mass from the left limit \(f^{(2)\mathrm{L}}(z, t^\ast)\) towards the evaluation \(f^{(2)}(z, t^\ast)\) and the right limit \(f^{(2)\mathrm{R}}(z, t^\ast)\). The pointwise ordering \(f^{(2)\mathrm{L}}(z, t^\ast) \le f^{(2)}(z, t^\ast) \le f^{(2)\mathrm{R}}(z, t^\ast)\) ensures this rightward weight shift gives \(g(z, \vartheta_1) \le g(z, \vartheta_2)\).
Thus, \(g_Z\) is non-decreasing. ◻
Let \(f \colon \mathcal{Z} \times \Theta \to \mathbb{R}\) be bounded, and define the lattice functions \(\overline{f^+}, \overline{f^-}, \overline{|f|} \colon \mathcal{Z} \to \mathbb{R}\) by \[\overline{f^+}(z) = \sup_{\theta\in\Theta} (f(z,\theta))^+, \qquad \overline{f^-}(z) = \sup_{\theta\in\Theta} (f(z,\theta))^-, \qquad \overline{|f|}(z) = \sup_{\theta\in\Theta} |f(z,\theta)|\,.\] Each is a bounded function on \(\mathcal{Z}\), hence a function on \(\mathcal{Z}\times\{\ast\}\) with singleton index set, and so has an effective sample size.
Lemma 34. The lattice functions have effective sample size at least that of \(f\); that is, \[\eta_{\overline{f^+}} \ge \eta_f, \qquad \eta_{\overline{f^-}} \ge \eta_f, \qquad \eta_{\overline{|f|}} \ge \eta_f \,.\]
Proof. For each coordinate \(i\), \[c_i(\overline{f^+}) \le c_i(f), \qquad c_i(\overline{f^-}) \le c_i(f), \qquad c_i(\overline{|f|}) \le c_i(f)\,,\] and the result follows. ◻
We use the following in the proof of Lemma 25 for the upper bound of the main theorem.
Corollary 6. Let \(f_Z \colon \Theta \to \mathbb{R}\) be a randomised function with finite \(\eta_f\), and let \(h \colon \Theta \to \mathbb{R}\) be a bounded deterministic function. If \[S(Z) := \sup_{\theta\in \Theta} |f_Z(\theta) - h(\theta)|\] is a measurable random variable, then \[\mathbb{P}\left\{ \sqrt{\eta_f}\left( S(Z) - \mathbb{E}\left[S(Z)\right] \right) \ge \varepsilon \right\} \le \mathrm{e}^{-2\varepsilon^2}, \qquad \text{for all } \varepsilon > 0 \,.\]
Proof. Write \(g(z,\theta) = f(z,\theta) - h(\theta)\). Because \(h\) is deterministic, \(c_i(g) = c_i(f)\) for all \(i\), so \(\eta_g = \eta_f\); and \(S = \overline{|g|}\), so Lemma 34 gives \(\eta_S \ge \eta_g = \eta_f\). If \(\eta_S = \infty\), then \(S(Z) = \mathbb{E}\left[S(Z)\right]\) surely and the bound holds. Otherwise \(\infty > \eta_S \ge \eta_f > 0\), and for all \(\varepsilon > 0\) McDiarmid’s inequality gives \[\mathbb{P}\left\{ S(Z) - \mathbb{E}\left[S(Z)\right] \ge \varepsilon \right\} \le \mathrm{e}^{-2\varepsilon^2\eta_S} \le \mathrm{e}^{-2\varepsilon^2\eta_f}\,.\] With \(\varepsilon/\sqrt{\eta_f}\) in place of \(\varepsilon\), the bound follows. ◻
Remark 35 ().
Consider the pointwise coefficient \[\frac{1}{\gamma_f} := \sup_{\theta\in \Theta}\sum_{i=1}^n c_i(f, \theta)^2, \qquad\text{where}\qquad c_i(f, \theta) := \sup_{(z^\ast_i,z)\in \mathcal{Z}_i\times \mathcal{Z}} \left|
f_\theta(z_i^\ast,z_{-i})-f_\theta(z) \right| \,.\] Because the supremum over \(\Theta\) is inside the sum rather than on each coordinate, \(\gamma_f \ge \eta_f\), and the pointwise
McDiarmid bounds, e.g. 1 of the introduction, hold with \(\gamma_f\) in place of \(\eta_f\); thus pointwise tighter. But \(\gamma_f\) does not satisfy the lattice inequalities of Lemma 34, and Corollary 6 fails for \(\gamma_f\). \(\diamond\)
Proof of Lemma 23. Positive deviation. Let \(S := \sup_{t\in[0,1]} (f(t)-t)^+ \ge \varepsilon\). Choose a sequence \((t_n) \subset [0,1]\) such that \(f(t_n)-t_n \to S\). By compactness, passing to a subsequence if necessary, assume \(t_n \to t^\ast \in [0,1]\). Because \(f \le 1\), \(t_n \le 1 - (f(t_n) - t_n)\). Taking limits gives \(t^\ast \le 1 - S \le 1 - \varepsilon\). Thus, \(t^\ast + \alpha \in [0,1]\) for all \(\alpha \in [0, \varepsilon]\).
Fix \(\alpha \in [0, \varepsilon]\). For sufficiently large \(n\), \(t_n \le t^\ast + \alpha\). Monotonicity of \(f\) implies \(f(t^\ast + \alpha) \ge f(t_n)\). Thus, \[\begin{align} f(t^\ast + \alpha) - (t^\ast + \alpha) & \ge f(t_n) - t^\ast - \alpha \\ & = (f(t_n) - t_n) + t_n - t^\ast - \alpha \,. \end{align}\] Letting \(n \to \infty\) gives us \(f(t^\ast + \alpha) - (t^\ast + \alpha) \ge S - \alpha \ge \varepsilon - \alpha \ge 0\).
Negative deviation. Let \(S := \sup_{t\in[0,1]} (t-f(t))^+ \ge \varepsilon\). Choose \((t_n) \subset [0,1]\) such that \(t_n - f(t_n) \to S\), and assume \(t_n \to t^\ast \in [0,1]\). Because \(f \ge 0\), \(t_n \ge t_n - f(t_n)\), giving \(t^\ast \ge S \ge \varepsilon\). Thus, \(t^\ast - \alpha \in [0,1]\) for all \(\alpha \in [0, \varepsilon]\).
Fix \(\alpha \in [0, \varepsilon]\). For large \(n\), \(t_n \ge t^\ast - \alpha\). Monotonicity of \(f\) gives \(f(t^\ast - \alpha) \le f(t_n)\), and therefore \[\begin{align} (t^\ast - \alpha) - f(t^\ast - \alpha) & \ge t^\ast - \alpha - f(t_n) \\ & = (t_n - f(t_n)) + t^\ast - t_n - \alpha \,. \end{align}\] Letting \(n \to \infty\) gives us \((t^\ast - \alpha) - f(t^\ast - \alpha) \ge S - \alpha \ge \varepsilon - \alpha \ge 0\). ◻
Proof of Lemma 27. Write \(c:=2\ln2\), and set \[\Phi(y):=\ln\left( 4\sqrt{2\pi}\,y+4\sqrt{y}\,\mathrm{e}^{-y/2} +\mathrm{e}^{-y} \right), \qquad \Lambda(y)=\frac{\Phi(y)}{y} \,.\] Recall that \(y_0\) solves \[\ln(\mathrm{e}+y)=\sqrt{y} \,.\] Numerically, \(y_0=3.102884\ldots\), so \(y_0>3\).
Reduction by squaring. Multiplying the claimed bound by \(\sqrt{c/y}>0\) and isolating the square root, it is equivalent to \[\sqrt{\left(1+\frac{\Phi(y)}{y}\right)^2+\frac{c}{y}} \le 1+\frac{2\sqrt{c}\,\ln(\mathrm{e}+y)}{y} \,.\] Both sides are positive. Squaring and multiplying by \(y^2\), the equivalent inequality is \[2y\Phi(y)+\Phi(y)^2+cy \le 4\sqrt{c}\,y\ln(\mathrm{e}+y) +4c\,\ln(\mathrm{e}+y)^2 \,.\]
Majorising \(\Phi\). The function \[y\mapsto 4\sqrt{y}\,\mathrm{e}^{-y/2}+\mathrm{e}^{-y}\] is decreasing for \(y\ge1\). Hence, for \(y\ge3\), \[4\sqrt{y}\,\mathrm{e}^{-y/2}+\mathrm{e}^{-y} \le 4\sqrt{3}\,\mathrm{e}^{-3/2}+\mathrm{e}^{-3} \le (11-4\sqrt{2\pi})\,3 \,.\] Because \(y\ge y_0>3\), we have \[4\sqrt{2\pi}\,y +4\sqrt{y}\,\mathrm{e}^{-y/2} +\mathrm{e}^{-y} \le 11y \,.\] Therefore \[0<\Phi(y)\le \ln(11y), \qquad y\ge y_0 \,.\] Because the map \(\xi\mapsto2y\xi+\xi^2\) is increasing on \(\xi\ge0\), it is enough to prove \[2y\ln(11y)+\ln(11y)^2+cy \le 4\sqrt c\,y\ln(\mathrm{e}+y) +4c\,\ln(\mathrm{e}+y)^2 \,.\] Define \[\Psi(y) := 4\sqrt c\,y\ln(\mathrm{e}+y) +4c\,\ln(\mathrm{e}+y)^2 -2y\ln(11y)-\ln(11y)^2-cy \,.\] We prove that \(\Psi(y)\ge0\) for \(y\ge y_0\).
Convexity of \(\Psi\). Differentiating, \[\Psi'(y) = 4\sqrt c\left(\ln(\mathrm{e}+y)+\frac{y}{\mathrm{e}+y}\right) +\frac{8c\ln(\mathrm{e}+y)}{\mathrm{e}+y} -2\ln(11y)-2-\frac{2\ln(11y)}{y}-c\] and \[\Psi''(y) = \frac{4\sqrt{c}}{\mathrm{e}+y} +\frac{4\sqrt c\,\mathrm{e}}{(\mathrm{e}+y)^2} +\frac{8c(1-\ln(\mathrm{e}+y))}{(\mathrm{e}+y)^2} -\frac{2}{y} +\frac{2(\ln(11y)-1)}{y^2} \,.\] We show that \(\Psi''(y)\ge0\) for \(y\ge y_0\).
The only negative terms in \(\Psi''\) are \(-2/y\) and the negative part of the term containing \(1-\ln(\mathrm{e}+y)\). We use \[\frac{8c(\ln(\mathrm{e}+y)-1)}{(\mathrm{e}+y)^2} \le \frac{4c}{\mathrm{e}^3} \,,\] because \[\max_{u>0}\frac{\ln(u)-1}{u^2}=\frac{1}{2\mathrm{e}^3} \,.\]
The tail (\(y\ge6\)). For \(y\ge 6\), discard the nonnegative terms \[\frac{4\sqrt c\,\mathrm{e}}{(\mathrm{e}+y)^2}, \qquad \frac{2(\ln(11y)-1)}{y^2}\] from \(\Psi''\). After multiplying by \((\mathrm{e}+y)^2\), it is enough to prove \[\left( \frac{4\sqrt c}{\mathrm{e}+y}-\frac{2}{y} \right)(\mathrm{e}+y)^2 -8c(\ln(\mathrm{e}+y)-1) \ge0 \,.\] The first term equals \[(4\sqrt c-2)y +(4\sqrt c-4)\mathrm{e} -\frac{2\mathrm{e}^2}{y} \,.\] For \(y\ge6\), \[\frac{2\mathrm{e}^2}{y}\le\frac{\mathrm{e}^2}{3} \,.\] Thus, it is enough to prove \(\Xi(y)\ge0\), where \[\Xi(y) := (4\sqrt c-2)y +(4\sqrt c-4)\mathrm{e} -\frac{\mathrm{e}^2}{3} -8c(\ln(\mathrm{e}+y)-1) \,.\] But \[\Xi'(y) = (4\sqrt c-2)-\frac{8c}{\mathrm{e}+y}>0, \qquad y\ge6 \,,\] and \[\Xi(6)>2.79 \,.\] Hence \(\Psi''(y)\ge0\) for \(y\ge6\).
The bounded range (\(y_0\le y\le6\)). Using the global bound above in the formula for \(\Psi''\), we have \[\Psi''(y) \ge \frac{4\sqrt c}{\mathrm{e}+y} +\frac{4\sqrt c\,\mathrm{e}}{(\mathrm{e}+y)^2} +\frac{2(\ln(11y)-1)}{y^2} -\frac{2}{y}-\frac{4c}{\mathrm{e}^3} \,.\] The function \[y\mapsto\frac{\ln(11y)-1}{y^2}\] is decreasing on \([y_0,\infty)\), because \[\frac{\mathrm d}{\mathrm dy} \left(\frac{\ln(11y)-1}{y^2}\right) = \frac{3-2\ln(11y)}{y^3}<0 \,.\] Hence, on \([y_0,4]\), \[\Psi''(y) \ge \frac{4\sqrt c}{4+\mathrm{e}} +\frac{4\sqrt c\,\mathrm{e}}{(4+\mathrm{e})^2} +\frac{2(\ln44-1)}{16} -\frac{2}{y_0} -\frac{4c}{\mathrm{e}^3} >0 \,.\] Similarly, on \([4,6]\), \[\Psi''(y) \ge \frac{4\sqrt c}{6+\mathrm{e}} +\frac{4\sqrt c\,\mathrm{e}}{(6+\mathrm{e})^2} +\frac{2(\ln44-1)}{36} -\frac{1}{2} -\frac{4c}{\mathrm{e}^3} >0 \,.\] Therefore \(\Psi''(y)\ge0\) for all \(y\ge y_0\).
Because \[\Psi'(y_0)=1.440209\ldots>0\] and \[\Psi(y_0)=4.275887\ldots>0 \,,\] the function \(\Psi\) is increasing and positive on \([y_0,\infty)\). This proves the majorised inequality, hence the squared inequality, and therefore the scalar envelope bound. ◻