Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

the details derivation of the local current? #25

Open
qqgu opened this issue Jan 27, 2022 · 3 comments
Open

the details derivation of the local current? #25

qqgu opened this issue Jan 27, 2022 · 3 comments

Comments

@qqgu
Copy link

qqgu commented Jan 27, 2022

Hi, in the transport example, you show the local current calculation formula as,
image

can you tell me the details of derivation for this formula? in the orthogonal basis set, the local current should be H\rho , where is the -S\epsilon term com form? and how to calculate the \epsilon?

Thanks a lot.

@awvwgk
Copy link
Member

awvwgk commented Jan 27, 2022

in the orthogonal basis set, the local current should be H\rho , where is the -S\epsilon term com form

Probably from exactly this. The basis set is not orthogonal, therefore the overlap of the basis functions has to be accounted for.

@qqgu
Copy link
Author

qqgu commented Jan 27, 2022

in the orthogonal basis set, the local current should be H\rho , where is the -S\epsilon term com form

Probably from exactly this. The basis set is not orthogonal, therefore the overlap of the basis functions has to be accounted for.

Thanks for the immediate response.
Yes, the non-orthogonal basis set is the reason of the additional terms -s\epsilon.
But i just don't know the derivation of this. \epsilon is referred to as energy-weighted density matrix, what exactly does it look like? something like $\int E*\rho(E) dE$ ?

@aradi
Copy link
Member

aradi commented Jan 27, 2022

When calculated with eigenvalues and eigenvectors, the density matrix is $P_{\mu,\nu} = \sum_i f_i c_{i,\mu} c_{i,\nu}^*$ and the energy weighted density matrix is $P^\epsilon_{\mu,\nu} = \sum_i f_i \epsilon_i c_{i,\mu} c_{i,\nu}^*$, with $\epsilon_i$ being the eigenvalue of the eigenstate $i$.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

3 participants