Temperature-Induced Reorganization of Supported Zn$_3$ Clusters on Cu(111): From Minimum-Energy Structures to Finite-Temperature Ensembles

기술 요약: Cu(111) 상에 지지된 $\text{Zn}_3$ 클러스터의 온도 유도 재구성

문제 제기
불균일 촉매, 특히 메탄올 합성에 사용되는 산업용 구리/산화아연/알루미나(CZA) 촉매에서, 원자 수준의 구조와 안정성은 여전히 논쟁의 대상이다. Cu 계면에서의 작은 Zn 기반 종들이 활성 부위로 제안되고 있으나, 이들의 식별은 전적으로 $0\text{ K}$ 밀도 범함수 이론(DFT) 최적화에 크게 의존하고 있다. 이러한 접근 방식은 $0\text{ K}$에서의 최소 에너지 구성이 작동 조건하에서도 활성 부위를 나타낸다고 가정한다. 그러나 촉매는 실제 작동 온도에서 상당한 구조적 유연성, 원자 이동, 그리고 흡착제와의 동적 상호작용을 나타내는 경우가 많다. 결과적으로, 열역학적으로 유의미한 활성 부위는 단일한 정적 구조가 아니라 작동 온도에서 형성되는 다양한 구성들의 앙상블(ensemble)일 수 있다. 제일원리 방법을 사용하여 유한 온도에서의 자유 에너지 지형을 명시적으로 샘플링하는 것은 필요한 시간 및 길이 척도로 인해 계산적으로 매우 부담스럽다.

방법론
저자들은 머신러닝 원자간 포텐셜(MLIP), 전역 최적화(global optimization), 분자 동역학(MD), 그리고 향상된 샘플링 자유 에너지 계산을 결합한 다중 스케일 계산 프레임워크를 채택하였다.

1. MLIP 개발: 순수 금속(Cu, Zn), 합금, 산화물 및 지지된 클러스터를 아우르는 44,269개의 구조 데이터셋을 바탕으로 딥 포텐셜(Deep Potential, DP) 모델을 훈련시켰다. 훈련 데이터는 PBE+D3 수준의 DFT를 통해 생성되었다. 결과적인 모델은 에너지 측면에서 0.005 eV/atom, 힘(force) 측면에서 0.056 eV/Å의 제곱평균제곱근 오차(RMSE)를 달축성하여, DFT에 근접한 정확도를 유지하면서도 시뮬레이션 속도를 크게 가속화하였다.
2. 시스템 모델: 연구는 메탄올 합성과 관련된 두 가지 클러스터 모델인 수산화된 $\text{Zn}_3(\text{OH})_3$와 포메이트를 포함한 $\text{Zn}_3(\text{OH})_2\text{CHOO}$를 조사하였다. 이들은 다양한 정도의 Zn 합금을 모사하기 위해 세 가지 서로 다른 Cu(111) 기반 표면에 지지되었다:
   • 순수 Cu(111) (0% Zn).
   • 단일 Zn 사이트 표면 (희박한 합금, ~1.56% Zn).
   • 합금화된 $\text{Cu}_3\text{Zn}$ 표면 (25% Zn).
3. 전역 최적화: DP 모델과 함께 유전 알고리즘(genetic algorithm)을 사용하여 각 표면 위 클러스터의 저에너지 삼각형 및 선형 구성을 식별하였으며, 이후 DFT로 검증하였다.
4. 유한 온도 샘플링:
   • 비편향 MD (Unbiased MD): 300 K와 450 K에서 20 ns 시뮬레이션을 수행하여 자발적인 구조 전이와 클러스터의 이동성을 관찰하였다.
   • 향상된 샘플링 (Enhanced Sampling): 온더플라이 확률 향상 샘플링(On-the-fly Probability Enhanced Sampling, OPES) 방법을 사용하여 150, 300, 450 K에서의 자유 에너지 지형을 재구성하였다. 두 가지 집단 변수(collective variables, CVs)가 사용되었다: 삼각형($\approx 60^\circ$)과 선형($\approx 170^\circ$) 기하 구조를 구분하기 위한 최대 $\text{Zn-Zn-Zn}$ 각도($\theta_{max}$), 그리고 전이를 가속화하기 위한 거리 기반 기술자($\tilde{d}_{max}$).
5. 이동성 분석: 클러스터 질량 중심의 평균 제곱 변위(MSD)와 사이트 상관 함수를 계산하여 공간적 안정성과 소결(sintering)에 대한 저항성을 평가하였다.

주요 결과

• 온도 유도 구성 재구성:

  • $0\text{ K}$에서, 순수 Cu(111) 위의 $\text{Zn}_3(\text{OH})_3$는 일반적으로 조밀한 삼각형 구성이 최소 에너지 구조이다.
  • 유한 온도 자유 에너지 계산 결과, 온도가 상승함에 따라 확장된 선형 구성으로의 뚜렷한 이동이 나타났다. 450 K에서는 선형 구성이 삼각형 구성과 열역학적으로 경쟁할 만큼 안정해진다.
  • 이러한 전이는 주로 엔트로피적 안정화에 의해 구동된다. 온도 변화에 따라 구성 간의 퍼텐셜 에너지 차이는 비교적 일정하게 유지되지만, 선형 구성은 더 높은 엔트로피를 가지며, 이는 온도가 높아짐에 따라 엔탈피적 불리함을 점진적으로 상쇄시킨다.

• 표면 합금화의 영향:

  • 순수 및 희박 표면: 순수 Cu(111) 및 단일 Zn 사이트 표면에서, 온도가 증가함에 따라 선형 앙상블이 점진적으로 안정화되어 450 K에서는 삼각형 구조와 선형 구조의 인구(population)가 비슷해진다.
  • 합금화된 표면: $\text{Cu}_3\text{Zn}$ 표면에서는 에너지 순위가 $0\text{ K}$에서도 역전되어, 선형 구성이 0.19 eV만큼 더 안정적이다. 광범위한 Zn 합금화는 전체 온도 범위에서 선형 구성을 안정화하고 경직시킨다.
  • 흡착제의 역할: 포메이트($\text{Zn}_3(\text{OH})_2\text{CHOO}$)의 존재는 선형 구성을 향한 에너지 균형을 더욱 심화시키지만, 온도 및 합금화에 따른 동일한 정성적 경향은 유지한다.

• 공간적 안정성 및 이동성:

  • 순수 Cu(111) 위에서 지지된 $\text{Zn}_3$ 클러스터는 상당한 횡방향 이동성을 보였으며, 이는 이동 및 응집(소결)에 대한 높은 성향을 나타낸다.
  • Zn 합금화는 확산을 강력하게 억제한다. 단일 Zn 사이트 표면의 경우, 구성적 평형은 순수 표면과 유사하지만, 고립된 표면 Zn 원자의 존재가 안정적인 $\text{Zn}_c\text{-adsorbate-}\text{Zn}_s$ 브릿징 모티프를 형성한다. 이러한 모티프는 앵커링 센터(anchoring center) 역할을 하여 클러스터의 체류 시간을 크게 연장시키고, $\text{Zn}_c\text{-Cu}$ 상호작용에 비해 이동을 억제한다.
  • 고도로 합금화된 $\text{Cu}_3\text{Zn}$ 표면에서는 확산이 거의 발생하지 않으며, 이는 고립된 계면 Zn 종을 소결으로부터 효과적으로 안정화한다.

의의 및 주장
본 논문은 지지된 Zn 기반 클러스터의 열역학적으로 유의미한 구조가 정적인 $0\text{ K}$ 최적화에 의해 예측된 것과 근본적으로 다르다는 것을 보여준다. 활성 부위 앙상블은 엔탈피와 엔트로피 효과 사이의 경쟁에 의해 결정되며, 엔트로피는 작동 조건하에서 조밀한 구조에서 확장된 기하 구조로의 변화를 유도할 수 있다.

저자들은 활성 부위를 단지 최소 에너지 구조에 기반하여 식별하는 것만으로는 불충분하다고 주장한다. 대신, 활성 부위는 온도, 흡착제 피복도 및 국부 조성에 따라 인구와 표면 이동성이 진화하는 동적인 앙상블로 간주되어야 한다. 본 연구는 정적 최적화가 삼각형 클러스터가 지배적이라고 예측할 수 있지만, 유한 온도 샘플링은 특히 작동 온도에서 선형 구조의 상당한 인구를 드러낸다는 점을 강조한다. 또한, 본 연구는 표면 조성이 구성적 평형뿐만 아니라 클러스터의 공간적 안정성까지 제어하며, Zn 합금화가 앵커링 모티프를 생성함으로써 소결을 방지하는 데 결정적인 역할을 한다는 점을 밝힌다.

이러한 결과는 실험적 관측(예: 삼각형 구조를 보여주는 저온 STM 데이터)을 해석하기 위한 틀을 제공하며, 운전 중인 촉매의 복잡하고 동적인 구조에 근접하기 위해 운전 온도 근처에서 나타나는 확장된 동적 앙상블을 타겟팅해야 함을 시사한다. 본 연구는 MLIP 가속 자유 에너지 샘플링이 이상적인 정적 모델과 실제 작동 중인 촉매의 복잡하고 동적인 구조 사이의 간극을 메우는 실질적인 경로임을 입증한다.