1原子核和强相互作用物质的相变张三,李四(湖南科技大学建筑与城乡规划学院,湘潭411201)摘要:简要回顾原子核和强相互作用物质的相结构及相变研究的现状。说明原子核和强相互作用物质的相结构和相变的研究是原子核物理、粒子物理、天体物理、宇宙学和统计物理等领域共同关心重要前沿领域,到目前为止已取得重大进展,但无论是具体实际问题还是研究方法等方面都需要系统深入的研究。关键词:原子核物理;强相互作用物质;相与相变1引言100年前,爱因斯坦通过分析充满空腔的辐射系统的熵与充满空腔的气体系统的熵,提出电磁辐射由光量子组成[1,2],从而建立了光子的概念,吹响了引导人们探索微观世界的冲锋号。进一步的深入研究表明,组成物质世界的粒子可以分为强子和轻子两类,粒子间相互作用可分为引力作用、电磁作用、弱作用和强作用4类。参与强相互作用的粒子或具有强相互作用的系统统称为强相互作用物质(包括强子物质、夸克物质等)及其特殊形式——原子核(由有限个强子组成的系统),对原子核和强相互作用系统的相结构及相变的研究,对于认识强相互作用系统的相结构、相变,了解宇宙的起源和演化至关重要,并且可能是有限系统的统计物理的检验平台。因此,近年来关于原子核和强相互作用系统的相变的研究不仅是原子核物理、天体物理、宇宙学及粒子物理等领域研究的重要前沿课题,还引起了有限量子多体系统领域和统计物理学界的极大关注。本文简要介绍原子核及强相互作用系统的相及相变研究的现状。2原子核的相及相变2.1原子核的单粒子运动与集体运动原子核是有限数目的强子组成的束缚系统,其中的核子(质子和中子)具有单粒子运动,并建立壳模型成功的描述原子核的相应性质。实验上对原子核的能谱和电磁跃迁等的研究表明,原子核还具有整体运动,并建立了原子核具有形状和振动、转动等集体运动模式的概念。人们通常利用将核半径按球谐函数展开来描述原子核的形状,并将相应的形变称为极形变(如图1所示)。已经观测到和已经预言的原子核形状多种多样[3,4],比较重要的是四极形变,实验上已经观测到的最高极形变是16极形变[3,4]。按照壳模型和集体模型的观点,幻数核多为球形,而偏离满壳的核形变则为形变核,形变核可以细分为长椭球形、扁椭球形、三轴不对称形、梨形、香蕉形、纺锤形等。这里是个图。注意:一个图必须在一页上画,不能分在两页上!图1原子核的极形变的形状示意图近年来的研究表明,在较高激发能和较高角动量情况下,原子核的集体能谱消失,即出现带终结现象[5],这表明发生了由集体运动到单粒子运动的相变。小二号黑体粗四号楷体五号宋体五号黑体五号楷体五号楷体五号黑体四号宋体加粗(下同)五号宋体(下同)小四号宋体(下同),加粗若有图,应写这样的说明图题,小五宋体加粗(下同)参考文献的标注—顺序码(下同)注意页边距:上为2.50厘米,下为2.00厘米左边2.8厘米右边2.0厘米22.2原子核的形状相变原子核形状的研究一直是原子核结构理论中一个重要的问题,这是因为原子核形状与原子核组成成分及其两种运动形式--集体运动和单粒子运动、中子质子比、角动量、激发态能量和核环境的温度等都密切相关。例如,集体模型中计算单粒子运动时常用的变形平均势就和核形状有关,不同形状原子核的集体运动模式各不相同[6];同时原子核的形状由所有核子的空间分布决定,而且随集体运动模式的不同而变化[7]。另一方面,原子核的形状和一定的动力学对称性相联系[4],核形状变化与原子核的动力学对称性的破缺相联系。原子核的形状发生变化表明其状态和性质发生了变化,即发生了相变。因此,原子核的形状相结构和相变的研究是原子核结构研究的重要内容。由于形状共存可能是单粒子运动和集体运动较强耦合的结果[7],因此形状共存也是核形状研究中关注的焦点[8]。早期对于原子核形状相变的研究大多集中在一系列同位素或同中子素的基态[4,9],基态核的形状相变普遍存在于各个质量区[3],近年来关于临界状态对称性和三相点的研究[10-16]以及对超重核的形变和形状共存的研究[17],极大地丰富了基态和形状相变的研究内容。另一方面,由于实验上g-射线探测器阵列技术的进步,使得我们不仅可以对原子核基态的形状进行研究,而且可以对激发态、尤其是高自旋态的核形状进行研究。激发态核的形变则更富含物理内容,如超形变带、回弯现象、同核异能态等都和形变直接相关;2003年观测到的沿Yrast带出现的集体振动模式到定轴转动模式的变化表明低激发态中可能存在转动(或角动量)驱动的由球形(振动)到长椭球形(定轴转动)的形状相变[18]。对原子核基态形状的研究通常采用的理论模型有集体模型[6]、相互作用玻色子模型(IBM)[4]、HFB方法[19],另外还可以使用热力学统计理论[20]。而对于原子核激发态的形状的研究则采用Landau相变理论[21]、有限温度推转HFB[22]、推转IBM[23]等。在这些方法中,集体模型有比较直观的几何图象,但是缺乏微观机制;而微观理论没有直接的几何图象。由于IBM既有较好的微观基础[24],又可以由相干态理论建立直观的几何图象[4],所以IBM理论在原子核的形状相变研究中得到了广泛的应用。早期利用IBM对原子核基态的形状相变的研究可以归纳为Casten三角形[4],近年来Iachello利用几何模型对原子核基态形状相变的研究将Casten三角形扩展到四面体[25],如图2所示。图2扩展的IBM的对称性间的演化图图中三个顶点对应IBM的U(5)、SU(3)、O(6)三种对称性极限,另一个顶点对应SU*(3)对称性(将SU(3)的生成元中的替换为)。由相干态理论知,U(5)、SU(3)、SU*(3)、O(6)对称性分别对应球形、轴对称长椭球形变、轴对称扁椭球形变、g-不稳定形变[4]。并且,沿球形到g-不稳定形变的相变为二级相变,临界点附近的核态具有E(5)对称性[10];从球形区到长椭球形变区的相变为一级相变,临界点附近的核态具有X(5)对称性[11];还存在球形、长椭球形和g-不稳定形变三相共存的三相点[15,16]。此外,长椭球形变与扁椭球形变之间的临界点附近的核态具有O(6)对称性[12]、Y(5)对称性[13],也有人认为长椭球与扁椭球形状相变临界点附近的核态还可能具有Z(5)对称性[14]。理论上发现形状共存和各种临界点对称性之后,很快就在实验上找到了对应的原子核。如152Sm可能有形状共存现象[26],与E(5)对称性对应的原子核有134Ba[27]、108Pd[28]、130Xe[29]等,与X(5)对称性对应的原子核有152Sm、和上面的要求相同3154Gd、156Dy和其他N=90的同中子素链[30],与Y(5)对称性对应的原子核有166,168Er[31]等,与Z(5)对称性相对应的原子核有194Pt等[25]。同时,类似Iachello四面体的工作很快被推广到区分质子玻色子和中子玻色子的IBM-2[32],同样成功的找到了各种极限对称性之间的相变。对于角动量变化可能引起的原子核形状相变,早期的研究主要基于液滴模型[19]。近年来,人们开始利用Landau相变理论[21]、有限温度推转HFB理论[22]、推转IBM[23]、推转无规位相近似[33]以及IBM框架下考虑角动量投影的相干态方法[34,35]进行研究,结果表明即使是核的低激发态也可能存在各种形状之间的相变,并说明低激发能谱中出现振动到定轴转动的相变的机制可能是,随着角动量升高,振动逐渐减弱,转动逐渐加强,临界点以后成为很好的定轴转动。另一方面,直接从核子层次对原子核形状相变的研究也已取得进展[36]。3强相互作用物质的相变3.1原子核的液气相变早在20世纪30年代,根据实验观测到的原子核的性质,人们就对原子核的结构提出了费米气体模型和液滴模型。这说明在某些条件下,原子核呈液相,或者说其某些性质表现为液相的性质;而在另一些方面,原子核表现为气相。在这一层次上,所谓的“液相”和“气相”只是作为原子核的不同性质的唯象表述,根本没有关心这两种相之间的演化。到20世纪90年代中期,随着中高能核核碰撞研究的深入,人们研究了核核碰撞形成的系统的温度与其中核子的激发能之间的关系,最早的由德国GSI报告的结果[37]如图3所示,这一关系显然与通常物质处于液相、气相及其间相变中温度与单粒子平均能量间的关系相同,从而说明发生了液气相变。由于相变通常由热力学函数和状态方程出发进行研究,原子核的液气相变自然成为研究核物质状态方程、进而研究核天体状态及其演化的突破口。于是,美国Brookhaven国家实验室、Lawrence国家实验室、Michigan州立大学、德州农机学院、俄罗斯的Dubna、德国的GSI、法国的GANIL和LNSSaclay、意大利的delSud国家实验室等国际大型实验室的核物理学家系统研究了中高能核核形成的系统的温度与单粒子激发能的关系、热容、高碎裂多重度、集体膨胀、有限尺寸及Fisher定律标度等[38~41],理论上发展了核玻尔兹曼方程[42]、有限系统费米子-分子动力学[43]、全反对称分子动力学[44]等方法、并利用渗渝理论[45]对这些系统进行研究,结果都表明,在一定的条件下,中高能核核碰撞形成的系统中都会出现液气相变,并说明该相变的机制是失稳分解。事实上,这些研究还都有待深化,尤其是相变的序参量、同位旋依赖性、相变的临界温度、对核天体的结构和演化的影响等都是目前研究关注的重要问题。图3核核碰撞形成的系统的温度与单核子能量的关系3.2强相互作用物质的相变强相互作用物质是由强子(包括重子和介子)组成的强子物质和由夸克、胶子组成的夸克物质的统称。因此,对强相互作用物质的组分、性质、相结构及相变的研究是当代原子核物理、粒子物理、天体物理和宇宙学等领域共同关注的重大课题。和上面的要求相同4我们已经知道,强子由夸克和胶子组成,并且可以形象地将之比喻为束缚有夸克和胶子的口袋,口袋内的夸克、胶子的相互作用与强相互作用真空内的作用之间的差异提供的袋常数常被用来描述束缚的强度。随着强子物质系统温度的升高,强子无规则运动的能量和其内部夸克、胶子无规则运动的能量都会升高,压强会增大;系统密度的增大也会引起压强增大,当系统的真空压不能平衡强子内部的压强时,强子将消失,夸克和胶子将成为夸克物质,也就是可以发生退禁闭相变。退禁闭形成的夸克物质可能以等离子体状态存在,从而形成夸克胶子等离子体(QGP)。另一方面,描述强相互作用的基本理论是量子色动力学(QCD),QCD具有渐近自由的性质(上述退禁闭相变正是渐近自由的结果和表现),并且零质量的费米子(夸克等)具有左旋和右旋的等价性,这种等价性称为手征对称性。然而,现实的强子世界处于低能区域,夸克是禁闭的、有质量的,并且不具有手征对称性。但当退禁闭相变发生以后,手征对称性可能恢复,从而发生手征恢复相变。再者,我们知道,由于电声作用的相互影响,声子可以为电子之间提供一个较弱的吸引力,从而形成电子库珀对,出现超导现象;由于夸克之间的特殊的相互作用道本来就是吸引的,因此夸克之间也可以形成夸克库珀对,由于夸克具有3种颜色,3种色混合或一种色与其反色混合形成无色的强子,但两个夸克形成的对却带有颜色,因此由夸克库珀对形成的凝聚状态称为色超导态[46]。根据色超导态的夸克库珀对的色味结构,色超导态具有两味色超导、色味锁定色超导等多种相(有时简单地统称之为色超导相)。目前的研究表明,强相互作用物质的相图如图4所示。图4强相互作用物质相图由于QCD具有渐近自由的性质,因此,对于高能区的场和粒子性质,可以利用微扰QCD进行研究,并得到了很好的结果。但对于低能区域,QCD的求解问题尚没有解决,于是人们发展了QCD因子化和重求和(硬热圈展开和硬密圈展开)方法[47],并利用QCD的非微扰有效场论模型方法和唯象模型方法(Dyson-Schwinger方程、瞬子模型、整体色对称模型、手征模型、孤立子模型、夸克介子耦合模型、NJL模型、袋模型)[48~54]等对强相互作用物质进行理论研究。近年来,随着对基本原理的扩展和计算方法的发展,