fatpig 发表于 2023-2-28 18:21 明白了,谢谢你 |
HJLCZ 发表于 2023-2-27 20:04 1. 就是建议至少把基组加大到6-31g(d,p),多加的p是给H的,或者至少指定给那个要反应的H。 可以看http://sobereva.com/336 2. 如果你第一点的初始猜测就错了,你算后面的有啥用...所以建议先做stable=opt;接着做freq,可以在gaussview里看到你第一点的振动方向,而你用calcfc是没法直接看到的。 |
chands 发表于 2023-2-27 21:18 好的谢谢老师 |
HJLCZ 发表于 2023-2-27 21:05 建议还是先查相关资料。注意到你是做气相反应,(气态)羧酸在空气中相当稳定。另外这个反应就算有,能垒也会非常高。 |
chands 发表于 2023-2-27 20:40 老师是这样的,我们专业不是化学专业,这个反应在本专业的一些文献当中有所体现,但是说的很笼统,文献中只是给出了反应能垒图,不知道本领域的审稿人有没有验证他这个反应的过渡态的真实性。我的目的是想把这个反应做的详细些,然后像老师后面说的,得到PH-CH2-COO·,分解生成PH-CH2·和CO2。 |
wzkchem5 发表于 2023-2-27 20:20 谢谢老师,博文中提到gview中的成键关系不影响量化计算的结果。那想请问一下老师,对于我上述反应,找过渡态一直不收敛的原因,老师能说一下可能是什么原因导致的吗? |
| 这个反应很奇怪,你确定这个反应会发生吗?你确定机理是这样吗?另外,Ph-CH2-COO·容易脱羧,生成Ph-CH2·和CO2。 |
| 仔细阅读http://sobereva.com/414,了解为什么双键不是导致不收敛的原因 |
本帖最后由 HJLCZ 于 2023-2-27 20:19 编辑 fatpig 发表于 2023-2-27 19:33 谢谢你的解答,对于您的回答有如下疑问: 1.因为刚接触计算不久,不太明白你说的第1点是什么意思,能麻烦您说的稍微详细点吗? 2.你的意思是对于我初始摆的过渡态的结构进行stabl=opt,看看波函数是否稳定,再对该结构进行freq操作,若振动方向与我预期相同,是不是可以直接读取力常数进行opt=ts的操作。想问一下这样操作和我帖子里面那样操作有什么不一样的地方吗? 因为刚接触量化计算,若有些常识性的东西问的不合适,还请多多包涵。 |
|
如果是我的话: 1. 既然是H transfer,至少加个p给H吧... 2. 做单点的stable=opt; 3. 读上一步的guess,单点做freq,看虚频的振动模式是否符合预期; 4. readfc读上一步的力常数找TS. |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-2-24 07:12 , Processed in 0.233797 second(s), 25 queries , Gzip On.