solid phase diffusivity as a function of Li+ concentration

[swhan22]

Hi PyBaMM team,

Just curious, is there a method to make solid phase diffusivity as a function of stoichiometry? For example, can I use 'measured_graphite_diffusivity_Ecker2015.csv' as a lookup table to determine diffusivity? I tried the following code but the solution never converges. I would appreciate any advice. I am new to Python. Thank you.

def D1pos(sto, T):
Tref = 25 + 273.15
E_D1pos = 30e3
arrhenius = pybamm.exp(E_D1pos / pybamm.constants.R * (1 / Tref - 1 / T))
D1pos_lookup = pybamm.Interpolant(D1posLookUp[:, 0], D1posLookUp[:, 1], sto, interpolator = 'cubic')
return D1pos_lookup * arrhenius

[swhan22]

I am trying to utilize solid diffusivity function estimated from GITT. Could anybody tell me whether the above function is capable in pybamm? Thanks so much.