I modeled a liquid crystal with GB potential under reduced temperature 1.9. It should be random structure according to several published papers. But what I got is perfect ordering structure.
About one week ago, I sent my GB parameters to the mail list and validated these parameters. If the GB parameters are right, the control of pressure or temperature has bugs maybe. Does anyone have the similar opinion with me? Are there some information missed in kinetic energy calculation for temperature or in Virial for pressure? So I tried a series of pressures ( style xyz) and temperatures. Where these bugs are ? thanks a million.