최근 실험적 연구는 분자 빔 에피택시로 Ag(111) 위에서 원자 수준 얇은 붕소 시트—보로핀—의 성공적 성장을 입증하였다. 두 가지 뚜렷한 표면 상이 관찰되었으나, 정확한 원자 배열과 열역학적 안정성은 여전히 논쟁 중이다. 본 연구에서는 밀도 범함수 이론을 사용하여 Ag(111) 위에서 지지된 가장 그럴듯한 보로핀 다형체—버클드 삼각형, β12, χ3—를 조사하였다. 계산 결과 세 가지 독립 시트 모두 열역학적으로 불안정하고 금속성이지만, Ag(111) 기판이 이들을 안정화시킴을 밝혔다. 시뮬레이션된 주사 터널링 현미경(STM) 이미지는 실험적으로 관찰된 스트라이프 상과 균질 상을 재현하며 근본적인 격자 구조를 명확히 식별한다.
배경
그래핀 발견 이후 2차원(2D) 물질은 독특한 전자적, 기계적, 촉매적 특성으로 집중적 연구 관심을 끌어왔다. 지난 10년간 2D 붕소 시트—총칭 보로핀—에 대한 다양한 이론적 제안이 보고되었으며, 육각형 공공과 버클링 패턴의 다양한 배열을 특징으로 한다. 광범위한 제1원리 연구는 주기적 육각형 구멍이 있는 평면 삼각형 붕소 격자가 에너지적으로 유리함을 시사했으나, 2D 형태는 여전히 벌크 붕소보다 에너지가 높아 3D 결정화로 향하는 열역학적 추진력이 있음을 암시한다. 따라서 3D 핵생성을 억제하고 원자 수준 얇은 층 형성을 촉진하기 위해 충분히 "끈적한" 기판이 필요하다.
최근 실험적 돌파구로 금속 표면에서 보로핀 구현이 이루어졌다. Mannix 등은 Ag(111)에 붕소를 증착하고 고해상도 STM으로 스트라이프 상과 균질 상 두 가지 뚜렷한 표면 상을 관찰하였다. Feng 등은 유사한 상을 보고하고 균질 구조를 χ3 격자로, 스트라이프 상을 버클드 삼각형 격자 또는 직사각형 β12 격자로 해석하였다. 이러한 상이한 해석은 가능한 원자 구성과 Ag(111) 위 안정성에 대한 체계적 이론적 조사의 필요성을 강조한다.
구조 식별을 넘어 보로핀 다형체는 놀라운 특성을 보인다. 버클드 삼각형 시트는 그래핀을 초과하는 영방향 영률과 강한 이방성 격자 열전도도를 가진 높은 이방성 금속으로 예측된다. β12와 χ3 보로핀에서는 MgB2 박막에 필적할 초전도성도 예측된다. 그러나 β12와 χ3의 열역학적 안정성은 논쟁적이며, 일부 연구는 Γ 점 근처에서 허수 포논 주파수를 보고한다.
계산 방법
모든 계산은 Vienna Ab-initio Simulation Package(VASP)를 사용한 밀도 범함수 이론(DFT)으로 수행하였다. 투영 증강파 법으로 전자-이온 상호작용을 기술하였고, 일반화 기울기 근사(GGA)와 Perdew-Burke-Ernzerhof(PBE) 범함수를 교환-상관에 사용하였다. 평면파 기저 세트는 500 eV에서 절단되었으며, 브릴루앵 존은 버클드 삼각형, β12, χ3 상에 각각 25×15×1, 15×9×1, 11×11×1 k-포인트 메쉬로 샘플링하였다. z 방향 최소 20 Å 진공 간격으로 주기적 이미지 간 상호작용을 최소화하였다. 기하학 최적화는 힘을 0.02 eV/Å 미만으로 수렴시켰으며, Ag 하부 두 층은 고정하였다. 포논 분산은 PHONOPY에 구현된 유한 변위 법으로 얻었다.
STM 이미지는 Tersoff-Hamann 프레임워크 내에서 시뮬레이션되었으며, 이는 터널링 전류를 시료의 국소 상태 밀도(LDOS)와 연관시킨다. 일정 전류 모드를 사용하여 다양한 바이어스 전압에서 이미지를 생성하여 실험 STM 데이터와 직접 비교하였다.
결과 및 논의
그림 1은 고려된 세 보로핀 다형체의 최적화된 기하학을 보여준다: 버클드 삼각형(a), β12(b), χ3(c). 버클드 삼각형 시트는 b-축을 따라 단일 아웃오브플레인 버클링을 보이며, β12와 χ3는 육각형 공공의 독특한 배열을 가진 평면 구조이다. 격자 상수—버클드 삼각형 1.613 Å × 2.866 Å, β12 2.916 Å × 5.075 Å, χ3 4.448 Å—는 이전 이론 및 실험 보고와 잘 일치한다.
버클드 삼각형(a), β12(b), χ3(c) 붕소 시트의 상단 및 측면 뷰. 녹색 구는 붕소 원자. 실선으로 둘러싸인 직사각형과 마름모는 단위 셀. 문자 a와 b는 격자 파라미터 나타냄
공공 농도 η는 β12에서 1/6, χ3에서 1/5로 χ3에서 더 높은 결함 정렬도를 나타낸다. 상대적 안정성은 붕소 원자당 평균 결합 에너지 EFB(=Eboron/n)로 평가하였다. 결과(표 2)는 독립 β12 시트가 가장 안정적이고 χ3가 가장 불안정하며 원자당 0.08 eV 차이임을 보여준다. 그러나 포논 분산 계산(그림 2)은 세 구조 모두 Γ 근처에서 허수 주파수를 보여 독립적으로 열역학적으로 불안정함을 확인한다.
a 버클드 삼각형, b β12, c χ3 붕소 시트의 포논 분산. 고대칭 점은 왼쪽 모서리에 표시
전자 밴드 구조(그림 3)는 세 다형체 모두 금속성임을 확인한다. 버클드 삼각형 시트는 강한 이방성을 보이며 S-Y와 G-X에서 밴드 교차, X-S(9.63 eV)와 Y-G(4.32 eV)에서 큰 밴드갭을 가진다. 이는 비버클드 a-축으로 제한된 방향성 전도도를 시사한다.
a 버클드 삼각형, b β12, c χ3 붕소 시트의 계산된 밴드 구조. 페르미 에너지는 0으로 설정. 고대칭 점은 왼쪽 모서리에 표시
Ag(111) 위에 놓았을 때 격자 부정합은 β12에서 최소(~1%), 버클드 삼각형에서 중간(~3%)이다. χ3 시트는 기판에서 두 가지 뚜렷한 레지스트리(χ3와 χ3')를 채택하며, 격자 상수는 각각 8.67 Å와 2.89 Å × 25.02 Å이다. 보로핀 층과 Ag 표면 사이 수직 간격은 2.4~3.3 Å로, β12, χ3, χ3'의 평면 특성을 보존하는 약한 물리흡착을 나타낸다. 이 층간 거리는 Ag(111) 위 보로핀의 실험적 두께 보고와 일치한다.
시뮬레이션된 STM 이미지(그림 5)는 실험적으로 관찰된 스트라이프 상(Mannix 등)과 균질 상(Feng 등)을 재현한다. 버클드 삼각형의 스트라이프 패턴은 방추형 밝은 능선으로 나타나며, β12는 밝은 모서리로 둘러싸인 타원형 돌출부로 나타난다. χ3 시트는 Feng 등이 보고한 S2 상과 일치하는 마름모꼴 덤벨 패턴을 보인다. 바이어스 의존 시뮬레이션(그림 6)은 두 스트라이프 구성을 추가로 구별한다: 버클드 삼각형 스트라이프는 양/음 바이어스에서 대비가 반전되는 반면, β12 스트라이프는 바이어스에 관계없이 타원형 형태를 유지하여 Mannix 스트라이프를 버클드 삼각형으로, Feng 스트라이프를 β12로 할당함을 뒷받침한다.
Ag(111) 표면 위 독립 및 에피택셜 붕소 시트의 시뮬레이션 STM 이미지. 독립 a 삼각형, b β12, c χ3 붕소 시트. 독립 d 삼각형, e β12, f χ3 붕소 시트의 부분 전하 밀도. g 버클드 삼각형, h β12, i χ3, j χ3' Ag(111) 위 붕소 시트. 바이어스 전압 1.0 V. 녹색 구는 붕소 원자. 실선 빨간 선으로 둘러싸인 직사각형과 마름모는 독립 및 성장 보로핀 시트의 단위 셀. 실험적으로 관찰된 k 스트라이프 상(Ref.[22]), l 스트라이프 상(Ref.[23]), m 균질 상(Ref.[23])
요약하면, 제1원리 계산은 독립 버클드 삼각형, β12, χ3 보로핀이 열역학적으로 불안정하고 금속적임을 입증한다. Ag(111) 기판이 이 시트들을 안정화하여 붕소 원자당 0.1~0.2 eV 에너지를 낮춘다. 시뮬레이션된 STM 이미지는 관찰된 표면 상과 근본적인 격자 구조 간의 결정적 매핑을 제공하여 Ag(111) 위 보로핀의 원자 구성에 대한 오랜 논쟁을 해결한다. 이러한 통찰은 보로핀 기반 소자의 합리적 설계와 금속 기판 위 미래 실험적 합성 지침의 기초를 마련한다.