Optimal Trade Execution with Instantaneous Price Impact and Stochastic Resilience

Similar documents
On continuous time contract theory

Nonlinear representation, backward SDEs, and application to the Principal-Agent problem

Singular Perturbations of Stochastic Control Problems with Unbounded Fast Variables

A class of globally solvable systems of BSDE and applications

A new approach for investment performance measurement. 3rd WCMF, Santa Barbara November 2009

Minimal Supersolutions of Backward Stochastic Differential Equations and Robust Hedging

Order book modeling and market making under uncertainty.

Numerical scheme for quadratic BSDEs and Dynamic Risk Measures

Order book resilience, price manipulation, and the positive portfolio problem

On the Multi-Dimensional Controller and Stopper Games

Thomas Knispel Leibniz Universität Hannover

Backward Stochastic Differential Equations with Infinite Time Horizon

Stochastic Calculus. Kevin Sinclair. August 2, 2016

Theoretical Tutorial Session 2

Economics 2010c: Lectures 9-10 Bellman Equation in Continuous Time

Uniformly Uniformly-ergodic Markov chains and BSDEs

An Uncertain Control Model with Application to. Production-Inventory System

HJB equations. Seminar in Stochastic Modelling in Economics and Finance January 10, 2011

Quadratic BSDE systems and applications

Random G -Expectations

Proving the Regularity of the Minimal Probability of Ruin via a Game of Stopping and Control

Mathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio ( )

Multi-dimensional Stochastic Singular Control Via Dynkin Game and Dirichlet Form

Shadow prices and well-posedness in the problem of optimal investment and consumption with transaction costs

Harnack Inequalities and Applications for Stochastic Equations

Prof. Erhan Bayraktar (University of Michigan)

Econ Lecture 3. Outline. 1. Metric Spaces and Normed Spaces 2. Convergence of Sequences in Metric Spaces 3. Sequences in R and R n

On Stopping Times and Impulse Control with Constraint

Markov Chain BSDEs and risk averse networks

B8.3 Mathematical Models for Financial Derivatives. Hilary Term Solution Sheet 2

Stochastic optimal control with rough paths

Homogenization with stochastic differential equations

An Introduction to Moral Hazard in Continuous Time

Control of McKean-Vlasov Dynamics versus Mean Field Games

Classical and Restricted Impulse Control for the Exchange Rate under a Stochastic Trend Model

Convergence of Particle Filtering Method for Nonlinear Estimation of Vortex Dynamics

Peter Kratz May 17, 2018

ON THE POLICY IMPROVEMENT ALGORITHM IN CONTINUOUS TIME

Non-Zero-Sum Stochastic Differential Games of Controls and St

Lecture 3. Econ August 12

Reflected Brownian Motion

1 Lyapunov theory of stability

Controlled Diffusions and Hamilton-Jacobi Bellman Equations

Dynamic Risk Measures and Nonlinear Expectations with Markov Chain noise

Conditional and Dynamic Preferences

56 4 Integration against rough paths

Econ Lecture 14. Outline

Uncertainty in Kalman Bucy Filtering

1. Nonlinear Equations. This lecture note excerpted parts from Michael Heath and Max Gunzburger. f(x) = 0

Weak solutions of mean-field stochastic differential equations

Variance Reduction Techniques for Monte Carlo Simulations with Stochastic Volatility Models

Optimal trade execution and absence of price manipulations in limit order book models

Interpolation and function representation with grids in stochastic control

Universität des Saarlandes. Fachrichtung 6.1 Mathematik

Continuous Time Finance

Introduction to Algorithmic Trading Strategies Lecture 4

Stabilization of persistently excited linear systems

Robust Mean-Variance Hedging via G-Expectation

Perturbative Approaches for Robust Intertemporal Optimal Portfolio Selection

The priority option: the value of being a leader in complete and incomplete markets.

SDE Coefficients. March 4, 2008

Deterministic Dynamic Programming

Paracontrolled KPZ equation

arxiv: v1 [math.ap] 20 Dec 2018

Brownian Motion. 1 Definition Brownian Motion Wiener measure... 3

Large deviations and averaging for systems of slow fast stochastic reaction diffusion equations.

Consumption. Consider a consumer with utility. v(c τ )e ρ(τ t) dτ.

Worst Case Portfolio Optimization and HJB-Systems

Vector Calculus. A primer

The first order quasi-linear PDEs

Optimal portfolio strategies under partial information with expert opinions

2tdt 1 y = t2 + C y = which implies C = 1 and the solution is y = 1

Set-Valued Risk Measures and Bellman s Principle

Strong uniqueness for stochastic evolution equations with possibly unbounded measurable drift term

General Intensity Shapes in Optimal Liquidation

Bessel-like SPDEs. Lorenzo Zambotti, Sorbonne Université (joint work with Henri Elad-Altman) 15th May 2018, Luminy

Stochastic Differential Equations

1. Stochastic Processes and filtrations

Question 1. The correct answers are: (a) (2) (b) (1) (c) (2) (d) (3) (e) (2) (f) (1) (g) (2) (h) (1)

2 A Model, Harmonic Map, Problem

A Stochastic Paradox for Reflected Brownian Motion?

Lecture 4: Ito s Stochastic Calculus and SDE. Seung Yeal Ha Dept of Mathematical Sciences Seoul National University

1. Stochastic Process

Short-time expansions for close-to-the-money options under a Lévy jump model with stochastic volatility

Mean-Field optimization problems and non-anticipative optimal transport. Beatrice Acciaio

On Chern-Simons-Schrödinger equations including a vortex point

The Azéma-Yor Embedding in Non-Singular Diffusions

Stability of Stochastic Differential Equations

A note on stabilization with saturating feedback. Thomas I. Seidman. Houshi Li

Lecture No 1 Introduction to Diffusion equations The heat equat

Viscosity Solutions of Path-dependent Integro-Differential Equations

Optimal Saving under Poisson Uncertainty: Corrigendum

Production and Relative Consumption

Bessel Functions Michael Taylor. Lecture Notes for Math 524

Functional Limit theorems for the quadratic variation of a continuous time random walk and for certain stochastic integrals

z x = f x (x, y, a, b), z y = f y (x, y, a, b). F(x, y, z, z x, z y ) = 0. This is a PDE for the unknown function of two independent variables.

OBSERVATION OPERATORS FOR EVOLUTIONARY INTEGRAL EQUATIONS

Dynamical systems with Gaussian and Levy noise: analytical and stochastic approaches

Some Aspects of Universal Portfolio

Backward martingale representation and endogenous completeness in finance

Transcription:

Optimal Trade Execution with Instantaneous Price Impact and Stochastic Resilience Ulrich Horst 1 Humboldt-Universität zu Berlin Department of Mathematics and School of Business and Economics Vienna, Nov. 18th, 2016 1 Based on joint work with Paulwin Graewe

Aims and Goals portfolio liquidation model with instantaneous price impact permanent price impact stochastic resilience only absolutely continuous trading strategies what s new? value functions can be described by a coupled BSDE system one component has a singular terminal condition solved using asymptotic expansion at the terminal time

Portfolio liquidation models models with instantaneous price impact LQ cost function absolutely continuous strategies value functions described by one-dimensional PDEs, B(S)PDEs models with permanent price impact and resilience LQ cost function absolutely continuous and block trades absolutely continuous trades: permanent impact block trades: instantaneous and permanent impact optimal strategy characterised in terms of BSDE (systems)

A blend... only absolutely continuous strategies (value function) persistent price impact with only a.c. strategies: trading rate adds drift to a benchmark price impact of past trades on current prices decreases instantaneous impact with only a.c. strategies: no instantaneous impact: trade infinitely fast weak instantaneous impact: trade fast initially large instantaneous impact: no permanent impact

Outline the control problem the BSDE system for the value function solving the BSDE system the verification argument

The control problem

The control problem For any initial state (t, x, y) [0, T ) R R the value function is [ T V t (x, y) = ess inf E ξ A (t,x) t 1 2 ηξ2 s + ξ s Y s + 1 2 λ sx 2 s ds ] F t where { dxs = ξ s ds, t s T ; X t = x; dy s = { ρ s Y s + γξ s } ds, t s T ; Y t = y. The process ξ L 2 F (t, T ; R) is the control and η, γ R + ; ρ, λ L F (0, T ; R +). Remark ρ 0 (y = 0): only temporary impact.

The control problem A control is called admissible if the terminal state constraint is satisfied a.s. and X T = 0 ξ L 2 F (t, T ; R). The set of all admissible controls is denoted A (t, x). Remark Since ξ L 2 F (t, T ; R) boundedness of the coefficients guarantees X, Y L 2 F (Ω; C([t, T ]; R)).

The control problem Recall that [ T V t (x, y) = ess inf E ξ A (t,x) t 1 2 ηξ2 s + ξ s Y s + 1 2 λ sx 2 s ds ] F t X t 1.0 0.8 0.6 Out[167]= 0.4 0.2 η = 10 η = 0.5 η = 0.01 Obizhaeva & Wang Almgren & Chriss 0.0 0.2 0.4 0.6 0.8 1.0 t

The HJB equation The stochastic HJB equation to our problem is dv t (x, y) = inf ξ R {... xv t (x, y)... y V t (x, y)...} dt Z t (x, y) dw t. Remark The HJB equation is given by the cost function and the state dynamics; the terminal condition by the liquidation constraint.

The HJB equation The LQ structure of the control problem suggest the ansatz V t (x, y) = 1 2 A tx 2 + B t xy + 1 2 C ty 2 Z t (x, y) = 1 2 Z t A x 2 + Zt B xy + 1 2 Z t C y 2 for the solution (V (x, y), Z(x, y)) to the HJB equation, where da t = { λ t η 1 (A t γb t ) 2} dt Zt A dw t db t = { ρ t B t + η 1 (γc t B t + 1)(A t γb t ) } dt Zt B dc t = { 2ρ t C t η 1 (γc t B t + 1) 2} dt Zt C dw t. dw t What is the terminal condition of this BSDE system?

The terminal condition we expect the trading rate ξ to tend to infinity for any non-trivial initial position as t T. we expect the resulting trading cost to dominate any resilience effect. we expect that V t (x, y) V ρ=0 t (x, y) as t T where V ρ=0 is the value function corresponding to ρ 0.

Lemma (The case ρ 0 (Graewe, H and Qiu (2015))) Vt ρ=0 (x, y) = 1 2 (Ãt + γ)x 2 + xy, where { d à t = λ t η 1 à 2 t dt Z t dw t à t in L as t T. Remark From this lemma we see why ρ 0 corresponds to a model with only temporary impact.

The terminal condition Since and in view of the ansatz Vt ρ=0 (x, y) = 1 2 (Ã t + γ)x 2 + xy, V t (x, y) = 1 2 A tx 2 + B t xy + 1 2 C ty 2 we expect the coefficients of the linear-quadratic ansatz to satisfy (A t, B t, C t ) (, 1, 0) in L as t T. Remark Our approach uses the precise rate of A t in L.

Theorem (Existence and uniqueness of solutions) The above BSDE system imposed with the above singular terminal condition admits at least one solution ((A, B, C), (Z A, Z B, Z C )) L F (Ω; C([0, T ]; R 3 )) L 2 F (0, T ; R 3 m ). Suppose a solution exists, then the value function is of the linear quadratic form and the optimal strategy is ξ t (x, y) = η 1 (A t γb t )x η 1 (γc t B t + 1)y. In particular, the BSDE systems admits at most one solution.

Remark Recall that We show that ξ t (x, y) = η 1 (A t γb t )x η 1 (γc t B t + 1)y. A t γb t > 0; γc t B t + 1 0 but we have no result on the sign of ξ along (X, Y ). We can not rule our price triggered round trips.

Solving the BSDE system

The BSDE system We need to solve a fully coupled BSDE system with singular terminal condition to which multi-dimensional comparison results don t apply The idea is to find an equivalent BSDE system with regular terminal condition but singular driver that can be solved in a suitable space: as t T : A t = η T t + H t (T t) 2, B t = 1 + G t T t, C t = P t where H t, B t, C t = O((T t) 2 ).

Asymptotics: some intuition Recall that da t = { λ t η 1 (A t γb t ) 2} dt Zt A dw t db t = { ρ t B t + η 1 (γc t B t + 1)(A t γb t ) } dt Zt B dc t = { 2ρ t C t η 1 (γc t B t + 1) 2} dt Zt C dw t. dw t B T = 1 A behaves as in the one-dimensional case rate of convergence of B is slower than that of C A t 1 T t B t 1 = O(T t)

A priori estimates Let ((A, B, C), (Z A, Z B, Z C )) denote any solution to our BSDE system that satisfies (A t, B t, C t ) (, 1, 0) in L as t T. It will be necessary to also consider the processes D := η 1 (A γb) and E := η 1 (γc B + 1) that appear in the characterization of the candidate strategy. Remark We need the full picture: value function and strategy. The BSDEs for D, E and the system for (B, D) allow for comparison.

A priori estimates Using multi-dimensional comparison we can show that a.s. A, D 0, B, γc, ηe [0, 1] Since ξt (x, y) = D t x E t y this does not guarantee that ξ 0.

A priori estimates using a comparison principle for quasi-monotone BSDE systems we obtain a priori estimates for B, D and E from these we conclude that the following asymptotic behaviors hold in L as t T : (T t)a t = η + O(T t), B t = 1 + O(T t), C t = O((T t) 3 ).

Existence The asymptotic behavior suggests the following ansatz: A t = η T t + H t (T t) 2, H t = O((T t) 2 ) B t = 1 + G t T t, G t = O((T t) 2 ) C t = P t, P t = O((T t) 2 ) where H, G and P are of the form { ( ) p } Yt dy t =... T t + coupling dt + Z t dw t.

Existence We have a BSDE system with singular driver, dy t = f (ω, t, Y t )dt Z t dw t We need to look for a solution in the right space, namely H = {Y L F (Ω; C([T δ, T ]; R3 )) : Y H < + } endowed with the norm Y H = (T ) 2 Y L F (Ω;C([T δ,t ];R 3 )). This space is complete and the driver f is locally Lipschitz: f (, Y ) f (, X ) H L Y X H Y, X B H (R).

Lemma There exists a short-time solution (Y, Z) H 2 L 2 F (T δ, T ; R3 m ) to the above BSDE with singular driver. In particular, there exists a small time solution to the BSDE system for (A, B, C) with singular terminal value. The small time solution for (A, B, C) can be extended to a global solution on [0, T ). Remark The extension uses that f is independent of Z, analogously to PDEs, when the inhomogeneity is independent of the gradient.

Verification

Verification we can apply Itô on [0, s] for all s < T to get: V t (x, y) E t [V s (X s, Y s )] [ s + E t t { 1 2 ηξ2 r + ξ r Y r + 1 2 λ r X 2 r } dr with equality for the candidate strategy ξ ].

Lemma the candidate is admissible and ξ, Y L F (Ω; C([t, T ]; R)). for every ξ A (t, x) it holds that E t [A s X 2 s + B s X s Y s + C s Y 2 s ] s T 0.

Verification By the Itô-Kunita formula we have (after stopping): [ s ] V t (x, y) E t [V s (X s, Y s )] + E t { 1 2 ηξ2 r + ξ r Y r + 1 2 λ r Xr 2 } dr. Since X, Y L 2 F (Ω; C([t, T ]; R)), and ξ L2 F (t, T ; R), we get [ T ] 1 V t (x, y) E t 2 ηξ2 r + ξ r Y r + 1 2 λ r Xr 2 dr. t t

Conclusion liquidation model with instantaneous and permanent impact main results: value function is characterised by a BSDE system BSDE system satisfies a singular terminal condition existence via an asymptotic expansion at the terminal time main limiting factors: a priori estimates comparison result some open problems: sign of the optimal trading strategy introduce dark pools random order book hight (liquidity) convergence as η 0...

Many Thanks!