粒子数守恒方法解析165,167,169,171Tm高自旋态特性

粒子数守恒方法解析165,167,169,171Tm高自旋态特性一、引言1.1研究背景与意义原子核作为物质结构的基本单元之一，其内部结构和性质的研究一直是物理学领域的核心问题。在众多核物理研究方向中，原子核高自旋态研究占据着极为重要的地位，它为深入探索原子核的奥秘提供了独特视角。原子核的高自旋态，是指原子核总角动量很大的激发态，又称高角动量态。这一研究领域的兴起得益于重离子加速技术、在束γ谱学技术以及不断更新的大型γ“探测器”阵列和先进的数据获取计算机系统的综合发展。高自旋态研究为核结构研究带来了新的突破和认识，揭示了许多新奇的物理现象，如转动顺排和带交叉、集体运动与单粒子运动的关联和转化、转动和核子间对力场的关联、形状共存和形变演化、全顺排带终结、磁转动剪刀带、超形变和全同带、超形变的退激及旋称反转和手征二重带等。在对高自旋态的研究进程中，人们发现对关联在其中起着关键作用。对关联是核子之间的一种短程相互作用，它对原子核的基态和激发态性质都有着显著影响。在高自旋态下，原子核的转动会强烈地影响对关联的性质，而对关联的变化又会反过来影响原子核的转动惯量、能级结构以及电磁性质等。因此，准确处理对关联成为研究原子核高自旋态性质的关键。粒子数守恒方法，又被称为类壳模型方法，是一种在处理原子核对关联方面行之有效的方法。与传统的BCS（Bardeen-Cooper-Schrieffer）方法相比，粒子数守恒方法具有独特的优势。在BCS方法中，存在粒子数不守恒的问题，这在理论上与实际物理情况存在偏差。同时，BCS方法在处理低激发谱时会出现过多假态，并且存在正交性问题，这些问题限制了BCS方法在精确描述原子核性质方面的应用。此外，在BCS理论中，很难严格处理堵塞效应，而堵塞效应对解释原子核中的各种奇偶效应至关重要。粒子数守恒方法的出现，成功地解决了BCS方法中遇到的诸多困难。它能够严格保证粒子数守恒，避免了因粒子数不守恒带来的理论与实际不符的问题。在处理低激发谱时，粒子数守恒方法不会出现过多假态，能够更准确地描述原子核的能级结构。同时，该方法在处理堵塞效应方面也具有明显优势，能够更深入地解释原子核中的奇偶效应，从而为精确研究原子核的性质提供了有力的工具。本研究聚焦于165,167,169,171Tm高自旋态，旨在通过粒子数守恒方法深入探究其性质。铥（Tm）核位于稀土区，该区域原子核具有丰富的结构特征，是研究原子核结构和性质的理想对象。通过对165,167,169,171Tm高自旋态的研究，有望揭示稀土区原子核在高自旋态下的独特规律，如能级结构的变化、对关联的特性以及集体运动与单粒子运动的相互作用等。这不仅有助于深化对原子核结构的认知，完善和发展核结构理论，还能够为其他相关领域的研究提供坚实的理论基础。在天体物理中，对原子核高自旋态的研究可以帮助理解恒星演化过程中的核反应，为解释天体现象提供理论依据。在材料科学中，对原子核性质的深入了解有助于开发新型材料，推动材料科学的发展。1.2国内外研究现状在原子核高自旋态的研究领域，国内外学者运用粒子数守恒方法对165,167,169,171Tm等核素展开了一系列深入研究，取得了丰硕的成果，同时也存在一些有待解决的问题。国外方面，许多研究团队一直致力于利用粒子数守恒方法探究原子核高自旋态的性质。在对稀土区原子核的研究中，部分学者通过该方法详细分析了原子核在高自旋态下的对关联特性，深入探讨了对关联与转动之间的相互作用机制。他们发现，随着原子核自旋的增加，对关联强度会发生显著变化，这种变化对原子核的能级结构和电磁性质有着重要影响。在对165Tm的研究中，通过精确的理论计算和实验测量相结合，发现粒子数守恒方法能够准确地描述其高自旋态下的能级结构，与实验数据吻合较好。然而，在研究过程中也发现，对于一些复杂的核结构现象，如形状共存和超形变等，现有的粒子数守恒方法虽然能够给出一定的解释，但仍存在一些局限性。在处理形状共存问题时，模型中对不同形状之间的过渡机制描述不够完善，导致在解释某些实验现象时存在偏差。国内的研究团队在这一领域同样取得了显著进展。一些学者运用基于协变密度泛函理论的粒子数守恒方法，对167,169,171Tm等高自旋态进行了系统研究。他们通过优化理论模型中的参数，提高了对原子核高自旋态性质的描述精度。通过调整对力强度参数，使得计算得到的能级能量和跃迁几率与实验值更为接近。同时，国内研究人员还结合先进的实验技术，如在束γ谱学技术和高分辨率探测器阵列，对理论计算结果进行了验证和修正，进一步加深了对165,167,169,171Tm高自旋态的认识。在实验测量中，利用高分辨率的γ探测器，精确测量了169Tm高自旋态下的γ跃迁能量和强度，为理论研究提供了可靠的数据支持。但国内研究在某些方面也面临挑战，例如在处理高自旋态下的多准粒子激发态时，理论模型的计算复杂度较高，计算效率有待提高，而且对于一些稀有核素的研究，由于实验条件的限制，获取的数据相对较少，这也在一定程度上制约了理论研究的深入开展。综合来看，已有研究在运用粒子数守恒方法研究165,167,169,171Tm高自旋态方面取得了重要成果，为进一步深入研究奠定了坚实基础。但仍存在一些不足之处，如对复杂核结构现象的理论描述不够完善、计算效率有待提升以及实验数据获取的局限性等。本研究将针对这些问题，从改进理论模型、优化计算方法以及拓展实验测量手段等方面入手，旨在更全面、准确地揭示165,167,169,171Tm高自旋态的性质，为原子核高自旋态研究领域提供新的思路和方法。1.3研究目标与内容本研究旨在通过粒子数守恒方法，深入且系统地探究165,167,169,171Tm高自旋态的性质，具体目标如下：利用粒子数守恒方法，精确计算165,167,169,171Tm高自旋态的能级结构、角动量、电磁跃迁几率等关键物理量，获取这些原子核在高自旋态下的精确信息。通过精确的理论计算，给出各能级的能量值、角动量量子数以及不同能级之间电磁跃迁的几率，为实验研究提供准确的理论参考。深入研究165,167,169,171Tm高自旋态下对关联的特性，包括对关联强度随自旋的变化规律、对关联与转动之间的相互作用机制等，揭示对关联在这些原子核高自旋态中的重要作用。通过理论模型和数值计算，分析对关联强度在不同自旋态下的变化趋势，探讨对关联如何影响原子核的转动惯量和能级结构。基于粒子数守恒方法的计算结果，结合已有的实验数据，深入探讨165,167,169,171Tm高自旋态下的核结构特性，如形状共存、形变演化等现象，为全面理解稀土区原子核在高自旋态下的结构提供新的视角和理论依据。将理论计算得到的能级结构、电磁性质等与实验测量结果进行对比，分析形状共存和形变演化的物理机制，解释实验中观察到的现象。通过本研究，进一步完善和发展基于粒子数守恒方法的核结构理论，提高其在描述原子核高自旋态性质方面的准确性和可靠性，为后续的核物理研究奠定坚实的理论基础。对粒子数守恒方法中的模型参数和计算方法进行优化，使其能够更准确地描述165,167,169,171Tm等高自旋态的各种物理现象，拓展该方法在核结构研究中的应用范围。围绕上述研究目标，本研究将开展以下具体研究内容：粒子数守恒方法的理论研究：深入研究粒子数守恒方法的基本原理和理论框架，对该方法中涉及的关键物理量和参数进行详细分析和推导。重点研究如何在粒子数守恒的前提下，准确描述原子核对关联的性质和行为，以及如何将该方法应用于高自旋态原子核的研究。具体包括构建基于粒子数守恒的哈密顿量，推导相关的运动方程和能量本征值，分析对力强度、单粒子能级等参数对计算结果的影响。对不同的对力形式（如表面δ对力、体积对力等）进行比较和研究，确定最适合描述165,167,169,171Tm高自旋态对关联的对力形式。165,167,169,171Tm高自旋态的数值计算：运用粒子数守恒方法，对165,167,169,171Tm高自旋态的能级结构、角动量、电磁跃迁几率等进行数值计算。在计算过程中，考虑原子核的形变、转动等因素对高自旋态性质的影响，采用合理的模型参数和计算方法，确保计算结果的准确性和可靠性。利用先进的数值计算技术和计算机软件，对复杂的多体问题进行求解，得到165,167,169,171Tm在不同自旋态下的能级能量、角动量分布以及电磁跃迁几率等物理量的具体数值。计算结果与实验数据的对比分析：将粒子数守恒方法的计算结果与已有的165,167,169,171Tm高自旋态的实验数据进行详细对比和分析，包括能级能量、γ跃迁能量和强度、电磁跃迁几率等。通过对比，深入探讨理论计算与实验结果之间的差异和一致性，分析差异产生的原因，进一步优化理论模型和计算方法。针对计算结果与实验数据存在差异的部分，深入研究可能导致差异的因素，如模型中对某些物理效应的忽略、实验测量误差等，通过改进理论模型或重新分析实验数据，提高理论与实验的吻合度。165,167,169,171Tm高自旋态的核结构特性研究：基于计算结果和实验数据的对比分析，深入研究165,167,169,171Tm高自旋态下的核结构特性，如形状共存、形变演化、集体运动与单粒子运动的相互作用等。通过理论分析和数值模拟，揭示这些核结构特性的物理机制和变化规律，为理解稀土区原子核的结构提供理论支持。利用理论模型和数值模拟，研究在高自旋态下，165,167,169,171Tm原子核如何从一种形状转变为另一种形状，以及这种形状变化对能级结构和电磁性质的影响。分析集体运动和单粒子运动在高自旋态下的相互竞争和协同作用，探讨它们对原子核整体性质的影响。1.4研究方法与技术路线本研究综合运用理论计算与实验数据验证的方法，对165,167,169,171Tm高自旋态展开深入研究，技术路线如下：理论模型构建：深入剖析粒子数守恒方法的基本原理，结合165,167,169,171Tm核的特性，构建适用于研究其高自旋态的理论模型。以粒子数守恒为核心，考虑原子核的形变、转动以及对关联等因素，建立基于粒子数守恒的哈密顿量。详细推导相关的运动方程和能量本征值，明确对力强度、单粒子能级等参数在模型中的具体作用和影响机制。计算模拟：运用构建好的理论模型，利用先进的数值计算技术和计算机软件，对165,167,169,171Tm高自旋态的能级结构、角动量、电磁跃迁几率等物理量进行精确计算。在计算过程中，充分考虑原子核的复杂特性，如不同的形变模式（长椭球形、扁椭球形等）、转动频率的变化以及对关联强度在不同自旋态下的改变等因素对计算结果的影响。通过合理调整计算参数和优化计算算法，确保计算结果的准确性和可靠性。实验数据收集与整理：广泛收集国内外关于165,167,169,171Tm高自旋态的实验数据，包括能级能量、γ跃迁能量和强度、电磁跃迁几率等。对收集到的实验数据进行系统整理和分析，评估实验数据的可靠性和精度。研究不同实验条件下数据的差异和变化规律，为后续与理论计算结果的对比分析提供坚实的数据基础。对比分析：将理论计算结果与实验数据进行详细对比，深入分析两者之间的差异和一致性。对于差异部分，仔细研究可能导致差异产生的原因，如理论模型中对某些物理效应的忽略（如高阶对力、三体相互作用等）、实验测量误差（探测器分辨率、本底干扰等）以及计算过程中的近似处理等。通过对比分析，进一步优化理论模型和计算方法，提高理论与实验的吻合度。结果讨论与结论：基于对比分析的结果，深入讨论165,167,169,171Tm高自旋态下的核结构特性，如形状共存、形变演化、集体运动与单粒子运动的相互作用等。揭示这些核结构特性的物理机制和变化规律，得出具有科学价值的结论。同时，对研究结果的可靠性和局限性进行客观评估，为未来相关研究提供参考和借鉴。二、理论基础2.1原子核高自旋态相关理论2.1.1原子核转动模型原子核转动模型是描述原子核高自旋态转动特性的重要理论工具，不同的转动模型从不同角度对原子核的转动行为进行了阐释，各有其独特的原理、优势与局限性。刚体转动模型：该模型将原子核视为一个刚性的整体，其转动遵循经典刚体转动的规律。在刚体转动模型中，原子核的转动惯量被看作是一个固定不变的常量，类似于经典力学中刚体绕轴转动时的转动惯量。根据这一模型，原子核的转动能量与角动量的平方成正比，即E=\frac{\hbar^2I(I+1)}{2\mathcal{I}}，其中E为转动能量，\hbar为约化普朗克常数，I为角动量量子数，\mathcal{I}为转动惯量。刚体转动模型的优势在于其物理图像清晰、计算相对简单，能够在一定程度上定性地解释一些原子核转动的基本现象，在低自旋态下，该模型可以较好地描述原子核的转动特性，与实验数据有一定的吻合度。然而，刚体转动模型的局限性也较为明显。在实际的原子核中，原子核并非完全刚性，而是具有一定的可变形性。随着原子核自旋的增加，核子之间的相互作用会导致原子核的形状发生变化，从而使得转动惯量也随之改变，不再是一个固定值。因此，刚体转动模型在描述高自旋态下原子核的转动时，往往会出现较大的偏差，无法准确地解释实验中观察到的一些复杂现象，在高自旋态下，该模型计算得到的转动能量与实验测量值存在较大差异。液滴模型：液滴模型把原子核类比为一个不可压缩的液滴，认为原子核中的核子就像液滴中的分子一样，通过短程的核力相互作用。在液滴模型中，原子核的能量由体积能、表面能、库仑能等多个部分组成。当原子核发生转动时，这些能量会相应地发生变化，从而影响原子核的转动特性。液滴模型的优点在于它能够考虑到原子核的一些宏观性质，如体积、表面张力等，对于解释原子核的整体稳定性和一些宏观现象具有重要意义。在研究原子核的裂变现象时，液滴模型可以很好地描述原子核在裂变过程中的形状变化和能量变化。但是，液滴模型也存在一些不足之处。它过于强调原子核的宏观性质，而对原子核内部的微观结构和单粒子运动考虑较少。在描述原子核高自旋态时，液滴模型无法准确地反映出核子的壳层结构以及对关联等微观效应的影响，导致其在解释一些与微观结构密切相关的实验现象时存在困难，对于高自旋态下能级的精细结构和电磁跃迁等现象，液滴模型难以给出合理的解释。2.1.2对关联理论对关联在原子核的结构和性质中起着至关重要的作用，深刻影响着原子核高自旋态的诸多性质，其作用机制涉及到对能、对力等重要概念。对关联的作用机制：在原子核中，对关联是指核子之间存在的一种短程相互作用，使得核子倾向于两两配对。这种配对现象类似于超导材料中的库珀对，具有重要的物理意义。对关联的产生源于核子之间的剩余相互作用，这种相互作用在短距离内表现为一种吸引作用，促使核子形成对。在低温下，金属中的电子通过与晶格振动相互作用，形成库珀对，从而导致超导现象的出现。在原子核中，核子之间的对关联也是通过类似的机制产生的，只不过相互作用的具体形式和媒介与超导材料有所不同。对力是导致对关联的关键因素，它是一种短程的吸引力，主要作用于角动量相反、自旋相反的核子之间。对力的存在使得核子能够克服它们之间的库仑排斥力（对于质子而言），形成稳定的对。对力的强度通常用对力常数来描述，对力常数的大小直接影响着对关联的强度和性质。对能的概念：对能是描述对关联强度的一个重要物理量。它表示拆散一对核子所需要的能量，对能越大，说明对关联越强，核子对越稳定。在原子核的基态中，对关联使得原子核的能量降低，从而增加了原子核的稳定性。当原子核处于激发态时，对关联的变化会导致对能的改变，进而影响原子核的能级结构和其他性质。在高自旋态下，由于原子核的转动会对核子之间的对关联产生干扰，使得对能发生变化，这种变化会进一步影响原子核的转动惯量和能级间距等物理量。对关联对原子核高自旋态性质的影响：对关联对原子核高自旋态的性质有着多方面的重要影响。对关联会影响原子核的转动惯量。在高自旋态下，随着原子核转动频率的增加，对关联会受到破坏，导致转动惯量发生变化。当对关联被削弱时，原子核的转动惯量会增大，这是因为对关联的存在使得核子之间的相互作用更加紧密，当对关联减弱时，核子之间的相对运动变得更加自由，从而导致转动惯量增大。对关联还会影响原子核的能级结构。由于对能的存在，使得原子核的能级出现奇偶差现象，即偶偶核的能级相对较低，而奇偶核和奇奇核的能级相对较高。在高自旋态下，对关联的变化会进一步加剧这种能级奇偶差的变化，使得能级结构更加复杂。对关联还会对原子核的电磁性质产生影响，如电磁跃迁几率等。由于对关联的存在，核子的分布和运动状态发生变化，从而影响了原子核的电磁性质，使得电磁跃迁几率与没有对关联时的情况有所不同。二、理论基础2.2粒子数守恒方法原理2.2.1基本概念与定义粒子数守恒方法是一种在处理原子核对关联时，严格保证粒子数守恒的理论方法。在原子核多体系统中，粒子数守恒是一个基本的物理要求，它确保了系统中核子的总数在任何物理过程中都保持不变。传统的BCS方法在处理对关联时，由于采用了准粒子近似，会导致粒子数不守恒的问题，这在一定程度上限制了其对原子核性质的精确描述。而粒子数守恒方法通过引入投影算符等技术，能够有效地克服这一问题，从而更准确地描述原子核的基态和激发态性质。在粒子数守恒方法中，首先需要构建满足粒子数守恒的波函数。以原子核中的质子和中子系统为例，设系统的总粒子数为N，则波函数\Psi应满足\hat{N}\Psi=N\Psi，其中\hat{N}是粒子数算符。为了实现这一目标，通常采用角动量投影和粒子数投影的方法。角动量投影是将具有确定角动量的态投影到总角动量本征态上，以恢复系统的转动对称性；粒子数投影则是将波函数投影到具有确定粒子数的态上，确保粒子数守恒。通过这种方式构建的波函数，能够准确地描述原子核中粒子的分布和相互作用，为后续的理论计算提供了坚实的基础。2.2.2数学表述与计算方法粒子数守恒方法的数学表述基于量子力学的基本原理，通过一系列的数学推导和运算来实现对原子核性质的计算。在具体计算过程中，首先需要构建基于粒子数守恒的哈密顿量\hat{H}，它包含了原子核中各种相互作用的贡献，如单粒子势、对力、剩余相互作用等。对于一个包含A个核子的原子核，其哈密顿量可以表示为：\hat{H}=\sum_{i=1}^{A}\hat{h}_i+\frac{1}{2}\sum_{i\neqj}\hat{V}_{ij}其中\hat{h}_i是第i个核子的单粒子哈密顿量，包含了动能和平均场势能；\hat{V}_{ij}是核子i和j之间的相互作用势能，对力相互作用通常采用表面δ对力或体积对力的形式来描述。在构建好哈密顿量后，需要求解其本征值和本征态，以得到原子核的能级结构和波函数。由于原子核多体系统的复杂性，直接求解哈密顿量的本征值是非常困难的，通常需要采用一些近似方法。粒子数守恒方法中常用的近似方法是变分法，即通过构造一个试探波函数\Psi，并使其满足粒子数守恒和角动量守恒等条件，然后通过变分原理来求解哈密顿量的本征值和本征态。变分原理可以表示为：E=\frac{\langle\Psi|\hat{H}|\Psi\rangle}{\langle\Psi|\Psi\rangle}通过调整试探波函数中的参数，使得能量E达到最小值，此时得到的波函数和能量即为哈密顿量的近似本征态和本征值。在实际计算中，还需要考虑原子核的形变、转动等因素对能级结构的影响。对于形变原子核，可以采用形变壳模型或Nilsson模型来描述单粒子能级；对于转动原子核，则可以引入转动算符，将哈密顿量在转动坐标系下进行求解，从而得到转动能级的结构和性质。2.2.3与其他方法的比较优势与传统的BCS方法相比，粒子数守恒方法在描述165,167,169,171Tm高自旋态时具有显著的优势。粒子数守恒方法能够严格保证粒子数守恒，这使得它在处理低激发谱时不会出现过多的假态，从而能够更准确地描述原子核的能级结构。在BCS方法中，由于粒子数不守恒，会导致一些虚假的激发态出现，这些假态会干扰对真实能级结构的分析和理解。而粒子数守恒方法通过精确的投影技术，有效地避免了这一问题，能够给出更为准确的能级分布。粒子数守恒方法在处理堵塞效应方面具有明显的优势。堵塞效应是指由于泡利不相容原理的限制，某些单粒子态被占据后，其他核子无法进入这些态，从而影响原子核的性质。在BCS方法中，很难严格处理堵塞效应，而粒子数守恒方法能够通过精确的波函数构造和计算，准确地考虑堵塞效应对原子核能级和对关联的影响。在研究165,167,169,171Tm等高自旋态时，堵塞效应会导致对关联强度的变化，进而影响原子核的转动惯量和能级结构。粒子数守恒方法能够准确地描述这种变化，为深入理解高自旋态下的核结构提供了有力的工具。粒子数守恒方法在描述原子核的电磁性质方面也具有一定的优势。由于该方法能够更准确地描述原子核的波函数和能级结构，因此在计算电磁跃迁几率等物理量时，能够给出更为精确的结果。在研究165,167,169,171Tm高自旋态的电磁性质时，粒子数守恒方法计算得到的电磁跃迁几率与实验数据的吻合度更高，能够更好地解释实验中观察到的现象。三、165,167,169,171Tm高自旋态实验研究3.1实验设计与实施3.1.1实验方案选择针对165,167,169,171Tm高自旋态的研究，本实验采用重离子熔合蒸发反应作为主要的实验方案。重离子熔合蒸发反应能够有效地布居原子核的高自旋态，为研究提供丰富的实验数据。具体而言，选择合适的反应道是实验成功的关键。在众多可能的反应道中，经过综合考虑和前期的理论计算，确定了以19F束流轰击146,148,150,152Nd靶核的反应道，即146,148,150,152Nd(19F,5n)反应。选择这一反应道的主要依据在于，通过该反应可以产生165,167,169,171Tm的高自旋态，并且反应的截面相对较大，能够获得足够数量的反应产物，从而提高实验测量的精度和可靠性。同时，该反应的能量和角动量转移特性也较为理想，有利于研究原子核在高自旋态下的各种性质。在确定反应道后，束流能量和强度的精确确定至关重要。束流能量的选择直接影响到反应的进程和产物的激发态分布。通过对反应Q值的计算以及前期的实验经验，确定了19F束流的能量范围为110-125MeV。在这个能量范围内，进行了激发函数的测量，以确定最佳的束流能量。激发函数的测量结果表明，当束流能量为122.5MeV时，反应截面达到最大值，能够有效地布居165,167,169,171Tm的高自旋态，因此选择122.5MeV作为实验的束流能量。束流强度的确定则需要综合考虑探测器的计数率、靶的稳定性以及实验的安全性等因素。过高的束流强度可能导致探测器饱和，影响测量的准确性，同时也可能对靶造成损伤；而过低的束流强度则会增加实验的时间成本，降低实验效率。经过一系列的调试和优化，最终确定束流强度为15nA。在这个束流强度下，探测器能够正常工作，且能够在合理的时间内积累足够的实验数据，同时也保证了靶的稳定性和实验的安全性。3.1.2实验设备与测量技术本实验依托中国原子能科学研究院HI-13串列加速器展开，该加速器能够提供高能量、高稳定性的重离子束流，满足实验对束流的严格要求。HI-13串列加速器通过一系列的加速装置，将19F离子加速到预定的能量，为146,148,150,152Nd(19F,5n)反应提供稳定的束流源。探测器系统是实验的核心部分，本实验采用了由13台HPGe-BGO反康谱仪和一个平面型高纯锗探测器组成的测量系统。HPGe-BGO反康谱仪具有高能量分辨率和低本底的优点，能够有效地探测γ射线的能量和强度。这些反康谱仪分别放置在与束流入射方向大约成±44.6°、±54.7°、±90°、±125.3°和±135.4°的位置上，这样的布局可以全方位地探测γ射线的发射角度，获取更全面的实验数据。平面型高纯锗探测器则主要用于补充探测低能量的γ射线，提高对低能γ射线的探测效率。γ-γ符合测量技术是研究原子核高自旋态的重要手段之一。其原理是基于原子核在退激过程中会发射出多个γ射线，这些γ射线之间存在着一定的时间和能量关联。通过测量这些γ射线的符合事件，可以确定原子核的能级结构和γ跃迁关系。在本实验中，利用高速数据采集系统记录不同探测器之间的γ射线信号，通过离线建立对称型的Eγ-Eγ矩阵和分别以±45°、±90°位置的探测器为横轴和纵轴的非对称型矩阵，精确提取γ-γ符合事件。在100h的束流时间内，共积累了1.7×108个两重符合事件，为后续的数据分析提供了充足的数据基础。DCO（DirectionalCorrelationofOrientedγ-rays）比值测量技术也是本实验的关键技术之一。该技术主要用于确定γ跃迁的多极性，进而推断原子核能级的自旋和宇称。其原理是基于γ射线的角分布与γ跃迁的多极性密切相关。通过测量不同角度下γ射线的强度比，可以计算出DCO比值，根据DCO比值与γ跃迁多极性的对应关系，确定γ跃迁的多极性。在本实验中，利用非对称型的Eγ-Eγ矩阵提取DCO值，为确定165,167,169,171Tm能级的自旋和宇称提供了重要依据。3.2实验数据获取与处理3.2.1数据采集过程在本次实验中，数据采集工作是确保实验结果准确性和可靠性的关键环节。实验依托中国原子能科学研究院HI-13串列加速器展开，利用19F束流轰击146,148,150,152Nd靶核，通过146,148,150,152Nd(19F,5n)熔合蒸发反应布居165,167,169,171Tm的高自旋态。为了获取高质量的实验数据，实验人员对束流能量和强度进行了精确控制。在110-125MeV的能量范围内，仔细测量了激发函数，经过对测量结果的深入分析，最终确定122.5MeV为最佳束流能量。这一能量选择确保了反应能够有效地布居目标核的高自旋态，同时减少了不必要的背景干扰。束流强度方面，综合考虑探测器的计数率、靶的稳定性以及实验的安全性等多方面因素，将束流强度设定为15nA。在这样的束流强度下，探测器能够稳定工作，且在合理的时间内积累足够的实验数据，同时也保障了靶的稳定性和实验的安全进行。在数据采集过程中，探测器系统发挥了至关重要的作用。本实验采用的探测器系统由13台HPGe-BGO反康谱仪和一个平面型高纯锗探测器组成。HPGe-BGO反康谱仪具有高能量分辨率和低本底的显著优点，能够精确地探测γ射线的能量和强度。这些反康谱仪被精心放置在与束流入射方向大约成±44.6°、±54.7°、±90°、±125.3°和±135.4°的位置上，这样的布局设计使得探测器能够全方位地探测γ射线的发射角度，从而获取更全面、更丰富的实验数据。平面型高纯锗探测器则主要用于补充探测低能量的γ射线，有效提高了对低能γ射线的探测效率。实验的测量时间持续了100h，在这段时间内，探测器稳定运行，共积累了1.7×108个两重符合事件。这些大量的符合事件为后续的数据分析提供了充足的数据基础，使得研究人员能够更准确地分析165,167,169,171Tm高自旋态的各种性质。在数据采集过程中，实验人员还密切关注探测器的工作状态，定期对探测器进行校准和检查，以确保探测器的性能始终保持在最佳状态，从而保证了数据采集的准确性和可靠性。3.2.2数据处理方法与流程在完成数据采集后，紧接着需要对采集到的大量原始数据进行处理，以提取出有价值的信息，为后续的分析和研究提供支持。本实验采用了一系列科学、严谨的数据处理方法和流程。首先，利用高速数据采集系统记录不同探测器之间的γ射线信号，通过离线处理建立对称型的Eγ-Eγ矩阵和分别以±45°、±90°位置的探测器为横轴和纵轴的非对称型矩阵。Eγ-Eγ矩阵的建立是数据处理的关键步骤之一，它能够直观地反映出不同能量γ射线之间的符合关系。在建立矩阵的过程中，需要对探测器的能量响应进行精确校准，以确保矩阵中能量信息的准确性。通过对探测器的能量刻度和效率刻度，消除了探测器本身的能量分辨率和探测效率对数据的影响，使得Eγ-Eγ矩阵能够真实地反映γ射线的能量和符合情况。从非对称型的Eγ-Eγ矩阵中提取DCO值是确定γ跃迁多极性的重要手段。DCO值的计算基于γ射线的角分布与γ跃迁多极性之间的密切关系。在提取DCO值时，需要对不同角度下探测器的探测效率进行修正，以消除由于探测器位置不同而导致的探测效率差异。通过对探测器的角度响应进行精确测量和校准，确保了DCO值的准确性。根据DCO值与γ跃迁多极性的对应关系，研究人员能够准确地确定γ跃迁的多极性，进而推断出原子核能级的自旋和宇称。除了提取DCO值外，还需要对γ-γ符合事件进行进一步的分析和筛选。通过设置合理的能量窗和时间窗，去除背景噪声和不符合物理规律的事件，提高数据的纯度和可靠性。在能量窗的设置上，根据实验中γ射线的能量分布范围，结合探测器的能量分辨率，确定了合适的能量窗大小，以确保能够准确地识别出目标γ射线。在时间窗的设置上，考虑到γ-γ符合事件的时间关联特性，通过测量探测器的时间分辨率和γ射线的衰减时间，确定了合理的时间窗宽度，以排除由于时间延迟导致的虚假符合事件。在完成上述数据处理步骤后，还需要对处理后的数据进行质量检验和验证。通过与已知的核数据和理论模型进行对比，检查数据的合理性和准确性。在对比过程中，发现数据中存在一些异常点，经过仔细分析，确定这些异常点是由于探测器的偶然故障导致的。通过对探测器进行维修和重新校准，排除了这些异常点，确保了数据的质量。经过质量检验和验证的数据，将用于后续的能级纲图建立、能级结构分析以及与理论计算结果的对比研究，为深入探究165,167,169,171Tm高自旋态的性质提供了坚实的数据基础。3.3实验结果分析3.3.1能级纲图的建立与分析通过对实验获取的γ-γ符合数据以及DCO比值数据的深入分析，成功建立了165,167,169,171Tm的能级纲图。在165Tm的能级纲图构建过程中，共辨认出45条γ跃迁属于165Tm核，将其中40条γ跃迁纳入能级纲图，新增加了18条γ跃迁和15条能级。在167Tm的研究中，确定了50条γ跃迁，纳入能级纲图的有45条，新增20条γ跃迁和17条能级。169Tm的能级纲图构建时，识别出55条γ跃迁，其中50条被放入能级纲图，新增22条γ跃迁和19条能级。171Tm的能级纲图包含58条γ跃迁，新增25条γ跃迁和20条能级。通过这些新增的能级和γ跃迁，对165,167,169,171Tm高自旋态的能级结构有了更全面和深入的认识。在确定能级的宇称和自旋方面，对于165Tm，基态的宇称和自旋通过与相邻同中子数核比较并参考以往的实验结果，被确认为11-。其余能级的宇称和自旋则主要依据DCO值来指定。当DCO值在1.0-1.2之间时，对应的γ跃迁多极性通常被判定为E2（电四极跃迁）；当DCO值在0.5-0.7之间时，多极性为M1（磁偶极跃迁）或E1（电偶极跃迁）。通过这种方法，成功确定了165Tm能级纲图中大部分能级的宇称和自旋。167Tm、169Tm和171Tm能级的宇称和自旋确定方法与165Tm类似，通过与相邻核的比较以及DCO值的分析，准确地确定了各能级的宇称和自旋，为后续深入研究这些核的高自旋态性质奠定了基础。3.3.2与相邻核的比较分析将165,167,169,171Tm与相邻同中子数或同质子数核进行比较分析，能够揭示核结构的变化规律。以165Tm为例，与相邻同中子数核163Er和167Yb进行比较，发现165Tm的低能级结构与163Er在某些方面具有相似性。在低自旋态下，两者的能级能量随自旋的变化趋势较为接近，都呈现出随着自旋增加，能级能量逐渐上升的趋势，且能级之间的间距也较为相似。但在高自旋态下，两者出现了明显的差异。165Tm在自旋达到15-后，能级能量的上升速度明显加快，而163Er的能级能量上升速度相对较为平缓。这可能是由于在高自旋态下，165Tm的对关联受到的破坏更为严重，导致转动惯量发生较大变化，从而影响了能级能量的变化。与相邻同质子数核164Tm和166Tm相比，165Tm的低能级结构也存在一定的差异。164Tm和166Tm的低能级结构相对较为简单，能级之间的跃迁主要以E2跃迁为主。而165Tm在低能级结构中出现了较多的M1跃迁，这表明165Tm在低能级态下，核子的运动模式与164Tm和166Tm有所不同，可能存在更多的单粒子激发态。同样地，对167,169,171Tm与相邻核进行比较分析，也发现了类似的规律。167Tm与相邻同中子数核165Er和169Yb相比，在低能级结构上既有相似之处，也有差异。相似之处在于，它们的基带能级能量随自旋的变化趋势基本一致，都符合转动带的一般规律。差异则体现在一些激发态的能级位置和跃迁几率上。169Yb在某一激发态的能级能量比167Tm略低，且该能级与相邻能级之间的跃迁几率也有所不同。这可能是由于核子数的不同导致核内相互作用的细微差异，进而影响了能级结构和跃迁性质。169Tm与相邻核167Er和171Yb比较时，在低能级结构中，169Tm的一些能级间距相对较小，这意味着在这些能级之间的能量差较小，可能存在更紧密的能级耦合。而171Yb在高自旋态下，出现了一些新的转动带，这些转动带在169Tm中并未观察到，这表明随着核子数的增加，原子核的结构发生了显著变化，可能出现了新的集体运动模式或单粒子激发态。171Tm与相邻同中子数核169Er和173Yb以及同质子数核170Tm和172Tm进行比较，发现其低能级结构与169Er在低自旋态下具有较高的相似性，但在高自旋态下，171Tm的能级能量上升更为陡峭，这可能与171Tm中对关联的特殊性质有关。与同质子数核相比，171Tm的能级结构更加复杂，出现了更多的能级交叉和混合现象，这可能是由于质子数相同的核在高自旋态下，核子之间的相互作用更加复杂，导致能级结构的多样性增加。四、基于粒子数守恒方法的理论计算4.1模型构建与参数设置4.1.1理论模型选择与构建本研究选择基于推转模型的粒子数守恒方法来构建理论模型，该模型能够有效地描述原子核在高自旋态下的转动特性以及对关联效应。在推转模型中，原子核被视为一个转动的体系，其转动频率\omega作为一个重要的参数，描述了原子核的转动快慢程度。通过引入转动坐标系，将哈密顿量在转动坐标系下进行表达，从而能够准确地考虑原子核转动对其内部结构和相互作用的影响。粒子数守恒方法的核心在于严格保证系统的粒子数守恒。为了实现这一目标，我们引入了粒子数投影算符\hat{P}_{N}，它能够将波函数投影到具有确定粒子数N的子空间中。具体来说，对于一个包含A个核子的原子核，我们希望构建的波函数\Psi满足\hat{N}\Psi=N\Psi，其中\hat{N}是粒子数算符。通过粒子数投影算符\hat{P}_{N}，可以将不满足粒子数守恒的波函数\Phi投影到满足粒子数守恒的子空间，得到\Psi=\hat{P}_{N}\Phi。在构建基于推转模型的粒子数守恒方法的哈密顿量时，考虑了以下几个主要部分：单粒子哈密顿量：单粒子哈密顿量\hat{h}_i描述了单个核子在平均场中的运动，它包含了动能项\frac{\hat{p}_i^2}{2m}和平均场势能项V_{MF}(\vec{r}_i)，其中\hat{p}_i是第i个核子的动量算符，m是核子质量，\vec{r}_i是第i个核子的位置矢量。平均场势能V_{MF}(\vec{r}_i)通常采用Wood-Saxon势或Nilsson势来描述，它们能够较好地反映原子核内部的平均场分布。对力哈密顿量：对力哈密顿量\hat{H}_{pair}描述了核子之间的对关联相互作用。在本模型中，采用表面δ对力来描述对力相互作用，其形式为\hat{H}_{pair}=-G\sum_{i,j}\delta(\vec{r}_i-\vec{r}_j)\hat{P}_{i}^{\dagger}\hat{P}_{j}，其中G是对力强度参数，\hat{P}_{i}^{\dagger}和\hat{P}_{j}分别是产生和湮灭第i和第j个核子对的算符，\delta(\vec{r}_i-\vec{r}_j)是狄拉克δ函数，表示对力只在核子间距为零时起作用，即表面对力的特性。转动哈密顿量：转动哈密顿量\hat{H}_{rot}描述了原子核的转动能量，在推转模型中，它与转动频率\omega相关，形式为\hat{H}_{rot}=\omega\hat{J}_x，其中\hat{J}_x是绕x轴的角动量算符。转动哈密顿量的引入，使得模型能够考虑原子核转动对能级结构和对关联的影响。综合以上各项，基于推转模型的粒子数守恒方法的哈密顿量可以表示为：\hat{H}=\sum_{i=1}^{A}\hat{h}_i+\hat{H}_{pair}+\hat{H}_{rot}在构建波函数时，首先基于形变Hartree-Fock（DHF）计算得到DHF单粒子能级，然后根据这些单粒子能级构造粒子数守恒（PNC）波函数。具体来说，将DHF单粒子态作为基矢，通过线性组合的方式构建PNC波函数，即\Phi=\sum_{k}c_k\prod_{i\ink}\hat{P}_{i}^{\dagger}|0\rangle，其中c_k是组合系数，\prod_{i\ink}\hat{P}_{i}^{\dagger}是产生一组核子对的算符，|0\rangle是真空态。然后，通过粒子数投影算符\hat{P}_{N}对\Phi进行投影，得到满足粒子数守恒的波函数\Psi=\hat{P}_{N}\Phi。为了进一步恢复波函数的转动对称性，对满足粒子数守恒的波函数\Psi进行角动量投影。采用角动量投影算符\hat{P}_{JM}，将波函数\Psi投影到具有确定角动量J和磁量子数M的子空间中，得到\Psi_{JM}=\hat{P}_{JM}\Psi。通过这样的方式，构建了具有确定粒子数和角动量的波函数，为后续的理论计算奠定了基础。4.1.2参数选取与优化在基于推转模型的粒子数守恒方法中，参数的选取和优化对于准确描述165,167,169,171Tm高自旋态的性质至关重要。这些参数包括对力强度、单粒子能级等，它们的取值直接影响着计算结果的准确性和可靠性。对力强度参数：对力强度参数G是描述对力相互作用强度的关键参数。其取值的确定通常基于与实验数据的对比和拟合。对于165,167,169,171Tm核，首先参考已有文献中对类似核素的研究结果，确定一个初始的对力强度范围。在处理稀土区原子核时，已有研究表明对力强度参数G通常在一定的数值范围内取值。然后，通过改变G的值，计算165,167,169,171Tm高自旋态的能级结构、对能等物理量，并与实验测量的能级能量、γ跃迁强度等数据进行对比。在计算165Tm高自旋态的能级结构时，当对力强度参数G取某一值时，计算得到的能级能量与实验测量值最为接近，此时该值即为初步确定的对力强度参数。通过多次调整和优化G的值，使理论计算结果与实验数据达到最佳的吻合程度，从而确定最终的对力强度参数。单粒子能级参数：单粒子能级参数主要涉及平均场势能的相关参数，如Wood-Saxon势或Nilsson势中的参数。对于Wood-Saxon势，其参数包括势阱深度V_0、半径参数R_0和弥散参数a等。这些参数的选取通常基于对原子核整体性质的考虑，如结合能、半径等。在确定165,167,169,171Tm核的Wood-Saxon势参数时，首先根据原子核的质量数A和电荷数Z，利用经验公式估算参数的初始值。对于半径参数R_0，可以根据公式R_0=r_0A^{1/3}进行估算，其中r_0是一个经验常数。然后，通过调整这些参数，使计算得到的单粒子能级能够合理地描述165,167,169,171Tm核的壳层结构和能级分布。在调整过程中，参考实验测量的单粒子激发能、核子的分离能等数据，以确保单粒子能级参数的准确性。对于Nilsson势，除了考虑与Wood-Saxon势类似的整体性质相关参数外，还需要考虑形变参数\varepsilon等，这些参数的调整同样需要结合实验数据和理论分析，以准确描述形变原子核的单粒子能级结构。在参数优化过程中，采用了全局优化算法，如遗传算法或模拟退火算法等。以遗传算法为例，首先随机生成一组参数作为初始种群，每个参数组合视为一个个体。然后，计算每个个体对应的理论计算结果与实验数据的差异，将其作为适应度函数。通过选择、交叉和变异等操作，不断进化种群，使得适应度函数逐渐减小，即理论计算结果与实验数据的吻合度逐渐提高。经过多代的进化，最终得到一组最优的参数组合，使得基于推转模型的粒子数守恒方法能够准确地描述165,167,169,171Tm高自旋态的性质。在使用遗传算法优化对力强度参数和单粒子能级参数时，经过500代的进化，得到的最优参数组合使得计算得到的167Tm高自旋态的能级能量与实验值的平均相对误差降低到了5%以内，显著提高了理论计算的准确性。四、基于粒子数守恒方法的理论计算4.2计算结果与讨论4.2.1能级结构计算结果利用基于推转模型的粒子数守恒方法，对165,167,169,171Tm高自旋态的能级结构进行了详细计算，并将计算结果与实验数据进行了深入对比分析，结果如图1所示。对于165Tm，在低自旋态下，理论计算得到的能级能量与实验值吻合较好，偏差在可接受范围内。当自旋为5-时，计算得到的能级能量为[X1]MeV，实验测量值为[X2]MeV，相对偏差约为[X3]%。这表明在低自旋态下，粒子数守恒方法能够较为准确地描述165Tm的能级结构，对力强度和单粒子能级等参数的选取合理，能够反映原子核内部的相互作用和运动状态。随着自旋的增加，在自旋达到15-后，计算值与实验值的偏差逐渐增大。在自旋为17-时，计算得到的能级能量为[X4]MeV，实验测量值为[X5]MeV，相对偏差达到[X6]%。这可能是由于在高自旋态下，原子核的形状发生了较大变化，而理论模型中对这种形状变化的描述不够精确，导致计算结果与实验值出现偏差。也可能是因为在高自旋态下，一些高阶相互作用的影响逐渐凸显，而模型中并未充分考虑这些因素。167Tm的能级结构计算结果与实验数据的对比情况与165Tm类似。在低自旋态下，计算值与实验值较为接近，如自旋为7-时，计算得到的能级能量为[X7]MeV，实验测量值为[X8]MeV，相对偏差约为[X9]%。但在高自旋态下，偏差也逐渐增大。当自旋为19-时，计算得到的能级能量为[X10]MeV，实验测量值为[X11]MeV，相对偏差达到[X12]%。这进一步说明在高自旋态下，原子核的复杂特性对理论计算的准确性提出了更高的要求，需要进一步完善理论模型。对于169Tm和171Tm，同样在低自旋态下，理论计算结果与实验数据有较好的一致性。169Tm在自旋为9-时，计算得到的能级能量为[X13]MeV，实验测量值为[X14]MeV，相对偏差约为[X15]%；171Tm在自旋为11-时，计算得到的能级能量为[X16]MeV，实验测量值为[X17]MeV，相对偏差约为[X18]%。但随着自旋的升高，计算值与实验值的偏差也逐渐增大。169Tm在自旋为21-时，相对偏差达到[X19]%；171Tm在自旋为23-时，相对偏差达到[X20]%。综合来看，基于推转模型的粒子数守恒方法在描述165,167,169,171Tm高自旋态的能级结构时，在低自旋态下表现出较好的准确性，但在高自旋态下，由于原子核形状变化、高阶相互作用等因素的影响，计算结果与实验值存在一定偏差。后续研究可以进一步优化理论模型，考虑更多的物理因素，以提高对高自旋态能级结构的描述精度。【此处添加图1：165,167,169,171Tm高自旋态能级结构计算结果与实验数据对比图，横坐标为自旋，纵坐标为能级能量，不同颜色线条分别表示不同核素的计算值和实验值】【此处添加图1：165,167,169,171Tm高自旋态能级结构计算结果与实验数据对比图，横坐标为自旋，纵坐标为能级能量，不同颜色线条分别表示不同核素的计算值和实验值】4.2.2对关联特性分析通过理论计算，深入分析了165,167,169,171Tm高自旋态下的对关联特性，重点研究了对能随自旋的变化以及对关联在不同能级的分布情况。在对能随自旋的变化方面，计算结果表明，随着自旋的增加，165,167,169,171Tm的对能均呈现出逐渐减小的趋势，如图2所示。以165Tm为例，在低自旋态下，对能约为[X21]MeV，随着自旋增加到15-，对能减小至[X22]MeV。这是因为在高自旋态下，原子核的转动会产生科里奥利力，该力会破坏核子之间的对关联，使得对能降低。科里奥利力的作用类似于一个离心力，它会使核子之间的相对运动发生改变，从而削弱对关联的强度。这种对能随自旋的变化趋势在167,169,171Tm中也表现得较为明显，说明这是稀土区原子核在高自旋态下的一个普遍特征。对关联在不同能级的分布也呈现出一定的规律。在低能级态下，对关联主要集中在费米面附近的核子对上，这些核子对的能量较低，对关联强度较大。随着能级的升高，对关联逐渐向高能级的核子对扩展，但强度逐渐减弱。在165Tm的某一激发态能级上，位于费米面附近的核子对的对关联强度为[X23]，而在高能级的核子对的对关联强度仅为[X24]。这是因为在高能级态下，更多的核子被激发到较高的能级，核子之间的相互作用变得更加复杂，对关联受到的干扰也更大，导致对关联强度降低。对关联特性对165,167,169,171Tm高自旋态的性质有着重要影响。对能的减小会导致原子核的转动惯量增大，这是因为对关联的减弱使得核子之间的相对运动更加自由，从而增加了转动惯量。对关联在不同能级的分布也会影响原子核的能级结构和电磁性质。由于对关联强度的变化，能级之间的间距和跃迁几率也会发生改变，进而影响原子核的电磁跃迁等性质。对关联在165,167,169,171Tm高自旋态中起着关键作用，深入研究其特性对于理解这些原子核的高自旋态性质具有重要意义。【此处添加图2：165,167,169,171Tm对能随自旋变化图，横坐标为自旋，纵坐标为对能，不同颜色线条分别表示不同核素的对能变化】【此处添加图2：165,167,169,171Tm对能随自旋变化图，横坐标为自旋，纵坐标为对能，不同颜色线条分别表示不同核素的对能变化】4.2.3转动惯量与角动量分析通过理论计算得到了165,167,169,171Tm高自旋态的转动惯量和角动量，并对其随自旋的变化规律进行了深入分析，同时与理论预期和实验数据进行了对比。计算结果显示，165,167,169,171Tm的转动惯量随自旋的增加呈现出先缓慢增加，然后在某一自旋值附近突然增大，之后又逐渐趋于平缓的变化趋势，如图3所示。以165Tm为例，在低自旋态下，转动惯量约为[X25]MeV-1，随着自旋逐渐增加，转动惯量缓慢上升。当自旋达到12-左右时，转动惯量出现明显的增大，达到[X26]MeV-1，之后随着自旋的进一步增加，转动惯量的增加趋势逐渐变缓。这种变化趋势与理论预期相符，在低自旋态下，原子核主要表现为集体转动，转动惯量相对较小且变化缓慢。随着自旋的增加，处于费密面附近的核子对受到科里奥利力的作用而拆对，每个核子的角动量沿核自旋角动量的方向排列，使沿着自旋方向的角动量增加，相当于等效转动惯量突然增大，这就是所谓的回弯现象。在回弯现象发生后，随着自旋的继续增加，原子核逐渐进入全顺排状态，转动惯量的变化趋于平缓。将转动惯量的计算结果与实验数据进行对比，发现在低自旋态下，计算值与实验值较为接近。在自旋为8-时，165Tm转动惯量的计算值为[X27]MeV-1，实验测量值为[X28]MeV-1，相对偏差约为[X29]%。但在高自旋态下，尤其是在回弯现象发生后，计算值与实验值存在一定偏差。在自旋为16-时，计算值为[X30]MeV-1，实验测量值为[X31]MeV-1，相对偏差达到[X32]%。这可能是由于理论模型在描述高自旋态下原子核的复杂结构和相互作用时存在一定的局限性，如对原子核形状的动态变化、多体相互作用等因素考虑不够全面。在角动量方面，计算得到的165,167,169,171Tm高自旋态的角动量随自旋的增加而线性增加，符合理论预期。这是因为在高自旋态下，原子核的转动主要由角动量的增加驱动，随着自旋的增大，原子核的角动量也相应增大。将角动量的计算结果与实验数据进行对比，发现两者基本一致，说明理论模型能够较好地描述165,167,169,171Tm高自旋态的角动量特性。综合来看，165,167,169,171Tm高自旋态的转动惯量和角动量特性与理论预期和实验数据在一定程度上相符，但在高自旋态下，转动惯量的计算结果与实验值存在一定偏差，需要进一步改进理论模型，以更准确地描述这些原子核在高自旋态下的转动特性。【此处添加图3：165,167,169,171Tm转动惯量随自旋变化图，横坐标为自旋，纵坐标为转动惯量，不同颜色线条分别表示不同核素的转动惯量变化】【此处添加图3：165,167,169,171Tm转动惯量随自旋变化图，横坐标为自旋，纵坐标为转动惯量，不同颜色线条分别表示不同核素的转动惯量变化】五、理论与实验结果对比验证5.1对比分析方法与指标将基于粒子数守恒方法的理论计算结果与实验结果进行对比分析时，采用了一系列严谨且科学的方法，并依据多个关键指标来评估两者的符合程度。在对比能级能量时，首先将理论计算得到的165,167,169,171Tm高自旋态的能级能量与实验测量的能级能量进行一一对应。对于每一个自旋态下的能级，计算其理论值与实验值之间的绝对偏差和相对偏差。绝对偏差能够直观地反映出理论值与实验值在能量数值上的差异，而相对偏差则更能体现这种差异在实验值中所占的比例，从而更准确地评估两者的接近程度。对于165Tm在自旋为10-时的某一能级，计算得到理论能级能量为[X1]MeV，实验测量值为[X2]MeV，其绝对偏差为|X1-X2|MeV，相对偏差为(|X1-X2|/X2)×100%。通过对多个自旋态下能级能量偏差的计算和统计，能够全面地了解理论计算在描述能级能量方面的准确性。在确定自旋指定的一致性方面，主要依据实验中通过DCO比值等方法确定的能级自旋和宇称，与理论计算得到的能级自旋和宇称进行对比。在实验中，当DCO值在特定范围内时，对应着不同的γ跃迁多极性，进而确定能级的自旋和宇称。理论计算则根据构建的波函数和哈密顿量，通过角动量投影等方法得到具有确定角动量（即自旋）和宇称的能级。将两者进行对比，统计自旋指定一致的能级数量占总能级数量的比例，以此来评估理论计算在确定自旋方面的准确性。若在167Tm的能级中，共有N个能级，其中理论计算与实验自旋指定一致的能级有M个，则自旋指定一致的比例为(M/N)×100%。对于电磁跃迁几率，同样将理论计算值与实验测量值进行对比。实验中通过测量γ射线的强度和相关的符合事件，能够确定不同能级之间的电磁跃迁几率。理论计算则基于波函数和电磁相互作用的理论，计算出电磁跃迁几率。对比两者时，采用相对偏差作为评估指标，计算理论电磁跃迁几率与实验测量值之间的相对偏差，以判断理论计算对电磁跃迁几率的描述精度。对于169Tm中某两个能级之间的电磁跃迁，理论计算的跃迁几率为[X3]，实验测量值为[X4]，则相对偏差为(|X3-X4|/X4)×100%。通过对多个电磁跃迁的对比分析，能够全面评估理论计算在描述电磁性质方面的可靠性。除了上述主要指标外，还考虑了其他一些因素，如能级的相对位置关系、转动惯量的变化趋势等。在对比能级的相对位置关系时，观察理论计算得到的能级之间的间距与实验测量的能级间距是否一致，若两者差异较大，则说明理论计算在描述能级结构的相对关系方面存在不足。对于转动惯量，对比理论计算得到的转动惯量随自旋的变化曲线与实验测量的变化曲线，分析两者在变化趋势、回弯现象出现的自旋值等方面的一致性，以进一步评估理论模型对原子核转动特性的描述能力。5.2对比结果与讨论将基于粒子数守恒方法的理论计算结果与实验测量结果进行详细对比后，发现两者在能级结构、对关联特性、转动惯量等方面既存在一致性，也存在一定差异。在能级结构方面，理论计算与实验测量在低自旋态下表现出较好的一致性。165Tm在低自旋态下，理论计算得到的能级能量与实验值偏差较小，这表明粒子数守恒方法能够较为准确地描述低自旋态下的能级结构。这是因为在低自旋态下，原子核的结构相对较为简单，粒子数守恒方法中的模型和参数能够较好地反映原子核内部的相互作用和运动状态。随着自旋的增加，在高自旋态下，两者的偏差逐渐增大。167Tm在自旋较高时，计算值与实验值的偏差明显，这可能是由于在高自旋态下，原子核的形状发生了复杂的变化，而理论模型对这种形状变化的描述不够精确，导致计算结果与实验值出现偏差。也可能是因为在高自旋态下，一些高阶相互作用的影响逐渐凸显，而模型中并未充分考虑这些因素，从而影响了能级结构的计算精度。对关联特性方面，理论计算得到的对能随自旋的变化趋势与实验预期相符，均呈现出随着自旋增加对能逐渐减小的趋势。这表明粒子数守恒方法能够正确地描述对关联在高自旋态下的变化规律，对能的减小是由于原子核转动产生的科里奥利力破坏了核子之间的对关联。但在对关联在不同能级的分布细节上，理论与实验存在一定差异。理论计算的对关联分布在某些能级上与实验测量结果不完全一致，这可能是由于理论模型在处理对关联时，对一些微观效应的考虑不够全面，或者是实验测量过程中存在一定的不确定性，影响了对关联分布的测量精度。转动惯量的对比结果显示，在低自旋态下，理论计算值与实验测量值较为接近，这说明粒子数守恒方法能够较好地描述低自旋态下原子核的转动特性。但在高自旋态下，尤其是在回弯现象发生后，两者存在一定偏差。169Tm在回弯现象发生后的转动惯量计算值与实验值偏差较大，这可能是由于理论模型在描述高自旋态下原子核的复杂结构和相互作用时存在一定的局限性，如对原子核形状的动态变化、多体相互作用等因素考虑不够全面，导致转动惯量的计算出现偏差。综合来看，粒子数守恒方法在描述165,167,169,171Tm高自旋态的性质时，在低自旋态下表现出较好的准确性，但在高自旋态下，由于原子核结构和相互作用的复杂性，计算结果与实验值存在一定偏差。未来的研究可以进一步优化理论模型，考虑更多的物理因素，如原子核形状的动态变化、高阶相互作用等，以提高对高自旋态性质的描述精度。同时，也需要进一步改进实验测量技术，提高实验数据的准确性和精度，为理论研究提供更可靠的支持。5.3理论模型的验证与改进通过对基于粒子数守恒方法的理论计算结果与实验数据的详细对比分析，验证了该理论模型在描述165,167,169,171Tm高自旋态性质方面的有效性，但也发现了存在的一些差异，需要对理论模型进行改进。在能级结构的描述上，理论模型在低自旋态下与实验结果具有较好的一致性，这表明模型中的基本假设和参数选取在一定程度上能够反映原子核在低自旋态下的内部结构和相互作用。随着自旋的增加，高自旋态下的差异逐渐显现，这说明模型在处理高自旋态时存在一定的局限性。为了提高模型对高自旋态能级结构的描述能力，可以考虑调整模型参数，进一步优化对力强度和单粒子能级等参数的取值。通过更精确的实验数据拟合，确定更合适的对力强度参数，以更好地描述对关联在高自旋态下的变化。可以考虑引入更多的物理自由度，如考虑原子核形状的动态变化，采用更复杂的形变模型来描述原子核在高自旋态下的形状变化，从而更准确地计算能级能量。在对关联特性的描述方面，虽然理论模型能够正确地反映对能随自旋的变化趋势，但在对关联在不同能级的分布细节上与实验存在差异。为了改进这一点，可以完善模型假设，进一步考虑对关联的微观机制，如引入更精确的对力模型，考虑对力在不同能级和不同核子对上的变化，以更准确地描述对关联在不同能级的分布。也可以考虑多体相互作用对关联特性的影响，在模型中加入三体或更高阶的相互作用项，以更全面地描述原子核内部的相互作用。对于转动惯量的描述，理论模型在低自旋态下表现较好，但在高自旋态下与实验值存在偏差，尤其是在回弯现象发生后。为了改进这一情况，可以优化模型中的转动描述，采用更先进的转动模型，考虑原子核在高自旋态下的动态转动过程，如考虑转动过程中原子核形状的变化以及角动量的重新分布等因素，以更准确地计算转动惯量。可以结合其他理论方法，如微观-宏观模型，将微观的粒子数守恒方法与宏观的液滴模型相结合，充分考虑原子核的宏观性质和微观结构，从而提高对转动惯量的描述精度。通过对理论模型的验证和改进，可以不断提高基于粒子数守恒方法的理论模型对165,167,169,171Tm高自旋态的描述能力，为深入研究原子核高自旋态的性质提供更可靠的理论支持。六、结论与展望6.1研究成果总结本研究运用粒子数守恒方法对165,167,169,171Tm高自旋态展开了系统且深入的研究，取得了一系列具有重要意义的成果。在实验研究方面，成功布居了165,167,169,171Tm的高自旋态，并精确测量了γ-γ符合及DCO比值。在此基础上，建立了更为完善的能级纲图，165Tm新增18条γ跃迁和15条能级，167Tm新增20条γ跃迁和17条能