shinchonproton.yonsei.ac.kr/pt_files/2020/shinchon_shlim... · 2020-04-20 · new results on...
TRANSCRIPT
![Page 1: SHINCHONproton.yonsei.ac.kr/pt_files/2020/SHINCHON_shlim... · 2020-04-20 · New results on Bàμflow • Newresults on charm and bottom muon flow (arXiv:2003.03565) • Non-zero](https://reader034.vdocument.in/reader034/viewer/2022042207/5ea98260df517c34df555e58/html5/thumbnails/1.jpg)
SHINCHON
Simulation for Heavy IoN Collision with Heavy-quark and ONia
Sanghoon LimPNU
2020/Apr/16
![Page 2: SHINCHONproton.yonsei.ac.kr/pt_files/2020/SHINCHON_shlim... · 2020-04-20 · New results on Bàμflow • Newresults on charm and bottom muon flow (arXiv:2003.03565) • Non-zero](https://reader034.vdocument.in/reader034/viewer/2022042207/5ea98260df517c34df555e58/html5/thumbnails/2.jpg)
MC simulation for heavy quarks
• Heavy quarks are mostly produced in the early stage of HI collisions– Inside full evolution of the QGP
• There are several theoretical frameworks for HQ considering radiative and collisional energy loss– Successfully describe D-meson RAA and v2 simultaneously in HI collisions
• MC for heavy quarks will be useful to compare with the data more differentially – Angular correlation, jet sub-structure, and others
2JHEP 10 (2018) 174
![Page 3: SHINCHONproton.yonsei.ac.kr/pt_files/2020/SHINCHON_shlim... · 2020-04-20 · New results on Bàμflow • Newresults on charm and bottom muon flow (arXiv:2003.03565) • Non-zero](https://reader034.vdocument.in/reader034/viewer/2022042207/5ea98260df517c34df555e58/html5/thumbnails/3.jpg)
New results on Bàμ flow
• New results on charm and bottom muon flow (arXiv:2003.03565)• Non-zero v2 for bottom muon• Smaller bottom muon v2 than charm muon v2 at lower pT but similar at higher pT
3Can be covered by ALICE
![Page 4: SHINCHONproton.yonsei.ac.kr/pt_files/2020/SHINCHON_shlim... · 2020-04-20 · New results on Bàμflow • Newresults on charm and bottom muon flow (arXiv:2003.03565) • Non-zero](https://reader034.vdocument.in/reader034/viewer/2022042207/5ea98260df517c34df555e58/html5/thumbnails/4.jpg)
MC simulation for heavy quarks
• Heavy quarks are mostly produced in the early stage of HI collisions– Inside full evolution of the QGP
• There are several theoretical frameworks for HQ considering radiative and collisional energy loss– Successfully describe D-meson RAA and v2 simultaneously in HI collisions
• MC for heavy quarks will be useful to compare with the data more differentially – Angular correlation, jet sub-structure, and others
• More things to understand in small collision systems…
4
Phys. Rev. Lett. 121, 082301 (2018)
JHEP 12 (2019) 092
![Page 5: SHINCHONproton.yonsei.ac.kr/pt_files/2020/SHINCHON_shlim... · 2020-04-20 · New results on Bàμflow • Newresults on charm and bottom muon flow (arXiv:2003.03565) • Non-zero](https://reader034.vdocument.in/reader034/viewer/2022042207/5ea98260df517c34df555e58/html5/thumbnails/5.jpg)
General procedure (e.g. JETSCAPE)
5
https://indico.cern.ch/event/792436/contributions/3570550/
Initial geometryHydro evolution
Medium responseTranslation or MC HadronizationHadronic
interaction
![Page 6: SHINCHONproton.yonsei.ac.kr/pt_files/2020/SHINCHON_shlim... · 2020-04-20 · New results on Bàμflow • Newresults on charm and bottom muon flow (arXiv:2003.03565) • Non-zero](https://reader034.vdocument.in/reader034/viewer/2022042207/5ea98260df517c34df555e58/html5/thumbnails/6.jpg)
General procedure (e.g. JETSCAPE)
6
https://indico.cern.ch/event/792436/contributions/3570550/
![Page 7: SHINCHONproton.yonsei.ac.kr/pt_files/2020/SHINCHON_shlim... · 2020-04-20 · New results on Bàμflow • Newresults on charm and bottom muon flow (arXiv:2003.03565) • Non-zero](https://reader034.vdocument.in/reader034/viewer/2022042207/5ea98260df517c34df555e58/html5/thumbnails/7.jpg)
Initial condition
• MC Glauber– 2D distribution of smeared wounded nucleons in transverse direction– For hydrodynamics simulation, should convert into energy density by multiplying a scale factor
• The scale factor is tuned to match the measured charged particle multiplicity
• IP-Glasma– Gluon field starting from MC Glauber initial condition– Time evolution considering weak coupling of gluons
• This is possibly related to the bulk viscosity
– Requires another set of scale factor to matched the measured charged particle multiplicity
7
MC Glauber Pb+Pb 6-8 fm MC Glauber p+Pb 0-2 fm IP-Glasma Au+Au 6-8 fm
![Page 8: SHINCHONproton.yonsei.ac.kr/pt_files/2020/SHINCHON_shlim... · 2020-04-20 · New results on Bàμflow • Newresults on charm and bottom muon flow (arXiv:2003.03565) • Non-zero](https://reader034.vdocument.in/reader034/viewer/2022042207/5ea98260df517c34df555e58/html5/thumbnails/8.jpg)
Hydrodynamic simulation
• SONIC– 2+1D hydrodynamic simulation– https://sites.google.com/site/revihy/home– Option for pre-hydrodynamic evolution based on AdS/CFT (called superSONIC)
– Interfaced with the B3D hadronic cascade code
8
Pb+Pb 6-8 fm p+Pb 0-2 fm
![Page 9: SHINCHONproton.yonsei.ac.kr/pt_files/2020/SHINCHON_shlim... · 2020-04-20 · New results on Bàμflow • Newresults on charm and bottom muon flow (arXiv:2003.03565) • Non-zero](https://reader034.vdocument.in/reader034/viewer/2022042207/5ea98260df517c34df555e58/html5/thumbnails/9.jpg)
Hydrodynamic simulation
• SONIC– 2+1D hydrodynamic simulation– https://sites.google.com/site/revihy/home– Option for pre-hydrodynamic evolution based on AdS/CFT (called superSONIC)
– Interfaced with the B3D hadronic cascade code
• MUSIC– 2+1D or 3+1D hydrodynamic simulation– http://www.physics.mcgill.ca/music/– Recommended to use IP-Glasma initial condition
• From the recent update, there is no longer an option with MC Glauber
– Hydrodynamic evolution works well but need additional work for validation
– Will be useful to test bulk viscosity effect and longitudinal expansion
9
![Page 10: SHINCHONproton.yonsei.ac.kr/pt_files/2020/SHINCHON_shlim... · 2020-04-20 · New results on Bàμflow • Newresults on charm and bottom muon flow (arXiv:2003.03565) • Non-zero](https://reader034.vdocument.in/reader034/viewer/2022042207/5ea98260df517c34df555e58/html5/thumbnails/10.jpg)
HQ diffusion
• Reference– Phys. Rev. C 90, 024911 (2014)
• Initial HQ distribution– Spatial distribution from the MC Glauber initial condition– pT distribution from FONLL/PYTHIA
• Medium response – Follow time-dependent temperature/energy density profile from hydrodynamic simulation– Drag and diffusion based on Langevin equation
• No coalescence
10
Pb+Pb 6-8 fm Charm quark
dp(t)
dt= �⌘D(p)p(t) + ⇠(t)
<latexit sha1_base64="hKDH3o/4i2B7Qdu98xcMx6zRbe8=">AAACKHicbVDLSgMxFM3UV62vUZdugkVoEctMEXQjFnThsoJ9QKeUTCbThmYeJHfEMsznuPFX3Igo0q1fYvpY2OqBwOGcc8m9x40FV2BZYyO3srq2vpHfLGxt7+zumfsHTRUlkrIGjUQk2y5RTPCQNYCDYO1YMhK4grXc4c3Ebz0yqXgUPsAoZt2A9EPuc0pASz3z2vEloannBAQGrp/GWQnKWepBhq/wmcOA9NLbrBSX8UICn2LniWvSM4tWxZoC/yX2nBTRHPWe+e54EU0CFgIVRKmObcXQTYkETgXLCk6iWEzokPRZR9OQBEx10+mhGT7Riof9SOoXAp6qvydSEig1ClydnGyrlr2J+J/XScC/7KY8jBNgIZ195CcCQ4QnrWGPS0ZBjDQhVHK9K6YDopsD3W1Bl2Avn/yXNKsV26rY9+fFWnVeRx4doWNUQja6QDV0h+qogSh6Rq/oA30aL8ab8WWMZ9GcMZ85RAswvn8AP8mlVA==</latexit><latexit sha1_base64="hKDH3o/4i2B7Qdu98xcMx6zRbe8=">AAACKHicbVDLSgMxFM3UV62vUZdugkVoEctMEXQjFnThsoJ9QKeUTCbThmYeJHfEMsznuPFX3Igo0q1fYvpY2OqBwOGcc8m9x40FV2BZYyO3srq2vpHfLGxt7+zumfsHTRUlkrIGjUQk2y5RTPCQNYCDYO1YMhK4grXc4c3Ebz0yqXgUPsAoZt2A9EPuc0pASz3z2vEloannBAQGrp/GWQnKWepBhq/wmcOA9NLbrBSX8UICn2LniWvSM4tWxZoC/yX2nBTRHPWe+e54EU0CFgIVRKmObcXQTYkETgXLCk6iWEzokPRZR9OQBEx10+mhGT7Riof9SOoXAp6qvydSEig1ClydnGyrlr2J+J/XScC/7KY8jBNgIZ195CcCQ4QnrWGPS0ZBjDQhVHK9K6YDopsD3W1Bl2Avn/yXNKsV26rY9+fFWnVeRx4doWNUQja6QDV0h+qogSh6Rq/oA30aL8ab8WWMZ9GcMZ85RAswvn8AP8mlVA==</latexit><latexit sha1_base64="hKDH3o/4i2B7Qdu98xcMx6zRbe8=">AAACKHicbVDLSgMxFM3UV62vUZdugkVoEctMEXQjFnThsoJ9QKeUTCbThmYeJHfEMsznuPFX3Igo0q1fYvpY2OqBwOGcc8m9x40FV2BZYyO3srq2vpHfLGxt7+zumfsHTRUlkrIGjUQk2y5RTPCQNYCDYO1YMhK4grXc4c3Ebz0yqXgUPsAoZt2A9EPuc0pASz3z2vEloannBAQGrp/GWQnKWepBhq/wmcOA9NLbrBSX8UICn2LniWvSM4tWxZoC/yX2nBTRHPWe+e54EU0CFgIVRKmObcXQTYkETgXLCk6iWEzokPRZR9OQBEx10+mhGT7Riof9SOoXAp6qvydSEig1ClydnGyrlr2J+J/XScC/7KY8jBNgIZ195CcCQ4QnrWGPS0ZBjDQhVHK9K6YDopsD3W1Bl2Avn/yXNKsV26rY9+fFWnVeRx4doWNUQja6QDV0h+qogSh6Rq/oA30aL8ab8WWMZ9GcMZ85RAswvn8AP8mlVA==</latexit><latexit sha1_base64="hKDH3o/4i2B7Qdu98xcMx6zRbe8=">AAACKHicbVDLSgMxFM3UV62vUZdugkVoEctMEXQjFnThsoJ9QKeUTCbThmYeJHfEMsznuPFX3Igo0q1fYvpY2OqBwOGcc8m9x40FV2BZYyO3srq2vpHfLGxt7+zumfsHTRUlkrIGjUQk2y5RTPCQNYCDYO1YMhK4grXc4c3Ebz0yqXgUPsAoZt2A9EPuc0pASz3z2vEloannBAQGrp/GWQnKWepBhq/wmcOA9NLbrBSX8UICn2LniWvSM4tWxZoC/yX2nBTRHPWe+e54EU0CFgIVRKmObcXQTYkETgXLCk6iWEzokPRZR9OQBEx10+mhGT7Riof9SOoXAp6qvydSEig1ClydnGyrlr2J+J/XScC/7KY8jBNgIZ195CcCQ4QnrWGPS0ZBjDQhVHK9K6YDopsD3W1Bl2Avn/yXNKsV26rY9+fFWnVeRx4doWNUQja6QDV0h+qogSh6Rq/oA30aL8ab8WWMZ9GcMZ85RAswvn8AP8mlVA==</latexit>
![Page 11: SHINCHONproton.yonsei.ac.kr/pt_files/2020/SHINCHON_shlim... · 2020-04-20 · New results on Bàμflow • Newresults on charm and bottom muon flow (arXiv:2003.03565) • Non-zero](https://reader034.vdocument.in/reader034/viewer/2022042207/5ea98260df517c34df555e58/html5/thumbnails/11.jpg)
HQ diffusion (charm quark)
11
Pb+Pb 6-8 fm
p+Pb 0-2 fm
Charm quark
Charm quark
![Page 12: SHINCHONproton.yonsei.ac.kr/pt_files/2020/SHINCHON_shlim... · 2020-04-20 · New results on Bàμflow • Newresults on charm and bottom muon flow (arXiv:2003.03565) • Non-zero](https://reader034.vdocument.in/reader034/viewer/2022042207/5ea98260df517c34df555e58/html5/thumbnails/12.jpg)
HQ diffusion (charm vs bottom)
12
Bottom quark
Charm quark
Pb+Pb 6-8 fm
Pb+Pb 6-8 fm
![Page 13: SHINCHONproton.yonsei.ac.kr/pt_files/2020/SHINCHON_shlim... · 2020-04-20 · New results on Bàμflow • Newresults on charm and bottom muon flow (arXiv:2003.03565) • Non-zero](https://reader034.vdocument.in/reader034/viewer/2022042207/5ea98260df517c34df555e58/html5/thumbnails/13.jpg)
SHINCHON School in 2020/Jan
• SHINCHON School at PNU– ~15 participants to learn basic software/code for MC simulation
MC Glauber, IP-Glasma, SONIC, MUSIC, HQ diffusion– indico page:
https://n-ext.inha.ac.kr/event/382/– Instruction slides:
https://docs.google.com/presentation/d/1qxnRDLZ_hsXlv39KJ8kSWfxR7fegsIeYogIX0ZzH6Qk/edit?usp=sharing– Git repository
https://github.com/sh-lim/SHINCHON– Slack workspace
mcshinchon.slack.com
13
![Page 14: SHINCHONproton.yonsei.ac.kr/pt_files/2020/SHINCHON_shlim... · 2020-04-20 · New results on Bàμflow • Newresults on charm and bottom muon flow (arXiv:2003.03565) • Non-zero](https://reader034.vdocument.in/reader034/viewer/2022042207/5ea98260df517c34df555e58/html5/thumbnails/14.jpg)
Discussion with J.H. Hong
• Currently, the theoretical calculation procedure is not suitable for even-by-event MC
14
Phys. Lett. B 801 (2020) 135147
![Page 15: SHINCHONproton.yonsei.ac.kr/pt_files/2020/SHINCHON_shlim... · 2020-04-20 · New results on Bàμflow • Newresults on charm and bottom muon flow (arXiv:2003.03565) • Non-zero](https://reader034.vdocument.in/reader034/viewer/2022042207/5ea98260df517c34df555e58/html5/thumbnails/15.jpg)
Discussion with J.H. Hong
• Currently, the theoretical calculation procedure is not suitable for even-by-event MC
15
Phys. Lett. B 801 (2020) 135147
![Page 16: SHINCHONproton.yonsei.ac.kr/pt_files/2020/SHINCHON_shlim... · 2020-04-20 · New results on Bàμflow • Newresults on charm and bottom muon flow (arXiv:2003.03565) • Non-zero](https://reader034.vdocument.in/reader034/viewer/2022042207/5ea98260df517c34df555e58/html5/thumbnails/16.jpg)
Discussion with J.H. Hong
• For initial study, plan to implement – Time-dependent temperature profile from 2+1D hydrodynamic simulation– Time-dependent b-quark dN/dpT from Langevin framework– Study in different centrality of different initial geometry, medium property, and b-quark density
– 3+1D with MUSIC?
16
From lecture slides at the school Phys. Lett. B 801 (2020) 135147
![Page 17: SHINCHONproton.yonsei.ac.kr/pt_files/2020/SHINCHON_shlim... · 2020-04-20 · New results on Bàμflow • Newresults on charm and bottom muon flow (arXiv:2003.03565) • Non-zero](https://reader034.vdocument.in/reader034/viewer/2022042207/5ea98260df517c34df555e58/html5/thumbnails/17.jpg)
Another approach
• Reproduce the nuclear modification factor in MC (medium from hydro simulation)– Sample temperature dependent thermal width based on hydro simulation– Calculate disassociation rate (RAA) at different time
17
![Page 18: SHINCHONproton.yonsei.ac.kr/pt_files/2020/SHINCHON_shlim... · 2020-04-20 · New results on Bàμflow • Newresults on charm and bottom muon flow (arXiv:2003.03565) • Non-zero](https://reader034.vdocument.in/reader034/viewer/2022042207/5ea98260df517c34df555e58/html5/thumbnails/18.jpg)
Current status
• Set up the general MC framework in the KISTI computing farm for ALICE– Tested sources have been added into the github repository
• Bi-weekly meeting for code development– Wed. 2:00 pm– Participants: 박재범(고려대), 김준이(전북대), 서진주(인하대)
• Working on tuning the MC configuration for Pb+Pb collision at the LHC energies
18scaler from MC Glauber to SONIC input Npart dependence in SONIC
Phys. Rev. Lett. 116, 222302 (2016)
Done by Junlee Kim