科研课堂(二)——平均场近似和Ising模型

第一节课主要是叙述了“相”的概念以及两种不同的相边界模型(diffusive interface vs sharp interface)。第二节课的主题是相变(phase transition)。

相变是一个非平衡过程,描述起来比平衡过程要复杂。我们从最简单的Ising模型入手,对相变进行刻画。虽然数学模型较为简单,但是好处是可以定性地分析系统的特征温度、相变时对应的宏观物理量的变化。是研究复杂相变过程的基础。

Ising 模型最主要的假设有两个:

  • 状态只有两种 +1 -1 (更为复杂的情况,状态可以描述为矢量等连续变化的量、甚至可以描述为波函数)
  • 假设粒子只与最近的粒子发生相互作用(即只存在近场相互作用,无长程作用)

基于这两点出发,可以给出自旋粒子的哈密顿量。从哈密顿量出发、结合Gibbs概率密度分布(概率依能量指数分配),可推倒出粒子自旋量的概率测度(tanh分布,slide-12)。得到自旋的概率分布,就可以推出一系列的宏观物理量。如net magnetization、能量、磁滞回线等。

因此问题的核心就在于如何解出自旋的概率分布?平均场近似假设:每个粒子所受到的外界作用(课堂上的例子是外界磁场作用)是周围所有粒子的平均。因此在公式(5)中,可以将\(q_i\)用\(<q>\)代替。经过这种方法处理,方程大为简化(不需要考虑周围粒子的反转等瞬态过程)由平均场近似可以得出一系列很好的结果:

  • 相变温度Tc
  • 相变磁场Hc
  • 磁滞回线
  • 区分一级相变二级相变
  • \(<q> ~ (T_c-T)^{0.5}\)

平均场近似是不精确的。Onsager给出了二维Ising模型的解析接,序参量的critical exponent值不是0.5,而是0.125 。[1]

给同学们留下的作业是

1、回去动手算一下\(\tanh (Ax+b) = x\)的解。我用的是直接迭代。如果用牛顿迭代是不是收敛性更好?可否用程序证明?

2、算一下\(T^* = 0\) 和\(T^* > 0)时候的磁滞回线,看一下与低温时候有何区别?

资料:

[第二节课课件] [计算不同温度下的] [计算磁滞回线的python程序]

python环境管理包:anaconda: https://www.anaconda.com/

python的IDE:spyder 可以在anaconda中下载

有用的包:

自行搜索手册、可用anaconda安装或者用pip安装

基本数学工具和数组:numpy

科学计算、求解方程: scipy (部分替代mathematica)

数据表: pandas (可替代excel)

绘图:matplotlib (可替代origin)

参考文献

[1] Onsager, Lars. “Crystal statistics. I. A two-dimensional model with an order-disorder transition.” Physical Review 65.3-4 (1944): 117.

理论力学驻课笔记(三)

日期:2021年9月14日

地点:一号楼


这节课进入第二章——刚体的平衡。

主要概念:

力偶,力对点的矩矢、力对轴的矩。其中在课本上我们注意到力对点的 矩矢 是矢量,用粗体表示,而力对轴的矩是带有正负的标量。后面我们注意,力对轴的矩,其实是力对点的矩矢在某一个方向上的投影。

问题的引入:从直升机的螺旋桨引入矩的概念。对于质点而言,力的作用点在质点上,因此是共点力系。而对于刚体而言,由于考虑到形状,力的作用点可以在不同的位置,因此不一定是共点力系。这时候就需要引入矩(moment)这一概念。矩在物理中是一个十分重要的概念。一般来说当一个物理量的作用效果与空间位置相关,就需要一个将空间包含的的物理量表示。将物理量位置矢量与物理量的乘积(product)作为矩。

例如:电学中,将电荷(\(q\))与电荷位置矢量(\(\mathbf{r}\))的乘积叫做电偶极矩。对于点电荷和多个独立的点电荷组成的集合而言,电偶极矩 \( \mathbf{P} =q\mathbf{r}\)。对于一个电荷密度分布为 \(\rho(\mathbf{r})\))而言,电偶极矩为体积分 \( \mathbf{P} =\int_v \rho(\mathbf{r}) \mathbf(r)/ d \text{v}\) ,是一个矢量。同时还存在高阶矩,如四极子(quadrapole) \( \mathbf{P} =\int_v \rho(\mathbf{r}) \otimes \mathbf(r) \mathbf(r)/ d \text{v}\) 为一个二阶张量。

力学中,位置与力的叉积也是一个矢量,被称为力矩(moment 或torque)。理论力学中,密度分布的零阶矩是质量大小,一阶矩刻画质心,二阶矩是转动惯量(moment of inertia,或者叫质量惯性矩)注意这里\(\\mathbf{r}^2\)取矢量点积。材料力学中对于某一截面可以对某一个轴取矩,叫做面积二次轴矩(second axial moment of area),又称面积惯性矩,或面积对某一轴的惯性矩。

动量对某一点取矩,叫做角动量(angular momentum)\(\mathbf{L}\times\mathbf{r}\),这里是叉积。


随后一个重要的概念是力对点的矩矢和力对轴的矩两者之间的关系,以及和坐标投影的关系。

为研究刚体的受力特点,引入一种特殊的力系:力偶。力偶对应的矩矢与参考点、或者说说坐标的选取无关。因此力偶可以一一对应一个力偶矩矢,而且这个力偶矩矢是一个自由矢量。力偶具有的性质:

  • 不能与力等效,力的矩矢与参考点相关,因此无法等效
  • 力偶可以在作用面随意转移或者移到平行平面,这是利用了力偶矩矢是自由矢量的特点
  • 保持力*力臂乘积不变,作用效果不变。此即表示力偶矩矢的模不变,效果不变

为了方便对作用在刚体上的力进行简化和合并,引入加减平衡原理(在原有力系基础上加上一组共轴反向力是等效力系)。最后可以推导出空间任意力系可以简化为主矢加主矩。主矢大小方向与简化中心无关,主矩大小与简化中心的选择有关。

我的思考:

(1)作用在刚体上的力,如果作用效果已知,需要几个未知量来确定这个力?或者说刚体上力的自由度是几?三维空间中的力,唯一的确定他需要6个参数:\((\mathbf{F},\mathbf{r})\)。

(2)力对轴的矩,是否可以选取与轴无关的任意坐标轴,写出解析表达式?

复变函数笔记(一)

时间:2021年9月9日

地点:学院路校区三号楼106

授课教师:赵寿根

内容:复变函数第一节课

首先赵老师通过三方面引入了复变函数课程

  • 大学生-硕士生-博士生-博士后
  • 工程-科学-数学-数学的重要性-工程数学的重要性
  • 复变函数的作用:量子力学

而后介绍了复变函数的发展历史:

Cardano:用复数表示某一类无法用实数表示出来的方程的根。

笛卡尔-高斯-柯西-雷曼发展了这些概念。

复变函数在中国的发展:熊庆来:熊庆来和华罗庚的故事。

然后进入第一章 复数和复变函数

虚数单位 \(i = \sqrt(-1)\)

复数的构造 \(z = x +y i\)

理论力学驻课笔记(二)

时间:2021年9月9日

地点:学院路校区一号楼206

内容:理论力学第二次课,主要讲述平衡问题的解法

第一部分讲约束和约束力。

什么是约束—> 引入自由体和非自由体的概念。

约束力:约束作用于非自由体的力。其他的力统称为主动力。约束力的特点是:约束力方向与限制运动的方向相反。

我的思考:约束的概念在拉格朗日力学里是核心概念。约束减少了了体系的自由度。约束可以用自由度作为变量的等式(完全约束)或者不等式(非完全)约束来表示。对于完全约束,对于求解的问题而言,可以作为拉格朗日乘子引入。

在工程中,我们会遇到各种各样的约束。这里介绍了最常见的集中约束,并介绍了他们的特点:

  • 柔索:只能承受拉力,方向沿着柔索,背离物体。
  • 光滑面约束:切线上自由。约束力方向在公切线
  • 光滑圆柱铰链约束:有销钉。分为固定铰支座(二维1自由度)连接铰链(二维1自由度,链接两个物体),径向轴承(限制垂直于轴承的面方向的位移),活动铰链(水平无约束)。
  • 光滑球铰链。我的思考:被光滑球铰连约束的物体,只能绕着一点进行三维旋转,只有两个自由度。
  • 二力构件:根据二力平衡条件,二力构件收力所在直线已定

受力分析实例。举了三个例子,在课件上有。注意:

  • 整体法,对整个物体进行受力分析,忽略内部相互作用
  • 利用作用力-反作用力大小相等方向相反
  • 合理利用几何法
  • 投影法(或解析法)需要选取合适的坐标系

我的思考:最后的例子,用正余弦定理也挺方便。

理论力学驻课笔记(一)

时间:2021年9月7日

地点:学院路校区一号楼206

内容:理论力学第一节课,绪论与静力学基本概念

第一节课的主要目的是对理论力学这门课程,给学生一个感性的认知。

(1)王琪老师从科学哲学的角度出发,引入了力学的概念。告诉学生们学西力学有什么意义。

例子:比萨斜塔与竹子的类比,长城的建造以及赵州桥。

个人理解:长城的建造,视频中显示用树枝和粘土分层灌注,其实可以引申出复合材料的概念。为何加上树枝之后更加稳固(取向)

(2)力学处理问题的一般方法:
力学模型 — 仿真 — 实验 –对比 的脉络。并用车辆、行走器进行举例。

(3)学习目的

最终目的是提高学习能力。

最后一点时间讲了力的概念。力的作用是使物体产生移动、转动、变形和流动。这里就可以引出理论力学与材料力学的区别。对于理论力学来说,之研究物体的移动转动,并不考虑物体的变形与流动。


下半节

第一章内容——质点的平衡

(1)平衡的概念。首先引入力学模型,即怎么把工程问题中的实物进行力学建模。是质点、质点系,还是刚体?质点——只考虑质量,也就是说别的影响因素相比于质量来说可以忽略不计。刚体——考虑了尺寸和姿态(我理解的姿态,旋转、矩等),而质点系则考虑系统各部分之间的相互作用。理论力学可对物体进行如上所述的抽象建模。具体应该选择何种抽象,应根据所关心的问题而定。即使是同一个系统、同一个物理过程,也会因为研究问题的不同,选择不同的模型。这一点是需要着重强调的。

为了更好地理解这一点,举了乒乓球、保龄球、铅球的例子。马格纳斯效应。

我的思考:保龄球是否可以说是马格纳斯效应的结果,还是跟地面摩擦的结果?

(2)力的效应:分为内效应和外效应。内效应,内力,会引起形变,是材料力学、弹性力学考虑的问题。

(3)一连串的概念,这是理论力学后续内容的基础。

  • 力系:{\(\mathbf{F}_1\), \(\mathbf{F}_2\)… }
  • 等效力系:效应相同的 力系
  • 合力:n=1的等效力系,注意不是矢量和,要考虑矩
  • 平衡力系:和0力系等效的力系
  • 共点力系:力矢的延长线共点
  • 汇交力系:力汇交
  • 以后还会介绍:力系=共点力系+附加力矩

力的合成,基础是二里合成的平行四边形法则,或者三角形法则

第一章只涉及共点力系的合成。介绍两种方法,几何法以及投影法。投影法即解析法。在采用解析法时候,采用不同的坐标轴的选取,可以简化问题。举例。

思考:非共点力如何合成?三个力平衡的时候如何进行合成?依据是什么?

共点力系的平衡条件:(1)封闭(2)\(F_{Ri} = 0\)

思考:空间共点力系的平衡方程的投影轴必须正交吗?

我的思考:这个问题让学生回忆线性代数的基础知识:什么是线性相关,什么是空间基?正交基有什么好处?内积?更远一点:协变基矢量和逆变基矢量等等。


总结:两小节的内容大概相当于理论力学静力学绪论以及第一章前两节的内容

分子动力学计算SrMoO3 and BaMoO3的晶体及热力学参数

Molecular dynamics simulation of crystal structure and heat capacity in perovskite-type molybdates SrMoO3 and BaMoO3

本项工作与TU Darmstadt化学系的Albert课题组合作,研究了SMO和BMO的晶格常数以及热力学参数(比热容)。

记得之前在洛桑参加ISAF的时候遇到某位做XRD的老教授,他告诉我现在专心分析某些体系的基础参数的工作越来越少。一方面是因为工作本身很trivial,另一方面是这类文章很难发在有影响力的杂志上,但是这项工作又很有用。总体来说属于吃力不讨好的事情,很少有人愿意去做。大概这项工作属于他提及的那种范畴。

这项工作的idea其实很简单。我们隔壁Alff组十几年一直在用PLD、BME生长金属氧化物薄膜。一般来说通过PLD和BME生常驻来的薄膜不会很厚。因为随着晶体的生长,不可避免地会存在缺陷,此外基体与薄膜的晶格参数也有所区别会产生misfit strain。而他们课题组的绝活是,他们长出的钼基钙钛矿金属氧化物薄膜特别厚,可以达到50微米以上(参见P. Salg et al. APL applied. 2019),但是其中的原理并没有研究的很深入。想研究清楚这些,第一步是对这些氧化物晶体的基本物理参数进行测定,首先是晶格常数。但是经过一番文献调研,我们发现,在不多的SMO和BMO晶格常数测定的文献中,不同组做的结果有所差别:

科研课堂(一)——初识相场模型

相场模型(Phase field model)主要是为了解决界面问题。思想来源最早可以追溯到Van der Waals在处理液体表面张力时候的非均匀界面假设,经过Landau&Lifshitz、Allen-Cahn、Cahn-Hillard发展,日趋完善。对于相场模型的表达式,在简单、低维、平衡的条件下,人们已经给出了解析解。但是对于复杂的边界条件以及非平衡状态下序参量的演化过程,计算量十分庞大。

随着计算机科学的不断发展,求解复杂边界条件下的相场方程成为可能。典型的例子是R. Kayabushi(1993)的晶枝生长、以及LQ Chen等的铁电畴变的算例。目前相场方法应用于众多领域,如

  • 铁电
  • 铁磁
  • 合金
  • 马氏体相变
  • 多相合金
  • 粉末冶金

等等。

那么,究竟什么是相场模型?对于一个界面,我们可以有两种处理方式:

diffusive interface v.s. sharp interface

对于气、液界面,Tolman [1] 和Ono Kondo [2] 证明,采用如下所示的自由能表示方式,最后得到的界面的序参量是间断变化的。

\(\psi(z) = \psi(z) [\rho(z), T]\) (1)

自由能取最低的条件为表面张力为0。


Van der Waals 提出,在界面处的自由能分布不是均匀的,而是一个跟组份比例有关的参数:

\(\Psi(z) = \psi(z) [\rho(z), T]\) – \psi(\rho^{A,B}, T) (2)

与公式(1)相比增加的一部分自由能被称为surface excess。注意这里自由能的分布是位置坐标的函数。

Landau采用相同的思想,将铁磁体自由能的分布写为局部序参量(如magnetization)梯度的函数。Cahn等将梯度引入任意非均匀系统。在此假设下,自由能密度可以写为

\(\psi(c,\nabla c) = \psi _0(c) + G|\nabla c|^{2} \) (3)

其中c为系统的序参量,\( \psi _0(c) \) 为局部自由能密度, \( |\nabla c|^{2} \)为因序参量在空间分布的不均匀所导致的自由能的增加。 公式(3)即为平衡状态对应的相场公式。

与其他处理非均匀系统的方法相比, 相场法有着以下优点::

  • 相场方法对于界面的追踪是自动的。通过将体系自由能朝着能量最小的方向演化,可以自动获得界面的变化。
  • 相场方法是处理非平衡状态过程有力工具。
  • 相场的相参数可以与其它外场,如温度场、电场等进行耦合。


参考文献

[1] Tolman, Richard C. “The effect of droplet size on surface tension.” The journal of chemical physics 17.3 (1949): 333-337.

[2] Ono, Syu, and Sohei Kondo. “Molecular theory of surface tension in liquids.” Structure of Liquids/Struktur der Flüssigkeiten. Springer, Berlin, Heidelberg, 1960. 134-280.

[3] Cahn, John W., and John E. Hilliard. “Free energy of a nonuniform system. I. Interfacial free energy.” The Journal of chemical physics 28.2 (1958): 258-267.