计算化学公社

 找回密码 Forget password
 注册 Register

谈谈约束性DFT (CDFT)

查看数: 44331 | 评论数: 23 | 收藏 Add to favorites 25
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2014-12-29 10:17

正文摘要:

注:本文只是简要介绍NWChem里CDFT的用法。在北京科音高级量子化学培训班(http://www.keinsci.com/workshop/KAQC_content.html)里有对CDFT原理的详细得多的介绍,并给出了NWChem做CDFT的更丰富的例子和明显更深入 ...

回复 Reply

zjxitcc 发表于 Post on 2023-5-6 21:21:04
sinap 发表于 2015-1-4 21:26
**** 作者被禁止或删除 内容自动屏蔽 ****

挖个坟。可以了,CP2K有限制性DFT,不仅有能量,连力都有(不过CDFT的力很少会被用到)。
granvia 发表于 Post on 2019-8-10 09:36:12
香港城市大学的Patrick Sit刚发表了控制氧化态的CDFT方法:https://doi.org/10.1021/acs.jctc.9b00281
yjr 发表于 Post on 2017-3-16 09:58:14
sobereva 发表于 2017-3-15 12:59
CDFT只能用于HF和DFT

谢谢您的及时解答!
sobereva 发表于 Post on 2017-3-15 12:59:28
yjr 发表于 2017-3-15 10:44
能问一下限制电荷片段的这个MP2优化怎么写?

CDFT只能用于HF和DFT
yjr 发表于 Post on 2017-3-15 10:44:38
能问一下限制电荷片段的这个MP2优化怎么写?
sobereva 发表于 Post on 2016-10-13 18:28:46
wswrpd 发表于 2016-10-13 11:36
有没有什么办法转化成multiwfn能分析的格式的方法?


unfortunately,据我所知没有办法,除非额外再写代码。要怪就怪nwchem不支持输出.molden、.wfn这些标准波函数文件格式
wswrpd 发表于 Post on 2016-10-13 11:36:43

有没有什么办法转化成multiwfn能分析的格式的方法?
sobereva 发表于 Post on 2016-10-13 04:33:26
wswrpd 发表于 2016-10-13 04:11
想偷懒问个问题,Multiwfn目前支持NWCHEM的波函数分析吗? 由于在下实在是懒人请不要嫌弃

不支持
wswrpd 发表于 Post on 2016-10-13 04:11:30
sobereva 发表于 2016-10-12 23:00
这只是个抽象的说法,可能指波函数什么的,具体问作者

想偷懒问个问题,Multiwfn目前支持NWCHEM的波函数分析吗? 由于在下实在是懒人请不要嫌弃
sobereva 发表于 Post on 2016-10-12 23:00:00
wswrpd 发表于 2016-10-12 17:04
叔丁基碳正离子 的例子 在NWCHEM-6.6下计算CDFT时,要在geometry 后面加 noautosym 变成 geometry noautosy ...

这只是个抽象的说法,可能指波函数什么的,具体问作者
wswrpd 发表于 Post on 2016-10-12 17:04:30
叔丁基碳正离子 的例子 在NWCHEM-6.6下计算CDFT时,要在geometry 后面加 noautosym 变成 geometry noautosym
否则会出现vectors were symmetry contaminated的问题。想请教sob大大,这个vectors are symmetry contaminated是什么意思?是什么vectors啊?
sobereva 发表于 Post on 2015-1-26 13:35:52
nkallwar 发表于 2015-1-26 12:48
选了一个小的结构尝试了一下,BS和UKS几乎是一样的

不知道作者是出于什么样的考虑


他也就是为了严谨些吧,毕竟BS设定下心里有谱,而直接UKS说不准收敛到对应于哪种情况。一般来说就直接UKS就行了,没必要这么复杂。

你说得对。

nkallwar 发表于 Post on 2015-1-26 12:48:24
sobereva 发表于 2015-1-17 15:16
1 也没什么影响,毕竟收敛后的波函数是相同的(或者说,即便确实有差异,也小到察觉不到)
2 如果确实差 ...


选了一个小的结构尝试了一下,BS和UKS几乎是一样的

不知道作者是出于什么样的考虑

paper里后面还有一段是这么写的:

UKS calculations for the S = 1 dication 6 underwent spontaneous symmetry breaking to yield a BS(3,1) solution.

Similarly, UKS calculations for S = 1/2 [Cr(PDI)2]1+ undergo spontaneous symmetry breaking to yield BS(3,2) solutions.

UKS calculations for neutral complex 2 (S = 1) converged exclusively to BS(3,1) solutions.

总之就是用UKS算出来和BS是一样的,于是认为他的BS片段划分是合理的, 可是这样的话,那他费这么大劲弄出来的片段划分到底有什么意义呢?

我弄的那个双核, 划成S=1,BS(4,2)(aaaa, bb) 与 S=1 BS(0,2) (abab,aa),以及直接UKS,S=1,结果完全是一样的,应该还是最后从布居分析中看实际更符合哪种“BS” 情况吧?

sobereva 发表于 Post on 2015-1-17 15:16:57
nkallwar 发表于 2015-1-17 12:43
嗯,从能量和几何结构参数上看,几乎就是一样的优化结果。

那么这带来两个问题: 1,文章中的BS策略 ...

1 也没什么影响,毕竟收敛后的波函数是相同的(或者说,即便确实有差异,也小到察觉不到)
2 如果确实差异大,那么谁的能量最低就用哪个。如果UKS明摆着能量比BS的高很多,作者却反倒用UKS,那只能问作者是否有什么特殊的考虑。

可以试试,但我估计和UKS的是一样的。如果能量极为相近,但电子结构差异却很大,那就太巧了,除非是两个状态由于什么原因能量是简并的。

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-13 10:34 , Processed in 0.178201 second(s), 26 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list