Debye-Einstein-模型拟合比热容Python脚本

固体比热模型中的德拜模型和爱因斯坦模型是固体物理学中用于估算固体热容的两种重要原子振动模型。

爱因斯坦模型基于三种假设：1.晶格中的每一个原子都是三维量子谐振子；2.原子不互相作用；3.所有的原子都以相同的频率振动（与德拜模型不同）。

在高温下，爱因斯坦模型与实验结果一致，特别是与杜隆-珀替定律相符‌。

德拜模型将晶体中的原子振动视为连续弹性介质中传播的弹性波。固体的热容主要由低频的声学支声子贡献，存在截止频率，并未考虑光学支声子的贡献。在低温区与实验结果高度一致。

基本设置

import numpy as npimport osimport matplotlib.pyplot as pltfrom scipy.optimize import curve_fitimport scipy.integrate as integrate
from scipy.integrate import quad
R = 8.3144  # unit: J/ (mol·K）N = 10   #number of atomsn = 0.5 #Debye/(Debye+Einstein)names = ["data.dat"]colors = 'rgbmpyckrgbmpyc'

数据读入

def readata(name):