2015-11-24 17:01:09 +01:00
|
|
|
subroutine build_fock_tmp(fock_diag_tmp,det_ref,Nint)
|
|
|
|
|
use bitmasks
|
|
|
|
|
implicit none
|
|
|
|
|
BEGIN_DOC
|
|
|
|
|
! Build the diagonal of the Fock matrix corresponding to a generator
|
|
|
|
|
! determinant. F_00 is <i|H|i> = E0.
|
|
|
|
|
END_DOC
|
|
|
|
|
integer, intent(in) :: Nint
|
|
|
|
|
integer(bit_kind), intent(in) :: det_ref(Nint,2)
|
|
|
|
|
double precision, intent(out) :: fock_diag_tmp(2,mo_tot_num+1)
|
|
|
|
|
|
|
|
|
|
integer :: occ(Nint*bit_kind_size,2)
|
|
|
|
|
integer :: ne(2), i, j, ii, jj
|
|
|
|
|
double precision :: E0
|
|
|
|
|
|
|
|
|
|
! Compute Fock matrix diagonal elements
|
|
|
|
|
call bitstring_to_list_ab(det_ref,occ,Ne,Nint)
|
|
|
|
|
|
|
|
|
|
fock_diag_tmp = 0.d0
|
|
|
|
|
E0 = 0.d0
|
|
|
|
|
|
2016-12-29 22:00:41 +01:00
|
|
|
if (Ne(1) /= elec_alpha_num) then
|
|
|
|
|
print *, 'Error in build_fock_tmp (alpha)', Ne(1), Ne(2)
|
2018-07-26 11:31:36 +02:00
|
|
|
call debug_det(det_ref,N_int)
|
|
|
|
|
print *, occ(:,1)
|
|
|
|
|
print *, occ(:,2)
|
|
|
|
|
do i=1,10000
|
|
|
|
|
occ(i,1) = fock_diag_tmp(1,mo_tot_num+i) !traceback
|
|
|
|
|
enddo
|
2016-12-29 22:00:41 +01:00
|
|
|
stop -1
|
|
|
|
|
endif
|
|
|
|
|
if (Ne(2) /= elec_beta_num) then
|
|
|
|
|
print *, 'Error in build_fock_tmp (beta)', Ne(1), Ne(2)
|
2018-07-26 11:31:36 +02:00
|
|
|
do i=1,10000
|
|
|
|
|
occ(i,1) = fock_diag_tmp(1,mo_tot_num+i) !traceback
|
|
|
|
|
enddo
|
2016-12-29 22:00:41 +01:00
|
|
|
stop -1
|
|
|
|
|
endif
|
|
|
|
|
|
2015-11-24 17:01:09 +01:00
|
|
|
! Occupied MOs
|
|
|
|
|
do ii=1,elec_alpha_num
|
|
|
|
|
i = occ(ii,1)
|
|
|
|
|
fock_diag_tmp(1,i) = fock_diag_tmp(1,i) + mo_mono_elec_integral(i,i)
|
|
|
|
|
E0 = E0 + mo_mono_elec_integral(i,i)
|
|
|
|
|
do jj=1,elec_alpha_num
|
|
|
|
|
j = occ(jj,1)
|
|
|
|
|
if (i==j) cycle
|
|
|
|
|
fock_diag_tmp(1,i) = fock_diag_tmp(1,i) + mo_bielec_integral_jj_anti(i,j)
|
|
|
|
|
E0 = E0 + 0.5d0*mo_bielec_integral_jj_anti(i,j)
|
|
|
|
|
enddo
|
|
|
|
|
do jj=1,elec_beta_num
|
|
|
|
|
j = occ(jj,2)
|
|
|
|
|
fock_diag_tmp(1,i) = fock_diag_tmp(1,i) + mo_bielec_integral_jj(i,j)
|
|
|
|
|
E0 = E0 + mo_bielec_integral_jj(i,j)
|
|
|
|
|
enddo
|
|
|
|
|
enddo
|
|
|
|
|
do ii=1,elec_beta_num
|
|
|
|
|
i = occ(ii,2)
|
|
|
|
|
fock_diag_tmp(2,i) = fock_diag_tmp(2,i) + mo_mono_elec_integral(i,i)
|
|
|
|
|
E0 = E0 + mo_mono_elec_integral(i,i)
|
|
|
|
|
do jj=1,elec_beta_num
|
|
|
|
|
j = occ(jj,2)
|
|
|
|
|
if (i==j) cycle
|
|
|
|
|
fock_diag_tmp(2,i) = fock_diag_tmp(2,i) + mo_bielec_integral_jj_anti(i,j)
|
|
|
|
|
E0 = E0 + 0.5d0*mo_bielec_integral_jj_anti(i,j)
|
|
|
|
|
enddo
|
|
|
|
|
do jj=1,elec_alpha_num
|
|
|
|
|
j = occ(jj,1)
|
|
|
|
|
fock_diag_tmp(2,i) = fock_diag_tmp(2,i) + mo_bielec_integral_jj(i,j)
|
|
|
|
|
enddo
|
|
|
|
|
enddo
|
|
|
|
|
|
|
|
|
|
! Virtual MOs
|
|
|
|
|
do i=1,mo_tot_num
|
|
|
|
|
if (fock_diag_tmp(1,i) /= 0.d0) cycle
|
|
|
|
|
fock_diag_tmp(1,i) = fock_diag_tmp(1,i) + mo_mono_elec_integral(i,i)
|
|
|
|
|
do jj=1,elec_alpha_num
|
|
|
|
|
j = occ(jj,1)
|
|
|
|
|
fock_diag_tmp(1,i) = fock_diag_tmp(1,i) + mo_bielec_integral_jj_anti(i,j)
|
|
|
|
|
enddo
|
|
|
|
|
do jj=1,elec_beta_num
|
|
|
|
|
j = occ(jj,2)
|
|
|
|
|
fock_diag_tmp(1,i) = fock_diag_tmp(1,i) + mo_bielec_integral_jj(i,j)
|
|
|
|
|
enddo
|
|
|
|
|
enddo
|
|
|
|
|
do i=1,mo_tot_num
|
|
|
|
|
if (fock_diag_tmp(2,i) /= 0.d0) cycle
|
|
|
|
|
fock_diag_tmp(2,i) = fock_diag_tmp(2,i) + mo_mono_elec_integral(i,i)
|
|
|
|
|
do jj=1,elec_beta_num
|
|
|
|
|
j = occ(jj,2)
|
|
|
|
|
fock_diag_tmp(2,i) = fock_diag_tmp(2,i) + mo_bielec_integral_jj_anti(i,j)
|
|
|
|
|
enddo
|
|
|
|
|
do jj=1,elec_alpha_num
|
|
|
|
|
j = occ(jj,1)
|
|
|
|
|
fock_diag_tmp(2,i) = fock_diag_tmp(2,i) + mo_bielec_integral_jj(i,j)
|
|
|
|
|
enddo
|
|
|
|
|
enddo
|
|
|
|
|
|
|
|
|
|
fock_diag_tmp(1,mo_tot_num+1) = E0
|
|
|
|
|
fock_diag_tmp(2,mo_tot_num+1) = E0
|
|
|
|
|
|
|
|
|
|
end
|
