From chemistry-request %-% at %-% server.ccl.net Mon Mar 18 08:00:31 2002 Received: from MAIL.AD.Berry.edu ([66.20.28.66]) by server.ccl.net (8.11.6/8.11.0) with ESMTP id g2ID0Vp31641 for ; Mon, 18 Mar 2002 08:00:31 -0500 Received: from [10.16.232.24] ([10.16.232.24]) by MAIL.AD.Berry.edu with Microsoft SMTPSVC(5.0.2195.4453); Mon, 18 Mar 2002 07:22:27 -0500 User-Agent: Microsoft-Outlook-Express-Macintosh-Edition/5.02.2022 Date: Mon, 18 Mar 2002 07:21:01 -0500 Subject: SUMMARY: A couple of questions about AIM analysis From: Gary Breton To: Message-ID: Mime-version: 1.0 Content-type: text/plain; charset="US-ASCII" Content-transfer-encoding: 7bit X-OriginalArrivalTime: 18 Mar 2002 12:22:27.0420 (UTC) FILETIME=[911E59C0:01C1CE77] Hello all, I thank everyone for answering my question on AIM analysis. I was asked by some to post a summary, and so I am doing so. It would appear that the overall concensus is that the Gaussian implementation of AIM is not the most robust and that more specialized software is more appropriate (i.e., MORPHY, AIM2000, etc.). ************************* First, my original question: I have implemented AIM analysis from with Gaussian98W successfully many times. The keywords I use are: SCF(conver=8) density=current AIM=bondorders. The level of theory in all cases has been RB3LYP/6-311+G(2d,p) which is the same level of theory as the other calculations involving these molecules (e.g., minimizations, frequency calculations, NBO analysis). However,when attempting the same calculation on one particular structure, the calculation fails. I am given the information below: WARNING: RMS ERROR HAS INCREASED WARNING: RMS ERROR HAS INCREASED WARNING: RMS ERROR HAS INCREASED WARNING: RMS ERROR HAS INCREASED NEWTON STEP FAILED FOR SURFACE SHEET 16 Error termination via Lnk1e in C:\G98W\l609.exe. Job cpu time: 0 days 4 hours 16 minutes 8.0 seconds. File lengths (MBytes): RWF= 48 Int= 0 D2E= 0 Chk= 2 Scr= 1 The Gaussian User's Reference suggests that the AIM code will fail at times and cites the particular instance when the message NEWTON-RAPHSON STEP FAILED is due to a molecule with rings that has an unusual topology (which may well be the case). Is there any cure for this problem? Am I able to tweak something to get the calculation to go to completion? ******************** Don't know how to force gaussian to work round difficult topologies, but our in house code (MORPHY) is available for a modest fee and is it is easy adapt the internal settings to get round most problems [it will also look at the topology of the laplacian, and elf, and soon other fields like lol]. If you want any more info let me know. Richard Baders code (AIMPAC) is freely available and should be flexible enough to do what you want but I don't have any direct experience of it. regards noj Malcom ******************* We have done quite a bit of AIM recently but not using Gaussian. There have been several comments about the G98 AIM package that suggest that it is not robust, and will fail in cases like yours where there is extreme curvature in the atomic surfaces. I have heard of no "fix" for these cases. I believe that the G98 AIM package is essentially a re-packaging of the original AIMPAC bundle of Bader et al. We have found that the AIM2000 package has given us very good results, although the MORPHY program is also claimed to be very good and may in fact be superior if its long-awaited next version ever comes out. In several cases, atomic integrations attempted with the PROAIMV program from the AIMPAC bundle failed no matter what parameters we tried, but integration with AIM2000 succeeded. We have also found it much easier to enhance the accuracy of integrals with AIM2000. Good Luck. Phil Hultin ***************** The Gaussian implementation of AIM is not reknowned for its robustness. Their love affair ended a long time ago and the aim in gaussian is an old and not maintained version. I suggest using aim2000 or aimpac or morphy from a gaussian output wfn file. Yours, Alexandre HOCQUET *************** There is no known cure for this error. The program is attempting to fit a surface so it can do integration within the various regions. There are a number of topological features which cause it problems. The code is doing a first order search and this is obviously not high enough order but no one has taken on the task of implementing a higher order search, linear search etc, similar to what is done for difficult geometry optimization. There are a couple of bond order estimates in the NBO package. If you use # ... POP=NBORead and then after the structure add lines, $NBO BNDIDX NLMO $END you will see these printed out. If you need an estimate of the bond order better than something like Mulliken that is my recommendation. Douglas J. Fox Technical Support Gaussian, Inc. ******************* Dear Professor Gary Breton, Gaussian has a multiple of optimization and minimization methods. I think you are using the most default method, but I think there are other methods available, like QC, GDIIS, NODIIS, etc. I am not exactly sure where the error occurred because you cut off the output a little too short, but there are other methods. I am curious. Why do you use SCF(COVER=8)? The default I think is 7. Is there are a need for the extra accuracy in the SCF cycle for AIM calculations? Yours, Ha Yeon Cheong