La PizarraMatemática.

Máster y doctorado

Ecuaciones dispersivas

Un paquete de ondas puede ensancharse y bajar de altura sin perder su norma total. Eso es posible porque sus frecuencias acumulan fases distintas. Vamos a separar velocidad de una cresta, velocidad del paquete y conservación antes de estudiar la parte no lineal.

1. La relación de dispersión dice cómo evoluciona cada frecuencia

Una onda plana tiene forma e^{i(kx-\omega t)}. k es la frecuencia espacial y ω la temporal.

La ecuación relaciona ambas mediante ω=ω(k).
Una cresta conserva kx−ωt; si k≠0, se mueve a velocidad ω/k.

A eso llamamos velocidad de fase.

Un paquete mezcla frecuencias cercanas a k₀. Al desarrollar ω cerca de k₀, el término lineal desplaza su envolvente a velocidad ω′(k₀), llamada velocidad de grupo. En varias dimensiones es ∇ω. La curvatura de ω ayuda a que frecuencias cercanas se separen; oscilar no implica por sí solo ensancharse.

Un modelo que traslada sin ensanchar

ω(k)=ck con c constante.

  1. Fase

    Para k≠0, ω/k=c.
  2. Grupo

    ω′(k)=c en todas las frecuencias.
  3. Comparar

    Todos los componentes viajan a la misma velocidad: un perfil puede trasladarse sin cambiar de forma.

Resultado

La linealidad de ω produce transporte uniforme.

Lo comprobamos

La aparición de una fase oscilatoria no basta para concluir dispersión.

Dos frecuencias con distinta velocidad de grupo

ω(k)=k², cerca de k=1 y cerca de k=2.

  1. Derivar

    ω′(k)=2k.
  2. Valores

    Las velocidades de grupo son dos y cuatro.

  3. Separación

    Los centros de paquetes estrechos en frecuencia pueden alejarse uno del otro.

Resultado

La velocidad depende de la frecuencia.

Lo comprobamos

El movimiento exacto de un paquete requiere también sus amplitudes y su anchura frecuencial.

2. En Schrödinger libre el signo fija la dirección de la fase

Fijamos i uₜ+uₓₓ=0.
Sustituir la onda plana da i uₜ=ωu y uₓₓ=−k²u, por lo que ω=k².
La velocidad de fase es k si k≠0 y la de grupo 2k.

Una frecuencia negativa tiene velocidad de grupo negativa, aunque ω sea no negativa.

La velocidad de grupo es la pendiente de la relación de dispersión
ω(k)=k²ω=k²k<0k>0En k=1: ω/k=1Pendiente ω′(1)=2Cresta y envolvente pueden ir distintas.
Gráfica de omega igual a k al cuadrado. La recta verde es la tangente en k igual a uno, con ecuación omega igual a dos k menos uno; se dibuja solo donde queda dentro del cuadro. La velocidad de fase en ese punto es omega sobre k igual a uno. La de grupo es la pendiente dos. El eje vertical corresponde a omega cero en la horizontal de origen; para k negativo la pendiente es negativa aunque omega sea no negativa. Se compara una cresta de onda con una envolvente de frecuencias cercanas.

Con transformada espacial de fase e elevado a −ixξ, la ecuación en Rⁿ es i∂ₜû−|ξ|²û=0.

Se resuelve como \widehat u(t,\xi)=e^{-it|\xi|^2}\widehat u_0(\xi).
Denotamos el propagador por U(t)=eⁱᵗΔ.

Una onda plana en toda Rⁿ no es L²; sirve para comprobar el símbolo como función suave o distribución.

Una frecuencia positiva

k=3.

  1. Frecuencia temporal

    ω=3²=9.
  2. Cresta

    v de fase=ω/k=9/3=3.
  3. Paquete

    v de grupo=ω′(3)=2·3=6.
Resultado

Frecuencia temporal, velocidad de fase y velocidad de grupo son tres magnitudes distintas.

Lo comprobamos

Las dos velocidades no se intercambian aunque salgan de la misma curva.

Una frecuencia negativa

k=−2.

  1. Frecuencia temporal

    ω=(−2)²=4.
  2. Fase

    ω/k=4/(−2)=−2.
  3. Grupo

    2k=−4.
Resultado

El grupo se mueve hacia valores decrecientes de x.

Lo comprobamos

La no negatividad de ω no obliga a que las velocidades sean positivas.

3. Dispersar no es disipar la norma L dos

El multiplicador e elevado a −it|ξ|² tiene módulo uno.

Plancherel da ||U(t)u₀||₂=||u₀||₂.

Se obtiene un grupo unitario, continuo en tiempo en L², con U(t+s)=U(t)U(s) y U(−t)=U(t)*. Para un dato H¹ se conserva también la norma del gradiente.

El calor tiene multiplicador e elevado a −t|ξ|²: amortigua las frecuencias altas y pierde norma L². En Schrödinger puede bajar la altura de la solución al distribuir su masa en una zona más amplia. La reversibilidad y la conservación de norma no implican que se mantenga la forma espacial.

Conservar una norma dada

||u₀||₂=5 en la evolución libre.

  1. Fourier

    Cada coeficiente se multiplica por un número de módulo uno.

  2. Plancherel

    La integral de su módulo al cuadrado no cambia.

  3. Norma

    Al tomar raíz, ||u(t)||₂=5 en todo tiempo.
Resultado

La norma se conserva.

Lo comprobamos

No hemos concluido que ||u(t)||∞ sea cinco ni que permanezca constante.

El contraste con el calor

Comparar una frecuencia ξ con |ξ|²=4 a tiempo t=1.

  1. Schrödinger

    El factor es e elevado a −4i y tiene módulo uno.

  2. Calor

    El factor es e⁻⁴, menor que uno.

  3. Interpretar

    El primero cambia fase; el segundo reduce amplitud de esa frecuencia.

Resultado

Cambiar fase y amortiguar son mecanismos diferentes.

Lo comprobamos

Una frecuencia aislada en la recta no es un dato L²; la comparación lee los multiplicadores.

4. Una gaussiana exacta permite medir altura, anchura y masa

Tomamos u₀(x)=exp(−x²/(2σ²)), σ>0.

La solución es u(x,t)=\bigl(1+2it/\sigma^2\bigr)^{-1/2}\exp\!\left(\frac{-x^2}{2\sigma^2(1+2it/\sigma^2)}\right), con la raíz continua que vale uno en t=0. Se puede comprobar derivando o calcular por Fourier gaussiana.

La densidad se ensancha sin perder masa L dos
Schrödinger libre: densidad |u|²u₀(x)=exp(−x²/2)t=0t=1Anchura: 1 → √5Altura: 1 → 1/√5∫ |u|² = √π en ambos tiempos.
Se representa el módulo al cuadrado de la solución gaussiana de i u t más u x x igual a cero con sigma igual a uno. En tiempo cero la densidad es exp de menos x al cuadrado. En tiempo uno es uno sobre raíz de cinco por exp de menos x al cuadrado sobre cinco. Las curvas comparten escalas; sus integrales en toda la recta son raíz de pi. La menor altura representa dispersión, no disipación de la norma L dos. El dibujo de densidad no representa la fase compleja.
Si D=√(1+4t²/σ⁴), su densidad es |u|²=D⁻¹exp[−x²/(σ²D²)].

Su anchura característica es σD, la altura de la densidad es D⁻¹ y su integral σ√π. Al multiplicar el dato por eⁱᵏ⁰ˣ, aparece además un desplazamiento del centro a 2k₀t, con su fase correspondiente.

Ensanchar sin cambiar el área

σ=1, comparar t=0 y t=1.

  1. Factor

    D pasa de uno a √5.

  2. Perfil

    La anchura pasa de uno a √5 y la altura de la densidad de uno a 1/√5.

  3. Integral

    El producto altura por anchura conserva el factor gaussiano: ∫|u|²=√π.

Resultado

La densidad se ensancha y baja sin perder masa L².

Lo comprobamos

La norma ||u||₂ es π¹ᐟ⁴; √π es su cuadrado.

Mover un paquete además de ensancharlo

Dato e elevado a 2ix por exp(−x²/2).

  1. Centro de frecuencia

    k₀=2.
  2. Grupo

    La velocidad del centro es 2k₀=4.
  3. Tiempo

    En t=1 su centro está en x=4; la anchura sigue creciendo como en la gaussiana sin modular.

Resultado

Traslación y ensanchamiento pueden ocurrir a la vez.

Lo comprobamos

La modulación tiene módulo uno, así que no cambia la densidad inicial.

5. Una cota puntual necesita integrabilidad adicional del dato

En Rⁿ y t≠0, el núcleo de U(t) es (4\pi it)^{-n/2}\exp\!\left(\frac{i|x-y|^2}{4t}\right), interpretado con la fase y regularización oscilatoria apropiadas. Su módulo es (4π|t|)⁻ⁿᐟ².

La desigualdad triangular da ||U(t)u₀||∞≤(4π|t|)⁻ⁿᐟ²||u₀||₁.
Se exige un dato L¹ para esta cota L¹→L∞.

El núcleo no es absolutamente integrable en toda la recta o espacio sin el dato. Interpolar con unitariedad produce otras cotas. En un toro las recurrencias y la geometría de frecuencias impiden copiar sin cambios este decaimiento global en tiempo.

Leer una tasa de decaimiento

En una dimensión y con ||u₀||₁=2.

  1. Cota

    ||u(t)||∞≤2/√(4π|t|).
  2. Dos tiempos

    Al multiplicar |t| por cuatro, el lado derecho se divide entre dos.

  3. Interpretar

    La cota mejora a tiempos grandes; no es una fórmula de igualdad para cada dato.

Resultado

La dimensión determina el exponente temporal.

Lo comprobamos

En t=0 la expresión diverge y no se usa como una cota inicial uniforme.

Un dato L dos puede no ser L uno

g(x)=(1+|x|)⁻³ᐟ⁴.

  1. Cuadrado

    |g|² tiene cola de potencia −3/2, integrable en la recta.

  2. Sin cuadrar

    La cola de g tiene potencia −3/4 y no es integrable.

  3. Resultado

    g∈L² pero g∉L¹: unitariedad sí se aplica, la cota anterior con ||g||₁ finita no.

Resultado

Cada estimación requiere su propio espacio de datos.

Lo comprobamos

No se concluye que la solución de este dato sea necesariamente no acotada en todo tiempo.

6. Strichartz controla juntas las variables de tiempo y espacio

La norma Lᑫ del tiempo con valores en Lʳ del espacio integra cuánto y durante cuánto se concentra la solución. Para Schrödinger libre en Rⁿ, los pares admisibles cumplen 2≤q,r≤∞ y 2/q+n/r=n/2, excluyendo (q,r,n)=(2,∞,2).

Entonces ||U(t)u₀||LᑫₜLʳₓ≤C||u₀||₂ en toda la recta temporal.

La relación de exponentes refleja la escala; el teorema usa dispersión y dualidad, y no se reduce a sumar cotas puntuales con un exponente arbitrario. El par (∞,2) es la conservación de energía. Las versiones forzadas de Duhamel emplean pares y exponentes duales con las restricciones correspondientes.

Comprobar un par en tres dimensiones

q=2 y r=6.

  1. Tiempo

    2/q=1.
  2. Espacio

    n/r=3/6=1/2.
  3. Sumar

    1+1/2=3/2=n/2; el triple no es el excluido.
Resultado

El par (2,6) es admisible en dimensión tres.

Lo comprobamos

La conclusión mezcla norma L² temporal con L⁶ espacial.

El extremo que la cuenta de escala no basta para aceptar

n=2, q=2 y r=∞.

  1. Escala

    2/q+n/r=1+0=1=n/2.
  2. Restricción

    Es precisamente el triple (2,∞,2), excluido del teorema.

  3. Alternativa

    En n=2, q=r=4 sí cumple la relación y no está excluido.
Resultado

Pasar la igualdad algebraica no elimina la excepción analítica.

Lo comprobamos

En dimensión uno el par (4,∞) sí es admisible para el modelo libre.

7. Otra fase produce otras velocidades y otras tasas

En Airy, uₜ+uₓₓₓ=0.
Para e elevado a i(kx−ωt), uₜ=−iωu y uₓₓₓ=−ik³u: la relación es ω=−k³ y la velocidad de grupo −3k².
Su escala dispersiva L¹→L∞ es |t|⁻¹ᐟ³.

KdV añade un término no lineal como 6uuₓ.

La ecuación de ondas tiene ramas ω=±|ξ|.

En una dimensión, cada rama puede trasladar un perfil sin el ensanchamiento de Schrödinger; en dimensiones mayores interviene la separación de direcciones. Hay que mirar fase, dimensión y norma en cada modelo.

Cuidar el signo de una tercera derivada

Airy con k=2.

  1. Derivadas

    (ik)³=−ik³, no +ik³.
  2. Relación

    ω=−2³=−8.
  3. Grupo

    ω′(2)=−3·2²=−12.
Resultado

Cambiar el signo de la ecuación cambiaría la fase.

Lo comprobamos

La tasa cúbica tampoco se reemplaza por la tasa cuadrática de Schrödinger.

Una onda que mantiene su forma

u(x,t)=g(x−t), g suave, para uₜₜ−uₓₓ=0.

  1. Tiempo

    uₜₜ=g″(x−t).
  2. Espacio

    uₓₓ=g″(x−t).
  3. Cancelar

    La ecuación se cumple y el perfil solo se desplaza.

Resultado

Una evolución ondulatoria no obliga a que cada perfil se ensanche.

Lo comprobamos

Este ejemplo es unidimensional y no afirma ausencia de decaimiento geométrico en dimensiones mayores.

8. La fórmula integral separa evolución libre e interacción no lineal

Para NLS fijamos i uₜ+Δu=μ|u|ᵖ⁻¹u con μ real y p>1.
Su fórmula integral es u(t)=U(t)u_0-i\mu\int_0^tU(t-s)(|u|^{p-1}u)(s)\,ds.

La integral registra lo que aporta la interacción en cada tiempo y después propaga ese aporte.

Buscar un punto fijo necesita espacios y estimaciones dependientes de n,p,s y de la regularidad de la potencia. Para soluciones suficientemente regulares se conservan masa y energía E=∫|∇u|²/2+μ∫|u|ᵖ⁺¹/(p+1). Si μ>0 ambos términos son no negativos; con μ<0 la energía puede dejar de controlar el gradiente.

Recuperar el caso libre

μ=0 en Duhamel.

  1. Integral

    El factor μ hace cero todo el término no lineal.

  2. Solución

    Queda u(t)=U(t)u₀.
  3. Comprobar

    Se recuperan grupo y conservación L² del modelo lineal.

Resultado

La fórmula conserva el signo −iμ que corresponde a la ecuación elegida.

Lo comprobamos

Con μ≠0 no basta aplicar el decaimiento libre directamente a toda la solución.

Una solución homogénea no lineal

En un toro, dato constante complejo A.

  1. Espacio

    Δu=0 para una función constante en x.
  2. Propuesta

    u(t)=A exp(−iμ|A|ᵖ⁻¹t).
  3. Sustituir

    i uₜ=μ|A|ᵖ⁻¹u y |u|=|A|; satisface la ecuación.
Resultado

Puede cambiar la fase sin cambiar el módulo.

Lo comprobamos

En toda Rⁿ una constante no nula no es L²; aquí se ha especificado el toro.

9. La escala y las cantidades conservadas no resuelven solas el problema global

La transformación NLS es u_\lambda(t,x)=\lambda^{2/(p-1)}u(\lambda^2t,\lambda x).

Su orden crítico es s꜀=n/2−2/(p−1).
La masa es crítica si s꜀=0; la energía, si s꜀=1.
Para n≥3, este último caso corresponde a p=1+4/(n−2).

Conservar una cantidad no garantiza por sí solo regularidad global en un régimen que no controla. En el caso focalizante pueden existir perfiles persistentes y explosión bajo hipótesis específicas. Scattering significa aproximarse, en una norma declarada, a U(t)u₊ al tender t a infinito; es una conclusión adicional, no un sinónimo de existir.

Comparar dos potencias en dimensión tres

NLS cúbica p=3 y quíntica p=5.

  1. Cúbica

    s꜀=3/2−2/2=1/2.
  2. Quíntica

    s꜀=3/2−2/4=1.
  3. Interpretar

    La quíntica es crítica para energía; las dos no tienen el mismo orden crítico.

Resultado

La dimensión y la potencia deben entrar juntas en la cuenta.

Lo comprobamos

Esta clasificación de escala no demuestra un resultado global para todos los datos.

Un perfil focalizante que no se ensancha

En una dimensión, μ=−1, p=3 y Q(x)=√2 sech(x).

  1. Identidad

    Q″−Q+Q³=0, como se comprueba derivando sech.
  2. Solución

    u=eⁱᵗQ da i uₜ+uₓₓ=eⁱᵗ(−Q+Q″)=−eⁱᵗQ³.
  3. Módulo

    |u(x,t)|=Q(x), localizado y fijo; no obedece el decaimiento puntual del propagador libre.

Resultado

Es una onda estacionaria no lineal de masa finita.

Lo comprobamos

No contradice la cota libre: la ecuación incluye una interacción focalizante.

Ahora lo resolvemos de principio a fin.

EJEMPLO RESUELTO 1

Una frecuencia

Para k=3 en Schrödinger libre, calcula ω y v_g.

  1. Ecuación

    Fijamos i uₜ+uₓₓ=0 y una onda de fase kx−ωt.
  2. Relación

    La sustitución da ω−k²=0, por tanto ω=k².
  3. Dato

    k=3.
  4. Frecuencia

    ω=3²=9.
  5. Grupo

    v de grupo=dω/dk=2k=6.
  6. Distinguir

    La velocidad de fase sería ω/k=3.
    El resultado pedido es ω=9 y v de grupo=6.
Resultado

ω=9 y v_g=6.

Lo comprobamos

La velocidad de fase sería 3.

EJEMPLO RESUELTO 2

Norma conservada

Si ||u₀||₂=5 en la evolución libre, calcula ||u(t)||₂.

  1. Dato

    ||u₀||₂=5 en la evolución libre sobre Rⁿ.
  2. Propagador

    û(t,ξ)=e elevado a −it|ξ|² por û₀(ξ).
  3. Módulo

    La exponencial tiene módulo uno para t y ξ reales.

  4. Cuadrado

    La integral de |û(t,ξ)|² es la misma que inicialmente.

  5. Plancherel

    Esto conserva la norma L² de u en el espacio.

  6. Resultado

    ||u(t)||₂=5 para todo t.

    No implica que la norma suprema ni la forma del perfil sean constantes.

Resultado

La norma sigue siendo 5.

Lo comprobamos

No implica que la norma supremo sea constante.

Los símbolos que hemos usado

u(t,x)
Función compleja en tiempo real t y posición x∈Rⁿ, salvo que se indique otro dominio.
e^i(x·ξ−ωt)
Onda plana que fija los signos de la relación de dispersión ω=ω(ξ).
v_f=ω(k)/k
Velocidad de fase en una dimensión para k≠0.
v_g=∇ξω
Velocidad de grupo de un paquete concentrado cerca de una frecuencia.
U(t)=e^itΔ
Propagador de i uₜ+Δu=0; su multiplicador es e^−it|ξ|².
||u||Lᑫ_tLʳ_x
Norma que toma Lʳ en espacio y después Lᑫ en tiempo; los exponentes se eligen según la ecuación.
μ, p
En NLS, i uₜ+Δu=μ|u|^(p−1)u; μ=1 es defocalizante y μ=−1 focalizante, con p>1.
s_c=n/2−2/(p−1)
Orden crítico de Sobolev homogéneo para la escala de esa NLS.

Ahora te toca a ti.

Tienes cinco ejercicios. Inténtalos; las pistas te ayudan si te atascas.

Practicar este tema →

Para ampliar