From chemistry-request: at :ccl.net Thu Nov  4 15:15:38 2004
Received: from smtp-b.ua.edu (smtp-b.ua.edu [130.160.4.39])
	by server.ccl.net (8.12.8/8.12.8) with ESMTP id iA4KFWsE018438
	for <chemistry^at^ccl.net>; Thu, 4 Nov 2004 15:15:32 -0500
Received: from [130.160.44.65] ([130.160.44.65])
	by smtp-b.ua.edu (8.12.10/8.12.10) with ESMTP id iA4KSf04021582
	for <chemistry^at^ccl.net>; Thu, 4 Nov 2004 14:28:41 -0600 (CST)
Message-ID: <418A9090.9000500^at^bama.ua.edu>
Date: Thu, 04 Nov 2004 14:26:56 -0600
From: "Keith E. Gutowski" <gutow001^at^bama.ua.edu>
User-Agent: Mozilla/5.0 (X11; U; Linux i686; en-US; rv:1.7.3) Gecko/20040916
X-Accept-Language: en-us, en
MIME-Version: 1.0
To: chemistry^at^ccl.net
Subject: Volume keyword in Gaussian: problem
Content-Type: text/plain; charset=ISO-8859-1; format=flowed
Content-Transfer-Encoding: 7bit
X-Spam-Status: No, hits=4.5 required=7.5 tests=FROM_ENDS_IN_NUMS,NO_RDNS2 
	autolearn=no version=2.61
X-Spam-Checker-Version: SpamAssassin 2.61 (1.212.2.1-2003-12-09-exp) on 
	servernd.ccl.net

I am trying to compute ionic volumes using the 'volume=tight' keyword in 
Gaussian03.  In the manual, it states for a volume calc:*

'Tight*
Requests an increased density of points for more accurate integration. 
By default, the volume is computed to an accuracy of about 10%. Use of 
this option is recommended if the computed molecular volume is needed 
more quantitatively.'

Now, when I caculate the volume of PF6- using B3LYP, DZVP2 basis set 
(opt + volume=tight), I get a strange result.  Here is the volume output 
for 6 different runs of the same input file (in cubic cm per mol):

69.815
57.036
65.763
60.465
68.880
54.854

Now, I would hope that Gaussian03 could generate better reproducibility 
that this, but am I wrong?  I have been trying to track down the source 
of this, but I am a bit stumped.  The only this I can come up with has 
to do with the error on the integrated density.  For each run (PF6- or 
whatever molecule I do), the integrated density is about 4.23D+01, while 
the error on the integrated density is about -2.75D+01.   I thought the 
error on the integrated density should be around XXXD-06.  This is 
printed out in the output along with the volume output results.

Anyone have any idea what may be going on?

Keith