Recent detection of the neutrino events IceCube-170922A, 13 muon-neutrino events observed in 2014-2015 and IceCube-141209A by IceCube observatory from the Blazars, namely TXS 0506+056, PKS 0502+049/TXS 0506+056 and GB6 J1040+0617 respectively in the state of enhanced gamma-ray emission, indicates the acceleration of cosmic rays in the blazar jets. The photo-meson ($p\gamma$) interaction cannot explain the IceCube observations of 13 neutrino events. The non-detection of broadline emission in the optical spectra of the IceCube blazars, however, question the hadronuclear (pp) interaction interpretation through relativistic jet meets with high density cloud. In this work, we investigate the proton blazar model in which the non-relativistic protons that come into existence under the charge neutrality condition of the blazar jet can offer sufficient target matter for $pp$ interaction with shock-accelerated protons, to describe the observed high-energy gamma-rays and neutrino signal from the said blazars. Our findings suggest that the model can explain consistently the observed electromagnetic spectrum in combination with appropriate number of neutrino events from the corresponding blazars.

This has been suggested that the dwarf satellite galaxies in the Milky Way may contain substantial amount of dark matter in them. These dark matters may undergo self-annihilation to produce $\gamma$-rays. The satellite borne $\gamma$-rays telescope such as Fermi-LAT reported the detection of $\gamma$-rays from around 45 Dwarf Spheroidals (dSphs) of Milky Way. In this work, we consider a particle dark matter model and after studying its phenomenology, we calculate the $\gamma$-ray fluxes of each of these 45 dSphs and compare our results with the upper bounds of mass vs annihilation cross-sections of dark matter provided by the Fermi-LAT collaboration. We calculate these fluxes by considering different dark matter density profiles and make a comparison of the results. We also repeat our analysis with another dark matter candidate namely Kaluza-Klein dark matter inspired by extra-dimensional models. We make a critical comparison between the results obtained for each of these models vis-a-vis the observational results.

Considering the ultrahigh energy (UHE) neutrino events reported by IceCube in the PeV regime to have originated from the decay of superheavy dark matter, the IceCube UHE neutrino events are analysed and the best fit values of the two parameters namely the mass of the superheavy dark matter and its decay lifetime are obtained. The theoretical astrophysical flux is also included in theanalysis. We find that while the neutrino events in the energy range $\sim$ 60 TeV-$\sim$ 120 TeV appears to have astrophysical origin, the events in the energy range $\sim 1.2 \times 10^5$ GeV - $\sim 5 \times 10^7$ GeV can be well described from the superheavy dark matter decay hypothesis. We also find that although hadronic decay channel of the superheavy dark matter can well explain the events in the energy range $\sim 1.2 \times 10^5$ GeV - $\sim 5 \times 10^6$ GeV, the higher energy regime higher than this range can be addressed only when the leptonic decay channel is considered.

We consider a dark energy scenario driven by a scalar field $\Phi$ with a pseudo Nambu Goldstone boson (pNGB) type potential $V(\phi)=\mu^4 \left( 1+ {\rm cos}(\phi/f) \right)$. The pNGB originates out of breaking of spontaneous symmetry at a scale $f$ close to Planck mass $M_{\rm{pl}}$. We consider two cases namely the quintessence dark energy and the other, where the standard pNGB action is modified by the terms related to Slotheon cosmology. We demonstrate that for this pNGB potential, high-$f$ problem is better addressed when interaction between dark matter and dark energy is taken into account and that Slotheon dark energy scenario works even better over quintessence in this respect. To this end, a mass limit for dark matter is also estimated.

In this thorough study we focus on the indirect detection of Dark Matter (DM) through the confrontation of unexplained galactic and extragalactic $\gamma$-ray signatures for a low mass DM model. For this, we consider a simple Higgs-portal DM model, namely, the inert Higgs doublet model (IHDM) where the Standard Model is extended with an additional complex SU(2)$_L$ doublet scalar. The stability of the DM candidate in this model, i.e., the lightest neutral scalar component of the extra doublet, is ensured by imposing discrete $\mathbb{Z}_2$ symmetry. The reduced-$\chi^2$ analysis with the theoretical, experimental and observational constraints suggests the best-fit value of DM mass in this model to be $\sim$ 63.5 GeV. We analyse the anomalous GeV $\gamma$-ray excess from Galactic Centre in light of the best-fit IHDM parameters. We further check the consistency of the best-fit IHDM parameters with the Fermi-LAT obtained limits on photon flux for 18 Milky Way dwarf spheroidal satellite galaxies (dSphs) known to be mostly dominated by DM. Also since the $\gamma$-ray signal from DM annihilation is assumed to be embedded within the extragalactic $\gamma$-ray background (EGB), the theoretical calculations of photon flux for the best-fit parameter point in the IHDM framework are compared with the Fermi-LAT results for diffuse and isotropic EGB for different extragalactic and astrophysical background parametrisations. We show that the low mass DM in IHDM framework can satisfactorily confront all the observed continuum $\gamma$-ray fluxes originated from galactic as well as extragalactic sources. The extensive analysis performed in this work is valid for any Higgs-portal model with DM mass in the ballpark of that considered in this work.

We promote the idea of multi-component Dark Matter (DM) to explain results from both direct and indirect detection experiments. In these models as contribution of each DM candidate to relic abundance is summed up to meet WMAP/Planck measurements of $\Omega_{\rm DM}$, these candidates have larger annihilation cross-sections compared to the single-component DM models. This results in larger $\gamma$-ray flux in indirect detection experiments of DM. We illustrate this fact by introducing an extra scalar to the popular single real scalar DM model. We also present detailed calculations for the vacuum stability bounds, perturbative unitarity and triviality constraints on this model. As direct detection experimental results still show some conflict, we kept our options open, discussing different scenarios with different DM mass zones. In the framework of our model we make an interesting observation: The existing direct detection experiments like CDMS II, CoGeNT, CRESST II, XENON 100 or LUX together with the observation of excess low energy $\gamma$-ray from Galactic Centre and Fermi Bubble by FGST already have the capability to distinguish between different DM halo profiles.

We propose a two parameter generalization for the dark energy equation of state (EOS) $w_X$ for thawing dark energy models which includes PNGB, CPL and Algebraic thawing models as limiting cases and confront our model with the latest observational data namely SNe Ia, OHD, CMB, BOSS data. Our analysis reveals that the phantom type of thawing dark energy is favoured upto $2 \sigma$ confidence level. These results also show that thawing dark energy EOS is not unique from observational point of view. Though different thawing dark energy models are not distinguishable from each other from best-fit values (upto $2\sigma$ C.L.s) of matter density parameter ($\Omega_m^0$) and hubble parameter ($H_0$) at present epoch, best-fit plots of linear growth of matter perturbation ($f$) and average deceleration parameter ($q_{\rm av}$); the difference indeed reflects in best-fit variations of thawing dark energy EOS, model-independent geometrical diagnostics like the statefinder pair $\{r,s\}$ and $Om3$ parameter. We are thus led to the conclusion that unlike the standard observables ($\Omega_m^0$, $H_0$, $f$, $q_{\rm av}$), the model-independent parameters ($r,s,Om3$) and the variations of EOS (in terms of $w_X-w_X'$ plots) serve as model discriminators for thawing dark energy models.

We consider the lightest supersymmetric particle (LSP), neutralino in minimal anomaly mediated supersymmetry breaking model (mAMSB) to be a possible candidate for weakly interacting massive particles (WIMP) or cold dark matter and investigate its direct and indirect detections. The theoretically allowed supersymmetric parametric space for such a model along with the recent bounds from LHC is constrained by the WMAP results for relic densities. The spin independent and spin dependent scattering cross sections for dark matter off nucleon are thus constrained from the WMAP results. They are found to be within the allowed regions of different ongoing direct detection experiments. The annihilation of such dark matter candidates at the galactic centre produce different standard model particles such as gamma rays, neutrinos etc. In this work, we calculate the possible fluxes of these $\gamma$-rays and neutrinos coming from the direction of the galactic centre (and its neighbourhood) at terrestrial or satellite borne detectors. The calcutated $\gamma$-ray flux is compared with the observational results of HESS experiment. The neutrino flux of different flavours from the galactic centre and at different locations away from the galactic centre produced by WIMP annihilation in this model are also obtained for four types of galactic dark matter halo profiles. The detection prospects of such $\nu_\mu$ coming from the direction of the galactic centre at the ANTARES under sea detector are discussed in terms of muon signal yield from these muon neutrinos. Both the gamma and neutrino signals are estimated for four different dark matter halo profiles.

In this paper, we examine the cosmological viability of a slow moving galileon field in a potential. The Lagrangian $\cal{L}$ $ =-\frac{1}{2}g^{\mu\nu}\pi_{;\mu}\pi_{;\nu} +\frac{G^{\mu\nu}}{2M^2}\pi_{;\mu}\pi_{;\nu}$ respects the galileon symmetry in curved space time. We carry out detailed investigations of the underlying dynamics of this Lagrangian with Einstein-Hilbert term and a potential. We demonstrate that the model can give rise to a viable ghost free late time acceleration of universe. Furthermore we study the cosmological perturbation of the model and see that the model gives different BBN constraints at early times. We also carry out the observational analysis of the model and use observational data from growth, Type Ia Supernovae (SNIa), Baryon Acoustic Oscillations (BAO) and Cosmic Microwave Background (CMB) to constrain the parameters of the theory.

From the analysis of Supernova Ia data alongwith Observational Hubble Data (OHD) and Baryon Acoustic Oscillation (BAO) data, we attempt to find out the nature of a scalar potential that may be responsible for the Dark Energy of the universe. We demonstrate that in order to explain the varying dark energy equation of state ($\omega_X(z)$) as obtained in a model independent way from the analyses of observational data, we need to invoke a quintom scalar field having a "quintessence" part for $\omega_X(z) > -1$ and a "phantom" part for $\omega_X(z) < -1$. We consider a Gaussian type potential for these scalar fields and compare the dark energy equation of state derived from such potential with the one computed from the data analysis.

We consider a simplest extension of Standard Model in which a real SM gauge singlet scalar with an additional discrete symmetry $Z_2$ is introduced to SM. This additional scalar can be a viable candidate of cold dark matter since the stability of $S$ is achieved by the application of $Z_2$ symmetry on $S$. Considering $S$ as a possible candidate of cold dark matter we have solved Boltzmann's equation to find the freeze out temperature and relic density of $S$ for Higgs mass 120 GeV in the scalar mass range 5 GeV to 1 TeV. As $HHSS$ coupling $\delta_2$ appearing in Lagrangian depends upon the value of scalar mass $m_S$ and Higgs mass $m_h$, we have constrained the $m_S - \delta_2$ parameter space by using the WMAP limit on the relic density of dark matter in the universe and the results of recent ongoing dark matter direct search experiments namely CDMS-II, CoGeNT, DAMA, EDELWEISS-II, XENON-10, XENON-100. From such analysis we find two distinct mass regions (a lower and higher mass domain) for such a dark matter candidate that satisfy both the WMAP limit and the experimental results considered here. We have estimated the possible differential direct detection rates and annual variation of total detection rates for this scalar dark matter candidate $S$ for two detector materials namely Ge, Xe. Finally we have calculated the $\gamma-$ray flux from the galactic centre due to annihilation of two 130 GeV scalar dark matter into two monoenergetic $\gamma-$rays.

Direction sensitive direct detection of Weakly Interacting Massive Particles (WIMPs) as dark matter would provide an unambiguous non-gravitational signature of dark matter (DM). The diurnal variation of DM signal due to earth's rotation around its own axis can be a significant signature for galactic WIMPs. Because of particular orientation of earth's axis of rotation with respect to WIMP wind direction, the apparent direction of WIMP wind as observed at a detector can alter widely over a day. In this work we calculate the directional detection rates with their daily and yearly modulations in earth-bound dark matter experiments considering detailed features of the geometry and dynamics of the earth-sun system along with the solar motion in galactic frame. A separate halo model namely the dark disc model other than the usual standard halo model for dark matter halo is also considered and the results for two models are compared. We demonstrate the results for two types of gas detectors namely DRIFT (target material CS2) and NEWAGE (target material CF4) that use Time Projection Chamber techniques for measuring directionality of the recoil nucleus. The WIMP mass and recoil energy dependence of the daily variation of event rates are computed for specific detector and the sensitive ranges of mass and recoil energies for the considered detector are probed.

We adopt a model independent method to reconstruct the dark energy equation of state by analyzing 5 sets of SNe Ia data along with Baryon Acoustic Oscillation (BAO) and Observational Hubble Data (OHD). The SNe Ia data sets include the most recent UNION2 data and other data compilations from the year 2007 to the present. We assume a closed form parametrization of the luminosity distance in terms of redshift and perform a $\chi^2$ analysis of the observational data. The matter density at the present epoch $\Omega_m^0$ is also taken to be a parameter in the analysis and its best-fit values are obtained for each of the data sets. We found a strong dependence of dark energy equation of state on the matter density in the present and earlier epoch. From the analysis, we also predict the lower limit of matter density parameter at an earlier epoch within 1$\sigma$ confidence level for a flat FRW universe. The dark energy equation of state appears to be a slow varying function of $z$. The variation of dark energy density parameter and the matter density parameter are also shown along with their 1$\sigma$ variations.

We consider the recent limits on dark matter - nucleon elastic scattering cross section from the analysis of CDMS II collaboration using the two signal events observed in CDMS experiment. With these limits we try to interpret the Super-Kamiokande (SK) bounds on the detection rates of up-going muons induced by the neutrinos that are produced in the sun from the decay of annihilation products of dark matter (WIMPs) captured in the solar core. Calculated rates of up-going muons for different annihilation channels at SK using CDMS bounds are found to be orders below the predicted upper limits of such up-going muon rates at SK. Thus there exists room for enhancement (boost) of the calculated rates using CDMS limits for interpreting SK bounds. Such a feature is expected to represent the PAMELA data with the current CDMS limits. We also show the dependence of such a possible enhancement factor (boost) on WIMP mass for different WIMP annihilation channels.

We investigate deviation from the tri-bimaximal mixing in the case of ultra high energy neutrino using ICECUBE detector. We consider the ratio of No. of muon tracks to the shower generated due to eletrons and hadrons. Our analysis shows for tri-bimaximal mixing the ratio comes out around 4.05. Keeping $\theta_12$ and $\theta_23$ fixed at tri-bimaximal value, we have varied the angle $\theta_{13}$ = $3^o$, $6^o$, $9^o$ and the value of the ratio gradually decreases. The variation of ratio lies 8% to 18% from the tri-bimaximal mixing value.

The half lives are calculated for the process of $\beta^\pm$ decay and electron capture for nuclei in mass range $\sim$ 65 - 100 relevant for the core of a massive star at late burning stage of steller evolution that leads to supernova explosion. The rates of electron capture process are also computed for nuclei relevant during presupernova stage of massive stars. These half lives and rates are calculated by expressing $\beta^\pm$ Gamow-Teller decay strengths in terms of smoothed bivariate strength densities. These strength densities are constructed in the framework of spectral averaging theory for two body nuclear Hamiltonian in a large nuclear shell model space. The method has a natural extension to weak interaction rates for r and rp-processes.

We re-examine MSW solutions of the solar neutrino problem in a two flavor scenario taking (a) the results on total rates and the electron energy spectrum from the 1117-day SuperKamiokande (SK) data and (b) those on total rates from the Chlorine and Gallium experiments. We find that the SMA solution gives the best fit to the total rates data from the different experiments. One new feature of our analysis is the use of the moments of the SK electron spectrum in a $\chi^2$ analysis. The best-fit to the moments is broadly in agreement with that obtained from a direct fit to the spectrum data and prefers a $\Delta m^2$ comparable to the SMA fit to the rates but the required mixing angle is larger. In the combined rate and spectrum analysis, apart from varying the normalization of the $^8$B flux as a free parameter and determining its best-fit value we also obtain the best-fit parameters when correlations between the rates and the spectrum data are included and the normalization of the $^8$B flux held fixed at its SSM value. We observe that the correlations between the rates and spectrum data are important and the goodness of fit worsens when these are included. In either case, the best-fit lies in the LMA region.

The SuperKamiokande collaboration has presented results on the observation of solar neutrinos. The Sudbury Neutrino Observatory (SNO) is also expected to go on-line in the near future. We propose several new variables, insensitive to the absolute flux of the initial solar neutrino beam, which probe the shape of the observed spectrum at these experiments and can sensitively signal neutrino oscillations. One class of such variables involves normalised moments of the distributions recorded at the two facilities while another variable, specific to SNO, depends on the integrated charged and neutral current signals. The utility of these variables in the context of supernova neutrinos, both from the collapse epoch and the post-bounce era, is also emphasised. It is shown that, notwithstanding the imprecise nature of the information about the initial neutrino spectra from a supernova, oscillations can be detected using these variables and it will be possible to distinguish between the alternatives of oscillation to a sequential neutrino {\em vis-a-vis} that to a sterile neutrino.

Results for solar neutrino detection from the SuperKamiokande collaboration have been presented recently while those from the Sudbury Neutrino Observatory are expected in the near future. These experiments are sensitive to the 8B neutrinos from the sun, the shape of whose spectrum is well-known but the normalisation is less certain. We propose several variables, insensitive to the absolute flux of the incident beam, which probe the shape of the observed spectrum and can sensitively signal neutrino oscillations. They provide methods to extract the neutrino mixing angle and mass splitting from the data and also to distinguish oscillation to sequential neutrinos from those to a sterile neutrino.

Neutrinos emitted during stellar core collapse leading to a supernova are primarily of the electron neutrino type at source which may undergo oscillation between flavor eigenstates during propagation to an earth-bound detector. Although the number of neutrinos emitted during the pre-bounce collapse phase is much smaller than that emitted in the post-bounce phase (in which all flavors of neutrinos are emitted), a nearby supernova event may nevertheless register a substantial number of detections from the pre-bounce phase at SuperKamiokande (SK) and the Sudbury Neutrino Observatory (SNO). The calorimetric measurement of the supernova neutrino fluence from this stage via the charge current and neutral current detection channels in SNO and the corresponding distortion of detected spectrum in SK over the no-oscillation spectrum, can probe information about neutrino mass difference and mixing which are illustrated here in terms of two- and three-flavor oscillation models.