Schlogl’s model

class pints.toy.stochastic.SchloglModel(initial_molecule_count=20)

Schlogl’s system of chemical reactions has a single type of molecules and starts with an initial count \(A(0)\). The evolution of the molecule count is defined through the rates \(k_1\), \(k_2\), \(k_3\) and \(k_4\) and the following equations:

..math::

2A xrightarrow{k_1} 3A 3A xrightarrow{k_2} 2A 0 xrightarrow{k_3} A A xrightarrow{k_4} 0

Extends pints.ForwardModel, pints.toy.ToyModel.

Parameters:

initial_molecule_count – The initial molecule count \(A(0)\).

interpolate_mol_counts(time, mol_count, output_times)

Takes raw times and inputs and mol counts and outputs interpolated values at output_times

n_outputs()

Returns the number of outputs this model has. The default is 1.

n_parameters()

See pints.ForwardModel.n_parameters().

simulate(parameters, times)

See pints.ForwardModel.simulate().

simulate_raw(rates, max_time)

Returns raw times, mol counts when reactions occur.

suggested_parameters()

See pints.toy.ToyModel.suggested_parameters().

suggested_times()

See “meth:pints.toy.ToyModel.suggested_times().