LCPQ/quantum_package
10
0
Fork 0
mirror of https://github.com/LCPQ/quantum_package synced 2025-01-08 20:33:26 +01:00
Code Issues Releases Wiki Activity

Added level shift in Hartree-Fock

Browse Source
This commit is contained in:
Anthony Scemama 2015-11-16 19:00:49 +01:00
parent 6a640567e8
commit 0c84015982
2 changed files with 15 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
plugins/Hartree_Fock/EZFIO.cfg
View File

@ -10,6 +10,12 @@ doc: Maximum number of SCF iterations
interface: ezfio,provider,ocaml interface: ezfio,provider,ocaml
default: 200 default: 200
[level_shift]
type: Positive_float
doc: Energy shift on the virtual MOs
interface: ezfio,provider,ocaml
default: 0.0
[mo_guess_type] [mo_guess_type]
type: MO_guess type: MO_guess
doc: Initial MO guess. Can be [ Huckel | HCore ] doc: Initial MO guess. Can be [ Huckel | HCore ]

11
plugins/Hartree_Fock/Fock_matrix.irp.f
View File

@ -73,8 +73,13 @@
enddo enddo
endif endif
! Introduce level shift here
do i = elec_alpha_num+1, mo_tot_num
Fock_matrix_mo(i,i) += level_shift
enddo
do i = 1, mo_tot_num do i = 1, mo_tot_num
Fock_matrix_diag_mo(i) = Fock_matrix_mo(i,i) Fock_matrix_diag_mo(i) = Fock_matrix_mo(i,i)
enddo enddo
END_PROVIDER END_PROVIDER
@ -315,7 +320,9 @@ BEGIN_PROVIDER [ double precision, Fock_matrix_ao, (ao_num_align, ao_num) ]
! Fock matrix in AO basis set ! Fock matrix in AO basis set
END_DOC END_DOC
if (elec_alpha_num == elec_beta_num) then if ( (elec_alpha_num == elec_beta_num).and. &
(level_shift == 0.) ) &
then
integer :: i,j integer :: i,j
do j=1,ao_num do j=1,ao_num
!DIR$ VECTOR ALIGNED !DIR$ VECTOR ALIGNED