24小时热门版块排行榜    

查看: 1284  |  回复: 20
当前主题已经存档。

xiaowandouer

木虫 (正式写手)

★ ★
csfn(金币+2,VIP+0):thanks,欢迎常来计算版~
看来能算多大的体系与计算机的配置、所用的计算水平等有很大关系
我算的体系重原子也就四五个,再加上一些其他的小分子
基组中等水平:对重原子lanl2dz,其他6-311+g*
因为有过渡金属,经常遇到l502难收敛,所以一般都加scf=qc,这个可能也是费时的原因,再有计算只是在双核2.4G酷瑞CPU上运算,没并行

和大家交流真好,学到了不少东西
11楼2007-12-25 15:44:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hairan

木虫 (著名写手)

瞌睡木虫

★ ★
csfn(金币+2,VIP+0):thanks,欢迎常来计算版~
引用回帖:
Originally posted by xiaowandouer at 2007-12-25 15:44:
看来能算多大的体系与计算机的配置、所用的计算水平等有很大关系
我算的体系重原子也就四五个,再加上一些其他的小分子
基组中等水平:对重原子lanl2dz,其他6-311+g*
因为有过渡金属,经常遇到l502难收敛,所以 ...

原来如此,直接用scf=qc有点太费时,你可以用这个关键词,一般可以取得较好的效果
scf=(maxcycle=100,xqc)
意思是,如果scf在100个循环之后还没有收敛的话,转换到qc计算
一般来说502不收敛都发生在80个循环之后,能收敛的话50个循环之前就收敛了,这样的关键词可以在计算时间和收敛之间取得较好的平衡。
你试试看,希望能对你有帮助
据权威报纸调查,84%的青少年犯罪者有玩电子游戏的经历,故应该禁止经营电子游戏。据我们所知,100%的青少年犯罪者有穿鞋的经历,故制鞋厂应该停业。
12楼2007-12-25 18:22:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaowandouer

木虫 (正式写手)

太感谢了!这将会节省我很多时间!
不过还有一点不清楚:scf=(maxcyc=100,xqc)中的x需要带入具体的数值吗?还是直接这样用就行?
13楼2007-12-25 20:48:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hairan

木虫 (著名写手)

瞌睡木虫

引用回帖:
Originally posted by xiaowandouer at 2007-12-25 20:48:
太感谢了!这将会节省我很多时间!
不过还有一点不清楚:scf=(maxcyc=100,xqc)中的x需要带入具体的数值吗?还是直接这样用就行?

直接用就行了,我也是从网上搜出来的解决办法,祝你顺利
据权威报纸调查,84%的青少年犯罪者有玩电子游戏的经历,故应该禁止经营电子游戏。据我们所知,100%的青少年犯罪者有穿鞋的经历,故制鞋厂应该停业。
14楼2007-12-25 22:19:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaowandouer

木虫 (正式写手)

多谢hairan的耐心指导!
也谢谢其他虫子们
15楼2007-12-26 09:30:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hairan

木虫 (著名写手)

瞌睡木虫

引用回帖:
Originally posted by xiaowandouer at 2007-12-25 15:44:
看来能算多大的体系与计算机的配置、所用的计算水平等有很大关系
我算的体系重原子也就四五个,再加上一些其他的小分子
基组中等水平:对重原子lanl2dz,其他6-311+g*
因为有过渡金属,经常遇到l502难收敛,所以 ...

另外,你算的重金属如果是第一过渡系的,也就是从Sc到Zn,不推荐用lanl2dz,有点太小了,我在计算中就经常发现这个基组的计算结果不太可靠,一般用它来搜索初始结构和粗优化,除非你能证明这个基组得到的结果和大基组的一样。否则审稿人可能会不认账。最差也要用到6-311+g*。第二第三过渡系的金属也有一些比lanl2dz大一些的基组,另外对Ag,Au,Cd,Hg,可以使用aug-cc-pvdz/aug-cc-pvtz,结果的可信度要高得多。

[ Last edited by hairan on 2007-12-26 at 10:48 ]
据权威报纸调查,84%的青少年犯罪者有玩电子游戏的经历,故应该禁止经营电子游戏。据我们所知,100%的青少年犯罪者有穿鞋的经历,故制鞋厂应该停业。
16楼2007-12-26 10:40:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaowandouer

木虫 (正式写手)


zzgyb(金币+1,VIP+0):谢谢您的参与,欢迎您常来计算模拟版解答问题。
我算的是第三过渡系的金属,文献上用lanl2dz挺多,再有,得出的结果还算可靠(与实验值相比,eg:IE,EA等),以后尝试下您说的基组

hairan对计算化学确实颇有研究,而且乐于奉献,向您致敬!
17楼2007-12-26 15:46:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lei0736

荣誉版主 (职业作家)

优秀版主


zzgyb(金币+1,VIP+0):谢谢您的参与,欢迎您常来计算模拟版解答问题。
文献报道的用EAM势计算的体系有几千万个原子 做位错 和马氏体相变的 呵呵
18楼2007-12-27 14:38:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaowandouer

木虫 (正式写手)

引用回帖:
Originally posted by lei0736 at 2007-12-27 14:38:
文献报道的用EAM势计算的体系有几千万个原子 做位错 和马氏体相变的 呵呵

长见识了,呵呵
19楼2007-12-28 09:03:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaowandouer

木虫 (正式写手)

to hairan:
     用您之前说的scf=(maxcyc=100,xqc),我把原先用scf=qc能正常优化结束的构型的输入文件又重新算了一遍,想看看到底能节省多少时间,结果发现构型变掉了。我想主要是因为收敛算法的不同导致的,亦或是初始猜测不够好,不知道是不是这样?
20楼2007-12-28 09:45:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiaowandouer 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见