We investigate the relation between the stationary probability distribution of chemical reaction systems and the convective field derived from the chemical Fokker-Planck equation (CFPE) . This convective field takes into account the drift term of the CFPE and the reaction bias introduced by the diffusion term. For one-dimensional systems, fixed points and bifurcations of the convective field correspond to extrema and phenomenological bifurcations of the stationary probability distribution as long as the CFPE is a good approximation to the stochastic dynamics. This provides an efficient way to calculate the effect of system size on the number and location of probability maxima and their phenomenological bifurcations in parameter space. For two-dimensional systems, we study models that have saddle-node and Hopf bifurcations in the macroscopic limit. Here, the existence of two stable fixed points of the convective field correlates either with two peaks of the stationary probability distribution, or with a peak and a shoulder. In contrast, the system-size dependence of the parameter values for a Hopf bifurcation of the convective field is not reflected by the onset of a crater-shaped probability distribution , as the expected crater tends to be washed out for smaller system sizes .
We investigate the relation between the stationary probability distribution of chemical reaction systems and the convective field derived from the chemical Fokker-Planck equation (CFPE) by comparing predictions of the convective field to the results of stochastic simulations based on Gillespie's algorithm. The convective field takes into account the drift term of the CFPE and the reaction bias introduced by the diffusion term. For one-dimensional systems, fixed points and bifurcations of the convective field correspond to extrema and phenomenological bifurcations of the stationary probability distribution whenever the CFPE is a good approximation to the stochastic dynamics. This provides an efficient way to calculate the effect of system size on the number and location of probability maxima and their phenomenological bifurcations in parameter space. For two-dimensional systems, we study models that have saddle-node and Hopf bifurcations in the macroscopic limit. Here, the existence of two stable fixed points of the convective field correlates either with two peaks of the stationary probability distribution, or with a peak and a shoulder. In contrast, a Hopf bifurcation that occurs in the convective field for decreasing system size is not accompanied by the onset of a crater-shaped probability distribution ; decreasing system size rather destroys craters and replaces them by local maxima .
R
. This
by comparing predictions of the convective field to the results of stochastic simulations based on Gillespie's algorithm. The
R
as long as
whenever
D
the system-size dependence of the parameter values for
R
of
that occurs in
R
is not reflected
for decreasing system size is not accompanied
R
, as the expected crater tends to be washed out for smaller system sizes
; decreasing system size rather destroys craters and replaces them by local maxima