---
title: CDM3Y-IVF1 Equations of State
url: https://www.emergentmind.com/topics/cdm3y-ivf1-equations-of-state
type: topic
---

# CDM3Y-IVF1 Equations of State

CDM3Y-IVF1 equations of state are a family of semi-microscopic nuclear equations of state derived from the M3Y-Paris nucleon-nucleon interaction within a non-relativistic Hartree-Fock framework and extended to cold, charge-neutral, $\beta$-equilibrated $npe\mu$ matter relevant to neutron-star cores. Their defining feature is a density-dependent isoscalar sector calibrated to symmetric nuclear matter saturation and a separately tuned isovector sector that produces a stiff symmetry energy at supra-saturation density. In recent neutron-star applications, this framework has been used to span very soft to extremely stiff nuclear matter through variations of the saturation incompressibility $K_0$, while simultaneously retaining a full hierarchy of higher-order symmetry-energy coefficients. That structure makes CDM3Y-IVF1 especially suitable for systematic studies of core-crust transition properties, crustal moment of inertia, direct Urca thresholds, adiabatic indices, compactness, and tidal observables [2507.18384].

## 1. Microscopic definition and interaction structure

The CDM3Y construction starts from the semi-realistic M3Y-Paris interaction, decomposed into central isoscalar and isovector components with direct and exchange terms. Medium effects are incorporated through explicit density-dependent factors multiplying the finite-range interaction, so that the effective interaction is written in isoscalar and isovector channels as
\[
v^{D,Ex}_{00(01)}(\rho,r)=F_{0(1)}(\rho)\,v^{D,Ex}_{00(01)}(r),
\]
with
\[
F_{0(1)}(\rho)=C_{0(1)}\left(1+\alpha_{0(1)}e^{-\beta_{0(1)}\rho}-\gamma_{0(1)}\rho\right).
\]
A given set of parameters $\{C_i,\alpha_i,\beta_i,\gamma_i\}$ defines one member of the CDM3Y-$K_0$ family [2507.18384].

In asymmetric nuclear matter with baryon density $\rho$ and proton fraction $x_p=\rho_p/\rho$, the Hartree-Fock energy per nucleon is
\[
\begin{aligned}
E_A(\rho,x_p) &=  \frac{3\hbar^2 k_F^2}{20m} \Big[(2-2x_p)^{5/3} + (2x_p)^{5/3}\Big] \\
&\quad + \frac{\rho}{2} \Big[ F_0(\rho) J^{D}_{00} + (1-2x_p)^2 F_1(\rho) J^{D}_{01} \Big] \\
&\quad + \frac{\rho}{8} \int d\vec r \Big[ F_0(\rho) v^{Ex}_{00}(r) B_0^2(r,x_p) + F_1(\rho) v^{Ex}_{01}(r) B_1^2(r,x_p) \Big],
\end{aligned}
\]
where $k_F=(3\pi^2\rho/2)^{1/3}$, $k_{Fn}=(3\pi^2\rho_n)^{1/3}$, $k_{Fp}=(3\pi^2\rho_p)^{1/3}$, and
\[
\begin{aligned}
B_0(r,x_p) &= 2\Big[(1-x_p)\hat J_1(k_{Fn}r)+x_p\hat J_1(k_{Fp}r)\Big],\\
B_1(r,x_p) &= 2\Big[(1-x_p)\hat J_1(k_{Fn}r)-x_p\hat J_1(k_{Fp}r)\Big],
\end{aligned}
\]
with $\hat J_1(x)=3j_1(x)/x$ [2507.18384].

Within the IVF1 parametrization, the isoscalar density dependence $F_0(\rho)$ is adjusted to reproduce symmetric nuclear matter saturation in Hartree-Fock, with $\rho_0\approx 0.16\,\mathrm{fm}^{-3}$, $E_A(\rho_0,0)\approx -16\,\mathrm{MeV}$, and a prescribed incompressibility
\[
K_0=150\text{--}330~\mathrm{MeV}.
\]
The isovector density dependence $F_1(\rho)$ is fitted independently by matching Brueckner-Hartree-Fock neutron optical potentials, which yields a stiff symmetry energy at supra-saturation density; for Paris-IVF1, the quoted saturation values are $E_{\text{sym}0}\simeq 29.10~\mathrm{MeV}$ and $L\simeq 48.2~\mathrm{MeV}$ [2507.18384].

The broader CDM3Y literature clarifies the structural meaning of this construction. The underlying finite-range M3Y radial interaction is retained, while the density dependence is applied separately in the isoscalar and isovector channels. In earlier CDM3Y-n studies, the stiff isovector sector was likewise tied to Brueckner-Hartree-Fock or JLM isovector optical-potential information, whereas soft variants were generated by taking the isovector density dependence proportional to the isoscalar one [1105.5222]. This establishes CDM3Y-IVF1 as part of the stiff CDM3Y branch rather than the soft CDM3Y-ns branch.

## 2. Symmetry-energy expansion and higher-order coefficients

A central feature of CDM3Y-IVF1 is that it retains not only the quadratic symmetry energy but also higher even powers in the isospin asymmetry
\[
\mathcal{I}=\frac{\rho_n-\rho_p}{\rho}=1-2x_p.
\]
The nucleonic energy per particle is expanded as
\[
E_A(\rho,\mathcal{I})=E_A(\rho,0)+E_{\text{sym}2}(\rho)\mathcal{I}^2+E_{\text{sym}4}(\rho)\mathcal{I}^4+E_{\text{sym}6}(\rho)\mathcal{I}^6+\cdots,
\]
with
\[
E_{\text{sym}n}(\rho)=\frac{1}{n!}\left.\frac{\partial^n E_A(\rho,\mathcal{I})}{\partial \mathcal{I}^n}\right|_{\mathcal{I}=0},
\qquad n=0,2,4,\ldots
\]
Here $E_{\text{sym}2}$ is the usual symmetry energy, while $E_{\text{sym}4}$ and $E_{\text{sym}6}$ encode higher-order isovector structure [2507.18384].

Around saturation density, the expansion is organized through the scaled variable
\[
\mathbb{X}=\frac{\rho-\rho_0}{3\rho_0}.
\]
The symmetric-matter sector is written as
\[
\begin{aligned}
E_A(\rho,0)&=E_A(\rho_0,0)+\frac{K_0}{2!}\mathbb{X}^2+\frac{Q_0}{3!}\mathbb{X}^3+\frac{I_0}{4!}\mathbb{X}^4\\
&\quad+\frac{H_0}{5!}\mathbb{X}^5+\frac{G_0}{6!}\mathbb{X}^6+\cdots,
\end{aligned}
\]
while the quadratic and quartic symmetry-energy contributions are expanded as
\[
\begin{aligned}
E_{\text{sym}2}(\rho)&=E_{\text{sym}}(\rho_0)+L\,\mathbb{X}+\frac{K_2}{2!}\mathbb{X}^2+\frac{Q_2}{3!}\mathbb{X}^3\\
&\quad+\frac{I_2}{4!}\mathbb{X}^4+\frac{H_2}{5!}\mathbb{X}^5+\frac{G_2}{6!}\mathbb{X}^6+\cdots,
\end{aligned}
\]
\[
\begin{aligned}
E_{\text{sym}4}(\rho)&=E_{\text{sym}4}(\rho_0)+L_4\,\mathbb{X}+\frac{K_4}{2!}\mathbb{X}^2+\frac{Q_4}{3!}\mathbb{X}^3\\
&\quad+\frac{I_4}{4!}\mathbb{X}^4+\frac{H_4}{5!}\mathbb{X}^5+\frac{G_4}{6!}\mathbb{X}^6+\cdots.
\end{aligned}
\]
The coefficients $L_n$, $K_n$, $Q_n$, $I_n$, $H_n$, and $G_n$ are, respectively, slope, curvature, skewness, kurtosis, fifth-order, and sixth-order density derivatives at $\rho_0$ [2507.18384].

The incompressibility of asymmetric matter at its own saturation density is written as
\[
K_{0\mathcal{I}}=K_0+K_{\tau2}\mathcal{I}^2+K_{\tau4}\mathcal{I}^4+K_{\tau6}\mathcal{I}^6+\cdots,
\]
where $K_{\tau2}$, $K_{\tau4}$, and $K_{\tau6}$ are isobaric incompressibility coefficients [2507.18384].

The recent higher-order analysis identifies two opposite classes of coefficients. EOS stiffening, in the sense of increasing high-density pressure, is correlated with increasing $K_0$ and $Q_4$, and with less negative values of $Q_0$, $H_0$, $K_4$, $K_{\tau4}$, $I_{2,4}$, and $G_2$. By contrast, increasing $I_0$, $G_0$, $Q_2$, $H_2$, $K_6$, and $K_{\tau6}$, or making $K_2$ and $K_{\tau2}$ less negative, produces systematically opposite neutron-star trends relative to increasing $K_0$ [2507.18384]. This hierarchy is one of the main reasons CDM3Y-IVF1 is not reducible to a single-parameter incompressibility model.

Earlier work on CDM3Y-type neutron-star matter often relied on the parabolic approximation
\[
E(n_{\rm b},\delta)=E(n_{\rm b},0)+S(n_{\rm b})\delta^2+O(\delta^4),
\]
and treated higher-order terms as small in many applications [1105.5222]. The more recent CDM3Y-IVF1 treatment shows that explicit fourth- and sixth-order sectors can systematically shift observable thresholds and crustal properties [2507.18384]. A common misconception is therefore that only $J$, $L$, and $K_0$ control neutron-star observables in this family; the recent results show that high-order isoscalar, isovector, and isobaric coefficients modify several trends in a non-negligible way.

## 3. Extension to cold $\beta$-equilibrated $npe\mu$ matter

For neutron-star cores, CDM3Y-IVF1 is extended to uniform cold matter composed of neutrons, protons, electrons, and muons. The total energy density is
\[
\varepsilon(\rho,x_p,x_e,x_\mu)=\varepsilon_b(\rho,x_p)+\sum_{\ell=e,\mu}\varepsilon_\ell(\rho,x_\ell),
\]
with
\[
\varepsilon_b(\rho,x_p)=\rho\Big[E_A(\rho,x_p)+x_p m_p c^2+(1-x_p)m_n c^2\Big].
\]
Electrons and muons are treated as free relativistic Fermi gases with
\[
\mu_\ell=\sqrt{\hbar^2 c^2 k_{F\ell}^2+m_\ell^2 c^4},
\qquad
k_{F\ell}=(3\pi^2\rho_\ell)^{1/3},
\]
and the total pressure is
\[
P(\rho,x_p,x_e,x_\mu)=\rho^2\frac{\partial E_A(\rho,x_p)}{\partial \rho}+\sum_{\ell=e,\mu}(\mu_\ell\rho x_\ell-\varepsilon_\ell).
\]
The equilibrium conditions at $T=0$ are
\[
\mu_n-\mu_p=\mu_e=\mu_\mu,
\qquad
\mu_n-\mu_p=-\frac{\partial E_A(\rho,x_p)}{\partial x_p},
\]
together with charge neutrality
\[
x_p=x_e+x_\mu.
\]
These relations lead to an implicit equation for the proton fraction $x_p(\rho)$ [2507.18384].

In this framework, the isovector sector directly controls the composition of the stellar core. For CDM3Y-Paris-IVF1, the stiff symmetry energy produces higher proton content at intermediate densities than IVF0, but at very high densities the proton fraction decreases with increasing $K_0$ [2507.18384]. Quantitatively, for a maximum-mass star with soft EOS $K_0=160$ MeV, the central proton fraction reaches $x_{pc}\approx 0.376$ at $M_{\max}\simeq 0.91\,M_\odot$, whereas for a very stiff EOS $K_0=330$ MeV, $x_{pc}\approx 0.077$ at $M_{\max}\simeq 2.4\,M_\odot$ [2507.18384].

For a fixed $K_0$, the central proton fraction increases with stellar mass; one quoted example rises from $\approx 0.067$ at $0.7\,M_\odot$ to $\approx 0.161$ at its own $M_{\max}$ [2507.18384]. Higher-order coefficients modify these trends: more negative $K_{\tau2}$ and more positive $K_{\tau4}$ decrease $x_{pc}$, whereas larger $I_0$, $G_0$, $Q_2$, $K_6$, $H_2$, $K_{\tau6}$ and less negative $K_2$ increase the central proton fraction, opposite to increasing $K_0$ [2507.18384].

The contrast with soft CDM3Y variants is important. In the older stiff-versus-soft classification, stiff CDM3Y-n interactions produce a symmetry energy that rises with density and correspondingly large proton fractions, while soft CDM3Y-ns variants can drive the proton fraction down to a few percent and eventually toward pure neutron matter at high density [1105.5222]. CDM3Y-IVF1 belongs to the stiff class in this sense, but the higher-order analysis shows that even within the stiff class the composition is not monotonic in all stiffness indicators.

## 4. Core-crust transition and crustal structure

The core-crust transition in CDM3Y-IVF1 is determined by the thermodynamic instability of uniform $\beta$-equilibrated matter. The relevant quantity is the generalized compressibility at fixed lepton chemical potential,
\[
K_{\mu_i}=\left(\frac{\partial P}{\partial \rho}\right)_{\mu_i},
\]
which can be expressed as
\[
K_{\mu_i}=\frac{K(\rho,x_p)}{9}-\frac{\left(\rho\,\frac{\partial^2 E_A}{\partial \rho\,\partial x_p}\right)^2}{\frac{\partial^2 E_A}{\partial x_p^2}}.
\]
Uniform matter is stable for $K_{\mu_i}>0$. The transition density $\rho_t$ and transition pressure $P_t$ are defined where $K_{\mu_i}$ first vanishes upon decreasing density [2507.18384].

Within the CDM3Y-Paris and CDM3Y-Reid families, $P_t$ increases with increasing symmetry energy at saturation $E_{\text{sym}0}$. For nearly equal $E_{\text{sym}0}$, $P_t$ decreases as $L$ increases, while $\rho_t$ and the transition proton fraction $x_{pt}$ vary only weakly with $E_{\text{sym}0}$ and $L$, at the level of $\sim 0.003\text{--}0.007$ in absolute fractions [2507.18384]. Along the CDM3Y-Paris-IVF1 sequence, increasing $K_0$ raises both $\rho_t$ and $P_t$, and also increases $x_{pt}$ more mildly [2507.18384].

The higher-order coefficients introduce a more differentiated pattern. Increasing $K_0$ or $Q_4$, and making $Q_0$, $H_0$, $K_4$, $K_{\tau4}$, $I_{2,4}$, and $G_2$ less negative, raises $\rho_t$ and $P_t$ and slightly increases $x_{pt}$. Conversely, increasing $I_0$, $G_0$, $Q_2$, $H_2$, $K_6$, $K_{\tau6}$, or making $K_2$ and $K_{\tau2}$ less negative, modifies $\rho_t$ and $P_t$ in the direction opposite to the change in $K_0$ [2507.18384]. This is one of the clearest demonstrations that high-order coefficients cannot be collapsed into a single effective stiffness label.

A frequent oversimplification is that a stiffer EOS always implies a thicker crust. The CDM3Y-IVF1 results are more specific. A higher transition pressure is associated with a thicker crust and larger radius for a given mass, but a larger transition density implies that the uniform core extends farther outward, which can reduce fractional crust thickness in very massive configurations [2507.18384]. This distinction becomes central once rotational observables are considered.

Earlier CDM3Y-based neutron-star studies combined the uniform core EOS with crust models such as the compressible liquid drop model or the Douchin-Haensel inner crust, mainly to assess global stellar properties and the impact of soft versus stiff symmetry energy [1105.5222; 2507.11379]. The recent CDM3Y-IVF1 analysis brings the transition itself into the foreground by connecting it directly to the high-order symmetry-energy hierarchy [2507.18384].

## 5. Moment of inertia, crustal fraction, and compactness

Static stellar configurations are obtained from the Tolman-Oppenheimer-Volkoff equations,
\[
\frac{dP}{dr}=-\frac{G\varepsilon(r)M(r)}{c^2r^2}\left(1+\frac{P}{\varepsilon}\right)\left(1+\frac{4\pi r^3P}{c^2M(r)}\right)\left(1-\frac{2GM(r)}{c^2r}\right)^{-1},
\]
\[
\frac{dM}{dr}=\frac{4\pi r^2\varepsilon(r)}{c^2}.
\]
For slowly rotating stars, the total moment of inertia $I$ is calculated in the Hartle-Thorne formalism through the frame-dragging function $\omega(r)$ or, equivalently, through the accumulated moment-of-inertia function $\tilde\kappa(r)=J(r)/\Omega$ [2507.18384].

The crust is defined as the region between the core radius $R_{\text{core}}$ at density $\rho_t$ and the stellar radius $R$. Its contribution to the moment of inertia is
\[
I_{\text{crust}}=\tilde\kappa_{\text{cr}}(R),
\qquad
\frac{\Delta I}{I}=\frac{I_{\text{crust}}}{I}.
\]
The same formalism yields the fractional crust thickness $\Delta R/R$ with $\Delta R=R-R_{\text{core}}$ [2507.18384].

For the CDM3Y-Paris-IVF1 family, the total moment of inertia $I$ increases with increasing $K_0$ at fixed stellar mass, and this sensitivity becomes stronger with increasing mass. The maximum mass $M_{\max}(K_0)$ also increases with stiffness; soft EOS with $K_0\lesssim 160$ MeV may fail to support $1\,M_\odot$ [2507.18384]. The crust thickness $\Delta R$ and its fractional value $\Delta R/R$ are strongly anti-correlated with mass, and their dependence on $K_0$ at fixed mass is weak, weakening further for very stiff EOS [2507.18384].

The behavior of the crustal moment-of-inertia fraction is more nuanced than the behavior of the total moment of inertia. For fixed mass, increasing $K_0$ generally increases $I$ and can also increase $\Delta I/I$. For the maximum-mass configuration supported by a given EOS, however, increasing $K_0$ yields a larger $M_{\max}$ but a very thin crust, causing both $\Delta R/R$ and $\Delta I/I$ to decrease [2507.18384]. This distinction directly addresses a common misconception: within CDM3Y-IVF1, EOS stiffening does not have a single universal effect on the crustal fraction; the sign of the trend depends on whether the comparison is made at fixed mass or along the maximum-mass sequence.

The glitch-motivated criterion
\[
\Delta I/I \gtrsim 0.014
\]
is satisfied for neutron-star masses up to $\sim 1.6\,M_\odot$ in this family [2507.18384]. The same work states that both total $I$ and crustal fraction $\Delta I/I$ show approximately linear decrease with compactness
\[
\mathpzc{C}=\frac{GM}{Rc^2},
\]
and that the glitch condition implies $\mathpzc{C}\lesssim 0.21$ and $M\lesssim 1.6\,M_\odot$ within the CDM3Y-Paris-IVF1 family [2507.18384].

The higher-order coefficient dependence again separates into two groups. Increasing $K_0$ or $Q_4$, and reducing the magnitude of negative $Q_0$, $H_0$, $K_4$, $K_{\tau4}$, $I_{2,4}$, and $G_2$, increases $I$ and $\Delta I/I$ for a given mass, but decreases $\Delta I/I$ and $\Delta R/R$ for the maximum-mass configuration. Increasing $I_0$, $G_0$, $Q_2$, $K_6$, $K_{\tau6}$, $H_2$, or making $K_2$ and $K_{\tau2}$ less negative, reduces $I$ and $\Delta I/I$ for fixed mass, yet increases $\Delta R/R$ and $\Delta I/I$ for the maximum-mass configuration [2507.18384]. The paper explicitly notes that these patterns are consistent with independent constraints from Skyrme-based meta-modeling and relativistic mean-field studies.

The broader semi-microscopic CDM3Y program gives compatible global scales. In an earlier stiff-versus-soft survey, the stiff CDM3Y3, CDM3Y4, and CDM3Y6 interactions yielded maximum masses of $1.61\,M_\odot$, $1.73\,M_\odot$, and $1.97\,M_\odot$, respectively, with radii near $10$ km and moments of inertia increasing with stiffness, whereas the soft CDM3Y3s, CDM3Y4s, and CDM3Y6s gave smaller maximum masses and systematically lower moments of inertia [1105.5222]. This older dichotomy provides the global background against which the finer IVF1 higher-order trends should be read.

## 6. Direct Urca thresholds, adiabatic response, and astrophysical constraints

The nucleonic direct Urca process,
\[
n\to p+e^-+\bar\nu_e,\qquad p+e^-\to n+\nu_e,
\]
requires the proton fraction to exceed the threshold
\[
x_{\text{DU}}=\frac{1}{1+(1+r_e^{1/3})^3},
\qquad
r_e=\frac{1}{1+x_\mu/x_e}.
\]
When muons are absent, $x_{\text{DU}}=1/9\simeq 0.111$; when muons are present, $x_{\text{DU}}>1/9$ [2507.18384].

For CDM3Y-Paris-IVF1 over $K_0=150\text{--}280$ MeV, the threshold proton fraction grows weakly from $\approx 0.138$ to $\approx 0.141$, the threshold density increases from $\approx 3.28\rho_0$ to $\approx 6.15\rho_0$, and the threshold pressure rises strongly from $\sim 37$ to $\sim 705$ MeV fm$^{-3}$. For still stiffer EOS with $K_0\gtrsim 280$ MeV, the proton fraction may never reach the direct-Urca threshold before the central density is exceeded, so direct Urca does not occur [2507.18384]. IVF1 yields lower $x_{\text{DU}}$, $\rho_{\text{DU}}$, and $P_{\text{DU}}$ than IVF0, while CDM3Y-Reid-IVF1 behaves similarly but typically with slightly higher DU thresholds than Paris-IVF1 [2507.18384].

The dependence on the higher-order hierarchy is parallel to the crustal analysis. Increasing $K_0$ or $Q_4$, and making $Q_0$, $H_0$, $K_4$, $K_{\tau4}$, $I_{2,4}$, and $G_2$ less negative, raises $x_{\text{DU}}$, $\rho_{\text{DU}}$, and $P_{\text{DU}}$. Increasing $I_0$, $G_0$, $Q_2$, $K_6$, $K_{\tau6}$, $H_2$, or making $K_2$ and $K_{\tau2}$ less negative, lowers these thresholds [2507.18384]. A plausible implication is that modest changes in high-order symmetry coefficients can move DU cooling from the domain of only very massive stars into densities realized by lighter objects.

The local stiffness of cold matter is characterized by the adiabatic index
\[
\Gamma=\frac{\varepsilon+P}{P}\frac{dP}{d\varepsilon}
=\left(\frac{\varepsilon+P}{P}\right)\left(\frac{v_s}{c}\right)^2.
\]
For non-rotating stars with CDM3Y-Paris-IVF1, the central adiabatic index $\Gamma_c$ decreases with increasing stellar mass. For $1\,M_\odot$ and $1.4\,M_\odot$ stars, $\Gamma_c(K_0)$ exhibits a minimum at intermediate $K_0$ and then increases with further stiffening. Along the maximum-mass sequence, $\Gamma_c$ decreases with $K_0$ up to about $240$ MeV and increases beyond this value [2507.18384]. Less negative $K_{\tau4}$ increases $\Gamma_c$, whereas less negative $K_{\tau2}$ or larger $K_{\tau6}$ decreases it, mirroring the two opposite higher-order coefficient classes identified earlier [2507.18384].

Recent CDM3Y studies have also connected the same semi-microscopic EOS family to tidal observables. For CDM3Y parameterizations with $K_0=230\text{--}330$ MeV, the predicted canonical radii are $R_{1.4}=11.97\pm0.54$ km for Paris and $R_{1.4}=12.03\pm0.78$ km for Reid, while the corresponding maximum masses are $M_{\max}\approx 2.07\pm0.34\,M_\odot$ and $2.13\pm0.25\,M_\odot$ [2507.11379]. The same work reports $\Lambda_{1.4}\approx 32$ to $413$ for Paris and $\approx 98$ to $428$ for Reid when $K_0$ varies from 200 to 330 MeV, and states that the $K_0\ge 230$ sets reproduce most NICER and GW-inferred mass-radius constraints [2507.11379]. It also gives an approximately EOS-insensitive relation for $M\ge 1\,M_\odot$,
\[
\ln \Lambda = 11.327 - 32.8\,\mathpzc{C},
\]
with coefficient of determination $R^2=0.992$ across the CDM3Y-Paris and Reid sets with $K_0=230\text{--}330$ MeV [2507.11379].

Taken together, these results place CDM3Y-IVF1-type EOSs in a constrained but nontrivial region of parameter space. They must be stiff enough to sustain heavy pulsars, yet their higher-order symmetry structure controls whether the same EOS produces a sufficiently large crustal moment-of-inertia fraction, permits or suppresses direct Urca cooling, and shifts tidal deformabilities within observationally allowed bands. The recurring theme is that the macroscopic neutron-star phenomenology of CDM3Y-IVF1 is governed not only by $K_0$ and $L$, but by the full coupled hierarchy of isoscalar, isovector, and isobaric coefficients carried by the Hartree-Fock interaction [2507.18384].

Source: https://www.emergentmind.com/topics/cdm3y-ivf1-equations-of-state