육방정계 질화붕소 삼각형 공극: 크기 의존적 안정성, 자성, 에지 쌍 결손 효과
초록
밀도범함수이론(DFT)을 사용하여 단층 육방정계 질화붕소(h-BN) 삼각형 공극의 원자적, 전자적, 자기적 거동 조사. 계산 결과 완화된 삼각형 공극 기하구조가 크기에 강하게 의존하며, 모든 안정한 구성은 N-종단 지그재그 에지 가짐. 공극 에지에서 붕소-질소(BN) 쌍 제거 시 결과 구조 안정성과 자기 모멘트는 결손 쌍 위치에 따라 변하며, 면 영역에서 쌍 제거 시 가장 유리한 구성으로 최소 자기 모멘트 생성.
배경
육방정계 질화붕소(h-BN)는 그래핀과 유사한 2차원 절연체이나 B와 N 원자 번갈아 구성. 합성 중 h-BN은 필연적으로 점 결함, 특히 공극 포함하며 이는 전자적·기계적 특성 변화시켜 h-BN 기반 나노소자 성능 영향[1,2]. 격자 두 원자 유형으로 공극 에지는 N-종단 또는 B-종단 가능; 전자가 에너지적으로 선호되며 지그재그 모티프 채택 빈번[3,4]. 선행 이론 연구는 삼각형 공극 안정성과 자성이 에지 종단과 공극 크기 의존함 보여줌[3-11]. 실험은 전자빔 조사 하 삼각형 공극 형성 확인, 공극 성장하며 삼각형 모양 유지[12-15].
최근 연구에서 h-BN 원자가 공극 에지에서 개별 아닌 다발로 방출됨 관찰, 공극 성장 협동 메커니즘 시사[15]. 본 논문은 삼각형 공극 구조, 크기 의존적 완화, BN 쌍 제거 영향에 대한 포괄적 DFT 분석 제공.
계산 방법
모든 계산은 VASP(비엔나 어비니시오 시뮬레이션 패키지)로 평면파 컷오프 400 eV, 투영 증강파(PAW) 슈도포텐셜 사용[18,19]. 교환-상관 에너지는 Perdew-Burke-Ernzerhof 일반화 기울기 근사(GGA)로 처리, Grimme DFT-D2 반데르발스 보정 추가로 분산 상호작용 포착[20-22]. 원자 기하구조 잔여 힘 0.01 eV/Å 미만까지 완화; 브릴루앙 존 샘플링은 (9×9), (15×15) 슈퍼셀에 Γ-포인트만 사용. 계산 격자 상수 2.56 Å 실험값 일치[23].
결과 및 논의
h-BN 삼각형 공극
N-종단 에지가 B-종단보다 안정적[3,4]이므로 여러 삼각형 공극 크기 조사, N-종단 에지 집중. N-종단 삼각형에 대해 두 가지 뚜렷한 완화 기하구조 등장: 정점 원자 거의 변화 없는 대칭 N-symm 구성, 각 정점에서 N-N 동종 결합 형성하는 왜곡 NN-bond 구성. 가장 작은 공극(V1B)은 N-symm 구조만 발견, 2.66 Å N-N 거리로 원시 격자(2.48 Å) 대비 약간 확장. 공극 크기 V3B+1N, V6B+3N 증가 시 N-symm과 NN-bond 구성 모두 국소적 안정, 선행 연구와 일치[6].
더 큰 공극(V10B+6N 이상)에서는 NN-bond 구성이 유일한 안정한 기하구조. 이 영역에서 구멍 주변 B-N 결합 본질적으로 변화 없으며 정점 N-N 결합 지속. 표 1(미수록)에 N-symm과 NN-bond 구조 간 에너지 차이가 공극 크기 증가에 따라 감소, 정점 재구성 전체 격자 영향 감소 반영.
이 공극들 자기 모멘트는 에지 원자 종단과 상관. N-symm 구조에서 각 댕글링 N 원자 1 μB 기여, 모멘트 = 에지 N 원자 수. 반면 NN-bond 구조는 N-N 결합이 댕글링 결합 특성 소멸 또는 보존 여부에 따라 자성 감소 또는 증가. 예: V3B+1N, V6B+3N, V10B+6N NN-bond 상태 자기 모멘트 각각 0, 3, 6 μB.
공극 에지 BN 쌍 제거
실험 관찰 공극 성장이 에지에서 BN 쌍 방출로 진행[14,15]. 모방 위해 큰 V21B+15N 공극(15×15 슈퍼셀) 에지 다양한 위치에서 BN 쌍 제거. 여섯 개 결손 위치 조사: 두 모서리, 두 모서리 인접 면, 두 중앙 면. 완화 후 최적화 기하구조 세 클래스-모서리 결손, 면 결손, 중앙 결손-각 특징적 재구성 보임.
모서리 결손 구조는 약간 왜곡 BN 링 외 최소 변화. 반면 면 결손 구조는 공극 인접 오각형 링과 N 이합체 형성, N 이합체 안정성과 자성 크게 영향. 상대 에너지는 삼각형 중심 가까운 결손 쌍 구성 에너지 유리함 보여줌(표 2). 자기 모멘트 모서리 결손 12 μB에서 특정 면 결손 9 μB 범위, 비짝전자 재분포 반영.
스핀 밀도 맵(그림 3) 비짝 스핀이 에지 N 원자, 특히 왜곡 BN 링과 N 이합체에 국소화 확인. 상태밀도 분석(그림 4) 원시 h-BN 밴드 갭 내 결함 상태, 페르미 준위 주변 스핀 비대칭 특징 보여줌. 세 대표 결손 쌍 구성 밴드 갭 각각 0.35, 0.24, 0.36 eV.
결론
제1원리 연구는 h-BN 삼각형 공극 안정성과 자성이 공극 크기, 에지 종단, BN 쌍 제거 정확 위치에 의해 지배됨 입증. 큰 공극은 정점 NN-bond 구성 선호하며 이는 자기 모멘트 결정. 큰 삼각형 구멍 에지에서 BN 쌍 제거 시 가장 안정한 결과는 N-N 결합 형성 면 결손 구조이며, 자기 특성과 국소 상태밀도는 결손 위치 민감. 이 통찰은 전자 및 스핀트로닉 소자용 결함 공학 h-BN 나노구조 조정에 미시적 기반 제공.