简单量子化学问题答疑专帖（旧版，已关帖）

安衎然 · 发表于 Post on 2017-10-23 09:41:23

sobereva 发表于 2017-10-22 15:10
guess=read，读取真空下收敛的波函数再试

其它的做法

请问老师除了guess=read，还需要加geom=checkpoint关键词么？

lily19830130 · 发表于 Post on 2017-10-23 11:07:17

有些轨道的电子占据数是2，有些怎么不是整数？它们是怎么进行计算得到的？

阳光总在风雨后 · 发表于 Post on 2017-10-23 11:25:31

让你变成回忆发表于 2017-10-20 10:08
可以采用PCM模型，设置介电常数。

这是我的S0，S1，T1构型的输入文件，算出来的能量S0的比S1的还高，是不是有什么问题啊

让你变成回忆 · 发表于 Post on 2017-10-23 14:02:28

阳光总在风雨后发表于 2017-10-23 11:25
这是我的S0，S1，T1构型的输入文件，算出来的能量S0的比S1的还高，是不是有什么问题啊

先不看别的我觉得你的输入文件有些是多余的。
td计算的时候，默认算的激发态就是3个，其次优化的时候默认root=1，所以你td括号里面写的没任何意义。可以不写。
另外我不清楚你是否需要优化激发态来查看结构，我很少很遇到优化激发态的情况，太费时。一般就是考虑激发就行。
这是我个人的看法。

niubaotian85 · 发表于 Post on 2017-10-23 14:16:28

第一次使用MP2方法，请问，如果用MP2做单点能的时候，读能量，是读HF对应的数据还是MP2的呢

sobereva · 发表于 Post on 2017-10-23 14:34:15

安衎然发表于 2017-10-23 09:25
大家上午好，请问guess=read是默认上一个CHK文件吗？怎么指定高斯读取很早以前的文件啊?

%chk或%oldchk指向哪个chk，就从哪个chk读

sobereva · 发表于 Post on 2017-10-23 14:34:21

niubaotian85 发表于 2017-10-23 14:16
第一次使用MP2方法，请问，如果用MP2做单点能的时候，读能量，是读HF对应的数据还是MP2的呢

MP2

sobereva · 发表于 Post on 2017-10-23 14:34:51

安衎然发表于 2017-10-23 09:41
请问老师除了guess=read，还需要加geom=checkpoint关键词么？

输入文件里直接写了chk里的坐标就不需要geom=check，没写的话就得geom=check从chk里读取

sobereva · 发表于 Post on 2017-10-23 14:35:52

终_焉发表于 2017-10-22 22:35
Sob老师，我在新配置的电脑(WIN10 64位)上安装Gaussian 09 D.01跑任务，全都以
Two-electron integral sym ...

光看这些信息没法判断。
先算一个小体系，比如B3LYP/6-31G*算乙醇单点，能够正常跑完就说明高斯没问题。如果这个都跑不了，没准是软件冲突，机子里的安全程序阻碍程序运行之类。

sobereva · 发表于 Post on 2017-10-23 14:36:35

lily19830130 发表于 2017-10-23 11:07
有些轨道的电子占据数是2，有些怎么不是整数？它们是怎么进行计算得到的？

非整数的是自然轨道

sobereva · 发表于 Post on 2017-10-23 14:37:22

安衎然发表于 2017-10-22 15:45
谢谢老师，请问初猜为真空下收敛的波函数的相关命令怎么写呢？如何保证read的是指定的真空下收敛的波函数 ...

真空下单点算完了就有个chk文件，之后加了scrf的时候用guess=read，并且%chk还指向那个chk文件，就会从其中读取收敛的波函数作为初猜

sobereva · 发表于 Post on 2017-10-23 14:39:29

dingniu2 发表于 2017-10-22 22:08
请问老师，比较相互作用强度，用BCP处的电子密度有加和性吗？一个复合物A与B相互作用发生在3个位点，形成了 ...

1 不能，不是同类键
2 电子密度没法加和比较强度。即便是同类键，作用强度和BCP的密度虽然一般是正比的，但也不是线性关系。至于氢键相互作用强度，由于作用强度约为BCP处势能密度的一半，这时候多个BCP势能密度加和讨论总强度还说得过去。

sobereva · 发表于 Post on 2017-10-23 14:43:03

computational 发表于 2017-10-22 16:34
请问sob老师，我计算一个过多金属氧化物（含5个Fe原子）团簇与一些小分子的弱相互作用，用PBE/6-311++G（d ...

你的老师明显误导你
MP2算过渡金属明显不合适，远不如恰当选取的DFT泛函结果好，而且昂贵得多
什么泛函算过渡金属合适看
简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272

MP2/6-311+（3df，3pd）这种组合更是超级恶心
用什么基组看具体情况，见
谈谈赝势基组的选用
http://sobereva.com/373
谈谈量子化学中基组的选择
http://sobereva.com/336

一般来说，此体系建议B3LYP-D3或TPSSh-D3或M06-D3结合6-311+G**优化和振动分析，算相互作用能用那些泛函结合ma-TZVP或更贵的jun-cc-pVTZ（如果用得动双杂化泛函，对于高斯可以用B2PLYPD3）。

直接用PBE算弱相互作用结果很垃圾。除非加了D3校正

安衎然 · 发表于 Post on 2017-10-23 14:47:27

sobereva 发表于 2017-10-23 14:37
真空下单点算完了就有个chk文件，之后加了scrf的时候用guess=read，并且%chk还指向那个chk文件，就会从其 ...

谢谢耐心的老师，嘻嘻，我再试下~

k的研究生生活 · 发表于 Post on 2017-10-23 14:49:45

老师您好我想要用molpro计算beo分子的两个1重xigema态，输入卡部分如下
{MULTI;
                     occ,7,2,2,1;
                     wf,12,1,0;state,2;
                     wf,12,1,2;state,1;
                     wf,12,2,0;state,1;
                     wf,12,3,0;state,1;
                     wf,12,2,2;state,1;
                     wf,12,3,2;state,1;
                     }
                     ecas1(i)=energy(1);
                     ecas2(i)=energy(2);
                     ecas3(i)=energy(3);
                     ecas4(i)=energy(4);
                     ecas5(i)=energy(5);
                     ecas6(i)=energy(6);
                     ecas7(i)=energy(7);

                  {mrci;
                     occ,7,2,2,1;
                  wf,12,1,0;state,1;save,7111.1;maxiter,2000,2000;}
                     eci1(i)=energy(1);
                     eda1(i)=energd(1);
                     dipz1(i)=dmz(1);
                  {mrci;
                     occ,7,2,2,1;
                  wf,12,1,0;state,1,2;save,7121.1;maxiter,2000,2000;}
                     eci2(i)=energy(1);
                     eda2(i)=energd(1);
                     dipz2(i)=dmz(1);
可是总是勉强过了第一个点就报错，报错信息是problem occours in readm，我明明已经把两个态的ci计算分开了，为什么还是会出现readm问题

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[综合交流] 简单量子化学问题答疑专帖（旧版，已关帖）