24小时热门版块排行榜    

查看: 345  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 准备明年的基金了 +7 Tide man 2026-07-27 7/350 2026-08-02 16:43 by gltch
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cBGhfJFuNMca 2026-08-02 15/750 2026-08-02 16:31 by jHcrHMhbKCFV
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 13/650 2026-08-02 16:27 by jHcrHMhbKCFV
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 cBGhfJFuNMca 2026-08-02 19/950 2026-08-02 16:24 by jHcrHMhbKCFV
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 19/950 2026-08-02 16:23 by jHcrHMhbKCFV
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 15/750 2026-08-02 16:23 by jHcrHMhbKCFV
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 19/950 2026-08-02 16:21 by jHcrHMhbKCFV
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +18 低垂的野花 2026-07-31 24/1200 2026-08-02 16:03 by swallowwh
[基金申请] 面上提前没消息,有中的吗 +9 archvillain 2026-08-02 11/550 2026-08-02 15:44 by alexon
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 cBGhfJFuNMca 2026-08-02 16/800 2026-08-02 13:34 by jHcrHMhbKCFV
[基金申请] 2026年国自然面上资助率 +15 布布和一二 2026-07-30 21/1050 2026-08-02 10:56 by iaeyuan
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 10/500 2026-08-02 10:00 by 锦衣卫寒战
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[有机交流] 一个有机合成实验室都需要哪些设备? 50+3 kf2781974 2026-07-31 7/350 2026-08-02 07:59 by 88817753
[基金申请] 微信指数没变化,科研之友没阅读 +14 wangze12014 2026-07-28 18/900 2026-07-31 18:35 by lzchaihuo
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
[有机交流] 产物和副产物价值比较 30+4 小汤02 2026-07-27 5/250 2026-07-30 12:44 by czyzsu
[基金申请] 这种情况还有戏吗 +5 drbart 2026-07-27 11/550 2026-07-29 05:50 by drbart
[基金申请] 同事接到电话了,我却没有 +4 1234567wang 2026-07-27 4/200 2026-07-28 10:11 by GOODLUCKER
信息提示
请填处理意见