Equation of State
=====================

The EoS input file **is required** if the options **use_beta_equilibrium** or **use_charge_neutrality** are enabled. 
It must contain the following columns, presented in order:

+----------------------------------------------------+------------------+
| **Variable**                                       | **Unit**         |
+====================================================+==================+
| Temperature (T)                                    | MeV              |
+----------------------------------------------------+------------------+
| Baryon chemical potential (:math:`\mu_B`)          | MeV              |
+----------------------------------------------------+------------------+
| Strange chemical potential (:math:`\mu_S`)         | MeV              |
+----------------------------------------------------+------------------+
| Electric charge chemical potential (:math:`\mu_Q`) | MeV              |
+----------------------------------------------------+------------------+
| Baryon density (:math:`n_B`)                       | fm\ :sup:`-3`    |
+----------------------------------------------------+------------------+
| Strange density (:math:`n_S`)                      | fm\ :sup:`-3`    |
+----------------------------------------------------+------------------+
| Electric charge density (:math:`n_Q`)              | fm\ :sup:`-3`    |
+----------------------------------------------------+------------------+
| Energy density (:math:`\varepsilon`)               | MeV/fm\ :sup:`3` |
+----------------------------------------------------+------------------+
| Pressure (:math:`P`)                               | MeV/fm\ :sup:`3` |
+----------------------------------------------------+------------------+
| Entropy density (:math:`s`)                        | 1/fm\ :sup:`3`   |
+----------------------------------------------------+------------------+


If your model distinguishes between the 'physical' baryon density (:math:`n_B= \partial P/\partial \mu_B`) and 
the total baryon density resulting sum of particle densities (:math:`n_B^{particles}= \sum_i B_i n_i`), please pass the physical baryon density in the fifth EoS file, and append the particle baryon density in the last (11th) column, in fm\ :sup:`-3`.

This happens, for example, in models with a potential that is an explicit functions of :math:`\mu_B`. This is crucial if neutrinos are used, as the code will calculate the lepton fraction from the particle baryon density (:math:`Y_L= n_L / n_B^{particles}`), if provided.