Understanding Gaussian cube output of total density using DPLOT

From NWChem

Viewed 2564 times, With a total of 1 Posts
Jump to: navigation, search

Clicked A Few Times
Threads 9
Posts 19
I am a little puzzled by the cube files that I get when I use the following command to produce a Gaussian cube file with the total electron density:

dplot
gaussian
title monomer_12ang_anion_charge_density33
output monomer_12ang_anion_charge_density33.cube
vectors C70.movecs
spin total
LimitXYZ
-7.5 7.5 150
-7.5 7.5 150
-16.5 16.5 330
end
task dplot

I've been using some MATLAB code to process the cube files and decided to sum all data points in one of these density files to confirm that it adds up to the number of electrons. It does not; for a system of 421 electrons the sum of all the data points in the cube file is about 58,000. If I use the same value of the range-separation parameter (I use range-separated functionals), the sum of the numbers in the cube file is always about 58,000, with variations of 0.1-0.2% from this number. If I change the value of the range-separation parameter, the sum changes to about 61,000.

I seem to be missing something, I would have thought that a charge density file should sum to the number of electrons. Could anyone explain to me what DPLOT is doing here?

Thanks in advance for any advice.

Forum Vet
Threads 9
Posts 1462
Integration volume
Kenley,
Does your MATLAB code include the integration volume?

In your case, the integration volume is

(15./150.)^2*(33./330.)/(.529177)^3=.0067483

If you multiply .0067483*60,000 you get 405, that is not so far off your 421 number of electrons.


Forum >> NWChem's corner >> General Topics



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


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