In this talk, we present the package GravitinoPack that calculates decays of unstable supersymmetric particles, involving gravitinos in the final or initial state. If the gravitino is the dark matter particle and therefore stable, the package calculates the decays of the lightest neutralino, and the lighter stau or stop NLSP into the gravitino LSP and one or two Standard Model particles. On the other hand, assuming that the gravitino is unstable, GravitinoPack calculates all its two-body and the three-body decay widths to the neutralino LSP and Standard Model particles. Since all these decays, involving the gravitino, are of gravitational nature, the lifetime of the decaying particle can be of the order of seconds are more, hence called "late decays". The precise knowledge of all these partial decay widths enables the user to apply accurately the relevant cosmological constraints.

We present results from a new calculation of two- and three-body decays of the gravitino including all possible Feynman graphs. The work is done in the R-parity conserving Minimal Supersymmetric Standard Model, assuming that the gravitino is unstable. We calculate the full width for the three-body decay process Gravitino -> \chi_1 X Y, where \chi_1 is the lightest neutralino that plays the role of the lightest supersymmetric particle, and X, Y are standard model particles. For calculating the squared amplitudes we use the packages FeynArts and FormCalc, with appropriate extensions to accommodate the gravitino with spin 3/2. Our code treats automatically the intermediate exchanged particles that are on their mass shell, using the narrow width approximation. This method considerably simplifies the complexity of the evaluation of the total gravitino width and stabilizes it numerically. At various points of the supersymmetric parameter space, we compare the full two-body and the three-body results with approximations used in the past. In addition, we discuss the relative sizes of the full two-body and the three-body gravitino decay widths. Our results are obtained using the package GravitinoPack, that will be publicly available soon, including various handy features.