Gibbs construction

Calculation Engine Support
,
1

Hello,

We are currently trying to obtain an EOS including quarks in the core with the Gibbs construction, and we did not manage to make it. We are not entirely sure which grids are required for the synthesis, whether the files coming from the CMF or lepton module, and whether those which are charge neutral, beta equilibrated or flavour equilibrated.

We tried with

  • CMF_output_for_Lepton_quarks.csv and CMF_output_for_Lepton_baryons.csv, as model_Grid_eos
  • two beta_equilibrium_eos.csv, for both quark and baryons, as model_BetaEq_eos*
  • lepton_eos.csv

and got this error message:

mapping values are not allowed here in “/scratch/9af78c44-abd0-408b-87f2-f1e362bd45f7/synthesis-1/outputs/opt/output/status.yaml”, line 2, column 49

Thanks a lot in advance,

Federico and Fernando

2

Hi Federico,

Sorry for the late reply. Making the Gibbs construction right now is indeed a bit involved. Below is an example configuration that hopefully helps you figure out the files that go into each module.

processes:
- name: cmf
  module: cmf_solver
  config:
    computational_parameters:
      solution_resolution: 1.e-12
      options:
        vector_potential: 4
        use_octet: true
        use_hyperons: true
        use_decuplet: false
        use_quarks: true
      variables:
        chemical_optical_potentials:
          muB_begin: 954.0
          muB_end: 1800.0
          muB_step: 5.0
          muQ_begin: -400.0
          muQ_end: 0.0
          muQ_step: 1.0
- name: lepton-baryons
  module: lepton
  config:
    global:
      run_name: ''
      use_beta_equilibrium: true
    output:
      output_particle_properties: true
    particles:
      use_electron: true
      use_muon: true
  pipes:
    input_eos:
      label: CMF_for_Lepton_baryons_only
      module: cmf_solver
      process: cmf
    input_particle_properties:
      label: CMF_particle_properties_baryons_only
      module: cmf_solver
      process: cmf                           
- name: lepton-quarks
  module: lepton
  config:
    global:
      run_name: ''
      use_beta_equilibrium: true
    output:
      output_particle_properties: true
    particles:
      use_electron: true
      use_muon: true
  pipes:
    input_eos:
      label: CMF_for_Lepton_quarks_only
      module: cmf_solver
      process: cmf
    input_particle_properties:
      label: CMF_particle_properties_quarks_only
      module: cmf_solver
      process: cmf
- name: lepton-eos
  module: lepton
  config:
    global:
      use_beta_equilibrium: false
      use_charge_neutrality: false
      use_lepton_eos: true
    output:
      output_particle_properties: true
    particles:
      use_electron: true
      use_muon: true
    lepton_eos_parameters:
      electron_cp_initial: 0.0
      electron_cp_final: 500.0
      electron_cp_step: 1.0
- name: synthesis
  module: synthesis
  config:
    global:
      synthesis_type: 'gibbs'
    output:
      output_derivatives: false
      output_particle_properties: true
  pipes:
    model1_BetaEq_eos:
      label: eos_beta_equilibrium
      module: lepton
      process: lepton-baryons
    model2_BetaEq_eos:
      label: eos_beta_equilibrium
      module: lepton
      process: lepton-quarks
    model1_Grid_eos:
      label: CMF_for_Lepton_baryons_only
      module: cmf_solver
      process: cmf
    model2_Grid_eos:
      label: CMF_for_Lepton_quarks_only
      module: cmf_solver
      process: cmf             
    lepton_eos:
      label: lepton_eos
      module: lepton
      process: lepton-eos
    model1_BetaEq_particle_properties:
      label: particle_properties_beta_equilibrium
      module: lepton
      process: lepton-baryons
    model2_BetaEq_particle_properties:
      label: particle_properties_beta_equilibrium
      module: lepton
      process: lepton-quarks
    model1_Grid_particle_properties:
      label: CMF_particle_properties_baryons_only
      module: cmf_solver
      process: cmf
    model2_Grid_particle_properties:
      label: CMF_particle_properties_quarks_only
      module: cmf_solver
      process: cmf
    lepton_particle_properties:
      label: lepton_particle_properties
      module: lepton
      process: lepton-eos           
components:
- type: chain
  name: cmf-gibbs
  sequence:
  - cmf
  - lepton-baryons
  - lepton-quarks
  - lepton-eos
  - synthesis