From chemistry-request@server.ccl.net  Fri May 19 19:08:23 2000
Received: from soinplx1.solutia.com (gatekeeper.solutia.com [205.191.176.56])
	by server.ccl.net (8.8.7/8.8.7) with ESMTP id TAA26066
	for <chemistry@ccl.net>; Fri, 19 May 2000 19:08:23 -0400
From: GGTANA@solutia.com
Received: by soinplx1.solutia.com (NPlex 2.0.123); Fri, 19 May 2000 19:06:11 -0400
X400-Received: by /PRMD=solutia/ADMD=mci/C=us; relayed; 19 May 2000 17:29:13 -0400
To: chemistry@ccl.net
cc: GGTANA@solutia.com
Subject: COSMO in MOPAC7
X400-MTS-Identifier: [/PRMD=solutia/ADMD=mci/C=us;XGW-000519190538-0400-29361]
X400-Content-Identifier: COSMO in MOPAC7
Message-Id: <"0519230538-COSMO in MOPAC7"@MHS>
Date: 19 May 2000 17:29:13 -0400


     Hi,
     
     I have a difficulty to calculate the heat of formation using the COSMO 
     model in MOPAC7 revision 2 (QCPE). Some modification was necessary in 
     some subroutines in order to avoid incompatible common blocks. I use 
     MOPAC on SGI Indigo 2 and could reproduce HOF without COSMO.
     
     There is no problem for small molecules, though test_cosmo yielded a 
     slightly different result than that given in the test_cosmo.dat 
     (60.059 vs 59.766 kcal/mol). 
     
     ESP (= 78.4) alone could not complete the MOPAC job for some 
     molecules. I had to assign smaller NSPA than the default 42. For one 
     molecule for a single point energy calculation (max # of segments = 
     384), it took only 10 sec. for NSPA = 31, while MOPAC runs hours 
     without results for NSPA=32. I couldn't get any for another one (max # 
     of segments = 400) even for NSPA = 4. There was no error message or 
     warning. When I killed the job, the last line of the output file was 
     "RHF CALCULATION,....".
     
     What am I missing? I appreciate your help.
     
     Currently I am not a CCL subscriber, I appreciate if you could reply 
     this e-mail address.
     
     
     Thanks in advance.
     
     
     Genzo