We explore the thermal light sterile neutrino situation from cosmological perspective in the $\Lambda \textrm{CDM} + r_{0.05} + N_{\textrm{eff}} + m^{\textrm{eff}}_{\textrm{s}}$ model using combinations of latest data sets available. Here, $r_{0.05}$ is the tensor-to-scalar ratio at the pivot scale of $k_*=0.05h$ Mpc$^{-1}$, $N_{\textrm{eff}}$ is the effective number of relativistic species during recombination, and $m^{\textrm{eff}}_{\textrm{s}}$ is the effective mass of the sterile neutrino. Among Cosmic Microwave Background (CMB) datasets, we use Planck 2015 temperature and low-$l$ ($l <$ 30) polarization data and the latest data release on the B-mode polarization up to and including 2014 from the BICEP2/Keck collaboration (BK14). We also use the latest BAO data from SDSS-III BOSS DR12, MGS, and 6dFS; and a Gaussian prior (HST) on the Hubble constant ($H_0 = 73.24 \pm 1.74$ km/sec/Mpc) from direct measurements. We find that inclusion of BK14 data makes the constraints on the effective mass of sterile neutrino ($m^{\textrm{eff}}_{\textrm{s}}$) slightly stronger by preferring higher $\sigma_8$ values. The bound of $m^{\textrm{eff}}_{\textrm{s}} <$ 0.46 eV (95\% C.L.) is found for the combination of Planck 2015, BAO and BK14 datasets, whereas the bound is $m^{\textrm{eff}}_{\textrm{s}} <$ 0.53 eV (95\% C.L.) without the BK14 data. Our most aggressive bound of $m^{\textrm{eff}}_{\textrm{s}} <$ 0.28 eV (95\% C.L.) is obtained with Planck 2015, HST and BK14. However, the HST prior also leads to very high $N_{\textrm{eff}}$ which might be in conflict with bounds from BBN. Our analysis indicates that fully thermalized sterile neutrinos with mass $\sim 1$ eV are slightly more disfavoured with the inclusion of BK14 data. It also seems to make the agreement between Planck 2015 and CFHTLenS (weak gravitational lensing data) worse due to the higher $\sigma_8$ values (abstract abridged).

The non-thermal dark matter (DM) production via the so-called freeze-in mechanism provides a simple alternative to the standard thermal WIMP scenario. In this work, we consider a popular $U(1)_{B-L}$ extension of the standard model (SM) in the context of inverse seesaw mechanism which has at least one (fermionic) FIMP DM candidate. Due to the added $\mathbb{Z}_{2}$ symmetry, a SM gauge singlet fermion, with mass of order keV, is stable and can be a warm DM candidate. Also, the same $\mathbb{Z}_{2}$ symmetry helps the lightest right-handed neutrino, with mass of order GeV, to be a stable or long-lived particle by making a corresponding Yukawa coupling very small. This provides a possibility of a two component DM scenario as well. Firstly, in the absence of a GeV DM component (i.e., without tuning its corresponding Yukawa coupling to be very small), we consider only a keV DM as a single component DM, which is produced by the freeze-in mechanism via the decay of the extra $Z'$ gauge boson associated to $U(1)_{B-L}$ and can consistently explain the DM relic density measurements. In contrast with most of the existing literature, we have found a reasonable DM production from the annihilation processes. After numerically studying the DM production, we show the dependence of the DM relic density as a function of its relevant free parameters. We use these results to obtain the parameter space regions that are compatible with the DM relic density bound. Secondly, we study a two component DM scenario and emphasize that the current DM relic density bound can be satisfied for a wide range of parameter space.

We present strong bounds on the sum of three active neutrino masses ($\sum m_{\nu}$) in various cosmological models. We use the following baseline datasets: CMB temperature data from Planck 2015, BAO measurements from SDSS-III BOSS DR12, the newly released SNe Ia dataset from Pantheon Sample, and a prior on the optical depth to reionization from 2016 Planck Intermediate results. We constrain cosmological parameters in $\Lambda CDM$ model with 3 massive active neutrinos. For this $\Lambda CDM+\sum m_{\nu}$ model we find a upper bound of $\sum m_{\nu} <$ 0.152 eV at 95$\%$ C.L. Adding the high-$l$ polarization data from Planck strengthens this bound to $\sum m_{\nu} <$ 0.118 eV, which is very close to the minimum required mass of $\sum m_{\nu} \simeq$ 0.1 eV for inverted hierarchy. This bound is reduced to $\sum m_{\nu} <$ 0.110 eV when we also vary r, the tensor to scalar ratio ($\Lambda CDM+r+\sum m_{\nu}$ model), and add an additional dataset, BK14, the latest data released from the Bicep-Keck collaboration. This bound is further reduced to $\sum m_{\nu} <$ 0.101 eV in a cosmology with non-phantom dynamical dark energy ($w_0 w_a CDM+\sum m_{\nu}$ model with $w(z)\geq -1$ for all $z$). Considering the $w_0 w_a CDM+r+\sum m_{\nu}$ model and adding the BK14 data again, the bound can be even further reduced to $\sum m_{\nu} <$ 0.093 eV. For the $w_0 w_a CDM+\sum m_{\nu}$ model without any constraint on $w(z)$, the bounds however relax to $\sum m_{\nu} <$ 0.276 eV. Adding a prior on the Hubble constant ($H_0 = 73.24\pm 1.74$ km/sec/Mpc) from Hubble Space Telescope (HST), the above mentioned bounds further improve to $\sum m_{\nu} <$ 0.117 eV, 0.091 eV, 0.085 eV, 0.082 eV, 0.078 eV and 0.247 eV respectively. This substantial improvement is mostly driven by a more than 3$\sigma$ tension between Planck 2015 and HST measurements of $H_0$ and should be taken cautiously. (abstract abridged)

The annihilation of Weakly Interactive Massive Particles (WIMP) in the centre of the sun could give rise to neutrino fluxes. We study the prospects of searching for these neutrinos at the upcoming Iron CALorimeter (ICAL) detector to be housed at the India-based Neutrino Observatory (INO). We perform ICAL simulations to obtain the detector efficiencies and resolutions in order to simulate muon events in ICAL due to neutrinos coming from annihilation of WIMP in the mass range $m_\chi = (3-100)$ GeV. The atmospheric neutrinos pose a major background for these indirect detection studies and can be reduced using the fact that the signal comes only from the direction of the sun. For a given WIMP mass, we find the opening angle $\theta_{90}$ such that 90 \% of the signal events are contained within this angle and use this cone-cut criteria to reduce the atmospheric neutrino background. The reduced background is then weighted by the solar exposure function at INO to obtain the final background spectrum for a given WIMP mass. We perform a $\chi^2$ analysis and present expected exclusion regions in the $\sigma_{SD}-m_\chi$ and $\sigma_{SI}-m_\chi$, where $\sigma_{SD}$ and $\sigma_{SI}$ are the WIMP-nucleon Spin-Dependent (SD) and Spin-Independent (SI) scattering cross-section, respectively. For a 10 years exposure and $m_\chi=25$ GeV, the expected 90 \% C.L. exclusion limit is found to be $\sigma_{SD} < 6.87\times 10^{-41}$ cm$^2$ and $\sigma_{SI} < 7.75\times 10^{-43}$ cm$^2$ for the $\tau^{+} \tau^{-}$ annihilation channel and $\sigma_{SD} < 1.14\times 10^{-39}$ cm$^2$ and $\sigma_{SI} < 1.30\times 10^{-41}$ cm$^2$ for the $b~\bar b $ channel, assuming 100 \% branching ratio for each of the WIMP annihilation channel.

We discuss inflation and dark matter in the inert doublet model coupled non-minimally to gravity where the inert doublet is the inflaton and the neutral scalar part of the doublet is the dark matter candidate. We calculate the various inflationary parameters like $n_s$, $r$ and $P_s$ and then proceed to the reheating phase where the inflaton decays into the Higgs and other gauge bosons which are non-relativistic owing to high effective masses. These bosons further decay or annihilate to give relativistic fermions which are finally responsible for reheating the universe. At the end of the reheating phase, the inert doublet which was the inflaton enters into thermal equilibrium with the rest of the plasma and its neutral component later freezes out as cold dark matter with a mass of about 2 TeV.

We investigate the impact of non-standard neutrino interactions (NSIs) on atmospheric neutrinos using the proposed PINGU experiment. In particular, we focus on the matter NSI parameters $\varepsilon_{\mu\tau}$ and $|\varepsilon_{\tau\tau} - \varepsilon_{\mu\mu}|$ that have previously been constrained by the Super-Kamiokande experiment. First, we present approximate analytical formulas for the difference of the muon neutrino survival probability with and without the above-mentioned NSI parameters. Second, we calculate the atmospheric neutrino events at PINGU in the energy range (2-100) GeV, which follow the trend outlined on probability level. Finally, we perform a statistical analysis of PINGU. Using three years of data, we obtain bounds from PINGU given by $-0.0043~(-0.0048) < \varepsilon_{\mu\tau} < 0.0047~(0.0046)$ and $-0.03~(-0.016) < \varepsilon_{\tau\tau} < 0.017~(0.032)$ at 90 % confidence level for normal (inverted) neutrino mass hierarchy, which improve the Super-Kamiokande bounds by one order of magnitude. In addition, we show the expected allowed contour region in the $\varepsilon_{\mu\tau}$-$\varepsilon_{\tau\tau}$ plane if NSIs exist in Nature and the result suggests that there is basically no correlation between $\varepsilon_{\mu\tau}$ and $\varepsilon_{\tau\tau}$.

We propose the liquid-scintillator detector LENA (Low Energy Neutrino Astronomy) as a next-generation neutrino observatory on the scale of 50 kt. The outstanding successes of the Borexino and KamLAND experiments demonstrate the large potential of liquid-scintillator detectors in low-energy neutrino physics. LENA's physics objectives comprise the observation of astrophysical and terrestrial neutrino sources as well as the investigation of neutrino oscillations. In the GeV energy range, the search for proton decay and long-baseline neutrino oscillation experiments complement the low-energy program. Based on the considerable expertise present in European and international research groups, the technical design is sufficiently mature to allow for an early start of detector realization.

Collective flavor oscillations are known to bring multiple splits in the supernova (SN) neutrino and antineutrino spectra. These spectral splits depend not only on the mass hierarchy of the neutrinos but also on the initial relative flux composition. Observation of spectral splits in a future galactic supernova signal is expected to throw light on the mass hierarchy pattern of the neutrinos. However, since the Diffuse Supernova Neutrino Background (DSNB) comprises of a superposition of neutrino fluxes from all past supernovae, and since different supernovae are expected to have slightly different initial fluxes, it is pertinent to check if the hierarchy dependent signature of collective oscillations can survive this averaging of the flux spectra. Since the actual distribution of SN with initial relative flux spectra of the neutrinos and antineutrinos is unknown, we assume a log-normal distribution for them. We study the dependence of the hierarchy sensitivity to the mean and variance of the log-normal distribution function. We find that the hierarchy sensitivity depends crucially on the mean value of the relative initial luminosity. The effect of the width is to reduce the hierarchy sensitivity for all values of the mean initial relative luminosity. We find that in the very small mixing angle ($\theta_{13}$) limit considering only statistical errors even for very moderate values of variance, there is almost no detectable hierarchy sensitivity if the mean relative luminosities of $\nu_e$ and $\bar\nu_e$ are greater than 1.

We calculate the expected galactic supernova neutrino signal at large next-generation underground detectors. At different epochs after the explosion, the primary fluxes can be quite different. For these primary neutrino fluxes, spectral splits induced by collective neutrino flavor transformations can arise for either mass hierarchy in both neutrino and antineutrino channels. We classify flux models according to the nature and number of these splits, and calculate the observable electron-neutrino and electron-antineutrino spectra at Earth, taking into account subsequent matter effects. We find that some of the spectral splits could occur sufficiently close to the peak energies to produce significant distortions in the observable SN neutrino signal. The most striking signature of this effect would be presence of peculiar energy dependent modulations associated with Earth matter crossing, present only in portions of the SN neutrino energy spectra demarcated by spectral splits. These signatures at proposed large water Cherenkov, scintillation, and liquid Argon detectors could give hints about the primary SN neutrino fluxes, as well as on the neutrino mass hierarchy and the mixing angle theta_{13}.

Neutrino-neutrino interactions inside core-collapse supernovae may give rise to collective flavor oscillations resulting in swap between flavors. These oscillations depend on the initial energy spectra, and relative fluxes or relative luminosities of the neutrinos. It has been observed that departure from energy equipartition among different flavors can give rise to one or more sharp spectral swap over energy, termed as splits. We study the occurrence of splits in the neutrino and antineutrino spectra, varying the initial relative fluxes for different models of initial energy spectrum, in both normal and inverted hierarchy. These initial relative flux variations give rise to several possible split patterns whereas variation over different models of energy spectra give similar results. We explore the effect of these spectral splits on the electron fraction, $Y_e$, that governs r-process nucleosynthesis inside supernovae. Since spectral splits modify the electron neutrino and antineutrino spectra in the region where r-process is postulated to happen, and since the pattern of spectral splits depends on the initial conditions of the spectra and the neutrino mass hierarchy, we show that the condition $Y_e < 0.5$ required for successful r-process nucleosynthesis will lead to constraints on the initial spectral conditions, for a given neutrino mass hierarchy.

We examine how light neutrinos coming from distant active galactic nuclei (AGN) and similar high energy sources may be used as tools to probe non-standard physics. In particular we discuss how studying the energy spectra of each neutrino flavour coming from such distant sources and their distortion relative to each other may serve as pointers to exotic physics such as neutrino decay, Lorentz symmetry violation, pseudo-Dirac effects, CP and CPT violation and quantum decoherence. This allows us to probe hitherto unexplored ranges of parameters for the above cases, for example lifetimes in the range $ 10^{-3}-10^{4} $ s/eV for the case of neutrino decay. We show that standard neutrino oscillations ensure that the different flavours arrive at the earth with similar shapes even if their flavour spectra at source may differ strongly in both shape and magnitude. As a result, observed differences between the spectra of various flavours at the detector would be signatures of non-standard physics altering neutrino fluxes during propagation rather than those arising during their production at source. Since detection of ultra-high energy (UHE) neutrinos is perhaps imminent, it is possible that such differences in spectral shapes will be tested in neutrino detectors in the near future. To that end, using the IceCube detector as an example, we show how our results translate to observable shower and muon-track event rates.

We study spectral distortions of diffuse ultra-high energy (UHE) neutrino flavour fluxes resulting due to physics beyond the Standard Model (SM). Even large spectral differences between flavours at the source are massaged into a common shape at earth by SM oscillations, thus, any significant observed spectral differences are an indicator of new physics present in the oscillation probability during propagation. Neutrino decay and Lorentz symmetry violation (LV) are examples, and result in significant distortion of the fluxes and the well-known bounds on them, which may allow UHE detectors to probe lifetimes, the mass hierarchy and LV parameters over a broad range.

We parameterize the initial flux composition of high energy astrophysical neutrinos as (Phi_e^0 : Phi_mu^0 : Phi_tau^0) = (1 : n : 0), where n characterizes the source. All usually assumed neutrino sources appear as limits of this simple parametrization. We investigate how precise neutrino telescopes can pin down the value of n. We furthermore show that there is a neutrino mixing scenario in which the ratio of muon neutrinos to the other neutrinos takes a constant value regardless of the initial flux composition. This occurs when the muon neutrino survival probability takes its minimal allowed value. The phenomenological consequences of this very predictive neutrino mixing scenario are given.

Collective flavor oscillations driven by neutrino-neutrino self interaction inside core-collapse supernovae have now been shown to bring drastic changes in the resultant neutrino fluxes. This would in turn significantly affect the diffuse supernova neutrino background (DSNB), created by all core-collapse supernovae that have exploded in the past. In view of these collective effects, we re-analyze the potential of detecting the DSNB in currently running and planned large-scale detectors meant for detecting both electron neutrinos and antineutrinos. The next generation detectors should be able to observe DSNB fluxes. Under certain conducive conditions, one could learn about neutrino parameters. For instance, it might be possible to determine the neutrino mass hierarchy, even if theta_{13} is almost zero.

Among all neutrino mixing parameters, the atmospheric neutrino mixing angle theta_{23} introduces the strongest variation on the flux ratios of ultra high energy neutrinos. We investigate the potential of these flux ratio measurements at neutrino telescopes to constrain theta_{23}. We consider astrophysical neutrinos originating from pion, muon-damped and neutron sources and make a comparative study of their sensitivity reach to theta_{23}. It is found that neutron sources are most favorable for testing deviations from maximal theta_{23}. Using a chi^2 analysis, we show in particular the power of combining (i) different flux ratios from the same type of source, and also (ii) combining flux ratios from different astrophysical sources. We include in our analysis ``impure'' sources, i.e., deviations from the usually assumed initial (1 : 2 : 0), (0 : 1 : 0) or (1 : 0 : 0) flux compositions.

The MiniBooNE results have still not been able to comprehensively rule out the oscillation interpretation of the LSND experiment. So far the so-called short baseline experiments with energy in the MeV range and baseline of few meters have been probing the existence of sterile neutrinos. We show how signatures of these extra sterile states could be obtained in TeV energy range atmospheric neutrinos travelling distances of thousands of kilometers. Atmospheric neutrinos in the TeV range would be detected by the upcoming neutrino telescopes. Of course vacuum oscillations of these neutrinos would be very small. However, we show that resonant matter effects inside the Earth could enhance these very tiny oscillations into near-maximal transitions, which should be hard to miss. We show that imprint of sterile neutrinos could be unambiguously obtained in this high energy atmospheric neutrino event sample. Not only would neutrino telescopes tell the presence of sterile neutrinos, it should also be possible for them to distinguish between the different possible mass and mixing scenarios with additional sterile states.

We expound the impact of extra sterile species on the ultra high energy neutrino fluxes in neutrino telescopes. We use three types of well-known flux ratios and compare the values of these flux ratios in presence of sterile neutrinos, with those predicted by deviation from the tribimaximal mixing scheme. We show that in the upcoming neutrino telescopes, its easy to confuse between the signature of sterile neutrinos with that of the deviation from tribimaximal mixing. We also show that if the measured flux ratios acquire a value well outside the range predicted by the standard scenario with three active neutrinos only, it might be possible to tell the presence of extra sterile neutrinos by observing ultra high energy neutrinos in future neutrino telescopes.

The implications of extra sterile neutrinos for the Double Chooz experiment is expounded. The so-called ``3+2'' mass spectrum with 2 sterile neutrinos mixed with the active ones, is still allowed by the global neutrino data including MiniBooNE. We probe its impact on the resultant reactor antineutrino signal at the near and far detector of the Double Chooz experiment. The oscillations driven by the additional mass squared difference due to the sterile states bring an energy independent constant suppression at both the near and far detectors. We study to what extent the measurement of $\theta_{13}$ would get affected due to the presence of sterile mixing. We also give the projected sensitivity that Double Chooz will have to constrain the extra mixing angles associated with the sterile states.

The time dependent neutrino oscillation signals due to the passage of a shock wave through the supernovae are analyzed for the case of three active neutrinos and also for the case that there are two additional sterile neutrinos. It is shown that, even without flavour identification and energy measurement, detailed information about the masses and mixing angles of the neutrinos may be obtained with a detector with excellent time resolution such as IceCube. Such a signal would also give important information about the nature of the shock wave within the supernovae.

We use the detection of neutrinos from a future galactic type-II supernova event in a water Cerenkov detector like Super-Kamiokande to constrain the possible violation of spin-statistics by neutrinos resulting in their obeying a mixed statistics instead of Fermi-Dirac.

Pioneering atmospheric muon neutrino experiments have demonstrated the near-maximal magnitude of the flavor mixing angle $\theta_{23}$. But the precise value of the deviation $D \equiv 1/2 - \sin^2 \theta_{23}$ from maximality (if nonzero) needs to be known, being of great interest -- especially to builders of neutrino mass and mixing models. We quantitatively investigate in a three generation framework the feasibility of determining $D$ in a statistically significant manner from studies of the atmospheric $\nu_\mu,\bar\nu_\mu$ survival probability including both vacuum oscillations and matter effects. We show how this determination will be sharpened by considering the up-down ratios of observed $\nu_\mu$- and $\bar\nu_\mu$-induced events and the differences of these ratios in specified energy and zenith angle bins. We consider 1 Megaton year of exposure to a magnetized iron calorimeter such as the proposed INO detector ICAL, taking into account both energy and zenith angle resolution functions. The sensitivity of such an exposure and the dependence of the determination of $D$ on the concerned oscillation parameters are discussed in detail. The vital use of matter effects in fixing the octant of $\theta_{23}$ is highlighted.

We discuss the possibilities of high precision measurement of the solar neutrino mixing angle $\theta_\odot \equiv \theta_{12}$ in solar and reactor neutrino experiments. The improvements in the determination of $\sin^2\theta_{12}$, which can be achieved with the expected increase of statistics and reduction of systematic errors in the currently operating solar and KamLAND experiments, are summarised. The potential of LowNu $\nu-e$ elastic scattering experiment, designed to measure the $pp$ solar neutrino flux, for high precision determination of $\sin^2\theta_{12}$, is investigated in detail. The accuracy in the measurement of $\sin^2\theta_{12}$, which can be achieved in a reactor experiment with a baseline $L \sim (50-70)$ km, corresponding to a Survival Probability MINimum (SPMIN), is thoroughly studied. We include the effect of the uncertainty in the value of $\sin^2\theta_{13}$ in the analyses. A LowNu measurement of the $pp$ neutrino flux with a 1% error would allow to determine $\sin^2\theta_{12}$ with an error of 14% (17%) at 3$\sigma$ from a two-generation (three-generation) analysis. The same parameter $\sin^2\theta_{12}$ can be measured with an uncertainty of 2% (6%) at 1$\sigma$ (3$\sigma$) in a reactor experiment with $L \sim60 $ km, statistics of $\sim$60 GWkTy and systematic error of 2%. For the same statistics, the increase of the systematic error from 2% to 5% leads to an increase in the uncertainty in $\sin^2\theta_{12}$ from 6% to 9% at 3$\sigma$. The inclusion of the $\sin^2\theta_{13}$ uncertainty in the analysis changes the error on $\sin^2\theta_{12}$ to 3% (9%). The effect of $\sin^2\theta_{13}$ uncertainty on the $\sin^2\theta_{12}$ measurement in both types of experiments is considerably smaller than naively expected.

We investigate the impact of the 766.3 Ty KamLAND spectrum data on the determination of the solar neutrino oscillation parameters. We show that the observed spectrum distortion in the KamLAND experiment firmly establishes $\Delta m^2_{21}$ to lie in the low-LMA solution region. The high-LMA solution is excluded at more than 4$\sigma$ by the global solar neutrino and KamLAND spectrum data. The maximal solar neutrino mixing is ruled out at $6\sigma$ level. The $3\sigma$ allowed region in the $\Delta m^2_{21}-\sin^2\theta_{12}$ plane is found to be remarkably stable with respect to leaving out the data from one of the solar neutrino experiments from the global analysis. We perform a three flavor neutrino oscillation analysis of the global solar neutrino and KamLAND spectrum data as well. The $3\sigma$ upper limit on $\sin^2\theta_{13}$ is found to be $ \sin^2\theta_{13} <0.055$. We derive predictions for the CC to NC event rate ratio and day-night (D-N) asymmetry in the CC event rate, measured in the SNO experiment, and for the suppression of the event rate in the BOREXINO and LowNu experiments. Prospective high precision measurements of the solar neutrino oscillation parameters are also discussed.

We explore the potential of measuring the solar neutrino oscillation parameters in the proposed gadolinium loaded Super-Kamiokande (SK-Gd) detector. Gadolinium dissolved in water can detect neutrons much more efficiently than pure water. This imparts the detector the ability to observe electron type antineutrinos, transforming Super-Kamiokande into a huge reactor antineutrino detector with an event rate approximately 43 times higher than that observed in KamLAND. We simulate the reactor antineutrino data expected in this high statistics detector. We use these prospective data to study the precision with which the solar neutrino oscillation parameters, $\Delta m^2_{\odot}$ and $\sin^2\theta_{\odot}$, can be determined i) with the SK-Gd detector, and ii) by combining the SK-Gd data with the global data on solar neutrino oscillations. For comparison and completeness the allowed regions of $\Delta m^2_{\odot}$ and $\sin^2\theta_{\odot}$, expected to be obtained from the data of the solar neutrino and KamLAND experiments, are also presented.

The physics implications of the just published salt phase data from the SNO experiment are examined. The effect of these data on the allowed ranges of the solar neutrino oscillation parameters, $\Delta_m^2_{21}$ and $\sin^2\theta_{12}$, are studied in the cases of two- and three- neutrino mixing. In the latter case we derive an upper limit on the angle $\theta_{13}$. Constraints on the solar $\nu_e$ transitions into a mixture of active and sterile neutrinos are also presented. Finally, we give predictions for the day-night asymmetry in the SNO experiment, for the event rate in the BOREXINO and LowNu experiments, and discuss briefly the constraints on the solar neutrino oscillation parameters which can be obtained with prospective KamLAND data.

A new reactor power plant Shika-2, with a power of approximately 4 GW and at a distance of about 88 km from the KamLAND detector is scheduled to start operating in March 2006. We study the impact of the $\bar\nu_e$ flux from this reactor on the sensitivity of the KamLAND experiment to the solar neutrino oscillation parameters. We present results on prospective determination of $\Delta m^2_\odot$ and $\sin^2\theta_\odot$ using the combined data from KamLAND and the solar neutrino experiments, including the effect of the Shika-2 contribution to the KamLAND signal and the latest data from the salt enriched phase of the SNO experiment. We find that contrary to the expectations, the addition of the Shika-2 reactor flux does not improve the $\sin^2\theta_\odot$ sensitivity of KamLAND, while the ambiguity in $\Delta m^2_\odot$ measurement may even increase, as a result of the averaging effect between Kashiwazaki and the Shika-2 reactor contributions to the KamLAND signal.

We investigate the potential of 3 kiloTon-years(kTy) of KamLAND data to further constrain the $\Delta m^2$ and $\tan^2\theta$ values compared to those presently allowed by existing KamLAND and global solar data. We study the extent, dependence and characteristics of this sensitivity in and around the two parts of the LMA region that are currently allowed. Our analysis with 3 kTy simulated spectra shows that KamLAND spectrum data by itself can constrain $\Delta m^2$ with high precision. Combining the spectrum with global solar data further tightens the constraints on allowed values of $\tan^2\theta$ and $\Delta m^2$. We also study the effects of future neutral current data with a total error of 7% from the Sudbury Neutrino Observatory. We find that these future measurements offer the potential of considerable precision in determining the oscillation parameters (specially the mass parameter).

We discuss the physics potential of intermediate $L \sim 20 \div 30$ km baseline experiments at reactor facilities, assuming that the solar neutrino oscillation parameters $\Delta m^2_{\odot}$ and $\theta_{\odot}$ lie in the high-LMA solution region. We show that such an intermediate baseline reactor experiment can determine both $\Delta m^2_{\odot}$ and $\theta_{\odot}$ with a remarkably high precision. We perform also a detailed study of the sensitivity of the indicated experiment to $\Delta m^2_{\rm atm}$, which drives the dominant atmospheric $\nu_{\mu}$ ($\bar{\nu}_{\mu}$) oscillations, and to $\theta$ - the neutrino mixing angle limited by the data from the CHOOZ and Palo Verde experiments. We find that this experiment can improve the bounds on $\sin^2\theta$. If the value of $\sin^2\theta$ is large enough, $\sin^2\theta \gtap 0.02$, the energy resolution of the detector is sufficiently good and if the statistics is relatively high, it can determine with extremely high precision the value of $\Delta m^2_{\rm atm}$. We also explore the potential of the intermediate baseline reactor neutrino experiment for determining the type of the neutrino mass spectrum, which can be with normal or inverted hierarchy. We show that the conditions under which the type of neutrino mass hierarchy can be determined are quite challenging, but are within the reach of the experiment under discussion.

However the allowed range of $\theta_{\odot}$ ($\theta_{12}$) did not change much with the inclusion of the KamLAND results. In this paper we probe if future data from KamLAND can increase the accuracy of the allowed range in $\theta_{\odot}$ and conclude that even after 3 kton-year of statistics and most optimistic error estimates, KamLAND {\it may} find it hard to significantly improve the bounds on the mixing angle obtained from the solar neutrino data. We discuss the $\theta_{12}$ sensitivity of the survival probabilities in matter (vacuum) as is relevant for the solar (KamLAND) experiments. We find that the presence of matter effects in the survival probabilities for $^8B$ neutrinos give the solar neutrino experiments SK and SNO an edge over KamLAND, as far as $\theta_{12}$ sensitivity is concerned, particularly near maximal mixing. Among solar neutrino experiments we identify SNO as the most promising candidate for constraining $\theta_{12}$ and make a projected sensitivity test for the mixing angle by reducing the error in the neutral current measurement at SNO. Finally we argue that the most accurate bounds on $\theta_{12}$ can be achieved in a reactor experiment, if the corresponding baseline and energy can be tuned to a minimum in the survival probability. We propose a new reactor experiment which can give the value of $\tan^2\theta_{12}$ to within 14%. We also discuss the future Borexino and LowNu experiments.

The first results from the KamLAND experiment have provided confirmational evidence for the Large Mixing Angle (LMA) MSW solution to the solar neutrino problem. We do a global analysis of solar and the recently announced KamLAND data (both rate and spectrum) and investigate its effect on the allowed region in the $\Delta m^2-\tan^2\theta$ plane. The best-fit from a combined analysis which uses the KamLAND rate plus global solar data comes at $\Delta m^2 = 6.06 \times 10^{-5}$ eV $^2$ and $\tan^2\theta=0.42$, very close to the global solar best-fit, leaving a large allowed region within the global solar LMA contour. The inclusion of the KamLAND spectral data in the global fit gives a best-fit $\Delta m^2 = 7.15 \times 10^{-5}$ eV $^2$ and $\tan^2\theta=0.42$ and constrains the allowed areas within LMA, leaving essentially two allowed zones. Maximal mixing though allowed by the KamLAND data alone is disfavored by the global solar data and remains disallowed at about $3\sigma$. The LOW solution is now ruled out at about 5$\sigma$ w.r.t. the LMA solution.

We calculate the time delays of neutrinos emitted in gamma ray bursts due to the effects of neutrino mass and quantum gravity using a time dependent Hubble constant which can significantly change the naive results presented hitherto in the literature for large redshifts, and gives some sensitivity to the details of dark energy. We show that the effects of neutrino mass, quantum gravity and dark energy may be disentangled by using low energy neutrinos to study neutrino mass, high energy neutrinos to study quantum gravity, and large redshifts to study dark energy. From low energy neutrinos one may obtain direct limits on neutrino masses of order 10^{-3} eV, and distinguish a neutrino mass hierarchy from an inverted mass hierarchy. From ultra-high energy neutrinos the sensitivity to the scale of quantum gravity can be pushed up to E_{QG} ~ 5 times 10^{30} GeV. By studying neutrinos from GRBs at large redshifts a cosmological constant could be distinguished from quintessence.

We perform model independent and model dependent analyses of solar neutrino data including the neutral current event rate from SNO. The inclusion of the first SNO NC data in the model independent analysis determines the allowed ranges of $^{8}{B}$ flux normalisation and the $\nu_e$ survival probability more precisely than what was possible from the SK and SNO CC combination. We perform global $\nu_e-\nu_{active}$ oscillation analyses of solar neutrino data using the NC rate instead of the SSM prediction for the $^{8}{B}$ flux, in view of the large uncertainty in the latter. The LMA gives the best solution, while the LOW solution is allowed only at the $3\sigma$ level.

We examine the solar neutrino problem in the context of the realistic three neutrino mixing scenario including the SNO charged current (CC) rate. The two independent mass squared differences $\Delta m^2_{21}$ and $\Delta m^2_{31} \approx \Delta m^2_{32}$ are taken to be in the solar and atmospheric ranges respectively. We incorporate the constraints on $\Delta$m$^2_{31}$ as obtained by the SuperKamiokande atmospheric neutrino data and determine the allowed values of $\Delta m^2_{21}$, $\theta_{12}$ and $\theta_{13}$ from a combined analysis of solar and CHOOZ data. Our aim is to probe the changes in the values of the mass and mixing parameters with the inclusion of the SNO data as well as the changes in the two-generation parameter region obtained from the solar neutrino analysis with the inclusion of the third generation. We find that the inclusion of the SNO CC rate in the combined solar + CHOOZ analysis puts a more restrictive bound on $\theta_{13}$. Since the allowed values of $\theta_{13}$ are constrained to very small values by the CHOOZ experiment there is no qualitative change over the two generation allowed regions in the $\Delta m^2_{21} - \tan^2 \theta_{12}$ plane. The best-fit comes in the LMA region and no allowed area is obtained in the SMA region at 3$\sigma$ level from combined solar and CHOOZ analysis.

In this thesis we explore the implications of neutrino oscillations in the context of the solar neutrino data, the atmospheric neutrino data and results from the terrestrial accelerator/reactor neutrino oscillation experiments. We perform comprehensive $\chi^2$ analysis for the global solar neutrino data including SNO and present our results for two flavor $\nu_e-\nu_{active}$ and $\nu_e-\nu_{sterile}$ oscillations. We explore the apparent conflict between the flat SK spectrum and the non-monotonic energy dependence of the solar neutrino survival probability inferred from the global data on total rates. We analyze the SK atmospheric neutrino data in the framework of two and three generations. For the three-generation analysis we simultaneously take into account the results of the terrestrial neutrino experiments. We consider the possibility of unstable neutrinos and explore the viability of this decay model as a solution to the atmospheric neutrino anomaly. We make quantitative predictions for the number of events recorded in the water Cerenkov detectors due a galactic supernova and examine the signatures of neutrino mass and mixing in the resultant neutrino signal.

A seismic model for the Sun calculated using the accurate helioseismic data predicts a lower $^{8}{B}$ neutrino flux as compared to the standard solar model (SSM). However, there persists a discrepancy between the predicted and measured neutrino fluxes and it seems necessary to invoke neutrino oscillations to explain the measurements. In this work, we have performed a global, unified oscillation analysis of the latest solar neutrino data (including the results of SNO charged current rate) using the seismic model fluxes as theoretical predictions. We determine the best-fit values of the neutrino oscillation parameters and the $\chi^2_{\mathrm min}$ for both $\nu_e-\nu_{\mathrm active}$ and $\nu_e -\nu_{\mathrm sterile}$ cases and present the allowed parameter regions in the $\Delta m^2 - \tan^2 \theta$ plane for $\nu_e-\nu_{\mathrm active}$ transition. The results are compared with those obtained using the latest SSM by Bahcall and his collaborators.

We investigate the implications of the SNO charged-current (CC) and electron scattering (ES) measurements of solar \br neutrino fluxes for neutrino mass and mixing parameters by performing a global and unified $\chi^2$ analysis of the solar neutrino data in the framework of two neutrino mixing. We consider both $\nu_e -\nu_{active}$ and $\nu_e -\nu_{sterile}$ solutions and perform (i) analysis of the total rates data of Cl, Ga, SK and SNO experiments and (ii) global analysis including the total rates data, the recoil electron spectrum data of SK and the CC spectrum observed at SNO. For the $\nu_e-\nu_{active}$ case the inclusion of the SNO results in the analysis of the total rates reduces(enhances) the goodness-of-fit (GOF) of the SMA(LMA) solution. The flat spectrum observed at SK further favours the LMA solution over the SMA solution and no allowed area is obtained in the SMA region at 3$\sigma$ level from the global analysis. For the $\nu_e -\nu_{sterile}$ case, with the inclusion of the SNO results, all the solutions are disfavoured with a probability of more than 99% from the total rates analysis while for the global analysis the GOF of these become much worse.

The neutrinos and antineutrinos of all the three flavors released from a galactic supernova will be detected in the water Cerenkov detectors. We show that even though the neutral current interaction is flavor blind, and hence neutrino flavor mixing cannot alter the total neutral current signal in the detector, it can have a non-trivial impact on the delay of massive neutrinos and alters the neutral current event rate as a function of time. We have suggested various variables of the neutral and charged current events that can be used to study this effect. In particular the ratio of charged to neutral current events can be used at early times while the ratio of the energy moments for the charged to the neutral current events can form useful diagnostic tools even at late times to study neutrino mass and mixing.

The neutrinos from galactic supernovae can be detected by the Sudbury Neutrino Observatory and the Super-Kamiokande. The effect of neutrino mass can show up in the observed neutrino signal by (i) delay in the time of flight and (ii) distortion of the neutrino energy spectrum due to flavor mixing. We discuss a combination of the two effects for both charged and neutral current processes in the detectors for realistic three flavor scenarios and show that neutrino flavor mixing can lead to non-trivial changes in the event rate as a function of time.