TREX/qmckl
1
0
Fork 0
mirror of https://github.com/TREX-CoE/qmckl.git synced 2025-01-03 18:16:28 +01:00
Code Issues Projects Releases Wiki Activity

Added doc for the three energies. #41

Browse Source
This commit is contained in:
v1j4y 2021-10-11 18:48:48 +02:00
parent 1be1317210
commit edc6e0a4f6
1 changed files with 41 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

41
org/qmckl_local_energy.org
View File

@ -185,6 +185,19 @@ typedef struct qmckl_local_energy_struct {
:FRetType: qmckl_exit_code
:END:
Where the kinetic energy is given as:
\[
KE = -\frac{1}{2}\frac{\bigtriangleup \Psi}{\Psi}
\]
The laplacian of the wavefunction in the single-determinant
case is given as follows:
\[
\frac{\bigtriangleup \Psi(r)}{\Psi(r)} = \sum_{j=1}^{N_e} \bigtriangleup \Phi_j(r_i) D_{ji}^{-1}(r)
\]
*** Get
#+begin_src c :comments org :tangle (eval h_func) :noweb yes
@ -513,6 +526,27 @@ end function qmckl_compute_kinetic_energy_f
:FRetType: qmckl_exit_code
:END:
The potential energy is the sum of all the following terms
\[
PE = \mathcal{V}_{ee} + \mathcal{V}_{en} + \mathcal{V}_{nn}
\]
The potential for is calculated as the sum of single electron
contributions.
\[
\mathcal{V}_{ee} = \sum_{i=1}^{N_e}\sum_{j<i} \frac{1}{r_{ij}}
\]
\[
\mathcal{V}_{en} = - \sum_{i=1}^{N_e}\sum_{A=1}^{N_n}\frac{Z_A}{r_{iA}}
\]
\[
\mathcal{V}_{nn} = \sum_{A=1}^{N_n}\sum_{B<A}\frac{Z_A Z_B}{r_{AB}}
\]
*** Get
#+begin_src c :comments org :tangle (eval h_func) :noweb yes
@ -756,6 +790,13 @@ end function qmckl_compute_potential_energy_f
:FRetType: qmckl_exit_code
:END:
The local energy is the sum of kinetic and potential energies.
\[
E_L = KE + PE
\]
*** Get
#+begin_src c :comments org :tangle (eval h_func) :noweb yes