베이즈로 푸는 μSR 자기 모멘트 추정법

초록

이 논문은 제로 필드 μSR 실험에서 뮤온 위치가 알려지지 않은 상황에서도, 계산된 쌍극자장 분포 f(ν/μ)를 이용해 베이즈 정리를 적용함으로써 관측된 프리시전 주파수 ν로부터 자성 이온의 순간 μ에 대한 사후 확률 g(μ|ν)를 추정하는 방법을 제시한다. 실제로 저자들은 저자기 모멘트를 갖는 오스마이트, 분자 자석, 그리고 철-비소계 초전도체 등 세 가지 물질에 적용해 유용성을 입증하였다.

상세 분석

본 연구는 μSR 실험에서 가장 큰 불확실성인 뮤온 정착 위치 문제를 통계적·물리적 모델링으로 해결하려는 시도이다. 먼저, 뮤온이 경험하는 로컬 자기장 B_loc은 외부장, 쌍극자장, 로렌츠장, 탈자장, 그리고 초미세 접촉 초극자장으로 구성되지만, 저자는 항상성(antiferromagnet)와 같은 경우 로렌츠·탈자장을 무시하고, 초극자장 역시 계산이 어려워 제외한다. 핵심은 쌍극자장 B_dip이 뮤온 위치 r_μ와 각 자성 이온의 순간 m_i에 의해 결정된다는 점이다. 이를 수식 (2)·(3) 형태로 표현하고, 무한 격자에 대해 Ewald 합이나 구형 절단법을 이용해 수치적으로 계산한다.

그 다음, 뮤온이 전하 밀도가 낮은 영역에 정착한다는 물리적 제약을 두고, 단위 셀 내에서 무작위로 수천·수만 개의 후보 위치를 생성한다. 각 위치마다 B_dip을 계산하고, γ_μ/2π = 135.5 MHz T⁻¹ 를 이용해 프리시전 주파수 ν = (γ_μ/2π)|B_dip| 로 변환한다. 이렇게 얻은 ν 값들의 히스토그램을 정규화하면 f(ν/μ)라는 확률밀도함수가 된다. 여기서 ν는 μ에 비례하므로, f(ν/μ) 자체가 μ에 대한 조건부 확률 P(ν|μ)와 직접 연결된다.

베이즈 정리 P(μ|ν) ∝ P(μ) P(ν|μ) 를 적용하는데, 저자는 μ에 대한 사전분포를 0 ~ μ_max 구간에서 균일하게 가정한다. 따라서 사후분포는 식 (11)와 같이
g(μ|ν) =