reaxff模拟阳极腐蚀问题求助

321321321

各位老师同学好
我需要做一些阳极腐蚀的研究。通过文献查阅发现reaxff方法再加入外电场很适合我的研究，并且找到了相关力场文件。
但是有一个问题，阳极是要接入电源正极做失电子氧化反应的，就是阳极区域整体是失电子的。而仅仅加入外电场的话求解域电荷肯定是平衡的。

我的研究目的就是对阳极腐蚀过程的趋势有一个大致的定性，能做出来金属表面电化学反应的效果就ok，以此作为实验结果分析的一部分。
目前想到的解决办法有两种，但不知道能不能行得通
1.周期性的直接设置最左面金属原子的电荷，使它们带正电。reaxff力场应该可以把这些正电荷向右传递到反应面。
2.建立下图这样的模型，边界条件ppp，也就是说上下金属板实际上是连起来的。在金属的区域设置很大的电场迫使其中的电荷转移，使一边带正电一边带负电。给溶液区域小电场来模拟电迁移。

这两种解决方法都基于一个假设，就是reaxff方法可以实现金属原子之间的电子转移，也就是导电。

在这里请教各位老师，这种办法可不可以实现，或者有没有可以做出阳极表面溶解效果的方法呀。之前听说过dft、量子化学、第一性原理这些概念，暂时还没有仔细研究它们。如果有方法能实现的话还请各位老师指个方向我去慢慢研究。
谢谢各位了

GZU_220

同学，搞得咋样了，我今年开学导师也让我做这个，也是reaxff方法再加入外电场。
学机械的，过来搞电化学加工，化学啥的之前也没学过，现在也是一头雾水，不知道咋搞

baoluyao

无法实现，ReaxFF中所有原子的电荷根据反应的进行而动态调整的。另外在LAMMPS中加电场的话动态电荷平衡要用acks2，但是这个方法力场参数很少

321321321

GZU_220 发表于 2025-3-7 14:55
同学，搞得咋样了，我今年开学导师也让我做这个，也是reaxff方法再加入外电场。
学机械的，过来搞电化学加 ...

一模一样，分子动力学研究电化学加工，都找不到几篇这方面的文献。现在还在调研，打算实在不行就还是用有限元吧。做个流场加点实验就基本能毕业了。

GZU_220

Reactive molecular dynamics study on carbon steel corrosion induced by chloride: Effects of applied potential and temperature，我目前还在打算先实现这篇文献的效果，其实我只要能搞出像这图中的氧化膜就行了。

chencp

GZU_220 发表于 2025-3-20 16:32
Reactive molecular dynamics study on carbon steel corrosion induced by chloride: Effects of applied  ...

您好，我也在做这个，可以交流一下嘛

chencp

GZU_220 发表于 2025-3-20 16:32
Reactive molecular dynamics study on carbon steel corrosion induced by chloride: Effects of applied  ...

我好像私聊不了你

Graphite

我本来打了大概300个字，考虑到说太多意义不大，还是删了，简而言之目前没有真正理论上又合理、工程实现上又靠谱的办法。

要做出5楼那样的图并不需要改电荷什么，只要调参反复试就行了。

321321321

GZU_220 发表于 2025-3-20 16:32
Reactive molecular dynamics study on carbon steel corrosion induced by chloride: Effects of applied  ...

我也看了一些文献，发现这种原子尺度的模拟方法没办法做出真正的电极，近几年发表的级别很高的相关论文也只是用一些不直观的方法近似出电极的效果。用reaxff方法模拟出来的效果实际上是的金属表面氧化过程而没有电极电子转移的效果，对应于现实就是单纯把材料丢到电解液里让它自然腐蚀，这样得到的结论很难用来解释实验结果。我还看到有论文用阴离子侵入金属的深度来表示加工速度，这样倒是很直观但是感觉原理上不太靠谱