We present supersymmetric scenarios with gravitino LSP and stau NLSP in the case of a non-standard model of cosmology with the addition of a dark component in the pre-BBN era. In the context of the standard model of cosmology, gravitino LSP has drawn quite some attention as it is a good candidate for dark matter. It is produced in scattering processes during reheating after inflation and from the decay of the stau. With a long lifetime, the stau decays during Big Bang Nucleosynthesis. It is strongly constrained by the abundance of light elements but can however address the known "BBN lithium problem". It requires fairly massive staus $\mstau \gtrsim 1 \unit{TeV}$ and puts an upper bound on the reheating temperature $T_R \simeq 10^7\unit{GeV}$ which does not satisfy the requirements for thermal leptogenesis. For the non-standard cosmological scenario, the reheating temperature bound can be strongly relaxed $T_R\gg 10^9 \unit{GeV}$ and the lithium-7 problem solved with a stau typical mass of $\mstau \sim 600-700 \unit{GeV}$ and down to $\sim 400 \unit{GeV}$ with a very important dark component that could enable possible production and detection at the LHC.

We re-consider the gravitino as dark matter in the framework of the Constrained MSSM. We include several recently suggested improvements on: (i) the thermal production of gravitinos, (ii) the calculation of the hadronic spectrum from NLSP decay and (iii) the BBN calculation including stau bound-state effects. In most cases we find an upper bound on the reheating temperature $\treh\lsim {\rm a few}\times 10^7\gev$ from over-production of $^6{\rm Li}$ from bound state effects. We also find an upper limit on the stau lifetime of $3\times10^4 \sec$, which is nearly an order of magnitude larger than the simple limit $5\times 10^3 \sec$ often used to avoid the effect of bound-state catalysis. The bound on $\treh$ is relaxed to $\lsim {\rm a few}\times 10^8\gev$ when we use a more conservative bound on ${^6}{\rm Li}/{^7}{\rm Li}$, in which case a new region at small stau mass at $\sim 100\gev$ and much longer lifetimes opens up. Such a low stau mass region can be easily tested at the LHC.

Supersymmetric extensions of the standard model of particle physics assuming the gravitino to be the lightest supersymmetric particle (LSP), and with the next-to-LSP decaying to the gravitino during Big Bang nucleosynthesis, are analyzed. Particular emphasis is laid on their potential to solve the "Li7 problem", an apparent factor 2-4 overproduction of Li7 in standard Big Bang nucleosynthesis (BBN), their production of cosmologically important amounts of Li6, as well as the resulting gravitino dark matter densities in these models. The study includes several improvements compared to prior studies. Heavy gravitinos in the constrained minimal supersymmetric standard model (CMMSM) are reanalyzed, whereas light gravitinos in gauge-mediated supersymmetry breaking scenarios (GMSB) are studied for the first time. It is confirmed that decays of NLSP staus to heavy gravitinos, while producing all the dark matter, may at the same time resolve the Li7 problem. For NLSP decay times ~ 1000 sec, such scenarios also lead to cosmologically important Li6 (and possibly Be9) abundances. However, as such scenarios require heavy > 1 TeV staus they are likely not testable at the LHC. It is found that decays of NLSP staus to light gravitinos may lead to significant Li6 (and Be9) abundances, whereas NLSP neutralinos decaying into light gravitinos may solve the Li7 problem. Though both scenarios are testable at the LHC they may not lead to the production of the bulk of the dark matter. A section of the paper outlines particle properties required to significantly reduce the Li7 abundance, and/or enhance the Li6 (and possibly Be9) abundances, by the decay of an arbitrary relic particle.