Optimization-Based Thickness Estimation and Depth-Resolved FIB-Induced Damage Characterization via Multislice Electron Ptychography

기술 요약: 멀티슬라이스 전자 프티코그래피를 이용한 최적화 기반 두께 추정 및 깊이 분해능 기반 FIB 유도 손상 특성 분석

문제 정의
정량적 전자 현미경 기술, 특히 투과 전자 현미경(TEM) 및 주사 투과 전자 현미경(STEM)은 시편 두께에 대한 정확한 지식에 크게 의존한다. 동역학적 산란(Dynamical scattering)은 이미지 대비와 회절 강도가 전파 거리(propagation distance)에 비선형적으로 의존하게 만들기 때문이다. 멀티슬라이스 전자 프티코그래피(MEP)에서 시편 두께는 순방향 모델(forward model)을 위한 필수 입력 파라미터이다. 그러나 집속 이온 빔(FIB) 밀링으로 제작된 특정 부위 시편의 경우, 국부적인 두께를 알 수 없는 경우가 많으며 공간적으로 변화한다(예: 쐐기 모양의 라멜라). 또한, 두께는 프로브 디포커스(probe defocus) 및 슬라이스 개수와 같은 다른 재구성 파라미터와 결합되어 있어, 재구성 자체로부터 독립적으로 결정하는 것이 쉽지 않다. 게다가 FIB로 제작된 라멜라는 일반적으로 이온 주입 및 손상으로 인해 형성된 비정질 표면층에 둘러싸인 결정질 벌크(crystalline bulk)로 구성된다. 전통적인 투영 이미지에서는 이러한 층들이 빔 방향을 따라 결정질 벌크와 중첩되므로, 별도의 단면 시편을 준비하지 않고는 이들의 기여도를 구별하거나 손상 깊이를 정량화하기 어렵다.

방법론
저자들은 단일 4D-STEM 데이터셋으로부터 시편 두께, 프로브 디포커스, 슬라이스 개수를 공동으로 추정하기 위해 멀티슬라이스 전자 프티코그래피를 베이지안 최적화(Bayesian optimization, BO)와 결합하는 프레임워크를 제안한다.

1. 시편 준비: 듀얼 빔 FIB-SEM을 사용하여 Si (001) 웨이퍼로부터 쐐기 모양의 실리콘 TEM 라멜라를 제작하였다. 각 라멜라는 최종 밀링 전압을 2, 5, 8 또는 30 kV로 달리하여 암모퍼스화(amorphization)의 정도를 체계적으로 변화시켰다.
2. 데이터 획득: 직접 전자 검출기를 사용하여 200 kV에서 4D-STEM 데이터셋을 획득하였다. 데이터셋은 30 mrad의 프로브 수렴 반각(convergence semi-angle)을 가지며 9.2 × 9.2 nm²의 시야를 커버한다.
3. 베이지안 최적화 (BO): 세 가지 파라미터인 총 시편 두께($t$), 프로브 디포커스($\Delta f$), 슬라이스 개수($N$)를 공동 최적화하기 위해 BO 엔진을 사용하였다. 목적 함수는 측정된 회절 진폭과 시뮬레이션된 회절 진폭 사이의 제곱평균제곱근(RMS) 차이의 음의 로그로 정의되었다. 최적화 과정에서 잘못된 두께가 다른 자유도에 의해 보상되는 것을 방지하여, 추정된 두께가 수치적 아티팩트가 아닌 물리적으로 의미 있는 값이 되도록 하기 위해 여러 프로브 상태와 직교 프로브 완화(orthogonal probe relaxation)를 비활성화하였다.
4. 깊이 분해능 기반 특성 분석: 최적화된 파라미터를 사용하여 전체 MEP 재구성을 수행하였다. 결과로 얻은 3D 볼륨을 분석하여 결정질 내부와 FIB 유도 비정질 표면층을 분리하였다. 비정질 층의 두께는 대칭적인 이중 오차 함수(symmetric double error function)를 사용하여 원자 열(atomic columns)의 위상 대 깊이 프로파일을 피팅함으로써 결정질-비정질 계면을 식별하여 정량화하였다.
5. 검증: BO-MEP 데리베이션된 두께 추정치를 로그 비율법(log-ratio method)을 이용한 전자 에너지 손실 분광법(EELS)의 독립적인 측정값과 비교하였다.

주요 결과

• 두께 추정: 공동 최적화는 두께, 디포커스 및 슬라이스 개수에 대해 잘 정의된 고유한 최소값을 산출하였다. 쐐기형 라멜라의 서로 다른 하부 영역에 적용했을 때, 최적화된 두께 값은 아래로 갈수록 두꺼워지는 예상되는 공간적 경향을 재현하였다.
• EELS와의 비교: BO-MEP 두께 추정치는 EELS 측정값에 비해 10–17%의 체계적인 편차(두께 과소 추정)를 보였다. 그러나 저자들은 EELS의 절대 두께가 약 ±20%의 정확도를 가진 가정된 비탄성 평균 자유 행로(inelastic mean free path)에 의존한다는 점을 언급하였다. 이러한 편차에도 불구하고, 이 방법이 재구성 설정과 독립적으로 국부적인 두께 변화를 재현할 수 있다는 점은 이 방법이 물리적 추정치로서 해석될 수 있음을 뒷받침한다.
• 강건성(Robustness): 최적화된 두께는 레이어 정규화 강도(layer regularization strength)의 변화에 대해 강건한 것으로 나타난 반면, 최적의 슬라이스 개수는 일관된 축 방향 샘플링을 유지하기 위해 이동하였다.
• 손상 특성 분석: 깊이 분별 재구성은 비정질 표면층을 결정질 벌크로부터 성공적으로 분리하였다. 정량화된 비정질 층의 두께는 밀링 전압에 따라 단조롭게 증가하였다: 3.2 nm (2 kV), 5.0 nm (5 kV), 7.0 nm (8 kV), 26.4 nm (30 kV). 이 값들은 기존에 보고된 단면 측정값과 합리적으로 일치하였으나, 30 kV 영역은 더 두꺼운 시편의 복잡성과 확장된 왜곡 영역으로 인해 더 큰 변동성을 보였다.

의의 및 주장
본 논문은 단일 4D-STEM 데이터셋에 대한 멀티슬라이스 전자 프티코그래피와 베이지안 최적화의 적용이 자기 일관적이고 두께를 인지하는 정량적 분석을 향한 경로를 제공한다고 주장한다. 주요 의의는 다음 두 가지 측면에 있다:

1. 물리적 해석 가능성: 본 연구는 공동 최적화를 통해 파라미터 결합을 해결함으로써, 최적화에 의해 선택된 두께가 단순히 재구성 오차를 최소화하는 수치적 파라미터가 아니라 의미 있는 국부적 물리 시편 두께의 추정치로 해석될 수 있음을 입증한다.
2. 통합된 손상 평가: 이 방법은 원자 해상도 분석에 사용된 라멜라로부터 국부 시편 두께와 깊이 분해능 기반의 FIB 유도 손상(비정질 층 두께)을 동시에 결정할 수 있게 한다. 이는 손상 깊이나 총 두께를 평가하기 위해 별도의 단면 시편을 준비할 필요를 없애준다.

저자들은 이 접근 방식이 두께에 대한 거친 사전 추정치(coarse prior estimate)만을 사용하여 탐색 범위를 정의하는 정밀화 도구(refinement tool)로 자리매김하고 있으며, 따라서 시뮬레이션을 제약하고 동역학적 산란의 영향을 많이 받는 측정값의 해석을 개선하기 위해 일상적인 MEP 워크플로우에 통합하기 적합하다고 보고 있다.