计算化学公社

标题: NaCl晶体摻杂+2的Rh计算epr的g因子偏小，为什么呢？ [打印本页]

作者
Author: wudazhuang509 时间: 2016-12-12 15:32
标题: NaCl晶体摻杂+2的Rh计算epr的g因子偏小，为什么呢？
NaCl晶体摻杂+2的Rh后，计算epr的g因子的结果总是偏小，虽然优化后的结构和文献给的结构很接近。各位大神们，能不能帮我看下我的输入文件有什么问题么？

作者
Author: wudazhuang509 时间: 2016-12-12 15:34
忘记加原始文献了，见附件

作者
Author: winterzen 时间: 2016-12-12 19:03
你为什么要在周期体系中加入一个正电荷，在我的印象里带电荷的周期体系是不能够直接进行DFT计算分析性质的，因为带电周期体系中每个image之中都有一个电荷，这些电荷间会产生一个非物理真实的静电相互作用。如果非要计算带电周期体系的话需要一些方法进行校正，比如这篇论文：Phys. Rev. Lett. 110, 095505。我并不了解EPR的计算，但感觉误差可能大部分来源于这里。

作者
Author: wudazhuang509 时间: 2016-12-12 21:41

winterzen 发表于 2016-12-12 19:03
你为什么要在周期体系中加入一个正电荷，在我的印象里带电荷的周期体系是不能够直接进行DFT计算分析性质的 ...

谢谢你的回复，因为我希望用一个带+2电荷的Rh替换一个带+1电荷的Na，这样整个体系应该呈现的是带+1个电荷。我是新手，不是很了解带电荷周期性体系不能直接进行DFT计。之前是让整个体系呈电中性，得到的结果也是偏小。而且不知道怎样体现替换的Rh是带+2电荷，你能帮我解答下不？

作者
Author: winterzen 时间: 2016-12-12 22:57
你所考察的问题的物理真实是什么样的？是一个体相NaCl晶体中掺杂有极少量的Rh 2+吗？真实材料中的残余电荷是如何被平衡的？
如果你使用了带电的周期体系计算，那么你的模型中就引入了掺杂浓度以及电荷浓度的概念，例如，你用NaCl单元晶胞掺杂1个Rh和用2x2x2的NaCl超胞去掺杂1个Rh是两种不同的Rh掺杂浓度，带上电荷的话也会给体系带来不同的电荷浓度，同时造成了不同的非物理真实的静电相互作用。你所考察的相关性质也应该会依赖于这个浓度变化。
如果你想得到的是无限稀释情况下的性质（如果你的物理真实就是无限稀释情况，如晶体带电点缺陷的形成能），你可以用两种方法：
1. 使用从小到大不同规模的超胞掺杂1个Rh计算性质，然后作出性质随超胞规模的变化曲线进行外推。
2. 使用那篇PRL论文里的方法或其他类似方法进行校正。

或者你还可以使用足够大的非周期体系进行计算，这样就不会有浓度概念在里面干扰

作者
Author: wudazhuang509 时间: 2016-12-13 13:33

winterzen 发表于 2016-12-12 22:57
你所考察的问题的物理真实是什么样的？是一个体相NaCl晶体中掺杂有极少量的Rh 2+吗？真实材料中的残余电荷 ...

谢谢您的回答，希望获得的是无限稀释情况下的系统的局域结构性质，计算资源有限，我先看下PRL文章，如果有问题在继续向您请教。谢谢

作者
Author: sigma-jlu 时间: 2023-8-3 11:49

winterzen 发表于 2016-12-12 22:57
你所考察的问题的物理真实是什么样的？是一个体相NaCl晶体中掺杂有极少量的Rh 2+吗？真实材料中的残余电荷 ...

winterzen提到的电荷问题、低浓度外推以及那篇prl中添加势函数补偿的方法，都很有启发性，谢谢。

关于楼主提到的inp文件，是不是还有基组选择的不妥？个人经验，DZVP-MOLOPT-SR-GTH用于精细计算之前的geo_opt还凑合，计算性能，尤其是内层电子结构敏感的性能，比如fermi contact，可能不适合。

Kaupp组推荐过def2-TZVP基组，试了也不容易收敛，不如Jensen开发的pcseg-2基组。

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