经典流体的计算机模拟试验——蒙特卡洛法和分子动力学法

<正> 物质是由分子、原子组成的。尽管我们往往在不考虑分子、原子的微观运动情况下,照样可以把握物质的宏观运动。如流体力学的基本方程,虽然不涉及流体本身是由什么组成,但我们仍然可以成功地运用它去解决实际问题。但最终要深刻认识宏观规律或进一步发现新规律,从物质的微观结构出发去探求其宏观运动规律,无疑是一个根本性的重要途径。正是由于这个原因,《流体力学年鉴》的主编在1980年特约G.E.Uhlenbeck写了一篇文章,专门论述了流体力学与统计物理学间的关系。然而,从微观角度研究物质

Probing balance energy using momentum- and isospin-dependent potential

Using the momentum- and isospin-dependent Boltmann-Uehling-Uhlenbeck (BUU) model, we investigate the transverse flow and balance energy in two isotopic colliding systems Ca-48+Fe-58 and Cr-48+Ni-58 by adopting different symmetry potentials. By comparing the results between the two colliding systems, we find that the difference between the balance energies of two isotopic systems can be considered as a sensitive probe to the density dependence of symmetry energy.

Effects of elastic and inelastic NN scattering cross sections on pi(-)/pi(+) ratios in heavy-ion collisions at intermediate energies

Based on the isospin-dependent Boltzmann-Uehling-Uhlenbeck transport model and the scaling model according to nucleon effective mass, effects of elastic and inelastic NN scattering cross sections on pi(-)/pi(+) in the neutron-rich reaction of Ca-48 + Ca-48 at a beam energy of 400 MeV/nucleon are studied. It is found that cross-section effects of both NN elastic and inelastic scatterings affect Delta(1232), pi(-) and pi(+) production, as well as the value of pi(-)/pi(+).

Can the nuclear symmetry potential at supra-saturation densities be negative?

In the framework of an isospin-dependent Boltzmann-Uehling-Uhlenbeck (IBUU) transport model, for the central Au-197 + Au-197 reaction at an incident beam energy of 400 MeV/nucleon, the effect of nuclear symmetry potential at supra-saturation densities on the preequilibrium clusters emission is studied. It is found that for the positive symmetry potential at supra-saturation densities the neutron-to-proton ratio of lighter clusters with mass number A less than or similar to 3 [(n/p)(A less than or similar to 3)] is larger than that of the heavier clusters with mass number A > 3 [(n/p)(A>3)], whereas for the negative symmetry potential at supra-saturation densities the (n/p)(A less than or similar to 3) is smaller than the (n/p)(A>3). This may be considered as a probe of the negative symmetry potential at supra-saturation densities.

非相对论重离子碰撞动力学的两体关联输运理论

论文系统地介绍和分析了当前几类描述重离子碰撞动力学过程输运理论的优点和缺点，并针对这些理论的不足作为建立我们理论模型的出发点。例如BUU（Boltzmann-Uehling-Uhlenbeck）系列的输运过程，可以通过不同方法从量子多体理论依照扩展的时间相关的Hartree-Fock（ETDHF）的基本思想推导出来，并对重离子碰撞过程中有关平均场物理量给出了合理的描述。但由于推导中采用了一些半经典近似和参数化近似，破坏了平均场合碰撞项之间的动力学自洽耦合。特别是数值计算中采用了实验粒子系综平均法，从而丢失了多体关联和涨落，使其无法直接描述重离子碰撞中基本的碎块形成和多重碎裂过程。而量子分子动力学（QMD）系列的理论能够给出重离子碰撞过程中碎块形成的动力学描述，但至今无法从量子多体理论推导出有关QMD的输运方程。碰撞项是在数值计算中通过Monte Carlo抽样技术人为地加入的。那么，如何从量子多体理论出发推导出描述重离子碰撞动力学过程的，可将时间相关的平均场，多体关联进行自洽耦合描述的量子输运理论就成了本论文工作的中心目的。 基于王顺金等人建立的多体关联动力学理论，选用时间相关的相干态单粒子基矢作为新理论的工作表象，对两体关联动力学中的一体密度矩阵和两体关联函数进行轨道展开，推导出了描述非相对论重离子碰撞动力学过程的两体关联输运理论TBCTT（Two-Body Correlation Transport Theory）。TBCTT是一组包括时间相关的平均场，两体关联和Pauli原理的自洽耦合的动力学方程组。其中时间相关的相干单粒子基矢是该理论的一个关键问题。其时间演化的动力学可由多种不同的方法得到，如时间相关的Hartree-Fock方法，时间相关的变分方法等。但作为建立TBCTT工作的第一步，为了计算简便，我们采用经典的Hamilton方程来描述相干单粒子基矢中相空间参数的时间演化，然后通过与两体关联动力学的耦合而恢复TBCCT基本的量子特征。 利用TBCTT对几组轻的碰撞系统进行了数值计算和分析。计算结果表明：TBCTT可以给出重离子碰撞过程中的有关物理量时间演化过程的合理描述，得到了在不同入射道条件下与QMD模型的可比性结果。同时也在组态空间的有限截断和两体关联函数不同的等级截断下均得到了碰撞系统总动量，总能量和总粒子数的近似守恒结果，特别是能量守恒，这是一般半经典输运理论中一个重要的困难问题。另外还得到了两体关联函数不同的等级截断近似对碰撞动力学方面不同的描述。所用这些数值计算结果充分表明：TBCTT是一个有希望和有发展前途的能够描述重离子碰撞动力学的量子输运理论。最后我们对当前计算中所采用的近似和存在的问题进行了分析和讨论，提出了进一步改进和完善TBCTT的途径和方案。 我们在开始TBCTT的研究之前还对QMD进行了仔细的研究和改进。通过在平均场中引入Pauli势和对称势并利用摩擦冷却方法构造原子核基态，得到了一种改进的量子分子动力学MQMD。利用这种MQMD研究了12C+12C反应多重碎裂过程中的核结构效应，得到了与AMD和实验数据基本一致的结果。对QMD的改进和应用为开展TBCTT的研究工作创造了必要条件和准备。