A bstract We compute correlation functions of the primordial density perturbations when they couple to a gapless, strongly coupled sector of spectator fields — “unparticles” — during inflation. We first derive a four-point function of conformally coupled scalars for all kinematic configurations in de Sitter, which exchanges an unparticle at tree level, by performing direct integration using the Mellin-Barnes method. To obtain inflationary bispectra and trispectra, we apply weight-shifting operators to the conformally coupled scalar correlator. We show that the correlators solve differential equations determined by the additional symmetries enjoyed by the unparticle propagator. Based on these differential equations, we are able to discuss the spinning-unparticle exchanges, focusing on two possible cases where the currents or the stress tensor of unparticles are coupled to inflatons, with the help of spin-raising operators. Finally, we study the phenomenology of the resulting shape functions. Depending on the value of the unparticle scaling dimension, we classify three characteristic shapes for the inflationary bispectra, including near-equilateral, near-orthogonal, and a novel shape which appears when the scaling dimensions are close to half-integers. More generally, we find that the leading order squeezed limits are insufficient to conclusively determine the detection of a light particle or unparticle. Only the full shapes of bispectra and trispectra can break this degeneracy.
Pimentel et al. (2026) studied this question.