|
|
Joined: Sep 2004
Posts: 57
Forum Member
|
OP
Forum Member
Joined: Sep 2004
Posts: 57 |
Hi,
I am intersted in getting the average structure of two conformers of same molecule giving first conformer more weight than the other. I am reading two coordinate sets and making a trajectory of these two. But using COOR DYNA calculates the average structure of this trajectory (actually two coordinate sets in this trajectory in my case ) assigns equal weighting to both the coordinate sets to give average structure. However, I want to get the average structure giving different weightings to each coordinate sets in the trajectory.
Please guide.
Best Wishes Baloo
|
|
|
|
Joined: Sep 2003
Posts: 8,629 Likes: 24
Forum Member
|
Forum Member
Joined: Sep 2003
Posts: 8,629 Likes: 24 |
|
|
|
|
Joined: Sep 2004
Posts: 57
Forum Member
|
OP
Forum Member
Joined: Sep 2004
Posts: 57 |
Dear Rick,
Thanks for the suggestion. However 'coor average' does not allow to assign different weights to respective structures to calculate the average structure from more than 2 structures.
Best Wishes, Baloo
|
|
|
|
Joined: Sep 2003
Posts: 8,629 Likes: 24
Forum Member
|
Forum Member
Joined: Sep 2003
Posts: 8,629 Likes: 24 |
Your initial post explicitly referred to two conformers. Perhaps you could use COOR DYNA to average each conformer, then COOR AVER to combine the two averages.
|
|
|
|
Joined: Sep 2003
Posts: 4,863 Likes: 10
Forum Member
|
Forum Member
Joined: Sep 2003
Posts: 4,863 Likes: 10 |
How about something like (corman.doc): read coor card name ... coor scale fact1 read coor card comp name ... coor scale fact2 comp coor add comp read coor card name ... coor scale fact3 coor add comp etc, etc, etc ...
You can of course do it in a loop if you have many structures. The weights can be read from a file with the get command (miscom.doc).
Lennart Nilsson Karolinska Institutet Stockholm, Sweden
|
|
|
|
|
|