1 of 24

Heavy-Ion Collision Study at Intermediate Energies with DJBUU

2026.08.24 8th International Workshop on Nuclear Dynamics in Heavy-ion Reactions (IWND2026)

Kyungil Kim (IRIS, IBS)

2 of 24

Bird’s-eye-view of RAON

Guest House

Waste Storage

HQ Office

Control Center

Utility Bd.

Electricity Bd.

Assembly Bd.

SRF Test Bd.

Cryogenics Bd.

Low Energy B

Low Energy A

High Energy B

High Energy

SCL2

Accelerator

SCL3

IF

ISOL

● Accelerator System

● RI production System

● Conventional Utilities

● Experimental System

Campus Area: 952,066 m2

2

3 of 24

Motivation – Heavy-ion Collisions

3

Initial condition

Detection in experiments

How can we understand the process of heavy-ion collisions?

Ex.) How can we know the impact parameter of an event?

129Xe+120Sn, Ebeam = 65 MeV/n, b=1 fm

  • Transport models (BUU, QMD, …)
  • Hydrodynamics
  • Time-dependent Hatree-Fock models
  • Multifragmentation models
  • Statistical models

4 of 24

DJBUU

4

DaeJeon Boltzmann-Uehling-Uhlenbeck (DJBUU)

  • Lagrangian (QHD)
  • Phase-space density
  • Equation of motion

Ref.) M.Kim et al. PRC 101 (2020) 064614

5 of 24

DJBUU – Initialization

5

Relativistic Thomas-Fermi model calculation with same Lagrangian

Check the stability of an initialized nuclei

  • 197Au
  • for 300 fm/c

Comparison of the density evolution in a heavy-ion collision

  • 197Au + 197Au
  • Beam energy 100 MeV/n
  • for 140 fm/c

6 of 24

DJBUU – Pion Production

6

Ref.) M.Kim et al. Universe 8 (2022) 564

NN->N𝚫 𝚫 -> N𝛑

Production of pion in DJBUU

Isospin-dependent in-medium cross-sections

7 of 24

Motivation

7

  • Bubble nucleus ?

Nucleus which has a much smaller central density than the maximum density

In many nuclei, these bubble structures are observed in the proton density.

Since these nuclei have different density distributions, total and isospin densities at the collision center would be changed. We expect that these changes make the differences in the momentum distribution and isospin-dependent observables.

  • Chiral symmetry restoration at saturation density

8 of 24

Bubble Nuclei – Initialization

8

Density as a function of radius for 206Hg from RCHB and RTF calculations

  • We use two density configurations for 206Hg for comparison.
  • The density distribution from the RCHB calculation has a bubble configuration.
  • The DF values for total and proton densities are 13% and 18%, respectively.

Baryon and proton density distributions that are used in DJBUU calculations

  • Because of the width of the test particle (here, 4.2fm), with the same configuration obtained from the RCHB calculation, it is hard to get the bubble configuration. So, we modify the configuration.
  • For similar reasons, it is difficult to obtain small bubble nuclei with the test particle method which uses finite width.

9 of 24

Bubble Nuclei – Initialization

9

10 of 24

Bubble Nuclei – v1 (Directed Flow)

10

  • v1 of free protons

  • We choose test particles that are in effectively low density as free nucleons.

  • The differences are observed at high pT, but still within the error bar.

  • It is hard to reduce the error at high pT since the number of free nucleons is small in this region.

11 of 24

Bubble Nuclei – v1 (Directed Flow)

11

12 of 24

Bubble Nuclei – v2 (Elliptic Flow)

12

13 of 24

DJBUU with Various Models

13

  • DJBUU with the extended Parity Doublet Model (J. Zhou et al. Symmetry 2025, 17, 2155)
  • DJBUU with the Quark-Meson Coupling Model (D. I. Kim et al. PRC 113, 024615 (2026))

14 of 24

DJBUU – Quark-Meson Coupling (QMC) Model

14

  • Density evolution in the collision center

<Effective mass>

<QMC Lagrangian>

15 of 24

DJBUU – QMC Model

15

  • Directed flow (v1) comparison with experiment

<Directed flow>

<Transverse flow>

16 of 24

DJBUU – QMC Model

16

  • Delta baryon and pion production
  • Pion production ratios

17 of 24

DJBUU – Parity-Doublet Model (PDM)

17

  • Density evolution in the collision center

<Effective mass>

18 of 24

DJBUU – PDM Model

18

  • Effective masses with PDM model

Sn108+Sn112

Sn132+Sn124

19 of 24

DJBUU – PDM Model

19

  • DJBUU with the extended Parity Doublet Model (J. Zhou et al. Symmetry 2025, 17, 2155)

20 of 24

Deformed Nuclei and Flow Observables

20

N.M.Fortier et al. PRC111 (2025) L011901

The deformation of nuclei causes elliptic flows at ultra-central collisions.

21 of 24

Deformed Nuclei and Flow Observables

21

P. Danielewics et al. Science 298, 1592

22 of 24

Deformed Nuclei

22

23 of 24

Summary

23

  • Heavy-ion collisions at intermediate beam energies, from tens to hundreds of MeV per nucleon, offer a wide range of opportunities to explore interesting physics.
  • The DaeJeon Boltzmann-Uehling-Uhlenbeck (DJBUU) model has been developed to simulate and analyze these heavy-ion collisions at RAON.
  • Since the DJBUU model is based on relativistic mean-field calculations, it provides an opportunity to investigate fundamental physics, such as chiral symmetry, by adopting various mean-field models.
  • Bubble configurations in nuclei lead to differences in the directed flow, but the effect is marginal.
  • We are interested in the impact of nuclear deformation on heavy-ion collisions, and flow observables provide a useful tool for investigating these effects.

24 of 24

Thank you for your attention!