화학공학소재연구정보센터
홈 로그인 로그아웃 연락처 사이트맵
HWAHAK KONGHAK, Vol.40, No.5, 558-564, October, 2002
Export Citation
분자동력학과 몬테카를로의 혼성 분자시뮬레이션에 의한 다가전해질 복잡유체의 미세구조 특성 연구
Study on Microstructural Properties of Polyelectrolyte Complex Fluids by Employing Hybrid Molecular Simulations of Molecular Dynamics and Monte Carlo Scheme
이현수, 전명석†
  • 한국과학기술연구원 복잡유체연구팀, 서울 130-650
Hyun Su Lee, and Myung-Suk Chun†
  • Complex Fluids Team, Korea Institute of Science and Technology, Seoul 130-650, Korea
E-mail:
초록
분자동력학(molecular dynamics)과 몬테카를로(Monte Carlo)의 혼성(hybrid) 알고리즘을 이용하여 harmonic spring 사슬과 하전된 단분자(monomer)들 사이의 Debye-Huckel potential로 표현된 다가전해질(polyelectrolyte)의 미세구조 특성을 규명하였다. Polyelectrolyte의 단일 사슬에 대해서 기존에 제안된 분자동력학과 몬테카를로의 혼성 알고리즘을 개선하여, polyelectrolyte의 conformation 특성인 끝단거리(end-to-end distance), 선회반경(radius of gyration), 그리고 구조인자(structure factor)를 한층 효율적으로 구할 수 있는 새로운 모사 기법을 개발하였다. 본 연구에서 새로이 개발된 알고리즘에 의한 모사 결과와 기존 문헌값과의 비교를 통해 해석 방법의 타당성과 정확성을 검증하였다. 끝단거리와 선회반경간 평균제곱의 비와 구조인자의 계산을 통해, polyelectrolyte 단분자의 하전율(change fraction)을 대변하는 Bjerrum length λB와 용액의 이온화세기에서 결정되어 용매의 screening 정도를 나타내는 Debye length k(-1)가 증가함에 따라 사슬이 신장되는 거동을 확인하였다.
The microstructural properties of charged polyelectrolytes described by a potential model regarding both the harmonic springs and Debye-Huckel interaction were investigated by employing a hybrid scheme of molecular dynamics(MD) and Monte Carlo(MC) simulations. Based on the previous hybrid scheme, a novel hybrid scheme has been developed in the present study, with which computational efforts are effectively reduced. We present the conformational properties such as end-to-end distance, radius of gyration and structure factor. It is evident that the simulation results of the present study agree well with the previously reported results. The elongation behavior of the polyelectrolyte chain can successfully be observed by the calculations of the characteristic mean square ratio of end-to-end distance to radius of gyration as well as the structure factor. As the Bjerrum length λB and the Debye length k(-1) increase, the polyelectrolyte chain becomes elongated.
Keywords:Molecular Dynamics;Monte Carlo Simulation;Polyelectrolyte;Complex Fluid;Structure Factor
  1. Kremer K, Binder K, Comput. Phys. Rep., 7, 259 (1988) 
  2. Binder K, Ciccotti G, "Monte Carlo and Molecular Dynamics of Condensed Matter Systems," Societa Italiana Di Fisica, Bologna (1996)
  3. Chun MS, Park OO, Yang SM, J. Colloid Interface Sci., 161, 247 (1993) 
  4. Park PJ, Chun MS, Kim JJ, Macromolecules, 33(23), 8850 (2000) 
  5. Micka U, Kremer K, Phys. Rev. E, 54, 2653 (1996) 
  6. Lyulin AV, Dunweg B, Borisov OV, Darinskii AA, Macromolecules, 32(10), 3264 (1999) 
  7. Khokhlov AR, J. Phys. A, 13, 979 (1980) 
  8. Dobrynin AV, Rubinstein M, Obukhov SP, Macromolecules, 29(8), 2974 (1996) 
  9. Micka U, Kremer K, Europhys. Lett., 38, 279 (1997) 
  10. Micka U, Holm C, Kremer K, Langmuir, 15(12), 4033 (1999) 
  11. Odijk T, J. Polym. Sci. B: Polym. Phys., 15, 477 (1977)
  12. Skolnick J, Fixman M, Macromolecules, 10, 944 (1977) 
  13. Barrat JL, Joanny JF, Europhys. Lett., 24, 333 (1993)
  14. Ha BY, Thirumalai D, Macromolecules, 28(2), 577 (1995) 
  15. Bratko D, Dawson KA, J. Chem. Phys., 99, 5352 (1993) 
  16. Sadus RJ, "Molecular Simulation of Fluids: Theory, Algorithm and Object-Orientation," Elsevier, New York (1999)
  17. Doi M, Edward SF, "The Theory of Polymer Dynamics," Oxford Univ. Press, New York (1986)
  18. Hunter RJ, "Foundations of Colloid Science," 2nd Ed., Oxford Univ. Press, Oxford (2001)