고차원 접촉수 상한을 위한 고정밀 반정밀계획법

논문은 먼저 접촉수(kissing number) 문제의 배경을 소개한다. n차원 유클리드 공간에서 단위 구가 중심 구와 접촉할 수 있는 최대 개수는 차원 1, 2, 3, 4, 8, 24에서만 정확히 알려져 있다. 기존에는 Delsarte 방법과 그 변형을 이용해 상한을 구했지만, 차원 5 이상에서는 큰 격차가 남아 있었다. Bachoc과 Vallentin은 Jacobi 다항식과 다항식 합의 제곱(SOS) 제약을 결합한 반정밀계획(SDP) 모델을 제시하였다. 이 모델은 매개변수 d에 따라 일련의 최적화 문제 s_d(n)을 정의하고, s_d(n)의 최소값이 접촉수의 상한이 된다. d를 크게 잡을수록 상한은 단조 감소하고 실제 접촉수에 근접한다. 저자들은 이 모델을 차원 3부터 24까지 적용했으며, 특히 d = 11, 12, 13, 14까지 계산했다. 계산 과정에서 가장 큰 난관은 수치적 불안정성이었다. 일반적인 SDP 솔버는 30자리 이하의 정밀도에서 수렴이 멈추었고, 특히 d > 10에서는 신뢰할 수 없는 결과가 나왔다. 이를 해결하기 위해 GNU Multiple Precision Arithmetic Library(GMP)를 이용한 SDPA‑GMP를 사용했으며, 200~300비트(≈60~90자리) 정밀도로 연산했다. 파라미터 λ_star를 10⁰~10⁴ 사이에서 조정하고, 종료 기준을 10⁻³⁰으로 설정해 5~10주 정도의 연산 시간을 투자했다. 표 1은 기존에 알려진 상한과 새로 계산된 상한을 비교한다. 차원 3, 4, 8, 24는 이미 정확한 값이 알려져 있어 12.000… 등으로 표시되었고, 차원 5부터 15까지는 기존 상한보다 약 0.1%~0.5% 정도 개선된 값을 제공한다. 예를 들어 차원 5의 경우 기존 상한이 45.06107293…이었으나, d = 14에서 44.99899685…로 낮춰졌다. 차원 6, 7, 9, 10 등에서도 유사한 개선이 관찰된다. 특히 차원 16에 대한 결과가 눈에 띈다. 기존에 알려진 상한은 7515.95264…(d = 11)에서 7355.809036…(d = 14)까지 감소하였다. 이 값은 Conway‑Sloane이 제시한 평균 세타 급수 1+7680q³+4320q⁴+27648q⁵+61440q⁶+…를 만족하는 주기적 점집합이 존재한다면 최소 7680이어야 한다는 조건과 모순된다. 따라서 해당 급수를 갖는 점집합은 존재하지 않으며, 이는 Conway‑Sloane의 추측을 완전히 증명한 것이다. 논문은 또한 차원 4에 대한 흥미로운 질문을 제기한다. s_d(4) 값이 d→∞일 때 24에 수렴한다면, D₄ 루트계가 24점 배치의 유일 최적 배치가 된다. 현재 s₁₅(4)=23.06274…이라는 결과는 수렴이 느리지만 가능성을 시사한다. 만약 이 가설이 증명된다면, D₄ 루트계는 보편 최적성(Universally Optimal) 논쟁에서도 중요한 역할을 할 것이다. 마지막으로 저자들은 계산에 사용된 하드웨어와 소프트웨어 환경을 상세히 기술한다. Intel Xeon 2코어와 Quad 프로세서를 이용해 각 d값당 5~10주 정도 소요됐으며, SDPA‑GMP의 128비트 정밀도 버전은 속도는 느리지만 정확도는 충분히 확보했다. 이후 128비트 버전이 개선된 128‑bit(quad) 버전이 출시되어 d = 13, 14 계산 시간을 1~2.5주로 단축했다. 결론적으로, 고정밀 SDP 계산은 접촉수 상한을 기존보다 현저히 개선할 수 있음을 보여준다. 차원 16에서의 비존재 결과는 격자 이론, 모듈러 형식, 그리고 구체적 포장 문제에 새로운 통찰을 제공한다. 향후 더 높은 차원이나 더 큰 d값에 대한 계산이 가능해지면, 접촉수 문제뿐 아니라 관련된 최적화 문제 전반에 큰 영향을 미칠 것으로 기대된다.