最小面跳跃法

📁 version1

文件

• model.nb：Mathematica脚本

• fssh_v1.py：Python脚本

运行方法

1. 使用Mathematica打开并运行 model.nb，生成以下文件：

   • Tully-model1-e1.dat

   • Tully-model1-e2.dat

   • Tully-model1-d12.dat

   • potential.pdf

2. 然后运行Python脚本：

   1	python fssh_v1.py

   脚本将加载上述 .dat 文件并执行表面跳跃模拟，生成：

   • single.jpg：势能和耦合图

   • model-1-prob.jpg：跳跃统计图

📁 version2

文件

• fssh_v2.py：整合后的Python脚本

运行方法

1

python fssh_v2.py

该脚本将：

• 直接计算电子势能面和耦合项

• 执行FSSH模拟

• 生成与version1相同的输出图像：

  • single.jpg

  • model-1-prob.jpg

运行结果

参考

• Tully J C. Molecular dynamics with electronic transitions[J]. The Journal of Chemical Physics, 1990, 93(2): 1061-1071.

• mudslide

• FSSH

• Jain A, Sindhu A. Pedagogical overview of the fewest switches surface hopping method[J]. ACS omega, 2022, 7(50): 45810-45824.