Re: [AMBER] reaidt.f problem

From: Premila Samuel Mohan Dass <ps958.iqb.rutgers.edu>
Date: Sun, 4 Mar 2018 20:54:10 +0000

Here is an example on how to prepare input for Amber’s resp:
$AMBERHOME/bin/resp -O -i 1u7r.in -t 1u7r.chg -p 1u7r.pch -o testresp.out -q 1u7r.qin -e esp_espgen.esp

You get the esp_espgen.esp file using Amber’s espgen on your gaussian log file. To prepare the .in file, you can use Amber’s respgen. The .in file is used to specify the atomic numbers and whether the charges are restrained for specific atoms. In line 9 of the file, the total charge of the system and total number of atoms are specified in the first and 2nd column. In the following lines, the first column has each atomic number corresponding to each of the atoms of residues prepared for Gaussian input. In the second column, if an atom charge is restrained to a specific value to be referred to a .qin file prepared, you give a value of -1. If the charge is allowed to be free-floating, you give a value of 0 in the second column, and if the charge derivation is constrained with another atom to be the same as in the case for methyl hydrogens, then the number of the first atom constrained is given.

You also prepare the qin file specifying the reference charges for each atom for beginning of the fitting. The qin file has a specific format of 8 atoms separated by spaces per line. The format is similar to .chg output you will get.

Some references for resp fitting are below and you can see how they prepared the files required for Amber’s resp:
http://ambermd.org/antechamber/ac.html#respgen
http://ambermd.org/tutorials/advanced/tutorial1/section1.htm

Using RED is probably simpler however.

Premila Samuel


On Mar 4, 2018, at 3:26 PM, Premila Samuel Mohan Dass <ps958.iqb.rutgers.edu<mailto:ps958.iqb.rutgers.edu>> wrote:

[This sender failed our fraud detection checks and may not be who they appear to be. Learn about spoofing at http://aka.ms/LearnAboutSpoofing]

You can use espgen, Amber program, to extract ESP fit center from the gaussian output file. Gaussian does not calculate RESP charges, have to use Amber resp program to do that.
For example, use: espgen -i gaussian.log -o esp_espgen.esp

Then prepare the other input files for use with Amber’s resp program with Amber’s respgen.

Premila P. Samuel, Ph.D.
Post-Doctoral Associate,
Dr. David Case research group,
Center of Proteomics,
Rutgers University,
174 Frelinghuysen Road,
Piscataway, NJ 08854

On Mar 4, 2018, at 3:19 PM, Barış Kurt <bariskurt007.gmail.com<mailto:bariskurt007.gmail.com><mailto:bariskurt007.gmail.com>> wrote:

++and also I calculated RESP charges with gaussian09 now I don't know how
to import them into AMBer as a dat files?

2018-03-04 23:11 GMT+03:00 Barış Kurt <bariskurt007.gmail.com<mailto:bariskurt007.gmail.com><mailto:bariskurt007.gmail.com>>:

Sir, of course first I compiled fortran (otherwise it would be impossible
to get "syntax error") . I don't want to use RED I have no time for
learning another program.

2018-03-04 16:41 GMT+03:00 David A Case <david.case.rutgers.edu<mailto:david.case.rutgers.edu><mailto:david.case.rutgers.edu>>:

On Sun, Mar 04, 2018, Barış Kurt wrote:

I found readit.f on amber web site and copied pasted it but when I try
to
run it I got

syntax error near unexpected token `('

What do you mean by "try to run it"? This is Fortran program, which
would need to be compiled first. Can you say exactly what you did.

But further, this file is 20 years old, and not a part of Amber anymore.
What do you expect it to do? What problem are you trying to solve?

If you want to do RESP calculations, I'd suggest trying the R.E.D. web
site (see pointers at ambermd.org<http://ambermd.org><https://na01.safelinks.protection.outlook.com/?url=http%3A%2F%2Fambermd.org&data=02%7C01%7Cps958%40iqb.rutgers.edu%7C7d1cbf2c654a4c3d269408d5820e4d45%7Cb92d2b234d35447093ff69aca6632ffe%7C1%7C0%7C636557920319893163&sdata=v%2FX7LlQn4GYwMK5P0HHiS8%2FJ5bHKKszYNED59e6E1Vk%3D&reserved=0>).

....dac


_______________________________________________
AMBER mailing list
AMBER.ambermd.org<mailto:AMBER.ambermd.org><mailto:AMBER.ambermd.org>
https://na01.safelinks.protection.outlook.com/?url=http%3A%2F%2Flists.ambermd.org%2Fmailman%2Flistinfo%2Famber&data=02%7C01%7Cps958%40iqb.rutgers.edu%7C3fca8593a53e4752bde008d5820d4fb7%7Cb92d2b234d35447093ff69aca6632ffe%7C1%7C0%7C636557916062348231&sdata=1fnDHf%2F%2Fm8e7QR3mxpeAbqwYITif%2B1LdgfbCpWSCUSM%3D&reserved=0



_______________________________________________
AMBER mailing list
AMBER.ambermd.org<mailto:AMBER.ambermd.org>
https://na01.safelinks.protection.outlook.com/?url=http%3A%2F%2Flists.ambermd.org%2Fmailman%2Flistinfo%2Famber&data=02%7C01%7Cps958%40iqb.rutgers.edu%7C3fca8593a53e4752bde008d5820d4fb7%7Cb92d2b234d35447093ff69aca6632ffe%7C1%7C0%7C636557916062348231&sdata=1fnDHf%2F%2Fm8e7QR3mxpeAbqwYITif%2B1LdgfbCpWSCUSM%3D&reserved=0

_______________________________________________
AMBER mailing list
AMBER.ambermd.org
https://na01.safelinks.protection.outlook.com/?url=http%3A%2F%2Flists.ambermd.org%2Fmailman%2Flistinfo%2Famber&data=02%7C01%7Cps958%40iqb.rutgers.edu%7C7d1cbf2c654a4c3d269408d5820e4d45%7Cb92d2b234d35447093ff69aca6632ffe%7C1%7C0%7C636557920319893163&sdata=yL14mrq3nIfLhZ5cQgQS3jPiLzfU2Dq3hD2VAiySDgE%3D&reserved=0

_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Sun Mar 04 2018 - 13:00:02 PST
Custom Search