포트란 기반 QCDMAPT 패키지: 분석적 QCD 계산을 위한 새로운 도구

초록

QCDMAPT_F는 기존에 Maple용으로 제공되던 QCDMAPT 패키지를 포트란(Fortran) 환경으로 옮긴 버전이다. 1‑4 루프 수준까지의 스펙트럼 함수에 대한 명시적 식을 포함하고, 필요한 적분을 수행하는 서브루틴을 제공한다. CERNLIB의 MathLib(D102)만 설치되어 있으면 어떠한 Fortran 77 컴파일러에서도 실행 가능하며, 입력 파라미터 파일을 통해 (M)APT 형태의 시간·공간 전개 함수를 손쉽게 계산한다.

상세 분석

본 논문은 QCD의 비섭동적 영역을 기술하기 위해 널리 사용되는 ‘분산(analytic) 접근법’에 기반한 QCDMAPT 패키지를 포트란 언어로 구현한 QCDMAPT_F 버전을 소개한다. 핵심 물리적 개념은 강한 상호작용의 러닝 커플링 α_s(Q²)를 복소 평면에 전개하여 얻어지는 스펙트럼 함수 ρ(σ)이며, 이는 Källén‑Lehmann 표현을 통해 물리량을 비섭동적으로 재구성한다. 1‑루프 수준에서는 α_s가 로그 형태로 단순히 표현되므로 ρ(σ)도 간단히 적분으로 얻을 수 있다. 그러나 2‑루프 이상에서는 베타 함수 계수가 추가되고, 로그와 다항식이 복합적으로 나타나면서 ρ(σ)의 명시적 형태가 급격히 복잡해진다. 기존 Maple 구현에서는 이러한 복잡한 식을 기호적으로 다루어 자동화했지만, 대규모 수치 계산에서는 Maple의 실행 속도가 제한 요인이었다.

QCDMAPT_F는 이러한 문제점을 해결하기 위해 4‑루프까지의 스펙트럼 함수를 미리 계산한 명시적 식을 제공하고, 이를 Fortran 서브루틴으로 구현한다. Fortran 77 표준을 따르면서도 외부 라이브러리인 CERNLIB MathLib(D102)의 적분 루틴을 활용함으로써, 고정밀 수치 적분을 효율적으로 수행한다. 프로그램 구조는 메인 파일(QCDMAPT_F.f)과 두 개의 입력 파일(QCDMAPT_F.i1, QCDMAPT_F.i2)로 이루어져 있으며, 메인 파일 안에 스펙트럼 함수 정의와 적분 서브루틴이 포함된다. 사용자는 입력 파일에 Λ_QCD, 활성 쿼크 수, 계산하고자 하는 Q² 구간 등을 지정하면, (M)APT의 시간‑유사(expansion) 및 공간‑유사 함수 값을 자동으로 출력한다. 옵션으로는 지정 구간에 대한 연속 데이터 파일을 생성하는 기능도 제공한다.

성능 측면에서 저자들은 샘플 입력(i2)으로 전체 계산을 수행했을 때 약 1분 정도 소요된다고 보고한다. 이는 현대 CPU 환경에서 충분히 실용적인 수준이며, 대규모 파라미터 스캔이나 Monte‑Carlo 기반의 QCD 분석에 바로 적용할 수 있다. 또한 Fortran 기반이므로 고성능 컴퓨팅 클러스터에서 병렬화가 용이하고, 기존 Fortran 기반 코드와의 연동도 자연스럽게 이루어진다.

프로그램의 제한점은 CERNLIB MathLib이 사전 설치되어야 한다는 점이다. 현재 CERNLIB는 유지보수가 중단된 상태이지만, 대부분의 리눅스 배포판에서 여전히 제공되며, 필요 시 소스 코드를 직접 컴파일할 수 있다. 또한, 4‑루프 수준까지의 스펙트럼 함수는 현재까지 알려진 베타 함수 계수와 매스 스키마에 의존하므로, 새로운 고차 루프 결과가 발표될 경우 패키지 업데이트가 필요하다.

전반적으로 QCDMAPT_F는 분석적 QCD 계산을 위한 실용적인 도구로, Maple 기반의 한계를 극복하고 Fortran 환경에서 빠르고 안정적인 수치를 제공한다는 점에서 학계와 실험 분석팀 모두에게 유용할 것으로 기대된다.