计算化学公社

 找回密码 Forget password
 注册 Register

简单量子化学问题答疑专帖

查看数: 567631 | 评论数: 5022 | 收藏 Add to favorites 31
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2022-11-28 01:57

正文摘要:

经常看到有些人在其它帖子后面跟帖提问和主贴无关的比较简单的量化问题,本论坛很不鼓励这样的做法,楼就歪了,会造成讨论秩序的严重混乱。但有些问题非常简单,单独发个帖子问也不值得(要不然论坛平均讨论质量会显 ...

回复 Reply

KCM@SHU 发表于 Post on 3 hour ago
wzkchem5 发表于 2026-3-9 19:04
如果一个周期性体系的电荷数不为0,说明体系结构或者电子数搞错了,例如多氢、漏氢、多电子、少电子等。 ...

老师那如果像他这样算MOF结构的时候,从里面选取高度对称的片段作为代表,是应该把一些漏氢的部分补上让其带正电来算,还是不修改,片段带什么电就设置什么电荷去算呢?
wzkchem5 发表于 Post on 3 hour ago
KCM@SHU 发表于 2026-3-9 13:33
各位老师我最近想算ZIF-8的吸附,然后我看到一篇文章他用ZIF-8的γ点取样法选取一个片段最后用Gaussian去做 ...

如果一个周期性体系的电荷数不为0,说明体系结构或者电子数搞错了,例如多氢、漏氢、多电子、少电子等。不难估计出,只要是DFT能算得动的大小的晶胞,晶粒大到能用肉眼看见,但凡每个晶胞带一个电,库伦排斥能都足够让晶体原地爆炸了,这样的晶体根本存在不了
KCM@SHU 发表于 Post on 9 hour ago
本帖最后由 KCM@SHU 于 2026-3-9 13:36 编辑

各位老师我最近想算ZIF-8的吸附,然后我看到一篇文章他用ZIF-8的γ点取样法选取一个片段最后用Gaussian去做计算优化,如下图,我想问一下这样从周期性结构里面截取代表性片段计算的时候需不需要考虑片段电荷呢?(印象里面大多数周期性结构似乎不考虑带电?),正常情况下咪唑上一个N会脱去氢带一个负电,但是按照他这个片段,有六个咪唑,而且都脱去氢了,应该有六个负电,然后两个锌离子带+4,是不是gaussian计算的时候应该设置电荷为-2比较合适?(文章计算细节里面没提电荷)原文DOI:10.1021/acs.est.3c10653

202603091331548259..png (511.62 KB, 下载次数 Times of downloads: 0)

202603091331548259..png
sobereva 发表于 Post on yesterday 09:23
KCM@SHU 发表于 2026-3-7 10:30
各位老师我想请问一下可以根据HOMO-LUMO的轨道匹配程度用来预测有机配体和金属离子结合能力大小吗?我看好 ...

诸如此类能轻易较准确定量计算的就没必要去预测,参考
谈谈分子间结合能的构成以及分解分析思想
http://sobereva.com/733http://bbs.keinsci.com/thread-50890-1-1.html


基于轨道的分析更适合拿来解释定量计算的结果
Stardust0831 发表于 Post on 前天 21:57
KCM@SHU 发表于 2026-3-7 10:30
各位老师我想请问一下可以根据HOMO-LUMO的轨道匹配程度用来预测有机配体和金属离子结合能力大小吗?我看好 ...

不要过于高估和神化前线轨道理论。
想算结合能力直接算个binding energy就足够直接了吧。
相互作用分析还可以参考卢天老师的这篇综述:《Angew. Chem.上发表了全面介绍各种共价和非共价相互作用可视化分析方法的综述
KCM@SHU 发表于 Post on 前天 10:30
各位老师我想请问一下可以根据HOMO-LUMO的轨道匹配程度用来预测有机配体和金属离子结合能力大小吗?我看好像目前HOMO-LUMO轨道目前都是用来预测催化反应活性,不知道这个是不是同样可以用来预测配位结合。另外还想问一下各位老师有没有其他推荐的参数用来预测有机配体对金属离子吸附能力预测,或者推荐一些相关的理论书籍文献
sobereva 发表于 Post on 2026-3-1 21:07:38
Huijue 发表于 2026-3-1 13:47
卢老师陈老师好,在thermo这部分,做练习题Fe2+与(C5H5)-反应生成重叠型二茂铁时,我建立的二茂铁初始结构 ...

和版本有关。诸如G09和G16默认的积分格点精度就足以导致这种程度的差异。并且如果初始结构不精确相同,由于收敛限不是无限严格的,也可能最终带来这样的差异

Huijue 发表于 Post on 2026-3-1 13:47:02

重叠型二茂铁优化结果

卢老师陈老师好,在thermo这部分,做练习题Fe2+与(C5H5)-反应生成重叠型二茂铁时,我建立的二茂铁初始结构与您的一样(重叠型,D5h的point group),Gaussian跑出来的结果,在GaussView中Results显示的vibrations却有小小的差异,您给的结果前两个vibrations是39cm-1, 169cm-1...,而我的是43cm-1, 172cm-1... 这种情况是错在哪里了呢?Ps,用的理论方法和基组都是您练习题上要求的。谢谢
3079987747 发表于 Post on 2026-2-18 16:00:19
sobereva 发表于 2026-2-18 06:37
这不属于“简单量子化学问题答疑专帖”的范畴,请单独发主题

感谢sob老师的指正与建议,得益于sob老师的指导,本人已经解决现在的问题,以后一定多利用现在的资源,多多学习,避免再出现相似的问题!
sobereva 发表于 Post on 2026-2-18 06:37:25
Xx_DerEnTh 发表于 2026-2-17 22:02
各位老师好,如果一个反应路径存在势能面交叉,例如最低反应路径是三重态的产物复合物经过单重态的过渡态生 ...

这不属于“简单量子化学问题答疑专帖”的范畴,请单独发主题
Xx_DerEnTh 发表于 Post on 2026-2-17 22:02:15
各位老师好,如果一个反应路径存在势能面交叉,例如最低反应路径是三重态的产物复合物经过单重态的过渡态生成单重态的产物,为了对比该路径中三个结构的能量,三重态结构活性空间的选择需要与单重态结构活性空间的选择保证一致吗(均是(m,n))?我发现好像无法实现,还是说此时三重态结构活性空间的选择不受单重态结构活性空间的影响,只需要保证两个单重态结构活性空间一致即可呢?
3079987747 发表于 Post on 2026-2-14 23:41:33

硫酸根离子配合双PANI分子 ,以L502报错,解决方法都试过然而无效

各位老师好,本人在计算硫酸根离子配合双PANI分子时,每次都会出现报错,前几次出现了各种各样的报错,本人已解决,但是现在一直出现L502,不收敛的报错,尝试各种解决方案均以失败告终,请问各位老师,是分子结构的问题还是其他方面的问题


PANI SO4 2.gjf

5.89 KB, 下载次数 Times of downloads: 1

PANI SO4 2.LOG

354.55 KB, 下载次数 Times of downloads: 2

sobereva 发表于 Post on 2026-2-8 05:33:35
moiumuio 发表于 2026-2-8 00:52
Hi teacher,
I just want to ask, should I use rigid PES (only do SP) or relaxed PES (Opt and SP)?
...

Performing rigid scan is the easiest.
Relaxed scan may be feasible using the GIC feature of Gaussian. GIC allows users to perform a relaxed scan for a user-defined coordinate (including dihredral), which may involve center of ring(s). But the GIC code is not very robust and may fail.

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-3-9 22:52 , Processed in 0.218321 second(s), 27 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list