We show that breaking of the conformal invariance of electromagnetic Lagrangian which is required for inflationary magnetogenesis arises naturally in the Poincar{\'e} Gauge Theory. We use the minimal coupling prescription to introduce the electromagnetic gauge fields as well as non-abelian gauge fields in this theory. Due to the addition of non-abelian gauge fields, we show that the solar constraints on this model can be naturally evaded. We find that in the minimal version of this model the generated magnetic field is too small to explain the observations. We discuss some generalizations of the gravitational action, including the Starobinsky model and a model with conformal invariance. We show that such generalizations naturally generate the kinetic energy terms required for magnetogenesis. We propose a generalization of the minimal model by adding a potential term, which is allowed within the framework of this model, and show that it leads to sufficiently large magnetic fields.

We study the alignment of radio galaxies axes using the FIRST catalogue. we impose several cuts in order to select the candidates which are most likely to be free of systematic bias. In our study we primarily focus on testing for alignment among sources within a certain angular separation from one another since for most sources redshift information is not available. We find a very significant effect for angular distances less than 1 degrees. The distance scale of alignment is found to be roughly 28 Mpc, in agreement with earlier estimates, assuming that these sources are dominantly at redshift of 0.8. However we are not able to entirely rule out the possibility of systematic bias in data. We also perform a full three dimensional analysis using a smaller data sample for which redshift information is available. In this case we only find a very weak signal at much larger distances.

Low energy ground-based cosmic ray air shower experiments generally have energy threshold in the range of a few tens to a few hundreds of TeV. The shower observables are measured indirectly with an array of detectors. The atmospheric absorption of low energy secondaries limits their detection frequencies at the Earth's surface. However, due to selection effects, a tiny fraction of low energy showers, which are produced in the lower atmosphere can reach the observational level. But, due to less information of shower observables, the reconstruction of these showers are arduous. Hence, it is believed that direct measurements by experiments aboard on satellites and balloon flights are more reliable at low energies. Despite having very small efficiency ($\sim$0.1%) at low energies, the large acceptance ($\sim$5 m$^2$sr) of GRAPES-3 experiment allows observing primary cosmic rays down below to $\sim$1 TeV and opens up the possibility to measure primary energy spectrum spanning from a few TeV to beyond cosmic ray knee (up to 10$^{16}$ eV), covering five orders of magnitude. The GRAPES-3 energy threshold for primary protons through Monte Carlo simulations are calculated, which gives reasonably good agreement with data. Furthermore, the total efficiencies and acceptance are also calculated for protons primaries. The ability of GRAPES-3 experiment to cover such a broader energy range may provide a unique handle to bridge the energy spectrum between direct measurements at low energies and indirect measurements at ultra-high energies.

We have analyzed AstroSat observations of the galactic micro-quasar system GRS 1915+105, when the system exhibited C-type Quasi-periodic Oscillations (QPOs) in the frequency range of 3.4-5.4 Hz. The broad band spectra (1-50 keV) obtained from simultaneous LAXPC and SXT can be well described by a dominant relativistic truncated accretion disk along with thermal Comptonization and reflection. We find that while the QPO frequency depends on the inner radii with a large scatter, a much tighter correlation is obtained when both the inner radii and accretion rate of the disk are taken into account. In fact, the frequency varies just as the dynamic frequency (i.e. the inverse of the sound crossing time) as predicted decades ago by the relativistic standard accretion disk theory for a black hole with spin parameter of ~0.9. We show that this identification has been possible due to the simultaneous broad band spectral coverage with temporal information as obtained from AstroSat.

We introduce a pixel space method to detect dipole modulation or hemispherical power asymmetry in the cosmic microwave background (CMB) polarization. The method relies on the use of squared total polarized flux whose ensemble average picks up a dipole due to the dipole modulation in the CMB polarization. The method is useful since it can be applied easily to partial sky. We define several statistics to characterize the amplitude of the detected signal. By simulations we show that the method can be used to reliably extract the signal at 2.7$\sigma$ level or higher in future CORE-like missions, assuming that the signal is present in the CMB polarization at the level detected by the Planck mission in the CMB temperature. An application of the method to the 2018 Planck data does not detect a significant effect, when taking into account the presence of correlated detector noise and residual systematics in data. Using the FFP10 we find the presence of a very strong bias which might be masking any real effect.

The angular resolution of an extensive air shower (EAS) array plays a critical role in determining its sensitivity for the detection of point $\gamma$-ray sources in the multi-TeV energy range. The GRAPES-3 an EAS array located at Ooty in India (11.4$^{\circ}$N, 76.7$^{\circ}$E, 2200 m altitude) is designed to study $\gamma$-rays in the TeV-PeV energy range. It comprises of a dense array of 400 plastic scintillators deployed over an area of 25000 m$^2$ and a large area (560 m$^2$) muon telescope. A new statistical method allowed real time determination of the propagation delay of each detector in the GRAPES-3 array. The shape of shower front is known to be curved and here the details of a new method developed for accurate measurement of the shower front curvature is presented. These two developments have led to a sizable improvement in the angular resolution of GRAPES-3 array. It is shown that the curvature depends on the size and age of an EAS. By employing two different techniques, namely, the odd-even and the left-right methods, independent estimates of the angular resolution are obtained. The odd-even method estimates the best achievable resolution of the array. For obtaining the angular resolution, the left-right method is used after implementing the size and age dependent curvature corrections. A comparison of the angular resolution as a function of EAS energy by these two methods shows them be virtually indistinguishable. The angular resolution of GRAPES-3 array is 47$^{\prime}$ for energies E$>$5 TeV and improves to 17$^{\prime}$ at E$>$100 TeV and finally approaching 10$^{\prime}$ at E$>$500 TeV.

We explore the large to moderate scale anisotropy in distant radio sources using the TGSS ADR1 catalog. We use different measures, i.e. number counts, sky brightness and flux per source, for this study. In agreement with earlier results, we report a significant excess of clustering signal above the angular scale of roughly $10$ degrees (i.e. $l\lessapprox 20-30$). We find that some survey areas have a systematically low/high flux and argue this may be the cause of the observed signal of excess power at low multipoles. With mocks we demonstrate the effect of such large scale flux systematics and recover TGSS like excess clustering signal by assuming $20\%$ flux uncertainties over $\sim 10^\circ \times 10^\circ$ size patches. We argue that that TGSS at this stage, i.e. TGSS ADR1, is not suitable for large scale clustering measurements. We find that the measure, flux per source, shows evidence of isotropy for all multipoles $l > 2$ despite the presence of systematics in the data.

We study precession dynamics of generic binary black holes in eccentric orbits using an effective potential based formalism derived in [M. Kesden et al., PRL {\bf 114}, 081103 (2015)]. This effective potential is used to classify binary black holes into three mutually exclusive spin morphologies. During the inspiral phase, binaries make transitions from one morphology to others. We evolve a population of binary black holes from an initial separation of $1000\mathbf{M}$ to a final separation of $10\mathbf{M}$ using post-Newtonian accurate evolution equations. We find that, given suitable initial conditions, a binary's eccentricity can follow one of three distinct evolutionary patterns: (i) eccentricity monotonically increasing until final separation, (ii) eccentricity rising after decaying to a minimum value, and (iii) eccentricity monotonically decreasing throughout the inspiral. The monotonic growth or growth after reaching a certain minimum of eccentricity is due to the effect of 2PN spin-spin coupling. Further, we investigate the morphology transitions in eccentric binaries and find that the probability of such binaries transiting from one to other is similar to those in circular orbits, implying that eccentricity plays a sub-dominant role in spin morphology evolution of a precessing binary black hole. We, hence, argue that the morphological classification of spin precession dynamics is a robust tool to constrain the formation channels of binaries with arbitrary eccentricity as well.

We investigate weak lensing effect on cosmic microwave background (CMB) in the presence of dipole anisotropy. The approach of flat-sky approximation is considered. We determine the functions $\sigma_0^2$ and $\sigma_2^2$ that appear in expressions of the lensed CMB power spectrum in the presence of a dipole anisotropy. We determine the correction to B-mode power spectrum which is found to be appreciable at low multipoles ($l$). However, the temperature and E-mode power spectrum are not altered significantly.

We test the isotropy of radio polarization angle orientations with a robust and reliable dual frequency polarimetric survey of active galactic nuclei (AGN). We find that the polarization orientations are consistent with the assumption of isotropy for scales larger than or equal to $\sim 800$ Mpc. This provides further evidence of isotropy at large distance scales and is likely to impose strong constraints on some of the physical mechanisms that may be invoked to explain past observations of alignment of radio and optical polarizations at large distance scales.

In this work, we present a $\sim$90 ks continuous monitoring of the Galactic micro-quasar GRS 1915+105 with AstroSat when the source undergoes a major transition from a non-variable, $\chi$ class (similar to radio-quiet $\chi$ class) to a structured, large amplitude, periodic heartbeat state (similar to $\rho$ class). We show that such transition takes place via an intermediate state when the large-amplitude, irregular variability of the order of hundreds of seconds in the soft X-ray band turned into a 100-150 sec regular, structured, nearly periodic flares. The properties of a strong low-frequency quasi-periodic oscillations (LF QPO) in the frequency range 3-5 Hz also evolve marginally during these variability transitions. We also study time-lag and rms spectra at the QPO and harmonic component and the dynamic power spectra. We note few important differences between the heartbeat state and the $\rho$ class. Interestingly, the time-averaged LF QPO properties in the hard X-ray band is relatively stable in three states when compared to the significant evolution observed in the slow variability properties at mHz frequencies. Such relative stability of LF QPOs implies the inner disk-corona coupled accretion flow which determines the LF QPO properties, may be uninterrupted by the launch of long, large-amplitude flares.

The balloon-borne HiCal radio-frequency (RF) transmitter, in concert with the ANITA radio-frequency receiver array, is designed to measure the Antarctic surface reflectivity in the RF wavelength regime. The amplitude of surface-reflected transmissions from HiCal, registered as triggered events by ANITA, can be compared with the direct transmissions preceding them by O(10) microseconds, to infer the surface power reflection coefficient $\cal{R}$. The first HiCal mission (HiCal-1, Jan. 2015) yielded a sample of 100 such pairs, resulting in estimates of $\cal{R}$ at highly-glancing angles (i.e., zenith angles approaching $90^\circ$), with measured reflectivity for those events which exceeded extant calculations. The HiCal-2 experiment, flying from Dec., 2016-Jan., 2017, provided an improvement by nearly two orders of magnitude in our event statistics, allowing a considerably more precise mapping of the reflectivity over a wider range of incidence angles. We find general agreement between the HiCal-2 reflectivity results and those obtained with the earlier HiCal-1 mission, as well as estimates from Solar reflections in the radio-frequency regime. In parallel, our calculations of expected reflectivity have matured; herein, we use a plane-wave expansion to estimate the reflectivity R from both a flat, smooth surface (and, in so doing, recover the Fresnel reflectivity equations) and also a curved surface. Multiplying our flat-smooth reflectivity by improved Earth curvature and surface roughness corrections now provides significantly better agreement between theory and the HiCal 2a/2b measurements.

We present X-ray spectral and timing behavior of Cyg X-3 as observed by AstroSat during the onset of a giant radio flare on 01-02 April 2017. Within a time-scale of few hours, the source shows a transition from the hypersoft state (HPS) to a more luminous state (we termed as the very high state) which coincides with the time of the steep rise in radio flux density by an order of magnitude. Modeling the SXT and LAXPC spectra jointly in 0.5-70.0 keV, we found that the first few hours of the observation is dominated by the HPS with no significant counts above 17 keV. Later, an additional flat powerlaw component suddenly appeared in the spectra which extends to very high energies with the powerlaw photon index of 1.49 +/- 0.04. Such a flat powerlaw component has never been reported from Cyg X-3. Interestingly the fitted powerlaw model in 25-70 keV, when extrapolated to the radio frequency, predicts the radio flux density consistent with the trend measured from RATAN-600 telescope at 11.2 GHz. This provides a direct evidence of the synchrotron origin of flat X-ray powerlaw component and the most extensive monitoring of the broadband X-ray behavior at the moment of decoupling the giant radio jet base from the compact object in Cyg X-3. Using SXT and LAXPC observations, we determine the giant flare ejection time as MJD 57845.34 +/- 0.08 when 11.2 GHz radio flux density increases from ~100 to ~478 mJy.

The hemispherical power asymmetry, observed in the CMBR data, has generally been interpreted in terms of the dipole modulation model for the temperature fluctuations. Here we point out that this model leads to several predictions, which can be directly tested in the current data. We suggest tests of the hemispherical power asymmetry both in real and multipole space. We find a significant signal of the dipole modulation model in WMAP and PLANCK data with our tests. The dipole amplitude and direction also agrees, within errors, with earlier results based on hemispherical analysis in multipole space. We also find evidence that the effective dipole modulation amplitude increases with the multipole l in the range l=2-64.

We propose a dipole modulation model for the Cosmic Microwave Background Radiation (CMBR) polarization field. We show that the model leads to correlations between l and l+1 multipoles, exactly as in the case of temperature. We obtain results for the case of TE, EE and BB correlations. An anisotropic or inhomogeneous model of primordial power spectrum which leads to such correlations in temperature field also predicts similar correlations in CMBR polarization. We analyze the CMBR temperature and polarization data in order to extract the signal of these correlation between l and l+1 multipoles. Our results for the case of temperature using the latest PLANCK data agree with those obtained by an earlier analysis. A detailed study of the correlation in the polarization data is not possible at present. Hence we restrict ourselves to a preliminary investigation in this case.

We consider an inhomogeneous model and independently an anisotropic model of primordial power spectrum in order to describe the observed hemispherical anisotropy in Cosmic Microwave Background Radiation. This anisotropy can be parametrized in terms of the dipole modulation model of the temperature field. Both the models lead to correlations between spherical harmonic coefficients corresponding to multipoles, l and l \pm 1. We obtain the model parameters by making a fit to TT correlations in CMBR data. Using these parameters we predict the signature of our models for correlations among different multipoles for the case of the TE and EE modes. These predictions can be used to test whether the observed hemispherical anisotropy can be correctly described in terms of a primordial power spectrum. Furthermore these may also allow us to distinguish between an inhomogeneous and an anisotropic model.

We argue that an anisotropic dipolar imaginary primordial power spectrum is possible within the framework of noncommutative space-times. We show that such a spectrum provides a good description of the observed dipole modulation in CMBR data. We extract the corresponding power spectrum from data. The dipole modulation is related to the observed hemispherical anisotropy in CMBR data, which might represent the first signature of quantum gravity.

The primary science goal of the NASA-sponsored ANITA project is measurement of ultra-high energy neutrinos and cosmic rays, observed via radio-frequency signals resulting from a neutrino- or cosmic ray- interaction with terrestrial matter (atmospheric or ice molecules, e.g.). Accurate inference of the energies of these cosmic rays requires understanding the transmission/reflection of radio wave signals across the ice-air boundary. Satellite-based measurements of Antarctic surface reflectivity, using a co-located transmitter and receiver, have been performed more-or-less continuously for the last few decades. Satellite-based reflectivity surveys, at frequencies ranging from 2--45 GHz and at near-normal incidence, yield generally consistent reflectivity maps across Antarctica. Using the Sun as an RF source, and the ANITA-3 balloon borne radio-frequency antenna array as the RF receiver, we have also measured the surface reflectivity over the interval 200-1000 MHz, at elevation angles of 12-30 degrees, finding agreement with the Fresnel equations within systematic errors. To probe low incidence angles, inaccessible to the Antarctic Solar technique and not probed by previous satellite surveys, a novel experimental approach ("HiCal-1") was devised. Unlike previous measurements, HiCal-ANITA constitute a bi-static transmitter-receiver pair separated by hundreds of kilometers. Data taken with HiCal, between 200--600 MHz shows a significant departure from the Fresnel equations, constant with frequency over that band, with the deficit increasing with obliquity of incidence, which we attribute to the combined effects of possible surface roughness, surface grain effects, radar clutter and/or shadowing of the reflection zone due to Earth curvature effects.

The cumulative number count $N$, of sources above a threshold is known to approximately follow a power law behaviour $N\propto S^{-x}$. We study the variation of spectral index $x$ across the sky in order to look for possible signals of violation of isotropy. We develop a rigorous algorithm of likelihood maximisation to accurately fit for the spectral index. We divide the sky into upper and lower hemispheres for a particular choice of $z$-axis and determine the difference $\Delta x$ between the best fit values of the spectral indices between the two hemispheres. The maximum value of this difference obtained by varying over the $z$-axis provides us with a measure of departure from isotropy. We find that the data support isotropy of the spectral index. The maximum difference is found to be 1.3\% of the full sky best fit value of $x$. The deviation is found to be significant only at 2$\sigma$ level which indicates a weak departure from isotropy. We also perform a dipole fit to the spectral index as a function of the angular coordinates. The result is found to be consistent with isotropy.

There currently exist many observations which are not consistent with the cosmological principle. We review these observations with a particular emphasis on those relevant for Square Kilometre Array (SKA). In particular, several different data sets indicate a preferred direction pointing approximately towards the Virgo cluster. We also observe a hemispherical anisotropy in the Cosmic Microwave Background Radiation (CMBR) temperature fluctuations. Although these inconsistencies may be attributed to systematic effects, there remains the possibility that they indicate new physics and various theories have been proposed to explain them. One possibility, which we discuss in this review, is the generation of perturbation modes during the early pre-inflationary epoch, when the Universe may not obey the cosmological principle. Better measurements will provide better constraints on these theories. In particular, we propose measurement of the dipole in number counts, sky brightness, polarized flux and polarization orientations of radio sources. We also suggest test of alignment of linear polarizations of sources as a function of their relative separation. Finally we propose measurement of hemispherical anisotropy or equivalently dipole modulation in radio sources.

We propose a generalized star product which deviates from the standard product when the fields at evaluated at different space-time points. This produces no changes in the standard Lagrangian density in non-commutative space-time but produces a change in the cosmological power spectrum. We show that the generalized star product leads to physically consistent results and can fit the observed data on hemispherical anisotropy in the cosmic microwave background radiation.

We explain the large scale correlations in radio polarization in terms of the correlations of galaxy cluster/supercluster magnetic field. Assuming that the polarization correlations closely follow the spatial correlations of the background magnetic field we recover the magnetic field spectral index as -2.74$\pm$0.04.This remarkably agrees with cluster magnetic field spectral index obtained in cosmological magneto-hydrodynamic simulations. We discuss possible physical scenarios in which the observed polarization alignment is plausible.

We consider a generalization of the standard model which respects quantum conformal invariance. This model leads to identically zero vacuum energy. We show how non-relativistic matter and dark energy arises in this model. Hence the model is shown to be consistent with observations.

We consider a conformal model involving two real scalar fields in which the conformal symmetry is broken by a soft mechanism and is not anomalous. One of these scalar fields is representative of the standard model Higgs. The model predicts exactly zero cosmological constant. In the simplest version of the model, some of the couplings need to be fine tuned to very small values. We formulate the problem of fine tuning of these couplings. We argue that the problem arises since we require a soft mechanism to break conformal symmetry. The symmetry breaking is possible only if the scalar fields do not evolve significantly over the time scale of the Universe. We present two solutions to this fine tuning problem. We argue that the problem is solved if the classical value of one of the scalar fields is super-Planckian, i.e. takes a value much larger than the Planck mass. The second solution involves introduction of a strongly coupled hidden sector that we call hypercolor. In this case the conformal invariance is broken dynamically and triggers the breakdown of the electroweak symmetry. We argue that our analysis applies also to the case of the standard model Higgs multiplet.

We relate the observed hemispherical anisotropy in the cosmic microwave radiation data to an inhomogeneous power spectrum model. The hemispherical anisotropy can be parameterized in terms of the dipole modulation model. This model leads to correlations between spherical harmonic coefficients corresponding to multipoles, l and l+1. We extract the $l$ dependence of the dipole modulation amplitude, A, by making a fit to the WMAP and PLANCK CMBR data. We propose an inhomogeneous power spectrum model and show that it also leads to correlations between multipoles, l and l+1. The model parameters are determined by making a fit to the data. The spectral index of the inhomogeneous power spectrum is found to be consistent with zero.

We study the dipole anisotropy in number counts and flux density weighted number counts {or sky brightness} in the NRAO VLA Sky Survey (NVSS) data. The dipole anisotropy is expected due to our local motion with respect to the CMBR rest frame. We analyse data with an improved fit to the number density, n(S), as a function of the flux density S, which allows deviation from a pure power law behaviour. We also impose more stringent cuts to remove the contribution due to clustering dipole. In agreement with earlier results, we find that the amplitude of anisotropy is significantly larger in comparison to the prediction based on CMBR measurements. The extracted speed is found to be roughly 3 times the speed corresponding to CMBR. The significance of deviation is smaller, roughly 2 sigma, in comparison to earlier estimates. For the cut, S>30 mJy, the speed is found to be $1110\pm370$ Km/s using the source count analysis. The direction of the dipole anisotropy is found to be approximately in agreement with CMBR. We find that the results are relatively insensitive to the lower as well as upper limit imposed on the flux density. Our results suggest that the Universe is intrinsically anisotropic with the axis of anisotropy axis pointing roughly towards the CMBR dipole direction. Finally we present a method which may allow an independent extraction of the local speed and an intrinsic dipole anisotropy, provided a larger data set becomes available in future.

We show that a theory with conformal invariance, which is explicitly broken by small terms, provides a solution to the fine tuning problem of the cosmological constant. In the absence of the symmetry breaking terms, the cosmological constant is zero. Its value in the full theory is controlled by the symmetry breaking terms. The symmetry breaking terms also provide the slow roll conditions, which may be useful in constructing a model of inflation.

We study the dipole anisotropy in integrated linearly polarization flux density in NRAO VLA Sky Survey (NVSS). We extract the anisotropy parameters in the number counts, number counts weighted by polarization observables, i.e. degree of polarization (p) and polarization flux (P). We also determine the anisotropy in the degree of polarization per source and polarization flux per source. We consider data with several different cuts on the flux density, S>10,20,30,40,50 mJy. For studies with degree of polarization we impose the additional cut 0.01<p<1. Similarly for polarization flux we impose, 0.5<P<100 mJy. We find a very significant signal of dipole, both in number counts and P (or p) weighted number counts. For the case of P weighted number counts the significance ranges from 8 \sigma to 5 \sigma for S>20 mJy to S>50 mJy. The corresponding direction parameters are found to be stable with the cut on flux density. The significance is even higher for the case of p weighted number counts. The observed anisotropy is found to be much larger in comparison to the CMBR expectations. We find that polarization observables show a much higher level of anisotropy in comparison to pure number counts or flux weighted number counts.

We present a detailed statistical analysis of the alignment of polarizations of radio sources at high redshift. We use the JVAS/CLASS 8.4-GHz surveys for our study. This study is motivated by the puzzling signal of alignment of polarizations from distant quasars at optical frequencies. We explore several different cuts on the polarization flux for our analysis. We find that the entire data shows a very significant signal of alignment on very large distance scales of order 500 Mpc. The alignment starts to decay only at much larger distances of order Gpc. If we only consider data with polarization flux greater than 1 mJy, we find alignment at distance scales less than 150 Mpc. We also find that data with polarization flux less than 0.5 mJy does not show significant alignment. Similar results are seen for data with degree of polarization less than 0.01, although here a mild signal of alignment is observed for a narrow range of angular separations. We argue that the signal cannot be explained in terms of bias due to systematic errors in removal of instrumental polarization. We also find that the degree of polarization shows a strong negative correlation with the total flux. The data appears to fall into two classes, one of which shows such a correlation. The remaining set, which has total flux greater than 100 mJy and degree of polarization lying between 0.01 and 0.1, shows a more random behaviour. The latter set is also found to show no alignment whereas the first set shows a very strong polarization alignment.

We study several anisotropic inflationary models and their implications for the observed violation of statistical isotropy in the CMBR data. In two of these models the anisotropy decays very quickly during the inflationary phase of expansion. We explicitly show that these models lead to violation of isotropy only for low l CMBR modes. Our primary aim is to fit the observed alignment of l=2,3 multipoles to the theoretical models. We use two measures, based on the power tensor, which contains information about the alignment of each multipole, to quantify the anisotropy in data. One of the measures uses the dispersion in eigenvalues of the power tensor. We also define another measure which tests the overall correlation between two different multipoles. We perturbatively compute these measures of anisotropy and fix the theoretical parameters by making a best fit to l=2,3 multipoles. We show that some of the models studied are able to consistently explain the observed violation of statistical isotropy.

We propose a model of the Universe based on Minkowski flat space-time metric. In this model the space-time does not evolve. Instead the matter evolves such that all the mass parameters increase with time. We construct a model based on unimodular gravity to show how this can be accomplished within the framework of flat space-time. We show that the model predicts the Hubble law if the masses increase with time. Furthermore we show that it fits the high z supernova data in a manner almost identical to the standard Big Bang model. Furthermore we show that at early times the Universe is dominated by radiative energy density. The phenomenon of recombination also arises in our model and hence predicts the existence of CMBR. However a major difference with the standard Big Bang is that the radiative temperature and energy density does not evolve in our model. Furthermore we argue that the basic motivation for inflation is absent in our model.

We consider a model of gravity and matter fields which is invariant only under unimodular general coordinate transformations (GCT). The determinant of the metric is treated as a separate field which transforms as a scalar under unimodular GCT. Furthermore we also demand that the theory is invariant under a new global symmetry which we call generalized conformal invariance. We study the cosmological implications of the resulting theory. We show that this theory gives a fit to the high-z supernova data which is identical to the standard Big Bang model. Hence we require some other cosmological observations to test the validity of this model. We also consider some models which do not obey the generalized conformal invariance. In these models we can fit the supernova data without introducing the standard cosmological constant term. Furthermore these models introduce only one dark component and hence solve the coincidence problem of dark matter and dark energy.

We present the results of three-dimensional simulations of quasar polarizations in the presence of pseudoscalar-photon mixing in the intergalactic medium. The intergalactic magnetic field is assumed to be uncorrelated in wave vector space but correlated in real space. Such a field may be obtained if its origin is primordial. Furthermore we assume that the quasars, located at cosmological distances, have negligible initial polarization. In the presence of pseudoscalar-photon mixing we show, through a direct comparison with observations, that this may explain the observed large scale alignments in quasar polarizations within the framework of big bang cosmology. We find that the simulation results give a reasonably good fit to the observed data.

We consider models of gravitation that are based on unimodular general coordinate transformations (GCT). These transformations include only those which do not change the determinant of the metric. We treat the determinant as a separate field which transforms as a scalar under unimodular GCT. We consider a class of such theories. In general, these theories do not transform covariantly under the full GCT. We characterize the violation of general coordinate invariance by introducing a new parameter. We show that the theory is consistent with observations for a wide range of this parameter. This parameter may serve as a test for possible violations of general coordinate invariance. We also consider the cosmic evolution within the framework of these models. We show that in general we do not obtain consistent cosmological solutions if we assume the standard cosmological constant or the standard form of non-relativistic matter. We propose a suitable generalization which is consistent with cosmology. We fit the resulting model to the high redshift supernova data. We find that we can obtain a good fit to this data even if include only a single component, either cosmological constant or non-relativistic matter.

We study the power asymmetry between even and odd multipoles in the multipolar expansion of CMB temperature data from WMAP, recently reported in the literature. We introduce an alternate statistic which probes this effect more sensitively. We find that the data is highly anomalous and consistently outside $2\sigma$ significance level in the whole multipole range $l=[2,101]$. We examine the possibility that this asymmetry may be caused by the foreground cleaning procedure or by residual foregrounds. By direct simulations, using the Planck Sky Model for foregrounds we rule out this possibility. We also examine several possible sub-dominant foregrounds, which might lead to such an asymmetry. However in all cases we are unable to explain the signal seen in data. We next examine cleaned maps, using procedures other than the one followed by the WMAP Science team. Specifically we analysed the maps cleaned by the IPSE procedure, Needlets and the harmonic ILC procedure. In all these cases we do not find a statistically significant signal of power asymmetry. This is in contrast to the result obtained by the WMAP best fit power spectrum as well as the ILC map. Finally, we test for the contribution of low-$l$ multipoles to the observed power asymmetry. We find that if we eliminate the first six multipoles, $l=[2,7]$, the significance falls below $2\sigma$ CL. Hence we find that the signal gets dominant contribution from low-$l$ modes.

We show that perturbations generated during the anisotropic pre-inflationary stage of cosmic evolution may affect cosmological observations today for a certain range of parameters. Due to the anisotropic nature of the universe during such early times, it might explain some of the observed signals of large scale anisotropy. In particular we argue that the alignment of CMB quadrupole and octopole may be explained by the Sachs-Wolfe effect due to the large scale anisotropic modes from very early times of cosmological evolution. We also comment on how the observed dipole modulation of CMB power may be explained within this framework.

We investigate the effects of pseudoscalar-photon mixing on electromagnetic radiation in the presence of correlated extragalactic magnetic fields. We model the Universe as a collection of magnetic domains and study the propagation of radiation through them. This leads to correlations between Stokes parameters over large scales and consistently explains the observed large-scale alignment of quasar polarizations at different redshifts within the framework of the big bang model.

We show that scale invariance provides a solution to the fine tuning problem of the cosmological constant. We construct a generalization of the standard model of particle physics which displays exact quantum scale invariance. The matter action is invariant under global scale transformations in arbitrary dimensions. However the gravitational action breaks scale invariance explicitly. The scale symmetry is broken spontaneously in the matter sector of the theory. We show that the contribution to the vacuum energy and hence the cosmological constant is identically zero from the matter sector within the full quantum theory. However the gravitational sector may give non-zero contributions to the cosmological constant at loop orders. No fine tuning may be required at loop orders since the matter sector gives zero contribution to the cosmological constant. We also show that we do not require full scale invariance in order to constrain the vacuum energy from the matter sector. We only require invariance under pseudoscale transformations. Using this idea and motivated by the concept of unimodular gravity we propose an alternative model. In this case also we show that matter sector gives exactly zero contribution to the vacuum energy.

We explore the possibility that a new physics interaction can provide an explanation for the knee just above $10^6$ GeV in the cosmic ray spectrum. We model the new physics modifications to the total proton-proton cross section with an incoherent term that allows for missing energy above the scale of new physics. We add the constraint that the new physics must also be consistent with published $pp$ cross section measurements, using cosmic ray observations, an order of magnitude and more above the knee. We find that the rise in cross section required at energies above the knee is radical. The increase in cross section suggests that it may be more appropriate to treat the scattering process in the black disc limit at such high energies. In this case there may be no clean separation between the standard model and new physics contributions to the total cross section. We model the missing energy in this limit and find a good fit to the Tibet III cosmic ray flux data. We comment on testing the new physics proposal for the cosmic ray knee at the Large Hadron Collider.

We analyze the effect of foregrounds on the observed alignment of CMBR quadrupole and octopole. The alignment between these multipoles is studied by using a symmetry based approach which assigns a principal eigenvector (PEV) or an axis with each multipole. We determine the significance of alignment between these multipoles by using the Internal Linear Combination (ILC) 5 and 7 year map s and also the maps obtained by using the Internal Power Spectrum Estimation (IPSE) procedure. The effect of foreground cleaning is studied in detail within the framework of the IPSE method both analytically and numerically. By using simulated CMBR data, we study how the PEVs of the pure simulated CMB map differ from those of the final cleaned map. We find that, in general, the shift in the PEVs is relatively small and in random directions. Due to the random nature of the shift we conclude that it can only lead to misalignment rather than alignment of multipoles. We also directly estimate the significance of alignment by using simulated cleaned maps. We find that the results in this case are identical to those obtained by simple analytic estimate or by using simulated pure CMB maps.

We compute the flux of axions from Active Galactic Nuclei (AGN). Axions can be produced in the accretion disk by the Compton, Bremsstrahlung and Primakoff processes. We find that the axion luminosity due to these processes is negligible in comparison to the photon luminosity from AGNs. We also compute the luminosity of a hypothetical pseudoscalar, with very small mass, from the AGN atmosphere due to the phenomenon of pseudoscalar-photon mixing in background magnetic field. In this case we find that for some parameter ranges, the pseudoscalar flux can exceed that of photons. We comment on the implications of this result on the observed large scale alignment of optical polarizations from AGNs.

We argue that scale invariance is not anomalous in quantum field theory, provided it is broken cosmologically. We consider a locally scale invariant extension of the Standard Model of particle physics and argue that it fits both the particle and cosmological observations. The model is scale invariant both classically and quantum mechanically. The scale invariance is broken cosmologically producing all the dimensionful parameters. The cosmological constant or dark energy is a prediction of the theory and can be calculated systematically order by order in perturbation theory. It is expected to be finite at all orders. The model does not suffer from the hierarchy problem due to absence of scalar particles, including the Higgs, from the physical spectrum.

We estimate CMB polarization and temperature power spectra using WMAP 5-year foreground contaminated maps. The power spectrum is estimated by using a model independent method, which does not utilize directly the diffuse foreground templates nor the detector noise model. The method essentially consists of two steps, (i) removal of diffuse foregrounds contamination by making linear combination of individual maps in harmonic space and (ii) cross-correlation of foreground cleaned maps to minimize detector noise bias. For temperature power spectrum we also estimate and subtract residual unresolved point source contamination in the cross-power spectrum using the point source model provided by the WMAP science team. Our $TT$, $TE$ and $EE$ power spectra are in good agreement with the published results of the WMAP science team. The error bars on the polarization power spectra, however, turn out to be smaller in comparison to what is obtained by the WMAP science team. We perform detailed numerical simulations to test for bias in our procedure. We find that the bias is small in all cases. A negative bias at low $l$ in $TT$ power spectrum has been pointed in an earlier publication. We find that the bias corrected quadrupole power ($l(l+1)C_l/2\pi)$ is 532 $\mu K^2$, approximately 2.5 times the estimate (213.4 $\mu K^2$) made by the WMAP team.

We compute the cosmological constant in a scale invariant scalar field theory. The gravitational action is also suitably modified to respect scale invariance. Due to scale invariance the theory does not admit a cosmological constant term. The scale invariance is broken by a recently introduced mechanism called cosmological symmetry breaking. This leads to a nonzero cosmological constant. We compute the one loop corrections to the cosmological constant and show that it is finite.

In this paper, we propose & implement on WMAP 5-year data, a model independent approach of foreground power spectrum estimation for multifrequency observations of CMB experiments. Recently a model independent approach of CMB power spectrum estimation was proposed by Saha et al. 2006. This methodology demonstrates that CMB power spectrum can be reliably estimated solely from WMAP data without assuming any template models for the foreground components. In the current paper, we extend this work to estimate the galactic foreground power spectrum using the WMAP 5 year maps following a self contained analysis. We apply the model independent method in harmonic basis to estimate the foreground power spectrum and frequency dependence of combined foregrounds. We also study the behaviour of synchrotron spectral index variation over different regions of the sky. We compare our results with those obtained from MEM foreground maps which are formed in pixel space. We find that relative to our model independent estimates MEM maps overestimates the foreground power close to galactic plane and underestimates it at high latitudes.

Recently a symmetry-based method to test for statistical isotropy of the cosmic microwave background was developed. We apply the method to template-cleaned 3-year and 5-year WMAP-$DA$ maps. We examine a wide range of angular multipoles from $2 < l < 300$. The analysis detects statistically signicant signals of anisotropy inconsistent with an isotropic CMB in some of the foreground cleaned maps. We are unable to resolve whether the anomalies have a cosmological, local astrophysical or instrumental origin. Assuming the anisotropy arises due to residual foreground contamination, we estimate the residual foreground power in the maps. For the W band maps, we also find a highly improbable degree of isotropy we cannot explain. We speculate that excess isotropy may be caused by faulty modeling of detector noise.

We study the possibility of evading astrophysical bound on light pseudoscalars. We argue that the solar bounds can be evaded if we have a sufficiently strong self-coupling of the pseudoscalars. The required coupling does not conflict with any known experimental bounds. We also show that it is possible to find a coupling range such that the recent PVLAS experiment is not in conflict with any astrophysical bounds.

Pseudoscalar particles, with almost zero mass and very weak coupling to the visible matter, arise in many extensions of the standard model of particle physics. Their mixing with photons in the presence of an external magnetic field leads to many interesting astrophysical and cosmological consequences. This mixing depends on the medium properties, the momentum of the photon and the background magnetic field. Here we give a general treatment of pseudoscalar-photon oscillations in a background magnetic field, taking the Faraday term into account. We give predictions valid in all regimes, under the assumption that the frequency of the wave is much higher than the plasma frequency of the medium. At sufficiently high frequencies, the Faraday effect is negligible and we reproduce the standard pseudoscalar-photon mixing phenomenon. However at low frequencies, where Faraday effect is important, the mixing formulae are considerably modified. We explicitly compute the contribution due to the longitudinal mode of the photon and show that it is negligible.

We consider a generalization of Einstein's general theory of relativity such that it respects local scale invariance. This requires the introduction of a scalar and a vector field in the action. We show that the theory naturally displays both dark energy and dark matter. We solve the resulting equations of motion assuming an FRW metric. The solutions are found to be almost identical to those corresponding to the standard $\Lambda$CDM model

We model the effect of photon and ultra-light pseudoscalar mixing on the propagation of electromagnetic radiation through the extragalactic medium. The medium is modelled as a large number of magnetic domains, uncorrelated with one another. We obtain an analytic expression for the different Stokes parameters in the limit of small mixing angle. The different Stokes parameters are found to increase linearly with the number of domains. We also verify this result by direct numerical simulations. We use this formalism to estimate the effect of pseudoscalar-photon mixing on the Cosmic Microwave Background (CMB) polarization. We impose limits on the model parameters by the CMB observations. We find that the currently allowed parameter range admits a CMB circular polarization up to order $10^{-7}$.

We generalize the standard model of particle physics such it displays global scale invariance. The gravitational action is also suitably modified such that it respects this symmetry. This model is interesting since the cosmological constant term is absent in the action. We find that the scale symmetry is broken by the recently introduced cosmological symmetry breaking mechanism. This simultaneously generates all the dimensionful parameters such as the Newton's gravitational constant, the particle masses and the vacuum or dark energy. We find that in its simplest version the model predicts the Higgs mass to be very small, which is ruled out experimentally. We further generalize the model such that it displays local scale invariance. In this case the Higgs particle disappears from the particle spectrum and instead we find a very massive vector boson. Hence the model gives a consistent description of particle physics phenomenology as well as fits the cosmological dark energy.

We introduce new symmetry-based methods to test for isotropy in cosmic microwave background radiation. Each angular multipole is factored into unique products of power eigenvectors, related multipoles and singular values that provide 2 new rotationally invariant measures mode by mode. The power entropy and directional entropy are new tests of randomness that are independent of the usual CMB power. Simulated galactic plane contamination is readily identified, and the new procedures mesh perfectly with linear transformations employed for windowed-sky analysis. The ILC -WMAP data maps show 7 axes well aligned with one another and the direction Virgo. Parameter free statistics find 12 independent cases of extraordinary axial alignment, low power entropy, or both having 5% probability or lower in an isotropic distribution. Isotropy of the ILC maps is ruled out to confidence levels of better than 99.9%, whether or not coincidences with other puzzles coming from the Virgo axis are included. Our work shows that anisotropy is not confined to the low l region, but extends over a much larger l range.

Radiation propagating over cosmological distances can probe light weakly interacting pseudoscalar (or scalar) particles. The existence of a spin-0 field changes the dynamical symmetries of electrodynamics. It predicts spontaneous generation of polarization of electromagnetic waves due to mode mixing in the presence of background magnetic field.We illustrate this by calculations of propagation in a uniform medium, as well as in a slowly varying background medium, and finally with resonant mixing. Highly complicated correlations between different Stokes parameters are predicted depending on the parameter regimes. The polarization of propagating waves shows interesting and complex dependence on frequency, the distance of propagation, coupling constants, and parameters of the background medium such as the plasma density and the magnetic field strength. For the first time we study the resonant mixing of electromagnetic waves with the scalar field, which occurs when the background plasma frequency becomes equal to the mass of the scalar field at some point along the path. Dynamical effects are found to be considerably enhanced in this case. We also formulate the condition under which the adiabatic approximation can be used consistently, and find caveats about comparing different frequency regimes.

In recent years the goal of estimating different cosmological parameters precisely has set new challenges in the effort to accurately measure the angular power spectrum of CMB. This has required removal of foreground contamination as well as detector noise bias with reliability and precision. Recently, a novel {\em model-independent} method for the estimation of CMB angular power spectrum solely from multi-frequency observations has been proposed and implemented on the first year WMAP data by Saha et al. 2006. All previous estimates of power spectrum of CMB are based upon foreground templates using data sets from different experiments. However our methodology demonstrates that {\em CMB angular spectrum can be reliably estimated with precision from a self contained analysis of the WMAP data}. In this work we provide a detailed description of this method. We also study and identify the biases present in our power spectrum estimate. We apply our methodoly to extract the power spectrum from the WMAP 1 year and 3 year data.

The energy density of the universe today may be dominated by the vacuum energy of a slowly rolling scalar field. Making a quantum expansion around such a time dependent solution is found to break fundamental symmetries of quantum field theory. We call this mechanism cosmological symmetry breaking and argue that it is different from the standard phenomenon of spontaneous symmetry breaking. We illustrate this with a toy scalar field theory, whose action displays a U(1) symmetry. We identify a symmetry, called pseudo-scale invariance, which sets the cosmological constant exactly equal to zero, both in classical and quantum theory. This symmetry is also broken cosmologically and leads to a nonzero vacuum or dark energy. The slow roll condition along with the observed value of dark energy leads to a value of the background scalar field of the order of Planck mass. We also consider a U(1) gauge symmetry model. Cosmological symmetry breaking, in this case, leads to a non zero mass for the vector field. We also show that a cosmologically broken pseudo-scale invariance can generate a wide range of masses.

Weakly coupled, almost massless, spin 0 particles have been predicted by many extensions of the standard model of particle physics. Recently, the PVLAS group observed a rotation of polarization of electromagnetic waves in vacuum in the presence of transverse magnetic field. This phenomenon is best explained by the existence of a weakly coupled light pseudoscalar particle. However, the coupling required by this experiment is much larger than the conventional astrophysical limits. Here we consider a hypothetical self-interacting pseudoscalar particle which couples weakly with visible matter. Assuming that these pseudoscalars pervade the galaxy, we show that the solar limits on the pseudoscalar-photon coupling can be evaded.

We examine the Type 1a supernova data in order to determine if it shows any signal of large scale anisotropy. The anisotropy is modelled by an extended G\"{o}del metric, which incorporates expansion along with rotation. The model is smoothly connected to the usual FRW type, while expressing anisotropic metric effects depending on certain parameters. We find no significant signal of anisotropy in the data. We obtain bounds on an anisotropic redshift versus magnitude relationship, and accompanying parameters of the G\"{o}del-Obukhov metric.

We revisit the classification of polarization observables of the cosmic microwave background. There exists a unified approach to the $3 \times 3$ density matrix by which intensity, linear and circular polarization are treated on an equal footing. The representations of the rotation group on the sphere contain certain right-left setup symmetries which have not been taken into account before. Left-right symmetries revise the construction of invariants and certain predictions based on symmetries. There are two true rotational invariants mode by mode in polarization data. Independent of models of cosmological perturbations, we emphasize methods to construct invariant distributions in order to test predictions of isotropy, axial or parity symmetries, and generate covariant and invariant statistics order by order.

We analyze the three-year WMAP temperature anisotropy data seeking to confirm the power spectrum and likelihoods published by the WMAP team. We apply five independent implementations of four algorithms to the power spectrum estimation and two implementations to the parameter estimation. Our single most important result is that we broadly confirm the WMAP power spectrum and analysis. Still, we do find two small but potentially important discrepancies: On large angular scales there is a small power excess in the WMAP spectrum (5-10% at l<~30) primarily due to likelihood approximation issues between 13 <= l <~30. On small angular scales there is a systematic difference between the V- and W-band spectra (few percent at l>~300). Recently, the latter discrepancy was explained by Huffenberger et al. (2006) in terms of over-subtraction of unresolved point sources. As far as the low-l bias is concerned, most parameters are affected by a few tenths of a sigma. The most important effect is seen in n_s. For the combination of WMAP, Acbar and BOOMERanG, the significance of n_s =/ 1 drops from ~2.7 sigma to ~2.3 sigma when correcting for this bias. We propose a few simple improvements to the low-l WMAP likelihood code, and introduce two important extensions to the Gibbs sampling method that allows for proper sampling of the low signal-to-noise regime. Finally, we make the products from the Gibbs sampling analysis publically available, thereby providing a fast and simple route to the exact likelihood without the need of expensive matrix inversions.

The remarkable improvement in the estimates of different cosmological parameters in recent years has been largely spearheaded by accurate measurements of the angular power spectrum of Cosmic Microwave Background (CMB) radiation. This has required removal of foreground contamination as well as detector noise bias with reliability and precision. Recently, a novel model-independent method for the estimation of CMB angular power spectrum from multi-frequency observations has been proposed and implemented on the first year WMAP (WMAP-1) data by Saha et al.~2006. We review the results from WMAP-1 and also present the new angular power spectrum based on three years of the WMAP data (WMAP-3). Previous estimates have depended on foreground templates built using extraneous observational input to remove foreground contamination. This is the first demonstration that the CMB angular spectrum can be reliably estimated with precision from a self contained analysis of the WMAP data. The primary product of WMAP are the observations of CMB in 10 independent difference assemblies (DA) distributed over 5 frequency bands that have uncorrelated noise. Our method utilizes maximum information available within WMAP data by linearly combining DA maps from different frequencies to remove foregrounds and estimating the power spectrum from the 24 cross power spectra of clean maps that have independent noise. An important merit of the method is that the expected residual power from unresolved point sources is significantly tempered to a constant offset at large multipoles (in contrast to the ~l^2 contribution expected from a Poisson distribution) leading to a small correction at large multipoles. Hence, the power spectrum estimates are less susceptible to uncertainties in the model of point sources.

We study the possibility of evading astrophysical bounds on light pseudoscalars. We argue that the solar bounds can be evaded if we have a sufficiently strong self coupling of the pseudoscalars. The required couplings do not conflict with any known experimental bounds. We show that it is possible to find a coupling range such that the results of the recent PVLAS experiment are not in conflict with any astrophysical bounds.

Accurate measurements of angular power spectrum of Cosmic Microwave Background (CMB) radiation has lead to marked improvement in the estimates of different cosmological parameters. This has required removal of foreground contamination as well as detector noise bias with reliability and precision. We present the estimation of CMB angular power spectrum from the multi-frequency observations of WMAP using a novel {\em model-independent} method. The primary product of WMAP are the observations of CMB in 10 independent difference assemblies (DA) that have uncorrelated noise. Our method utilizes maximum information available within WMAP data by linearly combining all the DA maps to remove foregrounds and estimating the power spectrum from cross power spectra of clean maps with independent noise. We compute 24 cross power spectra which are the basis of the final power spectrum. The binned average power matches with WMAP team's published power spectrum closely. A small systematic difference at large multipoles is accounted for by the correction for the expected residual power from unresolved point sources. The correction is small and significantly tempered. Previous estimates have depended on foreground templates built using extraneous observational input. This is the first demonstration that the CMB angular spectrum can be reliably estimated with precision from a *self contained analysis* of the WMAP data.

We show that distances of objects at cosmological distances can be measured directly using interferometry. Our approach to interferometric parallax comes from analysis of 4-point amplitude and intensity correlations that can be generated from pairs of well-separated detectors. The baseline required to measure cosmological distances of Gigaparsec order are within the reach of the next generation of space-borne detectors. The semi-classical interpretation of intensity correlations uses a notion of a single photon taking two paths simultaneously. Semi-classically a single photon can simultaneously enter four detectors separated by an astronomical unit, developing correlations feasible to measure with current technology.

Type 1a supernova magnitudes conventionally include an additive parameter called the extinction coefficient. We find that the extinction coefficients of a popular ``gold'' set are well correlated with the deviation of magnitudes from Hubble diagrams. If the effect is due to bias, extinctions have been overestimated, which makes supernovas appear more dim. The statistical significance of the extinction-acceleration correlation has a random chance probability of less than one in a million. The hypothesis that extinction coefficients should be corrected empirically provides greatly improved fits to both accelerating and non-accelerating models, with the independent feature of eliminating any significant correlation of residuals.

Origin of ultra high energy cosmic rays is an unsolved problem. Several proposals such as Z-burst, decay of super massive matter, susy particles as a primary, neutrino as a primary in extra dimension models exist in the literature which try to address this issue. Many of these proposals solve the problem of propagation of cosmic rays over cosmological distances by introducing new physics. However these do not explain the origin of such high energy cosmic rays. The possible astrophysics sites, such as active galactic nuclei, are highly constrained. Here we determine whether these cosmic rays originated from the decay of some exotic objects, such as primordial black holes (PBHs), present in the early universe. In contrast to the usual Top Down scenario we do not assume that this exotic object necessarily has to decay in our astrophysical neighbourhood since we assume a beyond the standard model scenario, where the propagation problem is absent. We consider the standard 4-dimension PBHs as well as the brane world PBHs. We find that in both cases it is unable to produce the observed ultra high energy cosmic ray flux.

We point out that a well known axion model with an explicit Z(N) symmetry breaking term predicts both dark energy and cold dark matter. We estimate the parameters of this model which fit the observed densities of the dark components of the universe. We find that the parameters do not conflict with any observations.

We reconsider analysis of data on the cosmic microwave background on the largest angular scales. Temperature multipoles of any order factor naturally into a direct product of axial quantities and cosets. Striking coincidences exist among the axes associated with the dipole, quadrupole, and octupole CMB moments. These axes also coincide well with two other axes independently determined from polarizations at radio and optical frequencies propagating on cosmological scales. The five coincident axes indicate physical correlation and anisotropic properties of the cosmic medium not predicted by the conventional Big Bang scenario. We consider various mechanisms, including foreground corrections, as candidates for the observed correlations. We also consider whether the propagation anomalies may be a signal of ``dark energy'' in the form of a condensed background field. Perhaps {\it light propagation} will prove to be an effective way to look for the effects of {\it dark energy}.

Radiation traversing the observable universe provides powerful ways to probe anisotropy of electromagnetic propagation. A controversial recent study claimed a signal of dipole character. Here we test a new and independent data set of 361 points under the null proposal of {\it statistical independence} of linear polarization alignments relative to galaxy axes, versus angular positions. The null hypothesis is tested via maximum likelihood analysis of best fits among numerous independent types of factored distributions. We also examine single-number correlations which are parameter free, invariant under coordinate transformations, and distributed very robustly. The statistics are shown explicitly not to depend on the uneven distribution of sources on the sky. We find that the null proposal is not supported at the level of less than 5% to less than 0.1% by several independent statistics. The signal of correlation violates parity, that is, symmetry under spatial inversion, and requires a statistic which transforms properly. The data indicate an axis of correlation, on the basis of likelihood determined to be $[{\rm R.A.}=(0^{\rm h},9^{\rm m}) \pm (1^{\rm h},0^{\rm m})$, ${\rm Decl.} = -1^o\pm 15^o]$.

We present a detailed statistical study of the observed anisotropy in radio polarizations from distant extragalactic objects. This anisotropy was earlier found by Birch (1982) and reconfirmed by Jain and Ralston (1999) in a larger data set. A very strong signal was seen after imposing the cut $|RM-\bar{RM}|>6$ rad/m$^2$, where $RM$ is the rotation measure and $\bar{RM}$ its mean value. In this paper we show that there are several indications that this anisotropy cannot be attributed to bias in the data. We also find that a generalized statistic shows a very strong signal in the entire data without imposing the RM dependent cut. Finally we argue that an anisotropic background pseudoscalar field can explain the observations.

We introduce several coordinate invariant statistical procedures in order to test for local alignment of polarizations. A large scale alignment of optical polarizations from distant QSOs has recently been observed by \huts and collaborators. The new statistical procedures are based on comparing polarizations at different angular coordinates by making a parallel transport. The results of these statistical procedures continue to support the existence of the large scale alignment effect in the QSO optical polarization data. The alignment is found to be much more pronounced in the data sample with low degrees of polarization $p\le 2$%. This suggests that the alignment may be attributed to some propagation effect. The distance scale over which the alignment effect is dominant is found to be of order 1 Gpc. We also find that a very large scale alignment is present in the large redshift, $z\ge 1$, data sample. Infact the data sample with $z\ge 1$ appears to be aligned over the entire celestial sphere. We discuss possible physical effects, such as extinction and pseudoscalar-photon mixing, which may be responsible for the observations.

Cosmic rays in the regime of 10^19-10^20 eV should not exist, according to the GZK bound. Conflicts between data and the bound have led to suggestions that laws of physics must be modified. We find that physical processes neglected in the early history terminate the LPM effect, important in calibrating the energy of electromagnetic and hadronic showers in the GZK energy regime. The processes neglected are the direct losses of electrons and positrons, overlooked by LPM and greatly augmented by ultra-high energy circumstances. We present an exact formula for direct energy losses of e^+e^- pairs in electromagnetic showers. Numerical integration in current hadronic models shows that actual energy losses are always substantially larger than LPM estimates, but the energy at which crossover occurs depends on models.

Cubic kilometer neutrino telescopes are capable of probing fundamental questions of ultra-high energy neutrino interactions. There is currently great interest in neutrino interactions caused by low-scale, extra dimension models. Above 1 PeV the cross section in low scale gravity models rises well above the total Standard Model cross section. We assess the observability of this effect in the 1 PeV - 100 PeV energy range of kilometer-scale detectors with several new points of emphasis that hinge on enhanced neutral current cross sections. A major point is the importance of ``feed-down'' regeneration of upward neutrino flux, driven by new-physics neutral current interactions in the flux evolution equations. Feed-down is far from negligible, and it is essential to include its effect. We then find that the angular distribution of events has high discriminating value in separating models. In particular the ``up-to-down'' ratio between upward and downward-moving neutrino fluxes is a practical diagnostic tool which can discriminate between models in the near future. The slope of the angular distribution, in the region of maximum detected flux, is also substantially different in low-scale gravity and the Standard Model. These observables are only weakly dependent on astrophysical flux uncertainties. We conclude that angular distributions can reveal a breakdown of the Standard Model and probe the new physics beyond, as soon as data become available.

We investigate the effect of axions on the polarization of electromagnetic waves as they propagate through astronomical distances. We analyze the change in the dispersion of the electromagnetic wave due to its mixing with axions. We find that this leads to a shift in polarization and turns out to be the dominant effect for a wide range of frequencies. We analyze whether this effect or the decay of photons into axions can explain the large scale anisotropies which have been observed in the polarizations of quasars and radio galaxies. We also comment on the possibility that the axion-photon mixing can explain the dimming of distant supernovae.

Low energy string theory predicts the existence of an axion field which can lead to cosmic birefringence. We solve the electromagnetic wave equations in the presence of such an axion and a dilaton field in order to determine their effect on the polarization of light. We find that the presence of dilaton field leads to a nontrivial modification of the final result. We comment on the possibility of discovering such an effect by observations of radio wave polarizations from distant radio galaxies and quasars. We have also determined the limits on the string theory parameters that are imposed by the current polarization data from distant quasars and radio galaxies.

In models based on large extra dimensions where massive spin 2 exchange can dominate at high energies, the neutrino-proton cross section can rise to typical hadronic values at energies above 10^20 eV. The neutrino then becomes a candidate for the primary that initiates the highest energy cosmic ray showers. We investigate characteristics of neutrino-induced showers compared to proton-induced showers. The comparison includes study of starting depth, profile with depth, lateral particle distribution at ground and muon lateral distribution at ground level. We find that for cross sections above 20 mb there are regions of parameter space where the two types of showers are essentially indistinguishable. We conclude that the neutrino candidate hypothesis cannot be ruled out on the basis of shower characteristics.

Angular correlations of ultra-high energy cosmic rays with cosmologically distant sources may provide clues to these mysterious events. We compare cosmic ray tracks with energies above $ 10^{20}eV$ to a compilation of radio-loud compact QSO positions. The statistical method emphasizes invariant quantities and a test of {\it statistical independence} of track and source distributions. Statistical independence is ruled out by several independent statistics at confidence levels of less that $10^{-3}$ ($99.9% $.), indicating that track directions and QSO source positions are correlated at a highly significant level.

We propose a new method for the extraction of Rotation Measure from spectral polarization data. The method is based on maximum likelihood analysis and takes into account the circular nature of the polarization data. The method is unbiased and statistically more efficient than the standard $\chi^2$ procedure. We also find that the method is computationally much faster than the standard $\chi^2$ procedure if the number of data points are very large.

Cosmic ray events above $10^{20}$ eV are on the verge of confronting fundamental particle physics. The neutrino is the only candidate primary among established particles capable of crossing 100 Mpc intergalactic distances unimpeded. The magnitude of $\nu N$ cross sections indicated by events, plus consistency with the Standard Model at low-energy, point to new physics of massive spin-2 exchange. In models based on extra dimensions, we find that the $\nu N$ cross section rises to typical hadronic values of between 1 and 100 mb at energies above $10^{20}$ eV. Our calculations take into account constraints of unitarity. We conclude that air-showers observed with energies above $10^{19}$ eV are consistent with neutrino primaries and extra-dimension models. An {\it upper bound} of 1-10 TeV on the mass scale at which graviton exchange becomes strong in current Kaluza-Klein models follows.

We report on the feasibility of using an isotropic flux of cosmic neutrinos in the energy range of 10 to 10,000 TeV to study the interior geophysical structure of the Earth. The angular distribution of events in a $\sim {\rm km}^3$-scale neutrino telescope can be inverted to yield the Earth's nucleon density distribution. The approach is independent of previous seismic methods. The energy spectrum of the neutrino primaries can also be determined from consistency with the angular distribution. Depending on neutrino flux, a reasonable density determination can be obtained with a few hundred receivers operating for a few years.

We study the feasibility of using an isotropic flux of cosmic neutrinos in the energy range of 10 to 10000 TeV to study the interior structure of Earth. The angular distribution of events in a $\sim {\rm km}^3$-scale neutrino telescope can be inverted to yield information on the Earth's mass distribution that is independent of other methods. The energy spectrum of the neutrino primaries is also determined from consistency with the angular distribution. It is possible to make a model independent determination of the density profile of Earth's interior, separate from the absolute normalization of the incident cosmic neutrinos.

Polarizations of electromagnetic waves from distant galaxies are known to be correlated with the source orientations. These quantities have been used to search for signals of cosmological birefringence. We review and classify transformation properties of the polarization and source orientation observables. The classifications give a firm foundation to certain practices which have sprung up informally in the literature. Transformations under parity play a central role, showing that parity violation in emission or in the subsequent propagation is an observable phenomenon. We also discuss statistical measures, correlations and distributions which transform properly and which can be used for systematic data analysis.

We discuss a new mechanism which could cause a rotation of polarization of electromagnetic waves due to magnetic fields on cosmological scales. The effect is due to the geometrical phase of Pancharatnam and Berry, and causes a corkscrew twisting of the plane of polarization. The new effect represents an additional tool that allows possible intergalactic and cosmological magnetic fields to be studied using radio propagation.