Dear Dr. Gale,
I am interested in optimizing structures that contain water molecules coordinated to the metal centers of a periodic framework. I started testing the usual models (tip3p, tip4p, and so on), that are rigid. Thus, my system will contain a flexible framework allowed to relax, as well as the position of the waters. However, waters’ bonded terms should be kept as defined by their force field.
I saw the new feature of ‘rigid molecule optimization with quarternions’ in GULP 6.0, but I could not set up correctly the input file. I’ve tried using the ‘molrigid’ definition, but I am not sure about how to properly define (Q1 Q2 Q3) values for the quarternions. Unfortunately, there is no example of an input using the ‘molrigid’ keyword. Besides, although I doubted it were good approaches, I have also tried setting different regions (cartesian/cartesian rigid; cartesian region 1/cartesian region 2 rigid); and defining harmonic potentials for the water bonded terms with really high k constants, unsuccessfully.
Considering my problem (optimization with flexible and rigid molecules):
a) what is the best/correct approach to declare this in the input? If it is the usage of ‘molrigid’ definitions, can you give me an example for properly defining it?
b) If we want to use a water model that contains dummy atoms (massless point charges), what is the correct definition of them? (I usually set an atom with an element different from the remaining ones, and change the values for the mass (to zero), as well as cova and vdw).
Any other comment about properly defining water molecules is really welcome.
Thanks in advance,