24小时热门版块排行榜    

Znn3bq.jpeg
查看: 741  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[有机交流] 如何实现卤原子转化 10+3 BT20230424 2026-05-15 5/250 2026-05-16 16:20 by czyzsu
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 cjf4bx70cj 2026-05-14 6/300 2026-05-16 16:17 by 0i5p09z61n
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k37jurhrau 2026-05-16 3/150 2026-05-16 13:57 by vcdazktkjx
[论文投稿] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 xx7gd5zq4e 2026-05-15 5/250 2026-05-16 12:24 by h3oerqvkv9
[硕博家园] 申请博士 +3 呃?呃 2026-05-15 3/150 2026-05-16 11:01 by a4742549
[文学芳草园] 裁员滚滚,退居二线 +4 J_wei 2026-05-10 4/200 2026-05-16 10:52 by zh10246
[考博] 2026博士还有哪些学校有名额 +5 小王求读研 2026-05-15 6/300 2026-05-16 10:44 by a4742549
[基金申请] 这年头没有找到涵评专家,还有中面上的可能吗 +11 dd921ww 2026-05-12 13/650 2026-05-16 09:16 by Howard28
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 l7k6xnh0yc 2026-05-14 4/200 2026-05-16 04:44 by k37jurhrau
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 l7k6xnh0yc 2026-05-14 5/250 2026-05-16 04:29 by k37jurhrau
[找工作] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 l7k6xnh0yc 2026-05-14 3/150 2026-05-16 01:49 by k37jurhrau
[考博] 西南大学考核制博士 +4 lijunjie84 2026-05-11 7/350 2026-05-15 23:20 by 同仁堂教主
[论文投稿] 有带发论文的吗 +4 山楂之术 2026-05-09 4/200 2026-05-15 15:40 by 妹子不好惹
[文学芳草园] 风把牡丹吹跑了 +5 myrtle 2026-05-12 9/450 2026-05-15 15:27 by myrtle
[基金申请] 精华III评审感受-评审感受-评审感受 +14 ferrarichen 2026-05-11 18/900 2026-05-15 11:12 by cmhchen
[基金申请] 重磅!青年科学基金项目(C类)资助增幅预计超过50% +5 水和泥不是水泥 2026-05-13 7/350 2026-05-14 20:57 by 水和泥不是水泥
[有机交流] 求助2,4-二氯-5-嘧啶甲醛的合成方法 20+3 光吃不拉 2026-05-14 5/250 2026-05-14 20:15 by 一切都是空工
[高分子] 本人最近太闲了,谁有问题可以提,每天会统一回复 +8 一切都是空工 2026-05-12 19/950 2026-05-14 20:03 by 一切都是空工
[考博] 材料类只有一篇综述能申博么 +4 乐逍遥谷 2026-05-13 4/200 2026-05-14 12:05 by zhyzzh
[考博] 现在不知道怎么办,感觉很痛苦 +4 qweww 2026-05-11 5/250 2026-05-11 20:23 by Oversize
信息提示
请填处理意见