Hello everyone, I’m working on matching two equations of state under a first-order phase transition. On one side I’m using a Walecka (RMF) model with self-interactions, isospin, and beta equilibrium; on the other, a BPS-type EOS, all at zero temperature. Is there a workflow I can follow to build my code? Thanks in advance!
Hi Victor, the synthesis module is supposed to bridge different EOS, but I think it does so by using a hyperbolic tangent function, so it probably cannot create a first-order phase transition for you. @mrpelicer can correct me if that’s wrong. There is something else we are building that may create a first-order phase transition, but that is still under development.
Hi @vbruno1997,
Synthesis can compute the first-order phase transition from your two EoS files. There is an example in this notebook. See the function cmf_maxwell_wf_config.
However, the example uses input files coming from MUSES modules, not from an upload. You’ll have to modify it to use your uploaded files. Please take a look at this example to see what needs to be changed.
Let us know if you have any issues.
I should mention that Synthesis can also compute a 1st-order transition with phase coexistence. See cmf_gibbs_wf_config in the same notebook I linked above.
@mrpelicer and @yumuy2 Thank you all for your help. I will do the tests and will soon come back with the feedback.
Dear all, i would like to sincerely thank you for your assistance. I am pleased to inform you that everything worked out successfully. Although I did not use the model exactly as provided, I used it as a foundation, which proved very helpful. I truly appreciate your valuable feedback and support.
Glad to hear it. If you learned something through this experience that you think could improve the Synthesis module and/or CE documentation to help the next person, we would welcome your contribution to the docs source, either as an issue you open with your idea or (bonus points!) a merge request with suggested changes.