| 查看: 1169 | 回复: 5 | ||
cyclone985银虫 (小有名气)
|
[求助]
DFT 编程问题 已有1人参与
|
|
最近做了一些密度泛函编程方面的尝试，为了验证程序的正确性计算He原子基态的能量，使用高斯基矢线性组合（CGF）的基矢，使用解析方法计算动能势能和库伦积分。使用数值方法计算交换泛函的积分，发现结果总是和网站上给出的不一样，ftp://ftp.dl.ac.uk/qcg/dft_library/index.html#reference_data，大概有5%的误差，采用不同的计算方法结果都不正确。后来注意到使用的氦原子的DZVP基矢是非归一的。基矢的轨道指数和组合系数取自https://bse.pnl.gov/bse/portal，难道基矢有问题吗？此外其它的基矢如6-31G，6-31G**的第一个基矢也是非归一的，而3-21G的都是归一的。 请问非归一的基矢对计算是否有影响。这些基矢的数据是否不准确。 |
回复此楼
2楼2015-11-26 10:49:44
cyclone985
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 157.8
- 红花: 1
- 帖子: 175
- 在线: 34.9小时
- 虫号: 516204
- 注册: 2008-03-02
- 专业: 凝聚态物性 II ：电子结构
3楼2015-11-28 21:44:52
15588511668
新虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 324.6
- 散金: 700
- 红花: 7
- 沙发: 8
- 帖子: 1775
- 在线: 24.2小时
- 虫号: 4250160
- 注册: 2015-11-27
- 性别: MM
- 专业: 化学工程及工业化学
【答案】应助回帖
★
liliangfang: 金币+1, 谢谢交流 2015-11-29 08:07:24
liliangfang: 金币+1, 谢谢交流 2015-11-29 08:07:24
|
机组是正规称呼吧！不皈依的话。你至少要把，即函数怎么给求出来。最终的积分要乘以这个魔。 发自小木虫Android客户端 |
4楼2015-11-28 22:09:27
5楼2015-11-28 22:15:47
匿名
用户注销 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 269.9
- 散金: 103
- 帖子: 136
- 在线: 30.4小时
- 虫号: 0
- 注册: 2018-11-05
- 专业: 凝聚态物性 II ：电子结构
6楼2019-04-17 00:00:39
20