점토광물은 대기부터 지하수에 이르는 크리티컬존(critical zone) 영역에서 금속과 탄소 순환을 결정짓는 역할을 한다. 계산광물학 연구방법 중에 하나인 분자동역학(molecular dynamics) 시뮬레이션 은 지구물질을 원자단위로 계산하기 때문에, 점토광물의 물리화학적 현상들에 대해 원자수준의 자세 한 정보를 제공할 수 있다. 이번 연구에서는 clayFF 힘 장(force field)을 사용한 분자동역학 시뮬레이 션을 이팔면체 점토광물인 깁사이트(gibbsite, Al(OH)3), 카올리나이트(kaolinite, Al2Si2O5(OH)4), 파이 로필라이트(pyrophyllite, Al2Si4O10(OH)2)에 적용하여 300 K, 1기압조건에서 각 광물이 가지는 격자상 수와 원자간 거리를 계산하고 실험결과와 비교하였다. 더불어 수산기의 방향성 및 수소결합의 양상 그리고 파워스펙트럼(power spectrum)을 추가적으로 계산하였다. 계산 결과, 격자상수는 기존의 실험 결과와 약 0.1~3.7% 미만의 오차율을 보였으며, 원자간 거리는 실험결과와 약 5% 미만의 차이를 가 졌다. 깁사이트나 카올리나이트의 팔면체층 표면에 존재하는 수산기가 가지는 신축진동파수(stretching vibrational wavenumber)는 실험값 보다 약 200-300 cm-1 높게 계산되었지만, 팔면체층 표면에 존재하 는 수산기들의 상대적 크기의 경향은 기존 실험 결과와 일치하였다. 팔면체층 표면의 수산기가 (001) 면과 이루는 각도도 기존 실험결과와 상당부분 일치한 반면에 내부 수산기의 경우는 다소 차이를 보 였다. ClayFF를 사용한 분자동역학 시뮬레이션 연구 방법은 이팔면체 점토광물 표면 내(층간) 금속이 온 흡착에 대한 수산기의 역할을 규명하는데 유용한 연구방법이 될 수 있음을 시사한다.

점토광물은 대기부터 지하수에 이르는 크리티컬존(critical zone) 영역에서 금속과 탄소 순환을 결정짓는 역할을 한다. 계산광물학 연구방법 중에 하나인 분자동역학(molecular dynamics) 시뮬레이션 은 지구물질을 원자단위로 계산하기 때문에, 점토광물의 물리화학적 현상들에 대해 원자수준의 자세 한 정보를 제공할 수 있다. 이번 연구에서는 clayFF 힘 장(force field)을 사용한 분자동역학 시뮬레이 션을 이팔면체 점토광물인 깁사이트(gibbsite, Al(OH)3), 카올리나이트(kaolinite, Al2Si2O5(OH)4), 파이 로필라이트(pyrophyllite, Al2Si4O10(OH)2)에 적용하여 300 K, 1기압조건에서 각 광물이 가지는 격자상 수와 원자간 거리를 계산하고 실험결과와 비교하였다. 더불어 수산기의 방향성 및 수소결합의 양상 그리고 파워스펙트럼(power spectrum)을 추가적으로 계산하였다. 계산 결과, 격자상수는 기존의 실험 결과와 약 0.1~3.7% 미만의 오차율을 보였으며, 원자간 거리는 실험결과와 약 5% 미만의 차이를 가 졌다. 깁사이트나 카올리나이트의 팔면체층 표면에 존재하는 수산기가 가지는 신축진동파수(stretching vibrational wavenumber)는 실험값 보다 약 200-300 cm-1 높게 계산되었지만, 팔면체층 표면에 존재하 는 수산기들의 상대적 크기의 경향은 기존 실험 결과와 일치하였다. 팔면체층 표면의 수산기가 (001) 면과 이루는 각도도 기존 실험결과와 상당부분 일치한 반면에 내부 수산기의 경우는 다소 차이를 보 였다. ClayFF를 사용한 분자동역학 시뮬레이션 연구 방법은 이팔면체 점토광물 표면 내(층간) 금속이 온 흡착에 대한 수산기의 역할을 규명하는데 유용한 연구방법이 될 수 있음을 시사한다.

FORTRAN program PHYLS was developed to model the structures of 2:1 1M and 2M1 phyllosilicates on the basis of geometrical analyses. Input to PHYLS requires the chemical composition and d(001) spacing of the mineral. The output from PHYLS consists of the coordinates of the crystallographically independent sites in the unit cell, and such structural parameters as the cell dimensions, interaxial angle, cell volume, interatomic distances, and deformation angles of the polyhedra. PHYLS can generate these structural details according to the user's choice of space group and cation configuration. User can choose one of such space groups as C2/m, C2,and C2/c and such cation configurations as random and ordered tetrachedral/octahedral cation configurations. PHYLS simulated the structures of dioctahedral and trioctahedral phyllosilicates having random tetrahedral cation configuration fairly close to the reported experimentally determined structures. In contrast, the simulated structures for ordered tetrahedral cation configurations showed greater deviation from the experimentally determined structures than those for random configurations. However, if the cations were partially ordered and the sizes of the tetrahedra became similar, the simulated PHYLS may be helpful in various investigations on the relationships between structures and physicochemical properties of the phyllosilicates.