搜索

文章查询

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

温稠密铝等离子体物态方程及其电离平衡研究

王天浩 王坤 张阅 姜林村

温稠密铝等离子体物态方程及其电离平衡研究

王天浩, 王坤, 张阅, 姜林村
PDF
HTML
导出引用
导出核心图
  • 温稠密物质的物性参数在惯性约束聚变能源、Z箍缩等高能量密度物理领域的实验结果分析和物理过程数值模拟等方面有着重要的应用价值. 本文应用部分电离等离子体模型, 在理想自由能的基础上考虑了库仑相互作用、排斥体积作用和极化作用等非理想特性, 开展了温稠密等离子体物态方程和电离平衡的研究. 计算了温稠密铝等离子体的压强等物态方程数据和在密度为1.0 × 10–4—3.0 g/cm3, 温度为1.0 × 104—3.0 × 104 K范围内的粒子组分. 计算结果显示, 铝等离子体的平均电离度在临界密度区域内随着密度的增加而突然增大. 根据非理想Saha方程中有效电离能这一关键参数, 分析了铝等离子体平均电离度在临界密度区域内随密度迅速增大的现象.
      通信作者: 王坤, kunwang@hebut.edu.cn
    • 基金项目: 国家自然科学基金(批准号: 51807050)、河北省自然科学基金(批准号: E2019202297)和河北省高等学校青年拔尖人才项目(批准号: BJ2017038)资助的课题
    [1]

    付志坚, 贾丽君, 夏继宏, 唐可, 李召红, 权伟龙, 陈其峰 2016 物理学报 65 065201

    Fu Z J, Jia L J, Xia J H, Tang K, Li Z H, Quan W L, Chen Q F 2016 Acta Phys. Sin. 65 065201

    [2]

    Wallace M S, Haque S, Neill P, Pereira N R, Presura R 2018 Rev. Sci. Instrum. 89 015106

    [3]

    Oreshkin V I, Artyomov A P, Chaikovsky S A, Oreshkin E V, Rousskikh A G 2017 Phys. Plasmas 24 012703

    [4]

    格拉齐亚尼F, 德斯贾莱斯M P, 雷德默R, 特里基S B 著 (陈其峰 译) 2018 温稠密物质研究的前沿和挑战 (北京: 原子能出版社) 第vii−x页

    Graziani F, Desjarlais M P, Redmer R, Trickey S B (translated by Chen Q F) 2018 Frontiers and Challenges in Warm Dense Matter (Beijing: Atomic Energy Press) ppvii−x (in Chinese)

    [5]

    徐锡申, 张万箱 1986 实用物态方程理论导引 (北京: 科学出版社) 第1−5页

    Xu X S, Zhang W X 1986 Introduction to the Theory of Equation of State (Beijing: Science Press) pp1−5 (in Chinese)

    [6]

    Shi Z Q, Wang K, Li Y, Shi Y J, Wu J, Jia S L 2014 Phys. Plasmas 21 032702

    [7]

    Eliezer S, Ghatak A, Hora H, Teller E 2002 Fundamentals of Equations of State (Singapore: World Scientific) pp153−164

    [8]

    陈其峰, 顾云军, 郑君, 李江涛, 李治国, 权伟龙, 付志坚, 李成军 2017 科学通报 62 812

    Chen Q F, Gu Y J, Zheng J, Li J T, Li Z G, Quan W L, Fu Z J, Li C J 2017 Chin. Sci. Bull. 62 812

    [9]

    于继东, 李平, 王文强, 吴强 2014 物理学报 63 116401

    Yu J D, Li P, Wang W Q, Wu Q 2014 Acta Phys. Sin. 63 116401

    [10]

    Danel J F, Kazandjian L, Zérah G 2008 Phys. Plasmas 15 072704

    [11]

    Fu Z J, Quan W L, Zhang W, Li Z G, Zheng J, Gu Y J, Chen Q F 2017 Phys. Plasmas 24 013303

    [12]

    汤文辉, 徐彬彬, 冉宪文, 徐志宏 2017 物理学报 66 030505

    Tang W H, Xu B B, Ran X W, Xu Z H 2017 Acta Phys. Sin. 66 030505

    [13]

    张颖, 陈其峰, 顾云军, 蔡灵仓, 卢铁城 2007 物理学报 56 1318

    Zhang Y, Chen Q F, Gu Y J, Cai L C, Lu T C 2007 Acta Phys. Sin. 56 1318

    [14]

    Chen Q F, Cai L C, Gu Y J, Gu Y 2009 Phys. Rev. E 79 016409

    [15]

    Chen Q F, Zheng J, Gu Y J, Chen Y L, Cai L C 2011 Phys. Plasmas 18 112704

    [16]

    Quan W L, Chen Q F, Fu Z J, Sun X W, Zheng J, Gu Y J 2015 Phys. Rev. E 91 023106

    [17]

    Apfelbaum E M 2015 Phys. Plasmas 22 092703

    [18]

    Apfelbaum E M 2017 High Temp. 55 1

    [19]

    Kuhlbrodt S, Holst B, Redmer R 2005 Contrib. Plasma Phys. 45 73

    [20]

    Moldabekov Z A, Groth S, Dornheim T, Bonitz M, Ramazanov T S 2017 Contrib. Plasma Phys. 57 532

    [21]

    Stolzmann W, Blöcker T 1996 Astron. Astrophys. 314 1024

    [22]

    Fortov V E, Altshuler L V, Trunin R F, Funtikov A I 2004 High-Pressure Shock Compression of Solids VII (New York: Springer) pp437−489

    [23]

    付志坚, 陈其峰, 陈向荣 2011 物理学报 60 055202

    Fu Z J, Chen Q F, Chen X R 2011 Acta Phys. Sin. 60 055202

    [24]

    Redmer R, Rother T, Schmidt K, Kraeft W D, Röpke G 1988 Contrib. Plasma Phys. 28 41

    [25]

    Apfelbaum E M 2011 Phys. Rev. E 84 066403

    [26]

    Karakhtanov V S, Redmer R, Reinholz H, Röpke G 2011 Contrib. Plasma Phys. 51 355

    [27]

    More R M, Warren K H, Young D A, Zimmerman G B 1988 Phys. Fluids 31 3059

    [28]

    Iyetomi H, Ichimaru S 1986 Phys. Rev. A 34 433

    [29]

    Kresse G, Hafner J 1993 Phys. Rev. B 47 558

    [30]

    Zaghloul M R 2018 High Energy Density Phys. 26 8

    [31]

    Wu B, Shin Y C 2006 Appl. Phys. Lett. 89 111902

    [32]

    Morel V, Bultel A, Chéron B G 2009 Int. J. Thermophys. 30 1853

    [33]

    Preston T R, Vinko S M, Ciricosta O, Chung H K, Lee R W, Wark J S 2013 High Energy Density Phys. 9 258

    [34]

    Stransky M 2016 Phys. Plasmas 23 012708

    [35]

    Kemp A J, Meyer-Ter-Vehn J 1998 Nucl. Instrum. Methods Phys. Res., Sect. A 415 674

    [36]

    Krisch I, Kunze H J 1998 Phys. Rev. E 58 6557

    [37]

    田庆云 2015 硕士学位论文 (长沙: 国防科学技术大学)

    Tian Q Y 2015 M. S. Thesis (Changsha: National University of Defense Technology) (in Chinese)

    [38]

    Son S K, Thiele R, Jurek Z, Ziaja B, Santra R 2014 Phys. Rev. X 4 031004

    [39]

    Ebeling W, Föster A, Fortov V E, Gryaznov V K, Polishuk A Y 1991 Thermophysical Properties of Hot Dense Plasmas (Stuttgar-Leipzig: Teubner Verlagsgesellshaft) pp39—42

  • 图 1  在密度为0.1 g/cm3时, 不同模型计算的铝等离子体压强随温度的变化

    Fig. 1.  The pressure of aluminum plasma calculated by different models as a function of temperature at density of 0.1 g/cm3.

    图 2  铝等离子体的平均电离度随密度与温度变化等值线分布

    Fig. 2.  Contour map of average ionization degree of aluminum plasma as a function of density and temperature.

    图 3  等离子体中不同粒子的相对粒子分数在密度为0.1 g/cm3时随温度的变化

    Fig. 3.  Dependence of relative particle fraction of different particles on temperature at density of 0.1 g/cm3.

    图 4  不同模型计算的铝等离子体平均电离度在不同温度下随密度的变化 (a) 10000 K; (b) 15000 K

    Fig. 4.  Average ionization degree of aluminum plasma calculated by different models as a function of density at different temperatures: (a) 10000 K; (b) 15000 K.

    图 5  不同温度下非理想效应自由能密度以及电子、原子相对粒子分数随密度的变化 (a) 库仑相互作用; (b) 排斥体积作用; (c) 极化作用; (d) 电子、原子相对粒子分数

    Fig. 5.  Free energy density of different non-ideal effects and relative particle fraction for electrons and atoms as a function of density at different temperatures: (a) Coulomb interaction; (b) excluded volume effect; (c) polarization effect; (d) relative particle fraction for electrons and atoms.

    图 6  温度为15000 K时等离子体中不同粒子化学势的非理想部分随密度的变化

    Fig. 6.  Dependence of non-ideal chemical potential of particles on density at temperature of 15000 K.

    图 7  温度为15000 K时不同模型计算的电离能的降低∆E随密度的变化(黑色曲线, 非理想Saha方程; 红色曲线, DmEK模型)

    Fig. 7.  Depression of ionization potential calculated by different models as a function of density at 15000 K. Black lines correspond to nonideal Saha equation; red lines correspond to DmEK model.

    表 1  物态方程及电离平衡模型参数列表

    Table 1.  Parameters of equation of state and ionization equilibrium model.

    参数参数参数
    R02.5714R60.4167E32.0909
    R12.3377rc1.7712E48.8192
    R22.1429mk24597E511.3059
    R30.4678αD46.281E614.0008
    R40.4494E10.4400
    R50.4324E21.3839
    下载: 导出CSV
  • [1]

    付志坚, 贾丽君, 夏继宏, 唐可, 李召红, 权伟龙, 陈其峰 2016 物理学报 65 065201

    Fu Z J, Jia L J, Xia J H, Tang K, Li Z H, Quan W L, Chen Q F 2016 Acta Phys. Sin. 65 065201

    [2]

    Wallace M S, Haque S, Neill P, Pereira N R, Presura R 2018 Rev. Sci. Instrum. 89 015106

    [3]

    Oreshkin V I, Artyomov A P, Chaikovsky S A, Oreshkin E V, Rousskikh A G 2017 Phys. Plasmas 24 012703

    [4]

    格拉齐亚尼F, 德斯贾莱斯M P, 雷德默R, 特里基S B 著 (陈其峰 译) 2018 温稠密物质研究的前沿和挑战 (北京: 原子能出版社) 第vii−x页

    Graziani F, Desjarlais M P, Redmer R, Trickey S B (translated by Chen Q F) 2018 Frontiers and Challenges in Warm Dense Matter (Beijing: Atomic Energy Press) ppvii−x (in Chinese)

    [5]

    徐锡申, 张万箱 1986 实用物态方程理论导引 (北京: 科学出版社) 第1−5页

    Xu X S, Zhang W X 1986 Introduction to the Theory of Equation of State (Beijing: Science Press) pp1−5 (in Chinese)

    [6]

    Shi Z Q, Wang K, Li Y, Shi Y J, Wu J, Jia S L 2014 Phys. Plasmas 21 032702

    [7]

    Eliezer S, Ghatak A, Hora H, Teller E 2002 Fundamentals of Equations of State (Singapore: World Scientific) pp153−164

    [8]

    陈其峰, 顾云军, 郑君, 李江涛, 李治国, 权伟龙, 付志坚, 李成军 2017 科学通报 62 812

    Chen Q F, Gu Y J, Zheng J, Li J T, Li Z G, Quan W L, Fu Z J, Li C J 2017 Chin. Sci. Bull. 62 812

    [9]

    于继东, 李平, 王文强, 吴强 2014 物理学报 63 116401

    Yu J D, Li P, Wang W Q, Wu Q 2014 Acta Phys. Sin. 63 116401

    [10]

    Danel J F, Kazandjian L, Zérah G 2008 Phys. Plasmas 15 072704

    [11]

    Fu Z J, Quan W L, Zhang W, Li Z G, Zheng J, Gu Y J, Chen Q F 2017 Phys. Plasmas 24 013303

    [12]

    汤文辉, 徐彬彬, 冉宪文, 徐志宏 2017 物理学报 66 030505

    Tang W H, Xu B B, Ran X W, Xu Z H 2017 Acta Phys. Sin. 66 030505

    [13]

    张颖, 陈其峰, 顾云军, 蔡灵仓, 卢铁城 2007 物理学报 56 1318

    Zhang Y, Chen Q F, Gu Y J, Cai L C, Lu T C 2007 Acta Phys. Sin. 56 1318

    [14]

    Chen Q F, Cai L C, Gu Y J, Gu Y 2009 Phys. Rev. E 79 016409

    [15]

    Chen Q F, Zheng J, Gu Y J, Chen Y L, Cai L C 2011 Phys. Plasmas 18 112704

    [16]

    Quan W L, Chen Q F, Fu Z J, Sun X W, Zheng J, Gu Y J 2015 Phys. Rev. E 91 023106

    [17]

    Apfelbaum E M 2015 Phys. Plasmas 22 092703

    [18]

    Apfelbaum E M 2017 High Temp. 55 1

    [19]

    Kuhlbrodt S, Holst B, Redmer R 2005 Contrib. Plasma Phys. 45 73

    [20]

    Moldabekov Z A, Groth S, Dornheim T, Bonitz M, Ramazanov T S 2017 Contrib. Plasma Phys. 57 532

    [21]

    Stolzmann W, Blöcker T 1996 Astron. Astrophys. 314 1024

    [22]

    Fortov V E, Altshuler L V, Trunin R F, Funtikov A I 2004 High-Pressure Shock Compression of Solids VII (New York: Springer) pp437−489

    [23]

    付志坚, 陈其峰, 陈向荣 2011 物理学报 60 055202

    Fu Z J, Chen Q F, Chen X R 2011 Acta Phys. Sin. 60 055202

    [24]

    Redmer R, Rother T, Schmidt K, Kraeft W D, Röpke G 1988 Contrib. Plasma Phys. 28 41

    [25]

    Apfelbaum E M 2011 Phys. Rev. E 84 066403

    [26]

    Karakhtanov V S, Redmer R, Reinholz H, Röpke G 2011 Contrib. Plasma Phys. 51 355

    [27]

    More R M, Warren K H, Young D A, Zimmerman G B 1988 Phys. Fluids 31 3059

    [28]

    Iyetomi H, Ichimaru S 1986 Phys. Rev. A 34 433

    [29]

    Kresse G, Hafner J 1993 Phys. Rev. B 47 558

    [30]

    Zaghloul M R 2018 High Energy Density Phys. 26 8

    [31]

    Wu B, Shin Y C 2006 Appl. Phys. Lett. 89 111902

    [32]

    Morel V, Bultel A, Chéron B G 2009 Int. J. Thermophys. 30 1853

    [33]

    Preston T R, Vinko S M, Ciricosta O, Chung H K, Lee R W, Wark J S 2013 High Energy Density Phys. 9 258

    [34]

    Stransky M 2016 Phys. Plasmas 23 012708

    [35]

    Kemp A J, Meyer-Ter-Vehn J 1998 Nucl. Instrum. Methods Phys. Res., Sect. A 415 674

    [36]

    Krisch I, Kunze H J 1998 Phys. Rev. E 58 6557

    [37]

    田庆云 2015 硕士学位论文 (长沙: 国防科学技术大学)

    Tian Q Y 2015 M. S. Thesis (Changsha: National University of Defense Technology) (in Chinese)

    [38]

    Son S K, Thiele R, Jurek Z, Ziaja B, Santra R 2014 Phys. Rev. X 4 031004

    [39]

    Ebeling W, Föster A, Fortov V E, Gryaznov V K, Polishuk A Y 1991 Thermophysical Properties of Hot Dense Plasmas (Stuttgar-Leipzig: Teubner Verlagsgesellshaft) pp39—42

  • [1] 马桂存, 张其黎, 宋红州, 李琼, 朱希睿, 孟续军. 温稠密物质物态方程的理论研究. 物理学报, 2017, 66(3): 036401. doi: 10.7498/aps.66.036401
    [2] 卢铁城, 陈其峰, 顾云军, 蔡灵仓, 张 颖. 部分电离稠密氦等离子体物态方程的自洽变分计算. 物理学报, 2007, 56(3): 1318-1324. doi: 10.7498/aps.56.1318
    [3] 田杨萌, 王彩霞, 程新路, 杨向东, 姜 明. 惰性物质等离子体物态方程研究. 物理学报, 2007, 56(10): 5698-5703. doi: 10.7498/aps.56.5698
    [4] 侯 永, 袁建民. 第一性原理对金的高压相变和零温物态方程的计算. 物理学报, 2007, 56(6): 3458-3463. doi: 10.7498/aps.56.3458
    [5] 蔡灵仓, 顾云军, 经福谦, 陈志云, 田春玲. 用多次冲击压缩方法研究稠密氢氦等摩尔混合气体的物态方程. 物理学报, 2007, 56(7): 4180-4186. doi: 10.7498/aps.56.4180
    [6] 李晓杰. 热膨胀型固体物态方程. 物理学报, 2002, 51(5): 1098-1102. doi: 10.7498/aps.51.1098
    [7] 赵艳红, 刘海风, 张弓木. 基于统计物理的爆轰产物物态方程研究. 物理学报, 2007, 56(8): 4791-4797. doi: 10.7498/aps.56.4791
    [8] 王彩霞, 田杨萌, 程新路, 杨向东, 姜 明, 孟川民. 一种计算氩等离子物态方程的简单模型. 物理学报, 2006, 55(11): 5784-5789. doi: 10.7498/aps.55.5784
    [9] 顾云军, 郑君, 陈志云, 陈其峰, 蔡灵仓. H2+He流体混合物在部分离解区的物态方程. 物理学报, 2010, 59(7): 4508-4513. doi: 10.7498/aps.59.4508
    [10] 武娜, 杨皎, 肖芬, 蔡灵仓, 田春玲. 固氪物态方程的关联量子化学计算. 物理学报, 2014, 63(14): 146102. doi: 10.7498/aps.63.146102
    [11] 吴强, 谭华, 耿华运. 热力学物态方程参数的统计力学表示. 物理学报, 2001, 50(7): 1334-1339. doi: 10.7498/aps.50.1334
    [12] 彭小娟, 刘福生. 液相及固/液混合相区金属物态方程Grover模型存在的问题及修正. 物理学报, 2012, 61(18): 186201. doi: 10.7498/aps.61.186201
    [13] 郑君, 顾云军, 陈其峰, 陈志云. 稀有气体在电离区压缩特性研究. 物理学报, 2010, 59(10): 7472-7477. doi: 10.7498/aps.59.7472
    [14] 华劲松, 经福谦, 谭 华. 一种获得剪切模量压强二阶偏导数G″P的方法. 物理学报, 2000, 49(12): 2443-2447. doi: 10.7498/aps.49.2443
    [15] 金康, 经光银. 双电层相互作用下主动粒子系统的压强. 物理学报, 2019, 68(17): 170501. doi: 10.7498/aps.68.20190435
    [16] 宋海峰, 刘海风. 金属铍热力学性质的理论研究. 物理学报, 2007, 56(5): 2833-2837. doi: 10.7498/aps.56.2833
    [17] 孟续军, 孙永盛, 王藩侯, 陈敬平, 周显明, 李西军, 经福谦. 冲击压缩产生的氩等离子体辐射不透明度研究. 物理学报, 2001, 50(7): 1308-1312. doi: 10.7498/aps.50.1308
    [18] 朱希睿, 孟续军. 改进的含温有界原子模型对金的电子物态方程的计算. 物理学报, 2011, 60(9): 093103. doi: 10.7498/aps.60.093103
    [19] 范小兵, 陈俊祥, 向士凯. 基于比热的完全物态方程. 物理学报, 2016, 65(23): 236401. doi: 10.7498/aps.65.236401
    [20] 于继东, 李平, 王文强, 吴强. 金属铝固液气完全物态方程研究. 物理学报, 2014, 63(11): 116401. doi: 10.7498/aps.63.116401
  • 引用本文:
    Citation:
计量
  • 文章访问数:  115
  • PDF下载量:  10
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-12-02
  • 修回日期:  2020-02-04
  • 刊出日期:  2020-05-01

温稠密铝等离子体物态方程及其电离平衡研究

  • 1. 河北工业大学, 省部共建电工装备可靠性与智能化国家重点实验室, 天津 300130
  • 2. 河北工业大学, 河北省电磁场与电器可靠性重点实验室, 天津 300130
  • 通信作者: 王坤, kunwang@hebut.edu.cn
    基金项目: 国家自然科学基金(批准号: 51807050)、河北省自然科学基金(批准号: E2019202297)和河北省高等学校青年拔尖人才项目(批准号: BJ2017038)资助的课题

摘要: 温稠密物质的物性参数在惯性约束聚变能源、Z箍缩等高能量密度物理领域的实验结果分析和物理过程数值模拟等方面有着重要的应用价值. 本文应用部分电离等离子体模型, 在理想自由能的基础上考虑了库仑相互作用、排斥体积作用和极化作用等非理想特性, 开展了温稠密等离子体物态方程和电离平衡的研究. 计算了温稠密铝等离子体的压强等物态方程数据和在密度为1.0 × 10–4—3.0 g/cm3, 温度为1.0 × 104—3.0 × 104 K范围内的粒子组分. 计算结果显示, 铝等离子体的平均电离度在临界密度区域内随着密度的增加而突然增大. 根据非理想Saha方程中有效电离能这一关键参数, 分析了铝等离子体平均电离度在临界密度区域内随密度迅速增大的现象.

English Abstract

    • 温稠密物质是指介于固态和理想等离子体态之间的物质形态[1], 常见于丝阵Z箍缩[2]、X射线源[3]、惯性约束聚变能源[4]等高能量密度物理领域. 物态方程和电离平衡模型描述一定热力学条件下物质的宏观性质和粒子组分[5,6], 是高能量密度物理领域相关实验结果分析和物理过程数值模拟的重要物性参数.

      物质在固态和理想等离子体态都有较为精确的物态方程模型. 基于晶格点阵运动理论的Grüneisen物态方程是比较完善的高压固体物态方程[7], 而理想等离子体满足经典Boltzmann统计, 可以采用理想气体物态方程进行描述. 温稠密物质处于部分简并、强耦合和非理想等离子体状态[8], 这给准确描述温稠密等离子体特性带来了极大的困难和挑战. 研究者通常将不同相区的物态方程模型进行衔接处理, 或采用半经验的方法来建立适用温度-密度范围宽的多项物态方程[9]. Thomas-Fermi (TF)模型及相关的改进模型以其适用温度-密度范围宽、数学处理简洁等优点常用于描述温稠密等离子体的性质[10], 然而TF模型只能估算等离子体的平均电离度, 不能得到不同价态的离子分布. 为了获得更为精确的物性参数, 从头算方法、化学模型等也应用于温稠密物质物态方程和输运参数研究[11]. 汤文辉等[12]根据离子耦合系数和电子简并度参数对等离子体进行了划分, 并对不同区域内等离子体物态方程理论进行了总结. 张颖等[13]结合了微扰理论和线性响应理论, 在位形相互作用中考虑了原子间的微扰自由能和量子力学修正自由能, 库仑相互作用中考虑了电子间的交换作用, 能够自洽地对电离能进行修正, 建立了稠密氦等离子体的物态方程, 计算了电离度、压强等物性参数, 得到了粒子间的非理想相互作用引起的电离能降低是出现压致电离现象的主要原因的重要结论, 并将模型进一步改进, 推广至氩、疝等惰性元素等离子体物态方程的计算中[14-16]. 俄罗斯高温联合研究院的Apfelbaum[17,18]在提出的广义化学模型中考虑了电子间交换相互作用的影响, 计算了镍、钛、锌等离子体在温稠密区域的物态方程. 罗斯托克大学的Kuhlbrodt等[19]利用考虑了库仑相互作用和极化作用的自由能模型计算了等离子体组分及输运参数. 考虑非理想特性的温稠密等离子体特性研究显示平均电离度在临界密度区域内随着密度增加突然增大, 而解释这一现象需要针对温稠密等离子体物态方程和电离平衡开展更为深入的研究.

      本文基于部分电离等离子体模型, 在理想自由能模型的基础上考虑了库仑相互作用、排斥体积作用和极化作用等非理想特性, 对温稠密铝等离子体物态方程和电离平衡开展了研究. 获得了宽温度-密度范围内的物态方程和等离子体组分数据, 分析了非理想特性对等离子体组分内不同价态离子电离能的影响. 通过自由能和电离能的变化对临界密度区域内温稠密铝等离子体平均电离度随密度增加突然上升的现象进行了分析.

    • 温稠密等离子体处于部分电离等离子体状态, 由电子、原子和不同价态的离子组成. 对于给定体积V、温度T的热力学系统, 其自由能为F, 本文采用自由能密度f = F/V来建立温稠密等离子体物态方程模型. 等离子体的自由能密度f可以分为理想自由能部分fid和非理想自由能部分fint, 理想自由能密度可以表示为$ f_{\rm id} = f_{\rm id}^{\rm i} + f_{\rm id}^{\rm e} $, 其中上标i, e分别表示重粒子和电子. 重粒子的理想自由能密度可以表示成如下形式[16,18]

      $f_{{\rm{id}}}^{\rm{i}} = {k_{\rm{B}}}T\sum\limits_{k = 0, 1, 2, \cdots } {{n_k}\left[ {\ln \left( {\frac{{{n_k}\lambda _k^3}}{{{Q_k}}}} \right) - 1} \right]} ,$

      其中kB是玻尔兹曼常数; nk是粒子数密度; Qk是粒子的内部电子配分函数; λk = h/(2πmkkBT)1/2是德布罗意波长, h是普朗克常数, mk是粒子质量. 电子的理想自由能密度可以表示为[16]

      $f_{{\rm{id}}}^{\rm{e}}={n_{\rm{e}}}{k_{\rm{B}}}T\left( {{\eta _{{\rm{id}}}} - \frac{{2{I_{{3 / 2}}}\left( {{\eta _{{\rm{id}}}}} \right)}}{{3{I_{{1 / 2}}}\left( {{\eta _{{\rm{id}}}}} \right)}}} \right),$

      其中ne是电子数密度; ηid = $\mu_{{\rm{id}}}^{\rm{e}} $/kBT, $\mu_{{\rm{id}}}^{\rm{e}} $是理想电子化学势; In是Fermi-Dirac积分.

      当等离子体处于温稠密区域时, 粒子间的相互作用不可忽略[20]. 在本文中, 非理想自由能密度fint主要考虑带电粒子之间的库仑相互作用fcoul、重粒子之间的排斥体积作用fhs以及中性原子与带电粒子之间的极化作用fpol的贡献, 如下式所示:

      ${f_{{\rm{int}}}} = {f_{{\rm{coul}}}} + {f_{{\rm{hs}}}} + {f_{{\rm{pol}}}}.$

      带电粒子之间的库仑相互作用可以分为四部分[21]

      ${f_{{\rm{coul}}}} = f_{{\rm{ee}}}^{\rm{c}} + f_{{\rm{ie}}}^{\rm{c}} + f_{{\rm{ii}}}^{\rm{c}} + f_{{\rm{ii}}}^{{\rm{cq}}},$

      其中$f_{{\rm{ee}}}^{\rm{c}} $表示电子间相互作用项, $f_{{\rm{ie}}}^{\rm{c}} $表示电子-离子相互作用项, $f_{{\rm{ii}}}^{\rm{c}} $表示离子间相互作用项, $f_{{\rm{ii}}}^{\rm{cq}} $表示离子间量子修正项. 等离子体物态方程中库仑相互作用项各部分采用Pade近似公式描述[17,22], 其数学表达式较为复杂, 列于附录A中.

      排斥体积作用对自由能密度的贡献可以表示为[19,23]

      ${f_{{\rm{hs}}}} \!=\! n{k_{\rm{B}}}T\Big[ {X\frac{\eta }{{{{\left( {1 - \eta } \right)}^2}}} \!+ \!3Y\frac{\eta }{{1 \!-\! \eta }} \!+\! \left( {X - 1} \right)\ln (1 \!-\! \eta )} \Big],$

      式中n是原子核数密度, $\eta = \displaystyle\sum\nolimits_k {\dfrac{{4{\text{π}}{n_k}R_{_k}^3}}{3}} $, X = $ d_2^3/d_3^2 $, Y = d2d1/d3, ${d_i} = \displaystyle\sum\nolimits_k {{2^i}R_k^i{\alpha _k}} $, Rkαk分别为不同电离度离子的半径和相对粒子分数.

      带电粒子和中性原子间的极化作用可以表示为[24,25]

      ${f_{{\rm{pol}}}} = 2{k_{\rm{B}}}T{n_0}\sum\limits_{j = {\rm e},0,1, \cdots } {{n_j}{B_{0j}}} ,$

      其中B0j是第二维里系数, 可以由极化势UPP计算, $ B_{0j} = \displaystyle\int U_{\rm PP}{\rm d}^3 r $, 极化势表达式为

      ${U_{{\rm{PP}}}} = \frac{{{e^2}{\alpha _{\rm{D}}}\exp \left( { - 2\kappa r} \right)}}{{2{{\left( {4{\text{π}}{\varepsilon _0}} \right)}^2}{{\left( {{r^2} + r_{\rm{c}}^2} \right)}^2}}}{\left( {1 + \kappa r} \right)^2},$

      其中ε0是真空介电常数, αD是偶极化率, κ是德拜长度, rc是截断半径. 本文采用部分电离等离子体模型来描述温稠密等离子体, 因此极化作用项中只计算了电子-中性原子相互作用. 通过上述物态方程模型可以获得温稠密等离子体压强、粒子化学势等热力学量.

    • 电离平衡模型对准确计算物态方程和表征等离子体粒子组分有着重要作用. 铝等离子体内电离过程如下式所示:

      $\begin{split} & {\rm{A}}{{\rm{l}}^k} \leftrightarrow {\rm{A}}{{\rm{l}}^{k + 1}} + {\rm e}, \\ & \quad k = 0,1,2, \cdots , \end{split} $

      式中k = 0表示中性原子. 若等离子体处于电离平衡状态, 不同电离态离子、中性原子和电子也处于相平衡状态, 由化学势表示的多粒子体系电离平衡对应的相平衡表达式为[26]

      ${\mu ^k} = {\mu ^{k + 1}} + {\mu ^{\rm{e}}} + {E_{k + 1}},$

      其中Ek + 1是电离能, μk是重粒子化学势, μe是电子化学势.

      温稠密等离子体物态方程中自由能密度划分为理想部分和非理想部分. 与之相对应, 等离子体中粒子化学势也可划分为理想部分和非理想部分, 即μ = μid + μint, 将该式代入(9)式可以得到

      $\begin{split} &\mu _{{\rm{int}}}^k - \mu _{{\rm{int}}}^{k + 1} - \mu _{{\rm{int}}}^{\rm{e}} \\ ={}& \mu _{{\rm{id}}}^{k + 1} - \mu _{{\rm{id}}}^k + \mu _{{\rm{id}}}^{\rm{e}} + {E_{k + 1}},\end{split} $

      等离子体中粒子化学势的理想部分及非理想部分可由自由能密度求解, 如下式所示:

      $\mu _{{\rm{id}}/{\rm{int}}}^k = \frac{{\partial {f_{{\rm{id}}/{\rm{int}}}}}}{{\partial {n_k}}}.$

      将(1)式和(11)式代入(10)式, 化简可得到非理想Saha方程

      $\frac{{{n_{k + 1}}}}{{{n_k}}} = \frac{{{Q_{k + 1}}}}{{{Q_k}}}{{\rm{e}}^{ - \frac{{\mu _{{\rm{id}}}^{\rm{e}} + E_{{\rm{eff}}}^{k + 1}}}{{{k_{\rm{B}}}T}}}},$

      其中$E_{{\rm{eff}}}^{k + 1}$ = Ek + 1 – ∆E是有效电离能, ∆E = $\mu_{{\rm{int}}}^{k} $$\mu_{{\rm{int}}}^{k + 1} $$\mu_{{\rm{int}}}^{\rm{e}} $是非理想效应导致的电离能下降.

      非理想Saha方程结合质量守恒定律和电荷守恒定律形成的封闭方程组, 通过与物态方程模型迭代求解可以得到等离子体的组分. 质量守恒定律和电荷守恒定律如(13)式和(14)式所示:

      $\sum\nolimits_{k = 0}^Z {{n_k} = n} ,$

      $\sum\nolimits_{k = 0}^Z {k{n_k} = {n_{\rm{e}}}} .$

      本文应用上述物态方程和电离平衡模型针对温稠密铝等离子体物性参数开展了研究, 模型中相关参数(里德伯单位制)详见表1.

      参数参数参数
      R02.5714R60.4167E32.0909
      R12.3377rc1.7712E48.8192
      R22.1429mk24597E511.3059
      R30.4678αD46.281E614.0008
      R40.4494E10.4400
      R50.4324E21.3839

      表 1  物态方程及电离平衡模型参数列表

      Table 1.  Parameters of equation of state and ionization equilibrium model.

    • 温稠密等离子体压强等热力学量可通过自由能求得P = –($\partial $F/$\partial $V)|T. 本文模型计算的铝等离子体压强在0.1 g/cm3时随温度的变化曲线见图1, 由SCAALP-FT, QEOS以及VASP等模型[27-29]计算的结果也列于图中对比. 本文得到的压强曲线与VASP, SCAALP-FT模型计算的结果较为接近, 而QEOS模型的结果明显高于其他模型, 并且差别随着温度升高而变大.

      图  1  在密度为0.1 g/cm3时, 不同模型计算的铝等离子体压强随温度的变化

      Figure 1.  The pressure of aluminum plasma calculated by different models as a function of temperature at density of 0.1 g/cm3.

      描述多粒子体系的物态方程与非理想Saha方程迭代求解可获得等离子体的粒子组分. 本文计算了铝等离子体在密度1.0 × 10–4 g/cm3 < ρ < 3.0 g/cm3, 温度1.0 × 104 K < T < 3.0 × 104 K范围内平均电离度的等值线分布, 结果如图2所示. 等离子体的电离机制在低密度、高温度区域由热电离主导, 而在高密度、低温度区域由压致电离主导. 临界点是气-液相界线的终点, 已有大量研究工作对铝的临界点参数进行了估算[30-32], 一般认为金属铝的临界温度处于8.0 × 103—1.0 × 104 K范围内, 临界密度处于0.4—1.0 g/cm3范围内. 在本文中, 将上述密度区域称为临界密度区域. 在临界密度区域内, 平均电离度呈现随密度增大而迅速上升趋势. 在高密度区域, 当温度小于1.2 × 104 K时, 铝等离子体主要发生一级电离和二级电离, 体系主要由Al1+, Al2+和电子组成, 平均电离度的最大值为2; 而当温度大于1.2 × 104 K时, 二级电离随密度增大而减弱, 铝等离子体发生三级电离, 等离子体中带电粒子主要由Al3+和电子组成, 平均电离度的最大值为3.

      图  2  铝等离子体的平均电离度随密度与温度变化等值线分布

      Figure 2.  Contour map of average ionization degree of aluminum plasma as a function of density and temperature.

      等离子体中不同粒子的相对粒子分数在密度为0.1 g/cm3时随温度的变化如图3所示. 从图3可以看出, 铝等离子体在温度为1.0 × 104—2.0 × 106 K的范围内电离形成的离子最高, 为六价铝离子. 铝原子相对粒子分数随着温度的升高不断降低, 当温度达到3.0 × 104 K时, 铝原子相对粒子分数已低于0.1. 等离子体中一价到五价离子的相对粒子分数随着温度的增加都呈现出先增大后减小的变化规律. 各价态离子的相对粒子分数峰值所对应的温度随着电离态的增大而逐渐升高. 在温度为1.0 × 104—3.0 × 104 K、密度为1.0 × 10–4—3.0 g/cm3的范围内, 等离子体中重粒子组分主要由Al0, Al1+, Al2+和Al3+组成, 更高电离态离子可忽略不计.

      图  3  等离子体中不同粒子的相对粒子分数在密度为0.1 g/cm3时随温度的变化

      Figure 3.  Dependence of relative particle fraction of different particles on temperature at density of 0.1 g/cm3.

      采用非理想Saha方程计算得到的铝等离子体在温度为10000和15000 K时平均电离度随密度的变化如图4所示. 德拜理论修正的Saha方程、Ecker-Kröll (DmEK)模型[33,34]、FEOS模型[35]和COMPTRA04程序[19]等计算的结果以及实验测量的平均电离度数据[36]也列于图4中进行对比. 从图4(a)可以看出, 本文计算结果、德拜理论修正的Saha方程、DmEK模型和COMPTRA04程序计算的平均电离度分别在密度为0.63, 0.02, 1.71和0.21 g/cm3时出现突然上升的现象. 德拜理论修正的Saha方程在耦合系数小于1的范围内计算结果处于比较合理的范围内, 但在耦合系数大于1的高密度区域, 其计算结果存在比较大的偏差. DmEK模型在小于临界密度区域内能够描述电离机制由热电离向压致电离的转变, 但平均电离度突变时的密度较高. COMPTRA04程序计算的平均电离度在高密度区域最高达到2, 本文计算结果显示平均电离度在10000 K时最高可达到2, 而在15000 K时可达到3. 平均电离度随密度的突然增加导致采用过小的体积步长计算压强等热力学量时出现波动, 因此在该温度、密度区域内采用稍大的体积步长来获得平滑的压强曲线.

      图  4  不同模型计算的铝等离子体平均电离度在不同温度下随密度的变化 (a) 10000 K; (b) 15000 K

      Figure 4.  Average ionization degree of aluminum plasma calculated by different models as a function of density at different temperatures: (a) 10000 K; (b) 15000 K.

      物态方程模型中非理想特性对非理想Saha方程准确描述等离子体电离机制随密度的增加由热电离到压致电离的转变有非常重要的影响. 图5是不同温度下铝等离子体物态方程中非理想部分的库仑相互作用、排斥体积作用、极化作用和等离子体中电子以及原子的相对粒子分数随密度的变化. 随着密度的增加, 等离子体由弱耦合区域过渡到强耦合区域, 非理想效应增强.

      图  5  不同温度下非理想效应自由能密度以及电子、原子相对粒子分数随密度的变化 (a) 库仑相互作用; (b) 排斥体积作用; (c) 极化作用; (d) 电子、原子相对粒子分数

      Figure 5.  Free energy density of different non-ideal effects and relative particle fraction for electrons and atoms as a function of density at different temperatures: (a) Coulomb interaction; (b) excluded volume effect; (c) polarization effect; (d) relative particle fraction for electrons and atoms.

      图5可以看出, 排斥体积作用和极化作用等非理想部分远小于库仑相互作用. 因此, 温稠密铝等离子体的非理想特性由库仑相互作用主导, 这与Apfelbaum[18]在其文章中所述结论一致. 由(11)式可知, 非理想自由能密度的急剧变化会导致粒子化学势的变化, 进而影响有效电离能和等离子体组分. 各非理想特性随密度、温度的变化趋势并不一致. 库仑相互作用随着密度的增加而增强, 但随温度的变化并不明显. 如图5(b)所示, 排斥体积作用在密度小于1 g/cm3的区域内随密度的增加呈现先增强后减弱的变化规律. 当温度为15000 K时, 排斥体积作用在0.54 g/cm3时达到最大值. 此时, 铝原子相对粒子分数已很低, 等离子体平均电离度接近1, 处于平均电离度迅速增大的密度区域. 极化作用在密度为0.42 g/cm3时最强, 随着密度进一步增加, 铝原子相对粒子分数降低导致极化作用减弱. 排斥体积作用和极化作用的峰值所对应的密度随着温度的升高稍有降低. 图5(d)显示在温度15000 K时铝原子相对粒子分数在0.03 g/cm3时达到最大值0.27. 随着温度升高, 铝原子相对粒子分数不断降低. 当温度为30000 K时, 铝原子相对粒子分数随密度变化的峰值低于0.1. 在平均电离度随密度迅速增大的高密度区域, 铝原子几乎全部被电离.

      稠密等离子体处于部分电离等离子体状态, 很高的重粒子数密度使得电离产生的自由电子可占据的自由空间减少, 而带电粒子也会使原子和离子自身的化学势发生改变, 重粒子之间的空间也因此具有较低的势能, 这会导致离子的能级结构发生变化, 出现电离能下降的现象[34,37]. 为了获得电离能的分布, 本文计算了15000 K时各粒子相互作用部分化学势随密度的分布, 如图6所示. 从图6可以看出, 粒子的非理想化学势随密度增加而增大. 当密度达到0.5 g/cm3时, 由于离子密度增加, 非理想效应增强, 价态较高的离子非理想化学势的数值增大趋势明显.

      图  6  温度为15000 K时等离子体中不同粒子化学势的非理想部分随密度的变化

      Figure 6.  Dependence of non-ideal chemical potential of particles on density at temperature of 15000 K.

      粒子的非理想部分化学势决定了电离能下降, 本文模型计算的铝等离子体在温度为15000 K时不同价态离子电离能降低随密度的变化如图7所示. DmEK模型计算结果也画于图中对比. 从图7可以看出, 不同价态离子电离能的降低均随密度增加而增大, 且价态越高的离子电离能下降程度越大, 这与Son等[38]在其文章中所述结论是一致的. 非理想效应的增强导致有效电离能减小, Al1+, Al2+和Al3+的有效电离能分别在密度为0.11, 0.53和0.54 g/cm3时接近0, 过小的电离能阈值会使原子和低价态离子进一步电离. 当密度高于0.54 g/cm3时, Al2+全部电离为Al3+, 而这导致铝等离子体的平均电离度在临界密度区域出现突然上升的趋势. DmEK模型计算的电离能在密度小于0.1 g/cm3的区域与本文计算结果符合较好, 但低估了各级粒子电离能的降低, 从而平均电离度在密度约为3 g/cm3时才出现增大的现象.

      图  7  温度为15000 K时不同模型计算的电离能的降低∆E随密度的变化(黑色曲线, 非理想Saha方程; 红色曲线, DmEK模型)

      Figure 7.  Depression of ionization potential calculated by different models as a function of density at 15000 K. Black lines correspond to nonideal Saha equation; red lines correspond to DmEK model.

    • 温稠密等离子体态是高能量密度物理领域中最为常见的物理状态, 准确的描述温稠密态物质的特性是非常复杂的物理问题. 本文在理想自由能的基础上考虑了库仑相互作用、排斥体积作用和极化作用等重要的非理想特性, 结合非理想Saha方程开展了温稠密铝等离子体物态方程和电离平衡研究, 获得了压强、等离子体粒子组分等数据, 并与其他模型计算结果及相关实验数据作了对比、分析. 不同温度、密度下的非理想自由能曲线表明, 库仑相互作用随着密度增大而增强, 但对温度并不敏感; 在小于1 g/cm3的密度范围内, 排斥体积作用和极化作用都随着密度增加呈现先增强后减弱的变化规律. 多粒子体系物态方程与非理想Saha方程迭代求解获得了等离子体中不同电离态离子分布. 在温度为10000 K时, 高密度等离子体平均电离度为2, 而在温度为15000 K时, 高密度等离子体平均电离度为3. 通过等离子体中粒子非理想化学势分析了由于非理想特性导致的不同价态离子的电离能下降. 计算结果显示, 在某些温度、密度下, 离子的电离能下降与其电离能相当, 从而导致有效电离能接近0. 在1.0 × 104 K < T < 3.0 × 104 K的温度范围内, 当等离子体密度达到0.54 g/cm3时, Al3+的有效电离能降低至接近0, 从而导致Al2+电离成为Al3+, 等离子体平均电离度在这一密度区域内随着密度增加迅速地增大. 本文模型能够计算密度为1.0 × 10–4—3.0 g/cm3范围内的铝等离子体物态方程和粒子组分, 当密度进一步增加时, 迭代求解计算量增加, 迭代收敛性逐渐变差, 电离平衡方程变得难以求解, 并且过小的体积步长也会在临界密度区域内对压强的计算结果造成影响. 今后将对上述问题开展进一步深入的研究.

    • 本文采用Pade近似公式描述库仑相互作用, 电子-电子相互作用$f_{{\rm{ee}}}^{\rm{c}} $、电子-离子相互作用$f_{{\rm{ie}}}^{\rm{c}} $和离子-离子相互作用$f_{{\rm{ii}}}^{\rm{c}} $的表达式如下[39]:

      $f_{{\rm{ee}}}^{\rm{c}} = - {n_{\rm{e}}}\frac{{{f_0}{{\bar n}^{{1 / 2}}} + {f_3}\bar n + {f_2}{{\bar n}^2}{\varepsilon _{\rm{c}}}\left( {{r_{\rm{s}}}} \right)}}{{1 + {f_1}{{\bar n}^{{1 / 2}}} + {f_2}{{\bar n}^2}}},\tag{A1}$

      $\begin{split} \; & f_{{\rm{ie}}}^{\rm{c}} = - {n_ + }{k_{\rm{B}}}T\\ & \times\frac{{{Q_0}{{\tilde n}^{{1 / 2}}} + {Q_2}{{\bar n}^{ {3 / 2}}}{e_2}\left( {\hat n,{r_{\rm{s}}},{k_{\rm{B}}}T} \right)}}{{1 \!+\! {Q_1}{{\left( {\tilde n{k_{\rm{B}}}T} \right)}^{{1 / 2}}} \!+\! {Q_2}{{\bar n}^{ {3 / 2}}} + {Q_4}{{\tilde n}^{{1 / 2}}}\ln \left( {1 + \dfrac{{{Q_5}}}{{{{\tilde n}^{{1 / 2}}}}}} \right)}},\end{split}\tag{A2}$

      $f_{{\rm{ii}}}^{\rm{c}} = - {n_ + }{k_{\rm{B}}}T\frac{{{q_0}{{\tilde n}^{{1 / 2}}} + {q_2}{{\hat n}^{{3 / 2}}}{e_1}\left( {\hat n} \right)}}{{1 + {q_1}{{\left( {\tilde n{k_{\rm{B}}}T} \right)}^{{1 / 2}}} + {q_2}{{\hat n}^{{3 / 2}}}}},\tag{A3}$

      其中f0—2, e1—2, q0—2以及Q0—5的表达式为

      ${f_0} = \frac{2}{3}{\left( {\frac{{{k_{\rm{B}}}T}}{{\text{π}}}} \right)^{{1 / 4}}},\tag{A4}$

      ${f_1} = \left[ {\frac{{\sqrt 2 }}{8} + \frac{1}{4}{{\left( {\frac{{{k_{\rm{B}}}T}}{{\text{π}}}} \right)}^{{1 / 2}}} - \frac{{{k_{\rm{B}}}T}}{{8\sqrt 2 }}} \right]\Big/{f_0},\tag{A5}$

      ${f_2} = 3,~~{f_3} = \frac{1}{4}{\left( {\frac{{{k_{\rm{B}}}T}}{{\text{π}}}} \right)^{{1 / 2}}},\tag{A6}$

      $\begin{split}{e_1}\left( x \right) =\; & 1.4474{x^{{1 / 3}}} - 4.2944{x^{{1 / {12}}}} + 0.6712{x^{{{ - 1} / {12}}}} \\ &+ 0.2726\ln x + 2.983,\end{split}\tag{A7}$

      $\begin{split}{e_2}\left( {x,y,z} \right) = \; & \frac{{0.8511y{x^{{1 / {12}}}}}}{{\left( {1 + 0.3135{x^{{{ - 1} / {12}}}}} \right)\left( {1 + 1.137y{z^{{1 / 2}}}} \right)}} \\ & + \frac{{y{x^{{1 / 3}}}\left( {0.0726 + 0.0161y} \right)}}{{1 + 0.0887{y^2}}},\end{split}\tag{A8}$

      ${q_0} = \frac{{2{{\text{π}}^{{1 / 2}}}}}{3}{\left\langle {{z^2}} \right\rangle ^{{3 / 2}}},\tag{A9}$

      ${q_1} = \left[ {\frac{{{{\text{π}}^{{3 / 2}}}}}{8}\sum\limits_i {\sum\limits_j {{\zeta _i}{\zeta _j}Z_i^2Z_j^2{{\left( {{\gamma _i} + {\gamma _j}} \right)}^{{1 / 2}}}} } } \right]\Big/{q_0},~~{q_2} = 1000,\tag{A10}$

      ${Q_0} = \frac{{2{{\text{π}}^{{1 / 2}}}}}{3}\left[ {{{\left( {\zeta + \left\langle {{z^2}} \right\rangle } \right)}^{{3 / 2}}} - {\zeta ^{{3 / 2}}} - {{\left\langle {{z^2}} \right\rangle }^{{3 / 2}}}} \right],\tag{A11}$

      ${Q_1} = {{\left( {\frac{{{{\text{π}}^{{3 / 2}}}}}{2}\zeta \left\langle {{z^2}} \right\rangle - {q_1}{q_0}} \right)} \Big/ {{Q_0}}},~{Q_2} = 1,\tag{A12}$

      $\begin{split}{Q_3} ={}& \frac{1}{Q_0} \Big\{0.5236 ( \langle z^3 \rangle - \zeta )\big[ 4\langle z^3\ln z \rangle \\ & + (\langle z^3 \rangle - \zeta ) \ln (29.09 \big(\langle z^2 \rangle + \zeta )\big)\big]\Big\}, \end{split} \tag{A13} $

      ${Q_4} = {{1.0472{{\left( {\left\langle {{z^3}} \right\rangle - \zeta } \right)}^2}}}/{{{Q_0}}},\tag{A14}$

      ${Q_5} = \exp \left( { - {{{Q_3}}}/{{{Q_4}}}} \right),\tag{A15}$

      其中$\left\langle {{z^p}} \right\rangle \!=\! n_ + ^{ - 1}\displaystyle\sum\limits_{i = 1}^Z {Z_i^p} {n_i}$, ${n_ + } \!=\! \displaystyle\sum\limits_i\! {{n_i}} $, $\overline n \!=\! {n_{\rm{e}}}\lambda _{\rm{e}}^{\rm{3}}$, $\tilde n \!=\!\! \dfrac{8{n_ + }} {(k_{\rm B}T)^3}$, $\widehat n = \widetilde n\left\langle z \right\rangle {\left\langle {{z^{5/3}}} \right\rangle ^3}$, $\zeta ={{{n_{\rm{e}}}} / {{n_ + }}}$, ${\zeta _i}={{{n_i}} / {{n_ + }}}$, ${\gamma _i} = {{m_{\rm{e}}}}/{{m_i}} $; εc为电子基态能量, rs = (3/4πnea03)1/3是电子间的平均距离, a0是玻尔半径.

      离子-离子间量子修正项fiicq表达式如下[21]:

      $f_{{\rm{ii}}}^{{\rm{cq}}} = 3{q_0}\left( t \right) + 3{z_4}{q_4}\left( t \right) - \frac{9}{8}\frac{1}{{\sqrt {{r_{\rm{s}}}} }}{z_5}{q_5}\left( t \right) + 3{z_6}{q_6}\left( t \right),\tag{A16}$

      其中z4 = 0.827702, z5 = 1.131, z6 = 0.55045, 函数$ q_i(t) $可以通过(A17)—(A20)式计算

      ${q_0}\left( t \right) = \ln \left( {{{\sinh t}}/{t}} \right),\tag{A17}$

      ${q_4}\left( t \right) = \frac{1}{4} - \frac{t}{8}\coth t - \frac{{{t^2}}}{8}\frac{1}{{{{\sinh }^2}t}},\tag{A18}$

      ${q_5}\left( t \right) = t\left( {\frac{1}{3}{{\coth }^2}t - \frac{3}{{2t}}\coth t + \frac{2}{{{t^2}}} - \frac{5}{6}\frac{1}{{{{\sinh }^2}t}}} \right),\tag{A19}$

      ${q_6}\left( t \right) = - \frac{1}{6} + \frac{t}{{16}}\coth t + \frac{{{t^2}}}{{16}}\frac{1}{{{{\sinh }^2}t}} + \frac{{{t^3}}}{{24}}\frac{{\coth t}}{{{{\sinh }^2}t}}.\tag{A20}$

参考文献 (39)

目录

    /

    返回文章
    返回