Small‐time,large‐time,and asymptotics for the Rough Heston model

We characterize the behavior of the Rough Heston model introduced by Jaisson and Rosenbaum (2016, Ann. Appl. Probab., 26, 2860–2882) in the small‐time, large‐time, and (i.e., ) limits. We show that the short‐maturity smile scales in qualitatively the same way as a general rough stochastic volatility model , and the rate function is equal to the Fenchel–Legendre transform of a simple transformation of the solution to the same Volterra integral equation (VIE) that appears in El Euch and Rosenbaum (2019, Math. Financ., 29, 3–38), but with the drift and mean reversion terms removed. The solution to this VIE satisfies a space–time scaling property which means we only need to solve this equation for the moment values of and so the rate function can be efficiently computed using an Adams scheme or a power series, and we compute a power series in the log‐moneyness variable for the asymptotic implied volatility which yields tractable expressions for the implied vol skew and convexity which is useful for calibration purposes. We later derive a formal saddle point approximation for call options in the Forde and Zhang (2017) large deviations regime which goes to higher order than previous works for rough models. Our higher‐order expansion captures the effect of both drift terms, and at leading order is of qualitatively the same form as the higher‐order expansion for a general model which appears in Friz et al. (2018, Math. Financ., 28, 962–988). The limiting asymptotic smile in the large‐maturity regime is obtained via a stability analysis of the fixed points of the VIE, and is the same as for the standard Heston model in Forde and Jacquier (2011, Finance Stoch., 15, 755–780). Finally, using Lévy's convergence theorem, we show that the log stock price tends weakly to a nonsymmetric random variable as (i.e., ) whose moment generating function (MGF) is also the solution to the Rough Heston VIE with , and we show that tends weakly to a nonsymmetric random variable as , which leads to a nonflat nonsymmetric asymptotic smile in the Edgeworth regime, where the log‐moneyness as , and we compute this asymptotic smile numerically. We also show that the third moment of the log stock price tends to a finite constant as (in contrast to the Rough Bergomi model discussed in Forde et al. (2020, Preprint) where the skew flattens or blows up) and the process converges on pathspace to a random tempered distribution which has the same law as the hyper‐rough Heston model, discussed in Jusselin and Rosenbaum (2020, Math. Finance, 30, 1309–1336) and Abi Jaber (2019, Ann. Appl. Probab., 29, 3155–3200).     

Hedging nontradable risks with transaction costs and price impact

A risk‐averse agent hedges her exposure to a nontradable risk factor U using a correlated traded asset S and accounts for the impact of her trades on both factors. The effect of the agent's trades on U is referred to as cross‐impact. By solving the agent's stochastic control problem, we obtain a closed‐form expression for the optimal strategy when the agent holds a linear position in U. When the exposure to the nontradable risk factor is nonlinear, we provide an approximation to the optimal strategy in closed‐form, and prove that the value function is correctly approximated by this strategy when cross‐impact and risk‐aversion are small. We further prove that when is nonlinear, the approximate optimal strategy can be written in terms of the optimal strategy for a linear exposure with the size of the position changing dynamically according to the exposure's “Delta” under a particular probability measure.     

OPTIMAL SELLING RULES FOR MONETARY INVARIANT CRITERIA: TRACKING THE MAXIMUM OF A PORTFOLIO WITH NEGATIVE DRIFT

Considering a positive portfolio diffusion X with negative drift, we investigate optimal stopping problems of the form where f is a nonincreasing function, τ is the next random time where the portfolio X crosses zero and θ is any stopping time smaller than τ. Hereby, our motivation is the obtention of an optimal selling strategy minimizing the relative distance between the liquidation value of the portfolio and its highest possible value before it reaches zero. This paper unifies optimal selling rules observed for the quadratic absolute distance criteria in this stationary framework with bang–bang type ones observed for monetary invariant criteria but in finite horizon. More precisely, we provide a verification result for the general stopping problem of interest and derive the exact solution for two classical criteria f of the literature. For the power utility criterion with , instantaneous selling is always optimal, which is consistent with previous observations for the Black‐Scholes model in finite observation. On the contrary, for a relative quadratic error criterion, , selling is optimal as soon as the process X crosses a specified function φ of its running maximum . These results reinforce the idea that optimal stopping problems of similar type lead easily to selling rules of very different nature. Nevertheless, our numerical experiments suggest that the practical optimal selling rule for the relative quadratic error criterion is in fact very close to immediate selling.     

ROBUST FUNDAMENTAL THEOREM FOR CONTINUOUS PROCESSES

We study a continuous‐time financial market with continuous price processes under model uncertainty, modeled via a family of possible physical measures. A robust notion of no‐arbitrage of the first kind is introduced; it postulates that a nonnegative, nonvanishing claim cannot be superhedged for free by using simple trading strategies. Our first main result is a version of the fundamental theorem of asset pricing: holds if and only if every admits a martingale measure that is equivalent up to a certain lifetime. The second main result provides the existence of optimal superhedging strategies for general contingent claims and a representation of the superhedging price in terms of martingale measures.     

Optimal portfolio under fractional stochastic environment

Rough stochastic volatility models have attracted a lot of attention recently, in particular for the linear option pricing problem. In this paper, starting with power utilities, we propose to use a martingale distortion representation of the optimal value function for the nonlinear asset allocation problem in a (non‐Markovian) fractional stochastic environment (for all values of the Hurst index ). We rigorously establish a first‐order approximation of the optimal value, when the return and volatility of the underlying asset are functions of a stationary slowly varying fractional Ornstein–Uhlenbeck process. We prove that this approximation can be also generated by a fixed zeroth‐ order trading strategy providing an explicit strategy which is asymptotically optimal in all admissible controls. Furthermore, we extend the discussion to general utility functions, and obtain the asymptotic optimality of this fixed strategy in a specific family of admissible strategies.     

INDIFFERENCE PRICING FOR CONTINGENT CLAIMS: LARGE DEVIATIONS EFFECTS

We study utility indifference prices and optimal purchasing quantities for a nontraded contingent claim in an incomplete semimartingale market with vanishing hedging errors. We make connections with the theory of large deviations. We concentrate on sequences of semicomplete markets where in the nth market, the claim admits the decomposition . Here, is replicable by trading in the underlying assets , but is independent of . Under broad conditions, we may assume that vanishes in accordance with a large deviations principle (LDP) as n grows. In this setting, for an exponential investor, we identify the limit of the average indifference price , for units of , as . We show that if , the limiting price typically differs from the price obtained by assuming bounded positions , and the difference is explicitly identifiable using large deviations theory. Furthermore, we show that optimal purchase quantities occur at the large deviations scaling, and hence large positions arise endogenously in this setting.     

Azéma martingales for Bessel and CIR processes and the pricing of Parisian zero‐coupon bonds

In this paper, we study the excursions of Bessel and Cox–Ingersoll–Ross (CIR) processes with dimensions . We obtain densities for the last passage times and meanders of the processes. Using these results, we prove a variation of the Azéma martingale for the Bessel and CIR processes based on excursion theory. Furthermore, we study their Parisian excursions, and generalize previous results on the Parisian stopping time of Brownian motion to that of the Bessel and CIR processes. We obtain explicit formulas and asymptotic results for the densities of the Parisian stopping times, and develop exact simulation algorithms to sample the Parisian stopping times of Bessel and CIR processes. We introduce a new type of bond, the zero‐coupon Parisian bond. The buyer of such a bond is betting against zero interest rates, while the seller is effectively hedging against a period where interest rates fluctuate around 0. Using our results, we propose two methods for pricing these bonds and provide numerical examples.     

Convergence of optimal expected utility for a sequence of discrete‐time markets

We examine Kreps' conjecture that optimal expected utility in the classic Black–Scholes–Merton (BSM) economy is the limit of optimal expected utility for a sequence of discrete‐time economies that “approach” the BSM economy in a natural sense: The nth discrete‐time economy is generated by a scaled n‐step random walk, based on an unscaled random variable ζ with mean 0, variance 1, and bounded support. We confirm Kreps' conjecture if the consumer's utility function U has asymptotic elasticity strictly less than one, and we provide a counterexample to the conjecture for a utility function U with asymptotic elasticity equal to 1, for ζ such that .     

MULTIVARIATE SUBORDINATION OF MARKOV PROCESSES WITH FINANCIAL APPLICATIONS

This paper develops the procedure of multivariate subordination for a collection of independent Markov processes with killing. Starting from d independent Markov processes with killing and an independent d‐dimensional time change , we construct a new process by time, changing each of the Markov processes with a coordinate . When is a d‐dimensional Lévy subordinator, the time changed process is a time‐homogeneous Markov process with state‐dependent jumps and killing in the product of the state spaces of . The dependence among jumps of its components is governed by the d‐dimensional Lévy measure of the subordinator. When is a d‐dimensional additive subordinator, Y is a time‐inhomogeneous Markov process. When with forming a multivariate Markov process, is a Markov process, where each plays a role of stochastic volatility of . This construction provides a rich modeling architecture for building multivariate models in finance with time‐ and state‐dependent jumps, stochastic volatility, and killing (default). The semigroup theory provides powerful analytical and computational tools for securities pricing in this framework. To illustrate, the paper considers applications to multiname unified credit‐equity models and correlated commodity models.     

OPTIMAL INVESTMENT UNDER RELATIVE PERFORMANCE CONCERNS

We consider the problem of optimal investment when agents take into account their relative performance by comparison to their peers. Given N interacting agents, we consider the following optimization problem for agent i, : where is the utility function of agent i, his portfolio, his wealth, the average wealth of his peers, and is the parameter of relative interest for agent i. Together with some mild technical conditions, we assume that the portfolio of each agent i is restricted in some subset . We show existence and uniqueness of a Nash equilibrium in the following situations:

• ‐ unconstrained agents,
• ‐ constrained agents with exponential utilities and Black–Scholes financial market.

We also investigate the limit when the number of agents N goes to infinity. Finally, when the constraints sets are vector spaces, we study the impact of the s on the risk of the market.     

On the C‐property and ‐representations of risk measures

We identify a large class of Orlicz spaces for which the topology fails the C‐property introduced by Biagini and Frittelli. We also establish a variant of the C‐property and use it to prove a ‐representation theorem for proper convex increasing functionals, satisfying a suitable version of Delbaen's Fatou property, on Orlicz spaces with . Our results apply, in particular, to risk measures on all Orlicz spaces other than .     

A variation of the Azéma martingale and drawdown options

In this paper, we derive a variation of the Azéma martingale using two approaches—a direct probabilistic method and another by projecting the Kennedy martingale onto the filtration generated by the drawdown duration. This martingale links the time elapsed since the last maximum of the Brownian motion with the maximum process itself. We derive explicit formulas for the joint densities of , which are the first time the drawdown period hits a prespecified duration, the value of the Brownian motion, and the maximum up to this time. We use the results to price a new type of drawdown option, which takes into account both dimensions of drawdown risk—the magnitude and the duration.     

ON THE MARTINGALE PROPERTY IN STOCHASTIC VOLATILITY MODELS BASED ON TIME‐HOMOGENEOUS DIFFUSIONS

Lions and Musiela give sufficient conditions to verify when a stochastic exponential of a continuous local martingale is a martingale or a uniformly integrable martingale. Blei and Engelbert and Mijatovi? and Urusov give necessary and sufficient conditions in the case of perfect correlation (). For financial applications, such as checking the martingale property of the stock price process in correlated stochastic volatility models, we extend their work to the arbitrary correlation case (). We give a complete classification of the convergence properties of both perpetual and capped integral functionals of time‐homogeneous diffusions and generalize results in Mijatovi? and Urusov with direct proofs avoiding the use of separating times (concept introduced by Cherny and Urusov and extensively used in the proofs of Mijatovi? and Urusov).     

EXPECTATIONS OF FUNCTIONS OF STOCHASTIC TIME WITH APPLICATION TO CREDIT RISK MODELING

We develop two novel approaches to solving for the Laplace transform of a time‐changed stochastic process. We discard the standard assumption that the background process () is Lévy. Maintaining the assumption that the business clock () and the background process are independent, we develop two different series solutions for the Laplace transform of the time‐changed process . In fact, our methods apply not only to Laplace transforms, but more generically to expectations of smooth functions of random time. We apply the methods to introduce stochastic time change to the standard class of default intensity models of credit risk, and show that stochastic time‐change has a very large effect on the pricing of deep out‐of‐the‐money options on credit default swaps.     

STABILITY OF THE EXPONENTIAL UTILITY MAXIMIZATION PROBLEM WITH RESPECT TO PREFERENCES

This paper studies stability of the exponential utility maximization when there are small variations on agent's utility function. Two settings are considered. First, in a general semimartingale model where random endowments are present, a sequence of utilities defined on converges to the exponential utility. Under a uniform condition on their marginal utilities, convergence of value functions, optimal payoffs, and optimal investment strategies are obtained, their rate of convergence is also determined. Stability of utility‐based pricing is studied as an application. Second, a sequence of utilities defined on converges to the exponential utility after shifting and scaling. Their associated optimal strategies, after appropriate scaling, converge to the optimal strategy for the exponential hedging problem. This complements Theorem 3.2 in [Nutz, M. (2012): Risk aversion asymptotics for power utility maximization. Probab. Theory & Relat. Fields 152, 703–749], which establishes the convergence for a sequence of power utilities.     

Optimal investment and pricing in the presence of defaults

We consider the optimal investment problem with random endowment in the presence of defaults. For an investor with constant absolute risk aversion, we identify the certainty equivalent, and compute prices for defaultable bonds and dynamic protection against default. This latter price is interpreted as the premium for a contingent credit default swap, and connects our work with earlier articles, where the investor is protected upon default. We consider a multiple risky asset model with a single default time, at which point each of the assets may jump in price. Investment opportunities are driven by a diffusion X taking values in an arbitrary region . We allow for stochastic volatility, correlation, and recovery; unbounded random endowments; and postdefault trading. We identify the certainty equivalent with a semilinear parabolic partial differential equation with quadratic growth in both function and gradient. Under minimal integrability assumptions, we show that the certainty equivalent is a classical solution. Numerical examples highlight the relationship between the factor process, market dynamics, utility‐based prices, and default insurance premium. In particular, we show that the holder of a defaultable bond has a strong incentive to short the underlying stock, even for very low default intensities.     

Hyperbolic normal stochastic volatility model

For option pricing models and heavy-tailed distributions, this study proposes a continuous-time stochastic volatility model based on an arithmetic Brownian motion: a one-parameter extension of the normal stochastic alpha-beta-rho (SABR) model. Using two generalized Bougerol's identities in the literature, the study shows that our model has a closed-form Monte Carlo simulation scheme and that the transition probability for one special case follows Johnson's distribution—a popular heavy-tailed distribution originally proposed without stochastic process. It is argued that the distribution serves as an analytically superior alternative to the normal SABR model because the two distributions are empirically similar.     

EFFICIENT PRICING OF BARRIER OPTIONS AND CREDIT DEFAULT SWAPS IN LÉVY MODELS WITH STOCHASTIC INTEREST RATE

Recently, advantages of conformal deformations of the contours of integration in pricing formulas for European options have been demonstrated in the context of wide classes of Lévy models, the Heston model, and other affine models. Similar deformations were used in one‐factor Lévy models to price options with barrier and lookback features and credit default swaps (CDSs). In the present paper, we generalize this approach to models, where the dynamics of the assets is modeled as , where X is a Lévy process, and the interest rate is stochastic. Assuming that X and r are independent, and , the infinitesimal generator of the pricing semigroup in the model for the short rate, satisfies weak regularity conditions, which hold for popular models of the short rate, we develop a variation of the pricing procedure for Lévy models which is almost as fast as in the case of the constant interest rate. Numerical examples show that about 0.15 second suffices to calculate prices of 8 options of same maturity in a two‐factor model with the error tolerance and less; in a three‐factor model, accuracy of order 0.001–0.005 is achieved in about 0.2 second. Similar results are obtained for quanto CDS, where an additional stochastic factor is the exchange rate. We suggest a class of Lévy models with the stochastic interest rate driven by 1–3 factors, which allows for fast calculations. This class can satisfy the current regulatory requirements for banks mandating sufficiently sophisticated credit risk models.     

MODEL‐INDEPENDENT LOWER BOUND ON VARIANCE SWAPS

It is well known that, under a continuity assumption on the price of a stock S, the realized variance of S for maturity T can be replicated by a portfolio of calls and puts maturing at T. This paper assumes that call prices on S maturing at T are known for all strikes but makes no continuity assumptions on S. We derive semiexplicit expressions for the supremum lower bound on the hedged payoff, at maturity T, of a long position in the realized variance of S. Equivalently, is the supremum strike K such that an investor with a long position in a variance swap with strike K can ensure a nonnegative payoff at T. We study examples with constant implied volatilities and with a volatility skew. In our examples, is close to the fair variance strike obtained under the continuity assumption.     

HIGH‐ORDER SHORT‐TIME EXPANSIONS FOR ATM OPTION PRICES OF EXPONENTIAL LÉVY MODELS

The short‐time asymptotic behavior of option prices for a variety of models with jumps has received much attention in recent years. In this work, a novel second‐order approximation for at‐the‐money (ATM) option prices is derived for a large class of exponential Lévy models with or without Brownian component. The results hereafter shed new light on the connection between both the volatility of the continuous component and the jump parameters and the behavior of ATM option prices near expiration. In the presence of a Brownian component, the second‐order term, in time‐t, is of the form , with d₂ only depending on Y, the degree of jump activity, on σ, the volatility of the continuous component, and on an additional parameter controlling the intensity of the “small” jumps (regardless of their signs). This extends the well‐known result that the leading first‐order term is . In contrast, under a pure‐jump model, the dependence on Y and on the separate intensities of negative and positive small jumps are already reflected in the leading term, which is of the form . The second‐order term is shown to be of the form and, therefore, its order of decay turns out to be independent of Y. The asymptotic behavior of the corresponding Black–Scholes implied volatilities is also addressed. Our method of proof is based on an integral representation of the option price involving the tail probability of the log‐return process under the share measure and a suitable change of probability measure under which the pure‐jump component of the log‐return process becomes a Y‐stable process. Our approach is sufficiently general to cover a wide class of Lévy processes, which satisfy the latter property and whose Lévy density can be closely approximated by a stable density near the origin. Our numerical results show that the first‐order term typically exhibits rather poor performance and that the second‐order term can significantly improve the approximation's accuracy, particularly in the absence of a Brownian component.