극한 도핑과 스트레인 구현: 에피택셜 단결정 실리콘의 초고농도 보론 도핑

초록

본 연구는 나노초 레이저 도핑(GILD)으로 에피택셜 실리콘에 보론을 8 at.%까지 도핑하고, 격자 변형을 3 %까지 유도함으로써 전하 농도와 구조적 변형을 정량적으로 제어한다. 단순한 이항 모델과 DFT 계산을 통해 인접 보론 원자들의 복합체 형성이 활성 전하 농도의 내재적 한계를 만든다는 메커니즘을 제시한다.

상세 분석

이 논문은 초고농도 보론(보론) 도핑을 실현하기 위해 가스 침투 레이저 도핑(Gas Immersion Laser Doping, GILD)이라는 비열평형 공정을 적용하였다. 308 nm 파장의 25 ns 펄스 레이저를 사용해 실리콘 표면을 국부적으로 용융시킨 뒤, 용융된 액상에서 보론이 빠르게 확산·혼합되어 고속(≈4 m/s)으로 재결정화되는 과정에서 보론이 완전 치환(substitutional) 형태로 격자에 포획된다. 이때 레이저 펄스 수(N)를 1700번 조절함으로써 도핑 농도를 정밀하게 제어했으며, SIMS와 STEM 분석을 통해 20315 nm 두께의 레이어 전단에 걸쳐 보론 농도가 13 ± 8 % 수준의 균일성을 유지함을 확인했다.

전기적 특성은 Hall 측정을 통해 평가했으며, 보론 농도 C_B가 8 at.%(≈4 × 10²¹ cm⁻³)까지 증가함에 따라 정공 농도 h도 거의 100 % 활성화 비율(γ≈0.7)로 동반 상승하였다. 그러나 일정 농도(≈5–7 × 10²¹ cm⁻³)에서 정공 농도가 포화되기 시작했으며, 이는 격자 내 인접 보론 원자들이 2~3개가 인접한 자리에서 복합체(다이머, 트리머)를 형성해 전자를 비활성화시키는 현상으로 해석된다.

저자들은 이러한 현상을 설명하기 위해 간단한 이항(베르누이) 모델을 도입하였다. 보론 원자 하나가 네 개의 주변 실리콘 원자와만 결합될 확률을 (1‑p)⁴(p=C_B/n_Si)로 두고, 단일 보론(mononer) 농도 B₁ = C_B·(1‑p)⁴ 로 정의한다. 이 모델은 실험 데이터와 놀라울 정도로 일치하여, 보론 농도가 3 × 10²¹ cm⁻³ 이하에서는 활성 정공 농도가 B₁에 거의 동일함을 보여준다. 농도가 포화에 가까워질수록 B₁ 외에 B₂(다이머)와 B₃(트리머) 복합체가 비활성 정공을 차지하게 되며, 비활성 보론 농도 C_NA = C_B – h는 C_B²에 비례하는 2차 법칙을 따른다.

첫 원리 계산(Density Functional Theory, DFT)에서는 1~3개의 보론이 인접 격자점에 존재할 때의 형성 에너지를 평가하였다. 계산 결과, 첫·두번째 근접 이웃(shell S₁)과 1.5번째 근접 이웃(shell S₁.₅)까지 고려했을 때, 다이머와 트리머 복합체가 비교적 낮은 형성 에너지를 가져 실제 실험에서 관찰된 비활성화 메커니즘을 뒷받침한다. 또한, 약간의 보론 이동성(첫·두번째 근접 이웃 사이의 전이)만 허용해도 모델과 실험 사이의 차이가 최소화됨을 확인했다.

구조적 측면에서는 XRD와 STEM을 통해 격자 팽창이 보론 농도와 거의 선형적으로 증가함을 확인했으며, 3 % 수준의 인장 변형이 발생한다. 이 변형은 활성 보론뿐 아니라 비활성 복합체가 격자에 미치는 기여를 포함한다. 층이 얇을수록(≈20 nm) 완전 의사정질(pseudomorphic) 상태를 유지하지만, 두께가 170 nm 이상이면 부분적인 이완(relaxation)과 셀룰러 붕괴(cellular breakdown) 현상이 나타나 B‑rich 컬럼형 결함이 형성된다.

결과적으로, 이 연구는 (1) GILD 공정을 통해 실리콘에 8 at.% 수준의 초고농도 보론을 균일하게 도핑하고, (2) 단순한 확률 모델과 DFT 계산을 결합해 인접 보론 복합체가 전하 활성화와 격자 변형을 제한하는 근본 메커니즘을 정량적으로 규명했으며, (3) 이러한 제한은 결함이 전혀 없는 완전 치환 구조에서도 불가피하게 나타나는 ‘기하학적(geometrical) 한계’임을 입증하였다. 이는 그룹 IV 반도체 전반에 적용 가능한 보편적 원칙을 제시함으로써, 초고도 전자·광학 소자 설계와 초전도성 등 새로운 물성 탐색에 중요한 기반을 제공한다.