AMBER Archive (2008)
Subject: AMBER: RESP GAUSSIAN AMBER
I am just beginning to study QM, so I am not an expert. Please be gentle!
I follow this procedure:
1) "antechamber -fi pdb -i ligand.pdb -fo gcrt -o ligand.com -pf no" to generate my input template
2) I edit the gaussian options to: #HF/6-31g* scf=tight test pop=mk iop(6/33=2) iop(6/42=6) opt
To my understanding: 6/42=6 specifies the density points in each layer, 6/33=2 (which is NOT documented in gaussian 03) prints out potential points and potentials. On http://ambermd.org/tutorials/advanced/tutorial/section1.htm is is also recommended to use iop(6/41=10) on top, which would add 10 concentric layers of points for each atom. Is this option recommended or is it a carryover?
3) I run Gaussian03 which outputs a ligand.log file
4) I run 'antechamber -fi gout -fo prepi -i ligand.log -o ligand.prepi -c resp -s 2 -pf no'
and voila the resp charges are found in the prepi file. These charges are similar to AM1-BCC (which I suppose they should be as AM1-BCC has been parametrised to reproduce RESP).
To get a feeling of whether this procedure is correct I had a look at /usr/local/amber9/examples/resp_charge_fit/water. I created a water molecule in molden and applied the same procedure (this time with and without option iop(6/41=10)). The obtained charge are basically identical:
According to /usr/local/amber9/examples/resp_charge_fit/water/* I should get:
Looking at the gaussian optimised geometry of the water molecule I find small deviatons in the angle H1-O1-H2 compared to the standart water angle. What is happening here? Can I trust the RESP charges for my own ligands? Where can I find
I apologize if the mailing list is not the right place to post. However, I am grateful for any pointers!
-- Sensationsangebot verlängert: GMX FreeDSL - Telefonanschluss + DSL für nur 16,37 Euro/mtl.!* http://dsl.gmx.de/?ac=OM.AD.PD003K1308T4569a ----------------------------------------------------------------------- The AMBER Mail Reflector To post, send mail to amber_at_scripps.edu To unsubscribe, send "unsubscribe amber" (in the *body* of the email) to majordomo_at_scripps.edu