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.     

Local volatility under rough volatility

Several asymptotic results for the implied volatility generated by a rough volatility model have been obtained in recent years (notably in the small-maturity regime), providing a better understanding of the shapes of the volatility surface induced by rough volatility models, supporting their calibration power to SP500 option data. Rough volatility models also generate a local volatility surface, via the so-called Markovian projection of the stochastic volatility. We complement the existing results on implied volatility by studying the asymptotic behavior of the local volatility surface generated by a class of rough stochastic volatility models, encompassing the rough Bergomi model. Notably, we observe that the celebrated “1/2 skew rule” linking the short-term at-the-money skew of the implied volatility to the short-term at-the-money skew of the local volatility, a consequence of the celebrated “harmonic mean formula” of [Berestycki et al. (2002). Quantitative Finance, 2, 61–69], is replaced by a new rule: the ratio of the at-the-money implied and local volatility skews tends to the constant $1/(H+3/2)$ (as opposed to the constant 1/2), where H is the regularity index of the underlying instantaneous volatility process.     

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.     

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).     

Regime‐dependent smile‐adjusted delta hedging

We introduce several regime‐dependent smile‐adjusted deltas and compare their efficiency with the smile‐adjusted deltas that are popular with option traders. Using years of daily option prices, out‐of‐sample hedging performance tests for options of all moneyness and maturities and daily, weekly, or fortnightly rebalancing show that even the simplest regime‐dependent smile‐adjustment consistently outperforms implied BSM delta hedging and local volatility and minimum variance smile‐adjustments. Markov‐switching deltas offer the best performance, with delta‐hedging errors often half the size of implied BSM hedging errors. During volatile markets risk reduction from regime‐dependent delta hedging is much greater than during tranquil periods.     

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 .     

Semitractability of optimal stopping problems via a weighted stochastic mesh algorithm

In this paper, we propose a Weighted Stochastic Mesh (WSM) algorithm for approximating the value of discrete‐ and continuous‐time optimal stopping problems. In this context, we consider tractability of such problems via a useful notion of semitractability and the introduction of a tractability index for a particular numerical solution algorithm. It is shown that in the discrete‐time case the WSM algorithm leads to semitractability of the corresponding optimal stopping problem in the sense that its complexity is bounded in order by with being the dimension of the underlying Markov chain. Furthermore, we study the WSM approach in the context of continuous‐time optimal stopping problems and derive the corresponding complexity bounds. Although we cannot prove semitractability in this case, our bounds turn out to be the tightest ones among the complexity bounds known in the literature. We illustrate our theoretical findings by a numerical example.     

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.     

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.     

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.     

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.     

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.     

Cover Image,Volume 30, Issue 2

The cover image is based on the Original Article Semistatic and Sparse Variance‐Optimal Hedging by Di Tella et al., https://doi.org/10.1111/mafi.12235 .

     

Volatility of volatility is (also) rough

Using high-frequency data for major volatility indexes, we compute the volatility of volatility and show that its logarithm follows a fractional Brownian motion with Hurst parameter smaller than 1/2 thereby extending to the volatility asset class the recent findings obtained for the equity index markets. The results confirm that the volatility of volatility is a rough process and it possesses the long memory property. We also show that the correlation between the volatility and the volatility of volatility is positive, consistent with observations in the volatility option market. Lastly, a robustness check using volatility futures confirms the findings.     

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.     

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.     

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.     

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.     

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.