ggdh 发表于 2017-2-3 17:39 感谢回答~~ 感觉后面更换基团的做法1 是不是脱离了波函数分析的范畴，像是有些文献里找一堆不同取代基的底物算一遍、以某个表征位阻的参数做横坐标、能垒做纵坐标画条线的那种研究。我思考一下做法2如何实现。 看你的回复忽然想，一般所说位阻（例如两个叔丁基之间接近）的本质是不是以vdW为主？如果是的话使用完全不能描述vdW的方法做起点是不是可以提取出位阻的影响。 |
|
本帖最后由 liyuanhe211 于 2017-2-4 04:06 编辑 我明白了~~谢谢Sob老师 。 我想想怎么方便的实现它 |
liyuanhe211 发表于 2017-2-3 10:03 是的 |
liyuanhe211 发表于 2017-2-3 09:46 ‘位阻’可能包括静电，vdw或者轨道作用。如果位阻不大，一般不会涉及到轨道作用。 如果这个位阻没有轨道作用。那么能量分解的方法可以直接把这种位阻作用剔除。 如果这个位阻包括轨道作用。那么可以单独计算位阻基团的轨道作用。然后从总的轨道作用中扣除。 有两个办法： 1. 把位阻基团去掉，然后剩下的部分做能量分解。比较去掉位阻基团前后的能量分解各项。如果把其中轨道项的差别认为是位阻作用中轨道项的贡献。（如果位阻基团对共轭体系电子结构影响大。就不推荐这种方法） 2. 单独考虑位阻基团。单独把位阻基团弄出来，保持几何结构不变做能量分解。看看其中轨道作用大不大。如果轨道作用大，可以把这个轨道作用能记录下来。然后从总的轨道作用能中扣除位阻基团的轨道作用能。剩下的就认为是secondary interaction 可以上面两个方法同时算，然后相互印证。 |
本帖最后由 wswrpd 于 2017-2-3 12:57 编辑 liyuanhe211 发表于 2017-2-3 10:03 对了，除了GMH；其他两种方法都是基于DFT，对你而言不知道有没有用。所以除了GMH我也想不出更好的方法算，如果用GMH的话就要直接计算整个Supermolecule的能级结构，然后进行分析，  另外，感觉sob老师说的和片段组合波函数的某些想法有点关联，但具体我也不懂啊 |
本帖最后由 wswrpd 于 2017-2-4 12:24 编辑 liyuanhe211 发表于 2017-2-3 10:03 我觉得是这个意思，问题是你知道具体怎么用Gaussian Script 来实现这个计算？ 你可以看我曾经发过的帖子，也是类似的，如何计算两个Diabatic Local Excited State之间的Coupling的问题，目前我觉得利用Gaussian的话只能用GMH方法来算了（用adiabatic算出的数据去反算diabatic的Coupling），更好的方法比如CDFT，或者frozen density-embedding (FDE) formalism都需要对K-S方程进行额外的约束，所以必须用Q-Chem或者ADF这种软件才有现成的程序。 以上的方法主要用来处理 Acceptor和Doner之间的Electron Transfer过程，和D-A反应有一定区别，还请留意 如果想出来怎么做一定要发帖介绍经验啊 |
liyuanhe211 发表于 2017-1-31 02:28 因为前线轨道理论本身就是对于反应初期，将分子间二阶微扰作用能进行简化来得到的。如果不只有一对轨道间相互作用很明显，那就最好直接通过二阶微扰公式来明确考察轨道间作用对体系能量降低的贡献并进行横向比较。对于TS结构，两个单体各自的特征已经很大程度被破坏了，片段间呈现共价作用，相互作用已经很强了，就不适合再通过微扰方式考察了，这时就需要做SCF了。 我说的考查方式和复合物做普通的MP2这种后HF计算是两回事。倒是和考察分子间电荷转移要计算的<A_i|H|B_a>这种积分对应，其中A_i是A的某占据轨道，B_a是B的某空轨道，H是复合物的单电子哈密顿算符（可以是HF或KS算符），对应于两个片段直接放在一起，电子密度还没驰豫时的情况（相当于片断组合波函数当复合物初猜时的第一轮SCF的哈密顿）。 |
本帖最后由 ggdh 于 2017-2-1 20:15 编辑 liyuanhe211 发表于 2017-1-31 02:15 想了一下，这个做cda分析没有意义，因为cda只能看片段轨道系数，但是片段轨道系数大不代表作用强（比如离得比较远的有简并轨道的片段）。而且删除primary interacting 原子上的系数也不会改变片段轨道系数。所以1号方法作废。 secondary interaction 这个概念本身不会有严格的定义。比如这里就涉及到轨道的实空间怎么划分的问题。 我觉得比较好的是3号办法。因为三号办法里面已经做了能量分解。相当于把轨道作用能单独抽提出来了。那么你所说的‘掺杂了除 secondary interaction 之外的其他成分’也只能是其他的轨道作用了。按照你画的图中表示，红色和红色轨道作用是primary，蓝色和蓝色是secondary，那么剩下的就是红色和蓝色的之间的 tertiary作用了。。。。不过反正你也是考虑endo选择性。可以把primary以外的轨道作用都归结为secondary。另外还可以考察VDW作用对endo选择性的贡献，说不定发现其实VDW贡献更大，然后颠覆一下textbook  。更完善的做法是从初始结构到过渡态的irc路径上每隔一段算一算vdw 和 secondary 的比例。但是这个做法的问题是不好控制primary作用原子的相对位置一样。 |
|
这种pmo的东西本身是为了解释实验而造出的模型 没有啥物理意义 想办法用计算去验证的思路...... |
sobereva 发表于 2017-1-31 00:10 谢谢Sob老师回答~~不过选取“IRC上更靠近TS的结构”没大明白，只能用“更靠近TS的结构”而不能用TS是因为那样会造成即将形成的共价键对分析的干扰么？ 另外对 A+B-->AB 反应，您后面说的计算是不是类似于：对复合物跑一个MP2，但这个MP2的轨道是 A 和 B 的 HF 轨道，且只考虑 A (occupied) --> B (virtual) 的激发、和 B (occupied) --> A (virtual) 的激发，而不考虑 A-->A, B-->B 的激发？" |
smutao 发表于 2017-1-30 22:18 我曾想过E2的事，不过 NBO 的 E2 中是否包含了“NBO 轨道很渣”造成的相互作用能？比如苯环的NBO E2实际体现了芳香性（优化后的分子轨道）与 Pi-NBO 轨道（与优化的轨道相差很远的形式）两种表示间的差异，而没有主要体现分子间的相互作用。 不知道能不能用真正的分子轨道做E2？大概意思是考虑，两个单体各干各的，相比一块干活差了多少。不过如果对过渡态做这个分析似乎又混合了即将形成的共价键的贡献。不知如何解决。 |
ggdh 发表于 2017-1-30 10:08 谢谢您指点，不过有些不太明白。第一条里“删除原子上的波函数”是什么含义？仅删除该原子的基函数？那样相当于电子的离域性突然大大加强了，分子间的轨道相互作用感觉会变得出奇巨大。或是什么其他含义？ 2 & 3 的问题都在如何选取一个“其他相互作用类似，只有 secondary interaction不同”的参考点。因为希望解释的反应是 secondary interaction（倾向endo） 与位阻（倾向exo）的竞争，所以按图上两种情况计算的△△G≠或其他量、都掺杂了除 secondary interaction 之外的其他成分，还需要被进一步拆分才能分析 4. 找的图确实画错了，应该有三个节面→_→ |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2025-8-16 09:49 , Processed in 0.160534 second(s), 25 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites