Janus单层Cr2SSe中的应变可调多压电效应与谷电子学

Janus单层Cr2SSe中的应变可调多压电效应与谷电子学

PHYS. REV. B 113, 224428 (2026)

Janus单层Cr2SSe中的应变可调多压电效应与谷电子学

Strain-Tunable Multipiezoeffects in Janus Monolayer Cr2SSe

导读 导读：交错磁体（Altermagnet）同时具备铁磁体的自旋劈裂和反铁磁体的零净磁矩，是谷电子学的理想平台。本文通过第一性原理计算，预测Janus单层Cr2SSe展现出应变可调的多压电效应（压电谷+压电+压磁），并实现选择性谷极化反转和单自旋通道反常谷霍尔效应。核心发现：+4%应变下谷极化达57.6 meV，超过大多数已知ferrovalley材料；-2%到-3%应变下实现价带谷极化的选择性反转，导带谷极化方向不变，为低功耗谷电子学器件设计提供了新范式。

一、前言背景

交错磁体中的谷电子学：从对称性到器件

交错磁体（Altermagnet, AM）是2022年以来凝聚态物理最具革命性的发现之一。其核心特征在于：实空间自旋反平行排列（净磁矩为零），但倒空间能带出现非相对论自旋劈裂。这种独特的对称性特征使得AM材料同时具备铁磁体的自旋极化输运和反铁磁体的超快动力学、无杂散场等优势。

谷电子学（Valleytronics）利用能带中不同谷（valley）的自由度来编码和传递信息。在传统TMD材料中，谷自由度由时间反演对称性保护，需要外磁场或磁掺杂来打破。而在AM体系中，谷简并由晶体对称性（如镜面对称）保护，形成C-paired valleys，使应变工程成为调控谷极化的自然手段。

本文研究的Janus单层Cr2SSe具有P4mm空间群对称性，打破了面外镜像对称（Janus结构），同时保留了面内对角镜面Mxy对称性。这一独特的对称性组合使得该材料同时展现出压电谷、压电和压磁三种多压电效应，成为多功能谷电子学器件的理想平台。

多压电效应：应变同时调控电、磁、谷三个自由度

核心问题：单轴应变如何同时调控Janus AM材料的谷极化、电偶极矩和净磁矩？这三种效应之间是否存在耦合？

方法体系：VASP（DFT+U, PBE+Ueff=3.5 eV）-> Phonopy（声子谱）-> AIMD（热稳定性）-> Wannier90（Berry曲率）-> Berry-phase方法（电极化）-> Bader电荷分析（压磁机制）。

关键发现：(1) +4%应变下谷极化达57.6 meV，超过大多数ferrovalley材料；(2) -2%到-3%应变下实现价带谷极化的选择性反转，导带谷极化方向不变；(3) -3%应变下实现单自旋通道反常谷霍尔效应（AVHE）；(4) 压电、压磁、压电谷三种效应在单一材料中同时可调。

Janus Cr2SSe多压电效应与谷电子学研究流程。VASP DFT+U（结构稳定性、磁基态、电子结构、MAE）-> 单轴应变工程（-4%到+4%，打破Mxy对称性）-> Wannier90（Berry曲率）-> 三种压电效应（压电谷、压电、压磁）-> 单自旋通道反常谷霍尔效应。

二、研究方法

VASP DFT+U：基态电子结构与磁性

计算设置：PBE泛函 + PAW赝势, ENCUT=500 eV, 13x13x1 Gamma-centered k点（弛豫和电子结构）, 真空层30 A。Cr原子Hubbard Ueff=3.5 eV（J=0），与同类Cr基AM材料Cr2SeO一致。收敛标准：能量10^-6 eV，力0.001 eV/A。

磁基态确定：对比FM、AFM1、AFM2、AFM3四种磁构型，AFM1为基态。关键：AFM1中两个Cr子晶格由C2和Mxy对称性连接而非PT对称性，确认为AM序。SOC计算表明带隙基本不变，确认自旋劈裂的非相对论起源。

稳定性验证三重检查：(1) 声子谱无虚频（Phonopy, 5x5x1超胞）；(2) AIMD 300K/5ps能量波动极小；(3) 弹性常数C11=62.9 N/m, C12=8.2 N/m满足Born-Huang判据。

应变工程 + Wannier90 + Berry曲率：谷电子学计算

单轴应变施加方式：沿x或y方向施加-4%到+4%的单轴应变，每次改变晶格常数后重新弛豫原子位置。应变打破Mxy对称性，使X和Y谷不再简并，产生谷极化。

Wannier90紧束缚模型：基于Cr-d轨道和Se/S-p轨道构建MLWFs，用于插值计算Berry曲率。Berry曲率通过Kubo公式计算，是AVHE的核心物理量。注意：原始态Mxy对称性保留时Berry曲率处处为零，只有施加应变后才出现非零分布。

电极化计算：Berry-phase方法（LCALCPOL=.TRUE.），计算面外自发极化。Bader电荷分析：用于解释压磁效应的微观机制--应变导致两个Cr子晶格电荷不等价，通过Hund规则产生不同局域磁矩。

谷极化定义：Delta_C(V)为导带（价带）X谷与Y谷的能量差，正值表示Y谷能量更低

Berry曲率Kubo公式：f_n(k)为费米分布函数，psi_nk为Bloch波函数，v_x/v_y为速度算符

形变势常数：量化能带边缘对应变的敏感度，D值越大表示该能带边缘对应变越敏感

反常速度：Berry曲率作为倒空间有效磁场，在电场E下产生横向速度，驱动AVHE

三、核心结果

图 1：Janus单层Cr2SSe的晶体结构和稳定性。(a) 俯视图（P4mm空间群，Mxy对称性）。(b) 侧视图（S-Cr-Se堆叠）。(c) 声子谱（无虚频，动态稳定）。(d) AIMD模拟（300K/5ps，热稳定）。(e) 无SOC能带结构（间接带隙0.82 eV，AM自旋劈裂）。(f) 含SOC能带结构（带隙基本不变，非相对论自旋劈裂）。

AM序与自旋-谷锁定：C-paired valleys的对称性根源

Cr2SSe的AM序特征：两个Cr子晶格的自旋方向相反，但由C2和Mxy连接而非PT对称性连接。这导致倒空间中出现d波形式的自旋劈裂：X和Y谷的导带和价带分别被相反自旋通道占据，形成自旋-谷锁定（spin-valley locking）。

与传统TMD谷电子学的本质区别：在MoS2等TMD中，K和K'谷由T对称性配对（T-paired valleys），打破谷简并需要外磁场或磁掺杂。在Cr2SSe中，X和Y谷由Mxy对称性配对（C-paired valleys），单轴应变即可天然打破谷简并，无需外部磁场。

这一对称性特征使AM材料成为低功耗谷电子学的理想平台：器件操作仅需机械应变，不依赖电流产生的磁场，功耗极低。

图 2：+2%单轴应变下的能带结构和谷极化。(a,c) x方向应变：价带谷极化-12.4 meV，导带谷极化+28.7 meV。(b,d) y方向应变：价带谷极化+12.4 meV，导带谷极化-28.7 meV。(e,f) 谷极化随应变（-4%到+4%）的演化：导带线性响应，价带非线性响应，x/y方向趋势相反。

压电谷效应：应变方向控制谷极化符号

谷极化量级：+4%应变下导带谷极化达57.6 meV，超过LaBr2（33 meV）、YCl2（22 meV）、CrOBr（44 meV）等已知ferrovalley材料。这一大谷极化值来源于AM材料中C-paired valleys对应变的特殊敏感性。

方向可控性：x方向应变产生正谷极化，y方向应变产生负谷极化。这一符号反转由Mxy对称性决定：Mxy操作交换X和Y谷，因此沿x和y方向的应变产生相反的谷极化演化趋势。

导带vs价带响应差异：导带谷极化线性依赖于应变，价带谷极化非线性响应。根源在于轨道成分差异：导带以Cr-d轨道为主（多方向电子分布，应变敏感），价带以Se-p轨道为主（哑铃型分布，应变相对不敏感）。

图 3：压磁效应和压电效应的应变依赖性。(a) 净磁矩随应变演化：压缩和拉伸应变均诱导微弱净磁矩。(b) Bader电荷分析：应变导致两个Cr子晶格电荷不等价，通过Hund规则产生不等价局域磁矩。(c) 带隙随应变演化。(d) 电极化随应变线性变化：自发极化4.13 pC/m，-4%应变增至4.37 pC/m。

多压电效应的协同机制：三重对称性破缺的统一图景

压电谷效应：由Mxy对称性破缺驱动。无应变时Mxy交换X和Y谷，强制执行谷简并。单轴应变打破Mxy，使X和Y谷不再等价，产生谷极化。

压磁效应：同样由Mxy对称性破缺驱动。Mxy将自旋向上Cr子晶格映射到自旋向下Cr子晶格，锁定零净磁矩。应变打破Mxy后，两个子晶格不再等价，Bader电荷差异导致不等价局域磁矩，净磁矩不再为零。

压电效应：由Janus结构固有的空间反演对称性破缺驱动，与Mxy无关。Cr2SSe本身缺乏面外镜像对称性，存在自发极化。应变仅通过晶格畸变调制极化大小。

三种效应虽然共享应变这一外部刺激，但依赖不同的对称性破缺机制，因此可以实现独立调控而不相互干扰。

图 4：-2%和-3%压缩应变下的选择性谷极化反转和AVHE。(a,b) 能带结构。(c,d) 谷极化：-2%时价带+2.2 meV，-3%时价带-1.6 meV（反转），导带保持负值。(e,f) Berry曲率分布。(g,h) AVHE示意图：-2%时自旋通道混合，-3%时单自旋通道。

选择性谷极化反转与单自旋通道AVHE

选择性反转机制：-2%到-3%应变区间，价带形变势常数发生交叉：X谷价带边缘在-2%时更敏感，Y谷价带边缘在-3%时更敏感，导致价带谷极化符号反转。而导带形变势常数始终保持X谷>Y谷，谷极化方向不变。

单自旋通道AVHE：-3%应变下，价带谷极化反转使空穴通道（Y谷附近）和电子通道（X谷附近）均关联spin-up，实现理想的自旋极化输运。这为低功耗自旋电子学器件提供了新机制。

实验验证建议：可通过自旋分辨ARPES直接观测应变诱导的谷极化反转，通过横向霍尔电压测量验证AVHE，通过偏振分辨PL和MOKE间接验证谷和自旋依赖的载流子响应。

DFT Tips

【DFT Tip 1】AM材料磁基态确定的正确流程

对于交错磁体，磁基态确定比传统FM/AFM材料更复杂。步骤：(1) 先做非磁态弛豫获得初始结构；(2) 根据空间群对称性枚举所有可能的共线磁构型；(3) 在INCAR中手动设置MAGMOM（每个原子指定自旋方向和初始值），不能依赖默认FM初始化；(4) 比较各磁构型总能量，能量最低者为基态。

关键验证：确认基态磁构型满足自旋空间群对称性。AM的特征是Cr子晶格由C2或Mxy等晶体旋转操作连接而非PT操作连接。如果两个子晶格由PT连接，则为传统AFM而非AM。

常见陷阱：直接用FM初始磁矩计算AM材料，VASP可能收敛到亚稳态的FM解而非AM基态。务必手动设置MAGMOM标签，且对不同磁构型逐个测试。

【DFT Tip 2】Cr3+的Hubbard U值选择：3.5 eV vs 4 eV

Cr2SSe使用Ueff=3.5 eV，而Cr2Se2O使用Ueff=4 eV。差值0.5 eV看似不大，但可能影响自旋劈裂大小和带隙。Cr3+（d3）的典型U值范围在3-5 eV之间。

选择依据：(1) 参考同类材料的文献值（Cr2SSe参考Cr2SeO的3.5 eV）；(2) 线性响应方法从头计算U值；(3) 在U+/-1 eV范围内做敏感性测试，确保定性结论不变。

注意：不同化学环境（S vs Se vs O配位）中Cr3+的U值可能不同。S和Se的配位场弱于O，d电子更离域，理论预测U值应略小于氧化物中的值。Cr2SSe使用3.5 eV是合理的，但建议在论文中补充U值敏感性分析。

【DFT Tip 3】2D材料应变计算中的Poisson效应

施加单轴应变时，必须考虑Poisson效应：沿x方向拉伸时，y方向会因Poisson收缩而改变晶格常数。本文似未明确说明是否考虑了Poisson效应，这对应变下的定量结果有影响。

正确处理方式：(1) 固定应变方向的晶格常数，弛豫垂直方向的晶格常数和所有原子位置；(2) 或同时施加双轴应变，但双轴应变不破坏Mxy对称性，无法产生谷极化。

常见错误：只改变一个方向的晶格常数，保持另一个方向不变且不弛豫原子。这会导致非物理的应力状态和错误的能量/电子结构。

【DFT Tip 4】Berry曲率计算的Wannier插值技巧

从DFT能带计算Berry曲率的标准流程：(1) 做自洽计算获取CHGCAR；(2) 非自洽计算生成MLWFs（Wannier90）；(3) 用Wannier插值在密集k网格上计算Berry曲率。对于AM材料，必须包含SOC才能得到非零Berry曲率（原始Mxy对称性保留时Berry曲率为零）。

关键参数：Wannier90中disentanglement窗口需覆盖目标能带（导带底和价带顶附近），inner窗口设置为感兴趣的能带范围。对于Cr2SSe，Cr-d和Se/S-p轨道是构建MLWFs的最小基组。

常见陷阱：Wannier拟合质量差（spread不收敛或能带拟合偏差大），导致Berry曲率在k空间震荡。解决方法是增加初始投影轨道数量、调整disentanglement窗口、或使用SCDM方法自动选择初始投影。

【DFT Tip 5】Berry-phase方法计算2D材料电极化

2D材料的电极化计算需特别注意：(1) 真空层方向不能使用Berry-phase，只能计算面内极化；(2) 面外极化通过平面平均静电势的差异估算（vacuum level difference），不是严格Berry-phase；(3) 本文的4.13 pC/m是面外极化，应通过静电势方法获得。

VASP设置：LCALCPOL=.TRUE.用于面内Berry-phase极化计算。对于面外极化，做LOCPOT计算后处理得到平面平均静电势，取真空层两侧势差。

常见陷阱：将2D材料的"极化"单位pC/m与3D铁电体的"极化"单位uC/cm^2混淆。2D极化是1D线密度，3D极化是2D面密度，两者物理意义不同，不能直接比较数值。

【DFT Tip 6】Bader电荷分析的正确使用与局限

Bader电荷分析将电子密度划分为原子区域，常用于解释电荷转移和氧化态。本文用它解释压磁效应：应变导致Cr(I)和Cr(II)子晶格电荷不等价，通过Hund规则产生不同局域磁矩。

注意：Bader电荷是拓扑划分，不是物理观测量。不同划分方法（Bader vs Hirshfeld vs Mulliken）给出的绝对值不同，只有相对变化趋势有意义。

对于磁性材料，Bader电荷与磁矩的关系不是简单的线性关系。Hund规则是定性指导，定量上磁矩还受配位场、共价性、自旋极化程度等因素影响。建议结合自旋密度图（本文Fig.S3）更直观地展示磁矩的空间分布。

【DFT Tip 7】Phonopy声子谱计算中的超胞收敛

声子谱计算是验证2D材料动力学稳定性的标准方法。本文使用5x5x1超胞（基于2x2磁超胞？），需要确认超胞是否足够大以收敛声子频率。

收敛检查：(1) 比较不同超胞大小（3x3x1, 4x4x1, 5x5x1）的声子谱；(2) 重点关注Gamma点附近声学支是否满足二次色散（面外ZA模）和线性色散（面内TA/LA模）；(3) 确保无虚频（或仅有可忽略的Gamma点附近微小虚频）。

常见陷阱：2D材料中ZA模（面外声学支）在Gamma点附近应呈现二次色散而非线性色散。如果ZA模出现线性色散，说明真空层不够大或数值精度不足。

【DFT Tip 8】形变势常数计算：连接应变与电子结构

形变势常数D = dE/d(epsilon)是量化能带边缘对应变响应的核心参数。计算步骤：(1) 对不同应变值计算能带结构；(2) 提取VBM和CBM在X和Y谷的能量；(3) 对能量-应变关系做线性（或二次）拟合，斜率即为形变势常数。

本文的关键发现：导带X谷形变势常数始终大于Y谷，而价带在-2%到-3%之间发生交叉。这一交叉是选择性谷极化反转的微观起源。

注意：形变势常数依赖于参考能级的选择。通常以真空能级或深芯能级为参考。对于2D材料，推荐使用真空能级作为绝对参考，避免价带顶或导带底的绝对位移带来的歧义。

【DFT Tip 9】单自旋通道AVHE的理论预测与实验差距

本文预测的单自旋通道AVHE是在理想条件下（0K、无缺陷、完美周期性、均匀应变）的理论结果。实际实验中，有限温度、缺陷散射、自旋混合、谷间散射等因素会降低自旋极化率和霍尔响应。

在论文中正确表述：明确指出这是"理想clean-limit band-edge regime"的理论预测，而非实验可达到的精确值。强调核心发现是"应变诱导的自旋分辨谷输运通道分离机制"，而非具体的自旋极化率数值。

实验验证方案：通过栅压分别实现电子和空穴掺杂，选择性探测导带和价带输运通道。结合自旋分辨ARPES直接观测应变下的能带自旋劈裂。

【DFT Tip 10】Janus结构2D材料的vdW修正

Cr2SSe是Janus结构（一侧S一侧Se），面内共价键为主，vdW相互作用相对较弱。但vdW修正对层间距和晶格常数仍有影响，特别是在应变计算中。

建议：(1) 使用optB86b-vdW或DFT-D3等vdW修正方法；(2) 对比有无vdW修正的晶格常数和弹性常数差异；(3) 对于应变计算，如果vdW修正影响较大，需在每种应变下重新启用vdW修正弛豫。

注意：本文未明确提及vdW修正方法。对于单层材料，vdW修正主要影响面外方向，对面内应变计算结果影响较小，但仍建议在方法部分明确说明。

知识扩展

【知识扩展 1】C-paired valleys vs T-paired valleys：谷电子学的对称性范式转变

【理论解释】传统谷电子学（如MoS2）中，K和K'谷由时间反演对称性T连接（T-paired valleys）。T操作同时翻转波矢和自旋，因此K^和K'v具有相同能量。打破谷简并需要破缺T对称性，通常通过外磁场或磁掺杂实现。

【C-paired valleys】在AM体系中，谷简并由晶体对称性（如镜面、旋转）而非T保护。例如Cr2SSe中X和Y谷由Mxy镜面对称性连接。打破这类对称性只需要单轴应变，无需磁场或掺杂，功耗极低。

【方法比较】T-paired方案：需要外磁场（~1-10T）或磁掺杂，器件集成困难，功耗高。C-paired方案：仅需机械应变，可与MEMS/NEMS技术集成，超低功耗。但C-paired valleys目前仅在AM材料中发现，材料选择有限。

【经典参考】Smejkal et al., PRX 12, 031042 (2022) -- AM理论框架；Xiao et al., PRL 108, 196802 (2012) -- TMD谷电子学经典；Mak et al., Science 344, 1489 (2014) -- 谷霍尔效应实验。

【迁移能力】C-paired valleys的概念适用于所有具有适当晶体对称性的AM材料，包括Fe2MoS4、Ti2Se2S、VPSe3等。

【知识扩展 2】多压电效应：从单一效应到多场耦合

【理论解释】传统的压电效应（1880年Curie发现）仅涉及应变-电极化的耦合。压磁效应（piezomagnetism）涉及应变-磁化的耦合，在AM材料中由于补偿磁序极为敏感。压电谷效应（piezovalley effect）是近年提出的新概念，涉及应变-谷极化的耦合。

【多压电效应（multipiezo effect）】指应变同时调控电、磁、谷三个自由度的协同效应。Cr2SSe是目前少数同时展现三种压电效应的单材料体系。其物理根源在于：AM序提供磁自由度，Janus结构提供电偶极矩，C-paired valleys提供谷自由度。

【历史发展】压电效应（1880）-> 压磁效应（1960s）-> 多铁性（2000s）-> 压电谷效应（2020s）-> 多压电效应（2025）。每一步都代表了对材料中多自由度耦合认识的深化。

【经典参考】Curie & Curie (1880) -- 压电效应发现；Ikhlas et al., Nat. Phys. 18, 1086 (2022) -- 反铁磁压磁效应；Jiang et al., APL 126, 053102 (2025) -- AM多压电效应。

【迁移能力】多压电效应概念可推广至其他Janus型AM材料（Cr2SeO、Cr2SO、V2Se2O等），以及具有类似对称性的二维多铁材料。

科研经验

【科研经验 1】应变工程计算中的系统误差控制

问题：应变计算中，不同应变值下的计算精度是否一致？应变范围是否合理？非线性效应的物理来源是什么？

原因：(1) 不同应变下晶格常数变化，平面波基组的等效截断能改变，导致Pulay应力；(2) 应变过大（>5%）时，谐波近似失效，需考虑非谐效应；(3) 2D材料在压缩应变下可能出现buckling，改变电子结构。

解决方案：(1) 对所有应变值使用相同的ENCUT（以最大晶格常数为基准），避免Pulay应力；(2) 应变范围控制在+/-5%以内（大多数2D材料在此范围内弹性响应线性）；(3) 对压缩应变检查是否出现面外buckling，如有则需在弛豫中允许原子沿z方向移动。

建议：在论文中展示不同应变下的结构参数（键长、键角、层间距），佐证结构变化的合理性。对于非线性响应（如价带谷极化），务必提供物理机制解释。

【科研经验 2】AM vs 传统AFM的DFT区分陷阱

问题：在DFT计算中，如何确定一个材料是AM而不是传统AFM？仅凭能带自旋劈裂是否足够？

原因：很多传统AFM在SOC下也会出现自旋劈裂，但这是相对论效应（~meV量级），而非AM的非相对论自旋劈裂（~100 meV量级）。仅看能带劈裂存在无法区分AM和AFM+SOC。

解决方案：(1) 分别计算无SOC和含SOC的能带，如果无SOC就出现显著自旋劈裂，则为AM特征；(2) 分析磁构型的对称性：两个子晶格是否由晶体旋转（而非PT或平移+反演）连接？(3) 检查自旋劈裂的k空间分布：AM的d波劈裂有节点，AFM+SOC的劈裂通常无节点。

建议：在论文中同时展示无SOC和有SOC的能带，并明确标注自旋劈裂的对称性（d波/g波/i波）。仅凭能带劈裂声称AM是不够严谨的。

【科研经验 3】Berry曲率符号与AVHE的对应关系

问题：Berry曲率在X和Y谷符号相反，如何判断AVHE的方向和幅值？负Berry曲率是否意味着负Hall电压？

原因：AVHE的Hall电压正比于Berry曲率在占据态上的积分，不仅取决于Berry曲率的符号，还取决于费米能级的位置（doping类型和浓度）。对于空穴掺杂（费米能级在价带），积分来自价带Berry曲率；对于电子掺杂（费米能级在导带），积分来自导带Berry曲率。

解决方案：(1) 分别计算空穴和电子掺杂下的AVHE；(2) 在论文中明确标注费米能级的位置；(3) 对于Berry曲率符号相反的谷，AVHE的贡献相互竞争，净效应取决于费米能级处两个谷的态密度权重。

建议：在方法部分提供详细的AVHE计算公式，说明是整数Hall电导（量子化）还是非量子化的Hall电导。对于非量子化的AVHE，需明确区分与量子反常霍尔效应（QAHE）的不同。

如果是我，我还会继续算

【继续算 1】HSE06杂化泛函验证带隙与自旋劈裂

为什么值得算：PBE带隙0.82 eV可能被低估30-50%，HSE06可给出更准确的带隙和自旋劈裂值。带隙的准确性直接影响谷极化作为器件工作温度上限的估计。

能回答的问题：HSE06下的自旋劈裂模式和谷极化是否与PBE定性一致？带隙修正后AVHE的能量窗口是否变化？

适合体系：所有半导体AM材料。输入：HSE06计算（ENCUT=500 eV, 减少k点到5x5x1）

【继续算 2】Monte Carlo模拟磁转变温度

为什么值得算：本文未计算Cr2SSe的磁转变温度（Neel温度）。对于实际器件应用，需要知道AM序能维持到多高温度。

能回答的问题：Cr2SSe的Neel温度是多少？应变对磁转变温度有何影响？压磁效应在接近T_N时是否增强？

适合体系：所有2D磁性材料。输入：DFT提取交换耦合J_ij，EspinS或VAMPIRE进行Monte Carlo模拟到中等。

【继续算 3】载流子掺杂效应：超越应变调控

为什么值得算：本文仅通过应变调控谷极化，但静电栅压（载流子掺杂）是另一种重要的调控手段。掺杂可以移动费米能级，改变AVHE的大小和符号。

能回答的问题：电子/空穴掺杂如何影响谷极化和Berry曲率？是否存在最优掺杂浓度使AVHE最大化？不同掺杂浓度下选择性谷极化反转的阈值是否改变？

适合体系：所有半导体AM材料。输入：在DFT中通过增减电子数（NELECT标签）模拟掺杂，或使用刚性能带模型。

【继续算 4】非线性光学响应：SHG与谷极化的光学探测

为什么值得算：Janus结构缺乏反演对称性，天然具有二次谐波产生（SHG）响应。SHG强度对谷极化和应变敏感，可作为谷极化的全光学探针。

能回答的问题：Cr2SSe的SHG系数有多大？应变如何调控SHG？SHG能否区分X和Y谷的极化状态？

适合体系：所有非中心对称2D材料。输入：VASP计算非线性光学极化率（LEPSILON=.TRUE. + 非线性模块）。

【继续算 5】双层/异质结中的谷电子学：层间耦合效应

为什么值得算：单层Cr2SSe的谷极化在器件中可能受衬底影响。构建双层或异质结（如Cr2SSe/BN、Cr2SSe/Cr2SeO）可以研究层间耦合对谷极化和AVHE的影响。

能回答的问题：双层Cr2SSe的层间磁耦合是什么类型？异质结中谷极化是否被保持？层间耦合是否产生新的拓扑相？

适合体系：所有2D vdW材料。输入：构建异质结超胞，DFT弛豫+电子结构计算到高。

【继续算 6】TB2J提取交换耦合 + 自旋波激发谱

为什么值得算：从DFT提取J参数可以构建Heisenberg模型，进一步计算自旋波谱（磁振子）。磁振子对AM材料的热输运和自旋输运有重要贡献。

能回答的问题：Cr2SSe的交换耦合参数是多少？AM序的自旋波谱与AFM有何不同？自旋波是否携带Berry曲率？

适合体系：所有磁性材料。输入：VASP + TB2J + SpinW/Sunny。

Shen, Tan, Yao, Liao, Dong | Phys. Rev. B 113, 224428 (2026) | 交错磁体 多压电效应 谷电子学 AVHE