Capítulo Nº2 Correlaciones angulares

Capítulo No2

Correlaciones angulares

En el presente Capítulo se realiza una descripción sintetizada de la teoría de las

Correlaciones Angulares Perturbadas Diferenciales en Tiempo. Se aborda el caso
particular de la interacción cuadrupolar eléctrica entre el momento cuadrupolar del
núcleo y el Gradiente de Campo Eléctrico del campo extranuclear y por último se
describe la metodología experimental utilizada en el presente trabajo. Para un análisis
más exhaustivo del tema se recomienda ir al Apénice I y al trabajo [Frauenfelder, 1968].

2-1 Función de correlación

La desexcitación de un isótopo radiactivo debida a la emisión consecutiva de dos
radiaciones  está caracterizada por tres estados nucleares: un estado inicial de momento
angular Ii, un estado intermedio I (el cual tiene una vida media  y un momento
cuadrupolar Q) y un estado final If. Los núcleos decaen del estado Ii al estado I
emitiendo la primera radiación 1 y acceden al estado final emitiendo una segunda
radiación 2 (ver, por ejemplo, la sección 2 del Apéndice I). Este proceso de denomina
“cascada -” y se esquematiza en la Figura 2.1.

Ii

1
I, , Q

2
If

Figura 2.1: Esquema de desexcitación

de un núcleo radiactivo a través de una
cascada -.

7
Tesis Doctoral Emiliano L. Muñoz

Entre las direcciones de emisión k1 y k2 de las radiaciones 1 y 2, respectivamente,

existe una correlación debido a la conservación del momento angular del sistema. La
orientación del espín nuclear del estado intermedio I está relacionada con la dirección k1
y puede expresarse en términos del número cuántico momento angular magnético m
respecto a un eje de cuantificación arbitrario que en nuestro caso será k1.
Debido a que una muestra radiactiva está compuesta por un gran número de núcleos
realizando transiciones idénticas y teniendo en cuenta que sus espines se encuentran
orientados aleatoriamente, entonces no habrá una dirección preferencial de emisión por
lo que el patrón deviene isótropo. Para “romper” esta isotropía debemos lograr que los
núcleos estén orientados en alguna dirección preferencial. La detección de 1 en una
dirección k1 selecciona un subconjunto de núcleos particularmente orientados, por lo
que la radiación 2 detectada en coincidencia con 1 tendrá una distribución angular
anisotrópica definida.
En presencia de campo extranuclear, por ejemplo el campo generado en un sólido
cristalino, los momentos nucleares interactuarán con dicho campo. Esta interacción se
conoce como interacción hiperfina y es capaz de producir (semi-clásicamente) la
precesión de los espines nucleares alrededor de alguna dirección preferencial,
cambiando el patrón de emisión. Si el campo cristalino solo tiene una contribución
eléctrica, entonces, el gradiente del campo eléctrico (GCE) interactúa con el momento
cuadrupolar eléctrico nuclear. Para una diferencia de tiempo nula entre las emisiones de
1 y 2, la diferencia de contaje entre las coincidencias medidas entre el detector 2
(donde se detecta 2) y el detector 1 (donde se detecta 1) posicionados a 90º y 180º nos
da la anisotropía (ver Figura 2.2). En el caso de no tener una emisión simultánea de las
radiaciones, la anisotropía será función del tiempo transcurrido entre la emisión de 1 y
de 2 debido a que durante este tiempo el espín del estado intermedio nuclear precesa
alrededor del eje de cuantificación. La anisotropía se verá recuperada para un tiempo
t=2/, donde  es la frecuencia de precesión de espín. Esta frecuencia es
proporcional a la componente principal del tensor GCE y está relacionada con la
separación de los subestados m del nivel intermedio I cuando el campo externo rompe
su degeneración.

8
Capítulo Nº2 Correlaciones angulares

Figura 2.2: Esquema semi-clásico del principio básico de la

determinación de la frecuencia de precesión del espín nuclear
debido a la presencia de un campo extranuclear.

La relación entre la frecuencia de precesión  y el desdoblamiento E se conoce

como precesión de Larmor y es  = E/ħ. La Figura 2.3 muestra el desdoblamiento de
un estado nuclear para el caso particular de I=5/2 y sus energía en función de la
frecuencia de Larmor.

Figura 2.3: Esquema de decaimiento nuclear de un

núcleo a través de la emisión consecutiva de dos
radiaciones 1 y 2 y desdoblamiento del nivel
intermedio I de la cascada -.

9
Tesis Doctoral Emiliano L. Muñoz

La probabilidad de que 2 sea emitido en la dirección k2 en una abertura d2 al

tiempo t+dt, si 1 ha sido emitido en la dirección k1 en el instante t en una abertura d1
es:

   
P ( k1 , k 2 , t )  W ( k1 , k 2 , t ).e  t /  d  1d  2 dt , (2.1)

 
donde W (k1 , k 2 , t ) se define como la función de correlación angular y es la probabilidad
por unidad de ángulo sólido d1 y d2, y por unidad de tiempo de que 1 se emita en la
dirección k1 al tiempo t en correlación con la emisión de 2 en la dirección k2 al tiempo
t+dt.
Como se demuestra en el Apéndice I, la función de correlación no perturbada puede
expresarse como función del ángulo entre las direcciones de emisión  como sigue:

K máx
(2.2)
W ( )  A
K 0
KK PK (cos( )) ,

donde K=0, 2, 4, …, Kmáx y Kmáx ={2I, 2L1, 2L2}, Li es el momento angular de i.
PK(cos()) es el K-ésimo polinomio de Legendre, y los coeficientes AKK, denominados
coeficientes direccionales o coeficientes de anisotropía, miden la desviación respecto
del caso isotrópico (W()=1). La normalización de esta expresión se realiza tomando
A00=1 y como puede apreciarse en (2.2), para obtener un patrón anisotrópico se
requiere I1. Cabe destacar que debido a que la probabilidad de transición entre estados
nucleares decrece con el incremento del momento angular de la radiación emitida, en
general Kmáx no es mayor que 4.

2-2 Factor de Perturbación

Si consideramos el efecto de los campos externos al núcleo, entonces tendrá lugar el

desdoblamiento de los subestados m degenerados del estado intermedio de espín I.
Desde un punto de vista semi-clásico la interacción hiperfina entre el tensor GCE y el
momento cuadrupolar Q de un núcleo produce un movimiento de precesión del espín
nuclear alrededor del eje de cuantificación sobre el cual cada subestado tiene una
proyección m dando lugar a una re-orientación nuclear. El desdoblamiento típico de los

10
Capítulo Nº2 Correlaciones angulares

níveles de energía E es del orden de 10-8-10-6 eV, dependiendo de la intensidad del

campo. Este desdoblamiento es varios ordenes de magnitud menor que las energías
involucradas en las transiciones nucleares que son del orden de 102 keV.
Finalmente, estos efectos perturbativos modifican la orientación del espín nuclear y
cuantitativamente también se modifica la función de correlación angular, siendo ahora
función del ángulo  y del tiempo t (ver sección 1 del Apéndice I), tomando la forma:

K máx
W ( , t )  A
K 0
KK GKK (t ) PK (cos( )), (2.3)

donde GKK(t) se denomina Factor de Perturbación y contiene toda la información de la

interacción hiperfina entre el tensor GCE del campo extranuclear y el momento
cuadrupolar del núcleo.
La forma explícita del factor de perturbación depende del tipo de interacción
hiperfina experimentada por el núcleo (por ejemplo interacción dipolar magnética o
interacción cuadrupolar eléctrica), también depende de la dependencia temporal de
dicha interacción (si es estática o dinámica), como así también del tipo de muestra
(sistemas en estado bulk, monocristales, nanopartículas, etc). A modo de ejemplo,
mostramos el caso del factor de perturbación para una muestra policristalina donde el
núcleo experimenta una interacción cuadrupolar eléctrica con un campo externo
estático, esto es, independiente del tiempo.

nmáx
G KK (t )   S Kn cos( n t ) . (2.4)
n 0

Los coeficientes SKn y las frecuencias n dependen del momento cuadrupolar Q del
núcleo sonda, del valor del espín nuclear I del estado intermedio y de las componentes
principales del tensor GCE observado en el sitio de la sonda. Debido a que el núcleo a
tiempo t=0 aún no ha sido perturbado (GKK(t=0)=1), entonces los coeficientes SKn
nmáx
deben satisfacer S
n 0
Kn  1.

11
Tesis Doctoral Emiliano L. Muñoz

2-3 Interacción cuadrupolar eléctrica

En esta sección vamos a analizar el caso particular de la interacción entre el núcleo

y el campo eléctrico del campo cristalino del sólido donde se encuentra alojada la
sonda, ya que es la interacción de interés en esta Tesis. Supongamos que el campo
magnético del campo extranuclear es lo suficientemente débil como para despreciar la
interacción magnética frente a la interacción eléctrica. Entonces, la principal interacción
entre el núcleo de carga Ze y el campo cristalino será de tipo Coulombiana,
cuantitativamente dada por ZeV(0) (donde se supone que el núcleo se sitúa en el origen).
Debido a que la distribución de carga nuclear no es esférica y teniendo en cuenta que el
campo eléctrico cristalino no es homogéneo aparecerán términos correctivos en la
energía, rompiéndose así la degeneración y produciendo el desdoblamiento de los
subestados m del nivel intermedio I. La energía electrostática del sistema debido a la
interacción entre el núcleo y el potencial Vext del campo cristalino es la convolución de
la densidad de carga nuclear con dicho potencial:

 
Eelec    n ( r )Vext ( r ) d 3 r . (2.5)

El potencial externo Vext(r) puede desarrollarse en serie de Taylor [Jackson 1964] de

modo que la ecuación anterior queda expresada en una serie donde cada término de
orden n se relaciona con la multipolaridad de la interacción:

 1  E j  0 3  
Eelec  qV  0  p.E  Qij    el  0 n  r  r 2dr , (2.6)
6 i j  xi 2


donde q es la carga total del núcleo y representa el momento monopolar, p es el
momento dipolar eléctrico, Qij es el momento cuadrupolar eléctrico nuclear de traza
nula, y el y n son las densidades de carga electrónica y nuclear, respectivamente.
En la ecuación (2.6) se puede apreciar que la carga nuclear se acopla con el
potencial electrostático; el momento dipolar eléctrico, con el campo eléctrico; el
momento cuadrupolar, con el GCE; y el último término es la contribución de los
electrones s que tienen una probabilidad no nula de penetrar en la región nuclear.

12
Capítulo Nº2 Correlaciones angulares

Debido a que la simetría de la densidad carga de electrones s es esférica, el término es

de orden cero y, al igual que la interacción entre el momento monopolar y el potencial
V(0), sólo produce un corrimiento en la energía electrostática de la interacción. La
interacción dipolar se anula debido a la simetría de las funciones de onda nucleares. Por
lo tanto, el único término que rompe la degeneración del estado intermedio I es el
correspondiente a la interacción entre el momento cuadrupolar eléctrico nuclear y el
tensor GCE.
Como ya se mencionó anteriormente, la interacción entre el GCE y el momento
cuadrupolar Qij del núcleo en el estado intermedio I, semiclásicamente provoca la
precesión del espín nuclear entorno a la dirección del eje de cuantificación.
El tensor GCE queda definido a partir de las nueve componentes
 2V 
Vij  , i, j  1,2,3 . V( r ) es continuo de orden dos, por lo que la matriz real de 3x3
xi x j

del tensor GCE es simétrica y por ende diagonalizable, reduciéndose a las tres
componentes principales V11, V22 y V33. Estas tres componentes se definen de manera
que V33  V22  V11 y debido a que satisfacen la ecuación de Laplace

(V11+V22+V33=0, tensor de traza nula), el tensor GCE diagonalizado queda

completamente definido por dos de sus elementos diagonales. Es conveniente definir un
parámetro adimensional, el parámetro de asimetría , que depende de las tres
componentes principales:

V11  V22
 (2.7)
V33

De este modo el GCE queda completamente determinado a partir de  y de la

componente mayor V33.
Este parámetro  expresa la desviación del tensor GCE respecto de aquel que posee
simetría axial ( = 0). A partir de la definición de los elementos diagonales del tensor y
de , se tiene 0    1 .
El hamiltoniano perturbativo que describe la interacción entre el GCE del campo
extranuclear y el momento cuadrupolar del núcleo en el estado intermedio I viene dado
por [Frauenfelder, 1968]:

13
Tesis Doctoral Emiliano L. Muñoz

  K  
He    QKm VKm , (2.8)
K  0 m  K

 
donde QKm es el operador momento cuadrupolar eléctrico nuclear y VKm es el operador
tensor gradiente de campo eléctrico del campo extranuclear. Debido a que trabajamos
en la base Km , ambos operadores se expresan en coordenadas esféricas, esto es

r Km  YmK ( ,  ) . De esta manera, se puede demostrar [Rentería, 1992] que los

autovalores del hamiltoniano perturbativo de interacción cuadrupolar eléctrica H e , para
I=5/2, toman la forma:

E 5 2  1
3
Q cos  13 arccos  

E 3 2   1
3
Q cos 13   arccos   (2.9)

E1 2   1
3
Q cos 13   arccos   ,

eQV33
donde Q  es la frecuencia cuadrupolar,   4 7 3   2  y
40
  1920 3 1   2  /  3 .

Debido a que los autovalores de H e no dependen del signo de m, esto es, del signo
de la proyección del espín nuclear I en el eje de cuantificación, entonces las frecuencias
Em  Em '
de interacción   , tampoco dependerán de dicho signo, por lo que solo habrá

tres frecuencias de interacción i:

1  1
 E  32 
 E  1 2   Q  sen  13 arccos  

2  1
 E 52 
 E 3 2  Q  sen  13   arccos    (2.10)

3  1
 E  52 
 E  1 2   Q  sen  13   arccos   

14
 Capítulo Nº2 Correlaciones angulares

En la Figura 2.4 se grafican los autovalores de energía, en unidades de Q , en

función del parámetro de asimetría .

E /  Q Observando las energías (2.9) del

12
m =  5 /2 desdoblamiento del estado intermedio nuclear de
8
espín I=5/2, vemos que la única información de
4
2
la interacción cuadrupolar eléctrica necesaria
0
m =  3 /2

3 para calcular dichos autovalores reside en  y

-4
m =  1 /2 1
-8
Q.
-1 2 Las frecuencias de interacción n suelen
0 .0 0 .2 0 .4 0 .6 0 .8 1 .0
expresarse como un producto de factores

independientes, esto es:
Figura 2.4: Autovalores de energía
E vs. . Se identifican las tres
frecuencias de interacción.
 n  g n ( ) Q  g n ( )Q (2.11)

con n=1, 2, 3.
En la Figura 2.5 se grafican las funciones gn, con n =1, 2, 3, en función de  y los
2 3
cocientes
1 y 1 , y también en función de :

gn() 3.5
g () 3.0
 /
3
3.0
3 1
2.5

2.5
g () 2.0
2
2.0
1.5
1.5  /
2 1
1.0
1.0 g ()
1
0.5
0.5

0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0

 
Figura 2.5: Gráfico de gn() (n=1,2,3) vs.  y de
los cocientes 2/1 y 3/1en función de .

15
Tesis Doctoral Emiliano L. Muñoz

A modo de ejemplo veamos que forma tienen los desdoblamientos de energías y las
frecuencias de transición para el caso de GCE axialmente simétrico, esto es =0. El
desdoblamiento de energía entre dos subestados arbitrarios m y m’ es de la forma:

( m 2  m' 2 )
E m, m '  E m  E m '  3eQV33 , (2.12)
40

y explícitamente:

 E  3 / 2 , 1 / 2  6 Q    1  1  6 Q

 E  5 / 2 , 3 / 2  12 Q     2   2  12 Q  21 (2.13)

 E  5 / 2 , 1 / 2  18 Q    3   3  18 Q  31 .

Para las sondas empleadas en este trabajo el coeficiente A44 (ver ecuación 2.3) es del
orden de 1% con respecto al coeficiente A22. Por lo que la función de correlación toma
la forma:

W ( )  1  A22 G22 (t ) P2 (cos( )) . (2.15)

El factor de perturbación G22(t) puede escribirse como:

3
G22 (t )  S 20 ( )   S 2 n ( ) cos(n ( , Q )t ) , (2.16)
n 1

3
donde los coeficientes S2n están normalizados (  S 2n 1) . De (2.16) vemos que G22(t)
n 0

depende únicamente de  y Q, es decir, de los parámetros que caracterizan el tensor

GCE. Finalmente, vía una medida experimental de la correlación angular (perturbada)
de las emisiones 1 y 2 en la desexcitación del núcleo de interés podremos determinar 

16
Capítulo Nº2 Correlaciones angulares

y Q, y si conocemos Q podremos obtener una medida de la componente principal del

tensor GCE,V33, experimentada por el núcleo elegido.

2-4 Método Experimental

La técnica de las Correlaciones Angulares Perturbadas Diferenciales en Tiempo
(PAC) se basa en la determinación de la interacción hiperfina entre un núcleo radiactivo
y el campo extranuclear del sistema actuante en el sitio donde se aloja el núcleo a partir
de la determinación de la correlación angular perturbada entre las radiaciones 1 y 2
emitidas por el núcleo radiactivo. Entonces, estos núcleos inmersos en un sólido
devienen en sondas adecuadas para estudiar propiedades del sistema huésped. Es por
esto que cada uno de estos núcleos suele llamarse núcleo-sonda.

2-4-1 Medida de las Correlaciones Angulares Perturbadas Diferenciales en Tiempo

Para determinar experimentalmente el tensor GCE, conocido el momento
cuadrupolar nuclear Q del núcleo-sonda, es necesario obtener el factor de Perturbación
G22(t) presentado en (2.16). Para obtenerlo necesitamos medir el número de
coincidencias 1-2 en función del tiempo t transcurrido entre ambas emisiones,
detectadas en dos posiciones angulares diferentes. Si llamamos Cexp(,t) a las
coincidencias medidas para un ángulo  entre las radiaciones emitidas y luego de haber
transcurrido un tiempo t entre dichas radiaciones, se tiene:

C exp ( , t )  AP1 P2  1 2  CW exp ( , t )e t /   C A , (2.17)

donde A es la actividad de la muestra, P1 es la probabilidad de emisión de 1 por

cada desintegración del núcleo sonda, P2 es la probabilidad de emisión de 2 habiéndose
emitido 1, 1 y 2 son las eficiencias absolutas de los cristales para detectar 1 y 2,
respectivamente, c es la eficiencia del circuito de coincidencias y  es la vida media del
estado intermedio de la cascada -. Las coincidencias accidentales CA son aquellas en
las que 1 y 2 provienen de diferentes núcleos. Dichas coincidencias son proporcionales
al cuadrado de la actividad de la muestra y a la resolución temporal del equipo. En la
Figura 2.6 se muestran a modo de ejemplo, espectros de coincidencias típicos. Las
cuentas a la izquierda del máximo son las coincidencias accidentales, las que se deben
descontar para determinar las coincidencias reales.

17
Tesis Doctoral Emiliano L. Muñoz

La dependencia temporal de la
20000 función de correlación experimental
Número de coincidencias 16000 W exp ( , t ) involucra tiempos de emisión
12000
 = 90° anteriores y posteriores a t, debido a que
8000 la resolución temporal del equipo
 = 180°
4000
Accidentales
experimental (r(t-t’)) impide medir la
2 /
0 correlación exactactamente en t. Por lo
0 10 20 30 40 50
Tiempo (ns) que dos eventos no simultáneos pueden

Figura 2.6: Espectros de coincidencias interpretarse como tales, o bien dos

1-2 para  = 90 y  = 180. Los sucesos consecutivos producidos con un
espectros corresponden a 181Ta en C-
retardo dado son detectados con un
Ho2O3 [Errico 2002].
retardo diferente. Entonces, la función
W exp ( , t ) se puede describir a partir de la convolución temporal de la función de
correlación perturbada teórica W ( , t ) (ver la expresión (2.15)) con la resolución

temporal finita r (t  t ' ) del equipo experimental y el decaimiento exponencial e t /  del

estado intermedio nuclear I:

 r (t  t ' )W ( , t ' )e
t ' / 
dt '
W exp ( , t )  0
. (2.18)

 r (t  t ' )e
t ' / 
dt '

A partir de la expresión (1.14) para la función de correlación obtenemos:

W exp ( , t )  1  A22 G22 (t ) P2 (cos( )) ,

exp exp
(2.19)

donde hemos definido el factor de perturbación experimental G22exp(t) como:

 r (t  t ' )G 22 (t ' )e t '/  dt '

(2.20)
exp
G22 (t )  0
.
 r (t  t ' )e
t '/ 
dt '

18
Capítulo Nº2 Correlaciones angulares

La resolución r(t-t’) puede determinarse experimentalmente midiendo la respuesta

temporal del equipo ante una excitación del tipo “delta”. En la práctica esto se puede
realizar, por ejemplo, detectando la coincidencia de los rayos  de 1100 keV y 1300 keV
de la cascada - del 60Co, donde la vida media del estado intermedio es de 0.73 ps, la
cual resulta despreciable frente a la resolución típica de un equipo PAC (10-1 ns). Un
método alternativo es medir las dos radiaciones  de 511 keV emitidas en la
aniquilación positrón-electrón posterior al decaimiento + del 22
Na. Debido a la
conservación del momento y la energía, son emitidos en forma simultánea dos fotones
de 511 keV en direcciones opuestas con una diferencia temporal nula [Leo, 1994].
A modo de ejemplo, en la Figura 2.7 se muestra el espectro de coincidencias de la
cascada - del 60Co medido en el espectrómetro PACar de La Plata [Rentería 2008] con
las ventanas de energía de la cascada - de la sonda (181Hf→)181Ta, el tiempo de
resolución del equipo 2res es el ancho a mitad de altura.

Figura 2.7: Espectro de coincidencias de la

cascada 1100 – 1300 keV del 60Co y
resolución temporal del equipo PAC
[Rentería, 2008].

De la expresión (2.17) puede apreciarse que para cada tiempo t puede obtenerse un
sistema de dos ecuaciones con dos incógnitas (A22exp y G22exp), de modo que midiendo
las coincidencias en dos posiciones diferentes de los detectores se puede determinar
G22exp(t). Comúnmente se eligen las posiciones =90º y 180º puesto que la distribución
angular de la radiación del núcleo-sonda presenta la mayor diferencia relativa de

19
Tesis Doctoral Emiliano L. Muñoz

intensidad en tales ángulos. Se define la función R(t), llamada espectro de precesión de

espín, relación de asimetría o simplemente espectro PAC, como:

C 180  , t   C  90  , t  (2.21)
R t   2 ,
C 180 , t   2C  90 , t 
 

donde C(,t)=Cexp(,t)-CA.
Reemplazando por la forma explícita de Wexp(,t) (ec. 2.19) obtenemos:

exp exp exp exp

1  A22 G22 (t ) P2 (1)  (1  A22 G22 (t ) P2 (0)) (2.22)
R(t )  2 exp exp exp exp .
1  A22 G22 (t ) P2 (1)  2(1  A22 G22 (t ) P2 (0))

1
A partir de la expresión del polinomio de Legendre de orden 2 ( P2 ( x)  (3x 2  1) ),
2

tenemos que P2(-1)=1 y P2(0)=-1/2, entonces:

3
R(t )  exp exp
A22 G 22 (t )  exp
A22 S
n 0
2 n ( ) cos( n t ) ,
(2.23)

donde 0=0.
El espectro PAC definido según (2.21) nos permite obtener G22exp(t), siendo
conocido A22. Cabe destacar que la forma del espectro R(t) es válida para arreglos
experimentales de dos, tres o cuatro detectores.
En la Figura 2.8 se simulan distintos espectros R(t) para I=5/2 y distintos valores de
, con sus respectivas transformadas de Fourier donde se pueden apreciar las
frecuencias de interacción. Puede observarse que la anisotropía puede determinarse
midiendo la correlación angular para GCE nulo en el sitio del núcleo sonda (Figura 2.8
(a) izquierda) o del espectro R(t) a t=0.

20
 Capítulo Nº2 Correlaciones angulares

(a)
-0.2 0.08
-0.1 -0.2 (a) 1
0.04
0.0 2
-0.1 0.04
0.00 3
-0.2 (b) 0.08 0.0
1
2 0.00
-0.1
0.04 3 -0.2 (b)
0.0
-0.1 0.04
0.00

Amplitud (U. A.)

Amplitud (U. A.)

-0.2 (c) 0.08 0.0
0.00

R(t)
R(t)

-0.1
0.04 -0.2 (c)
0.0
-0.1 0.04
0.00
-0.2 (d) 0.08 0.0
-0.1 0.00
0.04 -0.2 (d)
0.0
-0.1 0.04
0.00
-0.2 (e) 0.08 0.0
-0.1 0.00
0.04 0 10 20 30 40 50 0 1 2 3 4
0.0 t (ns)  (Grad/s)
0.00
0 10 20 30 40 50 0 1 2 3 4
t (ns)  (Grad/s)

Figura 2.8: Espectros R(t) simulados para una única interacción caracterizada
por: Izquierda: (a) GCE nulo; (b)  = 0; (c)  = 0.3; (d)  = 0.7; (e)  = 1.
Para los casos (b)-(e) se consideró Q = 100 Mrad/s. Derecha: (a)  = 0; (b) 
= 0.05; (c)  = 0.10; (d)  = 0.15. En todos los casos (menos izquierda (a)) se
consideró Q = 100 Mrad/s y  = 0.

Resumiendo lo que hemos deducido hasta aquí, hemos encontrado una expresión
cuantitativa que da cuenta de la interacción cuadrupolar eléctrica entre el tensor GCE
del campo cristalino del sólido donde se aloja la sonda y el momento cuadrupolar
nuclear del núcleo-sonda en su estado intermedio de espín I. De esta forma,
encontramos una expresión de naturaleza experimental que nos permite ajustar los
espectros de rotación de espín experimentales y así obtener el tensor GCE en el sitio de
la sonda. Ahora, esta expresión experimental fue hallada a partir de varias hipótesis que
es conveniente detallar: la muestra y los detectores se han considerado puntuales, hemos
considerado que el sólido bajo estudio constituye una red cristalina perfecta y que las
sondas ocupan solamente sitios cristalográficos (propios de dicha red) libres de
defectos. Debido a estas idealizaciones es necesario introducir en el factor de
perturbación determinadas aproximaciones que nos acerquen de una manera más precisa
al factor real:
- Como dijimos anteriormente, para el cálculo de la función de correlación
angular W(,t) hemos supuesto puntuales los cristales centelladores de los detectores y
la fuente radiactiva. Su tamaño finito y por ende la abertura de ángulo sólido no

21
Tesis Doctoral Emiliano L. Muñoz

infinitesimal subtendida desde la muestra hasta los detectores provoca la atenuación de

la anisotropía A22, por lo que Wexp(,t) tendrá una amplitud menor que la teórica. Para
reducir estas anomalías se utilizan fuentes lo más puntuales posible y se trata de reducir
la abertura del ángulo sólido aumentando la distancia muestra-detector, pero teniendo en
cuenta que el número de eventos detectados decrece con el cuadrado de dicha distancia.
- Un cristal real no puede considerarse como una red cristalina perfecta sin
imperfecciones como ser impurezas, dislocaciones de planos, defectos puntuales como
vacancias, átomos en los intersticios, daños por la radiación, etc. Estas imperfecciones
producen una distribución finita  en las frecuencias cuadrupolares n alrededor de sus
valores medios. Dicha distribución puede suponerse una distribución gaussiana o
lorentziana P(n). En la Figura 2.8 (derecha) se aprecia como estas distribuciones
atenúan las oscilaciones en el espectro R(t). Por lo tanto, el factor de perturbación
teórico debe ser pesado con las distribuciones P(n) en la forma:


G22 (t )   P( )G

22 (t )d . (2.24)

- Las sondas pueden ubicarse en sitios con entornos cristalográficos distintos, esto
es, en sitios no equivalentes, aún para el caso de una red cristalina perfecta. Entonces, el
factor de perturbación G22(t) más general que podamos tener será una combinación
lineal de la forma: G22 (t )  f 0   f i G22
i
(t ) , donde fi es la fracción de núcleos sonda
i

que experimentan una interacción hiperfina caracterizada por los parámetros Qi y i en
el sitio inequivalente i y f0 es la fracción de núcleos sonda que experimenta un GCE
nulo. Dichas fracciones deben cumplir f 0   fi  1
i

- Suponiendo que la respuesta temporal del equipo experimental (ver Figura 2.7)
es gaussiana, y teniendo en cuenta que una función gaussiana y su transformada de
Fourier tienen el mismo aspecto funcional, por lo que la convolución de W(,t) con
dicha respuesta nos da Wexp(,t) multiplicado por tal distribución. Esto nos lleva a
introducir factores reales del tipo:

1 2 2
  res n
S real
2n  S 2 n ( ).e 2 (2.25)

22
 Capítulo Nº2 Correlaciones angulares

Como ejemplo, en la Figura 2.9 se muestra como varía el espectro R(t) con la
resolución temporal, a través de dos espectros R(t) para la misma muestra, K2ZrF6
181
dopado con Ta, medida con dos detectores

R (t) a de resolución temporal muy diferente.

Claramente puede apreciarse que las
oscilaciones del espectro R(t) tienden a
t (n s) desaparecer cuanto menor sea la resolución
temporal.

R (t) b
Por lo tanto, considerando todas las
modificaciones introducidas tendientes a
acercar la expresión del R(t) a su forma real,
t (n s)
dicha expresión queda:
Figura 2.9: Espectros R(t) de 181Ta en
K2ZrF6 tomados con equipos con diferente
exp  
resolución temporal: a) CsF (~ 0,8 ns); b) R (t )  A22  f 0   f i G22
i
(t )  . (2.26)
NaI(Tl) (~ 3 ns) [Martinez, 1981].  i 

El factor de perturbación para cada sitio inequivalente i tiene la forma:

3
S  ( ni  Qi t ) p / p
i
G22 (t )  i i i
2 n ( ) cos( n t )e
(2.27)
n 0

El factor exponencial proviene de convolucionar G22(t) con la distribución de

frecuencias P(). Si dicha distribución es lorentziana (p=1), P() toma la forma:

 (2.28)
P ( )  ,
    
2

  
2

 2 
 

y si es gaussiana (p=2), se tiene:

1
P    . e     / 2
2 2
(2.29)
 2. .

23
Tesis Doctoral Emiliano L. Muñoz

2-4-2 Equipo experimental

La determinación experimental de G22(t) da como resultado la determinación de los
parámetros hiperfinos característicos de la interacción cuadrupolar eléctrica entre el
momento cuadrupolar Q del núcleo-sonda y el tensor GCE del campo extranuclear del
sólido. Para esto es necesario medir las coincidencias C(,t) de la cascada - del núcleo
para dos ángulos distintos. Es decir, medir 1 en el detector 1 en coincidencia con la
detección de 2 en el detector 2. Tal medición está caracterizada por la energía de ambas
radiaciones y por la determinación precisa de la diferencia temporal entre la emisión de
1 y 2.
La mayor parte de los experimentos PAC cuyos resultados se describen en el
Capítulo 5 fueron desarrollados durante mi estadía (Julio – Agosto de 2009) en el
Instituto de Pesquisas Energéticas y Nucleares (IPEN) de la Universidad de San Pablo
utilizando el espectrómetro PAC de cuatro detectores perteneciente al grupo de
Interacciones Hiperfinas del IPEN. El primero de los experimentos en SnO fue
realizado con el mismo espectrómetro 2 años antes. El exhaustivo análisis de sus
espectros y el desarrollo preliminar del modelo que describe estos resultados a partir de
cálculos de primeros principios precedió y preparó el desarrollo del resto de los
experimentos y la elección de los sistemas estudiados. Medidas preliminares de óxidos
binarios dopados con la sonda (111In→)111Cd se realizaron en el equipo PACar de lógica
rápida-rápida [Rentería 2008] (optimizado para ser utilizado con la sonda
(181Hf→)181Ta), pero la aparición de oscilaciones espurias en el equipo para las energías
de la cascada γ-γ del 111Cd (ver Fig. 2.10) imposibilitaba una correcta determinación de
los espectros. En consecuencia, decidimos utilizar un equipo de lógica rápida-lenta y de
111
mayor eficiencia, óptimo para ser utilizado con el isótopo In(→111Cd), el cual posee
una vida media mucho menor que el 181Hf(→181Ta). Finalmente, el análisis de todos los
espectros PAC fue realizado en La Plata.
Los detectores están compuestos por cristales centelladores de Floruro de Bario
(BaF2). La forma de los cristales del tipo cono truncado (ver. Fig. 2.9). La forma cónica
de los centelladores mejora la respuesta temporal de los detectores minimizando el
tiempo de colección de luz dentro del cristal centellador [Butz, 1989]. Cada centellador
está acoplado ópticamente a un fotomultiplicador. La señal de cada detector es
procesada de dos maneras diferentes independientes entre sí, en un circuito de
electrónica rápida-lenta como se muestra en la Figura 2.9. Por cada radiación gamma

24
Capítulo Nº2 Correlaciones angulares

detectada se producen dos pulsos de tensión, uno proveniente de la salida de dinodo del
fotomuliplicador (rama lenta) cuya altura del pulso es proporcional a la energía del
fotón; y el otro el otro proveniente de la salida de ánodo (rama rápida), saturado en
altura y de rápido crecimiento inicial, por lo que se lo utiliza para indicar el momento de
detección de la radiación.
En la rama lenta, los pulsos son conformados por el Pre-Amplificador (PA) y
amplificados por el Amplificador (AMP). Un analizador Mono-Canal (SCA) selecciona
las energías deseadas en cada detector, produciendo pulsos lógicos que alimentan el
Circuito de Coincidencias Lento representado por la Unidad de Direccionamiento
(Rounting) [Domienikan 2001]. Observar que los SCA están definidos con un número i
que depende del detector al cual se encuentra conectado (1, 2, 3 o 4), y también se
definen con una letra a o b que esquematiza si se selecciona 1 (a) o 2 (b).
En la rama rápida, sólo se analiza el tiempo entre los dos pulsos,
independientemente de su energía. A partir de un Discriminador a Fracción Constante
(CFD) se eliminan las señales espureas de los pulsos de ánodo y por medio de sus
tiempos de crecimiento se determina el instante en que han llegado al detector,
generando un pulso lógico. Dado un par de detectores que intervendrá en una
coincidencia, uno de ellos provee el pulso lógico de comienzo de la coincidencia (start)
y el otro el de finalización de la misma (stop). Estos pulsos alimentan al Convertidor de
Tiempo en Amplitud (TAC), el cual genera un pulso de altura proporcional a la
diferencia de tiempo de llegada de los pulsos start y stop.
El espectrómetro utilizado es capaz de analizar los doce posibles espectros de
coincidencias que se pueden obtener con un arreglo coplanar de cuatro detectores. Para
un mismo detector, la salida de cada SCA (rama lenta) está conectado con la salida del
contador CFD (rama rápida) del detector a una compuerta AND. Cada salida de ánodo
(rama rápida) de los detectores presenta un retardo (D) para producir el retraso
adecuado para que la señal de la rama rápida y la de la rama lenta lleguen
correlacionadas a la compuerta AND. De este modo se tienen ocho compuertas AND
debido a los cuatro detectores y a las dos radiaciones 1 y 2 detectadas en cada uno.
Las salidas de las compuertas AND correspondientes al 1 se conectan a una compuerta
OR (OR start) y las cuatro salidas AND del 2 se conectan a otra compuerta OR (OR
stop). La salida de cada compuerta OR se encuentra conectada a un retardo (D) que
tiene como función correr el cero de tiempos y permitir detectar las coincidencias

25
Tesis Doctoral Emiliano L. Muñoz

accidentales en una región pequeña de canales a la izquierda del cero de las

coincidencias verdaderas (“tiempos negativos”, correspondientes a algún 2 proveniente
de un núcleo diferente al que emitió el 1 y que llega antes que 1) para poder restarlas
tal como vimos en la sección 2-4-1.

Figura 2.9: Esquema del equipo PAC de cuatro detectores utilizado para
medir los sistemas estudiados.

26
Capítulo Nº2 Correlaciones angulares

Las salidas start (stop) de la compuerta OR start (OR stop) son direccionadas al
TAC. Por cada pulso start (stop) que arriba al TAC, este produce un pulso de validación
true start (true stop) que es enviado al Routing para validar los pulsos que llegan a este
desde la rama lenta. El routing entrega un pulso de compuerta (gate) cada vez que se
produce una coincidencia de las cuatro posibles. La salida del TAC posee un retardo
variable para garantizar que el gate esté presente cuando los pulsos lógicos del TAC
lleguen al Convertidor Analógico-Digital (ADC) para luego ser digitalizados. Los
pulsos digitalizados por el ADC en presencia del gate son enviados a la memoria
correspondiente de las doce existentes del Analizador Multicanal (MCA) según la
combinación (de las doce posibles) que ha detectado el Routing.

2-4-3 La sonda PAC (111In(EC)→) 111Cd

Las sondas radiactivas más utilizadas en espectroscopía PAC son (181Hf→)181Ta y

(111In→)111Cd. Ambas tienen el espín del estado intermedio I=5/2 y presentan
características compatibles para realizar medidas PAC en cualquier tipo de sistema, las
cuales se pueden resumir según:
- Parámetros nucleares conocidos.
- Adecuada anisotropía para la cascada - donde se desea realizar la medida PAC.
- Energía de las radiaciones  elevadas de modo de evitar la absorción en el material
donde se aloja la sonda o donde se sostiene la muestra.
- La energía entre 1 y 2 es mayor que la resolución en energía de los detectores
utilizados.
- Momento cuadrupolar nuclear Q alto de manera que la frecuencia Q sea observable
en el rango de tiempos dado por la vida media del estado intermedio.
- Vida media  del estado intermedio acotado por T=2/Q de manera tal que T< para
poder determinar con precisión Q.
- Frecuencias  limitadas por  y la resolución temporal res de manera que pueda
medirse  si cumple 2/5 <  < 2/res. El límite superior está dado por res y el
inferior por 5 (suponiendo que es el máximo tiempo que puede medirse
estadísticamente).

27
 Tesis Doctoral Emiliano L. Muñoz

En la Tabla 2.1 mostramos los parámetros nucleares de la sonda (111In→)111Cd que

es la utilizada en el presente trabajo, como así también los parámetros de la sonda
(111mCd→)111Cd debido que es el caso de la misma sonda producida desde diferente
padre radiactivo, esto tiene incidencia directa en los sistemas a estudiar como veremos
más adelante. En la Figura 2.10 mostramos el esquema de decaimiento del
(111In→)111Cd y del (111mCd→)111Cd.

Isótopo Vida Isótopo

1-2 (keV) I Q [b] (a)  (ns) A22 (%) (b) A44 (%) (c)
padre media hijo
111 111
In 3.6 d Cd 171-245 5/2 0.8313 84.1 -18.01 0.23
111m 111
Cd 48.6 m Cd 151-245 5/2 0.8313 84.1 17.51 -0.43

Tabla 2.1: Parámetros nucleares de las sondas (111mCd→)111Cd y (111In→)111Cd.

Figura 2.10: Esquemas de decaimiento de la sonda

(111In→)111Cd y de la sonda (111mCd→)111Cd.

(a)
[Herzog 1980].
(b)
Para ambas sondas [Butz 1983].
(c)
Para ambas sondas [Segeth 1987].

28
Capítulo Nº2 Correlaciones angulares

Referencias

[Jackson, 1964] J. D. Jackson, “Classical Electrodynamics”, John Wiley & Sons, New
York (1964).
[Frauenfelder, 1968] H. Frauenfelder y R. M. Steffen en “Alpha-Beta-and Gamma-Ray
Spectoscopy”, vol.2 (North- Holland, Amsterdam 1968) P.997.
[Rentería, 1992] M. Rentería, Tesis Doctoral, cap.2, Biblioteca del Departamento de
Física, Facultad de Ciencias Exactas, Universidad Nacional de La Plata (1992).
[Errico 2002] L.A. Errico Tesis Doctoral, cap.2, Biblioteca del Departamento de Física,
Facultad de Ciencias Exactas, Universidad Nacional de La Plata (2002).
[Leo, 1994] W. R. Leo, “Techniques for Nuclear and Particle Physics Experiments”, 2
edición, Springer Verlag, New York (1994).
[Rentería 2008] M. Rentería, A.G. Bibiloni, G.N. Darriba, L.A. Errico, E.L. Muñoz, D.
Richard y [Link], Hyp. Int. 181, 145 (2008).
[Martinez 1981] J.A. Martinez, A.M. Rodriguez, M.C. Caracoche, R.C. Mercader y A.
Lopez García, Hyp. Int. 11, 261 (1981).
[Butz, 1989] T. Butz, S. Saibene, Th. Fraenzke y M. Weber, Nuc. Instr. and Meth. A
284, 417 (1989).
[Domienikan 2001] C. Domienikan M. Sc. Thesis, IPEN, San Pablo (2001).
[Herzog 1980] P. Herzog, K. Freitag, M. Reuschenbach y H. Walitzki, Z. Phys. A 294,
13 (1980).
[Butz, 1983] T. Butz y A. Lerf, Phys. Lett. A. 97, 217, (1983).
[Segeth 1987] W. Segeth, Tesis Doctoral, Rijksuniversiteit Groningen, Holanda (1987).

29
Tesis Doctoral Emiliano L. Muñoz

30