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.

Ever since Cosmic Microwave Background (CMB) signal is being measured by various satellites based observations with increasing experimental accuracies there has been a parallel increase in the demand for a CMB reconstruction technique which can provide accurate estimates of CMB signal and the theoretical angular power spectrum along with reliable statistical error estimates associated with them. In this work, we estimate the joint posterior of CMB E mode signal (S) and corresponding theoretical angular power spectrum (C^E_l) over large angular scales given the simulated polarization observations of future generation COrE satellite mission. To generate samples from the joint distribution we employ the ILC technique with prior information of CMB E mode covariance matrix augmented by a Gibbs sampling technique. We estimate the marginalized densities of S and C^E_l using the samples from full-posterior. The best fit cleaned E mode map and the corresponding angular power spectrum agree well with the input E mode map and the sky power spectrum implying accurate reconstruction using COrE like observations. Using the samples C^E_l of all Gibbs chains we estimate the likelihood function P(C^E_l|D) of any arbitrary C^E_l given simulated observed maps (D) of COrE mission following Blackwell-Rao estimator. The likelihood function can be seamlessly integrated to the cosmological parameter estimation method. Apart from producing an accurate estimate of E mode signal over large angular scales our method also builds a connection between the component reconstruction and reliable cosmological parameter estimation using CMB E mode observations over large angular scales. The entire method does not assume any explicit models for E mode foreground components in order to remove them, which is an attractive property since foreground modelling uncertainty does not pose as a challenge in this case.

In the recent years, the internal-linear-combination (ILC) method was investigated extensively in the context of reconstruction of Cosmic Microwave Background (CMB) temperature anisotropy signal using observations obtained by WMAP and Planck satellite missions. In this article, we, for the first time, apply the ILC method to reconstruct the large scale CMB E mode polarization signal, which serves as the unique probe of ionization history of the Universe, using simulated observations of 15 frequency CMB polarization maps of future generation COrE satellite mission. We find that usual ILC cleaned E mode map is highly erroneous due to presence of a chance-correlation between CMB and astrophysical foreground components in the empirical covariance matrix which is used to estimate the weight factors. The cleaned angular power spectrum for E mode is strongly biased and erroneous due to these chance correlation factors. In order to address the issues of bias and errors we extend and improve the usual ILC method for CMB E mode reconstruction by incorporating prior information of theoretical E mode angular power spectrum while estimating the weights for linear combination of input maps. Using the E mode covariance matrix effectively suppresses the CMB-foreground chance correlation power leading to an accurate reconstruction of cleaned CMB E mode map and its angular power spectrum. We provide a comparative study of the performance of the usual ILC and the new method over large angular scales of the sky and show that the later produces significantly statistically improved results than the former. The new E mode CMB angular power spectrum contains neither any significant negative bias at the low multipoles nor any positive foreground bias at relatively higher mutlipoles. The error estimates of the cleaned spectrum agree very well with the cosmic variance induced error.

In the era of precision cosmology accurate estimation of cosmological parameters is based upon the implicit assumption of Gaussian nature of Cosmic Microwave Background (CMB) radiation. Therefore, an important scientific question to ask is whether the observed CMB map is consistent with Gaussian prediction. In this work we devise a new method where we use Rao's statistic based on sample arc-lengths to examine the validity of the hypothesis that the temperature field of the CMB is consistent with a Gaussian random field, by comprehensively testing correlations within a given $\ell$ mode and between nearby $\ell$ modes phases. This circular statistic is ordered and non-parametric. We performed our analysis on the scales limited by spherical harmonic modes $\le$ 128, to restrict ourselves to signal dominated region. To find the correlated sets of phases, we calculate the statistic for the data and 10000 Monte Carlo simulated random sets of phases and used 0.01 and 0.05 $\alpha$ levels to distinguish between statistically significant and highly significant detections. We apply our method on Planck satellite mission's final released CMB temperature anisotropy maps- COMMANDER, SMICA, NILC, and SEVEM along with WMAP 9 year released ILC map. We report that phases corresponding to some of the modes are non-uniform in these maps. We also report that most of the mode pairs are uncorrelated in each map, but a few are found to be correlated, which are different pairs in different maps. The detection of non-uniformity and correlations in the phases indicates presence of non Gaussian signals in the foreground minimized CMB maps.

Residual error in calibration coefficients corresponding to observed CMB maps is an important issue while estimating a pure CMB signal. A component separation method, if these errors in the input foreground contaminated CMB maps are not properly taken into account, may lead to bias in the cleaned CMB map and estimated CMB angular power spectrum. But the inability to exactly determine the calibration coefficients corresponding to each observed CMB map from any CMB experiment makes it very difficult to incorporate their exact and actual values in a component separation analysis. Hence the effect of any random and residual calibration error on the cleaned CMB map and its angular power spectrum of a component separation problem can only be understood by performing detailed Monte Carlo simulations. In this paper, we investigate the impact of using input foreground contaminated CMB maps with random calibration errors on posterior density of cleaned CMB map and theoretical CMB angular power spectrum over large angular scales of the sky following the Gibbs ILC method proposed by \cite{Sudevan:2018qyj}. By performing detailed Monte Carlo simulations of WMAP and Planck temperature anisotropy observations with calibration errors compatible with them we show that the best-fit map corresponding to posterior maximum is minimally biased in Gibbs ILC method by a CMB normalization bias and residual foreground bias. The bias in best-fit CMB angular power spectrum with respect to the case where no calibration error is present are $\sim 28 \mu K^2$ and $-4.7 \mu K^2$ respectively between $2 \le \ell \le 15$ and $16 \le \ell \le 32$. The calibration error induced error in best-fit power spectrum causes an overall $6\%$ increase of the net error when added in quadrature with the cosmic variance induced error.

The inflationary epoch and the late time acceleration of the expansion rate of universe can be explained by assuming a gravitationally coupled scalar field. In this article, we propose a new method of finding exact solutions in the background of flat Friedmann-Robertson-Walker (FRW) cosmological models by considering both scalar field and matter. We consider different forms of scalar field potential as a function of scale factor. Our method provides analytical expressions for equation of state parameter of scalar field, deceleration parameter and Hubble parameter. This method can be applied to various other forms of scalar field potential, to the early radiation dominated epoch and very early scalar field dominated inflationary dynamics. Since the method produces exact analytical expression for $H(a)$ (i.e., H(z) as well) one can use this result to constrain the model parameters using Hubble parameter data at different redshifts. As an extension of the method, we also consider the inverse problem of reconstructing scalar field potential energy by assuming any general analytical expression of scalar field equation of state parameter as a function of scale factor.

In this work, we formalize a new technique to investigate joint posterior density of Cosmic Microwave Background (CMB) signal and its theoretical angular power spectrum given the observed data, using the global internal-linear-combination (ILC) method first proposed by Sudevan & Saha (2017). We implement the method on low resolution CMB maps observed by WMAP and Planck satellite missions, using Gibbs sampling, assuming that the detector noise is negligible on large angular scales of the sky. The main products of our analysis are best fit CMB cleaned map and its theoretical angular power spectrum along with their error estimates. We validate the methodology by performing Monte Carlo simulations that includes realistic foreground models and noise levels consistent with WMAP and Planck observations. Our method has an unique advantage that the posterior density is obtained without any need to explicitly model foreground components. Secondly, the power spectrum results with the error estimates can be directly used for cosmological parameter estimations.

We propose a new internal linear combination (ILC) method in the pixel space, applicable on large angular scales of the sky, to estimate a foreground minimized Cosmic Microwave Background (CMB) temperature anisotropy map by incorporating prior knowledge about the theoretical CMB covariance matrix. Usual ILC method in pixel space, on the contrary, does not use any information about the underlying CMB covariance matrix. The new approach complements the usual pixel space ILC technique specifically at low multipole region, using global information available from theoretical CMB covariance matrix as well as from the data. Since we apply our method over the large scale on the sky containing low multipoles we perform foreground minimization globally. We apply our methods on low resolution Planck and WMAP foreground contaminated CMB maps and validate the methodology by performing detailed Monte-Carlo simulations. Our cleaned CMB map and its power spectrum have significantly less error than those obtained following usual ILC technique at low resolution that does not use CMB covariance information. Another very important advantage of our method is that the cleaned power spectrum does not have any negative bias at the low multipoles because of effective suppression of CMB-foreground chance correlations on large angular scales of the sky. Our cleaned CMB map and its power spectrum match well with those estimated by other research groups.

In this paper we estimate diffuse foreground minimized Cosmic Microwave Background (CMB) Stokes Q and U polarization maps based upon the fundamental concept of Gaussian nature of CMB and strong non-Gaussian nature of astrophysical polarized foregrounds using WMAP nine year published polarization maps. We excise regions of the sky that define position of the known point sources, regions that are strongly contaminated by either the detector noise or by the diffuse foregrounds or both, and then perform foreground minimizations over the surviving sky regions that constitute approximately $50\%$ of the full sky area. We critically evaluate performance of foreground minimizations in several ways and show that our foreground minimization method removes significant foregrounds from input maps. The cleaned Stokes \{Q, U\} polarization maps have less EE and BB power from relevant sky region compared to WMAP foreground-reduced Stokes \{Q, U\} polarization maps at different multipole ranges. We validate our methodology by performing detailed Monte Carlo simulations. The main driving machinery of our method is an internal-linear-combination (ILC) approach, however, unlike simple variance minimization performed in usual ILC method, the core of the method being dependent on the theoretically well motivated concept of Gaussianity of CMB polarization a direct connection is established between observations and models of inflation. Additionally, the method, like the usual ILC method, is independent on modeling uncertainties of polarized foregrounds. It will be useful to apply our method in future generation low-noise CMB polarization experiments.

We report an improved technique for diffuse foreground minimization from Cosmic Microwave Background (CMB) maps using a new multi-phase iterative internal-linear-combination (ILC) approach in harmonic space. The new procedure consists of two phases. In phase 1, a diffuse foreground cleaned map is obtained by performing a usual ILC operation in the harmonic space in a single iteration over the desired portion of the sky. In phase 2, we obtain the final foreground cleaned map using an iterative ILC approach also in the harmonic space, however, now, during each iteration of foreground minimization, some of the regions of the sky that are not being cleaned in the current iteration, are replaced by the corresponding cleaned portions of the phase 1 cleaned map. The new ILC method nullifies a foreground leakage signal that is otherwise inevitably present in the old and usual harmonic space iterative ILC method. The new method is flexible to handle input frequency maps, irrespective of whether or not they initially have the same instrumental and pixel resolution, by bringing them to a common and maximum possible beam and pixel resolution at the beginning of the analysis. This dramatically reduces data redundancy and hence memory usage and computational cost. During the ILC weight calculation it avoids any need to deconvolve partial sky spherical harmonic coefficients by the beam and pixel window functions, which in strict mathematical sense, is not well-defined for azimuthally symmetric window functions. Using WMAP 9-year and Planck-2015 published frequency maps we obtain a pair of foreground cleaned CMB maps and CMB angular power spectrum. Our power spectrum match well with Planck-2015 results, with some difference. Finally, we show that the weights for ILC foreground minimization have an intrinsic characteristic that it tends to produce a statistically isotropic CMB map as well.

In this paper, we implement a perturbative approach, first proposed by Bouchet & Gispert (1999), to estimate variation of spectral index of galactic polarized synchrotron emission, using linear combination of simulated Stokes Q polarization maps of selected frequency bands from WMAP and Planck observations on a region of sky dominated by the synchrotron Stokes Q signal. We find that, a first order perturbative analysis recovers input spectral index map well. Along with the spectral index variation map our method provides a fixed reference index, \hat \beta_{0s}, over the sky portion being analyzed. Using Monte Carlo simulations we find that, <\hat \beta_{0s}> = -2.84 \pm 0.01, which matches very closely with position of a peak at \beta_s(p) = -2.85, of empirical probability density function of input synchrotron indices, obtained from the same sky region. For thermal dust, mean recovered spectral index, <\hat \beta_d> = 2.00 \pm 0.004, from simulations, matches very well with spatially fixed input thermal dust spectral index \beta_d = 2.00. As accompanying results of the method we also reconstruct CMB, thermal dust and a synchrotron template component with fixed spectral indices over the {\it entire} sky region. We use full pixel-pixel noise covariance matrices of all frequency bands, estimated from the sky region being analyzed, to obtain reference spectral indices for synchrotron and thermal dust, spectral index variation map, CMB map, thermal dust and synchrotron template components. The perturbative technique as implemented in this work has the interesting property that it can build a model to describe the data with an arbitrary but enough degree of accuracy (and precession) as allowed by the data. We argue that, our method of reference spectral index determination, CMB map, thermal dust and synchrotron template component reconstruction is a maximum likelihood method.

Cyanoformaldehyde (HCOCN) molecule has recently been suspected towards the Sagittarius B2(N) by the Green Bank telescope, though a confirmation of this observation has not yet been made. In and around a star forming region, this molecule could be formed by the exothermic reaction between two abundant interstellar species, H2CO and CN. Till date, the reaction rate coefficient for the formation of this molecule is unknown. Educated guesses were used to explain the abundance of this molecule by chemical modeling. In this paper, we carried out quantum chemical calculations to find out empirical rate coefficients for the formation of HCOCN and different chemical properties during the formation of HCOCN molecules. Though HCOCN is stable against unimolecular decomposition, this gas phase molecule could be destroyed by many other means, like: ion-molecular reactions or by the effect of cosmic rays. Ion-molecular reaction rates are computed by using the capture theories. We have also included the obtained rate coefficients into our large gas-grain chemical network to study the chemical evolution of these species in various interstellar conditions. Formation of one of the isotopologue(DCOCN) of HCOCN is also studied. Our study predicts the possibility of finding HCOCN and DCOCN in the ice phase with a reasonably high abundance. In order to detect HCOCN or DCOCN in various interstellar environments, it is necessary to know the spectroscopic properties of these molecules. To this effect, we carried out quantum chemical calculations to find out different spectral parameters of HCOCN for the transition in electronic, infrared and rotational modes. We clearly show how the isotopic substitution (DCOCN) plays a part in the vibrational progressions of HCOCN.

We analyse WMAP 7-year temperature data, jointly modeling the cosmic microwave background (CMB) and Galactic foreground emission. We use the Commander code based on Gibbs sampling. Thus, from the WMAP7 data, we derive simultaneously the CMB and Galactic components on scales larger than 1deg with sensitivity improved relative to previous work. We conduct a detailed study of the low-frequency foreground with particular focus on the "microwave haze" emission around the Galactic center. We demonstrate improved performance in quantifying the diffuse galactic emission when Haslam 408MHz data are included together with WMAP7, and the spinning and thermal dust emission is modeled jointly. We also address the question of whether the hypothetical galactic haze can be explained by a spatial variation of the synchrotron spectral index. The excess of emission around the Galactic center appears stable with respect to variations of the foreground model that we study. Our results demonstrate that the new galactic foreground component - the microwave haze - is indeed present.

We describe the processing of the 336 billion raw data samples from the High Frequency Instrument (HFI) which we performed to produce six temperature maps from the first 295 days of Planck-HFI survey data. These maps provide an accurate rendition of the sky emission at 100, 143, 217, 353, 545 and 857 GHz with an angular resolution ranging from 9.9 to 4.4^2. The white noise level is around 1.5 {\mu}K degree or less in the 3 main CMB channels (100--217GHz). The photometric accuracy is better than 2% at frequencies between 100 and 353 GHz and around 7% at the two highest frequencies. The maps created by the HFI Data Processing Centre reach our goals in terms of sensitivity, resolution, and photometric accuracy. They are already sufficiently accurate and well-characterised to allow scientific analyses which are presented in an accompanying series of early papers. At this stage, HFI data appears to be of high quality and we expect that with further refinements of the data processing we should be able to achieve, or exceed, the science goals of the Planck project.

In this paper we present a new method to estimate a foreground cleaned Cosmic Microwave Background (CMB) map at a resolution of $1^\circ$ by minimizing the non-Gaussian properties of the cleaned map which arise dominantly due to diffuse foreground emission components from the Milky Way. We employ simple kurtosis statistic as the measure of non-Gaussian properties and perform a linear combination of 5 frequency maps provided by Wilkinson Microwave Anisotropy Probe (WMAP) in its 7 year data release in such a way that the cleaned map has a minimum kurtosis which leads to a non-Gaussianity minimized, foreground cleaned CMB map. We validate the method by performing Monte-Carlo simulations. To minimize any residual foreground contamination from the cleaned map we flag out region near the galactic plane based upon results from simulations. Outside the masked region our new estimate of CMB map matches well with the WMAP's ILC map. A simple pseudo-$C_l$ based CMB TT power spectrum derived from the non-gaussianity minimized map reproduces the earlier results of WMAP's power spectrum. {\it An important advantage of the method is that it does not introduce any negative bias in angular power spectrum in low multipole regime, unlike usual ILC method.} Comparing our results with the previously published results we argue that CMB results are robust with respect to specific foreground removal algorithms employed.

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.

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.

Over the last decade, measurements of the CMB anisotropy has spearheaded the remarkable transition of cosmology into a precision science. However, addressing the systematic effects in the increasingly sensitive, high resolution, `full' sky measurements from different CMB experiments pose a stiff challenge. The analysis techniques must not only be computationally fast to contend with the huge size of the data, but, the higher sensitivity also limits the simplifying assumptions which can then be invoked to achieve the desired speed without compromising the final precision goals. While maximum likelihood is desirable, the enormous computational cost makes the suboptimal method of power spectrum estimation using Pseudo-C_l unavoidable for high resolution data. We provide a (semi)analytic framework to estimate bias in the power spectrum due to the effect of beam non-circularity and non-uniform sky coverage including incomplete/masked sky maps and scan strategy. The approach is perturbative in the distortion of the beam from non-circularity, allowing for rapid computations when the beam is mildly non-circular. We advocate that it is computationally advantageous to employ `soft' azimuthally apodized masks whose spherical harmonic transform die down fast with m. We numerically implement our method for non-rotating beams. We present preliminary estimates of the computational cost to evaluate the bias for the upcoming CMB anisotropy probes l_max~3000, with angular resolution comparable to the Planck surveyor mission. We further show that this implementation and estimate is applicable for rotating beams on equal declination scans and possibly can be extended to simple approximations to other scan strategies.

Following some of the recent articles on hole super-conductivity and related phenomena by Hirsch \cite{H1,H2,H3}, a simple model is proposed to explain the observed low surface magnetic field of the expected quark stars. It is argued that the diamagnetic moments of the electrons circulating in the electro-sphere induce a magnetic field, which forces the existing quark star magnetic flux density to become dilute. We have also analysed the instability of normal-superconducting interface due to excess accumulation of magnetic flux lines, assuming an extremely slow growth of superconducting phase through a first order bubble nucleation type transition.

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.

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.

In the era of high precision CMB measurements, systematic effects are beginning to limit the ability to extract subtler cosmological information. The non-circularity of the experimental beam has become progressively important as CMB experiments strive to attain higher angular resolution and sensitivity. The effect of non-circular beam on the power spectrum is important at multipoles larger than the beam-width. For recent experiments with high angular resolution, optimal methods of power spectrum estimation are computationally prohibitive and sub-optimal approaches, such as the Pseudo-Cl method, are used. We provide an analytic framework for correcting the power spectrum for the effect of beam non-circularity and non-uniform sky coverage (including incomplete/masked sky maps). The approach is perturbative in the distortion of the beam from non-circularity allowing for rapid computations when the beam is mildly non-circular. When non-circular beam effect is important, we advocate that it is computationally advantageous to employ `soft' azimuthally apodized masks whose spherical harmonic transform die down fast with m.

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.

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.

Using a relativistic version of Landau theory of Fermi liquid with $\sigma-\omega$ and $\rho$ mesons exchange, we have obtained an equation of state for dense neutron star matter in presence of strong quantizing magnetic field. It is found that in this scenario the self energies of both neutron and proton components of dense neutron star matter become complex under certain physical conditions. To be more specific, it is observed that in the exchange diagrams of $\sigma$, $\omega$ and neutral $\rho$ transfer processes and in the direct interaction diagram with $\rho_\pm$ transfer reactions, the nucleon self-energies become complex in nature.