☞走走的个人资料

【文献翻译】密度泛函理论的半个世纪课程作业，顺手贴上来了。Markdown写的。手翻，没顾上校订。

有没有一种可能，gtk库搞崩了会连BIOS启动界面都出现不了？前段时间在自己电脑上重装了一下Linux系统。现在是CentOS 8 Stream和Windows 10 21H2的双系统。有独立显卡和主板的集显。起因是为了在Linux下因为某些原因dnf安装和删除了一些GTK相关的库，重启就发现显示器（注，某4k显示器，下同）没有任何显示信号，包括BIOS的启动界面也没有显示信号。但是过几秒钟应该是grub2的自动引导可以成功进Linux系统，因为过一会以后可以用别的电脑远程登录这台电脑。一点一点排查吧。显示器用同一根线连笔记本可以正常显示。换显卡上不同的输出口无效。主板跳线改用集显的VGA输出没信号。用一台老旧显示器（16:10的，分辨率很低）VGA转HDMI接独显上无显示。又换了一台普通的1k显示器，从独显接上。这次BIOS正常，也能进Windows，但是Linux下只能用命令行界面登录(Ctrl+Alt+F2)可以登录，init 5然后Ctrl+Alt+F1没有和以前一样打开GUI的登录界面，依然是命令行的登录界面。怀疑是因为GTK库崩了导致的。请问有没有可能，这个BIOS界面都没有显示是Linux的GTK库崩了导致的？如果可能的话，应该尝试怎么解决？（因为没有BIOS显示，没法插U盘重装系统。）我目前是打算找一台别的电脑看看装的GTK相关的库都有哪些版本，然后替换一下...

请教一个特别诡异的bash初始化报错系统是CentOS 7.9 我们是一个集群，所有节点的/opt和/home都是挂载的同一个硬盘阵列的同一个目录，也就是各个节点的/opt和/home以及下面的子目录都是完全一样的。现在某个节点我用用户（假设叫myuser），不管是登录模式还是非登录模式，启动bash都会有个报错/yum: no such file or directory。我想那就是我.bashrc之类的文件手抖写错了什么呗，可是再一想不对，所有节点的/home/myuser下都是挂载的同一个目录，所以要是一个节点的登录有报错，别的节点应该也有啊。可是只有这一个节点会有这个报错。用root身份使用这个节点，不管是登录模式还是非登录模式，都没有这个报错。各个节点的/root是不同的目录。额外挂载硬盘阵列的只有各个节点的/home和/opt以及子目录。所有.bashrc，.bash_profile，BASH_ENV之类的都检查了，没发现/yum相关的语句。请问这个报错可能是哪个配置文件写错了？

求佬们帮忙看下这个配置有没有严重问题主要是兼容性问题，比如主板是否兼容PCI-E 4.0（这个查了应该能），机箱里能不能塞得下之类的。谢谢大家主板：微星 B660M DDR-4 Wifi型 CPU：Intel i7-12700，散片的板U套装 CPU散热器：利民PA120 TL-c12双风扇6铜管内存：三星DDR4 3200MHz 16GB*2 硬盘：三星PM9A1 1T。之前有一块西部数据的500G SATA SSD做系统盘，一块东芝的7200rpm 64MB 2T CMR HDD，都用了不到一年，拆下来接着用。显卡：之前买的Nvidia Geforce RTX 2060 电源：振华 LEADEX HG650W 金牌全模组机箱：先马蜜獾 M1商务版+3*先马风扇显示器：AOC U27P2C 键鼠有音箱有旧的别的线材之类的是否要自行购买呢？手头有一根DP的线1.2 144Hz 4k数据线。

【公告】封禁三小基甲苯 1天，并警示其他人近期以@掀起一阵冷为首的一系列账号，在吧里到处发布无意义回复，水贴、低俗内容，且主要是在别人正常楼贴下面回复，严重干扰了他人的阅读体验，同时有拉帮结派无视吧规的嫌疑，并在群内被警告后依然没有明显改进，还带领诸多账号发表类似内容。经过吧务组慎重讨论，一致决定给予封禁一天处罚，并对【其他从犯】提出严重警告，视其等人后续表现可能不再追究或者删帖、封禁、循环封禁甚至永封加黑名单。化吧不是某些人个人的群，发帖回帖者有义务充分考虑他人的阅读体验。极个别主题水贴可以容忍，但是在他人的正经贴下面水贴是很恶劣的行为。大量无意义甚至低俗讨论请自建群讨论。再次强调，这里是公共空间。最近发现大量账号和被封禁的此人有大量无意义互动，到处刷此类回复，使得恶劣影响扩大了好几倍。在此同时对这些账号提出严重警告。以后一旦发现再有如此行为，将同样甚至加重处理（因为已经警告过）。楼下将会陆续贴出此人和相关人士水贴的证据。不做打码处理。同时严禁楼下进行人身攻击。类似诸如“三小基甲苯是帅哥哥”“老公抱抱”和各种无意义表情等是上面指代的“上述行为”的主要特征。相关吧规截图。吧里很久以前就禁用了很长的小尾巴，同时非常不鼓励在任何情况下使用小尾巴，严重影响他人阅读体验。感谢大家的理解配合、为维护公共秩序和促进本吧氛围环境正面发展做出的贡献。

【科普】臭氧的电子结构和开壳层单重态 1L防抽

【个人公告】发帖请注意您的语气和态度提问者请务必注意，提问时应摆清自身定位，清楚你是来问问题的不是来发工资的。没有人欠你的，没有人有义务给你回答，也没有人有义务回复你。对回答者应该给予充分尊重。诸如对他人的解答给予适当的回应，感谢，尤其不要在前面已经有人回复过的情况下连看都不看一眼就接着说“人呢？”“人呢？”，这是对回答者的严重不尊重。严禁使用督促、吆喝、命令等语气提问。这种我个人会毫不犹豫地被删掉，视恶劣程度我个人可能会给予封号处理。这种贴每天都要处理掉非常多。反面例子： “来个人给我解答一下” “大神呢” “人呢？” “这么大贴吧没人么” “来个活人” “来人啊来人啊急急急” “都s光了么”(真的有过这么回复的) 不要没完没了顶自己的贴，严重影响版面和他人阅读。想想有人本来感兴趣你的问题，发现突然多了几十回复，结果点进来发现全是一堆“顶”“不要沉啊啊啊啊啊啊啊”“人呢”“来个人啊”，严重污染视线，也会失去回复的欲望。而且真正有用的回复会被淹没让别人没法一起阅读，降低版面质量。公共论坛不是你自己家，发贴者应充分考虑他人的阅读体验。经验表明，反复强调自己急急急的，99%不是作弊就是根本不急。公共论坛你写急急急并不能然你比别人获得更高的优先权。如果实在很急需要顶贴，自己在一层楼的楼中楼里面顶贴，这样可以最小化对版面的污染，方便别人阅读后面的内容。如果你在一个贴里面多个楼无意义顶贴，别人看到一堆回复点进来找不到有价值回复，谁知道你的问题是不是已经被解决了？其他人也无法从中学到任何东西。尤其上大学的，都是成年人了，你卖萌自称宝宝我管不着，但是提问标题还动不动“宝宝心里苦”“宝宝好难啊”的就让人很不舒服。这类问题99%都是极为基础的内容，但凡看过一遍书都不应该问。学习是你自己的事，能不能不挂科能不能毕业能不能找到好工作都和别人无关。不要把自己的任务退给别人。这类人和上面"来个人啊”这类有高度重合，懒到这样就别读书了，回家抱孩子去(希望不祸害孩子)。更有甚者，当别人指出这个问题是彻彻底底伸手的时候，居然对指出者言语攻击，这种直接黑名单，您下个ID学聪明点。伸手不是指的难度的绝对大小，而是如果这个是你所处阶段的基本知识课本上就有的重点，或者基础的个人作业等等任务。常有人一篇作业直接拍照传上来，啥也不说，这种行为非常恶劣。别人不是你的保姆。具体哪道题不会再问哪道题，一篇卷子都不会你应该做的是回去好好看书而不是问问问。还有“你就告诉我吧，我真不会啊，我太难了”之类的话，比你难的人有的是，好多人读了研究生问的问题还不如普通高中竞赛生，你自己丢不丢人。大家都是一路克服困难学上来的，除非你有身体障碍，我表示同情，否则不是你偷懒的理由。有很多资源贴，发贴的人的时间也是有限的。从中获得了资源后，建议对提供者表示口头上的感谢。由于某度抽贴的行为严重缺乏素质，含有url的贴经常被抽，或者资源也可能失效。这使得资源发布者只好用截图或者各种符号插入来屏蔽检索。提取出其中的有效信息只是举手之劳。请求重发之前应检查贴中各楼层的回复，观察楼主或者其他人有没有补lj，不要一lj失效了就问还有没有。尤其大量的lj并没有过期，比如图片发布的，经常有人分不清大写字母O和小写字母o以及数字0，还有大写字母I小写字母l以及数字1。至于我个人通常会选择尽量方便辨认的程序员字体。所以遇到lj不存在应该优先考虑double check你是不是打错了，而不是失效，通常失效lj不会被显示lj不存在。别人搜索发布资源不容易，没时间也没义务没完么了给你们重发，你只需要几分钟把容易打错的字符挨个试一遍就行，不要给别人添麻烦。资源发布者有时候面对明明有效的lj一天到晚有人求重发真的很崩溃。这类行为会被视为对楼主的严重打扰，必要时候会采取删封措施。对这个学科有基本的尊重，不要抱着戏谑的语气和态度来，例如“我今天铝热炸到天花板了真铥吓死我了哈哈哈”之类，这类人你们不配做化学实验，也趁早别接触这个学科。危险实验贴，无意义实验贴等会被毫不犹豫处理。不知道的内容不要强答，误导是非常恶劣的行为。这类内容同样会严肃处理。有后续内容会补充。感谢大家的配合，希望大家共同营造一个和谐的论坛环境。

一个芳香性反直觉的体系：1，4-二氢吡啶首先声明，下面提到的概念可能较多，如果有兴趣欢迎讨论，将在讨论中解释，这里一个一个解释打字实在太多了。之前有人发贴问过这个分子，一开始没多想。刚看到这个分子，有人可能会说大概没有芳香性，因为共轭体系没有闭环。也有人可能说吡啶的3号位和5号位的碳虽然没有直接连接，但是离得不是太远，而加一起是6个π电子的体系，符合所谓的Huckel规则，也许有所谓同芳香性。这两种观点在没有计算结果的时候，看起来好像都符合一些经验规则。那么真实结果是不是这样呢？我们用计算来说明。首先在B3LYP-D3(BJ)/6-31G*下优化分子，结果优化后的几何构型如下：看到这个构型，我们就会觉得哪里不对劲了。首先六个原子基本在一平面上，但是氮上的氢明显偏离平面，更像氨气的构型而不是吡咯的构型。这也暗示了后面的结论。考虑中间环的归一化多中心键级，这是将多中心键级开环上原子数次方的结果，用来消除对环原子数的依赖性，为0.44，一般这个数大于0.6才会被认为有芳香性，而小于0.5一般认为完全没有芳香性或者反芳香性。抗磁化率的各向异性Δχ定义为抗磁化率χ的矩阵绝对值从大到小的三个本征值，用绝对值最大的减去另两个的平均值。这里结果是-47.6908，这个数越负芳香性越强，一般芳香性体系都有-100以上，这个数说明这个体系应该不是芳香性体系。下面具体研究到底是非芳香性还是反芳香性。取平面上下的NICS(1)_ZZ，这个值是在距离环平面中心，正交于环平面1Å处各向同性对外磁场的化学屏蔽值的垂直平面分量的相反数，是最有效、最方便、最普适用，最常用的芳香性和反芳香性衡量指标之一。然后对环两侧的值取平均，得到9.45，是个比较大的正数，这是典型的反芳香性的特征。像苯这样的芳香性体系，NICS都是负几十的水平的负数。接下来用AICD方法分析环电流了，这是根据屏蔽张量矩阵元定义了一种实空间函数，数值越大处电流密度越大。电流密度矢量分布如下：AICD离域很差，内圈还是反左手定则的，这都是没有芳香性甚至反芳香性的标准。最后我们计算分子周围各点的化学屏蔽值，绘制等化学屏蔽表面，结果如下：绿色是正值，蓝色是负值。这个分布和苯完全相反，而和环丁二烯非常相似，这进一步说明了这个体系是明显反芳香性。因此有些震惊的结论是，这个体系是反芳香性的，电子结构某种意义上更像一个氨加上一个环丁二烯。而不是少了一个亚甲基后的吡咯。做实验的跟我说他们观测到这个结构的取代物不稳定，容易被氧化，一定程度上也说明了这个结果。下面发一下苯和环丁二烯的AICD和ICSS分析，方便对比。

一个芳香性反直觉的体系：1，4-二氢吡啶这是之前有人提过的一个问题，这里稍微整理一下回复。首先声明，下面提到的概念可能较多，如果有兴趣欢迎讨论，将在讨论中解释，这里一个一个解释打字实在太多了。之前有人发贴问过这个分子，一开始没多想。刚看到这个分子，有人可能会说大概没有芳香性，因为共轭体系没有闭环。也有人可能说吡啶的3号位和5号位的碳虽然没有直接连接，但是离得不是太远，而加一起是6个π电子的体系，符合所谓的Huckel规则，也许有所谓同芳香性。这两种观点在没有计算结果的时候，看起来好像都符合一些经验规则。那么真实结果是不是这样呢？我们用计算来说明。首先在B3LYP-D3(BJ)/6-31G*下优化分子，结果优化后的几何构型如下：看到这个构型，我们就会觉得哪里不对劲了。首先六个原子基本在一平面上，但是氮上的氢明显偏离平面，更像氨气的构型而不是吡咯的构型。这也暗示了后面的结论。考虑中间环的归一化多中心键级，这是将多中心键级开环上原子数次方的结果，用来消除对环原子数的依赖性，为0.44，一般这个数大于0.6才会被认为有芳香性，而小于0.5一般认为完全没有芳香性或者反芳香性。抗磁化率的各向异性Δχ定义为抗磁化率χ的矩阵绝对值从大到小的三个本征值，用绝对值最大的减去另两个的平均值。这里结果是-47.6908，这个数越负芳香性越强，一般芳香性体系都有-100以上，这个数说明这个体系应该不是芳香性体系。下面具体研究到底是非芳香性还是反芳香性。取平面上下的NICS(1)_ZZ，这个值是在距离环平面中心，正交于环平面1Å处各向同性对外磁场的化学屏蔽值的垂直平面分量的相反数，是最有效、最方便、最普适用，最常用的芳香性和反芳香性衡量指标之一。然后对环两侧的值取平均，得到9.45，是个比较大的正数，这是典型的反芳香性的特征。像苯这样的芳香性体系，NICS都是负几十的水平的负数。接下来用AICD方法分析环电流了，这是根据屏蔽张量矩阵元定义了一种实空间函数，数值越大处电流密度越大。电流密度矢量分布如下：AICD离域很差，内圈还是反左手定则的，这都是没有芳香性甚至反芳香性的标准。最后我们计算分子周围各点的化学屏蔽值，绘制等化学屏蔽表面，结果如下：绿色是正值，蓝色是负值。这个分布和苯完全相反，而和环丁二烯非常相似，这进一步说明了这个体系是明显反芳香性。因此有些震惊的结论是，这个体系是反芳香性的，电子结构某种意义上更像一个氨加上一个环丁二烯。而不是少了一个亚甲基后的吡咯。做实验的跟我说他们观测到这个结构的取代物不稳定，容易被氧化，一定程度上也说明了这个结果。下面发一下苯和环丁二烯的AICD和ICSS分析，方便对比。

关于卤键的介绍卤键是一类通常连在大电负性基团上的卤原子和强电负性原子或烯烃等其他局部富电子体系接近时候的一种弱相互作用，有饱和性和方向性。一个典型的体系是三氟一溴甲烷的溴对着氨气的氮的复合物结构。如下：关于卤键的成因和本质现在已经被研究清楚，主要是静电作用。具体地说，是卤原子连在大电负性原子上，电子向成键方向明显转移，从而导致成键方向的背面，卤原子电子密度下降，形成局部正电区域（叫做σ-穴或σ-hole），与另一分子的负电密集区域的静电相互吸引。本文将验证这件事。使用的软件有Gaussian，Multiwfn，VMD和PSI4。在B3LYP-D3(BJ)/ma-Def2-TZVP(-f)水平下对复合物结构优化，并用得到的波函数做下面的波函数分析。首先是静电势分析，按照Bader的定义取0.001的电子密度等值面作为分子的van def Waals表面，将表面上各点静电势通过颜色映射上去，越蓝越负，越红越正。按照优化后复合物的几何结构分别提取两个单体结构的分析结果如下：第一张图是溴代三氟甲烷的。最右边的小黄球代表了极大值点，即σ-hole的位置。可以看到溴原子外侧有明显的静电势为正的区域，且在Br-C键的反向位置。下面是氨气的静电势，最左边的蓝色小球代表了氮的孤对电子对应的静电势极小值点，有强烈的负的静电势。两者在复合物中相互吸引。我们保持两个单体“还没有相互作用”时候各自的静电势，叠加在一起，如下图：这是对静电势的分析。下面我们对弱相互作用做Yang等人提出的约化梯度密度分析。中间带点蓝色的白色区域就是弱相互作用区域，蓝色通常代表吸引作用，这里不太强。具体的内容就不多提了，这是将电子密度乘以电子密度Hessian矩阵的第二大本征值的符号，作为颜色映射到RDG等值面上。相互作用能是多少呢？用高精度的DSD-PBEP86-D3(BJ)/jul-cc-pVTZ结合CounterPoise计算的相互作用能为-17.13 kJ/mol。这里采用了完整的BSSE。用对称性匹配微扰法(SAPT)的SAPT2+/aug-cc-pVDZ即“银”方法对相互作用进行能量分解，结果如下。静电作用： -38.26196510 kJ/mol 交换互斥：42.93 kJ/mol 诱导作用：-9.68 kJ/mol 色散吸引：-12.33 kJ/mol 总相互作用：-17.34 kJ/mol。这个总相互作用和上面高精度计算的值是很接近的，一定程度上说明了该能量分解的合理性。可见，卤键的本质作用主要是静电吸引。由于卤键是σ-hole形成的，而分子动力学等方法用的原子电荷模型把一个原子周围的电荷全都归到这个原子上，相当于做了球对称化，而σ-hole完全不满足这种对称性，因此原子电荷完全无法体现这种卤键作用，从而分子动力学常规模型完全不能描述卤键。为了描述卤键，需要对模型进行一些校正，如在σ-hole的位置引入虚拟的拟合点，或者引入多极矩模型等。

关于苯环π-π堆积的本质自己算了一下，对苯环二聚体在B3LYP-gCP-D3(BJ)/def2-TZVP下优化。优化后的结构是两个苯环是错开的，几乎正好错开了半个苯环，这也说明了两者之间的相互作用主要根本不是什么轨道相互作用，否则应该正对着重叠。ωB97M-V/aug-cc-pVTZ结合CounterPoise下计算了结合能，结果是-12.2 kJ/mol。 SAPT的“铜”方法，sSAPT0/jun-cc-pVDZ下做能量分解，得到结合能是-12.1 kJ/mol，和上面计算的结合能-12.2 kJ/mol完全吻合。(这个SAPT级别有点儿低，再高内存不够，不排除结果这么准有点儿巧合)。能量分解结果如下：静电作用 -7.7 kJ/mol 交换互斥 27.4 kJ/mol 诱导作用 -2.9 kJ/mol 色散吸引 -28.9 kJ/mol 可见本例中pi-pi堆积的主要本质就是色散吸引作用。

【分享】使用Gaussian进行超小体系高精度热力学计算的程序 1.原理对于极小体系真空下的热力学计算，要考虑的体系热力学能量可以分为电子能量和热力学校正两部分。前者对应势能面极小值点的能量。电子能量在用Hartree-Fock和post Hartree-Fock方法计算的时候，可以分为自洽场能量即HF能量和电子相关能校正两部分。电子能量一个维度上受电子相关作用考虑程度的影响，一般来说越高级的电子相关post HF方法此维度上会越精确，但是计算量也会越大。另一个维度上受所用的基函数组的完备性的影响，越完备误差越小，但是也越昂贵。对于Dunning的相关一致性基组cc-pVxZ，随着zeta数增加基组越来越大，自洽场能量和电子相关校正能量都有很好的收敛性。有人事先拟合了一些公式和参数，按照这些方法可以在较小基组下获得较大基组的近似结果，这称为基组外推。 CCSD(T)方法对动态电子相关考虑已经很好，但是计算量太大。MP2方法考虑一般，但是计算量较小。利用电子相关能近似可加性，可以用MP2下小基组和大基组的电子相关能能量差来估计CCSD(T)下小基组和大基组的电子相关能能量差，这样可以用MP2下大基组的结果和CCSD(T)下小基组的结果(后者计算过程中会顺便得到MP2的结果)来估计CCSD(T)下大基组的电子相关能。使用2-zeta，3-zeta和4-zeta进行两者的外推是比较合理的。本文的脚本没有弥散函数，如果对于阴离子体系，可以使用月份基组加上一定程度的弥散函数。热力学校正部分能量在0K下为零点能，即最低振动能级和势能面极小值点的能量差。此温度下焓和吉布斯自由能的校正都是0。在室温下，总的吉布斯自由能在零点能的基础上额外还要加上吉布斯自由能随着升温的变化。这可以由Gibbs-Helmholtz方程用升温导致的焓的校正和熵的校正计算出来。由于非谐振模型计算极为耗时，常用的是谐振模型，这会引入一定系统偏差。为此有人提出了频率校正因子，分别对基频、零点能、升温的焓校正以及熵的校正做了拟合，将频率分别乘以对应的校正因子之后再用配分函数等算出的这类热力学校正能量更接近实验值。所有计算前必须进行几何优化，几何优化和频率计算必须在同一水平。本程序采用B3LYP-D3(BJ)/6-31G*进行优化和频率计算。这个精度对于这两个目的已经足够，而且其频率校正因子十分常见。电子能量部分分别计算CCSD(T)/cc-pVTZ、CCSD(T)/cc-pVDZ以及MP2/cc-pVQZ的能量。注意MP2的能量计算会顺便给出该基组的HF能量，CCSD(T)的能量会顺便给出该基组的MP2能量和HF能量。本方法最大的不足之一在于对大量涉及氢振动的低频模式，缺乏隧道效应的校正，会有精度的损失。但是由于热力学校正部分能量影响本来占的比例就较小，本文目的也只是学习为主，所以问题不大。其实Gaussian有自带的更好的高精度热力学组合方法使用。本文只是实现了一个低配的热力学组合方法。注意化学研究的是变化，只有相对数值有意义。直接用此程序算一个分子的热力学能量什么都不能说明。但是对于两个元素相同的反应，比如丙二烯到丙炔的异构化，用该方法分别计算两个分子的能量，然后相减，用得到的吉布斯自由能差通过van't Hoff等温式换算成标准平衡常数，是有意义的。

请教下关于windows下gvim终端版小键盘数字错码的问题使用的gvim，图形界面的gvim没有问题，cmd终端下的vim.exe用外接小键盘输入数字的时候会出现错码，像这样。还没有找到解决方案，请问大家有建议么，谢谢！

画了一组原子轨道... 用的球谐型函数，画图等值面数值0.036。需要注意的是，Schrödinger方程解的时候通常选取Hamilton算符、角动量大小平方算符和角动量z轴分量算符共同的本征函数，这样的轨道除了磁量子数为0的以外，都有复数部分，没法直接在三维空间作图。按照一定规则把磁量子数不为0的轨道线性组合，可以得到非角动量z轴分量算符本征函数的实函数，方便在三维空间作图。但是这样对f轨道和以上轨道实数化的处理方法是不唯一的，不同的量子化学软件处理方法也是不一样的。以下是按照Gaussian的处理方式。绘图使用Multiwfn和VMD。

画了一组s和p的杂化轨道 Gaussian+Multiwfn+VMD 等值面数值设的0.03，再大的话sp杂化的小端就看不见了。图形大小请只关注相对大小，因为绝对大小把显示用的等值面数值一改可能变化很大。 sp杂化：sp2杂化：sp3杂化：

关于双原子分子的s-p混杂以前本科上课时候的PPT，觉得挺好的，分享两张，不要随意传播。

关于CrO5(五氧化铬)中铬的氧化数这个问题总有人提到，+6和+10都有人说。这里稍微半定量分析一下。首先要明确，氧化数没有唯一定义，因为其是不可观测量，只是假想模型，因此原理上合理即可，也没法去测量。其相当于把所有化学键都当成离子键时候的各离子价态。但是无机化学讨论的时候又习惯用这种概念，因此这里简单分析一下。本文使用Head-Gordon等人提出的LOBA方法计算氧化态，其原理是将正则分子轨道做定域化，对定域化轨道分析其中的原子成分，若某轨道中某个原子的成分大于一约50%的阈值，就按照上面的假设，把这个轨道完全当做这个原子贡献的，最后按照轨道的电子布居数算各个原子的电子数，和核电荷数做差即得到该定义下的氧化态数。通常结果和轨道定域化化方法关系不大，本文采用Pipek-Mezey的轨道定域化方法基于Mulliken布居分析。此方法显然只与占据轨道有关。

【老生常谈】关于CrO5(五氧化铬)中铬的氧化数这个问题总有人提到，+6和+10都有人说。这里稍微半定量分析一下。首先要明确，氧化数没有唯一定义，因为其是不可观测量，只是假想模型，因此原理上合理即可，也没法去测量。其相当于把所有化学键都当成离子键时候的各离子价态。但是无机化学讨论的时候又习惯用这种概念，因此这里简单分析一下。本文使用Head-Gordon等人提出的LOBA方法计算氧化态，其原理是将正则分子轨道做定域化，对定域化轨道分析其中的原子成分，若某轨道中某个原子的成分大于一约50%的阈值，就按照上面的假设，把这个轨道完全当做这个原子贡献的，最后按照轨道的电子布居数算各个原子的电子数，和核电荷数做差即得到该定义下的氧化态数。通常结果和轨道定域化化方法关系不大，本文采用Pipek-Mezey的轨道定域化方法基于Mulliken布居分析。此方法显然只与占据轨道有关。

请警惕滥竽充数的某些培训班刚看到一个广告，一会儿删了。顺便转：辨别山寨坑钱科研培训的关键九点 http://tieba.baidu.com/mo/q/checkurl?url=http%3A%2F%2Fbbs.keinsci.com%2Fthread-3622-1-1.html&urlrefer=4e738800db117ebd70673183c7705ffb

【课外阅读】关于一氧化碳的偶极矩不再用水贴做标题了... 学过基础无机化学的同学们都知道，实验上一氧化碳分子偶极矩很弱，而且是（正指向负）是氧指向碳的，就是说氧这边带微弱正电，碳那边带微弱负电。之前有同学问过，为什么自己计算出来的结果不符合，这里简单计算一下。计算使用最常见的量子化学计算软件Gaussian，首先要清楚，Gaussian输出的偶极矩方向是【负指向正】的，是与常规定义相反的，因此结果必须取相反数，这是有些同学数据处理错误的重要原因。一氧化碳为C∞v点群分子，程序的标准取向会将分子摆在z轴上并让分子的质心在坐标原点。因此只要读输出的偶极矩在z轴上分量即可，另两个分量一定是0。再强调一点，有些同学用计算出来的Mulliken电荷布居来判断偶极矩的大小和方向，这是不妥的，加了弥散函数的Mulliken电荷布居意义非常差，而且Mulliken电荷布居只是一个很随意的经验模型，应该结合电荷密度分布直接计算偶极矩更为妥当。计算采用的方法分别有Hartree-Fock(HF)，Moller-Plesset的二阶微扰(MP2)和不考虑三重激发的四阶微扰(MP4(SDQ))，单激发和双重激发的耦合簇(CCSD)，普通杂化泛函中普遍认为计算偶极矩非常好的PBE0和Gaussian支持的两种双杂化泛函B2PLYP-D3(BJ)和DSD-PBEP86-D3(BJ)。基组分别采用Dunning的3-zeta相关一致性基组(CC-pVTZ)和对其每个角动量加弥散函数的版本(Aug-CC-pVTZ)。首先简单介绍一下，计算量大致是HF≈普通杂化泛函<MP2<双杂化泛函<CCSD<MP4(SDQ)，对于特定问题选用合适泛函的精度大致是HF<MP2<普通杂化泛函<双杂化泛函<CCSD≈MP4(SDTQ)（比MP4(SDQ)多考虑了三重激发）≈CCSD(T)，理论上金标准通常是用CCSD(T)，但是Gaussian里没有解析导数，不支持后HF电子密度计算和解析的偶极矩，因此很可惜无法采用。同样理由没有用MP4(SDTQ)。基组的话一般认为计算偶极矩等必须要加弥散函数。

【课外阅读】关于一氧化碳的偶极矩学过基础无机化学的同学们都知道，实验上一氧化碳分子偶极矩很弱，而且是（正指向负）是氧指向碳的，就是说氧这边带微弱正电，碳那边带微弱负电。之前有同学问过，为什么自己计算出来的结果不符合，这里简单计算一下。计算使用最常见的量子化学计算软件Gaussian，首先要清楚，Gaussian输出的偶极矩方向是【负指向正】的，是与常规定义相反的，因此结果必须取相反数，这是有些同学数据处理错误的重要原因。一氧化碳为C∞v点群分子，程序的标准取向会将分子摆在z轴上并让分子的质心在坐标原点。因此只要读输出的偶极矩在z轴上分量即可，另两个分量一定是0。再强调一点，有些同学用计算出来的Mulliken电荷布居来判断偶极矩的大小和方向，这是不妥的，加了弥散函数的Mulliken电荷布居意义非常差，而且Mulliken电荷布居只是一个很随意的经验模型，应该结合电荷密度分布直接计算偶极矩更为妥当。计算采用的方法分别有Hartree-Fock(HF)，Moller-Plesset的二阶微扰(MP2)和不考虑三重激发的四阶微扰(MP4(SDQ))，单激发和双重激发的耦合簇(CCSD)，普通杂化泛函中普遍认为计算偶极矩非常好的PBE0和Gaussian支持的两种双杂化泛函B2PLYP-D3(BJ)和DSD-PBEP86-D3(BJ)。基组分别采用Dunning的3-zeta相关一致性基组(CC-pVTZ)和对其每个角动量加弥散函数的版本(Aug-CC-pVTZ)。首先简单介绍一下，计算量大致是HF≈普通杂化泛函<MP2<双杂化泛函<CCSD<MP4(SDQ)，对于特定问题选用合适泛函的精度大致是HF<MP2<普通杂化泛函<双杂化泛函<CCSD≈MP4(SDTQ)（比MP4(SDQ)多考虑了三重激发）≈CCSD(T)，理论上金标准通常是用CCSD(T)，但是Gaussian里没有解析导数，不支持后HF电子密度计算和解析的偶极矩，因此很可惜无法采用。同样理由没有用MP4(SDTQ)。基组的话一般认为计算偶极矩等必须要加弥散函数。

【水贴】关于类氢原子非相对论B-O近似下的波函数程序 RT 无聊时候写出来了，python的版本比较好用也比较好写。C的版本之后放出来。 python使用了符号计算的sympy库。 R是径向部分，S是与z轴正方向夹角对应的函数部分，T是在xOy平面投影与x轴正半轴夹角对应函数部分，Y是全部角度部分函数，psi是整个波函数。有两个常量，Z是有几个核电荷，对于氢就是1，为了节省计算量可以把python程序中的Z=1注释删掉，注释掉原来定义的Z。另一个是a0，就是Bohr半径，原子单位下为1.，也可以按照上面的方法改掉节省计算量。 r是径向部分变量，theta是z轴正方向夹角，phi是在xOy平面投影与x轴正半轴夹角。这三者可以代入定义的符号，也可以代入数值进行某一点的计算。 n是主量子数，l是角量子数，m是磁量子数。没有写错误处理（后续大概会补上），但是需要注意n应该是正整数，l应该是非负整数且小于n，m应该是整数且绝对值不大于l。已经经过归一化。由于sympy库的原因，有些输出表达式不是最简，另外虽然原理上是精确计算不是数值求解，但是有时候根号表达式还是会丢掉变成小数。代码在楼下贴。

一个不考虑立体异构的烷基和烷烃同分异构体计数的算法用程序实现了一下，基于群论、圈分解、Burnside引理、Pólya计数定理以及生成函数。有时间补充原理。之前化吧有人发过，但是很不幸原稿似乎被百度删了，而且作者疑似没写完...这是运行结果。

这个广告呐啧啧啧啧

【讨论】不要强答针对某些用户，没人逼着你回复，自己不会的问题不要强答。最近胡答一气，基本事实都错误的贴子越来越多。五分钟连着删了五个。下次遇到有机会直接赠送黑名单。

请问这个是艾扎力么？一方去找芳川玩梗的路上出现的路人第五卷，一方初遇LO，LO吃饭晕倒，一方去找芳川路上的截图。时间线照理说是对的上的，那天就是因为早上假海原光贵勾搭炮姐还拆楼然后被当麻揍了，学园都市警戒等级提升。然后傍晚暗咲逢魔抓小黄书强行翻墙进学园都市，导致警戒等级进一步提升，天井亚雄才出不去的。这个看着有点儿像海原光贵，如果是被教主揍了后给组织打电话说不回去了之类的是不是说得通

新人请问下小说新约艾扎力（海原光贵）就没出场过么结标也是第一卷开头露个脸就再也没出场过了？

本周开始连着3周每周出1集！

请问cmd下有类似linux 的time命令那种计时命令么？请问cmd下有类似linux 的time命令那种给接下来跑的命令计时的命令吗？谢谢！

从半个月一集改成一个月一集了下一集4月4日...

华纳兄弟电视官网上又没更新第12季第16集还没出，还没看完网站写的要等多久，反正这周应该没了

字幕组辛苦啦真心感动，大过年的，2月8日晚上查的时候，最新一期已经被翻译出来了

请问一下已经添加的元素可以单个删除么请问一下已经添加的元素可以单个删除么，比如清空所有的annotation之类的，或者删掉某一条曲线啥的。。。感觉好像不能，毕竟没有这些元素的引用？谢谢

为什么cd if for这些命令没有manual手册。。。学习if和for命令全靠必应，求问下系统内置对if和for的说明在哪里，谢谢！

求助一句字符匹配有关的vim脚本小弟不才，在写一个vim script的函数，里面要进行一个判断，需要这样一个子函数：返回当前行的 "{" 比 "}" 多几个请问这个应该怎么写呢？另外，match()函数只能匹配左边第一个，如果我想匹配从右边开始第一个找到的，应该怎么写呢？谢谢大家！！

新人请问下二阶魔方XLL法的第一层 R.T. 第二层的42个公式已经背下来了，但是第一层总是找不到什么思路，总感觉做的太慢请教一下各位，第一层大概是什么思路，平均多少步比较好... p.s. 题主属于那种手残+空间想象能力差的，也就记忆力还算马马虎虎

新人想问下这个3*4*5的不等阶魔方 RT 只有U和D面可以有意义地转90°，剩下方向只能转180° 请问这个还原有什么思路么，三阶的各种最基本操作由于不能转90°，基本都不能用第一次怎么也拼不回去最后还弄散了，拆了重拼花了1小时。求个还原公式/方法，提供思路也好谢谢大家！ 2L放图

求问有可以算题的软件么就是输入一道题，比如死活题吃掉某一部分子，或者官子题之类的，而且可以根据用户的应对计算当前局部范围内的最优策略回应的，谢谢！

请问下这个单炮单象单士单兵胜单象双士这个怎么走啊，总是吃不掉象。我只会走炮士胜双士。谢谢大家，新年快乐！

求助一个vim的tags文件 RT 自己ctags -R生成的文件太大，有接近一个G。我不太清楚哪些该要哪些不该要，经常补全出一大堆却没一个想用的，而且都很慢。请问哪位大神有配置好的好用的tags文件，能不能发给我一下，多谢！邮箱[email protected]，或者其他方式传输也行谢谢大家了，祝新年快乐~ p.s.我现在用的是cygwin下的vim，默认没有tags。主要是希望bash脚本，C语言还有Python的补全，这些用的最多。另外有没有那种...英语字典或者英语单词表的补全，经常用英语写文章233

求问脚本的解释器路径带空格怎么办 cygwin下的模拟linux环境脚本文件test.py第一行 #!/cygdrive/c/Program Files/Python36/python 中间有个空格，执行起来会报错但是【加引号】和【转义】都试了都不行... 求助，多谢！

求助一个残局题库里的双马士象全胜单马士象全 R.T. 红先软件算杀一直没算出来，题目给的答案看不懂思路求个思路也好谢谢大家

那个，请问下英文怎么读 /目录 /root目录这两个英文怎么读啊。。。我一般叫“根目录”或者叫“杠”和“杠root"233

请问下vim如何判断是否是新文件我想在.vimrc里加一句，如果是已存在文件就不改变换行符类型，如果是新文件就默认用dos的换行符\r\n 就像这样 if * set fileformat=dos endif 请问这个if后面的条件怎么写，谢谢！刚想了下可以这么搞 if filereadable("*") == 0 set fileformat=dos endif 但是现在问题是我不知道怎么把文件本身作为参数传递进字符串，windows下的%1之类的好像行不通...就是命令行vim 1.txt这样后面那个参数 :file查询么？貌似也不好...

求问一下三阶魔方标准六面配色，精确的那种 RT 在自制贴纸图案，想要一下六个面颜色具体的RGB值（白色肯定是0xFFFFFF），一直查不到，多谢大家！