From: Jiaguo Yu <jia...@ya...> - 2016-06-11 16:29:54
|
<div>Dear cclib users,</div><div> </div><div> From <a href="http://cclib.github.io/methods.html," target="_blank">http://cclib.github.io/methods.html,</a> cclib can perform Mulliken population analysis (MPA) and store the results in aoresults, fragresults and fragcharges. The fragcharges can be used to compute Mulliken atomic charges. Nonetheless, how can one compute Mulliken atomic spin densities (for example, from fragresults) ?</div><div> </div><div> The blunt snippet below shows that "alpha" and "beta" variables contains almost the same values...</div><div><pre>sum_alpha = 0.0<br />sum_beta = 0.0<br />for i in range(len(m.fragcharges)):<br /> alpha = 0.0<br /> beta = 0.0<br /> for j in range(len(m.fragresults[0])):<br /> alpha += m.fragresults[0][j][i]<br /> beta += m.fragresults[1][j][i]<br /> sum_alpha += alpha<br /> sum_beta += beta<br /> print("sum_alpha=%10.5f, sum_beta=%10.5f, alpha=%10.5f, beta=%10.5f" % (sum_alpha, sum_beta, alpha, beta))</pre></div><div> </div><div>Best regards,</div><div>Jiaguo Yu</div> |