Negative and Zero Linear Compressibility in MCN (M = Ag, Au, Cu): A First-Principles Study

기술 요약: MCN(M = Ag, Au, Cu)에서의 음의 선형 압축률 및 제로 선형 압축률

문제 제기
음의 선형 압축률(NLC)과 제로 선형 압축률(ZLC)은 정수압 하에서 결정축이 팽창하거나 일정하게 유지되는 희귀한 기계적 현상이다. 이러한 현상은 프레임워크 물질, 무기 산화물, 분자 결정에서 알려져 있으나, 후보 물질에서 넓은 압력 범위에 걸쳐 강하고 안정적인 NLC/ZLC를 달ieving하는 것은 여전히 도전적인 과제이다. "와인 랙(wine-rack)" 힌징, 다면체 기울임, 층 슬라이딩과 같은 기존 메커니즘은 잘 문서화되어 있으나 제한적이다. 금속 시안족 계열(MCN, M = Ag, Au, Cu)은 음의 열팽창과 특정 상(예: CuCN $C222_1$)에서의 거대한 NLC를 포함한 이상 거동을 보여왔으나, AgCN, AuCN, CuCN의 $R\bar{3}m$ 및 $P6mm$ 상에서 ZLC 또는 NLC가 나타날 가능성은 아직 충분히 탐구되지 않았다.

연구 방법론
본 연구는 VASP(Vienna Ab-initio Simulation Package)를 이용한 밀도 범함수 이론(DFT) 기반의 제일원리 계산을 채택하였다.

• 조사된 구조: 저자들은 실험적으로 보고된 상(AgCN $R\bar{3}m$, AuCN $P6mm$, CuCN $R\bar{3}m$)과 이전에 보고되지 않은 세 가지 상(AgCN $P6mm$, CuCN $P6mm$, AuCN $R\bar{3}m$)을 포함하여 총 6개의 구조를 분석하였다. 보고되지 않은 상들은 알려진 구조의 금속 사이트를 치환하여 생성되었다.
• 범함수: 신뢰성을 결정하기 위해 PBE, PBE-D3, PBEsol, PBEsol-D3, r2SCAN, r2SCAN-D3의 6가지 교환-상관 범함수를 테스트하였다. 실험 데이터와의 격자 상수 비교를 바탕으로, 상세한 물성 분석을 위해 PBE-D3, PBEsol, r2SCAN이 선택되었다.
• 계산: 전체 구조 완화(relaxation)를 수행하였다. 탄성 상수($C_{ij}$)는 유한 차분법을 통해 계산되었으며, 이를 통해 선형 압축률($\beta$), 영률($E$), 부피 탄성률($K$)을 도출하였다. 구조적 안정성과 반응 진화를 평가하기 위해 최대 13 GPa의 정수압을 가하였다. 결합 강도와 원자 진동을 조사하기 위해 포논 분산 곡선을 계산하였다.

주요 기여 및 결과

1. 새로운 상의 예측: 본 연구는 실험적으로 보고되지 않았던 에너지적으로 경쟁력 있는 상들, 즉 AgCN과 CuCN의 $P6mm$ 상, 그리고 AuCN의 $R\bar{3}m$ 상의 존재를 예측한다. 제안된 상들과 기존 실험적 대응 상 사이의 에너지 차이는 매우 작으며(종종 < 50 meV/f.u.), 이는 이들이 실행 가능한 후보임을 시사한다.
2. 극한의 탄성 이방성 및 이상 압축률: 조사된 6개 MCN 상 모두 $c$-방향 탄성 상수($C_{33}$)가 $ab$-평면($C_{11}$)보다 한 자릿수 더 큰 극심한 탄성 이방성을 보인다.
   • ZLC/NLC 예측: 이러한 이방성은 $c$-축 방향의 거의 제로 또는 음의 선형 압축률로 이어진다. 구체적으로, AgCN $R\bar{3}m$ 상은 압력 하에서 NLC(c-축 팽창)를 나타낸다. 다른 상들은 근접 ZLC 거동을 보인다.
   • 압력 안정성: 이러한 이상적인 반응은 구조적 불안정성 없이 넓은 압력 범위(최대 13 GPa)에 걸쳐 지속된다.
3. 새로운 메커니즘 식별: 본 논문은 전통적인 "와인 랙" 또는 힌징 모티프와 구별되는 메커니즘을 식별하였다. MCN 구조는 강직한 1차원 MCN 로드(rod)로 구성된 "대나무 숲(bamboo forest)" 기하 구조를 형성한다.
   • 강직한 로드: 로드는 가벼운 원자들(M-C-N) 사이의 매우 강한 결합에 의해 유지되며, 이는 높은 주파수의 포논 모드(최대 68 THz)에 의해 입증된다.
   • 희소한 패킹: 로드들은 수직 평면에서 약한 힘을 통해 상호작용하며, 이는 매우 낮은 주파수의 포논 모드(~1.2 THz)에 의해 입증된다.
   • 반응: 정수압 하에서 부피 감소는 로드 자체를 압축하기보다는 로드들이 서로 더 가까워지는 방식(ab-평면에서의 패킹 밀도 증가)에 의해 거의 전적으로 수용된다. 특정 사례(AgCN $R\$\bar{3}m$)에서는, 이러한 재배열이 결합 각도를 유지하기 위해 강직한 로드를 그 길이 방향($c$-축)으로 약간 팽창하게 만들어 NLC를 유발한다.

의의
저자들은 본 연구가 ZLC 또는 NLC를 나타내는 것으로 알려진 상대적으로 작은 물질군을 확장한다고 주장한다. 강직한 1D 체인의 희소한 패킹이 이상적인 압축률을 허용하는 "대나무 숲" 메커니즘을 식별함으로써, 본 연구는 이전에 보고된 메커니즘과는 구별되는 이러한 현상의 새로운 미시적 기원을 제공한다. 이 결과는 구조적 이방성과 기계적 반응 사이의 관계에 대한 더 깊은 통찰을 제공하며, 정밀 기계 부품 및 압력 센서를 위한 재료 설계에 잠재적인 지침을 제공할 수 있으나, 본 논문은 특정 장치 응용보다는 근본적인 이해에 초점을 맞추고 있다. 또한, 본 연구는 복잡한 시안족 시스템을 정확하게 모델링하는 데 있어 특정 DFT 범함수(r2SCAN 및 PBE-D3)의 유용성을 강조한다.