Ecuación de Langevin

En física, la ecuación de Langevin (llamada así por Paul Langevin) es una ecuación diferencial estocástica que describe la evolución temporal de un subconjunto de los grados de libertad. Estos grados de libertad son típicamente variables colectivas (macroscópicas) que cambian solo lentamente en comparación con las otras variables (microscópicas) del sistema. Las variables rápidas (microscópicas) son responsables de la naturaleza estocástica de la ecuación de Langevin.

Movimiento browniano como prototipo 
La ecuación original de Langevin describe el movimiento browniano, el movimiento aparentemente aleatorio de una partícula en un fluido debido a las colisiones con las moléculas del fluido:

 

El grado de libertad de interés aquí es la posición    de la partícula,  denota la masa de la partícula. La fuerza que actúa sobre la partícula se escribe como la suma de una fuerza viscosa proporcional a la velocidad de la partícula (ley de Stokes), y un término de ruido   el nombre dado en contextos físicos a términos en ecuaciones diferenciales estocásticas (procesos estocásticos) que representan el efecto de las colisiones con las moléculas del fluido. La fuerza   tiene una distribución de probabilidad gaussiana con función de correlación 

 

donde es la constante de Boltzmann,  es la temperatura y es el componente i del vector  . La forma de la δ-función corre de las correlaciones en el tiempo significa que se supone que la fuerza en un momento  no está correlacionada completamente con ella en ningún otro momento. Esto es una aproximación; la fuerza aleatoria real tiene un tiempo de correlación distinto de cero correspondiente al tiempo de colisión de las moléculas. Sin embargo, la ecuación de Langevin se usa para describir el movimiento de una partícula "macroscópica" en una escala de tiempo mucho más larga, y en este límite la correlación y la ecuación de Langevin se vuelven exactas. 

Otra característica prototípica de la ecuación de Langevin es la ocurrencia del coeficiente de amortiguamiento  en la función de correlación de la fuerza aleatoria, un hecho también conocido como relación de Einstein.

Aspectos matemáticos 
Una fuerza fluctuante correlacionada estrictamente  no es una función en el sentido matemático habitual e incluso la derivada  no está definida en este límite. La ecuación de Langevin en su forma actual requiere una interpretación basada en la Integral de Itō.

Ecuación de Langevin genérica 
Existe una derivación formal de una ecuación de Langevin genérica de la mecánica clásica. Esta ecuación genérica juega un papel central en la teoría de la dinámica crítica y en otras áreas de la mecánica estadística sin equilibrio. La ecuación para el movimiento browniano anterior es un caso especial.

Una condición esencial de la derivación es un criterio que divide los grados de libertad en categorías lentas y rápidas. Por ejemplo, el equilibrio termodinámico local en un líquido se alcanza dentro de unos pocos tiempos de colisión. Pero las densidades de cantidades conservadas, como la masa y la energía, tardan mucho más tiempo en relajarse hasta alcanzar el equilibrio. Las densidades de cantidades conservadas, y en particular sus componentes de longitud de onda larga, son, por tanto, candidatos de variable lenta. Técnicamente esta división se realiza con el operador de proyección Zwanzig, la herramienta esencial en la derivación. La derivación no es completamente rigurosa porque se basa en suposiciones (plausibles) similares a las suposiciones requeridas en otros lugares en la mecánica estadística básica.

Deje que denota las variables lentas. La ecuación de Langevin genérica entonces lee  

 

La fuerza fluctuante obedece a una distribución de probabilidad gaussiana con función de correlación 

 

Esto implica la relación de reciprocidad de Onsager  para los coeficientes de amortiguamiento . La dependencia  de  en  es despreciable en la mayoría de los casos. El símbolo \ denota el Hamiltoniano del sistema, donde  es la distribución de probabilidad de equilibrio de las variables . Finalmente, es la proyección del corchete de Poisson de las variables lentas  y  en el espacio de lento variables.

En el caso de movimiento browniano, uno tendría ,  o y. La ecuación de movimiento   para    es exacta, no hay una fuerza fluctuante  y no coeficiente de amortiguamiento

Ejemplos

Oscilador armónico en un fluido 
Un oscilador armónico no ideal se ve afectado por alguna forma de amortiguación, a partir de la cual se sigue a través del teorema de fluctuación-disipación de que debe haber algunas fluctuaciones en el sistema. El diagrama a la derecha muestra un retrato de la fase de la evolución temporal del momento, , vs. position,  de un oscilador armónico. El movimiento determinista seguiría a lo largo de las trayectorias elipsoidales que no pueden cruzarse entre sí sin cambiar la energía. La presencia de alguna forma de amortiguación, por ej. un entorno de fluido molecular (representado por los términos de difusión y amortiguación), agrega y elimina continuamente la energía cinética del sistema, lo que hace que se extienda un conjunto inicial de osciladores estocásticos (círculos de puntos), que eventualmente alcanzan el equilibrio térmico.

Ruido térmico en una resistencia eléctrica 
Existe una estrecha analogía entre la partícula Brownian paradigmática discutida anteriormente y el ruido de Johnson, el voltaje eléctrico generado por las fluctuaciones térmicas en cada resistencia. El diagrama a la derecha muestra un circuito eléctrico que consiste en una resistencia R y una capacitancia C. La variable lenta es la tensión U entre los extremos de la resistencia. El hamiltoniano lee , y la ecuación de Langevin se convierte en

 

Esta ecuación se puede utilizar para determinar la función de correlación

 

que se convierte en un ruido blanco (ruido de Johnson) cuando la capacitancia C se vuelve despreciablemente pequeña.

Dinámica crítica 
La dinámica del parámetro de orden  de una transición de fase de segundo orden se ralentiza cerca del punto crítico y se puede describir con una ecuación de Langevin. El caso más simple es la clase de universalidad "modelo A" con un parámetro de orden escalar no conservado, realizado, por ejemplo, en ferromagnetos axiales,

 

Otras clases de universalidad (la nomenclatura es "modelo A", ..., "modelo J") contienen un parámetro de orden de difusión, parámetros de orden con varios componentes, otras variables críticas y/o contribuciones de los corchetes de Poisson.

Recuperando las estadísticas de Boltzmann 
Las ecuaciones de Langevin deben reproducir la distribución de Boltzmann. El movimiento browniano sobredimensionado unidimensional es un ejemplo instructivo. El caso sobredimensionado se realiza cuando la inercia de la partícula es despreciable en comparación con la fuerza de amortiguación. La trayectoria  de la partícula en un potencial  se describe mediante la ecuación de Langevin:

 

donde el ruido se caracteriza por  y  es la constante de amortiguamiento. Nos gustaría calcular la distribución   de la posición de la partícula en el transcurso del tiempo. Una forma directa de determinar esta distribución es introducir una función de prueba , y observar el promedio de esta función en todas las realizaciones (promedio del conjunto)  

 

Si  permanece finito, entonces esta cantidad es nula. Además, utilizando la interpretación de Stratonovich, podemos deshacernos de la eta en el segundo término para terminar con

 

donde hacemos uso de la función de densidad de probabilidad . Esto se hace calculando explícitamente el promedio, 

 

donde el segundo término estaba integrado por partes (de ahí el signo negativo). Como esto es cierto para funciones arbitrarias , debemos tener: 

 

recuperando así la distribución de Boltzmann.

Técnicas equivalentes 
La solución de una ecuación de Langevin para una realización particular de la fuerza fluctuante no tiene ningún interés por sí misma, lo que sí interesa son las funciones de correlación de las variables lentas después de promediar la fuerza fluctuante. Dichas funciones de correlación también pueden determinarse con otras técnicas (equivalentes).

Ecuación de Fokker Planck 
Una ecuación de Fokker-Planck es una ecuación determinística para la densidad de probabilidad dependiente del tiempo  de las variables estocásticas . La ecuación de Fokker-Planck correspondiente a la ecuación de Langevin genérica anterior se puede derivar con técnicas estándar.

 

La distribución de equilibrio  es una solución estacionaria.

Integral de caminos 
Se puede obtener una integral de trayectoria equivalente a una ecuación de Langevin a partir de la ecuación de Fokker-Planck correspondiente o transformando la distribución de probabilidad gaussiana  de la fuerza fluctuante a una distribución de probabilidad de las variables lentas, esquemáticamente . El determinante funcional y las sutilezas matemáticas asociadas se eliminan si la ecuación de Langevin se discretiza de manera natural (causal), donde depende de  pero no en . Resulta conveniente introducir variables de respuesta auxiliares . La integral de trayectoria equivalente a la ecuación de Langevin genérica se lee  

 

donde es un factor de normalización y

 

La formulación integral del camino no agrega nada nuevo, pero sí permite el uso de herramientas de la teoría cuántica de campos; por ejemplo, métodos de grupo de perturbación y renormalización (si tienen sentido).

Véase también 

 Dinámica de Langevin

Referencias

Lectura recomendada 

 W. T. Coffey (Trinity College, Dublin, Ireland) and Yu P. Kalmykov (Université de Perpignan, France, The Langevin Equation: With Applications to Stochastic Problems in Physics, Chemistry and Electrical Engineering (Third edition), World Scientific Series in Contemporary Chemical Physics - Vol 27.
 Reif, F. Fundamentals of Statistical and Thermal Physics, McGraw Hill New York, 1965. See section 15.5 Langevin Equation
 R. Friedrich, J. Peinke and Ch. Renner.  How to Quantify Deterministic and Random Influences on the Statistics of the Foreign Exchange Market, Phys. Rev. Lett. 84, 5224 - 5227 (2000)
 L.C.G. Rogers and D. Williams. Diffusions, Markov Processes, and Martingales, Cambridge Mathematical Library, Cambridge University Press, Cambridge,  reprint of 2nd (1994) edition, 2000.

Ecuaciones diferenciales estocásticas
Mecánica estadística