UNIVERSIDAD DE GUANAJUATO
DIVISIÓN DE CIENCIAS NATURALES Y EXACTAS
MODELOS DE CÁLCULO ESTOCÁSTICO PARA EL
MOVIMIENTO BROWNIANO RELATIVISTA
T E S I S
QUE PARA OBTENER EL TÍTULO DE
LICENCIADO EN MATEMÁTICAS
P R E S E N T A
JOSÉ ULISES MÁRQUEZ URBINA
DIRECTOR DE TESIS
Dr. VÍCTOR MANUEL PÉREZ‐ABREU CARRIÓN
GUANAJUATO, GTO. AGOSTO DE 2010
Tesis presentada el día 20 de agosto de 2010.
Sinodales
Dr. Daniel Hernández Hernández
Dr. Alfredo Sandoval Villalbazo
Dr. Víctor Manuel Pérez‐Abreu Carrión
A mis padres y mi hermana...
Índice general
1. Introducción General 7
2. Preliminares de Cálculo Estocástico 15
2.1. Procesos de Wiener . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
2.2. Cálculo estocástico de Itô . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
2.2.1. Procesos estocásticos no-anticipantes . . . . . . . . . . . . . . . . . . . . . 17
2.2.2. De�nición de la integral de Itô . . . . . . . . . . . . . . . . . . . . . . . . 19
2.2.3. Las diferenciales estocásticas y la fórmula de Itô . . . . . . . . . . . . . . 23
2.2.4. Ecuaciones diferenciales estocásticas . . . . . . . . . . . . . . . . . . . . . 26
2.2.5. Ecuaciones diferenciales estocásticas y los procesos de difusión . . . . . . 30
2.2.6. Cambio de tiempo aleatorio . . . . . . . . . . . . . . . . . . . . . . . . . . 32
2.3. La integral de Stratonovich . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
2.4. La integral Backward-Itô . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
3. El Movimiento Browniano No-Relativista 41
3.1. Ecuaciones de Langevin y Fokker-Planck . . . . . . . . . . . . . . . . . . . . . . . 41
3.1.1. La ecuación de Langevin . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
3.1.2. La ecuación de Langevin y las integrales estocásticas . . . . . . . . . . . . 43
3.1.3. La ecuación de Fokker-Planck . . . . . . . . . . . . . . . . . . . . . . . . . 45
3.1.4. Teorema de �uctuación-disipación . . . . . . . . . . . . . . . . . . . . . . 46
3.1.5. El proceso de Ornstein-Uhlenbeck no-relativista . . . . . . . . . . . . . . . 48
3.2. Modelos microscópicos: el modelo de colisión binaria . . . . . . . . . . . . . . . . 51
3.2.1. El modelo de colisión binaria elástica . . . . . . . . . . . . . . . . . . . . . 53
1
4. El Movimiento Browniano Relativista en el Espacio Fase 63
4.1. Ecuaciones de Langevin y Fokker-Planck relativistas . . . . . . . . . . . . . . . . 64
4.1.1. La ecuación de Langevin relativista: principios generales de construcción . 64
4.1.2. Las ecuaciones de Fokker-Planck . . . . . . . . . . . . . . . . . . . . . . . 70
4.2. Movimientos libres en un baño isotrópico y las relaciones de Einstein . . . . . . . 74
5. Ejemplos Unidimensionales y el Desplazamiento Cuadrático Medio 79
5.1. Consideraciones generales y las ecuaciones de la energía y la velocidad . . . . . . 79
5.2. Desplazamiento cuadrático medio asintótico . . . . . . . . . . . . . . . . . . . . . 85
5.3. Ejemplos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 86
5.3.1. Amplitud constante . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 87
5.3.2. Coe�ciente de fricción constante en la ecuación de Langevin-backward . . 90
5.3.3. Coe�ciente de fricción constante en la ecuación de Langevin-Itô . . . . . . 94
6. Reparametrizaciones de la Ecuación de Langevin 99
6.1. Reparametrización en términos del tiempo propio . . . . . . . . . . . . . . . . . . 99
6.2. Observadores en movimiento . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 105
6.2.1. Transformación de Lorentz de la densidad de transición en el espacio fase 106
6.2.2. Transformación de Lorentz de la ecuación de Langevin . . . . . . . . . . . 106
7. El Modelo de Colisión Binaria Relativista 115
7.1. Cinemática de las colisiones relativistas . . . . . . . . . . . . . . . . . . . . . . . 116
7.2. La distribución del baño térmico y la fuerza de deriva . . . . . . . . . . . . . . . 118
7.3. Aproximación de Langevin . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 120
A. Procesos de Markov y Procesos de Difusión 123
B. Distribución de Maxwell 129
C. Conceptos Básicos de Relatividad 131
C.1. Notación y convenciones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 131
C.2. Transformaciones de Lorentz-Poncairé . . . . . . . . . . . . . . . . . . . . . . . . 134
2
D. Densidades de Probabilidad en Relatividad Especial 137
D.1. Invarianza de Lorentz de las densidades en el espacio fase . . . . . . . . . . . . . 137
D.2. La distribución de Jüttner . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 138
E. Notación 143
Bibliografía 145
3
Agradecimientos
En estas líneas deseo expresar mi sincero reconocimiento a todos aquellos que, directa o
indirectamente, forman parte de este logro.
Comienzo con mis padres, cuyo invaluable apoyo y dirección me llevaron a bien concluir una
carrera. ¡Gracias por todo!
Una mención muy especial merece mi asesor, el Dr. Víctor Manuel Pérez-Abreu Carrión,
por su tiempo, tolerancia (¡mucha!), entusiasmo y apoyo. También le agradezco el tema de esta
tesis, el cual me ha permitido entrar en contacto con temas e ideas que me han conmocionado.
Agradezco a mis sinodales, el Dr. Daniel Hernández Hernández y el Dr. Alfredo Sandoval
Villalvazo, por el tiempo que dedicaron a la revisión de esta tesis y por sus valiosos comentarios.
También aprecio enormemente la lectura minuciosa y los provechosos comentarios hechos por
el Dr. Constantin Tudor.
Me siento afortunado de poder haber escuchado de viva voz los comentarios del Dr. Leopol-
do García-Colín Scherer y Tomoi Koide, fueron de gran ayuda para entender las di�cultades
asociadas al movimiento Browniano relativista.
A todos mis profesores, amigos y compañeros les doy gracias por la convivencia y por
haberme compartido sus conocimientos. Aquí quiero reconocer la ayuda que recibí por parte
de mi compañero Guillermo Basulto en la elaboración de la portada de este trabajo y en la
disolución de algunos problemas con Latex.
Finalmente, agradezco a la Universidad de Guanajuato y al CIMAT por las becas y apoyos
que me otorgaron durante mis estudios de licenciatura, pero sobre todo por la eduación que me
proporcionaron.
5
Capítulo 1
Introducción General
En 1905 Albert Einstein (1879-1955) publicó cuatro artículos que, según algunos especia-
listas, cambiaron para siempre la Física y la forma como percibimos el mundo. Estos traba-
jos explicaron el efecto fotoeléctrico, establecieron los fundamentos de la teoría especial de
la relatividad y justi�caron teóricamente el fenómeno del movimiento Browniano clásico (no-
relativista). Eliminando los efectos gravitacionales, la relatividad especial ha probado ser el
marco correcto para describir procesos físicos en las escalas terrestres. Durante el siglo pasado
se hicieron grandes esfuerzos para adaptar las teorías no-relativistas como la termodinámica,
la mecánica cuántica o las teorías de campo a los requerimientos de la relatividad especial.
Siguiendo esta línea, se ha venido desarrollando lo propio para el movimiento Browniano rela-
tivista en el marco de la relatividad especial. La presente tesis analiza los aspectos matemáticos
de algunos progresos recientes hechos en este tema y que siguen el enfoque del cálculo estocástico
de Itô.
Históricamente, el término �movimiento Browniano�hace referencia a la dinámica irregular
que se observa en algunas partículas (e.g. granos de polen) que se encuentran en un medio
líquido. Este fenómeno físico, que ya había sido descrito en 1784 por el médico holandés Jan
Ingenhousz (1730-1799), fue observado en detalle por primera vez en 1827 por el botánico
inglés Robert Brown (1773-1858). Tuvieron que pasar alrededor de 80 años para que William
Sutherland (1859-1911), Albert Einstein y Marian von Smoluchowski (1872-1917) pudieran
dar una explicación física teórica de estas observaciones. Ellos propusieron que el movimiento
Browniano es causado por interacciones microscópicas cuasi-aleatorias con las partículas que
7
forman el líquido. En 1909 sus teorías fueron con�rmadas experimentalmente por Jean Baptiste
Perrin (1870-1942), lo cual proporciono evidencia adicional apoyando la estructura atómica de
la materia. La teoría física del movimiento Browniano fue elaborada más a fondo durante la
primera mitad del siglo XX por los físicos Paul Langevin (1872-1946), Adriaan Fokker (1887-
1972), Max Planck (1858-1947), Oskar Klein (1894-1977), George Eugene Uhlenbeck (1900-
1988), Leonard Salomon Ornstein (1880-1941) y Hendrik Anthony Kramers (1894-1952).
Desde el punto de vista matemático, el estudio del movimiento Browniano aparece en la
tesis doctoral de Bachelier (1870-1946) en 1900, quién propuso este proceso para estudiar la
teoría de especulación de la Bolsa de Valores en París. El estudio riguroso del movimiento Brow-
niano fue iniciado en la década de 1920 por Wiener (1894-1964). Igualmente pioneros son los
trabajos de Lévy (1886-1971), Kolmogorov (1903-1987), Doob (1910-2004) y Feller (1906-1970),
entre otros, quienes contribuyeron a las bases matemáticas de la teoría de procesos estocásticos.
Posteriormente, entre 1944 y 1970 Itô (1915-2008), Gihman (1918-1985), Fisk y Stratonovich
(1930-1997), introdujeron y caracterizaron diferentes tipos de integrales estocásticas y ecua-
ciones diferenciales estocásticas. Estos modelos son un herramienta muy útil y e�ciente para
describir procesos estocásticos que surgen en el mundo real y su estudio atrajo considerable
interés en las décadas pasadas; especialmente en el área de �nanzas [9, 27, 30] y biología [3, 59].
Cabe destacar que en el año 2006, la Unión Matemática Internacional otorgó a Itô el primer
Premio Gauss por Aplicaciones Matemáticas, como reconocimiento a que el cálculo estocástico
que él fundó es actualmente una rica, importante y exitosa rama de las matemáticas, con un
impacto formidable en la tecnología, los negocios y la vida diaria de la gente.
El tema central de este trabajo de tesis es el estudio de generalización de los modelos del
movimiento Browniano basados en ecuaciones diferenciales estocásticas al marco de la teoría
especial de la relatividad. En la literatura física, las ecuaciones diferenciales estocásticas a
menudo se conocen como ecuaciones de Langevin, sin embargo nosotros dejaremos este nombre
para un tipo de ecuaciones especí�cas: aquellas que modelan la dinámica de una partícula
Browniana. Así mismo, reservaremos el término �movimiento Browniano� para el fenómeno
físico y utilizaremos �proceso de Wiener�para el caso del movimiento Browniano en el campo
de la probabilidad.
Desde una perspectiva matemática, las ecuaciones diferenciales estocásticas determinan
8
modelos precisos y bien de�nidos de procesos estocásticos. Desde el punto de vista físico, su
utilidad para la descripción de un sistema realista es, a priori, un tema abierto. La derivación
de ecuaciones diferenciales estocásticas a partir de modelos microscópicos atrajo gran interés en
las décadas pasadas. Algunos esfuerzos en esta dirección ayudaron a clari�car la aplicabilidad
de las ecuaciones diferenciales estocásticas a problemas en física.
Cuando se intenta generalizar los conceptos del movimiento Browniano clásico a la rela-
tividad especial, varios elementos de la termodinámica relativista del equilibrio y la mecánica
estadística relativista juegan un papel importante. Más precisamente, algunos principios ter-
moestadísticos gobiernan el comportamiento estacionario de las partículas Brownianas y, así,
imponen restricciones en la estructura de la ecuación de Langevin relativista. Los primeros
artículos sobre termodinámica relativista fueron publicados en 1907 por Planck y Einstein. Un
objetivo principal de sus trabajos fue el identi�car las llamadas leyes de transformación de
Lorentz de las variables termodinámicas (temperatura, presión, étc.). Estos resultados fueron
cuestionados en 1963 por Ott [49], cuyo trabajo inició un intenso debate sobre el comportamien-
to correcto de la transformación de las cantidades termodinámicas en la relatividad especial.
Sin embargo, como clari�caron van Kampen (1921- ) [39] y Yuen [61] a �nales de la década
de los 60, la controversia que rodeaba a la termodinámica relativista puede resolverse al darse
cuenta de que las cantidades termodinámicas pueden ser de�nidas de una manera diferente e
igualmente consistente.
Mientras algunos autores consideraron a la termodinámica relativista como un teoría pura-
mente macroscópica, otros trataron de adoptar un planteamiento más fundamental al enfocarse
en la mecánica estadística relativista del equilibrio (térmico). Los trabajos pioneros en esta
última dirección fueron hechos por los estudiantes de Planck, von Mosengeil (1884-1906) [47]
y von Laue (1879-1960) [45], y su colaborador Jüttner (1878-?) [38], quién en 1911 derivó la
generalización relativista de la distribución de Maxwell para la velocidad.
El debate recurrente sobre la termoestadística relativista puede rastrearse a la di�cultad de
tratar interacciones que envuelven muchas partículas de una manera �relativistamente�consis-
tente. En la física no-relativista, las interacciones pueden propagarse a velocidades arbitraria-
mente grandes, por lo que es posible modelarlas mediante potenciales de interacción instantá-
neos los cuales entran en la función de Hamilton; a partir de este punto, la mecánica estadística
9
no-relativista surge sin mucha di�cultad. Desafortunadamente, la situación se vuelve más com-
plicada en el caso relativista: debido a su propagación acotada, las interacciones relativistas
deben ser modeladas por medio de campos que pueden intercambiar energía con las partícu-
las. Estos campos agregan un número in�nito de grados de libertad al sistema de partículas.
Así, en relatividad especial usualmente es muy complicado, o incluso imposible, desarrollar un
formalismo consistente de campo Hamiltoniano libre para las partículas interactuantes.
A pesar de las di�cultades que impiden un tratamiento riguroso de los sistemas relativistas
de muchas partículas, durante la última mitad del Siglo XX se hicieron progresos conside-
rables en la construcción de una teoría cinética relativista aproximada basada en la ecuación de
Boltzmann relativista para las densidades de probabilidad en el espacio fase de una partícula.
A partir de tal teoría cinética, solo hay un paso relativamente pequeño en la formulación de
una teoría relativista del movimiento Browniano en términos de ecuaciones de Fokker-Planck
y ecuaciones de Langevin. Mientras que la ecuación de Boltzmann relativista es una ecuación
integro-diferencial no-lineal de la densidad de probabilidad, las ecuaciones de Fokker-Planck
son ecuaciones diferenciales parciales lineales y, por lo tanto, pueden ser resueltas o analizadas
con mayor facilidad.
Esta tesis se enfoca en procesos estocásticos relativistas que están caracterizados por ecua-
ciones de evolución lineales para sus respectivas densidad de probabilidad en el espacio fase.
La correspondiente teoría fenomenológica del movimiento Browniano relativista experimentó
un progreso considerable durante la década pasada, con aplicaciones en varias áreas de la as-
trofísica [12, 8, 62, 13, 4] y la física de altas energías [55, 35, 60, 36, 52, 1, 2, 63]. Desde una
perspectiva general, los procesos estocásticos relativistas proporcionan un enfoque útil cuando
se tiene que modelar el comportamiento cuasi-aleatorio de partículas relativistas en un entorno
complejo. Por lo tanto, es de esperarse que los conceptos de movimiento Browniano relativista y
de difusión jugarán un papel cada vez más importante en investigaciones futuras de, por ejem-
plo, procesos de relajación y termalización en astrofísica o en los experimentos de colisiones a
altas energías.
De acuerdo a nuestro conocimiento, los primero estudios matemáticos detallados sobre pro-
cesos de difusión relativistas fueron realizados de manera independiente por Lopuszànski, Rud-
berg y Schay, entre 1953 y 1961. En las décadas de 1960 y 1970 su trabajo pionero fue comple-
10
mentado por Dudley quien publicó una serie de artículos [14, 15, 16, 17] que buscaban proveer
de un enfoque axiomático a los procesos de Markov que fuesen Lorentz invariantes en el espacio
fase. De manera independiente, Hakim persiguió el mismo objetivo [31, 32, 33, 34], cuyo análisis
ayudo a dilucidar las sutilezas conceptuales de los procesos estocásticos relativistas. Especí�-
camente, Dudley y Hakim probaron la inexistencia de procesos de Markov Lorentz invariantes
no triviales en el espacio-tiempo de Minkowski, como había sido sugerido por Lopuszànski.
Este resultado fundamental implica que es difícil encontrar generalizaciones relativistas de la
ecuación de difusión no-relativista@
@t% = Dr2%;
donde D > 0 es la constante de difusión y % (t;x) � 0 es la densidad de probabilidad para las
posiciones x 2 Rd de la partícula al tiempo t. Con el �n de evitar una confrontación con el
resultado de inexistencia de Dudley y Hakim, usualmente se adoptan alguna de las siguientes
estrategias:
Se consideran procesos de difusión no MarkovianosX (t) en el espacio-tiempo de Minkows-
ki.
Se construyen procesos de Markov en el espacio fase, aceptables desde el punto de vista de
la relatividad, al considerar no sólo la posición X (t) de la partícula difusiva sino también
su momento P (t).
En esta tesis se tratará exclusivamente la segunda estrategia. Algunos ejemplos típicos de
procesos de Markov relativistas en el espacio fase son procesos descritos por ecuaciones de
Fokker-Planck o por ecuaciones de Langevin. De forma similar a lo que se hace con la ecuación
de Boltzmann relativista, las ecuaciones de Fokker-Planck relativistas en el espacio fase pueden
utilizarse para describir fenómenos fuera del equilibrio térmico y fenómenos de relajación en
sistemas relativistas de muchas partículas. Las ecuaciones de Fokker-Planck pueden derivarse de
la ecuaciones de Langevin, como aproximaciones a ecuaciones lineales maestras más generales
o al aproximar las integrales de colisión en la ecuación no-lineal de Boltzmann por expresiones
diferenciales que contienen los coe�cientes de fricción efectiva y de difusión. En particular, este
último método fue aplicado exitosamente en diferentes áreas de la física en la décadas pasadas,
11
incluyendo la física de plasmas, la física de altas energías y la astrofísica. En las décadas de 1980
y 1990 este enfoque fue mejorado y se desarrollaron varios métodos numéricos para resolver la
ecuación de Fokker-Planck. Aplicaciones recientes incluyen la modelación de procesos de difusión
y termalización en plasmas quark-gluón, así como en la descripción de procesos complejos de
alta energía en astrofísica, por mencionar algunos.
Un enfoque complementario para el estudio de los procesos estocásticos relativistas en el
espacio fase parten de la ecuación de Langevin. Las ecuaciones diferenciales estocásticas del
tipo Langevin generan trayectorias muestrales explícitas para el movimiento de una partícula
Browniana relativista. Estas ecuaciones de Langevin pueden postularse como modelos feno-
menológicos u obtenerse de modelos microscópicos más precisos al imponer una serie de apro-
ximaciones. Comparado con el caso no-relativista, la última tarea se vuelve considerablemente
más complicada, debido a las di�cultades conceptuales y técnicas ya mencionadas en la formu-
lación de teorías relativistamente consistentes para sistemas de muchas partículas. El enfoque
fenomenológico de Langevin del movimiento Browniano relativista fue iniciado por Debbasch,
Mallick y Rivet [10], quienes propusieron una generalización relativista simple del proceso de
Ornstein-Uhlenbeck, que representa casos límite especiales de una clase mayor de procesos de
Langevin relativistas. Desde un punto de vista práctico, las ecuaciones de Langevin relativistas
proporcionan una herramienta útil para modelar la dinámica de partículas relativistas en un
medio aleatorio, debido a que estas ecuaciones pueden simularse usando técnicas Monte-Carlo
que son numéricamente e�cientes. Aplicaciones recientes incluyen el análisis de los efectos de la
termalización en los plasmas quark-gluón [36, 52] y en la colisión de haces de plasmas ultrare-
lativistas [13].
Algunos de los trabajos más recientes en la modelación del movimiento Browniano rela-
tivista utilizando el enfoque descrito en el párrafo anterior fueron desarrollados por los físicos
Jörn Dunkel y Peter Hänggi, en una serie de artículos de los últimos cinco años [18, 19, 20,
21, 22, 24, 25, 26]. En particular estos autores consideran una ecuación de Langevin rela-
tivista usando integrales estocásticas de Itô, de Stratonovich y Backward-Itô; lo cual simpli�ca
el tratamiento del modelo al tener libertad de usar una de las integrales según el problema
especí�co bajo estudio. Asimismo, estos autores propusieron otro dos modelos especí�cos del
movimiento Browniano relativista: el RBM y el RBM(I). Estos dos ejemplos de procesos es-
12
tocásticos relativistas, al igual que el proceso de Ornstein-Uhlenbeck relativista, logran exhibir
el comportamiento límite no-relativista adecuado. En los artículos de Dunkel y Hänggi también
se encuentra las reparametrizaciones de la ecuación de Langevin en términos del tiempo propio
y del tiempo de un observador en movimiento; además, se propone un modelo microscópico
heurístico para determinar ecuaciones de Langevin �apegadas a la realidad�.
El objetivo principal de esta tesis de licenciatura en matemáticas es exponer de manera
rigurosa los modelos de cálculo estocástico para el movimiento Browniano relativista estudia-
dos por Dunkel y Hänggi. En particular, se presentan de manera detallada las condiciones y
demostraciones precisas de los resultados para los procesos estocásticos relativistas de Langevin
expuestos en estos trabajos. La exposición de esta tesis pretende ser autocontenida y teniendo
en mente como lector a personas con interés en la matemática, y en especial en la probabilidad
y los procesos estocásticos.
La estructura del trabajo de tesis es la siguiente. En el Capítulo 2 se introduce la teoría
matemática necesaria para nuestro enfoque al tratamiento de los modelos del movimiento
Browniano relativista: el cálculo estocástico. Se presentan las integrales estocásticas de Itô,
Stratonovich y Backward con sus correspondientes fórmulas de cambio de variable. Asimismo
se presentan resultados clásicos sobre existencia y unicidad para ecuaciones diferenciales es-
tocásticas de Itô, las relaciones entre la integral de Itô y las otras dos integrales y los teoremas
de cambio de tiempo aleatorio. Este capítulo puede considerarse el punto de partida al análisis
del caso relativista, además de que ayuda a clari�car algunos de los objetivos del caso general.
Se recomienda leer el Apéndice de notación para entender alguna de la simbología que se utiliza
en el capítulo.
En el Capítulo 3 se analiza un modelo para el movimiento Browniano no-relativista uni-
dimensional: se establecen las ecuaciones de Langevin y un teorema de �uctuación-disipación,
se deducen las ecuaciones de Fokker-Planck no-relativistas, se analiza el modelo del proceso de
Ornstein-Uhlenbeck y se describe un modelo microscópico para encontrar coe�cientes de ruido y
fricción realistas. La teoría de Langevin para movimiento Brownianos relativistas se estudia en
el Capítulo 4; en donde se establece la ecuación de Langevin relativista, se encuentran las ecua-
ciones de Fokker-Planck correspondientes y se presenta un teorema de �uctuación-disipación
relativista. Para leer el Capítulo 4 y los siguientes, es altamente recomendable leer primero los
13
apéndices B, C, D y E, ésto con el �n de tener en cuenta las convenciones y la notación que se
toma en éste y los próximos apartados. En el Capítulo 5 se proporciona una forma de calcular
constantes de difusión asintóticas para algunos modelos de movimiento Browniano y se analizan
ejemplos especí�cos de procesos estocásticos relativistas. El Capítulo 6 trata las reparametriza-
ciones de Langevin en términos del tiempo propio y del tiempo de un observador en movimiento.
La contraparte relativista del modelo microscópico no-relativista expuesto en el Capítulo 3 se
trata en el Capítulo 7. Los apéndices contienen información y resultados que son útiles desde
el punto de vista matemático y, sobre todo, físico. Aquellas personas sin conocimiento de las
teorías físicas involucradas, se les recomienda comenzar la lectura de esta tesis con los Apéndices
B, C, D y E.
14
Capítulo 2
Preliminares de Cálculo Estocástico
El objetivo de este capítulo es presentar la teoría matemática necesaria para trabajar en la
modelación del movimiento Browniano relativista. Se presentan las ideas y resultados princi-
pales del cálculo estocástico indicando las referencias en donde se pueden ver demostraciones y
aspectos más generales.
2.1. Procesos de Wiener
En esta sección se presenta un breve resumen sobre procesos de Wiener, lo su�ciente para
tratar la integral estocástica con respecto a este proceso. Estos procesos constituyen el ejemplo
más simple y mejor estudiado de un proceso de difusión.
De�nición 1 Un proceso estocástico real fWt; t � 0g se llama un proceso de Wiener si
(i) W0 = 0 a.s.;
(ii) Wt �Ws se distribuye como una variable aleatoria normal N (0; t� s), para cualesquiera
t > s � 0;
(iii) para todos 0 < t1 < t2 < ::: < tn, las variables aleatorias Wt1, Wt2 �Wt1,..., Wtn �Wtn�1
son independientes;
(iv) Wt tiene trayectorias continuas a.s.
15
A los procesos de Wiener también se les denomina movimiento Browniano en el campo de
la probabilidad. Para evitar confusiones con el movimiento Browniano físico, el cual tratare-
mos ampliamente en este trabajo, solo nos referiremos a ellos como procesos de Wiener. Cabe
mencionar que es común identi�car al espacio de parámetros con el tiempo y referirse a la
condiciones (ii) y (iii) como incrementos normales independientes.
Los primeros momentos estadísticos del proceso de Wiener son fáciles de calcular.
Lema 2 Si fWt; t � 0g es un proceso de Wiener, entonces hWti = 0 y hWtWsi = m��n(s; t)
para t � 0, s � 0.
Ahora extenderemos la de�nición de proceso de Wiener al caso n-dimensional.
De�nición 3 Un proceso estocástico fWt = (W 1t ; :::;W
nt ); t � 0g en Rn es un proceso de
Wiener n-dimensional si, para cada k = 1; :::; n, se cumple que
(i) W k es un proceso de Wiener (1-dimensional); y
(ii) las �-álgebras Wk := �(W k(t); t � 0) son independientes.
Usando los momentos estadísticos del proceso de Wiener unidimensional, es fácil ver lo
siguiente.
Lema 4 Si fWt; t � 0g es un proceso de Wiener n-dimensional y �kl una delta de Kronecker,
entonces
(i)W kt W
ls
�= m��n(s; t)�kl para k; l = 1; :::n, t � 0, s � 0; y
(ii)(W k
t �W ks )(W
lt �W l
s)�= (t� s)�kl para k; l = 1; :::; n y t � s � 0.
La di�cultad de integrar sobre un proceso de Wiener es una consecuencia del próximo
resultado.
Teorema 5 Sea fWt; t � 0g es un proceso de Wiener n-dimensional. Entonces, para casi toda
!, la trayectoria muestral t 7�! Wt (!) es de variación in�nita en cada subintervalo y no es
diferenciable en ningún t � 0.
16
2.2. Cálculo estocástico de Itô
En esta sección se presentan las herramientas necesarias del cálculo estocástico de Itô. En
la primera parte se establecen los conceptos de familia y proceso no-anticipante. En la segunda
subsección se de�ne la integral de Itô y se establecen algunas de la propiedades más importantes
de esta. La tercera subsección trata las diferenciales estocásticas y presenta la fórmula de Itô. En
la cuarta parte se de�ne el concepto de ecuación diferencial estocástica (de Itô) y su solución,
se presenta un Teorema de existencia y unicidad para las soluciones y se establecen algunos
otros resultados relacionados. En la Subsección 5 se presentan las condiciones necesarias para
que la solución de una ecuación diferencial estocástica sea un proceso de Markov o un proceso
de difusión. En la Subsección 6 se establecen un par de teoremas relacionados con el cambio de
tiempo en una integral Itô.
Los resultados de esta sección fueron tomados de las referencias [6, 28, 46, 48, 29, 56]; ahí
mismo puede ser consultada la demostración de los resultados. Para profundizar en los temas
tratados en esta sección, se recomienda consultar, además de las referencias anteriores, los libros
de Protter [51] y Karatzas et al. [42].
2.2.1. Procesos estocásticos no-anticipantes
Los procesos estocásticos no-anticipantes son de gran importancia en la de�nición de la
integral estocástica de Itô. Bajo este concepto quedan asentadas las condiciones de medibilidad
que debe cumplir un proceso para poder de�nir su integral estocástica en el sentido de Itô. Antes
de presentar la de�nición de proceso no-anticipante, es necesario de�nir lo que se entiende por
familia de �-álgebras no-anticipante.
De�nición 6 Sea t0 � 0 un real y Wt un proceso de Wiener m-dimensional de�nido sobre el
espacio de probabilidad (;U ;P). Las siguientes �-álgebras resultan de interés:
(i) La �-álgebra W [t0; s;W] := �(Wt; t0 � t � s) es llamada la historia del proceso de
Wiener desde el tiempo t0 hasta (e incluyendo) el tiempo s.
(ii) La �-álgebra W+ [t;W] := �(Ws �Wt; t � s < 1) es el futuro del proceso de Wiener
desde el tiempo t.
17
Cuando no haya confusión sobre el proceso de Wiener al que se hace referencia, se utilizará
W [t0; s] y W+t , en lugar de W [t0; s;W] y W+ [t;W], respectivamente.
De�nición 7 Sea t0 � 0 un real y Wt un proceso de Wiener m-dimensional de�nido sobre el
espacio de probabilidad (;U ;P). Una familia (Ft)t0�t<1 de sub-�-álgebras de U es llamada
no-anticipante respecto a Wt si tiene las siguientes tres propiedades:
(i) Ft � Fs para cualesquiera t � s � t0;
(ii) Ft �W [t0; t] para todo t � t0;
(iii) Ft es independiente de W+t para todo t � t0.
Es usual solo hablar de familia no-anticipante, sin hacer referencia explícita al proceso de
Wiener. En ocasiones, se usará esta convención.
Ahora es posible escribir la de�nición de proceso no-anticipante, sólo hay que tener en cuenta
que Matd�m (R) denota al conjunto de las matrices con coe�cientes reales y dimensión d�m.
De�nición 8 Sea t0 � 0 un real y Wt un proceso de Wiener m-dimensional de�nido sobre
el espacio de probabilidad (;U ;P). Sea (Ft)t0�t<1 una familia de �-álgebras no-anticipante
respecto aWt. Sea G = G (s; !) proceso estocástico de�nido en [t0; t]�, que toma valores en
Matd�m (R) y medible en (s; !) : Se dice que G es no-anticipante con respecto a Ft si G (s; �)
es Fs-medible para todo s 2 [t0; t].
Cuando es claro sobre que familia se trabaja, se acostumbra hablar de procesos no-anticipantes
sin mencionar sobre que familia.
La siguiente proposición establece un propiedad muy útil de los procesos estocásticos no-
anticipantes.
Proposición 9 Sea t0 � 0 un real y (Ft)t0�t<1 una familia de �-álgebras no-anticipante. Sea
g : [t0; t]� Rd ! R una función medible y G un proceso estocástico no-anticipante de�nido en
[t0; t]� y que toma valores en Matd�m (R). Entonces g (t;Gt) es un proceso no-anticipante.
Para poder de�nir la integral de Itô de un proceso estocástico, no es su�ciente que este sea
no-anticipante. También es necesario que el proceso cumpla ciertas condiciones de integrabilidad
(no estocástica). En la de�nición del siguiente conjunto quedan establecidas tales condiciones.
18
De�nición 10 Sea t0 � 0 un real, (Ft)t0�t<1 una familia de �-álgebras no-anticipante y el
espacio muestral donde está de�nido el proceso de Wiener sobre el cual la familia anterior es no-
anticipante. Denotaremos porMd;m2 [t0; t;F ] al conjunto de procesos estocásticos no-anticipantes
G de�nidos en [t0; t] � t > t0 y que toman valores en Matd�m (R), los cuales cumplen con
probabilidad 1 Z t
t0
jG (s; !)j2 ds <1;
siendo jGj :=�Pd
i=1
Pmj=1
�Gij�2�1=2
=�trGGT
�1=2y donde se toma la integral de Lebesgue.
Conviene observar que G 2Md;m2 [t0; t;F ] si, y solo si, Gij 2M1;1
2 [t0; t;F ] para cualesquiera
i y j. Además, se usaráMd;m2 [t0; t], en lugar deM
d;m2 [t0; t;F ], cuando sea claro sobre cual sobre
cual familia de �-álgebras no-anticipante (o sobre cual proceso de Wiener) se trabaja.
Proposición 11 Sea t0 � 0 un real y (Ft)t0�t<1 una familia de �-álgebras no-anticipante.
Entonces, Md;m2 [t0; t] es un espacio vectorial.
2.2.2. De�nición de la integral de Itô
En esta subsección se de�ne la integral de Itô. Para lograr este objetivo, se comienza de�nien-
do la integral de Itô para procesos escalonados.
De�nición 12 Sea t0 � 0 un real y (Ft)t0�t<1 una familia de �-álgebras no-anticipante. Un
proceso G 2Md;m2 [t0; t] se llama proceso escalonado si existe una partición
P := ft0 = t0 < t1 < ::: < tm = tg tal que
Gs � Gtk para 8s 2 [tk; tk+1] , donde k = 0; 1; :::;m� 1.
De�nición 13 Sea t0 � 0 un real, Wt un proceso de Wiener m-dimensional y (Ft)t0�t<1una familia de �-álgebras no-anticipante. Sea G 2Md;m
2 [t0; t] un proceso estocástico como en la
de�nición anterior. La variable aleatoria
Z t
t0
Gs � dWs :=
m�1Xk=0
Gtk (W (tk+1)�W (tk))
es la integral de Itô de G sobre el intervalo [t0; t] respecto al proceso Wt.
19
Para de�nir la integral de Itô de procesos más generales, se aproximan dichos procesos por
funciones escalonadas. Los siguientes dos resultados tratan sobre esta aproximación.
Lema 14 Sea t0 � 0 un real y (Ft)t0�t<1 una familia de �-álgebras no-anticipante. Para cada
proceso G 2Md;m2 [t0; t], existe una sucesión fGng�Md;m
2 [t0; t] de procesos escalonados tal que
ac - l��mn!1
Z t
t0
jG (s)�Gn (s)j2 ds = 0:
Lema 15 Sea t0 � 0 un real y (Ft)t0�t<1 una familia de �-álgebras no-anticipante. Además,
sean G 2Md;m2 [t0; t] y fGng�Md;m
2 [t0; t] una sucesión de procesos escalonados tal que
st - l��mn!1
Z t
t0
jG (s)�Gn (s)j2 ds = 0:
Entonces, se cumple que
st - l��mn!1
Z t
t0
Gn (s) � dWs = I (G) ;
donde I (G) es una variable aleatoria que no depende de la elección de la sucesión fGng.
Del Lema anterior se puede inferir la de�nición de la integral de Itô para procesos más
generales que los escalonados.
De�nición 16 Sea t0 � 0 un real, Wt un proceso de Wiener m-dimensional y (Ft)t0�t<1una familia de �-álgebras no-anticipante respecto a Wt: Para todo proceso G 2Md;m
2 [t0; t], la
integral de Itô de G sobre el intervalo [t0; t] con respecto al proceso de Wiener Wt se de�ne
como Z t
t0
Gs � dWs = st - l��mn!1
Z t
t0
Gn (s) � dWs;
donde fGng � Md;m2 [t0; t] es una sucesión de procesos escalonados que aproximan a G en el
sentido
st - l��mn!1
Z t
t0
jG (s)�Gn (s)j2 ds = 0:
A continuación se enlistan algunas de las propiedades más importantes que cumple la integral
de Itô. La propiedad (v) es de particular importancia pues establece una condición para la
existencia de los primeros momentos estadísticos de la integral de Itô.
20
Teorema 17 Sea t0 � 0 un real, Wt un proceso de Wiener m-dimensional y (Ft)t0�t<1una familia de �-álgebras no-anticipante respecto a Wt: Además, tomemos G 2Md;m
2 [t0; t] y
fGng �Md;m2 [t0; t] sucesión. Entonces, la integral de Itô tiene las siguientes propiedades:
(i)R tt0(aG1 (s) + bG2 (s)) � dWs = a
R tt0G1 (s) � dWs+ b
R tt0G2 (s) � dWs; para cualesquiera
a; b 2 R.
(ii)R tt0Gs � dWs =
Pmk=1
0BBB@R tt0G1ks � dW k
s
...R tt0Gdks � dW k
s
1CCCA :
(iii) Para N > 0 y c > 0,
P�����Z t
t0
Gs � dWs
���� > c
�� N
c2+ P
�Z t
t0
jGsj2 ds > N
�:
(iv) La relación
st - l��mn!1
Z t
t0
jG (s)�Gn (s)j2 ds = 0
implica
st - l��mn!1
Z t
t0
Gn (s) � dWs =
Z t
t0
Gs � dWs:
(v) Si se cumple Z t
t0
DjGsj2
Eds <1;
entonces, tenemos que �Z t
t0
Gs � dWs
�= 0
y *�Z t
t0
Gs � dWs
��Z t
t0
Gs � dWs
�T+=
Z t
t0
DGs (Gs)
TEds:
La integral de Itô tiene una de�nición alternativa para procesos estocásticos continuos a.s.;
tal de�nición involucra sumas de Riemann, lo cual coincide con la de�nición esperada de integral.
Corolario 18 Sea t0 � 0 un real, Wt un proceso de Wiener m-dimensional y (Ft)t0�t<1una familia de �-álgebras no-anticipante respecto a Wt. Sean t > t0 y G 2Md;m
2 [t0; t] un
21
proceso continuo con probabilidad 1. Supongamos que Pn := ft0 = tn0 < tn1 < ::: < tnm = tg es
una sucesión de particiones de [t0; t] con �n = m�ax1�k�n
(tk � tk�1)! 0 cuando n!1. Entonces,
se cumple que Z t
t0
Gs � dWs := st - l��mn!1
nXk=1
Gtk�1
�Wtk�Wtk�1
�:
La integral de Itô esta relacionada con el concepto de martingala, el cual tiene la siguiente
de�nición.
De�nición 19 Sea t0 � 0 un real, (;U ;P) un espacio de probabilidad, (Ut)t0�t�T una familia
de creciente sub-�-álgebras de U y fXt; t0 � t � Tg un proceso estocástico de�nido en (;U ;P)
y que toma valores en Rd. Además, supongamos que Xt es Ut-medible e integrable para cualquier
t 2 [t0; T ]. El par (Xt;Ut) es una martingala si
hXt j Usi = Xs a.s.
para cualesquiera t0 � s � t � T .
Ahora es posible establecer el siguiente resultado, el cual nos proporciona las propiedades
de continuidad y medibilidad de la integral de Itô, así como las condiciones bajo las cuales es
una martingala.
Teorema 20 Sea t0 � 0 un real, Wt un proceso de Wiener m-dimensional y (Ft)t0�t<1 una
familia de �-álgebras no-anticipante respecto a Wt. Supongamos que G 2Md;m2 [t0; T ] y que
Xt =R tt0Gs � dWs para t0 � t � T . Entonces, se cumple lo siguiente:
(i) Xt es Ft-medible.
(ii) Si Z t
t0
DjGsj2
Eds <1 8t � T; (2.1)
entonces (Xt,Ft) ; para t 2 [t0; T ], es una martingala; esto es, para t0 � s � t � T;
hXt j Fsi = Xs:
22
Además, para t; s 2 [t0; T ] tenemos hXti = 0 y
DXt (Xs)
TE=
Z m��n(t;s)
t0
DGu (Gu)
TEdu:
(iii) Xt tiene trayectorias muestrales continuas con probabilidad 1.
Corolario 21 Sea t0 � 0 un real,Wt un proceso de Wiener m-dimensional y (Ft)t0�t<1 una
familia de �-álgebras no-anticipante respecto aWt. Supongamos que G;H 2Md;m2 [t0; T ] y queZ t
t0
DjGsj2
Eds <1,
Z t
t0
DjHsj2
Eds <1:
Entonces, para t; s 2 [t0; T ] se cumple*�Z t
t0
Gu � dWu
��Z s
t0
Hu � dWu
�T+=
Z m��n(t;s)
t0
DGu (Hu)
TEdu:
Observación 22 A pesar de que no se cumpla la condición (2.1), el proceso Xt =R tt0Gs�dWs
es una martingala local continua. La propiedad de martingala local es más general que la de
martingala y tiene la siguiente de�nición.
De�nición 23 Sean (;U ;P) un espacio de probabilidad, (Ft)0�t<1 una �ltración de U y
fXt; t � 0g un proceso estocástico F-adaptado en Rd. El proceso Xt es una martingala lo-
cal si existe una sucesión fTng de tiempos de paro crecientes con l��mTn = 1 a.s. y tal que
Xm��nft;Tng1fTn>0g es una F-martingala para cada n.
En la de�nición anterior, una �ltración es una familia creciente de sub-�-álgebras de U ,
un proceso estocástico fXt; t � 0g es F-adaptado si Xt es Ft-medible para todo t � 0 y una
variable aleatoria T : ! [0;1] es un tiempo de paro si el eventofT � tg 2 Ft para cada
t 2 [0;1).
2.2.3. Las diferenciales estocásticas y la fórmula de Itô
Antes de pasar al estudio del las ecuaciones diferenciales estocásticas, es necesario de�nir el
concepto de diferencial estocástica.
23
De�nición 24 Sea t0 � 0 un real,Wt un proceso de Wiener m-dimensional y (Ft)t0�t<1 una
familia de �-álgebras no-anticipante respecto aWt. Supongamos que fXt; t � t0g es un proceso
estocástico que satisface
Xt = Xt0 +
Z t
t0
fsds+
Z t
t0
Gs � dWs 8t 2 [t0; T ] ; a.s. (2.2)
junto con los siguientes supuestos:
(i) G 2Md;m2 [t0; T ] :
(ii) Xt0 es una variable aleatoria Ft0-medible. En particular, este es el caso cuando Xt0 es
determinística.
(iii) El proceso f toma valores en Rd, es medible en (s; t), es no-anticipante (es decir, f (t; �)
es Ft-medible 8t 2 [t0; T ]) y cumple con probabilidad 1Z T
t0
jf (s; !)j ds <1;
tomando la integral de Lebesgue.
Entonces decimos que el proceso X tiene la diferencial estocástica
dXt = ftdt+Gt � dWt: (2.3)
Observación 25 Las dos integrales en la ecuación (2.2) son procesos estocásticos continuos,
así que el proceso Xt toma valores en Rd y, con probabilidad 1, tiene trayectorias muestrales
continuas. Además, Xt es Ft-medible (por tanto, no-anticipante) y se tiene que
Xt = Xs +
Z t
sfudu+
Z t
sGu � dWu;
para todo s tal que t0 � s � t � T . Debido a esta última propiedad, podemos habla de diferen-
ciables que se cumplen solo en (t0; T ) sin necesidad de establecer la condición inicial Xt0.
Básicamente, las diferenciales estocásticas son una notación más compacta y simple para
expresar relaciones de la forma (2.2). Conviene notar que, desglosando expresiones, la diferencial
24
estocástica (2.3) se traduce en el sistema
dX it = f itdt+
mXj=1
Gijt dWjt para i = 1; ::; n,
siendo X =�X1; : : : ; Xd
�,W =
�W 1; : : : ;W d
�, f =
�f1; : : : ; fd
�y G =
�Gij�.
Para facilitar el tratamiento de las diferenciales estocásticas, cuando se presente una dife-
rencial de la forma (2.3), se asumirá implícitamente que W es un proceso de Wiener (de
alguna dimensión a especi�car) y que f y G cumplen las propiedades de la de�nición anterior
para alguna familia de �-álgebras no-anticipante (respecto a W). La dimensión d tiene que
especi�carse y el valor inicial t0 coincidirá con el mínimo del intervalo de de�nición (el cual
tiene que ser especi�cado) del parámetro temporal del proceso G. Supuestos similares aplicarán
cuando solo se tenga a la integralR tt0Gu � dWu.
De ahora en adelante, se toman 0 � t0 � T . A continuación se presenta un resultado
fundamental en nuestro tratamiento, el cual es análogo a la regla de la cadena.
Teorema 26 (Fórmula de Itô) Denotemos por u = u (t;x) a una función continua de�nida
sobre [t0; T ]� Rd que toma valores en Rk y que tiene derivadas parciales continuas
@
@tu (t;x) = ut (t;x) ;
@
@xiu (t;x) = uxi (t;x) ;
@2
@xi@xju (t;x) = uxixj (t;x) ;
para i; j = 1; : : : ; d. Sea fXt; t0 � t � Tg un proceso estocástico en Rd, el cual esta de�nido por
la diferencial estocástica
dXt = ftdt+Gt � dWt;
donde W es un proceso de Wiener m-dimensional. Entonces, el proceso en Rk dado por
Yt := u (t;Xt) 8t 2 [t0; T ] ;
25
tiene diferencial estocástica con respecto al mismo proceso de WienerWt y esta tiene la forma
dYt =
0@(ut (t;Xt))T + ux (t;Xt) ft + 12
dXi;j=1
(uxixj (t;Xt))T�Gt (Gt)
T�ij1A dt
+ux (t;Xt)Gt � dWt;
donde (uxixj (t;Xt))T y (ut (t;Xt))
T son vectores columna, ux = (ux1 ; :::; uxd) es una matriz de
dimensión k � d y�Gt (Gt)
T�ij
denota la entrada ij de la matriz�Gt (Gt)
T�.
La fórmula de Itô tiene expresiones mucho más simples para ciertos valores de m, k y d.
Estos corolarios pueden ser consultados en las referencias propuestas.
2.2.4. Ecuaciones diferenciales estocásticas
En esta sección, se exponen algunos resultados de la teoría de ecuaciones diferenciales es-
tocásticas de Itô. Comenzaremos con la de�nición de solución de una ecuación diferencial es-
tocástica, posteriormente trataremos un teorema de existencia y unicidad de las soluciones; y
culminaremos con un par de ejemplos.
Consideremos la diferencial estocástica
dXt = f (t;Xt) dt+G (t;Xt) � dWt; Xt0 = c; (2.4)
válida para t0 � t � T <1; o bien, su forma integral
Xt = c+
Z t
t0
f (s;Xs) ds+
Z t
t0
G (s;Xs) � dWs; (2.5)
la cual se cumple también para t0 � t � T < 1. En las dos ecuaciones anteriores, Xt es un
proceso estocástico (conocido) de�nido en [t0; T ] yWt es un proceso de Wiener m-dimensional.
Además, se toman las funciones f : [t0; T ] � Rd ! Rd y G : [t0; T ] � Rd ! Matd�m (R) ; las
cuales se suponen medibles en [t0; T ]�Rd y no son procesos estocásticos (i.e. son independientes
de ! 2 ).
En las ecuaciones (2.4) y (2.5), el proceso Xt se supone conocido y que cumple todas las
propiedades necesarias para que el lado derecho de la diferencial estocástica (2.4) este perfec-
26
tamente de�nido. En particular, Xt debe ser no-anticipante. Las ecuaciones de la forma (2.4)
también se pueden interpretar como la de�nición de un proceso estocástico desconocido Xt con
valor inicial Xt0 = c.
Con respecto a la familia Ft de �-álgebras inherente, se toma la siguiente convención.
Convenci�on : Con el propósito de tratar ecuaciones diferenciales estocásticas (de Itô)
sobre el intervalo [t0; T ], siempre es su�ciente elegir a Ft como la �-álgebra más pequeña con
respecto a la cual el valor inicial c y las variables aleatorias Ws s � t son medibles. En otras
palabras
Ft = � (c;Ws (s � t)) 8t � t0:
Cuando no se especi�que una familia (Ft)t0�t<1 para una diferencial estocástica como (2.4),
se supondrá que se toma de esta manera.
Ahora es posible de�nir lo que es una ecuación diferencial estocástica de Itô y la solución
de ella.
De�nición 27 Una ecuación de la forma (2.4) se denomina ecuación diferencial estocástica de
Itô. La variable aleatoria c se denomina el valor inicial al tiempo t0. Un proceso estocástico Xt
es una solución de la ecuación (2.4) sobre el intervalo [t0; T ] si tiene las siguientes propiedades:
(i) Xt es Ft-medible, es decir no-anticipante para t 2 [t0; T ].
(ii) Los procesos f (s; !) = f (s;Xs) yG (s;Xs (!)) = G (s;Xs (!)) satisfacen con probabilidad
1 Z T
t0
��f (s; !)�� ds <1, Z T
t0
��G (s; !)��2 ds <1(es decir, G 2Md;m
2 [t0; T ]). Así, el lado derecho de la diferencial estocástica (2.4) tiene
sentido.
(iii) La ecuación integral (2.5) se cumple para todo t 2 [t0; T ] con probabilidad 1.
Observación 28 Debido a que Xt es no-anticipante, los procesos f (s; !) y G (s;Xs (!)) son
no-anticipantes.
27
Conviene observar que la ecuación (2.5) es equivalente a
Xt �Xs = +Z t
sf (u;Xu) du+
Z t
sG (u;Xu) � dWu; Xt0 = c;
donde t0 � s � t � T . Además, siempre se considerarán soluciones continuas de las ecuaciones
diferenciales estocásticas (de acuerdo a la observación 25 es posible hacerlo) cuando estas tengan
solución.
A continuación presentamos un teorema de existencia y unicidad para las soluciones de una
ecuación diferencial estocástica.
Teorema 29 (Existencia y unicidad) Supongamos que tenemos una ecuación diferencial
estocástica
dXt = f (t;Xt) dt+G (t;Xt) � dWt; Xt0 = c; t0 � t � T <1; (2.6)
dondeWt es un proceso de Wiener m-dimensional y c es una variable aleatoria independiente de
Wt�Wt0 para t � t0: Supongamos que las funciones f : [t0; T ]�Rd ! Rd y G : [t0; T ]�Rd !
Matd�m (R) son medibles en [t0; T ] � Rd y que tienen las siguientes propiedades: Existe una
constante K > 0 tal que
(i) (Condición de Lipschitz) para cualesquiera t 2 [t0; T ], x 2 Rd, y 2 Rd,
jf (t;x)� f (t;y)j+ jG (t;x)�G (t;y)j � K jx� yj :
(ii) (Restricción de crecimiento) Para todo t 2 [t0; T ] y x 2 Rd,
jf (t;x)j2 + jG (t;x)j2 � K2�1 + jxj2
�:
Entonces, la ecuación (2.4) tiene una única solución Xt sobre [t0; T ] que satisface la condi-
ción inicial Xt0 = c y que es continua con probabilidad 1; es decir, si Xt y Yt son soluciones
28
continuas de (2.4) con la misma condición inicial c, entonces
P
"sup
t0�t�T:jXt �Ytj > 0
#= 0:
El teorema anterior establece condiciones que pueden ser difíciles o tediosas de comprobar.
Los siguientes dos corolarios establecen condiciones bajo las cuales el resultado del Teorema 29
sigue siendo cierto, pero que son más fáciles de comprobar.
Corolario 30 El Teorema 29 de existencia y unicidad permanece válido si reemplazamos la
condición de Lipschitz con la condición más general: 8N > 0, existe una constante KN tal que,
para cualesquiera t 2 [t0; T ], jyj � N y jxj � N ,
jf (t;x)� f (t;y)j+ jG (t;x)�G (t;y)j � KN jx� yj : (2.7)
Corolario 31 Para que la condición de Lipschitz en el Teorema 29 de existencia y unicidad
(o su generalización (2.7))se cumpla, es su�ciente que las funciones f (t;x) y G (t;x) tengan
derivadas parciales @@xif (t;x), @
@xiG (t;x) ; i = 1; :::; d continuas para todo t 2 [t0; T ] y que estén
acotadas en [t0; T ]�Rd (o, en el caso de la generalización, acotadas sobre [t0; T ]�fjxj � Ng).
En el tratamiento del movimiento Browniano físico, resultan de gran importancia las ecua-
ciones autónomas. Para este tipo de ecuaciones puede establecerse un resultado de existencia y
unicidad más simple.
De�nición 32 Consideremos la ecuación diferencial estocástica (2.6). Diremos que es esta
ecuación es autónoma si las funciones f (t;x) y G (t;x) no dependen del tiempo; es decir, si
f (t;x) � f (x) y G (t;x) � G (x).
Corolario 33 Consideremos la ecuación diferencial estocástica autónoma
dXt = f (Xt) dt+G (Xt) � dWt; Xt0 = c;
donde Wt es un proceso de Wiener m-dimensional, c una variable aleatoria independiente de
Wt �Wt0 para t � t0 y se toman funciones f : Rd ! Rd y G : Rd ! Matd�m (R). Si existe
29
una constante K > 0 tal que las funciones f y G cumplen la condición de Lipschitz global
jf (x)� f (y)j+ jG (x)�G (y)j � K jx� yj 8 (x;y) 2 R2d; (2.8)
entonces la ecuación diferencial estocástica tiene exactamente una solución global Xt sobre el
intervalo [t0;1) tal que Xt0 = c.
2.2.5. Ecuaciones diferenciales estocásticas y los procesos de difusión
En esta sección se presentan las condiciones bajo las cuales la solución de una ecuación
diferencial estocástica es un proceso difusión. El material sobre procesos de difusión que se
necesita para leer este trabajo se presenta en el Apéndice A de procesos de Markov.
Considérese la ecuación diferencial estocástica
dXt = f (t;Xt) dt+G (t;Xt) � dWt; Xt0 = c; (2.9)
sobre el intervalo [t0; T ]. En esta ecuación, Xt es un proceso estocástico en Rd, f toma valores
en Rd, G toma valores en Matd�d (R) y Wt es un proceso de Wiener m-dimensional. El valor
inicial c es una variable aleatoria independiente deWt �Wt0 para t � t0.
También considérese la restricción de la ecuación (2.9) al intervalo [s; T ] � [t0; T ], la cual
puede expresarse como
Xt = x+
Z t
sf (u;Xu) du+
Z t
sG (u;Xu) � dWu; t0 � s � t � T; (2.10)
siendo Xs = x 2Rd.
El siguiente Teorema permite relacionar la solución de una ecuación diferencial estocástica
con los procesos de Markov.
Teorema 34 Si la ecuación (2.9) satisface las condiciones del Teorema de existencia y unicidad
29, la solución Xt de la ecuación (2.9) es un proceso de Markov sobre el intervalo [t0; T ] cuya
distribución de probabilidad inicial al instante t0 es la distribución de c y cuyas probabilidades
30
de transición están dadas por
P (s;x; t; B) = P (Xt 2 B j Xs = x) = P [Xt (s;x) 2 B] ;
donde Xt (s;x) es la solución de la ecuación (2.10).
Dado el Teorema anterior, resulta conveniente tener condiciones para las cuales la solución
de una ecuación diferencial estocástica es un proceso de Markov homogéneo.
Teorema 35 Supongamos que las condiciones del Teorema de existencia y unicidad 29 se
satisfacen por la ecuación (2.9). También supongamos que los coe�cientes f (t;x) � f (x) y
G (t;x) � G (x) son independientes del tiempo t sobre el intervalo [t0; T ]. Entonces, la solución
Xt de la ecuación (2.9) es un proceso de Markov homogéneo con probabilidades de transición
(estacionarias)
P (Xt 2 B j Xt0 = x) = P (t� t0;x;B) = P [Xt (t0;x) 2 B] ;
donde Xt (t0;x) es la solución de la ecuación (2.9) con valor inicial x =c = Xt0. En particular,
la solución de una ecuación autónoma
dXt = f (Xt) dt+G (Xt) � dWt; t � t0;
es un proceso de Markov homogéneo de�nido para todo t � t0.
Observación 36 Cuándo es la solución Xt de la ecuación (2.9) un proceso de Markov esta-
cionario? De acuerdo la observación (105), Xt tiene que ser homogéneo, lo cual sucede en el
caso de una ecuación autónoma. Para la existencia de una distribución estacionaria P , es decir
una distribución con la propiedad
P (B) =
ZRdP (t;x; B) dP (x) ; B 2 B
�Rd�; t � 0;
existen algunas condiciones analíticas [véase, [50] págs 272, 273]. En muchos casos, la densidad
p (x) de la distribución estacionaria puede obtenerse de la forma estacionaria de la ecuación
31
(A.1)dXi=1
@
@xi�f i (x) p (x)
�� 12
dXi;j=1
@2
@yi@yj
��G (x) [G (x)]T
�ijp (x)
�= 0:
El siguiente teorema proporciona condiciones para que la solución de una ecuación diferencial
estocástica sea un proceso de difusión.
Teorema 37 Supongamos que las condiciones del Teorema 29 de existencia y unicidad se
satisfacen para la ecuación diferencial estocástica
dXt = f (t;Xt) dt+G (t;Xt) � dWt; Xt0 = c; t0 � t � T;
donde Xt y f (t;x) toman valores en Rd, G (t;x) toma valores en Matd�m (R) y Wt es un
proceso de Wiener m-dimensional. Además, supongamos que las funciones f y G son continuas
con respecto a t. Entonces, la solución Xt es un proceso de difusión d-dimensional sobre [t0; T ]
con vector de deriva f (t;x) y matriz de difusión B (t;x) = [G (t;x)] [G (t;x)]T; y las relaciones
límite de la de�nición 106 se cumplen uniformemente en s 2 [t0; T ). En particular, la solución
de una ecuación diferencial estocástica autónoma siempre es un proceso de difusión homogéneo
sobre [t0,1).
2.2.6. Cambio de tiempo aleatorio
En esta sección se establecen un par de resultados relacionados con el cambio del parámetro
temporal en integrales estocásticas de Itô. El primero de estos resultados da algunas condiciones
necesarias para que un cambio de tiempo sea válido.
Teorema 38 Sea fWt; t � 0g un proceso de Wiener unidimensional y fGt; t � 0g un pro-
ceso estocástico real no-anticipante tal que PhR t0 G
2sds <18t � 0
i= 1: Además, pongamos
t (t) :=R t0 G
2sds 8t � 0 y tomemos a t�1 como la inversa por la derecha de t, es decir,
t�1 (t) = m��n (s : t (s) = t) 8t < t (1). Entonces, cWt :=R t�1(t)0 Gs � dWs es un proceso de
Wiener unidimensional para todo t < t (1) :
Observación 39 Conviene tener en cuenta lo siguiente:
(i) La función t�1 es continua por la izquierda;
32
(ii) Si tomamos x (t) :=R t0 Gsds, la función t se conoce como el tiempo intrínseco de x.
El siguiente Teorema nos proporciona una relación explícita entre una integral estocástica
y sus cambios de tiempo.
Teorema 40 Sea fWt; t � 0g un proceso de Wiener unidimensional. Sean fGt; t � 0g y fFt; t � 0g
procesos estocásticos reales continuos a.s. y no-anticipantes, tales que 0 < F <1 y
P�Z t
0G2sds+
Z t
0F�2s ds <1;8t � 0
�= 1:
De�nimos t (t) :=R t0 F
�2s ds y x (t) :=
R t0 Gs � dWs. Entonces, se cumple que
x�t�1 (t)
�=
Z t
0Gt�1(s)Ft�1(s) � dcWs 8t < t (1) ;
o en términos de diferenciales
dx�t�1 (t)
�= Gt�1(t)Ft�1(t) � dcWt = Gt�1(t)
�d
dtt�1 (t)
�1=2� dcWt; 8t 2 [0; t (1)) ;
donde cWt es un proceso de Wiener dado por cWt :=R t�1(t)0 F�1s dWs.
Además, se establecen un par de resultados que se utilizan en el Capítulo 6 para demostrar
un cambio de tiempo similar al del Teorema 38 para procesos de Wiener d-dimensionales. Antes
debemos establecer la de�nición de proceso de covariación cuadrática.
De�nición 41 Sean fXt; t � 0g y fYt; t � 0g dos procesos estocástico en R. El proceso de co-
variación cuadrática [X;Y ]t de Xt � X (t) y Yt � Y (t) se de�ne como el límite
[X;Y ]t = l��mn!1
2n�1Xk=0
�X
�(k + 1) t
2n
��X
�kt
2n
���Y
�(k + 1) t
2n
�� Y
�kt
2n
��:
En el caso de un proceso de Wiener d-dimensional W =�W 1; : : : ;W d
�, la covariación
cuadrática�W i;W j
�puede expresarse como
�W i;W j
�t=
8<: t si i = j;
0 si i 6= j:
33
Ahora estamos en posibilidades de establecer dos teoremas muy importantes. El primero de
ellos proporciona una expresión para la covariación cuadrática de las integrales unidimensionales
de Itô y el segundo es una caracterización de los procesos de Wiener.
Teorema 42 SeaWt un proceso de Wiener m-dimensional, c una variable aleatoria indepen-
diente de Wt �W0 para t � 0 y Ft = � (c;Ws (s � t)). Además, sea fHt; t � 0g un proceso
estocástico continuo con probabilidad 1 tal que es H 2 Md;m2 [0; T ;F ] 8T > 0. Entonces, si
Y jt =R t0 Hs � dW
js , se cumple que
�Y i; Y j
�t=
Z t
0H2s d�W i;W j
�s
8t � 0:
Teorema 43 (Lévy) Sea X =�X1; : : : ; Xd
�una martingala local continua tal que
�Xi; Xj
�t=
8<: t, si i = j,
0, si i 6= j.
Entonces X es un movimiento Browniano n-dimensional.
2.3. La integral de Stratonovich
En esta sección se de�ne la integral de Stratonovich-Fisk y se establecen relaciones con la
integral de Itô. Para mayor información se recomienda consultar [51, 53].
De�nición 44 Sea Wt un proceso de Wiener m-dimensional y (Ft)0�t<1 una familia de �-
álgebras no-anticipante respecto a Wt. Sea G 2 Md;m2 [0; T ] un proceso continuo a.s. Supon-
gamos que Pn := f0 = tn0 < tn1 < ::: < tnm = Tg es una sucesión de particiones de [0; T ] con
�n = m�ax1�k�n
(tk � tk�1)! 0 cuando n!1. La integral de Stratonovich del proceso Gt sobre el
proceso de Wiener Wt en el intervalo [0; T ] es
Z T
0Gt � dWt := st - l��m
n!1
nXk=1
�Gtk +Gtk�1
2
��Wtk�Wtk�1
�:
En la de�nición anterior se toma como límite inferior a 0, sin embargo es posible extender
la de�nición de la integral de Stratonovich para límites inferiores no-negativos. Conviene ob-
34
servar que, mientras en la integral de Itô se usa Gtk en las sumas de Riemann, en la integral
de Stratonovich estas sumas utilizan el promedioGtk
+Gtk�12 . Esta �ligera�diferencia produce
procesos estocásticos no equivalentes con propiedades distintas.
De�niremos la diferenciales de estocásticas de Stratonovich de manera análoga al caso de
Itô.
De�nición 45 Sea Wt un proceso de Wiener m-dimensional y (Ft)0�t<1 una familia de �-
álgebras no-anticipante respecto a Wt. Supongamos que fXt; t � 0g es un proceso estocástico
que satisface con probabilidad 1
Xt = X0 +
Z t
0fsds+
Z t
0Gs � dWs
junto con los siguientes supuestos:
(i) G 2Md;m2 [0; T ] :
(ii) X0 es una variable aleatoria F0-medible.
(iii) El proceso f toma valores en Rd, es medible en (s; t), es no-anticipante (es decir, f (t; �)
es Ft-medible 8t 2 [0; T ]) y cumple con probabilidad 1Z T
0jf (s; !)j ds <1;
tomando la integral de Lebesgue.
Entonces decimos que el proceso X tiene la diferencial estocástica
dXt = ftdt+Gt � dWt:
El siguiente resultado es un muy útil pues permite cambiar de la integral de Stratonovich a
la integral de Itô.
Teorema 46 Sea Wt un proceso de Wiener m-dimensional y (Ft)0�t<1 una familia de �-
álgebras no-anticipante respecto aWt. Sean f : [0; T ]�Rd ! Rd y G : [0; T ]�Rd ! Matd�m (R)
35
funciones medibles en [0; T ]�Rd con G de clase C1. Supongamos que fXt; t � 0g es un proceso
estocástico que satisface
dXt = f (t;Xt) dt+G (t;Xt) � dWt; Xt0 = c;
para 0 � t � T <1. Entonces, Xt satisface la diferencial de Itô
dXt = [f (t;Xt) + g (t;x)] dt+G (t;Xt) � dWt; Xt0 = c;
donde
gi (t;x) =1
2
mXk=1
dXj=1
@Gik
@xj(t;x)Gjk (t;x) i = 1; :::; d:
Un resultado fundamental para el tratamiento del movimiento Browniano relativista se es-
tablece a continuación. Establece la regla de cambio de variable de la diferencial de Stratonovich.
Teorema 47 Sea u : [0; T ]�Rd ! Rk una función de clase C2. Además, supongamos que se
cumple
dXt = ftdt+Gt � dWt,
con Wt un proceso de Wiener m-dimensional. Si de�nimos
Yt := u (Xt; t) ;
entonces se cumple que
dYt = (ut (t;Xt)) + ux (t;Xt) ftdt+ ux (t;Xt)Gt � dWt,
donde ux = (ux1 ; :::; uxd) es una matriz de dimensión k � d.
Vale la pena notar que las reglas del cálculo diferencial ordinario se conservan con la integral
de Stratonovich.
36
2.4. La integral Backward-Itô
En esta sección se de�ne la integral Backward-Itô y se establecen relaciones con la integral de
Itô. Para mayor información sobre este tipo de integral y sus respectivas ecuaciones diferenciales,
se recomienda consultar la referencia [43].
De�nición 48 Sea Wt un proceso de Wiener m-dimensional y (Ft)0�t<1 una familia de �-
álgebras no-anticipante respecto a Wt. Sea G 2 Md;m2 [0; T ] un proceso continuo a.s. Supon-
gamos que Pn := f0 = tn0 < tn1 < ::: < tnm = Tg es una sucesión de particiones de [0; T ] con
�n = m�ax1�k�n
(tk � tk�1) ! 0 cuando n ! 1. La integral backward del proceso Gt sobre el
proceso de Wiener Wt en el intervalo [0; T ] es
Z T
0Gt � dWt := st - l��m
n!1
nXk=1
Gtk
�Wtk�Wtk�1
�:
Conviene observar que en la integral backward las sumas de Riemann se evalúan en los
extremos superiores tnk+1 de los intervalos�tnk ;
nk+1
�; mientras que en la integral de Itô las sumas
de Riemann se evalúan en los extremos inferiores tnk . Al igual que en el caso de Stratonovich,
esta diferencia produce procesos estocásticos no equivalentes con propiedades distintas. Además,
esta integral no produce procesos que sean martingala de un proceso de Markov.
La de�nición general de integral Backward puede hacerse a pesar de que el proceso Gt no
sea no-anticipante. Sin embargo, en esta tesis, solo nos importa el caso no-anticipante, ya que
nos interesa su relación con la integral de Itô.
Las diferenciales estocásticas backward se de�nen de manera similar al caso de Itô.
De�nición 49 Sea Wt un proceso de Wiener m-dimensional y (Ft)0�t<1 una familia de �-
álgebras no-anticipante respecto a Wt. Supongamos que fXt; t � 0g es un proceso estocástico
que satisface con probabilidad 1
XT = X0 +
Z T
0fsds+
Z T
0Gs � dWs
junto con los siguientes supuestos:
(i) G 2Md;m2 [0; T ] :
37
(ii) X0 es una variable aleatoria F0-medible.
(iii) El proceso f toma valores en Rd, es medible en (s; t), es no-anticipante (es decir, f (t; �)
es Ft-medible 8t 2 [0; T ]) y cumple con probabilidad 1Z T
0jf (s; !)j ds <1;
tomando la integral de Lebesgue.
Entonces decimos que el proceso X tiene la diferencial estocástica
dXt = ftdt+Gt � dWt:
Bajo las condiciones con las que se de�nió la integral backward, se tiene una relación entre
la integral backward y la integral de Itô.
Teorema 50 Sea Wt un proceso de Wiener m-dimensional y (Ft)0�t<1 una familia de �-
álgebras no-anticipante respecto aWt. Sean f : [0; T ]�Rd ! Rd y G : [0; T ]�Rd ! Matd�m (R)
funciones medibles en [0; T ]�Rd con G de clase C1. Supongamos que fXt; t � 0g es un proceso
estocástico que satisface
dXt = f (t;Xt) dt+G (t;Xt) � dWt; Xt0 = c;
para 0 � t � T <1. Entonces, Xt satisface la diferencial de Itô
dXt = [f (t;Xt) + g (t;x)] dt+G (t;Xt) � dWt; Xt0 = c;
donde
gi (t;x) =mXk=1
dXj=1
@Gik
@xj(t;x)Gjk (t;x) i = 1; :::; d:
La integral backward tiene una fórmula de cambio de variable muy similar a la de la integral
de Itô.
38
Teorema 51 Sea u : [0; T ]�Rd ! Rk una función de clase C2. Además, supongamos que se
cumple
dXt = ftdt+Gt � dWt,
con Wt un proceso de Wiener m-dimensional. Si de�nimos
Yt := u (Xt; t) ;
entonces se cumple que
dYt =�(ut (t;Xt))
T + ux (t;Xt) ft � h (t;Xt)�dt+ ux (t;Xt)Gt � dWt,
donde h (t;Xt) := 12
Pdi;j=1 (uxixj (t;Xt))
T�Gt (Gt)
T�ij, (uxixj (t;Xt))
T y (ut (t;Xt))T son
vectores columna, ux = (ux1 ; :::; uxd) es una matriz de dimensión k � d y�Gt (Gt)
T�ij
de-
nota la entrada ij de la matriz�Gt (Gt)
T�.
Vale la pena notar que la diferencia fundamental con el caso de Itô es el signo menos que
aparece delante del primer símbolo de suma.
39
Capítulo 3
El Movimiento Browniano
No-Relativista
El objetivo de este capítulo es presentar un modelo fenomenológico que describa la dinámica
de una partícula Browniana utilizando herramientas del cálculo estocástico y la física clásica.
Este modelo constituye nuestro punto de partida para tratar el movimiento Browniano rela-
tivista; en particular, las ecuaciones de Langevin y Fokker-Planck que se consideran en este
capítulo constituyen casos límite de la teoría relativista que se estudia a lo largo de la tesis.
En la Sección 1 se introducen las ecuaciones de Langevin y Fokker-Planck, se establece un
teorema de �uctuación-disipación y se analiza el proceso de Ornstein-Uhlenbeck. En la segunda
sección se presenta el modelo de colisión binaria.
3.1. Ecuaciones de Langevin y Fokker-Planck
Con �nes de claridad conceptual, nos restringiremos al caso más simple donde los movimien-
tos están con�nados a una sola dimensión espacial. La generalización a dimensiones mayores es
directa.
3.1.1. La ecuación de Langevin
Considérese el movimiento de una partícula Browniana puntual con masaM > 0, la cual está
inmersa en un baño térmico estacionario y homogéneo que consiste, por ejemplo, de partículas
41
líquidas más pequeñas (masa m � M) a temperatura constante T . Además, supóngase que
la partícula Browniana está en equilibrio térmico con el baño. El marco de referencia inercial
� del baño térmico se identi�ca como el marco de referencia del laboratorio. La posición, la
velocidad y el momento no-relativista de la partícula Browniana en � al tiempo t; se toman
como Xt, Vt := dXtdt y Pt := MVt, respectivamente. Debido al comportamiento aleatorio de la
partícula Browniana, se asume que fXt; t � 0g, fVt; t � 0g y fPt; t � 0g son procesos estocás-
ticos de�nidos en algún espacio de probabilidad.
De acuerdo a la idea de Langevin del movimiento Browniano, la dinámica estocástica de la
partícula Browniana (debida a la interacción con el baño térmico) en presencia de un campo de
fuerza conservativo externo F (t; x), puede ser descrita por la ecuación diferencial estocástica
dXt =PtMdt; (3.1a)
dPt = F (t;Xt) dt� � (Pt)Ptdt+ [2D (Pt)]1=2 � dWt; (3.1b)
complementándola con las condiciones iniciales deterministas X0 = x0 y P0 = p01. A la igualdad
(3.1) se le conoce como la ecuación de Langevin. En el desarrollo de esta tesis se hará especial
énfasis en el caso F � 0.
En la ecuación (3.1b), Wt denota un proceso de Wiener y el símbolo � señala el tipo
de integral estocástica usada. Solo se consideran los tres tipos de integral estocástica que se
mencionan en el Capítulo 2; así, �� = ��cuando se toma la integral de Itô, �� = ��para la
integral de Stratonovich y �� = ��cuando se trabaja con la integral backward. La elección de la
integral se reduce a un mero problema de conveniencia. Además, se toman F : (0;1)�R! R;
� : R! R y D : R! R+ � f0g de clase C1. Por ser F , � y D de clase C1, la ecuación (3.1)
tiene sentido matemático para los tres tipos de integral considerados.
El segundo término del lado derecho de la ecuación (3.1b) representa una fuerza de fricción
con coe�ciente � (p); esta fuerza depende de detalles microscópicos de la interacción baño-
partícula. El último término del lado derecho de la ecuación (3.1b), [2D (Pt)]1=2 � dWt, se
denomina fuerza estocástica de Langevin y en ella quedan re�ejadas las �uctuaciones en el
baño circundante. Al coe�ciente D > 0 de la fuerza de Langevin se le conoce como coe�ciente
1Sin perder generalidad, su puede �jar el tiempo inicial en t0 = 0:
42
de ruido y se interpreta como la amplitud de esta fuerza. Para un baño térmico no homogéneo,
las funciones � y D también dependen de la posición de la partícula; mientras que para un
baño térmico no estacionario, las mismas funciones dependen además del tiempo.
Conviene tener en cuenta los cuatro supuestos físicos implícitos en la ecuación (3.1):
El baño térmico es espacialmente homogéneo y estacionario, es decir que los proceso de
relajación en el baño térmico ocurren en escalas de tiempo mucho menores que las escalas
temporales relevantes asociadas al movimiento de la partícula Browniana pesada.
Las colisiones entre la partícula Browniana y las partículas del baño térmico suceden
virtualmente sin correlación.
A un nivel macroscópico, las interacciones entre la partícula Browniana y el baño térmico
son lo su�cientemente bien descritas por el coe�ciente de fricción � y la fuerza estocástica
de Langevin.
La ecuación de Langevin (3.1) se cumple en el marco de referencia � del laboratorio.
En la Sección 3.2 se estudia la forma en que pueden derivarse o, al menos, motivarse ecua-
ciones diferenciales estocásticas similares a la igualdad (3.1).
3.1.2. La ecuación de Langevin y las integrales estocásticas
En el Capítulo 2 se presentaron las condiciones bajo las cuales las ecuaciones diferenciales
de Itô, Stratonovich y Backward-Itô generan procesos estocásticos equivalentes. Debido a estas
condiciones, la elección de una integral en la ecuación de Langevin se reduce a un problema de
conveniencia. Para ilustrar esto, considérese la ecuación de Langevin con la integral backward-
Itô dada por
dXt =PtMdt; (3.2a)
dPt = [F (t;Xt)� � (Pt)Pt] dt+ [2D (Pt)]1=2 � dWt: (3.2b)
43
Para cada par de coe�cientes (� (p) ; D (p)), es posible encontrar nuevos coe�cientes de fricción
��; �� : R! R de tal manera que las igualdades8<: dXt =PtM dt;
dPt = [F (t;Xt)� �� (Pt)Pt] dt+ [2D (Pt)]1=2 � dWt:
8<: dXt =PtM dt;
dPt = [F (t;Xt)� �� (Pt)Pt] dt+ [2D (Pt)]1=2 � dWt:
describen la misma dinámica que la ecuación (3.2). A saber,
�� (p) � � (p)�D0 (p) =p; �� (p) � � (p)�D0 (p) = (2p) ;
para cualquier p 2 R�f0g.
Cada uno de los tres tipos de integral estocástica considerados tiene sus pros y sus contras,
entre los que podemos mencionar:
La integral de Itô (�) tiene un formalismo matemático excelentemente establecido y bas-
tante conocido; además de tener la calidad de martingala, uno de los conceptos fundamen-
tales de la teoría moderna de procesos estocásticos. Sin embargo, algunas de las cantidades
y funciones útiles en la modelación del movimiento Browniano pueden tener expresiones
más complicadas que con la integral backward.
La integral de Stratonovich tiene la ventaja de cumplir las reglas usuales del cálculo
diferencial, pero no tiene la propiedad de martingala y también algunas de las cantidades
y funciones importantes pueden tener expresiones más complicadas que con la integral
backward.
La integral backward facilita algunos cálculos relacionados en la modelación del movimien-
to Browniano, tiene la desventaja de ser una integral poco común y no ser martingala.
A lo largo de este Capítulo, e incluso posteriores, se trabaja con la integral que más convenga
en el estudio bajo consideración. Básicamente utilizaremos la integral de Itô y la backward; la
integral de Stratonovich solo se presenta en esta tesis debido a que es usual encontrarla en
44
aplicaciones.
3.1.3. La ecuación de Fokker-Planck
En Física, cuando se estudian ecuaciones diferenciales estocásticas como la ecuación de
Langevin (3.1), suele ser también de interés conocer la función de densidad g (t; x; p) asociada
al proceso (Xt; Pt). Esta densidad no siempre existe y algunas condiciones para su existencia
en dimensiones mayores resultan ser un tanto complicadas. Sin embargo, en este trabajo esto
no es tan trascendente pues se construyen procesos de Langevin a partir del hecho de que tal
densidad existe, o que sus densidades marginales existen.
Supóngase que la densidad asociada al proceso de Langevin (3.1) con la integral backward
existe y que la ecuación (3.1) cumple las condiciones para la existencia de la ecuación de
Fokker-Planck. Debido a que para dos integrales estocásticas distintas la integral de Langevin
no genera procesos estocásticos equivalentes, se tiene que la densidad del proceso (3.1) depende
de la integral usada. Denótese por f (t; x; p) la densidad asociada al proceso de Langevin (Xt; Pt)
con la integral backward, la cual supondremos que está de�nida en el conjunto (0;1)�V�M,
donde V es el conjunto abierto y conexo al cual esta restringido el movimiento de la partícula
Browniana y M es el conjunto abierto y conexo de posibles momentos (en realidad podemos
tomar f de�nida en (0;1)�V� R). La función f satisface la normalización
ZV�M
f (t; x; p) dxdp = 1 8t > 0; (3.3)
además de la ecuación de Fokker-Planck
@f
@t+
p
M
@f
@x+ F (t; x)
@f
@p=
@
@p
�� (p) pf +D (p)
@f
@p
�8 (t; x; p) 2 (0;1)�V�M. (3.4)
Las condiciones iniciales X0 = x0 y P0 = p0 de la ecuación de Langevin se convierten en las
condición de frontera
f (0; x; p) = � (p� p0) � (x� x0) 8 (x; p) 2 V�M. (3.5)
Debido a que (3.4) es una ecuación diferencial parcial de primer orden en el tiempo, se tiene
45
que la ecuación de Langevin-backward (3.1) describe un proceso Markoviano.
Físicamente, la cantidad in�nitesimal f (t; x; p) dxdp se interpreta como la probabilidad de
encontrar la partícula al tiempo t en el subconjunto in�nitesimal [x; x+ dx] � [p; p+ dp] del
espacio fase V � M. Esta interpretación y las ecuaciones (3.3) y (3.5) son las mismas para
cualquier integral, siempre y cuando la densidad de (Xt; Pt) exista.
3.1.4. Teorema de �uctuación-disipación
A partir de la ecuación de Fokker-Planck y utilizando la situación de equilibrio térmico es
posible encontrar igualdades que relacionen a los coe�cientes � y D de la ecuación (3.1) entre sí.
Tales ecuaciones son comúnmente denominadas relaciones de �uctuación-disipación de Einstein
y juegan un papel fundamental en el desarrollo de esta tesis.
A continuación se calcula la relación de Einstein para la ecuación de Langevin-backward
en ausencia de campos de fuerza externos, o sea F � 0. Sea � : (0;1) � V ! R+ � f0g la
densidad marginal del momento de f (siguiendo la notación de la sección anterior), es decir
� (t; p) =RV f (t; x; p) dx 8 (t; p) 2 (0;1) � M: En ausencia de fuerzas externas, la ecuación
(3.4) implica que � satisface la ecuación de Fokker-Planck
@�
@t(t; p) =
@
@p
�� (p) p� (t; p) +D (p)
@�
@p(t; p)
�8 (t; p) 2 (0;1)�M:
Esta ecuación tiene asociada una solución estacionaria �1; dada por
�1 (p) = N exp
��Z p
�p�q� (q)
D (q)dq
�; (3.6)
donde N es una constante de normalización y �p� 2M es un real arbitrario tal que la integral
dentro de la función exponencial existe.
La forma de la relación (3.6) implica que es posible generar distribuciones estacionarias
casi arbitrarias para el momento eligiendo a los coe�cientes � y D de manera adecuada. Para
ilustrar esto, considérese una densidad de probabilidad b� : M! R+ � f0g diferenciable. Se
desea encontrar una relación entre � y D de tal manera que la solución �1 coincida con la
función b�. Igualando b� (p) con �1 (p), de la ecuación (3.6) se sigue que
46
� (p)
D (p)p = � d
dplog b� (p) 8p 2M. (3.7)
La igualdad anterior constituye una relación general de �uctuación-disipación de Einstein.
Debido a que la partícula se encuentra en equilibrio térmico con el baño circundante y no
tiene restricciones en su movimiento (es decir, V = R), el momento de la partícula Browniana
tiene como densidad estacionaria a una función Maxwell de la forma
bgM (p) = (2�M=�)�1=2 exp���p2= (2M)
�; 8p 2 R; (3.8)
tomando � := (kBT )�1, donde T es la temperatura del baño y kB es la constante de Boltzmann.
Para que los coe�cientes � y D capturen la información del equilibrio térmico y la no restricción
del movimiento, de la ecuación (3.7) y (3.8) se sigue que � y D deben cumplir
D (p) = � (p)MkBT = � (p)M=� 8p 2 R: (3.9)
Esta ecuación es conocida como relación de �uctuación-disipación de Einstein.
Las relaciones de �uctuación-disipación (3.7) y (3.9) �jan solo uno de los dos coe�cientes
� (p) y D (p) : Así, existe la libertad de adaptar, por ejemplo, a la función � (p) para que el
proceso (3.1) exhiba el comportamiento de relajación correcto. Esta libertad es una de las
razones principales por las que el enfoque de Langevin es existosamente aplicado en un amplio
rango de procesos de termalización. Algunas expresiones físicamente razonables para � pueden
deducirse a partir de la teoría cinética o de modelos microscópicos hamiltonianos que tienen en
cuenta las interacciones al igual que las propiedades estadísticas del baño térmico. Ejemplos de
métodos para encontrar expresiones para � y D se discuten en la Sección 3.2.
Con miras a nuestro tratamiento del movimiento Browniano relativista, conviene tener en
cuenta que las ecuaciones de Langevin constituyen una herramienta útil para construir modelos
de movimiento Browniano con distribuciones arbitrarias para el momento y la velocidad de la
partícula Browniana.
47
3.1.5. El proceso de Ornstein-Uhlenbeck no-relativista
El proceso de Ornstein-Uhlenbeck clásico constituye el ejemplo estándar de un modelo
de movimiento Browniano (no-relativista) en ausencia de fuerzas externas. Este proceso esta
de�nido por la elección de coe�cientes constantes positivos en la ecuación de Langevin-Itô. Mas
precisamente, se eligen
� (p) � �0 > 0, D (p) � D0 > 0;
dando lugar a las siguientes ecuaciones diferenciales estocásticas para la posición y el momento
de la partícula Browniana
dXt =PtMdt; (3.10a)
dPt = ��0Ptdt+ (2D0)1=2 � dWt: (3.10b)
La ecuación diferencial anterior tiene solución explícita única.
Lema 52 La ecuación diferencial (3.10) tiene solución única, dada por
Xt = X0 +
Z t
0
PsMds;
Pt = P0e��0t + (2D0)
1=2 e��0tZ t
0e�0s � dWs;
para todo t 2 [0;1) :
Observación 53 En la literatura matemática, al proceso del lema anterior se le conoce como
el proceso de Ornstein-Uhlenbeck.
A partir de las propiedades básicas de la integral de Itô y teniendo en cuenta las condiciones
iniciales deterministas de la ecuación de Langevin, se encuentra fácilmente la media y el segundo
momento estadístico de los procesos de posición y momento.
Lema 54 Supongamos que la ecuación (3.10) se cumple junto con las condiciones iniciales
48
deterministas X0 = x0 y P0 = p0. Entonces, los procesos Pt y Xt cumplen que
hPti = P0e��0t;
P 2t�= P 20 e
�2�0t +D0�0
�1� e�2�0t
�;
hXt �X0i =P0�0M
�1� e��0t
�;D
[Xt �X0]2E
=2D0
(�0M)2 t+
�P0�0M
�2 �1� e��0t
�2+
D0�30M
2
��3 + 4e��0t � e�2�0t
�;
para todo t 2 [0;1) :
Asociada al desplazamiento cuadrático medio esta una cantidad de suma importancia en
cualquier teoría del movimiento Browniano: la constante de difusión asintótica. Esta constante
se de�ne para un proceso (Xt; Pt) de la forma (3.1) por
D1 (Xt; Pt) = l��mt!1
1
2t
D[Xt �X0]2
E:
La constante de difusión asintótica mide la capacidad efectiva que tiene la partícula Browniana
para desplazarse a través del medio; también permite estimar los procesos de relajación dentro
del sistema y aplicar resultados de la teoría de difusiones.
En el caso del proceso de Ornstein-Uhlenbeck la constante de difusión asintótica tiene la
siguiente expresión simple.
Lema 55 El proceso de Ornstein-Uhlenbeck (3.10) tiene asociada la constante de difusión asin-
tótica
D1 =D0
(�0M)2 :
Como los coe�cientes del proceso de Ornstein-Uhlenbeck son constantes, se cumple que los
coe�cientes del proceso equivalente con la integral backward no cambian, es decir el proceso
dPt = ��0Ptdt+ (2D0)1=2 � dWt;
es equivalente al proceso de Ornstein-Uhlenbeck (3.10b). Por lo tanto, como la partícula Brow-
niana está en equilibrio térmico con el baño circundante, la relación de Einstein (3.9) D0 =
49
�0MkBT se cumple también en este caso. Luego, la constante de difusión asintótica puede
escribirse como
D1 =kBT�0M
:
Debido a que la ecuación para el momento del proceso de Ornstein-Uhlenbeck tiene coe�-
cientes constantes, del Teorema 109 y el Capítulo 3.5 del libro de McKean [46] se sigue que la
densidad marginal del momento satisface la ecuación de Fokker-Planck
@�
@t(t; p) =
@
@p
��0p� (t; p) +D0
@�
@p(t; p)
�8 (t; p) 2 (0;1)� R: (3.11)
Considerando la condición inicial determinista � (0; p) = � (p0 � p), la solución de la ecuación
(3.11) esta dada por
� (t; p) =
��0
2�D0 (1� e�2�0t)
�1=2exp
"��0�p� p0e��0t
�22D0 (1� e�2�0t)
#8p 2 R.
En el límite t!1 la solución anterior se reduce a una densidad Gaussiana
�1 (p) =
��02�D0
�1=2exp
���0p
2
2D0
�8p 2 R;
es decir que la distribución estacionaria marginal de momentos del proceso de Ornstein-Uhlenbeck
corresponde a una distribución normal, lo cual nos dice que es un modelo de movimiento Brow-
niano que atrapa la información del equilibrio térmico, algo bastante deseable.
Dada una densidad marginal b� (t; p) para el momento de la partícula Browniana y teniendoen cuenta la relación entre velocidad y momento p =Mv, se tiene que la densidad b (t; v) quegobierna la velocidad de la partícula Browniana está dada por
b (t; v) := dp
dvb� (t; p (v)) =Mb� (t; p (v)) .
Entonces, en el caso del proceso de Ornstein-Uhlenbeck, la densidad (t; v) de la velocidad de
la partícula Browniana está dada por
(t; v) =
��0
2�D0 (1� e�2�0t)
�1=2exp
"��0�Mv � p0e��0t
�22D0 (1� e�2�0t)
#8 (t; v) 2 R:
50
Teniendo en cuenta la relación de Einstein (3.9), podemos ver que la densidad estacionaria para
la velocidad de la partícula Browniana se expresa como
M (v) =
�M
2�kBT
�1=2exp
�� Mv2
2kBT
�8v 2 R;
la cual es una distribución de Maxwell para la velocidad. Esto vuelve a con�rmar que el modelo
de movimiento Browniano atrapa la información del equilibrio térmico.
Los modelos de Ornstein-Uhlenbeck con in�uencia de campos de fuerza externos fueron es-
tudiados intensamente durante el siglo pasado desde un enfoque de Kolmogorov. Sin embargo,
para campos de fuerza arbitrarios F (t; x) generalmente es muy difícil, y en muchos casos im-
posible, encontrar soluciones a la ecuación de Fokker-Planck. En el caso más simple se tienen
campos de fuerza conservativos y estacionarios F (t; x) � F (x) con potencial � (x), es decir
F (x) = � @
@x� (x) :
En este caso es posible determinar soluciones estacionarias en el límite t!1. La ecuación de
Fokker-Planck para la densidad g (t; x; p) de un proceso de Ornstein-Uhlenbeck bajo la in�uencia
de un campo de fuerza F (t; x) está dada por
@g
@t+
p
M
@g
@x+ F (t; x)
@g
@p=
@
@p
��0pg +D0
@g
@p
�:
Imponiendo la relación de Einstein D0 = �0MkBT , la solución estacionaria de la ecuación
anterior está dada por la densidad de Maxwell-Boltzmann
g (x; p) = Z�1 exp����p2
2M+�(x)
��; � = (kBT )�1 ;
donde Z�1 es una constante de normalización.
3.2. Modelos microscópicos: el modelo de colisión binaria
Cuando se consideran ecuaciones de Langevin del tipo (3.1), en principio se pueden distinguir
entre los dos siguientes enfoques:
51
(a) Postular estas ecuaciones como modelos fenomenológicos, estudiar las consecuencias mate
máticas y comparar las predicciones resultantes con experimentos con el �n de validar o
no la teoría. Adoptando este enfoque, los parámetros y la forma explícita de las funciones
de ruido y fricción tienen que ser determinados de la información experimental.
(b) Alternativamente, se puede tratar de motivar y derivar las ecuaciones de Langevin a partir
de modelos microscópicos. Si se tiene éxito, este enfoque permite encontrar expresiones
explícitas para las funciones de ruido y fricción en términos de los parámetros del modelo
microscópico.
En esta sección se considera el segundo enfoque, el cual atrajo considerable interés en las úl-
timas décadas. Las ecuaciones de Langevin proveen una descripción aproximada de la dinámica
microscópica �exacta�. Por lo tanto, con el �n de derivar ecuaciones diferenciales estocásticas
a partir de modelos microscópicos, se tienen que imponer ciertas aproximaciones; las cuales
determinan el rango de aplicabilidad del enfoque de Langevin para el movimiento Browniano.
Generalmente, se pueden seguir al menos dos caminos distintos para derivar ecuaciones dife-
renciales estocásticas de la forma (3.1) de modelos más precisos:
(1) Empezando con una ecuación tipo Boltzmann o una ecuación maestra para la densidad
de probabilidad de la partícula Browniana, se puede tratar de reducir estas ecuaciones
integro-diferenciables a una igualdad de Fokker-Planck realizando aproximaciones ade-
cuadas. Una vez encontrada la ecuación de Fokker-Planck es sencillo escribir la ecuación
de Langevin asociada. La dinámica de las colisiones microscópicas queda implícita en la
ecuación de Boltzmann.
(2) Alternativamente, se puede partir de un modelo microscópico que describa la interacción
entre la partícula Browniana y el baño térmico. Al eliminar los grados de libertad del
baño térmico en las ecuaciones del movimiento de la partícula Browniana, se obtiene una
ecuación de Langevin generalizada, la cual puede reducirse a la forma (3.1) en ciertos
casos límite. Como subproducto, con este enfoque las relaciones de �uctuación-disipación
surgen de una manera más natural al asumir una distribución de probabilidad para la
con�guración (inicial) del baño.
52
En esta sección se desarrolla el segundo procedimiento con un modelo de colisión binaria, el
cual puede ser extendido al caso relativista. Cabe destacar que el método más empleado para
deducir expresiones para las funciones de ruido y fricción se deriva del modelo del oscilador
armónico; sin embargo éste no puede extenderse a la relatividad, razón por la cual no se expone
en esta tesis.
3.2.1. El modelo de colisión binaria elástica
El modelo de colisión binaria parte de la idea de que el movimiento estocástico de una
partícula Browniana se genera debido a las frecuentes colisiones elásticas con las partículas
del baño térmico. Básicamente, este modelo consiste en construir una ecuación de Langevin a
partir de la cinemática de las colisiones de la partícula Browniana y las partículas del baño. El
modelo de colisión binaria se usa para encontrar expresiones útiles y realistas para las funciones
de ruido y fricción, además de que permite dar una justi�cación heurística de la ecuación de
Langevin (3.1).
En la exposición de esta sección, se considera que el baño térmico consiste de N 2 N
partículas de masa m, tomando N � 1 y m � M . Además, se denota por xr (t) y pr (t) 2
a la posición y el momento al tiempo t de la partícula �r�, r 2 f1; : : : ; Ng del baño térmico,
respectivamente.
Cinemática de las colisiones
Una colisión elástica entre la partícula Browniana y la partícula �r�del baño térmico está
gobernada por las leyes de conservación
E + �r = bE +b�r, P + pr = bP + bpr;donde E (P ) := P 2=2M es la energía de cinética de la partícula Browniana, �r (p) := (p)2=2m
es la energía de la partícula �r�del baño y los símbolos con (b) se re�eren al estado posterior ala colisión. Considerando las ecuaciones anteriores, se encuentra que la ganancia de momento
2Resulta importante notar que xr (t) y pr (t) también son procesos estocásticos. Para ganar un poco declaridad, no usamos la notación con subínices acostumbrada para los procesos.
53
de la partícula Browniana en la colisión con la partícula �r�está dada por
�Pr (t) := bPt � Pt = � 2m
M +mPt +
2M
M +mpr (t) : (3.12)
Con objeto de justi�car heurísticamente la ecuación de Langevin (3.1) y de encontrar coe�-
cientes � y D realistas, conviene considerar el cambio de momento �P (t) := Pt+�t � Pt de la
partícula Browniana en el intervalo de tiempo �mesoscópico� [t; t+ �t], haciendo las siguientes
tres suposiciones:
las colisiones que ocurren en [t; t+ �t] pueden ser vistas como eventos independientes;
el salto de tiempo �t es lo su�cientemente pequeño para que se cumpla j�P (t) =Ptj � 1 y
para que ocurra a lo más sólo una colisión entre la partícula Browniana y una partícula
especí�ca del baño térmico;
el valor de �t es lo su�cientemente grande para que el número de colisiones sea mayor a
1.
Estos requisitos solo pueden cumplirse simultáneamente si m�M , lo cual fue un supuesto
ya hecho. Con los supuestos anteriores es posible aproximar �P (t) por
�P (t) �NXr=1
�Pr (t) Ir (t; �t) ; (3.13)
donde Ir (t; �t) es la función indicadora para una colisión entre la partícula Browniana y la
partícula �r�del baño térmico durante el intervalo [t; t+ �t], es decir
Ir (t; �t) =
8>>><>>>:1 si una colisión entre la partícula Browniana
y la partícula �r�ha ocurrido en [t; t+ �t] ;
0 en otro caso.
En el caso unidimensional, Ir (t; �t) puede expresarse como
Ir (t; �t) = � (Xt � xrt )��exrt � eXt�+�(xrt �Xt)�� eXt � exrt� ; (3.14)
54
siendo � : R! R la función Heaviside y
eXt = Xt + Vt�t, exrt = xr (t) + vr (t) �t (3.15)
las posiciones proyectadas al tiempo t+�t de las partículas en colisión. Así, podemos expresar la
función indicadora Ir en términos de las posiciones y las velocidades de las partícula Browniana
y partícula �r�del baño, algo que resulta útil en la determinación de los coe�cientes � y D.
La ecuación (3.14) puede justi�carse al tener en cuenta que dos partículas colisionan en el
intervalo [t; t+ �t] cuando se dirigen la una contra la otra y la distancia que las separa es lo
su�cientemente pequeña.
Combinando las igualdades (3.12) y (3.13) se llega a que
�P (t) � �2"NXr=1
m
M +mIr (t; �t)
#Pt + 2
NXr=1
m
M +mpr (t) Ir (t; �t) ; (3.16)
donde, adicionalmente, se supone que para cada colisión que ocurre en el intervalo [t; t+ �t] ;
el momento de la partícula Browniana antes de la colisión es aproximadamente igual al valor
�inicial�Pt. En vista de que m�M , la ecuación anterior puede simpli�carse como
�P (t) � �2"NXr=1
m
MIr (t; �t)
#Pt + 2
NXr=1
pr (t) Ir (t; �t) : (3.17)
Una comparación entre la ecuación (3.1) y la igualdad (3.17) sugiere que, heurísticamente,
el primer término del lado derecho de la ecuación (3.17) puede interpretarse como fricción,
mientras que el segundo término puede pensarse como ruido. Esto justi�ca, heurísticamente, el
porque tomar ecuaciones de Langevin de la forma (3.1) con F � 0 para modelar la dinámica
de la partícula Browniana en ausencia de fuerzas externas. A pesar de que es bastante similar
a la ecuación de Langevin (3.1), la ecuación (3.17) es considerablemente más complicada que
la ecuación de Langevin. Esto se debe a que las funciones indicadoras Ir de la ecuación (3.17)
dependen tanto de la posición y la velocidad de la partícula Browniana, como de la posición y
la velocidad de las partículas del baño. Sin embargo, es posible calcular las propiedades estadís-
ticas de los incrementos de momento �P (t) de la ecuación (3.17), suponiendo una distribución
especí�ca para las partículas del baño térmico.
55
Distribución del baño térmico
En principio, se puede utilizar las ecuaciones (3.16) y (3.17) para calcular los momentos
estadísticosD[�P (t)]j
Epara una densidad arbitraria fN (x1; p1; : : : ; xN ; pN ) de las partículas
del baño térmico. Aquí se trata la situación cuando el baño térmico (arbitrariamente grande)
está dado por un gas quasi-ideal, el cual se encuentra en equilibrio térmico con su ambiente.
En este caso, la densidad que modela la probabilidad de encontrar a la partícula �r�del baño
térmico en algún lugar del espacio fase, está dada por la función espacialmente homogénea de
Maxwell
f1 (xr; pr) = (2�mkBT )�1=2 L�1 exp�
�p2r(2mkBT )
�8 (x; p) 2 (0; L)�R; (3.18)
donde L es el volumen (unidimensional) contenedor del baño. Además, se supondrá que:
las partículas del baño térmico son independiente e idénticamente distribuidas;
las colisiones entre la partícula Browniana y las partículas del baño térmico no afectan la
distribución de las partículas.
Las dos suposiciones anteriores pueden justi�carse para baños lo su�cientemente grandes,
siempre y cuando las colisiones entre las partículas del baño restablecen rápidamente una dis-
tribución espacialmente homogénea.
La fuerza media de deriva
La deriva media (al tiempo t) se de�ne como el cambio promedio de momento sobre el
intervalo [t; t+ �t], dado el valor del momento P al tiempo t; es decir la deriva media al tiempo
t es la esperanza h�P (t) j Pt = pi. En el caso de la ecuación (3.17), se encuentra que
h�P (t) j Pt = pi = �2N�mM
�hIr (t; �t)i p+ 2N hprIr (t; �t)i .
Debido a que la densidad en el espacio fase de las partículas del baño es una función de la
forma (3.18) y a que la función indicadora Ir se puede expresar como (3.14), se tiene que para
56
cualquier función continua G : (0; L)� R! R y para todo r 2 f1; : : : ; Ng se cumple que
hG (xr; vr) Ir (t; �t)i =
Z L
0
Z 1
�1G (xr; vr)� (Xt � xrt )�
�exrt � eXt� (vr) dvrdxr +Z L
0
Z 1
�1G (xr; vr)� (x
rt �Xt)�
� eXt � exrt� (vr) dvrdxr; (3.19)donde : R! R+ es la función de densidad de probabilidad para la velocidad de una partícula
del baño térmico, la cual es una función de Maxwell de la forma
(vr) =�v2B�
��1=2exp
�v2rv2B
�; vB :=
�2kBTm
�1=2: (3.20)
Cabe destacar que las integrales en la igualdad (3.19) son difíciles de calcular, sino imposi-
bles. Sin embargo, existen algunas identi�caciones heurísticas y métodos que permiten encon-
trar expresiones útiles. Por ejemplo, puede utilizarse una expansión de Taylor con derivadas
ge-neralizadas. De las ecuaciones (3.14) y (3.15), se sigue que la función indicadora Ir puede
ser aproximada por su expansión de Taylor (generalizada) hasta primer orden
Ir (t; �t) ��t
2jvr (t)� Vtj � (xr (t)�Xt) ;
donde � (�) es una delta de Dirac. Usando la última expresión, es posible encontrar que
hG (xr; vr) Ir (t; �t)i =�t
2L
Z 1
�1G (X; vr) jvr (t)� Vtj (vr) dvr;
tomando como en (3.20). A partir de la igualdad anterior, se llega a que
h�P (t) j Pt = pi � �2nbkBT(��1=2
�p
pB
�exp
"��p
pB
�2#)�t
�2nbkBT("�
p
pB
�2+1
2
#erf
�p
pB
�)�t; (3.21)
donde nb = N=L es la densidad de las partículas del baño térmico, pB :=MvB =M (2kBT =m)1=2
57
es una escala característica térmica para el momento y erf (z) es la función error de�nida por
erf (z) =2p�
Z z
0e�x
2dx:
Si bien obtener una expresión clara (sin aproximaciones de ningún tipo) para la deriva
media no es factible, se pueden obtener expresiones útiles a través de métodos numéricos como
aproximación por splines, entre muchos otros procedimientos.
Asociada a la deriva media se encuentra asociada la fuerza media de deriva K (p), la cual se
de�ne como
K (p) :=��P (t)
�tj Pt = p
�:
Esta fuerza puede considerarse como una medida de los procesos de relajación de la partícula
Browniana y el baño térmico. Al realizar aproximaciones numéricas podemos ver que la fuerza
media de deriva crece linealmente para valores pequeños del momento y de manera cuadrática
para valores grandes del momento.
A continuación se estudia un procedimiento para aproximar el modelo de la igualdad (3.17)
por una ecuación diferencial estocástica no lineal del tipo (3.1).
Aproximación de Langevin
La ecuación para los incrementos del momento (3.17) no es aún una ecuación de Langevin.
Las ecuaciones diferenciales estocásticas del tipo
dPt = �� (Pt)Ptdt+ [2D (Pt)]1=2 � dWt; (3.22)
constituyen modelos fenomenológicos que proveen una descripción simpli�cada de la dinámica
microscópica, la cual en el caso del modelo de colisión binaria está descrita más precisamente
por la ecuación (3.17). Por lo tanto, para obtener un modelo de Langevin útil, se deben elegir
a los coe�cientes � y D de tal manera que produzcan la mejor aproximación posible en la clase
de ecuaciones diferenciales estocásticas de�nida por la ecuación (3.22). Es posible que existan
una in�nidad de criterios para de�nir una �buena aproximación�, en este caso los criterios que
deberá cumplir el proceso (3.22) para decir que es una buena aproximación son dos:
58
Aproximar la distribución estacionaria correcta;
Exhibir el comportamiento de relajación media correcto.
El primer criterio es equivalente a imponer la relación de �uctuación-disipación apropiada
sobre los coe�cientes � y D. Para el modelo de colisión binaria considerado aquí, la distribución
estacionaria para el momento de la partícula Browniana esta dada por una densidad de Maxwell
�1 (p) = (2�MkBT )�1=2 exp�
�p2(2MkBT )
�8p 2 R:
De acuerdo a lo hecho en la Sección 3.1, esto implica que � y D deben satisfacer la relación de
Einstein
D (p) = � (p)MkBT 8p 2 R: (3.23)
El segundo criterio puede expresarse matemáticamente como
�dP (t)
dtj Pt = p
�=
��P (t)
�tj Pt = p
�;
donde se toma heurísticamente3 dP (t)dt = � (Pt) + D0 (Pt) +
[2D(Pt)]1=2�dWt
dt y se supone queR t0 h2D (Ps)i ds < 1 8t � 0. Teniendo en cuenta la relación de Einstein (3.23) y la relación
heurística anterior, se llega a que el lado izquierdo de la ecuación anterior esta dado por
�dP (t)
dtj Pt = p
�� �
�� (p) p� �0 (p)MkBT
�:
Así, el segundo criterio es equivalente a la ecuación diferencial ordinaria (no-estocástica)
�� (p) p+ �0 (p)MkBT = K (p) : (3.24)
3Nótese que el proceso Pt no necesariamente es diferenciable, por lo que estrictamente hablando dPtdtno tiene
sentido. Sin embargo, podemos cuanti�car la variación del proceso Pt heurísticamente. La ecuación (3.22) esequivalente al proceso
dPt =�� (Pt)�D0 (Pt)
�dt+ [2D (Pt)]
1=2 dWt:
Al dividir por dt, el cual será pensado como un escalar, se obtiene la identi�cación heurística. Además, se suponeR t0h2D (Ps)i ds <1 8t � 0 para poder asegurar que la esperanza
DR t0[2D (Ps)]
1=2 � dWs
E= 0 existe 8t � 0.
59
Con respecto a los dos criterios formulados anteriormente, la solución a la ecuación diferen-
cial ordinaria (3.24) proporciona el coe�ciente de fricción � que genera la �mejor�aproximación
de Langevin para el modelo de colisión binaria; la condición inicial para � (p) debe especi�carse
de tal manera que se obtenga el comportamiento asintótico correcto. La información sobre el
modelo de colisión binaria y la distribución del baño térmico queda implícita en la fuerza media
de deriva K (p). Evidentemente, el procedimiento que permite llegar a la ecuación (3.24) puede
generalizarse para otros modelos de interacción y otras distribuciones para el baño térmico,
siempre y cuando se conozca la distribución estacionaria para el momento de la partícula Brow-
niana. Otros tipos de interacción (e.g. colisiones no elásticas) podrían producir otra función
K (p); mientras una distribución no Maxwelliana para el baño térmico podría afectar no sólo
el lado derecho de la ecuación (3.24), también modi�caría el lado izquierdo pues la relación de
�uctuación-disipación sería distinta.
Desafortunadamente, casi siempre es muy difícil, o incluso imposible, encontrar la solu-
ción analítica de la ecuación diferencial ordinaria (3.24) para una función K (p) realista. Para
propósitos prácticos, se pueden obtener aproximaciones útiles, por ejemplo, al considerar el
comportamiento asintótico para p ! 0 o p ! 1. En el caso del modelo de colisión binaria, se
utiliza la aproximación
� (p) ' �K (p)p
:= �1 (p) :
Esta aproximación se toma como una opción sencilla de calcular y puede justi�carse para
algunos baño térmicos y relaciones de �uctuación, sin embargo no es totalmente cierta en todos
lo casos.
Para el proceso de Ornstein-Uhlenbeck, el valor de la constante �0 puede aproximarse por
�0 := �1 (0); según la ecuación (3.21) se tiene
�0 := �1 (0) = �K0 (0) � nbm
M
�8kBT�m
�1=2:
Otro valor que podría tomarse para �0 es el valor mínimo o máximo (si existen) de �1 (p) :
Cabe mencionar que en el caso de la aproximación (3.21) el valor mínimo de �1 (p) coincide con
�1 (0). El valor de D0 se obtiene a partir de la relaciones de Einstein D0 � D (p) = �0MkBT .
Adoptando estas simpli�caciones adicionales, puede esperarse que el proceso de Ornstein-
60
Uhlenbeck correspondiente genere una descripción útil para partículas Brownianas a tempera-
turas T lo su�cientemente bajas.
De una perspectiva más general, los ejemplos anteriores ilustran los objetivos que no es
posible alcanzar por los modelos fenomenológicos de Langevin basados en procesos de Wiener.
Las ecuaciones de Langevin de este tipo, y sus respectivas ecuaciones de Fokker-Planck, proveen
una descripción simpli�cada de la dinámica microscópica inherente. Las funciones coe�ciente en
la ecuación de Langevin/Fokker-Planck permiten construir procesos estocásticos que exhiben
el mismo comportamiento de relajación y aproximan la misma distribución estacionaria que
el proceso físico real. Las distribuciones estacionarias a menudo pueden inferirse de principios
termoestadísticos (e.g. máxima entropía), mientras que el comportamiento de relajación debe
deducirse de la dinámica microscópica exacta. En muchos casos, los modelos estocásticos resul-
tantes son su�cientes para compararlos con la información experimental accesible, pero pueden
volverse inexactos para describir correlaciones de mayor orden y/o la dinámica de relajación no
próxima al estado asintótico.
61
Capítulo 4
El Movimiento Browniano
Relativista en el Espacio Fase
El objetivo de este capítulo es estudiar y establecer la teoría de Langevin para movimien-
tos Brownianos en el marco de la teoría especial de la relatividad. Se consideran ecuaciones
diferenciales estocásticas que describen procesos de Markov en el espacio fase de una partícula.
La razón principal para considerar procesos en el espacio fase, en lugar de procesos que solo
involucren a la posición Xt, se debe a que es imposible de�nir un proceso de Markov relativista
no-trivial Xt en el espacio-tiempo de Minkowski; como lo ha demostrado R. Dudley [14] en
1966. Los desarrollos de esta sección pueden considerarse la contraparte relativista de aquellos
tratados en el Capítulo 3.
Desde una perspectiva conceptual, pueden distinguirse dos enfoques complementarios en
la modelación de procesos estocásticos que retratan la dinámica de una partícula Browniana
relativista. El primero de ellos consiste en postular ecuaciones de evolución para la densidad de
las probabilidades de transición del proceso estocástico; el ejemplo más típico son ecuaciones
diferenciales parciales del tipo Fokker-Planck. En un marco relativista, la cota en la velocidad
de propagación de las partículas impone restricciones bastante estrictas en la estructura de
las ecuaciones de evolución. Alternativamente, un segundo enfoque es comenzar postulando
ecuaciones diferenciales estocásticas como modelos fenomenológicos y, posteriormente, derivar
ecuaciones de evolución. En el presente capítulo se expone esta última vía, la cual tiene como
63
ventaja que el origen físico de la mecánica estocástica está �menos escondido�en comparación
con otros enfoques.
En la primera sección de este capítulo se establece la ecuación de Langevin relativista, se
discuten algunas diferencias con el caso no-relativista y se calculan las ecuaciones de Fokker-
Planck para la densidad de la trayectoria en el espacio fase. En la sección 2 se trata el caso
de un baño térmico estacionario, isotrópico y homogéneo, el cual está en ausencia de fuerzas
externas; además se establece un teorema de �uctuación-disipación para este caso.
4.1. Ecuaciones de Langevin y Fokker-Planck relativistas
Cuando se consideran ecuaciones de Langevin como modelos del movimiento Browniano,
implícitamente se asume que es posible y razonable separar los grados de libertad de la partícula
Browniana de los del ambiente (baño térmico). Desde este punto de vista, se pueden especi�car
dos marcos de referencia distinguibles: el marco del laboratorio �, de�nido como el marco
de referencia inercial del baño térmico, y el marco inercial �� que se está comoviendo con la
partícula Browniana en un instante de tiempo dado.
Para comenzar se considera la cuestión de como modelar la dinámica de una partícula
Browniana relativista en el marco del laboratorio �, el cual se supondrá que tiene d 2 N
dimensiones espaciales.
4.1.1. La ecuación de Langevin relativista: principios generales de construc-
ción
Considérese el movimiento de una partícula Browniana relativista puntual en �. Supóngase
que tal partícula tiene masa en reposo M > 0 y que está inmersa en un baño térmico que
consiste, por ejemplo, de partículas líquidas más pequeñas (masa en reposo m � M) a tem-
peratura constante T . Las cantidades M , m y T corresponden a las medidas en el marco de
referencia inercial �. La posición, la velocidad, el momento relativista y la energía relativista
de la partícula Browniana en � para cualquier tiempo t de �, se toman como Xt, Vt := dXtdt ,
Pt := M (Vt)Vt y P 0t :=�M2 +P2t
�1=2, respectivamente. Debido al comportamiento aleato-rio de la partícula Browniana, se asume que fXt; t � 0g, fPt; t � 0g, fVt; t � 0g y
�P 0t ; t � 0
64
son procesos estocásticos en Rd de�nidos en algún espacio de probabilidad.
Todos los modelos estocásticos para la posición de la partícula Browniana relativista son
procesos estocásticos que satisfacen la siguiente de�nición.
De�nición 56 Sea fYt; t � 0g un proceso estocástico que toma valores en Rd y que es diferen-
ciable en todas partes respecto a la variable temporal. Diremos que Y es un proceso estocástico
relativista si satisface ����dYtdt
���� < c = 1 8t � 0.
Nótese que la posición Xt es un proceso estocástico relativista si, y sólo si,
jPtj�M2 +P2t
�1=2 = jPtjP 0t
= jVtj =����dXtdt
���� < 1 8t � 0; (4.1)
lo cual obviamente se cumple.
Observación 57 Bajo la de�nición de proceso estocástico relativista, queda establecida la propiedad
de que la partícula Browniana no puede ir más rápido que la luz.
De particular interés son los modelos del movimiento Browniano que pueden ser expresados
como ecuaciones diferenciales estocásticas similares a las ecuaciones de Langevin no-relativistas
(3.1). La idea básica para construir tales procesos consiste en asociar integrales estocásticas
solo al momento relativista P =�P i�de la partícula Browniana, el cual puede tomar valores
en todo Rd. Haciendo esto se asegura que la condición (4.1) se cumpla automáticamente.
Los modelos de movimiento Browniano relativista que se desarrollan en esta tesis son ecua-
ciones similares a la ecuación de Langevin no-relativista para la posición Xt y el momento Pt
de la partícula Browniana, parametrizadas en el tiempo de �. Antes de escribir explícitamente
ecuaciones diferenciales estocásticas para las componentes del momento relativista P =�P i�,
conviene dar una pequeña explicación del porque se decidió parametrizar en términos del tiempo
de �.
La elección del parámetro temporal: ¿Por qué el tiempo de �?
Una suposición (postulado) fundamental de la física Galileana no-relativista es la existencia
de un tiempo universal. Por consiguiente, en la teoría no-relativista de Langevin, es bastante
65
natural identi�car al tiempo universal t con el parámetro temporal del proceso estocástico
conductor, el cual se toma como un proceso de Wiener multi-dimensional Wt. En contraste,
en relatividad especial la noción de tiempo se vuelve marco-dependiente y, por lo tanto, es
importante especi�car cual parámetro de tiempo se utiliza para cuanti�car las �uctuaciones del
proceso estocástico conductor.
Cuando se considera el movimiento estocástico de una partícula Browniana, pueden distin-
guirse dos parámetros temporales característicos: el tiempo t del marco de referencia inercial
�, el cual puede ser interpretado como el tiempo propio del baño térmico, y el tiempo propio
� de la partícula Browniana. En principio, cualquiera de los dos tiempos puede usarse para
formular ecuaciones diferenciales estocásticas para las componentes espaciales del momento re-
lativista de la partícula, P =�P i�. Sin embargo, en el tratamiento convencional de Langevin
del movimiento Browniano, usualmente se considera a la fricción y al ruido como fuerzas ex-
ternamente impuestas que actúan sobre la partícula Browniana y re�ejan �uctuaciones en el
baño térmico. Por lo tanto, parece más natural caracterizar las propiedades (estadísticas) de la
fuerza de Langevin en términos del tiempo del laboratorio t. Sin embargo, como se discute en
el Capítulo 6, es posible reparametrizar una ecuación de Langevin dada en términos del tiempo
propio asociado.
La ecuación de Langevin en el tiempo del laboratorio
Supóngase que la partícula Browniana está en presencia de un campo de fuerza determi-
nístico externo. La dinámica estocástica de la partícula Browniana relativista en � puede ser
descrita por la ecuación diferencial estocástica parametrizada por el tiempo de �
dXit =
P itP 0t
dt; (4.2a)
dP it = F i (t;Xt;Pt) dt� aij (t;Xt;Pt)P jt dt+ cij (t;Xt;Pt)� dWjt , (4.2b)
donde:
a) Xit y P
it i = 1; : : : ; d, denotan las componentes espaciales en � de la posición y el momento
relativista;
b) � 2 f�; �; �g señala el tipo de integral estocástica considerado;
66
c) Wt =�W 1t ; ::;W
dt
�es un proceso de Wiener d-dimensional;
d) F i : [0;1) � V � Rd ! R, aij : [0;1) � V � Rd ! R y cij : [0;1) � V � Rd ! R
i; j = 1; : : : ; d son funciones C1, con V el conjunto abierto y conexo al cual esta restringido
el movimiento de la partícula Browniana;
e) F i denota el campo de fuerza un campo de fuerza externo determinístico,�aij (t;Xt;Pt)P jtes una fuerza de fricción fenomenológica y el término cij (t;Xt;Pt)� dW j
t representa una
fuerza de Langevin relativista.
A la igualdad (4.2) se le denomina ecuación de Langevin relativista (en el marco de referencia
�). En el desarrollo de la tesis, se trata sólo el caso�F i�� 0. Por ser F i, aij y cij de clase
C1, la ecuación (4.2) tiene sentido matemático para los tres tipos de integral considerados. La
existencia y unicidad de la solución de la ecuación de Langevin relativista (4.2) depende de la
forma de las funciones F i, aij y cij ; así como de la condición inicial considerada. En el Capítulo
2 se establecieron algunos resultados que tratan claramente este tema. A lo largo de esta tesis, se
supondrá que la ecuación de Langevin (4.2) tiene solución única; solo se demostrará la existencia
y unicidad para los ejemplos clásicos de modelos del movimiento Browniano relativista. Además
se supondrá que las condiciones inicialesX0 y P0 de la ecuación (4.2) son deterministas, a menos
que se indique lo contrario.
Observación 58 Si además de la condiciones ya mencionadas, la ecuación de Langevin (4.2)
cumple la restricción de crecimiento del Teorema 29 de existencia y unicidad y tiene una condi-
ción inicial (X0;P0) independiente de Wt 8t � 0, tal igualdad tiene solución única.
Una posible motivación de la ecuación de Langevin relativista (4.2) es la analogía con el
caso no-relativista, sin embargo, esta es una razón muy pobre. Una motivación más adecuada de
esta ecuación y una forma de encontrar coe�cientes aij y cij realistas se discute en el Capítulo
7.
Las ecuaciones (4.2) gobiernan a las componentes espaciales de los cuadrivectores�X0t ;Xt
�y�P 0t ;Pt
�. Además, es posible agregar una ecuación para la componente temporal X0 poniendo
dX0t =
P 0tP 0t
dt = dt;
67
lo cual es consistente con el hecho de que X0t = t. La ecuación para la componente de energía
P 0t puede derivarse de la ecuación (4.2b) al aplicar las fórmulas de cambio de variable para
ecuaciones diferenciales estocásticas al proceso P 0t :=�M2 +P2t
�1=2.Proposición 59 Supongamos que se cumple la ecuación (4.2). Entonces, el proceso P 0t :=�M2 +P2t
�1=2 satisface la ecuacióndP 0t =
((F i (t;Xt;Pt)� aij (t;Xt;Pt)P jt )
P itP 0t
+ ��Dij (t;Xt;Pt)
2
"�ijP 0t
� P itPjt�
P 0t�3#)
dt
+P itP 0t
cir (t;Xt;Pt)� dW rt , (4.3)
donde Dij (t;Xt;Pt) := cir (t;Xt;Pt) cjr (t;Xt;Pt), �ij es una delta de Kronecker y, depen-
diendo de la integral estocástica, �� = 1, �� = 0 y �� = �1.
Demostración. Para demostrar lo deseado, usaremos la fórmula de Itô (26) y sus análogos
de Stratonovich (47) y backward (51); lo demostraremos primero para la integral de Itô. De�-
namos u : Rd �Rd � [0;1)! R dada por u (x;p;t) =�M2 + p2
�1=2: Nótese que @u
@t ,@u@xi
y @u@pi
existen 8i 2 f1; : : : ; dg y que
@u
@t(x;p;t) = 0;
@u
@xi(x;p;t) = 0;
@u
@pi(x;p;t) =
pi
(M2 + p2)1=2;
para todo�x1; : : : ; xd; p1; : : : ; pd; t
�= (x;p;t) 2 Rd � Rd � [0;1).
Por otro lado, de lo anterior se sigue que @2u@pj@pi
, @2u@xj@xi
; @2u@xj@pi
y @2u@pj@xi
existen para i; j =
1; : : : ; d y que
@2u
@ (pi)2(x;p;t) =
M2 + p2 ��pi�2
(M2 + p2)3=2;
@2u
@pj@pi(x;p;t) = � pjpi
(M2 + p2)3=2;
@2u
@ (xi)2(x;p;t) =
@2u
@xj@xi(x;p;t) =
@2u
@xj@pi(x;p;t) =
@2u
@pi@xj(x;p;t) = 0;
68
para cualesquiera�x1; : : : ; xd; p1; : : : ; pd; t
�= (x;p;t) 2 Rd � Rd � [0;1) e i; j 2 f1; : : : ; dg,
i 6= j.
De la de�nición de u y lo hecho arriba, es claro que u es continua y que @u@t ,@u@xi, @u@pi, @2u@pj@pi
,
@2u@xj@xi
; @2u@xj@pi
y @2u@pj@xi
son continuas para i; j = 1; : : : ; d. Entonces podemos aplicar la fórmula
de Itô al proceso P 0t :=�M2 +P2t
�1=2= u (Xt;Pt; t). Mediante la fórmula de Itô y evitando el
criterio de sumación de Einstein para ganar claridad, llegamos a que
dP 0t =dXi=1
@u
@pi(Xt;Pt; t) dP
it +
1
2
dXi;j=1
@2u
@pj@pi(Xt;Pt; t)
dXr=1
cir (t;Xt;Pt) cjr (t;Xt;Pt) dt
=
dXi=1
P itP 0t
dP it +1
2
264 dXi=1
0B@�P 0t �2 � �P it �2�P 0t�3 �
dXj=1i6=j
P jt Pit�
P 0t�3Dij (t;Xt;Pt)
1CA375 dt
=dXi=1
24F i (t;Xt;Pt)� dXj=1
aij (t;Xt;Pt)Pjt
35 P itP 0t
dt+dXi=1
P itP 0t
dXj=1
cij (t;Xt;Pt) � dW jt
+1
2
24 d
P 0t�
dXi;j=1
P jt Pit�
P 0t�3Dij (t;Xt;Pt)
35 dt=
dXi;j=1
P itP 0t
cij (t;Xt;Pt) � dW jt +
dXi=1
0@F i (t;Xt;Pt)� dXj=1
aij (t;Xt;Pt)Pjt
1A P itP 0t
dt
+
dXi;j=1
Dij (t;Xt;Pt)
2
�ijP 0t
� P jt Pit�
P 0t�3!dt:
Aplicando el criterio de sumación de Einstein a la ecuación anterior, tenemos que
dP 0t =P itP 0t
cij (t;Xt;Pt) � dW jt +
�F i (t;Xt;Pt)� aij (t;Xt;Pt)P jt
� P itP 0t
dt
+Dij (t;Xt;Pt)
2
"�ijP 0t
� P jt Pit�
P 0t�3#dt
=
(�F i (t;Xt;Pt)� aij (t;Xt;Pt)P jt
� P itP 0t
+Dij (t;Xt;Pt)
2
"�ijP 0t
� P jt Pit�
P 0t�3#)
dt
+P itP 0t
cij (t;Xt;Pt) � dW jt ;
69
lo que prueba lo deseado para la integral de Itô. Las pruebas para el caso de Stratonovich y
Backward-Itô son análogas a esta demostración.
En el Capítulo 6 se discute la forma en que las ecuaciones de Langevin parametrizadas
por el tiempo del laboratorio pueden ser reparametrizadas en términos del tiempo propio de la
partícula y del tiempo de un observador en movimiento. La reparametrización en términos del
tiempo propio es de especial interés, pues este tiempo es un invariante en todos los marcos de
referencia. Una motivación de la ecuación de Langevin (4.2) y una forma de obtener coe�cientes
� y D �realistas�se discute en el Capítulo 7.
4.1.2. Las ecuaciones de Fokker-Planck
Las condiciones para la existencia de la densidad de transición del proceso de Langevin (4.2)
dependen de la forma que tomen las funciones coe�ciente involucradas F i, aij y cij . En esta
tesis no se hace énfasis en estas condiciones pues se construyen modelos para el movimiento
Browniano relativista a partir del supuesto de que el proceso (4.2) tiene densidad estacionaria,
lo cual es su�ciente para el análisis que se hace. Las propiedades que debe tener la ecuación
de Langevin (4.2) para que su densidad de transición (si la tiene) cumpla con la ecuación de
Fokker-Planck son mucho más sencillas, aunque tampoco se hace gran énfasis en ellas. Los
casos en los que la ecuación de Fokker-Planck es necesaria, corresponden todos a ecuaciones de
Langevin autónomas, en las cuales se supone de entrada la existencia de la densidad estacionaria
y solo se necesita utilizar la ecuación de Fokker-Planck estacionaria. En este último caso, dada
una condición inicial como en el teorema 29 para la ecuación (4.2), la validez de la ecuación de
Fokker-Planck queda asegurada debido a la existencia de la densidad estacionaria, la autonomía
de la ecuación de Langevin y los supuestos hechos sobre los coe�cientes de esta ecuación.
En el Capítulo 2 se dieron algunas condiciones para la validez de la ecuación de Fokker-
Planck y la existencia de la densidad de transición en el caso unidimensional. Considerando la
información anterior, es fácil ver la forma que tiene la ecuación de Fokker-Planck.
Proposición 60 Sean V � Rd y M � Rd conjuntos abiertos y conexos. Supongamos lo si-
guiente:
(a) La ecuación de Langevin (4.2) satisface la condición de Lipschitz (2.8) y es autónoma, es
70
decir que los coe�cientes F i, aij y cij no dependen del tiempo t.
(b) La ecuación (4.2) tiene condición inicial (X0;P0) = (x0;p0), siendo (x0;p0) una variable
aleatoria en R2d con distribución P0 e independiente de Wt 8t � 0;
(c) La densidad de transición p (s; z; r;y) del proceso (4.2) existe y está de�nida en [0;1)�
E� [s;1)� E, con E : = V�M;
(d) Dado (s; z) 2 [0;1)� E, la función p (s; z; �; �) es de clase C1 en el conjunto [s;1)� E.
Entonces, la densidad del proceso (Xt;Pt) satisface (según el tipo de integral):
1. En el caso de la integral de Itô �� = ��, la densidad f� (t;x;p) del proceso (Xt;Pt) cumple
la ecuación de Fokker-Planck
�@
@t+pi
p0@
@xi
�f� =
@
@pi
���F i + aijpj
�f� +
1
2
@
@pk
�circkrf�
��, (4.4a)
para (t;x;p) 2 [0;1)� V�M;
2. Cuando se toma la integral de Stratonovich �� = ��, la ecuación de Fokker-Planck para
la densidad f� (t;x;p) del proceso (Xt;Pt) es�@
@t+pi
p0@
@xi
�f� =
@
@pi
���F i + aijpj � 1
2ckr
@
@pkcir�f� +
1
2
@
@pk
�circkrf�
��, (4.4b)
para los puntos (t;x;p) 2 [0;1)� V�M;
3. En el caso de la integral de backward �� = ��, la densidad f� (t;x;p) del proceso (Xt;Pt)
satisface la ecuación de Fokker-Planck
�@
@t+pi
p0@
@xi
�f� =
@
@pi
���F i + aijpj � ckr @
@pkcir�f� +
1
2
@
@pk
�circkrf�
��, (4.4c)
en [0;1)� V�M.
Observación 61 Hay dos cosas que conviene tener en cuenta:
(i) Una ecuación diferencial estocástica autónoma que satisface las condiciones (a) tiene una
única solución global.
71
(ii) Si se reemplaza (a) por las condiciones del Teorema 29, el resultado sigue siendo cierto.
Sin embargo, solo necesitaremos la ecuación de Fokker-Planck para ecuaciones autónomas.
Demostración. Comenzaremos con el caso de la ecuación de Langevin-Itô. Conviene es-
tablecer algo de notación para facilitar el entendimiento de la prueba. De�namos g : [0;1) �
V � M ! Rd y b : [0;1) � V � M ! Mat2d�2d (R) dadas por gi (t;x;p) = pi
(M2+p2)1=2,
gi+d (t;x;p) = F i (t;x;p)� aij (t;x;p) pj , bij (t;x;p) = 0, bi+d;j (t;x;p) = 0, bi;j+d (t;x;p) = 0
y bi+d;j+d (t;x;p) = cij (t;x;p) para i; j 2 f1; : : : ; dg.
Debido a que F i, aij y cij i; j = 1; : : : ; d no dependen del tiempo t, de los Teoremas 37 y
33 se sigue que (Xt;Pt) es un proceso de difusión en [0;1) con deriva g (t;x;p) y matriz de
difusión B (t;x;p) = c (t;x;p) [c (t;x;p)]T , siendo c :=�cij�y (t;x;p) 2 [0;1) � V �M. El
Teorema 37 también nos permite asegurar que las condiciones de la de�nición 106 para (Xt;Pt)
se cumplen uniformemente en [0;1).
Por otro lado, debido a las condiciones impuestas en la ecuación de Langevin y la den-
sidad de transición, se tiene que, para s y z; �jos las derivadas @@yi
�gi (t;y) p (s; z; r;y)
�y
@2
@yi@yj
�Bij (t;y) p (s; z; r;y)
�existen y son continuas para cualesquiera y 2E y t � s.
Tomemos a f (t;x;p) = f� (t;x;p) como la densidad del proceso (Xt;Pt), la cual esta dada
por la Observación 108. Entonces, teniendo en cuenta los dos párrafos anteriores, del Teorema
107 y la Observación 108 se sigue que
0 =@f
@t+
dXi=1
�@gif
@xi+@gi+df
@pi
�
�12
dXi;j=1
"@2�Bijf
�@xi@xj
+@2�Bi;j+df
�@xi@pj
+@�2Bi+d;jf
�@pi@xj
+@2�Bi+d;j+df
�@pi@pj
#
=@f
@t+
dXi=1
"pi
(M2 + p2)1=2@f
@xi+@gi+df
@pi
#� 12
24 dXi;j=1
@2
@xi@xj
2dXr=1
birbjrf
!35+1
2
dXi;j=1
"@2
@xi@pj
2dXr=1
birbj+d;rf
!+
@2
@pi@xj
2dXr=1
bi+d;rbjrf
!+
@2
@pi@pj
2dXr=1
bi+d;rbj+d;rf
!#
=@f
@t+
dXi=1
�pi
p0@f
@xi+@gi+df
@pi
�� 12
dXi;j=1
2dXr=1
@2
@pi@pjbi+d;rbj+d;rf
72
=
@
@t�
dXi=1
pi
p0@
@xi
!f +
dXi=1
0@ @
@pi
24F i � dXj=1
aijpj
35 f1A� 1
2
dXi;j=1
dXr=1
@2
@pi@pjcircjrf
=
@
@t�
dXi=1
pi
p0@
@xi
!f +
dXi=1
@
@pi
0@24F i � dXj=1
aijpj
35 f � 12
dXj;r=1
@
@pjcircjrf
1A ;
donde las funciones están evaluadas en el conjunto [0;1) � E = [0;1) � V �M. La ecuación
anterior es equivalente a
@
@t�
dXi=1
pi
p0@
@xi
!f =
dXi=1
@
@pi
240@�F i + dXj=1
aijpj
1A f +1
2
dXj;r=1
@
@pjcircjrf
35 ;en [0;1)� V�M. Usando el criterio de sumación de Einstein, llegamos a que
�@
@t� pi
p0@
@xi
�f =
@
@pi
���F i + aijpj
�f +
1
2
@
@pkcikcjkf
�;
en [0;1)�V�M, lo que prueba lo deseado en el caso de la integral de Itô. Para las pruebas en
el caso de Stratonovich y backward, conviene encontrar la ecuación de Langevin-Itô equivalente
y luego proceder de una manera similar.
Observación 62 La hipótesis sobre la condición C1 de la densidad de transición, puede reem-
plazarse al solo pedir que, para s y z �jos, las derivadas @@yi
�gi (t;y) p (s; z; r;y)
�y
@2
@yi@yj
�Bij (t;y) p (s; z; r;y)
�existan y sean continuas para cualesquiera y 2E y t � s:
De manera similar al caso no-relativista, la cantidad f� (t;x;p) dtdx1 : : :dxddp1 : : :dpd se
interpreta como la probabilidad de encontrar la partícula Browniana al tiempo t de � en el
conjunto in�nitesimaldYi=1
�xi; xi + dxi
��
dYi=1
�pi; pi + dpi
�, donde x =
�xi�y p =
�pi�se corres-
ponden con la posición y el momento medidos en �.
Las condiciones iniciales deterministas X0 = x0 y P0 = p0 se convierten en la condición de
frontera �localizada�
f� (0;x;p) = � (x� x0) � (p� p0) 8 (x;p) 2 V�M.
Debido a la linealidad de las ecuaciones (4.4), es posible obtener soluciones más generales para
73
estas ecuaciones al considerar una distribución inicial f0 (x0;p0) no localizada. Sin embargo,
debe tenerse en cuenta que, en un marco relativista, es muy complicado determinar información
inicial �-isócrona por medio de mediciones experimentales.
4.2. Movimientos libres en un baño isotrópico y las relaciones
de Einstein
De manera similar al caso no relativista, algunas restricciones físicas (y, por tanto, matemáti-
cas) sobre las funciones aij y cij surgen al considerar propiedades de simetría en el baño térmico
y del requisito de que las ecuaciones (4.2) deben reproducir la distribución estacionaria y el com-
portamiento de relajación correctos. Por ejemplo, en ausencia de campos de fuerza externos y
en caso de un baño térmico estacionario, isotrópico y homogéneo en el marco �, las funciones
coe�ciente de la ecuación (4.2) toman la forma
F i = 0, aij (t;x;p) = � (p) �ij , cij (t;x;p) = [2D (p)]1=2 �ij
donde � : Rd ! R y D : Rd ! (0;1) son C1; y �ij es una delta de Kronecker. La esta-
cionariedad y la homogeneidad implican que las funciones aij y cij solo dependen del momento
(y no del tiempo y la posición); mientras que la isotropía causa la forma diagonal de cij (por
aquello de que se presentan las mismas propiedades, independientemente de la dirección en que
se midan).
Bajo los supuestos establecidos (baño estacionario, isotrópico y homogéneo en �), las fun-
ciones � y D dependen sólo del momento de las partícula Browniana. Sin embargo, es posible
encontrar funciones de ruido y fricción que dependan de la energía relativista de la partícula
Browniana. Defínanse b� : [M;1) ! Rd y bD : [M;1) ! (0;1) dadas por b� �p0� = � (p) ybD �p0� = D (p). Así, en este caso, la ecuación de Langevin relativista (4.2) se simpli�ca de la
siguiente forma
dX it =
�P itP 0t
�dt, (4.5a)
dP it = �b� �P 0t �P it dt+ h2 bD �P 0t �i1=2 � dW it , (4.5b)
74
mientras que las ecuaciones de Fokker-Planck (4.4) toman la forma
�@
@t+pi
p0@
@xi
�f� =
@
@pi
�b�pif� + @
@pi
� bDf��� , (4.6a)�@
@t+pi
p0@
@xi
�f� =
@
@pi
�b�pif� + bD 12@
@pi
� bD 12 f���, (4.6b)�
@
@t+pi
p0@
@xi
�f� =
@
@pi
�b�pif� + bD @
@pif�
�. (4.6c)
Conviene notar que, debido a que � y D son funciones C1, las funciones b� y bD son C1 en
(M;1).
Una restricción adicional sobre las funciones � (p) = b� �p0� y D (p) = bD �p0� surge deconsideraciones termoestadísticas: si el movimiento de la partícula Browniana relativista está
restringido a un volumen �nito V � Rd (abierto y conexo) y el baño térmico esta en equilibrio
térmico a temperatura T = (kB�)�1, entonces el proceso (Xt;Pt) debe tener como densidad
estacionaria f1 (x;p) a una función de Jüttner espacialmente homogénea1
f1 (x;p) = N exph���M2 + p2
� 12
iJ (x;V) , (4.7)
siendo J (x;V) la función indicadora del volumen accesible V y N una constante de norma-
lización. Insertando la ecuación (4.7) dentro de las ecuaciones de Fokker-Planck (4.4), se encuen-
tra que, dependiendo de la integral considerada, las funciones � (p) = b� �p0� y D (p) = bD �p0�satisfacen la siguientes relaciones de �uctuación-disipación.
Teorema 63 (Relaciones relativistas de Einstein) Sea f : Rd�Rd ! [0;1) una función
de densidad de Jüttner espacialmente homogénea en V � Rd, dada por
f (x;p) = N exph���M2 + p2
�1=2iJ (x;V) .Si la densidad anterior satisface las ecuaciones de Fokker-Planck (4.6), entonces los coe�cientesb� y bD se cumplen las siguientes relaciones:
1Cuando decimos �debe tener�, nos referimos a que es una condición para que la ecuación (4.2) exhiba elcomportamiento de equilibrio térmico.
75
(i) Para la integral de Itô �� = ��,
0 = b� �p0� p0 � � bD �p0�+ bD0 �p0� 8p0 2 (M;1) ; (4.8a)
(ii) En el caso de la integral de Stratonovich, �� = ��,
0 = b� �p0� p0 � � bD �p0�+ bD0 �p0�2
8p0 2 (M;1) ; (4.8b)
(iii) Para la integral Backward-Itô �� = ��,
0 = b� �p0� p0 � � bD �p0� 8p0 2 (M;1) . (4.8c)
Observación 64 Las relaciones relativistas de Einstein también se conocen como relaciones
de �uctuación-disipación relativistas.
Demostración. Empezaremos con el caso de Itô. De la ecuación (4.6a) se sigue que
b� �p0� pif (x;p) + @
@pi
� bD �p0� f (x;p)� = c 8 (x;p) 2 V� Rd;
donde c 2 R es una constante. Derivando en la última ecuación, se encuentra que c = 0 y que
0 = b� �p0� pif (x;p) + @
@pi
� bD �p0� f (x;p)� == b� �p0� pif (x;p) + pi
p0@ bD �p0�@p0
f (x;p)� � pi
p0bD �p0� f (x;p) 8 (x;p) 2 V� Rd:
De lo anterior se obtiene que
0 = b� �p0� p0 + @ bD �p0�@p0
� � bD �p0� 8p 2 Rd � f0g ;
o lo que es equivalente
0 = b� �p0� p0 + bD0 �p0�� � bD �p0� 8p0 2 (M;1) :
76
Esta última igualdad prueba lo deseado en el caso de Itô. En los otros dos casos se procede
análogamente.
Comparando las relaciones relativistas de Einstein (4.8) con la relación no-relativista (3.9),
puede observarse que en el caso relativista la masa ha sido reemplazada por la energía p0. Las
relaciones relativistas de Einstein proporcionan condiciones que deben satisfacer los coe�cientes
de ruido y fricción para que la ecuación de Langevin modele adecuadamente la dinámica es-
tocástica de una partícula Browniana que se encuentra en ausencia de fuerzas externas e inmersa
en un baño termalizado (equilibrio térmico), estacionario, isotrópico y homogéneo en �:
77
Capítulo 5
Ejemplos Unidimensionales y el
Desplazamiento Cuadrático Medio
En este Capítulo se consideran algunos ejemplos de modelos del movimiento Browniano re-
lativista. Presentamos los detalles y cálculos del caso unidimensional (d = 1). La generalización
a dimensiones mayores es directa aunque los cálculos pueden ser más tediosos. Se comienza
estableciendo las ecuaciones diferenciales estocásticas para el proceso energía P 0t y el proceso
velocidad Vt := PtP 0t
1. Se discuten resultados para la constante de difusión asintótica de algunos
ejemplos de movimiento Browniano tratados.
5.1. Consideraciones generales y las ecuaciones de la energía y
la velocidad
En este Capítulo, se considera que la partícula Browniana está inmersa en un baño térmico
estacionario, isotrópico y homogéneo, que tiene d = 1 dimensiones espaciales. También se supone
que el baño térmico está en estado de equilibrio a temperatura T y en ausencia de fuerzas
1Debido a que se trabajara en dimensión uno, utilizaremos Xt, Pt y Vt, en lugar de Xt, Pt y Vt, para denotarla posición, momento y velocidad, respectivamente, de la partícula Browniana relativista.
79
externas. Bajo estas condiciones, la ecuación de Langevin del proceso (Xt; Pt) toma la forma
dXt =PtP 0t
dt; (5.1a)
dPt = ��� (Pt)Ptdt+ [2D (Pt)]1=2 � dWt: (5.1b)
Dado que se quiere que la ecuación anterior genere procesos estocásticos equivalentes para
cada una de las integrales, se pide que
�� (p) = �� (p)�D0 (p)
2p; (5.2a)
�� (p) = �� (p)�D0 (p)
p; (5.2b)
para todo p 2 R�f0g. Estas igualdades son una consecuencia directa de las condiciones de
equivalencia entre ecuaciones diferenciales estocásticas (46) y (50) que se describieron en el
Capítulo 2. Hay otras formas de expresar la equivalencia de las ecuaciones diferenciales, sin
embargo se elige esta por conveniencia, como más adelante se verá.
Proposición 65 Supongamos que la ecuación (5.1) se cumple. Entonces, bajo las relaciones
(5.2), los siguientes procesos son equivalentes:
dPt = ��� (Pt)Ptdt+ [2D (Pt)]1=2 � dWt
dPt = ��� (Pt)Ptdt+ [2D (Pt)]1=2 � dWt
dPt = ��� (Pt)Ptdt+ [2D (Pt)]1=2 � dWt
Además, como el baño térmico está termalizado (en equilibrio térmico) a temperatura T ,
tenemos que la densidad estacionaria de momentos del proceso (5.1) debe ser una función de
Jüttner unidimensional de la forma
�1 (p) = N exp����M2 + p2
�1=2� 8p 2 R; (5.3)
siendo N una constante de normalización y T = (kB�)�1. Luego, los coe�cientes de ruido y
80
fricción de la ecuación (5.1) deben satisfacer las relaciones de Einstein (4.8), las cuales pueden
escribirse de manera alternativa como en la proposición que sigue.
Proposición 66 Supongamos la ecuación (5.1) se cumple y que su densidad estacionaria de
momentos � es una función de Jüttner unidimensional de la forma (5.3). Entonces,
0 � �� (p) p0 � �D (p) (5.4a)
0 � �� (p) p0 � �D (p) +D0 (p)
p0
2p(5.4b)
0 � �� (p) p0 � �D (p) +D0 (p)
p0
2p(5.4c)
para cualquier p 2 R�f0g.
Demostración. Comenzaremos demostrando que los coe�cientes ��, ��, �� y D satis-
facen las relaciones de Einstein (4.8). Posteriormente, utilizaremos este hecho para demostrar
lo deseado.
Tenemos que la densidad estacionaria del proceso (5.1) cumple, según la integral que se tome,
ecuaciones de Fokker-Planck análogas a (4.6). Entonces, la densidad estacionaria de momentos
� de la ecuación de Langevin (5.1) satisface, para todo p 2 R, las igualdades
0 =@
@p
�b�� �p0� p� (p) + @
@p
� bD �p0�� (p)�� ,0 =
@
@p
�b�� �p0� p� (p) + bD 12�p0� @@p
� bD 12�p0�� (p)
��,
0 =@
@p
�b�� �p0� p� (p) + bD �p0� @@p� (p)
�,
siendo bD �p0� = D (p) y b�� �p0� = �� (p) ; 8p 2 R y 8� 2 f�; �; �g. Esto implica que (ver la
demostración de 63)
0 = b�� �p0� p� (p) + @
@p
� bD �p0�� (p)�= b�� �p0� p� (p) + bD0 �p0� p
p0� (p)� � p
p0� (p) bD �p0� ,
0 = b�� �p0� p� (p) + bD 12�p0� @@p
� bD 12�p0�� (p)
�= b�� �p0� p� (p) + 1
2bD0 �p0� p
p0� (p)� � p
p0� (p) bD �p0� ,
81
0 = b�� �p0� p� (p) + bD �p0� @@p� (p) = b�� �p0� p� (p)� � p
p0� (p) bD �p0� ,
para todo p 2 R. Luego, por ser � (p) > 0 8p 2 R, se tiene
0 = b�� �p0� p0 � � bD �p0�+ bD0 �p0� ;0 = b�� �p0� p0 � � bD �p0�+ bD0 �p0� =2;0 = b�� �p0� p0 � � bD �p0� ;
para todo p 2 R�f0g, lo cual demuestra que los coe�cientes ��, ��, �� y D satisfacen las
relaciones de Einstein (4.8).
Por otro lado, debido a que D (p) = bD �p0 (p)� 8p 2 R, se cumple queD0 (p)
p0
p= bD0 �p0� 8p 2 R�f0g :
Entonces, sustituyendo la igualdad de arriba en la ecuación (5.5), llegamos a que
0 � �� (p) p0 � �D (p)
0 � �� (p) p0 � �D (p) +D0 (p)
p0
2p
0 � �� (p) p0 � �D (p) +D0 (p)
p0
2p
para todo p 2 R�f0g, lo que prueba lo deseado.
Ahora procederemos a establecer las ecuaciones de la energía y la velocidad. TomandobD �p0� = D (p) y b�� �p0� = �� (p) 8p 2 R y 8� 2 f�; �; �g, de la ecuación (4.3) se obtiene el
siguiente resultado para el proceso de energía P 0t :=�M2 + P 2
�1=2.Corolario 67 Supongamos que la ecuación (5.1) se cumple. Entonces, el proceso P 0t :=
�M2 + P 2
�1=2tiene asociada la ecuación diferencial estocástica
dP 0t =
(�b�� �P 0t �P 0t
"1�
�M
P 0t
�2#+ ��
bD �P 0t �P 0t
�M
P 0t
�2)dt
(2 bD �P 0t �
"1�
�M
P 0t
�2#)1=2� dWt;
82
donde �� = 1, �� = 0 y �� = �1.
Resta encontrar la ecuación de la velocidad, para lo cual se de�ne eD; e�� : (�1; 1)! R para
� 2 f�; �; �g, dadas por
e�� (V ) := �� (P (V )) ; eD (V ) := D (P (V ))
para P (V ) = MV�1� V 2
��1=2 (la fórmula del momento relativista). Conviene notar que porser �� y D de clase C1, también b�� y D lo son.
Las ecuaciones diferenciales estocásticas de la velocidad están dadas por el siguiente resul-
tado, el cuál utilizaremos repetidamente en los ejemplos.
Proposición 68 Supongamos que la ecuación (5.1) se cumple. Entonces, el proceso estocástico
Vt :=PtP 0tsatisface
dVt =
"�e�� (Vt) �1� V 2t �� ��
3 eD (Vt)M2
!�1� V 2t
�2#Vtdt
+
" 2 eD (Vt)M2
!�1� V 2t
�3#1=2 � dWt;
siendo �� = 1, �� = 0 y �� = �1.
Demostración. Para esta demostración utilizaremos las fórmulas de Itô (26), Stratonovich
(47) y backward-Itô (51). De�namos u : R� [0;1)! R dada por u (p; t) = p
(M2+p2)1=2. Nótese
que @u@p y
@u@t existen y
@u
@p(p; t) =
M2
(M2 + p2)3=2,
@u
@t(p; t) = 0;
para todo (p; t) 2 R� [0;1).
Por otro lado, de lo anterior se tiene que @2u@p2
existe y
@2u
@p2(p; t) =
�3M2p
(M2 + p2)5=28 (p; t) 2 R� [0;1) :
83
De la de�nición de u y lo hecho arriba, es claro que u es continua y que @u@p ,@u@t ,
@2u@p2
existen
y son continuas. Luego, podemos aplicar las fórmulas de Itô, Stratonovich y backward-Itô al
proceso Vt := Pt
(M2+P 2t )1=2 = u (Pt; t). Mediante la fórmula de Itô, para la ecuación (5.1) con la
integral de Itô obtenemos
dVt =@u
@t(Pt; t) dt+
@u
@p(Pt; t) dPt +
1
2[2D (Pt)]
@2u
@p2(Pt; t) dt
=M2�
M2 + P 2t�3=2dPt � 3M2Pt�
M2 + P 2t�5=2D (Pt) dt
=
�1� V 2t
��M2 + P 2t
�1=2dPt � �3D (Pt)M2
��1� V 2t
�2Vtdt
= ��1� V 2t
��M2 + P 2t
�1=2�� (Pt)Ptdt+�1� V 2t
��M2 + P 2t
�1=2 [2D (Pt)]1=2 � dWt ��3D (Pt)
M2
��1� V 2t
�2Vtdt
=
���� (Pt)
�1� V 2t
���3D (Pt)
M2
��1� V 2t
�2�Vtdt+
�1� V 2t
�3=2M
[2D (Pt)]1=2 � dWt
=
���� (Pt)
�1� V 2t
���3D (Pt)
M2
��1� V 2t
�2�Vtdt+
��2D (Pt)
M2
��1� V 2t
�3�1=2 � dWt
=
"�e�� (Vt) �1� V 2t ��
3 eD (Vt)M2
!�1� V 2t
�2#Vtdt+
" 2 eD (Vt)M2
!�1� V 2t
�3#1=2 � dWt.
Además, considerando el desarrollo anterior tenemos que, para la ecuación (5.1) con las
integrales de Stratonovich y backward-Itô, se cumple
dVt =@u
@t(Pt; t) dt+
@u
@p(Pt; t) dPt
= �e�� (Vt) �1� V 2t �Vtdt+" 2 eD (Vt)M2
!�1� V 2t
�3#1=2 � dWt; y
dVt =@u
@t(Pt; t) dt+
@u
@p(Pt; t) dPt �
1
2[2D (Pt)]
@2u
@p2(Pt; t) dt
=
"�e�� (Vt) �1� V 2t �+
3 eD (Vt)M2
!�1� V 2t
�2#Vtdt+
" 2 eD (Vt)M2
!�1� V 2t
�3#1=2 � dWt.
84
Entonces, resulta que
dVt =
"�e�� (Vt) �1� V 2t �� ��
3 eD (Vt)M2
!�1� V 2t
�2#Vtdt
+
" 2 eD (Vt)M2
!�1� V 2t
�3#1=2 � dWt,
para �� = 1, �� = 0 y �� = �1; lo cual demuestra lo deseado.
Como se ve en la Sección 5.3, la ecuación de la proposición anterior se usará para calcular
el desplazamiento cuadrático medio.
5.2. Desplazamiento cuadrático medio asintótico
El objetivo de este apartado es tratar la constante de difusión efectiva asociada a modelos
especí�cos de movimiento Browniano relativista con el objetivo de utilizar tal información en el
análisis de los ejemplos que se proponen en la siguiente sección. Para mayor claridad, conviene
retomar la de�nición de constante de difusión asintótica.
De�nición 69 Sea fYt; t � 0g un proceso estocástico relativista. La constante de difusión asin-
tótica de Y es
D1 (Y ) = l��mt!1
D(Yt � Y0)2
E2t
:
Aquí no se demostrará que la constante de difusión asintótica esta bien de�nida, una posible
demostración puede encontrarse en el artículo de Angst y Franchi [5]. Determinar la constante de
difusión asintótica es un objetivo básico en cualquier teoría y modelo del movimiento Browniano
pues permite estimar los procesos de relajación de la partícula Browniana y aplicar resultados
de la teoría de difusiones.
Físicamente, la constante de difusión asintótica se interpreta como una constante de difusión
efectiva. Es decir que, para el caso de la partícula Browniana considerada, D1 (X) es una
medida de la �facilidad�(real) con la que se mueve la partícula Browniana.
Dada la importancia de la constante de difusión asintótica y lo poco práctico de su de�nición,
es deseable encontrar una ecuación que permita calcular dicha constante de manera sencilla, al
85
menos para los ejemplos que se consideran en este trabajo. Un resultado que cumple con esta
meta fue propuesto por B. Lindner (ver págs. 69 y 70 del artículo [26] de Dunkel y Hänggi).
Teorema 70 (Lindner) Sea fYt; t � 0g un proceso estocástico relativista cuya velocidad Vt :=dYtdt es un proceso estacionario y cumple
dVt = �a� (Vt)Vt + [2b (Vt)]1=2 � dWt; (5.6)
con coe�cientes a�; b : (�1; 1)! (0;1) diferenciables y tales que a� (�w) = a� (w) y b (�w) =
b (w). Además, defínase U; �� : [0; 1)! R dadas por
��(w) = a� (w)w � b0 (w) = a� (w)w,
U (w) =
Z w
0
��(s)
b (s)ds:
Entonces,
D1 (Y ) '
R 10 e
U(y)hR 1yxe�U(x)
b(x) dxi2dyR 1
0e�U(z)b(z) dz
: (5.7)
Observación 71 Para el baño térmico supuesto en este capítulo, cuando la ecuación de Langevin
tiene solución, ésta es un proceso homogéneo. Por otra parte, en los modelos de movimiento
Browniano relativista que se tratan aquí, se supone de entrada la existencia de una densidad
estacionaria para el proceso de Langevin. Por lo tanto, según la Observación 105, el momen-
to y la velocidad son proceso estacionarios en los ejemplos de este capítulo; en particular, la
velocidad cumple una de las hipótesis del Teorema anterior.
Ahora se tiene todo lo necesario para un análisis completo de algunos ejemplos.
5.3. Ejemplos
Se consideraran tres ejemplos unidimensionales de modelos relativistas de Langevin, cuyas
densidades estacionarias del momento se distribuyen como funciones de Jüttner �J (p) /
exp[��p2 +M2
�1=2], con T =(kB�)�1 la temperatura del baño. En estos casos, la relación
86
de Einstein (5.4) implica que, en la ecuación (5.1), solo una de las funciones �� (p) ó D (p)
puede ser elegida arbitrariamente.
5.3.1. Amplitud constante
Como primer ejemplo consideraremos el llamado �Proceso Relativista de Ornstein-Uhlenbeck�
(ROUP). Este proceso está de�nido por la elección
�� (p) = �cM
p08p 2 R;
en la ecuación (5.1), donde �c > 0 es un parámetro de fricción constante. A partir de las
relaciones de Einstein (5.4) se encuentra2
D (p) � �� (p) p0��1 = �cM��1 =: Dc 8p 2 R;
por lo cual es claro que el ROUP corresponde al caso de amplitud constante. Por lo tanto, la
ecuación de Langevin asociada al ROUP queda como
dPt = ��cM
p0Ptdt+
�2�cM
�
�1=2� dWt (5.8a)
= ��cM
p0Ptdt+
�2�cM
�
�1=2� dWt: (5.8b)
En este caso, debido a que la amplitud del ruido Dc := �cM=� no depende del momento P ,
los coe�cientes en la ecuación de Langevin no varían según el tipo de integral estocástica. Sin
embargo, es importante tener en cuenta la clase de integral estocástica para poder escribir la
ecuación diferencial estocástica asociada al proceso de la velocidad Vt := PtP 0t.
La existencia y unicidad de la solución de la ecuación (5.8) es una consecuencia inmediata
del Corolario 33.
Proposición 72 Supongamos que (X0; P0) es una variable aleatoria independiente de Wt para
t � 0. Entonces, la ecuación (5.8) tiene una única solución global.
2Se esta suponiendo que la partícula Browniana esta en equilibrio térmico con el baño circundante, por lotanto los coe�cientes deben cumplir las relaciones de Einstein (5.4).
87
A continuación se encuentra la ecuación diferencial estocástica asociada a la velocidad Vt.
Proposición 73 Sea (Xt; Pt) como en (5.8). Entonces el proceso Vt := PtP 0ttiene la diferencial
dVt = ��c��1� V 2t
�3=2 � 3
�
�1� V 2t
�2�Vtdt+
�2�c�
�1� V 2t
�3�1=2 � dWt (5.9a)
= ��c��1� V 2t
�3=2+3
�
�1� V 2t
�2�Vtdt+
�2�c�
�1� V 2t
�3�1=2 � dWt; (5.9b)
donde � = �M .
Observación 74 Considerando la relación Vt := Pt(M2+P 2t )
, P 2t :=M2V 2t1�V 2t
, es claro que dada
una solución de (5.9) se tiene una solución de (5.8), y viceversa.
Demostración. Basta aplicar la Proposición 68. Observemos que las ecuaciones (5.8) son
ecuaciones de Langevin del tipo (5.1) con
�� (p) = �cM
p0= �c
�1� v2
�1=2= e�� (v) ;
�� (p) = �cM
p0= �c
�1� v2
�1=2= e�� (v) ;
D (p) =�cM
�=�cM
2
�= eD (v) ;
bajo la identidad v = pp0. Aplicando la Proposición 68 a las ecuaciones (5.8) se llega al resultado
deseado.
Ahora estamos en condiciones de calcular el coe�ciente de difusión asintótica D1 (X). Para
esto conviene tener en cuenta el siguiente resultado auxiliar sobre las funciones de Bessel. Para
mayores detalles sobre estas funciones se sugiere [57].
Lema 75 La función de Bessel modi�cada de segunda especie de orden 0, denotada por K0,
satisface la relación
(�1)� d�
d��K0 (�) =
Z 1
0exp
� �
(1� v2)1=2
!�1� v2
��(�+2)=2dv:
88
Proposición 76 La posición Xt asociada al ROUP (5.8), satisface
DROUP1 (�c) � D1 (X) = (�c�)
�1 =kBTM�c
;
con � = �M .
Observación 77 Es interesante notar que con el ROUP se obtiene la misma constante de di-
fusión asintótica que para el proceso de Ornstein-Uhlenbeck no-relativista, DROUP1 = (�c�)
�1 =
DOU1 ; para cualesquiera parámetros (�c; T ;M). Esta puede ser considerada la razón por la cual
el ROUP recibe tal nombre.
Demostración. (Proposición 76). Para encontrar lo pedido, utilizaremos el Teorema 70.
Notemos que la velocidad del ROUP (5.8) cumple la ecuación (5.9), la cual es de la forma (5.6)
con
a� (v) = �c
��1� v2
�3=2 � 3
�
�1� v2
�2�;
b (v) =�c�
�1� v2
�3;
para todo v 2 (�1; 1). Además, según la notación del teorema de Lindner,
��(v) = �c
��1� v2
�3=2+3
�
�1� v2
�2�v;
U (v) =
Z v
0
�c
h�1� w2
�3=2+ 3
�
�1� w2
�2iw
�c� (1� w2)
3 dw
= �
Z v
0
w
(1� w2)3=2dw + 3
Z v
0
w
(1� w2)dw
=�
(1� v2)1=2� �� 3
2log�1� v2
�;
para v 2 [0; 1).
Luego, aplicando la Observación 71, la ecuación (5.7) y el Lema 75, vemos que
D1 (X) =e��
�c
R 10 e
U(y)nR 1
y x exph���1� x2
��1=2i �1� x2
��3=2dxo2dyR 1
0 exph�� (1� z2)�1=2
i(1� z2)�3=2 dz
89
= � e��
�cdd�K0 (�)
Z 1
0eU(y)
(� 1�exp
h���1� x2
��1=2i����x=1x=y
)2dy
= � 1
�c�dd�K0 (�)
Z 1
0exp
h��1� y2
��1=2i �1� y2
��3=2exp
h�2�
�1� y2
��1=2idy
= � 1
�c�dd�K0 (�)
Z 1
0exp
h���1� y2
��1=2i �1� y2
��3=2dy
=
dd�K0 (�)
�c�dd�K0 (�)
= (�c�)�1 ;
lo que demuestra lo deseado.
5.3.2. Coe�ciente de fricción constante en la ecuación de Langevin-backward
El siguiente ejemplo de modelo de movimiento Browniano relativista que se considera cor-
responde al caso especial de un coe�ciente de fricción constante �� (p) � �y > 0 en la ecuación
de Langevin (5.1) con la integral backward. Abreviaremos este ejemplo de proceso como RBM.
En este caso, la relación de Einstein (5.4) implica
D(p) = �yp0��1 = �y�
�1 �M2 + p2�1=2 8p 2 R.
Entonces la ecuación de Langevin asociada a un proceso RBM queda como
dPt = ��yPtdt+�2�yP
0t
�
�1=2� dWt (5.10a)
= ��y��P 0t � 1�P 0t
�Ptdt+
�2�yP
0t
�
�1=2� dWt: (5.10b)
A continuación se demuestra que la ecuación anterior tiene solución única, aunque para esto
conviene tener en cuenta el siguiente resultado.
Lema 78 Sea g : R! R+ una función dada por g (P ) =�M2 + P 2
�1=2=: P 0. Entonces, se
cumple que ��P 0 �Q0�� = jg (P )� g (Q)j � jP �Qj 8 (P;Q) 2 R2:
90
Demostración. Para todo (P;Q) 2 R2 � f0g, sucede que
��P 0 �Q0�� = jg (P )� g (Q)j
=����M2 + P 2
�1=2 � �M2 +Q2�1=2���
=
��P 2 �Q2�����(M2 + P 2)1=2 + (M2 +Q2)1=2���
���P 2 �Q2��jjP j+ jQjj
� jjP j � jQjj
� jP �Qj ;
lo que prueba lo deseado.
Proposición 79 Sea (X0; P0) una variable aleatoria independiente deWt para t � 0. Entonces,
la ecuación (5.10) tiene una única solución global.
Demostración. Por el Corolario 33, es claro que basta probar que los coe�cientes de ruido
y fricción satisfacen una condición de Lipschitz de la forma (2.8).
De�namos f : R! R dada por f (P ) = PP 0. Claramente, se cumple que f 0 (P ) = M2
(M2+P 2)3=2�
1M 8P 2 R. Aplicando el Teorema del Valor Medio a la función anterior, se llega a que
�����y��P 0 � 1�P 0
�P � �y
��Q0 � 1�Q0
�Q
���� = �y
����P � P
�P 0�Q+ Q
�Q0
����� �y
�jP �Qj+ 1
�
���� PP 0 � Q
Q0
������ �y
�jP �Qj+ 1
�
����m�axx2Rf 0 (x)
���� jP �Qj�=
��y +
1
�M
�jP �Qj (5.11)
para cualquier (P;Q) 2 R2. Por otro lado, del Lema anterior se sigue que
������2�yP
0
�
�1=2��2�yQ
0
�
�1=2����� =�2�y�
�1=2 ����P 0�1=2 � �Q0�1=2���
91
=
�2�y�
�1=2 ��P 0 �Q0�����(P 0)1=2 + (Q0)1=2����
�2�y�
�1=2 jP �Qj���(M2 + P 2)1=4 + (M2 +Q2)1=4���
��2�y�
�1=2 jP �Qj��2M1=2��
=
��y2M�
�1=2jP �Qj (5.12)
para todo (P;Q) 2 R2.
De las ecuaciones (5.11) y (5.12) se sigue que la ecuación (5.10b) satisface una condición de
Lipschitz de la forma (2.8), lo cual prueba lo deseado.
La ecuación diferencial estocástica asociada a la velocidad permite calcular el coe�ciente de
difusión asintótica de la siguiente manera.
Proposición 80 Sea (Xt; Pt) como en (5.10). Entonces, el proceso Vt := PtP 0ttiene la diferencial
estocástica
dVt = ��y��1� V 2t
�� 3
�
�1� V 2t
�3=2�Vtdt+
�2�y�
�1� V 2t
�5=2�1=2 � dWt (5.13a)
= ��y��1� V 2t
�+2
�
�1� V 2t
�3=2�Vtdt+
�2�y�
�1� V 2t
�5=2�1=2 � dWt; (5.13b)
tomando � = �M .
Demostración. Haremos uso de la Proposición 68. Notemos que las ecuaciones (5.10) son
ecuaciones de Langevin del tipo (5.1) con
�� (p) = �y = e�� (v) ;�� (p) = �y
��p0 � 1�p0
�= �y
�1� 1
�
�1� v2
�1=2�= e�� (v) ;
D (p) =�yp
0
�=�yM
2
�
�1� v2
��1=2= eD (v) ;
y v = pp0. Aplicando la Proposición 68 a las ecuaciones (5.10) se llega al resultado deseado.
Para el coe�ciente de difusión asintótica correspondiente se obtiene lo siguiente.
92
Proposición 81 La posición Xt asociada al proceso RBM (5.10) tiene como constante de
difusión asintótica a
DRBM1 (�y) � D1 (X) = (�y�)
�1 K0 (�)
K1 (�); (5.14)
donde � = �M y Kn denota la función de Bessel modi�cada de segunda especie de orden n 2 N.
Demostración. Similar a la demostración de 76, utilizaremos el Teorema 70. Notemos que
la velocidad del proceso Xt satisface la ecuación (5.13), la cual es de la forma (5.6) con
a� (v) = �y
��1� v2
�� 3
�
�1� v2
�3=2�;
b (v) =�y�
�1� v2
�5=2;
para cualquier v 2 (�1; 1). Entonces, acorde a la notación del Teorema de Lindner 70,
��(v) = �y
��1� v2
�+2
�
�1� v2
�3=2�v;
U (v) =
Z v
0
�yh�1� w2
�+ 2
�
�1� w2
�3=2iw
�y� (1� w2)
5=2dw
= �
Z v
0
w
(1� w2)3=2dw + 2
Z v
0
w
(1� w2)dw
=�
(1� v2)1=2� �� log
�1� v2
�;
para v 2 [0; 1).
Aplicando el Teorema de Lindner 70, la Observación 71 y el lema 75 sobre la función de
Bessel, y considerando la demostración de la Proposición 76, tenemos que
DRBM1 (�y) � D1 (X) =
e��
�y
R 10 e
U(y)nR 1
y x exph���1� x2
��1=2i �1� x2
��3=2dxo2dyR 1
0 exph�� (1� z2)�1=2
i(1� z2)�3=2 dz
= � e�
�y�dd�K0 (�)
Z 1
0eU(y) exp
h�2�
�1� y2
��1=2idy
= � 1
�y�dd�K0 (�)
Z 1
0exp
h��1� y2
��1=2i �1� y2
��1exp
h�2�
�1� y2
��1=2idy
93
= � 1
�y�dd�K0 (�)
Z 1
0exp
h���1� y2
��1=2i �1� y2
��1dy
= � K0 (�)
�y�dd�K0 (�)
:
Luego, como la función K0 (�) cumple dd�K0 (�) = �K1 (�), de lo anterior se sigue
D1 (X) = (�y�)�1 K0 (�)
K1 (�);
con lo que se prueba lo deseado.
A bajas temperaturas � := (kBT )�1 !1 y la relación (5.14) se reduce al resultado clásico
DRBM1 � kBT
M�y. Mientras que para temperaturas muy altas, i.e. para �M � 1, se encuentra
DRBM1 = (�yM)
�1�� " + log
�2
�+O
h(�M)2
i��, (5.15)
donde " � 0;577216 es la constante de Euler. Debe tenerse en mente que, debido a los pro-
cesos de aniquilación/creación de partículas a altas energías, las teorías clásicas no-cuánticas
se vuelven inválidas a muy altas temperaturas (en el límite � = �M � 1), y, por lo tanto, la
igualdad (5.15) tiene un uso práctico limitado.
Comparando el ejemplo de esta subsección con el proceso de Ornstein-Uhlenbeck relativista,
se observa que DRBM1 (�y) � DROUP
1 (�c) para �c = �y. Intuitivamente, esto puede explicarse
debido a que para el proceso ROUP (5.8), el valor absoluto de la fuerza de fricción esta acotado
por �cM , mientras que la fuerza de fricción no esta acotada para el modelo RBM (5.10). De
este modo, la difusión espacial es más fuerte en el último caso.
5.3.3. Coe�ciente de fricción constante en la ecuación de Langevin-Itô
El proceso RBM tratado en el apartado anterior está caracterizado por un coe�ciente de
fricción constante �y cuando se adopta la integral backward en la ecuación de Langevin (5.1).
Otro modelo de movimiento Browniano relativista, que denotaremos por RBM(I), se obtiene al
considerar un coe�ciente de fricción constante �� (p) � �� > 0 en la ecuación de Langevin (5.1)
con la integral de Itô.
94
La relación relativista de Einstein (5.4) implica que
D (p) =��
�2�1 + �p0
�8p 2 R:
Así que la ecuación de Langevin asociada al RBM(I) es
dPt = ���Ptdt+�2��
�2�1 + �P 0t
��1=2� dWt; (5.16)
la cual tiene solución única.
Proposición 82 Supongamos que (X0; P0) es una variable aleatoria independiente de Wt para
t � 0. Entonces, la ecuación (5.16) tiene una única solución global.
Demostración. Procederemos como en la demostración de la existencia y unicidad de la
solución del proceso RBM. Del Lema 78 se sigue que, para cada (P;Q) 2 R2,
j��P � ��Qj+������2��
�2�1 + �P 0
��1=2��2��
�2�1 + �Q0
��1=2������ �� jP �Qj+
�2��
�2
�1=2 ����1 + �P 0�1=2 � �1 + �Q0�1=2���� �� jP �Qj+
�2���
�1=2 ��P 0 �Q0�����(1 + �P 0)1=2 + (1 + �Q0)1=2���� �� jP �Qj+
���2�
�1=2jP �Qj
=
�� +
���2�
�1=2!jP �Qj ;
lo que prueba lo deseado.
Por otro lado, la velocidad del RBM(I) dado por (5.16) queda determinada de la siguiente
manera.
Proposición 83 Sea (Xt; Pt) como en la relación (5.16). Entonces, la ecuación diferencial
95
asociada al proceso Vt = PtP 0tes
dVt = �����1� V 2t
�+3
�
�1� V 2t
�3=2+3
�2�1� V 2t
�2�Vtdt+ (5.17)�
2���
h�1� V 2t
�5=2+ ��1
�1� V 2t
�3i� � dWt;
con � = �M .
Demostración. Utilizaremos una vez más la Proposición 68. Notemos que la ecuación
(5.16) es del tipo (5.1) con
�� (p) = �� = e�� (v) ;D (p) =
��
�2�1 + �p0
�= ��
�1� v2
�1=2�2
+M2
�
!�1� v2
��1=2= eD (v) ;
bajo la identidad v = pp0. Aplicando la Proposición 68 a la ecuación (5.16) se llega al resultado
deseado.
Para el coe�ciente de difusión efectiva obtenemos la siguiente expresión.
Proposición 84 La posición Xt asociada al proceso RBM(I) dado por la ecuación (5.16), tiene
como constante de difusión asintótica a
D1 (X) = [��K1 (�)]�1Z 1
0
e��(1�v2)�1=2
� (1� v2) + (1� v2)3=2dv;
para � = �M y y K1 la función de Bessel modi�cada de segunda clase de orden 1.
Observación 85 Cabe señalar que la integral que aparece en esta proposición no tiene una
forma explícita y se evalúa de manera numérica.
Demostración. (Proposición 84). Procederemos como en las pruebas de las Proposiciones
76 y 81. Nótese que la velocidad del proceso RBM(I) satisface la ecuación (5.17), la cual tiene
96
como proceso equivalente a una ecuación de la forma (5.6) con
b (v) =���
h�1� v2
�5=2+ ��1
�1� v2
�3i;
a� (v) = ��
��1� v2
�+3
�
�1� v2
�3=2+3
�2�1� v2
�2�� �� � 5�
�1� v2
�3=2+6
�2�1� v2
�2�= ��
��1� v2
�� 2
�
�1� v2
�3=2 � 3
�2�1� v2
�2�;
para cualquier v 2 [�1; 1]. Entonces, según la notación del Teorema 70
�� (v) = ��
��1� v2
�+3
�
�1� v2
�3=2+3
�2�1� v2
�2�v;
U (v) = �
Z v
0
h�1� w2
�+ 3
�
�1� w2
�3=2+ 3
�2
�1� w2
�2iwh
(1� w2)5=2 + ��1 (1� w2)3i dw
= �
Z v
0
w
(1� w2)3=2
241 + 3�
�1� w2
�1=2+ 3
�2
�1� w2
�1 + ��1 (1� w2)1=2
35 dw= �
Z v
0
w
(1� w2)3=2
24 1
1 + ��1 (1� w2)1=2+
3�
�1� w2
�1=2+ 3
�2
�1� w2
�1 + ��1 (1� w2)1=2
35 dw= �
Z v
0
"w
(1� w2)3=2 + ��1 (1� w2)2+
3w
� (1� w2)
#dw
= �
Z v
0
w
(1� w2)3=2 + ��1 (1� w2)2dw � 3
2log�1� v2
�
= �
264�1� w2��1=2 + 1
2�log
0B@ 1� w2h1 + ��1 (1� w2)1=2
i21CA375w=v
w=0
� 32log�1� v2
�
=�
(1� v2)1=2+1
2log
0B@ 1� v2h1 + ��1 (1� v2)1=2
i21CA� �+ log��+ 1
�
�� 32log�1� v2
�
=�
(1� v2)1=2� �+ log
"1
(1� v2) + ��1 (1� v2)3=2
#+ log
��+ 1
�
�
= log
"�+ 1
� (1� v2) + (1� v2)3=2
#� �
h1�
�1� v2
��1=2i;
97
para v 2 [0; 1).
Considerando las demostraciones de las Proposiciones 76 y 81, y aplicando la ecuación (5.7),
la Observación 71 y el Lema 75, se llega a que
D1 (X) =e��2
�� (�+ 1)
R 10 e
U(y)nR 1
y x exph���1� x2
��1=2i �1� x2
��3=2dxo2dyR 1
0 exph�� (1� z2)�1=2
i(1� z2)�3=2 dz
= � e�
�� (�+ 1)dd�K0 (�)
Z 1
0eU(y) exp
h�2�
�1� y2
��1=2idy
=1
��K1 (�)
Z 1
0
exph���1� y2
��1=2i� (1� y2) + (1� y2)3=2
dy
= [��K1 (�)]�1Z 1
0
e��(1�v2)�1=2
� (1� v2) + (1� v2)3=2dv;
lo que demuestra lo deseado.
El ROUP y los dos modelos RBM son casos límites especiales de la ecuación de Langevin
(5.1) con coe�cientes de fricción elegidos arbitrariamente. Para sistemas realistas, la forma
exacta del coe�ciente de fricción se determina por las interacciones microscópicas. Por tanto
se puede obtener información sobre las fuerzas microscópicas inherentes al realizar mediciones
simultáneas de la temperatura y las constantes de difusión (asintótica). Más adelante, en el
Capítulo 7, se dará un método general para deducir coe�cientes de fricción más realistas a
partir de modelos microscópicos.
98
Capítulo 6
Reparametrizaciones de la Ecuación
de Langevin
En los capítulos anteriores sólo se han considerado ecuaciones de Langevin parametrizadas
en términos del tiempo t del marco de referencia del laboratorio �. En este Capítulo se estudian
las reparametrizaciones de la ecuación de Langevin en términos del tiempo propio de la partícula
y del tiempo de un observador en movimiento. Lo primero se hace en la Sección 1 mientras que
lo segundo en la Sección 2.
6.1. Reparametrización en términos del tiempo propio
En esta sección se estudia como la ecuación de Langevin, en términos del tiempo t del marco
de referencia �, puede ser reparametrizada en términos del tiempo propio � de la partícula.
Para este propósito, se considera la ecuación de Langevin d-dimensional con la integral de Itô
dX�t =
�P�tP 0t
�dt, (6.1a)
dP it = Ai (t;Xt;Pt) dt+ cij (t;Xt;Pt) � dW j
t ; (6.1b)
para i = 1; : : : ; d, � = 0; 1; : : : ; d, Xt =�X1t ; : : : ; X
dt
�y Pt =
�P 1t ; : : : ; P
dt
�. Además, en la
ecuación anterior,Wt =�W 1t ; : : : ;W
dt
�es un proceso de Wiener d-dimensional y las funciones
Ai; cij : [0;1)� Rd � Rd ! R son de clase C1.
99
A partir de la fórmula de Itô, se llega a que la ecuación para la energía relativista P 0t =�M2 +P2t
�1=2 del proceso (6.1) satisfacedP 0t =
(AitP itP 0t
+Dij
2
"�ijP 0t
� P itPjt�
P 0t�3#)
dt+P itP 0t
cir � dW rt , (6.2)
con Dij := circjr =Pdr=1 c
ircjr, en donde hemos usado la convención de Einstein. Además,
si la función de densidad de transición del proceso (6.1) existe1, la ecuación de Fokker-Planck
para la densidad del proceso (Xt;Pt) es�@
@t+pi
p0@
@xi
�f (t;x;p) =
@
@pi
��Ai (t;x;p) f (t;x;p) + 1
2
@
@pk
�Dik (t;x;p) f (t;x;p)
��,
para cualquier (t;x;p) 2 (0;1)� Rd � Rd.
Se desea reescribir las ecuaciones (6.1) en términos del tiempo propio � , el cual es un proceso
estocástico de�nido por
� (t) :=
Z t
0
�1�V2
s
�1=2ds =
Z t
0
M
P 0sds 8t � 0.
Para lograr nuestro objetivo, conviene observar que � (t) tiene inversa (para todo !), a la cual
denotaremos por bX0� = t (�). Esta función inversa representa el tiempo de la partícula en el
marco de referencia � parametrizado por el tiempo propio � . Entonces, se puede observar que,
una forma de lograr el objetivo perseguido es encontrando ecuaciones diferenciales estocásticas
para los procesos reparametrizados bX�� := X�
t(�) ybP�� := P�t(�) en �. Precisamente, esta es la
forma como se procederá. Estableceremos primero unos resultados técnicos preliminares.
Lema 86 Los procesos P 0t :=�M2 +P2t
�1=2 y cij (t;Xt;Pt) cumplen las siguientes propiedades:(i) 0 <
�P 0tM
�1=2<1 8t � 0;
(ii)�P 0tM
�1=2es un proceso no-anticipante;
1Más precisamente, la ecuación (6.1), además de tener densidad de transición, debe cumplir las condiciones delos Teoremas 29, 37 y 107 y que su condición inicial tenga distribución P0. Esto asegura la validez de la ecuaciónde Fokker-Planck.Esta misma observación aplica para todos los casos donde se hable de la ecuación de Fokker-Planck en este
capítulo.
100
(iii)R t0MP 0sds <1 8t � 0.
Demostración. Los incisos (i) y (ii) se cumplen debido a que el proceso fPt; t � 0g es
no-anticipante y toma valores en Rd. Para comprobar el inciso (iii), basta notar que MP 0t
=
M
(M2+P2t )1=2 � M
M = 1, lo que implica
Z t
0
M
P 0sds �
Z t
0ds = t <1 8t � 0:
Esto concluye la demostración.
Lema 87 Los procesos cij (t;Xt;Pt) cumplen la igualdad
Z t(�0)
0cij (s;Xs;Ps) dW
js =
Z �0
0cij� bX0
t(r);Xt(r);Pt(r)
� P 0t(r)M
!1=2dcW j
r 8�0 2 [0; � (1)) ;
donde cW jt :=
R t(t)0
�MP 0s
�1=2dW j
s ; para todo t 2 [0;1).
Observación 88 Una consecuencia inmediata del Teorema 38 es que, para cada j = 1; :::; d,cW jt :=
R t(t)0
�MP 0s
�1=2dW j
s de�ne un nuevo proceso de Wiener.
Demostración. Basaremos nuestra demostración en el Teorema 40. Del Lema 86 se sigue
que podemos aplicar el Teorema 40 a las integralesR t0 c
ij (s;Xs;Ps) dWjs y al tiempo intrínseco
� (t) : Entonces, por el Teorema 40 tenemos
Z t(�0)
0cij (s;Xs;Ps) dW
js =
dXj=1
Z t(�0)
0cij (s;Xs;Ps) dW
js
=
dXj=1
Z �0
0cij�t (r) ;Xt(r);Pt(r)
� P 0t(r)M
!1=2dcW j
r
=
dXj=1
Z �0
0cij�X0t(s);Xt(r);Pt(r)
� P 0t(r)M
!1=2dcW j
r
=
Z �0
0cij�X0t(s);Xt(r);Pt(r)
� P 0t(r)M
!1=2dcW j
r ;
101
para todo 0 � �0 < �(1), donde cW jt :=
R t(t)0
�MP 0s
�1=2dW j
s para todo t 2 [0;1). Esto último
prueba lo deseado.
A partir del Lema anterior es posible encontrar la reparametrización de la ecuación de
Langevin en términos del tiempo propio, sin embargo es necesario demostrar primero que cWt :=�cW 1t ; : : : ;cW d
t
�es un proceso de Wiener.
Lema 89 Sea cW jt :=
R t(t)0
�MP 0s
�1=2dW j
s para cada t � 0 y j 2 f1; : : : ; dg. Entonces, cWt :=�cW 1t ; : : : ;
cW dt
�8t � 0 es un proceso de Wiener d-dimensional.
Demostración. Para demostrar lo deseado, aplicaremos los Teoremas 42 y 43. De�namos
Y jt :=R t0
�MP 0s
�1=2dW j
s para cualesquiera t � 0 y j 2 f1; : : : ; dg. Del Teorema 42, se sigue que
�Y i; Y j
�t=
Z t
0
M
P 0sd�W i;W j
�s
8t � 0:
Luego, se cumple que
hcW i;cW jit=�Y i; Y j
�t(t)
=
Z t(t)
0
M
P 0sd�W i;W j
�s=
8<: t si i = j
0 si i 6= j;
para cualesquiera t � 0; i; j 2 f1; : : : ; dg e i 6= j. Entonces, la Observación 22 y el Teorema 43
implican que cWt :=�cW 1
t ; : : : ;cW dt
�es un proceso de Wiener d-dimensional, que es lo que se
quería demostrar.
El siguiente teorema es el resultado más importante de esta sección. Establece la reparame-
trización de la ecuación de Langevin (6.1) en términos del tiempo propio � .
Teorema 90 Supongamos que se cumple la ecuación (6.1). Entonces, los procesos bX0� :=
X0t(�),
bX� := Xt(�) y bP� := Pt(�) satisfacen la ecuaciónd bX�
� =
bP��M
!d�; (6.3a)
d bP i� = bAi � bX0� ;bX� ;bP�� d� + bcij � bX0
� ;bX� ;bP�� � dcW j
� ; (6.3b)
102
para i = 1; :::; d y � = 0; 1; :::; d; donde bAi; bcij y cW jt están dados por
bAi � bX0� ;bX� ;bP�� : =
bP 0�M
!Ai� bX0
� ;bX� ;bP�� ;
bcij � bX0� ;bX� ;bP�� : =
bP 0�M
!1=2cij� bX0
� ;bX� ;bP�� ;
cW jt : =
Z t(t)
0
�M
P 0s
�1=2dW j
s 8t 2 [0;1) :
Observación 91 Conviene notar que por el Lema 89 se tiene que cWt :=�cW 1
t ; : : : ;cW dt
�es un
proceso de Wiener d-dimensional.
Demostración. Basaremos nuestra demostración en el Lema 87. Aplicando el cambio de
variable s = t (r) a la ecuación (6.1a), se llega a que
bX��0 = X�
t(�0)= X�
0 +
Z t(�0)
0
�P�sP 0s
�ds
= X�0 +
Z �0
0
P�t(r)
P 0t(r)
! P 0t(r)
M
!dr
= X�0 +
Z �0
0
P�t(r)
Mdr
= bX�0 +
Z �0
0
bP�rM
dr; (6.4)
para todo 0 � �0 < � (1) : Además, del Lema 87 y el cambio de variable anterior, se sigue
bP i�0 = P it(�0) = P i0 +
Z t(�0)
0Ai (s;Xs;Ps) ds+
Z t(�0)
0cij (s;Xs;Ps) � dW j
s
= P i0 +
Z �0
0Ai�t (r) ;Xt(r);Pt(r)
� P 0t(r)M
!dr
+
Z �0
0cij�X0t(r);Xt(r);Pt(r)
� P 0t(r)M
!1=2dcW j
r
= bP i0 + Z �0
0
bP 0rM
!Ai� bX0
r ;bXr; bPr� dr
+
Z �0
0cij� bX0
r ;bXr;bPr� bP 0r
M
!1=2dcW j
r
103
= bP i0 + Z �0
0
bAi � bX0� ;bX� ;bP�� d� + Z �0
0bcij � bX0
� ;bX� ;bP�� dcW j
� , (6.5)
para cualquier 0 � �0 < � (1), donde bAi; bcij y cW jt se toman como en las hipótesis del Teorema.
De las ecuaciones (6.4) y (6.5) se siguen las igualdades (6.3a) y (6.3b), que es lo que
queríamos probar.
Si la función de densidad de transición de la ecuación (6.1) existe, tenemos que el proceso
(6.3) también tiene función de densidad de transición bf �� ; x0;x;p�, la cual satisface la ecuaciónde Fokker-Planck
�@
@�+p�
M
@
@x�
� bf �� ; x0;x;p� =@
@pi��Ai
�� ; x0;x;p
�f�� ; x0;x;p
�+1
2
@
@pk
�Dik
�� ; x0;x;p
�f�� ; x0;x;p
���,
para cualquier�� ; x0;x;p
�2 (0;1)� (0;1)� Rd � Rd. La cantidad bf �� ; x0;x;p� dx0ddxddp
se interpreta como la probabilidad de encontrar la partícula al tiempo propio � en el espacio
(2d+ 1)-dimensional [t; t+ dt]��Yd
i=1
�xi; xi + dxi
����Yd
i=1
�pi; pi + dpi
��del marco de
referencia inercial �.
Una consecuencia interesante de esta reparametrización, surge al considerar el movimiento
libre (sin la presencia de campos de fuerza) en un baño térmico termalizado, estacionario,
isotrópico y homogéneo en el marco �. Como ya se ha comentado, en este caso la ecuación de
Langevin toma la forma
dX�t =
P itP 0t
dt; (6.6a)
dP it = �b� �P 0t �P it + h bD �P 0t �i1=2 �ij � dW jt ; (6.6b)
donde los coe�cientes de ruido y fricción b� y bD dependen de la energía p0 =�M2 + p2
�1=2 ysatisfacen la relación de Einstein
0 � b� �p0� p0 + bD0 �p0�� �D �p0� 8p0 2 R� [�M;M ] .
En esta situación, la densidad estacionaria (marginal) del momento f1 (p) para el proceso (6.6)
104
es una función de Jüttner, es decir
f1 (p) := l��mt!1
f (t;p) /exp
���p0
�Zj
= �J (p) ;
siendo f la densidad marginal del momento del proceso (6.6).
Para algunos sistemas, es realista suponer que la velocidad de la partícula Browniana está
acotada por una cantidad menor que c = 1. En estos casos, � (1) =1 y la densidad estacionaria
marginal del momento bf1 (p) correspondiente a la reparametrización en términos de � delproceso (6.6), está dada por
bf1 (p) := l��m�!1
bf (� ;p) / exp���p0
�p0Zj
= �MJ (p) ;
donde bf es la densidad marginal del momento de la reparametrización del proceso (6.6).Así, el proceso (6.6) y su reparametrización en términos de � tienen funciones de densidad
de transición distintas, a pesar de que modelan el mismo fenómeno. Físicamente, la diferencia
entre f1 y bf1 se debe a que las mediciones en t = const y � = const no son equivalentes,
incluso cuando t; � !1.
Usualmente, la función �MJ recibe el nombre de densidad de Jüttner modi�cada. Como
puede observarse, �MJ di�ere de la función de Jüttner estándar �J por un factor proporcional
al inverso de la energía. La función de Jüttner modi�cada �MJ puede ser derivada de un
principio de entropía relativa, usando el hecho de que la densidad de transición es un invariante
de Lorentz.
6.2. Observadores en movimiento
Hasta ahora, solo se han tratado en este trabajo ecuaciones de Langevin que describen
el movimiento aleatorio de una partícula Browniana relativista en el marco de referencia del
laboratorio �. En esta sección, se determina como percibe un observador en movimiento (a
velocidad constante w) el desplazamiento de la partícula Browniana, suponiendo que se tiene
una ecuación de Langevin de la forma (5.1) en el marco �.
Hay, al menos, dos maneras distintas de abordar este problema. Una de ellas utiliza la inva-
105
rianza de Lorentz de la densidad de transición. La otra forma consiste en aplicar directamente
una transformación de Lorentz a la ecuación de Langevin medida en �.
6.2.1. Transformación de Lorentz de la densidad de transición en el espacio
fase
Un observador en reposo en el marco � puede medir una densidad f (t;x;p) ; asociada al
proceso (Xt;Pt) ; que satisface las ecuaciones de Fokker-Planck (4.4) ó (4.6). Mientras que un
observador que se mueve a través del marco � a velocidad constante w mide una densidad
f 0 (t0;x0;p0). Como la densidad f resulta ser un escalar de Lorentz (ver Apéndice D.1), se tiene
que
f 0�t0;x0;p0
�= f
�t�t0;x0
�;x�t0;x0
�;p�p0��
(6.7)
y
f (t;x;p) = f 0�t0 (t;x) ;x0 (t;x) ;p0 (p)
�; (6.8)
donde (t0;x0;p0) y (t;x;p) están relacionados por
x0� (t;x) = ��0t+ ��ixi;
p0i (p) = �i0�M2 + p2
�1=2+ ��jpj ;
siendo ��� una transformación de Lorentz del marco � al marco del observador en movimiento.
Para calcular f 0 es su�ciente resolver las ecuaciones de Fokker-Planck (4.4) que se cumplen
en el marco �, para una condición inicial t-simultánea f (0;x;p), e insertar tal solución en
(6.7). Obteniendo f 0 se logra determinar la forma en que el observador en movimiento (velocidad
constante w) percibe el movimiento de la partícula Browniana, con lo cual completamos nuestro
objetivo.
6.2.2. Transformación de Lorentz de la ecuación de Langevin
Considérese un marco de referencia inercial �0, el cual se encuentra en movimiento respecto
a �. En esta subsección se encuentran ecuaciones diferenciales estocásticas explícitas para el
movimiento Browniano en �0, suponiendo que la ecuación de Langevin (6.1) se cumple en el
106
marco �. Más precisamente, se obtienen ecuaciones diferenciales parametrizadas por el tiempo
t0 (de �0) para la posición X0t0 y el momento P0t0 medidos en �
0: Esto se logra aplicando una
transformación de Lorentz a las ecuaciones de Langevin que se cumplen en �.
Con objeto de establecer las ecuaciones diferenciales estocásticas deseadas para X0t0 y P0t0 ,
conviene comenzar estableciendo igualdades que relacionen a�t0;X0t0 ;P
0t0�y (t;Xt;Pt). Con este
propósito, considérese una transformación de Lorentz � = ��� del marco del laboratorio � al
marco �0. Por conveniencia, se excluyen transformaciones de Lorentz que envuelven inversión
espacial o reversión temporal, es decir, nos restringimos a transformaciones propias de Lorentz
con �00 > 0. Entonces, si se de�ne
Y 0�t := ���X�t , G0�t := �
��P�t ,
para � = 0; 1; : : : ; d y t � 0; se tiene que los procesos�X0t0 ;P
0t0�y (Xt;Pt) cumplen las relaciones
X 0�t00= Y 0�
t(t00), P 0�t00
= G0�t(t00)
;
para � = 0; 1; : : : d y donde se toma t := (t0)�1. Además, también se tiene que
dt0 (t) = dY 00t =�0�P�tP 0t
dt =G00tP 0t
dt =P 00t0(t)
(��1)0� P 0�t0(t)dt; (6.9)
donde ��1 es la inversa de la transformación de Lorentz.
Con lo anterior, procederemos a encontrar las ecuaciones diferenciales estocásticas para
X0t0 y P0t0 . La idea a seguir es muy similar a la utilizada para encontrar la reparametrización
en términos del tiempo propio. Se comienza encontrando ecuaciones diferenciales estocásticas
para los procesos�Y 0�t0�y�G0�t0�: Posteriormente, se utiliza el Teorema 40 para establecer las
ecuaciones para�X 0�t0�y�P 0�t0�.
A continuación se establecen ecuaciones diferenciales para los procesos�Y 0�t0�y�G0�t0�: Se
inicia con un resultado que simpli�ca el encontrar tales ecuaciones, el cual es una consecuencia
inmediata de la ecuación (6.2).
Lema 92 Supongamos que la ecuación (6.1) se cumple. Entonces, la energía relativista P 0t :=
107
�M2 +P2
�1=2 tiene asociada la diferencialdP 0t = A0 (t;Xt;Pt) dt+ c
0j (t;Xt;Pt) dWjt ;
donde
A0 (t;x;p) : = Ai (t;x;p)pi
p0+Dij (t;x;p)
2
��ijp0� pipj
(p0)3
�;
c0j (t;x;p) : =pi
p0cij (t;x;p) ;
para cualquier (t;x;p) 2 [0;1)� Rd � Rd.
Proposición 93 Supongamos que la ecuación (6.1) se cumple. Entonces, los procesos Y 0�t :=
���X�t y G
0�t := �
��P�t satisfacen las ecuaciones diferenciales
dY 0�t =G0�tP 0t
dt; (6.10a)
dG0�t = �i�A� (t;Xt;Pt) dt+ �i�c�j (t;Xt;Pt) dW
jt : (6.10b)
Demostración. Comenzaremos probando (6.10b); para demostrarlo aplicaremos la fórmula
de Itô. De�namos u : Rd+1 ! Rd+1 dada por u�p0;p
�= �(p�).
Por ser � una transformación lineal, se tiene que u satisface las hipótesis necesarias para
aplicarle la fórmula de Itô d+1-dimensional. Entonces, si fe1; :::; ed+1g forma la base canónica
de Rd+1, de la fórmula de Itô 26 y el Lema 92, se sigue que el proceso G0�t := ���P�t = u
�P 0t ;Pt
�satisface la ecuación
dG�t =dX�=0
�@
@x�u�P 0t ;Pt
�; e���dP�t +
1
2
dX�;�=0
�@2
@x�@x�u�P 0t ;Pt
�; e���dP �t dP
�t
= ���dP�t
= ���A� (t;Xt;Pt) dt+ ���c�j (t;Xt;Pt) dW
jt ;
lo que prueba (6.10b). Para demostrar (6.10a) se puede proceder de manera similar.
108
Utilizando el resultado anterior, se encuentran las ecuaciones para�X 0�t0�y�P 0�t0�: Se es-
tablecen primero unos resultados técnicos preliminares.
Lema 94 Los procesos P 0t :=�M2 +P2t
�1=2 y cij (t;Xt;Pt) cumplen las siguientes propiedades:(i) 0 <
�P 0t
�0�P�t
�1=2<1 8t � 0;
(ii)�
P 0t�0�P�t
�1=2es un proceso no-anticipante;
(iii)R t0�0�P�tP 0t
ds <1 8t � 0.
Demostración. El inciso (i) se cumple debido a que P2t ; (P0t0(t))
2 2 R+[f0g 8t � 0 y a queP 0t
�0�P�t=
P 0tP 00t0(t)
=(M2+P2t )
1=2�M2+
�P0t0(t)
�2�1=2 8t � 0: El inciso (ii) es una consecuencia de que el procesofPt; t � 0g sea no-anticipante y de que la transformación � : Rd+1 ! Rd+1 sea continua.
Resta demostrar el inciso (iii). Sea t0 � 0. Debido a que s 7! �0�P�sP 0s
es una composi-
ción de funciones continuas, tenemos que tal función esta acotada en [0; t0] 8!. Por lo tanto,R t00
�0�P�sP 0s
ds <1, lo cual implica queR t0�0�P�sP 0s
ds <1 8t � 0, lo que demuestra el inciso (iii).
Con esto concluye la demostración.
Proposición 95 El tiempo t0 en �0 satisface t0 (1) =1.
Demostración. Debido a que � = ��� es una transformación propia de Lorentz, por el
Apéndice C (ecuación (C.7)) sabemos que
� = BR;
donde R es una rotación de la forma (C.5) y B es un boost de Lorentz de la forma (C.6) con
=�1�w2
��1=2. Luego, tenemos que��X0t ;Xt
�= BR
0@ X0t
Xt
1A = B
0@ t
(Rij)Xt
1A :
109
De esto último y la desigualdad de Cauchy-Schwarz, se sigue que
t0 (t) � t0 (t;Xt) = t� h(Rij)Xt; wi
� t� k(Rij)Xtk k wk
= (t� kXtk kwk)
� (t� (t+X0) kwk)
= t (1� kwk)� X0 kwk : (6.11)
Entonces, debido a que consideramos a X0 una condición inicial determinista, de la ecuación
(6.11) es claro que t0 (1) =1, que prueba lo deseado.
Sin perder generalidad se puede suponer queX0 = 0. Esta suposición permite establecer una
ecuación diferencial estocástica �sencilla� para describir el movimiento Browniano relativista
en �0.
Lema 96 Supongamos que se cumplen las ecuaciones (6.1). Los procesos cij (t;Xt;Pt) cumplen
la relación
Z t(t00)
0�i�c�j (s;Xs;Ps) dW
js =
Z t00
0�i�c�j
�t (r) ;Xt(r);Pt(r)
� P 0t(r)G00t(r)
!1=2dW 0j
r 8t00 � 0;
donde W 0jt :=
R t(t)0
�G00tP 0t
�1=2dW j
s ; para todo t 2 [0;1).
Observación 97 Una consecuencia inmediata del Teorema 38 es que, para cada j = 1; :::; d,
W 0jt :=
R t(t)0
�G00tP 0t
�1=2dW j
s de�ne un nuevo proceso de Wiener.
Demostración. Basaremos nuestra demostración en el Teorema 40. Por ser X0 = 0, tene-
mos que t0 (0) = 0. Entonces, de las ecuaciones (6.9), se sigue que
t0 (t) =
Z t
0
�0�P�sP 0s
ds =
Z t
0
G00sP 0s
ds =
Z t
0
P 00t0(s)
(��1)0� P 0�t0(s)ds 8t � 0:
Teniendo en cuenta la ecuación anterior, del Lema 94 se sigue que podemos aplicar el
Teorema 40 a las integralesR t0 �
i�c�j (s;Xs;Ps) dWjs y al tiempo intrínseco t0 (t). Entonces, por
110
el Teorema 40 tenemos que
Z t(t00)
0�i�c�j (s;Xs;Ps) dW
js =
dXj=1
Z t(t00)
0�i�c�j (s;Xs;Ps) dW
js
=dXj=1
Z t00
0�i�c�j
�t (r) ;Xt(r);Pt(r)
� P 0t(r)G00t(r)
!1=2dW 0j
r
=
Z t00
0�i�c�j
�t (r) ;Xt(r);Pt(r)
� P 0t(r)G00t(r)
!1=2dW 0j
r ;
para todo 0 � t00 < t0 (1) = 1, donde W 0jt :=
R t(t)0
�G00tP 0t
�1=2dW j
s para todo t 2 [0;1). Esto
último prueba lo deseado.
A partir del Lema anterior es posible encontrar las ecuaciones que describen el movimiento
de la partícula Browniana en �; sin embargo, al igual que en la reparametrización en términos
del tiempo propio, es necesario demostrar queW0t :=
�W 01t ; : : : ;W
0dt
�es un proceso de Wiener
d-dimensional.
Lema 98 Sea W 0jt :=
R t(t)0
�G00tP 0t
�1=2dW j
s para cada t � 0 y j 2 f1; : : : ; dg. Entonces,�W0
t :=�W 01t ; : : : ;W
0dt
�; t � 0
es un proceso de Wiener d-dimensional.
Demostración. Procederemos como en la demostración del Lema 98. De�namos Y jt :=R t0
�G00tP 0t
�1=2dW j
s =R t0
��0�P�tP 0t
�1=2dW j
s para cualesquiera t � 0 y j 2 f1; : : : ; dg. Por el Teorema
42, tenemos que �Y i; Y j
�t=
Z t
0
G00tP 0t
d�W i;W j
�s
8t � 0:
Luego, como supusimos X0 = 0, se cumple que
�W 0i;W 0j�
t=�Y i; Y j
�t(t)
=
Z t(t)
0
G00tP 0t
d�W i;W j
�s=
8<: t si i = j
0 si i 6= j;
para cualesquiera t � 0; i; j 2 f1; : : : ; dg e i 6= j. Entonces, la Observación 22 y el Teorema 43
implican que W0t :=
�W 01t ; : : : ;W
0dt
�es un proceso de Wiener d-dimensional, que es lo que se
quería probar.
El siguiente teorema es el resultado más importante de esta subsección. Establece las ecua-
111
ciones que describen el movimiento de la partícula Browniana en �0.
Teorema 99 Supongamos que se cumple la ecuación (6.1). Entonces, los procesos X 0�t00= Y 0�
t(t00)
y P 0�t00= G0�
t(t00)satisfacen las ecuaciones
dX 0�t0 =
�P 0�t0M
�dt0; (6.12a)
dP 0it0 = A0i�t0;X0t0 ;P
0t0�d� + c0ij
�t0;X0t0 ;P
0t0�� dW 0j
t0 ; (6.12b)
para i = 1; : : : ; d y � = 0; 1; : : : ; d; donde A0i, c0ij y W 0jt0 están dados por
A0i�x00;x0;p0
�: =
"���1
�0�p0�
p00
#�i�A�
�����1
���x0��d�=0
;h���1
�j�p0�idj=1
�
c0ij�x00;x0;p0
�: =
"���1
�0�p0�
p00
#1=2�i�c�j
�����1
���x0��d�=0
;h���1
�j�p0�idj=1
�
W 0jt : =
Z t(t)
0
�G00tP 0t
�1=2dW j
s ; 8t 2 [0;1) :
Observación 100 Hay dos hechos implícitos importantes en la de�nición de A0i y c0ij:
1. Dado un vector�z0; z1; : : : ; zd
�2 Rd+1 y un número m 2 f0; 1; : : : ; dg, de�nimos
hzlidl=m
:=�zm; zm+1; : : : ; zd
�2 Rd+1�m:
Luego,�zl�dl=m
es un vector en Rd+1�m.
2. En un pequeño abuso de notación, estamos entendiendo
����1
���x0��d�=0
=
����1
�0�x0�;
h���1
�j�x0�idj=1
�:
Entonces, tiene sentido evaluar Ai y cij en���
��1���
x0��d�=0
;h���1
�j�p0�idj=1
�: Hace-
mos uso de esta notación para facilitar el tratamiento, así como para simpli�car y com-
pactar expresiones.
Demostración. (Teorema 99). Basaremos nuestra demostración en las Proposiciones 93 y
112
96. Aplicando el cambio de variable s = t (r) a la ecuación (6.10a), se llega a que
X 0�t00
= Y 0�t(t00)
= Y 0�0 +
Z t(t00)
0
G0�sP 0s
ds
= Y 0�0 +
Z t00
0
G0�t(r)P 0t(r)
P 0t(r)
G00t(r)dr
= Y 0�0 +
Z t00
0
G0�t(r)G00t(r)
dr
= X 0�0 +
Z t00
0
P 0�rP 00r
dr; (6.13)
para todo t00 � 0: Además, del Lema 96 y el cambio de variable anterior, se sigue
P 0it00= G0i
t(t00)= G0i0 +
Z t(t00)
0�i�A� (s;Xs;Ps) ds+
Z t(t00)
0�i�c�j (s;Xs;Ps) dW
js
= G0i0 +
Z t00
0�i�A�
�t (r) ;Xt(r);Pt(r)
� P 0t(r)G00t(r)
dr
+
Z t00
0�i�c�j
�t (r) ;Xt(r);Pt(r)
� P 0t(r)G00t(r)
!1=2dW 0j
r
= P 0i0 +
Z t00
0�i�A�
�����1
���x0��d�=0
;h���1
�j�p0�idj=1
� ���1
�0�P 0�r
P 00rdr
+
Z t00
0�i�c�j
�����1
���x0��d�=0
;h���1
�j�p0�idj=1
� ���1
�0�P 0�r
P 00r
!1=2dW 0j
r
= P 0i0 +
Z t00
0A0i�t;X0t;P
0t
�dr +
Z t00
0c0ij�t;X0t;P
0t
�dW 0j
r ; (6.14)
para cualquier t00 � 0, donde A0i, c0ij y W 0j se toman como en las hipótesis del problema.
De (6.13) y (6.14) se siguen las ecuaciones (6.12), que es lo que queríamos probar.
El Teorema anterior muestra que la trayectoria�X0t0 ;P
0t0�en �0 de la partícula Browniana,
satisface una ecuación diferencial de la forma
dX 0�t0 =
�P 0�t0
P 00t0
�dt0;
dP 0it0 = A0idt0 + c0ij � dW 0jt0 :
Además, se puede encontrar la ecuación diferencial estocástica para la energía relativista P 00t =
113
�M2 +P02t
�1=2 en �0 reemplazando las cantidades sin prima (0) por cantidades con prima en laecuación (6.2).
Debido a los prefactores�(��1)
0�p0�
p00
�en las de�niciones de A0i y c0ij , se tiene que los coe-
�cientes Ai y cij no constituyen tensores de Lorentz. Esta es, esencialmente, una consecuencia
de la parametrización temporal de las ecuaciones (6.1) y (6.12).
Si la función de densidad de transición de la ecuación (6.1) existe, tenemos que el proceso
(6.12) también tiene función de densidad de transición f 0 (t0;x0;p0), la cual satisface la ecuación
de Fokker-Planck
�@
@t0+p0i
p00@
@x0i
�f 0�t0;x0;p0
�=
@
@p0i��A0i
�t0;x0;p0
�f 0�t0;x0;p0
�+1
2
@
@p0k
�D0ik �t0;x0;p0� f 0 �t0;x0;p0��� ,
para cualquier (t0;x0;p0) 2 (0;1)� Rd � Rd, donde D0ij := c0irc0jr.
Cuando se considera el movimiento libre en un baño térmico termalizado, estacionario,
isotrópico y homogéneo en el marco �, la ecuación de Langevin en � toma la forma
dX�t =
P�tP 0t
dt;
dP it = �b� �P 0t �P it + h bD �P 0t �i1=2 �ij � dW jt ;
donde los coe�cientes de ruido y fricción b� y bD dependen de la energía p0 =�M2 + p2
�1=2 y, sila función de densidad estacionaria del proceso (6.6) existe, satisfacen la relación de Einstein
0 � b� �p0� p0 + bD0 �p0�� �D �p0� 8p0 2 R� [�M;M ] .
En esta situación, si la transformación de Lorentz es solo un boost, la densidad estacionaria
marginal del momento en el marco en movimiento �0; esta dada por una función de Jüttner
transformada de la forma
�0�p0�/ exp
���U 0�p0�
�;
donde U 0� es la cuadri-velocidad del baño térmico en �0.
114
Capítulo 7
El Modelo de Colisión Binaria
Relativista
Para que el enfoque de Langevin se apegue a la realidad, se deben conocer los coe�cientes
de ruido D (p) y fricción � (p) que son apropiados para el sistema bajo consideración. Algunos
modelos de fricción realistas pueden obtenerse, por ejemplo, derivando ecuaciones de Fokker-
Planck de las ecuaciones relativistas de Boltzmann. Sin embargo, los métodos conocidos en la
literatura resultan ser complicados y son poco entendidos.
En este capítulo bosquejaremos un procedimiento alternativo para obtener coe�cientes de
ruido y fricción a partir de un modelo de interacción microscópica simple y bajo el supuesto de
equilibrio térmico. Tal procedimiento puede considerarse como la generalización relativista del
modelo de colisión binaria presentado en el Capítulo 3, así como una justi�cación heurística de
la ecuación de Langevin relativista.
En la exposición de este capítulo se considera un sistema unidimensional que consiste de una
partícula Browniana pesada (masaM) inmersa en un baño térmico de partículas más pequeñas
(masa m�M , número total de partículas N � 1). El modelo que se desarrolla supone que el
movimiento aleatorio de una partícula Browniana surge debido a frecuentes colisiones elásticas
con las partículas del baño térmico circundante.
De manera similar al modelo de colisión binaria no-relativista, el interés primordial radica
en encontrar la mejor aproximación (posible) de la dinámica exacta en la clase de ecuaciones
115
diferenciales estocásticas de�nidas por la ecuación de Langevin-backward (5.1).
En la primera sección se desarrolla la cinemática de las colisiones relativistas, se construye
una ecuación del tipo Langevin teniendo en cuenta algunos supuestos y se comenta sobre una
posible justi�cación de la ecuación de Langevin (4.2) en ausencia de campos de fuerza externos.
La sección 2 trata la distribución del baño térmico considerada (aquella del equilibrio térmico)
y se realiza un par de comentarios sobre la fuerza media de deriva. Finalmente, en la sección
3 se estudia como aproximar, según algunos supuestos que se establecen, la �ecuación de tipo
Langevin�obtenida de la cinemática relativista por ecuaciones diferenciales estocásticas.
7.1. Cinemática de las colisiones relativistas
Considérese una (única) colisión elástica entre la partícula Browniana (momento Pt, energía
E (Pt) = P 0t ) y una partícula del baño térmico (momento pt, energía " (pt) = p0t ). La energía
relativista, el momento y la velocidad de las dos partículas antes de la colisión, están dadas por
Pt = MVt (Vt) ; E (Pt) =�M2 + P 2t
�1=2; (7.1a)
pt = mvt (vt) ; " (pt) =�m2 + p2t
�1=2; (7.1b)
donde (w) :=�1� w2
��1=2. Por tener colisiones elásticas, la cinemática de la colisión obedecelas leyes de conservación
cM =M; bm = m; E + " = bE + b"; Pt + pt = bPt + bpt; (7.2)
donde los símbolos con (b) representan el estado después de la colisión. Insertando las ecuaciones(7.1) en las leyes de conservación (7.2) y despejando el momento bPt de la partícula Brownianadespués de la colisión, se obtiene que
bPt = (ut)2 �2utE (Pt)� �1 + u2t �Pt� ;
116
donde ut es la velocidad del centro de masa (invariante bajo la colisión), el cual esta dado por
la fórmula
ut � u (pt; Pt) =Pt + pt
E (Pt) + " (pt):
Obsérvese que las partículas del baño térmico pueden enumerarse con números entre 1 y
N . Si denotamos por �P rt := bPt � Pt al cambio de momento de la partícula Browniana en unaúnica colisión con la partícula �r�(0 � r � N) del baño térmico, se tiene que
�P rt = �2 (u (prt ; Pt))2 " (prt )
E (Pt) + " (prt )Pt + 2 (u (p
rt ; Pt))
2 E (Pt)
E (Pt) + " (prt )prt ; (7.3)
donde prt es el momento de la partícula �r�al tiempo t � 0. En el caso límite no-relativista,
donde u (prt ; Pt)2 � 1, E 'M y " (prt ) ' m, la ecuación (7.3) se reduce a la ecuación (3.12).
Para construir una ecuación del tipo Langevin a partir de las ecuaciones anteriores, se
considera el cambio de momento total �P (t) := Pt+�t � Pt de la partícula Browniana en el
intervalo �mesoscópico�[t; t+ �t], asumiendo que:
las colisiones que ocurren en [t; t+ �t] pueden ser vistas como eventos independientes;
el salto de tiempo �t es lo su�cientemente pequeño para que se cumpla j�P (t) =P (t)j � 1
y para que ocurra a lo más solo una colisión entre la partícula Browniana y una partícula
especí�ca del baño térmico;
el valor de �t es lo su�cientemente grande para que el número de colisiones sea mayor que
1.
Los requisitos anteriores solo pueden cumplirse simultáneamente cuando m � M , lo cual
constituye uno de los supuestos iniciales de nuestro modelo.
Entonces se tiene que el cambio de momento �P (t) de la partícula Browniana, durante el
117
intervalo de tiempo [t; t+ �t] ; puede aproximarse por
�P (t) = Pt+�t � Pt �NXr=1
�P rt Ir (t; �t)
= �"2NXr=1
(u (prt ; Pt))2 " (prt )
E (Pt) + " (prt )Ir (t; �t)
#Pt
+2
NXr=1
(u (prt ; Pt))2 E (Pt)
E (Pt) + " (prt )prt Ir (t; �t) ; (7.4)
donde Ir (t; �t) es la función indicadora para una colisión entre la partícula Browniana y la
partícula �r�, es decir
Ir (t; �t) =
8>>><>>>:1 si una colisión entre la partícula Browniana
y la partícula r ha ocurrido en [t; t+ �t] ;
0 en otro caso.
Conviene notar que la ecuación (7.4) es la contraparte relativista de la ecuación (3.16). Heurís-
ticamente, el primer término de la ecuación (7.4) puede ser interpretado como fricción, mientras
que el segundo término puede pensarse como ruido. Esto justi�ca, de cierta manera, el porque
tomar ecuaciones de Langevin de la forma (4.2) con F � 0 para modelar la dinámica relativista
de la partícula Browniana en ausencia de campos de fuerza externos.
Por la sección 3.2.1., se tiene que heurísticamente
Ir (t; �t) ��t
2jvrt � Vtj � (xrt �Xt) ; (7.5)
siendo xrt y vrt la posición y la velocidad de la partícula �r�del baño térmico, respectivamente.
Sin embargo, en este caso, debemos usar las velocidades relativistas Vt = Pt
(M2+P 2t )1=2 y v
rt =
prt
(m2+(prt )2)1=2 .
7.2. La distribución del baño térmico y la fuerza de deriva
Como es de interés sólo el caso en el que hay equilibrio térmico, asumiremos que el baño está
termalizado. También es realista suponer que las partículas del baño térmico están con�nadas
118
a un volumen V = [0; L] � R �nito. Entonces, bajo los supuestos anteriores, la densidad que
modela la probabilidad de encontrar a una partícula especí�ca del baño térmico en algún lugar
del espacio fase, esta dada por la función de Jüttner espacialmente homogénea
f1b (x; p) = (2mK1 (�m)L)�1 exp
h���m2 + p2
�1=2i�(L� x)� (x) 8 (x; p) 2 V� R;
donde T = (�kB)�1 es la temperatura del baño, � es una función Heaviside y K1 denota la
función de Bessel modi�cada de segunda especie y orden 1.
Para completar la información de la distribución del baño térmico, también se asumirá que:
las partículas del baño térmico son independiente e idénticamente distribuidas;
las colisiones entre la partícula Browniana y las partículas del baño térmico no afectan la
distribución de las partículas.
Las dos suposiciones anteriores pueden justi�carse para baños lo su�cientemente grandes,
siempre y cuando las colisiones entre las partículas del baño restablecen rápidamente un dis-
tribución espacialmente homogénea.
Utilizaremos la distribución de las partículas del baño térmico para calcular la fuerza media
de deriva K, la cual se de�ne como
K (p) :=��P (t)
�tj Pt = p
�:
Al sustituir la ecuación (7.4) en la de�nición anterior, obtenemos
K (p) = �nb� (u (prt ; Pt))
2 " (prt )
E (Pt) + " (prt )LIr (t; �t) j Pt = p
�p
+nb
� (u (prt ; Pt))
2 E (Pt)
E (Pt) + " (prt )prtLIr (t; �t) j Pt = p
�;
donde nb = N=L es la densidad de las partículas del baño térmico.
Desafortunadamente, es muy complicado, o incluso imposible, evaluar analíticamente la es-
peranza anterior. Sin embargo, puede ser evaluada numéricamente haciendo uso de expansiones
de Taylor en el sentido de funciones generalizadas. Cabe destacar que, a diferencia del caso
119
no-relativista, la expansión de Taylor (7.5) no permite encontrar un expresión analítica para la
fuerza media de deriva K.
Conviene observar que se pueden hacer desarrollos similares si se considera una distribución
distinta para el baño térmico. Es posible que se obtengan expresiones diferentes para estos
casos, pero la forma de proceder es análoga.
7.3. Aproximación de Langevin
Para �nalizar, se estudia como, en principio, se pueden aproximar la ecuación (7.4) por
ecuaciones diferenciales estocásticas de la forma (141). Solo se estudiará a fondo el caso de la
ecuación de (5.1) con la integral backward, la cual tiene la expresión
dPt = �� (Pt)Ptdt+ [2D (Pt)]1=2 � dWt: (7.6)
Los procedimientos resultan análogos para las otras integrales estocásticas, aunque con ex-
presiones más complicadas. Se sugiere utilizar las relaciones entre las ecuaciones diferenciales
estocásticas de los Teoremas 46 y 50 para obtener ecuaciones diferenciales estocásticas con la
integral de Itô y la integral de Stratonovich que aproximen a (7.4), en lugar de aplicar el método
directo que a continuación se desarrollará.
Para obtener un modelo de Langevin útil y apegado a la realidad, los coe�cientes � y D
tienen que elegirse de tal manera que produzcan la mejor aproximación posible en la clase de
ecuaciones diferenciales estocásticas de�nida por la ecuación (7.6). Es posible que existan una
in�nidad de criterios para construir una �buena aproximación�, en este trabajo los criterios que
deberá cumplir el proceso (7.6) son dos:
Aproximar la distribución estacionaria correcta;
Exhibir el comportamiento de relajación media correcto, lo cual se puede expresar matemáti-
camente como1DdP (t)dt j Pt = p
E= K (p).
1Estamos usando la misma aproximación heurística para dP (t)dt
que en el modelo de colisión binaria no-relativista (ver pág. 48 de esta tesis). También se supone que
R t0h2D (Ps)i ds < 1 8t � 0, para tenerDR t
0[2D (Ps)]
1=2 dWs
E= 0 8t � 0.
120
El criterio anterior resulta plausible pues con él se logra captar información muy importante
del baño térmico.
Debido a que el baño térmico esta termalizado, la distribución estacionaria del momento de
la partícula Browniana en el modelo de colisión binaria, está dada por la función de Jüttner
�J (p) = (2MK1 (�M))�1 exp
h���p2 +M2
�1=2i 8p 2 R. (7.7)
Entonces para que la ecuación (7.6) cumpla el primer criterio de buena aproximación impuesto,
se pide que la ecuación (7.6) tenga como densidad estacionaria a la función (7.7). De las rela-
ciones de Einstein (4.8) se sigue que esto es equivalente a que las funciones � y D satisfagan
D (p) = ��1� (p)E (p) ; (7.8)
para cualquier p 2 R, donde E (p) =�m2 + p2
�1=2.Además, para que la ecuación (7.6) cumpla ambos criterios de buena aproximación, también
se pide que tal ecuación cumpla la relación
�dP (t)
dtj Pt = p
�= K (p) =
��P (t)
�tj Pt = p
�: (7.9)
Por las propiedades de la integral backward-Itô, se sabe que
�dP (t)
dtj Pt = p
�= �� (p) p+ d
dpD (p) :
Sustituyendo la relación de Einstein (7.8) en la ecuación anterior, se llega a que
�dP (t)
dtj Pt = p
�= �� (p) p+ ��1 d
dp[� (p)E (p)] :
Así que la condición (7.9) es equivalente a la ecuación diferencial
�� (p) p+ ��1 ddp[� (p)E (p)] = K (p) :
Cuando no se conoce una expresión para la función K (p), se puede, por ejemplo, tratar de
121
�jar una expresión analítica sencilla y, subsecuentemente, usar esta aproximación en la ecuación
anterior. Para esto existen diversas técnicas conocidas como aproximar por splines, entre otros.
122
Apéndice A
Procesos de Markov y Procesos de
Difusión
En esta sección se presenta un breve resumen de procesos de Markov y procesos de difusión.
Para mayor información sobre estos temas se recomienda consultar, por ejemplo, los libros
[6, 29, 50, 53, 56].
De�nición 101 Sean t0 � 0 y t0 < T <1. Sea fXt; t0 � t � Tg un proceso estocástico en Rd
de�nido sobre el espacio de probabilidad (;U ;P). El proceso Xt es un proceso de Markov si se
satisface la propiedad: para t0 � s � t � T y para todo B 2 B�Rd�, la ecuación
P (Xt 2 B j � (Xr; t0 � r � s)) = P (Xt 2 B j Xs)
se cumple con probabilidad 1.
Dado un proceso de Markov, el pasado y el futuro son estocásticamente independientes cuan-
do se conoce el presente. Este principio puede ser considerado como el principio de causalidad
de la física clásica llevado a sistemas dinámicos estocásticos.
Un concepto básico en los procesos de Markov lo constituye la probabilidad de transición.
De�nición 102 Sea fXt; t0 � t � Tg un proceso de Markov. Una función P (s;x; t; B) es lla-
mada una probabilidad de transición si tiene las siguientes propiedades:
123
(i) Para s � t y B 2 B�Rd��jos, se cumple con probabilidad 1
P (s;Xs; t; B) = P (Xt 2 B j Xs) :
(ii) P (s;x; t; �) es una medida de probabilidad sobre B�Rd�, para s � t y x 2 Rd �jos.
(iii) P (s; �; t; B) es B�Rd�-medible, para s � t y B 2 B
�Rd��jos.
(iv) Para t0 � s � u � t � T y B 2 B�Rd�y para todo x 2 Rd con la posible excepción de
un conjunto N � Rd tal que P (Xs 2 N) = 0, tenemos la llamada ecuación de Chapman-
Kolmogorov
P (s;x; t; B) =
ZRdP (u;y; t; B)P (s;x; u; dy) :
Si P (s; x; t; B) es una probabilidad de transición, la función P (s; x; t; B) también recibe el
nombre de probabilidad de transición del proceso de Markov Xt.
Para s; t 2 [t0; T ] �jos tales que s � t, P (s;x; t; B) esta únicamente de�nida como función
de x y B con la posible excepción de un conjunto N de valores de x (independiente de B) tal
que P (Xs 2 N) = 0. Además, en ocasiones se usa la notación
P (s;x; t; B) = P (Xt 2 B j Xs = x) ;
la cual es la probabilidad de que el proceso observado estará en el conjunto B al tiempo t si al
tiempo s, donde s � t, esta en el estado x.
Las probabilidades de transición, como toda distribución de una variable aleatoria, pueden
tener asociada una función de densidad. En el caso de los procesos de Markov, a la función de
densidad se le conoce como densidad de transición.
De�nición 103 Sea fXt; t0 � t � Tg un proceso de Markov y P (s;x; t; B) su probabilidad de
transición. Si existe una función p (s;x; t;y) que satisface
P (s;x; t; B) =
ZBp (s;x; t;y) dy
124
para cualesquiera t0 � s � t � T , x 2 Rd y B 2 B�Rd�, tal función recibe el nombre de
densidad de transición del proceso Xt.
Un tipo especial e importante de procesos de Markov son los procesos homogéneos.
De�nición 104 Un proceso de Markov fXt; t0 � t � Tg se dice homogéneo si su probabilidad
de transición P (s;x; t; B) es estacionaria, es decir si la condición
P (s+ u;x; t+ u;B) = P (s;x; t; B)
se cumple para cualesquiera t0 � s � t � T y t0 � s+ u � t+ u � T .
Cuando se tiene un proceso de Markov homogéneo, se acostumbra escribir
P (t� s;x; B) = P (s;x; t; B) ;
para 0 � t� t0 � T � t0.
Observación 105 ¿Cuándo un proceso de Markov fXt; t0 � t � Tg es también un proceso
estacionario? Una condición necesaria y su�ciente para ser estacionario es:
(i) Xt es homogéneo;
(ii) existe una distribución invariante P 0 en el espacio de estados, es decir
P 0 (B) =
ZRdP (t;x;B)P 0 (dx) para cualesquiera B 2 B
�Rd�, t 2 [0; T � t0] :
Si elegimos esta P 0 como distribución inicial para Xt0, entonces Xt es un proceso esta-
cionario. Además, si existe un invariante P 0, tenemos que para distribuciones iniciales arbi-
trarias y T =1;
l��mt!1
P [Xt 2 B] = P 0 (B)
para todo B 2 B�Rd�cuya frontera tiene P 0-medida 0; esto es, la distribución invariante
es una distribución estacionaria límite y es independiente de la distribución inicial. Existen
125
algunas condiciones bajo las cuales una función de transición homogénea P (t;x; B) tiene una
distribución invariante (véase [50]).
Los procesos de difusión son casos especiales de procesos de Markov con trayectorias mues-
trales continuas. Más precisamente, estos procesos quedan caracterizados bajo la siguiente
de�nición.
De�nición 106 Un proceso de Markov fXt; t0 � t � Tg en Rd con trayectorias muestrales
continuas a.s. es llamado un proceso de difusión si sus probabilidades de transición P (s;x; t; B)
satisfacen, para cualesquiera s 2 [t0; T ], x 2 Rd y " > 0, las siguientes tres condiciones :
(i) l��mt#s 1t�sRjx�yj>" P (s;x; t; dy) = 0;
(ii) existe una función f (s;x) que toma valores en Rd y
l��mt#s
1
t� s
Zjx�yj�"
(x� y)P (s;x; t; dy) = f (s;x) ;
(iii) existe una función B (s;x) que toma valores en Matd�d (R) tal que
l��mt#s
1
t� s
Zjx�yj�"
(x� y) (x� y)T P (s;x; t; dy) = B (s;x) :
Las funciones f y B son llamados los coe�cientes del proceso de difusión. En particular,
f es denominada vector de deriva y B es llamada matriz de difusión. B (s;x) es simétrica y
no-negativa de�nida.
El siguiente resultado establece una ecuación de evolución para la densidad de transición de
un proceso de difusión.
Teorema 107 Sea fXt; t0 � t � Tg un proceso de difusión en Rd; el cual posee una densidad
de transición p (s;x; t;y) y cumple uniformemente en s y x las relaciones límite de la de�nición
106. Si las derivadas @p@t ,@@yi
�f i (t;y) p
�y @2
@yi@yj
�Bij (t;y) p
�existen y son continuas, entonces,
para s y x �jos con s � t, la densidad de transición p (s;x; t;y) es una solución fundamental
de la ecuación de Fokker-Planck
@p
@t+
dXi=1
@
@yi�f i (t;y) p
�� 12
dXi;j=1
@2
@yi@yj�Bij (t;x) p
�= 0: (A.1)
126
Observación 108 Si se de�ne la distribución de Xt0 en términos de la probabilidad inicial
Pt0, la densidad de probabilidad p (t;y) del proceso Xt puede obtenerse por
p (t;y) =
ZRdp (s;x; t0;y)Pt0 (dx) :
Si en la ecuación de Fokker-Planck anterior integramos respecto a Pt0 (dx), obtenemos que
p (t;y) también satisface tal ecuación de Fokker-Planck.
La ecuación de Fokker-Planck también se conoce como la ecuación forward de Kolmogorov.
Algunas condiciones para la existencia de las densidades p (s;x; t;y) y p (t;y) se derivan de
las condiciones para que la ecuación (A.1) tenga solución. En general, la existencia de estas
densidades es un problema no trivial. El siguiente Teorema proporciona algunas condiciones
para la existencia de la densidad de transición de un proceso de difusión unidimensional.
Teorema 109 Sea fXt; 0 � t � Tg un proceso de difusión en R con deriva f (t; x) y coe�ciente
de difusión �2 (t; x) > 0, el cual satisface
(i) @@xf (t; x) y
@@x� (t; x) existen y son acotadas; las derivadas
@2
@x2� (t; x), @2
@x@t� (t; x),@2
@t2� (t; x)
y @@t� (t; x) existen;
(ii) g (t; x) esta de�nida por x =R g(t;x)0
dy�(t;y) ;
a (t; x) =
Z g(t;x)
0
@@t� (t; x)
�2 (t; y)dy +
f (t; g (t; x))
� (t; g (t; x))� 12
@
@x� (t; g (t; x)) ;
B (t; x) = �12a2 (t; x)��1
2
@
@xa (t; x)�
Z x
0
@
@ta (t; y) dy;
y B (t; x) cumple que
l��mjxj!1
1
1 + x2sup0�t�T
B (t; x) � 0:
Bajo estos supuestos la densidad transición del proceso existe.
Las condiciones del Teorema anterior incluso permiten encontrar una expresión explícita
para la densidad de transición (véase, [29] págs. 96-97).
127
Apéndice B
Distribución de Maxwell
La distribución de Maxwell modela a un gas (cuasi-ideal) diluido, cerca del equilibrio tér-
mico, con efectos cuánticos insigni�cantes y con velocidades no relativistas. Tal distribución
establece que la velocidad de una partícula de un gas con las propiedades anteriores obedece
una densidad M : R! R+ dada por
M (v;m;�; d) =
��m
2�
�d=2exp
���mv22
�; (B.1a)
donde m es la masa en reposo de una partícula del gas, v 2 Rd es la velocidad no-relativista y
T = (kB�)�1 es la temperatura. A partir de la ecuación anterior, se encuentra que la densidad
que gobierna el momento de una partícula del gas es
�M (p;m;�; d) =
��
2�m
�d=2exp
���mp22
�; (B.1b)
siendo p =mv 2 Rd el momento no-relativista y �M : R! R+.
Obviamente la densidad Maxwelliana para la velocidad no respeta el postulado relativista
fundamental de que nada puede viajar más rápido que la luz. La función de densidad adecuada
para trabajar en el caso relativista se expone en el Apéndice D.
Cuando se trabaje en el caso clásico, se asume que las partículas de cualquier baño en
equilibrio térmico obedecen una densidad de Maxwell de la forma (B.1). Lo mismo aplica para
una partícula Browniana que este en equilibrio térmico con el baño circundante.
129
Apéndice C
Conceptos Básicos de Relatividad
El objetivo de este apéndice es establecer la notación, las convenciones y los conceptos de
relatividad especial que se utilizan en el tratamiento del movimiento Browniano relativista.
Para profundizar en estos temas se recomienda consultar [37, 26, 54, 58].
C.1. Notación y convenciones
En relatividad especial, un marco de referencia inercial1 � corresponde a un sistema coorde-
nado Cartesiano global del espacio-tiempo. Un evento E del espacio-tiempo esta etiquetado por
un vector (1 + d)-dimensional x = (x�) = (ct;x) =�ct; x1; : : : ; xd
�en �, donde d es el número
de dimensiones espaciales y, adoptando unidades naturales, c = 1 es la velocidad de la luz. Los
super y subíndices griegos �; �; : : : toman los valores 0; 1; : : : ; d; mientras que se usan índices
latinos i; k; : : : 2 f1; : : : ; dg para las componentes espaciales. Los vectores con superíndices son
llamados contravariantes.
Se usa el convenio de sumación de Einstein, el cual permite simpli�car sumas como
y = cixi :=
dXi=1
cixi,
1Entendemos por marco de referencia inercial a un marco de referencia no acelerado.
131
para los índices latinos, o bien
y = c�x� :=dX�=0
c�x�,
para los índices griegos. Además, se usa la notación
x2 := jxj2 =�xi�2
para denotar la norma euclidiana del vector de componentes espaciales x =�x1; : : : ; xd
�.
Con respecto al marco de coordenadas Cartesiano �, las componentes ��� del tensor métrico
del espacio-tiempo de Minkowski se de�nen por
��� =
8>>><>>>:�1 � = � = 0
+1 � = � = 1; : : : ; d
0 � 6= �:
Por de�nición, las componentes de los vectores covariantes (x�) se obtienen al contraer el
vector contravariante (x�) con ��� , es decir
x� := ���x� ) (x�) = (�t;x) :
Los vectores (x�) y (x�) son llamados cuadrivectores, sin importar el número de dimensiones.
La distancia en el espacio-tiempo de Minkowski entre dos eventos xA = (x�A) = (tA;xA) y
xB = (x�B) = (tB;xB) se de�ne por
d (xA; xB)2 := ��� (x�A � x�B) (x�A � x�B) = � (tA � tB)
2 + (xA � xB)2 .
Por de�nición, la separación entre dos eventos es
Tipo tiempo, si d (xA; xB)2 < 0;
Tipo luminoso, si d (xA; xB)2 = 0;
Tipo espacio, si d (xA; xB)2 > 0:
Los eventos con separación tipo tiempo pueden estar causalmente conectados por (una serie
132
de) señales viajando a una velocidad menor o igual a la velocidad de la luz. Los eventos con
separación tipo luminosa pueden estar relacionados solo por señales ininterrumpidas que viajan
a la velocidad de la luz. Los eventos con separación tipo espacio son causalmente disconexos.
El movimiento clásico de una partícula masiva puntual a través del espacio tiempo, corres-
ponde a una curva en � lo su�cientemente suave del tipo tiempo. Tal curva recibe el nombre
de línea de universo. Considérese un observador estacionario O, el cual se encuentra en reposo
en �. Es natural que O parametrice el movimiento de la partícula usando el tiempo t de �, es
decir que O describe la línea de universo como�x0 (t) ; xi (t)
�con x0 (t) = t. En una vecindad
de cada punto (evento) sobre la línea de universo de la partícula, una diferencial del tiempo
propio in�nitesimal puede de�nirse por
d� := �����dx
�dx��1=2
=�dt2 � dx2
�1=2= dt
�1� v2
�1=2; (C.1)
donde v (t) = dx(t)dt es la velocidad de la partícula en �. De acuerdo a los postulados de la
relatividad especial, d� es el intervalo de tiempo medido por un reloj intrínseco comoviéndose
con la partícula, mientras que dt es un intervalo de tiempo medido por un reloj en reposo en �.
La cuadrivelocidad (u�) de una partícula masiva se de�ne como la derivada de la línea de
universo con respecto a su tiempo propio, es decir
u� :=dx�
d�) u�u
� = �1:
Para una partícula puntual con masa en reposo m > 0, el cuadrivector de energía-momento
(p�) =�p0; p1; : : : ; pd
�= (�;p) se de�ne por
p� := mu� ) p�p� = �m2: (C.2)
Al comparar con (C.1), se tiene que para una partícula con velocidad v en �
p0 = � = m (v) ; p = m (v)v; (v) =�1� v2
��1=2:
La entrada p0 = � es la energía relativista de la partícula.
133
C.2. Transformaciones de Lorentz-Poncairé
En relatividad especial, las transformaciones lineales a�nes de Lorentz-Poncairé de la forma
x0 = �x+ a () x0� = ���x� + a� (C.3)
describen la transición de un marco de referencia inercial � a otro marco inercial �0. Se usa
la notación (x0�) para los eventos en �0. El cuadrivector constante (a�) traslada el origen del
tiempo y el espacio, mientras que la matriz constante de Lorentz����
�puede dar cuenta de
una rotación, un cambio de orientación o una velocidad relativa entre los dos marcos � y �0.
Las componentes matriciales ��� están determinadas por la condición
d�x0A; x
0B
�2= d (xA; xB)
2 , �� ������ = � �; (C.4)
lo cual signi�ca que las relaciones causales se preservan durante transiciones entre sistemas
inerciales. Las transformaciones de Lorentz-Poncairé forman un grupo. De particular interés
para nuestro tratamiento es el subgrupo de las transformaciones propias de Lorentz, de�nido
al tomar a� = 0 y las restricciones adicionales
�00 � 1; det����
�= +1:
Los requisitos anteriores excluyen reversión temporal e inversión espacial. Ejemplos de trans-
formaciones propias de Lorentz son las rotaciones
�00 = 1; �i0 = �0i = 0; �ij = Rij ; (C.5)
donde�Rij�es una matriz de rotación (i.e., det
�Rij�= 1 y RijRkj = �ij con �ij una delta de
Kronecker). Otro ejemplo lo constituyen los boosts de Lorentz
�00 = ; �i0 = �0i = � wi; �ij = �ij +wiwj
w2( � 1) ; (C.6)
134
con velocidadw =�w1; : : : ; wd
�y factor de Lorentz :=
�1�w2
��1=2. Para ilustrar brevementeel efecto de un boost, considérese una partícula en reposo en el origen de � y, por lo tanto,
descrita por la línea de universo (x� (t)) � (t;0) en �. Aplicando el boost de Lorentz (C.6) a
(x�) = (t;0), se encuentra que
x00 = �00x0 = t = t0, x0i = �i0x0 = � wit = �wit0;
lo cual signi�ca que la partícula viaja a velocidad constate w0= �w a través de �0. La inversa
de la transformación (C.6) se obtiene al reemplazar w con �w.
Toda transformación propia de Lorentz puede escribirse como el producto de un boost con
una rotación. Más precisamente, si � es un transformación propia de Lorentz, existe un boost
�bst de la forma (C.6) y una rotación �rot de la forma (C.5) tales que
� = �bst�rot: (C.7)
De las ecuaciones (C.2), (C.3) y (C.4), se obtiene la ley de transformación de energía-
momento relativista
p0� = ���p�: (C.8)
Combinando la ecuación anterior y la igualdad (C.4), es posible veri�car la condición de �mass-
shell�
m2 = �2 � p2 = �02 � p02 = m02;
lo cual implica que la masa en reposo es un invariante de Lorentz. En particular, de la condición
anterior se sigue que las componentes espaciales de la ecuación (C.8) pueden escribirse como
p0i (p) = �i0�m2 + p2
�1=2+ �ijpj :
135
Apéndice D
Densidades de Probabilidad en
Relatividad Especial
En este apéndice se presentan algunos resultados para las distribuciones de probabilidad
en el espacio fase de una partícula. Este apartado esta dividido en dos secciones, la primera
de ellas trata la invarianza de Lorentz de las densidades en el espacio fase y la segunda habla
sobre la densidad de Jüttner. El primer tema se trata con mayor profundidad en las referencias
[26, 39, 40, 41], mientras que el segundo puede ser consultado en [26].
D.1. Invarianza de Lorentz de las densidades en el espacio fase
Considérense dos marcos de referencia inercial � y �0, de tal manera que �0 se mueve con
velocidad constante w 6= 0 relativa a �. También, considérense observadores O y O0 tales que
O se encuentra en estado de reposo en � y O0 se encuentra en reposo en el marco �0. Dada
un partícula en movimiento, los observadores medirán distintas densidades para la posición en
el espacio fase de la partícula. Denótese por f (t;x;p) y f 0 (t0;x0;p0) la densidades medidas
por O y O0, respectivamente. Naturalmente, surge la pregunta de como están relacionadas
f (t;x;p) y f 0 (t0;x0;p0). En la teoría no relativista, el cambio de un sistema inercial a otro no
afecta la coordenada temporal; de ahí que, en este caso, se puedan usar las leyes usuales de
transformación entre densidades. En contraste, la situación se vuelve más complicada en la teoría
relativista, pues la de�nición de f (t;x;p) y f 0 (t0;x0;p0) se basa en la noción de simultaneidad:
137
las mediciones de O y O0 se re�eren a diferentes hiperplanos �t = const�y �t0 = const�en el
espacio de Minkowski, respectivamente.
En un artículo publicado en 1969, van Kampen [40] probó que la densidad f (t;x;p) de la
posición de una partícula en el espacio fase se transforma como un escalar de Lorentz, es decir
f (t;x;p) = f 0�t0;x0;p0
�;
donde (t;x;p) y (t0;x0;p0) están conectados por transformaciones de Lorentz. Además, van
Kampen demostró que
Zf�t0;x0;p0
�ddx0ddp0 =
Zf (t;x;p) ddxddp;
lo que implica que la función f 0 satisface la condición de normalización t0-simultánea
Zf�t0;x0;p0
�ddx0ddp0 = 1:
En la referencia [26] se pueden encontrar varias consecuencias de esto.
D.2. La distribución de Jüttner
En el estudio del movimiento Browniano relativista, estamos especialmente interesados en
distribuciones que se cumplen en equilibrio térmico. Cuando se postulan ecuaciones relativistas
de Langevin, estas distribuciones deben conocerse de antemano para poder especi�car la relación
de �uctuación-disipación correcta. También las distribuciones del equilibrio térmico se usan para
derivar ecuaciones del tipo Langevin a partir de modelos microscópicos.
Al comienzo del siglo pasado, era comúnmente aceptado que un gas (cuasi-ideal) diluido en
equilibrio térmico estaba descrito por la densidad Maxwelliana para la velocidad M : R! R+
con
M (v;m;�; d) =
��m
2�
�d=2exp
���mv22
�; (D.1a)
138
o, equivalentemente, por la densidad para el momento �M : R! R+ con
�M (p;m;�; d) =
��
2�m
�d=2exp
���mp22
�; (D.1b)
dondem es la masa en reposo de una partícula del gas, v = p=m 2 Rd la velocidad no-relativista
y T = (kB�)�1 la temperatura. Después de que Einstein formuló su teoría de la relatividad
especial en 1905, Planck y otros inmediatamente se dieron cuenta de que la distribución (D.1)
estaba en con�icto con el postulado relativista fundamental de que nada puede ir más rápido
que la luz. Una solución para este problema fue encontrada por Jüttner en 1911, quien propuso
reemplazar la densidad de Maxwell por
J (v;m;�; d) = Z�1d md (v)2+d exp (��m (v))� (1� jvj) ; (D.2a)
donde � : R! R es una función Heaviside y J : (�1; 1)! R+ solo puede evaluarse en los
valores permitidos para la velocidad. Para el momento relativista p = m (v)v, la distribución
anterior da lugar a la densidad
�J (p;m;�; d) = Z�1d exph���m2 + p2
�1=2i; (D.2b)
donde �J : R! R+. De manera similar a la distribución de Maxwell (D.1), la distribución de
Jüttner (D.2) se re�ere al marco de referencia del laboratorio � donde el recipiente que contiene
al gas está en reposo. Para las dimensiones espaciales d = 1; 2; 3, la constante de normalización
Zd (m;�) =ZRexp
h���m2 + p2
�1=2iddp
puede expresarse como
Z1 = 2mK1 (�m) ;
Z2 = 2�m2 exp (��m) (1 + �m) (�m)�2 ;
Z3 = 4�m3K2 (�m) (�m)�1 ;
con Kn (z) denotando funciones de Bessel modi�cadas de segundo tipo. La energía promedio
139
por partícula se obtiene de la diferenciación
h�id = �@
@�logZd;
dando lugar a
h�i1 = mK0 (�m) +K2 (�m)
2K1 (�m);
h�i2 =2
�+
m2�
1 +m�;
h�i3 =3
�+m
K1 (�m)
K2 (�m);
y exhibiendo el comportamiento límite
l��m�!1
h�id = m;
l��m�!0
� h�id = d:
Además, se puede demostrar que, para dimensiones espaciales arbitrarias d, el valor medio
hp � vi es independiente de la masa; más precisamente
hp � vi =�p2
�
�=d
�;
lo cual permite considerar a hp � vi como un termómetro estadístico.
De particular relevancia en termodinámica del equilibrio son los sistemas con�nados que
pueden ser descritos por una densidad ' (x;p) estacionaria e isotrópica en un marco de re-
ferencia especi�co �. Un ejemplo típico lo constituye un gas en equilibrio enclaustrado en un
contenedor en reposo en el marco del laboratorio �. En la teoría estándar del movimiento
Browniano tales sistemas a menudo juegan el rol del baño térmico. En el marco �, un gas
espacialmente homogéneo está descrito por una función de densidad de la forma
' (x;p) = V �1J (x;V)�J (p) ; (D.3)
donde V es el volumen en reposo (medida de Lebesgue) de la región contenedor V � Rd (V
140
abierto y conexo) en � y J (x;V) denota la función indicadora
J (x;V) =
8<: 1; si x 2V;
0; si x 62V:
Más especí�camente, estaremos interesados en contenedores cúbicos de la forma V =dYi=1
[L=2; L=2],
donde L > 0.
En general, cuando un baño térmico relativista se encuentra en equilibrio térmico, se asume
que sus partículas cumplen una distribución de Jüttner (en el marco del laboratorio �) de la
forma (D.3) cuando se especi�que un volumen �nito V; o bien, de la forma (D.2) cuando solo
importe el momento. Lo mismo aplica para una partícula Browniana relativista cuando ésta se
encuentra en equilibrio térmico con el baño.
141
Apéndice E
Notación
Conviene establecer algo de notación con el �n de facilitar la exposición de este trabajo.
Notación 110 A continuación se presenta una lista de símbolos y convenciones:
1. Rd es el espacio euclidiano d-dimensional con la topología usual.
2. B�Rd�es la �-álgebra de Borel en Rd.
3.�xi�denota al vector x =
�x1; : : : ; xd
�, donde d tiene que especi�carse.
4. A =�Aij�es la matriz con entradas Aij, la dimensión a de ser especi�cada.
5. AT es la transpuesta de la matriz A (no excluye vectores).
6. Matd�m (R) es el conjunto de las matrices con coe�cientes reales y dimensión d�m.
7. trA es la traza de la matriz A.
8. Se toma como norma de una matriz A de dimensión d�m a jAj2 =Pdi=1
Pmj=1
�Aij�2=
tr�AAT
�9. �ij denota a un delta de Kronecker.
10. � (X) es la �-álgebra generada por X y � (Xt; t 2 B) es la �-álgebra más pequeña que
contiene a � (Xt) para todo t 2 B.
11. hXi denota la esperanza de la variable aleatoria X:
143
12. st-l��m es el límite estocástico, ac-l��m es el límite con probabilidad 1 y qm-l��m el límite
cuadrático medio.
13. kB es la constante de Boltzmann.
14. c = 1 es la velocidad de la luz.
144
Bibliografía
[1] M. Abdel-Aziz y S. Gavin. Causal di¤usion and the survival of charge �uctuations in
nuclear collisions. Phys. Rev. C, 70(3):034905, 2004.
[2] M. Abdel-Aziz y S. Gavin. Causal di¤usion at RHIC. J. Phys. G: Nucl. Phys., 31:S77�S84,
2005.
[3] L. J. S. Allen. Stochastic Processes with Applications to Biology. Prentice-Hall, Upper
Saddle River, NJ, 2002.
[4] R. Aloisio, V. Berezinsky, y A. Gazizov. Superluminal problem in di¤usion of relativistic
particles and its phenomenological solution. arXiv/0805.1867, 2008.
[5] J. Angst y J. Franchi. Central limit theorem for a class of relativistic di¤usions. J. Math.
Phys., 48(8):083101, 2007.
[6] L. Arnold. Stochastic Di¤erential Equations: Theory and Applications. John Wiley & Sons,
1974.
[7] C. Barbachoux, F. Debbasch y J. P. Rivet. The spatially one-dimensional relativistic
Ornstein-Uhlenbeck process in an arbitrary in-ertial frame. Eur. Phys. J. B, 19:37�47,
2001.
[8] E. Bertschinger. Brownian Motion of Stars, Dust, and Invisible Matter. AIP Conf. Proc.,
861(1):97�105, 2006.
[9] J. P. Bouchaud y M. Potters. Theory of Financial Risks: From Statistical Physics to Risk
Management. Cambridge University Press, Cambridge, 2001.
145
[10] F. Debbasch, K. Mallick, y J. P. Rivet. Relativistic Ornstein-Uhlenbeck Process. J. Stat.
Phys, 88:945�966, 1997.
[11] F. Debbasch y J. P. Rivet. A di¤usion equation from the relativistic Ornstein-Uhlenbeck
process. J. Stat. Phys, 90:1179�1199, 1998.
[12] C. D. Dermer. The production spectrum of a relativistic Maxwell-Boltzmann gas. Astro-
phys. J., 280:328�333, 1984.
[13] M. E. Dieckmann, L. O�C. Drury, y P. K. Shukla. On the ultrarelativistic two-stream
instability, electrostatic turbulence and Brownian motion. New J. Phys., 8:40, 2006.
[14] R. M. Dudley. Lorentz-invariant Markov processes in relativistic phase space. Ark. Mat.
Astron. Fys., 6:241�268, 1965.
[15] R. M. Dudley. A note on Lorentz-invariant Markov processes. Ark. Mat. Astron. Fys.,
6:575�581, 1967.
[16] R. M. Dudley. Asymptotics of Some Relativistic Markov Processes. Proc. Nat. Acad. Sci.
USA, 70(12):3551�3555, 1973.
[17] R. M. Dudley. Recession of some relativistic Markov processes. Rocky Mountain J. Math.,
4:401�406, 1967.
[18] J. Dunkel y P. Hänggi. Theory of the relativistic Brownian motion: The (1+1)-dimensional
case. Phys. Rev. E, 71:016124, 2005.
[19] J. Dunkel y P. Hänggi. Theory of the relativistic Brownian motion: The (1+3)-dimensional
case. Phys. Rev. E, 72:036106, 2005.
[20] J. Dunkel, S. Hilbert y P. Hänggi. Langevin-Gleichungen mit nichtlinearer Reibung. In T.
Pöschel, H. Malchow, and L. Schimansky-Geier, editors, Irreversible Prozesse und Selb-
storganisation, pages 11�21. Logos-Verlag, Berlin, 2006. ISBN 3-8325-1350-7.
[21] J. Dunkel y P. Hänggi. Relativistic Brownian motion: From a microscopic binary collision
model to the Langevin equation. Phys. Rev. E, 74:051106, 2006. Erratum, Phys. Rev. E,
74:069902(E), 2006.
146
[22] J. Dunkel y P. Hänggi. One-dimensional nonrelativistic and relativistic Brownian motions:
A microscopic collision model. Physica A, 374(2):559�572, 2007.
[23] J. Dunkel. Relativistic Brownian Motion and Di¤usion Processes. Ph.D. thesis, Universität
Augsburg, 2008.
[24] J. Dunkel, P. Talkner y P. Hänggi. Relativistic di¤usion processes and random walk models.
Phys. Rev. D, 75:043001, 2007.
[25] J. Dunkel, P. Hänggi y S. Weber. Time parameters and Lorentz transformations of rela-
tivistic stochastic processes. arXiv:0812.0466v1, 2008.
[26] J. Dunkel y P. Hänggi. Relativistic Brownian Motion. Phys. Rep., 471(1):1-73, 2009.
[27] R. J. Elliott y P. E. Kopp. Mathematics of Financial Markets. Springer-Verlag, New York,
1999.
[28] L. C. Evans. An Introduction to Stochastic Di¤erential Equations.
http://math.berkeley.edu/~evans, 2009.
[29] I. I. Gihman y A. V. Skorohod. Stochastic Di¤erential Equations. Springer-Verlag, Berlín,
Heidelberg, New York, 1972.
[30] P. Glasserman.Monte Carlo Methods in Financial Engineering. Number 53 in Applications
of Mathematics. Springer, New York, 2004.
[31] R. Hakim. A covariant theory of relativistic Brownian Motion. I. Local Equilibrium. J.
Math. Phys., 6(10):1482�1495, 1965.
[32] R. Hakim. Relativistic stochastic processes. J. Math. Phys., 9:1805�1818, 1968.
[33] R. Hakim. Remarks on Relativistic Statistical Mechanics. I. J. Math. Phys., 8(6):1315�
1344, 1965.
[34] R. Hakim. Remarks on Relativistic Statistical Mechanics. II. Hierarchies for the Reduced
Densities. J. Math. Phys., 8(7):1397�1400, 1965.
147
[35] H. van Hees y R. Rapp. Thermalization of heavy quarks in the quark-gluon plasma. Phys.
Rev. C, 71(3):034907, 2005.
[36] H. van Hees, V. Greco, y R. Rapp. Heavy-quark probes of the quark-gluon plasma and
interpretation of recent data taken at the BNL Relativistic Heavy Ion Collider. Phys. Rev.
C, 73:034913, 2006.
[37] O. D. Johns. Analytical Mechanics for Relativity and Quantum Mechanics. Oxford Univer-
sity Press, 2005.
[38] F. Jüttner. Das Maxwellsche Gesetz der Geschwindigkeitsverteilung in der Relativtheorie.
Ann. Phys. (Leipzig), 34(5):856�882, 1911.
[39] N. G. van Kampen. Relativistic thermodynamics of moving systems. Phys. Rev., 173:295�
301, 1968.
[40] N. G. van Kampen. Lorentz-invariance of the distribution in phase space. Physica, 43:244�
262, 1969.
[41] N. G. van Kampen. Relativistic thermodynamics. J. Phys. Soc. Jap. Suppl., 26:316�321,
1969.
[42] I. Karatzas y S. E. Shreve. Brownian Motion and Stochastic Calculus. Springer-Verlag,
Berlín, Heidelberg, New York, 1988.
[43] N. E. Karoui y L. Mazliak. Backward Stochastic Di¤erential Equations. Addison-Wesley,
1997.
[44] P. E. Kloeden y E. Platen. Numerical Solution of Stochastic Di¤erential Equations.
Springer-Verlag, Berlín, Heidelberg, New York, 1995.
[45] M. von Laue. Zur Dynamik der Relativitätstheorie. Ann. Phys. (Leipzig), 35:524�542, 1911.
[46] H. P. McKean Jr. Stochastic Integrals. Academic Press, New York, 1969.
[47] K. von Mosengeil. Theorie der stationären Strahlung in einem gleichförmig bewegten
Hohlraum. Ann. Phys. (Leipzig), 22:867�904, 1907.
148
[48] B. Oksendal. Stochastic Di¤erential Equations: An Introduction with Applications.
Springer-Verlag, Berlín, Heidelberg, New York, Tokio, 1989.
[49] H. Ott. Lorentz-Transformation der Wärme und der Temperatur. Z. Phys., 175:70�104,
1963.
[50] Yu. V. Prohorov y Yu. A. Rozanov. Probability Theory. Springer-Verlag, New York, Hei-
delberg, Berlín, 1969.
[51] P. E. Protter. Stochastic Integration and Di¤erential Equations. Springer-Verlag, Berlín,
Heidelberg, New York, 2003.
[52] R. Rapp, V. Greco, y H. van Hees. Heavy-Quark Spectra at RHIC and Resonances in the
QGP. Nucl. Phys. A, 774:685�688, 2006.
[53] L. C. G. Rogers y D. Williams. Di¤usions, Markov Processes and Martingales, Vol. 2.
Cambridge University Press, 2004.
[54] R. U. Sexl y H. K. Urbantke. Relativity, Groups, Particles. Springer Physics. Springer,
Wien, 2001.
[55] B. Svetitsky. Di¤usion of charmed quarks in the quark-gluon plasma. Phys. Rev. D,
37(9):2484�2491, 1988.
[56] C. Tudor. Procesos Estocásticos. Sociedad Matemática Mexicana. Aportaciones Matemáti-
cas,Textos Nivel Avanzado 2. 2002.
[57] G. N. Watson. A Treatise on the Theory of Bessel Functions. Merchant Books, 2008.
[58] S. Weinberg. Gravitation and Cosmology. John Wiley & Sons, 1972.
[59] D. Wilkinson. Stochastic Modelling for Systems Biology, volumen 11 de Mathematical and
Computational Biology Series. Chapman & Hall/CRC, London, 2006.
[60] G.Wolschin, M. Biyajima, T. Mizoguchi, y N. Suzuki. Time evolution of relativistic d+au
and au+au collisions. Ann. Phys., 15(6):369�378, 2006.
149
[61] C. K. Yuen. Lorentz transformation of thermodynamic quantities. Am. J. Phys., 38:246�
252, 1970.
[62] B. Wolfe y F. Melia. Covariant kinetic theory with an application to the coma cluster.
Astrophys. J., 638, 2006.
[63] G. Wolschin. Di¤usion and local decon�nement in relativistic systems. Phys. Rev. C,
69(2):024906, February 2004.
150