GHL 发表于 2022-10-31 00:11 糖苷的话，每根键都扭转就比较麻烦了，所以最好按我说的，宁可不穷尽所有排列组合，也要保证所有键都有扭转的机会，而不是只允许个别键扭转但是把这些键扭转的所有排列组合都搜索到。 |
GHL 发表于 2022-10-30 15:17 最好所有旋转都考虑，但不需要搜索所有可能的排列组合。也就是说，宁可让所有键都可以旋转、但不严格保证搜索到最优构象，也不要只让一部分键旋转、但是保证搜索到只让这些键旋转的情况下的最优构象。 当然话说回来，你的体系只有5根键可以旋转的话，用xtb等比较快的方法完全搜索所有的排列组合也是完全做得到的 |
GHL 发表于 2022-10-24 05:38 这个只是确定是局部极小，不能确定是全局极小。这是最最基础的数学问题，仔细看https://en.wikipedia.org/wiki/Second_partial_derivative_test |
GHL 发表于 2022-10-30 11:52 M06-2X方法在高斯里的关键词叫做M062X。那个横杠在输入文件里不能写，但是在文章里必须写。 |
wzkchem5 发表于 2022-10-21 16:52 您好，大佬，我后面采用的是confab+molculs 做的构象搜索，后期采用B3LYP 6-31**g做的分子优化，单点能采用的是同一水平的6-31**gd, 后期依照能量最小的选取（负值最大），计算的结果不是太理想，想采用之前您建议的M06-2X方法，发现gaussian09 中没有选项，不知道怎样处理？ |
|
The energies of all geometries of anthocyanin-derived pigments and their respective radicals were calculated in vacuum with the density functional theory (DFT) approach, at the unrestricted B3LYP38−40 hybrid density functional level (Becke−Slater−HF exchange with Lee−Yang−Parr correlation functional) using the 6-311+G(d) basis set as implemented in the Gaussian 09 package.41 The inclusion of polarization functions correctly takes into account intramolecular H bonds. The optimized structures were confirmed as true minima by vibrational analysis performed at the same level of theory 好像他们在计算前对结构进行振动分析确保结构最小，不知道这一步是怎样做到了，麻烦大佬了! |
GHL 发表于 2022-10-23 12:53 即使是算垂直电离势，严格的做法也涉及到构象搜索，不然你不知道你优化的电离前的分子的构象对不对。他们没提到构象搜索，更可能是因为早年没有CREST、molclus等专业构象搜索软件，也没有xtb等适合做构象搜索的半经验方法，做构象搜索不方便，所以很多人牺牲严谨程度换取便利性。但是现在这些工具有了以后，不做构象搜索就越来越不被审稿人接受了。 垂直电离势修改电荷、自旋多重度以后，只算单点能。注意结构优化并不是得到体系能量的唯一方式，对于计算垂直电离势而言，电离后的分子的能量只能从单点计算获得，而不能从结构优化获得。只要知道垂直电离势的定义就能看出为什么是这样的 |
wzkchem5 发表于 2022-10-20 21:55 刚看到垂直电离势是能量的差值，而绝热电离势是焓的差值 |
本帖最后由 GHL 于 2022-10-23 19:54 编辑 wzkchem5 发表于 2022-10-20 21:55 大佬，不好意思，再次打扰 我查阅很多英文文献，在计算IP时，结果均比较好，会不会是因为他人计算的是垂直电离势而我计算的是绝热电离势的原因。 因为如果是后者的话，肯定都有提及到构象搜索，但是都没有提及。 如果是前者的话我对于“垂直电离势“”只优化不修改，不再次优化”这个不是太清楚 是改变电荷，自旋多重度，只做结构优化吗？ 计算IP我都只做了OPT+Fre |
wzkchem5 发表于 2022-10-21 16:52 谢谢大佬关照，谢谢！！！ |
GHL 发表于 2022-10-21 09:07 基组不用加弥散，除非你之后还要研究这个分子带负电的状态（即使那样也只需要对重原子加弥散，不需要对氢加弥散）。泛函B3LYP怕有自相互作用误差问题，未必适用（但也未必不适用），可能M06-2X合适一点。其他没问题。 构象搜索和结构优化是两个概念，构象搜索要用http://bbs.keinsci.com/thread-577-1-1.html的方法做。仔细看这篇文章，了解为什么结构优化不足以保证得到最低能量构象 |
wzkchem5 发表于 2022-10-21 14:14 我想问下大佬，我的输入文件方法参数是否合理，计算IP时，采用的是离子态（电荷+2，自旋多重度2）-母体分子（电荷+1，自旋多重度1）. 构象搜索那块不知道如何实现，分子结构优化结果收敛，是不是已经将最小能量的构象记录？ 可能说的不对，别介意啊! |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-2-19 13:59 , Processed in 0.197827 second(s), 25 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites