[求助] 动力学跑出这样的结构，是什么原因啊？ 已有1人参与

RT，跑铁的超胞动力学出来一个这样的结果，感觉是错了，请各位大神分析一下是什么原因，不胜感激，，，    输入参数： md conp conse qok nomod pres spat pot efg ensemble npt tau_thermostat  0.1 ps tau_barostat  0.1 ps temperature 1500.000000 K pressure 0.000101 GPa equilibration 10 ps timestep 2 fs production 50 ps write 100 sample 100 title GULP calculation from Materials Studio for Fe9826 end cell 48.728800 48.728800 48.728800 90.000000 90.000000 90.000000 fractional Fe1 core 0.000000 0.000000 0.000000 0.000000 1.000000 0.0 1 1 1 Fe1 core 0.029412 0.029412 0.029412 0.000000 1.000000 0.0 1 1 1 输入的结果以几张附图所示，，， 结构.jpg constant of motion.jpg pressure.jpg temperature.jpg

回复此楼