Source author record

Matthew Lorig

Matthew Lorig appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

26works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

26 published item(s)

preprint2026arXiv

Optimal Liquidation of Perpetual Contracts

An agent holds a position in a perpetual contract with payoff function $ψ$ and attempts to liquidate the position while managing transaction costs, inventory risk, and funding rate payments. By solving the agent's stochastic control problem we obtain a closed-form expression for the optimal trading strategy when the payoff function is given by $ψ(s) = s$. When the payoff function is non-linear we provide approximations to the optimal strategy which apply when the funding rate parameter is small or when the length of the trading interval is small. We further prove that when $ψ$ is non-linear, the short time approximation can be written in terms of the closed-form trading strategy corresponding to the case of the identity payoff function.

preprint2022arXiv

A primer on perpetuals

We consider a continuous-time financial market with no arbitrage and no transactions costs. In this setting, we introduce two types of perpetual contracts, one in which the payoff to the long side is a fixed function of the underlyers and the long side pays a funding rate to the short side, the other in which the payoff to the long side is a fixed function of the underlyers times a discount factor that changes over time but no funding payments are required. Assuming asset prices are continuous and strictly positive, we derive model-free expressions for the funding rate and discount rate of these perpetual contracts as well as replication strategies for the short side. When asset prices can jump, we derive expressions for the funding and discount rates, which are semi-robust in the sense that they do not depend on the dynamics of the volatility process of the underlying risky assets, but do depend on the intensity of jumps under the market's pricing measure. When asset prices can jump and the volatility process is independent of the underlying risky assets, we derive an explicit replication strategy for the short side of a perpetual contract. Throughout the paper, we illustrate through examples how specific perpetual contracts relate to traditional financial instruments such as variance swaps and leveraged exchange traded funds.

preprint2022arXiv

Optimal times to buy and sell a home

We consider a financial market in which the risk-free rate of interest is modeled as a Markov diffusion. We suppose that home prices are set by a representative home-buyer, who can afford to pay only a fixed cash-flow per unit time for housing. The cash-flow is a fraction of the representative home-buyer's salary, which grows at a rate that is proportional to the risk-free rate of interest. As a result, in the long-run, higher interest rates lead to faster growth of home prices. The representative home-buyer finances the purchase of a home by taking out a mortgage. The mortgage rate paid by the home-buyer is fixed at the time of purchase and equal to the risk-free rate of interest plus a positive constant. As the home-buyer can only afford to pay a fixed cash-flow per unit time, a higher mortgage rate limits the size of the loan the home-buyer can take out. As a result, the short-term effect of higher interest rates is to lower the value of homes. In this setting, we consider an investor who wishes to buy and then sell a home in order to maximize his discounted expected profit. This leads to a nested optimal stopping problem. We use a nonnegative concave majorant approach to derive the investor's optimal buying and selling strategies. Additionally, we provide a detailed analytic and numerical study of the case in which the risk-free rate of interest is modeled by a Cox-Ingersoll-Ross (CIR) process. We also examine, in the case of CIR interest rates, the expected time that the investor waits before buying and then selling a home when following the optimal strategies.

preprint2022arXiv

Robust replication of barrier-style claims on price and volatility

We show how to price and replicate a variety of barrier-style claims written on the $\log$ price $X$ and quadratic variation $\langle X \rangle$ of a risky asset. Our framework assumes no arbitrage, frictionless markets and zero interest rates. We model the risky asset as a strictly positive continuous semimartingale with an independent volatility process. The volatility process may exhibit jumps and may be non-Markovian. As hedging instruments, we use only the underlying risky asset, zero-coupon bonds, and European calls and puts with the same maturity as the barrier-style claim. We consider knock-in, knock-out and rebate claims in single and double barrier varieties.

preprint2021arXiv

Optimal Bookmaking

We introduce a general framework for continuous-time betting markets, in which a bookmaker can dynamically control the prices of bets on outcomes of random events. In turn, the prices set by the bookmaker affect the rate or intensity of bets placed by gamblers. The bookmaker seeks a price process that maximizes his expected (utility of) terminal wealth. We obtain explicit solutions or characterizations to the bookmaker's optimal bookmaking problem in various interesting models.

preprint2020arXiv

Bond indifference prices and indifference yield curves

In a market with stochastic interest rates, we consider an investor who can either (i) invest all if his money in a savings account or (ii) purchase zero-coupon bonds and invest the remainder of his wealth in a savings account. The indifference price of the bond is the price for which the investor could achieve the same expected utility under both scenarios. In an affine term structure setting, under the assumption that an investor has a utility function in either exponential or power form, we show that the indifference price of a zero-coupon bond is the root of an integral expression. As an example, we compute bond indifference prices and the corresponding indifference yield curves in the Vasicek setting and interpret the results.

preprint2015arXiv

Indifference prices and implied volatilities

We consider a general local-stochastic volatility model and an investor with exponential utility. For a European-style contingent claim, whose payoff may depend on either a traded or non-traded asset, we derive an explicit approximation for both the buyer's and seller's indifference price. For European calls on a traded asset, we translate indifference prices into an explicit approximation of the buyer's and seller's implied volatility surface. For European claims on a non-traded asset, we establish rigorous error bounds for the indifference price approximation. Finally, we implement our indifference price and implied volatility approximations in two examples.

preprint2015arXiv

Leveraged {ETF} implied volatilities from {ETF} dynamics

The growth of the exhange-traded fund (ETF) industry has given rise to the trading of options written on ETFs and their leveraged counterparts {(LETFs)}. We study the relationship between the ETF and LETF implied volatility surfaces when the underlying ETF is modeled by a general class of local-stochastic volatility models. A closed-form approximation for prices is derived for European-style options whose payoff depends on the terminal value of the ETF and/or LETF. Rigorous error bounds for this pricing approximation are established. A closed-form approximation for implied volatilities is also derived. We also discuss a scaling procedure for comparing implied volatilities across leverage ratios. The implied volatility expansions and scalings are tested in three well-known settings: CEV, Heston and SABR.

preprint2015arXiv

Optimal Static Quadratic Hedging

We propose a flexible framework for hedging a contingent claim by holding static positions in vanilla European calls, puts, bonds, and forwards. A model-free expression is derived for the optimal static hedging strategy that minimizes the expected squared hedging error subject to a cost constraint. The optimal hedge involves computing a number of expectations that reflect the dependence among the contingent claim and the hedging assets. We provide a general method for approximating these expectations analytically in a general Markov diffusion market. To illustrate the versatility of our approach, we present several numerical examples, including hedging path-dependent options and options written on a correlated asset.

preprint2015arXiv

Portfolio Optimization under Local-Stochastic Volatility: Coefficient Taylor Series Approximations & Implied Sharpe Ratio

We study the finite horizon Merton portfolio optimization problem in a general local-stochastic volatility setting. Using model coefficient expansion techniques, we derive approximations for the both the value function and the optimal investment strategy. We also analyze the `implied Sharpe ratio' and derive a series approximation for this quantity. The zeroth-order approximation of the value function and optimal investment strategy correspond to those obtained by Merton (1969) when the risky asset follows a geometric Brownian motion. The first-order correction of the value function can, for general utility functions, be expressed as a differential operator acting on the zeroth-order term. For power utility functions, higher order terms can also be computed as a differential operator acting on the zeroth-order term. We give a rigorous accuracy bound for the higher order approximations in this case in pure stochastic volatility models. A number of examples are provided in order to demonstrate numerically the accuracy of our approximations.

preprint2015arXiv

Second Order Multiscale Stochastic Volatility Asymptotics: Stochastic Terminal Layer Analysis & Calibration

Multiscale stochastic volatility models have been developed as an efficient way to capture the principle effects on derivative pricing and portfolio optimization of randomly varying volatility. The recent book Fouque, Papanicolaou, Sircar and Sølna (2011, CUP) analyzes models in which the volatility of the underlying is driven by two diffusions -- one fast mean-reverting and one slow-varying, and provides a first order approximation for European option prices and for the implied volatility surface, which is calibrated to market data. Here, we present the full second order asymptotics, which are considerably more complicated due to a terminal layer near the option expiration time. We find that, to second order, the implied volatility approximation depends quadratically on log-moneyness, capturing the convexity of the implied volatility curve seen in data. We introduce a new probabilistic approach to the terminal layer analysis needed for the derivation of the second order singular perturbation term, and calibrate to S&P 500 options data.

preprint2014arXiv

Asymptotics for $d$-dimensional Lévy-type processes

We consider a general d-dimensional Levy-type process with killing. Combining the classical Dyson series approach with a novel polynomial expansion of the generator A(t) of the Levy-type process, we derive a family of asymptotic approximations for transition densities and European-style options prices. Examples of stochastic volatility models with jumps are provided in order to illustrate the numerical accuracy of our approach. The methods described in this paper extend the results from Corielli et al. (2010), Pagliarani and Pascucci (2013) and Lorig et al. (2013a) for Markov diffusions to Markov processes with jumps.

preprint2014arXiv

Explicit implied volatilities for multifactor local-stochastic volatility models

We consider an asset whose risk-neutral dynamics are described by a general class of local-stochastic volatility models and derive a family of asymptotic expansions for European-style option prices and implied volatilities. Our implied volatility expansions are explicit; they do not require any special functions nor do they require numerical integration. To illustrate the accuracy and versatility of our method, we implement it under five different model dynamics: CEV local volatility, quadratic local volatility, Heston stochastic volatility, $3/2$ stochastic volatility, and SABR local-stochastic volatility.

preprint2014arXiv

From characteristic functions to implied volatility expansions

For any strictly positive martingale $S = \exp(X)$ for which $X$ has a characteristic function, we provide an expansion for the implied volatility. This expansion is explicit in the sense that it involves no integrals, but only polynomials in the log strike. We illustrate the versatility of our expansion by computing the approximate implied volatility smile in three well-known martingale models: one finite activity exponential Lévy model (Merton), one infinite activity exponential Lévy model (Variance Gamma), and one stochastic volatility model (Heston). Finally, we illustrate how our expansion can be used to perform a model-free calibration of the empirically observed implied volatility surface.

preprint2014arXiv

Pricing approximations and error estimates for local Lévy-type models with default

We find approximate solutions of partial integro-differential equations, which arise in financial models when defaultable assets are described by general scalar Lévy-type stochastic processes. We derive rigorous error bounds for the approximate solutions. We also provide numerical examples illustrating the usefulness and versatility of our methods in a variety of financial settings.

preprint2013arXiv

A family of density expansions for Lévy-type processes

We consider a defaultable asset whose risk-neutral pricing dynamics are described by an exponential Levy-type martingale subject to default. This class of models allows for local volatility, local default intensity, and a locally dependent Levy measure. Generalizing and extending the novel adjoint expansion technique of Pagliarani, Pascucci, and Riga (2013), we derive a family of asymptotic expansions for the transition density of the underlying as well as for European-style option prices and defaultable bond prices. For the density expansion, we also provide error bounds for the truncated asymptotic series. Our method is numerically efficient; approximate transition densities and European option prices are computed via Fourier transforms; approximate bond prices are computed as finite series. Additionally, as in Pagliarani et al. (2013), for models with Gaussian-type jumps, approximate option prices can be computed in closed form. Sample Mathematica code is provided.

preprint2013arXiv

A Taylor series approach to pricing and implied vol for LSV models

Using classical Taylor series techniques, we develop a unified approach to pricing and implied volatility for European-style options in a general local-stochastic volatility setting. Our price approximations require only a normal CDF and our implied volatility approximations are fully explicit (ie, they require no special functions, no infinite series and no numerical integration). As such, approximate prices can be computed as efficiently as Black-Scholes prices, and approximate implied volatilities can be computed nearly instantaneously.

preprint2013arXiv

Exponential Lévy-type models with stochastic volatility and stochastic jump-intensity

We consider the problem of valuing a European option written on an asset whose dynamics are described by an exponential Lévy-type model. In our framework, both the volatility and jump-intensity are allowed to vary stochastically in time through common driving factors -- one fast-varying and one slow-varying. Using Fourier analysis we derive an explicit formula for the approximate price of any European-style derivative whose payoff has a generalized Fourier transform; in particular, this includes European calls and puts. From a theoretical perspective, our results extend the class of multiscale stochastic volatility models of \citet*{fpss} to models of the exponential Lévy type. From a financial perspective, the inclusion of jumps and stochastic volatility allow us to capture the term-structure of implied volatility. To illustrate the flexibility of our modeling framework we extend five exponential Lévy processes to include stochastic volatility and jump-intensity. For each of the extended models, using a single fast-varying factor of volatility and jump-intensity, we perform a calibration to the S&P500 implied volatility surface. Our results show decisively that the extended framework provides a significantly better fit to implied volatility than both the traditional exponential Lévy models and the fast mean-reverting stochastic volatility models of \citet{fpss}.

preprint2013arXiv

The Smile of certain Lévy-type Models

We consider a class of assets whose risk-neutral pricing dynamics are described by an exponential Lévy-type process subject to default. The class of processes we consider features locally-dependent drift, diffusion and default-intensity as well as a locally-dependent Lévy measure. Using techniques from regular perturbation theory and Fourier analysis, we derive a series expansion for the price of a European-style option. We also provide precise conditions under which this series expansion converges to the exact price. Additionally, for a certain subclass of assets in our modeling framework, we derive an expansion for the implied volatility induced by our option pricing formula. The implied volatility expansion is exact within its radius of convergence. As an example of our framework, we propose a class of CEV-like Lévy-type models. Within this class, approximate option prices can be computed by a single Fourier integral and approximate implied volatilities are explicit (i.e., no integration is required). Furthermore, the class of CEV-like Lévy-type models is shown to provide a tight fit to the implied volatility surface of S{&}P500 index options.

preprint2013arXiv

Variance Swaps on Defaultable Assets and Market Implied Time-Changes

We compute the value of a variance swap when the underlying is modeled as a Markov process time changed by a Lévy subordinator. In this framework, the underlying may exhibit jumps with a state-dependent Lévy measure, local stochastic volatility and have a local stochastic default intensity. Moreover, the Lévy subordinator that drives the underlying can be obtained directly by observing European call/put prices. To illustrate our general framework, we provide an explicit formula for the value of a variance swap when the underlying is modeled as (i) a Lévy subordinated geometric Brownian motion with default and (ii) a Lévy subordinated Jump-to-default CEV process (see \citet{carr-linetsky-1}). {In the latter example, we extend} the results of \cite{mendoza-carr-linetsky-1}, by allowing for joint valuation of credit and equity derivatives as well as variance swaps.

preprint2012arXiv

A Fast Mean-Reverting Correction to Heston's Stochastic Volatility Model

We propose a multi-scale stochastic volatility model in which a fast mean-reverting factor of volatility is built on top of the Heston stochastic volatility model. A singular pertubative expansion is then used to obtain an approximation for European option prices. The resulting pricing formulas are semi-analytic, in the sense that they can be expressed as integrals. Difficulties associated with the numerical evaluation of these integrals are discussed, and techniques for avoiding these difficulties are provided. Overall, it is shown that computational complexity for our model is comparable to the case of a pure Heston model, but our correction brings significant flexibility in terms of fitting to the implied volatility surface. This is illustrated numerically and with option data.

preprint2012arXiv

Pricing Derivatives on Multiscale Diffusions: an Eigenfunction Expansion Approach

Using tools from spectral analysis, singular and regular perturbation theory, we develop a systematic method for analytically computing the approximate price of a derivative-asset. The payoff of the derivative-asset may be path-dependent. Additionally, the process underlying the derivative may exhibit killing (i.e. jump to default) as well as combined local/nonlocal stochastic volatility. The nonlocal component of volatility is multiscale, in the sense that it is driven by one fast-varying and one slow-varying factor. The flexibility of our modeling framework is contrasted by the simplicity of our method. We reduce the derivative pricing problem to that of solving a single eigenvalue equation. Once the eigenvalue equation is solved, the approximate price of a derivative can be calculated formulaically. To illustrate our method, we calculate the approximate price of three derivative-assets: a vanilla option on a defaultable stock, a path-dependent option on a non-defaultable stock, and a bond in a short-rate model.

preprint2012arXiv

Spectral Decomposition of Option Prices in Fast Mean-Reverting Stochastic Volatility Models

Using spectral decomposition techniques and singular perturbation theory, we develop a systematic method to approximate the prices of a variety of options in a fast mean-reverting stochastic volatility setting. Four examples are provided in order to demonstrate the versatility of our method. These include: European options, up-and-out options, double-barrier knock-out options, and options which pay a rebate upon hitting a boundary. For European options, our method is shown to produce option price approximations which are equivalent to those developed in [5]. [5] Jean-Pierre Fouque, George Papanicolaou, and Sircar Ronnie. Derivatives in Financial Markets with Stochas- tic Volatility. Cambridge University Press, 2000.

preprint2012arXiv

Time-Changed Fast Mean-Reverting Stochastic Volatility Models

We introduce a class of randomly time-changed fast mean-reverting stochastic volatility models and, using spectral theory and singular perturbation techniques, we derive an approximation for the prices of European options in this setting. Three examples of random time-changes are provided and the implied volatility surfaces induced by these time-changes are examined as a function of the model parameters. Three key features of our framework are that we are able to incorporate jumps into the price process of the underlying asset, allow for the leverage effect, and accommodate multiple factors of volatility, which operate on different time-scales.