큐픽스 우주

#분자동역학

2건의 기사

천문·천체물리

응력이 유카와 결정 내 확산을 가속한다 — 중성자별 지각 크리프의 미시 메커니즘

결정화된 쿨롱(유카와) 플라스마에서 자기확산은 공공-격자간 원자 쌍(프렌켈 쌍)의 열적 활성화 생성과 이동을 통해 진행되는데, 연구진은 분자동역학으로 변형되지 않은 결정에서 밝혀진 아이링 법칙을 변형된 결정으로 확장했다. 전단된 유카와 결정 시뮬레이션 결과, 정적 탄성 변형은 프렌켈 쌍 생성에 참여하는 영역에 추가로 저장된 탄성 에너지 때문에 확산을 지수적으로 증가시키며, 이는 D(ε)=D₀(Γ*)exp[V*σ/k_BT] 형태로 기술된다(σ=½με²은 전단 변형 ε이 저장한 탄성 에너지 밀도, V*는 활성화 부피). 전형적인 천체물리학적 차폐 조건에서는 단일 활성화 부피 N_act = n_iV* ≃ 31로 충분하며 차폐가 강할수록 값이 줄어드는데, 이는 프렌켈 쌍 생성이 진동의 열적 최대점에서 최근접 이웃들의 짧은 협동만을 필요로 함을 시사한다. 응력에 의한 확산 증대는 변형된 중성자별 지각과 결정화 중인 백색왜성 중심부에서 일어나는 크리프의 미시물리적 메커니즘을 제공한다.

arXiv astro-ph.HE (고에너지)

천문·천체물리

따뜻한 중성자별 외각 지각의 상태방정식, 분자동역학 시뮬레이션으로 정밀화

연구진은 분자동역학(Molecular Dynamics) 기법의 미시적 다체 시뮬레이션으로 따뜻한 이온 플라즈마의 상태방정식(EoS)을 기술했으며, 전자 차폐와 이온을 유한 크기 가우스 분포로 모델링한 상호작용 퍼텐셜, 효율적인 에발트(Ewald) 에너지 합산 절차를 포함해 기존 계산을 개선했다. 이를 통해 중성자별 외각 지각 조건에서 바리온 수밀도 7.48×10⁻¹⁰~2.09×10⁻⁴ fm⁻³, 온도 k_BT 1~5 MeV 범위의 상태방정식 P(n_B,T)를 얻었다. 결과의 활용도를 높이기 위해 표로 정리한 데이터와 신경망 파라미터화를 Zenodo 저장소에 공개했다. 열적 단열지수 Γ_th 기반의 열적 유효 파라미터화로 기술할 경우, 중간 정도의 온도에서도 이온의 열적 효과가 내부 지각에 가까운 고밀도 영역에서 핵심적인 역할을 한다는 점을 확인했다.

arXiv astro-ph.HE (고에너지)