将Gaussian与ORCA联用搜索过渡态、产生IRC、做振动分析

将Gaussian与ORCA联用搜索过渡态、产生IRC、做振动分析 文/Sobereva @北京科音  2018-May-31   1 前言 在《将Gaussian与Grimme的xtb程序联用搜索过渡态、产生IRC、做振动分析》(http://sobereva.com/421)中笔者介绍了Gaussian的external功能的使用,以及如何利用这个功能,借用xtb程序产生的能量、受力、Hessian来进行过渡态搜索...

将Gaussian与Grimme的xtb程序联用搜索过渡态、产生IRC、做振动分析

注:本文的脚本支持的是2019-Feb-14更新的xtb,对更老版本不兼容,对之后的版本可能兼容也可能不兼容,发现不兼容时请自行修改脚本或者用之前的版本。 将Gaussian与Grimme的xtb程序联用搜索过渡态、产生IRC、做振动分析  文/Sobereva @北京科音 First release: 2018-May-29  Last update: 2019-Feb-15 0 前言 ...

2018年度计算化学公社杯最常用的量子化学程序和DFT泛函投票结果统计

2018年度计算化学公社杯最常用的量子化学程序和DFT泛函投票结果统计 文/Sobereva @北京科音  2018-May-27 在2018年4月16号,在著名的计算化学公社论坛(http://bbs.keinsci.com)开展了为期一个月的“你最常用的量子化学程序投票”(http://bbs.keinsci.com/thread-9713-1-1.html)和“你最常用的DFT泛函投票”(http://bbs.keinsci....

高性价比热力学组合方法G4(MP2)-6X的计算方法

高性价比热力学组合方法G4(MP2)-6X的计算方法 文/Sobereva @北京科音  2018-May-26 在2011年,J. Chem. Theory Comput., 7, 112-120中,几个非Gn系列方法的官方人员提出了G4(MP2)-6X热力学组合方法,号称耗时与G4(MP2)相仿佛,而精度很接近昂贵得多的G4。其相对于G4(MP2)最主要变化是把MP2换成了SCS-MP2,把CCSD(T)的CCSD和(T)部分的相...

利用布居分析判断基函数与原子轨道的对应关系

利用布居分析判断基函数与原子轨道的对应关系 文/Sobereva @北京科音  2018-May-26 0 前言 波函数分析程序Multiwfn (http://sobereva.com/multiwfn)的很多分析,比如计算分子轨道中某些原子轨道所占成分(见《谈谈轨道成份的计算方法》http://sobereva.com/131)、绘制PDOS图(见Multiwfn手册4.10节的例子)等等都需要在理解基函数与原子轨道的...

在大阪拍的“百合展2018”照片

PS: 本文图非常多,20多MB,请非常耐心等待加载... 在大阪拍的“百合展2018”照片 摄影/Chris  文/Sobereva  2018-May-14 之前看到一个百合展预告http://acg.178.com/201803/315211303529.html,在四月份大阪有个“百合展2018”,而且有很多知名作家会参展。虽然Sobereva非常想去,但由于诸事繁忙,抽身不得,于是托正巧要在日本一游的人去百合展把照片都拍了回来...

2018.4在大阪所购的二次元物品

PS: 本文图非常多,20MB左右,请非常耐心等待加载... 2018.4在大阪所购的二次元物品 文/Sobereva  2018-May-14 2018年4月托人去了趟大阪的百合展,展会照片见http://sobereva.com/417,在展会上搜刮了大量周边,在旁边的偶像大师灰姑娘女孩主题咖啡厅也搞了点限定周边。之后托人也顺便在animate购入不少物资,在这里展示一下。 百合展的宣传单 百合展限定T恤,耻度略大,据...

量子化学研究中切换泛函应当注意的问题

量子化学研究中切换泛函应当注意的问题 文/Sobereva @北京科音  2018-May-2 量子化学研究者经常会用一种做法,即计算一种问题用A泛函,计算另一种问题则切换为B泛函。笔者在线上和线下答疑的时候,都常看到有人不恰当地切换泛函,故此文专门说说这个问题。 一定要明白,如果你要切换泛函,必须必须必须有非常正当的理由。一般有这三种正当理由: (1) A泛函和B泛函适用的问题不同,因此在计算不同问题时用不同泛函。什么问题适合...

谈谈原子间是否成键的判断问题

谈谈原子间是否成键的判断问题 文/Sobereva @北京科音  2018-Apr-27 0 前言 经常有人在网上问类似这种问题:“我优化出的结构,在gview里怎么显示A原子和B原子没有成键啊,应该成键才对啊,我是不是算错啦?”这种问题每次都解释一遍比较麻烦,遂写此文专门谈一谈。 如果对于化学体系中的相互作用类型缺乏基本的了解,建议看一下《谈谈“计算时是否需要加DFT-D3色散校正?》(http://soberev...

谈谈“计算时是否需要加DFT-D3色散校正?”

谈谈“计算时是否需要加DFT-D3色散校正?” 文/Sobereva @北京科音 First release: 2018-Apr-21  Last update: 2019-Feb-14 0 前言 笔者经常被问及的一个问题就是“我要计算xxx,是否需要加DFT-D色散校正?”。关于DFT-D,在以下博文里已经有充分的介绍: 乱谈DFT-D http://sobereva.com/83 DFT-D色散校正的使用...

2018年4月完结动画简评

2018年4月完结动画简评 文/Sobereva  2018-Apr-15 本文是笔者完整观看过的2018年4月完结的19部全长或半长TV版动画和3部泡面番的简评。所有新番笔者都看过前几话,没追的要么是笔者对其题材有抗拒(比如恐怖或者胡搞类),要么就是制作质量太差、太低俗、太粗野、内容太老套、太平庸、太莫名其妙等作品自身问题。 通常评分区间是90~95分(低于90分一般不追,超过95的则是神作水准,极度罕见),不了解此系列文章评分传统的勿试图吐槽这点。评分依据...

使用Gaussian+PySOC在TDDFT下计算旋轨耦合矩阵元

2018-Dec-29注:有不少人表示PySOC运行不了。只要读者使用本文的操作系统、严格按本文的操作步骤计算,至少本文的例子是绝对能运行成功的。读者应当先确保本文例子能运行正常,由此确认安装方式和自己的操作过程无误后再去跑你自己的体系。如果你当前没有本文用的系统,可以用VMware虚拟机装,半小时就能装好。 2019-Feb-10注:笔者后来又写了《使用ORCA在TDDFT下计算旋轨耦合矩阵元和绘制旋轨耦合校正的UV-Vis光谱》(http://sobereva...

在VMD中将距离较近的分子质心连线的脚本

在VMD中将距离较近的分子质心连线的脚本 文/Sobereva @北京科音  2018-Mar-7 在计算化学公社论坛上有人问下图这种把距离较近的分子质心连线的图怎么绘制 其实这种问题用VMD的tcl脚本非常容易实现,下面是笔者写的绘制这种图的tcl脚本,20多行就实现了。可见VMD脚本稍微会一点就能解决很大问题。 draw delete all draw color orange set nres [llength [l...

root用户在用openmpi并行计算时避免加--allow-run-as-root的方法

root用户在用openmpi并行计算时避免加--allow-run-as-root的方法 文/Sobereva @北京科音   2018-Mar-4 OpenMPI从2.0开始有个变态要求,即如果用户是root,则通过mpirun来并行运行程序时,会提示 mpirun has detected an attempt to run as root. Running at root is *strongly* discouraged as any mistake ...

使用Multiwfn可视化分子孔洞并计算孔洞体积

使用Multiwfn可视化分子孔洞并计算孔洞体积 文/Sobereva @北京科音  First release: 2018-Feb-28  Last update: 2018-Aug-17     1 前言 前些天有个人在思想家公社群里问这种图怎么绘制,怎么计算孔洞体积 其实有很多程序和在线服务器都是专门用来显示孔洞、口袋并计算其体积的,但是这些程序绝大部分都是给生物分子用的,而对于小分子体系...