b****e 发帖数: 906 | 1 具体问题在附件中
看起来那个投影算符似乎符合CCSD所使用的未定算符形式,并且这里更简单因为那个系数在我这里是预先确定的,如果能够简单修改一下CCSD的代码就能严格计算小系统的能量期望值就太好了 | b****e 发帖数: 906 | 2 有一点觉得不行了,似乎CCSD并没有计算中间过程的波函数的期望值,而是直接滑向了
最小值,望谁能给我一个肯定得答复,谢谢
系数在我这里是预先确定的,如果能够简单修改一下CCSD的代码就能严格计算小系统的
能量期望值就太好了
【在 b****e 的大作中提到】 : 具体问题在附件中 : 看起来那个投影算符似乎符合CCSD所使用的未定算符形式,并且这里更简单因为那个系数在我这里是预先确定的,如果能够简单修改一下CCSD的代码就能严格计算小系统的能量期望值就太好了
| b****e 发帖数: 906 | | 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到收敛时达到系统的基态。所以基态能量 : 也许是严格的,任何中间值都只是过度的。这个结论请各位大大加以评论一下,谢谢
|
|