The [contact process](https://en.wikipedia.org/wiki/Contact_process_(mathematics)) is a well-studied Markov process. I'm just concerned with the one-dimensional nearest-neighbor version here.


The state space is $\eta\in\{0,1\}^\mathbb Z$, and for state $\eta$ at each site $x\in\mathbb Z$ the dynamics are: $0\rightarrow 1$ at rate $\lambda\sum\_{y:|x-y|=1}\eta(y)$ and $1\rightarrow0$ at rate $1$.


There is a critical rate $\lambda\_c$ such that for $\lambda>\lambda\_c$ we have a positve probability of survival for any initial distribution that does not place unit mass of the state of all zeros. We have convergence to a stationary $\nu\_\lambda$ such that 
$$\rho\_\lambda=\nu\_\lambda\left(\left\{\eta : \eta(x)=1\right\}\right)>0$$
for $\lambda>\lambda\_c$ and is independent of $x.$


What I am interested in is this process restricted to the half-line, $\mathbb Z^+\_0=\{0,1,2,\ldots\}$. In particular, I want to know what the stationary distribution is. It seems like we would have
$$u\_\lambda(x)=\nu\_\lambda\left(\left\{\eta : \eta(x)=1\right\}\right)>0$$
an increasing function of $x$ with $u\_\lambda(0)>0$ and $u\_\lambda(x)\rightarrow\rho\_\lambda$ as $x\rightarrow\infty$.


I can't seem to find any papers that discuss this. Does anyone have any relevant references or maybe can provide some relevant calculations that help to understand $u\_\lambda(x)?$


**Edit:** I'm specifically interested in the questions:  

1) Does $u\_\lambda(x)\rightarrow\rho\_\lambda$ as $x\rightarrow\infty?$  

2) Is $u\_\lambda(x)$ strictly increasing in $x?$  

3) What is the relationship between $u\_\lambda(0)$ and $\displaystyle\lim\_{x\rightarrow\infty}u\_\lambda(x)?$  



I'm fairly certain (1) and (2) are true, but I'd especially like quantitative results for (3). The more specific about $u\_\lambda(x)$ the better, however I'm not asking for anything close to an exact formula as that is likely not possible. 


You can assume I'm familiar with all the basics (monotonicity, coupling, graphical representation, duality, etc., whatever is in Liggett's two books at least).



What is the stationary distribution for the contact process on the half line?
For the first part of your question, you can do this by coupling:
Define two processes: the first is the contact process on $\{0,1,2,\ldots\}$
and the second is the contact process on $\{1,2,3,\ldots\}$ with the same parameters. Call the processes $\eta$ and $\eta'$ respectively.


They can be monotonically coupled: that is: there is a probabaility measure on the stochastic processes $(\eta,\eta')$ such that $\eta(x)\ge \eta'(x)$ for all $x$ such that the evolution of $\eta$ and $\eta'$ follow the correct laws. 


Now $u(x)=\nu(\eta(x)=1)=\nu'(\eta'(x+1)=1)\le \nu(\eta(x+1)=1)=u(x+1)$. 
The inequality here arises from the coupling: on the space where both processes are defined, the even $\{\eta'(x+1)=1\}$ is a subset of $\{\eta(x+1)=1\}$. 


**EDIT:** Some more details here, giving a (in principle) concrete estimate on the difference $u(1)-u(0)$. We'll use the description of the coupling coming from the graphic method. So in the $\eta$ process, there are 'lifelines' above each $x\in\mathbb Z\_0^+$ marked with "kill" symbols (a Poisson process with intensity 1), "left to right" arrows joining a pair of lifelines with intensity $\lambda$ and "right to left arrows" joining a pair of lifelines with intensity $\lambda$.
The $\eta'$ process has the same symbols in the same locations, except there is no lifeline above 0. Now consider the lifeline above 1. If the system is in equilibrium, the system should be in state 1 at each kill symbol with probability $u(0)$. At each kill symbol with no kill symbol in the previous unit of time, the system should be in state 1 with probability at least $u(0)$. Call this an isolated kill symbol. If in addition, there are no arrows from the 2 lifeline to the 1 lifeline for the following 2 units of time, call the kill symbol a *target*. By the Markov property, the $\eta'$ process above $x=1$ is in state 1 immediately below at least $u(0)$ proportion of targets. Now if a target has the property that in the $\eta$ process, there is a right-to-left arrow (from 1 to 0) in the unit of time preceding the target, no kill symbols in the 0 lifeline for 1 unit before and after the target and a left-to-right arrow (from 0 to 1) in the time unit after the target, then $\eta(1)$ is 1 while $\eta'(1)$ is 0 for at least one unit of time. One can then do the calculation to get a concrete estimate of the difference. Since $u(1)\le \rho\_\lambda$, this also gives a lower bound for $\rho\_\lambda-u(0)$. 


