origin6.0磁性拟合过程教程

origin6.0磁性拟合过程教程

origin软件的磁性拟合过程

应小木虫及晶体解析联合网站上部分虫友恳求，将origin软件的拟合过程做成视频文件，遗憾的是视频文件太大上传不了网站，故以截图的方式给各位同行粗略演示拟合过程。

1.先将数据导入origin软件，并建立温度（T）,磁化率（Xm）和磁化率温度乘积（XmT）各列。如先右击空白栏勾点add new column，出现空白栏，在空白栏的顶端右击勾点set column values,出现对话框（如下），在“Col(X)=”右边填写：col(温度栏)*相应式量/外场强度/样品质量，如col(b)*15.68/5000/0.021，――OK

2.挑选A 及D栏（ctrl键选择）再点plot—scatter做图

3.建立拟合公式（该例子为双核锰(S1=S2=5/2)

(A)双核锰(S1=S2=5/2)的公式：

见Molecular Magnetism---Olivier Hahn 114页

Ng222 e2x10 e6x28 e12x60 e20x110 e30xd 2x6x12x20x30xkT13 e5 e7 e9 e11 eOrigin 拟合过程：对

Mn2拟合mT曲线

符号：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(30*p2/x);

y=x*0.375*p1*p1*(A/B)

(B)考虑分子场近似：

Origin 拟合过程：对

Mn2拟合mT曲线

符号：p1---g；p2---J/kB；p3---zJ/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(30*p2/x); s=0.375*p1*p1*(A/B); y=s/(1-2*p3*s/(0.375*p1*p1))

(C)考虑顺磁杂质及温度无关项：

Origin 拟合过程：对

Mn2拟合mT曲线

符号：p1---g；p2---J/kB；p3---zJ/kB；p4---m’；p5---TIP

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

4.拷贝拟合式子（若同时考虑分子场近似、顺磁杂质及温度无关项）其式子如下：

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

（注意每行后边“；”号）

5.打开XmT-T曲线origin软件界面，点analysis-non linear curve fit 出来如下界面

6.先建立函数，（如果已建有该函数点More再选择）

若需创建新函数即点New 按钮，

Name,Form(选Y-script),Number of Param 个数，Definition（此栏就将从别处拷来的拟合式子粘贴），如下填写

Save保存，再点Basic Mode进入上级界面，

再点accept

7.选择数据列，点Select Dataset

预先记住温度及磁化率温度乘积两列的栏名（column name）,如A及D

再点对话框其中---->x indep ??? [1…0] 再从available 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,一直拟合直到拟合值及拟合诊断因子值(Chi^2)不变.如果Chi^2很小，同时拟合值也很合乎情理，估计就差不多了，如果二者任何一值不合理，就必须再行分析，是初值设的不好（重设），有些量没有纳入拟合（同时考虑），有些量不适合同时拟合（牺牲一些量另拟合另一些量），模型不对（重建模型）等。

如果觉得拟合到最佳程度，就点Done.

附注：该实例拟合还没完善（Chi^2诊断值偏大）。

谢谢大家

倏忽于2010年12月12日