Irina V. Lebedeva 与 Andrey M. Popov 在《The Journal of Physical Chemistry C》期刊上于2020年发表了一项理论研究，题为“Energetics and Structure of Domain Wall Networks in Minimally Twisted Bilayer Graphene Under Strain”。该研究隶属于凝聚态物理和材料科学领域，特别关注二维材料（如石墨烯）的结构与物性调控。研究背景在于，近年来实验上在微小扭转角度的双层石墨烯中观察到了由畴壁（domain walls）分隔开的三角形共度（commensurate）畴网络，这些畴壁网络对体系的电子、磁学和光学性质产生了显著影响，例如产生了拓扑保护的螺旋态，为“能谷电子学”（valleytronics）提供了新平台。然而，尽管其电子性质有趣且应用潜力巨大，对此类系统中畴壁网络的能量学和结构本身的基础研究却相对匮乏。先前的大多数理论研究集中于孤立的、不相交的畴壁，而对于在应变下形成的周期性网络，特别是当共度畴尺寸远大于畴壁宽度时，缺乏有效的理论模型进行描述。因此，本研究旨在填补这一空白，通过发展理论模型，系统研究在底部石墨烯层同时存在双轴拉伸和微小扭转的情况下，双层石墨烯中三角形畴壁网络的参数（如畴壁取向、网络周期等）如何随应变变化，并探索利用外部应变调控该网络从而改变材料电子性质的可能性。

本研究是一项纯粹的理论与计算研究，其工作流程主要包含以下几个相互关联的步骤：理论模型构建、关键参数确定、能量表达式推导、以及针对不同物理情景的数值分析与讨论。研究未涉及具体的实验对象或样本，而是以双层石墨烯这一体系为研究对象，通过物理建模和数学推导来探究其规律。

首先，在理论模型构建方面，研究者采用了双链Frenkel-Kontorova (FK) 模型。该模型是一种经典的模型，用于描述置于周期势场中的弹性链，在此被推广用于描述二维的双层石墨烯系统。模型的核心是平衡两层石墨烯之间的弹性形变能与层间相互作用能。层间相互作用能采用基于第一傅里叶谐波的势能面（potential energy surface, PES）来近似描述（见原文公式1），该势能面有两个简并的基态堆叠方式（AB和BA），它们之间的最小能量路径通过一个鞍点（SP）。畴壁被视为此路径上的“孤子”（soliton）解，即相对位移在空间中的平滑过渡。畴壁的特性由其伯格斯矢量（Burgers vector）与畴壁法线之间的夹角β决定：β=0°对应拉伸畴壁（tensile domain walls），沿锯齿形（zigzag）方向排列；β=90°对应剪切畴壁（shear domain walls），沿扶手椅（armchair）方向排列。通过求解FK模型，可以得到单个畴壁的形成能（δw）和特征宽度（位错宽度，dislocation width, ld）（见原文公式3-6）。

其次，研究者将单个畴壁的模型扩展到规则的三角形畴壁网络。他们假设在各项同性负载（如双轴拉伸或均匀扭转）下，基态可能对应一个由六个相同畴壁在每个位错节点（dislocation node）交汇而成的等边三角形网络（如图2所示）。通过几何考虑，网络的周期l、层间扭转角φ以及由网络引入的额外双轴应变ε0之间存在明确关系（见原文公式8-9）。研究者进而推导了具有畴壁网络的双层体系相对于共度状态（无畴壁，层间对准）的总相对能量δW_tot的表达式（见原文公式10-21）。这个表达式包含三部分贡献：1）由于网络引入的额外应变ε0导致的弹性能变化δW_el；2）所有畴壁自身的能量贡献δW_dw；3）位错节点（即畴壁交汇处）的能量贡献δW_dn。其中，位错节点被建模为具有特征尺寸ld（畴壁宽度）的六边形区域，假设在该区域内两层石墨烯简单地发生相对旋转和均匀拉伸，层间处于完全非共度（incommensurate）状态，其平均能量可由对整个势能面取平均得到。

在关键参数确定环节，研究严重依赖于第一性原理计算（DFT）的结果。其中最重要的参数是两层石墨烯相对滑移的能垒V_max。作者指出，由于描述范德华相互作用的困难，文献中该值存在较大范围（0.5-2.1 meV/原子）。他们通过比较多种范德华修正泛函对石墨烯和氮化硼一系列性质（剪切模量、体模量、剪切模式频率等）的预测与实验数据，认为vdW-DF2泛函对此类材料的势能面描述最佳，并采用其计算值V_max = 1.61 meV/原子，但同时承认该值可能存在高达40%的不确定性，这是模型预测精度的主要限制因素。其他参数如石墨烯的键长l、单轴应力下的弹性常数k（与杨氏模量和厚度相关）和泊松比ν也取自他们先前基于DFT的计算结果。

基于构建的模型和确定的参数，研究进入了详细的数值分析与讨论阶段，主要分为两种情景：

情景一：上层可自由旋转。此情景下，通过最小化总能量δW_tot对网络周期l和畴壁角β，研究者发现，当底部层双轴拉伸应变ε超过一个临界值ε_c0（计算约为3.0×10^-3）时，体系会发生从共度态到非共度态的相变（commensurate-incommensurate phase transition）。此时能量最优的构型对应β=0°的拉伸畴壁网络，且两层保持对准（φ=0）。这意味着如果上层可以自由转动，在拉伸底部层时，两层倾向于保持对齐，并形成规则的拉伸畴壁三角形网络。他们进一步证明，该相变为二级相变，网络的逆周期l_0^{-1}可作为序参量，其随应变的变化关系为l_0 ∝ (ε - ε_c0)^{-1}，而网络相对能量δW_0 ∝ -l_0^{-2}。

情景二：层间扭转角固定。此情景模拟了实验上常见的具有固定微小扭转角的双层石墨烯。当没有外加应变（ε=0）时，能量最低的网络由β=90°的剪切畴壁构成，其周期l_0与扭转角φ满足简单的几何关系：l_0 = l / (√3 sin φ) ≈ l / (√3 φ)（对于小φ）。例如，φ=0.1°时，l_0 ≈ 142 nm，与实验观测的120-140 nm吻合良好。当施加底部层拉伸应变（ε>0）时，情况变得复杂。通过数值求解能量最小化条件，研究者绘制了最优畴壁角β和网络周期l_0随应变ε和扭转角φ变化的相图（见图5）。结果显示，随着ε增大，β从90°单调减小至0°，即畴壁的取向从扶手椅方向逐渐转向锯齿形方向，特性从剪切型转变为拉伸型。同时，网络的周期l_0（即共度畴尺寸）随ε和φ的增大而减小。研究者还区分了“近似剪切畴壁”（ε很小）和“近似拉伸畴壁”（ε显著大于ε_c0且φ很小）两种极限情况，并给出了相应的解析近似表达式，清晰地描述了β和l_0的变化规律。

除了双轴拉伸，研究者在讨论部分还简要分析了其他类型的机械负载：压缩、弯曲和剪切应变。对于压缩，理论公式与拉伸相同，但需要考虑可能发生的面外屈曲（buckling），这取决于石墨烯与衬底的粘附能。对于弯曲，上下层会承受符号相反的应变，预计也会促进拉伸畴壁的形成，但需要修改FK模型以进行定量描述。对于剪切应变，他们特别分析了在固定扭转角的同时，对底部层施加沿扶手椅方向剪切的情况。此时，几何对称性被破坏，网络中的三角形变为等腰三角形（两条边相等）。他们通过纯几何论证（无需能量优化）给出了等腰三角形网络参数（非等边的长度l、底角α）与扭转角φ、剪切应变τ之间的直接关系：φ = √3 l_0 / l 且 α = arctan[√3 / (φ + τ)]，这展示了在特定负载下网络参数确定的简洁性。

本研究的结论系统而明确：1）若双层石墨烯的上层可自由旋转，在双轴拉伸底部层时，两层将保持对准，并在应变超过临界值后发生向拉伸畴壁三角形网络的共度-非共度相变。2）若层间存在固定扭转角，零应变下观察到剪切畴壁网络；施加拉伸应变会使畴壁逐渐转变为拉伸型，并缩小共度畴尺寸。3）位错节点的形成能估计为151 eV（拉伸壁）和35 eV（剪切壁）。4）在剪切应变和扭转共同作用下，等腰三角形畴壁网络的参数可由纯几何考虑确定。

该研究的科学价值在于，首次通过扩展的双链Frenkel-Kontorova模型，为应变和扭转共同调控下的双层石墨烯畴壁网络提供了系统的理论框架和定量预测。它连接了原子尺度的层间相互作用（通过DFT参数化）与介观尺度的网络结构，并揭示了应变与扭转在决定网络形态中的竞争与协同效应。其应用价值显著：一方面，实验上通过测量畴壁取向和网络周期（例如使用扫描隧道显微镜、透射电镜等技术），可以验证理论预测，进而反推出层间的相对应变和扭转角，作为一种应变/转角探测方法。另一方面，更重要的是，该研究明确指出外部应变可以作为一种有效工具来主动调控畴壁网络的参数（类型、周期），从而为调控与之紧密关联的电子性质（如拓扑通道的密度、导电性等）提供了理论基础和清晰路径。例如，拉伸可以使系统从具有软能隙的剪切畴壁网络转变为绝缘性更强的拉伸畴壁网络，从而可能改变其输运特性。

本研究的亮点突出体现在以下几个方面：首先是方法的有效性与创新性。在原子模拟难以处理大周期（小转角）网络的情况下，研究者成功地将经典的FK模型发展并应用于二维双层系统，建立了一个能够解析描述畴壁网络能量与结构的模型。该模型在畴壁宽度远小于网络周期的假设下成立，恰恰适用于大多数实验关注的“极小扭转角”情形，填补了理论研究空白。其次是结论的系统性与预测性。研究不仅分别处理了自由旋转和固定转角两种基本情形，还给出了详细的参数相图，明确了畴壁类型转变和网络周期变化的定量趋势，这些预测可直接与实验对比验证。第三是对位错节点能量的定量估计，这是以往研究中较少详细讨论的部分。第四是对多种应变类型（拉伸、压缩、弯曲、剪切）的全面考量，特别是对剪切应变下等腰三角形网络的纯几何分析，展现了物理图像的清晰性。最后，研究始终紧扣应用导向，明确指出了理论结果在传感（探测应变/转角）和功能调控（通过应变控制电子性质）两方面的潜在价值，将基础研究与可能的技术应用紧密联系起来。