Heavy-ion collision simulation with high performance computer

이 논문은 고성능 컴퓨터 (HPC) 를 활용하여 DJBUU 와 SQMD 모델을 기반으로 한 중이온 충돌 시뮬레이션에 대한 최근 예비 연구 결과를 보고합니다.

핵심

🌌 1. 왜 이런 실험을 할까요? (동기)

우리는 우주의 비밀을 풀기 위해 두 가지 방법을 쓰고 있습니다.

1. 천문학: 블랙홀이나 중성자별이 부딪힐 때 나오는 '중력파'를 관측합니다. (마치 우주의 진동을 듣는 것)
2. 지상 실험: 지상에서 무거운 원자핵들을 아주 빠르게 충돌시켜, 중성자별 내부처럼 압축된 핵물질을 만들어 봅니다.

한국에는 **RAON(라온)**이라는 거대한 가속기 시설이 지어지고 있습니다. 이곳에서는 평소에는 존재하지 않는 '희귀 동위원소'를 만들어 충돌 실험을 하려고 합니다. 이 실험 결과를 예측하고 이해하기 위해, 과학자들은 **초고속 컴퓨터 (HPC)**를 이용해 가상 실험을 하고 있는 것입니다.

🎮 2. 컴퓨터 시뮬레이션은 어떻게 작동할까요? (모델)

원자핵 충돌은 너무 복잡해서 손으로 계산할 수 없습니다. 그래서 과학자들은 두 가지 다른 '게임 엔진' 같은 모델을 개발했습니다.

• DJBUU (다제온 모델):

  • 비유: 마치 밀려다니는 구름을 보는 것과 같습니다.
  • 수만 개의 입자가 모여 있는 '밀도'와 '흐름'을 중점적으로 봅니다. 입자 하나하나를 세기보다는, 전체적인 유체 (물이나 공기) 의 흐름처럼 핵물질을 다룹니다.
  • 특징: 입자들이 서로 부딪히거나 퍼지는 과정을 '확률'과 '평균장'으로 계산합니다.

• SQMD (신동 모델):

  • 비유: 마치 수만 개의 공 (구슬) 이 서로 부딪히는 것을 보는 것과 같습니다.
  • 각 입자 (핵자) 를 하나의 작은 파동 덩어리로 보며, 하나하나의 운동 궤적을 추적합니다.
  • 특징: 입자들이 뭉쳐서 '조각 (Fragment)'이 되는 과정을 더 자세히 묘사합니다.

이 두 모델은 서로 다른 방식으로 세상을 바라보지만, 같은 현상 (원자핵 충돌) 을 시뮬레이션합니다.

🧪 3. 무엇을 발견했나요? (결과)

연구팀은 한국과학기술정보연구원 (KISTI) 의 슈퍼컴퓨터 '누리온'을 이용해 이 두 모델을 가동했습니다.

① 두 모델의 비교 (동일한 상황, 다른 결과)

• 안정된 원자핵 (납 + 칼슘 등) 을 충돌시켰을 때:
  • 에너지가 낮을 때는 두 모델이 거의 같은 결과를 냅니다. (두 게임 엔진이 같은 그림을 그리는 셈)
  • 하지만 에너지가 매우 높거나, 충돌 각도가 작을 때는 두 모델의 결과가 달라집니다. 이는 모델이 핵물질을 바라보는 '안경 (방정식)'이 다르기 때문입니다.
• 불안정한 원자핵 (나트륨 -20 등) 을 충돌시켰을 때:
  • 여기서 큰 차이가 나타났습니다. DJBUU 모델은 조각이 더 크게 뭉치는 반면, SQMD 모델은 더 작게 부서지는 것으로 나타났습니다.
  • 이유: 불안정한 원자핵은 마치 약한 구조의 성처럼 쉽게 무너집니다. 두 모델이 이 '약한 성'을 어떻게 초기에 설정하고, 충돌 전까지 어떻게 유지하는지에 따라 결과가 크게 달라진 것입니다.

② 새로운 이론 적용 (QMC 모델)

• 기존에 쓰던 이론 (QHD) 에 '쿼크 (입자의 구성 성분)'의 영향을 더한 새로운 이론 (QMC) 을 DJBUU 모델에 적용해 보았습니다.
• 결과: 충돌 순간의 밀도가 더 높아지는 것을 발견했습니다.
• 의미: 이는 핵물질이 얼마나 단단한지 (압축성) 를 알려주는 중요한 지표입니다. 밀도가 높으면 파이온 (입자) 이 더 많이 만들어질 수 있다는 뜻입니다.

🚀 4. 이 연구의 의미는 무엇인가요? (결론)

이 연구는 단순히 컴퓨터 게임을 하는 것이 아닙니다.

1. RAON 가속기의 나침반: 앞으로 RAON 에서 일어날 실험 결과를 미리 예측하고, 실험 설계에 도움을 줍니다. 특히 '불안정한 원자핵'을 다룰 때 어떤 모델이 더 정확한지 알려줍니다.
2. 우주 이해의 확장: 지상의 실험으로 알게 된 '핵물질의 성질'은 중성자별 내부의 비밀을 푸는 열쇠가 됩니다. 중력파 관측 (우주) 과 원자핵 충돌 실험 (지상) 이 서로의 결과를 검증하며 우주의 비밀을 밝혀내는 것입니다.

💡 한 줄 요약

**"슈퍼컴퓨터로 만든 두 가지 가상 실험실 (DJBUU, SQMD) 을 비교하며, 우주의 극한 상태 (중성자별) 를 이해하고, 한국형 거대 가속기 (RAON) 의 실험을 성공적으로 이끌기 위한 지도를 그리는 연구"**입니다.

이처럼 과학자들은 거대한 컴퓨터와 정교한 수학을 통해, 눈에 보이지 않는 우주의 가장 작은 입자부터 가장 거대한 별까지 연결하는 비밀을 풀어가고 있습니다.

기술 요약

논문 기술 요약: 고성능 컴퓨팅을 활용한 중이온 충돌 시뮬레이션 (DJBUU 및 SQMD 모델)

1. 문제 제기 (Motivation)

• 배경: 중이온 충돌 실험은 중성자별 내부와 같은 고밀도 핵물질의 특성을 이해하는 핵심 도구입니다. 최근 GW170817 중력파 관측과 NICER X-ray 관측을 통해 중성자별의 물리적 특성이 규명되면서, 천체물리학과 핵물리학의 교차점에서 중성자별 연구의 중요성이 부각되었습니다.
• 과제: 이러한 천체 관측 결과를 보완하고 검증하기 위해 지상에서의 중이온 충돌 실험이 필수적입니다. 특히 한국에서 건설 중인 희귀 동위원소 가속기 (RAON) 를 통해 불안정 동위원소 빔을 이용한 실험이 예정되어 있습니다.
• 난제: 중이온 충돌 시스템은 매우 복잡하며, 이를 이론적으로 모델링하고 실험 결과를 예측하기 위해서는 방대한 계산 자원이 필요합니다. 기존 모델들의 한계를 극복하고 RAON 실험을 정확히 예측하기 위해 고성능 컴퓨팅 (HPC) 을 활용한 정밀 시뮬레이션이 요구됩니다.

2. 방법론 (Methodology)

연구진은 한국과학기술정보연구원 (KISTI) 의 고성능 컴퓨팅 시스템 (NURION) 을 활용하여 두 가지 수송 모델 (Transport Models) 을 개발 및 적용했습니다.

• 사용된 모델:

  1. DJBUU (DaeJeon Boltzmann-Uehling-Uhlenbeck):
     • 상대론적 BUU 방정식을 기반으로 하며, 1-체 분포 함수의 역학을 다룹니다.
     • 핵심 특징: 테스트 입자 (Test Particle) 방법을 사용하며, 기존 가우스 또는 삼각형 프로파일 대신 다항식 (Polynomial) 프로파일을 사용하여 위상 공간 밀도를 계산합니다. 이는 유한한 끝점을 가지며 부드러운 형태를 제공하여 수치적 안정성을 높입니다.
     • 상호작용: 상대론적 평균장 (RMF) 이론을 기반으로 하며, 기존 QHD(Quantum Hadrodynamics) 모델 대신 쿼크 - 메손 결합 (QMC) 모델을 도입하여 Bag 상수 (Bag constant) 의 효과를 고려한 유효 핵자 질량을 계산합니다.
  2. SQMD (Sindong Quantum Molecular Dynamics):
     • QMD 접근법을 사용하여 n-체 역학을 다룹니다.
     • 핵심 특징: 핵자를 고정된 폭을 가진 가우스 파동 패킷으로 취급하며, 위상 공간 밀도는 Wigner 변환을 통해 유도됩니다.
     • 상호작용: 밀도 의존적 Skyrme 퍼텐셜을 사용하며, 충돌은 탄성 산란 단면적을 기반으로 처리됩니다.
     • 클러스터링: 최소 스패닝 트리 (Minimal Spanning Tree) 알고리즘을 사용하여 파편화 및 클러스터링을 분석합니다.

• 컴퓨팅 환경:

  • KISTI 의 NURION 시스템을 활용하여 대규모 시뮬레이션을 수행했습니다.
  • 예시: $^{196}\text{Au} + ^{196}\text{Au}$ 충돌 시 100 개의 테스트 입자, 700 개의 시간 단계, 100 만 개의 격자 셀 등을 처리하는 고부하 계산을 수행했습니다.

3. 주요 기여 (Key Contributions)

• RAON 실험 대응 모델 개발: 불안정 동위원소 빔을 다루기 위해 DJBUU 와 SQMD 모델을 고도화하고, KISTI 의 HPC 자원을 활용한 대규모 병렬 시뮬레이션 파이프라인을 구축했습니다.
• QMC 모델의 DJBUU 적용: 기존 DJBUU 가 사용하던 QHD 모델에서 Quark-Meson Coupling (QMC) 모델로 전환하여, 핵자의 유효 질량과 압축성 (Incompressibility, $K_0$) 을 더 정교하게 묘사할 수 있도록 모델을 개선했습니다.
• 모델 간 비교 분석 프레임워크: 두 가지 서로 다른 수송 모델 (BUU 형식 vs QMD 형식) 을 동일한 조건에서 비교하여, 고밀도 핵물질 연구에서의 모델 의존성을 체계적으로 분석했습니다.

4. 결과 (Results)

• 모델 간 비교 ($^{208}\text{Pb} + ^{40,48}\text{Ca}$):
  • 낮은 에너지 (50 AMeV): 두 모델 모두 가장 큰 파편 (Biggest Fragments, BFs) 의 크기가 유사하게 생성되었습니다.
  • 높은 에너지 (100 AMeV): 에너지가 높아지고 충돌 파라미터가 작아질수록 모델 간 차이가 두드러졌습니다. 이는 상태 방정식 (EOS) 과 모델의 고유한 역학적 차이 때문입니다.
• 불안정 동위원소 시뮬레이션 ($^{20}\text{Na} + ^{208}\text{Pb}$):
  • 불안정 핵 ($^{20}\text{Na}$) 을 표적으로 한 시뮬레이션에서 모델 간 편차가 더욱 커졌습니다. DJBUU 는 SQMD 에 비해 약 30% 더 큰 BFs 를 생성했습니다.
  • 원인 분석: 충돌 전 불안정 핵의 안정성 차이에서 기인했습니다. DJBUU 는 중심 밀도가 더 조밀해지는 경향을 보인 반면, SQMD 는 밀도가 감소하고 퍼지는 경향을 보였습니다. 이는 불안정 핵의 초기 조건 설정 및 진화 방식의 차이에서 비롯된 것으로 판단됩니다.
• QMC 모델 도입 효과:
  • DJBUU+QMC 모델은 기존 DJBUU+QHD 모델에 비해 충돌 중심부에서 더 높은 최대 핵자 밀도를 보였습니다.
  • 이는 핵물질의 상태 방정식에 민감한 관측량 (예: 파이온 생성률, 직접 흐름 등) 에 유의미한 영향을 미칠 것으로 예상됩니다.

5. 의의 및 향후 전망 (Significance & Future Work)

• RAON 실험 예측 능력 향상: 이번 연구는 RAON 에서 수행될 희귀 동위원소 빔 중이온 충돌 실험의 결과를 예측하고 해석하는 데 필수적인 도구로 작용할 것입니다. 특히 불안정 핵의 거동을 더 정확하게 묘사하기 위해 표면 항을 포함한 메손 장 방정식 솔버와 시간 의존적 파동 패킷 폭 도입 등을 통해 모델을 지속적으로 개선할 계획입니다.
• 다중 메신저 천체물리학과의 연계: 중이온 충돌 실험을 통해 얻은 고밀도 핵물질의 상태 방정식 (EOS) 정보는 중성자별 내부 구조 이해 및 중력파 관측 데이터 해석에 중요한 제약 조건을 제공합니다.
• 고성능 컴퓨팅의 활용: 복잡한 핵물리 현상을 규명하기 위해 HPC 를 활용한 대규모 시뮬레이션의 중요성을 입증했으며, 이는 향후 더 정밀한 핵물질 연구의 기반이 될 것입니다.

이 논문은 이론적 모델링, 고성능 컴퓨팅, 그리고 실험 물리학을 결합하여 고밀도 핵물질의 성질을 규명하려는 한국 연구진의 선구적인 노력을 보여줍니다.