다크스타스 코드 공개: 별 내부에서 WIMP 포획·소멸을 동시에 계산하는 최신 진화 패키지

초록

다크스타스(DarkStars)는 별 내부에서 WIMP(약한 상호작용 질량 입자)의 포획·소멸과 전통적인 별 진화 방정식을 동시에 풀어주는 공개 소프트웨어이다. 다양한 WIMP 속도 분포와 평균 자유행로에 기반한 복합 분포·전도 에너지 전달 모델을 포함해 현재 가장 포괄적인 WIMP 미세물리학을 제공한다. 코드 사용법과 입력 물리량을 간략히 소개하고, 은하 중심에서의 다크 스타 시뮬레이션 사례를 제시한다.

상세 분석

다크스타스는 별 내부에서 WIMP이 포획되고, 서로 만나 소멸하면서 발생하는 에너지 주입을 기존의 별 진화 방정식에 직접 결합한다는 점에서 혁신적이다. 기존 코드들은 보통 WIMP 포획률을 별도 후처리하거나, 단순한 에너지 주입 파라미터만 적용했지만, 다크스타스는 포획률을 별의 온도·밀도 프로파일과 연동된 미세한 속도 분포 함수(맥스웰, 비등방성, 스트림 등)로 계산한다. 또한, WIMP의 평균 자유행로(Mean Free Path, MFP)를 기준으로 두 가지 전도 모델을 자동 전환한다. MFP가 별의 규모보다 작을 경우 확산형 전도 방정식을, 크면 자유 입자 흐름 모델을 적용해 에너지 전달 효율을 정확히 추정한다. 이러한 복합 분포·전도 스킴은 별 중심부에서 WIMP이 축적될 때 발생하는 온도 상승과 압력 변화를 실시간으로 반영한다. 코드 구조는 MESA와 유사한 모듈식 설계로, 핵반응 네트워크, 방사전달, 대류 모델 등 기존 천체물리 모듈과 원활히 연동된다. 입력 파라미터에는 WIMP 질량, 단면적(스칼라·스핀-의존), 은하 중심의 밀도·속도 분포, 별의 초기 질량·조성 등이 포함되며, 사용자는 파라미터 파일을 수정해 다양한 시나리오를 손쉽게 탐색할 수 있다. 다크스타스는 또한 포획·소멸 균형을 자동으로 찾아주는 수렴 알고리즘을 제공해, 장시간 진화 시뮬레이션에서도 수치적 안정성을 유지한다. 논문에서 제시된 은하 중심 다크 스타 예시는, 높은 DM 밀도(∼10⁹ GeV cm⁻³)와 낮은 WIMP 속도 분산이 별의 핵융합을 억제하고, WIMP 소멸 에너지에 의해 팽창된 저밀도 구조를 형성함을 보여준다. 이는 기존 이론이 예측한 ‘다크 스타’의 특성과 일치하면서, 코드가 실제 물리적 현상을 정량적으로 재현할 수 있음을 입증한다. 전체적으로 다크스타스는 WIMP 물리와 별 진화의 복합 문제를 통합적으로 다루는 최초의 공개 패키지이며, 향후 DM 탐색, 은하 중심 천체물리, 초신성 전구 단계 연구 등에 광범위하게 활용될 전망이다.