Generic two-phase coexistence in a type-2 Schloegl model for autocatalysis on a square lattice: Analysis via heterogeneous master equations
- Ames Lab., Ames, IA (United States); Iowa State Univ., Ames, IA (United States)
- Ames Lab., Ames, IA (United States)
Schloegl's second model (also known as the quadratic contact process) on a square lattice involves spontaneous annihilation of particles at lattice sites at rate p, and their autocatalytic creation at unoccupied sites with n ≥ 2 occupied neighbors at rate kn. Here, kinetic Monte Carlo (KMC) simulation reveals that these models exhibit a nonequilibrium discontinuous phase transition with generic two-phase coexistence: the p value for equistability of coexisting populated and vacuum states, peq(S), depends on the orientation or slope, S, of a planar interface separating those phases. The vacuum state displaces the populated state for p > eq(S), and the opposite applies for p < peq(S) for 0 < S < ∞. The special “combinatorial” rate choice kn = n(n–1)/12 facilitates an appealing simplification of the exact master equations for the evolution of spatially heterogeneous states in the model, which aids analytic investigation of these equations via hierarchical truncation approximations. Truncation produces coupled sets of lattice differential equations which can describe orientation-dependent interface propagation and equistability. The pair approximation predicts that peq(max) = peq(S = 1) = 0.09645 and peq(min) = peq (S → ∞) = 0.08827, values deviating less than 15% from KMC predictions. In the pair approximation, a perfect vertical interface is stationary for all p < peq(S = ∞) = 0.08907, a value exceeding peq(S → ∞). One can regard an interface for large S → ∞ as a vertical interface decorated with isolated kinks. For p < peq(S = ∞), the kink can move in either direction along this otherwise stationary interface depending upon p, but for p = peq(min) the kink is also stationary.
- Research Organization:
- Ames Laboratory (AMES), Ames, IA (United States)
- Sponsoring Organization:
- USDOE Office of Science (SC), Basic Energy Sciences (BES). Chemical Sciences, Geosciences & Biosciences Division
- Grant/Contract Number:
- AC02-07CH11358
- OSTI ID:
- 1960546
- Report Number(s):
- IS-J 11,008; TRN: US2313228
- Journal Information:
- Physical Review. E, Vol. 107, Issue 3; ISSN 2470-0045
- Publisher:
- American Physical Society (APS)Copyright Statement
- Country of Publication:
- United States
- Language:
- English
Similar Records
Discontinuous non-equilibrium phase transition in a threshold Schloegl model for autocatalysis: Generic two-phase coexistence and metastability
Discontinuous non-equilibrium phase transition in a threshold Schloegl model for autocatalysis: Generic two-phase coexistence and metastability