Lower energy with larger basis set


Just Got Here
Hello,

I am running all-electron DFT calculations for single atoms and I am interested in the total energy. As DFT is variational one expects the larger basis sets to yield lower total energies, which is also the case in most of my calculations, except for the aug-cc-pv6z basis which yields a higher total energy ( -128.226194282698 a.u.) than the aug-cc-pv5z basis (-128.227026486337 a.u.). Can there be an error in the implementation of the basis, are some parameters like the grid for the integration not choseb appropriately or can this really happen from the way the aug-cc-pc6z basis is set up ?

The input file I use is

title "Ne"
echo
geometry units angstroms
Ne      0.0  0.0  0.0
end
basis
Ne library aug-cc-pv6Z 
end
dft
grid xfine
xc slater perdew81 
end
TASK DFT


Thanks for your answers

Volker

Forum Vet
Linear dependencies
Hi Volker,

What you are seeing is the effect of removing linear dependencies in the basis set. Only two are found for the aug-cc-pv5z, while 33 are found for the aug-cc-pv6z. We have a pretty aggressive threshold for removing linear dependencies from the basis set of 1.0D-5 in the overlap matrix.

Keywords to adjust the linear dependencies are:

set lindep:tol 1.0d-5 (sets linear dependency cutoff at 1.0d-5)

set lindep:n_dep 0 (number specifies how many linear dependent functions should be removed. By setting this to zero all functions are included)

I would recommend using "set lindep:tol 1.0d-6". This will get you the behavior you want. You could even use 1.0d-7.

While "set lindep:n_dep 0" is an option it will not work here. There are such strong linear dependencies in the basis set that it the DFT energy will simply not converge.

Bert


Forum >> NWChem's corner >> General Topics