{"id":"5ecfc2b1-fccb-4aec-b23d-b17a48b9c9b0","arxiv_id":"2608.04102","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A 2PI-Hartree effective potential for two mixing scalars is renormalized and used to show that self-consistent thermal resummation can substantially alter predicted phase transition strengths and gravitational wave spectra.","lead":"This paper builds a finite-temperature effective potential for two interacting scalar fields using the two-particle irreducible (2PI) formalism, which resums thermal masses self-consistently. The authors show that this Hartree-resummed potential can change the predicted strength of cosmological phase transitions and the associated gravitational wave signal by orders of magnitude compared with standard resummation schemes.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Renormalization of mixing 2PI-Hartree potential relies on unproven divergence structure (Eq. 3.21); additional subdivergences would invalidate Eqs. (4.5)/(5.3).","rationale":"The reader's weakest_assumption identifies exactly the point on which the entire formalism pivots: the completeness of the divergence structure (3.21) and the sufficiency of the cancellation conditions (3.23)-(3.25). I agree this is the most load-bearing concern. The paper's renormalization procedure is plausible and follows the single-field and O(N) literature, but the extension to two mixing fields is new, and the proof that no other subdivergences arise is not supplied. In particular, the effective potential (4.5) is obtained after assuming the local correlation function's divergence is fully captured by m^2 Delta_epsilon; if any additional divergence enters through the mixing or through the resummed vertex functions, the resulting potential would not be finite and all numerical results would be void. This is a correctness risk, not a disagreement with consensus: the standard 2PI renormalization literature supports the single-field case, but the two-field mixing extension is not independently verified. The concrete test—explicit pole cancellation in a fourth derivative of the potential at an arbitrary field point—would settle the issue. Other concerns (e.g., non-identical parameter sets in the two-step comparison, or the finite-temperature scale invariance) are secondary: the parameter mismatch affects only the strength of the demonstrated phenomenological difference, and the thermal scale invariance is likely to follow from the vacuum proof with a Q-independent thermal correction. Thus the reader's conditional verdict remains appropriate, with the proposed check as a necessary condition for acceptance.","tokens_in":42017,"tokens_out":26660,"duration_ms":266039,"concrete_test":"Compute the fourth-order field derivative d^4 V_HT/dphi^4 at an off-axis field point (phi,chi) using the analytical mass derivatives of Appendix C, in dimensional regularization d=4-2epsilon with the counterterms solved from (3.24)-(3.25). Check explicitly that all 1/epsilon and 1/epsilon^2 poles cancel to the same order as at the minimum (v,w), separately for the diagonal and off-diagonal mass sectors. A surviving pole would falsify the assumption (3.21) and invalidate the effective potential.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that Eq. (5.3) is a self-consistent, all-temperature, scale-invariant effective potential—depends entirely on the assertion in Sec. 3.2.1 that every UV divergence in the Hartree-truncated 2PI action has the local form m^2_{alpha beta}Delta_epsilon (Eq. 3.21), so that the counterterm conditions (3.23)-(3.25) cancel all subdivergences for arbitrary field configurations. This is stated as an assumption, not derived: the text says 'if the theory is to be renormalizable, the mass squared matrix elements must be finite...' and then takes (3.21) as the divergence structure. For the two-field mixing case, which is the paper's novel contribution, the off-diagonal propagators and the resummed vertex functions V_{alpha beta bar-gamma} of Appendix A introduce extra opportunities for subdivergences, including possible derivative terms (e.g., proportional to d_mu m^2 or d_mu phi) that are not present in the assumed structure. If such additional divergences exist, the cancellation conditions (3.23)-(3.25) would be incomplete, and the effective potential (4.5) and its finite-temperature extension (5.3) would retain 1/epsilon poles, invalidating all numerical predictions. No independent verification is provided that the assumed divergence form is complete for the two-field Hartree resummation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops the 2PI Hartree approximation for a model with two mixing real scalar fields, renormalizing the 2PI effective action through the cancellation-of-subdivergences method and connecting auxiliary MS parameters to physical p^2=0 masses and fourth derivatives of the potential. It derives a vacuum effective potential (4.5) and a finite-temperature version (5.3), argued to be scale-invariant and free of high-temperature approximations. Using benchmark points, it compares one- and two-step phase transitions and gravitational wave spectra with the ring and Parwani resummation schemes. The central claim is that the Hartree potential is self-consistent, valid for all temperatures, and can change predicted transition strengths and gravitational wave amplitudes by orders of magnitude relative to conventional schemes.","tokens_in":42409,"tokens_out":7140,"duration_ms":82088,"significance":"If the renormalization proof holds, this is a valuable methodological contribution: it extends 2PI renormalization to the multi-field mixing case, gives a concrete physical-parameter mapping, and provides a resummation scheme for cosmological phase transitions that does not rely on the high-temperature expansion. The explicit scale-invariance proof in Appendix B, the analytical derivative formulas in Appendix C, and the bounce/GW numerical implementation are substantial strengths. The main uncertainties are the completeness of the assumed divergence structure in Eq. (3.21), which is stated rather than proven, and the fact that the two-step numerical comparison uses different λφχ values for the two schemes. These issues prevent acceptance in the present form.","major_comments":[{"comment":"The UV finiteness of VHT rests on the assertion that the only divergence in the local correlation functions is m^2_{αβ}Δε. This is demonstrated for constant fields in Eqs. (3.17)–(3.20), and then asserted to hold for arbitrary field configurations and for the two-field mixing case. The text itself says 'if the theory is to be renormalizable' (Sec. 3.2), which indicates this is an assumption rather than a derivation. In a space-dependent background the resummed propagator can develop additional subdivergences not of the form m^2_{αβ}Δε, for instance terms involving derivatives of m^2 or of the fields, and the operator structure in Appendix A introduces further opportunities for such terms. If any such divergence exists, the cancellation conditions (3.23)–(3.25) are incomplete and the potential (4.5) and its thermal extension (5.3) would retain 1/ε poles. Please provide a proof of completeness of the divergence structure for the two-field Hartree case, or explicitly restrict the claims to constant backgrounds and qualify the broader renormalizability statement.","section":"Sec. 3.2.1, Eqs. (3.21)–(3.25)"},{"comment":"The statement that the Hartree potential (5.3) is 'valid for all temperatures' is stronger than what the Hartree truncation supports. The scheme avoids the high-temperature expansion and retains Boltzmann suppression, which is a genuine improvement. However, the Hartree approximation is the leading-order 2PI truncation, and the paper provides no estimate of the neglected higher-order 2PI diagrams or of the resulting accuracy of thermodynamic quantities. I recommend replacing 'valid for all temperatures' with a more guarded formulation, such as 'free of the high-temperature expansion within the Hartree truncation', unless an explicit error estimate or convergence test is added.","section":"Sec. 5.1 and Abstract"},{"comment":"The two-step comparison is not performed at the same physical input parameters. The Parwani case uses λφχ = 1.95 while the Hartree case uses λφχ = 1.92, and the text states that for λφχ = 1.92 the Parwani potential gives a one-step transition. Consequently, the differences in T_n, α*, β/H*, and the GW spectra in Figs. 8 and 9 may reflect the change in λφχ rather than the difference in resummation scheme. The authors acknowledge this in words, but the conclusion that the Hartree potential 'tends to predict a much stronger transition' in the two-step scenario is not established by a controlled comparison. Please either find overlapping parameter regions where both schemes give two-step transitions, or present the case as an illustration with different physical inputs and avoid drawing scheme-dependent conclusions from it.","section":"Sec. 6.2, Table 3"}],"minor_comments":[{"comment":"In Eq. (3.31b), the term displayed as (λ(4)φχ + δ(4)λφχ)ϕ^2 appears to be a typo for φ^2 or χ^2; please correct the notation for consistency with the surrounding equations.","section":"Sec. 3.3, Eq. (3.31b)"},{"comment":"The column headers list 'λ(0)(m1)' twice; the last entry should be explicitly labelled as ¯λ(0)φχ (or ¯λ(0)φχ evaluated at the indicated scale) to avoid ambiguity.","section":"Tables 1–3"},{"comment":"The numerical inversion of the hatted couplings via a guessed λ(4g) and repeated solution of (4.14c) would benefit from a statement of the convergence criterion and, if available, an estimate of the numerical uncertainty; the discontinuous jumps of the global minimum shown in Fig. 3 suggest that such robustness information is relevant.","section":"Sec. 4.1, around Eq. (4.15)"},{"comment":"The final paragraph speculates about field- and temperature-dependent wave-function renormalization factors that are not used elsewhere in the paper; consider moving this discussion to the conclusions or removing it to keep the appendix focused on the quantities actually employed.","section":"Appendix A, final paragraph"},{"comment":"The bounce action (5.12) uses constant Zα,2 factors evaluated at the vacuum minimum, while Appendix A derives them at that point; please state this explicitly in Sec. 5.3 so that the reader is not left to infer it.","section":"Secs. 5.3 and 6.1"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of JHEP and the authors are honest about several caveats, including the 'guess' nature of Eq. (3.34b) and the non-overlapping parameter sets in Sec. 6.2. The main technical risk is the unproven completeness of the divergence structure in Sec. 3.2.1; I would ask the authors to address that point directly, or to weaken the corresponding claims. There are no concerns about citation practices beyond the expected self-citations."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth a serious look. The genuinely new thing is the renormalized 2PI-Hartree effective potential for two explicitly mixing scalars, with a physical parameter matching that goes beyond the single-field and O(N) cases. The derivation is detailed and mostly self-contained: the split of local correlation functions, the cancellation of subdivergences, the running of the auxiliary MS couplings, and the proof of scale invariance are all laid out carefully. The connection to physical four-point functions via fourth derivatives of the potential is a practical step that will be useful for phenomenology. Appendix A, with the wave-function renormalization factors for the classical fields, is a real addition. The numerical comparison with ring and Parwani schemes makes the point that the choice of resummation can change PT strength and GW amplitude by orders of magnitude, which is worth knowing even if the model is a toy.\n\nThe soft spots are real but not fatal. First, the divergence structure in Eq. (3.21) — that all local divergences have the form m^2_alpha beta Delta_epsilon — is assumed rather than proved in this paper. The stress-test note asks whether additional subdivergences could appear for the two-field case. My reading is that this is a standard result for the 2PI-Hartree approximation in multicomponent scalar models, and the paper cites the relevant work (Fejos, Patkos, Szep). But because the paper's novelty rests on the mixing, a referee should ask for a more explicit argument or a precise citation that covers the two-field case with off-diagonal propagators. Second, there is no code or data release; the numerical results (gap equations, parameter matching, bounce actions) are described but not independently checkable. For a methods paper, that is a significant omission. Third, the two-step benchmark uses different lambda_phi_chi for the Hartree and Parwani runs because the parameter windows do not overlap, so the comparison is not one-to-one. The authors are upfront about this, but it weakens the quantitative conclusions. Finally, 'valid for all temperatures' is a bit strong for a Hartree truncation; within the truncation it's fine, but the paper should hedge.\n\nIf you are in the phase-transition or thermal-field-theory business, this is worth your time. I would send it to a competent referee — the technical core is substantial and the flaws are fixable in revision. I would not desk-reject it.","headline":"Solid 2PI-Hartree renormalization for two-field mixing with a careful physical parameter connection; the main open questions are about the assumed divergence structure and the lack of reproducible numerics.","tokens_in":42847,"tokens_out":3532,"would_cite":true,"duration_ms":32492,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper derives a thermal effective potential from renormalized 2PI actions that is valid at all temperatures and shows the choice of resummation can change predicted phase-transition strengths and gravitational wave spectra by orders…","keywords":["thermal resummation","two-particle irreducible effective action","Hartree approximation","finite-temperature effective potential","cosmological phase transitions","gravitational waves","renormalization","gap equation"],"falsifier":"Compute $V_{\\rm HT}(\\varphi,\\chi;T)$ for a fixed physical parameter set at two widely separated renormalization scales, say $Q=100$ GeV and $Q=10^6$ GeV; if the potential changes at any field value or temperature, the claimed exact scale invariance is false. A second, sharper test is to include the sunset diagram in the 2PI expansion and check whether a divergence not of the form $m^2_{\\alpha\\beta}\\Delta_\\epsilon$ appears; such a divergence would show that the Hartree-level cancellation does not generalize.","tokens_in":41825,"feed_emoji":"🌌","tokens_out":8521,"duration_ms":81741,"temperature":0.7,"pith_summary":"The paper sets out to give cosmological first-order phase transitions a thermal resummation that does not rely on the high-temperature expansion. Working with two mixing real scalar fields, it renormalizes the two-particle-irreducible (2PI) effective action in the Hartree approximation and connects the auxiliary renormalized parameters to physical masses and couplings. The result is the Hartree-resummed finite-temperature effective potential of Eq. (5.3), which is scale invariant and valid from high to low temperature. The authors then show that this potential can predict transition strengths, and hence gravitational wave amplitudes, that differ by orders of magnitude from the standard ring and Parwani resummations, in either direction depending on the benchmark.","feed_headline":"Thermal resummation choice changes gravitational wave forecasts","feed_subtitle":"Scale-invariant 2PI potential shifts predicted gravitational wave signals by orders of magnitude.","key_machinery":"The load-bearing object is the Hartree-resummed effective potential $V_{\\rm HT}(\\varphi,\\chi;T)$ of Eq. (5.3): a vacuum one-loop-type piece evaluated with thermally corrected gap masses plus the resummed thermal integrals $J$ and $I$. The gap equation (3.22) for the mass matrix $m^2_{\\alpha\\beta}$, solved self-consistently with the local correlation functions, is what makes the resummation consistent at all temperatures. The renormalization rests on the cancellation-of-subdivergences conditions (3.24)-(3.25), which assume the sole divergence is $m^2_{\\alpha\\beta}\\Delta_\\epsilon$ and thereby fix the counterterms. Finally, the connection between auxiliary MS-parameters and physical parameters is established by demanding that the second and fourth derivatives of the potential at the minimum reproduce the $p^2=0$ masses $\\hat m^2_{\\alpha\\beta}$ and zero-momentum couplings $\\hat\\lambda_{\\alpha\\beta}$.","core_discovery":"On the paper's own terms, the central discovery is that a fully renormalized 2PI Hartree treatment supplies a consistent, self-consistent thermal resummation for a multi-field scalar model: thermal masses are obtained as solutions of the gap equation (3.22), the effective potential (5.3) is built from those masses, and no high-temperature approximation enters. The same renormalization procedure that removes ultraviolet divergences in vacuum also removes them at finite temperature, because the only divergent part of the local correlation function is $m^2_{\\alpha\\beta}\\Delta_\\epsilon$ and the cancellation conditions (3.24)-(3.25) make every $n$-point function finite. The potential is exactly scale invariant, and the auxiliary MS-parameters are tied to the physical $p^2=0$ masses and to the zero-momentum four-point functions of the theory. When applied to one- and two-step phase transitions, the Hartree potential typically, but not always, predicts the strongest transitions, and the resulting gravitational wave spectrum can be larger by up to four orders of magnitude than the Parwani prediction.","pith_inferences":["A natural next test is to push the 2PI expansion beyond Hartree: if a next-to-leading-order sunset-type calculation shifts the $p^2=0$ masses by a large amount, the Hartree prediction of phase-transition strength would need revision.","Because the Hartree potential is scale invariant by construction, using it should reduce the renormalization-scale uncertainty that plagues standard perturbative phase-transition predictions; this could be checked directly by repeating the benchmarks at different renormalization scales $Q$.","The same framework could be used to compute bubble-wall velocities and baryogenesis observables in parameter regions where the high-temperature expansion is invalid, exactly where standard schemes are least trustworthy."],"forward_implications":["If $V_{\\rm HT}$ is the correct all-temperature potential, ring-resummed predictions for transition strength and gravitational waves are unreliable whenever $T$ is not much larger than the relevant masses.","Parwani-type resummation captures qualitatively similar physics because it also thermalizes the vacuum part, but its non-self-consistent masses introduce uncontrolled errors near $T\\sim m$.","A stronger, or weaker, Hartree transition directly translates into a larger, or smaller, gravitational wave amplitude, by up to four orders of magnitude in the benchmarks studied, so phenomenological scans should include this scheme.","The renormalization and parameter-matching procedure extends to models with more scalar fields, such as singlet extensions of the Standard Model, making the method applicable to realistic electroweak phase transitions."],"supporting_citations":[{"why":"Introduce the 2PI effective action as the resummation framework on which the whole calculation is built.","marker":"[52, 53]"},{"why":"Establishes that 2PI actions can be renormalized with local counterterms despite the mixing of loop orders.","marker":"[57]"},{"why":"Supply the cancellation-of-subdivergences method and the MS renormalization conditions used here.","marker":"[64, 65, 82]"},{"why":"Shows how to connect 2PI MS-parameters to physical parameters in a single-field model, the procedure extended here to two mixing scalars.","marker":"[56]"},{"why":"Provides the single-field example in which wave-function renormalization factors differ from unity and gives the renormalized equations of motion.","marker":"[76]"},{"why":"Define the ring resummation scheme against which the Hartree potential is compared.","marker":"[24, 25]"},{"why":"Defines the Parwani resummation scheme against which the Hartree potential is compared.","marker":"[26]"}],"fun_headline_variants":["2PI thermal resummation changes gravitational wave forecasts","2PI resummation shifts gravitational wave spectra by orders","Renormalized 2PI Hartree potential changes gravitational wave signal","Consistent 2PI resummation avoids high-T approximation in GW spectra","2PI thermal masses alter phase transitions and gravitational wave forecasts"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that every ultraviolet divergence in the local propagators has exactly the form $m^2_{\\alpha\\beta}\\Delta_\\epsilon$; if additional divergence structures exist beyond the Hartree level, the cancellation conditions (3.23)-(3.25) would not render the potential finite and the central construction collapses.","fun_headline_variants_meta":{"raw":{"variants":["2PI thermal resummation changes gravitational wave forecasts","2PI resummation shifts gravitational wave spectra by orders","Renormalized 2PI Hartree potential changes gravitational wave signal","Consistent 2PI resummation avoids high-T approximation in GW spectra","2PI thermal masses alter phase transitions and gravitational wave forecasts"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000841,"raw_usage":{"total_tokens":3641,"prompt_tokens":899,"completion_tokens":2742,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":515,"completion_tokens_details":{"reasoning_tokens":2656}},"tokens_in":515,"tokens_out":2742,"duration_ms":23095,"temperature":1.0,"reasoning_tokens":2656,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T00:35:05.062667+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute $V_{\\rm HT}(\\varphi,\\chi;T)$ for a fixed physical parameter set at two widely separated renormalization scales, say $Q=100$ GeV and $Q=10^6$ GeV; if the potential changes at any field value or temperature, the claimed exact scale invariance is false. A second, sharper test is to include the sunset diagram in the 2PI expansion and check whether a divergence not of the form $m^2_{\\alpha\\beta}\\Delta_\\epsilon$ appears; such a divergence would show that the Hartree-level cancellation does not generalize.","supporting_citations":[{"cited_title":"Nonperturbative renormalization for 2PI effective action techniques","cited_arxiv_id":"hep-ph/0503240","evidence_quote":"Establishes that 2PI actions can be renormalized with local counterterms despite the mixing of loop orders."},{"cited_title":"Renormalized equations of motions for scalars and fermions in the 2PI formalism","cited_arxiv_id":"2307.14983","evidence_quote":"Provides the single-field example in which wave-function renormalization factors differ from unity and gives the renormalized equations of motion."}],"review_version":1}