LISE++ Abrasion-Ablation




Contents :

1. Introduction: Heavy Ions reactions at intermediate and high energies

2. Geometrical Abrasion-Ablation model

Excitation energy

3. Ablation stage (Evaporation cascade)

4. Abrasion-Ablation & Evaporation calculator

5. Settings

Evaporation cascade
Prefragment Excitation energy
Selection of Abrasion-Ablation model to calculate cross sections
Plotting Abrasion-Ablation cross sections
Mass model choice

6. Abrasion-Ablation vs. Experimental data

"MATRIX" solution
Minimization

7. Abrasion-Ablation model development in LISE++

8. Next development steps




1. Introduction: Heavy Ions reactions at intermediate and high energies




 
Nuclear charge yields calcualted by LISE++
for different de-excitation channels
after 238U(1AGeV) abrasion on a Be-target

2. Geometrical Abrasion-Ablation model

Geometrical LISE++ Abrasion-Ablation Model is a simple theory of fragmentation based upon a two-step model [Wil87]. The abrasion process accounts for removal of nuclear matter in the overlap region of the colliding ions. An average transmission factor is used for the projectile and target nuclei at a given impact parameter to account for the finite mean-free path in nuclear matter. The ions are treated otherwise on a geometrical basis assuming uniform spheres. The surface distortion excitation energy of the projectile prefragment following abrasion of nucleons is calculated from the clean-cut abrasion formalism [Gos77].

[Gos77] J.Gosset et al., Phys.Rev. C16 (1977) 629.
[Wil87] J.W.Wilson, L.W.Towsend, F.F.Badavi, NIM B18 (1987) 225-231.

2.1. Excitation energy

Currently the code contains four methods to calculate a prefragment excitation energy (See the dialog below).

Menu ➝ Physics Models ➝ Excitation Energy of Prefragment

click on image to get a high resolution picture

Fig. The Prefragment excitation energy dialog (Click to zoom).

click on image to get a high resolution picture

Fig. The Prefragment excitation energy plots (Click to zoom).



3. : Ablation stage (Evaporation cascade)

The Ablation step is based on the fusion-evaporation model LisFus [Tar03] for fast analytical calculations of fusion residues cross sections, where the evaporation stage is treated in a macroscopic way on the basis of a master equation which leads to diffusion equations as proposed by Campi and Hűfner [Cam81], and reexamined lately by Gaimard and Schmidt [Gai91]. Level densities and decay widths are taken from the statistical analysis of Iljinov et al. [Ilj92]. The LISE evaporation model works with probability distributions as function of excitation energy taking into account initially eight possible parent and daughter channels (n, 2n, p, 2p, d, t, 3He, a), and fission and breakup de-excitation channels implemented later [Tar08f]. Analytical solution of the evaporation cascade was performed with the transport integral theory [Baz94] providing fast calculations, that allows to calculate the cross-section for nuclei far from stability which are not accessible with Monte Carlo technique.

[Baz94] D. Bazin, B. Sherrill, Phys. Rev. E 50 (1994) 4017.
[Gai91] J.-J.Gaimard, K.-H.Schmidt, Nucl.Phys. A531 (1991) 709-745.
[Ilj92] A.S. Iljinov et al., Nucl. Phys. A543 (1992) 517.
[Tar03] O.B. Tarasov, D. Bazin, NIM B 204 (2003) 174-178.
[Tar08f] O.B. Tarasov, A.C.C.Villari, NIM B 266 (2008) 4670-4673.

Menu ➝ Physics Models ➝ Prefragment Search & Evaporation options

click on image to get a high resolution picture

Fig. Prefragment Search & Evaporation options dialog (Click to zoom).

click on image to get a high resolution picture

Fig. De-excitation channels plot (Click to zoom).


4. Abrasion-Ablation / Evaporation calculator

The "Evaporation calculator" is available in the Calculations menu or via the icon in the toolbar. The "Evaporation Calculator" dialog in action shown in Fig.4.1 allows to calculate and plot production cross-sections of different nuclei as a result of deexcitation of an excited nucleus (see Fig.4.2), and visulaize the evolution of excitation functions (Fig.4.3) using plots of input and output excitation energy distributions for a given nucleus:

  • Input parent distributions
  • Nucleus de-excitation channels
  • Emitted particles spectra
  • Output daughter distributions

click on image to get a high resolution picture

Fig.4.1. The Evapoaration dialog in action (Click to zoom).

click on image to get a high resolution picture

Fig.4.2. Two-dimensional plot of calculated cross-sections. The data is represented by rectangles centered at N (neutrons) and Z (protons),
where the cross-section value depends on color and size of the rectangle.
The user can change the display modes of rectangles, and also change the scale (Click to zoom).


click on image to get a high resolution picture

Fig.4.3. Excitation functions of 40Ca produced in reaction of the 48Ca beam with a Be-target. (Click to zoom).


06/02/2001 v.4.18   Evaporation Calculator


5. Settings

Evaporation cascade

see Ablation stage (Evaporation cascade)

Prefragment Excitation energy

see Excitation energy

Selection of Abrasion-Ablation model to calculate cross sections

Menu ➝ Physics Models ➝ Production Mechanism ➝
"Settings" button (Projectile Fragmentation) ➝ "Cross Sections" tab



Plotting Abrasion-Ablation cross sections

Before plotting cross-sections with
Menu ➝ 1-D Plot ➝ Cross section distributions

preliminary it is necessary to check the "Show Abasion-Ablation in X-section plots" box in the Preference dialog:
Menu ➝ Options ➝ Preferences



Mass model choice

Menu ➝ Physics Models ➝ Production Mechanism ➝ "Database: Masses, Isomers" tab




6. Abrasion-Ablation vs. Experimental data

"MATRIX" solution

09/01/2005   v.7.05   User cross-section analysis using Abrasion-Ablation model (Matrix method)

Menu ➝ Utilities ➝ Reaction utilities ➝
User cross-sections analysis using Abrasion-Ablation model ➝ Calculation of E* vs. Sigma matrix



Fig. Dialog of User cross-section analysis using Abrasion-Ablation model (Matrix method)



Fig. Final plot combined on four plots (Local-chi2, Local-LoD, Global-chi2, Global-LoD) based on user weights.

Minimization

03/10/2019   v.11.0.45   Abrasion-Ablation minimization to describe user cross-sections

Menu ➝ Utilities ➝ Reaction utilities ➝
User cross-sections analysis using Abrasion-Ablation model ➝ Minimization of E* parameters

The new minimization utility recently developed in the LISE code allows to deduce Abrasion-Ablation model parameters from comparison of AA-calculation results with experimental cross-sections with selection one from 28 mass models distributed with the LISE++ suite. The utility is based on the levmar package using the Levenberg-Marquardt nonlinear least square algorithm [LevMar].

[LevMar] M.I.A. Lourakis, levmar: Levenberg–Marquardt nonlinear least squares algorithms in C/C++, http://www.ics.forth.gr/lourakis/levmar/, 2004.

click on image to get a high resolution picture

Fig. Dialog of User cross-section analysis using Abrasion-Ablation model (Minimization method) (Click to zoom).



Fig. Plot of logarithm ratios of Abrasion-Ablation calculations and experimental values.



Fig. Plot of chi-square values.


7. Abrasion-Ablation model development in LISE++

06/02/2001   v.4.18   Abrasion-Ablation v.1.0
06/25/2001   v.5.04   Abrasion-Ablation v.1.3
12/09/2003   v.6.04   Abrasion-Ablation revision

06/25/2001   v.5.03   Prefragment Excitation Energy
09/21/2001   v.5.08   Auto mode for the Abrasion-Ablation model
11/01/2002   v.6.01   LDM parameterizations and accuracy of cross-section calculations

09/16/2003   v.6.03   Level density calculations
09/16/2003   v.6.03   Fission channel in evaporation cascade
09/16/2003   v.6.03   Comparison between LISE calculations and experimental results : LISE Abrasion-Ablation model

09/01/2005   v.7.05   The break-up channel in the evaporation cascade
09/01/2005   v.7.05   Dissipation effects in fission
09/01/2005   v.7.05   Evaporation calculator modifications
09/01/2005   v.7.05   User cross-section analysis using Abrasion-Ablation model (Matrix method)

11/18/2005   v.7.08   Decay channel analysis
03/07/2013   v.9.04   Abrasion-Ablation update, Initial prefragment plot
03/10/2019   v.11.0.45   Abrasion-Ablation minimization to describe user cross-sections


8. LISE AA model next development plans

  • Three step Abrasion-Ablation model

  • Monte Carlo method
  • Purposes:
    ⋅ benchmarks of the "Distribtuion" method,
    ⋅ using more sophisticated event gates ,
    ⋅ application of angular momentum