24小时热门版块排行榜    

查看: 868  |  回复: 4

kent1022

金虫 (著名写手)

[求助] scan同时能否进行布局分析

对一体系进行势能曲线扫描,发现程序只计算了不同点的单点能。为了确定其离解极限,可不可以在scan的同时要求程序进行布局分析,即在每个不同核间距计算单点能是也计算电荷布局?如能,输入程序怎样写?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yanrding

金虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
zhou2009: 金币+3 2013-03-28 11:56:06
kent1022: 金币+10 2013-03-29 13:21:43
可以。
方法一是手动扫描,即计算许多个高斯任务,每次只算一个点的,当然可以写一个脚本做。
方法二是改IOP,具体可参考sob的文章。
Gaussian的Link、IOp与非标准计算路径 文/Sobereva 2010-MAR-9
原文地址: http://hi.baidu.com/sobereva/blo ... a71d14bf09e643.html
华北电力大学(北京)数理学院,计算化学,凝聚态物理。欢迎合作!招硕士/博士研究生,招老师。
2楼2013-03-28 11:03:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kent1022

金虫 (著名写手)

引用回帖:
2楼: Originally posted by yanrding at 2013-03-28 11:03:15
可以。
方法一是手动扫描,即计算许多个高斯任务,每次只算一个点的,当然可以写一个脚本做。
方法二是改IOP,具体可参考sob的文章。
Gaussian的Link、IOp与非标准计算路径 文/Sobereva 2010-MAR-9
原文地址:  ...

感谢回帖!
我计算的是一负分子离子体系,尝试过第一种方法,当然并不是对所有的点都进行了计算,而是间隔大概1~2埃时计算单点能。根据此时各点的电荷布局分析判断体系的离解极限。问题是这样分别计算的单点能其波函数均是各自分别猜测而得到,而用scan命令扫描时后一个点的初始波函是前一个点的计算所得满足收敛标准的波函,如此迭代。这两者应该是有所差别的,对结果的影响有多大?

对第二种方法,本人看得是云里雾里,不过还是稍有眉目。
scan及pop=full的计算路径输入文件分别为
----
scan
----
1/38=1,60=1/1,8;
2/12=2,17=6,18=5,29=3,40=1/2;
3/11=9,16=1,25=1,30=1/1,2,3;
4//1;
5/5=2,38=5/2;
6/7=2,8=2,9=2,10=2,18=1,28=1/1;
1/60=1/8(1);
99/9=1/99;
2/29=3/2;
3/11=9,16=1,25=1,30=1/1,2,3;
4/5=5,16=3/1;
5/5=2,38=5/2;
1/60=1/8(-4);
99/9=1/99;

--------
pop=full
--------
1/38=1/1;
2/12=2,17=6,18=5,40=1/2;
3/11=9,16=1,25=1,30=1/1,2,3;
4//1;
5/5=2,38=5/2;
6/7=3,28=1/1;
99/5=1,9=1/99;

怎样将两者结合起来,在scan的同时计算每个点的电荷布局,本人不知从何下手,请高手指点!
3楼2013-03-28 17:01:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yanrding

金虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
kent1022: 金币+10 2013-03-29 13:21:56
前一种方法,理论上如果都是收敛到正确的波函数,那么应该结果是一样的。
当然手工扫描的时候,也可以read前一次的chk啊,为什么非要“各自分别猜测”?

第二种方法我们在做modr软扫描时用过,所以scan肯定也可以,不过需要仔细琢磨那一段iop。sob写得很清楚,仔细看看吧。
华北电力大学(北京)数理学院,计算化学,凝聚态物理。欢迎合作!招硕士/博士研究生,招老师。
4楼2013-03-28 17:11:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kent1022

金虫 (著名写手)

引用回帖:
2楼: Originally posted by yanrding at 2013-03-28 11:03:15
可以。
方法一是手动扫描,即计算许多个高斯任务,每次只算一个点的,当然可以写一个脚本做。
方法二是改IOP,具体可参考sob的文章。
Gaussian的Link、IOp与非标准计算路径 文/Sobereva 2010-MAR-9
原文地址:  ...

感谢回帖!
我计算的是一负分子离子体系,尝试过第一种方法,当然并不是对所有的点都进行了计算,而是间隔大概1~2埃时计算单点能。根据此时各点的电荷布局分析判断体系的离解极限。问题是这样分别计算的单点能其波函数均是各自分别猜测而得到,而用scan命令扫描时后一个点的初始波函是前一个点的计算所得满足收敛标准的波函,如此迭代。这两者应该是有所差别的,对结果的影响有多大?

对第二种方法,本人看得是云里雾里,不过还是稍有眉目。
scan及pop=full的计算路径输入文件分别为
----
scan
----
1/38=1,60=1/1,8;
2/12=2,17=6,18=5,29=3,40=1/2;
3/11=9,16=1,25=1,30=1/1,2,3;
4//1;
5/5=2,38=5/2;
6/7=2,8=2,9=2,10=2,18=1,28=1/1;
1/60=1/8(1);
99/9=1/99;
2/29=3/2;
3/11=9,16=1,25=1,30=1/1,2,3;
4/5=5,16=3/1;
5/5=2,38=5/2;
1/60=1/8(-4);
99/9=1/99;

--------
pop=full
--------
1/38=1/1;
2/12=2,17=6,18=5,40=1/2;
3/11=9,16=1,25=1,30=1/1,2,3;
4//1;
5/5=2,38=5/2;
6/7=3,28=1/1;
99/5=1,9=1/99;

怎样将两者结合起来,在scan的同时计算每个点的电荷布局,本人不知从何下手,请高手指点!
5楼2013-03-28 22:16:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 kent1022 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 7K1CJE38xLG4 2026-08-25 4/200 2026-08-26 02:04 by cNXvBfCpiZOM
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 7K1CJE38xLG4 2026-08-25 4/200 2026-08-26 02:04 by cNXvBfCpiZOM
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +17 木水思豆 2026-08-25 20/1000 2026-08-26 01:26 by chengyan1220
[基金申请] 今天务委会开完了,明天出结果吗 +18 angus9576 2026-08-25 22/1100 2026-08-26 00:41 by merchancy
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 有没有大神帮我看看基金代码 28+4 1234567wang 2026-08-24 10/500 2026-08-25 19:15 by lfy8008
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 19/950 2026-08-25 17:57 by chick875
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +3 苏东坡二世 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +5 process2012 2026-08-23 7/350 2026-08-25 09:42 by huixian257
[基金申请] 只有每年这种时候来逛逛小木虫 +25 yaoyewhu2008 2026-08-20 27/1350 2026-08-25 08:01 by Equinoxhua
[基金申请] filecode,4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤?
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
信息提示
请填处理意见