删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

平行板受限下伊辛普适类临界卡西米尔力的一种计算方法

本站小编 Free考研考试/2021-12-25

陈东艳, 苗兵
中国科学院大学材料科学与光电技术学院, 北京 100049
2018年4月8日 收稿; 2018年4月27日 收修改稿
基金项目: 国家自然科学基金(21774131, 21544007)资助
通信作者: 苗兵,E-mail:bmiao@ucas.ac.cn

摘要: 研究统计力学中伊辛普适类的临界卡西米尔力。针对临界点上平行板几何约束下的伊辛普适类,基于欧拉-麦克劳林求和法发展了一套计算临界卡西米尔力的有效方法。结果表明,在平行板几何约束下,伊辛普适类临界卡西米尔力表现为以标度指数为-3的方式随板间距离衰减的长程吸引力。
关键词: 临界卡西米尔力伊辛普适类平行板受限欧拉-麦克劳林求和
An effective method for calculating the critical Casimir force between two parallel plates for the Ising universality class
CHEN Dongyan, MIAO Bing
College of Materials Science and Opto-Electronic Technology, University of Chinese Academy of Sciences, Beijing 100049, China


Abstract: Due to the constraints introduced by boundaries for the long range fluctuations in a massless field, a long-range fluctuation-induced effective force between boundaries is induced, and it is called the Casimir force. The classical Casimir force refers to the interaction between two conducting plates immersed in a vacuum, and it is originated from the quantum fluctuations of the electromagnetic field. In this work, we study the critical Casimir force for the Ising universality class in statistical mechanics. Specifically, we develop an effective method, namely, the Euler-Maclaurin summation method, to calculate the critical Casimir force between two parallel plates at the critical point of the Ising universality class. Our results demonstrate that, in this geometric set-up, the critical Casimir force behaves as a long-range attractive force with a scaling exponent of -3 with respect to the boundary distance.
Keywords: critical Casimir forceIsing universality classparallel plates confinementEuler-Maclaurin summation
卡西米尔效应是1948年由著名荷兰物理学家Casimir[1]发现的。Casimir通过计算预测置于真空中的两块无穷大平行导体板,由于对真空量子涨落的边界约束,从而板间会诱导产生一个长程吸引力,即卡西米尔力,该现象称为卡西米尔效应,是真空零点能的一种物理体现。
由于卡西米尔力是一个微弱的力,所以在实验上进行验证比较困难。对卡西米尔效应的实验验证,最早来自于Sparnaay[2]在1958年进行的一个实验,Sparnaay利用弹簧巧妙地构造了一个电容器,力的变化体现为容易测量的电容量的变化,从而定性地从实验上验证了卡西米尔力的存在。随着实验技术的进步,人们越来越追求实验的精度,美国华盛顿大学的Lamoreaux[3]在1997年率先以足够的精度从实验上验证了卡西米尔力,实现了卡西米尔力的定量测量。Lamoreaux测量的是一个镀金圆球透镜与一块金属板之间的卡西米尔力,金属板与一个精准的扭秤相连,测量精度在5%~10%。1998年,Mohideen和Roy[4]使用激光技术更精确地测量了镀铝平板和小球之间的卡西米尔力,精度达到1%。
对经典卡西米尔效应的研究持续了半个多世纪,研究越来越细致和深入,有大量的研究成果。在理论研究方面,人们主要集中在研究不同的几何构型和边界条件下的卡西米尔力,详见综述性文献[5-9]。
20世纪70年代,人们发现只要在存在长程关联的涨落场中引入边界约束,边界间就会产生一个由涨落所诱导的力,这种力被称为广义的卡西米尔力。由真空涨落所引起的卡西米尔力只是其中的一种,另外如统计力学临界点处的密度场等,这些具有长程涨落特性的场都可以满足产生卡西米尔力的条件。1978年,Fisher和De Gennes[10]首次将经典的卡西米尔效应推广到热卡西米尔效应或者临界卡西米尔效应,涨落背景场由量子涨落替换为临界点处的长程密度涨落,开辟了临界卡西米尔效应研究的新领域。
近年来,在对临界卡西米尔效应的研究中,一个研究重点是在统计力学伊辛普适类中平行板几何约束下的卡西米尔力,包括平衡态的临界卡西米尔力和非平衡态的卡西米尔力[11]
本文研究伊辛普适类在平行板受限下的临界卡西米尔力,不同于常用的场论方法[12-13],我们基于欧拉-麦克劳林求和法(Euler-Maclaurin summation),发展了一套计算两块平行板之间临界卡西米尔力的有效方法,这种方法得到的结果与用场论方法得到的结果基本一致,并且简单易操作。计算结果表明,在平行板几何约束下,伊辛普适类的临界卡西米尔力表现为一种长程吸引力,与边界间距的标度指数为-3,该标度指数与标度分析的结果一致。
1 模型描述如图 1所示,考虑在伊辛普适类中,板间距离为a的平行板受限体系中的临界卡西米尔效应。方便起见,取热能kBT=1。在边界上,采用Dirichlet边界条件,即要求边界上密度场为零。该边界条件使得垂直于平板方向的涨落模式离散化:(qz=nπ/a),其中n=1, 2, …。涨落波矢大小为q2=|q//|2+qz2,其中q//表示平行于平板方向的二维波矢。注意在本文中所有的长度都已由体系的微观相互作用力程约化为无量纲量。
Fig. 1
Download: JPG
larger image
图 1 板间距离为a的平行板受限体系 Fig. 1 Confinement set-up of parallel plates with the distance a
图 1 板间距离为a的平行板受限体系

Fig. 1 Confinement set-up of parallel plates with the distance a -->

由统计力学可知,伊辛普适类在临界点上的高斯涨落自由能写成如下形式
$\frac{F}{A} = \frac{1}{2}\int {\frac{{{{\rm{d}}^2}{\mathit{\boldsymbol{q}}_{//}}}}{{{{(2\pi )}^2}}}} \sum\limits_{n = 1}^\infty {\ln } \left[ {{{\left( {\frac{{n\pi }}{a}} \right)}^2} + q_{//}^2} \right], $ (1)
式中:A是平行板的面积。显然,这是一个无穷大的量,为了对其进行正规化处理,仿照量子场论的做法[14],引入正规化因子$ {{\rm{e}}^{ - \ln \left[ {{{\left( {\frac{{n\pi }}{a}} \right)}^2} + q_{//}^2} \right]/\Lambda }}$,这里,Λ是截断常数。完成所有计算后,将令Λ→∞。此时,自由能写成
$\frac{F}{A}=\lim\limits_{x \rightarrow 0} \frac{F(x)}{A}.$ (2)
这里定义x=1/Λ,并且
$ \begin{array}{*{20}{l}}\frac{{F(x)}}{A} = \frac{1}{2}\int {\frac{{{{\rm{d}}^2}{\mathit{\boldsymbol{q}}_{//}}}}{{{{(2\pi )}^2}}}} \sum\limits_{n = 1}^\infty {\ln } \left[ {{{\left( {\frac{{n\pi }}{a}} \right)}^2} + q_{//}^2} \right]{{\rm{e}}^{ - x\ln \left[ {{{\left( {\frac{{n\pi }}{a}} \right)}^2} + q_{//}^2} \right]}}\\ = - \frac{1}{2}\frac{\partial }{{\partial x}}\int {\frac{{{{\rm{d}}^2}{\mathit{\boldsymbol{q}}_{//}}}}{{{{(2\pi )}^2}}}} \sum\limits_{n = 1}^\infty {{{\rm{e}}^{ - x\ln \left[ {{{\left( {\frac{{n\pi }}{a}} \right)}^{2 + }}q_{//}^2} \right]}}} \\{ = - \frac{1}{2}\frac{\partial }{{\partial x}}\int {\frac{{{{\rm{d}}^2}{\mathit{\boldsymbol{q}}_{//}}}}{{{{(2\pi )}^2}}}} \sum\limits_{n = 1}^\infty {{{\left[ {{{\left( {\frac{{n\pi }}{a}} \right)}^2} + q_{//}^2} \right]}^{ - x}}} }\\{ = - \frac{1}{{8\pi }}\frac{\partial }{{\partial x}}\left[ {\frac{1}{{x - 1}}{{\left( {\frac{\pi }{a}} \right)}^{ - 2(x - 1)}}\sum\limits_{n = 1}^\infty {{n^{ - 2(x - 1)}}} } \right].}\end{array}$ (3)
其中,在对平行波矢的积分结果中,舍去了与无穷大有关的项,因为最后计算卡西米尔力的涨落自由能只需要保留有限大小的项。
2 欧拉-麦克劳林求和法为了计算上式中的无穷级数求和,采用欧拉-麦克劳林求和法。欧拉-麦克劳林求和公式是一个将积分与求和联系起来的式子,被数学家和物理学家广泛应用,其具体表达式(只取到3阶导数)如下
$\begin{aligned} \int_{a}^{b} f(t) \mathrm{d} t=&\left[\frac{f(a)}{2}+f(a+h)+\cdots+\frac{f(b)}{2}\right] h+\\ & \frac{1}{12}\left[f^{\prime}(a)-f^{\prime}(b)\right]-\\ & \frac{1}{720}\left[f^{\prime \prime\prime}(a)-f^{\prime \prime\prime}(b)\right] .\end{aligned}$ (4)
式中:f(t)表示一个可微函数,ab是积分限常数,h是离散步长。
定义函数f(t)=t-2(x-1), 取离散步长h=1,根据欧拉-麦克劳林公式有
$\begin{aligned} \int_{1}^{\infty} t^{-2(x-1)} \mathrm{d} t=& \frac{1^{-2(x-1)}}{2}+2^{-2(x-1)}+\cdots+\infty+\\ & \frac{1}{12} f^{\prime}(1)-\frac{1}{720} f^{\prime \prime\prime}(1) \\=& \sum\limits_{n=1}^{\infty} n^{-2(x-1)}-\frac{1}{2}+\frac{1}{12} f^{\prime}(1)-\\ & \frac{1}{720} f^{\prime\prime \prime}(1)+\infty. \end{aligned}$ (5)
其中,
$\begin{aligned} f^{\prime}(t) &=-2(x-1) t^{-2 x+1}, \\ f^{\prime \prime\prime}(t) &=-4 x(x-1)(2 x-1) t^{-2 x-1}, \\ f^{\prime}(1) &=-2(x-1), \\ f^{\prime\prime \prime}(1) &=-4 x(x-1)(2 x-1).\end{aligned}$ (6)
这里,将所有与无穷大有关的项归并,用∞表示。简单的运算可得
$\begin{aligned} \sum\limits_{n=1}^{\infty} n^{-2(x-1)}=& \int_{1}^{\infty} t^{-2(x-1)} \mathrm{d} t+\frac{1}{2}+\frac{1}{6}(x-1)-\\ & \frac{1}{180} x(x-1)(2 x-1) \\&= \infty+\frac{1}{2}+\frac{1}{6}(x-1)-\\ & \frac{1}{180} x(x-1)(2 x-1)+\frac{1}{2 x-3} \end{aligned}$ (7)
同样,所有与无穷大有关的项已经用∞表示。进一步,得到
$\mathop {\lim }\limits_{x \to 0} \sum\limits_{n = 1}^\infty {{n^{ - 2(x - 1)}}} = 0.$ (8)
这里,类似于量子场论中由无穷的真空零点能计算卡西米尔力的方法,已经舍去无穷大,只取了有限部分。
我们注意到黎曼函数的级数定义为
$\zeta (z) = \sum\limits_{n = 1}^\infty {{n^{ - z}}} , $ (9)
因此易知,需要求算的无穷级数可用黎曼函数表示如下
$\sum\limits_{n = 1}^\infty {{n^{ - 2(x - 1)}}} = \zeta (2(x - 1)).$ (10)
由黎曼函数的性质有
$\mathop {\lim }\limits_{x \to 0} \sum\limits_{n = 1}^\infty {{n^{ - 2(x - 1)}}} = \mathop {\lim }\limits_{x \to 0} \zeta (2(x - 1)) = \zeta ( - 2) = 0.$ (11)
显然,这一结果与我们用欧拉-麦克劳林求和法得到的结果(8)一致。
由式(7),进一步可得
$\begin{array}{*{20}{c}}{\frac{\partial }{{\partial x}}\sum\limits_{n = 1}^\infty {{n^{ - 2(x - 1)}}} = \frac{1}{6} - \frac{1}{{30}}\left( {{x^2} - x + \frac{1}{6}} \right) - }\\{\frac{1}{{2{{(x - 3/2)}^2}}}.}\end{array}$ (12)
同样,只保留有限部分。由此可得
$\mathop {\lim }\limits_{x \to 0} \frac{\partial }{{\partial x}}\sum\limits_{n = 1}^\infty {{n^{ - 2(x - 1)}}} = - \frac{{11}}{{180}}.$ (13)
综合以上结果,得到涨落自由能如下
$\begin{array}{l}\frac{F}{A} = \mathop {\lim }\limits_{x \to 0} \frac{{F(x)}}{A}\\ = \mathop {\lim }\limits_{x \to 0} \left[ {\frac{1}{{8\pi }}\left[ {\begin{array}{*{20}{c}}{\frac{1}{{{{(x - 1)}^2}}}{{\left( {\frac{\pi }{a}} \right)}^{ - 2(x - 1)}}}\\{ + \frac{2}{{x - 1}}{{\left( {\frac{\pi }{a}} \right)}^{ - 2(x - 1)}}\ln \left( {\frac{\pi }{a}} \right)}\end{array}} \right]\sum\limits_{n = 1}^\infty {{n^{ - 2(x - 1)}}} } \right] + \\\mathop {\lim }\limits_{x \to 0} \left[ { - \frac{1}{{8\pi }}\frac{1}{{x - 1}}{{\left( {\frac{\pi }{a}} \right)}^{ - 2(x - 1)}}\frac{\partial }{{\partial x}}\sum\limits_{n = 1}^\infty {{n^{ - 2(x - 1)}}} } \right]\\ = - \frac{{11\pi }}{{1440}}\frac{1}{{{a^2}}} \approx - 0.024{a^{ - 2}}.\end{array}$ (14)
该结果与场论结果[13]基本一致。
进一步,临界卡西米尔力由涨落自由能对边界距离的导数得到,即
$\begin{aligned} f(a) &=-\frac{\partial}{\partial a}(F / A)=-\frac{11 \pi}{720 a^{3}} \\ & \approx-0.048 a^{-3}. \end{aligned}$ (15)
由该结果可知,对于伊辛普适类,平行板受限下的临界卡西米尔力表现为一个长程吸引力,并且与板间距离的标度指数为-3。
3 标度分析本小节对临界卡西米尔力的标度指数做简单的标度分析。
这里定义的临界卡西米尔力表现为两平行板之间的分离压(disjoining pressure),因此有[f]=E/L3, 这里EL分别代表能量和长度。注意,能量已由热能约化。在伊辛普适类中,临界点上涨落的关联长度发散,因此不允许存在一个有限大小的关联长度进入临界卡西米尔力的表达式。体系中只有两平行板之间的距离是可以进入最后结果的特征长度,因此,立即得到fa-3, 这正是本文得到的标度指数。需要指出的是虽然标度指数可以通过标度分析得到,然而临界卡西米尔力在特定几何约束和边界条件下表现为吸引或排斥,只能通过具体计算得到。
4 结论本文应用欧拉-麦克劳林求和公式对卡西米尔效应进行研究,对平行板受限下伊辛普适类的临界卡西米尔力进行解析推导,基于欧拉-麦克劳林求和公式与黎曼函数的相关知识发展出一套计算临界卡西米尔力的有效方法。此方法对研究其他体系中的卡西米尔效应具有重要的借鉴意义。在未来的工作中,我们将尝试将该方法用于其他受限几何,以及其他统计力学普适类中临界卡西米尔力的计算。
参考文献
[1] Casimir H B G. On the attraction between two perfectly conducting plates[J]. Proc Kon Ned Akad Wet, 1948, 51: 793-795.
[2] Sparnaay M J. Measurements of attractive forces between flat plates[J]. Physica, 1958, 24(6-10): 751-764. DOI:10.1016/S0031-8914(58)80090-7
[3] Lamoreaux S K. Demonstration of the Casimir force in the 0.6 to 6μm range[J]. Physical Review Letters, 1997, 78(1): 5-8. DOI:10.1103/PhysRevLett.78.5
[4] Mohideen U, Roy A. Precision measurement of the Casimir force from 0.1 to 0.9μm[J]. Physical Review Letters, 1998, 81(21): 4549-4552. DOI:10.1103/PhysRevLett.81.4549
[5] Plunien G, Müller B, Greiner W. The casimir effect[J]. Physics Reports, 1986, 134(2/3): 87-193.
[6] Elizalde E, Romeo A. Essentials of the Casimir effect and its computation[J]. American Journal of Physics, 1991, 59: 711-719. DOI:10.1119/1.16749
[7] Kardar M, Golestanian R. The "friction" of vacuum, and other fluctuation-induced forces[J]. Reviews of Modern Physics, 1999, 71(4): 1233-1245. DOI:10.1103/RevModPhys.71.1233
[8] Bordag M, Klimchitskaya G L, Mohideen U, et al. Advances in the Casimir effect[M]. Oxford: OUP Oxford, 2009.
[9] Lamoreaux S K. The Casimir force:background, experiments, and applications[J]. Reports on Progress in Physics, 2004, 68(1): 201-236.
[10] Fisher M E, De Gennes P G. Wall phenomena in a critical binary mixture[J]. Comptes Rendus Hebdomadaires Des Seances De L Academie Des Sciences Serie B, 1978, 287(8): 207-209.
[11] Aminov A, Kafri Y, Kardar M. Fluctuation-induced forces in nonequilibrium diffusive dynamics[J]. Physical Review Letters, 2015, 114(23): 230602. DOI:10.1103/PhysRevLett.114.230602
[12] Uchida N. Casimir effect in fluids above the isotropic-lamellar transition[J]. Physical Review Letters, 2001, 87(21): 216101. DOI:10.1103/PhysRevLett.87.216101
[13] Li H, Kardar M. Fluctuation-induced forces between manifolds immersed in correlated fluids[J]. Physical Review A, 1992, 46(10): 6490-6500. DOI:10.1103/PhysRevA.46.6490
[14] Zee A. Quantum field theory in a nutshell[M]. Princeton: Princeton university press, 2010.


相关话题/计算 实验 测量 力学 统计

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 展向磁场作用下液态金属GaInSn多层膜流实验研究
    齐天煜1,阳倦成2,倪明玖11.中国科学院大学工程科学学院,北京100049;2.西安交通大学航天航空学院机械结构与强度国家重点实验室,西安7100492018年2月27日收稿;2018年4月10日收修改稿基金项目:国家自然科学基金(51406193,51636009)和科技部项目(2013GB11 ...
    本站小编 Free考研考试 2021-12-25
  • 地球椭球模型中太阳位置计算的改进
    李文1,2,3,赵永超1,2,31.中国科学院空间信息处理与应用系统技术重点实验室,北京100190;2.中国科学院电子学研究所,北京100190;3.中国科学院大学,北京1000492018年1月25日收稿;2018年3月29日收修改稿基金项目:国家重大科研仪器研制项目(41427805)资助通信 ...
    本站小编 Free考研考试 2021-12-25
  • 小麦炭化过程中质量和颜色变化的实验模拟
    苏鑫,李玉梅,谷永建中国科学院大学地球与行星科学学院,北京100049;中国科学院计算地球动力学重点实验室,北京100049;中国科学院地质与地球物理研究所,北京100029;中国科学院大学测试中心分子化石实验室,北京1000492018年2月9日收稿;2018年4月18日收修改稿基金项目:国家自然 ...
    本站小编 Free考研考试 2021-12-25
  • 参数可调椭圆滑动轴承转子系统动力学研究及抑振机理
    张磊1,裴世源1,徐华1,21.西安交通大学现代设计及转子轴承系统教育部重点实验室,西安710049;2.新疆大学机械工程学院,乌鲁木齐8300472017年8月31日收稿;2018年1月22日收修改稿基金项目:国家自然科学基金(51605367,51575421)资助通信作者:裴世源,E-mail ...
    本站小编 Free考研考试 2021-12-25
  • 强磁场对导电流体热毛细流动和换热影响的实验研究
    陈然,王增辉,倪明玖中国科学院大学工程科学学院,北京1000492017年12月8日收稿;2018年3月5日收修改稿基金项目:国家自然科学基金(51476162,51876201)资助通信作者:王增辉,E-mail:wzhawk@ucas.ac.cn摘要:核聚变装置限制器有效地屏蔽来自器壁的杂质,排 ...
    本站小编 Free考研考试 2021-12-25
  • 腹主动脉瘤形态学二维断面和三维立体测量的误差分析
    尉驰俊1,邱越1,袁丁2,郑庭辉11.四川大学建筑与环境学院力学系,成都610207;2.四川大学华西医院血管外科,成都6100002018年5月14日收稿;2018年7月10日收修改稿基金项目:国家自然基金(81500370)资助通信作者:郑庭辉,E-mail:tinghuizh@scu.edu. ...
    本站小编 Free考研考试 2021-12-25
  • 沉积环境下涡轮叶片前缘气膜冷却的实验研究*
    在燃气涡轮发动机运行的过程中,涡轮部件会受到污染物的沉积、腐蚀和侵蚀。当飞机在起降或低空飞行时,燃气涡轮发动机往往会吸入灰尘、砂砾或其他颗粒物,特别是在沙漠或火山活跃地区。此外,航空燃料在高温环境下生成的金属化合物的副产物,也是极易沉积的颗粒物。其中,钒燃烧的产物V2O3会在涡轮叶片上沉积,在过剩热 ...
    本站小编 Free考研考试 2021-12-25
  • 翼尖铰接复合飞行器动力学特性研究*
    作为一种特定的新概念飞行器系统,翼尖铰接复合飞行技术具体是以单架次分别进行起飞爬升和下降着陆,在巡航过程中通过特定设计的翼尖铰柱接合机构实现飞机复合飞行,这种技术能够有效降低飞机的油耗率,提升飞机的航程和航时等续航性能指标。对于远距离长航时飞机,这种技术的优点具体体现在:可以增加展弦比、翼型弯度,同 ...
    本站小编 Free考研考试 2021-12-25
  • 气冷涡轮导叶流热耦合计算及机理*
    涡轮内部的流场与温度场是高度耦合的,一方面涡轮内部流场存在着如激波、二次涡系等复杂流动结构,会影响到部件温度场分布及其换热特性,另一方面,固体域的温度分布及冷却射流与主流的相互作用也将反过来影响流场的流动特征[1]。因此,理解和掌握涡轮内部流场温度场相互作用机制,对于提高冷却效率和准确预测热负荷有着 ...
    本站小编 Free考研考试 2021-12-25
  • 轴向通流旋转盘腔流动换热稳态实验与数值模拟*
    当前,提高涡轮前温度是提升航空发动机综合性能的主要途径之一。为使涡轮等热端部件能够在愈发恶劣的环境下正常、持久地工作,需要从压气机级间抽取冷气对其进行冷却。冷却气体用量的增加可以提高冷却效果,但提高了主流损失,对整机性能影响巨大;冷气由压气机级间进入盘腔后会与压气机转盘和盘罩换热,形成复杂的流动结构 ...
    本站小编 Free考研考试 2021-12-25