计算化学公社

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

[群聊Q&A整理] 公社群聊Q&A整理:2015.06.21 20 ~ 2015.06.25 19 Mark

[复制链接 Copy URL]

3097

帖子

29

威望

1万

eV
积分
17098

Level 6 (一方通行)

本帖最后由 liyuanhe211 于 2017-3-6 16:09 编辑

本帖由 @978142355 整理

相关概况可阅读:http://bbs.keinsci.com/thread-2022-1-1.html
十分欢迎一起校读聊天记录片段并发布:http://bbs.keinsci.com/forum.php?mod=viewthread&tid=5208
用于检索聊天记录全文的小程序已发布:http://bbs.keinsci.com/forum.php?mod=viewthread&tid=5207

暂无全文附件

本次附件(附件内容与下文一致):
2015.06.21 20 ~ 2015.06.25 19 Mark.txt (59.68 KB, 下载次数 Times of downloads: 5)


  1. 2015.06.21 19:14:04
  2. Q:
  3.     请教一个水分子氢键的问题,文献上是这样定义的[图片],我想问一下这个angle 小于30度,是指角O-H-O的角吗?? 如果没表述清楚的化,那么请告知一下在编程中如何定义这个角! 谢谢

  4. A:
  5.     文中表达得不确切。应该是O-H---O大于150度,或者O-O-H小于30度

  6. ----------------------------------------------------
  7. 2015.06.21 20:34:35
  8. Q:
  9.     想请问下下图中的T1 对应Excited State   1:      Triplet-A      2.1985 eV(自旋多重度为1时)的值么? [图片]
  10.     这个T1 和直接用自旋多重度为3在单态下的计算有关系么?

  11. A:
  12.      only 0 1基态的SCF能量加上2.1985eV,以及你直接在0 3下得到的SCF能量,都对应于这个T1能量。

  13. Q:
  14.     好的。谢谢。
  15.     请问Sob,我做了一个简单的测试。分别计算了甲烷的三态的单点能,和单态下的三态激发能,分别用如下命令:#p  b3lyp/6-31g*;0,3 和#p td(triplet,root=1) b3lyp/6-31g*;0,1. 但是在三态下的能量是SCF Done:  E(UB3LYP) =  -40.0692977242,单态下T1的能量是Excited State   1:      Triplet-T,Total Energy, E(TD-HF/TD-KS) =  -40.0862074709。 这两个能量之间有关系么?

  16. A:
  17.     两种方法结果肯定不同,但原理上都对

  18. Q:
  19.     那哪个对应的才是T1的能量?

  20. A:
  21.     都是

  22. Q:
  23.     我用什么方法才能和上图中T1对比?

  24. A:
  25.     一个是delta-SCF方法算的,一个是TDDFT算的。哪个都能对比

  26. Q:
  27.     好的。那我就用TDDFT的结果对比吧。谢谢了。
  28.     在做激发态三态计算时,优化前几步激发能都是正值,可是优化了一段时间后激发能就变成了负值,这个问题该怎么客服?

  29. A:
  30.     输入输出文件都给他就成了
  31.     看错了,以为你想说“这个问题该怎么问客服”
  32.     用出现负值的结构,改变多重度设定看看激发能是否为正。即原来的基态自旋多重度在当前结构下可能已经变了。

  33. ----------------------------------------------------
  34. 2015.06.21 20:41:52
  35. Q:
  36.      S老师,amber软件里
  37.     伞形取样法计算的是自由能面(基于反应坐标)上的最低能量路径?
  38.     而NEB方法计算的是势能面(基于反应坐标)的最低能量路径?
  39.     是这样吗?

  40. A:
  41.     NEB是那样。但是伞状取样并不是得到自由能面的最低能量路径,应该说是得到给定的一个路径上的自由能变化。

  42. Q:
  43.     另有一问,伞形取样法得到的路径上某点的自由能,应该是这点的最低自由能吧?

  44. A:
  45.     没法这么说

  46. ----------------------------------------------------
  47. 2015.06.21 21:42:05
  48. Q:
  49.     不好意思打错字了。我没试过改变多重度。我是想优化T1下的稳定构型,可以改变多重度么?

  50. A:
  51.     你就用0 3优化

  52. ----------------------------------------------------
  53. 2015.06.21 22:05:30
  54. Q:
  55.     [图片]
  56.     请问,我要把CH3CH3改成ch3ch2d
  57.     应该怎么改这个输入文件?

  58. A:
  59.     最后一个H改为H(iso=2)

  60. Q:
  61.     [图片]
  62.     是这样么?

  63. A:
  64.      对

  65. Q:
  66.     嗯,懂了,谢谢大神
  67.     为什么CH3CH2D和CH3CH3优化出来的所有键长键角的结果都是一样的啊?

  68. A:
  69.     因为优化过程本来就不依赖于原子质量。只看核电荷、电子数

  70. Q:
  71.     哦,懂了,我查了一下频率,是不同的

  72. A:
  73.     只有振动平均结构能反映出同位素质量的影响

  74. Q:
  75.     了解,谢谢大神

  76. A:
  77.     看样子真是进入认真写论文的阶段了

  78. Q:
  79.     是哒,写得差不多了,分块写的,现在再拼一个完整的文章~
  80.     图和表格也差不多准备完了

  81. A:
  82.     頑張れ!

  83. Q:
  84.     好期待九月份博士开学~
  85.     谢谢[图片]
  86.     [Emoticon]打算在科研的道路上一直走下去

  87. A:
  88.     刮目相看了

  89. ----------------------------------------------------
  90. 2015.06.21 22:08:05
  91. Q:
  92.     sob老师,Multiwfn可以做这种图么
  93.     [图片]

  94. A:
  95.     不能

  96. Q:
  97.     嗯   这个我看介绍说是等值面填色图

  98. A:
  99.     结合VMD能作

  100. Q:
  101.     不知道做这个图有啥意义

  102. A:
  103.     使用Multiwfn结合VMD分析和绘制分子表面静电势分布
  104.     http://sobereva.com/196
  105.     在GaussView里绘制分子表面静电势填色图的过程
  106.     http://sobereva.com/253
  107.     幻灯片:Predicting reactive sites (反应位点的预测)
  108.     http://sobereva.com/234

  109. ----------------------------------------------------
  110. 2015.06.21 22:09:30
  111. Q:
  112.     嗯  我看看老师发的地址
  113.     上面那个图是用这个软件做的  Moliso  看起来也很简单
  114.     http://www.moliso.de/

  115. A:
  116.     真是稀奇

  117. Q:
  118.     老师也可以看看

  119. A:
  120.     moliso以前玩过,并不好用

  121. ----------------------------------------------------
  122. 2015.06.21 22:57:59
  123. Q:
  124.     As I explained in my original comments: unless the
  125.     same basis is used for PCA, the pictures of the PCA projections cannot
  126.     directly be compared by visual inspection.一个审稿人两次提到了这个问题,这是什么意思呢
  127.     针对pca图提出来的

  128. A:
  129.     除非PCA的PC的定义相同(产生PC的过程一致),否则PCA图没法直接对比

  130. Q:
  131.     我做的是针对一个系统突变前后分析的主成分
  132.     这不存在pc产生的不一致吧

  133. A:
  134.      那你应该把突变后的结果往突变前产生的PC上投影

  135. ----------------------------------------------------
  136. 2015.06.21 23:00:18
  137. Q:
  138.      记得卖萌君是男人么,怎么头像是这么萌的菇凉
  139.     玛丽卖萌是女的
  140.     [Emoticon]
  141.     可是上次咋看到卖萌君的男人照了呢
  142.     [Emoticon]
  143.     那卖萌君确实好菇凉,有对象了吗[Emoticon]
  144.     阿萌说还没有

  145. A:
  146.     有了

  147. Q:
  148.      你骗娃娃
  149.     [图片]
  150.     我想你们都疯了
  151.     [Emoticon]
  152.     竟然有了
  153.     我只想安静的写论文
  154.     大家排队
  155.     不要争抢
  156.     我们可以给你写论文(⊙o⊙)哦
  157.     都要做有素质的竞争者
  158.     明明我排在第一位
  159.     你们这么说,我家妹子会吃醋的

  160. A:
  161.     散了散了,别影响人家用功

  162. ----------------------------------------------------
  163. 2015.06.22 01:03:14
  164. Q:
  165.     Sob老师你好……使用ADMP的时候,输出文件有上万步……这种情况下用GaussView打开很慢,请问有什么好的办法吗?

  166. A:
  167.     将Gaussian从头算动力学轨迹转换为pdb轨迹的程序GauMDconv
  168.     http://sobereva.com/126
  169.    
  170.     转换完了你可以用VMD看,快得多

  171. ----------------------------------------------------
  172. 2015.06.22 05:44:35
  173. Q:
  174.     com输入格式正确和log正常结束也可以正常打开 但fchk却打不开[Emoticon]  是计算出错了?[Emoticon]
  175.     [图片]
  176.     [图片]
  177.     那样算不了 分子轨道了[Emoticon]
  178.     [图片]这个是输入文件  有什么关键字冲突吗?

  179. A:
  180.      #P前面都没换行

  181. ----------------------------------------------------
  182. 2015.06.22 16:34:07
  183. Q:
  184.     北京科音自然科学研究中心(www.keinsci.com)将于2015年8月于北京开办以下培训班:
  185.     (1) 8月5日至7日:2015北京科音量子化学基础暑期班,通知见
  186.     http://www.keinsci.com/workshop/2015SKBQC_ann1.html
  187.     (2) 8月10日至13日:第二届量子化学波函数分析暨Multiwfn程序使用培训班,通知见
  188.     http://www.keinsci.com/workshop/2WFN_ann1.html
  189.     对于量子化学初学者,两个培训班同时参加不仅能充分学习基础量子化学知识、掌握主流量子化学程序的使用,在学习波函数分析和Multiwfn使用后,还将使得量子化学研究水准明显进一步提升、将问题分析得加更深入透彻、令文章内容更为丰富充实。两个培训班同时参加可享受较多优惠,优惠信息见http://www.keinsci.com/workshop/2015SKBQC+2WFN.html。
  190.    
  191.     培训的全部课程将由本群群主之一、Sobereva网名主要使用者、Multiwfn波函数分析程序的开发者卢天博士讲授。

  192. A:
  193.     欢迎参加

  194. ----------------------------------------------------
  195. 2015.06.22 16:45:26
  196. Q:
  197.     [图片]gamess运行过程中产生的这些文件可以删除吗?

  198. A:
  199.     算完了再删

  200. Q:
  201.     [图片]sob老师,磁盘空间快被占满了,这种情况下,怎么解决呢?

  202. A:
  203.     腾空间呗

  204. Q:
  205.     就是上述的两个很大的文件占住了空间

  206. A:
  207.     如果硬盘都占满了还算不完,那就只能尝试用别的计算方法、设定

  208. Q:
  209.     哦 这样啊  谢谢sob老师
  210.     那你怎么解决的?

  211. A:
  212.     计算级别的选择要量力而行

  213. Q:
  214.     6-31g*  估计是体系太大了  200多个原子呢  看来确实不能算能量分解了

  215. A:
  216.     DIRSCF=.T.你用了么?

  217. Q:
  218.     在另一个体系中用了,但是算一会儿就会自动停止  不清楚为什么

  219. A:
  220.     看输出报错。如果程序把GTF积分储存到硬盘上,自然会大得很

  221. Q:
  222.     [图片]最后没有报错信息

  223. A:
  224.     Direct SCF既不会占多少内存也不会占多少硬盘

  225. Q:
  226.     手动删除硬盘上的GTF积分可行吗?  恩

  227. A:
  228.     必须用Direct SCF,这是不占硬盘的唯一办法
  229.     其它出错多尝试bypass

  230. ----------------------------------------------------
  231. 2015.06.22 17:14:34
  232. Q:
  233.     我之前也遇到这样的问题,算个qcisd(t)/tz结果,rwf写满了磁盘,算不下去了
  234.     好像劈裂了rwf也没用,因为分开的小文件如果删了其中一个,计算也就直接死掉了
  235.     算不下去啊~

  236. A:
  237.     计算级别的选择要量力而行
  238.     DIRSCF=.T.你用了么?
  239.     多尝试。
  240.     显然

  241. ----------------------------------------------------
  242. 2015.06.22 17:28:01
  243. Q:
  244.     求教各位  把MO6优化好的构型保存好,用M06/aug-cc-pVDZ做单点计算,直接报scf不收敛 这个是不是可以把scf收敛精度降下来

  245. A:
  246.     int=ultrafine scf=novaracc之类没辙再降收敛限

  247. ----------------------------------------------------
  248. 2015.06.22 17:35:11
  249. Q:
  250.     请教,使用在DFT中,选择GGA+U计算 内包金属富勒烯 这类分子可以吗,会有明显的缺陷吗?

  251. A:
  252.     这种体系就拿普通杂化泛函算就行了
  253.    
  254.    
  255.     达夫(707476975)  9:35:11
  256.     请教,使用在DFT中,选择GGA+U计算 内包金属富勒烯 这类分子可以吗,会有明显的缺陷吗?

  257. ----------------------------------------------------
  258. 2015.06.22 17:38:38
  259. Q:
  260.     gamess,lmoeda用b3lyp计算。按顺序会经历那些过程呢?[图片]这个又是进行到哪一步了呢?  如果有较小的.log文件,能否发来一个参考一下,谢谢啦~~

  261. A:
  262.     自身就带了例子文件

  263. ----------------------------------------------------
  264. 2015.06.22 18:30:02
  265. Q:
  266.     gamess中,做能量分解时,用dfttyp=b3lyp是不是不能再加scftyp=rhf?

  267. A:
  268.     不矛盾

  269. Q:
  270.     加不加scftyp=rhf看什么呢
  271.     看体系是闭合还是开合?

  272. A:
  273.     闭壳层都写这个,开壳层写uhf或rohf
  274.     这设定的是SCF是什么类型,和泛函选择是两回事

  275. ----------------------------------------------------
  276. 2015.06.22 21:22:36
  277. Q:
  278.     [图片]轨道分析出现这样的错误[Emoticon] 求指点[Emoticon]

  279. A:
  280.     用Multiwfn做

  281. ----------------------------------------------------
  282. 2015.06.22 21:47:37
  283. Q:
  284.      老师 我最近算了一个ts 加上了tight calcall得到的振动也凑合 但是做IRC时用stepsize=15,取点=10只能出5个点 我加大stepsize=20 取点=10,加上了LQA,但是得到的只有11个点 图形如下[图片]这该怎么办呢

  285. A:
  286.     再加上calcall
  287.     IRC里加

  288. ----------------------------------------------------
  289. 2015.06.23 05:24:52
  290. Q:
  291.      请问用Multiwfn可以得到金属团簇的势能或平均势能的径向分布函数吗?

  292. A:
  293.       不能

  294. ----------------------------------------------------
  295. 2015.06.23 05:30:53
  296. Q:
  297.     什么键不影响结果,不画都行。
  298.     请问Multiwfn能输出HF能量吗?

  299. A:
  300.      不能

  301. ----------------------------------------------------
  302. 2015.06.23 06:15:34
  303. Q:
  304.     请教一个问题,把小分子放进酶的活性口袋里,小分子的前线轨道形状会不会发生改变? 从而导致小分子的活性位点也发生改变?

  305. A:
  306.      首先此时分子轨道就是整体的了。若把周围环境当成外场考虑,小分子的轨道肯定会发生变化,但毕竟是弱相互作用,变得不会那么显著。

  307. ----------------------------------------------------
  308. 2015.06.23 15:19:33
  309. Q:
  310.        今天使徒没有来呢 好伤心

  311. A:
  312.     ?

  313. ----------------------------------------------------
  314. 2015.06.23 15:22:20
  315. Q:
  316.     eva的设定是20156.22使徒来袭。。。。

  317. A:
  318.     还以为你想说UBW 24话使出咖喱棒之后那个爆炸很有EVA风格但是使徒没来呢

  319. ----------------------------------------------------
  320. 2015.06.23 16:44:00
  321. Q:
  322.     老师,求指导,NCI是指什么分析?multiwfn程序有相关分析功能吗?谢谢老师 卢天

  323. A:
  324.      使用Multiwfn图形化研究弱相互作用
  325.     http://sobereva.com/68
  326.     Multiwfn支持的弱相互作用的分析方法概览
  327.     http://sobereva.com/252
  328.     使用Multiwfn研究分子动力学中的弱相互作用
  329.     http://sobereva.com/186
  330.     弱相互作用的分析在8月份的Multiwfn培训班中会详细介绍和演示

  331. ----------------------------------------------------
  332. 2015.06.23 17:10:17
  333. Q:
  334.     [图片]
  335.     我算DFT-D3的输入文件,
  336.     提示书写有错误

  337. A:
  338.     dispersion后面多写一个i

  339. Q:
  340.     那个是修改时候Insert的i
  341.     我看了DFT-D3的帖子
  342.     [图片]

  343. A:
  344.     用的哪个g09版本?

  345. Q:
  346.     D.01

  347. A:
  348.      把图截完整

  349. Q:
  350.     不对,这个版本是B.01
  351.     应该是这个问题啦?

  352. A:
  353.     D.01才行

  354. Q:
  355.     D.01才加入GD3?

  356. A:
  357.     sure

  358. ----------------------------------------------------
  359. 2015.06.23 17:11:42
  360. Q:
  361.       老师,susceptibility该翻译成极化率还是磁化率呢?

  362. A:
  363.     磁化率

  364. ----------------------------------------------------
  365. 2015.06.23 17:14:22
  366. Q:
  367.     请问,如何看激发态优化是否收敛呢@sobereva

  368. A:
  369.     收敛了任务自然就正常结束了

  370. Q:
  371.     现在。log文件都一万多行了

  372. A:
  373.     [图片]
  374.     [图片]

  375. ----------------------------------------------------
  376. 2015.06.23 17:30:17
  377. Q:
  378.     [图片]请问这个BSSE矫正是什么意思呢
  379.     有谁知道吗 谢谢

  380. A:
  381.     谈谈BSSE校正与Gaussian对它的处理
  382.     http://sobereva.com/46

  383. ----------------------------------------------------
  384. 2015.06.23 17:39:15
  385. Q:
  386.     [图片]
  387.     2000人不远了
  388.     估计到时间会清理一批僵尸

  389. A:
  390.     肯定的

  391. Q:
  392.     能看到那些人从来没说话?

  393. A:
  394.     当然。
  395.     所有人加入时间加入理由都有详细记录,最后发言日期也能看到,不会错杀的

  396. ----------------------------------------------------
  397. 2015.06.23 17:46:11
  398. Q:
  399.     活跃级别会被杀了不?好忧桑啊

  400. A:
  401.     一年内发过言就不会被踢

  402. ----------------------------------------------------
  403. 2015.06.23 17:47:40
  404. Q:
  405.     理论上能没能预测在一在分子上引入什么样的基团可以降低熔点。

  406. A:
  407.     根据QSPR关系,分子的性质,包括熔点、蒸发焓、密度等等,都是分子描述符的函数。分子描述符包括分子体积、质量、偶极矩、分子表面上静电势极大极小值/平均值/方差等等。
  408.     因此可以对分子进行各种修饰,然后计算描述符,通过QSPR方程预测出相应的熔点
  409.    
  410.     更精准的方式则是做动力学模拟。步骤更复杂,耗时也高。

  411. ----------------------------------------------------
  412. 2015.06.23 18:40:33
  413. Q:
  414.     对于比如C60这样的不含氢的分子,用6-31g(d,p),是不是这个p加上去是没有任何意义的啊,因为p是对氢加的极化函数,是这样吗
  415.     2d是什么意思?

  416. A:
  417.      是的。6-31G*就够了。

  418. ----------------------------------------------------
  419. 2015.06.23 18:57:01
  420. Q:
  421.     在用multiwfn做临界点BCP分析的时候 是不是只能用WFN格式的啊 能够用FCHK格式的了?

  422. A:
  423.      wfn/wfx/fch/molden格式都可以。凡是要用wfn格式的都能用wfx/fch/molden格式。

  424. ----------------------------------------------------
  425. 2015.06.23 18:57:41
  426. Q:
  427.     [图片]大家好。我想请教下,我想在水溶剂下做一个分子的NMR,这是我输入文件
  428.     对吗
  429.     [图片]提示的错误是这样的
  430.     只能 做气相中的NMR吗

  431. A:
  432.      别把外迭代和NMR搅在一起。你写的关键词几乎大多都是多余的。直接写
  433.     PBE1PBE/6-311++G** NMR SCRF=SMD
  434.     就够了

  435. ----------------------------------------------------
  436. 2015.06.23 19:20:39
  437. Q:
  438.     请问一下~VMD的这个键的半径,只能是0.1   0.2   这样的间隔么  不能0.15啥的么
  439.     [图片]

  440. A:
  441.      你选file-log tcl commands to console,然后你调节这个的时候屏幕上就会显示相应的指令,然后重新输入那条指令,但是把参数写成0.15即可。

  442. ----------------------------------------------------
  443. 2015.06.23 19:22:29
  444. Q:
  445.     [图片]
  446.     这个打不开

  447. A:
  448.      不要在本群宣扬工口或奇怪的网站,警告一次。

  449. ----------------------------------------------------
  450. 2015.06.23 19:49:33
  451. Q:
  452.     嗯  谢谢老师   再问一下老师,这个作用指认为C-H------pi可以么

  453. A:
  454.     可以

  455. ----------------------------------------------------
  456. 2015.06.23 19:51:24
  457. Q:
  458.     显示溶剂化与隐式溶剂化模型的各有什么好处

  459. A:
  460.     [图片]

  461. ----------------------------------------------------
  462. 2015.06.23 19:53:46
  463. Q:
  464.     PCM 和cpcm 是一样的?

  465. A:
  466.     不一样,但类似

  467. ----------------------------------------------------
  468. 2015.06.23 20:24:28
  469. Q:
  470.     优化不考虑对称性限制 是加nosymm吗

  471. A:
  472.     恩
  473.     但是如果初始结构严格满足对称性,即便加了这个往往对称性也不会被破坏

  474. Q:
  475.     老师 我不考虑对称性 为什么看不到电子态的不可约表示了
  476.     哦哦

  477. A:
  478.     此时相对于用C1点群,显然没有那些不可约表示

  479. Q:
  480.     除了c1点群 都可以显示不可约表示吗

  481. A:
  482.     C1不可约表示全是A,没什么意义
  483.     其它点群有哪些不可约表示看特征标表

  484. ----------------------------------------------------
  485. 2015.06.23 20:27:07
  486. Q:
  487.     gaussian有把结构往高对称去优化的倾向么

  488. A:
  489.     没有此倾向,纯粹取决于势能面

  490. ----------------------------------------------------
  491. 2015.06.23 21:35:39
  492. Q:
  493.      冒昧请教一个问题:
  494.     扫描分子的某个键长 资料得到的这个图的横坐标是coordinate[图片]
  495.     可是我得到的是(# B3LYP/6-311++G(d,p) opt(modredundant,maxcycles=500) iop(5/13=1) Test ):
  496.     [图片]
  497.     横坐标怎么就不一样了[Emoticon]?
  498.     (想知道有什么方法可以得到资料那样的横坐标)

  499. A:
  500.     真不希望视野里再出现maxcyc=几百,以及iop(5/13=1)这样的关键词

  501. ----------------------------------------------------
  502. 2015.06.23 21:43:59
  503. Q:
  504.     [图片]文献上说这个反应能自发进行,我觉得不太可能。

  505. A:
  506.     原子数都不一样,自然不可能反应

  507. ----------------------------------------------------
  508. 2015.06.23 22:09:55
  509. Q:
  510.      帮解答一下小草儿的问题~

  511. A:
  512.     hrtang已经解答了

  513. ----------------------------------------------------
  514. 2015.06.23 22:50:39
  515. Q:
  516.      咨询一下我用

  517. A:
  518.      操作过程不对

  519. ----------------------------------------------------
  520. 2015.06.23 22:52:57
  521. Q:
  522.      咨询一下我用multiwfn计算了一下ICSS中苯这个例子,out文件算完了,再用multiwfn计算,multiwfn软件就自动关闭了,这是什么原因啊,求指点!
  523.      用高斯算完了,怎么再用Multiwfn算到那个Isotropic选择完1之后Multiwfn软件就自懂关闭了
  524.      您知不知道那块操作出问题了呢

  525. A:
  526.      看此文开头提供的文件
  527.     http://sobereva.com/216

  528. ----------------------------------------------------
  529. 2015.06.23 22:59:39
  530. Q:
  531.     请教下分子的电子态如何判断啊?比如说Σ是如何判断的。有相关的材料可以阅读吗?或者帖子啥的?~拜谢[Emoticon]

  532. A:
  533.      看Levine的quantum chemistry的双原子分子部分

  534. ----------------------------------------------------
  535. 2015.06.23 23:09:42
  536. Q:
  537.     请教大家,有没有什么详细的理论资料介绍在DFT计算中加外电场?

  538. A:
  539.     [图片]

  540. Q:
  541.     谢谢sob,这是操作方面的,有没有介绍理论实现细节的材料?

  542. A:
  543.     没有什么值得说的,就是哈密顿算符里加一项外场和体系电荷的作用势

  544. Q:
  545.     谢谢sob,那比如说,如果我想在ab initio MD中加入一个交流变化的外电场,那是不是这个外电场就会同时产生磁场,这样外磁场也要考虑?目前有软件能实现这样的功能吗?

  546. A:
  547.     octopus

  548. ----------------------------------------------------
  549. 2015.06.23 23:15:22
  550. Q:
  551.     老师:您好!之前看您说
  552.     如果初始结构严格满足对称性,即便加了nosymm这个往往对称性也不会被破坏。请问nosymm和symm=veryloose、symm=loose,这几个有什么区别吗?使用的时候需要注意写什么吗?
  553.     使用时需要注意些什么?

  554. A:
  555.     这控制的是对点群判断的严格程度。平时没必要写veryloose、loose

  556. Q:
  557.     那在什么情况下可以写veryloose,loose呢?

  558. A:
  559.     程序没判断对实际点群的时候
  560.     尤其是体系对称性很高,如C60,默认设置下程序不容易判断对,需要把判断阈值设松一点

  561. Q:
  562.     在把分子对称性调高的时候加入,是这样理解的吗?

  563. A:
  564.     不是调高,而是没判断对当前的点群,勿把点群判断低了
  565.     [图片]

  566. ----------------------------------------------------
  567. 2015.06.23 23:34:46
  568. Q:
  569.      我之前的计算中断了,restart之后能不能把之前的rwf文件删了呢?

  570. A:
  571.     可以

  572. Q:
  573.     噢,那restart是不是只需要chk文件,其他的都可以删去?

  574. A:
  575.     恩

  576. ----------------------------------------------------
  577. 2015.06.23 23:46:31
  578. Q:
  579.      [图片]我下载了您算的结果可是还是不行,这是怎么回事啊

  580. A:
  581.     仔细阅读屏幕提示
  582.     文件名又不是benzene0001、benzene0002...

  583. ----------------------------------------------------
  584. 2015.06.24 00:17:51
  585. Q:
  586.     [图片]
  587.     好低调啊
  588.     黎乐民是南大的啊
  589.     黎叔华
  590.     黎叔很生气
  591.     人家是黎书华

  592. A:
  593.     很久很久很久以前,把叔本华和黎书华名字还搞混了

  594. Q:
  595.     刘文剑 vs 刘剑
  596.     Schleyer  vs Schlegel

  597. ----------------------------------------------------
  598. 2015.06.24 01:08:24
  599. Q:
  600.     我是新手求助高手:文献上是这么说的molecular orbital calculations were performed by using the hybrid B3LYP functional in combination with the triple zeta basis set。其中 triple zeta basis set我在计算的时候直接写TZ还是写LANL2TZ?谢谢

  601. A:
  602.      文献没明确写清楚,作者自己才知道

  603. ----------------------------------------------------
  604. 2015.06.24 01:09:50
  605. Q:
  606.     Sob老师:您的博文<Multiwfn支持的弱相互作用的分析方法概览> 中提到的:7.3 对单体和复合物状态下的原子电荷求差值得到变化量,可以再加和成片段或分子的变化量。这样一来,弱相互作用发生前后每个原子、片段或分子的电子增减了多少就非常清楚了。这个方法极为重要、常用!
  607.     这里的电荷是指NBO分析得到的NPA电荷吗[图片]

  608. A:
  609.      原子电荷定义方法千差万别,看群共享里的《原子电荷计算方法的对比》。

  610. ----------------------------------------------------
  611. 2015.06.24 04:09:10
  612. Q:
  613.     各位大神,泛函M06-X 与M06-2X有什么区别?哪一个算弱作用力更准确呢?!谢谢
  614.     [Emoticon][Emoticon][Emoticon]谢谢谢谢 @华电ding

  615. A:
  616.      根本没有M06-X

  617. ----------------------------------------------------
  618. 2015.06.24 04:19:45
  619. Q:
  620.     [Emoticon][Emoticon]谢谢谢谢,这两个泛函有什么区别呢?!能麻烦您帮忙解答一下吗?
  621.     各位大神,算分子间弱作用力,除了m06-2x外还有什么好的泛函可以使用吗?谢谢各位大神

  622. A:
  623.    
  624.     乱谈DFT-D
  625.     http://sobereva.com/83
  626.     DFT-D色散校正的使用
  627.     http://sobereva.com/210
  628.     大体系弱相互作用计算的解决之道
  629.     http://sobereva.com/214
  630.    
  631.     弱相互作用的计算在今年8月的科音量子化学暑期基础班里也会详细讲

  632. ----------------------------------------------------
  633. 2015.06.24 05:28:03
  634. Q:
  635.     哪位知道,那种免费的动力学模拟软件,比较好上手,windows版本,谢谢
  636.     好的,谢谢了。我看看这个软件
  637.     我就是跑20个原子左右的团簇体系,可以吗

  638. A:
  639.      gromacs也有windows版,只不过不是官方发布

  640. ----------------------------------------------------
  641. 2015.06.24 06:12:54
  642. Q:
  643.     请教一下,如果我们用高斯计算的时候,计算的结果error,原因是收敛不了的。有没有更好的方法让其收敛。(已经增加循环次数500了,还是收敛不了o(╯□╰)o)
  644.     换句话说就是输入文件怎么编命令?
  645.     老师的文章?请问您指的是……

  646. A:
  647.    
  648.     解决SCF不收敛问题的方法
  649.     http://sobereva.com/61
  650.     量子化学计算中帮助几何优化收敛的常用方法
  651.     http://sobereva.com/164
  652.    
  653.     增加循环次数到500只会白浪费时间。
  654.     无论是科音的量化培训班里、群里、论坛里、博文里,我都反复强调绝对不要用maxcyc=几百这种做法

  655. ----------------------------------------------------
  656. 2015.06.24 15:52:03
  657. Q:
  658.     sob老师:如果没机会去现场听您的课,可以做成视频,我们买课,类似于新东方的在线教学,打破时间空间的限制

  659. A:
  660.     暂无此计划

  661. Q:
  662.     可以考虑先录制,如果考虑好马上就可以推出

  663. A:
  664.     暂不打算这么做

  665. Q:
  666.     [Emoticon]
  667.     有录像就有流出啊

  668. A:
  669.     正解

  670. ----------------------------------------------------
  671. 2015.06.24 16:02:19
  672. Q:
  673.     请教一个问题:怎么用li+在C60里扫描稳定的位置

  674. A:
  675.     直接优化

  676. Q:
  677.     取不同位置优化,得到的是同一结构,能量也差不多。但我想把这个过程体现在文章里,所以想画个类似势能面曲线的图

  678. A:
  679.     那你就扫把

  680. Q:
  681.     取不同位置直接算单点能是不是也可以

  682. A:
  683.     可以

  684. ----------------------------------------------------
  685. 2015.06.24 16:38:06
  686. Q:
  687.     请问在计算一个共轭体系时使用lc-b3lyp/6-311++G* 时 301错误,怎么解决啊?谢谢

  688. A:
  689.     LC不能和杂化泛函一起用,只能用到GGA泛函上

  690. ----------------------------------------------------
  691. 2015.06.24 16:48:52
  692. Q:
  693.     [图片]

  694. A:
  695.     基础量化班培训内容如下:
  696.    
  697.     •        计算化学概览
  698.     •        量子力学概要(介绍所有量子化学学习需要知道的最基本的量子力学知识)
  699.     •        量子化学理论基础(Hartree-Fock、半经验、密度泛函理论、组态相互作用、微扰、耦合簇、电子相关的本质、电子积分、SCF迭代、分子轨道、自旋多重度、开壳层与闭壳层、自旋污染、导数计算等)
  700.     •        基组与赝势(介绍基函数和赝势的基本概念,盘点所有主流的基组和赝势,介绍如何根据具体情况选择最合适的基组)
  701.     •        分子结构建模(各种坐标系的介绍,演示构建各类体系的最佳、最快的途径)
  702.     •        势能面的概念
  703.     •        单点能计算(能量的读取、gap、轨道图形的观看、解决不收敛的方法、对称破缺处理、波函数稳定性测试等)
  704.     •        几何优化(几何优化的原理算法、限制性优化、加速收敛的处理等)
  705.     •        过渡态搜索与反应速率的计算(过渡态搜索的基本算法、操作技巧;反应速率计算的原理与操作)
  706.     •        IRC计算
  707.     •        势能面扫描(刚性扫描、柔性扫描、不同坐标下的扫描)
  708.     •        振动分析(谐振频率计算的原理、非谐振效应、频率校正因子等)
  709.     •        热力学性质与能量相关性质的计算(ZPE、焓、自由能、反应热、电离能、电子亲和能、电负性、键能等常见属性的计算,并介绍热力学组合算法、构象分布、平衡常数的计算)
  710.     •        红外、拉曼、VCD、ROA光谱的计算与绘制
  711.     •        NMR、核自旋耦合常数、磁化率的计算
  712.     •        极化率与超极化率的计算
  713.     •        弱相互作用的计算(弱相互作用的本质、计算方法的选择、BSSE的考虑、色散校正、弱相互作用的分析等)
  714.     •        溶剂效应(溶剂的影响、非平衡溶剂效应、溶剂化能的计算、pKa的计算等)
  715.     •        激发态与紫外、荧光、ECD光谱的计算和绘制(同时包括各类光化学过程、势能面交叉、旋轨耦合、激发态寿命、振动耦合等概念的讨论)
  716.     •        盘点主流计算化学软件和理论、计算化学书籍
  717.     •        其它(相对论效应、周期性计算简介、设置外电场、背景电荷、偶/多极矩的计算等)

  718. ----------------------------------------------------
  719. 2015.06.24 16:50:15
  720. Q:
  721.     求两篇文献,谢谢!
  722.     http://journals.aps.org/prl/abstract/10.1103/PhysRevLett.104.136102
  723.     http://journals.aps.org/prl/abstract/10.1103/PhysRevLett.87.216102

  724. A:
  725.     自行用sci-hub.org解决

  726. ----------------------------------------------------
  727. 2015.06.24 18:04:54
  728. Q:
  729.     请教 mmff94的力场参数文件在哪里能找到 ?

  730. A:
  731.     tinker支持MMFF94,可以看看

  732. ----------------------------------------------------
  733. 2015.06.24 18:13:06
  734. Q:
  735.     [图片]    大神么,能帮我看看我安装出现这个问题是怎么回事?

  736. A:
  737.      Multiwfn入门tips
  738.     http://sobereva.com/167

  739. ----------------------------------------------------
  740. 2015.06.24 18:13:38
  741. Q:
  742.     谁有CPMD的压缩文件,帮忙传一个,官方申请的都两天了,还没动静。谢谢啊!

  743. A:
  744.     计算化学公社-量子化学版-资源分享区里面就有CPMD

  745. ----------------------------------------------------
  746. 2015.06.24 18:21:01
  747. Q:
  748.     formchk显示内存过小,这个怎么处理

  749. A:
  750.     当你转换过大的chk文件的时候,需要在环境变量中设置:export GAUSS_MEMDEF=2000000000,即先尝试设置2GB

  751. ----------------------------------------------------
  752. 2015.06.24 18:33:11
  753. Q:
  754.     RDG分析中CH...HC之间的spike值小于-0.01  是否可解释为二氢键?

  755. A:
  756.     这类作用算不上什么二氢键
  757.     这类正负电荷吸引才算
  758.     [图片]

  759. Q:
  760.     那属于范德华?  原文不是说小于0.005的属于范德华吗?

  761. A:
  762.     如果是因为体系结构原因导致挨得太近才导致它们之间密度较高,这什么都不算。弱相互作用分析那个博文里也说了,相互作用的考察应当是那个相互作用处在平衡结构下。

  763. Q:
  764.     结构经过优化,是否可以说是处在平衡位置?  只是接触面积大、等值面比较多,这是其中的一个等值面,这样的情况是否属于“因为体系结构原因导致挨得太近”?

  765. A:
  766.     这就是我指的因为结构原因被迫挨得太近
  767.     [图片]
  768.     并非因为自发的成键而离得近,是结构迫使其挨得太近
  769.     一些文章瞎讨论一堆,说这C-H----H-C有某某某吸引/成键作用

  770. ----------------------------------------------------
  771. 2015.06.24 18:37:24
  772. Q:
  773.     sob老师 ,multiwfn搜索(3,-1)的BCPs时,能不能把共价键的CPs除去,只留下剩下的弱相互作用的CPs。

  774. A:
  775.     如果是因为体系结构原因导致挨得太近才导致它们之间密度较高,这什么都不算。弱相互作用分析那个博文里也说了,相互作用的考察应当是那个相互作用处在平衡结构下。
  776.     选
  777.     -4 Modify or export CPs (critical points)
  778.     再选
  779.     2 Delete some CPs

  780. Q:
  781.     您觉得可行不?老师

  782. A:
  783.     没有相应的选项控制这个

  784. Q:
  785.     恩,我也没有找到,那我用 2 Delete some CPs吧

  786. A:
  787.     特别需要这个功能也可以给你扩展一下程序

  788. Q:
  789.     哇塞,老师,您太棒了,我真需要这个功能,
  790.     我昨天就用 2 Delete some CPs,昨天一天才处理了三个结构,而且有漏的CP点

  791. A:
  792.      正在弄

  793. ----------------------------------------------------
  794. 2015.06.24 18:42:05
  795. Q:
  796.     因为我的体系比较大,成共价键的CPs比较多,找出来就很麻烦。我的想法是:能不能再搜索的过程中就不要搜索电子密度大于0.2的
  797.     sob老师,能不能帮忙扩展下程序,在搜索(3,-1)type CPs 的过程中可以选择不要搜索电子密度大于0.2的,[Emoticon]

  798. A:
  799.     已经在Multiwfn主页上更新了。你可以先用选项2搜索出NCP,然后选-1,再选9,输入0,0.2,然后再返回去用选项3搜索BCP,此时在CP搜索过程中,只有密度在0~0.2范围内的会被保留。
  800.    
  801.     wyy(1070840197)  10:51:49
  802.     sob老师,能不能帮忙扩展下程序,在搜索(3,-1)type CPs 的过程中可以选择不要搜索电子密度大于0.2的,[Emoticon]

  803. Q:
  804.     老师这速度 真是佩服
  805.     谢谢sob老师,这速度,佩服,太感谢了

  806. A:
  807.     要不是有其它事耽误,最多半个小时就能搞完

  808. Q:
  809.     我要是这么厉害就好了

  810. A:
  811.     关键是程序写得清楚,代码干净易读,这点功能总共也就加了十来行代码

  812. ----------------------------------------------------
  813. 2015.06.24 18:54:42
  814. Q:
  815.     [图片] magicmonk 分享文件 10:54:42
  816.     "JCR2015_all_.xls" 下载

  817. A:
  818.     群共享已经有这个了

  819. ----------------------------------------------------
  820. 2015.06.24 18:55:47
  821. Q:
  822.     请教一个问题:在小基组下优化的结构,要比较大基组下的稳定性,做了大机组的单点能计算,要加溶剂化效应,用大基组?还是小基组呢?哪个更好?

  823. A:
  824.     溶剂化能应该用拟合溶剂模型时的级别,比如SMD最好的级别是M052X/6-31G*,用更大基组结果反而差

  825. Q:
  826.     我用B3LYP/6-31G*优化的,现在用6-311+g(d,p)做了单点能,不知道用6-31G*、6-311+g(d,p),哪一个做溶剂化能呢?
  827.     希望SOB老师指点
  828.     那我就用6-31g*,就可以了?

  829. A:
  830.     [图片]

  831. ----------------------------------------------------
  832. 2015.06.24 18:56:48
  833. Q:
  834.     sob今年不分析一下计算方面期刊的趋势?

  835. A:
  836.     适合理论、计算化学投稿的期刊及其2014年影响因子(2015年公布)
  837.     http://sobereva.com/296

  838. ----------------------------------------------------
  839. 2015.06.24 19:03:29
  840. Q:
  841.     刚发现量化基础班和Multiwfn程序使用培训班中间隔了两天时间

  842. A:
  843.      场地方的限制,不得已如此。
  844.     不过可以趁这个时间再休息休息、复习复习(内容很多,幻灯片近1000张),没搞明白的话等波函数分析班的时候也可以一起解答。

  845. ----------------------------------------------------
  846. 2015.06.24 19:03:32
  847. Q:
  848.     请教sob老师,如果想把moclus用于Gaussian批量计算的话,能否实现呢?比如还想保存chk文件呢?
  849.     [图片]

  850. A:
  851.      molclus本身就会调用高斯批量计算,不明白你的意思
  852.     #的位置写错了
  853.     chk没法一一保留

  854. Q:
  855.     嗯,就是不知道能否保留批量计算过程中的详细信息,因为我看out文件的话缺少波函数分析的信息。

  856. A:
  857.     不能

  858. ----------------------------------------------------
  859. 2015.06.24 19:16:21
  860. Q:
  861.      Sob老师,我想对比氨分子参与到硫酸二聚体前后,两个硫酸之间所有氢键作用的变化,我这么操作是否合理?
  862.     1.先做一个氨分子不参与的纯硫酸二聚体的RDG图,找到spike位置
  863.     2.将优化好的(两个硫酸+氨分子)团簇中氨分子去掉,得到wfn文件再做RDG
  864.     3.由于两个RDG图的spike均只体现硫酸二聚体的弱相互作用,因此对比两个图中的spike位置。
  865.     [图片]

  866. A:
  867.     这不合理。不光是从几何结构上影响,氨的电子(比如导致电子的转移)也会影响硫酸之间的作用

  868. Q:
  869.      那我想要考虑硫酸二聚体那一部分的弱作用力,可以这么操作:
  870.     设定格点数据计算范围,让计算的盒子中心处在某个弱相互作用区域中央,适当调节盒子的延展范围,让计算的空间区域恰好囊括那个弱相互作用区域

  871. A:
  872.     对

  873. ----------------------------------------------------
  874. 2015.06.24 19:18:28
  875. Q:
  876.     在vmd中如何lable BCP index 而不是lable atom[图片]

  877. A:
  878.     在VMD中显示原子序号的方法
  879.     http://sobereva.com/197
  880.    
  881.     里面的id参数输入显示BCP对应的那个id就能只把BCP的标号显示出来

  882. ----------------------------------------------------
  883. 2015.06.24 20:09:26
  884. Q:
  885.      sob老师好,进行高斯MD时,我需要用时间作为步长,并且选择不同的温度,可是高斯若选择时间步长,只能GradientOnly,并输入某一温度下所有原子的初始速度,您知道如何能得到某一温度下的初始速度吗?Multiwfn是否可以给出某一温度的初始速度吗?谢谢

  886. A:
  887.     Multiwfn没这个功能。
  888.     温度反应的是平均动能,自行写个程序,让速度在一定程度内随机分布,而且平均动能对应你需要的温度即可。
  889.     如果不会写也可以参考分子动力学程序的源码,都可以根据设定的温度随机生成初始速度。

  890. Q:
  891.     好的,谢谢sob老师
  892.     那我有个简单的方法,不知道是否正确。
  893.     我用质量加权步长,做BOMD,输出文件里面会有以下两部分:
  894.    
  895.     Cartesian coordinates: (bohr)
  896.     I=    1 X=  -2.906865256658D+00 Y=   7.735784453974D-02 Z=   2.268174450444D-01
  897.     I=    2 X=  -3.507391366331D+00 Y=  -6.168637802479D-01 Z=  -1.224258313529D+00
  898.     I=    3 X=  -8.910749910962D-01 Y=   3.070240426046D-01 Z=   1.630056378567D-01
  899.     I=    4 X=   2.760162082140D+00 Y=  -7.902684073311D-02 Z=  -2.387267758202D-01
  900.     I=    5 X=   3.091158189998D+00 Y=  -9.303973272618D-01 Z=   1.424247840108D+00
  901.     I=    6 X=   3.635593998476D+00 Y=   1.266725245559D+00 Z=  -1.739854433188D-01
  902.     MW cartesian velocity: (sqrt(amu)*bohr/ sec)
  903.     I=    1 X=   4.523317494829D+13 Y=  -3.689497215427D+13 Z=   3.703257091831D+13
  904.     I=    2 X=   4.728199101524D+13 Y=   1.169636447502D+14 Z=  -2.041033595795D+13
  905.     I=    3 X=  -6.013600794576D+13 Y=   5.472232688398D+12 Z=   2.115175166759D+13
  906.     I=    4 X=  -2.060671197266D+13 Y=   7.883213056708D+12 Z=   1.102790135581D+13
  907.     I=    5 X=  -5.693315437223D+12 Y=   4.500653047233D+13 Z=  -1.404532560750D+14
  908.     I=    6 X=  -7.955976612752D+13 Y=  -5.186513194880D+13 Z=  -5.175185020973D+13
  909.    
  910.     第二部分MW cartesian velocity: (sqrt(amu)*bohr/ sec)
  911.     应该是质量加权的笛卡尔速度,单位也对着呢,
  912.     可是第一部分也像速度,可是单位是bohr,不是bohr/ sec.
  913.     您觉得第一部分会不会是原子单位的笛卡尔速度呢

  914. A:
  915.     换算一下就知道了,大抵是

  916. Q:
  917.     amu的平方根值如何算呢,我想算算对不对

  918. A:
  919.     就是原子的原子质量

  920. Q:
  921.     哦,这样呀,好的,谢谢sob老师了。
  922.     那麻烦再问下,我采用质量加权步长,设定了具体温度,可是我看到MD过程中,每一步的速度都不一样,那是不是都是在我设定温度下的模拟呢。温度定了,速度会在一定程度内随机分布?

  923. A:
  924.     坐标、速度怎么改变是根据牛顿运动方程来演变的

  925. Q:
  926.     那都是在我设定一定温度下变化的吗?
  927.     因为整个MD我都设定的是600K

  928. A:
  929.     模拟过程中总能量是保守的,动能和势能会发生一定交换,衡量温度的是动能

  930. Q:
  931.     那这么说,以上的速度我还不能用,因为如果动能和势能交换了,就不是我要的600K了

  932. A:
  933.     对
  934.    
  935.     保持600K严格来说应通过热浴来实现,通过调节速度来让对应的温度一直是你设的值,但高斯的AIMD没有。

  936. Q:
  937.     哦,这样。
  938.     那我设定温度是600K,模拟的第一个速度会不会就是对应600K的动能下的初始速度呢,用的高斯BOMD。如果是的话,第一个速度我也可以拿来用,后面的因为像您说的动能势能有交换,应该不是600K
  939.     第一个速度,如果没有交换,应该还是600K的动能,那对应的速度应该就是600K的
  940.     我是这样想的,不知道是否正确

  941. A:
  942.     对

  943. ----------------------------------------------------
  944. 2015.06.24 21:12:17
  945. Q:
  946.     算这一步时老是出错,怎么办?
  947.     Performing Berny Optimization
  948.     [图片]
  949.     求大神指教!

  950. A:
  951.     opt=cartesian

  952. ----------------------------------------------------
  953. 2015.06.24 21:34:34
  954. Q:
  955.     群里有老师用ACEMD这个软件的吗,我想知道这个软件是不是开源的,我下载了这个软件,但是不会看是不是开源,谢谢

  956. A:
  957.     显然不开源,否则没法卖了

  958. Q:
  959.     这个软件不是可以自由下载[图片]

  960. A:
  961.     完整版不会免费提供的

  962. ----------------------------------------------------
  963. 2015.06.24 22:17:18
  964. Q:
  965.      [图片]
  966.     激发态优化出现这样的情况,是继续算还是改一些参数呢

  967. A:
  968.     算起码15步再说

  969. Q:
  970.     这4步就1个月了
  971.     能不能降低一点收敛标准
  972.     你看如何降低收敛表针呢
  973.     收敛标准呢

  974. A:
  975.     说明你的计算模型、计算级别设定完全不合适。这样算下去没有未来

  976. Q:
  977.     你看就是gaussian计算荧光的第四步
  978.     降低收敛标准,怎么降低呢?

  979. A:
  980.     根本不是降低收敛标准的事

  981. Q:
  982.     那是,请指教。

  983. A:
  984.     Sobereva(190258442)  14:31:27
  985.     说明你的计算模型、计算级别设定完全不合适。这样算下去没有未来

  986. ----------------------------------------------------
  987. 2015.06.24 22:44:55
  988. Q:
  989.       #p wb97xd/6-31+G(d) TD=(NStates=130,Root=1) SCRF=(Solvent=Water) Geom= Modify Guess=Read Opt=ReadFC
  990.     你看可能减少一些收敛标准
  991.     或者减少一些激发态的数目

  992. A:
  993.     多少原子?

  994. Q:
  995.     89个原子
  996.     NStates=130  需要这么大的nstate吗?

  997. A:
  998.     计算时用了多少核?

  999. Q:
  1000.     Sobereva老师是回答的我的问题吗?只要能量最小化就可以了,是吗?
  1001.     92个原子其中80个C,12个氢

  1002. A:
  1003.     能量最小化也不是必须的

  1004. Q:
  1005.     i7的cpu,用了4个核心。算了接近1个半月了
  1006.     纯CH的体系不用130个态吧?
  1007.     需要。不然算不到200nm

  1008. A:
  1009.     你若是算荧光,写
  1010.     #p B3LYP/6-31G(d) TD SCRF opt

  1011. Q:
  1012.     紫外吸收是哪个态跃迁到200附件的?

  1013. A:
  1014.     是

  1015. Q:
  1016.     我有点糊涂了。上次你不是说算荧光要带+么,
  1017.     如果用白
  1018.     b3lyp得算到300个以上
  1019.     [Emoticon]

  1020. A:
  1021.     算荧光是算荧光,算紫外是算紫外
  1022.     选的泛函耗时高
  1023.     还带弥散
  1024.     优化第一激发态每步却居然算到130个态
  1025.     CPU只有四核
  1026.    
  1027.     算得动才怪

  1028. Q:
  1029.     多问两句,荧光不就是紫外发射么?

  1030. A:
  1031.     [图片]
  1032.     你先弄明白kasha's rule

  1033. ----------------------------------------------------
  1034. 2015.06.24 22:45:51
  1035. Q:
  1036.     刚刚开始学amber,我想用autoduck算一个小分子和蛋白质的对接,蛋白质分子用amber做能量最小化之后还需要动力学模拟吗,是在真空的条件下还是加水分子的条件下模拟呢?有哪位做生物体系的,能帮我解答下吗。

  1037. A:
  1038.     不用

  1039. ----------------------------------------------------
  1040. 2015.06.24 23:20:02
  1041. Q:
  1042.     为什么 分子氧 可以控制热异构化的快慢

  1043. A:
  1044.     字体最好别用太大

  1045. ----------------------------------------------------
  1046. 2015.06.24 23:20:35
  1047. Q:
  1048.     @Sobereva大神,想画有机环分子的结构,有点大,怕画出来也优化不了,有什么网站有现成的晶体结构,在它基础上我做opt,谢谢

  1049. A:
  1050.     没有那么恰好的网站正好能找到你需要的晶体结构

  1051. Q:
  1052.     哦,那画出来的有机分子的结构,优化不了怎么办

  1053. A:
  1054.     优化不动就用MOPAC半经验

  1055. Q:
  1056.     优化opt完后,需要用frequent看一下有没有虚频吗

  1057. A:
  1058.     不是必须的。有没有虚频往往凭直觉也能判断。有时有点虚频也不碍事

  1059. ----------------------------------------------------
  1060. 2015.06.24 23:23:48
  1061. Q:
  1062.     老师:加入int=ultrafine会对能量产生影响,这个影响是好还是不好呢?

  1063. A:
  1064.     横向比较时,都加int=ultrafine,相对能量会更精准一点

  1065. Q:
  1066.     恩,好的,马上改,横向比较?老师指的是?

  1067. A:
  1068.     诸如不同构型能量差
  1069.     反应能

  1070. Q:
  1071.     [图片]根据高斯view中的这个对话框,我调整对称性为C3v,但是优化的结果仍然是Cs。然后加入symm=veryloose.也没解决,老师您觉得还可以从哪个角度调整呢?对于该体系用Cs或者是C3v做NBO分析和WBIs和NICs分析和MO轨道分析有差别吗?

  1072. A:
  1073.     用的什么基组?

  1074. Q:
  1075.     6-311++G**和Sdd

  1076. A:
  1077.     把弥散去掉再优化

  1078. ----------------------------------------------------
  1079. 2015.06.24 23:41:29
  1080. Q:
  1081.     S老师,请问像三氟化硼和苯甲醚络合体系的优化是否需要对两个分子的势能进行扫描分析以确定最低能量?另外,如果我需要考虑硼同位素的话,请问您推荐用何种基组计算?

  1082. A:
  1083.      不用扫描,直接优化。
  1084.     几何优化不会体现同位素效应。

  1085. ----------------------------------------------------
  1086. 2015.06.24 23:53:46
  1087. Q:
  1088.     求教:为什么我做NBO分析时,可以正常得到.47文件,但是继续做NBO计算时,会提示 Structure accepted: No low occupancy Lewis orbitals
  1089.    
  1090.     4 low occupancy (<1.9990e) core orbitals found on Nb 1
  1091.    
  1092.     AO with unknown angular symmetry, l =   5
  1093.    
  1094.     AO with unknown angular symmetry found.
  1095.     从而不能正常得到NBO分析结果。同时还会出来个FILE.48文件。大家有了解的吗?~

  1096. A:
  1097.      用的什么基组?

  1098. Q:
  1099.     很大的。aug-cc-pwCVQZ-PP和aug-cc-pVQZ

  1100. A:
  1101.     把基组改小。
  1102.     NBO支持不了那么高角动量
  1103.     而且NBO分析对基组很不敏感,用那么大毫无意义

  1104. ----------------------------------------------------
  1105. 2015.06.25 00:02:19
  1106. Q:
  1107.     大家好有个问题想向大家请教,我是用的方法是B97-1,基组是cc-pwCVQZ-DK和cc-pVQZ,输入信息怎么写呀?这种基组是高斯03自带的吗?

  1108. A:
  1109.     EMSL上拷cc-pwCVQZ-DK定义。
  1110.     cc-pVQZ是自带的
  1111.     DFT用这类基组太浪费

  1112. ----------------------------------------------------
  1113. 2015.06.25 00:04:27
  1114. Q:
  1115.     谢谢S老师。我用[图片]计算出来的反应热是-21.33 KJ/mol,但实验测定值为-51 KJ/mol 左右,请问这个偏差算正常吗?

  1116. A:
  1117.     偏差太大。
  1118.     一般来说建议你改用M062X/def2TZVP来算
  1119.     而且要考虑对焓的热力学校正

  1120. Q:
  1121.     您说的是针对同位素体系采用这个算比较好吗?

  1122. A:
  1123.     和同位素没直接关系

  1124. Q:
  1125.     大家好,请教一下,用B3LYP方法DZVP(d)(opt+3f)基组优化过渡金属离子,基组信息在哪里可以查到
  1126.      老师,我的过渡态优化震荡,想调整步长上限,需要重新calcfc么?我改步长后,力常数算了20多小时还没搞定…

  1127. A:
  1128.     不要用这种非主流基组,没任何好处

  1129. Q:
  1130.     那我做NiBr+催化丁烷,用什么方法跟基组,能给个建议吗?
  1131.     请问我做的是铁原子裂解2分子水这个体系应该用什么方法和基组呢[Emoticon]

  1132. A:
  1133.     穷:TZVP
  1134.     富:def2TZVP

  1135. Q:
  1136.     老师是说我吗

  1137. A:
  1138.     都是

  1139. Q:
  1140.     老师我用小基组sto_3g优化得到中间体构型,到换有大基组6_31优化单点能却不收敛,是我找的中间体不对吗,我该怎么办啊

  1141. A:
  1142.     穷穷:def2SVP
  1143.     穷穷穷:lanl2DZ+6-31G**
  1144.     STO-3G根本别用
  1145.     6-31G*是下限

  1146. Q:
  1147.     哦,我是用的6-31G**大几组不收敛
  1148.     我们师妹跟我探讨我们到底是穷还是富[Emoticon]
  1149.     老师,6-31G*是下限是什么意思啊

  1150. A:
  1151.     比6-31G*还低,被拒风险就很大了

  1152. Q:
  1153.     请问文献中使用的基组和方法和我使用的不一样可以吗?
  1154.     壕:def2QZVPP
  1155.     神:aug-cc-pV5Z
  1156.     仙:aug-cc-pV6Z
  1157.     外星人:[图片]

  1158. A:
  1159.     高斯里用def2QZVPP是壕,ORCA里用def2QZVPP,对中小体系即便是屌丝也能享受

  1160. Q:
  1161.     ORCA开了挂,壕是不屑于用的
  1162.     aug-cc-pvtz和def2-TZVP是一个水平的吗

  1163. A:
  1164.     比后者大得多

  1165. Q:
  1166.     [图片]属于穷穷穷级别

  1167. A:
  1168.     cc-pVTZ和def2-TZVP一个级别,但对于DFT前者耗时略多一些,后者更为理想

  1169. ----------------------------------------------------
  1170. 2015.06.25 01:00:12
  1171. Q:
  1172.     我用SMD溶剂模型作优化,直接报错了,提示说[图片],#部分是opt=gdiis b3lyp/6-31+g(d,p) scrf=(solvent=methanol,smd) nosymm geom=
  1173.     connectivity。。是不是当前结构下不适合用SMD模型呢

  1174. A:
  1175.     数值巧合

  1176. Q:
  1177.     就是说在数值计算上巧合碰到不合适的地方导致失败是吗?

  1178. A:
  1179.     对

  1180. ----------------------------------------------------
  1181. 2015.06.25 01:27:47
  1182. Q:
  1183.     如果我只想计算总态密度和分态密度,是不是定义片段的时候不用定义1和2[图片],而是从3开始定义啊,还是怎么做?

  1184. A:
  1185.     作PDOS定义哪个无所谓。但只有定义1、2才有OPDOS
  1186.    
  1187.     这类问题just try和RTFM大法都能解决

  1188. ----------------------------------------------------
  1189. 2015.06.25 01:54:43
  1190. Q:
  1191.     帮忙下载一篇文献:Accounts of Chemical Research (2012), 45, (3), 404-413. Publisher: (American Chemical Society, ) CODEN:ACHRE4 ISSN:0001-4842.谢谢大伙

  1192. A:
  1193.     [图片]
  1194.     没看过群共享里的公社规章制度者都先看看

  1195. Q:
  1196.     Sob老师,我的两篇文献,学校没有权限,实在是下载不了,有一个连连接都没有,别人推荐的,我只能求助大家帮帮忙了,第一个:http://onlinelibrary.wiley.com/doi/10.1002/cphc.201402472/abstract 第二个:找不到连接: Acc Chem Res 2012, 45, 404,谢谢各位帮个忙

  1197. A:
  1198.     不可能找不到链接
  1199.     [图片]

  1200. ----------------------------------------------------
  1201. 2015.06.25 03:12:38
  1202. Q:
  1203.     【话唠】ruanyang 谢谢你,我刚吃饭去了,现在回来了,你说sci-hub可以下载那篇weili的,我打开了SCI-Hub,但是不知道如何使用,你能教教我吗?谢谢

  1204. A:
  1205.     看群共享里本群的规章制度,用法在里面写得很清楚

  1206. ----------------------------------------------------
  1207. 2015.06.25 03:16:32
  1208. Q:
  1209.     [图片]@Sobereva大神,中间有金属的多吡啶结构优化之后,金属和N原子的键,没有了,怎么回事?谢谢大神[Emoticon][Emoticon]

  1210. A:
  1211.     [图片]

  1212. ----------------------------------------------------
  1213. 2015.06.25 03:41:48
  1214. Q:
  1215.     [图片]其中碳键是哪种键?

  1216. A:
  1217.      Phys. Chem. Chem. Phys., 2013, 15, 14377--14383

  1218. ----------------------------------------------------
  1219. 2015.06.25 03:53:19
  1220. Q:
  1221.     请问有懂openCL的大神么?我现在有64核的AMD CPU和96核的Nvidia GPU,已经安装N卡驱动和cuda,请问怎么安装和设置OpenCL的环境才能使用OpenCL异构计算呢?

  1222. A:
  1223.      应该直接就能运行OpenCL程序。起码Windows里用OpenCL测试程序直接就能跑。

  1224. ----------------------------------------------------
  1225. 2015.06.25 04:16:24
  1226. Q:
  1227.     加上noysmm,就不能显示分子轨道?可以这样说吗?
  1228.     哦,好的,也就是说,算与光谱有关的性质时不要加nosymm.可以这样理解吗?
  1229.     手头没书,暂时查不来了...

  1230. A:
  1231.       数值巧合而已

  1232. ----------------------------------------------------
  1233. 2015.06.25 04:18:42
  1234. Q:
  1235.     ppt上有
  1236.     [图片]

  1237. A:
  1238.      幻灯片里说得根本不对。网上的计算化学资料净是严重误导。画一个D3h的NH3(在gv里对称化),写上nosymm后优化,得到的还是D3h的结构而得不到锥形。
  1239.    
  1240.     劝大家不要轻易相信网上的计算化学资料,尤其是中文的。

  1241. ----------------------------------------------------
  1242. 2015.06.25 05:06:56
  1243. Q:
  1244.     群主还有一篇DFT-D乱谈,没看见了,你可以尝试找找
  1245.     百度空间太坑,关了,现在找群主的精华贴子都不好找了,有些群主没有移植过来

  1246. A:
  1247.      计算化学相关的博文全都移植过来了,费了好大劲

  1248. ----------------------------------------------------
  1249. 2015.06.25 05:26:32
  1250. Q:
  1251.     oc1sq=(owx(i)-bounc1x)**2+ *
  1252.     (owy(i)-bounc1y)**2 + *
  1253.     (owz(i)-bounc1z)**2
  1254.     ocensq3=(owx(i)-cenx)**2+
  1255.     &                            (owy(i)-ceny)**2+
  1256.     &                            (owz(i)-cenz)**2
  1257.     把+放到上一行去 还出错?
  1258.     ocensq3=(owx(i)-cenx)**2+
  1259.     &                            (owy(i)-ceny)**2+
  1260.     &                            (owz(i)-cenz)**2
  1261.     用* 代替上面的&
  1262.     ocensq3=(owx(i)-cenx)**2+
  1263.     *                            (owy(i)-ceny)**2+
  1264.     *                            (owz(i)-cenz)**2
  1265.     差
  1266.     frak

  1267. A:
  1268.      对f90,末尾截断的位置写上&即可,比如
  1269.     write(*,*) "lovelive nico &
  1270.     maki nozomi"

  1271. ----------------------------------------------------
  1272. 2015.06.25 16:08:14
  1273. Q:
  1274.     请问一下,对于同一个分子使用相同的基组,分别使用b3lyp与mo62x计算得到的吸收光谱为什么会相差那么大,我该相信哪一种方法,是要相信与实验结果相近的吗?谢谢

  1275. A:
  1276.     乱谈激发态的计算方法
  1277.     http://sobereva.com/265

  1278. ----------------------------------------------------
  1279. 2015.06.25 16:56:15
  1280. Q:
  1281.      sob老师,为什么我计算NPA得到的电荷是0.9多
  1282.     ,而不是设定的-1。我是读取opt计算结果的chk文件来计算的npa。(抱歉,这个问题我昨天问了,但没看到回复,怕是自己错过了,所以再问一下)
  1283.     哦,是-0.9多

  1284. A:
  1285.     把问题发到计算化学公社波函数分析版,输出文件压缩一下也作为附件传上去
  1286.    
  1287.     我优先回答计算化学公社上的问题
  1288.     Sob回答问题的优先级顺序见群共享里群规章制度最后的部分

  1289. ----------------------------------------------------
  1290. 2015.06.25 17:02:46
  1291. Q:
  1292.     [图片]从这个图中可以看到的5,6在LUMO部分对总态密度的贡献很小,是怎么回事啊,还是咋分析啊?

  1293. A:
  1294.     实际就是这样,没有什么怎么回事

  1295. Q:
  1296.     5,6是受体部分的片段,按理来说应该对LUMO贡献比较大啊

  1297. A:
  1298.     都不知道你算的什么体系,上来就一张图,没法回答

  1299. Q:
  1300.     从分子轨道空间分布图可以看到LUMO位于5,6部分,

  1301. A:
  1302.     用的什么基组?

  1303. Q:
  1304.     是染料敏化剂,D-
  1305.     6-31g(d,p)

  1306. A:
  1307.     你先做轨道成分分析,看mulliken/SCPA方法得到的LUMO的轨道成份和图中是否定性一致。
  1308.     PDOS的贡献计算就是根据Mulliken/SCPA方法来的,轨道成份结果合理则PDOS必然合理

  1309. ----------------------------------------------------
  1310. 2015.06.25 17:21:01
  1311. Q:
  1312.     [图片]老师,你手头有没有做这种分子轨道图的文献或者能不能给我建议一下怎么做复杂分子分子轨道图谢谢老师

  1313. A:
  1314.      用Multiwfn的CDA模块作轨道相互作用图,然后根据实际情况适当修改调整
  1315.     使用Multiwfn做电荷分解分析(CDA)、绘制轨道相互作用图
  1316.     http://sobereva.com/166

  1317. ----------------------------------------------------
  1318. 2015.06.25 17:46:11
  1319. Q:
  1320.     科音有微信公众号吗?

  1321. A:
  1322.     没有

  1323. Q:
  1324.     运营微信公众号是不是比较麻烦?

  1325. A:
  1326.     不知道。
  1327.     寡人抵制微信

  1328. Q:
  1329.     也不是很麻烦,要看你的运营的微信公众平台的目的
  1330.     貌似现在扣扣都是工作,微信都是闲聊

  1331. A:
  1332.     闲聊可以在公社大杂烩群

  1333. Q:
  1334.     谢谢,呵呵 。思想家公社和计算化学公社有什么区别啊

  1335. A:
  1336.     计算化学公社是科音直属的高水平计算化学论坛。
  1337.     和本群互为补充

  1338. ----------------------------------------------------
  1339. 2015.06.25 17:48:23
  1340. Q:
  1341.     好吧,

  1342. A:
  1343.     不能。得自己手动标

  1344. Q:
  1345.     [图片]

  1346. A:
  1347.     要么你就把ADCH电荷输出贴到高斯输出文件里Mulliken电荷部分,然后在gv里显示Mulliken电荷即可

  1348. ----------------------------------------------------
  1349. 2015.06.25 18:07:27
  1350. Q:
  1351.     老师 请问BSSE矫正直接在做优化的时候加上Counterpoise=2吗
  1352.     我看过啦
  1353.     就是优化后总是出现问题

  1354. A:
  1355.     [图片]

  1356. Q:
  1357.     你好 这个是我的输入文件[图片]
  1358.     最后计算相互作用能跟我之前计算优化时候的能量有什么区别吗 我想知道A B 分子之间的弱相互作用

  1359. A:
  1360.     1 幻灯片里已经说了,优化和振动分析不要用counterpoise,完全划不来
  1361.     2 B971算弱相互作用糟糕透顶

  1362. Q:
  1363.     那就是将优化出来的结构做下一步counterpoise是吧 还有B971我是看文献中说比较好 要不就是MP2

  1364. A:
  1365.     是

  1366. ----------------------------------------------------
  1367. 2015.06.25 18:34:31
  1368. Q:
  1369.     各位老师们,请教大家个问题啊,初学者,不是很懂
  1370.     [图片]
  1371.     我是研究表面活性剂的,我想得到类似于这些图的结果,用的Gaussian 09,可以做出来这些吗?
  1372.     这必须VMD啊
  1373.     他说的可能是怎么去模拟吧

  1374. A:
  1375.     gromacs+VMD

  1376. Q:
  1377.     就是说高斯不能做到这些吗?
  1378.     模拟的话Gromacs, Lammps都行

  1379. A:
  1380.     计算化学常识问题
  1381.     gmx是首选

  1382. ----------------------------------------------------
  1383. 2015.06.25 18:37:02
  1384. Q:
  1385.     老师:请问怎样理解“NBO轨道定域性强(化学键往往能直接对应),而MO离域性强”这句话?

  1386. A:
  1387.      MO往往涉及很大空间范围,NBO只涉及一个或两个,至多3个原子

  1388. ----------------------------------------------------
  1389. 2015.06.25 18:38:57
  1390. Q:
  1391.     [图片]log文件中有mulliken电荷,这里却显示灰色,请问如何解决?

  1392. A:
  1393.     你先正常算完一个体系,比如水分子,做普通单点计算,用gv看Mulliken电荷

  1394. ----------------------------------------------------
  1395. 2015.06.25 18:40:13
  1396. Q:
  1397.      [图片]老师 我这样做IRC做出来还是只有3个点

  1398. A:
  1399.     步长再加大

  1400. ----------------------------------------------------
  1401. 2015.06.25 18:47:25
  1402. Q:
  1403.     [图片]老师,这个mulliken电荷在输出文件可以看到,用高斯view也可以直接看吗?

  1404. A:
  1405.     对

  1406. ----------------------------------------------------
  1407. 2015.06.25 18:49:31
  1408. Q:
  1409.     那老师们,就是我算一个体系,用B3LYP/6-31G** 去优化,正常结束后用B3LYP/6-31+G**去优化,显示的则为Restarting incremental Fock formation.
  1410.     >>>>>>>>>> Convergence criterion not met

  1411. A:
  1412.     多尝试,自然就明白了
  1413.     解决SCF不收敛问题的方法
  1414.     http://sobereva.com/61

  1415. Q:
  1416.     网上说的都是过了,就是没有解决

  1417. A:
  1418.     不要信网上别的地方的说法
  1419.     大多都是错的

  1420. Q:
  1421.     应该是初始结构有问题吧?

  1422. A:
  1423.     都是基础知识问题,掌握了基础理论和算法知识就不会有这些疑问

  1424. Q:
  1425.     但是初始结构也看不出有什么问题啊

  1426. A:
  1427.     你先好好看我的帖子再问

  1428. ----------------------------------------------------
  1429. 2015.06.25 18:49:52
  1430. Q:
  1431.     sob老师 ,vmd中怎么保存当前的所有操作,下次直接load,就不用一个个设置了

  1432. A:
  1433.     file-save visualization state

  1434. ----------------------------------------------------
  1435. 2015.06.25 18:55:56
  1436. Q:
  1437.      Typecho怎么能在首页只显示摘要而不是全文啊

  1438. A:
  1439.     不属于计算化学相关问题,在大杂烩群问吧

  1440. ----------------------------------------------------
复制代码


本版积分规则 Credits rule

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

GMT+8, 2025-8-17 09:16 , Processed in 0.235517 second(s), 23 queries , Gzip On.

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