In the future, our goal will be to apply the present method to begin assessing non-adiabatic dynamics near plasmonic surfaces, which is nowadays a rather "hot" topic in physical chemistry. B A useful identity for computing determinants using the Schur complement Suppose we divide a general square matrix into 2-by-2 blocks with the diagonal blocks being square matrices, say,