입체 랜덤 보로노이 테셀레이션: 격자에서 포아송까지의 연속 변환

본 논문은 3차원 유클리드 공간에서 세 가지 전형적인 입방 격자(SC, BCC, FCC)를 출발점으로 삼아, 각 격자점에 동일한 통계적 특성을 갖는 가우시안 잡음을 가함으로써 점군을 점진적으로 무작위화한다. 잡음 강도는 무차원 파라미터 α로 정의되며, α=0은 완전한 결정, α→∞는 포아송 점군에 해당한다. 연구는 크게 네 부분으로 구성된다. 1. **이론적 배경 및 방법론** - 보로노이 테셀레이션의 기본 정의와 3차원 다면체의 위상 관계(V+E+F=2) 등을 정리하고, 면수(f), 면적(A), 부피(V), 주변 길이(P), 등축비(Q) 등 주요 통계량을 정의한다. - 가우시안 잡음은 격자 상수 a₀에 대한 비율 σ/a₀=α로 설정하고, Monte‑Carlo 시뮬레이션을 통해 10⁶ 셀 규모의 샘플을 생성한다. 각 α에 대해 10⁴개의 독립 실현을 평균화해 통계적 오차를 최소화한다. 2. **위상적 안정성 분석** - BCC는 트렁케이티드 옥타헤드론(14면)이라는 비퇴화된 셀을 갖기 때문에, α가 아주 작아도 면수와 연결 구조가 변하지 않는다. 이는 “위상적 강인성”이라 부른다. - 반면 SC(정육면체, 6면)와 FCC(러미노드케이드, 12면)는 퇴화된 형태이므로, 아주 미세한 잡음(α≈0.05)만으로도 면수가 급격히 변한다. 특히 SC는 면수 6에서 8~10 사이로 빠르게 전이한다. 3. **계량적 특성 및 확률분포** - 평균 면적 ⟨A⟩와 평균 부피 ⟨V⟩는 α에 따라 서로 다른 스케일링을 보인다. BCC와 FCC는 α≲0.5 구간에서 ⟨A⟩∝α²(즉, 2차 스케일)이며, SC는 ⟨A⟩가 최소가 되는 최적 α≈0.3을 가진다. - α>0.5가 되면 세 구조 모두 평균 면수 ⟨f⟩≈15.5, 평균 면적 ⟨A⟩≈5.3, 평균 부피 ⟨V⟩≈1.0(격자 부피 단위) 수준으로 수렴한다. α>2에서는 포아송 보로노이 한계값(⟨f⟩≈15.54, ⟨Q⟩≈0.66 등)과 거의 일치한다. - 면수, 면적, 부피, 등축비 각각의 확률밀도함수는 두 파라미터 감마분포 Γ(k,θ)로 매우 높은 적합도를 보이며, 이는 기존 2‑파라미터 모델과 일치한다. 4. **비정상적 A–V 스케일링과 면수 의존성** - 전체 셀을 대상으로 하면 A와 V 사이의 관계가 ⟨A⟩∝⟨V⟩^{β} 형태를 띠며, β≈1.67(≈5/3)으로 전통적인 기하학적 기대치(β=2/3)와 크게 다르다. 이는 셀 형태가 크게 변동하고, 특히 면수가 많을수록 셀은 구형에 가까워지는 현상 때문이다. - 면수별로 데이터를 분리해 회귀하면 β는 1.5에 가까워지고, 면수가 많을수록 등축비 Q가 증가한다. 따라서 “면수-구형도-스케일링” 삼각관계가 존재한다는 결론에 도달한다. **주요 결론 및 의의** - BCC는 위상적으로 가장 강인하며, 작은 잡음에도 셀 구조가 유지된다. 이는 물질의 미세구조 설계에서 강인성을 요구하는 경우 유용하다. - FCC는 잡음에 대해 면적 최소화 특성을 보이며, 최적의 잡음 강도에서 평균 면적이 최소가 된다. 이는 재료의 표면 에너지 최소화와 연관될 수 있다. - SC는 잡음에 매우 민감해 위상 변이가 빠르게 일어나므로, 결정 구조가 쉽게 붕괴되는 경우를 모델링하는 데 적합하다. - 모든 경우에 2‑파라미터 감마분포가 통계량을 잘 설명하므로, 복잡한 3차원 구조의 확률적 모델링에 강력한 도구가 된다. - 비정상적 A–V 스케일링은 셀 형태 변동이 큰 시스템(예: 다공성 매체, 금속 폼)에서 물리적 특성(전도성, 강도)과 연관된 새로운 스케일링 법칙을 제시한다. 본 연구는 격자 → 무작위 전이 과정을 정량적으로 규명함으로써, 재료 과학, 물리학, 데이터 압축 등 다양한 분야에서 보로노이 테셀레이션을 활용하는 데 이론적·실용적 기반을 제공한다.