WFC3/UVIS Geometric Distortion -- Time Evolution of Linear Terms w.r.t Gaia

기술 요약: ISR 2024-15 – WFC3/UVIS 기하학적 왜곡: Gaia 대비 선형 항의 시간적 진화

문제 정의
정밀한 하늘 좌표, 시차(parallax), 고유 운동을 도출하기 위해서는 허블 우주 망원경(HST)의 WFC3에 대한 정확한 기하학적 왜곡 모델링이 매우 중요하다. 현재의 왜곡 솔루션(IDCTAB, NPOLFILE, D2IMFILE 참조 파일에 저장됨)은 4차 다항식과 비다항식 룩업 테이블을 활용하고 있으나, 이의 선형 항(이동, 회전, 스케일, 스큐)이 전체 운용 기간 동안 Gaia DR3 카탈로그를 기준으로 얼마나 안정적인지에 대한 종합적인 평가는 이루어지지 않았다. 선형 항에서 유의미한 시간적 진화가 발생할 경우, 드리즐링(drizzling)을 통해 다중 에포크(multi-epoch) 데이터를 결합할 때 소스의 흐림 현상이나 부정확한 천체 측정을 초래할 수 있다. 이전 연구들은 표본 크기가 작거나, 초기 Gaia 릴리스에서 충분한 신뢰도를 확보하지 못한 단일 대상($\omega$ Centauri)에 의존했다는 한계가 있었다.

방법론
저자들은 2009년부터 2022년 사이에 수집된 5개의 광대역 필터(F275W, F336W, F438W, F606W, F814W)에 걸친 7,491개의 WFC3/UVIS 풀 프레임 노출 데이터셋을 분석하였다. 방법론은 다음과 같다:

1. 소스 검출: 개별 FLC 파일을 hst1pass 루틴으로 실행하여 픽셀 위치와 밝기를 포함한 고정밀 소스 카탈로그를 생성한다.
2. 참조 매칭: 이미지 영역 내의 소스에 대해 Gaia DR3 카탈로그를 쿼리하고, 관측 에포크에 맞춰 고유 운동을 보정한 Gaia 위치를 적용한다.
3. 정렬: tweakwcs(및 tweakreg) 루틴을 사용하여 WFC3/UVIS 카탈로그를 Gaia 참조값에 정렬한다. 이 과정은 XYXYMatch 알고리즘을 사용하여 소스를 교차 매칭하고, 접평면(tangent plane)에서의 잔차를 최소화하는 아핀 변환(이동, 회전, 스케일, 스큐)을 계산한다.
4. 분석: 표준 IDCTAB을 사용하여 생성된 초기 WCS와 Gaia에 정렬된 WCS로부터 도출된 선형 변환 파라미터를 비교한다. 저자들은 참조 파일의 오차를 결정하기 위해 차이(또는 스케일의 경우 비율)를 계산하였다.
5. 필터링: RMS 오차(RMSE)가 0.2 픽셀 미만이고 매칭된 소스가 최소 50개 이상인 이미지만 포함하도록 필터링하였다. 선형 추세를 피팅하기 전 이상치를 제거하기 위해 시그마 클리핑(sigma clipping)을 적용하였다.

주요 기여

• 대규모 시계열 분석: 본 연구는 단일 대상을 사용하는 대신 약 7,500개의 다양한 이미지 샘플을 사용하여 13년의 베이스라인 동안 WFC3/UVIS 선형 왜곡 항에 대한 최초의 종합적인 평가를 제공한다.
• 필터 의존적 오프셋의 정량화: 현재의 단일 참조 파일 왜곡 모델에서 고려되지 않은 스케일 항의 체계적이고 필터 의존적인 오프셋을 식별하고 정량화하였다.
• 운용 권고 사항: 고정밀 천체 측정을 위해 이미지 재정렬의 필요성과 구체적인 피팅 기하 구조(fitting geometries)에 관한 구체적이고 실행 가능한 지침을 제공한다.

결과
분석 결과, 선형 항은 시간적 진화를 보이지만, 이러한 변화는 매끄럽고 결정론적인 드리프트보다는 본질적인 산포(scatter)와 포인팅 불확실성에 의해 주로 지배되는 것으로 나타났다.

• 이동(Shift): 이동 항은 상당한 산포(RMS $\approx$ 5 픽셀)를 보였으며, 이는 주로 망원경의 포인팅 및 가이드 스타 카탈로그(GSC) 좌표의 불확실성에 기인한다. 시간적 기울기는 작다 (X축 0.17 픽셀/년, Y축 -0.27 픽셀/년).
• 회전(Rotation): 2009년부터 2017년 중반까지 $0.003 \pm 0.004$ 도의 작고 안정적인 오프셋이 관찰되었다. 2017년 중반부터 2022년까지는 오프셋이 최대 0.01도까지 증가하며 큰 산포를 보였는데, 이는 HST 포인팅 지터(jitter)가 증가한 시기와 상관관계가 있다. 선형 기울기는 무시할 만한 수준이다 (-0.00011 도/년).
• 스케일(Scale): 13년 동안 약 0.15–0.2 픽셀의 작은 시간적 진화가 감지되었다. 그러나 더 중요한 발견은 스케일 항에서 나타나는 필터 의존적 오프셋으로, 모든 에포크에서 F275W와 F814W 필터 간에 최대 0.3 픽셀의 차이가 발생했다. 이는 현재의 IDCTAB이 필터별 기하학적 왜곡을 완전히 반영하지 못하고 있음을 시사한다.
• 스큐(Skew): 스큐 항은 13년 동안 약 0.2 $\pm$ 0.3 픽셀만큼 진화한다. 시간적 변화는 일반적으로 RMS 산포와 비슷하거나 그보다 작으며, 명확한 필터 의존적 오프셋은 발견되지 않았다.

의의 및 권고 사항
본 논문은 선형 왜곡 항의 시간 의존적 진화가 포인팅 불확실성 및 측정치의 본질적 산포에 비해 일반적으로 작다고 결론짓는다. 그러나 필터 의존적 스케일 오프셋과 회전/이동 항의 큰 산포는 사용자에 의한 능동적인 보정을 필요로 한다.

저자들은 관측자들에게 다음과 같이 권고한다:

1. 이미지 재정렬: FLT/FLC 이미지를 재정렬하기 위해 DrizzlePac의 tweakreg 루틴을 사용할 것.
2. 적절한 기하 구조 피팅:
   • 표준 정렬의 경우, 이동, 회전, 스케일을 포함하는 4-파라미터 피팅(fitgeometry='rscale')을 사용할 것.
   • 장기 베이스라인에 대한 고정밀 천체 측정을 수행할 경우, 스큐와 독립적인 x/y 스케일을 모두 고려하는 6-파라미터 피팅(fitgeometry='general')을 사용할 것.
3. MAST 프로세싱 확인: MAST 파이프라인 제품의 NMATCHES 키워드를 확인할 것. Gaia와의 매칭 수가 20개 미만인 경우, WCS를 "a priori" 솔루션으로 리셋한 후 수동 재정렬을 수행하는 것을 고려해야 한다.

이 보고서는 MAST 파이프라인이 정렬 단계를 수행하지만, 특히 서로 다른 필터나 에포크 간의 데이터를 결합할 때 최고 수준의 정밀도를 달 equilibrium 하기 위해서는 권장된 피팅 기하 구조를 사용하는 수동 개입이 종종 필요함을 강조한다. 향후 연구에서는 시간 의dependent 왜곡 솔루션 및 IDCTAB의 글로벌 스케일 보정 구현을 조사할 예정이다.