10
1
mirror of https://github.com/pfloos/quack synced 2024-12-25 13:53:41 +01:00
QuAcK/src/MBPT/print_G0T0.f90

63 lines
2.5 KiB
Fortran
Raw Normal View History

2021-10-16 15:34:34 +02:00
subroutine print_G0T0(nBas,nO,eHF,ENuc,ERHF,SigT,Z,eGT,EcRPA)
2019-10-16 18:14:47 +02:00
! Print one-electron energies and other stuff for G0T0
implicit none
include 'parameters.h'
integer,intent(in) :: nBas,nO
double precision,intent(in) :: ENuc
2019-10-20 08:18:58 +02:00
double precision,intent(in) :: ERHF
double precision,intent(in) :: EcRPA(nspin)
2021-10-16 15:34:34 +02:00
double precision,intent(in) :: eHF(nBas)
2019-10-28 16:34:09 +01:00
double precision,intent(in) :: SigT(nBas)
double precision,intent(in) :: Z(nBas)
2021-10-16 15:34:34 +02:00
double precision,intent(in) :: eGT(nBas)
2019-10-16 18:14:47 +02:00
2019-10-28 16:34:09 +01:00
integer :: p,HOMO,LUMO
2019-10-16 18:14:47 +02:00
double precision :: Gap
! HOMO and LUMO
HOMO = nO
LUMO = HOMO + 1
2020-03-21 16:31:39 +01:00
if(nBas >= LUMO) then
2021-10-16 15:34:34 +02:00
Gap = eGT(LUMO) - eGT(HOMO)
2020-03-21 16:31:39 +01:00
else
Gap = 0d0
end if
2019-10-16 18:14:47 +02:00
! Dump results
write(*,*)'-------------------------------------------------------------------------------'
write(*,*)' One-shot G0T0 calculation (T-matrix self-energy) '
write(*,*)'-------------------------------------------------------------------------------'
write(*,'(1X,A1,1X,A3,1X,A1,1X,A15,1X,A1,1X,A15,1X,A1,1X,A15,1X,A1,1X,A15,1X,A1,1X)') &
2019-10-17 22:07:41 +02:00
'|','#','|','e_HF (eV)','|','Sigma_T (eV)','|','Z','|','e_QP (eV)','|'
2019-10-16 18:14:47 +02:00
write(*,*)'-------------------------------------------------------------------------------'
2019-10-28 16:34:09 +01:00
do p=1,nBas
2019-10-16 18:14:47 +02:00
write(*,'(1X,A1,1X,I3,1X,A1,1X,F15.6,1X,A1,1X,F15.6,1X,A1,1X,F15.6,1X,A1,1X,F15.6,1X,A1,1X)') &
2021-10-16 15:34:34 +02:00
'|',p,'|',eHF(p)*HaToeV,'|',SigT(p)*HaToeV,'|',Z(p),'|',eGT(p)*HaToeV,'|'
2019-10-16 18:14:47 +02:00
enddo
write(*,*)'-------------------------------------------------------------------------------'
2021-10-16 15:34:34 +02:00
write(*,'(2X,A50,F15.6,A3)') 'G0T0 HOMO energy (eV) =',eGT(HOMO)*HaToeV,' eV'
write(*,'(2X,A50,F15.6,A3)') 'G0T0 LUMO energy (eV) =',eGT(LUMO)*HaToeV,' eV'
write(*,'(2X,A50,F15.6,A3)') 'G0T0 HOMO-LUMO gap (eV) =',Gap*HaToeV,' eV'
write(*,*)'-------------------------------------------------------------------------------'
write(*,'(2X,A50,F20.10,A3)') 'Tr@RPA@G0T0 correlation energy (singlet) =',EcRPA(1),' au'
write(*,'(2X,A50,F20.10,A3)') 'Tr@RPA@G0T0 correlation energy (triplet) =',EcRPA(2),' au'
write(*,'(2X,A50,F20.10,A3)') 'Tr@RPA@G0T0 correlation energy =',EcRPA(1) + EcRPA(2),' au'
write(*,'(2X,A50,F20.10,A3)') 'Tr@RPA@G0T0 total energy =',ENuc + ERHF + EcRPA(1) + EcRPA(2),' au'
2019-10-16 18:14:47 +02:00
write(*,*)'-------------------------------------------------------------------------------'
write(*,*)
end subroutine print_G0T0