CP2K固定核实时电子响应：三方向24fs完整张量与续算

1. 快速绘图：解压进入cp2k-rt-tensor，bash run.sh。这只处理资料包中鸢算科技集群真实计算的三方向24fs数据，不重新跑电子传播。results、plots为本次后处理输出；reference-plots是原示例图。
2. 短原生练习：在获分配CPU作业内，bash run.sh smoke。x方向从基态开始10步，再用成对的excitation-1.restart和excitation-RESTART.rtpwfn续到第20步。总共0.05fs，只学习完整重启，不用于吸收谱。已有native/smoke-x目录时拒绝覆盖，请另解压新目录。
3. 完整重新计算：sbatch job.slurm，三个方向依次各9600步、时间步0.0025fs，共每方向24fs。cpu-8173m按56物理核整节点分配和计费，脚本每次一个MPI、8个OpenMP线程，不使用GPU。它是CPU模板；请确认相应分区当前可用并遵守集群提交规则。
4. 完整运行后，在相同环境中python -B analyze.py --native，然后python -B plot.py --native，从自己新生成的native/full24-x、y、z/FULL-FIELDS.npz读取，写results-native和plots-native。默认不带--native始终读取资料包inputs中的原示例数据，避免把原图当成自己的新计算结果。

示例使用CP2K2026.2，PBE/DZVP-MOLOPT-SR-GTH/GTH-PBE、12Å孤立盒、600Ry/REL_CUTOFF50、固定水分子核坐标、ETRS/PADE。当前CPU构建oneAPI2024.2，绘图sisl0.16.4。原始24fs例子先由现有P100构建开始，后续CPU检查点延续前做完整10步、三分量偶极和能量对照；这不是新的CPU/GPU速度比较。

每次kick产生张量的一列：alpha[能量,响应方向,脉冲方向]。扣除真正的t=0原生偶极，Debye换成原子单位，实际电场冲量kappa=2*pi*0.0001/(12Å/Bohr)，约2.77076538e-5。CP2K2026.2实际施加exp(-i*kappa*r)，正频变换用exp(+i*omega*t)，电偶极响应除以正kappa；不能把打印的电子波函数相位方向或DELTA_PULSE_SCALE直接当作电场冲量。

采用0.3eV指数阻尼，5–30eV共1001点，直接复合梯形积分；17个真正FFT频率用独立均匀FFT对照。12fs与24fs的平均虚部相对L2约0.34924%，完整复张量约0.41924%。这只表明此固定设置和阻尼下的窗口稳定性，不代表基组、网格、泛函或全部材料均收敛。0.025eV画图间隔不是0.025eV物理分辨率；24fs基本能量尺度h/T约0.1723eV。0.3eV不是测量的激发态寿命。

水位于yz平面，非对角小项保留原值，不据数值残差推断物理对称性破缺。示例最大平均峰约24.775eV属于PBE有限基组教学结果，不当作实验峰预测。本篇没有全电子芯态、核运动或非绝热动力学计算。0.005fs与0.0025fs的旧4fs对照差6.824%未满足原5%目标；另做0.00125fs与0.0025fs得到1.731%，实际保留0.0025fs。所有失败和原生原件保留在内部，不通过平移或拟合改变这些差别。

本包仅此固定水体系。run_native.py会与此结构的原始参考核对，换自己的结构、方法或盒子应另建项目，重新定义自己的收敛与物理比较，不能沿用固定水参考作为其他体系的通用程序。仅含输入、数值数组与自己的脚本，无波函数/检查点二进制；运行时从现有模块读取基组和赝势。input中的electronic_density是原生带符号电子电荷-8，而不是负电子数。
