在J. Chem. Phys. 162, 110901 (2025)这个用CP2K跑水合电子AIMD模拟的工作中，表征水合电子分布的指标除回转半径外还有个“非球形度”(asphericity)，和http://sobereva.com/661提到的球形度只考虑表面积和体积不同，定义为<A>=<Tr^2 – 3M>/<Tr^2>（此文中公式错写为Tr^3了，对照先前文献应该是Tr^2），其中Tr和M都能用回转张量的本征值算出来。接近球对称的形状的<A>小得几近于0，而细长的形状的<A>接近1。

既然multiwfn已经能给出二阶矩本征值，输出Tr^2和M应该是小菜一碟；这样逐帧批处理时外部指令提取出各帧Tr^2 – 3M和Tr^2，再分别取时间平均然后相除，就能得到一条轨迹的过剩电子的非球形度了。

编辑：这帖提到的方法似乎也有点关系

社长威武！