> So what type you will choose?
I would either choose some other H type at random (it likely
won't matter), or use Li+. To be procedurally accurate I would
create my own type by using experimental data from neutron
scattering if available, plugging the first radial distribution
peak of H+ and water O into the vdw parameter equation derived
in Ross and Hardin (JACS 116, 6070 (1994)) to get R* for epsilon
of .1, then verify the 1st peak in a simulation with water.
> Aslo, I tried to add a H+ ion by add 3 lines in pdb:
>
> TER
> ATOM 1232 H+ H+ 77 45.303 29.749 7.497 1.00 0.00
> TER
>
> But after I load this pdb and check the charge, it's still 0.
The charge comes from the residue definition. You would need
to define a H+ residue, see amber7/dat/leap/lib/ions94.cmd.
> While a line
> ATOM 1232 Li+ Li+ 77 45.303 29.749 7.497 1.00 0.00
> can change the total charge to +1.
Try setting it to 2 - the charge will remain 1.
Bill
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Tue Feb 17 2004 - 21:53:00 PST
