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

sobereva · 发表于 Post on 2017-9-22 15:49:58

wangjianmin 发表于 2017-9-22 10:49
老师，又来麻烦您了，orca提示 Cannot open input file,是文件路径不对吗？

没写对路径

shangfang · 发表于 Post on 2017-9-23 09:26:17

sob老师，你好，请问如何计算自旋多重度是2的激发能，我查看高斯说明书发现不管TD还是eom-ccsd都只有singlets 和triplets的关键词，如果这个物质是带电荷的二重态如何计算激发能呢。

dingniu2 · 发表于 Post on 2017-9-23 10:20:27

请问老师，氢键键能估算用BCP势能密度的一半只适用于X-H----O (X=C,N,O)型氢键吗？如果是其它类型，如X-H----N (X=C,N,O)等适用吗？

computational · 发表于 Post on 2017-9-23 15:54:27

请问sob老师，（1）如果mayer键级在0.5左右算成键吗？
（2）对于分析成键情况，RDG与ELF相比是不是没有任何优势？
（3）文献中ELF比LOL更常见，是不是分析成键ELF更好？

renlangfen · 发表于 Post on 2017-9-23 19:00:21

各位大神，求问如果计算溶剂化，高斯手册中没有该溶剂，该如何处理，感激不尽。

happyknighthawk · 发表于 Post on 2017-9-23 19:06:42

请问Sob老师
因为要计算权重光谱或核磁，需要提取各个构象frequency分析后的热力学数据，现在想通过一个脚本把Sum of electronic and thermal Free Energies这个数值批量提取出来，所有按照【log文件名 + gibbs自由能】的格式输出到一个文件里，方便后续转到ECXEL里面进行波尔兹曼分布的计算。仿照公社里之前有位楼主的帖子写了如下脚本：
#!/bin/bash
for m in *.log
do
echo searching ${m} ...
name=${m%.*}
grep 'Sum of electronic and thermal Free Energies=' ${m} > ${name}.txt
echo ${m} completed
echo
done
但是都是一个log输出一个txt。想请问您，如何把所有log里面的能量按照特定格式输出到一个文件里，要使用什么命令？
谢谢您

xylz6188 · 发表于 Post on 2017-9-23 19:24:24

自定义，小木虫有，Gaussian手册里也有讲解
http://muchong.com/t-11007750-1，本论坛也有，我前两天也看到过，自己再查查

renlangfen · 发表于 Post on 2017-9-23 20:10:05

xylz6188 发表于 2017-9-23 19:24
自定义，小木虫有，Gaussian手册里也有讲解
http://muchong.com/t-11007750-1，本论坛也有，我前两天也看 ...

环氧丙烷找了好久哦

renlangfen · 发表于 Post on 2017-9-23 20:38:06

xylz6188 发表于 2017-9-23 19:24
自定义，小木虫有，Gaussian手册里也有讲解
http://muchong.com/t-11007750-1，本论坛也有，我前两天也看 ...

您好，我想知道化合物的ε在哪边查，本人菜鸟，望大神解答

renlangfen · 发表于 Post on 2017-9-23 20:58:03

由于高斯自带溶剂中没有环氧丙烷，小弟翻遍百度都找到环氧丙烷介电常数，求大神指点迷津，计算溶剂化，感激不尽。

sobereva · 发表于 Post on 2017-9-23 21:29:40

renlangfen 发表于 2017-9-23 20:58
由于高斯自带溶剂中没有环氧丙烷，小弟翻遍百度都找到环氧丙烷介电常数，求大神指点迷津，计算溶剂化，感激 ...

去搜CRC Handbook of Chemistry and Physics、兰氏化学手册、Google，搜不到就用环氧乙烷凑合，或者自行做动力学模拟计算

sobereva · 发表于 Post on 2017-9-23 21:30:33

computational 发表于 2017-9-23 15:54
请问sob老师，（1）如果mayer键级在0.5左右算成键吗？
（2）对于分析成键情况，RDG与ELF相比是不是没有任 ...

算
一般来说是
未必，有的体系用LOL的图更清楚

sobereva · 发表于 Post on 2017-9-23 21:32:16

dingniu2 发表于 2017-9-23 10:20
请问老师，氢键键能估算用BCP势能密度的一半只适用于X-H----O (X=C,N,O)型氢键吗？如果是其它类型，如X-H-- ...

其它类型未必适用，除非自行找体系测试一下证明可用性

sobereva · 发表于 Post on 2017-9-23 21:32:58

shangfang 发表于 2017-9-23 09:26
sob老师，你好，请问如何计算自旋多重度是2的激发能，我查看高斯说明书发现不管TD还是eom-ccsd都只有single ...

直接写TD就完了

sobereva · 发表于 Post on 2017-9-23 21:39:08

happyknighthawk 发表于 2017-9-23 19:06
请问Sob老师
因为要计算权重光谱或核磁，需要提取各个构象frequency分析后的热力学数据，现在想通过一个脚 ...

这种问题我都是直接写Fortran程序做...

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

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

各位大神，求问如果计算溶剂化，高斯手册中没有该溶剂，该如何处理

想要计算溶剂化，小弟翻遍百度都没找到环氧丙烷介电常数，求大神指点.