配合物磁性拟合过程(共10页)

1、精选优质文档-倾情为你奉上网上原有的拟合过程：回应小木虫及晶体解析联合网站上部分虫友恳求，将origin软件的拟合过程做成视频文件，遗憾的是视频文件太大上传不了网站，故以截图的方式给各位同行粗略演示拟合过程。1.先将数据导入origin软件，并建立温度（T）,磁化率（Xm）和磁化率温度乘积（XmT）各列。如先右击空白栏勾点add new column，出现空白栏，在空白栏的顶端右击勾点set column values,出现对话框（如下），在“Col(X)=”右边填写：col(温度栏)*相应式量/外场强度/样品质量，如col(b)*1589.68/5000/0.02641，OK2.挑选A 及D

2、栏（ctrl键选择）再点plotscatter做图3.建立拟合公式（该例子为双核锰(S1=S2=5/2)(A)双核锰(S1=S2=5/2)的公式：见Molecular Magnetism-Olivier Hahn 114页Origin 拟合过程：对Mn2拟合cmT曲线符号：p1-g; p2-J/kB A=2*exp(2*p2/x)+10*exp(6*p2/x)+28*exp(12*p2/x)+60*exp(20*p2/x)+110*exp(30*p2/x);B=1+3*exp(2*p2/x)+5*exp(6*p2/x)+7*exp(12*p2/x)+9*exp(20*p2/x)+11*exp(

3、30*p2/x);y=x*0.375*p1*p1*(A/B)(B)考虑分子场近似：Origin 拟合过程：对Mn2拟合cmT曲线符号：p1-g；p2-J/kB；p3-zJ/kBA=2*exp(2*p2/x)+10*exp(6*p2/x)+28*exp(12*p2/x)+60*exp(20*p2/x)+110*exp(30*p2/x);B=1+3*exp(2*p2/x)+5*exp(6*p2/x)+7*exp(12*p2/x)+9*exp(20*p2/x)+11*exp(30*p2/x);s=0.375*p1*p1*(A/B);y=s/(1-2*p3*s/(0.375*p1*p1)(C)考虑顺磁

4、杂质及温度无关项：Origin 拟合过程：对Mn2拟合cmT曲线符号：p1-g；p2-J/kB；p3-zJ/kB；p4-cm；p5-cTIPA=2*exp(2*p2/x)+10*exp(6*p2/x)+28*exp(12*p2/x)+60*exp(20*p2/x)+110*exp(30*p2/x);B=1+3*exp(2*p2/x)+5*exp(6*p2/x)+7*exp(12*p2/x)+9*exp(20*p2/x)+11*exp(30*p2/x);s=0.375*p1*p1*(A/B)/x;t=s/(1-2*p3*s/(0.375*p1*p1)*(1-p4)+0.125*2.5*3.5/x

5、*p4+p5;y=t*x4.拷贝拟合式子（若同时考虑分子场近似、顺磁杂质及温度无关项）其式子如下：A=2*exp(2*p2/x)+10*exp(6*p2/x)+28*exp(12*p2/x)+60*exp(20*p2/x)+110*exp(30*p2/x);B=1+3*exp(2*p2/x)+5*exp(6*p2/x)+7*exp(12*p2/x)+9*exp(20*p2/x)+11*exp(30*p2/x);s=0.375*p1*p1*(A/B)/x;t=s/(1-2*p3*s/(0.375*p1*p1)*(1-p4)+0.125*2.5*3.5/x*p4+p5;y=t*x（注意每行后边“；

6、”号）5.打开XmT-T曲线origin软件界面，点analysis-non linear curve fit出来如下界面6.先建立函数，（如果已建有该函数点More再选择）若需创建新函数即点New 按钮，Name,Form(选Y-script),Number of Param 个数，Definition（此栏就将从别处拷来的拟合式子粘贴），如下填写Save保存，再点Basic Mode进入上级界面，再点accept7.选择数据列，点Select Dataset预先记住温度及磁化率温度乘积两列的栏名（column name）,如A及D再点对话框其中-x indep ? 10再从availabl

7、e datasets中选后边是a的那一行，如a9129mt300k_a再点Assign,则自变量x的值就变为A列了，同时把因变量y的值变为D列。当然可以改动温度范围，改动=Row=两侧的值就行。8.好，现在就开始拟合了，点Start Fitting记住P1、 P2、 P3、 P4-各代表什么量，估计其值并填入value栏中，有些量不想参与拟合或定值的，可以在Value值改为0或某个定值，Vary处的勾去掉；此处P1-g,先给个初值2.0；P2-J为二核锰簇内的交换作用从曲线知很弱，可以设个小的负值，如-1.0; P3-zJ二核锰簇间的磁交换作用。如果无法同时拟合就先固定一些量来拟合。再点10 lter,一直拟合直到拟合值及拟合诊断因子值(Chi2)不变.如果Chi2很小，同时拟合值也很合乎情理，估计就差不多了，如果二者任何一值不合理，就必须再行分析，是初值设的不好（重设），有些量没有纳入