We installed two sets of Astronomical Site Monitoring System(ASMS) at Lijiang Observatory(GMG), for the running of the 2.4-meter Lijiang optical telescope(LJT) and the 1.6-meter Multi-channel Photometric Survey Telescope (Mephisto). The Mephistro is under construction. ASMS has been running on robotic mode since 2017. The core instruments: Cloud Sensor, All-Sky Camera and Autonomous-DIMM that are developed by our group, together with the commercial Meteorological Station and Sky Quality Meter, are combined into the astronomical optical site monitoring system. The new Cloud Sensor's Cloud-Clear Relationship is presented for the first time, which is used to calculate the All-Sky cloud cover. We designed the Autonomous-DIMM located on a tower, with the same height as LJT. The seeing data have been observed for a full year. ASMS's data for the year 2019 are also analysed in detail, which are valuable to observers.

Cosmological simulations play an important role in the interpretation of astronomical data, in particular in comparing observed data to our theoretical expectations. However, to compare data with these simulations, the simulations in principle need to include gravity, magneto-hydrodyanmics, radiative transfer, etc. These ideal large-volume simulations (gravo-magneto-hydrodynamical) are incredibly computationally expensive which can cost tens of millions of CPU hours to run. In this paper, we propose a deep learning approach to map from the dark-matter-only simulation (computationally cheaper) to the galaxy distribution (from the much costlier cosmological simulation). The main challenge of this task is the high sparsity in the target galaxy distribution: space is mainly empty. We propose a cascade architecture composed of a classification filter followed by a regression procedure. We show that our result outperforms a state-of-the-art model used in the astronomical community, and provides a good trade-off between computational cost and prediction accuracy.

Measuring the sum of the three active neutrino masses, $M_\nu$, is one of the most important challenges in modern cosmology. Massive neutrinos imprint characteristic signatures on several cosmological observables in particular on the large-scale structure of the Universe. In order to maximize the information that can be retrieved from galaxy surveys, accurate theoretical predictions in the non-linear regime are needed. Currently, one way to achieve those predictions is by running cosmological numerical simulations. Unfortunately, producing those simulations requires high computational resources -- seven hundred CPU hours for each neutrino mass case. In this work, we propose a new method, based on a deep learning network (U-Net), to quickly generate simulations with massive neutrinos from standard $\Lambda$CDM simulations without neutrinos. We computed multiple relevant statistical measures of deep-learning generated simulations, and conclude that our method accurately reproduces the 3-dimensional spatial distribution of matter down to non-linear scales: $k < 0.7$ h/Mpc. Finally, our method allows us to generate massive neutrino simulations 10,000 times faster than the traditional methods.

The Quijote simulations are a set of 43100 full N-body simulations spanning more than 7000 cosmological models in the $\{\Omega_{\rm m}, \Omega_{\rm b}, h, n_s, \sigma_8, M_\nu, w \}$ hyperplane. At a single redshift the simulations contain more than 8.5 trillions of particles over a combined volume of 43100 $(h^{-1}{\rm Gpc})^3$. Billions of dark matter halos and cosmic voids have been identified in the simulations, whose runs required more than 35 million core hours. The Quijote simulations have been designed for two main purposes: 1) to quantify the information content on cosmological observables, and 2) to provide enough data to train machine learning algorithms. In this paper we describe the simulations and show a few of their applications. We also release the Petabyte of data generated, comprising hundreds of thousands of simulation snapshots at multiple redshifts, halo and void catalogs, together with millions of summary statistics such as power spectra, bispectra, correlation functions, marked power spectra, and estimated probability density functions.

We present Period-Luminosity and Period-Luminosity-Color relations at maximum-light for Mira variables in the Magellanic Clouds using time-series data from the Optical Gravitational Lensing Experiment (OGLE-III) and {\it Gaia} data release 2. The maximum-light relations exhibit a scatter typically up to $\sim 30\%$ smaller than their mean-light counterparts. The apparent magnitudes of Oxygen-rich Miras at maximum-light display significantly smaller cycle-to-cycle variations than at minimum-light. High-precision photometric data for Kepler Mira candidates also exhibit stable magnitude variations at the brightest epochs while their multi-epoch spectra display strong Balmer emission lines and weak molecular absorption at maximum-light. The stability of maximum-light magnitudes for Miras possibly occurs due to the decrease in the sensitivity to molecular bands at their warmest phase. At near-infrared wavelengths, the Period-Luminosity relations of Miras display similar dispersion at mean and maximum-light with limited time-series data in the Magellanic Clouds. A kink in the Oxygen-rich Mira Period-Luminosity relations is found at 300 days in the $VI$-bands which shifts to longer-periods ($\sim 350$~days) at near-infrared wavelengths. Oxygen-rich Mira Period-Luminosity relations at maximum-light provide a relative distance modulus, $\Delta \mu = 0.48\pm0.08$~mag, between the Magellanic Clouds with a smaller statistical uncertainty than the mean-light relations. The maximum-light properties of Miras can be very useful for stellar atmosphere modeling and distance scale studies provided their stability and the universality can be established in other stellar environments in the era of extremely large telescopes.

Matter evolved under influence of gravity from minuscule density fluctuations. Non-perturbative structure formed hierarchically over all scales, and developed non-Gaussian features in the Universe, known as the Cosmic Web. To fully understand the structure formation of the Universe is one of the holy grails of modern astrophysics. Astrophysicists survey large volumes of the Universe and employ a large ensemble of computer simulations to compare with the observed data in order to extract the full information of our own Universe. However, to evolve trillions of galaxies over billions of years even with the simplest physics is a daunting task. We build a deep neural network, the Deep Density Displacement Model (hereafter D$^3$M), to predict the non-linear structure formation of the Universe from simple linear perturbation theory. Our extensive analysis, demonstrates that D$^3$M outperforms the second order perturbation theory (hereafter 2LPT), the commonly used fast approximate simulation method, in point-wise comparison, 2-point correlation, and 3-point correlation. We also show that D$^3$M is able to accurately extrapolate far beyond its training data, and predict structure formation for significantly different cosmological parameters. Our study proves, for the first time, that deep learning is a practical and accurate alternative to approximate simulations of the gravitational structure formation of the Universe.

Lijiang 2.4-meter Telescope(LJT), the largest common-purpose optical telescope in China, has been applied to the world-wide astronomers since 2008. It is located at Gaomeigu site, Lijiang Observatory(LJO), the southwest of China. The site has very good observational conditions. Since 10-year operation, several instruments have been equipped on the LJT. Astronomers can perform both photometric and spectral observations. The main scientific goals of LJT include photometric and spectral evolution of supernova, reverberation mapping of active galactic nucleus, physical properties of binary star and near-earth object(comet and asteroid), identification of exoplanet, and all kinds of transients. Until now, the masses of 41 high accretion rate black holes have been measured, and more than 168 supernova have been identified by the LJT. More than 190 papers related to the LJT have been published. In this paper, the general observation condition of the Gaomeigu site is introduced at first. Then, the LJT structure is described in detail, including the optical, mechanical, motion and control system. The specification of all the instruments, and some detailed parameters of the YFOSC is also presented. Finally, some important scientific results and future expectations are summarized.

One of the most promising ways to observe the Universe is by detecting the 21cm emission from cosmic neutral hydrogen (HI) through radio-telescopes. Those observations can shed light on fundamental astrophysical questions only if accurate theoretical predictions are available. In order to maximize the scientific return of these surveys, those predictions need to include different observables and be precise on non-linear scales. Currently, one of the best ways to achieve this is via cosmological hydrodynamic simulations; however, the computational cost of these simulations is high -- tens of millions of CPU hours. In this work, we use Wasserstein Generative Adversarial Networks (WGANs) to generate new high-resolution ($35~h^{-1}{\rm kpc}$) 3D realizations of cosmic HI at $z=5$. We do so by sampling from a 100-dimension manifold, learned by the generator, that characterizes the fully non-linear abundance and clustering of cosmic HI from the state-of-the-art simulation IllustrisTNG. We show that different statistical properties of the produced samples -- 1D PDF, power spectrum, bispectrum, and void size function -- match very well those of IllustrisTNG, and outperform state-of-the-art models such as Halo Occupation Distributions (HODs). Our WGAN samples reproduce the abundance of HI across 9 orders of magnitude, from the Ly$\alpha$ forest to Damped Lyman Absorbers. WGAN can produce new samples orders of magnitude faster than hydrodynamic simulations.

Cosmological surveys aim at answering fundamental questions about our Universe, including the nature of dark matter or the reason of unexpected accelerated expansion of the Universe. In order to answer these questions, two important ingredients are needed: 1) data from observations and 2) a theoretical model that allows fast comparison between observation and theory. Most of the cosmological surveys observe galaxies, which are very difficult to model theoretically due to the complicated physics involved in their formation and evolution; modeling realistic galaxies over cosmological volumes requires running computationally expensive hydrodynamic simulations that can cost millions of CPU hours. In this paper, we propose to use deep learning to establish a mapping between the 3D galaxy distribution in hydrodynamic simulations and its underlying dark matter distribution. One of the major challenges in this pursuit is the very high sparsity in the predicted galaxy distribution. To this end, we develop a two-phase convolutional neural network architecture to generate fast galaxy catalogues, and compare our results against a standard cosmological technique. We find that our proposed approach either outperforms or is competitive with traditional cosmological techniques. Compared to the common methods used in cosmology, our approach also provides a nice trade-off between time-consumption (comparable to fastest benchmark in the literature) and the quality and accuracy of the predicted simulation. In combination with current and upcoming data from cosmological observations, our method has the potential to answer fundamental questions about our Universe with the highest accuracy.

Galaxy redshift surveys, such as 2dF, SDSS, 6df, GAMA and VIPERS, have shown that the spatial distribution of matter forms a rich web, known as the cosmic web. The majority of galaxy survey analyses measure the amplitude of galaxy clustering as a function of scale, ignoring information beyond a small number of summary statistics. Since the matter density field becomes highly non-Gaussian as structure evolves under gravity, we expect other statistical descriptions of the field to provide us with additional information. One way to study the non-Gaussianity is to study filaments, which evolve non-linearly from the initial density fluctuations produced in the primordial Universe. In our study, we report the first detection of CMB (Cosmic Microwave Background) lensing by filaments and we apply a null test to confirm our detection. Furthermore, we propose a phenomenological model to interpret the detected signal and we measure how filaments trace the matter distribution on large scales through filament bias, which we measure to be around 1.5. Our study provides a new scope to understand the environmental dependence of galaxy formation. In the future, the joint analysis of lensing and Sunyaev-Zel'dovich observations might reveal the properties of `missing baryons', the vast majority of the gas which resides in the intergalactic medium and has so far evaded most observations.

Previous studies have shown the filamentary structures in the cosmic web influence the alignments of nearby galaxies. We study this effect in the LOWZ sample of the Sloan Digital Sky Survey using the "Cosmic Web Reconstruction" filament catalogue. We find that LOWZ galaxies exhibit a small but statistically significant alignment in the direction parallel to the orientation of nearby filaments. This effect is detectable even in the absence of nearby galaxy clusters, which suggests it is an effect from the matter distribution in the filament. A nonparametric regression model suggests that the alignment effect with filaments extends over separations of 30-40 Mpc. We find that galaxies that are bright and early-forming align more strongly with the directions of nearby filaments than those that are faint and late-forming; however, trends with stellar mass are less statistically significant, within the narrow range of stellar mass of this sample.

Deep learning is a promising tool to determine the physical model that describes our universe. To handle the considerable computational cost of this problem, we present CosmoFlow: a highly scalable deep learning application built on top of the TensorFlow framework. CosmoFlow uses efficient implementations of 3D convolution and pooling primitives, together with improvements in threading for many element-wise operations, to improve training performance on Intel(C) Xeon Phi(TM) processors. We also utilize the Cray PE Machine Learning Plugin for efficient scaling to multiple nodes. We demonstrate fully synchronous data-parallel training on 8192 nodes of Cori with 77% parallel efficiency, achieving 3.5 Pflop/s sustained performance. To our knowledge, this is the first large-scale science application of the TensorFlow framework at supercomputer scale with fully-synchronous training. These enhancements enable us to process large 3D dark matter distribution and predict the cosmological parameters $\Omega_M$, $\sigma_8$ and n$_s$ with unprecedented accuracy.

With growing data from ongoing and future supernova surveys it is possible to empirically quantify the shapes of SNIa light curves in more detail, and to quantitatively relate the shape parameters with the intrinsic properties of SNIa. Building such relationship is critical in controlling systematic errors associated with supernova cosmology. Based on a collection of well-observed SNIa samples accumulated in the past years, we construct an empirical SNIa light curve model using a statistical method called the functional principal component analysis (FPCA) for sparse and irregularly sampled functional data. Using this method, the entire light curve of an SNIa is represented by a linear combination of principal component functions, and the SNIa is represented by a few numbers called principal component scores. These scores are used to establish relations between light curve shapes and physical quantities such as intrinsic color, interstellar dust reddening, spectral line strength, and spectral classes. These relations allow for descriptions of some critical physical quantities based purely on light curve shape parameters. Our study shows that some important spectral feature information is being encoded in the broad band light curves, for instance, we find that the light curve shapes are correlated with the velocity and velocity gradient of the Si II $\lambda$6355 line. This is important for supernova surveys, e.g., LSST and WFIRST. Moreover, the FPCA light curve model is used to construct the entire light curve shape, which in turn is used in a functional linear form to adjust intrinsic luminosity when fitting distance models.

We report here the results of searching for inelastic scattering of dark matter (initial and final state dark matter particles differ by a small mass splitting) with nucleon with the first 79.6-day of PandaX-II data (Run 9). We set the upper limits for the spin independent WIMP-nucleon scattering cross section up to a mass splitting of 300 keV/c$^2$ at two benchmark dark matter masses of 1 and 10 TeV/c$^2$.

We report a new search of weakly interacting massive particles (WIMPs) using the combined low background data sets in 2016 and 2017 from the PandaX-II experiment in China. The latest data set contains a new exposure of 77.1 live day, with the background reduced to a level of 0.8$\times10^{-3}$ evt/kg/day, improved by a factor of 2.5 in comparison to the previous run in 2016. No excess events were found above the expected background. With a total exposure of 5.4$\times10^4$ kg day, the most stringent upper limit on spin-independent WIMP-nucleon cross section was set for a WIMP with mass larger than 100 GeV/c$^2$, with the lowest exclusion at 8.6$\times10^{-47}$ cm$^2$ at 40 GeV/c$^2$.

We study the near-infrared properties of 690 Mira candidates in the central region of the Large Magellanic Cloud, based on time-series observations at JHKs. We use densely-sampled I-band observations from the OGLE project to generate template light curves in the near infrared and derive robust mean magnitudes at those wavelengths. We obtain near-infrared Period-Luminosity relations for Oxygen-rich Miras with a scatter as low as 0.12 mag at Ks. We study the Period-Luminosity-Color relations and the color excesses of Carbon-rich Miras, which show evidence for a substantially different reddening law.

We present the discovery of 1847 Mira candidates in the Local Group galaxy M33 using a novel semi-parametric periodogram technique coupled with a Random Forest classifier. The algorithms were applied to ~2.4x10^5 I-band light curves previously obtained by the M33 Synoptic Stellar Survey. We derive preliminary Period-Luminosity relations at optical, near- & mid-infrared wavelengths and compare them to the corresponding relations in the Large Magellanic Cloud.

In a two-component jet model, the emissions are the sum of the core and extended emissions: $S^{\rm ob}=S_{\rm core}^{\rm ob}+S_{\rm ext}^{\rm ob}$, with the core emissions, $S_{\rm core}^{\rm ob}= f S_{\rm ext}^{\rm ob}\delta^{q}$, being a function of the Doppler factor, $\delta$, the extended emission, $S_{\rm ext}^{\rm ob}$, jet type dependent factor, $q$, and the ratio of the core to the extended emissions in the comoving frame, $f$. The $f$ is an unobservable but important parameter. Following our previous work, we collect 65 blazars with available Doppler factor, $\delta$, superluminal velocity, $\beta_{app}$, and core-dominance parameter, $R$, calculate the ratio, $f$, and peform statistical analyses. We find that the ratio, $f$, in BL Lacs is on average larger than that in FSRQs. We suggest that the difference of the ratio $f$ between FSRQs and BL Lacs is one of the possible reasons that cause the difference of other observed properties between them. We also find some significant correlations between $\log f$ and other parameters, including intrinsic (de-beamed) peak frequency, $\log \nu _{\rm p}^{\rm in}$, intrinsic polarization, $\log P^{\rm in}$, and core-dominance parameter, $\log R$, for the whole sample. In addition, we show that the ratio, $f$, can be estimated by $R$.

We develop a non-linear semi-parametric Gaussian process model to estimate periods of Miras with sparsely-sampled light curves. The model uses a sinusoidal basis for the periodic variation and a Gaussian process for the stochastic changes. We use maximum likelihood to estimate the period and the parameters of the Gaussian process, while integrating out the effects of other nuisance parameters in the model with respect to a suitable prior distribution obtained from earlier studies. Since the likelihood is highly multimodal for period, we implement a hybrid method that applies the quasi-Newton algorithm for Gaussian process parameters and search the period/frequency parameter over a dense grid. A large-scale, high-fidelity simulation is conducted to mimic the sampling quality of Mira light curves obtained by the M33 Synoptic Stellar Survey. The simulated data set is publicly available and can serve as a testbed for future evaluation of different period estimation methods. The semi-parametric model outperforms an existing algorithm on this simulated test data set as measured by period recovery rate and quality of the resulting Period-Luminosity relations.

We investigate anew the distribution of absolute carbon abundance, $A$(C) $= \log\,\epsilon $(C), for carbon-enhanced metal-poor (CEMP) stars in the halo of the Milky Way, based on high-resolution spectroscopic data for a total sample of 305 CEMP stars. The sample includes 147 CEMP-$s$ (and CEMP-r/s) stars, 127 CEMP-no stars, and 31 CEMP stars that are unclassified, based on the currently employed [Ba/Fe] criterion. We confirm previous claims that the distribution of $A$(C) for CEMP stars is (at least) bimodal, with newly determined peaks centered on $A$(C)$=7.96$ (the high-C region) and $A$(C)$ =6.28$ (the low-C region). A very high fraction of CEMP-$s$ (and CEMP-r/s) stars belong to the high-C region, while the great majority of CEMP-no stars reside in the low-C region. However, there exists complexity in the morphology of the $A$(C)-[Fe/H] space for the CEMP-no stars, a first indication that more than one class of first-generation stellar progenitors may be required to account for their observed abundances. The two groups of CEMP-no stars we identify exhibit clearly different locations in the $A$(Na)-$A$(C) and $A$(Mg)-$A$(C) spaces, also suggesting multiple progenitors. The clear distinction in $A$(C) between the CEMP-$s$ (and CEMP-$r/s$) stars and the CEMP-no stars appears to be $as\ successful$, and $likely\ more\ astrophysically\ fundamental$, for the separation of these sub-classes as the previously recommended criterion based on [Ba/Fe] (and [Ba/Eu]) abundance ratios. This result opens the window for its application to present and future large-scale low- and medium-resolution spectroscopic surveys.

We demonstrate a new method to constrain gravity on the largest cosmological scales by combining measurements of cosmic microwave background (CMB) lensing and the galaxy velocity field. $E_G$ is a statistic, constructed from a gravitational lensing tracer and a measure of velocities such as redshift-space distortions (RSD), that can discriminate between gravity models while being independent of clustering bias and $\sigma_8$. While traditionally, the lensing field for $E_G$ has been probed through galaxy lensing, CMB lensing has been proposed as a more robust tracer of the lensing field for $E_G$ at higher redshifts while avoiding intrinsic alignments. We perform the largest-scale measurement of $E_G$ ever, up to 150 Mpc/$h$, by cross-correlating the Planck CMB lensing map with the Sloan Digital Sky Survey III (SDSS-III) CMASS galaxy sample and combining this with our measurement of the CMASS auto-power spectrum and the RSD parameter $\beta$. We report $E_G(z=0.57)=0.243\pm0.060$ (stat) $\pm0.013$ (sys), a measurement in tension with the general relativity prediction at a level of 2.6$\sigma$. Note that our $E_G$ measurement deviates from GR only at scales greater than 80 Mpc/$h$, scales which have not been probed by previous $E_G$ tests. Upcoming surveys, which will provide an order-of-magnitude reduction in statistical errors, can significantly constrain alternative gravity models when combined with better control of systematics.

Optical and ultraviolet observations for the nearby type II-plateau supernova (SN IIP) 2013am in the nearby spiral galaxy M65 are presented in this paper. The early spectra are characterized by relatively narrow P-Cygni features, with ejecta velocities much lower than observed in normal SNe IIP (i.e., $\sim$2000 km s$^{-1}$ vs. $\sim$5000 km $^{-1}$ in the middle of the plateau phase). Moreover, prominent Ca II absorptions are also detected in SN 2013am at relatively early phases. These spectral features are reminiscent of those seen in the low-velocity and low-luminosity SN IIP 2005cs. However, SN 2013am exhibits different photometric properties, having shorter plateau phases and brighter light-curve tails if compared to SN 2005cs. Adopting $R_{V}$=3.1 and a mean value of total reddening derived from the photometric and spectroscopic methods(i.e., $E(B-V)=0.55\pm$0.19 mag), we find that SN 2013am may have reached an absolute $V$-band peak magnitude of $-15.83\pm0.71$ mag, and produced a $^{56}$Ni mass of $0.016^{+0.010}_{-0.006} solar masses in the explosion. These parameters are close to those derived for SN 2008in and SN 2009N which have been regarded as "gap-filler" objects linking the faint SNe IIP to the normal ones. This indicates that some low-velocity SNe IIP may not necessarily result from the low-energetic explosions, and the low expansion velocities could be due to a lower metallicity of the progenitor stars, a larger envelope mass ejected in the explosion, or that is was observed at an angle that is away from the polar direction.

Recently a covariant entropy conjecture has been proposed for dynamical horizons. We apply this conjecture to concordance cosmological models, namely, those cosmological models filled with perfect fluids, in the presence of a positive cosmological constant. As a result, we find this conjecture has a severe constraint power. Not only does this conjecture rule out those cosmological models disfavored by the anthropic principle, but also it imposes an upper bound $10^{-60}$ on the cosmological constant for our own universe, which thus provides an alternative macroscopic perspective for understanding the long-standing cosmological constant problem.

As a compelling pattern for the holographic principle, our covariant entropy bound conjecture is proposed for more general dynamical horizons. Then we apply our conjecture to $\Lambda$CDM cosmological models, where we find it imposes a novel upper bound $10^{-90}$ on the cosmological constant for our own universe by taking into account the dominant entropy contribution from super-massive black holes, which thus provides an alternative macroscopic perspective to understand the longstanding cosmological constant problem. As an intriguing implication of this conjecture, we also discuss the possible profound relation between the present cosmological constant, the origin of mass, and the anthropic principle.

The generalized second law of thermodynamics for a system containing a black hole dynamical horizon is proposed in a covariant way. Its validity is also tested in case of adiabatically collapsing thick light shells.

Taking into account that Sandage-Loeb test is unique in its coverage of the redshift desert and available in the near future, we explore the cosmic time evolution behavior of the source redshift for holographic dark energy model, an important competing cosmological model. As a result, we find that Sandage-Loeb test can provide a extremely strong bound on $\Omega^0_m$, while its constraint on another dimensionless parameter $\lambda$ is weak. In addition, it is proposed here for the first time that we can also constrain various cosmological model by measuring the value of $z_{max}$ at which the peak of redshift velocity occurs. Combining this new proposed method with the traditional Sandage-Loeb test, we should be able to provide a better constraint on $\lambda$, at least from the theoretical perspective.

We here propose a covariant entropy conjecture on cosmological dynamical horizon. After the formulation of our conjecture, we test its validity in adiabatically expanding universes with open, flat and closed spatial geometry, where our conjecture can also be viewed as a cosmological version of the generalized second law of thermodynamics in some sense.