The Effect of Lattice Topology on Benzyl Alcohol Adsorption on Kaolinite Surfaces: Quantum Chemical Calculations of Mulliken Charges and Magnetic Shielding Tensor

In order to have better insights into adsorption of organic molecules on kaolinite surfaces, we performed quantum chemical calculations of interaction between three different model clusters of kaolinite siloxane surfaces and benzyl alcohol, with emphasis on the effect of size and lattice topology of the cluster on the variation of electron density and magnetic shielding tensor. Model cluster 1 is an ideal silicate tetrahedral surface that consists of 7 hexagonal rings, and model cluster 2 is composed of 7 ditrigonal siloxane rings with crystallographically distinct basal oxygen atoms in the cluster, and finally model cluster 3 has both tetrahedral and octahedral layers. The Mulliken charge analysis shows that siloxane surface of model cluster 3 undergoes the largest electron density transfer after the benzyl alcohol adsorption and that of model cluster 1 is apparently larger than that of model cluster 2. The difference of Mulliken charges of basal oxygen atoms before and after the adsorption is positively correlated with hydrogen bond strength. NMR chemical shielding tensor calculation of clusters without benryl alcohol shows that three different basal oxygen atoms (O3, O4, and O5) in model cluster 2 have the isotropic magnetic shielding tensor as $228.2{\pm}3.9,\;228.9{\pm}3.4,\;and\;222.3{\pm}3.0ppm$, respectively. After the adsorption, the difference of isotropic chemical shift varies from 1 to 5.5 ppm fer model cluster 1 and 2 while model cluster 2 apparently shows larger changes in isotropic chemical shift. The chemical shift of oxygen atoms is also positively correlated with electron density transfer. The current results show that the adsorption of benzyl alcohol on the kaolinite siloxane surfaces can largely be dominated by a weak hydrogen bonding and electrostatic force (charge-charge interaction) and demonstrate the importance of the cluster site and the lattice topology of surfaces on the adsorption behavior of the organic molecules on clay surfaces.

캐올리나이트 규산염 층과 벤질알코올의 반응에 대한 양자화학계산에서 결정학적 위상이 멀리켄 전하와 자기 차폐 텐서에 미치는 영향

본 연구에서는 점토광물 표면 클러스터의 크기와 결정학적 위상이 전자 밀도와 자기 차폐 텐서에 미치는 영향을 살펴보기 위하여 캐올리나이트 규산염 층을 대표하는 세 개의 서로 다른 위상의 모델 클러스터와 벤질 알코올과의 상호작용에 대해 다양한 수준의 양자화학 계산을 수행하였다. 모델 클러스터 1은 단순화된 7개의 규산염 고리로 이루어졌고, 모델 클러스터 2는 결정학적 위상을 가진 7개의 규산염 고리로 이루어졌으며, 모델 클러스터 3은 세 개의 규산염 고리와 팔면체 고리로 이루어져 있다. 멀리켄 전하 계산 결과 벤질 알코올과의 반응 후의 상대적인 전자 밀도 이동의 크기는 모델클러스터 3의 사면체 쪽 > 모델 클러스터 1 > 모델 클러스터 2 > 모델 클러스터 3의 팔면체 쪽의 순으로 계산되었다. 또한 벤질 알코올과 강한 수소 결합을 하는 원자들의 전자 밀도 이동이 상대적으로 크다 벤질 알코올 흡착 전에 대한 사면체 표면 원자들의 자기 차폐 텐서 결과는 결정학적 위상을 고려하지 않은 경우 표면 중심으로부터의 거리가 비슷한 산소들끼리 유사한 등방 자기 차폐 텐서 값들을 갖고, 결정학적 위상을 고려한 경우는 결정학적으로 서로 다른 산소 자리(O3, O4, O5)에 대해 각각 $228.2{\pm}3.9,\;228.9{\pm}3.4,\;222.3{\pm}3.0ppm$으로 계산되었다. 흡착 전후의 산소 원자의 화학 차폐의 차이는 알코올과 근접한 산소들에서 약 $1{\sim}5.5ppm$ 정도의 변화가 나타나며 이러한 변화는 최근의 고분해능 이차원 핵자기공명분광 분석을 이용하면 실험으로 관찰할 수 있을 것으로 예상된다. 또한 모델 클러스터 2의 화학 차폐의 변화는 모델 클러스터 1보다 상대적으로 큰 특징을 보인다. 전자밀도 이동과 화학 차폐의 변화는 약한 양의 상관관계를 가진다. 이러한 결과들은 캐올리나이트 규산염 사면체 층과 벤질 알코올이 약한 수소 결합과 벤젠 고리와 규산염 층 산소 원자들의 약한 정전기적 힘에 의해 흡착되고 있음을 보여준다. 본 연구는 점토광물과 유기물에 대한 양자 화학 계산에서 클러스터 크기와 결정학적 위상이 고려되어야 함을 제시한다.