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

dingniu2

sob老师，你好：一个归一化的波函数乘以绝对值不为1的常数就不是归一化的了，请问老师波函数乘以-1所代表的量子态与原波函数表示的量子态是否完全一样？为什么？谢谢老师！

ccccn345

请问老师，要得到HOMO和LUMO是不是应该计算分子的基态，计算激发态得到的HOMO和LUMO有没有意义？

jingruoyoulan

请问sob老师，我想构建K2O（110）裸簇模型，首先在MS里做完110切面（1层），添加好真空层，将切面原子移动到真空层内，构建2*2*1的超胞，然后导出为.cif文件，用VESTA转换为pdb文件，VESTA转换步骤：edit--bonds--Do not search atoms beyond the boundary……接着用高斯做SP计算能正常结束，opt计算SCF不收敛，K2O（110）裸簇模型可以拿来直接计算吗？还是必须要经过进一步加工？请sob老师指点，谢谢~~~

yangjianfeng

sob老师  量子化学零基础的人（研一，其他系调剂到化学系的）  学量化该用哪本教材呢

Wind_runner

yangjianfeng 发表于 2017-11-7 19:24
sob老师  量子化学零基础的人（研一，其他系调剂到化学系的）  学量化该用哪本教材呢

levine的第七版量子化学

Wind_runner

计算寻找过渡态出现l502错误   求解哪里出了错
%nprocs=4
%mem=1GB
%chk=ts.chk
# UHF/6-31G(d) Opt=QST2 test

r

0 2
C                 -0.00536900   -0.61780300    0.00000000
H                 -0.48144300   -1.00815500    0.88911300
H                 -0.48144300   -1.00815500   -0.88911300
H                  1.03805900   -0.91696800    0.00000000
O                 -0.00536900    0.83001200    0.00000000

p

0 2
C                 -0.69842700   -0.03120700    0.06942400
H                 -1.13557900   -0.97393800   -0.20083500
H                 -1.22962000    0.88369300   -0.07652600
O                  0.67273300    0.11869800   -0.03296100
H                  1.17389300   -0.67209500    0.12450200

错误
Annihilation of the first spin contaminant:
S**2 before annihilation     1.1923,   after     0.7599
Convergence failure -- run terminated.
Error termination via Lnk1e in /share/apps/gaussian/g09D01/g09/l502.exe at Tue Nov  7 19:17:31 2017.
Job cpu time:       0 days  0 hours  0 minutes  7.4 seconds.
File lengths (MBytes):  RWF=      5 Int=      0 D2E=      0 Chk=      1 Scr=      1

wrtgcn

Wind_runner 发表于 2017-11-7 20:12
计算寻找过渡态出现l502错误   求解哪里出了错
%nprocs=4
%mem=1GB

http://bbs.keinsci.com/forum.php ... 829&highlight=F

wrtgcn

Sob老师，最近有几个想不通的问题，还得请教您一下：
1.石墨烯及C3N4体系掺杂或改性的计算类文献中，往往是直接将杂原子替换进去。但是譬如将N原子替入石墨烯，应该会造成周围C原子的不饱和（类似的还有掺杂引起断键，也会造成相邻原子不饱和）；而实际情况下，存在不饱和原子的体系应该是不稳定的。那这样的计算体系是否合理？
2.同样是这类体系下，往往是3重态比单重态更稳定（验stable），是否可以理解为不饱和原子处的电子难以与系统内其它电子配对所致？
3.人为补足不饱和原子（比如加上H原子），与保持原体系不变，哪个更为合理？
之前一直是照着文献进行的计算，突然想到了这些，难以得出结论。还请老师指点。

yangjianfeng

sob 老师 麻烦了 请问下①由图一i中3.31怎样就得到了图二3.32
                                ②图二中的3.32后2项只是σ位置变了  怎么就能直接写成3.33式中的复共轭了

steven

sob老师好。想请教您一下怎么把理论计算的结果与真实反应条件关联起来。
一.例如我如果通过理论计算找出了一条反应途径，怎么知道：现实中是否成立，试剂浓度和怎么选择溶剂，以及每一步基元反应（或者总包反应）的条件？
二.怎么体现反应可以进行到什么程度？
在回家的火车上，听完这次课觉得很兴奋。问题有点多，谢谢老师了！

卡开发发

dingniu2 发表于 2017-11-6 09:32
sob老师，你好：一个归一化的波函数乘以绝对值不为1的常数就不是归一化的了，请问老师波函数乘以-1所代表的 ...

一样。比如一个系统的态ψ在r1和r2两处的相对概率并不会因为乘以常数而发生变化，但按照波函数统计解释，一般会要求概率归一化。

卡开发发

jingruoyoulan 发表于 2017-11-6 22:48
请问sob老师，我想构建K2O（110）裸簇模型，首先在MS里做完110切面（1层），添加好真空层，将切面原子移动 ...

你可以参考这两篇帖子以及下面的讨论：
http://bbs.keinsci.com/forum.php ... mp;highlight=%B4%D8
http://bbs.keinsci.com/forum.php?mod=viewthread&tid=6755
包含的内容已经比较丰富了。

卡开发发

yangjianfeng 发表于 2017-11-8 15:09
sob 老师 麻烦了 请问下①由图一i中3.31怎样就得到了图二3.32
                                ②图二中 ...

下次发图有条件的话还是尽量转过来吧，要不然大家看着都累。

详细过程我写在下面，但建议最好自己推导一遍，后续推导Hartree-Fock近似和Kohn-Sham方程的手段照葫芦画瓢。

有可能会有笔误，哪怕是书上的也会有，比如3.33的求和的[]括错了地方，所以请多包涵。

meatball1982

关于compliance对角线的输出。
根据：
http://sobereva.com/364
我希望通过 compliance得到一个分子中的若干键的力常数。



帖子中说，各个键的力常数是通过对角线的元素的倒数得到的。

我的问题是：
这个表没法复制，自己一个一个输入倒是行。有点麻烦。
有没有能导出这个表格的方法。

谢谢大家。

yangjianfeng

卡开发发 发表于 2017-11-9 02:24
下次发图有条件的话还是尽量转过来吧，要不然大家看着都累。

详细过程我写在下面，但建议最好自己推导 ...

以后发图会注意的，非常感谢 写的不能更仔细了。    另外3中的共轭能简单科普下嘛（没学过）或者说推荐点书或者博文