Transcript SOLUCIONES I NICAS BIDIMENSIONALES: UN ESTUDIO DE DIN MICA MOLECULAR

SOLUCIONES IÓNICAS BIDIMENSIONALES: UN ESTUDIO DE
DINÁMICA MOLECULAR
G. A. Méndez Maldonado1, M. González Melchor2, H. Ruiz Estrada1.
de Ciencias Físico Matemáticas, 2 Instituto de Física Luis Rivera Terrazas, Benemérita Universidad Autónoma
de Puebla.
E-mail: [email protected], [email protected]
1Facultad
Resumen
Presentamos un estudio numérico de una mezcla binaria bidimensional de discos cargados suspendidos en agua
en equilibrio termodinámico. Usamos dinámica molecular [1] y una modificación del modelo primitivo
restringido para obtener las funciones de distribución radial y la energía interna del sistema.
 Sistema: N iones inmersos en agua en equilibrio
termodinámico a temperatura T.
qi  Z i e
 Z
i
Interacción electrostática (método de Ewald 2D) [2]
N 1 N
U  
l
qi q j erfc( rij )
rij
i 1 j 1
Interacción de corto alcance
i
0
A
i

N
 
qi q j erfc(G / 2  )
qi2
  
exp(

i
G

r
)



ij
A i 1 j 1 G  0

G
i 1

N
N
N 1 N
 
U s   Aij  ij 
i 1 j i
 r 
n
 Resultados
Dinámica molecular
♦ Algoritmo de velocidades de Verlet
fi riU
♦ Simulaciones a número N, área A y
temperatura T constantes
♦ Cálculo de funciones de distribución radial
♦ Propiedades termodinámicas
Energía potencial electrostática
u* weis
[3]
u*
ρ*=0.3
ρ*=0.6
T*=0.1
6
T*=0.1
6
-2.9375
± 0.0029
-2.93
-3.5052
± 0.0041
-3.47
T=117.52 K
ρ*=6.6 x10 18 par/m2
ρ*=3.3 x10 18 par/m2
Distancia
media
d=3.89 A
d=5.5 A
g(r) para una densidad ρ*=0.3 a T=0.166 con u* =-3.5052 ± 0.004
g(r) para una densidad ρ*=0.6 a T=0.1 con u* =-6.1108 ±
Bibliografia
0.0038
[1] M.P. Allen and D.J. Tildesley, Computer Simulation of liquids, Clarendon Press Oxford, (1987).
[2] G. Gao, X. Zeng, Vapor-Liquid coexistence of quasi-two-dimensional Stockmayer fluids, J.Chem.Phys. 106, 3311 (1997).
[3] J.J.Weis, D.Levesque, Restricted primitive model of an ionic solution confined to a plane, J.Chem.Phys.109,7486 (1998).