喜讯！PWmat新增体材料非绝热动力学模块

非绝热动力学（NAMD）是描述载流子冷却过程的重要方法。龙讯旷腾PWmat的NAMD方法已经被用在不同的纳米材料、缺陷以及界面等系统上。但以前的这些NAMD计算，都是基于基态的AIMD模拟，而不是显性的电声耦合计算。但对有的问题来说，需要算很大的系统，这尤其表现在对周期性材料的问题上。虽然PWmat已经有一个跨k点的基于AIMD以及NAMD的方法，但此方法还是需要通过超胞来计算，对于纯体材料也许并不是最有效的方法。

对于周期体系，另一个途径是通过Wannier函数，以及声子模式来计算电声耦合，这些功能都已经在PWmat的其他模块中实现了。如果有这样的电声耦合系数，我们就可以在k空间来执行NAMD。我们可以根据声子模式在Harmonic oscillator的近似下来描述原子的轨迹，而利用电声耦合系数来描述电子哈密顿量的变化。理论上讲，我们可以通过以前的P-matrix来做相应的NAMD计算。但基于k-point的电子态数目有可能是巨大的，比如20x20x20个k-mesh就有8000个k-point，而每个k-point上有可能有10个能带，这样总共就有8万个电子态。很多情况下20x20x20 k-mesh可能还不够。这种情况下，P-matrix计算就很慢了。这时我们宁可用单波函数计算。这需要一套stochastic dynamics方法。文献中我们有DISH方法， 但我们试用下来并不那么满意。

特此我们开发了一个新的Single Wave Function Boltzman方法。此方法基于波函数的塌陷，并通过施加Boltzmann factor来达到细致平恒。另一方面，我们采用了两套k-mesh格点的方法来加速NAMD计算。在粗格点上，我们利用电声耦合系数来描述k点之间的coupling，但在细格点之间，我们利用deformation potential来描述带内的耦合。我们发现，利用这种方法，我们可以大大加快NAMD的计算。