掺杂模型里，一个原子替位常被当成“低浓度掺杂”。但在周期性边界下，它其实代表一个无限重复的有序掺杂阵列。超胞越小，掺杂原子之间的周期性距离越近；超胞越大，才更接近稀释掺杂。掺杂浓度不是输入参数里的小细节，而是模型回答问题的边界。
一个掺杂原子、两个掺杂原子和高浓度模型回答的是三类问题。单掺杂看孤立杂质的局域效应；双掺杂看杂质之间相互吸引、排斥或配对；高浓度模型则接近合金化、固溶体或有序相。把它们放在一起只按形成能最低排序，很容易把浓度效应误写成位点效应。
一个掺杂原子也有周期性邻居
若在 2 × 2 表面超胞里替换一个原子，浓度可能已经很高；若在 4 × 4 超胞里替换一个原子，周期性镜像距离增大，局域畸变和电子态更接近孤立杂质。两者不能只叫“一个掺杂原子”，还要写清超胞大小和被替换位点总数。
c = Ndopant / Nhost-site 式(1)
对于表面 slab，还要说明 c 是按表层位点、全部可替换位点，还是整个超胞原子数来算。若不说明分母，同一个“6.25%”可能代表完全不同的物理模型。掺杂浓度必须和超胞、层位、替换位点一起报告。
图1 ZnO 中双 Li 掺杂不同构型的形成能，展示近邻与远离构型会改变能量。DOI：10.1039/c8ra04811h
低浓度和高浓度关注点不同
低浓度模型主要问：一个杂质会不会稳定存在，会不会引入局域态，会不会改变附近吸附位。高浓度模型主要问：多个杂质是否形成有序分布，晶格是否整体畸变，材料电子结构是否整体改写。低浓度强调局域，浓度升高后则要关注杂质-杂质相互作用。
如果一个低浓度模型显示吸附能改善，不能直接推出高浓度也更好。高浓度下杂质可能相互抵消电荷调节，也可能形成新相或聚集。此时应把“单杂质效应”和“浓度演化”分开写。
近邻、远邻和成对构型要分组
两个掺杂原子可以相邻，也可以隔得很远；可以在同一层，也可以在表层和次表层；可以围住一个吸附位，也可以分别影响不同区域。它们的能量差反映的是杂质之间相互作用，而不是单个掺杂效应的简单加倍。
Epair = E2D − 2E1D + Ehost 式(2)
若 Epair pair > 0，更倾向分散。这个量能帮助判断高浓度模型是否会形成局域富集。双掺杂要先看相互作用，再谈吸附活性。
图2 Li 掺杂 ZnO 的相稳定区域，说明化学势和构型会共同影响掺杂稳定性。DOI：10.1039/c8ra04811h
同浓度也可能有很多局域环境
同样是 12.5% 或 25% 掺杂，杂质分布可以均匀、团聚或平面富集。局域环境不同，吸附位点附近的配位、键长、应变和电子态都会不同。若只取一个构型，就可能把构型选择当成浓度趋势。
比较时建议同一浓度下先枚举多个分布，再选低能构型和代表性高能构型。低能构型回答稳定分布，高能构型帮助判断趋势敏感性。
图3 12.5% Y 掺杂体系中的均匀、团聚和平面分布构型，展示同浓度下局域环境仍可不同。DOI：10.1021/acs.jpcc.0c04594
从缺陷模型变成有序相模型
当掺杂浓度很高时，体系不再只是一个缺陷扰动，而可能接近合金、固溶体或有序化合物。此时形成能、吸附能和电子结构都受到整体晶格变化影响。若仍按“一个杂质改变一个活性位”的语言解释，就会低估浓度带来的结构重排。
高浓度还会改变应变状态。多个掺杂原子可能共同拉伸或压缩 M-O、M-S、M-N 键，使局域畸变从单点扰动变成全局晶格响应。此时要看晶格参数、键长分布和态密度整体变化，而不是只盯着一个吸附位。
图4 25% Y 掺杂下不同 Y-Y 构型的相对能量，说明高浓度需要考虑杂质排布。DOI：10.1021/acs.jpcc.0c04594
形成能要按同一化学势口径比较
掺杂形成能常包含宿主原子和掺杂原子的化学势。不同浓度比较时，化学势边界要一致；若一个模型用富金属条件，另一个用富氧条件，能量差就不是浓度差。
Eform = Edoped − Ehost − nμD + nμH 式(3)
这里 D 是掺杂元素，H 是被替换的宿主元素。若替换多个原子，n 要对应实际替换数。多个掺杂原子的形成能还可以除以 n 得到平均形成能，但平均值会掩盖配对效应，因此最好同时报告总形成能、平均形成能和配对能。
图5 不同 Y 掺杂浓度下金属-氧和金属-金属距离变化，展示浓度会改变局域结构。DOI：10.1021/acs.jpcc.0c04594
先做浓度阶梯，再做位点枚举
建议先选一套浓度阶梯：低浓度单掺杂、中等浓度双掺杂、高浓度多掺杂。每个浓度内部再枚举不同层位和相对距离。这样可以先看浓度趋势，再看局域构型差异。若顺序反过来，很容易在大量构型里看不清主变量。
对催化吸附计算，最好只把同一浓度、同一超胞大小、同一覆盖度下的吸附能放在一起比较。跨浓度比较可以讨论趋势，但不能把最低吸附能直接当成最优掺杂浓度，因为高浓度模型也可能更难形成或更容易相分离。
图6 Li 掺杂 ZnO 中不同缺陷构型之间的转化能垒，提示浓度和构型稳定性还可能涉及迁移路径。DOI：10.1039/c8ra04811h
结果表达要把三类模型分开
单掺杂模型写“孤立杂质效应”，双掺杂模型写“杂质相互作用”，高浓度模型写“集体结构调节”。如果某个结论只在高浓度下成立，就不要反推到稀释掺杂；如果某个吸附位只在双掺杂成对时出现，也不要写成单个掺杂原子的普遍作用。
一个清楚的结论可以这样组织：低浓度 D 改变邻近金属位吸附；双掺杂显示 D-D 近邻更低能或更高能；高浓度下晶格整体畸变使吸附自由能进一步移动。这样能把浓度、分布和催化活性之间的关系拆开。
对二维材料或表面 slab，还要额外注意真空方向和表面覆盖度。一个掺杂原子如果只替换表层位点，它的浓度不能和体相全部位点的替换浓度混写；如果上下表面都暴露，单面掺杂和双面掺杂也不是同一类模型。对于吸附计算，掺杂浓度和吸附覆盖度会同时改变局部环境，因此最好固定其中一个变量，再考察另一个变量。
若目标是筛选催化剂，可以先用较低浓度判断杂质是否有局域调节能力，再用双掺杂或高浓度模型检查这种调节是否会被杂质相互作用放大、抵消或改写。这样得到的结论会比直接挑最低吸附能更清楚。
如果不同浓度下最稳定构型来自不同层位，也不要急着说浓度越高越稳定。表层掺杂、次表层掺杂和体相掺杂面对的化学环境不同，形成能变化可能来自层位差异，而不是浓度本身。比较浓度趋势时，最好先固定层位；比较层位趋势时，再固定浓度。
最终写结论时，可以把浓度作用拆成三句话：低浓度提供局域电子调节，双掺杂体现杂质相互作用，高浓度带来整体晶格和能带改变。这样一来，即使三个模型都显示吸附能改善，也能说明改善来源是否相同。
要点一：一个掺杂原子也代表周期性阵列，必须写清超胞和位点分母。
要点二：双掺杂不是单掺杂翻倍，要看近邻、远邻、层位和配对能。
要点三：高浓度模型可能从缺陷问题变成有序相或固溶体问题。
要点四：形成能、平均形成能和配对能回答不同问题，不能只报一个最低值。
要点五：吸附能比较要绑定浓度、超胞、局域构型和覆盖度。
