请问用MS中Forcite模块计算缓蚀剂分子在Fe(110)的吸附能出现不收敛，该怎么解决

暴躁兔子Zz_

请问用MS中Forcite模块计算缓蚀剂分子在Fe(110)的吸附出现不收敛，该怎么解决。
在MS软件中引入铁晶胞，晶格参数为a=b=c=2.8664 Å，α=β=γ=90°，沿Fe晶胞（110）表面切割，切割厚度为21.20 Å，构建了一个11×11×11的超晶胞，超晶胞的长度和宽度为27.30 Å，选择Amorphous Cell模块建立缓蚀剂溶液，缓蚀剂溶液中放入1000个水分子，5个HCl分子，HCl分子的存在形式为H3O+和Cl-。

用Build Layers 构建最后的模型由铁晶胞和缓蚀剂溶液层组成，模型顶部留有40 Å真空层
      

选择几何优化模块，选择Smart算法优化。收敛准则设定为2×10-5 kcal/mol，最大优化步长为20000步。固定铁原子层。用基于原子的方法和Ewald方法计算了范德华作用和静电作用。截断距离为18.5 Å

   


采用正则系综（NVT）进行模拟，温度为298K，温度由Andersen恒温器控制。总模拟时间为5000 ps，模拟积分步长为1 fs。每10000 fs为一帧。用速度方程计算牛顿运动方程。分别用基于原子的方法和Ewald方法计算了范德华作用和静电作用，缓蚀剂分子没有在Fe(110)表面进行吸附。



吸附能

lixiaoguai

您这个问题解决了吗，我想要建模的跟您的类似，我想要让fe与上面的溶液进行一个化学反应，您有什么思路提供的吗

暴躁兔子Zz_

lixiaoguai 发表于 2024-12-25 21:09
您这个问题解决了吗，我想要建模的跟您的类似，我想要让fe与上面的溶液进行一个化学反应，您有什么思路提供 ...

最后解决了，边界条件不收敛

暴躁兔子Zz_

lixiaoguai 发表于 2024-12-25 21:09
您这个问题解决了吗，我想要建模的跟您的类似，我想要让fe与上面的溶液进行一个化学反应，您有什么思路提供 ...

截断距离为12.5 Å，就解决了后续问题