由买买提看人间百态

boards

本页内容为未名空间相应帖子的节选和存档,一周内的贴子最多显示50字,超过一周显示500字 访问原贴
Chemistry版 - 能否稍微修改一下CCSD的代码来计算我的这个问题?
相关主题
问一个gaussian计算的问题这里有没有Phi Lambda Upsilon Honorary Chemical Society (转载)
量化疑问求助文獻一篇 (转载)
一个关于DFT的初级问题:求助文獻 (转载)
问个electronic coupling的问题吧paper help,Thanks.----- PCCP
文献求助 Rev. Mod. Phys. 78, 973 (2006) [17 pages]paper help
paper needed, help.求chem phys lett 文献两篇,谢谢!!!
paper help, thanks文献求助
求一篇Ann. Rev. Phys. Chem(已求到,谢谢kenny514paper help please (links inside) thank you!
相关话题的讨论汇总
话题: ccsd话题: 计算话题: 代码话题: 修改话题: 能量
进入Chemistry版参与讨论
1 (共1页)
b****e
发帖数: 906
1
具体问题在附件中
看起来那个投影算符似乎符合CCSD所使用的未定算符形式,并且这里更简单因为那个系数在我这里是预先确定的,如果能够简单修改一下CCSD的代码就能严格计算小系统的能量期望值就太好了
b****e
发帖数: 906
2
有一点觉得不行了,似乎CCSD并没有计算中间过程的波函数的期望值,而是直接滑向了
最小值,望谁能给我一个肯定得答复,谢谢

系数在我这里是预先确定的,如果能够简单修改一下CCSD的代码就能严格计算小系统的
能量期望值就太好了

【在 b****e 的大作中提到】
: 具体问题在附件中
: 看起来那个投影算符似乎符合CCSD所使用的未定算符形式,并且这里更简单因为那个系数在我这里是预先确定的,如果能够简单修改一下CCSD的代码就能严格计算小系统的能量期望值就太好了

b****e
发帖数: 906
3
en, please help ...
a****o
发帖数: 6612
4
其实你没有必要解释Hamiltonian, 或者分子轨道, 或者第一和第二量子化。能看懂你
问题的肯定有这些背景。但是你对那个投影算符没有解释清楚。好奇一下,你的公式(1
) 求和符号下面为什么没有j ? 如果你用CCSD, 你是基于原子轨道积分(alpha, beta)
还是基于分子轨道积分(比如i, j,说实在的,我不懂你图片里i, j的含义)。 我的
理解是对于多电子体系,能量算法涉及到四个电子轨道 (i, j, k, l).
如果真的对应CCSD的话,你可以尝试用的你参数(平方?可能有phase factor的问题需
要解决。)替换amplitudes t_ir, t_ijrs (i, j是占据轨道,r, s是空轨道)。
商业化的软件里CCSD的解不是"variational". Variational CC 的理论有,但是计算公
式展开有无穷多项,而现在通用的方法是有限多项。
即使你能够代入计算,也得不到你倒过来第二个公式的能量。还有,我估计你用参数算
出来的能量,误差太大,没有实际用途。
b****e
发帖数: 906
5
谢谢aaddo的回复。没想到仔细查过的note还是有这么多的typo,不好意思,公式(1)
中的j是个typo,应该为i,同时它的spin也错了,应该为down。之所以写这么多是因为
询问一个化学系的大牛他要求这么多的解释,然后就没有了下文,所以希望到这里能找
到答案。现在我倾向于觉得CCSD没法计算一个给定算符的期望值,CCSD本身并不严格计
算一个算符的期望值,而是利用了某个trick到收敛时达到系统的基态。所以基态能量
也许是严格的,任何中间值都只是过度的。这个结论请各位大大加以评论一下,谢谢
a****o
发帖数: 6612
6
有人做过CC各项展开的计算,不过他们的程序并不在公开的程序里面。见
“Alternative single-reference coupled cluster approaches for multireference
problems: The simpler, the better”, F. A. Evangelista, J. Chem. Phys. 134,
224102 (2011)
(link: http://jcp.aip.org/resource/1/jcpsa6/v134/i22/p224102_s1?isAuthorized=no)

【在 b****e 的大作中提到】
: 谢谢aaddo的回复。没想到仔细查过的note还是有这么多的typo,不好意思,公式(1)
: 中的j是个typo,应该为i,同时它的spin也错了,应该为down。之所以写这么多是因为
: 询问一个化学系的大牛他要求这么多的解释,然后就没有了下文,所以希望到这里能找
: 到答案。现在我倾向于觉得CCSD没法计算一个给定算符的期望值,CCSD本身并不严格计
: 算一个算符的期望值,而是利用了某个trick到收敛时达到系统的基态。所以基态能量
: 也许是严格的,任何中间值都只是过度的。这个结论请各位大大加以评论一下,谢谢

1 (共1页)
进入Chemistry版参与讨论
相关主题
paper help please (links inside) thank you!文献求助 Rev. Mod. Phys. 78, 973 (2006) [17 pages]
paper helppaper needed, help.
请问哪里能看到被J. Chem. Phys.接收了的文章?paper help, thanks
paper help, thanks!求一篇Ann. Rev. Phys. Chem(已求到,谢谢kenny514
问一个gaussian计算的问题这里有没有Phi Lambda Upsilon Honorary Chemical Society (转载)
量化疑问求助文獻一篇 (转载)
一个关于DFT的初级问题:求助文獻 (转载)
问个electronic coupling的问题吧paper help,Thanks.----- PCCP
相关话题的讨论汇总
话题: ccsd话题: 计算话题: 代码话题: 修改话题: 能量