24小时热门版块排行榜    

查看: 1671  |  回复: 10

yafang2012

木虫 (正式写手)

[求助] cas计算问题 已有1人参与

我想要用CAS算freq,算出来opt的结构了,直接读opt后的chk,为什么一开始算上就死了呢?log文件看不懂啊,,今天算freq总是出问题,opt可以,就是算不错freq。。郁闷死了。
本来是用的大基组6-31g(d),算不出来,就改用小基组算了。有没有大神支支招的。。多谢多谢~~
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : 1-freq.log
  • 2015-04-22 21:07:58, 19.53 K

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

经验

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

virtualzx

木虫 (著名写手)

★ ★ ★ ★
gmy1990: 金币+2, 谢谢交流 2015-04-24 15:49:43
gmy1990: 金币+2, 谢谢交流 2015-04-24 15:50:09
引用回帖:
5楼: Originally posted by yafang2012 at 2015-04-22 20:10:14
这个是opt的输出文件。轨道我不是很会选,我计算的是一个反应的中间体,因为之前有分子计算是选的这个活性空间,为了一致,所以只有这样选了。而且这个分子有硝基,好像有N原子的选轨道比较复杂,,我不太懂。。像 ...

看你的输出是需要11.4G内存。进行MCSCF计算需要很了解需要选取哪些轨道。如果不想考虑这一点的话,就不要用CAS了。硝基完全可以用SCF处理。如果涉及激发态、自由基、成断键、电子组态变化等才需要MCSCF,而这些相对应的轨道放在活性空间里。不要试图用增加MCSCF活性空间来处理电子的动态相关,这样很昂贵,而且不准确,有时甚至可以比单组态方法得到的结果更差。动态相关用DFT、微扰理论、CC或CI来处理。
6楼2015-04-23 09:23:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

virtualzx

木虫 (著名写手)

引用回帖:
5楼: Originally posted by yafang2012 at 2015-04-22 20:10:14
这个是opt的输出文件。轨道我不是很会选,我计算的是一个反应的中间体,因为之前有分子计算是选的这个活性空间,为了一致,所以只有这样选了。而且这个分子有硝基,好像有N原子的选轨道比较复杂,,我不太懂。。像 ...

如果对活性空间有疑问,建议看这篇文章。看完基本就是专家了。

http://www.annualreviews.org/doi ... v.physchem.49.1.233
7楼2015-04-23 09:39:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

virtualzx

木虫 (著名写手)


gmy1990: 金币+1, 谢谢交流 2015-04-24 15:49:12
(10,10)CAS很大啊。freq需要计算解析二阶导数,需要的内存比opt大得多。你真的需要这么大的活性空间吗?CAS不像DFT,不能随便设定活性空间,你需要很小心地想好想要那些轨道在活性空间里。6-31g(d)对于波函数相关方法来说是蛮小的基组。活性空间比基组对计算量的影响还要大得多。

[ 发自手机版 http://muchong.com/3g ]

» 本帖已获得的红花(最新10朵)

2楼2015-04-23 03:23:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

枪下游魂

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
yafang2012: 金币+30, ★★★★★最佳答案, 谢谢~ 2015-04-23 11:05:09
gmy1990: 金币+2, 谢谢交流 2015-04-24 15:49:30
你没有贴出结构优化的输出文件,所以我不知道你在结构优化上的设置,诚如楼上所说,这种规模的分子,真的需要casscf(10,10)这么大的活性空间么?另外你在结构优化时是否对活性空间内的轨道进行了选择呢?看最后问题是内存不够了。如果一定要用cas,可以先从最简单的4,4开始。

» 本帖已获得的红花(最新10朵)

3楼2015-04-23 08:16:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yafang2012

木虫 (正式写手)

送红花一朵
引用回帖:
2楼: Originally posted by virtualzx at 2015-04-23 03:23:41
(10,10)CAS很大啊。freq需要计算解析二阶导数,需要的内存比opt大得多。你真的需要这么大的活性空间吗?CAS不像DFT,不能随便设定活性空间,你需要很小心地想好想要那些轨道在活性空间里。6-31g(d)对于波函数相关方 ...

因为我算的是几个反应的中间体,反应过程中可能会涉及到很多轨道,老师的意思是尽可能多的把反应轨道包含进去。所以她不同意我改小活化空间。轨道其实我不是很会选,我看他之前算过一个分子选的的(10,10)。这个是因为内存不足所以死掉的吗?我菜鸟一个额,,我还有一个文件,我觉得算的很奇怪,freq好像算过了。出现4个yes,但是最后还是死了。 下面这个命令好像算了两次,我不太懂。你有时间帮我看下。多谢多谢~~
#N Geom=AllCheck Guess=TCheck SCRF=Check GenChk CASSCF/STO-3G Freq
4楼2015-04-23 09:03:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yafang2012

木虫 (正式写手)

送红花一朵
引用回帖:
3楼: Originally posted by 枪下游魂 at 2015-04-23 08:16:18
你没有贴出结构优化的输出文件,所以我不知道你在结构优化上的设置,诚如楼上所说,这种规模的分子,真的需要casscf(10,10)这么大的活性空间么?另外你在结构优化时是否对活性空间内的轨道进行了选择呢?看最后问 ...

这个是opt的输出文件。轨道我不是很会选,我计算的是一个反应的中间体,因为之前有分子计算是选的这个活性空间,为了一致,所以只有这样选了。而且这个分子有硝基,好像有N原子的选轨道比较复杂,,我不太懂。。像这样的分子在这个活性空间下算freq时候需要内存是多大?我能不能再改大一点内存,看可不可以?
5楼2015-04-23 09:10:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yafang2012

木虫 (正式写手)

引用回帖:
6楼: Originally posted by virtualzx at 2015-04-23 09:23:34
看你的输出是需要11.4G内存。进行MCSCF计算需要很了解需要选取哪些轨道。如果不想考虑这一点的话,就不要用CAS了。硝基完全可以用SCF处理。如果涉及激发态、自由基、成断键、电子组态变化等才需要MCSCF,而这些相对 ...

谢谢,大神啊,膜拜一下~~虽然有些地方我还是不太懂,嘿嘿。老师不同意改活性空间,我可能要试试其他的方法了,
8楼2015-04-23 10:55:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

virtualzx

木虫 (著名写手)

引用回帖:
8楼: Originally posted by yafang2012 at 2015-04-22 21:55:46
谢谢,大神啊,膜拜一下~~虽然有些地方我还是不太懂,嘿嘿。老师不同意改活性空间,我可能要试试其他的方法了,...

STO-3G都需要这么大,可能是有点困难。STO-3G太小,是没法发表的。如果一定要用这么大的活性空间,建议不要用高斯。molpro,molcas或者GAMESS会好一些。
9楼2015-04-23 12:49:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yafang2012

木虫 (正式写手)

引用回帖:
9楼: Originally posted by virtualzx at 2015-04-23 12:49:20
STO-3G都需要这么大,可能是有点困难。STO-3G太小,是没法发表的。如果一定要用这么大的活性空间,建议不要用高斯。molpro,molcas或者GAMESS会好一些。...

恩恩,我准备用molcas算呢,有没有molcas的入门资料?我什么都不懂呀,,,
10楼2015-04-23 13:22:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yafang2012 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jRl3mE6ddZGq 2026-08-03 8/400 2026-08-04 06:15 by k7OM8YghWbkC
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 cu9Nq1xK233Z 2026-08-03 6/300 2026-08-04 06:10 by k7OM8YghWbkC
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 cu9Nq1xK233Z 2026-08-03 9/450 2026-08-04 06:07 by k7OM8YghWbkC
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 jRl3mE6ddZGq 2026-08-03 4/200 2026-08-04 04:54 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 6/300 2026-08-04 04:47 by k7OM8YghWbkC
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 04:47 by k7OM8YghWbkC
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 7/350 2026-08-04 04:46 by k7OM8YghWbkC
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 04:46 by k7OM8YghWbkC
[有机交流] 一个有机合成实验室都需要哪些设备? 50+3 kf2781974 2026-07-31 9/450 2026-08-04 03:56 by kf2781974
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 03:21 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 cu9Nq1xK233Z 2026-08-03 12/600 2026-08-04 03:17 by k7OM8YghWbkC
[高分子] UV压敏胶开发 +3 ichall 2026-07-30 5/250 2026-08-03 14:40 by Sunrisepay
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 16/800 2026-08-03 12:23 by fuweiguochen
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +19 低垂的野花 2026-07-31 25/1250 2026-08-03 09:25 by gy116024
[基金申请] 2026年国自然面上资助率 +16 布布和一二 2026-07-30 22/1100 2026-08-03 09:20 by gy116024
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 12/600 2026-08-03 08:22 by syl200707
[基金申请] 微信指数没变化,科研之友没阅读 +15 wangze12014 2026-07-28 19/950 2026-08-02 20:05 by 蔡棒棒菂
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[高分子] 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
信息提示
请填处理意见