Search

We consider the portfolio optimisation problem where the terminal function is an S-shaped utility applied at the difference between the wealth and a random benchmark process. We develop several numerical methods for solving the problem using deep learning and duality methods. We use deep learning methods to solve the associated Hamilton-Jacobi-Bellman equation for both the primal and dual problems, and the adjoint equation arising from the stochastic maximum principle. We compare the solution of this non-concave problem to that of concavified utility, a random function depending on the benchmark, in both complete and incomplete markets. We give some numerical results for power and log utilities to show the accuracy of the suggested algorithms.

In this paper, we revisit the multistage spectral risk minimization models proposed by Philpott et al.~\cite{PdF13} and Guigues and R\"omisch \cite{GuR12} but with some new focuses. We consider a situation where the decision maker's (DM's) risk preferences may be state-dependent or even inconsistent at some states, and consequently there is not a single deterministic spectral risk measure (SRM) which can be used to represent the DM's preferences at each stage. We adopt the recently introduced average randomized SRM (ARSRM) (in \cite{li2022randomization}) to describe the DM's overall risk preference at each stage. To solve the resulting multistage ARSRM (MARSRM) problem, we apply the well-known stochastic dual dynamic programming (SDDP) method which generates a sequence of lower and upper bounds in an iterative manner. Under some moderate conditions, we prove that the optimal solution can be found in a finite number of iterations. The MARSRM model generalizes the one-stage ARSRM and simplifies the existing multistage state-dependent preference robust model \cite{liu2021multistage}, while also encompassing the mainstream multistage risk-neutral and risk-averse optimization models \cite{GuR12,PdF13}. In the absence of complete information on the probability distribution of the DM's random preferences, we propose to use distributionally robust ARSRM (DR-ARSRM) to describe the DM's preferences at each stage. We detail computational schemes for solving both MARSRM and DR-MARSRM. Finally, we examine the performance of MARSRM and DR-MARSRM by applying them to an asset allocation problem with transaction costs and compare them with standard risk neutral and risk averse multistage linear stochastic programming (MLSP) models., Comment: 33 pages, 4 figures and 3 tables

The paper investigates the consumption-investment problem for an investor with Epstein-Zin utility in an incomplete market. A non-Markovian environment with unbounded parameters is considered, which is more realistic in practical financial scenarios compared to the Markovian setting. The optimal consumption and investment strategies are derived using the martingale optimal principle and quadratic backward stochastic differential equations (BSDEs) whose solutions admit some exponential moment. This integrability property plays a crucial role in establishing a key martingale argument. In addition, the paper also examines the associated dual problem and several models within the specified parameter framework., Comment: 25 pages

This article studies linear-quadratic Stackelberg games between two dominating players (or equivalently, leaders) and a large group of followers, each of whom interacts under a mean field game (MFG) framework. Unlike the conventional major-minor player game, the mean field term herein is endogenously affected by the two leaders simultaneously. These homogeneous followers are non-cooperative, whereas the two leaders can either compete or cooperate with each other, which are respectively formulated as a Nash and a Pareto game. The complete solutions of the leader-follower game can be expressed in terms of the solutions of some non-symmetric Riccati equations. Notably, our analysis suggests that both modes of interactions between leaders has their own merits and neither of them is always more favourable to the community of followers. In our knowledge, a comparative study of the effect of different modes of governance on the society is relatively rare in the existing literature, we here provide its first preliminary quantitative analysis; under a broad class of practically relevant models, we provide sufficient conditions to decide whether cooperation or competition between leaders is more favourable to the followers. Being in common with modern folklore, the relative merits of the two Stackelberg games depend on whether the interests between the two leaders and the followers align among themselves. Representative numerical examples are also supplemented.

In this paper we study a general optimal liquidation problem with a control-dependent stopping time which is the first time the stock holding becomes zero or a fixed terminal time, whichever comes first. We prove a stochastic maximum principle (SMP) which is markedly different in its Hamiltonian condition from that of the standard SMP with fixed terminal time. We present a simple example in which the optimal solution satisfies the SMP in this paper but fails the standard SMP in the literature., Comment: 23 pages, 2 figures

We establish a rigorous duality theory, under No Unbounded Profit with Bounded Risk, for an infinite horizon problem of optimal consumption in the presence of an income stream that can terminate randomly at an exponentially distributed time, independent of the asset prices. We thus close a duality gap encountered by Vellekoop and Davis in a version of this problem in a Black-Scholes market. Many of the classical tenets of duality theory hold, with the notable exception that marginal utility at zero initial wealth is finite. We use as dual variables a class of supermartingale deflators such that deflated wealth plus cumulative deflated consumption in excess of income is a supermartingale. We show that the space of discounted local martingale deflators is dense in our dual domain, so that the dual problem can also be expressed as an infimum over the discounted local martingale deflators. We characterise the optimal wealth process, showing that optimal deflated wealth is a potential decaying to zero, while deflated wealth plus cumulative deflated consumption over income is a uniformly integrable martingale at the optimum. We apply the analysis to the Vellekoop and Davis example and give a numerical solution., Comment: Major revision of earlier paper entitled "Duality and deep learning for optimal consumption with randomly terminating income". arXiv admin note: text overlap with arXiv:2006.04687, arXiv:2009.00972

In this work we study a finite horizon optimal liquidation problem with multiplicative price impact in algorithmic trading, using market orders. We analyze the case when an agent is trading on a market with two financial assets, whose difference of log-prices is modelled with a mean-reverting process. The agent's task is to liquidate an initial position of shares of one of the two financial assets, without having the possibility of trading the other stock. The criterion to be optimized consists in maximising the expected final value of the agent, with a running inventory penalty. The main result of this paper consists in finding a classical solution of the Hamilton-Jacobi-Bellman (HJB) equation associated to this problem, which is proved to not coincide with the value function. However, we find the value function as a solution to the forward-backward stochastic differential equation (FBSDE) associated to the problem. We provide numerical tests showing that the HJB and FBSDE solutions are close to each other and analysing performance of the described model. We also prove a verification theorem and a comparison principle for the viscosity solution to the HJB equation.

This paper proposes two algorithms for solving stochastic control problems with deep learning, with a focus on the utility maximisation problem. The first algorithm solves Markovian problems via the Hamilton Jacobi Bellman (HJB) equation. We solve this highly nonlinear partial differential equation (PDE) with a second order backward stochastic differential equation (2BSDE) formulation. The convex structure of the problem allows us to describe a dual problem that can either verify the original primal approach or bypass some of the complexity. The second algorithm utilises the full power of the duality method to solve non-Markovian problems, which are often beyond the scope of stochastic control solvers in the existing literature. We solve an adjoint BSDE that satisfies the dual optimality conditions. We apply these algorithms to problems with power, log and non-HARA utilities in the Black-Scholes, the Heston stochastic volatility, and path dependent volatility models. Numerical experiments show highly accurate results with low computational cost, supporting our proposed algorithms.

This paper studies an optimal investment and consumption problem with heterogeneous consumption of basic and luxury goods, together with the choice of time for retirement. The utility for luxury goods is not necessarily a concave function. The optimal heterogeneous consumption strategies for a class of non-homothetic utility maximizer are shown to consume only basic goods when the wealth is small, to consume basic goods and make savings when the wealth is intermediate, and to consume almost all in luxury goods when the wealth is large. The optimal retirement policy is shown to be both universal, in the sense that all individuals should retire at the same level of marginal utility that is determined only by income, labor cost, discount factor as well as market parameters, and not universal, in the sense that all individuals can achieve the same marginal utility with different utility and wealth. It is also shown that individuals prefer to retire as time goes by if the marginal labor cost increases faster than that of income. The main tools used in analyzing the problem are from PDE and stochastic control theory including variational inequality and dual transformation. We finally conduct the simulation analysis for the featured model parameters to investigate practical and economic implications by providing their figures.

A speculative agent with Prospect Theory preference chooses the optimal time to purchase and then to sell an indivisible risky asset to maximize the expected utility of the round-trip profit net of transaction costs. The optimization problem is formulated as a sequential optimal stopping problem and we provide a complete characterization of the solution. Depending on the preference and market parameters, the optimal strategy can be "buy and hold", "buy low sell high", "buy high sell higher" or "no trading". Behavioral preference and market friction interact in a subtle way which yields surprising implications on the agent's trading patterns. For example, increasing the market entry fee does not necessarily curb speculative trading, but instead it may induce a higher reference point under which the agent becomes more risk-seeking and in turn is more likely to trade.

In this paper we study a utility maximization problem with both optimal control and optimal stopping in a finite time horizon. The value function can be characterized by a variational equation that involves a free boundary problem of a fully nonlinear partial differential equation. Using the dual control method, we derive the asymptotic properties of the dual value function and the associated dual free boundary for a class of utility functions, including power and non-HARA utilities. We construct a global closed-form approximation to the dual free boundary, which greatly reduces the computational cost. Using the duality relation, we find the approximate formulas for the optimal value function, trading strategy, and exercise boundary for the optimal investment stopping problem. Numerical examples show the approximation is robust, accurate and fast., Comment: 30 pages, 10 figures

We discuss the turnpike property for optimal investment and consumption problems. We find there exists a threshold value that determines the turnpike property for investment policy. The threshold value only depends on the Sharpe ratio, the riskless interest rate and the discount rate. We show that if utilities behave asymptotically like power utilities and satisfy some simple relations with the threshold value, then the turnpike property for investment holds. There is in general no turnpike property for consumption policy. We also provide the rate of convergence and illustrate the main results with examples of power and non-HARA utilities and numerical tests., Comment: 24 pages

The aim of this paper is to study the fast computation of the lower and upper bounds on the value function for utility maximization under the Heston stochastic volatility model with general utility functions. It is well known there is a closed form solution of the HJB equation for power utility due to its homothetic property. It is not possible to get closed form solution for general utilities and there is little literature on the numerical scheme to solve the HJB equation for the Heston model. In this paper we propose an efficient dual control Monte Carlo method for computing tight lower and upper bounds of the value function. We identify a particular form of the dual control which leads to the closed form upper bound for a class of utility functions, including power, non-HARA and Yarri utilities. Finally, we perform some numerical tests to see the efficiency, accuracy, and robustness of the method. The numerical results support strongly our proposed scheme., Comment: 25 pages, 2 figures

In this paper we consider a utility maximization problem with defaultable stocks and looping contagion risk. We assume that the default intensity of one company depends on the stock prices of itself and other companies, and the default of the company induces immediate drops in the stock prices of the surviving companies. We prove that the value function is the unique viscosity solution of the HJB equation. We also perform some numerical tests to compare and analyse the statistical distributions of the terminal wealth of log utility and power utility based on two strategies, one using the full information of intensity process and the other a proxy constant intensity process., Comment: 26 pages, 6 figures. This paper has been presented on 10th Bachelier Conference in Dublin

In this paper, we study a constrained utility maximization problem following the convex duality approach. After formulating the primal and dual problems, we construct the necessary and sufficient conditions for both the primal and dual problems in terms of FBSDEs plus additional conditions. Such formulation then allows us to explicitly characterize the primal optimal control as a function of the adjoint process coming from the dual FBSDEs in a dynamic fashion and vice versa. Moreover, we also find that the optimal primal wealth process coincides with the adjoint process of the dual problem and vice versa. Finally we solve three constrained utility maximization problems, which contrasts the simplicity of the duality approach we propose and the technical complexity of solving the primal problems directly., Comment: 23 pages

In this paper we study a continuous-time stochastic linear quadratic control problem arising from mathematical finance. We model the asset dynamics with random market coefficients and portfolio strategies with convex constraints. Following the convex duality approach, we show that the necessary and sufficient optimality conditions for both the primal and dual problems can be written in terms of processes satisfying a system of FBSDEs together with other conditions. We characterise explicitly the optimal wealth and portfolio processes as functions of adjoint processes from the dual FBSDEs in a dynamic fashion and vice versa. We apply the results to solve quadratic risk minimization problems with cone-constraints and derive the explicit representations of solutions to the extended stochastic Riccati equations for such problems., Comment: 22 pages

We investigate the computational aspects of the basket CDS pricing with counterparty risk under a credit contagion model of multinames. This model enables us to capture the systematic volatility increases in the market triggered by a particular bankruptcy. The drawback of this problem is its analytical complication due to its path-dependent functional, which bears a potential failure in its convergence of numerical approximation under standing assumptions. In this paper we find sufficient conditions for the desired convergence of the functionals associated with a class of path-dependent stochastic differential equations. The main ingredient is to identify the weak convergence of the approximated solution to the underlying path-dependent stochastic differential equation.

In this paper we aim to address two questions faced by a long-term investor with a power-type utility at high levels of wealth: one is whether the turnpike property still holds for a general utility that is not necessarily differentiable or strictly concave, the other is whether the error and the convergence rate of the turnpike property can be estimated. We give positive answers to both questions. To achieve these results, we first show that there is a classical solution to the HJB equation and give a representation of the solution in terms of the dual function of the solution to the dual HJB equation. We demonstrate the usefulness of that representation with some nontrivial examples that would be difficult to solve with the trial and error method. We then combine the dual method and the partial differential equation method to give a direct proof to the turnpike property and to estimate the error and the convergence rate of the optimal policy when the utility function is continuously differentiable and strictly concave. We finally relax the conditions of the utility function and provide some sufficient conditions that guarantee the turnpike property and the convergence rate in terms of both primal and dual utility functions., Comment: 29 pages

In this paper, we study a continuous time structural asset value model for two correlated firms using a two-dimensional Brownian motion. We consider the situation of incomplete information, where the information set available to the market participants includes the default time of each firm and the periodic asset value reports. In this situation, the default time of each firm becomes a totally inaccessible stopping time to the market participants. The original structural model is first transformed to a reduced-form model. Then the conditional distribution of the default time together with the asset value of each name are derived. We prove the existence of the intensity processes of default times and also give the explicit form of the intensity processes. Numerical studies on the intensities of the two correlated names are conducted for some special cases. We also indicate the possible future research extension into three names case by considering a special correlation structure., Comment: 24 pages, 6 figures

In this paper we propose a new robust algorithm to find the optimal static replicating portfolios for general nonlinear payoff functions and give the estimate of the rate of convergence that is absent in the literature. We choose the static replication by minimizing the error bound between the nonlinear payoff function and the linear spline approximation and derive the equidistribution equation for selecting the optimal strike prices. The numerical tests for variance swaps and swaptions and also for the static quadratic replication and the model with counterparty risk show that the proposed algorithm is simple, fast and accurate. The paper has generalized and improved the results of the static replication and approximation in the literature., Comment: 15 pages

In this paper we discuss a credit risk model with a pure jump L\'evy process for the asset value and an unobservable random barrier. The default time is the first time when the asset value falls below the barrier. Using the indistinguishability of the intensity process and the likelihood process, we prove the existence of the intensity process of the default time and find its explicit representation in terms of the distance between the asset value and its running minimal value. We apply the result to find the instantaneous credit spread process and illustrate it with a numerical example., Comment: 15 pages, 2 figures

In the aftermath of the global financial crisis, much attention has been paid to investigating the appropriateness of the current practice of default risk modeling in banking, finance and insurance industries. A recent empirical study by Guo et al.(2008) shows that the time difference between the economic and recorded default dates has a significant impact on recovery rate estimates. Guo et al.(2011) develop a theoretical structural firm asset value model for a firm default process that embeds the distinction of these two default times. To be more consistent with the practice, in this paper, we assume the market participants cannot observe the firm asset value directly and developed a reduced-form model to characterize the economic and recorded default times. We derive the probability distribution of these two default times. The numerical study on the difference between these two shows that our proposed model can both capture the features and fit the empirical data., Comment: arXiv admin note: text overlap with arXiv:1012.0843 by other authors

In this paper, we propose a two-sector Markovian infectious model, which is an extension of Greenwood's model. The central idea of this model is that the causality of defaults of two sectors is in both direction, which enrich dependence dynamics. The Bayesian Information Criterion is adopted to compare the proposed model with the two-sector model in credit literature using the real data. We find that the newly proposed model is statistically better than the model in past literature. We also introduce two measures: CRES and CRVaR to give risk evaluation of our model.

Corporate defaults may be triggered by some major market news or events such as financial crises or collapses of major banks or financial institutions. With a view to develop a more realistic model for credit risk analysis, we introduce a new type of reduced-form intensity-based model that can incorporate the impacts of both observable "trigger" events and economic environment on corporate defaults. The key idea of the model is to augment a Cox process with trigger events. Both single-default and multiple-default cases are considered in this paper. In the former case, a simple expression for the distribution of the default time is obtained. Applications of the proposed model to price defaultable bonds and multi-name Credit Default Swaps (CDSs) are provided.

In this paper we derive an easily computed approximation to European basket call prices for a local volatility jump-diffusion model. We apply the asymptotic expansion method to find the approximate value of the lower bound of European basket call prices. If the local volatility function is time independent then there is a closed-form expression for the approximation. Numerical tests show that the suggested approximation is fast and accurate in comparison with the Monte Carlo and other approximation methods in the literature., Comment: 12 pages

In this paper we continue the study of Bian-Miao-Zheng (2011) and extend the results there to a more general class of utility functions which may be bounded and non-strictly-concave and show that there is a classical solution to the HJB equation with the dual control method. We then apply the results to study the efficient frontier of wealth and conditional VaR (CVaR) problem and the turnpike property problem. For the former we construct explicitly the optimal control and discuss the choice of the optimal threadshold level and illustrate that the wealth and the CVaR are positively correlated. For the latter we give a simple proof to the turnpike property of the optimal policy of long-run investors and generalize the results of Huang-Zariphopoulou (1999)., Comment: 20 pages

In this paper we discuss the optimal liquidation over a finite time horizon until the exit time. The drift and diffusion terms of the asset price are general functions depending on all variables including control and market regime. There is also a local nonlinear transaction cost associated to the liquidation. The model deals with both the permanent impact and the temporary impact in a regime switching framework. The problem can be solved with the dynamic programming principle. The optimal value function is the unique continuous viscosity solution to the HJB equation and can be computed with the finite difference method., Comment: 17 pages

In this paper we prove a weak necessary and sufficient maximum principle for Markovian regime switching stochastic optimal control problems. Instead of insisting on the maximum condition of the Hamiltonian, we show that 0 belongs to the sum of Clarke's generalized gradient of the Hamiltonian and Clarke's normal cone of the control constraint set at the optimal control. Under a joint concavity condition on the Hamiltonian and a convexity condition on the terminal objective function, the necessary condition becomes sufficient. We give four examples to demonstrate the weak stochastic maximum principle., Comment: 30 pages

In this paper we propose a simple and efficient method to compute the ordered default time distributions in both the homogeneous case and the two-group heterogeneous case under the interacting intensity default contagion model. We give the analytical expressions for the ordered default time distributions with recursive formulas for the coefficients, which makes the calculation fast and efficient in finding rates of basket CDSs. In the homogeneous case, we explore the ordered default time in limiting case and further include the exponential decay and the multistate stochastic intensity process. The numerical study indicates that, in the valuation of the swap rates and their sensitivities with respect to underlying parameters, our proposed model outperforms the Monte Carlo method.

We maximize the expected utility of terminal wealth in an incomplete market where there are cone constraints on the investor's portfolio process and the utility function is not assumed to be strictly concave or differentiable. We establish the existence of the optimal solutions to the primal and dual problems and their dual relationship. We simplify the present proofs in this area and extend the existing duality theory to the constrained nonsmooth setting., Comment: 20 pages

In this paper we propose a copula contagion mixture model for correlated default times. The model includes the well known factor, copula, and contagion models as its special cases. The key advantage of such a model is that we can study the interaction of different models and their pricing impact. Specifically, we model the marginal default times to follow some contagion intensity processes coupled with copula dependence structure. We apply the total hazard construction method to generate ordered default times and numerically compare the pricing impact of different models on basket CDSs and CDOs in the presence of exponential decay and counterparty risk., Comment: 12 pages

In this paper we prove that there exists a smooth classical solution to the HJB equation for a large class of constrained problems with utility functions that are not necessarily differentiable or strictly concave. The value function is smooth if admissible controls satisfy an integrability condition or if it is continuous on the closure of its domain. The key idea is to work on the dual control problem and the dual HJB equation. We construct a smooth, strictly convex solution to the dual HJB equation and show that its conjugate function is a smooth, strictly concave solution to the primal HJB equation satisfying the terminal and boundary conditions., Comment: 18 pages

In this paper we discuss the basket options valuation for a jump-diffusion model. The underlying asset prices follow some correlated local volatility diffusion processes with systematic jumps. We derive a forward partial integral differential equation (PIDE) for general stochastic processes and use the asymptotic expansion method to approximate the conditional expectation of the stochastic variance associated with the basket value process. The numerical tests show that the suggested method is fast and accurate in comparison with the Monte Carlo and other methods in most cases., Comment: 16 pages, 4 tables