Arxiv:1809.04341V1 [Quant-Ph] 12 Sep 2018 Progress, See [18])

Total Page:16

File Type:pdf, Size:1020Kb

Arxiv:1809.04341V1 [Quant-Ph] 12 Sep 2018 Progress, See [18]) Stochastic Integral Representation for the Dynamics of Disordered Systems Ivana Kureciˇ c´1 and Tobias J. Osborne2 1Max Planck Institute of Quantum Optics, Hans-Kopfermann-Str. 1, 85748 Garching, Germany 2Institut fur¨ Theoretische Physik, Leibniz Universitat¨ Hannover, Appelstr. 2, 30167 Hannover, Germany The dynamics of interacting quantum systems in the presence of disorder is studied and an exact rep- resentation for disorder-averaged quantities via Itoˆ stochastic calculus is obtained. The stochastic integral representation affords many advantages, including amenability to analytic approximation, applicability to interacting systems, and compatibility with existing tensor network methods. The integral may be ex- panded to produce a series of approximations, the first of which already includes all diffusive corrections and, further, is manifestly completely positive. The addition of fluctuations leads to a convergent series of systematic corrections. As examples, expressions for the density of states, spectral form factor, and out-of-time-order correlators for the Anderson model are obtained. Since Anderson’s disovery of the localisation phe- averaged propagator of an arbitrary quantum system. The nomenon [1], quantum systems in the presence of disor- derivation of this representation is reminiscent, in parts, of der have attracted considerable interest. The localisation field-theoretic approaches to disordered systems and also phenomenon is a result of quantum coherent effects, thus separately to recent calculations [22, 23] of Prosen and presenting many challenges and exhibiting rich physics. coworkers. However, there are crucial differences. For ex- During the past 60 years there has been excellent progress ample, there is no mapping to the nonlinear sigma model. in the study of disordered systems in the single-particle Also, the representation is provably exact and not an ap- regime, culminating in a solid understanding of these phys- proximation. We exploit the integral to produce a series ical phenomena (see, e.g., the recent review [2]). These ad- of approximations, the first of which already includes all vances are the result of several powerful techniques which diffusive corrections and, further, is manifestly the result have been developed to quantitatively study disordered sys- of a completely positive evolution, thus conserving proba- tems, including numerical methods [3], the renormalisa- bilities. The addition of fluctuations leads to a convergent tion group [4], and field-theoretic methods [5] such as the series of systematic corrections. Using this expansion we replica trick [6, 7], the supersymmetric method [8, 9], and compute the density of states and out-of-time-order corre- the Keldysh formalism [10–12]. lations for the Anderson model. With the recent discovery of the intriguing phenomenon Disordered quantum systems, Brownian motions, and of many body localisation (MBL) [13] (see also the recent stochastic calculus.—Here we introduce the systems un- reviews [14–17]), renewed interest in disordered systems der consideration and give a brief notational summary of has emerged. As a collective many particle effect, MBL stochastic calculus. We consider quantum systems with a is proportionally more challenging to study than its single- Hamiltonian of the form m particle counterpart. The most progress in this field has X come from perturbation theory and direct numerical sim- H(x) = H0 + xjDj; (1) ulation via Monte Carlo sampling. In general, these tech- j=1 niques are limited either to weak disorder or finite system where H0 is a fixed Hamiltonian (typically a kinetic energy sizes (of the order of 20 spins), and few samples. Ex- term), D represent disordered terms (e.g., a local mag- ploiting the most powerful field-theoretic approaches to j netic field or potential energy terms), and xj are random study the strongly interacting quantum spin systems ex- variables drawn from the Gaussian distribution [24] with hibiting MBL appears to be a difficult task (for some recent the probability density function arXiv:1809.04341v1 [quant-ph] 12 Sep 2018 progress, see [18]). 2 − x Some intriguing new directions — emerging from holo- e 2γ2 gγ (x) ≡ p ; graphic arguments in high energy physics [19] — are also 2πγ being explored in the study of the complex dynamics of quantum systems such as the Anderson model [20]. Here a where γ is its standard deviation. This class of models is central role is played by out-of-time-order correlation func- sufficiently general to describe a diverse variety of models tions (OTOCs) [21] as signatures of quantum chaos. The from single impurity models, the Anderson model in arbi- behaviour of OTOCs for complex disordered systems is a trary dimensions, as well as MBL systems. We illustrate central goal in the study of new phases of disordered mat- the results in terms of the one-dimensional Anderson tight- ter. binding model with periodic boundary conditions, where Motivated by the twin challenges of MBL and calculat- N X ing OTOCs for disordered systems, in this Letter we intro- H0 ≡ 2I − jj + 1ihjj + jjihj + 1j; (2) duce a stochastic integral representation for the disorder- j=1 2 and Dj ≡ jjihjj. We study disorder-averaged dynamical where the product is taken from right to left and δ(f(x)) = quantities, such as the disorder-averaged propagator: δ(x2 −x1) ··· δ(xn −xn−1). The second step is to use the itH identity S(t) = Ex[e ]; (3) 1 Z 1 δ(x) ≡ eikx dk; (9) whose Fourier transform yields the density of states (DOS), 2π the density operator: −∞ and to carry out the integral over the xj variables. This −itH itH ρ(t) = Ex[e j 0ih 0je ]; (4) leaves an integral over k1; k2; : : : ; kn−1: 1 Z nγ2 Pn 2 spectral form factors: − 2 j=1(kj−1−kj ) S(t) = 1 e Fe(k) dk; n− 2 (−β−it)H 2k (2π) γ Ex[j tr(e )j ]; (5) (10) and OTOCs: where we’ve introduced two additional k variables, k0 = kn = 0, as well as the function Ex[hA(0)B(t)C(0)D(t)i]: (6) n Y These quantities are intimately interrelated (see, e.g., [25]), Fe(k) ≡ Fej(kj−1 − kj); (11) and can all be calculated in terms of S(t) by taking tensor j=1 copies of the Hamiltonian [26]. where Z 1 p 2 − p 1 x −i n γy it We make use of stochastic calculus, but it is not at 2nγ j 2 (A+xj B) Fej(y) = e e n dxj: (12) all necessary to be familiar with this formalism, as all −∞ derivations can be understood with little more than a pass- The third step is to expand the exponents of the operator ex- ing familiarity with Gaussian integrals and a tolerance ponentials in F s, collect terms to O(1=n2), and then carry for lengthy derivations using discretisations [27]. At this e out the integrals over x : stage it is sufficient to comment that a Brownian motion or j Z n Wiener process Wt is characterised by the following four Y F (kj−1 − kj) dµ, (13) properties: (1) W0 = 0; (2) Wt is (almost surely) continu- j=1 ous; (3) Wt has independent increments; and (4) Wt − Ws is distributed according to the normal distribution with a where mean of zero and variance t − s for 0 ≤ s ≤ t. We will it t2 F (y) ≡ + (A + iγ2nyB) − (nγ2 − n2γ4y2)B2 also encounter stochastic differential equations (SDE) of I n 2n2 the form dy = f ds + g dW . At this stage, it is sufficient (14) to regard these as equations that, upon discretisation, define and s new random variables y in terms of W . For further details 2 n−1 2 n(nγ ) − nγ kT Mk on stochastic calculus see, e.g., [28, 29]. dµ ≡ e 2 dk1 ··· dkn−1; (15) (2π)n−1 A stochastic integral representation.— In this section we summarise the salient features of the derivation of our in- with tegral representation, with an emphasis on the physical 0 2 −1 0 0 ··· 01 foundations of the argument. As explained previously, for B−1 2 −1 0 ··· 0C Hamiltonians of the form (1), it is sufficient to restrict our B 0 −1 2 −1 0C M = B C = 2I − Pn−1: (16) attention to the study of the quantity B . C @ . .. A itH(x) S(t) ≡ Ex[e ]: (7) 0 ··· 0 −1 2 For simplicity, we focus on Hamiltonians of the form In physical terms, µ is the (discretisation of the) equilib- rium probability distribution for a free particle with the H(x) = A+xB, where x is a single Gaussian-distributed 2 random variable, and A and B are N × N matrices. (The Hamiltonian H = −∇ , which is diffusively moving on extension to Hamiltonians of the form from Eq. (1) is the interval [0; 1] with vanishing Dirichlet boundary condi- entirely straightforward and requires no additional tech- tions. niques.) We derive our representation in four steps. The By employing the approximation itH it 2 t2 2 2 first step is to break the propagator e into n small pieces n (A+inγ (kj−1−kj )B)− 2n γ B n F (y) ≈ e ; (17) it H e n , introduce n independent Gaussian-distributed valid to O(1=n2), we have already arrived at an approxi- variables xj, j = 1; 2; : : : ; n, and enforce equality via mation of great utility in tensor-network simulations (this delta functions: will be the subject of a forthcoming paper): n n Z kxk2 Z 2 1 − Y it (A+x B) Y it A−tγ2(k −k )B− t γ2B2 p S(t) = p e 2nγ2 δ(f(x)) e n j dx; S(t) ≈ e n j−1 j 2n dµ+O(1= n): 2πγ j=1 j=1 (8) (18) 3 The final step is to recognise the integral measure in the The stochastic Dyson series: diffusions and limit of n ! 1 as the path measure for the Brownian fluctuations.—In this section we describe a Dyson se- bridge [28], which is a continuous-time stochastic process ries procedure to develop the integral in a power series with the same conditional probability distribution as the in the disorder parameter γ.
Recommended publications
  • Markovian Bridges: Weak Continuity and Pathwise Constructions
    The Annals of Probability 2011, Vol. 39, No. 2, 609–647 DOI: 10.1214/10-AOP562 c Institute of Mathematical Statistics, 2011 MARKOVIAN BRIDGES: WEAK CONTINUITY AND PATHWISE CONSTRUCTIONS1 By Lo¨ıc Chaumont and Geronimo´ Uribe Bravo2,3 Universit´ed’Angers and Universidad Nacional Aut´onoma de M´exico A Markovian bridge is a probability measure taken from a disin- tegration of the law of an initial part of the path of a Markov process given its terminal value. As such, Markovian bridges admit a natural parameterization in terms of the state space of the process. In the context of Feller processes with continuous transition densities, we construct by weak convergence considerations the only versions of Markovian bridges which are weakly continuous with respect to their parameter. We use this weakly continuous construction to provide an extension of the strong Markov property in which the flow of time is reversed. In the context of self-similar Feller process, the last result is shown to be useful in the construction of Markovian bridges out of the trajectories of the original process. 1. Introduction and main results. 1.1. Motivation. The aim of this article is to study Markov processes on [0,t], starting at x, conditioned to arrive at y at time t. Historically, the first example of such a conditional law is given by Paul L´evy’s construction of the Brownian bridge: given a Brownian motion B starting at zero, let s s bx,y,t = x + B B + (y x) . s s − t t − t arXiv:0905.2155v3 [math.PR] 14 Mar 2011 Received May 2009; revised April 2010.
    [Show full text]
  • Patterns in Random Walks and Brownian Motion
    Patterns in Random Walks and Brownian Motion Jim Pitman and Wenpin Tang Abstract We ask if it is possible to find some particular continuous paths of unit length in linear Brownian motion. Beginning with a discrete version of the problem, we derive the asymptotics of the expected waiting time for several interesting patterns. These suggest corresponding results on the existence/non-existence of continuous paths embedded in Brownian motion. With further effort we are able to prove some of these existence and non-existence results by various stochastic analysis arguments. A list of open problems is presented. AMS 2010 Mathematics Subject Classification: 60C05, 60G17, 60J65. 1 Introduction and Main Results We are interested in the question of embedding some continuous-time stochastic processes .Zu;0Ä u Ä 1/ into a Brownian path .BtI t 0/, without time-change or scaling, just by a random translation of origin in spacetime. More precisely, we ask the following: Question 1 Given some distribution of a process Z with continuous paths, does there exist a random time T such that .BTCu BT I 0 Ä u Ä 1/ has the same distribution as .Zu;0Ä u Ä 1/? The question of whether external randomization is allowed to construct such a random time T, is of no importance here. In fact, we can simply ignore Brownian J. Pitman ()•W.Tang Department of Statistics, University of California, 367 Evans Hall, Berkeley, CA 94720-3860, USA e-mail: [email protected]; [email protected] © Springer International Publishing Switzerland 2015 49 C. Donati-Martin et al.
    [Show full text]
  • Constructing a Sequence of Random Walks Strongly Converging to Brownian Motion Philippe Marchal
    Constructing a sequence of random walks strongly converging to Brownian motion Philippe Marchal To cite this version: Philippe Marchal. Constructing a sequence of random walks strongly converging to Brownian motion. Discrete Random Walks, DRW’03, 2003, Paris, France. pp.181-190. hal-01183930 HAL Id: hal-01183930 https://hal.inria.fr/hal-01183930 Submitted on 12 Aug 2015 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Discrete Mathematics and Theoretical Computer Science AC, 2003, 181–190 Constructing a sequence of random walks strongly converging to Brownian motion Philippe Marchal CNRS and Ecole´ normale superieur´ e, 45 rue d’Ulm, 75005 Paris, France [email protected] We give an algorithm which constructs recursively a sequence of simple random walks on converging almost surely to a Brownian motion. One obtains by the same method conditional versions of the simple random walk converging to the excursion, the bridge, the meander or the normalized pseudobridge. Keywords: strong convergence, simple random walk, Brownian motion 1 Introduction It is one of the most basic facts in probability theory that random walks, after proper rescaling, converge to Brownian motion. However, Donsker’s classical theorem [Don51] only states a convergence in law.
    [Show full text]
  • Arxiv:2004.05394V1 [Math.PR] 11 Apr 2020 Not Equal
    SIZE AND SHAPE OF TRACKED BROWNIAN BRIDGES ABDULRAHMAN ALSOLAMI, JAMES BURRIDGE AND MICHALGNACIK Abstract. We investigate the typical sizes and shapes of sets of points obtained by irregularly tracking two-dimensional Brownian bridges. The tracking process consists of observing the path location at the arrival times of a non-homogeneous Poisson process on a finite time interval. The time varying intensity of this observation process is the tracking strategy. By analysing the gyration tensor of tracked points we prove two theorems which relate the tracking strategy to the average gyration radius, and to the asphericity { a measure of how non-spherical the point set is. The act of tracking may be interpreted either as a process of observation, or as process of depositing time decaying \evidence" such as scent, environmental disturbance, or disease particles. We present examples of different strategies, and explore by simulation the effects of varying the total number of tracking points. 1. Introduction Understanding the statistical properties of human and animal movement processes is of interest to ecologists [1, 2, 3], epidemiologists [4, 5, 6], criminologists [7], physicists and mathematicians [8, 9, 10, 11], including those interested in the evolution of human culture and language [12, 13, 14]. Advances in information and communication technologies have allowed automated collection of large numbers of human and animal trajectories [15, 16], allowing real movement patterns to be studied in detail and compared to idealised mathematical models. Beyond academic study, movement data has important practical applications, for example in controlling the spread of disease through contact tracing [6]. Due to the growing availability and applications of tracking information, it is useful to possess a greater analytical understanding of the typical shape and size characteristics of trajectories which are observed, or otherwise emit information.
    [Show full text]
  • On Conditionally Heteroscedastic Ar Models with Thresholds
    Statistica Sinica 24 (2014), 625-652 doi:http://dx.doi.org/10.5705/ss.2012.185 ON CONDITIONALLY HETEROSCEDASTIC AR MODELS WITH THRESHOLDS Kung-Sik Chan, Dong Li, Shiqing Ling and Howell Tong University of Iowa, Tsinghua University, Hong Kong University of Science & Technology and London School of Economics & Political Science Abstract: Conditional heteroscedasticity is prevalent in many time series. By view- ing conditional heteroscedasticity as the consequence of a dynamic mixture of in- dependent random variables, we develop a simple yet versatile observable mixing function, leading to the conditionally heteroscedastic AR model with thresholds, or a T-CHARM for short. We demonstrate its many attributes and provide com- prehensive theoretical underpinnings with efficient computational procedures and algorithms. We compare, via simulation, the performance of T-CHARM with the GARCH model. We report some experiences using data from economics, biology, and geoscience. Key words and phrases: Compound Poisson process, conditional variance, heavy tail, heteroscedasticity, limiting distribution, quasi-maximum likelihood estimation, random field, score test, T-CHARM, threshold model, volatility. 1. Introduction We can often model a time series as the sum of a conditional mean function, the drift or trend, and a conditional variance function, the diffusion. See, e.g., Tong (1990). The drift attracted attention from very early days, although the importance of the diffusion did not go entirely unnoticed, with an example in ecological populations as early as Moran (1953); a systematic modelling of the diffusion did not seem to attract serious attention before the 1980s. For discrete-time cases, our focus, as far as we are aware it is in the econo- metric and finance literature that the modelling of the conditional variance has been treated seriously, although Tong and Lim (1980) did include conditional heteroscedasticity.
    [Show full text]
  • Arxiv:1910.05067V4 [Physics.Comp-Ph] 1 Jan 2021
    A FINITE-VOLUME METHOD FOR FLUCTUATING DYNAMICAL DENSITY FUNCTIONAL THEORY ANTONIO RUSSO†, SERGIO P. PEREZ†, MIGUEL A. DURAN-OLIVENCIA,´ PETER YATSYSHIN, JOSE´ A. CARRILLO, AND SERAFIM KALLIADASIS Abstract. We introduce a finite-volume numerical scheme for solving stochastic gradient-flow equations. Such equations are of crucial importance within the framework of fluctuating hydrodynamics and dynamic density func- tional theory. Our proposed scheme deals with general free-energy functionals, including, for instance, external fields or interaction potentials. This allows us to simulate a range of physical phenomena where thermal fluctuations play a crucial role, such as nucleation and other energy-barrier crossing transitions. A positivity-preserving algorithm for the density is derived based on a hybrid space discretization of the deterministic and the stochastic terms and different implicit and explicit time integrators. We show through numerous applications that not only our scheme is able to accurately reproduce the statistical properties (structure factor and correlations) of the physical system, but, because of the multiplicative noise, it allows us to simulate energy barrier crossing dynamics, which cannot be captured by mean-field approaches. 1. Introduction The study of fluid dynamics encounters major challenges due to the inherently multiscale nature of fluids. Not surprisingly, fluid dynamics has been one of the main arenas of activity for numerical analysis and fluids are commonly studied via numerical simulations, either at molecular scale, by using molecular dynamics (MD) or Monte Carlo (MC) simulations; or at macro scale, by utilising deterministic models based on the conservation of fundamental quantities, namely mass, momentum and energy. While atomistic simulations take into account thermal fluctuations, they come with an important drawback, the enormous computational cost of having to resolve at least three degrees of freedom per particle.
    [Show full text]
  • Stochastic Flows Associated to Coalescent Processes II
    Stochastic flows associated to coalescent processes II: Stochastic differential equations Dedicated to the memory of Paul-Andr´e MEYER Jean Bertoin(1) and Jean-Fran¸cois Le Gall(2) (1) Laboratoire de Probabilit´es et Mod`eles Al´eatoires and Institut universitaire de France, Universit´e Pierre et Marie Curie, 175, rue du Chevaleret, F-75013 Paris, France. (2) DMA, Ecole normale sup´erieure, 45, rue d'Ulm, F-75005 Paris, France. Summary. We obtain precise information about the stochastic flows of bridges that are associated with the so-called Λ-coalescents. When the measure Λ gives no mass to 0, we prove that the flow of bridges is generated by a stochastic differential equation driven by a Poisson point process. On the other hand, the case Λ = δ0 of the Kingman coalescent gives rise to a flow of coalescing diffusions on the interval [0; 1]. We also discuss a remarkable Brownian flow on the circle which has close connections with the Kingman coalescent Titre. Flots stochastiques associ´es a` des processus de coalescence II : ´equations diff´erentielles stochastiques. R´esum´e. Nous ´etudions les flots de ponts associ´es aux processus de coagulation appel´es Λ-coalescents. Quand la mesure Λ ne charge pas 0, nous montrons que le flot de ponts est engendr´e par une ´equation diff´erentielle stochastique conduite par un processus de Poisson ponctuel. Au contraire, le cas Λ = δ0 du coalescent de Kingman fait apparaitre un flot de diffusions coalescentes sur l'intervalle [0; 1]. Nous ´etudions aussi un flot brownien remarquable sur le cercle, qui est ´etroitement li´e au coalescent de Kingman.
    [Show full text]
  • A Acceleration, 199 Adjacency Matrix, 115 AIC, 212 Akaike Information
    Index A Bayesian sampling, 392 Acceleration, 199 Behavioral dynamics, 218 Adjacency matrix, 115 Bergstrom’s statistic, 345 AIC, 212 Between-and within-person changes, 259 Akaike information criterion, 208 BFGS secant updates, 402 Aliasing, 11 Bias, 107 Aliasing problem, 329 BIC, 212 Ambulatory assessment, 245 Bifurcation, 390 Amplitude, 293 Boundaries, 225 Approximate discrete model, 11 Boundary effect or edge effect, 255–256 ARMAX model, 2 Bowman-Shenton test, 314 Assumption of invertibility, 325 Box-Ljung statistic, 314 Attractor force, 67 Box-Pierce test, 153 Attribute evolution, 117–118 BPS, 67–68 Augmented state space form, 328 Brownian bridge, 412 Autocorrelation, 29 Brownian motion, 59 Auto-covariance function, 372 Business cycle fluctuations, 283 Autoeffects, 38, 184 Business cycle theory, 290 Autoregression, 180 Autoregression coefficients, 169 Autoregressive and moving average (ARMA) C model, 2 CARIMA, 346 CARMA, 1 Causal inference, 106 B Chaotic Lorenz model, 389 Backward operator, 405 Chronometric factor models, 266 Backward pass, 310 Classification and regression tree (CART), 264 Base and input processes, 102 Closed-form solution, 319 Base process, 88 Coevolution model, 120 Bayes factor, 74 Coevolution of networks and behaviors, 112 Bayes formula, 392 COGITO study, 273 Bayesian estimation, 56 Coherence, 370 Bayesian framework, 55–76 Cointegration, 325 Bayesian inference, 68 Competitive equilibrium, 285 Bayesian information criterion, 208 Complete mediation, 197–198 © Springer International Publishing AG, part of Springer Nature
    [Show full text]
  • A Guide to Brownian Motion and Related Stochastic Processes
    Vol. 0 (0000) A guide to Brownian motion and related stochastic processes Jim Pitman and Marc Yor Dept. Statistics, University of California, 367 Evans Hall # 3860, Berkeley, CA 94720-3860, USA e-mail: [email protected] Abstract: This is a guide to the mathematical theory of Brownian mo- tion and related stochastic processes, with indications of how this theory is related to other branches of mathematics, most notably the classical the- ory of partial differential equations associated with the Laplace and heat operators, and various generalizations thereof. As a typical reader, we have in mind a student, familiar with the basic concepts of probability based on measure theory, at the level of the graduate texts of Billingsley [43] and Durrett [106], and who wants a broader perspective on the theory of Brow- nian motion and related stochastic processes than can be found in these texts. Keywords and phrases: Markov process, random walk, martingale, Gaus- sian process, L´evy process, diffusion. AMS 2000 subject classifications: Primary 60J65. Contents 1 Introduction................................. 3 1.1 History ................................ 3 1.2 Definitions............................... 4 2 BM as a limit of random walks . 5 3 BMasaGaussianprocess......................... 7 3.1 Elementarytransformations . 8 3.2 Quadratic variation . 8 3.3 Paley-Wiener integrals . 8 3.4 Brownianbridges........................... 10 3.5 FinestructureofBrownianpaths . 10 arXiv:1802.09679v1 [math.PR] 27 Feb 2018 3.6 Generalizations . 10 3.6.1 FractionalBM ........................ 10 3.6.2 L´evy’s BM . 11 3.6.3 Browniansheets ....................... 11 3.7 References............................... 11 4 BMasaMarkovprocess.......................... 12 4.1 Markovprocessesandtheirsemigroups. 12 4.2 ThestrongMarkovproperty. 14 4.3 Generators .............................
    [Show full text]
  • Convergence of a Randomised Change Point Estimator in GARCH Models
    Journal of Mathematical Finance, 2021, 11, 234-245 https://www.scirp.org/journal/jmf ISSN Online: 2162-2442 ISSN Print: 2162-2434 Convergence of a Randomised Change Point Estimator in GARCH Models George Awiakye-Marfo1, Joseph Mung’atu2, Patrick Weke3 1Department of Mathematics and Statistics, Pan African University, Nairobi, Kenya 2Department of Statistics and Actuarial Science, JKUAT, Nairobi, Kenya 3School of Mathematics, University of Nairobi, Nairobi, Kenya How to cite this paper: Awiakye-Marfo, Abstract G., Mung’atu, J. and Weke, P. (2021) Con- vergence of a Randomised Change Point In this paper, the randomised pseudolikelihood ratio change point estimator Estimator in GARCH Models. Journal of for GARCH models in [1] is employed and its limiting distribution is derived Mathematical Finance, 11, 234-245. as the supremum of a standard Brownian bridge. Data analysis to validate the https://doi.org/10.4236/jmf.2021.112013 estimator is carried out using the United states dollar (USD)-Ghana cedi Received: February 7, 2021 (GHS) daily exchange rate data. The randomised estimator is able to detect Accepted: May 9, 2021 and estimate a single change in the variance structure of the data and pro- Published: May 12, 2021 vides a reference point for historic data analysis. Copyright © 2021 by author(s) and Scientific Research Publishing Inc. Keywords This work is licensed under the Creative GARCH, Randomised, Limiting Distribution, Brownian Bridge, Volatility, Commons Attribution International License (CC BY 4.0). CUSUM, IGARCH, Supremum http://creativecommons.org/licenses/by/4.0/ Open Access 1. Introduction Volatility models are becoming increasingly important due to their role in asset pricing and risk management.
    [Show full text]
  • Multiscale Stochastic Optimization: Modeling Aspects and Scenario Generation
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by International Institute for Applied Systems Analysis (IIASA) Computational Optimization and Applications https://doi.org/10.1007/s10589-019-00135-4 Multiscale stochastic optimization: modeling aspects and scenario generation Martin Glanzer1 · Georg Ch. Pflug1,2 Received: 4 March 2019 © The Author(s) 2019 Abstract Real-world multistage stochastic optimization problems are often characterized by the fact that the decision maker may take actions only at specific points in time, even if relevant data can be observed much more frequently. In such a case there are not only multiple decision stages present but also several observation periods between consec- utive decisions, where profits/costs occur contingent on the stochastic evolution of some uncertainty factors. We refer to such multistage decision problems with encap- sulated multiperiod random costs, as multiscale stochastic optimization problems. In this article, we present a tailor-made modeling framework for such problems, which allows for a computational solution. We first establish new results related to the gen- eration of scenario lattices and then incorporate the multiscale feature by leveraging the theory of stochastic bridge processes. All necessary ingredients to our proposed modeling framework are elaborated explicitly for various popular examples, including both diffusion and jump models. Keywords Stochastic programming · Scenario generation · Bridge process · Stochastic bridge · Diffusion bridge · Lévy bridge · Compound Poisson bridge · Simulation of stochastic bridge · Multiple time scales · Multi-horizon · Multistage stochastic optimization 1 Introduction Optimization models over a large time frame can be classified into two types: B Martin Glanzer [email protected] 1 Department of Statistics and Operations Research (DSOR), University of Vienna, Oskar-Morgenstern-Platz 1, 1090 Vienna, Austria 2 International Institute for Applied Systems Analysis (IIASA), Laxenburg, Austria 123 M.
    [Show full text]
  • Path Transformations in Lévy Processes
    Path transformations in L´evy processes Lo¨ıc Chaumont Laboratoire de Probabilit´es et Mod`eles Al´eatoires, Universit´eParis6 4, Place Jussieu, 75252, Paris cedex 05 E-mail: [email protected] 1. Introduction Let D be the space of c.`a.d. l.`a.g. functions with finite or infinite life time. We are interested in maps from D onto itself which connect the paths of stochastic processes related to L´evy processes. The existence of such path transformations is mainly due to the special properties of the increments of L´evy processes and their main interest is to provide and explain some identities in law between remarkable functionals. In most of the cases such identities are first noticed by pure calculation in particular cases like Brownian motion or Cauchy processes. Then, it sometimes arises that they have pathwise explanation which can be spread to a large class of L´evy processes. Moreover, these path transformations usually lead to more complete identities in law than the initial one. To illustrate this topic, we chose to present two examples. The first one is known as Ver- vaat’s transformation and its aim was to explain the identity in law between the amplitude of the Brownian bridge and the maximum of the normalized Brownian excursion. The second one is due to Bertoin, (1993). It provides a pathwise explanation of the identity in law between the time spent under 0 by a L´evy process up to time 1 and the time at which this process first hits its minimum before 1.
    [Show full text]