Controlling HER activity and stability of $γ$- and 6,6,12-Graphyne through engineered B-N doping: DFT and Reactive MD simulations

기술 요약: 설계된 B-N 도핑을 통한 $\gamma$- 및 6,6,12-그래파인의 HER 활성 및 안정성 제어

문제 정의
수소 발생 반응(HER)은 지속 가능한 수소 생산에 필수적이지만, 실제 응용을 위해서는 높은 활성과 장기적인 안정성 사이의 균형을 맞추는 촉매가 필요하다. 귀금속이 표준이지만, 지구가 풍부한 금속 불포함 대안에 대한 필요성이 존재한다. 2차원 탄소 동소체, 특히 그래파인(sp/sp$^2$ 혼성 탄소 시트)은 조절 가능한 전자 구조를 제공한다. 그러나 두 가지 특정 과제가 그래파인 기반 HER 촉매에 남아 있다: (1) 붕소-질소(B-N) 공동 도핑의 상대적 기하 구조(meta, ortho, para)가 결함 열역학 및 페르미 준위($E_F$) 근처의 전자 상태를 어떻게 동시에 제어하는지 이해하는 것, 그리고 (2) 0 K에서 수소 흡착 자유 에너지($\Delta G_{ads}$)를 최적화하도록 예측된 모티프가 유한한 온도에서 수소화에 의한 결합 절단에 대해 운동학적으로 견고한지를 결정하는 것이다.

방법론
본 연구는 순수, B-도핑, N-도핑 및 B-N 공동 도핑된 $\gamma$-그래파인과 6,6,12-그래파인을 평가하기 위해 밀도 범함수 이론(DFT)과 반응 분자 역학(ReaxFF) 시뮬레이션을 결합한 다중 스케일 계산 접근법을 채택한다.

• DFT 계산: GGA-PBE 범함수와 DZP 기저 집합을 사용하는 SIESTA 코드를 사용하여 수행되었다. 본 연구는 전자 밴드 구조, 결함 형성 에너지($E_{form}$), 계산적 수소 전극 프레임워크 내에서의 수소 흡착 깁스 자유 에너지($\Delta G_{ads}$)를 계산하였다. 스핀 편극 계산을 통해 해당 구조들이 비자성임을 확인하였다.
• 반응 분자 역학(MD) 시뮬레이션: 300 K에서 그래파인 시트와 원자 수소 간의 상호작용을 시뮬레이션하기 위해 ReaxFF 힘의 장(force field)을 사용하는 LAMMPS를 사용하여 수행되었다. 이 시뮬레이션은 결합의 끊어짐과 형성을 명시적으로 모델링하며, 유한 온도에서의 구조적 안정성과 수소 흡수 역학을 평가하기 위해 NVT 앙상블에서 500 ps 동안 실행되었다.

주요 기여 및 결과

1. 전자 구조 및 도핑 효과:

   • 순수 $\gamma$-그래파인은 반도체이며, 6,6,12-그래파인은 이방성 디락 유사 준금속 분산을 보인다.
   • 단일 B 또는 N 치환은 단순한 리지드 밴드 이동(rigid band shift)이 아닌 도펀트-$\pi$ 혼성을 통해 $E_F$ 근처의 상태를 재구성한다. B는 수용체(acceptor)로, N은 공여체(donor)로 작용한다.
   • B-N 공동 도핑은 기하 구조 의존적인 섭동을 도입한다. B-N ortho 구성이 공여체-수용체 전하 보상으로 인해 가장 열역학적으로 유리한 것으로 확인되었으며, 그 뒤를 meta와 para가 잇는다.

2. 열역학적 활성 ($\Delta G_{ads}$):

   • 순수 격자는 수소 결합력이 약하여 HER에 최적화되지 못한다.
   • 도핑은 수소 흡착을 위한 "핫스팟"을 생성한다. 결정적으로, 가장 활발한 사이트는 도펀트 원자 자체가 아니라 인접한 탄소 사이트, 특히 sp 혼성 아세틸렌 사슬 내의 탄소들이다.
   • N-도핑 단독으로는 N의 높은 전기 음성도로 인해 직접적인 H 결합에 대해 열역학적으로 불리하다. 대신, N은 주변 골격을 편극시킨다.
   • B-N 공동 도핑은 $\Delta G_{ads}$를 사바티에 최적값(Sabatier optimum, $\approx 0$ eV)으로 조율한다. 6,6,12-그래파인의 경우, ortho 치환은 특정 sp-근접 사이트를 열중립 영역($\Delta G_{ads} \approx 0.03$ eV) 내에 직접 위치시킨다.

3. 운동학적 안정성 및 활성-안정성 트레이드오프:

   • 300 K에서의 반응 MD는 결정적인 트레이드오프를 보여준다: 전자적 환경을 강력하게 활성화하는 구성은 과잉 수소화에 의한 파괴적인 결합 절단을 유발할 수도 있다.
   • $\gamma$-그래파인: B-N ortho 구성은 구조적 퇴화 없이 제어된 수소 흡수를 유일하게 유지한다. 반면, B-N meta 및 para 치환은 결합 절단을 초래한다. 단일 B 또는 N 도핑은 구조적으로 온전하게 유지된다.
   • 6,6,12-그래파인: 이 격자는 일반적으로 과잉 수소화 및 결합 끊어짐에 더 취약하다. 순수 시트조차 초기 단계의 퇴화를 보이며, B-N 공동 도핑은 특히 meta 및 para 구성에서 심각한 네트워크 변형과 결합 절단을 유발한다.

의의 및 주장
본 논문은 B-N 기하 구조(특히 도펀트의 상대적 위치)가 전하 재분배와 구조적 안정성을 공동 제어하는 그래파인 기반 HER 촉매의 주요 설계 변수임을 입증한다.

저자들은 근사 열중립적 $\Delta G_{ads}$를 달ieving하는 것이 필요한 열역학적 기준이지만, 그것만으로는 불충분하다고 주장한다. 실용적인 촉매는 수소화에 의한 퇴화에 대해 운동학적 견고성 또한 만족해야 한다. 본 연구는 적절한 흡착 에너지와 유한 온도 안정성 사이의 균형을 맞추는 유망한 후보로 B-N ortho-도핑된 $\gamma$-그래파인을 식별하였다. 반대로, 결과는 6,6,12-그래파인이 공동 도핑될 때 특히 수소가 풍부한 조건에서 본질적으로 덜 안정적임을 시사한다. 이 연구는 촉매적으로 유의미한 모티프가 헤테로 원자 자체가 아니라 "BN-편극된 탄소 환경"(특히 인접한 sp-탄소)임을 강조하며, 전자적 활성화와 구조적 무결성 사이의 상호작용이 동시에 최적화되어야 함을 보여준다.