物理化学学报 >> 2020, Vol. 36 >> Issue (11): 1908035.doi: 10.3866/PKU.WHXB201908035

论文 上一篇    下一篇

镎(Ⅳ)、镅(Ⅲ)、锔(Ⅲ)的AMBER力场参数化及评估

刘子义1,2, 夏苗仁1,2, 柴之芳1,3, 王东琪1,*()   

  1. 1 中国科学院高能物理研究所多学科中心, 北京 100049
    2 中国科学院大学化学学院, 北京 100049
    3 苏州大学放射医学与辐射防护国家重点实验室, 苏州大学放射医学与防护学院, 江苏 苏州 215123
  • 收稿日期:2019-08-29 录用日期:2019-09-24 发布日期:2019-09-29
  • 通讯作者: 王东琪 E-mail:dwang@ihep.ac.cn
  • 基金资助:
    国家自然科学基金(91026000);国家自然科学基金(21473206);国家自然科学基金(91226105);中国科学院百人计划项目(Y2291810S3)

Parameterization and Validation of AMBER Force Field for Np4+, Am3+, and Cm3+

Ziyi Liu1,2, Miaoren Xia1,2, Zhifang Chai1,3, Dongqi Wang1,*()   

  1. 1 Multidisciplinary Initiative Center, Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 100049, P. R. China
    2 University of Chinese Academy of Sciences, Beijing 100049, P. R. China
    3 State Key Laboratory of Radiation Medicine and Protection, and School of Radiation Medicine and Interdisciplinary Sciences (RAD-X), Soochow University, Suzhou 215123, Jiangsu Province, P. R. China
  • Received:2019-08-29 Accepted:2019-09-24 Published:2019-09-29
  • Contact: Dongqi Wang E-mail:dwang@ihep.ac.cn
  • Supported by:
    the National Natural Science Foundation of China(91026000);the National Natural Science Foundation of China(21473206);the National Natural Science Foundation of China(91226105);CAS Hundred Talents Program(Y2291810S3)

摘要:

:锕系元素具有高放射性和毒性,为开展锕系水溶液化学的实验研究带来了困难与挑战。得益于理论计算方法和计算能力的快速发展,运用分子动力学模拟水溶液环境并研究锕系在水溶液当中的化学行为成为一种可替代的方法。本文拟合了Np4+、Am3+、Cm3+三种金属离子的力场参数并进行了评估,进而将其用于溶液配位化学和动力学研究。为了研究的系统性,也开展了Th4+、U4+、Pu4+三种锕系离子的分子模拟研究。基于对上述六种离子的分子动力学研究,系统报道了该六种锕系离子的第一、二溶剂化层结构、性质、立体化学及驻留时间等方面的共性与特性,并从能量方面探究这些内在差异的原因。结果显示,本文报道的参数适用于锕系离子的配位结构和结合自由能的研究,但对于动力学性质的研究,宜谨慎用于相同价态的锕系离子溶液体系的定性理解。本文也探究了六种锕系离子在水溶液中与Cl-、NO3-和CO32-的相互作用,分析了各配合物配位模式、成键规律及结合强弱。综上所述,本项工作发展了三种关键锕系离子的力场参数,有助于开展基于AMBER力场的动力学方法的An3+/4+溶液化学的研究,丰富了对An3+/4+的水溶液中配位化学的认识。

关键词: An3+, An4+, AMBER力场

Abstract:

The radioactivity and toxicity of actinides impede experimental investigation into their chemical properties in the condensed phase. The rapid development of computational methods and computational facilities allows for alternative experimental methods, including the use of a molecular force field, to gain insight into the coordination chemistry and dynamics of actinides. The key to this method is the force fieild parameters. In the present work, we report the development and validation of the AMBER (Assisted Model Building with Energy Refinement) force field parameters for Np4+, Am3+, and Cm3+ based on the experimentally determined ion-oxygen distance (IOD). The parameter set, together with that reported for Th4+, U4+, and Pu4+, was then applied to investigate the coordination chemistry and dynamics of these six actinide ions in the aqueous phase, in the absence and presence of counterions Cl-, NO3-, and CO32-. The simulations showed a shorter An-Ow coordination length for An4+ than for An3+, and for higher atomic numbers of ions in the same valence state. Th4+ preferentially existed in a 10-coordinated state, adopting a BCASP (bicapped square antiprism) conformation, while the other ions tended to be 9-coordinated with a CASP (capped square antiprism) conformation. The only exception was Cm3+, which adopted a TCTP (tricapped trigonal prism) conformation. The results also showed that the water molecules around An4+ were more ordered than those around An3+, as indicated by the smaller angles between the An-Ow vector and the dipole direction of the water ligand. This highly ordered structure of coordinated water affected their translation and rotation, i.e., the diffusion coefficient and rotational relaxation time of the water molecules around An4+ were smaller than those in the case of An3+ due to the stronger electrostatic interaction between An4+ and ligating water. The hydration free energies of the targeted actinide ions were also calculated by the FEP (free energy perturbation) method. An4+ underwent a greater degree of stabilization than did An3+ upon hydration; among the ions in the same oxidation states, those with a higher atomic number were better stabilized. In summary, the results of the simulations were consistent with the literature data in terms of the hydration structure, coordination of counterions, and hydration free energy of the actinide ions. The ability of the parameter set to describe the dynamics of water in the vicinity of actinides remains to be verified due to the lack of reference data. We tentatively propose that it may be used to investigate the coordination chemistry of actinides both in conformational analysis and binding strength, while special care should be taken when studying the kinetics of the solvated system. This work is expected to enrich our understanding of the solution behavior of An3+/An4+ at the force field level.

Key words: An3+, An4+, AMBER force field