24小时热门版块排行榜    

查看: 1280  |  回复: 6

ldzh0531

木虫 (小有名气)

[求助] 过渡态优化

请教用TS QST2 QST3优化过渡态时的输入格式怎么写啊,我看着好像要2行title。另外,我用赝势基组做TS优化怎么老是出错,构型读入就错误了。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beyondjun

铁杆木虫 (著名写手)

【答案】应助回帖

ldzh0531(金币+2): 多解释解释 qst2、qst3他们的关键词不会用啊 2011-11-01 19:21:16
引用回帖:
1楼: Originally posted by ldzh0531 at 2011-10-31 12:35:00:
请教用TS QST2 QST3优化过渡态时的输入格式怎么写啊,我看着好像要2行title。另外,我用赝势基组做TS优化怎么老是出错,构型读入就错误了。

建议贴出输入文件前几行命令行,好让其他愿意帮助你的人好诊断啊。我也在做过渡态优化,还没结束也倒还没出错

%d2e=/....../3tsopt.d2e
%int=/....../3tsopt.int
%rwf=/....../3tsopt.rwf
%scr=/....../3tsopt.scr
%mem=1800mb
%nprocshared=2  
%Chk=3tsopt.chk
#p opt=(calcfc,ts) ub3lyp/6-31g(d,p) scrf=(cpcm,solvent=cyclohexane) guess(,mix,always)

QST2 和QST3还不会,等待高手……
2楼2011-10-31 18:32:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

scottfreedom

金虫 (正式写手)

【答案】应助回帖

ldzh0531(金币+1): 主要是qst2、qst3他们的关键词不会用啊 2011-11-01 19:22:12
引用回帖:
1楼: Originally posted by ldzh0531 at 2011-10-31 12:35:00:
请教用TS QST2 QST3优化过渡态时的输入格式怎么写啊,我看着好像要2行title。另外,我用赝势基组做TS优化怎么老是出错,构型读入就错误了。

建议楼主可以在论坛里搜搜高斯教程,这样更有针对性!
Terse
3楼2011-10-31 22:37:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dreamyeye

木虫 (著名写手)

引用回帖:
1楼: Originally posted by ldzh0531 at 2011-10-31 12:35:00:
请教用TS QST2 QST3优化过渡态时的输入格式怎么写啊,我看着好像要2行title。另外,我用赝势基组做TS优化怎么老是出错,构型读入就错误了。

高斯自己应该有相应的例子
4楼2011-11-01 00:09:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ldzh0531

木虫 (小有名气)

引用回帖:
2楼: Originally posted by beyondjun at 2011-10-31 18:32:17:
建议贴出输入文件前几行命令行,好让其他愿意帮助你的人好诊断啊。我也在做过渡态优化,还没结束也倒还没出错

%d2e=/....../3tsopt.d2e
%int=/....../3tsopt.int
%rwf=/....../3tsopt.rwf
%scr=/....../ ...

%d2e=/....../3tsopt.d2e


%scr=/....../3tsopt.scr

%nprocshared=2  

scrf=(cpcm,solvent=cyclohexane) guess(,mix,always)
这些都是什么意思啊?
我用的是pc,好像用你的方法做不了啊
5楼2011-11-01 19:20:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beyondjun

铁杆木虫 (著名写手)

【答案】应助回帖

ldzh0531(金币+2): 谢谢,我也试过qst2 不过刚运行电脑就显示需要3个结构,所以我总是用qst3了,ts也试过,不过运行了一天又死了,不敢用了 2011-11-02 12:56:28
引用回帖:
5楼: Originally posted by ldzh0531 at 2011-11-01 19:20:03:
%d2e=/....../3tsopt.d2e


%scr=/....../3tsopt.scr

%nprocshared=2  

scrf=(cpcm,solvent=cyclohexane) guess(,mix,always)
这些都是什么意思啊?
我用的是pc,好像用你的方法做不了啊

哦,那前面那堆%就别管了
#p opt=(calcfc,ts) ub3lyp/6-31g(d,p) scrf=(cpcm,solvent=cyclohexane)  guess(,mix,always)
这句命令行中scrf=(cpcm,solvent=cyclohexane) 你可以不用,这是加入溶剂效应的命令。所以你的命令行应该写作


#p opt=(calcfc,ts) ub3lyp/6-31g(d,p) guess(,mix,always)

title
电荷 多重度
(猜测的过渡态原子坐标)


当然了,基组和方法你按照自己的需要稍作调整。#p opt=(calcfc,ts)是必须的。
QTS2据说不大靠谱,过渡态不好算,还是猜一个过渡态用TS好,在瞎猜方面,电脑不如人脑,呵呵……
QTS2格式如下:
# Opt=QST2 ub3lyp/6-31g Test

title
0,2
structure for reactant
0,2
structure for product
6楼2011-11-02 12:29:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beyondjun

铁杆木虫 (著名写手)

引用回帖:
1楼: Originally posted by ldzh0531 at 2011-10-31 12:35:00:
请教用TS QST2 QST3优化过渡态时的输入格式怎么写啊,我看着好像要2行title。另外,我用赝势基组做TS优化怎么老是出错,构型读入就错误了。

算过渡态不死才怪呢。建议你把高斯出错的提示看懂再对症下药,一般是9999错误,计算不能收敛。得加些命令,避免高斯钻牛角尖,你应该赶快贴出你的命令行和输出文件的出错提示,否则谁知道你的问题在哪里
7楼2011-11-02 14:16:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ldzh0531 的主题更新
信息提示
请填处理意见