Simulación de Difusión Mediante El Método de Montecarlo Cinético

1 Laboratorio Ing.
Superficies e Intercaras 2019 – 2020
Simulación de difusión mediante el Método de Montecarlo Cinético (KMC)

▪ Introducción teórica migración, que es la energía necesaria para mover una
La difusión o movilidad atómica se hace a través de un vacante en un cristal dado; E* es una propiedad de cada
determinado mecanismo atómico que requiere, para su cristal y se relaciona principalmente con la intensidad del
desarrollo, temperatura y tiempo. La velocidad de la enlace presente.
difusión atómica, que es la velocidad a la cual el sistema
se podrá acercar a la situación de equilibrio, se relaciona
con la cinética química. Si la temperatura fuese
suficientemente baja, la cinética será prácticamente nula
y no se podrá llegar al equilibrio.
En cristales geométricamente densos, como lo son
aquellos de los metales puros y de sus soluciones sólidas
de sustitución, el mecanismo más importante de difusión
atómica se asocia a las vacantes.
El aumento en la temperatura T de un cristal tiene dos

efectos que conjugados permiten que los átomos se
muevan más rápido aprovechando las vacantes. El
primero está relacionado con el equilibrio químico y el
segundo con la cinética química.
El primer efecto depende del Equilibrio Químico (EQ) y,
Ilustración 1. Energía que debe vencer un átomo para saltar a una
por ende, NO depende del tiempo. Este efecto se traduce vacante vecina.
en que como a mayor temperatura T [K] hay más energía
cinética promedio por temperatura (EC= kBT, donde k es ▪ Técnicas y resultados experimentales
la constante de Boltzmann) en el cristal, también podrá Para realizar la simulación se ha utilizado el llamado
haber más desorden. En particular, a T= 0K un cristal al EQ Método de Montecarlo Cinético (KMC).
debiese estar libre de defectos cristalinos y, por ende, En el método de simulación Monte Carlo se intenta seguir
teóricamente, el cristal debiese ser perfecto. la ‘dependencia del tiempo’ de un modelo para el cual su
cambio o evolución no se desarrolla de una manera
Sin embargo, la energía necesaria para formar una rigurosamente predefinida, sino más bien de una manera
vacante en un cristal dado, Ev, es relativamente estocástica y que depende de una secuencia de números
moderada. Tal energía Ev es una propiedad del material y aleatorios generada durante la simulación. Con una
depende de la intensidad del enlace presente. Así, a una segunda secuencia diferente de números aleatorios, la
temperatura muy alta, próxima a la temperatura de simulación no arroja resultados idénticos, pero sí valores
fusión de cristal, por cada 1000 átomos podrá haber, en que están acordes con aquellos obtenidos de la primera
forma aproximadamente, 1 vacante al equilibrio por cada secuencia, generándose eso si un pequeño ‘error
mil átomos, Cv*≈ 1/1000 [número de vacancias / número estadístico’ pero conservándose las tendencias en los
de átomos]. resultados. En el método se habla también de los pasos
de Monte Carlo, pero existe una buena libertad para
El segundo efecto se relaciona con la Cinética, lo cual diseñar un paso de Monte Carlo. Este puede contener
implica considerar el tiempo. Esta cinética se refiere a la usualmente, N pasos elementales donde N es el número
velocidad promedio con que se mueven los átomos, Ra, y de partículas en el sistema; o bien, puede ser tomando
las vacancias Rv. Para estas velocidades las unidades son como ejemplo las simulaciones de la física del estado
saltos atómicos por unidad de tiempo (1s). Así, por sólido, un cambio en los primeros vecinos de un átomo o
ejemplo, si una vacante cambia de sitio atómico 10 veces vacancia en una red cristalina. Dentro de las simulaciones
en un segundo, entonces: Rv= 10 [saltos atómicos /s]. con Monte Carlo se encuentra una buena parte dedicada
Para que una vacante se mueva un sitio atómico vecino, al crecimiento de películas o recubrimientos, y se pueden
se requiere que el átomo de ese sitio salte a la vacante, distinguir algunos métodos como deposición balística que
ver Ilustración 1. Para ello tal átomo debe superar una en términos computacionales es la más sencilla de
barrera de potencial E*, denominada energía de implementar y que los resultados son fácilmente
1
2 Laboratorio Ing. Superficies e Intercaras 2019 – 2020
interpretados, también se habla del método de

sedimentación que es básicamente una versión mejorada
del método de deposición balística pero añadiendo el
proceso de difusión de cada átomo que llega al sustrato
con la inclusión también de reglas para el salto en la
difusión.
Ilustración 2. Esquema comparativo entre Monte Carlo (MC) y Monte

Carlo Cinético (KMC).
Recientemente la atención se ha centrado en el estudio

Ilustración 3. Evolución del sistema para KT=0.1 e intervalos de 1.
de una variante del método Monte Carlo en equilibrio
llamada Monte Carlo Cinético (KMC) y que se utiliza para Para KBT=0.3 e intervalos semejantes el resultado es muy
simular la técnica de crecimiento llamada “Molecular similar, aunque en este caso podemos ver que se han
Beam Epitaxy” o crecimiento epitaxial por haces formado un número menor de islas que en la situación
moleculares. El Monte Carlo cinético puede ser aplicado anterior dado que el aumento de temperatura y, por
a sistemas en condiciones de equilibrio y fuera del tanto, el aumento de energía hace que el desorden sea
equilibrio y para crecimiento epitaxial por haces algo mayor.
moleculares trata tanto los procesos de deposición como
La evolución del sistema la podemos observar en la
los de difusión; la principal diferencia con el Monte Carlo
tradicional aplicado a la difusión de átomos en la Ilustración 4.
superficie es que debe calcularse la energía de activación
del sistema y tenerla en cuenta al realizar un paso de
Monte Carlo entre dos sitios con energía E1 y E2 y aceptar
una configuración nueva del sistema (E2) solo sí existe la
energía de activación necesaria EA. Esta relación se
observa en la Ilustración 2.
Resultados experimentales
Para bajas temperaturas (KBT=0.1) e intervalos bajos
entorno a τ=1, podemos observar que se van formando
islas de vacantes por diferentes capas. Esto se debe a que
la baja energía del sistema hace que cuando se produce
la coalescencia de vacantes, no haya suficiente energía
para que estas sigan difundiendo por la superficie, es
decir, la energía cinética del sistema es muy baja.
Sin embargo, si aumentamos la energía del sistema a

KBT=0.9 para un mismo intervalo τ=1, la alta energía del
sistema se traduce en un aumento de la difusión de las
2
3 Laboratorio Ing. Superficies e Intercaras 2019 – 2020
vacantes haciendo que la inestabilidad del sistema ▪ Referencias

aumente dando lugar a un gas de vacantes. El - Desarrollo de un Código Kinetic Monte Carlo para
movimiento de las vacantes es constante por la su Aplicación en Ciencia de Materiales. TFM.
superficie. Solo encontramos vacantes en la primera Universidad de Cantabria. Pablo Martín García.
superficie y no se aprecia formación de islas dado que la
cinética es alta. - Simulación Monte Carlo cinético de la difusión
atómica en la aleación FeAl. Universidad Nacional
Esta evolución del sistema la podemos comprobar en la
Mayor de San Marcos. Tesis Erich Víctor Manrique
siguiente Ilustración 5.
Castillo.
- Fundamentos de simulación de materiales por

medio del método de Monte Carlo. Giovany Orozco
Hernández1, Jhon Jairo Olaya Flórez2, Elisabeth
Restrepo Parra3.
▪ Conclusiones
Por tanto, como conclusión podemos decir que un
aumento de la temperatura provoca un aumento de la
energía cinética del sistema, lo que se traduce en una
aceleración de la difusión de las vacantes que tienen
energía suficiente para moverse con facilidad. De este
modo, altas energías no favorecerán la formación de islas
de vacantes ni la erosión capa a capa ya que la
inestabilidad es alta.
Por el contrario, bajas temperaturas, conllevan energías
inferiores que resultan en la formación de islas de
vacantes y en la erosión capa a capa ya que de este modo
el sistema se encuentra en una mayor estabilidad.
Menores temperaturas también favorecen que el Autores:
aumento de las islas de vacantes sea mayor; a medida Ye Yuan
que elevamos la temperatura este tamaño de islas se irá Jesús García-Arisco Rivera
reduciendo por el aumento de la inestabilidad del sistema
al crecer la movilidad de las vacantes.

Simulación de Difusión Mediante El Método de Montecarlo Cinético

Cargado por

Copyright:

Formatos disponibles

Simulación de Difusión Mediante El Método de Montecarlo Cinético

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Simulación de Difusión Mediante El Método de Montecarlo Cinético

Cargado por

Copyright:

Formatos disponibles

1 Laboratorio Ing.

Superficies e Intercaras 2019 – 2020

Simulación de difusión mediante el Método de Montecarlo Cinético (KMC)

El aumento en la temperatura T de un cristal tiene dos

interpretados, también se habla del método de

Ilustración 2. Esquema comparativo entre Monte Carlo (MC) y Monte

Recientemente la atención se ha centrado en el estudio

Ilustración 4. Evolución del sistema para KT=0.3 e intervalos de 1.

Sin embargo, si aumentamos la energía del sistema a

vacantes haciendo que la inestabilidad del sistema ▪ Referencias

- Fundamentos de simulación de materiales por

Ilustración 5. Evolución del sistema para KT=0.9 e intervalos de 1.

También podría gustarte