红红的太阳 发表于 2025-11-13 09:22
1.俺寻思哈，这个两性离子吸附到石墨烯表面的构型应该不是你这样的，即使你的石墨烯片带正电荷，也有概率和 ...

感谢您的回答~
1、我的想法是初始模型可以先这样构建，然后在跑动力学的过程中，I3-逐步吸附到N+上，然后醚链卷曲，就像您说的，醚和石墨烯产生了相互作用。
2、石墨烯是正极材料，为了尽可能还原实际电池体系，所以考虑是需要加石墨烯的。
3、我会参考您关于盒子调整的建议，再次感谢！但是如果体系中有石墨烯的话，就不能再做NPT了吧，因为需要保持石墨烯的周期性，盒子大小不能改变，唯一能改变的只有Z方向的长度。
4、我参考的场强来源于http://bbs.keinsci.com/thread-22552-1-1.html.由于对所加场强没有太明确的概念，所以本着试一试的想法选取的0.5V/nm场强。您提供的这两篇文献我会认真阅读并联系作者尝试要到源代码，非常感谢您提供的帮助

1.俺寻思哈，这个两性离子吸附到石墨烯表面的构型应该不是你这样的，即使你的石墨烯片带正电荷，也有概率和醚那边产生相互作用吧？要不要先拿石墨烯和两性离子跑一个初始构型出来？
2.有必要加入石墨烯吗？如果只是看两性离子和碘以及I3-阴离子的相互作用，直接拿这几个物质做盒子不就好了..没必要实验里面有什么就加什么吧
3.溶剂跑到上层有可能是盒子太大，分子太少，调整一下盒子到你电解液密度1.2倍左右的体积，然后加分子做NPT NVT应该就可以了
4. 俺寻思哈，你已经用平板电极去近似多孔电极了，再按照实验的电流密度换算模型的面电流密度，不靠谱..
别处参考的场强来源建议给一下。另外，如果目的是想知道不同电极电势下电极表面分子的吸附情况，我倒是见过通过更改石墨原子电荷来近似的方法(constant potential method (CPM))：10.1021/nl4031083；10.1021/ja104273r；不过呢，有文献报道认为使用上述方法来推断多孔电极的双电层结构会存在误差较大的问题，见10.1038/s41467-022-29761-z所以现在比较"流行"constant charge method (CCM) 方法，但需要写代码自己实现。或者可能的话，可以尝试联系联系文章作者看看能否要到源码