Reflecting Brownian Motion in Two Dimensions: Exact Asymptotics for the Stationary Distribution

Abstract

We consider a two-dimensional semimartingale reflecting Brownian motion (SRBM) in the nonnegative quadrant. The data of the SRBM consists of a two-dimensional drift vector, a 2 ×2 positive definite covariance matrix, and a 2×2 reflection matrix. Assuming the SRBM is positive recurrent, we are interested in tail asymptotic of its marginal stationary distribution along each direction in the quadrant. For a given direction, the marginal tail distribution has the exact asymptotic of the form bx κ exp(−αx) as x goes to infinity, where α and b are positive constants and κ takes one of the values −3/2, −1/2, 0, or 1; both the decay rate α and the power κ can be computed explicitly from the given direction and the SRBM data. A key tool in our proof is a relationship governing the moment generating function of the two-dimensional stationary distribution and two moment generating functions of the associated one-dimensional boundary measures. This relationship allows us to characterize the convergence domain of the two-dimensional moment generating function. For a given direction c, the line in this direction intersects the boundary of the convergence domain at one point, and that point uniquely determines the decay rate α. The one-dimensional moment generating function of the marginal distribution along direction c has a singularity at α. Using analytic extension in complex analysis, we characterize the precise nature of the singularity there. Using that characterization and complex inversion techniques, we obtain the exact asymptotic of the marginal tail distribution.

Conjectures on tail asymptotics of the marginal stationary distribution for a multidimensional SRBM

“…This is equivalent to the fixed point equation (2.12) and (2.13) in [4], and therefore all the conjectures are valid for d = 2. Moreover, exact asymptotics are obtained in [4].…”

“…Furthermore, we may find exact asymptotics including polynomial prefactor to the exponential main term as studied in [4,16]. Our conjectures may suggest how the optimal path looks like in the variational problem for the large deviations principle.…”

“…However, this does not answer what geometrical expression it has. For d = 2, this is answered using extreme points in [4]. For extending this idea, we introduce the following fixed point equations for vectors in R d .…”

“…In this case, J = {1, 2} and A = {1}, {2} for A ∈ 2 J \ {∅, J }. Conjecture 3.1 is proved in [4]. (3.7) can be written as…”

“…A basic idea is the decay rates are determined by the convergence domain of the moment generating function. For d = 2, this approach has been successfully performed by solving a certain fixed point equation in [4]. Similar ideas have been used for the special case of the two dimensional SRBM in [16] and the two dimensional skip free reflecting random walk in [14].…”