Showing 1–48 of 48 results for author: Saha, R

Foreground removal and angular power spectrum estimation of 21 cm signal using harmonic space ILC method

Authors: Albin Joseph, Rajib Saha

Abstract: Mapping the distribution of neutral atomic hydrogen (HI) in the Universe through its 21 cm emission line provides a powerful cosmological probe to map the large-scale structures and shed light on various cosmological phenomena. The Baryon Acoustic Oscillations at low redshifts can potentially be probed by sensitive HI intensity mapping experiments and constrain the properties of dark energy. Howev… ▽ More Mapping the distribution of neutral atomic hydrogen (HI) in the Universe through its 21 cm emission line provides a powerful cosmological probe to map the large-scale structures and shed light on various cosmological phenomena. The Baryon Acoustic Oscillations at low redshifts can potentially be probed by sensitive HI intensity mapping experiments and constrain the properties of dark energy. However, the 21 cm signal detection faces formidable challenges due to the dominance of various astrophysical foregrounds, which can be several orders of magnitude stronger. Our current work introduces a novel and model-independent Internal Linear Combination (ILC) method in harmonic space using the principal components of the 21 cm signal for accurate foreground removal and power spectrum estimation. We estimate the principal components by incorporating prior knowledge of the theoretical 21 cm covariance matrix. We test our methodology by detailed simulations of radio observations, incorporating synchrotron emission, free-free radiation, extragalactic point sources, and thermal noise. We estimate the full sky 21 cm angular power spectrum after application of a mask on the full sky cleaned 21 cm signal by using the mode-mode coupling matrix. These full sky estimates of angular spectra can be directly used to measure the cosmological parameters. For the first time, we demonstrate the effectiveness of a foreground model-independent ILC method in harmonic space to reconstruct the 21 cm signal. △ Less

Submitted 4 May, 2024; originally announced May 2024.

Comments: 20 pages, 15 figures

Accurate and Unbiased Reconstruction of CMB B Mode using Deep Learning

Authors: Srikanta Pal, Sarvesh Kumar Yadav, Rajib Saha, Tarun Souradeep

Abstract: An ingeniously designed autoencoder (PrimeNet) using simulated observations of future generation ECHO satellite mission recovers CMB B mode map, angular spectrum for multipoles $\ell \lesssim 9$ and tensor to scalar ratio $r$ {\it limited only by cosmic variance down to $r= 0.0001$ and below}. We use diverse, realistically complex and detailed foreground models. PrimeNet predicts accurate results… ▽ More An ingeniously designed autoencoder (PrimeNet) using simulated observations of future generation ECHO satellite mission recovers CMB B mode map, angular spectrum for multipoles $\ell \lesssim 9$ and tensor to scalar ratio $r$ {\it limited only by cosmic variance down to $r= 0.0001$ and below}. We use diverse, realistically complex and detailed foreground models. PrimeNet predicts accurate results even when data with $r=0$ are tested which were not used in training, implying robust and efficient predictive power. The work eliminates a major bottleneck of weak CMB B mode reconstruction and takes a leap forward for understanding fundamental physics of the primordial Universe. △ Less

Submitted 28 April, 2024; originally announced April 2024.

The Massive and Distant Clusters of WISE Survey 2: Equatorial First Data Release

Authors: Khunanon Thongkham, Anthony H. Gonzalez, Mark Brodwin, Ariane Trudeau, Ripon Saha, Peter Eisenhardt, S. A. Stanford, Emily Moravec, Thomas Connor, Daniel Stern

Abstract: The Massive and Distant Clusters of WISE Survey 2 (MaDCoWS2) is a new survey designed as the successor of the original MaDCoWS survey. MaDCoWS2 improves upon its predecessor by using deeper optical and infrared data and a more powerful detection algorithm (PZWav). As input to the search, we use grz photometry from DECaLS in combination with W1 and W2 photometry from the CatWISE2020 catalog to deri… ▽ More The Massive and Distant Clusters of WISE Survey 2 (MaDCoWS2) is a new survey designed as the successor of the original MaDCoWS survey. MaDCoWS2 improves upon its predecessor by using deeper optical and infrared data and a more powerful detection algorithm (PZWav). As input to the search, we use grz photometry from DECaLS in combination with W1 and W2 photometry from the CatWISE2020 catalog to derive the photometric redshifts with full redshift probability distribution functions for WISE-selected galaxies. Cluster candidates are then detected using the PZWav algorithm to find three-dimensional galaxy overdensities from the sky positions and photometric redshifts. This paper provides the first MaDCoWS2 data release, covering 1461 (1838 without masking) deg^2 centered on the Hyper-SuprimeCam Subaru Strategic Program equatorial fields. Within this region, we derive a catalog of 22,970 galaxy cluster candidates detected at S/N>5. These clusters span the redshift range 0.1<z<2, including 1312 candidates at z>1.5. We compare MaDCoWS2 to six existing catalogs in the area. We rediscover 60%-92% of the clusters in these surveys at S/N>5. The medians of the absolute redshift offset are <0.02 relative to these surveys, while the standard deviations are less than 0.06. The median offsets between the detection position from MaDCoWS2 and other surveys are less than 0.25 Mpc. We quantify the relation between S/N and gas mass, total mass, luminosity, and richness from other surveys using a redshift-dependent power law relation. We find that the S/N-richness relation exhibits the lowest scatter. △ Less

Submitted 11 April, 2024; v1 submitted 25 March, 2024; originally announced March 2024.

Comments: 27 pages, 7 figures. Typo corrected. Accepted for publication in ApJ

ParamANN: A Neural Network to Estimate Cosmological Parameters for $Λ$CDM Universe using Hubble Measurements

Authors: Srikanta Pal, Rajib Saha

Abstract: In this article, we employ a machine learning (ML) approach for the estimations of four fundamental parameters, namely, the Hubble constant ($H_0$), matter ($Ω_{0m}$), curvature ($Ω_{0k}$) and vacuum ($Ω_{0Λ}$) densities of non-flat $Λ$CDM model. We use $53$ Hubble parameter values measured by differential ages (DA) and baryon acoustic oscillations (BAO) techniques in the redshift interval… ▽ More In this article, we employ a machine learning (ML) approach for the estimations of four fundamental parameters, namely, the Hubble constant ($H_0$), matter ($Ω_{0m}$), curvature ($Ω_{0k}$) and vacuum ($Ω_{0Λ}$) densities of non-flat $Λ$CDM model. We use $53$ Hubble parameter values measured by differential ages (DA) and baryon acoustic oscillations (BAO) techniques in the redshift interval $0.07 \leq z \leq 2.36$. We create an artificial neural network (called ParamANN) and train it with simulated values of $H(z)$ using various sets of $H_0$, $Ω_{0m}$, $Ω_{0k}$, $Ω_{0Λ}$ parameters chosen from different and sufficiently wide prior intervals. We use a correlated noise model in the analysis. We demonstrate accurate validation and prediction by ParamANN. ParamANN provides an excellent cross-check for the validity of the $Λ$CDM model and alleviates the Hubble tension problem which has been reported earlier in the literature. We obtain $H_0 = 66.11 \pm 2.59$ $\rm{kmMpc^{-1}sec^{-1}}$, $Ω_{0m} = 0.3359 \pm 0.0814$, $Ω_{0k} = 0.0237 \pm 0.1248$ and $Ω_{0Λ} = 0.6405 \pm 0.0861$ by using the trained network. These parameter values agree very well with the results of global CMB observations of Planck collaboration. △ Less

Submitted 26 September, 2023; originally announced September 2023.

Comments: 21 pages, 5 figures, 3 tables

Detection of Dipole Modulation in CMB Temperature Anisotropy Maps from WMAP and Planck using Artificial Intelligence

Authors: Md Ishaque Khan, Rajib Saha

Abstract: Breakdown of rotational invariance of the primordial power spectrum manifests in the statistical anisotropy of the observed Cosmic Microwave Background (CMB) radiation. Hemispherical power asymmetry in the CMB may be caused due to a dipolar modulation, indicating the presence of a preferred direction. Appropriately re-scaled local variance maps of the CMB temperature anisotropy data effectively en… ▽ More Breakdown of rotational invariance of the primordial power spectrum manifests in the statistical anisotropy of the observed Cosmic Microwave Background (CMB) radiation. Hemispherical power asymmetry in the CMB may be caused due to a dipolar modulation, indicating the presence of a preferred direction. Appropriately re-scaled local variance maps of the CMB temperature anisotropy data effectively encapsulate this dipolar pattern. As a first-of-its-kind method, we train Artificial Neural Networks (ANNs) with such local variances as input features to distinguish statistically isotropic CMB maps from dipole modulated ones. Our trained ANNs are able to predict components of the amplitude times the unit vector of the preferred direction for mixed sets of modulated and unmodulated maps, with goodness of fit ($R^2$) scores $>0.97$ for full sky, and $>0.96$ for partial sky coverage. On all observed foreground-cleaned CMB maps, the ANNs detect the dipolar modulation signal with overall consistent values of amplitudes and directions. This detection is significant at $97.21\%-99.38\%$ C.L. for all full sky maps, and at $98.34\%-100\%$ C.L. for all partial sky maps. Robustness of the signal holds across full and partial skies, various foreground cleaning methods, inpainting algorithms, instruments and all the different periods of observation for Planck and WMAP satellites. The significant and robust detection of the signal, in addition to the consistency of values of amplitude and directions, as found independent of any pre-existing methods, further mitigates the criticisms of look-elsewhere effects and a posteriori inferences for the preferred dipole direction in the CMB. △ Less

Submitted 22 May, 2023; v1 submitted 8 December, 2022; originally announced December 2022.

Comments: 21 pages, 16 figures

Journal ref: 2023 ApJ 947 47

Reconstruction of full sky CMB $\bf{E}$ and $\bf{B}$ modes spectra removing $\bf{E}$-to-$\bf{B}$ leakage from partial sky using deep learning

Authors: Srikanta Pal, Rajib Saha

Abstract: Incomplete sky analysis of cosmic microwave background (CMB) polarization spectra poses a major problem of leakage between $E$- and $B$-modes. We present a machine learning approach to remove this $E$-to-$B$ leakage using a convolutional neural network (CNN) in presence of detector noise. The CNN predicts the full sky $E$- and $B$-modes spectra for multipoles $2 \leq \ell \leq 384$ from the partia… ▽ More Incomplete sky analysis of cosmic microwave background (CMB) polarization spectra poses a major problem of leakage between $E$- and $B$-modes. We present a machine learning approach to remove this $E$-to-$B$ leakage using a convolutional neural network (CNN) in presence of detector noise. The CNN predicts the full sky $E$- and $B$-modes spectra for multipoles $2 \leq \ell \leq 384$ from the partial sky spectra for $N_{\rm{side}} = 256$. We use tensor-to-scalar ratio $r=0.001$ to simulate the CMB polarization maps. We train our CNN using $10^5$ full sky target spectra and an equal number of noise contaminated partial sky spectra obtained from the simulated maps. The CNN works well for two masks covering the sky area of $\sim 80\%$ and $\sim 10\%$ respectively after training separately for each mask. For the assumed theoretical $E$- and $B$-modes spectra, predicted full sky $E$- and $B$-modes spectra agree well with the corresponding target spectra and their means agree with theoretical spectra. The CNN preserves the cosmic variances at each multipole, effectively removes correlations of the partial sky $E$- and $B$-modes spectra, and retains the entire statistical properties of the targets avoiding the problem of so-called $E$-to-$B$ leakage for the chosen theoretical model. △ Less

Submitted 26 September, 2023; v1 submitted 16 November, 2022; originally announced November 2022.

Comments: 24 pages, 13 figures

Journal ref: J. Astrophys. Astron. 44 (2023) 84

A Foreground Model Independent Bayesian CMB Temperature and Polarization Signal Reconstruction and Cosmological Parameter Estimation over Large Angular Scales

Authors: Albin Joseph, Ujjal Purkayastha, Rajib Saha

Abstract: Recent CMB observations have resulted in very precise observational data. A robust and reliable CMB reconstruction technique can lead to efficient estimation of the cosmological parameters. We demonstrate the performance of our methodology using simulated temperature and polarization observations using cosmic variance limited future generation PRISM satellite mission. We generate samples from the… ▽ More Recent CMB observations have resulted in very precise observational data. A robust and reliable CMB reconstruction technique can lead to efficient estimation of the cosmological parameters. We demonstrate the performance of our methodology using simulated temperature and polarization observations using cosmic variance limited future generation PRISM satellite mission. We generate samples from the joint distribution by implementing the CMB inverse covariance weighted internal-linear-combination (ILC) with the Gibbs sampling technique. We use the Python Sky Model (PySM), d4f1s1 to generate the realistic foreground templates. The synchrotron emission is parametrized by a spatially varying spectral index, whereas the thermal dust emission is described as a two-component dust model. We estimate the marginalized densities of CMB signal and theoretical angular power spectrum utilizing the samples from the entire posterior distribution. The best-fit cleaned CMB map and the corresponding angular power spectrum are consistent with the CMB realization and the sky angular power spectrum, implying an efficient foreground minimized reconstruction. The likelihood function estimated by making use of the Blackwell-Rao estimator is used for the estimation of the cosmological parameters. Our methodology can estimate the tensor to scalar ratio $r\ge 0.0075$ for the chosen foreground models and the instrumental noise levels. Our current work demonstrates an analysis pipeline starting from the reliable estimation of CMB signal and its angular power spectrum to the case of cosmological parameter estimation using the foreground model independent Gibbs-ILC method. △ Less

Submitted 30 January, 2023; v1 submitted 15 September, 2022; originally announced September 2022.

Forecast Analysis on Interacting Dark Energy Models from Future Generation PICO and DESI Missions

Authors: Albin Joseph, Rajib Saha

Abstract: The next-generation CMB satellite missions are expected to provide robust constraints on a wide range of cosmological parameters with unprecedented precision. But these constraints on the parameters could weaken if we do not attribute dark energy to a cosmological constant. The cosmological models involving interaction between dark energy and dark matter can give rise to comparable energy densitie… ▽ More The next-generation CMB satellite missions are expected to provide robust constraints on a wide range of cosmological parameters with unprecedented precision. But these constraints on the parameters could weaken if we do not attribute dark energy to a cosmological constant. The cosmological models involving interaction between dark energy and dark matter can give rise to comparable energy densities at the present epoch, thereby alleviating the so-called cosmic coincidence problem. In the present paper, we perform a forecast analysis to test the ability of the future generation high-sensitive Cosmic Microwave Background (CMB), and Baryon Acoustic Oscillation (BAO) experiments to constrain phenomenological interacting dark energy models. We consider cosmic variance limited future CMB polarization experiment PICO along with BAO information from the DESI experiment to constrain the parameters of the interacting dark sector. Based on the stability of the cosmological perturbations, we consider two possibilities for the interaction scenario. We investigate the impact of both the coupling constant and the equation of state parameter of dark energy on the CMB temperature power spectrum, matter power spectrum, and $fσ_8$. We have used simulated temperature and polarization data from PICO within the multipole ranges ($\ell = 2 - 4000$), and as expected, we do see PICO alone produces better constraints than Planck on the $Λ$CDM parameters. With the integration of the PICO and DESI missions, we observe a significant improvement in the constraints on several cosmological parameters, especially the equation of state parameter of dark energy. However, we note that additional data is required to constrain a small positive coupling constant. △ Less

Submitted 5 December, 2022; v1 submitted 15 September, 2022; originally announced September 2022.

On Direct Estimation of Density Parameters Alleviating Hubble Tension for $Λ$CDM Universe using Hubble Measurements

Authors: Srikanta Pal, Rajib Saha

Abstract: The set of cosmological density parameters ($Ω_{0m}h_{0}^{2}$, $Ω_{0k}h_{0}^{2}$, $Ω_{0Λ}h_{0}^{2}$) and Hubble constant ($h_{0}$) are useful for fundamental understanding of the universe from many perspectives. In this article, we propose a new procedure to estimate these parameters for $Λ$ cold dark matter ($Λ$CDM) universe in the Friedmann-Robertson-Walker (FRW) background. We generalize the tw… ▽ More The set of cosmological density parameters ($Ω_{0m}h_{0}^{2}$, $Ω_{0k}h_{0}^{2}$, $Ω_{0Λ}h_{0}^{2}$) and Hubble constant ($h_{0}$) are useful for fundamental understanding of the universe from many perspectives. In this article, we propose a new procedure to estimate these parameters for $Λ$ cold dark matter ($Λ$CDM) universe in the Friedmann-Robertson-Walker (FRW) background. We generalize the two-point statistic first proposed by Sahni et al. (2008) to the three point case and estimate the parameters using currently available Hubble parameter ($H(z)$) values in the redshift range $0.07 \leq z \leq 2.36$ measured by differential age (DA) and baryon acoustic oscillation (BAO) techniques. All the parameters are estimated assuming the general case of non-flat universe. Using both DA and BAO data we obtain $Ω_{0m}h_{0}^{2}=0.1485 \pm 0.0065$, $Ω_{0k}h_{0}^{2}=-0.0137 \pm 0.017$, $Ω_{0Λ}h_{0}^{2}=0.3126 \pm 0.0145$ and $h_{0}=0.6689 \pm 0.0021$. These results are in satisfactory agreement with the Planck results. An important advantage of our method is that to estimate value of any one of the independent cosmological parameters one does not need to use values for rest of them. Each parameter is obtained solely from the measured values of Hubble parameters at different redshifts without any need to use values of other parameters. Such a method is expected to be less susceptible to the undesired effects of degeneracy issues between cosmological parameters during their estimations. Moreover, there is no requirement of assuming spatial flatness in our method. △ Less

Submitted 26 September, 2023; v1 submitted 14 April, 2022; originally announced April 2022.

Comments: 15 pages, 6 figures

Estimation of Full Sky Power Spectrum between Intermediate to Large Angular Scales from Partial Sky CMB Anisotropies using Artificial Neural Network

Authors: Srikanta Pal, Pallav Chanda, Rajib Saha

Abstract: Reliable extraction of cosmological information from observed cosmic microwave background (CMB) maps may require removal of strongly foreground contaminated regions from the analysis. In this article, we employ an artificial neural network (ANN) to predict the full sky CMB angular power spectrum between intermediate to large angular scales from the partial sky spectrum obtained from masked CMB tem… ▽ More Reliable extraction of cosmological information from observed cosmic microwave background (CMB) maps may require removal of strongly foreground contaminated regions from the analysis. In this article, we employ an artificial neural network (ANN) to predict the full sky CMB angular power spectrum between intermediate to large angular scales from the partial sky spectrum obtained from masked CMB temperature anisotropy map. We use a simple ANN architecture with one hidden layer containing $895$ neurons. Using $1.2 \times 10^{5}$ training samples of full sky and corresponding partial sky CMB angular power spectra at Healpix pixel resolution parameter $N_{side} = 256$, we show that predicted spectrum by our ANN agrees well with the target spectrum at each realization for the multipole range $2 \leq l \leq 512$. The predicted spectra are statistically unbiased and they preserve the cosmic variance accurately. Statistically, the differences between the mean predicted and underlying theoretical spectra are within approximately $3σ$. Moreover, the probability densities obtained from predicted angular power spectra agree very well with those obtained from `actual' full sky CMB angular power spectra for each multipole. Interestingly, our work shows that the significant correlations in input cut-sky spectra, due to mode-mode coupling introduced on the partial sky, are effectively removed since the ANN learns the hidden pattern between the partial sky and full sky spectra preserving the entire statistical properties. The excellent agreement of statistical properties between the predicted and the ground-truth demonstrates the importance of using artificial intelligence systems in cosmological analysis more widely. △ Less

Submitted 23 February, 2023; v1 submitted 26 March, 2022; originally announced March 2022.

Comments: 17 pages, 15 figures

Journal ref: ApJ 945 (2023) 77

The Massive and Distant Clusters of WISE Survey XI: Stellar Mass Fractions and Luminosity Functions of MaDCoWS Clusters at $z \sim 1$

Authors: Bandon Decker, Mark Brodwin, Ripon Saha, Thomas Connor, Peter R. M. Eisenhardt, Anthony H. Gonzalez, Emily Moravec, Mustafa Muhibullah, S. Adam Stanford, Daniel Stern, Khunanon Thongkham, Dominika Wylezalek, Simon R. Dicker, Brian Mason, Tony Mroczkowski, Charles E. Romero, Florian Ruppin

Abstract: We present stellar mass fractions and composite luminosity functions (LFs) for a sample of \Ncl\ clusters from the Massive and Distant Clusters of WISE Survey (MaDCoWS) at a redshift range of $0.951 \leq z \leq 1.43$. Using SED fitting of optical and deep mid-infrared photometry, we establish the membership of objects along the lines-of-sight to these clusters and calculate the stellar masses of m… ▽ More We present stellar mass fractions and composite luminosity functions (LFs) for a sample of \Ncl\ clusters from the Massive and Distant Clusters of WISE Survey (MaDCoWS) at a redshift range of $0.951 \leq z \leq 1.43$. Using SED fitting of optical and deep mid-infrared photometry, we establish the membership of objects along the lines-of-sight to these clusters and calculate the stellar masses of member galaxies. We find stellar mass fractions for these clusters largely consistent with previous works, including appearing to display a negative correlation with total cluster mass. We measure a composite $3.6~\mathrm{μm}$ LF down to $m^*+2.5$ for all 12 clusters. Fitting a Schechter function to the LF, we find a characteristic $3.6~\mathrm{μm}$ magnitude of $m^*=19.83\pm0.12$ and faint-end slope of $α=-0.81\pm0.10$ for the full sample at a mean redshift of $\bar{z} = 1.18$. We also divide the clusters into high- and low-redshift bins at $\bar{z}=1.29$ and $\bar{z}=1.06$ respectively and measure a composite LF for each bin. We see a small, but statistically significant evolution in $m^*$ and $α$ -- consistent with passive evolution -- when we study the joint fit to the two parameters, which is probing the evolution of faint cluster galaxies at $z\sim1$. This highlights the importance of deep IR data in studying the evolution of cluster galaxy populations at high-redshift. △ Less

Submitted 22 December, 2021; originally announced December 2021.

Comments: 13 pages, 11 figures, Submitted to ApJ

Isotropy statistics of CMB hot and cold spots

Authors: Md Ishaque Khan, Rajib Saha

Abstract: Statistical Isotropy of the Cosmic Microwave Background (CMB) radiation has been studied and debated extensively in recent years. Under this assumption, the hot spots and cold spots of the CMB are expected to be uniformly distributed over a 2-sphere. We use the orientation matrix, first proposed by Watson (1965) and Scheidegger (1965) and associated shape and strength parameters (Woodcock, 1977) t… ▽ More Statistical Isotropy of the Cosmic Microwave Background (CMB) radiation has been studied and debated extensively in recent years. Under this assumption, the hot spots and cold spots of the CMB are expected to be uniformly distributed over a 2-sphere. We use the orientation matrix, first proposed by Watson (1965) and Scheidegger (1965) and associated shape and strength parameters (Woodcock, 1977) to analyse whether the hot and cold spots of the observed CMB temperature anisotropy field are uniformly placed. We demonstrate the usefulness of our estimators by using simulated toy models containing non-uniform data. We apply our method on several foreground minimized CMB maps observed by WMAP and Planck over large angular scales. The shape and strength parameters constrain geometric features of possible deviations from uniformity (isotropy) and the power of the anomalous signal. We find that distributions of hot or cold spots in cleaned maps show no unusual signature of clustering or girdling. Instead, we notice a strikingly uniform distribution of hot spots over the full sky. The signal remains robust with respect to the four cleaned maps used and presence or absence of the non-Gaussian cold spot (NGCS). On the partial sky with WMAP KQ75 and Planck U73 masks we find anomalously low strength of non-uniformity for cold spots which is found to be robust with respect to various cleaning methods, masks applied, instruments, frequencies, and the presence or absence of the NGCS. Interestingly we find that the signal of anomalously weak non-uniformity could be due to contributions from the quadrupole and octupole and may be related with the low CMB temperature variance anomaly. △ Less

Submitted 11 November, 2022; v1 submitted 10 November, 2021; originally announced November 2021.

Comments: 25 pages, 11 figures

Journal ref: JCAP06 (2022) 006

A Partial-Sky Gibbs ILC Approach for the Estimation of CMB Posterior over Large Angular Scales of the Sky

Authors: Vipin Sudevan, Ujjal Purkayastha, Rajib Saha

Abstract: In this article we present a formalism to incorporate the partial-sky maps to the Gibbs ILC algorithm to estimate the joint posterior density of the Cosmic Microwave Background (CMB) signal and the theoretical CMB angular power spectrum given the observed CMB maps. In order to generate the partial-sky maps we mask all the observed CMB maps provided by WMAP and Planck satellite mission using a Gaus… ▽ More In this article we present a formalism to incorporate the partial-sky maps to the Gibbs ILC algorithm to estimate the joint posterior density of the Cosmic Microwave Background (CMB) signal and the theoretical CMB angular power spectrum given the observed CMB maps. In order to generate the partial-sky maps we mask all the observed CMB maps provided by WMAP and Planck satellite mission using a Gaussian smoothed mask formed on the basis of thermal dust emissions in Planck 353 GHz map. The central galactic region from all the input maps is removed after the application of the smoothed mask. While implementing the Gibbs ILC method on the partial-sky maps, we convert the partial-sky cleaned angular power spectrum to the full-sky angular power spectrum using the mode-mode coupling matrix estimated from the smoothed mask. The main products of our analysis are partial-sky cleaned best-fit CMB map and an estimate of the underlying full-sky theoretical CMB angular power spectrum along with their error estimates. We validate the methodology by performing detailed Monte Carlo simulations after using realistic models of foregrounds and detector noise consistent with the WMAP and Planck frequency channels used in our analysis. We can estimate the posterior density and full-sky theoretical CMB angular power spectrum, without any need to explicitly model the foreground components, from partial-sky maps using our method. Another important feature of this method is that the power spectrum results along with the error estimates can be directly used for cosmological parameter estimations. △ Less

Submitted 15 October, 2021; originally announced October 2021.

Dark energy with oscillatory tracking potential: Observational Constraints and Perturbative effects

Authors: Albin Joseph, Rajib Saha

Abstract: The cosmological models exhibiting tracker properties have great significance in the context of dark energy as they can reach the present value of dark energy density from a wide range of initial conditions, thereby alleviating both the fine-tuning and the cosmic coincidence problem. The $α$-attractors, which are originally discussed in the context of inflation, can exhibit the properties of dark… ▽ More The cosmological models exhibiting tracker properties have great significance in the context of dark energy as they can reach the present value of dark energy density from a wide range of initial conditions, thereby alleviating both the fine-tuning and the cosmic coincidence problem. The $α$-attractors, which are originally discussed in the context of inflation, can exhibit the properties of dark energy as they can behave like cosmological trackers at early times and show the late time behaviour of a cosmological constant. In the present paper, we study the Oscillatory Tracker Model (OTM), which belongs to the family of $α$-attractor dark energy models. Using the current observational data sets like Cosmic Microwave Background (CMB), Baryon Acoustic Oscillation (BAO) and type 1a supernova data (Pantheon compilation), we constrain the parameters of the model and estimate both the mean and best-fit values. Although the oscillatory tracker model contains a larger set of parameters than the usual LCDM model, the common set of parameters of both agree within $1\, σ$ error limits. Our observations using both high redshift and low redshift data supports Hubble parameter value $H_0 = 67.4$ Kms$^{-1}$Mpc$^{-1}$. We study the effect of the OTM on the CMB temperature and polarization power spectra, matter power spectrum and $f σ_8$. Our analysis of the CMB power spectrum and matter power spectrum suggests that the oscillatory tracker dark energy model has noticeable differences from usual LCDM predictions. Yet, in most cases, the agreement is very close. △ Less

Submitted 27 January, 2022; v1 submitted 1 October, 2021; originally announced October 2021.

Application of Genetic Algorithm to Estimate the Large Angular Scale Features of Cosmic Microwave Background

Authors: Parth Nayak, Rajib Saha

Abstract: Genetic Algorithm (GA) -- motivated by natural evolution -- is a robust method to estimate the global optimal solutions of problems involving multiple objective functions. In this article, for the first time, we apply GA to reconstruct the CMB temperature anisotropy map over large angular scales of the sky using (internal) linear combination (ILC) of the final-year WMAP and Planck satellite observ… ▽ More Genetic Algorithm (GA) -- motivated by natural evolution -- is a robust method to estimate the global optimal solutions of problems involving multiple objective functions. In this article, for the first time, we apply GA to reconstruct the CMB temperature anisotropy map over large angular scales of the sky using (internal) linear combination (ILC) of the final-year WMAP and Planck satellite observations. To avoid getting trapped into a local minimum, we implement the GA with generous diversity in the populations by selecting pairs with diverse fitness coefficients and by introducing a small but significant amount of mutation of genes. We find that the new GA-ILC method produces a clean map which agrees very well with that obtained using the exact analytical expression of weights in ILC. By performing extensive Monte Carlo simulations of the CMB reconstruction using the GA-ILC algorithm, we find that residual foregrounds in the cleaned map are minimal and tend to occupy localized regions along the galactic plane. The CMB angular power spectrum shows no indication of any bias in the entire multipole range $2 \leq \ell \leq 32$ studied in this work. The error in the CMB angular power spectrum is also minimal and given entirely by the cosmic-variance-induced error. Our results agree well with those obtained by various other reconstruction methods by different research groups. This problem-independent robust GA-ILC method provides a flexible way towards the complex and challenging task of CMB component reconstruction in cosmology. △ Less

Submitted 27 March, 2022; v1 submitted 12 February, 2021; originally announced February 2021.

Journal ref: MNRAS 510-2(2022) 2173-2185

An Unbiased Estimator of the Full-sky CMB Angular Power Spectrum at Large Scales using Neural Networks

Authors: Pallav Chanda, Rajib Saha

Abstract: Accurate estimation of the Cosmic Microwave Background (CMB) angular power spectrum is enticing due to the prospect for precision cosmology it presents. Galactic foreground emissions, however, contaminate the CMB signal and need to be subtracted reliably in order to lessen systematic errors on the CMB temperature estimates. Typically bright foregrounds in a region lead to further uncertainty in te… ▽ More Accurate estimation of the Cosmic Microwave Background (CMB) angular power spectrum is enticing due to the prospect for precision cosmology it presents. Galactic foreground emissions, however, contaminate the CMB signal and need to be subtracted reliably in order to lessen systematic errors on the CMB temperature estimates. Typically bright foregrounds in a region lead to further uncertainty in temperature estimates in the area even after some foreground removal technique is performed and hence determining the underlying full-sky angular power spectrum poses a challenge. We explore the feasibility of utilizing artificial neural networks to predict the angular power spectrum of the full sky CMB temperature maps from the observed angular power spectrum of the partial sky in which CMB temperatures in some bright foreground regions are masked. We present our analysis at large angular scales with two different masks. We produce unbiased predictions of the full-sky angular power spectrum and recover the underlying theoretical power spectrum using neural networks. Our predictions are also uncorrelated to a large extent. We further show that the multipole-space covariances of the predictions of full-sky spectra made by the ANNs are much smaller than those of the estimates obtained using the pseudo-$C_\ell$ method. △ Less

Submitted 22 September, 2021; v1 submitted 8 February, 2021; originally announced February 2021.

Comments: 10 pages, 11 figures; altered methodology, added links to references, updated analysis using latest available data, modified the write-up accordingly

Level correlations of the CMB temperature angular power spectrum

Authors: Md Ishaque Khan, Rajib Saha

Abstract: If the principle of statistical isotropy is valid, then the angular power spectrum (APS) of cosmic microwave background (CMB) radiation is uncorrelated between different multipoles. We propose a novel technique to analyse any possible correlations of the foreground cleaned CMB temperature angular power spectrum (APS) measures ($C_{\ell}$ and $\mathcal{D}_\ell=\frac{\ell(\ell+1)}{2π}C_\ell$). This… ▽ More If the principle of statistical isotropy is valid, then the angular power spectrum (APS) of cosmic microwave background (CMB) radiation is uncorrelated between different multipoles. We propose a novel technique to analyse any possible correlations of the foreground cleaned CMB temperature angular power spectrum (APS) measures ($C_{\ell}$ and $\mathcal{D}_\ell=\frac{\ell(\ell+1)}{2π}C_\ell$). This is motivated by the behaviour of level spacings between random matrix eigenvalues. The method helps distinguish uncorrelated statistically isotropic CMB APS from correlated APS, where the latter can arise due to breakdown of isotropy or presence of some residual systematics in the foreground cleaned CMB maps. Spacings of statistically isotropic CMB $C_\ell$'s and $\mathcal{D}_\ell$'s are seen to closely obey Poisson statistics and introduction of correlations changes the distribution to appropriate Wigner-Dyson statistics. For foreground cleaned CMB, we employ the average spacing of consecutive multipole APS for multipoles $\in [2,31]$. This estimator is sensitive to departures from the null hypothesis of zero correlations between the APS measures of statistically isotropic CMB. We study full sky WMAP 9 year ILC and 2018 Planck foreground cleaned maps (Commander, NILC and SMICA). Sans parity distinctions, average spacings are in good agreement with theoretical expectation. With parity distinctions, even multipoles indicate unusually low average spacings for both $C_\ell$'s (at $\geq 98.86\%$ C.L.) and $\mathcal{D}_\ell$'s (at $\geq 95.07\%$ C.L.). We use an inpainting method based on constrained Gaussian realisations and show that for the Planck $U73$ and WMAP $KQ75$ masks, all the foreground cleaned inpainted CMB maps robustly confirm the existence of such unusually low average spacings of even multipole APS. In addition, this signal is independent of the non-Gaussian cold spot. △ Less

Submitted 21 November, 2022; v1 submitted 17 January, 2021; originally announced January 2021.

Comments: 16 pages, 9 figures

Journal ref: J Astrophys Astron 43, 100 (2022)

Comparing Classification Models on Kepler Data

Authors: Rohan Saha

Abstract: Even though the original Kepler mission ended due to mechanical failures, the Kepler satellite continues to collect data. Using classification models, we can understand the features exoplanets possess and then use those features to investigate further for any more information on the candidate planet. Based on the classification model, the idea is to find out the probability of the planet under obs… ▽ More Even though the original Kepler mission ended due to mechanical failures, the Kepler satellite continues to collect data. Using classification models, we can understand the features exoplanets possess and then use those features to investigate further for any more information on the candidate planet. Based on the classification model, the idea is to find out the probability of the planet under observation being a candidate for an exoplanet or a false positive. If the model predicts that the observation is a candidate for being an exoplanet, then the further investigation can be conducted. From the model, we can narrow down the features that might explain the difference between a candidate and a false-positive which ultimately helps us to increase the efficiency of any model and fine-tune the model and ultimately the process of searching for any future exoplanets. The model comparison is supported by McNemar's test for checking significance. △ Less

Submitted 6 January, 2021; v1 submitted 6 January, 2021; originally announced January 2021.

A Bayesian ILC method for CMB B-mode posterior estimation and reconstruction of primordial gravity wave signal

Authors: Sarvesh Kumar Yadav, Rajib Saha

Abstract: The Cosmic Microwave Background (CMB) radiation B mode polarization signal contains the unique signature of primordial metric perturbations produced during the inflation. The separation of the weak CMB B-mode signal from strong foreground contamination in observed maps is a complex task, and proposed new generation low noise satellite missions compete with the weak signal level of this gravitation… ▽ More The Cosmic Microwave Background (CMB) radiation B mode polarization signal contains the unique signature of primordial metric perturbations produced during the inflation. The separation of the weak CMB B-mode signal from strong foreground contamination in observed maps is a complex task, and proposed new generation low noise satellite missions compete with the weak signal level of this gravitational background. In this article, for the first time, we employ a foreground model-independent internal linear combination (ILC) method to reconstruct the CMB B mode signal using simulated observations over large angular scales of the sky of 6 frequency bands of future generation CMB mission Probe of Inflation and Cosmic Origins (PICO). We estimate the joint CMB B mode posterior density following the interleaving Gibbs steps of B mode angular power spectrum and cleaned map samples using the ILC method. We extend and improve the earlier reported Bayesian ILC method to analyze weak CMB B mode reconstruction by introducing noise bias corrections at two stages during the ILC weight estimation. By performing $200$ Monte Carlo simulations of the Bayesian ILC method, we find that our method can reconstruct the CMB signals and the joint posterior density accurately over large angular scales of the sky. We estimate Blackwell-Rao statistics of the marginal density of CMB B mode angular power spectrum and use them to estimate the joint density of scalar to tensor ratio $r$ and a lensing power spectrum amplitude $A^{\textrm{lens}}$. Using $200$ Monte Carlo simulations of the delensing approach, we find that our method can achieve an unbiased detection of the primordial gravitational wave signal $r$ with more than 8$σ$ significance for levels of $r \geqslant 0.01$. △ Less

Submitted 30 September, 2020; originally announced September 2020.

Comments: 12 pages,13 figures,1 table

Euclid preparation: X. The Euclid photometric-redshift challenge

Authors: Euclid Collaboration, G. Desprez, S. Paltani, J. Coupon, I. Almosallam, A. Alvarez-Ayllon, V. Amaro, M. Brescia, M. Brodwin, S. Cavuoti, J. De Vicente-Albendea, S. Fotopoulou, P. W. Hatfield, W. G. Hartley, O. Ilbert, M. J. Jarvis, G. Longo, R. Saha, J. S. Speagle, A. Tramacere, M. Castellano, F. Dubath, A. Galametz, M. Kuemmel, C. Laigle , et al. (148 additional authors not shown)

Abstract: Forthcoming large photometric surveys for cosmology require precise and accurate photometric redshift (photo-z) measurements for the success of their main science objectives. However, to date, no method has been able to produce photo-$z$s at the required accuracy using only the broad-band photometry that those surveys will provide. An assessment of the strengths and weaknesses of current methods i… ▽ More Forthcoming large photometric surveys for cosmology require precise and accurate photometric redshift (photo-z) measurements for the success of their main science objectives. However, to date, no method has been able to produce photo-$z$s at the required accuracy using only the broad-band photometry that those surveys will provide. An assessment of the strengths and weaknesses of current methods is a crucial step in the eventual development of an approach to meet this challenge. We report on the performance of 13 photometric redshift code single value redshift estimates and redshift probability distributions (PDZs) on a common set of data, focusing particularly on the 0.2--2.6 redshift range that the Euclid mission will probe. We design a challenge using emulated Euclid data drawn from three photometric surveys of the COSMOS field. The data are divided into two samples: one calibration sample for which photometry and redshifts are provided to the participants; and the validation sample, containing only the photometry, to ensure a blinded test of the methods. Participants were invited to provide a redshift single value estimate and a PDZ for each source in the validation sample, along with a rejection flag that indicates sources they consider unfit for use in cosmological analyses. The performance of each method is assessed through a set of informative metrics, using cross-matched spectroscopic and highly-accurate photometric redshifts as the ground truth. We show that the rejection criteria set by participants are efficient in removing strong outliers, sources for which the photo-z deviates by more than 0.15(1+z) from the spectroscopic-redshift (spec-z). We also show that, while all methods are able to provide reliable single value estimates, several machine-learning methods do not manage to produce useful PDZs. [abridged] △ Less

Submitted 18 November, 2020; v1 submitted 25 September, 2020; originally announced September 2020.

Comments: Accepted in A&A, 25 pages, 13 figures, 7 tables

Journal ref: A&A 644, A31 (2020)

EinsteinPy: A Community Python Package for General Relativity

Authors: Shreyas Bapat, Ritwik Saha, Bhavya Bhatt, Shilpi Jain, Akshita Jain, Sofía Ortín Vela, Priyanshu Khandelwal, Jyotirmaya Shivottam, Jialin Ma, Gim Seng Ng, Pratyush Kerhalkar, Hrishikesh Sudam Sarode, Rishi Sharma, Manvi Gupta, Divya Gupta, Tushar Tyagi, Tanmay Rustagi, Varun Singh, Saurabh Bansal, Naman Tayal, Abhijeet Manhas, Raphael Reyna, Gaurav Kumar, Govind Dixit, Ratin Kumar , et al. (24 additional authors not shown)

Abstract: 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 increas… ▽ More 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. △ Less

Submitted 22 May, 2020; originally announced May 2020.

Comments: We plan to submit it to A&A. 9 pages

MSC Class: 83-04

A Foreground Model Independent Estimation of Joint Posterior of CMB E mode Polarization over Large Angular Scales

Authors: Ujjal Purkayastha, Vipin Sudevan, Rajib Saha

Abstract: 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. I… ▽ More 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. △ Less

Submitted 18 April, 2020; originally announced April 2020.

Comments: 11 figures and 14 pages

An Accurate Reconstruction of CMB E Mode Signal over Large Angular Scales using Prior Information of CMB Covariance Matrix in ILC Algorithm

Authors: Ujjal Purkayastha, Vipin Sudevan, Rajib Saha

Abstract: 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… ▽ More 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. △ Less

Submitted 30 March, 2020; originally announced March 2020.

Comments: 29 pages, 10 figures

Investigating non-Gaussianity in Cosmic Microwave Background Temperature Maps using Spherical Harmonic Phases

Authors: Sarvesh Kumar Yadav, Rajib Saha

Abstract: In this article, we extend previous studies based on CMB spherical harmonic phases (SHP) to examine the validity of the hypothesis that the temperature field of the CMB is consistent with a Gaussian random field (GRF). The null hypothesis is that the corresponding CMB SHP are independent and identically distributed in terms of a uniform distribution in the interval [0, 2$π$] \citep{1986ApJ...304..… ▽ More In this article, we extend previous studies based on CMB spherical harmonic phases (SHP) to examine the validity of the hypothesis that the temperature field of the CMB is consistent with a Gaussian random field (GRF). The null hypothesis is that the corresponding CMB SHP are independent and identically distributed in terms of a uniform distribution in the interval [0, 2$π$] \citep{1986ApJ...304...15B,2013rossmanith}. We devise a new model-independent method where we use ordered and non-parametric Rao's statistic, based on sample arc-lengths to comprehensively test uniformity and independence of SHP. We performed our analysis on the scales limited by spherical harmonic modes $\le$ 128, to restrict ourselves to signal-dominated regions. To find the non-uniform or dependent sets of SHP, we calculate the statistic for the data and 10000 Monte Carlo simulated uniformly random sets of SHP and use 0.05 and 0.001 $α$ levels to distinguish between statistically significant and highly significant detections. We first establish the performance of our method using simulated Gaussian, non-Gaussian CMB temperature maps, along with observed non-Gaussian 100 and 143 GHz Planck channel maps. We find that our method performs efficiently and accurately in detecting phase correlations generated in all of the non-Gaussian simulations and observed foreground contaminated 100 and 143 GHz Planck channel temperature maps. We apply our method on the 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 SHP corresponding to some of the $m$-modes is non-uniform, some of the $\ell$ mode SHP and neighboring mode pair SHP are correlated in cleaned CMB maps. The detection of non-uniformity or correlation in the SHP indicates the presence of non-Gaussian signals in the foreground minimized CMB maps. △ Less

Submitted 10 January, 2022; v1 submitted 29 January, 2020; originally announced January 2020.

Comments: 42 pages, 23 figures, 6 tables

On the Impact of Random Residual Calibration Error on the Gibbs ILC Estimates over Large Angular Scales

Authors: Vipin Sudevan, Rajib Saha

Abstract: 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 calibrati… ▽ More 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 μK^2$ and $-4.7 μ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. △ Less

Submitted 9 January, 2020; originally announced January 2020.

Comments: 8 pages, 4 figures

Exact Solutions and Constraints on the Dark Energy Model in FRW Universe

Authors: Albin Joseph, Rajib Saha

Abstract: 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 where the scalar field potential is a function of th… ▽ More 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 where the scalar field potential is a function of the 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), we then constrain the model with currents data sets, which includes-Baryon Acoustic Oscillations, Hubble parameter data and Type 1a Supernova data (Pantheon Dataset). 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. △ Less

Submitted 16 September, 2022; v1 submitted 13 December, 2019; originally announced December 2019.

Comments: 21 pages, 7 figures

An Application of Global ILC Algorithm over Large Angular Scales to Estimate CMB Posterior Using Gibbs Sampling

Authors: Vipin Sudevan, Rajib Saha

Abstract: 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 Gi… ▽ More 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. △ Less

Submitted 20 October, 2018; originally announced October 2018.

Comments: 10 pages, 8 figures

A Global ILC Approach in Pixel Space over Large Angular Scales of the Sky using CMB Covariance Matrix

Authors: Vipin Sudevan, Rajib Saha

Abstract: 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 C… ▽ More 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. △ Less

Submitted 28 December, 2017; originally announced December 2017.

Comments: 11 pages, 12 figures

Separating CMB Stokes Q and U polarization signals from Non-Gaussian Emissions

Authors: Ujjal Purkayastha, Rajib Saha

Abstract: 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 st… ▽ More 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. △ Less

Submitted 6 July, 2017; originally announced July 2017.

Comments: 19 pages, 14 figures

An Improved Diffuse Foreground Subtraction by ILC method: CMB Map and Angular Power Spectrum using Planck and WMAP Observations

Authors: Vipin Sudevan, Pavan K. Aluri, Sarvesh Kumar Yadav, Rajib Saha, Tarun Souradeep

Abstract: 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 desi… ▽ More 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. △ Less

Submitted 11 December, 2016; originally announced December 2016.

Comments: 17 pages, 17 figures

Journal ref: The Astrophysical Journal, 842, 62 (2017)

A Perturbative Analysis of Synchrotron Spectral Index Variation over Microwave Sky

Authors: Rajib Saha, Pavan Kumar Aluri

Abstract: 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… ▽ More 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 β_{0s}, over the sky portion being analyzed. Using Monte Carlo simulations we find that, <\hat β_{0s}> = -2.84 \pm 0.01, which matches very closely with position of a peak at β_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 β_d> = 2.00 \pm 0.004, from simulations, matches very well with spatially fixed input thermal dust spectral index β_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. △ Less

Submitted 11 September, 2015; originally announced September 2015.

Comments: 11 pages, 6 figures

Journal ref: The Astrophysical Journal, 829, 113 (2016)

Formation of Cyanoformaldehyde in the interstellar space

Authors: Ankan Das, Liton Majumdar, Sandip K. Chakrabarti, Rajdeep Saha, Sonali Chakrabarti

Abstract: 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 t… ▽ More 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. △ Less

Submitted 17 June, 2013; v1 submitted 7 June, 2013; originally announced June 2013.

Comments: 20 pages, 7 figures

Analysis of WMAP 7-year Temperature Data: Astrophysics of the Galactic Haze

Authors: Davide Pietrobon, Krzysztof M. Gorski, James Bartlett, Anthony J. Banday, Gregory Dobler, Loris P. L. Colombo, Luca Pagano, Graca Rocha, Rajib Saha, Jeffrey B. Jewell, Sergi R. Hildebrandt, Hans Kristian Eriksen, Charles R. Lawrence

Abstract: 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 foregr… ▽ More 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. △ Less

Submitted 30 July, 2012; v1 submitted 25 October, 2011; originally announced October 2011.

Comments: 16 pages, 16 figures, Published on ApJ

Journal ref: ApJ, 755, 69 (2012)

A Foreground Cleaned CMB Map from Non-Gaussianity Measurement

Authors: Rajib Saha

Abstract: 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 o… ▽ More 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. △ Less

Submitted 31 May, 2011; originally announced May 2011.

Comments: 5 pages, 5 figures

Planck early results. VI. The High Frequency Instrument data processing

Authors: Planck HFI Core Team, P. A. R. Ade, N. Aghanim, R. Ansari, M. Arnaud, M. Ashdown, J. Aumont, A. J. Banday, M. Bartelmann, J. G. Bartlett, E. Battaner, K. Benabed, A. Benoît, J. -P. Bernard, M. Bersanelli, J. J. Bock, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger, T. Bradshaw, M. Bucher, J. -F. Cardoso, G. Castex, A. Catalano , et al. (141 additional authors not shown)

Abstract: 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 μK… ▽ More 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 μ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. △ Less

Submitted 23 December, 2011; v1 submitted 11 January, 2011; originally announced January 2011.

Comments: Replaced by the accepted version for publication, as part of a package of papers describing first results of the Planck mission The paper with figures at full resolution and full color tables can also be downloaded from the ESA site http://www.rssd.esa.int/Planck

Report number: Planck2011-1.6

Journal ref: A&A, Volume 536, December 2011

CMB Polarization and Temperature Power Spectra Estimation using Linear Combination of WMAP 5-year Maps

Authors: Pramoda Kumar Samal, Rajib Saha, Jacques Delabrouille, Simon Prunet, Pankaj Jain, Tarun Souradeep

Abstract: 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… ▽ More 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π)$ is 532 $μK^2$, approximately 2.5 times the estimate (213.4 $μK^2$) made by the WMAP team. △ Less

Submitted 20 March, 2009; originally announced March 2009.

Comments: 28 pages 13 figures

Journal ref: Astrophys.J.714:840-851, 2010

Model Independent Foreground Power Spectrum Estimation using WMAP 5-year Data

Authors: Tuhin Ghosh, Rajib Saha, Pankaj Jain, Tarun Souradeep

Abstract: 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 ass… ▽ More 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. △ Less

Submitted 12 January, 2009; originally announced January 2009.

Comments: 12 pages, 4 figures

Journal ref: Phys.Rev.D79:123011,2009

Signals of Statistical Anisotropy in WMAP Foreground-Cleaned Maps

Authors: Pramoda Kumar Samal, Rajib Saha, Pankaj Jain, John P. Ralston

Abstract: 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.… ▽ More 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. △ Less

Submitted 27 November, 2008; v1 submitted 11 November, 2008; originally announced November 2008.

Comments: 32 pages 12 figures

Journal ref: Mon.Not.Roy.Astron.Soc.396:511,2009

On the Low Surface Magnetic Field Structure of Quark Stars

Authors: Nandini Nag, Sutapa Ghosh, Roni Saha, Somenath Chakrabarty

Abstract: 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… ▽ More 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. △ Less

Submitted 27 October, 2008; originally announced October 2008.

Comments: 24 pages REVTEX, one .eps figure, psfig.sty is included

Testing Isotropy of Cosmic Microwave Background Radiation

Authors: Pramoda Kumar Samal, Rajib Saha, Pankaj Jain, John P. Ralston

Abstract: 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 powe… ▽ More 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. △ Less

Submitted 21 August, 2007; originally announced August 2007.

Comments: 40 pages 15 figures

Journal ref: Mon.Not.Roy.Astron.Soc.385:1718,2008

CMB anisotropy power spectrum using linear combinations of WMAP maps

Authors: Rajib Saha, Simon Prunet, Pankaj Jain, Tarun Souradeep

Abstract: 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… ▽ More 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. △ Less

Submitted 25 June, 2007; originally announced June 2007.

Comments: 38 pages, 17 figures

Journal ref: Phys.Rev.D78:023003,2008

CMB power spectrum estimation with non-circular beam and incomplete sky coverage

Authors: Sanjit Mitra, Anand S. Sengupta, Subharthi Ray, Rajib Saha, Tarun Souradeep

Abstract: 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 siz… ▽ More 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. △ Less

Submitted 20 November, 2008; v1 submitted 5 February, 2007; originally announced February 2007.

Comments: 22 pages, 7 figures. Revised presentation to highlight significance of extended results. Matches version accepted to the MNRAS

Report number: IUCAA PREPRINT 06/2007

Journal ref: Mon.Not.Roy.Astron.Soc.394:1419-1439,2009

Non-Circular beam correction to the CMB power spectrum

Authors: Tarun Souradeep, Sanjit Mitra, Anand Sengupta, Subharthi Ray, Rajib Saha

Abstract: 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… ▽ More 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. △ Less

Submitted 23 August, 2006; originally announced August 2006.

Comments: 12 pages, 2 figures; Proceedings of the Fundamental Physics With CMB workshop, UC Irvine, March 23-25, 2006, to be published in New Astronomy Reviews

Report number: IUCAA-41/06

Journal ref: New Astron.Rev.50:1030-1035,2006

Angular power spectrum of CMB anisotropy from WMAP

Authors: Tarun Souradeep, Rajib Saha, Pankaj Jain

Abstract: 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 e… ▽ More 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. △ Less

Submitted 9 August, 2006; originally announced August 2006.

Comments: 12 pages, 4 figures; Proceedings of the Fundamental Physics With CMB workshop, UC Irvine, March 23-25, 2006, to be published in New Astronomy Reviews

Report number: IUCAA-37/06

Journal ref: New Astron.Rev.50:854-860,2006

A re-analysis of the three-year WMAP temperature power spectrum and likelihood

Authors: H. K. Eriksen, Greg Huey, R. Saha, F. K. Hansen, J. Dick, A. J. Banday, K. M. Gorski, P. Jain, J. B. Jewell, L. Knox, D. L. Larson, I. J. O'Dwyer, T. Souradeep, B. D. Wandelt

Abstract: 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… ▽ More 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. △ Less

Submitted 2 October, 2006; v1 submitted 5 June, 2006; originally announced June 2006.

Comments: 14 pages, 7 figures. Accepted for publication in ApJ. Numerical results unchanged, but interpretation sharpened: Likelihood approximation issues at l=13-30 far more important than potential foreground issues at l <= 12. Gibbs products (spectrum and sky samples, and "easy-to-use" likelihood module) available from http://www.astro.uio.no/~hke/ under "Research"

Journal ref: Astrophys.J.656:641-652,2007

Neutron Star Equation of State and the Possibility of Complex Self-Energy in Landau Theory of Fermi Liquid in Presence of Strong Quantizing Magnetic Field

Authors: Soma Mandal, Roni Saha, Sutapa Ghosh, Somenath Chakrabarty

Abstract: Using a relativistic version of Landau theory of Fermi liquid with $σ-ω$ and $ρ$ 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 spe… ▽ More Using a relativistic version of Landau theory of Fermi liquid with $σ-ω$ and $ρ$ 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 $σ$, $ω$ and neutral $ρ$ transfer processes and in the direct interaction diagram with $ρ_\pm$ transfer reactions, the nucleon self-energies become complex in nature. △ Less

Submitted 15 May, 2006; v1 submitted 19 September, 2005; originally announced September 2005.

Comments: REVTEX 14 Pages, no figure

Journal ref: Phys. Rev. C 74 , 015801 (2006)

A blind estimation of the power spectrum of CMB anisotropy from WMAP

Authors: Rajib Saha, Pankaj Jain, Tarun Souradeep

Abstract: 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 WMA… ▽ More 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. △ Less

Submitted 25 June, 2006; v1 submitted 17 August, 2005; originally announced August 2005.

Comments: 4 pages, 3 figures, Matches version accepted to ApJ. Lett. Minor editing; Level of bias at low multipoles mentioned. Suppl. details at http://meghnad.iucaa.ernet.in/~tarun/sup0508383

Report number: IUCAA-34/2005

Journal ref: Astrophys.J.645:L89-L92,2006

Probing Light Pseudoscalars with Light: Propagation, Resonance and Spontaneous Polarization

Authors: Sudeep Das, Pankaj Jain, John P. Ralston, Rajib Saha

Abstract: 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… ▽ More 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. △ Less

Submitted 9 August, 2007; v1 submitted 18 August, 2004; originally announced August 2004.

Comments: 41 pages, 7 figures, revised version

Journal ref: JCAP 0506 (2005) 002