# Understanding Gaussian cube output of total density using DPLOT

Viewed 2678 times, With a total of 1 Posts

 Kpelzer Member Profile Send PM
 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 Profile Send PM
 Forum Vet Threads 9 Posts 1492
 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.

 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