Source author record

Ivan Corwin

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

42works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

42 published item(s)

preprint2022arXiv

Exceptional times when the KPZ fixed point violates Johansson's conjecture on maximizer uniqueness

In 2002, Johansson conjectured that the maximum of the Airy$_2$ process minus the parabola $x^2$ is almost surely achieved at a unique location. This result was proved a decade later by Corwin and Hammond; Moreno Flores, Quastel and Remenik; and Pimentel. Up to scaling, the Airy$_2$ process minus the parabola $x^2$ arises as the fixed time spatial marginal of the KPZ fixed point when started from narrow wedge initial data. We extend this maximizer uniqueness result to the fixed time spatial marginal of the KPZ fixed point when begun from any element of a very broad class of initial data. None of these results rules out the possibility that at random times, the KPZ fixed point spatial marginal violates maximizer uniqueness. To understand this possibility, we study the probability that the KPZ fixed point has, at a given time, two or more locations where its value is close to the maximum, obtaining quantitative upper and lower bounds in terms of the degree of closeness for a very broad class of initial data. We also compute a quantity akin to the joint density of the locations of two maximizers and the maximum value. As a consequence, the set of times of maximizer non-uniqueness almost surely has Hausdorff dimension at most two-thirds. Our analysis relies on the exact formula for the distribution function of the KPZ fixed point obtained by Matetski, Quastel and Remenik, the variational formula for the KPZ fixed point involving the Airy sheet constructed by Dauvergne, Ortmann and Virág, and the Brownian Gibbs property for the Airy$_2$ process minus the parabola $x^2$ demonstrated by Corwin and Hammond.

preprint2022arXiv

Harold Widom's work in random matrix theory

This is a survey of Harold Widom's work in random matrices. We start with his pioneering papers on the sine-kernel determinant, continue with his and Craig Tracy's groundbreaking results concerning the distribution functions of random matrix theory, touch on the remarkable universality of the Tracy-Widom distributions in mathematics and physics, and close with Tracy and Widom's remarkable work on the asymmetric simple exclusion process.

preprint2022arXiv

Some recent progress on the stationary measure for the open KPZ equation

This note, dedicated in Harold Widom's memory, is an expanded version of a lecture I gave in fall 2021 at the MSRI program "Universality and Integrability in Random Matrices and Interacting Particle Systems". I will focus on the behavior of the stationary measure for the open KPZ equation, a paradigmatic model for interface growth in contact with boundaries. Much of this will review elements of my joint work with A. Knizel as well as with H. Shen, as well as subsequent works of W. Bryc, A. Kuznetsov, Y. Wang, and J. Wesolowski and of G. Barraquand and P. Le Doussal. The basis for this advance is fundamental work of B. Derrida, M. Evans, V. Hakim and V. Pasquier from 1993, of T. Sasamoto, M. Uchiyama and M. Wadati from 2003, and of W. Bryc and J. Wesolowski from 2010 and 2017. I will try to explain how all of this fits together, without laboring details for the sake of exposition. Though this work does not directly follow from Harold Widom's own work, it (and a great deal of my research) is very much inspired by his and Craig Tracy's work on ASEP.

preprint2022arXiv

Spatial tightness at the edge of Gibbsian line ensembles

Consider a sequence of Gibbsian line ensembles, whose lowest labeled curves (i.e., the edge) have tight one-point marginals. Then, given certain technical assumptions on the nature of the Gibbs property and underlying random walk measure, we prove that the entire spatial process of the edge is tight. We then apply this black-box theory to the log-gamma polymer Gibbsian line ensemble which we construct. The edge of this line ensemble is the transversal free energy process for the polymer, and our theorem implies tightness with the ubiquitous KPZ class $2/3$ exponent, as well as Brownian absolute continuity of all the subsequential limits. A key technical innovation which fuels our general result is the construction of a continuous grand monotone coupling of Gibbsian line ensembles with respect to their boundary data (entrance and exit values, and bounding curves). {\em Continuous} means that the Gibbs measure varies continuously with respect to varying the boundary data, {\em grand} means that all uncountably many boundary data measures are coupled to the same probability space, and {\em monotone} means that raising the values of the boundary data likewise raises the associated measure. This result applies to a general class of Gibbsian line ensembles where the underlying random walk measure is discrete time, continuous valued and log-convex, and the interaction Hamiltonian is nearest neighbor and convex.

preprint2022arXiv

The ASEP speed process

For ASEP with step initial data and a second class particle started at the origin we prove that as time goes to infinity the second class particle almost surely achieves a velocity that is uniformly distributed on $[-1,1]$. This positively resolves Conjecture 1.9 and 1.10 of [Amir, Angel and Valko, "The TASEP speed process", Annals of Probability 39, 1205--1242, 2011] and allows us to construct the ASEP speed process.

preprint2021arXiv

Stochastic PDE limit of the dynamic ASEP

We study a stochastic PDE limit of the height function of the dynamic asymmetric simple exclusion process (dynamic ASEP). A degeneration of the stochastic Interaction Round-a-Face (IRF) model of arXiv:1701.05239, dynamic ASEP has a jump parameter $q\in (0,1)$ and a dynamical parameter $α>0$. It degenerates to the standard ASEP height function when $α$ goes to $0$ or $\infty$. We consider very weakly asymmetric scaling, i.e., for $\varepsilon$ tending to zero we set $q=e^{-\varepsilon}$ and look at fluctuations, space and time in the scales $\varepsilon^{-1}$, $\varepsilon^{-2}$ and $\varepsilon^{-4}$. We show that under such scaling the height function of the dynamic ASEP converges to the solution of the space-time Ornstein-Uhlenbeck process. We also introduce the dynamic ASEP on a ring with generalized rate functions. Under the very weakly asymmetric scaling, we show that the dynamic ASEP (with generalized jump rates) on a ring also converges to the solution of the space-time Ornstein-Uhlenbeck process on $[0,1]$ with periodic boundary conditions.

preprint2020arXiv

Half-space Macdonald processes

Macdonald processes are measures on sequences of integer partitions built using the Cauchy summation identity for Macdonald symmetric functions. These measures are a useful tool to uncover the integrability of many probabilistic systems, including the Kardar-Parisi-Zhang (KPZ) equation and a number of other models in its universality class. In this paper we develop the structural theory behind half-space variants of these models and the corresponding half-space Macdonald processes. These processes are built using a Littlewood summation identity instead of the Cauchy identity, and their analysis is considerably harder than their full-space counterparts. We compute moments and Laplace transforms of observables for general half-space Macdonald measures. Introducing new dynamics preserving this class of measures, we relate them to various stochastic processes, in particular the log-gamma polymer in a half-quadrant (they are also related to the stochastic six-vertex model in a half-quadrant and the half-space ASEP). For the polymer model, we provide explicit integral formulas for the Laplace transform of the partition function. Non-rigorous saddle point asymptotics yield convergence of the directed polymer free energy to either the Tracy-Widom GOE, GSE or the Gaussian distribution depending on the average size of weights on the boundary.

preprint2020arXiv

Invariance of polymer partition functions under the geometric RSK correspondence

We prove that the values of discrete directed polymer partition functions involving multiple non-intersecting paths remain invariant under replacing the background weights by their images under the geometric RSK correspondence. This result is inspired by a recent and remarkable identity proved by Dauvergne, Orthmann and Virág which is recovered as the zero-temperature, semi-discrete limit of our main result.

preprint2020arXiv

KPZ equation correlations in time

We consider the narrow wedge solution to the Kardar-Parisi-Zhang stochastic PDE under the characteristic $3:2:1$ scaling of time, space and fluctuations. We study the correlation of fluctuations at two different times. We show that when the times are close to each other, the correlation approaches one at a power-law rate with exponent $2/3$, while when the two times are remote from each other, the correlation tends to zero at a power-law rate with exponent $-1/3$. We also prove exponential-type tail bounds for differences of the solution at two space-time points. Three main tools are pivotal to proving these results: 1) a representation for the two-time distribution in terms of two independent narrow wedge solutions; 2) the Brownian Gibbs property of the KPZ line ensemble; and 3) recently proved one-point tail bounds on the narrow wedge solution.

preprint2020arXiv

KPZ equation tails for general initial data

We consider the upper and lower tail probabilities for the centered (by time$/24$) and scaled (according to KPZ time$^{1/3}$ scaling) one-point distribution of the Cole-Hopf solution of the KPZ equation when started with initial data drawn from a very general class. For the lower tail, we prove an upper bound which demonstrates a crossover from super-exponential decay with exponent $3$ in the shallow tail to an exponent $5/2$ in the deep tail. For the upper tail, we prove super-exponential decay bounds with exponent $3/2$ at all depth in the tail.

preprint2020arXiv

KPZ line ensemble

For each $t\geq 1$ we construct an $\mathbf{N}$-indexed ensemble of random continuous curves with three properties: 1. The lowest indexed curve is distributed as the time $t$ Hopf-Cole solution to the Kardar-Parisi-Zhang (KPZ) stochastic PDE with narrow wedge initial data; 2. The entire ensemble satisfies a resampling invariance which we call the $\mathbf{H}$-Brownian Gibbs property (with $\mathbf{H}(x)=e^{x}$); 3. Increments of the lowest indexed curve, when centered by $-t/24$ and scaled down vertically by $t^{1/3}$ and horizontally by $t^{2/3}$, remain uniformly absolutely continuous (i.e. have tight Radon-Nikodym derivatives) with respect to Brownian bridges as time $t$ goes to infinity. This construction uses as inputs the diffusion that O'Connell discovered in relation to the O'Connell-Yor semi-discrete Brownian polymer, the convergence result of Nica of the lowest indexed curve of that diffusion to the solution of the KPZ equation with narrow wedge initial data, and the one-point distribution formula proved by Amir-Corwin-Quastel for the solution of the KPZ equation with narrow wedge initial data. We provide four main applications of this construction: 1. Uniform (as $t$ goes to infinity) Brownian absolute continuity of the time $t$ solution to the KPZ equation with narrow wedge initial data, even when scaled vertically by $t^{1/3}$ and horizontally by $t^{2/3}$; 2. Universality of the $t^{1/3}$ one-point (vertical) fluctuation scale for the solution of the KPZ equation with general initial data; 3. Concentration in the $t^{2/3}$ scale for the endpoint of the continuum directed random polymer; 4.Exponential upper and lower tail bounds for the solution at fixed time of the KPZ equation with general initial data.

preprint2020arXiv

SPDE Limit of Weakly Inhomogeneous ASEP

We study ASEP in a spatially inhomogeneous environment on a torus $ \mathbb{T}^{(N)} = \mathbb{Z}/N\mathbb{Z} $ of $ N $ sites. A given inhomogeneity $ \widetilde{a}(x) \in (0,\infty) $, $x \in \mathbb{T} $, perturbs the overall asymmetric jumping rates $ r < \ell \in (0,1) $ at bonds, so that particles jump from site $ x $ to $ x+1 $ with rate $ r\widetilde{a}(x) $ and from $x+1$ to $x$ with rate $ \ell\widetilde{a}(x) $ (subject to the exclusion rule in both cases). Under the limit $ N\to\infty $, we suitably tune the asymmetry $ (\ell-r) $ to zero like $N^{-\frac{1}{2}}$ and the inhomogeneity $ \widetilde{a}(x) $ to unity, so that the two compete on equal footing. At the level of the Gärtner (or microscopic Hopf--Cole) transform, we show convergence to a new SPDE -- the Stochastic Heat Equation with a mix of spatial and spacetime multiplicative noise. Equivalently, at the level of the height function we show convergence to the Kardar--Parisi--Zhang equation with a mix of spatial and spacetime additive noise. Our method applies to a general class of $ \widetilde{a}(x) $, which, in particular, includes i.i.d., long-range correlated, and periodic inhomogeneities. The key technical component of our analysis consists of a host of estimates on the \emph{kernel} of the semigroup $ \mathcal{Q}(t):=e^{t\mathcal{H}} $ for a Hill-type operator $ \mathcal{H}:= \frac12\partial_{xx} + \mathcal{A}'(x) $, and its discrete analog, where $ \mathcal{A} $ (and its discrete analog) is a generic Hölder continuous function.

preprint2020arXiv

Stochastic six-vertex model in a half-quadrant and half-line open ASEP

We consider the asymmetric simple exclusion process (ASEP) on the positive integers with an open boundary condition. We show that, when starting devoid of particles and for a certain boundary condition, the height function at the origin fluctuates asymptotically (in large time $τ$) according to the Tracy-Widom GOE distribution on the $τ^{1/3}$ scale. This is the first example of KPZ asymptotics for a half-space system outside the class of free-fermionic/determinantal/Pfaffian models. Our main tool in this analysis is a new class of probability measures on Young diagrams that we call half-space Macdonald processes, as well as two surprising relations. The first relates a special (Hall-Littlewood) case of these measures to the half-space stochastic six-vertex model (which further limits to ASEP) using a Yang-Baxter graphical argument. The second relates certain averages under these measures to their half-space (or Pfaffian) Schur process analogs via a refined Littlewood identity.

preprint2019arXiv

Limit Shape of Subpartition Maximizing Partitions

This is an expository note answering a question posed to us by Richard Stanley, in which we prove a limit shape theorem for partitions of $n$ which maximize the number of subpartitions. The limit shape and the growth rate of the number of subpartitions are explicit. The key ideas are to use large deviations estimates for random walks, together with convex analysis and the Hardy-Ramanujan asymptotics. Our limit shape coincides with Vershik's limit shape for uniform random partitions.

preprint2015arXiv

Directed random polymers via nested contour integrals

We study the partition function of two versions of the continuum directed polymer in 1+1 dimension. In the full-space version, the polymer starts at the origin and is free to move transversally in the reals, and in the half-space version, the polymer starts at the origin but is reflected at the origin and stays in the negative reals. The partition functions solves the stochastic heat equation in full-space or half-space with mixed boundary condition at the origin; or equivalently the free energy satisfies the Kardar-Parisi-Zhang equation. We derive exact formulas for the Laplace transforms of the partition functions. In the full-space this is expressed as a Fredholm determinant while in the half-space this is expressed as a Fredholm Pfaffian. Taking long-time asymptotics we show that the limiting free energy fluctuations scale with exponent 1/3 and are given by the GUE and GSE Tracy-Widom distributions. These formulas come from summing divergent moment generating functions, hence are not mathematically justified. The primary purpose of this work is to present a mathematical perspective on the polymer replica method which is used to derive these results. In contrast to other replica method work, we do not appeal directly to the Bethe ansatz for the Lieb-Liniger model but rather utilize nested contour integral formulas for moments as well as their residue expansions.

preprint2015arXiv

Fluctuations of TASEP and LPP with general initial data

We prove Airy process variational formulas for the one-point probability distribution of (discrete time parallel update) TASEP with general initial data, as well as last passage percolation from a general lattice path to a point. We also consider variants of last passage percolation with inhomogeneous parameter geometric weights and provide variational formulas of a similar nature. This proves one aspect of the conjectural description of the renormalization fixed point of the Kardar-Parisi-Zhang universality class.

preprint2015arXiv

Stationary measure of the driven two-dimensional q-Whittaker particle system on the torus

We consider a q-deformed version of the uniform Gibbs measure on dimers on the periodized hexagonal lattice (equivalently, on interlacing particle configurations, if vertical dimers are seen as particles) and show that it is invariant under a certain irreversible q-Whittaker dynamic. Thereby we provide a new non-trivial example of driven interacting two-dimensional particle system, or of (2+1)-dimensional stochastic growth model, with explicit stationary measure. We emphasize that this measure is far from being a product Bernoulli measure. These Gibbs measures and dynamics both arose earlier in the theory of Macdonald processes. The q=0 degeneration of the Gibbs measures reduce to the usual uniform dimer measures with given tilt, the degeneration of the dynamics originate in the study of Schur processes and the degeneration of the results contained herein were recently treated in work of the second author.

preprint2015arXiv

The q-PushASEP: A New Integrable Model for Traffic in 1+1 Dimension

We introduce a new interacting (stochastic) particle system q-PushASEP which interpolates between the q-TASEP introduced by Borodin and Corwin (see arXiv:1111.4408, and also arXiv:1207.5035; arXiv:1305.2972; arXiv:1212.6716) and the q-PushTASEP introduced recently by Borodin and Petrov (arXiv:1305.5501). In the q-PushASEP, particles can jump to the left or to the right, and there is a certain partially asymmetric pushing mechanism present. This particle system has a nice interpretation as a model of traffic on a one-lane highway in which cars are able to accelerate or slow down. Using the quantum many body system approach, we explicitly compute the expectations of a large family of observables for this system in terms of nested contour integrals. We also discuss relevant Fredholm determinantal formulas for the distribution of the location of each particle, and connections of the model with a certain two-sided version of Macdonald processes and with the semi-discrete stochastic heat equation.

preprint2015arXiv

The strict-weak lattice polymer

We introduce the strict-weak polymer model, and show the KPZ universality of the free energy fluctuation of this model for a certain range of parameters. Our proof relies on the observation that the discrete time geometric q-TASEP model, studied earlier by A. Borodin and I. Corwin, scales to this polymer model in the limit q->1. This allows us to exploit the exact results for geometric q-TASEP to derive a Fredholm determinant formula for the strict-weak polymer, and in turn perform rigorous asymptotic analysis to show KPZ scaling and GUE Tracy-Widom limit for the free energy fluctuations. We also derive moments formulae for the polymer partition function directly by Bethe ansatz, and identify the limit of the free energy using a stationary version of the polymer model.

preprint2014arXiv

From duality to determinants for q-TASEP and ASEP

We prove duality relations for two interacting particle systems: the $q$-deformed totally asymmetric simple exclusion process ($q$-TASEP) and the asymmetric simple exclusion process (ASEP). Expectations of the duality functionals correspond to certain joint moments of particle locations or integrated currents, respectively. Duality implies that they solve systems of ODEs. These systems are integrable and for particular step and half-stationary initial data we use a nested contour integral ansatz to provide explicit formulas for the systems' solutions, and hence also the moments. We form Laplace transform-like generating functions of these moments and via residue calculus we compute two different types of Fredholm determinant formulas for such generating functions. For ASEP, the first type of formula is new and readily lends itself to asymptotic analysis (as necessary to reprove GUE Tracy--Widom distribution fluctuations for ASEP), while the second type of formula is recognizable as closely related to Tracy and Widom's ASEP formula [Comm. Math. Phys. 279 (2008) 815--844, J. Stat. Phys. 132 (2008) 291--300, Comm. Math. Phys. 290 (2009) 129--154, J. Stat. Phys. 140 (2010) 619--634]. For $q$-TASEP, both formulas coincide with those computed via Borodin and Corwin's Macdonald processes [Probab. Theory Related Fields (2014) 158 225--400]. Both $q$-TASEP and ASEP have limit transitions to the free energy of the continuum directed polymer, the logarithm of the solution of the stochastic heat equation or the Hopf--Cole solution to the Kardar--Parisi--Zhang equation. Thus, $q$-TASEP and ASEP are integrable discretizations of these continuum objects; the systems of ODEs associated to their dualities are deformed discrete quantum delta Bose gases; and the procedure through which we pass from expectations of their duality functionals to characterizing generating functions is a rigorous version of the replica trick in physics.

preprint2014arXiv

Macdonald processes, quantum integrable systems and the Kardar-Parisi-Zhang universality class

Integrable probability has emerged as an active area of research at the interface of probability/mathematical physics/statistical mechanics on the one hand, and representation theory/integrable systems on the other. Informally, integrable probabilistic systems have two properties: 1) It is possible to write down concise and exact formulas for expectations of a variety of interesting observables (or functions) of the system. 2) Asymptotics of the system and associated exact formulas provide access to exact descriptions of the properties and statistics of large universality classes and universal scaling limits for disordered systems. We focus here on examples of integrable probabilistic systems related to the Kardar-Parisi-Zhang (KPZ) universality class and explain how their integrability stems from connections with symmetric function theory and quantum integrable systems.

preprint2014arXiv

Moments and Lyapunov exponents for the parabolic Anderson model

We study the parabolic Anderson model in $(1+1)$ dimensions with nearest neighbor jumps and space-time white noise (discrete space/continuous time). We prove a contour integral formula for the second moment and compute the second moment Lyapunov exponent. For the model with only jumps to the right, we prove a contour integral formula for all moments and compute moment Lyapunov exponents of all orders.

preprint2014arXiv

Stochastic six-vertex model

We study the asymmetric six-vertex model in the quadrant with parameters on the stochastic line. We show that the random height function of the model converges to an explicit deterministic limit shape as the mesh size tends to 0. We further prove that the one-point fluctuations around the limit shape are asymptotically governed by the GUE Tracy-Widom distribution. We also explain an equivalent formulation of our model as an interacting particle system, which can be viewed as a discrete time generalization of ASEP started from the step initial condition. Our results confirm an earlier prediction of Gwa and Spohn (1992) that this system belongs to the KPZ universality class.

preprint2014arXiv

The $(q,μ,ν)$-Boson process and $(q,μ,ν)$-TASEP

We prove a intertwining relation (or Markov duality) between the $(q,μ,ν)$-Boson process and $(q,μ,ν)$-TASEP, two discrete time Markov chains introduced by Povolotsky. Using this and a variant of the coordinate Bethe ansatz we compute nested contour integral formulas for expectations of a family of observables of the $(q,μ,ν)$-TASEP when started from step initial data. We then utilize these to prove a Fredholm determinant formula for distribution of the location of any given particle.

preprint2013arXiv

Brownian Gibbs property for Airy line ensembles

Consider N Brownian bridges B_i:[-N,N] -> R, B_i(-N) = B_i(N) = 0, 1 <= i <= N, conditioned not to intersect. The edge-scaling limit of this system is obtained by taking a limit as N -> infinity of these curves scaled around (0,2^{1/2} N) horizontally by a factor of N^{2/3} and vertically by N^{1/3}. If a parabola is added to each limit curve, an x-translation invariant process sometimes called the multi-line Airy process is obtained. We prove the existence of a version of this process (which we call the Airy line ensemble) in which curves are a.s. everywhere continuous and non-intersecting. This process naturally arises in the study of growth processes and random matrix ensembles, as do related processes with "wanderers" and "outliers". We formulate our results to treat these relatives as well. Note that the law of the finite collection of Brownian bridges above has the property -- called the Brownian Gibbs property -- of being invariant under the following action. Select an index 1 <= k <= N and erase B_k on a fixed time interval (a,b) subset of (-N,N); then replace this erased curve with a new curve on (a,b) according to the law of a Brownian bridge between the two existing endpoints (a,B_k(a)) and (b,B_k(b)), conditioned to intersect neither the curve above nor the one below. We show that this property is preserved under the edge-scaling limit and thus establish that the Airy line ensemble has the Brownian Gibbs property. An immediate consequence is a proof of M. Prahofer and H. Spohn's prediction that the lines of the Airy line ensemble are locally absolutely continuous with respect to Brownian motion. We also prove the conjecture of K. Johansson that the top line of the Airy line ensemble minus a parabola attains its maximum at a unique point, thus establishing the asymptotic law of the transversal fluctuation of last passage percolation with geometric weights.

preprint2013arXiv

Crossover distributions at the edge of the rarefaction fan

We consider the weakly asymmetric limit of simple exclusion process with drift to the left, starting from step Bernoulli initial data with $ρ_-<ρ_+$ so that macroscopically one has a rarefaction fan. We study the fluctuations of the process observed along slopes in the fan, which are given by the Hopf--Cole solution of the Kardar-Parisi-Zhang (KPZ) equation, with appropriate initial data. For slopes strictly inside the fan, the initial data is a Dirac delta function and the one point distribution functions have been computed in [Comm. Pure Appl. Math. 64 (2011) 466-537] and [Nuclear Phys. B 834 (2010) 523-542]. At the edge of the rarefaction fan, the initial data is one-sided Brownian. We obtain a new family of crossover distributions giving the exact one-point distributions of this process, which converge, as $T\nearrow\infty$ to those of the Airy $\mathcal{A}_{2\to \mathrm{BM}}$ process. As an application, we prove moment and large deviation estimates for the equilibrium Hopf-Cole solution of KPZ. These bounds rely on the apparently new observation that the FKG inequality holds for the stochastic heat equation. Finally, via a Feynman-Kac path integral, the KPZ equation also governs the free energy of the continuum directed polymer, and thus our formula may also be interpreted in those terms.

preprint2013arXiv

Discrete time q-TASEPs

We introduce two new exactly solvable (stochastic) interacting particle systems which are discrete time versions of q-TASEP. We call these geometric and Bernoulli discrete time q-TASEP. We obtain concise formulas for expectations of a large enough class of observables of the systems to completely characterize their fixed time distributions when started from step initial condition. We then extract Fredholm determinant formulas for the marginal distribution of the location of any given particle. Underlying this work is the fact that these expectations solve closed systems of difference equations which can be rewritten as free evolution equations with k-1 two-body boundary conditions -- discrete q-deformed versions of the quantum delta Bose gas. These can be solved via a nested contour integral ansatz. The same solutions also arise in the study of Macdonald processes, and we show how the systems of equations our expectations solve are equivalent to certain commutation relations involving the Macdonald first difference operator.

preprint2013arXiv

Free energy fluctuations for directed polymers in random media in 1+1 dimension

We consider two models for directed polymers in space-time independent random media (the O'Connell-Yor semi-discrete directed polymer and the continuum directed random polymer) at positive temperature and prove their KPZ universality via asymptotic analysis of exact Fredholm determinant formulas for the Laplace transform of their partition functions. In particular, we show that for large time tau, the probability distributions for the free energy fluctuations, when rescaled by tau^{1/3}, converges to the GUE Tracy-Widom distribution. We also consider the effect of boundary perturbations to the quenched random media on the limiting free energy statistics. For the semi-discrete directed polymer, when the drifts of a finite number of the Brownian motions forming the quenched random media are critically tuned, the statistics are instead governed by the limiting Baik-Ben Arous-Peche distributions from spiked random matrix theory. For the continuum polymer, the boundary perturbations correspond to choosing the initial data for the stochastic heat equation from a particular class, and likewise for its logarithm -- the Kardar-Parisi-Zhang equation. The Laplace transform formula we prove can be inverted to give the one-point probability distribution of the solution to these stochastic PDEs for the class of initial data.

preprint2013arXiv

Macdonald processes

Macdonald processes are probability measures on sequences of partitions defined in terms of nonnegative specializations of the Macdonald symmetric functions and two Macdonald parameters q,t in [0,1). We prove several results about these processes, which include the following. (1) We explicitly evaluate expectations of a rich family of observables for these processes. (2) In the case t=0, we find a Fredholm determinant formula for a q-Laplace transform of the distribution of the last part of the Macdonald-random partition. (3) We introduce Markov dynamics that preserve the class of Macdonald processes and lead to new "integrable" 2d and 1d interacting particle systems. (4) In a large time limit transition, and as q goes to 1, the particles of these systems crystallize on a lattice, and fluctuations around the lattice converge to O'Connell's Whittaker process that describe semi-discrete Brownian directed polymers. (5) This yields a Fredholm determinant for the Laplace transform of the polymer partition function, and taking its asymptotics we prove KPZ universality for the polymer (free energy fluctuation exponent 1/3 and Tracy-Widom GUE limit law). (6) Under intermediate disorder scaling, we recover the Laplace transform of the solution of the KPZ equation with narrow wedge initial data. (7) We provide contour integral formulas for a wide array of polymer moments. (8) This results in a new ansatz for solving quantum many body systems such as the delta Bose gas.

preprint2013arXiv

Tropical Combinatorics and Whittaker functions

We establish a fundamental connection between the geometric RSK correspondence and GL(N,R)-Whittaker functions, analogous to the well known relationship between the RSK correspondence and Schur functions. This gives rise to a natural family of measures associated with GL(N,R)-Whittaker functions which are the analogues in this setting of the Schur measures on integer partitions. The corresponding analogue of the Cauchy-Littlewood identity can be seen as a generalisation of an integral identity for GL(N,R)-Whittaker functions due to Bump and Stade. As an application, we obtain an explicit integral formula for the Laplace transform of the law of the partition function associated with a one-dimensional directed polymer model with log-gamma weights recently introduced by one of the authors (TS).

preprint2013arXiv

Two ways to solve ASEP

The purpose of this article is to describe the two approaches to compute exact formulas (which are amenable to asymptotic analysis) for the probability distribution of the current of particles past a given site in the asymmetric simple exclusion process (ASEP) with step initial data. The first approach is via a variant of the coordinate Bethe ansatz and was developed in work of Tracy and Widom in 2008-2009, while the second approach is via a rigorous version of the replica trick and was developed in work of Borodin, Sasamoto and the author in 2012.

preprint2011arXiv

The Kardar-Parisi-Zhang equation and universality class

Brownian motion is a continuum scaling limit for a wide class of random processes, and there has been great success in developing a theory for its properties (such as distribution functions or regularity) and expanding the breadth of its universality class. Over the past twenty five years a new universality class has emerged to describe a host of important physical and probabilistic models (including one dimensional interface growth processes, interacting particle systems and polymers in random environments) which display characteristic, though unusual, scalings and new statistics. This class is called the Kardar-Parisi-Zhang (KPZ) universality class and underlying it is, again, a continuum object -- a non-linear stochastic partial differential equation -- known as the KPZ equation. The purpose of this survey is to explain the context for, as well as the content of a number of mathematical breakthroughs which have culminated in the derivation of the exact formula for the distribution function of the KPZ equation started with {\it narrow wedge} initial data. In particular we emphasize three topics: (1) The approximation of the KPZ equation through the weakly asymmetric simple exclusion process; (2) The derivation of the exact one-point distribution of the solution to the KPZ equation with narrow wedge initial data; (3) Connections with directed polymers in random media. As the purpose of this article is to survey and review, we make precise statements but provide only heuristic arguments with indications of the technical complexities necessary to make such arguments mathematically rigorous.

preprint2011arXiv

Universality of slow decorrelation in KPZ growth

There has been much success in describing the limiting spatial fluctuations of growth models in the Kardar-Parisi-Zhang (KPZ) universality class. A proper rescaling of time should introduce a non-trivial temporal dimension to these limiting fluctuations. In one-dimension, the KPZ class has the dynamical scaling exponent $z=3/2$, that means one should find a universal space-time limiting process under the scaling of time as $t\,T$, space like $t^{2/3} X$ and fluctuations like $t^{1/3}$ as $t\to\infty$. In this paper we provide evidence for this belief. We prove that under certain hypotheses, growth models display temporal slow decorrelation. That is to say that in the scalings above, the limiting spatial process for times $t\, T$ and $t\, T+t^ν$ are identical, for any $ν<1$. The hypotheses are known to be satisfied for certain last passage percolation models, the polynuclear growth model, and the totally / partially asymmetric simple exclusion process. Using slow decorrelation we may extend known fluctuation limit results to space-time regions where correlation functions are unknown. The approach we develop requires the minimal expected hypotheses for slow decorrelation to hold and provides a simple and intuitive proof which applied to a wide variety of models.

preprint2010arXiv

Current fluctuations for TASEP: A proof of the Prähofer--Spohn conjecture

We consider the family of two-sided Bernoulli initial conditions for TASEP which, as the left and right densities ($ρ_-,ρ_+$) are varied, give rise to shock waves and rarefaction fans---the two phenomena which are typical to TASEP. We provide a proof of Conjecture 7.1 of [Progr. Probab. 51 (2002) 185--204] which characterizes the order of and scaling functions for the fluctuations of the height function of two-sided TASEP in terms of the two densities $ρ_-,ρ_+$ and the speed $y$ around which the height is observed. In proving this theorem for TASEP, we also prove a fluctuation theorem for a class of corner growth processes with external sources, or equivalently for the last passage time in a directed last passage percolation model with two-sided boundary conditions: $ρ_-$ and $1-ρ_+$. We provide a complete characterization of the order of and the scaling functions for the fluctuations of this model's last passage time $L(N,M)$ as a function of three parameters: the two boundary/source rates $ρ_-$ and $1-ρ_+$, and the scaling ratio $γ^2=M/N$. The proof of this theorem draws on the results of [Comm. Math. Phys. 265 (2006) 1--44] and extensively on the work of [Ann. Probab. 33 (2005) 1643--1697] on finite rank perturbations of Wishart ensembles in random matrix theory.

preprint2010arXiv

Limit processes for TASEP with shocks and rarefaction fans

We consider the totally asymmetric simple exclusion process (TASEP) with two-sided Bernoulli initial condition, i.e., with left density rho_- and right density rho_+. We consider the associated height function, whose discrete gradient is given by the particle occurrences. Macroscopically one has a deterministic limit shape with a shock or a rarefaction fan depending on the values of rho_{+/-}. We characterize the large time scaling limit of the fluctuations as a function of the densities rho_{+/-} and of the different macroscopic regions. Moreover, using a slow decorrelation phenomena, the results are extended from fixed time to the whole space-time, except along the some directions (the characteristic solutions of the related Burgers equation) where the problem is still open. On the way to proving the results for TASEP, we obtain the limit processes for the fluctuations in a class of corner growth processes with external sources, of equivalently for the last passage time in a directed percolation model with two-sided boundary conditions. Additionally, we provide analogous results for eigenvalues of perturbed complex Wishart (sample covariance) matrices.

preprint2010arXiv

Probability Distribution of the Free Energy of the Continuum Directed Random Polymer in 1+1 dimensions

We consider the solution of the stochastic heat equation \partial_T \mathcal{Z} = 1/2 \partial_X^2 \mathcal{Z} - \mathcal{Z} \dot{\mathscr{W}} with delta function initial condition \mathcal{Z} (T=0)= δ_0 whose logarithm, with appropriate normalizations, is the free energy of the continuum directed polymer, or the solution of the Kardar-Parisi-Zhang equation with narrow wedge initial conditions. We obtain explicit formulas for the one-dimensional marginal distributions -- the {\it crossover distributions} -- which interpolate between a standard Gaussian distribution (small time) and the GUE Tracy-Widom distribution (large time). The proof is via a rigorous steepest descent analysis of the Tracy-Widom formula for the asymmetric simple exclusion with anti-shock initial data, which is shown to converge to the continuum equations in an appropriate weakly asymmetric limit. The limit also describes the crossover behaviour between the symmetric and asymmetric exclusion processes.