审稿意见.docVIP

意见如下： 介绍中因说明镎这种元素的应用背景和产生的过程。 融合极限通过混合自由能曲线的最小值来判断来进行判断，但是仅仅选用了三个点，既镎的含量分别为0、50%、100%，选用的点有点少。 计算融合能时候，没有给出所有位置的能量、键强数据，计算初始模型的构建没有说明清楚，此外，反应物和生成物的物相（固体、液体，气体）对计算结果产生较大的影响，每个反应的计算过程中，选择什么样的状态没有给出说明，对计算结果的评估将造成一定的影响。 计算过程中，不能仅仅考虑取代阳离子的位置，NP的位置也有两种。 在一个反应过程中采用CASTEP和DMOL3相结合方法进行计算获得的能量，是否有其它的例子来验证这种方法的可行性。此外，带电体系用CASTEP计算的主要问题是什么？ 镎与铀酰根的融合：量子力学评估 摘要： 镎会融合到腐蚀的UO2变化的相中，实验表明层间带有正电阳离子的铀酰比层间不带阳离子的更容易融合镎，这种融合机制尚未明确。用密度泛函理论来比较镎融合到黄硅钾铀矿–键合的环境变化，提供电荷平衡耦合取代机制的一种直观的认识。在这项研究中，采用量子力学计算，比较了各种机制的Np+5基态的融合量。 为什么要研究镎这种元素？镎的应用背景和产生背景？？ 2 黄硅钾铀矿:发生和结构： 一种天然铀（VI）硅酸盐，二氧化铀实验室腐蚀研究的中间产物。铀硅酸盐空间群为P21 / m，铀五角双锥和硅氧四面体，其中的U（VI）与两氧组成UO22+分子，五个氧原子形成一个五边形，阳离子四面体的硅氧和OH -形成SiO3OH。形成的K /钠boltwoodite固溶。为了简化计算， boltwoodite夹层是由水和K组成（图1）。 计算方法： 3.1计算方法： CASTEP，平面波基组，周期性边界条件，超软赝势描述心电子与价电子的相互作用，GGA（PBE）用来近似电子交换和相关性。此外，自旋极化近似考虑两个未成对的5F电子 （Np5+）电子，Cutoff为800， K-point为0.07 ? - 1。总能量收敛性为1.0.10 -5 eV /原子，几何收敛为2.10 -5eV /原子。 3.2能量定义 定量结果（焓，融合能）必须加以比较，结合实验结果来验证计算准确性。实验上，焓或形成能是固体测量获得的。但是，量子力学计算获得的能量不是形成能，根据不同的计算方法，是晶格能（使用经验力场）或零价气相元素的反应能（如CASTEP计算）。图2流程图给出了三种化学反应的描述，从铀的形成（一）自然零价，（二）气体零价元素，（三）离子，分别对应（一）形成能，（二）CASTEP能源，以及（三）晶格能。 形成能由UO2和反应物（U和气态氧）能量差计算获得，其中气态氧气的氧分子为计算 在10 × 10 × 10一个盒子。 对铀的晶格能也可以计算，主要是在CASTEP中加入电离能（U0变成U4+的）和电子亲和力（O0变成O2-）。对电离能和电子亲和的值可采用高斯程序中量子力学计算或从实验测量。 Elattice = ECASTEP(UO2) – [(ECASTEP(U0) + IE)+2*(ECASTEP(O0) + EA)] 融合能的计算 融合能为NP取代U的反应过程的能量的变化。 K2(UO2)2(SiO3OH)2(H2O)3 + ?Np2O5 + ?O2 ? K2(UO2)(NpO2)(SiO3OH)2(H2O)3 + UO3 源是NP2O5，水复合物是UO3，生成物与反应物之间的能量差即为融合能，计算可以采用CASTEP计算获得的能量，也可以是形成能或是晶格能，只要每种化合物采用的标准一致即可。采用形成能或是晶格能同时提供了源与sink的相的信息。 融合能由依赖于融入机制，如H+增加或是层间取代，取代离子的位置以及取代反应发生时离子的相状态。比较了不同电荷平衡机制下的融入机制及离子取代位置；比较了H+融入机制下氧化物及水复合物的区别，比较不同层间的替代机制下氧化物和硅酸盐相的情况。 融合能高低的判断？应该是低的好吧？ 3.3 评估NP的融入极限 实验是测量不到NP融入铀酰的能量的，实验上能测量到的是NP的融入量，如果能计算出NP的融入量对实验有很大的帮助。 如何从理论上来评估融入极限没太明白？？从混合自由能曲线的最小值来判断。（11页） 4结果 4.1 取代机制1：Np5+ + H+ ? U6+ 计算了四中H的取代位置：1）硅氧四面体与铀多面体间的O上；2）硅氧四面体顶部的O上；3）层间水分子形成水合氢；4）镎酰的氧上。（最稳定位置） K2(UO2)2(SiO3OH)2(H2O)3 + ?Np2O5 + ?H2O ?HK2(UO2)(NpO2)(SiO3OH)2(H2O)3 + UO3 ΔEincorporation = 0.79 e