알루미늄 규산염(Al2SiO5) 광물은 온도와 압력 환경에 따라 남정석, 홍주석, 규선석으로 상전이가 일어나는 동질이상(polymorph)으로 변성암의 변성정도를 유추하는 데 사용되는 중요한 광물이다. 이번 연구에서는 고전분자동력학 시뮬레이션(classical molecular dynamic simulations)과 양자역학 계산방법인 밀도범함수이론(density functional theory)을 이용하여 압력에 따른 알루미늄 규산염 광물의결정구조와 엔탈피를 계산하고, 상대적인 안정성을 평가하였다. 격자상수 계산결과, 분자동력학과 밀도범함수이론 계산 모두 압력에 따라 부피가 줄어드는 기존의 실험결과와 동일한 경향을 보였다. 특히, 밀도범함수이론으로 얻어진 격자상수는 실험과 약 1% 이내의 오차로 매우 정확한 결과를 보였다. 그러나 엔탈피 계산 결과, 분자동력학에서는 압력에 따른 엔탈피의 변화가 거의 없어 광물 간 안정성이 역전되는 상전이 압력을 구할 수 없었다. 밀도범함수이론 계산 결과는 실험과 동일한 경향을 보여주었지만, 전자의 교환-상관 관계를 나타내는 범함수에 따라 상전이 압력이 크게 다른 결과를 보여주었다. 밀도범함수이론 계산 결과는 결정구조와 엔탈피에 대해서 모두 높은 수준의 정확도를 보여주지만, 동질이상의 상도표 작성에는 정량적으로 큰 오차를 보여주었다.

점토광물의 물리화학적 특성에 대한 분자 또는 원자 스케일의 연구 중요성이 강조되고 있다. 그러나 실험만으로는 광물의 미시적 현상을 이해하기 어려운 경우가 많다. 특히 2:1 점토광물 팔면체에 존재하는 수산기(hydroxyl)가 금속 양이온 흡착과정에 큰 역할을 한다는 가정은 X-ray를 이용하는 실험만으로는 명확하게 테스트하기 어렵다. 이번 논문에서는 점토광물 내의 수산기 연구에 대한 전산광물학(computational mineralogy) 이용 가능성에 대하여 조사하였다. 점토광물의 기본구조인 팔면체 층만으로 구성된 광물, 1:1 구조를 갖는 광물, 2:1 구조를 갖는 광물 중 대표적인 이팔면체 및 삼팔면체 층상규산염 광물을 선별하여 구조최적화를 실시하였다. 분자역학적(molecular mechanics) 계산과 양자역학적(quantum mechanical) 계산 모두 실험값의 격자상수(lattice parameters)를 잘 재현할 수 있었다. 그러나, 사면체층과 팔면체의 구조적 뒤틀림(structural distortion) 등 광물 내부구조를 기존 실험결과와 비교했을 때, 양자역학적 계산결과가 분자역학적 방법을 이용한 결과 보다 더 낮은 오차를 보였다. 파이로필라이트(pyrophyllite) 수산기가 (001)면과 이루는 각은, 수산기의 H(proton)과 사면체의 Si 양이온 간의 척력으로 결정되는데, 양자역학적 방법은 약 25-26˚로 예측하였고, 분자역학적 방법은 약 35˚ 정도로 양자역학계산 결과와 무려 10˚의 큰 차이를 보였다. 전산광물학은 점토광물 구조연구에 신뢰성이 매우 높은 연구방법으로 양이온 흡착과정 중 수산기의 역할 규명에 사용될 수 있다.