计算化学公社

标题: 求助,抗氧剂分子晶体失去四个H自由基后的THz光谱计算 [打印本页]

作者
Author:
xjtuglr    时间: 2025-5-19 16:48
标题: 求助,抗氧剂分子晶体失去四个H自由基后的THz光谱计算
本帖最后由 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?
谢谢大家
(, 下载次数 Times of downloads: 4)
抗氧剂1010
(, 下载次数 Times of downloads: 0) 抗氧剂1010cif
(, 下载次数 Times of downloads: 0) 去掉H原子以后的cif

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

作者
Author:
sobereva    时间: 2025-5-20 06:18
1 取决于是否要考虑H失去导致的结构弛豫效应,这取决于自己的研究目的

2 搜论坛上之前关于CP2K算磁性体系讨论的帖子,我回答过无数遍。如果参加过北京科音CP2K第一性原理计算培训班(http://www.keinsci.com/KFP),复习里面的“磁性体系的计算”部分,该讲的都讲得极为透彻了。当前体系的磁性和反铁磁性状态都可以算,恰当设置带单电子的氧的MAGNETIZATION(例如都设+1对应铁磁性状态,一半+1一半-1对应反铁磁性状态)以及这个体系的自旋多重度就完了
作者
Author:
xjtuglr    时间: 2025-5-20 11:31
sobereva 发表于 2025-5-20 06:18
1 取决于是否要考虑H失去导致的结构弛豫效应,这取决于自己的研究目的

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

好的谢谢老师




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3