ldatea 发表于 2020-6-30 08:37 嗯嗯好的，十分感谢您的解答 |
j5888xm 发表于 2020-6-29 23:35 1.α-H的氢键比较弱的，但具体多弱我不清楚，你找个模型算一下就知道了。 2.一般来说，问题不大，构象平均以后降低的自由能不多。前提是你必须找准最低能量的构象，这个找不准误差就很大了。 |
ldatea 发表于 2020-6-29 20:33 谢谢解答，构象的问题我会多加注意，位阻作用可能确实对苯环共轭有不小影响。 如果用醛或醇作溶剂，会不会和阿尔法氢有相互作用而促进离解，还是说只有氢键并且主要还是在氧上。 还有，图省事都采用能量最低的构象，应该不会有什么不妥吧（手动狗头）。 哈哈，再次谢谢指点 |
本帖最后由 ldatea 于 2020-6-29 21:17 编辑 j5888xm 发表于 2020-6-29 20:17 一般来说，确实会回到同一个结构。但这个体系有特殊性，这两个苯环根本不可能共平面，至少有一个苯环要扭转。这种情况下阴离子可能不只一种构象。 ps:你这个体系确实一看就是不适合做反应的（除非你要在氧上反应，位阻的影响就会小很多）。 你要算pKa,HA=H++A-，那么这里面三个物种的能量一个都不能少，溶剂化氢离子能量不用算，一般有数据可查。其他两个都得自己算。 如果只关心烯醇负离子哪个更多，反应物的结构和能量其实根本不用算，因为你关心的是两种烯醇负离子的比例，这只与烯醇负离子的自由能有关。算pKa的思路是对的，但是没有必要。因为你只需要关心两种电离方式的pKa之比就可以了，这个pKa之比完全由两种烯醇负离子的自由能差决定。 但烯醇负离子可能存在构象问题，最简单的方法是分别找两种烯醇负离子中能量最低的构象，能量低的比例就更高。要算的再准一点就是考虑构象平均。 你算异丙基一侧烯醇负离子的能量，那就必须考虑苯环一侧哪种构象能量更低。 算苯环这一侧烯醇负离子的能量，也要考虑苯环这一侧由于有一个苯环不能共平面，是否有多种构象。 |
ldatea 发表于 2020-6-29 16:45 好的，谢谢纠正，十分感谢。 其实我拿这个分子举例，主要还是想考察双苯基共轭对烯醇化取向的影响，氧和氢异侧时能量较低，是不是简单拿这种构象的能量和另一侧的直接对比有失偏颇，因为这个构象并不适合外界夺氢。或者说我试着计算一下两种构象的pKa值看看有无明显差别。但是两种构象最后的阴离子应该是一样的，直接把两侧的阴离子能量对比是不是就能简明判断出结果。 再次感谢您的解答，辛苦了。 |
本帖最后由 ldatea 于 2020-6-29 19:28 编辑 j5888xm 发表于 2020-6-29 12:30 你说的二聚体，是指一个溶质分子+一个溶剂分子 吗？这种一般不叫二聚体，二聚体通常指同一种分子形成的复合物。 如果你还是用DMSO溶剂，隐式溶剂模型就可以了，DMSO是典型偶极溶剂，偶极矩较大，且正电中心在硫上，藏在内部，负电中心露在外面（氧），只考虑静电作用的时候，阳离子在DMSO中的溶剂化作用强，阴离子溶剂化作用弱。现在你的溶质也可以近似看成一个正电荷中心藏在内部，负电荷中心露在外面的偶极子。 如果用能和酮形成氢键的溶剂(比如水，位阻小的醇，胺，氯仿等），最最最少得带一个溶剂（体现氢键），然后对这个复合物用隐式溶剂模型。具体细节请看sob老师的博文(《谈谈隐式溶剂模型下溶解自由能和体系自由能的计算》)http://sobereva.com/327 我其实现在都不敢确定你这个体系到底会有几个构象，目前只是优化出两个极小值点而已。 我这边算出来的结果： 6-31g*(D3BJ)几何优化，然后用ORCA的wB97M-V/ma-def2-TZVP（RIJCOSX加速） 计算单点（这里用弥散函数的目的是为了和之后要算的阴离子体系用同一个基组，本身中性分子加弥散函数是不划算的，还有一个次要目的是估算后面算阴离子的时间，就是看wB97M-V/ma-def2-TZVP能不能算得动） O-C-C-H二面角168.106的构象，wB97M-V/ma-def2-TZVP单点-733.841302485153Hartree,M05-2X/6-31g*气相单点-733.793625308Hartree，SMD溶剂模型下M05-2X/6-31g*单点-733.812249149Hartree O-C-C-H二面角8.964的构象，wB97M-V/ma-def2-TZVP单点-733.835717944143Hartree,M05-2X/6-31g*气相单点-733.787252648Hartree，SMD溶剂模型下M05-2X/6-31g*单点-733.803927949Hartree第三个减第二个就是溶剂化自由能(状态是默认的298.15K，1atm,溶液标态1mol\L)。 |
sobereva 发表于 2020-6-29 00:00 十分感谢sob老师解答，我试着重新计算一下电荷 |
ldatea 发表于 2020-6-28 23:27 十分感谢，这个构象问题之前并没有注意到。 是不是需要和溶剂二聚体继续优化才能得出结论。 |
Accelerator 发表于 2020-6-28 17:27 谢谢Accelerator老师。 |
本帖最后由 ldatea 于 2020-6-29 09:41 编辑 sobereva 发表于 2020-6-29 00:00 谢谢sob老师。 |
ldatea 发表于 2020-6-28 15:45 B3LYP对这种情况完全没问题 当前体系用CHELPG不是什么好的选择，因为当前体系较为柔性，结果受构象影响很大，而且有些原子部分被包埋，拟合静电势电荷质量低。用Multiwfn算ADCH、Hirshfeld-I等电荷更有意义。 |
| 参与人数Participants 1 | eV +2 | 收起 理由Reason |
|---|---|---|
|
| + 2 |
本帖最后由 ldatea 于 2020-6-29 09:55 编辑 j5888xm 发表于 2020-6-28 22:31 我刚算了一下你的这个分子，发现气相下，苯环这一侧的CO-CHPh2单键旋转有两个构象，B3LYP/6-31g*（D3BJ）优化，O和H在C-C异侧的时候单点能更低。 但计算级别比较粗糙，没振动分析也没溶剂化，不好说谁更优势。 没有位阻和分子内弱相互作用的时候，酮的优势构象是O和α-H在同一侧的。说实话我也没料到这个体系这么坑爹，可能有两个能量差不多的构象，这两种构象的去质子会有明显差别，O和α-H的O-C-C-H在90度附近有利于较快地去质子（具体哪个角度最有利，跟酮和碱结构都有关，没法讲）。 |
| 参与人数Participants 1 | eV +2 | 收起 理由Reason |
|---|---|---|
|
| + 2 |
本帖最后由 j5888xm 于 2020-6-29 19:11 编辑 Accelerator 发表于 2020-6-28 17:27 谢谢解答，是不是说不考虑位阻时，苄位的酸性还是要强一些，但是考虑位阻以后就不一定了 |
ldatea 发表于 2020-6-28 15:45 谢谢解答，关键词部分以后会注意 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-2-20 04:26 , Processed in 0.227593 second(s), 26 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites