クォーク・核物質研究部⾨報告
Quark-Nuclear Matters【拠点・施設】
ALICE/LHC, STAR/RHIC
J-PARC, RIKEN/RIBF
計算科学研究センター
宇宙史研究センター

・格⼦QCD 研究
・宇宙元素合成研究
・QGP/臨界点研究

TCHoU strearing committee meeting, QNP report, 18/Mar/2022

ShinIchi Esumi, Quark-Nuclear Matters, TCHoU

1【代表的な論⽂】

Lattice QCD group

•
••

QCD•
•
•
!

NJLモデルTricritical point

2nd【講演、賞など】•! =0!!

1st"

Figure 1. Schematic view of expected phase diagram of the NJL model on the T -µ plane. Solidand broken curves represent the first- and second-order phase transitions, respectively. Closed circle
denotes the tricritical point where the first-order phase transition line terminates.

“Application of tensor renormalization group to Nambu-Jona-Lasinio model and Hubbard model”
Y. Kuramashi (招待講演),whereTensor
Networks in Many Body and Lattice Field (Online, July 26-30, 2021).
n = (n , n , n , n )(∈ Z ) specifies a position in the lattice Λ, with the lattice spacinga. χ(n) and χ̄(n) are Grassmann-valued
the Dirac structure.
Since Theories”
they“Application of tensor renormalizationgroupfields
towithoutQuantum
Field
describe the Kogut-Susskind fermions, χ(n) and χ̄(n) are single-component Grassmann
Y. Kuramashi (招待講演),variables.
DWQ@25，(Online,Dec.
18,
η (n) is the staggered sign function definedby η 13-Dec.
(n) = (−1)
with2022).
η (n) = 1. The partition function is defined in the ordinal manner:
“Tensor renormalization group approach! to$ (1+1)-dimensional
Hubbard model”
Z= dχ(n)dχ̄(n) e .
(2.2)
S. Akiyama, 熱場の量⼦論とその応⽤2021，(オンライン, Aug. 30-Sep. 1, 2021).
For vanishing mass m, eq.
(2.1) is invariant under the following continuous chiral transfor秋⼭進⼀郎(D3) 筑波⼤学学⽣表彰
2022年3⽉mation:12344

ν•JHEP01(2021)121! ≠0!!

•“Calculation of derivative of nucleon form factors in Nf=2+1 lattice QCD at Mπ=138 MeV on a (5.5 fm)3 volume”,
PACS Collaboration: K.-I. Ishikawa et al., Phys.Rev.D104(2021)074504“Tensor renormalization group approach to (1+1)-dimensional Hubbard model”,
S. Akiyama and Y. Kuramashi, Phys.Rev.D104(2021)014504“Metal-insulator transition in (2+1)-dimensional Hubbard model with tensor renormalization group”,
S. Akiyama, Y. Kuramashi and T. Yamashita, PTEP2022(2022)023“Phase transition of four-dimensional lattice ϕ4 theory with tensor renormalization group”,
S. Akiyama, Y. Kuramashi, and Y. Yoshimura, Phys.Rev.D104(2021)034507“Charm and beauty in the deconfined plasma from quenched lattice QCD”,
H.-T. Ding, O. Kaczmarek, A.-L. Lorenz, H. Ohno and H.-T. Shu, Phys.Rev.D104(2021)114508“Finite-size scaling around the critical point in the heavy quark region of QCD”,
A. Kiyohara, M. Kitazawa, S. Ejiri, and K. Kanaya, Phys.Rev.D104(2021)114509

νn1 +···+nν−11−Sn∈Λ•χ(n) → eiα#(n) χ(n),(2.3)χ̄(n) →χ̄(n)e4 talks in the 38th International Symposium on
Lattice
Field Theory(2.4)(Lattice 2021)，(Online, July 26-30, 2021).
with α ∈ R and $(n) = (−1).
•“Restoration of chiral symmetry
in cold anddense Nambu–-Jona-Lasinio model with tensor renormalization group”, Y. Kuramashi
2.2 Tensor network representation•“Tensor renormalization Wegroup
approach to (1+1)-dimensional Hubbard model”, S. Akiyama
introduce the tensor network representation for eq. (2.2) in a similar way with refs. [10,11]. Hereafter, we set a = 1 for simplicity. Firstly, we expand the local Boltzmann weights•“Critical endpoints in (2+1)and 4-flavor QCD with Wilson-Clover fermions”, H. Ohno•
“The lower moments of nucleon structure functions in lattice QCD with physical quark masses”, R. Tsuji for PACS Collaborationiα#(n)

n1 +n2 +n3 +n433See ref. [25] for a different TRG approach with the Kogut-Susskind fermion, where the TRG procedure
is applied to the Schwinger model after integrating out the fermion fields analytically.

TCHoU strearing committee meeting, QNP report,– 3 – 18/Mar/2022

ShinIchi Esumi, Quark-Nuclear Matters, TCHoU

2

宇宙元素合成の研究
Ni領域の質量測定実験
2018年

重元素の⽣成過程を探る
重元素合成仮説(Rプロセス)
110Sn, 80Zr  

等時性度%%((2018 12⽉2021年    ))Sn領域
2018 11⽉、2021 11⽉

Ni領域2018 11⽉、2021 4⽉
稀少RIリング@RIBF（理研）

TCHoU strearing committee meeting, QNP report, 18/Mar/2022

33

ker magnets system was successfully
T0125Ag
all of 20201) and extraction timing
ded from 100ns to 400ns, leading towhere is the velocity of124Pd
the particle of interest relative
ored events in a single extraction.
t were produced at RIBF by impingto the speed of the light in vacuum. Revolution time
ium beam on a 5-mm thick Be tarspectrum of all injected nuclei is shown in Fig. 1 (b) (deidentified at BigRIPS by energy loss
tails of determination of the revolution time in R3 can behamber (IC) placed at F3 and theirF) from F3 to F5. After injection
found in [19]). Since the isochronous condition of the ringorage for about 1 ms, equivalent tois optimized for the reference particle, T0 is independent
the particles were extracted. Fig. 1of the momentum. To determine the mass, the velocity
3 of all events produced at BigRIPS
Time-of-Flight F3 - F5 [ns]events from R3, which are indicatedneeds to be determined event-by-event from the timeles. These events were selected byof-flightalong
thebeamline
from F3 to S0 (T OF3S0 ) byFig. 1. Particle Identification(PID)
showingthe energy
TOF/dE gate method that is necloss in theIonizarion(IC) atF3 and ToF
usingtheChamber
following
equation,vidual injection method2) . Due tofrom F3 to F5. Extracted particles from R3 are shownuring the beam transport from Biin black empty circles.
circulation inside the ring, the total
Length3S0=.
(2)
low compared to previous experi(TOF+ T OFof f set )3S01, we show the preliminary numbers
Table 1. Preliminary total events produced at BigRIPS
s compared to the produced eventsand extracted event after the ring.
e PID shown in Fig. 1. The estiThe average path length from F3 to S0 (Length3S0 )bers were performed by a PID gate
Isotope # Produced # Extracted
128and the TOFof f set caused by the electronics and them F3 to F5 asTCHoU
well as a strearing
gate for en- committee meeting,
Sn
1,380
2
QNP
report, 18/Mar/2022
127
energyloss
on the beamline, are deIn
76,270in the detectors
108. Unfortunately, no events of 124 Pd
126125
Cd
137,059
wever, sufficient events ofAg andtermined
via Eq.(2) 266
by using known masses of 124 Ag125
Ag
23,81963
were extracted.
127
124andSn.803 The parametersthat could reproduce thePd0re, we will focus our e↵orts on solvknown m/q values are Length3S0 = 84.859(2) m andes encountered during the experin to continue mass measurements ofT OFof f set = 325.47(1) ns. The mass is then determined2) Y. Abe et al., RIKEN Accel. Prog. Rep. 52 (2019)as well as neutron-deficient nuclei.

粒⼦識別（今回の実験で稀少RIリングに⼊射できた粒⼦を⿊丸で表記）

forviaRep.
Eq.(1).
3) S. Naimi etal., each
RIKENevent
Accel. Prog.
52 (2019)Additional systematic uncer4) A. Ozawa tainties,
et al., RIKEN Accel.
(2019), Prog.
dueRep.
to 52thedetermination of parameterssyssuch as Length3S0 , T OFof f set and T0 were estimated

TABLE I. Mass excess from literature and the mass excessof nuclei measured in this work are shown in the second andthird column, respectively. Total uncertainties are shown aswell as thecontribution from the reference mass uncertainty
2018年11⽉に⾏ったSn領域m0 and the statistical uncertainty stat . The systematic uncertainty の実験の結果をPRLに投稿
sys is estimated from the uncertainty of T0 and thefit parameters Length3S0 and T OFof f set of Eq.(2).
Nucleus MEAM E20[keV]
126
125
123

InCd
Pd

MER3 total m0 stat sys[keV] [keV] [keV] [keV] [keV]-77809(4) -77707 269-73348.1(29) -73237 320-60430(790) -60282 265

254
252
248

65
192
86

6240401000

Δ ME(R3 - AME20) keV

Energy Loss at F3 IC [ch]TOF in the R3 was measured by the E-MCP detector atandthea plasticscintillator
asurement of neutron-rich nuclei S0with
Rare-RI
Ring detector placed at ELC after
the ejection from the R3. Another IC was installed at⇤1⇤4maguchi,⇤1 Y. Nagata,⇤2 D. Nagae,⇤2 Y. Abe,⇤1 F. Suzaki,
M. Wakasugi,
N. Kaname,
ELC,⇤1⇤3wherean additional
PID
was performed. Finally,⇤5
2021年11⽉:２度⽬のSn領域
Yano,⇤4 T. Moriguchi,⇤4 A. Ozawa,⇤4 M. Kanda,
Y. Koizumi,⇤5 K. Okubo,⇤5 M. Otsu,⇤5 K.
intheNaI scintillator detector⇤5⇤5⇤5⇤6 particles⇤6 were stopped⇤6⇤6ki, N. Shinozaki, T. Yamaguchi, T. Chillery, C. Hanai, Y. Hijitaka, J. Li, R.
⇤6⇤1 IC at ELC.
behindtheDozono,⇤6 S. Ota,⇤6⇤7 S. Michimasa,の質量測定実験
H. Baba,⇤1 H. placed
Suzuki,⇤1(実験データは解析中)
N. Fukuda,
Y. Shimizu,⇤1 H.
⇤1Takeda,⇤1 M. Yashimoto,⇤1 H.F. Li,⇤8 and T. TheUesaka,
mass-to-charge
ratio (m/q) of the particle of interest with a revolution time T is determined relative21, we conducted an experiment atto a reference particle with m0 /q0 and T0 by using the(R3) aiming to measure masses offollowing
formula [13, 18]:(Ag) and palladium (Pd) isotopes.
128Sn
ei were measured before at R3, how127Inv
125124es ofAg andPd showed large
m126Cd m0 T uu 1⇣ 2 ⌘inties that were suspected to be due=,(1)traction timing of the storage ring.
T 2
qq0 T0 1t al., RIKEN Accel. Prog. Rep. 54

ShinIchi Esumi, Quark-Nuclear Matters, TCHoU

AME20 uncertainty
R3500

0−500−1000

123Pd の質量誤差を改善
127Sn126In125Cd124Ag123

Pd

FIG. 2. Mass excess values of nuclei measured at R3 comparedto literature values from AME2020 [20]12 Esumi,
3 Quark-Nuclear
45ShinIchi
Matters,TCHoU

masses at A= 112 124 was thought to be due to theshell quenching at N = 82 [23]. However, it has been
recently demonstrated from mass measurement of 132 Cdthat this shell gap, although reduced, is not quenched
for Z < 50 [24]. Another suggested solution to this issue
was the increase of the two-neutron separation energies(S2n ) that might be associated with a sudden transition
from deformed shape to spherical shape for nuclei with4

1年延⻑

ALICE Run3 with LS2 upgrades

O2 : fast/continuous readout

TPC with GEM readout

TCHoU strearing committee meeting, QNP report, 18/Mar/2022

ShinIchi Esumi, Quark-Nuclear Matters, TCHoU

Partonic Energy Loss
Jet-quenching••

5

ALICE Run3&4 and future

ALICE internal-jetstructure
Jet broadening
or narrowing

FoCal-E pad 検出器の読出エレキHGCROC試験

ALICE-Focal
Tsukuba group• 東北⼤学 ELPH のGeVーγ 施設 (7⽉)
• CERN 研究所 SPS 加速器 (9⽉−10⽉）
ALICE 遠隔データ収集シフト (9⽉− 10⽉)
研究員 Abderrahmane Ghimouz⽒着任 (11⽉)

ALICE-3 : Letter of Intentfor Run5 and future
all silicon detectors
TCHoU strearing committee meeting, QNP report, 18/Mar/2022

ShinIchi Esumi, Quark-Nuclear Matters, TCHoU

6

Beam Energy Scan II programat RHIC-STAR (COL / FXT modes)

fluctuation

200GeV d+Au, 200GeV O+O during run21

polarization

directed flow

TCHoU strearing committee meeting, QNP report, 18/Mar/2022

elliptic flow

ShinIchi Esumi, Quark-Nuclear Matters, TCHoU

7ShinIchi Esumi, Quark-Nuclear Matters, TCHoU

8Isobar analysis for CME search and beyond

TCHoU strearing committee meeting, QNP report, 18/Mar/2022

STAR setup
Fixed Target modeeRHICfuture•••••

STAR-BES-II (and beyond) data analysis
including FXT mode setup
FAIR-CBM experiment preparation towards
high-density QCD phase-diagram
Neutron detector development for J-PARC-HI
with phenix TOF detector
HIMAC H355/H447 experiment with Riken-group
MRPC detector development for J-PARC E16

FAIR (ドイツ)
NICA (ロシア)HIAF (中国)
J-PARC-HI (⽇本)

TCHoU strearing committee meeting, QNP report, 18/Mar/2022

ShinIchi Esumi, Quark-Nuclear Matters, TCHoU

9外部資⾦・科研費など
2019 − 2023 基盤研究(S) 江⾓ 晋⼀
⾼次ゆらぎと粒⼦相関による⾼密度クォーク核物質の１次相転移と臨界点観測への挑戦
2020 − 2024 基盤研究(S) 中條 達也
LHC 超前⽅光⼦測定によるグルーオン飽和とQGP⽣成起源
2020 − 2024 基盤研究(S) ⼩沢 恭⼀郎（KEKクロスアポ）
⾼輝度陽⼦ビームによる原⼦核中での明確な中間⼦質量変化の実験的確⽴
2020 − 2023 基盤研究(C) 森⼝ 哲朗
中性⼦星の構造解明に向けた重い中性⼦過剰核のスキン厚測定
2021 − 2022 挑戦的研究(萌芽) 笹 公和
⻑半減期放射性セシウム135を加速器質量分析法により超⾼感度で検出する試み
2017 − 2022 伊藤科学振興会/物理学分野研究助成⾦ 江⾓ 晋⼀
原⼦核衝突ビームエネルギー⾛査による⾼密度クォーク核物質の研究

TCHoU strearing committee meeting, QNP report, 18/Mar/2022

ShinIchi Esumi, Quark-Nuclear Matters, TCHoU

10

TCHoU Workshop, Photon & Particle Detector
Mar 22 (Tue), 2022, 10:00 – 12:20

クォーク・核物質 まとめ
・格⼦QCD 研究
クォーク物質の熱⼒学的特性
１次相転移、臨界終点、プラズマ中の重クォーク

・

宇宙元素合成研究
不安定核⽤蓄積リング
陽⼦過剰領域実験
Ni,Sn領域質量測定、寿命測定

・

QGP/臨界点研究
LHC-ALICE実験（ジェット、光⼦、HF、FoCal検出器）
RHIC-STAR実験（エネルギー⾛査、臨界点、固定標的）
J-PARC-HI/E16実験（MRPC、中性⼦検出器開発）
FAIR-CBM実験（⾼密度QCD物質研究）
⾼温・⾼密度領域における将来実験計画

TCHoU strearing committee meeting, QNP report, 18/Mar/2022

TCHoU Workshop, Quark-Nuclear MattersMar 24 (Thr), 2022, 10:30 – 18:00

ShinIchi Esumi, Quark-Nuclear Matters, TCHoU

11