기술적 요약: 리만 제타 함수에 대한 단순하고 정확한 근사

문제 제기

허수부 $t = \text{Im}(s)$가 큰 값에 대한 리만 제타 함수 $\zeta(s)$의 계산은 상당한 어려움을 안고 있습니다. 오일러 - 마클로린 합계 방법은 엄밀한 오차 한계를 제공하지만, 그 계산 복잡도가 $O(|s|)$로 증가하여 큰 $t$ 값에서는 실용적이지 않습니다. 큰 $t$ 값에 대한 표준적인 접근법은 리만 - 지겔 (Riemann-Siegel) 공식으로, 주 합과 나머지 항에 대한 점근 전개식을 활용하여 복잡도를 $O(t^{1/2})$로 줄입니다. 그러나 리만 - 지겔 공식의 고차 보정항은 특수 함수 (구체적으로 함수 $F(z)$) 의 복잡한 미분항을 포함하며, 이는 특이점과 복잡한 공식으로 인해 수치적으로 평가하기 어렵습니다. 또한, 갈웨이 (Galway) 의 것과 같은 기존 고정밀 알고리즘은 최적의 정확도를 달성하기 위해 여러 매개변수를 조정해야 하는 경우가 많습니다.

본 논문은 고차 리만 - 지겔 보정에 따른 구현의 어려움 없이 임의의 수직 띠 (임계 띠 포함) 에 걸쳐 계산 효율성 (적은 항 수 요구) 과 높은 정확도를 모두 갖춘 근사법의 필요성을 다룹니다.

방법론

저자는 리만 - 지겔 공식의 주 합과 나머지 항에 대한 새로운 근사를 결합한 $\zeta_p(s)$로 표시되는 근사 함수 군을 개발했습니다. 핵심 방법론은 다음 단계들에 기반합니다:

1. 적분 표현: 정확한 리만 - 지겔 공식의 나머지 항은 $\theta = e^{-i\pi/4}$인 복소수 측정치 $\eta(dx) = \frac{e^{-\pi x^2}}{\cosh(\pi \theta x)} dx$를 포함하는 적분 $I_M(s)$로 표현됩니다.
2. 가우스 구적법: 점근 급수 전개를 사용하는 대신, 저자는 이 적분을 가우스 구적법 규칙을 사용하여 근사합니다. 목표는 $4p+2$개의 특정 시험 함수에 대해 구적법이 정확하도록 노드 집합 $\{x_{p,j}\}$과 가중치 집합 $\{\omega_{p,j}\}$을 찾는 것입니다.
3. 계수 유도: 계수 $\{\omega_{p,j}\}$와 $\{\lambda_{p,j}\}$ (여기서 $\lambda_{p,j} = \theta x_{p,j}$) 는 선형 방정식 시스템을 풀어서 유도됩니다. 이 시스템은 $[-1, 1]$ 내의 등간격 점 $y_{p,k}$에 대해 $f_k(x) = e^{-\pi x^2 + 2\pi \theta y_{p,k} x}$ 형태의 함수에 대해 구적법이 정확하도록 요구함으로써 구성됩니다.
   • 문제는 모르델 (Mordell) 적분 $H(y)$의 모멘트에 의해 정의된 선형 범함수 $L$에 대한 가우스 구적법의 노드와 가중치를 찾는 문제로 변환됩니다.
   • 이러한 노드와 가중치는 3 항 재귀 관계를 통해 생성된 직교 다항식을 사용하여 계산됩니다.
4. 근사 공식: 최종 근사는 다음과 같이 정의됩니다:
   $$\zeta_p(s) = \sum_{n=1}^{N_t} n^{-s} + \chi(s) \sum_{n=1}^{N_t} n^{s-1} - \frac{(-1)^{N_t}}{2} \left[ I_{N_t+1/2, p}(s) + \chi(s) I_{N_t+1/2, p}(1-s) \right]$$
   여기서 $I_{M,p}(s)$는 지수함수, 로그함수와 같은 기본 함수와 미리 계산된 복소수 계수만을 포함하는 적분항의 구적법 근사입니다.

주요 기여

• 간소화된 구현: 제안된 근사는 단일 정수 매개변수 $p$와 미리 계산된 계수 집합 $\{\omega_{p,j}\}$ 및 $\{\lambda_{p,j}\}$에 의존합니다. 고차항을 위해 특수 함수의 복잡한 미분항을 평가해야 하는 리만 - 지겔 공식과 달리, $\zeta_p(s)$는 기본 함수와 감마 함수만 요구합니다.
• 미리 계산된 계수: 저자는 $1 \le p \le 30$에 대한 계수와 $p=150$까지의 선택된 값에 대한 계수를 제공합니다. 이러한 계수는 중간 크기를 가지며 급격히 감소하는 경향을 보여 안정적인 수치 평가를 용이하게 합니다.
• 미분 근사: 논문은 $N_t$의 바닥 함수를 미분하는 것과 관련된 불연속성을 피하기 위해 고정된 $N_t$에서 근사 공식을 $s$에 대해 미분함으로써 $\zeta'(s)$를 근사하는 방법론을 확장합니다.
• 하이브리드 알고리즘: $s$의 영역에 따라 직접 합계, 오일러 - 마클로린 합계, 그리고 $\zeta_p(s)$ 근사 사이를 전환하여 전체 복소평면에서 $\zeta(s)$를 계산하는 전략이 제안됩니다.

결과

논문은 $\zeta_p(s)$의 정확도에 관한 광범위한 수치적 증거를 제시합니다:

• 임계 띠 내 정확도: 임계 띠 ($0 \le \sigma \le 1$) 에서의 수치 실험은 $p$가 증가함에 따라 오차 $|\zeta_p(s) - \zeta(s)|$가 현저히 감소함을 보여줍니다.
  • $p=10$인 경우, $t > 250$에서 오차는 $10^{-15}$ 미만입니다.
  • $p=50$인 경우, $t > 4000$에서 오차는 $10^{-100}$ 미만입니다.
  • $p=150$인 경우, $t > 6900$에서 오차는 $10^{-300}$ 미만입니다.
• 오차 패턴: 오차는 구간 $[t_n, t_{n+1}]$ (여기서 $t_n = 2\pi n^2$) 의 끝점에서 그리고 이러한 구간 내의 $4p$개의 등간격 점에서 소멸하는 뚜렷한 패턴을 보입니다. 이 행동은 지배적인 오차 항이 구적법 노드에서 0 이 되는 $H_p(B(t)) - H(B(t))$에 비례한다는 사실에 의해 분석적으로 설명됩니다.
• 리만 - 지겔과의 비교: 중간 크기 $t$ (예: $10^2 \le t \le 10^4$) 의 경우, 근사 $\zeta_7(s)$는 $K=5$인 리만 - 지겔 근사 ($\zeta^{RS}_5$) 보다 더 정확한 것으로 나타납니다. $\zeta^{RS}_5$는 $\zeta_7$ ($O(t^{-1/4})$) 에 비해 더 빠른 점근 감소율 ($O(t^{-13/4})$) 을 가지지만, $\zeta_7$은 고정밀도를 위해 훨씬 적은 미리 계산된 계수 (15 개 대 약 100~170 개) 만 필요로 하며 구현이 더 쉽습니다.
• 정밀도 한계: 배정밀도 (double-precision) 구현에서는 $t > 200$에서 반올림 오차가 근사 오차를 지배합니다. 그러나 사중정밀도 (quadruple precision) 에서는 근사가 훨씬 더 큰 $t$ (예: $t > 10^4$에서 오차 $< 10^{-29}$) 까지 높은 정확도를 유지합니다.
• 성능: 하이브리드 알고리즘의 MATLAB 구현은 내장된 MATLAB zeta 함수보다 수백 배에서 수천 배 더 빠른 것으로 보고됩니다.

중요성과 주장

본 논문은 이러한 근사법이 큰 $t$에 대해 임계 띠 및 기타 수직 띠에서 $\zeta(s)$와 $\zeta'(s)$의 고정밀도 계산을 가능하게 한다고 주장합니다. 그 중요성은 구현 복잡도와 정확도 사이의 절충점에 있습니다. 이 방법은 고차 보정항의 계산 부담을 피하면서도 광범위한 $t$에 대해 경쟁력 있거나 우수한 정확도를 유지하는 리만 - 지겔 공식에 대한 "단순한" 대안을 제공합니다.

저자는 이 작업에서 $\zeta_p(s)$에 대한 엄밀한 오차 한계를 제공하지는 않았지만, 수치적 증거가 광범위하다고 지적합니다. 이 방법은 메르텐스 추측의 반증과 같은 수학적 발견에 필수적인 고정밀도 계산을 위한 실용적인 도구로 제시됩니다. 논문은 연구 커뮤니티의 즉시 사용을 용이하게 하기 위해 MATLAB, Python, Fortran 으로 구현된 오픈 소스 코드를 제공하며 마무리합니다.