Thermodynamic Databases ======================= .. image:: https://img.shields.io/badge/Reaktoro%20Tutorial-Interactive!-E66581.svg?logo=data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAFkAAABZCAMAAABi1XidAAAB8lBMVEX///9XmsrmZYH1olJXmsr1olJXmsrmZYH1olJXmsr1olJXmsrmZYH1olL1olJXmsr1olJXmsrmZYH1olL1olJXmsrmZYH1olJXmsr1olL1olJXmsrmZYH1olL1olJXmsrmZYH1olL1olL0nFf1olJXmsrmZYH1olJXmsq8dZb1olJXmsrmZYH1olJXmspXmspXmsr1olL1olJXmsrmZYH1olJXmsr1olL1olJXmsrmZYH1olL1olLeaIVXmsrmZYH1olL1olL1olJXmsrmZYH1olLna31Xmsr1olJXmsr1olJXmsrmZYH1olLqoVr1olJXmsr1olJXmsrmZYH1olL1olKkfaPobXvviGabgadXmsqThKuofKHmZ4Dobnr1olJXmsr1olJXmspXmsr1olJXmsrfZ4TuhWn1olL1olJXmsqBi7X1olJXmspZmslbmMhbmsdemsVfl8ZgmsNim8Jpk8F0m7R4m7F5nLB6jbh7jbiDirOEibOGnKaMhq+PnaCVg6qWg6qegKaff6WhnpKofKGtnomxeZy3noG6dZi+n3vCcpPDcpPGn3bLb4/Mb47UbIrVa4rYoGjdaIbeaIXhoWHmZYHobXvpcHjqdHXreHLroVrsfG/uhGnuh2bwj2Hxk17yl1vzmljzm1j0nlX1olL3AJXWAAAAbXRSTlMAEBAQHx8gICAuLjAwMDw9PUBAQEpQUFBXV1hgYGBkcHBwcXl8gICAgoiIkJCQlJicnJ2goKCmqK+wsLC4usDAwMjP0NDQ1NbW3Nzg4ODi5+3v8PDw8/T09PX29vb39/f5+fr7+/z8/Pz9/v7+zczCxgAABC5JREFUeAHN1ul3k0UUBvCb1CTVpmpaitAGSLSpSuKCLWpbTKNJFGlcSMAFF63iUmRccNG6gLbuxkXU66JAUef/9LSpmXnyLr3T5AO/rzl5zj137p136BISy44fKJXuGN/d19PUfYeO67Znqtf2KH33Id1psXoFdW30sPZ1sMvs2D060AHqws4FHeJojLZqnw53cmfvg+XR8mC0OEjuxrXEkX5ydeVJLVIlV0e10PXk5k7dYeHu7Cj1j+49uKg7uLU61tGLw1lq27ugQYlclHC4bgv7VQ+TAyj5Zc/UjsPvs1sd5cWryWObtvWT2EPa4rtnWW3JkpjggEpbOsPr7F7EyNewtpBIslA7p43HCsnwooXTEc3UmPmCNn5lrqTJxy6nRmcavGZVt/3Da2pD5NHvsOHJCrdc1G2r3DITpU7yic7w/7Rxnjc0kt5GC4djiv2Sz3Fb2iEZg41/ddsFDoyuYrIkmFehz0HR2thPgQqMyQYb2OtB0WxsZ3BeG3+wpRb1vzl2UYBog8FfGhttFKjtAclnZYrRo9ryG9uG/FZQU4AEg8ZE9LjGMzTmqKXPLnlWVnIlQQTvxJf8ip7VgjZjyVPrjw1te5otM7RmP7xm+sK2Gv9I8Gi++BRbEkR9EBw8zRUcKxwp73xkaLiqQb+kGduJTNHG72zcW9LoJgqQxpP3/Tj//c3yB0tqzaml05/+orHLksVO+95kX7/7qgJvnjlrfr2Ggsyx0eoy9uPzN5SPd86aXggOsEKW2Prz7du3VID3/tzs/sSRs2w7ovVHKtjrX2pd7ZMlTxAYfBAL9jiDwfLkq55Tm7ifhMlTGPyCAs7RFRhn47JnlcB9RM5T97ASuZXIcVNuUDIndpDbdsfrqsOppeXl5Y+XVKdjFCTh+zGaVuj0d9zy05PPK3QzBamxdwtTCrzyg/2Rvf2EstUjordGwa/kx9mSJLr8mLLtCW8HHGJc2R5hS219IiF6PnTusOqcMl57gm0Z8kanKMAQg0qSyuZfn7zItsbGyO9QlnxY0eCuD1XL2ys/MsrQhltE7Ug0uFOzufJFE2PxBo/YAx8XPPdDwWN0MrDRYIZF0mSMKCNHgaIVFoBbNoLJ7tEQDKxGF0kcLQimojCZopv0OkNOyWCCg9XMVAi7ARJzQdM2QUh0gmBozjc3Skg6dSBRqDGYSUOu66Zg+I2fNZs/M3/f/Grl/XnyF1Gw3VKCez0PN5IUfFLqvgUN4C0qNqYs5YhPL+aVZYDE4IpUk57oSFnJm4FyCqqOE0jhY2SMyLFoo56zyo6becOS5UVDdj7Vih0zp+tcMhwRpBeLyqtIjlJKAIZSbI8SGSF3k0pA3mR5tHuwPFoa7N7reoq2bqCsAk1HqCu5uvI1n6JuRXI+S1Mco54YmYTwcn6Aeic+kssXi8XpXC4V3t7/ADuTNKaQJdScAAAAAElFTkSuQmCC :target: https://mybinder.org/v2/gh/reaktoro/reaktoro-jupyter/master?filepath=tutorial%2Foverview.ipynb Thermodynamic databases allow us to define and model a chemically reactive system by providing the means for the computation of necessary thermodynamic properties (e.g., standard Gibbs energies of species, equilibrium constants of reactions). In such databases, we find a collection of chemical species and/or reactions and their accompanying data, which include substance's name and chemical formula, reaction equations, thermodynamic data and/or model parameters. .. attention:: There are many thermodynamic databases available in the literature and they are in general very different from each other. At the moment, there is no standard on how thermodynamic databases should be formatted. Some databases are based on chemical reactions and their equilibrium constants (e.g., PHREEQC databases), while others are based on substances and their model parameters for evaluation of their standard thermodynamic properties at temperature and pressure of interest (e.g., SUPCRT92 databases). Reaktoro currently supports the following thermodynamic databases: * SUPCRT92 databases; * PHREEQC databases; and * GEMS databases. SUPCRT92 Databases ------------------ .. |supcrt98| replace:: :download:`supcrt98.xml <../../databases/supcrt/supcrt98.xml>` .. |supcrt07| replace:: :download:`supcrt07.xml <../../databases/supcrt/supcrt07.xml>` .. |supcrt98-organics| replace:: :download:`supcrt98-organics.xml <../../databases/supcrt/supcrt98-organics.xml>` .. |supcrt07-organics| replace:: :download:`supcrt07-organics.xml <../../databases/supcrt/supcrt07-organics.xml>` .. |slop98| replace:: :download:`slop98.dat <../../databases/supcrt/slop98.dat>` .. |slop07| replace:: :download:`slop07.dat <../../databases/supcrt/slop07.dat>` .. _supcrt98.xml: :download:`<../../databases/supcrt/supcrt98.xml>` .. _slop98.dat: :download:`<../../databases/supcrt/slop98.dat>` .. _slop07.dat: :download:`<../../databases/supcrt/slop07.dat>` .. _supcrt07.xml: :download:`<../../databases/supcrt/supcrt07.xml>` .. _supcrt98-organics.xml: databases/supcrt/supcrt98-organics.xml .. _supcrt07-organics.xml: databases/supcrt/supcrt07-organics.xml .. sidebar:: SUPCRT92 Database Files | |supcrt98| | |supcrt98-organics| | |supcrt07| | |supcrt07-organics| The SUPCRT92 thermodynamic databases supported in Reaktoro are presented next. They contain parameters for the calculation of standard thermodynamic properties of aqueous species, gases, and minerals for temperatures 0-1000 °C and pressures 1-5000 bar. The standard properties of aqueous species are calculated using the revised *Helgeson-Kirkham-Flowers* (HKF) equations of state and, for the gases and minerals, a thermodynamic model based on Maier--Kelly heat capacity polynomial equation. .. note:: The thermodynamic databases ``supcrt98.xml`` and ``supcrt07.xml``, in XML format, were derived, respectively, from the original SUPCRT92 database files |slop98| and |slop07|. In the process, **all organic aqueous species were removed**! If you need them in your modeling problem, you should then use instead ``supcrt98-organics.xml`` and ``supcrt07-organics.xml``. .. tip:: If your problem requires an aqueous phase without organic species and you are using an automatic initialization scheme for its construction (e.g., creating an aqueous phase with all species in the database whose elements are H, O, C, or Ca), then make sure you are using one of the SUPCRT92 databases **without organic species**! Otherwise, you might end up with an aqueous phase containing an extremely long list of organic species that will only serve to decrease the performance of the calculations. .. note:: The equation of state of Wagner and Pruss (2002) is used to calculate the thermodynamic properties of water and its temperature and pressure derivatives. PHREEQC Databases ----------------- .. |Amm| replace:: :download:`Amm.dat <../../databases/phreeqc/Amm.dat>` .. |frezchem| replace:: :download:`frezchem.dat <../../databases/phreeqc/frezchem.dat>` .. |iso| replace:: :download:`iso.dat <../../databases/phreeqc/iso.dat>` .. |llnl| replace:: :download:`llnl.dat <../../databases/phreeqc/llnl.dat>` .. |minteq| replace:: :download:`minteq.dat <../../databases/phreeqc/minteq.dat>` .. |minteq.v4| replace:: :download:`minteq.v4.dat <../../databases/phreeqc/minteq.v4.dat>` .. |phreeqc| replace:: :download:`phreeqc.dat <../../databases/phreeqc/phreeqc.dat>` .. |pitzer| replace:: :download:`pitzer.dat <../../databases/phreeqc/pitzer.dat>` .. |sit| replace:: :download:`sit.dat <../../databases/phreeqc/sit.dat>` .. |wateq4f| replace:: :download:`wateq4f.dat <../../databases/phreeqc/wateq4f.dat>` Reaktoro can use PHREEQC_ as a :ref:`thermodynamic backend`, which permits us to take advantage of the rich collection of PHREEQC thermodynamic databases that are listed next. .. sidebar:: PHREEQC Database Files | |Amm| | |frezchem| | |iso| | |llnl| | |minteq| | |minteq.v4| | |phreeqc| | |pitzer| | |sit| | |wateq4f| .. _PHREEQC: https://wwwbrr.cr.usgs.gov/projects/GWC_coupled/phreeqc/ GEMS Databases -------------- Reaktoro can also use GEMS as a :ref:`thermodynamic backend` and take advantage of its databases. .. todo:: Write about GEMS databases.