Despite its successful predictions, the Standard Model (SM) suffers from a major drawback as it contains massless neutrinos to all orders in the perturbation theory. Even after including the non-perturbative effects, this problem persists. This is in contradiction with the experimental results that suggest that the neutrinos have non-zero masses and oscillations. The best-fit values at error level for these neutrino oscillation parameters in the three-flavor framework are summarised as follows 
Thus, the oscillation experiments indicate that at least some neutrinos must be massive. However, the above relations do not provide the overall scale of masses which means that other methods must be employed to understand the neutrino mass spectrum. One way to obtain meaningful bounds on the absolute scale for the neutrinos is to look for kinematic effects that can be present as a con- sequence of their non-zero masses in the tritium -decay
. Two groups from Mainz  and from Troitsk  , respectively, have reported on the bounds of and . Also, the upcoming KATRIN experiment  is expected to produce results at a sensitivity of about 0.3 eV, which will further narrow down the scale of the neutrino spectrum. Another way to probe the neutrino mass scale is via studies of the lepton number violating neutrinoless double -decay
 . Several groups such as Heidelberg-Moscow  and IGEX   collaborations conducted experiments with , while the more recent CUORICINO experiment  used to test the lepton number conservation. The best upper bounds on the decay lifetimes are presently provided by CUORICINO (which is still running), whose results are translated to
for the neutrino mass. Note that the large range is due to the uncertainty in the nuclear matrix elements. Upcoming experiments like CUORE    , GERDA   and MAJORANA   are expected to further improve these results with projected sensitivity of about 0.05 eV. Finally, it must be mentioned that some of the strongest bounds on the overall scale for neutrino masses come from cosmology. The studies of the data from the Wilkinson Microwave Anisotropy Probe (WMAP) and the Sloan Digital Sky Survey (SDSS) have deduced that the sum of neutrino masses (three species assumed) is constrained by   and 1.6 eV   .
On the other hand, the LhCb reported recently a deviation of of the measured ratios of the branching fractions in the individual lepton flavour model with respect to the SM in the low invariant mass region given by  . In this range, is defined as
Thus, the experimental results from  put new numerical constraints on the scalar and pseudoscalar couplings. As observed in   , the low invariant mass range excludes the resonant regions and thus improving the theoretical predictions. After these interesting results, the CMS collaboration published an intriguing deviation from the SM in the channel, in the mass region . No significant deviation was observed in the channel   . The ATLAS measured an excess1
with respect to SM predictions in the production of di-electroweak gauge bosons VV (where V = W; Z) that decay hadronically  . The ATLAS and CMS collaborations have recently presented the results of di-photon resonance searches2 in early Run II of data     .
Beside the neutrino masses and the results from Atlas, LhCb and CMS, there are other features of the SM that require an explanation from a more funda- mental point of view:
1) The coupling constants do not meet at a single definite value  .
2) The hierarchy problem  .
3) The naturalness or fine tuning problem  .
4) The large number of parameters    .
Also, it is expected that in a fundamental theory, the gravity take a natural place alongside the other three fundamental interactions.
One promising class of theories that could solve the problems of the SM is formed by the supersymmetric extensions of the SM based on a postulated fundamental symmetry between the bosons and the fermions. The model from this class that contains a minimum number of physical states and interactions is the Minimal Supersymmetric Standard Model (MSSM)       3. The MSSM suffers from the -problem which is the generation of a coupling in the term of the order of the electro-weak scale. The - problem is solved by SSM proposed in  4 which represents a modifi- cation of the MSSM by introducing new Yukawa interactions that generates light neutrino masses, as we will present at Sec. (3.1).
The aim of the present paper is to propose a modification SSM that can explain the recent data from Atlas, CMS and LHCb by introducing new interactions among the leptons with gauginos and higgsinos while the masses of charginos and neutralinos are left unchanged. Since the new model has an explicit broken -parity and lepton number, there are flat directions of the superpotential that can generate the cosmological inflation. Also, the matter anti-matter asymmetry could be obtained from the letogenesis mechanism.
This paper is organized as follows. In order to make the paper self-contained, we review in Section 2 the MSSM and establish our notations. In Section 3 we present a model that generalizes the MSSM. Next, we calculate all flat directions of this model and show that it can generate a viable leptogenesis mechanism5. Also, we show how this model can explain the data from CMS and LHCb. The last section is devoted to conclusions.
2. Review of the SMM
The Minimal Supersymmetric Standard Model (MSSM) is the supersymmetric extension of the SM that contains a minimal number of states and interactions    . It aims at providing a general frame for the solving of the hierarchy problem, for the stabilization of the weak scale, for the unification of the coupling constants and for addressing the dark matter issues, among other things. The model has the gauge symmetry extended by the supersymmetry to include the supersymmetric partners of the SM fields which have spins that differ by as required by the supersymmetric algebra     . Since the SM fermions are left-handed and right-handed and they transform differently under , and groups, the particles of the MSSM must belong to chiral or gauge supermultiplets. The degrees of freedom are grouped in gauge superfields for gauge bosons and left-handed chiral superfields for spinors.
The chiral supermultiplet    contains three families of left-handed (right-handed) quarks and , . Here, the numbers in parenthesis refers to the quantum numbers, respectively and refers to the generation index (or flavor indices) and we neglected the color indices. Also, we use the notation in the anti-right-chiral superfield according to  . The model contains three families of leptons , respectively. The Higgs boson has spin 0, therefore it must belong to a chiral supermultiplet. However, in this case a single Higgs boson cannot provide mass for all quarks that have different weak isospin charges . Therefore, the MSSM contains two left-handed chiral superfields for Higgs fields      . The particle content of each chiral superfield given above is presented in the Table 1 and Table 2 below.
The gauge supermultiplets are described by three vector superfilds
, where , with and
. The particle content in each vector superfield is presented in the Table 3.
The supersymetric Lagrangian of the MSSM is given by
Table 1. Particle content in the left-chiral superfields in MSSM, is flavour index
Table 2. Particle content in the anti-right-chiral superfields in MSSM, is flavour index .
Table 3. Particle content in the vector superfields in MSSM.
The Lagrangian defined in the Equation (4) contains contributions from all sectors of the model
and the terms have the following explicit form
here, and (with ) are the generators of
and where (with ) are the generators of . As usual, , and are the gauge couplings for the , and groups, respectively, as shown in the Table 3. The action in the lepton and Higgs sectors are defined by the following Lagrangians
The last two terms define the superpotential of the MSSM as where
The supersymmetric parameter is a complex numbers and the terms are elements of the complex 3 × 3 Yukawa coupling matrices in the family space. The color indices on the triplet (antitriplet) superfield contract trivially, and have been suppressed. The second terms of the Lagrangian defined by the Equation (4) is given by the following equation
The gauge superfields have the following explicit form
where is the covariant derivative and it is given by 
and is a spinorial index.
In principle, one could add to the Lagrangian defined by the Equation (5) other terms that, even if they break the baryon number and the lepton number conservation laws, are still allowed by the supersymmetry. However, no physical process with this property has been discovered so far. This phenomenological fact suggest imposing a symmetry that rules out such terms called -parity which is defined in terms of the following operators
where and are the baryon and lepton numbers, respectively, and is the spin for a given state. The -parity of the Lagrangian implies that the usual particles of the SM have while their superpartners have . The terms that break -parity are
Here, we have suppressed the indices and is the antisymmetric tensor. Some of the coupling constants in the Equation (13) should be set to zero in order to avoid a too fast proton decay and neutron-anti-neutron oscillation     . The choice of the -parity violation couplings with and 3, are constrained from various low energy observables such as: 1) charge-current universality, 2) universality, 3) atomic parity violation. The bounds on the product can be obtained from the charged -meson decay mixing , and the transition as observed in  .
The experimental evidence suggests that the supersymmetry is not an exact symmetry. Therefore, supersymmetry breaking terms should be added to the Lagrangian defined by the Equation (5). One possibility is by requiring that the divergences cancel at all orders of the perturbation theory. The most general soft supersymmetry breaking terms, which do not induce quadratic divergence, where described by Girardello and Grisaru  . They found that the allowed terms can be categorized as follows: a scalar field with mass terms
a fermion field gaugino with mass terms
and finally trilinear scalar interaction terms
The terms in this case are similar with the terms allowed in the superpotential of the model we are going to consider next.
Taken all this information into account, we can add the following soft supersymmetry breaking terms to the MSSM
where the scalar mass term is given by the following relation
The matrices and are hermitian and and are real. The gaugino mass term is written as
Here, and are complex. Finally, there is an interaction term , see the Equation (13), of the form
The matrices and matrices are complex.
The total Lagrangian of the MSSM is obtained by adding all Lagrangians above
see the Equations (4), (17). The MSSM contains 124 free parameters  and the symmetry breaking parameters are completely arbitrary  . The main goal in the SUSY phenomenology is to find some approximation about the way we can break SUSY in order to have a drastic reduction in the number of these parameters6. Many phenomenological analyses adopt the universality hypothesis at the scale GeV:
The assumptions that the MSSM is valid between the weak scale and GUT scale, and that the “boundary conditions”, defined by the Equation (22) hold, are often referred to as mSUGRA, or minimal supergravity model. The mSUGRA model is completely specified by the parameter set  
The new free parameter is defined in the following way
where is the vev of while is the vev of the Higgs in the doublet representation of group. Due the fact that and are both positive, it imples that .
In the context of the MSSM, it is possible to give mass to all charged fermions. With this superpotential we can explain the mass hierarchy in the charged fermion masses as showed in   . On the other hand, give mass to the gauge bosons: the charged ones and the neutral ( and get a massless foton but the neutrinos remain massless. Due to this fact, it is generated a spectrum that contains five physical Higgs bosons, two neutral scalar , one neutral pseudoscalar , and a pair of charged Higgs particles . At the level of tree level, we can write the following relations hold in the Higgs sector   :
Therefore, the light scalar has a mass smaller than the gauge boson at the tree level. This implies that one has to consider the one-loop corrections which lead to the following result 
In the MSSM there are four neutralinos ( with ) and two charginos ( with )   .
The mass matrix of neutrinos from this model was studied in  -  . The mass matrix has two zero eigenvalues. Thus, there are two neutrinos , which are massless at the tree level. More realistic neutrino masses require radiative corrections     . The neutrinos are Majorana particles, therefore the neutrinoless double beta decay must be observed. Neutrinoless double beta decay is a sensitive probe of physics beyond the SM since it violates the conservation of the lepton number  -  . The nucleon level process of decay, , can be obtained via the lepton number violating sub-process .
The supersymmetric mechanism of decay was first suggested by Mohapatra  and further studied in   . In  , the -parity violating Yukawa coupling of the first generation is strongly bounded by due to the gluino exchange -decay. Babu and Mohapatra  have latter implemented another contribution comparable with that via the gluino exchange. This set stringent bounds on the products of -parity violating Yukawa couplings of th generation index 
On the other hand, the confrontation of the experimental results with the predictions of the MSSM set a phenomenological constraint on the magnitude of . Indeed, the mass of Higgssino from the Equation (8) and the terms from the , given by the Equation (17), are of order of the electro- weak scale of while the natural cut-off scale is the Planck scale . The MSSM does not provide any mechanism to explain the difference between the two scales. This is know as the problem.
In order to address this problem, the Next-to-the-Minimal Supersymmetric Standard-Model (NMSSM)   was developed within the framework of the Grand Unification Theory (GUTs) as well as the superstring theorie    7. The NMSSM is characterized by the a new singlet field introduced in the following chiral superfield8 
where is the scalar in the singlet and its vacuum expectation value is given by . Its superpartner is known as the singlino. The rest of the particle content of this model is the same as of the MSSM given above in the Tables 1-3. The superpotential of the NMSSM model has the following form
where is defined by the Equation (9). The way in which the -problem is solved in the NMSSM is by generating dynamically the term in the superpotential through with a dimensionless coupling and the vacuum expectation value of the Higgs singlet. Another essential feature of the NMSSM is the fact that the mass bounds for the Higgs bosons and neutralinos are weakened. For more details about the scalar sector of this model see  . We summarize them in the Table 4 below.
Note that the neutralino sector is extended to a 5 × 5 mass matrix. If the following vector basis for fields is adopted (see, e.g.  )
Table 4. The physical Higgs states of the NMSSM  , the parameter is defined by the Equation (24).
We note that the singlino does not mix directly with the gauginos but it can mix with the neutral higgsinos  . The neutrinos are massless. However, the term can produce an effective term when the sneutrinos get vev as we will show later on. This would allow us to solve the problem  , without having to introduce an extra singlet superfield as we have done in the NMSSM. This new model is called “ from ” Supersymmetric Standard Model ( SSM). The field content of the SSM is the same as MSSM supplemented by three neutrino superfields  and is given by the following equation given by:
If the terms like are considered, then the term is induced when the right handed sneutrinos acquires a vev. By adding right handed neutrinos to the model, one can choose only the terms that break the lepton number conservation instead of the ones that break the baryon number conservation. The vev of these models are
In this case, all the neutrinos of the model can get mass at the tree level. Then the double beta decay can occur and the nucleon is stabilized.
The symmetry generates the following transformation of each chiral superfield
where is an entire number. The superpotential of this model can be obtained by requiring that it be -symmetric invariant. As a consequence, it takes the following form
where is defined by the Equation (9). It turns out that symmetry forbids all the bilinear terms in the superpotential. The expression obtained in the Equation (36) is consistent with the phenomenological models derived from the superstring theory that generate only trilinear couplings.
In the present context the string theory is relevant to the unification of all interactions, including gravity. The term proportional to gives an effective Majorana mass term to neutrinos, while the coupling generates Dirac mass term to neutrinos.
When the scalar components of the superfields , denoted by , acquire vev’s of the order of the electroweak scale, an effective interaction is generated with the effective coupling given by
In the same situation, the term can be generated with
the contribution of to the minimization conditions for the left- handed neutrinos . That provides an explanation for the neutrino’s masses in MSSM  .
In this model the -parity (and also the lepton number conservation) is broken explicitly. One of the candidates for the dark matter in NMSSM is the gravitino. Recently, some experimental bounds on gravitino masses were presented in see  . For an analysis of the gravitino as dark matter without -parity see      . Other possibilities that LSP be the axino were presented in  . The mass spectrum of this model can be found in  and ths spectrum is consistent with the experimental values obtained for both masses and mixing. The nice phenomenological aspects of this model were discussed in  . There are some works in SSM that consider gravitino as dark matter    .
3. A “New” SSM
In this section we propose a generalization of the SSM by adding new interaction terms that explicitly break the -parity and the lepton number symmetries, respectively. Therefore, the new model has potentially interesting cosmological consequences such as flat directions that provide a mechanism for the cosmological inflation and leptogenesis which explains the asymmetry between the matter and the anti-matter. We determine the flat directions of the generalized SSM and explain the leptogenesis mechanism. Then we show how our proposal addresses the recent experimental results from CMS and LHCb obtained in   and discussed in  .
3.1. New Terms in the Superpotential of SSM Model
The superpotential of the SSM model given by the Equation (36) can be generalized as follows
In the above equation, we have introduced two new terms that explicitly break the -parity and the lepton number symmetry. There is a new parameter that generates one more contribution to the mixing between the usual leptons with higgsinos. The usual techniques allow to determine from the superpotential the following mass matrix elements
The terms that describe the mixing between the usual leptons with the gauginos are the same as in the MSSM. In our notation, they are given by the following relations
The mixing between the usual leptons with the higgsinos is given by the equation
Beside the new mixing sectors given above, there are interactions between gauginos and higgsinos given by the same terms as in the MSSM. One can calculate the mass matrices of the charged leptons following the reference  . The result in the basis is given by the matrix
Here, the winos (superpartners of the -boson) defined as
See also the Equation (34). In the neutralino sector we use the basis
Then the mass matrices take the following form
where is the neutralino mass matrix presented in  . However, one should emphasize that the neutrinos are Majorana particles. Therefore, the double beta decay without neutrinos can occur in this model. This process is permitted by the new first term from the Equation (39). The second term alone will not induce neither the fast proton decay, nor the neutron anti-neutron oscillation  .
It is important to note that the last term in the superpotential will modify the down quarks masses. Some algebra shows that they are given by the following relation
The new superpotential given by the Equation (39) has all properties required in    . It has the advantage that it induces the double beta decay while maintaining the nucleon stability  as was discussed in the previous section.
3.2. Flat Direction of SSM Model
One of the most remarkable aspects of the supersymmetric gauge theories is that they have a vacuum degeneracy at the classical level. It is a well established fact that the renormalizable scalar potential is a sum of squares of -terms and - terms which implies that it can vanish identically along certain flat directions10 in the space of fields11. The properties of the space of flat directions of a supersymmetric model are crucial for making realistic considerations in cosmology and whenever the behaviour of the theory at large field strengths is an issue.
In the MSSM the flat directions can give rise to a host of cosmologically interesting dynamics (for a review, see  ). These include Affleck-Dine baryogenesis    , the cosmological formation and fragmentation of the MSSM flat direction condensate and subsequent -ball formation   -  , reheating the Universe with -ball evaporation   , generation of baryon isocurvature density perturbations    , as well as curvaton scenarios where MSSM flat directions reheat the Universe and generate adiabatic density perturbations  -  . Adiabatic density pertur- bations induced by fluctuating inflaton-MSSM flat direction coupling has also been discussed in   .
We can calculate all flat directions in the MSSM using the prescription given in   . It turns out that a MSSM flat direction is some linear combination of the MSSM scalars and can be thought of as a trajectory in the moduli space described by a single scalar degree of freedom.
Using the same technique from   we can calculate the flat directions in SSM model. The Table 5 presents the computed flat directions when the only terms taken into account are the renormalizable terms. In this model gives a particular flat direction. It follows that the field
generates a similar inflanton scenario as the one discussed in    . Another flat direction is given by the interaction term from which is generated the following field
It will be interesting to compare the cosmological consequences of the fields and given above.
3.3. Leptogenesis in SSM Model
The mechanisms to create a baryon asymmetry from an initially symmetric state must in general satisfy the three basic conditions for baryogenesis as pointed out by Sakharov  :
1) Violate baryon number, , conservation.
2) Violate and conservation.
3) To be out of thermal equilibrium.
It is found that the violation observed in the quark sector  , e.g. in or mesons system, is far too small to give rise to the observed baryon asymmetry  . Therefore, these conditions should be extended to include the lepton number violation processes.
The classic leptogenesis scenario of Fukugita and Yanagida  described in   can occur in the SSM model  . The Yukawa coupling can then induce heavy right handed neutrino decays via the following two channels:
Table 5. Flat direction of the model SSM.
that violate the lepton number by one unit. In the new SSM model given by the superpotential (39), it is one of the heavies neutralinos that is responsible for the right handed neutrion decay according to the Equation (46). All Sakharov’s conditions for leptogenesis are satisfied if these decays violate and go out of equilibrium at some stage during the evolution of the early universe. The requirement for violation means that the coupling matrix must be complex and the mass of must be greater than the combined mass of and , so that the interferences between the tree-level processes and the one-loop corrections with on-shell intermediate states will be non-zero   . Since is the scalar field of the SM, the usual Higgs can suffer the following decays
Note that none of these decays violate the lepton number conservation. Nevertheless, in this model the fields have both chiralities. Therefore, they will induce the followings decays
Thus, both violate the lepton number conservation. On the other hand, we note that there are scattering processes that can alter the abundance of the neutrino flavour in the -channel and -channel
besides the tree-level interaction . In addition to these, there are also scattering processes mediated by which can be important for the evolution of . Also, if we consider the couplings and to be complex, we can generate the leptogenesis in this model as shown in  by inducing decays as .
It is interesting to note that the superpotential from the Equation (39) induces the following processes    
1) New contributions to the neutrals and Systems.
2) New contributions to the muon decay.
3) Leptonic Decays of Heavy Quarks Hadrons such as .
4) Rare Leptonic Decays of Mesons like .
5) Hadronic Meson Decay Asymmetries.
Also, ir gives the following direct decays of the lightest neutralinos
and for the lightest charginos
These decays are similar to the ones from the MSSM when -Parity violating scenarios are taken into account. Therefore, we expect that the missing energy plus jets be the main experimental signal in the “new” SSM as is in the MSSM. These decays violate only the lepton number conservation but they conserve the baryon number.
As we have seen above, all necessary conditions to generate a viable leptogenesis mechanism from the SSM model are present in the “new” SSM model  as well as the CP violation processes. Also, this model could contain an invisible axion. These properties deserve a deeper study. Another interesting phenomenological avenue is to analyse the total cross section of the Dark Matter-Nucleon (DM-N) elastic scattering process.
4. Explanation of the Data from ATLAS, CMS and LHCb in SSM Model
One possible explanation to the excess of electrons is given if the following processes are considered  
Neglecting finite width effects, the color and spin-averaged parton total cross section of a single slepton production is  
where is the partonic center of mass energy, and is the mass of the resonant slepton. Including the effects of the parton distribution functions, we find the total cross section
to a good approximation in the parameter region of interest.
As was discussed in  , these processes represent one of the possibilities to explain the data of CMS   if the selectron mass is fixed to and the lightest neutralino mass is taken to be in the range from 400 GeV up to 1 TeV. The measurement can be consistent with the new physics arising from the electron or muon sector of the SM and it was shown in  that if we consider the muon sector in the MSSM with -parity violation scenarios, the can also account for both data arising from CMS and LHCb. In the “new” SSM model we have both terms present. With respect with the di-boson data, there is a similar explanation. Indeed, in the case of there is the single production of smuons  , while in the case of di-photons the stau is produced  . Due this fact, we expect that our model fit the new data coming from ATLAS  , CMS  and from LHCb   . To confirm that this is the true mechanism employed, the double beta decay must be detected in experiments like CUORE  , GERDA  and MAJORANA  and no proton decay must occur in the neutron anti-neutron oscillation.
In this article we have reviewed some of the basic properties of the MSSM, NMSSM and SSM essential to the cosmological applications. Also, in order to incorporate the recent data from the CMS and LHCb into this class of models, we have proposed a “new” SSM model characterized by the superpotential given in the Equation (39). The terms added to the of the SSM in order to obtain the modified model, explicitly break the -parity and the lepton number conservation. This makes the model attractive for cosmological applications as it presents flat directions that represent a possibility to generate inflation and a viable leptogenesis mechanism that is necessary to generate the matter anti-matter asymmetry. These properties make the model interesting for further investigations on which we hope to report in the near future.
M. C. R would like to thanks to Laboratório de Fsica Experimental at Centro Brasileiro de Pesquisas F sicas (LAFEX-CBPF) for their nice hospitality and special thanks to Professores J. A. Helayël-Neto and A. J. Accioly. Both authors acknowledge R. Rosenfeld for hospitality at ICTP-SAIFR where part of this work was accomplished. We acknowledge P. S. Bhupal Dev for information on the latest results from ATLAS and CMS in di-bosons, and D. E. Lopez-Fogliani for useful correspondence on the gravitino in the SSM.
1The excess was at a di-boson invariant masses in the range from 1.3 to 3.0 TeV.
2The ressonance appear at around 750 GeV in the di-photon invariant mass.
3About the history of MSSM, see e.g.   .
4The term generate the term when the sneutrino get its values expectation values.
5We recall that the flat directions provide a viable mechanism to generate the cosmological inflation and the leptogenesis is important to explaining the asymmetry between the matter and the anti-matter.
6Different assumptions result in different version of the Constrained Minimal Supersymmetric Model (CMSSM).
7References to the original work on the NMSSM may be found in the reviews   .
8 , where are the three Pauli matrices plus the the identity matrix.
9Where we have defined .
10The flat directions are noncompact lines and surfaces in the space of scalars fields along which the scalar potential vanishes. The present flat direction is an accidental feature of the classical potential and gets removed by quantum corrections.
11In quantum field theories, the possible vacua are usually labelled by the vacuum expectation values of scalar fields, as Lorentz invariance forces the vacuum expectation values of any higher spin fields to vanish. These vacuum expectation values can take any value for which the potential function is a minimum. Consequently, when the potential function has continuous families of global minima, the space of vacua for the quantum field theory is a manifold (or orbifold), usually called the vacuum manifold. This manifold is often called the moduli space of vacua, or just the moduli space.
 Klapdor-Kleingrothaus, H.V., et al. (2001) Latest Results from the HEIDELBERG-MOSCOW Double Beta Decay Experiment. The European Physical Journal A, 12, 147-154.
 Aalseth, C.E., et al. (1999) Neutrinoless Double-β Decay of 76 Ge: First Results from the International Germanium Experiment (IGEX) with Six Isotopically Enriched Detectors. Physical Review C, 59, 2108.
 Avignone, F.T. (2008) The Majorana 76Ge Neutrino Less Double-Beta Decay Project: A Brief Update. Journal of Physics: Conference Series, 120, Article ID: 052059.
 Spergel, D.N., et al. (2007) Three-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Beam Profiles, Data Processing, Radiometer Characterization, and Systematic Error Limits. The Astrophysical Journal Supplement Series, 170, 377.
 Seljak, U., et al. (2005) Cosmological Parameter Analysis Including SDSS Ly α Forest and Galaxy Bias: Constraints on the Primordial Spectrum of Fluctuations, Neutrino mass, and Dark Energy. Physical Review D, 71, Article ID: 103515.
 Aad, G., et al. (2015) Search for High-Mass Diboson Resonances with Boson-Tagged Jets in Proton-Proton Collisions at √ s = 8 TeV with the ATLAS Detector. Journal of High Energy Physics, 2015: 55. [arXiv:1506.00962 [hep-ex]]
 CMS Collaboration (2016) Search for New Physics in High Mass Diphoton Events in 3.3 fb-13.3 fb-1 of Proton-Proton Collisions at s√=13 TeVs and Combined Interpretation of Searches at 8 TeV and 13 TeV.
 Amaldi, U., de Boer, W. and Fürstenau, H. (1991) Comparison of Grand Unified Theories with Electroweak and Strong Coupling Constants Measured at LEP. Physics Letters B, 260, 447.
 Allanach, B.C., Kom, C.H. and Pas, H. (2009) Large Hadron Collider Probe of Supersymmetric Neutrinoless Double-Beta-Decay Mechanism. Physical Review Letters, 103, Article ID: 091801.
 Haber, H.E. and Hempfling, R. (1991) Can the Mass of the Lightest Higgs Boson of the Minimal Supersymmetric Model be Larger Than mZ. Physical Review Letters, 66, 1815.
 Borzumati, F., Hamaguchi, K. and Yanagida, T. (2001) Supersymmetric See Saw Model for the (1 + 3)-Scheme of Neutrino Masses. Physics Letters B, 497, 259-264.
 Baudis, L., et al. (1997) The Heidelberg-Moscow Experiment: Improved Sensitivity for 76Ge Neutrinoless Double Beta Decay. Physics Letters B, 407, 219-224.
 Klapdor-Kleingrothaus, H.V., Dietz, A., Harney, H.L. and Krivosheina, I.V. (2001) Evidence for Neutrinoless Double Beta Decay. Modern Physics Letters A, 16, 2409.
 Chemtob, M. (2005) Phenomenological Constraints on Broken R Parity Symmetry in Supersymmetry Models. Progress in Particle and Nuclear Physics, 54, 71-191.
 Hirsch, M., Klapdor-Kleingrothaus, H.V. and Kovalenko, S.G. (1995) New Contributions to Supersymmetric Mechanism of Neutrinoless Double Beta Decay. Physics Letters B, 352, 1-7.
 Hirsch, M., Klapdor-Kleingrothaus, H.V. and Kovalenko, S.G. (1995) New Constraints on R-Parity-Broken Supersymmetry from Neutrinoless Double Beta Decay. Physical Review Letters, 75, 17.
 Babu, K.S. and Mohapatra, R.N. (1995) New Vector-Scalar Contributions to Neutrinoless Double Beta Decay and Constraints on R-Parity Violation. Physical Review Letters, 75, 2276.
 Ellwanger, U., Rausch de Traubenberg, M. and Savoy, C.A. (1993) Particle Spectrum in Supersymmetric Models with a Gauge Singlet. Physics Letters B, 315, 331-337.
 Ananthanarayan, B. and Pandita, P.N. (1997) The Nonminimal Supersymmetric Standard Model at Large Tan β. International Journal of Modern Physics A, 12, 2321.
 Choi, K.Y., Lopez-Fogliani, D.E., Munoz, C. and de Austri, R.R. (2010) Gamma-Ray Detection from Gravitino Dark Matter Decay in the Mu Nu SSM. Journal of Cosmology and Astroparticle Physics, 1003, 028.
 Gomez-Vargas, G.A., Fornasa, M., Zandanel, F., Cuesta, A.J., Munoz, C., Prada, F. and Yepes, G. (2012) CLUES on Fermi-LAT Prospects for the Extragalactic Detection of Mu Nu SSM Gravitino Dark Matter. Journal of Cosmology and Astroparticle Physics, 02, 001.
 Albert, A., et al. (2014) Search for 100 MeV to 10 GeV Gamma-Ray Lines in the Fermi-LAT Data and Implications for Gravitino Dark Matter in the μνSSM. Journal of Cosmology and Astroparticle Physics, 10, 023.
 Enqvist, K., Kasuya, S. and Mazumdar, A. (2003) Adiabatic Density Perturbations and Matter Generation from the Minimal Supersymmetric Standard Model. Physical Review Letters, 90, Article ID: 091302.
 Postma, M. and Mazumdar, A. Resonant Decay of Flat Directions: Applications to Curvaton Scenarios, Affleck-Dine Baryogenesis, and Leptogenesis from a Sneutrino Condensate. arXiv: hep-ph/0304246.
 Gherghetta, T., Kolda, C.F. and Martin, S.P. (1996) Flat Directions in the Scalar Potential of the Supersymmetric Standard Model. Nuclear Physics B, 468, 37-58.
 Kajiyama, Y., Khalil, S. and Raidal, M. (2009) Electron EDM and Soft Leptogenesis in Supersymmetric B-L Extension of the Standard Model. Nuclear Physics B, 820, 75.
 Allanach, B.C., Dev, P.S.B., Renner, S.A. and Sakurai, K. (2015) Di-Photon Excess Explained by a Resonant Sneutrino in R-Parity Violating Supersymmetry. Physical Review D, 93, Article ID: 115022.