2D GeP₃ 반도체 조절: 도펀트 원자가 전자가 전자적 특성 형성

초록

밀도 함수 이론 내 제1원리 계산으로 2D GeP₃ 단일층에 Ⅲ-Ⅵ족 도펀트 치환 도핑 시 전자적 특성과 안정성 조사. 도전 특성이 도핑 부위와 도펀트 원자가 전자 수에 의해 극적으로 수정됨. 구체적으로 Ge 자리 치환 시 도펀트 원자가 전자 수 함수로 금속-반도체 진동 발생, P 자리 치환 시 완전 역전. 공도핑 경우 'AND' 논리 현상 발견: 공도핑 GeP₃ 도전 특성이 단일 도핑 결과로 단순 논리 관계 통해 추론 가능. 또한 도펀트 형성 에너지 연구로 전자-정공, 정공-정공 공도핑 시스템이 쿨롱 인력으로 훨씬 에너지적으로 유리함 확인. 본 연구는 2D 도핑 현상 포괄적 이해 제시와 2D 이진 반도체 전자적 특성 조절 흥미로운 경로 제안.

서론

그래핀 발견 이래 TMD, 실리신, 저마넨, 포스포린, 텔루렌 등 2D 결정 군이 독특한 전기적, 광학적, 자기적 특성으로 큰 주목. 그래핀은 무질량 디랙 페르미온 거동으로 초고이동도 보이며 양자 스핀 홀 효과, 향상된 열전효과, 초전도성, 양자 이방성 홀 효과 가능성. 결정 구조 데이터베이스 증가와 현대 계산 도구로 미발견 2D 물질 탐구 활발. 현재 1000개 이상 2D 물질 예측, 일부 실험 제작되어 물리, 화학, 재료 과학 흥미 분야. 2D 물질 기초 연구는 센싱 분야 잠재 응용도 촉진.

최근 Jing 등이 새로운 2D 물질 GeP₃ 단일층 보고: BP 단일층보다 높은 화학적 안정성, 우수한 전자적, 광학적 특성. 강한 층간 양자 구속으로 반도체 특성, 약 0.55eV 조절 가능 밴드갭. 높은 용량과 양호한 순환 안정성으로 GeP₃ 박막이 리튬이온 배터리 유망 양극 제안. Li 등 GeP₃ 나노리본 연구로 폭 증가에 따른 밴드갭 짝-홀 진동 발견.

도핑은 2D 층상 물질 전자적, 자기적 특성 근본 조절 실용 전략. 단일 물질 한계 극복해 다분야, 소자 응용 확장. 2D 단일층 반도체에서 현저히 향상된 전자-전자 상호작용이 큰 밴드갭 재규격화와 엑시톤 생성 유도(다체 이론 계산과 실험 입증). 벌크 반도체 도핑 대비 2D 도핑도 강한 전자 구속 효과로 이상 거동 기대: 붕소/질소 도핑 그래핀 디랙 포인트 미세 밴드갭 개방, 작은 BN 도메인 도입으로 K(K') 포인트 주변 밴드갭 효과적 개방. 흑린 원자가 전자 수 짝/홀 다른 원소 도핑 시 밴드갭 진동 거동. 본 연구는 2D 이진 GeP₃ 단일층 반도체에서 Ⅳ-Ⅴ족 도핑 원소 조사 확장.

여기서 Ⅲ-Ⅵ족 도펀트로 치환 도핑된 GeP₃ 단일층 체계적 연구 수행. 도핑 시스템 전자적 특성이 도펀트 원자가 전자 수와 도핑 부위 둘 다에 의해 극적 영향. 핵심: (1) 단일 도펀트 결과가 치환 부위와 민감하게 의존, 두 도핑 부위 치환 시 완전히 역전 결과. (2) 공도핑 도전 특성이 단일 도펀트 결과로 논리 연산자 통해 추론 가능. 추가로 다양한 도핑 유형 형성 에너지 계산으로 일부가 열적 변동 대비 고도로 에너지 유리함 제시.

계산 방법

일반 기울기 근사 내 모든 밀도 함수 이론 계산은 비엔나 압시니오 시뮬레이션 패키지(VASP)로 수행. 교환 상관 항은 PBE 범함수, 투영 증강 파동 포텐셜로 전자-이온 상호작용 기술. 도핑 GeP₃ 단일층을 32원자 주기적 2×2 슈퍼셀로 모델링, 3×3 더 큰 슈퍼셀로 결과 검증. z 방향 약 20Å 진공 공간으로 인접 층 간 상호작용 제거. 단일 도핑 시 Ge 또는 P 원자 하나를 Ⅲ(Ⅳ, Ⅴ, Ⅵ)족 도펀트로 치환. 기하 구조는 격자 상수와 호스트 GeP₃ 단일층 전자적 특성 등 보고 결과 비교로 결정. 도핑 시스템에서 슈퍼셀 모든 원자 헬만-파인만 힘 0.02eV/Å 미만까지 완화, 표면 셀 격자 상수는 원자 완화 중 고정. 운동 에너지 컷오프 약 600eV, 6×6×1 k-메시 사용.

GeP₃ 단일층 내 도펀트 가용성 확인 위해 도펀트 X(X=Ⅲ-Ⅵ족) 형성 에너지(E_f) 다음 두 공식으로 계산. 단일 도펀트:

E_f(GeP₃:X) = E(GeP₃:X) - E(GeP₃) - E_X + E_i (1)

공도핑 시스템 유사 공식:

E_f(GeP₃:XY) = E(GeP₃:XY) - E(GeP₃) - E_X - E_Y + E_i + E_j (2)

여기서 E_f(GeP₃:X), E(GeP₃) 같은 슈퍼셀 X-도핑 및 고유 GeP₃ 단일층 총 에너지. E(GeP₃:XY) XY 공도핑 시스템 총 에너지. E_X, E_Y 도펀트 X, Y 원자 에너지(해당 벌크 구조 기준). E_i, E_j 치환 원자 에너지(i, j는 Ge 또는 P 원자 지시).

결과 및 고찰

단일 원소 도핑 시스템 짝-홀 진동

그림 1a는 GeP₃ 2×2 슈퍼셀 위/옆면 구조, 그림 1b 해당 2D 브릴루앙 영역. GeP₃ 단일층 최적화 격자 상수 a=b=6.96Å, 계산 밴드갭 약 0.26eV로 다른 이론 계산과 잘 일치.

먼저 Ge 원자 치환 단일 원소 도핑 GeP₃ 단일층 밴드 구조 플롯(도펀트: B, C, N, O, Al, Si, P, S, Ga, As, Se). 결과 그림 2a-l. 페르미 준위가 Ⅴ족(N, P, As) 도펀트에서 전자 하나 더해 전도대 건너며 상승, Ⅲ족(B, Al, Ga) 전자 하나 적어 원자가대 건너며 하강. 예: 그림 2f, j 원자가대 최대가 그림 2e, i 부분 점유 밴드와 대응. Ⅳ족(C, Si, Ge)과 Ⅵ족(O, S, Se) 도펀트는 Ge 원자와 같거나 두 개 더 많은 전자로 반도체 특성 유지. 반도체-금속 전이 조절은 원자가 전자 수 점유에 기인: 홀수(짝수) 원자가 전자 점유가 금속(반도체) 특성 유도.

PBE 범함수 유도 결과 타당성 확인 위해 최상단 행 도핑 시스템에 하이브리드 밀도 범함수(HSE06)도 사용. PBE가 밴드갭 과소평가로 밴드갭 오류 야기함 명확. 그러나 연구 시스템 모두 상당한 갭 가져 PBE 오류로 금속/반도체 특성 실수 드묾(일부 작은 밴드갭 반도체에서 PBE가 전도/금속 특성 실수 유발). 본 연구 관심사는 특정 밴드갭 값 아닌 금속/반도체 특성. HSE06 갭이 PBE 대비 명확히 확대되나 금속-반도체 진동 유지. 따라서 PBE 기반 핵심 결론 신뢰성 있음.

그러나 P 원자 치환 같은 도펀트 경우 완전히 역전(그림 3a-l). Ⅴ족(N, As)과 Ⅲ족(B, Al, Ga) 도펀트에서 반도체 특성 유지, Ⅳ족(C, Si, Ge)과 Ⅵ족(O, S, Se)에서 금속 특성 변화. 이는 원자가 전자가 Ⅴ족(Ⅲ족) 도펀트에서 고유 GeP₃와 같음(두 개 적음), Ⅳ족(Ⅵ족)에서 하나 적음(많음) 때문.

반도체-금속 특성 전이 진동 더 잘 제시 위해 도펀트별 밴드갭 변화 추세 플롯(그림 4a,b). 반도체-금속 전이 급격히 역전됨 명확. 구체적으로 Ge 자리 치환 시 금속(반도체)-반도체(금속) 진동, P 자리 치환 시 역전. 또한 도펀트가 Ge 원자와 같은 원자가 전자 수 가질 때 밴드갭 크기 고유 GeP₃와 거의 같음. 두 개 더 많은 전자 가질 때 밴드갭 변화 상대적 큼. 반면 P 자리 도펀트에서는 원자가 전자 수 무관하게 밴드갭 변화 항상 상대적 큼. 이는 원자 반경과 가용 원자가 전자 공동 효과 이해 가능: Ge 원자와 거의 같은(작거나 큰) 반경과 원자가 전자 가진 도펀트가 전자적 특성(밴드갭 등)에 상대적으로 작은(큰) 효과. 즉 적절한 도펀트와 다른 도핑 부위 선택으로 반도체-금속 전이 진동뿐 아니라 밴드갭 크기도 조절 가능.

GeP₃ 단일층 내 다른 도펀트 전자 구조 변화 이해 위해 고유 및 Ⅳ-Ⅴ족 도핑 GeP₃ 부분 상태 밀도(PDOS) 플롯(그림 5a-d). GeP₃ 단일층 원자가대 최대(VBM)와 전도대 최소(CBM) 주로 Ge와 P p 오비탈 기원. Ge 원자와 같은 원자가 전자 수 도펀트(C, Si) 사용 시 p 오비탈 에너지 준위가 P 원자보다 높아 고유 GeP₃ VBM 바로 위에 불순물 상태 위치(그림 5a). 따라서 전도 특성 유지, 밴드갭 변화 크기 상대적 작음. 그러나 도펀트가 Ge보다 하나 더 많은 전자(N, P, As) 가질 때 밴드갭 내 불순물 상태 존재, 불순물 상태는 분열 CBM(주도적)과 도펀트 상태 기원(그림 5b).

반면 P 자리 도핑 시 도펀트가 P 원자보다 하나 적은 원자가 전자(Ⅳ족) 가질 때 페르미 준위 가로지르는 불순물 상태, 불순물 상태는 분열 VBM(주도적)과 도펀트 상태 구성. 도펀트가 P와 같은 원자가 전자(Ⅴ족) 가질 때 반도체 특성 유지(그림 5c). 격자 상수 큰 불일치로 밴드갭 고유 GeP₃보다 상대적 큼. 또한 P 자리 치환 시 N 도펀트 밴드갭이 As 도펀트보다 큰 것 관측: As 원자 p 오비탈 에너지 준위가 N보다 높아 더 높은 에너지 준위 p 오비탈일수록 불순물 상태가 VBM에서 더 업시프트(그림 5d).

공도핑 시스템 논리 관계

앞선 단일 도펀트 발견 기반으로 원하는 전자적 특성 만족 공도핑 시스템 설계 가능. 여기 B, C, N, O 예시로 공도핑 효과 설명, 결론은 다른 선택 도펀트에도 강건. 예: Ge 자리 공도핑 시 두 도펀트 모두 하나 적은 원자가 전자면 자연스럽게 반도체 특성, 하나 적고 하나 많은 원자가 전자 도펀트면 공도핑 시스템도 반도체 특성 가능.

그러나 하나 적음(많음)과 같은 원자가 전자 두 도펀트면 여전히 금속 특성 유지(하나 적음(많음) 원자가 전자 도펀트 결과 특성). 단순화: 이 아이디어가 추가 DFT 계산으로 정확히 확인됨, 그림 6a-l B, C, N, O 공도핑 GeP₃ 단일층 밴드 구조 결과 참조.

이제 금속 특성 'M', 반도체 특성 'S'로 설정해 논리적 'AND' 연산 예시: M AND M = S, S AND S = S, M AND S = M 정의. 위 발견이 이런 논리적 관계 준수. 예: 하나 많고 하나 적은 원자가 전자 도펀트 금속 특성 유도, 하지만 두 도펀트 공도핑 시 Ge 자리 B, N(그림 7a,b) 공도핑 시스템 반도체 특성 됨(그림 6c). B-C 공도핑 GeP₃ 단일층 시스템 선택 시 M AND S 경우 금속 특성(그림 4a,b). C-N, N-O, B-O 원자 GeP₃ 공도핑도 두 Ge, 두 P, 또는 Ge-P 치환 시 동일(그림 7c-f).

마지막으로 단일 도핑, 공도핑 시스템 안정성 확인 위해 실험 추가 실현 가능성 검토. 형성 에너지 식 (1), (2)로 단일, 공도핑 계산. 결과 그림 8a-e. 그림 8a,b에서 Ge 자리 단일 도펀트 E_f 모두 GeP₃ 단일층(기준점 0 설정) 근접, C, N, S 원자 제외. B, O, P, Ge, Se 원자 도펀트 형성 에너지 훨씬 작아 실험 도핑 용이. P 자리 도펀트 중 B, O, P, Ge 원자 상대적 작은 형성 에너지로 도핑 용이. C, N, Al, Ga 도핑 어려움.

공도핑 형성 에너지 관해 그림 8c-e는 두 도펀트 위치 Ge-Ge, P-P, Ge-P 사이트 공도핑 형성 에너지. B, C, N, O 도펀트 예시만 보임. Ge-Ge, P-P 사이트 공도핑 형성 에너지는 단일 원소 도핑 형성 에너지 평균으로 대략 추정 가능. 명확히 Ge-Ge 사이트 BB, BO, OO 공도핑, P-P 사이트 BB, BO, NO, OO 공도핑 형성 에너지 상대적 작아 실험 실현 용이. 그러나 Ge-Ge 사이트 CC, CN, NN, P-P 사이트 CC, CN 공도핑 형성 에너지 상대적 커 도핑 어려움. Ge-P 사이트 공도핑(그림 8e)은 도펀트 간 전하 이동으로 Ge-Ge, P-P보다 복잡. 어쨌든 BB, BO, OO 공도핑 작은 형성 에너지, CC, CN, NN 큰 형성 에너지. 일반적으로 형성 에너지는 도펀트 원자가 전자 수에 고도로 의존. 구체적으로 치환 원자보다 하나 적음(많음) 전자 두 도펀트면 공도핑 시스템 형성 에너지가 해당 단일 도펀트 BB(NN) 공도핑 Ge 자리보다 낮음(높음). 이는 감소(증가) 전자 에너지 이득과 쿨롱 반발 경쟁 때문. 정공-정공 공도핑 경우 전자 이득이 후자보다 훨씬 커 공도핑 시스템 형성 에너지 현저 감소(BB 등). 전자-전자 공도핑 경우 둘 다 높은 형성 에너지 유도(NN 등). 그러나 정공-전자 공도핑 시스템(BN 공도핑 Ge 자리 등) 형성 에너지가 단일 도핑 대비 극적 감소. 이런 공도핑 시스템에서 순 추가/감소 전자 에너지 이득 없고 쿨롱 상호작용이 공도핑 도펀트 형성 결정적 역할. 종합하면 흑린 원소 도핑 선행 연구 고려 시 현재 연구 일정 보편성 가지며 BN, MoS₂ 등 다른 2D 반도체 단일층 적용 기대.

마지막으로 위 도핑 시스템 비도핑 대비 안정성 확인 위해 AIMD(ab initio 분자 동역학)로 에너지-시간 플롯(그림 9a-f). 시간이 충분히 길어짐(~4ps) 진동 진폭 수렴, C 도핑 GeP₃(그림 9a) 300K까지 열적 변동 대비 붕괴 안함 의미. 가장 활성 C 원자 도핑 GeP₃도 극한 온도 300K까지(그림 9c). 금속 Al Ge 자리 치환 예시로 그림 9e,f 계산 결과: 에너지 진동 진폭 시간 경과에 따라 점차 감소, 시간 충분히 길면 에너지 수렴 의미, 도핑 시스템 구조 300K 열적 변동 대비 열적 안정. 따라서 고품질 GeP₃ 단일층 제조 시 이런 도핑 시스템 추가 실험 실현 기대.

끝으로 본 연구 신뢰성 논의. 여기 제시 결론 이론적 예측이나 고신뢰도. 호스트 물질 이미 보고됨, 벌크 층상 GeP₃ 이미 존재[26]. 연구 도핑 유도 현상 단일층 GeP₃ 실현 후 실험 추가 확인 필요. 단순화하면 단일층 GeP₃ 전자/정공 도핑 일부 분자 흡착으로 실현 가능.

결론

요약: 2D GeP₃ 단일층 Ⅲ-Ⅵ족 도펀트 전자적 특성 조사. Ge 자리 치환 시 도펀트 원자가 전자 수 함수로 금속-반도체 진동, P 자리 치환 시 역전. 단일 도펀트 결과 기반 GeP₃ 공도핑 도전 특성 단순 논리 연산으로 획득 가능. 마지막으로 다양한 도펀트 형성 에너지 계산으로 일부 공도핑 시스템, 특히 전자-정공, 정공-정공 공도핑이 쿨롱 인력으로 훨씬 에너지 유리함 확인. 본 발견 새 현상 제시와 2D 이진 반도체 전자적 특성 조절 흥미로운 경로 제안.

데이터 및 재료 가용성

현재 연구 중 생성/분석된 데이터셋은 요청 시 교신 저자로부터 제공 가능.

약어

1D

1차원

2D

2차원

AIMD

어비니시오 분자 동역학

BP

흑린

CBM

전도대 최소

DFT

밀도 함수 이론

HSE06

하이브리드 밀도 범함수

PBE

퍼듀-버크-에른처호프

PDOS

부분 상태 밀도

VBM

원자가대 최대