24小时热门版块排行榜    

查看: 250  |  回复: 2
当前主题已经存档。

wuy069

铁杆木虫 (正式写手)

[交流] 【求助】如何任意选取格点大小

The grid size of the density cube may be chosen arbitrarily,
thus balancing the computational cost of the coupling calculation
with the errors introduced by the graining of space within the
density cube. We find the Gaussian 94 default grid size to be
quite conservative and the resulting coupling calculations
lengthy. For the work here, we employ a coarser graining in
which each cube element contains a volume of about 0.23 bohr^3,
which is roughly twice the volume of the Gaussian 94 default。
在Gaussian中格点选取:
命令行加cube=(potential,cards)
分子坐标后空行
a.cube               输出文件名称(可要可不要)
IFlag, X0, Y0, Z0 输出单元编号和网格原点。
N1, X1, Y1, Z1    X-方向的点数和点的间隔。
N2, X2, Y2, Z2    Y-方向的点数和点的间隔。
N3, X3, Y3, Z3    Z-方向的点数和点的间隔。
这样能任意选取格点吗?这篇文章思想是把跃迁密度cube分割成很多cube元,如何任意选取格点大小,怎样做到使cube元的体积大小为0.23bohr^3?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主

★ ★ ★ ★
wuy069(金币+4,VIP+0):谢谢 12-8 10:16
以前我也问过这个问题
无人回答,后来我自己解决了它
解决过程在http://www.mdbbs.org/thread-19342-1-1.html
2楼2009-12-05 10:39:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)

★ ★ ★ ★ ★ ★ ★ ★ ★
yjcmwgk(金币+5,VIP+0): 12-7 17:53
wuy069(金币+4,VIP+0):谢谢 12-8 10:17
X0, Y0, Z0是cube最左下方的起点,通常是三个相等的负数。

X1、Y2、Z3为格点的步长。如果你要求cube元的体积大小为0.23bohr^3,那你对X1、Y2、Z3都取值0.23即可。

N1、N2、N3是三个坐标方向的格点数目。显然N*0.23就是空间一个轴的长度,它是你cube空间的大小。cube空间的大小要能装下你的分子。

N*0.23就是空间一个轴的长度,也决定了X0, Y0, Z0的取值。通常把坐标原点取在图的中心。

还要注意,如果IFlag取正值才是bohr单位,如果取负值,则单位采用埃。

[ Last edited by zhou2009 on 2009-12-8 at 11:28 ]
3楼2009-12-07 15:38:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wuy069 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见