quinta-feira, 20 de maio de 2010

Densidade da água correta

Estou fazendo simulações da água SPC/E à 298 K e densidade 0.998 g/cm^3, ensemble NVT. Quanto às regras de combinação, estou usando o Lorentz-Berthelot, que é o utilizado pelo Charmm. Parece que é o mais utilizado pelas simulações em água. Engraçado o OPLS ter optado por uma regra de combinação diferente. O coeficiente de difusão parece melhor agora, mas ainda abaixo do esperado. Vamos ver se o resultado fica bom depois que eu termalizar bem, e depois ainda com as simulações NVE.

rho = 3,33894 x 10^22 molec de água/cm^3 = 3,33894 x 10^-2 molec de água/Angs^3

equivale a

29,95 Angs^3/molec de água

Nenhum comentário: