Grupo - Lennard Jones: mudanças entre as edições

De Física Computacional
Ir para navegação Ir para pesquisar
Sem resumo de edição
Sem resumo de edição
Linha 42: Linha 42:


Dessa forma, utiliza-se amostragem aleatória massiva para estimar <math> \langle f(x) \rangle </math>, que é a média da função no intervalo de interesse.
Dessa forma, utiliza-se amostragem aleatória massiva para estimar <math> \langle f(x) \rangle </math>, que é a média da função no intervalo de interesse.


=== Amostragem por importância ===
=== Amostragem por importância ===


Se queremos encontrar
Um problema da amostragem simples é que ela utiliza uma distribuição uniforme para
 




Linha 58: Linha 56:
Dado uma amostra com <math>N</math> partículas, a abordagem introduzida por Metropolis segue o seguinte esquema:
Dado uma amostra com <math>N</math> partículas, a abordagem introduzida por Metropolis segue o seguinte esquema:


(1) Selecionar uma partícula aleatóriamente, e calcular sua energia <math>U(r)</math>;
(1) Selecionar uma partícula aleatóriamente, e calcular sua energia <math>U(r)</math>;


(2) Dado o deslocamento <math>r_n = r + \Delta</math>, calcular <math>U(r_n)</math>;
(2) Dado o deslocamento <math>r_n = r + \Delta</math>, calcular <math>U(r_n)</math>;


(2) Aceitar o movimento <math>r \rightarrow r_n</math> com probabilidade <math>p = min\{1; \exp[-\beta (U(r_n) - U(r))]\}</math>
(3) Aceitar o movimento <math>r \rightarrow r_n</math> com probabilidade <math>p = min\{1; \exp[-\beta (U(r_n) - U(r))]\}</math>


=== Estimadores no Equilíbrio ===
=== Estimadores no Equilíbrio ===

Edição das 13h23min de 13 de janeiro de 2018

O potencial devido a interação entre duas partículas separadas por uma distância Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle r'} pode ser modelado pelo potencial de Lennard Jones:

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle U(r') = 4 \epsilon \left [ \left ( \frac{\sigma}{r'} \right )^{12} + \left ( \frac{\sigma}{r'} \right )^{6} \right ].}

Posto em unidades reduzidas (Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle r \equiv r'/ \sigma} e Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle U \equiv U' / \epsilon} ), o potencial reduz-se a:

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle U(r) = 4 \left [ \left ( \frac{1}{r} \right )^{12} + \left ( \frac{1}{r} \right )^{6} \right ].}

Trabalha-se, por conveniência, com o seguintes sistema de unidades básicas:

Grandeza Comprimento Tempo Massa Temperatura Energia Pressão Densidade
Unidade Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \sigma} Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \sigma \sqrt{m_p / \epsilon}} Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle m_p} Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \epsilon/k_B} Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \epsilon} Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \epsilon / \sigma^3} Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle 1 / \sigma^{3}}

onde Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle m_p} é a massa da partícula e Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle k_B} é a constante de Boltzmann. .

Método Monte Carlo

Denomina-se método de Monte Carlo métodos estatísticos que se baseiam em amostragem aleatória massiva para cálculo numérico.

Amostragem simples

Pode-se querer calcular uma integral numericamente utilizando Monte Carlo. Uma forma de fazer isso parte de que uma integral pode ser reescrita como

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle F = \int_a^b{f(x) dx} = (b - a)\langle f(x) \rangle }

Dessa forma, utiliza-se amostragem aleatória massiva para estimar Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \langle f(x) \rangle } , que é a média da função no intervalo de interesse.

Amostragem por importância

Um problema da amostragem simples é que ela utiliza uma distribuição uniforme para


A integração Monte Carlo por amostragem aleatória com distribuição de probabilidade uniforme consiste em fazer o uma média da função a ser integrada e multiplicar pela largura 

pode-se ter problemas por realizar muitas vezes

Algoritmo de Metropolis

Dado uma amostra com Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle N} partículas, a abordagem introduzida por Metropolis segue o seguinte esquema: 

(1) Selecionar uma partícula aleatóriamente, e calcular sua energia Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle U(r)}
;

(2) Dado o deslocamento Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle r_n = r + \Delta}
, calcular Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle U(r_n)}
;

(3) Aceitar o movimento Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle r \rightarrow r_n}
 com probabilidade Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle p = min\{1; \exp[-\beta (U(r_n) - U(r))]\}}

Estimadores no Equilíbrio

Detalhes Técnicos

Condições de Contorno

Truncagem nas interações

Translação

Diagramas de fase

Referências

Disponível em “http://fiscomp.if.ufrgs.br/index.php?title=Grupo_-_Lennard_Jones&oldid=1233”