This version includes new aspects that improve the computation of the counting efficiency for each one of the three available atomic rearrangement detection models (i.e., KLM, KL 1L 2L 3M and KLMN). The first modification involves a correction algorithm that simulates the non-linear response of the detector to photoionization for low-energy X-ray photons. Although this correction has the inconvenience of substantially increasing the number of atomic rearrangement detection pathways, the computed counting efficiency for low- Z nuclides is reduced by 2% for moderate quenching in agreement with experiment. The program also simulates how the addition of extra components, such as a quencher or aqueous solutions, affects the counting efficiency. Since the CIEMAT/NIST method requires identical ionization quench functions for the electron-capture nuclide and the tracer, the computation of the counting efficiency for 3H, the low-energy beta-ray emitter commonly used as tracer, is included in the program as an option. Program summary Title of program:EMILIA Catalogue identifier:ADWK Program summary URL: http://cpc.cs.qub.ac.uk/summaries/ADWK Program obtainable from: CPC Program Library, Queen's University of Belfast, N. Ireland Licensing previsions: none Computers: revisions: any IBM PC compatible with 80386 or higher Intel processors Operating systems under which the program has been tested:MS-DOS and higher systems Programming language used:FORTRAN 77 Memory required to execute with typical data: 253 kword No. of bits in a word: 32 No. of lines in distributed program, including test data, etc.:7147 No. of bytes in distributed program, including test data, etc.:74 776 Distribution format:tar.gz Nature of the physical problem: The determination of radioactivity in liquid samples of electron-capture nuclides is demanded in radiation physics, radiation protection, dosimetry, radiobiology and nuclear medicine. The CIEMAT/NIST method has proved to be suitable for radionuclide standardizations when the counting efficiency of the liquid-scintillation spectrometer is sufficiently high. Although the method has widely been applied to beta-ray nuclides, its applicability to electron-capture nuclides nowadays has not the required degree of accuracy. The inaccuracies of the method are mainly induced by the huge number of low-energy electrons and X-ray photons emitted by the atomic rearrangement cascade after the electron-capture process, which are efficiently detected by liquid-scintillation counting, but are of difficult modeling due to the inherent complexity of the atomic rearrangement process and the non-linear response of the spectrometer in the low-energy range. Solution method: A detailed simulation of the non-linear response of the spectrometer to photoionization must include the radiation emitted by the atomic rearrangement cascade. However, a model considering all possible scintillator de-excitations at atomic level increases exponentially the number of atomic rearrangement detection pathways subsequent to capture. Since the contribution of the non-linear effects to the counting efficiency are only corrective, we can approximate the reduced energy involved in the photoionization process to a sum of only three terms: the photoelectron energy, the mean energy of the KXY Auger electrons and the global energy contribution of the remaining radiation (electrons and X-rays) emitted by the atomic rearrangement cascade. The value of each term depends on the nature of the atom in which photoabsorption is produced and on the atom shell from where the photoelectron is ejected. The non-linear correction required to simulate low-energy X-ray photoabsorption is basically important when the scintillator cocktail contains elements of high atomic numbers. For heavy atoms, the K- and L-shell binding energies can be slightly less than the energy of the colliding photon. For such cases, the non-linear effects can play an important role. Restrictions on the complexity of the problem: The simulation of all possible detection pathways of atomic rearrangement that follow to photoionization complicates the problem unnecessary. To correct the non-linear effects we consider only significative photoelectric interactions with the K- and L-shells. Also we assume that K-shell photoionizations only generate three types of entities: the photoelectron itself, KXY Auger electrons and a radiation group that includes the remaining emitted particles. For L-shell photoelectrons the radiation emission subsequent to photoionization is considered as a whole.