I. INTRODUCCIÓN
Las nanoestructuras magnéticas unidimensionales, como nanohilos (NHs) y nanotubos (NTs), se caracterizan por tener una gran relación de aspecto y área superficial, por lo que estas estructuras a escala nanométrica adquieren propiedades magnéticas distintas a las que se observan a gran escala.1,2 Esto les otorga características únicas de interés para potenciales aplicaciones tecnológicas. En biomedicina, por ejemplo, estas nanoestructuras se utilizan como agentes de dispersión para la obtención de imágenes mediante resonancia magnética3,4 o también en suministro de fármacos,5 tecnología de la información para el desarrollo de memorias espintrónicas de estado sólido,6,7 entre otras. Las nanoestructuras magnéticas pueden estar constituidas por varios materiales. Entre las aleaciones más prometedoras para aplicaciones tecnológicas se encuentra el Permalloy (Py),8,9 un material ferromagnético blando, el cual está compuesto en un 80% por níquel y en un 20% por hierro.
Los NHs pueden ser caracterizados a través de su diámetro externo y su longitud (L), mientras que los NTs exhiben un parámetro adicional que es el espesor de la pared del tubo , tal y como puede observarse en la Fig. 1. Al aplicar un campo magnético externo en direcciones paralela y perpendicular al eje fácil (z) de las nanoestructuras (Fig. 1), es posible observar ciertos cambios en sus propiedades magnéticas, variando la magnetización de remanencia , la coercitividad y la magnetización de saturación . También, pueden estudiarse los distintos modos de reversión de la magnetización, la cual puede darse mediante rotación coherente de los momentos magnéticos o a través de la nucleación y propagación de paredes de dominio.10 Los más usuales para nanoestructuras cilíndricas son: coherente (C), transversal (T) y vórtice (V). Así mismo, recientemente se ha encontrado el modo wave (W) para nanotubos con ciertos parámetros geométricos y para la situación en la cual el campo magnético externo se aplica de manera perpendicular al eje z. Estos modos de reversión se observan en la Fig. 2.11,12 También se encuentra en la literatura el modo pseudo-coherente, en el cual casi todos los momentos magnéticos revierten de forma coherente.11 Estos modos de reversión dependen del ángulo de aplicación del campo magnético externo con respecto al eje fácil y, además, del tipo y variables geométricas de las nanoestructuras.
Consecuentemente, resulta interesante estudiar cómo influyen en las propiedades magnéticas cuasi-estáticas, tanto el tipo de nanoestructura como sus parámetros geométricos, con el fin de controlar sus comportamientos en potenciales aplicaciones. Es por ello que las simulaciones micromagnéticas constituyen una herramienta fundamental para describir, predecir y comprender el efecto de cada factor en las propiedades magnéticas.
En este trabajo, se presentan los resultados obtenidos de simulaciones micromagnéticas de un NH y de un NT de Permalloy de 10 m de longitud en función de sus diámetros externos, cuando se aplica un campo magnético externo en la dirección paralela y perpendicular al eje de fácil magnetización de estas nanoestructuras. Además, se discute en forma comparativa estos resultados con los obtenidos para nanoestructuras análogas de menor longitud, simuladas bajo las mismas condiciones.
II. SIMULACIONES MICROMAGNÉTICAS
La ecuación de Landau-Lifshitz-Gilbert (LLG) describe la evolución de la magnetización con respecto al tiempo13
Donde es la magnetización, es la magnetización de saturación, γ es la razón giromagnética, 𝛼 es el coeficiente de amortiguamiento y es el campo magnético efectivo que incluye el campo magnético externo, el de desmagnetización, el de anisotropía y el de intercambio.14
Las simulaciones micromagnéticas fueron realizadas utilizando el software libre Object Oriented Micromagnetic Framework (OOMMF),15 el cual resuelve la ecuación de LLG de forma iterativa para cada celda de una malla seleccionada.
Se estudiaron las propiedades magnéticas de dos configuraciones, un NT y un NH de Py, con diámetros externos nm y 60 nm y con m; en particular, para el NT se utilizó un nm. La constante de intercambio para el Permalloy fue J/m, la magnetización de saturación A/m y el coeficiente de amortiguamiento, , igual a 0,5.12,16 Las nanoestructuras se consideraron policristalinas, por lo que se despreció la anisotropía magnetocristalina.10 Todas las nanoestructuras se discretizaron tomando en consideración un tamaño de celda de nm3. Dado que se simularon NHs y NTs muy largos, se eligió un tamaño de celda suficientemente grande a lo largo del eje fácil (z) para completar las simulaciones en un tiempo razonable, pero lo suficientemente pequeño en el plano xy, con valores menores que su longitud de intercambio (5.3 nm) para reproducir la geometría cilíndrica de estas nanoestructuras. Este criterio de discretización ha sido elaborado para reducir la rugosidad de los bordes provocada por la simulación de un cilindro con celdas cúbicas.
Con el fin de observar la variación de las propiedades magnéticas de las nanoestructuras estudiadas en función del campo magnético externo, H, se aplicó éste en la dirección paralela y perpendicular al eje de fácil magnetización de estas nanoestructuras, es decir, con y respectivamente, como se observa en la Fig. .
III. RESULTADOS Y ANÁLISIS
En las Figs. 3 y 4 se visualizan las diferentes curvas de histéresis para los NHs y NTs simulados, respectivamente, en donde la magnetización se encuentra normalizada a la magnetización de saturación.
Se observa una marcada cuadratura para , lo que indica que esta es la dirección de fácil magnetización. Esto significa que, si el campo magnético aplicado es paralelo al eje fácil, éste debe ser lo suficientemente intenso para revertir la magnetización de la nanoestructura. Por otro lado, a , cuando el campo magnético externo es perpendicular al eje de fácil magnetización, se observa un comportamiento con coercitividad nula, permitiendo visualizar un cambio significativo de las propiedades magnéticas con respecto al campo magnético externo. Por ello, un análisis comparativo de las coercitividades de ambos sistemas (Tablas 1 y 2) indica que es necesario una mayor intensidad de campo magnético externo a que a para llevar la imanación a cero. Este comportamiento se debe a la gran anisotropía de forma de estas nanoestructuras, independientemente de si es un NH o un NT. De la misma manera, la remanencia de ambas estructuras es mayor para que la observada para .17
Por otro lado, se observa que para un diámetro mayor de NH (o NT), el campo coercitivo disminuye (ver Tablas 1 y 2). Esto se debe a que la interacción dipolar interna se vuelve más importante, produciendo dominios de cierre.10,18 Tal efecto facilita la nucleación y propagación de las paredes de dominio, dando lugar a una menor coercitividad. Un análisis comparativo entre ambas nanoestructuras evidencia que la coercitividad es mayor en un NH que en un NT. Este resultado se interpreta en función de que es más sencillo nuclear una pared de vórtice en los extremos de un NT que en los de un NH, por lo que se necesita un campo magnético menor para nuclear y propagar las paredes de dominio, responsables del proceso de reversión de la magnetización en estas nanoestructuras que tienen una relación de aspecto importante.
Asimismo, se han señalado en las Figs. 3 y 4 algunos puntos correspondientes a instantáneas de la magnetización, para los cuales pueden observarse los diferentes modos de reversión para ambas estructuras y diámetros estudiados, en función de la dirección de campo aplicado con respecto al eje fácil (Fig. 5 y 6). Para , la reversión de la magnetización se presenta como la propagación de paredes de dominio tipo vórtice, mientras que para se observa una reversión de la magnetización en modo pseudo-coherente. Además, al comparar estos valores para los mismos diámetros en función de la longitud, observamos que los valores de y prácticamente no varían de manera significativa12,17 (ver Tablas 1 y 2). Sin embargo, un resultado notable se evidencia para el NT con nm y . En efecto, se observa un cambio en el modo de reversión, respecto a lo observado para la misma nanoestructura con m.16 Así, en estas condiciones, se obtiene que para una longitud más corta (1 m), el NT revierte su magnetización mediante el modo wave, recientemente propuesto.12,16 Por ello, es importante notar cómo puede cambiar la forma en la que los momentos magnéticos de un material revierten su magnetización de forma distinta influenciada por un factor geométrico como la longitud.
IV. CONCLUSIONES
Se analizaron en forma comparativa las curvas de histéresis y modos de reversión para un NH y un NT de 10 de longitud, obtenidos mediante simulaciones micromagnéticas. Se observó que para , el proceso de reversión de la magnetización ocurre mediante la nucleación y propagación de paredes de dominio tipo vórtice, mientras que para se identificó el eje de difícil magnetización, con modo de reversión de la magnetización pseudo-coherente. Además, se observó que a un mismo diámetro ( nm), pero bajo dos longitudes distintas, la reversión de la magnetización es notablemente diferente cuando el campo externo se aplica de manera perpendicular al eje de fácil magnetización.
Los resultados obtenidos nos han permitido describir el efecto de los parámetros geométricos en las propiedades magnéticas de nanoestructuras unidimensionales de Py, lo que evidencia el rol fundamental de tales factores, así como el aporte de las simulaciones al entendimiento del comportamiento de tales nanoestructuras, de interés para el diseño, predicción y control de su conducta magnética en diversas aplicaciones nanotecnológicas.
AGRADECIMIENTOS
Los autores agradecen al Centro de Computación de Alto Desempeño (CCAD) de la UNC, por el acceso al cluster Mendieta.
J.E. agradece a Fondecyt 1200302 y al proyecto Basal AFB180001.
CRÉDITOS DE AUTOR
Selene Levis y Angelo Villegas: Investigación del problema. Procesamiento e interpretación de datos. Escritura del primer borrador.
Juan Escrig: Investigación del problema. Interpretación de datos.
Noelia Bajales Luna: Investigación del problema. Diseño y conceptualización. Financiamiento.
Diana Arciniegas Jaimes. Investigación del problema. Diseño y conceptualización. Orientadora del proyecto.