Spanish (pdf)
Article in xml format
Article references
Automatic translation
Send this article by e-mail
Cited by SciELO 
Cited by Google 
Similars in
SciELO 
Similars in Google 
Permalink
Introducción
Compuestos binarios III-arsénico exhiben propiedades que los hacen atractivos para aplicaciones comerciales y tecnológicas, entre ellos se encuentra el GaAs, usado en diferentes dispositivos electrónicos y optoelectrónicos tales como diodos emisores de luz, foto detectores, láseres y circuitos integrados, entre otros, debido a que posee gab directo, (Ahmed, et al. , 2007).
Los semiconductores que no poseen propiedades magnéticas, se pueden dopar con elementos magnéticos, variando así las propiedades magnéticas del semiconductor, éstos se conocen como semiconductores magnéticos diluidos (DMS Diluted Magnetic Semiconductors) (He & HaiMing, 2008), de amplia aplicación en espintrónica. Diversos estudios realizados con el GaAs confirman que al adicionar impurezas de algunos metales de transición 3d, este compuesto exhibe un momento magnético considerable.
Para analizar las modificaciones de las propiedades electrónicas, estructurales y magnéticas del Ga1-x Cr x As, se varió la concentración del Cr, se trabajó con celdas 1x1x1, 1x1x2 y 2x2x2, que corresponden a concentraciones de Cr de 0.25, 0.125 y 0.03125 respectivamente. Concentraciones que están de acuerdo con lo reportado por Shirai (2001), Baykov, et al. (2007) y Sandratskii & Bruno (2003).
Optimización de parámetros
Los cálculos se realizaron con el software quantum-espresso (Giannozzi, et al. , 2009), se trabajó con el formalismo de la teoría del funcional de la densidad DFT, se empleó el método de ondas planas y pseudopotenciales atómicos ultra-suaves. Se trabajó con la aproximación de Gradiente generalizado (GGA) en la parametrización de Perdew, Burke y Ernzerhof (PBE), J. P. Perdew, et al. , 1996. El GaAs se estudió en la estructura zincblenda, con una constante de red optimizada de 5.66 Å (10.85 Bohr), se empleó una malla de puntos k de 6x6x6.
Se realizó optimización de la energía de corte Ecutoff y la malla de puntos k con respecto al volumen de equilibrio y a la constante de red, como se muestra en el Tabla 1.
Resultados y discusiones
Teniendo en cuenta los resultados que se presentan en el cuadro 1, se grafica EvsV para cada una de las energías de corte (Figura 1). El valor de la constante de red que hace mínima la energía se encuentra a partir de una malla de puntos k de 6x6x6, para valores de E cutoff de 30 Ry, 35 Ry y 40 Ry, de igual manera para las energías de corte mencionadas, con una malla de puntos k de 10x10x10 el valor de la constante de red es similar al encontrado con la malla de 6x6x6. Teniendo en cuenta lo mencionado y el coste computacional al correr los cálculos con los parámetros descritos, se selecciona la E cutoff de 30 Ry con una malla de puntos k de 6x6x6, para el valor de a que minimiza la energía correspondiente al ao(10.85 bohr).
Figura 1 Energía Vs Volumen para el GaAs en su estructura zincblenda. Para una energía de corte de (a) 25 Ry, (b) 30 Ry, (c) 35 Ry y (d) 40 Ry.
En el Tabla 2, se compara el valor de la constante de red y el módulo de volumen calculado, para el GaAs en la estructura zincblenda, con otros datos reportados en la literatura. Se relacionan los porcentajes de error para cada uno de los parámetros escogidos.
Los valores reportados por Ahmed, et al. , (2007) son de carácter teórico, calculados mediante el método FP-LDA y FP-GGA por sus siglas en inglés Full-potential linear augmented plane wave (FP-LAPW) con la aproximación de densidad local (LDA - Local density aproximation) y aproximación de gradiente generalizado (GGA -generalized gradiente approximation), respectivamente, cálculos realizados mediante código WIEN2K; los datos reportados por Hellwege, et al. , (1982) y Fillipi, et al. , (1994) son de carácter experimental.
La diferencia en los valores obtenidos se debe a los métodos empleados para hallar dicho parámetro.
Propiedades electrónicas y estructurales
Las gráficas de bandas de energía que se presentan, se hicieron para el camino Γ, M, X, Γ, R, X, M, R; camino de alta simetría de la Primera Zona de Brillouin, lo enunciado se trabajó para la celda zincblenda y la supercelda 2x2x2. Para la supercelda 1x1x2, se tuvo en cuenta el camino de alta simetría Γ, M, X, Γ, A, M, X, A, en la PZB.
A continuación se presentan los resultados para el GaAs en su estructura zincblenda.
De acuerdo con la Figura 2 el GaAs posee gap directo en Γ, con un valor de 0.25 eV, dicho valor permite ver que el compuesto es un semiconductor; éste resultado concuerda con otros datos teóricos reportados por Ahmed, et al., (2007); Grivickas, et al., (2009); Ihm & Joannopoulus (1981); Sze & Ng (2007); Wang & Klein (1981). Resultados similares se encuentran para las superceldas 1x1x2 y 2x2x2.
Al sustituir un átomo de Ga por uno de Cr en la red, se encuentra que las bandas de energía difieren de la estructura pura. Se encuentra que para una orientación de espín presenta un comportamiento conductor, y para la otra, un comportamiento aislante R.A.Groot, et al. (1983); Sandratskii & Bruno (2003), (Figura 3), es decir el sistema exhibe un comportamiento half-metallic.
Figura 3 Bandas de energía para el GaAs 1x1x1 spin up (a) y spin down (b), al sustituir un átomo de Ga por uno de Cr.
En la Figura 3(a) se observa un traslapamiento de bandas de valencia y conducción, hecho que evidencia un carácter metálico del Ga1-x Cr x As (x=0.25). Cerca al nivel de Fermi y hasta -4 eV, la contribución de orbital s corresponde a orbitales 4p del Ga, 4p del As y 3d del Cr.
En la Figura 3(b) se observan tres regiones por debajo del nivel de Fermi. La región más cerca al nivel de Fermi, entre -0.7 eV y -4 eV, tiene la mayor contribución en electrones 4p del As y 4p del Ga, la contribución de electrones 3d no es mayoritaria como en el caso de las bandas con spin up, se observa la presencia de un gap de energía.
Al analizar las gráficas de bandas de energía para el Ga Cr As (graficas 3 (a) y 4 (a)), sustituyendo un átomo 1-x x de Ga por uno de Cr en las estructuras 1x1x1 y 1x1x2 respectivamente, se confirma el comportamiento “half-metallic” del sistema.
Figura 4 Bandas de energía para el GaAs 1x1x2 spin up (a) y spin down (b), reemplazando un átomo de Ga por uno de Cr.
En la Figura 3(b) se observan tres regiones por debajo del nivel de Fermi. La región más cerca al nivel de Fermi, entre -0.7 eV y -4 eV, tiene la mayor contribución en electrones 4p del As y 4p del Ga, la contribución de electrones 3d no es mayoritaria como en el caso de las bandas con spin up, se observa la presencia de un gap de energía.
Al analizar las gráficas de bandas de energía para el Ga1-x Cr x As ((gráficas 3 (a) y 4(a)), sustituyendo un átomo de Ga por uno de Cr en las estructuras 1x1x1 y 1x1x2 respectivamente, se confirma el comportamiento “half- metallic” del sistema.
Al hacer la sustitución de un átomo de Ga por uno de Cr en la supercelda 2x2x2, la concentración del dopante es de 0.03125. Al comparar las Figuras 3, 4 y 5, se encuentra que la contribución de orbitales es similar para las estructuras 1x1x1, 1x1x2 y 2x2x2 como se ha descrito en párrafos anteriores, al sustituir un átomo de galio por uno de cromo, el compuesto presenta magnetización.
Figura 5 Bandas de energía para el GaAs 2x2x2 spin up (a), spin down, al hacer una sustitución de un átomo de Ga por uno de Cr.
En el Tabla 3 se presenta la variación de la constante de red en función de la concentración, para las estructuras estudiadas.
Tabla 3 Constante de red para el GaAs con una concentración de Cr de 0.25, 0.125 y 0.03125, correspondientes a la estructura zincblenda, y a las superceldas 1x1x2 y 2x2x2 respectivamente.
Al hacer la sustitución de un átomo de Ga por uno de Cr, se encuentra que la constante de red no varía significativamente (analizando el caso en cada estructura), esto se debe a que la distribución de carga es similar para el átomo de Ga y el de Cr.
En el Tabla 4 se presenta el momento magnético en Ga Cr As, para las tres concentraciones trabajadas x=0.25, 1-x x x=0.125 y x=0.03125. Se muestra el momento magnético de la impureza 3d, quien induce momento en los átomos vecinos de As. Igualmente se presenta el momento magnético total de la celda.
El momento magnético total en los casos presentados al hacer la sustitución de un átomo de Ga por uno de Cr en las estructuras consideradas es de 3 µB, independiente de la concentración del dopante, esto está de acuerdo con lo reportado por Sandratskii (Sandratskii & Bruno, 2003), y por De-ming Ma (De-min Ma & Shi, 2016); la magnetización se presenta por la hibridación de los orbitales As-4p y Cr-3d.
Conclusiones
El GaAs en su estructura zincblenda, posee gap directo. Este hecho lo convierte en un material atractivo para su uso en celdas solares, heteroestructuras para emisores de luz, entre otros. Cerca al nivel de Fermi para este material se encuentra que la contribución electrónica corresponde a electrones 4p del As y 4p del Ga, y estando aun en el core, entre -4 eV y -7 eV la contribución corresponde a electrones 4p del As y 4s del Ga. Dicho material no presenta magnetización. Al hacer la sustitución de un átomo de Ga por uno de Cr para diferentes concentraciones de dopante, el material presenta un comportamiento “half-metallic”, dicho material puede ser usado como un semiconductor magnético diluido, como inyector de espín, también es un material que se puede usar en aplicaciones en espintrónica.
La magnetización se presenta por la hibridacion de orbitales As-4p y Cr-3d.
