24小时热门版块排行榜    

查看: 899  |  回复: 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,我:8O.551.O54,科目全,可十急 +3 40ms4Wlo4umh 2026-09-14 3/150 2026-09-15 10:02 by VMDqgnlKimtZ
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:26 by C79jjtjAKjEn
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:14 by C79jjtjAKjEn
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 07:58 by C79jjtjAKjEn
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-15 07:46 by C79jjtjAKjEn
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-15 07:34 by C79jjtjAKjEn
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 06:50 by 5BDX0d0WFp7t
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-15 06:14 by 5BDX0d0WFp7t
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 3n8v2C8RimXI 2026-09-13 4/200 2026-09-15 06:01 by 5BDX0d0WFp7t
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:50 by 5BDX0d0WFp7t
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:38 by 5BDX0d0WFp7t
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:38 by 5BDX0d0WFp7t
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 QUjhNVAcOSff 2026-09-13 7/350 2026-09-15 05:22 by 5BDX0d0WFp7t
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 QUjhNVAcOSff 2026-09-13 8/400 2026-09-15 05:01 by mibUvS8DDCwf
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-15 01:21 by mibUvS8DDCwf
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-14 22:45 by vZfe6xYu34yj
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-14 22:25 by vZfe6xYu34yj
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +4 两三岁ss 2026-09-14 4/200 2026-09-14 20:29 by 班主任的同桌
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
[教师之家] 浙江师范大学是怎么坑我的 +10 recruit123 2026-09-09 15/750 2026-09-13 10:45 by seaskyy
信息提示
请填处理意见