검증 및 문헌 근거 (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과 함께일 때만 유효).

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

반경 이온 플럭스 비균일도

이온 플럭스 반경 프로파일의 평탄도를 나타내며, 웨이퍼 전면에서 공정이 얼마나 균일하게 진행되는지를 나타내는 핵심 지표입니다. 값이 높으면 웨이퍼 중심과 가장자리 간에 식각 깊이, 막 두께, 임계 치수(CD) 편차가 커질 수 있습니다. 시뮬레이터는 반경 프로파일에서 비균일도(%) = (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

식각/증착 속도 상대 프록시

이 값은 절대적인 식각 속도(nm/min)가 아닙니다. 이온 보조 식각(ion-assisted etching)에서 속도는 이온 플럭스와 이온 에너지에 의존한다는 물리적 근거(Coburn & Winters 1979)에 기반한 상대적 비교 프록시입니다. 시뮬레이터는 각 반경 위치에서 이온 플럭스 × √(이온 에너지)를 계산한 후 반경 평균을 취합니다. 표면 화학 반응, 스퍼터 수율, 반응 부산물, 마스크 선택비는 전혀 모델링하지 않으므로, 이 값은 동일 조건 내에서의 상대 경향 비교에만 사용해야 합니다.

• Coburn & Winters, J. Appl. Phys. 50, 3189 (1979) — ion-assisted etching: ion flux and energy dependence of etch rate [verify DOI] [verify]
• Lieberman & Lichtenberg 2005, Ch. 15–16 — etch and deposition rate physics