Advanced
Study of Excess Gibbs Energy for a Lattice Solution by Random Number Simulation
Study of Excess Gibbs Energy for a Lattice Solution by Random Number Simulation
Journal of the Korean Chemical Society. 2007. Aug, 51(4): 312-317
Copyright © 2007, The Korean Chemical Society
  • Received : May 12, 2007
  • Published : August 20, 2007
Download
PDF
e-PUB
PubReader
PPT
Export by style
Article
Author
Metrics
Cited by
TagCloud
About the Authors
해영 정

Abstract
난수 모의실험을 수행하여 격자 위에 분자를 무작위하게 배열하는 경우의 수의 분포를 다른 분자간의 최근린 상호작용수 N 12 에 대한 정규분포로 근사하였다. 이 분포로부터 논랜덤 혼합 격자용액의 과잉깁스에너지 GE 에 대한 근사식을 유도하였다. 이를 이용하여 여러 이성분용액의 액체-증기 상평형 계산을 하였고 기존의 식들의 계산 결과와 비교하여 보았다.
Keywords
서 론
액체용액에 대한 열역학적 설명을 하기 위하여 cell, hole, free volume 또는 격자(lattice)의 개념에 근거를 둔 여러 모델들이 1 제안되어 왔다. 특히 격자 모델의 경우 취급과 이해가 용이하여 많은 연구가 진행되어 왔으며 현재도 격자에 근거를 둔 이론들이 폭넓게 사용되고 있다.
실제용액에서는 이상적인 경우와는 달리 분자들의 혼합이 랜덤하게 이루어 지지 않으며, 이러한 논랜덤 혼합은 용액의 열역학적 성질에 많은 영향을 주고 있다. 용액에서 일어나는 분자간의 상호작용에 의한 논랜덤 혼합효과를 설명하는 대표적인 이론중의 하나가 격자의 개념에 근거를 둔 quasi-chemical 이론 2 이다. 이론의 결과가 화학평형의 관점에서 직관적으로 유추가 되는 까닭에 quasi-chemical 이론이라 불리우고 있다. 이를 바탕으로 하여 여러 모델 3 이 제시되어 왔으며, Wilson 4 이 제안한 수정된 식은 이성분액체-증기 상평형의 실험적인 값들을 잘 계산해 주고 있다. 이 후 Wilson식에 근거를 둔 여러 이론들이 제시되어 왔으며 UNIQUAC 이론 3 은 공학적으로 상당한 성과를 거두었으나 이론적인 관점에서 논랜덤 혼합효과를 정확하게 설명하고 있지는 못하다.
본 연구에서는 난수발생기를 이용하여 5가지 종류의 격자에서 입자를 무작위하게 배열하는 모의실험을 수행하였다. 이렇게 해서 입자를 격자에 배열하는 경우의 수의 분포를 다른 분자간의 최근린 상호작용수 N 12 에 대한 정규분포로 근사하였고, 논랜덤혼합 격자용액에서의 과잉깁스에너지 GE 에 대한 근사식을 유도하였다. 이 식을 사용하여 10가지 이성분계의 액체-증기에 대한 상평형을 계산하였고, Redlich-Kister식, Van Laar식, Wilson식 등 3 의 계산 결과와 비교하여 보았다.
- 난수 모의실험을 통한 N12분포
모의실험은 난수발생기를 통하여 생성된 난수를 이용하여 N 1 개의 입자-1과 N 2 개의 입자-2를 N 1 + N 2 개의 격자에 무작위로 배열하는 것이다. 1차원 선형격자, 2차원 평방격자, 입방격자, 체심방격자, 면심방격자의 5가지 격자에 대하여 모의실험이 이루어 졌다. 각 격자에서 최근린 분자수 z는 2, 4, 6, 8, 12이다. 1차원 선형격자에 대하여서는 정확한 해가 나와 있으며 비교를 위하여 실험에 포함시켰다. 1차원 선형격자의 경우 N 이 매우 클 때 배열에 대한 정확한 확률밀도함수를 입자-1과 입자-2의 최근린 상호작용수 N 12 에 대하여 쓰면 다음과 같다. 5
PPT Slide
Lager Image
식 (1)의 확률분포는 초기하분포에 해당하여 X에 대한 평균과 표준편차는 다음과 같다. 6
PPT Slide
Lager Image
PPT Slide
Lager Image
모의실험은 입자의 개수 N 은 대략 10000~80000 범위, 입자-1의 몰분율 x 1 = N 1 /N은 0.1, 0.2, 0.3, 0.4, 0.5에 대하여 수행하였다. 무작위 배열에 대한 시행횟수는 모두 10 6 으로 하였다. . 1 에서 보는 바와 같이 입자의 수가 클 경우, X 에 관한 확률분포는 정규분포에 수렴하게 되는 데, 이는 표본의 수가 매우 큰 이산분포의 경우 중심극한정리 6 에 의하여 정규분포에 수렴한다는 사실에 부합되는 것이다. 그리고 . 2 에서 보는 바와 같이 X 의 평균과 편차, x 1 사이에는 다음과 같은 비례관계가 성립함을 알 수 있다.
PPT Slide
Lager Image
Probability distribution of X/N2 from 106 trials when x1=0.3, N=10648 for a cubic lattice.
PPT Slide
Lager Image
, X/N2 vs x1 from 106 trials when N=10648 for a cubic lattice.
PPT Slide
Lager Image
PPT Slide
Lager Image
모의실험을 통해 얻어진 식 (4), (5)의 비례상수는 1 , 2 에 수록되어 있다. 1 , 2 에서 Linear, Square, Cubic, BCC, FCC 는 각각 1차원 선형격자, 2차원 평방격자, 입방격자, 체심입방격자, 면심입방격자를 나타내며 입자의 개수 N 은 1차원 선형격자를 제외하고는 대략적인 값이다. 1 , 2 의 값과 1차원 선형격자에 대한 정확한 식 (2), (3)과의 비교를 통하여 식 (4), (5)의 비례상수는 격자의 차원이나 구조와 무관하게 각각 z /2,
PPT Slide
Lager Image
임을 알 수 있다. 따라서, N 이 매우 커지면 무작위 배열에 대한 X 의 확률분포는 다음의 평균과 편차를 갖는 정규분포에 수렴하게 된다.
Slope of vsx1for 5 different lattice structures
PPT Slide
Lager Image
N in the 1st column are approximate values except for a linear lattice
Slope ofvsx1for 5 different lattice structures
PPT Slide
Lager Image
N in the 1st column are approximate values except for a linear lattice
PPT Slide
Lager Image
PPT Slide
Lager Image
식 (6), (7)은 N 이 매우 클 때 무작위 배열의 경우의 수에 대한 확률분포는 격자의 차원이나 구조와는 무관하게 몰분율과 최근린 입자수 z 에 의존한다는 것을 보여 주고 있다.
- 정규분포를 이용한 과잉깁스에너지 GE의 근사
본 연구에서는 X 에 대한 확률분포를 정규분포로 근사하기로 한다.
그러면, 입자배열에 대한 경우의 수 Ω는 다음과 같이 된다.
PPT Slide
Lager Image
격자용액에 대한 통계역학적인 분배함수 Q 는 다음과 같다. 5
PPT Slide
Lager Image
식 (9) 우변에서 총에너지 E ( X )는 입자- i , j 간의 상호작용에너지 ε ij 를 사용하여 다음과 같이 표시된다.
PPT Slide
Lager Image
PPT Slide
Lager Image
Q 는 다음과 같이 식 (9) 우변의 최대항으로 근사할 수 있다. 5
PPT Slide
Lager Image
식 (12)에서 X *는 다음을 만족하는 X 값이다.
PPT Slide
Lager Image
식 (8), (10), (13)으로부터
PPT Slide
Lager Image
식 (8), (12), (14)으로부터
PPT Slide
Lager Image
N 1 , N 2 가 매우 클 경우 식 (15)은 다음과 같이 근사할 수 있다.
PPT Slide
Lager Image
식 (16)에 식 (6), (7)을 대입하여 정리하면
PPT Slide
Lager Image
격자용액에서 혼합부피변화 Δ Vmix 를 무시하고 통계열역학의 공식 3 을 사용하여 식 (17)으로부터 혼합깁스에너지 Δ Gmix GE 는 다음과 같이 표시된다.
PPT Slide
Lager Image
PPT Slide
Lager Image
식 (18), (19)에서
PPT Slide
Lager Image
식 (20)에서 최근린 분자수 z 를 12로 고정시키면 식 (19)은 1-매개변수의 식이 됨을 알 수 있다. B 는 랜덤혼합의 경우에는 나타나지 않으므로 논랜덤효과를 나타내는 항이라는 것을 알 수 있다. 그리고 식 (19)은 다음의 3-매개변수 Redlich-Kister식에서 두 번째 매개변수 B '=0인 경우에 해당된다는 것을 알 수 있다.
PPT Slide
Lager Image
식 (19)을 보면 x 1 =0.5에 대하여 대칭적인 관계를 나타내는 데 이것은 분자의 크기,구조등이 고려되지 않았기 때문이다. 그래서 본 연구에서는 몰분율 x 1 대신 분자의 구조, 크기등을 고려한 식 (22)의 분율을 사용하기로 한다.
PPT Slide
Lager Image
식 (22)에서 r 1 , r 2 는 고분자용액을 설명하는 Flory-Huggins 격자이론 7 에서 고분자-1과 2의 사슬길이를 나타내는 매개변수로 저분자용액에서는 분자 크기, 구조등을 고려한 실험적 매개변수의 의미를 가지고 있다.
식 (19), (22)으로부터
PPT Slide
Lager Image
즉 식 (23)은 z 를 12로 놓아 A , r 의 2 매개변수를 갖는 식이 된다. A 는 에너지, B 는 논랜덤혼합에 의한 효과를 나타내며, r 은 분자의 크기, 구조등에 의한 효과를 나타내는 매개변수라 할 수 있다. 그런데 식 (23)은 근사적인 식 (8)으로부터 유도되었으므로 A , B 를 각각 독립적인 반실험적 매개변수로 생각하여 A , B , r 의 3-매개변수 식으로 취급하는 것이 더 적합하다고 볼 수 있다.
계산 및 결과
본 연구에서는 CH 3 NO 2 +CCl 4 8 , EtOH +chloroform 9 , MeOH+benzene 10 , benzene+methylcyclohexane 11 , cyclohexane+toluene 11 , cyclopentane+benzene 11 , propanol+benzene 12 , EtOH+benzene 12 , methylcyclopentane+benzene 11 , methylcyclopentane+toluene 11 의 10가지 이성분계 용액에 대한 상평형 계산을 하였다. 계산에는 식 (23)을 포함하여 Redlich-Kister식, Van Laar식, Wilson식을 사용하였다. Redlich-Kister식, Van Laar식, Wilson식에 대하여는 참고문헌 3 에 잘 정리되어 있다. 그리고 각 식에 나타나는 매개변수를 최적화하기 위해서 증기압에 대한 실험값과 계산값의 근평균평방오차 식 (24)을 최소화시키는 Barker 13 의 방법을 사용하였다.
PPT Slide
Lager Image
식 (24)에서 증기압 계산값 Pcalc 은 다음의 식을 사용하여 계산한다.
PPT Slide
Lager Image
γ 1 , γ 2 는 활동도 계수로 과잉깁스에너지 GE 로부터 계산된다. 3
식 (25)에서
relative error % in vapor pressure calculated for 10 binary mixtures
PPT Slide
Lager Image
relative error % in vapor pressure calculated for 10 binary mixtures
Δycalculated for 10 binary mixtures
PPT Slide
Lager Image
Δy calculated for 10 binary mixtures
PPT Slide
Lager Image
PPT Slide
Lager Image
식 (26)에서 v 1 , v 2 는 액체성분의 몰부피로 modified Rackett식 14 을 사용하였으며, P 1 s , P 2 s 는 각 성분의 평형증기압으로 참고문헌14에 수록되어 있다. B 11 , B 22 , B 12 는 각 성분의 비리얼 계수로 Hayden과 O’Connell 15 의 방법을 사용하여 상당히 정확히 계산되며 계산에 필요한 식과 매개변수들은 참고문헌 16 에 수록되어 있다.
계산 결과는 3 , 4 에 나타내었다. 3 에는 증기압에 대한 상대오차%, 4 에는 증기상 몰분율에 대한 근평균평방오차 Δ y 가 수록되어 있다. 3 , 4 에서 RK-2, RK-3는 각각 매개변수가 2개, 3개인 Redlich-Kister식을 나타내며 Eqn. (23)-B는 식 (23)에서 B를 독립된 매개변수로 취급하여 총 3개의 매개변수를 사용한 경우를 나타낸다. 3 , 4 에서 보는 바와 같이 3개의 매개변수를 사용하는 식 (23)이 가장 좋은 결과를 보여주고 있으며, 전반적으로는 본 연구의 식 (23)은 Redlich-Kister식보다 정량적으로 약간 개선된 결과를 보여 주고 있다. 매개변수를 2개만 사용하는 경우에는 Wilson식이 가장 좋은 결과를 보여 주고 있음을 알 수 있다.
결 론
본 연구에서는 난수 모의실험을 수행하여 격자용액에서 분자의 무작위 배열에 대한 경우의 수의 분포를 N 12 에 대한 정규분포로 근사하였다. 이로부터 논랜덤 혼합 격자용액의 과잉깁스에너지 GE 에 대한 반실험적인 공식을 구하였으며, 식에서 나타나는 3개의 매개변수는 각각 에너지, 논랜덤혼합, 분자 구조 및 크기의 효과를 나타내는 정성적인 의미를 가지고 있다. 계산 결과는 기존의 널리 사용되던 Redlich-Kister식 보다 다소 개선된 결과를 보여줌을 알 수 있었다. 그런데 서로 다른 분자 간에 발생할 수 있는 특정상호작용등에 의한 효과를 고려하지 않았으므로 이런 부분을 고려하게 되면 보다 개선된 결과를 얻을 수 있을 것으로 생각된다.
본 연구는 2006년도 덕성여자대학교 교내연구비 지원으로 수행되었음.
References
Hirschfelder J. O. , Curtiss C. F. , Bird R. B. 1954 The Molecular Theory of Gases and Liquids John Wiley & Sons
Guggenheim E. A. 1952 Mixtures Clarendon Press Oxford
Prausnitz J. M. , Lichtenthaler R. N. , de Azevedo E. G. 1986 Molecular Thermodynamics of Fluid-Phase Equilibria 2nd Ed. Prentice-Hall
Wilson J. M. 1964 J. Am. Chem. Soc. 86 127 -    DOI : 10.1021/ja01056a002
Hill T. L. 1987 An Introduction to Statistical Thermodynamics Dover Publications
Mendenhall W. , Sincich T. 1995 Statistics for Engineering and the Sciences 4th Ed. Prentice-Hall
Flory P. J. 1953 Principles of Polymer Chemistry Cornell Univ. Press
Brown I , Smith F. 1957 Aust. J. Chem. 10 423 -    DOI : 10.1071/CH9570423
Scatchard G. , Raymond C. 1938 J. Am. Chem. Soc. 60 1279 -
Scatchard G. , Wood S. , Mochel J. 1946 J. Am. Chem. Soc. 68 1957 -    DOI : 10.1021/ja01214a024
Myers H. S. 1956 Ind. Eng. Chem. 48 1104 -    DOI : 10.1021/ie50558a041
Wehe A. , Coates J. 1955 AICHE J. 1 241 -    DOI : 10.1002/aic.690010218
Barker J. A. 1953 Aust. J. Chem. 6 207 -    DOI : 10.1071/CH9530207
Reid R. C. , Prausnitz J. M. , Sherwood T. K. 1988 The Properties of Gases and Liquids 4th Ed. McGraw-Hill
Hayden J. G. , O’Connell J. P. 1975 Ind. Eng. Chem., Process Des. Dev. 14 209 -    DOI : 10.1021/i260055a003
Prausnitz J. M. 1980 Computer calculations for multicomponent vapor-liquid and liquid-liquid equilibria Prentice-Hall