eagle
EMPIRE-II

What is new? (ChangeLog)


EMPIRE-2.18
Mondovi,
FORTRAN sources,
parameter lib.,
GUI and manual


Installation script
(use shift+click)
Bash script for
installing and setting
up the whole
EMPIRE-II package

EXFOR library
size 62 Mb,
needed only for
comparison plots)

HF-BCS lev. dens.
size 43 Mb; needed
only if tabulated
HF-BCS level densities
are to be used

Manual (PS)
Manual (pdf)

Installation

See also:
Screen-shots
Tsukuba contr.
CGS11 contr.
Test cases

EMPIRE home

contact:
M. Herman


Changes in version 2.18 with respect to 2.17.1
Vienna, October 2002
  1. Bug fixed in reading levels file when preparing collective levels for ECIS (R. Capote).
  2. Bug fixed in reading Coulomb parameters of o.m.p. from the RIPL-2 file (affects charged particle channels)(R. Capote).
  3. Relativistic correction in tl.f removed (gamma=1). It is not needed since SCAT2 and ECIS take care of it (R. Capote).
  4. Bug in EMPEND, which caused skipping the last but one energy point fixed (A. Trkov).
  5. EMPEND formats branching ratios in the ENDF file (A. Trkov).
  6. xterm closes when calculations are done as it was in versions 2.17 and earlier (p.3.iii below removed) since it was inconvenient for running piece-wise and multiple calculations(M. Herman).
  7. gamma-strength function uses a temperature consistent with the current level densities for the final state (not for the initial one as erroneously coded before) (M. Herman).
  8. Beta version of the new GUI (Xrun.tcl) added (M. Herman).
  9. Corrected renormalization of the absorption cross section in DEGAS (E. Betak).
  10. Corrected way of restricting exciton configurations used in DEGAS (E. Betak).
  11. Maximum number of excitons used in DEGAS set to 5 (3 configurations) (M. Herman).
  12. PLOTC4 updated by A. Trkov to produce gamma-spectra (continuum contribution from the capture reaction has not yet been included).
  13. Bugs in PLOTC4 fixed by A. Trkov. PLOTC4 tested on Linux with g77 and Absoft90 ver. 7.0 compilers and on Ms Windows with Lahey 75.
  14. If discrete levels are inconsistent with Gilbert-Cameron level densities and LEVFIT is non-negative calculations are continued using nuclear temperature taken from the systematics. Appropriate massage is printed and cumulative plot displayed (without blocking further calculations) (M. Herman).
  15. Bug in empire-specific level densities at very high energies (close to 200 MeV) when level density exceeds capacity of the 32-bit word has been fixed by ignoring fit to the discrete levels (unfortunately in all nuclei) (M. Herman).
  16. File Ko in tl.f is rewound to avoid messages that particular omp is not found in the local file (M. Herman).
  17. Number of iterations when fitting discrete levels with Gilbert-Cameron level densities limited to 300 to avoid the possibility of infinite loop (M. Herman).
  18. Z=105 name changed to 'Db'; thanks to Erik Strub for pointing it out (M. Herman).
  19. Problem with the discontinuity of neutron capture cross section when incident energy becomes lower than the integration bin width has been solved. Capture cross sections, even at very low energies (<10 keV) can now be calculated even with as few as 50 energy bins (M. Herman).
  20. All collective levels are used in ECIS calculations with DIRECT=1,2, even if they can not be excited because of the too low incident energy. In the previous versions only open channels were considered and coupling to the closed ones was ignored (R. Capote).
  21. Minor corrections of nuclear masses in EMPEND avoid potentially infinite loop while the ENDF formatted file is processed by SIXTAB (A. Trkov).
  22. Thresholds for particle emission channels corrected in EMPEND (A. Trkov).
  23. Discrete level library updated to the RIPL-2 version of September 24, 2002).
  24. LSTTAB corrected to allow for plotting spectra at 90deg with ZVView (A.Trkov).
Known bugs:
  • Fission is wrong when HRTW is active; will be corrected along with the upgrade of the fission channel
  • Slight inconsistency with DIRECT=1,2,3 option - last coll. level(s) is included in the fusion renormalization but may be excluded from actual calculations as being above the last accepted level.
To do:
  • improve fission channel (under way)
  • adjust energy for the ddx spectra in EMPEND to account for the recoils
  • clean mixed character and numerical data in commons in tl.f
Changes in version 2.17.1 with respect to 2.17
Vienna, April 2002
  1. EXFOR library converted from DOS to UNIX format.
  2. Bug in reading RIPL-2 discrete levels fixed by R. Capote.
  3. All scripts modified: (i) all links made symbolic to allow for using network drives, (ii) ./ added in front of calls (thus ./ is no longer needed in the PATH), (iii) xterm remains open after calculations are completed so that runtime messages can be inspected (M. Herman).
  4. GUI background color changed to gray (M. Herman).
  5. guizvv.tcl for graphical comparison of up to three calculations/evaluations (already present in 2.17) modified (nonstandard Tcl/Tk widget replaced) (M. Herman).
Changes in version 2.17 with respect to 2.17beta
Vienna, February 2002
  1. Bug in adding CN contribution to angular distributions for discrete levels corrected. Gammas were added to CSAlev(.,.,0) in HF-comp.f with 0 being out of dimension. Adding of gammas is now blocked. Thanks to V. Plujko for pointing it out.
  2. PLOTC4 updated by A. Trkov to produce spectra (including double-differential ones) for outgoing protons and alphas. Structure of the PostScript files improved.
  3. RIPL files placed in the empire/RIPL-2 directory following RIPL-2 structure and naming convention (compatibility with the RIPL-2 CD-ROM).
  4. Optical model database again updated to the preliminary version form RIPL-2
  5. EXFOR database updated to the version of 18.05.2001 (NOTE: previous EXFOR database will NOT work due to a different directory structure).
  6. File with levels (14) closed at the end.
  7. The new 'setup-emp' script assists installation procedure and allows to use EXFOR library and/or HF-BCS level densities directly from the CD-ROM without copying them onto the hard disk.
  8. EMPIRE manual accessible from the GUI ("Help -> EMPIRE manual" button)
Changes in version 2.17beta with respect to 2.16.2
Vienna, November 2001
  1. Dispersive optical model potential implemented (SCAT2, ECIS, ORION3) (R. Capote)
  2. Optical model database updated to the preliminary version form RIPL-2 (including 4 dispersive omp) (R. Capote)
  3. Discrete levels database updated to the preliminary version of RIPL-2 (R. Capote)
  4. Number of discrete levels constituting a complete scheme taken from the preliminary version of RIPL-2 (R. Capote)
  5. In all cases elastic is calculated on the 2.5 deg grid while all the rest on the 10 deg grid. Previous inconsistency in elastic transfer when DIRECT>0 was used removed (M. Herman)
  6. GUI starts with gvim editor selected by default (M. Herman)
  7. Two compilation warnings in tl.f fixed (M. Herman)
  8. Cumulative plots of discrete levels created as a PostScript file rather than being dumped onto the screen. Stored as *-cum.ps and available through GUI (M. Herman)
  9. Damping of the rotational level density enhancement below BCS critical energy made consistent with the one above critical energy. Removes discontinuity at critical energy for perfectly spherical nuclei (def=0) (M. Herman)
  10. List of warnings available through GUI (M. Herman)
  11. Tabulated HFBCS levels densities as provided to RIPL-2 by S. Goriely included as additional level density option (M. Herman)

Changes in version 2.16.2 with respect to 2.16
Vienna, 21 August 2001

  1. Implementation of the HMS model completed (needs testing)
  2. New version of V. Zerkin viewer (zvv94l.exe) provides eps format and interaction with plot titles (use A + right shift)
  3. Double use of file 33 removed (36 employed in lev-dens.f)

Changes in version 2.16 with respect to 2.15
Vienna, July 2001

  1. EXFOR retrieval using FORTRAN instead of UNIX 'grep'
  2. Fission barrier set to 1000 MeV for Z<19
  3. GDR second hump energy set to 0 for spherical nuclei
  4. Temperature in the generalized Lorenzian protected against 0 excitation energy
  5. Plotc4.f improved (A. Trkov) to allow for up to 23 data sets to be drawn on a single plot. Plotc4 input modified to include whole energy range.
  6. Orion calculations performed up to 99% of the maximum energy loss (instead of 90% that happen to cause problems with extrapolation to low outgoing energies)
  7. Plotc4.f improved (A. Trkov) to allow for energy spectra and angular distributions plots.
  8. Legend and lsttab added for future improvements but not fully implemented yet.
  9. ENDF=2 option added (lumped channel representation MT=10 or 5). SO FAR DDX ARE NOT SUPPORTED BY EMPEND!!!!
  10. Exact treatment of recoil spectra with ENDF=2 option.
  11. Call to Monte Carlo preequilibrium model HMS coded by M.B. Chadwick. NO TRANSFER OF THE RESULTS!!!!
  12. Exciton model code DEGAS implemented by P. Oblozinsky
  13. Fixed bug in flux conservation when GST option selected but MSC gamma channel closed.
  14. New version of SCAT2 introduced by R. Capote
  15. Link to ECIS for calculation of the elastic, absorption and cross sections to discrete collective levels within rotational and vibrational CC model. Includes automatic selection of collective levels (R. Capote).
  16. Link to the omp segment of RIPL introduced by R. Capote. NOTE: change in the format and name of the local file with internal omp. New name is *-omp.int and two columns are added in definition of om potentials (just before radii). RIPL potentials are stored in the file *-omp.ripl.
  17. New version of V. Zerkin viewer (zvv93l.exe)
  18. List of RIPL omp added to the GUI under Help menu
  19. DWBA option with ECIS added by R. Capote
  20. pol20, pol21 and pol22 variables in SCAT2 initialized with 0.
  21. Bug in ACCUMSD when MSD transitions to continuum were energetically closed fixed (seems had no effect on the results).
  22. Total cross section added to EXFOR retrieval and together with fission to the MT list for ZVView plotting.
  23. 'zvpl' script added for plotting excitation function for any cross section from the main EMPIRE output (*.lst)
  24. run-piece-wise' script added for running up to 3 different inputs for the same case in the 3 non-overlapping energy ranges (incident energies inside and among the 3 ranges must increase monotonically).
  25. Division by 0 in fitting field parameters in MSD-tristan.f protected.
  26. GUI modified to include 'zvpl' script and allow for merging ZVV plots.
  27. Form factor for the l=0 transitions in MSD controlled from the optional input (default: standard surface form factor)
  28. Provision for use of combined preequilibrium models in a single run.
  29. 'zvcomb' script added to combine various existing ZVView plots.
  30. Utility code manuals and Frequently Asked Questions (FAQ) added to the GUI under Help menu
  31. Utility code inputs accessible from the GUI Options menu
  32. Utility code c4sort added to sort *.c4 file in ascending energy. Actually, it is disabled (commented) in run and runE scripts as it takes long for large files.

Changes in version 2.15 with respect to 2.14.1
Vienna, November 2000

  1. Plotting of double differential cross sections for neutron production using PLOTC4 (modified by A. Trkov) added. This involves modification of EXFOR retrieval (new REAC_SIG.TXT file, changes in the 'sel' script and input.f), as well as modifications in the util/x4toc4/reaction file.
  2. Printout of total and shape elastic cross sections, shape elastic angular distributions, strength functions, and scattering radius from SCAT-2 activated.
  3. Minimum Tl set to 1.0E-10 to avoid underflow in the calculations. This limit might be changed in the line 761 of tl.f

Changes in version 2.14.1 with respect to 2.14
Vienna, October 2000

  1. Bug in flux conservation at very low incident energies (in HRTW) fixed
  2. Mismatch of elastic channels for negative parity targets (in HRTW) fixed
  3. Division by zero 1p-0h level density in MSC gamma emission fixed
  4. Undefined variables in MSC-NVWY.f and MSD-tristan.f defined (KASE in MSD-orion.f still remains undefined, Z1 and Z2 do not make any problem since unused)
  5. MATIN1 replaced by MTXINV in MSC gamma emission
  6. VMS specific statements marked with *-IF VMS in input.f and io.h
  7. Print of C.M. incident energy to standard output added in main.f

Changes in version 2.14 with respect to 2.13
Vienna, September 2000

  1. Second chance preequilibrium emission applying M.B. Chadwick model to the result of the MSD
  2. Bug in clearing the population of discrete levels when ENDF option used fixed
  3. Bug in writing optical model parameters for the incident channel on the ompar file fixed.
  4. SHC(2)=SHC(1) bug in main.f fixed
  5. Miss-print for 95-Zr in data/ldp.dat(7676) fixed by setting the value to 0
  6. Redundant FLOAT in ORION2 removed
  7. Collective level densities with a=A/LEVDEN used whenever LEVDEN.GT.2.0
  8. Star format for recoil mass replaced by the fixed one
  9. Alphanumeric branching ratios in orsi.liv replaced by numbers and input.f changed appropriately (BCDNUM removed)
  10. HRTW formulation of the statistical model included (width fluctuation correction).
  11. Bug causing use of Weisskopf estimates for gamma strength at all incident energies except the first one fixed.
  12. Bug in passing GDR parameters through input removed.
  13. Zerkin's zvd plotting capability added (through zvd script and Tcl/Tk GUI).
  14. EXFOR entries for neutron capture coded as SIG,,AV (averaged cross sections) included in plots.

Changes in version 2.13 with respect to 2.12.2
Vienna, March, 2000

  1. General clean-up of the FORTRAN source
  2. New organization of modules that reflects physical contents
  3. Check whether MSD cross section is not larger than the fusion cross section
  4. Gamma cascade in the first CN is default if CN excitation energy is below 20 MeV and can be controlled with GCASC in optional input
  5. l=0 transfer strength in MSD-orion set to the self-consistent value by default (rather than to GMR energy as before)
  6. Fusion cross sections read from the file name.fus are treated as fusion cross sections at subsequent l values (starting with l=0) rather than fusion cross sections to a given spin of CN. Transmission coefficients are calculated from the read values and a proper angular momentum and parity coupling is performed to obtain Compound Nucleus spin distribution.
  7. Error in BCS blocking in TRISTAN removed
  8. BCS blocking in TRISTAN made automatic (no input required)
  9. pipe.c replaced with the one by R. Capote
  10. Automatic retrieval of EXFOR data
  11. Link to EMPEND to create ENDF formatted file
  12. Link to X4TOC4 and PLOTC4 to produce plots at the end of the run
  13. Fission barriers for Z>102 according to Myers & Swiatecki Phys. Rev. C60 1999
  14. FITLEV creates a whole set of cumulative plots without interaction
  15. Fusion barrier can be read from input (BFUS)
  16. Bass option for fusion disabled
  17. Distributed fusion barrier model now relies on the CCFUS fusion barrier
  18. OMPAR (.omp) file with optical model parameters produced and used for input

Changes in version 2.12.2 with respect to 2.12.1
Vienna, 11 April 1999

  1. Energy shift used to fit low energy discrete levels gradually decreased between Ecut and neutron binding to recover agreement with Dobs.
  2. Entire code checked with FTNCHEK and some errors corrected (undefined variables, unreferenced ones removed).
  3. All DO loops given separate endings (END DO, no labels)
  4. Entire code consistently indented.

Changes in version 2.12.1 with respect to 2.12
Vienna, 07 Dec. 1998

  1. source directory added with explicit double precision source (no compiler -r8 option required).
  2. roign.f routine replaced by roemp.f. It fits dynamical level densities (BCS approach below critical energy) to discrete levels by applying an energy shift. ATTENTION: for the time being these shifts destroy agreement with Dobs at neutron binding.

Changes in version 2.12 with respect to 2.11
Vienna, 15 Sept. 1998

  1. Error in the determination of pairing corrections in readnix.f (in the case of Gilbert-Cameron) removed.
  2. Errors in coding of the dynamical level densities removed.
  3. New, EMPIRE specific, systematics of level density parameter 'a' for dynamical level densities introduced and set as default if dynamical level density chosen. This systematics fits Ilijnov-Mebel D-observed using EMPIRE formulas for l.d. If an experimental value for a given nucleus is present in the ldp.dat file it is taken. Average normalization factor is calculated for these cases and applied for the remaining ones for which there is no experimental value. All nuclei are treated as deformed as far as the collective enhancement is concerned. Temporarily damping of the collective effects has been removed, i.e. level densities are calculated in the adiabatic approach.