第十八届中国国际模型博览会的照片 Part2

提示:本文图片极多,一百多张,20多兆,请耐心等待缓冲完毕。另外一百多张图在本文Part 1部分,见http://sobereva.com/374 第十八届中国国际模型博览会的照片 Part2 文/摄影 Sobereva  2017-Apr-25 我颇喜欢舰娘里的俾斯麦 辽宁号 我超喜欢的大和号 中间这坦克的纹理我很喜欢,有数码感 演示3D打印 展厅外头的一部车,我挺喜欢 ...

第十八届中国国际模型博览会的照片 Part1

提示:本文图片极多,一百多张,20多兆,请耐心等待缓冲完毕。另外一百多张图在本文Part 2部分,见http://sobereva.com/375 第十八届中国国际模型博览会的照片 Part1 文/摄影 Sobereva  2017-Apr-25 第十八届中国国际模型博览会于2017年4月21日~23日在北京展览馆举办,洒家去了第一天,由于是周五工作日,人不算太多,估计主要是商务上的合作洽谈。本展主要是航模、车模、船模、无人机等,还有各种人物模型(主要是欧美、军...

谈谈赝势基组的选用

谈谈赝势基组的选用 文/Sobereva First release: 2017-Apr-18   Last update: 2017-Apr-19 1 前言 “算xxx元素该用什么赝势基组”是网上最常问的问题之一。本文说说该怎么选用赝势基组。看本文之前一定要非常仔细看《谈谈量子化学中基组的选择》(http://sobereva.com/336),本文假定读者已经充分...

Molcas的计算双原子分子光谱常数的模块vibrot使用简介

Molcas的计算双原子分子光谱常数的模块vibrot使用简介 文/Sobereva  2017-Apr-17 计算双原子分子的转动-振动谱涉及的光谱常数是量子化学研究双原子分子体系常涉及的。Molcas程序自带了一个不错的vibrot模块,专门用来计算这些量,使用时只需要直接提供离散的势能面扫描的数据即可,完全不依赖于Molcas的其它模块。另外它还可以计算其它信息,包括振动波函数、转-振态间的跃迁性质、温度平均下的可观测量(需要提供事先算好的不同键长时对应的...

使用Multiwfn考察轨道间重叠程度和质心距离

使用Multiwfn考察轨道间重叠程度和质心距离 文/Sobereva   2017-Apr-16 很早以前版本的Multiwfn的主功能100里的子功能11(简写为100-11)就能计算轨道间重叠程度。前一阵子有人问怎么衡量轨道间距离,虽然直接用Multiwfn的格点数据处理功能可以算,但是步骤略繁琐,笔者遂对100-11功能进行了扩展,现可以计算轨道间重叠程度、轨道质心位置和轨道间质心距离,原理和用法在此文简单说一下。这个功能对于分析很多问题比较有...

写计算化学文章时引用理论方法、基组、程序时应注意的问题

写计算化学文章时引用理论方法、基组、程序时应注意的问题 文/Sobereva  2017-Apr-15 写文章时恰当引用是非常重要的,这不仅是学术严谨性、对读者负责的体现,也是对原作者的尊重。如今,文章被他引的数目是衡量一个科研人员功绩最重要的指标之一,重要性远远大于文章篇数和IF加和(Google学术搜索、ResearchGate等站点的科研人员个人档案也都着重体现被引次数而不体现IF),因此只有当某个人的成果以正确的方式被引用,此引用才容易被有效地记入他的总...

用NBO计算原子轨道杂化后的能量变化

用NBO计算原子轨道杂化后的能量变化 文/Sobereva   2017-Apr-13 今天在思想家公社1群里有人问这种图里面的能量怎么得到 可以通过做价键计算来得到,但是比较麻烦。一种简单粗略的得到方法是用NBO。 用过NBO的人都知道,NBO分析可以直接给出NAO和NBO的能量。比如计算甲烷,有这样的输出    NAO  Atom  No  lang   T...

用NWChem做SODFT在DFT计算中考虑旋轨耦合效应

用NWChem做SODFT在DFT计算中考虑旋轨耦合效应 文/Sobereva  2017-Apr-12 1 前言 众所周知,对重原子构成的体系,尤其是从第五周期开始(H,He算第一周期,下同),需要考虑相对论效应,否则结果可能严重错误。相对论效应分为标量和旋轨耦合两部分。标量相对论效应好考虑,用相对论赝势时即可等效体现,多数主流量化程序如Gaussian也支持标量相对论哈密顿,结合专为标量相对论计算优化的全电子基组即可(见《在赝势下做波函数分析的...

使用DORI函数同时考察共价和非共价相互作用

使用DORI函数同时考察共价和非共价相互作用 文/Sobereva  2017-Apr-7 RDG(约化密度梯度)和ELF(电子定域化函数)大家都熟悉,前者适合用来考察非共价相互作用,见《使用Multiwfn图形化研究弱相互作用》(http://sobereva.com/68);后者可以考察共价作用,见《电子定域性的图形分析》(http://sobereva.com/63)。如果想把共价和非共价相互作用一起考察,一个显而易见的做法是把ELF和RDG结合使用,同时...

2017年4月完结动画简评

2017年4月完结动画简评 文/Sobereva  2017-Apr-6 本文是笔者完整观看过的2017年4月完结的18部全长TV版动画和4部泡面番的简评。注意其中多多少少有一些剧透的成份,怕被剧透又想了解哪些值得追的话就看评分就好了。所有新番笔者都看过前几话,没追的要么是笔者对其题材有抗拒(比如恐怖或者胡搞类),要么就是制作质量太差、太低俗、太粗野、内容太老套、太平庸、太莫名其妙等作品自身问题。 正在追的还没完结的有《虎面人W》、《BanG Dream!》和...

将GROMACS的原子电荷信息读入VMD的方法

将GROMACS的原子电荷信息读入VMD的方法 文/Sobereva  2017-Mar-27 VMD可以载入GROMACS的gro、trr、xtc文件,但不会载入速度、受力、电荷。以前写过怎么让VMD把GROMACS产生的速度和受力载入并绘制出来的博文: 《使VMD实时显示gromacs轨迹中原子的受力》(http://sobereva.com/36) 《使VMD读入Gromacs产生的trr轨迹中速度信息的方法》(http://sobereva.com/11...

通过柔性力常数考察键的强度

通过柔性力常数考察键的强度 文/Sobereva   2017-Mar-1 1 关于键强的衡量 键的强度是化学研究者十分关心的问题。但键的强度不是一个确切的可观测量,只是一个概念,具体衡量方式多种多样,常见的有 (1)键解离能(BDE) (2)键级、AIM理论定义的BCP上的性质 (3)键长 (4)键方向上的力常数 其中,BDE定义为H(片段1)+H(片段2)-H(整体),这里H指焓,整体和片段结构都是分别优化的。BDE意义很明确,而且实...

使用NWChem做分数占据数的DFT计算

使用NWChem做分数占据数的DFT计算 文/Sobereva   2017-Feb-25 1 前言 分数占据数(Fractional occupation number, FON)计算是指被计算的体系有的轨道占据数不为整数,也因此允许整个体系的电子数不为整数。这看似有点超现实,实际中体系电子数不可能为非整数,毕竟电子不可分割,但对于DFT(或HF)计算是可以实现的,而且在程序实现上超简单。大家知道,常规DFT计算的时候都是已知电子数,根据...

使用Multiwfn通过LOBA方法计算氧化态

使用Multiwfn通过LOBA方法计算氧化态 文/Sobereva  2017-Feb-22 1 原理 讨论化学体系中的原子,特别是配合物中的金属的时候经常用到氧化态(oxidation state)这个概念。这是完全人为的概念,不是可观测的。用氧化态这个概念时相当于假定所有键都是纯粹的离子键,电子在原子间转移量是精确的整数。显然,这个假定对绝大多数情况都是极为糟糕的。懂量子化学的人不爱用氧化态这个虚构的概念,因为原子电荷能明显能更真实客观地描述...

谈谈计算第一超极化率的双能级公式

谈谈计算第一超极化率的双能级公式 文/Sobereva  2017-Feb-21 SOS是计算(超)极化率常用方法之一,介绍见: 使用Multiwfn基于完全态求和(SOS)方法计算极化率和超极化率 http://sobereva.com/232 双能级公式是对SOS简化而来的。研究D-A体系的第一超极化率(β)问题时,研究者总喜欢用双能级公式来解释β值的内在本质,以及解释对体系做不同修改(如替换不同原子或基团)时β的变化趋势。双能级公式表示为: ...