Equação de Langevin

De Física Computacional
Revisão de 11h07min de 17 de outubro de 2022 por Arturuf (discussão | contribs) (Criou página com 'O método numérico escolhido para realizar a integração da equação é conhecido como BAOAB, desenvolvido por Leimkuhler e Mattews <ref> Leimkuhler, B., & Matthews, C. (20...')
(dif) ← Edição anterior | Revisão atual (dif) | Versão posterior → (dif)
Ir para navegação Ir para pesquisar

O método numérico escolhido para realizar a integração da equação é conhecido como BAOAB, desenvolvido por Leimkuhler e Mattews [1] utilizado para resolver equações diferenciais estocásticas.

Ele faz o uso de um método de separação das equações entre as denominadas A, B e O, respectivamente representadas:

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 \vec{p'}\left(t + \frac{\Delta t}{2}\right) = exp(-\gamma \Delta t)\vec{p}\left(t + \frac{\Delta t}{2}\right) + \sqrt{1-exp(-2\gamma \Delta t)}\sqrt{mk_{B}T}\vec{G} }

O 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 \vec{G}} aqui representa um número aleatório Gaussiano que faz o papel da força estocástica.

A equação "A" realiza meio passo no tempo da distância, a "B" realiza um meio passo para o momentum e o "O" contabiliza a contribuição estocástica equação.

Essas equações podem formar vários algoritmos de integração mas o utilizado nesse trabalho será o BAOAB:


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 \vec{p}\left(t + \frac{\Delta t}{2}\right) = \vec{p}(t) + \frac{\Delta t}{2}\vec{f}(t) \quad (1) }
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 \vec{r}\left(t+\frac{\Delta t}{2}\right) = \vec{r}(t) +\frac{\Delta t}{2}\vec{p}\left ( t +\frac{\Delta t}{2}\right )\frac{1}{m} \quad (2) }
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 \vec{p'}\left(t + \frac{\Delta t}{2}\right) = exp(-\gamma \Delta t)\vec{p}\left(t + \frac{\Delta t}{2}\right) + \sqrt{1-exp(-2\gamma \Delta t)}\sqrt{mk_{B}T}\vec{G} \quad (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 \vec{r}(t + \Delta t) = \vec{r}\left(t + \frac{\Delta t}{2}\right) +\frac{\Delta t}{2}\vec{p'}\left ( t +\frac{\Delta t}{2}\right )\frac{1}{m} \quad (4) }
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 \vec{p}(t + \Delta t) = \vec{p'}\left(t + \frac{\Delta t}{2}\right) + \frac{\Delta t}{2}\vec{f}(t) \quad (5) }

É importante lembrar que entre os dois últimos passos é necessário atualizar o termo 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 \vec{f}} , já que ele pode depender de termos já atualizados 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 \vec{p}} ou 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 \vec{r}} .

Referências

  1. Leimkuhler, B., & Matthews, C. (2015). Molecular Dynamics: With Deterministic and Stochastic Numerical Methods. (Interdisciplinary Applied Mathematics; Vol. 39). Springer. https://doi.org/10.1007/978-3-319-16375-8