In [ ]:

    

%matplotlib inline
import matplotlib as mpl
import matplotlib.pyplot as plt
import numpy as np

from rmgpy.kinetics.arrhenius import Arrhenius


    

In [ ]:

    

kinetics = [
    Arrhenius(A=(3.95e+16,'s^-1'), n=-0.98, Ea=(285,'kJ/mol'), T0=(298,'K')),
    Arrhenius(A=(1.98e+11,'s^-1'), n=0, Ea=(192,'kJ/mol'), T0=(1,'K')),
    Arrhenius(A=(8.59e+14,'s^-1'), n=0, Ea=(273,'kJ/mol'), T0=(1,'K')),
    Arrhenius(A=(1e+13,'s^-1'), n=0, Ea=(226,'kJ/mol'), T0=(1,'K')),
]


    

In [ ]:

    

pressure = 1e5 # Pa
temperature = np.linspace(298, 2000, 20)


    

In [ ]:

    

plt.style.use('seaborn-talk')

fig = plt.figure()

for i, rate in enumerate(kinetics):
    # Evaluate kinetics
    k = []
    for t in temperature:
        k.append(rate.getRateCoefficient(t, pressure))
    
    x = 1000 / temperature
    y = np.log10(k)
    plt.plot(x, y)
    
plt.xlabel('1000/T (K)')
plt.ylabel('log(k)')
plt.legend(range(len(kinetics)))


    

In [ ]: