From the Virgo interferometer calibration to the bias and uncertainty of the h(t) detector strain during the O4 run

기술 요약: Virgo 간섭계 교정 및 O4에서의 $h(t)$ 변형률 재구성

문제 정의
2015년 첫 중력파 검출 이후 지상 기반 간섭계의 민감도는 크게 향상되었으며, 이에 따라 검출기 변형률 신호인 $h(t)$를 정확하게 재구성하기 위해 매우 정밀한 교정이 필요해졌다. Virgo 간섭계는 현수된 테스트 질량 사이의 극미한 길이 변화(약 $10^{-18}$ m 단위)를 측정한다. 이러한 신호를 검출하기 위해 간섭계는 피드백 루프와 전자기(EM) 액추에이터를 사용하여 특정 관측 구성을 유지해야 한다. 다크 프린지(dark fringe) 신호로부터 $h(t)$를 재구성하려면 모든 제어 신호의 기여도와 거울 액추에이터의 주파수 응답 모델을 고려해야 한다. 이러한 모델은 거울 변위에 대한 독립적인 참조를 필요로 한다. 주요 과제는 관측 대역(10 Hz ~ 2 kHz) 전체에 걸쳐 정밀하고 주파수 종속적인 교정을 달성하는 것과, 저지연(low-latency) 및 오프라인 데이터 분석을 지원하기 위해 재구성된 변형률의 편향(bias)과 불확실성을 정확하게 정량화하는 것이다.

방법론
본 논문은 O4 런(특히 2024년 4월에 시작된 O4b)을 위해 구현된 교정 절차 및 재구성 모니터링 전략을 설명한다.

1. 액추에이터 교정 참조:

   • 뉴턴 교정기 (Newtonian Calibrator, NCal): 회전하는 질량을 사용하여 알려진 중력장을 유도하며, 0.12%의 불확실도로 거울 변위에 대한 절대적 참조 역할을 한다. 그러나 150 Hz까지의 주파수로 제한되며 단일 주파수 동작만 가능하다.
   • 광자 교정기 (Photon Calibrator, PCal): 보조 레이저의 복사압을 사용하여 변위를 유도한다. 이는 10 Hz에서 2 kHz 대역에서 NCal을 보완하며 0.48%의 정밀도를 가진다. 변위 $\Delta x_{PCal}$은 기계적 응답 모델, 레이저 입사각, 반사 전력 및 센서 체인 응답을 사용하여 계산된다.
   • 교정 과정: PCal은 EM 액추에이터를 교정하는 데 사용된다. 이는 참조 액추에이터(PCal)의 주입에 대한 간섭계의 응답과 테스트 대상 액추에이터의 응답을 비교함으로써 수행되며, 이를 통해 새로운 액추에이터의 응답 모델을 추출할 수 있다.

2. 변형률 재구성 ($h(t)$):

   • 재구성 알고리즘인 Hrec는 교정된 액추에이터 응답을 사용하여 다크 프린지 신호로부터 종방향 제어 신호의 기여분을 차감한다.
   • 모니터링: 재구성 품질은 전달 함수 $h_{raw}/h_{inj}$를 통해 모니터링된다. 여기서 $h_{raw}$는 재구성된 변형률이고, $h_{inj}$는 액추에이터의 교정 모델만을 사용하여 계산된 변형률이다. 크기(modulus)가 1에서 벗어나거나 위상(phase)이 0에서 벗어나는 현상은 재구성 편향을 나타낸다.
   • 주입 신호: 모니터링을 위해 두 가지 유형의 신호가 사용된다.
     • 영구 라인 (Permanent Lines): 지속적으로 주입되는 16개의 사인파 신호 (16 Hz ~ 1 kHz).
     • 주간 라인 (Weekly Lines): 매주 1분 동안 주입되는 27개의 사인파 신호 (18 Hz ~ 1158 Hz).

3. 편향 및 불확실성 계산:

   • **온라인 언바이싱(unbiasing)**을 가능하게 하기 위해 주파수 종속적 편향 및 불확실성을 계산하는 새로운 방법이 개발되었다.
   • 편향 계산: 주간 라인 데이터를 사용한다. 크기와 위상 값은 특정 기간(통상 한 달) 동안 평균화된다. 각 라인에 대해 가우시안 분포가 생성되며, 주파수 범위(10 Hz ~ 10 kHz)에 걸쳐 선형 보간을 수행하기 위해 무작위 점 집합이 선택된다. 이를 통해 $h_{raw}$를 $h_{unbias}$로 보정하는 데 사용되는 주파수 종속적 편향 추정치를 얻는다.
   • 불확실성 계산: 영구 라인을 적용하여 언바이어스된 변형률($h_{unbias}$)에 대해 수행된다. 이 접근 방식은 단기적인 변동(예: 정렬 또는 광학 응답 변화 < 1시간)에서 발생하는 계통적 불확실성을 포착한다.
   • 총 불확실성: 재구성 불확실성은 교정 불확실성(크기에 대한 보수적인 0.6% PCal 불확실성 및 0.5% 액추에이터 모델 불확실성, 그리고 위상에 대한 5 mrad 포함)과 이차적으로 결합된다. 50 Hz 및 150 Hz 근처의 추가적인 오차(전원 공급 장치 고조파)도 포함된다.

주요 결과

• 정밀도: 10 Hz ~ 2 kHz 대역에서 재구성된 변형률은 크기에서 2%, 위상에서 30 mrad의 정밀도를 달행한다.
• 언바이싱 성능: 온라인 언바이싱 과정은 10 Hz ~ 1 kHz 구간에서 잔류 편향을 **크기 기준 < 2%**로, 30 ~ 500 Hz 구간에서 **< 1%**로 감소시킨다. 잔류 위상 편향은 10 ~ 1 kHz 구간에서 < 15 mrad이다.
• 안정성: 재구성 편향은 시간 경과에 따라 매우 안정적인 것으로 나타났으며, 16개월 동안 단 두 번의 업데이트만을 필요로 했다.
• 전체 불확실성: 총 변형률 재구성 불확실성은 10 Hz ~ 1 kHz 대역에서 크기 기준 2.5% 미만이다. 위상의 경우, 200 Hz 미만에서는 25 mrad 이내를 유지하며, 1 kHz에서는 100 mrad까지 증가한다.

의의 및 주장
본 논문의 주요 의의는 Virgo 검출기에 대한 온라인 언바이싱의 구현과 주파수 종속적 불확실성 계산 방법의 개발에 있다고 주장한다.

• 온라인 교정이 수행되지 않았던 O3 런과 달리, O4는 저지연 분석 파이프라인으로 배포되기 전 약 10초의 지연 시간을 가진 언바이어스된 $h(t)$ 변형률을 특징으로 한다.
• 이 방법론은 주파수 독립적인 불확실성 추정치를 10 Hz ~ 2 kHz 대역에 걸친 보수적인 주파수 종속적 재구성 불확실성으로 대체함으로써 이전의 O3 관행을 개선한다.
• 이 시스템은 잔류 편향 및 불확실성을 월 단위로 모니터링할 수 있게 하여, 데이터 품질을 추적하고 품질이 낮은 기간(변형률 $h(t)$를 0으로 설정하는 기간)을 식별할 수 있도록 한다.
• 본 연구는 정확한 중력파 소스 특성화를 위해 필수적인 고충실도 변형률 데이터를 LIGO-Virgo-KAGRA(LVK) 협력단에 제공하기 위한 견고한 프레임워크를 구축한다.