Source author record

Kavita Ramanan

Kavita Ramanan 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

24works
6topics
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

24 published item(s)

preprint2023arXiv

Interacting stochastic processes on sparse random graphs

Large ensembles of stochastically evolving interacting particles describe phenomena in diverse fields including statistical physics, neuroscience, biology, and engineering. In such systems, the infinitesimal evolution of each particle depends only on its own state (or history) and the states (or histories) of neighboring particles with respect to an underlying, possibly random, interaction graph. While these high-dimensional processes are typically too complex to be amenable to exact analysis, their dynamics are quite well understood when the interaction graph is the complete graph. In this case, classical theorems show that in the limit as the number of particles goes to infinity, the dynamics of the empirical measure and the law of a typical particle coincide and can be characterized in terms of a much more tractable dynamical system of reduced dimension called the mean-field limit. In contrast, until recently not much was known about corresponding convergence results in the complementary case when the interaction graph is sparse (i.e., with uniformly bounded average degree). This article provides a brief survey of classical work and then describes recent progress on the sparse regime that relies on a combination of techniques from random graph theory, Markov random fields, and stochastic analysis. The article concludes by discussing ramifications for applications and posing several open problems.

preprint2022arXiv

Local weak convergence for sparse networks of interacting processes

We study the limiting behavior of interacting particle systems indexed by large sparse graphs, which evolve either according to a discrete time Markov chain or a diffusion, in which particles interact directly only with their nearest neighbors in the graph. To encode sparsity we work in the framework of local weak convergence of marked (random) graphs. We show that the joint law of the particle system varies continuously with respect to local weak convergence of the underlying graph marked with the initial conditions. In addition, we show that the global empirical measure converges to a non-random limit for a large class of graph sequences including sparse Erdös-Rényi graphs and configuration models, whereas the empirical measure of the connected component of a uniformly random vertex converges to a random limit. Along the way, we develop some related results on the time-propagation of ergodicity and empirical field convergence, as well as some general results on local weak convergence of Gibbs measures in the uniqueness regime which appear to be new. The results obtained here are also useful for obtaining autonomous descriptions of marginal dynamics of interacting diffusions and Markov chains on sparse graphs. While limits of interacting particle systems on dense graphs have been extensively studied, there are relatively few works that have studied the sparse regime in generality.

preprint2020arXiv

Invariant states of hydrodynamic limits of randomized load balancing networks

Randomized load-balancing algorithms play an important role in improving performance in large-scale networks at relatively low computational cost. A common model of such a system is a network of $N$ parallel queues in which incoming jobs with independent and identically distributed service times are routed on arrival using the join-the-shortest-of-$d$-queues routing algorithm. Under fairly general conditions, it was shown by Aghajani and Ramanan that as $N\rightarrow\infty$, the state dynamics converges to the unique solution of a countable system of coupled deterministic measure-valued equations called the hydrodynamic equations. In this article, a characterization of invariant states of these hydrodynamic equations is obtained and, when $d=2$, used to construct a numerical algorithm to compute the queue length distribution and mean virtual waiting time in the invariant state. Additionally, it is also shown that under a suitable tail condition on the service distribution, the queue length distribution of the invariant state exhibits a doubly exponential tail decay, thus demonstrating a vast improvement in performance over the case $d=1$, which corresponds to random routing, when the tail decay could even be polynomial. Furthermore, numerical evidence is provided to support the conjecture that the invariant state is the limit of the steady-state distributions of the $N$-server models. The proof methodology, which entails analysis of a coupled system of measure-valued equations, can potentially be applied to other many-server systems with general service distributions, where measure-valued representations are useful.

preprint2020arXiv

Large deviations for configurations generated by Gibbs distributions with energy functionals consisting of singular interaction and weakly confining potentials

We establish large deviation principles (LDPs) for empirical measures associated with a sequence of Gibbs distributions on $n$-particle configurations, each of which is defined in terms of an inverse temperature $% β_n$ and an energy functional consisting of a (possibly singular) interaction potential and a (possibly weakly) confining potential. Under fairly general assumptions on the potentials, we use a common framework to establish LDPs both with speeds $β_n/n \rightarrow \infty$, in which case the rate function is expressed in terms of a functional involving the potentials, and with speed $β_n =n$, when the rate function contains an additional entropic term. Such LDPs are motivated by questions arising in random matrix theory, sampling, simulated annealing and asymptotic convex geometry. Our approach, which uses the weak convergence method developed by Dupuis and Ellis, establishes LDPs with respect to stronger Wasserstein-type topologies. Our results address several interesting examples not covered by previous works, including the case of a weakly confining potential, which allows for rate functions with minimizers that do not have compact support, thus resolving several open questions raised in a work of Chafa\"ı et al.

preprint2016arXiv

A Skorokhod Map on Measure-Valued Paths with Applications to Priority Queues

The Skorokhod map on the half-line has proved to be a useful tool for studying processes with non-negativity constraints. In this work we introduce a measure-valued analog of this map that transforms each element $ζ$ of a certain class of càdlàg paths that take values in the space of signed measures on the half-line to a càdlàg path that takes values in the space of non-negative measures on $[0,\infty)$ in such a way that for each $x > 0$, the path $t \mapsto ζ_t[0,x]$ is transformed via a Skorokhod map on the half-line, and the regulating functions for different $x > 0$ are coupled. We establish regularity properties of this map and show that the map provides a convenient tool for studying queueing systems in which tasks are prioritized according to a continuous parameter. Three such well known models are the earliest-deadline-first, the shortest-job-first and the shortest-remaining-processing-time scheduling policies. For these applications, we show how the map provides a unified framework within which to form fluid model equations, prove uniqueness of solutions to these equations and establish convergence of scaled state processes to the fluid model. In particular, for these models, we obtain new convergence results in time-inhomogeneous settings, which appear to fall outside the purview of existing approaches.

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

Large Deviation Principle For Finite-State Mean Field Interacting Particle Systems

We establish a large deviation principle for the empirical measure process associated with a general class of finite-state mean field interacting particle systems with Lipschitz continuous transition rates that satisfy a certain ergodicity condition. The approach is based on a variational representation for functionals of a Poisson random measure. Under an appropriate strengthening of the ergodicity condition, we also prove a locally uniform large deviation principle. The main novelty is that more than one particle is allowed to change its state simultaneously, and so a standard approach to the proof based on a change of measure with respect to a system of independent particles is not possible. The result is shown to be applicable to a wide range of models arising from statistical physics, queueing systems and communication networks. Along the way, we establish a large deviation principle for a class of jump Markov processes on the simplex, whose rates decay to zero as they approach the boundary of the domain. This result may be of independent interest.

preprint2016arXiv

On directional derivatives of Skorokhod maps in convex polyhedral domains

The study of both sensitivity analysis and differentiability of the stochastic flow of a reflected process in a convex polyhedral domain is challenging because the dynamics are discontinuous at the boundary of the domain and the boundary of the domain is not smooth. These difficulties can be addressed by studying directional derivatives of an associated extended Skorokhod map, which is a deterministic mapping that takes an unconstrained path to a suitably reflected version. In this work we develop an axiomatic framework for the analysis of directional derivatives of a large class of Lipschitz continuous extended Skorokhod maps in convex polyhedral domains with oblique directions of reflection. We establish existence of directional derivatives at a path whose reflected version satisfies a certain boundary jitter property, and also show that the right-continuous regularization of such a directional derivative can be characterized as the unique solution to a Skorokhod-type problem, where both the domain and directions of reflection vary (discontinuously) with time. A key ingredient in the proof is establishing certain contraction properties for a family of (oblique) derivative projection operators. As an application, we establish pathwise differentiability of reflected Brownian motion in the nonnegative quadrant with respect to the initial condition, drift vector, dispersion matrix and directions of reflection. The results of this paper are also used in subsequent work to establish pathwise differentiability of a much larger class of reflected diffusions in convex polyhedral domains.

preprint2015arXiv

Cramér's theorem is atypical

The empirical mean of $n$ independent and identically distributed (i.i.d.) random variables $(X_1,\dots,X_n)$ can be viewed as a suitably normalized scalar projection of the $n$-dimensional random vector $X^{(n)}\doteq(X_1,\dots,X_n)$ in the direction of the unit vector $n^{-1/2}(1,1,\dots,1) \in \mathbb{S}^{n-1}$. The large deviation principle (LDP) for such projections as $n\rightarrow\infty$ is given by the classical Cramér's theorem. We prove an LDP for the sequence of normalized scalar projections of $X^{(n)}$ in the direction of a generic unit vector $θ^{(n)} \in \mathbb{S}^{n-1}$, as $n\rightarrow\infty$. This LDP holds under fairly general conditions on the distribution of $X_1$, and for "almost every" sequence of directions $(θ^{(n)})_{n\in\mathbb{N}}$. The associated rate function is "universal" in the sense that it does not depend on the particular sequence of directions. Moreover, under mild additional conditions on the law of $X_1$, we show that the universal rate function differs from the Cramér rate function, thus showing that the sequence of directions $n^{-1/2}(1,1,\dots,1) \in \mathbb{S}^{n-1},$ $n \in \mathbb{N}$, corresponding to Cramér's theorem is atypical.

preprint2015arXiv

Large deviations for random projections of $\ell^p$ balls

Let $p\in[1,\infty]$. Consider the projection of a uniform random vector from a suitably normalized $\ell^p$ ball in $\mathbb{R}^n$ onto an independent random vector from the unit sphere. We show that sequences of such random projections, when suitably normalized, satisfy a large deviation principle (LDP) as the dimension $n$ goes to $\infty$, which can be viewed as an annealed LDP. We also establish a quenched LDP (conditioned on a fixed sequence of projection directions) and show that for $p\in(1,\infty]$ (but not for $p=1$), the corresponding rate function is "universal", in the sense that it coincides for "almost every" sequence of projection directions. We also analyze some exceptional sequences of directions in the "measure zero" set, including the directions corresponding to the classical Cramér's theorem, and show that those directions yield LDPs with rate functions that are distinct from the universal rate function of the quenched LDP. Lastly, we identify a variational formula that relates the annealed and quenched LDPs, and analyze the minimizer of this variational formula. These large deviation results complement the central limit theorem for convex sets, specialized to the case of sequences of $\ell^p$ balls.

preprint2015arXiv

Limits of relative entropies associated with weakly interacting particle systems

The limits of scaled relative entropies between probability distributions associated with N-particle weakly interacting Markov processes are considered. The convergence of such scaled relative entropies is established in various settings. The analysis is motivated by the role relative entropy plays as a Lyapunov function for the (linear) Kolmogorov forward equation associated with an ergodic Markov process, and Lyapunov function properties of these scaling limits with respect to nonlinear finite-state Markov processes are studied in the companion paper [6].

preprint2015arXiv

Local stability of Kolmogorov forward equations for finite state nonlinear Markov processes

The focus of this work is on local stability of a class of nonlinear ordinary differential equations (ODE) that describe limits of empirical measures associated with finite-state weakly interacting N-particle systems. Local Lyapunov functions are identified for several classes of such ODE, including those associated with systems with slow adaptation and Gibbs systems. Using results from [5] and large deviations heuristics, a partial differential equation (PDE) associated with the nonlinear ODE is introduced and it is shown that positive definite subsolutions of this PDE serve as local Lyapunov functions for the ODE. This PDE characterization is used to construct explicit Lyapunov functions for a broad class of models called locally Gibbs systems. This class of models is significantly larger than the family of Gibbs systems and several examples of such systems are presented, including models with nearest neighbor jumps and models with simultaneous jumps that arise in applications.

preprint2015arXiv

Obliquely reflected Brownian motion in non-smooth planar domains

We construct obliquely reflected Brownian motions in all bounded simply connected planar domains, including non-smooth domains, with general reflection vector fields on the boundary. Conformal mappings and excursion theory are our main technical tools. A key intermediate step, which may be of independent interest, is an alternative characterization of reflected Brownian motions in smooth bounded planar domains with a given field of angles of oblique reflection on the boundary in terms of a pair of quantities, namely an integrable positive harmonic function, which represents the stationary distribution of the process, and a real number that represents, in a suitable sense, the asymptotic rate of rotation of the process around a reference point in the domain. Furthermore, we also show that any obliquely reflected Brownian motion in a simply connected Jordan domain can be obtained as a suitable limit of obliquely reflected Brownian motions in smooth domains.

preprint2014arXiv

Large Deviations for Weighted Sums of Stretched Exponential Random Variables

We consider the probability that a weighted sum of $n$ i.i.d. random variables $X_j$, $j = 1, . . ., n$, with stretched exponential tails is larger than its expectation and determine the rate of its decay, under suitable conditions on the weights. We show that the decay is subexponential, and identify the rate function in terms of the tails of $X_j$ and the weights. Our result generalizes the large deviation principle given by Kiesel and Stadtmüller [8] as well as the tail asymptotics for sums of i.i.d. random variables provided by Nagaev [10, 11]. As an application of our result, motivated by random projections of high-dimensional vectors, we consider the case of random, self-normalized weights that are independent of the sequence $\{X_j\}_{j \in \mathbb N}$, identify the decay rate for both the quenched and annealed large deviations in this case, and show that they coincide. As another example we consider weights derived from kernel functions that arise in non-parametric regression.

preprint2014arXiv

On the submartingale problem for reflected diffusions in domains with piecewise smooth boundaries

Two frameworks that have been used to characterize reflected diffusions include stochastic differential equations with reflection and the so-called submartingale problem. We introduce a general formulation of the submartingale problem for (obliquely) reflected diffusions in domains with piecewise C^2 boundaries and piecewise continuous reflection vector fields. Under suitable assumptions, we show that well-posedness of the submartingale problem is equivalent to existence and uniqueness in law of weak solutions to the corresponding stochastic differential equation with reflection. Our result generalizes to the case of reflecting diffusions a classical result due to Stroock and Varadhan on the equivalence of well-posedness of martingale problems and well-posedness of weak solutions of stochastic differential equations in d-dimensional Euclidean space. The analysis in the case of reflected diffusions in domains with non-smooth boundaries is considerably more subtle and requires a careful analysis of the behavior of the reflected diffusion on the boundary of the domain. In particular, the equivalence can fail to hold when our assumptions are not satisfied. The equivalence we establish allows one to transfer results on reflected diffusions characterized by one approach to reflected diffusions analyzed by the other approach. As an application, we provide a characterization of stationary distributions of a large class of reflected diffusions in convex polyhedral domains.

preprint2012arXiv

Asymptotic approximations for stationary distributions of many-server queues with abandonment

A many-server queueing system is considered in which customers arrive according to a renewal process and have service and patience times that are drawn from two independent sequences of independent, identically distributed random variables. Customers enter service in the order of arrival and are assumed to abandon the queue if the waiting time in queue exceeds the patience time. The state of the system with $N$ servers is represented by a four-component process that consists of the forward recurrence time of the arrival process, a pair of measure-valued processes, one that keeps track of the waiting times of customers in queue and the other that keeps track of the amounts of time customers present in the system have been in service and a real-valued process that represents the total number of customers in the system. Under general assumptions, it is shown that the state process is a Feller process, admits a stationary distribution and is ergodic. It is also shown that the associated sequence of scaled stationary distributions is tight, and that any subsequence converges to an invariant state for the fluid limit. In particular, this implies that when the associated fluid limit has a unique invariant state, then the sequence of stationary distributions converges, as $N\rightarrow \infty$, to the invariant state. In addition, a simple example is given to illustrate that, both in the presence and absence of abandonments, the $N\rightarrow \infty$ and $t\rightarrow \infty$ limits cannot always be interchanged.

preprint2012arXiv

Characterization of stationary distributions of reflected diffusions

Given a domain G, a reflection vector field d(.) on the boundary of G, and drift and dispersion coefficients b(.) and σ(.), let L be the usual second-order elliptic operator associated with b(.) and σ(.). Under suitable assumptions that, in particular, ensure that the associated submartingale problem is well posed, it is shown that a probability measure $π$ on \bar{G} is a stationary distribution for the corresponding reflected diffusion if and only if $π(\partial G) = 0$ and $\int_{\bar{G}} L f (x) π(dx) \leq 0$ for every f in a certain class of test functions. Moreover, the assumptions are shown to be satisfied by a large class of reflected diffusions in piecewise smooth multi-dimensional domains with possibly oblique reflection.

preprint2012arXiv

Distributed Parameter Estimation in Sensor Networks: Nonlinear Observation Models and Imperfect Communication

The paper studies distributed static parameter (vector) estimation in sensor networks with nonlinear observation models and noisy inter-sensor communication. It introduces \emph{separably estimable} observation models that generalize the observability condition in linear centralized estimation to nonlinear distributed estimation. It studies two distributed estimation algorithms in separably estimable models, the $\mathcal{NU}$ (with its linear counterpart $\mathcal{LU}$) and the $\mathcal{NLU}$. Their update rule combines a \emph{consensus} step (where each sensor updates the state by weight averaging it with its neighbors' states) and an \emph{innovation} step (where each sensor processes its local current observation.) This makes the three algorithms of the \textit{consensus + innovations} type, very different from traditional consensus. The paper proves consistency (all sensors reach consensus almost surely and converge to the true parameter value,) efficiency, and asymptotic unbiasedness. For $\mathcal{LU}$ and $\mathcal{NU}$, it proves asymptotic normality and provides convergence rate guarantees. The three algorithms are characterized by appropriately chosen decaying weight sequences. Algorithms $\mathcal{LU}$ and $\mathcal{NU}$ are analyzed in the framework of stochastic approximation theory; algorithm $\mathcal{NLU}$ exhibits mixed time-scale behavior and biased perturbations, and its analysis requires a different approach that is developed in the paper.

preprint2011arXiv

Heavy traffic analysis for EDF queues with reneging

This paper presents a heavy-traffic analysis of the behavior of a single-server queue under an Earliest-Deadline-First (EDF) scheduling policy in which customers have deadlines and are served only until their deadlines elapse. The performance of the system is measured by the fraction of reneged work (the residual work lost due to elapsed deadlines) which is shown to be minimized by the EDF policy. The evolution of the lead time distribution of customers in queue is described by a measure-valued process. The heavy traffic limit of this (properly scaled) process is shown to be a deterministic function of the limit of the scaled workload process which, in turn, is identified to be a doubly reflected Brownian motion. This paper complements previous work by Doytchinov, Lehoczky and Shreve on the EDF discipline in which customers are served to completion even after their deadlines elapse. The fraction of reneged work in a heavily loaded system and the fraction of late work in the corresponding system without reneging are compared using explicit formulas based on the heavy traffic approximations. The formulas are validated by simulation results.

preprint2010arXiv

A Dirichlet process characterization of a class of reflected diffusions

For a class of stochastic differential equations with reflection for which a certain ${\mathbb{L}}^p$ continuity condition holds with $p>1$, it is shown that any weak solution that is a strong Markov process can be decomposed into the sum of a local martingale and a continuous, adapted process of zero $p$-variation. When $p=2$, this implies that the reflected diffusion is a Dirichlet process. Two examples are provided to motivate such a characterization. The first example is a class of multidimensional reflected diffusions in polyhedral conical domains that arise as approximations of certain stochastic networks, and the second example is a family of two-dimensional reflected diffusions in curved domains. In both cases, the reflected diffusions are shown to be Dirichlet processes, but not semimartingales.

preprint2010arXiv

Fluid limits of many-server queues with reneging

This work considers a many-server queueing system in which impatient customers with i.i.d., generally distributed service times and i.i.d., generally distributed patience times enter service in the order of arrival and abandon the queue if the time before possible entry into service exceeds the patience time. The dynamics of the system is represented in terms of a pair of measure-valued processes, one that keeps track of the waiting times of the customers in queue and the other that keeps track of the amounts of time each customer being served has been in service. Under mild assumptions, essentially only requiring that the service and reneging distributions have densities, as both the arrival rate and the number of servers go to infinity, a law of large numbers (or fluid) limit is established for this pair of processes. The limit is shown to be the unique solution of a coupled pair of deterministic integral equations that admits an explicit representation. In addition, a fluid limit for the virtual waiting time process is also established. This paper extends previous work by Kaspi and Ramanan, which analyzed the model in the absence of reneging. A strong motivation for understanding performance in the presence of reneging arises from models of call centers.

preprint2010arXiv

SPDE Limits of Many Server Queues

A many-server queueing system is considered in which customers with independent and identically distributed service times enter service in the order of arrival. The state of the system is represented by a process that describes the total number of customers in the system, as well as a measure-valued process that keeps track of the ages of customers in service, leading to a Markovian description of the dynamics. Under suitable assumptions, a functional central limit theorem is established for the sequence of (centered and scaled) state processes as the number of servers goes to infinity. The limit process describing the total number in system is shown to be an Ito diffusion with a constant diffusion coefficient that is insensitive to the service distribution. The limit of the sequence of (centered and scaled) age processes is shown to be a Hilbert space valued diffusion that can also be characterized as the unique solution of a stochastic partial differential equation that is coupled with the Ito diffusion. Furthermore, the limit processes are shown to be semimartingales and to possess a strong Markov property.

preprint2010arXiv

The multi-state hard core model on a regular tree

The classical hard core model from statistical physics, with activity $λ> 0$ and capacity $C=1$, on a graph $G$, concerns a probability measure on the set ${\mathcal I}(G)$ of independent sets of $G$, with the measure of each independent set $I \in {\mathcal I}(G)$ being proportional to $λ^{|I|}$. Ramanan et al. proposed a generalization of the hard core model as an idealized model of multicasting in communication networks. In this generalization, the {\em multi-state} hard core model, the capacity $C$ is allowed to be a positive integer, and a configuration in the model is an assignment of states from $\{0,\ldots,C\}$ to $V(G)$ (the set of nodes of $G$) subject to the constraint that the states of adjacent nodes may not sum to more than $C$. The activity associated to state $i$ is $λ^{i}$, so that the probability of a configuration $σ:V(G)\rightarrow \{0,\ldots, C\}$ is proportional to $λ^{\sum_{v \in V(G)} σ(v)}$. In this work, we consider this generalization when $G$ is an infinite rooted $b$-ary tree and prove rigorously some of the conjectures made by Ramanan et al. In particular, we show that the $C=2$ model exhibits a (first-order) phase transition at a larger value of $λ$ than the $C=1$ model exhibits its (second-order) phase transition. In addition, for large $b$ we identify a short interval of values for $λ$ above which the model exhibits phase co-existence and below which there is phase uniqueness. For odd $C$, this transition occurs in the region of $λ= (e/b)^{1/\ceil{C/2}}$, while for even $C$, it occurs around $λ=(\log b/b(C+2))^{2/(C+2)}$. In the latter case, the transition is first-order.

preprint2007arXiv

An explicit formula for the Skorokhod map on $[0,a]$

The Skorokhod map is a convenient tool for constructing solutions to stochastic differential equations with reflecting boundary conditions. In this work, an explicit formula for the Skorokhod map $Γ_{0,a}$ on $[0,a]$ for any $a>0$ is derived. Specifically, it is shown that on the space $\mathcal{D}[0,\infty)$ of right-continuous functions with left limits taking values in $\mathbb{R}$, $Γ_{0,a}=Λ_a\circ Γ_0$, where $Λ_a:\mathcal{D}[0,\infty)\to\mathcal{D}[0,\infty)$ is defined by \[Λ_a(ϕ)(t)=ϕ(t)-\sup_{s\in[0,t]}\biggl[\bigl(\ phi(s)-a\bigr)^+\wedge\inf_{u\in[s,t]}ϕ(u)\biggr]\] and $Γ_0:\mathcal{D}[0,\infty)\to\mathcal{D}[0,\infty)$ is the Skorokhod map on $[0,\infty)$, which is given explicitly by \[Γ_0(ψ)(t)=ψ(t)+\sup_{s\in[0,t]}[-ψ(s)]^+.\] In addition, properties of $Λ_a$ are developed and comparison properties of $Γ_{0,a}$ are established.