Partially reactive force field for the UiO-66 metal-organic framework

본 논문은 UiO-66 금속 - 유기 골격체의 구조적, 기계적 및 결함 특성을 정확하게 모델링하고 그 용매열 합성 및 자기 조립 메커니즘에 대한 분자 동역학 시뮬레이션을 가능하게 하는 새로운 부분 반응성 힘장인 nb-UiO-FF 를 소개한다.

핵심

UiO-66와 같은 금속 - 유기 골격체 (MOF) 를 거대하고 미세한 3 차원 퍼즐로 상상해 보세요. 이는 모서리 역할을 하는 금속 '노드'(예: 지르코늄 클러스터) 와 연결봉 역할을 하는 유기 '링크어'(예: 벤젠 고리) 로 구성되어 있습니다. 과학자들은 이러한 퍼즐을 매우 튼튼하며, 기체를 포획하거나, 약물을 전달하거나, 화학 반응을 가속화하도록 조정할 수 있기 때문에 좋아합니다.

그러나 문제가 하나 있습니다. 바로 이러한 퍼즐이 처음에 어떻게 조립되는지 정확히 파악하는 것이 빛의 속도로 재생되는 영화를 보는 것과 같다는 점입니다. 금속과 링크어를 함께 묶는 화학 결합은 컴퓨터에서 시뮬레이션하기가 까다롭습니다. 대부분의 컴퓨터 모델은 이러한 연결을 영구적인 접착제로 간주합니다. 즉, 접착제가 도포되는 과정, 조각들이 조립되는 모습, 혹은 조각이 하나라도 빠졌을 때의 상황을 보여줄 수 없습니다.

해결책: "스마트 접착제" 힘장 (Force Field)
이 논문에서 저자들은 nb-UiO-FF라는 새로운 컴퓨터 도구를 소개합니다. 이는 퍼즐 조각들이 '부분적으로 반응성'을 가질 수 있도록 하는 시뮬레이션 게임의 새로운 규칙 세트라고 생각하시면 됩니다.

다음은 몇 가지 간단한 비유를 사용하여 그들이 이를 어떻게 작동시켰는지 설명한 것입니다:

• '더미' 원자 (보이지 않는 손):
  실제 세계에서는 지르코늄 금속 노드가 복잡한 전하를 가지고 있어 링크어를 특정 방향으로 끌어당깁니다. 표준 컴퓨터 모델은 이를 messy 해지지 않고 모방하는 데 어려움을 겪습니다. 저자들은 금속 노드에 보이지 않는 '더미' 원자 (작은 자석 더미와 같은) 를 부착함으로써 이를 해결했습니다. 이러한 더미는 링크어를 올바른 모양과 방향으로 고정하는 보이지 않는 손처럼 작용하여, 복잡하고 무거운 계산 없이도 퍼즐이 올바르게 조립되도록 합니다.

• '모스 퍼텐셜' (신축성 있는 스프링):
  일반적으로 컴퓨터 모델은 결합을 뻣뻣한 막대기로 취급합니다. 당기면 즉시 끊어집니다. 저자들은 이러한 뻣뻣한 막대기를 모스 퍼텐셜로 대체했는데, 이는 더 신축성 있는 스프링처럼 작동합니다. 이를 통해 시뮬레이션은 금속과 링크어가 늘어나고, 흔들리며, 역동적으로 조립되거나 분리되는 모습을 보여줄 수 있습니다. 이는 물질의 '탄생'을 관찰하는 데 결정적입니다.

테스트 내용
저자들은 단순히 도구를 개발한 것에 그치지 않고, 신뢰성을 확보하기 위해 엄격한 스트레스 테스트를 거쳤습니다:

1. 완벽한 퍼즐: 도구가 완벽한 UiO-66 결정의 정확한 모양을 재현할 수 있는지 확인했습니다. 결과는 실제 측정값과 거의 완벽하게 일치했습니다 (수십 분의 일 퍼센트 이내).
2. 젖은 퍼즐: 도구를 결정 제조에 사용되는 두 가지 액체 (DMF 및 에탄올) 에 적신 결정으로 테스트했습니다. 모델은 결정이 젖었을 때에도 강하게 유지되며 부서지지 않음을 보여주었습니다.
3. 깨진 퍼즐: 실제 세계의 결정에는 종종 결손 (결함) 이 있습니다. 저자들은 시뮬레이션에서 의도적으로 링크어나 전체 노드를 제거했습니다. 도구는 이러한 구멍이 있더라도 실제 물질이 그러하듯 결정이 여전히 형태를 유지할 수 있음을 성공적으로 보여주었습니다.
4. 튕기는 퍼즐: 결정이 변형되기 전까지 얼마나 강하게 누를 수 있는지 테스트했습니다. 결과는 고급 물리 계산과 일치하여, 이 모델이 물질의 강도를 이해하고 있음을 입증했습니다.
5. 사촌 퍼즐: 약간 더 큰 버전의 퍼즐 (UiO-67) 에 도구를 적용해 보았으며, 그곳에서도 작동하여 규칙이 유연함을 입증했습니다.

마술이 일어나는 모습 관찰
이 논문에서 가장 흥미로운 부분은 이 새로운 도구를 사용하여 자기 조립 과정을 관찰한 것입니다. 모든 퍼즐 조각 (금속 노드와 링크어) 과 액체 용매를 상자에 넣고 '재생' 버튼을 누르는 상황을 상상해 보세요.

시뮬레이션은 조각들이 떠다니다가 서서히 서로 달라붙기 시작하는 모습을 보여주었습니다.

• 금속 노드와 링크어가 서로를 찾아 초기 빌딩 블록을 형성하는 것을 목격했습니다.
• 때로는 조각들이 '잘못된' 위치에 걸려 멈추는 (운동학적 함정) 것을 발견했는데, 이는 퍼즐 조각이 헐겁게 맞지만 정확히 맞지 않는 경우와 같습니다.
• 과정이 느리다는 것을 관찰했습니다. 조각들이 무겁고 둔하게 움직이기 때문에, 시뮬레이션한 시간 동안 퍼즐이 완전히 조립되지는 않았습니다.

결론
이 논문은 UiO-66 물질에 대한 '스마트 현미경'처럼 작동하는 새로운 고정밀 컴퓨터 모델을 제시합니다. 이 모델은 물질의 구조, 강도, 그리고 결함을 처리하는 능력을 시뮬레이션할 수 있습니다. 가장 중요한 점은, 이 도구가 물질이 처음부터 스스로 조립되는 역동적인 과정을 현실적으로 시뮬레이션할 수 있는 최초의 도구라는 것입니다. 이를 통해 과학자들은 이러한 놀라운 물질들이 어떻게 태어나고 어떻게 그 결함을 제어할 수 있는지 이해하는 데 도움을 받을 수 있습니다.

기술 요약

기술 요약: UiO-66 금속 - 유기 골격체를 위한 부분적 반응성 힘장

문제 제기
금속 - 유기 골격체 (MOF), 특히 UiO-66 은 구조적 조절 가능성과 안정성으로 인해 광범위하게 연구되고 있습니다. 그러나 반응성 과정을 시뮬레이션할 수 있는 계산 모델의 부재로 인해 합성 메커니즘과 결함 형성을 이해하는 것은 여전히 어렵습니다. 기존 힘장 (예: UFF4MOF, BTW-FF, MOF-FF) 은 일반적으로 사전 정의된 연결성에 의존하여 결합 끊김과 형성을 모델링할 수 없으며, 이는 MOF 자기 조립 및 결함 공학 연구에 필수적입니다. Reax-FF 와 같은 반응성 힘장과 머신러닝 포텐셜이 존재하지만, MOF 에 대한 정확도, 계산 비용, 또는 용매와 조절제를 포함하도록 확장하는 어려움과 관련된 과제를 안고 있습니다. 또한 UiO-66 을 위한 특정 모델들은 종종 수산화된 형태를 정확하게 표현하는 데 어려움을 겪거나, 계산 비용이 많이 드는 더미 원자 구성 (예: 노드당 48 개의 더미 원자) 을 필요로 합니다.

방법론
저자들은 UiO-66 골격체와 그 등형체 (isoreticular) 인 UiO-67 을 모델링하도록 설계된 부분적 반응성 힘장인 nb-UiO-FF를 개발했습니다. 이 방법론은 이전에 확립된 ZIF-8 접근법을 기반으로 하며, 비결합 상호작용과 양이온성 더미 원자 (CDA) 의 조합을 활용합니다.

• 힘장 구성: 금속 노드와 리간드 사이의 경직된 결합 포텐셜을 비결합 모스 (Morse) 포텐셜로 대체하여 동적인 결합 끊김과 형성을 가능하게 합니다.
• 더미 원자 (CDA): 결합 항 없이 Zr 원자의 이방성 전하 분포를 재현하고 올바른 배위 기하학을 유지하기 위해, 각 Zr 중심당 4 개의 더미 원자 (D) 가 도입됩니다. 이들은 강한 조화 결합과 각도를 통해 Zr 원자에 연결됩니다. Zr 원자의 +1.968 전하는 이 네 개의 더미 원자 (+0.492 씩) 에 완전히 분배되어 Zr 원자가 전기적으로 중성이 됩니다.
• 매개변수화: 힘장 매개변수는 Yang 등에게서 유래된 비반응성 힘장을 기반으로 하여, 수산화된 Zr6O4(OH)4 노드에 맞게 조정되었습니다. 금속 클러스터의 결합 매개변수는 보고된 모델에서 가져왔으며, 평형 결합 길이와 각도는 Valenzano 등의 DFT 결과에서 출처를 얻었습니다. Zr–리간드 상호작용에 대한 모스 포텐셜 매개변수 ($D_0$ 및 $\alpha$) 는 DFT 및 실험 참조 대비 수밀도와 국소 결합 길이 (특히 Zr–리간드 방사 분포 함수) 의 정확도를 균형 있게 맞추기 위해 매개변수 공간을 체계적으로 탐색하여 최적화되었습니다.
• 시뮬레이션 프로토콜: 모든 시뮬레이션은 LAMMPS 를 사용하여 수행되었습니다. 검증에는 구조 분석, 용매 (DMF 및 에탄올) 내 안정성 테스트, 그리고 기계적 특성 계산이 포함되었습니다. 자기 조립 연구는 16 개의 Zr 클러스터, 96 개의 BDC 리간드, 2136 개의 DMF 분자로 구성된 시스템을 사용하여 NPT 앙상블에서 6 ns 동안 시뮬레이션되었습니다.

주요 결과

• 구조적 정확도: nb-UiO-FF 는 순수한 UiO-66 의 실험적 격자 상수를 약 0.046% 의 편차 (20.71 Å 대 20.7004 Å) 로 재현합니다. 주요 원자 쌍 (Zr–$\mu_3$O, Zr–$\mu_3$OH, Zr–Zr, Zr–리간드) 에 대한 방사 분포 함수 (RDF) 는 DFT 및 실험 데이터와 합리적인 일치를 보이며, 편차는 일반적으로 10% 이내입니다.
• 용매 및 결함 안정성: 이 힘장은 DMF 및 에탄올 용매의 존재 하에서, 그리고 결손 리간드와 결손 노드가 포함된 결함 구조에서도 골격체의 무결성을 유지합니다. 결함 모델에서는 하전 불균형 중심을 염화 이온 또는 단치 BDC 조각으로 캡핑하여 전하 중성을 유지했습니다.
• 기계적 특성: 계산된 탄성 상수 ($C_{11}$, $C_{12}$, $C_{44}$) 및 유도된 기계적 특성 (체적 탄성률: 44 GPa, 영률: 40 GPa) 은 DFT 벤치마크와 잘 일치하지만, Zr 과 더미 원자 사이의 강한 결합 상호작용으로 인한 강성으로 인해 약간 더 높은 값을 보입니다.
• 전이성: 이 모델은 등형체인 UiO-67 골격체의 구조를 성공적으로 재현하여 서로 다른 리간드 길이 전반에 걸친 전이성을 입증했습니다.
• 자기 조립 역학: 동적 시뮬레이션에서 이 힘장은 금속 클러스터와 유기 리간드 사이의 결합 형성의 점진적 과정을 성공적으로 포착했습니다. 연구는 초기 2 차 구조 단위 (SBU) 와 비이상적 결합 구성 (동역학적으로 갇힌 상태) 을 포함한 일시적 구조 모티프를 확인했습니다. 그러나 6 ns 시뮬레이션 창 내에서 배위수의 완전한 수렴과 대형 이온성 응집체의 형성은 달성되지 않았으며, 이는 더 긴 시간 규모를 위한 향상된 샘플링 방법이 필요함을 시사합니다.

의의 및 주장
이 논문은 nb-UiO-FF 가 정적 구조 모델링과 반응성 역학 사이의 간극을 연결하며, UiO-66 과 UiO-67 에 대한 일관되고 물리적으로 신뢰할 수 있는 설명을 제공한다고 주장합니다. 이전 모델에 비해 더미 원자 수를 줄임 (Zr 당 4 개) 으로써, 저자들은 모델 효율성과 반응성 과정 시뮬레이션 능력 사이의 균형을 달성했습니다. 이 힘장은 전통적인 고정 연결성 힘장이 접근할 수 없는 자기 조립 메커니즘과 결함 형성 연구를 가능하게 합니다. 저자들은 이 작업을 MOF 합성의 초기 단계 조사와 점 결함 공학의 기초로 위치시키면서, 동역학적 함정을 완전히 극복하고 완전한 자기 조립을 모델링하기 위해서는 메타동역학과 같은 향상된 샘플링 계획을 통한 추가 작업이 필요함을 인정합니다.