24小时热门版块排行榜    

查看: 1478  |  回复: 3

buxiangshuo

铜虫 (初入文坛)

[求助] 同样用DFT模拟,castep或Dmol与gaussian有什么差别?

我会gaussian,但不懂castep或Dmol,同样是DFT模拟的情况下,不明白后者为何要设置截止能和边界条件,因为高斯并不用设置;根本原因是什么呢,是因为它们采用的近似算法不同吗,还是基组的不同或其他原因导致的?谢谢!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

★ ★ ★ ★ ★
zhou2009: 2016-07-04 10:42:44
zhou2009: 金币+5 2016-07-04 10:43:10
castep是做周期体系的平面波+赝势程序。周期体系需要用边界条件,gaussian的pbc也要做这种设置。平面波难于描述芯轨道,截断能越低需要的平面波就越多,计算量也就越大,但结果未必更好。

dmol和gaussian、castep有相似的地方,但是用赝势+数值基组
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
2楼2016-07-04 06:16:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

buxiangshuo

铜虫 (初入文坛)

引用回帖:
2楼: Originally posted by beefly at 2016-07-04 06:16:13
castep是做周期体系的平面波+赝势程序。周期体系需要用边界条件,gaussian的pbc也要做这种设置。平面波难于描述芯轨道,截断能越低需要的平面波就越多,计算量也就越大,但结果未必更好。

dmol和gaussian、caste ...

谢谢大神!这么说来,主要是因为使用基组函数的不同,导致它们这些参数设置和计算特点上有差别吗?因为我的理解是castep用平面波/赝势函数如同gaussian计算小分子体系时使用斯莱特或高斯基组函数,我的理解对吗?
3楼2016-07-04 09:04:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

引入周期性之后,根据Bloch定理会引入一个Bloch vector,所以实际周期性的计算要多设置一个k网格,dmol可以不设置周期边界,这种情况下也不必(程序也不允许)设置k-point。

大致从图像上说,如果将孤立的分子在有限距离上平移复制无限个,就能够得到周期体系;反过来说,把周期体系的边界拉开到无穷大就会回到孤立体系,清楚这一点我们就能够知道大概怎么进行计算(比如怎么用castep算孤立的分子)。

dmol和gaussian采用的基函数比较接近,同属于LCAO形式的,控制基函数的尺寸主要是看原子轨道的zeta数目、扩展基组等(因为劈裂价基组还涉及到GTO的数目),计算几乎在实空间完成;castep则是采用的平面波,控制基函数尺寸是通过能量截断的大小来完成的,计算中动能项和赝势的一部分的积分是在倒空间完成(某些的势能项可能在实空间完成)。同样都是周期性的话,这两种不同的形式存在变换关系。

平面波方法是肯定需要边界条件的,因为在无限的空间中基函数没办法归一化。还有一点,波函数在原子核附近有岐点效应(可以认为是振荡强烈),这将会导致需要非常高频(高能)的平面波才能够比较好的描述,所以才考虑赝势或者缀加原子轨道。

平面波的方法量子化学的书上说的少,一般固体理论方法之类或者第一性原理的书,如R M Martin的书其中会提到,谢希德先生的《固体能带理论》里面也有。
不一定挂在论坛,计算问题问题欢迎留言。
4楼2016-07-05 18:28:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 buxiangshuo 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 329求调剂 +3 钮恩雪 2026-03-25 3/150 2026-03-25 14:43 by 糖加冰
[考研] 293求调剂 +7 加一一九 2026-03-24 7/350 2026-03-25 12:02 by userper
[考研] 化学调剂 +6 yzysaa 2026-03-21 6/300 2026-03-25 09:27 by aa331100
[考研] 085601求调剂总分293英一数二 +3 钢铁大炮 2026-03-24 3/150 2026-03-24 22:03 by bingxueer79
[考研] 08工学调剂 +12 用户573181 2026-03-20 17/850 2026-03-24 21:45 by 番茄炖了西红柿
[考研] 【考研调剂】化学专业 281分,一志愿四川大学,诚心求调剂 +15 吃吃吃才有意义 2026-03-19 16/800 2026-03-24 19:51 by 了了了了。。
[考研] 一志愿北化315 求调剂 +3 akrrain 2026-03-24 3/150 2026-03-24 19:35 by 了了了了。。
[考研] 一志愿吉大化学322求调剂 +4 17501029541 2026-03-23 6/300 2026-03-24 10:21 by 戴围脖的小蚊子
[考研] 求调剂 +7 十三加油 2026-03-21 7/350 2026-03-23 23:48 by 热情沙漠
[考研] 350求调剂 +6 weudhdk 2026-03-19 6/300 2026-03-23 15:47 by tangyuan0840221
[考研] 280分求调剂 一志愿085802 +4 PUMPT 2026-03-22 7/350 2026-03-22 22:13 by 星空星月
[考研] 寻找调剂 +4 倔强芒? 2026-03-21 4/200 2026-03-22 16:14 by 木托莫露露
[考研] 311求调剂 +3 26研0 2026-03-20 3/150 2026-03-22 14:46 by ColorlessPI
[考研] 初试 317 +7 半拉月丙 2026-03-20 7/350 2026-03-21 22:26 by peike
[考研] 一志愿南大,0703化学,分数336,求调剂 +3 收到VS 2026-03-21 3/150 2026-03-21 18:42 by 学员8dgXkO
[考研] 材料与化工(0856)304求 B区 调剂 +3 邱gl 2026-03-21 3/150 2026-03-21 13:47 by lature00
[考研] 330求调剂0854 +3 assdll 2026-03-21 3/150 2026-03-21 13:01 by 搏击518
[考研] 求调剂 +6 Mqqqqqq 2026-03-19 6/300 2026-03-21 08:04 by JourneyLucky
[考研] 一志愿中海洋材料工程专硕330分求调剂 +8 小材化本科 2026-03-18 8/400 2026-03-20 23:16 by JourneyLucky
[考研] 求调剂 +3 @taotao 2026-03-20 3/150 2026-03-20 19:35 by JourneyLucky
信息提示
请填处理意见