We compute the dark matter relic densities of neutralinos and axions in a supersymmetric model with a gauged anomalous U(1) symmetry, kinetically mixed with $U(1)_Y$ of hypercharge. The model is a variant of the USSM (the U(1) extended NMSSM), containing an extra U(1) symmetry and an extra singlet in the superpotential respect to the MSSM, where gauge invariance is restored by Peccei-Quinn interactions using a Stuckelberg multiplet. This approach introduces an axion (Im b) and a saxion (Re b) in the spectrum and generates an axino component for the neutralino. The Stuckelberg axion (Im b) develops a physical component (the gauged axion) after electroweak symmetry breaking. We classify all the interactions of the Lagrangian and perform a complete simulation study of the spectrum, determining the neutralino relic densities using micrOMEGAs. We discuss the phenomenological implications of the model analyzing mass values for the axion from the milli-eV to the MeV region. The possible scenarios that we analyze are significantly constrained by a combination of WMAP data, the exclusion limits from direct axion searches and the veto on late entropy release at the time of nucleosynthesis.

We illustrate the structure and the main phenomenological features of a supersymmetric model (the USSM-A) built following a bottom-up approach and containing an anomalous abelian gauge symmetry. This model supports a gauged axion in its spectrum and provides a generalization of the global (supersymmetric) Peccei-Quinn construction. Complete simulations of the neutralino relic density are performed. Bounds from CAST and WMAP, combined with dark matter simulations, provide significant constraints on the scale of the interactions between the axion and the gauge fields.

We analyze the most salient cosmological features of axions in extensions of the Standard Model with a gauged anomalous extra U(1) symmetry. The model is built by imposing the constraint of gauge invariance in the anomalous effective action, which is extended with Wess-Zumino counterterms. These generate axion-like interactions of the axions to the gauge fields and a gauged shift symmetry. The scalar sector is assumed to acquire a non-perturbative potential after inflation, at the electroweak phase transition, which induces a mixing of the Stuckelberg field of the model with the scalars of the electroweak sector, and at the QCD phase transition. We discuss the possible mechanisms of sequential misalignments which could affect the axions of these models, and generated, in this case, at both transitions. We compute the contribution of these particles to dark matter, quantifying their relic densities as a function of the Stuckelberg mass. We also show that models with a single anomalous U(1) in general do not account for the dark energy, due to the presence of mixed U(1)-SU(3) anomalies.