Smooth Surrogates for CVaR Optimization

• The paper introduces smooth surrogates for CVaR, demonstrating that EVaR and DD-GPCE-Kriging offer scalable, differentiable, and computationally efficient alternatives for risk optimization.
• The methodology employs convex programming and multifidelity sampling to accurately estimate tail risk while reducing computational burdens in high-dimensional settings.
• Key results show significant speedups and enhanced portfolio performance, validating the surrogate approaches as robust substitutes for traditional CVaR formulations.

Smooth surrogates for Conditional Value-at-Risk (CVaR) are a family of mathematical constructs and computational methodologies designed to overcome limitations inherent in the canonical definition of CVaR, especially its nonsmoothness with respect to underlying stochastic and optimization variables. These surrogates facilitate scalable and differentiable CVaR approximation or replacement within high-dimensional, nonsmooth, or sample-based risk quantification and optimization tasks. Two principal approaches have been recently advanced: the entropic value-at-risk (EVaR), a variationally tight, coherent, strongly monotone upper bound for CVaR that is infinitely differentiable; and high-dimensional, Gaussian-process–augmented polynomial chaos surrogates such as DD-GPCE-Kriging, whose global smoothness and multifidelity sampling enable efficient, accurate tail risk estimation.

1. Mathematical Foundations of CVaR and Smooth Surrogates

The Conditional Value-at-Risk at probability level $\beta\in(0,1)$, for a random variable $Y=y(\mathbf{X})$, is defined as the tail expectation beyond the $\beta$-quantile: $$\mathrm{CVaR}_\beta[Y] = \min_{\eta\in\mathbb{R}} \left\{ \eta + \frac{1}{1-\beta}\,\mathbb{E}[(Y-\eta)_+] \right\}$$ or, when the cumulative distribution $F_Y$ is continuous at $\mathrm{VaR}_\beta$,

$$\mathrm{CVaR}_\beta[Y] = \frac{1}{1-\beta}\,\mathbb{E}[\,Y\,\mathbb{I}_{\{Y\geq \mathrm{VaR}_\beta[Y]\}}\,].$$

The mapping $(y,\eta)\mapsto (y-\eta)_+$ is non-differentiable in $y$ and $\eta$, which propagates nonsmoothness into sample-based or optimization-based estimators and creates challenges for gradient-based optimization and surrogate modeling.

Smooth surrogates address this by either tightly upper-bounding CVaR with a variationally defined, smooth, coherent risk measure (e.g., EVaR) or by constructing globally smooth functional approximations (e.g., DD-GPCE-Kriging) that can replace or facilitate differentiation through the risk mapping (Ahmadi-Javid et al., 2017, Lee et al., 2022).

2. Entropic Value-at-Risk (EVaR): A Strongly Smooth Convex Surrogate

The entropic value-at-risk (EVaR) is defined for a real-valued loss $Y=y(\mathbf{X})$0 with moment-generating function $Y=y(\mathbf{X})$1 and tail-probability parameter $Y=y(\mathbf{X})$2 as: $Y=y(\mathbf{X})$3 or, equivalently, for $Y=y(\mathbf{X})$4,

$Y=y(\mathbf{X})$5

This is the tightest exponential (Chernoff) upper bound on VaR and CVaR available by Markov's inequality. Importantly, EVaR is coherent, strongly monotone, strictly monotone for all continuous distributions, and is $Y=y(\mathbf{X})$6 in both its arguments and all underlying statistical parameters—unlike CVaR itself, which lacks these monotonicity and smoothness properties (Ahmadi-Javid et al., 2017).

As $Y=y(\mathbf{X})$7, the optimal $Y=y(\mathbf{X})$8 and an analytic expansion shows $Y=y(\mathbf{X})$9, so EVaR converges to CVaR for deep-tail regimes in continuous distributions.

3. Convexity, Smoothness, and Computational Structures

EVaR admits a differentiable convex program structure. For $\beta$0 samples $\beta$1 (portfolio returns), weights $\beta$2, and a linear portfolio $\beta$3, define

$\beta$4

as the empirical EVaR objective. Its derivatives are:

• Gradient w.r.t. $\beta$5: $\beta$6
• Gradient w.r.t. $\beta$7: $\beta$8
• Hessian blocks: Provided in block notation for $\beta$9 and are all continuous and finite for $$\mathrm{CVaR}_\beta[Y] = \min_{\eta\in\mathbb{R}} \left\{ \eta + \frac{1}{1-\beta}\,\mathbb{E}[(Y-\eta)_+] \right\}$$0.

This yields a strictly convex, twice-differentiable optimization problem for portfolio weights $$\mathrm{CVaR}_\beta[Y] = \min_{\eta\in\mathbb{R}} \left\{ \eta + \frac{1}{1-\beta}\,\mathbb{E}[(Y-\eta)_+] \right\}$$1 and entropic parameter $$\mathrm{CVaR}_\beta[Y] = \min_{\eta\in\mathbb{R}} \left\{ \eta + \frac{1}{1-\beta}\,\mathbb{E}[(Y-\eta)_+] \right\}$$2. The number of optimization variables is $$\mathrm{CVaR}_\beta[Y] = \min_{\eta\in\mathbb{R}} \left\{ \eta + \frac{1}{1-\beta}\,\mathbb{E}[(Y-\eta)_+] \right\}$$3 plus user-imposed convex constraints, independent of $$\mathrm{CVaR}_\beta[Y] = \min_{\eta\in\mathbb{R}} \left\{ \eta + \frac{1}{1-\beta}\,\mathbb{E}[(Y-\eta)_+] \right\}$$4 (the sample size) (Ahmadi-Javid et al., 2017).

By contrast, the canonical CVaR LP reformulation introduces $$\mathrm{CVaR}_\beta[Y] = \min_{\eta\in\mathbb{R}} \left\{ \eta + \frac{1}{1-\beta}\,\mathbb{E}[(Y-\eta)_+] \right\}$$5 auxiliary variables and constraints, making large-$$\mathrm{CVaR}_\beta[Y] = \min_{\eta\in\mathbb{R}} \left\{ \eta + \frac{1}{1-\beta}\,\mathbb{E}[(Y-\eta)_+] \right\}$$6 optimization intractable by general-purpose solvers and yielding considerable overhead or memory exhaustion.

4. DD-GPCE-Kriging Surrogates for CVaR Estimation in High Dimensions

An alternate surrogate paradigm approximates the loss function $$\mathrm{CVaR}_\beta[Y] = \min_{\eta\in\mathbb{R}} \left\{ \eta + \frac{1}{1-\beta}\,\mathbb{E}[(Y-\eta)_+] \right\}$$7 by a globally smooth, high-dimensional surrogate: $$\mathrm{CVaR}_\beta[Y] = \min_{\eta\in\mathbb{R}} \left\{ \eta + \frac{1}{1-\beta}\,\mathbb{E}[(Y-\eta)_+] \right\}$$8 where

• $$\mathrm{CVaR}_\beta[Y] = \min_{\eta\in\mathbb{R}} \left\{ \eta + \frac{1}{1-\beta}\,\mathbb{E}[(Y-\eta)_+] \right\}$$9 is a vector of multivariate orthonormal polynomials on the support of the inputs, truncated by small interaction degree $F_Y$0 and total order $F_Y$1;
• $F_Y$2 are fitted coefficients;
• $F_Y$3 is a stationary, mean-zero Gaussian process (Kriging) used to correct local discrepancies;
• $F_Y$4 is a smooth, positive-definite correlation kernel (e.g., Gaussian or exponential).

The full surrogate is rendered infinitely differentiable (globally smooth) wherever $F_Y$5 is, and is constructed by regression from (potentially expensive) function evaluations of $F_Y$6 (Lee et al., 2022).

5. Algorithms for Efficient CVaR Estimation Using Smooth Surrogates

5.1 Sample-based Convex Optimization via EVaR

The EVaR convex program may be solved by primal-dual interior-point methods, which only require operations in the low-dimensional $F_Y$7 variable space and aggregate all $F_Y$8 samples via a log-sum-exp term. The algorithm iterates Newton steps with $F_Y$9 system size independent of the sample count, supporting millions of samples efficiently. Empirical results demonstrate orders-of-magnitude speedup over CVaR-LP formulations, which become prohibitive for large $\mathrm{VaR}_\beta$0 (Ahmadi-Javid et al., 2017).

5.2 Surrogate-based Monte Carlo and Multifidelity Importance Sampling

DD-GPCE-Kriging surrogates enable two major estimation strategies:

• Surrogate MCS: Fast MC is performed with the cheap surrogate $\mathrm{VaR}_\beta$1 as the loss proxy, increasing tail-sample efficiency but introducing surrogate bias of order $\mathrm{VaR}_\beta$2.
• Multifidelity Importance Sampling (MFIS): The surrogate is used solely to design a risk-region–biased sampling density. High-fidelity evaluations are then made only in these tail areas, and weighted by likelihood ratios to guarantee unbiased CVaR estimation. This hybrid approach combines the computational speedup of the surrogate with the statistical fidelity of the true $\mathrm{VaR}_\beta$3 in the relevant region.

Empirically, MFIS using DD-GPCE-Kriging achieves up to 104$\mathrm{VaR}_\beta$4 CPU speedup for composite finite-element models of dimension 20–28 and correlated inputs, with CVaR errors under 1–2% (Lee et al., 2022).

6. Comparative Performance and Practical Implications

A direct comparison between EVaR-based and standard CVaR-based optimization shows that, for sufficiently large $\mathrm{VaR}_\beta$5, the EVaR program is significantly faster and more memory-efficient, yielding nearly identical or superior portfolios in terms of expected return and tail risk. In a 20-asset S&P 500 study, EVaR-optimized portfolios outperformed CVaR portfolios in mean return (by up to +40% at $\mathrm{VaR}_\beta$6) and improved high-confidence VaR by +20% at only a marginal increase in standard deviation (5–15%). EVaR's strong and strict monotonicity is credited for mitigating both deep-tail and moderate losses (Ahmadi-Javid et al., 2017).

For nonsmooth high-dimensional outputs, DD-GPCE-Kriging with MFIS acutely reduces estimator bias relative to naive surrogate MC, and is scalable to complex systems with moderate sample budgets.

A summary of salient properties:

Surrogate Class	Differentiability	Variable/constraint count	Scalability to $\mathrm{VaR}_\beta$7
EVaR	$\mathrm{VaR}_\beta$8	$\mathrm{VaR}_\beta$9 ($$\mathrm{CVaR}_\beta[Y] = \frac{1}{1-\beta}\,\mathbb{E}[\,Y\,\mathbb{I}_{\{Y\geq \mathrm{VaR}_\beta[Y]\}}\,].$$0 constraints)	$$\mathrm{CVaR}_\beta[Y] = \frac{1}{1-\beta}\,\mathbb{E}[\,Y\,\mathbb{I}_{\{Y\geq \mathrm{VaR}_\beta[Y]\}}\,].$$1 function eval., variable dimension $$\mathrm{CVaR}_\beta[Y] = \frac{1}{1-\beta}\,\mathbb{E}[\,Y\,\mathbb{I}_{\{Y\geq \mathrm{VaR}_\beta[Y]\}}\,].$$2
CVaR (LP)	Piecewise-linear	$$\mathrm{CVaR}_\beta[Y] = \frac{1}{1-\beta}\,\mathbb{E}[\,Y\,\mathbb{I}_{\{Y\geq \mathrm{VaR}_\beta[Y]\}}\,].$$3	Limited by $$\mathrm{CVaR}_\beta[Y] = \frac{1}{1-\beta}\,\mathbb{E}[\,Y\,\mathbb{I}_{\{Y\geq \mathrm{VaR}_\beta[Y]\}}\,].$$4
DD-GPCE-Kriging	$$\mathrm{CVaR}_\beta[Y] = \frac{1}{1-\beta}\,\mathbb{E}[\,Y\,\mathbb{I}_{\{Y\geq \mathrm{VaR}_\beta[Y]\}}\,].$$5 (kernel)	Surrogate only	Scalable to $$\mathrm{CVaR}_\beta[Y] = \frac{1}{1-\beta}\,\mathbb{E}[\,Y\,\mathbb{I}_{\{Y\geq \mathrm{VaR}_\beta[Y]\}}\,].$$6 via cheap MC or MFIS

7. Broader Context and Outlook

Smooth surrogates for CVaR underpin a systematic shift in risk-aware modeling from nonsmooth, memory-intensive, sample-based formulations to compact, scalable, and differentiable paradigms. The EVaR construction offers a variationally tight, strongly monotone, and coherent substitute for CVaR, especially suited to convex portfolio optimization and large-scale sample regimes. DD-GPCE-Kriging enables smooth probabilistic risk estimation in high-dimensional, possibly dependent, and nonsmooth scenarios through multifidelity surrogate methodology.

A plausible implication is that as system dimension, input dependence, and sample count increase, smooth surrogates will become increasingly dominant in both risk estimation and optimization workflows, enabling robust risk management across finance and engineering domains without the computational costs historically associated with CVaR approaches (Ahmadi-Javid et al., 2017, Lee et al., 2022).