From szilagyi@miat0.vein.hu Wed Jan 12 04:18:43 1994
Received: from sztaki.hu  for szilagyi@miat0.vein.hu
	by www.ccl.net (8.6.4/930601.1506) id EAA04840; Wed, 12 Jan 1994 04:16:48 -0500
From: 
Received: by sztaki.hu with SMTP
	(5.67a8/SZTAKI-3.13) id Ba12098; Wed, 12 Jan 1994 10:01:28 +0100
Received: by miat0.vein.hu with SMTP
	(5.61++/Vein-1.0) id AA22620; Wed, 12 Jan 94 10:01:12 +0100
Message-Id: <9401120901.AA22620@miat0.vein.hu>
To: chemistry@ccl.net
Cc: szilagyi@miat0.vein.hu
Subject: METMOD1 force field
Date: Wed, 12 Jan 94 10:01:11 +0100
X-Mts: smtp
X-Charset: US
X-Char-Esc: 0
Status: RO

Dear Netters,

	I am sending now, the first part of the METMOD1 force field.
Anyone can use it who has the PCMODEL V4.0. If someone has only the
MM2 (or similar force field MMX, MM+, etc.), he/she can very easy
install the additional parameters (marked with C). Do not need to 
reoptimize it.

Part I. Modelling of WCl4CR1R2 complexes

PARAM.NEW file:

-------------------------------------------------------------------------------
           MMX (1986)
             (THE PARAMETERS INDICATED WITH AN ASTERISK ARE NOT FINAL.)
         0 ATOMS INDICATE GENERALIZED PARAMETERS
             TORSIONAL PARAMETERS
                        ANGLE          V1          V2          V3

TOR           1      1- 1- 1- 1       0.200       0.270       0.093
        ---------- rest deleted ------------

Some instruction in use

	I marked with a letter C the new (additional) parameters.
Lets see the general formula of the complex:


	Cla		The structure is very close to the
     Cle |   R1		trigonal bipyramidal coordination.
	>W=Cc<		(but there some significant distortion)
     Cle |   R2
        Cla

Atom types:

	W 	tunsgten 	31
	Cle	equatorial Cl	12
	Cla	axial Cl	59
	Cc	carbene C	 2

		anything else is MMX default.

--------------------------------------------
      stuff deleted
--------------------------------------------

The average and RMS fitting error are less than 0.05 A !!
(Large and significant difference is in the W-Cla bond
length. From the ab initio calcns. it seems to be longer
then the equatorial one. In the literature we did not
find any similar effect, so at this time we did not make
any difference between these two bond lenght)

There is a possibility to manage with alpha hydrogen
coordination to metal atom. The user needs redefine
the alpha carbon and hydrogen and the carbene group will
rotate to minimise the distance between W and alpha-H.
(In this particular application it is not really important
but it will be in other systems.)

Test files:

--------------------------------------------
      stuff deleted
--------------------------------------------

Please find the rotational barriers, the precise description of the
parameter optimization in the published papers.

If you want to minimized the structures you need to retype the atoms
and if you want to display it you need to "backtype" the atoms, of
course.

If something is missing or not clear please write me

						Rob