원자 중력계의 응답 함수와 시스템 해석

초록

본 논문은 원자 중력계를 선형 시불변(LTI) 시스템으로 모델링하고, 가속도·속도·변위에 대한 임펄스 응답 함수를 각각 삼각형, 메앤더, 디랙 콤 형태로 도출한다. 이를 이용해 중력 구배와 자기 인력 등 비균일 중력장의 영향을 정량화하고, 펌프-오실레이션을 적용한 장치가 1차 저역통과 필터와 동일함을 밝힌다.

상세 분석

논문은 먼저 3‑펄스 원자 중력계가 자유 낙하 시간 T에 따라 위상 차 Δφ = k g T² + φ₁−2φ₂+φ₃ 로 표현된다는 점을 강조한다. 여기서 k는 레이저 파수벡터이며, Δφ를 0으로 만드는 주파수 증가율 α를 측정함으로써 중력 가속도 g = 2π α/k 를 얻는다. 저자들은 이 관계를 시간에 따라 변하는 가속도 g(t)와의 선형 결합 형태로 재작성하고, 측정 결과가 입력 신호와 컨볼루션으로 연결되는 LTI 시스템으로 해석한다.

두 번의 적분 연산을 갖는 시스템의 임펄스 응답은 기본적으로 단위 계단 함수와 단위 램프 함수의 조합으로 나타난다. 구체적으로 가속도에 대한 가중 함수 w_g(t)는 0≤t≤T 구간에서는 (t/T²)·(2T−t) 형태의 삼각형 곡선이며, T≤t≤2T 구간에서는 (2T−t)/T²·t 형태로 대칭된다. 이는 가속도 입력이 시간 구간 전체에 걸쳐 평균화되는 효과를 의미한다.

속도와 변위에 대한 가중 함수는 각각 w_V(t)=−dw_g/dt, w_z(t)=d²w_g/dt² 로 정의되며, 앞서 도출된 w_g(t)의 1차·2차 미분으로 얻어진다. 속도 가중 함수는 메앤더 형태(양·음 구간이 교대로 나타나는 사각형 파형)이며, 변위 가중 함수는 디랙 콤(δ(t)−2δ(t−T)+δ(t−2T)) 형태로 세 점에서만 신호를 샘플링한다. 이러한 특성은 가속도·속도·변위 각각에 대한 시스템의 민감도가 어떻게 다르게 나타나는지를 명확히 보여준다.

시스템의 순간 응답을 이용해 비균일 중력장, 특히 상수 구배 γ와 선형 구배(γ₁+γ₂z)에 의한 오차를 정량화한다. 상수 구배의 경우, 측정된 중력은 g₀±γ·h_eff 로 표현되며, h_eff는 실제 궤적의 1/6(분수형 분출형) 또는 7/24(방출형) 높이에 해당한다. 이는 전통적인 코너 큐브 중력계와 달리 측정 고도가 궤적 상단보다 낮은 위치에 존재함을 의미한다. 선형 구배의 경우, 가속도 신호를 다항식으로 전개하고 가중 함수의 모멘트 C_n을 이용해 각 차수별 기여를 계산한다.

자기 인력과 같은 복잡한 교란은 고차 다항식 근사와 잔차 함수를 이용해 평가한다. 저자들은 12차까지의 모멘트를 사용해 −1.27 µGal 수준의 교정값을 얻었으며, 근사 오차가 최대 잔차의 절반 이하가 되도록 설계함으로써 교정 불확실성을 제한한다.

마지막으로, 블록 오실레이션을 이용해 원자들을 정지시킨 장치는 적분 구간이 T, T+T₀, 2T+T₀ 로 확장되며, 임펄스 응답이 거의 균일한 사각형 형태가 된다. 따라서 이 장치는 1차 저역통과 필터와 동등한 주파수 응답을 보이며, 기존 3‑펄스 장치는 2차 저역통과 필터, 코너 큐브 장치는 3차 저역통과 필터와 동일한 특성을 가진다. 이는 진동 억제 설계 시 필터 차수에 기반한 직관적인 해석을 가능하게 한다.

전반적으로 논문은 원자 중력계를 LTI 시스템으로 모델링함으로써 복잡한 시스템atics를 간단한 수학적 도구(임펄스 응답, 가중 함수, 모멘트)로 정량화하고, 다양한 설계(분출형, 방출형, 블록 오실레이션형)의 성능 차이를 명확히 제시한다.