应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » scan同时能否进行布局分析
51/1返回列表
查看: 918  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

kent1022

金虫 (著名写手)

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

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

» 收录本帖的淘帖专辑推荐

过渡态计算

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-03-28 09:55:09
已阅   回复此楼   关注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的回帖
查看全部 5 个回答

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的回帖

yanrding

金虫 (正式写手)

【答案】应助回帖

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

第二种方法我们在做modr软扫描时用过，所以scan肯定也可以，不过需要仔细琢磨那一段iop。sob写得很清楚，仔细看看吧。
赞一下(1人) 回复此楼
华北电力大学（北京）数理学院，计算化学，凝聚态物理。欢迎合作！招硕士/博士研究生，招老师。
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的回帖
查看全部 5 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[有机交流] OpenAI又攻克数学难题了，离AI做有机合成还有多远？ +8 wangyl_123 2026-10-10 12/600 2026-10-13 20:40 by Merryya
[硕博家园] 极限拉扯，如何抉择 +9 Isak0320 2026-10-12 13/650 2026-10-13 11:57 by Isak0320
[教师之家] 课题组会帮你通过考核么？ +4 好运来～9527 2026-10-10 7/350 2026-10-13 08:41 by xli1984
[基金申请] 入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷！
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬，路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[有机交流] 我应该选哪条合成路线让实验成功率高一些？ +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
关闭
请填处理意见
关闭