Simulação de dinâmica molecular de vapor de água?

8

Estou tentando fazer MD com vapor de água. Como eu sei, existem alguns modelos de água para água líquida, como SPC, SPC / E, TIP3P, mas eles também se aplicam ao estado de vapor da água? E qual é a diferença de simulações de água líquida e de vapor?

dofine
fonte
Oi dofine. Bem-vindo ao Scicomp! :) Deixe-me ver se entendi sua pergunta corretamente. Você está procurando um modelo de MD para em várias fases simultaneamente ou apenas um que o explique em uma fase gasosa? H2O
Paul
@ Paul Obrigado pelo seu comentário. Acho que estou procurando um modelo que se encaixe na água em várias fases, ou pelo menos na fase líquida e gasosa simultaneamente. Porque estou tentando simular a coexistência de vapor líquido. Embora muitos trabalhos tenham discutido sobre esse tópico, ainda me sinto intrigado com a forma como eles tratam o modelo (por exemplo, CEP / E) na fase gasosa. O modelo original não é adequado apenas para a fase líquida?
Dofine 18/04

Respostas:

6

Certamente, é verdade que a maioria dos modelos de água padrão é "sintonizada" de acordo com as propriedades da fase líquida - eles devem, afinal, ser usados ​​como modelos para água em massa , em vez de modelar propriedades interfaciais.

Dito isto, você pode usar esses modelos de fase líquida para simular o equilíbrio vapor-líquido, desde que esteja disposto a conviver com as inadequações substanciais associadas a eles. Alguns anos atrás, estudei a tensão superficial de modelos de água comuns (link arXiv) e constatei que nenhum deles apresentava tensões superficiais superiores a 62 mN / m (o valor experimental é mais próximo de 72 mN / m). Mesmo o melhor modelo de água atualmente disponível no mercado, o modelo TIP4P / 2005, fornece apenas valores de aproximadamente 65 mN / m (de acordo com as simulações Gromacs, NAMD e LAMMPS).

A principal coisa a lembrar sobre as fases de vapor é que, a baixas temperaturas, existem muito poucas moléculas de água na fase de vapor, uma vez que a densidade é muito baixa. Isso significa que você precisará ter um sistema grande ou um longo tempo de amostragem para obter estatísticas suficientes para cálculos precisos.

aeismail
fonte
2
Métodos avançados de amostragem, como amostragem de guarda-chuva ou metadinâmica, também podem ser úteis para obter dados sobre eventos raros ou transições de energia mais altas.
Aesin
3

Geralmente, você deve escolher o campo de força que reproduza com maior precisão a quantidade de maior interesse para você (constante dielétrica, momento dipolar, tensão superficial, etc.). Este é um excelente guia para todos os diferentes modelos de água e seus respectivos pontos fortes e fracos.

Existem várias dúzias de modelos de água diferentes por aí, e todos sugam de uma maneira especial; então, eventualmente, você só precisa fechar os olhos e escolher um. Minha experiência é que a maioria dos não especialistas acaba usando o padrão defacto, ou seja, o TIP3P / TIP4P, que provavelmente é adequado para a maioria das aplicações.

Quanto à diferença na precisão dos campos de força nas fases gasosa e líquida, uma das fraquezas de todos os campos de força MM é que eles são muito sensíveis ao contexto. Um campo de força que seja preciso para a fase gasosa será sempre menos preciso para a fase líquida e vice-versa. Alguns recursos do campo de força, como a polarização, permitem algum grau de resposta dinâmica às mudanças no contexto molecular, mas mesmo assim nunca haverá um campo de força MM que tenha precisão química em duas fases. Eu acho que a maioria dos simuladores aprende a conviver com um certo nível desses tipos de imprecisões.

Se você deseja uma simulação multifásica muito precisa, precisará realmente começar a estudar os métodos quânticos .

tel
fonte
1
Por outro lado, encontrar um quantummétodo (por exemplo, DFT funcional) que funcione bem em todas as fases da água é um desafio.
Toon Verstraelen
2
@Toon Verstraelen Querido doce Jesus, um desafio impossível e com certeza. Estou iniciando um doutorado em biofísica no outono em grande parte porque quero encontrar maneiras precisas de simular proteínas inteiras , e o fato de que as simulações mais precisas e computacionalmente caras não conseguem nem acertar a água pura, mas me deprimem sem fim .
tel