|
|
本帖最后由 YZH_jesse 于 2022-11-18 22:17 编辑
各位老师：
我在用VASP计算一个带有季铵根的有机分子和碘分子 (I2) 的结合能，在优化过程中碘—碘键长逐渐拉长至330 pm，超过文献中报道的c.a. 270 pm (直观的感觉是I2分子散架了)
以下是我的初始结构和优化后的结构 (我在计算过程中观察到不正常的碘—碘键长后就取消了后续的计算，已运行255个离子步)
图1，初始结构，碘—碘键长为266 pm。
图2，优化后结构，碘—碘键长约330 pm。
以下是我的INCAR文件：
- ISTART = 0 # do not read wavecar information
- ICHARG = 2 # generate electronic information by initial guessing
- ENCUT = 450 # cutoff energy
- ISMEAR = 0
- SIGMA = 0.1
- ALGO = Fast # mixture of Davidson and RMM-DllS algorithms
- EDIFF = 1E-4 # the global break condition for the electronic SC-loop
- GGA = PE # pseudopotential PAW-GGA-PBE
- ISIF = 0 # cell shape: no; cell volume: no; positions: yes
- EDIFFG = -0.02 # the relaxation is stopped when the norms of all the forces are smaller than |EDIFFG|
- NSW = 500 # maximum number of ionic steps (maximum optimization steps)
- IBRION = 2 # ionic relaxation (conjugate gradient algorithm)
- IVDW = 11 # Dispersion correction DFT-D3
- LWAVE = .FALSE.
- LCHARG = .FALSE.
- LREAL = Auto
- NCORE = 16 # vasp strongly urge to set /# of cores per compute node
PS: 1. 当我将I2分子替换为I3-和I5-时，后两个分子会出现更加离谱的结构崩溃。
2. 已尝试过考虑自旋极化 (ISPIN = 2)；尝试过取消色散校正 (去除DFT-D3，虽然这样做完全没有道理)；将模拟的BOX尺寸增大以尽可能减小周期性影响；更换碘的赝势文件 (I_GW, I_sv_GW, 前者无影响，后者计算耗时太高)；尝试过引入隐式溶剂模型来做结构优化 (LSOL = .TRUE. )；尝试过用MD率先进行1 ns的动力学计算以尽可能优化初始结构；但所有的这些尝试都不能解决碘及碘的复合物结构崩溃的问题。
想问一下各位老师，碘及碘的复合物分子/离子的结构崩溃是什么原因引起的，我该如何解决这个问题呢？
|
|
发表于 Post on 2022-11-18 21:53:57
|
只看该作者 Only view this author
|只看大图 View big picture
|倒序浏览 Reverse view
|阅读模式 Reading model
收藏 Add to favorites
楼主 Author