격자 몬테카를로 확산 시뮬레이션 정확도 최적화

초록

본 논문은 격자 기반 몬테카를로(Lattice Monte Carlo) 알고리즘이 연속 확산 방정식을 근사하는 과정에서, 시간 간격(Δt)과 격자 간격(Δx) 중 하나를 고정했을 때 다른 파라미터의 최적값이 존재함을 보인다. 1차원에서는 최근접 이웃 이동만으로도 두 번째와 네 번째 모멘트를 정확히 재현하고, 첫 통과 시간 분포의 첫 두 모멘트를 보존한다. 2·3차원에서는 대각선 이동을 허용해야 같은 정확도를 얻으며, 경계와의 충돌 시 단순히 거부하는 대신 경계에 평행하게 투사(projection)하는 방식이 필요하다. 흡수 경계에 대한 처리 방법도 제시한다.

상세 분석

논문은 먼저 Lattice Monte Carlo(LMC) 알고리즘을 연속 확산 방정식의 유한 차분 형태와 비교함으로써, Δt와 Δx가 동시에 0으로 수렴할 때 수렴성을 보장한다는 일반적인 전제를 확인한다. 그러나 실제 시뮬레이션에서는 계산 비용과 메모리 제한으로 두 파라미터 중 하나를 고정하는 경우가 빈번하다. 저자들은 이러한 상황에서 “최적”의 Δt(또는 Δx)값이 존재한다는 점을 수학적으로 증명한다. 1차원에서는 입자가 현재 격자점에서 왼쪽 혹은 오른쪽 인접 격자점으로 이동하는 확률 p=½를 유지하면서, 이동 간격 Δx와 시간 간격 Δt 사이에 다음 관계식
Δt = (Δx)² / (2D)·(1‑α)
(α≈0.5) 를 만족하도록 선택하면, 확산 계수 D를 정확히 재현한다. 이때 확산 과정의 두 번째 모멘트 ⟨x²⟩와 네 번째 모멘트 ⟨x⁴⟩가 연속 해와 일치하고, 고차 모멘트에 대한 오차는 시간에 비례해 감소한다. 또한 첫 통과 시간(first‑passage time, FPT) 분포의 평균과 분산이 정확히 보존되며, 이는 반사 경계와 흡수 경계 모두에서 중요한 검증 지표가 된다.

2차원·3차원으로 차원을 확장할 경우, 단순히 각 축에 대해 독립적인 1차원 이동만 허용하면 네 번째 모멘트까지 정확히 맞추기가 불가능함을 보인다. 따라서 “대각선 이동”(예: (±Δx,±Δx) 혹은 (±Δx,±Δx,±Δx) 형태)을 포함시켜 이동 집합을 확장해야 한다. 대각선 이동을 포함하면 전이 확률을 적절히 조정해 전체 확산 텐서가 등방성을 유지하도록 할 수 있다. 특히 경계 근처에서 대각선 이동이 경계면을 관통하려 할 경우, 이동을 무조건 거부(reject)하는 전통적인 방법은 확산 속도를 인위적으로 감소시켜 오차를 크게 만든다. 대신 이동 벡터를 경계에 평행하게 투사(projection)함으로써 입자는 경계에 머물면서도 올바른 확산 통계량을 유지한다. 이 투사 방식은 반사 경계와 흡수 경계 모두에 적용 가능하며, 흡수 경계에서는 투사된 위치가 경계면에 도달하면 즉시 소멸시켜 첫 통과 시간 통계가 정확히 재현된다.

수치 실험에서는 다양한 Δx와 Δt 조합에 대해 평균 제곱 변위, 네 번째 모멘트, 그리고 FPT 평균·분산을 측정하였다. 최적 파라미터 집합에서는 오차가 최소화되고, 특히 장시간 시뮬레이션에서 고차 모멘트 오차가 선형적으로 감소함을 확인했다. 또한 대각선 이동을 포함한 경우와 단순 거부 방식을 비교했을 때, 전자는 경계 근처에서의 확산 속도가 실제 연속 모델과 거의 일치함을 보여준다.

결론적으로, LMC 알고리즘의 정확도를 극대화하려면 (1) 시간·격자 간격 사이에 이론적으로 도출된 최적 비율을 유지하고, (2) 다차원에서는 대각선 이동을 포함시키며, (3) 경계와의 충돌 시 투사 방식을 적용해야 한다는 실용적인 가이드라인을 제공한다. 이러한 설계 원칙은 복잡한 형상이나 변동 경계가 존재하는 물리·생물 시스템의 확산 시뮬레이션에 바로 적용 가능하다.