AMBER Archive (2008)Subject: AMBER: How to add radius of an atom
From: saurabh agrawal (imsam100_at_yahoo.co.in)
Date: Mon Apr 21 2008 - 00:45:43 CDT
Dear Amber Users,
The molecule I am using for simulation
contains two Selenium atoms (viz. Se1 and Se2). For setting “Se” as an atom type I
have used following commands:
*****************************************
>addAtomTypes { {"Se1"
"Se" "sp2"} }
addAtomTypes: type Se1 - max length is
2
> addAtomTypes { {"Se2"
"Se" "sp2"} }
addAtomTypes: type Se2 - max length is
2
****************************************
Than to set partial charges for the atom which I have these commands:
****************************************
> s1=createAtom Se1 Se 0.37744 > s2=createAtom Se2 Se 0.32396 ****************************************
and to set element following commands **************************************** > set s1 element Se > set s2 element Se **************************************** My questions are-
1. Is this a correct way to add an atom?
2. How to set atomic radius for Se atom? Because when I add solvent in the system It is giving following output:
> solvateoct c TIP3PBOX 8 (using default radius 1.500000 for Se1) (using default radius 1.500000 for Se2) Scaling up box by a factor of 1.159796 to meet diagonal cut criterion Solute vdw bounding box: 44.319 22.848 42.645 Total bounding box for atom centers: 62.876 62.876 62.876 (box expansion for 'iso' is 56.0%) Solvent unit
box: 18.774 18.774 18.774 (using default radius 1.500000 for Se1) (using default radius 1.500000 for Se2) Volume: 130310.325 A^3 (oct) Total mass 66695.892 amu, Density 0.850 g/cc Added 3240 residues.
3.What is the correct value for Se radius?
|
Explore your hobbies and interests. Click here to begin.
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber@scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo@scripps.edu
|