Physics Overview
The thermodynamic properties of the Leptons are calculated from the Lagrangian formalism of a free Fermi gas:
where \(\psi_i\) and \(\nu_{Li}\) are the electron, muon, tau, and neutrino fields, respectively. :math m_i is the mass of the lepton and \(\gamma_{\mu}\) are the Dirac matrices.
The thermodynamic properties are calculated at zero temperature, so the Fermi-Dirac distribution simplifies to a step function and the integrations are made analytically. The density, energy density and pressure are given by:
and
where \(\gamma_i=2\) for electrons and \(\gamma_i=1\) for neutrinos, and \(k_F\) is the Fermi momentum, defined from \(\mu_i=\sqrt{k_F^2 + m_i^2}\).
The three principal functionalities of the module are:
Generate a Leptonic Equation of State (EoS):
The user can generate a gridded EoS exclusively with leptons by setting the option use_lepton_eos to true in the configuration. The table is generated assuming equilibrium between massive leptons, such that the chemical potentials of the electron, muon and tau are equal to \(\mu_e\). The same is true for the neutrino flavors, which have a chemical potential \(\mu_\nu\). The EoS is generated by looping over the chemical potentials these chemical potentials (\(\mu_e, \mu_\nu\)) and their range is flexible and can also be adjusted in the input configuration.
Make the nuclear grid electrically neutral:
For charge neutral matter, the code takes each row in the input EoS file and uses a root-solver from the Ceres library to find the chemical potential(s) (\(\mu_e\) and \(\mu_\nu\)) that ensure a zero net charge density and, if neutrinos are enabled, that satisfy the lepton fraction that was given in the input configuration. If there are rows in the input EoS file for which there is no solution, the row is excluded from the output (for example, when if nuclear charge density is negative at zero temperature). To enable the charge neutrality option, the user must set the option use_charge_neutrality to true. The lepton charge fraction and the lepton fraction are defined as:
and
where \(n_i\) is the number density of particle species \(i\), \(n_{Q, l}\) is the lepton charge density and \(n_B\) is the baryon number density.
Find the beta-equilibrated EoS:
This functionality employs the same `Ceres <http://ceres-solver.org/>`_root-solver to locate points satisfying the bemath:ta-equilibrium condition for chemical potentials. This entails ensuring:
If neutrinos are not present, then \(\mu_\nu=0\), and the electric charge potential of the baryons becomes the negative of the proton \(\mu_Q = -\mu_e\). \(n_{Q, b}\) is the baryon charge density. To enable the beta-equilibrium calculation, the user must set the option use_beta_equilibrium to true.
For the calculation, the module checks if the baryon chemical potential (\(\mu_B\)) is equally spaced. If so, it will secctionalize the 2d input data into vectors of 1d data of fixed \(\mu_B\) and attempt to find beta-equilibrium for each value. By default, if \(\mu_B\) is not equally spaced, it will asume that the data can be sectionalized in the baryon density. If irregular grids are anticipated, the user can enable the use_multidimensional_interpolator option.
Throughout the functionalities above, the user has the freedom to choose the types of leptons to include—electrons, muons, taus, and their corresponding neutrinos. Be aware that if neutrinos are chosen, the lepton fraction must be specified in the input configuration.