From NWChem
Viewed 1186 times, With a total of 1 Posts
|
|
|
Clicked A Few Times
Threads 9
Posts 19
|
|
| 2:24:47 PM PDT - Wed, Oct 16th 2013 |
|
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.
|
|
|
-
Edoapra Forum:Admin, Forum:Mod, bureaucrat, sysop
|
|
|
Forum Vet
Threads 3
Posts 855
|
|
| 3:21:23 PM PDT - Thu, Oct 17th 2013 |
|
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.
|
|
|
AWC's:
2.5.10 MediaWiki - Stand Alone Forum Extension
Forum theme style by: AWC
