High-dimensional Laplace asymptotics up to the concentration threshold
Abstract
We study high-dimensional Laplace-type integrals I(λ):=(λ/2π)d/2∫ Rd g(x)e-λ f(x)dx in the regime where both d and λ are large. Existing rigorous Laplace-expansion results in growing dimension are largely confined to the "Gaussian-approximation" regime d2/λ0, which excludes many practically relevant settings that lie beyond this threshold but still satisfy the concentration condition d/λ0. We close this gap by deriving an explicit asymptotic expansion for I(λ) with quantitative remainder bounds that remain valid throughout this intermediate region, arbitrarily close to the concentration threshold. Fix L1 and assume that, in a neighborhood of the global minimizer of f, the operator norms of derivatives of f and g are bounded independently of d,λ up to orders 2L+2 and 2L, respectively. Assuming also some mild global growth conditions, we prove I(λ)=Σk=1L-1 bk(f,g)λ-k+O(dL+1/λL), dL+1/λL0, with coefficients satisfying bk(f,g)=O(dk+1). Moreover, the bk(f,g) coincide with the coefficients from the formal cumulant expansion of I(λ). We also study computation for concentrating densities π(x) e-λ f(x). For smooth observables g, our expansion yields closed-form, analytic approximations of EXπ[g(X)]. For sampling, we construct explicit polynomial transports xL such that πL:=(xL)\# N(0,λ-1Id) satisfies TV(π,πL) dL+1/λL for L=1,2,3,…, yielding an accurate procedure arbitrarily close to the concentration threshold d=o(λ).