Gradiente de Campo Elctrico en metales no cbicos

Gradiente de Campo Eléctrico en metales no cúbicos En este trabajo se comparan dos enfoques teóricos al estudio del Gradiente de Campo Eléctrico (GCE) en metales no cúbicos. El primer estudio fue publicado en 1975 y realiza un análisis sistemático teóricoexperimental basado en el modelo de cargas puntuales (MCP) y con datos reportados hasta esa fecha [1]. El segundo trabajo fue publicado en 1988 y versa sobre una comparación de resultados experimentales con cálculos de primeros principios [2]. EMILIANO L. MUNOZ

Outline Ø Modelo de Cargas Puntuales (MCP) Ø Método “Linearized Augmented Plane Wave” (LAPW) Ø Correlaciones Angulares Perturbadas (PAC) Ø Correlación Universal del GCE electrónico en metales nó cúbicos [1] iónico y Ø Cálculos de primeros principios del GCE en metales no cúbicos [2] Ø Comentarios finales

Modelo de Cargas Puntuales Sólido arreglo 3 D de iones puntuales en los sitios atómicos de red. Las componentes del tensor GCE se determinan según: Densidad de carga PCM: el tensor GCE en PCM toma la forma:

Para calcular el GCE a partir del modelo MCP debemos conocer los parámetros a, b, c; los ángulos , , y las posiciones de los átomos en la celda unidad. O Al Ø Generación de la red por traslación de celda hasta un radio R. Ø Cálculo de Vij en un contribuciones de c/ ión dado sitio la Ta de red, sumando Ø Diagonalización del tensor por medio de una subrutina, determinación de los ejes principales, Vii y . Contribución del carozo del átomo sonda: se supone proporcional a Vijred vía el factor de Sternheimer, el cual depende de cada átomo sonda y es calculado para átomos libres.

Cálculos ab initio: Teoría de la Funcional Densidad (DFT) Teoremas de Hohemberg y Kohn* E = E[ (r)] (r) real minimiza la E del sistema El éxito de la aplicación de esta teoría depende de la buena elección de (r) Aproximación de Born-Oppenheimer Energía de interacción Coulombiana electrón de Hartree Energía de intercambio de los electrones relacionada con la antisimetría de las funciones de onda electrónicas

Método Full-Potential Linearized Augmented Plane Wave (FP-LAPW) II I: funciones esféricas II: ondas planas I I II

Si (r) minimiza E[ (r)] Teoría variacional Ecuaciones de Kohn y Sham dada (r) Ciclo autoconsistente

Método experimental Técnica PAC: 1 I=5/2 Vzz Ii I 2 If 1 2 Detector 1 Orientación de espín al tiempo t I/ t Vzz Detector 2 Orientación de espín al tiempo t+dt

Tratamiento de datos: Correlaciones angulares Probabilidad de detección en coincidencia para correlación angular no perturbada Y en el caso de que la correlación sea perturbada, se tiene: es el ángulo entre las direcciones de emisión de los dos W( , t) es la Función de Correlación y tiene la forma: Z 1 A 22 es la Anisotropía y G 22(t) es la función de Perturbación y contiene toda la información de la interacción hiperfina entre el nucleo sonda y el campo extranuclear k 1 X k 2 2 Sistema de referencia del laboratorio Y

Para policristales el factor de perturbación viene dado por: Parámetros hiperfinos que se extraen de la medida PAC: donde Q es el momento cuadrupolar eléctrico del núcleo sonda. Parámetro de asimetría Frecuencia cuadrupolar Existe una relación unívoca entre , Q y las componentes principales del GCE Vii

Correlación angular, la medida: C( , t)exp es la coincidencia de detectar 1 y 2 formando un ángulo . Para determinar el espectro R(t) se miden las coincidencias a 90 o y 180 o, C(90 o, t) y C(180 o, t) respectivamente. El relación R(t) se define según: La función que se utiliza para fitear el espectro R(t)exp tiene la forma:

Correlación Universal del GCE Iónico y Electrónico en Metales no Cúbicos [1]

* eqlatt(1 - ) es la contribución iónica amplificada por el factor del átomo interactuante donde se calcula el gradiente. El segundo término es la contribución de los electrones de conducción, el cual es el efecto dominante y de carácter local. Tal término es corregido por la interacción con el core a partir del factor de apantallamiento 1 -R. Este paper plantea la cuestión: “El parámetro (1 -R) es adecuado para tener en cuenta la influencia del carozo electrónico de la sonda sobre la distribución espacial de los electrones de conducción en un metal? ” * R. E. Watson et. al. PRB 140, A 375 (1965).

En la siguiente figura se presenta una sistemática teóricoexperimental del GCE extraiónico VZZ’=Vzz - Vzzred (1 - ) en función del GCE iónico Vzzred (1 - ). Esta sistemática muestra una correlación universal VZZ’ con el GCE iónico y presenta las siguientes características: Ø La mayoría de los gradientes extrraiónicos están correlacionados sus contrapartes iónicas, siendo de signos opuestos entre si. Ø VZZ’ es proporcional al GCE iónico, Vzz’ ≈ -K Vzzred (1 - ) con K≈3; salvo para valores grandes del GCE iónico donde K se hace progresivamente más chico indicando la saturación del GCE extraiónico. Ø Para salvar desviaciones particulares de impurezas, redes huésped o combinaciones de ambas de la correlación universal es necesario adicionar un término de correción a la ecuación derivada de la sistemática. Vzz = Vzzred (1 - )(1 -K) + <Vzz>

Cálculos de primeros principios del GCE en metales no cúbicos [2]

Desarrollo del potencial II I II V( ) FP-LAPW I I II V( ) muffin-tin

El primer término considera el efecto de la densidad de carga dentro de la esfera atómica, mientras que el segundo y tercer término pesan la influencia de la densidad de carga fuera de la esfera atómica. El primer término se define como GCE de valencia por estar originado por la densidad electrónica no esférica de los electrones de valencia y semi-core. Los otros dos términos corresponden al GCE de red, el mismo incluye una sumatoria de multipolos. Esta suma es exacta pero es función de la elección del radio de la esfera. Tal elección debe satisfacer la no superposición de carga, esto es, de esferas atómicas.

Aspectos que afectan el cálculo del GCE: Ø Limitación numérica en el programa de cálculo de estructura de banda. El número del puntos k es el parámetro crítico del “volumen” de cálculo. Ø El tratamiento de los electrones de core. Se particionaron en estados de semi-core los estados 3 s y 3 p, los de la serie 3 d; 3 s, 4 s y 4 p, los de la serie 4 d. Estos son tratados como estados de bandas (no esféricos). Los estados de core son tratados como atómicos.

En el modelo de la ref. [1] se tiene: Vzz = Vzzred (1 - )(1 -K) + <Vzz> Este modelo fue utilizado para explicar el GCE en metales no cúbicos a partir de la “Correlación Universal” [1] basicamente a partir de la única determinación de Vzzred. La contribución de red es más pequeña que la local ya que está entre -10 y -80. El efecto de los electrones de conducción es mayor que el GCE iónico (K>1). En la tabla II se muestran todas las contribuciones, encontrándose que la mayor contribución a la local es la de valencia de carácter no esférico, mientras que la contribución de la parte esférica (core) es despreciable. En la tabla I se comparan valores del GCE experimentales y ab initio para metales no cúbicos.

Comentarios finales De la correlación universal del GCE extraiónico V’ZZ vs. VZZred(1 - ), donde V’ZZ = VZZ-VZZred(1 - ), se aprecia que para valores pequeños del GCE iónico, V’ZZ es proporcional (con pendiente K≈3) al GCE iónico y de signo opuesto. Para GCE’s iónicos grandes, K se hace progresivamente más pequeña. A esta correlación universal hace falta imponerle un término para salvar desviaciones particulares de alguna impureza y/o huésped (o combinación de ambos): VZZ = VZZred(1 - )(1 -K) + VZZ Por lo expresado en esta correlación universal concluimos que MCP vía este modelo considera la contribución electrónica al GCE a partir de un factor externo independiente de la simetría e intensidad del campo externo de la sonda. Además, debido a que se calcula para átomos aislados, no tiene en cuenta la interacción de la estructura electrónica de la sonda átomos de la red.

Ø El análisis ab initio sugiere fuertemente que el GCE está determinado por la densidad electrónica no esférica, predominantemente por los electrones p cercanos al núcleo. Ø Para obtener la densidad de carga fue necesario un cálculo del potencial en el marco del método LAPW teniendo en cuenta la contribución de todos los electrones del sistema. Es por esto que los métodos de seudopotenciales no deben ser utilizados para calcular el GCE en metales, támpoco es conveniente el cálculo vía potenciales de muffin-tin.

Ø Debido a que el GCE escalea con r -3, donde r es la distancia entre la carga que genera el GCE y el sitio sonda, el conocimiento exacto de la densidad de electrones cerca del núcleo sonda se torna impresindible en el cálculo del GCE. Ø En sistemas complejos como ser con vacancias, impurezas o superficies no debemos tratar de explicar un observable como el GCE a partir de MCP y sumas red, si en cambio, con estructura de bandas electrónicas, que pueden realizarse con un nivel altamente sofisticado no solo para sistemas bulk, sino también en superficies o sistemas con impurezas o vacancias.

Referencias [1] R. S. Raghavan, E. N. Kaumann and P. Raghavan, Phys. Rev. Let. 34 1280 (1975). [2] P. Blaha, K. Schwarz and P. H. Dederichs, Phys. Rev. B 37 2792 (1988). [3] E. N. Kaufmann and R. J. Vianden, Rev. Mod. Phys. 51 161 (1979).

Fín Obrigado