Ni ANO-RCC

From NWChem

Viewed 1005 times, With a total of 4 Posts
Jump to: navigation, search

Just Got Here
Threads 1
Posts 4
A strange HF energy of -1357 appears for Nickel atom with ANO-RCC, while it should be -1520.

Forum Vet
Threads 7
Posts 1356
Just got an energy of -1518.8

Total DFT energy = -1518.847822293797

This is the input file I used
start ni_ano
geometry 
ni 0 0 0
symmetry oh
end
memory stack 900 mb heap 200 mb global 700 mb noverify
BASIS "ao basis" spherical noprint
#BASIS SET: (21s,15p,10d,6f,4g,2h) -> [10s,9p,8d,6f,4g,2h]
Ni    S
....
end
dft
xc hfexch
end
relativistic
 douglas-kroll on dkh
 clight 137.036d0
end
task dft

Just Got Here
Threads 1
Posts 4
a singlet d10
That's my output: http://www.primefan.ru/stuff/chem/ni.txt

Just Got Here
Threads 1
Posts 4
the bug is located
It appears that the bug goes from Gaussian's routine for basis orthogonalization. NoBasisTransform option solves the problem, but significantly increases the job time:
http://www.primefan.ru/stuff/chem/ni_nobasistransform.txt

Surprisingly, for gold all seems to be OK. Orthogonalization routine leaves only 168 of 177 basis functions as it does for nickel, but the total energy do not change.

So I think some kind of linear dependency problems do occur in the ANO-RCC basis set.

Just Got Here
Threads 1
Posts 4
Quote:PrimeFan Sep 11th 10:05 pm
So I think some kind of linear dependency problems do occur in the ANO-RCC basis set.
I was wrong. Rev.D01 does not reproduce the bug. All is OK with the basis set. Sorry for the inconvenience.


Forum >> BSE: Basis Set Exchange >> General BSE Topics



Who's here now Members 0 Guests 1 Bots/Crawler 0


AWC's: 2.5.10 MediaWiki - Stand Alone Forum Extension
Forum theme style by: AWC