칼슘탄산염 나노클러스터의 에너지와 구조 변천 연구

초록

본 연구는 진공 상태에서 CaCO₃ 클러스터를 12000 화학식 단위까지 구성하고, GULP 기반의 Rohl 전위와 DFT(SIESTA) 검증을 통해 에너지와 구조 변화를 조사하였다. 작은 클러스터(≤12 화학식)는 무정형 형태가 가장 안정적이며, 중간 규모(14160)에서는 부분적인 결정질 질서가 나타난다. 170개 이상에서는 초기 구조가 유지되지만 표면 재배열이 일어나며, 모든 클러스터는 대량 결정보다 약 25~30 kJ·mol⁻¹ 높은 에너지를 가진다. 이러한 결과는 고전적 핵생성 이론이 원자 수준에서의 핵 형성을 설명하기에 부족함을 시사한다.

상세 분석

이 논문은 CaCO₃ 다형체(칼사이트, 아라고나이트, 바테라이트)의 이상적인 원자 배열을 출발점으로 삼아, 진공에서 고립된 클러스터의 구조 최적화를 수행한다. 짧은 거리 상호작용은 Buckingham·Morse·스프링·삼체·평면 외력 등 12개의 파라미터로 구성된 Rohl 전위 집합을 사용했으며, 전하 분극을 묘사하기 위해 산소 이온에 쉘 모델을 적용하였다. 전위의 신뢰성을 확보하기 위해 10개 이하의 소규모 클러스터에 대해 SIESTA 기반 DFT 계산을 병행했으며, 에너지와 최종 좌표가 GULP 결과와 일치함을 확인하였다.

에너지 분석은 ‘비이완(unrelaxed)’과 ‘이완(relaxed)’ 두 단계로 나뉘어 진행되었다. 비이완 클러스터는 초기 결정 구조를 그대로 유지하므로, 표면 원자들의 배위가 완전하지 않아 에너지 감소가 급격히 나타난다. 특히 3034 화학식(칼사이트)이나 1216 화학식(아라고나이트) 구간에서 비연속적인 안정성 피크가 관찰되었다. 이완 단계에서는 구조 최적화가 진행되어 표면의 CO₃²⁻ 그룹이 평면에 접하도록 회전하고, Ca²⁺의 배위 다각형이 부분적으로 복구된다. 그 결과 에너지 감소 폭은 비이완보다 완만해지지만, 2000 화학식 클러스터조차도 대체 결정(칼사이트, 아라고나이트)보다 약 25~30 kJ·mol⁻¹ 높은 에너지를 유지한다.

클러스터 크기에 따라 세 가지 전형적인 군이 제시된다. (1) 12 이하의 초소형 클러스터는 전부 표면 원자이며, Ca 배위가 완전하지 않아 무정형 구조가 가장 안정적이다. XRD 시뮬레이션에서도 뚜렷한 회절 피크가 나타나지 않는다. (2) 14~160 사이의 중간 규모 클러스터는 내부에 칼사이트·아라고나이트 특유의 Ca–CO₃ 층이 보존되지만, 표면은 크게 재배열된다. 이때는 초기 구조에 따라 에너지 차이가 발생하고, 일부 아라고나이트 기반 클러스터가 칼사이트보다 낮은 에너지를 보인다. (3) 170 이상 대형 클러스터는 내부 구조가 거의 완전한 결정과 동일해지며, 표면 재배열만이 남는다. 이 경우 초기 구조가 칼사이트인 경우가 전반적으로 더 낮은 에너지를 갖는다.

XRD 패턴을 Debye 식으로 계산한 결과, 중·대형 클러스터는 초기 다형체의 주요 회절선이 남아 있어 실험적 검출이 가능하지만, 소형 클러스터는 회절 신호가 거의 사라진다. 이러한 구조·에너지 관계는 핵생성 과정에서 ‘전구상(precursor) 무정형’ 혹은 ‘아라고나이트’ 단계가 존재할 수 있음을 원자 수준에서 뒷받침한다. 또한, 고전적 핵생성 이론이 가정하는 일정한 표면 에너지와 구형 핵 모델이 실제 원자 배열에서는 적용되기 어렵다는 점을 강조한다.