Nonmonotonic EOS when including hyperons in CMF

Calculation Engine Support
,
1

Hi All,

I’ve submitted some jobs using CMF-lepton-QLIMR, as well as Crust DFT, CEFT, CMF- lepton- QLIMR with hyperons enabled (`use_hyperons: true`) and nonzero values for m_3H. However, these runs produce non-monotonically increasing EOS.

I’ve tried setting m_3H to the default value (0.85914584), and also tested values{0.05, 0.5, 1}, but the issue persists in all of them.

Has anyone encountered this before or found a potential cause of this issue?

Thanks,

Ania

|461px;x285px;

2

Thanks for using the CE and seeking help here. I’m guessing that @carlosc7 and @cnc6 may have some ideas, as they are the developers of QLIMR and CMF, respectively.

3

Hi Ania, how are you connecting these 3 EoS’s? 1st order phase transition, mixed phase, a particular smoothing function? Maybe @mrpelicer can help.

4

Hi,

I’m quite certain the issue is that the CMF table goes down to muB<945 MeV, reaching the liquid-gas phase transition, and then, when Lepton computes beta-equilibrium by muB it finds the wrong solution due to the interpolator being built with multiple phases at a single muB.

Turning on the following options in Lepton should fix the problem:

config:
  multidimensional_interpolator:
    use_multidimensional_interpolator: true
    number_of_points: 100 #change to your preference

This way the interpolator runs in (nB, muQ), and it avoids having multiple phases per muB.

5

It looks like I need to perform some cluster maintenance today, which may cause intermittent connection issues and may abort running jobs. If so, accept my apologies in advance and be prepared to launch the jobs again until they complete properly.

6

Yes, it solves the problem. Thank you!

7

@aliu I just enabled the “Mark as solved” plugin for this category. Would you like the honor of marking the first solution? You should see a checkbox on the posts now that you can use to select Mateus’s reply post.

8

Done :smiley: Thanks!