求教，使用vasp软件如何计算跃迁偶极矩？

Shana

各位老师好，小弟最近处于兴趣，在研究偶极近似下的辐射跃迁过程。查阅资料发现，辐射速率正比于跃迁偶极矩的平方。跃迁偶极矩μ=<i|e·r|f>，i和f是初末态。小弟检索到一些计算方法，想请各位老师指点一下：
1. vaspkit 713 能够计算跃迁偶极矩。但是小弟不太懂代码，不知道vaspkit计算的是哪一部分？是μ？是<i|r|f>?还是<i|p|f>？vaspkit输出的单位是德拜，应该是<i|e·r|f>？
2. 小弟关心的问题是载流子在缺陷上辐射复合的过程？以GaN中C取代N (C_N) 为例，假设是导带电子到C_N^0的复合过程，那么初态是C_N^0+e^-1，末态是C_N^-1。使用vaspkit计算跃迁偶极矩时，指定的初态band是导带，末态是C_N^-1的band?
如果是-1价的C_N和空穴复合的过程，那么初态是C_N^-1+h^+，末态是C_N^0。使用vaspkit计算跃迁偶极矩时，指定的初态band是C_N^-1的band?末态是价带中的空穴band?
3. 小弟在文献中看到，涉及自旋态的跃迁偶极矩没法计算。这是指，自旋三重态到自旋单态的跃迁偶极矩不好计算吗？小弟对束缚激子的复合比较感兴趣，求老师们指点。
4. vaspkit只能读取std计算出的WAVECAR，有没有能够直接计算跃迁偶极矩的软件呢？
请各位老师指点！