where r f,j and r r,j are the rates of reaction j in the forward and reverse direction, → k f,j and ← k r,j the respective kinetic constants, and y i refers to θ i or P i , depending on the phase of species i. The initial conditions assume an empty surface exposed to the gas mixture defined by the user: