Efecto Zeeman

Las líneas espectrales de la lámpara de vapor de mercurio en la longitud de onda 546.1 nm, mostrando un efecto Zeeman anómalo. (A) Sin campo magnético. (B) Con el campo magnético, las líneas espectrales se dividen como un efecto Zeeman transversal. (C) Con campo magnético, dividido como efecto Zeeman longitudinal. Las líneas espectrales se obtuvieron mediante un interferómetro de Fabry-Pérot.

[[Archivo:Breit-rabi-Zeeman-en.svg|miniaturadeimagen|420x420px| División de Zeeman del nivel 5 de [[Rubidio|Plantilla:ExpRb]], incluida la división de estructura fina y estructura hiperfina. Aquí F = J + I, donde I es el espín nuclear (para Plantilla:ExpRb, I = Plantilla:Frac)]] Archivo:Explanation of how the magnetic field on a star affects the light emitted.webm El efecto Zeeman, que lleva el nombre del físico neerlandés Pieter Zeeman, es el efecto de la división de una línea espectral en varios componentes en presencia de un campo magnético estático. Es análogo al efecto Stark, la división de una línea espectral en varios componentes en presencia de un campo eléctrico. También similar al efecto Stark, las transiciones entre diferentes componentes tienen, en general, diferentes intensidades, estando algunas totalmente prohibidas (en la aproximación dipolar), según se rige por las reglas de selección.

Dado que la distancia entre los subniveles de Zeeman es una función de la intensidad del campo magnético, este efecto se puede utilizar para medir la intensidad del campo magnético, por ejemplo, el del Sol y otras estrellas o en plasmas de laboratorio. El efecto Zeeman es muy importante en aplicaciones como la espectroscopia de resonancia magnética nuclear, la espectroscopia de resonancia de espín electrónico, la formación de imágenes por resonancia magnética (MRI) y la espectroscopia de Mössbauer. También se puede utilizar para mejorar la precisión en la espectroscopia de absorción atómica. Una teoría sobre el sentido magnético de las aves asume que una proteína en la retina cambia debido al efecto Zeeman.

Cuando las líneas espectrales son líneas de absorción, el efecto se denomina efecto Zeeman inverso.

Nomenclatura

Históricamente, se distingue entre el efecto Zeeman normal y anómalo (descubierto por Thomas Preston en Dublín, Irlanda). El efecto anómalo aparece en las transiciones en las que el spin neto de los electrones no es cero. Se le llamó "anómalo" porque el espín del electrón aún no se había descubierto, por lo que no había una buena explicación en el momento en que Zeeman observó el efecto.

A mayor intensidad de campo magnético, el efecto deja de ser lineal. A una fuerza de campo aún mayor, cuando la fuerza del campo externo es comparable a la fuerza del campo interno del átomo, el acoplamiento de electrones se altera y las líneas espectrales se reorganizan. Esto se llama efecto Paschen-Back.

En la literatura científica moderna, estos términos se utilizan raramente, con una tendencia a utilizar solo "efecto Zeeman".

Presentación teórica

El hamiltoniano total de un átomo en un campo magnético es

$H=H_0+V_{\mathrm{M}},$

donde $H_0$ es el hamiltoniano imperturbable del átomo, y $V_{\mathrm{M}}$ es la perturbación debida al campo magnético:

$V_{\mathrm{M}}=-\overrightarrow{\mu }\cdot \overrightarrow{B},$

donde $\overrightarrow{\mu }$ es el momento magnético del átomo. El momento magnético consta de las partes electrónicas y nucleares; sin embargo, este último es muchos órdenes de magnitud más pequeño y se pasará por alto aquí. Por lo tanto,

$\overrightarrow{\mu }\approx -\frac{{\mu }_{\mathrm{B}}g\overrightarrow{J}}{\hslash },$

donde ${\mu }_{\mathrm{B}}$ es el magnetón de Bohr, $\overrightarrow{J}$ es el momento angular electrónico total, y $g$ es el factor g de Landé. Un enfoque más preciso es tener en cuenta que el operador del momento magnético de un electrón es una suma de las contribuciones del momento angular orbital $\overrightarrow{L}$ y el momento angular de giro $\overrightarrow{S}$, con cada uno multiplicado por la proporción giromagnética apropiada:

$\overrightarrow{\mu }=-\frac{{\mu }_{\mathrm{B}}(g_l\overrightarrow{L}+g_s\overrightarrow{S})}{\hslash },$

donde $g_l=1$ y $g_s\approx 2.0023192$ (este último se llama relación giromagnética anómala; la desviación del valor de 2 se debe a los efectos de la electrodinámica cuántica). En el caso del acoplamiento LS, se pueden sumar todos los electrones del átomo:

$g\overrightarrow{J}=⟨\sum _i(g_l\overrightarrow{l_i}+g_s\overrightarrow{s_i})⟩=⟨(g_l\overrightarrow{L}+g_s\overrightarrow{S})⟩,$

donde $\overrightarrow{L}$ y $\overrightarrow{S}$ son el momento orbital total y el spin del átomo, y el promedio se realiza sobre un estado con un valor dado del momento angular total.

Si el término de interacción $V_M$ es pequeño (menos que la estructura fina), puede tratarse como una perturbación; este es el efecto Zeeman propiamente dicho. En el efecto Paschen-Back, que se describe a continuación, $V_M$ excede el acoplamiento LS significativamente (pero aún es pequeño en comparación con $H_0$). En campos magnéticos ultrafuertes, la interacción del campo magnético puede exceder $H_0$, en cuyo caso el átomo ya no puede existir en su significado normal, y en su lugar se habla de niveles de Landau. Hay casos intermedios que son más complejos que estos casos límite.

Campo débil (efecto Zeeman)

Si la interacción espín-órbita domina sobre el efecto del campo magnético externo, $\overrightarrow{L}$ y $\overrightarrow{S}$ no se conservan por separado, solo el momento angular total $\overrightarrow{J}=\overrightarrow{L}+\overrightarrow{S}$ se conserva. Se puede pensar que los vectores de momento angular orbital y de espín precesan alrededor del vector de momento angular total (fijo) $\overrightarrow{J}$. El vector de giro (tiempo-) "promediado" es entonces la proyección del giro en la dirección de $\overrightarrow{J}$:

${\overrightarrow{S}}_{\mathrm{a}\mathrm{v}\mathrm{g}}=\frac{(\overrightarrow{S}\cdot \overrightarrow{J})}{J^2}\overrightarrow{J}$

y para el vector orbital (tiempo -) "promediado":

${\overrightarrow{L}}_{\mathrm{a}\mathrm{v}\mathrm{g}}=\frac{(\overrightarrow{L}\cdot \overrightarrow{J})}{J^2}\overrightarrow{J}.$

Por lo tanto,

$⟨V_{\mathrm{M}}⟩=\frac{{\mu }_{\mathrm{B}}}{\hslash }\overrightarrow{J}(g_L\frac{\overrightarrow{L}\cdot \overrightarrow{J}}{J^2}+g_S\frac{\overrightarrow{S}\cdot \overrightarrow{J}}{J^2})\cdot \overrightarrow{B}.$

Usando $\overrightarrow{L}=\overrightarrow{J}-\overrightarrow{S}$ y cuadrando ambos lados, obtenemos

$\overrightarrow{S}\cdot \overrightarrow{J}=\frac{1}{2}(J^2+S^2-L^2)=\frac{{\hslash }^2}{2}[j(j+1)-l(l+1)+s(s+1)],$

y: usando $\overrightarrow{S}=\overrightarrow{J}-\overrightarrow{L}$ y cuadrando ambos lados, obtenemos

$\overrightarrow{L}\cdot \overrightarrow{J}=\frac{1}{2}(J^2-S^2+L^2)=\frac{{\hslash }^2}{2}[j(j+1)+l(l+1)-s(s+1)].$

Combinando todo y tomando $J_z=\hslash m_j$, obtenemos la energía potencial magnética del átomo en el campo magnético externo aplicado,

$\begin{array}{cc}{V}_{\mathrm{M}}& ={\mu }_{\mathrm{B}}B{m}_j[g_L\frac{j(j+1)+l(l+1)-s(s+1)}{2j(j+1)}+g_S\frac{j(j+1)-l(l+1)+s(s+1)}{2j(j+1)}]\\ & ={\mu }_{\mathrm{B}}B{m}_j[1+(g_S-1)\frac{j(j+1)-l(l+1)+s(s+1)}{2j(j+1)}],\\ & ={\mu }_{\mathrm{B}}B{m}_j{g}_j\end{array}$

donde la cantidad entre corchetes es el factor g de Landé g_J del átomo ($g_L=1$ y $g_S\approx 2$) y $m_j$ es la componente z del momento angular total. Por un solo electrón por encima de las capas llenas $s=1/2$ y $j=l\pm s$, el factor g de Landé se puede simplificar en:

$g_j=1\pm \frac{g_S-1}{2l+1}$

Tomando $V_m$ para ser la perturbación, la corrección de Zeeman a la energía es

$\begin{array}{c}{E}_{\mathrm{Z}}^{(1)}=⟨nlj{m}_j|H_{\mathrm{Z}}^{\text{'}}|nlj{m}_j⟩=⟨V_M{⟩}_{\mathrm{\Psi }}={\mu }_{\mathrm{B}}{g}_J{B}_{\mathrm{e}\mathrm{x}\mathrm{t}}{m}_j\end{array}$

Ejemplo: transición Lyman-alfa en hidrógeno

La transición Lyman-alfa en hidrógeno en presencia de la interacción espín-órbita implica las transiciones

$2P_{1/2}\to 1S_{1/2}$ y $2P_{3/2}\to 1S_{1/2}.$

En presencia de un campo magnético externo, el efecto Zeeman de campo débil divide los niveles 1S_1/2 y 2P_1/2 en 2 estados cada uno ($m_j=1/2,-1/2$) y el nivel 2P_3/2 en 4 estados ($m_j=3/2,1/2,-1/2,-3/2$). Los factores g de Landé para los tres niveles son:

$g_J=2$ por $1S_{1/2}$ (j = 1/2, l = 0)

$g_J=2/3$ por $2P_{1/2}$ (j = 1/2, l = 1)

$g_J=4/3$ por $2P_{3/2}$ (j = 3/2, l = 1).

Tenga en cuenta en particular que el tamaño de la división de energía es diferente para los diferentes orbitales, porque los valores de g_J son diferentes. A la izquierda, se muestra una fina estructura dividida. Esta división ocurre incluso en ausencia de un campo magnético, ya que se debe al acoplamiento espín-órbita. Representado a la derecha es la división adicional de Zeeman, que ocurre en presencia de campos magnéticos.

Posibles transiciones para el efecto Zeeman débil

Estado inicial ( $n=2,l=1$ ) $\mid j,m_j⟩$	Estado final ( $n=1,l=0$ ) $\mid j,m_j⟩$	Perturbación energética
$|\frac{1}{2},\pm \frac{1}{2}⟩$	$|\frac{1}{2},\pm \frac{1}{2}⟩$	$\mp \frac{2}{3}{\mu }_{\mathrm{B}}B$
$|\frac{1}{2},\pm \frac{1}{2}⟩$	$|\frac{1}{2},\mp \frac{1}{2}⟩$	$\pm \frac{4}{3}{\mu }_{\mathrm{B}}B$
$|\frac{3}{2},\pm \frac{3}{2}⟩$	$|\frac{1}{2},\pm \frac{1}{2}⟩$	$\pm {\mu }_{\mathrm{B}}B$
$|\frac{3}{2},\pm \frac{1}{2}⟩$	$|\frac{1}{2},\pm \frac{1}{2}⟩$	$\mp \frac{1}{3}{\mu }_{\mathrm{B}}B$
$|\frac{3}{2},\pm \frac{1}{2}⟩$	$|\frac{1}{2},\mp \frac{1}{2}⟩$	$\pm \frac{5}{3}{\mu }_{\mathrm{B}}B$

Campo fuerte (efecto Paschen-Back)

El efecto Paschen-Back es la división de los niveles de energía atómica en presencia de un fuerte campo magnético. Esto ocurre cuando un campo magnético externo es lo suficientemente fuerte como para interrumpir el acoplamiento entre orbitales ($\overrightarrow{L}$) y girar ($\overrightarrow{S}$) momentos angulares. Este efecto es el límite de campo fuerte del efecto Zeeman. Cuando $s=0$, los dos efectos son equivalentes. El efecto lleva el nombre de los físicos alemanes Friedrich Paschen y Ernst E. A. Back.

Cuando la perturbación del campo magnético excede significativamente la interacción espín-órbita, se puede asumir con seguridad $[H_0,S]=0$. Esto permite que los valores esperados de $L_z$ y $S_z$ para ser fácilmente evaluado por un estado $|\psi ⟩$. Las energías son simplemente

$E_z=⟨\psi |H_0+\frac{B_z{\mu }_{\mathrm{B}}}{\hslash }(L_z+g_s{S}_z)|\psi ⟩=E_0+B_z{\mu }_{\mathrm{B}}(m_l+g_s{m}_s).$

Lo anterior puede interpretarse como implicando que el acoplamiento LS está completamente roto por el campo externo. Sin embargo $m_l$ y $m_s$ siguen siendo números cuánticos "buenos". Junto con las reglas de selección para una transición de dipolo eléctrico, es decir, $\mathrm{\Delta }s=0,\mathrm{\Delta }{m}_s=0,\mathrm{\Delta }l=\pm 1,\mathrm{\Delta }{m}_l=0,\pm 1$ esto permite ignorar por completo el grado de libertad de giro. Como resultado, solo tres líneas espectrales serán visibles, correspondientes a la $\mathrm{\Delta }{m}_l=0,\pm 1$ regla de selección. La división $\mathrm{\Delta }E=B{\mu }_{\mathrm{B}}\mathrm{\Delta }{m}_l$ es independiente de las energías no perturbadas y las configuraciones electrónicas de los niveles considerados. En general (si $s\ne 0$), estos tres componentes son en realidad grupos de varias transiciones cada uno, debido al acoplamiento de espín-órbita residual.

En general, ahora se debe agregar el acoplamiento de espín-órbita y las correcciones relativistas (que son del mismo orden, conocidas como "estructura fina") como una perturbación a estos niveles "imperturbables". La teoría de perturbación de primer orden con estas correcciones de estructura fina produce la siguiente fórmula para el átomo de hidrógeno en el límite de Paschen – Back:

$E_{z+fs}=E_z+\frac{m_e{c}^2{\alpha }^4}{2n^3}\{\frac{3}{4n}-\left[\frac{l(l+1)-m_l{m}_s}{l(l+1/2)(l+1)}\right]\}$

Posibles transiciones Lyman-alfa para el régimen fuerte

Estado inicial ( $n=2,l=1$ ) $\mid m_l,m_s⟩$	Perturbación energética inicial	Estado final ( $n=1,l=0$ ) $\mid m_l,m_s⟩$
$|1,\frac{1}{2}⟩$	$\pm 2{\mu }_{\mathrm{B}}{B}_z$	$|0,\frac{1}{2}⟩$
$|0,\frac{1}{2}⟩$	$+{\mu }_{\mathrm{B}}{B}_z$	$|0,\frac{1}{2}⟩$
$|1,-\frac{1}{2}⟩$	$0$	$|0,-\frac{1}{2}⟩$
$|-1,\frac{1}{2}⟩$	$0$	$|0,\frac{1}{2}⟩$
$|0,-\frac{1}{2}⟩$	$-{\mu }_{\mathrm{B}}{B}_z$	$|0,-\frac{1}{2}⟩$
$|-1,-\frac{1}{2}⟩$	$-2{\mu }_{\mathrm{B}}{B}_z$	$|0,-\frac{1}{2}⟩$

Campo intermedio para j = 1/2

En la aproximación del dipolo magnético, el hamiltoniano que incluye las interacciones hiperfina y de Zeeman es

$H=hA\overrightarrow{I}\cdot \overrightarrow{J}-\overrightarrow{\mu }\cdot \overrightarrow{B}$

$H=hA\overrightarrow{I}\cdot \overrightarrow{J}+({\mu }_{\mathrm{B}}{g}_J\overrightarrow{J}+{\mu }_{\mathrm{N}}{g}_I\overrightarrow{I})\cdot \overrightarrow{\mathrm{B}}$

donde $A$ es la división hiperfina (en Hz) en el campo magnético aplicado cero, ${\mu }_{\mathrm{B}}$ y ${\mu }_{\mathrm{N}}$ son el magnetón de Bohr y el magnetón nuclear respectivamente, $\overrightarrow{J}$ y $\overrightarrow{I}$ son los operadores de momento angular de electrones y nucleares y $g_J$ es el factor g de Landé:

$g_J=g_L\frac{J(J+1)+L(L+1)-S(S+1)}{2J(J+1)}+g_S\frac{J(J+1)-L(L+1)+S(S+1)}{2J(J+1)}$ .

En el caso de campos magnéticos débiles, la interacción Zeeman puede tratarse como una perturbación del $|F,m_f⟩$ base. En el régimen de campo alto, el campo magnético se vuelve tan fuerte que dominará el efecto Zeeman, y uno debe usar una base más completa de $|I,J,m_I,m_J⟩$ o solo $|m_I,m_J⟩$, ya que $I$ y $J$ será constante dentro de un nivel dado.

Para obtener una imagen completa, incluidas las intensidades de campo intermedias, debemos considerar los estados propios que son superposiciones de $|F,m_F⟩$ y $|m_I,m_J⟩$ estados base. Para $J=1/2$, el hamiltoniano se puede resolver analíticamente, dando como resultado la fórmula de Breit-Rabi. En particular, la interacción cuadrupolo eléctrico es cero para $L=0$ ($J=1/2$), por lo que esta fórmula es bastante precisa.

Ahora utilizamos operadores de escalera mecánica cuántica, que se definen para un operador de momento angular general $L$ como

$L_{\pm }\equiv L_x\pm i{L}_y$

Estos operadores de escalera tienen la propiedad

$L_{\pm }|L_{,}{m}_L⟩=\sqrt{(L\mp m_L)(L\pm m_L+1)}|L,m_L\pm 1⟩$

siempre y cuando $m_L$ se encuentra en el rango $-L,\dots ...,L$ (de lo contrario, devuelven cero). Usando operadores de escalera $J_{\pm }$ y $I_{\pm }$ podemos reescribir el hamiltoniano como

$H=hA{I}_z{J}_z+\frac{hA}{2}(J_{+}{I}_{-}+J_{-}{I}_{+})+{\mu }_{\mathrm{B}}B{g}_J{J}_z+{\mu }_{\mathrm{N}}B{g}_I{I}_z$

Ahora podemos ver que en todo momento, la proyección del momento angular total $m_F=m_J+m_I$ se conservará. Esto es porque ambos $J_z$ y $I_z$ dejar estados con definidas $m_J$ y $m_I$ sin cambios, mientras $J_{+}{I}_{-}$ y $J_{-}{I}_{+}$ o bien aumentar $m_J$ y disminuir $m_I$ o viceversa, por lo que la suma siempre no se ve afectada. Además, dado que $J=1/2$ solo hay dos valores posibles de $m_J$ las cuales son $\pm 1/2$. Por lo tanto, para cada valor de $m_F$ solo hay dos estados posibles, y podemos definirlos como base:

$|\pm ⟩\equiv |m_J=\pm 1/2,m_I=m_F\mp 1/2⟩$

Este par de estados es un sistema mecánico cuántico de dos niveles. Ahora podemos determinar los elementos de la matriz del hamiltoniano:

$⟨\pm |H|\pm ⟩=-\frac{1}{4}hA+{\mu }_{\mathrm{N}}B{g}_I{m}_F\pm \frac{1}{2}(hA{m}_F+{\mu }_{\mathrm{B}}B{g}_J-{\mu }_{\mathrm{N}}B{g}_I))$

$⟨\pm |H|\mp ⟩=\frac{1}{2}hA\sqrt{(I+1/2{)}^2-m_F^2}$

Resolviendo los valores propios de esta matriz, (como se puede hacer a mano o más fácilmente, con un sistema de álgebra por computadora) llegamos a los cambios de energía:

$\mathrm{\Delta }{E}_{F=I\pm 1/2}=-\frac{h\mathrm{\Delta }W}{2(2I+1)}+{\mu }_{\mathrm{N}}{g}_I{m}_{F}B\pm \frac{h\mathrm{\Delta }W}{2}\sqrt{1+\frac{2m_{F}x}{I+1/2}+x^2}$

$x\equiv \frac{B({\mu }_{\mathrm{B}}{g}_J-{\mu }_{\mathrm{N}}{g}_I)}{h\mathrm{\Delta }W}\mathrm{\Delta }W=A(I+\frac{1}{2})$

donde $\mathrm{\Delta }W$ es la división (en unidades de Hz) entre dos subniveles hiperfinos en ausencia de campo magnético $B$, $x$ se conoce como el parámetro de intensidad de campo (Nota: para $m_F=\pm (I+1/2)$ la expresión debajo de la raíz cuadrada es un cuadrado exacto, por lo que el último término debe reemplazarse por $+\frac{h\mathrm{\Delta }W}{2}(1\pm x)$). Esta ecuación se conoce como fórmula de Breit-Rabi y es útil para sistemas con un electrón de valencia en un nivel $s$ ($J=1/2$).

Tenga en cuenta que el índice $F$ en $\mathrm{\Delta }{E}_{F=I\pm 1/2}$ debe considerarse no como momento angular total del átomo, sino como momento angular total asintótico. Es igual al momento angular total solo si $B=0$ de lo contrario, los vectores propios correspondientes a diferentes valores propios del hamiltoniano son las superposiciones de estados con diferentes $F$ pero igual $m_F$ (las únicas excepciones son $|F=I+1/2,m_F=\pm F⟩$).

Aplicaciones

Astrofísica

Efecto Zeeman en una línea espectral de manchas solares

George Ellery Hale fue el primero en notar el efecto Zeeman en los espectros solares, lo que indica la existencia de fuertes campos magnéticos en las manchas solares. Estos campos pueden ser bastante altos, del orden de 0,1 tesla o más. El efecto Zeeman se utiliza para producir magnetogramas que muestran la variación del campo magnético del sol.

Refrigeración por láser

El efecto Zeeman se utiliza en muchas aplicaciones de enfriamiento láser, como una trampa magnetoóptica y el Ralentizador Zeeman.

Acoplamiento de espín y movimientos orbitales mediado por la energía de Zeeman

La interacción spin-órbita en los cristales generalmente se atribuye al acoplamiento de matrices de Pauli $\mathit{\sigma }$ al impulso de los electrones $𝒌$ que existe incluso en ausencia de campo magnético $𝑩$ . Sin embargo, bajo las condiciones del efecto Zeeman, cuando $𝑩\ne 0$, se puede lograr una interacción similar acoplando $\mathit{\sigma }$ a la coordenada del electrón $𝒓$ a través de la espacialmente no homogénea Zeeman Hamiltoniana

$H_{\mathrm{Z}}=\frac{1}{2}(𝑩\hat{g}\mathit{\sigma })$ ,

donde $\hat{g}$ es un factor g tensorial de Landé y $𝑩=𝑩(𝒓)$ o $\hat{g}=\hat{g}(𝒓)$, o ambos, dependen de la coordenada del electrón $𝒓$. Semejante $𝒓$ -dependiente al Hamiltoniano de Zeeman $H_{\mathrm{Z}}(𝒓)$ de las parejas de espín de electrones $\mathit{\sigma }$ al operador $𝒓$ que representa el movimiento orbital del electrón. El campo no homogéneo $𝑩(𝒓)$ puede ser un campo suave de fuentes externas o un campo magnético microscópico de oscilación rápida en antiferromagnetos. El acoplamiento spin-órbita a través de un campo macroscópicamente no homogéneo $𝑩(𝒓)$ de nanoimanes se utiliza para la operación eléctrica de espines de electrones en puntos cuánticos a través de la resonancia de espines dipolares eléctricos, y la conducción de espines mediante un campo eléctrico debido a $\hat{g}(𝒓)$ , también se ha demostrado.

Véase también

• Efecto Kerr magneto-óptico
• Efecto Voigt
• Efecto Faraday
• Efecto Cotton-Mouton
• Energía de Zeeman
• Efecto Stark
• Efecto Lamb
  • La configuración electrónica dice que en la subcapa p (l = 1), hay 3 niveles de energía ml = -1,0,1, pero solo vemos dos p1/2 y p3/2, para la subcapa s (l = 0), solo hay 1 nivel de energía (ml = 0), pero aquí tenemos 2. l correspondiente a estructura fina, ml correspondiente a estructura hiperfina.

Referencias

Histórico

• Plantilla:Cita libro (Chapter 16 provides a comprehensive treatment, as of 1935.)
• Plantilla:Cita publicación
• Plantilla:Cita publicación
• Plantilla:Cita publicación
• Plantilla:Cita publicación
• Plantilla:Cita publicación

Moderno

• Plantilla:Cita libro
• Plantilla:Cita publicación
• Plantilla:Cita libro
• Plantilla:Cita libro
• Plantilla:Cita libro
• Plantilla:Cita libro