应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】过渡态问题 分子结构输入部分如何得来？
141/1返回列表
查看: 2021  |  回复: 27
查看全部回帖 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

zhaoyxcas

捐助贵宾 (正式写手)

[交流] 【求助】过渡态问题 分子结构输入部分如何得来？

%mem=64mb
#p RHF/6-31G(d) Opt=QST2

SiH2 + H2 --> SiH4 Reactants

0,1
Si
X 1 1.0
H 1 1.48 2 55.0
H 1 1.48 2 55.0 3 180.0
H 1 R 2 A1 3 90.0
H 1 R 5 A2 2 180.0

R=2.0
A1=80.0
A2=22.0

SiH2 + H2 --> SiH4 Products

0,1
Si
X 1 1.0
H 1 1.48 2 55.0
H 1 1.48 2 55.0 3 180.0
H 1 R 2 A1 3 90.0
H 1 R 5 A2 2 180.0

R=1.48
A1=125.2
A2=109.5
分子结构输入部分怎么得来的？ 不是原子坐标？？费解，谢谢！！
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2010-11-11 11:10:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhaoyxcas

捐助贵宾 (正式写手)
引用回帖:
Originally posted by zhangmt at 2010-11-11 13:49:43:
当然是坐标，只不过不是你通常看到的笛卡尔坐标，而是内坐标。也就是用键长、键角、二面角来描述的分子坐标。
由于是用QST2方式找过渡态的，所以分别给了反应物和产物两套坐标，让gaussian根据提供的情况去猜测、 ...

嗯，谢谢，还没算过过渡态，直接由gjf产生的坐标可以吗。尝试中谢谢指导！
赞一下 回复此楼
4楼2010-11-11 15:51:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhaoyxcas

捐助贵宾 (正式写手)
引用回帖:
Originally posted by manson1998 at 2010-11-11 11:40:15:
貌似是98的格式？

不太清楚，分子很大的话也适合用内坐标表示吗？看着很复杂！
赞一下 回复此楼
5楼2010-11-11 16:03:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhaoyxcas

捐助贵宾 (正式写手)
引用回帖:
Originally posted by zhangmt at 2010-11-11 13:49:43:
当然是坐标，只不过不是你通常看到的笛卡尔坐标，而是内坐标。也就是用键长、键角、二面角来描述的分子坐标。
由于是用QST2方式找过渡态的，所以分别给了反应物和产物两套坐标，让gaussian根据提供的情况去猜测、 ...

反应物有两个内坐标可以写在一起吗？
赞一下 回复此楼
6楼2010-11-11 16:05:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhaoyxcas

捐助贵宾 (正式写手)
引用回帖:
Originally posted by manson1998 at 2010-11-11 17:03:54:


qst2是由反应物和产物来计算过渡态的，所以输入文件里当然包含了两个分子的坐标。例子里分的很清楚用不同的title card隔开了。至于内坐标 或者 笛卡尔坐标，用谁都可以 无所谓。

我是说的反应物有两个，SiH2 和H2，这两个均在反应物中，就直接把坐标写在一起吗，是不是由结构直接转换过来的坐标直接紧挨着放在一起。中间不用分割或其他标注。感谢+麻烦！！
赞一下 回复此楼
9楼2010-11-12 11:46:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhaoyxcas

捐助贵宾 (正式写手)
引用回帖:
Originally posted by dx357 at 2010-11-11 16:39:20:
同惑。。。

我惑很多，希望你能弄明白，帮我解惑！！
赞一下 回复此楼
10楼2010-11-12 11:48:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhaoyxcas

捐助贵宾 (正式写手)
引用回帖:
Originally posted by manson1998 at 2010-11-12 12:25:53:


建议你google一下qst2在gv下的实现操作，或者看一下gv的手册。gv是直接可以保存为qst2型输入文件的。
当然你也可以将你认为符合qst2要求的两个分子的坐标手动添加到文本里，格式就是
#

title

0 1
co ...

嗯，非常感谢。呵呵是自己太笨了，看了有时也弄不明白。好好学习争取自己能算出过渡态。有困难.....想办法要解决地。
赞一下 回复此楼
12楼2010-11-12 22:25:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhaoyxcas

捐助贵宾 (正式写手)
引用回帖:
Originally posted by dx357 at 2010-11-15 15:53:19:



hehe，我也是刚学，问题也很多。就知道反应物中有两个分子，用高斯view构造时，把反应物中两个分子放的足够远，再优化。至于优化后的坐标两个分子就一个内坐标表示的。产物也是这样。

收益了！您的意思是用GUSSVIEW构筑好两个反应物，放的足够远，优化后应保证还是两个不相连的分子，然后将两个分子作为一个分子找出原子坐标作为起始物。（那直接用笛卡尔坐标不用内坐标也可以吧？）产物优化后。然后将其进行过渡态的优化？？说得有些啰嗦不知我表达的意思对不对？？
赞一下 回复此楼
16楼2010-11-15 17:10:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhaoyxcas

捐助贵宾 (正式写手)
引用回帖:
Originally posted by manson1998 at 2010-11-12 12:25:53:


建议你google一下qst2在gv下的实现操作，或者看一下gv的手册。gv是直接可以保存为qst2型输入文件的。
当然你也可以将你认为符合qst2要求的两个分子的坐标手动添加到文本里，格式就是
#

title

0 1
co ...

顺道，问个其它问题，homo及lumo的能量如log文件中的Alpha  occ. eigenvalues --   -0.22223  -0.21668  -0.16484  -0.15357
Alpha virt. eigenvalues --   -0.05149  -0.01632  -0.00178   0.00071   0.00507
在后面这么多数值，文献报道的一般指其后的哪一个数，以前在小木虫上问过不过没找到答案。谢谢帮忙！！
赞一下 回复此楼
18楼2010-11-15 17:14:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhaoyxcas

捐助贵宾 (正式写手)
引用回帖:
Originally posted by fatpig8832 at 2010-11-15 17:39:41:


莫非不是occ的最后一个和vir的第一个...

呵呵...是的，刚弄明白。
赞一下 回复此楼
20楼2010-11-15 18:41:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhaoyxcas

捐助贵宾 (正式写手)
引用回帖:
Originally posted by dx357 at 2010-11-16 10:01:49:


前面的说的都对，就是到优化后保证还是两个不相连的分子。然后分别要将优化后的反应物和产物（不是有两个输出文件吗，）然后把输出文件的坐标（两个）中间的题目之间用空行隔开形成一新的输入文件，关键词是q ...

反应物是两个，一起优化，优化后找出分子坐标，再和优化后的产物坐标在一起做为一个输出文件，我这样理解你的意思不知对不对？？期待回复。
赞一下 回复此楼
22楼2010-11-16 11:40:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhaoyxcas

捐助贵宾 (正式写手)
引用回帖:
Originally posted by fatpig8832 at 2010-11-15 17:39:41:


莫非不是occ的最后一个和vir的第一个...

nbo计算完成后出来的             Natural  -----------------------------------------------
    Atom  No    Charge         Core      Valence    Rydberg      Total
-----------------------------------------------------------------------
     Al    1    1.11105      

是说Al的电荷为+1.11105au?  这个数值与其价态有什么关系？？
谢谢！！！
没人回复，直接问你吧，谢谢了！！
赞一下 回复此楼
23楼2010-11-16 11:47:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhaoyxcas

捐助贵宾 (正式写手)
引用回帖:
Originally posted by dx357 at 2010-11-16 17:39:36:


嗯，是的。是一起作为输入文件的（注意两坐标之间题目前后要有空行）。给你个意见，不用找分子坐标。用高斯view新建，然后复制output分子构型到新建。然后点view file就ok。不知道我表达到位了没。因为是菜鸟 ...

嗯，已经很明白了，我试试看吧，很感谢，谢谢指教。
赞一下 回复此楼
27楼2010-11-16 19:56:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhaoyxcas 的主题更新
141/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭