[论文解读] Density dependence of the "symmetry energy" in the lattice gas model

Q. M. Su, Y. G.+4 more|arXiv (Cornell University)|Nov 1, 2007

Statistical Mechanics and Entropy被引用 5

一句话总结

本研究利用同位旋依赖的格点气体模型，研究了非对称核物质中对称能的密度依赖性。通过分析具有不同同位旋不对称性、温度和冻结-out密度的热化源的同位旋标度参数，作者提取了对称能，并发现其与密度呈幂律依赖关系，对于特定的核子-核子势，C_sym = 30(ρ/ρ₀)^0.62，与近期的实验估计一致。

ABSTRACT

Isoscaling behavior of the statistical emission fragments from the equilibrated sources with $Z$ = 30 and $N$ = 30, 33, 36 and 39, resepectively, is investigated in the framework of isospin dependent lattice gas model. The dependences of isoscaling parameters $α$ on source isospin asymmetry, temperature and freeze-out density are studied and the "symmetry energy" is deduced from isoscaling parameters. Results show that "symmetry energy" $C_{sym}$ is insensitive to the change of temperature but follows the power-law dependence on the freeze-out density $ρ$. The later gives $C_{sym}$ = 30$(ρ/ρ_0)^{0.62}$ if the suitable asymmetric nucleon-nucleon potential is taken. The effect of strength of asymmetry of nucleon-nucleon interaction potential on the density dependence of the "symmetry energy" is dicussed.

研究动机与目标

• 利用格点气体模型研究非对称核物质中对称能在温度和密度上的依赖性。
• 研究同位旋标度参数 α 和 β 如何随源的同位旋不对称性、温度和冻结-out密度变化。
• 从同位旋标度行为中提取对称能 C_sym，并评估其对核子-核子相互作用势的敏感性。
• 评估核子-核子势中不对称性对对称能密度依赖性的影响。
• 将模型推导出的对称能与近期的实验和理论估计进行比较。

提出的方法

• 同位旋依赖的格点气体模型模拟了 A = 60、Z = 30 以及 N = 30、33、36、39 的平衡核源，覆盖不同的温度和冻结-out密度。
• 蒙特卡罗梅特罗波利斯采样将核子放置在三维立方格点上，指定同位旋（s_i = ±1 表示质子/中子，0 表示空位）。
• 哈密顿量包含同位旋依赖的相互作用：ε_pn = -5.33 MeV（吸引力），通过调节 ε_pp 和 ε_nn 来调节势能不对称性。
• 通过动能阈值判断碎片形成：若相对动能不足以破坏键，则形成簇。
• 从同位素产额比 R21(N,Z) = Y2(N,Z)/Y1(N,Z) = C exp(αN + βZ) 中提取同位旋标度参数 α 和 β。
• 从 α 推导出对称能 C_sym，并通过拟合 C_sym = C₀(ρ/ρ₀)^γ 来提取 C₀ 和 γ，以研究其密度依赖性。

实验结果

研究问题

• RQ1 在格点气体模型中，对称能 C_sym 如何随冻结-out密度 ρ 变化？
• RQ2 对称能对温度和源同位旋不对称性的变化有多敏感？
• RQ3 核子-核子势的不对称性（ε_pn - ε_pp）如何影响对称能的密度依赖性？
• RQ4 该模型在多大程度上再现了经验上建议的 C_sym ∝ (ρ/ρ₀)^0.62 形式？
• RQ5 同位旋标度参数能否可靠地从热化、平衡的核源中提取对称能？

主要发现

• 对称能 C_sym 对温度的依赖性较弱，但对冻结-out密度表现出强烈的幂律依赖性，形式为 C_sym = C₀(ρ/ρ₀)^γ。
• 对于势参数组（ε_pn = -7 MeV，ε_pp = 0 MeV），拟合结果为 C₀ ≈ 30 MeV 且 γ ≈ 0.62，与近期实验估计一致。
• 随着核子-核子势不对称性增加（即 |ε_pn - ε_pp| 越大），C₀ 值随之增大，表明对称能状态方程变得更硬。
• 所有测试的势参数配置下，同位旋标度行为均保持一致，证实了该方法的鲁棒性。
• 模型在不同温度和密度下均再现了观测到的同位旋标度定律 R21(N,Z) = C exp(αN + βZ)。
• 推导出的对称能对核子-核子相互作用参数的选择敏感，尤其是 ε_pn 与 ε_pp 之间的差异。