Francesco Costa
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Research Subjects
Francesco Costa has a total of 22 co-authors. Here are the top co-authors by number of shared publications:
- Jinsu Kim
University of Shanghai for Science and Technology
5 shared publications - Oleg Lebedev
Department of Physics and Helsinki Institute of Physics, Gustaf Hällströmin katu 2a, FI-00014 Helsinki, Finland
4 shared publications - Sarif Khan
Aliah University
4 shared publications - Catarina Cosme
Ottawa-Carleton Institute for Physics Carleton University 1125 Colonel By Drive K1S 5B6 Ottawa Ontario Canada
2 shared publications - Jaime Hoefken Zink
National Centre for Nuclear Research
2 shared publications - Michele Lucente
Università di Bologna
2 shared publications
Latest Publications
- Open Access
- 11 September 2026
- Published by Springer Nature in Journal of High Energy Physics
- Vol. 2026 (9) , 150
Decaying spin-3/2 dark matter from baryon number violation
A bstract: We explore an uncharted corner of dark matter phenomenology: non-supersymmetric spin-3/2 dark matter with baryon number violating interactions. In an effective field theory description, we identify the leading baryonic portal between the spin-3/2 state and Standard Model quarks and show that it can account for the observed dark matter abundance through UV freeze-in and Boltzmann-suppressed freeze-in, while the freeze-out region is completely excluded. The resulting phenomenology is distinctive, with relic production controlled by the competition between baryon-violating single-particle processes and baryon-conserving pair production. We map the viable parameter space against indirect detection, direct detection, and LHC monojet bounds, finding strong complementarity between these probes and especially stringent limits when production and decay are tied to the same operator. We also present a dark QCD-like ultraviolet completion in which the spin-3/2 particle arises as a composite baryon, naturally generating the effective interactions and mitigating the main theoretical obstacles of elementary higher-spin states. This framework opens a novel and testable connection between baryonic portals, confining dark sectors, and higher-spin dark matter.
- 11 September 2026
- Published by Springer Nature in Journal of High Energy Physics
- Vol. 2026 (9) , 150
Open AccessDecaying spin-3/2 dark matter from baryon number violation
A bstract: We explore an uncharted corner of dark matter phenomenology: non-supersymmetric spin-3/2 dark matter with baryon number violating interactions. In an effective field theory description, we identify the leading baryonic portal between the spin-3/2 state and Standard Model quarks and show that it can account for the observed dark matter abundance through UV freeze-in and Boltzmann-suppressed freeze-in, while the freeze-out region is completely excluded. The resulting phenomenology is distinctive, with relic production controlled by the competition between baryon-violating single-particle processes and baryon-conserving pair production. We map the viable parameter space against indirect detection, direct detection, and LHC monojet bounds, finding strong complementarity between these probes and especially stringent limits when production and decay are tied to the same operator. We also present a dark QCD-like ultraviolet completion in which the spin-3/2 particle arises as a composite baryon, naturally generating the effective interactions and mitigating the main theoretical obstacles of elementary higher-spin states. This framework opens a novel and testable connection between baryonic portals, confining dark sectors, and higher-spin dark matter.
- Open Access
- 22 July 2026
- Published by Springer Nature in The European Physical Journal C
- Vol. 86 (7) , 851
ELENA: a software for fast and precise computation of first order phase transitions and gravitational waves production in particle physics models
We present (EvaLuator of tunnElliNg Actions), an open-source Python package designed to compute the full evolution of first-order phase transitions in the early Universe generated by particle physics models, taking into account several refinements that go beyond commonly assumed simplifications. The core of is based on a vectorized implementation of the tunnelling potential formalism, making it the first and only public code to compute finite-temperature tunnelling actions within this framework. This allows for a fast computation of the finite-temperature tunnelling action. This, in turn, enables the sampling of the full range of temperatures where two phases coexist and the use of integral expressions that track the complete evolution of the transition, providing a comprehensive picture of it. In addition, provides all the tools to compute the resulting stochastic gravitational waves spectrum, allowing for the full chain of computations – from the Lagrangian parameter inputs to the final gravitational waves spectrum – in a fast and self-contained implementation.
- 22 July 2026
- Published by Springer Nature in The European Physical Journal C
- Vol. 86 (7) , 851
Open AccessELENA: a software for fast and precise computation of first order phase transitions and gravitational waves production in particle physics models
We present (EvaLuator of tunnElliNg Actions), an open-source Python package designed to compute the full evolution of first-order phase transitions in the early Universe generated by particle physics models, taking into account several refinements that go beyond commonly assumed simplifications. The core of is based on a vectorized implementation of the tunnelling potential formalism, making it the first and only public code to compute finite-temperature tunnelling actions within this framework. This allows for a fast computation of the finite-temperature tunnelling action. This, in turn, enables the sampling of the full range of temperatures where two phases coexist and the use of integral expressions that track the complete evolution of the transition, providing a comprehensive picture of it. In addition, provides all the tools to compute the resulting stochastic gravitational waves spectrum, allowing for the full chain of computations – from the Lagrangian parameter inputs to the final gravitational waves spectrum – in a fast and self-contained implementation.
- Open Access
- 19 June 2026
- Published by Springer Nature in Journal of High Energy Physics
- Vol. 2026 (6)
Gravitational scalar production with a generic reheating scenario
A bstract: Gravitational production of decoupled scalars during inflationary and post-inflationary phases is efficient and can lead to over-production. We study this production with various reheating scenarios such as a generic power-law inflaton potential V inf ∝ ϕ k as well as a multi-stage reheating scenario. We derive constraints on the scalar self-interaction coupling λ s , the mass m s , and coefficients of quantum gravity-induced operators. We find that the constraints depend sensitively on the reheating dynamics. Our analysis demonstrates that universal gravity effects do not necessarily spoil the predictivity of non-thermal dark matter scenarios with k < 4 and low reheating temperatures, as an extended reheating phase dilutes gravitationally-produced relics. For k > 4, on the other hand, the relic abundance is enhanced during the reheating phase, leading to stringent constraints on the scalar. In multi-stage reheating, we show that the enhancement/dilution effect of subsequent reheating phases factorises.
- 19 June 2026
- Published by Springer Nature in Journal of High Energy Physics
- Vol. 2026 (6)
Open AccessGravitational scalar production with a generic reheating scenario
A bstract: Gravitational production of decoupled scalars during inflationary and post-inflationary phases is efficient and can lead to over-production. We study this production with various reheating scenarios such as a generic power-law inflaton potential V inf ∝ ϕ k as well as a multi-stage reheating scenario. We derive constraints on the scalar self-interaction coupling λ s , the mass m s , and coefficients of quantum gravity-induced operators. We find that the constraints depend sensitively on the reheating dynamics. Our analysis demonstrates that universal gravity effects do not necessarily spoil the predictivity of non-thermal dark matter scenarios with k < 4 and low reheating temperatures, as an extended reheating phase dilutes gravitationally-produced relics. For k > 4, on the other hand, the relic abundance is enhanced during the reheating phase, leading to stringent constraints on the scalar. In multi-stage reheating, we show that the enhancement/dilution effect of subsequent reheating phases factorises.
- Open Access
- 1 October 2025
- Published by Elsevier in Nuclear Physics B
- Vol. 1020
From oversimplified to overlooked: The case for exploring rich dark sectors
- 1 October 2025
- Published by Elsevier in Nuclear Physics B
- Vol. 1020
Open AccessFrom oversimplified to overlooked: The case for exploring rich dark sectors
- 1 December 2024
- Published by IOP Publishing in Journal of Cosmology and Astroparticle Physics
- Vol. 2024 (12)
High Frequency Gravitational Wave bounds from galactic neutron stars
High-Frequency Gravitational Waves (HFGWs) constitute a unique window on the early Universe as well as exotic astrophysical objects. While the current gravitational wave experiments are more dedicated to the low frequency regime, the graviton conversion into photons in a strong magnetic field constitutes a powerful tool to probe HFGWs. In this paper, we show that neutron stars, due to their extreme magnetic field, are a perfect laboratory to study the conversion of HFGWs into photons. Using realistic models for the galactic neutron star population, we calculate for the first time the expected photon flux induced by the conversion of an isotropic stochastic gravitational wave background in the magnetosphere of the ensemble of neutron stars present in the Milky Way. We compare this photon flux to the observed one from several telescopes and derive upper limits on the stochastic gravitational wave background in the frequency range 108 Hz–1025 Hz. We find our limits to be competitive in the frequency range 108 Hz–1012 Hz.
- 1 December 2024
- Published by IOP Publishing in Journal of Cosmology and Astroparticle Physics
- Vol. 2024 (12)
High Frequency Gravitational Wave bounds from galactic neutron stars
High-Frequency Gravitational Waves (HFGWs) constitute a unique window on the early Universe as well as exotic astrophysical objects. While the current gravitational wave experiments are more dedicated to the low frequency regime, the graviton conversion into photons in a strong magnetic field constitutes a powerful tool to probe HFGWs. In this paper, we show that neutron stars, due to their extreme magnetic field, are a perfect laboratory to study the conversion of HFGWs into photons. Using realistic models for the galactic neutron star population, we calculate for the first time the expected photon flux induced by the conversion of an isotropic stochastic gravitational wave background in the magnetosphere of the ensemble of neutron stars present in the Milky Way. We compare this photon flux to the observed one from several telescopes and derive upper limits on the stochastic gravitational wave background in the frequency range 108 Hz–1025 Hz. We find our limits to be competitive in the frequency range 108 Hz–1012 Hz.
- Open Access
- 4 July 2024
- Published by Springer Nature in Journal of High Energy Physics
- Vol. 2024 (7) , 1-19
Higgs portal dark matter freeze-in at stronger coupling: observational benchmarks
We study freeze-in production of Higgs portal dark matter (DM) at temperatures far below the dark matter mass. The temperature of the Standard Model (SM) thermal bath may have never been high such that dark matter production via thermal emission has been Boltzmann-suppressed. This allows for a significant coupling between the Higgs field and DM, which is being probed by the direct DM detection experiments and invisible Higgs decay searches at the LHC. We delineate the corresponding parameter space in the Higgs portal framework with dark matter of spin 0, 1/2 and 1.
- 4 July 2024
- Published by Springer Nature in Journal of High Energy Physics
- Vol. 2024 (7) , 1-19
Open AccessHiggs portal dark matter freeze-in at stronger coupling: observational benchmarks
We study freeze-in production of Higgs portal dark matter (DM) at temperatures far below the dark matter mass. The temperature of the Standard Model (SM) thermal bath may have never been high such that dark matter production via thermal emission has been Boltzmann-suppressed. This allows for a significant coupling between the Higgs field and DM, which is being probed by the direct DM detection experiments and invisible Higgs decay searches at the LHC. We delineate the corresponding parameter space in the Higgs portal framework with dark matter of spin 0, 1/2 and 1.
- Open Access
- 1 June 2024
- Published by IOP Publishing in Journal of Cosmology and Astroparticle Physics
- Vol. 2024 (06)
Temperature evolution in the Early Universe and freeze-in at stronger coupling
Dark matter freeze-in at stronger coupling is operative when the Standard Model (SM) bath temperature never exceeds the dark matter mass. An attractive feature of this scenario is that it can be probed by direct detection experiments as well as at the LHC. In this work, we show how the mechanism can be realized in a simple UV complete framework, emphasizing the role of the maximal temperature of the SM thermal bath. We demonstrate that the maximal temperature can coincide with the reheating temperature or be close to it such that dark matter production is always Boltzmann-suppressed. This possibility is realized, for example, if the inflaton decays primarily into feebly interacting right-handed neutrinos, which subsequently generate the SM thermal bath. In this case, the SM sector temperature remains constant over cosmological times prior to reheating.
- 1 June 2024
- Published by IOP Publishing in Journal of Cosmology and Astroparticle Physics
- Vol. 2024 (06)
Open AccessTemperature evolution in the Early Universe and freeze-in at stronger coupling
Dark matter freeze-in at stronger coupling is operative when the Standard Model (SM) bath temperature never exceeds the dark matter mass. An attractive feature of this scenario is that it can be probed by direct detection experiments as well as at the LHC. In this work, we show how the mechanism can be realized in a simple UV complete framework, emphasizing the role of the maximal temperature of the SM thermal bath. We demonstrate that the maximal temperature can coincide with the reheating temperature or be close to it such that dark matter production is always Boltzmann-suppressed. This possibility is realized, for example, if the inflaton decays primarily into feebly interacting right-handed neutrinos, which subsequently generate the SM thermal bath. In this case, the SM sector temperature remains constant over cosmological times prior to reheating.
- Open Access
Freeze-in at stronger coupling
Predictivity of many non-thermal dark matter (DM) models is marred by the gravitational production background. This problem is ameliorated in models with lower reheating temperature $T_R$, which allows for dilution of gravitationally produced relics. We study the freeze-in dark matter production mechanism in the thermal bath with the electroweak scale temperature. The process is Boltzmann-suppressed if the dark matter mass is above $T_R$. In this case, the coupling to the thermal bath has to be significant to account for the observed dark matter relic density. As a result, the direct DM detection experiments already probe such freeze-in models, excluding significant parts of parameter space. The forthcoming experiments will explore this framework further, extending to lower couplings and higher reheating temperatures.Open AccessFreeze-in at stronger coupling
Predictivity of many non-thermal dark matter (DM) models is marred by the gravitational production background. This problem is ameliorated in models with lower reheating temperature $T_R$, which allows for dilution of gravitationally produced relics. We study the freeze-in dark matter production mechanism in the thermal bath with the electroweak scale temperature. The process is Boltzmann-suppressed if the dark matter mass is above $T_R$. In this case, the coupling to the thermal bath has to be significant to account for the observed dark matter relic density. As a result, the direct DM detection experiments already probe such freeze-in models, excluding significant parts of parameter space. The forthcoming experiments will explore this framework further, extending to lower couplings and higher reheating temperatures. - Open Access
- 22 April 2024
- Published by American Physical Society (APS) in Physical Review D
- Vol. 109 (7) , 075038
Freeze-in at stronger coupling
The predictivity of many nonthermal dark matter (DM) models is marred by the gravitational production background. Nevertheless, if the reheating temperature is low, the gravitationally produced relics can be diluted. We study the freeze-in dark matter production mechanism at temperatures below the dark matter mass. In this case, the coupling to the thermal bath has to be significant to account for the observed dark matter relic density. As a result, the direct DM detection experiments already probe such freeze-in models, excluding significant parts of parameter space.- 22 April 2024
- Published by American Physical Society (APS) in Physical Review D
- Vol. 109 (7) , 075038
Open AccessFreeze-in at stronger coupling
The predictivity of many nonthermal dark matter (DM) models is marred by the gravitational production background. Nevertheless, if the reheating temperature is low, the gravitationally produced relics can be diluted. We study the freeze-in dark matter production mechanism at temperatures below the dark matter mass. In this case, the coupling to the thermal bath has to be significant to account for the observed dark matter relic density. As a result, the direct DM detection experiments already probe such freeze-in models, excluding significant parts of parameter space. - Open Access
Neutrino mixing sum rules and the Littlest Seesaw
In this work, we study the neutrino mixing sum rules arising from discrete symmetries, and the class of Littlest Seesaw (LS) neutrino models. These symmetry based approaches all offer predictions for the cosine of the leptonic CP phase $\cos δ$ in terms of the mixing angles, $θ_{13}$, $θ_{12}$, $θ_{23}$, while the LS models also predict the sine of the leptonic CP phase $\sin δ$ as well as making other predictions. In particular we study the \textit{solar} neutrino mixing sum rules, arising from charged lepton corrections to Tri-bimaximal (TB), Bi-maximal (BM), Golden Ratios (GRs) and Hexagonal (HEX) neutrino mixing, and \textit{atmospheric} neutrino mixing sum rules, arising from preserving one of the columns of these types of mixing, for example the first or second column of the TB mixing matrix (TM1 or TM2), and confront them with an up-to-date global fit of the neutrino oscillation data. We show that some mixing sum rules, for example an \textit{atmospheric} neutrino mixing sum rule arising from a version of neutrino Golden Ratio mixing (GRa1), are already excluded at 3$σ$, and determine the remaining models allowed by the data. We also consider the more predictive LS models (which obey the TM1 sum rules and offer further predictions) based on constrained sequential dominance CSD($n$) with $n\approx 3$. We compare for the first time the three cases $n=2.5$, $n=3$ and $n=1+\sqrt{6}\approx 3.45$ which are favoured by theoretical models, using a new type of analysis to accurately predict the observables $θ_{12}$, $θ_{23}$ and $δ$. We study all the above approaches, \textit{solar} and \textit{atmospheric} mixing sum rules and LS models, together so that they may be compared, and to give an up to date analysis of the predictions of all of these possibilities, when confronted with the most recent global fits.Open AccessNeutrino mixing sum rules and the Littlest Seesaw
In this work, we study the neutrino mixing sum rules arising from discrete symmetries, and the class of Littlest Seesaw (LS) neutrino models. These symmetry based approaches all offer predictions for the cosine of the leptonic CP phase $\cos δ$ in terms of the mixing angles, $θ_{13}$, $θ_{12}$, $θ_{23}$, while the LS models also predict the sine of the leptonic CP phase $\sin δ$ as well as making other predictions. In particular we study the \textit{solar} neutrino mixing sum rules, arising from charged lepton corrections to Tri-bimaximal (TB), Bi-maximal (BM), Golden Ratios (GRs) and Hexagonal (HEX) neutrino mixing, and \textit{atmospheric} neutrino mixing sum rules, arising from preserving one of the columns of these types of mixing, for example the first or second column of the TB mixing matrix (TM1 or TM2), and confront them with an up-to-date global fit of the neutrino oscillation data. We show that some mixing sum rules, for example an \textit{atmospheric} neutrino mixing sum rule arising from a version of neutrino Golden Ratio mixing (GRa1), are already excluded at 3$σ$, and determine the remaining models allowed by the data. We also consider the more predictive LS models (which obey the TM1 sum rules and offer further predictions) based on constrained sequential dominance CSD($n$) with $n\approx 3$. We compare for the first time the three cases $n=2.5$, $n=3$ and $n=1+\sqrt{6}\approx 3.45$ which are favoured by theoretical models, using a new type of analysis to accurately predict the observables $θ_{12}$, $θ_{23}$ and $δ$. We study all the above approaches, \textit{solar} and \textit{atmospheric} mixing sum rules and LS models, together so that they may be compared, and to give an up to date analysis of the predictions of all of these possibilities, when confronted with the most recent global fits.