나트륨 클러스터 Naₙ (n 39‑55)의 전자구조·평형구조·온도 의존성 연구

초록

밀도범함수 이론을 이용해 n=39∼55 범위의 나트륨 클러스터 Naₙ의 구조적 변화를 조사하였다. 초기 구간( n=39‑43)에서는 대칭적인 부분 이코사헤드럴(두 겹) 구조가 나타나며, 이후 n=44‑52 구간에서는 이코사헤드럴 핵이 사라진 무질서한 형태가 지배한다. n>52에서는 세 겹 이코사헤드럴 구조가 급격히 재등장한다. 또한, 밀도범함수 분자동역학을 통해 대표적인 크기 n=43, 45, 48, 52에 대한 비열 곡선을 계산했으며, 기존에 보고된 n=40, 50, 55에 대한 열역학 데이터와 비교하였다. 전체 시리즈에 걸친 열용량 곡선은 각 클러스터의 기하학적 변천과 강하게 상관관계를 보이며, 나트륨 클러스터의 열용량이 크기에 따라 민감하게 달라짐을 확인하였다.

상세 요약

본 연구는 나트륨 금속 클러스터가 원자 수에 따라 어떻게 구조적 전이를 겪는지를 체계적으로 밝힌 점에서 의미가 크다. 먼저, n=39‑43 구간에서 관찰된 ‘부분 이코사헤드럴(두 겹)’ 구조는 전통적인 금속 클러스터 성장 모델과 일치한다. 이 구간에서는 중심에 13개의 원자가 이루는 icosahedral 핵을 두고, 그 주위에 추가 원자들이 규칙적으로 배열되어 두 번째 껍질을 형성한다. 이러한 대칭성은 전자구조 계산에서 높은 HOMO‑LUMO 간격과 안정적인 전자 구름을 제공함으로써 클러스터의 총 에너지를 최소화한다는 점에서 물리적으로 타당하다.

반면 n=44‑52 구간에서는 이코사헤드럴 핵이 사라지고 ‘무질서’한 구조가 지배한다는 점이 눈에 띈다. 이는 원자 수가 증가함에 따라 기존의 두 겹 구조를 유지하기 위한 기하학적 제약이 커져, 스트레스를 최소화하기 위해 비대칭적인 배열을 채택하게 되는 현상으로 해석할 수 있다. DFT‑MD 시뮬레이션 결과, 이 구간의 클러스터는 다수의 저에너지 메타스테이블 구조가 존재하며, 에너지 풍경이 평탄해지는 경향을 보인다. 이러한 구조적 다중성은 열용량 곡선에서 넓은 ‘용융’ 구간과 낮은 피크 강도를 초래한다.

n>52에서 급격히 나타나는 ‘세 겹 이코사헤드럴’ 구조는 다시 대칭성을 회복하는 전형적인 ‘재결정’ 현상이다. 추가 원자들이 기존 두 겹 구조에 균일하게 배치되면서 새로운 외부 껍질을 형성하고, 이는 전자구조적으로도 안정적인 전자 구름을 재구성한다. 특히 n=55와 같은 완전한 세 겹 이코사헤드럴은 높은 전자 구속 에너지와 뚜렷한 HOMO‑LUMO 간격을 보여, 열용량 곡선에서 뚜렷한 피크와 급격한 상승을 야기한다.

열역학적 측면에서 저자들은 DFT‑MD 기반의 특정 열용량(Cv) 곡선을 n=43, 45, 48, 52에 대해 계산하고, 기존 n=40, 50, 55 데이터와 비교하였다. 결과는 구조 전이와 열용량 피크 형태 사이에 강한 상관관계를 보여준다. 예를 들어, 부분 이코사헤드럴(두 겹) 구조를 가진 n=40, 43은 비교적 뾰족하고 높은 피크를 보이며, 이는 ‘고체‑액체 전이’가 급격히 일어남을 의미한다. 반면 무질서 클러스터인 n=45, 48, 52는 피크가 넓고 낮으며, 전이 온도가 분산되어 있음을 나타낸다. 최종적으로 세 겹 이코사헤드럴인 n=55는 다시 뚜렷한 피크를 보이며, 구조적 대칭성이 회복될 때 열용량도 급격히 변한다는 점을 확인한다.

이러한 결과는 금속 클러스터의 ‘크기 의존적 열역학’이 단순히 원자 수에 비례하는 것이 아니라, 내부 기하학적 재배열에 의해 크게 좌우된다는 중요한 교훈을 제공한다. 향후 촉매, 나노전자소자, 그리고 양자점 설계에 있어 클러스터의 정확한 크기와 구조를 제어하는 것이 열적 안정성 및 전자적 특성을 최적화하는 핵심 전략이 될 것으로 기대된다.