计算化学公社

标题: 请问用MS中Forcite模块计算缓蚀剂分子在Fe(110)的吸附能出现不收敛，该怎么解决 [打印本页]

作者
Author: 暴躁兔子Zz_ 时间: 2023-10-19 10:34
标题: 请问用MS中Forcite模块计算缓蚀剂分子在Fe(110)的吸附能出现不收敛，该怎么解决
请问用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 Å

(, 下载次数 Times of downloads: 23)
(, 下载次数 Times of downloads: 17)

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

吸附能
(, 下载次数 Times of downloads: 18)

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

作者
Author: 暴躁兔子Zz_ 时间: 2025-5-12 17:03

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

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

作者
Author: 暴躁兔子Zz_ 时间: 2025-5-12 17:04

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

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)