고압 환경에서 규산염 용융체의 원자 구조에 대한 정보는 지구 내부 마그마의 열전도율이나 주변 암석과의 원소 분배계수와 같은 이동 물성을 이해하는 단서를 제공한다. 규소의 전자 구조는 규산염 다면체 주 변의 산소 원자 분포와 연관성을 가질 것으로 예상되나, 이 사이의 상관관계가 명확하게 밝혀져 있지 않다. 본 연구는 SiO2의 고밀도화에 따른 규소의 전자 구조 변화의 미시적인 기원을 규명하기 위해 SiO2 동질이상 의 규소 부분 상태 밀도와 L3-edge X-선 흡수분광분석(X-ray absorption spectroscopy; XAS) 스펙트럼을 계산 하였다. 규소의 전도 띠 영역에서 전자 구조는 결정 구조에 따라서 변화하였다. 특히 d-오비탈은 108, 130 eV 영역에서 배위 환경에 따른 뚜렷한 차이를 보였다. 계산된 XAS 스펙트럼은 규소 전도 띠의 s,d-오비탈에서 기인하는 피크를 보였으며, 결정 구조에 따라 s,d-오비탈과 유사한 양상으로 변화했다. 계산된 석영의 XAS 스펙트럼은 SiO2 유리의 X R S 실험 결과와 유사하였으며 규소 주변 원자 환경이 비슷하기 때문으로 생각된 다. XAS 스펙트럼을 수치화한 무게 중심 값은 Si-O 결합 거리와 밀접한 상관관계를 가지며 이로 인하여 고 밀도화 과정에서 체계적으로 변화한다. 본 연구의 결과는 Si-O 결합 거리에 민감한 규소 L2,3-edge XRS가 규 산염 유리 및 용융체의 고밀도화 기작을 규명하는 과정에서 유용하게 적용될 수 있음을 지시한다.
SiO2는 지각과 맨틀을 구성하는 풍부한 물질로 고압 상태의 SiO2 원자구조를 결정짓는 전자구 조적 특성에 대한 상세한 이해는 지구 내부의 탄성과 열역학적 성질에 대한 통찰을 제공한다. SiO2처 럼 경원소(low-z)로 이루어진 지구 물질의 고압상 전자구조는 in situ 고압 XRS (x-ray Raman scattering) 실험을 통해 연구되어 왔다. 하지만 기존의 고압 실험 방법으로는 물질의 국소 원자구조와 XRS 스펙트럼 간 상관관계를 밝히는데 한계가 있다. 이를 극복하고 더 높은 압력에서 존재하는 SiO2 에 대한 XRS 정보를 얻기 위해 밀도 범함수 이론(density functional theory; DFT)에 기반을 둔 제1원 리(ab initio) 계산법을 이용한 XRS 스펙트럼 계산 연구들이 진행되고 있다. 비탄성 X-선 산란에 의하 여 원자핵 주변 1s 오비탈에 만들어지는 전자-정공(core-hole)은 경원소 물질의 국소 전자구조에 크게 영향을 미치기 때문에 O K-edge XRS 스펙트럼 형태를 계산할 때 중요하게 고려해야 한다. 본 연구에서는 온-퍼텐셜 선형보충파(full-potential linearized augmented plane wave; FP-LAPW) 방법 론에 기반하는 WIEN2k 프로그램을 사용하여 α-quartz, α-cristobalite 그리고 CaCl2-구조를 갖는 SiO2 에 대한 O 원자 전자 오비탈의 부분 상태밀도(partial density of states; PDOS)와 O K-edge XRS 스펙 트럼을 계산하였다. 또한, CaCl2-구조를 갖는 SiO2의 O 원자 PDOS의 전자-정공 효과의 적용 여부에 따른 차이를 비교하여, 원자핵 부근 전자구조 변화에 따른 PDOS의 피크 세기와 위치 변화가 크게 나 타났다는 사실을 확인할 수 있었다. 또한 계산된 각 SiO2 구조의 O K-edge XRS 스펙트럼이 각 SiO2 구조에서 계산된 O 원자의 p* 오비탈의 PDOS 결과와 매우 유사한 형태를 갖고 있음을 확인하였다. 이는 O K-edge XRS 스펙트럼이 갖는 대부분의 특징적인 피크들이 O 원자의 점유 1s 오비탈에서 2p* 오비탈로의 전자전이에 기인하기 때문이다. 본 연구의 결과는 SiO2에 대한 정확한 O K-edge XRS 스 펙트럼을 계산하는데 있어 전자-정공 효과를 고려해야 한다는 사실을 보여준다. 또한, 실험적으로는 재현이 어려운 고압 환경에 존재하는 CaCl2-구조를 갖는 SiO2 (~63 GPa)에 대한 O K-edge XRS 스펙 트럼 계산을 통해, 제1원리 계산이 고압상 물질의 물성 연구에 이용될 수 있다는 사실을 보여준다.
지구시스템 이해에 중요한 지구 내부 맨틀 물질의 거시적인 성질을 이해하기 위해서는 고압상태의 Mg-규산염 결정질 및 비정질 물질에 대한 원자구조와 그에 수반하는 전자구조에 대한 이해가 필요하다. 근래에 in-situ 고압 실험의 어려움을 피하여 고압환경에 존재하는 지구물질의 원자구조와 그 전자구조를 규명하기 위한 방법론으로서 밀도 범함수 이론에 기반을 둔 양자화학계산이 많이 이용되고 있다. 본 연구에서는 FP-LAPW (full-potential linearized augmented plane wave) 방법론을 이용하는 WIEN2k 프로그램을 통하여 25 GPa와 120 GPa의 MgSiO3 페로브스카이트(Pv)의 전자 오비탈의 PDOS (partial density of states)와 O 원자 K-전자껍질 ELNES (energy-loss near-edge structure) 스펙트럼을 계산하였다. 두 압력 조건의 MgSiO3 Pv에 대하여 계산된 전자 오비탈의 PDOS와 O원자 K-전자껍질 ELNES 스펙트럼은 뚜렷한 차이를 보이고 있었다. 이와 같은 결과는 MgSiO3 Pv에서 압력 증가에 의한 Si 원자 배위수의 변화가 나타나지 않더라도 Si-O 결합거리, O-O거리, Mg-O거리와 같은 O 원자 주변 국소 원자구조의 변화가 O원자 주변 전자구조에 뚜렷한 영향을 미칠 수 있음을 의미한다. 본 연구의 결과는 MgSiO3 결정질 및 비정질 물질의 압력에 의한 전자구조 변화의 미시적 기원을 이해하고 더욱 나아가 다양한 지구물질의 압력에 의한 원자구조 변화와 그에 수반되는 전자구조 변화의 관계를 이해하는데 많은 도움을 줄 수 있을 것이다.