🔬 材料科学

Two-dimensional structures formed by pair potentials: An analytical and numerical study

本研究通过解析与数值方法证明，由于长程原子间相互作用的存在，薄膜过剩能量的 $1/N$ 标度律在单层极限下失效，从而证明仅利用对势即可形成非范德华力二维材料。

原作者： Shota Ono

发布于 2026-10-22

📖 1 分钟阅读☕ 轻松阅读

原作者： Shota Ono

原始论文采用 CC BY 4.0 许可（http://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写或认可的。如需技术准确性，请参阅原始论文。阅读完整免责声明

在材料科学领域，科学家们不断寻找创造极其薄的、二维物质薄层的方法。像石墨烯这样的某些材料，由自然堆叠在一起的层组成，这些层由微弱且温和的力量维系。这些被称为范德华材料，并且相对容易被剥离成单层薄片。然而，许多其他有用的材料，例如金属或某些半导体，并不具备这种天然的层状结构。它们的原子通过强烈的、方向性的键结合在一起，形成一个固体的三维块体。长期以来，人们一直不清楚是否可能从这些非层状结构中分离出单一且稳定的薄层。预测这种材料是否能作为独立薄层存在的关键，在于理解其能量随着厚度变薄而发生的变化。如果能量在移除层数的过程中表现出可预测的行为，那么该材料可能无法以单层形式存在。但如果能量在极薄极限处出现意料之外的偏移，则表明一个稳定的二维版本是有可能存在的。

来自日本室兰工业技术学院的一位研究人员通过剥离现实世界化学中的复杂细节，转而关注原子间吸引与排斥的基本规则，对这一问题进行了全新的审视。该研究并未使用通常模拟特定金属或化合物时所需的沉重且复杂的计算，而是采用了一种基于“对势”（pair potentials）的简化方法。这种方法假设原子之间的能量仅取决于任意两个原子之间的距离，忽略了当三个或更多原子参与时发生的更复杂的相互作用。通过构建一个由这些简单的相互作用原子组成的薄膜数学模型，研究人员能够通过解析方式，精确追踪薄膜的能量如何随着层数从厚块体减少到单层而发生变化。

研究首先建立了一个适用于大多数厚膜的基础规则。在这些较厚的堆叠体中，由于表面产生的额外能量遵循一种非常规律的模式：当层数翻倍时，每层的额外能量减半。这是一种标准的标度律，当原子间的力是短程的（即一个原子只能感受到其直接相邻原子的拉力或推力）时，该规律表现良好。在这种情况下，数学预测即使是单层也会遵循这一平滑趋势。然而，研究人员发现，当原子间的力是长程的，这一规则就会失效。在这些场景下，一个原子可以感受到数层之外的邻居的影响，而不仅仅是与其接触的邻居。在这种情况下，单层的能量并不符合厚膜所预测的平滑曲线；相反，它比预期值更低。

为了证实这种理论上的失效是真实的而非仅仅是数学上的巧合，研究人员使用两种著名的原子相互作用模型——Lennard-Jones势和Morse势——进行了计算机模拟。这些模型就像是描述原子行为的不同规则集，其参数可以通过调整来使原子间的力变为短程或长程。当模拟使用短程力时，结果完美符合标准规则，即使对于单层也是如此。但当参数被调整以产生长程力时，模拟显示出明显的偏差。单层的能量显著低于标准规则所预测的数值。这种偏差通过一个特定的数值来量化，该数值衡量了能量曲线在单层极限处的弯曲程度。对于短程力，该值接近 -1，表明完全符合标准规则。对于长程力，该值偏移到了如 -0.8 或 -0.47 之类的数字，标志着对预期模式的明显突破。

研究结果表明，从非层状块体结构中形成稳定的二维材料的能力，并不一定需要复杂的、多原子的化学键。相反，只要这些相互作用足够长程，就可以仅利用简单的两体相互作用来创造这些结构。这提供了一种更简单的新视角，用以思考为什么像铅烯（plumbene，单层铅）或金烯（goldene，单层金）这类材料可以存在，而其他材料却不能。该研究暗示，将原子结合在一起的“粘性”（即这些力的强度），特别是这种粘性延伸的距离，才是关键因素。如果力延伸得足够远，能够连接起跨越多个层级的原子，材料就能将自己稳定为单层薄片。虽然这项研究并未考虑环境效应（例如基底产生的应变或其它物质的存在），但它提供了基础性的见解。它表明，发现新型二维材料的路径，或许可以通过寻找那些原子相互作用能延伸至紧邻邻居之外的系统来实现，从而让单层的微妙平衡能够抵御来自块体的拉力。

技术摘要：由对势形成的二维结构

问题陈述
从非层状体相结构中衍生出的二维（2D）材料（即非范德华力（non-vdW）二维材料，例如硅烯、铅烯、金烯）的形成，在弥合理论预测与实验合成之间的差距方面仍面临挑战。预测此类材料稳定性的一个关键指标是薄膜相对于体相的过剩能量（ $\Delta E$ ）对厚度的依赖性。前人的工作已经确立，在单层极限下， $1/N$ 标度律（其中 $N$ 为层数）的失效是识别非范德华力二维材料的一个有用描述符。然而，尽管第一性原理计算可以准确预测这些能量，但支配这种标度律失效背后的物理机制在解析层面上尚未得到充分理解。具体而言，目前尚不清楚仅靠对势（pair potentials）是否足以解释在实际材料中观察到的偏离 $1/N$ 标度律的现象，或者是否必须引入多体相互作用。

方法论
本研究采用结合解析与数值的方法，研究由对势形成的 $N$ 层薄膜的过剩能量。

解析模型： 作者假设一个由单原子晶体组成的紧密堆积结构（沿表面法线方向堆叠的三角晶格），其原子间相互作用由对势 $v(R)$ 控制。 $N$ 层薄膜的总能量被推导为层内相互作用能与层间相互作用能之和。过剩能量 $\Delta E(N)$ 被定义为薄膜每层能量与体相能量之差。分析区分了短程与长程原子间相互作用，以确定它们对 $N=1$ 时标度行为的影响。
数值模拟： 为了验证解析理论，研究使用原子模拟环境（ASE）进行了数值模拟。测试了两种特定的对势：Lennard-Jones (LJ) 势和 Morse 势。通过改变 LJ 势的截断半径（ $R_c$ ）和 Morse 势的指数前因子（ $\rho_0$ ）来系统地调节这些相互作用的范围。研究考察了源自面心立方结构 (111) 面（ABC 堆叠）的 $N$ 层薄膜，且未进行几何弛豫，以隔离相互作用范围的影响。

主要贡献与结果

标度失效的解析推导： 本文从解析上证明，过剩能量 $1/N$ 标度律的有效性关键取决于原子间相互作用的范围。
- 对于短程相互作用，即只有相邻层发生显著相互作用时，过剩能量即使在单层极限（ $N=1$ ）下也遵循 $1/N$ 标度律。在这种情况下，外推值 $\Delta E_{ext}^1$ 与实际单层能量 $\Delta E(1)$ 相匹配。
- 对于长程相互作用，即相互作用延伸至最近邻以外时， $1/N$ 标度律会失效。解析推导表明，在 $N=1$ 时的实际过剩能量低于从 $1/N$ 趋势中外推的值。这导致了一个“向下偏差”，由参数 $p_1 > -1$ 来量化（其中 $p_1$ 是 $\ln \Delta E$ 对 $\ln N$ 在 $N=1$ 处的斜率）。
数值确认： 使用 LJ 和 Morse 势的模拟证实了解析预测。
- 当相互作用范围较短时（例如 LJ 势的 $R_c = 2\sigma$ ，Morse 势的 $\rho_0 = 8$ ），系统遵循 $1/N$ 标度律（ $p_1 \approx -1$ ）。
- 当相互作用范围增加时（例如 LJ 势的 $R_c = 10\sigma$ ，Morse 势的 $\rho_0 = 3$ ），观察到了明显的向下偏差（分别为 $p_1 = -0.8$ 和 $-0.47$）。
- 研究进一步指出，这种偏差即使在没有结构弛豫的情况下也会发生，表明几何弛豫对标度律本身的影响较小。

意义与主张
这项工作的核心意义在于确立了非范德华力二维结构在理论上可以仅由对势形成，前提是原子间相互作用具有足够的长程性。 $1/N$ 标度律在单层极限下的失效并不完全是复杂多体效应（如共价键或化学键）的结果，而是可以从对相互作用的范围这一基本属性中产生的。

作者得出结论，层间相互作用的范围是调节过剩能量标度的关键因素。虽然本研究并不声称能取代针对特定材料预测的第一性原理方法，但它提供了一个能够捕捉标度律失效背后本质物理过程的极简解析模型。研究表明，在这一简化框架内，长程对势是形成稳定的非范德华力二维材料的必要条件，这为理解为什么某些非层状体相材料可以剥离成稳定的二维单层提供了理论基础。作者明确指出，超越对势（多体相互作用）的影响以及环境因素（衬底、应变）仍是未来使用密度泛函理论进行研究的课题。

技术摘要：由对势形成的二维结构

类似论文