CCL:G: AMSOL Question



Three issues. Two answered as an AMSOL developer, one as a sometime solvation meddler...

1)  AMSOL's Z-matrix format is that used by AMPAC and MOPAC. It's described in lots of places, e.g., http://www.chem.cmu.edu/courses/09-560/docs/msi/modenv/D_Files.html#944883 (under subheading MOPAC).

2)  One does not need to specify atomic partial charges in AMSOL unless one is using the EXTCM (external charge model) keyword to calculate generalized Born solvation free energies with user supplied charges. It is quite possible, however, that example input decks contain such charges. That is because a convenient way to generate a new input geometry is to cut and paste from the .log file of a previously completed calculation. In that case, the partial atomic charges from that prior calculation are indeed found at the end of each atom's line, but that information is not read in the absence of EXTCM.

3)  Jim Kress suggested that Mark Gordon's effective fragment potential (EFP) as an option for a supermolecular solute, where the supermolecule is a cluster of solvent molecules about a central non-solvent-molecule solute. EFP is certainly designed specifically for the purpose of modeling such clusters, but it does NOT include charge transfer to the solute molecules when they are represented by the EFP. Only fully quantum solvent molecules can transfer charge to/from the solute. The EFP may be thought of as a tremendously sophisticated molecular mechanics description (including exchange repulsion, polarizability, etc.) of the explicit solvent molecules, but they do not carry partial charges (unless Mark and his co-workers have an extension that I haven't yet seen). It could certainly be added, in principle, but it introduces some trickiness associated with partial occupation numbers in the quantum MOs.

Chris Cramer


On Jan 28, 2006, at 6:59 PM, Jim Kress ccl_nospam[A]kressworks.com wrote:

Sent to CCL by: "Jim Kress" [ccl_nospam|*|kressworks.com]
Why don't you try the Fragment Method in GAMESS/ PCGAMESS?  It was
specifically designed for system such as yours. 

Jim

-----Original Message-----
From: Mark Zottola mzottola * gmail.com 
Sent: Saturday, January 28, 2006 11:49 AM
To: Kress, Jim 
Subject: CCL:G: AMSOL Question

Sent to CCL by: Mark Zottola [mzottola||gmail.com]
------=_Part_11700_8950891.1138462191404
Content-Type: text/plain; charset=ISO-8859-1
Content-Transfer-Encoding: quoted-printable
Content-Disposition: inline

I am exploring the behaviour of a modestly-sized molecule in 
water.  As the PCM models for water are known to fail to 
account for charge transfer, the other ab initio method would 
be to include a solvation shell.  As this molecule would 
require approximately 10 waters, I would like to get an 
initial answer much quicker than ab initio methods on such a 
system would allow.  Since AMSOL is parameterized to account 
for charge transfer, etc., it seems the obvious choice to 
answer my questions.

To that end, I have a few unresolved questions after reading the AMSOL
manual:

1)  How does one build the z-matrix format that is used in 
AMSOL?  This is not the traditional z-matrix form used in 
Gaussian.  I can build a z-matrix
- the question is how does one build an AMSOL z-matrix?

2)  There appears to be a vicious circle in the creation of 
an AMSOL input deck.  The input deck requires charges from an 
AMSOL run in order to run th= e calculation.  How does one 
get AMSOL charges before running an AMSOL job?  = I am sure I 
am misunderstanding the situation, but not sure where my 
mis-understanding arises.

Thanks for your help!


Mark

------=_Part_11700_8950891.1138462191404
Content-Type: text/html; charset=ISO-8859-1
Content-Transfer-Encoding: quoted-printable
Content-Disposition: inline

<div>I am exploring the behaviour of a modestly-sized 
molecule in water.&nb= sp; As the PCM models for water are 
known to fail to account for charge tra= nsfer, the other ab 
initio method would be to include a solvation shell.&nb= sp; 
As this molecule would require approximately 10 waters, I 
would like to=  get an initial answer much quicker than ab 
initio methods on such a system=  would allow.&nbsp; Since 
AMSOL is parameterized to account for charge tran= sfer, 
etc., it seems the obvious choice&nbsp;to answer my questions.
</div>
<div>&nbsp;</div>
<div>To that end, I have a few unresolved questions after 
reading the AMSOL=  manual:</div> <div>&nbsp;</div> 
<div>1)&nbsp; How does one build the z-matrix format that is 
used in AMSOL?= &nbsp; This is not the traditional z-matrix 
form used in Gaussian.&nbsp; I = can build a z-matrix - the 
question is how does one build an AMSOL z-matrix= ?</div> 
<div>&nbsp;</div> <div>2)&nbsp; There appears to be a vicious 
circle in the creation of an AM= SOL input deck.&nbsp; The 
input deck requires charges from an AMSOL run in = order to 
run the calculation.&nbsp; How does one get AMSOL charges 
before r= unning an AMSOL job?&nbsp; I am sure I am 
misunderstanding the situation, b= ut not sure where my 
mis-understanding arises.&nbsp;=20 </div> <div>&nbsp;</div> 
<div>Thanks for your help!</div> <div>&nbsp;</div> 
<div>&nbsp;</div> <div>Mark</div>

------=_Part_11700_8950891.1138462191404--



-= This is automatically added to each message by the mailing 
script =- To recover the email address of the author of the 
message, please change the strange characters on the top line 
to the ^ sign. You can also look up the X-Original-From: line 
in the mail header.> Conferences: 

Search Messages: http://www.ccl.net/htdig  (login: ccl, 
Password: search)> 
-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
-+-+-+-+-+









-= This is automatically added to each message by the mailing script =-
To recover the email address of the author of the message, please change
the strange characters on the top line to the |-| sign. You can also
look up the X-Original-From: line in the mail header.

E-mail to subscribers: CHEMISTRY|-|ccl.net or use:

E-mail to administrators: CHEMISTRY-REQUEST|-|ccl.net or use

Subscribe/Unsubscribe: 

Before posting, check wait time at: http://www.ccl.net


Search Messages: http://www.ccl.net/htdig  (login: ccl, Password: search)

If your mail bounces from CCL with 5.7.1 error, check:






--


Christopher J. Cramer

University of Minnesota

Department of Chemistry

207 Pleasant St. SE

Minneapolis, MN 55455-0431

--------------------------

Phone:  (612) 624-0859 || FAX:  (612) 626-2006

Mobile: (952) 297-2575

cramer|-|pollux.chem.umn.edu

http://pollux.chem.umn.edu/~cramer

(website includes information about the textbook "Essentials

    of Computational Chemistry:  Theories and Models, 2nd Edition")