Acessibilidade / Reportar erro

Quatro abordagens para o movimento browniano

Four approaches to the brownian motion

Resumos

A "dança" aleatória de pequenas partículas em suspensão num líquido, um fenômeno conhecido como movimento browniano, foi primeiramente explicado por Einstein na sua famosa tese de doutorado. Seguindo uma perspectiva histórica, mostramos como este fenômeno pode ser descrito de quatro maneiras distintas, a saber: o tratamento difusivo de Einstein, a variante estocástica ou de força flutuante proposta por Paul Langevin, a abordagem via equação de Fokker-Planck, e finalmente, as caminhadas aleatórias de Mark Kac. Discutiremos também as limitações presentes na abordagem difusiva. Em particular, mostramos que a equação parabólica na qual Einstein baseou sua explicação deve ser substituída por uma equação do tipo hiperbólica que também surge naturalmente no tratamento via caminhadas aleatórias. A solução geral dessa equação é obtida e comparada com o resultado padrão. Para tempos curtos, comparados com as escalas de tempo características do sistema, o movimento das partículas segue um comportamento ondulatório.

Einstein; movimento browniano; Paul Langevin; Fokker-Planck


The random trajectory of small particles in suspension within a liquid, a phenomenon known as brownian motion, was firstly explained by Einstein in his famous doctorate thesis. From a historical perspective, we show how such a phenomenon can be described in four different ways, namely: the Einstein diffusive treatment, the stochastic variant or fluctuating force proposed by Paul Langevin, the approach through the Fokker-Planck equation, and, finally, the random walks by Mark Kac. Some limitations present in the standard diffusive approach are also discussed. In particular, we show that the parabolic equation in which Einstein based his explanation should be replaced by a hyperbolic equation of motion which also appears naturally in the treatment of random walks. The general solution of the generalized diffusion equation is obtained and compared to the standard result. For short times, in comparison with the characteristic time scales of the system, the motion of the particles is described by a wave behavior.

Einstein; brownian motion; Paul Langevin; Fokker-Planck


ARTIGOS GERAIS

Quatro abordagens para o movimento browniano

Four approaches to the brownian motion

J.M. SilvaI,1 1 E-mail: jmsilva@dfte.ufrn.br. ; J.A.S. LimaII

IDepartamento de Física, Universidade Federal do Rio Grande no Norte, Natal, RN, Brasil

IIDepartamento de Astronomia, Universidade de São Paulo, São Paulo, SP, Brazil

RESUMO

A "dança" aleatória de pequenas partículas em suspensão num líquido, um fenômeno conhecido como movimento browniano, foi primeiramente explicado por Einstein na sua famosa tese de doutorado. Seguindo uma perspectiva histórica, mostramos como este fenômeno pode ser descrito de quatro maneiras distintas, a saber: o tratamento difusivo de Einstein, a variante estocástica ou de força flutuante proposta por Paul Langevin, a abordagem via equação de Fokker-Planck, e finalmente, as caminhadas aleatórias de Mark Kac. Discutiremos também as limitações presentes na abordagem difusiva. Em particular, mostramos que a equação parabólica na qual Einstein baseou sua explicação deve ser substituída por uma equação do tipo hiperbólica que também surge naturalmente no tratamento via caminhadas aleatórias. A solução geral dessa equação é obtida e comparada com o resultado padrão. Para tempos curtos, comparados com as escalas de tempo características do sistema, o movimento das partículas segue um comportamento ondulatório.

Palavras-chave: Einstein, movimento browniano, Paul Langevin, Fokker-Planck.

ABSTRACT

The random trajectory of small particles in suspension within a liquid, a phenomenon known as brownian motion, was firstly explained by Einstein in his famous doctorate thesis. From a historical perspective, we show how such a phenomenon can be described in four different ways, namely: the Einstein diffusive treatment, the stochastic variant or fluctuating force proposed by Paul Langevin, the approach through the Fokker-Planck equation, and, finally, the random walks by Mark Kac. Some limitations present in the standard diffusive approach are also discussed. In particular, we show that the parabolic equation in which Einstein based his explanation should be replaced by a hyperbolic equation of motion which also appears naturally in the treatment of random walks. The general solution of the generalized diffusion equation is obtained and compared to the standard result. For short times, in comparison with the characteristic time scales of the system, the motion of the particles is described by a wave behavior.

Keywords: Einstein, brownian motion, Paul Langevin, Fokker-Planck.

1. Introdução

O movimento irregular de pequenas partículas imersas numa solução foi originalmente observado em 1828 pelo botânico inglês Robert Brown [1]. Ele notou que as partículas em suspensão adquiriam uma espécie de movimento errático que posteriormente ficaria popularmente conhecido pelo nome de movimento browniano (MB).

Nas décadas seguintes, inúmeras tentativas foram realizadas para desvendar a natureza do movimento browniano. Experimentos de laboratório mostraram que o movimento fica mais intenso quando se reduz a viscosidade do meio ou o tamanho das partículas brownianas, e também quando se eleva a temperatura da solução. Muitas causas possíveis foram aos poucos sendo eliminadas, tais como: atrações ou repulsões entre as partículas suspensas, ações capilares ou higrométricas, bolhas temporárias de ar, correntes de conveção no interior da solução, gradientes de temperatura ou algum tipo de perturbação mecânica, além de outros tipos de instabilidades no fluido.

Somente a partir de 1860 começou a tomar forma o ponto de vista moderno de que o "zigue-zague" das partículas brownianas poderia ser devido às colisões com as moléculas do fluido. Verificou-se que suas trajetórias não apresentavam tangentes (ou seja, as curvas não seriam diferenciáveis), e também que o movimento randômico aparentemente nunca cessava. No entanto, a verdadeira causa do fenômeno permaneceu um mistério até 1905, quando finalmente foi elucidado por Einstein no seu artigo de 1905.

O tratamento de Einstein para o MB é um dos trabalhos intelectuais mais notáveis de todos os tempos. Sua solução representou um grande avanço científico nos campos da química e da física, tornando a teoria atômico-molecular uma parte fundamental da estrutura da matéria. Como uma espécie de bônus extra, o tratamento de Einstein também forneceu uma estimativa do número de Avogadro que foi verificada, com grande precisão, nos experimentos efetuados por Jean Perrin [3].

Posteriormente, um esforço considerável foi canalizado por muitos investigadores para generalizar e compreender o tratamento de Einstein do movimento browniano. Importantes contribuições foram dadas por Smoluchowski [4], Langevin [5], Fokker [6], Burger [7], Fürther [8], Ornstein [9], Planck [10], Kac [11] e muitos outros.

Atualmente, o movimento browniano permanece na fronteira da pesquisa como um exemplo importante de processo estocástico, e constitui uma ferramenta fundamental para o estudo de sistemas físicos de não equilíbrio. Tais sistemas são encontrados em diferentes áreas da física, desde o nível microscópico, como verificado na difusão de partículas num solvente, até escalas de ordem astronômica, tal como observado em sistemas estelares [12]. Um exemplo interessante desse último tipo é representado por um buraco negro (BN) no centro de um sistema estelar denso. Teoricamente, quando sua massa é muito grande, o BN pode adquirir um movimento que é semelhante ao de uma partícula em suspensão num líquido ou num gás [13]. Em cosmologia, movimentos brownianos com barreiras fixas ou móveis são também bastante empregados para estudar os processos de formação da estrutura de larga escala, tais como galáxias, aglomerados de galáxias e vazios [14]. Mais recentemente, outros tipos de contribuições foram obtidas na investigação de sistemas com memória, objetivando estabelecer relações entre os regimes de difusão anômala e normal [15].

Neste contexto, no presente trabalho mostraremos como é possível abordar o movimento browniano de quatro maneiras distintas, a saber: o tratamento difusivo de Einstein, o procedimento estocástico ou de força flutuante proposto por Paul Langevin, a abordagem via equação de Fokker-Planck, e finalmente, as caminhadas aleatórias de Mark Kac [11]. Discutiremos também com bastante detalhe, as limitações presentes na abordagem difusiva. Em particular, mostraremos que a equação parabólica na qual Einstein baseou sua explicação deve ser substituída por uma equação do tipo hiperbólica que também surge naturalmente no tratamento via caminhadas aleatórias.

2. MB e equação de difusão: O tratamento de Einstein

Para estudar o comportamento irregular das partículas em suspensão que surge devido aos movimentos moleculares térmicos, suporemos que cada partícula executa um movimento completamente independente das outras partículas. Como veremos, essa hipótese é válida somente se os intervalos de tempos considerados não são demasiadamente pequenos. Seguindo Einstein [2], consideraremos um intervalo de tempo t, que é pequeno em comparação com o tempo de observação, porém suficientemente longo, para que os movimentos executados por diferentes partículas no neste intervalo de tempo possam ser considerados eventos independentes.

Suponhamos que existam N partículas em suspensão no líquido. No intervalo de tempo t, as coordenadas x das partículas variam de Dx = m, onde m pode assumir valores diferentes (positivo ou negativo) para cada partícula. Uma determinada lei de distribuição de probabilidades deve ser satisfeita pela variável m: A fração de partículas que sofre um deslocamento entre x e x+m no intervalo de tempo t, pode ser expressa por uma equação da forma [2]

com a distribuição f(m) satisfazendo a condição de normalização

onde f(m) é uma função par, f(m) = f(-m), suposta diferente de zero apenas para pequenos valores de m. Considere também que h(x,t) é o número de partículas por unidade de comprimento, e calculemos a distribuição de partículas no instante t+t, a partir da distribuição delas no instante t. Pela definição da função f(m), o número de partículas que no instante t+t se encontram entre x e x+m, é dado por

Como t é muito pequeno, podemos fazer uma expansão temporal de h até segunda ordem2 2 Einstein obteve seus resultados fazendo a expansão no tempo somente até primeira ordem [2]. Por razões que serão discutidas adiante, consideraremos termos até segunda ordem em t na expansão da funcão h( x,t).

e como m também é pequeno, para sermos consistentes devemos desenvolver h(x+m,t) em potências até segunda ordem em m

Inserindo os resultados acima na Eq. (3) obtemos

No lado direito dessa equação, o segundo termo é identicamente nulo uma vez que f(m) = f(-m). Logo, considerando a Eq. (2), vemos que h satisfaz a seguinte equação diferencial

onde definimos

A Eq. (7) representa uma espécie de difusão generalizada. A quantidade h(x,t) é a concentração de partículas por unidade de comprimento em torno de x num instante arbitrário e a constante D é o coeficiente de difusão. No limite

a Eq. (7) se reduz para

que é a forma padrão da equação de difusão, na qual Einstein baseou a sua explicação do MB. A Eq. (7) é do tipo hiperbólica e generaliza a equação de difusão usual que é do tipo parabólica.3 3 Uma classificacão das equacões diferenciais parciais pode ser vista na Ref. [16]. Na seção final desse trabalho analisaremos a solução analítica da equação de difusão generalizada. Por enquanto, prosseguiremos com a descrição einsteiniana do movimento browniano.

Como um exemplo para ilustrar esse tratamento, vamos obter a solução da Eq. (10) quando o processo difusivo satisfaz a seguinte condição inicial

onde N é o número total de partículas e d denota a função delta de Dirac. Como seria esperado, tal condição implica que

A solução da Eq. (10) pode ser facilmente obtida pela técnica das integrais de Fourier. De acordo com esse método, a concentração pode ser definida como

onde os coeficientes da expansão, hk(t), são determinados pela transformada inversa

Calculando as derivadas temporal e espacial de h(x,t) e substituindo suas expressões na equação de difusão (10), obtemos a seguinte forma integral

Como a equação acima é válida para todo instante, seu integrando deve ser identicamente nulo, ou seja,

cuja solução é da forma

Com este resultado, a definição (13) pode ser escrita como

Por outro lado, considerando que

temos para a transformada inversa

e da Eq. (18) podemos escrever

Note que a segunda parcela na expressão acima é igual a zero, pois se trata do produto de uma função par por uma função ímpar, com a equação se reduzindo para

A integração deste resultado é mais facilmente obtida considerando as seguintes mudanças de variáveis: k = y, m = x-x' e a = Dt, com a expressão (22) tomando a seguinte forma

Finalmente, observando que a condição (11), implica que as partículas estão inicialmente localizadas na origem, ou seja, h(x',0) = Nd(x'), a concentração pode ser escrita como

O resultado acima nos mostra que as partículas se comportam como num processo gaussiano difusivo. A função h(x,t) inicialmente representa uma delta centrada em torno da origem x = 0. No entanto, à medida que o tempo passa a distribuição evolui como uma gaussiana de largura variável (ver Fig. 2).


Tendo calculado a função h(x,t), é interessante determinar a distribuição de probabilidade de que uma partícula da amostra ocupe a posição entre x e x+m, quando em t = 0, iniciou seu movimento da posição x0 com velocidade inicial v0. O conhecimento de tal função é de fundamental importância para se calcular quantidades de interesse físico, tais como o deslocamento quadrático médio e a variância. A distribuição de probabilidade pode ser obtida dividindo-se a concentração pelo número total de partículas. Ou seja,

Comparando essa equação com a distribuição de probabilidades gaussiana

vemos que < x > = 0, enquanto que para a variância temos s2 = 2Dt.

Este resultado significa que na teoria do movimento browniano, as grandezas fisicamente relevantes estão diretamente relacionadas com os primeiros e os segundos momentos da distribuição, que é uma propriedade geral da função gaussiana [16]. Tais momentos podem ser calculados da relação

Utilizando a função distribuição (26), o valor de < x > e s2 podem ser obtidos diretamente por cálculos algébricos considerando a expressão geral acima. O primeiro momento é o deslocamento médio4 4 O integrando de (28) é composto pelo produto de uma funcão par por uma funcão ímpar.

Seguindo a mesma prescrição, o segundo momento da distribuição é o deslocamento quadrático médio5 5 Note que a dx = para n par.

que na teoria do MB é conhecida como relação de Einstein. O coeficiente de difusão D na Eq. (29) deve ser uma função da temperatura e da geometria das partículas. Einstein mostrou que para partículas esféricas de raio a, o coeficiente D pode ser calculado a partir da mobilidade b e da temperatura do meio no qual a partícula se encontra. O parâmetro b pode ser obtido da fluidodinâmica, mais precisamente a partir da lei de Stokes [19]. A relação satisfeita por D é

onde kB é a constante de Boltzmann, T é a temperatura, g representa o coeficiente de viscosidade do meio e b = 1/6pga. Inserindo a Eq. (30) na Eq. (29), temos para o deslocamento quadrático médio no MB

Note que para obtermos a forma originalmente deduzida por Einstein [2], utilizamos o fato de que a constante de Boltzmann kB pode ser escrita como kB = R/Na, onde R é a constante universal dos gases e Na é o número de Avogadro.

Portanto, vemos que a partícula se comporta como um processo difusivo com < x2 > µ t. Toda essa formulação unidimensional pode ser consistentemente generalizada para três dimensões. Neste caso, não é difícil demonstrar que a Eq. (31) assume a seguinte forma [19]

É importante também mencionar que o resultado de Einstein (31), ou equivalentemente, (32), foi um dos primeiros exemplos de uma relação onde uma flutuação quadrática média está associada com um processo dissipativo (descrito pelo coeficiente de viscosidade g). Além disso, como os valores das quantidades < r2 > , t, g e a são diretamente mensuráveis, isto significa que o número de Avogadro pode ser estimado se tivermos um bom cronômetro e um microscópio [3]. Seguindo esse procedimento, Jean Perrin [20] obteve valores experimentais do desvio quadrático médio que permitiram uma determinação mais precisa do número de Avogadro [20,21]. Tais resultados também contribuíram significativamente para que a hipótese atômico-molecular tivesse aceitação geral como uma descrição realista da matéria.

Posteriormente, Einstein observou que seus resultados apresentavam inconsistências para tempos curtos comparados aos tempos característicos do sistema. Uma forma simples de perceber tais dificuldades é calculando a "velocidade média" da partícula usando a Eq. (31)

Vemos que no limite t ® 0, a velocidade v ® ¥, sendo esta a raiz da dificuldade. Outra maneira fácil de compreender este problema pode ser vista na Fig. 2. Note que, para tempos próximos de zero temos uma função d de Dirac centrada em x = 0, mostrando que inicialmente todas as partículas estão localizadas na origem. Por outro lado, para intervalos de tempos tão pequenos quanto se queira (t = 0+), a curva é uma gaussiana que se estende a todo espaço, indicando que as partículas se difundiram com uma velocidade infinita. Portanto, fica claro que os resultados de Einstein só permanecem válidos para um regime de tempo suficientemente longo em comparação à escala de tempo característica do sistema.

Para corrigir tais dificuldades, precisaríamos considerar o termo de derivada segunda com respeito ao tempo na equação de difusão (10). Em outras palavras, é importante considerar a solução analítica da Eq. (7), já que ela incorpora naturalmente o termo ¶2h/¶t2, sugerindo que para tempos curtos teremos uma descrição ondulatória. Discutiremos alguns detalhes dessa abordagem na seção final. Por enquanto, vamos proseguir examinando as diversas variantes da teoria do movimento browniano.

3. MB e forças flutuantes: a visão de Langevin

Poucos anos após o trabalho de Einstein, o físico francês Paul Langevin [5], posteriormente seguido por Fürther [8], Ornstein [9] e outros mais, iniciaram uma série de estudos tentando uma possível generalização daqueles resultados. Tal abordagem, comumente conhecida como tratamento de Langevin, será examinada com detalhe nesta seção.

Segundo Langevin, o MB de uma partícula na ausência de um campo de força conservativo pode ser entendido com base numa equação diferencial estocástica, agora popularmente conhecida como equação de Langevin [5,12,18,22]

onde v denota a velocidade da partícula. Nesta equação, a influência do meio sobre o movimento da partícula é decomposta em duas partes. Em primeiro lugar, existe uma força que varia lentamente, F = -gv, representando uma fricção dinâmica sobre o movimento da partícula, onde g é o coeficiente de viscosidade do meio. Existe também uma força aleatória, x(t), que varia rapidamente em comparação com os tempos de observação. Em outras palavras, x(t) é uma força flutuante que é uma característica básica de uma equação diferencial estocástica. Langevin definiu as propriedades dessa função por duas condições (G é uma constante)

que caracterizam o chamado ruído branco.6 6 O ruído é branco ("white noise") se o espectro de potência S(w) da funcão correlacão < x( t)x(0) > é independente da frequência, sendo S(w) = < x( t)x(0) > dt. No tratamento de Langevin, < x( t)x(0) > = Gd( t), temos S(w) = G.

Para determinar a solução analítica da Eq. (34), vamos primeiramente supor uma equação de Langevin geral escrita na seguinte forma

onde f(t) é uma função arbitrária. Definindo

sendo g(t) também arbitrária e (t) sua derivada, a Eq. (37) pode ser reescrita como

Note que a Eq. (34) é recuperada para g(t) = e g t. Portanto, a equação acima, ou equivalentemente, a Eq. (37), pode ser reescrita na forma

que pela mudança de variável, u = ve g t, se reduz a forma elementar

com solução

Retornando para a antiga variável v, vemos que a solução geral da equação de Langevin é dada por

O valor médio e a variância na velocidade devem ser calculados através das propriedades da função x(t). Utilizando a condição (35) temos

A variância é mais facilmente obtida calculando-se primeiramente a diferença v- < v > , de onde obtemos

ou ainda

Tomando a média, utilizando a condição (36) e efetuando a integração, obtemos facilmente

onde (Dv)2 = < v2 > - < v > 2 é a variância. Para calcular a constante G, observemos que o regime estacionário é obtido para tempos longos em comparação com os tempos de flutuação da função x(t), indicando que < v > se anula na Eq. (44), e da Eq. (47) temos que

Por outro lado, o teorema da equipartição garante que a energia cinética média de uma partícula em movimento corresponde a

kBT para cada grau de liberdade, mais precisamente

onde kB é a constante de Boltzmann. Combinando as duas últimas expressões, obtemos a relação exata entre G e a temperatura do meio externo

Uma vez determinada a variância da velocidade é conveniente calcular o deslocamento quadrático médio, já que este é uma grandeza, experimentalmente mensurável (mais detalhes nessa abordagem pode ser vista em na Ref. [22])

onde x0 é a posição da partícula em t = 0. Substituindo na integral acima o valor de v(t) dado pela expressão (43), segue que

Desta equação obtemos o deslocamento médio

sendo o deslocamento quadrático médio obtido calculando-se primeiramente a diferença

de onde obtemos

Na expressão acima, tomando a média, usando a condição (36) e efetuando as integrais obtemos facilmente

Observe que no limite de tempos longos o termo dominante é o primeiro, mais precisamente

ou equivalentemente,

que é a relação de Einstein (ver Eq. (29)). Vemos portanto, que no regime de tempos longos a abordagem de Langevin é equivalente à descrição de Einstein. Neste limite também pode ser mostrado [22] que a distribuição de probabilidades relativa à variável v obedece a uma distribuição maxweliana de velocidades

4. A equação de Fokker-Planck

Como vimos na seção 3, a equação de Langevin na forma (34) descreve o movimento de uma partícula de massa m imersa num fluido com coeficiente de viscosidade g. Este mesmo sistema pode ser descrito por uma equação de movimento que governa a evolução temporal de uma distribuição de probabilidade. Tal equação é comumente conhecida como equação de Fokker-Planck e constitui o objeto de investigação desta seção. A equação de Fokker-Planck é um tipo especial de equação mestra [17,22], freqüentemente usada como uma boa aproximação para descrever processos markovianos mais gerais.

Considere uma equação do tipo Langevin da seguinte forma

onde a variável x denota uma coordenada generalizada que, em princípio, pode ser a posição ou velocidade. Para esta variável independente, a equação de Fokker-Planck dependente do tempo é comumente escrita como [17,22]

onde f(x) relaciona a natureza da força atuando na Eq. (60) e h(x,t) representa a distribuição de probabilidade de encontrar a partícula no intervalo entre x e x+m. A equação acima também pode ser reescrita como

que representa uma equação de continuidade para a densidade de probabilidade h(x,t), na qual a quantidade S(x,t) deve ser interpretada como uma corrente de probabilidade definida por

A integração da Eq. (62) com x assumindo valores no intervalo [a,b] nos fornece

e como

segue que

nos mostrando que a conservação da probabilidade total é, uma conseqüência direta das condições de contorno.

Vamos determinar a solução da equação de Fokker-Planck na forma (62) para o caso estacionário, considerando que os valores extremos S(x = a,t) e S(x = b,t) são nulos. Nestas condições, segue da Eq. (63) que

cuja solução é

onde a constante A é fixada pela condição de normalização de h(x). Para o caso de uma força viscosa, f = -gv e a constante G dada pela Eq. (50), a solução acima assume a seguinte forma

que é a distribuição maxwelliana de velocidades.

A solução não estacionária é obtida diretamente da Eq. (61). Utilizando a mesma força viscosa do exemplo acima, tal equação pode ser representada como

com solução dada por [17,22,25]

Comparando a expressão acima com a distribuição gaussiana (veja Eq. (26)), vemos que os valores da média < v > e da variância (Dv)2 são respectivamente

e

que são os mesmos valores obtidos no tratamento de Langevin (cf. Eqs. (44) e (47). Como seria esperado, vemos também de (71) que para tempos suficientemente longos o sistema relaxa para o estado de equilíbrio, pois a distribuição de probabilidades se reduz à distribuição de velocidades maxwelliana.

5. Caminhadas aleatórias: o tratamento de M. Kac

O problema do caminhante aleatório, é dotado de um caráter bastante universal em física. No magnetismo, por exemplo, um átomo de spin 1/2 tem um momento magnético e de acordo com a mecânica quântica, o spin pode está "up" ou "down", com respeito a uma dada direção. Se essas possibilidades são igualmente prováveis, então qual o momento magnético médio < > para uma amostra contendo N átomos? Um outro problema bastante familiar, corresponde à difusão de partículas num meio intermolecular. Suponha que uma partícula percorre uma distância média l entre duas colisões sucessivas com as moléculas do meio. Qual será a distância percorrida após N colisões?

A solução para o problema da caminhada aleatória, na sua forma mais geral, é facilmente entendido considerando-se a versão mais simples do problema em uma dimensão, tal como originalmente investigado por M. Kac [11]. Suponha que um caminhante aleatório partindo da origem e se deslocando em linha reta, realiza n1 passos de comprimento fixo l para a direita com probabilidade p e n2 passos para a esquerda com probabilidade q = 1-p, de modo que p+q = 1. Além do mais, estamos considerando que os passos são eventos mutuamente independentes. O problema é determinar qual a probabilidade PN(m) de encontrar o caminhante na posição x = ml, onde -N < m < N, depois de ter dado N passos. O número total de passos é

sendo m a grandeza que parametriza a distância líquida percorrida, isto é,

e como cada passo tem comprimento l, a distância que o caminhante percorre a partir da origem é dada por

Considerando que os passos são estatisticamente independentes, de probabilidades p e q, a probabilidade de realizar n1 passos para a direita e n2 passos para a esquerda é independente da seqüência de passos e pode ser escrita como [26]

Existem várias maneiras de arranjar os N passos de forma que n1 seja o número de passos para a direita e n2 seja o número de passos para a esquerda. Na verdade, descobrir o número de maneiras de arranjar os n1 e n2 passos, é descobrir de quantas maneiras distintas podem ser arranjados n1 +n2 objetos, sendo que a permutação de qualquer um dos objetos (n1 +n2) é irrelevante. Tal fato significa que o número de possibilidades distintas é exatamente [16]

e que a probabilidade total, PN(n1), de realizar n1 passos para a direita e n2 para a esquerda num total de N passos, em qualquer ordem, é dada pelo produto

pois todas as seqüências são independentes. Como vemos, o valor de PN(n1) é uma distribuição binomial. Lembrando que a expansão binomial de (p+q)N, onde p e q são dois números quaisquer, é dada por

segue que a distribuição PN(n1) é normalizada, ou seja,

Vamos determinar a probabilidade PN(m) do caminhante se encontrar na posição x = ml. Das Eqs. (74) e (75), temos

Substituindo esses resultados na Eq. (79), pode ser visto facilmente que a distribuição PN(m) tem a forma

ou, equivalentemente,

Para estabelecer uma conexão com o fenômeno de difusão, é necessário descrever o problema do caminhante aleatório por meio de uma equação diferencial envolvendo variáveis contínuas [11,22,26]. Suponha que t seja o tempo necessário para realizar um passo, então PN(m) dado pela Eq. (84) é a probabilidade da partícula se encontrar na posição x = ml no tempo Nt. Somente uma partícula que esteja em x = (m-1)l ou x = (m+1)l no tempo t = (N-1)t poderá atingir a posição x = ml. No passo seguinte, a probabilidade PN(m) obedece a seguinte relação de recorrência [11]

que representa um exemplo típico de um Processo Markoviano7 7 Nos chamados processos markovianos não existem efeitos de memória, ou seja, a probabilidade condicional relativa a cada variável aleatória qi( t) de uma partícula, só depende do valor de qi = num instante anterior t 0 [17,18]. . Equações estocásticas dessa natureza, nas quais os detalhes da dinâmica de um sistema físico são substituídos por leis probabilísticas, desempenham um papel extremamente importante no estudo de sistemas fora do equilíbrio.

Conforme visto anteriormente, se N é suficientemente grande, a função discreta PN(m) pode ser substituída por uma função contínua h(Nt,ml) = h(t,x). Reescrevendo a relação de recorrência (85) para h(t,x), temos

Substituindo esses resultados na Eq. (85) e expandindo ambos os lados em série de Taylor até segunda ordem, obtemos

Considerando que a probabilidade total satisfaz p+q = 1, a equação acima se reduz para

que representa uma equação generalizada para a caminhada aleatória. Por se tratar de uma equação diferencial do tipo hiperbólica, devemos esperar que sua solução seja válida também no regime de pequenos tempos, já que ela incorpora naturalmente, uma derivada segunda com respeito ao tempo na função h(x,t). Como veremos na seção seguinte, esse fato é de fundamental importância para corrigir as inconsistências presentes na descrição de Einstein.

Algumas aproximações interessantes devem ser discutidas na Eq. (90). Primeiramente, observamos que a conexão direta com o movimento browniano difusivo é estabelecida quando assumimos que p = q = 1/2. Neste caso, definindo

a Eq. (90) se reduz a

que é precisamente a Eq. (7). Novamente, a equação de difusão que serviu de base para o tratamento de Einstein é recuperada quando fazemos o mesmo tipo de aproximação (veja a Eq. (9)), ou seja,

Portanto, a conexão com o contínuo é estabelecida de maneira consistente, de modo que todo o tratamento posterior, em particular, o cálculo dos valores médios das grandezas fisicamente relevantes, permanece idêntico ao das seções 2 e 3.

6. A Equação de difusão generalizada

Como vimos, a equação comumente utilizada para descrever transmissão de calor e difusão de partículas, constitui na verdade, um modelo aproximado, ou seja, uma descrição menos rigorosa de tais fenômenos. Um argumento favorável a essa visão se baseia na idéia de que equações parabólicas do tipo da Eq. (10) transmitem (em alguns regimes) sinais com velocidades infinitas. Naturalmente, tal resultado é inconsistente já que a velocidade máxima com a qual uma perturbação se propaga num fluido ou meio elástico deve ser da ordem da velocidade do som.

Se considerarmos que em cada intervalo de tempo t uma partícula se desloca aleatoriamente com velocidade v = l/t, vemos que a Eq. (92) pode ser reescrita como

que representa uma equação de onda amortecida para a caminhada aleatória.

Para estudar a influência do termo adicional na equação de movimento, vamos considerar uma onda plana se deslocando num meio infinito. Em x = 0 supomos que h(0,t) = ei w t, onde w é a freqüência de vibração da onda. Escrevendo a solução geral de (94) na forma

onde A e B são constantes, obtemos

e

com a velocidade de propagação da onda escrita como

Para o caso em que w << v2/D, ou equivalentemente, ¶2h/¶t2 << ¶h/¶t as relações anteriores se reduzem a

que são os resultados obtidos da equação de difusão usual. Por outro lado, para o caso em que w >> v2/D, os resultados são também fsicamente consistentes, pois a velocidade de propagação da onda tem como limite a velocidade das partículas. De fato, a freqüência de vibração de uma onda se deslocando num meio difusivo não deve exceder a freqüência de colisão das partículas do meio.

A solução da Eq. (94) para as condições gerais h(x,0) = Nd(x) e (¶h/¶t)t = 0 = 0, válida para |x | < vt pode ser escrita como [27]

sendo J0 e J1 funções de Bessel de primeira espécie.

Para o caso | x | > vt, a solução de d'Alembert para uma onda plana amortecida se deslocando na direção x é recuperada

Da expressão acima vemos também que a velocidade de propagação da onda nunca excede a velocidade das partículas. Como o produto vt é da ordem do livre caminho médio l, o argumento das funções J0 e J1 cresce rapidamente quando | x | é muito menor que vt. Neste caso, a expansão assintótica para as funções de Bessel fornecem

Portanto, a Eq. (101) pode ser reescrita como

que representa a solução geral da equação de onda modificada para a caminhada aleatória, sendo y = x/vt < 1. Note que a expressão acima é composta de duas partes. O primeiro termo relaciona a solução de onda de d'Alembert que rapidamente se torna desprezível, enquanto que o segundo se refere a difusão das partículas.

No limite y << 1, ou equivalentemente x << vt (tempos longos), o segundo termo da solução acima tende para

que é precisamente a solução da equação de difusão usual (Cf. Eq. (24)). Note que D foi reintroduzido pela definição (91).

Concluindo, vemos que a equação ondulatória hiperbólica (7), ou equivalentemente (94), resolve o problema difusivo para tempos curtos, cuja existência foi reconhecida pelo próprio Einstein ao propor sua teoria do MB. Nesse aspecto, é importante ressaltar que muitos livros-texto que tratam o problema difusivo não discutem o problema de tempos curtos, ou equivalentemente, se a propagação de uma perturbação com velocidade infinita num meio contínuo é conceitualmente correta.

Agradecimentos

Os autores agradecem a Mário José de Oliveira por sua leitura e comentários na versão preliminar do manuscrito. Este trabalho foi parcialmente financiado pelo CNPq e CAPES.

Referências

[1] R. Brown, Phil. Mag. 4, 161 (1828).

[2] A. Einstein, Ann. d. Phys. 17, 549 (1905). Uma versão para o português encontra-se na Ref. [21].

[3] J. Perrin, Ann. Chem. Phys. 18, 1 (1909); veja também Les Atomes (Libraire Alcan, Paris, 1914), 4 ed. e Atoms (Van Nostrand, Nova York, 1916).

[4] M.V. Smoluchowski, ADP 21, 756 (1906); Phys. Zeits 17, 557 (1916).

[5] P. Langevin, Compt. Rend. 146, 530 (1908).

[6] A. Fokker, Tese de doutorado, Leiden, (1913).

[7] H.C. Burger, Versl. Acad. Amst. 25, 1482 (1917).

[8] R. Fürther, Ann. d. Physik 53, 177 (1917).

[9] L.S. Ornstein, Versl. Acad. Amst. 26, 1005 (1917).

[10] M. Planck, Berl. Ber., 324 (1927).

[11] M. Kac, Am. Math. Month. 54, 369 (1947).

[12] S. Chandrasekhar, Rev. Mod. Phys. 15, 1 (1943).

[13] P. Chatterjee, L. Hernquist and A. Loeb, Phys. Rev. Lett. 88, 121103 (2002).

[14] Veja, por exemplo, J. Zhang and L. Hui, astro-ph/0508384 e referências citadas.

[15] R. Morgado and F. Oliveira, Phys. Rev. Lett. 89, 100601 (2001).

[16] J. Mathews and R. Walker, Mathematical Methods of Physics (Addison Wesley, Nova Iorque, 1973), 2 ed.

[17] H. Risken, The Fokker-Planck Equation: Methods of Solutions and Applications (Springer, Nova Iorque, 1996), Series in Synergetics, v. 18.

[18] F. Reif, Fundamentals of Statistical and Thermal Physics (McGraw-Hill, Nova Iorque, 1965).

[19] L.D. Landau and E.M. Lifshitz, Fluid Mechanics (Pergamon Press, Nova Iorque, 1959).

[20] A. Pais, Sutil é o Senhor: Vida e Pensamento de Albert Einstein (Editora Nova Fronteira, Rio de Janeiro, 1982).

[21] J. Stachel, O Ano Miraculoso de Einstein: Cinco Artigos que Mudaram a Face da Física (Editora UFRJ, Rio de Janeiro, 2001).

[22] T. Tomé e M.J. Oliveira, Dinâmica Estocástica e Irreversibilidade (Editora da Universidade de São Paulo, São Paulo, 2001).

[23] Y.B. Zel'dovich, A.A. Ruzmaikin, and D.D. Sokoloff, The Almighty Chance (World Scientific, Singapore, 1990).

[24] N.G.V. Kamper, Stochastic Process in Physics and Chemistry (North-Holland, Amsterdan, 1992).

[25] G.E. Uhlenbeck and L.S. Ornstein, Phys. Rev. 36, 823 (1930).

[26] S.R.A. Salinas, Introdução à Física Estatística (Editora da Universidade de São Paulo, São Paulo, 1999).

[27] H.D. Weymann, Am. J. Phys. 35, 488 (1965).

Recebido em 31/8/2006; Aceito em 26/9/2006

  • [1] R. Brown, Phil. Mag. 4, 161 (1828).
  • [2] A. Einstein, Ann. d. Phys. 17, 549 (1905). Uma versăo para o portuguęs encontra-se na Ref. [21].
  • [3] J. Perrin, Ann. Chem. Phys. 18, 1 (1909);
  • veja também Les Atomes (Libraire Alcan, Paris, 1914), 4 ed.
  • e Atoms (Van Nostrand, Nova York, 1916).
  • [4] M.V. Smoluchowski, ADP 21, 756 (1906);
  • Phys. Zeits 17, 557 (1916).
  • [5] P. Langevin, Compt. Rend. 146, 530 (1908).
  • [6] A. Fokker, Tese de doutorado, Leiden, (1913).
  • [7] H.C. Burger, Versl. Acad. Amst. 25, 1482 (1917).
  • [8] R. Fürther, Ann. d. Physik 53, 177 (1917).
  • [9] L.S. Ornstein, Versl. Acad. Amst. 26, 1005 (1917).
  • [10] M. Planck, Berl. Ber., 324 (1927).
  • [11] M. Kac, Am. Math. Month. 54, 369 (1947).
  • [12] S. Chandrasekhar, Rev. Mod. Phys. 15, 1 (1943).
  • [13] P. Chatterjee, L. Hernquist and A. Loeb, Phys. Rev. Lett. 88, 121103 (2002).
  • [14] Veja, por exemplo, J. Zhang and L. Hui, astro-ph/0508384 e referęncias citadas.
  • [15] R. Morgado and F. Oliveira, Phys. Rev. Lett. 89, 100601 (2001).
  • [16] J. Mathews and R. Walker, Mathematical Methods of Physics (Addison Wesley, Nova Iorque, 1973), 2 ed.
  • [17] H. Risken, The Fokker-Planck Equation: Methods of Solutions and Applications (Springer, Nova Iorque, 1996), Series in Synergetics, v. 18.
  • [18] F. Reif, Fundamentals of Statistical and Thermal Physics (McGraw-Hill, Nova Iorque, 1965).
  • [19] L.D. Landau and E.M. Lifshitz, Fluid Mechanics (Pergamon Press, Nova Iorque, 1959).
  • [20] A. Pais, Sutil é o Senhor: Vida e Pensamento de Albert Einstein (Editora Nova Fronteira, Rio de Janeiro, 1982).
  • [21] J. Stachel, O Ano Miraculoso de Einstein: Cinco Artigos que Mudaram a Face da Física (Editora UFRJ, Rio de Janeiro, 2001).
  • [22] T. Tomé e M.J. Oliveira, Dinâmica Estocástica e Irreversibilidade (Editora da Universidade de Săo Paulo, Săo Paulo, 2001).
  • [23] Y.B. Zel'dovich, A.A. Ruzmaikin, and D.D. Sokoloff, The Almighty Chance (World Scientific, Singapore, 1990).
  • [24] N.G.V. Kamper, Stochastic Process in Physics and Chemistry (North-Holland, Amsterdan, 1992).
  • [25] G.E. Uhlenbeck and L.S. Ornstein, Phys. Rev. 36, 823 (1930).
  • [26] S.R.A. Salinas, Introduçăo ŕ Física Estatística (Editora da Universidade de Săo Paulo, Săo Paulo, 1999).
  • [27] H.D. Weymann, Am. J. Phys. 35, 488 (1965).
  • 1
    E-mail:
  • 2
    Einstein obteve seus resultados fazendo a expansão no tempo somente até primeira ordem [2]. Por razões que serão discutidas adiante, consideraremos termos até segunda ordem em t na expansão da funcão h(
    x,t).
  • 3
    Uma classificacão das equacões diferenciais parciais pode ser vista na Ref. [16].
  • 4
    O integrando de (28) é composto pelo produto de uma funcão par por uma funcão ímpar.
  • 5
    Note que a
    dx =
    para
    n par.
  • 6
    O ruído é branco ("white noise") se o espectro de potência
    S(w) da funcão correlacão < x(
    t)x(0) > é independente da frequência, sendo
    S(w) =
    < x(
    t)x(0) >
    dt. No tratamento de Langevin, < x(
    t)x(0) > = Gd(
    t), temos
    S(w) = G.
  • 7
    Nos chamados processos markovianos não existem efeitos de memória, ou seja, a probabilidade condicional relativa a cada variável aleatória
    qi(
    t) de uma partícula, só depende do valor de
    qi =
    num instante anterior
    t
    0 [17,18].
  • Datas de Publicação

    • Publicação nesta coleção
      15 Jun 2007
    • Data do Fascículo
      2007

    Histórico

    • Aceito
      26 Set 2006
    • Recebido
      31 Ago 2006
    Sociedade Brasileira de Física Caixa Postal 66328, 05389-970 São Paulo SP - Brazil - São Paulo - SP - Brazil
    E-mail: marcio@sbfisica.org.br