REVIEW 1 major objections 3 minor 13 references
Expected improvement attains minimax-optimal simple-regret rates in Matérn and squared-exponential RKHSs.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · deepseek-v4-flash
2026-08-03 10:52 UTC pith:2UVUNY2O
load-bearing objection Fixed-prior EI finally gets the Matérn minimax rate and exponential SE rates; the main thing to verify is the one-dimensional SE upper bound that leans on Yarotsky. the 1 major comments →
Simple-regret rates and minimax optimality of fixed-prior expected improvement in Mat\'ern and squared-exponential RKHSs
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central claim is that fixed-prior expected improvement, used without randomized exploration, has essentially optimal simple-regret behavior on RKHS balls. For a Matérn kernel of smoothness ν on a compact subset of R^d, every weak-EI trajectory—one whose selected EI is at least a fixed fraction of the maximum—satisfies r_T = O(N^{-ν/d}) after N post-initial queries, uniformly over f in the RKHS ball of radius B. For the squared-exponential kernel, weak EI satisfies r_T = O(exp[-c min{N, N^{1/d} log(eN)}]), and exact EI maximization gives r_T = O(exp[-c N^{1/d} log(eN)]). If X has nonempty interior and B>0, exact EI is minimax-rate optimal among deterministic methods for Matérn kernels, an
What carries the argument
The load-bearing objects are the sequential separation radius and the Kolmogorov width. The innovation at a query point is the component of k(x,·) orthogonal to the span of previous evaluation representers; its norm equals the normalized posterior standard deviation. The separation radius δ_m(k,X) is the supremum over ordered m-point sequences of the smallest such innovation norm, and it bounds every ranked selected-point innovation over any finite horizon. Kolmogorov widths—the best uniform error in approximating the evaluation representers by an m-dimensional subspace—control the separation radii through Gram-determinant and eigenvalue estimates. The one-step regret inequality r_{n+1} ≤ q
Load-bearing premise
The load-bearing premise is that the best approximation of the kernel's evaluation functions by m basis functions decays at exactly the assumed rates (m^{-ν/d} for Matérn, e^{-c m^{1/d} log m} for squared-exponential) with no hidden logarithmic losses; the paper imports these scattered-data estimates rather than proving them.
What would settle it
Compute the sequential separation radius δ_m directly for the squared-exponential kernel in one dimension by greedy farthest-point sampling; the paper's Proposition E.5 predicts limsup_{m→∞} log δ_m/(m log m) ≤ -1/2, so a value strictly above -1/2 would disprove the refined exponential regret bound. For the Matérn claim, check whether the best m-point interpolation error on a cube really decays as O(m^{-ν/d}) uniformly over compact domains; any slower decay would invalidate the N^{-ν/d} upper bound.
If this is right
- On Matérn RKHS balls of radius B, exact EI (and any weak-EI policy) achieves worst-case simple regret O(N^{-ν/d}) after N post-initial evaluations, matching the deterministic minimax lower bound for every smoothness ν.
- On squared-exponential RKHS balls, exact EI achieves O(exp[-c N^{1/d} log(eN)]) simple regret, and no deterministic method can beat that exponent by more than a constant factor when the domain has interior.
- The weak-EI condition means the rates hold for any implementation that keeps a fixed positive fraction of the maximum EI, so approximate EI maximizers inherit the guarantees.
- The upper bounds are uniform over the RKHS ball and extend to a known continuous prior mean, so the conclusions apply to shifted function classes as well.
Where Pith is reading between the lines
- Editorial inference: Because the proof reduces regret to Kolmogorov-width decay, the same machinery should yield EI regret rates for other kernels with known width decay—for instance, Sobolev-type kernels of finite smoothness—though the paper only states the Matérn and squared-exponential cases.
- Editorial inference: The arbitrariness of the fraction η in the weak-EI condition suggests the rates are robust to coarse EI maximization; a practitioner who maximizes EI on a finite candidate grid should still attain them as long as the grid resolves a constant fraction of the true EI maximum—an empirically testable claim.
- Editorial inference: The minimax optimality is for deterministic, noiseless settings with a fixed prior. Noisy evaluations or data-driven kernel hyperparameters, which the paper explicitly leaves out, could change the rates; the feature-space mechanism alone does not predict by how much.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies the expected improvement (EI) policy for minimizing a deterministic function f in the RKHS of a continuous kernel on a compact set, under exact evaluations and a fixed zero-mean GP prior with covariance σ²k. The main claims are worst-case simple-regret bounds for weak-EI policies: O(N^{-ν/d}) for isotropic Matérn kernels and exp(-b' min{N, N^{1/d} log(eN)}) for the squared-exponential kernel; with exact EI maximization it claims the sharper SE bound exp(-c N^{1/d} log(eN)) and, on domains with nonempty interior, minimax-rate optimality for Matérn and optimality up to constants in the exponent for squared-exponential kernels. The proof methodology is a geometric reduction: posterior standard deviations are identified with feature-space innovation norms, their ranked values are bounded by sequential separation radii, width decay is transferred to separation-radius decay via Kolmogorov widths, and a one-step EI regret inequality is iterated through a finite-budget recursion.
Significance. If the concerns below are addressed, this is a substantial paper. It would resolve a real open question: fixed-prior EI, without randomized space-filling perturbations, attains the Matérn minimax rate for every smoothness ν>0, beyond the saturation in Bull (2011), and gives near-optimal squared-exponential rates. The paper is particularly strong in its general architecture: the innovation-norm viewpoint is elegant, the finite-budget bound of Theorem 4.4 is a genuinely useful tool, the constants are explicit and uniform over the RKHS ball, and no free parameters or fitting steps are introduced. The treatment of measurable exact-EI policies (Proposition F.1) and the use of external minimax lower bounds are careful. The Matérn upper bound and the d≥2 SE upper bound do not rely on the questionable step discussed below.
major comments (1)
- [§E.5, Proposition E.5; used in Theorem 6.6(2) and Theorem 7.4 (d=1)] The uniform super-exponential power-function bound (E.18) is the only source of the N log(eN) exponent in the d=1 SE upper bound; without it, the best available bound is the exponential (6.10), which is exponentially worse than the minimax lower bound exp(-C N log(eN)) in Theorem 7.4. The proof of (E.18) is a single citation: 'Yarotsky (2013, Theorem 2 and Equation (14)) ... gives (E.21)'. The cited theorem's hypotheses, direction, and domain of validity are not stated. Since (E.21) is then converted to (E.23) by replacing each factor |x̃-x̃_i|² with 4, the reader cannot verify that the external inequality is truly an upper bound valid for arbitrary distinct sites in [-1,1] rather than a lower bound or a bound requiring separated designs or additional smoothness conditions. This is load-bearing: if (E.21) does not hold in exactly the stated form, the d=1 part of the squared-exponential m
minor comments (3)
- [Eq. (E.22)] The expansion Fγ(K) = -K log K + O(K) is central to the exponent in (E.23). A short derivation of this expansion would help the reader check the constants and the uniformity in K.
- [Theorem 6.6(2)] The statement says '(6.11) holds for every N≥1'. As written this is only true after choosing A'' large enough to absorb a finite initial block; the clarity would be improved by saying 'for all N sufficiently large, with A'' adjusted for finitely many small N'.
- [Sections 6 and 7] The symbol N denotes the post-initial budget in Section 6 but the total number of evaluations in the minimax theorems of Section 7. The switch is not announced; please add a sentence in Section 7 to avoid confusion.
Circularity Check
No significant circularity: the EI regret bounds are genuine derivations from RKHS geometry and independent external approximation estimates; authors' prior self-citations are explicitly non-load-bearing.
full rationale
The derivation chain is self-contained in the sense required by the circularity check. The upper bounds follow from (i) the exact identity s_n(x)=dist(k(·,x),V_n) (Eq. 2.6), (ii) the one-step regret inequality (3.8) derived from the EI comparison Lemma 3.1, (iii) ranked-innovation bounds via Gram determinants and Kolmogorov widths (Lemmas 4.2, 5.1, 5.2), and (iv) external scattered-data width estimates (Schaback–Wendland 2002; Wendland 2004; Yarotsky 2013). No parameter is fitted to the target rates; the constants q_B, C_B depend only on B, sigma, eta. The authors' own Vazquez–Bect results are explicitly quarantined: 'The consistency and density results discussed in this subsection are recorded for completeness and insight. They are not used in the quantitative bounds that follow' (Section 2.5), and NEB/density is not needed for any of the rate theorems. The minimax lower bounds come from Bull (2011), Kühn (2011), and Xu et al. (2024), i.e., independent external work; they are not derived from EI, so the upper/lower comparison is not a self-comparison. The one genuinely fragile link, Proposition E.5's uniform super-exponential power-function bound imported from Yarotsky (2013), is a load-bearing external correctness risk, not a circular reduction: the paper does not define its conclusion into that citation, and the result is not an EI-specific fitted quantity. Thus central claims have independent content; only minor non-load-bearing self-citations are present.
Axiom & Free-Parameter Ledger
axioms (17)
- standard math RKHS theory: continuous PSD kernel k on compact X, reproducing property, restriction/extension theorems (Aronszajn 1950; Wendland 2004 Thms 10.46–10.47)
- standard math Gram–Schmidt identity: det K_x = product of squared successive innovation norms
- standard math Eckart–Young theorem for rank-i approximation in Hilbert–Schmidt norm
- standard math Greedy approximation product inequalities (Binev et al. 2011; DeVore et al. 2013; Wenzel et al. 2023)
- standard math Matérn power-function/fill-distance estimate (Schaback–Wendland 2002, Thm 3.2)
- standard math Squared-exponential power-function estimate on cubes (Wendland 2004, Thm 11.22 and Eq. (11.11))
- standard math Yarotsky 2013, Thm 2 conditional-variance bound for 1D squared-exponential kernel
- standard math Kühn 2011, Thm 3 metric entropy of squared-exponential RKHS ball
- standard math Xu et al. 2024, Thm 5.1 evaluation lower bound for noiseless optimization
- standard math Bull 2011, Thm 1 Matérn minimax lower bound
- domain assumption Objective lies in RKHS ball: f ∈ H_k, ||f||_{H_k} ≤ B
- domain assumption Exact noiseless observations of f at queried points
- domain assumption Fixed zero-mean GP prior with covariance σ²k, σ>0 fixed independently of data
- domain assumption Weak-EI condition: EI at selected point is at least η times the maximum EI, η∈(0,1]
- domain assumption Initial design fixed independently of f and of budget; kernel diagonal normalized to ≤1
- domain assumption For minimax comparisons: X has nonempty interior and B>0
- domain assumption Deterministic algorithms with Borel measurable final recommendation; randomized wrapper for lower bounds
read the original abstract
We study the expected improvement (EI) policy for minimizing a deterministic objective function $f$ on a nonempty compact set $\mathcal X \subset\mathbb R^d$. We assume that $f$ belongs to the RKHS $\mathcal H_k$ of a continuous positive-semidefinite kernel $k$ on $\mathcal X$. Function values are observed exactly, and EI is computed from a fixed zero-mean Gaussian-process model with covariance $\sigma^2k$. After an initial design, the policy queries a point whose EI is at least a fixed positive fraction of its maximum. We identify the normalized posterior standard deviation at a candidate point $x$ with the norm of the corresponding innovation in the canonical feature space, namely the component of $k(x,\cdot)$ orthogonal to the span of the preceding evaluation representers. Sequential separation radii bound the ranked innovation norms along arbitrary query sequences. We estimate these radii using Gram determinants and Kolmogorov widths for subspaces of different dimensions, then combine the estimates with a one-step regret inequality to obtain finite-budget bounds for simple regret. After $N$ post-initial queries, simple regret is $O(N^{-\nu/d})$ for isotropic Mat\'ern kernels of smoothness $\nu>0$. For the isotropic squared-exponential kernel, simple regret is $O(\exp[-c_1\min\{N, N^{1/d}\log(eN)\}])$ for some $c_1>0$. With exact EI maximization, it is $O(\exp[-c_2N^{1/d} \log(eN)])$ for some $c_2>0$. For every fixed $B\geq0$, these bounds are uniform over the RKHS ball of radius $B$. If $\mathcal X$ has nonempty interior and $B>0$, then, among deterministic methods whose final recommendation may be any point of $\mathcal X$, the exact EI policy is minimax-rate optimal over the RKHS ball of radius $B$ for Mat\'ern kernels and minimax-rate optimal up to constants in the exponent for squared-exponential kernels.
Reference graph
Works this paper leans on
-
[1]
N. Aronszajn. Theory of reproducing kernels.Transactions of the American Mathematical Society, 68(3):337–404, 1950.doi:10.1090/S0002-9947-1950-0051437-7. P. Binev, A. Cohen, W. Dahmen, R. DeVore, G. Petrova, and P. Wojtaszczyk. Convergence rates for greedy algorithms in reduced basis methods.SIAM Journal on Mathematical Analysis, 43(3):1457–1472, 2011.doi...
arXiv 1950
-
[13]
URL https://arxiv.org/abs/2501.09262v2. H. Wendland.Scattered Data Approximation, volume 17 ofCambridge Monographs on Applied and Computational Mathematics. Cambridge University Press, Cambridge,
-
[1936]
doi:10.2307/1968691. T. K¨ uhn. Covering numbers of Gaussian reproducing kernel Hilbert spaces.Journal of Complexity, 27(5):489–499, 2011.doi:10.1016/j.jco.2011.01.005. H. J. Kushner. A new method of locating the maximum point of an arbitrary multipeak curve in the presence of noise.Journal of Basic Engineering, 86(1):97–106,
-
[1962]
J. Moˇ ckus. On Bayesian methods for seeking the extremum. In G. I. Marchuk, editor, Optimization Techniques: IFIP Technical Conference, Novosibirsk, July 1–7, 1974, volume 27 ofLecture Notes in Computer Science, pages 400–404. Springer, Berlin, Heidelberg, 1975.doi:10.1007/3-540-07165-2_55. J. Moˇ ckus, V. Tiesis, and A.ˇZilinskas. The application of Bay...
-
[1964]
doi:10.1115/1.3653121. Y. Li and J. W. Siegel. Entropy-based convergence rates of greedy algorithms. Mathematical Models and Methods in Applied Sciences, 34(5):779–802,
-
[1978]
G. Santin, T. Wenzel, and B. Haasdonk. On the optimality of target-data-dependent kernel greedy interpolation in Sobolev reproducing kernel Hilbert spaces.SIAM Journal on Numerical Analysis, 62(5):2249–2275, 2024.doi:10.1137/23M1587956. R. Schaback and H. Wendland. Inverse and saturation theorems for radial basis function interpolation.Mathematics of Comp...
-
[1998]
doi:10.1023/A:1008306431147. A. N. Kolmogorov. ¨Uber die beste Ann¨ aherung von Funktionen einer gegebenen Funktionenklasse.Annals of Mathematics, 37(1):107–110,
-
[2002]
doi:10.1090/S0025-5718-01-01383-7. H. Tran-The, S. Gupta, S. Rana, and S. Venkatesh. Regret bounds for expected improvement algorithms in Gaussian process bandit optimization. InProceedings of the 25th International Conference on Artificial Intelligence and Statistics, volume 151 of Proceedings of Machine Learning Research, pages 8715–8737. PMLR,
-
[2004]
doi:10.1017/CBO9780511617539. T. Wenzel, G. Santin, and B. Haasdonk. Analysis of target data-dependent greedy kernel algorithms: convergence rates forf-,f·P-, andf /P-greedy.Constructive Approximation, 57(1):45–74, 2023.doi:10.1007/s00365-022-09592-3. W. Xu, Y. Jiang, E. T. Maddalena, and C. N. Jones. Lower bounds on the noiseless worst-case complexity of...
-
[2013]
doi:10.1007/s00365-013-9186-2. C. Eckart and G. Young. The approximation of one matrix by another of lower rank. Psychometrika, 1(3):211–218, 1936.doi:10.1007/BF02288367. R. Garnett.Bayesian Optimization. Cambridge University Press, Cambridge,
-
[2023]
ISBN 978-1-108-42578-0.doi:10.1017/9781108348973. S. Iwazaki. Gaussian process upper confidence bound achieves nearly-optimal regret in noise-free Gaussian process bandits. InAdvances in Neural Information Processing Systems, volume 38, pages 65863–65886. Curran Associates, Inc.,
-
[2024]
doi:10.1142/S0218202524500143. G. Matheron and P. Formery. Recherche d’optimum dans la reconnaissance et la mise en exploitation des gisements miniers. Note g´ eostatistique 41, Bureau de recherches g´ eologiques et mini` eres, Nov
-
[2025]
URL https://proceedings.neurips.cc/paper_files/paper/2025/file/ 5f1eee2509599faeeb3570a887016a64-Paper-Conference.pdf. 41 D. R. Jones, M. Schonlau, and W. J. Welch. Efficient global optimization of expensive black-box functions.Journal of Global Optimization, 13(4):455–492,
2025
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.