理论计算新手一个,还望各位大佬不惜赐教。如题,我在复现一篇文献的时候,文献中提到了使用72个显式水分子、反应中间体和表面原子进行5ps的AIMD以获得合适的结构,再去掉顶层的36个显式水分子并由隐式溶剂模型代替。而在其Support Information中提到的是使用Ag原子、48个水分子来准备其反应物。(而且这文章里面正文说slab是3x3x4,Support Information中说是3x4x3....)想请教一下各位师兄师姐和老师,他这里的带有显式水分子的隐式溶剂模型中到底有48个还是36个水分子,AIMD跑完以后应该删去几个?如果我用这样的模型去算吸附能、自由能的话,体系里面这么多水分子,计算量是不是会特别大?我是不是应该对每一个反应物比如CO2,H2O,*H,**COOH都进行一下AIMD然后再加入隐式溶剂然后再对其进行计算呢?
还有一个问题时关于M06-L泛函的,这个泛函我在用于气相CO结构优化时不收敛,这是为什么,是因为它不可以用于结构优化嘛?
Singh, M. R., Goodpaster, J. D., Weber, A. Z., Head-Gordon, M., & Bell, A. T. (2017). Mechanistic insights into electrochemical reduction of CO2over Ag using density functional theory and transport models. Proceedings of the National Academy of Sciences, 114(42), E8812–E8821.
doi:10.1073/pnas.1713164114