计算化学公社

 找回密码 Forget password
 注册 Register
Views: 9470|回复 Reply: 0
打印 Print 上一主题 Last thread 下一主题 Next thread

[群聊Q&A整理] 公社群聊Q&A整理:2016.07.12 19 ~ 2016.07.20 20 Concate

[复制链接 Copy URL]

3097

帖子

29

威望

1万

eV
积分
17098

Level 6 (一方通行)

跳转到指定楼层 Go to specific reply
楼主
相关概况可阅读:http://bbs.keinsci.com/thread-2022-1-1.html


暂无全文附件

本次附件(附件内容与下文一致):
2016.07.12 19 ~ 2016.07.20 20 Concate.txt (53.51 KB, 下载次数 Times of downloads: 4)

  1. 2016.07.12 18:10:33
  2. Q:
  3.     HF exchange percentage (HF%) 这个是怎么算的

  4. A:
  5.     不同DFT泛函的HF成份一览
  6.     http://sobereva.com/282

  7. ----------------------------------------------------
  8. 2016.07.12 18:25:00
  9. Q:
  10.     各位老师们,我想请问:
  11.     在一个多步反应路径中,其中有一步的吉布斯能变大于零,整体的吉布斯能变远远小于零,这条路径还能发生么?
  12.     还有一个路径有两步且每一步的吉布斯能变都小于零,整体也小于零但比三步的那个能变大,可以得出三步路径比两步路径更容易发生么?
  13.     各位老师们,我想请问:
  14.     在一个多步反应路径中,其中有一步的吉布斯能变大于零,整体的吉布斯能变远远小于零,这条路径还能发生么?
  15.     还有一个路径有两步且每一步的吉布斯能变都小于零,整体也小于零但比三步的那个能变大,可以得出三步路径比两步路径更容易发生么?

  16. A:
  17.     1 能发生
  18.     2 意义不明

  19. Q:
  20.       “2 意义不明” 意思是不能去直接比较这两个路径么?

  21. A:
  22.     完全没道理三步的更容易发生

  23. Q:
  24.     哦
  25.       老师,“1 能发生”  文献报道都是说不发生,您支持发生是从哪方面考虑的呢?直接因为总吉布斯能小于零么?

  26. A:
  27.     凭什么不能发生,根本没有依据
  28.     除非势垒特别高

  29. ----------------------------------------------------
  30. 2016.07.12 19:03:01
  31. Q:
  32.     哪位老师或同学可以给我一些关于自旋密度的定义,分析之类的资料吗??
  33.     哪位老师或同学可以给我一些关于自旋密度的定义,分析之类的资料吗??

  34. A:
  35.      这里面有定义和例子
  36.     幻灯片:Multiwfn与波函数分析简介
  37.     http://sobereva.com/239

  38. ----------------------------------------------------
  39. 2016.07.12 21:44:34
  40. Q:
  41.     Sob老师,您好!我在使用您编译的dushin-linux计算分子重组能的时候出现了如下报错,[图片],请问这可能是什么原因造成的?@Sobereva

  42. A:
  43.     不好说

  44. Q:
  45.     会跟基组有关系吗?@Sobereva

  46. A:
  47.     不像

  48. ----------------------------------------------------
  49. 2016.07.12 22:04:07
  50. Q:
  51.     请教各位前辈,混合基组在gamess中怎么写?有没有相关资料可以推荐一下?

  52. A:
  53.      http://bbs.keinsci.com/forum.php?mod=viewthread&tid=745

  54. ----------------------------------------------------
  55. 2016.07.12 22:13:07
  56. Q:
  57.     [图片]大家好,请问在高斯软件中,怎样设置键变成黑色呢?[表情][表情]

  58. A:
  59.      file-preference[图片]

  60. ----------------------------------------------------
  61. 2016.07.12 22:23:04
  62. Q:
  63.     [图片]
  64.     循环很久了,不出结果,怎么办?
  65.     [图片]
  66.     [图片]
  67.     那为什么循环不出来

  68. A:
  69.     明显是SCF收敛困难
  70.     解决SCF不收敛问题的方法
  71.     http://sobereva.com/61

  72. Q:
  73.     泛函不适合?

  74. A:
  75.     此文一定要看
  76.     常见的多余的和被滥用的Gaussian关键词
  77.     http://sobereva.com/331

  78. ----------------------------------------------------
  79. 2016.07.12 23:10:46
  80. Q:
  81.     [图片]计算E(adia)就是激发态几何优化的HF减去基态的几何优化HF
  82.     [图片]计算E(adia)就是激发态几何优化的HF减去基态的几何优化HF?请教一下

  83. A:
  84.      不对。激发态几何优化给出的HF=后面的能量是参考态(基态能量)。应当用E(TD-HF/TD-KS) 后面的,这才是激发态能量

  85. ----------------------------------------------------
  86. 2016.07.12 23:24:30
  87. Q:
  88.     老师:您好!请问这种DOS图,[图片],感觉是分解成了轨道贡献,请问这种图怎么做呢?谢谢!我参照了您手册中的DOS做法,没找到 ,请您指点指点  谢谢···[表情]

  89. A:
  90.      Multiwfn做DOS时候自定义片段,让片段对应要考察的轨道,仔细看手册4.10节例子

  91. ----------------------------------------------------
  92. 2016.07.12 23:43:21
  93. Q:
  94.     [图片]
  95.     黄色的是硫原子,数值是MULLIKEN电荷,当H远离S原子时,它应该是以质子形式的远离的。但把H逐渐远离S原子,算MULLIKEN电荷时,H总是是原子形式,不怎么带电荷。距离越远,越是中性。
  96.     这是怎么回事?
  97.     [图片]
  98.     黄色的是硫原子,数值是MULLIKEN电荷,当H远离S原子时,它应该是以质子形式远离的。但把H逐渐远离S原子,算MULLIKEN电荷时,H总是是原子形式,不怎么带电荷。距离越远,越是中性。
  99.     这是怎么回事?
  100.     [图片]
  101.     黄色的是硫原子,灰色的是碳原子,白色的是H,,数值是MULLIKEN电荷,当H从S上解离时,应该是质子形式。取S-H不同的距离,计算单点,却发现:H逐渐远离S原子时,H总是是原子形式,不怎么带电荷。距离越远,越是中性。
  102.     这是怎么回事?
  103.     怎么计算模拟H以质子形式从S上解离这个过程?

  104. A:
  105.      首先确保没用弥散函数
  106.     如果没用弥散函数,当前解离相当于S-H均裂,可以考察拉得更远的情况,如果实际上就是C2H5S-和H+这种状态能量最低,也会自然而然收敛到那种态
  107.     基态是A- B+还是A·  ·B距离有关,在某个距离左右可能发生这两种透热态的交叉
  108.     基态是A- B+还是A·  ·B跟AB距离有关,在某个距离左右可能发生这两种透热态的交叉

  109. Q:
  110.     改变体系的自旋多重态可以模拟这种异裂吗?
  111.     以氢负离子形式离去,以氢质子形式离去。分别如何模拟(量化水平)?

  112. A:
  113.     和自旋多重度没关系。
  114.     有的距离范围基态可能是均裂,有的距离范围基态可能是异裂。
  115.     假设某个距离下,基态不是异裂而是均裂,但你又非得算异裂,就只能用CDFT
  116.     [图片]

  117. Q:
  118.     高斯还不能算CDFT?

  119. A:
  120.     只有nwchem和q-chem支持
  121.     谈谈约束性DFT (CDFT)
  122.     http://sobereva.com/271

  123. ----------------------------------------------------
  124. 2016.07.13 00:14:38
  125. Q:
  126.     各位下午好,想请教各位大神个问题:Jarzynski’s equality 这个名词如何理解呢?

  127. A:
  128.      Jarzynski equality:从平衡态A到某态B的自由能变ΔG=G_B-G_A<=W,其中W是外力做功,只有过程是平衡过程等号才成立。Jarzynski equality将不等式化为等式,exp(-ΔG/(kT))=<exp(-W/(kT))>,T代表在A时的温度,<>代表是以各种方式干这个活儿的平均值。如果过程快,非平衡特征强,就得需要做很多次取平均才精确。如果过程慢,接近平衡态,只需要很少次数即可。比如通过SMD将受体-配体拉开,如果拉得十分慢,得到的"做功-拉伸距离"曲线就是"结合自由能-拉伸距离"曲线,如果拉得快,就要拉好多次把结果取平均。

  129. ----------------------------------------------------
  130. 2016.07.13 00:18:33
  131. Q:
  132.     请问下在做弱相互作用时,由于之前用的B3LYP/6-31G**,现想加入D3矫正,只算个单点可以吗?还是需要做下重新优化。

  133. A:
  134.      色散作用主导时,优化时候就得考虑

  135. ----------------------------------------------------
  136. 2016.07.13 00:29:31
  137. Q:
  138.     请问gauss view画出来的静电势图如何导出成tif格式呢?

  139. A:
  140.      保存图片时候直接就能选tif啊
  141.      [图片]

  142. ----------------------------------------------------
  143. 2016.07.13 01:02:21
  144. Q:
  145.     请教老师,分析AdNDP时,这个错误要怎么解决[图片]谢谢

  146. A:
  147.     我只回答Multiwfn做adndp的问题。其它程序谁开发的问谁

  148. ----------------------------------------------------
  149. 2016.07.13 01:20:39
  150. Q:
  151.     请问Marcus 理论算的是分子内还是分子间的电子转移

  152. A:
  153.     一般用来讨论分子间的

  154. ----------------------------------------------------
  155. 2016.07.13 03:24:52
  156. Q:
  157.     我记得S老师写过一个关于高斯操作的贴子,比如,一行里面写多个g09 命令。这个贴子名字是什么呀?

  158. A:
  159.     使用Gaussian时的几个实用脚本和命令
  160.     http://sobereva.com/258

  161. ----------------------------------------------------
  162. 2016.07.13 03:38:36
  163. Q:
  164.     想计算ADCH电荷,wfn和fch文件是不是都可以用作输入文件?还是必须是fch文件?ADCH电荷是否可用于求算激发态分子的原子电荷? 谢谢老师!

  165. A:
  166.     都行
  167.     可
  168.     激发态时候最好用.wfn文件,否则用fch还得额外多做一步

  169. Q:
  170.     老师,计算电荷的方法用HF/6-31G*还是B3LYP/6-31G*?是不是用HF更好些? 麻烦老师了,总是打扰您。

  171. A:
  172.     若无特殊情况用B3LYP

  173. ----------------------------------------------------
  174. 2016.07.13 04:14:56
  175. Q:
  176.     我突然发现内存对高斯计算速度也有很大作用的影响
  177.     192G的比64G的就是算得快
  178.     内存设置均约为总内存 90%

  179. A:
  180.     肯定的

  181. Q:
  182.     嗯,就是说内存如果不是完全充足,高斯会降低速度来满足内存设置

  183. A:
  184.     程序实现同一个目的往往有多种算法,有耗内存高但是速度快的,也有耗内存低但速度慢的,程序会根据分配的内存量来确定最优的选择。

  185. Q:
  186.     [表情]
  187.     MS有这个选项, 我一直以为高斯是占硬盘

  188. A:
  189.     还比如,内存大的情况一个事情可能分三批算完,内存小的情况需要分7批算完,前者耗时明显低

  190. Q:
  191.     没想到也是占内存

  192. A:
  193.     可不


  194. ----------------------------------------------------
  195. 2016.07.13 04:14:56
  196. Q:
  197.     老师,它已经停了,没有继续算

  198. A:
  199.     重算一遍,还不行就换个版本,或者换关键词
  200.        
  201. ----------------------------------------------------
  202. 2016.07.13 06:43:19
  203. Q:
  204.     各位老师:通过NBO看轨道间转移的电荷数从哪看呀
  205.     [图片]我以前都是通过二阶微扰能看趋势

  206. A:
  207.     自行根据公式算
  208.     [图片]

  209. ----------------------------------------------------
  210. 2016.07.13 18:05:38
  211. Q:
  212.     sob,老师,什么样的体系适合Anndp分析呢@Sobereva
  213.     我觉得你应该问 从AdNDP分析 能得到什么信息

  214. A:
  215.     AdNDP常用在有显著的局部多中心离域的体系

  216. ----------------------------------------------------
  217. 2016.07.13 18:19:55
  218. Q:
  219.     请问一下gromacs中的键能是怎样提取的?
  220.     算出来的是coul和lj
  221.     没有键能

  222. A:
  223.     本来力场就不可能算出键能

  224. Q:
  225.     老师,力场里面不是有键参数吗?

  226. A:
  227.     力场成键项的能量和键能是完全两码事

  228. Q:
  229.     我想算rna二面角的扭转能
  230.     也就是不能计算bonded energy?
  231.     只能提取coul 和lj吗?

  232. A:
  233.     记得是
  234.     把键长、二面角变化数据搞出来,自己代进公式一算就有了

  235. ----------------------------------------------------
  236. 2016.07.13 19:01:22
  237. Q:
  238.     老师,我刚才我试了一下,6-31G**可以算铜,但是不能算Mo,我还要算巯基乙酸与Mo的化合物,该怎么选择基组 ,谢谢

  239. A:
  240.     用lanl2TZ,NBO分析还不行就SDD,还不行就def2SVP,还不行就全电子相对论计算(到时候再说)

  241. Q:
  242.     用lanl2TZ提示Subroutine NAOANL could not find a S-type valence orbital on atom Cu   11.
  243.     IVAL :  1  1  1  0     M :  1     LA :  1
  244.     Error@Sobereva
  245.     请问,Sob老师说的是不是,Cu用SDD,其他原子还是用6-31G**啊?

  246. A:
  247.     现在说的是Mo
  248.     Cu不都解决了

  249. Q:
  250.     /抱拳
  251.     您的意思是对Mo适用的基组对Cu也适用,对吧?@Sobereva
  252.     好的,谢啦

  253. A:
  254.     你对Cu用SDD没准还碰上那个报错
  255.     不是用的赝势基组合理不合理的问题,是NBO 3.1本身就有bug,你得想办法绕过那个bug
  256.     全电子基组100%没问题

  257. ----------------------------------------------------
  258. 2016.07.13 19:07:52
  259. Q:
  260.     [图片]CBS-QB3算到这一步就慢的一动不动的,24核也算不动

  261. A:
  262.     算不动就甭愣算
  263.     降低点档次

  264. ----------------------------------------------------
  265. 2016.07.13 19:12:03
  266. Q:
  267.      sob老师,td 计算triplet的时候,用ub3lyp 还是b3lyp更好些?

  268. A:
  269.     基态是闭壳层单重态就不要写U。而且用了U再用triplet关键词毫无效果

  270. ----------------------------------------------------
  271. 2016.07.13 19:14:10
  272. Q:
  273.     请教各位 类似这样的图用什么软件做最方便 我只想做CPK和线性的组合那部分
  274.     [图片]

  275. A:
  276.     VMD

  277. Q:
  278.     [图片]请问两种显示方式叠合之后要突出某一部分 除了删除部分原子 还有别的办法么

  279. A:
  280.     VDW方式显示的representation的selected atom只选要凸显的原子范围就完了,比如serial 1-30
  281.     serial 1 to 30

  282. ----------------------------------------------------
  283. 2016.07.13 19:14:39
  284. Q:
  285.     我要看反应过程的能量以及S和Cu之间电子转移,老师有啥建议吗?不用NBO有其他方法算吗?@Sobereva

  286. A:
  287.     能量和NBO本身就没直接关系
  288.    
  289.     电子转移算原子电荷、做密度差图、做CDA什么的方法多了去了

  290. ----------------------------------------------------
  291. 2016.07.13 19:17:53
  292. Q:
  293.     vmd能看高斯的输出吗

  294. A:
  295.     不能

  296. ----------------------------------------------------
  297. 2016.07.13 19:18:21
  298. Q:
  299.     [图片]CBS-QB3中CCSD(T)/6-31+G(d') 算20几个原子算不动。用老师上课讲的较好级别的基组比如:CCSD(T)/cc-PVTZ,一般能算到多少个原子?

  300. A:
  301.     不用DLPNO之类黑科技顶多也就20个原子到头了

  302. Q:
  303.     DLPNO-CCSD(T)只有用ORCA?

  304. A:
  305.     恩

  306. ----------------------------------------------------
  307. 2016.07.13 19:22:01
  308. Q:
  309.     老师,有一个问题,这个是我用multiwfn看的等值面,想保存这个图片,但是右击图片,图片就没了,怎么可以保存呢?要截图吗?

  310. A:
  311.     菜单上有选项保存图片

  312. ----------------------------------------------------
  313. 2016.07.13 19:39:20
  314. Q:
  315.     [图片]
  316.     请问老师这怎么破

  317. A:
  318.     opt=cartesian

  319. ----------------------------------------------------
  320. 2016.07.13 20:57:06
  321. Q:
  322.     请教各位,请教6-31g++(d)不收敛但是6-31g++(d,p)可以收敛的原因

  323. A:
  324.     数值巧合

  325. ----------------------------------------------------
  326. 2016.07.13 21:29:06
  327. Q:
  328.     群主不知道可有好的工具生成碳120和240这样的工具?

  329. A:
  330.     CaGe可以试试
  331.     生成富勒烯的螺旋算法简介以及使CaGe中的编号与Fowler-Manolopoulos编号相符的方法
  332.     http://sobereva.com/104

  333. ----------------------------------------------------
  334. 2016.07.13 23:41:04
  335. Q:
  336.     --------Markdown Markdown--------

  337. A:
  338.     --------Markdown Markdown--------
  339.      做变分过渡态理论计算需要IRC上一系列点的自由能

  340. ----------------------------------------------------
  341. 2016.07.13 23:41:10
  342. Q:
  343.     获得零点能,然后获得势能面,寻找范德瓦尔斯complex的位置
  344.     #p irc=(maxpoints=50,calcfc,maxcycle=100) bhandhlyp/6-311g(d,p) geom=connectivity freq=noraman我这样算的,但是好像没有输出freq

  345. A:
  346.    
  347.     可以用这个得到.xyz轨迹,然后用molclus批量执行
  348.     将Gaussian的IRC任务输出转换为.xyz轨迹文件的工具
  349.     http://sobereva.com/285
  350.    
  351.     或者借助此文的IRCsplit生成一批.gjf
  352.     产生Gaussian09 IRC和SCAN任务每个点的波函数文件的工具
  353.     http://sobereva.com/199

  354. ----------------------------------------------------
  355. 2016.07.13 23:43:54
  356. Q:
  357.     --------Markdown Markdown--------

  358. A:
  359.     --------Markdown Markdown--------
  360.      y
  361.          初猜给得不好,就不容易收敛

  362. ----------------------------------------------------
  363. 2016.07.14 00:01:18
  364. Q:
  365.      老师您好,multiwfn里计算出的分子电场只是点电荷电场吗?包括偶极场吗?

  366. A:
  367.     不一回事。
  368.     每个点都有静电势值,公式就是前面我幻灯片里贴出来的,电场就是静电势的梯度

  369. Q:
  370.      那对[图片]求梯度,就可以认为是相当于点电荷产生的电场?我加了电场后有了极化,有偶极场,其实这时,multiwfn里再计算出电场也还是没包含偶极场。老师,是这样吗?

  371. A:
  372.     图中公式用的是连续的电子密度,而不是把体系当成个点多极矩来考虑,牵扯不到什么偶极场

  373. ----------------------------------------------------
  374. 2016.07.14 00:04:58
  375. Q:
  376.      老师,您好,如果做NTO计算的话是不是就不能进行electron-hole 分析了?

  377. A:
  378.     如果你用saveNTO把NTO存到chk里就不能用这个chk做electron-hole分析了
  379.     electron-hole分析是基于MO来做的

  380. ----------------------------------------------------
  381. 2016.07.14 00:52:33
  382. Q:
  383.      sob 老师,我用ub3lyp  opt    自旋0  3优化得到的t1 和用 td(triplet,root=1)  自旋 0  1得到的能量差别不小啊。

  384. A:
  385.     这种设定并没什么意义,收敛完了都一样

  386. Q:
  387.      老师,td 优化s1或者t1 最后HF=这个能量就是他本身能量吧。真奇怪,和直接计算得到的有所差别

  388. A:
  389.      那是那个结构下基态DFT能量

  390. Q:
  391.     这个基态dft能量,t1下应该就是t1的能量吧。
  392.     linux 格式,还有空格

  393. A:
  394.      TDDFT算T1时基态能量当然不是T1能量

  395. Q:
  396.     明白了,谢谢Total Energy, E(TD-HF/TD-KS)  这个才是
  397.     这个基态能量只是输入结构的能量
  398.     http://bbs.keinsci.com/forum.php?mod=viewthread&tid=3163&extra=page%3D1 计算t1-s1 必须要这样改么?,把IF (DE .lt. TDE) THEN改成IF (.true.) THEN

  399. A:
  400.     改它干嘛
  401.     本来就能交叉
  402.     不想要“能量差必须很小”这个条件才那么改

  403. Q:
  404.     请问老师,用DFT方法计算激发态和TDDFT方法计算激发态关键词有什么不同呢?

  405. A:
  406.     Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法
  407.     http://sobereva.com/314
  408.        
  409. ----------------------------------------------------
  410. 2016.07.14 01:57:59
  411. Q:
  412.     请问我计算一个约50个原子以上的相互作用体系,其中4个氧原子,其他均为碳原子和氢原子。使用4核计算机单机计算,使用高斯计算。需计算20个以上构型,需1到1个半月算完。使用M06-2X方法,使用6-31+G**进行几何构型优化,使用6-311+G**计算单点能。请问各位老师同学这样算可以吗?

  413. A:
  414.     具体什么级别最合适还得看体系具体结构,不是光原子数就能决定的
  415.     优化可以不用弥散。即便色散作用明显,用6-311G**也比6-31+G**效率更高结果更好

  416. ----------------------------------------------------
  417. 2016.07.14 04:20:19
  418. Q:
  419.     sob老师,multiwfn主功能18里的电子空穴积分对于HOMO--LUMO主导(95.9%)的跃迁是不是基本等于min{HOMO^2,LUMO^2}的积分呢?我用主功能13-11-22,算出来是0.0477,而电子空穴积分是0.11,差别有点大,该如何解决呢?

  420. A:
  421.     如果13-11-22你是对hole和electron的格点数据做的,只可能是格点设定不一样所致。
  422.     没必要自行用13-11-22去搞

  423. Q:
  424.     恩,我只是想确认下我的理解有没有错。13-11-22是对HOMO^2和LUMO^2做的。全都是用的medium格点精度啊,还有什么参数会影响吗?另外,刨去手动操作的时间,感觉18-1的计算速度相比于13-11-22要慢些呢

  425. A:
  426.     18-1要考虑各种轨道产生格点数据,所以显然慢

  427. Q:
  428.     感觉既然HOMO--LUMO主导(95.9%),那么两种方法算出来结果应该差别不大啊

  429. A:
  430.     差那么百分之几也可能会有影响
  431.     你还不如直接用18-1产生的electron和hole的cube文件,再用13-11-22去计算min的积分,对比18-1直接输出的

  432. Q:
  433.     这个是一样的,我试过了
  434.     差百分之几也会有这么大的影响?老师能不能详细解释一下原因呢?非常感谢/西瓜
  435.     差百分之几也会有这么大的影响?老师能不能详细解释一下原因呢?非常感谢[表情]

  436. A:
  437.     那就是说明考虑所有轨道跃迁的electron-hole和只考虑HOMO->LUMO的差异造成那些结果差异
  438.     你可以先从图形上对比electron和LUMO^2、hole和HOMO^2,应该图形上看不出太大差异,但是定量到min的积分数值上会有那些影响

  439. ----------------------------------------------------
  440. 2016.07.14 05:39:26
  441. Q:
  442.     群主可能也不喜欢这个[表情]表情

  443. A:
  444.     是的
  445.     很莫名奇妙

  446. ----------------------------------------------------
  447. 2016.07.14 05:47:19
  448. Q:
  449.      老师好,研究有质子转移的反应,是不是没必要去刻意指定H为质子?整个体系在同一个电荷数和自旋多重度下计算就行了?
  450.     计算结果转移过程中H是质子,那么就是质子转移反应,如果不是,那就不是质子转移。

  451. A:
  452.     对

  453. Q:
  454.     质子转移是一个古老而又富有活力的课题,质子转移是最从木和最重要也是最有发展前景的生物化学过程之一
  455.     [图片]
  456.     在转移过程中,H的电荷变化不大,始终在0.47左右。这算是质子转移吗?

  457. A:
  458.      勉强算
  459.     实际体系中不会有完全质子形态存在的情况,电子终究会有一部分离域上去

  460. ----------------------------------------------------
  461. 2016.07.14 05:52:40
  462. Q:
  463.     话说ORCA怎么这么大

  464. A:
  465.     因为用了很多静态库

  466. ----------------------------------------------------
  467. 2016.07.14 05:52:41
  468. Q:
  469.     [图片]
  470.     请问老师,为什么做出来的图一半完整,一半有窟窿呢?

  471. A:
  472.     用Multiwfn+VMD做RDG分析时的一些要点和常见问题
  473.     http://sobereva.com/291
  474.    
  475.     看2.12节

  476. Q:
  477.     老师,我把RDG_maxrho设为0,中间小圆片的窟窿没有了,但两边的原子又出现了我不想要的形状,该怎样修改呀?
  478.     这是修改前[图片]
  479.     修改后[图片]

  480. A:
  481.     RDG_maxrho设成合适的值,或者先设成0,算完了之后看散点图,再决定屏蔽掉指定的sign(lambda2)rho的范围。手册4.100.1节Interactively set RDG value where sign(l2)r is in specific range提到了这点

  482. ----------------------------------------------------
  483. 2016.07.14 06:10:53
  484. Q:
  485.      sob老师,我将一个未优化的构象转化成wfn文件,然后用Multiwfn算CHELPG,算完发现分子坐标变了,但以前我把优化后的构象这样操作,坐标是不会变的,请问老师,中间那一步导致构象会发生变化,谢谢
  486.     有没有方法让坐标保持不变吗?

  487. A:
  488.     是高斯自动把它调整到了标准朝向,看
  489.     谈谈Gaussian中的对称性与nosymm关键词的使用
  490.     http://sobereva.com/297

  491. ----------------------------------------------------
  492. 2016.07.14 16:56:09
  493. Q:
  494.     请问各位高人,GsGrid程序和Sigmaploth程序在哪里可以下载到?

  495. A:
  496.      gsgrid功能早就融入到Multiwfn主功能13了,不要再考虑gsgrid

  497. ----------------------------------------------------
  498. 2016.07.14 17:47:39
  499. Q:
  500.     请问sob老师,PM6-D3H4算弱相互作用时怎么样在输入文件写@Sobereva

  501. A:
  502.     MOPAC里直接写PM6-D3H4就完了

  503. Q:
  504.     老师,我用Gaussian算的。

  505. A:
  506.     高斯不支持

  507. ----------------------------------------------------
  508. 2016.07.14 17:51:59
  509. Q:
  510.      老师您好,用gaussian计算,在m062x/6-31+g(d)下进行结构优化,然后在m062x/6-311+g(d)进行BSSE校正,这样处理可以吗?这个精度发文章会有问题吗?谢谢老师

  511. A:
  512.     没问题
  513.     优化最好改用6-311g(d)
  514.     若牵扯氢键要对氢加极化

  515. ----------------------------------------------------
  516. 2016.07.14 18:00:29
  517. Q:
  518.     我用multiwfn保存ELF的cube文件,然后用VMD观看的时候,一直在出现invalid command name "  "啊,求问老师怎么回事

  519. A:
  520.     怎么操作的?
  521.     直接拖进VMD main窗口就完了,不要拖进控制台

  522. Q:
  523.     我是把文件放到vmd那,然后load visersion state 的

  524. A:
  525.     自行在图形界面设定显示方式
  526.     你那个visualization state的文件未必适合当前情况

  527. Q:
  528.     sob老师,拖进VMD之后,并没有出现ELF等值面图啊,您手册上有关于这方面的操作吗?

  529. A:
  530.     自己在display-representation里面设定
  531.     [图片]

  532. Q:
  533.      你说的是Graphics-representations里面吧?

  534. A:
  535.     y
  536.     不会操作参考这里面VMD的操作步骤
  537.     使用Multiwfn观看分子轨道
  538.     http://sobereva.com/269
  539.     要么chimera显示,专有个博文
  540.     在Chimera中用不同的颜色显示不同的等值面
  541.     http://sobereva.com/168

  542. ----------------------------------------------------
  543. 2016.07.14 18:02:16
  544. Q:
  545.      老师,为了得到偶极矩,我用SVPD基组对只含C、H、O共54个原子的化合物做优化时,由于电脑配置低,参考“解决SCF不收敛问题的方法”,加入 scf=(novarace,conver=6)或scf=(nodiis,fermi),运行约36h后都出现了收敛失败。那么想实现收敛用scf=(xqc,conver=6)会好点吗?也请教一下有没有技巧快速得到可靠的偶极矩值?[表情]谢谢

  546. A:
  547.     就是在6-31G**基础上把极化函数的指数给改了
  548.     你优化时候不用带弥散,算偶极矩再带弥散就行了
  549.     我一直强烈反对在优化时候加弥散,除非有特殊情况

  550. Q:
  551.      老师,我一般是从优化的输出文件中看偶极矩值,是不是优化两次,先不带弥散优化,再带弥散优化一次?请指正

  552. A:
  553.     不带弥散优化,带弥散算单点

  554. Q:
  555.      老师,细问一下,为了得到偶极矩值,那么用6-31g(d)优化,再用def2-SVPD算单点能,由单点能的log文件读取偶极矩值,可以吗

  556. A:
  557.     def2-SVP优化,def2-SVPD算单点
  558.     就是那么读

  559. ----------------------------------------------------
  560. 2016.07.14 18:08:50
  561. Q:
  562.     服务器是不是没有必要配固态硬盘

  563. A:
  564.     看算什么。硬盘读写厉害的任务用固态还是速度有提升

  565. ----------------------------------------------------
  566. 2016.07.14 18:11:25
  567. Q:
  568.     5D 0,5D-2中x^i*y^j*z^k*exp(-ar^2)的i,j,k分别是多少?

  569. A:
  570.     那是球谐函数的情况,此时没法写成这种笛卡尔GTF的形式,必须得做线性组合
  571.     谈谈5d、6d型d壳层基函数与它们在Gaussian中的标识
  572.     http://sobereva.com/51
  573.     球谐型与笛卡尔型Gauss函数的转换关系
  574.     http://sobereva.com/97

  575. ----------------------------------------------------
  576. 2016.07.14 18:25:36
  577. Q:
  578.     各位老师,我想问下氢键和π-π作用一般相互作用能大小有多大的差值,还是不能根据能量差分析这个,要做rdg图看哪

  579. A:
  580.     不是作用类型决定的,而是具体体系决定的,后者也可以比前者大得多

  581. ----------------------------------------------------
  582. 2016.07.14 18:31:28
  583. Q:
  584.     请问老师,用过渡态的关键词得到的结构没有虚频是什么原因
  585.     请问老师,用过渡态的关键词得到的结构没有虚频是什么原因

  586. A:
  587.     过渡态结构不对

  588. ----------------------------------------------------
  589. 2016.07.14 18:32:28
  590. Q:
  591.      老师,我算如下涉及自由基的反应,请问对过渡态结构做自旋密度布居图有意义吗,如果有意义说明什么呢,说明这时单电子已同时分布在N和O原子上吗?
  592.     [图片]

  593. A:
  594.     这是自旋密度图,不是自旋密度布居图
  595.     物理意义就是你说的

  596. ----------------------------------------------------
  597. 2016.07.14 18:38:24
  598. Q:
  599.     老师 文献中的这个最高占据自然轨道的占据数[图片]怎么查找, 谢谢

  600. A:
  601.     pop=NO输出自然轨道和占据数

  602. Q:
  603.     [图片]老师这样写输入文件可以吗

  604. A:
  605.     可

  606. Q:
  607.     我是这样写的,计算完了之后输入,我输入NO 结果查找不到,我应该输入什么关键词可以查找

  608. A:
  609.     那可能正好搜到,仔细看输出文件
  610.     太想当然了

  611. ----------------------------------------------------
  612. 2016.07.14 19:24:32
  613. Q:
  614.     键级两种都算出来,mayer 和Laplacian 放到同一个表里有意义么 我主要想说明从IM TS P过程中关键的几个键的变化

  615. A:
  616.      可以,有

  617. ----------------------------------------------------
  618. 2016.07.14 19:39:34
  619. Q:
  620.     请教老师,怎么感觉MS功能不如高斯强大?
  621.     [表情],专业做计算的 还用什么界面
  622.     错, 专业做计算的才用界面
  623.     比如结构优化, 时时观测

  624. A:
  625.     目的不一样
  626.     纯粹量化范畴的计算高斯更强

  627. ----------------------------------------------------
  628. 2016.07.14 19:46:14
  629. Q:
  630.     sob老师您好,我显示富勒烯碳笼和内部cluster之间的轨道重叠部分的时候,我格点用的是cluster@Cage,然后cage作为第一部分,cluster作为第二部分,我选轨道的时候选的是HOMO重叠,但是显示出来没有重叠啊

  631. A:
  632.     肯定是操作问题
  633.     先把博文里那个例子完整重复一遍
  634.     或者等值面数值调小点再看看,如果怎么调小也显示不出来一定是操作过程不对

  635. ----------------------------------------------------
  636. 2016.07.14 19:48:26
  637. Q:
  638.     [图片] 这东西安装会麻烦?
  639.     还要配置好多环境变量
  640.     不如MS直接把一个破解文件拷贝进去
  641.     就破解了
  642.     medeA也破解了?
  643.     我现在已经被免费软件的环境变量搞疯了
  644.     免费的软件你就不要挑了,不收费让你用,就知足吧。

  645. A:
  646.     一个新东西,而且没有压倒性的优势,想要流行开来,有盗版几乎是必须的

  647. ----------------------------------------------------
  648. 2016.07.14 20:20:22
  649. Q:
  650.     [图片]
  651.     Amber14+Ambertools15 尝试做有机小分子的动力学模拟,运行 tleap -s -f leaprc.ff10 时提示找不到 leaprc.ff10 文件,搜索找到其在oldff中如上,是否需要将它从 oldff 中提到上一层目录?或选用哪个新的参数文件比较合适?

  652. A:
  653.     14SB

  654. ----------------------------------------------------
  655. 2016.07.14 20:48:00
  656. Q:
  657.      老师,刚刚问您mayer 键级和laplacian键级都算出来放一个表里,老师您说有意义,那意义是否在于:两种不同定义的键级都表明了键级(强)的变化趋势?

  658. A:
  659.      键级特征不同,侧重点不同。看拉普拉斯键级原文的讨论

  660. ----------------------------------------------------
  661. 2016.07.14 22:10:26
  662. Q:
  663.     Counterpoise corrected energy =   -1164.078839135480
  664.     BSSE energy =       0.000659415350
  665.     sum of monomers =   -1164.075358838885
  666.     complexation energy =      -2.60 kcal/mole (raw)
  667.     complexation energy =      -2.18 kcal/mole (corrected)    这里的配位能表示什么? @Sobereva

  668. A:
  669.     相互作用能

  670. Q:
  671.     是两分子之间的?还是

  672. A:
  673.     当然是你设的两个片段之间的
  674.        
  675. Q:
  676.     Counterpoise corrected energy =   -1164.078839135480
  677.     BSSE energy =       0.000659415350
  678.     sum of monomers =   -1164.075358838885
  679.     complexation energy =      -2.60 kcal/mole (raw)
  680.     complexation energy =      -2.18 kcal/mole (corrected)                                        老师,那两个分子的总能量是哪一项??  谢谢老师

  681. A:
  682.     Counterpoise corrected energy
  683.        
  684. ----------------------------------------------------
  685. 2016.07.14 22:22:59
  686. Q:
  687.     大家有知道怎样用gaussian能够得到吸收光谱图么?
  688.     类似这样的[图片]

  689. A:
  690.      使用Multiwfn绘制红外、拉曼、UV-Vis、ECD和VCD光谱图
  691.     http://sobereva.com/224

  692. Q:
  693.     好的老师,我用的关键词是:# freq td=singlets b3lyp/3-21g 只能得到紫外可见光谱。吸收光谱怎样到处啊?

  694. A:
  695.     Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法
  696.     http://sobereva.com/314
  697.     紫外可见光谱本身就是吸收光谱的一种

  698. Q:
  699.     老师,是不是我要得到吸收光谱,是不是需要分别对体系的基态和激发态进行计算,然后得到他们的差值?

  700. A:
  701.     按博文去做就完了

  702. ----------------------------------------------------
  703. 2016.07.14 23:36:30
  704. Q:
  705.       老师 您好 通过紫外计算 怎么判断吸收光谱是哪个基团产生的 谢谢

  706. A:
  707.      图解电子激发的分类
  708.     http://sobereva.com/284
  709.    
  710.     做NTO、electron-hole分析等方式可考察,看博文

  711. ----------------------------------------------------
  712. 2016.07.14 23:50:15
  713. Q:
  714.     sob老师,请问GROMACS的培训什么时候能开始呢?

  715. A:
  716.      下学期

  717. ----------------------------------------------------
  718. 2016.07.15 01:52:19
  719. Q:
  720.     请问,Gaussian的NTO分析结果,除了最高NTO对儿的本征值之外,后面的布居分析结果(mulliken电荷,偶极矩等)是不是没么子意义呀?

  721. A:
  722.      对

  723. ----------------------------------------------------
  724. 2016.07.15 04:25:24
  725. Q:   
  726.         我觉得如果要得到激发前后的电荷密度差,激发态计算应该写上td density,才能得到激发态的电子密度fchk。然后用这个fchk减去基态的fchk,
  727.     或者减去td不写density的fchk
  728.     那算不出来?

  729. A:
  730.     光用td density得到的fch里存的还是基态轨道,必须按照multiwfn手册第四节开头的说明把自然轨道转存到fch,Multiwfn里给出的才是激发态密度。
  731.    
  732.     用.wfn就省事多了,直接用td density产生的.wfn文件在Multiwfn里面得到的就是激发态密度
  733.        
  734. Q:
  735.     对了,老师我fch文件是按照说明书写得出的,可不可以利用multiwfn像wfn一样得出电荷转移量?

  736. A:
  737.     可以
  738.        
  739. Q:
  740.     前提是已经得到了fchk文件,这是其输入文件[图片]
  741.     ,想得到第一激发态到基态的电荷转移量,用fchk与out,请问老师可以吗?

  742. A:
  743.     具体怎么定义转移量?

  744. Q:
  745.     从正电荷中心到负电荷中心,设得是0.0004,密度差的正值部分积分,怪我没说清楚....

  746. A:
  747.     这种输入文件不行。直接
  748.     cam-b3lyp/genecp td density out=wfn
  749.     用产生的wfn文件和基态的wfn文件求差就完了
  750.     非要用激发态的fch文件,如前所述必须把激发态自然轨道转存到fch里,否则fch里的轨道还是基态DFT轨道

  751. Q:
  752.     好的,谢谢老师了
  753.     明天我直接用wfn,再算一遍
  754.     [图片]
  755.     定义完转移量后,电荷转移量直接就是这个multiwfn数据对吧?

  756. A:
  757.     “定义完转移量后”意义不明
  758.         具体怎么定义转移量?
  759.     激发态与基态的密度差的正值或负值部分积分?
  760.        
  761. Q:
  762.     好的,谢谢老师了
  763.     明天我直接用wfn,再算一遍
  764.     [图片]
  765.     定义完转移量后,电荷转移量直接就是这个multiwfn数据对吧?
  766.        
  767. A:
  768.     就取这个
  769.        
  770. ----------------------------------------------------
  771. 2016.07.15 04:25:24
  772. Q:

  773.     iop(6/8=3)输出的是个方阵,multiwfn绘图的时候也是方阵。但是nto分析里却是用Nocc*Nvirtual的来做奇异值分解...为什么应用中有这么个差异呢?
  774.     尤其想问,如果拿iop6/8=3的方阵去做特征值分解,和nto结果会一样吗?
  775.     integral of hole是电子空穴分析的结果啊
  776.     老师来了[表情] 原来这样啊
  777.     sob老师麻烦你看看我刚才的问题啊

  778. A:
  779.     啥问题?

  780. Q:
  781.     iop(6/8=3)输出的是个方阵,multiwfn绘图的时候也是方阵。但是nto分析里却是用Nocc*Nvirtual的来做奇异值分解...为什么应用中有这么个差异呢?
  782.     尤其想问,如果拿iop6/8=3的方阵去做特征值分解,和nto结果会一样吗?
  783.     [表情]
  784.     [图片]请问用IOp(6/8=3)输出的这个density matrix是这个博文里说的跃迁密度矩阵吗?好像维度对不上

  785. A:
  786.     NTO本来就不是直接简单地对跃迁密度矩阵做个对角化
  787.     跃迁密度矩阵就是Nbasis*Nbasis维的

  788. Q:
  789.     还是有点糊涂..那写成Nocc*Nvirtual的好处是什么呢?直接对角化那个方阵会是什么情形?

  790. A:
  791.     直接对角化大抵什么也不是,需要做特殊的变换,有兴趣看原文琢磨

  792. ----------------------------------------------------
  793. 2016.07.15 04:30:25
  794. Q:
  795.     Sob老师您好,我想问一下两根键的electron density ellipicity along bond path不同,我该从哪方面去分析这种不同呢?因为我发现拉普拉斯算符沿着这个相同的键径变化趋势是几乎相同的。

  796. A:
  797.      键的pi特征不同。pi键特征越强,椭率越大(只有一个pi键的情况,两个pi键还是会接近轴对称的)

  798. Q:
  799.     谢谢老师。老师,我想问一下,能否接着研究造成这个Pi键特征不同的原因在哪?可以从哪方面研究?老师可否提点一下。

  800. A:
  801.     不结合实际问题很难说
  802.     先说清楚前提,什么过程的

  803. ----------------------------------------------------
  804. 2016.07.15 11:22:17
  805. Q:
  806.     请教下6-311G,3-21G这种split valence基组,例如2s分裂成2s'和2s'',被分裂的基函数是不是在线性组合成分子轨道时,2s'与2s''都参与,而不是通过2s=a*2s'+b*2s'', 用2s来参与分子轨道的

  807. A:
  808.     对

  809. ----------------------------------------------------
  810. 2016.07.15 17:01:12
  811. Q:
  812.     [图片][图片]请问大家,当前设置的半宽高是0.46ev, 我应该怎么设置把峰的最大值降下来
  813.     [图片]请问,怎么导出图右边的数据?
  814.     这是Multiwfn展的,貌似只能导出左边

  815. A:
  816.      spectrum_line.txt绘制出来的离散竖线就是对应右边的轴。选项-1也可以把每个跃迁的信息都输出出来,从屏幕上直接拷就完了

  817. Q:
  818.     老师我的意思是曲线上的值,以能量为横坐标,振子强度为纵坐标,导出一组这种组合方式的点@Sobereva

  819. A:
  820.      选项-1直接就输出了

  821. Q:
  822.     [图片]-1输出的是这个,我是想要(nm,f)以能量为横坐标,振子强度为纵坐标,这样的一系列的很多很多的点

  823. A:
  824.     根本就没有那样的东西,振子强度是对于一个跃迁来说的,哪可能有很多点
  825.     要么离散的竖线,要么连续的曲线

  826. Q:
  827.     对对,我就是想要连续的曲线,(nm,f)以能量为横坐标,振子强度为纵坐标,导出很多很多点就可以做连续曲线了

  828. A:
  829.     连续曲线的纵坐标是epsilon,不是振子强度!
  830.     仔细看此文了解是怎么从量化算的数据展宽成曲线的
  831.     使用Multiwfn绘制红外、拉曼、UV-Vis、ECD和VCD光谱图
  832.     http://sobereva.com/224
  833.     (nm,epsilon):连续曲线
  834.     (nm,f):离散的值

  835. ----------------------------------------------------
  836. 2016.07.15 17:52:01
  837. Q:
  838.     Sob老师,您好,请问在做电荷分析时,使用gaussian做NBO分析和使用VASP做bader charge分析,哪个会更好一些,假设分析的体系使用gaussian 和VASP均可以进行优化。谢谢sob老师@Sobereva

  839. A:
  840.      什么体系?分子体系当然不要用VASP,周期性体系想用Gaussian也很费劲

  841. Q:
  842.     谢谢Sob老师,我计算的体系是一个分子体系,在进行结构优化时均是在gaussian中进行的,可是在使用NBO分析电荷时,我做的体系总会出现和理论值相悖 的结果,同时又试了mulliken和APT分析,结果都和NBO分析的一样,可是使用bader charge分析时却可以很好的吻合~我该怎么选择呢?再次谢谢sob老师~

  843. A:
  844.      尝试ADCH或Hirshfeld电荷,看群共享《原子电荷计算方法的对比》
  845.     所谓bader charge就是AIM电荷,一般都叫AIM电荷,要计算AIM电荷甭用其它程序,用Multiwfn(看手册4.17.1节的例子)或者AIMALL
  846.     VASP+Bader程序算出来严重不准

  847. ----------------------------------------------------
  848. 2016.07.15 18:22:24
  849. Q:
  850.     各位老师好,为什么我安装GV508W是出现
  851.     Failed tolocate .ini file in Gaussian directory

  852. A:
  853.     全都卸干净,先安装高斯,再把gview装到相同目录
  854.     目录不得有中文、空格

  855. ----------------------------------------------------
  856. 2016.07.15 18:30:52
  857. Q:
  858.     [表情]
  859.     能不能录个视频呢,淘宝上卖,更多人可以学习

  860. A:
  861.     不能

  862. ----------------------------------------------------
  863. 2016.07.15 18:35:03
  864. Q:
  865.     sob老师,给体和受体的质心距离越大是不是可以说明它的CT态更多一些

  866. A:
  867.     没法说“CT态更多”,语焉不详

  868. ----------------------------------------------------
  869. 2016.07.15 18:35:10
  870. Q:
  871.     老师给推荐本入门dft化学书吧

  872. A:
  873.     levine的quantum chemistry

  874. Q:
  875.     好的,谢谢,那本书没看完,以前没注意到有dft

  876. A:
  877.     起码把目录看一遍啊

  878. ----------------------------------------------------
  879. 2016.07.15 18:50:45
  880. Q:
  881.     sob老师请问下计算荧光时候,opt正常结束了,但是计算频率还是有虚频,这个是不是应该重新优化?

  882. A:
  883.     看距离收敛限还差多少
  884.     Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法
  885.     http://sobereva.com/278

  886. Q:
  887.     好的,是不是如果差的少就可以忽略?

  888. A:
  889.      y

  890. ----------------------------------------------------
  891. 2016.07.15 18:56:03
  892. Q:
  893.     [图片]
  894.     请教一下各位老师,优化一个这样的配合物并计算他的NBO用什么基组比较好
  895.     高斯小白,希望大家不要笑话~

  896. A:
  897.     6-31G*

  898. Q:
  899.     不是说有F这样电负性大的原子要加弥散函数吗?

  900. A:
  901.     根本不用!
  902.     至多3-zeta也已经描述得很好

  903. Q:
  904.     好的,谢谢老师,我之前用的6-31+g*优化的结果只能得到两个yes是什么原因呢

  905. A:
  906.     弥散函数容易有碍收敛
  907.     量子化学计算中帮助几何优化收敛的常用方法
  908.     http://sobereva.com/164

  909. ----------------------------------------------------
  910. 2016.07.16 00:12:12
  911. Q:
  912.      老师,请问用Multiwfn做AdNDP分析多中心键,如果共86个价电子,做完分析,剩8.16个电子,是否合理?

  913. A:
  914.     有的体系整体离域性强,注定会剩下很多电子不好归属
  915.     没法光看这个判断

  916. ----------------------------------------------------
  917. 2016.07.16 00:13:00
  918. Q:
  919.     请问老师,IRC图上的total energy指的是什么能量,跟zero-point energy有关系吗

  920. A:
  921.     就是电子能量
  922.     无关

  923. ----------------------------------------------------
  924. 2016.07.16 01:22:59
  925. Q:
  926.     请问各位,哪里有multiwfn程序下载及help?

  927. A:
  928.     [图片]

  929. Q:
  930.     [表情]

  931. A:
  932.     [图片]

  933. ----------------------------------------------------
  934. 2016.07.16 03:35:51
  935. Q:
  936.     --------Markdown Markdown--------

  937. A:
  938.     --------Markdown Markdown--------
  939.     公告:本月16~19日群主Sobereva因为讲科音基础量化培训班的课十分忙碌,这期间不回复群里任何学术方面的提问,私聊更不回,请相互转告。

  940. Q:
  941.     今天是最后一天啊
  942.     抓紧时间问

  943. A:
  944.     因为要备课,今天晚上也几乎不回
  945.        
  946. ----------------------------------------------------
  947. 2016.07.16 03:36:46
  948. Q:
  949.     老师我们一般研究NBO时大多看LP (2) F 21→  BD*(  1) C  14 - F  20  这种孤对电子向反键轨道转移电子的二阶稳定化能。但是我现在研究的体系,要断裂的键没有出现类似情况只有两个反键轨道的二阶稳定化能 BD*(   1) S  10 - O  13 →BD*(   1) O  13 - C  15 ,(这里O13-C15键断裂)@Sobereva

  950. A:
  951.     没意义。甭试图用NBO讨论断键问题
  952.     较新版本NBO应该都不会输出BD*->BD*这种无意义的鬼东西了

  953. Q:
  954.     那断键成新键用什么手段比较好呢,键级什么的都用了@Sobereva

  955. A:
  956.     ELF、键级、自旋密度等等都可以讨论

  957. ----------------------------------------------------
  958. 2016.07.16 04:48:11
  959. Q:
  960.      老师您好,分析弱作用时,用RDG函数,分析出来需要对比的几个做出来都是绿色的,不好进行对比。请教一下老师在,这样的情况我应该如何处理呢?谢谢老师

  961. A:
  962.     说明本身作用强度就相仿佛

  963. ----------------------------------------------------
  964. 2016.07.18 05:10:27
  965. Q:
  966.     请教老师一个问题,ONIOM方法中,MM区对QM区的极化作用,非键合原子之间也能体现出来吗?

  967. A:
  968.     用电子嵌入就可以体现

  969. ----------------------------------------------------
  970. 2016.07.18 06:26:35
  971. Q:
  972.     请问,使用内存为16G的计算机在windows系统下进行计算。使用关键词%mem=400mw,按理说内存应该够用,为什么还是出现2070内存不够错误。

  973. A:
  974.     因为你用的高斯是32bit的

  975. ----------------------------------------------------
  976. 2016.07.19 04:23:07
  977. Q:
  978.     大家好,最近被含时的密度泛函所干扰,请问可以用gauss计算么?还有计算的时候命令怎么写啊,自学中,求指教

  979. A:
  980.     Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法
  981.     http://sobereva.com/314

  982. ----------------------------------------------------
  983. 2016.07.19 07:00:40
  984. Q:
  985.    
  986.     请教一个问题, 安装gromacs-5.1.2, 命令如下
  987.     cmake .. -DGMX_BUILD_OWN_FFTW=ON -DCMAKE_INSTALL_PREFIX=/home/x_qinli/Programs/Gromacs512/
  988.     make -j 8 ; make install
  989.    
  990.     安装一切顺利, 只是在文件夹bin中只有如下几个文件:
  991.    
  992.     demux.pl  gmx  gmx-completion.bash  gmx-completion-gmx.bash  GMXRC  GMXRC.bash  GMXRC.csh  GMXRC.zsh  xplor2gmx.pl
  993.    
  994.     其它的可执行工具都没有, 试过多次, 还是如此! 而gromacs-4.6.7就可以顺利安装, 请问有人知道如何解决这个问题吗?  十分感谢!

  995. A:
  996.     本来就这样

  997. ----------------------------------------------------
  998. 2016.07.20 16:50:53
  999. Q:
  1000.     请教一下:我想将基态优化和激发态计算通过一个gjf进行计算。
  1001.     设置是这样的:
  1002.     %chk=DMAC-DPS-link.chk
  1003.     %mem=20GB
  1004.     %nprocshared=24
  1005.     # opt freq b3lyp/6-31g(d)
  1006.    
  1007.     s0 opt
  1008.    
  1009.     0 1
  1010.     --link1--
  1011.     %oldchk=DMAC-DPS-link.chk
  1012.     %chk=DMAC-DPS-link-S1-T1.chk
  1013.     %mem=20GB
  1014.     %nprocshared=24
  1015.     #p td=(50-50,nstates=6) b3lyp/6-31g(d) geom=check guess=read
  1016.    
  1017.     S1 T1 calculation
  1018.    
  1019.     0 1
  1020.     但是第一步的基态优化成功,接着进行第二步计算的时候出现下面的语法错误提示:%oldchk=DMAC-DPS-link.chk
  1021.     Error parsing percent card:
  1022.     QPErr --- A syntax error was detected in the input line.
  1023.     %oldchk=DMAC-DPS-link.chk
  1024.     '
  1025.     Last state=""
  1026.     TCursr=  184 LCursr=    1
  1027.     Error termination via Lnk1e in /home/customer/g09/l1.exe at Thu Jun  9

  1028. A:
  1029.      老的g09版本不支持oldchk写法

  1030. Q:
  1031.     那sober老师我直接这么写可不可以?--link1--
  1032.     %chk=DMAC-DPS-link.chk
  1033.     %mem=20GB
  1034.     %nprocshared=24
  1035.     #p td=(50-50,nstates=6) b3lyp/6-31g(d) geom=check guess=read
  1036.     @Sobereva

  1037. A:
  1038.     可
  1039.        
  1040. Q:
  1041.     谢谢老师计算没有问题 ,请问这样计算的.chk文件是否保存了S0 opt+freq的信息和S1/T1的信息,还是只保存了S1/T1的计算信息,谢谢老师

  1042. A:
  1043.     具体指什么信息?

  1044. Q:
  1045.     具体指的是S0的HOMO/LUMO信息@Sobereva

  1046. A:
  1047.     你做td=(50-50,nstates=6) b3lyp/6-31g(d)计算本来就会产生 b3lyp/6-31g(d)的轨道并保存在chk里
  1048.     这就是S0的
  1049.     TDDFT没有所谓激发态轨道
  1050.        
  1051. ----------------------------------------------------
  1052. 2016.07.20 17:47:34
  1053. Q:
  1054.      请问老师Multiwfn 可以做出来这样的图吗[图片]

  1055. A:
  1056.      使用Multiwfn观看分子轨道
  1057.     http://sobereva.com/269

  1058. Q:
  1059.     老师,可以用多功能波函数做出二聚体的HOMO图吗?

  1060. A:
  1061.     当然可以
  1062.     那叫Multiwfn不叫多功能波函数

  1063. Q:
  1064.     那老师该怎么看呢

  1065. A:
  1066.     Sobereva(190258442)  9:48:36
  1067.      使用Multiwfn观看分子轨道
  1068.     http://sobereva.com/269
  1069.     Sobereva(190258442)  9:48:36
  1070.      使用Multiwfn观看分子轨道
  1071.     http://sobereva.com/269
  1072.     Sobereva(190258442)  9:48:36
  1073.      使用Multiwfn观看分子轨道
  1074.     http://sobereva.com/269

  1075. ----------------------------------------------------
  1076. 2016.07.20 18:12:08
  1077. Q:
  1078.     我计算了一个Ir配合物,一直不收敛,尝试了Sobereva老师解决SCF不收敛帖子中所有方法。 哪位老师有空可以帮我看一下我的输入文件?
  1079.     [图片]

  1080. A:
  1081.     genecp啊
  1082.     否则结构不塌到一起才怪

  1083. Q:
  1084.     试过 genecp

  1085. A:
  1086.     不是试,是必须用!!!!!!!!!!!!!!1
  1087.     否则即便碰巧收敛了结果又有何意义!?

  1088. Q:
  1089.     [表情]第一次见社长大人用这么多叹号

  1090. A:
  1091.     赝势基组哪能不用赝势?

  1092. Q:
  1093.     最开始的时候 勇的是 genecp ,看到有人没有加 ecp 也收敛了
  1094.     所以我把所有的能找到的资料全部尝试了一边

  1095. A:
  1096.     那个人算出来的结果是翔

  1097. Q:
  1098.     全部没有收敛
  1099.     学到了~
  1100.     可以我用genecp 也没有收敛, 您帖子上的十几种方法都用过了

  1101. A:
  1102.     博文里的方式你用的方式不对,或者结构、自旋多重度之类不合理

  1103. Q:
  1104.      也就是说控制行 修改 genecp 后,在保证结构正确都情况下,尝试您帖子上的方法 有很大可能可以解决收敛问题,对吧?

  1105. A:
  1106.      对

  1107. ----------------------------------------------------
  1108. 2016.07.20 18:25:21
  1109. Q:
  1110.     请问,有没有办法对cube文件进行插值。有没有这个程序?计算出平滑的cube太耗时了

  1111. A:
  1112.     记得网上有不少现成的spline三维插值的C语言代码可以用

  1113. ----------------------------------------------------
  1114. 2016.07.20 18:42:55
  1115. Q:
  1116.     [图片]各位老师好,我想问下这些5f(0),5f(c1)等这些轨道具体代表哪些轨道?比如2px,2py这样的,在wfn手册没找到,谢谢

  1117. A:
  1118.     这是NBO程序里用的表示符号,这里说了
  1119.     球谐型与笛卡尔型Gauss函数的转换关系
  1120.     http://sobereva.com/97
  1121.     [图片]

  1122. ----------------------------------------------------
  1123. 2016.07.20 18:47:37
  1124. Q:
  1125.     TDDFT计算是否不能给出n-π*跃迁?我做了苯甲醛的计算,只有发现有π-π*的跃迁。求大侠指点

  1126. A:
  1127.     当然能给出
  1128.     算的态数不够,或者没看仔细

  1129. ----------------------------------------------------
  1130. 2016.07.20 18:51:19
  1131. Q:
  1132.      大社长,纯B团簇吸附碱金属,用PBE0是不是比B3LYP要好?此前用B3LYP优化无虚频,但是PBE0就有虚频

  1133. A:
  1134.      没听说,通常都是半斤八两

  1135. Q:
  1136.     [图片]
  1137.     这句话应该可以认为PBE0好于B3LYP?

  1138. A:
  1139.     >=的关系

  1140. ----------------------------------------------------
  1141. 2016.07.20 18:51:39
  1142. Q:
  1143.     请问下,高斯里,设置电场时,如field=z+400.。。。优化过程中会对原子的坐标有限制吗

  1144. A:
  1145.     没限制

  1146. ----------------------------------------------------
  1147. 2016.07.20 19:05:34
  1148. Q:
  1149.      请问可否通过multiwfn得到前线轨道的具体轨道组成?比如主要是pi还是pi*等?

  1150. A:
  1151.     一看轨道图形就知道了

  1152. ----------------------------------------------------
  1153. 2016.07.20 19:06:02
  1154. Q:
  1155.      大群主我开始追LL了

  1156. A:
  1157.     sunshine追了?

  1158. ----------------------------------------------------
  1159. 2016.07.20 19:08:22
  1160. Q:
  1161.     求助Sob老师,我们老师说让我尝试一下PBE-D和PBE0-D3方法计算我的体系?我查了PBE0就是PBE1PBE,那PBE是?它们的关键词这么输入正确吗?PBE-D:PBE/EmpiricalDispersion=GD3; PBE1PBE/EmpiricalDispersion=GD3BJ

  1162. A:
  1163.     PBE=PBEPBE

  1164. Q:
  1165.     哦哦,好的,明白了

  1166. A:
  1167.     PBEPBE/def2TZVP EmpiricalDispersion=GD3BJ

  1168. Q:
  1169.     老师,那-D和-D3的区别是在后面EmpiricalDispersion=GD3(BJ)是吗

  1170. A:
  1171.     DFT-D色散校正的使用
  1172.     http://sobereva.com/210

  1173. ----------------------------------------------------
  1174. 2016.07.20 19:10:24
  1175. Q:
  1176.     恩恩,谢谢,也就是如果multiwfn进行处理,如果要用到fchk文件,我将这个计算完成的fchk文件放入即可,不管是基态还是激发态。如果需要基态的log文件,我只需要将这个一次计算得到的log文件中的激发态计算信息删除就可了;如果需要激发态的log文件,我只需要将这个一次计算得到的log文件中的基态计算信息删除就可了?谢谢@Sobereva

  1177. A:
  1178.     什么叫基态信息?
  1179.     你不写TD算的就是基态

  1180. ----------------------------------------------------
  1181. 2016.07.20 19:18:08
  1182. A:
  1183.     这文章真够省事
  1184.     http://www.orientjchem.org/pdf/vol32no3/OJC_Vol32_No3_p_1589-1600.pdf
  1185.    
  1186.     直接把Multiwfn文本窗口截图贴上去当表格
  1187.     [图片]

  1188. Q:
  1189.     感觉他里面的很多图片都是截图的感觉

  1190. A:
  1191.     然
  1192.     伊朗人现在计算化学水文甚多

  1193. Q:
  1194.     是因为没条件做实验么
  1195.     差不多是
  1196.     国内也有很多这种情况,实验条件有限,就搞理论计算去了。

  1197. A:
  1198.     大抵和天朝一样,竞争激烈,不得不多灌水

  1199. Q:
  1200.     要是总灌水,参加个什么会,都不好意思见人。[表情]
  1201.     这个都看不清楚图 怎么会接收发表出来的呢

  1202. A:
  1203.     毕竟是个非SCI刊物

  1204. Q:
  1205.     奥,要求不高啊
  1206.     750镑一篇文章…
  1207.     就rsc advance
  1208.     这叫开放获取吧,开源什么鬼
  1209.     那收钱是什么意思

  1210. A:
  1211.     一般的期刊是向读者要钱,不向作者要钱
  1212.     免费观看的期刊是向作者要钱,不向读者要钱

  1213. Q:
  1214.     [表情]那谁还投这个期刊了?
  1215.     如果nature管你要两万块钱,给你发一篇你发不发

  1216. A:
  1217.     开放访问是期刊的趋势
  1218.     作者有经费,大多也不差这点钱。当然对民科是残忍的
  1219.     20万都照样有人发

  1220. Q:
  1221.     国内大学体制运作的好20万一篇nature主刊出来马上就回本了

  1222. A:
  1223.     自筹费用
  1224.     要么发别的

  1225. ----------------------------------------------------
  1226. 2016.07.20 19:19:10
  1227. Q:
  1228.     用高斯结构优化,虚频有的在40-50左右,还有个一百左右的。有的有两到三个虚频。。请问怎么调合适。试了一下那个start  animation,感觉不怎么管用。求助各位大神。

  1229. A:
  1230.     Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法
  1231.     http://sobereva.com/278

  1232. ----------------------------------------------------
复制代码

评分 Rate

参与人数
Participants 1
eV +8 收起 理由
Reason
sobereva + 8

查看全部评分 View all ratings

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-15 22:23 , Processed in 0.168283 second(s), 25 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list