Nonequilibrium steady-state dynamics of Markov processes on graphs

We propose an analytic approach for the steady-state dynamics of Markov processes on locally tree-like graphs. It is based on time-translation invariant probability distributions for edge trajectories, which we encode in terms of infinite matrix products. For homogeneous ensembles on regular graphs, the distribution is parametrized by a single d× d× r^2 <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"> <mml:mrow> <mml:mi>d</mml:mi> <mml:mo>×</mml:mo> <mml:mi>d</mml:mi> <mml:mo>×</mml:mo> <mml:msup> <mml:mi>r</mml:mi> <mml:mn>2</mml:mn> </mml:msup> </mml:mrow> </mml:math> tensor, where r <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"> <mml:mi>r</mml:mi> </mml:math> is the number of states per variable, and d <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"> <mml:mi>d</mml:mi> </mml:math> is the matrix-product bond dimension. While the method becomes exact in the large- d <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"> <mml:mi>d</mml:mi> </mml:math> limit, it typically provides highly accurate results even for small bond dimensions d <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"> <mml:mi>d</mml:mi> </mml:math> . The d^2r^2 <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"> <mml:mrow> <mml:msup> <mml:mi>d</mml:mi> <mml:mn>2</mml:mn> </mml:msup> <mml:msup> <mml:mi>r</mml:mi> <mml:mn>2</mml:mn> </mml:msup> </mml:mrow> </mml:math> parameters are determined by solving a fixed point equation, for which we provide an efficient belief-propagation procedure. We apply this approach to a variety of models, including Ising-Glauber dynamics with symmetric and asymmetric couplings, as well as the SIS model. Even for small d <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"> <mml:mi>d</mml:mi> </mml:math> , the results are compatible with Monte Carlo estimates and accurately reproduce known exact solutions. The method provides access to precise temporal correlations, which, in some regimes, would be virtually impossible to estimate by sampling.