A time-dependent Hartree-Fock study of triple-alpha dynamics

J. L. Willerton; P. D. Stevenson

REVIEW 3 major objections 5 minor 16 references

A time-dependent Hartree-Fock study of triple-alpha dynamics

T0 review · 3 major / 5 minor · reviewed 2026-08-14 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read A time-dependent Hartree-Fock simulation of the triple-alpha reaction finds that carbon-12 forms either a linear-chain vibrational state near 9 MeV or a triangular vibrational state near 4 MeV, depending on the orientation of the incoming…

desk verdict Honest, exploratory TDHF study of triple-alpha dynamics whose qualitative mode structure is interesting but whose quoted 9 and 4 MeV energies are model-dependent numbers, not physical predictions. read the letter →

arxiv 1909.01924 v2 pith:ARVZRYSB submitted 2019-09-04 nucl-th

P. D. Stevenson , J. L. Willerton This is my paper · ORCID

classification nucl-th

keywords triple-alphareactiontime-dependentHartree-Fockcarbon-12clusterstateslinearchainstatetriangularvibrationberyllium-8Skyrmeinteractionquadrupoleoscillations

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper follows the full triple-alpha path in time-dependent Hartree-Fock: two helium-4 nuclei fuse into beryllium-8, and a third helium-4 fuses with that beryllium-8 to make carbon-12. Its central claim is that the shape of the resulting carbon-12 depends on how the deformed beryllium-8 is oriented when the third alpha arrives: tip-on collisions produce a linear chain that vibrates near 9 MeV, side-on collisions produce a compact triangular vibration near 4 MeV, and the chain state can decay into the triangular state. This matters because the triple-alpha reaction is the stellar bottleneck that turns helium into carbon, and the calculation provides a fully dynamical mean-field picture of the cluster states involved. The authors note that the chosen interaction underbinds helium-4 by 40 percent, so the absolute energies are not expected to be final.

What carries the argument

The carrying mechanism is time-dependent Hartree-Fock (TDHF), a mean-field approximation in which each nucleon moves in the self-consistent average field generated by all the others, evolved on a spatial grid. The primary observable is the mass quadrupole moment Q of the total density, defined in the paper's Eq. (1); its time series is Fourier-transformed to read off vibrational energies of the compound nucleus. The decisive control parameter is the orientation of the deformed 8Be* nucleus ('tip' versus 'side') at the moment the third alpha arrives, together with the impact parameter and the snapshot of the beryllium internal state.

What would settle it

Repeat the same two-stage TDHF calculation with a Skyrme energy functional that reproduces the experimental helium-4 binding energy of 28.30 MeV, and inspect the Fourier spectrum of the carbon-12 quadrupole moment. If the peak near 9 MeV and the peak near 4 MeV do not survive, or if no chain-to-triangle transition is observed, the central claim is refuted. A cheaper check is to see whether the 8Be* peaks at 8.5 and 14.7 MeV move toward the known 3.03 and 11.35 MeV levels under a better-binding force; if they shift while the 12C peaks do not, the two-alpha and three-alpha interpretations are inconsistent.

Watch

Extended reading notes

Core claim

The paper reports that in two-stage TDHF simulations of 4He+4He -> 8Be* followed by 4He+8Be* -> 12C, the fused carbon nucleus supports identifiable large-amplitude cluster vibrations rather than a featureless mean-field state. Head-on tip collisions of the third alpha with the deformed beryllium leave the density in a linear chain of three alphas oscillating with a Fourier peak around 9 MeV, while side collisions form a compact triangular arrangement oscillating around 4 MeV. In a tip collision with impact parameter b = 1 fm, the chain state decays into the triangular state within about 2 zs. In the two-alpha stage, the calculation extracts beryllium-8 vibrational peaks at 8.5 and 14.7 MeV, which the paper compares with the known levels at 3.03 and 11.35 MeV and attributes to the same interaction deficiency.

Load-bearing premise

The results depend on the chosen nuclear force binding helium-4 at 17.67 MeV instead of the experimental 28.30 MeV, a 40 percent underbinding that the authors concede could strongly influence the energies and even the qualitative vibrational structure they report.

Editorial extensions

If this is right

If the central claim is right, the triple-alpha reaction has at least two distinct doorway configurations in 12C, selected by the geometry of the 8Be+alpha encounter.
The chain vibrational state at about 9 MeV is not the endpoint: it decays into the compact triangular state on a timescale of about 2 zs, so the longer-lived product of the fusion is the triangular configuration.
The extracted energies of 9 and 4 MeV give concrete benchmark values that any future microscopic model of 12C cluster states would need to reproduce, even though the absolute scale is hostage to the interaction's alpha underbinding.
The failure of the two-alpha calculation to reproduce the known 3.03 MeV 2+ state in 8Be implies that, with this interaction, the mean-field reaction mechanism either does not excite that state or places it at an overestimated energy.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

If the orientation dependence is real, the effective triple-alpha rate in a stellar plasma may depend on the alignment distribution of 8Be at the moment the third alpha approaches, a dependence not present in standard equilibrium rate formulas.
The systematic 40 percent alpha underbinding most likely shifts all extracted vibrational energies, so the qualitative chain-versus-triangle structure is the more robust claim than the specific 9 and 4 MeV values.
A natural next calculation is a systematic scan of the 8Be internal vibrational phase at impact: the paper samples only two snapshots (2000 and 4000 iterations), and the chain-to-triangle branching ratio could vary across the full beryllium oscillation cycle.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

You should know two things about this one. It's a short proceedings paper that does something genuinely new in the TDHF context—two-step 4He+4He→8Be* then 4He+8Be*, with explicit tip vs side orientation dependence—and it is refreshingly upfront about its own limits. The headline numbers, though, are not established as physical: with SLy4d the alpha is underbound by 40% and the same method puts the 8Be first 2+ at 8.5 MeV instead of 3.03 MeV, so the ~9 and ~4 MeV 12C "vibrations" are at best model-dependent frequencies.

Previous TDHF work (Umar et al. 2010) already saw chain and triangular configurations; the authors cite it and say their results qualitatively agree. What they add is the two-step collision protocol (saving the 8Be* wave function at two chosen times), the orientation study (tip vs side), and the Fourier analysis of the quadrupole moment, which lets them identify a chain mode at ~9 MeV decaying into a triangular mode at ~4 MeV in a b=1 tip collision. The computations are standard Sky3D with SLy4d, and the paper honestly flags the underbinding, the arbitrary choice of starting configurations, and the lack of angular momentum projection. For a proceedings contribution, that's decent scientific practice.

The soft spots are real but proportionate. The abstract states the 9 and 4 MeV energies without caveats; that's overreach given the 8Be benchmark already misses by 5.5 MeV. The spectral analysis uses a single Slater determinant trajectory and a single observable (Q), so peak assignment to distinct chain/triangular states is reasonable but not rigorous; the authors acknowledge this. The parameter scan is thin (two orientations, b up to 1 fm, Ecm=2 MeV, two 8Be* snapshots), so one cannot assess how robust the modes are. These are limitations typical of a first exploratory study, not fatal flaws. The claim that a force fitting alpha binding is needed is worth taking seriously; the authors themselves say so.

Who is this for? People working on TDHF and alpha clustering will want to see it as evidence that the two-step protocol and orientation dependence matter. It does not claim to resolve the physical 12C spectrum, and the reader's conditional verdict is fair. I'd send it to peer review as a proceedings contribution; with revision, mainly to soften the abstract and add a benchmark comparison to 8Be, it would be a solid contribution. The citation pattern looks fine; the self-references to Sky3D and related Surrey work are appropriate for a methods paper.

Referee Report

3 major / 5 minor

Summary. The paper reports time-dependent Hartree-Fock (TDHF) calculations of the triple-alpha reaction, modeled as a two-step process: first 4He+4He fusion to form 8Be*, then 4He+8Be* collisions with different orientations and impact parameters. Using the SLy4d Skyrme interaction and the Sky3D code, the authors extract Fourier power spectra of the quadrupole moment of the fused 12C system and identify two vibrational modes: a linear chain state at about 9 MeV and a triangular state at about 4 MeV, with transitions from the chain to the triangular configuration observed in some trajectories. The paper explicitly acknowledges the 40% underbinding of 4He by SLy4d and the lack of angular momentum projection.

Significance. If the results were quantitatively reliable, the identification of distinct chain and triangular cluster vibrational modes in a TDHF treatment of the triple-alpha reaction would be a useful contribution to the study of cluster structures in 12C, complementing prior work by Umar et al. and others. The two-step collision setup and the use of Fourier analysis of a dynamical observable are methodologically interesting. However, the significance is substantially tempered by the model's known deficiencies: the interaction underbinds 4He by about 40%, and the internal calibration against the 8Be spectrum shows a ~5.5 MeV discrepancy for the first excited state. These issues undermine the quantitative claim of vibrational energies at ~9 and ~4 MeV, leaving the paper as a qualitative exploratory study rather than a quantitative prediction. The authors are honest about the limitations, which is commendable, but the abstract and conclusion present the energies without those caveats.

major comments (3)

[Section 3, first paragraph; Section 3.1; Abstract] The central quantitative claim that the chain and triangular vibrations occur at ~9 and ~4 MeV is not supported by the model's internal consistency. The SLy4d interaction binds the alpha particle at 17.67 MeV against the experimental 28.30 MeV, a 40% underbinding that the authors themselves state 'could clearly have a strong influence on the results.' The one available internal control, the 2-alpha calculation in Section 3.1, extracts the first excited state of 8Be at 8.5 MeV, whereas the known 2+ state is at 3.03 MeV, a discrepancy of 5.5 MeV. Since the 12C peaks are extracted by the same Fourier method with the same interaction, there is no reason to expect them to be accurate to better than a similar scale. The abstract's unconditional wording ('occur at ~9 and 4 MeV') is therefore premature. The authors should either provide evidence that the peak energies are robust to the interaction choice (e.g., by testing a force that reproduces alpha binding) or explicitly reframe the claims as model-dependent values within the SLy4d TDHF framework.
[Section 3.2, especially the final paragraph; Figure 2 insets; Conclusion] The identification of the two Fourier peaks with distinct chain and triangular states of 12C is not established. The assignment rests on visual inspection of density snapshots and on the time evolution of a single collective observable, the quadrupole moment Q. The authors explicitly defer angular momentum projection and state that 'a more sophisticated treatment than mixing via Fourier analysis would be needed to obtain definite spins for each state.' Without such projection or another quantum-number assignment, the peaks should be described as 'vibrational modes of the TDHF trajectory' rather than 'states of 12C' in the abstract and conclusion. As written, the conclusion repeats the unqualified claim that 'identifiable chain and triangular vibrational states at around 9 and 4 MeV respectively are found.'
[Section 3.2, first paragraph; Figure 2] The generality of the orientation-dependent claim is limited by the small and partly arbitrary parameter sampling. The two 8Be* starting configurations are 'somewhat arbitrarily chosen' (iterations 2000 and 4000), only two orientations ('tip' and 'side') are considered, the impact parameter is restricted to b=0 and b=1 fm, and the center-of-mass energy is fixed at Ecm=2 MeV. While the paper appropriately calls for a fuller study, the abstract's statement that 'depending on the orientation of the initial state' one finds chain or triangular states is a strong generalization from this limited set of trajectories. The authors should quantify the sensitivity of the claimed modes to these choices, or restrict the claim to the specific initial conditions studied.

minor comments (5)

[Section 2, Methodology] There is a typo: 'as obtained form the static Hartree-Fock calculation' should be 'as obtained from the static Hartree-Fock calculation.'
[Section 3.1, paragraph after Figure 1] The sentence 'The second observed excited state is 8Be is a broad 4+ resonance' is ungrammatical; it should be 'The second observed excited state of 8Be is a broad 4+ resonance.'
[Figure 2 and its insets] The inset spectra are difficult to read at the printed size, and the yellow and pink dotted lines may be indistinguishable in grayscale. Plotting the spectra on a common energy axis with labeled peaks would improve clarity.
[Section 3.2, paragraph beginning 'Mixing of the mean-field...'] The phrase 'mixing of the mean-field TDHF configurations via a Fourier spectrum analysis' is imprecise; a Fourier transform of the quadrupole time series is not a mixing of configurations. Suggest rewording to 'spectral analysis of the quadrupole moment time series.'
[Abstract and Conclusion] The abstract uses '~9 and 4 MeV' while the conclusion uses 'around 9 and 4 MeV' and Section 3.2 uses 'around 9 MeV' and 'around 4 MeV'; the wording should be made consistent, ideally with the caveat that these are model-dependent values.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the 9 and 4 MeV vibrational energies are extracted from Fourier power spectra of the TDHF quadrupole-moment time series, not fitted or defined into existence.

full rationale

The paper's derivation chain runs from static Hartree-Fock 4He ground states, through TDHF 4He+4He fusion, saved 8Be* configurations, TDHF 4He+8Be* collisions, and finally Fourier analysis of the quadrupole-moment time series Q(t) to obtain the reported ~9 and ~4 MeV peaks. No parameter is fitted to those energies; the peak positions emerge from the dynamics, and the assignment to chain versus triangular motion is made from density snapshots and the time evolution of the spectra. The self-citations present (Sky3D code references [10,11], the Skyrme-interaction review [9], the time-reversal validation [7], and the continuum discretisation discussion [14]) are methodological or validation references and are not load-bearing for the central claim. The paper's own stated limitations--SLy4d underbinding 4He by about 40%, the 8Be first excited state appearing at 8.5 MeV instead of 3.03 MeV, and the absence of angular momentum projection--are accuracy and interpretation concerns, not circularity. The central result is therefore not equivalent to its inputs by construction.

Assumptions & free parameters 5 free parameters · 4 assumptions · 0 invented entities

The paper's central results rest on a set of tuned numerical choices (collision energies, impact parameters, selected 8Be configurations, and grid resolution) and on domain assumptions about TDHF and the SLy4d interaction. No new physical entities are introduced. The underbinding of the alpha particle is the most consequential admitted assumption, as it directly shifts the energy scale of all reported vibrational states.

free parameters (5)

8Be* starting configuration = iterations 2000 and 4000 of the 2-alpha TDHF trajectory
Chosen 'somewhat arbitrarily' (Sec. 3.2) to represent the compound beryllium nucleus; one near a maximum and one near a minimum of the quadrupole moment. The results depend on this choice.
Triple-alpha collision energy Ecm = 2.0 MeV
Fixed close to the Coulomb barrier between 1 and 2 MeV to favor fusion; the paper defers study of other energies (Sec. 3.2).
Impact parameter b = 0 and 1 fm
Only these two values are sampled; the full parameter space (Ecm, b, orientation) is conceded to need mapping (Sec. 3.2).
2-alpha collision energy Ecm = 1.0 MeV
Single energy used to form the 8Be*; chosen to keep the two alphas fused (Sec. 3.1).
Grid spacing = 1 fm
Numerical resolution in Sky3D; peak fragmentation is attributed to box discretization (Sec. 3.1), so resolution affects the extracted energies.

assumptions (4)

domain assumption TDHF mean-field dynamics with a Skyrme interaction captures the essential cluster structure and reaction dynamics of the triple-alpha process.
The entire analysis relies on TDHF as implemented in Sky3D (Sec. 2); no beyond-mean-field corrections are used.
domain assumption The SLy4d interaction is adequate for this study despite underbinding 4He by 40 percent.
Section 3 states the underbinding 'could clearly have a strong influence on the results', yet the calculation proceeds with this force.
domain assumption Fourier analysis of the quadrupole moment time series yields the energies of the chain and triangular vibrational states.
Sec. 3.2 uses the power spectrum to assign 9 and 4 MeV; the paper notes spins and rigorous spectrum require projection.
ad hoc to paper Beryllium-8 configurations at iterations 2000 and 4000 are representative starting points for the triple-alpha collision.
Sec. 3.2 labels these 'somewhat arbitrarily chosen'; there is no systematic scan or justification that they sample the physically relevant distribution.

how reviews work

0 comments

Cite this review

Pith. "Pith review of A time-dependent Hartree-Fock study of triple-alpha dynamics." pith.science (2026). https://pith.science/paper/ARVZRYSB

@misc{pith2026190901924,
  author       = {Pith},
  title        = {Pith review of: A time-dependent Hartree-Fock study of triple-alpha dynamics},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/ARVZRYSB}},
  note         = {Machine review of arXiv:1909.01924}
}

read the original abstract

Time-dependent Hartree-Fock calculations have been performed for fusion reactions of He-4 + He-4 -> Be*-8, followed by He-4 + Be*-8 . Depending on the orientation of the initial state, a linear chain vibrational state or a triangular vibration is found in 12C, with transitions between these states observed. The vibrations of the linear chain state and the triangular state occur at ~9 and 4 MeV respectively.

Figures

Figures reproduced from arXiv: 1909.01924 by the authors.

**Figure 2.** Quadrupole moment of the matter density during collisions of [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗

**Figure 3.** Snapshots of the time evolution of the total density for four of the trajectories [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

16 extracted references · 14 canonical work pages

[1]

E. M. Burbidge, G. R. Burbidge, W . A. Fowler and F . Hoyle,Synthesis of the elements in stars , Rev. Mod. Phys. 29, 547 (1957), doi:10.1103 /RevModPhys.29.547

work page 1957
[2]

Freer and H

M. Freer and H. Fynbo, The Hoyle state in 12C, Progress in Particle and Nuclear Physics 78, 1 (2014), doi:10.1016 /j.ppnp.2014.06.001

work page 2014
[3]

Freer, H

M. Freer, H. Horiuchi, Y. Kanada-En’Yo, D. Lee and U. G. Meißner,Microscopic clustering in light nuclei, Rev. Mod. Phys. 90(3), 35004 (2018), doi:10.1103 /RevModPhys.90.035004

work page 2018
[4]

A. S. Umar, J. A. Maruhn, N. Itagaki and V . E. Oberacker, Microscopic Study of the Triple- α Reaction, Phys. Rev. Lett. 104(21), 212503 (2010), doi:10.1103 /PhysRevLett.104.212503

work page 2010
[5]

Ichikawa, N

T . Ichikawa, N. Itagaki, N. Loebl, J. A. Maruhn, V . E. Oberacker, S. Ohkubo, B. Schuetrumpf and A. S. Umar, Static and Dynamic Chain Structures in the Mean-Field Theory , EPJ Web Conf. 17, 07002 (2011), doi:10.1051 /epjconf/20111707002

work page arXiv 2011
[6]

A. S. Umar and M. R. Strayer, Nuclear shape-isomeric vibrations , Phys. Lett. B 171(4), 353 (1986), doi:10.1016 /0370-2693(86)91419-X

work page 1986
[7]

Iwata and P

Y. Iwata and P . Stevenson, Conditional recovery of time-reversal symmetry in many nucleus systems, New Journal of Physics 21(4), 043010 (2019), doi:10.1088 /1367-2630/ab0e58

work page 2019
[8]

Simenel, Nuclear quantum many-body dynamics , Eur

C. Simenel, Nuclear quantum many-body dynamics , Eur. Phys. J. A 48(11), 152 (2012), doi:10.1140/epja/i2012-12152-0

work page doi:10.1140/epja/i2012-12152-0 2012

Show all 16 references

[9]

P . D. Stevenson and M. C. Barton, Low-energy heavy-ion reactions and the Skyrme effective interaction , Progress in Particle and Nuclear Physics 104, 142 (2019), doi:10.1016/j.ppnp.2018.09.002

2019 doi
[10]

J. A. Maruhn, P .-G. Reinhard, P . D. Stevenson and A. S. Umar,The TDHF code Sky3D, Comput. Phys. Commun. 185(7), 2195 (2014), doi:10.1016 /j.cpc.2014.04.008

2014
[11]

Schuetrumpf, P .-G

B. Schuetrumpf, P .-G. Reinhard, P . D. Stevenson, A. S. Umar and J. A. Maruhn, The TDHF code Sky3D version 1.1 , Comput. Phys. Commun. 229, 211 (2018), doi:10.1016/j.cpc.2018.03.012

2018 doi
[12]

K.-H. Kim, T . Otsuka and P . Bonche, Three-dimensional TDHF calculations for reactions of unstable nuclei, J. Phys. G 23(10), 1267 (1997), doi:10.1088 /0954-3899/23/10/014

1997
[13]

M. Wang, G. Audi, F . G. Kondev, W . J. Huang, S. Naimi and X. Xu, The AME2016 atomic mass evaluation (II) T ables, graphs and references , Chin. Phys. C 41(3), 030003 (2017), doi:10.1088/1674-1137/41/3/030003

2017 doi
[14]

Pardi, P

C. Pardi, P . Stevenson and K. Xu, Extension of the continuum time-dependent Hartree-Fock method to proton states , Physical Review E - Statistical, Nonlinear, and Soft Matter Physics 89(3), 033312 (2014), doi:10.1103 /PhysRevE.89.033312

2014
[15]

D. R. Tilley , J. H. Kelley , J. L. Godwin, D. J. Millener, J. E. Purcell, C. G. Sheu and H. R. Weller, Energy levels of light nuclei A = 8,9,10, Nuclear Physics A 745(3-4), 155 (2004), doi:10.1016/j.nuclphysa.2004.09.059. 6 SciPost Physics Proceedings Submission

2004 doi
[16]

R. Imai, T . Tada and M. Kimura,Real-time evolution method and its application to the 3α cluster system, Physical Review C 99(6), 064327 (2019), doi:10.1103 /PhysRevC.99.064327. 7

2019

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.

[1] [1]

E. M. Burbidge, G. R. Burbidge, W . A. Fowler and F . Hoyle,Synthesis of the elements in stars , Rev. Mod. Phys. 29, 547 (1957), doi:10.1103 /RevModPhys.29.547

work page 1957

[2] [2]

Freer and H

M. Freer and H. Fynbo, The Hoyle state in 12C, Progress in Particle and Nuclear Physics 78, 1 (2014), doi:10.1016 /j.ppnp.2014.06.001

work page 2014

[3] [3]

Freer, H

M. Freer, H. Horiuchi, Y. Kanada-En’Yo, D. Lee and U. G. Meißner,Microscopic clustering in light nuclei, Rev. Mod. Phys. 90(3), 35004 (2018), doi:10.1103 /RevModPhys.90.035004

work page 2018

[4] [4]

A. S. Umar, J. A. Maruhn, N. Itagaki and V . E. Oberacker, Microscopic Study of the Triple- α Reaction, Phys. Rev. Lett. 104(21), 212503 (2010), doi:10.1103 /PhysRevLett.104.212503

work page 2010

[5] [5]

Ichikawa, N

T . Ichikawa, N. Itagaki, N. Loebl, J. A. Maruhn, V . E. Oberacker, S. Ohkubo, B. Schuetrumpf and A. S. Umar, Static and Dynamic Chain Structures in the Mean-Field Theory , EPJ Web Conf. 17, 07002 (2011), doi:10.1051 /epjconf/20111707002

work page arXiv 2011

[6] [6]

A. S. Umar and M. R. Strayer, Nuclear shape-isomeric vibrations , Phys. Lett. B 171(4), 353 (1986), doi:10.1016 /0370-2693(86)91419-X

work page 1986

[7] [7]

Iwata and P

Y. Iwata and P . Stevenson, Conditional recovery of time-reversal symmetry in many nucleus systems, New Journal of Physics 21(4), 043010 (2019), doi:10.1088 /1367-2630/ab0e58

work page 2019

[8] [8]

Simenel, Nuclear quantum many-body dynamics , Eur

C. Simenel, Nuclear quantum many-body dynamics , Eur. Phys. J. A 48(11), 152 (2012), doi:10.1140/epja/i2012-12152-0

work page doi:10.1140/epja/i2012-12152-0 2012