1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
|
! Copyright (C) 2020 J. K. Dewhurst and S. Sharma.
! This file is distributed under the terms of the GNU General Public License.
! See the file COPYING for license details.
subroutine energytd
use modmain
use modtddft
implicit none
! local variables
integer is,ias
real(8) ca,engya,t1,t2
! external functions
real(8), external :: rfinp
! Coulomb potential energy
engyvcl=rfinp(rhomt,rhoir,vclmt,vclir)
! Madelung term
engymad=0.d0
do ias=1,natmtot
is=idxis(ias)
engymad=engymad+0.5d0*spzn(is)*(vclmt(1,ias)-vcln(1,is))*y00
end do
! exchange and correlation energy
engyx=rfinp(rhomt,rhoir,exmt,exir)
engyc=rfinp(rhomt,rhoir,ecmt,ecir)
! vector potential contributions to energy
ca=-1.d0/solsc
! coupling term -1/c A(t)⋅J(t)
engya=ca*dot_product(afieldt(1:3,itimes),jtot(1:3))
! constant term 1/2c² A²(t)Q where Q is the 'dynamic' charge
ca=0.5d0/solsc**2
t1=sum(afieldt(1:3,itimes)**2)
t2=sum(chgstot(1:3))/3.d0
engya=engya+ca*t1*(chgtot-t2)
! total energy
engytot=engykn+0.5d0*engyvcl+engymad+engyx+engyc+engya
end subroutine
|
