With telescope apertures becoming larger and larger, the deployment of large-field telescopes is becoming increasingly popular. However, optical path calibration is necessary to ensure the image quality of large-field and large-diameter telescopes. In particular, focal plane attitude calibration is an essential optical path calibration technique that has a direct impact on image quality. In this paper, a focal plane attitude detection method using eight acquisition cameras is proposed based on the calibration requirements of the wide-field telescope, LAMOST. Comparison of simulation and experimental results shows that the detection accuracy of the proposed method can reach 30 arcsec. With additional testing and verification, this method could be used to facilitate regular focal plane attitude calibration for LAMOST as well as other large-field telescopes.

We present the Data Release 4&5 quasar catalog from the quasar survey by Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST), which includes quasars observed between September 2015 and June 2017. There are a total of 19,253 quasars identified by visual inspections of the spectra. Among them, 11,458 are independently discovered by LAMOST, in which 3296 were reported by SDSS DR12 and DR14 quasar catalog after our survey began, while the rest 8162 are new discoveries of LAMOST. We provide the emission line measurements for the Halpha, Hbeta, MgII and/or CIV for 18100 quasars. Since LAMOST does not have absolute flux calibration information, we obtain the monochromatic continuum luminosities by fitting the SDSS photometric data using the quasar spectra, and then estimate the black hole masses. The catalog and spectra for these quasars are available online. This is the third installment in the series of LAMOST quasar survey which has released spectra for totally ~43,000 quasars hitherto. There are 24,772 independently discovered quasars, 17,128 of which are newly discovered. In addition to the great supplement to the new quasar discoveries, LAMOST has also provided a large database (overlapped with SDSS) for investigating the quasar spectral variability and discovering unusual quasars, including changing-look quasars, with ongoing and upcoming large surveys.

Chebyshev coefficients of a coordinate representation can be used to form the corresponding velocity representation. One way is to directly apply them to the derivatives of Chebyshev polynomials, another is to compute from them the Chebyshev coefficients of the velocity representation. The advantages of the latter over the former ways are illustrated. Also, the approach of generating Chebyshev coefficients developed by Newhall (1989) is extended such that coordinate, velocity and acceleration are consistently treated. The resulting representations are all continuous.

The rare case of changing-look (CL) AGNs, with the appearance or disappearance of broad Balmer emission lines within a few years, challenges our understanding of the AGN unified model. We present a sample of 21 new CL AGNs at $0.08<z<0.58$, which doubles the number of such objects known to date. These new CL AGNs were discovered by several ways, from (1) repeat spectra in the SDSS, (2) repeat spectra in the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) and SDSS, and (3) photometric variability and new spectroscopic observations. We use the photometric data from surveys, including the SDSS imaging survey, the Pan-STARRS1, the DESI Legacy imaging survey, the Wide-field Infrared Survey Explorer (WISE), the Catalina Real-time Transient Survey, and the Palomar Transient Factory. The estimated upper limits of transition timescale of the CL AGNs in this sample spans from 0.9 to 13 years in the rest frame. The continuum flux in the optical and mid-infrared becomes brighter when the CL AGNs turn on, or vice versa. Variations of more than 0.2 mag in $W1$ band were detected in 15 CL AGNs during the transition. The optical and mid-infrared variability is not consistent with the scenario of variable obscuration in 10 CL AGNs at more than $3\sigma$ confidence level. We confirm a bluer-when-brighter trend in the optical. However, the mid-infrared WISE colors $W1-W2$ become redder when the objects become brighter in the $W1$ band, possibly due to a stronger hot dust contribution in the $W2$ band when the AGN activity becomes stronger. The physical mechanism of type transition is important for understanding the evolution of AGNs.

In this paper we illustrate the technique used by the New Vacuum Solar Telescope to increase the spatial resolution of two-dimensional (2D) solar spectroscopy observation involving two dimensions of space and one of wavelength. Without an image stabilizer at the NVST, a large scale wobble motion is present during the spatial scanning, whose instantaneous amplitude could reach up to 1.3" due to the earth's atmosphere and the precision of the telescope guiding system, and seriously decreases the spatial resolution of 2D spatial maps composed with the scanning spectra. We make the following effort to resolve this problem: the imaging system (e.g., the TiO-band) is used to record and detect the displacement vectors of solar image motion during the raster scan, in both the slit and scanning directions. The spectral data (e.g., the Ha line) which are originally obtained in time sequence are corrected and re-arranged in space according to those displacement vectors. Raster scans are carried out in several active regions with different seeing conditions (two rasters are illustrated in this paper). Given a certain spatial sample and temporal resolution, the spatial resolution of the composed 2D map could be close to that of the slit-jaw image. The resulting quality after correction is quantitatively evaluated with two methods. Two-dimensional physical quantity, such as the line-of-sight velocities in multi-layer of the solar atmosphere, is also inferred demonstrating the effect of this technique.

We present a new algorithm to estimate quasar photometric redshifts (photo-$z$s), by considering the asymmetries in the relative flux distributions of quasars. The relative flux models are built with multivariate Skew-t distributions in the multi-dimensional space of relative fluxes as a function of redshift and magnitude. For 151,392 quasars in the SDSS, we achieve a photo-$z$ accuracy, defined as the fraction of quasars with the difference between the photo-$z$ $z_p$ and the spectroscopic redshift $z_s$, $|\Delta z| = |z_s-z_p|/(1+z_s)$ within 0.1, of 74%. Combining the WISE W1 and W2 infrared data with the SDSS data, the photo-$z$ accuracy is enhanced to 87%. Using the Pan-STARRS1 or DECaLS photometry with WISE W1 and W2 data, the photo-$z$ accuracies are 79% and 72%, respectively. The prior probabilities as a function of magnitude for quasars, stars and galaxies are calculated respectively based on (1) the quasar luminosity function; (2) the Milky Way synthetic simulation with the Besan\c{c}on model; (3) the Bayesian Galaxy Photometric Redshift estimation. The relative fluxes of stars are obtained with the Padova isochrones, and the relative fluxes of galaxies are modeled through galaxy templates. We test our classification method to select quasars using the DECaLS $g$, $r$, $z$, and WISE W1 and W2 photometry. The quasar selection completeness is higher than 70% for a wide redshift range $0.5<z<4.5$, and a wide magnitude range $18<r<21.5$ mag. Our photo-$z$ regression and classification method has the potential to extend to future surveys. The photo-$z$ code will be publicly available.

We present a precise and complete procedure for processing spectral data observed by the 1-meter New Vacuum Solar Telescope (NVST). The procedure is suitable for both the sit-and-stare and raster-scan spectra. In this work, the geometric distortions of the spectra are firstly corrected for subsequent processes. Then, considering the temporal changes and the remnants of spectral lines in the flat-field, the original flat-field matrix is split into four independent components for ensuring a high precision flat-fielding correction, consisting of the continuum gradient matrix, slit non-uniform matrix, CCD dust matrix, and interference fringe matrix. Subsequently, the spectral line drifts and intensity fluctuations of the science data are further corrected. After precise reduction with this procedure, the measuring accuracies of the Doppler velocities for different spectral lines and of the oscillation curves of the chromosphere and photosphere are measured. The results show that the highest measuring accuracy of the Doppler velocity is within 100 ms-1, which indicates that the characteristics of the photosphere and chromosphere can be studied co-spatially and co-temporally with the reduced spectra of NVST.

Over short time intervals planetary ephemerides have been traditionally represented in analytical form as finite sums of periodic terms or sums of Poisson terms that are periodic terms with polynomial amplitudes. Nevertheless, this representation is not well adapted for the evolution of the planetary orbits in the solar system over million of years as they present drifts in their main frequencies, due to the chaotic nature of their dynamics. The aim of the present paper is to develop a numerical algorithm for slowly diffusing solutions of a perturbed integrable Hamiltonian system that will apply to the representation of the chaotic planetary motions with varying frequencies. By simple analytical considerations, we first argue that it is possible to recover exactly a single varying frequency. Then, a function basis involving time-dependent fundamental frequencies is formulated in a semi-analytical way. Finally, starting from a numerical solution, a recursive algorithm is used to numerically decompose the solution on the significant elements of the function basis. Simple examples show that this algorithm can be used to give compact representations of different types of slowly diffusing solutions. As a test example, we show how this algorithm can be successfully applied to obtain a very compact approximation of the La2004 solution of the orbital motion of the Earth over 40 Myr ([-35Myr,5Myr]). This example has been chosen as this solution is widely used for the reconstruction of the climates of the past.

We study the effect of newly emerged solar active regions (ARs) on the large-scale magnetic environment of pre-existing ARs (PEARs). We first present a theoretical approach to quantify the "interaction energy" between new ARs and PEARs as the difference between (i) the summed magnetic energies of their individual potential fields and (ii) the energy of their superposed potential fields. We expect that this interaction energy can, depending upon the relative arrangements of newly emerged and PEAR magnetic flux, indicate the existence of "topological" free magnetic energy in the global coronal field that is independent of any "internal" free magnetic energy due to coronal electric currents flowing within the newly emerged and PEAR flux systems. We then examine the interaction energy in two well-studied cases of flux emergence, but find that the predicted energetic perturbation is relatively small compared to energies released in large solar flares. Next, we present an observational study on the influence of the emergence of new ARs on flare statistics in PEARs, using NOAA's Solar Region Summary and GOES flare databases. As part of an effort to precisely determine the emergence time of ARs in a large event sample, we find that emergence in about half of these regions exhibits a two-stage behavior, with an initial gradual phase followed by a more rapid phase. Regarding flaring, we find that the emergence of new ARs is associated with a significant increase in the occurrence rate of X- and M-class flares in PEARs. This effect tends to be more significant when PEARs and new emerging ARs are closer. Given the relative weakness of the interaction energy, this effect suggests that perturbations in the large-scale magnetic field, such as topology changes invoked in the "breakout" model of coronal mass ejections, might play a significant role in the occurrence of some flares.

For most hierarchical triple stars, the classical double two-body model of zeroth-order cannot describe the motions of the components under the current observational accuracy. In this paper, Marchal's first-order analytical solution is implemented and a more efficient simplified version is applied to real triple stars. The results show that, for most triple stars, the proposed first-order model is preferable to the zeroth-order model both in fitting observational data and in predicting component positions.

We report on observations of recurrent jets by instruments onboard the Interface Region Imaging Spectrograph (IRIS), Solar Dynamics Observatory (SDO) and Hinode spacecrafts. Over a 4-hour period on July 21st 2013, recurrent coronal jets were observed to emanate from NOAA Active Region 11793. FUV spectra probing plasma at transition region temperatures show evidence of oppositely directed flows with components reaching Doppler velocities of +/- 100 km/s. Raster Doppler maps using a Si IV transition region line show all four jets to have helical motion of the same sense. Simultaneous observations of the region by SDO and Hinode show that the jets emanate from a source region comprising a pore embedded in the interior of a supergranule. The parasitic pore has opposite polarity flux compared to the surrounding network field. This leads to a spine-fan magnetic topology in the coronal field that is amenable to jet formation. Time-dependent data-driven simulations are used to investigate the underlying drivers for the jets. These numerical experiments show that the emergence of current-carrying magnetic field in the vicinity of the pore supplies the magnetic twist needed for recurrent helical jet formation.

Embedded clusters are ideal laboratories to understand the early phase of the dynamical evolution of clusters as well as the massive star formation. An interesting observational phenomenon is that some of the embedded clusters show mass segregation, i.e., the most massive stars are preferentially found near the cluster center. In this paper, we develop a new approach to describe mass segregation. Using this approach and the Two Micron All Sky Survey Point Source Catalog (2MASS PSC), we analyze eighteen embedded clusters in the Galaxy. We find that eleven of them are mass-segregated and that the others are non-mass-segregated. No inversely mass-segregated cluster is found.

This work is devoted to improving empirical mass-luminosity relations and mass-metallicity-luminosity relation for low mass stars. For these stars, observational data in the mass-luminosity plane or the mass-metallicity-luminosity space subject to non-negligible errors in all coordinates with different dimensions. Thus a reasonable weight assigning scheme is needed for obtaining more reliable results. Such a scheme is developed, with which each data point can have its own due contribution. Previous studies have shown that there exists a plateau feature in the mass-luminosity relation. Taking into account the constraints from the observational luminosity function, we find by fitting the observational data using our weight assigning scheme that the plateau spans from 0.28 to 0.50 solar mass. Three-piecewise continuous improved mass-luminosity relations in K, J, H and V bands, respectively, are obtained. The visual mass-metallicity-luminosity relation is also improved based on our K band mass-luminosity relation and the available observational metallicity data.

We simulate the bulge formation in very late-type dwarf galaxies from circumnuclear super star clusters (SSCs) moving in a responding cuspy dark matter halo (DMH). The simulations show that (1) the response of DMH to sinking of SSCs is detectable only in the region interior to about 200 pc. The mean logarithmic slope of the responding DM density profile over that area displays two different phases: the very early descent followed by ascent till approaching to 1.2 at the age of 2 Gyrs. (2) the detectable feedbacks of the DMH response on the bulge formation turned out to be very small, in the sense that the formed bulges and their paired nuclear cusps in the fixed and the responding DMH are basically the same, both are consistent with $HST$ observations. (3) the yielded mass correlation of bulges to their nuclear (stellar) cusps and the time evolution of cusps' mass are accordance with recent findings on relevant relations. In combination with the consistent effective radii of nuclear cusps with observed quantities of nuclear clusters, we believe that the bulge formation scenario that we proposed could be a very promising mechanism to form nuclear clusters.

In this paper a prescription for generating an equilibrium spherical system depicted with cuspy density profiles is extended to a core density profile, the Burkert profile, which is observationally more suitable to dwarfs. By using a time-saving Monte Carlo method instead of N-body simulations, we show that the Burkert halo initialized with the distribution function that depends only on energy is in a practically stable equilibrium. The one generated using the local Maxwellian approximation is unstable, where the flat core density structure tends to steepen. This is fundamentally different from the previous study on a halo with cuspy density profile. The deviation of the unstable "Burkert" from the initial Burkert profile is found to be closely related to taking a Gaussian as the velocity distribution at any given point. The significance of not using the local Maxwellian approximation is further demonstrated by exploring the dynamical evolution of compact super star clusters in the Burkert halo. In particular, the local Maxwellian approximation will result in underestimating the dynamical friction and overestimating sinking time scale. Accordingly, it will lead to lower probability of forming massive bulges and young nuclear star clusters in bulgeless galaxies.

The dynamical evolution of super star clusters has been investigated in dark matter halos depicted with a cuspy- or soft-core density profile. The simulations show that (1) exponential bulges with central cusps form in both cases;(2) distinctive bulge formation rates are derived for cuspy- and soft-core profile; (3) masses of bulges and nuclear clusters are more heavier in case of cuspy-core profile; (4) massive, nuclear star clusters possibly form with no discernible bulges at different ages in cuspy- or soft-core cases.

The dynamical evolution of super star clusters (SSCs) moving in the background of a dark matter halo has been investigated as a possible event responsible for the formation of bulges in late-type spirals. The underlying physical processes include sinking of SSCs due to the dynamical friction and stripping of SSCs on their way to the center. Based on the assumption of a universal density profile for the dark matter halo, and an isothermal model for the SSCs, our simulations have yielded bulges that are similar in many aspects to the observational ones. In particular, the derived surface density profiles can be well fitted by an exponential structure with nuclear cusps, which is consistent with HST observations.The preliminary simulations with the Burkert density profile yield very instructive predictions on the bulge formation, which is surely worth further investigating.

The dynamical evolution of super star clusters (SSCs) moving in the background of a dark matter halo has been investigated as a possible event responsible for the formation of bulges in late-type spirals. The underlying physical processes include sinking of SSCs due to the dynamical friction and stripping of SSCs on their way to the center. Our model calculations show that only sinking of circumnuclear SSCs contributes to the formation of galactic bulges at the early stage. Based on the assumption of a universal density profile for the dark matter halo, and an isothermal model for the SSCs, our simulations have yielded bulges that are similar in many aspects to the observational ones. In particular, the derived surface density profiles can be well fitted by anexponential structure with nuclear cusps, which is consistent with HST observations.

The dynamical evolution of super star clusters (SSCs) moving in the background of dark matter halo has been investigated as a possible event causing the observed correlation between the mass of galactic bulge, $M_{bulge}$, and the mass of its central black hole, $M_{BH}$. The involved physical processes are the sinking of SSCs due to the dynamical friction, and the stripping of SSCs on their way to the center. Model calculations show that only sinking of circumnuclear SSCs contribute to both the growth of the central object and the formation of the galactic bulge at the early stage. On the assumption of a universal density profile for the dark matter halo, and an isothermal model for the SSCs, our simulations have yielded the mass ratio of the central objects to the bulges formed this way to be about a few times $10^{-4}$, less than the observed median value for early type galaxies. It is, however, consistent with the observed mass ratio for disk spirals, implying that the proposed scenario might be a possible event for the formation of bulges and central black holes of late type galaxies and for $M_{BH} - M_{bulge}$ correlation of disk galaxies.