Output

YAML output file

The user can name this file in the YAML node “yaml_output_file” of the input file. This file is stored in the repo home directory and shows info of the format of the output files.

EoS output files

The user can specify the folder name where they want to store the EoS output files in the yaml node “output_path” of the input file. By default, the folder is named “output”.

Inside the folder, there are several different output files.

phi0Phi1_lines.csv: This table is directly calculated from the initial conditions given by phi0 Phi1 values.

eos.csv: This table is interpolated based on the output file phi0Phi1_lines.csv and thus has a uniform grid in \(T\) and \(\mu_B\). It has all the phases: stable, unstable, and meta phases.

stable_eos.csv: This table is obtained by performing a Maxwell construction on the output file eos.csv, and thus only has the stable phase.

transition_line.csv: This table records the first-order transition line and relevant thermodynamics given by the EoS. The calculation is not precise enough yet, and the first-order transition line does not go all the way to the critical point.

spinodal_lines.csv: This table records the spinodal lines and relevant thermodynamics.

Transport and energy loss output

When transport_coefficient: true is set in transport_options, transport quantities are appended as extra columns to phi0Phi1_lines.csv, spinodal_lines.csv, eos.csv, and (after Maxwell construction) stable_eos.csv. No separate transport file is generated. Points where transport values are unphysical (NaN or negative intermediate results) are output as NaN in all tables.

Column metadata files

Every output table is written together with a companion metadata file named <table>_metadata.yaml (for example, eos.csv is accompanied by eos_metadata.yaml). The metadata file lists the columns of the table in the order in which they appear in the data file, and for each column gives its description, data type, and units:

T:
  description: temperature
  type: number
  x-unit: MeV
mu_B:
  description: baryon chemical potential
  type: number
  x-unit: MeV
# ... one entry per column, in column order

This is the metadata layout used by the MUSES Porter library, so a “muses-synthesis” table (comma-separated, no header row) can be loaded directly with:

from muses_porter import Porter

porter = Porter()
porter.import_table(
    "output/eos.csv",
    extension="CSV",
    filename_metadata="output/eos_metadata.yaml",
    delimiter=",",
)

The metadata files are generated for both the “muses-synthesis” and the “legacy” output formats.

Output file format

If one selects the table format to be “muses-synthesis” (the default format), the equation of state tables follow the following format:

Variable	Unit
Temperature (T)	MeV
Baryon chemical potential (\(\mu_B\))	MeV
Strange chemical potential (\(\mu_S\))	MeV
Electric charge chemical potential (\(\mu_Q\))	MeV
Baryon density (\(n_B\))	MeV³
Strange density (\(n_S\))	MeV³
Electric charge density (\(n_Q\))	MeV³
Energy density (\(\varepsilon\))	MeV⁴
Pressure (\(P\))	MeV⁴
Entropy density (\(s\))	MeV³

Since the model currently only allows one charge, columns related to Q and S charges will be set to 0. If one specifies muses_extra_columns, those columns will be appended at last.

phi0Phi1_lines, transition_line, and spinodal_lines are not equation of state tables, so they will not follow the same format. When transport_coefficient: true, transport columns (qhat, sigma_B, zeta, and per-velocity drag/kappa columns) are appended after the standard columns in phi0Phi1_lines and spinodal_lines. The format for each file is listed below:

For phi0Phi1_lines and spinodal_lines:

Variable	Unit
Temperature (T)	MeV
Baryon chemical potential (\(\mu_B\))	MeV
Entropy density (\(s\))	MeV³
Baryon density (\(n_B\))	MeV³
Energy density (\(\varepsilon\))	MeV⁴
Pressure (\(P\))	MeV⁴
Trace Anomaly (\(I\))	MeV⁴
Phase	N/A
Baryon susceptibility (\(chi_2\))	N/A
(\(\frac{\partial s}{\partial T}\))	N/A
(\(\frac{\partial n_B}{\partial \mu_B}\))	N/A
(\(\frac{\partial s}{\partial \mu_B}\))	N/A
(\(\phi_0\))	N/A
(\(\Phi_1\))	N/A
precision	N/A

For transition_line:

Variable	Unit
Baryon chemical potential (\(\mu_B\))	MeV
Temperature (T)	MeV
Pressure (\(P\))	MeV⁴
Latent heat (\(L\))	MeV⁴

When transport_coefficient: true, the following columns are appended to the EoS output (eos.csv and stable_eos.csv). One set of drag and Langevin columns is added per velocity value specified in transport_options:

Variable	Unit
Jet quenching parameter (\(\hat{q}\))	N/A
Baryon conductivity (\(\sigma_B/T\))	N/A
Bulk viscosity (\(\zeta/s\))	N/A
Drag force (\(F_\mathrm{drag}\)) per velocity	N/A
Langevin coefficient (\(\kappa_\perp\)) per velocity	N/A
Langevin coefficient (\(\kappa_\parallel\)) per velocity\| N/A