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

초록

SC, BCC, FCC 격자점을 가우시안 잡음으로 변위시켜 3차원 보로노이 셀의 면수, 부피·면적 통계와 구형도(등축비)를 조사하였다. BCC는 작은 잡음에도 위상 안정성을 유지하지만 SC와 FCC는 위상이 급변한다. 잡음 강도가 중간(α≈0.5) 이상이면 세 구조의 통계가 거의 동일해지고, 강한 잡음(α>2)에서는 포아송 보로노이 한계에 수렴한다. 면적·부피 관계는 지수 1.67 정도의 비정상적 스케일링을 보이며, 면수별로 분석하면 이 효과가 크게 감소한다. 모든 측정량은 2‑파라미터 감마 분포로 잘 설명된다.

상세 분석

본 연구는 3차원 격자 구조인 단순 입방체(SC), 체심 입방체(BCC), 면심 입방체(FCC)를 기준점으로 삼아, 각 격자점에 동일한 표준편차 σ를 갖는 가우시안 잡음을 독립적으로 부여한다. 잡음 강도는 무차원 파라미터 α=σ/a₀ (a₀는 격자 상수)로 정의되며, α=0은 완전한 결정, α→∞는 완전 무작위 포아송 점군에 해당한다. Monte‑Carlo 시뮬레이션을 10⁶ 셀 규모로 수행해 각 α에 대해 면수(f), 면적(A), 부피(V), 주변 길이(P) 및 등축비(Q)=36πV²/A³ 등을 통계적으로 측정하였다.

위상적 관점에서 BCC는 트렁케이티드 옥타헤드론(14면)이라는 비퇴화 다면체를 기본 셀로 갖는다. 이는 작은 변위에도 셀의 면수와 연결 구조가 변하지 않아 위상적 안정성을 보인다. 반면 SC(6면)와 FCC(12면)는 각각 정육면체와 러미노드케이드(12면)라는 퇴화된 형태이므로, 미세한 잡음만으로도 면수 변동이 급격히 일어나며, 평균 면수는 α≈0.1부터 급격히 증가한다.

계량적 결과는 다음과 같다. α≲0.5 구간에서 BCC와 FCC의 평균 면적 ⟨A⟩은 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‑파라미터 모델과 일치하며, 복잡한 다변량 상관관계를 단순화하는 강력한 도구임을 확인한다.

특이한 스케일링 현상으로, A와 V 사이의 관계가 ⟨A⟩∝⟨V⟩^{β} 형태를 보이며, β≈1.67(≈5/3)으로 전통적인 기하학적 기대치(β=2/3)와 크게 다르다. 이는 셀 형태의 변동성, 특히 면수와 등축비가 크게 변동하면서 발생한다. 면수별로 데이터를 분리해 회귀하면 β는 1.5에 가까워지며, 면수가 많을수록 셀은 더 구형에 가까워져 스케일링이 정상화된다.

결론적으로, 격자 구조에 대한 점진적 무작위화는 위상·계량적 특성을 연속적으로 변화시키며, BCC는 위상적 강인성을, FCC는 면적 최소화를, SC는 잡음에 대한 민감성을 보여준다. 또한, 2‑파라미터 감마분포와 비정상적 A–V 스케일링은 복합 구조 물질의 미세구조 설계와 데이터 압축 등에 활용될 수 있는 중요한 통계적 법칙임을 제시한다.