24小时热门版块排行榜

返回列表

shuaizhaoz

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 2
帖子: 11
在线: 2.8小时
虫号: 27500565
注册: 2021-10-10
专业: 固体力学

[求助] Multiwfn画DFT的电子密度图时，格点（grid point）处的值是如何得到的？已有1人参与

新手使用multiwfn后，有一些疑问希望能够搞清楚，grid point（default: 200,200）是可调的参数，那么很有可能是通过插值的办法来得到，那么fchk文件中是否有电子密度的参考格点的值？
如果有，要如何在fchk中查看呢？
如果没有，是如何计算得到的呢？
我采用的泛函是B3LYP，基组是6-31g，发现DFT积分格点是一个可调的参数，例如int=ultrafine或superfine等，但是不太清楚具体空间的格点数量/密度，不像VASP可以通过NGX,NGY,NGZ直接指定，所以有些困惑，希望大神来解答一下。

回复此楼

» 猜你喜欢

书籍求助：汽车市场营销理论与实务（电子版）——章小平已经有0人回复
东方理-中科大联合博士生招聘已经有0人回复
物理化学论文润色/翻译怎么收费? 已经有229人回复
钴酸锂半电池小倍率容量上不去已经有1人回复
吉林大学材料物理本科生求问调剂信息已经有23人回复
2026第二届光电子与半导体器件前沿技术研讨会——光电子赋能·半导体创芯！已经有1人回复
光电子赋能·半导体创芯！已经有1人回复
光电子赋能·半导体创芯！已经有0人回复
Chemical Bonding at Surfaces and Interfaces,最经典的一本表面上化学相互作用教材已经有0人回复
QE计算电声耦合的时候报错Error in routine lambda (100)wrong or too many modes 已经有2人回复
分子动力学模拟合作需求已经有0人回复

1楼 2023-05-04 04:40:20

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

paramecium86

版主 (著名写手)

QC强帖: 5
应助: 587 (博士)
贵宾: 0.01
金币: 26356.1
红花: 348
帖子: 2024
在线: 1441.6小时
虫号: 2466443
注册: 2013-05-16
性别: GG
专业: 理论和计算化学
管辖: 量子化学

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
shuaizhaoz: 金币+26 2023-05-05 13:34:29

高斯的积分格点可以直接用数字来控制. 不过并不是通过笛卡尔的形式而是极坐标形式。也就是定义径向（radial）和角度（angular）的格点。用int 关键词来控制。  int(Grid=99302) 就代表了径向有99个而角度是302个格点。总格点数=99*302 个也就是差不多3万个。而默认的int=ultrafine 实际上是 int(Grid=99590).
高斯的有一些固定的格点级别可以用关键词来调用从低到高

Coarse（33,110）< SG1 (50,194)< Fine (99,302)<UltraFine (99,590)< SuperFine (250,974)，它对H和He原子是（175,974）

不想用上述的默认设置好的格点关键词，就用int(Grid=xxxxxx) 来自己设置。

除此之外，高斯在做freq的时候会用到CPHF，这里也需要格点。CPHF的格点在进行普通freq的时候会默认比Int的格点低两个级别比如默认 int=ultrafine 的时候对于freq的CPHF来说格点就是 SG1 级别的，某些特殊任务下的cphf会高一级达到fine 比如计算非谐振之类的。而设置CPHF格点使用  CPHF(Grid=xxxxxx)

赞一下(3人)