Revisiting the Extraction of Coupling Strength for Polaron Hopping from $ab~initio$ Approach

Accurately determining the coupling strength between polaron states is essential to describe charge-hopping transport in materials. In this work, we revisit methodologies to extract coupling strengths using ab initio approaches. Our findings underscore the critical role of incorporating anharmonic effects in the model Hamiltonian when analyzing total energy variations along reaction coordinates. Furthermore, we demonstrate that coupling strength extraction based on total energy calculated from ab initio approaches is fundamentally limited due to the stabilization of diabatic states, which do not involve coupling strength. We demonstrate that such limitation exists in both DFT+U and HSE hybrid functionals, which are widely used in the study of polaron transport. Instead, we suggest extracting coupling strength directly from the electronic structure of a system in neutral conditions. The neutral condition avoids overestimating coupling strength due to additional energy splitting between bonding and anti-bonding states from the charging energy. This study highlights the limitations of existing methods and introduces a robust framework for accurately extracting coupling parameters, paving the way for improved modeling of charge transport in complex materials.