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

gangyecp · 发表于 Post on 2018-6-27 08:41:47

sobereva 发表于 2018-6-26 22:57
基于过渡态理论计算反应速率常数的Excel表格
http://sobereva.com/310
使用KiSThelP结合Gaussian基于过 ...

谢谢老师。老师还有个问题，就是利用老师您提供的表格计算单分子反应速率时，随着温度的升高，计算出的反应速率增大；而计算双分子反应速率时，随着温度的升高，计算出的反应速率降低。在实际反应中，反应速率应随着温度的升高而增大，也是按照老师提供的表格将自己求得的数据代入的，不知道什么原因，双分子的时候得到的与实际情况相反。还希望老师给予指导下，谢谢老师。

菩城 · 发表于 Post on 2018-6-27 11:19:28

sob老师好，您在回答一个动力学半径的问题提到：“动力学直径本身是实验上测的量，要想理论计算，就得做分子模拟来重现实验过程，也就是搭建不同孔径的分子筛通过模拟看分子能否过去”。我的疑问是：我计算的孔材料尺寸如何跟动力学直径关联起来？比如CO2的动力学直径是3.30埃，我怎么判断我的孔材料的尺寸是否大于CO2的动力学直径呢？我在文献上看到有人自己定义一个尺寸来比较，是否有比较普适的定义方法呢？谢谢！

sobereva · 发表于 Post on 2018-6-27 14:13:33

菩城发表于 2018-6-27 11:19
sob老师好，您在回答一个动力学半径的问题提到：“动力学直径本身是实验上测的量，要想理论计算，就得做分 ...

这不属于量化问题，请单独发帖

sobereva · 发表于 Post on 2018-6-27 14:16:01

gangyecp 发表于 2018-6-27 08:41
谢谢老师。老师还有个问题，就是利用老师您提供的表格计算单分子反应速率时，随着温度的升高，计算出的反 ...

没有具体信息没法回答。
假设不考虑透射系数，随温度升高，用我那个excel表格来算k一定是和温度正相关的
倘若你考虑了透射系数才出现你说的问题，那可能是因为温度越低隧道效应越明显，kappa越大所致

菩城 · 发表于 Post on 2018-6-27 15:03:28

sobereva 发表于 2018-6-27 14:13
这不属于量化问题，请单独发帖

好的，谢谢老师！

gangyecp · 发表于 Post on 2018-6-27 16:41:28

sobereva 发表于 2018-6-27 14:16
没有具体信息没法回答。
假设不考虑透射系数，随温度升高，用我那个excel表格来算k一定是和温度正相关 ...

好的，谢谢老师，明白了

she · 发表于 Post on 2018-6-27 16:52:26

请问各位老师，gaussian里面算偶极矩是怎么算的？个人理解是先求马利肯电荷分布，再求正负电荷中心之类，这样算的。不知道对不对。
先谢谢回答

sobereva · 发表于 Post on 2018-6-27 17:46:15

she 发表于 2018-6-27 16:52
请问各位老师，gaussian里面算偶极矩是怎么算的？个人理解是先求马利肯电荷分布，再求正负电荷中心之类，这 ...

完全不对
Mulliken电荷对偶极矩的重现性极差，只有我提出的ADCH原子电荷才能精确重现基于波函数得到的偶极矩

Gaussian算完了（不管是何种任务），输出文件末尾直接就给你当前结构下的偶极矩了，是基于波函数算的，用这个就完了

qianbao · 发表于 Post on 2018-6-28 09:09:40

我在做优化（图opt）、算激发态（图td）、以及激发态的优化（图tdopt）时，发现其输出构型有很大的变化，优化激发态的时候构型怎么会变这么多呢？请问老师，这个可能是由于色散力太大吗？

329678001 · 发表于 Post on 2018-6-28 09:11:08

本帖最后由 329678001 于 2018-6-28 09:12 编辑

sob老师您好，1. 我们算气体吸附能时候，应该是用的各自基态来做运算的。计算金属吸附能时，我导师说公式中金属的价态是由体系的价态决定，这个可以么？有相关的文献么，我看文章里没有这么详细的计算细节。
2. 我最近在研究金属价态的不同，所导致的对同一气体吸附能影响，不知道有没有一定相关性，想让您帮我推荐2篇文献。

CCCC · 发表于 Post on 2018-6-28 09:53:52

Requested convergence on RMS density matrix=1.00D-08 within 128 cycles.
Requested convergence on MAX density matrix=1.00D-06.
Requested convergence on          energy=1.00D-06.
No special actions if energy rises.
Integral accuracy reduced to 1.0D-05 until final iterations.
Problem detected with inexpensive integrals.
Switching to full accuracy and repeating last cycle.
Spurious integrated density or basis function:
NE=  456 NElCor= 0 El error=1.29D+00 rel=1.33D-03 Tolerance=1.00D-03
Shell 709    absolute error=2.75D-04             Tolerance=1.20D-02
Shell 709    signed error=2.75D-04             Tolerance=1.00D-01
Inaccurate quadrature in CalDSu.
Error termination via Lnk1e in /home/ftcc/g09B01/l502.exe at Wed Jun 27 22:50:27 2018.

she · 发表于 Post on 2018-6-28 09:55:16

sobereva 发表于 2018-6-27 17:46
完全不对
Mulliken电荷对偶极矩的重现性极差，只有我提出的ADCH原子电荷才能精确重现基于波函数得到的 ...

先非常感谢sob老师回答。
sob老师，就是想算一对阴阳离子之间的偶极矩，在用高斯# opt b3lyp/6-31+g（2df,p） geom=connectivity算出来的偶极矩是将这个体系看做一个整体再做电荷分布求得的偶极矩吗？
如果想求阴离子和阳离子之间偶极矩的矢量和，这样怎么算？
或者这二者有区别吗？

app · 发表于 Post on 2018-6-28 10:00:49

老师，您好，请教您一个问题，我想看一个复合物中，一个分子的带电情况，在gaussview中打开输出文件之后可以怎么操作，就可以看到局部分子所带的总电荷是多少了呢。谢谢您！

zjxitcc · 发表于 Post on 2018-6-28 10:34:16

唉，这些论坛上早就有了，自己搜索一下“Inaccurate quadrature in CalDSu”立即找到答案，何必发帖。
看《Gaussian FAQ；新手求助报错时的注意事项》http://bbs.keinsci.com/thread-4829-1-1.html

CCCC · 发表于 Post on 2018-6-28 10:49:28

zjxitcc 发表于 2018-6-28 10:34
唉，这些论坛上早就有了，自己搜索一下“Inaccurate quadrature in CalDSu”立即找到答案，何必发帖。
看 ...

好的谢谢

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

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

算能量的时候出现L502错误