xxzj 发表于 2019-6-1 11:00 是 你哪能用绝对能量跟相对能量比 跟是否考虑溶剂没必然关系。当然溶剂下的情况应当考虑溶解自由能 |
本帖最后由 xxzj 于 2019-6-1 11:04 编辑 sobereva 发表于 2019-6-1 08:53 老师，您不要生气，原谅我没有认真看帖子，我认识到错误了， 我重新修改自己的问题：从opt freq得到的Thermal correction to Gibbs Free Energy值为0.191661，然后DFT单点计算HF=后面值为-708.2937197，两者相加得到-708.1020587，就是自由能值吗，但是单位是不是Hartree，然后单位换算结果为-4.44*(10^5)kcal/mol,与实际相差太远了，老师帮看一下哪里出现错误了，还是我计算错误了，是不是还需要计算溶液的，然后按照上述方法计算，最终两个自由能结果相减，才是最终的自由能值，然后通过最后的这个值比较稳定性，辛苦老师了 |
xxzj 发表于 2019-5-31 09:17 你就是不看我贴的文章，还用“HF值”这种说法，我拒绝回复 |
本帖最后由 xxzj 于 2019-5-31 21:10 编辑 sobereva 发表于 2019-5-31 07:04 老师，从opt freq得到的Thermal correction to Gibbs Free Energy值为0.191661，然后单点计算的HF值为-708.2937197，两者相加得到-708.1020587，就是自由能值吗，但是单位是不是Hartree，然后单位换算结果为-4.44*(10^5)kcal/mol,与实际相差太远了，老师帮看一下哪里出现错误了，还是我计算错误了，是不是还需要计算溶液的，然后按照上述方法计算，最终两个自由能结果相减，才是最终的自由能值，然后通过最后的这个值比较稳定性，辛苦老师 |
xxzj 发表于 2019-5-30 09:17 仔细看此文 谈谈该从Gaussian输出文件中的什么地方读电子能量 http://sobereva.com/488 用自由能对比 |
本帖最后由 xxzj 于 2019-5-30 18:58 编辑 sobereva 发表于 2019-5-30 08:46 老师，那比较两个同分异构体的稳定性时，是不是可以用气相自由能对比得出结论或者简单的做过单点进行能量比较就可以，辛苦老师了，老师我现在用B3LYP-D3/6-31G*进行opt freq计算，然后比较最后HF值的大小，哪个小说明哪个分子能量稳定，可以吗 |
xxzj 发表于 2019-5-29 19:23 1 算的是形成氢键过程的自由能变 2 显然不正确。你单点和几何优化用相同的级别，你做单点还有什么意义？仔细看 谈谈隐式溶剂模型下溶解自由能和体系自由能的计算 http://sobereva.com/327（http://bbs.keinsci.com/thread-3345-1-1.html） 3 跟D3根本就没直接关系。显然不能 仔细看 量子化学计算中帮助几何优化收敛的常用方法 http://sobereva.com/164 谈谈“计算时是否需要加DFT-D3色散校正？” http://sobereva.com/413（http://bbs.keinsci.com/thread-9772-1-1.html） |
本帖最后由 xxzj 于 2019-5-29 19:28 编辑 sobereva 发表于 2019-5-29 06:08 谢谢老师的耐心回答，我还有两个问题想请教老师， 1）为了比较分子内存在氢键的两个分子的稳定性，是不是可以通过计算它们气相下自由能值就能得到结果；2）在计算气相下自由能时，输入文件#p m062x/6-31+G* empiricaldispersion=gd3 opt freq ，然后读取Thermal correction to Gibbs Free Energy的数值，然后用m062x/6-31+G** empiricaldispersion=gd3 sp 计算单点，读取输出文件中最后的HF能量，然后最后与Thermal correction to Gibbs Free Energy的数值相加，作为气相下自由能的结果，然后比较自由能大小，哪个能量低哪个稳定，这样做是正确的吗，因为需要修稿，是实验理论相结合的，怕出现错误，辛苦老师帮我分析一下，我上面的计算和描述是否合理，如果不合理，是否有其他的方法 3）用m062x/6-31+G** empiricaldispersion=gd3进行优化频率计算时，不能收敛，D3校正是不是可以不加，然后计算单点时加入，辛苦老师了 |
xxzj 发表于 2019-5-28 20:27 1 键连关系根本不影响量化计算结果 2 应当算自由能对比 3 用这个级别100%文章发不出去，只会被内行审稿人骂一顿然后拒掉。仔细看完以下文章再做弱相互作用相关计算，否则算出的数据毫无意义 谈谈“计算时是否需要加DFT-D3色散校正？” http://sobereva.com/413 DFT-D色散校正的使用 http://sobereva.com/210 乱谈DFT-D http://sobereva.com/83 谈谈BSSE校正与Gaussian对它的处理 http://sobereva.com/46 4 仔细看此文 谈谈原子间是否成键的判断问题 http://sobereva.com/414（http://bbs.keinsci.com/thread-9840-1-1.html） |
本帖最后由 xxzj 于 2019-5-28 20:38 编辑 小范范1989 发表于 2016-11-20 19:21 老师，我有几个问题有些不清楚，请老师帮我解答一下： 1）如果分子内一个片段上的氧原子与另一个片段上的氢原子形成具有氢键的六圆环，那么在gaussview中建模时，O与和她形成氢键的H原子之间是需要用虚线，单键还是二者之间不用任何链接，然后进行基态几何构型的优化 2）比较两个同分异构体的稳定性时，是否可以之间取基态几何构型优化的out文件中的HF能量进行对比，看哪个分子能量低则哪个分子稳定 3）含有分子内氢键时，用B3LYP/6-31G*进行优化计算合适吗，还是选择B3LYP/6-311G**的方法； 4）当分子内存在氢键时，但是我在建模的时候，O与另一个片段上的H原子之间没有任何链接，但是优化的构型中，原来另一个分子片段与H原子之间的键突然没有了，是不是说明O与H之间氢键较强，然后由于O原子的吸引，所以导致H原子与之前连接的分子键断裂；问题有点多，辛苦老师了，学的基础有些不扎实，辛苦老师了  |
sobereva 发表于 2016-11-21 12:23 恩，谢谢sob老师指点。谢谢老师 |
小范范1989 发表于 2016-11-21 08:04 1 是。如果你不做freq只用单点能求差，就当成是对焓变的近似讨论。 2 是，但这是最懒、最不精确的做法。 |
本帖最后由 小范范1989 于 2016-11-21 08:24 编辑 sobereva 发表于 2016-11-21 06:19 谢谢sob老师指点。 1：也就是说，需要计算阴离子和H+的freq喽？都要做freq分析才能找到焓变。能不能对阴离子不做freq计算，近似采用它的文件最后的HF能量呢？ 因为，我要看的是激发态的解离能，这个如果对阴离子激发态做频率，有点慢。不知道这样近似可以不？或者有什么其他更好的办法， 2：如果计算了频率，焓变是不是看输出中的： sun of electronic and thermal Enthaipies 这个能量呢？ 谢谢sob老师指点 |
小范范1989 发表于 2016-11-20 19:21 解离之后形成阴离子和质子，显然是两个产物的焓减去一个反应物的焓，不可能光是阴离子减中性状态。 用经典统计力学，H+的焓是5/2*RT，可以不用算，你用高斯做freq也会给你这个值 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-2-25 01:55 , Processed in 0.177766 second(s), 25 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites