下载地址: download
相关软件: related software
Study on the proton-induced spallation reactions and the equation of state for asymmetry nuclear matter(1765)重离子熔合反应和强阻尼反应的输运理论研究(266)Au+Au三分裂反应的输运理论研究(646)重离子碰撞中的同位旋及近阈K产生(460)质子引起的散裂反应及介质内的核子-核子有效相互作用和两体散射截面(1083)
软件简介: brief introduction
改进的量子分子动力学(ImQMD)程序
用途: 模拟近垒以及中能重离子碰撞, 研究核反应的微观动力学过程.
This program (ImQMD) is to simulate heavy-ion collisions at low and intermediate energies,
based on the improved quantum molecular dynamics (ImQMD) model.
version:2.2
Nov. 22, 2015
Corresponding authors: Ning Wang, Zhuxia Li
Emails: wangning@gxnu.edu.cn (N. Wang)
lizwux@ciae.ac.cn (Z. Li)
Website: http://www.ImQMD.com/code/
References:
Ning Wang, Kai Zhao, and Zhuxia Li, Sci. China Ser. G - Phys. Mech. Astron. 58 (2015) 112001.
Ning Wang, Li Ou, Yingxun Zhang, and Zhuxia Li, Phys. Rev. C 89 (2014) 064601.
Ning Wang, Zhuxia Li, Xizhen Wu, Junlong Tian, Yingxun Zhang, and Min Liu, Phys. Rev. C 69 (2004) 034608.
attached two exe files "tt.exe" and "pot2.1.exe" are to see the movements of particles and to calculate the density distribution.
As an example, "fusion_16O92Zr.dat" presents the calculated fusion cross sections for 16O+92Zr with "ImQMD2.2.exe"
under the given input file (with Windows Operation System).
INPUT FILE:
================================
The first line: IFERMI,IFLAG1,IFLAG2,IFLAG3
---------------------------------------------------
iFermi : Fermi Constraint
---> 0 : without FermiC
---> 1 : with FermiC
iflag1 : collsion control
---> 0 : without collision
---> 1 : reduced collision
---> 2 : full collision
---> 3 : Pauli-block from PaPa
iflag2 : reaction type
---> 0 : test
---> 1 : fusion
---> 2 : multifragmentation
---> 3 : Qausi-fission
---> 4 : Strongly damped
iflag3 : record xyz and pxyz
---> 0 : DO NOT record x,y,z and px,py,pz
---> N : call output at each N-time-steps, record x,y,z and px,py,pz
================================
The second line: A2, Z2, A1, Z1 (mass and charge numbers of target and projectile)
================================
The third line: NT,DT, nrun
---------------------------------------------------
NT: total time steps to simulate
DT: each time step (default value: 1fm/c)
nrun: number of simulation events at each energy and each impact parameter
================================
The fourth line: bmin,bmax,db
---------------------------------------------------
bmin: minimal impact parameter
bmax: maximal impact parameter
db: step size for impact parameter
default values for fusion reaction: bmin=1, bmax=12, db=1
================================
The fifth line: EcmL,EcmH,dEcm
---------------------------------------------------
EcmL: minimal incident energy at center-of-mass frame
EcmH: maximal incident energy at center-of-mass frame
dEcm: step size for incident energy
================================
The sixth line: ALPHA,BETA,GAMMA,symtry,GRADI,gtau,eta,cks,rho0,c0,c1,gtaus, ct1, ct2, ct3
---------------------------------------------------
model parameters adopted for Hamiltonian
gtaus: the coefficient of the symmetry potential for the non-linear density dependence part (default value 0)
ct1, ct2, ct3: the corresponding cofficients of Thomas-Fermi kinetic energy, which don't influence the simulation process.
================================
The seventh line: d00, iseed
---------------------------------------------------
d00: initial distance at beam direction (z-axis)
(default value 30fm for light and intermediate systems and 40fm for heavy systems)
iseed: intial random seed (iseed<0, to repeat selected events, with the file seed.dat)
================================
The eighth line: Nrepeat, FCC
---------------------------------------------------
Nrepeat: The maximal times to re-perform the two-body elastic scattering in Fermi constraint at each time
FCC: Transfer factor for inelastic scattering adopted in Fermi constraint
(For heavy-ion collisions, the value of FCC should be adjusted to guarantee the conservation of energy of the system)
-------------------------------------------------------------------------------
Example I: for fusion reaction 16O+92Zr
1, 0, 1, 0
92, 40, 16, 8
2000, 1, 100
1, 12, 1
60, 40, -2
-0.207, 0.138, 1.16666, 0.034, 0.0165, 0.0140, 1.6666, 0.4, 0.165, 0.94, 0.020, 0, 0.0124, 0.001 , 0 IQ3a
30 , 0
4, 5E-6
The fusion cross sections are recorded in "fusion.dat"
------------------------------------------------------------------------------
Example II: for multifragmentation 197Au+197Au at 35AMeV and b=1fm
1, 0, 2, 500
197, 79, 197, 79
3000, 1, 1000
1, 1, 1
3447.5, 3447.5, -5
-0.207, 0.138, 1.16666, 0.034, 0.0165, 0.0140, 1.6666, 0.4, 0.165, 0.94, 0.020, 0, 0.0124, 0.001 , 0 IQ3a
80 , 2000
4, 5E-6
The charge and mass distributions are recorded in "cluster.dat"
The isotopic distributions are recorded in "isotope.dat"
------------------------------------------------------------------------------
Example III: to test the stability of an individual nucleus 92Zr at t=2000fm/c
1, 0, 0, 5
92, 40, 0, 0
2000, 1, 200
1, 1, 2
80, 80, -5
-0.207, 0.138, 1.16666, 0.034, 0.0165, 0.0140, 1.6666, 0.4, 0.165, 0.94, 0.020, 0, 0.0124, 0.001 , 0 IQ3a
40 , 0
4, 5E-6
The binding energies per particle and rms radii are given in "rms.dat"
网友评论: comments 发表评论