计算化学公社

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

[群聊Q&A整理] 公社群聊Q&A整理:2016.06.23 06 ~ 2016.06.26 19 Concate

[复制链接 Copy URL]

3097

帖子

29

威望

1万

eV
积分
17098

Level 6 (一方通行)



相关概况可阅读:http://bbs.keinsci.com/thread-2022-1-1.html


暂无全文附件

本次附件(附件内容与下文一致):
2016.06.23 06 ~ 2016.06.26 19 Concate.txt (54.65 KB, 下载次数 Times of downloads: 3)
  1. 2016.06.23 05:11:31
  2. Q:
  3.     opt=ts freq联合计算
  4.     然后出来,link9999[图片]请问这个是ts的问题,还是freq的问题呢?应该在输出文件的哪里去判断呢?

  5. A:
  6.     几何优化没收敛

  7. ----------------------------------------------------
  8. 2016.06.23 05:12:01
  9. Q:
  10.     请教个问题:linux 下高斯运行提示:段错误  (核心已转储) g09 XX.gjf 这该怎么解决呢???
  11.     请教个问题:linux 下高斯运行提示:段错误  (核心已转储) g09 XX.gjf 这该怎么解决呢???求指导
  12.     谢谢 !输出提示g09权限过高Files in the Gaussian directory are world accessible.
  13.     This must be fixed.
  14.     我用chmod 755 g09之后,仍然出现相同的错误提示。。。

  15. A:
  16.     750

  17. Q:
  18.     谢谢 !尝试了770和750依然不行。。。,还有一个让我费解的现象:高斯的输入(XX.gjf)/输出(XX.log)文件全是绿色的。。。

  19. A:
  20.     管它什么颜色
  21.     Gaussian目录下运行chmod 750 -R *
  22.     心里没数就ls -l再看看权限到底改没改
  23.     750绝对行

  24. ----------------------------------------------------
  25. 2016.06.23 05:28:30
  26. Q:
  27.      老师,我在计算激发态的时候遇到一个问题[图片],是激发态优化的那一步,不仅振子强度=0,而且波长已经上万,而实验值只有400。请问这种情况我应该怎样处理呢

  28. A:
  29.      什么关键词?

  30. Q:
  31.     抱歉抱歉, 老师,我的关键词是:#p opt td b3lyp/genecp fchk=all,然后基组是[图片]
  32.     那个只是输出fchk文件的
  33.      谢谢老师了,学生不知道为什么会得到这种结果,请您指点

  34. A:
  35.     用M062X再试试

  36. ----------------------------------------------------
  37. 2016.06.23 06:06:51
  38. Q:
  39.      老师,看文献上好多从头算动力学都是平衡几个或几十ps就拿来说事,而经典动力学大都是要平衡几十几百ns才行,是因为AIMD太耗时还是因为它的能量比较准,ps级就够了,还是因为AIMD大都涉及化学反应,而化学反应大都在很短时间就完成了(如0.1-1ps)。不太懂,请您指教

  40. A:
  41.     太耗时

  42. ----------------------------------------------------
  43. 2016.06.23 07:28:12
  44. Q:
  45.      老师,计算基态与某个激发态的偶极矩用什么泛函比较好

  46. A:
  47.     B3LYP就可以

  48. ----------------------------------------------------
  49. 2016.06.23 16:24:22
  50. Q:
  51.     请教下各位,多个分子在VMD中显示,怎么找到其中的一个想要的分子呢?谢谢指点

  52. A:
  53.     @阳光肉眼一看就知道。如果分子太多,建立新的representation用一些选择语句将某些部分高亮以便于观看

  54. ----------------------------------------------------
  55. 2016.06.23 16:53:32
  56. Q:
  57.      老师,您好,怎么通过垂直激发,看哪些是禁阻,哪些是允许的?
  58.     各位,如何通过垂直激发,看哪些是禁阻?哪些是允许?

  59. A:
  60.      看振子强度

  61. ----------------------------------------------------
  62. 2016.06.23 17:19:40
  63. Q:
  64.     请问,[图片]这个版本是32位点的还是64位的?

  65. A:
  66.      64
  67.          
  68. ----------------------------------------------------
  69. 2016.06.23 17:19:40
  70. Q:
  71.     gview上面不可输入def-TZVP基组?
  72.     [图片]请问,这样输入和上面的关键字等价吗?
  73.     [图片]和[图片]是等价的吗?

  74. A:
  75.      不等价
  76.      def-TZVP要写成TZVP

  77. ----------------------------------------------------
  78. 2016.06.23 17:27:01
  79. Q:
  80.     [图片][图片][图片]  esp图为什么会这样老师

  81. A:
  82.      操作流程不对

  83. ----------------------------------------------------
  84. 2016.06.23 18:06:26
  85. Q:
  86.      sober老师,我用guassview打开的某个fchk文件,想知道怎么来寻找哪个是我算的这个碱基和HIS+之间形成氢键的分子轨道?
  87.     [图片]
  88.     [图片]

  89. A:
  90.      不要试图用分子轨道分析氢键,毫无意义。氢键怎么分析看
  91.     Multiwfn支持的弱相互作用的分析方法概览
  92.     http://sobereva.com/252

  93. ----------------------------------------------------
  94. 2016.06.23 18:13:06
  95. Q:
  96.     [图片]
  97.     联系molclus 报错,请问该怎么解决?@Sobereva

  98. A:
  99.      图裂了

  100. Q:
  101.     单独运行mopac 没问题

  102. A:
  103.     看MOPAC输出文件找原因
  104.     那就怪了

  105. ----------------------------------------------------
  106. 2016.06.23 18:19:02
  107. Q:
  108.     sob,我想请教一下,在用一些片段做分析的时候,要保持和原来结构的原子序号一致,但是经常拆开成两个片段以后,序号顺序就变了,有什么办法可以重新编号吗?只能手动的调整吗?@Sobereva

  109. A:
  110.     @ D调 gview可以。看http://bbs.keinsci.com/forum.php?mod=viewthread&tid=3585

  111. ----------------------------------------------------
  112. 2016.06.23 18:26:55
  113. Q:
  114.     各位,用cis优化出来的激发态和用tddft优化出来的最后结果会有很大差别吗?
  115.     40个原子的CHO构型tddft优化激发态要一个星期正常吗?

  116. A:
  117.      激发能差很多
  118.      正常与否看输出文件判断。双路16核机子正常情况不需要这么久

  119. ----------------------------------------------------
  120. 2016.06.23 19:05:15
  121. Q:
  122.     [图片]
  123.     安装好出现图片错误怎么解决啊

  124. A:
  125.      机子没装gv,或者安装流程不对(先安高斯,再把gv安到相同目录下)

  126. ----------------------------------------------------
  127. 2016.06.23 19:15:17
  128. Q:
  129.     [图片]请问大家出现这种错误怎么改呢?

  130. A:
  131.      CCSD(T)没解析梯度,把要优化的坐标用变量表示。

  132. ----------------------------------------------------
  133. 2016.06.23 23:32:12
  134. Q:
  135.     MS做动力学模拟怎样确定盒子大小啊?
  136.     没有实验数据怎么办
  137.     是不是主要是要一个密度?

  138. A:
  139.     做NPT模拟,通过控压可自发调整盒子

  140. ----------------------------------------------------
  141. 2016.06.23 23:37:16
  142. Q:
  143.     老师,计算解离能的那个公式,配体为未优化的,那配体的能量我直接从配合物上拆下来计算吗?这样合理吗

  144. A:
  145.      算解离能时候片段要单独优化

  146. Q:
  147.     老师,那结合能计算的时候,配体的能量怎么计算呢?可以将优化完成的配合物分别减去拆开的配体和金属离子吗?这样对吗

  148. A:
  149.     整体能量减去金属和配体能量即是结合能
  150.     直接优化配体,算个单点,就是配体能量

  151. Q:
  152.     老师,那算解离能的时候片段要单独优化,和算结合能时配体的优化差别在哪里呢?没理清楚

  153. A:
  154.     结合能没有严格的计算方式的定义
  155.    
  156.     无非就是片段优化还是不优化。
  157.     解离能是一定要优化片段的,而理论上研究键能则一般不优化片段
  158.        
  159. ----------------------------------------------------
  160. 2016.06.23 23:38:04
  161. Q:
  162.      老师,我通过单重激发态,算20个激发态,然后通过垂直激发,看哪些禁阻?哪些允许,这样输入可以吗?# p opt td=(singlets,nstates=20,direct)

  163. A:
  164.      算垂直激发能不要写opt。direct、singlets都多余,是默认的

  165. ----------------------------------------------------
  166. 2016.06.23 23:43:30
  167. Q:
  168.     --------Markdown Markdown--------

  169. A:
  170.     --------Markdown Markdown--------
  171.      赝势定义的格式不对

  172. ----------------------------------------------------
  173. 2016.06.23 23:49:41
  174. Q:
  175.     请问哪位手中还有有关Gaussian计算零点能,热力学能,焓和自由能数据的解释图片,就是给每个数据后面加个解释的那张图片,之前在咱们群中看过多次,但是现在没有搜到[表情]

  176. A:
  177.     [图片]

  178. ----------------------------------------------------
  179. 2016.06.23 23:50:17
  180. Q:
  181.     请问,溶剂的光学介电常数是怎样得来的?

  182. A:
  183.     实验可以测,也可以根据折射率间接估计
  184.     感兴趣可以看看atkin的物理化学,里面有讨论

  185. ----------------------------------------------------
  186. 2016.06.24 00:14:34
  187. Q:
  188.      老师,这个GV通过fchk做ESP流程我按照老师你的博文做了,好多次,也在win下跑了chk用GV做,都不成功,到底是什么情况[图片]
  189.     为什么在GV里面调取HOMO 和LUMO 轨道是只有两个即 [图片]而用CHEM3D却能查看好多比如HOME+1等  @华电ding
  190.     [图片]

  191. A:
  192.      把计算时用的关键词列出来

  193. ----------------------------------------------------
  194. 2016.06.24 00:20:49
  195. Q:
  196.     请问gauss是否可以计算材料化学上面的紫外可见漫反射光谱(DRS)?

  197. A:
  198.      算不了

  199. ----------------------------------------------------
  200. 2016.06.24 00:34:33
  201. Q:
  202.     各位老师,请问对于多重态,β电子α电子是分开计算的,取轨道某一能量时,是取平均值还是取α轨道就可以了?

  203. A:
  204.      非限制性开壳层时,取平均值没意义。此时“某一轨道”指的只能是某一自旋轨道,你必须考虑自旋

  205. ----------------------------------------------------
  206. 2016.06.24 02:02:43
  207. Q:
  208.     宏观的磁场会不会对微观的芳香化合物环电流产生影响,进而改变反应结果呢?

  209. A:
  210.     有可能

  211. Q:
  212.     可以用量化来计算磁场对环电流的影响么?

  213. A:
  214.     Utilizing AICD and GIMIC programs to study magnetically induced current density (幻灯片)
  215.     http://sobereva.com/148
  216.     使用AICD程序研究电子离域性和磁感应电流密度
  217.     http://sobereva.com/147

  218. ----------------------------------------------------
  219. 2016.06.24 02:02:43
  220. Q:
  221.     请问,我想把名称为DRG resid为1的分子显示为drawing method=VDW,但是不知道怎么写tcl
  222.     我写了第一行,drawing method在display或者draw关键词里面找不到
  223.     atomselect top "name DRG"

  224. A:
  225.      file-log tcl command to console,然后每一步操作对应的tcl语句都在控制台输出。自行在图形界面里切换显示方式,就知道对应的tcl语句了

  226. ----------------------------------------------------
  227. 2016.06.24 04:18:14
  228. Q:
  229.     [图片]老师,请问怎样利用Multiwfn得到图表中的每个原子的FED的数值

  230. A:
  231.      做轨道成分分析
  232.     谈谈轨道成份的计算方法
  233.     http://sobereva.com/131
  234.     PS:我讨厌FED这个词,完全没有意义地去创造一个新词

  235. ----------------------------------------------------
  236. 2016.06.24 04:39:37
  237. Q:
  238.     [图片](甲醛分子)
  239.     请问这个是啥轨道?

  240. A:
  241.      sigma

  242. ----------------------------------------------------
  243. 2016.06.24 05:00:04
  244. Q:
  245.     [图片]
  246.     这个配置3800亏不亏

  247. A:
  248.      peace 性价比不高。内存太小

  249. ----------------------------------------------------
  250. 2016.06.24 05:09:21
  251. Q:
  252.     sob老师,molclus 中的isostat,是不是有原子数限制啊
  253.     [图片]
  254.     我用的是免费版的,应该上限是30个原子吧

  255. A:
  256.      isostat没限制。molclus免费版上限是30个

  257. Q:
  258.     [图片]
  259.       再问一遍,在练习molclus, 运行mopac 时报错,Mopac 单独可以运行,可能的原因是什么
  260.     我更新到Mopac2016了

  261. A:
  262.     看MOPAC输出文件找原因

  263. Q:
  264.     看文件说是密码问题,但是单独运行Mopac 不存在这个问题,所以就不知道怎么办了

  265. A:
  266.      你单独运行的mopac和你用molclus调用的MOPAC确信是同一个么?
  267.     mopac会读取其所在目录下的license信息,单独能运行没理由molclus调用时无法运行

  268. Q:
  269.       恩,是用的同一路径的mopac, 我觉得跟mopac版本有关,你有没有试过用moleclus 运行mopac2016? 现在我ChemDraw 也没法连到Mopac了, 我已经给mopac 发邮件了,看看他们怎么回吧。

  270. A:
  271.     2016没问题
  272.     你也可以直接从dos/linux下,在不同目录下调用mopac2016看看能否运行,能的话没理由molclus没法调用
  273.     2016和2012在用法上没任何区别
  274.        
  275. ----------------------------------------------------
  276. 2016.06.24 06:29:28
  277. Q:
  278.     大家好,我是华师小硕,计算方面的新手。最近在算玻尔兹曼分布,需要用到吉布斯自由能,但不知道自由能怎么算,求助大神
  279.     看不懂[表情]

  280. A:
  281.     谈谈隐式溶剂模型下溶解自由能和体系自由能的计算
  282.     http://sobereva.com/327
  283.    
  284.     看第三节

  285. Q:
  286.     我看了这个贴
  287.     我算的是分子自由能

  288. A:
  289.     怎么算都已经很清楚得写了

  290. Q:
  291.     看不懂

  292. A:
  293.     没救了
  294.     [图片]
  295.    
  296.     还怎么写得更清楚?

  297. ----------------------------------------------------
  298. 2016.06.24 06:41:41
  299. Q:
  300.     本人也是新手。我有一个小问题: 我用IQmol 或者Avogadro 去画分子的结构,然后生成内坐标, 然后 做几何优化。  我可以用优化后的 分子内坐标,或者z-matrix 坐标, 去构建一个分子结构么,就像我刚才用avogadro 画出来的那样, 同时显示键长,键交。  我知道我们可以通过 z-matrix 坐标读出来 键长 键角。 但是想省事,我试了保存xyz 格式 或者gzmat 格式,然后用软件都读不出来。 有没有别的什么办法?
  301.     谢谢

  302. A:
  303.     你用gview画多好,最方便
  304.     内坐标、笛卡尔坐标都直接定义了分子结构

  305. Q:
  306.     我用的是q-chem。  gview 是gaussian  的软件么? 能混用么?

  307. A:
  308.     用gview建出来的结构可以给任何第三方量化程序用

  309. ----------------------------------------------------
  310. 2016.06.24 06:46:53
  311. Q:
  312.     老师,g09e.01用哪个版本的redhat安装?装了一个5,没法运行。

  313. A:
  314.     6.x、7

  315. ----------------------------------------------------
  316. 2016.06.24 07:00:07
  317. Q:
  318.     老师,气相下不是标况下自由能怎么计算,分子不大,想用G4:温度20K,气压10^-3 pa,

  319. A:
  320.     用temperature、pressure设定温度和压力即可

  321. Q:
  322.     需要设定同位素值吗?

  323. A:
  324.     不

  325. ----------------------------------------------------
  326. 2016.06.24 07:05:12
  327. Q:
  328.     请问各位老师,能垒为负(-1.19kcal/mol)会是哪一步出了问题?各关键词如下:
  329.     TS:opt=(calcall,ts,noeigen,tight) m062x/6-311+g** geom=connectivity nosymm
  330.     频率分析:freq m062x/6-311+g** scale=0.970 geom=connectivity
  331.     单点能:b2plyp/ma-tzvp
  332.     IRC:irc=(LQA,maxpoints=50,calcfc,stepsize=5) m062x/6-311g(d,p) geom=connectivity
  333.     复合物OPT:opt=tight m062x/6-311+g**

  334. A:
  335.      你是指电子能量为负?
  336.     走IRC也得6-311+g**,否则和找TS级别不一致

  337. Q:
  338.     是电子能量+zpe为负,老师

  339. A:
  340.     电子能量本身呢?

  341. Q:
  342.     电子能量本身也是负的,-0.002091kcal。。。
  343.     好的老师,我重新算一下IRC。 是IRC基组不同导致复合物找不准确吗?[表情]
  344.      我学习机理不太久,所以对很多数据的意义还把握的不太好。。[表情]
  345.    
  346.     这个barrier应该是没问题的吧?是羟基自由基攻击呋喃上的一个C。

  347. A:
  348.      基于热力学形式的TST算k的时候才必须自由能求差得到自由能垒
  349.    
  350.     m062x/6-311+g**级别下反应物和TS谁单点能高?
  351.     如果还是反应物高,100%算错了。因为IRC两边的点必然能量比TS低,拿一侧的点优化,能量继续降低,注定不可能反倒比TS还高了。
  352.     算得动的话,可以反应物、TS优化都用M062X/ma-TZVP,再看看B2PLYP/ma-TZVP算出来势垒还是不是负的
  353.     B2PLYP/ma-TZVP和M062X/6-311+G**有2kcal/mol的差异很正常

  354. Q:
  355.     刚刚才找到[表情],正如 所说ts高1.2kcal,谢谢。好的老师,我换这个基组试试。那IRC是不是也得更换成ma-tzvp?@Sobereva
  356.     这个反应能垒太低了……

  357. A:
  358.     是,IRC、振动分析用的级别必须与几何优化严格一致

  359. Q:
  360.      为什么这个能垒不准呢?我分别算了电子能量,但是自由能没算,这个情况要用自由能来评价能垒吗?

  361. A:
  362.     [图片]
  363.     [图片]

  364. Q:
  365.     老师,我今天也在看这个ppt,有点疑惑。什么时候适合用电子能量求差,什么时候适合用自由能求差么?@Sobereva
  366.      您的意思是说这种情况下用电子能量算的能垒低估了熵损失,所以要用自由能能垒?那啥时候熵损失会被低估呀?/衰
  367.      您的意思是说这种情况下用电子能量算的能垒低估了熵损失,所以要用自由能能垒?那啥时候熵损失会被低估呀?[表情]

  368. A:
  369.      自由能跟电子能量,最关键是两点差异
  370.     1 ZPE
  371.     2 熵效应
  372.    
  373.     自由能求差获得自由能垒主要用来根据TST算k
  374.     纯理论研究一般用电子能量差,也算是作为自由能差的近似(因为不用做振动分析)

  375. Q:
  376.     谢谢sob老师~之前做了一个路径里一些自由能差,发现有的比电子能量差高,有的比电子能量差底,差别还挺大的(约有10kcal吧),这样的差异电子能量差为什么还能作为自由能差的近似呢?

  377. A:
  378.     电子能量占大头

  379. Q:
  380.     /衰那不会出现有的反应1比反应2电子能垒低,但是反应2比反应1自由能能垒低的情况吗?万一出现了应该信哪一个呢?

  381. A:
  382.     可能会
  383.     没有什么“信哪个”的说法,只不过是给出两种不同的量,都有意义,看你讨论什么

  384. Q:
  385.     噢~原来是这样子。那如果讨论哪一条反应路径更可行用电子能量就行了吗(无论反应的温度压力条件)?

  386. A:
  387.     自由能垒

  388. Q:
  389.     哦哦[表情],那电子能量能垒的意义是什么呢?

  390. A:
  391.     答案在我发的幻灯片里

  392. ----------------------------------------------------
  393. 2016.06.24 07:07:26
  394. Q:
  395.     用m062x优化,再用b2plyp 计算能量?
  396.     这不是扯吗?b2plyp有m06-2x准吗?计算能量

  397. A:
  398.     b2plyp当然比m06-2x准

  399. Q:
  400.       有测试过?

  401. A:
  402.      双杂化泛函怎么可能精度比M062X还低?这样的话谁还用它?
  403.     B2PLYP在双杂化泛函中再烂它也是双杂化泛函,跟普通泛函不是一个档
  404.     瘦死的骆驼比马大
  405.     不用测,属于常识

  406. Q:
  407.     b2plyp用的人多吗?怎么都是见到m06-2x

  408. A:
  409.     多了去了
  410.     更昂贵当然用的人少

  411. ----------------------------------------------------
  412. 2016.06.24 07:19:48
  413. Q:
  414.     sob老师,NMODES是什么程序啊
  415.     [图片]

  416. A:
  417.     不知道。如果做PED可以看
  418.     PED分析软件VEDA4简介
  419.     http://sobereva.com/183
  420.     Potential energy distribution(PED)计算软件GAR2PED使用简介
  421.     http://sobereva.com/75
  422.     以及
  423.     Gaussian中分析振动模式成份的关键词freq=intmodes
  424.     http://sobereva.com/106

  425. Q:
  426.     哦哦,感谢sob老师我就要最后一个

  427. ----------------------------------------------------
  428. 2016.06.24 07:23:02
  429. Q:
  430.     --------Markdown Markdown--------

  431. A:
  432.     --------Markdown Markdown--------
  433.      发射过程是从已经占据了的pi*返回空着的n

  434. ----------------------------------------------------
  435. 2016.06.24 07:27:07
  436. Q:
  437.     Multiwfn不可以做振动模式成分分析么

  438. A:
  439.     不能,和波函数分析八竿子打不着

  440. ----------------------------------------------------
  441. 2016.06.24 07:38:03
  442. Q:
  443.     请问 Multiwfn 软件的 manual 有中文版的吗?

  444. A:
  445.     无
  446.     相关帖子一堆,好几十万字,这里有汇总
  447.     Multiwfn入门tips
  448.     http://sobereva.com/167

  449. ----------------------------------------------------
  450. 2016.06.24 07:44:44
  451. Q:
  452.     还有计算计算熵时,把小于100的频率都认为是100,可以减小误差。
  453.     b2plyp 和 b2plyp(full) 有什么区别?

  454. A:
  455.     后者把内核电子的相关作用也考虑了,前者用冻核近似

  456. Q:
  457.     算了一个水分子,给出的 scf done E 都是一样的

  458. A:
  459.     scf部分当然一样,差异在后面算MP2部分时候才体现

  460. Q:
  461.     计算能垒不是取 scf done E 吗?

  462. A:
  463.     双杂化泛函取最后头输出的能量

  464. Q:
  465.     是 E2(...)=.....    E(B2PLYP)=.....  这一行的 E(B2PLYP) ?

  466. A:
  467.     y

  468. ----------------------------------------------------
  469. 2016.06.24 10:03:18
  470. Q:
  471.     创腾李明宪那个castep 班怎么那么贵
  472.     网上不都有教程吗?
  473.     还有个教学网站

  474. A:
  475.     创腾的培训历来就很贵
  476.     再加上讲课的人来头大,出场费估计不少给

  477. Q:
  478.     网上也有网站的
  479.     搞不懂还要开培训做啥
  480.     [图片]
  481.     6000太多了,黑心
  482.     我知道有个老师去创腾的班总共吃住花了4万
  483.     国家基金都被这么用了

  484. A:
  485.     唉

  486. Q:
  487.     我们想搞点钱,好难[表情]
  488.     算了,反正商业公司的培训再也不去了
  489.     有钱就支持科音吧
  490.     以前被商业公司培训忽悠过几次,啥也没学到

  491. A:
  492.     定价一人6000块钱,来15个人,和定价一人1000块钱,来60个人,商业公司会选前者,我肯定选后者,一方面能造福更多研究者,一方面也对得起自己良心。

  493. ----------------------------------------------------
  494. 2016.06.24 16:33:12
  495. Q:
  496.     [图片][图片]老师这就是昨天做ESP图出错的命令行,ESP的步骤是按照博文做的,好多次都这样,

  497. A:
  498.      fch文件给我

  499. ----------------------------------------------------
  500. 2016.06.24 16:46:41
  501. Q:
  502.     请问研究石墨烯的层间,是否可以在gauss上固定两层之间的的距离?从而模拟实际情况?

  503. A:
  504.      可以

  505. ----------------------------------------------------
  506. 2016.06.24 17:24:20
  507. Q:
  508.     老师,你好,问一下,MS-display style-(Origin defult)格式的区别是什么?是不是模型计算的时候一般采用Origin模式
  509.     刚刚接触计算,有好些地方不明白
  510.     [表情]

  511. A:
  512.      随便选,你看着效果好看就行。不明白的看help

  513. ----------------------------------------------------
  514. 2016.06.24 17:34:47
  515. Q:
  516.      老师 我在做一个COSMO的文件时 采用这样的输入命令 总是算到一半就停下来了 log文件也没有报错 我不知道哪里出问题了 请老师帮我看看 谢谢您了
  517.     [图片]
  518.     [图片]输出文件的最后一页

  519. A:
  520.      就一个原子你优化什么?

  521. ----------------------------------------------------
  522. 2016.06.24 17:40:14
  523. Q:
  524.     请问,M062X算单点,最后输出的[图片]和[图片]都是M062X算出的电子能量呢?

  525. A:
  526.      都是m062x能量

  527. ----------------------------------------------------
  528. 2016.06.24 18:55:46
  529. Q:
  530.     利用orbital info.-show all 为什么不是全部,而是从300-628?
  531.     利用orbital info.-show all 为什么不是全部,而是从300-628? @LiYuanhe

  532. A:
  533.      show all就是全部,如果用wfn/wfx,不包含空轨道,才只会给你占据轨道,这里写了
  534.     使用Multiwfn观看分子轨道
  535.     http://sobereva.com/269

  536. ----------------------------------------------------
  537. 2016.06.24 18:58:42
  538. Q:
  539.     [图片]
  540.     优化至过渡时,碰到这种情况怎么办? 即,能量最高点不是一个峰,而是一个平台。
  541.     [图片]
  542.     优化至过渡态时,碰到这种情况怎么办? 即,能量最高点不是一个峰,而是一个平台。
  543.     [图片]

  544. A:
  545.      如果当前是实际的IRC曲线,说明有三个过渡态,之间有两个很浅的极小点。都要找到

  546. Q:
  547.     irc一般不都是能量依次减小的平滑去西曲线吗

  548. A:
  549.     它这个是复合反应的总IRC
  550.     当然前提是实际就是如此而不是随手画的
  551.        
  552. Q:
  553.     恩,我指的应该是基元反应的一步IRC
  554.     那就说明,他有两个极小在势能面非常平滑的地方,可能不是很好找

  555. A:
  556.     y
  557.        
  558. ----------------------------------------------------
  559. 2016.06.24 19:39:24
  560. Q:
  561.     DFT方法算一个势能面上的鞍点能量,有的加int=ultrafine,有的不加可以么
  562.     Gaussian输出文件的计算时间不对是怎么回事呢

  563. A:
  564.      要加必须都加

  565. Q:
  566.      老师,就是说同一个势能面上不同鞍点可以有的加int=ultrafine,有的不加,对么?

  567. A:
  568.     比较能量当然要加都要加
  569.     用了int=ultrafine会影响势能面

  570. Q:
  571.     [图片] 李老师这个不再优化了,加int=ultrafine跑个单点行不行?@Sobereva

  572. A:
  573.     int=ultrafine对于结构影响很小,所以可以基于没带int=ultrafine的结构带着int=ultrafine只算个单点
  574.      int=ultrafine对于结构影响很小,所以可以基于没带int=ultrafine的结构带着int=ultrafine只算个单点
  575.          
  576. ----------------------------------------------------
  577. 2016.06.24 19:47:51
  578. Q:
  579.     vmd里面如何去掉盒子?我记得昨天有人做过?求指导

  580. A:
  581.      选上
  582.     [图片]

  583. ----------------------------------------------------
  584. 2016.06.24 21:00:48
  585. Q:
  586.     [图片]
  587.     [图片]
  588.     这是什么意思

  589. A:
  590.      即四个核以SMP方式单机并行,即曰:%nproc=4

  591. ----------------------------------------------------
  592. 2016.06.24 21:02:14
  593. A:
  594.      多中心键级、NICS等,哪些能用视环的具体特征而定

  595. Q:
  596.     [图片]我想做这种大环的不同稳定构象芳香性对比,您觉得什么指标比较好?NICS(1)zz,AICD的磁感应电流密度?可以吗?还有什么比较好的参数
  597.     [图片]

  598. A:
  599.      这个不好考察芳香性,而且这么大环,电子整体离域可能性也很低,没什么好考察的
  600.          
  601. ----------------------------------------------------
  602. 2016.06.24 21:29:54
  603. Q:
  604.     请问要求273K条件下的反应焓, 是不是直接将temperture=273 加入关键词就行了?

  605. A:
  606.      y

  607. ----------------------------------------------------
  608. 2016.06.24 21:39:19
  609. Q:
  610.     --------Markdown Markdown--------

  611. A:
  612.     --------Markdown Markdown--------
  613.     也有可能

  614. ----------------------------------------------------
  615. 2016.06.24 22:14:59
  616. A:
  617.     那个lammps培训的公司居然还打电话来,要挟要删掉计算化学公社那个打假帖子
  618.     还死不承认是跑到本QQ群来获取成员qq邮箱地址
  619.     实在太猖獗了

  620. Q:
  621.     特别气人
  622.     死不承认自己的行为

  623. A:
  624.     一定不要去那个培训,一定要向同行相互转告
  625.     一看就是骗子被揭穿,狗急跳墙

  626. Q:
  627.     那个公司的?

  628. A:
  629.     就是老往群成员qq邮箱发培训邮件的那个公司

  630. Q:
  631.     不学lammps。。
  632.     [图片]屏蔽了吗?

  633. A:
  634.     对
  635.     要不然丧心病狂的人指不定干出什么来

  636. Q:
  637.     被谁屏蔽了?

  638. A:
  639.     我屏蔽的

  640. Q:
  641.     被他们威胁了?

  642. A:
  643.     对

  644. Q:
  645.     卧槽[表情]
  646.     是不是堵了他们的财路
  647.     毕竟断人财路如同杀人父母啊

  648. A:
  649.     显然
  650.     猖獗至极

  651. Q:
  652.     无耻
  653.     这么无耻,不过消息还是挺灵通的
  654.     Lammps看看视频,看看命令网页就可以了

  655. A:
  656.     骗子消息当然灵通

  657. Q:
  658.     [图片]
  659.     当骗子还这么嚣张
  660.     把优酷那个视频看完了
  661.     估计就可以了
  662.     没有含金量的培训没用
  663.     可惜我也被骗过
  664.     莫要理会便好,要是再嚣张便举报他们
  665.     举报110
  666.     骗子太多,傻子不够用
  667.     这也有一门学科,从网络安全发展出来的-》 社会工程学
  668.     反正科音光明磊落,扫尽一切邪魔外道
  669.     我想说的是那个公司确实在5月份在苏州办了个lammps培训

  670. A:
  671.     不是苏州那个

  672. Q:
  673.     不是那家吗,

  674. A:
  675.     泰科那个lammps培训是真的(虽然质量怎么样我不知道),假的是之前往大家QQ邮箱里偷偷摸摸发漏洞百出的北京lammps培训通知的那家

  676. Q:
  677.     代理maps的那个公司还是正规的,别的就不清楚了!乱入了啊
  678.     界面也是个鸡肋,到处卖没人要
  679.     [图片]
  680.     说的是这个吧

  681. A:
  682.     就是这个
  683.     联合抵制

  684. Q:
  685.     这个是有问题的
  686.     [图片]搞得还一副挺有情怀的样子。。。
  687.     叫做 北京博宏科睿企业管理顾问有限公司
  688.     企业管理和培训感觉不太搭啊

  689. A:
  690.     稍有点常识的人看了通知都知道是假的
  691.     还说什么他们的学员有固定客户,有固定客户还跑来偷偷摸摸往大家qq邮箱发信息?

  692. Q:
  693.     这样的管理公司怎么也对分子模拟感兴趣
  694.     估计也骗了一些钱了
  695.     所以敢嚣张
  696.     [图片]
  697.     感觉这个公司貌似有点真[表情]

  698. A:
  699.     主页一看就和那个骗子公司不是一个档次的

  700. ----------------------------------------------------
  701. 2016.06.24 22:41:36
  702. Q:
  703.      你的博客里面,**** 自己装FFTW的情况,"注意--enable-float代表编译单精度版本FFTW。如果是给双精度Gromacs用的,应该把--enable-float去掉。"
  704.     我该如何知道我的Gromacs是单精度还是双精度呢?

  705. A:
  706.     默认是单精度

  707. Q:
  708.     [图片]
  709.     那那句,我就直接./configure --prefix=/sob/fftw333 --enable-sse2 --enable-float --enable-shared
  710.     这样就行了?不需要去掉--enable-float?

  711. A:
  712.     编译双精度时应去掉

  713. Q:
  714.     实际上,我默认单精度就够了吧?SOB老师

  715. A:
  716.     对
  717.        
  718. Q:
  719.     [图片]单机并行gromacs的话,这个openmpi不用装,是吧

  720. A:
  721.     不用
  722.        
  723. Q:
  724.     [表情]
  725.      自己装FFTW的情况,make install这一步末尾出现:
  726.     [图片]
  727.     是怎么回事?
  728.     为什么我装个GROMACS这么坎坷啊,晕死

  729. A:
  730.     看更上头

  731. Q:
  732.     [图片]

  733. A:
  734.     [图片]
  735.     编译命令有问题,牵扯到了/sob

  736. Q:
  737.     [图片]
  738.     绝对路径
  739.     那我的如下:./configure --prefix=/sob/fftw333 --enable-sse2 --enable-float --enable-shared
  740.     该怎么写?
  741.     [图片]
  742.     sob老师,我直接把./configure --prefix=/sob/fftw333 --enable-sse2 --enable-float --enable-shared
  743.     这里的sob去掉吗?还是改成liuzeyu,或者是liuzeyu/fftw-3.3.4?我还是没搞清楚,不好意思
  744.     刚才忘了不是root装的,这个有关系?

  745. A:
  746.     肯定得改成自己的路径啊

  747. Q:
  748.     难道这个路径不是自己自动生成吗?我还是不明白,就拿我这个例子,那一部分怎么写?现在的绝对路径是:/home/luzeyu/fftw-3.3.

  749. A:
  750.     /sob/fftw333改成比如/home/luzeyu/fftw33

  751. Q:
  752.     谢谢,我试试哈
  753.     要用管理员装吧

  754. A:
  755.     你装在home你自己目录下面当然不需要管理员

  756. Q:
  757.     那比如我装在跟liuzeyu平行的目录下,就需要管理员吗?

  758. A:
  759.     什么叫平行目录?如果装在liuzeyu目录下面由你自己权限就够了,装在home下面其它用户目录下必须管理员或者那个人的权限

  760. Q:
  761.     那接下来,我要装cmake和gromacs,如果都在我自己的目录下,就都不用管理员权限,是吧

  762. A:
  763.     y

  764. ----------------------------------------------------
  765. 2016.06.24 23:41:31
  766. A:
  767.     pdb、xyz等VMD都能打开

  768. ----------------------------------------------------
  769. 2016.06.24 23:41:31
  770. Q:
  771.     VM中的linux是直接调用内存还是通过vm向windows调用呢?
  772.     window显示内存用了9G多,但是vm只用了30M而已,win7本身也用不了多少内存

  773. A:
  774.     通过vm

  775. ----------------------------------------------------
  776. 2016.06.25 00:03:44
  777. Q:
  778.      老师问一下,在Gaussian 09计算,中氧气是三重态稳定,那氧气吸附到一个三重态的团簇结构上时,算这个吸附结构的自旋多重度是几呢?

  779. A:
  780.     3
  781.     1的情况也可以一起算算

  782. ----------------------------------------------------
  783. 2016.06.25 00:10:17
  784. Q:
  785.     sob老师请教一下,我设置了赝势,但是在输出文件中会有这个提示,是为什么?。。。
  786.     [图片]
  787.     [图片]

  788. A:
  789.      还没报错呢,继续算

  790. ----------------------------------------------------
  791. 2016.06.25 00:38:10
  792. Q:
  793.     谢谢sob耐心!
  794.     [图片]之前出现这个错误,我将输入微调还是一样的,是不是改成内坐标输入会解决?

  795. A:
  796.     opt=cartesian

  797. Q:
  798.     OK
  799.     为什么局部结构微调解决不了?这个问题到底如何引起的

  800. A:
  801.     内坐标本身就有限制,比如角度不能是180度等。而且优化算法内部还牵扯到内坐标与笛卡尔坐标转换之类,也可能转换过程出错。

  802. ----------------------------------------------------
  803. 2016.06.25 00:42:47
  804. Q:
  805.     请问各位老师,VASP算原子电荷有什么比较好的方法吗,用了Bader Charge Analysis感觉不是很好,限制性很差,中性的东西算出来电荷将近-0.2,他们论坛上也有人正面质疑过这个问题,官方人员只是回答调整INCAR里的一个参数调得到你desired的数字。。。这我觉得有点虚。。。有没有像是比较靠谱的Hirhfeld之类的电荷在第一性原理中有什么软件可以实现的吗?

  806. A:
  807.     可以尝试这个
  808.     http://sourceforge.net/projects/ddec/
  809.     可以算DDEC电荷,和Hirsfheld类似,支持VASP

  810. ----------------------------------------------------
  811. 2016.06.25 01:43:56
  812. Q:
  813.     请问本征值波函数模的平方就是电荷密度吗?

  814. A:
  815.      叙述不清
  816.    
  817.     可以说轨道波函数的模方就是这个轨道的电子密度函数

  818. ----------------------------------------------------
  819. 2016.06.25 01:59:11
  820. Q:
  821.     可不可以用m06-2x/6-311+G(d,p)的频率校正因子来替代m06/6-311+G(d,p)的啊?因为没有找到后者的校正因子…… 老师

  822. A:
  823.     不能
  824.     HF成份对校正因子影响极大
  825.     找不到就用其它相似基组的

  826. ----------------------------------------------------
  827. 2016.06.25 03:22:04
  828. Q:
  829.     老师,请教您一个问题。算出来的ELF 盆的电子的多少  能说明哪些问题? 谢谢

  830. A:
  831.      体现有多少电子参与那个电子结构特征区域,比如共价键、孤对电子
  832.     对于键的话,对讨论键的强度、是几重键有用

  833. ----------------------------------------------------
  834. 2016.06.25 04:28:05
  835. Q:
  836.     [图片]各位老师,这个错误到底怎么可以解决调呢,最近一直是这个错误

  837. A:
  838.      没具体信息不好说,看看当前结构有没有诡异之处,或者换换关键词试试

  839. ----------------------------------------------------
  840. 2016.06.25 04:32:31
  841. Q:
  842.     [图片]请问老师用QST2法找过渡态时一开始就出现这个问题是什么意思?

  843. A:
  844.      opt里写cartesian

  845. ----------------------------------------------------
  846. 2016.06.25 04:33:19
  847. Q:
  848.      [图片]缺了一块啊老师

  849. A:
  850.      盒子设得太小。当等值面数值取得太小,等值面很大,就会超过盒子,出现截断
  851.     @天堂鸟

  852. ----------------------------------------------------
  853. 2016.06.25 04:33:19
  854. Q:
  855.     请问大家ccsd方法能否用于寻找过度态呢?

  856. A:
  857.      可以,但一般没必要用这个级别,太耗时

  858. ----------------------------------------------------
  859. 2016.06.25 05:12:05
  860. Q:
  861.     热力学组合方法算能量时,如何获得不同温度的能量,关键字如何写?
  862.      有没有办法一次计算一系列的温度
  863.     常规的方法Sob在课上讲过, 我想用CBS-QB3或者G4MP2计算,遇到了这样的问题
  864.     文献中有用CBS-QB3 计算300-2500k的热力学数据的,20多个温度的热力学数据

  865. A:
  866.      有办法获得其他温度的,手册里有,
  867.    
  868.     [图片]

  869. ----------------------------------------------------
  870. 2016.06.25 06:27:26
  871. Q:
  872.     请问下,为何我这个图显得那么虚化?
  873.     [图片]VMD作图

  874. A:
  875.      用tachyon渲染一下看看结果

  876. ----------------------------------------------------
  877. 2016.06.25 06:47:08
  878. Q:
  879.     用DFT计算第一过渡系金属配合物的单点能,6-311G**和赝势基组lanl2dz哪个更合适些?

  880. A:
  881.      用lanl2TZ或lanl2TZ(f)多好

  882. ----------------------------------------------------
  883. 2016.06.25 07:08:34
  884. Q:
  885.     老师,您的意思是tz比dz更准确?这两种基组计算量差别大吗?

  886. A:
  887.     当然更准,计算量大一些,lanl2DZ算得动TZ也能算

  888. Q:
  889.     /抱拳谢谢老师,又学到了新东西,lanl2tz基组在g09说明书里都找不到
  890.     [表情]谢谢老师,又学到了新东西,lanl2tz基组在g09说明书里都找不到

  891. A:
  892.     EMSL上下载定义

  893. ----------------------------------------------------
  894. 2016.06.25 07:34:58
  895. Q:
  896.      老师,我将计算得到的电子、空穴中心的坐标文件复制到垂直激发的输入中,以Bq,X来代替,但为什么其电子、空穴中心这么偏?[图片]

  897. A:
  898.     高斯会自动把输入坐标调整到标准朝向

  899. Q:
  900.     意思是输出的电子、空穴中心是以 调整后的构型来计算得到的?

  901. A:
  902.     你把输出文件或fch拖到gv,看到的就是标准朝向的,保存成的新的输入文件后,你再把质心写进去。

  903. Q:
  904.     哦,好的

  905. A:
  906.     fch、wfn文件在你没用nosymm都是在标准朝向下的

  907. Q:
  908.     还有个问题,我输入空穴和电子的cub文件,在VMD中显示,原理上应该iso=0.002是空穴的,-0.002是电子的,但我只有显示0.002才能显示电子的图形(-0.002啥都没有),这样电子-空穴交叠cub就只能显示一种颜色,为啥?

  909. A:
  910.     空穴和电子是两个cube文件,都是正值,都得用比如0.002才能同时显示

  911. Q:
  912.      "空穴和电子是两个cube文件,都是正值,都得用比如0.002才能同时显示"
  913.     那holeeleovlp.cub应该能够既显示空穴又显示电子啊,若只有一个值,那不就只有一种颜色吗?

  914. A:
  915.     不可能靠holeeleovlp.cub同时显示空穴和电子,这取的是二者在各个格点的最小值

  916. ----------------------------------------------------
  917. 2016.06.25 11:33:09
  918. Q:
  919.     群里有谁用过wB97XD这个泛函吗?我最近在计算激发态的时候用这个,看文献需要优化其w的大小,不知道怎么弄,看高斯说明书也没有看到对应的关键词

  920. A:
  921.     没事甭去试图优化w

  922. Q:
  923.     但是我外导要我去优化这个
  924.     这两个帖子我都看了

  925. A:
  926.     不知道怎么优化就甭优化,乱优化只不过是去凑数据而已,毫无意义
  927.     他非让你优化你就问他想怎么去优化

  928. Q:
  929.     优化的过程我是知道的,我只想知道这个在高斯输入文件里面怎么体现
  930.     我看官网都找不到

  931. A:
  932.     IOp(3/107) 设定短程和长程HF交换项。0时使用标准HF交换项,MMMMMNNNNN时短程用NNNNN/10000的HF交换而长程用MMMMM/10000的HF交换。这两个数对应于输出文件中HFx后面的wShort和wLong。
  933.     IOp(3/108) 设定短程和长程DFT交换项。0时使用标准DFT交换项,MMMMMNNNNN时短程用NNNNN/10000的DFT交换而长程用MMMMM/10000的DFT交换。这两个数对应于输出文件中DFx后面的wShort和wLong。
  934.     (注:上面两个选项字面上的意思有点误导性。实际上,将3/107和3/108都设为MMMMM00000,就相当于用了ω参数为MMMMM/10000的范围分离泛函。对GGA加LC-时ω默认是0.47/Bohr,等价于0470000000。wB97X默认是ω=0.3,等价于0300000000,但注意wB97X短程是杂化泛函,其杂化比例有参数另设)

  935. Q:
  936.     将3/107和3/108都设为0300000000,就相当于默认的wB97XD了?

  937. A:
  938.     wB97XD和wB97X的参数不一样
  939.     看原文

  940. Q:
  941.     Sob老师,那我wB97XD的omega参数的默认值是多少?我看了iop设置,有两项需要设置:
  942.     [图片]
  943.     这两项都要设置吗?

  944. A:
  945.     iop(3/107=0200000000,3/108=0200000000)

  946. Q:
  947.     Sob老师,还有一个问题,iop(3/107,3/108) 这两项设定是一起变的吧?两个的值在变化的过程中是一致的吧?

  948. A:
  949.     一致
  950.     简谈量子化学计算中DFT泛函的选择
  951.     http://sobereva.com/272
  952.     乱谈激发态的计算方法
  953.     http://sobereva.com/265
  954.    
  955.     这里都提到了

  956. Q:
  957.     109项不用设定啊?

  958. A:
  959.     不设

  960. Q:
  961.     iop(3/107=0200000000,3/108=0200000000) 这两个是wB97XD的默认值吧(0.2)?最后一个问题,如果我要优化其omega,你觉得我变化的间隔多少合适?0.05?或者再小一点?

  962. A:
  963.      是。看你计算能力,可以诸如先0.04,如果间隔太大,不好拟合,到时候再改成0.02,补算没算过的那批点

  964. ----------------------------------------------------
  965. 2016.06.25 11:51:34
  966. Q:
  967.     大家有几个夜猫子啊
  968.    
  969.     不知道有机磷农药分子与乙酰胆碱酯酶相互作用可用gauss算吗

  970. A:
  971.     可能能算

  972. Q:
  973.     哦,就是说理论上可以算,但有可能是计算不动吧

  974. A:
  975.      对。

  976. ----------------------------------------------------
  977. 2016.06.25 13:32:45
  978. Q:
  979.     [图片]测试了一下甲烷,可以算出不同温度下的能量,老师您看一下这个gjf有没有问题

  980. A:
  981.     可以

  982. ----------------------------------------------------
  983. 2016.06.25 15:14:58
  984. Q:
  985.     sob老师,利用multiwfn画DOS图的时候,竖直虚线是HOMO,请问sob老师,这条虚线一定经过DOS虚线图的一个波峰吗?还有,LUMO也过DOS图的一个波峰吗?谢谢。

  986. A:
  987.     不一定经过峰,因为HOMO/LUMO的峰可能淹没在相邻的峰中

  988. Q:
  989.     我能不能通过改变FWHM的值,让其通过波峰?试了几次不行。

  990. A:
  991.     FWHM改小了肯定会通过波峰

  992. Q:
  993.     可是,sob老师,我改变FWHM的值后,发现,出现波峰的位置对应的不是HOMO的能量值。

  994. A:
  995.     上图

  996. Q:
  997.     [图片]
  998.     sob老师,按理说,红色的HOMO-LUMO的值应该小于黑色的值。
  999.     可是,我按照对应的波峰画却大于。
  1000.     返了。

  1001. A:
  1002.     看不明白图

  1003. Q:
  1004.     红色的HOMO-LUMO的值应该大于黑色的值。
  1005.     [图片]

  1006. A:
  1007.     这跟DOS都没直接关系,一看算出来的gap就知道谁大谁小

  1008. Q:
  1009.     这是原图。sob老师,审稿人的意思是改变FWHM的值,让其HOMO和LUMO对应的虚线通过波峰。

  1010. A:
  1011.     改呗

  1012. Q:
  1013.     是的,改变了FWHM的值,

  1014. A:
  1015.     前面的图已经穿过波峰了
  1016.     不管FWHM设大设小,gap就是个固定的值

  1017. Q:
  1018.     是的。可是,你发现没,前面的图中,红色的间距小于了黑色的间距。
  1019.     其实,应该是大于的。2.96>2.04

  1020. A:
  1021.     这大小毫无确切意义
  1022.     实际就是个gap
  1023.     FWHM只是影响图像效果
  1024.     间距终究还是对应的gap

  1025. Q:
  1026.     sob老师,比如,我利用结果计算得到HOMO(-0.2414)的值,我通过改变FWHM的值,希望通过-0.2414这个点的竖直虚线通过一个波峰。

  1027. A:
  1028.     FWHM改小了自然就通过了

  1029. Q:
  1030.     可是,画图后发现,通过改变FWHM的值,通过-0.2414这个点的竖直虚线不是一个波峰。
  1031.     你的意思是还变小。?

  1032. A:
  1033.     FWHM越大展得越宽,越容易被淹没

  1034. Q:
  1035.     我让FWHM=0.02了

  1036. A:
  1037.     FWHM无限小就相当于分立的竖线

  1038. Q:
  1039.     还得变小?
  1040.     我感觉FWHM=0.02已经分裂了很多的峰了。

  1041. A:
  1042.     反复试就知道了

  1043. Q:
  1044.     一般情况小FWHM最小到多少?

  1045. A:
  1046.     无限小

  1047. ----------------------------------------------------
  1048. 2016.06.25 16:09:29
  1049. Q:
  1050.     片段组合波函数初猜的时候出现 ntrerr called from fileio 错误,之前用b3lyp进行初猜没问题,现在改用b3lyp*就出现这个错误了,其他什么都没动

  1051. A:
  1052.     干脆用B3LYP的收敛波函数当成B3LYP*的初猜

  1053. Q:
  1054.     老师,是不是初猜的时候用b3lyp,然后计算的时候用b3lpy*,guess=read读取初猜得到的chk文件?

  1055. A:
  1056.     y

  1057. Q:
  1058.     [表情]谢谢老师
  1059.     这个是不是您提到的Gaussian A02修订版的bug问题

  1060. A:
  1061.     试试E01就知道了

  1062. ----------------------------------------------------
  1063. 2016.06.25 16:59:39
  1064. Q:
  1065.     sob老师,BLYP是的交换想关部分分别是什么?

  1066. A:
  1067.     Becke88, LYP

  1068. ----------------------------------------------------
  1069. 2016.06.25 17:02:38
  1070. Q:
  1071.      老师,gaussian计算的frequency是什么温度下的?可以自己设置计算某一个温度下的frequency吗?

  1072. A:
  1073.     频率和温度无关
  1074.     y

  1075. ----------------------------------------------------
  1076. 2016.06.25 17:03:01
  1077. Q:
  1078.      [图片]
  1079.     老师,我想算石墨烯表面氧官能团对NO2氧化成NO3-的反应,体系电荷设置成0合适还是设置成-1合适呢?设置成0,优化后C=O与NO2不作用(左图),设置成-1,生成了NO3(右图)

  1080. A:
  1081.     带氧的石墨烯是中性,NO2也是中性,整体应该是中性。

  1082. Q:
  1083.     哦,那如果能生成NO3-的话,也是因为基底电子向NO3基团移动,导致NO3带负电(所谓硝酸根),基底这时带正电,但体系整体还是电中性的,对吗?
  1084.     谢谢老师!还有个昨天的问题想请教您一下:
  1085.     [图片]
  1086.     我设计算一个O2分子在缺陷石墨烯表面的裂解过程,不知是否合理。第一步:对三重态体系进行优化,结果是O2与缺陷处发生化学吸附;第二步:将该结构按单重态算单点,发现能量升高18.7kcal/mol,可以认为吸附态氧气分子在常温下可以到单重激发态吗?第三步:对该结构进行几何优化,发现O-O间距从1.20A增大到了1.45A,这时体系能量降到-35.34,可以认为氧气分子发生了裂解,这个过程是无能垒的吗?

  1087. A:
  1088.      也可能氧分子打过来之前就处于单重态。如果放石墨烯上优化O-O键自动伸长了,这个过程就算没势垒。

  1089. Q:
  1090.     [图片]
  1091.     那这个过程算室温下能进行吗,我算单重态和三重态氧气的能量差大约39kcal/mol,但化学吸附后的氧气单重态和三重态差了18.7kcal/mol,大大降低,可以认为这个吸附态的自旋翻转的过程更加容易,甚至室温下就能激发到单重态了吗?
  1092.     [图片]
  1093.     前半部分的过程@Sobereva

  1094. A:
  1095.     氧气单-三重态能量差为22.64 kcal/mol,你算得不准
  1096.     不太大

  1097. Q:
  1098.     [表情]恩,我再算下,但是18.7也比22.6降低明显,是否能说明化学吸附使三重态向单重态的激发更容易在常温环境下更容易进行?

  1099. A:
  1100.     还指不定因计算精度关系,未必真的降低明显

  1101. ----------------------------------------------------
  1102. 2016.06.25 17:16:58
  1103. Q:
  1104.     老师,G4MP2计算到CCSDT的时候出错,是不是也需要加opt=tight int=ultrafine关键词帮助收敛?

  1105. A:
  1106.     ultrafine没用,那是对DFT用的
  1107.     高斯不支持CCSDT
  1108.     看具体什么报错,跟优化无关的用tight毫无意义

  1109. Q:
  1110.     [图片]

  1111. A:
  1112.     哪个g09版本?

  1113. Q:
  1114.     g09w e.01

  1115. A:
  1116.     bug

  1117. ----------------------------------------------------
  1118. 2016.06.25 17:17:33
  1119. Q:
  1120.     sob老师,请问我高斯结构中的HOMO=-0.24214,可是,我利用multiwfn画的结果中的HOMO虚线对应的值确不是这个值,为什么?[图片]

  1121. A:
  1122.     估计你绘制的是alpha的

  1123. ----------------------------------------------------
  1124. 2016.06.25 18:53:46
  1125. Q:
  1126.     老师,在Gv里能直接做两个分子之间HOMO-HOMO的轨道重叠吗?我刚才试了一下,总是退出来。

  1127. A:
  1128.      分子间轨道重叠的图形显示和计算
  1129.     http://sobereva.com/163

  1130. Q:
  1131.     老师,我是想做这种图[图片]

  1132. A:
  1133.      博文里写得很清楚了

  1134. Q:
  1135.     好的,老师,我正在画图呢
  1136.     谢谢老师
  1137.     Summing up all value and multiply differential element老师,我到最后没出现这行数信息啊?

  1138. A:
  1139.     仔细看屏幕,肯定有
  1140.     [图片]

  1141. ----------------------------------------------------
  1142. 2016.06.25 20:40:34
  1143. Q:
  1144.      老师,昨天的问题你还没回答我呢
  1145.     【传说】娃娃-首席卖萌官(979170845) 19:31:17
  1146.      老师,画超极化率密度图的时候,需要计算单点的水平要和计算超极化率的水平一致吗?

  1147. A:
  1148.      需要

  1149. ----------------------------------------------------
  1150. 2016.06.26 01:05:30
  1151. Q:
  1152.     我想请教一下   高斯能不能算出PMF曲线
  1153.     之前一直在做密度泛函    没有弄过动力学的东西

  1154. A:
  1155.      不能。PMF是分子动力学程序做的

  1156. ----------------------------------------------------
  1157. 2016.06.26 05:07:11
  1158. Q:
  1159.     [图片]
  1160.     骗子

  1161. A:
  1162.     骗子
  1163.     正规培训没有这么偷偷摸摸发通知的,肯定有正规的培训通知网页

  1164. ----------------------------------------------------
  1165. 2016.06.26 05:10:41
  1166. Q:
  1167.     请教老师,chemoffice 没有序列号啊,怎么注册呢

  1168. A:
  1169.     有专门注册机
  1170.     破解版的要么带注册机,要么带破解文件,即便压缩包里没明确告诉你怎么破解,稍微一动脑也知道怎么用注册机或破解文件

  1171. ----------------------------------------------------
  1172. 2016.06.26 05:34:58
  1173. Q:
  1174.     [图片]
  1175.     这个是HOMO轨道,我们能从这个图上得到哪些信息呢

  1176. A:
  1177.         那图没用
  1178.     光有个轨道图什么也说明不了

  1179. Q:
  1180.     回复 Sobereva:我们这个图上能读出哪些信息?比如电荷啊什么的?

  1181. A:
  1182.     只给HOMO分布而已,不结合具体问题什么也说明不了
  1183.     计算化学不是先算出一些数据然后才去想办法给他找意义

  1184. Q:
  1185.     我希望能考察反应活性,不知道这样是否合理?回复 Sobereva:

  1186. A:
  1187.     幻灯片:Predicting reactive sites (反应位点的预测)
  1188.     http://sobereva.com/234
  1189.     《亲电取代反应中活性位点预测方法的比较》doi: 10.3866/PKU.WHXB201401211
  1190.    
  1191.     看这些

  1192. ----------------------------------------------------
  1193. 2016.06.26 06:13:26
  1194. Q:
  1195.     [图片]field=Z+60  这个是静态电场的强度吗?

  1196. A:
  1197.     y

  1198. Q:
  1199.     那如果不加,是默认的?是多少呢?

  1200. A:
  1201.     0

  1202. Q:
  1203.     默认是不加电场的?

  1204. A:
  1205.     默认加什么电场

  1206. Q:
  1207.     #p BMK/gen polar=(DCSHG) int=(ultrafine) scrf=(smd,solvent=toluene)
  1208.     那我这样的输入,是加了多大的电场呢?方向呢?搞糊涂了
  1209.     # HF/aug-cc-pVTZ field=Z+60 out=wfx
  1210.    
  1211.     B3LYP/TZVP opted
  1212.    
  1213.     0 1
  1214.     C                  0.00000000    0.00000000   -0.52710800
  1215.     H                  0.00000000    0.93885600   -1.11413900
  1216.     H                  0.00000000   -0.93885600   -1.11413900
  1217.     O                  0.00000000    0.00000000    0.67386600
  1218.    
  1219.     C:\H2CO-2FZ.wfx
  1220.     ************
  1221.     #p PBE1PBE/aug-cc-pVTZ polar CPHF=RdFreq
  1222.    
  1223.     test
  1224.    
  1225.     0 1
  1226.     C                 -0.55391731    0.43227932    0.14513431
  1227.     H                 -1.29461381   -0.13844703    0.71478821
  1228.     H                 -0.40071988    1.37810560    0.67504812
  1229.     H                 -0.99858078    0.66864399   -0.83680924
  1230.     N                  0.70741247   -0.31205725    0.11165308
  1231.     H                  0.57479710   -1.18908235   -0.38659422
  1232.     H                  1.39864521    0.20863272   -0.42322026
  1233.    
  1234.     0.07 0.1
  1235.      这两个加外场的方法有什么区别?sob老师

  1236. A:
  1237.     polar算的是(超)极化率,这是分子内在属性(外场为0时候能量对电场的导数),和field无关,除非你自行用有限场来计算
  1238.     field=Z+60算的是带静态外场时候的情况,不是算(超)极化率
  1239.    
  1240.     后面是计算含频(超)极化率

  1241. Q:
  1242.     我再理解一下,谢谢Sob!
  1243.     DCSHG 对 dc-SHG(直流二次谐波生成)超极化率执行额外的含频 CPHF 计算。该选项 还意味着 CPHF=RdFreq。
  1244.     也就是说:polar=(DCSHG)和polar=(DCSHG) CPHF=RdFreq是一样的,写了polar=(DCSHG)就意味着CPHF=RdFreq,是吗?

  1245. A:
  1246.     y

  1247. Q:
  1248.     [图片]末尾的0.0是否可以省略?

  1249. A:
  1250.     try
  1251.        
  1252. ----------------------------------------------------
  1253. 2016.06.26 06:33:36
  1254. Q:
  1255.     各位前辈好!能不能请教一个问题: 优化中间体出现Small interatomic distances encountered:
  1256.     2    1
  1257.    
  1258.     3    1
  1259.    
  1260.     3    2
  1261.    
  1262.     4    1
  1263.    
  1264.     4    2
  1265.    
  1266.     4    3
  1267.    
  1268.     5    1
  1269.    
  1270.     5    2
  1271.    
  1272.     5    3
  1273.    
  1274.     5    4
  1275.    
  1276.     6    1
  1277.    
  1278.     6    2
  1279.    
  1280.     6    3
  1281.    
  1282.     6    4
  1283.    
  1284.     6    5
  1285.    
  1286.     7    1
  1287.    
  1288.     7    2
  1289.    
  1290.     7    3
  1291.    
  1292.     7    4
  1293.    
  1294.     7    5
  1295.     Problem with the distance matrix.
  1296.     Error termination via Lnk1e in /bin/g03e/g03/l202.exe at Fri Jun 24 21:12:06 2016.
  1297.     Job cpu time:  1 days  0 hours 12 minutes 46.4 seconds.
  1298.     但是出错的上个结构是对的,不知道为什么距离会突然变小,该如何解决呀?
  1299.     网上说加nosymm,不知道原因是什么

  1300. A:
  1301.     跟nosymm根本无关
  1302.     用g09

  1303. ----------------------------------------------------
  1304. 2016.06.26 06:47:00
  1305. Q:
  1306.     sob老师,我计算两个甲胺吸收一个二氧化碳的过渡态时,总是出现l9999不收敛的情况,也按照您帖子里给出的方法尝试了,还是不收敛,怎么办啊?opt=(calcfc,cartesian,ts,noeigen) int=ultrafine freq b3lyp/6-31g(d)

  1307. A:
  1308.     改小步长上限
  1309.     尝试gdiis
  1310.     甭用cartesian
  1311.     还不行就calcall
  1312.     没必要用ultrafine

  1313. Q:
  1314.     cartesian是在过渡态中有三个分子程180度的时候才用是吗

  1315. A:
  1316.     没报错就甭用它,只会让收敛更难

  1317. ----------------------------------------------------
  1318. 2016.06.26 07:43:36
  1319. Q:
  1320.     老师,怎么给高斯添加新的基组?比如lanl2tz,def2-tzvp

  1321. A:
  1322.     详解Gaussian中混合基组、自定义基组和赝势基组的输入
  1323.     http://sobereva.com/60
  1324.    
  1325.     def2TZVP本来就是内置的

  1326. ----------------------------------------------------
  1327. 2016.06.26 07:59:08
  1328. Q:
  1329.     wb97xd/TZVP
  1330.     这样写方法基组对吗?

  1331. A:
  1332.     对

  1333. ----------------------------------------------------
  1334. 2016.06.26 17:22:55
  1335. Q:
  1336.     我宁可以后我的办公室的电脑都装红帽
  1337.     那样科研的潜移默化更快
  1338.     乌班图不好么
  1339.     社长不推荐
  1340.     [表情]如果是单机,乌班图可以一试,不过我最近也发现,越来越花哨的乌班图对于科研没什么卵用

  1341. A:
  1342.     ubuntu太轻浮

  1343. Q:
  1344.     [表情]不过追求华丽,可以试试
  1345.     Soga,看来要换红帽了
  1346.     红帽不是收费吗?免费的是CentOs吧
  1347.     嗯嗯。不过网上有红帽的破解版
  1348.     6.5版的
  1349.     这都可以破解~~

  1350. A:
  1351.     红帽用不着破解,只不过不注册没法用官方的源,借用CentOS的即可

  1352. Q:
  1353.     只是个安装码

  1354. A:
  1355.     安装都不要密码

  1356. ----------------------------------------------------
  1357. 2016.06.26 17:52:36
  1358. Q:
  1359.     请问,能不能直接将高斯优化的结果文件比如log或者chk文件转换为wfn文件呢
  1360.     加关键字的话意味着需要重新计算吧
  1361.     有没有命令转换两个文件的

  1362. A:
  1363.     chk转换为fch就直接能给multiwfn用了
  1364.     log不可能转换为wfn
  1365.     multiwfn本身也有fch->wfn转换功能

  1366. ----------------------------------------------------
  1367. 2016.06.26 17:56:25
  1368. Q:
  1369.     GTBas 1,2,GTMP2largeXP,GFHFB3,4都是什么类型的基组
  1370.     GTBas 1,2,GTMP2largeXP,GFHFB3,4都是什么类型的基组

  1371. A:
  1372.     甭用这些烂七八糟的,净是热力学组合方法里定义的稀奇古怪的。
  1373.     不适合一般目的计算

  1374. ----------------------------------------------------
  1375. 2016.06.26 18:08:34
  1376. Q:
  1377.     Sob老师,您好。 我在 《International Journal of Oil, Gas and Coal Technology》期刊投了一篇文章,名为《Weak interaction between water molecules and different rank coals: A DFT-D3 study》,我已推荐您为审稿人,希望得到您的指点与建议。@Sobereva

  1378. A:
  1379.     恩

  1380. Q:
  1381.     /握手
  1382.     怎么不给SOB发个红包表示诚意?!

  1383. A:
  1384.     岂不成了贿赂审稿人

  1385. ----------------------------------------------------
复制代码



评分 Rate

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

查看全部评分 View all ratings

本版积分规则 Credits rule

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

GMT+8, 2025-8-15 22:01 , Processed in 0.210470 second(s), 24 queries , Gzip On.

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