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

基于蒙特卡罗-离散纵标方法的氘氚激光等离子体聚变反应率数值模拟

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

摘要:惯性约束聚变(ICF)是实现受控热核聚变可能途径之一. 聚变中子源项是氘氚激光等离子体物理设计与分析的重要参数之一, 其准确性直接影响分析结果的可靠性. 目前国内外对于ICF氘氚聚变反应产生的中子源项研究主要基于解析公式法, 在温度和反应类型等方面适用范围有限. 本文采用粒子云概念对氘、氚粒子云团开展了随机抽样与时空网格划分, 然后基于麦克斯韦速率分布律对氘氚聚变反应开展了多普勒能量展宽效应分析与微分截面温度修正工作, 耦合蒙特卡罗方法和离散纵标方法, 开展了激光等离子体中D-T粒子云团聚变反应率的数值模拟工作. 研究结果显示, 与原核数据库截面相比, D-T, D-D, T-D截面经修正后多普勒温度效应显著. 在20—100 keV的等离子体温度范围内, 相较传统的解析公式法, 本文模拟结果更符合最新的ENDF核数据库的氘氚反应截面数据, 且与解析公式法结果在低能区存在较大误差, 可能是计算方法不同与低温截面差异过大原因导致.
关键词: 激光等离子体/
氘氚聚变反应/
微分截面温度修正/
蒙特卡罗方法/
离散纵标方法

English Abstract


--> --> -->
惯性约束聚变(ICF)是实现受控热核聚变可能途径之一[1,2]. 点火是研究惯性约束核聚变的关键, 所获得的聚变中子源项是激光等离子体物理设计与分析的重要参数之一, 其准确性直接影响分析结果的可靠性[314]. 目前国内外对于ICF氘氚聚变反应产生的中子源项研究主要基于解析公式法, 如对于氘氚各占50%密度的情况, 其产生率计算公式为
$\begin{split} & S({n_i} ,{T_i}) = n_i^2\overline {\sigma {v_{{\rm{DT}}}}} ({T_i})\\ & (0 \leqslant T \leqslant 100\;{\rm{keV}},\;{n_i} = {{10}^{14}}{\rm{c}}{{\rm{m}}^{ - 3}})\;,\end{split}$
其中, S(ni, Ti)为中子产生率, ni为氘或者氚粒子密度, Ti为等离子体温度, σ为反应截面, vDT为麦克斯韦分布下的氘氚粒子相对速度, $\overline {\sigma {v_{{\rm{DT}}}}} \;$表示两者的加权平均[15].
目前大量使用的解析公式法主要针对氘氚粒子密度相同且温度不大于100 keV, 在温度和反应类型等方面适用范围有限. 本文通过结合蒙特卡罗(MC)方法和离散纵标(SN)方法, 开展了激光等离子体中D-T粒子云团聚变反应率的数值模拟算法程序设计.
2
2.1.物理模型与碰撞物理过程
-->对于所给定的氘氚等离子体集团, 在进行模拟时, 由于数据量过于庞大, 不可能把每个粒子都进行聚变反应的模拟, 因此本文首先引入了粒子云概念, 即多个位置、速度相同的同类粒子的集合, 且认为粒子云所发生的作用均为同一反应. 其次使用了网格模型, 即将所模拟的时空进行了网格划分, 其中整个等离子体空间经网格化后构成一个三维的空间坐标, 一个空间网格中可包含多个粒子云.
氘氚等离子体为大量粒子组成, 满足麦克斯韦速率分布律, 即
$ {f(v) = 4{\text{π}}{{\left( {\frac{m}{{2{\text{π}}kT}}} \right)}^{\textstyle\frac{3}{2}}} \cdot {{\rm{e}}^{ - \textstyle\frac{{m{v^{\rm{2}}}}}{{{\rm{2}}kT}}}}{v^{\rm{2}}},} $
其中, m为粒子质量, k为玻尔兹曼常数, T为等离子团温度.
本文选取了ENDF/B-VI和ENDF/B-VII的氘氚聚变反应数据库, 其所给出的粒子能量为系统的质心系能量, 因此需要将实验室系下入射粒子与靶粒子碰撞时的对应能量转换为质心系能量. 等离子体中, 氘、氚粒子由于速度和质量不同, 用动量来进行矢量合成. 根据动量守恒
$\overline {{P}}_1 + \overline {{P}}_2 = \overline {{P}} ,$
其中, $\overline {{P}}_1 $为入射粒子的动量, $\overline {{P}}_2$为靶粒子的动量, $\overline {{P}} \;$为系统的质心系动量. 固定$\overline {{P}}_1$始终处于v轴方向, 进行矢量分解, 如图1所示.
图 1 动量矢量分解
Figure1. Momentum vector decomposition.

将(3)式写成标量形式:
${{P_1} + {P_2}{\rm{cos}}\varphi = {P_x},}$
${{P_2}{\rm{sin}}\varphi = {P_y},} $
${{P^2} = P_x^2 + P_y^2,}$
其中, φ为两矢量夹角, PxPy分别为质心系动量的x轴和y轴分量. φ角即为需要离散化处理的角度分布.
所模拟的物理过程为: 某一随机位置处, 具有某一随机速率的入射粒子云沿某一随机方向发射, 则根据核反应率的计算原理, 沿途路径上的所有网格, 均认为入射粒子云与其中所有粒子发生反应, 且按照截面大小发生核反应. 由此, 在一个时间网格内, 氘氚聚变反应的中子产生率为
${S = n_{\rm{i}} \cdot {n_{\rm{t}}} \cdot \;\sigma ' \cdot l/t,}$
其中, S为中子产生率; n i为入射粒子的数密度, n t为靶粒子的数密度, 当氘氚粒子各占50%时, n i = n t; $\sigma '$为修正后的氘氚聚变反应截面; t为模拟的时间步长; l为入射粒子云在一个时间网格内走过的路程.
2
2.2.时空网格划分与误差控制
-->氘氚等离子体发生聚变, 并实现可持续惯性热核聚变燃烧, 必须满足以下3个基本条件:
1)劳森判据条件$ {n_i\tau \geqslant {{10}^{14}}} \;{\rm{s}}/{\rm{cm}}^3$;
2)燃料等离子体温度条件$ {{T_h} \geqslant 5 \times {{10}^7}}\;{\rm{K}} $;
3) ρr乘积条件$ {{\rho _m}r \geqslant 3} \;{\rm{g}}/{\rm{cm}}^2$;
其中, ni为热核燃料等离子体密度, τ为由惯性维持该离子数密度不变的时间间隔, ρm为被压缩的低温热核燃料质密度, r为预压缩到高质密度热核燃料小球的半径[16].
为减少氘氚聚变反应率计算误差及计算量, 激光等离子体的时空网格划分极为重要. 从误差产生源来看, 假设入射粒子云所经过的网格, 均要计算其内部的粒子云与入射粒子云的反应, 但如果入射粒子云只经过网格边缘部分, 则其内部大部分粒子云并没有与入射粒子云发生碰撞, 由此导致误差. 以二维网格为例, 如图2所示.
图 2 二维网格误差分析
Figure2. Two-dimensional grid error analysis.

图2中的灰色网格代表着仅其边缘部分被入射粒子云穿过, 可以看出, 若将该类网格均计入氘氚聚变反应, 模拟结果将由此偏大. 同样, 入射粒子云的初始点和终点位置也会造成误差, 如图3所示.
图 3 二维网格大小造成的误差
Figure3. Error caused by two-dimensional mesh size.

图3给出了两个入射粒子云, 由于终点位置的不同, 导致实际路径长度差别接近一个网格. 因此, 网格划分越细, 则误差越小. 而网格大小最主要的判定条件取决于入射粒子所走的距离长短, 若网格过大, 每个时间网格(步长)所通过的路程则有可能在一两个网格内, 由此会导致巨大的误差. 根据上述的误差来源, 需要考虑两个条件: 一是时间网格大小即时间步长, 二是氘氚等离子体热运动的速率. 根据公式
${L = vt}$
计算所走路程L, 以此推出合适的网格大小设置.
对于惯性约束核聚变装置, 如美国国家点火装置, 其实际发生核反应的时间大约在10–10—10–11 s. 根据流体力学中判断计算收敛的柯朗-弗里德里奇-列维条件, 计算并设定时间步长为t = 1 ns.
对于氘氚等离子体热运动速率, 应以模拟中的最小速率为下限, 选择模拟所设定的最小温度的最概然速率作为估算网格速率的参考速率值, 则有
${{v_{\rm{m}}} = \sqrt {\frac{{2kT}}{m}},}$
其中, vmT温度下麦克斯韦速率分布的最概然速率. 根据(9)式计算, 等离子体温度为20 keV时, 氘粒子入射, 其最概然速率为1.9575 × 106 m/s; 氚粒子入射, 其最概然速率为1.1325 × 106 m/s. 估算时, 只考虑这个速率的数量级106, 由(8)式可得
${L = vt = {{10}^6}\;{\rm{m}}/{\rm{s}} \times {{10}^{ - 9}}\;{\rm{ s}} = {{10}^{ - 3}}\;{\rm{ m}},}$
即, 入射粒子云所走的路径长度在10–3 m量级以上, 为保证10%以下的误差, 网格大小应至少在10–4 m以下, 故设定网格大小为10–4 m.
2
2.3.不同反应类型截面修正
-->本文所用数据来自ENDF/B-VI和ENDF/B-VII的氘氚聚变反应数据库, 其数据由美国洛斯阿拉莫斯国家实验室(LANL)利用EDA-R矩阵码对氘氚聚变反应进行研究所生成的评价中子数据. 其中, D-T的反应数据是ENDF/B-VI和ENDF/B-VII的官方原数据, D-D的反应数据是洛斯阿拉莫斯国家实验室的一个初步结果, 有ENDF格式、HTML格式和PDF格式三种[17].
本文模拟的是氘氚等离子体聚变反应的中子产生率, 只需要考虑3个反应, 分别是
${{\rm{D}} + {\rm{D}} = {}^{{\rm{ 3}}}{\rm{He}} + {\rm{n}} + 3.27\;{\rm{MeV}},}$
${{\rm{D}} + {\rm{T}} = {}^{\rm{4}}{\rm{He}} + {\rm{n}} + 17.6\;{\rm{MeV}},}$
${{\rm{D}} + {\rm{T}} = {}^{\rm{4}}{\rm{H}}{{\rm{e}}^*} + {\rm{n}}.}$
对于给定温度, 氘氚等离子体在满足麦克斯韦速率分布的条件下, 在任何速度、任何角度的位置上都存在粒子, 因此在进行截面修正时, 应同时考虑这两个物理量对微分截面的影响. 本文修正截面的目的在于解决大数据量模拟计算时计算空间不足的问题, 将所有满足同一麦克斯韦速率分布的粒子的截面归一化为一个截面值, 即当粒子数量足够多时, 无论入射粒子与靶粒子的夹角有多大, 靶粒子的速率大小有多大, 可使用同一微分截面值.
按照麦克斯韦速率分布函数的概率做离散化处理, 即
$\overline v = \sum\limits_{v_{{\rm{h1}}}}^{{v_{{\rm{h2}}}}} vf(v)\;.$
由于不需要考虑出射后的动量方向, 所以只需要将入射粒子和靶粒子的动量矢量夹角进行归一化, 即仰角归一化. 同时, 当仰角确定时, 整个方向角的圆周上, 都应有粒子存在, 所以每一仰角所占概率, 应为仰角所对应的球带占整个球表面积的比例, 如图4所示.
图 4 仰角对应球带(球坐标)
Figure4. Elevation corresponding to the ball belt (ball coordinate).

其概率为
$\begin{split} \frac{{\Delta S}}{S} \!={}&\! \frac1{{4{\text{π}}{r^2}}} {{\displaystyle\int_\phi ^{\phi + \Delta \phi }\!{2{\text{π}}r{\rm{sin}}\varphi \cdot r{\rm{d}}\varphi } }} \\ ={}& \frac{1}{2}({\rm{cos}}\varphi \!-\! {\rm{cos}}(\varphi \!+\! \Delta \varphi )), \end{split}$
其中, φ为仰角, ?φ为离散化的角度间隔. 则归一化的仰角为
$ {\overline \varphi = \sum\limits_0^{\text{π}} {\varphi \cdot } \frac{1}{2}({\rm{cos}}\varphi - {\rm{cos}}(\varphi + \Delta \varphi )\;).} $
对于D-D, D-T, T-D反应截面, 首先计算出某一温度下粒子能够具有的速率区间[vh1, vh2], 计算入射粒子动量P1和靶粒子动量P2, 计算速率所对应的麦克斯韦速率分布概率f. 以仰角为循环变量, 从0到π, 计算碰撞后的动量P, 以之计算质心系能量, 并计算仰角对应的概率${\Delta S}/{S}$. 得到质心系能量后, 从数据库中读取对应截面σ, 利用公式
$ {\sigma ' = \sum\limits_{{v_{{\rm{h1}}}}}^{{v_{{\rm{h2}}}}} {\sum\limits_{\rm{0}}^{\text{π}} {f \times } } \frac{{\Delta S}}{S} \cdot \sigma,} $
计算出修正后的截面值, 获得对应不同入射粒子能量的矩阵.
2
2.4.路径函数与随机函数设计
-->本文所模拟的等离子体只包含氘离子和氚粒子, 且认为这些粒子均匀分布. 因为氘氚聚变反应率需要通过路径长度计算, 因此计算入射粒子云的运动路径(路径函数)较为关键. 这里对路径函数做一简要说明: 首先输入入射粒子云的初始位置$\left( {x,\; y,\; z} \right)$和入射粒子云速率的三个坐标分量$\left( {{v_x},\;{v_y},\;{v_z}} \right)$. 然后建立一个n × 3的矩阵, 其中n代表着入射粒子云走过的网格数, 3即为三个坐标值. 为避免漏算网格导致误差, 将原时间步长1 ns细分成时间网格10–12 s, 计算每一个时间网格内入射粒子云走过的网格数, 进而计算路径长度.
本文中, 除等离子体温度和氘氚粒子数密度已知外, 对于氘氚粒子云的位置、速率、方向等均随机产生. 本文采用了全部随机方法, 即无论是入射粒子云的种类、位置, 或是速率, 还是方向, 均按照一定概率进行随机抽样. 由于入射粒子云不同, 对应的麦克斯韦速率分布并不完全相同, 所以要先判断是哪种粒子, 分两种情况进行后续的部分. 对于入射氘粒子云, 先求出氘粒子云在此温度下的最概然速率vm, 并求出两个半高宽的速率范围[vh1, vh2], 然后求解麦克斯韦速率分布函数, 离散化后利用函数依概率进行随机, 得到一个速率值. 对极角也做同样的离散化处理. 对于入射氚粒子云, 与氚粒子云入射算法相同, 更改质量参数即可.
2
2.5.算法流程图
-->总的算法流程见图5.
图 5 氘氚等离子体聚变反应模拟流程图
Figure5. Simulation flow of deuterium-tritium plasma fusion reaction.

2
3.1.初始边界条件
-->为与解析公式(1)进行比较, 本文中, ni = 1014 cm–3, 氘氚粒子各占50%, 并满足劳森判据. 惯性约束等离子体时间步长通常在纳秒量级[18], 因此本文设置时间步长为1 ns. 根据2.2节的工作, 设定网格大小为10–4 m. 模拟的温度取值为20, 40, 60, 80, 100 keV. 为确保精确度, 对于每一个所模拟的温度, 粒子云速率离散为1000个网格, 极角离散为180个网格, 方位角离散为360个网格. 采用MC方法对105个网格单元进行采样, 采样网格单元作为入射粒子云网格来控制统计误差.
2
3.2.模拟结果与讨论
-->三种聚变反应截面的温度修正结果如图6所示(1 barn = 10–24 cm2).
图 6 氘氚等离子体聚变反应修正 (a) D-D反应截面修正; (b) D-T反应截面修正; (c) T-D反应截面修正
Figure6. Correction of fusion reaction of deuterium-tritium plasma: (a) D-D fusion cross section correction; (b) D-T fusion cross section correction; (c) T-D fusion cross section correction.

图6(b)图6(c)可明显地看到多普勒效应. 这是由于氘氚等离子体中, 靶粒子服从麦克斯韦速率分布, 因此在截面温度修正后, 能量将有所展宽, 且温度越高, 麦克斯韦速率分布范围越宽, 由此靶粒子展宽越大, 同时峰值截面也逐渐减小. 经截面修正, 将所有满足同一麦克斯韦速率分布的粒子的截面归一化为一个截面值, 解决了大量靶粒子热运动所带来的大数据量模拟计算时计算空间不足的问题, 大大节约了计算时间.
氘氚等离子体聚变反应模拟结果如图7所示.
图 7 氘氚等离子体聚变反应中子产生率
Figure7. Neutron production rate of deuterium-tritium plasma fusion reaction.

图7可以看出, 根据解析公式法, 随等离子体温度的上升, 中子产生率首先增加, 在60 keV左右达到峰值, 而后开始降低; 而数值模拟结果则是中子产生率随等离子体温度上升而一直增加. 根据三种聚变反应截面的温度修正即图6, 分析可知氘-氘反应、氘-氚反应、氚-氘反应的截面在20—100 keV的范围内呈递增趋势, 在此范围内, 本文采用的数值模拟法相较传统解析公式法, 更符合ENDF/B-VI和ENDF/B-VII的氘氚聚变反应数据库的截面数据趋势.
下面分析等离子体处于低温时数值模拟结果与解析公式法结果误差过大的原因. 首先, 对于氘氚聚变反应率, 解析公式法和本文采用了不同的计算方法, 具体可参见本文第2节所述. 另一方面, 解析方法所用截面来自核物理分析方法, 本文则采用了最新的ENDF核数据库. 根据图7, 解析方法所用的D-T聚变截面(聚变反应最重要的过程)在60 keV左右达到峰值, 所获得的中子产生率也在60 keV左右达到峰值. 而本文所用数值模拟方法所获得的中子产生率和ENDF/B-VII数据库的聚变截面则在20—100 keV的范围内呈递增趋势. 两种方法相比, 在约20 keV, 本文方法给出的聚变反应率远低于解析方法, 差异约70%. 图8对比分析了两种方法所对应的聚变反应截面. 根据该图, D-T聚变截面在低于60 keV的范围内出现明显差异. 就20 keV而言, 解析法使用截面为0.4077 barn, 本文方法使用截面为0.0597 barn, 差异约85%. 这可能是导致这两种方法之间出现显著差异的根本原因.
图 8 两种方法聚变反应截面对比[19,20]
Figure8. Comparison of fusion cross sections of the two methods[19,20].

本文针对惯性约束的激光等离子体, 采用MC方法和SN方法, 开展了D-T粒子云团聚变反应率的数值模拟工作. 得到以下结论: 1)与原核数据库截面相比, D-T, D-D, T-D截面经修正后多普勒温度效应显著; 2)在等离子体温度20—100 keV范围内, 本文数值模拟结果较解析公式法更符合ENDF/B-VI和ENDF/B-VII的氘氚聚变反应数据库的截面数据; 3)本文数值模拟结果与与解析公式法结果在低能区存在较大误差, 可能是计算方法不同与低温截面差异过大原因导致.
相关话题/计算 数据 数据库 激光 概率

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 结构光照明技术在二维激光诱导荧光成像去杂散光中的应用
    摘要:背景杂散光信号的干扰制约了激光片光成像技术的发展,本文将结构光照明技术应用到激光片光成像测量中来消除杂散光的干扰.介绍了基于结构光照明技术的工作原理和实验测量系统,基于Matlab软件理论分析了相位移动结构光照明技术具有完全消除杂散光的作用,并针对稳定罗丹明B溶液进行了二维激光诱导荧光成像实验 ...
    本站小编 Free考研考试 2021-12-29
  • 垂直腔面发射激光器与异质结双极型晶体管集成结构的设计和模拟
    摘要:垂直腔面发射激光器(verticalcavitysurfaceemittinglasers,VCSELs)和异质结双极型晶体管(heterojunctionbipolartransistor,HBT)都是纵向电流器件,可以集成在同一外延片上,通过HBT基极电流调制VCSELs的输出光功率.本文 ...
    本站小编 Free考研考试 2021-12-29
  • 光频梳频域干涉测距主要参数分析及一种改进的数据处理方法
    摘要:本文对光学频率梳频域干涉测距中的测距范围、分辨力、非模糊范围等的影响因素进行了分析,并说明了传统傅里叶变换法的局限性和系统误差产生原因;提出了一种等频率间隔重采样数据处理方法,该方法基于三次样条插值,修正了傅里叶变换法因频率量不等间隔造成的误差;在此基础上提出峰值位置拟合算法,解决了包络随距离 ...
    本站小编 Free考研考试 2021-12-29
  • 垂直腔面发射激光器低温光电特性
    摘要:垂直腔面发射激光器通常被用作常温下850nm波段短波长短距离光互连领域的激光光源,多在室温下进行测试和使用.在低温环境下垂直腔面发射激光器工作状态的表征是本文的研究重点.我们表征了在不同温度下直流驱动垂直腔面发射激光器的发光光谱和10%占空比脉冲电流驱动垂直腔面发射激光器的发光光谱和功率-电流 ...
    本站小编 Free考研考试 2021-12-29
  • 拉曼激光边带效应对冷原子重力仪测量精度的影响
    摘要:电光调制技术是产生拉曼光的几种方法之一,其优点是系统简单、易搭建且环境适应性强.然而,这种调制技术会产生额外的边带光,并影响冷原子干涉绝对重力仪的测量精度.本文利用自行研制的可移动冷原子重力仪,研究了边带效应对冷原子重力仪测量精度的影响.详细分析了拉曼反射镜的位置、拉曼脉冲的作用时刻及其间隔、 ...
    本站小编 Free考研考试 2021-12-29
  • 激光3D纳米打印温度敏感的微球激光器
    摘要:耳语回廊模式(WGM)微腔具有品质因子高、模式体积小、制作工艺简单多样、同时对周围环境敏感性极高等优点,已被广泛应用于传感和检测.然而,真正的尺寸可控的三维微腔却少有报道.本文报道了一种有源回音壁模式微球腔,由商业光刻胶SU-8作为腔体材料,有机染料罗丹明B作为增益介质.利用飞秒激光双光子聚合 ...
    本站小编 Free考研考试 2021-12-29
  • 不同散射介质对飞秒脉冲激光传输特性影响研究
    摘要:基于分层传输模型和Mie散射理论,在粒子散射模型中充分考虑了谱分布特征,数值模拟了800nm飞秒激光在冰云、水云、雾、气溶胶和降雨环境中的传输特性.结果表明,谱分布和粒子相态对光丝传输特性有较大的影响.雨滴的粒径较大,光丝在降雨环境中传输时,由于散射导致的能量衰减最强,产生的光丝峰值光强和能量 ...
    本站小编 Free考研考试 2021-12-29
  • 相对论皮秒激光在低密度等离子体中直接加速的电子束的横向分布特征研究
    摘要:相对论皮秒激光与低密度等离子体作用可以通过“激光直接加速”机制获得超有质动力定标率的高能电子,且电荷量可以达到百nC级,在伽马射线产生、正电子产生等方面具有重要应用.然而激光直接加速电子束相比激光尾场加速电子束具有更大的发散角,同时实验观测的横向束分布也不均匀,但是其中的物理机制研究较少.本文 ...
    本站小编 Free考研考试 2021-12-29
  • 基于二维纳米材料的超快脉冲激光器
    摘要:石墨烯以其独特的光电特性打开了二维纳米材料的大门,随后拓扑绝缘体、过渡金属硫化物、黑磷等二维材料相继被报道,这些材料由于具有良好的非线性光学特性,可用作被动饱和吸收体来产生脉冲激光.本文总结了近年来基于二维材料的光纤激光器和固体激光器的研究状况,从激光器的中心波长、脉宽、重复频率、脉冲能量和输 ...
    本站小编 Free考研考试 2021-12-29
  • 激光冷却KCl<sup>–</sup>阴离子的理论研究
    摘要:本文采用多组态相互作用方法计算了KCl–阴离子前两个离解极限K(2Sg)+Cl–(1Sg)和K(2Pu)+Cl–(1Sg)所对应的3个电子态(X2Σ+,A2Π和B2Σ+)的电子结构.在计算中考虑了Davidson修正,核-价电子关联效应及自旋轨道耦合效应.预测了KCl–阴离子的光谱常数和跃迁性 ...
    本站小编 Free考研考试 2021-12-29