【独家|#杭州金钱豹逃亡一个月#：老猎人遇挫 救援队员“死要见尸”】三只未成年的金钱豹逃离动物园，令有着千万人口的杭州陷入某种慌张。这是一场疏忽、瞒报引发的系列风波。直到五一假期结束，村民居民连续报警，动物园才向官方坦白实情，早在4月19日，三只豹子已经丢失。直到今天，逃亡第二十九日，还有一只豹子生死未卜。豹子们不会想到，自己的出走会引发一场足以载入历史的大围捕。杭州几乎动用了全部力量，公安、特警、护林员、猎人，各路民间救援队，数千人进入山林。无人机、直升机和猎犬也纷纷加入战斗。正如一位搜捕队员所说，“反正是活要见豹，死要见尸。”
59岁的茶农祝财松第一个见到了豹子。5月2日，他独自在茶山上施肥，用手机拍下了豹子模糊的身影。饭桌上，这件事被他当作奇闻分享。然而，在他生活的龙门坎村，没人细究豹子来源，人们只是互相叮嘱上山要注意安全。直到5月6日早上，同村的萍姐才第一个发出报警讯息。显然，对当地公安来说，抓豹子也是项陌生工作。萍姐回忆，做笔录的民警一时想不起“豹”字怎么写，只好在纸上记下：“包子”。
祝财松没再见到豹子。但5月7日，二十公里外，别墅区金苑山庄的看家犬小黑看到了，还和它纠缠了一番。金苑山庄和祝财松所在的龙门坎村位于午潮山森林的两端，这里有山林，小溪，大白天可能偶遇在小区闲逛的黄鼠狼，碰到蛇也不稀奇。但顶级捕食者的出现还是让人恐慌。狗和豹子打架的时候，女主人正在家。她拨打了电话报警。
杭州已经许多年没有野生豹子了。人们很容易想到，豹子可能来自杭州市野生动物世界（以下称动物园），这里同处午潮山森林脚，与金苑山庄只隔了一个山头。但多家媒体找上门时，却被园方否认。动物园甚至对官方撒了谎。直到第二次报警后的5月7日深夜，在杭州森保中心专业人员追问下，动物园管理者才坦白：该园逃出三只豹子，已经自行捕回一只。
5月8日白天，民兵来了，特警来了。多路民间救援队也带着成群猎犬赶到动物园门口。再后来，直升机开始在金苑山庄上空低低地盘旋，隆隆作响。很快，下午5时许，传来第二只豹子落网的消息，就在金苑小区背后的山林。媒体立即和现场的狼群救援队龙羊猎人分队队长通上电话。这位队长开口就震惊了所有人，“（豹子）已经……除掉了。”
“已经死掉了是吗？”有记者问。“对对对，危害到人的生命安全了。”“是在搜救过程中导致的吗？”队长支吾着，没有回复。没一会儿，直升机的轰鸣声中，手机里弹出新闻——“第二只外逃豹子已经找到，搜救队长证实：被狗咬死了”。原本分散的记者在慌乱和诧异中聚到一起，有人反复确认通话录音，“没人说被狗咬死啊”。
现场始终处于混乱的状况。突然，一个记者向着西边小跑，紧接着两个、三个，一群记者跟在后面奔跑起来，“怎么啦？”没人知道发生了什么。领头者发觉后，转过身说，他只是去拿盒饭。
就当人们相信豹子死亡是事实时，媒体记者们得到了官方辟谣，第二只豹子并没死，而是麻醉捕获。广场上的记者再次奔跑起来——在动物园服务中心，官方安排狼群救援队队员邵琦接受采访。“被麻醉的豹子是通过猎犬找到的。狗没有咬，是在叫。豹子没有受伤，狗受伤了。”邵琦介绍追捕情况。“狗脸上的血是怎么回事儿？”“就是豹子抓的。”“为什么用烈性犬来搜救？””我无权回答。”不到两分钟，邵琦被工作人员拉走了。
但在网络上，关于第二只豹子的讨论远没有结束。人们从央视拍摄画面发现，第二只豹子疑似缺失一只后掌。专家在现场说，这是趾甲丢失。又有专业人士科普，豹子行动时脚趾着地，趾甲丢失并非皮外伤。杭州富阳区宣传部工作人员则称不存在脚断的情况。直到现在，没看到更清晰特写，人们不清楚豹子伤情究竟如何。也无人解释，豹子为何受伤。更不会有人知道，失踪的大半个月，它经历了什么。
那些曾经看到过野猴、狐狸的人们则开始怀疑，这些动物同样逃自动物园。动物园一路之隔的住户说，自己捡到过跑出来的孔雀、猴子，还了回去。还有参与过动物园建设的人记得，这里跑出过马鹿。人们也想知道，为什么豹子能轻易逃出动物园。闭园后，有记者从山中出发，依旧能轻易翻过围墙进入动物园，最矮的墙大概1.6米高，上面铁丝圈已经损坏。
搜捕一直在继续。很快，第三只豹子也有了消息。何家村附近的山林里，热成像捕捉到了豹子踪迹，5月9日凌晨1点，搜寻人员带着十多只搜救犬，采取合围的方式从外围慢慢靠近。但也许是察觉到什么，小豹子速度飞快地跑掉了。不过，一支救援队队长说：“我们队员和麻醉师上山确认。那只是两只野猫之恋，变长了身体骗过了热成像仪。”
何家村从未如此喧哗。警力将周边山林团团围住，最前线的山林搜寻，依旧主要由各个民间救援队负责。在这之前，这些救援队主要成就在于处置疯狗，摘除马蜂窝，抓野猪。老人走丢，有人溺水，也会找到他们。一支救援队的成员神情真挚地举例，某次为了找寻失踪老人，上山的救援队整整三天三夜没睡。但对所有人来说，找豹子还是头一回。
如果没有瞒报，搜寻显然不会如此困难。现在，没人知道豹子还在不在这座森林。一位搜救队员说，“跑到安徽去都有可能”。没有经验可以借鉴，人们只能自己想办法。随着时间推移，搜索工作越来越像大海捞针，失望的情绪正在弥漫开来。5月12日，经历四天搜寻后，一位救援队队员说，自己每天凌晨回家，身体已经疲惫，他觉得希望不大了。
最大的干扰，还是人太多。野生动物对气味敏感，以前打山鸡、野猪，猎人都单独行动，出发前不抽烟，不喝酒。“上山打猎必须安静，”一个老猎人说，“这回搜豹子这么多人一起上山，豹子肯定警觉，不好找的。”参与了两三天搜救之后，老猎人不再继续，接力棒交给了其他猎户。
在午潮山，这场战斗没有轻易结束。豹子一天没找到，有茶农一天不敢上山。原本计划在5月16日举行的杭州越野赛也延期了。一些救援队的参与人数明显变少，但仍在坚持，一位队员说，“反正是活要见豹，死要见尸”。
某天晚上，公狼救援队又上了山，十几名队员带着两只狗，一只比特猎犬，一只萨摩耶。手电筒的白光射向幽静漆黑的森林，寻找可能反光的动物眼睛，再用红外热成像仪比对。深夜的半山腰，只能听见人的喘气和说话声，“小豹子，回家咯”，有人冲着远方呼喊。爬了半个小时，萨摩耶率先崩溃了，赖在草堆里拼命喝矿泉水，说什么都不愿再向上挪一步。
三个小时的跋涉，终于传来令人兴奋的线索。 在刚下过雨没多久的山路上，救援队发现了一个清晰的、疑似猫科动物的脚印。消息传回，指挥部瞬间启动，几分钟里，公安、特警，拿着相机的人，带着动物园工作证的专家就坐车赶到。两架价值几十万的大型无人机在夜空升起，像星星一样闪着光，盘旋在高处，监测整座山林。
珍贵的脚印被人们围了起来，闪光灯不断，专家细看后皱了眉，“不像啊”，和豹子的肉垫相比，脚印中间一块太小。但还是不敢怠慢，专家让动物园饲养员发来出生一个月小豹的小脚。救援队的人嘀咕，不会是狗的脚印吧？一旁的白色萨摩耶被拉到泥地里，有人抓着它的脚在泥里按了按，也不像。
就在这时，一架无人机突然失去了平衡，从高空中掉落，一只螺旋桨碎成了三片，像只巨大的蜘蛛，静静停在黄梅坞水库的桥上，熄火了。#洞见计划# https://t.cn/A6V4y52G
59岁的茶农祝财松第一个见到了豹子。5月2日，他独自在茶山上施肥，用手机拍下了豹子模糊的身影。饭桌上，这件事被他当作奇闻分享。然而，在他生活的龙门坎村，没人细究豹子来源，人们只是互相叮嘱上山要注意安全。直到5月6日早上，同村的萍姐才第一个发出报警讯息。显然，对当地公安来说，抓豹子也是项陌生工作。萍姐回忆，做笔录的民警一时想不起“豹”字怎么写，只好在纸上记下：“包子”。
祝财松没再见到豹子。但5月7日，二十公里外，别墅区金苑山庄的看家犬小黑看到了，还和它纠缠了一番。金苑山庄和祝财松所在的龙门坎村位于午潮山森林的两端，这里有山林，小溪，大白天可能偶遇在小区闲逛的黄鼠狼，碰到蛇也不稀奇。但顶级捕食者的出现还是让人恐慌。狗和豹子打架的时候，女主人正在家。她拨打了电话报警。
杭州已经许多年没有野生豹子了。人们很容易想到，豹子可能来自杭州市野生动物世界（以下称动物园），这里同处午潮山森林脚，与金苑山庄只隔了一个山头。但多家媒体找上门时，却被园方否认。动物园甚至对官方撒了谎。直到第二次报警后的5月7日深夜，在杭州森保中心专业人员追问下，动物园管理者才坦白：该园逃出三只豹子，已经自行捕回一只。
5月8日白天，民兵来了，特警来了。多路民间救援队也带着成群猎犬赶到动物园门口。再后来，直升机开始在金苑山庄上空低低地盘旋，隆隆作响。很快，下午5时许，传来第二只豹子落网的消息，就在金苑小区背后的山林。媒体立即和现场的狼群救援队龙羊猎人分队队长通上电话。这位队长开口就震惊了所有人，“（豹子）已经……除掉了。”
“已经死掉了是吗？”有记者问。“对对对，危害到人的生命安全了。”“是在搜救过程中导致的吗？”队长支吾着，没有回复。没一会儿，直升机的轰鸣声中，手机里弹出新闻——“第二只外逃豹子已经找到，搜救队长证实：被狗咬死了”。原本分散的记者在慌乱和诧异中聚到一起，有人反复确认通话录音，“没人说被狗咬死啊”。
现场始终处于混乱的状况。突然，一个记者向着西边小跑，紧接着两个、三个，一群记者跟在后面奔跑起来，“怎么啦？”没人知道发生了什么。领头者发觉后，转过身说，他只是去拿盒饭。
就当人们相信豹子死亡是事实时，媒体记者们得到了官方辟谣，第二只豹子并没死，而是麻醉捕获。广场上的记者再次奔跑起来——在动物园服务中心，官方安排狼群救援队队员邵琦接受采访。“被麻醉的豹子是通过猎犬找到的。狗没有咬，是在叫。豹子没有受伤，狗受伤了。”邵琦介绍追捕情况。“狗脸上的血是怎么回事儿？”“就是豹子抓的。”“为什么用烈性犬来搜救？””我无权回答。”不到两分钟，邵琦被工作人员拉走了。
但在网络上，关于第二只豹子的讨论远没有结束。人们从央视拍摄画面发现，第二只豹子疑似缺失一只后掌。专家在现场说，这是趾甲丢失。又有专业人士科普，豹子行动时脚趾着地，趾甲丢失并非皮外伤。杭州富阳区宣传部工作人员则称不存在脚断的情况。直到现在，没看到更清晰特写，人们不清楚豹子伤情究竟如何。也无人解释，豹子为何受伤。更不会有人知道，失踪的大半个月，它经历了什么。
那些曾经看到过野猴、狐狸的人们则开始怀疑，这些动物同样逃自动物园。动物园一路之隔的住户说，自己捡到过跑出来的孔雀、猴子，还了回去。还有参与过动物园建设的人记得，这里跑出过马鹿。人们也想知道，为什么豹子能轻易逃出动物园。闭园后，有记者从山中出发，依旧能轻易翻过围墙进入动物园，最矮的墙大概1.6米高，上面铁丝圈已经损坏。
搜捕一直在继续。很快，第三只豹子也有了消息。何家村附近的山林里，热成像捕捉到了豹子踪迹，5月9日凌晨1点，搜寻人员带着十多只搜救犬，采取合围的方式从外围慢慢靠近。但也许是察觉到什么，小豹子速度飞快地跑掉了。不过，一支救援队队长说：“我们队员和麻醉师上山确认。那只是两只野猫之恋，变长了身体骗过了热成像仪。”
何家村从未如此喧哗。警力将周边山林团团围住，最前线的山林搜寻，依旧主要由各个民间救援队负责。在这之前，这些救援队主要成就在于处置疯狗，摘除马蜂窝，抓野猪。老人走丢，有人溺水，也会找到他们。一支救援队的成员神情真挚地举例，某次为了找寻失踪老人，上山的救援队整整三天三夜没睡。但对所有人来说，找豹子还是头一回。
如果没有瞒报，搜寻显然不会如此困难。现在，没人知道豹子还在不在这座森林。一位搜救队员说，“跑到安徽去都有可能”。没有经验可以借鉴，人们只能自己想办法。随着时间推移，搜索工作越来越像大海捞针，失望的情绪正在弥漫开来。5月12日，经历四天搜寻后，一位救援队队员说，自己每天凌晨回家，身体已经疲惫，他觉得希望不大了。
最大的干扰，还是人太多。野生动物对气味敏感，以前打山鸡、野猪，猎人都单独行动，出发前不抽烟，不喝酒。“上山打猎必须安静，”一个老猎人说，“这回搜豹子这么多人一起上山，豹子肯定警觉，不好找的。”参与了两三天搜救之后，老猎人不再继续，接力棒交给了其他猎户。
在午潮山，这场战斗没有轻易结束。豹子一天没找到，有茶农一天不敢上山。原本计划在5月16日举行的杭州越野赛也延期了。一些救援队的参与人数明显变少，但仍在坚持，一位队员说，“反正是活要见豹，死要见尸”。
某天晚上，公狼救援队又上了山，十几名队员带着两只狗，一只比特猎犬，一只萨摩耶。手电筒的白光射向幽静漆黑的森林，寻找可能反光的动物眼睛，再用红外热成像仪比对。深夜的半山腰，只能听见人的喘气和说话声，“小豹子，回家咯”，有人冲着远方呼喊。爬了半个小时，萨摩耶率先崩溃了，赖在草堆里拼命喝矿泉水，说什么都不愿再向上挪一步。
三个小时的跋涉，终于传来令人兴奋的线索。 在刚下过雨没多久的山路上，救援队发现了一个清晰的、疑似猫科动物的脚印。消息传回，指挥部瞬间启动，几分钟里，公安、特警，拿着相机的人，带着动物园工作证的专家就坐车赶到。两架价值几十万的大型无人机在夜空升起，像星星一样闪着光，盘旋在高处，监测整座山林。
珍贵的脚印被人们围了起来，闪光灯不断，专家细看后皱了眉，“不像啊”，和豹子的肉垫相比，脚印中间一块太小。但还是不敢怠慢，专家让动物园饲养员发来出生一个月小豹的小脚。救援队的人嘀咕，不会是狗的脚印吧？一旁的白色萨摩耶被拉到泥地里，有人抓着它的脚在泥里按了按，也不像。
就在这时，一架无人机突然失去了平衡，从高空中掉落，一只螺旋桨碎成了三片，像只巨大的蜘蛛，静静停在黄梅坞水库的桥上，熄火了。#洞见计划# https://t.cn/A6V4y52G
希望咸鱼收物的妹妹们可以明白只是出物的地方。可以讨价还价，给拍实物图只是让确认有无瑕疵和具体状态，但不是买家秀让你们种草拔草的，想看裙子颜色款式好不好看满不满意请自行搜种草图。
浪费我看山河令直播的时间给你拍照录视频，最后来句觉得太白怕显黑，问题是详情的官方图就挺白的（有色差），而且我拍的实物整体是生成色（米黄），东西没问题就是被当成买家秀拔草了[拜拜]
浪费我看山河令直播的时间给你拍照录视频，最后来句觉得太白怕显黑，问题是详情的官方图就挺白的（有色差），而且我拍的实物整体是生成色（米黄），东西没问题就是被当成买家秀拔草了[拜拜]
结构优化和单点能计算的能量为何不同
结构优化OPT和单点能计算SP是第一性原理计算最最基础的计算任务，计算结束时都会输出能量。细心的朋友可能发现同一个结构、用相同泛函和精度，这两个能量并不相等。那么，到底该用哪一个？鉴于不少同学问我，就干脆写个帖子解释一下。为严肃起见，我专门做了一个金属银的体系作为测试，并附录我测试的数据。以下为方便起见，下面用到的例子都基于密度泛函方法，而且采用PBE泛函，精度相关的设置并不高，仅为示例。
首先把我用到的模型我提供下图，包含三个：A是元胞1×1×1，然后我沿着X方向扩大获得B与C，分别是2×1×1和3×1×1, 这些数字表示我沿着X、Y、Z扩胞时采用的重复单元的数量。之所以选择不同大小的超胞，是下面会看到这个计算误差与这个尺寸也是有关的。我的测试采用八核工作站，CASTEP软件，每个测试花费时间在2分钟以内，感兴趣的朋友可以自行对比重复。具体操作就是先做结构优化（OPT），结束后把最后输出的能量收集起来（成为E1），然后用得到的结构重新做一次单点能计算（SP），收集能量并标为E2。我关注的是OPT和SP的能量不一致问题，所以研究ΔE=E2-E1。
测试结果放在下图，其中左图代表ΔE随模型大小的变化，而右图则是随计算精度的变化。可以总结出两个重要信息：（1）元胞1×1×1，ΔE很小，大约为10-4 eV,但一旦扩胞，ΔE迅速变大到0.1 eV量级，说明ΔE对体系大小非常敏感（参看左图，A与BC差1000倍）；（2）ΔE数值对优化精度不敏感（自洽精度从10-3 eV（low）提升10-5 eV（medium）到10-6 eV(high)，ΔE变化在10-5eV水平）。简单而言，OPT和SP肯定有差别，这个差别随体系尺寸变化会很显著。
从测试来看，显然OPT相比SP可能少了点什么。要理解这一差异，得明白第一性原理真正计算是如何运行的。这得追溯到多体粒子的薛定谔方程，物理化学方面一般指多电子体系，要解析求解太困难，所以一开始就采用了三大近似，分别为非相对论近似、波恩-奥本海默近似（又称绝热近似）以及轨道近似。按照波恩-奥本海默近似，原子核因体量太大总跟不上电子的运动，有点类似大象和蜜蜂，所以就把原子核和电子运动分离，同时要让二者关联，所以每次大象走一步都让会停下来让蜜蜂调整位置跟上，这样蜜蜂始终在大象周围，但坐标又不需要同步更新。这样处理肯定有误差，但大大简化了问题，直接效果就是求解电子状态（蜜蜂）时可以冻结原子核（大象）的坐标，所以一切与原子核相关的作用可以处理为常数。
这和今天的问题有什么关系呢？关系就在于结构优化OPT针对的是原子坐标（大象的位置）而没有精准定位电子的分布（蜜蜂），但是体系的能量是二者共同贡献的，也就是说，从能量角度看结构优化只是完成了一部分任务，蜜蜂完全可能在大象的鼻子或尾巴上，究竟在哪里体系能量最低呢？这需要在原子核位置确定下来后单独对电子状态进行优化，这正是单点能计算的意义所在。所以，当我们需要获得除结构参数以外的、和电子相关的准确信息时（例如能量、电子分布、电荷等）都需要基于优化良好的电子状态（即电子波函数）的基础上进行，否则是会存在人为误差的。这里需要提醒一些初学者，不要为了图简单或省时间，就直接使用结构优化的能量输出来做精确计算。这个误差会随体系电子数增加而增加，通过简单提高收敛精度并不能消除（参看上图）。
最后问一个更深入的问题：是不是每次SP的能量就一定比OPT要低？答案是不一定！这个不确定性来自于优化电子波函数（即SP步骤）的初始波函数是哪里来的。如果我们从前一步结构优化那里读取（前提是你的计算和软件都支持你输出波函数供后续计算分析，这个文件一般较大），那么我们一定是会得到更低的能量（哪怕只低一点点）。但很多初学者没有很强保留或者分析波函数的意识，而只读取了结构坐标，其本质是原子核的坐标，那么程序计算时往往是基于原子的基态波函数（一般就是原子轨道）出发来优化电子分布和具体波函数，完全有可能得到和几何优化有轻微差别的结果和能量本征值。其原因就是此时程序是看不到你前面已经获得的波函数，因此它进行自洽判断的依据是重新优化波函数前后步骤能量差，达到收敛标准就退出。这样就很可能得到的能量比几何优化输出的能量还略高。原则上讲，这并没有找到体系真正的电子基态。正是因为这个原因，一般我会建议做单点能计算都尽可能读取结构优化时输出的波函数和电荷分布，这不仅能大大加快计算进度，而且能获得更精准的基态波函数。
对量子力学和密度泛函理论不熟悉的朋友，可能对于三大近似没有概念，或难以理解波函数优化。建议阅读量子力学的相关教材，我也会在后面写个初浅的帖子，分别解释三大近似以及带来的问题，进而介绍一些更精确的修正方案。如果您觉得有用，请常回来，我会不时更新。这个帖子，放在了我的公众号和B站（请搜 计算老司机），请关注。本帖由计算老司机原创，转载请联系我，请勿用以商途。
结构优化OPT和单点能计算SP是第一性原理计算最最基础的计算任务，计算结束时都会输出能量。细心的朋友可能发现同一个结构、用相同泛函和精度，这两个能量并不相等。那么，到底该用哪一个？鉴于不少同学问我，就干脆写个帖子解释一下。为严肃起见，我专门做了一个金属银的体系作为测试，并附录我测试的数据。以下为方便起见，下面用到的例子都基于密度泛函方法，而且采用PBE泛函，精度相关的设置并不高，仅为示例。
首先把我用到的模型我提供下图，包含三个：A是元胞1×1×1，然后我沿着X方向扩大获得B与C，分别是2×1×1和3×1×1, 这些数字表示我沿着X、Y、Z扩胞时采用的重复单元的数量。之所以选择不同大小的超胞，是下面会看到这个计算误差与这个尺寸也是有关的。我的测试采用八核工作站，CASTEP软件，每个测试花费时间在2分钟以内，感兴趣的朋友可以自行对比重复。具体操作就是先做结构优化（OPT），结束后把最后输出的能量收集起来（成为E1），然后用得到的结构重新做一次单点能计算（SP），收集能量并标为E2。我关注的是OPT和SP的能量不一致问题，所以研究ΔE=E2-E1。
测试结果放在下图，其中左图代表ΔE随模型大小的变化，而右图则是随计算精度的变化。可以总结出两个重要信息：（1）元胞1×1×1，ΔE很小，大约为10-4 eV,但一旦扩胞，ΔE迅速变大到0.1 eV量级，说明ΔE对体系大小非常敏感（参看左图，A与BC差1000倍）；（2）ΔE数值对优化精度不敏感（自洽精度从10-3 eV（low）提升10-5 eV（medium）到10-6 eV(high)，ΔE变化在10-5eV水平）。简单而言，OPT和SP肯定有差别，这个差别随体系尺寸变化会很显著。
从测试来看，显然OPT相比SP可能少了点什么。要理解这一差异，得明白第一性原理真正计算是如何运行的。这得追溯到多体粒子的薛定谔方程，物理化学方面一般指多电子体系，要解析求解太困难，所以一开始就采用了三大近似，分别为非相对论近似、波恩-奥本海默近似（又称绝热近似）以及轨道近似。按照波恩-奥本海默近似，原子核因体量太大总跟不上电子的运动，有点类似大象和蜜蜂，所以就把原子核和电子运动分离，同时要让二者关联，所以每次大象走一步都让会停下来让蜜蜂调整位置跟上，这样蜜蜂始终在大象周围，但坐标又不需要同步更新。这样处理肯定有误差，但大大简化了问题，直接效果就是求解电子状态（蜜蜂）时可以冻结原子核（大象）的坐标，所以一切与原子核相关的作用可以处理为常数。
这和今天的问题有什么关系呢？关系就在于结构优化OPT针对的是原子坐标（大象的位置）而没有精准定位电子的分布（蜜蜂），但是体系的能量是二者共同贡献的，也就是说，从能量角度看结构优化只是完成了一部分任务，蜜蜂完全可能在大象的鼻子或尾巴上，究竟在哪里体系能量最低呢？这需要在原子核位置确定下来后单独对电子状态进行优化，这正是单点能计算的意义所在。所以，当我们需要获得除结构参数以外的、和电子相关的准确信息时（例如能量、电子分布、电荷等）都需要基于优化良好的电子状态（即电子波函数）的基础上进行，否则是会存在人为误差的。这里需要提醒一些初学者，不要为了图简单或省时间，就直接使用结构优化的能量输出来做精确计算。这个误差会随体系电子数增加而增加，通过简单提高收敛精度并不能消除（参看上图）。
最后问一个更深入的问题：是不是每次SP的能量就一定比OPT要低？答案是不一定！这个不确定性来自于优化电子波函数（即SP步骤）的初始波函数是哪里来的。如果我们从前一步结构优化那里读取（前提是你的计算和软件都支持你输出波函数供后续计算分析，这个文件一般较大），那么我们一定是会得到更低的能量（哪怕只低一点点）。但很多初学者没有很强保留或者分析波函数的意识，而只读取了结构坐标，其本质是原子核的坐标，那么程序计算时往往是基于原子的基态波函数（一般就是原子轨道）出发来优化电子分布和具体波函数，完全有可能得到和几何优化有轻微差别的结果和能量本征值。其原因就是此时程序是看不到你前面已经获得的波函数，因此它进行自洽判断的依据是重新优化波函数前后步骤能量差，达到收敛标准就退出。这样就很可能得到的能量比几何优化输出的能量还略高。原则上讲，这并没有找到体系真正的电子基态。正是因为这个原因，一般我会建议做单点能计算都尽可能读取结构优化时输出的波函数和电荷分布，这不仅能大大加快计算进度，而且能获得更精准的基态波函数。
对量子力学和密度泛函理论不熟悉的朋友，可能对于三大近似没有概念，或难以理解波函数优化。建议阅读量子力学的相关教材，我也会在后面写个初浅的帖子，分别解释三大近似以及带来的问题，进而介绍一些更精确的修正方案。如果您觉得有用，请常回来，我会不时更新。这个帖子，放在了我的公众号和B站（请搜 计算老司机），请关注。本帖由计算老司机原创，转载请联系我，请勿用以商途。
✋热门推荐