pith. sign in

arxiv: 2405.02157 · v1 · submitted 2024-05-03 · ⚛️ nucl-th

Dynamics of dilute nuclear matter with light clusters and in-medium effects

Rui Wang , Stefano Burrello , Maria Colonna , Francesco Matera This is my paper

Pith reviewed 2026-05-24 01:12 UTC · model grok-4.3

classification ⚛️ nucl-th
keywords dilute nuclear matterlight clustersspinodal instabilitiesin-medium effectsMott effectslinear responseheavy-ion collisionsfragment formation
0
0 comments X

The pith

Spinodal instabilities and growth rates in dilute nuclear matter change sharply when light clusters and in-medium Mott effects are included.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper studies the dynamics of dilute nuclear systems made of nucleons plus light clusters by means of a linear response analysis. In-medium effects on the clusters are introduced through a density-dependent momentum cut-off that represents Mott dissolution. The central result is that both the location of spinodal instabilities and the rates at which they grow become strongly altered by the clusters themselves and by how the in-medium modifications are handled. A reader would care because these instabilities control how nuclear matter fragments in low-density conditions. The work therefore supplies a more realistic description of the early stages of breakup in heavy-ion collisions and in astrophysical dilute environments.

Core claim

Within a linear response treatment of dilute matter containing nucleons and light clusters, the inclusion of in-medium Mott effects through a density-dependent momentum cut-off demonstrates that spinodal instabilities and the associated growth rates are severely affected by the presence of the clusters and by the particular treatment of those in-medium effects.

What carries the argument

Linear response approach with a density-dependent momentum cut-off that encodes the in-medium Mott dissolution of light clusters.

If this is right

  • Spinodal instabilities shift when light clusters are present.
  • Instability growth rates depend on the chosen treatment of in-medium effects.
  • Fragment formation in heavy-ion collisions is expected to be influenced by these changes.
  • The results carry over to astrophysical settings that involve dilute nuclear matter.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Models of heavy-ion collisions that omit clusters may therefore predict different fragmentation timescales and sizes.
  • The same mechanism could modify the equation of state used in supernova or neutron-star crust simulations.
  • The density-dependent cut-off would benefit from direct checks against microscopic calculations that include additional many-body corrections.

Load-bearing premise

The linear response framework remains valid once light clusters and the density-dependent momentum cut-off are added to the dilute regime.

What would settle it

A direct comparison, in the same linear response calculation, of spinodal growth rates computed with and without the light clusters plus the density-dependent cut-off that shows no measurable difference would falsify the claim.

Figures

Figures reproduced from arXiv: 2405.02157 by Francesco Matera, Maria Colonna, Rui Wang, Stefano Burrello.

Figure 1
Figure 1. Figure 1: Spinodal border in the (ρb, T ) plane in three cases: 1) pure nucleonic matter (SNM, black); 2) for nuclear mat￾ter with deuterons, including in-medium effects (Rd ≫ Rδρb , red); 3) for nuclear matter with deuterons, neglecting (Rd ≪ Rδρb , cyan) in-medium effects in the dynamics. The result from a “hybrid” case is also included (green, see text for de￾tails). The inset shows the ρb dependence of the deute… view at source ↗
Figure 3
Figure 3. Figure 3: Growth rate of the instability, Im(ω), as a function of the wave number k, for the same cases as in Figs. 1 and 2, at various density and temperature values. The dependence of the growth rate on the wave num￾ber, for two typical density values lying inside the spin￾odal region (ρb = 0.02 fm−3 and ρb = 0.06 fm−3 ) is displayed in [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 2
Figure 2. Figure 2: Growth rate of the imaginary sound velocity, [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 4
Figure 4. Figure 4: The quantity ∆ (see text) as a function of the total baryon density ρb for nuclear matter with deuterons, neglecting (cyan) or including (red) in-medium effects in the dynamics, for three temperature values. Lines are drawn only inside the spinodal region. An in-depth insight into the direction of the unstable modes in the space of density fluctuations is provided by the (δρS/δρd) ratio, where ρS = ρn + ρp… view at source ↗
Figure 5
Figure 5. Figure 5: Mott momentum of deuteron Λd obtained by Eq. (9) as a function of the cubic root of the total baryon density, (ρb/ρ0) 1/3 with ρ0 = 0.16 fm−3 , for three temperature values [PITH_FULL_IMAGE:figures/full_fig_p009_5.png] view at source ↗
read the original abstract

We investigate the dynamics of dilute systems composed of nucleons and light clusters within a linear response approach, taking into account the in-medium Mott effects on cluster appearance, through a density-dependent momentum cut-off. We find that spinodal instabilities and associated growth rates are severely affected by the presence of light clusters and, in particular, by the treatment of in-medium effects, foreshadowing intriguing consequences for fragment formation in heavy-ion collisions and in the broader astrophysical context.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 2 minor

Summary. The manuscript investigates the dynamics of dilute nuclear matter consisting of nucleons and light clusters using a linear response approach. In-medium Mott effects on cluster appearance are incorporated via a density-dependent momentum cut-off. The central claim is that spinodal instabilities and their associated growth rates are severely affected by the presence of light clusters and, in particular, by the treatment of in-medium effects, with implications for fragment formation in heavy-ion collisions and astrophysical contexts.

Significance. If the central claim holds after addressing the validity of the framework, the result would underscore the importance of cluster degrees of freedom and in-medium modifications in models of dilute nuclear systems. This could influence interpretations of heavy-ion collision data and equations of state relevant to astrophysics, extending standard linear response techniques to clustered regimes.

major comments (1)
  1. [Linear response and dispersion relation section] The linear response framework (Section describing the dispersion relation and response function) assumes that cluster degrees of freedom remain perturbative and that the density-dependent momentum cut-off does not introduce additional poles or non-analyticities invalidating the linearization. No explicit verification is provided that the extracted eigenvalues remain small enough for the approximation to hold at the densities where clusters appear, nor is there a comparison against a non-linear or kinetic-theory treatment. This is load-bearing for the growth-rate claims.
minor comments (2)
  1. [Abstract] The abstract would benefit from a brief mention of the specific cluster model or interaction employed to allow readers to assess the scope immediately.
  2. [Formalism] Notation for the density-dependent cut-off parameters should be introduced with an explicit equation reference in the formalism section for clarity.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for raising this important point about the validity of the linear response framework. We address the comment in detail below.

read point-by-point responses

  1. Referee: The linear response framework (Section describing the dispersion relation and response function) assumes that cluster degrees of freedom remain perturbative and that the density-dependent momentum cut-off does not introduce additional poles or non-analyticities invalidating the linearization. No explicit verification is provided that the extracted eigenvalues remain small enough for the approximation to hold at the densities where clusters appear, nor is there a comparison against a non-linear or kinetic-theory treatment. This is load-bearing for the growth-rate claims.

    Authors: The linear response approach is the standard method for identifying the onset and growth rates of spinodal instabilities in nuclear matter, as it directly yields the dispersion relation whose imaginary part gives the growth rate of unstable modes. The density-dependent momentum cut-off is constructed as a smooth function of density that suppresses cluster contributions above the Mott density; within the density and wave-number range explored, it does not generate additional poles or non-analyticities in the response function. In the revised manuscript we have added an explicit check (new paragraph and supplementary figure) demonstrating that the extracted growth rates remain small relative to the real part of the frequency (typically < 0.1 in natural units) precisely in the density window where light clusters are present, thereby confirming the perturbative character of the linearization. A direct comparison with non-linear hydrodynamic or kinetic-theory simulations lies outside the scope of the present work, which is focused on the linear regime that provides the initial growth rates seeding subsequent non-linear evolution; such a comparison would be a natural extension for future studies. revision: yes

Circularity Check

0 steps flagged

Derivation chain shows no significant circularity

full rationale

The paper applies a standard linear response framework to a model of nucleons plus light clusters, introducing a density-dependent momentum cut-off as an explicit ansatz to encode Mott dissolution. The reported effects on spinodal instabilities and growth rates follow from solving the resulting dispersion relation; no equation reduces a claimed prediction to a fitted parameter by construction, no uniqueness theorem is imported via self-citation, and no known empirical pattern is merely renamed. The derivation therefore remains self-contained against external benchmarks and receives the default non-finding.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

Only the abstract is available, so the ledger is necessarily incomplete. The density-dependent momentum cut-off is expected to introduce at least one free parameter whose value is chosen to reproduce cluster dissolution; linear response itself rests on the domain assumption that fluctuations remain small.

free parameters (1)
  • density-dependent momentum cut-off parameters
    Introduced to model in-medium Mott effects on cluster appearance; value(s) must be chosen or fitted.
axioms (1)
  • domain assumption Linear response theory remains valid once clusters and the density-dependent cut-off are added
    The method chosen in the abstract.

pith-pipeline@v0.9.0 · 5597 in / 1148 out tokens · 22478 ms · 2026-05-24T01:12:11.604725+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

59 extracted references · 59 canonical work pages

  1. [1]

    2 These temperature values mainly lie beyond the critical one for the functional F does not include any contribution from boson condensation

    Then 1 A unitary value is adopted for the speed of light. 2 These temperature values mainly lie beyond the critical one for the functional F does not include any contribution from boson condensation. Our aim here is to account, within a unified theoret- ical framework, for the formation of heavy fragments driven by the volume instabilities and the presence...

    work page

  2. [2]

    admits plane-wave solutions δfj, periodic in time with frequency ω and wave vector k, such as δfj ∼ ∑ k δf k j ei(k·r− ωt). Following a standard Landau procedure [ 44], one can derive a system of three coupled equations for neutron, proton and the light cluster species considered, which can be expressed in the following com- pact form (see Supplemental Ma...

    work page

  3. [3]

    The result from a “hybrid” case is also included (green, see text for de- tails)

    pure nucleonic matter (SNM, black); 2) for nuclear mat- ter with deuterons, including in-medium effects ( Rd ≫ Rδρb , red); 3) for nuclear matter with deuterons, neglecting ( Rd ≪ Rδρb , cyan) in-medium effects in the dynamics. The result from a “hybrid” case is also included (green, see text for de- tails). The inset shows the ρb dependence of the deuteron...

    work page

  4. [4]

    The main panel of Fig

    vanishes for ω = 0 (χ q = χ d = 1) [18, 47]. The main panel of Fig. 1 displays the spinodal border in the (ρb, T ) plane, as obtained by taking into account the lo- cal density-dependence of the cut-off in the dynamics (red line). This represents a suitable choice, considering that at the spinodal boundary one has ω = 0 , therefore the condition Rd ≫ Rδρb ...

    work page

  5. [5]

    002 fm− 3). Moreover, in the case of the full calcula- tions, a further escape and re-entrance from the region of spinodal instabilities is seen at a larger density, in som e analogy with the findings of recent works [ 48, 49], where the emergence of such a meta-stable region was discussed. /s48/s46/s48/s48 /s48/s46/s48/s50 /s48/s46/s48/s52 /s48/s46/s48/s5...

    work page

  6. [6]

    From this condition, one ex- tracts the dispersion relation, connecting the frequency ω to the wave number k

    to be equal to zero. From this condition, one ex- tracts the dispersion relation, connecting the frequency ω to the wave number k. Inside the spinodal region, a pure imaginary ω will be obtained in a certain k interval, lead- ing to an unstable growth of those modes. In Fig. 2, we show the imaginary sound velocity, Im(ω )/k (for k → 0), as a function of t...

    work page

  7. [7]

    This suppression contrasts with the increase that would occur if the variation of the cut-off had been ignored during the fragmentation pro- cess (Rd ≪ Rδρb )

    We notice that in-medium effects associated with the density de- pendence of the cut-off along the dynamics ( Rd ≫ Rδρb ) typically induce a suppression of the growth rate of the instability with respect to the case of pure nucleonic mat- ter, thus slowing down fragment formation, especially at the largest temperature considered, where deuterons sur- vive u...

    work page

  8. [8]

    C. J. Horowitz and A. Schwenk, Nucl. Phys. A 776, 55 (2006)

    work page 2006

  9. [9]

    Typel, G

    S. Typel, G. Röpke, T. Klähn, D. Blaschke, and H. H. Wolter, Phys. Rev. C 81, 015803 (2010)

    work page 2010

  10. [10]

    Gulminelli and A

    F. Gulminelli and A. R. Raduta, Phys. Rev. C 92, 055803 (2015)

    work page 2015

  11. [11]

    Bauswein, N.-U

    A. Bauswein, N.-U. F. Bastian, D. B. Blaschke, K. Chatziioannou, J. A. Clark, T. Fischer, and M. Oertel, Phys. Rev. Lett. 122, 061102 (2019)

    work page 2019

  12. [12]

    Burgio, H.-J

    G. Burgio, H.-J. Schulze, I. Vidaña, and J.-B. Wei, Prog. Part. Nucl. Phys. 120, 103879 (2021)

    work page 2021

  13. [13]

    Sedrakian, J

    A. Sedrakian, J. J. Li, and F. Weber, Prog. Part. Nucl. Phys. 131, 104041 (2023)

    work page 2023

  14. [14]

    B. Zhou, Y. Funaki, H. Horiuchi, Y.-G. Ma, G. Röpke, P. Schuck, A. Tohsaki, and T. Yamada, Nat. Commun. 14, 8206 (2023)

    work page 2023

  15. [15]

    Burrello and M

    S. Burrello and M. Grasso, Eur. Phys. J. A 58, 1 (2022)

    work page 2022

  16. [16]

    Keller, K

    J. Keller, K. Hebeler, and A. Schwenk, Phys. Rev. Lett. 130, 072701 (2023)

    work page 2023

  17. [17]

    Typel, Phys

    S. Typel, Phys. Rev. C 89, 064321 (2014)

    work page 2014

  18. [18]

    Burrello, J

    S. Burrello, J. Bonnard, and M. Grasso, Phys. Rev. C 103, 064317 (2021)

    work page 2021

  19. [19]

    I. Tews, J. Margueron, and S. Reddy, Eur. Phys. J. A 55, 97 (2019)

    work page 2019

  20. [20]

    Oertel, A

    M. Oertel, A. Pascal, M. Mancini, and J. Novak, Phys. Rev. C 102, 035802 (2020)

    work page 2020

  21. [21]

    Burrello, F

    S. Burrello, F. Gulminelli, F. Aymard, M. Colonna, and A. R. Raduta, Phys. Rev. C 92, 055804 (2015)

    work page 2015

  22. [22]

    H. D. Thi, T. Carreau, A. F. Fantina, and F. Gulminelli, Astron. Astrophys. 654, A114 (2021) . 6

    work page 2021

  23. [23]

    Grams, J

    G. Grams, J. Margueron, R. Somasundaram, and S. Reddy, Eur. Phys. J. A 58, 56 (2022)

    work page 2022

  24. [24]

    A. Rios, A. Polls, and W. Dickhoff, Phys. Rev. C 89, 044303 (2014)

    work page 2014

  25. [25]

    Chomaz, M

    P. Chomaz, M. Colonna, and J. Randrup, Phys. Rept. 389, 263 (2004)

    work page 2004

  26. [26]

    Röpke, Phys

    G. Röpke, Phys. Rev. C 92, 054001 (2015)

    work page 2015

  27. [27]

    Hempel, T

    M. Hempel, T. Fischer, J. Schaffner-Bielich, and M. Liebendörfer, Astrophys. J. 748, 70 (2012)

    work page 2012

  28. [28]

    H. Pais, F. Gulminelli, C. Providência, and G. Röpke, Phys. Rev. C 99, 055806 (2019)

    work page 2019

  29. [29]

    Röpke, Phys

    G. Röpke, Phys. Rev. C 101, 064310 (2020)

    work page 2020

  30. [30]

    Baran, M

    V. Baran, M. Colonna, V. Greco, and M. Di Toro, Phys. Rep. 410, 335 (2005)

    work page 2005

  31. [31]

    Li, L.-W

    B.-A. Li, L.-W. Chen, and C. M. Ko, Phys. Rep. 464, 113 (2008)

    work page 2008

  32. [32]

    Colonna, Prog

    M. Colonna, Prog. Part. Nucl. Phys. 113, 103775 (2020)

    work page 2020

  33. [33]

    Sorensen, K

    A. Sorensen, K. Agarwal, K. W. Brown, Z. Cha- jęcki, P. Danielewicz, C. Drischler, S. Gandolfi, J. W. Holt, M. Kaminski, C.-M. Ko, et al., Prog. Part. Nucl. Phys. 134, 104080 (2024)

    work page 2024

  34. [34]

    Wolter et al

    H. Wolter et al. (TMEP), Prog. Part. Nucl. Phys. 125, 103962 (2022)

    work page 2022

  35. [35]

    Xu, Prog

    J. Xu, Prog. Part. Nucl. Phys. 106, 312 (2019)

    work page 2019

  36. [36]

    Danielewicz and G

    P. Danielewicz and G. F. Bertsch, Nucl. Phys. A 533, 712 (1991)

    work page 1991

  37. [37]

    A. Ono, J. Phys.: Conf. Ser. 420, 012103 (2013)

    work page 2013

  38. [38]

    Oliinychenko, L.-G

    D. Oliinychenko, L.-G. Pang, H. Elfner, and V. Koch, Phys. Rev. C 99, 044907 (2019)

    work page 2019

  39. [39]

    Ono, Prog

    A. Ono, Prog. Part. Nucl. Phys. 105, 139 (2019)

    work page 2019

  40. [40]

    Wang, Y.-G

    R. Wang, Y.-G. Ma, L.-W. Chen, C. M. Ko, K.-J. Sun, and Z. Zhang, Phys. Rev. C 108, L031601 (2023)

    work page 2023

  41. [41]

    G. Coci, S. Gläßel, V. Kireyeu, J. Aichelin, C. Blume, E. Bratkovskaya, V. Kolesnikov, and V. Voronyuk, Phys. Rev. C 108, 014902 (2023)

    work page 2023

  42. [42]

    Cheng and Z.-Q

    H.-G. Cheng and Z.-Q. Feng, Phys. Rev. C 109, L021602 (2024)

    work page 2024

  43. [43]

    K.-J. Sun, R. Wang, C. M. Ko, Y.-G. Ma, and C. Shen, Nat. Commun. 15, 1074 (2024)

    work page 2024

  44. [44]

    Napolitani and M

    P. Napolitani and M. Colonna, Phys. Lett. B 726, 382 (2013)

    work page 2013

  45. [45]

    Lin and P

    H. Lin and P. Danielewicz, Phys. Rev. C 99, 024612 (2019)

    work page 2019

  46. [46]

    Röpke, Nucl

    G. Röpke, Nucl. Phys. A 867, 66 (2011)

    work page 2011

  47. [47]

    Burrello and S

    S. Burrello and S. Typel, Eur. Phys. J. A 58, 120 (2022)

    work page 2022

  48. [48]

    Wu, S.-B

    X.-H. Wu, S.-B. Wang, A. Sedrakian, and G. Röpke, J. Low Temp. Phys. 189, 133 (2017)

    work page 2017

  49. [49]

    Colonna and F

    M. Colonna and F. Matera, Phys. Rev. C 77, 064606 (2008)

    work page 2008

  50. [50]

    Burrello, M

    S. Burrello, M. Colonna, G. Colò, D. Lacroix, X. Roca-Maza, G. Scamps, and H. Zheng, Phys. Rev. C 99, 054314 (2019)

    work page 2019

  51. [51]

    Landau, Sov

    L. Landau, Sov. Phys. JETP 8, 70 (1959)

    work page 1959

  52. [52]

    Qin et al., Phys

    L. Qin et al., Phys. Rev. Lett. 108, 172701 (2012)

    work page 2012

  53. [53]

    Kuhrts, M

    C. Kuhrts, M. Beyer, P. Danielewicz, and G. Röpke, Phys. Rev. C 63, 034605 (2001)

    work page 2001

  54. [54]

    Burrello, M

    S. Burrello, M. Colonna, and F. Matera, Phys. Rev. C 89, 057604 (2014)

    work page 2014

  55. [55]

    Röpke, D

    G. Röpke, D. Voskresensky, I. Kryukov, and D. Blaschke, Nucl. Phys. A 970, 224 (2018)

    work page 2018

  56. [56]

    Voskresensky,Prog

    D. Voskresensky,Prog. Part. Nucl. Phys. 130, 104030 (2023)

    work page 2023

  57. [57]

    Borderie, N

    B. Borderie, N. Le Neindre, M. Rivet, P. Dés- esquelles, E. Bonnet, R. Bougault, A. Chbihi, D. Dell’Aquila, Q. Fable, J. Frankland, et al., Phys. Lett. B 782, 291 (2018)

    work page 2018

  58. [58]

    Huth et al., Nature 606, 276 (2022)

    S. Huth et al., Nature 606, 276 (2022)

    work page 2022

  59. [59]

    C. Y. Tsang, M. B. Tsang, W. G. Lynch, R. Kumar, and C. J. Horowitz, Nat. Astron. 8, 328 (2024) . 7 Supplemental Material I. ENERGY DENSITY FUNCTIONAL AND SINGLE P AR TICLE ENERGY In the expression of the free energy density functional F = E − T S, the energy density functional E = K + U is the sum of the kinetic energy density K = ∑ j gj ∫ Λj dp (2π ℏ)3 ...

    work page 2024