💡简介 Orbital Viewer是一个主要由JavaScript与python写成的简易轨道查看器，用于提供直观的三维界面来展示和分析轨道，并快速整理计算结果，令使用者可以大致掌握分子的基本信息。 系统要求：目前支持Windows 10及以上版本，Windows 7未测试。Linux理论上可兼容，但作者没有具有图形界面的linux机子来测试。 网络要求：可以本地运行，未来可...

去年在量子化学公众号上看到了一种用Gaussian 16中的GIC功能实现同时扫描多个坐标的方法，当时没需求就收藏起来了。今天有个水传质子的过渡态猜测是协同的，想到了这个方法，于是研究了一下。 原理 假设需要计算一个乙烯与水的加成反应，需要扫描C4-O7’与C1-H9’同时缩短。 则定义冗余内坐标： RC4O7=R(4,7) RC1H9=R(1,9) 初始结构是随便摆的，RC4O7=1...

存点自用配置 gview风格配色 vmd可以按照如下格式自定义颜色，第四列是颜色ID，后面三列应该是RGB值/255 color change rgb 101 0.8000 0.8000 0.8000 color change rgb 102 0.8471 1.0000 1.0000 color change rgb 103 0.8000 0.4863 1.0000 color c...

态跟踪原理笔者在SDNTO中提过，参考了J. Comput. Chem. 2019, 40, 1420–1428。利用这个原理，可以对一条轨迹进行透热势能面跟踪，这对于追踪CT态与LE态的交叉十分有用。本文记录利用baneSNTO程序进行透热态跟踪的流程，选择的模型是对硝基苯胺。 ...

假期刷到初中的一个实验的视频，氯气和氢气在光照下化合。想起当时很不能理解明明室内开了灯，但是不能发生反应，非要镁条点一下才行。曾经问老师得到的解释是灯光强度不够，要Mg点燃的时候那个晃眼的强度才行。现在掌握的知识多了些，忽然有兴致随便算算玩玩。 首先分析一下。反应的关键是氯氯键需要被光打断，解离成氯自由基。这个过程涉及到氯气的激发和断键，应该是吸收到S1然后过MECI断键。既然涉及到吸收，...

记录一些origin偷懒小玩具 %Y：用户文件夹路径，可以把自己搓的小玩意放在这里，labtalk脚本里%Y调用 小玩具 raw_data 操作表格示例，把荧光仪测试的数据删掉前12行，第一列注释，第二列X，后续每隔两列删一列： import PyOrigin as po # 获取当前活动工作表 wks = po.ActiveLayer() # 获取工作表的列数 num_colu...

在计算光化学中，优化分子的激发态几何结构是理解光化学反应机理的关键步骤。然而，这个看似简单的任务实际上充满挑战。想象一下这样的场景：想要优化某个分子的S2，但在优化过程中，随着分子几何构型的变化，原本的S2态可能会变成S1，甚至S3。如果不加以追踪，优化算法很容易丢失目标，最终得到的不是你想要的那个电子态的最小值，而是另一个完全不同性质的态。这就是激发态追踪（state tracking）要...

BaneTask 2.1.0：把量子化学批量任务、条件分支与结果查询放进同一套工作流 BaneTask 2.1.0 是一个面向量子化学批量计算的任务控制与结果整理工具。它使用 .bt / .projbt 描述文件统一组织任务、依赖、脚本、结果快照与数据库，把一批分子、一套模板、一组后处理步骤沉淀为可复用、可查询、可比较的工作流。 它的设计目标很直接：让日常算化项目里的“铺输入、串上下游、...

已在公社发帖介绍，那里写的比较详细 Multiwfn以其强大的波函数分析功能而闻名。然而，数字交互的运行方式使得其输入文件可能较为难读。比如，笔者想导出前线轨道，则需要这一串抽象的数字作为输入： 200 3 h 2 3 l 2 0 q 虽然可以存为模板等待调用，但而假若下次想微调一些流程，就不能像QM计算的输入文件那样给opt加个calcfc这种方式来实现，这是笔者一直以来比较想解决的问...

理论背景 量子化学的核心目标是求解分子的薛定谔方程: [\hat{H}\psi = E\psi] 对于一个包含 $N(N>1)$ 个电子的分子，哈密顿算符 $\hat{H}$ 包含电子动能、电子-核吸引能、电子-电子排斥能。但这个多体问题太复杂，无法精确求解。 Hartree-Fock方法 为避免直接求解多电子波函数，一个经典的思路是Hartree-Fock方法。HF的核心思想是...