Astrophysics

Current archives of the LAMOST telescope contain millions of pipeline-processed spectra that have probably never been seen by human eyes. Most of the rare objects with interesting physical properties, however, can only be identified by visual analysis of their characteristic spectral features. A proper combination of interactive visualisation with modern machine learning techniques opens new ways to discover such objects. We apply active learning classification supported by deep convolutional networks to automatically identify complex emission-line shapes in multi-million spectra archives. We used the pool-based uncertainty sampling active learning driven by a custom-designed deep convolutional neural network with 12 layers inspired by VGGNet, AlexNet, and ZFNet, but adapted for one-dimensional feature vectors. The unlabelled pool set is represented by 4.1 million spectra from the LAMOST DR2 survey. The initial training of the network was performed on a labelled set of about 13000 spectra obtained in the region around H α by the 2m Perek telescope of the Ondřejov observatory, which mostly contains spectra of Be and related early-type stars. The differences between the Ondřejov intermediate-resolution and the LAMOST low-resolution spectrographs were compensated for by Gaussian blurring. After several iterations, the network was able to successfully identify emission-line stars with an error smaller than 6.5%. Using the technology of the Virtual Observatory to visualise the results, we discovered 1013 spectra of 948 new candidates of emission-line objects in addition to 664 spectra of 549 objects that are listed in SIMBAD and 2644 spectra of 2291 objects identified in an earlier paper of a Chinese group led by Wen Hou. The most interesting objects with unusual spectral properties are discussed in detail.

This project outlines the complete development of a variable star classification algorithm methodology. With the advent of Big-Data in astronomy, professional astronomers are left with the problem of how to manage large amounts of data, and how this deluge of information can be studied in order to improve our understanding of the universe. While our focus will be on the development of machine learning methodologies for the identification of variable star type based on light curve data and associated information, one of the goals of this work is the acknowledgment that the development of a true machine learning methodology must include not only study of what goes into the service (features, optimization methods) but a study on how we understand what comes out of the service (performance analysis). The complete development of a beginning-to-end system development strategy is presented as the following individual developments (simulation, training, feature extraction, detection, classification, and performance analysis). We propose that a complete machine learning strategy for use in the upcoming era of big data from the next generation of big telescopes, such as LSST, must consider this type of design integration.

A close partnership between people and partially autonomous machines has enabled decades of space exploration. But to further expand our horizons, our systems must become more capable. Increasing the nature and degree of autonomy - allowing our systems to make and act on their own decisions as directed by mission teams - enables new science capabilities and enhances science return. The 2011 Planetary Science Decadal Survey (PSDS) and on-going pre-Decadal mission studies have identified increased autonomy as a core technology required for future missions. However, even as scientific discovery has necessitated the development of autonomous systems and past flight demonstrations have been successful, institutional barriers have limited its maturation and infusion on existing planetary missions. Consequently, the authors and endorsers of this paper recommend that new programmatic pathways be developed to infuse autonomy, infrastructure for support autonomous systems be invested in, new practices be adopted, and the cost-saving value of autonomy for operations be studied.

We present the first version of the ALeRCE (Automatic Learning for the Rapid Classification of Events) broker light curve classifier. ALeRCE is currently processing the Zwicky Transient Facility (ZTF) alert stream, in preparation for the Vera C. Rubin Observatory. The ALeRCE light curve classifier uses variability features computed from the ZTF alert stream, and colors obtained from AllWISE and ZTF photometry. We apply a Balanced Random Forest algorithm with a two-level scheme, where the top level classifies each source as periodic, stochastic, or transient, and the bottom level further resolves each of these hierarchical classes, amongst 15 total classes. This classifier corresponds to the first attempt to classify multiple classes of stochastic variables (including core- and host-dominated active galactic nuclei, blazars, young stellar objects, and cataclysmic variables) in addition to different classes of periodic and transient sources, using real data. We created a labeled set using various public catalogs (such as the Catalina Surveys and {\em Gaia} DR2 variable stars catalogs, and the Million Quasars catalog), and we classify all objects with ≥6 g -band or ≥6 r -band detections in ZTF (868,371 sources as of 2020/06/09), providing updated classifications for sources with new alerts every day. For the top level we obtain macro-averaged precision and recall scores of 0.96 and 0.99, respectively, and for the bottom level we obtain macro-averaged precision and recall scores of 0.57 and 0.76, respectively. Updated classifications from the light curve classifier can be found at the \href{http://alerce.online}{ALeRCE Explorer website}.

We present a real-time stamp classifier of astronomical events for the ALeRCE (Automatic Learning for the Rapid Classification of Events) broker. The classifier is based on a convolutional neural network, trained on alerts ingested from the Zwicky Transient Facility (ZTF). Using only the \textit{science, reference} and \textit{difference} images of the first detection as inputs, along with the metadata of the alert as features, the classifier is able to correctly classify alerts from active galactic nuclei, supernovae (SNe), variable stars, asteroids and bogus classes, with high accuracy ( ∼ 94\%) in a balanced test set. In order to find and analyze SN candidates selected by our classifier from the ZTF alert stream, we designed and deployed a visualization tool called SN Hunter, where relevant information about each possible SN is displayed for the experts to choose among candidates to report to the Transient Name Server database. From June 26th 2019 to February 28th 2021, we have reported 6846 SN candidates to date (11.8 candidates per day on average), of which 971 have been confirmed spectroscopically. Our ability to report objects using only a single detection means that 70\% of the reported SNe occurred within one day after the first detection. ALeRCE has only reported candidates not otherwise detected or selected by other groups, therefore adding new early transients to the bulk of objects available for early follow-up. Our work represents an important milestone toward rapid alert classifications with the next generation of large etendue telescopes, such as the Vera C. Rubin Observatory.

The GWAC-N is an observation network composed of multi-aperture and multi-field of view robotic optical telescopes. The main instruments are the GWAC-A. Besides, several robotic optical telescopes with narrower field of views provide fast follow-up multi-band capabilities to the GWAC-N. The primary scientific goal of the GWAC-N is to search for the optical counterparts of GRB that will be detected by the SVOM. The GWAC-N performs many other observing tasks including the follow-ups of ToO and both the detection and the monitoring of variable/periodic objects as well as optical transients. To handle all of those scientific cases, we designed 10 observation modes and 175 observation strategies, especially, a joint observation strategy with multiple telescopes of the GWAC-N for the follow-up of GW events. To perform these observations, we thus develop an AOM system in charge of the object management, the dynamic scheduling of the observation plan and its automatic broadcasting to the network management and finally the image management. The AOM combines the individual telescopes into a network and smoothly organizes all the associated operations. The system completely meets the requirements of the GWAC-N on all its science objectives. With its good portability, the AOM is scientifically and technically qualified for other general purposed telescope networks. As the GWAC-N extends and evolves, the AOM will greatly enhance the discovery potential for the GWAC-N. In the first paper of a series of publications, we present the scientific goals of the GWAC-N as well as the hardware, the software and the strategy setup to achieve the scientific objectives. The structure, the technical design, the implementation and performances of the AOM will be also described in details. In the end, we summarize the current status of the GWAC-N and prospect for the development plan in the near future.

We describe a new algorithm to solve the time dependent, frequency integrated radiation transport (RT) equation implicitly, which is coupled to an explicit solver for equations of magnetohydrodynamics (MHD) using {\sf Athena++}. The radiation filed is represented by specific intensities along discrete rays, which are evolved using a conservative finite volume approach for both cartesian and curvilinear coordinate systems. All the terms for spatial transport of photons and interactions between gas and radiation are calculated implicitly together. An efficient Jacobi-like iteration scheme is used to solve the implicit equations. This removes any time step constrain due to the speed of light in RT. We evolve the specific intensities in the lab frame to simplify the transport step. The lab-frame specific intensities are transformed to the co-moving frame via Lorentz transformation when the source term is calculated. Therefore, the scheme does not need any expansion in terms of v/c . The radiation energy and momentum source terms for the gas are calculated via direct quadrature in the angular space. The time step for the whole scheme is determined by the normal Courant -- Friedrichs -- Lewy condition in the MHD module. We provide a variety of test problems for this algorithm including both optically thick and thin regimes, and for both gas and radiation pressure dominated flows to demonstrate its accuracy and efficiency.

An alternative scheme is described to estimate the layer 1 LAXPC 20 background for faint sources where the source contribution to the 50-80 keV count rate is less than 0.25 counts/sec (15 milli-crabs or 6? 10 ??1 ergs/s/cm 2 ). We consider 12 blank sky observations and based on their 50-80 keV count rate in 100 second time-bins, generate four template spectra which are then used to estimate the background spectrum and lightcurve for a given faint source observation. The variance of the estimated background subtracted spectra for the 12 blank sky observations is taken as the energy dependent systematic uncertainty which will dominate over the statistical one for exposures longer than 5 ksecs. The estimated 100 second time bin background lightcurve in the 4-20 keV band with a 3\% systematic error matches with the blank sky ones. The 4-20 keV spectrum can be constrained for a source with flux �? milli-crab. Fractional r.m.s variability of 10\% can be determined for a ?? milli-crab source lightcurve binned at 100 seconds. To illustrate the scheme, the lightcurves, and spectra of three different blank sky observations, three AGN sources (Mrk 0926, Mrk 110, NGC 4593), and LMC X-1 are shown.

The effective observation time of Imaging Air Cherenkov Telescopes (IACTs) plays an important role in the detection of gamma-ray sources, especially when the expected flux is low. This time is strongly limited by the atmospheric conditions. Significant extinction of Cherenkov light caused by the presence of clouds reduces the photon detection rate and also complicates or even makes impossible proper data analysis. However, for clouds with relatively high atmospheric transmission, high energy showers can still produce enough Cherenkov photons to allow their detection by IACTs. In this paper, we study the degradation of the detection capability of an array of small-sized telescopes for different cloud transmissions. We show the expected changes of the energy bias, energy and angular resolution and the effective collection area caused by absorption layers located at 2.5 and 4.5 km above the observation level. We demonstrate simple correction methods for reconstructed energy and effective collection area. As a result, the source flux that is observed during the presence of clouds is determined with a systematic error of < 20%. Finally, we show that the proposed correction method can be used for clouds at altitudes higher than 5 km a.s.l.. As a result, the analysis of data taken under certain cloudy conditions will not require additional time-consuming Monte Carlo simulations.

FAST is the largest single-dish aperture telescope with a cable-driven parallel robot introduced to achieve the highest sensitivity in the world. However, to realize the high-precision, mechanical equations of such robot are always complicated, so that it is difficult to achieve real-time control by the traditional iterative method. In this regard, this paper proposes an approximately analytical solution method, which uses the approximately linear relationship between the main parameters of FAST to bypass some iterations. With the coefficients of the relationship extracted, static or quasi-static mechanical equations can be analytically solved. In this paper's example, this method saves at least 90% of the calculating time and the calculated values are consistent with the experimental data. With such huge efficiency improvements, real-time and high-precision control of FAST will no longer be a difficult work. Besides, all the work in this paper is expected to be used in the FAST.