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

liyuanhe211 · 发表于 Post on 2016-3-19 12:25:12

alexwpch 发表于 2016-3-19 11:33
请教Sob,
我尝试使用ORCA，照网上模板套的输入文件，但出现“ERROR: expect a '$', '!', '%', '*' or '[' ...

经测试没有返回你说的错误。
不过你这原子坐标的元素符号全都是“X”是什么鬼

alexwpch · 发表于 Post on 2016-3-19 14:46:50

sobereva 发表于 2016-3-19 12:21
些简单点，比如

! PBE0 def2-TZVP miniprint TightSCF pal4

目前版本不能? 为嘛？以后还支持IGLO吗？那有没有IGLO package? 我问的有点多哦。。。。

alexwpch · 发表于 Post on 2016-3-19 14:47:46

liyuanhe211 发表于 2016-3-19 12:25
经测试没有返回你说的错误。
不过你这原子坐标的元素符号全都是“X”是什么鬼

那能帮我算下吗，X改成N就行，谢谢。

498746012 · 发表于 Post on 2016-3-19 21:59:05

本帖最后由 498746012 于 2016-4-8 14:40 编辑

sob老师：
您好！

sobereva · 发表于 Post on 2016-3-20 02:33:52

498746012 发表于 2016-3-19 21:59
sob老师：
您好！能不能请教您两个问题：一是，下图是Ti和N之间的分子轨道（用B3LYP计算）HOMO,HOM ...

1 用Multiwfn计算Mayer键级就行了，看Multiwfn手册4.9.1节的例子
2 算原子电荷，比如高斯里可以用pop=NPA得到NPA原子电荷，是很常用的

lei234 · 发表于 Post on 2016-3-20 02:50:32

sobereva 发表于 2016-3-19 12:08
SMD里极性部分用的IEFPCM，高斯里PCM就等于IEFPCM，直接用PCM没有考虑非极性部分的贡献

谢谢您的回答~
那么请问使用SMD的时候关键词是这样写的吗：scrf=(solvent=water,smd) ？，读自由能的时候找这项是吗？Sum of electronic and thermal Free Energies，直接读不需要任何的数据处理吗？

498746012 · 发表于 Post on 2016-3-20 08:53:32

sobereva 发表于 2016-3-20 02:33
1 用Multiwfn计算Mayer键级就行了，看Multiwfn手册4.9.1节的例子
2 算原子电荷，比如高斯里可以用pop=NP ...

谢谢sob老师您在深夜里还能解答我的问题

sobereva · 发表于 Post on 2016-3-20 11:05:54

lei234 发表于 2016-3-20 02:50
谢谢您的回答~
那么请问使用SMD的时候关键词是这样写的吗：scrf=(solvent=water,smd) ？，读自由能的时 ...

是
这是比较粗略的溶液下的自由能，要是想准确考虑，要把高精度单点能、合适级别下计算的溶解自由能、一般级别下考虑频率校正因子得到的自由能校正能加起来。

15whyhq · 发表于 Post on 2016-3-20 16:36:54

Symmetry consistency failure.
Error termination via Lnk1e in /public/software/Gaussian09/g09B01/l202.exe

请问高斯输出文件显示l202错误该怎么解决呀？谢谢！

sobereva · 发表于 Post on 2016-3-20 18:15:09

15whyhq 发表于 2016-3-20 16:36
Symmetry consistency failure.
Error termination via ...

加上nosymm看能否解决，不行的话把输入输出文件传上来

lei234 · 发表于 Post on 2016-3-21 01:31:07

sobereva 发表于 2016-3-20 11:05
是
这是比较粗略的溶液下的自由能，要是想准确考虑，要把高精度单点能、合适级别下计算的溶解自由能、 ...

我还是粗略得了，非常谢谢您~

lei234 · 发表于 Post on 2016-3-21 02:54:00

请问，加了溶剂效应的自由能用单点计算后发现有很多虚频，是不是这种情况的能量不能用？是不是所有的单点计算都有虚频？

lei234 · 发表于 Post on 2016-3-21 02:58:39

我加了freq在算单点的时候，因为要读Sum of electronic and thermal Free Energies

lei234 · 发表于 Post on 2016-3-21 03:09:12

sobereva 发表于 2016-3-20 11:05
是
这是比较粗略的溶液下的自由能，要是想准确考虑，要把高精度单点能、合适级别下计算的溶解自由能、 ...

勒个。。。我要得到溶剂效应下的吉布斯自由能好像只能按照你说的这个了，因为我算单点的时候加了freq发现有好多虚频，看见网上说算单点不能加freq，加了也不对。。所以，请问您说的这三项是不是这样的？
高精度单点能--就是加了溶剂效应后，单点计算最后得到的HF？
合适级别下计算的溶解自由能--这个是什么？怎么算的？读的那个能量？
一般级别下考虑频率校正因子得到的自由能校正能--请问是不是气相下的这个值：Thermal correction to Gibbs Free Energy？
还有，您说的一般级别和合适级别是指的不同基组吗？我在气相下用的M062X/6-31+g**,液相下也是这个，只是加了溶剂效应，如果要对比气相和液相下能垒的高低，不是应该保持方法和基组一样吗？

期待您的回复，谢谢~

498746012 · 发表于 Post on 2016-3-21 08:44:01

sob老师：
您好！进行NPA原子电荷计算，是不是只在#行加上pop=NPA,其他地方和溶剂校正的一样？如何在gview中看NPA？如果中心金属的NPA是2左右那么就是+2价，如何是4左右，就是+4价，这样理解对吗？希望老师您指教

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

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

高斯输出文件显示l202错误怎么解决

浏览过的版块