RESNUCLEi

Summary

Scores stopping nuclei on a region basis. ​

​

See also DCYSCORE, DCYTIMES, IRRPROFIle, RADDECAY)

WHAT(1) :

type of products to be scored

= 1.0:

spallation products​ (all inelastic interactions​ except those induced by neutrons below the threshold for multigroup treatment)

= 2.0:

low-energy neutron products​, i.e. those produced by neutrons below the threshold for multigroup treatment (provided the information is available, see Note 1).

= 3.0:

all residual nuclei are scored (if available, see above)

= 0.0:

resets the default (= 1.0)

Default

= 1.0  (only spallation products are scored)

WHAT(2) =

logical output unit

> 0.0:

formatted data are written on the WHAT(2) unit

< 0.0:
unformatted data are written on the |WHAT(2)| unit
Values of |WHAT(2)| < 21.0 should be avoided
(with the exception of WHAT(2) = +11.0 = standard output​).
= 0.0:

resets the default (= 11.0)

Default

= 11.0  (standard output unit)

WHAT(3) =

maximum atomic number Z of the residual nuclei distribution

Default

:  according to the Z of the element(s) of the material assigned to the scoring region

WHAT(4) =

maximum M = N - Z - (NMZ)\(_{min}\) of the residual nuclei distribution (see Notes 2 and 3)

Default

:  according to the mass and atomic number A, Z of the element(s) of the material assigned to the scoring region

WHAT(5) =

scoring region number or name

= -1.0:

all regions (see Note 9)​

Default

= 1.0

WHAT(6) =

volume of the region in cm\(^3\) (or, more in general, a normalization factor by which the scoring shall be divided)

Default

= 1.0

SDUM

= any character string identifying the detector (max. 10 characters)


Notes

  1. Elements or isotopes for which the FLUKA low-energy neutron cross sections contain information on the production of residual nuclei are indicated by “Yes” in column 5 (“Residual nuclei”) of Table 10.3 where the components of the neutron cross section library are listed.

    The same information can be obtained by setting the printing flag in the LOW–NEUT option (WHAT(4)> 0.0). If such data are available for a given nuclide, the following message is printed on standard output:

    (RESIDUAL NUCLEI INFORMATIONS AVAILABLE)

  2. To minimise storage, nuclei are indexed by Z (with Z\(_{min}\)= 1) and NMZ = N - Z (with (NMZ)\(_{min}\)= -5). The parameter M is defined as M = NMZ - (NMZ)\(_{min}\): therefore M\(_{min}\)= 1. The following relations can also be useful:

    N - Z = M + (NMZ)\(_{min}\) N = M + Z + (NMZ)\(_{min}\)

  3. In the case of heavy ion projectiles the default NMZ, based on the region material, is not necessarily sufficient to score all the residual nuclei, which could include possible ion fragments.

  4. In order to achieve reasonable results for residual nuclei production the new evaporation​ module must be activated (it is currently the default) and heavy fragment evaporation should also be activated (it is not the default beacuse of the related large CPU penalty). Coalescence should also be activated (see option PHYSICS for all these settings). The old evaporation​ is still available, mostly because of historical reasons, but it does not produce meaningful results for residuals. The new evaporation, available since 1997, is far more sophisticated in this respect, while differences in emitted particle spectra and multiplicities are weak.

  5. Starting with Fluka2006.3 protons are scored together with \(^2\)H, \(^3\)H, \(^3\)He, \(^4\)He, at the end of their path, if transported​ (see option IONTRANS,​). This is a change with respect to previous versions where protons were not scored.

  6. All residual nuclei are scored when they have been fully de-excited down to their ground​ or isomeric state​.

  7. Radioactive decay​ of residual nuclei can be performed by FLUKA in the same run: see commands DCYTIMES, DCYSCORE, IRRPROFIle and RADDECAY, or can be done off-line by a user-written code (see for instance the program Usrsuwev available with the normal FLUKA distribution). If command IRRPROFI has been issued, RESNUCLEi results provided by detectors associated to a cooling time index via DCYSCORE will be expressed in Bq.

  8. An example on how to read RESNUCLEi unformatted output is shown below. An explanation of the meaning of the different variables is given in the comments at the beginning of the program. The program lists the Z and A of the produced nuclei, followed by the corresponding amount per unit volume.

    A more complex program Usrsuw​, which allows to compute also standard deviations over several runs, is available with the normal FLUKA code distribution in directory $FLUPRO/flutil.
    A special version of the same program, Usrsuwev, provides in addition a calculation of induced activity and of its evolution in time.

  9. Setting WHAT(5) = -1 will provide the sum of the residual nuclei in all regions, divided by the value set by the user for WHAT(6).

      PROGRAM RDRESN
*---------------------------------------------------------------------*
*     Up to MXRSNC user-defined track or coll are allowed             *
*            izrhgh = maximum Z of the scoring (minimum Z: 1)         *
*            imrhgh = maximum M=N-Z-NMZ_min of the scoring            *
*                    (minimum M: 1). Note:                            *
*                     N-Z = M + NMZ_min, N = M + Z + NMZ_min          *
*            itursn = type of binning: 1 = spallation products,       *
*                     2 = low-energy neutrons products,               *
*                     3 = all products                                *
*            nrursn = region                                          *
*            vursnc = volume (cm**3) of the detector                  *
*            tiursn = scoring name                                    *
*---------------------------------------------------------------------*
      PARAMETER ( MXRSNC = 400 )
      CHARACTER*10 TIURSN
      CHARACTER RUNTIT*80, RUNTIM*32, FILNAM*80

      DIMENSION TIURSN(MXRSNC), ITURSN(MXRSNC), NRURSN(MXRSNC),
     &          VURSNC(MXRSNC), IMRHGH(MXRSNC), IZRHGH(MXRSNC)
      DIMENSION RNDATA(MXRSNC,100,260)

      WRITE(*,*)' Type the name of the input file:'
      READ (*,'(A)') FILNAM
      LQ = INDEX(FILNAM,' ') - 1
      OPEN (UNIT=1, FILE=FILNAM, STATUS='OLD', FORM='UNFORMATTED')
      OPEN (UNIT=2, FILE=FILNAM(1:LQ)//'.txt', STATUS='UNKNOWN')
*----------- read and write 1st record ---------------------------------
      READ (1) RUNTIT, RUNTIM, WEIPRI, NCASE
      WRITE(2,100) RUNTIT, RUNTIM, NCASE, WEIPRI
*------- loop on residual nuclei detector data in the present file -----
      DO 1 IRN = 1, MXRSNC
         READ (1, END=1000) NRN, TIURSN(NRN), ITURSN(NRN),
     &   NRURSN(NRN), VURSNC(NRN), IMRHGH(NRN), IZRHGH(NRN), K
         IF (ABS(ITURSN(NRN)) .LE. 1) THEN
            WRITE (2,200) NRN, TIURSN(NRN), NRURSN(NRN),
     &        VURSNC(NRN), IZRHGH(NRN), IMRHGH(NRN) + K, K + 1
         ELSE IF (ABS(ITURSN(NRN)) .LE. 2) THEN
            WRITE (2,300) NRN, TIURSN(NRN), NRURSN(NRN),
     &        VURSNC(NRN), IZRHGH(NRN), IMRHGH(NRN) + K, K + 1
         ELSE
            WRITE (2,400) NRN, TIURSN(NRN), NRURSN(NRN),
     &        VURSNC(NRN), IZRHGH(NRN), IMRHGH(NRN) + K, K + 1
         END IF
         WRITE(2,'(/,A)') '   Z   A   Residual nuclei'
         WRITE(2,'(A,/)') '         per cm**3 per primary'
         READ (1) ((RNDATA(NRN,I,J), I=1,IZRHGH(NRN)), J=1,IMRHGH(NRN))
         DO 2 I = 1, IZRHGH(NRN)
            DO 3 J = 1, IMRHGH(NRN)
               IF(RNDATA(NRN,I,J) .GT. 0.)
     &         WRITE(2,'(2I4,1P, G15.6)') I, J+K+2*I, RNDATA(NRN,I,J)
  3         CONTINUE
  2      CONTINUE
  1   CONTINUE
 1000 CONTINUE
 100  FORMAT(/,1X,'*****',2X,A80,2X,'*****',/,/,10X,A32,/,/,
     &        10X,'Total number of particles followed ',I9,', for a ',
     &        'total weight of ',1P,E15.8,/)
 200  FORMAT (/,3X,'Res. nuclei n. ',I3,'  "',A10,
     &        '" , "high" energy products, region n. ',I5,
     &        /,6X,'detector volume: ',1P,E11.4,' cm**3',/
     &          6X,'Max. Z: ',I3,', Max. N-Z: ',I3,' Min. N-Z:',I3)
 300  FORMAT (/,3X,'Res. nuclei n. ',I3,'  "',A10,
     &        '" , "low"  energy products, region n. ',I5,
     &        /,6X,'detector volume: ',1P,E11.4,' cm**3',/
     &          6X,'Max. Z: ',I3,', Max. N-Z: ',I3,' Min. N-Z:',I3)
 400  FORMAT (/,3X,'Res. nuclei n. ',I3,'  "',A10,
     &        '" ,  all         products, region n. ',I5,
     &        /,6X,'detector volume: ',1P,E11.4,' cm**3',/
     &          6X,'Max. Z: ',I3,', Max. N-Z: ',I3,' Min. N-Z:',I3)
      END

Example:

* Calculate residual nuclei produced in an iron slab (region 6) and in a zinc
* vessel (region 10). Heavy recoils are transported (option IONTRANS) and scored
* at the point where they stop. The new evaporation model is activated to ensure
* a better quality of the results. For iron, all residual nuclei are scored. For
* zinc, no data are available for low-energy neutrons, so only nuclei produced
* by spallation/evaporation are scored. Results are written (formatted) on
* logical unit 22 and 23, respectively.
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
MATERIAL        26.0       0.0      7.87       11.       0.0       0. IRON
MATERIAL        30.0       0.0     7.133       12.       0.0       0. ZINC
ASSIGNMAT       11.0       6.0       9.0       0.0   !  Four Fe slabs
ASSIGNMAT       12.0      10.0       0.0       0.0   !  Zn vessel
IONTRANS        -2.0
PHYSICS          2.0       0.0       0.0       0.0       0.0       0. EVAPORAT
RESNUCLEI        3.0      22.0       0.0       0.0       6.0       0. FirstFe
RESNUCLEI        1.0      23.0       0.0       0.0      10.0       0. Znvessel