We developed a 3-D elastic Boundary Element Method (BEM) computer code, called AstroSeis, to model seismic wavefields in a body with an arbitrary shape, such as an asteroid. Besides the AstroSeis can handle arbitrary surface topography, it can deal with a liquid core in an asteroid model. Both the solid and liquid domains are homogenous in our current code. For seismic sources, we can use single forces or moment tensors. The AstroSeis is implemented in the frequency domain and the frequency-dependent Q can be readily incorporated. The code is in MATLAB and it is straightforward to set up the model to run our code. The frequency-domain calculation is advantageous to study the long-term elastic response of a celestial body due to a cyclic force such as the tidal force with no numerical dispersion issue suffered by many other methods requiring volume meshing. Our AstroSeis has been benchmarked with other methods such as normal modes summation and the direct solution method (DSM). This open-source AstroSeis will be a useful tool to study the interior and surface processes of asteroids.

The radial acceleration relation (RAR) in galaxies describes a tight empirical scaling law between the total acceleration $g_\mathrm{tot}(r)=GM_\mathrm{tot}(<r)/r^2$ observed in galaxies and that expected from their baryonic mass $g_\mathrm{bar}(r)=GM_\mathrm{bar}(<r)/r^2$, with a characteristic acceleration scale of $g_\dagger\simeq 1.2\times 10^{-10}$ms$^{-2}$. Here, we examine if such a correlation exists in galaxy clusters using weak-lensing, strong-lensing, and X-ray data sets available for 20 high-mass clusters targeted by the CLASH survey. By combining our CLASH data with stellar mass estimates for the brightest cluster galaxies (BCGs) and accounting for the stellar baryonic component in clusters, we determine, for the first time, an RAR on BCG--cluster scales. The resulting RAR is well described by a tight power-law relation, $g_\mathrm{tot}\propto g_\mathrm{bar}^{0.51^{+0.04}_{-0.05}}$, with lognormal intrinsic scatter of $14.7^{+2.9}_{-2.8}\%$. The slope is consistent with the low acceleration limit of the RAR in galaxies, $g_\mathrm{tot}=\sqrt{g_\dagger\,g_\mathrm{bar}}$, whereas the intercept implies a much higher acceleration scale of $g_\ddagger = (2.02\pm0.11)\times 10^{-9}$ms$^{-2}$, indicating that there is no universal RAR that holds on all scales from galaxies to clusters. We find that the observed RAR in CLASH clusters is consistent with predictions from a semi-analytical model developed in the standard $\Lambda$CDM framework. Our results also predict the presence of a baryonic Faber--Jackson relation ($\sigma_v^4\propto M_\mathrm{bar}$) on cluster scales.

The Large Optical/infrared Telescope (LOT) is a ground-based 12m diameter optical/infrared telescope which is proposed to be built in the western part of China in the next decade. Based on satellite remote sensing data, along with geographical, logistical and political considerations, three candidate sites were chosen for ground-based astronomical performance monitoring. These sites include: Ali in Tibet, Daocheng in Sichuan, and Muztagh Ata in Xinjiang. Up until now, all three sites have continuously collected data for two years. In this paper, we will introduce this site testing campaign, and present its monitoring results obtained during the period between March 2017 and March 2019.

The Gamma-Ray Integrated Detectors (GRID) is a space mission concept dedicated to monitoring the transient gamma-ray sky in the energy range from 10 keV to 2 MeV using scintillation detectors onboard CubeSats in low Earth orbits. The primary targets of GRID are the gamma-ray bursts (GRBs) in the local universe. The scientific goal of GRID is, in synergy with ground-based gravitational wave (GW) detectors such as LIGO and VIRGO, to accumulate a sample of GRBs associated with the merger of two compact stars and study jets and related physics of those objects. It also involves observing and studying other gamma-ray transients such as long GRBs, soft gamma-ray repeaters, terrestrial gamma-ray flashes, and solar flares. With multiple CubeSats in various orbits, GRID is unaffected by the Earth occultation and serves as a full-time and all-sky monitor. Assuming a horizon of 200 Mpc for ground-based GW detectors, we expect to see a few associated GW-GRB events per year. With about 10 CubeSats in operation, GRID is capable of localizing a faint GRB like 170817A with a 90% error radius of about 10 degrees, through triangulation and flux modulation. GRID is proposed and developed by students, with considerable contribution from undergraduate students, and will remain operated as a student project in the future. The current GRID collaboration involves more than 20 institutes and keeps growing. On August 29th, the first GRID detector onboard a CubeSat was launched into a Sun-synchronous orbit and is currently under test.

Recently, from the new Spitzer Photometry and Accurate Rotation Curves (SPARC) data, McGaugh et al. (2016) reported a tight Radial Acceleration Relation (RAR) between the observed total acceleration and the acceleration produced by baryons in spiral galaxies. The relation can be fitted by different functions. However, these functions can be discerned if we express the data in the form of halo acceleration relation (HAR). The data reveals a maximum in the halo acceleration. We examined NFW (cusp) and Burkert (core) profiles in the context of dark matter and different parameter families of the interpolating function in the framework of Modified Newtonian Dynamics (MOND).

Tidal force plays an important role in the evolution of the planet-moon system. The tidal force of a moon can excite seismic waves in the planet it is orbiting. A tidal-seismic resonance is expected when a tidal force frequency matches a free-oscillation frequency of the planet. Here we show that when the moon is close to the planet, the tidal-seismic resonance can cause large-amplitude seismic waves, which can change the shape of the planet and in turn exert a negative torque on the moon to cause it to fall rapidly toward the planet. We postulate that the tidal-seismic resonance may be an important mechanism which can accelerate planet accretion process. On the other hand, tidal-seismic resonance effect can also be used to interrogate planet interior by long term tracking of the orbital change of the moon.

A study on cosmogenic activation in germanium was carried out to evaluate the cosmogenic background level of natural and $^{70}$Ge depleted germanium detectors. The production rates of long-lived radionuclides were calculated with Geant4 and CRY. Results were validated by comparing the simulated and experimental spectra of CDEX-1B detector. Based on the validated codes, the cosmogenic background level was predicted for further tonne-scale CDEX experiment. The suppression of cosmogenic background level could be achieved by underground germanium crystal growth and high-purity germanium detector fabrication to reach the sensitivity requirement for direct detection of dark matter. With the low cosmogenic background, new physics channels, such as solar neutrino research and neutrinoless double-beta decay experiments, were opened and the corresponding simulations and evaluations were carried out.

The Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) is the largest existing spectroscopic survey telescope, having 32 scientific charge-coupled-device (CCD) cameras for receiving spectra. Stability and automation of the camera-control software are essential, but cannot be provided by the existing system. The Remote Telescope System 2nd Version (RTS2) is an open-source and automatic observatory-control system. However, all previous RTS2 applications have concerned small telescopes. This paper focuses on implementation of an RTS2-based camera-control system for the 32 CCDs of LAMOST. A virtual camera module inherited from the RTS2 camera module is built as a device component working on the RTS2 framework. To improve the controllability and robustness, a virtualized layer is designed using the master-slave software paradigm, and the virtual camera module is mapped to the 32 real cameras of LAMOST. The new system is deployed in the actual environment and experimentally tested. Finally, multiple observations are conducted using this new RTS2-framework-based control system. The new camera-control system is found to satisfy the requirements for LAMOST automatic camera control. This is the first time that RTS2 has been applied to a large telescope, and provides a referential solution for full RTS2 introduction to the LAMOST observatory control system.

We study the Mass Discrepancy-Acceleration Relation (MDAR) of 57 elliptical galaxies by their Einstein rings from the Sloan Lens ACS Survey (SLACS). The mass discrepancy between the lensing mass and the baryonic mass derived from population synthesis is larger when the acceleration of the elliptical galaxy lenses is smaller. The MDAR is also related to surface mass density discrepancy. At the Einstein ring, these lenses belong to high-surface-mass density galaxies. Similarly, we find that the discrepancy between the lensing and stellar surface mass density is small. It is consistent with the recent discovery of dynamical surface mass density discrepancy in disk galaxies where the discrepancy is smaller when surface density is larger. We also find relativistic modified Newtonian dynamics (MOND) can naturally explain the MDAR and surface mass density discrepancy in 57 Einstein rings. Moreover, the lensing mass, the dynamical mass and the stellar mass of these galaxies are consistent with each other in relativistic MOND.

The dynamics of an elliptical galaxy within a couple of effective radii can be probed effectively by stars. However, at larger distances planetary nebulae (PNe) replace stars as the tracer of the dynamics. Making use of the motion of PNe, Romanowsky et al. (2003) measured the dynamics of three luminous elliptical galaxies (NGC821, NGC3379, and NGC4494) at large distances from the galactic center. They found that little dark matter is needed up to 6 effective radii. Milgrom & Sanders (2003) showed that this result can be understood in the framework of MOdified Newtonian Dynamics (MOND). As more data are available in the past decade, we revisit this problem. We combine PNe data (up to 6{8 effective radii) and stellar data from SAURON of 7 elliptical galaxies, including those 3 galaxies in Romanowsky et al. (2003) with updated data and 4 other galaxies which have not been analyzed before. We conclude that the dynamics of these galaxies can be well explained by MOND.

The Large sky Area Multi-Object Spectroscopic Telescope (LAMOST) General Survey is a spectroscopic survey that will eventually cover approximately half of the celestial sphere and collect 10 million spectra of stars, galaxies and QSOs. Objects both in the pilot survey and the first year general survey are included in the LAMOST First Data Release (DR1). The pilot survey started in October 2011 and ended in June 2012, and the data have been released to the public as the LAMOST Pilot Data Release in August 2012. The general survey started in September 2012, and completed its first year of operation in June 2013. The LAMOST DR1 includes a total of 1202 plates containing 2,955,336 spectra, of which 1,790,879 spectra have observed signal-to-noise S/N >10. All data with S/N>2 are formally released as LAMOST DR1 under the LAMOST data policy. This data release contains a total of 2,204,696 spectra, of which 1,944,329 are stellar spectra, 12,082 are galaxy spectra and 5,017 are quasars. The DR1 includes not only spectra, but also three stellar catalogues with measured parameters: AFGK-type stars with high quality spectra (1,061,918 entries), A-type stars (100,073 entries), and M stars (121,522 entries). This paper introduces the survey design, the observational and instrumental limitations, data reduction and analysis, and some caveats. Description of the FITS structure of spectral files and parameter catalogues is also provided.

The BICEP2 experiment confirms the existence of primordial gravitational wave with the tensor-to-scalar ratio $r=0$ ruled out at $7\sigma$ level. The consistency of this large value of $r$ with the {\em Planck} data requires a large negative running $n'_s$ of the scalar spectral index. Herein we propose two types of the single field inflation models with simple potentials to study the possibility of the consistency of the models with the BICEP2 and {\em Planck} observations. One type of model suggested herein is realized in the supergravity model building. These models fail to provide the needed $n'_s$ even though both can fit the tensor-to-scalar ratio and spectral index.

Time delay in galaxy gravitational lensing systems has been used to determine the value of Hubble constant. As in other dynamical phenomena at the scale of galaxy, dark matter is often invoked in gravitational lensing to account for the "missing mass" (the apparent discrepancy between the dynamical mass and the luminous mass). Alternatively, modified gravity can be used to explain the discrepancy. In this paper we adopt the Tensor-Vector-Scalar gravity (TeVeS), a relativistic version of MOdified Newtonian Dynamics (MOND), to study gravitational lensing phenomena and derive the formulae needed to evaluate the Hubble constant. We test our method on quasar lensing by elliptical galaxies in the literature. We focus on double-image systems with time delay measurement. Three candidates are suitable for our study: HE 2149-2745, FBQ J0951+2635 and SBS 0909+532. The Hubble constant obtained is consistent with the value used in fitting the CMB result in neutrino cosmological model.

We present new quasars discovered in the vicinity of the Andromeda and Triangulum galaxies with the LAMOST during the 2010 and 2011 observational seasons. Quasar candidates are selected based on the available SDSS, KPNO 4 m telescope, XSTPS optical, and WISE near infrared photometric data. We present 509 new quasars discovered in a stripe of ~135 sq. deg from M31 to M33 along the Giant Stellar Stream in the 2011 pilot survey datasets, and also 17 new quasars discovered in an area of ~100 sq. deg that covers the central region and the southeastern halo of M31 in the 2010 commissioning datasets. These 526 new quasars have i magnitudes ranging from 15.5 to 20.0, redshifts from 0.1 to 3.2. They represent a significant increase of the number of identified quasars in the vicinity of M31 and M33. There are now 26, 62 and 139 known quasars in this region of the sky with i magnitudes brighter than 17.0, 17.5 and 18.0 respectively, of which 5, 20 and 75 are newly-discovered. These bright quasars provide an invaluable collection with which to probe the kinematics and chemistry of the ISM/IGM in the Local Group of galaxies. A total of 93 quasars are now known with locations within 2.5 deg of M31, of which 73 are newly discovered. Tens of quasars are now known to be located behind the Giant Stellar Stream, and hundreds behind the extended halo and its associated substructures of M31. The much enlarged sample of known quasars in the vicinity of M31 and M33 can potentially be utilized to construct a perfect astrometric reference frame to measure the minute PMs of M31 and M33, along with the PMs of substructures associated with the Local Group of galaxies. Those PMs are some of the most fundamental properties of the Local Group.

Using strong lensing data Milgrom's MOdified Newtonian Dynamics (MOND) or its covariant TeVeS (Tensor-Vector-Scalar Theory) is being examined here as an alternative to the conventional $\Lambda$CDM paradigm. We examine 10 double-image gravitational lensing systems, in which the lens masses have been estimated by stellar population synthesis models. While mild deviations exist, we do not find out that strong cases for outliers to the TeVeS theory.

Modified Newtonian Dynamics (MOND) and its relativistic version - TeVeS offer us an alternative perspective to understand the universe without the demand of the elusive cold dark matter. This MONDian paradigm is not only competitive with the conventional CDM in a large range of scales, but also even more successful in the galactic scale. Recently, by studying 6 lensing systems, Ferreras et al. (2008) claimed that MOND still needs dark matter even in galactic scales. When we study the same systems, however, we yield an opposite conclusion. In this contribution, we report our result and conclude that MOND does not need dark matter in galactic lensing systems. Furthermore, we extend our study to 22 SLACS (Sloan Lens ACS Survey) lenses, and obtain the same conclusion as well, i.e., no dark matter is needed in elliptical galaxies.

Since Bekenstein's (2004) creation of his Tensor-Vector-Scalar theory (TeVeS), the Modified Newtonian dynamics (MOND) paradigm has been redeemed from the embarrassment of lacking a relativistic version. One primary success of TeVeS is that it provides an enhancement of gravitational lensing, which could not be achieved by other MONDian theories. Following Bekenstein's work, we investigate the phenomena of gravitational lensing including deflection angles, lens equations and time delay. We find that the deflection angle would maintain its value while the distance of closest approach vary in the MOND regime, this coincides with the conclusion of Mortlock and Turner's (2001) intuitional approach. Moreover, the scalar field, which is introduced to enhance the deflection angle in TeVeS, contributes a negative effect on the potential time delay. Unfortunately this phenomenon is unmeasurable in lensing systems where we can only observe the time delay between two images for a given source. However, this measurable time delay offers another constraint on the mass ratio of the DM and MOND scenarios, which in general differs from that given by the deflection angle. In other words, for a lensing system, if two masses, m_gN and m_gM, are mutually alternatives for the deflection angles in their own paradigm, regarding the time delay they are in general in an exclusive relation.