| 查看: 1754 | 回复: 4 | ||||
584931986木虫 (小有名气)
|
[求助]
Siesta基组表达式是什么样子的？
|
|
我们知道在高斯程序包中，基组基本表示形式是高斯函数R（a, r）=r^(n-1)*exp(-ar^2),这是由于描述原子轨道的slater函数不易进行多中心分子积分，故而采用的高斯函数拟合slater函数的方法，达到描述原子轨道的目的。 但是在siesta中，我们只说是用数值原子轨道作为基组，进行求解，那所谓的数值原子轨道是什么样子的呢，有没有解析解呢？是角动量L和m的数值形式么？又是怎么得到的呢？ 而实际操作中，我们常采用SZP或者DZP作为基组，这其中的Zeta函数代表的是什么？它有没有基本表达式？而且zeta可以类比于高斯计算中的slater函数（表征原子轨道的波函数）中的变量吗？ 另外，在siesta计算中，基组采用了截断半径Rc，也就是说基组表示的波函数中R>Rc时，波函数为零。那么Rc表征的波函数是谁呢？是什么样的表达形式呢？是不是上边所说的L和m表征的数值原子轨道的表达形式呢？ 总之，就是不明白siesta的基组表示是什么样子的，真的是迷糊了，请大家帮忙解答。 [ Last edited by 584931986 on 2013-11-19 at 12:11 ] |
回复此楼» 猜你喜欢
让我中一个面上吧！
已经有15人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有5人回复
-
跳槽后在研项目怎么办？
已经有10人回复
-
什么时候开奖？
已经有10人回复
-
2026国自然函评费到账
已经有16人回复
-
今天放榜吗？
已经有16人回复
-
93BebMhtakh前后11位开头都是大写
已经有8人回复
-
范进中举一文的中心思想
已经有4人回复
-
只有每年这种时候来逛逛小木虫
已经有26人回复
-
filecode，4个jtjc了
已经有16人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教如何加快siesta的收敛速度
已经有9人回复
- Siesta从官网下载赝势计算出错
已经有7人回复
- 在gaussian计算中，先用小基组计算，在改成大基组计算出现的问题
已经有23人回复
- 求助siesta收敛问题
已经有3人回复
- 分子中电荷变化种种
已经有82人回复
- 请大家帮忙看看siesta输入文件有哪些错误
已经有8人回复
- SIESTA计算能级不准怎么解决
已经有10人回复
- 紧急求助：Gaussian输入问题，两个原子同时使用不同的赝势基组
已经有5人回复
- siesta计算dipole moment（电偶极矩）
已经有7人回复
- 请问siesta计算BornCharge时为什么输出的张量不满足对称性要求呢？谢谢
已经有15人回复
- 求助赝势基组的问题
已经有9人回复
- 为什么只有前四周期有全电子基组呢？
已经有3人回复
- 关于SIESTA中如何计算功函数的一种做法探讨
已经有8人回复
- 【分享】SIESTA中Bader analysis
已经有43人回复
- 【求助】SIESTA 中ATOM生成赝势测试后如何看赝势的好坏
已经有4人回复
- 【求助】如何在siesta中设置ghost state原子的计算
已经有16人回复
- SIESTA scf不收敛，怎么办？
已经有5人回复
- 【求助】siesta 计算DOS应该注意什么问题
已经有3人回复
- 【求助】请教在ATK中怎样使用自己的或者Siesta的赝势
已经有4人回复
- 【求助】siesta输入文件中赝势文件的作用？
已经有8人回复
- 【求助】siesta运行需要.psf文件吗
已经有4人回复
- 【求助】siesta中transiesta模块的收敛标准是啥~~
已经有20人回复
- 【求助】siesta中关于几何结构优化的问题【完结】
已经有17人回复
- 【讨论】关于赝势基组
已经有6人回复
- 【求助】求gamess 输入使用赝势基组与分裂基组的例子
已经有6人回复
wiseguy2008
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 752.5
- 帖子: 25
- 在线: 199.1小时
- 虫号: 1556693
- 注册: 2011-12-30
- 专业: 凝聚态物性 II ：电子结构
2楼2014-03-05 10:58:37
卡开发发
专家顾问 (著名写手)
Ab Initio Amateur
-
专家经验: +224 - 1ST强帖: 2
- 应助: 47 (小学生)
- 金币: 3566.8
- 散金: 14
- 红花: 223
- 帖子: 2477
- 在线: 1254小时
- 虫号: 1369423
- 注册: 2011-08-16
- 性别: GG
- 专业: 金属材料的磨损与磨蚀
- 管辖: 第一性原理
★ ★ ★
franch: 金币+3, 谢谢回帖交流。。。 2014-05-12 22:21:33
franch: 金币+3, 谢谢回帖交流。。。 2014-05-12 22:21:33
|
数值轨道意味着径向部分那些类似于gaussian的函数(GTO)都是格点的、数值的，而不是像gaussian函数那样解析的。这样的好处就是可以直接通过样条插值等直接进行数值积分，速度快得多。 这些zeta和和slater的应该也是类似的，就是所谓的“Slater指数”，即“有效电荷”，但是对于数值计算而言，尤其是DZ型以上的基组这个概念就模糊了。 Rc指的是指定原子轨道径向上的截断，原子轨道要通过基函数展开，这个Rcut自然也就是基函数的截断，对于siesta就是数值的离散的那些基函数。 具体更细节的部分我也并不了解，也许上面的有些地方还需前辈多指教。 因为Dmol3也是数值原子轨道的方法，这个问题我是参考了Dmol3的help文档。 |
3楼2014-03-05 23:13:24
超级月亮
木虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 2057.7
- 红花: 7
- 帖子: 1080
- 在线: 152.3小时
- 虫号: 2818437
- 注册: 2013-11-21
- 性别: GG
- 专业: 半导体光电子器件
4楼2014-05-07 15:21:13
卡开发发
专家顾问 (著名写手)
Ab Initio Amateur
-
专家经验: +224 - 1ST强帖: 2
- 应助: 47 (小学生)
- 金币: 3566.8
- 散金: 14
- 红花: 223
- 帖子: 2477
- 在线: 1254小时
- 虫号: 1369423
- 注册: 2011-08-16
- 性别: GG
- 专业: 金属材料的磨损与磨蚀
- 管辖: 第一性原理
5楼2014-05-07 15:46:55
