Showing 1–50 of 161 results for author: Kumar, S

First Search for Kaluza-Klein Gravitons and Radion Using Planck Data

Authors: Alexander P. Cassem, Soubhik Kumar

Abstract: Heavy moduli and Kaluza-Klein (KK) gravitons from extra dimensions may evade terrestrial probes but can be produced during inflation, generating primordial non-Gaussianity (NG) through unavoidable couplings to density perturbations. In a warped five-dimensional (5D) model, we compute the full radion- and KK-graviton-mediated bispectra and perform the first search for these signals using Planck 201… ▽ More Heavy moduli and Kaluza-Klein (KK) gravitons from extra dimensions may evade terrestrial probes but can be produced during inflation, generating primordial non-Gaussianity (NG) through unavoidable couplings to density perturbations. In a warped five-dimensional (5D) model, we compute the full radion- and KK-graviton-mediated bispectra and perform the first search for these signals using Planck 2018 temperature and polarization data. We find no significant evidence for NG, with the maximum significance being $1.8σ$ for $m_{\rm KK}\approx 1.6H$. We also identify a 5D setup which naturally generates NG with $f_{\rm NL}\sim 1-50$, within the reach of future surveys. △ Less

Submitted 2 July, 2026; originally announced July 2026.

Comments: 20 pages, 12 figures

Scalars at the Cosmological Collider: Full Shapes of Tree Diagrams and Bispectrum Searches using Planck Data

Authors: Soubhik Kumar, Qianshu Lu, Zhong-Zhi Xianyu, Yisong Zhang

Abstract: The Cosmological Collider (CC) provides a unique opportunity to probe the particle spectrum and fundamental interactions at extremely high energies. Massive particles, via their decay into inflaton quanta, can induce a non-analytic, oscillatory, primordial non-Gaussianity (NG), including the bispectrum. At tree level, three classes of such processes contribute to the bispectrum: 'single exchange',… ▽ More The Cosmological Collider (CC) provides a unique opportunity to probe the particle spectrum and fundamental interactions at extremely high energies. Massive particles, via their decay into inflaton quanta, can induce a non-analytic, oscillatory, primordial non-Gaussianity (NG), including the bispectrum. At tree level, three classes of such processes contribute to the bispectrum: 'single exchange', 'double exchange', and 'triple exchange', depending on the number of massive particle propagators. We provide a unified evaluation of all three diagrams and derive the explicit shape functions for the bispectrum, valid across the entire kinematic space. We perform a search for these three processes with the Planck data, finding no evidence for NG. We also consider simple extensions of the minimal scenario that can counter the exponential suppression of the non-analytic signature, and produce on-shell particles with masses $M\gg H$, the Hubble scale during inflation. In particular, we focus on the 'scalar chemical potential' mechanism and extend our previous search to a wider range of chemical potential ($ω$) and $M$, finding global 1.5$σ$ evidence for non-zero NG for the parameter space $ω- M \simeq 3H$. △ Less

Submitted 8 April, 2026; originally announced April 2026.

Comments: 22 pages, 21 figures. Companion paper to the letter arXiv:2603.15728

Cosmological Collider Searches beyond the Hubble Scale with Planck Data

Authors: Soubhik Kumar, Qianshu Lu, Zhong-Zhi Xianyu, Yisong Zhang

Abstract: Searches for primordial non-Gaussianity (NG) has the potential to not only reveal the physics of cosmic inflation, but also the structure of fundamental interactions at the highest energies. The cosmological collider (CC) physics program exemplifies this possibility and demonstrates how searches for oscillatory NG can lead to mass-spin spectroscopy of extremely heavy states. Adopting an effective… ▽ More Searches for primordial non-Gaussianity (NG) has the potential to not only reveal the physics of cosmic inflation, but also the structure of fundamental interactions at the highest energies. The cosmological collider (CC) physics program exemplifies this possibility and demonstrates how searches for oscillatory NG can lead to mass-spin spectroscopy of extremely heavy states. Adopting an effective field theory approach, we find the class of Feynman diagrams that can give the largest NG mediated by a heavy scalar particle with mass $M\sim H$, the inflationary Hubble scale. We compute the full shape of the NG and perform the first search for this shape using Planck data, finding no evidence for NG. This search loses its sensitivity as $M\gg H$ since quantum vacuum fluctuations cannot efficiently produce such heavier particles. We then focus on a mechanism where a chemical potential excites on-shell scalar particles with mass $M\gg H$. Computing the full shapes, we perform the first CC search for particles parametrically heavier than $H$ using Planck data. For a range of chemical potential $ω$ and $M$ satisfying $ω-M \simeq 3H$, we find a global $1.7σ$ evidence for non-zero NG, after taking into account the look-elsewhere effect. △ Less

Submitted 16 March, 2026; originally announced March 2026.

Comments: 8 pages, 3 figures

Cosmological and Astrophysical Constraints on Late First-Order Phase Transitions

Authors: Kylar Greene, Daven Wei Ren Ho, Soubhik Kumar, Yuhsin Tsai

Abstract: First-order cosmological phase transitions (PT) can take place in a dark sector at relatively late times between the big-bang nucleosynthesis and recombination epochs. Because bubble nucleation is stochastic, the PT completes at different times in different regions of the Universe. This fluctuation sources a curvature perturbation whose (dimensionless) power spectrum ${\cal P}_ζ(k)$ features a uni… ▽ More First-order cosmological phase transitions (PT) can take place in a dark sector at relatively late times between the big-bang nucleosynthesis and recombination epochs. Because bubble nucleation is stochastic, the PT completes at different times in different regions of the Universe. This fluctuation sources a curvature perturbation whose (dimensionless) power spectrum ${\cal P}_ζ(k)$ features a universal infrared tail, independent of the microscopic details of the PT. Even in the absence of any non-gravitational interaction between the dark sector and the Standard Model, these additional curvature perturbations at small scales impact a variety of observables. We derive new constraints on dark sector phase transitions using {\it Planck}, baryon acoustic oscillation (BAO), Lyman-$α$ observations, spectral distortion limits from FIRAS, constraints on early reionization, and the existence of ultra-faint dwarf galaxies. △ Less

Submitted 27 February, 2026; originally announced March 2026.

Comments: 32 pages, 6 figures, 2 tables

Freezing-in the Axiverse

Authors: Christopher Dessert, Soubhik Kumar, Joshua T. Ruderman

Abstract: The presence of multiple light axions in the infrared is a generic feature of many ultraviolet (UV) scenarios. In many cases the number of axions ${\cal N}$ is ${\cal O}(10-100)$ or more. Even in the scenario where these axions interact very weakly with the Standard Model (SM), the presence of ${\cal N}$ light axions poses a challenge to the stringent constraint on the number of relativistic degre… ▽ More The presence of multiple light axions in the infrared is a generic feature of many ultraviolet (UV) scenarios. In many cases the number of axions ${\cal N}$ is ${\cal O}(10-100)$ or more. Even in the scenario where these axions interact very weakly with the Standard Model (SM), the presence of ${\cal N}$ light axions poses a challenge to the stringent constraint on the number of relativistic degrees of freedom $N_{\rm eff}$. In order to remain agnostic about the UV, we adopt an effective field theory (EFT) approach, and parametrize the interactions of ${\cal N}$ axions with the SM to quantify the contribution to $N_{\rm eff}$. We consider operators up to dimension six, uncovering one previously-unconsidered charge radius operator, and pay particular attention to the flavor structure of the axion-SM fermion couplings and consider EFTs based on anarchy, textures, and minimal flavor violation. For various choices of such EFTs, we identify the discovery space for current and future cosmic microwave background surveys, including the Simons Observatory and CMB-HD. We show this discovery space depends sensitively on the flavor structure and exhibits a rich interplay with terrestrial and astrophysical probes. △ Less

Submitted 11 February, 2026; v1 submitted 12 November, 2025; originally announced November 2025.

Comments: 38 + 11 pages, 13 figures, 1 table

Warped Dimensions at the Cosmological Collider

Authors: Soubhik Kumar, Michael Nee

Abstract: Extra dimensions are present in many beyond the Standard Model scenarios, most notably in string theory. However, direct signatures of extra dimensions are difficult to observe in many cases. This is the situation, for example, if the energy scales associated with extra dimensions are close to the string or Grand Unification scale. The energetic early universe provides an exciting opportunity to o… ▽ More Extra dimensions are present in many beyond the Standard Model scenarios, most notably in string theory. However, direct signatures of extra dimensions are difficult to observe in many cases. This is the situation, for example, if the energy scales associated with extra dimensions are close to the string or Grand Unification scale. The energetic early universe provides an exciting opportunity to overcome this challenge, since the heavy states associated with high-scale extra dimensions, such as scalar moduli and Kaluza-Klein (KK) gravitons, could have been produced on-shell at early epochs. In this work, we illustrate this by focusing on how such states can be produced during inflation and leave signatures in primordial non-Gaussianity (NG). Specifically, we consider a 5D spacetime with a warped extra dimension that remains stabilized as inflation proceeds in the four non-compact dimensions. By discussing an explicit stabilization mechanism, we compute the masses and couplings of the radion modulus and the KK graviton modes. Being gravitational degrees of freedom, these unavoidably couple to the field(s) generating curvature perturbation, and can lead to observable NG with a distinctive oscillatory shape and characteristic angular dependence. We give example benchmarks which can already be probed by the Planck data and identify targets for the future. Our study shows that cosmological surveys have the potential to observe on-shell imprints of extra dimensions in the coming years. △ Less

Submitted 11 November, 2025; v1 submitted 22 October, 2025; originally announced October 2025.

Comments: 23+10 pages, 4 figures, 1 table, Journal Version

Constraints on Axion-Like Particles from VERITAS Observations of a Flaring Radio Galaxy in the Perseus Cluster

Authors: C. B. Adams, A. Archer, P. Bangale, J. T. Bartkoske, W. Benbow, Y. Chen, J. L. Christiansen, A. J. Chromey, A. Duerr, M. Errando, M. Escobar Godoy, J. Escudero Pedrosa, S. Feldman, Q. Feng, S. Filbert, L. Fortson, A. Furniss, W. Hanlon, O. Hervet, C. E. Hinrichs, J. Holder, Z. Hughes, T. B. Humensky, M. Iskakova, W. Jin , et al. (40 additional authors not shown)

Abstract: Background: Axion-like particles (ALPs) are hypothetical particles that emerge in numerous theoretical extensions to the Standard Model. Their coupling to electromagnetic field implies that ALPs would mix with photons in the presence of external magnetic fields. As ALP phenomenology is governed by the mass and strength of its coupling, there is a subset of this parameter space in which this mixing… ▽ More Background: Axion-like particles (ALPs) are hypothetical particles that emerge in numerous theoretical extensions to the Standard Model. Their coupling to electromagnetic field implies that ALPs would mix with photons in the presence of external magnetic fields. As ALP phenomenology is governed by the mass and strength of its coupling, there is a subset of this parameter space in which this mixing would be expected to leave an imprint on the spectra of TeV gamma-ray sources. Data: In 2017, the VERITAS gamma-ray observatory recorded the second day of a dramatic flare of the radio galaxy NGC 1275, embedded at the center of the Perseus galaxy cluster. This serendipitous locale provides a spatially-extended magnetic field of strength O(10$μ$G) through which escaping photons traverse, making it an excellent target to study ALPs. Methods: We analyze the VERITAS data of NGC 1275's 2017 flare with the gammapy analysis package. Extensive fitting and modeling are performed to ultimately conduct a likelihood analysis used to search for any evidence of a preference for ALPs and to explore the confidence with which constraints can be set. We adopt the CLs method for this study for its conservative approach to setting limits in regimes where the search has limited sensitivity. Results: No evidence for the existence of ALPs is found, and no combination of mass and coupling strength can be excluded at or above 95% confidence level. We provide a map showing the strength of our exclusions in the mass and coupling parameter space. The strongest exclusions are found in the mass range $2 \times 10^{-7}$eV $\lesssim m_a \lesssim 4 \times 10^{-7}$eV and at the coupling strength of $g_{aγ} \gtrsim 3 \times 10^{-11}$ GeV$^{-1}$ up to 80% confidence level, which are consistent with previous studies. Conclusions: We find the CLs method to be a trustworthy approach, and advocate for its... △ Less

Submitted 21 October, 2025; originally announced October 2025.

Comments: 18 pages, 6 figures, accepted for publication in the Physical Review D (PRD)

Cosmological Constraints on Secluded Dark Radiation

Authors: Jae Hyeok Chang, Peizhi Du, Subhajit Ghosh, Soubhik Kumar

Abstract: Dark radiation (DR) is ubiquitous in physics beyond the Standard Model (SM), and its interactions with the SM and dark matter (DM) lead to a variety of interesting effects on cosmological observables. However, even in scenarios where DR is 'secluded', i.e., only gravitationally interacting with SM and DM, it can leave discernible signatures. We present a comprehensive study of four different types… ▽ More Dark radiation (DR) is ubiquitous in physics beyond the Standard Model (SM), and its interactions with the SM and dark matter (DM) lead to a variety of interesting effects on cosmological observables. However, even in scenarios where DR is 'secluded', i.e., only gravitationally interacting with SM and DM, it can leave discernible signatures. We present a comprehensive study of four different types of DR: free-streaming, self-interacting (coupled), decoupling, and recoupling DR, and vary initial conditions to include both adiabatic and isocurvature perturbations. In addition to these properties, we also vary neutrino energy density, DR energy density, and the SM neutrino masses to perform a general analysis and study degeneracies among neutrino and DR properties. We derive constraints using the cosmic microwave background, large-scale structure, and supernova datasets. We find no significant preference for physics beyond the $Λ$CDM model, but data exhibit interesting interplays between different physical quantities. When the neutrino energy density is allowed to vary, we find that the cosmological dataset prefers massless free-streaming DR over massive neutrinos, leading to a significant relaxation of the neutrino mass bound. Although we do not find any evidence of DR isocurvature, the data show support for a strong blue tilt of the isocurvature power spectrum. Our analysis also highlights the degeneracy of various DR parameters with the Hubble constant $H_0$ resulting in a mild relaxation of the $H_0$ tension. △ Less

Submitted 1 October, 2025; originally announced October 2025.

Comments: 40 pages, 13 figures

Report number: FERMILAB-PUB-25-0710-T-V, UT-WI-31-3025

Heavy QCD Axions at High-Energy Muon Colliders

Authors: Ravneet Bedi, Tony Gherghetta, Soubhik Kumar, Peiran Li, Zhen Liu

Abstract: We study the physics potential of heavy QCD axions at high-energy muon colliders. Unlike typical axion-like particles, heavy QCD axions solve the strong CP problem with phenomenology driven by the anomalous gluon ($aG\widetilde G$) couplings. Several ultraviolet scenarios are presented in which QCD axions with TeV-scale masses and decay constants arise consistently with a solution to both the stro… ▽ More We study the physics potential of heavy QCD axions at high-energy muon colliders. Unlike typical axion-like particles, heavy QCD axions solve the strong CP problem with phenomenology driven by the anomalous gluon ($aG\widetilde G$) couplings. Several ultraviolet scenarios are presented in which QCD axions with TeV-scale masses and decay constants arise consistently with a solution to both the strong CP problem and the axion quality problem. We perform a detailed collider analysis for both a 3 and 10~TeV muon collider, focusing on hadronic axion decays that gives rise to a dijet-resonance signature. Our projections for the axion discovery reach in the multi-TeV mass range demonstrate that a muon collider can significantly extend sensitivity to heavy QCD axions compared to existing experiments. △ Less

Submitted 12 September, 2025; originally announced September 2025.

Comments: 30 pages, 11 figures

Love beyond Einstein: Metric reconstruction and Love number in quadratic gravity using WEFT

Authors: Arpan Bhattacharyya, Saptaswa Ghosh, Naman Kumar, Shailesh Kumar, Sounak Pal

Abstract: We study tidal Love numbers of static black holes in four-dimensional quadratic theory of gravity, extending the result of GR. We use worldline effective field theory (WEFT) methods to compute metric perturbations from one-point functions, treating the higher-derivative terms perturbatively. We show that insertions of scalar fields on the worldline induce non-zero tidal tails, and the correspondin… ▽ More We study tidal Love numbers of static black holes in four-dimensional quadratic theory of gravity, extending the result of GR. We use worldline effective field theory (WEFT) methods to compute metric perturbations from one-point functions, treating the higher-derivative terms perturbatively. We show that insertions of scalar fields on the worldline induce non-zero tidal tails, and the corresponding Love number displays no RG running. The same conclusion holds for the insertions of tensor fields. Furthermore, for scalar dipole perturbations, we derive a Yukawa-deformed Frobenius solution and match the asymptotic behavior to fix the UV charge, finding agreement with EFT predictions of Wilson coefficients. Our work demonstrates that quadratic higher-curvature corrections induce non-zero but scale-independent tidal responses, offering a robust EFT framework to test deviations from GR in gravitational wave observations. △ Less

Submitted 9 October, 2025; v1 submitted 4 August, 2025; originally announced August 2025.

Comments: 59 pages, 7 figures, version to appear in JHEP

Journal ref: JHEP 11 (2025) 155

Not-quite-primordial black holes

Authors: Wenzer Qin, Soubhik Kumar, Priyamvada Natarajan, Neal Weiner

Abstract: We propose a new mechanism for the formation of seeds of supermassive black holes at early cosmic epochs. Enhanced density fluctuations with amplitudes that are not large enough to form primordial black holes post-inflation can still lead to collapsed dark matter halos at very early times. For halos forming prior to $1+z \approx 200$, the Cosmic Microwave Background (CMB) is energetic enough to su… ▽ More We propose a new mechanism for the formation of seeds of supermassive black holes at early cosmic epochs. Enhanced density fluctuations with amplitudes that are not large enough to form primordial black holes post-inflation can still lead to collapsed dark matter halos at very early times. For halos forming prior to $1+z \approx 200$, the Cosmic Microwave Background (CMB) is energetic enough to suppress the formation of molecular hydrogen, hence preventing cooling and fragmentation, as a consequence of which baryons falling into the potential well of the halo may undergo "direct collapse" into a black hole. We show using a few illustrative models how this mechanism may account for the abundance of high-redshift black holes inferred from observations by the James Webb Space Telescope while remaining consistent with current limits from CMB spectral distortions. Limits on the primordial power spectrum are also derived by requiring that the universe not reionize too early. △ Less

Submitted 16 June, 2025; originally announced June 2025.

Comments: 9 pages, 5 figures

The CosmoVerse White Paper: Addressing observational tensions in cosmology with systematics and fundamental physics

Authors: Eleonora Di Valentino, Jackson Levi Said, Adam Riess, Agnieszka Pollo, Vivian Poulin, Adrià Gómez-Valent, Amanda Weltman, Antonella Palmese, Caroline D. Huang, Carsten van de Bruck, Chandra Shekhar Saraf, Cheng-Yu Kuo, Cora Uhlemann, Daniela Grandón, Dante Paz, Dominique Eckert, Elsa M. Teixeira, Emmanuel N. Saridakis, Eoin Ó Colgáin, Florian Beutler, Florian Niedermann, Francesco Bajardi, Gabriela Barenboim, Giulia Gubitosi, Ilaria Musella , et al. (516 additional authors not shown)

Abstract: The standard model of cosmology has provided a good phenomenological description of a wide range of observations both at astrophysical and cosmological scales for several decades. This concordance model is constructed by a universal cosmological constant and supported by a matter sector described by the standard model of particle physics and a cold dark matter contribution, as well as very early-t… ▽ More The standard model of cosmology has provided a good phenomenological description of a wide range of observations both at astrophysical and cosmological scales for several decades. This concordance model is constructed by a universal cosmological constant and supported by a matter sector described by the standard model of particle physics and a cold dark matter contribution, as well as very early-time inflationary physics, and underpinned by gravitation through general relativity. There have always been open questions about the soundness of the foundations of the standard model. However, recent years have shown that there may also be questions from the observational sector with the emergence of differences between certain cosmological probes. In this White Paper, we identify the key objectives that need to be addressed over the coming decade together with the core science projects that aim to meet these challenges. These discordances primarily rest on the divergence in the measurement of core cosmological parameters with varying levels of statistical confidence. These possible statistical tensions may be partially accounted for by systematics in various measurements or cosmological probes but there is also a growing indication of potential new physics beyond the standard model. After reviewing the principal probes used in the measurement of cosmological parameters, as well as potential systematics, we discuss the most promising array of potential new physics that may be observable in upcoming surveys. We also discuss the growing set of novel data analysis approaches that go beyond traditional methods to test physical models. [Abridged] △ Less

Submitted 4 August, 2025; v1 submitted 2 April, 2025; originally announced April 2025.

Comments: 416 pages, 81 figures

Journal ref: Phys. Dark Univ. 49 (2025) 101965

New constraints on interacting dark energy from DESI DR2 BAO observations

Authors: Emanuelly Silva, Miguel A. Sabogal, Mateus Scherer, Rafael C. Nunes, Eleonora Di Valentino, Suresh Kumar

Abstract: In its second data release (DR2), the Dark Energy Spectroscopic Instrument (DESI) publicly released measurements of Baryon Acoustic Oscillations (BAO) from over 13.1 million galaxies and 1.6 million quasars, covering the redshift range $0.295 \leq z \leq 2.330$. In this work, we investigate the impact of this new dataset on dark sector interaction models, which are motivated by non-gravitational i… ▽ More In its second data release (DR2), the Dark Energy Spectroscopic Instrument (DESI) publicly released measurements of Baryon Acoustic Oscillations (BAO) from over 13.1 million galaxies and 1.6 million quasars, covering the redshift range $0.295 \leq z \leq 2.330$. In this work, we investigate the impact of this new dataset on dark sector interaction models, which are motivated by non-gravitational interactions between dark energy (DE) and dark matter (DM), commonly referred to as interacting dark energy models (IDE). We focus on two frameworks: the traditional IDE model and the recently proposed sign-switching Interacting model (S-IDE), aiming to derive new and robust constraints on both scenarios. After carefully selecting the sample for the joint analysis, ensuring compatibility among the data without significant tension, our main results indicate that both models can alleviate the $H_0$ tension, reducing it to moderate tension approximately $2.7σ$. The IDE model shows compatibility with the latest $S_8$ constraints from cosmic shear surveys, while the S-IDE model predicts lower values of $S_8$, which align with alternative perspectives on the $S_8$ tension. For the traditional IDE model, we derive new bounds for the coupling parameter, marking the strongest constraints to date through geometric measurements. This highlights the crucial role that supernova samples can play in refining these constraints. For the S-IDE model, we find mild evidence (over $2σ$) for a non-zero coupling, once the PantheonPlus dataset calibrated with Cepheid-based magnitude measurements is included in the analysis. △ Less

Submitted 11 June, 2025; v1 submitted 29 March, 2025; originally announced March 2025.

Comments: 13 pages, 5 figures, 4 tables. Minor changes and updated references

Journal ref: Phys. Rev. D 111, 123511 (2025)

Early Galaxies from Rare Inflationary Processes and JWST Observations

Authors: Soubhik Kumar, Neal Weiner

Abstract: Rare Poisson processes (PP) during cosmic inflation can lead to signatures that are localized in position space and are not well captured by the standard two- or higher-point correlation functions of primordial density perturbations. As an example, PP can lead to localized overdense regions that are far denser than the ones produced through standard inflationary fluctuations. As a result, such ove… ▽ More Rare Poisson processes (PP) during cosmic inflation can lead to signatures that are localized in position space and are not well captured by the standard two- or higher-point correlation functions of primordial density perturbations. As an example, PP can lead to localized overdense regions that are far denser than the ones produced through standard inflationary fluctuations. As a result, such overdense regions collapse earlier than expected based on the standard $Λ$CDM model and would host anomalously high-redshift galaxies. We describe some general aspects of such PP and consider a particular realization in the context of inflationary particle production. We then show that the masses and redshifts of the resulting galaxies can lie in a range discoverable by the James Webb Space Telescope (JWST) and future surveys, while being consistent with existing constraints on the matter power spectrum and UV luminosity functions at lower redshifts. △ Less

Submitted 12 February, 2025; originally announced February 2025.

Comments: 14 pages, 9 figures

Visions in Quantum Gravity

Authors: Luca Buoninfante, Benjamin Knorr, K. Sravan Kumar, Alessia Platania, Damiano Anselmi, Ivano Basile, N. Emil J. Bjerrum-Bohr, Robert Brandenberger, Mariana Carrillo González, Anne-Christine Davis, Bianca Dittrich, Paolo Di Vecchia, John F. Donoghue, Fay Dowker, Gia Dvali, Astrid Eichhorn, Steven B. Giddings, Alessandra Gnecchi, Giulia Gubitosi, Lavinia Heisenberg, Renata Kallosh, Alexey S. Koshelev, Stefano Liberati, Renate Loll, Leonardo Modesto , et al. (14 additional authors not shown)

Abstract: To deepen our understanding of Quantum Gravity and its connections with black holes and cosmology, building a common language and exchanging ideas across different approaches is crucial. The Nordita Program "Quantum Gravity: from gravitational effective field theories to ultraviolet complete approaches" created a platform for extensive discussions, aimed at pinpointing both common grounds and sour… ▽ More To deepen our understanding of Quantum Gravity and its connections with black holes and cosmology, building a common language and exchanging ideas across different approaches is crucial. The Nordita Program "Quantum Gravity: from gravitational effective field theories to ultraviolet complete approaches" created a platform for extensive discussions, aimed at pinpointing both common grounds and sources of disagreements, with the hope of generating ideas and driving progress in the field. This contribution summarizes the twelve topical discussions held during the program and collects individual thoughts of speakers and panelists on the future of the field in light of these discussions. △ Less

Submitted 8 October, 2025; v1 submitted 11 December, 2024; originally announced December 2024.

Comments: Collection of summaries of twelve topical panel discussions and individual thoughts of speakers and panelists, Nordita Scientific Program "Quantum Gravity: from gravitational EFTs to UV complete approaches". 63 pages + references, no figures - V2: typos fixed, references and one contribution added - V3: matches published version

Journal ref: SciPost Phys. Comm. Rep. 11 (2025)

Towards a Complete Treatment of Scalar-induced Gravitational Waves with Early Matter Domination

Authors: Soubhik Kumar, Hanwen Tai, Lian-Tao Wang

Abstract: Large curvature perturbations can source an observable amount of stochastic gravitational wave background (SGWB). We consider several scenarios where small-scale curvature perturbations are naturally enhanced due to the presence of additional spectator fields during inflation. The same spectator fields can lead to a period of early matter domination (EMD) after inflation. We compute the inflationa… ▽ More Large curvature perturbations can source an observable amount of stochastic gravitational wave background (SGWB). We consider several scenarios where small-scale curvature perturbations are naturally enhanced due to the presence of additional spectator fields during inflation. The same spectator fields can lead to a period of early matter domination (EMD) after inflation. We compute the inflationary spectrum of curvature perturbation and determine its evolution at later times, taking into account both the onset and the end of the EMD epoch, and also the impact of relative velocity perturbation between matter and radiation. The feature that the same field is responsible for both enhanced perturbations and the EMD era, leads to a predictive framework within which the full frequency dependence of SGWB can be computed. The SGWB can be observed in several detectors, including those focused on the nano-Hz regime. Our numerical framework can also be used to study other non-standard cosmological histories. △ Less

Submitted 21 October, 2024; originally announced October 2024.

Comments: 38 pages, 14 figures

An indirect search for dark matter with a combined analysis of dwarf spheroidal galaxies from VERITAS

Authors: A. Acharyya, C. B. Adams, P. Bangale, J. T. Bartkoske, P. Batista, W. Benbow, J. L. Christiansen, A. J. Chromey, A. Duerr, M. Errando, A. Falcone, Q. Feng, G. M. Foote, L. Fortson, A. Furniss, W. Hanlon, D. Hanna, O. Hervet, C. E. Hinrichs, J. Holder, T. B. Humensky, W. Jin, M. N. Johnson, P. Kaaret, M. Kertzman , et al. (37 additional authors not shown)

Abstract: Understanding the nature and identity of dark matter is a key goal in the physics community. In the case that TeV-scale dark matter particles decay or annihilate into standard model particles, very-high-energy (VHE) gamma rays (greater than 100 GeV) will be present in the final state. The Very Energetic Radiation Imaging Telescope Array System (VERITAS) is an imaging atmospheric Cherenkov telescop… ▽ More Understanding the nature and identity of dark matter is a key goal in the physics community. In the case that TeV-scale dark matter particles decay or annihilate into standard model particles, very-high-energy (VHE) gamma rays (greater than 100 GeV) will be present in the final state. The Very Energetic Radiation Imaging Telescope Array System (VERITAS) is an imaging atmospheric Cherenkov telescope array that can indirectly detect VHE gamma rays in an energy range of 100 GeV to > 30 TeV. Dwarf spheroidal galaxies (dSphs) are ideal candidates in the search for dark matter due to their high dark matter content, high mass-to-light ratios, and their low gamma-ray fluxes from astrophysical processes. This study uses a legacy data set of 638 hours collected on 17 dSphs, built over 11 years with an observing strategy optimized according to the dark matter content of the targets. The study addresses a broad dark matter particle mass range, extending from 200 GeV to 30 PeV. In the absence of a detection, we set the upper limits on the dark matter velocity-weighted annihilation cross section. △ Less

Submitted 6 August, 2024; v1 submitted 23 July, 2024; originally announced July 2024.

Comments: 18 pages, 7 figures, 3 tables, accepted in PRD

$Λ_{\rm s}$CDM cosmology: Alleviating major cosmological tensions by predicting standard neutrino properties

Authors: Anita Yadav, Suresh Kumar, Cihad Kibris, Ozgur Akarsu

Abstract: We investigate a two-parameter extension of the $Λ_{\rm s}$CDM model by allowing variations in the effective number of neutrino species $N_{\rm eff}$ and their total mass $\sum m_ν$. Our motivation is twofold: (i) to examine whether $Λ_{\rm s}$CDM retains its success in fitting the data and addressing major cosmological tensions, without suggesting a need for a deviation from the standard model of… ▽ More We investigate a two-parameter extension of the $Λ_{\rm s}$CDM model by allowing variations in the effective number of neutrino species $N_{\rm eff}$ and their total mass $\sum m_ν$. Our motivation is twofold: (i) to examine whether $Λ_{\rm s}$CDM retains its success in fitting the data and addressing major cosmological tensions, without suggesting a need for a deviation from the standard model of particle physics, and (ii) to determine whether the data indicate new physics that could potentially address cosmological tensions, either in the post-recombination universe through the late-time mirror AdS-dS transition, or in the pre-recombination universe through modifications in the standard values of $N_{\rm eff}$ and $\sum m_ν$, or both. Within the extended $Λ_{\rm s}$CDM model, referred to as $Λ_{\rm s}$CDM+$N_{\rm eff}$+$\sum m_{\rm ν}$, we find no significant tension when considering the Planck-alone analysis. We observe that incorporating BAO data limits the further success of the $Λ_{\rm s}$CDM extension. However, the weakly model-dependent BAOtr data, along with Planck and Planck+PP\&SH0ES, favor $H_0\sim 73\,{\rm km\, s^{-1}\, Mpc^{-1}}$. In cases where BAOtr dataset is used, the mirror AdS-dS transition is very effective in providing enhanced $H_0$ values, and thus the model requires no significant deviation from the standard value of $N_{\rm eff} = 3.044$. Both the $H_0$ and $S_8$ tensions are effectively addressed, with some compromise in the case of the Planck+BAO dataset. Finally, the upper bounds obtained on $\sum m_ν\lesssim 0.5$~eV are fully compatible with neutrino oscillation experiments. Our findings provide evidence that late-time physics beyond $Λ$CDM, such as $Λ_{\rm s}$CDM, without altering the standard pre-recombination universe, can suffice to alleviate the major cosmological tensions. △ Less

Submitted 13 January, 2025; v1 submitted 26 June, 2024; originally announced June 2024.

Comments: 24 pages, 6 figures, 3 tables, matches the version published in JCAP

Journal ref: JCAP 01(2025)042

Cosmological constraints on $Λ_{\rm s}$CDM scenario in a type II minimally modified gravity

Authors: Özgür Akarsu, Antonio De Felice, Eleonora Di Valentino, Suresh Kumar, Rafael C. Nunes, Emre Özülker, J. Alberto Vazquez, Anita Yadav

Abstract: The idea of a rapid sign-switching cosmological constant (mirror AdS-dS transition) in the late universe at $z\sim1.7$, known as the $Λ_{\rm s}$CDM model, has significantly improved the fit to observational data and provides a promising scenario for alleviating major cosmological tensions, such as the $H_0$ and $S_8$ tensions. However, in the absence of a fully predictive model, implementing this… ▽ More The idea of a rapid sign-switching cosmological constant (mirror AdS-dS transition) in the late universe at $z\sim1.7$, known as the $Λ_{\rm s}$CDM model, has significantly improved the fit to observational data and provides a promising scenario for alleviating major cosmological tensions, such as the $H_0$ and $S_8$ tensions. However, in the absence of a fully predictive model, implementing this fit required conjecturing that the dynamics of the linear perturbations are governed by general relativity. Recent work embedding the $Λ_{\rm s}$CDM model with the Lagrangian of a type II minimally modified gravity known as VCDM has propelled $Λ_{\rm s}$CDM to a fully predictive model, removing the uncertainty related to the aforementioned assumption; we call this new model $Λ_{\rm s}$VCDM. In this work, we demonstrate that not only does $Λ_{\rm s}$CDM fit the data better than the standard $Λ$CDM model, but the new model, $Λ_{\rm s}$VCDM, performs even better in alleviating cosmological tensions while also providing a better fit to the data, including CMB, BAO, SNe Ia, and cosmic shear measurements. Our findings highlight the $Λ_{\rm s}$CDM framework, particularly the $Λ_{\rm s}$VCDM model, as a compelling alternative to the standard $Λ$CDM model, especially by successfully alleviating the $H_0$ tension. Additionally, these models predict higher values for $σ_8$, indicating enhanced structuring, albeit with lower present-day matter density parameter values and consequently reduced $S_8$ values, alleviating the $S_8$ tension as well. This demonstrates that the data are well fit by a combination of background and linear perturbations, both having dynamics differing from those of $Λ$CDM. This paves the way for further exploration of new ways for embedding the sign-switching cosmological constant into other models. △ Less

Submitted 25 November, 2024; v1 submitted 11 June, 2024; originally announced June 2024.

Comments: 17 pages, 6 figures, 2 tables; matches the version published in Physical Review D

Report number: YITP-24-57

Journal ref: Phys. Rev. D 110, 103527 (2024)

Supernova electron-neutrino interactions with xenon in the nEXO detector

Authors: nEXO Collaboration, S. Hedges, S. Al Kharusi, E. Angelico, J. P. Brodsky, G. Richardson, S. Wilde, A. Amy, A. Anker, I. J. Arnquist, P. Arsenault, A. Atencio, I. Badhrees, J. Bane, V. Belov, E. P. Bernard, T. Bhatta, A. Bolotnikov, J. Breslin, P. A. Breur, E. Brown, T. Brunner, E. Caden, G. F. Cao, L. Q. Cao , et al. (122 additional authors not shown)

Abstract: Electron-neutrino charged-current interactions with xenon nuclei were modeled in the nEXO neutrinoless double-$β$ decay detector (~5 metric ton, 90% ${}^{136}$Xe, 10% ${}^{134}$Xe) to evaluate its sensitivity to supernova neutrinos. Predictions for event rates and detectable signatures were modeled using the Model of Argon Reaction Low Energy Yields (MARLEY) event generator. We find good agreement… ▽ More Electron-neutrino charged-current interactions with xenon nuclei were modeled in the nEXO neutrinoless double-$β$ decay detector (~5 metric ton, 90% ${}^{136}$Xe, 10% ${}^{134}$Xe) to evaluate its sensitivity to supernova neutrinos. Predictions for event rates and detectable signatures were modeled using the Model of Argon Reaction Low Energy Yields (MARLEY) event generator. We find good agreement between MARLEY's predictions and existing theoretical calculations of the inclusive cross sections at supernova neutrino energies. The interactions modeled by MARLEY were simulated within the nEXO simulation framework and were run through an example reconstruction algorithm to determine the detector's efficiency for reconstructing these events. The simulated data, incorporating the detector response, were used to study the ability of nEXO to reconstruct the incident electron-neutrino spectrum and these results were extended to a larger xenon detector of the same isotope enrichment. We estimate that nEXO will be able to observe electron-neutrino interactions with xenon from supernovae as far as 5-8 kpc from Earth, while the ability to reconstruct incident electron-neutrino spectrum parameters from observed interactions in nEXO is limited to closer supernovae. △ Less

Submitted 29 November, 2024; v1 submitted 29 May, 2024; originally announced May 2024.

Comments: 17 pages, 16 figures

Report number: LLNL-JRNL-864783

Journal ref: Phys. Rev. D 110, 093002 (2024)

Probing Dark Matter Isocurvature with Primordial Non-Gaussianity

Authors: Michael Geller, Soubhik Kumar, Lian-Tao Wang

Abstract: Multiple fields can become dynamical during the inflationary epoch. We consider an example where a light field acquires isocurvature fluctuations during inflation and contributes to the dark matter abundance at late times. Interactions between the light field and the adiabatic sector contribute to mixed adiabatic-isocurvature non-Gaussianity (NG). We show the resulting form of NG has a different k… ▽ More Multiple fields can become dynamical during the inflationary epoch. We consider an example where a light field acquires isocurvature fluctuations during inflation and contributes to the dark matter abundance at late times. Interactions between the light field and the adiabatic sector contribute to mixed adiabatic-isocurvature non-Gaussianity (NG). We show the resulting form of NG has a different kinematic dependence than the 'local shape' commonly considered, and highlight the parameter space where a dedicated search is expected to significantly improve the current $\textit{Planck}$ sensitivity. We interpret our results in the context of the QCD axion and illustrate how the proposed NG searches can improve upon the existing searches for isocurvature power spectrum and bispectrum. △ Less

Submitted 15 May, 2024; originally announced May 2024.

Comments: 17 pages, 4 figures

Journal ref: JCAP 12 (2024) 018

Too Hot to Handle: Searching for Inflationary Particle Production in Planck Data

Authors: Oliver H. E. Philcox, Soubhik Kumar, J. Colin Hill

Abstract: Non-adiabatic production of massive particles is a generic feature of many inflationary mechanisms. If sufficiently massive, these particles can leave features in the cosmic microwave background (CMB) that are not well-captured by traditional correlation function analyses. We consider a scenario in which particle production occurs only in a narrow time-interval during inflation, eventually leading… ▽ More Non-adiabatic production of massive particles is a generic feature of many inflationary mechanisms. If sufficiently massive, these particles can leave features in the cosmic microwave background (CMB) that are not well-captured by traditional correlation function analyses. We consider a scenario in which particle production occurs only in a narrow time-interval during inflation, eventually leading to CMB hot- or coldspots with characteristic shapes and sizes. Searching for such features in CMB data is analogous to searching for late-Universe hot- or coldspots, such as those due to the thermal Sunyaev-Zel'dovich (tSZ) effect. Exploiting this data-analysis parallel, we perform a search for particle-production hotspots in the Planck PR4 temperature dataset, which we implement via a matched-filter analysis. Our pipeline is validated on synthetic observations and found to yield unbiased constraints on sufficiently large hotspots across $\approx 60\%$ of the sky. After removing point sources and tSZ clusters, we find no evidence for new physics and place novel bounds on the coupling between the inflaton and massive particles. These bounds are strongest for larger hotspots, produced early in inflation, whilst sensitivity to smaller hotspots is limited by noise and beam effects. Through such methods we can constrain particles with masses $\mathcal{O}(100)$ times larger than the inflationary Hubble scale, which represents possibly the highest energies ever directly probed with observational data. △ Less

Submitted 26 April, 2025; v1 submitted 6 May, 2024; originally announced May 2024.

Comments: 12 pages, 7 figures, accepted by Phys. Rev. D

On Lepton Flavor Violation and Dark Matter in Scotogenic model with Trimaximal Mixing

Authors: Tapender, Surender Verma, Sanjeev Kumar

Abstract: We examine the Scotogenic model employing the TM$_2$ mixing matrix, $U_{\text{TM}_2}$, for neutrinos and parameterize the Yukawa coupling matrix $y$ based on the diagonalization condition for the neutrino mass matrix, $m_ν$. Our investigation centers on analyzing the relic density of cold dark matter ($Ωh^2$) and possible lepton flavor violation (LFV) in the model. In particular, we study coannihi… ▽ More We examine the Scotogenic model employing the TM$_2$ mixing matrix, $U_{\text{TM}_2}$, for neutrinos and parameterize the Yukawa coupling matrix $y$ based on the diagonalization condition for the neutrino mass matrix, $m_ν$. Our investigation centers on analyzing the relic density of cold dark matter ($Ωh^2$) and possible lepton flavor violation (LFV) in the model. In particular, we study coannihilation dynamics and LFV, in the model, considering various coannihilation scenarios including non-zero mass splitting between lightest sterile neutrinos. While analyzing, we have taken into consideration respective experimental constraints on $Ωh^2$ and LFV alongside neutrino oscillation data. Our study reveals that in both normal and inverted hierarchy of neutrino masses, splitting between masses of $N_1$ and $N_2$ can be up to $\approx 15\%$ for the model to be in consonance with the above constraints. In the second part, we have extended the analysis incorporating extended magic symmetry in $m_ν$ enabling us to completely determine Yukawa coupling matrix ($y$). We observe a notable exclusion of the effective Majorana mass $|m_{ee}|$ parameter space by cosmological bound on sum of neutrino masses, particularly in the normal hierarchy while inverted hierarchy scenario is excluded due to constraints coming from extended magic symmetry. These findings shed light on the interplay among the Scotogenic model, TM$_2$ mixing, and extended magic symmetry, offering insights into the permitted parameter space and hierarchy exclusion. △ Less

Submitted 26 June, 2024; v1 submitted 26 February, 2024; originally announced February 2024.

Comments: 40 pages, 11 figures

$Λ_{\rm s}$CDM cosmology from a type-II minimally modified gravity

Authors: Özgür Akarsu, Antonio De Felice, Eleonora Di Valentino, Suresh Kumar, Rafael C. Nunes, Emre Özülker, J. Alberto Vazquez, Anita Yadav

Abstract: We integrate $Λ_{\rm s}$CDM, a promising scenario for alleviating cosmological tensions, into VCDM, a type-II minimally modified gravity. This promotes the scenario to a fully predictive model (dubbed $Λ_{\rm s}$VCDM) that specifies the cosmological evolution self-consistently, including through the late-time AdS-to-dS transition epoch. In this theory, an auxiliary scalar field generates an effect… ▽ More We integrate $Λ_{\rm s}$CDM, a promising scenario for alleviating cosmological tensions, into VCDM, a type-II minimally modified gravity. This promotes the scenario to a fully predictive model (dubbed $Λ_{\rm s}$VCDM) that specifies the cosmological evolution self-consistently, including through the late-time AdS-to-dS transition epoch. In this theory, an auxiliary scalar field generates an effective cosmological constant with either a constant or a linear potential. This allows an abrupt mirror AdS-to-dS transition via a piecewise-linear potential with a sudden slope change. To remove the associated sudden singularity and ensure stable evolution, we smooth the junction using a blended sigmoid interpolant, obtaining rapid but continuous transitions. We identify two qualitatively distinct smooth mirror AdS-to-dS realisations of $Λ_{\rm s}$: (i) an agitated transition, in which the potential interpolates between equal-magnitude AdS and dS plateaus and $Λ_{\rm s}$ develops a central bump; and (ii) a quiescent transition, in which the potential remains continuous but changes slope across the transition layer, so that $Λ_{\rm s}(a)$ can remain monotone, with possible shallow shoulders, and a central bump is not automatic. Depending on type and sharpness, a finite-width transition can induce a transient accelerated-expansion interval ($\ddot a>0$) around $z\sim 1.5-2$, in addition to present-day acceleration, and, if the background enters a region where $V_{,φφ}>2/3$, a nested super-acceleration episode. These distinct transient histories can imprint signatures on background and perturbation evolution. Our construction enables a self-consistent observational assessment of smooth $Λ_{\rm s}$CDM realisations and motivates multi-probe analyses to test transition dynamics and reassess cosmological tensions. △ Less

Submitted 28 April, 2026; v1 submitted 12 February, 2024; originally announced February 2024.

Comments: 28 pages (including appendices), 4 figures, and no tables; matches the version published in MNRAS

Report number: YITP-24-18

Journal ref: Mon.Not.Roy.Astron.Soc. 546 (2026) 1, staf2276

An Effective Cosmological Collider

Authors: Nathaniel Craig, Soubhik Kumar, Amara McCune

Abstract: Effective field theories (EFTs) of heavy particles coupled to the inflaton are rife with operator redundancies, frequently obscured by sensitivity to both boundary terms and field redefinitions. We initiate a systematic study of these redundancies by establishing a minimal operator basis for an archetypal example, the abelian gauge-Higgs-inflaton EFT. Working up to dimension 9, we show that certai… ▽ More Effective field theories (EFTs) of heavy particles coupled to the inflaton are rife with operator redundancies, frequently obscured by sensitivity to both boundary terms and field redefinitions. We initiate a systematic study of these redundancies by establishing a minimal operator basis for an archetypal example, the abelian gauge-Higgs-inflaton EFT. Working up to dimension 9, we show that certain low-dimensional operators are entirely redundant and identify new non-redundant operators with potentially interesting cosmological collider signals. Our methods generalize straightforwardly to other EFTs of heavy particles coupled to the inflaton. △ Less

Submitted 19 January, 2024; originally announced January 2024.

Comments: 31 pages, 4 figures

Journal ref: JHEP 07 (2024) 108

The Cosmological Dynamics of String Theory Axion Strings

Authors: Joshua N. Benabou, Quentin Bonnefoy, Malte Buschmann, Soubhik Kumar, Benjamin R. Safdi

Abstract: The quantum chromodynamics (QCD) axion may solve the strong CP problem and explain the dark matter (DM) abundance of our Universe. The axion was originally proposed to arise as the pseudo-Nambu Goldstone boson of global $\mathrm{U}(1)_{\rm PQ}$ Peccei-Quinn (PQ) symmetry breaking, but axions also arise generically in string theory as zero modes of higher-dimensional gauge fields. In this work we s… ▽ More The quantum chromodynamics (QCD) axion may solve the strong CP problem and explain the dark matter (DM) abundance of our Universe. The axion was originally proposed to arise as the pseudo-Nambu Goldstone boson of global $\mathrm{U}(1)_{\rm PQ}$ Peccei-Quinn (PQ) symmetry breaking, but axions also arise generically in string theory as zero modes of higher-dimensional gauge fields. In this work we show that string theory axions behave fundamentally differently from field theory axions in the early Universe. Field theory axions may form axion strings if the PQ phase transition takes place after inflation. In contrast, we show that string theory axions do not generically form axion strings. In special inflationary paradigms, such as D-brane inflation, string theory axion strings may form; however, their tension is parametrically larger than that of field theory axion strings. We then show that such QCD axion strings overproduce the DM abundance for all allowed QCD axion masses and are thus ruled out, except in scenarios with large warping. A loop-hole to this conclusion arises in the axiverse, where an axion string could be composed of multiple different axion mass eigenstates; a heavier eigenstate could collapse the network earlier, allowing for the QCD axion to produce the correct DM abundance and also generating observable gravitational wave signals. △ Less

Submitted 13 December, 2023; originally announced December 2023.

Comments: 36 pages, 6 figures, video abstract at https://youtu.be/LnN9W9Vkxt4

Neutrino Mass Matrices with Generalized CP Symmetries and Texture Zeros

Authors: Sanjeev Kumar, Radha Raman Gautam

Abstract: We investigate the properties of neutrino mass matrices that incorporate texture zeros and generalized CP symmetries associated with tribimaximal mixing. By combining these approaches, we derive predictive neutrino mass matrices and explore their implications for mass hierarchies, mixing angles, and CP-violating phases. We find that the three angles defining the generalized CP symmetries have narr… ▽ More We investigate the properties of neutrino mass matrices that incorporate texture zeros and generalized CP symmetries associated with tribimaximal mixing. By combining these approaches, we derive predictive neutrino mass matrices and explore their implications for mass hierarchies, mixing angles, and CP-violating phases. We find that the three angles defining the generalized CP symmetries have narrow allowed ranges. We also obtain distinct correlations between the three mixing angles and the CP-violating phases that distinguish the various texture patterns from one another. Moreover, we compute the effective neutrino mass for neutrinoless double beta decay and the sum of neutrino masses. Our results highlight the predictability and testability of neutrino mass matrices with generalized CP symmetry. △ Less

Submitted 26 March, 2024; v1 submitted 12 December, 2023; originally announced December 2023.

Comments: 22 pages, 5 figures, 4 tables

Journal ref: Nucl. Phys. B 1001 (2024) 116520

Finite Bubble Statistics Constrain Late Cosmological Phase Transitions

Authors: Gilly Elor, Ryusuke Jinno, Soubhik Kumar, Robert McGehee, Yuhsin Tsai

Abstract: We consider first order cosmological phase transitions (PT) happening at late times, below Standard Model (SM) temperatures $T_{\rm PT} \lesssim$ GeV. The inherently stochastic nature of bubble nucleation and the finite number of bubbles associated with a late-time PT lead to superhorizon fluctuations in the PT completion time. We compute how such fluctuations eventually source curvature fluctuati… ▽ More We consider first order cosmological phase transitions (PT) happening at late times, below Standard Model (SM) temperatures $T_{\rm PT} \lesssim$ GeV. The inherently stochastic nature of bubble nucleation and the finite number of bubbles associated with a late-time PT lead to superhorizon fluctuations in the PT completion time. We compute how such fluctuations eventually source curvature fluctuations with universal properties, independent of the microphysics of the PT dynamics. Using Cosmic Microwave Background (CMB) and Large Scale Structure (LSS) measurements, we constrain the energy released in a dark-sector PT. For 0.1 eV $\lesssim T_{\rm PT} \lesssim$ keV this constraint is stronger than both the current bound from additional neutrino species $ΔN_{\rm eff}$, and in some cases, even CMB-S4 projections. Future measurements of CMB spectral distortions and pulsar timing arrays will also provide competitive sensitivity for keV $\lesssim T_{\rm PT} \lesssim$ GeV. △ Less

Submitted 27 November, 2023; originally announced November 2023.

Comments: 5+7 pages, 2+5 figures

Report number: FTPI-MINN-23-20, UTWI-39-2023

Neutrino Phenomenology in a Model with Generalized CP symmetry within Type-I seesaw framework

Authors: Tapender, Sanjeev Kumar, Surender Verma

Abstract: We investigate the consequences of generalized CP (GCP) symmetry within the context of the two Higgs doublet model (2HDM), specifically focusing on the lepton sector. Utilizing the Type-I seesaw framework, we study an intriguing connection between the Dirac Yukawa couplings originating from both Higgs fields, leading to a reduction in the number of independent Yukawa couplings and simplifying the… ▽ More We investigate the consequences of generalized CP (GCP) symmetry within the context of the two Higgs doublet model (2HDM), specifically focusing on the lepton sector. Utilizing the Type-I seesaw framework, we study an intriguing connection between the Dirac Yukawa couplings originating from both Higgs fields, leading to a reduction in the number of independent Yukawa couplings and simplifying the scalar and Yukawa sectors when compared to the general 2HDM. The CP3 constraint results in two right-handed neutrinos having equal masses and leads to a diagonal right-handed Majorana neutrino mass matrix. Notably, CP symmetry experiences a soft break due to the phase associated with the vacuum expectation value of the second Higgs doublet. The model aligns well with observed charged lepton masses and neutrino oscillation data, explaining both masses and mixing angles, and yields distinct predictions for normal and inverted neutrino mass hierarchies. It features a novel interplay between atmospheric mixing angle $θ_{23}$ and neutrino mass hierarchy: the angle $θ_{23}$ is below maximal for the normal hierarchy and above maximal for inverted hierarchy. Another interesting feature of the model is inherent CP violation for the inverted hierarchy. △ Less

Submitted 8 September, 2023; originally announced September 2023.

Comments: 18 pages, 6 figures

Signatures of Primordial Energy Injection from Axion Strings

Authors: Joshua N. Benabou, Malte Buschmann, Soubhik Kumar, Yujin Park, Benjamin R. Safdi

Abstract: Axion strings are horizon-size topological defects that may be produced in the early Universe. Ultra-light axion-like particles may form strings that persist to temperatures below that of big bang nucleosynthesis. Such strings have been considered previously as sources of gravitational waves and cosmic microwave background (CMB) polarization rotation. In this work we show, through analytic argumen… ▽ More Axion strings are horizon-size topological defects that may be produced in the early Universe. Ultra-light axion-like particles may form strings that persist to temperatures below that of big bang nucleosynthesis. Such strings have been considered previously as sources of gravitational waves and cosmic microwave background (CMB) polarization rotation. In this work we show, through analytic arguments and dedicated adaptive mesh refinement cosmological simulations, that axion strings deposit a sub-dominant fraction of their energy into high-energy Standard Model (SM) final states, for example, by the direct production of heavy radial modes that subsequently decay to SM particles. This high-energy SM radiation is absorbed by the primordial plasma, leading to novel signatures in precision big bang nucleosynthesis, the CMB power spectrum, and gamma-ray surveys. In particular, we show that CMB power spectrum data constrains axion strings with decay constants $f_a \lesssim 10^{12}$ GeV, up to model dependence on the ultraviolet completion, for axion masses $m_a \lesssim 10^{-29}$ eV; future CMB surveys could find striking evidence of axion strings with lower decay constants. △ Less

Submitted 2 August, 2023; originally announced August 2023.

Comments: 30 pages, 20 figures, Supplementary Animations at https://goo.by/qHk9d, Video Abstract at https://youtu.be/Ae4AR46qKnU

Unraveling Cosmological Anisotropies within Stochastic Gravitational Wave Backgrounds

Authors: Yanou Cui, Soubhik Kumar, Raman Sundrum, Yuhsin Tsai

Abstract: Identifying the anisotropies in a cosmologically sourced stochastic gravitational wave background (SGWB) would be of significance in shedding light on the nature of primordial inhomogeneities. For example, if SGWB carries isocurvature fluctuations, it would provide evidence for a multi-field inflationary origin of these inhomogeneities. However, this is challenging in practice due to finite detect… ▽ More Identifying the anisotropies in a cosmologically sourced stochastic gravitational wave background (SGWB) would be of significance in shedding light on the nature of primordial inhomogeneities. For example, if SGWB carries isocurvature fluctuations, it would provide evidence for a multi-field inflationary origin of these inhomogeneities. However, this is challenging in practice due to finite detector sensitivity and also the presence of the astrophysical foregrounds that can compete with the cosmological signal. In this work, we explore the prospects for measuring cosmological SGWB anisotropies in the presence of an astrophysical counterpart and detector noise. To illustrate the main idea, we perform a Fisher analysis using a well-motivated cosmological SGWB template corresponding to a first order phase transition, and an astrophysical SGWB template corresponding to extra-galactic binary mergers, and compute the uncertainty with which various parameters characterizing the isotropic and anisotropic components can be extracted. We also discuss some subtleties and caveats involving shot noise in the astrophysical foreground. Overall, we show that upcoming experiments, e.g., LISA, Taiji, Einstein Telescope, Cosmic Explorer, and BBO, can all be effective in discovering plausible anisotropic cosmological SGWBs. △ Less

Submitted 19 July, 2023; originally announced July 2023.

Comments: 24 pages, 9 figures

Journal ref: JCAP10(2023)064

Gravitational Waves from Stochastic Scalar Fluctuations

Authors: Reza Ebadi, Soubhik Kumar, Amara McCune, Hanwen Tai, Lian-Tao Wang

Abstract: We present a novel mechanism for gravitational wave generation in the early Universe. Light spectator scalar fields during inflation can acquire a blue-tilted power spectrum due to stochastic effects. We show that this effect can lead to large curvature perturbations at small scales (induced by the spectator field fluctuations) while maintaining the observed, slightly red-tilted curvature perturba… ▽ More We present a novel mechanism for gravitational wave generation in the early Universe. Light spectator scalar fields during inflation can acquire a blue-tilted power spectrum due to stochastic effects. We show that this effect can lead to large curvature perturbations at small scales (induced by the spectator field fluctuations) while maintaining the observed, slightly red-tilted curvature perturbations at large cosmological scales (induced by the inflaton fluctuations). Along with other observational signatures, such as enhanced dark matter substructure, large curvature perturbations can induce a stochastic gravitational wave background (SGWB). The predicted strength of SGWB in our scenario, $Ω_{\rm GW}h^2 \simeq 10^{-20} - 10^{-15}$, can be observed with future detectors, operating between $10^{-5}$ Hz and 10 Hz. We note that, in order to accommodate the newly reported NANOGrav observation, one could consider the same class of spectator models. At the same time, one would need to go beyond the simple benchmark considered here and consider a regime in which a misalignment contribution is also important. △ Less

Submitted 3 July, 2023; originally announced July 2023.

Comments: 12 pages, 5 figures

Dissipative Inflation via Scalar Production

Authors: Paolo Creminelli, Soubhik Kumar, Borna Salehian, Luca Santoni

Abstract: We describe a new mechanism that gives rise to dissipation during cosmic inflation. In the simplest implementation, the mechanism requires the presence of a massive scalar field with a softly-broken global $U(1)$ symmetry, along with the inflaton field. Particle production in this scenario takes place on parametrically sub-horizon scales, at variance with the case of dissipation into gauge fields.… ▽ More We describe a new mechanism that gives rise to dissipation during cosmic inflation. In the simplest implementation, the mechanism requires the presence of a massive scalar field with a softly-broken global $U(1)$ symmetry, along with the inflaton field. Particle production in this scenario takes place on parametrically sub-horizon scales, at variance with the case of dissipation into gauge fields. Consequently, the backreaction of the produced particles on the inflationary dynamics can be treated in a \textit{local} manner, allowing us to compute their effects analytically. We determine the parametric dependence of the power spectrum which deviates from the usual slow-roll expression. Non-Gaussianities are always sizeable whenever perturbations are generated by the noise induced by dissipation: $f_{\rm NL}^{\rm eq} \gtrsim {O}(10)$. △ Less

Submitted 12 May, 2023; originally announced May 2023.

Comments: 31 pages + appendices, 8 figures

Leptogenesis and Neutrinoless Double Beta Decay in the Scotogenic Hybrid Textures of Neutrino Mass Matrix

Authors: Ankush, Rishu Verma, Sahil Kumar, B. C. Chauhan

Abstract: In our recent work we identify the hybrid textures of neutrino mass matrix which simultaneously account for dark matter (DM) and neutrinoless double beta decay ($0νββ$). We also obtained the bounds on dark matter mass and effective Majorana mass $|M_{ee}|$. In this work we look for those hybrid textures which altogether accounts for DM, $0νββ$ and leptogenesis. We have found correlation of baryon… ▽ More In our recent work we identify the hybrid textures of neutrino mass matrix which simultaneously account for dark matter (DM) and neutrinoless double beta decay ($0νββ$). We also obtained the bounds on dark matter mass and effective Majorana mass $|M_{ee}|$. In this work we look for those hybrid textures which altogether accounts for DM, $0νββ$ and leptogenesis. We have found correlation of baryon asymmetry of universe $Y$ with dark matter mass $M_1$ and effective Majorana mass $|M_{ee}|$. We use experimental bounds on relic density of dark matter ($Ωh^2$) and baryon asymmetry of universe to identify the hybrid textures. We found that out of five hybrid textures which simultaneously satisfies the physics observations of the DM and $0νββ$ only three hybrid textures altogether satisfy the DM, $0νββ$ and leptogenesis. It is interesting to note that these three hybrid textures gives lower bound to the effective Majorana mass $|M_{ee}|$ which can be probed in current and future experiments like SuperNEMO, KamLAND-Zen, NEXT, and nEXO (5 year) have sensitivity reaches of 0.05 eV, 0.045 eV, 0.03 eV, and 0.015 eV, respectively. △ Less

Submitted 21 March, 2023; v1 submitted 20 March, 2023; originally announced March 2023.

Comments: 19 pages, 6 figures, 4 tables, section 2, 3 and 5 restructured

Probing Cosmological Particle Production and Pairwise Hotspots with Deep Neural Networks

Authors: Taegyun Kim, Jeong Han Kim, Soubhik Kumar, Adam Martin, Moritz Münchmeyer, Yuhsin Tsai

Abstract: Particles with masses much larger than the inflationary Hubble scale, $H_I$, can be pair-produced non-adiabatically during inflation. Due to their large masses, the produced particles modify the curvature perturbation around their locations. These localized perturbations eventually give rise to localized signatures on the Cosmic Microwave Background (CMB), in particular, pairwise hotspots (PHS). I… ▽ More Particles with masses much larger than the inflationary Hubble scale, $H_I$, can be pair-produced non-adiabatically during inflation. Due to their large masses, the produced particles modify the curvature perturbation around their locations. These localized perturbations eventually give rise to localized signatures on the Cosmic Microwave Background (CMB), in particular, pairwise hotspots (PHS). In this work, we show that Convolutional Neural Networks (CNN) provide a powerful tool for identifying PHS on the CMB. While for a given hotspot profile a traditional Matched Filter Analysis is known to be optimal, a Neural Network learns to effectively detect the large variety of shapes that can arise in realistic models of particle production. Considering an idealized situation where the dominant background to the PHS signal comes from the standard CMB fluctuations, we show that a CNN can isolate the PHS with $\mathcal{O}(10)\%$ efficiency even if the hotspot temperature is $\mathcal{O}(10)$ times smaller than the average CMB fluctuations. Overall, the CNN search is sensitive to heavy particle masses $M_0/H_I=\mathcal{O}(200)$, and constitutes one of the unique probes of very high energy particle physics. △ Less

Submitted 15 March, 2023; originally announced March 2023.

The Present and Future of QCD

Authors: P. Achenbach, D. Adhikari, A. Afanasev, F. Afzal, C. A. Aidala, A. Al-bataineh, D. K. Almaalol, M. Amaryan, D. Androić, W. R. Armstrong, M. Arratia, J. Arrington, A. Asaturyan, E. C. Aschenauer, H. Atac, H. Avakian, T. Averett, C. Ayerbe Gayoso, X. Bai, K. N. Barish, N. Barnea, G. Basar, M. Battaglieri, A. A. Baty, I. Bautista , et al. (378 additional authors not shown)

Abstract: This White Paper presents the community inputs and scientific conclusions from the Hot and Cold QCD Town Meeting that took place September 23-25, 2022 at MIT, as part of the Nuclear Science Advisory Committee (NSAC) 2023 Long Range Planning process. A total of 424 physicists registered for the meeting. The meeting highlighted progress in Quantum Chromodynamics (QCD) nuclear physics since the 2015… ▽ More This White Paper presents the community inputs and scientific conclusions from the Hot and Cold QCD Town Meeting that took place September 23-25, 2022 at MIT, as part of the Nuclear Science Advisory Committee (NSAC) 2023 Long Range Planning process. A total of 424 physicists registered for the meeting. The meeting highlighted progress in Quantum Chromodynamics (QCD) nuclear physics since the 2015 LRP (LRP15) and identified key questions and plausible paths to obtaining answers to those questions, defining priorities for our research over the coming decade. In defining the priority of outstanding physics opportunities for the future, both prospects for the short (~ 5 years) and longer term (5-10 years and beyond) are identified together with the facilities, personnel and other resources needed to maximize the discovery potential and maintain United States leadership in QCD physics worldwide. This White Paper is organized as follows: In the Executive Summary, we detail the Recommendations and Initiatives that were presented and discussed at the Town Meeting, and their supporting rationales. Section 2 highlights major progress and accomplishments of the past seven years. It is followed, in Section 3, by an overview of the physics opportunities for the immediate future, and in relation with the next QCD frontier: the EIC. Section 4 provides an overview of the physics motivations and goals associated with the EIC. Section 5 is devoted to the workforce development and support of diversity, equity and inclusion. This is followed by a dedicated section on computing in Section 6. Section 7 describes the national need for nuclear data science and the relevance to QCD research. △ Less

Submitted 4 March, 2023; originally announced March 2023.

Comments: QCD Town Meeting White Paper, as submitted to 2023 NSAC LRP committee on Feb. 28, 2023

Journal ref: Nucl.Phys.A 1047 (2024) 122874

Measuring Electroweak Quantum Numbers of Color Sextet Resonances at the LHC

Authors: Soubhik Kumar, Rafiqul Rahaman, Ritesh K. Singh

Abstract: We study the prospect of measuring the electroweak quantum numbers of beyond the Standard Model (SM) color sextet particles that decay into same-sign top quark pairs. Among these particles, the color sextet scalars give rise to top quarks with the same chirality, while the top quarks coming from the color sextet vector would have opposite chirality. This difference gets encoded in the angular dist… ▽ More We study the prospect of measuring the electroweak quantum numbers of beyond the Standard Model (SM) color sextet particles that decay into same-sign top quark pairs. Among these particles, the color sextet scalars give rise to top quarks with the same chirality, while the top quarks coming from the color sextet vector would have opposite chirality. This difference gets encoded in the angular distributions of the bottom quarks and leptons originating from the decays of the top quarks. We utilize this feature and the energy distributions of the final state jets and leptons to distinguish among the three possible color sextet resonances, taking into account various SM background processes at the $13$ TeV LHC. △ Less

Submitted 20 September, 2023; v1 submitted 22 February, 2023; originally announced February 2023.

Comments: 9 pages, 5 captioned figures, and 5 tables

Report number: HRI-RECAPP-2023-01

Journal ref: Physical Review D 108, 055006 (2023)

Exploring the $H_0$ tension and the evidence of dark sector interaction from 2D BAO measurements

Authors: Armando Bernui, Eleonora Di Valentino, William Giarè, Suresh Kumar, Rafael C. Nunes

Abstract: We explore observational constraints on a cosmological model with an interaction between dark energy (DE) and dark matter (DM), using a compilation of 15 measurements of the 2D BAO (i.e., transversal) scale in combination with Planck-CMB data, to explore the parametric space of a class of interacting DE models. We find that 2D BAO measurements can generate different observational constraints compa… ▽ More We explore observational constraints on a cosmological model with an interaction between dark energy (DE) and dark matter (DM), using a compilation of 15 measurements of the 2D BAO (i.e., transversal) scale in combination with Planck-CMB data, to explore the parametric space of a class of interacting DE models. We find that 2D BAO measurements can generate different observational constraints compared to the traditional approach of studying the matter clustering in the 3D BAO measurements. Contrary to the observations for the $Λ$CDM and IDE models when analyzed with Planck-CMB + 3D BAO data, we note that Planck-CMB + 2D BAO data favor high values of the Hubble constant $H_0$. From the joint analysis with Planck-CMB + 2D BAO + Gaussian prior on $H_0$, we find $H_0 = 73.4 \pm 0.88$ km/s/Mpc. We conclude that the $H_0$ tension is solved in the IDE model with strong statistical evidence (more than 3$σ$) for the IDE cosmologies. △ Less

Submitted 1 June, 2023; v1 submitted 15 January, 2023; originally announced January 2023.

Comments: 7 Pages, 2 Figures, 2 Tables. Matches the version published in PRD

Journal ref: Phys. Rev. D 107, 103531 (2023)

Multifield Positivity Bounds for Inflation

Authors: Marat Freytsis, Soubhik Kumar, Grant N. Remmen, Nicholas L. Rodd

Abstract: Positivity bounds represent nontrivial limitations on effective field theories (EFTs) if those EFTs are to be completed into a Lorentz-invariant, causal, local, and unitary framework. While such positivity bounds have been applied in a wide array of physical contexts to obtain useful constraints, their application to inflationary EFTs is subtle since Lorentz invariance is spontaneously broken duri… ▽ More Positivity bounds represent nontrivial limitations on effective field theories (EFTs) if those EFTs are to be completed into a Lorentz-invariant, causal, local, and unitary framework. While such positivity bounds have been applied in a wide array of physical contexts to obtain useful constraints, their application to inflationary EFTs is subtle since Lorentz invariance is spontaneously broken during cosmic inflation. One path forward is to employ a $\textit{Breit parameterization}$ to ensure a crossing-symmetric and analytic S-matrix in theories with broken boosts. We extend this approach to a theory with multiple fields, and uncover a fundamental obstruction that arises unless all fields obey a dispersion relation that is approximately lightlike. We then apply the formalism to various classes of inflationary EFTs, with and without isocurvature perturbations, and employ this parameterization to derive new positivity bounds on such EFTs. For multifield inflation, we also consider bounds originating from the generalized optical theorem and demonstrate how these can give rise to stronger constraints on EFTs compared to constraints from traditional elastic positivity bounds alone. We compute various shapes of non-Gaussianity (NG), involving both adiabatic and isocurvature perturbations, and show how the observational parameter space controlling the strength of NG can be constrained by our bounds. △ Less

Submitted 10 September, 2023; v1 submitted 19 October, 2022; originally announced October 2022.

Comments: 49 pages, 4 figures; v2: clarifications added, journal version

Report number: CERN-TH-2022-160

Journal ref: JHEP 09 (2023) 041

Enhancing Searches for Heavy QCD Axions via Dimuon Final States

Authors: Raymond T. Co, Soubhik Kumar, Zhen Liu

Abstract: Heavy QCD axions are well-motivated extensions of the QCD axion that address the quality problem while still solving the strong CP problem. Owing to the gluon coupling, critical for solving the strong CP problem, these axions can be produced in significant numbers in beam dump and collider environments for axion decay constants as large as PeV, relevant for addressing the axion quality problem. In… ▽ More Heavy QCD axions are well-motivated extensions of the QCD axion that address the quality problem while still solving the strong CP problem. Owing to the gluon coupling, critical for solving the strong CP problem, these axions can be produced in significant numbers in beam dump and collider environments for axion decay constants as large as PeV, relevant for addressing the axion quality problem. In addition, if these axions have leptonic couplings, they can give rise to long-lived decay into lepton pairs, in particular, dominantly into muons above the dimuon threshold and below the GeV scale in a broad class of axion models. Considering existing constraints, primarily from rare meson decays, we demonstrate that current and future neutrino facilities and long-lived particle searches have the potential to probe significant parts of the heavy QCD axion parameter space via dimuon final states. △ Less

Submitted 5 October, 2022; originally announced October 2022.

Comments: 19 pages + references, 6 figures

Report number: UMN-TH-4201/22, FTPI-MINN-22/25

Higher-Order Cumulants and Correlation Functions of Proton Multiplicity Distributions in $\sqrt{s_{\mathrm{NN}}}$ = 3 GeV Au+Au Collisions at the RHIC STAR Experiment

Authors: STAR Collaboration, M. S. Abdallah, B. E. Aboona, J. Adam, L. Adamczyk, J. R. Adams, J. K. Adkins, I. Aggarwal, M. M. Aggarwal, Z. Ahammed, D. M. Anderson, E. C. Aschenauer, J. Atchison, V. Bairathi, W. Baker, J. G. Ball Cap, K. Barish, R. Bellwied, P. Bhagat, A. Bhasin, S. Bhatta, J. Bielcik, J. Bielcikova, J. D. Brandenburg, X. Z. Cai , et al. (349 additional authors not shown)

Abstract: We report a measurement of cumulants and correlation functions of event-by-event proton multiplicity distributions from fixed-target Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 3 GeV measured by the STAR experiment. Protons are identified within the rapidity ($y$) and transverse momentum ($p_{\rm T}$) region $-0.9 < y<0$ and $0.4 < p_{\rm T} <2.0 $ GeV/$c$ in the center-of-mass frame. A systematic a… ▽ More We report a measurement of cumulants and correlation functions of event-by-event proton multiplicity distributions from fixed-target Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 3 GeV measured by the STAR experiment. Protons are identified within the rapidity ($y$) and transverse momentum ($p_{\rm T}$) region $-0.9 < y<0$ and $0.4 < p_{\rm T} <2.0 $ GeV/$c$ in the center-of-mass frame. A systematic analysis of the proton cumulants and correlation functions up to sixth-order as well as the corresponding ratios as a function of the collision centrality, $p_{\rm T}$, and $y$ are presented. The effect of pileup and initial volume fluctuations on these observables and the respective corrections are discussed in detail. The results are compared to calculations from the hadronic transport UrQMD model as well as a hydrodynamic model. In the most central 5\% collisions, the value of proton cumulant ratio $C_4/C_2$ is negative, drastically different from the values observed in Au+Au collisions at higher energies. Compared to model calculations including Lattice QCD, a hadronic transport model, and a hydrodynamic model, the strong suppression in the ratio of $C_4/C_2$ at 3 GeV Au+Au collisions indicates an energy regime dominated by hadronic interactions. △ Less

Submitted 22 February, 2023; v1 submitted 24 September, 2022; originally announced September 2022.

Comments: 25 pages, 20 figures, 4 tables

Journal ref: Phys. Rev. C 107, 024908(2023)

Report of the Topical Group on Cosmic Probes of Fundamental Physics for for Snowmass 2021

Authors: Rana X. Adhikari, Luis A. Anchordoqui, Ke Fang, B. S. Sathyaprakash, Kirsten Tollefson, Tiffany R. Lewis, Kristi Engel, Amin Aboubrahim, Ozgur Akarsu, Yashar Akrami, Roberto Aloisio, Rafael Alves Batista, Mario Ballardini, Stefan W. Ballmer, Ellen Bechtol, David Benisty, Emanuele Berti, Simon Birrer, Alexander Bonilla, Richard Brito, Mauricio Bustamante, Robert Caldwell, Vitor Cardoso, Sukanya Chakrabarti, Thomas Y. Chen , et al. (96 additional authors not shown)

Abstract: Cosmic Probes of Fundamental Physics take two primary forms: Very high energy particles (cosmic rays, neutrinos, and gamma rays) and gravitational waves. Already today, these probes give access to fundamental physics not available by any other means, helping elucidate the underlying theory that completes the Standard Model. The last decade has witnessed a revolution of exciting discoveries such as… ▽ More Cosmic Probes of Fundamental Physics take two primary forms: Very high energy particles (cosmic rays, neutrinos, and gamma rays) and gravitational waves. Already today, these probes give access to fundamental physics not available by any other means, helping elucidate the underlying theory that completes the Standard Model. The last decade has witnessed a revolution of exciting discoveries such as the detection of high-energy neutrinos and gravitational waves. The scope for major developments in the next decades is dramatic, as we detail in this report. △ Less

Submitted 23 September, 2022; originally announced September 2022.

Comments: Report of theTopical Group on Cosmic Probes of Fundamental Physics, for the U.S. decadal Particle Physics Planning Exercise (Snowmass 2021)

Beam Energy Dependence of Triton Production and Yield Ratio ($\mathrm{N}_t \times \mathrm{N}_p/\mathrm{N}_d^2$) in Au+Au Collisions at RHIC

Authors: STAR Collaboration, M. I. Abdulhamid, B. E. Aboona, J. Adam, J. R. Adams, G. Agakishiev, I. Aggarwal, M. M. Aggarwal, Z. Ahammed, A. Aitbaev, I. Alekseev, D. M. Anderson, A. Aparin, S. Aslam, J. Atchison, G. S. Averichev, V. Bairathi, W. Baker, J. G. Ball Cap, K. Barish, P. Bhagat, A. Bhasin, S. Bhatta, I. G. Bordyuzhin, J. D. Brandenburg , et al. (333 additional authors not shown)

Abstract: We report the triton ($t$) production in mid-rapidity ($|y| <$ 0.5) Au+Au collisions at $\sqrt{s_\mathrm{NN}}$= 7.7--200 GeV measured by the STAR experiment from the first phase of the beam energy scan at the Relativistic Heavy Ion Collider (RHIC). The nuclear compound yield ratio ($\mathrm{N}_t \times \mathrm{N}_p/\mathrm{N}_d^2$), which is predicted to be sensitive to the fluctuation of local ne… ▽ More We report the triton ($t$) production in mid-rapidity ($|y| <$ 0.5) Au+Au collisions at $\sqrt{s_\mathrm{NN}}$= 7.7--200 GeV measured by the STAR experiment from the first phase of the beam energy scan at the Relativistic Heavy Ion Collider (RHIC). The nuclear compound yield ratio ($\mathrm{N}_t \times \mathrm{N}_p/\mathrm{N}_d^2$), which is predicted to be sensitive to the fluctuation of local neutron density, is observed to decrease monotonically with increasing charged-particle multiplicity ($dN_{ch}/dη$) and follows a scaling behavior. The $dN_{ch}/dη$ dependence of the yield ratio is compared to calculations from coalescence and thermal models. Enhancements in the yield ratios relative to the coalescence baseline are observed in the 0\%-10\% most central collisions at 19.6 and 27 GeV, with a significance of 2.3$σ$ and 3.4$σ$, respectively, giving a combined significance of 4.1$σ$. The enhancements are not observed in peripheral collisions or model calculations without critical fluctuation, and decreases with a smaller $p_{T}$ acceptance. The physics implications of these results on the QCD phase structure and the production mechanism of light nuclei in heavy-ion collisions are discussed. △ Less

Submitted 18 May, 2023; v1 submitted 16 September, 2022; originally announced September 2022.

Comments: 6 pages, 4 figures, Supplemental Material: http://link.aps.org/supplemental/10.1103/PhysRevLett.130.202301

Journal ref: Phys. Rev. Lett. 130, 202301 (2023)

Search for the Chiral Magnetic Effect in Au+Au collisions at $\sqrt{s_{_{\rm{NN}}}}=27$ GeV with the STAR forward Event Plane Detectors

Authors: STAR Collaboration, B. E. Aboona, J. Adam, L. Adamczyk, J. R. Adams, I. Aggarwal, M. M. Aggarwal, Z. Ahammed, D. M. Anderson, E. C. Aschenauer, J. Atchison, V. Bairathi, W. Baker, J. G. Ball Cap, K. Barish, R. Bellwied, P. Bhagat, A. Bhasin, S. Bhatta, J. Bielcik, J. Bielcikova, J. D. Brandenburg, X. Z. Cai, H. Caines, M. Calderón de la Barca Sánchez , et al. (347 additional authors not shown)

Abstract: A decisive experimental test of the Chiral Magnetic Effect (CME) is considered one of the major scientific goals at the Relativistic Heavy-Ion Collider (RHIC) towards understanding the nontrivial topological fluctuations of the Quantum Chromodynamics vacuum. In heavy-ion collisions, the CME is expected to result in a charge separation phenomenon across the reaction plane, whose strength could be s… ▽ More A decisive experimental test of the Chiral Magnetic Effect (CME) is considered one of the major scientific goals at the Relativistic Heavy-Ion Collider (RHIC) towards understanding the nontrivial topological fluctuations of the Quantum Chromodynamics vacuum. In heavy-ion collisions, the CME is expected to result in a charge separation phenomenon across the reaction plane, whose strength could be strongly energy dependent. The previous CME searches have been focused on top RHIC energy collisions. In this Letter, we present a low energy search for the CME in Au+Au collisions at $\sqrt{s_{_{\rm{NN}}}}=27$ GeV. We measure elliptic flow scaled charge-dependent correlators relative to the event planes that are defined at both mid-rapidity $|η|<1.0$ and at forward rapidity $2.1 < |η|<5.1$. We compare the results based on the directed flow plane ($Ψ_1$) at forward rapidity and the elliptic flow plane ($Ψ_2$) at both central and forward rapidity. The CME scenario is expected to result in a larger correlation relative to $Ψ_1$ than to $Ψ_2$, while a flow driven background scenario would lead to a consistent result for both event planes. In 10-50\% centrality, results using three different event planes are found to be consistent within experimental uncertainties, suggesting a flow driven background scenario dominating the measurement. We obtain an upper limit on the deviation from a flow driven background scenario at the 95\% confidence level. This work opens up a possible road map towards future CME search with the high statistics data from the RHIC Beam Energy Scan Phase-II. △ Less

Submitted 19 April, 2023; v1 submitted 7 September, 2022; originally announced September 2022.

Comments: main: 16 pages, 5 figures; supplementary material: 2 pages, 1 figure

Exploring dark $Z_d$-boson in future Large Hadron-electron collider

Authors: Ashok Goyal, Mukesh Kumar, Satendra Kumar, Rafiqul Rahaman

Abstract: The interaction between the dark $U(1)_d$ sector with the visible Standard Model (SM) sector takes place through the kinetic mixing between the dark photon $U(1)_d$ field $Z_d^μ$ and the SM $U(1)_Y$ gauge field $B_μ$. After the electroweak and $U(1)_d$ symmetry breaking, the dark photon $Z_d^μ$ acquires a mass and mixes with the SM neutral vector boson $Z_μ$. This mixing leads to parity-violating… ▽ More The interaction between the dark $U(1)_d$ sector with the visible Standard Model (SM) sector takes place through the kinetic mixing between the dark photon $U(1)_d$ field $Z_d^μ$ and the SM $U(1)_Y$ gauge field $B_μ$. After the electroweak and $U(1)_d$ symmetry breaking, the dark photon $Z_d^μ$ acquires a mass and mixes with the SM neutral vector boson $Z_μ$. This mixing leads to parity-violating coupling between the $Z_d^μ$ and SM. The coupling between the dark photon and SM can be explored in low energy phenomenology as well as in collider experiments. The Lorentz structure of dark photon interaction with SM fermions is explored in the proposed high energy future Large Hadron-electron collider, which would provide efficient energy and a clean environment using cross-section and asymmetries associated with polarisation observable of the dark photon in leptons decay. A $χ^2$-analysis is performed to compare the strength of various variables for both the charge- and neutral-current processes. Based on this analysis, $90\%$ confidence level (C.L.) contours in the $ε$-$m_{Z_d}$ and $ε$-$g_V$ plane are obtained to put limits on the $Z_d^μ$ mass up to $100$ GeV, coupling strength $ε$ and on the Lorentz structure of dark photon coupling with the SM fermions ($g_V$) at $\sqrt{s} \approx 1.3$ TeV. △ Less

Submitted 13 February, 2023; v1 submitted 7 September, 2022; originally announced September 2022.

Comments: published version in EPJC

Journal ref: Eur. Phys. J. C (2023) 83:132

Dark Grand Unification in the Axiverse: Decaying Axion Dark Matter and Spontaneous Baryogenesis

Authors: Joshua W. Foster, Soubhik Kumar, Benjamin R. Safdi, Yotam Soreq

Abstract: The quantum chromodynamics axion with a decay constant near the Grand Unification (GUT) scale has an ultralight mass near a neV. We show, however, that axion-like particles with masses near the keV - PeV range with GUT-scale decay constants are also well motivated in that they naturally arise from axiverse theories with dark non-abelian gauge groups. We demonstrate that the correct dark matter abu… ▽ More The quantum chromodynamics axion with a decay constant near the Grand Unification (GUT) scale has an ultralight mass near a neV. We show, however, that axion-like particles with masses near the keV - PeV range with GUT-scale decay constants are also well motivated in that they naturally arise from axiverse theories with dark non-abelian gauge groups. We demonstrate that the correct dark matter abundance may be achieved by the heavy axions in these models through the misalignment mechanism in combination with a period of early matter domination from the long-lived dark glueballs of the same gauge group. Heavy axion dark matter may decay to two photons, yielding mono-energetic electromagnetic signatures that may be detectable by current or next-generation space-based telescopes. We project the sensitivity of next-generation telescopes including $\textit {Athena,}$ AMEGO, and e-ASTROGAM to such decaying axion dark matter. If the dark sector contains multiple confining gauge groups, then the observed primordial baryon asymmetry may also be achieved in this scenario through spontaneous baryogenesis. We present explicit orbifold constructions where the dark gauge groups unify with the SM at the GUT scale and axions emerge as the fifth components of dark gauge fields with bulk Chern-Simons terms. △ Less

Submitted 22 August, 2022; originally announced August 2022.

Comments: 52 pages, 6 figures

Report number: MIT-CTP/5458

Pion, kaon, and (anti-)proton production in U+U Collisions at $\sqrt{s_{NN}}$ = 193 GeV measured with the STAR detector

Authors: STAR Collaboration, M. S. Abdallah, B. E. Aboona, J. Adam, J. R. Adams, J. K. Adkins, G. Agakishiev, I. Aggarwal, M. M. Aggarwal, Z. Ahammed, A. Aitbaev, I. Alekseev, D. M. Anderson, A. Aparin, J. Atchison, G. S. Averichev, V. Bairathi, W. Baker, J. G. Ball Cap, K. Barish, P. Bhagat, A. Bhasin, S. Bhatta, I. G. Bordyuzhin, J. D. Brandenburg , et al. (330 additional authors not shown)

Abstract: We present the first measurements of transverse momentum spectra of $π^{\pm}$, $K^{\pm}$, $p(\bar{p})$ at midrapidity ($|y| < 0.1$) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV with the STAR detector at the Relativistic Heavy Ion Collider (RHIC). The centrality dependence of particle yields, average transverse momenta, particle ratios and kinetic freeze-out parameters are discussed. The results… ▽ More We present the first measurements of transverse momentum spectra of $π^{\pm}$, $K^{\pm}$, $p(\bar{p})$ at midrapidity ($|y| < 0.1$) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV with the STAR detector at the Relativistic Heavy Ion Collider (RHIC). The centrality dependence of particle yields, average transverse momenta, particle ratios and kinetic freeze-out parameters are discussed. The results are compared with the published results from Au+Au collisions at $\sqrt{s_{NN}} =$ 200 GeV in STAR. The results are also compared to those from A Multi Phase Transport (AMPT) model. △ Less

Submitted 11 February, 2023; v1 submitted 1 August, 2022; originally announced August 2022.

Comments: 17 pages, 14 figures and 7 tables; Replaced with the updated version published in Physical Review C

Journal ref: Phys. Rev. C 107 (2023) 024901

First Constraints on Heavy QCD Axions with a Liquid Argon Time Projection Chamber using the ArgoNeuT Experiment

Authors: ArgoNeuT Collaboration, R. Acciarri, C. Adams, B. Baller, V. Basque, F. Cavanna, R. T. Co, R. S. Fitzpatrick, B. Fleming, P. Green, R. Harnik, K. J. Kelly, S. Kumar, K. Lang, I. Lepetic, Z. Liu, X. Luo, K. F. Lyu, O. Palamara, G. Scanavini, M. Soderberg, J. Spitz, A. M. Szelc, W. Wu, T. Yang

Abstract: We present the results of a search for heavy QCD axions performed by the ArgoNeuT experiment at Fermilab. We search for heavy axions produced in the NuMI neutrino beam target and absorber decaying into dimuon pairs, which can be identified using the unique capabilities of ArgoNeuT and the MINOS near detector. This decay channel is motivated by a broad class of heavy QCD axion models that address t… ▽ More We present the results of a search for heavy QCD axions performed by the ArgoNeuT experiment at Fermilab. We search for heavy axions produced in the NuMI neutrino beam target and absorber decaying into dimuon pairs, which can be identified using the unique capabilities of ArgoNeuT and the MINOS near detector. This decay channel is motivated by a broad class of heavy QCD axion models that address the strong CP and axion quality problems with axion masses above the dimuon threshold. We obtain new constraints at a 95\% confidence level for heavy axions in the previously unexplored mass range between 0.2-0.9 GeV, for axion decay constants around tens of TeV. △ Less

Submitted 24 April, 2023; v1 submitted 18 July, 2022; originally announced July 2022.

Comments: 6 pages, 4 figures. Version accepted by PRL

Report number: FERMILAB-PUB-22-527-ND-T

Journal ref: Phys. Rev. Lett. 130 (2023) 22, 221802

Constraining mass of the graviton with GW170817

Authors: Andrey A. Shoom, Sumit Kumar, N. V. Krishnendu

Abstract: We consider the massive graviton phenomenological model based on the graviton's dispersion terms included into phase of gravitational wave's waveform. Such model was already considered in many works but it was based on a single leading-order dispersion term only. Here we derive a relation between relativistic gravitons emission and absorption time intervals computed up to ${\cal O}(γ^{-6})$, where… ▽ More We consider the massive graviton phenomenological model based on the graviton's dispersion terms included into phase of gravitational wave's waveform. Such model was already considered in many works but it was based on a single leading-order dispersion term only. Here we derive a relation between relativistic gravitons emission and absorption time intervals computed up to ${\cal O}(γ^{-6})$, where $γ$ is the Lorentz factor. Including the dispersion terms into the phase of gravitational wave's waveform results in two non-GR parameters of the $1^{st}$ and the $-2^{nd}$ post-Newtonian orders whose posteriors are used to put a constraint on the graviton's rest mass. We use the TaylorF2 waveform model to analyse the event GW170817 and report the following $95\%$-confidence upper bounds on the graviton's rest mass: $m^{Low\,Spin}_{g}\leq1.305\times10^{-54}$g and $m^{High\,Spin}_{g}\leq2.996\times10^{-54}$g for the high and low spin priors. △ Less

Submitted 30 May, 2022; originally announced May 2022.

Comments: 9 pages, 6 figures

Faint light of old neutron stars and detectability at the James Webb Space Telescope

Authors: Shiuli Chatterjee, Raghuveer Garani, Rajeev Kumar Jain, Brijesh Kanodia, M. S. N. Kumar, Sudhir K. Vempati

Abstract: Isolated ideal neutron stars (NS) of age $>10^9$ yrs exhaust thermal and rotational energies and cool down to temperatures below $\mathcal{O}(100)$ K. Accretion of particle dark matter (DM) by such NS can heat them up through kinetic and annihilation processes. This increases the NS surface temperature to a maximum of $\sim 2550$ K in the best case scenario. The maximum accretion rate depends on t… ▽ More Isolated ideal neutron stars (NS) of age $>10^9$ yrs exhaust thermal and rotational energies and cool down to temperatures below $\mathcal{O}(100)$ K. Accretion of particle dark matter (DM) by such NS can heat them up through kinetic and annihilation processes. This increases the NS surface temperature to a maximum of $\sim 2550$ K in the best case scenario. The maximum accretion rate depends on the DM ambient density and velocity dispersion, and on the NS equation of state and their velocity distributions. Upon scanning over these variables, we find that the effective surface temperature varies at most by $\sim 40\%$. Black body spectrum of such warm NS peak at near infrared wavelengths with magnitudes in the range potentially detectable by the James Webb Space Telescope (JWST). Using the JWST exposure time calculator, we demonstrate that NS with surface temperatures $\gtrsim 2400$ K, located at a distance of 10\,pc can be detected through the F150W2 (F322W2) filters of the NIRCAM instrument at SNR\,$\gtrsim 10$ (5) within 24 hours of exposure time. Independently of DM, an observation of NS with surface temperatures $\gtrsim 2500$ K will be a formative step towards testing the minimal cooling paradigm during late evolutionary stages. △ Less

Submitted 13 July, 2023; v1 submitted 10 May, 2022; originally announced May 2022.

Comments: v2: 7+3 pages, 2+2 figures; discussion on background modelling added. Matches journal version

Journal ref: Phys.Rev.D 108 (2023) 2, L021301