Nonzero-sum stochastic differential games with impulse controls offer a realistic and far-reaching modelling framework for applications within finance, energy markets, and other areas, but the difficulty in solving such problems has hindered their proliferation. Semi-analytical approaches make strong assumptions pertaining to very particular cases. To the author's best knowledge, the only numerical method in the literature is the heuristic one we put forward to solve an underlying system of quasi-variational inequalities. Focusing on symmetric games, this paper presents a simpler, more precise and efficient fixed-point policy-iteration-type algorithm which removes the strong dependence on the initial guess and the relaxation scheme of the previous method. A rigorous convergence analysis is undertaken with natural assumptions on the players strategies, which admit graph-theoretic interpretations in the context of weakly chained diagonally dominant matrices. A novel provably convergent single-player impulse control solver is also provided. The main algorithm is used to compute with high precision equilibrium payoffs and Nash equilibria of otherwise very challenging problems, and even some which go beyond the scope of the currently available theory.

In this paper, we present a backward deep BSDE method applied to Forward Backward Stochastic Differential Equations (FBSDE) with given terminal condition at maturity that time-steps the BSDE backwards. We present an application of this method to a nonlinear pricing problem - the differential rates problem. To time-step the BSDE backward, one needs to solve a nonlinear problem. For the differential rates problem, we derive an exact solution of this time-step problem and a Taylor-based approximation. Previously backward deep BSDE methods only treated zero or linear generators. While a Taylor approach for nonlinear generators was previously mentioned, it had not been implemented or applied, while we apply our method to nonlinear generators and derive details and present results. Likewise, previously backward deep BSDE methods were presented for fixed initial risk factor values $X_0$ only, while we present a version with random $X_0$ and a version that learns portfolio values at intermediate times as well. The method is able to solve nonlinear FBSDE problems in high dimensions.

During the COVID-19 pandemic of 2019/2020, authorities have used temporary ad-hoc policy measures, such as lockdowns and mass quarantines, to slow its transmission. However, the consequences of widespread use of these unprecedented measures are poorly understood. To contribute to the understanding of the economic and human consequences of such policy measures, we therefore construct a mathematical model of an economy under the impact of a pandemic, select parameter values to represent the global economy under the impact of COVID-19, and perform numerical experiments by simulating a large number of possible policy responses. By varying the starting date of the policy intervention in the simulated scenarios, we find that the most effective policy intervention occurs around the time when the number of active infections is growing at its highest rate -- that is, the results suggest that the most severe measures should only be implemented when the disease is sufficiently spread. The intensity of the intervention, above a certain threshold, does not appear to have a great impact on the outcomes in our simulations, due to the strongly concave relationship that we identify between production shortfall and infection rate reductions. Our experiments further suggest that the intervention should last until after the peak established by the reduced infection rate, which implies that stricter policies should last longer. The model and its implementation, along with the general insights from our policy experiments, may help policymakers design effective emergency policy responses in the face of a serious pandemic, and contribute to our understanding of the relationship between the economic growth and the spread of infectious diseases.

The considerable cessation of human activities during the COVID-19 pandemic has affected global energy use and CO2 emissions. Here we show the unprecedented decrease in global fossil CO2 emissions from January to April 2020 was of 7.8% (938 Mt CO2 with a +6.8% of 2-{\sigma} uncertainty) when compared with the period last year. In addition other emerging estimates of COVID impacts based on monthly energy supply or estimated parameters, this study contributes to another step that constructed the near-real-time daily CO2 emission inventories based on activity from power generation (for 29 countries), industry (for 73 countries), road transportation (for 406 cities), aviation and maritime transportation and commercial and residential sectors emissions (for 206 countries). The estimates distinguished the decline of CO2 due to COVID-19 from the daily, weekly and seasonal variations as well as the holiday events. The COVID-related decreases in CO2 emissions in road transportation (340.4 Mt CO2, -15.5%), power (292.5 Mt CO2, -6.4% compared to 2019), industry (136.2 Mt CO2, -4.4%), aviation (92.8 Mt CO2, -28.9%), residential (43.4 Mt CO2, -2.7%), and international shipping (35.9Mt CO2, -15%). Regionally, decreases in China were the largest and earliest (234.5 Mt CO2,-6.9%), followed by Europe (EU-27 & UK) (138.3 Mt CO2, -12.0%) and the U.S. (162.4 Mt CO2, -9.5%). The declines of CO2 are consistent with regional nitrogen oxides concentrations observed by satellites and ground-based networks, but the calculated signal of emissions decreases (about 1Gt CO2) will have little impacts (less than 0.13ppm by April 30, 2020) on the overserved global CO2 concertation. However, with observed fast CO2 recovery in China and partial re-opening globally, our findings suggest the longer-term effects on CO2 emissions are unknown and should be carefully monitored using multiple measures.

We constructed a near-real-time daily CO2 emission dataset, namely the Carbon Monitor, to monitor the variations of CO2 emissions from fossil fuel combustion and cement production since January 1st 2019 at national level with near-global coverage on a daily basis, with the potential to be frequently updated. Daily CO2 emissions are estimated from a diverse range of activity data, including: hourly to daily electrical power generation data of 29 countries, monthly production data and production indices of industry processes of 62 countries/regions, daily mobility data and mobility indices of road transportation of 416 cities worldwide. Individual flight location data and monthly data were utilised for aviation and maritime transportation sectors estimates. In addition, monthly fuel consumption data that corrected for daily air temperature of 206 countries were used for estimating the emissions from commercial and residential buildings. This Carbon Monitor dataset manifests the dynamic nature of CO2 emissions through daily, weekly and seasonal variations as influenced by workdays and holidays, as well as the unfolding impacts of the COVID-19 pandemic. The Carbon Monitor near-real-time CO2 emission dataset shows a 7.8% decline of CO2 emission globally from Jan 1st to Apr 30th in 2020 when compared with the same period in 2019, and detects a re-growth of CO2 emissions by late April which are mainly attributed to the recovery of economy activities in China and partial easing of lockdowns in other countries. Further, this daily updated CO2 emission dataset could offer a range of opportunities for related scientific research and policy making.

We develop a methodology which replicates in great accuracy the FTSE Russell indexes reconstitutions, including the quarterly rebalancings due to new initial public offerings (IPOs). While using only data available in the CRSP US Stock database for our index reconstruction, we demonstrate the accuracy of this methodology by comparing it to the original Russell US indexes for the time period between 1989 to 2019. A python package that generates the replicated indexes is also provided.

As an application, we use our index reconstruction protocol to compute the permanent and temporary price impact on the Russell 3000 annual additions and deletions, and on the quarterly additions of new IPOs . We find that the index portfolios following the Russell 3000 index and rebalanced on an annual basis are overall more crowded than those following the index on a quarterly basis. This phenomenon implies that transaction costs of indexing strategies could be significantly reduced by buying new IPOs additions in proximity to quarterly rebalance dates.

In this paper, an approximate version of the Barndorff-Nielsen and Shephard model, driven by a Brownian motion and a L\'evy subordinator, is formulated. The first-exit time of the log-return process for this model is analyzed. It is shown that with certain probability, the first-exit time process of the log-return is decomposable into the sum of the first exit time of the Brownian motion with drift, and the first exit time of a L\'evy subordinator with drift. Subsequently, the probability density functions of the first exit time of some specific L\'evy subordinators, connected to stationary, self-decomposable variance processes, are studied. Analytical expressions of the probability density function of the first-exit time of three such L\'evy subordinators are obtained in terms of various special functions. The results are implemented to empirical S&P 500 dataset.

Green hydrogen can help to decarbonize transportation, but its power sector interactions are not well understood. It may contribute to integrating variable renewable energy sources if production is sufficiently flexible in time. Using an open-source co-optimization model of the power sector and four options for supplying hydrogen at German filling stations, we find a trade-off between energy efficiency and temporal flexibility: for lower shares of renewables and hydrogen, more energy-efficient and less flexible small-scale on-site electrolysis is optimal. For higher shares of renewables and/or hydrogen, more flexible but less energy-efficient large-scale hydrogen supply chains gain importance as they allow disentangling hydrogen production from demand via storage. Liquid hydrogen emerges as particularly beneficial, followed by liquid organic hydrogen carriers and gaseous hydrogen. Large-scale hydrogen supply chains can deliver substantial power sector benefits, mainly through reduced renewable surplus generation. Energy modelers and system planners should consider the distinct flexibility characteristics of hydrogen supply chains in more detail when assessing the role of green hydrogen in future energy transition scenarios.

A Hidden Markov Model for intraday momentum trading is presented which specifies a latent momentum state responsible for generating the observed securities' noisy returns. Existing momentum trading models suffer from time-lagging caused by the delayed frequency response of digital filters. Time-lagging results in a momentum signal of the wrong sign, when the market changes trend direction. A key feature of this state space formulation, is no such lagging occurs, allowing for accurate shifts in signal sign at market change points. The number of latent states in the model is estimated using three techniques, cross validation, penalized likelihood criteria and simulation-based model selection for the marginal likelihood. All three techniques suggest either 2 or 3 hidden states. Model parameters are then found using Baum-Welch and Markov Chain Monte Carlo, whilst assuming a single (discretized) univariate Gaussian distribution for the emission matrix. Often a momentum trader will want to condition their trading signals on additional information. To reflect this, learning is also carried out in the presence of side information. Two sets of side information are considered, namely a ratio of realized volatilities and intraday seasonality. It is shown that splines can be used to capture statistically significant relationships from this information, allowing returns to be predicted. An Input Output Hidden Markov Model is used to incorporate these univariate predictive signals into the transition matrix, presenting a possible solution for dealing with the signal combination problem. Bayesian inference is then carried out to predict the securities $t+1$ return using the forward algorithm. Simple modifications to the current framework allow for a fully non-parametric model with asynchronous prediction.

This paper aims to make a new contribution to the study of lifetime ruin problem by considering investment in two hedge funds with high-watermark fees and drift uncertainty. Due to multi-dimensional performance fees that are charged whenever each fund profit exceeds its historical maximum, the value function is expected to be multi-dimensional. New mathematical challenges arise as the standard dimension reduction cannot be applied, and the convexity of the value function and Isaacs condition may not hold in our ruin probability minimization problem with drift uncertainty. We propose to employ the stochastic Perron's method to characterize the value function as the unique viscosity solution to the associated Hamilton Jacobi Bellman (HJB) equation without resorting to the proof of dynamic programming principle. The required comparison principle is also established in our setting to close the loop of stochastic Perron's method.

A major part of the balance sheets of the largest US banks consists of credit card portfolios. Hence, managing the charge-off rates is a vital task for the profitability of the credit card industry. Different macroeconomic conditions affect individuals' behavior in paying down their debts. In this paper, we propose an expert system for loss forecasting in the credit card industry using macroeconomic indicators. We select the indicators based on a thorough review of the literature and experts' opinions covering all aspects of the economy, consumer, business, and government sectors. The state of the art machine learning models are used to develop the proposed expert system framework. We develop two versions of the forecasting expert system, which utilize different approaches to select between the lags added to each indicator. Among 19 macroeconomic indicators that were used as the input, six were used in the model with optimal lags, and seven indicators were selected by the model using all lags. The features that were selected by each of these models covered all three sectors of the economy. Using the charge-off data for the top 100 US banks ranked by assets from the first quarter of 1985 to the second quarter of 2019, we achieve mean squared error values of 1.15E-03 and 1.04E-03 using the model with optimal lags and the model with all lags, respectively. The proposed expert system gives a holistic view of the economy to the practitioners in the credit card industry and helps them to see the impact of different macroeconomic conditions on their future loss.

We study an $N$-player and a mean field exponential utility game. Each player manages two stocks; one is driven by an individual shock and the other is driven by a common shock. Moreover, each player is concerned not only with her own terminal wealth but also with the relative performance of her competitors. We use the probabilistic approach to study these two games. We show the unique equilibrium of the $N$-player game and the mean field game can be characterized by a novel multi-dimensional FBSDE with quadratic growth and a novel mean-field FBSDEs, respectively. The well-posedness result and the convergence result are established.

Family history is usually seen as a significant factor insurance companies look at when applying for a life insurance policy. Where it is used, family history of cardiovascular diseases, death by cancer, or family history of high blood pressure and diabetes could result in higher premiums or no coverage at all. In this article, we use massive (historical) data to study dependencies between life length within families. If joint life contracts (between a husband and a wife) have been long studied in actuarial literature, little is known about child and parents dependencies. We illustrate those dependencies using 19th century family trees in France, and quantify implications in annuities computations. For parents and children, we observe a modest but significant positive association between life lengths. It yields different estimates for remaining life expectancy, present values of annuities, or whole life insurance guarantee, given information about the parents (such as the number of parents alive). A similar but weaker pattern is observed when using information on grandparents.

We propose a novel model for describing the spreading processes, in particular, epidemics. Our model is an extension of the SIQR (susceptible-infected-quarantine-recovered) and SIRP (susceptible-infected-recovered-pathogen) models used earlier to describe various scenarios of epidemic spread. As compared to the basic SIR model, our model takes into account two possible routes of virus transmission: direct from the infected compartment to the susceptible compartment and indirect via some intermediate medium or fomites. The transmission rates are estimated in terms of the average distances between the individuals in selected social environments and characteristic relaxation times. We also introduce a resource activation function that reflects the load of the epidemics on economics and the limited capacity of the medical infrastructure. Our model brings an advantage of building various control strategies to minimize the effect of the epidemic and can be applied to modeling the recent COVID-19 outbreak.

We provide a general framework for no-arbitrage concepts in topological vector lattices, which covers many of the well-known no-arbitrage concepts as particular cases. The main structural condition which we impose is that the outcomes of trading strategies with initial wealth zero and those with positive initial wealth have the structure of a convex cone. As one consequence of our approach, the concepts NUPBR, NAA$_1$ and NA$_1$ may fail to be equivalent in our general setting. Furthermore, we derive abstract versions of the fundamental theorem of asset pricing. We also consider a financial market with semimartingales which does not need to have a num\'{e}raire, and derive results which show the links between the no-arbitrage concepts by only using the theory of topological vector lattices and well-known results from stochastic analysis in a sequence of short proofs.

In this paper we develop numerical pricing methodologies for European style Exchange Options written on a pair of correlated assets, in a market with finite liquidity. In contrast to the standard multi-asset Black-Scholes framework, trading in our market model has a direct impact on the asset's price. The price impact is incorporated into the dynamics of the first asset through a specific trading strategy, as in large trader liquidity model. Two-dimensional Milstein scheme is implemented to simulate the pair of assets prices. The option value is numerically estimated by Monte Carlo with the Margrabe option as controlled variate. Time complexity of these numerical schemes are included. Finally, we provide a deep learning framework to implement this model effectively in a production environment.

This paper studies of the multifractal dynamics in 84 cryptocurrencies. It fills an important gap in the literature, by studying this market using two alternative multi-scaling methodologies. We find compelling evidence that cryptocurrencies have different degree of long range dependence, and --more importantly -- follow different stochastic processes. Some of them follow models closer to monofractal fractional Gaussian noises, while others exhibit complex multifractal dynamics. Regarding the source of multifractality, our results are mixed. Time series shuffling produces a reduction in the level of multifractality, but not enough to offset it. We find an association of kurtosis with multifractality.

This article presents a new method for forecasting Value at Risk. Convolutional neural networks can do time series forecasting, since they can learn local patterns in time. A simple modification enables them to forecast not the mean, but arbitrary quantiles of the distribution, and thus allows them to be applied to VaR-forecasting. The proposed model can learn from the price history of different assets, and it seems to produce fairly accurate forecasts.

This paper begins with a study on the dual representations of risk and regret measures and their impact on modeling multistage decision making under uncertainty. A relationship between risk envelopes and regret envelopes is established by using the Lagrangian duality theory. Such a relationship opens a door to a decomposition scheme, called progressive hedging, for solving multistage risk minimization and regret minimization problems. In particular, the classical progressive hedging algorithm is modified in order to handle a new class of linkage constraints that arises from reformulations and other applications of risk and regret minimization problems. Numerical results are provided to show the efficiency of the progressive hedging algorithms.

This paper studies a robust portfolio optimization problem under the multi-factor volatility model introduced by Christoffersen et al. (2009). The optimal strategy is derived analytically under the worst-case scenario with or without derivative trading. To illustrate the effects of ambiguity, we compare our optimal robust strategy with some strategies that ignore the information of uncertainty, and provide the corresponding welfare analysis. The effects of derivative trading to the optimal portfolio selection are also discussed by considering alternative strategies. Our study is further extended to the cases with jump risks in asset price and correlated volatility factors, respectively. Numerical experiments are provided to demonstrate the behavior of the optimal portfolio and utility loss.

Fire sales are among the major drivers of market instability in modern financial systems. Due to iterated distressed selling and the associated price impact, initial shocks to some institutions can be amplified dramatically through the network induced by portfolio overlaps. In this paper we develop models that allow us to investigate central characteristics that drive or hinder the propagation of distress. We investigate single systems as well as ensembles of systems that are alike, where similarity is measured in terms of the empirical distribution of all defining properties of a system. This asymptotic approach ensures a great deal of robustness to statistical uncertainty and temporal fluctuations, and we give various applications. A natural characterization of systems resilient to fire sales emerges, and we provide explicit criteria that regulators may exploit in order to assess the stability of any system. Moreover, we propose risk management guidelines in form of minimal capital requirements, and we investigate the effect of portfolio diversification and portfolio overlap. We test our results by Monte Carlo simulations for exemplary configurations.

In this paper, we introduce a new time series model having the stochastic exponential tail. This model is constructed based on the Normal Tempered Stable distribution with a time varying parameter. The model captures the stochastic exponential tail which generates the volatility smile effect and volatility term structure in option pricing. Moreover, the model describes the time varying volatility of volatility. We empirically show the stochastic skewness and stochastic kurtosis by applying the model to analyze S\&P 500 index return data. We present Monte-Carlo simulation technique for a parameter calibration of the model for the S\&P 500 option prices. By the calibration, we can see that the stochastic exponential tail makes the model better to analyze the market option prices.

By 2030 Austria aims to meet \SI{100}{\percent} of its electricity demand from domestic renewable sources, predominantly from wind and solar energy. While wind power reduces \COO emissions, it is also connected to negative impacts at the local level. such as interference with landscape aesthetics. Nevertheless, wind power comes at lower system integration cost than solar power, so that it effectively reduces system cost. We quantify the opportunity cost of replacing wind turbines with solar power, using the power system model medea. Our findings suggest that these cost of undisturbed landscapes are considerable, particularly when PV is not entirely rolled out as utility scale, open space installations. The opportunity cost is likely high enough to allow for significant compensation of the ones affected by local wind turbine externalities.

Especially in the insurance industry interest rate models play a crucial role e.g. to calculate the insurance company's liabilities, performance scenarios or risk measures. A prominant candidate is the 2-Additive-Factor Gaussian Model (Gauss2++) - in a different representation also known as the 2-Factor Hull-White model. In this paper, we propose a framework to estimate the model such that it can be applied under the risk neutral and the real world measure in a consistent manner. We first show that any progressive and square-integrable function can be used to specify the change of measure without loosing the analytic tractability of e.g. zero-coupon bond prices in both worlds. We further propose two time dependent candidates, which are easy to calibrate: a step and a linear function. They represent two variants of our framework and distinguish between a short and a long term risk premium, which allows to regularize the interest rates in the long horizon. We apply both variants to historical data and show that they indeed produce realistic and much more stable long term interest rate forecast than the usage of a constant function. This stability over time would translate to performance scenarios of e.g. interest rate sensitive fonds and risk measures.

The dependence of world GDP on current energy consumption and total energy produced over the previous period and materialized in the form of production infrastructure is studied. The dependence describes empirical data with high accuracy over the entire observation interval 1965-2018.

We introduce time-inhomogeneous stochastic volatility models, in which the volatility is described by a nonnegative function of a Volterra type continuous Gaussian process that may have extremely rough sample paths. The drift function and the volatility function are assumed to be time-dependent and locally $\omega$-continuous for some modulus of continuity $\omega$. The main results obtained in the paper are sample path and small-noise large deviation principles for the log-price process in a Gaussian model under very mild restrictions. We use these results to study the asymptotic behavior of binary up-and-in barrier options and binary call options.

Financial markets across all asset classes are known to exhibit trends. These trends have been exploited by traders for decades. Here, we empirically measure when trends revert, based on 30 years of daily futures prices for equity indices, interest rates, currencies and commodities. We find that trends tend to revert once they reach a critical level of statistical significance. Based on polynomial regression, we carefully measure this critical level. We find that it is universal across asset classes and has a universal scaling behavior, as the trend's time horizon runs from a few days to several years. The corresponding regression coefficients are small, but statistically highly significant, as confirmed by bootstrapping and out-of-sample testing. Our results signal to investors when to exit a trend. They also reveal how markets have become more efficient over the decades. Moreover, they point towards a potential deep analogy between financial markets and critical phenomena: our analysis supports the conjecture that financial markets can be modeled as statistical mechanical ensembles of Buy/Sell orders near critical points. In this analogy, the trend strength plays the role of an order parameter, whose dynamcis is described by a Langevin equation.

We discuss the impact of a Covid-like shock on a simple toy economy, described by the Mark-0 Agent-Based Model that we developed and discussed in a series of previous papers. We consider a mixed supply and demand shock, and show that depending on the shock parameters (amplitude and duration), our toy economy can display V-shaped, U-shaped or W-shaped recoveries, and even an L-shaped output curve with permanent output loss. This is due to the existence of a self-sustained "bad" state of the economy. We then discuss two policies that attempt to moderate the impact of the shock: giving easy credit to firms, and the so-called helicopter money, i.e. injecting new money into the households savings. We find that both policies are effective if strong enough, and we highlight the potential danger of terminating these policies too early. While we only discuss a limited number of scenarios, our model is flexible and versatile enough to allow for a much wider exploration, thus serving as a useful tool for the qualitative understanding of post-Covid recovery.