[논문 리뷰] Perturbed distribution functions with accurate action estimates for the Galactic disc

H. Al Kazwini, Q. Agobert외 9명|arXiv (Cornell University)|2022. 02. 01.

Stellar, planetary, and galactic studies참고 문헌 37인용 수 4

한 줄 요약

이 논문은 은하판에서 왜소된 분포함수를 모델링하기 위한 정교화된 이론적 프레임워크를 제시한다. 여기서 에피시클릭 근사법을 AGAMA를 통해 계산된 정확한 작용-각도 변수로 대체함으로써 3차원 위상공간에서의 공명 정밀 식별이 가능해진다. 또한 성장하는 바 또는 일시적인 나선 패턴을 모델링하기 위해 시간에 따라 변화하는 외부 힘을 도입하여, 고도에 따라 속도 공간에서 공명이 이동함을 밝혀내었으며, 이는 Gaia DR3 데이터를 이용한 은하의 3차원 잠재력에 대한 새로운 제약 조건을 제공한다.

ABSTRACT

In the <i>Gaia<i/> era, understanding the effects of the perturbations of the Galactic disc is of major importance in the context of dynamical modelling. In this theoretical paper we extend previous work in which, making use of the epicyclic approximation, the linearized Boltzmann equation had been used to explicitly compute, away from resonances, the perturbed distribution function of a Galactic thin-disc population in the presence of a non-axisymmetric perturbation of constant amplitude. Here we improve this theoretical framework in two distinct ways in the new code that we present. First, we use better estimates for the action-angle variables away from quasi-circular orbits, computed from the AGAMA software, and we present an efficient routine to numerically re-express any perturbing potential in these coordinates with a typical accuracy at the per cent level. The use of more accurate action estimates allows us to identify resonances such as the outer 1:1 bar resonance at higher azimuthal velocities than the outer Lindblad resonance, and to extend our previous theoretical results well above the Galactic plane, where we explicitly show how they differ from the epicyclic approximation. In particular, the displacement of resonances in velocity space as a function of height can in principle constrain the 3D structure of the Galactic potential. Second, we allow the perturbation to be time dependent, thereby allowing us to model the effect of transient spiral arms or a growing bar. The theoretical framework and tools presented here will be useful for a thorough analytical dynamical modelling of the complex velocity distribution of disc stars as measured by past and upcoming <i>Gaia<i/> data releases.

연구 동기 및 목표

• 비원형 궤도 및 은하판 상단에서 높은 수직 진동수를 보이는 경우에 대해 에피시클릭 근사법을 초월한 작용-각도 변수 추정의 정확도를 향상시키기 위해.
• 성장하는 바나 일시적인 나선 패턴과 같은 시간에 따라 변화하는 외부 힘을 포함한 선형 외부 힘 이론을 확장하기 위해.
• 특히 은하 평면 위아래 영역에서 3차원 속도 공간에서의 공명을 정밀하게 식별하기 위해.
• Gaia 관측된 항성 속도 분포를 분석하기 위한 분석적 천체역학 모델링을 위한 계산 효율적인 프레임워크를 제공하기 위해.
• 공명 영역에서 보다 고급 방법과 선형 외부 힘 이론을 융합하여 향후 자율적이고 일관된 모델링을 위한 기초를 마련하기 위해.

제안 방법

• 연구는 AGAMA 소프트웨어를 사용하여 원형 궤도에서 벗어나고 은하판 상단에서 높은 고도에 위치한 궤도에 대해 정확한 작용-각도 변수를 계산한다.
• 어느 외부 힘의 잠재력도 새로운 작용-각도 기저에서 약 1%의 정확도로 표현할 수 있는 수치적 루틴을 개발한다.
• 새로운 좌표계에서 선형화된 볼츠만 방정식을 적용하여 왜소된 분포함수를 계산하고, 이를 통해 에피시클릭 근사법을 대체한다.
• 성장하는 바의 진폭 변화를 모델링하기 위해 시간에 따라 변화하는 성장 함수를 도입한다.
• 이 방법은 새로운 작용-각도 기저에서 외부 힘의 푸리에 전개를 가능하게 하여 체계적인 공명 분석이 가능하다.
• 결과는 에피시클릭 근사법과의 비교 및 성장 외부 힘에 따른 분포함수의 시간적 진화 분석을 통해 검증된다.

실험 결과

연구 질문

• RQ1 정확한 작용-각도 변수를 에피시클릭 근사법 대신 사용할 경우, 외부 리드블라드 공명과 1:1 바 공명과 같은 공명이 고도가 높아질수록 속도 공간에서 어떻게 이동하는가?
• RQ2 은하 잠재력의 3차원적 구조는 항성 속도 공간에서 공명 특성의 공간적 및 운동학적 분포에 어느 정도의 영향을 미치는가?
• RQ3 성장하는 바와 같은 외부 힘의 시간적 진화는 작용-속도 공간에서 왜소된 분포함수의 형태와 구조에 어떤 영향을 미치는가?
• RQ4 선형 외부 힘 이론을 사용할 때, 외부 힘의 시간적 진화와 디스크의 정적 반응을 분리할 수 있는가?
• RQ5 속도 복원에 Stäckel fudge를 사용할 경우, 더 정밀한 방법에 비해 결과 속도장의 정확도에 어떤 영향을 미치는가?

주요 결과

• AGAMA 기반의 작용-각도 변수 사용으로 인해 속도 공간에서 공명의 위치가 이동한다: 고도가 높을수록 외부 1:1 바 공명이 외부 리드블라드 공명보다 낮은 순환 속도에서 나타난다.
• 공명은 평면에서 높이 1킬로파섹 증가할 때마다, 공명점에서의 속도 이동이 각각 코로레이션 8 km s⁻¹ kpc⁻¹, 외부 리드블라드 공명 6 km s⁻¹ kpc⁻¹, 1:1 공명 4 km s⁻¹ kpc⁻¹의 속도 이동을 보인다.
• 고도에 따라 속도 공간에서 공명이 이동하는 현상은 은하 잠재력의 3차원적 구조에 대한 직접적인 관측적 제약 조건을 제공한다.
• 해석적 진폭 진화를 가진 성장 바에 대해선 선형 응답이 매우 빠르게 진전되며, 바가 최종 진폭의 절반에 도달할 무렵에는 정적 경우와 거의 유사한 형태를 띤다.
• 시간에 따라 변화하는 처리 방식은 성장 중인 바와 일정 진폭을 가진 완성된 바 사이의 차이가 특히 패턴 속도 변화가 없을 경우 속도 분포상에서 항상 명확하지 않음을 보여준다.
• 현재 접근법의 한계로는 속도 복원에 Stäckel fudge를 사용하는 것과 속도 축에 대한 근사화가 필요로 하며, 향후 딥러닝(예: ACTIONFINDER)과 더 정밀한 변환 방법을 통해 향상 가능성이 있음을 밝혔다.