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