24小时热门版块排行榜    

CyRhmU.jpeg
查看: 796  |  回复: 8
当前主题已经存档。
【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 sandy_too 的 5 个金币

sandy_too

木虫 (著名写手)

忽忽家族之忽忽猪

[交流] 【求助】计算溶剂效应onsager模型的输入命令

最近在计算溶剂效应时的时候用SCIPCM模型时一直都有虚频,而在气相时的计算结果都没有虚频,开始以为是结构错误,然后换了不同的对称性,不同的结构,还是算不出来。
后来老板说是模型错误,用onsager模型计算,结果还是死掉……
  我的输入命令是:B3LYP 6-311+G**  opt scrf=(dipole,sovent=THF) freq
结果是死掉,在l502就死掉了
在论坛上看到onsager模型还要输入a0(A),请问这个参数在哪里可以查到?
谢谢高手!
回复此楼
岁月就象一条河,左岸是无法忘却的回忆,右岸是值得把握的青春年华,中间飞快流淌的,是年轻隐隐的伤感。世间有许多美好的东西,但真正属于自己的却并不多。看
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

digilab

木虫 (小有名气)

★ ★
sandy_too(金币+2,VIP+0):谢谢你的解答,但是怎么计算volume呢?在高斯说明书里找不到…… 11-21 10:08
计算volume得到a0
2楼2009-11-21 09:52:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

alexander8568

木虫 (正式写手)

★ ★ ★ ★ ★
sandy_too(金币+5,VIP+0):谢谢解答,我再试试看! 11-21 10:55
频率和优化不能同时进行吧,你先优化然后才算吧!试试看
3楼2009-11-21 10:47:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjp01

木虫 (小有名气)

★ ★ ★
sandy_too(金币+3,VIP+0):谢谢解答,pcm模型和onsager模型都不适用,还有什么模型呢? 11-21 10:56
onsager模型可能不太适合你需要计算的这个物质,你可以考虑使用其他的模型进行计算。
4楼2009-11-21 10:48:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

digilab

木虫 (小有名气)


sandy_too(金币+1,VIP+0):谢谢了! 11-21 17:57
高斯说明书里有volume关键词

volume就跟freq opt是同类的啊
5楼2009-11-21 15:02:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dummy1

木虫 (正式写手)


yjcmwgk(金币+1,VIP+0): 11-22 21:09
你可以尝试一下别的PCM方法,里面很多或许有的可以,或者实在不行就调整一下构型,再就是调一下收敛的算法什么的。
引用回帖:
Originally posted by sandy_too at 2009-11-20 13:54:
最近在计算溶剂效应时的时候用SCIPCM模型时一直都有虚频,而在气相时的计算结果都没有虚频,开始以为是结构错误,然后换了不同的对称性,不同的结构,还是算不出来。
后来老板说是模型错误,用onsager模型计算, ...

6楼2009-11-22 09:27:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

swyswy

铜虫 (小有名气)


yjcmwgk(金币+1,VIP+0): 11-22 21:09
通过vol的多次计算得出a。,再进行溶剂化计算,其中需要输入a。和所用溶剂的介电常数
7楼2009-11-22 09:31:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

digilab

木虫 (小有名气)

★ ★ ★
sandy_too(金币+3,VIP+0):谢谢解答,我在试试看! 11-23 11:01
优化好结构以后用volume关键字计算
搜索Recommended a0 for SCRF calculation就可以了
8楼2009-11-23 09:34:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hzfish

金虫 (正式写手)


sandy_too(金币+1,VIP+0):谢谢,辛苦了! 11-24 09:38
GAUSSIAN自带的例子
10_01b.com
%Chk=cl2_hf
#T RHF/6-31+G(D) Test SCRF=Dipole

GAUCHE-1,2-DICHLOROETHANE

0,1
C
C,1,RCC
CL,1,RCL,2,A1
CL,2,RCL,1,A1,3,B1
H,1,RCH1,2,A2,4,B2
H,1,RCH2,2,A3,4,B3
H,2,RCH1,1,A2,3,B2
H,2,RCH2,1,A3,3,B3

RCC=1.51490
RCL=1.78594
RCH1=1.08019
RCH2=1.07802
A1=112.88944
A2=108.96431
A3=111.38813
B1=70.51710
B2=188.89555
B3=-50.42251

3.65 10.1

3.65=>A0
10.1=>epson
9楼2009-11-23 15:49:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 sandy_too 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见