VIRIAL
VIRIAL
VIRIAL
u(r)
r 0
Potencial de Pozo Rectangular
es el modelo ms sencillo que tiene en cuenta a las fuerzas
atractivas y repulsivas en forma simultnea
consiste en una esfera repulsiva de dimetro rodeada de un foso
atractivo de profundidad que se extiende desde hasta g
este potencial ha sido ampliamente usado para el clculo del
segundo coeficiente virial
es el modelo ms sencillo que tiene en cuenta a las fuerzas
atractivas y repulsivas en forma simultnea
consiste en una esfera repulsiva de dimetro rodeada de un foso
atractivo de profundidad que se extiende desde hasta g
este potencial ha sido ampliamente usado para el clculo del
segundo coeficiente virial
Expresin matemtica:
u(r) = r <
u(r) = - < r < g
u(r) = 0 r > g
Forma:
u(r)
u(r)
0 r
Potencial de Pozo Triangular
Expresin matemtica:
u(r) = r <
< r < g
u(r) = 0 r > g
Forma:
u r
g
g
r
( ) =
|
]
|
c
o 1
Expresin matemtica:
u(r) = r <
< r < g
u(r) = 0 r > g
Forma: u(r)
r
g
Modelos con Campo Central de Fuerzas
Algunos de estos modelos son:
Lennard-Jones Rydberg
Potencial e
-6
Singer
Potencial e
-6
, e
-8
Guggenheim-McGlashan
Carr-Konowalow Morse
Potencial 12-6-4 Potencial generalizado
Vamos a tratar: Lennard-Jones
Algunos de estos modelos son:
Lennard-Jones Rydberg
Potencial e
-6
Singer
Potencial e
-6
, e
-8
Guggenheim-McGlashan
Carr-Konowalow Morse
Potencial 12-6-4 Potencial generalizado
Vamos a tratar: Lennard-Jones
Potencial de Lennard-Jones
es uno de los potenciales ms representativos de este tipo y
uno de los ms utilizados
est constituido por un trmino de largo alcance (r
-m
) y otro de
corto alcance (r
-n
)
generalmente m = 12 y n = 6
es fcil de integrar para obtener el segundo coeficiente virial
(da valores muy cercanos a la realidad)
es uno de los potenciales ms representativos de este tipo y
uno de los ms utilizados
est constituido por un trmino de largo alcance (r
-m
) y otro de
corto alcance (r
-n
)
generalmente m = 12 y n = 6
es fcil de integrar para obtener el segundo coeficiente virial
(da valores muy cercanos a la realidad)
Expresin matemtica:
Forma:
u(r) 4
r r
12 6
=
'
|
.
|
'
|
.
|
|
]
|
u(r)
r
r
0
|
.
|
'
|
.
|
|
]
|
|
+
) )
4c
o o
b 3
u r
r r
u
con a
a b
( ) ( ; ) =
'
|
.
|
'
|
.
|
|
]
|
|
+
) )
4c
o o
b 3
donde los son los momentos dipolares y u(,) es la
contribucin debida al dipolo-dipolo
Modelos de Fuerzas Pseudocentrales
Estos potenciales eliminan la dificultad de uso de los potenciales
con fuerzas que dependen de la orientacin
Entre este tipo de modelos podemos encontrar
potencial de cuatro centros
potenciales cuasiesfricos
potencial de Kihara
Vamos a tratar: potencial de Kihara
Estos potenciales eliminan la dificultad de uso de los potenciales
con fuerzas que dependen de la orientacin
Entre este tipo de modelos podemos encontrar
potencial de cuatro centros
potenciales cuasiesfricos
potencial de Kihara
Vamos a tratar: potencial de Kihara
Potencial de Kihara
puede aplicarse a molculas no esfricas (ms realista que Lennard-
Jones)
este modelo reemplaza a la molcula por un cuerpo rgido convexo
rodeado de un campo de fuerzas
supone que la energa de interaccin entre dos ncleos moleculares
slo depende de la distancia ms corta entre las superficies de los
dos ncleos moleculares
su expresin analtica es igual a la de Lennard-Jones:
no es la distancia entre los centros de las molculas, sino
la mnima distancia a la que pueden acercarse dos cuerpos rgidos
o ncleos moleculares.
puede aplicarse a molculas no esfricas (ms realista que Lennard-
Jones)
este modelo reemplaza a la molcula por un cuerpo rgido convexo
rodeado de un campo de fuerzas
supone que la energa de interaccin entre dos ncleos moleculares
slo depende de la distancia ms corta entre las superficies de los
dos ncleos moleculares
su expresin analtica es igual a la de Lennard-Jones:
no es la distancia entre los centros de las molculas, sino
la mnima distancia a la que pueden acercarse dos cuerpos rgidos
o ncleos moleculares.
u( ) p c
o
p
o
p
=
'
|
.
|
'
|
.
|
|
]
|
|
4
12 6
Ecuaciones de Estado
Ecuacin de estado ms sencilla: la del gas ideal
P.v = R.T
buenos resultados para todos los gases a bajas presiones
desviaciones a altas presiones
Este ltimo inconveniente se puede resolver con varias ecuaciones de
estado. Una de stas es la del desarrollo virial, que tiene una base
terica dada por la termodinmica estadstica
Ecuacin de estado ms sencilla: la del gas ideal
P.v = R.T
buenos resultados para todos los gases a bajas presiones
desviaciones a altas presiones
Este ltimo inconveniente se puede resolver con varias ecuaciones de
estado. Una de stas es la del desarrollo virial, que tiene una base
terica dada por la termodinmica estadstica
P v
R T
1
B
v
C
v
D
v
2 3
= + + + +
Clculo del Segundo Coeficiente Virial
El segundo coeficiente virial es la magnitud ms usada para
calcular parmetros intermoleculares
Para un gas real, la termodinmica estadstica da la siguiente
expresin:
donde u(r) es el potencial de interaccin y k la constante de
Boltzmann (k = 1,380658.10
-23
J/K)
Vamos a usar las siguientes expresiones de u(r): esferas
rgidas, pozo cuadrado, pozo triangular, Lennard-Jonnes
(con y sin orientacin), y Kihara
B 2 e 1 r dr
2
u(r)
kT
0
2
=
|
]
|
'
|
.
|
}
El segundo coeficiente virial es la magnitud ms usada para
calcular parmetros intermoleculares
Para un gas real, la termodinmica estadstica da la siguiente
expresin:
donde u(r) es el potencial de interaccin y k la constante de
Boltzmann (k = 1,380658.10
-23
J/K)
Vamos a usar las siguientes expresiones de u(r): esferas
rgidas, pozo cuadrado, pozo triangular, Lennard-Jonnes
(con y sin orientacin), y Kihara
B 2 e 1 r dr
2
u(r)
kT
0
2
=
|
]
|
'
|
.
|
}
Modelo de Esferas Rgidas
B
2
3
2
3
=
t
o
es el dimetro de la esfera
Este potencial se usa como referencia de otros potenciales a
travs del coeficiente virial reducido
B
B
B
ER
2
2
=
Modelo de Pozo Cuadrado
B B e
ER
kT
2
3
1 1 1 =
'
|
.
|
|
|
]
|
|
( )
c
B
B
B
ER
2
2
=
es el rango del pozo atractivo (vara entre 1,5 y 2)
Coeficiente B
2
(T,x) reducido
) ( 1 exp
) , (
) , (
3
2
) ( 1 exp
3
2
3 *
0
*
0
2 * *
2
3
0
3 *
0
*
3
2
*
*
r d
T b
x T B
x T B
b si
r d
T
B
kT
T
r
r
}
}
|
]
|
|
.
|
'
= =
=
|
]
|
|
.
|
'
=
=
=
) ( 1 exp
) , (
) , (
3
2
) ( 1 exp
3
2
3 *
0
*
0
2 * *
2
3
0
3 *
0
*
3
2
*
*
r d
T b
x T B
x T B
b si
r d
T
B
kT
T
r
r
}
}
|
]
|
|
.
|
'
= =
=
|
]
|
|
.
|
'
=
=
=
|
]
|
t t o
o
o
, ) , )
A
g
kT g
=
c c
o 1
y B=
kT g - 1
donde
Modelo de Kihara
B 2 e 1 r dr
2
4
kT
12 6
0
2
=
|
]
|
|
|
|
'
|
.
|
'
|
.
|
'
|
.
|
|
]
|
|
}
t
c o
p
o
p
resolucin numrica
Modelo de Lennard-Jones
B 2 e 1 r dr
2
4
kT r r
0
2
12 6
=
|
]
|
|
'
|
.
|
'
|
.
|
'
|
.
|
|
]
|
}
resolucin numrica
Modelo de Lennard-Jones con orientacin Modelo de Lennard-Jones con orientacin
B 2 e 1 r dr
con a b 3
2
4
kT r r
0
2
a b
=
|
]
|
|
) )
'
|
.
|
'
|
.
|
'
|
.
|
|
]
|+
}
u( ; )
resolucin numrica
Clculo de las Propiedades Termodinmicas Residuales
Relaciones fundamentales:
dU = T dS - P dV (1)
dH = T dS + V dP (2)
dA = -P dV - S dT (3)
dG = V dP - S dT (4)
G es funcin de variables medibles: temperatura y presin
Aplicando el diferencial total a G/(RT):
(5)
Reemplazando a dG por (4) y a G por G = H - T S se llega a:
Relaciones fundamentales:
dU = T dS - P dV (1)
dH = T dS + V dP (2)
dA = -P dV - S dT (3)
dG = V dP - S dT (4)
G es funcin de variables medibles: temperatura y presin
Aplicando el diferencial total a G/(RT):
(5)
Reemplazando a dG por (4) y a G por G = H - T S se llega a:
d
G
RT RT
dG
G
RT
dT
'
|
.
|
=
1
2
d
G
RT
V
RT
dP
H
RT
dT
'
|
.
|
=
2
(6)
A partir de la ecuacin (6) se pueden obtener las expresiones a T y a
P constantes:
c
c
c
c
G
RT
P
V
RT
G
RT
T
H
RT
T cte
P cte
'
|
.
|
'
|
.
|
|
|
|
=
'
|
.
|
'
|
.
|
|
|
|
=
=
=
.
.
2
(7)
(8)
c
c
c
c
G
RT
P
V
RT
G
RT
T
H
RT
T cte
P cte
'
|
.
|
'
|
.
|
|
|
|
=
'
|
.
|
'
|
.
|
|
|
|
=
=
=
.
.
2
Cuando se conoce G/(RT) en funcin de T y P, V/(RT) y H/(RT) se
obtienen por simple diferenciacin. Las propiedades restantes se
calculan con las ecuaciones que las definen, o sea:
y
S
R
H
RT
G
RT
=
U
RT
H
RT
PV
RT
=
(9) (10)
(8)
INCONVENIENTE: no hay mtodos experimentales adecuados para
determinar valores numricos de G/(RT). Por eso el tratamiento de las
propiedades termodinmicas se suele hacer mediante propiedades
residuales. stas se definen como:
M
R
= M- M
gi
(11)
donde, M:valor real de cualquier propiedad termodinmica extensiva
(tal como G, U, H, S o V)
M
gi
:valor correspondiente al gas ideal en iguales condiciones
M
R
:valor residual de la propiedad
INCONVENIENTE: no hay mtodos experimentales adecuados para
determinar valores numricos de G/(RT). Por eso el tratamiento de las
propiedades termodinmicas se suele hacer mediante propiedades
residuales. stas se definen como:
M
R
= M- M
gi
(11)
donde, M:valor real de cualquier propiedad termodinmica extensiva
(tal como G, U, H, S o V)
M
gi
:valor correspondiente al gas ideal en iguales condiciones
M
R
:valor residual de la propiedad
Segn (11), la energa residual de Gibbs es:
G
R
= G - G
gi
(12)
Anlogamente, para el volumen se tiene:
V
R
= V - V
gi
(13)
(14)
Puesto que
(15)
la relacin entre el volumen residual y el factor de compresibilidad ser:
(16)
V V
RT
P
R
=
Segn (11), la energa residual de Gibbs es:
G
R
= G - G
gi
(12)
Anlogamente, para el volumen se tiene:
V
R
= V - V
gi
(13)
(14)
Puesto que
(15)
la relacin entre el volumen residual y el factor de compresibilidad ser:
(16)
V V
RT
P
R
=
, )
V
RT
P
Z
R
= 1
V
ZRT
P
=
Para el caso especial de un gas ideal, la ecuacin (6) es:
d
G
RT
V
RT
dP
H
RT
dT
gi gi gi
'
|
.
| =
2
(17)
Restando esta ecuacin de la (6) se obtiene:
d
G
RT
V
RT
dP
H
RT
dT
R R R
'
|
.
| =
2
(18)
A partir de ella se obtienen: A partir de ella se obtienen:
c
c
c
c
G
RT
P
V
RT
G
RT
T
H
RT
R
T cte
R
R
P cte
R
'
|
.
|
'
|
.
|
|
|
|
|
=
'
|
.
|
'
|
.
|
|
|
|
|
=
=
=
.
.
2
(19)
(20)
Haciendo la resta de G = H - T S y G
gi
= H
gi
- T S
gi
se obtiene
G
R
= H
R
- T S
R
(21)
de la cual se puede despejar la entropa residual:
(22)
La energa interna de un sistema viene dada por la siguiente relacin:
S
R
H
RT
G
RT
R R R
=
Haciendo la resta de G = H - T S y G
gi
= H
gi
- T S
gi
se obtiene
G
R
= H
R
- T S
R
(21)
de la cual se puede despejar la entropa residual:
(22)
La energa interna de un sistema viene dada por la siguiente relacin:
U
RT
H
RT
PV
RT
=
Para un gas real esta expresin se convierte en:
U
RT
H
RT
Z =
(23)
(24)
mientras que para un gas ideal queda:
U
RT
H
RT
gi gi
= 1
(25)
Restando (24) y (25) tendremos la energa interna residual:
, )
U
RT
H
RT
Z
R R
= 1
(26)
Para un sistema a temperatura constante, la ecuacin (19) puede
escribirse como:
Para un sistema a temperatura constante, la ecuacin (19) puede
escribirse como:
d
G
RT
V
RT
dP
R R
'
|
.
| = (27)
G
RT
V
RT
dP (T
R R P
=
}
constante)
0
Integrando (27) desde P = 0 (G
R
/(RT) = 0 por ser gas ideal) hasta una
presin arbitraria P,
(28)
Sustituyendo (16) en (28) llegamos a:
G
RT
R P
=
}
(Z - 1)
dP
P
(T constante)
0
Si para hallar el valor de Z usamos un desarrollo virial del tipo
( )
(
Z
B T)P
RT
= 1 ( )
(
Z
B T)P
RT
= 1
donde B(T) es una funcin de la temperatura, tendremos que:
G
RT
B T P
RT
R
=
( )
(29)
A partir de (29) podemos obtener:
c
c
c
c
G
RT
T
B T P
RT
T
R
P cte
'
|
.
|
'
|
.
|
|
|
|
|
=
'
|
.
|
= .
( )
c
c
c
c
G
RT
T
B T
T
P
RT
B T P
RT
R
P cte
'
|
.
|
'
|
.
|
|
|
|
|
=
= .
( ) ( )
2
Entonces,
Para obtener una expresin de la entalpa residual hay que reemplazar
esta ltima ecuacin en (20):
Para obtener una expresin de la entalpa residual hay que reemplazar
esta ltima ecuacin en (20):
H
RT
B T)P
RT
B T)
T
P
R
R
=
( ( c
c
U
RT
B T)
T
P
R
R
=
c
c
(
S
R
B T)
T
P
R
R
=
c
c
(
Luego, podemos calcular la entropa residual con (22) y la energa
interna residual con (26):
Entonces, con la expresin del segundo coeficiente virial deducida para
alguno de los potenciales y su derivada, se pueden calcular las
propiedades termodinmicas residuales.
Se usaron dos potenciales: pozo cuadrado y pozo triangular
Sus expresiones son:
, )
c
c
c
c
B
T
B
kT
e
ER
kT 2 3
2
1 =
'
|
.
|
pozo cuadrado
pozo tringular
, )
c
c
c
c
B
T
B
kT
e
ER
kT 2 3
2
1 =
'
|
.
|
, ) , )
c
c
t
c
c
c
c
t
c
c
B
T
e
A
T
B
T
r
r
B
r
B B
e
r
B
r
B B
B
T
A Br A Br
2
2
2 3
2
2 3 4
2
2 2
2
4 6
=
'
|
.
|
|
]
|
+ + +
|
]
|
, )
c
c
c A
T
g
kT g
=
2
1 , )
c
c
c
o
B
T
kT g
=
2
1
pozo cuadrado
pozo tringular
donde y
EOS V2.00 EOS V2.00 -- DESCRIPCIN DE LAS DESCRIPCIN DE LAS
DISTINTAS OPCIONES DISTINTAS OPCIONES
EOS V2.00 EOS V2.00 -- DESCRIPCIN DE LAS DESCRIPCIN DE LAS
DISTINTAS OPCIONES DISTINTAS OPCIONES
OPCIN 1 - Lennard-Jones 6-12
Datos requeridos: temperatura [K] y presin [ATM]
Potencial empleado: Lennard-Jones 6-12 (con o sin orientacin)
Clculo del segundo coeficiente virial
Clculo del factor de compresibilidad
Grfico del segundo coeficiente virial reducido Vs. temperatura
Impresin de resultados
Datos requeridos: temperatura [K] y presin [ATM]
Potencial empleado: Lennard-Jones 6-12 (con o sin orientacin)
Clculo del segundo coeficiente virial
Clculo del factor de compresibilidad
Grfico del segundo coeficiente virial reducido Vs. temperatura
Impresin de resultados
OPCIN 2 - Lennard-Jones generalizado
Datos requeridos: temperatura [K] y presin [ATM]
Potencial empleado: Lennard-Jones n-m (con o sin orientacin)
Clculo del segundo coeficiente virial
Clculo del factor de compresibilidad
Grfico del segundo coeficiente virial reducido Vs. temperatura
Impresin de resultados
Datos requeridos: temperatura [K] y presin [ATM]
Potencial empleado: Lennard-Jones n-m (con o sin orientacin)
Clculo del segundo coeficiente virial
Clculo del factor de compresibilidad
Grfico del segundo coeficiente virial reducido Vs. temperatura
Impresin de resultados
OPCIN 3 - Esferas rgidas
Datos requeridos: temperatura [K] y presin [ATM]
Potencial empleado: esferas rgidas
Clculo del segundo al sptimo coeficiente virial
Clculo del factor de compresibilidad
Impresin de resultados
Datos requeridos: temperatura [K] y presin [ATM]
Potencial empleado: esferas rgidas
Clculo del segundo al sptimo coeficiente virial
Clculo del factor de compresibilidad
Impresin de resultados
OPCIN 4 - Pozo cuadrado
Datos requeridos: temperatura [K] y presin [ATM]
Potencial empleado: pozo cuadrado
Clculo del segundo coeficiente virial
Clculo del factor de compresibilidad
Grfico del segundo coeficiente virial reducido Vs. temperatura
Impresin de resultados
Datos requeridos: temperatura [K] y presin [ATM]
Potencial empleado: pozo cuadrado
Clculo del segundo coeficiente virial
Clculo del factor de compresibilidad
Grfico del segundo coeficiente virial reducido Vs. temperatura
Impresin de resultados
OPCIN 5 - Pozo triangular
Datos requeridos: temperatura [K] y presin [ATM]
Potencial empleado: pozo triangular
Clculo del segundo coeficiente virial
Clculo del factor de compresibilidad
Grfico del segundo coeficiente virial reducido Vs. temperatura
Impresin de resultados
Datos requeridos: temperatura [K] y presin [ATM]
Potencial empleado: pozo triangular
Clculo del segundo coeficiente virial
Clculo del factor de compresibilidad
Grfico del segundo coeficiente virial reducido Vs. temperatura
Impresin de resultados
OPCIN 6 - Kihara
Datos requeridos: temperatura [K] y presin [ATM]
Potencial empleado: Kihara
Clculo del segundo coeficiente virial
Clculo del factor de compresibilidad
Grfico del segundo coeficiente virial reducido Vs. temperatura
Impresin de resultados
Datos requeridos: temperatura [K] y presin [ATM]
Potencial empleado: Kihara
Clculo del segundo coeficiente virial
Clculo del factor de compresibilidad
Grfico del segundo coeficiente virial reducido Vs. temperatura
Impresin de resultados
OPCIN 7 - Propiedades termodinmicas
Datos requeridos: temperatura [K] y presin [ATM]
Potencial empleado: pozo cuadrado o pozo triangular
Clculo de las propiedades termodinmicas reducidas (U, S, H,
V y )
Impresin de resultados
Datos requeridos: temperatura [K] y presin [ATM]
Potencial empleado: pozo cuadrado o pozo triangular
Clculo de las propiedades termodinmicas reducidas (U, S, H,
V y )
Impresin de resultados
OPCIN 8 - Agregar compuestos
Datos requeridos: parmetros del compuesto para ese potencial
OPCIN 9 - Eliminar compuestos
Datos requeridos: cdigo del compuesto para ese potencial Datos requeridos: cdigo del compuesto para ese potencial
OPCIN 10 - SALIR
81 100 400 600
Temperatura [K]
0
50
-50
-100
-150
-200
-250
-300
B2 [cm3/mol]
Potencial
Esferas rgidas
Pozo cuadrado
Lennard-Jones 6-12
Kihara
Experimental
81 100 400 600
Temperatura [K]
0
50
-50
-100
-150
-200
-250
-300
B2 [cm3/mol]
Potencial
Esferas rgidas
Pozo cuadrado
Lennard-Jones 6-12
Kihara
Experimental
110 170 400 600
Temperatura [K]
0
100
-100
-200
-300
-400
B2 [cm3/mol]
Potencial
Esferas rgidas
Pozo cuadrado
Lennard-Jones 6-12
Kihara
Experimental
110 170 400 600
Temperatura [K]
0
100
-100
-200
-300
-400
B2 [cm3/mol]
Potencial
Esferas rgidas
Pozo cuadrado
Lennard-Jones 6-12
Kihara
Experimental
75 150 300 400 600
Temperatura [K]
0
50
-50
-100
-150
-200
-250
-300
B2 [cm3/mol]
Potencial
Esferas rgidas
Pozo cuadrado
Lennard-Jones 6-12
Experimental
75 150 300 400 600
Temperatura [K]
0
50
-50
-100
-150
-200
-250
-300
B2 [cm3/mol]
Potencial
Esferas rgidas
Pozo cuadrado
Lennard-Jones 6-12
Experimental
110 200 300 400 600
Temperatura [K]
0
100
-100
-200
-300
-400
B2 [cm3/mol]
Potencial
Esferas rgidas
Pozo cuadrado
Lennard-Jones 6-12
Experimental
110 200 300 400 600
Temperatura [K]
0
100
-100
-200
-300
-400
B2 [cm3/mol]
Potencial
Esferas rgidas
Pozo cuadrado
Lennard-Jones 6-12
Experimental