Chess1357 发表于 2025-12-27 16:37 原理上对于任意分子都可以。但是巨大的生物分子体系基于波函数算静电势算不动,只能基于原子电荷算,如下文 基于原子电荷极快速绘制超大体系的分子表面静电势图 http://sobereva.com/639(http://bbs.keinsci.com/thread-28913-1-1.html) |
Chess1357 发表于 2025-12-27 16:36 Multiwfn主功能12对静电势分析完后,后处理菜单选择导出表面顶点成为vtx.pqr,放到VMD里显示,选择语句按照charge属性(当前对应于静电势)设置。诸如选择语句用charge < -0.03 代表只显示表面静电势值< -0.03 a.u.的表面顶点部分。 |
Chess1357 发表于 2025-12-27 16:36 另外生物大分子比如核酸蛋白质,也可以参照文中的方法进行MPI计算吗?期待指点,谢谢!!! |
| 真是一篇非常出色的介绍!现在我想知道运用MPI方法衡量分子极性时,能否将极性表面和非极性表面进行可视化呢按照文中的定义?谢谢sob老师!!!! |
sobereva 发表于 2025-11-17 20:12 啊我想明白了,谢谢老师指点 ![]() |
sgwzq 发表于 2025-11-17 16:26 最低阶非零多极矩没有原点依赖性 甲烷最低阶非零多极矩是四极矩,所以不依赖于原点。而偶极矩不为0的分子的四极矩就有原点依赖性 |
| sob老师您好,我有一个问题:电四极矩的定义是电荷量,乘以直角坐标平方,再求和,如果坐标选取不同,电四极矩的值应该不一样吧?如果建模的时候把甲烷质心放在(0.0,0.0,0.0),和把它放在(1.0, 0.0, 0.0),Gaussian应该会给出不一样的值吧?这样比较是不是就没有意义了? |
EmperorXiang 发表于 2024-3-9 21:05 不是导出,而是需要额外计算 后HF波函数不是后HF能量的直接的副产物,就如同你计算能量也不会免费得到能量的导数 |
sobereva 发表于 2024-3-9 05:44 不过为什么post-HF波函数的导出会花费那么长的时间呢? |
sobereva 发表于 2024-3-9 05:44 明白了,感谢社长 |
wzkchem5 发表于 2024-3-9 00:14 但测试时用CCSD和CCSD(T)单点算出来的偶极矩和HF是一样的。是我的设置出了问题吗?比如我用CCSD(T)/aug-cc-pVQZ//CCSD/cc-pVTZ计算得到的水分子的偶极矩是1.99D,而HF/aug-cc-pVQZ//CCSD/cc-pVTZ得到的也是这个值。 |
EmperorXiang 发表于 2024-3-8 16:39 足够高阶的post-HF方法不会高估偶极矩,比如CCSD(T)一般问题不大 |
sobereva 发表于 2024-3-8 05:10 确实有注意到这一点,似乎连带一众post-HF都高估偶极矩。那应该使用什么方法? |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-1-23 20:13 , Processed in 0.177270 second(s), 26 queries , Gzip On.