CN105675711A

CN105675711A - 一种批量获取磁性材料m-t、δs-t曲线的方法

Info

Publication number: CN105675711A
Application number: CN201610008914.1A
Authority: CN
Inventors: 尹海宏; 杨静; 王志亮; 宋长青; 史敏
Original assignee: Nantong University
Current assignee: Ganzhou Xinli Electronics Co ltd
Priority date: 2016-01-05
Filing date: 2016-01-05
Publication date: 2016-06-15
Anticipated expiration: 2036-01-05
Also published as: CN105675711B

Abstract

一种批量获取磁性材料M-T、ΔS-T曲线的方法，建议分类号为G01N？27/72。本发明既不是对现有热力学理论的改进，也不是对现有磁熵变的计算机计算方法作出的算法改进，而是专门针对利用PPMS系统进行磁性材料的测量得到磁矩—磁场—温度测试数据并进行技术处理而得到测量磁性材料的磁熵变-温度特性的技术方案，不仅解决了现有技术中PPMS测量得到的数据难以评估磁性材料的磁熵变—温度特性的难题，其效率大大提升，极大地降低了对计算机内存的需求。

Description

一种批量获取磁性材料M-T、ΔS-T曲线的方法

技术领域

本发明涉及一种磁熵变测量方法，具体地说是一种专门适用于PPMS的铁磁性材料的磁熵变测量方法。

背景技术

本发明的建议分类号为G01N27/72。

本发明既不是对现有热力学理论的改进，也不是对现有磁熵变的计算机计算方法作出的算法改进，而是专门针对利用PPMS系统进行磁性材料的测量得到磁矩—磁场—温度测试数据并进行技术处理而得到测量磁性材料的磁熵变-温度特性的技术方案。

为便于理解本发明，以下对现有技术及现有技术作扼要的介绍。

早在1907年郎杰斐(P.Langevin)就注意到：顺磁体绝热去磁过程中，其温度会降低。从机理上说，固体磁性物质在受磁场作用磁化时，系统的磁有序度加强(磁熵减小)，对外放出热量；再将其去磁，则磁有序度下降(磁熵增大)，又要从外界吸收热量。这种磁性系统在磁场施加与除去过程中所出现的热现象称为磁热效应，又称为磁卡效应。1927年德贝(Debye)和杰克(Giauque)预言了可以利用此效应制冷。1933年杰克实现了绝热去磁制冷。从此，在极低温领域(mK级至16K范围)磁制冷发挥了很大作用。利用磁卡效应可制成固态磁制冷器。

磁熵是磁矩混乱程度的量度，以S_m表示。磁熵S_m是温度T和外加磁场H的函数。磁矩排列愈混乱即无序度愈大，磁熵就愈大。为了研制高效率的磁制冷器，必须研究各种铁磁性材料的磁熵变。

为便捷、可靠地、自动化地完成相关测试，本领域技术人员常常采用PPMS来完成测量过程。PPMS——PhysicalPropertyMeasurementSystem(综合物性测量系统)是美国QuantumDesign公司开发的一款材料的各种物理性质进行基础研究的仪器。PPMS系统的设计理念是在一个完美控制的极低温(50mK)和强磁场(±16T)平台上，集成全自动的磁学、电学、热学和形貌，甚至铁电和介电等各种物性测量手段。一个PPMS系统由基本系统和各种测量和拓展功能选件构成：基本系统提供低温和强磁场的环境，以及整个系统的软硬件控制中心。

根据热力学中Maxwell方程可很容易地得到：当外加磁场从H₀变化到H₁时，磁性材料在这一过程中的磁熵变计算公式为：

Δ S = S (H_{1}, T) - S (H_{0}, T) = {&Integral;}_{H_{0}}^{H_{1}} {(\frac{\partial M (H, T)}{\partial T})}_{H} d H - - - (1)

上式中，S为磁性材料的熵，H为外加磁场，T为热力学温度，M为磁化强度，下标0、1分别代表磁场变化前、后的状态。

通常可以通过如下途径测量得到磁性材料的磁熵变：

(A)在给定温度下给磁性材料施加一随时间变化的外磁场H，H以等步长增加或减小，测量在所述温度下的第一条M-H曲线(也可以称为M-H磁滞回线)；上述测试完成后，使温度增加或减少一固定值，即温度等步长地增大或减小，重复进行测量步骤(A)，得到若干条在不同温度下的M-H曲线；在每次测量过程中温度保持恒定不变；通过上述测量步骤，可以得到所测磁性材料的M(H，T)的三维数组；上述过程可以概括为：测试不同温度下的M-H曲线，即固定温度扫描磁场，每个温度设置的外磁场的序列要一致，这样才比较容易得到M(H，T)的三维数组；

(B)根据公式

Δ S = S (H_{1}, T) - S (H_{0}, T) = {&Integral;}_{H_{0}}^{H_{1}} {(\frac{\partial M (H, T)}{\partial T})}_{H} d H - - - (2)

计算得到磁熵变。

而在测试过程中得到的是一系列离散的测试数据，进行具体数值计算时，通常采用如下两种步骤来计算：

(一)按照如下步骤来进行计算：

(1)将测试得到的按温度组合(M-H)_T(每列数据T相等)的实测数据转换为按磁场强度重新组合的数据(M-T)_H(每列数据H相等)；

(2)对上一步重组的每列数据M对温度T求偏导；

(3)然后先微分再积分运算得到所测材料的磁熵变。

由于数据量巨大，上述步骤通常利用软件来完成，例如Matlab、Origin、Excel等商业软件；

为了减小计算过程中的误差，通常将外磁场变化的步长、温度变化的步长设置小一些，使得数据点足够密集，M-H曲线足够平滑。

(二)根据在测试过程中得到的是一系列离散的测试数据，进行具体计算时，公式(1)改写为采用如下公式来进行计算：

- Δ S = \underset{i}{Σ} \frac{1}{T_{i + 1} - T_{i}} (M_{i} - M_{i + 1}) {ΔH}_{i} - - - (3)

在上式中，i表示测试过程中第i次测试；利用公式(3)也可直接计算得到指定参数条件下的磁性材料的磁熵变。

如果上述外磁场的序列不一致，则根本无法得到M(H，T)的三维数组，也就无从根据上述方法(一)和(二)来计算磁熵变。

在利用PPMS完成相关测量并导出测量数据时，通常得到的数据，每一条M-H曲线中，①：M的变化步长并不是固定值，而是具有一定随机性的值，且②：每一条M-H曲线中的数据点数量也是变化不定的，而且，③：由于所导出的记录数值中，温度值也是由温度传感器所测量得到的实测值，是带有轻微波动的温度值，而不是测量前程序设定的温度序列(确定值)，例如某次测量得到的部分数据为：

Temperature(K),MagneticField(Oe),Moment(emu)

10.0140295,8825.666,4.77920531340975

10.0114298,9997.748,5.35753168743939

10.0107646,11606.298,6.1806593824613

……

9.967885,5850.0365,3.42836713837576

9.9699335,4077.1785,2.44207219325932

9.9713125,2630.387,1.65372870264995

……

12.0178871,8729.4175,0.546468356803494

12.0119467,10598.8385,4.62828899150997

12.0118856,12381.02,5.66277344040436

……

此外，④：而在不同的设定温度值下测量得到的M-H数据点的数量则为3340，3268，3421，3319，3486，3511，3403，……

由此可见，根据PPMS测试导出的实验数据根本无法得到为了计算磁熵变所需的所测磁性材料的M(H，T)的三维数组。

由于上述几点原因，导致无法使用前述计算方法来根据PPMS系统测得的实验数据来得到所测材料的磁熵变。目前可通过人工辅助以其他数据处理程序如Excel、Origin等来完成数据的三维数组，然而，即使借助于这些计算机程序，处理过程往往也需要几个星期甚至几个月的时间，工作繁杂、枯燥且极易出错。

利用方法(一)进行处理时，涉及的变量有T、H、M、ΔS，在计算机数值处理时这些变量构成了一个五维数组，通常测试得到T的数量、计算的ΔS在100左右，H、M、的数量通常在10000-100000之间，因此这样的一个五维数组的数据量通常在10¹⁶-10²⁰之间，若每个数值采用为双精度存储、运算，那么该五维数值所需要的内存(RAM)空间大约在10¹⁰-10¹⁷MB之间，即使降低测量过程中H的数量，即增大H的间隔，使H、M、的数量控制在100-1000之间，该五维数组所需要的内存(RAM)空间也大约在10⁴-10¹¹MB之间，再进一步降低H的数量，尽管可以再降低内存的需求量，但其计算精确度也会大大降低，将导致结果严重不可信。由此可见，利用方法(一)进行数据处理时，其内存需求量极大，通常需要大型计算机来进行处理。这就严重制约了其应用场合。

另外，在研究磁性材料的磁卡效应时，除了分析磁性材料的M-H曲线，科研人员往往还需要分析材料的M-T曲线、曲线、ΔS-T曲线，参见说明书附图1、2；利用方法(二)只能计算得到指定参数条件下的磁性材料的磁熵变，而不能得到M-T曲线、曲线、ΔS-T曲线，使得对磁性材料的分析研究大打折扣。

综上所述，现有技术中各种方法都存在着重大缺点，且无法应用于PPMS测量的磁性材料的数据处理。

发明内容

为了解决现有技术问题，本发明的目的在于提供一种专门适用于PPMS的铁磁性材料的磁学性质测量方法，可批量获得M-T特性曲线、特性曲线、ΔS-T特性曲线。

实现本发明目的具体技术方案是：

A)按照符合PPMS操作规范的流程，称取一定质量的将待测试的磁性材料，记录下该质量值，将待测试的磁性材料装入PPMS测量系统中；

B)在PPMS的控制计算机上设定测试的温度序列值，各个测试温度成等间距序列(间隔为ΔT)，ΔT可以取1K、2K、3K、4K、5K、……、10K；设定测试的磁场强度，这里所述的磁场强度是指磁滞回线中的磁场强度的绝对值的最大值，磁场从0增加到正向最大值，再逐渐变化到负向最大值，然后再逐渐回到正向最大值；在每一个设定的测试温度下，磁性材料均历经同样的所述磁场扫描，所述T指PPMS系统导出的数据文件中的物理量Temperature；

C)在给定温度下给磁性材料施加一随时间变化的外磁场H，测量在所述温度下的第一条M-H曲线(也可以称为M-H磁滞回线、磁化曲线)；上述测试完成后，使温度增加或减少一固定值，即温度等步长地增大或减小，重复进行上述测量步骤；在每次测量过程中温度保持恒定不变；依次逐次测量磁性材料的M-H曲线，直至完成所设定的所有温度序列下的M-H曲线；所述H指PPMS系统导出的数据文件中的物理量MagneticField，所述M指PPMS系统导出的数据文件中的物理量Moment；

依照PPMS操作规范的流程，导出所测试材料的所有数据，除标题行外，每一行数据中由温度值、所施加的外磁场强度、所测磁性材料的磁矩构成；根据测量数据量的多少，导出的测试数据文件可能是一个或多个数据文件，不同的温度、不同的H值下的测试数据合并存储在所述一个或多个数据文件中，即使导出的数据文件为多个，在每个单独的数据文件中，也包含有不同的温度、不同的H值下的测试数据，且这些H值的间隔(ΔH)是不均匀的；

D)将步骤C)测试得到的实测数据分割为按温度组合(M-H)_T(每列数据T相等，但H间隔不均匀)的实测数据：逐行读取所导出的测试数据，利用计算机按照设定的测试温度进行数据分割，将步骤C)中导出的数据文件分割成各个按照设定测量温度值来独立存储M-H数据的文件，将在每个测试温度下测量得到的(M-H)_T数据分别存储在各自独立的(M-H)_T文件中，这些分割得到的(M-H)_T文件存储在外存储器中；在每个独立的(M-H)_T文件中，只存储在测量得到的从正向最大H值扫描到负向最大H值的过程的数据；当步骤D)中得到多个数据文件时，则对多个数据文件轮次进行逐行读取；具体步骤为：

D-a)从外部存储器(如磁盘等存储器，出于读写效率的考虑，不包括磁带)中逐行读取所导出的测试数据到计算机内存(RAM)或CPU的高速缓存(Cache)中；当步骤C)中得到多个数据文件时，对多个数据文件轮次进行逐行读取；当读取第一行(即标题行)时，暂存第一行(即标题行)的文本内容；

D-b)当逐行读取测试数据文件的第二行至最后一行时，以PPMS数据导出时所采用的数据分隔符为不同物理量的分隔判别标志，以步骤C)中所设定的分隔符来区分开不同物理量(Temperature、MagneticField、Moment等)，将每一行的原始数据存储于一个临时数组中；将所读取到的温度实测值取为测试前设定的温度序列中最接近的设定温度值，并在外部存储器中建立一个包含有该设定温度值为文件名的数据文件，其扩展名可以任意设定，如.txt，.dat，.bin等等；在所建立的数据文件中，以追加的方式写入步骤D-a)中暂存的第一行即标题行的内容；

在上述步骤D-b)中，还包括：

逐行读取第三行、第四行……时，每一次读取到的H(MagneticField)值均与上一次读取到的H(MagneticField)值进行比较，以获取在整个测试过程中所施加的外磁场H(MagneticField)的最大值；当读取到的H(MagneticField)值比上一次读取到的H(MagneticField)值小时，即将上一次读取到的H(MagneticField)值记为H最大值，以一个指定的变量来表示这一最大值；

当读取到外磁场H(MagneticField)的最大值的一行测量数据时，开始向位于外存储器中的、在步骤D-b)中建立的数据文件中以追加的方式逐行记入测量数据，直至读取到下一不同设定温度值或导出数据文件的结尾为止；

当读取到的温度值所对应的设定温度值与上一行的设定温度值不同时，在外存储器中建立新的包含有该新的设定温度值为文件名的数据文件，文件命名规则前后保持一致；

在上述步骤D-b)关于写入数据文件的操作中，温度值既可以以温度实测值即原始记录数据写入，也可以以前述所取的测试前设定的温度序列中最接近的设定温度值写入；每一个物理量的数据均以自定义的分隔符进行隔离；

D-c)在步骤C)中导出的测试数据文件为多个时，继续遍历其他所有导出的测试数据文件，重复执行步骤D-b)；

D-d)关闭上述过程中所有打开的数据文件；

E)以计算机外存储器为数据暂存介质，将步骤D得到的按各个温度T单独存储、且磁场强度H间隔不均匀的(M-H)_T实测数据转换、存储为按指定磁场强度H序列单独存储、且H值间隔均匀的(M-H)_T数据文件：

获得指定等间距H序列下的一系列的(M-T)_H数据文件，构成了(H，M，T)三维数组：根据步骤E)分割得的代表不同温度下测量得到的(M-H)_T数据文件，利用计算机按照指定的H变化步长ΔH(即对指定H序列)进行数据采样，获得指定H序列下的一系列的(M-T)_H数据文件，所述指定H序列；H变化步长过大会导致后续计算过程中误差较大，而H变化步长过小则会导致计算效率低下且并不能切实提高计算精确度；

步骤E)具体步骤如下：

E-a)在外部存储器中建立一个用于存储按指定磁场强度H序列、且H值间隔均匀的(M-H)_T数据文件的文件夹；

指定一系列的H序列值，所述H序列值以恒定步长ΔH递增，ΔH依据所需计算精度而定，通常可选择100Oe-1000Oe；以包含指定H序列的H值为文件名，在外部存储器的上述新建的文件夹中建立一系列代表各自H值的独立的数据文件，其扩展名可以任意设定；

在外部存储器所建立的数据文件中，以追加的方式写入标题行，标题行包含M、T两个物理量，如“Temperature(K),Moment(emu)”；

E-b)获取指定H序列中第一个指定H值下的(M-T)_H测试数据序列，并将所述(M-T)_H测试数据序列以追加的方式写入在步骤E-a)中所建立的对应H值的数据文件中：

依照温度从低到高或从高到低的次序，遍历读取步骤E)中存储在外存储器中的、代表在每个测试温度下测量得到的、各自独立的(M-H)_T数据文件；每个文件的读取方式均为逐行读取；

在逐行读取每一个文件时，均比较所读取的每一行的H值，并将所读取到的H值与所指定的H值求差值后取绝对值，当所得到的绝对值最小时，取该行所在的M值作为指定H值下所对应的M值，将所读取的数据文件所代表的温度值、该M值作为一行数据以追加写入的方式写入在步骤E-a)中所建立的对应H值的数据文件中；

当遍历读取完所有的遍历读取步骤E)中存储在外存储器中的、代表在每个测试温度下测量得到的、各自独立的(M-H)_T数据文件后，关闭步骤E-b1)中所有读、写过的文件，获得指定H序列中第一个指定H值下的(M-T)_H测试数据文件；

E-c)更改指定H值为步骤E-a)中指定的H序列的第二、第三、……直至最后一个H值，重复步骤E-b)，分别获取指定的H序列的第二、第三、……直至最后一个H值下的(M-T)_H测试数据文件；

这样便得到了按指定磁场强度H序列单独存储、且H值间隔均匀的(M-H)_T数据文件，可以构成(H，M，T)三维数组；

F)依次逐行读取步骤E)得到的H值间隔均匀的(M-H)_T数据文件，求解每一个H值下的M对温度T的偏导数具体步骤为：

在外部存储器中建立一个用于存储的数据的文件夹；

逐行读取步骤E)得到的H值间隔均匀的(M-H)_T数据文件，除去第一行标题行外，将读取到的每一行的数据存储于一个临时数组中；

在外部存储器的用于存储的数据的文件夹中建立一个包含有该H值为文件名的数据文件，其扩展名可以任意设定，如.txt，.dat，.bin等等；在所建立的数据文件中，以追加的方式写入标题行的内容，标题行包括磁场强度H和偏导数两个物理量及分隔符；

从每个(M-H)_T数据文件的第三行起，将读取到的该行的M值减去上一行读取到的M值，记为ΔM；将读取到的该行的H值减去上一行读取到的H值，记为ΔT；利用下面的公式来计算每一个指定H值下的M对温度T的偏导数：

{(\frac{\partial M (H, T)}{\partial T})}_{H} = \frac{Δ M}{Δ T};

将T值及相应的偏导数以逐行追加的方式写入步骤F)建立的文件夹下的数据文件中；

重复上述过程，遍历其余所有的(M-H)_T数据文件，获得全部指定H序列下的数据文件；

G)指定T序列，T序列可以是步骤B)中所设定的T序列，或者是步骤B)中所设定的T序列的子集；设定积分公式中的积分下限和上限；依次逐行读取步骤F)中建立的数据文件，对步骤F)中得到的一系列数据曲线在指定条件下求积分，得到指定T序列下的一系列的ΔS-T数据，具体步骤如下：

在外部存储器中建立一个用于存储的ΔS-T数据文件的文件夹；

依此遍历并逐行读取步骤F)得到的数据文件，除去第一行标题行外，将读取到的每一行的数据存储于一个临时数组中；

在外部存储器的用于存储的ΔS-T数据的文件夹中建立一个包含有步骤G)所设定的该积分上限H值为文件名的数据文件，其扩展名可以任意设定，如.txt，.dat，.bin等等；在所建立的数据文件中，以追加的方式写入标题行的内容，标题行包括温度T和磁熵变两个物理量及分隔符，例如：“TemperatureΔS”；

逐行读取步骤F)得到的数据文件，从第二行开始，将读取到的每一行的数据存储于一个临时数组中，当读取到T值等于步骤G)中指定T序列中的T值时，以一临时变量来存储该行的值；终止继续读取该数据文件，继续读取下一个数据文件，重复本段所述的操作，将所有读取到的指定T值所对应的进行累加，直至遍历完全部的数据文件，得到其累加值

Σ {(\frac{\partial M (H, T)}{\partial T})}_{H};

将上述累加值乘以H步长，即其乘积即为指定T值、积分下限、上限条件下的磁熵变ΔS，其中ΔH为相邻两行的H值的差值，因为在步骤E)中已经指定H序列为等间隔的序列，因此ΔH是一个固定不变的值；

将所得到的温度T、磁熵变ΔS以及分隔符，以逐行追加的方式写入步骤G)所建立的文件下的ΔS-T数据文件中；

重复上述步骤G)所述过程，遍历其余所有指定T序列下的一系列的数据文件，获得全部指定条件下的一系列的ΔS-T数据；

保持积分下限不变，更改积分上限，重复G步骤)中上述步骤，获得新的指定磁场变化条件下的全部指定T序列下的一系列的ΔS-T数据；

H)关闭上述所有步骤中读、写过的数据文件。

在步骤E-b1)中，所采取的M取值的依据在于：由于在一条M-H曲线中通常都包含有成千上万个数据点，所指定的H值与最邻近的实际测量H值之间的差值与整条M-H曲线中的H变化范围相比是微不足道的，因此，即使取该行所在的M值作为指定H值下所对应的M值，所导致的误差也是很小的，通常，由此带来的误差小于0.07％。

当计算精度要求较高时，可采取步骤E-b2)中所述的差值的方法来获得更为准确的M值。

优选地，在步骤D)中，所读取到的M值除以将待测试的磁性材料的质量，其结果替代原实测的M值，写入步骤D)建立的数据文件中。

优选地，步骤C)中得到的数据中，过长位数的数据并不能带来计算结果准确度的提高，为了为了减小数据量和运算资源消耗，在步骤D中将M、H的数值统一截取为固定长度的数值，例如小数点后取6位、7位或8位。

本发明的有益效果：

通过采用本发明的方法，可以彻底解决现有技术在利用PPMS测量磁性材料过程中获取的M-T曲线、ΔS-T曲线难题，利用计算机来完成相关操作，方便、快捷。

附图说明：

图1：某种磁性材料的M-T曲线；

图2：某种磁性材料的ΔS-T曲线；

具体实施方式

为便于理解本发明，下面结合实例来具体介绍本发明的技术方案。

A)按照符合PPMS操作规范的流程，称量5g的某种磁性材料粉末，严格按照PPMS操作规范将待测试的磁性材料装入PPMS测量系统中；

以下所有操作步骤均利用计算机来完成：

B)在PPMS的控制计算机上设定测试的温度序列值：2K，4K，6K，8K，10K，……48K，50K，52K，54K，……96K，98K，100K，各个测试温度成等间距序列，间隔为2K；设定测试的磁场强度为70000Oe，磁场从0Oe增加到70000Oe，再逐渐减小到-70000Oe，然后回到70000Oe；在每一个设定的测试温度下，磁性材料均历经同样的所述磁场扫描；所述T指PPMS系统导出的数据文件中的物理量Temperature；

C)开始测试：系统温度在控温系统的作用下达到设定的第一温度值2K并稳定后，测量第一M-H曲线；然后，系统温度在控温系统的作用下达到设定的第二温度值4K并稳定后，测量第二M-H曲线；然后，系统温度在控温系统的作用下达到设定的第三温度值6K并稳定后，测量第三M-H曲线；……重复这一过程，逐次测量磁性材料的M-H曲线，直至完成100K温度下的M-H曲线；在每次测量过程中温度保持恒定不变；依次逐次测量磁性材料的M-H曲线，直至完成所设定的所有温度序列下的M-H曲线；所述H指PPMS系统导出的数据文件中的物理量MagneticField，所述M指PPMS系统导出的数据文件中的物理量Moment；

依照PPMS操作规范的流程，导出所测试材料的所有数据，除标题行外，每一行数据中由温度值、所施加的外磁场强度、所测磁性材料的磁矩、误差项构成；根据测量数据量的多少，导出的测试数据文件得到的数据文件有7个，分别保存为MCE1.dat，MCE2.dat，MCE3.dat，MCE4.dat，MCE5.dat，MCE6.dat，MCE7.dat；不同的温度、不同的H值下的测试数据合并存储在所述一个或多个数据文件中，即使导出的数据文件为多个，在每个单独的数据文件中，也包含有不同的温度、不同的H值下的测试数据，且这些H值的间隔(ΔH)是不均匀的；

例如：MCE4.dat的文件部分内容如下：

Temperature(K),MagneticField(Oe),Moment(emu),M.Std.Err.(emu)

32.00202755,-0.118,-0.00126688856058318,4.51899702302814E-6

32.0013313,7.0295,-0.00127597555348906,4.86920603808076E-6

32.00174335,565.58,0.0422035493506414,0.000402803915835685

32.0015945,1025.09,0.0729658843502942,0.000458270810619051

32.001915,1255.996,0.0928575299552316,0.000451647444187443

32.0020218,1451.7385,0.112409561976028,0.000444269277498134

……

31.99976255,61045.3245,4.15532584682656,0.000495170438107461

31.9995384,60795.461,4.14521563817148,0.000600448386830866

31.999567,60635.5585,4.13403371695014,0.000654675443479942

31.9996605,60448.9255,4.12348915203937,0.000595515738766516

31.9995708,60272.258,4.11842844701878,0.000676258986960682

31.9993477,60103.1605,4.10819758191176,0.000710948878591925

……

34.0042763,7076.996,0.527432757413177,0.000347193892695242

34.0044746,7269.762,0.542291854429929,0.000342030768021547

34.0042934,7361.835,0.557095205406481,0.000343129104902078

34.0041885,7575.8715,0.564461725900955,0.000331726538447699

34.00409885,7648.111,0.579335477010551,0.000345789072590493

34.0038986,7906.3135,0.594013092437098,0.000347164011560761

34.0036354,8121.3095,0.608767264647458,0.000339320358035022

……

36.0040131,24387.041,1.71385973148928,0.000361950939258065

36.0034981,24582.1975,1.72698549024281,0.000356809845964897

36.003273,24802.459,1.74071538042355,0.000372159199191297

36.0034294,24997.042,1.75377960350089,0.0003968064715188

36.0033073,25199.7705,1.76690363483609,0.000361884393402058

36.0032425,25385.827,1.7807557086519,0.000414432388799402

36.003231,25607.715,1.79401333551337,0.000398482166958347

36.0033607,25800.865,1.80709933818233,0.000397470439631694

36.0035381,26074.1085,1.82711095863772,0.000376102430731429

……

38.00268745,36812.916,2.42817117789424,0.000443950779835792

38.0027237,37032.891,2.43384123902899,0.000430941893335765

38.0025177,37227.664,2.44604049048833,0.00043964961855296

38.0028267,37413.3395,2.45793496722591,0.000523089852077352

38.00296595,37488.2615,2.47011950372595,0.000507111775134114

38.00333215,37628.906,2.47642474163063,0.00043463581691392

38.00377655,37756.8065,2.48184218842265,0.000441226548974584

38.00416185,37977.6405,2.49401960102861,0.000448238522768833

……

40.0026207,69337.508,4.04897113030786,0.000630537000910383

40.0025177,69412.051,4.05820400138428,0.000568074112521435

40.0022697,69635.313,4.06209889195388,0.000543928762905382

40.00183485,69789.746,4.07179846876265,0.00071221520587325

40.0017433,69901.8945,4.07808757914927,0.000505949420696191

40.0022182,69999.922,4.08278275651664,0.000600963729794102

很明显地，单个MCE4.dat文件中合并存储了32K、34K、36K、38K、40K等五个温度下测试得到的M-H数据，并且这些H值的间隔(ΔH)是不均匀的；

D)在硬盘中建立文件夹C:\Data\M-H\，分别逐行读取所导出的测试数据MCE1.dat，MCE2.dat，MCE3.dat，MCE4.dat，MCE5.dat，MCE6.dat，MCE7.dat，利用计算机按照测试温度进行数据分割，将在每个测试温度(2K，4K，6K，……，98K，100K)下测量得到的(M-H)_T数据分别存储在各自独立的(M-H)_T文件中，这些分割得到的(M-H)_T文件(分别命名为2K.dat，4K.dat，6K.dat，……50K.dat，52K.dat，……98K.dat，100K.dat)存储在硬盘建立文件夹C:\Data\M-H\文件夹中；在每个独立的(M-H)_T文件(2K.dat，4K.dat，6K.dat，……，98K.dat，100K.dat)中，只存储在测量得到的从正向最大H值扫描到负向最大H值的过程的数据，即从70000Oe到-70000Oe的测试数据；具体步骤为：

D-a)从硬盘中逐行读取所导出的测试数据MCE1.dat到计算机内存(RAM)中；当读取第一行(即标题行)时，暂存第一行(即标题行)的文本内容，内容为“Temperature(K),MagneticField(Oe),Moment(emu)”；

D-b)当逐行读取测试数据文件的第二行至最后一行时，以PPMS数据导出时所采用的数据分隔符为不同物理量的分隔判别标志，以步骤C)中所设定的分隔符来区分开不同物理量(Temperature、MagneticField、Moment等)，将每一行的原始数据存储于一个临时数组中；

MCE1.dat的部分内容如下：

Temperature(K),MagneticField(Oe),Moment(emu),M.Std.Err.(emu)——第1行

2.00786255,-0.062,-0.0018420076453778,4.68147230612165E-6——第2行

2.0080564,54.536,-0.00160844591554442,0.000247550499783577——第3行

2.0114515,455.301,0.0390593777601995,0.000711834956395969——第4行

2.0107373,525.147,0.054575639945636,0.000707314880304041

2.01013325,643.76,0.0699908666254729,0.000700355705591364

2.0099548,762.4715,0.0852845765830705,0.000694930050706897

2.0094968,891.817,0.100431437809785,0.000691586807257692

2.0088098,1035.1515,0.115400953689561,0.00068140220226271

2.00849115,1162.775,0.130234228395757,0.000673690201895222

2.00874615,1263.1835,0.144902554721851,0.000664242728502225

……

1.99570285,69982.0585,6.47701266154109,0.000910002725693407

1.9964498,70000.418,6.47737046047771,0.00100083032970079——第741行

1.997007,69974.453,6.47764781770063,0.000882885617488435——第742行

1.99819445,69895.484,6.47590669194502,0.00112676914363876——第743行

……

1.9975571,-70014.594,-6.47646090206977,0.0010909447795878

1.9979685,-70045.5355,-6.47722052937364,0.00110168935393343

1.9983413,-70056.3205,-6.47672913871526,0.00117970710047581——第2121行

1.99813915,-70000.57,-6.47673318183019,0.00101478431642557——第2122行

1.9984763,-69949.7855,-6.47683668464876,0.00104135098506153

……

4.0043831,-0.335,0.00401174413532849,4.68827127114532E-5——第3474行

4.003768,65.656,0.00376830558528498,0.000854753231569799——第3475行

4.0111647,446.867,0.0442976647184307,0.000574437113725152——第3476行

4.0035796,920.939,0.108339487673647,0.000711857730897711

4.00412775,993.0845,0.123853303763718,0.000705167335263082

4.00524405,1110.071,0.139273263828744,0.000705910455217505

4.00522045,1210.5775,0.154591409474921,0.000703423510718862

4.004051,1325.4535,0.16976336326451,0.000689500239741493

4.0041368,1454.0315,0.184904992382575,0.000693296538700884

4.00412225,1577.531,0.199936441375344,0.000688937041088548

4.003361,1700.839,0.214808299881983,0.000681304686754516

4.00367215,1823.7655,0.229638867059882,0.000671341384018587

4.00299145,2260.854,0.273674767571859,0.000663858971864902

……

10.0039368,-0.193,-0.00126743023287644,4.67445883867269E-6——第13906行

10.002696,64.1845,-0.00128526790687863,4.56030727999837E-6——第13907行

10.0027289,92.4845,0.00617136531883082,0.00121020852470981——第13908行

10.0064507,631.5895,0.0848930971419456,0.00088703653307265

10.00640635,693.3865,0.104340095631526,0.000884500313919431

10.00629905,799.9285,0.123644325266753,0.000877777808707441

10.0061498,928.0285,0.142794671249241,0.000874610713522662

10.00592045,1056.32,0.161861059302663,0.000871550237515841

10.0056553,1183.941,0.180684210640686,0.000858980501435077

……

则，读取到第2行时，实测温度值2.00786255K，与设定的温度序列(2K，4K，6K，……，98K，100K)中最为接近的是2K，因此温度值记为2K，在硬盘上文件夹C:\Data\M-H\文件夹中建立一个以2K.dat为文件名的数据文件；在所建立的数据文件2K.dat中，以追加的方式写入步骤a)中暂存的第一行即标题行的内容“Temperature(K),MagneticField(Oe),Moment(emu)”(不储存误差项)；在一临时数组中记录H、M的实测值分别为-0.062Oe和-0.0018420076453778emu；分隔符为逗号“,”；

读取到第3行时，实测温度值2.0080564K，与设定的温度序列(2K，4K，6K，……，98K，100K)中最为接近的还是2K，即设定的温度值没有改变，在一临时数组中记录H、M的实测值分别为54.536Oe和-0.00160844591554442emu；将本次读取到的H值与上一次读取到的H值进行比较，本次读取到的H值的绝对值大于上一次读取到的H值的绝对值，不作任何操作，继续读取下一行；

读取到第4行时，实测温度值2.0114515K，与设定的温度序列(2K，4K，6K，……，98K，100K)中最为接近的还是2K，即设定的温度值没有改变，在一临时数组中记录H、M的实测值分别为455.301Oe和0.0390593777601995emu；将本次读取到的H值与上一次读取到的H值进行比较，本次读取到的H值的绝对值大于上一次读取到的H值的绝对值，不作任何操作，继续读取下一行；

……

读取到第741行时，实测温度值1.9964498K，与设定的温度序列(2K，4K，6K，……，98K，100K)中最为接近的还是2K，即设定的温度值没有改变，在一临时数组中记录H、M的实测值分别为70000.418Oe和6.47737046047771emu；将本次读取到的H值与上一次读取到的H值进行比较，本次读取到的H值的绝对值大于上一次读取到的H值的绝对值，不作任何操作，继续读取下一行；

读取到第742行时，实测温度值1.997007K，与设定的温度序列(2K，4K，6K，……，98K，100K)中最为接近的还是2K，即设定的温度值没有改变，在一临时数组中记录H、M的实测值分别为69974.453Oe和6.47764781770063emu；将本次读取到的H值与上一次读取到的H值进行比较，本次读取到的H值的绝对值小于上一次读取到的H值的绝对值，即意味着上一行读取到的H值为2K温度下测量M-H曲线中的正向最大值，在C:\Data\M-H\2K.dat文件中以追加的方式写入上一行读取到的H、M数值70000.418Oe和6.47737046047771emu；H、M数值之间的分隔符为逗号“,”；

读取到第743行时，实测温度值1.99819445K，与设定的温度序列(2K，4K，6K，……，98K，100K)中最为接近的还是2K，即设定的温度值没有改变，在一临时数组中记录H、M的实测值分别为69895.484Oe和6.47590669194502emu；将本次读取到的H值与上一次读取到的H值进行比较，本次读取到的H值的绝对值小于上一次读取到的H值的绝对值，在C:\Data\M-H\2K.dat文件中以追加的方式写入上一行读取到的H、M数值69974.453Oe和6.47764781770063emu；H、M数值之间的分隔符为逗号“,”；

……

读取到第2121行时，实测温度值1.9983413K，与设定的温度序列(2K，4K，6K，……，98K，100K)中最为接近的还是2K，即设定的温度值没有改变，在一临时数组中记录H、M的实测值分别为-70056.3205Oe和-6.47672913871526emu；将本次读取到的H值与上一次读取到的H值进行比较，本次读取到的H值的绝对值大于上一次读取到的H值的绝对值，即意味着上一行读取到的H值为2K温度下测量M-H曲线中的负向最大值，在C:\Data\M-H\2K.dat文件中以追加的方式写入上一行读取到的H、M数值-70045.5355Oe和-6.47722052937364emu；H、M数值之间的分隔符为逗号“,”；

读取到第2122行时，实测温度值1.99813915K，与设定的温度序列(2K，4K，6K，……，98K，100K)中最为接近的还是2K，即设定的温度值没有改变，在一临时数组中记录H、M的实测值分别为-70000.57Oe和-6.47673318183019emu；将本次读取到的H值与上一次读取到的H值进行比较，本次读取到的H值的绝对值小于上一次读取到的H值的绝对值，在C:\Data\M-H\2K.dat文件中以追加的方式写入上一行读取到的H、M数值-70056.3205Oe和-6.47672913871526emu；H、M数值之间的分隔符为逗号“,”；

读取到第2123行时，实测温度值1.9984763K，与设定的温度序列(2K，4K，6K，……，98K，100K)中最为接近的还是2K，即设定的温度值没有改变，不作任何操作；

……

读取到第3474行时，实测温度值4.0043831K，与设定的温度序列(2K，4K，6K，……，98K，100K)中最为接近的是4K，即设定的温度值发生了改变，因此温度值记为4K，在硬盘上文件夹C:\Data\M-H\文件夹中建立一个以4K.dat为文件名的数据文件；以追加的方式写入步骤a)中暂存的第一行即标题行的内容“Temperature(K),MagneticField(Oe),Moment(emu)”(不储存误差项)；

其后操作与上面的类似；

……

读取到第13906行时，实测温度值10.0039368K，与设定的温度序列(2K，4K，6K，……，98K，100K)中最为接近的是10K，即设定的温度值发生了改变，因此温度值记为10K，在硬盘上文件夹C:\Data\M-H\文件夹中建立一个以10K.dat为文件名的数据文件；以追加的方式写入步骤a)中暂存的第一行即标题行的内容“Temperature(K),MagneticField(Oe),Moment(emu)”(不储存误差项)；

其后操作与上面的类似；

……

直至MCE1.dat的最后一行结束；

D-c)依次类推，重复执行步骤E-b)，直至遍历完其余的导出的测量数据文件MCE2.dat，MCE3.dat，MCE4.dat，MCE5.dat，MCE6.dat，MCE7.dat；

D-d)关闭上述过程中所有打开的数据文件MCE1.dat，MCE2.dat，MCE3.dat，MCE4.dat，MCE5.dat，MCE6.dat，MCE7.dat，以及2K.dat，4K.dat，6K.dat，……50K.dat，52K.dat，……98K.dat，100K.dat；

E)获得指定等间距H序列(70000Oe，69800Oe，69600Oe，69400Oe，……0Oe，……-69600Oe，-69800Oe，-70000Oe，即H间隔为200Oe)下的一系列的(M-T)_H数据文件，构成了(H，M，T)三维数组：根据步骤E)分割得的代表不同温度下测量得到的(M-H)_T数据文件2K.dat，4K.dat，6K.dat，……50K.dat，52K.dat，……98K.dat，100K.dat，利用计算机按照指定的H变化步长(即对指定H序列)进行数据采样，获得指定H序列下的一系列的(M-T)_H数据文件；步骤E)具体步骤如下：

E-a)指定一系列的H序列值：70000Oe，69800Oe，69600Oe，69400Oe，……0Oe，……-69600Oe，-69800Oe，-70000Oe，所述H序列值以恒定步长ΔH＝200Oe递增；以包含指定H序列的H值为文件名，在硬盘中建立文件夹C:\Data\M-T\，在该文件夹中建立一系列代表各自H值下测量得到的独立的数据文件：70000Oe.dat，69800Oe.dat，69600Oe.dat，69400Oe.dat，……0Oe.dat，……-69600Oe.dat，-69800Oe.dat，-70000Oe.dat；

在70000Oe.dat，69800Oe.dat，69600Oe.dat，69400Oe.dat，……0Oe.dat，……-69600Oe.dat，-69800Oe.dat，-70000Oe.dat数据文件中，以追加的方式写入标题行，标题行包含M、T两个物理量，如“Temperature(K),Moment(emu)”；

E-b)获取指定H序列中第一个指定H值70000Oe下的(M-T)_H测试数据序列，并将所述(M-T)_H测试数据序列以追加的方式写入在步骤F-a)中所建立的对应H值的数据文件中70000Oe.dat：

依照温度从低到高的次序，遍历读取步骤D)中存储在硬盘文件夹C:\Data\M-H\中、代表在每个测试温度下测量得到的、各自独立的(M-H)_T数据文件：2K.dat，4K.dat，6K.dat，……50K.dat，52K.dat，……98K.dat，100K.dat；每个文件的读取方式均为逐行读取；

例如：2K.dat的部分数据如下：

MagneticField(Oe),Moment(emu)——第1行

70000.418,6.47737046047771——第2行

69974.453,6.47764781770063——第3行

69895.484,6.47590669194502——第4行

69792.5035,6.4736558963573——第5行

69694.109,6.47083457293903——第6行

69577.418,6.46978183022505——第7行

69474.805,6.46605271136606——第8行

69365.7735,6.4648762127737——第9行

69258.1835,6.4626036312054——第10行

69151.3595,6.4609293819902——第11行

69055.266,6.45822662541325——第12行

68960.801,6.45632540543144——第13行

68854.2615,6.4553284188999——第14行

68750.8825,6.4529345096850——第15行

……

4K.dat的部分数据如下：

MagneticField(Oe),Moment(emu)——第1行

70000.6285,6.45364032917115——第2行

69947.031,6.451612849781050——第3行

69891.5625,6.45188523236246——第4行

69789.1485,6.44881227197867——第5行

69680.2185,6.44725088629155——第6行

69584.316,6.44350190352285,——第7行

69481.5115,6.44309404080617——第8行

69342.7815,6.44011894070729——第9行

69254.352,6.438656685704360——第10行

69176.844,6.435399066070250——第11行

69052.965,6.433575860754200——第12行

68948.1525,6.43176163955381——第13行

68862.5,6.43122044221395000——第14行

……

当逐行读取2K.dat并记录每一行中的H值，将所读取到的H值与指定H值70000Oe进行比较，得出2K.dat中第2行即数据记录的第1行中实际测量得到的H值70000.418Oe与指定值70000Oe最接近，因此以实际测量得到的H值70000.418Oe所对应的M值6.47737046047771emu作为指定H值70000Oe下所对应的M值，并将对应的温度值2K、M值6.47737046047771emu以追加的方式写入数据文件70000Oe.dat中；

继续逐行读取4K.dat并记录每一行中的H值，将所读取到的H值与指定H值70000Oe进行比较，得出4K.dat中第2行即数据记录的第1行中实际测量得到的H值70000.6285Oe与指定值70000Oe最接近，因此以实际测量得到的H值70000.6285Oe所对应的M值6.45364032917115emu作为指定H值70000Oe下所对应的M值，并将对应的温度值4K、M值6.45364032917115emu以追加的方式写入数据文件70000Oe.dat中；

依照上述操作方式，遍历所有其余的数据文件6K.dat，……50K.dat，52K.dat，……98K.dat，100K.dat，便得到70000Oe下的(M-T)_H数据文件70000Oe.dat；

E-c)获取指定H序列中第二个指定H值69800Oe下的(M-T)_H测试数据序列，并将所述(M-T)_H测试数据序列以逐行追加的方式写入在步骤E-a)中所建立的对应H值的数据文件中69800Oe.dat：

当逐行读取2K.dat并记录每一行中的H值，将所读取到的H值与指定H值69800Oe进行比较，得出2K.dat中第5行即数据记录的第4行中实际测量得到的H值69792.5035Oe与指定值69800Oe最接近，因此以实际测量得到的H值69792.5035Oe所对应的M值6.4736558963573emu作为指定H值69800Oe下所对应的M值，并将对应的温度值2K、M值6.4736558963573emu以追加的方式写入数据文件69800Oe.dat中；

继续逐行读取4K.dat并记录每一行中的H值，将所读取到的H值与指定H值69800Oe进行比较，得出4K.dat中第5行即数据记录的第4行中实际测量得到的H值69789.1485Oe与指定值69800Oe最接近，因此以实际测量得到的H值69789.1485Oe所对应的M值6.44881227197867emu作为指定H值70000Oe下所对应的M值，并将对应的温度值4K、M值6.44881227197867emu以追加的方式写入数据文件69800Oe.dat中；

依照上述操作方式，遍历所有其余的数据文件6K.dat，……50K.dat，52K.dat，……98K.dat，100K.dat，便得到70000Oe下的(M-T)_H数据文件69800Oe.dat；

依照同样的方法，获取指定H序列中69600Oe、69400Oe、……0Oe、……、-69600Oe、-69800Oe、-70000Oe下的(M-T)_H测试数据序列69600Oe.dat，69400Oe.dat，……0Oe.dat，……-69600Oe.dat，-69800Oe.dat，-70000Oe.dat；

关闭步骤E-b)中所有读、写过的文件，获得指定H序列中第一个指定H值下的(M-T)_H测试数据文件；

这样便得到了获得指定等间距H序列下的一系列的(M-T)_H数据文件，构成了(H，M，T)三维数组；

F)依次逐行读取步骤E)得到的(M-H)_T数据文件，求解每一个H值下的M对温度T的偏导数具体步骤为：

在硬盘上建立一个用于存储的数据的文件夹C:\Data\dM-dT\；

在外部存储器的用于存储的数据的文件夹中建立一个包含有该H值为文件名的数据文件200Oe.txt等等；在所建立的数据文件200Oe.txt中，以追加的方式写入标题行的内容，标题行为“Temperature(K)Slope(emu/(g*K))”；

{(\frac{\partial M (H, T)}{\partial T})}_{H} = \frac{Δ M}{Δ T};

例如，在步骤E)得到的200Oe下的M-T文件内容如下：

Temperature(K)Moment(emu)M.Std.Err.(emu)

25.05231864850.000912099112341484

44.92286511290.000848345823308769

64.70764960220.000951939432121166

84.43557178510.00104602122157555

104.15401647070.000756497790415489

123.90869667600.000759289863431564

143.65791397450.00065282260535558

163.43580357220.000712153496413619

183.23363921520.000641743252684032

203.05421860570.000711178731830154

222.88750547400.000633583014363018

242.73639440390.000612027386979506

262.60672229680.000582255503098482

282.48583058500.00052096015783399

302.37073827870.000531041995268534

322.27086963940.000527704246626589

342.17221313490.000468627209191499

362.08564902320.000424715399090915

382.00703910060.000444781136608427

401.93469572990.00036176508505969

421.86875007680.000374445379170067

441.80827090580.00032804006324117

461.74592611050.000339148252897741

481.69603536590.000326136824685095

501.64605452530.000354802492701038

551.53583709770.000322637398867676

601.44111776090.000280788103370659

651.36024897930.000271314849759564

701.28718586370.000256295679790413

751.22053431490.000230965052151243

801.16399353990.000242077424031515

851.10902789770.00023282699142299

901.05917765310.000240851206211287

951.01551616480.000161831296804004

1000.97332717680.000194509617496906

第3行的M值减去第2行读取到的M值，ΔM＝-0.129453536emu；

第3行的T值减去第2行读取到的T值，ΔT＝2K；

则

{(\frac{\partial M (H, T)}{\partial T})}_{H} = \frac{Δ M}{Δ T} = - 0.064726768 e m u / K

重复上述过程，遍历其余所有的(M-H)T数据文件，获得全部指定H序列下的数据文件：0Oe.txt、200Oe.txt、400Oe.txt、600Oe.txt、……、69400Oe.txt、69600Oe.txt、69800Oe.txt、70000Oe.txt；

举例而言，得到的36000Oe下的数据结果如下：

Temperature(K)PartialDerivative(emu/(g*K))

4-0.0655024472

6-0.1003210984

8-0.1221701159

10-0.1236801023

12-0.1196835623

14-0.1221636892

16-0.1158554408

18-0.1031784337

20-0.0940045455

22-0.0852988782

24-0.0791938393

26-0.0729467498

28-0.0657919604

30-0.0608403727

32-0.0555074351

34-0.0511001996

36-0.0477874164

38-0.0443303011

40-0.0407991128

42-0.0380759010

44-0.0336589478

46-0.0308700606

48-0.0293266307

50-0.0260036567

55-0.0240959068

60-0.0212309865

65-0.0182615201

70-0.0165310821

75-0.0151322573

80-0.0133145686

85-0.0120863879

90-0.0110094887

95-0.0098515785

特别指出的是，这里所获得的各个数据文件的名称与步骤E)所获得的各个数据文件的名称分别相同，但并不会造成混淆，因为这两个步骤中产生的数据文件存放于各自的文件夹中；

G)指定T序列：4K、6K、8K、10K、……、46K、48K、50K、……94K、96K、98K、100K，设定积分公式中的积分下限和上限分别为0Oe、600Oe；依次逐行读取步骤F)中建立的数据文件，对步骤F)中得到的一系列数据曲线在指定条件下求积分，得到指定T序列下的一系列的ΔS-T数据，具体步骤如下：

在外部存储器中建立一个用于存储的ΔS-T数据文件的文件夹C:\Data\ΔS-T\；

依此遍历并逐行读取步骤F)得到的数据文件:0Oe.txt、200Oe.txt、400Oe.txt、600Oe.txt、……、69400Oe.txt、69600Oe.txt、69800Oe.txt、70000Oe.txt，除去第一行标题行外，将读取到的每一行的数据存储于一个临时数组中；

在外部存储器的用于存储的ΔS-T数据的文件夹中建立一个包含有步骤G)所设定的该积分上限H值为文件名的数据文件：0Oe.txt、200Oe.txt、400Oe.txt、600Oe.txt、……、69400Oe.txt、69600Oe.txt、69800Oe.txt、70000Oe.txt；在所建立的数据文件中，以追加的方式写入标题行的内容，标题行包括温度T和磁熵变两个物理量及分隔符，例如：“TemperatureΔS”；

Σ {(\frac{\partial M (H, T)}{\partial T})}_{H};

例如：计算10K温度下，该材料在0Oe至600Oe之间的磁熵变：

步骤F)得到的0Oe.txt、200Oe.txt、400Oe.txt、600Oe.txt、800Oe部分内容如下：

0Oe.txt部分内容：

Temperature(K)Slope(emu/(g*K))syd4182163.com200syd4182163.com

4-0.0040352674——第2行

6-0.0020341544——第3行

8-0.0003398646——第4行

10-0.0008626960——第5行

120.0001786028

……

200Oe.txt部分内容：

Temperature(K)Slope(emu/(g*K))syd4182163.com200syd4182163.com

40.0030047162——第2行

60.0023498100——第3行

80.0004700273——第4行

10-0.007406361——第5行

12-0.0080658803

……

400Oe.txt部分内容：

Temperature(K)Slope(emu/(g*K))syd4182163.com200syd4182163.com

4-0.0015049657——第2行

60.0008345383——第3行

8-0.0009918612——第4行

10-0.0017225218——第5行

12-0.0010413388

……

600Oe.txt部分内容：

Temperature(K)Slope(emu/(g*K))syd4182163.com200syd4182163.com

4-0.0011649709——第2行

60.0019778708——第3行

8-0.0006583236——第4行

10-0.0008542947——第5行

12-0.0013106532

……

800Oe.txt部分内容：

Temperature(K)Slope(emu/(g*K))syd4182163.com200syd4182163.com

40.0001960065——第2行

60.0043273019——第3行

8-0.0009248119——第4行

10-0.0020094518——第5行

12-0.0020280588

……

逐行读取上述数据文件0Oe.txt、200Oe.txt、400Oe.txt、600Oe.txt，由于数据文件600Oe.txt所代表的测试H值已经达到所指定的积分上限，因此800Oe.txt、1000Oe.txt、1200Oe.txt、……、69800Oe.txt、70000Oe.txt等数据文件则不再读取；对应于10K温度的数据均位于各自数据文件的第5行，将这些文件的第5行中记录的偏导数的值进行累加，得而H的步长ΔH在步骤E)中已经设定为200Oe，因此，在10K温度下，该材料在0Oe至600Oe之间的磁熵变为：

Δ S = Σ {(\frac{\partial M (H, T)}{\partial T})}_{H} \times Δ H = - 2.1691747 (e m u * O e / (g * K)

重复上述步骤G)所述过程，遍历其余所有指定T序列(6K、8K、10K、12K、……、96K、98K、100K)下的一系列的数据文件，获得全部指定条件下的一系列的ΔS-T数据；

保持积分下限不变，更改积分上限为12000Oe、16000Oe、20000Oe、24000Oe、28000Oe、35000Oe、45000Oe、55000Oe、65000Oe、70000Oe，重复G步骤)中上述步骤，获得新的指定磁场变化条件下的全部指定T序列下的一系列的ΔS-T数据：ΔS-T12000Oe.txt、ΔS-T16000Oe.txt、ΔS-T20000Oe.txt、ΔS-T24000Oe.txt、ΔS-T28000Oe.txt、ΔS-T35000Oe.txt、ΔS-T45000Oe.txt、ΔS-T55000Oe.txt、ΔS-T65000Oe.txt、ΔS-T70000Oe.txt；

举例而言，ΔS-T55000Oe.txt的内容如下：

Temperature(K)ΔS(J/kgK)

4-0.2467038758

6-0.3696258054

8-0.4703001883

10-0.4802541481

12-0.4824572476

14-0.4888180513

16-0.4544233692

18-0.4184906003

20-0.3841392399

22-0.3540794608

24-0.3222721539

26-0.2936000779

28-0.2698646268

30-0.2503259984

32-0.2309927611

34-0.2140408844

36-0.1997334306

38-0.1833204023

40-0.1670496952

42-0.1545068372

44-0.1416488434

46-0.1309885549

48-0.1239820804

50-0.1098914411

55-0.1000167777

60-0.0879078722

65-0.0773784541

70-0.0692550325

75-0.0619101999

80-0.0560074570

85-0.0510289580

90-0.0462624511

95-0.0426979831

可以将上述数据导入Origin等数据分析软件中进行处理、分析；

H)关闭上述所有步骤中读、写过的数据文件。

本发明的处理方法，在IntelCore2DuoP8700、2GBDDR21333、7200转500GB硬盘+Windows7professional软硬件平台上处理上述实施例的数据时，以硬盘暂存中间数据文件的技术方案，整个过程耗时2min37s。由此可见，本发明的技术方案，不仅解决了现有技术中PPMS测量得到的数据难以评估磁性材料的磁熵变—温度特性的难题，且其效率大大提升，极大地降低了科研人员的人力成本。

由于本发明的方法在步骤D)-步骤G)中需要多次地对硬盘读写零碎的数据，每次读写的数据只有一行，仅约数十字节，而数据运算则极为简单，因此整个过程中的耗时基本取决于硬盘读写操作。相对于现有技术中将全部数据导入计算机内存、现有技术中以多维数组进行处理的方法，本发明的方法耗时大于现有技术。

然而，本发明的技术方案对计算机内存的需求则远远小于现有技术，经实际测量，在数据处理的过程中，数据+代码运行所占据的全部内存不到4MB，。而现有技术中若达到同样的计算精度，则五维数组所需要的内存(RAM)空间大约在10¹⁰-10¹⁷MB之间，即使将测量过程中各个物理量的数量均降低1～2个数量级，五维数组所需要的内存(RAM)空间也仍然在10⁴-10¹¹MB之间，且结果精确度也会大大降低，将导致结果严重不可信。

因此即便本发明耗时较长，也仅仅在1-3min而已，这一缺点在科研过程中是微不足道的。

Claims

1.一种专用于PPMS的批量获取磁性材料的M-T曲线、ΔS-T曲线的方法，其特征在于：

该方法仅用于以美国QuantumDesign公司的PPMS(PhysicalPropertyMeasurementSystem，综合物性测量系统)进行磁性材料的磁学性质测量；且，

包括如下步骤：

A)按照符合PPMS操作规范的流程，称量将待测试的磁性材料的质量，记录下该质量值，将待测试的磁性材料装入PPMS测量系统中；

以下所有操作步骤均通过计算机来完成：

B)在PPMS的控制计算机上设定测试的温度序列值，各个测试温度成等间距序列(间隔为ΔT)，ΔT可以取1K、2K、3K、4K、5K、……、10K；设定测试的磁场强度，这里所述的磁场强度是指磁滞回线中的磁场强度的绝对值的最大值，磁场从0增加到正向最大值，再逐渐变化到负向最大值，然后再逐渐回到正向最大值；在每一个设定的测试温度下，磁性材料均历经同样的所述磁场扫描；

C)在给定温度下给磁性材料施加一随时间变化的外磁场H，测量在所述温度下的第一条M-H曲线(也可以称为M-H磁滞回线、磁化曲线)；上述测试完成后，使温度增加或减少一固定值，即温度等步长地增大或减小，重复进行上述测量步骤；在每次测量过程中温度保持恒定不变；依次逐次测量磁性材料的M-H曲线，直至完成所设定的所有温度序列下的M-H曲线；

D-a)从外部存储器中逐行读取所导出的测试数据到计算机内存(RAM)或CPU的高速缓存(Cache)中；当步骤C)中得到多个数据文件时，对多个数据文件轮次进行逐行读取；当读取第一行(即标题行)时，暂存第一行(即标题行)的文本内容；

D-b)当逐行读取测试数据文件的第二行至最后一行时，以PPMS数据导出时所采用的数据分隔符为不同物理量的分隔判别标志，以步骤C)中所设定的分隔符来区分开不同物理量，将每一行的原始数据存储于一个临时数组中；将所读取到的温度实测值取为测试前设定的温度序列中最接近的设定温度值，并在外部存储器中建立一个包含有该设定温度值为文件名的数据文件；在所建立的数据文件中，以追加的方式写入步骤D-a)中暂存的第一行即标题行的内容；

在上述步骤D-b)中，还包括：

D-d)关闭上述过程中所有打开的数据文件；

获得指定等间距H序列下的一系列的(M-T)_H数据文件，构成了(H，M，T)三维数组：根据步骤E)分割得的代表不同温度下测量得到的(M-H)_T数据文件，利用计算机按照指定的H变化步长ΔH(即对指定H序列)进行数据采样，获得指定H序列下的一系列的(M-T)_H数据文件；

F)依次逐行读取步骤E)得到的H值间隔均匀的(M-H)_T数据文件，求解每一个H值下的M对温度T的偏导数

G)指定T序列，T序列可以是步骤B)中所设定的T序列，或者是步骤B)中所设定的T序列的子集；设定积分公式dH中的积分下限和上限；依次逐行读取步骤F)中建立的数据文件，对步骤F)中得到的一系列数据曲线在指定条件下求积分，得到指定T序列下的一系列的ΔS-T数据；

H)关闭上述所有步骤中读、写过的数据文件；

所述H指PPMS系统导出的数据文件中的物理量MagneticField，即外加磁场强度；所述M指PPMS系统导出的数据文件中的物理量Moment，即磁性材料的磁矩；所述T指PPMS系统导出的数据文件中的物理量Temperature，即测试温度；所述S指磁性材料的磁熵；

2.一种如权利要求1所述的批量获取磁性材料的M-T曲线、ΔS-T曲线的方法，其特征在于：

步骤E)具体步骤如下：

在外部存储器所建立的数据文件中，以追加的方式写入标题行，标题行包含M、T两个物理量；

E-b)获取指定H序列中第一个指定H值下的(M-T)_H测试数据序列，并将所述(M-T)_H测试数据序列以追加的方式写入在步骤E-a)中所建立的对应H值的数据文件中，该操作具体步骤如下：

E-c)更改指定H值为步骤E-a)中指定的H序列的第二、第三、……直至最后一个H值，重复步骤E-b)，分别获取指定的H序列的第二、第三、……直至最后一个H值下的(M-T)_H测试数据文件，得到按指定磁场强度H序列单独存储、且H值间隔均匀的(M-H)_T数据文件。

3.一种如权利要求1所述的批量获取磁性材料的M-T曲线、ΔS-T曲线的方法，其特征在于：

步骤G具体步骤如下：

在外部存储器的用于存储的ΔS-T数据的文件夹中建立一个包含有步骤G)所设定的该积分上限H值为文件名的数据文件；在所建立的数据文件中，以追加的方式写入标题行的内容，标题行包括温度T和磁熵变两个物理量及分隔符；

将上述累加值乘以H步长，即其乘积即为指定T值、积分下限、上限条件下的磁熵变ΔS，其中ΔH是一个固定不变的值；

将所得到的温度T、磁熵变ΔS以及分隔符，以逐行追加的方式写入步骤G所建立的文件下的ΔS-T数据文件中；

重复上述步骤G所述过程，遍历其余所有指定T序列下的一系列的数据文件，获得全部指定条件下的一系列的ΔS-T数据；

4.一种如权利要求1或2所述的测量磁性材料的磁熵变-温度特性的方法，其特征在于：

步骤F具体步骤为：

在外部存储器中建立一个用于存储的数据的文件夹；

逐行读取步骤E得到的H值间隔均匀的(M-H)_T数据文件，除去第一行标题行外，将读取到的每一行的数据存储于一个临时数组中；

在外部存储器的用于存储的数据的文件夹中建立一个包含有该H值为文件名的数据文件；在所建立的数据文件中，以追加的方式写入标题行的内容，标题行包括磁场强度H和偏导数两个物理量及分隔符；

{(\frac{\partial M (H, T)}{\partial T})}_{H} = \frac{Δ M}{Δ T};

重复上述过程，遍历其余所有的(M-H)_T数据文件，获得全部指定H序列下的数据文件。

5.一种如权利要求1或2所述的测量磁性材料的磁熵变-温度特性的方法，其特征在于：

在步骤D中将M、H的数值统一截取为固定的5位或6位或7位或8位长度的数值。