cs.LGNov 17, 2025

On the Gradient Complexity of Private Optimization with Private Oracles

Authors: Michael MenartAleksandar Nikolov

Organizations: Department of Computer Science, University of Toronto, Vector Institute · Department of Computer Science, University of Toronto

Abstract

We study the running time, in terms of first order oracle queries, of differentially private empirical/population risk minimization of Lipschitz convex losses. We first consider the setting where the loss is non-smooth and the optimizer interacts with a private proxy oracle, which sends only private messages about a minibatch of gradients. In this setting, we show that expected running time Ω(min{dα2,dlog(1/α)})Ω(\min\{\frac{\sqrt{d}}{α^2}, \frac{d}{\log(1/α)}\}) is necessary to achieve αα excess risk on problems of dimension dd when d1/α2d \geq 1/α^2. Upper bounds via DP-SGD show these results are tight when d>Ω~(1/α4)d>\tildeΩ(1/α^4). We further show our lower bound can be strengthened to Ω(min{dmˉα2,dlog(1/α)})Ω(\min\{\frac{d}{\bar{m}α^2}, \frac{d}{\log(1/α)} \}) for algorithms which use minibatches of size at most mˉ<d\bar{m} < \sqrt{d}. We next consider smooth losses, where we relax the private oracle assumption and give lower bounds under only the condition that the optimizer is private. Here, we lower bound the expected number of first order oracle calls by Ω~(dα+min{1α2,n})\tildeΩ\big(\frac{\sqrt{d}}α + \min\{\frac{1}{α^2}, n\}\big), where nn is the size of the dataset. Modifications to existing algorithms show this bound is nearly tight. Compared to non-private lower bounds, our results show that differentially private optimizers pay a dimension dependent runtime penalty. Finally, as a natural extension of our proof technique, we show lower bounds in the non-smooth setting for optimizers interacting with information limited oracles. Specifically, if the proxy oracle transmits at most ΓΓ-bits of information about the gradients in the minibatch, then Ω(min{dα2Γ,dlog(1/α)})Ω\big(\min\{\frac{d}{α^2Γ}, \frac{d}{\log(1/α)}\}\big) oracle calls are needed. This result shows fundamental limitations of gradient quantization techniques in optimization.

Explore similar work

Jul 21, 2026cs.DS

Stronger Memory-Query Tradeoffs for Convex Optimization: The Limitations of Subquadratic Memory

We prove two lower bounds for the first order oracle complexity of minimizing a dd-dimensional 11-Lipschitz convex function over the unit ball with mm bits of memory. We first show that any such (possibly randomized) algorithm must make Ω~(d2m)\tildeΩ(\frac{d^2}{\sqrt{m}}) oracle queries. For deterministic optimization algorithms, we show that Ω~(min{d1.6,d8/3m2/3})\tildeΩ(\min\{d^{1.6},\frac{d^{8/3}}{m^{2/3}}\}) queries are required. For all memory regimes of interest, these improves upon the previous best known lower bounds of Ω~(max{d8/3m4/3,d4/3m1/6})\tildeΩ(\max\{\frac{d^{8/3}}{m^{4/3}},\frac{d^{4/3}}{m^{1/6}}\}) and Ω~(d5/3m1/3)\tildeΩ(\frac{d^{5/3}}{m^{1/3}}) for randomized and deterministic algorithms respectively. Notably, due to existing upper bounds, our lower bound for deterministic algorithms is the first to show a sharp oracle complexity phase transition around md2m\approx d^2, where a polylogarithmic change in memory leads to a poly(d)\mathsf{poly}(d) change in the number of required oracle calls. Further, when the suboptimality is polynomially small in dd, our lower bound randomized algorithms is the first to show that Ω~(d2)\tildeΩ(d^2) memory is necessary to nearly match the optimal query complexity among algorithms without memory constraints. Previously, such a result was only known for the regime where the suboptimality is quasipolynomially small in dd.
Michael Menart, Aleksandar Nikolov, Ohad Shamir
May 28, 2025cs.LG

Private Rate-Constrained Optimization with Applications to Fair Learning

Many problems in trustworthy ML can be expressed as constraints on prediction rates across subpopulations, including group fairness constraints (demographic parity, equalized odds, etc.). In this work, we study such constrained minimization problems under differential privacy (DP). Standard DP optimization techniques like DP-SGD rely on objectives that decompose over individual examples, enabling per-example gradient clipping and noise addition. Rate constraints, however, depend on aggregate statistics across groups, creating inter-sample dependencies that violate this decomposability. To address this, we develop RaCO-DP, a DP variant of Stochastic Gradient Descent-Ascent (SGDA) that solves the Lagrangian formulation of rate constraint problems. Through careful design, the extra privacy cost incurred by incorporating these constraints in our approach is limited to that of privately estimating a histogram over each mini-batch at every step. We prove the convergence of our algorithm through a novel analysis of SGDA that leverages the linear structure of the dual parameter. Empirical results show that our method Pareto-dominates existing private learning approaches under group fairness constraints and also achieves strong privacy-utility-fairness performance on neural networks.
Mohammad Yaghini, Tudor Cebere, Michael Menart +2
Jun 18, 2026cs.LG

On the Oracle Complexity of Interpolation-Based Gradient Descent

Recent work on first-order optimizers for empirical risk minimization (ERM) has suggested that smoothness of ERM loss functions in the training data, rather than in the optimization parameters, can be leveraged to improve the oracle complexity of gradient descent (GD) methods. In this paper, we propose an inexact gradient method, piecewise polynomial interpolation-based gradient descent (PPI-GD), which approximates the full gradient in each iteration by querying the first-order oracle at equidistant points in the data domain to construct polynomial interpolants of the resulting gradient samples over appropriately sized patches of the data domain. We analyze the oracle complexity of PPI-GD for strongly convex and non-convex loss functions when the data space dimension is bounded by a polylogarithmic function of the number of training samples, and find it to outperform several GD variants in key regimes when the loss function is sufficiently smooth. Furthermore, our analysis extends several techniques from the error analysis of bicubic spline interpolants to the setting of dd-variate tensor product polynomial interpolants which may be of independent interest in interpolation analysis.
Dongmin Lee, William Lu, Anuran Makur