Análise de Energias DM Microcanônico: mudanças entre as edições

De Física Computacional
Ir para navegação Ir para pesquisar
Linha 7: Linha 7:
<gallery>
<gallery>


Ficheiro:g0.001.jpg|MSD para uma densidade de 5%. É possível perceber o comportamento balístico inicial e então o comportamento linear. O ajuste linear foi realizado no intervalo [10, 50].
Ficheiro:g0.001.jpg|Função de distribuição radial (g(r)) para um dt de 0.001.






Ficheiro:var.jpg| Variância da energia total do sistema versus o delta t associado na simulação.


</gallery>


O potencial de Lennard-Jones tem uma contribuição bastante pequena após determinado ponto. Este ponto onde desconsideramos essa contribuição chamamos Raio de corte(Rc). Em muitos livros o raio de corte do potencial de lennard-jones é usado Rc = 2,5. Mas o ponto que quero levantar neste tópico é que mesmo após Rc = 2,5 ele ainda é relevante para o calculo da energia total, potencial e cinética do sistema. Para demonstrar, usei uma densidade constante e 1000 partículas para montar o gráfico abaixo onde é possível notar a contribuição que esse potencial tem sobre o sistema após o Rc = 2,5.
[[Arquivo:rc.jpg|300px|thumb|right|]]
</gallery>
== Número de Partículas ==
== Número de Partículas ==
== Análise do Raio de Corte ==
== Análise do Raio de Corte ==

Edição das 14h24min de 15 de junho de 2016

Em dinâmica molecular é necessário que se tenha ideia do erro relacionado à limitação computacional que temos em relação à quantidade de partículas que podemos simular em um sistema, ao passo de tempo associado à simulação (dt) e no caso do ensemble microcanônico farei uma análise de como o raio de corte do potencial de Lennard-Jones também pode afetar a energia total do sistema em DM.

Delta t

Para poder mostrar como o comportamento da simulação pode mudar de acordo com o dt associado mostrarei a variância da energia total do sistema (que supostamente deveria 0, pois a energia deveria ser constante). Juntamente com a estrutura do sistema que será mostrado pela g(r).


Número de Partículas

Análise do Raio de Corte

O potencial de Lennard-Jones tem uma contribuição bastante pequena após determinado ponto. Este ponto onde desconsideramos essa contribuição chamamos Raio de corte(Rc). Em muitos livros o raio de corte do potencial de lennard-jones é usado Rc = 2,5. Mas o ponto que quero levantar neste tópico é que mesmo após Rc = 2,5 ele ainda é relevante para o calculo da energia total, potencial e cinética do sistema. Para demonstrar, usei uma densidade constante e 1000 partículas para montar o gráfico abaixo onde é possível notar a contribuição que esse potencial tem sobre o sistema após o Rc = 2,5.

Rc.jpg