Class implementing cosmological equation of state as a mixture of HRG model, ideal gases of leptons and photons. More...

#include <CosmicEoS.h>

Classes
class	BroydenEquationsCosmology
	Class implementing the Broyden equations for cosmology. More...

Public Member Functions
	CosmicEoS (ThermalModelBase *THMbase, bool pionsinteract=false)
	Constructor.

virtual	~CosmicEoS (void)
	Destructor.

ThermalModelBase *	HRGModel () const
	Gets the pointer to the HRG model object.

const std::vector< double > &	ChemicalPotentials () const
	Gets the chemical potentials (baryon number, electric charge, three lepton charges).

virtual void	SetTemperature (double T)
	Set the temperature.

double	Temperature () const
	Gets the current temperature.

virtual void	SetBaryonChemicalPotential (double muB)
	Set the baryon chemical potential.

double	BaryonChemicalPotential () const
	Gets the baryon chemical potential.

virtual void	SetElectricChemicalPotential (double muQ)
	Set the electric chemical potential.

double	ElectricChemicalPotential () const
	Gets the current electric chemical potential.

virtual void	SetLeptonChemicalPotential (LeptonFlavor::Name flavor, double muL)
	Set the lepton chemical potential.

double	LeptonChemicalPotential (LeptonFlavor::Name flavor) const
	Gets the current lepton chemical potential for a specific flavor.

virtual void	SetBaryonAsymmetry (double b)
	Sets the baryon number asymmetry (baryon over entropy density).

double	BaryonAsymmetry () const
	Gets the current baryon asymmetry.

virtual void	SetChargeAsymmetry (double q)
	Sets the electric charge asymmetry (charge over entropy density).

double	ChargeAsymmetry () const
	Gets the current electric charge asymmetry.

virtual void	SetLeptonAsymmetry (LeptonFlavor::Name flavor, double l)
	Sets the lepton flavor asymmetry (lepton flavor over entropy density).

double	LeptonAsymmetry (LeptonFlavor::Name flavor) const
	Gets the current lepton flavor asymmetry.

void	SetAsymmetries (const std::vector< double > &asymmetries)
	Sets all asymmetries at once.

void	CalculatePrimordialDensities ()
	Calculates number densities of all particle species.

double	EntropyDensity ()
	Calculates the total entropy density.

double	EntropyDensityHRG ()
	Calculates the entropy density of the HRG part.

double	Pressure ()
	Calculates the total pressure.

double	PressureHRG ()
	Calculates the partial pressure of the HRG part.

double	EnergyDensity ()
	Calculates the total energy density.

double	EnergyDensityHRG ()
	Calculates the energy density of the HRG part.

double	NetDensityChargedLepton (int iL)
	Calculates the net density of charged lepton flavor.

double	PressureChargedLepton (int iL)
	Calculates the partial pressure of charged lepton flavor.

double	EnergyDensityChargedLepton (int iL)
	Calculates the energy density of charged lepton flavor.

double	BaryonDensity (bool absolute=false)
	Calculates the total baryon density.

double	ElectricChargeDensity (bool absolute=false)
	Calculates the electric charge density.

double	ElectricChargeDensityHRG (bool absolute=false)
	Calculates the electric charge density of the HRG part.

double	IsospinChargeDensity (bool absolute=false)
	Calculates the isospin charge density.

double	LeptonFlavorDensity (LeptonFlavor::Name flavor, bool absolute=false)
	Calculates the lepton flavor density.

std::vector< double >	SolveChemicalPotentials (double T, const std::vector< double > &muInit=std::vector< double >())
	Calculates the values of the chemical potential (B,Q,{L}) that satisfy the given asymmetry constraints at given temperature.

double	GetPionMass () const
	Gets the mass of pi+.

void	SetPionsInteracting (bool pionsinteract=true, double fpiChPT=fpi)
	Sets whether to include pion interactions via effective mass model.

bool	InteractingPions () const
	Checks if pions are interacting in the model.

bool	InPionCondensedPhase () const
	Checks if the system has non-zero BEC of pions.

const ThermalParticle &	PhotonParticle () const
	Gets the ThermalParticle object instance corresponding to photons.

int	NumberOfElectroWeakSpecies () const
	Gets the number of electroweak species (charged leptons and neutrinos) in the model.

std::string	GetSpeciesName (int id) const
	Gets the name of particle species of given id.

double	GetDensity (int id) const
	Gets the number density for given species.

Static Public Attributes
static double	fpi = 0.133
	Pion decay constant for pion interactions a la ChPT.

Protected Member Functions
void	ClearEMMs ()
	Clears the effective mass models.

Protected Attributes
ThermalModelBase *	m_modelHRG
	Pointer to an HRG model object.

ThermalParticle	m_Photon
	Photons.

std::vector< ThermalParticle >	m_ChargedLeptons
	Charged leptons.

std::vector< ThermalParticle >	m_Neutrinos
	Neutrinos.

bool	m_InteractingPions
	Whether to include pion interactions.

bool	m_IsCalculated
	Whether the system has been calculated.

double	m_T
	Temperature in GeV.

std::vector< double >	m_ChemCurrent
	Vector of chemical potentials (baryon, charge, lepton flavors)

std::vector< double >	m_Asymmetries
	Vector of asymmetries (baryon, charge, lepton flavors)

Detailed Description

Class implementing cosmological equation of state as a mixture of HRG model, ideal gases of leptons and photons.

Examples: CosmicTrajectory.cpp.

Definition at line 56 of file CosmicEoS.h.

Constructor & Destructor Documentation

◆ CosmicEoS()

thermalfist::CosmicEoS::CosmicEoS	(	ThermalModelBase *	THMbase,
		bool	pionsinteract = false )

Constructor.

Parameters

THMbase	Pointer to the HRG model object.
pionsinteract	Specifies whether to incorporate pions interactions through effective mass model.

Definition at line 25 of file CosmicEoS.cpp.

◆ ~CosmicEoS()

virtual thermalfist::CosmicEoS::~CosmicEoS ( void )

inlinevirtual

Destructor.

Definition at line 71 of file CosmicEoS.h.

Member Function Documentation

◆ BaryonAsymmetry()

double thermalfist::CosmicEoS::BaryonAsymmetry ( ) const

inline

Gets the current baryon asymmetry.

Returns: The baryon asymmetry value.

Definition at line 157 of file CosmicEoS.h.

◆ BaryonChemicalPotential()

double thermalfist::CosmicEoS::BaryonChemicalPotential ( ) const

inline

Gets the baryon chemical potential.

Returns: The baryon chemical potential value in GeV.

Definition at line 113 of file CosmicEoS.h.

◆ BaryonDensity()

double thermalfist::CosmicEoS::BaryonDensity ( bool absolute = false )

Calculates the total baryon density.

Parameters

absolute If true, calculates absolute baryon density.

Returns: The baryon density value.

Definition at line 221 of file CosmicEoS.cpp.

◆ CalculatePrimordialDensities()

void thermalfist::CosmicEoS::CalculatePrimordialDensities ( )

Calculates number densities of all particle species.

Definition at line 111 of file CosmicEoS.cpp.

◆ ChargeAsymmetry()

double thermalfist::CosmicEoS::ChargeAsymmetry ( ) const

inline

Gets the current electric charge asymmetry.

Returns: The charge asymmetry value.

Definition at line 171 of file CosmicEoS.h.

◆ ChemicalPotentials()

const std::vector< double > & thermalfist::CosmicEoS::ChemicalPotentials ( ) const

inline

Gets the chemical potentials (baryon number, electric charge, three lepton charges).

Returns: Vector of chemical potentials.

Definition at line 85 of file CosmicEoS.h.

◆ ClearEMMs()

void thermalfist::CosmicEoS::ClearEMMs ( )

protected

Clears the effective mass models.

Definition at line 401 of file CosmicEoS.cpp.

◆ ElectricChargeDensity()

double thermalfist::CosmicEoS::ElectricChargeDensity ( bool absolute = false )

Calculates the electric charge density.

Parameters

absolute If true, calculates absolute charge density.

Returns: The electric charge density value.

Definition at line 227 of file CosmicEoS.cpp.

◆ ElectricChargeDensityHRG()

double thermalfist::CosmicEoS::ElectricChargeDensityHRG ( bool absolute = false )

Calculates the electric charge density of the HRG part.

Parameters

absolute If true, calculates absolute charge density.

Returns: The HRG electric charge density value.

Examples: CosmicTrajectory.cpp.

Definition at line 245 of file CosmicEoS.cpp.

◆ ElectricChemicalPotential()

double thermalfist::CosmicEoS::ElectricChemicalPotential ( ) const

inline

Gets the current electric chemical potential.

Returns: The electric chemical potential value in GeV.

Definition at line 127 of file CosmicEoS.h.

◆ EnergyDensity()

double thermalfist::CosmicEoS::EnergyDensity ( )

Calculates the total energy density.

Returns: The total energy density value.

Examples: CosmicTrajectory.cpp.

Definition at line 167 of file CosmicEoS.cpp.

◆ EnergyDensityChargedLepton()

double thermalfist::CosmicEoS::EnergyDensityChargedLepton ( int iL )

Calculates the energy density of charged lepton flavor.

Parameters

iL	The lepton flavor index.

Returns: The energy density value.

Examples: CosmicTrajectory.cpp.

Definition at line 212 of file CosmicEoS.cpp.

◆ EnergyDensityHRG()

double thermalfist::CosmicEoS::EnergyDensityHRG ( )

Calculates the energy density of the HRG part.

Returns: The HRG energy density value.

Examples: CosmicTrajectory.cpp.

Definition at line 185 of file CosmicEoS.cpp.

◆ EntropyDensity()

double thermalfist::CosmicEoS::EntropyDensity ( )

Calculates the total entropy density.

Returns: The total entropy density value.

Examples: CosmicTrajectory.cpp.

Definition at line 116 of file CosmicEoS.cpp.

◆ EntropyDensityHRG()

double thermalfist::CosmicEoS::EntropyDensityHRG ( )

Calculates the entropy density of the HRG part.

Returns: The HRG entropy density value.

Examples: CosmicTrajectory.cpp.

Definition at line 134 of file CosmicEoS.cpp.

◆ GetDensity()

double thermalfist::CosmicEoS::GetDensity ( int id ) const

Gets the number density for given species.

Parameters

id	0 - photon, 1 - e+, 2 - e-, 3 - mu+, ..., 7 - nu_e, 8 - anti-nu_e, ...

Returns: The number density value.

Definition at line 375 of file CosmicEoS.cpp.

◆ GetPionMass()

double thermalfist::CosmicEoS::GetPionMass ( ) const

Gets the mass of pi+.

Used to determine if Bose condensation reached in ideal gas.

Returns: The mass of pi+ in GeV.

Definition at line 307 of file CosmicEoS.cpp.

◆ GetSpeciesName()

std::string thermalfist::CosmicEoS::GetSpeciesName ( int id ) const

Gets the name of particle species of given id.

Parameters

id	The particle species id.

Returns: The name of the particle species.

Definition at line 350 of file CosmicEoS.cpp.

◆ HRGModel()

ThermalModelBase * thermalfist::CosmicEoS::HRGModel ( ) const

inline

Gets the pointer to the HRG model object.

Returns: Pointer to the HRG model object.

Examples: CosmicTrajectory.cpp.

Definition at line 78 of file CosmicEoS.h.

◆ InPionCondensedPhase()

bool thermalfist::CosmicEoS::InPionCondensedPhase ( ) const

Checks if the system has non-zero BEC of pions.

Returns: True if the system is in a pion condensed phase, false otherwise.

Examples: CosmicTrajectory.cpp.

Definition at line 338 of file CosmicEoS.cpp.

◆ InteractingPions()

bool thermalfist::CosmicEoS::InteractingPions ( ) const

inline

Checks if pions are interacting in the model.

Returns: True if pions are interacting, false otherwise.

Definition at line 345 of file CosmicEoS.h.

◆ IsospinChargeDensity()

double thermalfist::CosmicEoS::IsospinChargeDensity ( bool absolute = false )

Calculates the isospin charge density.

Parameters

absolute If true, calculates absolute isospin density.

Returns: The isospin charge density value.

Examples: CosmicTrajectory.cpp.

Definition at line 254 of file CosmicEoS.cpp.

◆ LeptonAsymmetry()

double thermalfist::CosmicEoS::LeptonAsymmetry ( LeptonFlavor::Name flavor ) const

inline

Gets the current lepton flavor asymmetry.

Parameters

flavor The lepton flavor.

Returns: The lepton asymmetry value.

Definition at line 187 of file CosmicEoS.h.

◆ LeptonChemicalPotential()

double thermalfist::CosmicEoS::LeptonChemicalPotential ( LeptonFlavor::Name flavor ) const

inline

Gets the current lepton chemical potential for a specific flavor.

Parameters

flavor The lepton flavor.

Returns: The lepton chemical potential value in GeV.

Definition at line 143 of file CosmicEoS.h.

◆ LeptonFlavorDensity()

double thermalfist::CosmicEoS::LeptonFlavorDensity	(	LeptonFlavor::Name	flavor,
		bool	absolute = false )

Calculates the lepton flavor density.

Parameters

flavor	The lepton flavor.
absolute	If true, calculates absolute lepton flavor density.

Returns: The lepton flavor density value.

Definition at line 267 of file CosmicEoS.cpp.

◆ NetDensityChargedLepton()

double thermalfist::CosmicEoS::NetDensityChargedLepton ( int iL )

Calculates the net density of charged lepton flavor.

Parameters

iL	The lepton flavor index.

Returns: The net density value.

Examples: CosmicTrajectory.cpp.

Definition at line 192 of file CosmicEoS.cpp.

◆ NumberOfElectroWeakSpecies()

int thermalfist::CosmicEoS::NumberOfElectroWeakSpecies ( ) const

inline

Gets the number of electroweak species (charged leptons and neutrinos) in the model.

Returns: The number of electroweak species.

Definition at line 369 of file CosmicEoS.h.

◆ PhotonParticle()

const ThermalParticle & thermalfist::CosmicEoS::PhotonParticle ( ) const

inline

Gets the ThermalParticle object instance corresponding to photons.

Returns: Reference to the photon particle object.

Definition at line 362 of file CosmicEoS.h.

◆ Pressure()

double thermalfist::CosmicEoS::Pressure ( )

Calculates the total pressure.

Returns: The total pressure value.

Examples: CosmicTrajectory.cpp.

Definition at line 141 of file CosmicEoS.cpp.

◆ PressureChargedLepton()

double thermalfist::CosmicEoS::PressureChargedLepton ( int iL )

Calculates the partial pressure of charged lepton flavor.

Parameters

iL	The lepton flavor index.

Returns: The partial pressure value.

Examples: CosmicTrajectory.cpp.

Definition at line 203 of file CosmicEoS.cpp.

◆ PressureHRG()

double thermalfist::CosmicEoS::PressureHRG ( )

Calculates the partial pressure of the HRG part.

Returns: The HRG pressure value.

Examples: CosmicTrajectory.cpp.

Definition at line 160 of file CosmicEoS.cpp.

◆ SetAsymmetries()

void thermalfist::CosmicEoS::SetAsymmetries ( const std::vector< double > & asymmetries )

inline

Sets all asymmetries at once.

Parameters

asymmetries Vector of asymmetries (baryon, charge, lepton flavors).

Examples: CosmicTrajectory.cpp.

Definition at line 194 of file CosmicEoS.h.

◆ SetBaryonAsymmetry()

virtual void thermalfist::CosmicEoS::SetBaryonAsymmetry ( double b )

inlinevirtual

Sets the baryon number asymmetry (baryon over entropy density).

Parameters

b	The baryon asymmetry value.

Definition at line 150 of file CosmicEoS.h.

◆ SetBaryonChemicalPotential()

virtual void thermalfist::CosmicEoS::SetBaryonChemicalPotential ( double muB )

inlinevirtual

Set the baryon chemical potential.

Parameters

muB	Baryon chemical potential (GeV)

Definition at line 106 of file CosmicEoS.h.

◆ SetChargeAsymmetry()

virtual void thermalfist::CosmicEoS::SetChargeAsymmetry ( double q )

inlinevirtual

Sets the electric charge asymmetry (charge over entropy density).

Parameters

q	The charge asymmetry value.

Definition at line 164 of file CosmicEoS.h.

◆ SetElectricChemicalPotential()

void thermalfist::CosmicEoS::SetElectricChemicalPotential ( double muQ )

virtual

Set the electric chemical potential.

Parameters

muQ	Electric chemical potential (GeV)

Definition at line 105 of file CosmicEoS.cpp.

◆ SetLeptonAsymmetry()

virtual void thermalfist::CosmicEoS::SetLeptonAsymmetry	(	LeptonFlavor::Name	flavor,
		double	l )

inlinevirtual

Sets the lepton flavor asymmetry (lepton flavor over entropy density).

Parameters

flavor	The lepton flavor.
l	The lepton asymmetry value.

Definition at line 179 of file CosmicEoS.h.

◆ SetLeptonChemicalPotential()

virtual void thermalfist::CosmicEoS::SetLeptonChemicalPotential	(	LeptonFlavor::Name	flavor,
		double	muL )

inlinevirtual

Set the lepton chemical potential.

Parameters

muL	Lepton chemical potential (GeV)
flavor	Lepton flavor

Definition at line 135 of file CosmicEoS.h.

◆ SetPionsInteracting()

void thermalfist::CosmicEoS::SetPionsInteracting	(	bool	pionsinteract = true,
		double	fpiChPT = fpi )

Sets whether to include pion interactions via effective mass model.

Parameters

pionsinteract	If true, pion interactions are included.
fpiChPT	The pion decay constant (default: fpi).

Definition at line 315 of file CosmicEoS.cpp.

◆ SetTemperature()

void thermalfist::CosmicEoS::SetTemperature ( double T )

virtual

Set the temperature.

Parameters

T	Temperature (GeV)

Definition at line 99 of file CosmicEoS.cpp.

◆ SolveChemicalPotentials()

std::vector< double > thermalfist::CosmicEoS::SolveChemicalPotentials	(	double	T,
		const std::vector< double > &	muInit = std::vector<double>() )

Calculates the values of the chemical potential (B,Q,{L}) that satisfy the given asymmetry constraints at given temperature.

Parameters

T	Temperature in GeV.
muInit	Initial guesses of the chemical potentials.

Returns: Vector of calculated chemical potentials.

Examples: CosmicTrajectory.cpp.

Definition at line 287 of file CosmicEoS.cpp.

◆ Temperature()

double thermalfist::CosmicEoS::Temperature ( ) const

inline

Gets the current temperature.

Returns: The temperature value in GeV.

Definition at line 99 of file CosmicEoS.h.

Member Data Documentation

◆ fpi

double thermalfist::CosmicEoS::fpi = 0.133

static

Pion decay constant for pion interactions a la ChPT.

Pion decay constant.

Definition at line 321 of file CosmicEoS.h.

◆ m_Asymmetries

std::vector<double> thermalfist::CosmicEoS::m_Asymmetries

protected

Vector of asymmetries (baryon, charge, lepton flavors)

Definition at line 416 of file CosmicEoS.h.

◆ m_ChargedLeptons

std::vector<ThermalParticle> thermalfist::CosmicEoS::m_ChargedLeptons

protected

Charged leptons.

Definition at line 395 of file CosmicEoS.h.

◆ m_ChemCurrent

std::vector<double> thermalfist::CosmicEoS::m_ChemCurrent

protected

Vector of chemical potentials (baryon, charge, lepton flavors)

Definition at line 413 of file CosmicEoS.h.

◆ m_InteractingPions

bool thermalfist::CosmicEoS::m_InteractingPions

protected

Whether to include pion interactions.

Definition at line 404 of file CosmicEoS.h.

◆ m_IsCalculated

bool thermalfist::CosmicEoS::m_IsCalculated

protected

Whether the system has been calculated.

Definition at line 407 of file CosmicEoS.h.

◆ m_modelHRG

ThermalModelBase* thermalfist::CosmicEoS::m_modelHRG

protected

Pointer to an HRG model object.

Definition at line 389 of file CosmicEoS.h.

◆ m_Neutrinos

std::vector<ThermalParticle> thermalfist::CosmicEoS::m_Neutrinos

protected

Neutrinos.

Definition at line 398 of file CosmicEoS.h.

◆ m_Photon

ThermalParticle thermalfist::CosmicEoS::m_Photon

protected

Photons.

Definition at line 392 of file CosmicEoS.h.

◆ m_T

double thermalfist::CosmicEoS::m_T

protected

Temperature in GeV.

Definition at line 410 of file CosmicEoS.h.

The documentation for this class was generated from the following files:

include/CosmicEos/CosmicEoS.h
src/library/CosmicEos/CosmicEoS.cpp

Classes

Public Member Functions

Static Public Attributes

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ CosmicEoS()

◆ ~CosmicEoS()

Member Function Documentation

◆ BaryonAsymmetry()

◆ BaryonChemicalPotential()

◆ BaryonDensity()

◆ CalculatePrimordialDensities()

◆ ChargeAsymmetry()

◆ ChemicalPotentials()

◆ ClearEMMs()

◆ ElectricChargeDensity()

◆ ElectricChargeDensityHRG()

◆ ElectricChemicalPotential()

◆ EnergyDensity()

◆ EnergyDensityChargedLepton()

◆ EnergyDensityHRG()

◆ EntropyDensity()

◆ EntropyDensityHRG()

◆ GetDensity()

◆ GetPionMass()

◆ GetSpeciesName()

◆ HRGModel()

◆ InPionCondensedPhase()

◆ InteractingPions()

◆ IsospinChargeDensity()

◆ LeptonAsymmetry()

◆ LeptonChemicalPotential()

◆ LeptonFlavorDensity()

◆ NetDensityChargedLepton()

◆ NumberOfElectroWeakSpecies()

◆ PhotonParticle()

◆ Pressure()

◆ PressureChargedLepton()

◆ PressureHRG()

◆ SetAsymmetries()

◆ SetBaryonAsymmetry()

◆ SetBaryonChemicalPotential()

◆ SetChargeAsymmetry()

◆ SetElectricChemicalPotential()

◆ SetLeptonAsymmetry()

◆ SetLeptonChemicalPotential()

◆ SetPionsInteracting()

◆ SetTemperature()

◆ SolveChemicalPotentials()

◆ Temperature()

Member Data Documentation

◆ fpi

◆ m_Asymmetries

◆ m_ChargedLeptons

◆ m_ChemCurrent

◆ m_InteractingPions

◆ m_IsCalculated

◆ m_modelHRG

◆ m_Neutrinos

◆ m_Photon

◆ m_T