[논문 리뷰] MARS Code Developments

N. Mokhov, S. Striganov외 4명|ArXiv.org|1998. 12. 15.

Particle physics theoretical and experimental studies참고 문헌 14인용 수 32

한 줄 요약

이 논문은 1 MeV–100 TeV 에너지 범위에서 복잡한 매질 내 강입자 및 전자기 캐스케이드를 시뮬레이션하기 위한 MARS 몽테카를로 코드에서의 주요 발전을 제시한다. 신규 핵 반응단면적 라이브러리, 파이온 생성 및 뮤온 흡수 모델의 향상, 전자기 상호작용 처리의 개선을 포함한다. 업데이트된 mars13(98) 코드는 광범위한 에너지 범위에서 실험 데이터와 강한 일치를 보이며, 가속기 차폐 및 검출기 설계 응용 분야에서의 신뢰성을 입증한다.

ABSTRACT

Recent developments in the physical model of 1 MeV to 100 TeV hadron and lepton interactions with nuclei and atoms are described. These include a new nuclear cross section library, a model for soft pion production, the cascade-exciton model, the dual parton model, deuteron-nucleus and neutrino-nucleus interaction models, detailed description of muon, pion, and antiproton absorption and a unified treatment of muon and charged hadron electromagnetic interactions with matter. New algorithms are implemented into the MARS13(98)} Monte Carlo code and benchmarked against experimental data. The code capabilities to simulate cascades and generate a variety of results in complex media have been also enhanced.

연구 동기 및 목표

• 1 MeV–100 TeV의 넓은 에너지 범위에서 강입자 및 렙톤의 물질과의 상호작용 물리 모델링을 향상시키기 위해.
• 새로운 집합 및 이론적 프레임워크를 사용하여 강입자-핵 및 광자-핵 상호작용에 대한 핵 반응단면적의 정확도를 향상시키기 위해.
• 고에너지 강입자 상호작용에서 소프트 파이온 생성, 캐스케이드-에크스아론, 이중 파르톤 역학에 대한 고도화된 모델을 개발하고 구현하기 위해.
• 뮤온 및 고전하 강입자에 대한 전자기 상호작용 처리를 통합하고 개선하여 에너지 손실 및 브레머스트랄링을 포함하기 위해.
• 복잡한 공학 응용 분야를 위한 스코어링, 기하 구조 처리, 그리고 MCNP 및 ANSYS와 같은 다른 시뮬레이션 도구와의 통합 능력을 확장하기 위해.

제안 방법

• 1 MeV에서 100 TeV까지의 강입자(p, n, π±, K±, p̄)에 대해 캐스케이드-에크스아론 모델(Cascade-Exciton Model, CEM)의 개선된 매개변수화 및 데이터 집합을 사용하여 새로운 핵 반응단면적 라이브러리를 구현하였다.
• 고에너지에서의 총, 비탄성, 탄성 강입자-핵 단면적을 계산하기 위해, 핵자 밀도 분포(Fermi형)를 사용한 Glauber 다중산란 이론을 적용하였다.
• 입사 입자 및 탈진 입자 운동량에 대해 약하게 의존하는 전단면적 비율 함수 R^{pA→π±X}를 사용한 새로운 현상학적 모델을 도입하여 pA 상호작용에서의 파이온 생성을 기술하였다.
• 정확한 공식과 1% 이내에서 총 단면적 및 dE/dx를 일치시키는 단순화된 알고리즘으로 복잡한 뮤온 브레머스트랄링 형식을 대체하였다.
• Bhabha 공식을 통한 δ 전자 에너지 샘플링 및 제한된 에너지 손실 계산을 통한 전자기 상호작용 모델링을 향상시켰다.
• 저에너지 중성자 및 광자 운반을 위한 자동 MCNP 인터페이스를 통합하였으며, 유연한 스코어링을 위한 객체 지향 기하 구조 및 메쉬 기반 히스토그램화를 도입하였다.

실험 결과

연구 질문

• RQ1 1 MeV–100 TeV 에너지 범위에서 강입자-핵 총 및 미분 단면적의 정확도를 어떻게 향상시킬 수 있는가?
• RQ2 pA 상호작용에서의 새로운 현상학적 파이온 생성 모델이 전방 및 후방 반구에서 실험 데이터를 얼마나 잘 재현하는가?
• RQ3 뮤온 및 고전하 강입자에 대한 업데이트된 전자기 상호작용 모델이 실험적 에너지 손실 및 각도 분포 데이터와 비교해 어떻게 성능을 보이는가?
• RQ4 새로운 MARS13(98) 코드가 개선된 스코어링 및 물질 조성 처리 능력을 통해 복잡하고 이질적인 매질 내 캐스케이드를 신뢰성 있게 시뮬레이션할 수 있는가?
• RQ5 neutron 감쇠 및 선량 계산을 포함한 기준 시나리오에서 MARS13(98)의 성능이 GEANT, MCNP, FLUKA 등의 기존 코드들과 비교해 어떻게 되는가?

주요 결과

• pA 상호작용에서의 새로운 파이온 생성 모델은 전체 운동량 범위에서 실험 데이터와 양호한 일치를 보이며, 특히 파이온 운동량 0.5–5 GeV/c의 도전적인 영역에서 뛰어난 성능을 나타낸다.
• 8 GeV/c의 프로톤 비틀림에서 두꺼운 lead 타겟에서의 파이온 스펙트럼을 MARS13(98)으로 시뮬레이션한 결과, 낮은 운동량 영역에서 GEANT보다 우수한 일치를 보였다.
• 50 cm 및 150 cm 두께의 리튬 흡수체를 통과한 뒤의 뮤온 에너지 스펙트럼과 각도 분포는 GEANT 및 ICool의 결과와 매우 유사하여, 새로운 전자기 상호작용 모델의 타당성을 검증하였다.
• 단순화된 뮤온 브레머스트랄링 알고리즘은 총 단면적과 dE/dx 모두에서 정확한 공식과 1% 이내에서 일치하여, 낮은 계산 비용으로도 높은 정확도를 확보하였다.
• SATIF-4 데이터와의 비교에서 MARS13(98)은 콘crete에서 중성자 감쇠(100–400 MeV) 및 조직 등가 필라멘트에서의 흡수 선량(0.1–10 GeV)에 대해 HETC, ANISN, MCNP 등 다른 코드들과 우수한 일치를 보였다.
• 자동 MCNP 인터페이스를 통해 저에너지 중성자 및 광자 운반을 정확하게 시뮬레이션할 수 있었으며, 새로운 물질 처리 방식은 동질화가 필요 없어져 시뮬레이션의 정밀도를 향상시켰다.