We explore the cosmological multi-tracer synergies between a DESI-like Emission Line Galaxy (ELG) distribution and a Tianlai-like 21cm intensity map. We use simulated maps generated from a particle simulation in the light-cone volume (Horizon Run 4), sky-trimmed and including the effects of foreground contamination, its removal, and instrument noise. We first validate how the foreground residual affects the recovered 21cm signal by putting different levels of foreground contamination into the 21cm maps. We find that the contamination cannot be ignored in the angular auto-power spectra of HI even when it is small, but has no influence on the accuracy of the angular cross-power spectra between HI and galaxies. In the foreground-cleaned map case, as information is lost in the cleaning procedure, there is also a bias in the cross-power spectrum. However, we found that the bias from the cross-power spectrum is scale-independent, which is easily parameterised as part of the model, while the offset in the HI auto-power spectrum is non-linear. In particular, we tested that the cross power also benefits from the cancellation of the bias in the power spectrum measurement that is induced by the instrument noise, which changes the shape of the auto-power spectra but leaves the cross-power unaffected. We then modelled the angular cross-power spectra to fit the BAO feature in broadband shape of the angular cross power spectrum, including contamination from the residual foreground and the effect of instrument noise. We forecast a constraint on the angular diameter distance $D_\mathrm{A}$ for the Tianlai Pathfinder redshift $0.775<z<1.03$, giving a distance measurement with a precision of 3.3% at that redshift.

To reach the full potential for the next generation of weak lensing surveys, it is necessary to mitigate the contamination of intrinsic alignments (IA) of galaxies in the observed cosmic shear signal. The self calibration (SC) of intrinsic alignments provides an independent method to measure the IA signal from the survey data and the photometric redshift information. It operates differently from the marginalization method based on IA modeling. In this work, we present the first application of SC to the KiDS450 data and the KV450 data, to split directly the intrinsic shape - galaxy density (Ig) correlation signal and the gravitational shear - galaxy density (Gg) correlation signal, using the information from photometric redshift (photo-z). We achieved a clear separation of the two signals and performed several validation tests. Our measured signals are found to be in general agreement with the KiDS450 cosmic shear best-fit cosmology, for both lensing and IA measurements. For KV450, we use partial (high-z) data, and our lensing measurements are also in good agreement with KV450 cosmic shear best-fit, however, our IA signal suggests a larger IA amplitude. We discussed the impact of photo-z quality on IA detection and several other potential systematic biases. Finally, we discuss the potential application of the information extracted for both the lensing signal and the IA signal in future surveys.

Galaxy intrinsic alignment (IA) is both a source of systematic contamination to cosmic shear measurement and its cosmological applications, and a source of valuable information on the large scale structure of the universe and galaxy formation. The self-calibration (SC) method \citep{SC2008} was designed to separate IA from cosmic shear, free of IA modeling. It was first successfully applied to the KiDS450 and KV450 data \citep{Yao2019}. We apply the SC method to the DECaLS DR3 shear + photo-z catalog and significantly improve the IA detection to $\sim 20\sigma$. We find a strong dependence of IA on galaxy color, with strong IA signal for red galaxies, and weak/vanishing IA signal for blue galaxies. The detected IA for red galaxies is in reasonable agreement with the non-linear tidal alignment model and the inferred IA amplitude increases with redshift. We discuss various caveats and possible improvements in the measurement, theory and parameter fitting that will be addressed in future works.

The distribution of cosmological neutral hydrogen will provide a new window into the large-scale structure of the Universe with the next generation of radio telescopes and surveys. The observation of this material, through 21cm line emission, will be confused by foreground emission in the same frequencies. Even after these foregrounds are removed, the reconstructed map may not exactly match the original cosmological signal, which will introduce systematic errors and offset into the measured correlations. In this paper, we simulate future surveys of neutral hydrogen using the Horizon Run 4 (HR4) cosmological N-body simulation. We generate HI intensity maps from the HR4 halo catalogue, and combine with foreground radio emission maps from the Global Sky Model, to create accurate simulations over the entire sky. We simulate the HI sky for the frequency range 700-800 MHz, matching the sensitivity of the Tianlai pathfinder. We test the accuracy of the fastICA, PCA and log-polynomial fitting foreground removal methods to recover the input cosmological angular power spectrum and measure the parameters. We show the effect of survey noise levels and beam sizes on the recovered the cosmological constraints. We find that while the reconstruction removes power from the cosmological 21cm distribution on large-scales, we can correct for this and recover the input parameters in the noise-free case. However, the effect of noise and beam size of the Tianlai pathfinder prevents accurate recovery of the cosmological parameters when using only intensity mapping information.

Weak gravitational lensing is one of the most promising cosmological probes to constrain dark matter, dark energy and the nature of gravity at cosmic scales. Intrinsic alignments (IA) of galaxies have been recognized as one of the most serious systematic effects facing gravitational lensing. Such alignments must be isolated and removed to obtain a pure lensing signal. Furthermore, the alignments are related to the processes of galaxy formation, so their extracted signal can help in understanding such formation processes and improving their theoretical modeling. We report in this letter the first detection of the gravitational shear--intrinsic shape (GI) correlation and the intrinsic shape--galaxy density (Ig) correlation in a photometric redshift survey using the self-calibration method. These direct measurements are made from the KiDS-450 photometric galaxy survey with a significance of 2.74$\sigma$ in the third bin for the Ig correlation, and 2.73$\sigma$ for the GI cross-correlation between the third and fourth bins. The self-calibration method uses the information available from photometric surveys without needing to specify an IA model and will play an important role in validating IA models and IA mitigation in future surveys such as LSST, Euclid and WFIRST.

The Analytical Blind Separation (ABS) method is dedicated to extract the CMB signal from its contaminants without any prior information about foregrounds. In this study, we apply the ABS method to simulated polarization sky maps in order to test its ability in the recovery of the CMB $E$- and $B$-mode power spectra. We simulate multi-frequency $Q$ and $U$ microwave sky maps in the presence of two polarized foreground components (synchrotron and thermal dust), together with an uncorrelated Gaussian distributed instrumental noise, measured with a hypothetical PICO-like experiment operating in 10-frequency bands in the range of 30 to 321 GHz. Both full and partial sky surveys are taken into account in our analysis. For the full sky case with the tensor-to-scalar ratios of $r=0$ and $r=0.05$, the ABS can accurately recover the $E$-mode and $B$-mode power spectra with a relative error below $20\%$ with respect to the input CMB power spectra for the full multipole range ($2 \leq \ell \leq 1050$), close to theoretical expectations. For the partial sky analysis with a E/B separation proposed by \cite{ebmixture10} to avoid "$E$-to-$B$" leakage, the results are comparable to the full sky case. We find that the ABS is able to estimate both CMB $E$-mode and $B$-mode power spectra within 1-$\sigma$ at all scales for full sky, and for most of the scales ($50 \leq \ell \leq 1050$) for partial sky. For low-$\ell$ multipoles ($\ell \leq 50$), a noticeable divergence occurs in the partial sky analysis, which might be due to uncertainties originated from both the ABS and the E/B separation. A more detailed analysis for the partial sky case is thus needed for estimating primordial $B$-modes unbiasedly and accurately.

In this study, we apply the Analytical method of Blind Separation (ABS) of the cosmic microwave background (CMB) from foregrounds to estimate the CMB temperature power spectrum from multi-frequency microwave maps. We test the robustness of the ABS estimator and assess the accuracy of the power spectrum recovery by using realistic simulations based on the seven-frequency Planck data, including various frequency-dependent and spatially-varying foreground components (synchrotron, free-free, thermal dust and anomalous microwave emission), as well as an uncorrelated Gaussian-distributed instrumental noise. Considering no prior information about the foregrounds, the ABS estimator can analytically recover the CMB power spectrum over almost all scales with less than $0.5\%$ error for maps where the Galactic plane region ($|b|<10^{\circ}$) is masked out. To further test the flexibility and effectiveness of the ABS approach in a variety of situations, we apply the ABS to the simulated Planck maps in three cases: (1) without any mask, (2) imposing a two-times-stronger synchrotron emission and (3) including only the Galactic plane region ($|b|<10^{\circ}$) in the analysis. In such extreme cases, the ABS approach can still provide an unbiased estimate of band powers at the level of 1 $\mu\rm{K}^2$ on average over all $\ell$ range, and the recovered powers are consistent with the input values within 1-$\sigma$ for most $\ell$ bins.

We introduce a self-calibration method that can be applied to the intrinsic ellipticity--intrinsic ellipticity (II) and gravitational shear -- intrinsic ellipticity (GI) types of intrinsic alignment of galaxies. The method combines previous self-calibration techniques with modifications to one of them in order to use auto-spectra in addition to cross-spectra between redshift bins. This allows one to use the self-calibration while preserving all the constraining power of cosmic shear from surveys. We show that the new method provides more flexibility in using various redshift bin widths. We perform cosmological parameter constraint forecast when this method is applied to the Large Synoptic Survey Telescope (LSST). Compared to the original self-calibration, we find that the new method provides further significant reduction in any residual shift in the cosmological parameters (e.g. factors of $2-4$ for the dark energy equation of state) which is promising for accurate cosmology.

The 9.7$\mu$m interstellar spectral feature, arising from the Si--O stretch of amorphous silicate dust, is the strongest extinction feature in the infrared (IR). In principle, the spectral profile of this feature could allow one to diagnose the mineralogical composition of interstellar silicate material. However, observationally, the 9.7$\mu$m interstellar silicate extinction profile is not well determined. Here we utilize the Spitzer/IRS spectra of five early-type (one O- and four B-type) stars and compare them with that of unreddened stars of the same spectral type to probe the interstellar extinction of silicate dust around 9.7$\mu$m. We find that, while the silicate extinction profiles all peak at ~9.7$\mu$m, two stars exhibit a narrow feature of FWHM ~2.0$\mu$m and three stars display a broad feature of FWHM ~3.0$\mu$m. We also find that the width of the 9.7$\mu$m extinction feature does not show any environmental dependence. With a FWHM of ~2.2$\mu$m, the mean 9.7\mu m extinction profile, obtained by averaging over our five stars, closely resembles that of the prototypical diffuse interstellar medium along the lines of sight toward Cyg OB2 No.12 and WR 98a. Finally, an analytical formula is presented to parameterize the interstellar extinction in the IR at $0.9\mu {\rm m} < \lambda < 15\mu {\rm m}$.

Intrinsic alignments (IA) of galaxies have been recognized as one of the most serious contaminants to weak lensing. These systematics need to be isolated and mitigated in order for ongoing and future lensing surveys to reach their full potential. The IA self-calibration (SC) method was shown in previous studies to be able to reduce the GI contamination by up to a factor of 10 for the 2-point and 3-point correlations. The SC method does not require to assume an IA model in its working and can extract the GI signal from the same photo-z survey offering the possibility to test and understand structure formation scenarios and their relationship to IA models. In this paper, we study the effects of the IA SC mitigation method on the precision and accuracy of cosmological parameter constraints from future cosmic shear surveys LSST, WFIRST and Euclid. We perform analytical and numerical calculations to estimate the loss of precision and the residual bias in the best fit cosmological parameters after the self-calibration is performed. We take into account uncertainties from photometric redshifts and the galaxy bias. We find that the confidence contours are slightly inflated from applying the SC method itself while a significant increase is due to the inclusion of the photo-z uncertainties. The bias of cosmological parameters is reduced from several-$\sigma$, when IA is not corrected for, to below 1-$\sigma$ after SC is applied. These numbers are comparable to those resulting from applying the method of marginalizing over IA model parameters despite the fact that the two methods operate very differently. We conclude that implementing the SC for these future cosmic-shear surveys will not only allow one to efficiently mitigate the GI contaminant but also help to understand their modeling and link to structure formation.

We derive the Sun's offset from the local mean Galactic plane($z_\odot$) using the observed $z$ distribution of young pulsars. Pulsar distances are obtained from measurements of annual parallax, HI absorption spectra or associations where available and otherwise from the observed pulsar dispersion and a model for the distribution of free electrons in the Galaxy. We fit the cumulative distribution function for a ${\rm sech}^2(z)$ distribution function, representing an isothermal self-gravitating disk, with uncertainties being estimated using the bootstrap method. We take pulsars having characteristic age $\tau_c<10^{6.5}$~yr and located within 4.5~kpc of the Sun, omitting those within the local spiral arm and those significantly affected by the Galactic warp, and solve for $z_\odot$ and the scale height, $H$, for different cutoffs in $\tau_c$. We compute these quantities using just the independently determined distances, and these together with DM-based distances separately using the YMW16 and NE2001 Galactic electron density models. We find that an age cutoff at $10^{5.75}$~yr with YMW16 DM-distances gives the best results with a minimum uncertainty in $z_\odot$ and an asymptotically stable value for $H$ showing that, at this age and below, the observed pulsar $z$-distribution is dominated by the dispersion in their birth locations. From this sample of 115 pulsars, we obtain $z_\odot=13.4\pm$4.4~pc and $H=56.9\pm$6.5~pc, similar to estimated scale heights for OB stars and open clusters. Consistent results are obtained using the independent-only distances and using the NE2001 model for the DM-based distances.

We present a new model for the distribution of free electrons in the Galaxy, the Magellanic Clouds and the intergalactic medium (IGM) that can be used to estimate distances to real or simulated pulsars and fast radio bursts (FRBs) based on their dispersion measure (DM). The Galactic model has an extended thick disk representing the so-called warm interstellar medium, a thin disk representing the Galactic molecular ring, spiral arms based on a recent fit to Galactic HII regions, a Galactic Center disk and seven local features including the Gum Nebula, Galactic Loop I and the Local Bubble. An offset of the Sun from the Galactic plane and a warp of the outer Galactic disk are included in the model. Parameters of the Galactic model are determined by fitting to 189 pulsars with independently determined distances and DMs. Simple models are used for the Magellanic Clouds and the IGM. Galactic model distances are within the uncertainty range for 86 of the 189 independently determined distances and within 20\% of the nearest limit for a further 38 pulsars. We estimate that 95\% of predicted Galactic pulsar distances will have a relative error of less than a factor of 0.9. The predictions of YMW16 are compared to those of the TC93 and NE2001 models showing that YMW16 performs significantly better on all measures. Timescales for pulse broadening due to interstellar scattering are estimated for (real or simulated) Galactic and Magellanic Cloud pulsars and FRBs.

We present a three dimensional extinction map in $r$ band. The map has a spatial angular resolution, depending on latitude, between 3 -- 9\,arcmin and covers the entire XSTPS-GAC survey area of over 6,000\,$\rm deg^2$ for Galactic longitude $\rm 140 \leq$ $l$ $ \leq 220\deg$ and latitude $\rm -40\leq$ $b$ $ \leq 40\deg$. By cross-matching the photometric catalog of the Xuyi Schmidt Telescope Photometric Survey of the Galactic Anticentre (XSTPS-GAC) with those of 2MASS and WISE, we have built a multi-band photometric stellar sample of about 30 million stars and applied spectral energy distribution (SED) fitting to the sample. By combining photometric data from the optical to the near-infrared, we are able to break the degeneracy between the intrinsic stellar colours and the amounts of extinction by dust grains for stars with high photometric accuracy, and trace the extinction as a function of distance for low Galactic latitude and thus highly extincted regions. This has allowed us to derive the best-fit extinction and distance information of more than 13 million stars, which are used to construct the three dimensional extinction map. We have also applied a Rayleigh-Jeans colour excess (RJCE) method to the data using the 2MASS and WISE colour $(H-W2)$. The resulting RJCE extinction map is consistent with the integrated two dimensional map deduced using the best-fit SED algorithm. However for individual stars, the amounts of extinction yielded by the RJCE method suffer from larger errors than those given by the best-fit SED algorithm.

We present a catalog of RR Lyrae stars (RRLs) observed by the Xuyi Schmidt Telescope Photometric Survey (XDSS). The area we consider is located in the North Galactic Cap, covering 376.75 sq deg at RA $\approx$ 150 deg and Dec $\approx$ 27 deg down to a magnitude limit of i $\approx$ 19. Using the variability information afforded by the multi-epoch nature of our XDSS data, combined with colors from the Sloan Digital Sky Survey, we are able to identify candidate RRLs. We find 318 candidates, derive distances to them and estimate the detection efficiency. The majority of our candidates have more than 12 observations and for these we are able to calculate periods. These also allows us to estimate our contamination level, which we predict is between 30% to 40%. Finally we use the sample to probe the halo density profile in the 9-49 kpc range and find that it can be well fitted by a double power law. We find good agreement between this model and the models derived for the South Galactic Cap using the Watkins et al. (2009) and Sesar et al. (2010) RRL data-sets, after accounting for possible contamination in our data-set from Sagittarius stream members. We consider non-spherical double power law models of the halo density profile and again find agreement with literature data-sets, although we have limited power to constrain the flattening due to our small survey area. Much tighter constraints will be placed by current and future wide-area surveys, most notably ESA's astrometric Gaia mission. Our analysis demonstrates that surveys with a limited number of epochs can effectively be mined for RRLs. Our complete sample is provided as accompanying online material.

We present measurements of the optical broadband atmospheric extinction coefficients and the night sky brightness at the Xuyi Observational Station of Purple Mountain Observatory (PMO). The measurements are based on CCD imaging data taken in the Sloan Digital Sky Survey g, r and i bands with the Xuyi 1.04/1.20m Schmidt Telescope for the Xuyi Schmidt Telescope Photometric Survey of the Galactic Anti-center (XSTPS-GAC), the photometric part of the Digital Sky Survey of the Galactic Anti-center (DSS-GAC). The data were collected in more than 130 winter nights from 2009 to 2011. We find that the atmospheric extinction coefficients for the g, r and i bands are 0.70, 0.55 and 0.38 mag/airmass, respectively, based on observations taken in several photometric nights. The night sky brightness determined from images of good quality has median val- ues of 21.7, 20.8 and 20.0 mag/arcsec2 and reaches 22.1, 21.2 and 20.4 mag/arcsec2 under the best observing conditions for the g, r and i bands, respectively. The relatively large extinction coefficients compared with other good astronomical observing sites are mainly due to the relatively low elevation (i.e. 180 m) and high humidity of the Station.

The Damocloids are a group of unusual asteroids, recently enrolling a new member of 2010 EJ104. The dynamical evolution for the Damocloids may uncover a connection passage from the Main Belt, the Kuiper Belt and the scattered disk beyond. According to our simulations, two regions may be considered as possible origin of the Damocloids: the scattered disk, or a part of Oort cloud which will be perturbed to a transient region locating between 700 AU and 1000 AU. Based on the potential origin, the Damocloids can be classified into two types, with relation to their semi-major axes, and about 65.5% Damocloids is classified into type I which mainly originate from Oort cloud. Whether the Damocloids is inactive nuclei of Halley Family Comets may rely on their origin.

The issue of the influence of coronal holes (CHs) on coronal mass ejections (CMEs) in causing solar energetic particle (SEP) events is revisited. It is a continuation and extension of our previous work (Shen et al., 2006), in which no evident effect of CHs on CMEs in generating SEPs were found by statistically investigating 56 CME events. This result is consistent with the conclusion obtained by Kahler in 2004. In this paper, we extrapolate the coronal magnetic field, define CHs as the regions consisting of only open magnetic field lines and perform a similar analysis on this issue for totally 76 events by extending the study interval to the end of 2008. Three key parameters, CH proximity, CH area and CH relative position, are involved in the analysis. The new result confirms the previous conclusion that CHs did not show any evident effect on CMEs in causing SEP events.

The Smart Focal Planes (SmartFP) activity is a European Joint Research Activity funded to develop novel optical technologies for future large telescope instrumentation. In this paper, we will discuss the image slicer developments being carried out as part of this initiative. Image slicing technique s have many applications in the plans for instrumentation on Extremely Large Telescopes and will be central to the delivery of the science case. A study of a virtual "multi-object multi-ifu spectrograph and imager" (MOMSI) for a hypothetical OWL-class telescope reveals the need for focal plane splitting, deployable imagers and very small beam steering elements like deployable IFUs. The image slicer workpackage, lead from Durham University in collaboration with LFM Bremen, TNO Delft, UKATC Edinburgh, CRAL Lyon, LAM Marseille, Padua University and REFLEX Prague, is evaluating technologies for manufacturing micro optics in large numbers to enable multi-object integral field spectroscopy.

The Durham Microoptics Programme was established to develop key components to be used for integral field spectrographs for upcoming instrumentation projects, focussing on currently existing telescopes as well as on the next generation of ELTs. These activities include monolithic multi-optics machining and grinding, optical surface improvement using various post polishing techniques and replication of micro-optical components. While these developments have mostly slicer-type IFUs in mind, also new types of microlens arrays are in development for fiber based high contrast IFU systems.