Insights on molecular P implantation for scalable spin-qubit arrays

기술 요약: 확장 가능한 스핀 큐비트 어레이를 위한 분자 P 이온 주입에 관한 통찰

문제 정의
실리콘 기판 내 인 도너($^{31}$P)의 스핀을 활용하는 고체 상태 스핀 큐비트는 양자 컴퓨팅을 위한 유력한 후보이다. 스핀 결맞음(coherence)은 크게 개선되었으나, 결정적인 과제는 실리콘 격자 내에 도너 원자를 정밀하고 확장 가능하게 배치하는 것이다. 최적의 큐비트 작동을 위해 도너는 표면으로부터 7~20 nm 아래에 위치해야 한다. 현재의 이온 빔 주입 기술은 확장성은 제공하지만, 배치 불확정성이 주사 터널링 현미경(STM) 리소그래피보다 수 자릿수 더 높다는 문제가 있다.

이를 완화하기 위해, 최근에는 단원자 인 이온 대신 분자 이온(특히 $PF_2$)을 사용하는 방안이 제안되었다. 가설은 분자 이온이 검출을 위한 높은 전자 신호를 생성하면서도, 9 keV 단원자 P 이온과 유사한 배치 불확정성을 유지할 수 있다는 것이다. 그러나 분자가 충돌 즉시 분해된다고 가정하는 기존의 이진 충돌 근사(BCA) 방식은 다체 효과(many-body effects), 충돌 중 분자의 무결성, 전자 신호의 정확한 성격, 그리고 분자 주입과 단원자 주입에 의해 발생하는 구체적인 방사선 손상 프로파일을 포착하지 못한다.

방법론
저자들은 단원자 및 분자($PF_2$) 이온의 실리콘 조사 차이를 조사하기 위해 대규모 원자/분자 병렬 시뮬레이터(LAMMPS)를 이용한 분자 역학(MD) 시뮬레이션을 채택하였다.

• 시뮬레이션 설정: 침투 역학(채널링 포함)을 추적하기 위해 결정질 실리콘(c-Si) 구조를 가진 대형 시뮬레이션 셀($43.5 \times 43.5 \times 65.7$ nm$^3$)을 사용하였다. 또한 웨이퍼에서 발견되는 천연 산화물 층을 모사하기 위해 5 nm 두께의 비정질 실리카($a$-$SiO_2$) 층이 포함된 작은 셀을 사용하였다.
• 포텐셜: 상호작용은 Si-Si 상호작용을 위한 Stillinger-Weber(SW) 포텐셜(계산 효율성과 머신러닝 포텐셜에 필적하는 정확성을 위해 선택됨), 이온-격자 상호작용을 위한 ZBL 포텐셜, 그리고 분자 내 상호작용을 위한 Lennard-Jones 포텐셜을 사용하여 모델링되었다.
• 조건: 시뮬레이션은 300 K에서 수행되었다. 구성 성분인 P 원자의 에너지가 9 keV가 되도록 $PF_2$ 분자 이온에 총 20.041 keV의 에너지를 할당하였다. 채널링 확률을 최소화하기 위해 입사각은 7도로 설정하였으며, 산화물 층의 존재 여부를 변화시켜 에너지 손실 및 분자 무결성에 미치는 영향을 평가하였다.
• 분석: 저자들은 에너지 손실, 전자 저지능(electronic stopping power), 결함 형성(Wigner-Seitz 분석 사용), 그리고 결함 클러스터링을 추적하였다. 검출에 사용되는 전자-정공 쌍의 생성을 시뮬레이션하기 위해 전자 저지능을 모델링하였다.

주요 결과

1. 분자 무결성 및 에너지 손실:
   분자가 표면 충돌 즉시 분해된다는 BCA의 가정과 달리, MD 시뮬레이션 결과 $PF_2$ 분자는 기판 내부로 몇 나노미터 동안 온전하거나 부분적으로 온전한 상태($PF$ 복합체 형태)를 유지하는 경우가 많음을 밝혔다. 순수 c-Si 내에서는 40% 이상의 경우에서 6 nm 깊이에 도달하기 전에 하나의 불소(F) 원자가 탈착되었으나, 나머지 하나는 인(P) 근처에 머물러 있었다. $a$-$SiO_2$ 층의 존재는 충돌 빈도를 높여 더 이른 해리를 유발한다.

• 에너지 퍼짐: P 이온은 $a$-$SiO_2$ 층을 통과할 때 F 이온보다 약 1 keV 더 많은 에너지를 잃는다.
• 주입 프로파일: 분자 주입은 단원자 주입에 비해 더 평탄한 깊이 프로파일을 보이며, P 원자가 더 깊게 밀려 들어가는 경향이 있다. 이는 분자가 해리될 때까지 하나의 더 무거운 투사체로서 작형하기 때문이다. 특히, 분자 주입은 채널링된 이온의 "긴 꼬리(long tail)"를 더 강하게 생성하는데, 이는 산화물을 통과한 후 P 원자가 더 높은 유효 에너지를 가지고 c-Si에 진입하기 때문인 것으로 보인다.

2. 전자 신호 및 채널링 검출:
   연구는 전자 저지 신호(주입 깊이 검출에 사용됨)와 P 원자의 최종 위치 사이의 상관관계를 분석하였다.

• 신호 상관관계: 전자 신호 강도와 P 원자 단독의 최종 침투 깊이 사이에는 약한 상관관계(피어슨 계수 $r \approx 0.5$)가 발견되었다. 그러나 신호는 세 가지 분자 구성 성분(P + 2F)의 평균 침투 깊이와 강하게 상관된다.
• 채널링의 모호성: 강한 전자 신호가 반드시 P 원자의 채널링을 의미하는 것은 아니다. 신호는 $PF_2$ 분자의 세 원자 중 어느 것에 의해서도 생성될 수 있기 때문이다. 반대로, 저자들은 낮은 신호는 모든 입자에 대해 채널링이 발생하지 않았음을 의미한다고 언급하였다. 즉, 강한 신호는 모호할 수 있지만(예: P는 적절히 배치되었으나 F가 채널링되어 신호가 강한 경우 등), 낮은 신호는 어떤 입자에 대해서도 채널링이 일어나지 않았음을 확실히 나타낸다.
• 검출 한계: 신호 임계값을 사용하여 채널링 이벤트를 구분하려는 시도는 상당한 불확실성을 초래한다. 약 절반의 이벤트는 신호 강도만으로는 채널링 여부를 제대로 라벨링할 수 없다. 다만, 채널링 위험을 피하기 위해 강한 신호를 보이는 이벤트를 안전하게 폐기할 수는 있으나, 이 과정에서 올바르게 주입된 P 원자까지 버리게 될 가능성이 있다.

3. 방사선 손상 및 결함 진화:

• 시너지적 손상: 분자 주입은 세 개의 독립적인 단원자 이온의 합보다 더 많은 결함을 생성하며, 결함 수의 분포도 더 넓다. 이는 서로 가깝게 배치된 분자 원자들의 겹쳐진 충돌 캐스케이드가 더 광범위한 손상을 만들어내는 "시너지적 손상 효과"를 나타낸다.
• 클러스터링: 분자 주입 시 결함은 P 원자를 중심으로 하는 더 큰 클러스터(일부는 500개 이상의 원자 포함)를 형성하는 경향이 있다. 이러한 큰 클러스터들은 종종 공공(vacancy)이 풍부한 비정질 포켓 형태를 띤다.
• 산화물 층의 역할: $a$-$SiO_2$ 층은 원자의 진입 에너지를 감소시켜 가장 큰 결함 클러스터의 크기를 제한하며, 이는 격자 회복을 위한 어닐링 요구 사항을 줄일 수 있다.
• 어닐링 시사점: 공공이 풍부한 큰 클러스터와 특정 결함 쌍(E-center 등)의 존재는 전기적 활성화 및 양자 특성에 위험 요소가 된다. 이러한 결함들은 회복을 위해 고온 어닐링이 필요하지만, 분자 이온의 특정 손상 프로파일은 해로운 원자 확산을 억란하기 위해 단원자 이온을 위한 표준 어닐링 프로토콜의 조정이 필요함을 시사한다.

의의 및 주장
본 논문은 현재의 BCA 기반 모델 및 실험적 해석에 내재된 여러 가정을 반박하며, 분자 $PF_2$ 주입에 대한 포괄적인 분자 역학 기반 평가를 제공한다고 주장한다.

• 즉각적 해리 가설 반박: 연구는 분자가 표면에서 즉시 분해된다는 가정이 부정확함을 입증하였다. 분자는 해리되기 전 몇 나노미터를 관통하는 경우가 많으며, 이는 겹쳐진 캐스케이드와 시너지적 손상을 초래한다.
• 신호 해석: 저자들은 높은 전자 신호와 P 원자의 채널링을 직접 연관시키는 것에 주의를 촉구한다. 분자 주입에서는 신호가 어떤 구성 원자로부터도 발생할 수 있으므로, 이를 P 도너의 깊이나 채널링 상태를 나타내는 신뢰할 수 있는 단독 지표로 사용하는 것은 부적절하다. 반면, 낮은 신호는 모든 입자에 대한 비채널링의 신뢰할 수 있는 지표이다.
• 손상 프로파일: 본 연구는 분자 주입이 단원자 주입에 비해 더 큰 클러스터 형태의 결함을 특징으로 하는 독특한 손상 지형을 생성함을 강조한다. 이는 도너가 잔류 변형이나 큐비트 성능을 저하시키는 결함 없이 전기적으로 활성화되도록 하기 위해 사후 주입 어닐링 전략의 재평가가 필요함을 의미한다.

요약하자면, 분자 이온은 확장 가능한 큐비트 생산을 위한 경로를 제공하지만, 본 논문은 그 사용이 배치 검출 및 손상 회복에 관한 복잡성을 수반하며, 특히 전자 신호의 해석과 열적 어닐링 프로세스의 최적화 측면에서 세심한 관리가 필요하다고 결론짓는다.