From chemistry-request-!at!-server.ccl.net  Tue Jun 27 01:31:09 2000
Received: from uvaix7e1.comp.UVic.CA (root:~at~:uvaix7e1.comp.UVic.CA [142.104.5.107])
	by server.ccl.net (8.8.7/8.8.7) with ESMTP id BAA14935
	for <CHEMISTRY {*at*} ccl.net>; Tue, 27 Jun 2000 01:31:08 -0400
Received: from p17-111.dialup.UVic.CA (p17-111.dialup.UVic.CA [142.104.17.111])
	by uvaix7e1.comp.UVic.CA (8.9.3/8.9.3) with SMTP id WAA117580
	for <CHEMISTRY %-% at %-% ccl.net>; Mon, 26 Jun 2000 22:31:00 -0700
From: Roy Jensen <royj;at;uvic.ca>
To: CHEMISTRY _-at-_)ccl.net
Subject: Help on transition metal--chlorine optimization
Date: Mon, 26 Jun 2000 22:30:57 -0700
Message-ID: <vvegls84taoqsh17ip1snv00i6ugf9f6j0 "-at-" 4ax.com>
X-Mailer: Forte Agent 1.8/32.548
MIME-Version: 1.0
Content-Type: text/plain; charset=us-ascii
Content-Transfer-Encoding: 8bit
X-MIME-Autoconverted: from quoted-printable to 8bit by server.ccl.net id BAA14937

All

I posted a message to the CCL a few weeks ago; a summary was just
posted. DFT was the choice for geometry optimizations on transition
metal--chlorine clusters (neutral and ionic). 

I am using Gaussian 94W and 98W however I cannot get the simplest job
to run. I have tried several methods and various basis sets. Almost
all fail to converge on the first step. A sample input and output file
are given below. 

Thanks,
Roy Jensen


INPUT FILE
=========
# UB3LYP LANL2DZ opt test
 
Optimized geometry for FeCl.
 
0 2
Fe
Cl 1 R1

R1=2.19


OUTPUT FILE
===========
...
 <S**2> of initial guess=  .7500
 Requested convergence on RMS density matrix=1.00E-08 within  64
cycles.
 Requested convergence on MAX density matrix=1.00E-06.
 Keep R1 and R2 integrals in memory in canonical form, NReq=
727936.
 Integral accuracy reduced to 1.0E-05 until final iterations.
 Problem detected with inexpensive integrals.
 Switching to full accuracy and repeating last cycle.
 >>>>>>>>>> Convergence criterion not met.
 SCF Done:  E(UB+HF-LYP) =  -137.386151674     A.U. after   65 cycles
             Convg  =     .4226E-01             -V/T =  2.6548
             S**2   =    .8263
 Annihilation of the first spin contaminant:
 S**2 before annihilation      .8263,   after      .7500
 Convergence failure -- run terminated.
 Error termination via Lnk1e in D:\G94W\l502.exe.
 Job cpu time:  0 days  0 hours 11 minutes 38.0 seconds.
 File lengths (MBytes):  RWF=    5 Int=    0 D2E=    0 Chk=    1 Scr=
1


(increasing 'maxcyc' and using 'novaracc' does not help)