Method Article

극자외선 리소그래피에서의 주석 이물질 제어를 위한 통합 다중 방법 시뮬레이션 프레임워크

DOI:

10.3791/69818

March 27th, 2026

In This Article

Summary

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

이 프로토콜은 사용자가 극자외선(EUV) 및 신흥 Blue-X 리소그래피에서 주석 이물질 제어를 달성하기 위한 통합 시뮬레이션 프레임워크를 통해 안내하고, 운동학 모델링, 볼츠만 수송 방정식(BTE), 밀도 함수 이론(DFT) 기반 방법을 통합하여 이온 상호작용과 수소 보조 세척을 평가하는 것을 목표로 합니다.

Abstract

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

이 프로토콜은 대표적인 결과를 보여주는 개념적이고 통합된 모델링 프레임워크로, 사용자에게 볼츠만 수송 방정식(BTE), 입자-인-셀(PIC), 그리고 운동 시뮬레이션을 결합하여 극자외선(EUV) 리소그래피에서 주석(Sn) 잔해 완화를 연구하는 방법을 안내합니다. 프로토콜에는 Mo/Si 다층 미러(MLM)의 반사율, 스퍼터 수율, 이식 깊이, 운동학 모델링, BTE 계산이 포함됩니다. BTE 및 PIC 시뮬레이션은 수소 플라즈마의 전자 에너지 분포 함수(EEDF)를 해소하고, 다양한 플라즈마 조건에서 고에너지 Sn 이온의 생성 및 가속을 분석하는 데 사용됩니다. 수소 흐름이 이온 감속과 방사선 효율에 미치는 영향도 정량화되었습니다. SnxHy 종의 이온화 단면과 해리 채널을 바탕으로, Sn-H 충돌의 상호작용 전위는 밀도 함수 이론(DFT) 방법을 사용하여 이식 깊이를 계산합니다. 또한, MLM의 Sn 이물질과 MLM 내 Ru 코팅 간 상호작용에서 발생하는 MLM 반사율과 스퍼터 수율은 반경험적 공식을 사용하여 계산됩니다. 이 프로토콜을 따르면 사용자는 스퍼터 수율, 주입 깊이, MLM 반사율, 다양한 수소 플라즈마 EEDF 하에서의 SH4 형성 등 Sn 파편 제어와 관련된 주요 물리적 매개변수를 얻을 수 있습니다. 이러한 결과물은 EUV 리소그래피 시스템에서 오염, 세척, 검출 공정을 체계적으로 평가할 수 있게 합니다.

Introduction

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

극자외선 리소그래피(EUVL)는 집적 회로 소형화를 진전시키는 최첨단 기술로, 2nm 미만의 특징 파운팅을 가능하게 합니다. 일반적인 EUV 소스에서는 주석(Sn) 미세방울이 먼저 Nd:YAG 레이저의 프리펄스로 기화되고 이온화되며, 생성된 플라즈마 구름은 10.6 μm에서 작동하는CO2 레이저에 의해 재가열되어 EUV 방사선을 생성하며, 이는 Mo/Si 다층 거울(MLM)1,2에 의해 수집됩니다. ASML이 개발한 상용 시스템의 경우, 전원 전력은 대량 생산에 충분한 수준에 도달했습니다. 그럼에도 불구하고, 특히 중국에서의 연구는CO2 레이저 생성 Sn 플라즈마의 효율 향상에 계속 집중하고 있습니다.

EUV 광원의 주요 과제 중 하나는 고에너지 Sn 이온의 생성입니다. 고강도CO2 레이저 펄스로 Sn 방울을 조사하면 keV 범위의 에너지를 가진 이온이 생성되며, 이는 다층 거울(MLM)을 손상시키고 시스템 수명을 단축시킬 수 있습니다. 이온에 의한 손상을 줄이기 위해 수소(H-2)가 완충 가스로 널리 사용됩니다. 충돌 감속을 통해H2는 Sn 이온 수송을 완화하고 광학 부품에 도달하는 이물질을 줄입니다. 따라서 신뢰할 수 있는 제동 데이터와 Sn–H 상호작용의 정확한 모델은 소스 효율성과 내구성 최적화에 매우 중요합니다 5,6,7.

또 다른 중요한 문제는 진공 챔버 내 표면, 특히 플라즈마 근처에 위치한 집전기 거울에 Sn 조각이 침착되는 것과 관련이 있습니다. 얇은 Sn 코팅도 EUV 반사율을 낮추고 광학 성능과 작동 안정성을 저하시킵니다 8,9,10. 실용적인 산업 해결책으로는H2를 배경 가스11로 연속적으로 주입하는 방법이 있습니다. 이 방법에서는 수소 라디칼이 다음 발열 반응을 통해 Sn 코팅을 식각하여 휘발성 스탄난(SnH4)을 생성하고, 이를 펌핑으로 제거합니다.

Sn(s) + 4H(g) → SnH4(g),

이 방법은 Sn 제거를 향상시키는 데 효과적이지만, 새로운 합병증을 초래합니다. H-2 플라즈마 해리에서 생성된 수소 라디칼은 SnH4의 사슬 분해를 유도하여 Sn을 재생하고 2차 오염을 유발할 수 있습니다. 이러한 공정은 세척 효율을 떨어뜨리고 거울 안정성과 광학 성능을 저해할 수 있습니다. 따라서 주석 수소화물의 형성, 분해, 표면 상호작용에 대한 상세한 이해는 수소 기반 세척 방법 개선에 필수적입니다. 최근 표면 연구는 주석 수소화물과 그 중간체를 특성화하여 오염 경로를 정확히 식별하고 Sn재침착을 억제하는 데 매우 중요함을 강조합니다.

이러한 노력에도 불구하고, Sn–H 플라즈마 화학의 핵심 측면은 여전히 충분히 특성화되지 못하고 있습니다. 특히, EUV 관련 플라즈마 조건 하에서 Sn-H 종(예: Sn2H2 및 SnHx)의 구조, 반응성, 단편화, 형성/해리 속도는 직접적인 실험적 검증이 부족하다13. 더불어, Sn-H 플라즈마의 부작용 경로, 그 발생 확률을 결정하는 요인, 부작용의 임계 임계값, 장기 작동 안정성에 대해서는 체계적으로 조사되지않았습니다.

이 문제들을 종합해 보면, 원자물리학 및 분자 물리학, 플라즈마 물리학, 양자화학 관점에서 플라즈마-표면 상호작용에 대한 근본적인 연구가 필요함을 강조합니다. 기존 모델링 접근법은 일반적으로 Sn 이온 생성,H-2 의 고에너지 Sn 이온에 대한 차단력, 이온-표면 상호작용 등 Sn 이물질 제어의 고립된 측면만 다루므로, 오염-청소-검출 전체 과정을 포착하지 못합니다. 이러한 한계를 해결하기 위해 입자-인-셀(PIC) 시뮬레이션, 볼츠만 수송 방정식(BTE) 분석, 밀도 함수 이론(DFT), 운동 모델링을 결합한 통합 시뮬레이션 프로토콜을 개발하고자 했습니다. 극자외선(EUV) 광원 연구는 레이저-물방울, 레이저-플라즈마, 플라즈마-플라즈마, 플라즈마-가스 상호작용 등 여러 결합 과정을 포함합니다. 이 프로토콜은 유체역학, 입자-인-셀(PIC), 밀도 함수 이론(DFT) 방법을 결합한 통합 시뮬레이션 프레임워크를 설명하여 주석(Sn) 잔해 완화와 수소 정화를 모델링합니다. 이 프로토콜은 Sn 잔해 생성, 운송, 표면 상호작용 및 수소 보조 완화를 위한 통합되고 재현 가능한 워크플로우를 제공합니다. 다음 절에서는 이 방법론의 단계별 구현을 자세히 설명합니다.

Access restricted. Please log in or start a trial to view this content.

Protocol

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

참고: 유체, 운동학, 양자화학 접근법의 통합을 포함한 전체 워크플로우를 포함합니다. 작업 흐름은 그림 1 (빨간 박스에 강조됨)에 나와 있습니다.

figure-protocol-1
그림 1. 극자외선 리소그래피용 통합 시뮬레이션 프레임워크 개략도. 약어: MLM = 다층 거울; PIC = 입자-인-셀; BTE = 볼츠만 수송 방정식; EEDF = 전자 에너지 분포 함수. 이 그림의 더 큰 버전을 보시려면 여기를 클릭해 주세요.

1. MLM 반사율 시뮬레이션

  1. 다층 파라미터를 설정하세요. EUV 소스에서 Mo/Si MLM을 수집자로 활용하세요. 다음 층 두께로 Mo/Si 다층 거울(MLM) 구조를 정의합니다: Mo (1.950 nm), Mo-on-Si (0.806 nm), Si (3.843 nm), Si-on-Mo (0.386 nm)15.
  2. 표면 보호 재료를 평가하세요. Mo/Si 표면은 산화와 카바이드 형성에 취약하여 시간이 지남에 따라 광학 성능을 저하시키므로, 산화 및 카바이드 저항성을 평가하기 위해 Ru, RuO2, ZrO2, TiO2 코팅을 포함합니다.
  3. MLM 반사율을 계산하세요. 굴절률 데이터를 사용하여 Ru 캡핑층을 가진 Mo/Si 다층의 반사율을 평가하여 보호와 광학 효율 간의 절충관계를 정량적으로 평가할 수 있습니다.
    figure-protocol-2
    figure-protocol-3
    figure-protocol-4
    참고: 다양한 재료에 대한 δβ 값은 로렌스 버클리 국립연구소17 X-선 광학 센터에서 확인할 수 있습니다.
  4. 루 캡핑층에 대한 MLM 반사율: 굴절률을 사용하여 캡핑층 두께에 따른 반사율 변화를 계산합니다. 결과를 비교하여 광학 효율과 내구성 간의 트레이드오프를 파악합니다(그림 2).
  5. 출력 및 재현성 체크포인트: 그림 2 또는 Liu 등이 보고한 참조값에 따라 13.5 nm에서 반사율-두께 곡선을 생성하여 이 섹션의 성공 수행을 확인한다.

figure-protocol-5
그림 2. Ru 캡핑층의 두께가 다양한 Mo/Si 다층의 반사율. 이 그림의 더 큰 버전을 보시려면 여기를 클릭해 주세요.

2. 스퍼터 수율 계산

  1. 야마무라 공식을 적용하세요. Yamamura 등이 제안한 공식을 사용하여 스퍼터 수율(Y)을 계산합니다.18figure-protocol-6
  2. 정지 단면적을 계산하세요. 핵(Sn) 및 전자(Se) 정지 단면적을 방정식을 사용하여 평가합니다. (3)–(4).
    figure-protocol-7
    그리고figure-protocol-8
  3. 상수를 결정하세요. 경험적 상수 K를 식(5)을 사용하여 계산합니다.
    figure-protocol-9
    여기서 Z1 과 Z2는 각각 입사한 투사체와 표적 물질의 원자 번호를 나타냅니다; M1과 M2는 각각 입사 탄환과 표적 물질의 질량을 나타냅니다. ErEth 는 각각 환원 에너지와 임계 에너지이며, Es 는 표적 물질18의 표면 결합 에너지입니다.
  4. 실행 단계: 그림 3에 표시된 Python 스크립트를 실행하여 스퍼터 수율을 계산합니다. 그림 4에 보이는 Python 스크립트를 사용하여 야마무라 공식을 구현합니다. 컴퓨터에 Python 3와 NumPy 라이브러리가 장착되어 있는지 확인하세요. 그림 3 에 보이는 Python 스크립트를 실행하면 계산된 스퍼터 산출량을 담은 2열 텍스트 yield.dat파일이 생성되며, 이는 그림 5에 나타난 것입니다.
  5. 재현성 체크포인트: Ru에 영향을 미치는 Sn 이온에 대해 스퍼터링 수율 대 입사 에너지 곡선을 생성하여 이 섹션의 성공적인 실행을 확인한다(그림 5). Ru에 대한 Ar의 계산된 스퍼터 수율이 발표된 실험 데이터와 ±30% 이내로 일치하는지 확인하여 보정 검증 역할을 합니다.

figure-protocol-10
그림 3. 스퍼터 수율 계산을 위한 파이썬 스크립트입니다. 이 그림의 더 큰 버전을 보시려면 여기를 클릭해 주세요.

figure-protocol-11
그림 4. 야마무라 공식용 파이썬 스크립트. 이 그림의 더 큰 버전을 보시려면 여기를 클릭해 주세요.

figure-protocol-12
그림 5. Ru에서 Ar, Ru에서 Sn이 계산된 퍼터 수확량입니다. 왼쪽: 루; 오른쪽: Sn in Ru. Yamamura 등의 2.1단계에서 설명한 공식이 사용되었다. 현재 시뮬레이션과 Wu 등(26 ), Laegreid 등(등)의 시뮬레이션 간 비교가 수행되었다. 이 그림의 더 큰 버전을 보시려면 여기를 클릭해 주세요.

3. 이식 깊이 시뮬레이션

  1. 잠재적 모델을 선택하세요. RustBCA 코드19 의 KrC 퍼텐셜을 이온-고체 상호작용에 사용하세요:figure-protocol-13
  2. 스크리닝 기능을 정의하세요. Φ(r/a)를 지수 항들의 합으로 구현하라:
    figure-protocol-14
    1. 1의 다른 매개변수 cidi를 사용하여 KrC 퍼텐셜에 대해 a의 값을 다음 식과 같이 표현하라.figure-protocol-15
  3. 실행 단계: 그림 6에 표시된 Python 스크립트를 실행하여 RustBCA 실행 명령어가 스크립트에 통합되어 있어 이식 깊이를 계산합니다:
    1. 명령어 = "cargo run --release 1D "+ 입력파일 입력
    2. 그다음 os.system(command)를 입력하세요
  4. 그림 6에 표시된 Python 스크립트를 열고 스크립트에 따라 매개변수를 설정한 후 실행하면 계산된 이식 깊이가 포함된 depth.dat라는 두 열 텍스트 파일을 얻습니다.
  5. 재현성 체크포인트: 평균 이식 깊이 Sn을 생성하여 이 섹션의 성공적인 실행을 확인합니다(그림 7).
C1C2C3d1d2D3
0.190950.473670.335380.278540.637171.91925

표 1: KrC 퍼텐셜에 관여하는 매개변수 cidi.

figure-protocol-16
그림 6. 이식 깊이 계산을 위한 파이썬 스크립트입니다. 이 그림의 더 큰 버전을 보시려면 여기를 클릭해 주세요.

figure-protocol-17
그림 7. Ru-Mo-Si 다층 거울에 Sn 이온의 주입 깊이를 계산했습니다. 왼쪽: 2.0 keV(노란색)와 3.0 keV(파란색) 두 입사 에너지에서 10,000개의 입사 Sn 이온의 주입 깊이 분포; 오른쪽: 프로토콜 3.1단계에서 설명된 RustBCA에서 구현된 KrC 전위로 계산된 평균 Sn의 이식 깊이. 이 그림의 더 큰 버전을 보시려면 여기를 클릭해 주세요.

4. 제동력 계산

  1. 수소를 완충 가스로 모델링하세요. keV Sn 이온이 MLM에 미치는 손상을 줄이기 위해 완충 가스로 수소를 도입하세요.
    참고: 따라서 수소 및 MLM 표면 존재 시 keV Sn 이온의 정지력과 스퍼터링은 여전히 중요한 문제로 남아 있습니다.
  2. DFT 기반 퍼텐셜을 사용하세요. 수소-금속 시스템의 계산된 원자 간 퍼텐셜을 지글러–비어삭–리트마크(ZBL)와 모스 퍼텐셜 형식 모두에 적합시킨다.
    참고: 최근 연구20에서는 밀도 함수 이론(DFT) 계산에 기반한 수소-금속 시스템의 원자 간 퍼텐셜이 개발되었습니다.
  3. 재현성 체크포인트: SRIM 시뮬레이션 및 발표된 실험 데이터셋에서 얻은 참조 데이터와 에너지 의존적 정지 곡선을 비교하여 수소 내 Sn 이온의 계산된 차단력을 검증합니다.
    참고: 이 데이터는 Feng 등의 그림 6과 비교해야 합니다.
  4. 섹션 1–4의 출력(MLM 반사율, 스퍼터 수율, 이식 깊이, 정착 파워)을 결합하여 Sn 이온 노출 하에서 Mo/Si 다층 거울의 상대적 수명을 추정합니다.
    참고: 표면 거칠기 진화, 거울 기하학, 광선 추적과 같은 효과는 현재 프로토콜에 포함되지 않으며 향후 확장에 포함되어야 합니다.
  5. 동일한 작업 흐름을 Blue-X 리소그래피와 같은 대체 파장 영역에도 적용하여 광학 상수와 이온 에너지 분포를 조정합니다.

5. SnH4 형성 및 분해

참고: SnH4의 형성 및 분해에 대한 상세한 동역학적 연구는 여러 단면과 반응 속도 간의 Sn-H 간의 측정이 필요합니다. 이전에는 스탄인21의 전자 충격 이온화와 파편화, 반응률 XH4+H→XH3+H2,SnH 4+SnH→Sn2H3+H2, SnH4+SnH→Sn2H522,23이 보고된 바 있습니다. 그러나 SnH4의 플라즈마상 형성과 다양한 물질과의 상호작용 및 반응 메커니즘은 아직 완전히 규명되거나 이해되지 않았습니다. 따라서 스탄난 화학과 관련 분해 경로에 대한 실험 연구는 여전히 부족해12,24 추가 연구의 필요성을 강조한다.

  1. DFT 및 TST 계산: 밀도 함수 이론(DFT)과 가우시안 16에 구현된 전이 상태 이론(TST)을 결합하여 놓친 반응률을 계산합니다.
    참고: 이러한 계산 접근법은 반응 에너지학, 전이 상태, 속도 상수를 계산할 수 있게 하여 플라즈마 조건에서 스탄난 형성에 대한 상세한 메커니즘적 이해를 제공합니다.
  2. 반응 경로를 정의하세요. 여기에는 SnH4 형성으로 이어지는 두 개의 연속된 반응 경로가 포함되어 있습니다.
    (1) Sn+H2→SnH2
    (2) SnH2+H2→SnH4
  3. DFT 및 TST 계산을 수행하세요. 두 반응의 반응 에너지, 전이 상태, 속도 상수(k)를 계산하며, 결과는 그림 8그림 9에 나와 있습니다. 반응 열역학은 표 2표 4 에, 아레니우스 매개변수는 표 3표 5에 요약하라.
  4. 출력 및 재현성 체크포인트: 그림 8그림 9에 표시된 온도 의존 속도 곡선, 또는 보고된 값22,23을 재현하여 계산된 반응률 상수를 검증합니다.
  5. 검증된 속도 상수를 표 또는 기계 판독 가능한 형식(예: CSV 또는 TXT)으로 내보내 Sn–H 플라즈마 화학의 후속 동역학 모델링에 입력 매개변수로 직접 사용합니다.

figure-protocol-18
그림 8. Sn+H2→SnH2의 반응 속도와 에너지 장벽. 왼쪽: 반응률 상수 Sn+H2→SnH2; 오른쪽: 반응 경로를 위한 에너지 장벽(모든 회색 원자는 수소를, 파란색 원자는 Sn을 나타냅니다). 계산은 가우시안 16으로 수행됩니다. 이 그림의 더 큰 버전을 보시려면 여기를 클릭해 주세요.

반응제품ΔHΔGΔE
Sn+H2→SnH2SnH2-24.71-19.1317.87

표 2: 298.15 K 및 1 기압에서 세 반응 채널의 반응 엔탈피(H), 깁스 자유 에너지(G), 퍼텐셜 장벽(E)(kcal/mol).

아레니우스 매개변수방법반응
Sn+H2→SnH2
ATST2.50×10-13
TST/위그너1.13×10-13
TST/에카르트1.45×10-29
nTST0.85
TST/위그너0.93
TST/에카르트5.56
Ea(kJ/mol)TST68.99
TST/위그너65.3
TST/에카르트30.4
k(298K)(cm3 mol-1 sec-1)TST2.72×10-23
TST/위그너8.94×10-23
TST/에카르트1.03×10-21

표 3: 180에서 2000 K 온도 범위 내 Sn+H2→SnH2 반응의 아레니우스 매개변수.

figure-protocol-19
그림 9. SnH2+H2→SnH4의 반응 속도와 에너지 장벽. 왼쪽: 반응 속도 상수 SnH2+H2→SnH4; 오른쪽: 반응 경로를 위한 에너지 장벽(모든 회색 원자는 수소를, 파란색 원자는 Sn을 나타냅니다). 계산은 가우시안 16으로 수행됩니다. 이 그림의 더 큰 버전을 보시려면 여기를 클릭해 주세요.

반응제품ΔHΔGΔE
SnH2+H2→SnH4SnH4-26.5-32.8126.26

표 4: 298.15 K 및 1 기압에서 세 반응 채널의 반응 엔탈피(H), 깁스 자유 에너지(G), 퍼텐셜 장벽(E)(kcal/mol).

아레니우스 매개변수방법반응
SnH2+H2→SnH4
ATST3.73×10-17
TST/위그너1.23×10-17
TST/에카르트1.29×10-37
nTST1.55
TST/위그너1.67
TST/에카르트7.5
Ea(kJ/mol)TST136.39
TST/위그너132.94
TST/에카르트90.83
k(298K)(cm3 mol-1 sec-1)TST3.39×10-37
TST/위그너9.33×10-37
TST/에카르트6.56×10-36

표 5: 180에서 2,000 K 온도 범위 내 SnH2+H2→SnH4 반응의 아레니우스 매개변수.

6. 전자 에너지 분포 함수(EEDF) 계산

참고: 볼츠만 수송 방정식

이온화된 기체 내 전자 앙상블에 대한 볼츠만 방정식은 다음과 같습니다.

figure-protocol-20

여기서 f 는 6차원 위상 공간 내 전자 분포, v 는 속도 좌표, e 는 기본 전하, m 은 전자 질량(9.10956 × 10-31kg ), E 는 전기장, figure-protocol-21 는 속도-기울기 연산자, C 는 충돌로 인한 f 의 변화율을 나타냅니다.

  1. BOLSIG+ 해법을 실행하여 수소 플라즈마25의 볼츠만 수송 방정식을 풀기 위해 2항 근사법을 실행하세요.
  2. 실행 단계: BOLSIG+는 그래픽 창입니다.
    1. 그림 10A에 표시된 충돌 기 버튼을 클릭하면 H2의 단면적 데이터를 읽을 수 있습니다.
    2. 그림 10B에 표시된 "conditions" 파일에서 계산 매개변수를 선택하세요.
    3. 마지막으로, 그림 10C에 나타난 대로 EEDF 도면 버튼을 클릭하여 EEDF 이미지를 그립니다.
  3. 출력 및 재현성 체크포인트: 지정된 축소 전기장(E/N) 범위 내 수소 플라즈마에 대한 전자 에너지 분포 함수(EEDF)를 생성하여 BOLSIG+ 해법기의 성공적인 실행을 확인합니다. 그림 11로 EEDF를 검증하세요.
  4. 최종 EEDF 데이터를 표 형태(예: ASCII 또는 CSV 형식)로 내보내 Sn–H 플라즈마 화학의 동역학 모델링에 직접 입력할 수 있습니다.

figure-protocol-22
그림 10. BOLSIG+ 소프트웨어의 그래픽 인터페이스. 이 그림의 더 큰 버전을 보시려면 여기를 클릭해 주세요.

7. Sn–H 플라즈마 화학의 운동 운동 모델링

  1. PIC 시뮬레이션에서 플라즈마 파라미터를 가져오세요. 유체 시뮬레이션에서 전자 밀도와 플라즈마 온도를 포함한 플라즈마 매개변수를 추출합니다. 이 매개변수들을 PIC 시뮬레이션의 초기 조건으로 사용하여 Sn 이온의 시공간 분포와 에너지 스펙트럼을 얻습니다.
  2. 운동 시뮬레이션을 수행하세요. PIC에서 유도된 이온 에너지 분포와 DFT/TST에서 유래한 반응 속도를 입력으로 사용하여 Sn, SnHx 및 관련 중간체의 결합 속도 방정식을 풀어. EUV 소스 작동과 관련된 수소 플라즈마 조건 하에서 종 밀도의 시간적 변화를 추적합니다.
  3. 운동 출력과 표면 상호작용 모델을 결합하세요. 2–4절에서 얻은 제동력, 스퍼터 수율, 이식 깊이 분포와 함께 동역학적 결과를 결합합니다. 이러한 결합된 출력을 사용하여 분해 메커니즘을 평가하고 Mo/Si MLM의 유효 수명을 추정합니다.

Access restricted. Please log in or start a trial to view this content.

Results

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

스퍼터링 수율 보정 및 검증
보정 단계로서 Ru 내 Ar 원자의 스퍼터 생성율을 계산하세요. 이 스퍼터 수율은 프로토콜 단계 2.1(야마무라 모델)의 출력을 나타냅니다. 결과는 그림 5 (왼쪽)에 나와 있습니다. Wu 등26 과 Laegreid 등27 이 보고한 실험 데이터는 대체로 일관성이 있다. 본 모델의 이론적 결과는 80 eV 이상의 입사 에너지에서의 실험 측정과 좋은 일치를 보입니다. 저입사 에너지(<40 eV)에서는 Wu 등의 데이터와의 상대적 편차가 약 60%에 달한다. 이 비교는 프로토콜 2.4의 올바른 실행을 확인하는 검증 체크포인트 역할을 합니다.

이 보정을 바탕으로 Ru에 입사한 Sn 이온의 스퍼터 수율이 그림

Access restricted. Please log in or start a trial to view this content.

Discussion

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

볼츠만 수송 방정식(BTE), 입자-인-셀(PIC), 운동 시뮬레이션을 결합한 통합 방법론은 극자외선(EUV) 리소그래피에서 주석(Sn) 잔해 저감을 조사하는 통합된 틀을 구축합니다. 구체적으로, 유체 시뮬레이션은 플라즈마 매개변수—밀도와 온도—를 산출하며, 이를 PIC 프로그램에 통합하여 SnxHy 분자의 시공간적 분포를 얻을 수 있습니다. 이 PIC 결과와 DFT/TST 계산에서 얻은 반응률을 결합함으로써 Sn-H 종의 완전한 동역학적 조사가 가능하다. 이 다중 규모 접근법은 플라즈마 속도론, 이온 수송, 표면 반응의 동시 모델링을 가능하게 하며, 이는 원천 효율성과 내구성 향상에 필수적입니다. 정확한 단면과 반응 속도 수집에 대한 최근 진전은 Sn–수소 시스템의 현실적인 동역학 시뮬레이션에 견고한 토대를 제공합니다.

여러 대표적인 결과가 이 프로토콜의 신뢰성을 확인시켜 줍니다...

Access restricted. Please log in or start a trial to view this content.

Disclosures

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

저자들은 공개할 이해 상충이 없습니다.

Acknowledgements

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

중국국가자연과학기금(National Natural Science Foundation of China)의 지원금(Grant No.12374231)에 감사드립니다.

Access restricted. Please log in or start a trial to view this content.

Materials

List of materials used in this article
NameCompanyCatalog NumberComments
볼시그+플라즈마 및 에너지 변환 실험실, 폴 사바티에 대학교2025년 4월 24일에 업데이트된 버전입니다
가우시안가우시안 주식회사가우시안 16
러스트BCA일리노이 대학교 어배너-샴페인 캠퍼스 원자력, 플라즈마 및 방사선공학과1.2.0

References

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,
  1. O’Sullivan, G., et al. Spectroscopy of highly charged ions and its relevance to EUV and soft X-ray source development. J Phys B At Mol Opt Phys. 48, 144025(2015).
  2. Versolato, O. O. Physics of laser-driven tin plasma sources of EUV radiation for nanolithography. Plasma Sources Sci Technol. 28, 083001(2019).
  3. Bayerle, A., et al. Sn ion energy distributions of ns- and ps-laser produced plasmas. Plasma Sources Sci Technol. 27 (4), 045001(2018).
  4. Rai, S., et al. Evidence of production of keV Sn+ ions in the H2 buffer gas surrounding an Sn-plasma EUV source. Plasma Sources Sci Technol. 32 (3), 035006(2023).
  5. Spatial separation of EUV emission and energetic ions by use of double-laser-pulse irradiation. Sugiura, T., et al. Proc SPIE 13177 Photomask Japan XXX Symposium on Photomask and Next-Generation Lithography Mask Technology, , 1317715(2024).
  6. Fleur, V. Energy loss and scattering of energetic Sn ions interacting with H2: prospects of time-of-flight investigations. [PhD thesis]. , University of Groningen. (2021).
  7. Rai, S. Ionic interactions around EUV generating tin plasma. [PhD thesis]. , University of Groningen. (2023).
  8. Mertens, B., et al. Progress in EUV optics lifetime expectations. Microelectron Eng. 73-74, 16-22 (2004).
  9. Ugur, D., Storm, A. J., Verberk, R., Brouwer, J. C., Sloof, W. G. Decomposition of SnH4 molecules on metal and metal–oxide surfaces. Appl Surf Sci. 288, 673-676 (2014).
  10. Elg, D. T., et al. Removal of tin from extreme ultraviolet collector optics by in-situ hydrogen plasma etching. Plasma Chem Plasma Process. 38, 223-245 (2018).
  11. van Herpen, M. M. J. W., Klunder, D. J. W., Soer, W. A., Moors, R., Banine, V. Sn etching with hydrogen radicals to clean EUV optics. Chem Phys Lett. 484 (4-6), 197-199 (2010).
  12. Garza, R., et al. Stannane in extreme ultraviolet lithography and vacuum technology: synthesis and characterization. J Vac Sci Technol A. 41 (6), 063209(2023).
  13. Biggerstaff, S., et al. Comparative study of neutral and cationic Sn2H2: toward laboratory detection of the cation. J Phys Chem A. 128, 7090-7104 (2024).
  14. Xiao, Z., et al. Multiple states and roles of hydrogen in p-type SnS semiconductors. Phys Chem Chem Phys. 20, 20952-20960 (2018).
  15. Liu, X., et al. Comparative study on microstructure of Mo/Si multilayers deposited on large curved mirror with and without the shadow mask. Micromachines (Basel). 14 (3), 526(2023).
  16. Yao, D., et al. Fabrication and characterization of TiO2 and SiO2 as protective coating for Mo/Si multilayer by ion beam sputtering. Vacuum. 238, 114287(2025).
  17. Henke, B. L., Gullikson, E. M., Davis, J. C. X-ray interactions: photoabsorption, scattering, transmission, and reflection at E = 50-30000 eV, Z = 1-92. At Data Nucl Data Tables. 54, 181-342 (1993).
  18. Yamamura, Y., Tawara, H. Energy dependence of ion-induced sputtering yields from monatomic solids at normal incidence. At Data Nucl Data Tables. 62, 149-253 (1996).
  19. Drobny, J. T., Curreli, D. RustBCA: a high-performance binary-collision-approximation code for ion-material interactions. J Open Source Softw. 6, 3298(2021).
  20. Feng, X., Song, Y., Ma, Y., Li, B. DFT-based interatomic potentials for hydrogen-metal systems: improved stopping power modeling. Nucl Instrum Methods Phys Res B. 572, 165999(2026).
  21. Song, Y., Ma, Y., Li, B., Chen, X. Decomposition of electron ionization mass spectra and calculation of electron ionization cross sections of SnxHy for extreme ultraviolet lithography. Phys Scr. 100, 045405(2025).
  22. Ma, Y., Li, B. A comparative study of kinetic and thermodynamic mechanisms of XH4 + H → XH3 + H2 reaction (X = Si, Ge, and Sn). AIP Adv. 15, 075138(2025).
  23. Ma, Y., Li, B. Reaction pathways between SnH4 and SnH relevant to EUV lithography: a DFT and TST study. Plasma Chem Plasma Process. 46, 40(2026).
  24. Rieger, J., Benter, T., Kersten, H. High-resolution electron ionization mass spectrometry of stannane: deconvolution of superimposed fragmentation patterns. J Am Soc Mass Spectrom. 35, 1523-1532 (2024).
  25. Hagelaar, G. J. M., Pitchford, L. C. Solving the Boltzmann equation to obtain electron transport coefficients and rate coefficients for fluid models. Plasma Sources Sci Technol. 14, 722-733 (2005).
  26. Wu, S. Sputtering yields of Ru, Mo, and Si under low energy. J Appl Phys. 106, 054902(2009).
  27. Laegreid, N., Wehner, G. K. Sputtering yields of metals for Ar+ and Ne+ ions with energies from 50 to 600 eV. J Appl Phys. 32, 365-369 (1961).

Access restricted. Please log in or start a trial to view this content.

Reprints and Permissions

Request permission to reuse the text or figures of this JoVE article

Request Permission

Tags

EUV LithographyTin Debris ControlBoltzmann Transport EquationParticle In Cell SimulationKinetic ModelingHydrogen PlasmaSputtering YieldImplantation DepthMLM ReflectivityDensity Functional Theory

Related Articles