24小时热门版块排行榜    

查看: 611  |  回复: 4
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

kent1022

金虫 (著名写手)

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

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

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

过渡态计算

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注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写得很清楚,仔细看看吧。
华北电力大学(北京)数理学院,计算化学,凝聚态物理。欢迎合作!招硕士/博士研究生,招老师。
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的回帖
信息提示
请填处理意见