Abstract: We propose a minimal gauged U(1)(B-L) extension of the minimal supersymmetric Standard Model (MSSM) which resolves the cosmological moduli problem via thermal inflation, and realizes late-time Affleck-Dine leptogensis so as to generate the right amount of baryon asymmetry at the end of thermal inflation. The present relic density of dark matter can be explained by sneutrinos, MSSM neutralinos, axinos, or axions. Cosmic strings from U(1)(B-L) breaking are very thick, and so the expected stochastic gravitational wave background from cosmic string loops has a spectrum different from the one in the conventional Abelian-Higgs model, as would be distinguishable at least at LISA and DECIGO. The characteristic spectrum is due to a flat potential, and may be regarded as a hint of supersymmetry. Combined with the resolution of moduli problem, the expected signal of gravitational waves constrains the U(1)(B-L) breaking scale to be O(10(12-13)) GeV. Interestingly, our model provides a natural possibility for explaining the observed ultra-high-energy cosmic rays thanks to the fact that the core width of strings in our scenario is very large, allowing a large enhancement of particle emissions from the cusps of string loops. Condensation of LHu flat-direction inside of string cores arises inevitably and can also be the main source of the ultra-high-energy cosmic rays accompanied by ultra-high-energy lightest supersymmetric particles.

Abstract: We use the latest constraints on the population of stellar origin binary black holes (SOBBH) from LIGO/Virgo/KAGRA (LVK) observations, to estimate the stochastic gravi-tational wave background (SGWB) they generate in the frequency band of LISA. In order to account for the faint and distant binaries, which contribute the most to the SGWB, we extend the merger rate at high redshift assuming that it tracks the star formation rate. We adopt different methods to compute the SGWB signal: we perform an analytical evaluation, we use Monte Carlo sums over the SOBBH population realisations, and we account for the role of the detector by simulating LISA data and iteratively removing the resolvable signals until only the confusion noise is left. The last method allows the extraction of both the expected SGWB and the number of resolvable SOBBHs. Since the latter are few for signal-to-noise ratio thresholds larger than five, we confirm that the spectral shape of the SGWB in the LISA band agrees with the analytical prediction of a single power law. We infer the probability dis-tribution of the SGWB amplitude from the LVK GWTC-3 posterior of the binary population model: at the reference frequency of 0.003 Hz it has an interquartile range of h(2 Omega)GW(f = 3 x 10(-3) Hz) is an element of [5.65, 11.5] x 10(-13), in agreement with most previous estimates. We then perform a MC analysis to assess LISA's capability to detect and characterise this signal. Ac-counting for both the instrumental noise and the galactic binaries foreground, with four years of data, LISA will be able to detect the SOBBH SGWB with percent accuracy, narrowing down the uncertainty on the amplitude by one order of magnitude with respect to the range of possible amplitudes inferred from the population model. A measurement of this signal by LISA will help to break the degeneracy among some of the population parameters, and pro-vide interesting constraints, in particular on the redshift evolution of the SOBBH merger rate.

Abstract: Cosmological simulations play a pivotal role in understanding the properties of the dark matter (DM) distribution in both galactic and galaxy -cluster environments. The characterization of DM structures is crucial for informing indirect DM searches, aiming at the detection of the annihilation (or decay) products of DM particles. A fundamental quantity in these analyses is the astrophysical J -factor. In the DM phenomenology community, J -factors are typically computed through the semi -analytical modelling of the DM mass distribution, which is affected by large uncertainties. With the scope of addressing and possibly reducing these uncertainties, we present NAJADS, a self-contained framework to derive the DM J -factor directly from the raw simulations data. We show how this framework can be used to compute all -sky maps of the J -factor, automatically accounting for the complex 3D structure of the simulated halos and for the boosting of the signal due to the density fluctuations along the line of sight. After validating our code, we present a proof -of -concept application of NAJADS to a realistic halo from the IllustrisTNG suite, and exploit it to make a thorough comparison between our numerical approach and traditional semi -analytical methods. JCAP02(2024)001