求助，抗氧剂分子晶体失去四个H自由基后的THz光谱计算

xjtuglr · 发表于 Post on 2025-5-19 16:48:57

本帖最后由 xjtuglr 于 2025-5-19 19:31 编辑

请问各位老师，我想计算抗氧剂1010分子晶体失去4个H原子后酚羟基变成酚氧自由基之后的THz光谱，之前从COD官网下载了抗氧剂1010分子的cif文件并计算了THz光谱，现在遇到了问题求助一下大家
1.如果要计算失去H自由基后的结构，是用去掉H原子之前的结构完成opt之后去掉对应的H原子，还是直接从下载的cif文件去掉H原子再重新进行opt计算？
2.下载的cif文件我通过multiwfn去掉了H原子之后，该如何设置指定氧原子变为氧自由基？氧自由基在做opt计算的时候要需要添加什么语句吗，如何通过multiwfn设置MAGNETIZATION并开UKS？
谢谢大家

抗氧剂1010

Anti1010.cif (10.16 KB, 下载次数 Times of downloads: 0) 抗氧剂1010cif

noH.cif (9.92 KB, 下载次数 Times of downloads: 0) 去掉H原子以后的cif

1010noH.inp (47.46 KB, 下载次数 Times of downloads: 1) 修改以后的inp文件，麻烦大家看看这样可以了吗？没有改ESP_SCF的值，想先跑出一个初步结果再变严格

sobereva · 发表于 Post on 2025-5-20 06:18:41

1 取决于是否要考虑H失去导致的结构弛豫效应，这取决于自己的研究目的

2 搜论坛上之前关于CP2K算磁性体系讨论的帖子，我回答过无数遍。如果参加过北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/KFP），复习里面的“磁性体系的计算”部分，该讲的都讲得极为透彻了。当前体系的磁性和反铁磁性状态都可以算，恰当设置带单电子的氧的MAGNETIZATION（例如都设+1对应铁磁性状态，一半+1一半-1对应反铁磁性状态）以及这个体系的自旋多重度就完了

xjtuglr · 发表于 Post on 2025-5-20 11:31:39

sobereva 发表于 2025-5-20 06:18
1 取决于是否要考虑H失去导致的结构弛豫效应，这取决于自己的研究目的

2 搜论坛上之前关于CP2K算磁性体 ...

好的谢谢老师

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[CP2K] 求助，抗氧剂分子晶体失去四个H自由基后的THz光谱计算

浏览过的版块