• The ICAL Collaboration: Shakeel Ahmed, M. Sajjad Athar, Rashid Hasan, Mohammad Salim, S. K. Singh, S. S. R. Inbanathan (The American College), Venktesh Singh, V. S. Subrahmanyam, Shiba Prasad Behera, Vinay B. Chandratre, Nitali Dash, Vivek M. Datar, V. K. S. Kashyap, Ajit K. Mohanty, Lalit M. Pant, Animesh Chatterjee, Sandhya Choubey, Raj Gandhi, Anushree Ghosh, Deepak Tiwari (HRI, Allahabad), Ali Ajmi, S. Uma Sankar, Prafulla Behera, Aleena Chacko, Sadiq Jafer, James Libby, K. Raveendrababu, K. R. Rebin, D. Indumathi, K. Meghna, S. M. Lakshmi, M. V. N. Murthy, Sumanta Pal, G. Rajasekaran, Nita Sinha, Sanjib Kumar Agarwalla, Amina Khatun, Poonam Mehta, Vipin Bhatnagar, R. Kanishka, A. Kumar, J. S. Shahi, J. B. Singh, Monojit Ghosh, Pomita Ghoshal, Srubabati Goswami, Chandan Gupta, Sushant Raut, Sudeb Bhattacharya, Suvendu Bose, Ambar Ghosal, Abhik Jash, Kamalesh Kar, Debasish Majumdar, Nayana Majumdar, Supratik Mukhopadhyay, Satyajit Saha (Saha Inst. Nucl. Phys., Kolkata), B. S. Acharya, Sudeshna Banerjee, Kolahal Bhattacharya, Sudeshna Dasgupta, Moon Moon Devi, Amol Dighe, Gobinda Majumder, Naba K. Mondal, Asmita Redij, Deepak Samuel, B. Satyanarayana, Tarak Thakore, C. D. Ravikumar, A. M. Vinodkumar, Gautam Gangopadhyay, Amitava Raychaudhuri, Brajesh C. Choudhary, Ankit Gaur, Daljeet Kaur, Ashok Kumar, Sanjeev Kumar, Md. Naimuddin, Waseem Bari, Manzoor A. Malik, S. Krishnaveni, H. B. Ravikumar, C. Ranganathaiah, Saikat Biswas, Subhasis Chattopadhyay, Rajesh Ganai, Tapasi Ghosh, Y. P. Viyogi
    May 9, 2017 hep-ph, hep-ex, physics.ins-det
    The upcoming 50 kt magnetized iron calorimeter (ICAL) detector at the India-based 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 multi-GeV 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 multi-GeV 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 saddle-point 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 next-to-leading 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 leading-order (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 self-interacting dark matter, as well as a related model of unstable long-lived 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 India-based 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 (s-wave) 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 inverse-squared radial potential with a lower cut-off radius. The lowest-order WKB quantization, including the Langer correction, is shown to reproduce the geometric scaling of the energy spectrum. The (WKB) mean-squared 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 two-body 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 two-body counterpart. This allows a parameter-free 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 kinetic-energy 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 kinetic-energy 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 one-dimensional 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 non-periodic 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 pair-wise by an inverse square potential in one dimension (Calogero-Sutherland-Moser model). When trapped harmonically, its classical canonical partition function for the repulsive regime is known in the literature. We start by presenting a concise re-derivation 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 re-derive the classical virial expansion of the homogeneous gas using this distribution function.
  • We first give a short review of the ``local-current 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 quantum-mechanical mean-field 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 semi-quantitative agreement with microscopic time-dependent 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 many-particle 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 Hardy-Ramanujan 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 power-law spectrum yields the partitioning formula for p^s(n), the latter being the number of partitions of n into a sum of s-th 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 Erdos-Lehner 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 semi-classical Thomas-Fermi (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 leading-order 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 leading-order oscillating terms.
  • It is well-known 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 two-body interaction. In two dimensions, similar results are obtained using a self-consistent Thomas-Fermi (TF) model for a repulsive zero-range interaction. Both models may be mapped on to a system of non-interacting particles obeying the Haldane-Wu 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 mass-squared difference and large mixing angle solution (the so-called 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 non-vanishing masses and mixings among neutrino flavours on the detection of neutrinos from stellar collapse by a water Cerenkov detector. We consider a frame-work 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 electron-antineutrinos changes dramatically for some allowed mixing parameters. Furthermore, contributions from charged-current 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 non-vanishing 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 electron-antineutrinos may change dramatically for some allowed mixing parameters. Furthermore, contributions from charged-current 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 backward-peaked. Hence, mixing has a non-trivial 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 low-energy 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 two-body 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 three-body interaction is included, albeit in a particular model.
  • We report on a study of a finite system of classical confined particles in two-dimensions in the presence of a uniform magnetic field and interacting via a two-body 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 many-body Hamiltonian with two- and three-body interactions in two dimensions introduced recently by Murthy, Bhaduri and Sen. Apart from an analysis of some exact solutions in the many-body system, we analyze in detail the two-body problem which is completely solvable. We show that the solution of the two-body problem reduces to solving a known differential equation due to Heun. We show that the two-body 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 inverse-square potential properly and discuss its implications for our model.
  • We construct a new many-body Hamiltonian with two- and three-body interactions in two space dimensions and obtain its exact many-body 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 Calogero-Sutherland model in one dimension. The model reduces to an analog of the well-known trigonometric Sutherland model when projected on to a circular ring.