2023-09-05 17:23:40 +02:00
|
|
|
subroutine ROHF_fock_matrix(nBas,nOa,nOb,S,c,Fa,Fb,F)
|
|
|
|
|
|
|
|
! Construct the ROHF Fock matrix
|
|
|
|
|
|
|
|
implicit none
|
|
|
|
include 'parameters.h'
|
|
|
|
|
|
|
|
! Input variables
|
|
|
|
|
|
|
|
integer,intent(in) :: nBas
|
|
|
|
integer,intent(in) :: nOa
|
|
|
|
integer,intent(in) :: nOb
|
|
|
|
|
|
|
|
double precision,intent(in) :: S(nBas,nBas)
|
|
|
|
double precision,intent(in) :: c(nBas,nBas)
|
|
|
|
double precision,intent(inout):: Fa(nBas,nBas)
|
|
|
|
double precision,intent(inout):: Fb(nBas,nBas)
|
|
|
|
|
|
|
|
! Local variables
|
|
|
|
|
|
|
|
double precision :: aC,bC
|
|
|
|
double precision :: aO,bO
|
|
|
|
double precision :: aV,bV
|
|
|
|
|
|
|
|
integer :: nC
|
|
|
|
integer :: nO
|
|
|
|
integer :: nV
|
|
|
|
|
|
|
|
! Output variables
|
|
|
|
|
|
|
|
double precision,intent(out) :: F(nBas,nBas)
|
|
|
|
|
|
|
|
! Roothan canonicalization parameters
|
|
|
|
|
|
|
|
aC = -0.5d0
|
|
|
|
bC = +1.5d0
|
|
|
|
|
|
|
|
aO = +0.5d0
|
|
|
|
bO = +0.5d0
|
|
|
|
|
|
|
|
aV = +1.5d0
|
|
|
|
bV = -0.5d0
|
|
|
|
|
|
|
|
! Number of closed, open, and virtual orbitals
|
|
|
|
|
|
|
|
nC = min(nOa,nOb)
|
|
|
|
nO = abs(nOa - nOb)
|
|
|
|
nV = nBas - nC - nO
|
|
|
|
|
|
|
|
! Block-by-block Fock matrix
|
|
|
|
|
2023-09-06 09:50:29 +02:00
|
|
|
! call AOtoMO_transform(nBas,c,Fa)
|
|
|
|
! call AOtoMO_transform(nBas,c,Fb)
|
2023-09-05 17:23:40 +02:00
|
|
|
|
|
|
|
F(1:nC, 1:nC ) = aC*Fa(1:nC, 1:nC ) + bC*Fb(1:nC, 1:nC )
|
|
|
|
F(1:nC, nC+1:nC+nO ) = Fb(1:nC, nC+1:nC+nO )
|
|
|
|
F(1:nC,nO+nC+1:nC+nO+nV) = 0.5d0*Fa(1:nC,nO+nC+1:nC+nO+nV) + 0.5d0*Fb(1:nC,nO+nC+1:nC+nO+nV)
|
|
|
|
|
|
|
|
F(nC+1:nO+nC, 1:nC ) = Fb(nC+1:nO+nC, 1:nC )
|
2023-09-06 09:50:29 +02:00
|
|
|
F(nC+1:nO+nC, nC+1:nC+nO ) = aO*Fa(nC+1:nO+nC, nC+1:nC+nO ) + bO*Fb(nC+1:nO+nC,nC+1:nC+nO)
|
|
|
|
F(nC+1:nO+nC,nO+nC+1:nC+nO+nV) = Fa(nC+1:nO+nC,nO+nC+1:nC+nO+nV)
|
2023-09-05 17:23:40 +02:00
|
|
|
|
|
|
|
F(nO+nC+1:nC+nO+nV, 1:nC ) = 0.5d0*Fa(nO+nC+1:nC+nO+nV, 1:nC ) + 0.5d0*Fb(nO+nC+1:nC+nO+nV, 1:nC )
|
|
|
|
F(nO+nC+1:nC+nO+nV, nC+1:nC+nO ) = Fa(nO+nC+1:nC+nO+nV, nC+1:nC+nO )
|
|
|
|
F(nO+nC+1:nC+nO+nV,nO+nC+1:nC+nO+nV) = aV*Fa(nO+nC+1:nC+nO+nV,nO+nC+1:nC+nO+nV) + bV*Fb(nO+nC+1:nC+nO+nV,nO+nC+1:nC+nO+nV)
|
|
|
|
|
2023-09-06 09:50:29 +02:00
|
|
|
! call MOtoAO_transform(nBas,S,c,F)
|
|
|
|
! call MOtoAO_transform(nBas,S,c,Fa)
|
|
|
|
! call MOtoAO_transform(nBas,S,c,Fb)
|
2023-09-05 17:23:40 +02:00
|
|
|
|
|
|
|
end subroutine
|