arxiv: 2604.11058 · v1 · submitted 2026-04-13 · ⚛️ nucl-th · hep-ph· nucl-ex

Recognition: unknown

Nonlinear response of flow harmonics in Gubser flow with participant-reaction planes mismatch

Xiang Ren , Jin-Yu Hu , Hao-jie Xu , Shi Pu

Authors on Pith no claims yet

Pith reviewed 2026-05-10 16:16 UTC · model grok-4.3

classification ⚛️ nucl-th hep-phnucl-ex

keywords nonlinear responseflow harmonicsGubser floweccentricitiesparticipant plane mismatchheavy ion collisionsanalytic relationshydrodynamics

0 comments

The pith

A mismatch between participant and reaction planes in Gubser flow can change both the strength and sign of nonlinear response coefficients for flow harmonics v2 and v4.

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

The paper extends perturbative Gubser flow solutions to derive analytic nonlinear response relations that connect initial eccentricities e2 and e4 to observed flow harmonics v2 and v4. It shows that including the mismatch between participant and reaction planes adds an angle-dependent factor to the conventional nonlinear response coefficients. This factor, usually treated as random noise from fluctuations, can instead systematically modify the effective coefficients and even reverse their sign depending on the initial nuclear configuration. A sympathetic reader would care because this provides a concrete analytical handle on how initial geometry influences collective flow in heavy-ion collisions beyond statistical averaging.

Core claim

By extending the perturbative solutions of Gubser flow, analytic nonlinear response relations are derived connecting e2, e4 to v2, v4, reproducing v4/v2^2 to 1/2 at large pT, and the participant-reaction plane mismatch introduces an additional factor from the participant-plane angles that modifies the strength and can even change the sign of the effective nonlinear response coefficient.

What carries the argument

Perturbative extension of Gubser flow solutions incorporating participant-reaction plane mismatch angles in the nonlinear response relations.

If this is right

The ratio v4/v2^2 approaches 1/2 in the large transverse momentum limit.
The conventional nonlinear response coefficients are multiplied by an angle-dependent factor from the mismatch.
This factor makes the effective response sensitive to the specific initial configuration of the colliding nuclei.
The mismatch is not merely statistical noise but can have deterministic effects on observed flows.

Where Pith is reading between the lines

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

This suggests that experimental analyses of flow harmonics should account for event-by-event plane alignments rather than averaging them out.
Similar mismatch effects might appear in other hydrodynamic models or at different collision energies.
Future work could test these analytic relations against full numerical simulations of viscous hydrodynamics.

Load-bearing premise

The perturbative extension of Gubser flow solutions remains valid for the eccentricities and mismatch angles considered in the derivation.

What would settle it

Observation of the effective nonlinear response coefficient v4/v2^2 deviating from the predicted angle-dependent modification in collisions with controlled or measured participant plane orientations.

Figures

Figures reproduced from arXiv: 2604.11058 by Hao-jie Xu, Jin-Yu Hu, Shi Pu, Xiang Ren.

**Figure 2.** Figure 2: FIG. 2. The ratio [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗

**Figure 3.** Figure 3: FIG. 3. Angle [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗

**Figure 4.** Figure 4: FIG. 4. Illustration of Eq. ( [PITH_FULL_IMAGE:figures/full_fig_p010_4.png] view at source ↗

read the original abstract

We investigate the nonlinear response of flow harmonics $v_2,v_4$ to initial-state eccentricities $\epsilon_2,\epsilon_4$ within the Gubser-flow framework. By extending the perturbative solutions of Gubser flow, we derive analytic nonlinear response relations connecting the eccentricities $\epsilon_2,\epsilon_4$ to the flow harmonics $v_2,v_4$. Our results reproduce the well-known result $v_4/v_2^2 \to 1/2$ in large transverse momentum $p_T$ limit. Furthermore, we study the effects of a mismatch between the participant and reaction planes. We find that the conventional nonlinear response coefficients acquire an additional factor determined by the participant-plane angles, which is often approximated as statistical noise driven by event-by-event fluctuations. This factor can modify both the strength but even the sign of the effective nonlinear response coefficient, making it sensitive to the initial configuration of the colliding nuclei. Our study provides new analytical insight into the origin of collective phenomena in relativistic heavy-ion collisions.

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.

Desk Editor's Note private letter to a colleague

The paper adds an explicit participant-reaction plane mismatch factor to analytic nonlinear response coefficients in Gubser flow, which can flip their sign, while recovering the standard v4/v2² limit.

read the letter

This paper derives analytic nonlinear response relations for v2 and v4 from initial eccentricities inside the Gubser flow framework, with the main addition being a trigonometric factor from the fixed mismatch angle between participant and reaction planes. That factor multiplies the usual coefficients and can change both magnitude and sign depending on the initial geometry. The authors treat the mismatch as a deterministic input rather than averaging it away as event-by-event noise, which is the concrete step beyond prior work on Gubser perturbations. They also recover the known high-pT limit where v4 over v2 squared approaches one half, which serves as a useful check that the extension is consistent with established results. The derivation begins from the hydrodynamic equations in Gubser coordinates with the angle difference inserted directly, so the algebra follows from the ideal, conformal, boost-invariant assumptions without extra fitting. The calculation is therefore transparent within its stated limits. The main limitation is that the second-order truncation is not tested for convergence when the mismatch angle makes the effective nonlinear source term order-one. In that regime the omitted higher-order contributions could be sizable, and the reported sign flip might shift once those terms are kept. Gubser flow also fixes exact boost invariance and a particular radial profile, so the result is tied to those idealizations and will likely change once viscosity or longitudinal fluctuations are restored. This work is aimed at theorists who use analytic hydrodynamic solutions to interpret flow data in heavy-ion collisions. It gives a compact way to see how specific initial-state geometry details enter the nonlinear mixing without needing full numerical runs. The paper engages the standard literature on nonlinear responses and Gubser solutions in a direct way. I would send it for peer review. The analytic extension is worth a detailed check by people familiar with these methods, even if the final applicability remains narrow.

Referee Report

1 major / 1 minor

Summary. The paper claims to derive analytic nonlinear response relations connecting initial eccentricities ε₂, ε₄ to flow harmonics v₂, v₄ by extending perturbative Gubser flow solutions. It reproduces the known high-p_T limit v₄/v₂² → 1/2 and shows that a mismatch between participant and reaction planes introduces an additional trigonometric factor that can modify both the magnitude and sign of the conventional nonlinear response coefficients.

Significance. If the derivations hold, this work supplies parameter-free analytic expressions for nonlinear hydrodynamic response within an exactly solvable boost-invariant conformal model, which is a clear strength for benchmarking. The reproduction of the standard v₄/v₂² limit serves as an important consistency check, while the mismatch-angle factor offers a deterministic mechanism that could explain event-by-event variations in observed nonlinear coefficients beyond pure statistical fluctuations.

major comments (1)

[§4] §4 (nonlinear response with plane mismatch): the derivation of the additional angle-dependent factor multiplying the nonlinear coefficients assumes the second-order perturbative truncation remains accurate even when the mismatch angle renders the effective azimuthal source non-small. No estimate of the size of omitted O(ε³) terms or viscous corrections is supplied for the regime where the sign change is predicted; this is load-bearing for the central claim that the factor can reverse the sign of the effective response.

minor comments (1)

[Abstract] Abstract: the phrasing 'modify both the strength but even the sign' is grammatically awkward and should be revised to 'modify both the strength and even the sign' for clarity.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive feedback. We address the single major comment point by point below.

read point-by-point responses

Referee: [§4] §4 (nonlinear response with plane mismatch): the derivation of the additional angle-dependent factor multiplying the nonlinear coefficients assumes the second-order perturbative truncation remains accurate even when the mismatch angle renders the effective azimuthal source non-small. No estimate of the size of omitted O(ε³) terms or viscous corrections is supplied for the regime where the sign change is predicted; this is load-bearing for the central claim that the factor can reverse the sign of the effective response.

Authors: We agree that the validity of the second-order truncation is central and that explicit error estimates would strengthen the presentation. In the derivation the eccentricities ε₂ and ε₄ remain the small expansion parameters; the participant-reaction-plane mismatch angle ψ enters only through geometric projection factors (e.g., cos(4ψ)) that modulate the coefficient but do not alter the magnitude of the source term itself. Consequently the size of the omitted O(ε³) contributions continues to be controlled by |ε| ≪ 1 independently of ψ. Because the underlying Gubser solution is conformal and ideal, viscous corrections lie outside the present analytic framework. We will add a short paragraph in the revised §4 that (i) reiterates the small-ε assumption, (ii) notes that the sign-reversal feature is a prediction within this controlled approximation, and (iii) suggests numerical hydrodynamic tests for regimes where |ε| is no longer small. This constitutes a partial revision that directly responds to the referee’s concern without requiring new third-order calculations. revision: partial

Circularity Check

0 steps flagged

Derivation from Gubser hydrodynamic equations is self-contained with no circular reduction

full rationale

The paper derives the nonlinear response relations by extending the known perturbative Gubser flow solutions (ideal, boost-invariant, conformal hydrodynamics) to second order in eccentricities ε2, ε4 while treating participant-reaction plane mismatch angles as fixed external parameters. The resulting analytic expressions for v2, v4 and the modified nonlinear coefficients follow directly from the hydrodynamic equations in Gubser coordinates; the reproduction of the known large-pT limit v4/v2² → 1/2 functions as an independent consistency check rather than an input. No step reduces by construction to a fitted parameter, self-referential definition, or load-bearing self-citation chain, and the mismatch trigonometric factor emerges from the azimuthal expansion without presupposing the final sign-flip result. The framework remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Limited information available from abstract alone; relies on standard assumptions of ideal hydrodynamics and perturbative expansion around Gubser flow without introducing new free parameters or entities in the summary.

axioms (1)

domain assumption Gubser flow provides an exact analytic solution to ideal hydrodynamic equations under boost-invariant and azimuthally symmetric conditions
Framework invoked for extending perturbative solutions to nonlinear response.

pith-pipeline@v0.9.0 · 5494 in / 1239 out tokens · 28151 ms · 2026-05-10T16:16:54.176409+00:00 · methodology

discussion (0)

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

Investigation of Nonlinear Collective Dynamics in Relativistic Heavy-Ion Collisions Using A Multi-Phase Transport Model
nucl-th 2026-04 unverdicted novelty 4.0

AMPT simulations show the ratio of nonlinear response coefficients between U+U and Au+Au collisions remains stable across evolution stages, isolating initial geometric correlations.

Reference graph

Works this paper leans on

33 extracted references · 30 canonical work pages · cited by 1 Pith paper

[1]

0.5 1. 1.5 2. 2.5 3. 0.5 1. 1.5 4(ψ2-ψ4)/rad ξ*/rad FIG. 3. Angleξ ∗ as a function of 4(ψ 2 −ψ 4). The parameters are chosen to be the same as those in Fig. 1, withϵ 2 = 0.04 and ϵ4 = 0.01. reaction-plane angle Ψ n. To better quantify anisotropic flow, we therefore need to account for the mismatch between ψn and Ψn and its impact on the relation betweenϵ ...
[2]

Experimental and Theoretical Challenges in the Search for the Quark Gluon Plasma: The STAR Collaboration's Critical Assessment of the Evidence from RHIC Collisions

STAR, J. Adamset al., Nucl. Phys. A757, 102 (2005), nucl-ex/0501009

work page Pith review arXiv 2005
[3]

Adcoxet al.(PHENIX), Nucl

PHENIX, K. Adcoxet al., Nucl. Phys. A757, 184 (2005), nucl-ex/0410003

work page arXiv 2005
[4]

Arsene et al

BRAHMS, I. Arseneet al., Nucl. Phys. A757, 1 (2005), nucl-ex/0410020

work page arXiv 2005
[5]

Aamodtet al., JINST3, S08002 (2008)

ALICE, K. Aamodtet al., JINST3, S08002 (2008)

2008
[6]

Chatrchyanet al., Phys

CMS, S. Chatrchyanet al., Phys. Rev. C89, 044906 (2014), 1310.8651

work page arXiv 2014
[7]

Aadet al., Phys

ATLAS, G. Aadet al., Phys. Rev. C90, 024905 (2014), 1403.0489

work page arXiv 2014
[8]

Adareet al., Phys

PHENIX, A. Adareet al., Phys. Rev. C93, 051902 (2016), 1412.1038

work page arXiv 2016
[9]

Alver and G

B. Alver and G. Roland, Phys. Rev. C81, 054905 (2010), 1003.0194, [Erratum: Phys.Rev.C 82, 039903 (2010)]

work page arXiv 2010
[10]

Alveret al., Phys

PHOBOS, B. Alveret al., Phys. Rev. Lett.98, 242302 (2007), nucl-ex/0610037. 16

work page arXiv 2007
[11]

Qiu and U

Z. Qiu and U. W. Heinz, Phys. Rev. C84, 024911 (2011), 1104.0650

work page arXiv 2011
[12]

F. G. Gardim, Y. Hama, and F. Grassi, Prog. Theor. Phys. Suppl.193, 319 (2012), 1110.5658

work page arXiv 2012
[13]

Yan and J.-Y

L. Yan and J.-Y. Ollitrault, Phys. Lett. B744, 82 (2015), 1502.02502

work page arXiv 2015
[14]

Zhao, H.-j

S. Zhao, H.-j. Xu, Y.-X. Liu, and H. Song, Phys. Lett. B839, 137838 (2023), 2204.02387

work page arXiv 2023
[15]

H.-j. Xu, J. Zhao, and F. Wang, Phys. Rev. Lett.132, 262301 (2024), 2402.16550

work page arXiv 2024
[16]

Z. Wang, J. Chen, H.-j. Xu, and J. Zhao, Phys. Rev. C110, 034907 (2024), 2405.09329

work page arXiv 2024
[17]

S. S. Gubser, Phys. Rev. D82, 085027 (2010), 1006.0006

work page arXiv 2010
[18]

S. S. Gubser and A. Yarom, Nucl. Phys. B846, 469 (2011), 1012.1314

work page arXiv 2011
[19]

Hatta and B.-W

Y. Hatta and B.-W. Xiao, Phys. Lett. B736, 180 (2014), 1405.1984

work page arXiv 2014
[20]

3n(1−r 2) g p 1 +g 2 4(g+t) +λ(g)(1 +g 2 −r 2) − −4 r2 ¯H0g−1/3 − 4r2 + (1 +g 2) 1 +g 2 1 + ¯H0g2/3 2 n(1−r 2) 2g+t # , δA2,⊥ = n t2r(1 +g 2) ϵnAn(r) cosnϕ

Y. Hatta, J. Noronha, G. Torrieri, and B.-W. Xiao, Phys. Rev. D90, 074026 (2014), 1407.5952

work page arXiv 2014
[21]

Hatta, A

Y. Hatta, A. Monnai, and B.-W. Xiao, Phys. Rev. D92, 114010 (2015), 1505.04226

work page arXiv 2015
[22]

Hatta, A

Y. Hatta, A. Monnai, and B.-W. Xiao, Nucl. Phys. A947, 155 (2016), 1507.04690

work page arXiv 2016
[23]

CERN Report No., , 2012 (unpublished)

2012
[24]

Borghini and J.-Y

N. Borghini and J.-Y. Ollitrault, Phys. Lett. B642, 227 (2006), nucl-th/0506045

work page arXiv 2006
[25]

Teaney and L

D. Teaney and L. Yan, Phys. Rev. C86, 044908 (2012), 1206.1905

work page arXiv 2012
[26]

Lang and N

C. Lang and N. Borghini, Eur. Phys. J. C74, 2955 (2014), 1312.7763

work page arXiv 2014
[27]

Gombeaud and J.-Y

C. Gombeaud and J.-Y. Ollitrault, Phys. Rev. C81, 014901 (2010), 0907.4664

work page arXiv 2010
[28]

F. G. Gardim, J. Noronha-Hostler, M. Luzum, and F. Grassi, Phys. Rev. C91, 034902 (2015), 1411.2574

work page arXiv 2015
[29]

Luzum, J

M. Luzum, J. Phys. G38, 124026 (2011), 1107.0592

work page arXiv 2011
[30]

Aamodtet al.(ALICE), Phys

ALICE, K. Aamodtet al., Phys. Rev. Lett.107, 032301 (2011), 1105.3865

work page arXiv 2011
[31]

Adareet al., Phys

PHENIX, A. Adareet al., Phys. Rev. Lett.107, 252301 (2011), 1105.3928

work page arXiv 2011
[32]

L. V. Bravinaet al., Phys. Rev. C89, 024909 (2014), 1311.0747

work page arXiv 2014
[33]

L. V. Bravinaet al., Eur. Phys. J. C74, 2807 (2014), 1311.7054

work page arXiv 2014