公社群聊Q&A整理：2016.06.16 02 ~ 2016.06.18 00 Concate - 计算化学公社

2016.06.16 01:13:37
Q:
老师您好请教一个问题，我想用Truhlar组的M08-HX这个泛函，但是g09里面没有，而他们的网站http://comp.chem.umn.edu/info/DFT.htm上只有MPW1K等的IOp设置似乎没有M08-HX的。请问我要怎样才能在gaussian或者qchem里面使用这个泛函呢？
A:
那就没辙了。
GAMESS-US里有
Q:
没有办法根据他们文献里面的描述进行IOp设置是吗？
A:
没法
----------------------------------------------------
2016.06.16 01:27:24
Q:
化学反应路径自动搜索软件GRRM（全称为Global reaction Route Map），有人用过这软件吗？有这种自动搜索路径的靠谱软件？
A:
寻找势能面上全部极小点、过渡态的算法GRRM简介
http://sobereva.com/31
Q:
您给的链接不是介绍GRRM的呀
A:
怎么不是
Q:
【链接】寻找势能面上全部极小点、过渡态的方法GRR
http://sobereva.com/319
少了个9
有了这个，做小分子反应的岂不是很幸福，啥都不用搞了
A:
做全维势能面扫描，做小分子反应的岂不是很幸福，极小点、过渡态和IRC都有了[图片]
----------------------------------------------------
2016.06.16 01:40:03
Q:
请问下谁有packmol安装包
谢谢
A:
http://www.ime.unicamp.br/~martinez/packmol/home.shtml
----------------------------------------------------
2016.06.16 03:05:26
Q:
md输入内容：
[图片]
A:
这么高温度分子本来就会跑散
如果你是指组成分子的原子跑散了，dt降到0.001，temp0降到0，tempi降到300再试
Q:
有两个甲基跑出了分子老远
A:
参数按我说的改，如果还这样，说明力场参数有问题
Q:
这次没跑散，也没报错，不过，这样的温度设定是否限制了其越过势能垒，导致总是在某个或某几个势阱徘徊？
A:
可以温度再调高点，逐渐试
温度太高，运动速度太快，2fs容易不稳定。现在1fs应该好多乐
Q:
变哪个温度？temp0还是tempi?
A:
temp0
Q:
谢谢sob老师，1fs确实好多了。tempi保持300，temp0逐渐升高，是这个意思吧，我试试
A:
y
----------------------------------------------------
2016.06.16 04:29:34
Q:
各位前辈好：我用b3lyp/6-31g**计算含有两个茂环和一个羰基的钛配合物，开壳层单重态的能量为-1349.8708，闭壳层单重态的能量为-1349.8704，三重态的能量为-1349.8703，这样开壳层单重态为基态，为其他两个态与基态能量几乎没差别，但是开壳层单重态的LUMO-HOMO为63.6kcal/mol，应该很小才对呀（应为gap小才会有近简并情况）。是不是这样体系不能用b3lyp优化？希望前辈给希望建议
而其他两个态与基态能量几乎没差别
但是他们的自旋确实不同，应该是三个钛
态
A:
B3LYP没问题，你可以看看开壳层单重态的自旋分布
开壳层单重态和gap小并没直接关系
Q:
sob老师，您说的是自旋密度吗？
A:
y
Q:
sob老师，如下是自旋密度 Mulliken atomic spin densities:
1
1 Ti -0.025787
2 O 0.010453
3 C 0.012264
4 C 0.005356
5 C -0.015162
6 C 0.012271
7 C 0.005378
8 C -0.015154
9 C 0.018626
10 H -0.000821
11 H 0.000653
12 H -0.000342
13 H -0.000820
14 H 0.000652
15 H -0.000343
16 C 0.006830
17 C -0.010958
18 C 0.006808
19 C -0.010974
20 H -0.000047
21 H 0.000581
22 H -0.000046
23 H 0.000582
A:
开壳层单重态本身也没法定义HOMO、LUMO。如果你说闭壳层单重态时候gap小那倒有点理
建议你也用BP86算一次看看
Q:
我算的是alfa轨道的GAP
A:
单重态双自由基gap大是很正常的
Q:
sob老师您好！闭壳层单重态的gap也有61.9kcal/mol，这样合理吗？
A:
也不能说不合理
本身你的体系双自由基特征也没那么强，自旋极化不是很厉害
----------------------------------------------------
2016.06.16 04:33:56
Q:
sob老师，我报名了培训班
A:
[表情]
----------------------------------------------------
2016.06.16 04:58:20
Q:
[图片][图片]？
？？？
A:
不用
----------------------------------------------------
2016.06.16 04:59:36
Q:
sob老师，请问用gaussian03计算溶剂化能DeltaG (solv)中是否要加命令scfvac？
A:
溶剂下算一次气相下算一次求差
Q:
不能一次性算完么？
A:
不要投机取巧
一次算完就自己写个--link1--
----------------------------------------------------
2016.06.16 05:07:51
Q:
sob老师，您有没有Guassian算例的教程，我刚接触，自己在网上也查了相关教程，可是很乱，不详细，不知道您有没有比较详细的教程之类的？@Sobereva
A:
最好学习途径就是参加7月份的我讲的量化培训班
其它方式学，花费的时间起码比这高一个数量级。网上所有相关教程有用的东西加起来也没培训讲的1/3多。
想自己找零零碎碎的教程，去群共享的公社专用邮箱索引里面搜
Q:
我是在上海太远了@Sobereva
A:
有不少从香港、广州那边过来的
Q:
老师暑假的培训班能不能提供电脑？我目前还没有买笔记本
A:
没法提供，可以租或借一个。不带笔记本也问题不大
Q:
租一个得多少钱呀？
有见过除了第一作者，后面全是星星的么
A:
我们不统一提供租赁服务
Q:
Sob老师，surface pro性能够不够
A:
对这没兴趣，很鄙视现在微软搞的一切东西
Q:
我是说培训班带这个够不够啊
I7版本的
A:
对性能没要求，程序能启动就够了
----------------------------------------------------
2016.06.16 05:26:15
Q:
见过全都是通讯作者的
A:
[图片]
----------------------------------------------------
2016.06.16 05:35:00
Q:
有没有现在比较流行的用于gaussian计算和分子动力学模拟(主要MS和AMBER)的小型工作站配置？
A:
2* XEON 2640v3
主板X10DRL-i
DDR4 64GB
2T企业级SATA
Q:
价格？
A:
2640v3，8核2.6，5200每颗
X10DRL-i（C612，支持v3，集显） 1950
三星DDR4 ECC 16GB*4 600*4
希捷企业级2TB 800
Q:
英特尔E5-2640V3（2.6G 8核）*2 64G内存 500G固态硬盘+2TB SATA企业盘 2U+导轨+600W电源
这个是别人给我的，您看呢？
A:
可以
Q:
[图片]你给的这个配置，价格大概多少
A:
[图片]
----------------------------------------------------
2016.06.16 05:57:03
Q:
求助大神Wiberg bond indices的原始文献！同时也求助各路大神，跪谢！！！
A:
[图片]
定义在脚注
[图片]
Q:
找了半天没找到，只能求助了
----------------------------------------------------
2016.06.16 06:18:33
Q:
! LANL08d EMSL Basis Set Exchange Library 6/15/16 7:15 AM
! Elements References
! -------- ----------
! Si P S Cl Ge As Se Br Sn Sb Te I Pb Bi: P.J. Hay and W.R. Wadt, J. Chem. Phys. 82, 284 (1985).
! Si P S Cl Ge As Se Br Sn Sb Te I Pb Bi: L.E. Roy, P.J. Hay, and R.L. Martin, J. Chem. Theory Comput. 4, 1029 (2008).
! Si P S Cl Ge As Se Br Sn Sb Te I Pb Bi: C.E. Check, T.O. Faust, J.M. Bailey, B.J. Wright, T.M. Gilbert and L.S. Sunderlin, J Phys. Chem. A 105, 8111 (2001).
!
****
I 0
S 1 1.00
0.7242000 1.0000000
S
1 1.00
0.4653000 1.0000000
S 1 1.00
0.1336000 1.0000000
P 1 1.00
1.2900000 1.0000000
P 1 1.00
0.3180000 1.0000000
P 1 1.00
0.1053000 1.0000000
P 1 1.00
0.0308000 1.0000000
D 1 1.00
0.2940000 1.0000000
****
! Elements References
! -------- ----------
! Na - Hg: P. J. Hay and W. R. Wadt, J. Chem. Phys. 82, 270 (1985).
! P. J. Hay and W. R. Wadt, J. Chem. Phys. 82, 284 (1985).
! P. J. Hay and W. R. Wadt, J. Chem. Phys. 82, 299
(1985).
!
I 0
I-ECP 3 46
f-ul potential
5
0 1.0715702 -0.0747621
1 44.1936028 -30.0811224
2 12.9367609 -75.3722721
2 3.1956412 -22.0563758
2 0.8589806 -1.6979585
s-ul potential
5
0 127.9202670 2.9380036
1 78.6211465 41.2471267
2 36.5146237 287.8680095
2 9.9065681 114.3758506
2 1.9420086 37.6547714
p-ul potential
5
0 13.0035304 2.2222630
1 76.0331404 39.4090831
2 24.1961684
! LANL08d EMSL Basis Set Exchange Library 6/15/16 7:15 AM
! Elements References
! -------- ----------
! Si P S Cl Ge As Se Br Sn Sb Te I Pb Bi: P.J. Hay and W.R. Wadt, J. Chem. Phys. 82, 284 (1985).
! Si P S Cl Ge As Se Br Sn Sb Te I Pb Bi: L.E. Roy, P.J. Hay, and R.L. Martin, J. Chem. Theory Comput. 4, 1029 (2008).
! Si P S Cl Ge As Se Br Sn Sb Te I Pb Bi: C.E. Check, T.O. Faust, J.M. Bailey, B.J. Wright, T.M. Gilbert and L.S. Sunderlin, J Phys. Chem. A 105, 8111 (2001).
!
****
I 0
S 1 1.00
0.7242000 1.0000000
S 1 1.00
0.4653000 1.0000000
S 1 1.00
0.1336000 1.0000000
P 1 1.00
1.2900000 1.0000000
P 1 1.00
0.3180000 1.0000000
P 1 1.00
0.1053000 1.0000000
P 1 1.00
0.0308000 1.0000000
D 1 1.00
0.2940000 1.0000000
****
! Elements References
! -------- ----------
! Na - Hg: P. J. Hay and W. R. Wadt, J. Chem. Phys. 82, 270 (1985).
! P. J. Hay and W. R. Wadt, J. Chem. Phys. 82, 284 (1985).
! P. J. Hay and W. R. Wadt, J. Chem. Phys. 82, 299 (1985).
!
I 0
I-ECP 3 46
f-ul potential
5
0 1.0715702 -0.0747621
1 44.1936028 -30.0811224
2 12.9367609 -75.3722721
2 3.1956412 -22.0563758
2 0.8589806 -1.6979585
s-ul potential
5
0 127.9202670 2.9380036
1 78.6211465 41.2471267
2 36.5146237 287.8680095
2 9.9065681 114.3758506
2 1.9420086 37.6547714
p-ul potential
5
0 13.0035304 2.2222630
1 76.0331404 39.4090831
2 24.1961684 177.4075002
2 6.4053433 77.9889462
2 1.5851786 25.7547641
d-ul potential
5
0 40.4278108 7.0524360
1 28.9084375 33.3041635
2 15.6268936 186.9453875
2 4.1442856 71.9688361
2 0.9377235 9.3630657
恕我愚钝，我还是没搞清I原子的这个基组是该用前一部分
****
I 0
S 1 1.00
0.7242000 1.0000000
S 1 1.00
0.4653000 1.0000000
S 1 1.00
0.1336000 1.0000000
P 1 1.00
1.2900000 1.0000000
P 1 1.00
0.3180000 1.0000000
P 1 1.00
0.1053000 1.0000000
P 1 1.00
0.0308000 1.0000000
D 1 1.00
0.2940000 1.0000000
****
还是后一部分
I 0
I-ECP 3 46
f-ul potential
5
0 1.0715702 -0.0747621
1 44.1936028 -30.0811224
2 12.9367609 -75.3722721
2 3.1956412 -22.0563758
2 0.8589806 -1.6979585
s-ul potential
5
0 127.9202670 2.9380036
1 78.6211465 41.2471267
2 36.5146237 287.8680095
2 9.9065681 114.3758506
2 1.9420086 37.6547714
p-ul potential
5
0 13.0035304 2.2222630
1 76.0331404 39.4090831
2 24.1961684 177.4075002
2 6.4053433 77.9889462
2 1.5851786 25.7547641
d-ul potential
5
0 40.4278108 7.0524360
1 28.9084375 33.3041635
2 15.6268936 186.9453875
2 4.1442856 71.9688361
2 0.9377235 9.3630657
这两部分又是什么关系？
177.4075002
2 6.4053433 77.9889462
2 1.5851786 25.7547641
d-ul potential
5
0 40.4278108 7.0524360
1 28.9084375 33.3041635
2 15.6268936 186.9453875
2 4.1442856 71.9688361
2 0.9377235 9.3630657
恕我愚钝，我还是没搞清I原子的这个基组是该用前一部分
****
I 0
S 1 1.00
0.7242000 1.0000000
S 1 1.00
0.4653000 1.0000000
S 1 1.00
0.1336000 1.0000000
P 1 1.
00
1.2900000 1.0000000
P 1 1.00
0.3180000 1.0000000
P 1 1.00
0.1053000 1.0000000
P 1 1.00
0.0308000 1.0000000
D 1 1.00
0.2940000 1.0000000
****
还是后一部分
I 0
I-ECP 3 46
f-ul potential
5
0 1.0715702 -0.0747621
1 44.1936028 -30.0811224
2 12.9367609 -75.3722721
2 3.1956412 -22.0563758
2 0.8589806 -1.6979585
s-ul potential
5
0 127.9202670 2.9380036
1 78.6211465 41.
2471267
2 36.5146237 287.8680095
2 9.9065681 114.3758506
2 1.9420086 37.6547714
p-ul potential
5
0 13.0035304 2.2222630
1 76.0331404 39.4090831
2 24.1961684 177.4075002
2 6.4053433 77.9889462
2 1.5851786 25.7547641
d-ul potential
5
0 40.4278108 7.0524360
1 28.9084375 33.3041635
2 15.6268936 186.9453875
2 4.1442856 71.9688361
2 0.9377235 9.3630657
这两部分又是什么关系？
A:
I 0
S 1 1.00
0.7242000 1.0000000
S 1 1.00
...
这一堆是赝势基组
I 0
I-ECP 3 46
f-ul potential
5
0 1.0715702 -0.0747621
1 44.1936028 -30.0811224
...
这一堆是赝势
两个你都得用
赝势基组必须结合赝势用
Q:
明白了
等我写完给你看
其他元素 0
其他元素的基组
****
I 0
I的基组
I 0
I的赝势
****
大概应该是这样，对吧
A:
其他元素 0
其他元素的基组
****
I 0
I的基组
****
I 0
I的赝势
Q:
坐标
空格
H C N 0
6-31g(d)
****
I 0
S 1 1.00
0.7242000 1.0000000
S 1 1.00
0.4653000 1.0000000
S 1 1.00
0.1336000 1.0000000
P 1 1.00
1.2900000 1.0000000
P 1 1.00
0.3180000 1.0000000
P 1 1.00
0.1053000 1.0000000
P 1 1.00
0.0308000 1.0000000
D 1 1.00
0.2940000 1.0000000
****
I 0
I-ECP 3 46
f-ul potential
5
0 1.0715702 -0.0747621
1 44.1936028 -30.0811224
2 12.9367609 -75.3722721
2 3.1956412 -22.0563758
2 0.8589806 -1.6979585
s-ul potential
5
0 127.9202670 2.9380036
1 78.6211465 41.2471267
2 36.5146237 287.8680095
2 9.9065681 114.3758506
2 1.9420086 37.6547714
p-ul potential
5
0 13.0035304 2.2222630
1 76.0331404 39.4090831
2 24.1961684 177.4075002
2 6.4053433 77.9889462
2 1.5851786 25.7547641
d-ul potential
5
0 40.42
78108 7.0524360
1 28.9084375 33.3041635
2 15.6268936 186.9453875
2 4.1442856 71.9688361
2 0.9377235 9.3630657
坐标
空格
H C N 0
6-31g(d)
****
I 0
S 1 1.00
0.7242000 1.0000000
S 1 1.00
0.4653000 1.0000000
S 1 1.00
0.1336000 1.0000000
P 1 1.00
1.2900000 1.0000000
P 1 1.00
0.3180000 1.0000000
P 1 1.00
0.1053000 1.0000000
P 1 1.00
0.0308000 1.0000000
D 1 1.00
0.2940000 1.0000000
****
I 0
I-ECP 3 46
f-ul potential
5
0 1.0715702 -0.0747621
1 44.1936028 -30.0811224
2 12.9367609 -75.3722721
2 3.1956412 -22.0563758
2 0.8589806 -1.6979585
s-ul potential
5
0 127.9202670 2.9380036
1 78.6211465 41.2471267
2 36.5146237 287.8680095
2 9.9065681 114.3758506
2 1.9420086 37.6547714
p-ul potential
5
0 13.0035304 2.2222630
1 76.0331404 39.4090831
2 24.1961684 177.4075002
2 6.4053433 77.9889462
2 1.5851786 25.7547641
d-ul potential
5
0 40.4278108 7.0524360
1 28.9084375 33.3041635
2 15.6268936 186.9453875
2 4.1442856 71.9688361
2 0.9377235 9.3630657
最下面俩空格，哦了不？
A:
对
Q:
突然有个想法，能混合基组，那能不能混合方法啊，呵呵
用QM/MM不就混合咯/偷笑
不分层
A:
混合方法必须分层
----------------------------------------------------
2016.06.16 16:47:15
Q:
[图片]请问下mulliken 、CHELPG、npa、Hirshfeld、ADCH电荷分别有什么区别呀？
A:
看http://www.whxb.pku.edu.cn/CN/abstract/abstract27818.shtml
----------------------------------------------------
2016.06.16 16:51:01
Q:
请问MP2对应DFT-D中的什么水平呢？@Sobereva
有没有M06-2X-D3好呢？
A:
没有
一般问题连B3LYP大多情况都不如
Q:
我之前的那个体系，MnCl3与管之间的弱相互作用为主的相互作用能，还是定性错误，我想看看后HF法和DFT-D哪一个好一点
基组已经加到了ma-def2-tzvp
A:
DFT-D
Q:
[表情]相较于之前的低水平的DFT_D+低水平基组，已经在小了不少，但是还是定性错误，我该继续增大水平和基组吗？还是考虑用MOPAC的半经验色散校正考察一下？
A:
没必要再增加
大抵是计算方式有问题
比ma-def2-tzvp再增大纯属浪费时间
Q:
计算方式？先低精度优化，高精度算单点，计算方式该怎么修改呢？
A:
电荷、自旋多重度之类问题
Q:
我修改电荷或者自旋吗？
A:
probably
----------------------------------------------------
2016.06.16 16:58:47
Q:
[图片]各位大神请看看我linux计算出来的log文件win下的gview打不开哇
A:
log打不开就打开fch
Q:
sob老师，我用fch打开的话我计算的预共振拉曼怎么看结果啊？
输出文件里可以看到，但是gview里找不到
A:
尝试改用multiwfn绘制
使用Multiwfn绘制红外、拉曼、UV-Vis、ECD和VCD光谱图
http://sobereva.com/224
Q:
其实multiwfn还不太熟，所以才准备参加培训
很简单，一个上午让你学会
A:
绘光谱图，一分钟学会
----------------------------------------------------
2016.06.16 17:21:54
Q:
[图片]
请问，实在找不见cu的极化函数怎么弄
是有固定值吗
A:
你直接到EMSL上复制Lanl2DZ(f)或Lanl2TZ(f)定义就完了
Q:
老师EMSL在哪里找到啊？
A:
在线基组和赝势数据库一览
http://sobereva.com/309
详解Gaussian中混合基组、自定义基组和赝势基组的输入
http://sobereva.com/60
----------------------------------------------------
2016.06.16 17:29:06
Q:
sob老师，请问SAOP potential是怎么计算得到？
A:
看J. Chem. Phys. 2000, 112, 1344; 2001, 114, 652
----------------------------------------------------
2016.06.16 17:30:02
Q:
我们有一些同学要参加，请问可以现在报名，7月15现场交现金吗
A:
建议汇款。现金的话培训期间估计给不了发票，得等到培训后再寄发票。如果觉得这样不碍事的话也可以，报名时注明付现金即可。
----------------------------------------------------
2016.06.16 17:30:49
Q:
M06-2X应引用哪篇lunwen？
论文
[图片]
TCA是吧
A:
Theor Chem Account (2008) 120:215–241
Q:
MBO键级和NPA电荷不引用也问题不大吗？
A:
应当引用
----------------------------------------------------
2016.06.16 17:30:49
Q:
请教各位老师，基态优化时，收敛标准是下面这样它就输出了，这个结果可靠吗？
.--------------------.
----------------------|Geometry convergence|---------------------
Item value Tolerance Converged
-----------------------------------------------------------------
Energy change 0.00000374 0.00000500 YES
RMS gradient 0.00016193 0.00010000 NO
MAX gradient 0.00084315 0.00030000 NO
RMS step 0.00011906 0.00200000 YES
MAX step 0.00059617 0.00400000 YES
....................................................
Max(Bonds) 0.0000 Max(Angles) 0.01
Max(Dihed) 0.03 Max(Improp) 0.00
-----------------------------------------------------------------
The step convergence is overachieved with
reasonable convergence on the gradient
Convergence will therefore be signaled now
[图片]
A:
单纯几何优化目的已经足够了
----------------------------------------------------
2016.06.16 17:44:48
Q:
Gaussian可不可算出分子某个原子的核外电子排布@Sobereva
A:
可
Q:
Gaussian可不可算出分子某个原子的核外电子排布@LiYuanhe
A:
做Mulliken、SCPA或者NPA分析
Q:
Mulliken charges默认就算了，但是这不是电子密度吗，原子的电子排布结果在哪呢@Sobereva
[图片]HNbN负离子的极小点是0.13 eV，并且出现在Nb原子附近；C2H4的极小点为9.26 eV在双键的位置。请问可以说明这两种物质发生反应时，Nb与双键相连接吗？
A:
pop=NPA就会输出啊
[图片]
Q:
NPA电子排布结果如这样，[core]5S(0.40)4d(9.98)5p(0.35)6p(0.01)，是什么意思@Sobereva
A:
说明不同原子轨道壳层有多少电子
字号不要超过12
----------------------------------------------------
2016.06.16 17:46:18
Q:
老师，我想咨询一下，我柔性扫面二面角，有两个能量极小点，那这两个结构有没有可能都存在呢
A:
可能
Q:
谢谢sob老师
都优化了，正在做频率分析
我的意思是在现实环境中会不会存在这种情况，怕投稿被质疑
A:
当然可能都存在
根据Boltzmann分布计算分子不同构象所占比例
http://sobereva.com/165
----------------------------------------------------
2016.06.16 17:54:06
Q:
MBO引用Mayer,I. Chem.Phys.Lett. 1983, 97,270
NPA引用Reed,A.E.;Weinstock,R.B.;Weinhold,F. J.Chem.Phys.
1985, 83, 735.
这样吗？
A:
对
对
对
Q:
请教sob老师如果我想计算活化能，上面那个只有三项收敛的优化结果，可以直接拿来用吗？
A:
问题不大
----------------------------------------------------
2016.06.16 18:26:13
Q:
sob老师，请教一个问题。我的一个结构用比b3lyp/6-31g(d)优化收敛，频率分析没有虚频；但拿结果用M062X/def2TZVP优化收敛，但结构变化很大，从原来的反式变到了顺式。结构优化结果是否以高精度的方法为标准？不是说几何优化对方法和基组的敏感性低吗？
A:
用B3LYP优化完的结构作为初猜，M062X优化时减小步长上限
----------------------------------------------------
2016.06.16 18:38:36
Q:
请问高斯计算分子尺寸结果中的a0
是分子半径还是直径啊
求大神解释之前说过但我忘了
A:
半径
----------------------------------------------------
2016.06.16 18:48:07
Q:
sob老师，请问配合物分子CH3COOCu算CHELPG 用HF/3-21g合适吗？
A:
这么点体系用3-21G做什么，起码TZVP
----------------------------------------------------
2016.06.16 18:51:01
Q:
[图片]请教一下各位老师，rwf文件设置中，设置为0的话是不是意味着rwf文件会保持一定的大小，一直正常计算，谢谢
A:
y
----------------------------------------------------
2016.06.16 19:20:48
Q:
请问高斯能计算反应速率吗
A:
使用KiSThelP结合Gaussian基于过渡态理论计算反应速率常数
http://sobereva.com/246
----------------------------------------------------
2016.06.16 19:40:35
Q:
咨询一下sob老师，在做背景电荷的时候，为什么用CHELPG电荷呢？原因是什么，或者有这样的参考文献？因为我怕后面写的时候审稿人问，就想先说一下为什么用CHELPG电荷。谢谢老师。
A:
看群共享《原子电荷计算方法的对比》，拟合静电势电荷对静电势重现性最好，能最充分表现静电相互作用
----------------------------------------------------
2016.06.16 19:42:19
Q:
请教个gaussian的问题，CCSD（T）和CCSD时，怎么改最大循环步数？我试了scf（maxcycle=200），但是还是每次到默认的最大循环50步就出错误信息，怎么改不了这个50步呢？
[图片]
input 文件
下面是output 文件[图片]
A:
在CCSD里面写maxcyc
你这控制的是SCF的，当然不生效
Q:
哦，那么这个output的cycle不是scf的cycle？
请问这是什么cycle？
一直以为这个是scf呢
A:
CCSD本身也是迭代求解的
完全不是一码事
Q:
那我能认为scf每次cycle是update波函数吗，这个ccsd的cycle是update什么？singlet excitation和doublet的excitation的那部分吗？
A:
[图片]
SCF迭代解的是LCAO系数，CCSD迭代解的是slater行列式的组合系数
----------------------------------------------------
2016.06.16 20:05:57
Q:
sob老师，β(-2ω;ω,ω)是SHG形式，那么β(-ω;ω,0)是什么形式呢。
A:
dc-Pockels形式
----------------------------------------------------
2016.06.16 20:22:52
Q:
北京科音自然科学研究中心将于2016年7月16日至19日于北京开办“2016北京科音量子化学基础暑期班”，以使得量子化学研究者以最好、最快的方式掌握量子化学应用性研究所需要的各方面理论知识和计算技巧，并熟练掌握主流量子化学程序的使用。此培训讲授者为卢天博士（网名Sobereva，北京科音自然科学研究中心主任、思想家公社群群主、计算化学公社论坛创始人、Multiwfn程序开发者）。本次培训详情、内容和报名方式参见http://www.keinsci.com/workshop/2016SKBQC.html。报名于7月12日截止，人数有限，欢迎踊跃报名并转发消息。
这次和5月份的培训讲的一样吗
A:
通知上有写
Q:
通知说了，是一样的，有些删减，有些改进
A:
本次培训内容与2016年5月举办的第二届北京科音基础量子化学培训班相同，但培训时间比5月份的培训缩短一天，将一些相对小众、重要性不是很高的部分略去不讲，而更突出重点。所用讲义内容和5月份的培训一致，但做了进一步的修订和完善。
Q:
是不是跟梁山一样，按登山先后排座次？
A:
随便坐
Q:
动力学的培训会不会直接连在这个后面呢？
A:
不会
Q:
我们学校出台的新政策原则上也不让学生飞机
A:
高铁、动车也挺好
Q:
原则归原则，实际上具体操作的时候写一个五六十字的说明就可以了
A:
原则就是用来违反的
Q:
第一次去帝都有点紧张啊
。。
帝都夏天干么？我怎么觉着桑拿天那个湿热
A:
时干时湿
Q:
sob老师，到时候培训完了撸串不？[表情]
A:
不
Q:
一起撸
Sob撸串的画面，想象不出来
A:
最后一天讲完了就赶紧回家歇着去了
Q:
撸点串，睡的香
边撸边讲
娃娃太胖，夏天怕热[表情]
[表情]
娃娃怕大伙把他撸吃了
我瘦了很多了
A:
娃娃已经学差不多了
Q:
对了sob老师，资料都会发么？
A:
会
Q:
电子的么
来兰州告诉我
A:
impossible
----------------------------------------------------
2016.06.16 20:59:11
Q:
hcth407这个泛函怎么样？
A:
不怎么样
Q:
哦哦
几乎不怎么用吗
A:
没什么人用，一堆参数，靠使劲拟合参数搞出来的。
Q:
写程序的话，可以用来写的比较好的泛函有哪些？
A:
水平够高都能写出来
Q:
xalpha写出来了，想加上别的
[表情]
我看了一下sob老师你写的关于密度泛函库的博文
A:
建议按顺序写出这些
SVWN->BLYP->PBE->B3LYP->PBE0
有余力再把M06-2X写出来
----------------------------------------------------
2016.06.16 21:04:16
Q:
--------Markdown Markdown--------
A:
--------Markdown Markdown--------
那个...
----------------------------------------------------
2016.06.16 21:10:29
Q:
<i|V_xc|j>，这个东西的计算比较耗时，需要先储存起来吗？貌似这个东西的量不小
A:
存，也不占什么内存
----------------------------------------------------
2016.06.16 21:21:46
Q:
老师我想请问下packmol里面比如我知道盒子体积3*3*8与浓度1mol/L 加入的钠离子和氯离子数目我知道怎么确定那如何确定加入的水分子数谢谢
A:
浓度知道了，离子和水的比例一算就知道了
----------------------------------------------------
2016.06.16 21:46:53
Q:
请问free rotor，harmonic rotor 和 hindered rotor有什么区别？
前两者不是一样地吗?
A:
[图片]
----------------------------------------------------
2016.06.16 22:37:42
Q:
各位，谁有关于高斯计算 exploring chemistry with electronic structure method 这本书里的练习的电子文件啊，网上只找到了这本书，没找到里面练习的com文件。
A:
第二版的
http://bbs.keinsci.com/forum.php?mod=viewthread&tid=96&highlight=exploring
第三版的http://expchem3.com
----------------------------------------------------
2016.06.16 23:00:48
Q:
同一输入文件，m06-2x/TZVP 与 B3LYP/6-311+g**优化出的过渡态有根本性不同，都有合理性……哪个更靠谱？
A:
前者
----------------------------------------------------
2016.06.16 23:55:13
Q:
老师：您好！文章中的这个图。[图片]，可以用multiwfn做吗？如果绘制的话，是用哪个部分呢？谢谢··
A:
可以，看手册4.4.6节的例子
Q:
恩老师，我正在绘制，可是绘制的图线比较细，而且[图片]，用multiwfn只能调等值线的粗细，可是试了一下，好像改不了粗细，只能是颜色
A:
保存出图片就能看到粗细效果了
----------------------------------------------------
2016.06.16 23:58:05
Q:
AIM的分子图？
说起来为什么很少有见做这个图的？
A:
因为你不搞这方面
AIM学习资料和重要文献合集
http://bbs.keinsci.com/forum.php?mod=viewthread&tid=362
这类图太多了
----------------------------------------------------
2016.06.16 23:59:33
Q:
咨询一下老师，这个ONIOM优化完之后，如何找这个high层的能量？谢谢老师。
ONIOM: gridpoint 1 method: low system: model energy: 0.243576548440
ONIOM: gridpoint 2 method: high system: model energy: -1302.687378502348
ONIOM: gridpoint 3 method: low system: real energy: 4.127407810362
是这个-1302.687378502348吗？谢谢
A:
model就是指high层区域，有low和high两个级别算的结果
----------------------------------------------------
2016.06.17 00:20:24
Q:
请教各位，[图片][图片]用的txt文件去展锋，但是展完后什么都没有，可以问一下我是哪一步操作错了吗？
A:
输入文件大抵有问题
Q:
[图片]
这样难道不对？
A:
你用-1看看是否正确载入了信息
Q:
好，谢谢老师，我试试
A:
之前忘了提供设定粗细的选项，刚在主页更新了，用最新的3.3.9(dev)版本你用这个选项在设定颜色后就会问你粗细设多少
Q:
老师，是没载入，我是直接将文件拖进去，然后有一个路径，回车进入主页面，我再试试看是哪儿出错了
A:
载入之后选
11
3
0
看到的应该是这样的
[图片]
----------------------------------------------------
2016.06.17 00:28:50
Q:
MPW1B95：31% ； PBE0-1/3：33.33%这两个泛函做激发常用吗？
A:
不常用
Q:
常用的pbe0,pbe38,cam-b3lyp,MO6-2X……是吗
ωB97XD
A:
是0不是O
是
----------------------------------------------------
2016.06.17 00:31:38
Q:
若单体和二聚体的某个原子的NPA电荷变化只有0.009，是否相当于没有发生电子转移？
A:
几乎没转
----------------------------------------------------
2016.06.17 00:36:15
Q:
一直有点疑惑~AIM中的闭壳层作用和我们平时说的闭壳层不是一会事吧？
A:
闭壳层作用就是指非共价作用
----------------------------------------------------
2016.06.17 00:36:52
Q:
[图片]，谢谢老师，我已经成功了，老师，比如我想把几种方法算的数据展在一张图上，我没看到manual上有相关介绍的，是不是得把得到的所有spectrum_curve.txt文件导入到origin去作图？
A:
是
----------------------------------------------------
2016.06.17 01:18:36
Q:
请问分子直径如何计算吧
A:
谈谈分子半径的计算和分子形状的描述
http://sobereva.com/190
----------------------------------------------------
2016.06.17 01:37:05
Q:
sob老师，请问酱紫的图是用什么软件画的呀
[图片]
就是不知道用了这样的画图软件，在paper里需不需要文献引用
请问是不是acknowledgement里面提一下就好了
A:
不用引
如果你买了chemcraft就引网站
----------------------------------------------------
2016.06.17 01:39:51
Q:
Sob老师请问开壳层体系是不是不能用ALIE分析呢？
A:
开壳层也能用ALIE
----------------------------------------------------
2016.06.17 02:19:18
Q:
[图片]请问这个下面的donor，acceptor以及overlap的图怎么做出来的？
A:
使用Multiwfn绘制NBO及相关轨道
http://sobereva.com/134
----------------------------------------------------
2016.06.17 02:24:03
A:
鉴于经常用人问基组选用的问题，专门写了个博文
http://sobereva.com/336
鉴于经常用人问基组选用的问题，专门写了个博文
http://sobereva.com/336
Q:
顶
里面全是各种高频问题
----------------------------------------------------
2016.06.17 02:25:36
Q:
[图片]
[表情]
我感觉做激发态计算的时候还是跟体系有关，我的研究体系，用TZVP作基组，泛函用了很多种，最后发现B3LYP/TZVP才是与实验值比较相符的组合
A:
泛函选择显然和体系有关，小分子和大共轭体系用的泛函明显不同
----------------------------------------------------
2016.06.17 02:38:29
Q:
bmp多少位色质量最高？
多少色位
A:
24bit
----------------------------------------------------
2016.06.17 02:38:29
Q:
老师，有没有关于能级轨道和电子特征的帖子？就是关于成键轨道和反键轨道
A:
看结构化学
----------------------------------------------------
2016.06.17 03:32:44
Q:
老师，gaussian计算得到的FCHK文件里是不是不包含外场信息呢？比如我现在加了电场，那么加了电场后的FCHK里就没有把相应的外电场信息包含？这样的话，想看加外电场后分子的一些信息，如静电势，就不能得到了？@Sobereva
A:
能得到啊
加了外场时得到的波函数会相应地改变，都体现在fch里了
Q:
老师，既然fchk可以反应外电场信息，那我可以先用加完电场的FCHK在MULTIWFN里画静电势的梯度图，然后用没加电场时的FCHK画静电势的梯度图，最后用前者减去后者，得到的可以说是外电场引起的分子的电场强度的变化吗？（还有，在MULTIWFN里用没加电场和加电场的FCHK去画分子平面电场图，二者是一样的，看不出有差别,这样如何解释呢？）老师，其实我现在就想看外电场下分子电场强度。@Sobereva
A:
说明加的电场太小，从图上看不出显著差异。
你可以先做个密度差图看看
----------------------------------------------------
2016.06.17 03:54:16
Q:
[图片]我也是醉了
又是培训广告
A:
这种培训，交了钱不给发票或给假发票，甚至去了都找不到人都有可能
Q:
[图片]
有些公司的培训竟然开了yy
[图片]
本来还庆幸没收到，结果lammps广告还是发到这里了
A:
人人喊打
Q:
有这个群，不容易上当
难道真是群里出内鬼了
那个培训不靠谱啊?
lammps
[图片]请问各位大神这个是咋回事啊?谢谢
那个培训不值得去啊?还准备和老板申请去呢?
A:
坑钱的培训，这两天群里都在批
Q:
价格就已经非常坑了
A:
这种山寨味儿这么重的培训一看就不能去
Q:
唉刚还把邮件转了
A:
赶紧通知他们千万别去
现在有一些骗子混入本群挨个给群成员qq邮箱发坑钱的培训、会议通知，一定不能信
Q:
怪不得垃圾增多
群主可以排除一下
后面加的可能大一些
A:
问题是不知道是谁
揪出来有赏
Q:
1210个人群主排除难度有点大啊
A:
1998个人，在线1211个
Q:
看加入时间
A:
从验证信息上看不出来
而且可能已经潜伏了很久了
Q:
后面加入的可能性大一些
我估计不会，他们性急
潜伏不是他们的特征
很狡猾的，故意用查不到的邮箱发
A:
没法查，没证据
Q:
肯定是有计划有组织的
写到群通知上就可以了吧，或者在群共享里面加个‘已认定垃圾培训清单’
A:
群通知字数很有限，还得放科音的培训通知
放群共享里没人看
只要大家相互转告，平时不用qq邮箱的人如果qq邮箱收到培训、会议的信息，一律不能去
Q:
几天前收到了这个
我也收到了
这个收到过好几次
也是骗子
A:
GRRM这个倒不是骗子
不过如果嫌老收到邮件烦可以拒收
Q:
老板刚给申请的所里的
这个没事吧
A:
[图片]
Q:
当然，发邮件前咨询一下群主比较好。
A:
我是极反感也不允许私自给群成员qq邮箱发广告邮件的行为的
即便科音办培训都没给每个群成员发邮件，免得被人嫌扰民
A:
谨防各种骗钱的科研、软件培训和学术会议
http://bbs.keinsci.com/forum.php?mod=viewthread&tid=3547
Q:
哈哈，看了社长义愤填膺的鞭笞骗子的文章。。。科音培训的口碑肯定会越来越好，美誉度肯定也会越来越高，时间会见证。
A:
tks
A:
鉴于某公司还在不停地往本群群成员的QQ邮箱发坑钱的培训邮件，已成公害，有必要再次发一下这个链接
谨防各种骗钱的科研、软件培训和学术会议
http://bbs.keinsci.com/forum.php?mod=viewthread&tid=3547
此公司的人极有可能潜伏在本群内，获取成员qq邮箱地址，用以大批量发骗钱培训的邮件。若能提供证据揪出此公司的人潜伏在本群内的QQ号，请私信管理员sobereva，核实后奖励100元！
顺带一提，北京科音自然科学研究中心举办的暑期量子化学培训将在7月16日于北京举行，本人主讲，欢迎参加，详见科音官网www.keinsci.com。
顺便再说一句，本群管理员从没授权过给任何第三方公司往成员QQ邮箱里发宣传邮件的权利！
Q:
[图片]我这邮箱是不是被黑了经常有这些培训的
A:
见http://bbs.keinsci.com/forum.php?mod=viewthread&tid=3547
一些不法分子将本群当成获取科研人员联络方式的渠道进行诈骗
Q:
估计不止本群....
A:
肯定的
本群人数最多、影响力最大，最容易被盯上
----------------------------------------------------
2016.06.17 03:54:16
Q:
http://fcode.cn/
这个网站咋样？
A:
网站不错
----------------------------------------------------
2016.06.17 03:54:16
Q:
[图片] sob老师，这本书主要涵盖了哪些方面的量化呢？是不是从头开始讲（基础班）？
A:
是
----------------------------------------------------
2016.06.17 04:32:26
Q:
[图片]老师，这个图我只想要几个关键的原子标签，有什么方法可以把其他的原子标签隐掉呢@Sobereva
[图片]
A:
可以
Q:
老师，选哪个选项啊
A:
[图片]
Q:
这样全部都隐掉了，老师
我只想隐掉部分的
A:
自行ps
后处理菜单有选项关掉标签
----------------------------------------------------
2016.06.17 05:16:49
Q:
请教一下，NBO分析结果中，标出的红色部分是什么意思呀[图片]
A:
看此文http://sobereva.com/97
[图片]
----------------------------------------------------
2016.06.17 05:21:39
Q:
请问杂志投稿的时候要求推荐3个审稿人，找不到怎么办
怎么找
[图片]
-你是说在文献上找3个大牛就可以了吗
就是说在文献上随便找3个就可以了吗
那理由怎么说呢
[图片]就名字和邮箱是必填的
请问sob老师的名字和邮箱在哪，有相关的文献可以搜索吗
谢谢
那sob老师的名字是什么
A:
如果文章用了Multiwfn，可能会收到审稿请求
如果文章没用Multiwfn，我收到审稿请求的可能性比较低
----------------------------------------------------
2016.06.17 05:28:25
Q:
老师，我想问一下计算NBO时，出现l607错误，怎么解决啊。如果在出现错误的时候已经生成47文件，生成的47文件是否正确呢
[图片]
输入文件是这样的
A:
看报错提示
Q:
[图片]
老师报错是这样的
A:
有47文件就够了，没问题
用了saveNBO导致的
Q:
老师，那就不用继续修改错误继续计算，是吗？
A:
y
----------------------------------------------------
2016.06.17 06:02:21
Q:
请教一下，轨道成分分析只能看出每个原子所占的比例，除了用chk 或者fchk自己判断具体属于dxy或者fxyz等等，还有什么方法能得到具体属于哪个轨道类型吗
A:
用Multiwfn做轨道成分分析，原子轨道的贡献也能得到
输出信息一目了然
----------------------------------------------------
2016.06.17 06:07:26
Q:
请问一下，用[图片]算出的核磁数据，为什么同一个C上的H的化学位移不一致呢？比如[图片]，C12原子上两个H的位移分别是[图片]，相差的特别多，请问是为什么呢
但是实验测出的峰虽然有分裂，但是化学位移不至于相差到1.6ppm
A:
计算级别严重不合适
B972/pcSseg-1是首选
用什么合适这里写了
简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272
谈谈量子化学中基组的选择
http://sobereva.com/336
Q:
那也要用这个基组重新优化结构，还是直接可以用这个B972/pcSseg-1基组计算核磁就可以了呢
A:
仔细仔细仔细看完博文你就懂了
----------------------------------------------------
2016.06.17 06:12:18
Q:
http://pubs.acs.org/doi/pdf/10.1021/jp1117135麻烦哪位帮我下载一下，谢谢
A:
看群共享里群规，自行用scihub解决
Q:
scihub？
A:
求sci-hub能下文献一次，违规已记录，扣6分，扣满12分将被移出本群，详见群共享中的群规（含sci-hub用法）
----------------------------------------------------
2016.06.17 06:19:10
Q:
请问gauss可以把化学键的长度标出来吗？
A:
不能直接标图上
----------------------------------------------------
2016.06.17 06:19:52
Q:
老师，我计算的几个过渡态结构，加关键词 int=ultrafine 与不加这个关键词所优化出来的过渡态结构之间的能量差差别很大，我应该怎么处理？
A:
当然用int=ultrafine的
Q:
那就是说，我得把这篇工作里面所有需要计算的结构都在 int=ultrafine 下计算？
但是有几个过渡态结构在我加了这个关键词之后反而出现了多个虚频，怎么破？
A:
那倒无所谓
Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法
http://sobereva.com/278
优化用这个算频率也必须用这个
Q:
算频率也需要加上这个关键词？
我的结构优化和频率计算是分开的
A:
显然啊
加不加这个势能面都不一样
[图片]
这种问题培训班全都会讲
Q:
多谢老师提醒
可是对于体系很大的结构，如果结构优化和频率计算一起的时候，计算频率会死掉，所以我都分开计算，为了避免死掉，会在计算频率的时候加上关键词 int=acc2e=11，这样也是不对的吧？
A:
别写这个，没用
分开算的时候freq用的关键词和优化一致就完了
Q:
如何判断一个结构是否需要加上int =ultrafine，我不大可能把所有的结构都在这两种情况下算一遍呀
A:
为保险就都加上
----------------------------------------------------
2016.06.17 06:41:37
Q:
各位老师，请问方法和基组可以随便搭配吗（在查文献后的得到的方法和基组之间）？谢谢，或者说有些方法不能搭配某些基组，有这样的说法吗
A:
谈谈量子化学中基组的选择
http://sobereva.com/336
----------------------------------------------------
2016.06.17 06:55:44
Q:
请问Sob老师一个问题，基态和激发态优化的结构显示激发态的结构二面角更小，但审稿人提出一般来说激发态的构型更扭曲，这应该怎么解释呢？
A:
反例很多啊，比如联苯
[图片]
Q:
非常感谢，请问是否有相关文献？
A:
Chem. Soc. Rev., 2013, 42, 845--856
----------------------------------------------------
2016.06.17 07:01:30
Q:
[图片][图片]老师我跑EM时候就显示inconsistent 是不是我的mdp里面算法有问题
top里面也有问题
A:
初始结构不合理、拓扑文件不合理、模拟参数不合理都有可能
和最小化算法本身没关系
----------------------------------------------------
2016.06.17 07:11:57
Q:
请问下，我在linux自己用户下，设置环境变量LD_LIBRARY_PATH后，所有命令都用不了了，显示Segmentation fault。。。请问下我应该怎么解决。因为现在vi都用不了，不能直接修改.bashrc文件
A:
如果你是root，用安装盘进恢复模式后操作
如果不是root，让root给你改
Q:
好的，谢谢
如果，我不是root。。。是不是必要要让root权限的人改啊。。。不然就解决不了
A:
那你试试把.bashrc拷出来，改完了放回去覆盖
比如通过ssh传到windows下
Q:
我试了。。。外部的Filezilla软件都链接不上它了。。。然后我进入root下的其它用户，想把文件copy进去也不行,但可以copy出来。。用scp
用绝对路径运行/usr/bin/vim 试试？
也不行，刚刚试的
或者echo "" >/.bashrc 覆盖掉bashrc
再重写
~/.bashrc
我也有两次这情况了…
A:
echo都已经把你的.bashrc重写了。你登出再登进来，看看原先命令能用了没
----------------------------------------------------
2016.06.17 07:14:50
Q:
看来我们老板还是很好的，她问我谁讲，我说卢天大神，她说好的去吧，这个我听过。PS.我老板搞光传感的
A:
/微笑
Q:
https://www.icms2016.org/?lang=zh，这个会议应该是正规的吧，有去的么？[表情]
A:
这是正规的，一看网站制作质量就知道
----------------------------------------------------
2016.06.17 17:05:27
Q:
ccsd（t）太吃硬盘了。有什么办法么？
A:
要么尝试DLPNO-CCSD(T)，要么用双杂化代替，要么减小点基组
----------------------------------------------------
2016.06.17 17:15:45
Q:
过渡态几何优化一直不能收敛，加上int=ultrafine scf=verytight后收敛了。算能垒的话是不是也要把反应物加上int=ultrafine scf=verytight重新优化下？
我把收敛后的结构取出，不加int=ultrafine scf=verytight, 结果又收敛不了
A:
/aini int=ultrafine就够了
反应物没必要非得也用ultrafine优化
----------------------------------------------------
2016.06.17 18:14:39
Q:
同志们，请教一个问题，我已经把一个结构的ELF算出来了，保存了它的.cub文件，用什么软件作图比较好看一些？谢谢
A:
VMD或chimera
例
在Chimera中用不同的颜色显示不同的等值面
http://sobereva.com/168
----------------------------------------------------
2016.06.17 18:42:25
Q:
想问一下比较吸附能力强弱的时候无论在几何数据还是AIM电子密度上都是A>B 可是为什么偏偏做BSSE校正算吸附能的时候出现B>A呢
A:
不考虑BSSE，直接算相互作用能谁大？
Q:
还是A>B
A:
也许没有正确使用counterpoise，或数据读取不对
Q:
[图片]
这是我的输入文件跟数据我计算直接有最后优化后的二聚体减去未优化的片段1 2
哦错了那个输入文件多重度是-1 ，1 -1，1 0，1
上面输入错了
但是结论还是错的
没啊正负直接1.577啊
我也觉得纳闷怎么可能吸附那么大
难道是电荷不应该设为-1 而是0 么
A:
甭用counterpoise就完了，省得不熟悉细节还把数据报错了
这么点体系，索性直接用jul-cc-pVTZ，弥散函数已足够多，此时根本不需要考虑BSSE问题。何况体系又带电荷，静电相互作用占主导，更不需要考虑BSSE。再加上，才B3LYP这么烂的泛函，算出的精度水平根本不配考虑BSSE问题
看
乱谈DFT-D
http://sobereva.com/83
谈谈量子化学中基组的选择
http://sobereva.com/336
----------------------------------------------------
2016.06.17 19:17:26
Q:
[图片]
[图片]是什么方法？
[图片]
[图片]是不是相当于把这个里面的QCISD换成CCSD?
,
[图片][图片]是不是相当于把这个里面的QCISD换成CCSD?
A:
这里写了
谈谈量子化学中基组的选择
http://sobereva.com/336
谈谈能量的基组外推
http://sobereva.com/172
UCCSD(T)就是基于UHF做CCSD(T)
Q:
，您好！
请问 ,
[图片][图片]是不是相当于把这个里面的QCISD(T)换成CCSD（T）?
A:
显然是
Q:
CCSD（T）和QCISD（T）哪个做单点好一些？
A:
前者
[图片]
Q:
相同的基组用时差不多吧？
好的，知道了，谢谢!
A:
曾经用QCI的人早已纷纷倒戈向CC
----------------------------------------------------
2016.06.17 19:50:02
Q:
[图片]
请问，gaussview这个截面，怎么保存数据啊
A:
如果你要获得平面上每个点的数据，没辙，必须用Multiwfn
----------------------------------------------------
2016.06.17 20:23:23
Q:
如果理论计算跟实验相比，吸收光谱的两个峰，一个红移30nm，一个蓝移30nm，是不是没有什么太好的解决办法？
A:
恩
Q:
30nm的差距是不是很大了？
A:
还行
非得峰位置算那么准其实没必要
光谱曲线形状能对上更重要
Q:
问题是我的这个体系曲线也不太好，我怀疑是优化结构的问题
因为与之相关的一个体系，非常好（给你看过）
[图片]
A:
也很有可能
----------------------------------------------------
2016.06.17 20:58:15
Q:
sob老师，化学反应需要的自由能是通过这个计算的吧 Thermal correction to Gibbs Free Energy=
A:
看此文第三节
谈谈隐式溶剂模型下溶解自由能和体系自由能的计算
http://sobereva.com/327
----------------------------------------------------
2016.06.17 21:05:16
Q:
[图片]求助一个输入问题：我用6-31g*优化好的结构然后用6-311++g**计算一个单点能量用来最为计算结合能的数值，是不是输入文件这样写就可以了呢？
还是直接拿6-31g*优化还的结构直接再用# b3lyp/6-311++G**这样的形式来计算各单点能量就完事了呢？
A:
不要写opt，用//直接就自动用opt了
----------------------------------------------------
2016.06.17 21:06:05
Q:
sob老师，一个分子含有三个金属，其中两个都是 4 个单电子，另外一个没有单电子，用Gaussian算单点能，多重度指认 9。算完看单电子占据的 8 个轨道，其中两个轨道是在那个没有单电子的金属附近，这个结果是不是就有问题了？能不能通过片段的方法分别指认多重度解决？或者有没有其他方法？
A:
算出来是什么样就是什么样
最好直接看自旋布居数或者自旋密度而不是看轨道
Q:
out里面自旋密度是不是[图片]这个？
[图片]
[图片]
[图片]
A:
这叫自旋布居
量子化学中的一些常见不良写法和用词
http://sobereva.com/298
虽然高斯可以用片段初猜，但是该是什么样收敛完了还是什么样
----------------------------------------------------
2016.06.17 21:14:16
Q:
//的意思就已经表示用后面的6-31g*优化用6-311++g**计算单点能量了吗？老师
A:
y
----------------------------------------------------
2016.06.17 21:56:43
Q:
高斯com文件：
[图片]
log文件最后：
体系49个原子。请问这个出错信息说磁盘空间不够，请问该怎么解决呢？
[图片]
组里集群上算的
除了这个原因，还有没有可能是精度高，体系大导致的，能用maxdisk解决吗？
A:
当然可能，算不了就甭愣算
要么用orca
----------------------------------------------------
2016.06.17 22:26:35
Q:
记得群主有篇博文讲硬件配置
A:
单路计算化学攒机配置推荐
http://sobereva.com/243
----------------------------------------------------
2016.06.17 22:27:16
Q:
http://www.sciencedirect.com/science/bookseries/00653276/69
麻烦大家帮我下载一下这篇文献吧
谢谢！
A:
[图片]
自行下载
----------------------------------------------------
2016.06.17 22:30:31
Q:
请问荧光寿命能计算得到吗？
A:
Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法
http://sobereva.com/314
这里写了
培训班里也会讲
Q:
中国量化培训，我只相信sob [表情]
三重态和单重态之间是自旋禁阻的，仅当考虑旋轨耦合时振子强度才不为零，但数值依然很小，所以磷光寿命远比荧光寿命长。
仅提到这一句荧光寿命
骗子太多，傻子不够用
具体能不能计算？如何计算？
培训班回去参加，但是现在您先给个简要的回答呗
A:
文中给了寿命的公式啊
Q:
呵呵，这样说会让sobereva压力很大的
A:
[图片]
----------------------------------------------------
2016.06.17 22:51:04
Q:
--------Markdown Markdown--------
A:
--------Markdown Markdown--------
Q:
我给大家发一个IRC曲线让大家开开眼
[图片]
CPU时间白费了……
将近29天……
40
stepsize=40
A:
这是粗暴跑IRC，不震荡才怪
步长越小，或精确重算一次Hessian所隔步数越小，IRC曲线月平滑
Q:
默认的stepsize感觉太磨蹭，点过于密了，于是……
嗯，谢sob老师，谢谢楼上各位
A:
关键看走IRC是什么目的，对于确认TS无误的目的，粗暴跑出来，确认是往极小点走就够了
Q:
嗯，老师，我就想确认一下过渡态，是不是连接产物反应物，这么震荡，也无所谓吗
A:
无所谓
Q:
好的。谢谢老师
老兄，你那算法是研究什么的呢
l9999要加noeigentest。几乎总是要加。
判断体系自旋多重度要根据相应结构化学知识
A:
无非又是震荡不收敛
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164
Q:
我在优化过程中出现这个错误，请问可能是什么？
[图片]
A:
Sobereva(190258442) 15:11:26
无非又是震荡不收敛
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164
Q:
那是不是我把步长设长点，可能会收敛？
A:
看完博文再问
----------------------------------------------------
2016.06.17 23:06:24
Q:
遇到了一个问题NO2+CO=NO+CO2反应，找过渡态，opt+freq,总是l9999，无法得到频率信息。试过各种方法都不行。这里NO2有单电子，自旋多重度=2,这么做对吗？谢谢！
A:
无非又是震荡不收敛
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164
Q:
好，我再试试，谢谢！
A:
[图片]
很简单就找到了
Q:
有什么好的方法可以判断体系自旋多重度？只能一个个试吗？
A:
自旋多重度就是2
----------------------------------------------------
2016.06.17 23:19:26
Q:
sob老师，G=H - TS 这个H就是H=E + RT 这个E应该是单点加ZPE矫正的能量吧？
A:
H=U+RT
U是单点能加上内能热校正量
ZPE只是内能热校正量在0K时的值
----------------------------------------------------
2016.06.17 23:20:40
Q:
[图片]
[图片]
请问一下，这种轨道图就是把结构拆开算单点得到的吗？》
[图片]sob, 他每个片段都给出来一个电子的数目,括号里面的值，是怎么搞出来滴？
A:
@ D调图中根本就没说清楚，不要去纠结
A:
@ D调就是对两个反应物分别计算分子轨道，然后和到一张图上
Q:
应该不需要优化吧，就把整个分子优化好了以后，拆开算一下单点就好了吧？
A:
结构看不清楚。
应该就是把反应物分子分别优化分别计算单点
图片描述不清
貌似就是个复合物，算的片段电荷得到电荷转移量
----------------------------------------------------
2016.06.17 23:24:05
Q:
Please submit your 3D molecular models in one of the following formats: PDB, PSE or MOL/MOL2.请问图片格式可以转化吗
A:
把分子结构存成这些格式之一提交
一般常用pdb
Q:
sob老师，我的图片格式，可以转为PDB格式吗
A:
不能
----------------------------------------------------

复制代码