寫了個小程式，與Gaussian聯用實現quasiclassical direct dynamics。 程式本身毫無難度，寫它的唯一目的是不用看滿是bug慘不忍睹的Singleton組的程式碼。

Quasiclassical

direct dynamics

原理見 doi：10。1021/jacs。6b12017。簡單來說就是對過渡態的每個振動維度賦其零點能來近似描述quantum effect。施加此種近似之後，就可以簡單地用NVE下的BOMD來描述體系。私以為這種魔性的近似方法能廣為流傳，主要還是因為大家都不想用Path Integral。

專案地址：

https：//

github。com/WangXinyan94

0/quasidyn