Machine learning lies at the heart of new possibilities for scientific discovery, knowledge generation, and artificial intelligence. Its potential benefits to these fields requires going beyond predictive accuracy and focusing on interpretability. In particular, many scientific problems require interpretations in a domain-specific interpretable feature space (e.g. the frequency domain) whereas attributions to the raw features (e.g. the pixel space) may be unintelligible or even misleading. To address this challenge, we propose TRIM (TRansformation IMportance), a novel approach which attributes importances to features in a transformed space and can be applied post-hoc to a fully trained model. TRIM is motivated by a cosmological parameter estimation problem using deep neural networks (DNNs) on simulated data, but it is generally applicable across domains/models and can be combined with any local interpretation method. In our cosmology example, combining TRIM with contextual decomposition shows promising results for identifying which frequencies a DNN uses, helping cosmologists to understand and validate that the model learns appropriate physical features rather than simulation artifacts.

By combining the photometric, spectroscopic, and astrometric information of the stars in the sightline of SNRs, the distances to and the extinctions of 32 Galactic supernova remnants (SNRs) are investigated. The stellar atmospheric parameters are from the SDSS$-$DR14$/$APOGEE and LAMOST$-$DR5$/$LEGUE spectroscopic surveys. The multi-band photometry, from optical to infrared, are collected from the {\it Gaia}, APASS, Pan--STARRS1, 2MASS, and {\it WISE} surveys. With the calibrated {\it Gaia} distances of individual stars, the distances to 15 of 32 SNRs are well determined from their produced extinction and association with molecular clouds. The upper limits of distance are derived for 3 SNRs. The color excess ratios $E(g_{\rm P1}-\lambda) / E(g_{\rm P1}-r_{\rm P1})$ of 32 SNRs are calculated, and their variation with wavebands is fitted by a simple dust model. The inferred dust grain size distribution bifurcates: while the graphite grains have comparable size to the average ISM dust, the silicate grains are generally larger. Along the way, the average extinction law from optical to near-infrared of the Milky Way is derived from the 1.3 million star sample and found to agree with the CCM89 law with $R_{\rm V}=3.15$.

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

Outflows from starburst galaxies can be driven by thermal pressure, radiation and cosmic rays. We present an analytic phenomenological model that accounts for these contributions simultaneously to investigate their effects on the hydrodynamical properties of outflows. We assess the impact of energy injection, wind opacity, magnetic field strength and the mass of the host galaxy on flow velocity, temperature, density and pressure profiles. For an M82-like wind, a thermally-dominated driving mechanism is found to deliver the fastest and hottest wind. Radiation-driven winds in typical starburst-galaxy configurations are unable to attain the higher flow velocities and temperatures associated with thermal and cosmic ray-driven systems, leading to higher wind densities which would be more susceptible to cooling and fragmentation at lower altitudes. High opacity winds are more sensitive to radiative driving, but terminal flow velocities are still lower than those achieved by other driving mechanisms at realistic opacities. We demonstrate that variations in the outflow magnetic field can influence its coupling with cosmic rays, where stronger fields enable greater streaming but less driving near the base of the flow, instead with cosmic rays redirecting their driving impact to higher altitudes. The gravitational potential is less important in M82-like wind configurations, and substantial variations in the flow profiles only emerge at high altitude in massive haloes. This model offers a more generalised approach to examine the large scale hydrodynamical properties for a wide variety of starburst galaxies.

Statistical inference of analytically non-tractable posteriors is a difficult problem because of marginalization of correlated variables and stochastic methods such as MCMC and VI are commonly used. We argue that stochastic KL divergence minimization used by MCMC and VI is noisy, and we propose instead EL_2O, expectation optimization of L_2 distance squared between the approximate log posterior q and the un-normalized log posterior of p. When sampling from q the solutions agree with stochastic KL divergence minimization based VI in the large sample limit, however EL_2O method is free of sampling noise, has better optimization properties, and requires only as many sample evaluations as the number of parameters we are optimizing if q covers p. As a consequence, increasing the expressivity of q improves both the quality of results and the convergence rate, allowing EL_2O to approach exact inference. Use of automatic differentiation methods enables us to develop Hessian, gradient and gradient free versions of the method, which can determine M(M+2)/2+1, M+1 and 1 parameter(s) of q with a single sample, respectively. EL_2O provides a reliable estimate of the quality of the approximating posterior, and converges rapidly on full rank gaussian approximation for q and extensions beyond it, such as nonlinear transformations and gaussian mixtures. These can handle general posteriors, while still allowing fast analytic marginalizations. We test it on several examples, including a realistic 13 dimensional galaxy clustering analysis, showing that it is several orders of magnitude faster than MCMC, while giving smooth and accurate non-gaussian posteriors, often requiring a few to a few dozen of iterations only.

The Simons Observatory (SO) is a ground-based cosmic microwave background (CMB) experiment sited on Cerro Toco in the Atacama Desert in Chile that promises to provide breakthrough discoveries in fundamental physics, cosmology, and astrophysics. Supported by the Simons Foundation, the Heising-Simons Foundation, and with contributions from collaborating institutions, SO will see first light in 2021 and start a five year survey in 2022. SO has 287 collaborators from 12 countries and 53 institutions, including 85 students and 90 postdocs. The SO experiment in its currently funded form ('SO-Nominal') consists of three 0.4 m Small Aperture Telescopes (SATs) and one 6 m Large Aperture Telescope (LAT). Optimized for minimizing systematic errors in polarization measurements at large angular scales, the SATs will perform a deep, degree-scale survey of 10% of the sky to search for the signature of primordial gravitational waves. The LAT will survey 40% of the sky with arc-minute resolution. These observations will measure (or limit) the sum of neutrino masses, search for light relics, measure the early behavior of Dark Energy, and refine our understanding of the intergalactic medium, clusters and the role of feedback in galaxy formation. With up to ten times the sensitivity and five times the angular resolution of the Planck satellite, and roughly an order of magnitude increase in mapping speed over currently operating ("Stage 3") experiments, SO will measure the CMB temperature and polarization fluctuations to exquisite precision in six frequency bands from 27 to 280 GHz. SO will rapidly advance CMB science while informing the design of future observatories such as CMB-S4.

We present three dimensional (3D) dust mapping of 12 supernova remnants (SNRs) in the Galactic anti-center (Galactic longitude $l$ between 150\degr\ and 210\degr) based on a recent 3D interstellar extinction map. The dust distribution of the regions which cover the full extents in the radio continuum for the individual SNRs are discussed. Four SNRs show significant spatial coincidences between molecular clouds (MCs) revealed from the 3D extinction mapping and the corresponding radio features. The results confirm the interactions between these SNRs and their surrounding MCs. Based on these correlations, we provide new distance estimates of the four SNRs, G189.1+3.0 (IC443, $d=1729^{+116}_{-94} \rm \,pc$), G190.9-2.2 ($d=1036^{+17}_{-81} \rm \,pc$), G205.5+0.5 ($d=941^{+96}_{-94}$ or $1257^{+92}_{-101} \rm \,pc$) and G213.0-0.6 ($d=1146^{+79}_{-80} \rm \,pc$). In addition, we find indirect evidences of potential interactions between SNRs and MCs for three other SNRs. New distance constraints are also given for these three SNRs.

The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at $2 < z < 5$, together with the ability to efficiently target high-$z$ galaxies with known techniques, enables large gains in the study of Inflation and Dark Energy. A future spectroscopic survey can test the Gaussianity of the initial conditions up to a factor of ~50 better than our current bounds, crossing the crucial theoretical threshold of $\sigma(f_{NL}^{\rm local})$ of order unity that separates single field and multi-field models. Simultaneously, it can measure the fraction of Dark Energy at the percent level up to $z = 5$, thus serving as an unprecedented test of the standard model and opening up a tremendous discovery space.

The Simons Observatory (SO) is a new cosmic microwave background experiment being built on Cerro Toco in Chile, due to begin observations in the early 2020s. We describe the scientific goals of the experiment, motivate the design, and forecast its performance. SO will measure the temperature and polarization anisotropy of the cosmic microwave background in six frequency bands: 27, 39, 93, 145, 225 and 280 GHz. The initial configuration of SO will have three small-aperture 0.5-m telescopes (SATs) and one large-aperture 6-m telescope (LAT), with a total of 60,000 cryogenic bolometers. Our key science goals are to characterize the primordial perturbations, measure the number of relativistic species and the mass of neutrinos, test for deviations from a cosmological constant, improve our understanding of galaxy evolution, and constrain the duration of reionization. The SATs will target the largest angular scales observable from Chile, mapping ~10% of the sky to a white noise level of 2 $\mu$K-arcmin in combined 93 and 145 GHz bands, to measure the primordial tensor-to-scalar ratio, $r$, at a target level of $\sigma(r)=0.003$. The LAT will map ~40% of the sky at arcminute angular resolution to an expected white noise level of 6 $\mu$K-arcmin in combined 93 and 145 GHz bands, overlapping with the majority of the LSST sky region and partially with DESI. With up to an order of magnitude lower polarization noise than maps from the Planck satellite, the high-resolution sky maps will constrain cosmological parameters derived from the damping tail, gravitational lensing of the microwave background, the primordial bispectrum, and the thermal and kinematic Sunyaev-Zel'dovich effects, and will aid in delensing the large-angle polarization signal to measure the tensor-to-scalar ratio. The survey will also provide a legacy catalog of 16,000 galaxy clusters and more than 20,000 extragalactic sources.

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

With low redshift probes reaching unprecedented precision, uncertainty of the CMB optical depth is expected to be the limiting factor for future cosmological neutrino mass constraints. In this paper, we discuss to what extent combinations of CMB lensing and galaxy surveys measurements at low redshifts $z\sim 0.5-5$ will be able to make competitive neutrino mass measurements without relying on any optical depth constraints. We find that the combination of LSST galaxies and CMB-S4 lensing should be able to achieve constraints on the neutrino mass sum of 25meV without optical depth information, an independent measurement that is competitive with or slightly better than the constraint of 30meV possible with CMB-S4 and present-day optical depth measurements. These constraints originate both in structure growth probed by cross-correlation tomography over a wide redshift range as well as, most importantly, the shape of the galaxy power spectrum measured over a large volume. We caution that possible complications such as higher-order biasing and systematic errors in the analysis of high redshift galaxy clustering are only briefly discussed and may be non-negligible. Nevertheless, our results show that new kinds of high-precision neutrino mass measurements at and beyond the present-day optical depth limit may be possible.

Double neutron star (DNS) merger events are promosing candidates of short Gamma-ray Burst (sGRB) progenitors as well as high-frequecy gravitational wave (GW) emitters. On August 17, 2017, such a coinciding event was detected by both the LIGO-Virgo gravitational wave detector network as GW170817 and Gamma-Ray Monitor on board NASA's {\it Fermi} Space Telescope as GRB 170817A. Here we show that the fluence and spectral peak energy of this sGRB fall into the lower portion of the distributions of known sGRBs. Its peak isotropic luminosity is abnormally low. The estimated event rate density above this luminosity is at least $190^{+440}_{-160} {\rm Gpc^{-3} \ yr^{-1}}$, which is close to but still below the DNS merger event rate density. This event likely originates from a structured jet viewed from a large viewing angle. There are similar faint soft GRBs in the {\it Fermi} archival data, a small fraction of which might belong to this new population of nearby, low-luminosity sGRBs.

Delensing, the removal of the limiting lensing B-mode background, is crucial for the success of future cosmic microwave background (CMB) surveys in constraining inflationary gravitational waves (IGWs). In recent work, delensing with large-scale structure tracers has emerged as a promising method both for improving constraints on IGWs and for testing delensing methods for future use. However, the delensing fractions (i.e., the fraction of the lensing-B mode power removed) achieved by recent efforts have been only $20-30\%$. In this work, we provide a detailed characterization of a full-sky, dust-cleaned cosmic infrared background (CIB) map for delensing and construct a further-improved delensing template by adding additional tracers to increase delensing performance. In particular, we build a multitracer delensing template by combining the dust-cleaned Planck CIB map with a reconstructed CMB lensing map from Planck and a galaxy number density map from the Wide-field Infrared Survey Explorer (WISE) satellite. For this combination, we calculate the relevant weightings by fitting smooth templates to measurements of all the cross- and auto-spectra of these maps. On a large fraction of the sky ($f_\mathrm{sky}=0.43$), we demonstrate that our maps are capable of providing a delensing factor of $43 \pm 1\%$; using a more restrictive mask ($f_\mathrm{sky}=0.11$), the delensing factor reaches $48 \pm 1\%$. For low-noise surveys, our delensing maps, which cover much of the sky, can thus improve constraints on the tensor-to-scalar ratio ($r$) by nearly a factor of 2. The delensing tracer maps are made publicly available, and we encourage their use in ongoing and upcoming B-mode surveys.

We investigate the impact of a variety of analysis assumptions that influence cluster identification and location on the kSZ pairwise momentum signal and covariance estimation. Photometric and spectroscopic galaxy tracers from SDSS, WISE, and DECaLs, spanning redshifts $0.05<z<0.7$, are considered in combination with CMB data from Planck and WMAP. With two complementary techniques, analytic offset modeling and direct comparisons of redMaPPer brightest and central catalog samples, we find that miscentering uncertainties average to $0.4-0.7\sigma$ for the Planck kSZ statistical error budget obtained with a jackknife estimator. We also find that jackknife covariance estimates are significantly more conservative than those obtained by CMB rotation methods. Using redMaPPer data, we concurrently compare the impact of photometric redshift errors and miscentering. At separations $<\sim 50$ Mpc, where the kSZ signal is largest, miscentering uncertainties can be comparable to JK errors, while photometric redshifts are lower but still significant. For the next generation of CMB and LSS surveys the statistical and photometric errors will shrink markedly. Our results demonstrate that uncertainties introduced through using galaxy proxies for cluster locations will need to be fully incorporated, and actively mitigated, for the kSZ to reach its full potential as a cosmological constraining tool for dark energy and neutrino physics.

Weakly interacting massive particles are a widely well-probed dark matter candidate by the dark matter direct detection experiments. Theoretically, there are a large number of ultraviolet completed models that consist of a weakly interacting massive particle dark matter. The variety of models makes the comparison with the direct detection data complicated and often non-trivial. To overcome this, in the non-relativistic limit, the effective theory was developed in the literature which works very well to significantly reduce the complexity of dark matter-nucleon interactions and to better study the nuclear response functions. In the effective theory framework for a spin-1/2 dark matter, we combine three independent likelihood functions from the latest PandaX, LUX, and XENON1T data, and give a joint limit on each effective coupling. The astrophysical uncertainties of the dark matter distribution are also included in the likelihood. We further discuss the isospin violating cases of the interactions. Finally, for both dimension-five and dimension-six effective theories above the electroweak scale, we give updated limits of the new physics mass scales.

We present a three dimensional (3D) extinction analysis in the region toward the supernova remnant (SNR) S147 (G180.0-1.7) using multi-band photometric data from the Xuyi Schmidt Telescope Photometric Survey of the Galactic Anticentre (XSTPS-GAC), 2MASS and WISE. We isolate a previously unrecognised dust structure likely to be associated with SNR S147. The structure, which we term as "S147 dust cloud", is estimated to have a distance $d$ = 1.22 $\pm$ 0.21 kpc, consistent with the conjecture that S147 is associated with pulsar PSR J0538 + 2817. The cloud includes several dense clumps of relatively high extinction that locate on the radio shell of S147 and coincide spatially with the CO and gamma-ray emission features. We conclude that the usage of CO measurements to trace the SNR associated MCs is unavoidably limited by the detection threshold, dust depletion, and the difficulty of distance estimates in the outer Galaxy. 3D dust extinction mapping may provide a better way to identify and study SNR-MC interactions.

We carry out a systematical study of the spectral lag properties of 50 single-pulsed Gamma-Ray Bursts (GRBs) detected by Fermi/GBM. By dividing the light curves into multiple consecutive energy channels we provide a new measurement of the spectral lag which is independent on energy channel selections. We perform a detailed statistical study of our new measurements. We find two similar power-law energy dependencies of both the pulse arrival time and pulse width. Our new results on the power-law indices would favor the relativistic geometric effects for the origin of spectral lag. However, a complete theoretical framework that can fully account for the diverse energy dependencies of both arrival time and pulse width revealed in this work is still missing. We also study the spectral evolution behaviors of the GRB pulses. We find that the GRB pulse with negligible spectral lag would usually have a shorter pulse duration and would appear to have a "hardness-intensity tracking" (HIT) behavior and the GRB pulse with a significant spectral lag would usually have a longer pulse duration and would appear to have a "hard-to-soft" (HTS) behavior.

The coupling of a scalar field to electromagnetic field via the Chern-Simons term will rotate the polarization directions of the cosmic microwave background radiation. The rotation angle which relies on the distribution of the scalar field on the CMB sky is direction dependent. Such anisotropies will give rise to new distortions to the power spectra of CMB polarization and it can be used to probe the detailed physics of the scalar field. In this paper we use the updated observational data to constrain the anisotropic rotation angle in a model independent way. We find that the dominant effect of the anisotropic rotation on CMB comes from its variance and it is constrained tightly by the current data.

The luminescent properties of CsI(Na) crystals are studied in this report. By using a TDS3054C oscilloscope with a sampling frequency of 5 GS/s, we find out that nuclear recoil signals are dominated by very fast light pulse with a decay time of ~20 ns, while {\gamma}-ray signals have a decay time of ~600 ns. The wavelength of nuclear recoil and {\gamma}-ray signals are also different. The study of n/{\gamma} separation shows that the rejection factor can reach an order of 10-7 with signal efficiency more than 80% at an equivalent electron recoil energy of 20 keV or more. Such a property makes CsI(Na) an ideal candidate for dark matter searches.

In this paper, we study the motion of photons around a Kehagias-Sfetsos (KS) black hole and obtain constraints on IR modified Ho$\check{r}$ava gravity without cosmological constant ($\sim \Lambda_{W}$). An analytic formula for the light deflection angle is obtained. For a propagating photon, the deflection angle $\delta \phi$ increases with large values of the Ho$\check{r}$ava gravity parameter $\omega$. Under the UV limit $\omega \longrightarrow \infty$, deflection angle reduces to the result of usual Schwarzschild case, $4GM/R$. It is also found that with increasing scale of astronomical observation system the Ho$\check{r}$ava-Lifshitz gravity should satisfy $|\omega M^2|>1.1725 \times10^{-16}$ with 12% precision for Earth system, $|\omega M^2| > 8.27649 \times 10^{-17}$ with 17% precision for Jupiter system and $|\omega M^2| > 8.27650\times 10^{-15}$ with 0.17% precision for solar system.

The 28 October 2003 flare gave us the unique opportunity to compare the acceleration time of high-energy protons with the escaping time of those particles which have been measured onboard spacecraft and by neutron monitors network as GLE event. High-energy emission time scale and shock wave height and velocity time dependencies were also studied.

In this paper we generalize the flux ratio method Bailey et al. (2009) to the case of two luminosity indicators and search the optimal luminosity-flux ratio relations on a set of spectra whose phases are around not only the date of bright light but also other time. With these relations, a new method is proposed to constrain the host galaxy extinction of SN Ia and its distance. It is first applied to the low redshift supernovas and then to the high redshift ones. The results of the low redshift supernovas indicate that the flux ratio method can indeed give well constraint on the host galaxy extinction parameter E(B-V), but weaker constraints on R_{V}. The high redshift supernova spectra are processed by the same method as the low redshift ones besides some differences due to their high redshift. Among 16 high redshift supernovas, 15 are fitted very well except 03D1gt. Based on these distances, Hubble diagram is drew and the contents of the Universe are analyzed. It supports an acceleration behavior in the late Universe. Therefore, the flux ratio method can give constraints on the host galaxy extinction and supernova distance independently. We believe, through further studies, it may provide a precise tool to probe the acceleration of the Universe than before.

The large scatters of luminosity relations of gamma-ray bursts (GRBs) have been one of the most important reasons that prevent the extensive applications of GRBs in cosmology. In this paper, we extend the two-dimensional (2D) luminosity relations with $\tau_{\mathrm{lag}}$, $V$, $E_{\mathrm{peak}}$, and $\tau_{\mathrm{RT}}$ as the luminosity indicators to three dimensions (3D) using the same set of luminosity indicators to explore the possibility of decreasing the intrinsic scatters. We find that, for the 3D luminosity relations between the luminosity and an energy scale ($E_{\mathrm{peak}}$) and a time scale ($\tau_{\mathrm{lag}}$ or $\tau_{\mathrm{RT}}$), their intrinsic scatters are considerably smaller than those of corresponding 2D luminosity relations. Enlightened by the result and the definition of the luminosity (energy released in units of time), we discussed possible reasons behind, which may give us helpful suggestions on seeking more precise luminosity relations for GRBs in the future.

In this paper, we revisit the weak gravitational lensing of CMB. The key point is that the deflection angle of CMB photon varies with the redshift. From this we obtain a set of modified power spectra. Compared with the ones calculated before, the temperature, the polarization E mode and the cross power spectra correction are of the order $0.01 %$, while the polarization B mode correction is about $0.4 %$ (to $l=2000$).

The usual inflationary scenario predicts a Gaussian random primordial density fluctuation, different Fourier modes of which do not correlate with each other. In this paper we propose a correlation between these different modes. A simple case is that these different Fourier modes correlate with each other following a Gaussian function. For such a primordial density fluctuation we calculate the CMB angular power spectrum and find that its amplitude decreases but the decrease is different for different $l$. This feature can be used to constrain the the correlation strength from the real data.

Since the inhomogeneous instrument noise can produce extra non-Gaussianity in the CMB anisotropy, its effect should be carefully subtracted in the primordial non-Gaussianity estimation. We calculate the probability distribution function of the CMB anisotropy for local type of non-Gaussianity, from which the optimal estimator in the general case (inhomogeneous noise and cut sky) is obtained. The new estimator obtained here is different from the popular one, since the inhomogeneous noise and cut sky effects are completely accounted. The CMB anisotropy in the new estimator is noise weighted. The noise weight is different from that used by WMAP Group in their 5-year data analysis. Although it is still difficult to calculate the new estimator rigorously, for the case of the slightly inhomogeneous noise, there exists a series expansion method to compute the new estimator. Each order in the series is suppressed by two factors, $(\frac{\sigma^2}{\sigma^2_{i}}-1)$ and $\frac{C_l}{C^{tot}_{l}}$, which make the method feasible. Through the Edgeworth expansion we can generalize our discussion to other types of non-Gaussianity.