低SIS能区重离子碰撞中的非平衡动力学研究

同位旋效应对直接流消失的影响李庆峰合作者：沈彩万，郭琛琛，王永佳，李祝霞，J.Lukasik，W.Trautmann主要内容取自：PHYSICALREVIEWC83,044617(2011);SCIENCECHINAPhysics55252(2012);SCIENCECHINAPhysics即将投稿。2020/5/192012超重核合成与性质研讨会（兰州）QF-2主要内容动机模型简介及相关改进集体流、流消失现象及研究现状动力学过程对集体流和核阻止的影响INDRA能区直接流激发函数的理论描述流消失现象对对称能密度依赖程度的探测流消失的同位旋效应依赖的多系统研究小结及展望2020/5/192012超重核合成与性质研讨会（兰州）QF-3主要动机①希望扩展我们的输运模型UrQMD至低SIS能区，从而在一个最大的能量范围内（从SIS到LHC）研究核反应的动力学性质及所存在的可能的相变情况；②在低SIS能区，最近对一些观测量---如集体流和核阻止效应的---的激发函数的实验数据和初步理论研究的比较表明，现在的理论解释还存在较大的不自洽性，需要进一步完善。③INDRA和MSU两实验组对SIS能区的直接流消失能量的测定有差别。④INDRA实验给出的流消失能量是出射粒子种类依赖的。2020/5/192012超重核合成与性质研讨会（兰州）QF-4模型简介•UrQMD=Ultra-relativisticQuantumMolecularDynamics极端相对论量子分子动力学模型，是一个研究较大能量区域、较多反应系统的微观非平衡动力学模型；•新粒子产生通过两种方式获得：共振态衰变和弦成块（模拟新相的运动学）•粒子的相互作用考虑了平均场和两体碰撞；•具体更新情况（cascade模式的）请见•虽然叫极端相对论量子分子动力学模型，但是，在低能端，为了简化问题，其继承了QMD对平均场的处理方式（而不是RQMD的协变动力学处理方式），因此，在低能端，UrQMD除了一些技术上的特殊处理外，所涵盖的物理内容和QMD完全相同。我们利用UrQMD对LHC能区软物理研究的最新成果•研究由HBT两π介子关联所抽取的火球时空信息在LHC能区的的集体行为、膨胀行为和标度行为，具体请见：•PHYSICALREVIEWC85,034908(2012)•PHYSICALREVIEWC85,044901(2012)•J.Phys.G39(2012)0651012020/5/192012超重核合成与性质研讨会（兰州）QF-52020/5/192012超重核合成与性质研讨会（兰州）QF-6UrQMD用于低能运算的模型改进•势修正部分：①平均场势or②Skyrme势（王永佳同学将给出一些初步计算）•碰撞项部分：①Pauli-blocking阻塞的再处理；②两体碰撞截面的介质修正的仔细考虑。•这些处理成立的前提：初始核合理的束缚能及方均根半径在一个时间段（150fm/c）内的相对稳定。2020/5/192012超重核合成与性质研讨会（兰州）QF-7UrQMD中现在已经采用的平均场势项•除了通常的动能项以外，UrQMD内现已考虑的势项：I.已具备：II.新加入：介子库仑势部分也考虑了。JPG31_1359(2005);32_151(2006);36_015111(2009)2020/5/192012超重核合成与性质研讨会（兰州）QF-8Pauli-blocking阻塞的再处理•默认版本的UrQMD中为了加快运算速度，使用了一个等效的拟合公式用于PB阻止：•但是，此公式在低能区被证明不够准确：这会允许过多的碰撞发生。2020/5/192012超重核合成与性质研讨会（兰州）QF-9PB阻止续•因此，我们采用了传统的处理PB的方式，即满足如下两条件的某次碰撞的同类粒子就可以发生碰撞：①PB：②UU：自旋兼并因子两粒子相空间分布函数随机数碰撞阻塞因子时时满足2020/5/192012超重核合成与性质研讨会（兰州）QF-10两体碰撞截面的介质修正核子-核子两体弹性散射截面的密度、动量修正是众所周知的，许多理论组都对此作较为详细的研究，如：PRC69_017601(2004)我们最早研究了质量劈裂效应在截面中的影响PRC62_014606(2000)我们最早研究了截面密度修正中的同位旋依赖性。2020/5/192012超重核合成与性质研讨会（兰州）QF-11截面介质修正在输运模型中的实现•暂不考虑非弹部分：JPG32_407(2006)2020/5/192012超重核合成与性质研讨会（兰州）QF-12反应系统、末态处理、观测量在INDRA能区，广泛选择不同的碰撞系统（见结果部分）及碰撞参数,事件停止时间150fm/c，每一计算用10万以上事件分析以保证统计的可靠性。反应停止后，采用相空间判断的Coalescence模型的相空间判断，P0和R0参数分别为0.2GeV/c和2.8fm.•集体流：•核阻止参数vartl：集体流v1的产生及时间演化举例2020/5/192012超重核合成与性质研讨会（兰州）QF-13流消失现象2020/5/192012超重核合成与性质研讨会（兰州）QF-14流消失成因：吸引与排斥2020/5/192012超重核合成与性质研讨会（兰州）QF-15碰撞前后碰撞时MSU和INDRA的实验2020/5/192012超重核合成与性质研讨会（兰州）QF-16PRC,61,021602（R）（2000）PRC,65,044604（2002）INDRA：更好的反应平面确定，更好的粒子分辨流消失观测量2020/5/192012超重核合成与性质研讨会（兰州）QF-17directedtransversemomentumwithpx(i)beingthetransversemomentumoftheithparticleandwhereyc.m.(i)istheassociatedrapidityoftheithparticle.Intransversemomentum,allrapiditybinsaretakenintoaccount.“Inseveralexperiments,itwasshownthatthebalanceenergyEbalisnearlyindependentofthenatureoftheparticle.Thebalanceenergyisnearlythesamefornucleons,light,andheavyfragments.“一个典型例子：PHYSICALREVIEWC583494（1998）MSU实验-流消失的初始核同位旋依赖2020/5/192012超重核合成与性质研讨会（兰州）QF-18PRC78，1022（1997）INDRA实验-流消失的粒子种类依赖性2020/5/192012超重核合成与性质研讨会（兰州）QF-19Eur.Phys.J.A30,31(2006)2020/5/192012超重核合成与性质研讨会（兰州）QF-20改变初始化、势修正、碰撞项后的14个UrQMD测试版本P-B:Pauli-blockingEoS：状态方程Fu：截面密度依赖因子Fp：截面动量修正因子Ini.:初始化方式H-S：硬球模式W-S：Woods-Saxon整体而言，UrQMD中的HS模式比WS模式模拟的核要小一点点（2.5%）反应系统的密度时间演化图2020/5/192012超重核合成与性质研讨会（兰州）QF-21反应系统的同位旋不对称度的时间演化图2020/5/192012超重核合成与性质研讨会（兰州）QF-222020/5/192012超重核合成与性质研讨会（兰州）QF-23Pauli-Blocking效应对观测量的影响2020/5/192012超重核合成与性质研讨会（兰州）QF-24集体现象的产生原因分析I:初始化、平均场、碰撞项无碰撞无势修正无库仑、对称势HSWS2020/5/192012超重核合成与性质研讨会（兰州）QF-25集体现象的产生原因分析II:平均场和截面的动量依赖效应UrQMD-ISMFP2UrQMD-IISMFP1UrQMD-VIIISMNopNNlimitUrQMD-IXSFP1*均为FU12020/5/192012超重核合成与性质研讨会（兰州）QF-26核阻止此反应条件下实验值：0.800.032020/5/192012超重核合成与性质研讨会（兰州）QF-27自洽性检验如上对不同入射能量下碰撞的集体流和核阻止观测量的研究表明：UrQMD中不同贡献项对它们的影响顺序是相同的，它们之间具有高度的结果自洽性。因此，接下来的研究中我们将只选定其中一个观测量用于和实验值的比较。（注意：实验结果的自洽性检验也是十分重要的）2020/5/192012超重核合成与性质研讨会（兰州）QF-28截面中的密度修正和动量修正效应密度压低效应增强无碰撞项动量恢复效应增强两边修正均无动量依赖两边修正均无动量依赖2020/5/192012超重核合成与性质研讨会（兰州）QF-29直接流在INDRA能区的激发函数转变能低能核反应更细致考虑INDRA和FOPI实验结果有不同下面关注一个次级量：对称能2020/5/192012超重核合成与性质研讨会（兰州）QF-3018setsFromB.AlexBrown,PRL85,5296(2000)放射性束流设备现状Q.LiforAPCTPinPohang31CoolingStorageRing(CSR)FacilityatHIRFL/LanzhouinChina(2008)upto500MeV/Afor238UBeijingRadioactiveIonFacility(BRIF-II)atCIAEinChina(2012)RadioactiveIonBeamFactory(RIBF)atRIKENinJapan(2007)FacilityforAntiprotonandIonResearch(FAIR)/GSIinGermany(2016)upto2GeV/Afor132Sn(NUSTAR-NUclearSTructure,AstrophysicsandReactions)SPIRAL2/GANILinFrance(2013)SelectiveProductionofExoticSpecies(SPES)/INFNinItaly(2015)FacilityforRareIsotopeBeams(FRIB)/MSUinUSA(2018)upto400(200)MeV/Afor132SnTheKoreanRareIsotopeAccelerator(KoRIA)(2016?)upto250MeV/Afor132SnQ.LiforAPCTPinPohang32理论上研究IEOS的方法列表MicroscopicMany-BodyApproachesNon-relativisticBrueckner-Bethe-Goldstone(BBG)TheoryRelativisticDirac-Brueckner-Hartree-Fock(DBHF)approachSelf-ConsistentGreen’sFunction(SCGF)TheoryVariationalMany-Body(VMB)approachGreen’sFunctionMonteCarloCalculationVlowk+Renormaliz