Capítulo 4

Funciones de operadores

1. ¿Por qué?
De nuevo recurrimos a la Mecánica Cuántica para justificar la necesidad de este
estudio en Física: escribir objetos tan importantes como la evolución de u sistema
físico necesitará del cálculo de la exponencial de un operador (la energía, en este
caso). En muchas otras áreas, como la mecánica clásica, serán también importantes
para el estudio de simetrías, o en el estudio de la estabilidad de sistemas dinámicos.

2. Operadores nilpotentes
Ejemplo 4.1. Vamos a considerar una aplicación lineal representada en una cierta
base en la forma:
 
0 1 2
A = 0 0 3 .
0 0 0
Se trata de una matriz que ya se encuentra representada en una forma triangular,
pero en la que las entradas de la diagonal son nulas.
Es sencillo verificar que si calculamos las potencias de la matriz:
   
0 0 3 0 0 0
A2 = 0 0 0 A3 = 0 0 0 .
0 0 0 0 0 0

Vemos pues que la composición de la matriz consigo misma un cierto número de veces
nos conduce a la matriz nula.

Es evidente también que esta no es una propiedad de todas las apliaciones linea-
les, pues basta considerar el endomorfismo identidad para darse cuenta de que hay
aplicaciones para las que no existe una potencia que las anule.
Definición 4.1. Sea f : V → V un endomorfismo. Diremos que f es un endo-
morfismo nilpotente si existe un número natural N ∈ N que verifica que
N
� �� �
f N = f ◦ . . . ◦ f = 0.
Análogamente diremos que una matriz es nilpotente cuando existe un número natural
N ∈ N que verifica que
N
N � �� �
A = A·...·A = 0
En ambos casos denominaremos índice de nilpotencia al exponente N.
111
112 4. FUNCIONES DE OPERADORES

El interés de los operadores nilpotentes se entenderá en las secciones siguientes.

Pero vamos a presentar un par de resultados previos que proporcionan información
interesante:

Lema 4.1. Sea f un operador nilpotente con índice de nilpotencia N. El único

autovalor de f es λ = 0.

Demostración. Sea |v � un autovector con autovalor λ. Entonces,

f |v � = λ|v � ⇒ f k |v � = λk |v � ∀k.

Entonces, sabemos que

�
f N |v � = λN |v �
f N |v � = 0 por ser nilpotente de orden N

Esto prueba el resultado. �

La conclusión del resultado anterior es por tanto que el polinomio característico

de un endomorfismo nilpotente tendrá la forma:

Pf (λ) = λn . (196)

Y dado que sabemos que el polinomio mínimo debe dividir a Pf , sabemos que se
debe cumplir que

mf (λ) = λk k ≤ n.

Pero sabemos que por ser f nilpotente con orden N, el orden más bajo para el que se
anula la ecuación

fk =0

es k = N. Concluimos por tanto:

Lema 4.2. Sea f : V → V un endomorfismo nilpotente de orden N definido

en un espacio vectorial V de dimensión n. Entonces su polinomio característico y su
polinomio mínimo corresponden a:

Pf (λ) = λn mf (λ) = λN . (197)

Consideremos entonces la base asociada B a la forma canónica de Jordan que

vimos en la Sección 9 del capítulo anterior. Respecto a ella, la expresión matricial del
 3. FUNCIONES DE OPERADORES 113

endomorfismo adopta una estructura diagonal por cajas, en la forma:

 
0 1 0 ... 0
 0 0 1 ... 0 
 
 .. .. .. .. . 
 .
 . . . .. 0 ... 0 

 0 0 0 ... 1 
 
 0 0 0 ... 0 
 
 0 1 0 ... 0 
 
 0 0 1 ... 0 
 
 .. .. .. .. .. 
 0 
 . . . . . ... 0 
f |B →   (198)
 0 0 0 ... 1 
 
 0 0 0 ... 0 
 
 .. 
 0 0 . 0 
 
 0 1 0 ... 0 
 
 0 0 1 ... 0 
 
 .. .. .. . . .. 
 0 0 ... . . . . . 
 
 0 0 0 ... 1 
0 0 0 ... 0
El tamaño relativo de cada caja queda determinado por la longitud de las cadenas de
vectores de la base de la forma canónica. Sabemos, sin embargo que todas ellas son
de longitud menor o igual que N (el índice de nilpotencia) por ser éste el índice de la
raíz λ = 0 en el polinomio mínimo
Pero con ello hemos determinado:
Lema 4.3. Sea f : V → V un endomorfismo nilpotente. Entonces, existe una
base de V en la que la matriz que representa a f es una matriz estrictamente trian-
gular superior, es decir, una matriz triangular superior con todas las entradas sobre la
diagonal nulas.
Ejercicio 4.1. Demostrar que la parte fuera de la diagonal de la forma canónica
de Jordan representa siempre un operador nilpotente.

3. Funciones de operadores
3.1. Definición general. En muchas ocasiones en Física va a ser preciso definir
magnitudes definidas a partir de otras, que se corresponden con operadores asociados a
magnitudes simples. El ejemplo más conocido lo da el operador evolución en Mecánica
Cuántica, que se define como la exponencial del operador que representa la energía
del sistema (el operador Hamiltoniano).
En esta sección vamos a ocuparnos de estas definiciones funciones de operadores.
Nos restringiremos no obstante a aquellas funciones para las cuales dispongamos de
una serie de Taylor convergente en un entorno alrededor de cada punto de su dominio,
que es lo que se denomina función analítica. Una función de este tipo se escribirá
entonces, para cualquier punto x ∈ D donde D es el abierto de K n donde esté definida:
∞
� f (k) (x0 )
f (x) = (x − x0 )k
k!
k=0
114 4. FUNCIONES DE OPERADORES

Si consideramos las series en torno al origen, tendremos:

∞
� f (k) (0)
f (x) = xk
k!
k=0

Es esta la expresión que vamos a considerar para la definición de la función de un

endomorfismo.
Definición 4.2. Sea f : K → K una función analítica definida sobre el cuerpo K, y
sea H : V → V un endomorfismo definido sobre un K–espacio vectorial V . Entonces,
definimos la función del endomorfismo H como el K–endomorfismo definido por la
serie:
∞ k
� f (k) (0) k � �� �
f (H) = H H k := H ◦ . . . ◦ H (199)
k!
k=0

La definición permite definir operaciones sobre el conjunto End(V ), mucho más

generales que los polinomios que habíamos introducido al considerar el polinomio mí-
nimo. Se trata de objetos que, por lo general, resultan muy complicados de calcular y
manejar. No obstante, vamos a ver que las definiciones que hicimos de formas canó-
nicas nos permiten construir, sobre las bases correspondientes, expresiones sencillas
de manejar.

3.2. Operadores diagonalizables. Consideremos en primer lugar que el endo-

morfismo H es diagonalizable, esto es, existe una descomposición en subespacios pro-
pios y una base de vectores propios, en la que la expresión de la matriz que representa
a H es
 
λ1 ... 0
d  .. .. .. 
AH =  . . . (200)
0 ... λn
Y dada una expresión matricial AH en una base arbitraria, sabemos que existe una
aplicación lineal invertible, que tiene una representación matricial dada por una matriz
S que verifica que
AH = SAdH S −1
En estas circunstancias, en inmediato probar que:
Lema 4.4. Sea H : V → V un endomorfismo diagonalizable. Entonces, la repre-
sentación matricial de cualquier potencia del endomorfismo está dada por:
 k 
λ1 . . . 0
. .
H k → S  .. . . . ..  S −1 , (201)
0 ... λkn
donde S es la matriz del cambio de base anterior.
Demostración. En una base arbitraria, la expresión de las potencias toma la forma:
k
H k = SAdH S −1 ����
. . . SAdH S −1 = S(AdH )k S −1 .
�
 3. FUNCIONES DE OPERADORES 115

Esto nos permite extender de forma sencilla el cálculo de cualquier función del
operador H en función unicamente de los valores sobre la diagonal:
Teorema 4.1. Sea f : Kn → K una función escalar sobre el cuerpo K y H : V → V
un endomorfismo diagonalizable. Entonces, la expresión matricial de la función f (H)
corresponde a  
f (λ1 ) . . . 0
 . .. ..  −1
f (H) → S  .. . . S (202)
0 ... f (λn )

Demostración. Basta con emplear el Lema 4.4 y escribir:

∞
�∞ �
� f (k) (0) � f (k) (0)
f (H) → (SAdH S −1 )k = S (AdH )k S −1 = Sf (AdH )S −1
k! k!
k=0 k=0

Vamos a analizar ahora algunas de las funciones de operadores que con más
frecuencia encontraremos en Física:
3.2.1. La función exponencial. En este caso,
∞
� 1 k
f (x) = exp x = x .
k!
k=0

Por tanto, tenemos que, dada un endomorfismo H, diagonalizable, con expresión

diagonal
 
λ1 . . . 0
. .
AdH =  .. . . . ..  ;
0 ... λn
y una expresión matricial en una base arbitraria conjugada de la anterior por una matriz
invertible S como
AH = SAdH S −1 ,
Definición 4.3. Sea H : V → V un endomorfismo definido sobre un espacio
vectorial complejo. Entonces, la expresión de su exponencial corresponderá a:
 λ1 
e ... 0
 . .. ..  −1
Aexp H = S  .. . . S .
0 ... e λn

Ejercicio 4.2. Sea H : V → V un endomorfismo diagonalizable. Probar que

det exp H = exp(TrH)
Ejercicio 4.3. Sea H : V → V un endomorfismo diagonalizable y P : V → V otro
endomorfismo invertible. Entonces
exp(Y ◦ H ◦ Y −1 ) = Y ◦ exp H ◦ Y −1
Ejercicio 4.4. Sea H : V → V un endomorfismo diagonalizable. Probar que
(exp H)−1 = exp(−H)
116 4. FUNCIONES DE OPERADORES

3.2.2. La función logaritmo. La función logaritmo real sabemos que está defini-
da unicamente para el conjunto de puntos no negativos de la recta real. En el caso del
campo complejo, la función logaritmo tiene particularidades que la hacen de especial
dificultad respecto a su análogo real, pues se trata de una función multivaluada:

Figura 1. Expresión de la función logaritmo compleja

Es este un dato a tener en cuenta, si bien no afecta a la definición del logaritmo

de un endomorfismo sobre un espacio vectorial complejo, salvo en lo que se refiere
a la interpretación de la expresión matricial resultante. Y tanto en el caso real como
en el caso complejo es necesario recordar la existencia de una singularidad sobre el
origen. Con estos condicionantes, consideraremos:
Definición 4.4. Sea H : V → V un endomorfismo definido sobre un espacio
vectorial complejo. Entonces, si H es diagonalizable y no tiene autovalores nulos, la
expresión del matricial operador logaritmo log H se escribirá:
 
log λ1 . . . 0
 . .. ..  −1
Alog H = S  .. . . S
0 ... log λn
3.2.3. La función raiz cuadrada. La raíz cuadrada es una función bien definida
para valores reales no negativos cuando se considera como función real y bien definida
en cualquier punto del plano complejo cuando se considera como función compleja.
Es una función analítica en ambos casos.
Se trata pues de determinar cuándo podemos obtener un operador tal que elevado
al cuadrado sea igual a uno dado:
bien sobre el cuerpo de los reales
Definición 4.5. Sea H : V → V un endomorfismo definido sobre un espacio
vectorial real. Si H es diagonalizable y todos sus valores propios son no negativos,
entonces existe su raíz cuadrada que corresponde al endomorfismo que tiene como
expresión matricial
 1/2 
λ1 ... 0
 . .. ..  −1
H 1/2 → AH1/2 = S  .. . . S
1/2
0 . . . λn
 3. FUNCIONES DE OPERADORES 117

bien sobre el cuerpo de los complejos:

Definición 4.6. Sea H : V → V un endomorfismo definido sobre un espacio

vectorial complejo. Entonces, si H es diagonalizable la expresión del matricial operador
raiz cuadrada H 1/2 se escribirá:
 1/2 
λ1 ... 0
 . .. ..  −1
H 1/2 → AH1/2 = S  .. . . S
1/2
0 . . . λn
.

3.3. Funciones de operadores no diagonalizables. La definición formal en el

caso de endomorfismos no diagonalizables es completamente análoga, aunque desde
un punto de vista operativo los resultados san muy diferentes. Por ejemplo, podemos
comenzar analizando qué ocurre con operadores nilpotentes: es evidente que, al anu-
larse a partir de una determinada potencia, las series infinitas que definen las funciones
se convierten en sumas finitas, más o menos sencillas de manejar.

Ejercicio 4.5. Sea N : V → V un endomorfismo nilpotente de orden q. Entonces,

N2 N3 N q−1
exp N = Id + N + + + ... + .
2 6 (q − 1)!

Evidentemente, este ejemplo no proporciona más información que la referida a

operadores nilpotentes. Sin embargo, recordemos que la forma canónica de Jordan
estaba disponible para todos los endomorfismos, y, precisamente, se compone de la
suma de un operador diagonal y un operador nilpotente. Es por ello natural pensar
en usar la forma de Jordan, tal como en la sección anterior hemos usado la forma
diagonal.

3.4. La función exponencial. Consideremos entonces un endomorfismo cual-

quiera f : V → V en una base arbitraria. Sabemos que siempre existe una base en la
que la representación matricial de f corresponde a la forma canónica de Jordan, es
decir, sabemos que:

Af = SJS −1
118 4. FUNCIONES DE OPERADORES

donde
 
λ1 1 0 ... 0
 0 λ1 1 ... 0 
 

 .. .. .. .. .. 

. . . . . 0 ... 0
 
 0 0 0 ... 1 
 
 
 0 0 0 ... λ1 
 
 λ2 1 0 ... 0 
 
 0 λ2 1 ... 0 
 
 .. .. .. .. .. 
 0 . . . . . ... 0 
J=



 0 0 0 ... 1 
 0 0 0 ... λ2 
 
 .. 
 . 
 0 0 0 
 
 λk 1 0 ... 0 
 
 0 λk 1 ... 0 
 
 .. .. .. .. .. 
 0 0 ... . . . . . 
 
 0 0 0 ... 1 
0 0 0 ... λk

La expresión de la función exponencial de Af se escribe por tanto:

∞
� ∞
1 k � 1 1
exp Af = Af = (SJS −1 )k = S( J k )S −1
k! k! k!
k=0 k=0

Al tener la matriz estructura por bloques, es sencillo entender que cada bloque es
independiente de los demás, de forma que se puede escribir:
 
exp J1 0 ... 0
 0 exp J2 ... 0 
 
exp J =  .. 
 0 0 . 0 
0 0 ... exp Jk
Podemos escribir las matrices diagonal y nilpotente que cada bloque definen la forma
de Jordan:
Jj = Λ + N
con
 
λj ... 0
 .. .. .. 
Λ=. . .
0 ... λj
y
 
0 1 ... 0
 .. .. . . .
. . . .. 
N=

;

0 0 . . . 1
0 0 ... 0
 3. FUNCIONES DE OPERADORES 119

con las potencias de la matriz nilpotente correspondiendo a:

j +1
↓
 
0 ... 0 1 ... 0
 .. .. .. .. . . .. 
 .
 . . . . . 

Nj = 
 0 . . . 0 0 . . . 1  ← nj − j.

 . .. .. .. .. .. 
 .. . . . . . 
0 ... 0 0 ... 0

Con esta expresión podemos escribir las expresiones de las potencias de la forma
de Jordan:
� �k �
k k
Jj = (Λ + N) = Λp N k−p .
p
p

Un sencillo cálculo nos conduce a:

p+1
↓
 
0 ... 0 1 ... 0
 .. .. .. .. . . . 
 . . . . . .. 
�k 
k−p 

Jjk = k! 
p=0 p!(k−p)! λj  0 . . . 0 0 . . . 1  ← np − p.
 . . . . . . 
 .. .. .. .. .. .. 
0 ... 0 0 ... 0

Pero entonces, en cada entrada de la matriz anterior, vamos a encontrar el desa-

rrollo de la exponencial del autovalor, dividido por un número factorial (el coeficiente
p!):
 k 
λj 1 k−1 1 k−2 1 k−nj +1
λ λ . . . λ
 k! 1!(k−1)!k j 2!(k−2)! j (nj −1)!(k−nj +1)! j 
 λj 1 k−1 1 k−nj +2 
 0 k! 1!(k−1)! λ j . . . (nj −2)!(k−nj +2)! λ j 
1 k  .. .. .. .. ..

Jj = 



k!  . . . . . 
 0 ... 0 ... 1 k−1 
 1!(k−1)! λj 
λkj
0 ... 0 ... k!

Por consiguiente, si sumamos todos los elementos de la serie, vamos a recuperar la

exponencial del autovalor, pero siempre con el factorial extra:
 λj λj λj

e λj e1! e2! . . . (nje−1)!
 
 0 e λj e λ j . . . e λj 
 1! (nj −2)! 
 . .. .. .. .. 
exp Jj =  .  (203)
 . . . . . 
 e λj 
 0 ... 0 ... 1!

0 ... 0 ... e λj
120 4. FUNCIONES DE OPERADORES

3.5. La función logaritmo. EL cálculo de la función logaritmo es similar al

anterior, pero debemos emplear un desarrollo en serie diferente. Tal como ocurría en
el caso de la exponencial, el proceso se base en dos etapas:
consideramos la forma canónica de Jordan, es decir, escribimos la matriz en
una base arbitraria como un producto
Af = SJS −1 ,
donde
 
J1 0 ... 0
 0 J2 ... 0 
 
J= .. 
 0 0 . 0 
0 0 ... Jk
y
 
λi 1 0 ... 0
0 λi 1 ... 0
 
. .. .. . . .
Ji =  .. . . . .. 
 
0 0 0 ... 1
0 0 0 . . . λi
Esta descomposición permite escribir
log Af = S log JS −1
Por conveniencia, vamos a escribir la matriz J de una forma alternativa,
separándola como un producto, más adecuado para el cálculo del logaritmo:
J = Λ(Id + K)
donde
   
Λ1 0 ... 0 K1 0 ... 0
 0 Λ2 ... 0   0 K2 ... 0 
   
Λ= .. ; K= .. 
 0 0 . 0   0 0 . 0 
0 0 . . . Λk 0 0 ... Kk
y
   
λi 0 0 ... 0 0 λ−1
i 0 ... 0
0 λi 0 ... 0 0 0 λ−1 ... 0 
   i 
. .. .. . . ..  . .. .. .. .. 
Λi =  .. . . . . Ki =  .. . . . . 
   
0 0 0 ... 0 0 0 0 ... λ −1 
i
0 0 0 . . . λi 0 0 0 ... 0
Escribiremos entonces, haciendo uso de las propiedades del logaritmo que:
log J = log Λ + log(I + K).
El primer sumando es el logaritmo de una matriz diagonal, que ya hemos
aprendido a hacer en la sección anterior. Para el caso del segundo sumando,
procedemos como sigue:
 3. FUNCIONES DE OPERADORES 121

Cada uno de los bloques no se ve afectado por los productos de los demás, de
forma que podemos pensar en el logaritmo de cada uno de forma individual.
Pensaremos entonces en desarrollar los términos de la forma
log(I + Kj ).
En este punto emplearemos el desarrollo de Taylor de la función
x2 x3 x4 xk
log(1 + x) = x − + − + . . . + (−1)k+1 + ...
2 3 4 k
Afirmamos por tanto que
Kj2 Kj3 Kj4 Kjk
log(I + Kj ) = Kj − + − + . . . + (−1)k+1 + ...
2 3 4 k
La expresión anterior puede formularse para cualquier matriz, pero en general,
nunca converge. En el caso que nos ocupa, sin embargo, las matrices K j son nilpo-
tentes en todos los casos. Esto implica que en vez de series infinitas tendremos una
suma finita (y por tanto necesariamente convergente):
q
Kj2 Kj3 Kj4 Kj j
log(I + Kj ) = Kj − + − + . . . + (−1)k+1 ;
2 3 4 qj
donde qj es el orden de nilpotencia de la matriz Kj .
Concluimos por consiguiente que para cada uno de los bloques de la matriz J
encontramos la expresión
Lema 4.5. Para cada una de los bloques de la forma de Jordan tenemos:
q
Kj2 Kj3 Kj4 Kj j
log J = log Λj + Kj − + − + . . . + (−1)k+1 (204)
2 3 4 qj
Ejercicio 4.6. Determinar la forma de las matrices Kji , para cualquier valor de i .
La combinación de estas matrices en una única matriz que represente la forma
de Jordan completa es difícil de escribir de forma compacta, dado que cada uno de
los bloques puede tener distinto orden de nilpotencia. Es recomendable, no obstante:
Ejercicio 4.7. Determinar la forma del logaritmo de la matriz completa en el
caso de dos bloques, asumiendo que las dimensiones de los bloques son uno y dos,
respectivamente.
Ejercicio 4.8. Verificar que la expresión del logaritmo anterior coincide con la
expresión obtenida en el caso de una matriz diagonalizable.
122 4. FUNCIONES DE OPERADORES

Problemas
1. Ejercicio 4.1.
2. Ejercicio 4.2.
3. Ejercicio 4.3.
4. Ejercicio 4.4.
5. Consideremos las matrices de Pauli, definidas sobre C2 :
� � � � � �
0 1 0 −i 1 0
σ1 = σ2 = σ3 = . (205)
1 0 i 0 0 −1
Determinar la expresión de
a) exp(i� πσ1 ) �
b) exp −i π4 σ2
c) exp(i ασ3 ), con α ∈ (−π, π].
6. Consideremos V = R2 y la matriz
� �
0 −1
R= .
1 0
Determinar exp R y exp(i R).
7. Consideremos las matrices sobre R3 :
     
0 0 0 0 0 1 0 −1 0
L1 = 0 0 −1 L2 =  0 0 0 L3 = 1 0 0 . (206)
0 1 0 −1 0 0 0 0 0
Calcular
2k+1
a) L2k
1 y L1 para cualquier k ∈ N
b) Idem para L2 y L3
c) exp(i αL3 ), con α ∈ (−π, π].
8. Sea H la matriz definida como
� σ ≡ B1 σ1 + B2 σ2 + B3 σ3 .
H = B�
Esta es la matriz que representa la energía de un espín acoplado a un campo
magnético B,� definida sobre un espacio de estados que es igual a C 2 .
a) Determinar los autovalores y autovectores de H en función del valor del
campo magnético B.
b) Determinar la evolución de un estado
� �
ψ1
|ψ� =
ψ2
sabiendo que está dada por la curva |ψ(t)� definida en C2 por la expresión
|ψ(t)� = exp(iHt)|ψ(0)�.
9. Ejercicio 4.5.
10. Ejercicio 4.6
11. Ejercicio 4.7.
12. Ejercicio 4.8.
13. Sea V = C2 . Consideremos un elemento arbitrario de V
� �
ψ1
ψ= |ψ1 |2 + |ψ2 |2 = 1
ψ2