Rolando Belardinelli, Sergio Manzi y Víctor Pereyra

1. Universidad Nacional de San Luis Facultad de Ciencias Físico-Matemáticas y Naturales Departamento de Física Instituto Nacional de Física Aplicada El algoritmo 1/t, un método eficiente para el calculo de la función de estado en la física de superficie Rolando Belardinelli, Sergio Manzi y Víctor Pereyra Simposio sobre Adsorción Adsorbentes y sus Aplicaciones (SAASA). San Luis 2009.

2. En el marco de la mecánica estadística toda magnitud en equilibrio puede ser obtenida a partir de la función de partición: • donde  es el estado del sistema, E es la energía correspondiente a , kB es la constante de Bolzman y T es la temperatura. La función de partición la podemos re-escribir agrupando los sistemas  que tienen la misma energía: • donde g(E) es el numero de estados con energía E.

3. Es una tarea dificultosa calcular directamente la g(E) con precisión en sistemas relativamente grandes. En el 2001, Wang y Landau presentaron un método de MC el cual calcula la g(E) en tamaños grandes, y con reducido esfuerzo computacional. Este método permite acceder directamente a la energía libre y entropía, es independiente de temperatura, y es eficaz en el estudio de transiciones de fase de primer y segundo orden. F. Wang and D. P. Landau, ‘‘Efficient, multiple-range random walk algorithm to calculate the density of states,’’ Phys. Rev. Lett. 86, 2050–2053

4. El algoritmo se ha usado en diferentes sistemas tales como, modelo de Ising, modelo de Potts, sistemas de espines aleatorios, sistemas cuánticos, fluidos, vidrios binarios de Lennard-Jones, cristales líquidos, polímeros, proteínas, sistemas moleculares, racimos atómicos, problemas de optimización, y teoría de numeración combinatoria... La g(E) calculada no converge a la verdadera densidad de estados. Algoritmo 1/t R. E. Belardinelli, and V. D. Pereyra.“Fast algorithm to calculate density of states”, Physical Review E 75, 046701 (2007). R. E. Belardinelli, and V. D. Pereyra.“Wang-Landau algorithm: a theoretical analysis of the saturation of the error”, Journal of Chemical Physics. 127, 184105 (2007).

5. Algoritmo 1/t 1) Se elige una configuración al azar del sistema; se inicializa S(E)=ln[g(E)]=0; se fija Fo=1 y también Ffinal= 10 -8. 2) Se sortea una partícula adsorbida y se la intenta mover a un lugar vacío de acuerdo a la probabilidad, 3) Se incrementa S(E)S(E)+Fk. 4) Luego de 1000 pasos de MC, se verifica que el programa haya visitado todas las energía “E” accesibles del sistema con el mismo Fk, si visito todos entonces se reduce Fk+1 Fk/2. 5) Si Fk+1 1/t, entonces Fk+1=F(t)=1/t. El paso 4) no se usa mas. 6) Si Fk<Ffinal el computo termina, sino se vuelve a 2).

6. Monómeros en Redes Triangulares Aplicamos el algoritmo 1/t al modelo de gas de red para la adsorción de partículas sobre una red triangular en el régimen de monocapa. Nos concentraremos en la interacción adsorbato-adsorbato (NN) repulsiva, ya que provoca dos fases bien definidas en la red triangular, una a cubrimiento 1/3 y la otra a 2/3,

7. Cada sitio de la red podrá ser ocupado por solo un monómero. Sea c1,c2,c3,....,cM la ocupación del sitio 1,2,3...,M respectivamente. Cada ci puede ser 0 o 1 de acuerdo si el sito esta vacío o ocupado por una partícula. El número de sitios ocupados es n, y el cubrimiento es definido como =n/L2. Suponemos solo interacción a primeros vecinos. Definimos el número de interacciones como, El Hamiltoniano del sistema puede ser escrito como, donde w es la energía de interacción a primeros vecinos, siendo positiva w>0 para interacciones repulsivas y negativas w<0 para atractivas

8. Estudiaremos la temperatura crítica de la fase 2/3 en la asamblea canónica, es decir a cubrimiento fijo =2/3.

10. La Tc es bastante consistente respecto a la Tc reportada por Pasinetti et al., Tc(∞)=0.3354  0.0001. P. M. Pasinetti, F. Romá, A. J. Ramirez-Pastor and J. L. Riccardo,” Critical behavior of repulsive linear k-mers on triangular lattices” cond-mat/0606391. Phys. Rev. B 74, 155418 (2006).

11. Dímeros en Redes Rectangulares Cada sitio de la red podrá ser ocupado por una unidad del dímero. Cada ci puede ser 0 o 1 de acuerdo si el sito esta vacío o ocupado por una partícula. El número de sitios ocupados es n, y el cubrimiento es definido como =n/L2. Suponemos solo interacción a primeros vecinos. Definimos el número de interacciones como, El Hamiltoniano del sistema puede ser escrito como, donde w es la energía de interacción a primeros vecinos, siendo positiva w>0 para interacciones repulsivas y negativas w<0 para atractivas

13. ADSORCIÓN DE DÍMEROS EN SUBSTRATOS BIDIMENSIONALES CON INTERACCIÓN ANISOTRÓPICAS La Energía de interacción es anisotrópica, es decir que depende de la forma direccional de contacto entre dímeros. Como indica la figura tenemos tres tipos de interacción, una Longitudinal JL, transversal JT y mixta JM =(JT+JL)/2. El Hamiltoniano del sistema es Y la gran función de partición,

14. La isoterma de adsorción la podemos calcular a partir de, Red 36x36 JT=5 JL=5

16. El valor de la Energía por sitio desde la gran función de partición: JT=5 JL=5

17. JL=5 JT=5

18. El parámetro de orden lo definimos como  = |(NV-NH)/N, donde NV es el numero de dímeros orientados verticalmente, NH horizontalmente y N es la cantidad total de dímeros sobre la red (N = NH+NV).

19. Conclusiones • Los métodos de simulación en donde la probabilidades de transición dependen de la temperatura son poco eficiente. • El calculo de la función densidad de estados es una de las herramientas alternativas para mejorar la eficiencia en el calculo de los observables en equilibrio pues es independiente de la temperatura. • Todo algoritmo donde el parámetro de refinamiento decrezca exponencialmente, es no convergente. En contraste el algoritmo 1/t es convergente. • El “Algoritmo 1/t” muestra ser el mas eficiente, simple y preciso en el cálculo de la densidad de estado tanto en sistemas discretos como aproximaciones al continuo.