求助拟合二面角QM-MM数值比文献中高该如何解决

zwli · 发表于 Post on 2025-8-24 14:34:51

最近在用gromacs模拟聚合物的时候，发现聚合物的HOMO能级实验CV测量的高很多，不知道是什么原因。想请教一下各位老师。
我在一篇文献中发现易老师组有做相关聚合物能级的研究。
标题为:Toward Low Energetic Disorder in Organic Solar Cells: The Critical Role of Polymer Donors
https://doi.org/10.1021/acs.jpclett.5c00132
因此我想学习他们的方法，从二面角拟合开始方法验证。针对图中的结构，运用柔性扫描还原了DFT（黑色曲线），文献中的黄色曲线是QM-MM的差值。
因此进一步将0到180°的结构，运用cat进行了拼接，执行md.mdp动力学模拟，以及gmx_mpi energy，获取了这些构型的单点能量。但是，得到的图像结果会比图中数值高得多。
在itp中，我们将转动的二面角力常数设为了0，其他键长参数是采用b3lyp/6-311g(d,p)优化，通过sobtop生成，用def2TZVP基组获取RESP电荷。想请教一下各位老师，这种情况是因为mdp有问题，还是因为非键参数有问题呢？或者是别的其他什么原因？
该如何解决呢？

sobereva · 发表于 Post on 2025-8-28 03:40:33

zwli 发表于 2025-8-27 13:15
感谢社长的回复，我一开始的理解也是觉得这些体系可能在HOMO能级上，并不能很好的和实验的CV吻合。但是在 ...

拿HOMO能级负值去近似CV得到的溶剂下的IP的做法本来就是不当的，在这个不当的前提下做什么其它事都没有意义，我也不想做任何评论

zwli · 发表于 Post on 2025-8-27 13:15:24

sobereva 发表于 2025-8-25 08:19
仔细把此文看了
正确地认识分子的能隙(gap)、HOMO和LUMO
http://sobereva.com/543（http://bbs.keinsci.c ...

感谢社长的回复，我一开始的理解也是觉得这些体系可能在HOMO能级上，并不能很好的和实验的CV吻合。但是在易老师的文章中，以PM6为例，我发现他们将IP等同于-HOMO后，能级与实验的CV还比较吻合。所以我学习他们的二面角拟合方式进行了拟合，但是在这个拟合上与他们的结果有些差距，您觉得我的拟合过程是否存在问题呢？

zwli · 发表于 Post on 2025-8-25 19:47:28

本帖最后由 zwli 于 2025-8-27 13:15 编辑

感谢社长的回复。

sobereva · 发表于 Post on 2025-8-25 08:19:10

仔细把此文看了
正确地认识分子的能隙(gap)、HOMO和LUMO
http://sobereva.com/543（http://bbs.keinsci.com/thread-16758-1-1.html）

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

[GROMACS] 求助拟合二面角QM-MM数值比文献中高该如何解决

浏览过的版块