Greetings, I am trying to do electrostatics calculations for a protein system which requires knowledge of the average dipole vector for an experimental probe. I've previously obtained this average dipole vector by using umbrella sampling with WHAM over the chi_2 dihedral angle of this probe, however, it's become clear that the chi_1 orientations are important enough to the average dipole vector that being in a local min for chi_1 and biasing the single degree of freedom is not sufficient. Rather than doing the simulations for 2 DOF, it was suggested that I bias the dipole vector direction relative to the moment of inertia tensor's eigenvectors. Gromacs has a function for orientation restraints relative to the rotation matrix for each frame, however, this is presented in the context of NMR spectroscopy and observables. Would it be possible to use the orientation restraints function for umbrella sampling of dipole vector space and is there a straightforward way to interpret the biasing orientation region (which I believe is referred to as the observable in the manual) as spherical polar coordinates phi and psi?
Thanks, Andrew Ritchie The University of Texas at Austin-- gmx-users mailing list gmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/Search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org. Can't post? Read http://www.gromacs.org/Support/Mailing_Lists