Friday, 23 October 2015

computational chemistry - Is it possible for the kinetic energy integral to be negative?



Is it possible for the kinetic-energy integral, Tij, to be negative? I was messing around with some HF code and found that the integral became negative on some off-diagonal terms. (This could also be a bug with the code.)


If the kinetic energy integral is negative, what would the physical meaning be?



Answer



Generally, the off-diagonal elements of matrices, especially in physical systems, are interpreted as the coupling between whatever the $i$ and $j$ elements correspond to. So, in this case, a negative element of $T_{ij}$ corresponds to a negative kinetic coupling between atomic orbital basis functions. In plain language, this means that basis functions $i$ and $j$ tend to mutually lower the kinetic energy of an electron placed in one of those orbitals. I am not sure it really makes sense to give much of an interpretation to this, however, because what one is really interested in is diagonalizing the Fock matrix, of which the kinetic energy is only one part.


Also, it is always possible to choose a basis where the kinetic energy matrix is diagonal, but again, this would not get you anywhere as in solving the HF problem, you will diagonalize the Fock Matrix which will surely put you in a basis in which the kinetic energy matrix is non-diagonal. What might be more interesting would be to take the actual molecular orbitals which are part of the solutions of the Roothan-Hall equations and re-compute the kinetic energy matrix over these orbitals. Looking at the coupling between these orbitals may be more easily interpretable for e.g. aromatic $\pi$-systems.


No comments:

Post a Comment

readings - Appending 内 to a company name is read ない or うち?

For example, if I say マイクロソフト内のパートナーシップは強いです, is the 内 here read as うち or ない? Answer 「内」 in the form: 「Proper Noun + 内」 is always read 「ない...