mirror of
https://github.com/QuantumPackage/qp2.git
synced 2024-12-22 03:23:29 +01:00
Improved huckel guess
This commit is contained in:
parent
21d42dc7d7
commit
008d1b9e8f
@ -15,7 +15,7 @@ function run() {
|
|||||||
qp run scf
|
qp run scf
|
||||||
# qp set_frozen_core
|
# qp set_frozen_core
|
||||||
energy="$(ezfio get hartree_fock energy)"
|
energy="$(ezfio get hartree_fock energy)"
|
||||||
eq $2 $energy $thresh
|
eq $energy $2 $thresh
|
||||||
}
|
}
|
||||||
|
|
||||||
|
|
||||||
|
@ -15,7 +15,8 @@ subroutine huckel_guess
|
|||||||
A = 0.d0
|
A = 0.d0
|
||||||
do j=1,ao_num
|
do j=1,ao_num
|
||||||
do i=1,ao_num
|
do i=1,ao_num
|
||||||
A(i,j) = c * ao_overlap(i,j) * (ao_one_e_integrals_diag(i) + ao_one_e_integrals_diag(j))
|
A(i,j) = c * ao_overlap(i,j)/(ao_overlap(i,i)*ao_overlap(j,j)) * &
|
||||||
|
(ao_one_e_integrals_diag(i) + ao_one_e_integrals_diag(j))
|
||||||
enddo
|
enddo
|
||||||
A(j,j) = ao_one_e_integrals_diag(j) + ao_two_e_integral_alpha(j,j)
|
A(j,j) = ao_one_e_integrals_diag(j) + ao_two_e_integral_alpha(j,j)
|
||||||
enddo
|
enddo
|
||||||
|
Loading…
Reference in New Issue
Block a user