Next Article in Journal
Generalized Multiplicative Risk Apportionment
Next Article in Special Issue
Risk Factor Evolution for Counterparty Credit Risk under a Hidden Markov Model
Previous Article in Journal
A Renewal Shot Noise Process with Subexponential Shot Marks
Previous Article in Special Issue
The Determinants of Market-Implied Recovery Rates
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Default Ambiguity

by
Tolulope Fadina
1 and
Thorsten Schmidt
1,2,3,*
1
Department of Mathematical Stochastics, University of Freiburg, Ernst-Zermelo Str.1, 79104 Freiburg im Breisgau, Germany
2
Freiburg Research Institute of Advanced Studies (FRIAS), 79104 Freiburg im Breisgau, Germany
3
University of Strasbourg Institute for Advanced Study (USIAS), 67081 Strasbourg, France
*
Author to whom correspondence should be addressed.
Submission received: 29 March 2019 / Revised: 17 May 2019 / Accepted: 17 May 2019 / Published: 10 June 2019
(This article belongs to the Special Issue Advances in Credit Risk Modeling and Management)

Abstract

:
This paper discusses ambiguity in the context of single-name credit risk. We focus on uncertainty in the default intensity but also discuss uncertainty in the recovery in a fractional recovery of the market value. This approach is a first step towards integrating uncertainty in credit-risky term structure models and can profit from its simplicity. We derive drift conditions in a Heath–Jarrow–Morton forward rate setting in the case of ambiguous default intensity in combination with zero recovery, and in the case of ambiguous fractional recovery of the market value.

1. Introduction

Recently, an increasing amount of literature focuses on uncertainty as it relates to financial markets. The problem is that, the probability distribution of randomness in these markets is unknown. Typically, the unknown distribution is either estimated by statistical methods or calibrated to given market data by means of a model for the financial market. For example, in credit risk, the default probability is not observed, hence, have to be estimated from observable data. These methods introduce a large model risk.
Already, Knight (1921) pointed towards a formulation of risk which is able to treat such challenges in a systematic way. He was followed by Ellsberg (1961), who called random variables with known probability distribution certain, and those where the probability distribution is not known uncertain.
In this paper, we address these problems by constructing a model such that the parameters are characterised by uncertainty. Then, a single probability measure in a classical model is replaced by a family of probability measures, that is, a full class of models.
Following the modern literature in the area, we will call the feature that the probability distribution is not entirely fixed or, cannot be modelled by a single probability measure, ambiguity. This area has recently renewed the attention of researchers in mathematical finance to fundamental subjects such as arbitrage conditions, pricing mechanisms, and super-hedging. In equity markets, volatility uncertainty plays a crucial role, and has been extensively investigated, see for example, Avellaneda et al. (1995); Denis and Martini (2006); Lyons (1995); Vorbrink (2014). A major difficulty in this setting is that volatility uncertainty is characterized (at least in continuous time) by probabilities measures being mutually singular. Thus, the classical fundamental theorem of asset pricing fails to justify the no-arbitrage conditions, and new techniques are demanded, see Bayraktar and Zhang (2013); Bouchard and Nutz (2015); Burzoni et al. (2017), and Biagini et al. (2017).
In this paper, we introduce the concept of ambiguity to defaultable term structure models. The starting point for term structure models are typically bond prices of the form
P ( t , T ) = e t T f ( t , u ) d u , 0 t T ,
where ( f ( t , T ) ) 0 t T is the instantaneous forward rate and T is the maturity time. This follows the seminal approach proposed in Heath et al. (1992). The presence of credit risk1 in the model introduces an additional factor known as the default time. In this setting, bond prices are assumed to be absolutely continuous with respect to the maturity of the bond. This assumption is typically justified by the argument that, in practice, only a finite number of bonds are liquidly traded and the full term structure is obtained by interpolation, thus is smooth. There are two classical approaches to model market default risk: the structural approach Merton (1974) and the reduced-form approach (see for example, Artzner and Delbaen (1995); Duffie et al. (1996); Lando (1994) for some of the first works in this direction).
Structural models of credit risk describe the modelling of credit events specific to a particular corporate firm. Here, the underlying state is the value of a firm’s assets which is observable. Default time is define as the first time the firm’s asset value falls below a certain barrier level (for example, its liabilities). Hence, default is not a surprise. The approach links the default events to the firm’s economic fundamentals, consequently, default time is endogenously within the model. However, the assumption that the firm’s asset value is observable is often too strong in real applications, see Duffie and Lando (2001) and the survey article Frey and Schmidt (2011). Structural models in credit risk have been studied under many different viewpoints, see Black and Cox (1976); Frey and Schmidt (2009, 2012); Gehmlich and Schmidt (2018); Geske (1977); Kim et al. (1993); Leland and Toft (1996); Merton (1974).
In comparison to the structural approach, reduced form approaches take a less stringent viewpoint regarding the mechanisms leading to default and model default events via an additional random structure. This additional degree of freedom together with their high tractability led to a tremendous success of this model class in credit risk modelling. For more details on the reduced form approach, we refer to Bielecki and Rutkowski (2002), and the references therein. Structural models can be embedded into (generalized) reduced form models as pointed out in Bélanger et al. (2004) and Fontana and Schmidt (2018).
Reduced-form models typically postulate that default time is totally inaccessible and, consequently, bond prices are absolutely continuous with respect to the maturity. Under the assumption of zero recovery2, this implies that credit risky bond prices P ( t , T ) are given by
P ( t , T ) = 𝟙 { τ > t } e t T f ( t , u ) d u
with τ denoting the random default time. Since τ is totally inaccessible, it has an intensity λ . For example, if the intensity is constant, the default time is the first jump of a Poisson process with constant arrival rate λ . More generally, λ t may be viewed as the conditional rate of occurrence of default at time t, given information up to that time. In a situation where the owner of a defaultable claim recovers part of its initial investment upon default, the associated survival process 𝟙 { τ > t } in (2), is replaced by a semimartingale. The quantity of the investment recovered is the so-called recovery.
Under ambiguity, we suggest that there is some prior information at hand which gives a upper and lower bounds on the intensity. The implicit assumption that the probability distribution of default is known is quite restrictive. Thus, we analyse our problem in a multiple priors model which describe uncertainty about the “true probability distribution". By means of the Girsanov theorem, we construct the set of priors from the reference measure. The assumption is that all priors are equivalent.
In view of our framework, it is only important to acknowledge that a rating class provides an estimate of the one-year default probability in terms of a confidence interval. Also estimates for 3-, and 5-year default probabilities can be obtained from the rating migration matrix. Thus, leading to a certain amount of model risk.
The aim of this paper is to incorporate ambiguity into the context of single-name credit risk. We focus on ambiguity on the default intensity, and also discuss ambiguity on the recovery.
The main results are as follows: we obtain a necessary and sufficient condition for a reference probability measure to be a local martingale measure for credit risky bond markets under default ambiguity, thereby ensuring the absence of arbitrage in a sense to be precisely specified below. Furthermore, we consider the case where we have partial information on the amount that the owner of a defaulted claim receives upon default.
The next section of this paper introduces homogeneous ambiguity, and its example. Section 3 introduces the fundamental theorem of asset pricing (FTAP) under homogeneous ambiguity. In Section 4 and Section 6, we derive the no-arbitrage conditions for defaultable term structure models with zero-recovery, and fractional recovery of the market value, in our framework. We conclude in Section 7.

2. Intensity-Based Models

Intensity-based models are the most used model class for modeling credit risk (see (Bielecki and Rutkowski 2002, Chapter 8) for an overview of relevant literature). The default intensity, however, is difficult to estimate and therefore naturally carries a lot of uncertainty. This has led to the emergence of rating agencies which, since the early 20th century, estimate bond’s credit worthiness3.
Modeling of credit risk has up to now incorporated uncertainty in the default intensity in a systematic way. On the other side, a number of Bayesian approaches exist, utilizing filtering technologies (see, for example, Duffie and Lando (2001); Frey and Schmidt (2009), among many others).
Here, we introduce an alternative treatment of the lack of precise knowledge of the default intensity based on the concept of ambiguity following the seminal ideas from Frank Knight in Knight (1921).
Uncertainty in our setting will be captured through a family of probability measures P replacing the single probability measure P in classical approaches. Intuitively, each P represents a model and the family P collects models which we consider equally likely.
In this spirit, working with a single P , or with a set P = { P } which contains only one element, is in a one-to-one correspondence to assuming that the parameters of the underlying processes are exactly known. In financial markets, this is certainly not the case and ambiguity helps to incorporate this uncertainty into the used models.
We consider throughout a fixed finite time horizon T * > 0 . In light of our discussion above, let ( Ω , F ) be a measurable space and P be a set of probability measures on the measurable space ( Ω , F ) . In particular, there is no fixed and known measure P (except in the special case where P contains only one element which we treat en passant).
Intensity-based default models correspond to the case where the ambiguity is homogeneous, i.e., there is a measure P such that P P for all P P . Here, P P means that P and P are equivalent, that is, they have the same nullsets. The reference measure P only has the role of fixing events of measure zero for all probability measures under consideration. Intuitively, this means there is no ambiguity on these events of measure zero. In the following, E is the expectation with respect to the reference measure P .
Remark 1.
As a consequence of the equivalence of all probability measures in P , all equalities and inequalities will hold almost-surely with respect to any probability measure P P , or, respectively, to P .

Ambiguity in Intensity-Based Models

In this section, we introduce ambiguity in intensity-based models. Our goal is not the most general approach in this setting: we rather focus on simpler, but still practically highly relevant cases. For a more general treatment, we refer the reader to Biagini and Zhang (2017). The main mathematical tool we use here is enlargement of filtrations and we refer the reader to Aksamit and Jeanblanc (2017) for further details and a guide to the literature.
Assume that under P we have a d-dimensional Brownian motion W with canonical and augmented4 filtration F = ( F t ) 0 t T * and a standard exponential random variable τ , independent of F T * , that is, P ( t < τ | F t ) = exp ( t ) , 0 t T * . The Brownian motion W has the role of modelling market movements and general information, excluding default information. We therefore call F the market filtration in the following. The filtration G = ( G t ) 0 t T * includes default information and is obtained by a progressive enlargement of F with τ , i.e.,
G t = ϵ > 0 σ ( 𝟙 { t τ } , W s : 0 s t + ϵ ) , 0 t T * .
To finalize our setup, we assume that F = G T * .
Note that up to now, everything has been specified under the reference measure P and nothing was said about the concrete models we are interested in (except about the nullsets). These models will now be introduced using the Girsanov theorem, i.e., by changing from P to the measures we are interested in.
Consequently, the next step is to construct measures P λ with appropriate processes λ – under P λ , the default time τ will have the intensity λ . More precisely, assume that λ is some positive process which is predictable with respect to the market filtration, F . Define the density process Z λ by
Z t λ : = exp 0 t ( 1 λ s ) d s , t < τ λ τ exp 0 τ ( 1 λ s ) d s t τ .
Note that Z λ is a G -local martingale and corresponds to a Girsanov-type change of measure (see Theorem VI.2.2 in Brémaud (1981)). If E [ Z T * λ ] = 1 we obtain an equivalent measure P λ P via
P λ ( A ) : = E [ 𝟙 A Z T * λ ] A F .
Under the measure P λ , τ has intensity λ : more precisely, this means that the process
M t λ : = 𝟙 { t τ } 0 t τ λ s d s , 0 t T * ,
is a P λ -martingale.
Now we introduce a precise definition of ambiguity on the default intensity which is very much in spirit of the G-Brownian motion in Peng (2010): we consider an interval [ λ ̲ , λ ¯ ] ( 0 , ) where λ ̲ and λ ¯ denote lower (upper) bounds in the default intensity. Intuitively, we include all possible intensities lying in these bounds in our family of models P . More precisely, we define the set of density generators H by
H : = { λ : λ is F - predictable and P ( λ ̲ λ t λ ¯ , t [ 0 , T * ] ) = 1 } .
Ambiguity on the default intensity is now covered by considering the concrete family of probability measures
P : = { P λ : λ H } .
In the following, we will always consider this P . First, we observe that this set is convex.
Lemma 1.
P is a convex set.
Proof. 
Consider P λ , P λ P and α ( 0 , 1 ) . Then,
α P λ ( A ) + ( 1 α ) P λ ( A ) = E 𝟙 A ( α Z T * λ + ( 1 α ) Z T * λ ) .
Now consider the (well-defined) intensity λ , given by
0 t λ s d s : = t log α e 0 t ( 1 λ s ) d s + ( 1 α ) e 0 t ( 1 λ s ) d s ,
0 t T * . Then,
α Z T * λ + ( 1 α ) Z T * λ = Z T * λ
such that by (4), P λ P refers to an equivalent change of measure. Finally, we have to check that λ H , which means that λ [ λ ̲ , λ ¯ ] , 0 t T * : note that
t log α e 0 t ( 1 λ s ) d s + ( 1 α ) e 0 t ( 1 λ s ) d s t log α e 0 t ( 1 λ ¯ ) d s + ( 1 α ) e 0 t ( 1 λ ¯ ) d s t t ( 1 λ ¯ ) = λ ¯ t ,
and λ t λ ¯ follows. Similarly, λ ̲ λ t and the claim follows since t was arbitrary. □
Remark 2.
Intuitively, the requirement λ ̲ > 0 states that there is always a positive risk of experiencing a default, which is economically reasonable. Technically it has the appealing consequence that all considered measures in P are equivalent.
It turns out that the set of possible densities will play an important role in connection with measure changes. In this regard, we define admissible measure changes with respect to P by
A : = { λ * : λ * is positive , F - predictable and E [ Z T * λ * ] = 1 } .
The associated Radon-Nikodym derivatives Z T * λ * for λ * A are the possible Radon-Nikodym derivatives for equivalent measure changes.
Remark 3.
It is of course possible to consider an ambiguity setting more general than the specific one in (6). One possibility is to consider only a subset of P . Another possibility is to allow the bounds λ ̲ and λ ¯ to depend on time, or even on the state of the process – this latter case is important for considering affine processes under uncertainty and we refer to Fadina et al. (2019) for further details. In Section 5, we consider indeed such a more general setting.

3. Absence of Arbitrage under Ambiguity

Absence of arbitrage and the respective generalizations, no free lunch (NFL), and no free lunch with vanishing risk (NFLVR), are well established concepts when the underlying probability measure P is known and fixed. Here, we give a small set of sufficient conditions for absence of arbitrage extended to the setting with ambiguity. In this regard, consider, a fixed set P of probability measures on the measurable space ( Ω , F ) . In addition, let G = ( G t ) 0 t T * be a right-continuous filtration.
Discounted price processes of the traded assets are given by a finite dimensional G -semimartingale X = ( X t ) 0 t T * . The semimartingale property holds equivalently in any of the filtration G + or the augmentation of G + , see (Neufeld and Nutz 2014, Proposition 2.2). It is well known that then X is a semimartingale for all P P .
A self-financing trading strategy is a predictable and X-integrable process Φ and the associated discounted gains process is given by the stochastic integral of Φ with respect to X,
( Φ · X ) t = 0 t Φ u d X u , 0 t T * .
Intuitively, an arbitrage is an admissible self-financing trading strategy which starts from zero initial wealth, has non-negative pay-off under all possible future scenarios, hence for all P P , there is at least one P P , such that the pay-off is positive. This is formalized in the following definition, compare for example Vorbrink (2014). As usual a trading strategy is called a-admissible, if ( Φ · X ) t a for all 0 t T * .
Definition 1.
A self-financing trading strategy Φ is called P -arbitrage if it is a-admissible for some a > 0 and
(i) 
for every P P we have that ( Φ · X ) T * 0 , P -almost surely, and
(ii) 
for at least one P P it holds that P ( ( Φ · X ) T * > 0 ) > 0 .
Since all probability measures P P are considered as possible, a P -arbitrage is a riskless trading strategy for all possible models (i.e., for all P P ) while it is a profitable strategy for at least one scenario (i.e., for at least one P P ).
The main tool for classifying arbitrage free markets will be local martingale measures, even in the setting with ambiguity. In this regard, we call a probability measure Q a local martingale measure if X is a Q -local martingale.
It is well-known that no arbitrage or, more precisely, no free lunch with vanishing risk (NFLVR) in a market where discounted price processes are locally bounded semimartingales is equivalent to the existence of an equivalent local martingale measure (ELMM), see Delbaen and Schachermayer (1994, 1998). The technically difficult part of this result is to show that a precise criterion of absence of arbitrage implies the existence of an ELMM. In the following, we will not aim at such a deep result under ambiguity, but utilize the easy direction, namely that existence of an ELMM implies the absence of arbitrage as formulated below.
From the classical fundamental theorem of asset pricing (FTAP), the following result follows easily.
Theorem 1.
If, for every P P there exists an equivalent local martingale measure Q ( P ) , then there is no arbitrage in the sense of Definition 1.
Proof. 
Indeed, assume on the contrary that there is an arbitrage Φ with respect to some measure P P which we fix for the remainder of the proof. If there exists an ELMM Q ( P ) then Φ would be an arbitrage strategy together with an ELMM, a contradiction to the classical FTAP. □
This (sufficient) condition directly corresponds to the existing results in the literature (see, for example, Biagini et al. (2017)) where arbitrages of the first kind are studied under the additional assumption of continuity for the traded assets.

4. Ambiguity on the Default Intensity

Our aim is to discuss dynamic term structure models under default risk with ambiguity on the default intensity. The relevance of this issue has, for example, already been reported in Riedel (2015). Here, we take this as motivation to propose a precise framework taking ambiguity on the default intensity into account.

4.1. Dynamic Defaultable Term Structures

We specialize the considerations of absence of arbitrage in Section 3 to defaultable bond markets. Recall that we have a filtration G at hand and that τ is the G -stopping time at which the company defaults. We define the default indicator process H by
H t = 𝟙 { t τ } , 0 t T * .
The associated survival process is 1 H . A credit risky bond with a maturity time T T * is a contingent claim promising to pay one unit of currency at T. We denote the price of such a bond at time t T by P ( t , T ) . If no default occurs prior to T, P ( T , T ) = 1 . We will first consider zero recovery, i.e., assume that the bond loses its total value at default. Then P ( t , T ) = 0 on { t τ } .
Besides zero recovery, we only make the weak assumption that bond-prices prior to default are positive and absolutely continuous with respect to maturity T. This follows the well-established approach by Heath et al. (1992). More formally, we assume that
P ( t , T ) = 𝟙 { τ > t } exp t T f ( t , u ) d u 0 t T .
The initial forward curve T f ( 0 , T ) is then assumed to be sufficiently integrable and the forward rate processes f ( · , T ) are assumed to follow Itô processes satisfying
f ( t , T ) = f ( 0 , T ) + 0 t a ( s , T ) d s + 0 t b ( s , T ) d W s ,
for 0 t T . Recall that W was chosen to be a Brownian motion. We denote by O the optional σ -algebra and by B the Borel σ -algebra.
Assumption 1.
We require the following technical assumptions:
(i) 
the initial forward curve is measurable, and integrable on [ 0 , T * ] :
0 T * | f ( 0 , u ) | d u < ,
(ii) 
the drift parameter a ( ω , s , t ) is R -valued O B -measurable and integrable on [ 0 , T * ] :
0 T * 0 T * | a ( s , t ) | d s d t < ,
(iii) 
the volatility parameter b ( ω , s , t ) is R d -valued, O B -measurable, and
sup s , t T * b ( s , t ) < .
(iv) 
Let r t be the short rate process at time t, for 0 t T * . With probability one, it holds that
0 < f ( t , t ) r t , 0 t T * .
Set for 0 t T T * ,
a ¯ ( t , T ) = t T a ( t , u ) d u , b ¯ ( t , T ) = t T b ( t , u ) d u .
Lemma 2.
Under Assumption 1 it holds that,
t T f ( t , u ) d u = 0 T f ( 0 , u ) d u + 0 t a ¯ ( · , u ) d u + 0 t b ¯ ( · , u ) d W u 0 t f ( u , u ) d u
for 0 t T T * , almost surely.
This follows as in Heath et al. (1992): for the case W is a Brownian motion, this is Lemma 6.1 in Filipović (2009). This result could also be generalized where W is replaced by a semimartingale with absolutely continuous characteristics, see Proposition 5.2 in Björk et al. (1997). Note that the strong condition (iii) of uniform boundedness of b in Assumption 1 is needed for the application of the stochastic Fubini theorem.

4.2. Absence of Arbitrage without Ambiguity on the Default Intensity

The first step towards the study of term structure models with default ambiguity is the study of absence of arbitrage in (classical) intensity-based dynamic term structure models. Consider λ = ( λ t ) 0 t T * A and the probability measure P λ . Then, the dual predictable projection H p of H is given by H t p = 0 t τ λ s d s (under P λ ) . Moreover, the Doob-Meyer decomposition yields that
M λ : = H 0 · τ λ s d s
is P λ -martingale, compare equation (5).
For discounting, we use the bank account. Its value is given by a stochastic process starting with 1 which is then upcounted by the short rate r, i.e., the value process of the bank account is γ ( t ) = exp ( 0 t r s d s ) with an G -predictable process r.
In the bond market context considered here, a measure Q is called local martingale measure if, for any maturity T ( 0 , T * ] , the discounted bond price process for the bond with maturity T is a Q -local martingale. Then, we obtain the following result.
Proposition 1.
Assume that Assumption 1 holds. Consider a measure Q on ( Ω , F ) , such that M λ is a Q -martingale, that W is a Q -Brownian motion and that Q ( 0 T * | r s | d s < ) = 1 . Then Q is a local martingale measure if and only if
(i) 
f ( t , t ) = r t + λ t ,
(ii) 
the drift condition
a ¯ ( t , T ) = 1 2 b ¯ ( t , T ) 2 ,
holds d t d Q -almost surely for 0 t T T * on { τ > t } .
Proof. 
We set E = 1 H and F ( t , T ) = exp t T f ( t , u ) d u . Then (7) can be written as P ( t , T ) = E ( t ) F ( t , T ) . Integrating by part yields
d P ( t , T ) = F ( t , T ) d E ( t ) + E ( t ) d F ( t , T ) + d [ E , F ( · , T ) ] t .
For { t < τ } ,
d P ( t , T ) = P ( t , T ) λ t d t + f ( t , t ) + 1 2 b ¯ ( t , T ) 2 a ¯ ( t , T ) d t P ( t , T ) d M λ + b ¯ ( t , T ) d W t .
The discounted bond price process is a local martingale if and only if the predictable part in the semimartingale decomposition vanishes, i.e.,
f ( t , t ) r t λ t a ¯ ( t , T ) + 1 2 b ¯ ( t , T ) 2 = 0 .
Letting T = t we obtain (i) and (ii) and the result follows. □

4.3. Absence of Arbitrage with Ambiguity on the Default Intensity

Next, we derive the no-arbitrage conditions for the forward rate in term of the intensity and the short rate, and also the conditions for the drift and volatility parameters, under ambiguity on the default intensity. In this regard, we require a bit more structure: we assume that the setting detailed in Section 2 holds, in particular, we consider the family of probability measures P constructed in equation (6). Recall that for all P λ P , W is a Brownian motion and that P λ P . We may, for the moment, safely assume that the market filtration F satisfies the usual conditions under P .
For a generic real-valued, F -progressive process θ = ( θ t ) t 0 , let the process z θ = ( z t θ ) 0 t T * be given as the unique strong solution of
d z t θ = θ t z t θ d W t , z 0 θ = 1 .
Then, z θ is a continuous local martingale. If E [ z T * θ ] = 1 , we can define a probability measure P ˜ θ by letting
P ˜ θ ( A ) : = E [ 𝟙 A z T * θ ] , A F ,
just as in equation (4). Under P ˜ θ the process W ˜ = W 0 · θ s d s is a Brownian motion, i.e., W itself became a Brownian motion with drift θ , see Theorem 5.1 in Chapter 3 of Karatzas and Shreve (1998).
Moreover, set λ ˜ t : = ( f ( t , t ) r t ) · λ t 1 , t [ 0 , T * ] . Note that under Assumption 1, λ ˜ is positive (which is necessary for an equivalent change of measure). The associated density is abbreviated by
Z T * * = Z T * λ ˜ .
Theorem 2.
Consider P λ P . Under Assumption 1, there exists an ELMM to P λ , if there exists an F -progressive process θ * such that
(i) 
E P λ [ Z T * * z T * θ * ] = 1 ,
(ii) 
the drift condition
a ¯ ( t , T ) = 1 2 b ¯ ( t , T ) 2 b ¯ ( t , T ) θ t * , 0 t T T *
holds d t d P λ -almost surely on { t < τ } .
Intuitively, the theorem states that for the probability measure P λ , we find an ELMM if we are able to perform an equivalent change of measure (condition (i)) in such a way that under the new measure the drift condition holds for the Brownian motion with drift θ * (condition (ii)).
Proof. 
We start from some P λ P and fix this measure in the following. This means that, under P λ , W is a Brownian motion and τ has intensity λ . In the search for an ELMM we are looking for an equivalent measure P * which satisfies the conditions of Proposition 1.
In this regard, note the following: by its definition, (11), together with condition (i), z θ * is a density process for a change of measure via the Girsanov theorem for Itô processes, see Theorem 5.1 in Chapter 3 of Karatzas and Shreve (1998). Moreover, by (3) together with (i), Z * = Z λ ˜ is the density for the change in intensity from λ under P λ to the intensity given by λ t * : = ( f ( t , t ) r t ) 0 t T * , see (Brémaud 1981, Theorem VI.2.T2). We set
d P * : = Z T * * z T * θ * d P λ .
According to Theorem 3.40 in Chapter III of Jacod and Shiryaev (2003), this refers to a Girsanov-type (and equivalent) change of measure. Moreover, W t * = W t 0 t θ s * d s , 0 t T * , is a P * -Brownian motion and M λ * is a P * -martingale.
We now show that P * is also a local martingale measure. Recall from (9) that, under P λ ,
d P ( t , T ) P ( t , T ) = f ( t , t ) + 1 2 b ¯ ( t , T ) 2 a ¯ ( t , T ) d t d H t b ¯ ( t , T ) d W t ,
for { t < τ } ; here H t = 𝟙 { t τ } is the default indicator. We introduce the martingales W * and M λ * into this representation: note that
d P ( t , T ) P ( t , T ) = λ * ( t ) + f ( t , t ) + 1 2 b ¯ ( t , T ) 2 a ¯ ( t , T ) b ¯ ( t , T ) θ t * d t d M t λ * b ¯ ( t , T ) d W t * ,
again on { t < τ } . Since, by assumption, λ t * + f ( t , t ) = r t , together with the drift condition (ii), we obtain for the discounted bond price process P ˜ ( t , T ) = P ( t , T ) / γ t ,
d P ˜ ( t , T ) = P ˜ ( t , T ) · d M t λ * b ¯ ( t , T ) d W t * ,
which is a P * -local martingale and the proof is finished. □

5. Examples

The setting proposed above can, in dimension one, be directly linked to a special case of the non-linear affine processes introduced in Fadina et al. (2019). Indeed, note that for a progressive process λ , the integral
X t : = x + 0 t λ s d s , 0 t T *
is a special semimartingale. Moreover, there are affine bounds on drift and volatility (the bound of the volatility is simply zero) since
λ ̲ λ t λ ¯ ,
such that X is a non-linear affine process.
The major advantage of this setting is that numerical methods via non-linear PDE come into reach. More precisely, Theorem 4.1 in Fadina et al. (2019) shows that whenever ψ is Lipschitz, the non-linear expectation
E [ ψ ( X T * ) ] : = sup P P E P [ ψ ( X T * ) ]
can be expressed as viscosity solution of the fully non-linear PDE
t v ( t , x ) G x , x v ( t , x ) = 0 on [ 0 , T * ) × [ λ ̲ , λ ¯ ] , t v ( t , x ) G ( x , v ( T , x ) = ψ ( x ) x [ λ ̲ , λ ¯ ] ,
where G is defined by
G ( x , p ) : = sup λ [ λ ̲ , λ ¯ ] { λ p }
and v ( 0 , x ) = E [ ψ ( X T * ) ] (the dependency on x arises through X 0 = x ).
Clearly, when p is either strictly positive (hence x v ( t , x ) ) or negative, then the supremum in (15) is immediate and the PDE (14) can be solved using standard methods. This means that the solution to the non-linear expectation is obtained simply by the upper bound λ ¯ (or the lower bound λ ̲ , respectively). Such a condition holds if ψ is monotone. The more general case has to be solved using numerical methods and we provide a simple example now.
Example 1.
Consider a butterfly on X T * , i.e., the derivative with the payoff
ψ ( x ) = ( x K 1 ) + 2 ( x K 2 ) + + ( x K 3 ) + ,
where we choose K 1 = 0.2 , K 2 = 0.3 , and K 3 = 0.8 . Moreover, let λ ̲ = 0.1 and λ ¯ = 0.5 . Then the upper and lower price bounds for the butterfly are shown in Figure 1 (the upper prices are given by the nonlinear expectation in equation (13), while the lower prices are obtained by replacing the supremum in (13) by an infimum).

6. Ambiguity on the Recovery

A detailed study of bond markets beyond zero recovery is often neglected, the high degree of uncertainty about the recovery mechanism being a prime reason for this. This motivates us to take some time for developing a deeper understanding of a suitable recovery model under ambiguity.
We start from the observation that intensity-based models always need certain recovery assumptions, as for example, zero recovery, fractional recovery of treasury, and fractional recovery of par value, see (Bielecki and Rutkowski 2002, Chapter 8). We have so far considered the case where the credit risky bond becomes worthless at default (zero recovery). In the following, we will consider fractional recovery of market value where the credit risky bond looses a fraction of its market value upon default. Other recovery models can be treated in a similar fashion.

6.1. Fractional Recovery without Ambiguity

Fractional recovery of market value (RMV) is specified through a marked point process ( T n , R n ) n 1 where the stopping times ( T n ) denote the default times and R n ( 0 , 1 ] denotes the associated fractional recovery. Let
R t = T n t R n , 0 t T * ,
denote the recovery process. Then, R is non-increasing, positive with R 0 = 1 . The recovery process replaces the default indicator in (7). More precisely, we assume that the family of defaultable bond prices under RMV satisfy
P R ( t , T ) = R t exp t T f ( t , u ) d u , 0 t T T * .
Remark 4.
If a default occurs at t = T n , the bond looses a random fraction q t = 1 R T n [ 0 , 1 ) of its pre-default value. Thus, the value ( 1 q t ) P ( t , T ) is immediately available to the bond owner at default. It is still subject to default risk because of the possible future defaults occurring at T n + 1 , T n + 2 , .
First, we state a generalization of Proposition 1 to this setting. To this end, we require more structure and continue in the setting of the Section 2. Assume that the marked point process ( T n , R n ) n 1 is independent from W and standard in the following sense: the random times ( T n ) are the jumping times from a Poisson process with intensity one, and the recovery values ( R n ) are independent from ( T n ) and W, and uniformly distributed in [ r ̲ , r ¯ ] ( 0 , 1 ] .
The filtration G = ( G t ) 0 t T * is obtained by a progressive enlargement of F with default and recovery information (given by R), i.e.,
G t = ϵ > 0 σ ( R s , W s : 0 s t + ϵ ) , 0 t T * .
We assume again F = G T * . As next step, we introduce measure changes for the marked point process. Let
Φ t = T n t R n , 0 t T * .
This implies that the defaultable bond prices under RMV, P R ( 0 , t ) , at maturity time t of this bond, pays out Φ t , the accumulated fractional default losses. Then, Φ is a special semimartingale w.r.t. G . Let
μ Φ ( d t , d x ) = n 1 δ ( T n , R n ) ( d t , d x )
denote the associated jump measure and let ν Φ ( d t , d x ) denote its compensator, see Chapter II.1 in Jacod and Shiryaev (2003) or Chapter VIII.1 in Brémaud (1981). Note that ν Φ ( d t , d x ) = 𝟙 { x [ r ̲ , r ¯ ] } ( r ¯ r ̲ ) 1 d x d t .
We introduce the densities
L μ , h = T n T * μ T n h ( T n , R n ) · e 0 T * r ̲ r ¯ ( 1 μ t h ( t , x ) ) ( r ¯ r ̲ ) 1 d x d t ,
where the predictable process μ is positive and, for any x [ r ̲ , r ¯ ] , the G -predictable process ( h ( t , x ) ) 0 t T * is also positive. If E [ L T * ] = 1 , we can define the equivalent measure P μ , h by
d P μ , h = L T * d P .
By A * we denote all pairs ( μ , h ) which satisfy the above properties. Then, the compensator of the jump measure μ Φ under P μ , h is
μ t h ( t , x ) ν Φ ( d t , d x ) = μ t h ( t , x ) 𝟙 { x [ r ̲ , r ¯ ] } ( r ¯ r ̲ ) 1 d x d t = : K t μ , h ( d x ) d t ,
see T10 in Section VIII.3 of Brémaud (1981). Next, we compute the compensator of R. We obtain from (16) that
R t 0 t R s ( x 1 ) K s μ , h ( d x ) d s , 0 t T *
is a P μ , h -martingale. For a G -progressive process g, we denote
M g = R + 0 · R s g s d s .
Proposition 2.
Assume that Assumption 1 holds and let g be a positive and G -predictable process. Consider a measure Q on ( Ω , F ) , such that M g is a Q -martingale, W is a Q -Brownian motion, and Q ( 0 T * | r s | d s < ) = 1 . Then Q is a local martingale measure if and only if
(i) 
f ( t , t ) = r t + g t ,
(ii) 
the drift condition
a ¯ ( t , T ) = 1 2 b ¯ ( t , T ) 2 , 0 t T T * ,
holds d t d Q -almost surely.
Proof. 
We generalize the proof of Proposition 1 to the case of RMV. To this end, let F ( t , T ) = exp t T f ( t , u ) d u . Then (17) reads P R ( t , T ) = R ( t ) F ( t , T ) . Integrating by part yields
d P R ( t , T ) = F ( t , T ) d R ( t ) + R ( t ) d F ( t , T ) + d [ R , F ( · , T ) ] t .
Note that, by assumption,
M t g = R t + 0 t R s g s d s , 0 t T * ,
is a Q-martingale and that [ R , F ( . , T ) ] = 0 since R is of finite variation and F ( . , T ) is continuous. Hence, by Lemma 2,
d P R ( t , T ) = P R ( t , T ) g t + f ( t , t ) + 1 2 b ¯ ( t , T ) 2 a ¯ ( t , T ) d t + P R ( t , T ) r t d t + d M g b ¯ ( t , T ) d W t ,
and we obtain the result as in the proof of Proposition 1. □
Example 2.
A classical example is when the defaults ( T n ) arrive at rate λ > 0 , and the recovery values ( R n ) are i.i.d. Then, ( x 1 ) ν Φ ( d t , d x ) = λ E [ R 1 1 ] d t . We obtain that the instantaneous forward rate of the defaultable bond f ( t , t ) equals r t + λ E [ 1 R 1 ] . In the case of zero recovery, we recover f ( t , t ) = r t + λ , and, in the case of full recovery (the case without default risk), f ( t , t ) = r t .

6.2. Fractional Recovery with Ambiguity

We introduce ambiguity in this setting by changing from the standardized measure P to various appropriate measures via the Girsanov theorem. We also generalize the setting for ambiguity from the quite specific P to a general set of probability measures P * here, see Remark 3. The reason for this is also economic: while bounding the intensity from above and below seems to be quite plausible, an upper/lower bound on the recovery (i.e., on ( R n ) ) sounds too strong for some applications.
Recall that A * was the set of all candidates ( μ , h ) which induce the measure changes via (19). Ambiguity is introduced by the set P * of probability measures satisfying
P * { P μ , h : ( μ , h ) A * } .
If P * contains only one probability measure, we are in the classical setting, otherwise there is ambiguity in the market. Measure changes from P μ , h to a new measure are done via the density L μ * , h * (see (18)) where, as above, μ * , h * ( . , x ) , x [ r ̲ , r ¯ ] are positive and progressive. Recall the definition of the density z θ in (11).
Theorem 3.
Let g t * : = f ( t , t ) r t , t [ 0 , T * ] , and assume that Assumption 1 holds. Then there exists an ELMM for P μ , h P * if
(i) 
there exists an F -progressive θ * such that E [ z T * θ * ] = 1 ,
(ii) 
there exist μ * and h * ( t , x ) , such that E [ L μ * , h * ] = 1 and
g t * = ( x 1 ) μ t * h * ( t , x ) K t μ , h ( d x ) , 0 t T * ,
d t d P -almost surely, and
(iii) 
the drift condition
a ¯ ( t , T ) = 1 2 b ¯ ( t , T ) 2 b ¯ ( t , T ) θ t * , 0 t T T * ,
holds d t d P -almost surely.
Absence of arbitrage in this general ambiguity setting can now be classified, thanks to Theorem 1 as follows: if an ELMM exists for each P μ , h P * , then the market is free of arbitrage in the sense of Definition 1.
Proof. 
Fix P μ , h P * . We can define an equivalent measure P * P μ , h by
d P * : = L μ * , h * z T * θ * d P μ , h ,
with μ * and h * as in (ii). According to Theorem 3.40 in Chapter III of Jacod and Shiryaev (2003), this refers to a Girsanov-type (and equivalent) change of measure. Moreover, W * = W 0 · θ s * d s is a P * -Brownian motion. Next, note that the compensator of the jump measure μ Φ under P * computes, according to T10 in Section VIII.3 in Brémaud (1981), to
ν * ( d t , d x ) : = μ t * h * ( t , x ) K t μ , h ( d x ) d t
with K t μ , h from Equation (20). This implies that
M t g * = R s ( x 1 ) ν * ( d s , d x ) = R t 0 t R s g s * d s
is a P * -martingale.
Now, we show that P * is indeed a martingale measure: from (22) we obtain that
d P R ( t , T ) P R ( t , T ) = f ( t , t ) + 1 2 b ¯ ( t , T ) 2 a ¯ ( t , T ) d t + d R t b ¯ ( t , T ) d W t .
It follows that
d P R ( t , T ) P R ( t , T ) = f ( t , t ) + 1 2 b ¯ ( t , T ) 2 a ¯ ( t , T ) b ¯ ( t , T ) θ t * d t + d M t g b ¯ ( t , T ) d W t * = r t d t + d M t g b ¯ ( t , T ) d W t * ,
by the definition of g * and the drift condition (iii). Hence, discounted bond prices are P * -local martingales and the proof is finished. □
Remark 5.
We can view zero recovery in the above setting by assuming that P ( R 1 = 0 ) = 1 and letting τ = T 1 . Note that this case is excluded in RMV setting, since, under this assumption, at the first default all prices drop to zero and further defaults can not occur.

7. Conclusions

This paper provides a first step towards including ambiguity in intensity based models for credit risk. Many research questions are still open: first, the extension of constant boundaries λ ̲ , λ ¯ to time-dependent, or, as in Fadina et al. (2019), state-dependent boundaries. Second, the extension to two or more defaultable assets, where default dependence comes into play. Third, estimation and determination of the bounds by statistical methods or calibration to market data. We hope that our paper provides the foundation for future works in these directions.

Author Contributions

All authors contributed equally to the paper. The author appearance is in alphabetical order.

Funding

Financial support by Carl-Zeiss-Stiftung, and German Research Foundation (DFG) via CRC 1283 is gratefully acknowledged. We also thank the Freiburg Institute of Advanced Studies (FRIAS) for its hospitality and financial support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Aksamit, Anna, and Monique Jeanblanc. 2017. Enlargements of Filtration with Finance in View. SpringerBriefs in Quantitative Finance. Berlin: Springer. [Google Scholar]
  2. Artzner, Philippe, and Freddy Delbaen. 1995. Default risk insurance and incomplete markets. Mathematical Finance 5: 187–95. [Google Scholar] [CrossRef]
  3. Avellaneda, Marco, Arnon Levy, and Antonio Parás. 1995. Pricing and hedging derivative securities in markets with uncertain volatilities. Appl. Math. Finance 2: 73–88. [Google Scholar] [CrossRef]
  4. Bayraktar, Erhan, and Yuchong Zhang. 2013. Fundamental theorem of asset pricing under transaction costs and model uncertainty. arXiv arXiv:1309.1420v2. [Google Scholar] [CrossRef]
  5. Bélanger, Alain, Steven E. Shreve, and Dennis Wong. 2004. A general framework for pricing credit risk. Mathematical Finance 14: 317–50. [Google Scholar] [CrossRef]
  6. Biagini, Francesca, and Yinglin Zhang. 2017. Reduced-form framework under model uncertainty. arXiv arXiv:1707.04475. [Google Scholar]
  7. Biagini, Sara, Bruno Bouchard, Constantinos Kardaras, and Marcel Nutz. 2017. Robust fundamental theorem for continuous processes. Mathematical Finance 27: 963–87. [Google Scholar] [CrossRef]
  8. Bielecki, Tomasz, and Marek Rutkowski. 2002. Credit Risk: Modeling, Valuation and Hedging. Berlin, Heidelberg and New York: Springer. [Google Scholar]
  9. Björk, Tomas, Giovanni Di Masi, Yuri Kabanov, and Wolfgang Runggaldier. 1997. Towards a general theory of bond markets. Finance and Stochastics 1: 141–74. [Google Scholar] [CrossRef]
  10. Black, Fischer, and John C. Cox. 1976. Valuing corporate securities: some effects of bond indenture provisions. Journal of Finance 31: 351–67. [Google Scholar] [CrossRef]
  11. Bouchard, Bruno, and Marcel Nutz. 2015. Arbitrage and duality in nondominated discrete-time models. The Annals of Applied Probability 25: 823–59. [Google Scholar] [CrossRef]
  12. Brémaud, Pierre. 1981. Point Processes and Queues. Berlin, Heidelberg and New York: Springer. [Google Scholar]
  13. Burzoni, Matteo, Frank Riedel, and Mete H. Soner. 2017. Viability and arbitrage under knightian uncertainty. arXiv arXiv:1707.03335. [Google Scholar] [CrossRef]
  14. Delbaen, Freddy, and Walter Schachermayer. 1994. A general version of the fundamental theorem of asset pricing. Mathematische Annalen 300: 463–520. [Google Scholar] [CrossRef]
  15. Delbaen, Freddy, and Walter Schachermayer. 1998. The fundamental theorem of asset pricing for unbounded stochastic processes. Mathematische Annalen 312: 215–50. [Google Scholar] [CrossRef]
  16. Denis, Laurent, and Claude Martini. 2006. A theoretical framework for the pricing of contingent claims in the presence of model uncertainty. The Annals of Applied Probability 16: 827–52. [Google Scholar] [CrossRef]
  17. Duffie, Darrell, and David Lando. 2001. Term structures of credit spreads with incomplete accounting information. Econometrica 69: 633–64. [Google Scholar] [CrossRef]
  18. Duffie, Darrell, Mark Schroder, and Costis Skiadas. 1996. Recursive valuation of defaultable securities and the timing of resolution of uncertainty. The Annals of Applied Probability 6: 1075–90. [Google Scholar] [CrossRef]
  19. Ellsberg, Daniel. 1961. Risk, ambiguity, and the Savage axioms. The Quarterly Journal of Economics 75: 643–69. [Google Scholar] [CrossRef]
  20. Fadina, Tolulope, Ariel Neufeld, and Thorsten Schmidt. 2019. Affine processes under parameter uncertainty. Probability, Uncertainty and Quantitative Risk 4: 1–35. [Google Scholar] [CrossRef]
  21. Filipović, Damir. 2009. Term Structure Models: A Graduate Course. Berlin, Heidelberg and New York: Springer. [Google Scholar]
  22. Fontana, Claudio, and Thorsten Schmidt. 2018. General dynamic term structures under default risk. Stochastic Processes and their Applications 128: 3353–86. [Google Scholar] [CrossRef]
  23. Frey, Rüdiger, and Thorsten Schmidt. 2009. Pricing corporate securities under noisy asset information. Mathematical Finance 19: 403–21. [Google Scholar] [CrossRef]
  24. Frey, Rüdiger, and Thorsten Schmidt. 2011. Filtering and Incomplete Information. In Credit Risk Frontiers. Edited by T. Bielecki, D. Brigo and F. Patras. Hoboken: Wiley. [Google Scholar]
  25. Frey, Rüdiger, and Thorsten Schmidt. 2012. Pricing and hedging of credit derivatives via the innovations approach to nonlinear filtering. Finance and Stochastics 16: 105–133. [Google Scholar] [CrossRef]
  26. Gehmlich, Frank, and Thorsten Schmidt. 2018. Dynamic defaultable term structure modelling beyond the intensity paradigm. Mathematical Finance 28: 211–39. [Google Scholar] [CrossRef]
  27. Geske, Robert. 1977. The valuation of corporate liabilities as compound options. Journal of Financial and Quantitative Analysis 12: 541–52. [Google Scholar] [CrossRef]
  28. Heath, David, Robert Jarrow, and Andrew Morton. 1992. Bond pricing and the term structure of interest rates. Econometrica 60: 77–105. [Google Scholar] [CrossRef]
  29. Jacod, Jean, and Albert N Shiryaev. 2003. Limit Theorems for Stochastic Processes, 2nd ed. Berlin: Springer. [Google Scholar]
  30. Karatzas, Ioannis, and Steven E Shreve. 1998. Methods of Mathematical Finance. New York: Springer. [Google Scholar]
  31. Kim, In Joon, Krishna Ramaswamy, and Suresh Sundaresan. 1993. The valuation of corporate fixed income securities. Working paper. Philadelphia: Wharton School, University of Pennsylvania. [Google Scholar]
  32. Knight, Frank H. 1921. Risk, Uncertainty and Profit. Boston: Houghton Mifflin. [Google Scholar]
  33. Lando, David. 1994. Three Essays on Contingent Claim Pricing. Ph.D. thesis, Cornell University, Ithaca, NY, USA. [Google Scholar]
  34. Leland, Hayne E., and Klaus Bjerre Toft. 1996. Optimal capital structure, endogenous bankruptcy, and the term structure of credit spreads. Journal of Finance 51: 987–1019. [Google Scholar] [CrossRef]
  35. Lyons, Terry. J. 1995. Uncertain volatility and the risk-free synthesis of derivatives. Appl. Math. Finance 2: 117–133. [Google Scholar] [CrossRef]
  36. Merton, Robert. 1974. On the pricing of corporate debt: the risk structure of interest rates. Journal of Finance 29: 449–70. [Google Scholar]
  37. Neufeld, Ariel, and Marcel Nutz. 2014. Measurability of semimartingale characteristics with respect to the probability law. Stochastic Processes and their Applications 124: 3819–45. [Google Scholar] [CrossRef]
  38. Peng, Shige. 2010. Nonlinear expectations and stochastic calculus under uncertainty. arXiv arXiv:1002.4546. [Google Scholar]
  39. Riedel, Frank. 2015. Financial economics without probabilistic prior assumptions. Decisions in Economics and Finance 38: 75–91. [Google Scholar] [CrossRef]
  40. Sylla, Richard. 2002. An Historical Primer on the Business of Credit Rating. Boston: Springer, pp. 19–40. [Google Scholar]
  41. Vorbrink, Jörg. 2014. Financial markets with volatility uncertainty. Journal of Mathematical Economics 53: 64–78. [Google Scholar] [CrossRef] [Green Version]
1
The risk that an agent fails to fulfil contractual obligations. Example of an instrument bearing credit risk is a corporate bond.
2
The amount that the owner of a defaulted claim receives upon default.
3
For a historical account, see Sylla (2002): John Moody founded the first rating agency in 1909, in the United States.
4
Augmentation can be done in a standard fashion with respect to P .
Figure 1. This figure shows the solution of the nonlinear PDE in Equation (14) with boundary condition ψ ( y ) = ( y K 1 + x ) + 2 ( y K 2 + x ) + + ( y K 3 + x ) + , K 1 = 0.2 , K 2 = 0.3 , K 3 = 0.8 , and x [ 0.5 , 0.7 ] is depicted on the x-axis of the plot. The dashed lines show the solution for the lower bound (upper bound, respectively), i.e., for the constants λ ̲ = 0.1 and λ ¯ = 0.5 . The upper and lower solid lines show the upper and lower price bounds.
Figure 1. This figure shows the solution of the nonlinear PDE in Equation (14) with boundary condition ψ ( y ) = ( y K 1 + x ) + 2 ( y K 2 + x ) + + ( y K 3 + x ) + , K 1 = 0.2 , K 2 = 0.3 , K 3 = 0.8 , and x [ 0.5 , 0.7 ] is depicted on the x-axis of the plot. The dashed lines show the solution for the lower bound (upper bound, respectively), i.e., for the constants λ ̲ = 0.1 and λ ¯ = 0.5 . The upper and lower solid lines show the upper and lower price bounds.
Risks 07 00064 g001

Share and Cite

MDPI and ACS Style

Fadina, T.; Schmidt, T. Default Ambiguity. Risks 2019, 7, 64. https://0-doi-org.brum.beds.ac.uk/10.3390/risks7020064

AMA Style

Fadina T, Schmidt T. Default Ambiguity. Risks. 2019; 7(2):64. https://0-doi-org.brum.beds.ac.uk/10.3390/risks7020064

Chicago/Turabian Style

Fadina, Tolulope, and Thorsten Schmidt. 2019. "Default Ambiguity" Risks 7, no. 2: 64. https://0-doi-org.brum.beds.ac.uk/10.3390/risks7020064

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop