{"id":"c5f3238e-1a70-4b62-9164-50b17604451d","arxiv_id":"2504.13725","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"The conventional kinetic equation of motion for bubble walls is incomplete, missing a condensate self-energy term that produces additional friction from particle production, mixing, and transition radiation.","lead":"This paper derives the bubble wall equation of motion from nonequilibrium quantum field theory and finds that the commonly used kinetic equation misses a non-local self-energy term. The missing term captures friction from particle production by the wall's background field, which matters for gravitational wave and dark matter predictions.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The leading-order truncation of the 2PI effective action in Eq. (21) is asserted without a power-counting proof; if neglected higher-loop or VIA-subleading terms carry an extra log(γw) with only coupling suppression, the quantitative friction estimates in Secs. 3.3–4 are not controlled.","rationale":"The central derivation is internally consistent: the non-local self-energy term follows from the exact operator EoM and the 2PI effective action, and the recovery of known kick-approach results (Bödeker–Moore, Azatov–Vanvlasselaer, Ai's logarithmic friction) provides external anchoring. The concern raised here is about quantitative control, not about the existence of the new term. The paper asserts—without proof—that the two-loop diagrams with two condensate insertions are the leading dissipative contribution and that subleading terms only give higher-order corrections. This is load-bearing because the quantitative friction formulas in Sections 3.3 and 4 depend on the truncation, and the VIA expansion used for mixing and transition radiation has an uncontrolled expansion parameter for realistic wall amplitudes. The paper's own Appendix A.3 demonstrates that a naively organized cumulant expansion misses an O(1) constant, which had to be fixed by hand in Eq. (A190); this is a concrete symptom that the truncation is not fully systematic. A three-loop check would settle whether large log γw enhancements invalidate the leading-order estimate. Given that the conceptual claim is well-supported but the quantitative accuracy is not fully established, the CONDITIONAL verdict from the reader remains appropriate; no verdict change is needed.","tokens_in":49316,"tokens_out":18707,"duration_ms":183029,"concrete_test":"In the toy model (4), compute the O(g^4) three-loop contribution to the retarded condensate self-energy Π^R_φ—e.g., the two-loop sunset diagram with a one-loop χ self-energy insertion—and evaluate its contribution to P_vertex via Eq. (51) in the ultrarelativistic limit. If the result is g² times a constant times the leading term, the truncation is safe; if it is g² log(γw) times the leading term, then at γw values with g² log γw ~ 1 the leading-order friction estimate is uncontrolled and the quantitative claims in Secs. 3.3–4 need revision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim is that the non-local self-energy term in Eq. (30) captures the leading dissipative friction from condensate-dependent vertices. This rests on the assertion that the two-loop 2PI diagrams of Eq. (21) are the leading contributions and that sub-leading local terms only give higher-order corrections. No systematic power counting is provided: the expansion parameter for the VIA in Sec. 4 (e.g., κφ/mχ² for mixing, g2φ/mA for transition radiation) is not bounded for realistic wall amplitudes v_b/T ~ O(1–10), and the loop expansion can mix with large log γw enhancements. Concretely, the pair-production pressure in Eq. (65) grows like log(γw T/(2π Lw mχ²)); a next-order diagram with the same log but suppressed only by g² would be comparable when g² log γw ~ O(1), which can occur for large γw. The paper's own Appendix A.3 shows the cumulant expansion misses an O(1) constant, requiring a fit (Eq. A190), indicating the truncation is not fully under control. If higher-order terms modify the γw-scaling or the coefficient, the central claim that all kick-approach friction is captured quantitatively is weakened, even though the existence of the new term is likely correct.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper derives the bubble wall equation of motion from the closed-time-path (CTP) formalism and the two-particle-irreducible (2PI) effective action. The central claim is that the conventional kinetic equation for the condensate, Eq. (1a), is incomplete: it misses a non-local condensate self-energy term of the form ∫d⁴x' Π^R_φ(x,x')φ(x'), shown in Eq. (30). This term arises from two-loop 2PI vacuum diagrams with two condensate insertions, Eq. (21), and describes particle production from condensate-dependent vertices such as φϕχ². The paper identifies this as a new source of dissipative friction, P_vertex in Eq. (37), and shows how the kinetic approach can also accommodate, via the VEV insertion approximation (VIA), the previously kick-approach-only processes of 1-to-1 mixing and 1-to-2 transition radiation. It closes with a localization procedure, Eq. (134), that makes the non-local term usable in numerical computations. The derivation of the exact operator equation, Eq. (8), and the 2PI-based identification of the missing term are internally consistent, but several quantitative steps rely on control assumptions that are not fully established.","tokens_in":49588,"tokens_out":5235,"duration_ms":54174,"significance":"If the central claim holds, this paper fills a genuine conceptual gap between the kinetic and kick approaches to bubble wall dynamics: it provides a first-principles derivation that reproduces the existing kick-approach friction formulas for particle production, mixing, and transition radiation from a single nonequilibrium field-theoretic framework. The strength of the paper is the clean derivation of Eq. (30), the explicit cutting-rule interpretation of Im Π^R_φ, and the verification in Appendix B that the VIA expansions satisfy the Dyson-Schwinger equations. The paper also delivers a locally approximated equation of motion that can be implemented in existing codes such as WallGo. However, the quantitative content is uneven: the pair-production formula has an undetermined O(1) constant fixed only by fitting to numerics, the mixing result is stated to match the literature only up to a factor of two, and the transition-radiation analysis yields a scaling law rather than a controlled coefficient. These issues do not undermine the existence of the new term, but they do affect the claim that the kinetic approach quantitatively reproduces all kick-approach friction processes.","major_comments":[{"comment":"The reduction of the condensate self-energy to the two-loop 2PI vacuum diagrams with two condensate insertions is asserted rather than derived. No power-counting argument is given for neglecting higher-loop diagrams or for the VIA expansion parameters κφ/m_χ² and g₂φ/m_A used in Section 4, and these parameters are not bounded for realistic wall amplitudes. Because the pressure formulas in Eqs. (65), (84), and (128) rely on this truncation, and because large log γ_w enhancements could in principle compensate coupling suppression, the quantitative estimates are not controlled. Please provide a systematic power-counting estimate, or state explicitly the regime in which the truncation is valid and soften the quantitative claims accordingly.","section":"Section 3.1, Eq. (21)"},{"comment":"The paper itself shows that the cumulant expansion of Eq. (A170) does not fix an O(1) constant: Eq. (A185) contains the factor √e σ in the logarithm, and Eq. (A190) drops this factor by hand to match numerics. The final formula quoted in the main text, Eq. (65), therefore contains a fitted constant. This is a legitimate diagnostic, but it means that the prefactor in Eq. (65) has an undetermined order-one uncertainty. Since Eq. (65) is used in the phenomenological comparison of Section 3.3, the uncertainty should be stated explicitly and propagated into the dominance condition g²v_b² ≳ 32π²Δm_φ².","section":"Appendix A.3, Eq. (A190)"},{"comment":"The text says that Eq. (84) 'matches' the result of Ref. [92] but 'up to a factor of two'. A factor of two is not a benign normalization difference for a quantitative recovery claim, and the source of this mismatch is not identified. Please determine whether it originates in the definition of κ versus 2B, in the Fourier transform convention, in the treatment of both chiralities/spins, or in an actual discrepancy, and either resolve the factor or explicitly state that the kinetic approach reproduces the mixing pressure only up to an unexplained factor of two.","section":"Section 4.1.1, Eq. (84)"},{"comment":"The transition-radiation result is derived only at the level of a scaling law, P_TR ∝ γ_w T³ v_b, under several uncontrolled approximations: the soft limit x≪1, the replacement |φ̃₂(Δp_z)|² by v_b⁴/(Δp_z)², and the use of a thermal screening mass as the IR cutoff. The paper also acknowledges differences from Ref. [89], including a factor of 1/2 and the replacement of 1/k⊥⁴ by 1/(k⊥²+m_A²)². This is sufficient to show that the kinetic approach contains the transition-radiation process, but it does not quantitatively reproduce the existing kick-approach result. The text should distinguish 'captured in principle' from 'quantitatively reproduced', otherwise the central claim is overstated.","section":"Section 4.2.2, Eqs. (123)-(128)"}],"minor_comments":[{"comment":"The notation Δ² = 4m_χ² − m_φ² is introduced but the radicand in Eq. (63) and the derivation in Appendix A.1 would be clearer if the dimensionless variables x, y, z defined in Eq. (A153) were used consistently in the main text as well.","section":"Eq. (63) and Appendix A.1"},{"comment":"The right panel of Fig. 6, showing the ratio of numerical to analytic pressure, uses an inset with small axis labels; the figure should be reproduced at larger size or with a separate panel so that the claimed factor-of-two discrepancy is legible.","section":"Fig. 6"},{"comment":"The phenomenological paragraph discusses the electroweak phase transition and processes such as h→hh, h→WW, and Z→hZ, but it does not give the corresponding formula or estimate for these light-particle cases. A brief indication of which terms in Eq. (64) would dominate, or a reference to future numerical work, would help the reader assess the phenomenological relevance.","section":"Section 3.3, after Eq. (64)"},{"comment":"The localization procedure identifies the last term of Eq. (132) as 'the standard form for a friction term', but the sign convention relative to the frictional pressure P_vertex in Eq. (37) is not discussed; a short comment on the sign would avoid confusion when the local equation is used in practice.","section":"Section 5, Eq. (132)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript makes a substantive and mostly convincing conceptual advance, and the simultaneous posting with Ref. [106] is acknowledged appropriately. The main risk is not circularity but quantitative control: the paper's own Appendix A.3 demonstrates that the analytic pair-production formula requires a fitted constant, and the factor-of-two mismatch in the mixing channel is unexplained. I recommend major revision, with the expectation that the authors either prove or explicitly bound the neglected terms in the 2PI/VIA truncations, and that they clearly separate conceptual recovery from quantitative reproduction in their claims. Once those points are addressed, the paper would be a strong candidate for acceptance."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper is likely right about the conceptual gap. The commonly used kinetic EoM (1a) misses a nonlocal self-energy term coming from condensate-dependent vertices, and the derivation from the exact operator EoM plus the 2PI effective action is genuine. The kick-approach pair production, mixing, and transition radiation results being recovered from the VIA is real external anchoring. This is worth refereeing.\n\nWhat is new: Eq. (30) adds the nonlocal self-energy term, derived rather than postulated, and the exact operator EoM in Eq. (8) plus the 2PI derivation are internally consistent. Appendix B verifying Dyson-Schwinger consistency is useful. The authors also document honestly in Appendix A.3 that the cumulant expansion cannot fix an O(1) constant and fit it to numerics; I take that as credit, not a hidden flaw.\n\nSoft spots, in order of importance. First, the two-loop truncation of Gamma_2 in Eq. (21) is asserted without a power-counting proof. The stress-test note is on target: the VIA expansion parameters like kappa*phi/m_chi^2 and g2*phi/m_A are not bounded for realistic v_b/T ~ 1-10, and higher orders could carry log(gamma_w) enhancements with only coupling suppression. The pair production pressure itself grows like log(gamma_w T/(2 pi L_w m_chi^2)), so g^2 log(gamma_w) ~ O(1) is possible. This leaves the quantitative size of the new friction less controlled than the existence of the term.\n\nSecond, the fitted O(1) constant is not a small blemish. Figure 6 shows the analytic formula underestimates the numerical pair production pressure by a factor of two even at gamma_w ~ 250, after the fit. So the asymptotic formulas in Sec. 3.3 are order-of-magnitude estimates, not predictions.\n\nThird, the transition radiation result is parametric: Eq. (128) P_TR ~ gamma_w T^3 v_b is a scaling statement with an unquantified IR cutoff, acknowledged in the text. That is fine for the paper's conceptual claim but limits direct phenomenological use.\n\nThe central claim, that the kinetic approach can capture all four kick processes, holds up at the level of the derivation. What is not yet established is whether the new friction is numerically important in realistic models. A revision should add a systematic power-counting argument or at least a numerical check of higher-order diagrams in a simple model, and state the parametric uncertainty in the transition radiation pressure.\n\nThis paper deserves a serious referee despite the soft spots. It directly addresses a contradiction with Ref. [105], gives a first-principles derivation, and reproduces external benchmarks. I would send it to review with a request for major revision, not desk reject.","headline":"Likely right about the missing nonlocal friction term in the kinetic bubble-wall EoM, but the quantitative friction estimates are under-controlled and need revision before the paper is used for phenomenology.","tokens_in":50133,"tokens_out":1927,"would_cite":true,"duration_ms":17667,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The standard bubble-wall equation of motion is incomplete, and the missing term adds friction from particle production.","keywords":["bubble wall velocity","first-order phase transitions","nonequilibrium quantum field theory","closed-time-path formalism","two-particle-irreducible effective action","condensate self-energy","bubble wall friction","relativistic bubble walls"],"falsifier":"Compute the next-order, three-loop, three-condensate-insertion contribution to the condensate self-energy for the scalar model with a $\\varphi\\phi\\chi^2$ vertex and evaluate its friction pressure in the ultrarelativistic limit; if that pressure grows with $\\gamma_w$ as fast as the two-loop result, the truncation is not the leading dissipative effect. Equivalently, solve the full coupled equations for the one- and two-point functions numerically for a single planar wall and compare the terminal velocity with the localised equation's prediction.","tokens_in":49074,"feed_emoji":"🫧","tokens_out":13409,"duration_ms":108571,"temperature":0.7,"pith_summary":"This paper derives the coupled motion of a first-order phase-transition bubble wall and the plasma from nonequilibrium quantum field theory, using the closed-time-path (real-time) formalism and the two-particle-irreducible effective action, a functional whose variation gives coupled equations for the wall and the particle propagators. It claims that the standard equation of motion of the wall, in which the plasma acts only through field-dependent particle masses, is incomplete. The missing piece is a non-local condensate self-energy term, $\\int d^4x'\\, \\Pi^R_\\varphi(x,x')\\varphi(x')$, generated by condensate-particle vertices such as $\\varphi\\phi\\chi^2$. This term describes particle production from the wall itself and adds a friction pressure that earlier kinetic-approach studies overlooked. If the derivation is right, the kinetic approach can reproduce every process that the microscopic kick picture counts, namely mass-gain transmission, 1-to-1 mixing, 1-to-2 transition radiation, and 1-to-2 pair production, and it does so for walls at any speed, not only ultrarelativistic ones.","feed_headline":"One missing term changes bubble-wall friction","feed_subtitle":"Adding the nonlocal condensate self-energy recovers pair production, mixing, and transition radiation in the kinetic approach.","key_machinery":"The central object is the retarded condensate self-energy, $\\Pi^R_\\varphi(x,x') \\equiv \\Pi^{++}_\\varphi(x,x') - \\Pi^{+-}_\\varphi(x,x')$, defined by varying the two-particle-irreducible effective action twice with respect to the background field. In a planar-wall frame with translational symmetry along the wall and in time, its Fourier transform with respect to $z-z'$ enters the vertex friction as an integral of $q_z |\\tilde\\varphi(q_z)|^2 \\operatorname{Im}\\tilde\\Pi^R_\\varphi(q_z)$. The imaginary part behaves like a collision term, so the wall equation gains dissipative dynamics through the same object that describes particle production. The VEV insertion approximation, in which $\\varphi$-dependent masses are treated as insertions of the background field, converts mixing and transition-radiation effects into induced self-energies of the same non-local form. A phase-space transform and gradient expansion then localise the term into a field-dependent mass correction $\\Delta m^2_{\\Pi_\\varphi}\\varphi$ and a damping term $-(\\partial_\\mu\\varphi)\\lim_{q\\to 0}\\operatorname{Im}\\Pi^R_\\varphi(q,x)/q_\\mu$.","core_discovery":"The paper shows that the retarded condensate self-energy, $\\Pi^R_\\varphi = \\Pi^{++}_\\varphi - \\Pi^{+-}_\\varphi$, originating from two-loop two-particle-irreducible vacuum diagrams with two insertions of the background field, belongs in the wall equation of motion. Its imaginary part is a collision term: cutting the diagrams gives particle-production processes in which a plasma particle extracts momentum from the wall. Integrating the self-energy against the wall profile yields a vertex-induced friction $P_{\\rm vertex}$ distinct from the familiar mass-change friction $P_{\\rm mass}$. In the ultrarelativistic limit the formalism recovers the known logarithmic growth of pair-production pressure with $\\gamma_w$, and with a VEV insertion approximation, which expands propagators in powers of the background field, it also reproduces the pressure from scalar and fermion mixing and from gauge-boson transition radiation, including the linear growth with $\\gamma_w$ that prevents runaway walls. A gradient expansion localises the non-local term into a mass correction plus a standard dissipative damping term, giving an equation of motion suitable for numerical wall-velocity computations.","pith_inferences":["The same condensate self-energy should also appear in the collision terms of the Boltzmann equations for the fluctuations, not only in the wall equation; the paper notes this possibility but leaves it for future work.","If the two-loop truncation is robust, the framework implies a testable interpolation for friction at intermediate wall velocities, where the kick approach has no prediction and local-thermal-equilibrium estimates break down.","Comparing the localised equation's terminal velocity against a direct numerical solution of the full coupled equations for the one- and two-point functions in a simple scalar model would quantify the error from dropping subleading local terms.","For electroweak-scale transitions, the same formalism should generate new collision terms involving the Higgs self-interactions $h \\to hh$ and $h \\to W^+W^-,ZZ$; the paper lists these processes as natural next applications."],"forward_implications":["Conventional wall-velocity calculations that use only mass-dependent friction undercount the stopping force whenever the theory has condensate-particle vertices such as $\\varphi\\phi\\chi^2$.","The kinetic approach can recover the kick-approach results, namely pair production with logarithmic growth in $\\gamma_w$, mixing friction, and transition-radiation friction linear in $\\gamma_w$, without assuming ballistic particle motion.","In dark-sector phase transitions where the light-field mass change is loop-suppressed, the new pair-production pressure can exceed the mass-gain pressure and dominate the wall's terminal velocity.","The localised equation of motion provides a concrete dissipative damping term that can be inserted into numerical bubble-wall calculations, extending them to non-ultrarelativistic walls.","Transition radiation tends to prevent runaway bubble walls in gauged phase transitions because its pressure grows with the Lorentz factor."],"supporting_citations":[{"why":"Defines the two-particle-irreducible effective action whose variation yields the condensate equation of motion.","marker":"[102]"},{"why":"Shows how the two-particle-irreducible equations reduce to Boltzmann equations, grounding the kinetic approach in the same formalism.","marker":"[107]"},{"why":"Gives the standard ultrarelativistic mass-gain friction pressure that the new vertex term supplements.","marker":"[88]"},{"why":"Derives the transition-radiation friction in the kick approach that Section 4.2 reproduces from the closed-time-path formalism.","marker":"[89]"},{"why":"Provides the 1-to-1 mixing friction in the kick approach that Section 4.1 recovers via the VEV insertion approximation.","marker":"[92]"},{"why":"Establishes the logarithmically growing pair-production friction that Section 3.3 reproduces and generalises.","marker":"[95]"},{"why":"Identifies condensate-particle pair production as a production process for heavy dark matter, here shown to also contribute to wall friction.","marker":"[30]"},{"why":"Argues that the conventional equation of motion is complete; the paper's central claim directly addresses this opposing argument.","marker":"[105]"}],"fun_headline_variants":["Missing self-energy term rewrites bubble wall friction","Extra friction term from particle production in bubble walls","Self-energy collision term adds new bubble wall damping","Missing vertex term fixes bubble wall friction","Bubble wall friction gets new particle production term"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation depends on assuming that the simplest set of loop diagrams considered gives the dominant friction and that all neglected diagrams are only small corrections; the paper states this but does not prove it, and the quantitative friction results rely on it.","fun_headline_variants_meta":{"raw":{"variants":["Missing self-energy term rewrites bubble wall friction","Extra friction term from particle production in bubble walls","Self-energy collision term adds new bubble wall damping","Missing vertex term fixes bubble wall friction","Bubble wall friction gets new particle production term"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000573,"raw_usage":{"total_tokens":2701,"prompt_tokens":934,"completion_tokens":1767,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":550,"completion_tokens_details":{"reasoning_tokens":1698}},"tokens_in":550,"tokens_out":1767,"duration_ms":10702,"temperature":1.0,"reasoning_tokens":1698,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T12:01:23.248820+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the next-order, three-loop, three-condensate-insertion contribution to the condensate self-energy for the scalar model with a $\\varphi\\phi\\chi^2$ vertex and evaluate its friction pressure in the ultrarelativistic limit; if that pressure grows with $\\gamma_w$ as fast as the two-loop result, the truncation is not the leading dissipative effect. Equivalently, solve the full coupled equations for the one- and two-point functions numerically for a single planar wall and compare the terminal velocity with the localised equation's prediction.","supporting_citations":[],"review_version":1}