Skip to content

Python-based and open source data reduction tool for the major, minor and trace element analysis of minerals as well as of fluid and melt inclusions.

License

Notifications You must be signed in to change notification settings

MABeeskow/PySILLS

Folders and files

NameName
Last commit message
Last commit date

Latest commit

Β 
Β 
Β 
Β 
Β 
Β 
Β 
Β 
Β 
Β 
Β 
Β 
Β 
Β 
Β 
Β 
Β 
Β 
Β 
Β 
Β 
Β 
Β 
Β 
Β 

Repository files navigation

Python 3.12 Β  License: LGPL v3 Β  DOI

developed by Maximilian Alexander Beeskow

PySILLS is a newly developed, Python-based open source tool for a modern data reduction of LA-ICP-MS experiments. It is focused on the compositional analysis of major, minor and trace elements of minerals (and glasses) as well as of fluid and melt inclusions. PySILLS which was initially part of a M.Sc. thesis project, is developed by Maximilian Alexander Beeskow in the work group of Prof. Dr. Thomas Wagner and Dr. Tobias Fusswinkel at RWTH Aachen University. PySILLS was inspired conceptionally by the widely-used data reduction tool SILLS.

Top Features

The following list shows some of the main features that differentiate PySILLS from alternative data reduction tools.

  • works on all common computer systems that can run Python
  • use of multiple standard reference materials in one project file
  • use of multiple internal standards in one project file
  • consideration of isotope-specific standard reference materials
  • assemblage definition
  • file-specific quick analysis
  • intuitive, fast and flexible workflow
  • multiple check-up possibilities
  • export of processed LA-ICP-MS data (e.g. intensity ratios, analytical sensitivities, etc.) for external calculations
  • many quality-of-life features

Planned Features

The following features extend PySILLS and do not replace already existing features.

  • more outlier detection algorithms
  • replacement of scattered intensity values by regression curves
  • extended language support
  • in-built geothermometry analysis
  • Jupyter notebooks for a browser-based data reduction of LA-ICP-MS experiments
  • production of a video course on YouTube

If you like this project, please consider to follow (click on the star icon above) it.


⚠️ Disclaimer

I am aware of the fact that there are still some bugs and some of them can also cause a crash of the program. But, those bugs are the reason why I have decided to publish this pre-release version, because, as a developer, I see everywhere construction sites in the code, but I would like to know which ones are the most important. So, it is possible that different users have different ideas how a mineral or fluid inclusion analysis project should be done. I have developed PySILLS with my workflow in my mind. Now, the goal is to make it also applicable for your workflow.

Let's catch the bugs and smash them! πŸ‘Š

Attention

It was necessary to apply some changes concerning the saving/loading algorithm of a project. It is possible that an imported project does not contain any information about the acquisition time of the single measurement files. It is necessary to add them manually, to save the project again and to restart it finally again. If the user does not restart the project, it is possible that further bugs occur.

πŸš€ Installation

A detailed manual about the installation of PySILLS can be found here.

PySILLS can be easily installed via the following command:

pip install PySILLS

Alternatively, it is possible to run PySILLS by a manual installation of Python and - if not already installed - the following packages:

  • numpy
  • scipy
  • pandas
  • matplotlib
  • sympy
  • (tkinter)

In some cases, it will be necessary to install tkinter manually although it should be actually already part of your python version.

Python Β  NumPy Β  SciPy Β  Pandas Β  Matplotlib Sympy Β 

Windows Β  macOS Β  Linux

Running PySILLS

PySILLS can be started by running the command python pysills_app.py . If possible, use the newest version of PySILLS.

Attention (only for Mac OS users)

There is unfortunately a bug concerning the GUI library tkinter and default Python version that is installed on a Mac. You will notice this bug for example when you click on a button and for some cases, it will work directly, and for other cases, you have to click several times on the button or after moving the window a little bit. In order to avoid this bug, there is the following workaround until Apple will fix this one day:

  1. Install the newest version of PySILLS (it will work with Python 3.12). There are several ways to do this but I would recommend to install it via Homebrew which is a package manager for Mac OS.
  2. Install all needed packages explicitly for this newly installed Python version, for example via pip. Example command for the installation of numpy for Python 3.12: python3.12 -m pip install numpy .
  3. If you run PySILLS now by the command python3.12 pysills_app.py , the previously described bug should not appear anymore!

πŸ’» Resources

PySILLS manual (installation)
PySILLS manual (mineral analysis)
PySILLS manual (fluid inclusion analysis)
PySILLS manual (melt inclusion analysis)

PySILLS on ReadTheDocs
PySILLS on Blogger
PySILLS on YouTube
PySILLS on PyPi

Blogger Β  YouTube

πŸ’­ Citing PySILLS

Coming soon ...

DOI


πŸ’Ž Mineral Analysis

PySILLS allows the major, minor and trace element analysis of minerals and glasses. The calculations are based on the measured intensity signals and sensitivities of the ICP-MS instrument.

Short Step-by-Step Manual

In addition to a more detailed manual, I would like to describe briefly here which steps are necessary for a complete data reduction of a mineral analysis project.

Data import

  1. Select the correct analysis mode. In this case, select "Mineral Analysis" below "Select Mode".
  2. Import the measurement files by clicking on "Add" below "Standard Files" and "Sample Files".

Project setup

  1. Click on "Settings". A new window will be opened.

On the left side of the settings window, you can define some default parameters and features that influence the whole dataset, for example a spike elimination.

  1. Select a default standard reference material (SRM) for the standard files and isotopes below the header "Standard Reference Material (SRM)".
  2. Define the concentration of the internal standard. For this purpose, it is possible to calculate the concentration value based on the amount of a specific oxide, for example 100 % of SiO2, on the amount of an element, by importing a csv-file containing the needed information or by defining the internal standard and its concentration manually.
  3. Click on "Run" behind "Auto-Detection" below "Default Time Window (Background)" and "Default Time Window (Matrix)" in order to get automatically detected calculation intervals for the background and matrix signal. If you are very sure when the laser started and ended, you can also define manually the time limits of the calculation windows.
  4. Click on "Apply to all" for the standard and sample files in order to apply a spike elimination on all files.

In theory, all mandatory settings were defined now, but it is of course possible to change some parameters if this is necessary. We also recommend to check if the automatically set calculation intervals for the background and matrix signal were set correctly.

Project results

After all settings were defined, it is now finally possible to calculate the results of this mineral analysis project.

  1. Click on "Results" below "Mineral Analysis" in the main window of PySILLS. A new window will be opened.

The user can now specify which results should be displayed in the table, for example the concentration values for the smoothed (spike eliminated) sample files. Of course, the final step is to create a report file that contains all results.

  1. Click on "Export Results" in order to collect and export all data in a created report file. This file contains all possible values that can be displayed in the table.

Screenshots

Coming soon ...


πŸ’Ž Fluid Inclusion Analysis

PySILLS does not only allow the major, minor and trace element analysis of (homogenous) solid phases like minerals and glasses but also from fluid inclusions.

Short Step-by-Step Manual

In addition to a more detailed manual, I would like to describe briefly here which steps are necessary for a complete data reduction of a fluid inclusion analysis project.

Data import

  1. Select the correct analysis mode. In this case, select "Fluid Inclusion Analysis" below "Select Mode".
  2. Import the measurement files by clicking on "Add" below "Standard Files" and "Sample Files".

Project setup

  1. Click on "Settings". A new window will be opened.

On the left side of the settings window, you can define some default parameters and features that influence the whole dataset, for example a spike elimination.

  1. Select a default standard reference material (SRM) for the standard files and isotopes below the header "Standard Reference Material (SRM)".
  2. Define the matrix by setting a concentration value for the internal standard of the matrix which will be used for example if the "Matrix-only Tracer" quantification method was selected.
  3. Define the settings for the quantification method, for example the "Matrix-only Tracer" algorithm.
  4. Click on "Run" behind "Auto-Detection" below "Default Time Window (Background)" in order to get automatically detected calculation intervals for the background signal. If you are very sure when the laser started and ended, you can also define manually the time limits of the calculation windows.
  5. Click on "Apply to all" for the standard and sample files in order to apply a spike elimination on all files.
  6. Open every standard and sample file and define the calculation intervals for the matrix and inclusion signal.
  7. Define the concentration values of the internal standard that was selected for the sample files (this will be used for the fluid inclusion analysis). There are several options available, for example a mass balance calculation or a calculation based on PyPitzer.

Now, all mandatory settings were defined, but it is of course possible to change some parameters if this is necessary.

Project results

After all settings were defined, it is now finally possible to calculate the results of this fluid inclusion analysis project.

  1. Click on "Results" below "Fluid Inclusion Analysis" in the main window of PySILLS. A new window will be opened.

The user can now specify which results should be displayed in the table, for example the concentration values for the smoothed (spike eliminated) sample files. Of course, the final step is to create a report file that contains all results.

  1. Click on "Export Results" in order to collect and export all data in a created report file. This file contains all possible values that can be displayed in the table.

Screenshots

Coming soon ...

πŸ“¦ PyPitzer

PyPitzer allows thermodynamic modeling of fluid inclusion systems based on the Pitzer model. The quantification of those systems requires microthermometric data as well as element/Na ratios from LA-ICP-MS experiments. PyPitzer is able to calculate precise concentrations of fluid inclusion compositions for also complex multi-element systems. It is possible to use an already fully implemented version of PyPitzer in PySILLS or to use it externally which means that the user can export the necessary LA-ICP-MS data and import then the calculated Na concentrations.

PyPitzer


πŸ’Ž Melt Inclusion Analysis

PySILLS does not only allow the major, minor and trace element analysis of (homogenous) solid phases and fluid inclusions but also from melt inclusions.

Short Step-by-Step Manual

In addition to a more detailed manual, I would like to describe briefly here which steps are necessary for a complete data reduction of a melt inclusion analysis project.

Data import

  1. Select the correct analysis mode. In this case, select "Melt Inclusion Analysis" below "Select Mode".
  2. Import the measurement files by clicking on "Add" below "Standard Files" and "Sample Files".

Project setup

  1. Click on "Settings". A new window will be opened.

On the left side of the settings window, you can define some default parameters and features that influence the whole dataset, for example a spike elimination.

  1. Select a default standard reference material (SRM) for the standard files and isotopes below the header "Standard Reference Material (SRM)".
  2. Define the matrix by setting a concentration value for the internal standard of the matrix which will be used for example if the "Matrix-only Tracer" quantification method was selected.
  3. Define the settings for the quantification method, for example the "Matrix-only Tracer" algorithm.
  4. Click on "Run" behind "Auto-Detection" below "Default Time Window (Background)" and "Default Time Window (Matrix)" in order to get automatically detected calculation intervals for the background/matrix signal. If you are very sure when the laser started and ended, you can also define manually the time limits of the calculation windows.
  5. Click on "Apply to all" for the standard and sample files in order to apply a spike elimination on all files.
  6. Open every standard and sample file and define the calculation intervals for the missing matrix and inclusion signal.
  7. Define the concentration values of the internal standard that was selected for the sample files (this will be used for the melt inclusion analysis). There are several options available, for example the quantification based on a 100 wt.% oxides normalization.

Now, all mandatory settings were defined, but it is of course possible to change some parameters if this is necessary.

Project results

After all settings were defined, it is now finally possible to calculate the results of this fluid inclusion analysis project.

  1. Click on "Results" below "Melt Inclusion Analysis" in the main window of PySILLS. A new window will be opened.

The user can now specify which results should be displayed in the table, for example the concentration values for the smoothed (spike eliminated) sample files. Of course, the final step is to create a report file that contains all results.

  1. Click on "Export Results" in order to collect and export all data in a created report file. This file contains all possible values that can be displayed in the table.

πŸ“š References

  • Heinrich, C.A., et al., 2003, "Quantitative multi-element analysis of minerals, fluid and melt inclusions by laser-ablation inductively-coupled-plasma mass-spectrometry", Geochimica et Cosmochimica Acta, 67, pp. 3473-3496, Link
  • Guillong, M., et al., 2008, "SILLS: A MATLAB-based program for the reduction of Laser ablation ICP-MS data of homogenous materials and inclusions", Mineralogical Association of Canada Short Course, 40, pp. 323-333, Link
  • Liu, Y., et al., 2024, "An integrated approach for quantifiying fluid inclusion data combining microthermometry, LA-ICP-MS, and thermodynamic modeling", Chemical Geology, 644, pp. 1-18, Link
  • Longerich, H.P., et al., 1996, "Laser Ablation Inductively Coupled Plasma Mass Spectrometric Transient Signal Data Acquisition and Analyte Concentration Calculation", Journal of Analytical Atomic Spectrometry, 11, pp. 899-904, Link
  • Pettke, T., et al., 2012, "Recent developments in element concentration and isotope ratio analysis of individual fluid inclusions by laser ablation single and multiple collector ICP-MS", Ore Geology Reviews, 44, pp. 10-38, Link

Version: Pre-Release
Last updated: 27.06.2024