---
title: Barycentric Quadratic Cost in Wasserstein Barycenters
url: https://www.emergentmind.com/topics/barycentric-quadratic-cost
type: topic
---

# Barycentric Quadratic Cost in Wasserstein Barycenters

The barycentric quadratic cost arises in the definition, computation, and analysis of Wasserstein-2 barycenters, which are average distributions between multiple probability measures under the quadratic cost optimal transport metric. This concept plays a foundational role in statistical inference, geometric measure theory, and computational optimal transport, as well as in modern data science applications that depend on meaningful “interpolations” or “averaging” of distributions. The barycentric quadratic cost provides both the theoretical underpinning and algorithmic structure for algorithms that compute barycenters, such as fixed-point iterations and quadratic programs. Its study involves convex geometry, functional analysis, Riemannian geometry on measure spaces, and numerical computation.

## 1. Definition and Variational Principles

Given probability measures $\mu_1,\dots,\mu_N$ on $\mathbb{R}^d$ with finite second moment and positive weights $\lambda_1,\dots,\lambda_N$ summing to 1, the quadratic barycentric problem is to find $\mu^* \in \mathcal{P}_2(\mathbb{R}^d)$ minimizing the weighted sum of squared 2-Wasserstein distances:
$$
\mu^* = \arg\min_{\mu \in \mathcal P_2(\mathbb{R}^d)} \sum_{i=1}^N \lambda_i W_2^2(\mu,\mu_i),
$$
where $W_2^2(\mu, \nu) = \min_{\pi \in \Pi(\mu, \nu)} \int_{\mathbb{R}^d \times \mathbb{R}^d} \|x - y\|^2 d\pi(x, y)$, with $\Pi(\mu, \nu)$ the set of couplings with marginals $\mu$ and $\nu$.

This variational principle directly encodes the barycentric quadratic cost as the functional whose minimizer is the Wasserstein-2 barycenter [2501.04016].

## 2. Equivalent Formulations and Geometric Structure

The barycentric quadratic cost admits several equivalent formulations central to existence, uniqueness, and computational analysis.

- **Two-marginal formulation (C2M):** The barycenter $\bar\nu$ minimizes $\sum_i\lambda_i W_2^2(\mu_i, \nu)$ over all $\nu$.
- **Multi-marginal formulation (MM):** Introduces a coupling $\gamma$ on $(\mathbb{R}^d)^N$,
  $$
  J(\gamma) = \int \sum_{i=1}^N \lambda_i |x_i - B(x_1, \dots, x_N)|^2 d\gamma,
  $$
  where $B(x_1,\dots,x_N) = \sum_i \lambda_i x_i$ is the weighted Euclidean barycenter. The pushforward $B_{\#}\gamma$ yields the barycenter [2402.13176].

Under standard regularity (e.g., one $\mu_i$ absolutely continuous), the barycenter is unique, the optimal multi-marginal plan is Monge (graph over the first marginal), and the quadratic cost structure guarantees strict convexity with respect to both measure and coupling [2402.13176].

## 3. Algorithmic Computation: Fixed-Point and Quadratic Programs

The barycentric quadratic cost structure enables efficient computational schemes:

- **Fixed-point iteration:** Following the Fréchet mean paradigm, the iterates
  $$
  \mu^{(k+1)} = \left( \sum_{i=1}^N \lambda_i T_i^{(k)} \right)_{\#} \mu^{(k)},
  $$
  where $T_i^{(k)}$ are optimal transport maps from $\mu^{(k)}$ to $\mu_i$, converge under standard hypotheses to the unique barycenter. The energy $V(\mu) = \sum_i \lambda_i W_2^2(\mu, \mu_i)$ decreases at every step, and weak limit points satisfy the fixed-point property. Uniqueness ensures full sequence convergence [2501.04016].

- **Quadratic program for barycentric coordinates:** In the barycentric coding model, determining the weights $\lambda^*$ so that a measure $\mu_0$ is a barycenter for measures $\mu_1,\dots,\mu_p$ reduces to minimizing
  $$
  \lambda^* = \arg\min_{\lambda \in \Delta^p} \lambda^\top A \lambda,
  $$
  where $A_{ij} = \int \langle T_i(x) - x, T_j(x) - x \rangle d\mu_0(x)$, and $T_i$ are optimal maps from $\mu_0$ to $\mu_i$. Under appropriate conditions, the unique minimizer yields the true barycentric coordinates [2201.12195].

- For Gaussian inputs, the barycentric quadratic cost program reduces to a closed-form QP in terms of means and covariances [2201.12195].

## 4. Statistical, Geometric, and Stability Properties

The barycentric quadratic cost induces fundamental geometric structure:

- **Riemannian geometry:** The tangent space at $\mu_0$ in $\mathcal{P}_2(\mathbb{R}^d)$ is represented by optimal displacement fields $\Delta_i(x) = T_i(x) - x$; the Gram matrix $A$ encodes inner products $\langle \Delta_i, \Delta_j \rangle_{\mu_0}$. The cost function $\lambda^\top A \lambda$ measures the squared norm of the barycentric displacement [2201.12195].
- **Injectivity and stability:** The barycenter map $B$ is globally invertible (Lipschitz inverse) and the quadratic cost enforces a strong injectivity estimate: $|B(y) - B(\tilde y)| \geq \Lambda |y - \tilde y|$, with $\Lambda = \min_i \lambda_i > 0$, on the support of the optimal plan. This produces absolute continuity and stability of barycenters [2402.13176].

Statistical guarantees are available via sample-based estimators. The quadratic cost enables convergence rates for barycentric coordinate estimation that are dimension- and smoothness-dependent, and the quadratic program is robust under entropic regularization, regularity, and point cloud approximations [2201.12195].

## 5. Extensions: Regularization and General Costs

The barycentric quadratic cost paradigm serves as the baseline for more general transport cost models:

- **Entropic regularization:** Replaces the cost functional $T_{c_k}$ with its $\epsilon$-regularized counterpart:
  $$
  T_{c_k,\epsilon}(\mu, \nu_k) = \min_{\pi \in \Pi(\mu,\nu_k)} \int c_k\,d\pi + \epsilon\,\mathrm{KL}(\pi\mid \mu\otimes\nu_k),
  $$
  with convergence and decrease properties preserved due to strict convexity of KL [2501.04016].

- **Generic cost functions:** The fixed-point method, multi-coupling construction, and corresponding decrease–compactness arguments generalize to arbitrary continuous costs $c_k(x, y)$ under suitable “ground barycentre” hypotheses [2501.04016].

- **Entropic and energy-based algorithmic methods:** Modern approaches for approximating quadratic-cost barycenters in high-dimensional or structured domains utilize dual formulations, neural parameterizations (potential functions), stochastic optimization, and MCMC sampling, leveraging the smoothness and geometry provided by the quadratic cost [2310.01105].

## 6. Applications and Computational Aspects

The barycentric quadratic cost underlies both theoretical and practical algorithms for measure averaging, coding, and estimation:

- **Covariance estimation—including closed-form solutions for Gaussian measures**—image processing, and natural language processing are highlighted application domains [2201.12195].
- **Computation:** Efficient Sinkhorn-based and energy-guided stochastic gradient methods exploit the quadratic cost structure for scalable high-dimensional barycenter estimation [2310.01105].
- **Sparsity and uniqueness:** The quadratic cost enforces sparse optimizers (Monge-type plans), crucial for both interpretability and efficiency [2402.13176].

## 7. Theoretical Guarantees and Open Directions

The barycentric quadratic cost confers the following guarantees in the computation of Wasserstein barycenters:

- Unique barycenter and optimizer under absolute continuity and strict convexity.
- Energy descent and statistical consistency of estimators; sample complexity $O(n^{-1/2})$ under empirical optimization [2201.12195, 2310.01105].
- Universal approximation by neural-network parameterized potentials for the entropic barycenter problem [2310.01105].
- Extension to more general cost functions, with the quadratic case serving as the canonical example by virtue of its analytical tractability and rich geometric interpretation [2501.04016, 2402.13176].

The study of barycentric quadratic cost thus forms the backbone of both foundational and algorithmic developments in Wasserstein barycenter computation and its generalizations.

Source: https://www.emergentmind.com/topics/barycentric-quadratic-cost