검증 및 문헌 근거 (V&V)

이 시뮬레이터는 RF/CCP 플라즈마의 빠른 상대 트렌드 분석을 위한 축소 유체 공식을 사용합니다. 플라즈마 밀도, 이온 플럭스, RF 임피던스, 시스 응답 및 파라미터 민감도에서 정성적·반정량적 트렌드를 추정합니다. 완전 검증된 PIC/MCC, 상용 다유체 시뮬레이션 또는 실험을 대체하지 않습니다.

• 레시피 트렌드 비교 (같은 가스 조건에서 A vs B)
• 압력·전력·간격·주파수 민감도 분석
• 플라즈마 밀도 상대 트렌드 추정
• 이온 플럭스 상대 트렌드 추정
• RF 임피던스 및 매칭 경향 탐색
• 교육용 시각화 및 물리적 직관 개발
• 초기 공정 직관 형성 (설계 공간 선별)

아래 용도에는 이 시뮬레이터를 사용하지 마십시오. 결과가 오해를 불러일으킬 수 있습니다.

• 보정 없는 절대 밀도 수치
• 최종 공정 자격 또는 양산 적합성 판단
• 상세 래디컬 화학 반응 분석
• 저압 비국소 운동 분석
• 확률론적 시스 가열 분석
• 고주파 전자기 파동 효과 분석
• 하드웨어 설계 최종 확인
• 실험·PIC-MCC·검증된 상용 솔버의 대체

보정 없는 전자 밀도 불확도는 약 2–5배입니다. 이온 플럭스 절대값은 ~1 오더 수준입니다. 실험 데이터로 보정하면 특정 운전 조건 근처의 상대적 트렌드 정확도는 향상되지만, 절대 정확도는 시스템 의존적으로 유지됩니다. 더 넓은 파라미터 범위에 걸친 절대적 신뢰도를 요구하는 용도에는 전체 PIC/MCC 또는 검증된 다유체 상용 솔버를 사용하십시오.

아래 업그레이드는 두 그룹으로 나뉩니다. (a) 항상 적용 — 별도 설정 없이 모든 실행에 반영됩니다. (b) 플래그 게이트 — 환경변수로 제어하며, ION_SOLVER·RF_SHEATH·TE_PDE는 기본 ON("0"으로 비활성), GUMMEL만 기본 OFF(OFF일 때 bit-exact)입니다. 이 모델은 여전히 추세 수준 환산 유체 모델이며 절대 정확도를 보장하지 않습니다.

(a) 항상 적용 (기본 활성화)

(b) 플래그 게이트 (ION_SOLVER·RF_SHEATH·TE_PDE 기본 ON, GUMMEL만 기본 OFF)

ION_SOLVER·RF_SHEATH·TE_PDE는 이제 기본 ON이며 환경변수를 "0"으로 설정해 끌 수 있습니다. GUMMEL 외부 루프만 여전히 기본 OFF(opt-in)이고, OFF일 때 결과는 기존 파이프라인과 bit-exact입니다.

1. 공지된 GEC 셀 조건(Hargis et al. 1994)에서 Ar 기준 시뮬레이션을 실행합니다.
2. 계산된 전자 밀도와 이온 플럭스 트렌드를 동일 조건에서 보고된 PIC/MCC 또는 유체 벤치마크 결과와 비교합니다.
3. 전압·압력·간격 스윕에서 트렌드 방향(증가/감소)이 일치하는지 확인합니다.
4. 가능하면 실험 데이터(Langmuir 프로브, VI 프로브, 매칭)로 보정 인수를 도출합니다.
5. 보정 후 가장 가까운 운전점 근처에서만 절대 수치를 신뢰하십시오.
6. 결론 도출 전 반드시 범위와 한계 섹션을 검토하십시오.

• GEC 참조 셀 (Hargis et al. 1994) — 잘 알려진 평행판 13.56 MHz Ar CCP 기준 실험.
• Ar CCP 전자 밀도 및 이온 플럭스 — 문헌 값과 직접 비교 (밀도 ~10⁹–10¹¹ cm⁻³, 이온 플럭스 ~10¹⁵–10¹⁶ m⁻²s⁻¹).
• 압력·전력·간격·주파수 스윕 트렌드 — 방향성 동의 여부를 Boeuf & Pitchford (1995) 및 Surendra & Dalvie (1993) 결과와 비교.
• VI 프로브 및 매칭 네트워크 임피던스 — RF 임피던스 트렌드를 실험실 측정과 비교.
• COMSOL Plasma Module — 검증된 상용 다유체 솔버와의 트렌드 비교 기준.
• Zapdos/MOOSE — MOOSE 기반 오픈소스 플라즈마 유한요소 코드로 독립 검토 가능.

아래 표는 각 모델 출력을 공개 문헌 값에 앵커한 검증 체크입니다. 각 항목은 프리셋 딥링크(클릭하면 시뮬레이터에 바로 로드), 읽을 모델 출력, 문헌 근거, 통과 기준으로 구성됩니다. 이는 절대값이 아니라 오더/방향 수준의 검증입니다.

명시적 제외 (문헌 검증 대상 아님)

아래 항목은 의도적으로 검증 대상에서 제외됩니다 (physics-model.md §12 반영).

• 공정 결과 전반 — 식각률·증착률·선택비·프로파일 (라디칼 화학·표면 반응 미모델; 플라즈마 물리량만 보고하며 상대 프록시도 제공하지 않음)
• EEDF / 비-Maxwellian 상세
• space-charge 시스 전위 프로파일 (의도적 기각 — 미해상 λ_De)
• 절대 IEDF eV 스케일 (전력흡수·|E|는 이제 물리 단위 W/m³·V/m, 요청 전력에 파워-밸런스로 정규화)
• IEDF 이중 피크 구조 (RF 사이클 미분해)
• 확률론적 시스 가열 (Drude 옴 가열만)
• 절대 임피던스 (옴)
• 완전 CFD 중성 유동 (이제 showerhead→pump 포텐셜 유동 — 면-선택적 방출 + 제트 컬럼 방향 — 이지만 Navier–Stokes CFD는 아님)

아래 매트릭스는 대표 벤치마크 케이스별로 기대 트렌드와 판정 기준을 정리합니다. 이는 방향성·크기 규모(비율·오더) 검증이며 절대값 정확도가 아닙니다. 각 판정은 방향/비율 경계로 표현됩니다. B13(문헌 앵커 V1–V10)과 상호 보완적입니다. 상태 칩: 검증됨(초록)=판정 기준 전체가 tests/test_validation_matrix.py의 자동 테스트로 확인됨, 부분(파랑)=방향성은 CI로 확인되나 기준의 일부 조항은 미확인, 예상(중립)=문헌 근거는 있으나 회귀 체크 미고정, 예정(주황)=커버리지 없음.

명시적 범위 밖 (일치하지 않음) 목록은 위 우리가 주장하지 않는 것 →에 정리되어 있습니다.

아래 문헌은 모델 구성 요소별로 그룹화되었습니다. 이미 위 섹션에서 인용된 문헌은 재인용하지 않습니다.

드리프트-확산 / 유체 모델 (일반)

• Lieberman & Lichtenberg, Principles of Plasma Discharges and Materials Processing, 2nd ed., Wiley 2005 — drift-diffusion, global model, sheath (이미 섹션 1에 인용)
• Boeuf & Pitchford, Phys. Rev. E 51, 1376 (1995) — RF fluid CCP model, Poisson-coupled fluid + power absorption (이미 섹션 1에 인용)
• Scharfetter & Gummel, IEEE Trans. Electron Dev. 16, 64 (1969) — exponential-fit flux discretization (이미 섹션 3에 인용)

전자 속도 계수 LUT / EEDF (A1)

• Hagelaar & Pitchford, Plasma Sources Sci. Technol. 14, 722 (2005) — BOLSIG+ Boltzmann solver for electron transport coefficients and rate coefficients
• BOLSIG+ official site (LAPLACE, Toulouse) — free electron Boltzmann solver for generating LUT data
• LXCat cross-section database — source data for A1 LUT / BOLSIG+ inputs (electron-neutral cross sections)

0D 전자 에너지/입자 수지 (A2)

• Lee & Lieberman, J. Vac. Sci. Technol. A 13, 368 (1995) — global model 0D energy and particle balance for high-density plasmas

시스 / Bohm 기준 / RF 시스 (A3)

• Riemann, J. Phys. D 24, 493 (1991) — Bohm criterion and sheath formation (이미 섹션 4에 인용)
• Lieberman, IEEE Trans. Plasma Sci. 16, 638 (1988) — analytical RF sheath model (이미 섹션 4에 인용) [verify IEEE Xplore]
• Godyak & Sternberg, Phys. Rev. A 42, 2299 (1990) — dynamic model of electrode sheaths in symmetrically driven RF discharges

GEC 벤치마크 셀 (B8)

• Hargis et al., Rev. Sci. Instrum. 65, 140 (1994) — GEC RF reference cell: defined parallel-plate 13.56 MHz system for benchmark studies

비교 / 오픈소스 솔버 (B8)

• COMSOL Plasma Module — validated commercial multi-fluid plasma solver for benchmark comparison
• Zapdos (shannon-lab / MOOSE) — free and open-source MOOSE-based finite-element plasma simulation code

자기장(B) 입력이 활성화되면 시뮬레이터는 트렌드 수준의 자화 전자 구속 모델을 적용합니다. 핵심 원리: B 자기장이 전자를 자화하여 Hall 파라미터 β = ω_ce/ν_m을 형성합니다 (ω_ce = eB/m_e: 전자 사이클로트론 주파수, ν_m: 전자-중성 충돌 주파수). β가 크면 수직 방향(cross-field) 전자 수송이 1/(1+β²) 인수만큼 감소하며, 이를 방향에 따라 aniso_r, aniso_z로 적용합니다. 그 결과 전자 구속이 강화되고 밀도가 상승하는 경향(mag_gain)이 나타납니다.

이 자기장 모델은 다음을 포함하지 않습니다: Hall 텐서 기반 완전 이방성 수송 / ExB 드리프트 회전 / 전자 사이클로트론 공명(ECR) 가열 / 자화된 시스 물리. 방향별 Hall 파라미터 이방성 수송 + 경량 ExB 순수 트렌드 이득에 기반한 트렌드 지향 근사이며, 절대 자화 플라즈마 결과로 사용하지 마십시오.

• Lieberman & Lichtenberg 2005, Ch. 4 — electron transport in magnetic fields, Hall parameter
• Raizer 1991, Ch. 9 — magnetized electron transport and drift-diffusion in magnetic fields

아래 표는 이 축소 유체 모델이 파라미터 스윕을 통해 검증한 상대적 트렌드 방향을 정리한 것입니다. 화살표(↑ ↓ ~)는 출력이 움직이는 방향을 나타내며, 보정되지 않은 절대 크기가 아닙니다. 절대 사용을 위해서는 실험 데이터로 보정이 필요합니다.

↑ = 증가, ↓ = 감소, ~ = 거의 변화 없음. 크기는 트렌드 수준이며 절대 사용에는 실험 보정이 필요합니다.

Te(r,z) 맵 (fields.te_eV)

fields.te_eV 히트맵은 기본적으로 휴리스틱 에너지 맵을 보정된 스칼라 Te에 평균-정합한 것입니다 (휴리스틱과 0D 입자 밸런스의 기하 블렌드, [1,9] eV 클램프). PLASMACCP_TE_PDE를 켜면 공간 형상은 RF 흡수 맵으로 가열되는 screened-diffusion PDE에서 나옵니다.

nₑ 관측값 (fields.ne)

표시되는 fields.ne = v/(v+n_sat) ∈ [0,1] 포화 변환입니다 (요청마다 고정된 곡선이라 A/B 범례가 비교 가능하게 유지됨). 절대 밀도 스케일은 density_ref / 0D 블록에 있습니다.

전극 갭

갭 유도량(Te 갭 인자, Child-Langmuir 클램프, Λ=gap/π)은 격자에서 도출한 방전 갭 — 파워드 블록과 접지 블록 사이 플라즈마 셀 스팬의 중앙값 — 을 사용합니다. 도메인 높이 z_max는 폴백 전용입니다.

보정 노브 (calibration 블록)

중립 기본값(1.0 / 0.0) = 기존 결과와 bit-exact.

위 B1–B14가 무엇을 모델링하는지를 다뤘다면, 이 절은 그것을 실제로 어떻게 계산하는지를 정리합니다 — 격자, 이산화, 선형대수 솔버, 그리고 그 위에 얹힌 검증 게이트입니다.

수치적 정직성: 목표는 절대값이 아니라 트렌드입니다. CI가 매 커밋 고정하는 것 — 모든 벡터화 경로의 스칼라 대비 패리티, 모든 신규 플래그의 off=bit-exact 회귀, 방향성 어서션(파워↑→ne↑ 등), r-alternation(체커보드) 지표 가드, M1–M8 벤치마크(무수정 통과 필수).

bit-exact 탈출구: PLASMACCP_PLANAR_RADIAL_OPERATOR=1(구 planar 연산자 재현), PLASMACCP_TG_FEEDBACK=0, calibration 중립값(1.0/0.0), RF_SHEATH="0"(Child 폴백 단독). native C 어셈블러는 고종횡비 checkerboard M-matrix 버그에 더해 BC2 이후 원통 가중을 받지 못해 물리적으로도 구식(planar)이므로 prod 기본 비활성이며, ENABLE=1 단독 지정은 정합 가드가 무시하고 경고합니다(PLANAR_RADIAL_OPERATOR=1과 함께일 때만 유효).

V1 엔진에는 절대 전자 밀도를 구하는 두 가지 방식이 있습니다. 표준 모드(오늘의 기본값)는 보정된 경험적 밀도 스케일을 그대로 씁니다. 물리 모드(베타)는 대신 전극 사이의 벽 전력균형 — 흡수 전력을 이온화가 벽에서 잃는 에너지로 나눈 값 — 에서 밀도 수준을 직접 계산합니다. 둘 다 같은 2D 밀도 분포 솔버 위에서 돌아가며, 물리 모드는 그 분포의 절대 크기를 바꾸는 스위치입니다.

메커니즘

밀도는 n = P_abs / (e·u_B·(E_c(Te)+2Te+V_sheath)·A_eff) 형태의 면별 전력균형에서 풀립니다. 이온화 소스는 전력이 인가된 전극과 그 맞은편(위에 웨이퍼가 있으면 웨이퍼 뒤 도체 면까지 포함)의 구동 시스에 놓이고, 전자 에너지 완화 길이 λ_ε(Tsendin 2010)로 공간에 평활됩니다. 시스 전압은 같은 전력균형과 맞물려 자기일관적으로(반복 수렴) 다시 풀립니다. 마지막으로 Bohm 자기일관성 검사(R1)가 도는데, 계산된 벌크 밀도·Te에서 기대되는 이온 전류밀도와 실제로 보고된 이온 전류밀도를 각 채널(전력 전극·벽)마다 대조합니다 — 배율 κ=4를 벗어나면 그 정도에 비례해 로그 공간에서 경험적 스케일 쪽으로 되섞습니다.

R1은 자기일관성만 봅니다(그 채널 자신의 밀도·Te가 그 채널 자신의 전류와 맞는지) — 문헌의 절대 기준값과 대조하지 않으므로 절대 밀도가 맞다는 증거가 될 수 없습니다.

Te도 모든 기체·혼합물에서 입자균형으로 풉니다 — 3항 Lee–Lieberman h-factor와 기체별 이온 이동도(Ar·N2·O2·H2·He는 실측값, 나머지는 추정값이며 어느 쪽인지 실행 노트에 표시)를 씁니다. 화면에 보이는 Te는 밀도 계산에 실제로 쓰인 그 Te입니다. 다중 주파수 요청도 물리 밀도를 받습니다 — 여러 톤이 총 전력으로 한 플라즈마를 구동하고, 전력균형의 시스 전압은 톤마다 자기 전력으로 풀어 더합니다. 부착 손실은 실제로 부착하는 기체에만 걸립니다.

기체 계수(전리율 k_iz, 이온화 1건당 에너지 손실 E_c)는 기체마다 다른 근거를 씁니다: Ar·He·N2·O2·H2·SiH4·CF4·Cl2·N2O 9종은 전용 인용 fit(N2O는 Song et al. 2019 단면적을 직접 적분), NH3·SF6·CH4·C2H2는 NIST 공개 BEB 전리 단면 적분, Ne·Kr·Xe·NF3·C4F8·CHF3·HBr 7종은 화학적으로 유사한 기체의 전용 fit 형태에 그 기체 자신의 인용된 이온화 전위만 넣어 재조정한 값입니다(새로 맞춘 계수 아님). 등재되지 않은 기체는 실패하지 않고 아르곤 계수를 쓰되, 실행 노트에 그 사실이 명시됩니다. 혼합물의 에너지 손실 E_c는 성분별 값을 단순 몰분율이 아니라 각 성분이 실제로 기여하는 이온화율로 가중합니다.

다중 주파수 요청에서는 프로세스 모니터가 톤별 시스 전압·이온 전력을 따로 보여줍니다(2톤일 때 Lieberman의 이중주파 이온에너지식 적용).

결과 배지의 세 상태

검증 상태 — 정직하게

• 완주(넓은 범위): 26종의 기체·혼합물(신규 11종 포함)에 걸쳐 무작위 조합(LHC)과 단일 요인 스윕을 합친 950셀(표준·물리 각 475)이 전부 완주했습니다 — 예외·거부 0, 전 보고값 유한, 런타임은 표준의 3배 이내(중앙값 0.78배). 물리 밀도를 그대로 채택(위 표의 벽 전력균형)한 것이 85.9%, 자기일관성 검사(R1)가 경계 밖이라 경험 기준값으로 완전히 되돌아간 것이 14.1%입니다 — 이는 실패가 아니라 저하(degraded)로 집계되며, 대부분 극단적인 무작위 조합에서 나타났습니다. R1은 자기일관성 검사일 뿐이라 이 수치가 절대값이 맞다는 뜻은 아닙니다. (이전 라운드의 좁은 268셀 스윕 — 기체 10종 × 압력만 — 은 아래 절대값·모양·전력/압력 의존성 항목의 근거로 남아 있습니다.)
• 절대값 — 아직 검증되지 않음: 흡수 전력이 명시적으로 정의된 공개 기준인 Turner et al. (2013) He PIC 벤치마크와 대조했습니다. 간극을 Turner와 같게(67 mm) 맞춘 평판 챔버에서 중심 밀도는 물리 모드가 Turner 값의 0.85–4.8배(4개 압력 중 3개가 0.85–1.5배), 표준 모드가 0.06–0.62배였습니다. 작은 연구용 셀(간극 25 mm, 이전 측정)에서는 물리 모드 0.11–1.32배, 표준 모드 0.11–0.62배입니다. 참값은 두 모드 사이에 있는 것으로 보이며, 어느 쪽도 절대값이 검증되었다고 말할 수 없습니다.
• 분포 모양: 이전 라운드에서 물리 모드의 밀도 분포가 간극 중앙이 오목하게 파이는 비물리적 결함을 가진 적이 있었습니다(부피 손실항의 중복 계산). 이를 제거한 뒤로는 Turner의 4개 벤치마크 사례 전부에서 Turner Fig. 7–10과 같은 중앙-최대·평탄 상단 모양이 나옵니다.
• 전력 의존성: GEC 조건(26.7–53.3 Pa)에서 물리 모드의 전력 지수는 0.85–0.90으로, 실측(Land et al. 2009, ≈0.98–1.08)에 가깝습니다. 표준 모드는 0.086으로 이 전력 구간에서 사실상 반응하지 않습니다.
• 압력 의존성: Turner 벤치마크(30→1000 mTorr, ×18.4)에 대해 물리 모드는 ×6.6–14.5, 표준 모드는 ×2.6–3.3으로 물리 모드가 더 가깝지만 둘 다 완만합니다. GEC 조건(26.7→53.3 Pa)에서 두 모드 모두 압력과 함께 늘어납니다(물리 지수 약 +1.1). Land et al.(2009)의 측정은 줄어들지만, 그 전력은 위상을 뺀 피상 전력(V·I)이고 저자 스스로 이 결과가 문헌과 다르다고 적었습니다. 흡수 전력이 (거의) 고정된 공개 기준 — Turner 2013 PIC, Schulenberg et al. 2021, NIST GEC 특집(Hopkins·Overzet 1995) — 은 모두 증가를 보입니다. 가열 모드 전환(Godyak & Piejak 1990, 약 0.4–0.5 Torr)은 V1에 해당 물리가 없어 어느 모드에서도 나타나지 않습니다.

언제 쓰고 언제 쓰지 말아야 하는가

• 쓸 곳: 전력을 바꿀 때 밀도가 어떻게 반응하는지 보는 트렌드 탐색 — 물리 모드의 전력 지수(~0.9)가 표준 모드(~0.09)보다 실측에 훨씬 가깝습니다.
• 쓰지 말 곳: 절대 밀도 숫자가 필요한 어떤 용도든(공정 자격, 장비 스펙 비교) — 흡수 전력이 정의된 벤치마크에서 물리 모드는 조건에 따라 최대 약 5배 과대이고, 이온 에너지(시스 전압)는 크게 과소입니다. 압력 지수의 정확한 값과 가열 모드 전환 영역은 미검증입니다.
• 표준 모드는 여전히 기본값이고, 물리 모드는 opt-in 베타입니다. 배지에 뜨는 밀도 산출 방식(위 표)을 항상 함께 읽으십시오 — "혼합"이나 "경험 기준값"이 뜨면 그 결과의 절대 밀도는 물리 계산이 아니라 표준 모드와 같은 경로에서 온 것입니다.
• 알려진 한계 — 명시적으로: (1) 전기음성 기체(Cl2·CF4·SF6 등)는 여전히 양이온으로 풀립니다 — 음이온을 만드는 부착은 계수에 반영되지만, 음이온 자체를 별도 전하운반자로 추적하지는 않습니다. (2) 위 950셀 스윕에서 물리 밀도가 경험 기준값으로 완전히 되돌아가는 비율은 14.1%입니다 — 이 조건들에서는 물리 모드가 표준 모드와 다를 바 없는 숫자를 냅니다. (3) He와 NH3는 조회 가능한(라이선스 문제 없는) 실측 단면 자료를 구하지 못해, 전리율 이외의 항(He의 들뜸·탄성 항, NH3의 에너지 손실·부착)은 여전히 대용치이거나 다른 기체에서 빌려 온 값입니다.

• Turner, Hutchinson, Kushner et al., Phys. Plasmas 20, 013507 (2013) — arXiv:1211.5246 (open access); multi-code PIC-MCC benchmark, He, absorbed power explicitly defined
• Lieberman & Lichtenberg, Principles of Plasma Discharges and Materials Processing, 2nd ed., Wiley 2005 — wall power balance / global model (already cited above)
• Tsendin, Phys.-Usp. 53, 133 (2010) — electron energy relaxation length λ_ε, used to smooth the fixed ionization source

공정 요약(Process summary) 탭은 웨이퍼 수준의 트렌드 지표를 보여 줍니다. 각 지표가 무엇을 의미하는지, 공정 엔지니어가 왜 주목해야 하는지, 그리고 시뮬레이터가 어떻게 계산하는지 아래에 설명합니다.

반경 이온 플럭스 비균일도

이온 플럭스 반경 프로파일의 평탄도를 나타냅니다. 값이 높다는 것은 웨이퍼 중심과 가장자리에 도달하는 이온 플럭스 자체가 그만큼 다르다는 뜻입니다. 시뮬레이터는 반경 프로파일에서 비균일도(%) = (max − min) / (2 × mean) × 100으로 계산하며, 낮을수록 좋습니다. 제공된 웨이퍼 구간에 0 또는 비유한 표본이 있으면 누락된 표면 지지와 실제 0을 구분할 수 없으므로 숫자를 숨기고 사용 열/전체 열 coverage를 표시합니다. 실제 장비의 균일도는 가스 유량, 척 온도, 포커스 링 상태 등 여기서 모델링하지 않는 요소에도 크게 의존합니다.

• Lieberman & Lichtenberg, Principles of Plasma Discharges and Materials Processing, 2nd ed., Wiley 2005 — Ch. 10–11 (uniformity, transport)
• Boeuf & Pitchford, Phys. Rev. E 51, 1376 (1995) — 2D CCP fluid model with radial uniformity analysis

이온 에너지 / IEDF

기판에 입사하는 이온의 에너지는 시스 전압 강하로 결정되며, 이온 하나가 표면에 전달하는 운동에너지의 크기를 정합니다. 이온 에너지 분포 함수(IEDF)의 형태와 평균값이 그 분포를 기술합니다. 시뮬레이터는 Child-Langmuir 스케일링에 기반한 시간 평균 이온 에너지 프록시를 각 반경 위치에서 계산하며, 절대 에너지(eV)가 아닌 상대 트렌드를 나타냅니다. 이중 주파수 구성에서는 고주파(HF)가 이온 플럭스를, 저주파(LF)가 이온 에너지를 주로 제어합니다(Boyle et al. 2004). 실제 IEDF는 RF 시스의 위상 평균화로 인해 이중 피크 구조가 나타나며(Godyak et al. 1991), 이 시뮬레이터는 그 구조를 직접 분해하지는 않습니다.

• Lieberman & Lichtenberg 2005, Ch. 11 — RF sheath and IEDF theory
• Boyle, Ellingboe & Turner, J. Phys. D 37, 697 (2004) — HF controls ion flux, LF controls ion energy in dual-frequency CCP
• Godyak, Piejak & Alexandrovich, J. Appl. Phys. 69, 3455 (1991) — measured ion energy distributions and flux in RF discharges
• Lieberman, IEEE Trans. Plasma Sci. 16, 638 (1988) — analytical RF sheath model [verify IEEE Xplore] [verify]

시스 두께

시스(Sheath)는 전극 근처에 형성되는 양전하 공간 전하 영역으로, 이온이 기판을 향해 가속되는 구간입니다. 시스 두께는 인가 전압, 플라즈마 밀도, 주파수에 따라 결정되며(Child-Langmuir 법칙), 이온 충돌 에너지와 충돌 각도 분포에 직접 영향을 미칩니다. 시스가 두꺼울수록 이온 가속 거리가 길어지고 에너지가 높아지는 경향이 있습니다. 시뮬레이터는 Bohm 속도 조건과 Child-Langmuir 스케일링으로 각 반경 위치의 시스 두께를 추정하며, 이는 트렌드 신호로 해석해야 합니다.

• Lieberman, IEEE Trans. Plasma Sci. 16, 638 (1988) — analytical RF sheath model [verify IEEE Xplore] [verify]
• Riemann, J. Phys. D 24, 493 (1991) — Bohm criterion and sheath formation
• Lieberman & Lichtenberg 2005, Ch. 2 & 11 — sheath physics and Child-Langmuir scaling

이온 플럭스·에너지 지수

이온 플럭스와 이온 에너지라는 두 플라즈마 물리량을 하나의 무차원 지수로 묶은 값입니다. 반경 위치마다 I(r) = 이온 플럭스 × 이온 에너지 가중을 계산하고(에너지 가중은 절대 에너지 스케일이 있을 때 임계값 스퍼터 수율 √(max(E − E_th, 0)), 없을 때 √(상대 에너지) 또는 가중 없음), 면적 가중 평균이 1이 되도록 정규화합니다. 지수는 표면에 도달하는 이온 충격의 반경 분포를 기술합니다.

이 지수는 공정 결과가 아닙니다. 이 시뮬레이터는 식각률·증착률·선택비·프로파일을 예측하지 않습니다 — 그 양들은 라디칼 화학과 표면 반응이 지배하는데 이 솔버는 그것을 풀지 않기 때문입니다. 지수의 에너지 가중이 기술하는 물리적 스퍼터링은 실제 공정에서 여러 항 중 하나일 뿐이며, 지수를 무엇으로 나누어도 공정 결과가 나오지 않습니다.

• Coburn & Winters, J. Appl. Phys. 50, 3189 (1979) — ion/neutral synergy at a surface: the measurement showing why an ion-only index cannot stand in for a surface outcome [verify DOI] [verify]
• Lieberman & Lichtenberg 2005, Ch. 11 — RF sheath physics setting the ion flux and ion energy at the electrode