计算化学公社

标题: 计算富勒烯球心的NICS的步骤是否存在问题 [打印本页]

作者
Author: rorrim 时间: 2026-1-25 14:14
标题: 计算富勒烯球心的NICS的步骤是否存在问题
笔者想算（闭合的笼）富勒烯球心的NICS值，但不确定操作的方法是否正确，以下是操作步骤：
1.根据http://sobereva.com/643中开笼富勒烯的例子，用multiwfn计算球心，依次输入（以C60为例）cav，1-60，1，1，得到球心坐标
X/Y/Z of initial geometry center are -0.000000 -0.000000 -0.000000 Angstrom
Initial sphere radius is 1.849864 Angstrom

Step 1
Current coordinate: -0.000000 -0.000000 -0.000000 Angstrom
Gradient:       0.000059 0.000022 -0.000004  Norm 0.000063
Displacement: 0.000016 0.000006 -0.000001  Norm 0.000017 Angstrom
Goal: displacement norm <  0.01000000 Angstrom

Converged after    1 iterations

Final X/Y/Z of sphere center: 0.000016 0.000006 -0.000001 Angstrom
Radius is 1.849864 Angstrom
Diameter is 3.699728 Angstrom
Volume is    26.516010 Angstrom^3

2.在gjf文件中添加Bq 0.000016 0.000006 -0.000001，关键词为#p nmr=giao b3lyp/6-31g(d)
3.输出文件中Isotropic的负值就是NICS值
请问上述步骤是否存在问题？

作者
Author: sobereva 时间: 2026-1-25 16:12
通常就是在笼的几何中心计算，你这么搞多此一举。把结构文件载入Multiwfn，然后输入
100
21
[输入富勒烯上的原子序号；如果整个体系就是富勒烯就直接按回车]

马上屏幕上就得到了几何中心，即Geometry center (X/Y/Z)后面的

然后作为Bq原子写进输入文件就完了。=giao是多余的，本来这就是默认的

算完后取isotropic的负值

下文介绍NICS部分里也给了算NICS的例子
衡量芳香性的方法以及在Multiwfn中的计算
http://sobereva.com/176

作者
Author: rorrim 时间: 2026-1-25 16:22

sobereva 发表于 2026-1-25 16:12
通常就是在笼的几何中心计算，你这么搞多此一举。把结构文件载入Multiwfn，然后输入
100
21

非常感谢您的指点

欢迎光临计算化学公社 (http://bbs.keinsci.com/)