matlab在热学中的应用
198 MATLAB 及其在大学物理课程中的应用4-5 固体的热力学性质本节利用 MATLAB来处理固体热容量的三种模型、顺磁性固体及负温度状态。4.5.1 固体热容量的三种模型热容量是热力学系统的一个重要响应函数。经典理论曾用能量均分定理讨论了晶体在高温情况下的热容量，成功地解释了杜隆珀替定律。但是，经典理论不能说明低温下热容量随温度的降低而减小，以及它是系统特征量这两个实验事实。1907年，爱因斯坦应用量子概念处理晶体振动，定性地说明了固体的热容量随温度降低而趋于零的规律。1917 年，德拜修改了爱因斯坦模型，导出了3T定律，使固体热容量理论在定量上与实验结果相符合。(1)固体热容量的经典模型杜隆-珀替定律按照经典理论，由N 个原子或离子组成的固体可视为3N 个相互独立的经典线性谐振子的集合。由能量均分定理，每个线性简谐振子的能量为kT，固体的内能为U3NkT，热容量为3VCNk此即杜隆 - 珀替定律。 题目 ex4511 应用玻尔兹曼统计求经典固体的定容热容量。 解题分析经典固体可视为3N 个相互独立的经典线性谐振子的集合,每个经典线性谐振子的能量为222212rprem wm其中 ,212rpm是两原子相对运动的动能，1212m mmmm为约化质量 , r 是两原子间的距离， 为振动的圆频率。振动配分函数为第 4 章 MATLAB在热物理学中的应用199 22222 11ddrprv rZeprhbm wm- 蝌222 221ddrprreperhbbmw m-骣骣桫桫蝌求出配分函数后，再利用热力学公式13lnUNZ，V VUCT可求得经典固体的热容量。 程序 ex4511 syms V h beta N k T mu omiga r p; dbeta/2/mu;ebeta*mu*omiga2/2; zp2/d1/2*intexp-p2,0,inf; zr2/e1/2*intexp-r2,0,inf; Zv zp*zr/h, 振动配分函数Uv-3*N*difflogZv,beta; beta1/k/T; Uv1simplifyUv, 内能CvdiffUv1,T 热容量运行结果为Zv 2/beta/mu1/2*pi/beta*mu*omiga21/2/h Uv1 3*N*k*T Cv 3*N*k 实验表明， 杜隆 - 玻替定律在固体的温度较高时与测量结果符合，但在常温和低温下与实验结果严重不符。事实上，固体热容量是与温度和固体特性有关的量，并非该定律所描述的那样是与二者无关的常量。杜隆 - 玻替定律与实验事实偏离是对经典热力学理论的严重挑战。(2)爱因斯坦模型爱因斯坦将量子观点应用于固体热容量的研究，把固体看作由3N 个频率相同的，近独立的量子线性谐振子所组成的系统，应用玻尔兹曼统计得到了固体的内能和热容量表达式，这是继普朗克辐射理论之后，利用量子理论处理问题的第二个成功范例。200 MATLAB 及其在大学物理课程中的应用 题目 ex4512 应用玻尔兹曼统计求爱因斯坦固体的内能和定容热容量。 解题分析量子线性谐振子的能量为12nn0,1,2,3,.n谐振子的配分函数为1 2 1 0n Ze固体的内能和热容量分别为13lnUNZ, V VUCT 程序 ex4512 clear syms Z1 beta n hbar w U N k T Cv; Z1simplifysymsumexp-beta*hbar*w*n1/2,n,0,inf Usimplify-3*N*difflogZ1,beta beta1/k/T; U1subsU CvsimplifydiffU1,T 运行结果Z1 1/-1expbeta*hbar*w*exp1/2*beta*hbar*w U 3/2*N*hbar*w*expbeta*hbar*w1/-1expbeta*hbar*w U1 3/2*N*hbar*w*exp1/k/T*hbar*w1/-1exp1/k/T*hbar*w Cv 3*N*hbar2*w2*exp1/k/T*hbar*w/-1exp1/k/T*hbar*w2/k/T2 即1 21ee1Z, 3e12e1kTkTNU, 22e3e1kTVkTCNkkT( 3)德拜模型1917 年,德拜完成了他的固体热容量理论，他把固体看成连续介质，认为原子的 振动形成各种简正频率的弹性驻波，而把整个固体原子的微振动看作这些弹性驻波第 4 章 MATLAB在热物理学中的应用201 的叠加，每一个简正频率的弹性波的能量与同一频率简谐振子的能量是一样的。而弹性波又可分为纵波和横波，并且纵波和横波的波速均为一常数。根据这一思想，德拜从固体中原子振动的频率着手，得出固体的内能和定容热容量分别为3/30d9e1DTx DTxxUNkT；4/320ed9e1DxTVx DTxxCNk。其中，DxkTT, D称为德拜频率。 德拜的理论在低温区与实验符合得相当好，与实验发现的低温下热容量与T 3成正比的规律相一致，因此被称为德拜T 3律。 题目 ex4513 绘制杜隆珀替定律、 爱因斯坦模型和德拜模型的固体热容量随温度变化曲线，并讨论其在高、低温两端的性质。 解题分析 杜隆 -玻替定律113VCyNk 爱因斯坦模型令ExkTT，可将爱因斯坦固体热容量表达式改写为22 222ee3e1e1EExT VExTCyxNkT 德拜模型令Dxk TT 将德拜理论中热容量的表达式34/20ed9e1DxTVx DTxxCNk改写为3344/32200ed1ed333e1e1DyyTxV yy DCTyyyyyNkx202 MATLAB 及其在大学物理课程中的应用下面，采用数值方法计算上述积分， 程序 ex4513 clear,clf x00.011.3; y11; 杜隆珀替定律y21./x.2.*exp1./x./exp1./x-1.2; 爱因斯坦模型的热容量i0; 以下采用循环语句计算德拜模型的数值积分for x10.76920.5100 ii1; aiquadlexpy.*y.4./expy-1.2,0.001,x1; 德拜模型的热容量y3iai.*3./x1.3; end x10.76920.5100; plotx,y1,k-,x,y2,.r-,1./x1,y3,-bo axis0,1.3,0,1.1, xlabelT/theta, ylabelCv/3Nk legend杜隆珀替定律 ,爱因斯坦模型 ,德拜模型 从图 4-5-1 可知，在高温端，爱因斯坦模型和德拜模型的曲线都趋近于杜隆 -玻替定律，说明经典理论是量子理论的高温(或低频)近似。 实验表明， 在低温端， 爱因斯坦的热容量曲线比实验曲线要平缓一些 , 而德拜模型的热容量在低温端随温度的变化要比爱因斯坦模型来的快，与温度的三次方成比例,因此比爱因斯坦模型更符合实验结果。4.5.2 顺磁性固体的热力学性质顺磁性固体的理论模型是，磁性离子定域在晶体的特定格点上，认为离子间彼此相距甚远，相互作用可略去不计。因此，顺磁性固体是由定域、近独立的磁性图 4-5-1 固体热容量三种理论结果的比较第 4 章 MATLAB在热物理学中的应用203 离子组成的系统，遵从玻耳兹曼分布。1 顺磁体的热力学性质 题目 ex4521 计算顺磁体的磁化强度、内能和熵。 解题分析假定磁性离子的总角动量量子数为12，磁矩大小为2em 其中， 在外场中的能量的可能值为 B(磁矩沿外磁场方向)和 B(磁矩逆外磁场方向)， B 为外磁场的磁感应强度。由此，磁性离子的能量为BB离子的配分函数为1eeeBBZ磁化强度1lnNmZBB内能1lnUNZ熵11lnlnSNkZZ 程序 ex4521 ex45211磁化强度syms Z1 beta T k mu N B Z1expbeta*mu*Bexp-beta*mu*B; msimplifyN./beta.*difflogZ1,B 运行结果mN*mu*expbeta*mu*B-exp-beta*mu*B/expbeta*mu*Bexp-beta*mu*B 图 4-5-2 磁化强度曲线204 MATLAB 及其在大学物理课程中的应用即eeeeBBBBmN令 x B，y1 m / N，绘制 x-y1曲线。clf x-30.013; y1expx-exp-x./expxexp-x; plotx,y1,r-, xlabelbetamuB,ylabelm/muN grid on 运行结果如图4-5-2 所示。 ex45212内能Z1expbeta*mu*Bexp-beta*mu*B; U-N.*difflogZ1,beta 运行结果U-N*mu*B*expbeta*mu*B-mu*B*exp-beta*mu*B/expbeta*mu*B exp-beta*mu*B 即eeeeBBBBUNB令 x B，y2 U / NkT，绘制 x-y2曲线。syms x y; x-30.013; y2-x.*expx-exp-x./expxexp-x; plotx,y2,r- xlabelvB/kT ylabelU/NkT grid on 运行结果如图4-5-3 所示。 ex45213熵syms N k beta mu B z1expbeia*mu*Bexp-beia*mu*B; 图 4-5-3 内能与外磁场中能量B 的关系第 4 章 MATLAB在热物理学中的应用205 SN*k*logz1-beta*difflogz1,beta SN*k*logexpbeta*mu*Bexp-beta*mu*B-beta*mu*B*expbeta*mu*B- mu*B*exp-beta*mu*B/expbeta*mu*Bexp-beta*mu*B 即eelneeeeBB BB BBSNkNkB令 x B，y3 S / Nk，绘制x-y3曲线。x-30.013; y3logexpxexp-x-x.* expx-exp-x./expxexp-x; plotx,y3 xlabelvB/kT ylabelS/Nk grid 读者可根据玻耳兹曼关系来分析图 4-5-4 的物理意义。4.5.3 负温度状态 题目 ex4531 设核自旋量子数为12,在外磁场 B 下由于磁矩可与外磁场逆向或同向，其能量有两个可能值 2BeM，简记为。以 N 表示系统所含的总核磁矩数，x 和 y 分别表示能量为和的核磁矩数。 求系统的熵和内能之间的关系并绘制U-S曲线，由此讨论负温度状态及其意义。 解题分析由热力学知下述关系成立1xS TU在一般系统中熵随内能单调增加，因此T 为正。但是，也存在一些系统，当内能增加时熵反而减小，此时系统处在负温度状态，本题目就是这样的
如果觉得《ktt算法 约化_matlab在热学中的应用》对你有帮助，请点赞、收藏，并留下你的观点哦！
