Is there any way where we can calculate parameter covariance matrix instead of correlation matrix in python lmfit package

Question

I have re-parameterize Arrhenius equation of form k = kref*exp(-E/R((1/T)-(1/Tref)) and i wanted to estimate the values of parameters E & kref which i