Wrong Way Risk corrections to CVA in CIR reduced-form models

In this paper we provide an efficient methodology to compute the credit value adjustment of a European contingent claim subject to some default event concerning the issuer solvability, when the underlying and the default event are correlated. In particular, in a Black and Scholes market/CIR intensity-default model, we consider a second order expansion around the origin of a vulnerable call option with respect to a correlation parameter ρ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho$$\end{document}, which may be used to describe the wrong way risk of the contract, measuring the dependence between the underlying asset price and the option’s issuer default intensity. Numerical implementations of this approach are compared with the benchmark Monte Carlo simulations.


Introduction
As a consequence of the exponential growth of the over-the-counter (OTC) financial market, the valuation of contingent claims subject to possible counterparty default, named counterparty credit risk (CCR), has gained interest in the last decade, after the past global financial crisis, both from the practitioner and the academic perspective.Indeed, according to the Basel III Accord (2010), to 46 Page 2 of 28 compensate for those risks a financial institution must charge a premium, generally known as credit value adjustment (CVA), defined as the (positive) difference between the default-free contract value and its value under default risks, or equivalently as the risk-neutral expectation of the residual discounted net present contract value (or expected exposure) at time of default (see e.g.Brigo et al. 2013;Gregory 2012).When speaking of options, these are called defaultable (or vulnerable), and if the default-free option value can be evaluated in closed or semiclosed form, an efficient option evaluation is tantamount to efficiently determining the contract's CVA.
The two main approaches to CCR are the structural and the reduced-form ones which both use the time of default to represent the default event.Johnson and Stulz (1987) firstly used the structural approach to price a vulnerable option, and this was later extended by Klein (1996) to more general liability structures.In these papers, default could happen only at maturity.More recently, Ballotta et al. (2019) proposed a more complex structural approach, based on correlated Lèvy processes.
In the structural approach, the default time is typically a stopping time with respect to the observable market filtration F .On the contrary, in the reduced-form (or intensity) approach, it is not measurable with respect to F , to incorporate also risk factors exogenous to the market.It is assumed that its distribution conditioned to the market filtration is represented by a regular hazard process with stochastic intensity.As a general reference on the topic, we quote the extensive work by Bielecki and Rutkowski (2002) and Bielecki et al. (2011) just to mention some.This is the framework we are going to use in the present paper.
Although the representation of CVA in terms of a risk-neutral expectation enables to employ many derivatives valuation techniques (from Monte Carlo to PDE based and/or hybrid techniques, see e.g.Goudenège et al. 2020), the actual CVA computation can be challenging when dependence is allowed in the underlying stochastic model.This problem is particularly relevant where statistical correlation between the counterparty default event and the investor exposure is possible; that is when wrong way risk (WWR) might occur, which means that a decrease in the counterparty's credit quality produces a higher exposure for the derivative's holder.A classic example of WWR (see Zhu and Pykhtin 2007) comes up when a bank (the investor I) enters into a swap with an oil producer (the counterparty C), receiving the fixed and paying the floating, hence lower oil prices have the double effect of worsening the counterparty credit quality, even though increasing the swap value for the bank.Another example of WWR is as follows: investor I buys a put option from bank C (the counterparty), written on another bank B, whose investment portfolio may be highly and negatively affected by defaulting of C (for example, both banks belong to the same financial system).Hence the higher the investor's exposure, the greater the probability of failure.
Of course, the opposite effect, known as right way risk (RWR), can also occur when exposure is likely to be reduced if the counterparty defaults.Both effects are essential from a risk management point of view; neglecting such dependencies between counterparty default and investor exposure leads to overestimation or underestimation of the risk and mis-pricing.
Spurred by these considerations, a wide literature to include RWR/WWR in derivative evaluation developed.An extensive presentation of the correlation issues concerning Credit Risk can be found in Brigo et al. (2013), Glasserman and Yang (2016) discussed bounds on WWR using linear programming techniques.Resampling methods, proposed in Pykthin and Rosen (2011), applied again in Rosen and Saunders (2012) and Cherubini (2013), are based on the use of a copula function to impose a dependence between exposure and credit, leaving fixed the marginal distributions, a method that proved to be particularly appealing to practitioners, for its simplicity and computational efficiency.When using the reduced-form approach, the hazard (or the intensity) process is typically taken as correlated to the option's underlying, see e.g.Alòs et al. (2021), Antonelli et al. (2021), Baviera et al. (2016), Brigo and Vrins (2018), Brigo et al. (2018), Feng and Osterlee (2017) and Hull and White (2012).
Following those ideas, we consider a vulnerable European plain vanilla option (or equivalently its unilateral CVA) in a Black and Scholes market model and a default intensity represented by a Cox-Ingersoll-Ross (CIR) diffusion, with correlated driving processes.Even in this very classical model, the presence of WWR makes the problem much more complex from a computational point of view.We derive a handy approximation formula based on a second-order expansion of the price representation formula with respect to the correlation coefficient, as an alternative to the benchmark Monte Carlo approach.Monte-Carlo methods (MC) are most commonly used in this setting (see e.g.Hull and White 2012;Rosen and Saunders 2012).Our method (stemmed from the works Antonelli and Scarlatti 2009;Antonelli et al. 2021 and extending those results), being analytical, leads to a faster methodology than MC simulations, nevertheless exhibiting comparable accuracy.In particular, the correlation expansion presented in Antonelli et al. (2021) involves a representation of the coefficients as solutions of a chain of parabolic PDE problems, whose solutions are characterized by the Feynman Kaĉ formulas.The procedure implies evaluating multiple integrals whose dimension grows with the order of approximation.In the current work, we succeed in simplifying the coefficients representation, using intermediate analytic derivative price formulas, obtained by conditioning and change-ofmeasure techniques.
Recently, a similar methodology has been applied to deal with CVA in stochastic volatility market models (see Alòs et al. 2021Alòs et al. , 2022)).Along the years, other value adjustments have been introduced and an evaluation theory developed, an extension of the current method to that context is to be found in Antonelli et al. (2022).We remark that this methodology not only provides representation formulas with handy approximations (due to expansion coefficients written at correlation zero), but it is also easily extendable to include multiple risk factors in the model, as shown in the quoted papers.In general, we refer the interested reader to the volume (Glau et al. 2016) for an updated overview on the topic.
The paper is organized as follows.Section 2 introduces the pricing problem for a vulnerable option and the related CVA using the intensity-based approach to write a general representation formula.In Sect.3, we present the market/stochastic intensity model and we specialize the representation formula to this case.In Sect.4, we introduce a second-order approximation for the price of the vulnerable option and 46 Page 4 of 28 its CVA and we propose a method to compute the expansion coefficients.Finally we test the efficiency and accuracy of our approximation procedure, comparing the results with those obtained by implementing enhanced Monte-Carlo simulations.

CVA: the intensity-based approach
In this section, we briefly recall the main setting and definition of the intensity approach to vulnerable options.
We consider a finite time interval [0, T] and a complete probability space (Ω, F, Q) , endowed with a filtration = {F t } t∈[0,T] , augmented with the Q-null sets and made right continuous.We assume that all processes have a cádlág version.
The market is described by the money market account denoted by D −1 (t, s) = e r(s−t) , r being the constant risk-free rate, and by a process {X t } t∈[0,T] rep- resenting an asset log-price (whose dynamics will be specified later).We assume to be in absence of arbitrage and that the given probability Q is a risk neutral measure, already selected by some criterion.In this market, a defaultable European option paying Φ(X T ) ≥ 0 at maturity, is traded.We denote by the default time of the issuer of the claim and by {Z t } t∈[0,T] a bounded F t -adapted recovery process.We model by means of the canonical construction of the default time as in Bielecki et al. (2011) section 3.2.2,therefore we assume there exists a differentiable adapted -hazard process for some non negative -adapted intensity process { t } t∈[0,T] , and we define for some uniform random variable U independent of .To properly evaluate this type of derivative, we need to include the information generated by the default time.Let us define the process and let ℍ = {H t } t∈[0,T] be the filtration it generates.For t ∈ [0, T] , we set the pro- gressively enlarged filtration G t = F t ∨ H t , which is the smallest filtration making a stopping time.We uniquely extend the risk neutral probability to and we keep denoting it by Q.
Following the canonical construction, the so-called H-hypothesis (see e.g.Bielecki and Rutkowski 2002;Bielecki et al. 2011) is automatically satisfied (Bielecki et al. 2011-Remark 3.2.1)and every -martingale remains a -martingale under Q, which implies that D(t, s)e X s , for s ≥ t , remains a -martingale.
In , for any given time t ∈ [0, T] , the price of a defaultable call is given by no arbitrage pricing by We apply this result and we denote by c d (t, T) the -adapted projection of c d, (t, T).
Consequently we have an -adapted CVA projection given by coinciding with CVA (t, T) on { > t}.
By using the canonical construction, we obtain a more treatable expression of c d (t, T) and of CVA(t, T).Indeed, let us introduce the conditional distribution of the default time given F t .Then F t () < 1 a.s.for all t > 0 and we have By the consequences of the Key-Lemma (see again Antonelli et al. 2021), we may rewrite (1) on { > t} as recovering formulas (3.1) and (3.3) in Lando (1998).For the sake of exposition we shall assume the recovery process to be identically null, leading to the simpler expression and consequently (4) becomes 46 Page 6 of 28 From now on, we will be working in , keeping in mind that all equalities are meant on the set { > t}.
Remark 1 In case of non zero recovery, the process {Z t } t∈[0,T] in (7) should be specified.The two most common choices are either Z t = Rc(t, T) (risk-free close- out hypothesis) or Z t = Rc d (t, T) (replacement close-out hypothesis) for a constant 0 ≤ R ≤ 1 .When characterizing the price by means of a Backward Stochastic Dif- ferential Equation (BSDE), the first choice leads to the expression (see e.g.Antonelli et al. 2021) while the second, solving the equation, gives implying or respectively.Indeed, both cases are straightforward extensions of (8), so considering R = 0 is not restrictive.
Remark 2 We recall that introducing the survival process and the (deterministic) survival function and it is possible to rewrite (for details, see Brigo and Vrins 2018).If the intensity process is independent of the underlying, the expectation in the CVA factorizes, with [ t ] = 1 , so we simply get ( 9) which is the basic representation of the (unilateral) CVA on which the regulatory adjustment formula (see Basel III) is built (see Gregory 2010).As a matter of fact, when the interest rates follow a deterministic process, by taking t = 0 (hence S 0 = 1 ) and defining EE(u) = [c(u, T)|F 0 ] , we have (having chosen the risk-free close-out hypothesis).Banks with an internal model method must compute the CVA according to the following formula where 0 = t 0 < t 1 < ⋯ < t N = T is a given time grid, the s k 's are the correspond- ing counterparty credit spreads, and LGD MKT is the market loss-given-default of the counterparty, assumed equal to LGD (see Hull 2012), it is immediately seen that the expression in ( 12) is essentially the trapezoidal quadrature rule approximation of the Riemann-Stieltjes integral (11), and, as such, it converges there when max k (t k − t k−1 ) → 0 (see e.g.Dragomir 2011).According to the Basel Commitee requirements, WWR is taken into account by simply multiplying the previous formula by a factor  > 1.2.

A representation formula under correlation
In this section, we introduce a representation formula for the price of a vulnerable plain vanilla call option, and the related CVA, in the case of correlated risk factors, assuming without loss of generality (see Remark 1) R ≡ 0.
The setting here is the same as in Antonelli et al. ( 2021), but we are going to use a different approach to evaluate the expansion coefficients in a handier and more efficient way, which allows pushing the expansion up to the second order in a quite straightforward manner.
We write our market model, in flow notation, for x ∈ ℝ and y > 0 to be where (B 1 , B 2 ) is a two dimensional standard Brownian motion and the constants , , , ∈ (0, +∞) .We also assume the Feller condition 2 >  2 to be satisfied, so  s > 0, ∀s ∈ [0, T] , Q-almost surely.The pair defining the market is jointly Marko- vian, thus, from (8), the price of a defaultable call option with maturity T and strike price e can be rewritten as a deterministic function of the state variables at time t, on the event { > t} , as for X t,x t = x and t t,y = y .The intensity is a Markov process also individually and since it is an affine process (see e.g.Alfonsi 2015 for an overview on affine processes), bond pricing theory implies that the expected survival process, given by P(y, t, T) ∶= Q (e − ∫ T t t,y s ds ) can be written as the latter being a continuous positive function.
Correlation in the dynamical model ( 13) induces a dependence between the derivative contract's value and the counterparty default probability.Indeed, for a call option, a growth of the underlying asset price increases the contract value.When correlation is positive, the corresponding increase in the intensity t raises the probability of default.
Employing a CIR process meets the desirable requirement that the intensity be positive, but it makes quite impossible to compute c d (x, y, t, T; ) in closed form, when ≠ 0 .Therefore, we aim at developing a handy representation to construct an approximation formula of polynomial order by power series expansion in around 0, if enough regularity is achieved.
Without loss of generality, for ease of exposition, we set t = 0 , r = 0 and we simplify the notation of 0,y s , P(y, 0, T) and c d (x, y, 0, T; ) to s , P(0, T) and c d (x, y, T; ).
Conditioning internally with respect to , where N denotes the Gaussian distribution, so the inner expectation gives a conditioned Black & Scholes formula  Bjork 2009 for the method), we have We notice here that the martingale which defines the new measure is different from the one proposed in Brigo and Vrins (2018), allowing a factorization of each expectation in ( 22) different from the factorization in Brigo and Vrins (2018) (see Sect. 5 for more details).By Girsanov theorem, we know that the process (23) is a Q T -Brownian motion, and consequently we rewrite the representation formula for the price of the vulnerable option under the T-forward measure We remark that for = 0 , the above formulas simplify and setting we have where c BS (x, T) denotes the Call option Black and Scholes pricing function.

The second order approximation
In this section, we compute a second order approximation for (25) and later we compare its accuracy with a benchmark Monte Carlo estimates of the prices given by ( 14).The first order expansion was already considered in Antonelli et al. (2021), where the change-of-numéraire technique was not employed.As shown before, this allows to isolate a bond-like expression coming from the default intensity and consequently a simpler handling of each term.First, we remark that (25) shows that c d is a regular function in and we may approximate it by a Taylor polynomial around 0 of any order and estimate, at least numerically, the error.As we are going to see, the coefficients computations are a bit involved, hence we focus on a second-order approximation formula, which gives very satisfying results in terms of efficiency and accuracy.This represents an improvement with respect to the first order approximation obtained in Antonelli et al. (2021), as it catches better the price behavior with respect to the correlation parameter , as we will comment later in the numerical section.Also, this extension keeps the advantage of being rather simple to implement and computationally fast compared to any Monte Carlo approach.
Thus, for k ≥ 1 , we set (25) and, by recalling that c BS (x, T) = e x F T (0) − e G T (0) , we aim at computing the fol- lowing second order approximation of the vulnerable option price Correspondingly, we have the second order expansion for the CVA (see (3)): In (28), we immediately recognize that the first contribution c BS (x, T)(1 − P(0, T)) corresponds to the CVA under independence.The remaining terms are therefore the correction due to RWR/WWR.By differentiating and taking = 0 , after some lengthy calculations (see "Appen- dix 1"), we arrive at (by virtue of ( 26)) where To make the previous approximation implementable, it remains to compute m(T, y) and s 2 (T, y).
(27) c d,(2) (x, y, T; ) = P(0, T) c BS (x, T) + e x g 1,1 (x, y, T) − e g 2,1 (x, y, T) e x g 1,2 (x, y, T) − e g 2,2 (x, y, T) . (28) 46 Page 12 of 28 Approximation of m(T, y).In Grzelak and Oosterlee (2011), it is given the following efficient approximation for the standard CIR model: Since the dynamic of the process s under the new measure Q T is where h(s) = 1 + 2 b(s, T) , we consider the approximation (33) with modified parameters: We tested the quality of the approximation ( 35) by comparing it with the direct application of Δ-method and the Monte Carlo estimation, for several parameter choices.Results are reported in "Appendix 2".Approximation of s 2 (T, y) .By Itô's integration by parts ( is a continuous, square integrable, finite variation process), we have s , (under the Feller's condition the drift term is well-defined), again by Itô's integration by parts we get Thus, and by approximating the process s by its expectation, Q T [ s ] , the second term gives zero contribution and we obtain the approximation Finally, we approximate Q T [( T ) 2 ] ≈ m(T, y) 2 , and we conclude We tested the quality of the approximations ( 35) and ( 36) on a randomly chosen sets of CIR model parameters: results are reported in "Appendix 2".
Remark 3 A theoretical estimation of the error made by stopping the Taylor expansion at the second order can be done by taking the third derivative of ( 25) with respect to and trying to bound it uniformly in [0, 1).This is not possible since that derivative will contain terms of the type , for n, m = 0, 1, 2, … , which can be bounded uniformly only in an inter- val [0, 0 ] for some  0 < 1.
If such an interval is fixed, then exploiting that a CIR process has bounded moments of any order, that the process is positive and so is , the moment generating function of the Gaussian random variable B1 T and the derivatives of the Gaussian density, one may prove that Remark 4 Once the representation ( 25) is obtained, it is reasonable trying to determine a hedging strategy, which is not straightforward since we are in an incomplete market.As discussed in Antonelli et al. (2022), if the market is rich enough, a hedging strategy might be devised by investing on the underlying (delta hedging) and on a bond that has as its rate of return, by taking derivatives of the approximation formula with respect to those two assets.Nevertheless, this hedging will not be perfect, leaving out a component in the price represented by a martingale orthogonal to the space generated by the two assets (underlying and bond).This martingale will depend on the risk neutral probability determining the price, which is not unique and that needs being selected by some criterion.What criterion to use (minimal martingale measure, minimal variance, etc.) is a delicate issue and it deserves a careful analysis that goes beyond the aim of the present paper and that we hope to address in future work.

The Brigo-Vrins approach
In a recent paper, Brigo and Vrins (2018) consider a general methodology based on a change of measures to address the problem of computing CVA under WWR in a stochastic-intensity default setup.Their starting point is formula (10) applied 46 Page 14 of 28 to a general portfolio price process V t .In such a framework, the EPE (expected positive exposure) under WWR is the function (see Remark 2 for notation) and Girsanov theorem is used to factorize the EPE.Indeed, they define an equivalent martingale measure Q C F,t ∼ Q according to from which they prove that The so defined measure Q C F,t is called wrong-way measure and it is associated with the numéraire C F,t ⋅ = D(0, ⋅)M t ⋅ .Hence, the dynamics of V t under the measure Q C F,t has to be obtained to com- pute the previous expression.Assuming V t to be a diffusion under Q, the change of measure results in a drift adjustment (see Brigo and Vrins 2018 for the full details).
In Brigo et al. (2018), this result was applied to CVA for a call option under WWR in the market model described by ( 13).Since the risk free rate is assumed to be constant, we have [D(0, t) t ] = −e −rt .Moreover, the explicit expres- sion of the new drift is with the functions a(⋅, ⋅) and b(⋅, ⋅) given respectively by ( 16) and (17).To compute the expectations in (37), it is necessary to replace the process { t } {t∈[0,T]} with a deterministic proxy { (t)} {t∈[0,T]} in (38), so that (37) can be evaluated analytically leading to the approximation (see Brigo et al. 2018): where Two proxies are considered for (t) , [ t ] and C F,t [ t ] : while the first is analytically known, the second requires a further approximation step (see Brigo et al. 2018).CVA under WWR is then obtained through a numerical integration procedure by inserting ( 39) in ( 10), specialized for the CIR intensity process, that is by taking J(t) = P(y, 0, t) (see ( 15)).

Numerical results
In this section, we aim at assessing the performance of the second-order approximation obtained by inserting the approximations ( 35) and ( 36) in ( 28), that we finally denote with where and ĝi,j being the coefficients evaluated with m , ŝ2 .Furthermore, we intend to size the contribution to CVA due to the correlation .We recall that Wrong Way Risk corresponds to  > 0 : the asset value and the probability of default move in the same direction.
We choose the intensity process parameters (Table 1) as in Brigo et al. (2018), Brigo and Alfonsi (2005): in particular, Set a was exogeneously given and it corresponds to the scenario with the highest (1 year) default probability; Set b and Set c come from a calibration process on CDS market quotes.In the first set of experiments, the strike price was fixed to K = e = 100 and we three maturi- ties, T = 0.5 , T = 1 and T = 5 : the underlying volatility was = 10% and without loss of generality we also set the risk-free rate r = 0 and t = 0 .The log-asset was set to 4.6052.All the algorithms were implemented in MatLab © (R2019b), and ran on an Intel Core i7 2.40 GHZ with 8 GB RAM.
The benchmark Monte Carlo method was implemented by using the Euler discretization scheme with full truncation for the intensity process { t } {t∈[0,T]} (see Lord et al. 2010) and the exact simulation of the Brownian motion for the underlying {X t } {t∈[0,T]} .In order to reduce the variance of the estimator, a control variates tech- nique was considered, by using the default-free price as control: in the considered ( 40) h 1 (T) ∶= P(0, T) e x ĝ1,1 (x, y, T) − e  ĝ2,1 (x, y, T) h 2 (T) ∶= P(0, T) e x ĝ1,2 (x, y, T) − e  ĝ2,2 (x, y, T) , cases, this reduced the length of the confidence interval by at least one order of magnitude.In our numerical experiments we generated M = 10 6 sample paths with a time step discretization equal to 1.0e−03 for all the maturities.The length of the 95% confidence interval varies between 1.0e−03 and 1.0e−04.The implemented algorithm took about 52 s for the computation of each CVA value, ĈVA (MC) (T, ).The results of the numerical implementation are reported in Figs. 1, 2 and 3 and they show the quality of the proposed second order approximation.We also observed that the use of the different approximations of m(T, y), as reported in "Appendix 2", does not particularly significantly impact the quality of the estimator (40).A comparison with the first order approximation, i.e. and the Brigo-Vrins (BV) method are also displayed (see Sect. 5).In particular, the latter required to evaluate numerically the integral (10) and the integral defining the function Θ(t) for the computation of EPE(t), 0 < t ≤ T .Both integrations were implemented by using a trapezoidal routine: notice that the function Θ(t) must be computed for every value of the correlation parameter .The implemented procedure takes approximately 4.5e−02 s to calculate each CVA value, ĈVA (BV) (T, ).As it can be seen from Tables 2, 3 and 4 | , worsens for increasing maturity.The first and second-order approximations behave similarly for small , while for larger values, the second order is constantly better, especially in the most extended maturity scenario.We observe that the first-order approximation and the BV method produce in many instances similar relative error patterns, except for the parameter Set b, where the BV method seems to get worse for small T, as correlation decreases, although the error remains in the order of 1.0e−03.In contrast, the second-order approximation is systematically better at the cost of a minimal increase in computational complexity.As a matter of fact, we stress that it is sufficient to calculate the (approximated) coefficients ĝi,j only once to obtain an entire CVA curve as a function of , unlike the   Monte Carlo and the BV methods, which require a new run of simulations and computations for each choice of the correlation parameter.The computation of the first and second order coefficients h 1 and h 2 in our implementation took about 5.0e−04 and 4.0e−03 s, respectively: the overall evaluation of the our second order approximation formula took approximately 2.0e−03 s.The quality of ĈVA (2) is particularly appreciable in the case of Set b, where the convexity of the CVA figure seems most evident (see Fig. 2).
We further remark that the first-order approximation ĈVA (1) yields results comparable to that proposed in Antonelli et al. (2021).The coefficient of this approximation in the Taylor expansion is, in fact, the same but it has a more explicit representation in the present paper.
In the second set of experiments, we report the values of the (approximated) coefficients of the second-order expansion, i.e. the credit value adjustment under independence, CVA ind ∶= c BS (x, T) 1 − P(0, T) , and the terms h 1 and h 2 (see ( 40)), in Tables 5, 6 and 7 referring to Set a, Set b and Set c, respectively.
The coefficients h i , i = 1, 2 in formula ( 40) are found to be negative.This result is widely expected since it can be supposed that the Credit Value Adjustment (CVA) would increase when the parameter is positive.This is because, with a positive , the probability of the counterparty defaulting increases as the value of the contract increases.As a result, a higher credit adjustment is required to account for this risk.
Moreover, h 2 is found to be typically one or two times smaller than h 1 , except in Set b, where it is often comparable (or even larger) in magnitude to h 1 .This finding is consistent with the first onset of experiments (see Fig. 2) where it is observed that the dependence of the CVA on the parameter exhibits non-negligible convexity, which is reflected in the magnitude of the corresponding coefficient.Consequently, the coefficient introduces a non-negligible correction term, and it highlights the importance of considering the impact of parameter dependencies on CVA calculations.
Overall, these findings emphasize the significance of the parameter involved in CVA calculations to account for WWR/RWR, and they suggest the need for a comprehensive analysis of their impact on the resulting derivative values.

Conclusions
This paper addresses the issue of estimating CVA within the framework of intensity models, while accounting for wrong-way risk (WWR), which refers to the presence of a correlation between the derivative contract's value and the counterparty's default probability.We propose a technique based on approximating the CVA value as a function of the correlation parameter using first or second-order Taylor polynomials.The coefficients of this approximation are obtained using a change-of-numeraire approach.We demonstrate the effectiveness of our method in approximating the coefficient of the Taylor polynomials for a plain vanilla call option, where the underlying follows a GBM dynamic and the intensity of default is a correlated CIR process.Our approach proves to be both accurate and efficient, as shown by a series of numerical experiments compared with the results obtained by Monte Carlo estimations and the method in Brigo 46 Page 20 of 28

Table 5
The values of the coefficients of the second-order approximation for the CVA, for different time-to-maturity TtM

Table 7
The values of the coefficients of the second-order approximation for the CVA, for different time-to-maturity TtM  (2018).Additionally, our technique enucleates the WWR contribution of negative correlation providing a quantitative measure for the impact of such a risk.The approach in the present paper opens to further research in several directions.Firstly, it could be extended to consider other risk drivers, such as stochastic interest rates, which can significantly affect CVA estimates.Secondly, the method could be applied to other types of underlying/intensity models, including those with jumps or mean-reverting features.Lastly, it could be worth trying to apply it on other types of derivatives, such as forward starting options or path-dependent products, which are more complex and can pose greater challenges for CVA estimation.These extensions could help to further validate the proposed method and broaden its applicability in practical settings.
Besides approximation (45), in Grzelak and Oosterlee (2011) Grzelak and Oosterlee proposed to use, for the CIR model, the following function where the coefficients C 1 , C 2 and C 3 are set to fit as much as possible to the function Λ(t) (see Sect. 4).We tested the quality of the Δ-method approximation of Q T [ √ t ] (45) and its functional approximation (49) for a set of parameters , , and 0 : some results are displayed in Fig. 4, where they are compared with a Monte Carlo approximation, obtained from simulation of the process (46) through the Euler discretization scheme with full truncation, with time step 1.0e−03, and 1.0e+05 trajectories (Fig. 5).
, the absolute relative error with respect to the benchmark Monte Carlo estimation as a function of the correlation parameter ,

Fig. 1
Fig.1CVA approximations for the Set a: blue circles are MC estimates, magenta diamonds and red stars are the first and second order approximations, respectively, and green squares are the results of the BV method (color figure online)

Fig. 2
Fig. 2 CVA approximations for the Set b: blue circles are MC estimates, magenta diamonds and red stars are the first and second order approximations, respectively, and green squares are the results of the BV method (color figure online)

Fig. 4
Fig. 4 Approximation ofQ T [ √ t ].The dotted lines represent the Monte Carlo estimation, the squared dash-dotted and the star dash-dotted are the Δ and its functional approximation, respectively here 1 A denotes the indicator function of a set A), being the standard normal distribution function, where for i = 1, 2 , we have denoted the random variables Hence, from (19), we obtain where The above representation can be made a little more explicit by an appropriate change of Numéraire.Introducing the F -adapted Q-martingale

Table 1
Parameter sets for the CIR default intensity The 1-year survival probability for Set a and Set b is 96.9% , 98.7% and 98.5% , respectively

Table 3
Set b: the absolute relative error with respect to the Monte Carlo estimate