Pith. Sign in

REVIEW

Strongly mixed cosmological collider at unequal sound speeds

PaperPDF from arXiv

data files

ReviewThe objections and verdict, as JSON Paper breakdownThe authors' claims, conditions and measurements, as JSON Signed receiptProof this review has not been altered
X Bluesky LinkedIn Reddit Hacker News
Javier Huenupi, Gonzalo A. Palma, Spyros Sypsas

Cosmology and Nongalactic Astrophysics · Submitted 1 Oct 2026 · 87 references

Pith's verdict: Extends the equal-speed strong-mixing solutions to unequal sound speeds with an exact Heun-function spectrum and explicit hierarchy laws; the reader's keystone worry is mostly defused because identity (A.3) is actually derived, not just numerically checked.

T0 · AI review

How this review was made: An AI referee (deepseek/deepseek-v4-flash) read the paper on October 2, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 6 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv ID2610.02138 (version 1)
Pith IDGKXCHGTR, Pith's permanent record of this paper and its review
CitationGonzalo A. Palma, Javier Huenupi, Spyros Sypsas. "Strongly mixed cosmological collider at unequal sound speeds." arXiv:2610.02138 (2026). https://pith.science/paper/GKXCHGTR
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report

How Pith checks papers → · Full audit of this paper →

The paper

The paper claims

This paper derives the exact late-time curvature power spectrum of a two-field inflationary system at arbitrary quadratic mixing and unequal sound speeds, showing it reduces to a sum of two squared Heun functions whose parameters are…

The reading

This paper claims that the two-field inflationary system with constant quadratic mixing $\lambda$ and unequal sound speeds $c_\varphi$, $c_\sigma$ admits an exact, closed-form late-time curvature power spectrum for any positive entropy mass $\mu^2 > 0$, with no expansion in the mixing or the sound-speed ratio. The spectrum reduces to a sum of two squared Heun functions whose parameters are fixed by the mass, the mixing, and the ratio $r_c = c_\sigma/c_\varphi$, recovering the unmixed result at $\lambda=0$ and the equal-speed gamma-function result as $r_c\to 1$. At fixed mixing and $r_c\ll 1$ the power is enhanced as $r_c^{-2\sqrt{\nu^2+\lambda^2}}$, grows as $\ln^2(1/r_c)$ at vanishing exponent, or oscillates in $\ln r_c$, while for $r_c\gg 1$ the correction to the unmixed spectrum vanishes as $\lambda^2\ln r_c/r_c^2$. The same normalized modes give tree-level squeezed bispectra for three cubic operators, with the heavy-field oscillation frequency set strictly by the final entropy mass. This matters because strong mixing can bring inflationary collider signals within observational reach, and unequal sound speeds change both the primordial power and the collider amplitude and phase.

What carries the argument

The Laplace representation: the coupled field equations are recast through contour integrals over plane waves, reducing the dynamics to a single second-order differential equation in a spectral variable $u$, which for unequal sound speeds is a general Heun equation with four regular singular points at $u = \pm1$ and $u = \pm r_c$. Contour choices select the positive-frequency incoming curvature and entropy waves, with Bunch\textendash Davies conditions and the conserved inner product fixing the two normalized mode functions $Y_1$ and $Y_2$. The exact spectrum follows once the frozen amplitudes are reduced to Heun values at the origin through the integral identity (A.3), which connects the even imaginary-axis spectral solution to the endpoint value at the entropy singularity; the identity is derived in Appendix A.1 and verified numerically to about $2\times10^{-11}$ in relative power at six representative points.

Extended reading notes

Core claim

The central claim is Eq. (1.8): for $\mu^2 > 0$, the exact late-time curvature power spectrum of the two-field system at arbitrary quadratic mixing $\lambda$ and sound-speed ratio $r_c$ is a sum of two squared Heun functions, $\mathrm{HeunG}(1-r_c^{-2}, Q_1; \tfrac14-\tfrac{\nu}{2}, \tfrac14+\tfrac{\nu}{2}, 1, \tfrac12; 1)$ and $\mathrm{HeunG}(1-r_c^2, Q_2; \tfrac54-\tfrac{\nu}{2}, \tfrac54+\tfrac{\nu}{2}, 2, \tfrac12; 1)$, times fixed gamma-function factors, with accessory parameters $Q_1$, $Q_2$ defined in Eq. (1.9). As the mixing vanishes the spectrum reduces continuously to the unmixed result, and as $r_c\to 1$ it tends to the equal-speed gamma-function expression of the earlier exact solutions. The associated hierarchy laws state that at fixed nonzero mixing the small-$r_c$ regime is controlled by the combination $\chi = \sqrt{\nu^2+\lambda^2}$: the entropy contribution grows as $\lambda^2|C(\lambda,\nu)|^2 r_c^{-2\chi}$, grows as $\ln^2(1/r_c)$ at $\chi=0$, and oscillates as a bounded function of $\ln r_c$ for imaginary $\chi$, while for $r_c\gg 1$ the leading correction to the unmixed spectrum is $\lambda^2[2\ln r_c + B(\nu)]/r_c^2$ plus a subleading $r_c^{-3}$ term, a decoupling limit that holds even at large mixing. For the squeezed bispectrum of the three cubic operators the paper obtains coefficients $b_{I,\pm}$, which for a heavy field with $\nu = i\rho$ combine into an oscillation of logarithmic frequency $\rho$ fixed by the final entropy mass; unequal sound speeds modify the amplitude and phase, with the amplitude non-monotonic in $r_c$, but preserve the frequency and the $q^{3/2}$ squeezed scaling.

Load-bearing premise

The entire exact spectrum rests on the integral identity (A.3), which converts the contour integral into a closed-form Heun expression; the paper establishes this identity only through numerical verification at six points, with agreement below $2\times10^{-11}$ in relative power, and an unanticipated branch or contour subtlety would invalidate the central formula. A secondary fragile premise is that the long mode has reached its late-time form during the short-mode interaction, quantified by Eq. (7.27) and roughly $q \ll \min(r_c, 1/r_c)$, which the paper states as the applicability condition of the squeezed bispectrum.

Editorial extensions

If this is right

  • The two-field power spectrum at any mixing and sound-speed ratio is available in closed form, so template building for CMB and large-scale-structure searches need not expand in either parameter.
  • At slow entropy speed ($r_c \ll 1$) the spectrum can be strongly enhanced as $r_c^{-2\sqrt{\nu^2+\lambda^2}}$, grow as $\ln^2(1/r_c)$, or oscillate in $\ln r_c$, depending only on the combination $\nu^2+\lambda^2$.
  • At fast entropy speed ($r_c \gg 1$) the system decouples: the correction to the unmixed spectrum vanishes as $\lambda^2[2\ln r_c + B(\nu)]/r_c^2$ even for order-one mixing.
  • For heavy entropy fields the squeezed bispectrum is an oscillation in $\ln(k_S/k_L)$ with frequency fixed by $\mu/H$; the sound-speed hierarchy changes the amplitude and phase, and can delay the approach to the squeezed regime.

Reading between the lines

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

  • If the integral identity (A.3) admits a closed-form proof, it would constitute a new connection formula for the general Heun equation, a special-function result applicable well beyond cosmology.
  • Because the heavy-field oscillation frequency is robust to $r_c$ while the amplitude and phase vary, inferring the entropy mass from collider detections should marginalize over the sound-speed ratio, which otherwise enters only amplitude and phase.
  • The non-monotonic amplitude with $r_c$ implies that templates fixed at $r_c = 1$ could substantially misestimate collider amplitudes; full-shape surveys would need to scan $r_c$ alongside mass and mixing.
  • The two-sound-horizon picture suggests the $\chi = \sqrt{\nu^2+\lambda^2}$ enhancement is generic to any multi-field setup whose extra field crosses its sound horizon before curvature, so the hierarchy laws may transfer to models with derivative couplings beyond the constant-mixing case.
Original abstract

We derive the canonically normalized modes of a two-field inflationary system with arbitrary constant quadratic mixing and unequal sound speeds $c_\varphi$ and $c_\sigma$ for the curvature and isocurvature perturbations, respectively. A Laplace representation reduces the coupled dynamics to a second-order Heun equation and gives the exact late-time curvature power spectrum for arbitrary values of the mixing parameter and entropy mass. The sound-speed ratio $r_c=c_\sigma/c_\varphi$ controls the separation between the two sound horizons. At fixed nonzero mixing and $r_c\ll1$, the power exhibits power-law enhancement, growth as $\ln^2(1/r_c)$, or bounded oscillations in $\ln r_c$, depending on the mass and mixing. For $r_c\gg1$, the correction to the unmixed spectrum vanishes as $\ln r_c/r_c^2$ at fixed mass and mixing. Using the same normalized modes, we compute the tree-level bispectrum for three representative cubic interactions, with the mixing parameter treated exactly. For heavy entropy fields, unequal speeds modify the collider amplitude and phase while preserving the strict squeezed frequency set by the final entropy mass. The power-normalized amplitude can vary nonmonotonically with $r_c$, and a sound-speed hierarchy can delay the approach to this asymptotic regime.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13.0 via TypeSafe API, question bank abstract-p1-20260923, 2026-10-02: 1 request, 6,381 input tokens.

What the abstract itself states, as probabilities: claims first 0.03 · states limitation 0.22 · human participants 0.01 · controlled comparison 0.22 · randomized assignment 0.01 · states quantitative result 0.76 · compares with prior results 0.51 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; .

We derive the canonically normalized modes of a two-field inflationary system with arbitrary constant quadratic mixing and unequal sound speeds $c_\varphi$ and $c_\sigma$ for the curvature and isocurvature perturbations, respectively. A Laplace representation reduces the coupled dynamics to a second-order Heun equation and gives the exact late-time curvature power spectrum for arbitrary values of the mixing parameter and entropy mass. The sound-speed ratio $r_c=c_\sigma/c_\varphi$ controls the separation between the two sound horizons. At fixed nonzero mixing and $r_c\ll1$, the power exhibits power-law enhancement, growth as $\ln^2(1/r_c)$, or bounded oscillations in $\ln r_c$, depending on the mass and mixing. For $r_c\gg1$, the correction to the unmixed spectrum vanishes as $\ln r_c/r_c^2$ at fixed mass and mixing. Using the same normalized modes, we compute the tree-level bispectrum for three representative cubic interactions, with the mixing parameter treated exactly. For heavy entropy fields, unequal speeds modify the collider amplitude and phase while preserving the strict squeezed frequency set by the final entropy mass. The power-normalized amplitude can vary nonmonotonically with $r_c$, and a sound-speed hierarchy can delay the approach to this asymptotic regime.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 We derive the canonically normalized modes of a two-field inflationary system with arbitrary constant quadratic mixing and unequal sound speeds $c_\varphi$ and $c_\sigma$ for the curvature and isocurvature perturbations, respectively. contribution 0.93 0.96 0.07 asserted 1.00
s1 A Laplace representation reduces the coupled dynamics to a second-order Heun equation and gives the exact late-time curvature power spectrum for arbitrary values of the mixing parameter and entropy mass. method 0.72 0.94 0.22 asserted 1.00
s2 The sound-speed ratio $r_c=c_\sigma/c_\varphi$ controls the separation between the two sound horizons. result 0.64 0.59 0.08 asserted 1.00
s3 At fixed nonzero mixing and $r_c\ll1$, the power exhibits power-law enhancement, growth as $\ln^2(1/r_c)$, or bounded oscillations in $\ln r_c$, depending on the mass and mixing. result 1.00 0.91 0.21 asserted 0.99
s4 For $r_c\gg1$, the correction to the unmixed spectrum vanishes as $\ln r_c/r_c^2$ at fixed mass and mixing. result 1.00 0.91 0.08 asserted 1.00
s5 Using the same normalized modes, we compute the tree-level bispectrum for three representative cubic interactions, with the mixing parameter treated exactly. contribution 0.82 0.95 0.07 asserted 1.00
s6 For heavy entropy fields, unequal speeds modify the collider amplitude and phase while preserving the strict squeezed frequency set by the final entropy mass. result 0.99 0.89 0.19 asserted 1.00
s7 The power-normalized amplitude can vary nonmonotonically with $r_c$, and a sound-speed hierarchy can delay the approach to this asymptotic regime. result 0.86 0.90 0.12 asserted 0.81

What kind of paper the classifier reports this is: theoretical result 1.00 · evidence from proof or derivation 1.00.

What was examined, and what was not

sources read
  • arXiv listing abstract v1
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. Amplitude and phase of the squeezed signal from the σ 3 interaction as a function of rc = cσ/cφ, at fixed γ = 1 and ν = i, for three quadratic mixing strengths. The left panel shows |bγ,+|/2, while the coefficient of q 3/2 multiplying the cosine in Eq. (1.13) is 2|bγ,+|. The right panel shows arg bγ,+, continued smoothly from its principal value at rc = 1. Each marker is an independent numerical evaluation, and the lines guide the eye; the vertical dotted line marks equal sound speeds. All curves have the same final entropy mass µ/H = √ 13/2 and logarithmic frequency ρ = 1, while…
Figure 2. Equal-speed spectral contours in the complex u-plane. The cuts extend vertically upwards from u = 1 and u = −1, and the Hankel contours Γ+ and Γ− turn counterclockwise around the corresponding branch points. The kernel is analytically continued along each contour, whose two infinite ends generally lie on different sheets. For real τ < 0, the exponential exp(−ikcτu) decays along both upward tails. The endpoints supply the early-time factors exp(−ikcτ ) and exp(ikcτ ), selecting the positive- and negative-frequency sectors, respectively. 2.4 Choosing the contour integrals With the…
Figure 3. The connected positive-frequency contour Γ+ for the two possible orderings of u = 1 and u = rc. All four branch cuts extend upward. The blue contour descends to the left of the two positive endpoints, passes below them, and returns to positive imaginary infinity on their right, with counterclockwise orientation. Both infinite tails extend into the region where exp(−ikcφτu) is damped for real τ < 0. The negative-endpoint cuts, shown in gray, remain outside the contour. As the positive endpoints merge at rc = 1, this path becomes a contour around the single positive cut of the…
Figure 4. Equivalent contour representations for 0 < rc < 1. Left: the connected positive-frequency contour Γ+ passes below both positive endpoints. Right: its integral is decomposed into contributions from the counterclockwise contours Γrc and Γ1, surrounding the individual positive cuts. The kernel on each contour is continued on the sheet inherited from Γ+. All branch cuts extend upward, and the infinite contour tails lie in the region where exp(−ikcφτu) is damped for real τ < 0. The negative￾endpoint cuts, shown in gray, remain outside the integration paths. For rc > 1, the positive…
Figure 5. Contour deformation for the frozen curvature amplitude, shown for 0 < rc < 1. The limit τ → 0 − is taken first, leaving the integrand Yb(u)/(1 − u 2 ). (a) The connected contour Γ+ passes below both positive endpoints. (b) Its left leg is moved onto the imaginary axis and extended to −iR, while the horizontal connector ends at R − iR and the right leg rises along Re u = R, with R > max(1, rc). The contour remains open throughout the deformation and crosses no branch cut. (c) As R → ∞, the horizontal connector and right leg give vanishing contributions for Re ν < 3/2, leaving the…
Figure 6. Exact curvature power spectrum (solid curves) and independent solutions of the coupled 1 time-domain equations (open circles), for two light entropy fields and one heavy example with ν = i. At rc = 1, the spectral reference is the equal-speed gamma-function expression. All three examples have µ 2 > 0 and ν 2 + λ 2 > 0. The dotted lines mark equal sound speeds and the unmixed curvature amplitude. 6.1 Small sound-speed ratio Consider first rc ≪ 1. For the first Heun function in Eq. (5.15), the singularity parameter 1 − r −2 c tends to −∞. The ratio Q1/(1 − r −2 c ) approaches 1/16 −…
Figure 7. Evolution of the curvature and entropy mode sums for λ = 1 and ν = 0.8, shown by the solid and dashed curves, respectively, in units of N 2 k . Writing x = −kcφτ , the vertical lines mark the bare sound-horizon scales x = 1 and x = 1/rc. Time increases to the right. For rc = 0.1, the entropy fluctuation crosses its sound horizon first and transfers an enhanced amplitude to curvature. For rc = 10, the curvature fluctuation crosses first and the correction to its final power is small. more general multifield systems [75, 76]. When the entropy mode can be eliminated, this behavior can…

Figures reproduced from arXiv:2610.02138 by the authors.

References

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.