搜索

x

留言板

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

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

耗散响应理论及其在开放系统中的应用

陈宇

张国帅, 尹超, 王兆繁, 陈泽, 毛世峰, 叶民友. 中子辐照诱导钨再结晶的模拟研究. 物理学报, 2023, 72(16): 162801. doi: 10.7498/aps.72.20230531
引用本文: 张国帅, 尹超, 王兆繁, 陈泽, 毛世峰, 叶民友. 中子辐照诱导钨再结晶的模拟研究. 物理学报, 2023, 72(16): 162801. doi: 10.7498/aps.72.20230531
Zhang Guo-Shuai, Yin Chao, Wang Zhao-Fan, Chen Ze, Mao Shi-Feng, Ye Min-You. Simulation of neutron irradiation-induced recrystallization of tungsten. Acta Phys. Sin., 2023, 72(16): 162801. doi: 10.7498/aps.72.20230531
Citation: Zhang Guo-Shuai, Yin Chao, Wang Zhao-Fan, Chen Ze, Mao Shi-Feng, Ye Min-You. Simulation of neutron irradiation-induced recrystallization of tungsten. Acta Phys. Sin., 2023, 72(16): 162801. doi: 10.7498/aps.72.20230531

耗散响应理论及其在开放系统中的应用

陈宇

Dissipative linear response theory and its appications in open quantum systems

Chen Yu
Article Text (iFLYTEK Translation)
PDF
HTML
导出引用
  • 近些年来, 随着实验技术的进步, 对量子多体系统的耗散控制能力得到了增强, 同时耗散动力学过程表征技术方面的实验也有了较大进展. 实验上的进展驱使我们在理论上建立量子多体系统的耗散动力学计算体系. 最近我们发现, 通过把系统和环境之间的相互作用看成对系统的一个微扰, 可以得到一般性的耗散响应理论. 通过这一响应理论, 可以回答物理可观测量以及熵在耗散下一定时间尺度内的动力学演化的问题. 本文建立了非Markov 环境下的一般理论, 并讨论了何时可以取到 Markov 近似, 同时综述了这种方法在计算强关联体系的耗散动力学、强相互作用开放体系的熵的动力学演化等方面的应用.
    With the recent development of experimental technology, the ability to control the dissipation in quantum many-body system is greatly enhanced. Meanwhile, many new breakthroughs are achieved in detecting the quantum states and others. All these advances make it necessary to establish a new theory for calculating the dissipative dynamics in strongly correlated sstems. Very recently, we found that by taking the interactions between the system and the bath as a perturbation, a systematic dissipative response theory can be established. In this new approach, the calculation of dissipative dynamics for any physical observables and the entropies can be converted into the calculation of certain correlation functions in initial states. Then we discuss how Markovian approximation at low temperature limit and at high temperature limit can be reached Also, we review the progress of the dissipative dynamics in open Bose-Hubbard model. In the fourth section, we review recent progress of entropy dynamics of quench dynamics of an open quantum system. Finally, we draw a conclusion and discuss possible development in the future.
      通信作者: 陈宇, ychen@gscaep.ac.cn
    • 基金项目: 北京市自然科学基金重点专项(批准号: Z180013)和国家自然科学基金重点项目(批准号: 11734010)资助的课题
      Corresponding author: Chen Yu, ychen@gscaep.ac.cn
    • Funds: Project supported by the Key Program of the Natural Science Foundation of Beijing, China (Grant No. Z180013) and the Key Program of the National Natural Science Foundation of China (Grant No. 11734010).

    钨由于具有高熔点、高热导、高溅射阈值、良好的抗热冲击性能及较低的氘氚滞留等优点, 被确定为国际热核聚变实验堆(ITER)的偏滤器靶板材料, 并且是聚变示范堆(DEMO)面向等离子体部件(plasma-facing component, PFC)中铠甲和靶板的候选材料[1-4]. 然而, 聚变堆严苛的服役环境, 会造成钨作为PFC服役性能的退化, 进而减短部件寿命[5]. 钨在聚变装置中的工作温度涵盖了500—1200 ℃的范围[6]. 研究表明塑性变形加工后的钨在高温下会发生退火现象, 包括回复、再结晶和晶粒生长过程[7-10]. 当钨发生再结晶时, 其转脆温度将从200 ℃上升至500 ℃, 从而在较高的工作温度下可能发生脆性断裂[11]. 因此, 分析并预测钨的再结晶过程对钨作为PFC的应用和评估具有重要意义. 根据未辐照等温退火实验外推预测, 钨在1000 ℃下服役两年后其再结晶分数约为50%[8,12]. 然而, 最近高通量同位素反应(high flux isotope reactor, HFIR)堆裂变中子辐照钨的实验显示, 钨在850 ℃下辐照24 d后的再结晶分数就已达到50%[13], 这表明中子辐照会诱导钨发生再结晶, 增加脆断风险.

    中子辐照会对钨造成辐照损伤, 由于缺陷之间的相互反应, 这些初级损伤产生的缺陷会进一步演化成尺寸更大的缺陷, 包括位错环和空洞[14-17]. 这些非热平衡的缺陷会增大钨内部储能[18], 从而提高了再结晶驱动力. Mannheim等[19,20]通过耦合团簇动力学(cluster dynamic, CD)模型和包含辐照增强再结晶驱动力效应的平均场再结晶模型, 从物理的角度模拟了再结晶过程中晶粒的成核和生长. 该模型预测在1000 ℃辐照温度下辐照近700 h后, 钨几乎不发生再结晶现象, 再结晶速率远低于HFIR堆实验结果, 这可能是该模型未考虑辐照增强晶界迁移效应所致.

    基于此, 本文在Johnson-Mehl-Avrami-Kolmogorov (JMAK)再结晶动力学模型的基础上, 同时考虑辐照增强再结晶驱动力效应和辐照增强晶界迁移效应, 建立了辐照诱导再结晶(irradiation-induced recrystallization, IIR)动力学模型.

    利用CD模型模拟钨在HFIR堆中子辐照下不同尺寸空位(vacancy, V)和自间隙原子(self-interstitial atom, I)团簇随时间的演化过程, 同时将模拟得到的缺陷尺寸与密度输入IIR模型以模拟高温中子辐照过程中钨再结晶分数的变化. 模拟的辐照温度区间为750—1300 ℃, 使用Matlab刚性求解器ODE15s进行数值求解, 相对容差设置为1×10–3, 并且设置每个解为非负值.

    如(1)式所示, 平均场CD模型基于扩散捕获速率理论, 通过一组常微分控制方程(ordinary differential equations, ODEs)来描述不同尺寸缺陷的密度随时间的演化[21]:

    dCθdt=Gθ+θJ+θ,θθJθ,θLθ,
    (1)

    其中Cθ为缺陷密度, 下标θ表示n个V组成的团簇(Vn)或n个I组成的团簇(In); θ为可移动缺陷, t为时间, Gθ为源项, J+θ,θ为吸收反应通量, 用于描述缺陷团簇吸收可移动缺陷的反应; Jθ,θ为发射反应通量, 用于描述缺陷团簇发射可移动缺陷的反应; Lθ为汇项, 用于描述缺陷阱对可移动缺陷的吸收. 除源项外, 吸收和发射反应项、汇项都是直接使用Mannheim等[19]建立和整理的反应式, 相应的参数在见补充材料中给出.

    为了获得源项, 首先通过SPECTRA-PKA程序计算特定中子能谱辐照钨所产生的初级反冲原子(primary knock-on atom, PKA)能谱[22,23], 再将PKA能谱叠加到PKA级联模拟结果[24]以确定级联存活损伤速率(Sdpa), 最后将Sdpa转换成单位体积内的级联存活损伤速率(Gtot):

    Gtot=SdpaρWNAMW,
    (2)

    其中ρW为钨的密度, NA为阿伏伽德罗常数, MW为钨的摩尔质量, 具体数值见补充材料. 另外, SPECTRA-PKA也给出了按照Norgett-Robinson-Torrens (NRT)模型得到的平均每离位原子速率(NRTdpa).

    表1为基于HFIR中子能谱[25]得到的NRTdpa, SdpaGtot. 这与文献[16]中根据NRT模型计算得出的平均每离位原子速率2.3×10–7 dpa/s基本吻合. 同时计算得到的级联存活效率(SdpaNRTdpa的比值)为0.30, 这与文献[26, 27]中MD模拟级联效率的结果一致. 由于级联内缺陷的相互反应, 使得级联内缺陷密度与尺寸呈幂律分布. 按照Vrielink等[28]的方法, 级联尺寸分布由下式给出:

    表 1  CD模型模拟HFIR堆中子辐照钨的源项相关参数
    Table 1.  Parameters related to source term of neutron irradiated tungsten in HFIR reactor simulated by CD model.
    参数NRTdpa/(10–7 dpa·s–1)Sdpa/(10–8 dpa·s–1)Gtot/(1021 m–3·s–1)
    数值2.166.414.06
    下载: 导出CSV 
    | 显示表格
    Gθn={(1fcθ)Gtot,n = 1,Aθ(1n1Nc)Sθ, 2nNc,0, otherwise,
    (3)

    其中fcθSθ分别为温度相关的缺陷团聚分数和幂律指数[19]. 钨低温自离子辐照实验结果表明级联产生的最大In尺寸约为1000[29], 并且MD级联模拟结果表明级联产生的最大Vn尺寸约为最大In尺寸的2倍[30]. 因此本文级联产生的最大Vn尺寸(Nc, V)取2000, 最大In尺寸(Nc, I)取1000.

    本文由于考虑的最大缺陷尺寸大于105, 为了提高计算效率, 在考虑大尺寸缺陷时, 采用Ghoniem和Sharafat[31]提出的Fokker-Plank近似方法. 此外, 测试的算例结果表明大尺寸缺陷(n > 100)的演化忽略I2反应项对缺陷的平均尺寸与密度影响不大, 这与Li等[32]的CD模拟结果一致. 因此本文中, 当缺陷尺寸大于100时, 通过二阶泰勒展开将(1)式转化成尺寸空间连续的形式[31]:

    C[x]t=G[x]x[FθC[x]x(DθC[x])]L[x],
    (4)

    其中, C[x]表示取x整数部分对应尺寸的缺陷密度, Fθ表示缺陷的净通量, Dθ表示缺陷的平均扩散通量, FθDθ分别可以写成[31]:

    Fθ=kθλθ,   Dθ=12(kθ+λθ),
    (5)

    其中kθλθ分别为相应缺陷团簇的生长速率和衰减速率[31].

    在尺寸空间x内按照Li等[32]的方法, 采取非均匀化的网格, 表达式如下:

    x1=1,    xi=xi1+Δxi,  2iNm,Δxi={1,2i100,Δxi1exp(ε),100iNm.
    (6)

    其中ε是远小于1的常数, Nm为考虑的最大缺陷尺寸, 且Nm需要满足边界条件CNm0或小于本研究考虑的缺陷的最大尺寸. 通过Fokker-Plank近似, 本文建立的CD模型可以考虑更大尺寸缺陷的演化, 其中最大的Vn尺寸约为3×105, 对应的空洞直径约20 nm, 最大的In尺寸约为6 × 103, 对应的位错环直径约20 nm.

    为了简化分析, 做出如下假设.

    1)认为整个空间缺陷密度分布是均匀的, 且缺陷密度足够低, 不影响基体钨的性质.

    2)认为只有V, I和I2是可移动的, 其他固定缺陷通过捕获可移动缺陷演化.

    3)所有Vn都是球体: 对于In, 小于7的是球体; 大于7的是位错环.

    4)无论是吸收反应还是反射反应, 只考虑吸收或发射一个可移动性缺陷, 即反应逐次进行.

    5)汇项只考虑位错对缺陷的吸收, 且位错不会因吸收缺陷而消失或改变密度.

    IIR模型通过未辐照下JMAK模型[9,33-35], 建立了中子辐照诱导再结晶模型, 假设中子辐照将会影响钨未发生再结晶的初始晶粒内再结晶驱动力P和晶界迁移率M. 温轧中等变形量钨的等温退火实验表明其再结晶过程符合位点饱和情况下的再结晶[7], 同时中子辐照产生的缺陷会提高钨自体扩散系数及内部储能, 应会加速新晶粒的成核过程, 因此本研究假设辐照下钨再结晶过程均为位点饱和情况. 对于位点饱和成核的JMAK模型, 再结晶分数(X)随时间的变化可以写成下列形式[33]:

    X=1exp[U(PMt)na],
    (7)

    其中U为成核数相关的材料常数; na为Avrami指数, 反映了晶粒的生长维度; 而中子辐照时的P可以写为

    P=EVold+EVolV+EVolI,
    (8)

    其中EVold表示位错造成的储能[33]:

    EVold=αρGb2,
    (9)

    其中α为常数, 取0.5; ρ表示位错线密度, G表示剪切模量, b是位错的伯格斯矢量; 具体数值见补充材料. 由于本研究假定初始晶粒的位错密度(7×1013 m–2[13])在再结晶过程中不变, 因此EVold为常数; EVolV表示空洞造成的储能; EVolI表示I位错环造成的储能. 钨的自离子辐照实验表明1/2111型I位错环占主导[36,37], 且1/2111型I位错环占比还会随温度的升高而增大[38]. 因此本研究中, 假设I位错环全部都是1/2111型. 根据Ma等[18]的模拟结果, 尺寸为n的空洞形成能(Ev(n))和尺寸为n的1/2111; 型I位错环形成能(El-1/2111(n))可以写成下列形式:

    Ev(n)=4πr2ˉγ=4π(3nΩ4π)2/3ˉγ,
    (10)
    El-1/2111(n)=2aE111nlnn+ε1n,
    (11)

    其中, a为晶格常数, r为空洞半径, n表示缺陷团簇包含的单体数, ˉγ表示空洞平均表面能密度, E111表示位错环平均线能量前对数因子, ε1取9.08 eV/[18](1 Å = 10–10 m), 具体数值见补充材料. 再将CD模型得到的缺陷密度信息结合(10)式和(11)式及小缺陷的形成能, 就可以求得任意时刻空位缺陷、自间隙缺陷的体积储能:

    EVolV=CVEfV+CV2EfV2+Nvmn=3CVnEv(n)TSv,
    (12)
    EVolI=6n=1CInEfIn+Nim7CInEl1/2111(n)TSi,
    (13)

    其中T为温度, SvSi分别是V缺陷和I缺陷的构型熵, 可利用Mannheim等[19]的方法得到.

    由于晶界迁移本质上是晶粒间原子的迁移, 因此假设晶界迁移率和自体扩散系数成正比. 在此定义辐照增强晶界迁移因子(R)为

    R=Drad/DradDtheDthe,
    (14)

    式中, R反映中子辐照增强晶界迁移效应, 其中DtheDrad分别是未辐照下热激活的自体扩散系数和中子辐照下的自体扩散系数. DtheDrad的表达式为[13,39]

    Dthe=DICeqI+DVCeqV,
    (15)
    Drad=DICI+DI2CI2+DVCV,
    (16)

    其中DI, DI2, DV分别是I, I2, V的扩散系数, 由(17)式给出. CeqICeqV分别是I和V的热平衡密度, 由(18)式给出. CI, CI2CV分别是辐照下I, I2和V的密度. D0θ, EmθEfθ分别为材料相关的扩散系数前置因子、可移动缺陷迁移能和可移动缺陷形成能, 具体数值见补充材料. Ω=a3/2为体心立方-钨的原子体积:

    Dθ=D0θexp(EmθkT),
    (17)
    Ceqθ=exp(EfθkT)/Ω.
    (18)

    未辐照钨等温退火实验数据拟合结果表明钨的晶界迁移率M满足阿伦尼乌斯的形式[7]. 通过引入辐照增强扩散因子, 中子辐照下晶界迁移率可以写成下列形式[40,41]:

    M=RM0Texp(QmkT),
    (19)

    其中M0为材料常数, Qm为再结晶激活能, T为温度, k为玻尔兹曼常数. 对(7)式进行微分得到辐照下再结晶分数随时间演化的微分表达式:

    dXdt=naURM0Texp[Qm/(kT)]naP(1X)×[ln(1X)](11/na).
    (20)

    为得到钨再结晶过程中的材料常数并与HFIR堆中子辐照钨(80%热轧变形量, 记为W80)实验结果进行比较, 通过对比文献[7, 8, 13]中W67(67%温轧变形量, 记为W67)和W80的显微组织和维氏硬度, 发现二者初始微观结构和力学性质相近, 因此这里可认为二者具有相似的再结晶动力学过程. 图1为使用JMAK模型拟合未辐照W67等温退火实验数据得到的结果[7]. 从图1可以看出, 除了再结晶的初始阶段, 拟合结果和实验结果基本符合. 拟合得到的材料常数如下:

    图 1 使用JMAK模型拟合钨等温退火实验再结晶分数的演变, 拟合用实验数据取自Lopez[7]\r\nFig. 1. Using the JMAK model to fit the evolution of the recrystallization fraction in the isothermal annealing experiment of pure tungsten, the experimental data used for the fitting were taken from Lopez[7].
    图 1  使用JMAK模型拟合钨等温退火实验再结晶分数的演变, 拟合用实验数据取自Lopez[7]
    Fig. 1.  Using the JMAK model to fit the evolution of the recrystallization fraction in the isothermal annealing experiment of pure tungsten, the experimental data used for the fitting were taken from Lopez[7].
    naUM0=5.02×1015Km3/J,
    EVold=6.04×105J/m3, Qm=6.08eV, na2.55.

    CD模型以HFIR堆中子辐照为例, 模拟了纯钨内部缺陷的演化. 图2为不同温度下缺陷演化的结果, 表明Vn和In的尺寸和密度对温度和时间存在不同的依赖关系. 从图2(a), (c), (e), (g), (i), (k)可以看出随着温度的上升, Vn的最大尺寸起初也随之增大, 在1100 ℃左右达到峰值, 之后最大尺寸急剧下降, 但其整体密度一直是呈降低的趋势. Vn随时间的演变与温度密切相关, 在低于1100 ℃时, Vn随时间一直生长; 高于1100 ℃, Vn很难生长起来. 这主要是由于V的迁移能(1.66 eV)较大的缘故. 根据扩散捕获速率理论, Vn与V发生吸收和发射反应的速率系数如下[32]:

    图 2 不同中子辐照温度下的缺陷团簇尺寸与密度随辐照时间的演变 (a), (c), (e), (g), (i), (k)分别为V团簇在750, 850, 950, 1100, 1200, 1300 ℃下的演变; (b), (d), (f), (h), (j), (l) 分别为I团簇在750, 850, 950, 1100, 1200, 1300 ℃下的演变\r\nFig. 2. Evolution of defect cluster size and density with irradiation time at different neutron irradiation temperatures: (a), (c), (e), (g), (i), (k) Evolution of V cluster at 750, 850, 950, 1100, 1200, 1300 ℃; (b), (d), (f), (h), (j), (l) I clusters evolution of cluster at 750, 850, 950, 1100, 1200, 1300 ℃.
    图 2  不同中子辐照温度下的缺陷团簇尺寸与密度随辐照时间的演变 (a), (c), (e), (g), (i), (k)分别为V团簇在750, 850, 950, 1100, 1200, 1300 ℃下的演变; (b), (d), (f), (h), (j), (l) 分别为I团簇在750, 850, 950, 1100, 1200, 1300 ℃下的演变
    Fig. 2.  Evolution of defect cluster size and density with irradiation time at different neutron irradiation temperatures: (a), (c), (e), (g), (i), (k) Evolution of V cluster at 750, 850, 950, 1100, 1200, 1300 ℃; (b), (d), (f), (h), (j), (l) I clusters evolution of cluster at 750, 850, 950, 1100, 1200, 1300 ℃.
    γ+n=4πrVnDV,
    (21)
    γn+1=γ+nΩexp(EbVn-VkT),
    (22)

    其中rVn为Vn的捕获半径; EbVn-V为Vn与V反应的结合能, 具体数值见补充材料.

    图3为由(21)式和(22)式得到的不同温度下Vn与V反应的速率系数, 而反应通量J±Vn, V正比于相应的速率系数. 从图3可知, 温度高于800 ℃, 吸收反应的速率系数γ+Vn就已经很大. 在750—1100 ℃之间, 随着温度的升高, Vn与V反应的吸收通量增大, 且吸收通量大于相应的发射反应通量, 导致Vn的生长速率与温度呈正相关. 由(22)式可知, 随着温度的升高, 发射反应的速率系数γVnγ+Vn增大得更快, 当温度超过1100 ℃后, Vn的发射反应通量占主导, 导致其平均尺寸和密度随温度的升高而急剧下降.

    图 3 不同温度下Vn与V发生反应的速率系数 (a)不同温度下Vn吸收V反应的速率系数; (b)不同温度下Vn+1发射V反应的速率系数\r\nFig. 3. Rate coefficients of the reaction between Vn and V at different temperatures: (a) Rate coefficients of Vn absorption V reactions at different temperatures; (b) rate coefficients of Vn+1 emission V reactions at different temperatures.
    图 3  不同温度下Vn与V发生反应的速率系数 (a)不同温度下Vn吸收V反应的速率系数; (b)不同温度下Vn+1发射V反应的速率系数
    Fig. 3.  Rate coefficients of the reaction between Vn and V at different temperatures: (a) Rate coefficients of Vn absorption V reactions at different temperatures; (b) rate coefficients of Vn+1 emission V reactions at different temperatures.

    有别于Vn, 图2(b), (d), (f), (h), (j)表明In的平均尺寸与整体密度都是随温度的上升而下降, 且在整个温度区间内其尺寸分布很快达到平衡且尺寸较小. 这主要是由于I的迁移能(0.013 eV)远低于V的迁移能. 在所考虑的温度区间内, I的扩散系数足够大, 导致In的分布很快达到平衡. 相较于In与I反应的吸收反应通量, 相应的发射反应通量随温度增大得更快, 导致其平均尺寸和密度随温度的升高而降低.

    上述结果与Klimenkov等[42]在裂变中子辐照钨实验上观察到的结果相比, 实验观察到空洞与析出物在1100 ℃时达到最大尺寸, 且空洞密度也基本随温度上升而下降, 这与本模拟结果一致; 但当辐照温度为1200 ℃时, 空洞密度不降反升, 另外1100 ℃时的位错环尺寸大于1000 ℃和1200 ℃, 这有别于本研究的模拟结果. 这可能是因为本研究没有考虑嬗变元素对缺陷演化的影响, 已有的实验和模拟研究表明不同密度的Re和Os会改变空洞与位错环的密度[17,43,44]. 此外, Vn并非只形成三维空洞, 部分Vn也会形成V位错环[18,45]. 因此, 模拟结果和实验有所出入.

    图4为基于3.1节CD模拟结果所得到的辐照增强P和辐照增强M, 其中考虑的最长时间为完全再结晶或辐照两年的时间. 图4(a)为不同温度下P随时间的演变, 可看出P与辐照温度呈负相关. 这是由于随着温度的上升, I和V扩散系数增大, 一方面使缺陷间复合增加, 另一方面使大尺寸Vn越容易生长, 进而吸收更多的点缺陷使得缺陷团簇密度下降, 而P与缺陷团簇密度成正相关, 因此当温度上升, 缺陷密度下降, P也跟着降低. 此外, 1100 ℃下P随着辐照时间成指数增大, 这与图2(g)中Vn的演变有直接关系. 为了说明这一点, 图4(b)展示了在所截取的时间点处, 不同缺陷类型对P的贡献占比. 可以明显看出低于1100 ℃下P主要由空洞贡献, 高于1100 ℃下P主要位错线贡献. 这是由于当辐照温度在1100 ℃以下, Vn一直生长, 所形成的空洞对P的贡献占比随时间增大; 而当辐照温度升高至1200 ℃及以上, Vn的平均尺寸和密度开始急剧下降, 且In的平均尺寸和密度也会降低, 因此P变成几乎全部都由位错线贡献. 值得一提的是, Vn在1100 ℃下生长得最快, 尤其是在演化后期出现更大尺寸的缺陷团簇, 虽然密度有所下降, 但尺寸相对于In高出2个数量级, 因此P在1100 ℃下的演化后期急剧上升.

    图 4 (a)不同辐照温度下P的演变; (b)在截取的时间点处不同辐照温度下P的占比; (c)不同辐照温度下R的演变; (d)不同辐照温度下P×M的演变\r\nFig. 4. (a) Evolution of driving force P at different irradiation temperatures; (b) proportion of defects contribution to driving force at different irradiation temperatures and time; (c) evolution of R at different irradiation temperatures; (d) evolution of product of driving force and grain boundary mobility (P×M) at different irradiation temperatures.
    图 4  (a)不同辐照温度下P的演变; (b)在截取的时间点处不同辐照温度下P的占比; (c)不同辐照温度下R的演变; (d)不同辐照温度下P×M的演变
    Fig. 4.  (a) Evolution of driving force P at different irradiation temperatures; (b) proportion of defects contribution to driving force at different irradiation temperatures and time; (c) evolution of R at different irradiation temperatures; (d) evolution of product of driving force and grain boundary mobility (P×M) at different irradiation temperatures.

    图4(c)所示, 与P随辐照温度变化的趋势相同, 中子辐照下的R与温度一直成负相关. 这是因为低温时点缺陷热平衡密度远低于辐照产生的点缺陷密度, 使得Dthe远低于Drad, 进而温度越低R越大. 因此可认为辐照温度越低, M越容易受到辐照增强自体扩散效应的影响[13]. 但随着温度逐渐升高, 点缺陷热平衡密度上升, 导致Dthe快速增大. 模拟结果表明辐照下Drad也会随温度的上升而上升, 但其上升的幅度比Dthe上升的幅度小, 造成R逐渐降低. 且在辐照温度升至1300 ℃时, R更是接近于1, 因此1300 ℃下辐照的M接近于未辐照的M.

    图4(d)为不同辐照温度下P×M的演变, 反映了二者对中子辐照诱导再结晶过程的综合影响. 钨的晶界迁移由于是热激活过程((19)式), 其随温度的上升而指数上升, 因此P×M基本上是随温度的上升而增大. 但在1100 ℃下, 由于P在演化后期急剧上升, 因此出现P×M > 1200 ℃的现象.

    图5为只考虑辐照增强P效应以及辐照增强P效应和M效应都考虑时再结晶分数随时间的演变曲线. 可以看出, 若只考虑辐照增强P, 当辐照温度低于850 ℃时, 两年内都不会发生明显的再结晶现象, 且1200 ℃下需要约171 h达到一半再结晶分数. 相比之下, 若中子辐照增强PM效应都考虑下, 850 ℃时纯钨的tX=0.5只有约503 h, 而1200 ℃下只需约27 h就能达到一半再结晶分数, 这比只考虑辐照增强P条件下快了144 h. 根据Lopez[7]的未辐照下等温退火实验结果, 1200 ℃纯钨的tX=0.5约为300 h, 因此, 可以认为中子辐照缺陷导致PM增大将使钨在1200 ℃下的tX=0.5缩短约273 h.

    图 5 只考虑辐照增强P及同时考虑辐照增强P和M (P+M)的再结晶分数演变曲线\r\nFig. 5. Recrystallization fraction (X) evolution curve considering only irradiation enhancement on driving force (P) and both irradiation enhancement driving force and grain boundary mobility (P+M).
    图 5  只考虑辐照增强P及同时考虑辐照增强PM (P+M)的再结晶分数演变曲线
    Fig. 5.  Recrystallization fraction (X) evolution curve considering only irradiation enhancement on driving force (P) and both irradiation enhancement driving force and grain boundary mobility (P+M).

    当将辐照增强PM都考虑下850 ℃的模拟结果与HFIR堆850 ℃下中子辐照钨24 d (576 h)达到一半再结晶分数的结果[13]进行对比可以发现, 模拟得到的tX=0.5 (503 h)与实验基本符合但略小些. 实验的tX=0.5略长的原因初步推测是HFIR堆高热中子通量辐照下产生的嬗变元素(Re和Os含量分别达到3.3%和1.9%原子分数[13]). 这些嬗变元素会形成纳米析出颗粒, 从会钉扎住晶界, 产生Zener阻力, 降低晶界迁移速率[33,46]. 因此, HFIR堆850 ℃下的实验结果相较于模拟结果而言有较长的tX=0.5.

    同时辐照增强PM时, 虽然M随着温度指数增长, 但1100 ℃和1200 ℃的完全再结晶时间却几乎相同. 这是由于在1100 ℃辐照温度下的演化后期, PM的乘积大幅增大(图4(d)), 导致后期再结晶速率加快. 另一个值得关注的现象是在1300 ℃的辐照温度下, 只考虑辐照增强P和辐照增强PM都考虑下的再结晶曲线几乎重合. 这是由于高温下热激活扩散效应占主导, 使R趋近为1, 因此二者的再结晶速率差别不大.

    图6为不同温度下基于热激活JMAK模型得到的tX=0.5和IIR模型得到的tX=0.5, 可以明显看出温度越高辐照与未辐照的tX=0.5量级相差越小. 相较于未辐照下的热激活再结晶过程, 750 ℃中子辐照下的tX=0.5缩短了8个量级, 1000 ℃下的缩短了5个量级, 而1300 ℃下的则相差无几. 这是由于随着温度的升高, 热激活主导了晶界迁移, 且温度高于1100 ℃后辐照缺陷密度和尺寸大幅下降, 使得辐照诱导再结晶现象基本消失.

    图 6 中子辐照下与未辐照下钨的半再结晶时间(${t_{X = 0.5}}$)随温度的演变\r\nFig. 6. Evolution of semi-recrystallization time (${t_{X = 0.5}}$) of tungsten under neutron irradiation and non-irradiation with temperature.
    图 6  中子辐照下与未辐照下钨的半再结晶时间(tX=0.5)随温度的演变
    Fig. 6.  Evolution of semi-recrystallization time (tX=0.5) of tungsten under neutron irradiation and non-irradiation with temperature.

    本研究通过建立CD模型和IIR模型模拟了不同温度下中子辐照钨缺陷团簇和再结晶分数的演化, 其中模拟的辐照缺陷尺寸和时间尺度分别达到十几纳米和年的量级, 结果包括以下3点.

    1) CD模拟结果表明不同温度下I团簇和V团簇随时间演变的趋势不同, 这主要是由于二者的迁移性质不同所导致的. 其中最突出的特征表现为V团簇尺寸在辐照温度为1100 ℃左右时达到极值, 这与中子辐照实验观察到的结果一致. 综合辐照增强PM的效应, 导致1100 ℃下完全再结晶时间和1200 ℃下的几乎相同.

    2) 本文在辐照诱导驱动力的基础上, 通过引入辐照增强晶界迁移因子R建立了新的IIR模型. 模拟结果与最近HFIR堆850 ℃中子辐照实验相符, 表明辐照增强晶界迁移是影响辐照诱导再结晶过程的重要因素. 模拟计算的850 ℃下中子辐照的半再结晶时间(tX=0.5)虽然与实验结果几乎相同, 但模拟预测的tX=0.5稍短于实验结果, 这可能是由于模型忽略了嬗变析出物对缺陷演化和晶界迁移影响的缘故. 由于是首次实验上观察到钨的辐照诱导再结晶现象, 模型需要进一步的实验数据进行完善.

    3) 通过IIR模型计算了不同辐照温度下钨再结晶分数曲线, 给出了偏滤器可能运行温度区间内辐照诱导再结晶动力学的完整图像. 同时, 模型比较了辐照和未辐照下半再结晶时间. 结果表明辐照增强效应和热激活效应对温度不同的依赖关系, 低温下辐照增强效应占主导, 高温下则是由热激活效应占主导.

    另外, 由于目前的IIR模型还并不完备, 在未来将会纳入以下因素以进一步的完善模型.

    1) 纳入钨微观结构的变化和嬗变元素对缺陷演化及晶界迁移的影响, 包括位错和晶界的演变以及晶界迁移过程中, 晶界与缺陷的相互作用.

    2) 将再结晶开始前的回复阶段纳入模型. 未辐照纯钨等温退火实验1200 ℃下都观察到了明显的再结晶成核孕育阶段.

    3) 本研究假设发生的再结晶过程是位点饱和成核, 即在短时间内再结晶晶核就已经饱和, 这在高温或辐照情况下可能是如此, 但还需要进一步探究成核过程以更确切的评估中子辐照下钨再结晶动力学过程.

    4) 改进CD模型, 考虑大尺寸位错环也具有一定的移动性以及位错环与位错线的相互作用.

    [1]

    Mahan G D 1981 Many Particle Physics (New York and London: Plenum Press)

    [2]

    Pan L, Chen X, Chen Y, Zhai H 2020 Nat. Phys. 16 767Google Scholar

    [3]

    Chen Y 2021 JHEP 04 215

    [4]

    Feynman R P, Vernon Jr F L 1963 Ann. Phys. 24 118Google Scholar

    [5]

    Caldeira A O, Leggett A J 1983 Phys. A 121 587Google Scholar

    [6]

    Bouganne R, Aguilera M B, Ghermaoui A, Beugnon J, Gerbier F 2020 Nat. Phys. 16 21Google Scholar

    [7]

    Syassen N, Bauer D M, Lettner M, Volz T, Dietze D, GarcÍa- ripoll J J, Cirac J I, Rempe G, Dürr S 2008 Science 320 1329

    [8]

    Barontini G, Labouvie R, Stubenrauch F, Vogler A, Guarrera V, Ott H 2013 Phys. Rev. Lett. 110 035302Google Scholar

    [9]

    Tomita T, Nakajima S, Danshita I, Takasu Y, Takahashi Y 2017 Sci. Adv. 99 e1701513

    [10]

    Sponselee K, Freystatzky L, Abeln B, Diem M, Hundt B, Kochanke A, Ponath T, Santral B, Mathey L, Sengstock K, and Becker C 2018 Quantum Sci. Technol. 4 014002Google Scholar

    [11]

    Li J, Harter A K, Liu J, de Moel L, Joglekar Y N, Luo L 2019 Nat. Commun. 10 855Google Scholar

    [12]

    Tomita T, Nakajima S, Takasu Y, Takahashi Y 2019 Phys. Rev. A 99 031601Google Scholar

    [13]

    Cai Z, Barthel T 2013 Phys. Rev. Lett. 111 150403Google Scholar

    [14]

    Poletti D, Barmettler P, Georges A, Kollath C 2013 Phys. Rev. Lett. 111 195301Google Scholar

    [15]

    Sachdev S 1999 Quantum Phase Transitions (Cambridge: Cambridge University Press)

    [16]

    Bardeen J M, Carter B, Hawking S 1973 Commun. Math. Phys. 31 161Google Scholar

    [17]

    Hawking S 1971 Phys. Rev. Lett. 26 1344Google Scholar

    [18]

    Hawking S 1974 Nature 248 30Google Scholar

    [19]

    Hawking S 1975 Commun. Math. Phys. 43 199Google Scholar

    [20]

    Christodoulou D 1970 Phys. Rev. Lett. 25 1596Google Scholar

    [21]

    Christodoulou D, Ruffini R 1971 Phys. Rev. D 4 3552Google Scholar

    [22]

    Bekenstein J D 1972 Lett. Nuovo Cimento 4 737Google Scholar

    [23]

    Bekenstein J D 1973 Phys. Rev. D 7 2333

    [24]

    Hawking S 1976 Phys. Rev. D 14 2460Google Scholar

    [25]

    Page D N 1993 Phys. Rev. Lett. 71 3743Google Scholar

    [26]

    Strominger A, Vafa C 1996 Phys. Lett. B 379 99Google Scholar

    [27]

    Dabholkar A, Gomes J, Murthy S 2015 JHEP 03 074

    [28]

    Gubser S S, Klebanov I R, Polyakov A M 1998 Phys. Lett. B 428 105Google Scholar

    [29]

    Witten E 1998 Adv. Theor. Math. Phys. 2 253Google Scholar

    [30]

    Maldacena J 1999 Int. J. Theor. Phys. 38 1113Google Scholar

    [31]

    Kitaev A talk at KITP 2015 http://online.kitp.ucsb.edu/online/entangled15/kitaev/ [2021-9-1]

    [32]

    Kitaev A talk at KITP 2015 http://online.kitp.ucsb.edu/online/entangled15/kitaev2/ [2021-9-1]

    [33]

    Ye J, Sachdev S 1993 Phys. Rev. Lett. 70 3339Google Scholar

    [34]

    Maldacena J, Stanford D 2016 Phys. Rev. D 94 106002Google Scholar

    [35]

    Kitaev A, Suh S J 2018 JHEP 05 183

    [36]

    Maldacena J, Stanford D, Yang Z 2016 PTEP 2016 12C104

    [37]

    Jakiw R 1985 Nucl. Phys. B 252 343Google Scholar

    [38]

    Teitelboim C 1983 Phys. Lett. B 126 41Google Scholar

    [39]

    Gross D J, Rosenhaus V 2017 JHEP 05 092

    [40]

    Gu Y, Lucas A, Qi X L 2017 JHEP 09 120

    [41]

    Huang Y, Gu Y 2017 Phys. Rev. D 100 041901

    [42]

    Liu C, Chen X, Balents L 2018 Phys. Rev. B 97 245126Google Scholar

    [43]

    Zhang P 2019 Phys. Rev. B 100 245104Google Scholar

    [44]

    Chen Y, Qi X L, Zhang P 2020 JHEP 06 121

    [45]

    Zhang P, Liu C, Chen X 2020 SciPost Phys. 8 094Google Scholar

    [46]

    Qi X L, Zhang P 2020 JHEP 05 129

    [47]

    Haldar A, Bera S, Banerjee S 2020 Phys. Rev. Res. 2 033505Google Scholar

    [48]

    Zhang P 2020 JHEP 06 143

    [49]

    Penington G 2020 JHEP 09 002

    [50]

    Almheiri A, Engelhardt N, Marolf D, Maxfield H 2019 JHEP 12 063

    [51]

    Almheiri A, Mahajan R, Maldacena J, Zhao Y 2020 JHEP 03 149

    [52]

    Dadras P, Kitaev A 2021 JHEP 03 198

    [53]

    Su K, Zhang P, Zhai H 2021 JHEP 06 156

    期刊类型引用(1)

    1. 李旭,廖晨伦,谢亚红,曹小岗,洪兵,梁立振,顾玉明,胡纯栋,孟献才. 基于全永磁微波离子源钛自成靶的D-D中子源设计研究. 真空与低温. 2024(03): 336-342 . 百度学术

    其他类型引用(0)

  • 图 1  耗散费曼图图形规则演示图

    Fig. 1.  Illustrations of the diagram rules of the dissipative Feynman diagrams.

    图 2  (10)式的图形表达. 这一图形法则可以用于高阶图的展开

    Fig. 2.  Diagram expressions of Eq. (10).

    图 3  第二Renyi熵的指数的耗散费曼图

    Fig. 3.  Dissipative diagrams of the exponential of the second Renyi entropy.

    图 4  (a)红色为0动量粒子数衰变的数据, 蓝色为动量空间中的粒子数宽度随时间变化的曲线; (b)拟合的不同晶格强度下的Bose-Hubbard模型中的0动量粒子数衰变曲线中的参数 η. 这里的反常维度应该在量子临界点处最小. 图(b)的小图里两个箭头所在的晶格强度就是量子临界点所在的位置[2]

    Fig. 4.  (a) Red curve shows the decay of zero momentum particle occupation. The blue curve shows how the width of particle momentum distribution evolutes over time. The solid line is our theoretical prediction. (b) Theoretical curve with experimental data for zero momentum particle decay for different parameters. In the inset figure, we shows the anomalous dimension eta extracted from experimental data and we can see its minimal being around quantum critical region[2].

    图 5  当环境比SYK系统温度高时系统的熵随时间的变化 (a)环境和系统的谱函数; (b)系统的熵随时间的演化, 可以看到经过初期的平方增长后变为线性增大. 系统的温度β=10, 相互作用强度 J=4; 环境温度βE=2, 相互作用强度 J=8. 该图引用自文献[3]

    Fig. 5.  (a)Spectral functions of the environment and the system; (b) entropy dynamics after the quench interaction between the system and the environment. Here the temperature of the environment is higher than the system’s initial temperature. The entropy dynamics shows a typical thermalization case. Cited from Ref. [3].

    图 6  (a)系统的谱和环境的谱; (b)“黑洞蒸发”的熵的变化图. 系统的温度β=2, 相互作用强度 J=2; 环境的温度βE=20, 相互作用强度 J=2. 该图引自文献[3]

    Fig. 6.  (a) Spectral functions of the environment and the system; (b) entropy dynamics after the quench interaction between the system and the environment. Here the temperature of the environment is lower than the system’s initial temperature. The entropy dynamics shows a typical cooling case. Here it can be compared with “Black Hole Evaporation” and the entropy dynamics looks like a Page curve. Cited from Ref. [3].

  • [1]

    Mahan G D 1981 Many Particle Physics (New York and London: Plenum Press)

    [2]

    Pan L, Chen X, Chen Y, Zhai H 2020 Nat. Phys. 16 767Google Scholar

    [3]

    Chen Y 2021 JHEP 04 215

    [4]

    Feynman R P, Vernon Jr F L 1963 Ann. Phys. 24 118Google Scholar

    [5]

    Caldeira A O, Leggett A J 1983 Phys. A 121 587Google Scholar

    [6]

    Bouganne R, Aguilera M B, Ghermaoui A, Beugnon J, Gerbier F 2020 Nat. Phys. 16 21Google Scholar

    [7]

    Syassen N, Bauer D M, Lettner M, Volz T, Dietze D, GarcÍa- ripoll J J, Cirac J I, Rempe G, Dürr S 2008 Science 320 1329

    [8]

    Barontini G, Labouvie R, Stubenrauch F, Vogler A, Guarrera V, Ott H 2013 Phys. Rev. Lett. 110 035302Google Scholar

    [9]

    Tomita T, Nakajima S, Danshita I, Takasu Y, Takahashi Y 2017 Sci. Adv. 99 e1701513

    [10]

    Sponselee K, Freystatzky L, Abeln B, Diem M, Hundt B, Kochanke A, Ponath T, Santral B, Mathey L, Sengstock K, and Becker C 2018 Quantum Sci. Technol. 4 014002Google Scholar

    [11]

    Li J, Harter A K, Liu J, de Moel L, Joglekar Y N, Luo L 2019 Nat. Commun. 10 855Google Scholar

    [12]

    Tomita T, Nakajima S, Takasu Y, Takahashi Y 2019 Phys. Rev. A 99 031601Google Scholar

    [13]

    Cai Z, Barthel T 2013 Phys. Rev. Lett. 111 150403Google Scholar

    [14]

    Poletti D, Barmettler P, Georges A, Kollath C 2013 Phys. Rev. Lett. 111 195301Google Scholar

    [15]

    Sachdev S 1999 Quantum Phase Transitions (Cambridge: Cambridge University Press)

    [16]

    Bardeen J M, Carter B, Hawking S 1973 Commun. Math. Phys. 31 161Google Scholar

    [17]

    Hawking S 1971 Phys. Rev. Lett. 26 1344Google Scholar

    [18]

    Hawking S 1974 Nature 248 30Google Scholar

    [19]

    Hawking S 1975 Commun. Math. Phys. 43 199Google Scholar

    [20]

    Christodoulou D 1970 Phys. Rev. Lett. 25 1596Google Scholar

    [21]

    Christodoulou D, Ruffini R 1971 Phys. Rev. D 4 3552Google Scholar

    [22]

    Bekenstein J D 1972 Lett. Nuovo Cimento 4 737Google Scholar

    [23]

    Bekenstein J D 1973 Phys. Rev. D 7 2333

    [24]

    Hawking S 1976 Phys. Rev. D 14 2460Google Scholar

    [25]

    Page D N 1993 Phys. Rev. Lett. 71 3743Google Scholar

    [26]

    Strominger A, Vafa C 1996 Phys. Lett. B 379 99Google Scholar

    [27]

    Dabholkar A, Gomes J, Murthy S 2015 JHEP 03 074

    [28]

    Gubser S S, Klebanov I R, Polyakov A M 1998 Phys. Lett. B 428 105Google Scholar

    [29]

    Witten E 1998 Adv. Theor. Math. Phys. 2 253Google Scholar

    [30]

    Maldacena J 1999 Int. J. Theor. Phys. 38 1113Google Scholar

    [31]

    Kitaev A talk at KITP 2015 http://online.kitp.ucsb.edu/online/entangled15/kitaev/ [2021-9-1]

    [32]

    Kitaev A talk at KITP 2015 http://online.kitp.ucsb.edu/online/entangled15/kitaev2/ [2021-9-1]

    [33]

    Ye J, Sachdev S 1993 Phys. Rev. Lett. 70 3339Google Scholar

    [34]

    Maldacena J, Stanford D 2016 Phys. Rev. D 94 106002Google Scholar

    [35]

    Kitaev A, Suh S J 2018 JHEP 05 183

    [36]

    Maldacena J, Stanford D, Yang Z 2016 PTEP 2016 12C104

    [37]

    Jakiw R 1985 Nucl. Phys. B 252 343Google Scholar

    [38]

    Teitelboim C 1983 Phys. Lett. B 126 41Google Scholar

    [39]

    Gross D J, Rosenhaus V 2017 JHEP 05 092

    [40]

    Gu Y, Lucas A, Qi X L 2017 JHEP 09 120

    [41]

    Huang Y, Gu Y 2017 Phys. Rev. D 100 041901

    [42]

    Liu C, Chen X, Balents L 2018 Phys. Rev. B 97 245126Google Scholar

    [43]

    Zhang P 2019 Phys. Rev. B 100 245104Google Scholar

    [44]

    Chen Y, Qi X L, Zhang P 2020 JHEP 06 121

    [45]

    Zhang P, Liu C, Chen X 2020 SciPost Phys. 8 094Google Scholar

    [46]

    Qi X L, Zhang P 2020 JHEP 05 129

    [47]

    Haldar A, Bera S, Banerjee S 2020 Phys. Rev. Res. 2 033505Google Scholar

    [48]

    Zhang P 2020 JHEP 06 143

    [49]

    Penington G 2020 JHEP 09 002

    [50]

    Almheiri A, Engelhardt N, Marolf D, Maxfield H 2019 JHEP 12 063

    [51]

    Almheiri A, Mahajan R, Maldacena J, Zhao Y 2020 JHEP 03 149

    [52]

    Dadras P, Kitaev A 2021 JHEP 03 198

    [53]

    Su K, Zhang P, Zhai H 2021 JHEP 06 156

  • [1] 陈若凡. 时间演化矩阵乘积算符方法及其在量子开放系统中的应用. 物理学报, 2023, 72(12): 120201. doi: 10.7498/aps.72.20222267
    [2] 李蕊轩, 张勇. 熵在非晶材料合成中的作用. 物理学报, 2017, 66(17): 177101. doi: 10.7498/aps.66.177101
    [3] 徐红梅, 金永镐, 郭树旭. 电压控制不连续导电模式DC-DC变换器的熵特性研究. 物理学报, 2013, 62(24): 248401. doi: 10.7498/aps.62.248401
    [4] 冯维, 丁辉, 林昊, 罗辽复. λ噬菌体溶源/裂解转换调控与定态熵. 物理学报, 2012, 61(16): 168701. doi: 10.7498/aps.61.168701
    [5] 潘欣裕, 赵鹤鸣. Logistic混沌系统的熵特性研究. 物理学报, 2012, 61(20): 200504. doi: 10.7498/aps.61.200504
    [6] 杨 波. 一般加速带电带磁的动态黑洞中标量场的熵. 物理学报, 2008, 57(4): 2614-2620. doi: 10.7498/aps.57.2614
    [7] 杨 波. 变加速直线运动黑洞的温度和Dirac场的熵. 物理学报, 2007, 56(11): 6772-6776. doi: 10.7498/aps.56.6772
    [8] 韩亦文, 洪 云, 杨树政. 广义不确定关系与整体单极黑洞Dirac场的熵. 物理学报, 2007, 56(1): 10-14. doi: 10.7498/aps.56.10
    [9] 郑元强. 球对称动态黑洞Dirac场的熵的再讨论. 物理学报, 2007, 56(3): 1266-1270. doi: 10.7498/aps.56.1266
    [10] 郑元强. 动态广义球对称含荷黑洞Dirac场的熵. 物理学报, 2006, 55(7): 3272-3276. doi: 10.7498/aps.55.3272
    [11] 牛振风, 刘文彪. 新Tortoise坐标变换与任意加速带电动态黑洞熵. 物理学报, 2005, 54(1): 475-480. doi: 10.7498/aps.54.475
    [12] 孙学锋, 景 玲, 刘文彪. 黑洞熵无截断薄层模型的改进与推广. 物理学报, 2004, 53(11): 4002-4006. doi: 10.7498/aps.53.4002
    [13] 强丽娥, 高新芹, 赵 峥. 动态黑洞温度和熵的再讨论. 物理学报, 2004, 53(10): 3619-3626. doi: 10.7498/aps.53.3619
    [14] 韩亦文, 洪 云. Schwarzschild-de-Sitter黑洞宇宙视界量子态的熵. 物理学报, 2004, 53(10): 3270-3273. doi: 10.7498/aps.53.3270
    [15] 王波波. 环面黑洞背景下量子场的熵. 物理学报, 2004, 53(7): 2401-2406. doi: 10.7498/aps.53.2401
    [16] 孙鸣超. 起源于引力场的Vaidya-Bonner-de Sitter黑洞的量子熵. 物理学报, 2003, 52(6): 1350-1353. doi: 10.7498/aps.52.1350
    [17] 宋太平, 侯晨霞, 黄金书. 一般球对称带电蒸发黑洞的熵. 物理学报, 2002, 51(8): 1901-1906. doi: 10.7498/aps.51.1901
    [18] 张靖仪, 赵峥. 直线加速动态黑洞Dirac场的熵. 物理学报, 2002, 51(10): 2399-2406. doi: 10.7498/aps.51.2399
    [19] 宋太平, 侯晨霞, 史旺林. Vaidya-Bonner黑洞的熵. 物理学报, 2002, 51(6): 1398-1402. doi: 10.7498/aps.51.1398
    [20] 魏志勇, 段利敏, 吴和宇, 靳根明, 李祖玉, 诸永泰, 郗洪飞, 沈文庆, 肖志刚, 王宏伟, 张保国, 柳永英, 王素芳, 胡荣江. 35MeV/u40Ar+197Au中的熵产生. 物理学报, 2001, 50(4): 649-654. doi: 10.7498/aps.50.649
  • 期刊类型引用(1)

    1. 李旭,廖晨伦,谢亚红,曹小岗,洪兵,梁立振,顾玉明,胡纯栋,孟献才. 基于全永磁微波离子源钛自成靶的D-D中子源设计研究. 真空与低温. 2024(03): 336-342 . 百度学术

    其他类型引用(0)

计量
  • 文章访问数:  6425
  • PDF下载量:  324
  • 被引次数: 1
出版历程
  • 收稿日期:  2021-09-10
  • 修回日期:  2021-11-01
  • 刊出日期:  2021-12-05

/

返回文章
返回