다중성으로 강화된 옥트리 코드 Treecode2 정적 테스트

초록

**
Treecode2는 매 시간 단계마다 루트 셀을 무작위로 이동·회전·크기조정하여 여러 좌표계에서 평균을 취함으로써 옥트리 기반 중력 계산의 평행·회전 불변성을 통계적으로 회복한다. 본 논문에서는 정적 입자 배치를 이용해 평균화가 포텐셜·가속도 오차를 크게 감소시키는지를 검증하고, 기존 알고리즘 대비 에너지·운동량·각운동량 보존이 향상될 가능성을 제시한다.

**

상세 분석

**
Treecode2가 제안하는 핵심 아이디어는 ‘루트 셀의 무작위 변환(이동, 회전, 스케일)’을 매 타임스텝마다 적용하고, 여러 변환에 대해 얻어진 힘을 평균함으로써 옥트리 구조가 갖는 고정된 좌표계 의존성을 없애는 것이다. 전통적인 옥트리(Barnes‑Hut)에서는 루트 셀을 고정된 위치와 방향에 두기 때문에 전체 시스템을 평행 이동하거나 회전했을 때 계산된 중력장이 달라지는, 즉 뉴턴의 제3법칙 위반과 선형·각운동량 비보존 문제가 발생한다. 이러한 비대칭성은 특히 비구형(예: 얇은 원반) 구조에서 크게 나타나며, Fig. 1에서 보듯이 z축을 따라 작은 변위만으로도 수십 배 차이의 ‘벌크 힘’이 생성된다.

Treecode2는 변환 행렬 Rₖ, 스케일 Sₖ, 이동벡터 Tₖ를 각각 무작위로 선택한다. Rₖ는 정규화된 회전 행렬, Sₖ는 로그 균등분포( |log Sₖ| ≤ log Sₘₐₓ ), Tₖ는 구형 영역 내 균등분포(반경 Tₘₐₓ)에서 추출한다. 입자 좌표 rₚ는 변환 후 r′ₚ = Rₖ Sₖ (rₚ − Tₖ) 로 옮겨진 뒤 트리를 구축한다. 루트 셀은 변환된 좌표계에서 원점에 위치하고, 가장 작은 2의 거듭제곱 크기로 모든 입자를 포함하도록 정의된다. 이렇게 하면 셀 경계와 중심이 부동 소수점 표현에 최적화되고, Sₘₐₓ = √2 정도면 모든 가능한 셀 크기를 포괄한다는 장점이 있다.

각 트리마다 질량, 무게중심, 쿼드러플(moment) 등을 기존 방식과 동일하게 계산한다. 쿼드러플은 트리 셀 내부 입자들의 상대 위치 sₚ = rₚ − r_c 를 이용해 Q_c = ∑ₚ mₚ( sₚ⊗sₚ − sₚ² I ) 로 정의되며, 이는 2차 다중극 전개에 사용된다. 소프트닝 길이 ε를 포함한 개선된 포텐셜식(식 6)은 기존 1차식보다 높은 정확도를 제공한다.

셀을 열어야 할지 결정하는 개방 기준은 전통적인 BH86의 θ-기준( s < ℓ_c/θ )에 셀 중심과 무게중심 사이 거리 d_c 를 추가한 식 9( s < ℓ_c/θ + d_c ) 로 보강한다. 이는 셀 내부 질량이 비대칭적으로 분포될 때 발생하는 ‘detonating satellite’ 현상을 방지하고, θ ≈ 1 수준에서도 안정적인 계산을 보장한다.

무작위 변환을 연속적으로 적용하면 개별 타임스텝에서 발생하는 편향이 서로 상쇄된다. 통계적으로는 트리 구조가 모든 가능한 회전·이동·스케일에 대해 평균화되므로, 기대값으로는 완전한 평행·회전 불변성을 회복한다. 실험에서는 동일한 입자 배치를 10⁴ ~ 10⁵ 번의 무작위 변환으로 평균화했을 때, 포텐셜 오차는 평균 0.5 % 이하, 가속도 오차는 1 % 미만으로 크게 감소함을 확인했다. 특히 얇은 원반에서 나타나는 대규모 벌크 힘이 거의 사라지고, 구형 구름에서는 오차가 거의 무시 수준으로 줄어든다.

이러한 결과는 정적 테스트에 국한되지만, 동적 시뮬레이션에서도 에너지·운동량·각운동량 보존이 크게 개선될 것으로 기대된다. 왜냐하면 평균화된 힘은 시간에 따라 독립적인 오차를 갖고, 오차가 서로 상쇄되면서 누적 오차가 억제되기 때문이다. 또한, 무작위 변환 자체는 O(N log N) 복잡도에 추가적인 연산이 거의 없으며, 회전·스케일·이동을 위한 행렬 연산과 변환만 수행하면 되므로 전체 계산 비용은 거의 변하지 않는다.

결론적으로 Treecode2는 기존 옥트리 코드가 갖는 좌표계 의존성을 ‘다중성(Pluralism)’이라는 통계적 평균화 기법으로 해결하고, 정밀도와 효율성을 동시에 향상시키는 혁신적인 접근법이다. 향후 동적 테스트와 다양한 천체물리 시뮬레이션에 적용해 보존 법칙 개선 효과를 정량화할 필요가 있다.

**