24小时热门版块排行榜    

查看: 1482  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 2JOx3r2CYEgw 2026-08-21 9/450 2026-08-23 18:46 by lXgU0RvNTDC9
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 h4CP7TrQR8Lg 2026-08-22 3/150 2026-08-23 16:30 by lXgU0RvNTDC9
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 14:57 by KM5EcsNQRBPn
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 6/300 2026-08-23 14:21 by KM5EcsNQRBPn
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +3 process2012 2026-08-23 4/200 2026-08-23 12:54 by huixian257
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 h4CP7TrQR8Lg 2026-08-22 5/250 2026-08-23 12:09 by LR9qGULyN2ew
[基金申请] 什么时候开奖? +9 CrisMessi 2026-08-18 10/500 2026-08-23 12:03 by 丶昵称占用
[基金申请] 2026国自然函评费到账 +15 羊腰板 2026-08-21 16/800 2026-08-23 10:45 by process2012
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 9/450 2026-08-23 10:44 by IXZuIJ2Q7OVy
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 10/500 2026-08-23 07:17 by OEbVnUOu01ol
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 03:52 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QTy3jDtz1uLt 2026-08-21 9/450 2026-08-23 03:04 by OEbVnUOu01ol
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 DpSrDtM079iu 2026-08-22 3/150 2026-08-23 00:55 by OEbVnUOu01ol
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 建议基金发布提前给出明确的时间点 +10 kulium 2026-08-21 13/650 2026-08-21 21:58 by alongwaytogo
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
信息提示
请填处理意见