Dark Sector Research Articles

Shedding light on dark sectors with gravitational waves

The nature of dark matter remains one of the greatest unsolved mysteries in elementary particle physics. It might well be that the dark matter particle belongs to a dark sector completely secluded or extremely weakly coupled to the visible sector. We demonstrate that gravitational waves arising from first-order phase transitions in the early Universe can be used to look for signatures of dark sector models connected to neutron physics. This introduces a new connection between gravitational-wave physics and nuclear physics experiments. Focusing on two particular extensions of the Standard Model with dark U(1) and SU(2) gauge groups constructed to address the neutron lifetime puzzle, we show how those signatures can be searched for in future gravitational-wave and astrometry experiments. Published by the American Physical Society 2024

Fractional Einstein–Gauss–Bonnet Scalar Field Cosmology

Our paper introduces a new theoretical framework called the Fractional Einstein–Gauss–Bonnet scalar field cosmology, which has important physical implications. Using fractional calculus to modify the gravitational action integral, we derived a modified Friedmann equation and a modified Klein–Gordon equation. Our research reveals non-trivial solutions associated with exponential potential, exponential couplings to the Gauss–Bonnet term, and a logarithmic scalar field, which are dependent on two cosmological parameters, m and α0=t0H0 and the fractional derivative order μ. By employing linear stability theory, we reveal the phase space structure and analyze the dynamic effects of the Gauss–Bonnet couplings. The scaling behavior at some equilibrium points reveals that the geometric corrections in the coupling to the Gauss–Bonnet scalar can mimic the behavior of the dark sector in modified gravity. Using data from cosmic chronometers, type Ia supernovae, supermassive Black Hole Shadows, and strong gravitational lensing, we estimated the values of m and α0, indicating that the solution is consistent with an accelerated expansion at late times with the values α0=1.38±0.05, m=1.44±0.05, and μ=1.48±0.17 (consistent with Ωm,0=0.311±0.016 and h=0.712±0.007), resulting in an age of the Universe t0=19.0±0.7 [Gyr] at 1σ CL. Ultimately, we obtained late-time accelerating power-law solutions supported by the most recent cosmological data, and we proposed an alternative explanation for the origin of cosmic acceleration other than ΛCDM. Our results generalize and significantly improve previous achievements in the literature, highlighting the practical implications of fractional calculus in cosmology.

Heavy neutrino as dark matter in a neutrinophilic U(1) model

We study the prospect of heavy singlet neutrinos as a dark matter (DM) candidate within a neutrinophilic U(1) model, where the Standard Model (SM) is extended with a U(1) gauge symmetry, and neutrino mass and oscillation parameters are explained through an inverse see-saw mechanism. The lightest of the heavy neutrinos plays the role of the DM, while the newly introduced scalars and the extra gauge boson Z′\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$Z'$$\\end{document} act as mediators between the dark sector and the SM sector. We show the range of model parameters where this DM candidate can be accommodated in the Weakly Interacting Massive Particle (WIMP) or Feebly Interacting Massive Particle (FIMP) scenario. The observed DM relic density (DMRD) is achieved via the new gauge boson and singlet scalar portals in the WIMP scenario, whereas within the FIMP scenario, these two particles assume a distinct yet pivotal role in generating the observed DMRD.

Muon-induced baryon number violation

The search for charged-lepton flavor violation in muon capture on nuclei is a powerful probe of heavy new physics. A smoking gun signal for μ→e conversion is a monochromatic electron with energy almost equal to the muon mass. We show that light new physics can mimic this signature and that it can also lead to electrons above the μ→e signal peak. A concrete example of such light new physics is μ−-nucleon annihilation into a light dark sector, which can produce an energetic e− as well as e+e− by-products. Due to the size of the muon mass, the exotic muon capture process can be kinematically allowed, while the otherwise stringent constraints, e.g., from proton decay, are kinematically forbidden. We also discuss other relevant constraints, including those from the stability of nuclei and muon capture in the interior of neutron stars. Published by the American Physical Society 2024

Interacting Dark Energy and Its Implications for Unified Dark Sector

Interacting Dark Energy and Its Implications for Unified Dark Sector

Probing interacting dark sector with the next generation of gravitational-wave detectors

Probing interacting dark sector with the next generation of gravitational-wave detectors

Generalising axion-like particle as the curvaton: sourcing primordial density perturbation and non-Gaussianities

We investigate the non-perturbatively generated axion-like particle (ALP) potential, involving fermions in the dark sector that couple to the ALP, in an early universe cosmological inflationary stage with the ALP being a spectator field. The potential here deviates from the standard cosine nature due to the presence of the two dark sector fermion masses mu\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$m_u$$\\end{document} and md\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$m_d$$\\end{document} which couple to the ALP. The ALP is a spectator field during inflation but it starts to oscillate and dominates the energy density of the universe after inflation ends, thereby sourcing isocurvature perturbations, while standard curvature fluctuations from the inflaton are assumed to be sub-dominant. Subsequently the ALP decays converting the isocurvature perturbations to adiabatic perturbations thereby acting as the origin of the primordial density perturbations. We identify the parameter space involving the axion decay constant fϕ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$f_{\\phi }$$\\end{document}, scale of confinement Λ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\Lambda $$\\end{document}, ALP mass m and the masses of the fermions, mu\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$m_u$$\\end{document} and md\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$m_d$$\\end{document} where it can satisfactorily behave as the curvaton and source the observed primordial density perturbation. We also predict local non-Gaussianity signals for bi-spectrum and tri-spectrum fNL\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$f_{NL}$$\\end{document} and gNL\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$g_{NL}$$\\end{document}, as a function of the ratio mu/md\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$m_u/m_d$$\\end{document}, which are within the allowed range in the latest Planck observations and are detectable with future observations. Particularly we find that the value of fNL\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$f_{NL}$$\\end{document} and gNL\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$g_{NL}$$\\end{document} are dependent on the ratio of mu\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$m_u$$\\end{document} and md\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$m_d$$\\end{document}: fNL\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$f_{NL}$$\\end{document} is more or less positive for all scenarios except mu=md\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$m_u = m_d$$\\end{document} and gNL\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$g_{NL}$$\\end{document} is always positive irrespective of the ratio between mu\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$m_u$$\\end{document} and md\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$m_d$$\\end{document}. The results of our analysis in the limit mu=md\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$m_u = m_d$$\\end{document} resembles vanilla curvaton scenario while in the limit mu≫md\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$m_u \\gg m_d$$\\end{document} resembles pure axion cosine potential. This is because the model has the unique feature where the extra matter in the dark sector determines the nature and shape of the potential.

Illuminating the dark: mono-γ signals at NA62

Dipole interactions between dark sector states or between a Standard Model particle and a dark state can efficiently be searched for via high-intensity fixed-target facilities. We propose to look for the associated mono-γ signature at the NA62 experiment running in beam dump mode. Focusing on models of dipole inelastic Dark Matter and active-sterile neutrino dipole interactions, we compute the corresponding expected sensitivities finding promising prospects for discovery already with ~ 1017 proton-on-target, corresponding to the present accumulated dataset.

Singlet–doublet Dirac Fermion dark matter from Peccei–Quinn symmetry

Weakly interacting massive particles (WIMPs) and axions are arguably the most compelling dark matter (DM) candidates in the literature. Here, we consider a model where the Peccei–Quinn (PQ) symmetry solves the strong CP problem, generates radiatively Dirac neutrino masses and gives origin to a multicomponent dark sector. Specifically, scotogenic Dirac neutrino masses arise at one-loop level. The lightest fermionic mediator acts as the second DM candidate due to a residual [Formula: see text] symmetry resulting from the PQ symmetry breaking. The WIMP DM component resembles the well-known singlet–doublet fermion DM. While the lower WIMP dark mass region is usually excluded, our model reopens that portion of the parameter space (for DM masses below [Formula: see text][Formula: see text]GeV). Therefore, we perform a phenomenological analysis that addresses the constraints from direct searches of DM, neutrino oscillation data and charged lepton flavor violating processes. The model can be tested in future facilities where neutrino telescopes search for DM annihilation into Standard Model particles.

Redshift-space distortions corner interacting dark energy

Despite the fact that the ΛCDM model has been highly successful over the last few decades in providing an accurate fit to a broad range of cosmological and astrophysical observations, different intriguing tensions and anomalies emerged at various statistical levels. Given the fact that the dark energy and the dark matter sectors remain unexplored, the answer to some of the tensions may rely on modifications of these two dark sectors. This manuscript explores the important role of the growth of structure in constraining non-standard cosmologies. In particular, we focus on the interacting dark energy (IDE) scenario, where dark matter and dark energy interact non-gravitationally. We aim to place constraints on the phenomenological parameters of these alternative models, by considering different datasets related to a number of cosmological measurements, to achieve a complementary analysis. A special emphasis is devoted to redshift-space distortion measurements (RSD), whose role in constraining beyond the standard paradigm models has not been recently highlighted. These observations indeed have a strong constraining power, rendering all parameters to their ΛCDM canonical values, and therefore leaving little room for the IDE models explored here.

Cold darkogenesis: Dark matter and baryon asymmetry in light of the PTA signal

We build upon the intriguing possibility that the recently reported nano-Hz gravitational wave signal by Pulsar Timing Array (PTA) experiments is sourced by a strong first-order phase transition from a nearly conformal dark sector. The phase transition has to be strongly supercooled to explain the signal amplitude, while the critical temperature has to be in the O(GeV) range, as dictated by the peak frequency of the gravitational wave spectrum. However, the resulting strong supercooling exponentially dilutes away any pre-existing baryon asymmetry and dark matter, calling for a new paradigm of their productions. We then develop a mechanism of cold darkogenesis that generates a dark asymmetry during the phase transition from the textured dark SU(2)D Higgs field. This dark asymmetry is transferred to the visible sector via neutron portal interactions, resulting in the observed baryon asymmetry. Furthermore, the mechanism naturally leads to the correct abundance of asymmetric dark matter, with self-interaction of the scale that is of the right order to solve the diversity problem in galactic rotation curves. Collider searches for mono-jets and dark matter direct detection experiments can dictate the viability of the model.

Constraints on light dark sector particles from lifetime difference of heavy neutral mesons

Heavy meson decays with missing energy in the final state offer interesting avenues to search for light invisible new physics such as dark matter (DM). In this context, we show that such new physics (NP) interactions also affect lifetime difference in neutral meson-antimeson mixing. We consider general dimension-six effective quark interactions involving a pair of DM particles and calculate their contributions to lifetime difference in beauty and charm meson systems. We use the latest data on mixing observables to constrain the relevant effective operators. We find that lifetime differences provide novel and complementary flavor constraints compared to those obtained from heavy meson decays. Published by the American Physical Society 2024

Dark-Sector Search via Pion-Produced η and η^{'} Mesons Decaying Invisibly in the NA64h Detector.

We present the first results from a proof-of-concept search for dark sectors via invisible decays of pseudoscalar η and η^{'} mesons in the NA64h experiment at the CERN SPS. Our novel technique uses the charge-exchange reaction of 50GeV π^{-} on nuclei of an active target as the source of neutral mesons. The η,η^{'}→invisible events would exhibit themselves via a striking signature-the complete disappearance of the incoming beam energy in the detector. No evidence for such events has been found with 2.9×10^{9} pions on target accumulated during one day of data taking. This allows us to set a stringent limit on the branching ratio Br(η^{'}→invisible)<2.1×10^{-4} improving the current bound by a factor of ≃3. We also set a limit on Br(η→invisible)<1.1×10^{-4} comparable with the existing one. These results demonstrate the great potential of our approach and provide clear guidance on how to enhance and extend the sensitivity for dark sector physics from future searches for invisible neutral meson decays.

Unveiling desert region in inert doublet model assisted by Peccei-Quinn symmetry

The Inert Higgs Doublet model (IDM), assisted by Peccei-Quinn (PQ) symmetry, offers a simple but natural framework of a dark sector that accommodates Weakly Interacting Massive Particle (WIMP) and axion as dark matter components. Spontaneous breaking of U(1)PQ symmetry, which was originally proposed as an elegant solution to the strong charge-parity (CP) problem, also ensures the stability of WIMP through a residual ℤ2 symmetry. Interestingly, additional fields necessitated by PQ symmetry further enrich the dark sector. These include a scalar field proprietor for axion DM and a vector-like quark (VLQ) that acts as a portal for the dark sector through Yukawa interactions. Moreover, this combination of the axion and WIMP components satisfies the observed DM relic density and reopens the phenomenologically exciting region of the IDM parameter space where the WIMP mass falls between 100 - 550 GeV. We investigate the model-independent pair production of VLQs exploring this region at the Large Hadron Collider (LHC), incorporating the effects of next-to-leading order (NLO) QCD corrections. After production, each VLQ decays into a top or bottom quark accompanied by an inert scalar, a consequence of the residual ℤ2 symmetry. Utilising relevant observables with a leptonic search channel and employing multivariate analysis, we demonstrate the ability of this analysis to exclude a significant portion of the parameter space with an integrated luminosity of 300 fb−1.

Constraints on Metastable Dark Energy Decaying into Dark Matter

We revisit the proposal that an energy transfer from dark energy into dark matter can be described in field theory by a first order phase transition. We analyze a metastable dark energy model proposed in the literature, using updated constraints on the decay time of a metastable dark energy from recent data. The results of our analysis show no prospects for potentially observable signals that could distinguish this scenario from the ΛCDM. We analyze, for the first time, the process of bubble nucleation in this model, showing that such model would not drive a complete transition to a dark matter dominated phase even in a distant future. Nevertheless, the model is not excluded by the latest data and we confirm that the mass of the dark matter particle that would result from such a process corresponds to the mass of an axion-like particle, which is currently one of the best motivated dark matter candidates. We argue that extensions to this model, possibly with additional couplings, still deserve further attention as it could provide an interesting and viable description for an interacting dark sector scenario based in a single scalar field.

Relic density and temperature evolution of a light dark sector

We have developed a set of four fully coupled Boltzmann equations to precisely determine the relic density and temperature of dark matter by including three distinct sectors: dark matter, light scalar, and standard model sectors. The intricacies of heat transfer between dark matter (DM) and the standard model sector through a light scalar particle are explored, inspired by stringent experimental constraints on the scalar-Higgs mixing angle and the DM-scalar coupling. Three distinct sectors emerge prior to DM freeze-out, requiring fully coupled Boltzmann equations to accurately compute relic density. Investigation of forbidden, resonance, and secluded DM scenarios demonstrates significant deviations between established methods and the novel approach with fully coupled Boltzmann equations. Despite increased computational demands, this emphasizes the need for improved precision in relic density calculations, underlining the importance of incorporating these equations in comprehensive analyses. Published by the American Physical Society 2024

No νs is Good News

The baryon acoustic oscillation (BAO) analysis from the first year of data from the Dark Energy Spectroscopic Instrument (DESI), when combined with data from the cosmic microwave background (CMB), has placed an upper-limit on the sum of neutrino masses, ∑mν< 70 meV (95%). In addition to excluding the minimum sum associated with the inverted hierarchy, the posterior is peaked at ∑mν = 0 and is close to excluding even the minumum sum, 58 meV at 2σ. In this paper, we explore the implications of this data for cosmology and particle physics. The sum of neutrino mass is determined in cosmology from the suppression of clustering in the late universe. Allowing the clustering to be enhanced, we extended the DESI analysis to ∑mν< 0 and find ∑mν =160±90 meV (68%), and that the suppression of power from the minimum sum of neutrino masses is excluded at 99% confidence. We show this preference for negative masses makes it challenging to explain the result by a shift of cosmic parameters, such as the optical depth or matter density. We then show how a result of ∑mν = 0 could arise from new physics in the neutrino sector, including decay, cooling, and/or time-dependent masses. These models are consistent with current observations but imply new physics that is accessible in a wide range of experiments. In addition, we discuss how an apparent signal with ∑mν< 0 can arise from new long range forces in the dark sector or from a primordial trispectrum that resembles the signal of CMB lensing.

One-point statistics in various cosmic environments in the presence of massive neutrinos

Abstract Studying the structures (halos and galaxies) within the cosmic environments (void, sheet, filament, and node) where they reside is an ongoing attempt in cosmological studies. The link between the properties of structures and the cosmic environments may help to unravel the nature of the dark sector of the Universe. In this paper, we study the cosmic web environments from the spatial pattern perspective in the context of ΛCDM and νΛCDM as an example of an extension to the vanilla model. To do this, we use the T-web classification method and classify the cosmic environments for the catalogues from the gevolution N-body simulations for ΛCDM and νΛCDM cosmology. Then, we compute the first nearest neighbour cumulative distribution function, spherical contact cumulative distribution function, and J-function for every cosmic environment. In the context of the standard model, the results indicate that these functions can differentiate the various cosmic environments. In association with distinguishing between extensions of the standard model of cosmologies, these functions within the cosmic environment seem beneficial as a complementary probe.

Composite dark matter with forbidden annihilation

A dark matter model based on QCD-like SU(Nc) gauge theory with electroweakly interacting dark quarks is discussed. Assuming the dark quark mass m is smaller than the dynamical scale Λd ~ 4πfd, the main component of the dark matter is the lightest G-parity odd dark pion associated with chiral symmetry breaking in the dark sector. We show that nonzero dark quark mass induces the universal mass contribution to both G-parity odd and even pions, and their masses tend to be degenerate. As a result, dark pion annihilation into heavier G-parity even dark pion also affects the dark matter relic abundance. Thus, our setup naturally accommodates forbidden dark matter scenario and realizes heavy dark matter whose mass is O\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{O} $$\\end{document}(1–100) TeV, which is different from conventional electroweakly interacting dark matter such as minimal dark matter. We also discuss CP-violation from the θ-term in the dark gauge sector and find that the predicted size of the electron electric dipole moment can be as large as ∼ 10−32e cm.

Conformal Killing cosmology: Geometry, dark sector, growth of structures, and a big rip

Conformal Killing cosmology: Geometry, dark sector, growth of structures, and a big rip