求助，关于振动弛豫计算工具的思路与问题讨论

不会扣篮的后卫

大家好，请教一个振动弛豫计算的问题。因为目前市面上没有可以计算出振动弛豫的软件，所以我想自己写一个。

代码的理论基础是需要根据激发态势能面解核的薛定谔方程，来求解出最低一些振动能级的振动态波函数，然后根据Fermi's golden rule去计算。
在写的过程中，我发现如果这些底层的计算有现成的代码可以算，只是写个弛豫的接口感觉还是比较容易的。

请问大家是否有推荐的工具/模块可以辅助完成这类底层计算？或者有没有相关经验可以分享？非常感谢！