
The upcoming 50 kt magnetized iron calorimeter (ICAL) detector at the
Indiabased Neutrino Observatory (INO) is designed to study the atmospheric
neutrinos and antineutrinos separately over a wide range of energies and path
lengths. The primary focus of this experiment is to explore the Earth matter
effects by observing the energy and zenith angle dependence of the atmospheric
neutrinos in the multiGeV range. This study will be crucial to address some of
the outstanding issues in neutrino oscillation physics, including the
fundamental issue of neutrino mass hierarchy. In this document, we present the
physics potential of the detector as obtained from realistic detector
simulations. We describe the simulation framework, the neutrino interactions in
the detector, and the expected response of the detector to particles traversing
it. The ICAL detector can determine the energy and direction of the muons to a
high precision, and in addition, its sensitivity to multiGeV hadrons increases
its physics reach substantially. Its charge identification capability, and
hence its ability to distinguish neutrinos from antineutrinos, makes it an
efficient detector for determining the neutrino mass hierarchy. In this report,
we outline the analyses carried out for the determination of neutrino mass
hierarchy and precision measurements of atmospheric neutrino mixing parameters
at ICAL, and give the expected physics reach of the detector with 10 years of
runtime. We also explore the potential of ICAL for probing new physics
scenarios like CPT violation and the presence of magnetic monopoles.

In this paper, we discuss P(n), the number of ways in which a given integer n
may be written as a sum of primes. In particular, an asymptotic form P_as(n)
valid for n towards infinity is obtained analytically using standard techniques
of quantum statistical mechanics. First, the bosonic partition function of
primes, or the generating function of unrestricted prime partitions in number
theory, is constructed. Next, the density of states is obtained using the
saddlepoint method for Laplace inversion of the partition function in the
limit of large n. This directly gives the asymptotic number of prime partitions
P_as(n). The leading term in the asymptotic expression grows exponentially as
sqrt[n/ln(n)] and agrees with previous estimates. We calculate the
nexttoleading order term in the exponent, porportional to ln[ln(n)]/ln(n),
and show that an earlier result in the literature for its coefficient is
incorrect. Furthermore, we also calculate the next higher order correction,
proportional to 1/ln(n) and given in Eq.(43), which so far has not been
available in the literature. Finally, we compare our analytical results with
the exact numerical values of P(n) up to n \sim 8 10^6. For the highest values,
the remaining error between the exact P(n) and our P_as(n) is only about half
of that obtained with the leadingorder (LO) approximation. But we also show
that, unlike for other types of partitions, the asymptotic limit for the prime
partitions is still quite far from being reached even for n \sim 10^7.

We consider the case of light dark matter ($\sim 10$ GeV). We discuss a
simple $Z_2$ model of scalar selfinteracting dark matter, as well as a related
model of unstable longlived dark matter which can explain the anomalous Kolar
events observed decades ago.

The results of a Monte Carlo simulation study of the hadron energy response
for the magnetized Iron CALorimeter detector, ICAL, proposed to be located at
the Indiabased Neutrino Observatory (INO) is presented. Using a GEANT4
modeling of the detector ICAL, interactions of atmospheric neutrinos with
target nuclei are simulated. The detector response to hadrons propagating
through it is investigated using the hadron hit multiplicity in the active
detector elements. The detector response to charged pions of fixed energy is
studied first, followed by the average response to the hadrons produced in
atmospheric neutrino interactions using events simulated with the NUANCE event
generator. The shape of the hit distribution is observed to fit the Vavilov
distribution, which reduces to a Gaussian at high energies. In terms of the
parameters of this distribution, we present the hadron energy resolution as a
function of hadron energy, and the calibration of hadron energy as a function
of the hit multiplicity. The energy resolution for hadrons is found to be in
the range 85% (for 1GeV)  36% (for 15 GeV).

We demonstrate that the (swave) geometric spectrum of the Efimov energy
levels in the unitary limit is generated by the radial motion of a primitive
periodic orbit (and its harmonics) of the corresponding classical system. The
action of the primitive orbit depends logarithmically on the energy. It is
shown to be consistent with an inversesquared radial potential with a lower
cutoff radius. The lowestorder WKB quantization, including the Langer
correction, is shown to reproduce the geometric scaling of the energy spectrum.
The (WKB) meansquared radii of the Efimov states scale geometrically like the
inverse of their energies. The WKB wavefunctions, regularized near the
classical turning point by Langer's generalized connection formula, are
practically indistinguishable from the exact wave functions even for the lowest
($n=0$) state, apart from a tiny shift of its zeros that remains constant for
large $n$.

It is suggested that for a fermi gas at unitarity, the twobody bond plays a
special role. We propose an equation of state using an ansatz relating the
interaction part of the $l$body cluster to its twobody counterpart. This
allows a parameterfree comparison with the recently measured equation of state
by the ENS group. The agreement between the two over a range of fugacity ($z<5$
for a homogeneous gas, and $z<10$ for the trapped gas) leads us to perform the
calculations of more sensitive quantities measured recently by the MIT group.

We investigate the particle and kineticenergy densities for a system of $N$
fermions confined in a potential $V(\bfr)$. In an earlier paper [J. Phys. A:
Math. Gen. {\bf 36}, 1111 (2003)], some exact and asymptotic relations
involving the particle density and the kineticenergy density locally, i.e. at
any given point $\bfr$, were derived for isotropic harmonic oscillators in
arbitrary dimensions. In this paper we show that these {\it local virial
theorems} (LVT) also hold exactly for linear potentials in arbitrary dimensions
and for the onedimensional box. We also investigate the validity of these LVTs
when they are applied to arbitrary smooth potentials. We formulate generalized
LVTs that are supported by a semiclassical theory which relates the density
oscillations to the closed nonperiodic orbits of the classical system. We test
the validity of these generalized theorems numerically for various local
potentials. Although they formally are only valid asymptotically for large
particle numbers $N$, we show that they practically are surprisingly accurate
also for moderate values of $N$.

We consider N particles interacting pairwise by an inverse square potential
in one dimension (CalogeroSutherlandMoser model). When trapped harmonically,
its classical canonical partition function for the repulsive regime is known in
the literature. We start by presenting a concise rederivation of this result.
The equation of state is then calculated both for the trapped and the
homogeneous gas. Finally, the classical limit of Wu's distribution function for
fractional exclusion statistics is obtained and we rederive the classical
virial expansion of the homogeneous gas using this distribution function.

We first give a short review of the ``localcurrent approximation'' (LCA),
derived from a general variation principle, which serves as a semiclassical
description of strongly collective excitations in finite fermion systems
starting from their quantummechanical meanfield ground state. We illustrate
it for the example of coupled translational and compressional dipole
excitations in metal clusters. We then discuss collective electronic dipole
excitations in C$_{60}$ molecules (Buckminster fullerenes). We show that the
coupling of the pure translational mode (``surface plasmon'') with
compressional volume modes in the semiclasscial LCA yields semiquantitative
agreement with microscopic timedependent density functional (TDLDA)
calculations, while both theories yield qualitative agreement with the recent
experimental observation of a ``volume plasmon''.

This paper exploits the connection between the quantum manyparticle density
of states and the partitioning of an integer in number theory. For $N$ bosons
in a one dimensional harmonic oscillator potential, it is well known that the
asymptotic (N > infinity) density of states is identical to the
HardyRamanujan formula for the partitions p(n), of a number n into a sum of
integers. We show that the same statistical mechanics technique for the density
of states of bosons in a powerlaw spectrum yields the partitioning formula for
p^s(n), the latter being the number of partitions of n into a sum of sth
powers of a set of integers. By making an appropriate modification of the
statistical technique, we are also able to obtain d^s(n) for distinct
partitions. We find that the distinct square partitions d^2(n) show pronounced
oscillations as a function of n about the smooth curve derived by us. The
origin of these oscillations from the quantum point of view is discussed. After
deriving the ErdosLehner formula for restricted partitions for the $s=1$ case
by our method, we generalize it to obtain a new formula for distinct restricted
partitions.

We investigate the particle and kinetic energy densities of harmonically
trapped fermion gases at zero temperature in arbitrary dimensions. We derive
analytically a differential equation connecting these densities, which so far
have been proven only in one or two dimensions, and give other interesting
relations involving several densities or the particle density alone. We show
that in the asymptotic limit of large particle numbers, the densities go over
into the semiclassical ThomasFermi (TF) densities. Hereby the Fermi energy to
be used in the TF densities is identified uniquely. We derive an analytical
expansion for the remaining oscillating parts and obtain very simple closed
forms for the leadingorder oscillating densities. Finally, we show that the
simple TF functional relation $\tau_{TF}[\rho]$ between kinetic and particle
density is fulfilled also for the asymptotic quantum densities $\tau(r)$ and
$\rho(r)$ including their leadingorder oscillating terms.

It is wellknown that the number fluctuation in the grand canonical ensemble,
which is directly proportional to the compressibility, diverges for an ideal
bose gas as T > 0. We show that this divergence is removed when the atoms
interact in one dimension through an inverse square twobody interaction. In
two dimensions, similar results are obtained using a selfconsistent
ThomasFermi (TF) model for a repulsive zerorange interaction. Both models may
be mapped on to a system of noninteracting particles obeying the HaldaneWu
exclusion statistics. We also calculate the number fluctuation from the ground
state of the gas in these interacting models, and compare the grand canonical
results with those obtained from the canonical ensemble.

We have analysed the recent results from the observation of charged current
\nu_e d \to e^ p p events from solar neutrinos by the Sudbury Neutrino
Observatory SNO assuming neutrino oscillations with three active flavours. The
data seem to prefer a low masssquared difference and large mixing angle
solution (the socalled LOW solution) in (12) parameter space. However, when
combined with the Gallium charged current interaction data from Gallex and GNO,
distinct (1\sigma) allowed regions corresponding to the large mixing angle
(LMA) and small mixing angle (SMA) appear while the LOW solution is disfavoured
upto 3\sigma standard deviation. The physical electron neutrino survival
probability corresponding to these best fit solutions are then determined and
analysed for their energy dependence.

We have identified some important and worthwhile physics opportunitites with
a possible neutrino detector located in India. Particular emphasis is placed on
the geographical advantage with a stress on the complimentary aspects with
respect to other neutrino detectors already in operation.

Signatures of neutrino and antineutrino signals from stellar collapse in
heavy water detectors are contrasted with those in water detectors. The effects
of mixing, especially due to the highly dense matter in the supernova core, are
studied. The mixing parameters used are those sets allowed by current
understanding of available neutrino data: from solar, atmospheric and
laboratory neutrino experiments. Signals at a heavy water detector, especially
the dominant charged current reactions on deuteron, are very sensitive to some
of these sets of allowed mixing parameters. Theoretical uncertainties on
supernova neutrino spectra notwithstanding, a combination of supernova
measurements with water and heavy water detectors may be able to distinguish
many of these mixing possibilities and thus help in ruling out many of them.

We consider a small and fixed number of fermions (bosons) in a trap. The
ground state of the system is defined at T=0. For a given excitation energy,
there are several ways of exciting the particles from this ground state. We
formulate a method for calculating the number fluctuation in the ground state
using microcanonical counting, and implement it for small systems of
noninteracting fermions as well as bosons in harmonic confinement. This exact
calculation for fluctuation, when compared with canonical ensemble averaging,
gives considerably different results, specially for fermions. This difference
is expected to persist at low excitation even when the fermion number in the
trap is large.

We study the effect of nonvanishing masses and mixings among neutrino
flavours on the detection of neutrinos from stellar collapse by a water
Cerenkov detector. We consider a framework in which there are four neutrino
flavours, including a sterile neutrino, whose mass squared differences and
mixings are constrained by the present understanding of solar, atmospheric and
laboratory neutrino detection. We also include the effects of high density
matter within the supernova core. Unlike in the three flavour scenario, we find
that the number of events due to the dominant process involving
electronantineutrinos changes dramatically for some allowed mixing parameters.
Furthermore, contributions from chargedcurrent scattering off oxygen nuclei in
the detector can be considerably enhanced due to flavour mixing. We also
present a comparison between the two possible scenarios, namely, when only
three active neutrino flavours are present and when they are accompanied by a
fourth sterile neutrino.

We study the effect of nonvanishing masses and mixings among neutrino
flavours on the detection of neutrinos from stellar collapse by a water
Cerenkov detector. We consider a realistic framework in which there are three
neutrino flavours whose mass squared differences and mixings are constrained by
the present understanding of solar and atmospheric neutrinos. We also include
the effects of high dense matter within the supernova core. We find that the
number of events due to the dominant process involving electronantineutrinos
may change dramatically for some allowed mixing parameters. Furthermore,
contributions from chargedcurrent scattering off oxygen atoms in the detector
can be considerably enhanced due to flavour mixing; such events have a distinct
experimental signature since they are backwardpeaked. Hence, mixing has a
nontrivial effect on the signature of neutrinos (and antineutrinos) from
stellar collapse.

We give a formulation of the single particle occupation probabilities for a
system of identical particles obeying fractional exclusion statistics of
Haldane. We first derive a set of constraints using an exactly solvable model
which describes an ideal exclusion statistics system and deduce the general
counting rules for occupancy of states obeyed by these particles. We show that
the problem of negative probabilities may be avoided with these new counting
rules.

Periodic orbit expressions for the density of states lead to spurious results
when directly used to calculate quantities of thermodynamic interest. This is
because the trace formula is usually valid only for large energies while the
calculations make use of it at all energies. We present a prescription for
circumventing this problem by isolating contributions which arise from the
inaccurate lowenergy behavior, and are spurious, from other subdominant
contributions which are physical and are not contained in the Weyl series. The
method is tested by analyzing the fermionic grand partition function for a disk
billiard and may readily be extended to other dynamical systems.

The classical ground state of a D dimensional many body system with two and
three body interactions is studied as a function of the strength of the three
body interaction. We prove exactly that beyond a critical strength of the three
body interaction, the classical ground state of the system is one in which all
the particles are on a line. The positions of the particles in this string
configuration are uniquely determined by the zeros of the Hermite polynomials.

We report on a study of a classical, finite system of confined particles in
two dimensions with a twobody repulsive interaction. We first develop a simple
analytical method to obtain equilibrium configurations and energies for few
particles. When the confinement is harmonic, we prove that the first transition
from a single shell occurs when the number of particles changes from five to
six. The shell structure in the case of an arbitrary number of particles is
shown to be independent of the strength of the interaction but dependent only
on its functional form. It is also independent of the magnetic field strength
when included. We further study the effect of the functional form of the
confinement potential on the shell structure. Finally we report some
interesting results when a threebody interaction is included, albeit in a
particular model.

We report on a study of a finite system of classical confined particles in
twodimensions in the presence of a uniform magnetic field and interacting via
a twobody repulsive potential. We develop a simple analytical method of
analysis to obtain ground state energies and configurations. We prove
analytically the minimum energy configurations are independent of the nature of
two body interaction and the magnetic field. In particular we prove that the
first transition from a single shell occurs when the number of particles
changes from five to six. These results are exact.

We discuss a manybody Hamiltonian with two and threebody interactions in
two dimensions introduced recently by Murthy, Bhaduri and Sen. Apart from an
analysis of some exact solutions in the manybody system, we analyze in detail
the twobody problem which is completely solvable. We show that the solution of
the twobody problem reduces to solving a known differential equation due to
Heun. We show that the twobody spectrum becomes remarkably simple for large
interaction strength and the level structure resembles that of the Landau
Levels. We also clarify the "ultraviolet" regularization which is needed to
define an inversesquare potential properly and discuss its implications for
our model.

We construct a new manybody Hamiltonian with two and threebody
interactions in two space dimensions and obtain its exact manybody ground
state for an arbitrary number of particles. This ground state has a novel
pairwise correlation. A class of exact solutions for the excited states is also
found. These excited states display an energy spectrum similar to the
CalogeroSutherland model in one dimension. The model reduces to an analog of
the wellknown trigonometric Sutherland model when projected on to a circular
ring.