Accounting for non-ideal mixing effects in the hydrogen-helium equation of state

The equation of state for hydrogen and helium is fundamental for studying stars and giant planets. It has been shown that because of interactions at atomic and molecular levels, the behaviour of a mixture of hydrogen and helium cannot be accurately represented by considering these elements separately. This paper aims at providing a simple method to account for interactions between hydrogen and helium in interior and evolution models of giant planets. Using on the one hand ab initio simulations that involve a system of interacting hydrogen and helium particles and pure equations of state for hydrogen and helium on the other, we derived the contributions in density and entropy of the interactions between hydrogen and helium particles. We show that relative variations of up to 15% in density and entropy arise when non-ideal mixing is accounted for. These non-ideal mixing effects must be considered in interior models of giant planets based on accurate gravity field measurements, particularly in the context of variations in the helium-to-hydrogen ratio. They also affect the mass-radius relation of exoplanets. We provide a table that contains the volume and entropy of mixing as a function of pressure and temperature. This table is to be combined with pure hydrogen and pure helium equations of state to obtain an equation of state that self-consistently includes mixing effects for any hydrogen and helium mixing ratio and may be used to model the interior structure and evolution of giant planets to brown dwarfs. Non-linear mixing must be included in accurate calculations of the equations of state of hydrogen and helium. Uncertainties on the equation of state still exist, however. Ab initio calculations of the behaviour of the hydrogen-helium mixture in the megabar regime for various compositions should be performed in order to gain accuracy.


Introduction
The hydrogen and helium equation of state (hereafter H-He EOS) is crucial in many contexts, but is particularly challenging to model in a high-pressure, low-temperature regime where interactions between molecules, atoms and ions are substantial, near a megabar and at temperatures lower than 10 4 K (e.g., Helled 2018).These conditions concern giant planets (Guillot 2005;Fortney et al. 2010;Militzer et al. 2016) and brown dwarfs in particular (Baraffe 2014).With pressures reaching up to terabar (for the most massive and cold brown dwarfs), these objects span a range of thermodynamical conditions that includes low-pressure regimes for which experimental data are available and high-pressure regimes that are reasonably well understood (Helled et al. 2020), but, with temperatures in the megabar regime ranging from about 5000 K to 20 000 K (Guillot & Gautier 2015), they cross this difficult and uncertain regime in parameter space.Understanding their interior structure and evolution requires an accurate EOS relating physical parameters such as pressure, temperature, density and entropy.
Except for the specific case of the hydrogen-helium phase separation (e.g., Brygoo et al. 2021), most of the effort has been invested so far in the analysis of pure systems of either hydrogen Table B.1 is only available at the CDS via anonymous ftp to cdsarc.cds.unistra.fr(130.79.128.5) or via https:// cdsarc.cds.unistra.fr/viz-bin/cat/J/A+A/672/L1or helium.This has been realised by means of laboratory experiments with anvil cells and laser-driven shocks, and theoretically, with approaches such as density functional theory (DFT), path integral Monte Carlo (PIMC) and quantum Monte Carlo (QMC) based on first-principles simulations (see Helled et al. 2020 for more details).The pure EOSs are usually coupled with the linear mixing rule (also known as the additive volume rule), whose principle is to consider that for extensive variables such as volume and entropy, the value of the mixture is equal to the sum of the pure species weighted by their respective abundances (Saumon et al. 1995).However, Vorberger et al. (2007) pointed out that in the regime in which hydrogen transitions from a molecular to an atomic fluid, the presence of helium shortens and strengthens the molecular hydrogen bonds, leading to nonlinear effects and potentially non-negligible deviations from the linear mixing hypothesis.This implies that in addition to pure EOSs, the behaviour of mixtures should be calculated for all values of the abundances considered.Some DFT simulations have been conducted to study the conditions for a phase separation of hydrogen and helium, but only for a limited range of parameters (Morales et al. 2013;Schöttler & Redmer 2018a,b, hereafter SR18).An extensive DFT calculation for a H-He mixture has been carried out by Militzer & Hubbard (2013) to provide an EOS (the so-called MH13 EOS) that fully accounts for the interactions between hydrogen and helium particles for temperatures between 1000 and 80 000 K and pressures between 0.1 and 300 Mbar.However, it is available only for one specific composition, namely a mixture of 220 hydrogen atoms for 18 helium atoms, corresponding to a helium mass-mixing ratio Y = 0.245.Extensions of this table to cover a wider range of pressures and temperatures have been made available by Miguel et al. (2016) and Chabrier & Debras (2021, hereafter CD21).However, while these works provide EOSs that in principle can be applied to any mixture of hydrogen and helium, the treatment of non-linear effects remains fixed and equal to that calculated by Militzer & Hubbard (2013) at Y = 0.245.
Our work aims to obtain a H-He EOS that includes the hydrogen and helium interactions and remains valid for any hydrogento-helium ratio.This EOS should recover the CD21 EOS as well as both the pure end members CMS19-H and CMS19-He (Chabrier et al. 2019), when there is only hydrogen or only helium in the system, respectively.Section 2 describes how the table is built that contains the contribution of the interactions between hydrogen and helium.Sections 3 and 4 are dedicated to applications to Solar System planets and exoplanets, respectively.

Derivation of the non-ideal mixing effects
In order to account for the interactions between hydrogen and helium, we make use of the linear mixing rule, but add a term ∆V or ∆S that corresponds to an excess or deficit of volume or entropy, respectively, where X and Y are the mass fractions of hydrogen and helium, ρ H , ρ He and ρ H−He are the densities of pure hydrogen, pure helium, and the hydrogen-helium mixture, respectively, and S H , S He , and S H−He are the specific entropies of pure hydrogen, pure helium, and the hydrogen-helium mixture, respectively.The volume of mixing ∆V and entropy of mixing ∆S depend on X and Y, namely the composition of the mixture.All quantities are evaluated at a given pressure P and temperature T .A non-zero volume of mixing results from interactions between hydrogen and helium.The entropy of mixing is the sum of a non-ideal part and an ideal part ∆S ideal , arising even without these interactions.Following Chabrier et al. (2019) (see also Saumon et al. 1995), we write where k B refers to the Boltzmann constant, x i is the number fraction of component i, m H is the atomic mass unit, and A = x i A i , A i being the mass number.
Given the known EOSs for pure hydrogen, for pure helium, and for one mixture of hydrogen and helium with mass mixing ratios X * = 1 − Y * and Y * , respectively, a natural choice is to assume as a first-order approximation, V mix and S mix are then quantities that are independent of X and Y (but vary as a function of P and T ).Using Eqs. ( 1) and (2), we can hence calculate V mix and S mix ,  Eq. ( 6)) with its first-order approximation, assumed proportional to XY, the product of the mass-mixing ratios of hydrogen and helium (orange curve), calculated as ∆S ideal ( Ỹ) XY/( X Ỹ), with Ỹ = 0.245.
Here we have included a complication that arises when the table for the mixture is calculated for a helium mass-mixing ratio Y * , but the mixing terms have been evaluated for a different mixing ratio Ỹ.In the case that we consider, the CD21 H-He table is given for a helium mass-mixing ratio Y * = 0.275, but includes the non-ideal mixing effects from Militzer & Hubbard (2013), who evaluated the H-He interactions for a mixture with Ỹ = 0.245.Equation.( 6) also accounts for the fact that the CD21 EOS evaluates the ideal entropy of mixing at Y * = 0.275, but the non-ideal entropy of mixing at Ỹ = 0.245.We note that Eq. ( 3) yields a value of ∆S ideal that is not strictly proportional to XY, in contrast to the assumption used in Eq. ( 4).However, Fig. 1 shows that the error made by approximating the ideal entropy of mixing as if it were proportional to XY is small: it is always smaller than 0.1 k B /baryon, and when we limit ourselves to compositions such that Y < 0.4, the error made is smaller than 0.01 k B /baryon.This thus supports our assumption that mixing terms should to first order be proportional to XY.
Using Eqs. ( 5) and ( 6), the H-He table from Chabrier & Debras (2021), and the pure H and He tables from Chabrier et al. (2019), which all use the same P − T grid (pressures from 10 −11 to 10 11 Mbar and temperatures from 10 2 to 10 8 K), we can calculate a new table for V mix and S mix for the same set of pressures and temperatures.We stress that this table is valid in the same domain as the EOSs from Chabrier et al. (2019) and Chabrier & Debras (2021), whose limitations concern the regions in which molecular hydrogen and ionised hydrogen become solid as well as the region in which ion quantum effects become important, namely at low T and high P (see Figs. 1 and 16 of Chabrier et al. (2019) for the precise locations in the phase diagrams).The table was cleaned as described in Appendix B to avoid spurious numerical effects and is available for download.Combined with the CMS19 tables, it can be used to calculate the internal structure of giant planets and brown dwarfs of variable compositions and including non-ideal mixing effects.The EOS for any mixture may then be evaluated from the EOSs for the pure elements, Eqs. ( 1) and ( 2), ∆V = XYV mix and ∆S = XYS mix .
L1, page 2 of 7 Figure 2 shows the relative amplitude of the mixing quantities ∆V and ∆S in a mixture with Y = 0.245.A typical Jupiter adiabat passes a region between 0.01 Mbar and 0.04 Mbar, where the non-ideal mixing effects are strongest.In this region, the density of the hydrogen and helium mixture inside Jupiter is denser by almost 10% compared to the density obtained without accounting for the H-He interactions.Around 0.06 Mbar, the sign of ∆V changes and becomes positive (the mixture becomes less dense).Then, ∆V continues to increase to approximately 0.2 Mbar and starts to decrease afterwards.
The variations in ∆V obtained in the 0.1−1 Mbar pressure range for Jupiter conditions are in line with those obtained by Vorberger et al. (2007).However, we point out that we obtain unanticipated negative contributions (i.e.non-linear mixing terms yielding higher densities) at lower pressures.We note that this low-pressure region, around 0.01 Mbar, is a sensitive area in which EOS tables are combined, possibly yielding interpolation errors.
The bottom panel of Fig. 2 shows that the presence of helium yields a mixing term that is small in the classical regime at low pressures and high-enough temperatures, but reaches up to 10% away from it.On the Jupiter adiabat, this increase in ∆S occurs at pressures close to 0.1 Mbar, peaks in the megabar region, and slowly decreases at higher pressures.For temperatures that are about five times lower than on the Jupiter adiabat, we enter the invalidity regime of the EOSs.Results obtained in this regime should not be trusted.
Figure 3 shows that in the classical regime, in which T 10 6 K (P/Mbar), the entropy excess is equal to the ideal value from Eq. (3).At lower temperatures, the mixing entropy term increases to up to almost three times the ideal value, except in a small region centred on 10 −2 Mbar and 2000 K, where we see a decrease to only about 0.5 times the ideal value.Again, this decrease may be due to interpolation issues inherent to the original CMS19 and CD21 tables in a region that combines output from different works.For a typical Jupiter adiabat, ∆S is equal to the classical value up to 5 kbar and then increases to between 1.5 and 2.5 times this value, with a maximum around 0.5 Mbar.It is important to realise that uncertainties on both the pure and mixed EOSs still exist, as shown by a comparison to the SR18 calculations in Appendix A.

Application to Jupiter and Saturn
After deriving a table containing the mixing terms, we applied it to the cases of Jupiter and Saturn.We first compared adiabats for pure hydrogen-helium mixtures just by integrating the adiabatic gradient starting from 1 bar, 166.1 K (i.e.adapted to Jupiter conditions).We used a simple homogeneous model, L1, page 3 of 7 We note that all the new EOSs are much denser than SCvH95 (Saumon et al. 1995) in a wide region, from pressures of about 0.1 to 10 Mbar.The CMS19 (Chabrier et al. 2019) and MLS22 (Mazevet et al. 2022) EOSs do not include the interactions between hydrogen and helium at all.They were obtained using Eqs.( 1) and (2) with the mixing terms ∆V and ∆S both equal to 0 and combining a pure hydrogen table (CMS19-H or MLS22-H) with a pure helium table (CMS19-He or SCvH95-He).They are particularly less dense between 0.01 and 0.1 Mbar but denser at depth than other EOSs.The MGF16+MH13 (Miguel et al. 2016) and CD21 (Chabrier & Debras 2021) EOSs include the non-ideal mixing effects from Militzer & Hubbard (2013), who evaluated the H-He interactions for a mixture with Ỹ = 0.245 and thus are very close to MH13* (Militzer & Hubbard 2013).They are also based on Eqs. ( 1) and ( 2) with ∆V and ∆S both equal to 0; but the non-ideal mixing effects evaluated at Ỹ = 0.245 are included in the pure effective hydrogen table (MGF16+MH13-H or CD21-H; see Chabrier & Debras 2021 for more explanations on how to derive an effective pure H table that accounts for the H-He interactions).While MGF16+MH13-H is combined with SCvH95-He, CD21-H is combined with CMS19-He.Then, based on the table derived in Sect.2, we added the mixing terms to the CMS19 and MLS22 EOSs and obtained adiabats (CMS19 w/NIE and MLS22 w/NIE) that are closer to the other EOSs that include the non-ideal mixing effects (MGF16+MH13, MH13*, and CD21).
An EOS that does not properly include the interactions between hydrogen and helium tends to underestimate the total amount of heavy elements in Jupiter.A quick test that consists of estimating the amount of heavy elements needed to match the equatorial radius of Jupiter, using a static model, can give a rough estimate of the impact of taking the non-ideal mixing effects into account.Focusing on a model similar to the one used for the adiabats of Fig. 4, with Y = 0.245, we ran two models: one model with CMS19, and one model with CMS19, including the non-ideal mixing effects.In the first case, we needed to add a compact core of 14.1 M ⊕ to fit the equatorial radius of Jupiter, while we needed to add a compact core of 20.8 M ⊕ in the second case.This shows that using the CMS19 EOS (see Ni 2019;Nettelmann et al. 2021;Miguel et al. 2022) leads to an underestimation of the amount of heavy elements in Jupiter.
Another issue arises when the helium mixing ratio departs from the value Ỹ = 0.245, and particularly, as in Jupiter and Saturn, when the helium abundance varies from the protosolar value due to H-He phase separation (Stevenson & Salpeter 1977).In this case, even EOSs based on Militzer & Hubbard (2013) such as MGF16+MH13 and CD21 continue to evaluate the non-ideal mixing effects at Ỹ = 0.245.In order to evaluate the magnitude of this effect, we proceeded as follows.First, we calculated some adiabats (at 1 bar conditions for Jupiter) for several values of Y and for the three EOSs CMS19, CD21, and our new EOS (CMS19 w/NIE), which is based on the combination of CMS19-H, CMS19-He, and the mixing terms derived in Sect.2, using Eqs.( 1) and (2).We then searched for the maximum difference in density between the adiabats derived for each EOS at every value of Y. Figure 5 shows that the maximum deviation between our new EOS and CMS19 reaches a peak near Y = 0.5 and is equal to 0 at Y = 0 and Y = 1 because our new EOS is constructed to recover the pure end members CMS19-H and CMS19-He.With CD21, the maximum deviation reaches a minimum near Y = 0.245, which corresponds to the value where the non-ideal mixing effects were estimated.(The minimum is not exactly zero because of the difference between Ỹ = 0.245 and Y * = 0.275.)In the case of Jupiter, a relative difference of 12% can exist (for Y close to 0.3) between our new EOS and CMS19.The comparison with CD21 yields a 2.5% difference in density at most.For Saturn, where the range of Y is wider (Mankovich & Fuller 2021), the effects of the H-He interactions are even stronger.Using CMS19 can lead to a 15% error in density (for Y = 0.5), while using CD21 can lead to an error of 8.5%.

Application to exoplanets
So far, we focused on static models to analyse internal structures at a given time.Now, we use CEPAM and a non-grey atmosphere L1, page 4 of 7 Fig. 6.Evolution models of giant planets for three EOSs: SCvH95, CMS19, and our new EOS (CMS19 w/non-ideal effect), with T eq = 1500 K, taken at 4.5 Gyr.The models consist of a pure rocky core overlaid with a H-He envelope of solar composition.The coloured area for each EOS is bounded by evolution models without a core and with a 20 M ⊕ core.(Guillot et al. 2006;Parmentier et al. 2015) to model the evolution of planets and determine the influence of the H-He interactions during their lifetime.As seen in Sect.3, when we do not include non-ideal mixing effects, static models lead to an underestimated amount of heavy elements, thus a denser mixture of hydrogen and helium overall.Hence, static models (with similar quantities of heavy elements) will yield a larger radius when accounting for these non-ideal mixing effects.Assuming simple structures consisting of a central rocky core and a surrounding H-He envelope of solar composition (Y = 0.3), we calculated evolution models for giant planets with T eq = 1500 K.For each EOS, Fig. 6 shows the mass-radius relation obtained at 4.5 Gyr, bounded by homogeneous models without a core (on top) and with a core of 20 M ⊕ (at the bottom).Static models (for fixed boundary conditions) yield larger radii when non-ideal mixing effects are included, by up to 2% for planets of 1-3 M Jup .When we consider the evolution, the present-day radius results both from the mass-radius relation for static models and from the history of the planet in the P − T diagram.Figure 6 shows that the effects are generally in the opposite sense, leading to a change in radius due to non-ideal mixing that is smaller than for purely static models.We find larger radii for planets that are more massive than Jupiter when we use our new EOS (CMS19 w/NIE) compared to CMS19; and in contrast to static models, we find smaller radii (by up to 2%) for planets that are less massive than Jupiter.The difference can reach 6% at most, at younger ages.Furthermore, we stress that the radii differ significantly (by up to 8%) from those obtained with the SCvH95 EOS, which is commonly used (Fortney et al. 2007).Because of the accuracy in masses and radii of exoplanets, more careful modelling is needed and improved knowledge of the H-He EOS is important to characterise exoplanets.

Conclusion
We provide a simple way to account for the interactions between hydrogen and helium in the EOS.We have derived mixing terms (including the non-ideal contribution) to be added to the clas-sically used linear mixing law for density and entropy, respectively.We built a table of these mixing terms as a function of pressure and temperature.This table is to be combined with pure EOSs to derive the properties of any hydrogen and helium mixture.We emphasise that omitting these mixing effects can lead to a relative error in density of up to 15% on the H-He EOS for Jupiter and Saturn, resulting in a serious underestimation of the heavy element reservoir and invalidating the fit to the observed gravitational moments.Even models including these non-ideal mixing effects evaluated for a helium mass-mixing ratio Y = 0.245 (Debras & Chabrier 2019;Miguel et al. 2022;Militzer et al. 2022;Mankovich & Fuller 2021) have intrinsic errors on the density evaluated to up to 2.5% for Jupiter and 8.5% for Saturn.In addition, when the non-ideal mixing effects are not included, a change of up to 6% in the calculated radii of exoplanets may result.Furthermore, given the differences between EOSs, we stress that uncertainties on the H-He EOS still exist.
Our approach provides a way to obtain a self-consistent EOS for different helium abundances.We stress that it relies on EOSs calculated for three compositions only (for Y = 0, Y = 0.245 and Y = 1).Ab initio calculations and high-pressure experiments for P = 10 −2 −10 Mbar and T = 1000−20 000 K with other hydrogen-helium compositions, but also including heavy elements, would be highly desirable.1).For T = 10000 K, r s = 1.4 a 0 and 1.25 a 0 , and various helium fractions x He , we calculated the total entropy (using Eq. ( 2), with CMS19-H, CMS19-He, and our table of the mixing terms to compute S H , S He and S mix , respectively).We also show the entropy resulting from the ideal mixing of only the pure terms (XS H + YS He ).The difference between the two curves represents the entropy of mixing ∆S .First, we note some important differences on the pure H and He EOSs between CMS19 and the work of SR18.Typically, for r s = 1.4 a 0 , the entropy of CMS19-H is 0.2 k b /atom lower than the pure H EOS from SR18, while the entropy of CMS19-He is 0.53 k b /atom higher than the pure He EOS from SR18.In H-He mixtures, the entropy we calculate (where the mixing entropy is initially based on the MH13 EOS) is always higher than the entropy from SR18 (0.58 k b /atom greater at most).  1 of their supplemental material) for T = 10000 K and various helium fractions x He .The two panels show results for values of Wigner-Seitz radii of r s = 1.4 a 0 and r s = 1.25 a 0 respectively.Each panel displays our calculations (red) of the total entropy (computed using Eq. 2 with CMS19-H, CMS19-He and our table for ∆S ) and also the entropy without including the entropy of mixing ∆S (hence corresponding to XS H + YS He ).Results from SR18 are shown in black.

Appendix A: Comparisons of the entropy of mixing
These intrinsic differences in the pure EOSs and also in the entropies yielded for H-He mixtures lead to discrepancies between the entropy of mixing that we obtain and the one from SR18.In this parameter space (T = 10000 K, r s = 1.4 a 0 or 1.25 a 0 ), while the entropy of mixing of SR18 ranges from 0.23 to 0.66 k b /atom, our entropy of mixing ranges from 0.26 to 1.2 k b /atom.The entropy of mixing that we obtain is hence 1.1 to 2.2 times higher than the entropy of mixing of SR18.(Similar results would be obtained using the Morales et al. (2009) calculations because their entropies of mixing are similar to those of SR18).This highlights the fact that the EOSs for pure hydrogen, pure helium, and mixtures of hydrogen and helium obtained by the different groups still differ; this is particu-larly obvious for the pure helium EOS used here (CMS19-He) and the one obtained by SR18.While CMS19 have used the Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional in their simulations, SR18 applied the van der Waals density functional (vdW-DF) exchange-correlation functional.Both exchange-correlation functionals support different experiments and have their own benefits and limitations (see SR18 for a more detailed discussion).A consistent, homogeneous exploration of EOSs in this regime of pressures and temperature and for different compositions is crucial.
In spite of these uncertainties, we believe that using the table of non-ideal mixing effects that we provide in addition to the CMS19 EOS is an improvement for two reasons: (i) It includes non-ideal mixing effects that are otherwise ignored, and for low helium abundances (Y <0.3), it yields total entropies that differ from those of SR18 by less than 4%.(ii) It enables calculations in the high-temperature regions spanned by massive planets and brown dwarfs, which are otherwise not possible with the CD21 EOS.

Fig. 1 .
Fig. 1.Comparison of the ideal entropy of mixing (blue curve and

Fig. 2 .
Fig. 2. Relative contributions of the mixing terms ∆V and ∆S to the total specific volume 1/ρ H−He and entropy S H−He , respectively, as a function of pressure and temperature.The solid black line corresponds to the Jupiter adiabat from Miguel et al. (2022).The values of ∆V × ρ H−He and ∆S /S H−He are calculated for X * = 0.725 and Y * = 0.275.The hashed areas correspond to the region in which the mixing volume and entropy are invalid (see Chabrier et al. 2019, for more details).

Fig. 3 .
Fig. 3. Ratio of the mixing entropy ∆S and the ideal entropy of mixing ∆S ideal for a mixture of Y = 0.245.The hashed area corresponds to the region in which the entropy of mixing is invalid (see Chabrier et al. 2019, for more details).

Fig. 4 .
Fig. 4. Adiabats obtained from different EOSs and corresponding to a homogeneous model, without a compact core, with Y = 0.245.NIE stands for non-ideal effects.

Fig. 5 .
Fig. 5. Comparison in density of our new EOS (this work) based on the non-ideal mixing effects derived in Sect. 2 with both CD21 and CMS19, according to the helium mass fraction Y.The shaded red and grey areas correspond to the range of values of Y in Jupiter (taken from models of Miguel et al. 2022) and Saturn (taken from Mankovich & Fuller 2021), respectively.

Figure
Figure A.1 shows a comparison of the entropies obtained from our work and by SR18 (SM Fig.1).For T = 10000 K, r s = 1.4 a 0 and 1.25 a 0 , and various helium fractions x He , we calculated the total entropy (using Eq. (2), with CMS19-H, CMS19-He, and our table of the mixing terms to compute S H , S He and S mix , respectively).We also show the entropy resulting from the ideal mixing of only the pure terms (XS H + YS He ).The difference between the two curves represents the entropy of mixing ∆S .First, we note some important differences on the pure H and He EOSs between CMS19 and the work of SR18.Typically, for r s = 1.4 a 0 , the entropy of CMS19-H is 0.2 k b /atom lower than the pure H EOS from SR18, while the entropy of CMS19-He is 0.53 k b /atom higher than the pure He EOS from SR18.In H-He mixtures, the entropy we calculate (where the mixing entropy is initially based on the MH13 EOS) is always higher than the entropy from SR18 (0.58 k b /atom greater at most).
Fig. A.1.Comparison of our entropy with the one calculated by SR18 (from Fig.1of their supplemental material) for T = 10000 K and various helium fractions x He .The two panels show results for values of Wigner-Seitz radii of r s = 1.4 a 0 and r s = 1.25 a 0 respectively.Each panel displays our calculations (red) of the total entropy (computed using Eq. 2 with CMS19-H, CMS19-He and our table for ∆S ) and also the entropy without including the entropy of mixing ∆S (hence corresponding to XS H + YS He ).Results from SR18 are shown in black.