无序缩短晶格相干长度而非将其关闭
论文追问的是:一旦存在随机分布的杂质,晶格组织还能在多大范围内存活,以及这个距离如何随无序增强而收缩。库仑排斥倾向于三角排列,而杂质倾向于不规则的密度图案,两者在同一批电子上相互竞争。这里报告的答案是一个渐变过程,而不是一条边界:无序增强会削弱晶格关联、展宽结构因子,并把测得的相干长度压缩向原始晶格间距。
这一表述之所以重要,是因为它排除了那种直觉图像,即存在一条有限无序的界线把完美晶体与非晶固体分开。作者把他们的数值演化与 Imry-Ma 和 Larkin 图像联系起来:任意弱的随机钉扎都会在足够长的尺度上破坏原始平移有序,而无序越强,畸变尺度越短。因此弱钉扎的畴在有限观测窗口内仍可能看起来像晶体,这也是论文把这些图视为非晶化程度、而非相边界的原因。
两个微观模型检验同一个结构问题
论文并行运行了两套计算,各自拥有独立的哈密顿量、干净参考态、无序控制参数和数值检验。量子研究使用连续 Hartree-Fock 理论,追踪短程随机势中的自洽电子密度,其超胞包含九个电子,在三角晶体的每个原胞方向上跨越三个原胞格矢。一个无序和离散化完全相同的纯直接 Hartree 计算充当对照,用来分离交换作用,因为干净的纯直接 Hartree 密度是均匀的,而所选的 Hartree-Fock 态在 r_s = 3、4、8 和 16 时是晶态。
经典研究使用分子动力学和零温弛豫,追踪具有库仑相互作用和开放静电边界的点状电子,它们处在随机分布的吸引势阱中。它解析粒子构型、配位缺陷、位置与取向关联以及笼逃逸。两种无序方案刻意互补:Hartree-Fock 增大一个固定高斯随机势的幅度,而分子动力学在一个嵌套的泊松杂质场中增加吸引势阱的浓度,同时保持其物理深度和半径不变。两套计算都不对无序实现做平均,因此两个无序轴是分别定义的,比较针对的是共同的物理机制,而不是两个无序轴之间的数值换算。
量子密度权重离开原始布拉格矢量
在量子计算中,非零超胞倒格矢处的弹性密度强度衡量自洽密度的静态空间组织。它是由密度算符期望值构建的相干贡献,带有每个原胞的归一化,而不是完整的等时量子结构因子,后者还包含本工作未计算的连通交换与涨落项。在这一界定下,无序增强会把权重从原始的六重壳层转移到其他超胞倒格矢,而周期性边界条件使这些谱线保持离散,未施加任何人为的峰展宽。
论文还报告了一个操作性的胞内相干长度。在 r_s = 4 时,它从弱无序下大于 1.5a 的未分辨值,变为 W/E_F = 2 和 4 时的约 1.1a 和 0.9a。超胞在每个原胞方向上只跨越三个原始晶格间距,因此密度相干诊断是在该范围的一半内评估的,并被明确区别于渐近关联长度。量子图描述的是一个重复原胞内计算得到的密度,并不能确立独立无序原胞中关联的热力学衰减。
钉扎密度上升使经典结构因子展宽
在电子密度 n_e = n_0 的经典计算中,随着杂质密度增加,结构因子从六个尖锐峰演变为弥散散射,位置关联和取向关联都随之减弱。测得的平移相干性从弱钉扎时的超过 5a,下降到杂质浓度 c = 1/8、1/4 和 4 时的约 2.7a、1.2a 和 0.8a。强无序端点是机械稳定的非晶构型,只保留局部堆积有序,而带宽检验确立的是量级层面的结果,而不是一个低于一个晶格常数的精确距离。
在固定势阱半径和深度下的双密度图把这一结果扩展到不同电子密度,包含五十个计算点,晶态由同时满足取向有序阈值和六重分数阈值来定义。作者强调这些等值线是有限样本的结构渐变。他们还指出,布拉格分数和第一壳层相干度量都具有取向选择性,因此一个旋转或重构的晶体也可能在固定参考矢量处损失权重,而并未失去全部有序。
更强的钉扎提高重排温度
在固定杂质构型下加热,可以分离两种容易被混为一谈的效应。主要笼逃逸区间从干净系统中的 k_B T/E_0 = 0.005-0.007,上升到 c = 4 时的 0.020-0.030,而在 c = 4 时晶态有序在零温下已经消失。随机势阱破坏位置和取向相干性,缩小可辨认晶体的稳定范围,但它们也制造了粒子重排的势垒,使已经不规则的构型在更高温度下仍被笼住。
论文对这类迁移率渐变是什么保持谨慎。它是滞后 100 和 400 t_0 的有限观测时间渐变,不是扩散系数,也不是热力学玻璃转变;作者指出,在哈密顿量中存在显式外无序时并不存在玻璃,因为哈密顿量缺乏平移不变性。所报告的误差棒跨越采样得到的交叉区间,不是统计置信区间。采样检验在这里很重要:一个单独的干净 N = 572 计算给出的表观阈值不同于主要的 N = 288 方案,因此可辩护的干净尺度大致是 theta = 0.005-0.010,而不是一个高精度数值。
Larkin 估计与测得的相干尺度相符
论文使用干净的弹性刚度和已知的无序协方差,为两个模型推导了弱钉扎 Larkin 估计,未做数值无序平均。对于固定 r_s 的 Hartree-Fock 幅度扫描,弱钉扎标度为 R_c 正比于 W^-1,在 r_s = 3、4、8 和 16 时的估计约为 3.2a、2.5a、1.9a 和 2.3a。在 r_s = 4 且 eta = 0.5 时估计约为 5a,已经大于模拟的量子胞;在 eta = 4 时形式值约为 0.6a,作者将其解读为集体晶体描述的失效,而不是分辨出一个亚晶格弹性畴。
对于经典浓度扫描,原始三角库仑晶格的剪切模量为 mu_cl = 0.245065 (e^2/4 pi epsilon) n_e^{3/2},而 q = 1 时的紧致势阱满足 2R_p < a,因此独立的泊松杂质在不同原始晶格位点上给出不相关的力。代入实际的势阱深度和半径,得到 R_c/a 约为 0.63/sqrt(c),线性高斯估计约为 0.87/sqrt(c)。在 c = 1/8 时,估计为 R_c 约 1.8a、高斯值约 2.5a,而测得的相干性约为 2.7a。作者称这是无拟合参数量级上的一致,而不是对前因子的验证,并指出每个 c 下 Larkin 圆盘只容纳约 1.5 个势阱,因此 R_c 只是跟随平均势阱间距。两种扫描具有不同的幂律,因为它们的随机力方差分别按 W^2 和按钉扎密度乘以势阱深度平方标度。
大密度调制并非晶态有序的证据
论文最尖锐的观点之一涉及一个容易被误用的诊断量。在量子研究中,密度对比在晶格组织丧失的同时增大,因此大调制不能充当晶态有序的替代指标。晶体和非晶钉扎集合体都可以具有强烈调制的密度,这也是作者主张空间关联及其傅里叶谱检验晶格组织、而迁移率测量检验重排的原因。
这一区分也推动了论文针对文献中某个具体主张的收尾论证。作者得出结论:说外无序增强 Wigner 晶化并降低临界 r_s 是误导性的:无序可能增强某些短程有序,但晶体特有的长程有序总是被无序抑制,系统变为非晶固体。他们给出的理由是,如果无序增强会促进晶化,那么只需向系统添加杂质就能任意降低临界 r_s,这说不通,因为高杂质密度必然导致非晶结构。论文自身对此的证据是:两个模型中相干长度都在收缩,而经典模型中重排温度在上升,二者共同描述的是一个被钉扎的非晶固体,而不是一个更好的晶体。
这些结果确立了什么、未确立什么
量子计算使用的是平均场 Hartree-Fock 固体,作者承认关联计算把干净晶化放在比 Hartree-Fock 低得多的密度上,不同蒙特卡洛研究报告的临界 r_s 值为 37 加减 5、31(不确定度为 1)和 35(不确定度为 1)。他们认为平均场技术保持了物理的透明性,即使实际关注的 r_s 区间发生移动,而他们认为这对于定性理解相互作用与无序的竞争是可以接受的。经典计算限于 q = 1;在足够高的电子密度下,必须保留近邻位点的力关联,因为 2R_p < a 不再成立。
若干报告的量是操作性的,而非热力学的。胞内相干诊断是在超胞范围的一半内评估的,经典结构等值线是有限样本渐变,迁移率渐变依赖观测时间,强无序相干尺度是量级层面的结果。论文还指出,弱无序下的重入有序区间曾被 Nelson 预言,并在更早的工作中被发现,但这些轨迹并未分辨出重入有序区间或精确的熔化温度。因此最有用的结论是一个比较,而不是一个阈值:晶格相干性和抗运动能力是不同的性质,在这篇论文中无序抑制前者而增强后者,所以一个被钉扎的二维电子固体可以看起来更不像晶体,却仍然更难重排。
本文引用的文献
Disorder-driven loss of crystalline coherence in two-dimensional electron solids: Quantum Hartree-Fock and classical molecular dynamics
在 Hartree-Fock 和经典分子动力学两种描述中,淬火无序增强都会削弱晶格关联,并把测得的相干长度压缩向原始晶格间距;而更强的钉扎即使在晶态有序消失之后仍会提高粒子重排的温度。
