Penelope/Penepma/Penfluor/Fanal Monte-Carlo Software and Physics

[Probeman]

It appears that the NBTH and NBPH parameters are based on the Penepma 2012 "bilayer.in" file which is the default input file for thin film calculations. Apparently the boys in Barcelona must have changed this parameter in Penepma 2016.

So the best solution is to simply edit the bilayer.in file for the new parameter. I've attached a new version below that should work with Penepma 2016 and will upload a new Penepma16.zip file with this modification.

All I did was remove these two lines:

NBTH   45                     [No. of bins for the polar angle THETA]
NBPH   30                   [No. of bins for the azimuthal angle PHI]

And replace them with this line:

NBANGL 45 30              [Nos. of bins for the angles THETA and PHI]

That should work.

[John Donovan]

OK, we've uploaded a ZIP file containing executable and support files for Penepma 2016. You can access the link on this page:

https://probesoftware.com/Technical.html

or you can just grab it with this direct link:

https://probesoftware.com/download/PENEPMA16.ZIP

Here are the short instructions:

Penepma 2016 is almost identical with Penepma 2012, except that it provides support for emission lines of Li, Be and B. However, Penepma 2016 is not distributed with CalcZAF/Standard because the physics database files for Penepma 2016 are *not* compatible with the Penfluor/Fanal secondary fluorescence application files. Therefore Penepma 2016 is a separate download (here) which can be unzipped and placed in a folder named UserData\Penepma16, if desired for Li, Be and B emission calculations. To utilize the Penepma GUI in the Standard application one should edit the following enties in their Probewin.ini file:

    PENDBASE_Path="C:\Userdata\Penepma16\Pendbase"
    PENEPMA_Path="C:\Userdata\Penepma16\Penepma"
    PENEPMA_Root="C:\Userdata\Penepma16"
    PENEPMA_PAR_Path="C:\Userdata\Penepma12\Penfluor"

The above path to the Penfluor files isn't changed because the Penfluor and Fanal programs won't work with the Penepma 2016 physics database files. 

We probably should re-do the code relating to the Penfluor path in CalcZAF, Standard and PFE, but right now if you want to run the Penfluor/Fanal GUI as described here:

https://probesoftware.com/smf/index.php?topic=58.0

for secondary fluorescence calculations at boundaries, you'll need to re-edit the Penepma paths in the Probewin.ini file back to Penepma12.

In the meantime this Penepma16.ZIP file contains sub folders of output files for pure Li, Be and B at 5, 10, 15, 20 and 25 keV. See the Penepma16\penepma\Batch folder.

[jrminter]

Many thanks, John. You put a lot of work into this...

[Probeman]

Many thanks, John. You put a lot of work into this...

Thank-you John.

I should also mention that in order to calculate intensities down to 50 eV (necessary for Li ka), the number of secondaries generated are too large for the "forcing factor" arrays as specified by default in the Penepma GUI.

For example, here are the default forcing factors generated by selecting the "Optimize Production of Characteristic X-rays" option (from the Cu_Cha.in file):

       >>>>>>>> Interaction forcing.
IFORCE 1 1 4 -10    0.1 1.0           [KB,KPAR,ICOL,FORCER,WLOW,WHIG]
IFORCE 1 1 5 -400   0.1 1.0           [KB,KPAR,ICOL,FORCER,WLOW,WHIG]

It's the FORCER factor in red (-400) that causes the problem, so if you intend to calculate characteristic spectra down to 50 eV or so you should edit this number to -200 as seen here:

       >>>>>>>> Interaction forcing.
IFORCE 1 1 4 -10    0.1 1.0           [KB,KPAR,ICOL,FORCER,WLOW,WHIG]
IFORCE 1 1 5 -200   0.1 1.0           [KB,KPAR,ICOL,FORCER,WLOW,WHIG]


Just FYI.
john