Abstract
We propose a structural default model for portfolio-wide valuation adjustments (xVAs) and represent it as a system of coupled backward stochastic differential equations. The framework is divided into four layers, each capturing a key component: (i) clean values, (ii) initial margin and Collateral Valuation Adjustment (ColVA), (iii) Credit/Debit Valuation Adjustments (CVA/DVA) together with Margin Valuation Adjustment (MVA), and (iv) Funding Valuation Adjustment (FVA). Since these layers depend on one another through collateral and default effects, a naive Monte Carlo approach would require deeply nested simulations, making the problem computationally intractable. To address this challenge, we use an iterative deep BSDE approach, handling each layer sequentially so that earlier outputs serve as inputs to the subsequent layers. Initial margin is computed via deep quantile regression to reflect margin requirements over the Margin Period of Risk. We also adopt a change-of-measure method that highlights rare but significant defaults of the bank or counterparty, ensuring that these events are accurately captured in the training process. We further extend Han and Long’s (2020) a posteriori error analysis to BSDEs on bounded domains. Due to the random exit from the domain, we obtain an order of convergence of O(h^{1/4-epsilon}) rather than the usual O(h^{1/2}). Numerical experiments illustrate that this method drastically reduces computational demands and successfully scales to high-dimensional, non-symmetric portfolios. The results confirm its effectiveness and accuracy, offering a practical alternative to nested Monte Carlo simulations in multi-counterparty xVA analyses.
Summary
The paper formulates portfolio-wide xVA computation as a stochastic control problem via a system of coupled BSDEs within a structural credit model framework. The bank and counterparty asset values follow diffusion processes, and default is triggered when these processes hit pre-specified barriers, making default times stopping times of the filtration. The adjusted portfolio value is defined implicitly as a conditional expectation of five types of discounted cashflows (contractual, collateral, initial margin, default, and funding), yielding a fixed-point equation that is well-posed by Banach’s fixed-point theorem. The key insight is to decompose this single high-dimensional BSDE into a layered system: Layer 1 solves clean values, Layer 2 computes initial margin via deep quantile regression, Layer 3 handles ColVA/CVA/DVA/MVA, and Layer 4 solves FVA. Each layer’s output feeds forward as input to subsequent layers.
The numerical solution employs a multi-layer deep BSDE solver that extends the original deep BSDE method of Han, Jentzen, and E (2018) in three directions: simultaneous solution of multiple coupled BSDEs sharing a single neural network function approximator, a Girsanov-based change of measure that tilts the forward SDE drift to increase default probabilities during training without altering the PDE solution, and incorporation of McKean-Vlasov BSDE features arising from the risk measure in the driver. The measure-change technique is particularly important because, under the original measure, default events are so rare that neural networks cannot learn the correct control behaviour near default boundaries. By simulating under a tilted measure and absorbing the reweighting into the BSDE generator, the algorithm avoids the large likelihood-ratio variance issues of traditional importance sampling.
On the theoretical side, the authors extend the a posteriori error analysis of Han and Long (2020) to BSDEs on bounded domains with random stopping times, obtaining a reduced convergence rate of O(h^{1/4-epsilon}) instead of the standard O(h^{1/2}). This is the first such result for deep BSDE methods under domain restrictions. Numerical experiments on a portfolio of 33 European basket call options on 5 underlying assets (a 93-dimensional problem in total) demonstrate that the method scales effectively and produces accurate results compared to both closed-form solutions and nested Monte Carlo benchmarks.
Key Contributions
- Formulation of the full xVA problem (CVA, DVA, FVA, ColVA, MVA) as a four-layer system of coupled BSDEs within a structural credit model, enabling Brownian-driver-only numerics
- Extension of the deep BSDE method to solve multiple BSDEs simultaneously with a shared neural network, reducing computational cost
- A Girsanov-based measure change that tilts forward SDE drift to increase default frequency during training, incorporated directly into the BSDE generator rather than via importance sampling reweighting
- Computation of initial margin via deep quantile regression using the check-function loss, producing VaR-based margin estimates at each time step
- First a posteriori error analysis for deep BSDE methods on bounded domains with stopping times, yielding O(h^{1/4-epsilon}) convergence
- Numerical demonstration on a 93-dimensional, non-symmetric portfolio of 33 basket options with 5 risk factors plus 2 defaultable entities
Methodology
The methodology combines several components arranged in a hierarchical pipeline:
- Forward SDE: A (d+2)-dimensional diffusion process X merging non-defaultable risk factors with bank and counterparty asset processes, with correlated Brownian drivers capturing wrong-way risk
- BSDE decomposition: The adjusted portfolio value BSDE is decomposed into individual xVA BSDEs (one per adjustment type) using the discounting cashflow approach of Brigo et al., enabling separate computation of each adjustment
- Layer 1: Clean value BSDEs solved simultaneously for all P derivatives via the deep BSDE method, with a single shared network Z: [0,T] x R^d → R^{sum d_j}
- Layer 2: Deep quantile regression for initial margin, training separate networks at each time step to estimate conditional VaR of portfolio value changes over the margin period of risk
- Layer 3: ColVA, CVA, DVA, MVA BSDEs solved with measure-tilted forward processes; the drift shift q is chosen to increase the default probability of the relevant party
- Layer 4: FVA BSDE, depending on all preceding layers
- Temporal discretization: Euler-Maruyama for the forward SDE, forward-propagated discretization for the BSDEs, with the unknown initial condition and control process Z approximated by neural networks
- Error analysis: Extension of Han-Long a posteriori bounds using Bouchard-Menozzi results for BSDEs on domains, yielding recursive error propagation across layers
Key Findings
- The multi-layer approach successfully scales to a 93-dimensional problem (33 basket options on 5 underlyings + 2 default processes) with 201 time steps
- The measure-change technique is essential: without it, the neural network learns Z approximately equal to zero near default boundaries, producing poor path-wise accuracy even when terminal condition errors appear comparable
- Clean values (Layer 1) match closed-form solutions with high accuracy across mean, 1st, and 99th percentiles
- Initial margin (Layer 2) approximations closely track nested Monte Carlo VaR benchmarks on representative paths
- ColVA, MVA, and CVA (Layer 3) approximations match nested Monte Carlo reference solutions on representative trajectories and show tight empirical distributions of terminal condition errors
- Under a stress test reducing counterparty default probability to approximately 0.6%, the measure-change technique maintains accuracy while the no-measure-change variant fails to capture fluctuations near default
- The O(h^{1/4-epsilon}) convergence rate for stopping-time BSDEs is slower than the standard O(h^{1/2}) rate, reflecting the inherent difficulty of accurately handling boundary exits
Important References
- Han, Jentzen, and E (2018) - “Solving high-dimensional partial differential equations using deep learning” (1940 citations). Foundational deep BSDE solver that this paper extends.
- Gnoatto, Reisinger, and Picarelli (2023) - “Deep xVA solver - A neural network based counterparty credit risk management framework”. The direct predecessor to this paper, which introduced the deep xVA solver in a reduced-form credit model setting.
- Han and Long (2020) - “Convergence of the deep BSDE method for coupled FBSDEs” (181 citations). Provides the a posteriori error analysis framework that this paper extends to bounded domains with stopping times.
Atomic Notes
- multi-layer deep BSDE solver
- change of measure for rare defaults in deep BSDE
- structural credit model
- deep quantile regression
- wrong-way risk in xVA
- xVA BSDE decomposition