VMD-L Mailing List
From: John Stone (johns_at_ks.uiuc.edu)
Date: Fri Sep 25 2015 - 15:02:07 CDT
- Next message: surya narayanan chandrasekaran: "Re: REG: adding hydrogens using TCL"
- Previous message: John Stone: "Re: Segmentation fault error only when executing command within script"
- In reply to: surya narayanan chandrasekaran: "REG: adding hydrogens using TCL"
- Next in thread: surya narayanan chandrasekaran: "Re: REG: adding hydrogens using TCL"
- Reply: surya narayanan chandrasekaran: "Re: REG: adding hydrogens using TCL"
- Messages sorted by: [ date ] [ thread ] [ subject ] [ author ] [ attachment ]
Hi,
To build structures and add missing atoms, you're best options would be to
use something like psfgen or topotools, to run molefacture, or to use
external tools and then pull the modified structure into psfgen/topotools.
While the VMD atom selectin system allows removal of structure easily,
adding is not as easy, and so that's why special purpose tools like
topotools or psfgen exist.
Cheers,
John Stone
vmd_at_ks.uiuc.edu
On Fri, Sep 25, 2015 at 06:11:41PM +0200, surya narayanan chandrasekaran wrote:
> Hi everyone,
> I have a system containing pyrrole like pigment for which
> i am doing QM/MM calculation for predicting excitation energy along the
> trajectory. The problem is the system don't have hydrogen on nitrogen
> according to the amber ff parametrized for it. so on truncating the pigment
> from protein and when given as XYZ coordinates for the QM calculation. The
> QM calculation crashes, because the Multiplicity is 1 apparently odd number
> of electrons. Probably the input coordinate is not recognized as resonance
> structure. So i am planning to add hydrogen to the nitrogen while going to
> the QM calculation using TCL script.
>
> I know how to use the truncation procedure...
>
> set TruncSel [atomselect top {resname XXX and name "CM[A-D]" CAA CBB CAC
> CED}]
>
> $TruncSel set name HHH
> $TruncSel set type HAN
>
> But i don't know how to add Hydrogen using TCL. Any help is appreciated.
>
> Note: the pigment is phycoerythrobillin so don't imagine it's only pyrrole.
-- NIH Center for Macromolecular Modeling and Bioinformatics Beckman Institute for Advanced Science and Technology University of Illinois, 405 N. Mathews Ave, Urbana, IL 61801 http://www.ks.uiuc.edu/~johns/ Phone: 217-244-3349 http://www.ks.uiuc.edu/Research/vmd/
- Next message: surya narayanan chandrasekaran: "Re: REG: adding hydrogens using TCL"
- Previous message: John Stone: "Re: Segmentation fault error only when executing command within script"
- In reply to: surya narayanan chandrasekaran: "REG: adding hydrogens using TCL"
- Next in thread: surya narayanan chandrasekaran: "Re: REG: adding hydrogens using TCL"
- Reply: surya narayanan chandrasekaran: "Re: REG: adding hydrogens using TCL"
- Messages sorted by: [ date ] [ thread ] [ subject ] [ author ] [ attachment ]