Adsorción de H y O sobre Cu(100): DFT, simulaciones Monte Carlo y percolación.

ELIZABETH DEL VALLE GÓMEZ; LUCÍA B. AVALLE; M. CECILIA GIMÉNEZ

La Plata, Buenos Aires, Argentina.

Congreso; 102a Reunión de la Asociación Física Argentina; 2017

UNLP.

p { margin-bottom: 0.25cm; line-height: 120%; }Elestudio del hidrógeno y oxígeno adsorbido sobre superficiesmetálicas representa un tema atractivo para la ciencia de lassuperficies desde hace tiempo, por lo que viene siendo investigado através de diferentes métodos.Sepresenta un trabajobasado en el empleode simulacionesMonte Carlo en el ensamble canónico ycálculos de primeros principios, paraobtener unconocimiento mas detallado de laadsorción de átomos de hidrógeno y oxígeno sobre Cu(100).Los cálculos DFT se llevan a cabo utilizando el paqueteQuantum-ESPRESSO. Se modela una celda p(3x3) con 5 láminas del metalpara describir la adsorción de las especies atómicas (H u O) condiferente cantidad de primeros vecinos (1, 2, 3 y 4). Así obtenemosla energía de adsorción de un átomo de H u O en función de lacantidad de primeros vecinos de su misma especie.Las simulaciones Monte Carlo del presente estudio se complementan conlos resultados obtenidos a partir de los cálculos DFT, donde seconsideran los diferentes entornos posibles, dependiendo de lapresencia de otros átomos de H u O en los sitios de adsorciónvecinos. Finalmente, se estudia la percolación de H y de O sobre lasuperficie de Cu(100). Para ello de realizan simulaciones de MonteCarlo en el ensamble canónico, de acuerdo a la dinámica de Kawasakiy, para cada grado de cubrimiento, se calcula la probabilidad de queel sistema percole, R(\theta), para luego graficar R en función de\theta. La probabilidad de percolación, R, se calcula contando elnúmero de configuraciones en las cuales existe un cluster depercolación y dividiendo este número por el número deconfiguraciones totales consideradas. Estas curvas se realizan paradiferentes tamaños de sistema, con el objeto de extrapolar el umbralde percolación para un sistema de tamaño infinito, obteniendo así\theta_P para las condiciones dadas. Este estudio se realiza paradiferentes temperaturas, T, con el objeto de obtener un diagrama defase de \theta_P en función de T, es decir una curva que separe lasregiones percolantes de las no percolantes en el plano \theta ? T.