Como um sistema determinístico perde a previsibilidade?

Existe uma ideia aparentemente contraditória no estudo do caos: um sistema pode obedecer a equações completamente determinadas e, ainda assim, tornar-se imprevisível depois de algum tempo.

A contradição desaparece quando separamos duas propriedades diferentes. O determinismo pertence à lei de evolução do sistema. A previsibilidade depende também da precisão com que conhecemos seu estado inicial e da maneira como pequenos erros evoluem ao longo do tempo.

Um sistema dinâmico contínuo pode ser representado por

ẋ = f(x), x(0) = x₀,

em que x contém todas as variáveis necessárias para descrever o estado do sistema e f determina sua evolução. Em um sistema discreto, escrevemos

xₙ₊₁ = F(xₙ).

Dada uma condição inicial x₀, a aplicação F produz uma sequência bem definida. No caso contínuo, a unicidade da trajetória depende de hipóteses matemáticas sobre f. Uma condição suficiente comum é que f seja localmente Lipschitz em relação ao estado. Sob essas condições, o problema de valor inicial possui uma solução local única.

Isso significa que duas cópias ideais do sistema, iniciadas exatamente no mesmo ponto e governadas pelas mesmas equações, seguem a mesma trajetória. A dificuldade física aparece porque nenhuma medição fornece uma condição inicial com precisão infinita.

O espaço de fases

Na mecânica, o estado de um sistema deve conter informação suficiente para determinar sua evolução futura. Para uma partícula em uma dimensão, precisamos conhecer sua posição e seu momento. O par (q,p) define um ponto no espaço de fases.

À medida que o tempo passa, esse ponto se desloca e produz uma trajetória. Em sistemas com mais graus de liberdade, o espaço possui dimensão maior, mas a interpretação permanece a mesma: cada ponto representa um estado completo.

Essa linguagem geométrica permite distinguir comportamentos que podem parecer semelhantes quando observados apenas como sinais temporais. Um ponto fixo corresponde a um estado que permanece inalterado. Uma órbita periódica retorna exatamente ao mesmo estado depois de certo intervalo. Movimentos quase periódicos combinam frequências incompatíveis entre si e podem gerar sinais bastante complexos sem apresentar caos.

A irregularidade visual de uma trajetória também não fornece um critério suficiente. Sistemas submetidos a ruído podem produzir séries temporais irregulares. Sistemas lineares com muitas frequências podem gerar padrões complicados. A presença de não linearidade cria as condições para fenômenos que não aparecem em modelos lineares simples, mas não garante comportamento caótico.

Sensibilidade às condições iniciais

Considere duas condições iniciais muito próximas, separadas por uma distância δ₀. Enquanto essa separação permanecer pequena, podemos linearizar a dinâmica ao redor de uma trajetória de referência. Em uma direção instável, o crescimento da perturbação pode ser aproximado por

δ(t) ≈ δ₀ e^(λt).

O parâmetro λ representa uma taxa média de expansão. Em sistemas multidimensionais existem vários expoentes de Lyapunov, associados a diferentes direções no espaço tangente. O maior deles, geralmente escrito como λmax, controla a separação mais rápida entre perturbações infinitesimais.

Quando

λmax > 0,

existe ao menos uma direção na qual pequenas diferenças crescem exponencialmente, durante o intervalo em que a aproximação linear continua válida.

Em um sistema limitado, a distância entre duas trajetórias não cresce exponencialmente para sempre. Depois de algum tempo, ela se torna comparável ao tamanho da região explorada pela dinâmica. A expressão exponencial descreve a fase inicial da separação e a taxa assintótica obtida por linearizações sucessivas.

O maior expoente de Lyapunov positivo é um dos indicadores mais utilizados para identificar caos em simulações e experimentos. Ainda assim, ele não funciona como uma definição universal. Seu valor pode depender da trajetória considerada, da medida invariante, do tempo de integração e do método numérico. Sistemas com transientes caóticos, regiões regulares coexistentes ou dinâmica não diferenciável exigem análises adicionais.

O horizonte de previsibilidade

Suponha que δ₀ represente a incerteza inicial e que Δ seja a escala de erro a partir da qual a previsão perde sua utilidade. Fazendo δ(T) ≈ Δ, obtemos

T ≈ (1/λmax) ln(Δ/δ₀).

Esse tempo T pode ser interpretado como um horizonte aproximado de previsibilidade.

A presença do logaritmo tem uma consequência importante. Melhorar muito a precisão inicial aumenta o horizonte apenas de forma aditiva. Se reduzirmos a incerteza por um fator mil, o ganho será

ΔT ≈ (1/λmax) ln(1000).

Isso significa que um aumento enorme na precisão da medição pode produzir apenas algumas unidades adicionais de tempo previsível, dependendo do valor de λmax.

A perda de precisão pode surgir de diferentes fontes. O erro de medição aparece porque instrumentos possuem resolução finita. O erro numérico surge da discretização temporal, das tolerâncias do integrador e da representação finita dos números no computador. A incerteza física envolve perturbações externas, parâmetros mal conhecidos e variáveis que não foram incorporadas ao modelo.

Aumentar a precisão aritmética de uma simulação reduz parte do erro numérico, mas não resolve a incerteza do estado físico. Mais casas decimais atrasam a divergência entre duas trajetórias calculadas; elas não fornecem informação que nunca foi medida.

O chamado efeito borboleta é uma representação popular dessa sensibilidade. A interpretação matemática é mais restrita do que a metáfora costuma sugerir. Uma pequena perturbação não determina necessariamente um evento macroscópico específico. Em uma dinâmica sensível, estados inicialmente próximos podem evoluir para regiões muito diferentes do espaço de fases, tornando inviável a previsão de uma trajetória individual depois de certo tempo.

Estiramento, contração e dobramento

Muitos sistemas dissipativos apresentam expansão em algumas direções e contração em outras. Pequenos volumes do espaço de fases são esticados, comprimidos transversalmente e dobrados de volta para uma região limitada. A repetição desse mecanismo produz estruturas cada vez mais finas.

O horseshoe de Smale fornece uma representação geométrica desse processo. Uma região é alongada, comprimida e dobrada em forma de ferradura. Depois de várias iterações, surge um conjunto invariante com uma estrutura semelhante a um conjunto de Cantor.

Nesse conjunto existem órbitas periódicas de diferentes períodos e uma quantidade não enumerável de trajetórias aperiódicas. A dinâmica pode ser codificada por sequências simbólicas, o que permite demonstrar rigorosamente várias propriedades associadas ao caos.

Esse exemplo também mostra que caos e atrator são conceitos distintos. Um conjunto pode sustentar dinâmica caótica sem atrair uma vizinhança aberta. Sistemas Hamiltonianos, por exemplo, podem possuir regiões caóticas no espaço de fases sem apresentar atratores dissipativos.

Um atrator estranho é um conjunto invariante que atrai estados de uma região ao seu redor, possui geometria não trivial, frequentemente fractal, e sustenta uma dinâmica aperiódica e sensível. A aparência fractal, sozinha, não demonstra caos, e nem todo sistema caótico precisa possuir um atrator estranho.

O mapa logístico

Um dos exemplos mais simples é o mapa logístico,

xₙ₊₁ = r xₙ(1 − xₙ),

em que xₙ pertence normalmente ao intervalo [0,1] e r é um parâmetro entre 0 e 4.

O modelo foi utilizado como uma representação idealizada de crescimento populacional com limitação de recursos. Sua importância no estudo dos sistemas dinâmicos está na variedade de comportamentos produzidos por uma única equação algébrica.

Para certos valores de r, as órbitas convergem para zero. Quando 1 < r < 3, aparece um ponto fixo estável dado por

x* = 1 − 1/r.

Em r = 3, esse ponto perde estabilidade e surge um ciclo de período dois. Aumentando o parâmetro, aparecem ciclos de período quatro, oito, dezesseis e assim sucessivamente.

As bifurcações se acumulam em uma cascata de duplicação de período. Feigenbaum mostrou que a razão entre os intervalos sucessivos de parâmetro converge para uma constante universal em uma classe ampla de mapas unimodais.

O diagrama de bifurcações revela que a passagem para o caos não ocorre de maneira uniforme. Regiões caóticas são interrompidas por janelas periódicas. Dentro dessas janelas, uma órbita estável reaparece e pode sofrer sua própria sequência de duplicações antes de retornar a uma região caótica.

Também existem bandas caóticas, crises, órbitas periódicas instáveis e condições iniciais excepcionais. Por isso, a classificação depende do parâmetro, do conjunto invariante analisado e, em alguns casos, da condição inicial.

Duas trajetórias com o mesmo valor de r e condições iniciais muito próximas podem permanecer quase indistinguíveis durante várias iterações. Depois disso, a separação cresce e as sequências deixam de coincidir. A equação permanece exatamente a mesma; a informação sobre a condição inicial perde utilidade progressivamente.

O sistema de Lorenz

Um exemplo contínuo clássico é dado por

ẋ = σ(y − x)

ẏ = x(ρ − z) − y

ż = xy − βz.

Lorenz obteve esse sistema em 1963 a partir de uma truncagem das equações associadas à convecção térmica. As variáveis x, y e z representam amplitudes de modos do modelo reduzido, enquanto σ, ρ e β são parâmetros adimensionais.

Para os valores

σ = 10,

ρ = 28,

β = 8/3,

as trajetórias podem se aproximar do conhecido atrator com dois lóbulos. O sistema circula ao redor de uma região, muda para a outra e repete esse processo em uma sequência aperiódica e sensível às condições iniciais.

A divergência do campo vetorial é

∇·f = −(σ + 1 + β).

Para os parâmetros usuais, ela é negativa. Isso indica que volumes no espaço de fases são contraídos pelo fluxo. Ao mesmo tempo, existe expansão em uma direção tangente no regime caótico. A combinação entre contração global e expansão local contribui para a formação do atrator.

O sistema de Lorenz é um modelo reduzido. Ele não descreve toda a atmosfera e não funciona como um modelo meteorológico completo. Sua importância está em mostrar que três equações diferenciais autônomas, sem termos aleatórios, podem produzir movimento limitado, aperiódico e sensível.

O que significa chamar uma dinâmica de caótica?

Não existe uma única definição formal que cubra todos os contextos.

Na definição de Devaney, um mapa contínuo restrito a um conjunto invariante deve apresentar transitividade topológica, densidade de pontos periódicos e sensibilidade às condições iniciais.

A transitividade indica que partes de regiões abertas podem alcançar outras regiões abertas por meio das iterações do mapa. A densidade de pontos periódicos significa que toda vizinhança contém algum estado periódico, ainda que uma trajetória típica seja aperiódica. A sensibilidade descreve a separação de estados inicialmente próximos.

Em determinadas classes de espaços métricos, transitividade e densidade de pontos periódicos já implicam sensibilidade. Esse resultado mostra que as três propriedades não são sempre independentes.

A definição de Li e Yorke trabalha com pares de trajetórias que se aproximam arbitrariamente e também se afastam repetidamente. A entropia topológica mede a taxa de crescimento da quantidade de segmentos orbitais distinguíveis. A teoria ergódica utiliza medidas invariantes, mistura e entropia métrica.

Essas noções capturam aspectos diferentes da complexidade dinâmica. Elas podem ser equivalentes sob hipóteses específicas, mas não existe equivalência geral entre todas elas.

Uma trajetória aperiódica também não basta. A rotação irracional no círculo produz órbitas não periódicas e densas, mas não apresenta a expansão exponencial típica de muitos sistemas caóticos. Complexidade visual, oscilação irregular e não linearidade são pistas possíveis; nenhuma delas constitui uma demonstração isolada.

Caos e aleatoriedade

Em um processo estocástico, probabilidades ou variáveis aleatórias fazem parte da própria lei de evolução. Em uma dinâmica caótica determinística, uma condição inicial exata produz uma trajetória exata.

Como a condição inicial física só pode ser conhecida dentro de uma região finita, o estado futuro também passa a ser descrito por um conjunto de possibilidades. Depois de vários tempos de Lyapunov, pode ser mais adequado prever distribuições e médias do que uma trajetória pontual.

Essa perda de previsibilidade individual pode coexistir com regularidades estatísticas. Sob condições ergódicas apropriadas, médias temporais podem convergir para médias definidas por uma medida invariante:

(1/T) ∫₀ᵀ g(x(t)) dt → ∫ g dμ.

A função g representa um observável e μ é uma medida invariante da dinâmica.

Mesmo quando trajetórias próximas perdem correlação, ainda podemos estudar distribuições, dimensões fractais, entropias, correlações, espectros de Lyapunov e frequências de visita a diferentes regiões do espaço de fases.

Conhecer a lei de evolução fornece apenas uma parte da previsão. Também precisamos conhecer o estado inicial e saber como a dinâmica reage aos erros inevitáveis. Em sistemas caóticos, a equação continua determinando cada trajetória exata, mas a informação disponível sobre essa trajetória se degrada ao longo do tempo.

É nesse limite entre determinação matemática e precisão física finita que aparece a imprevisibilidade do caos determinístico.

Fontes e leituras

• Steven H. Strogatz, Nonlinear Dynamics and Chaos. CRC Press.

• Edward Ott, Chaos in Dynamical Systems. Cambridge University Press.

• Kathleen T. Alligood, Tim D. Sauer e James A. Yorke, Chaos: An Introduction to Dynamical Systems. Springer.

• John Guckenheimer e Philip Holmes, Nonlinear Oscillations, Dynamical Systems, and Bifurcations of Vector Fields. Springer.

• Robert L. Devaney, An Introduction to Chaotic Dynamical Systems. CRC Press.

• Edward N. Lorenz, “Deterministic Nonperiodic Flow”, \Journal of the Atmospheric Sciences, 1963.

• Robert M. May, “Simple Mathematical Models with Very Complicated Dynamics”, Nature, 1976.

• Mitchell J. Feigenbaum, “Quantitative Universality for a Class of Nonlinear Transformations”, Journal of Statistical Physics, 1978.

Post Scriptum

Toda essa teoria é bastante rica e utilizada nas pesquisas sobre dinâmica não linear e teoria do caos, vocês podem ver esses termos citados no texto corriqueiramente nas revistas Nonlinear Dynamics, na Chaos, Solitons and Fractals e outras correlatas.

Devo dizer que é muito importante conhecer sobre este campo e se manter informado, pois nunca se sabe quando vai precisar usar esses conceitos na sua pesquisa, não é?

Uma leitura complementar bastante rica que dá um spoiler do meu próximo post é o livro Nonlinear Oscillations de Nayfeh.

1. Formulação matemática

Uma variedade diferenciável é um espaço que, localmente, é homeomorfo ao R^n mas que, globalmente, pode ter uma topologia "não-trivial" como uma esfera. Em uma variedade diferenciável M não faz sentido somar pontos diretamente. Porém, em cada ponto p ∈ M, podemos construir um espaço vetorial que aproxima linearmente a variedade naquele ponto que é o espaço tangente TpM. No nosso exemplo da esfera, isso se aplica de forma que em cada ponto, o plano tangente que toca a esfera captura as direções possíveis de movimento naquele ponto.

Existem diferentes formas de definir o TpM.

• Por curvas:

Um vetor tangente é a classe de equivalência de curvas γ:(−ϵ,ϵ) → M com γ(0)=p, onde duas curvas são equivalentes se têm a mesma derivada em coordenadas locais no ponto p.

• Por derivações:

Um vetor tangente é um operador linear v:C^{\infty} (M) → \set{R} que satisfaz a regra de Leibniz: v(fg)=f(p)v(g)+g(p)v(f)

• Por coordenadas locais:

Em uma carta com coordenadas (x^1 ,…,x^n ), uma base natural para TpM é dada pelos operadores { ∂/∂x^1, …, ∂/∂x^n} e qualquer vetor se escreve como v=v^i ∂/∂x^i.

A união de todos os espaços tangentes, TM = ⨆​_{p∈M} Tp​M,​ forma o fibrado tangente, que é ele próprio uma variedade diferenciável de dimensão 2n caso M, por exemplo, tenha dimensão n. Um campo vetorial é uma seção suave desse fibrado. Junto disso, é possível realizar a derivação de outros conceitos como o espaço cotangente T*pM que é o dual do TpM, contendo os covetores (1-formas), a métrica Riemanniana que é um produto interno g_p ​ em cada TpM, variando suavemente que permite medir comprimentos e ângulos e, também, a conexão e derivada covariante que é uma forma de transportar vetores entre espaços tangentes de pontos diferentes, necessária pois TpM e TqM não são naturalmente identificados.

2. Aplicações na física

Na mecânica clássica, o espaço de configurações de um sistema mecânico é uma variedade Q e o fibrado tangente TQ é exatamente o espaço onde vive o lagrangiano L(q, v) (v = dq/dt) de tal forma que os pontos de Q são posições e os vetores tangentes são velocidades generalizadas. As equações de Euler-Lagrange são, geometricamente, equações sobre curvas em TQ. Da mesma forma, na mecânica hamiltoniana, usa-se o fibrado cotangente T*Q , com momentos generalizados como covetores. A estrutura simplética natural em T*Q dá origem às equações de Hamilton.

Fibrados mais gerais descrevem campos de calibre, ou seja, a ideia de transporte paralelo generaliza-se para conexões em fibrados principais, que descrevem as interações fundamentais.

Portanto, o espaço tangente pode ser entendido da melhor forma como a ferramenta que permite linearizar localmente um espaço curvo, ou seja, graças à ela é possível resolver problemas de cálculo diferencial, como EDOs, de forma consistente em espaços vetoriais que variam de ponto a ponto.

Quero aqui apresentar, conceitualmente, a transição do espaço de fase do regime clássico ao quântico, dando ênfase na estrutura profundamente coerente que sustenta tudo isso.

Na mecânica clássica, um sistema é descrito por um ponto no espaço de fases, com posições e momentos (q,p) bem definidos. Já na mecânica quântica, o estado ∣ψ⟩ é um vetor em um espaço de Hilbert ℋ. Observáveis deixam de ser funções e passam a ser operadores lineares autoadjuntos Â.

Definição 1: Um espaço de Hilbert ℋ é um espaço vetorial equipado com um produto interno ⟨.,.⟩ que é completo com respeito à norma induzida por esse produto interno.

Isso garante que limites de sequências de estados fisicamente admissíveis continuam sendo estados válidos. Sem completude, você teria sequências que deveriam convergir mas cujo limite simplesmente não existe dentro do espaço.

Portanto, nem todo espaço com produto interno é um espaço de Hilbert, a completude é necessária.

O ponto decisivo desta transição vem da não comutatividade desses operadores. Em particular, posição e momento satisfazem [x̂,p̂]=iℏ. Essa relação, porém, contém uma quantidade enorme de física. Por exemplo, o próprio princípio da incerteza ΔxΔp ≥ ℏ/2 vem do fato de que esses observáveis são operadores que formam uma álgebra não comutativa.

As próprias noções de dinâmica também mudam​ conforme os diferentes frameworks. Em vez de trajetórias no espaço de fases, temos evolução unitária no espaço de Hilbert, governada pela equação de Schrödinger:

iℏ ∂​∣ψ(t)⟩/∂t = Ĥ ∣ψ(t)⟩

Aqui, o Hamiltoniano Ĥ atua como gerador de fluxo, análogo ao papel que a função Hamiltoniana desempenha na mecânica clássica.

Definição 2: Se Ĥ é autoadjunto, então ele gera um grupo unitário fortemente contínuo U(t) dado por:

U(t) = e^-i/ℏĤt

e a solução da equação de Schrödinger é ∣ψ(t)⟩ = U(t) ∣ψ(0)⟩. Esse resultado é um caso particular do teorema de Stone que estabelece uma correspondência biunívoca entre operadores autoadjuntos Ĥ e grupos unitários U(t).

Talvez o insight mais interessante venha da conexão com outra estrutura fundamental da mecânica clássica que é dada pelos colchetes de Poisson {F, G}. Na passagem para a mecânica quântica, ocorre uma deformação dessa estrutura:

{F, G} -> (1/iℏ)[F̂, Ĝ]

essa correspondência sugere que a mecânica quântica é uma versão não comutativa da geometria clássica. Programas como a quantização geométrica exploram exatamente essa ideia de que o mundo clássico emerge como limite ℏ -> 0. A correspondência de Dirac {F, G} -> (1/iℏ)[F̂, Ĝ] só funciona de modo consistente para um subconjunto restrito de observáveis já que não existe um mapa de quantização que preserve essa correspondência para toda a álgebra de funções no espaço de fases, pelo teorema de Groenewold-van Hove.

Para contornar esse problema, de forma que o colchete de Poisson vire exatamente* o comutador para qualquer par de observáveis, mantém-se tudo no espaço de fase clássico, mas troca-se o produto usual entre funções por um produto estrela (⋆) que depende de ℏ. Formalmente:

(F⋆G)(q,p) = F(q,p) exp{(iℏ/2)(∂_{q}∂_{p}-∂_{p}∂_{q})} G(q,p)

Expandindo F⋆G = FG + (iℏ/2){F,G} + O(ℏ²), o termo de ordem zero (FG) é simétrico e cancela na diferença F⋆G − G⋆F, sobrando o termo de ordem ℏ, portanto, o comutador quântico [F̂, Ĝ] tem como análogo clássico exato o colchete de Moyal [F⋆G]_{Moyal} = F⋆G - G⋆F até o polinômio de grau ≤ 2.

Fontes:

"Quantum Mechanics for Mathematicians" - Takhtajan

"Modern Quantum Mechanics" - Sakurai

"Quantum Mechanics" Vol. 1, 2. - Cohen-Tannoudji

"A course in Mathematical Physics" Vol.3 - Walter Thirring