Source author record

Mykhaylo Shkolnikov

Mykhaylo Shkolnikov 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

30works
12topics
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

30 published item(s)

preprint2022arXiv

Convergence of a time-stepping scheme to the free boundary in the supercooled Stefan problem

The supercooled Stefan problem and its variants describe the freezing of a supercooled liquid in physics, as well as the large system limits of systemic risk models in finance and of integrate-and-fire models in neuroscience. Adopting the physics terminology, the supercooled Stefan problem is known to feature a finite-time blow-up of the freezing rate for a wide range of initial temperature distributions in the liquid. Such a blow-up can result in a discontinuity of the liquid-solid boundary. In this paper, we prove that the natural Euler time-stepping scheme applied to a probabilistic formulation of the supercooled Stefan problem converges to the liquid-solid boundary of its physical solution globally in time, in the Skorokhod M1 topology. In the course of the proof, we give an explicit bound on the rate of local convergence for the time-stepping scheme. We also run numerical tests to compare our theoretical results to the practically observed convergence behavior.

preprint2022arXiv

Stefan problem with surface tension: global existence of physical solutions under radial symmetry

We consider the Stefan problem with surface tension, also known as the Stefan-Gibbs-Thomson problem, in an ambient space of arbitrary dimension. Assuming the radial symmetry of the initial data we introduce a novel "probabilistic" notion of solution, which can accommodate the discontinuities in time (of the radius) of the evolving aggregate. Our main result establishes the global existence of a probabilistic solution satisfying the natural upper bound on the sizes of the discontinuities. Moreover, we prove that the upper bound is sharp in dimensions d>2, in the sense that none of the discontinuities in the solution can be decreased in magnitude. The detailed analysis of the discontinuities, via appropriate stochastic representations, differentiates this work from the previous literature on weak solutions to the Stefan problem with surface tension.

preprint2021arXiv

Scaling limits of external multi-particle DLA on the plane and the supercooled Stefan problem

We consider (a variant of) the external multi-particle diffusion-limited aggregation (MDLA) process of Rosenstock and Marquardt on the plane. Based on the recent findings of [11], [10] in one space dimension it is natural to conjecture that the scaling limit of the growing aggregate in such a model is given by the growing solid phase in a suitable "probabilistic" formulation of the single-phase supercooled Stefan problem for the heat equation. To address this conjecture, we extend the probabilistic formulation from [10] to multiple space dimensions. We then show that the equation that characterizes the growth rate of the solid phase in the supercooled Stefan problem is satisfied by the scaling limit of the external MDLA process with an inequality, which can be strict in general. In the course of the proof, we establish two additional results interesting in their own right: (i) the stability of a "crossing property" of planar Brownian motion and (ii) a rigorous connection between the probabilistic solutions to the supercooled Stefan problem and its classical and weak solutions.

preprint2020arXiv

Global solutions to the supercooled Stefan problem with blow-ups: regularity and uniqueness

We consider the supercooled Stefan problem, which captures the freezing of a supercooled liquid, in one space dimension. A probabilistic reformulation of the problem allows to define global solutions, even in the presence of blow-ups of the freezing rate. We provide a complete description of such solutions, by relating the temperature distribution in the liquid to the regularity of the ice growth process. The latter is shown to transition between (i) continuous differentiability, (ii) Holder continuity, and (iii) discontinuity. In particular, in the second regime we rediscover the square root behavior of the growth process pointed out by Stefan in his seminal paper [Ste89] from 1889 for the ordinary Stefan problem. In our second main theorem, we establish the uniqueness of the global solutions, a first result of this kind in the context of growth processes with singular self-excitation when blow-ups are present.

preprint2020arXiv

Superposition and mimicking theorems for conditional McKean-Vlasov equations

We consider conditional McKean-Vlasov stochastic differential equations (SDEs), such as the ones arising in the large-system limit of mean field games and particle systems with mean field interactions when common noise is present. The conditional time-marginals of the solutions to these SDEs satisfy non-linear stochastic partial differential equations (SPDEs) of the second order, whereas the laws of the conditional time-marginals follow Fokker-Planck equations on the space of probability measures. We prove two superposition principles: The first establishes that any solution of the SPDE can be lifted to a solution of the conditional McKean-Vlasov SDE, and the second guarantees that any solution of the Fokker-Planck equation on the space of probability measures can be lifted to a solution of the SPDE. We use these results to obtain a mimicking theorem which shows that the conditional time-marginals of an Ito process can be emulated by those of a solution to a conditional McKean-Vlasov SDE with Markovian coefficients. This yields, in particular, a tool for converting open-loop controls into Markovian ones in the context of controlled McKean-Vlasov dynamics.

preprint2020arXiv

Zero kinetic undercooling limit in the supercooled Stefan problem

We study the solutions of the one-phase supercooled Stefan problem with kinetic undercooling, which describes the freezing of a supercooled liquid, in one spatial dimension. Assuming that the initial temperature lies between the equilibrium freezing point and the characteristic invariant temperature throughout the liquid our main theorem shows that, as the kinetic undercooling parameter tends to zero, the free boundary converges to the (possibly irregular) free boundary in the supercooled Stefan problem without kinetic undercooling, whose uniqueness has been recently established in [DNS19], [LS18b]. The key tools in the proof are a Feynman-Kac formula, which expresses the free boundary in the problem with kinetic undercooling through a local time of a reflected process, and a resulting comparison principle for the free boundaries with different kinetic undercooling parameters.

preprint2016arXiv

Intertwinings of beta-Dyson Brownian motions of different dimensions

We show that for all positive beta the semigroups of beta-Dyson Brownian motions of different dimensions are intertwined. The proof relates beta-Dyson Brownian motions directly to Jack symmetric polynomials and omits an approximation of the former by discrete space Markov chains, thereby disposing of the technical assumption beta>1 in [GS]. The corresponding results for beta-Dyson Ornstein-Uhlenbeck processes are also presented.

preprint2016arXiv

SPDE limit of the global fluctuations in rank-based models

We consider systems of diffusion processes ("particles") interacting through their ranks (also referred to as "rank-based models" in the mathematical finance literature). We show that, as the number of particles becomes large, the process of fluctuations of the empirical cumulative distribution functions converges to the solution of a linear parabolic SPDE with additive noise. The coefficients in the limiting SPDE are determined by the hydrodynamic limit of the particle system which, in turn, can be described by the porous medium PDE. The result opens the door to a thorough investigation of large equity markets and investment therein. In the course of the proof we also derive quantitative propagation of chaos estimates for the particle system.

preprint2016arXiv

Stochastic Airy semigroup through tridiagonal matrices

We determine the operator limit for large powers of random tridiagonal matrices as the size of the matrix grows. The result provides a novel expression in terms of functionals of Brownian motions for the Laplace transform of the Airy$_β$ process, which describes the largest eigenvalues in the $β$ ensembles of random matrix theory. Another consequence is a Feynman-Kac formula for the stochastic Airy operator of Ramírez, Rider, and Virág. As a side result, we find that the difference between the area underneath a standard Brownian excursion and one half of the integral of its squared local times is a Gaussian random variable.

preprint2015arXiv

A construction of infinite Brownian particle systems

The paper identifies families of quasi-stationary initial conditions for infinite Brownian particle systems within a large class and provides a construction of the particle systems themselves started from such initial conditions. Examples of particle systems falling into our framework include Brownian versions of TASEP-like processes such as the diffusive scaling limit of the q-TASEP process. In this context the spacings between consecutive particles form infinite-dimensional versions of the softly reflected Brownian motions recently introduced in the finite-dimensional setting by O'Connell and Ortmann and are of independent interest. The proof of the main result is based on intertwining relations satisfied by the particle systems involved which can be regarded as infinite-dimensional analogues of the suitably generalized Burke's Theorem.

preprint2015arXiv

Asymptotic analysis of forward performance processes in incomplete markets and their ill-posed HJB equations

We consider the problem of optimal portfolio selection under forward investment performance criteria in an incomplete market. The dynamics of the prices of the traded assets depend on a pair of stochastic factors, namely, a slow factor (e.g. a macroeconomic indicator) and a fast factor (e.g. stochastic volatility). We analyze the associated forward performance SPDE and provide explicit formulae for the leading order and first order correction terms for the forward investment process and the optimal feedback portfolios. They both depend on the investor's initial preferences and the dynamically changing investment opportunities. The leading order terms resemble their time-monotone counterparts, but with the appropriate stochastic time changes resulting from averaging phenomena. The first-order terms compile the reaction of the investor to both the changes in the market input and his recent performance. Our analysis is based on an expansion of the underlying ill-posed HJB equation, and it is justified by means of an appropriate remainder estimate.

preprint2015arXiv

Multidimensional sticky Brownian motions as limits of exclusion processes

We study exclusion processes on the integer lattice in which particles change their velocities due to stickiness. Specifically, whenever two or more particles occupy adjacent sites, they stick together for an extended period of time, and the entire particle system is slowed down until the ``collision'' is resolved. We show that under diffusive scaling of space and time such processes converge to what one might refer to as a sticky reflected Brownian motion in the wedge. The latter behaves as a Brownian motion with constant drift vector and diffusion matrix in the interior of the wedge, and reflects at the boundary of the wedge after spending an instant of time there. In particular, this leads to a natural multidimensional generalization of sticky Brownian motion on the half-line, which is of interest in both queuing theory and stochastic portfolio theory. For instance, this can model a market, which experiences a slowdown due to a major event (such as a court trial between some of the largest firms in the market) deciding about the new market leader.

preprint2014arXiv

Interacting particle systems at the edge of multilevel Dyson Brownian motions

We study the joint asymptotic behavior of spacings between particles at the edge of multilevel Dyson Brownian motions, when the number of levels tends to infinity. Despite the global interactions between particles in multilevel Dyson Brownian motions, we observe a decoupling phenomenon in the limit: the global interactions become negligible and only the local interactions remain. The resulting limiting objects are interacting particle systems which can be described as Brownian versions of certain totally asymmetric exclusion processes. This is the first appearance of a particle system with local interactions in the context of general $β$ random matrix models.

preprint2014arXiv

Multilevel Dyson Brownian motions via Jack polynomials

We introduce multilevel versions of Dyson Brownian motions of arbitrary parameter $β>0$, generalizing the interlacing reflected Brownian motions of Warren for $β=2$. Such processes unify $β$ corners processes and Dyson Brownian motions in a single object. Our approach is based on the approximation by certain multilevel discrete Markov chains of independent interest, which are defined by means of Jack symmetric polynomials. In particular, this approach allows to show that the levels in a multilevel Dyson Brownian motion are intertwined (at least for $β\ge 1$) and to give the corresponding link explicitly.

preprint2013arXiv

Large deviations for interacting Bessel-like processes and applications to systemic risk

We establish a process level large deviation principle for systems of interacting Bessel-like diffusion processes. By establishing weak uniqueness for the limiting non-local SDE of McKean-Vlasov type, we conclude that the latter describes the process level hydrodynamic limit of such systems and obtain a propagation of chaos result. This is the first instance of results of this type in the context of interacting diffusion processes under explicit assumptions on the coefficients, where the diffusion coefficients are allowed to be both non-Lipschitz and degenerate. In the second part of the paper, we explain how systems of this type naturally arise in the study of stability of the interbank lending system and describe some financial implications of our results.

preprint2013arXiv

Time-reversal of reflected Brownian motions in the orthant

We determine the processes obtained from a large class of reflected Brownian motions (RBMs) in the nonnegative orthant by means of time reversal. The class of RBMs we deal with includes, but is not limited to, RBMs in the so-called Harrison-Reiman class [4] having diagonal covariance matrices. For such RBMs our main result resolves the long-standing open problem of determining the time reversal of RBMs beyond the skew-symmetric case treated by R.J. Williams in [16]. In general, the time-reversed process itself is no longer a RBM, but its distribution is absolutely continuous with respect to a certain auxiliary RBM. In the course of the proof we introduce a novel discrete approximation scheme for the class of RBMs described above, and use it to determine the semigroups dual to the semigroups of such RBMs.

preprint2012arXiv

Fluctuations of Martingales and Winning Probabilities of Game Contestants

Within a contest there is some probability M_i(t) that contestant i will be the winner, given information available at time t, and M_i(t) must be a martingale in t. Assume continuous paths, to capture the idea that relevant information is acquired slowly. Provided each contestant's initial winning probability is at most b, one can easily calculate, without needing further model specification, the expectations of the random variables N_b = number of contestants whose winning probability ever exceeds b, and D_{ab} = total number of downcrossings of the martingales over an interval [a,b]. The distributions of N_b and D_{ab} do depend on further model details, and we study how concentrated or spread out the distributions can be. The extremal models for N_b correspond to two contrasting intuitively natural methods for determining a winner: progressively shorten a list of remaining candidates, or sequentially examine candidates to be declared winner or eliminated. We give less precise bounds on the variability of D_{ab}. We formalize the setting of infinitely many contestants each with infinitesimally small chance of winning, in which the explicit results are more elegant. A canonical process in this setting is the Wright-Fisher diffusion associated with an infinite population of initially distinct alleles; we show how this process fits our setting and raise the problem of finding the distributions of N_b and D_{ab} for this process.

preprint2012arXiv

Limits of Multilevel TASEP and similar processes

We study the asymptotic behavior of a class of stochastic dynamics on interlacing particle configurations (also known as Gelfand-Tsetlin patterns). Examples of such dynamics include, in particular, a multi-layer extension of TASEP and particle dynamics related to the shuffling algorithm for domino tilings of the Aztec diamond. We prove that the process of reflected interlacing Brownian motions introduced by Warren in \cite{W} serves as a universal scaling limit for such dynamics.

preprint2012arXiv

On a non-linear transformation between Brownian martingales

The paper studies a non-linear transformation between Brownian martingales, which is given by the inverse of the pricing operator in the mathematical finance terminology. Subsequently, the solvability of systems of equations corresponding to such transformations is investigated. The latter give rise to novel monotone pathwise couplings of an arbitrary number of certain diffusion processes with varying diffusion coefficients. In the case that there is an uncountable number of these diffusion processes and that the index set is an interval such couplings can be viewed as models for the growth of one-dimensional random surfaces. With this motivation in mind, we derive the appropriate stochastic partial differential equations for the growth of such surfaces.

preprint2012arXiv

Small time central limit theorems for semimartingales with applications

We give conditions under which the normalized marginal distribution of a semimartingale converges to a Gaussian limit law as time tends to zero. In particular, our result is applicable to solutions of stochastic differential equations with locally bounded and continuous coefficients. The limit theorems are subsequently extended to functional central limit theorems on the process level. We present two applications of the results in the field of mathematical finance: to the pricing of at-the-money digital options with short maturities and short time implied volatility skews.

preprint2012arXiv

Some universal estimates for reversible Markov chains

We obtain universal estimates on the convergence to equilibrium and the times of coupling for continuous time irreducible reversible finite-state Markov chains, both in the total variation and in the L^2 norms. The estimates in total variation norm are obtained using a novel identity relating the convergence to equilibrium of a reversible Markov chain to the increase in the entropy of its one-dimensional distributions. In addition, we propose a universal way of defining the ultrametric partition structure on the state space of such Markov chains. Finally, for chains reversible with respect to the uniform measure, we show how the global convergence to equilibrium can be controlled using the entropy accumulated by the chain.

preprint2012arXiv

Systems of Brownian particles with asymmetric collisions

We study systems of Brownian particles on the real line, which interact by splitting the local times of collisions among themselves in an asymmetric manner. We prove the strong existence and uniqueness of such processes and identify them with the collections of ordered processes in a Brownian particle system, in which the drift coefficients, the diffusion coefficients, and the collision local times for the individual particles are assigned according to their ranks. These Brownian systems can be viewed as generalizations of those arising in first-order models for equity markets in the context of stochastic portfolio theory, and are able to correct for several shortcomings of such models while being equally amenable to computations. We also show that, in addition to being of interest in their own right, such systems of Brownian particles arise as universal scaling limits of systems of jump processes on the integer lattice with local interactions. A key step in the proof is the analysis of a generalization of Skorokhod maps which include `local times' at the intersection of faces of the nonnegative orthant. The result extends the convergence of TASEP to its continuous analogue. Finally, we identify those among the Brownian particle systems which have a probabilistic structure of determinantal type.

preprint2012arXiv

Two Models of Stochastic Loss Given Default

We propose two structural models for stochastic losses given default which allow to model the credit losses of a portfolio of defaultable financial instruments. The credit losses are integrated into a structural model of default events accounting for correlations between the default events and the associated losses. We show how the models can be calibrated and analyze the impact of correlations between the occurrences of defaults and recoveries by testing our models for a representative sample portfolio.

preprint2011arXiv

Competing particle systems evolving by interacting Lévy processes

We consider finite and infinite systems of particles on the real line and half-line evolving in continuous time. Hereby, the particles are driven by i.i.d. Lévy processes endowed with rank-dependent drift and diffusion coefficients. In the finite systems we show that the processes of gaps in the respective particle configurations possess unique invariant distributions and prove the convergence of the gap processes to the latter in the total variation distance, assuming a bound on the jumps of the Lévy processes. In the infinite case we show that the gap process of the particle system on the half-line is tight for appropriate initial conditions and same drift and diffusion coefficients for all particles. Applications of such processes include the modeling of capital distributions among the ranked participants in a financial market, the stability of certain stochastic queueing and storage networks and the study of the Sherrington--Kirkpatrick model of spin glasses.

preprint2011arXiv

Convergence rates for rank-based models with applications to portfolio theory

We determine rates of convergence of rank-based interacting diffusions and semimartingale reflecting Brownian motions to equilibrium. Convergence rate for the total variation metric is derived using Lyapunov functions. Sharp fluctuations of additive functionals are obtained using Transportation Cost-Information inequalities for Markov processes. We work out various applications to the rank-based abstract equity markets used in Stochastic Portfolio Theory. For example, we produce quantitative bounds, including constants, for fluctuations of market weights and occupation times of various ranks for individual coordinates. Another important application is the comparison of performance between symmetric functionally generated portfolios and the market portfolio. This produces estimates of probabilities of "beating the market".

preprint2011arXiv

Large volatility-stabilized markets

We investigate the behavior of systems of interacting diffusion processes, known as volatility-stabilized market models in the mathematical finance literature, when the number of diffusions tends to infinity. We show that, after an appropriate rescaling of the time parameter, the empirical measure of the system converges to the solution of a degenerate parabolic partial differential equation. A stochastic representation of the latter in terms of one-dimensional distributions of a time-changed squared Bessel process allows us to give an explicit description of the limit.

preprint2011arXiv

On the one-sided Tanaka equation with drift

We study questions of existence and uniqueness of weak and strong solutions for a one-sided Tanaka equation with constant drift λ. We observe a dichotomy in terms of the values of the drift parameter: for λ\leq 0, there exists a strong solution which is pathwise unique, thus also unique in distribution; whereas for λ>0, the equation has a unique in distribution weak solution, but no strong solution (and not even a weak solution that spends zero time at the origin). We also show that strength and pathwise uniqueness are restored to the equation via suitable "Brownian perturbations".

preprint2011arXiv

Strong solutions of stochastic equations with rank-based coefficients

We study finite and countably infinite systems of stochastic differential equations, in which the drift and diffusion coefficients of each component (particle) are determined by its rank in the vector of all components of the solution. We show that strong existence and uniqueness hold until the first time three particles collide. Motivated by this result, we improve significantly the existing conditions for the absence of such triple collisions in the case of finite-dimensional systems, and provide the first condition of this type for systems with a countable infinity of particles.

preprint2010arXiv

Concentration of measure for systems of Brownian particles interacting through their ranks

We consider a finite or countable collection of one-dimensional Brownian particles whose dynamics at any point in time is determined by their rank in the entire particle system. Using Transportation Cost Inequalities for stochastic processes we provide uniform fluctuation bounds for the ordered particles, their local time of collisions, and various associated statistics over intervals of time. For example, such processes, when exponentiated and rescaled, exhibit power law decay under stationarity; we derive concentration bounds for the empirical estimates of the index of the power law over large intervals of time. A key ingredient in our proofs is a novel upper bound on the Lipschitz constant of the Skorokhod map that transforms a multidimensional Brownian path to a path which is constrained not to leave the positive orthant.

preprint2010arXiv

Large systems of diffusions interacting through their ranks

We study the limiting behaviour of the empirical measure of a system of diffusions interacting through their ranks when the number of diffusions tends to infinity. We prove that the limiting dynamics is given by a McKean-Vlasov evolution equation. Moreover, we show that in a wide range of cases the evolution of the cumulative distribution function under the limiting dynamics is governed by the generalized porous medium equation with convection. The uniqueness theory for the latter is used to establish the uniqueness of solutions of the limiting McKean-Vlasov equation and the law of large numbers for the corresponding systems of interacting diffusions. The implications of the results for rank-based models of capital distributions in financial markets are also explained.