전력법칙 디스크를 포함한 일반화 광중력 체르미니크‑유사 문제의 평형점과 선형 안정성

초록

본 논문은 반지름이 a·b인 원판이 전력법칙(ρ∝r³)으로 질량을 분포하고, 회전 각속도 n>1을 갖는 일반화된 광중력 체르미니크‑유사 제한 삼체 문제를 다룬다. 방사압, 소원시의 편평도, 디스크 질량의 영향을 포함한 운동 방정식을 도출하고, 평형점(5개 기존 라그랑주점 + 1개 새로운 점)과 제로‑속도 곡면을 구한다. 선형 안정성 분석을 통해 L₂·L₃는 디스크 내·외반경에 따라 안정될 수 있음을, 모든 공선점은 불안정, L₄는 Routh 임계 질량비 이하에서 조건부 안정임을 확인한다.

상세 분석

논문은 먼저 G(m₁+m₂)=1, 거리 단위=1, 시간 단위=ω⁻¹으로 정규화한 좌표계에서, 큰 1차원은 복사압을 받는 방사원, 작은 2차원은 편평한 타원체(편평도 A₂)로 가정한다. 디스크는 ρ(r)=c r³ 형태의 전력법칙을 따르며, 내부반경 a와 외부반경 b를 갖는다. 디스크의 중력 퍼텐셜 V와 방사압·편평도·디스크 중력에 의한 교정항을 포함한 유효 퍼텐셜 Ω를 정의하고, 코리올리 항을 포함한 운동방정식(1)-(2)을 얻는다. Jacobi 상수 C=−(ẋ²+ẏ²)+2Ω를 도출하여 제로‑속도 곡면을 구한다.

평형점은 Ωₓ=Ω_y=0을 만족하는 (x,y)로 정의한다. y=0인 공선점은 K(x)=0이라는 10차 다항식 형태로 전개되며, x 구간을 (−∞,−µ), (−µ,0), (0,1−µ), (1−µ,∞) 로 나누어 각각의 절대값 처리와 부호 변화를 고려한다. µ가 작은 경우, ρ를 µ¹⁄⁴의 멱급수로 전개하여 근사해를 구하고, 각 구간마다 계수 B_i, C_i 등을 통해 고차 방정식(13), (17), (20), (23) 등을 얻는다. 이러한 전개는 기존 라그랑주점 L₁, L₂, L₃와 새로운 점 Lₙ(두 원점 사이) 존재 조건을 명시한다.

선형 안정성은 작은 섭동 ξ,η에 대해 2차 미분 방정식 ξ¨−2nη˙=Ωₓₓ ξ+Ωₓ_y η, η¨+2nξ˙=Ω_yₓ ξ+Ω_y_y η 로 전환하고, 행렬식 |λ²I−A|=0을 풀어 특성다항식 λ⁴+(4n²−Ωₓₓ−Ω_y_y)λ²+ (ΩₓₓΩ_y_y−Ωₓ_y²)=0을 얻는다. 판별식 Δ>0와 실근의 부호를 통해 안정성을 판단한다. 결과적으로 L₂와 L₃는 디스크의 내·외반경 a, b가 특정 구간에 있을 때 Δ>0와 λ²<0을 만족해 선형적으로 안정함을 보인다. 반면 L₁과 기존 공선점은 언제나 Δ<0이므로 불안정이다. 삼각점 L₄와 L₅는 µ<µ_Routh≈0.03852일 때만 조건부 안정성을 유지하고, 방사압 β와 편평도 A₂가 증가하면 안정 영역이 축소된다.

마지막으로, 파라미터 β, A₂, a, b의 변화를 수치적으로 시뮬레이션하여 제로‑속도 곡면과 평형점 위치가 어떻게 변하는지 시각화한다. 특히 디스크 질량이 클수록 n이 커져 라그랑주점이 원점에서 멀어지고, 새로운 평형점이 두 원점 사이에 나타나는 현상이 강조된다.