本文档旨在说明“单晶 XRD 定向衍射模拟与晶格反演计算器”的核心功能、物理理论背景以及操作规范。该程序主要用于材料晶体学研究，建立从微观晶胞结构到宏观 X 射线衍射（XRD）图谱的数值映射，并通过实验数据的非线性最小二乘法优化，实现晶格参数的逆向推导。
一、 核心功能概述
本计算器具备以下三大独立计算模块：
- 晶体结构解析：读取标准 VASP POSCAR 或 CONTCAR 文件格式，解析元素种类、分数坐标及晶格矢量，并转换为标准空间矩阵。
- 定向衍射模拟：基于指定的密勒指数族（如 00l 或 h00），计算倒易点阵间距，通过布拉格定律与结构因子公式，预测理论 2θ 峰位及系统消光现象。
- 晶格参数反演：读入实验衍射峰位数据，利用基于自适应步长的数值梯度下降算法，对晶格参数 a、c 及测角仪零点偏移量（Zero Shift）进行联合优化，得到最小二乘解。
二、 物理背景与理论公式
本程序的计算引擎严格遵循运动学衍射理论，各项物理量的计算逻辑如下：
1. 倒易晶格与晶面间距 (d-spacing)
由实空间晶格基矢 a, b, c，定义单胞体积 V = a · (b × c)。依据晶体学惯例（不含 2π 因子），构造倒易空间基矢：
a* = (b × c) / V
b* = (c × a) / V
c* = (a × b) / V
对于指定的衍射面 (hkl)，其倒易空间矢量 Ghkl 及其对应的晶面间距 dhkl 满足：
Ghkl = ha* + kb* + lc*
dhkl = 1 / |Ghkl|
2. 布拉格定律 (Bragg's Law)
当入射 X 射线波长为 λ 时，满足相长干涉的衍射角 θ 受布拉格方程约束：
2dhkl sin(θ) = λ
若 λ > 2dhkl，该晶面在几何上无法发生衍射。
3. 结构因子与系统消光 (Structure Factor)
衍射峰的本征强度不仅由几何间距决定，更取决于晶胞内原子的相对位置所导致的干涉叠加。复数结构因子 F(hkl) 定义为：
F(hkl) = Σ fj(s) exp[ 2πi (hxj + kyj + lzj) ]
其中，(xj, yj, zj) 为原子的分数坐标，fj(s) 为原子散射因子。本程序采用 Cromer-Mann 参数化模型计算随散射矢量 s = sin(θ)/λ 衰减的 fj(s)。衍射强度正比于 |F(hkl)|2。若 |F(hkl)|2 趋近于 0，则判定为系统消光（Forbidden）。
4. 误差函数与晶格反演模型
仪器测量中不可避免地存在测角仪零点偏移（Δ2θ）。程序利用实验峰位 (2θexp) 反推晶格参数，构建的最小二乘目标损失函数 χ2 为：
χ2(a, c, Δ) = Σ [ 2θcalc, i(a, c) + Δ - 2θexp, i ]2
三、 物理参数层级关系
程序内部计算严格遵循单一流向的数据映射，明确区分几何可达性与干涉消光特性：
[POSCAR / 晶格常数]
↓ 决定
[倒易空间矩阵]
↓ 决定
[晶面间距 dhkl]
↓ 结合波长 λ (Bragg)
[理论衍射角 2θ] （决定峰的横坐标位置）
--------------------------
[POSCAR 原子分数坐标] + [原子种类]
↓ 结合各原子 Cromer-Mann 参数
[结构因子 F(hkl)]
↓ 决定
[本征衍射强度 / 系统消光] （决定峰是否存在及纵坐标相对强度）
四、 标准操作规范
步骤 1：导入结构模型
在“数据输入”页面，粘贴标准的 VASP POSCAR 内容。程序兼容 Fractional (Direct) 与 Cartesian 坐标，并可正确识别缩放因子（Scale Factor）。解析完毕后，可根据需要设定 X 射线靶材波长。
步骤 2：正向衍射模拟
切换至“反射模拟”页面。选择扫描类型（如二维材料常用的 00l 层间反射，或输入自定义密勒指数）。点击计算后，程序将输出允许衍射的峰位图谱及包含结构因子实部、虚部的详细数据表，并排除消光晶面。
步骤 3：实验数据逆向拟合
确保步骤 2 理论数据已生成。切换至“实验峰拟合”页面，粘贴实验测得的峰位数据（支持单列 2θ 或包含强度的双列格式，程序将自动提取角度）。勾选需要优化的参数（a, c 晶格常数或 Zero shift），点击执行。程序将自动匹配实验与理论峰位，利用梯度下降算法输出残差最小化的真实晶格参数。
五、 物理模型局限性声明
使用本计算器时，需明确其基于理想运动学衍射理论的科学边界：
- 相对强度误差： 程序计算的强度正比于 |F(hkl)|2。实际测量的绝对强度受洛伦兹-偏振因子（Lorentz-Polarization）、吸收效应（Absorption）、择优取向（Preferred orientation）及 Debye-Waller 因子强烈影响。本程序结果不可直接作为粉末衍射法（含多重性因子）强度的定量依据。
- 衍射几何： 模型假设严格的单晶定向布拉格衍射几何，未引入粉末晶粒取向统计模型。
- 反演收敛性： 非线性最小二乘法对初值存在依赖。若输入 POSCAR 的初始晶格参数与实验真实值偏差过大，可能导致算法收敛于局部极小值。
