VFISV: 초고속 태양광자기장 역전 코드

초록

VFISV는 SDO/HMI에서 매 10분마다 전면 4096×4096 픽셀의 전자기장 벡터를 제공하도록 설계된 초고속 스톡스 벡터 역전 코드이다. 밀레니엄-에디슨 모델을 기반으로 파라미터 수를 최소화하고, 레벤버그‑마콰르트 최적화와 특수 가중치, 정규화 기법을 적용해 16 백만 픽셀을 10분 안에 처리한다. 코드 구현 세부와 히노데 SP 데이터 검증 결과도 제시한다.

상세 분석

VFISV는 태양광학에서 가장 널리 사용되는 밀레니엄‑에디슨(ME) 대기 모델을 전제로 하여, Stokes I, Q, U, V 네 개의 스펙트럼 라인을 동시에 역전한다. 핵심 아이디어는 물리적 파라미터(자기장 세기 B, 경사 θ, 방위 φ, 속도 v, 온도·압력 비율 η₀, 도플러 폭 Δλ_D, 소스 함수 기울기 β 등)를 최소한의 자유도로 압축하고, 레벤버그‑마콰르트(LM) 알고리즘을 변형해 수렴 속도를 극대화하는 것이다.

첫 번째 최적화는 초기 추정값을 관측된 Stokes V의 대칭성 및 Q/U의 비대칭성으로부터 빠르게 도출하는 ‘quick‑guess’ 절차이다. 이는 전통적인 전역 탐색보다 훨씬 적은 반복 횟수로 파라미터 공간에 진입하게 해준다. 두 번째는 파라미터 스케일링과 정규화이다. 각 파라미터에 대해 물리적 범위와 민감도를 고려해 스케일 팩터를 부여함으로써 LM 행렬의 조건수를 크게 낮추고, 수치적 불안정을 방지한다.

또한 VFISV는 가중치 행렬 W를 스펙트럼 라인 중심과 주변부에 차등 적용한다. 중심부는 신호‑대‑노이즈 비가 높아 높은 가중치를 부여하고, 외곽은 낮은 가중치로 잡아 불필요한 잡음에 의한 오버피팅을 억제한다. 이와 더불어 ‘regularization term’으로 알려진 제약 함수를 도입해, 비물리적 해를 억제하고 특히 B와 θ 사이의 상관관계를 완화한다.

코드 구조는 Fortran 90 기반이며, OpenMP를 이용한 다중 스레드 병렬화를 적용했다. 각 픽셀에 대한 역전은 독립적이므로, 16 CPU 코어에서 16 백만 픽셀을 10 분 이내에 처리할 수 있다. 메모리 사용량도 최소화하기 위해 한 번에 한 라인(또는 블록)씩 읽어 들이고, 중간 결과를 즉시 버퍼에 기록한다.

성능 검증을 위해 히노데(SOT/SP) 데이터와 비교했으며, VFISV는 평균 절대 차이가 50 G 이하, 경사 차이가 5° 이하인 높은 정확도를 보였다. 특히 강한 자기장 영역에서의 속도와 소스 함수 파라미터는 기존 코드보다 더 일관된 값을 산출했다. 처리 속도는 동일 데이터셋에 대해 기존 MERLIN 코드 대비 5배 이상 빠르며, 이는 실시간 전역 관측 파이프라인에 필수적인 요소다.

종합하면 VFISV는 물리적 모델 단순화, 수치 최적화, 병렬 처리라는 세 축을 결합해, 대용량 전면 관측을 실시간으로 지원하는 혁신적인 역전 솔루션이다. 향후 고해상도 관측기와 결합하면, 태양 자기장 동역학 연구에 새로운 가능성을 열어줄 것으로 기대된다.