根据上面的这个公式,我使用gromacs中的gmx current 输出了dsp文件。里面包含了公式中dy/dt中的y函数,文件给出了y的单位是S/(ms·ps),按照文献上说的,只需要拟合数据,斜率就是电导率.
我查到了和我比较类型的一篇文献,和我使用相同量级和相同方法。
这是他的数据:
我把文献上的图片导入到origin里面进行拟合,结果拟合出来是1.78686E-12,但是文献上0.5M拟合出来表中是1.710.
我的疑问是:是不是从y(单位是:s·m-1·ps-1)对t(ps)求导,再求极限得到s/m有量级上的变化,而我没有乘呢?
我试着推导这一过程单位的变化,但是没有正确求出来。
谢谢各位老师!
|