What is the default radial grid for DFT in DALTON?

Find answers or ask questions regarding Dalton calculations.
Please upload an output file showing the problem, if applicable.
(It is not necessary to upload input files, they can be found in the output file.)

Post Reply
rob
Posts: 27
Joined: 15 Oct 2014, 13:43
First name(s): Roberto
Last name(s): Di Remigio
Affiliation: CTCC
Country: Norway

What is the default radial grid for DFT in DALTON?

Post by rob » 21 Jul 2015, 09:10

Hi all,
I've been trying to compare some PSI4 DFT calculations with DALTON results.
From the DALTON manual is extremely clear how to modify all the radial and angular grid parameters (page 272, under the "*DFT INPUT" section)
However, the default radial grid and default partitioning scheme are not specified. The information on the grid generation strategy (both radial and angular)
is also nowhere to be found within the output file.
So, as from the subject: What is the default radial grid for DFT in DALTON?
Also, how many radial points are used in the radial grid when .RADINT has its default value of 1.0e-13?
Thanks a lot in advance!

bast
Posts: 1210
Joined: 26 Aug 2013, 13:22
First name(s): Radovan
Last name(s): Bast
Affiliation: none
Country: Germany

Re: What is the default radial grid for DFT in DALTON?

Post by bast » 21 Jul 2015, 09:23

dear Roberto,

this is the radial grid in DIRAC and AFAIK also the default
in Dalton (i think there are several options within Dalton):
http://dx.doi.org/10.1007/s002140100263

the radial grid is basis dependent so the number of points not only
depends on the parameter but also on the basis set range and angular momentum.
in other words if you add steep or diffuse functions you will get more points.

the partitioning is AFAIK this one: http://dx.doi.org/10.1063/1.454033

i have implemented a standalone lib for numerical grids:
https://github.com/bast/numgrid
i am not meaning to place a shameless plug here but it may be clearer
from that code what is going on since it is isolated.
it uses the above mentioned radial and partitioning schemes.
the other thing that you could do is to plug it into Psi4 and hopefully
get exactly the Dalton grid (the lib has a Python interface).

good luck!
radovan

Post Reply

Who is online

Users browsing this forum: No registered users and 9 guests