24小时热门版块排行榜    

查看: 1662  |  回复: 8
当前主题已经存档。

nyzhaoyin

新虫 (正式写手)

[交流] 【分享】消除虚频的技巧

坛子里也有(搜索zzgyb的帖子)但还是拿上来再说一说!

如果输出文件有一个或两个虚频,数值大概在-40以内,最有效的方法就是把产生虚频的振动坐标直接与原坐标相加 然后继续优化

昨天晚上好几个结构都是有虚频,试了其他几种方法还是不行,就让坐标相加,优化之后计算频率  已经没有虚频了

所以感觉坐标相加这个方法很不错  我觉得首选就应该选这个  这个不行再尝试其他


回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chemniu

金虫 (文坛精英)

沙发,学习一下!
2楼2009-03-21 08:37:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyr3365754

木虫 (著名写手)

我有问题:
我每次把虚频坐标加到原子坐标上去以后,老是报错!
还有如果有两个虚频,那么虚频坐标就有两组,那么加到原子坐标上去的时候,是不是把两组坐标 XX 和 YY 相加后再加到原子坐标的 X 和 Y 上去!

谢谢楼主解答!
3楼2009-03-21 08:43:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nyzhaoyin

新虫 (正式写手)

引用回帖:
Originally posted by zyr3365754 at 2009-3-21 08:43:
我有问题:
我每次把虚频坐标加到原子坐标上去以后,老是报错!
还有如果有两个虚频,那么虚频坐标就有两组,那么加到原子坐标上去的时候,是不是把两组坐标 XX 和 YY 相加后再加到原子坐标的 X 和 Y 上去!

...

是的  有两组了就把两组都加上 然后优化

ps:你出错  错误是什么?
4楼2009-03-21 08:46:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

半支烟5591

银虫 (正式写手)

专家


引用回帖:
Originally posted by nyzhaoyin at 2009-3-21 08:46:

是的  有两组了就把两组都加上 然后优化

ps:你出错  错误是什么?

不必把两组都加上,只要选最负的虚频对应的坐标加上那个就可以!
5楼2009-03-21 10:11:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nyzhaoyin

新虫 (正式写手)

引用回帖:
Originally posted by 半支烟5591 at 2009-3-21 10:11:

不必把两组都加上,只要选最负的虚频对应的坐标加上那个就可以!

哦  学习了  不过两个加上也没错
6楼2009-03-21 10:51:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chenxf122

银虫 (正式写手)

先加那个最负的,然后计算,结果还有虚频的话就再加另一个,如果加了后报错的话就在计算路径中加入 iop(5/13=1)就可以了
7楼2009-03-23 09:10:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyr3365754

木虫 (著名写手)

这是我的输入文件,加了虚频坐标后的,运行后就报错。麻烦给看看!谢谢!
%chk=44.chk
%mem=6MW
%nproc=1
# opt=tight freq ub3lyp/lanl2dz scrf=(solvent=water,check) guess=read
geom=(connectivity,allcheck) genchk

Au2Pd2

0 1
Pd                 0.96188900   -0.14443900    0.00028900
Pd                -4.24098800    0.45255900    0.00001800
Au                 3.58099600    0.15298000   -0.00007900
Au                -1.67500200   -0.32599900   -0.00010000

1 3 1.0 4 1.0
2 4 1.0
3
4
8楼2009-03-23 16:58:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nyzhaoyin

新虫 (正式写手)

引用回帖:
Originally posted by zyr3365754 at 2009-3-23 16:58:
这是我的输入文件,加了虚频坐标后的,运行后就报错。麻烦给看看!谢谢!
%chk=44.chk
%mem=6MW
%nproc=1
# opt=tight freq ub3lyp/lanl2dz scrf=(solvent=water,check) guess=read
geom=(connectivity,allc ...

把那些没用的关键词全去掉 条件太多了  你就不能先简单的优化 优化完不成了再加其他条件吗?  一开始限制条件太多了。。。汗

还有 内存能改大点吗?  你机子内存多大啊?
9楼2009-03-24 13:33:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 nyzhaoyin 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 cBGhfJFuNMca 2026-08-02 19/950 2026-08-02 16:24 by jHcrHMhbKCFV
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 15/750 2026-08-02 16:23 by jHcrHMhbKCFV
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 19/950 2026-08-02 16:23 by jHcrHMhbKCFV
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 19/950 2026-08-02 16:21 by jHcrHMhbKCFV
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +18 低垂的野花 2026-07-31 24/1200 2026-08-02 16:03 by swallowwh
[基金申请] 面上提前没消息,有中的吗 +9 archvillain 2026-08-02 11/550 2026-08-02 15:44 by alexon
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 8/400 2026-08-02 15:01 by jHcrHMhbKCFV
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 cBGhfJFuNMca 2026-08-02 14/700 2026-08-02 14:58 by jHcrHMhbKCFV
[找工作] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 12/600 2026-08-02 13:43 by jHcrHMhbKCFV
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cBGhfJFuNMca 2026-08-02 14/700 2026-08-02 13:34 by jHcrHMhbKCFV
[基金申请] 2026年国自然面上资助率 +15 布布和一二 2026-07-30 21/1050 2026-08-02 10:56 by iaeyuan
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 10/500 2026-08-02 10:00 by 锦衣卫寒战
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[基金申请] 微信指数没变化,科研之友没阅读 +14 wangze12014 2026-07-28 18/900 2026-07-31 18:35 by lzchaihuo
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
[基金申请] 这种情况还有戏吗 +5 drbart 2026-07-27 11/550 2026-07-29 05:50 by drbart
[基金申请] 同事接到电话了,我却没有 +4 1234567wang 2026-07-27 4/200 2026-07-28 10:11 by GOODLUCKER
[基金申请] 准备明年的基金了 +5 Tide man 2026-07-27 5/250 2026-07-28 09:25 by newfuzzy1
信息提示
请填处理意见