将信息论与密度泛函中的一些概念相结合，在讨论化学问题中已有不少应用，特别是近年来北卡罗来纳大学&湖南师大的刘述斌教授在这方面做了大量工作，通过信息论定义的量（information-theoretic quantities）研究了诸多化学问题，获得了丰硕的成果，陆续发表了几十篇应用性或纯理论文章。之前我写过的介绍文章《用Hirshfeld电荷和信息增益研究亲电性和亲核性》（http://sobereva.com/257）就属于这个范畴。2019年末，刘述斌教授等人还在理论...

使用键级密度(BOD)和自然适应性轨道(NAdO)图形化研究化学键 文/Sobereva@北京科音  First release: 2020-Mar-2  Last update: 2020-Nov-28 1 BOD和NAdO简介 键级是十分重要的定量考察化学键特征的手段，这在《Multiwfn支持的分析化学键的方法一览》（http://sobereva.com/471）中已经充分介绍过。广义来说，离域化指数...

无需导数的局部极小化算法NEWUOA在Fortran中的使用简例 文/Sobereva@北京科音  2020-Mar-1 1 NEWUOA简介 数学上对多维函数做局部极小化的算法非常多，诸如simplex、BFGS、最陡下降法、共轭梯度法等等。NEWUOA（NEW Unconstrained Optimization Algorithm）是已故的剑桥大学教授Powell于2004年提出的一种不需要导数信息的非约束性局部极小化算法...

谈谈VMD可视化程序的连接关系的判断和设置问题 文/Sobereva@北京科音  2020-Feb-25 化学体系可视化程序VMD的用户越来越多，经常有人在网上问我怎么在VMD里修改键连关系，这里就专门说一下这个问题。本文内容对应于撰文时最新的VMD正式版1.9.3，对其它版本可能适用也可能不适用。 当VMD载入pdb、xyz等格式的结构文件时，程序就会自动判断连接关系。判断规则和Multiwfn、GaussView等程序类似，...

通过轨道权重福井函数和轨道权重双描述符预测亲核和亲电反应位点 文/Sobereva@北京科音  2020-Feb-24 0 前言 福井函数（Fukui function）与双描述符（dual descriptor）是概念密度泛函理论框架下定义的非常流行的用于预测反应位点的方法，笔者之前写过很多相关文章，连同大量相关综述都汇总提供在了这里：《概念密度泛函综述和重要文献合集》（http://bbs.keinsci.com/threa...

懒得在这里再发一遍了，本文直接在计算化学公社论坛看： http://bbs.keinsci.com/thread-16255-1-1.html 此文目的是介绍怎么用molclus基于xtb做的分子动力学轨迹作为初始构型，搜索治疗武汉肺炎的特效药瑞德西韦的在水溶液下的能量最低构象，最终结果如下。此文的做法是对较大体系构象搜索的极佳、普适做法。 后注：鉴于有人在病毒名称上特别敏感，笔者在此进行如下专门澄清，以后不要再试图给笔者发邮件在肺炎的用...

RESP2原子电荷的思想以及在Multiwfn中的计算 文/Sobereva@北京科音  2020-Feb-1 1 前言 RESP是几乎最适合用于分子动力学模拟的原子电荷计算方法，在《RESP拟合静电势电荷的原理以及在Multiwfn中的计算》（http://sobereva.com/441）中笔者对RESP电荷有非常详细的介绍。RESP电荷的具体数值明显依赖于计算电子波函数用的量子化学方法和基组，以及对溶剂效应的考虑。有不少人...

一键把所有gjf文件转成xyz文件、把所有Gaussian输出文件转成gjf文件的脚本 文/Sobereva@北京科音 First release: 2020-Jan-28  Last update: 2020-Feb-5 本文提供笔者写的两个Linux下的Bash shell脚本，对实际研究有一定价值，都利用了Multiwfn程序进行格式转换，因此Multiwfn必须已恰当安装从而能通过Multiwfn命令直接启动。Mult...

使用Multiwfn便利地查看所有激发态中的主要轨道跃迁贡献 文/Sobereva@北京科音 First release: 2020-Jan-28  Last update: 2020-Jul-7 之前笔者写过一篇文章《电子激发任务中轨道跃迁贡献的计算》（http://sobereva.com/230），介绍了如何对TDDFT及类似方法计算的激发态计算轨道跃迁的贡献。在很多电子激发的研究文章中都会给出各个激发态中各种主要MO跃...

2019年下半年观看的非TV版动画简评 文/Sobereva  2020-Jan-20 这里对从2019年7月1日开始截止到2019年12月31日观看的非TV版动画进行简单评论。还有些此类动画看完后忘写了，就也不再补了。由于这类动画不好对比和打分，我也就不给出评分了。以往的此系列文章和TV版动画评论见http://sobereva.com右方的“动画评论”分类。 关于我转生变成史莱姆这档事 OAD01 长度：23m16s 讲述...

2019年下半年优秀动画音乐推荐 文/Sobereva  2020-Jan-19 这是第11次写“优秀动画音乐推荐”系列文章，此系列从2014年起不断发布，致力于尽自己微薄之力，令最杰出、最美妙的音乐能有出头之地，而不被数目是它们20倍以上的海量的平庸、凑数的音乐埋没。之前的动画音乐推荐见 《2019年上半年优秀动画音乐推荐》（http://sobereva.com/496） 《2018年下半年优秀动画音乐推荐》（http://sobere...

使用PSI4做对称匹配微扰理论(SAPT)能量分解计算 文/Sobereva@北京科音 First release: 2019-Dec-24  Last update: 2020-Oct-15 1 前言 能量分解方法之前我在《Multiwfn支持的弱相互作用的分析方法概览》（http://sobereva.com/252）中简单介绍过。能量分解就是将片段间相互作用能分解成不同物理成分，从而能从能量角度更深入地了解相...

通过轨道离域指数(ODI)衡量轨道的空间离域程度 文/Sobereva@北京科音 First release: 2019-Dec-22  Last update: 2020-Mar-22 1 前言 轨道的离域程度就是指轨道分布的范围是比较窄还是比较广，分布范围越大可以说离域性越强。有人在思想家公社QQ群里发了几个轨道图片，问怎么衡量他这几个轨道的离域程度，即谁离域相对比较强一些。有很多情况直接看轨道图形就可以直接比...

后记：后来笔者在此文的研究的基础上做了大量扩展和深化，发表了一系列论文，并且还做了很多其它的与18碳环及其类似物有关的研究，汇总见http://sobereva.com/carbon_ring.html。非常欢迎阅读和引用其中的文章！ 一篇最全面、系统的研究新颖独特的18碳环的理论文章 Sobereva@北京科音 First release: 2019-Dec-15  Last update: 2020-Jan-3 ...

将Multiwfn计算的键级直接标注在分子结构图上的方法 文/Sobereva@北京科音  2019-Dec-4 有人在Multiwfn中文论坛上问怎么把Multiwfn计算的键级直接标注在分子结构图上，以便于直观观看。将Multiwfn与GaussView结合使用就可以很方便地实现这一点，在本文介绍一下。读者必须使用2019-Dec-4及以后更新的Multiwfn版本，可以在其官网http://sobereva.com/mult...