三元体系的共掺杂，掺杂形成能的计算该如何处理

灰飞的旋律

请问下各位做过掺杂的大佬像这种三元体系的共掺杂，掺杂形成能的计算该如何处理呢，比如AXY体系，A是金属组分，X Y 都是非金属组分，然后N P共掺杂分别替换X Y。我查了下文献，下面这篇NaNbO3共掺杂N和Mo的，他俩分别取代O和Nb。
文中的计算过程如下

我现在是有个疑问，他在富氧条件下， ΔμO=0时，ΔμNb的化学势如何根据公式 5 6 7 8 确定范围的，主要是这个时候Na的化学势如何选取的？？？
还有就是，他讨论的是富氧和富Nb条件下的化学势变化，那我要是替换AXY中X和Y元素，那应该就讨论富X和富Y的条件下化学势变化范围就好是吧，如果是的话，问题同样就是这两种情况下金属元素A的化学势分别怎么选择？？？

灰飞的旋律

上面的问题想明白了，新的问题是像这张图里这这种图形可以用origin做出来吗，这种图叫什么名字呢，没有搜到相关教程？？？
图4

它是定义了x轴和y轴的区间范围，然后填色是根据一个特定的函数公式来填的。很明显这个公式，也就是掺杂形成能的公式是一个二元一次公式。
这种图在origin里怎么填色画图呢？？

灰飞的旋律

灰飞的旋律 发表于 2019-5-22 21:32
上面的问题想明白了，新的问题是像这张图里这这种图形可以用origin做出来吗，这种图叫什么名字呢，没有搜到 ...

搞定，需要设置边界条件

granvia

个人感觉这些利用边界条件来估算化学势的方法非常粗糙，我曾在论坛提到过一篇JPCL的文章，国内搞物理的算的，他们就（很武断地）取取值范围窗口的中心点为某物质的化学势，非常无语。很可能编辑和审稿人都不是量化专业的吧，反正现在第一性原理模拟给我的印象就是，没有不能算的东西，至于精度和可靠性吗，就不用太较真了，跟实验符合就ok了。个人看法，勿喷

zmlly

楼主，能不能分享一下这篇文献？