Fast radio bursts (FRBs) are bright, millisecond-duration radio transients originating from extragalactic distances. Their origin is unknown. Some FRB sources emit repeat bursts, ruling out cataclysmic origins for those events. Despite searches for periodicity in repeat burst arrival times on time scales from milliseconds to many days, these bursts have hitherto been observed to appear sporadically, and though clustered, without a regular pattern. Here we report the detection of a $16.35\pm0.15$ day periodicity (or possibly a higher-frequency alias of that periodicity) from a repeating FRB 180916.J0158+65 detected by the Canadian Hydrogen Intensity Mapping Experiment Fast Radio Burst Project (CHIME/FRB). In 38 bursts recorded from September 16th, 2018 through February 4th, 2020, we find that all bursts arrive in a 5-day phase window, and 50% of the bursts arrive in a 0.6-day phase window. Our results suggest a mechanism for periodic modulation either of the burst emission itself, or through external amplification or absorption, and disfavour models invoking purely sporadic processes.

Magnetars are highly magnetized young neutron stars that occasionally produce enormous bursts and flares of X-rays and gamma-rays. Of the approximately thirty magnetars currently known in our Galaxy and Magellanic Clouds, five have exhibited transient radio pulsations. Fast radio bursts (FRBs) are millisecond-duration bursts of radio waves arriving from cosmological distances. Some have been seen to repeat. A leading model for repeating FRBs is that they are extragalactic magnetars, powered by their intense magnetic fields. However, a challenge to this model has been that FRBs must have radio luminosities many orders of magnitude larger than those seen from known Galactic magnetars. Here we report the detection of an extremely intense radio burst from the Galactic magnetar SGR 1935+2154 using the Canadian Hydrogen Intensity Mapping Experiment (CHIME) FRB project. The fluence of this two-component bright radio burst and the estimated distance to SGR 1935+2154 together imply a 400-800 MHz burst energy of $\sim 3 \times 10^{34}$ erg, which is three orders of magnitude brighter than those of any radio-emitting magnetar detected thus far. Such a burst coming from a nearby galaxy would be indistinguishable from a typical FRB. This event thus bridges a large fraction of the radio energy gap between the population of Galactic magnetars and FRBs, strongly supporting the notion that magnetars are the origin of at least some FRBs.

In experiments searching for axionic dark matter, the use of the standard threshold-based data analysis discards valuable information. We present a Bayesian analysis framework that builds on an existing processing protocol to extract more information from the data of coherent axion detectors such as operating haloscopes. The analysis avoids logical subtleties that accompany the standard analysis framework and enables greater experimental flexibility on future data runs. Performing this analysis on the existing data from the HAYSTAC experiment, we find improved constraints on the axion-photon coupling $g_\gamma$ while also identifying the most promising regions of parameter space within the $23.15$--$24.0$ $\mu$eV mass range. A comparison with the standard threshold analysis suggests a $36\%$ improvement in scan rate from our analysis, demonstrating the utility of this framework for future axion haloscope analyses.

This paper presents EinsteinPy (version 0.3), a community-developed Python package for gravitational and relativistic astrophysics. Python is a free, easy to use a high-level programming language which has seen a huge expansion in the number of its users and developers in recent years. Specifically, a lot of recent studies show that the use of Python in Astrophysics and general physics has increased exponentially. We aim to provide a very high level of abstraction, an easy to use interface and pleasing user experience. EinsteinPy is developed keeping in mind the state of a theoretical gravitational physicist with little or no background in computer programming and trying to work in the field of numerical relativity or trying to use simulations in their research. Currently, EinsteinPy supports simulation of time-like and null geodesics and calculates trajectories in different background geometries some of which are Schwarzschild, Kerr, and KerrNewmann along with coordinate inter-conversion pipeline. It has a partially developed pipeline for plotting and visualization with dependencies on libraries like Plotly, matplotlib, etc. One of the unique features of EinsteinPy is a sufficiently developed symbolic tensor manipulation utilities which are a great tool in itself for teaching yourself tensor algebra which for many beginner students can be overwhelmingly tricky. EinsteinPy also provides few utility functions for hypersurface embedding of Schwarzschild spacetime which further will be extended to model gravitational lensing simulation.

We introduce DAYENU, a linear, spectral filter for HI intensity mapping that achieves the desirable foreground mitigation and error minimization properties of inverse co-variance weighting with minimal modeling of the underlying data. Beyond 21 cm power-spectrum estimation, our filter is suitable for any analysis where high dynamic-range removal of spectrally smooth foregrounds in irregularly (or regularly) sampled data is required, something required by many other intensity mapping techniques. We show that DAYENU enables the access of large-scale line-of-sight modes that are inaccessible to tapered DFT estimators. Since these modes have the largest SNRs, DAYENU significantly increases the sensitivity of 21 cm analyses over tapered Fourier transforms. Slight modifications allow us to use DAYENU as a linear replacement for iterative delay CLEANing (DAYENUREST). An interactive jupyter tutorial on using DAYENU can be found at https://github.com/HERA-Team/uvtools/blob/master/examples/linear_clean_demo.ipynb DAYENU's source code can be found at https://github.com/HERA-Team/uvtools/blob/master/uvtools/dspec.py

Fluctuation at the neutron star center gives rise to a small deconfined quark core very close to the star center. The density discontinuity at the quark-hadron boundary initiates a shock wave, which propagates outwards of the star. The shock is strong enough to combust nuclear matter to 2-flavor quark matter in the star. The 2-flavor quark matter is not stable and there is an excess of down quarks. The 2-flavor matter settles to a stable 3-flavor matter in the weakly interacting time-scale. In this paper particularly study the conversion of 2-flavor matter to 3-flavor matter is carried out. We set up a differential equation for the conversion of the excess of down quarks to strange quarks involving weak reaction and diffusion of quarks. Calculating the reaction rate we solve the differential equation to find the velocity of the conversion front. The conversion velocity is about $0.002$ times of the speed of light and the time taken for the conversion is about $1.96$ millisecond. Once the front velocity is known we find the change in density profile of the star. As the conversion front moves out of the star the density profile changes bringing about a change in the quadrupole moment of the star. The change in the quadrupole moment is reflected in the GW amplitude which is of the order of $10^{-24}-10^{-23}$. The power spectrum has a peak frequencies of around $10-20$ kHz. Such amplitude and peak frequency is a unique signature of shock-induced first-order phase transition.

We present B-, V-, R-, and I-bands linear polarimetric observations of 73 stars in the direction of open star cluster Casado Alessi 1 (hereafter Alessi 1). We aim to use polarimetry as a tool to investigate the properties and distribution of dust grains toward the direction of the cluster. The polarimetric observations were carried out using the ARIES IMaging POLarimeter mounted at the 104 cm telescope of ARIES, Nainital (India). Using the Gaia photometric data the age and distance of the cluster are estimated to be $0.8\pm0.1$ Gyr and $673\pm98$ pc, respectively. A total of 66 stars with a 26 arcmin radius from the cluster are identified as members of the cluster using the astrometric approach. Out of these 66 members, 15 stars were observed polarimetrically and found to have the same value of polarization. The majority of the stars in the region follow the general law of the polarization for the interstellar medium, indicating that polarization toward the cluster Alessi 1 is dominated by foreground dust grains. The average values of the maximum polarization ($P_{max}$) and the wavelength corresponding to the maximum polarization ($\lambda_{max}$) toward the cluster are found to be $0.83\pm0.03$% and $0.59\pm0.04$ $\mu$m, respectively. Also, dust grains toward the cluster appear to be aligned, possibly due to the galactic magnetic field.

Fundamental plane of elliptical galaxies can be used to predict the intrinsic size of galaxies and has a number of plausible application to study cosmology and galaxy physics. We present a detailed analysis of the fundamental plane of the SDSS-III BOSS LOWZ and CMASS galaxies. For the standard fundamental plane, we find a strong redshift evolution for the mean residual and show that it is primarily driven by the redshift evolution of the surface brightness of the galaxies. After correcting for the redshift evolution, the FP residuals are strongly correlated with the galaxy properties and some observational systematics. We show that the variations in the FP between the central and satellite galaxies, that have been observed in the literature, can primarily be explained by the correlation of the FP with the galaxy luminosity. We also measure the cross correlations of the FP residuals with the galaxy density field. The amplitude of the cross correlations depends on the galaxy properties and environment with brighter and redder galaxies showing stronger correlation. In general, galaxies in denser environments (higher galaxy bias ) show stronger correlations. We also compare FP amplitude with the amplitudes of intrinsic alignments of galaxy shapes (IA), finding the two to be correlated. Finally, using the FP residuals we also study the impact of intrinsic alignments on the constraint of growth rate using redshift space distortions. We do not observe any significant trends in measurements of the growth rate $f$ as function of the amplitude of FP-density correlations, resulting in null detection of the effects of IA on the RSD measurements.

We present cosmological parameter constraints based on a joint modeling of galaxy-lensing cross correlations and galaxy clustering measurements in the SDSS, marginalizing over small-scale modeling uncertainties using mock galaxy catalogs, without explicit modeling of galaxy bias. We show that our modeling method is robust to the impact of different choices for how galaxies occupy dark matter halos and to the impact of baryonic physics (at the $\sim2\%$ level in cosmological parameters) and test for the impact of covariance on the likelihood analysis and of the survey window function on the theory computations. Applying our results to the measurements using galaxy samples from BOSS and lensing measurements using shear from SDSS galaxies and CMB lensing from Planck, with conservative scale cuts, we obtain $S_8\equiv\left(\frac{\sigma_8}{0.8228}\right)^{0.8}\left(\frac{\Omega_m}{0.307}\right)^{0.6}=0.85\pm0.05$ (stat.) using LOWZ $\times$ SDSS galaxy lensing, and $S_8=0.91\pm0.1$ (stat.) using combination of LOWZ and CMASS $\times$ Planck CMB lensing. We estimate the systematic uncertainty in the galaxy-galaxy lensing measurements to be $\sim6\%$ (dominated by photometric redshift uncertainties) and in the galaxy-CMB lensing measurements to be $\sim3\%$, from small scale modeling uncertainties including baryonic physics.

We describe our nonlinear emulation (i.e., interpolation) framework that combines the halo occupation distribution (HOD) galaxy bias model with $N$-body simulations of nonlinear structure formation, designed to accurately predict the projected clustering and galaxy-galaxy lensing signals from luminous red galaxies (LRGs) in the redshift range $0.16 < z < 0.36$ on comoving scales $0.6 < r_p < 30$ \hMpc. The interpolation accuracy is $\lesssim 1-2$ per cent across the entire physically plausible range of parameters for all scales considered. We correctly recover the true value of the cosmological parameter $S_8 = ({\sigma_8}/{0.8228}) ({\Omega_{\text{m}}}/{0.3107})^{0.6}$ from mock measurements produced via subhalo abundance matching (SHAM)-based lightcones designed to approximately match the properties of the SDSS LOWZ galaxy sample. Applying our model to Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 14 (DR14) LOWZ galaxy clustering and galaxy-shear cross-correlation measurements made with Sloan Digital Sky Survey (SDSS) Data Release 8 (DR8) imaging, we perform a prototype cosmological analysis marginalizing over $w$CDM cosmological parameters and galaxy HOD parameters. We obtain a 4.4 per cent measurement of $S_8 = 0.847 \pm 0.037$, in $3.5\sigma$ tension with the Planck cosmological results of $1.00 \pm 0.02$. We discuss the possibility of underestimated systematic uncertainties or astrophysical effects that could explain this discrepancy.

We explore the viability of laboratory-scale mechanical resonators as detectors for ultralight scalar dark matter. The signal we investigate is an atomic strain due to modulation of the fine structure constant and the lepton mass at the Compton frequency of dark matter particles. The resulting stress can drive an elastic body with acoustic breathing modes, producing displacements that are accessible with opto- or electromechanical readout techniques. To address the unknown mass of dark matter particles (which determines their Compton frequency), we consider various resonator designs operating at kHz to MHz frequencies, corresponding to $10^{-12}-10^{-5}$ eV particle mass. Current resonant-mass gravitational wave detectors that have been repurposed as dark matter detectors weigh $\sim \! 10^3$ kg. We find that a large unexplored parameter space can be accessed with ultra-high-$Q$, cryogenically-cooled, cm-scale mechanical resonators possessing $\sim \! 10^7$ times smaller mass.

The splashback radius ($R_{\rm sp}$) of dark matter halos has recently been detected using weak gravitational lensing and cross-correlations with galaxies. However, different methods have been used to measure $R_{\rm sp}$ and to assess the significance of its detection. In this paper, we use simulations to study the precision and accuracy to which we can detect the splashback radius with 3D density, 3D subhalo, and weak lensing profiles. We study how well various methods and tracers recover $R_{\rm sp}$ by comparing it with the value measured directly from particle dynamics. We show that estimates of $R_{\rm sp}$ from density and subhalo profiles correspond to different percentiles of the underlying $R_{\rm sp}$ distribution of particle orbits. At low accretion rates, a second caustic appears and can bias results. Finally, we show that upcoming lensing surveys may be able to constrain the splashback-accretion rate relation directly.

This white paper submitted for the 2020 Canadian Long-Range Planning process (LRP2020) presents the prospects for Canada and the Square Kilometre Array (SKA) from 2020-2030, focussing on the first phase of the project (SKA1) scheduled to begin construction early in the next decade. SKA1 will make transformational advances in our understanding of the Universe across a wide range of fields, and Canadians are poised to play leadership roles in several. Canadian key SKA technologies will ensure a good return on capital investment in addition to strong scientific returns, positioning Canadian astronomy for future opportunities well beyond 2030. We therefore advocate for Canada's continued scientific and technological engagement in the SKA from 2020-2030 through participation in the construction and operations phases of SKA1.

We report on the discovery of eight repeating fast radio burst (FRB) sources found using the Canadian Hydrogen Intensity Mapping Experiment (CHIME) telescope. These sources span a dispersion measure (DM) range of 103.5 to 1281 pc cm$^{-3}$. They display varying degrees of activity: six sources were detected twice, another three times, and one ten times. These eight repeating FRBs likely represent the bright and/or high-rate end of a distribution of infrequently repeating sources. For all sources, we determine sky coordinates with uncertainties of $\sim$10$^\prime$. FRB 180916.J0158+65 has a burst-averaged DM = $349.2 \pm 0.3$ pc cm$^{-3}$ and a low DM excess over the modelled Galactic maximum (as low as $\sim$20 pc cm$^{-3}$); this source also has a Faraday rotation measure (RM) of $-114.6 \pm 0.6$ rad m$^{-2}$, much lower than the RM measured for FRB 121102. FRB 181030.J1054+73 has the lowest DM for a repeater, $103.5 \pm 0.3$ pc cm$^{-3}$, with a DM excess of $\sim$ 70 pc cm$^{-3}$. Both sources are interesting targets for multi-wavelength follow-up due to their apparent proximity. The DM distribution of our repeater sample is statistically indistinguishable from that of the first 12 CHIME/FRB sources that have not repeated. We find, with 4$\sigma$ significance, that repeater bursts are generally wider than those of CHIME/FRB bursts that have not repeated, suggesting different emission mechanisms. Our repeater events show complex morphologies that are reminiscent of the first two discovered repeating FRBs. The repetitive behavior of these sources will enable interferometric localizations and subsequent host galaxy identifications.

Generating mocks for future sky surveys requires large volumes and high resolutions, which is computationally expensive even for fast simulations. In this work we try to develop numerical schemes to calibrate various halo and matter statistics in fast low resolution simulations compared to high resolution N-body and hydrodynamic simulations. For the halos, we improve the initial condition accuracy and develop a halo finder "relaxed-FOF", where we allow different linking length for different halo mass and velocity dispersions. We show that our relaxed-FoF halo finder improves the common statistics, such as halo bias, halo mass function, halo auto power spectrum in real space and in redshift space, cross correlation coefficient with the reference halo catalog, and halo-matter cross power spectrum. We also incorporate the potential gradient descent (PGD) method into fast simulations to improve the matter distribution at nonlinear scale. By building a lightcone output, we show that the PGD method significantly improves the weak lensing convergence tomographic power spectrum. With these improvements FastPM is comparable to the high resolution full N-body simulation of the same mass resolution, with two orders of magnitude fewer time steps. These techniques can be used to improve the halo and matter statistics of FastPM simulations for mock catalogs of future surveys such as DESI and LSST.

The EDGES collaboration reported the finding of an unexpectedly deep absorption in the radio background at 78 MHz and interpreted the dip as a first detection of redshifted 21-cm from Cosmic Dawn. We have attempted an alternate analysis, adopting a maximally smooth function approach to model the foreground. A joint fit to the spectrum using such a function together with a flattened absorption profile yields a best fit absorption amplitude of $921 \pm 35$ mK. The depth of the 21-cm absorption inferred by the EDGES analysis required invoking non-standard cosmology or new physics or new sources at Cosmic Dawn and this tension with accepted models is compounded by our analysis that suggests absorption of greater depth. Alternatively, the measured spectrum may be equally-well fit assuming that there exists a residual unmodeled systematic sinusoidal feature and we explore this possibility further by examining for any additional 21-cm signal. The data then favors an absorption with Gaussian model parameters of amplitude $133 \pm 60$ mK, best width at half-power $9 \pm 3$ MHz and center frequency $72.5 \pm 0.8$ MHz. We also examine the consistency of the measured spectrum with plausible redshifted 21-cm models: a set of 3 of the 265 profiles in the global 21-cm atlas of Cohen et al. 2017 are favored by the spectrum. We conclude that the EDGES data may be consistent with standard cosmology and astrophysics, without invoking excess radio backgrounds or baryon-dark matter interactions.

The advent of new deep+wide photometric lensing surveys will open up the possibility of direct measurements of the dark matter halos of dwarf galaxies. The HSC wide survey will be the first with the statistical capability of measuring the lensing signal with high signal-to-noise at log(M*)=8. At this same mass scale, LSST will have the most overall constraining power with a predicted signal-to-noise for the galaxy-galaxy lensing signal around dwarfs of SN=200. WFIRST and LSST will have the greatest potential to push below the log(M*) = 7 mass scale thanks to the depth of their imaging data. Studies of the dark matter halos of dwarf galaxies at z=0.1 with gravitational lensing are soon within reach. However, further work will be required to develop optimized strategies for extracting dwarfs samples from these surveys, determining redshifts, and accurately measuring lensing on small radial scales. Dwarf lensing will be a new and powerful tool to constrain the halo masses and inner density slopes of dwarf galaxies and to distinguish between baryonic feedback and modified dark matter scenarios.

We study the significance of non-Gaussianity in the likelihood of weak lensing shear two-point correlation functions, detecting significantly non-zero skewness and kurtosis in one-dimensional marginal distributions of shear two-point correlation functions in simulated weak lensing data though the full multivariate distributions are relatively more Gaussian. We examine the implications in the context of future surveys, in particular LSST, with derivations of how the non-Gaussianity scales with survey area. We show that there is no significant bias in one-dimensional posteriors of $\Omega_{\rm m}$ and $\sigma_{\rm 8}$ due to the non-Gaussian likelihood distributions of shear correlations functions using the mock data ($100$ deg$^{2}$). We also present a systematic approach to constructing an approximate multivariate likelihood function by decorrelating the data points using principal component analysis (PCA). When using a subset of the PCA components that account for the majority of the cosmological signal as a data vector, the one-dimensional marginal likelihood distributions of those components exhibit less skewness and kurtosis than the original shear correlation functions. We further demonstrate that the difference in cosmological parameter constraints between the multivariate Gaussian likelihood model and more complex non-Gaussian likelihood models would be even smaller for an LSST-like survey due to the area effect. In addition, the PCA approach automatically serves as a data compression method, enabling the retention of the majority of the cosmological information while reducing the dimensionality of the data vector by a factor of $\sim$5.

The Core Cosmology Library (CCL) provides routines to compute basic cosmological observables to a high degree of accuracy, which have been verified with an extensive suite of validation tests. Predictions are provided for many cosmological quantities, including distances, angular power spectra, correlation functions, halo bias and the halo mass function through state-of-the-art modeling prescriptions available in the literature. Fiducial specifications for the expected galaxy distributions for the Large Synoptic Survey Telescope (LSST) are also included, together with the capability of computing redshift distributions for a user-defined photometric redshift model. A rigorous validation procedure, based on comparisons between CCL and independent software packages, allows us to establish a well-defined numerical accuracy for each predicted quantity. As a result, predictions for correlation functions of galaxy clustering, galaxy-galaxy lensing and cosmic shear are demonstrated to be within a fraction of the expected statistical uncertainty of the observables for the models and in the range of scales of interest to LSST. CCL is an open source software package written in C, with a python interface and publicly available at https://github.com/LSSTDESC/CCL.

Community access to deep (i ~ 25), highly-multiplexed optical and near-infrared multi-object spectroscopy (MOS) on 8-40m telescopes would greatly improve measurements of cosmological parameters from LSST. The largest gain would come from improvements to LSST photometric redshifts, which are employed directly or indirectly for every major LSST cosmological probe; deep spectroscopic datasets will enable reduced uncertainties in the redshifts of individual objects via optimized training. Such spectroscopy will also determine the relationship of galaxy SEDs to their environments, key observables for studies of galaxy evolution. The resulting data will also constrain the impact of blending on photo-z's. Focused spectroscopic campaigns can also improve weak lensing cosmology by constraining the intrinsic alignments between the orientations of galaxies. Galaxy cluster studies can be enhanced by measuring motions of galaxies in and around clusters and by testing photo-z performance in regions of high density. Photometric redshift and intrinsic alignment studies are best-suited to instruments on large-aperture telescopes with wider fields of view (e.g., Subaru/PFS, MSE, or GMT/MANIFEST) but cluster investigations can be pursued with smaller-field instruments (e.g., Gemini/GMOS, Keck/DEIMOS, or TMT/WFOS), so deep MOS work can be distributed amongst a variety of telescopes. However, community access to large amounts of nights for surveys will still be needed to accomplish this work. In two companion white papers we present gains from shallower, wide-area MOS and from single-target imaging and spectroscopy.

LSST will open new vistas for cosmology in the next decade, but it cannot reach its full potential without data from other telescopes. Cosmological constraints can be greatly enhanced using wide-field ($>20$ deg$^2$ total survey area), highly-multiplexed optical and near-infrared multi-object spectroscopy (MOS) on 4-15m telescopes. This could come in the form of suitably-designed large surveys and/or community access to add new targets to existing projects. First, photometric redshifts can be calibrated with high precision using cross-correlations of photometric samples against spectroscopic samples at $0 < z < 3$ that span thousands of sq. deg. Cross-correlations of faint LSST objects and lensing maps with these spectroscopic samples can also improve weak lensing cosmology by constraining intrinsic alignment systematics, and will also provide new tests of modified gravity theories. Large samples of LSST strong lens systems and supernovae can be studied most efficiently by piggybacking on spectroscopic surveys covering as much of the LSST extragalactic footprint as possible (up to $\sim20,000$ square degrees). Finally, redshifts can be measured efficiently for a high fraction of the supernovae in the LSST Deep Drilling Fields (DDFs) by targeting their hosts with wide-field spectrographs. Targeting distant galaxies, supernovae, and strong lens systems over wide areas in extended surveys with (e.g.) DESI or MSE in the northern portion of the LSST footprint or 4MOST in the south could realize many of these gains; DESI, 4MOST, Subaru/PFS, or MSE would all be well-suited for DDF surveys. The most efficient solution would be a new wide-field, highly-multiplexed spectroscopic instrument in the southern hemisphere with $>6$m aperture. In two companion white papers we present gains from deep, small-area MOS and from single-target imaging and spectroscopy.

Previous studies have shown the filamentary structures in the cosmic web influence the alignments of nearby galaxies. We study this effect in the LOWZ sample of the Sloan Digital Sky Survey using the "Cosmic Web Reconstruction" filament catalogue. We find that LOWZ galaxies exhibit a small but statistically significant alignment in the direction parallel to the orientation of nearby filaments. This effect is detectable even in the absence of nearby galaxy clusters, which suggests it is an effect from the matter distribution in the filament. A nonparametric regression model suggests that the alignment effect with filaments extends over separations of 30-40 Mpc. We find that galaxies that are bright and early-forming align more strongly with the directions of nearby filaments than those that are faint and late-forming; however, trends with stellar mass are less statistically significant, within the narrow range of stellar mass of this sample.

Optimal analyses using the 2-point functions of large-scale structure probes require accurate covariance matrices. A covariance matrix of the 2-point function comprises the disconnected part and the connected part. While the connected covariance only becomes important on small scales, the disconnected covariance is dominant on large scales, where the survey window has a significant impact. In this work, we develop an analytical method to compute the disconnected covariance, accounting for the window effect. Derived under the flat-sky approximation, our formalism is applicable to wide surveys by swapping in the curved-sky window functions. Our method works for both the power spectrum and the correlation function, and applies to the covariances of various probes including the multipoles and the wedges of 3D clustering, the angular and the projected statistics of clustering and shear, as well as the cross covariances between different probes. We verify the analytic covariance against the sample covariance from the galaxy mock simulations in two test cases: (1) the power spectrum multipole covariance, and (2) the joint covariance of the projected correlation function and the correlation function multipoles. Our method achieve good agreement with the mocks, while at a negligible computational cost. Unlike mocks, our analytic covariance is free of sampling noise, which often leads to numerical problems and the need to inflate the errors. In addition, our method can use the best-fit power spectrum as input, in contrast to the standard procedure of using a fiducial model that may deviate significantly from the truth. We also show that a naive diagonal power spectrum covariance underestimates the signal-to-noise ratio compared to our analytic covariance. The code that accompanies this paper is available at https://github.com/eelregit/covdisc.

We present measurements of $E_G$, a probe of gravity from large-scale structure, using BOSS LOWZ and CMASS spectroscopic samples, with lensing measurements from SDSS (galaxy lensing) and Planck (CMB lensing). Using SDSS lensing and the BOSS LOWZ sample, we measure $\langle{E_G}\rangle=0.40^{+0.05}_{-0.04}$ (stat), $\pm 0.026$ (systematic), consistent with the predicted value from the Planck $\Lambda$CDM model, $E_G=0.46$. Using CMB lensing, we measure $\langle{E_G}\rangle=0.46^{+0.08}_{-0.09}$ (stat) for LOWZ (statistically consistent with galaxy lensing and Planck predictions) and $\langle{E_G}\rangle=0.39^{+0.05}_{-0.05}$ (stat) for the CMASS sample, consistent with the Planck prediction of $E_G=0.40$ given the higher redshift of the sample. We also study the redshift evolution of $E_G$ by splitting the LOWZ sample into two samples based on redshift, with results being consistent with model predictions. We estimate systematic uncertainties on the above $\langle{E_G}\rangle$ numbers to be $\sim 6$% (when using galaxy-galaxy lensing) or $\sim 3$% (when using CMB lensing), subdominant to the quoted statistical errors. These systematic error budgets are dominated by observational systematics in galaxy-galaxy lensing and by theoretical modeling uncertainties, respectively. We do not estimate observational systematics in galaxy-CMB lensing cross correlations.

We propose and study a semi-classical cosmological system akin to the Newton-Schr\"odinger equation where matter field evolution is determined by time dependent Schr\"odinger equation. The resulting dynamics is one where the scale factor self-consistently informs the quantum evolution of the scalar field wave function via the Friedmann equation. We study various potentials, and show that for each case, inflation arises naturally, without fine-tuning, for arbitrary initial wave functions for the scalar field. Furthermore, due to zero-point quantum fluctuations of the scalar field, the equation of state $P=-\rho$ arises dynamically at late times for all initial states and potentials.

In this work, we have studied the effect of higher order perturbations, particularly the second order in details, on the sonic horizon. We have considered two different schemes of perturbations which are velocity potential perturbation and mass acceleration rate perturbation. These two schemes give us qualitatively similar behaviour. We have found that the analogue gravity formalism also holds for the higher order perturbations.

Spectral distortions in the cosmic microwave background over the 40--200~MHz band are imprinted by neutral hydrogen in the intergalactic medium prior to the end of reionization. This signal, produced in the redshift range $z = 6-34$ at the rest frame wavelength of 21 cm, has not been detected yet; and poor understanding of high redshift astrophysics results in a large uncertainty in the expected spectrum. The SARAS~2 radiometer was purposely designed to detect the sky-averaged 21-cm signal. The instrument, deployed at the Timbaktu Collective (Southern India) in April--June 2017, collected 63~hr of science data, which were examined for the presence of the cosmological 21-cm signal. In our previous work the first-light data from SARAS~2 radiometer were analyzed with Bayesian likelihood-ratio tests using $264$ plausible astrophysical scenarios. In this paper we re-examine the data using an improved analysis based on the frequentist approach and forward modeling. We show that SARAS~2 data rejects 27 models, out of which 25 are rejected at a significance $>5\sigma$. All the rejected models share the scenario of inefficient heating of the primordial gas by the first population of X-ray sources along with rapid reionization.

The global 21 cm signal from Cosmic Dawn (CD) and the Epoch of Reionization (EoR), at redshifts $z \sim 6-30$, probes the nature of first sources of radiation as well as physics of the Inter-Galactic Medium (IGM). Given that the signal is predicted to be extremely weak, of wide fractional bandwidth, and lies in a frequency range that is dominated by Galactic and Extragalactic foregrounds as well as Radio Frequency Interference, detection of the signal is a daunting task. Critical to the experiment is the manner in which the sky signal is represented through the instrument. It is of utmost importance to design a system whose spectral bandpass and additive spurious can be well calibrated and any calibration residual does not mimic the signal. SARAS is an ongoing experiment that aims to detect the global 21 cm signal. Here we present the design philosophy of the SARAS 2 system and discuss its performance and limitations based on laboratory and field measurements. Laboratory tests with the antenna replaced with a variety of terminations, including a network model for the antenna impedance, show that the gain calibration and modeling of internal additives leave no residuals with Fourier amplitudes exceeding 2~mK, or residual Gaussians of 25 MHz width with amplitudes exceeding 2~mK. Thus, even accounting for reflection and radiation efficiency losses in the antenna, the SARAS~2 system is capable of detection of complex 21-cm profiles at the level predicted by currently favoured models for thermal baryon evolution.

We investigate particle production \`a la Schwinger mechanism in an expanding, flat de Sitter patch as is relevant for the inflationary epoch of our universe. Defining states and particle content in curved spacetime is certainly not a unique process. There being different prescriptions on how that can be done, we have used the Schr\"odinger formalism to define instantaneous particle content of the state etc. This allows us to go past the adiabatic regime to which the effect has been restricted in the previous studies and bring out its multifaceted nature in different settings. Each of these settings gives rise to contrasting features and behaviour as per the effect of electric field and expansion rate on the instantaneous mean particle number. We also quantify the degree of classicality of the process during its evolution using a "classicality parameter" constructed out of parameters of the Wigner function to obtain information about the quantum to classical transition in this case.

We study the covariance properties of real space correlation function estimators -- primarily galaxy-shear correlations, or galaxy-galaxy lensing -- using SDSS data for both shear catalogs and lenses (specifically the BOSS LOWZ sample). Using mock catalogs of lenses and sources, we disentangle the various contributions to the covariance matrix and compare them with a simple analytical model. We show that not subtracting the lensing measurement around random points from the measurement around the lens sample is equivalent to performing the measurement using the lens density field instead of the lens over-density field. While the measurement using the lens density field is unbiased (in the absence of systematics), its error is significantly larger due to an additional term in the covariance. Therefore, this subtraction should be performed regardless of its beneficial effects on systematics. Comparing the error estimates from data and mocks for estimators that involve the over-density, we find that the errors are dominated by the shape noise and lens clustering, that empirically estimated covariances (jackknife and standard deviation across mocks) are consistent with theoretical estimates, and that both the connected parts of the 4-point function and the super-sample covariance can be neglected for the current levels of noise. While the trade-off between different terms in the covariance depends on the survey configuration (area, source number density), the diagnostics that we use in this work should be useful for future works to test their empirically-determined covariances.

Long wavelength spectral distortions in the Cosmic Microwave Background arising from the 21-cm transition in neutral Hydrogen are a key probe of Cosmic Dawn and the Epoch of Reionization. These features may reveal the nature of the first stars and ultra-faint galaxies that transformed the spin temperature and ionization state of the primordial gas. SARAS~2 is a spectral radiometer purposely designed for precision measurement of these monopole or all-sky global 21-cm spectral distortions. We use 63~hr night time observing of the radio background in the frequency band 110-200~MHz with the radiometer deployed at the Timbaktu Collective in Southern India to derive likelihoods for plausible redshifted 21-cm signals predicted by theoretical models. First light with SARAS 2 disfavors the class of models that feature weak X-ray heating (with $f_X \leq 0.1$) and rapid reionization (with peak $\frac{dT_b}{dz} \geq 120~\textrm{mK per unit redshift interval}$ ).

We present results from cross-correlating Planck CMB lensing maps with the Sloan Digital Sky Survey (SDSS) galaxy lensing shape catalog and BOSS galaxy catalogs. For galaxy position vs. CMB lensing cross-correlations, we measure the convergence signal around the galaxies in configuration space, using the BOSS LOWZ ($z\sim0.30$) and CMASS ($z\sim0.57$) samples. With fixed Planck 2015 cosmology, doing a joint fit with the galaxy clustering measurement, for the LOWZ (CMASS) sample we find a galaxy bias $b_g=1.75\pm0.04$ ($1.95\pm 0.02$) and galaxy-matter cross-correlation coefficient $r_{cc}=1.0\pm0.2$ ($0.8\pm 0.1$) using $20<r_p<70h^{-1}$Mpc, consistent with results from galaxy-galaxy lensing. Using the same scales and including the galaxy-galaxy lensing measurements, we constrain $\Omega_m=0.284\pm0.024$ and relative calibration bias between the CMB lensing and galaxy lensing to be $b_\gamma=0.82^{+0.15}_{-0.14}$. The combination of galaxy lensing and CMB lensing also allows us to measure the cosmological distance ratios (with $z_l\sim0.3$, $z_s\sim0.5$) $\mathcal R=\frac{D_s D_{l,*}}{D_* D_{l,s}}=2.68\pm0.29$, consistent with predictions from the Planck 2015 cosmology ($\mathcal R=2.35$). We detect the galaxy position-CMB convergence cross-correlation at small scales, $r_p<1h^{-1}$Mpc, and find consistency with lensing by NFW halos of mass $M_h\sim10^{13}h^{-1}M_\odot$. Finally, we measure the CMB lensing-galaxy shear cross-correlation, finding an amplitude of $A=0.76\pm0.23$ ($z_\text{eff}=0.35$, $\theta<2^\circ$) with respect to Planck 2015 $\Lambda$CDM predictions ($1\sigma$-level consistency). We do not find evidence for relative systematics between the CMB and SDSS galaxy lensing.

We establish the existence and find the necessary and sufficient conditions for a new class of solutions of Locally Rotationally Symmetric spacetimes that have non vanishing rotation and spatial twist simultaneously. We transparently show that the existence of such solutions demand non vanishing and bounded heat flux and these solutions are self similar. We provide a brief algorithm indicating how to solve the system of field equations with the given Cauchy data on an initial spacelike Cauchy surface. Finally we argue that these solutions can be used as a first approximation from spherical symmetry to study rotating, inhomogeneous, dynamic and radiating astrophysical stars.

Inflation has by-far set itself as one of the prime ideas in the current cosmological models that seemingly has an answer for every observed phenomenon in cosmology. More importantly, it serves as a bridge between the early quantum fluctuations and the present-day classical structures. Although the transition from quantum to classical is still not completely understood till date, there are two assumptions made in the inflationary paradigm in this regard: (i) the modes (metric perturbations or fluctuations) behave classically once they are well outside the Hubble radius and, (ii) once they become classical they stay classical and hence can be described by standard perturbation theory after they re-enter the Hubble radius. We critically examine these assumptions for the tensor modes of (linear) metric perturbations in a toy three stage universe with (i) inflation, (ii) radiation-dominated and (iii) late-time accelerated phases. The quantum-to-classical transition for these modes is evident from the evolution of Wigner function in phase space and its peaking on the classical trajectory. However, a better approach to quantify the degree of classicality and study its evolution was given by Mahajan and Padmanabhan [1] using a classicality parameter constructed from the parameters of the Wigner function. We study the evolution of the classicality parameter across the three phases and it turns out that the first assumption holds true, there is emergence of classicality on Hubble exit, however the latter assumption of "once classical, always classical" seems to lie on a shaky ground.

Direct detection of gravitational waves is opening a new window onto our universe. Here, we study the sensitivity to continuous-wave strain fields of a kg-scale optomechanical system formed by the acoustic motion of superfluid helium-4 parametrically coupled to a superconducting microwave cavity. This narrowband detection scheme can operate at very high $Q$-factors, while the resonant frequency is tunable through pressurization of the helium in the 0.1-1.5 kHz range. The detector can therefore be tuned to a variety of astrophysical sources and can remain sensitive to a particular source over a long period of time. For reasonable experimental parameters, we find that strain fields on the order of $h\sim 10^{-23} /\sqrt{\rm Hz}$ are detectable. We show that the proposed system can significantly improve the limits on gravitational wave strain from nearby pulsars within a few months of integration time.

We study the background cosmological dynamics with a three component source content: a radiation fluid, a barotropic fluid to mimic the matter sector and a single scalar field which can act as dark energy giving rise to the late-time accelerated phase. Using the well-known dimensionless variables, we cast the dynamical equations into an autonomous system of ordinary differential equations (ASODE), which are studied by computing the fixed points and the conditions for their stability. The matter fluid and the scalar field are taken to be uncoupled at first and later, we consider a coupling between the two of the form $Q = \sqrt{2/3}\kappa\beta\rho_m\dot{\phi}$ where $\rho_m$ is the barotropic fluid density. The key point of our analysis is that for the closure of ASODE, we only demand that the jerk, $\Gamma = V V"/V'^2$ is a function of acceleration, $z = - M_p V'/ V$, that is, $\Gamma = 1+ f(z)$. In this way, we are able to accommodate a large class of potentials that goes beyond the simple exponential potentials. The analysis is completely generic and \emph{independent} of the form of the potential for the scalar field. As an illustration and confirmation of the analysis, we consider $f(z)$ of the forms $\mu/z^2$, $\mu/z$, $(\mu-z)/z^2$ and $(\mu-z)$ to numerically compute the evolution of cosmological parameters with and without coupling. Implications of the approach and the results are discussed.

Measurements of intrinsic alignments of galaxy shapes with the large-scale density field, and the inferred intrinsic alignments model parameters, are sensitive to the shape measurement methods used. In this paper we measure the intrinsic alignments of the Sloan Digital Sky Survey-III (SDSS-III) Baryon Oscillation Spectroscopic Survey (BOSS) LOWZ galaxies using three different shape measurement methods (re-Gaussianization, isophotal, and de Vaucouleurs), identifying a variation in the inferred intrinsic alignments amplitude at the 40% level between these methods, independent of the galaxy luminosity or other properties. We also carry out a suite of systematics tests on the shapes and their two-point correlation functions, identifying a pronounced contribution from additive PSF systematics in the de Vaucouleurs shapes. Since different methods measure galaxy shapes at different effective radii, the trends we identify in the intrinsic alignments amplitude are consistent with the interpretation that the outer regions of galaxy shapes are more responsive to tidal fields, resulting in isophote twisting and stronger alignments for isophotal shapes. We observe environment dependence of ellipticity, with brightest galaxies in groups being rounder on average compared to satellite and field galaxies. We also study the anisotropy in intrinsic alignments measurements introduced by projected shapes, finding effects consistent with predictions of the nonlinear alignment model and hydrodynamic simulations. The large variations seen using the different shape measurement methods have important implications for intrinsic alignments forecasting and mitigation with future surveys.

The fundamental plane (FP) is a widely used tool to investigate the properties of early-type galaxies, and the tight relation between its parameters has spawned several cosmological applications, including its use as a distance indicator for peculiar velocity surveys and as a means to suppress intrinsic noise in cosmic size magnification measurements. Systematic trends with the large-scale structure across the FP could cause serious biases for these cosmological probes, but may also yield new insights into the early-type population. Here we report the first detection of spatial correlations among offsets in galaxy size from an FP that explicitly accounts for redshift trends, using a sample of about $95,000$ elliptical galaxies from the Sloan Digital Sky Survey. We show that these offsets correlate with the density field out to at least $10h^{-1}$Mpc at $4\sigma$ significance in a way that cannot be explained by systematic errors in galaxy size estimates. We propose a physical explanation for the correlations by dividing the sample into central, satellite, and field galaxies, identifying trends for each galaxy type separately. Central (satellite) galaxies lie on average above (below) the FP, which we argue could be due to a higher (lower) than average mass-to-light ratio. We fit a simple model to the correlations of FP residuals and use it to predict the impact on peculiar velocity power spectra, finding a contamination larger than $10\,\%$ for $k>0.04\,h/$Mpc. Moreover, cosmic magnification measurements based on an FP could be severely contaminated over a wide range of scales by the intrinsic FP correlations.

Detection of the global redshifted 21-cm signal is an excellent means of deciphering the physical processes during the Dark Ages and subsequent Epoch of Reionization (EoR). However, detection of this faint monopole is challenging due to high precision required in instrumental calibration and modeling of substantially brighter foregrounds and instrumental systematics. In particular, modeling of receiver noise with mK accuracy and its separation remains a formidable task in experiments aiming to detect the global signal using single-element spectral radiometers. Interferometers do not respond to receiver noise; therefore, we explore here the theory of the response of interferometers to global signals. In other words, we discuss the spatial coherence in the electric field arising from the monopole component of the 21-cm signal and methods for its detection using sensor arrays. We proceed by first deriving the response to uniform sky of two-element interferometers made of unit dipole and resonant loop antennas, then extend the analysis to interferometers made of 1-D arrays and also consider 2-D aperture antennas. Finally, we describe methods by which the coherence might be enhanced so that the interferometer measurements yield improved sensitivity to the monopole component. We conclude that (a) it is indeed possible to measure the global 21-cm from EoR using interferometers, (b) a practically useful configuration is with omnidirectional antennas as the interferometer elements, and (c) that the spatial coherence may be enhanced and detectability of the global EoR signal may be smoothened using, for example, a space beam splitter between the interferometer elements.

Intrinsic alignments (IA) of galaxies, i.e. correlations of galaxy shapes with each other or with the density field, are a major astrophysical source of contamination for weak lensing surveys. We present the results of IA measurements of galaxies on 0.1- 200 Mpc/h scales using the SDSS-III BOSS LOWZ sample, in the redshift range 0.16<z<0.36. We extend the existing IA measurements for spectroscopic LRGs to lower luminosities, and show that the luminosity dependence of large-scale IA can be well-described by a power law. Within the limited redshift and color range of our sample, we observe no significant redshift or color dependence of IA. We measure the halo mass of LOWZ galaxies using galaxy-galaxy lensing, and show that the mass dependence of large-scale IA is also well described by a power law. We detect variation in the scale dependence of IA with mass and luminosity, which underscores the need to use flexible templates in order to remove the IA signal. We also study the environment dependence of IA by splitting the sample into field and group galaxies, which are further split into satellite and central galaxies. We show that group central galaxies are aligned with their halos at small scales and also are aligned with the tidal fields out to large scales. We also detect the radial alignments of satellite galaxies within groups, which results in a null detection of large-scale intrinsic alignments for satellites. These results can be used to construct better intrinsic alignment models for removal of this contaminant to the weak lensing signal.

The intrinsic alignment of galaxies with the large-scale density field is an important astrophysical contaminant in upcoming weak lensing surveys. We present detailed measurements of the galaxy intrinsic alignments and associated ellipticity-direction (ED) and projected shape ($w_{g+}$) correlation functions for galaxies in the cosmological hydrodynamic MassiveBlack-II (MB-II) simulation. We carefully assess the effects on galaxy shapes, misalignment of the stellar component with the dark matter shape and two-point statistics of iterative weighted (by mass and luminosity) definitions of the (reduced and unreduced) inertia tensor. We find that iterative procedures must be adopted for a reliable measurement of the reduced tensor but that luminosity versus mass weighting has only negligible effects. Both ED and $w_{g+}$ correlations increase in amplitude with subhalo mass (in the range of $10^{10} - 6.0\times 10^{14}h^{-1}M_{\odot}$), with a weak redshift dependence (from $z=1$ to $z=0.06$) at fixed mass. At $z \sim 0.3$, we predict a $w_{g+}$ that is in reasonable agreement with SDSS LRG measurements and that decreases in amplitude by a factor of $\sim 5$--18 for galaxies in the LSST survey. We also compared the intrinsic alignments of centrals and satellites, with clear detection of satellite radial alignments within their host halos. Finally, we show that $w_{g+}$ (using subhalos as tracers of density) and $w_{\delta+}$ (using dark matter density) predictions from the simulations agree with that of non-linear alignment models (NLA) at scales where the 2-halo term dominates in the correlations (and tabulate associated NLA fitting parameters). The 1-halo term induces a scale dependent bias at small scales which is not modeled in the NLA model.

Using the nuclear equation of states for a large variety of relativistic and non-relativistic force parameters, we calculate the static and rotating masses and radii of neutron stars. From these equation of states, we also evaluate the properties of rotating neutron stars, such as rotational and gravitational frequencies, moment of inertia, quadrupole deformation parameter, rotational ellipcity and gravitational wave strain amplitude. The estimated gravitational wave strain amplitude of the star is found to be $\sim 10^{-23}$.

We study the evolution of a quantum scalar field in a toy universe which has three stages of evolution, viz., (i) an early (inflationary) de Sitter phase (ii) radiation dominated phase and (iii) late-time (cosmological constant dominated) de Sitter phase. Using Schr\"odinger picture, the scalar field equations are solved separately for the three stages and matched at the transition points. The boundary conditions are chosen so that field modes in the early de Sitter evolves from Bunch-Davies vacuum state. We determine the (time-dependent) particle content of this quantum state for the entire evolution of the universe and describe the various features both numerically and analytically. We also describe the quantum to classical transition in terms of a {\it classicality parameter} which tracks the particle creation and its effect on phase space correlation of the quantum field.

We investigate the anomalies in the Earth - Moon system using ancient eclipse data. We identify nine groups of anomalous eclipses between 400 and 1800 AD recorded in parts of India that should have completely missed the subcontinent as per NASA simulations (Espenak and Meeus, 2011). We show that the typical correction to the lunar location required to reconcile the anomalous eclipses is relatively small and consistent with the fluctuations in the length of day that are observed in recent periods. We then investigate the change in Earth's moment of inertia due to differential acceleration of land and water that can account for this discrepancy. We show that 80 percent of these discrepancies occur when the Moon is at declinations greater than 10 deg and closer to its major standstill of 28 deg while it spends 46 percent of the time in this region. We simulate the differential interaction of the Moon's gravity with landmass and water using finite element method to account for landmass and water mass. We show that the results of eclipse error are consistent with the estimate of a small differential acceleration when the Moon is over land at high latitudes. However, we encounter some examples where the results from simulations studies cannot explain the phenomena. Hence we propose that dT corrections have to be coupled with some other mechanism possibly a small vertical oscillation in the Moon's rotational plane with a period of the order of a few hundred years to achieve the required adjustment in the Eclipse maps.

The U(1) gauge dynamics on a D4-brane is revisited, with a two form, to construct an effective curvature theory in a second order formalism. We exploit the local degrees in a two form, and modify its dynamics in a gauge invariant way, to incorporate a non-perturbative metric fluctuation in an effective D4-brane. Interestingly, the near horizon D4-brane is shown to describe an asymptotic Anti de Sitter (AdS) in a semi-classical regime. Using Weyl scaling(s), we obtain the emergent rotating geometries leading to primordial de Sitter (dS) and AdS vacua in a quantum regime. Under a discrete transformation, we re-arrange the mixed dS patches to describe a Schwazschild-like dS (SdS) and a topological-like dS (TdS) black holes. We analyze SdS vacuum for Hawking radiations to arrive at Nariai geometry, where a discrete torsion forms a condensate. We perform thermal analysis to identify Nariai vacuum with a TdS. Investigation reveals an AdS patch within a thermal dS brane, which may provide a clue to unfold dS/CFT. In addition, the role of dark energy, sourced by a discrete torsion, in the dS vacua is investigated using Painleve geometries. It is argued that a D-instanton pair is created by a discrete torsion, with a Big Bang/Crunch, at the past horizon in a pure dS. Nucleation, of brane/anti-brane pair(s), is qualitatively analyzed to construct an effective space-time on a D4-brane and its anti brane. Analysis re-assures the significant role played by a non-zero mode, of NS-NS two form, to generalize the notion of branes within a brane.

ARIES Baker-Nunn Schmidt telescope project is converting a Baker-Nunn satellite tracking camera for Astronomical research. Original Baker-Nunn camera produces an extremely large (5X30 degree) curved focal plane at the prime focus for photographic imaging. We present here the re-designing of the camera produces a wide (4 X 4 degree) flat field of view for CCD imaging observations, which have many scientific potentials in Astronomy. Imaging performance of the CCD camera is also estimated.

We constrain the parameters of the variable Chaplygin gas model, using the location of peaks of the CMBR spectrum and SNe Ia ``gold '' data set. Equation of state of the model is $P=-A(a)/\rho$, where $A(a)=A_0 a^{-n}$ is a positive function of the cosmological scale factor $a$, $A_0$ and $n$ being constants. The variable Chaplygin gas interpolates from dust-dominated era to quintessence dominated era. The model is found to be compatible with current type Ia Supernovae data and location of first peak if the values of $\Omega_m$ and $n$ lie in the interval $[0.017,~0.117]$ and $[-1.3,~2.6]$ respectively.

We use the two-dimensional Kolmogorov-Smirnov (KS) test to study the correlation between the 60 cosmic ray events above 4x10^19 eV from the AGASA experiment and the positions of infrared luminous galaxies from the IRAS PSCz catalog. These galaxies are expected to be hosts to gamma ray bursts (GRB) and magnetars, both of which are associated with core collapse supernovae and have been proposed as possible acceleration sites for ultra high energy cosmic rays. We find consistency between the models and the AGASA events to have been drawn from the same underlying distribution of positions on the sky with KS probabilities ~50%. Application of the same test to the 11 highest AGASA events above 10^20 eV, however, yields a KS probability of < 0.5%, rejecting the models at >99.5% significance level. Taken at face value, these highest energy results suggest that the existing cosmic ray events above 10^20 eV do not owe their origin to long burst GRBs, rapidly rotating magnetars, or any other events associated with core collapse supernovae. The larger data set expected from the AUGER experiment will test whether this conclusion is real or is a statistical fluke that we estimate to be at the 2 sigma level.

We develop a method based on the collisionless Boltzmann equation to calculate the gravitational clustering of relic neutrinos in realistic cosmological models dominated by cold dark matter (CDM) and the cosmological constant. This method can be used to estimate the phase-space distribution of any light particles in CDM halos. We find that neutrinos with masses > 0.3 eV cluster appreciably in dark matter halos above the galactic size. The resulting neutrino overdensity above the cosmic mean neutrino density increases with both the neutrino mass and the halo mass, ranging from about 10 for 0.3 eV neutrinos in 10^13 M_sun halos, to about 1500 for 1.8 eV neutrinos in 10^15 M_sun halos. We examine the implications of neutrino clustering for the Z-burst model of ultra high energy cosmic rays (UHECR), which interprets the observed events at E > 4 x 10^19 eV as decay products of Z-bosons from the resonant scattering between relic and high energy neutrinos and anti-neutrinos. We estimate the UHECR energy spectrum for various neutrino masses towards five of the most massive clusters in the local universe (within 100 Mpc): Virgo, Perseus-Pisces, Hydra, Centaurus, and Coma. The UHECR flux in the Z-burst model is expected to be significantly higher towards these clusters if m_nu > 0.3 eV and nearly isotropic otherwise.

We study the effect of inclusion of muons and the muon neutrinos on the phase transition from nuclear to quark matter in a magnetised proto-neutron star and compare our results with those obtained by us without the muons. We find that the inclusion of muons changes slightly the nuclear density at which transition occurs.However the dependence of this transition density on various chemical potentials, temperature and the magnetic field remains quantitatively the same.

We report on the study of the mass-radius (M-R) relation and the radial oscillations of proto strange stars. For the quark matter we have employed the well known density dependent quark mass model and its very recent modification, the temperature and density dependent quark mass model. We find that the maximum mass the star can support increases significantly with the temperature of the star in this model which implies that transition to a black hole at the early stage of formation of the star is inhibited. As for the neutrinos, we find, contrary to the expectation that the M-R and oscillation frequencies are almost independent of the neutrino chemical potentials.

Studies of the growth of cosmic perturbations are typically focused on galactic scales and above. In this paper we investigate the evolution of perturbations in baryons, photons, and dark matter for masses below 10^6 M_\odot (or wavenumbers above 100 Mpc^-1). Fluctuations on these scales are of interest and importance because they grow to become the earliest collapsed objects and provide the first light sources in the so called dark ages. We investigate both the linear evolution and the second-order nonlinear effects arising from the coupling of large-scale velocity fields to small-scale perturbations in the baryon density and the electron ionization fraction. We find that this second order nonlinear coupling dominates the growth of perturbations with masses \la 10^3 M_\odot immediately after recombination, enhancing the baryon fluctuation amplitudes by a factor of ~ 5, but the nonlinear effect does not persist at late times.

The eigen frequencies of radial pulsations of neutron stars are calculated in a strong magnetic field. At low densities we use the magnetic BPS equation of state(EOS) similar to that obtained by Lai and Shapiro while at high densities the EOS obtained from the relativistic nuclear mean field theory is taken and extended to include strong magnetic field. It is found that magnetised neutron stars support higher maximum mass where as the effect of magnetic field on radial stability for observed neutron star masses is minimal.

In this paper we study radial oscillations of rotating strange stars in strong magnetic fields in the Density Dependent Quark Mass (DDQM) model. We see that increase of frequency i.e. difference in frequency of rotating and non-rotating stars is more for higher magnetic fields. The change is small for low mass stars but it increases with the mass of the star. This change of frequency is significant for maximum mass whereas it is marginal for a 1.4 solar mass star.

We have studied phase transition from hadron matter to quark matter in the presence of high magnetic fields incorporating the trapped electron neutrinos at finite temperatures. We have used the density dependent quark mass (DDQM) model for the quark phase while the hadron phase is treated in the frame-work of relativistic mean field theory. It is seen that the nuclear energy at phase transition decreases with both magnetic field and temperature. A brief discussion of the effect of magnetic field in supernova explosions and proto-neutron star evolution is given.

We study the effect of magnetic field on the bulk viscosity of nuclear matter in neutron stars. We employ the framework of relativistic mean field theory to observe the dense nuclear matter in neutron stars. The effects are first studied for the case when the magnetic field does not exceed the critical value to confine the electrons to the lowest Landau levels. We then consider the case of intense magnetic field to evaluate viscosity for the URCA processes and show that the inequality $p_{F}(e)+p_{F}(p)\geq p_{F}(n)$ is no longer required to be satisfied for the URCA processes to proceed.

Cosmological models involving shear and rotation are considered, first in the General Relat ivistic and then in the Newtonian framework with the aim of investigating singularities in them by using numerical and analytical techniques. The dynamics of these rotating models ar e studied. It is shown that singularities are unavoidable in such models and that the centr ifugal force arising due to rotation can never overcome the gravitational and shearing forc e over a length of time.