We use observational measurements of thermal and kinetic Sunyaev-Zel'dovich effects, as well as soft X-ray emission associated with galaxy groups to constrain the gas density and temperature in the circumgalactic medium (CGM) for dark matter halos with masses above $10^{12.5}M_{\odot}$. A number of generic models are used together with a Bayesian scheme to make model inferences. We find that gas with a single temperature component cannot provide a consistent model to match the observational data. A simple two-phase model assuming a hot component and an ionized warm component can accommodate all the three observations. The total amount of the gas in individual halos is inferred to be comparable to the universal baryon fraction corresponding to the halo mass. The inferred temperature of the hot component is comparable to the halo virial temperature. The fraction of the hot component increases from $(15-40)\%$ for $10^{12.5}{M}_\odot$ halos to $(40-60)\%$ for $10^{14.5}{M}_\odot$ halos, where the ranges reflect uncertainties in the assumed gas density profile. Our results suggest that a significant fraction of the halo gas is in a non-thermalized component with temperature much lower than the virial temperature.

Energetic electrons accelerated by solar flares often give rise to type III radio bursts at a broad waveband and even interplanetary type III bursts (IT3) if the wavelength extends to decameter-kilometer. In this Letter, we investigate the probability of the flares that produce IT3, based on the sample of 2272 flares above M-class observed from 1996 to 2016. It is found that only 49.6% of the flares are detected to be accompanied with IT3. The duration, peak flux, and fluence of the flares with and without IT3 both present power-law distributions in the frequency domain, but the corresponding spectral indices for the former (2.06$\pm$0.17, 2.04$\pm$0.18, and 1.55$\pm$0.09) are obviously smaller than that for the latter (2.82$\pm$0.22, 2.51$\pm$0.19, and 2.40$\pm$0.09), showing that the flares with IT3 have longer durations and higher peak fluxes. We further examine the relevance of coronal mass ejections (CMEs) to the two groups of flares. It is found that 58% (655 of 1127) of the flares with IT3 but only 19% (200 of 1078) of the flares without IT3 are associated with CMEs, and that the associated CMEs for the flares with IT3 are inclined to be wider and faster. This indicates that CMEs may also play a role in producing IT3, speculatively facilitating the escape of accelerated electrons from the low corona to the interplanetary space.

The Galactic magnetar SGR1935+2154 has been reported to produce the first known example of a bright millisecond duration radio burst (FRB 200428) similar to the cosmological population of fast radio bursts (FRBs), bolstering the association of FRBs to active magnetars. The detection of a coincident bright X-ray burst has revealed the first observed multi-wavelength counterpart of a FRB. However, the search for similar emission at optical wavelengths has been hampered by the high inferred extinction on the line of sight. Here, we present results from the first search for second-timescale emission from the source at near-infrared wavelengths using the Palomar Gattini-IR observing system in J-band, made possible by a recently implemented detector read-out mode that allowed for short exposure times of 0.84 s with 99.9% observing efficiency. With a total observing time of 12 hours (47728 images) on source, we place median $3\,\sigma$ limits on the second-timescale emission of $< 20$ mJy (13.1 AB mag). We present non-detection limits from epochs of four simultaneous X-ray bursts detected by the Insight-{\it HXMT} and {\it NuSTAR} telescopes during our observing campaign. The limits translate to an extinction corrected fluence limit of $< 125$ Jy ms for an estimated extinction of $A_J = 2.0$ mag. These limits provide the most stringent constraints to date on the fluence of flares at frequencies of $\sim 10^{14}$ Hz, and constrain the ratio of the near-infrared (NIR) fluence to that of coincident X-ray bursts to $R_{\rm NIR} < 2.5 \times 10^{-2}$. Our observations were sensitive enough to easily detect a near-infrared counterpart of FRB 200428 if the NIR emission falls on the same power law as that observed across its radio to X-ray spectrum. The non-detection of NIR emission around the coincident X-ray bursts constrains the fluence index of the brightest burst to be steeper than $0.35$.

A relationship between dust-reprocessed light from recent star formation and the amount of star-forming gas in a galaxy produces a correlation between WISE 12$\mu$m emission and CO line emission. Here we explore this correlation on kiloparsec scales with CO(1-0) maps from EDGE-CALIFA matched in resolution to WISE 12$\mu$m images. We find strong CO-12$\mu$m correlations within each galaxy (median Pearson $r = 0.85$) and we show that the scatter in the global CO-12$\mu$m correlation is largely driven by differences from galaxy to galaxy. The correlation is stronger than that between star formation rate and H$_2$ surface densities over the same set of pixels (median $r=0.71$). We explore multi-variable regression to predict $\Sigma(\mathrm{H_2})$ surface density using the WISE 12$\mu$m data combined with global and resolved galaxy properties, and provide the fit parameters for the best estimators. We find that $\Sigma(\mathrm{H_2})$ estimators that include $\Sigma(\mathrm{12\mu m})$ are able to predict $\Sigma(\mathrm{H_2})$ with $>10$% better accuracy than estimators that include resolved optical properties ($\Sigma(\mathrm{SFR}),\Sigma(\mathrm{M_*}),A_V$ and $12+\log\mathrm{O/H}$) instead of $\Sigma(\mathrm{12\mu m})$. The best single-property estimator is $\log \frac{\Sigma(\mathrm{H_2})}{\mathrm{M_\odot\>pc^{-2}}} = (0.49 \pm 0.01) + (0.71 \pm 0.01)\log \frac{\Sigma(\mathrm{12\mu m})}{\mathrm{L_\odot\>pc^{-2}}}$, with an average predictive accuracy of 0.19 dex per pixel, and intrinsic scatter of 0.17 dex. This correlation can be used to efficiently estimate H$_2$ surface densities down to at least $1 \> M_\odot \> \mathrm{pc^{-2}}$ on small spatial scales within nearby galaxies. This correlation may prove useful to probe even lower gas densities with the better mid-infrared sensitivities expected from the James Webb Space Telescope.

Historically, photometry has been largely used to identify stellar populations (MPs) in Globular Clusters (GCs) by using diagrams that are based on colours and magnitudes that are mostly sensitive to stars with different metallicities or different abundances of helium, carbon, nitrogen and oxygen. In particular, the pseudo two-colour diagram called chromosome map (ChM), allowed the identification and the characterization of MPs in about 70 GCs by using appropriate filters of the Hubble Space Telescope (HST) that are sensitive to the stellar content of He, C, N, O and Fe. We use here high-precision HST photometry from F275W, F280N, F343N, F373N, and F814W images of Omega Centauri to investigate its MPs. We introduce a new ChM whose abscissa and ordinate are mostly sensitive to stellar populations with different magnesium and nitrogen, respectively, in monometallic GCs. This ChM is effective in disentangling the MPs based on their Mg chemical abundances, allowing us to explore, for the first time, possible relations between the production of these elemental species for large samples of stars in GCs. By comparing the colours of the distinct stellar populations with the colours obtained from appropriate synthetic spectra we provide 'photometric-like' estimates of the chemical composition of each population. Our results show that, in addition to first generation (1G) stars, the metal-poor population of Omega Centauri hosts four groups of second-generation stars with different [N/Fe], namely, 2GA--D. 2GA stars share nearly the same [Mg/Fe] as the 1G, whereas 2GB, 2GC and 2GD stars are Mg depleted by ~0.15, ~0.25 and ~0.45 dex, respectively. We provide evidence that the metal-intermediate populations host stars with depleted [Mg/Fe].

Global non-hydrostatic atmospheric models are becoming increasingly important for studying the climates of planets and exoplanets. However, such models suffer from computational difficulties due to the large aspect ratio between the horizontal and vertical directions. To overcome this problem, we developed a global model using a vertically-implicit-correction (VIC) scheme in which the integration time step is no longer limited by the propagation of acoustic waves in the vertical. We proved that our model, based on the $\rm Athena^{++}$ framework and its extension for planetary atmospheres - SNAP (Simulating Non-hydrostatic Atmosphere on Planets), rigorously conserves mass and energy in finite volume simulations. We found that traditional numerical stabilizers such as hyper-viscosity and divergence damping are not needed when using the VIC scheme, which greatly simplifies the numerical implementation and improves stability. We present simulation results ranging from 1D linear waves to 3D global circulations with and without the VIC scheme. These tests demonstrate that our formulation correctly tracks local turbulent motions, produces Kelvin-Helmholtz instability, and generates a super-rotating jet on hot Jupiters. Employing this VIC scheme improves the computational efficiency of global simulations by more than two orders of magnitude compared to an explicit model and facilitates the capability of simulating a wide range of planetary atmospheres both regionally and globally.

Identifying galaxy groups from redshift surveys of galaxies plays an important role in connecting galaxies with the underlying dark matter distribution. Current and future high-$z$ spectroscopic surveys, usually incomplete in redshift sampling, present both opportunities and challenges to identifying groups in the high-$z$ Universe. We develop a group finder that is based on incomplete redshift samples combined with photometric data, using a machine learning method to assign halo masses to identified groups. Test using realistic mock catalogs shows that $\gtrsim 90\%$ of true groups with halo masses $\rm M_h \gtrsim 10^{12} M_{\odot}/h$ are successfully identified, and that the fraction of contaminants is smaller than $10\%$. The standard deviation in the halo mass estimation is smaller than 0.25 dex at all masses. We apply our group finder to zCOSMOS-bright and describe basic properties of the group catalog obtained.

Scaling laws of dust, HI gas and metal mass with stellar mass, specific star formation rate and metallicity are crucial to our understanding of the buildup of galaxies through their enrichment with metals and dust. In this work, we analyse how the dust and metal content varies with specific gas mass ($M_{\text{HI}}$/$M_{\star}$) across a diverse sample of 423 nearby galaxies. The observed trends are interpreted with a set of Dust and Element evolUtion modelS (DEUS) - incluidng stellar dust production, grain growth, and dust destruction - within a Bayesian framework to enable a rigorous search of the multi-dimensional parameter space. We find that these scaling laws for galaxies with $-1.0\lesssim \log M_{\text{HI}}$/$M_{\star}\lesssim0$ can be reproduced using closed-box models with high fractions (37-89$\%$) of supernova dust surviving a reverse shock, relatively low grain growth efficiencies ($\epsilon$=30-40), and long dus lifetimes (1-2\,Gyr). The models have present-day dust masses with similar contributions from stellar sources (50-80\,$\%$) and grain growth (20-50\,$\%$). Over the entire lifetime of these galaxies, the contribution from stardust ($>$90\,$\%$) outweighs the fraction of dust grown in the interstellar medium ($<$10$\%$). Our results provide an alternative for the chemical evolution models that require extremely low supernova dust production efficiencies and short grain growth timescales to reproduce local scaling laws, and could help solving the conundrum on whether or not grains can grow efficiently in the interstellar medium.

Fast radio bursts (FRBs) are mysterious millisecond-duration radio transients observed from extragalactic distances. Magnetars have been long speculated as the possible engine to power repeating bursts from FRB sources, but no convincing evidence has been collected so far. Recently, a Galactic magnetar dubbed Soft Gamma-ray Repeater (SGR) J1935+2154 entered an active phase by emitting intense soft gamma-ray bursts. One fast radio burst with two peaks (hereafter FRB 200428) and a luminosity comparable to those of extragalactic FRBs was detected from the source in association with a soft $\gamma$-ray / hard X-ray flare. Here we report an eight-hour targeted observational campaign with the Five-hundred-meter Aperture Spherical radio Telescope (FAST) in four observing sessions, assisted by multi-wavelength (optical and hard X-rays) observations. During the third session, twenty-nine SGR bursts were detected by the {\it Fermi} Gamma-ray Burst Monitor (GBM). Throughout the observing period, especially during the scrutinized epochs when SGR bursts reached the Earth, no single dispersed pulsed emission was detected by FAST. This places a fluence upper limit eight orders of magnitude deeper than the fluence of FRB 200428. Our results suggest that FRB -- SGR burst associations are rather rare. FRBs may be highly relativistic and geometrically beamed or may carry narrow spectra with the characteristic frequencies of the most bursts far outside the FAST band. More intriguingly, it is also possible that the physical condition to achieve coherent radiation in SGR bursts is contrived. Only under extreme conditions could an FRB be made in association with an SGR burst.

Fast radio bursts (FRBs) are short pulses observed in radio band from cosmological distances, some of which emit repeating bursts. The physical origins of these mysterious events have been subject to wide speculations and heated debates. One class of models invoke soft gamma-ray repeaters (SGRs), or magnetars, as the sources of FRBs. Magnetars are rotating neutron stars with extremely strong magnetic field and can sporadically emit bursts from X-ray ($\sim$keV) to soft gamma-ray ($\sim$sub-MeV) with duration from $10^{-2}$~s to $10^2$~s. However, even though some bright radio bursts have been observed from some magnetars, no FRB-like events had been detected to be associated with any magnetar burst, including one giant flare, and no radio burst has been associated with any X-ray event from any magnetar. Therefore, there is still no observational evidence for magnetar-FRB association up to today. Recently, a pair of FRB-like bursts (FRB~200428 hereafter) separated by 30 milliseconds (ms) were detected from the general direction of the Galactic magnetar SGR~J1935+2154. Here we report the detection of a non-thermal X-ray burst in the 1-250 keV energy band with the Insight-HXMT satellite, which we identify as emitted from SGR~J1935+2154. The burst showed two hard peaks with a separation of $\sim30$ ms, consistent with the separation between the two bursts in FRB~200428. The delay time between the double radio and X-ray peaks is $\sim8.57$ s, fully consistent with the dispersion delay of FRB~200428. We thus identify the non-thermal X-ray burst is associated with FRB~200428 whose high energy counterpart is the two hard peaks in X-ray. Our results suggest that the non-thermal X-ray burst and FRB~200428 share the same physical origin in an explosive event from SGR~J1935+2154.

We study the equivalence principle, regarded as the cornerstone of general relativity, by analyzing the deformation observable of black hole shadows. Such deformation can arise from new physics and may be expressed as a phenomenological violation of the equivalence principle. Specifically, we assume that there is an additional background vector field that couples to the photons around the supermassive black hole. This type of coupling yields impact on the way the system depends on initial conditions, and affects the black hole shadow at different wavelengths by a different amount, and therefore observations of the shadow in different wavelengths could constrain such couplings. This can be tested by future multi-band observations. Adopting a specific form of the vector field, we obtain constraints on model parameters from Event Horizon Telescope observations and measurements of gas/stellar orbits.

In the era of big data astronomy, next generation telescopes and large sky surveys produce data sets at the TB or even PB level. Due to their large data volumes, these astronomical data sets are extremely difficult to transfer and analyze using personal computers or small clusters. In order to offer better access to data, data centers now generally provide online science platforms that enable analysis close to the data. The Chinese Virtual Observatory (China-VO) is one of the member projects in the International Virtual Observatory Alliance and it is dedicated to providing a research and education environment where globally distributed astronomy archives are simple to find, access, and interoperate. In this study, we summarize highlights of the work conducted at the China-VO, as well the experiences and lessons learned during the full life-cycle management of astronomical data. Finally, We discuss the challenges and future trends for astronomical science platforms.

We analyse two-dimensional maps and radial profiles of the equivalent width of H$\alpha$ emission (EW(H$\alpha$)), H$\delta$ absorption (EW(H$\delta_A$)), and the 4000\AA\ break (D$_n$(4000)) of low-redshift galaxies using integral field spectroscopy from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey. Out of $\approx2400$ nearly face-on galaxies with a redshift $z<0.05$ from the MaNGA Products Launch 7, we identify 124 "turnover" galaxies that each have a central upturn in EW(H$\alpha$), EW(H$\delta_A$) and/or a central drop in D$_n$(4000), indicative of ongoing/recent star formation (SF). The turnover features are found mostly in galaxies with a stellar mass above $\sim$10$^{10}$ M$_{\odot}$ and NUV-$r$ colour less than $\approx5$. The majority of the turnover galaxies are barred, with a bar fraction of 87$\pm$3\%. Furthermore, for barred galaxies the radius of the central turnover region ($R_{\rm{t}}$) is found to tightly correlate with the bar length ($R_{\rm{bar}}$): $R_{\rm{t}}$ = (0.35$\pm$0.02)$R_{\rm{bar}}$ + (0.61$\pm$0.13). Comparing the observed and the inward extrapolated star formation rate surface density ($\Sigma_{\rm SFR}$), we estimate that the central SF rate have been enhanced by an order of magnitude. Conversely, only half of the barred galaxies in our sample have a central turnover feature, implying that the presence of a bar is not sufficient to lead to a central SF enhancement. We further examined the SF enhancement in paired galaxies, as well as the density of the local environment, finding no relation. This implies that environment is not a driving factor for central turnover and SF enhancement in our sample galaxies.

Wolf-Rayet (WR) galaxies are a rare population of galaxies that host living high-mass stars during their WR phase (i.e. WR stars) and are thus expected to provide interesting constraints on the stellar Initial Mass Function, massive star formation, stellar evolution models, etc. Spatially resolved spectroscopy should in principle provide a more efficient way of identifying WR galaxies than single-fiber surveys of galactic centers such as SDSS-I & II, as WR stars should be more preferentially found in discs. Using IFU data from the ongoing SDSS-IV MaNGA survey, we have performed a thorough search for WR galaxies. We first identify H II regions in each datacube and carry out full spectral fitting to the stacked spectra. We then visually inspect the residual spectrum of each H II region and identify WR regions that present a significant "blue bump" at 4600-4750 A. The resulting WR catalog includes 267 WR regions of ~500pc (radius) sizes, distributed in 90 galaxies from the current sample of MaNGA (MaNGA Product Launch 7). We find WR regions are exclusively found in galaxies that show bluest colors and highest star formation rates for their mass. Most WR galaxies have late-type morphologies and show relatively large asymmetry in their images, implying that WR regions are more preferentially found in interacting/merging galaxies. We estimate the stellar mass function of WR galaxies and the mass-dependent detection rate. The detection rate of WR galaxies is typically ~2%, with weak dependence on stellar mass. This detection rate is about 40 times higher than previous studies with SDSS single fiber data, and by a factor of 2 lower than the CALIFA-based WR catalog. We make comparisons with SDSS and CALIFA studies, and conclude that different detection rates can be explained mainly by three factors: spatial coverage, spectral signal-to-noise ratio, and redshift ranges of the parent sample.

We study the mass distribution and kinematics of the giant elliptical galaxy M87 (NGC 4486) using discrete chemo-dynamical, axisymmetric Jeans equation modelling. Our catalogue comprises 894 globular clusters (GCs) extending to a projected radius of $\sim 430$ kpc with line-of-sight velocities and colours, and Multi Unit Spectroscopic Explorer (MUSE) integral field unit data within the central $2.4$ kpc of the main galaxy. The gravitational potential for our models is a combination of a luminous matter potential with a varying mass-to-light ratio for the main galaxy, a supermassive black hole and a dark matter (DM) potential with a cusped or cored DM halo. The best-fitting models with either a cusped or a cored DM halo show no significant differences and both are acceptable. We obtain a total mass of $(2.16 \pm 0.38) \times 10^{13} M_{\odot}$ within $\sim$ 400 kpc. By including the stellar mass-to-light ratio gradient, the DM fraction increases from $\sim$ 26 percent (with no gradient) to $\sim$ 73 percent within $1\,R_e^{\rm maj}$ (major axis of half-light isophote, 14.2 kpc), and from $\sim$ 84 percent to $\sim$ 94 percent within $5\,R_e^{\rm maj}$ (71.2 kpc). Red GCs have moderate rotation with $V_{\rm max}/\sigma \sim$ 0.4, and blue GCs have weak rotation with $V_{\rm max}/\sigma \sim$ 0.1. Red GCs have tangential velocity dispersion anisotropy, while blue GCs are consistent with being nearly isotropic. Our results suggest that red GCs are more likely to be born in-situ, while blue GCs are more likely to be accreted.

We present multi-wavelength measurements of the thermal, chemical, and cloud contrasts associated with the visibly dark formations (also known as 5-$\mu$m hot spots) and intervening bright plumes on the boundary between Jupiter's Equatorial Zone (EZ) and North Equatorial Belt (NEB). Observations made by the TEXES 5-20 $\mu$m spectrometer at the Gemini North Telescope in March 2017 reveal the upper-tropospheric properties of 12 hot spots, which are directly compared to measurements by Juno using the Microwave Radiometer (MWR), JIRAM at 5 $\mu$m, and JunoCam visible images. MWR and thermal-infrared spectroscopic results are consistent near 0.7 bar. Mid-infrared-derived aerosol opacity is consistent with that inferred from visible-albedo and 5-$\mu$m opacity maps. Aerosol contrasts, the defining characteristics of the cloudy plumes and aerosol-depleted hot spots, are not a good proxy for microwave brightness. The hot spots are neither uniformly warmer nor ammonia-depleted compared to their surroundings at $p<1$ bar. At 0.7 bar, the microwave brightness at the edges of hot spots is comparable to other features within the NEB. Conversely, hot spots are brighter at 1.5 bar, signifying either warm temperatures and/or depleted NH$_3$ at depth. Temperatures and ammonia are spatially variable within the hot spots, so the precise location of the observations matters to their interpretation. Reflective plumes sometimes have enhanced NH$_3$, cold temperatures, and elevated aerosol opacity, but each plume appears different. Neither plumes nor hot spots had microwave signatures in channels sensing $p>10$ bars, suggesting that the hot-spot/plume wave is a relatively shallow feature.

This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

We develop a method to estimate the dust attenuation curve of galaxies from full spectral fitting of their optical spectra. Motivated from previous studies, we separate the small-scale features from the large-scale spectral shape, by performing a moving average method to both the observed spectrum and the simple stellar population model spectra. The intrinsic dust-free model spectrum is then derived by fitting the observed ratio of the small-scale to large-scale (S/L) components with the S/L ratios of the SSP models. The selective dust attenuation curve is then determined by comparing the observed spectrum with the dust-free model spectrum. One important advantage of this method is that the estimated dust attenuation curve is independent of the shape of theoretical dust attenuation curves. We have done a series of tests on a set of mock spectra covering wide ranges of stellar age and metallicity. We show that our method is able to recover the input dust attenuation curve accurately, although the accuracy depends slightly on signal-to-noise ratio of the spectra. We have applied our method to a number of edge-on galaxies with obvious dust lanes from the ongoing MaNGA survey, deriving their dust attenuation curves and $E(B-V)$ maps, as well as dust-free images in $g$, $r$, and $i$ bands. These galaxies show obvious dust lane features in their original images, which largely disappear after we have corrected the effect of dust attenuation. The vertical brightness profiles of these galaxies become axis-symmetric and can well be fitted by a simple model proposed for the disk vertical structure. Comparing the estimated dust attenuation curve with the three commonly-adopted model curves, we find that the Calzetti curve provides the best description of the estimated curves for the inner region of galaxies, while the Milky Way and SMC curves work better for the outer region.

The CsI detectors of the High Energy X-ray Telescope of the Hard X-ray Modulation Telescope (HXMT/CsI) can be used for gamma-ray all sky monitoring and searching for the electromagnetic counterpart of gravitational wave source. The instrumental responses are mainly obtained by Monte Carlo simulation with the Geant4 tool and the mass model of both the satellite and all the payloads, which is updated and tested with the Crab pulse emission in various incident directions. Both the Energy-Channel relationship and the energy resolution are calibrated in two working modes (Normal-Gain mode & Low-Gain Mode) with the different detection energy ranges. The simulative spectral analyses show that HXMT/CsI can constrain the spectral parameters much better in the high energy band than that in the low energy band. The joint spectral analyses are performed to ten bright GRBs observed simultaneously with HXMT/CsI and other instruments (Fermi/GBM, Swift/BAT, Konus-Wind), and the results show that the GRB flux given by HXMT/CsI is systematically higher by $7.0\pm8.8\%$ than those given by the other instruments. The HXMT/CsI-Fermi/GBM joint fittings also show that the high energy spectral parameter can be constrained much better as the HXMT/CsI data are used in the joint fittings.

We use the In data collected during the 2019 outburst from X-ray pulsar 4U 1901+03 to complement the orbital parameters reported by Fermi/GBM. Using the Insight-HXMT, we examine the correlation between the derivative of the intrinsic spin frequency and bolometric flux based on accretion torque models. It was found that the pulse profiles significantly evolve during the outburst. The existence of two types of the profile's pattern discovered in the Insight-HXMT data indicates that this source experienced transition between a super-critical and a sub-critical accretion regime during its 2019 outburst. Based on the evolution of the pulse profiles and the torque model, we derive the distance to 4U 1901+03 as 12.4+-0.2 kpc.

We study the evolution of the temporal properties of MAXI 1820+070 during the 2018 outburst in its hard state from MJD 58190 to 58289 with Insight-HXMT in a broad energy band 1-150 keV. We find different behaviors of the hardness ratio, the fractional rms and time lag before and after MJD 58257, suggesting a transition occurred at around this point. The observed time lags between the soft photons in the 1-5 keV band and the hard photons in higher energy bands, up to 150 keV, are frequency-dependent: the time lags in the low-frequency range, 2-10 mHz, are both soft and hard lags with a timescale of dozens of seconds but without a clear trend along the outburst; the time lags in the high-frequency range, 1-10 Hz, are only hard lags with a timescale of tens of milliseconds; first increase until around MJD 58257 and decrease after this date. The high-frequency time lags are significantly correlated to the photon index derived from the fit to the quasi-simultaneous NICER spectrum in the 1-10 keV band. This result is qualitatively consistent with a model in which the high-frequency time lags are produced by Comptonization in a jet.

The 21~cm intensity mapping experiments promise to measure the baryon acoustic oscillations (BAO) through the mapping of H{\small I} gas at the post-reionization epoch. Its success depends on the assumption that H{\small I} gas traces the matter density distribution {\it linearly} on large scales. However, both nonlinear halo clustering and nonlinear effects modulating H{\small I} gas in halos may spoil this. We employ three approaches to generate the mock H{\small I} density from a large-scale N-body simulation at low redshifts, and demonstrate that the assumption of H{\small I} linearity is valid at the scale corresponding to the first BAO peak, but breaks down at $k \gtrsim 0.1\,h\, {\rm Mpc}^{-1}$. The nonlinear effects of halo clustering and H{\small I} content modulation counteract each other at small scales, and their competition results in a model-dependent "sweet-spot" redshift near $z$=1 where the H{\small I} bias is scale-independent down to small scales. We also find that the linear H{\small I} bias scales approximately linearly with redshift for $z\le 3$.

We investigate the co-evolution of dark matter (DM) density perturbation and metric perturbation in the WIMP paradigm. Instead of adopting the thermal equilibrium assumption as conventional, we propose a simple phase of DM production for the WIMP paradigm and extend our analysis to this phase. Being free from the envelop of thermal equilibrium, an amplified perturbative resonance between DM density perturbation and scalar modes of metric perturbation takes place during the DM production phase, and consequently results in a suppression of the tensor-to-scalar ratio of metric perturbation. By specifying the cosmic background with a typical realization of cosmic reheating, we establish a new relation between DM particle mass $m_\chi$ and the tensor-to-scalar ratio $r$ in the WIMP paradigm, which also contains two reheating parameters, the reheating temperature $\Gamma_0$ and the dissipative constant $T_{R_f}$. Notably, for a sizeable parameter region of WIMP candidate and cosmic reheating, this relation predicts a smaller value of $r$ in comparing with the conventional expectation obtained by ignoring the process of DM production. Once the suppression of $r$ is measured in future observations of primordial gravitational wave in CMB, this relation can be used to constrain $m_\chi$, $T_{R_f}$ and $\Gamma_0$ in principle.

Data of $^{12}$CO/$^{13}$CO/C$^{18}$O $J=1\to0$ emission toward the Galactic plane region of $l=35^\circ$ to $45^\circ$ and $b= -5^\circ$ to $+5^\circ$ are available with the Milky Way Imaging Scroll Painting (MWISP) project. Using the data, we found a giant molecular filament (GMF) around $l\approx38\sim42^\circ$, $b\approx-3.5\sim0^\circ$, $V_{LSR} \approx 27 \sim 40$ km~s$^{-1}$, named the GMF MWISP G041-01. At a distance of 1.7 kpc, the GMF is about 160 pc long. With a median excitation temperature about 7.5 K and a median column density about $10^{21}$ cm$^{-2}$, this GMF is very cold and very diffuse compared to known GMFs. Using the morphology in the data cube, the GMF is divided into four components among which three show filamentary structure. Masses of the components are $ 10^3 \sim 10^4 M_\odot$, with a total mass for the whole filament being about $7\times10^4 M_\odot$ from the LTE method. $^{13}$CO cores inside each component are searched. Virial parameters are about 2.5 for these cores and have a power-law index of -0.34 against the mass. The mass fraction of dense cores traced by $^{13}$CO to the diffuse clouds traced by $^{12}$CO are about 7% for all components of the GMF. We found signatures of possible large scale filament-filament collision in the GMF.

Glitches correspond to sudden jumps of rotation frequency ($\nu$) and its derivative ($\dot{\nu}$) of pulsars, the origin of which remains not well understood yet, partly because the jump processes of most glitches are not well time-resolved. There are three large glitches of the Crab pulsar, detected in 1989, 1996 and 2017, which were found to have delayed spin-up processes before the normal recovery processes. Here we report two additional glitches of the Crab pulsar occurred in 2004 and 2011 for which we discovered delayed spin up processes, and present refined parameters of the largest glitch occurred in 2017. The initial rising time of the glitch is determined as $<0.48$ hour. We also carried out a statistical study of these five glitches with observed spin-up processes. The two glitches occurred in 2004 and 2011 have delayed spin-up time scales ($\tau_{1}$) of $1.7\pm0.8$\,days and $1.6\pm0.4$\,days, respectively. We find that the $\Delta{\nu}$ vs. $|\Delta{\dot\nu}|$ relation of these five glitches is similar to those with no detected delayed spin-up process, indicating that they are similar to the others in nature except that they have larger amplitudes. For these five glitches, the amplitudes of the delayed spin-up process ($|\Delta{\nu}_{\rm d1}|$) and recovery process ($\Delta{\nu}_{\rm d2}$), their time scales ($\tau_{1}$, $\tau_{2}$), and permanent changes in spin frequency ($\Delta{\nu}_{\rm p}$) and total frequency step ($\Delta{\nu}_{\rm g}$) have positive correlations. From these correlations, we suggest that the delayed spin-up processes are common for all glitches, but are too short and thus difficult to be detected for most glitches.

Identifying the sources of the highest energy cosmic rays requires understanding how they are deflected by the stochastic, spatially intermittent intergalactic magnetic field. Here we report measurements of energetic charged-particle propagation through a laser-produced magnetized plasma with these properties. We characterize the diffusive transport of the particles experimentally. The results show that the transport is diffusive and that, for the regime of interest for the highest-energy cosmic rays, the diffusion coefficient is unaffected by the spatial intermittency of the magnetic field.

No robust detection of prompt electromagnetic counterparts to fast radio bursts (FRBs) has yet been obtained, in spite of several multi-wavelength searches carried out so far. Specifically, X/gamma-ray counterparts are predicted by some models. We planned on searching for prompt gamma-ray counterparts in the Insight-Hard X-ray Modulation Telescope (Insight-HXMT) data, taking advantage of the unique combination of large effective area in the keV-MeV energy range and of sub-ms time resolution. We selected 39 FRBs that were promptly visible from the High-Energy (HE) instrument aboard Insight-HXMT. After calculating the expected arrival times at the location of the spacecraft, we searched for a significant excess in both individual and cumulative time profiles over a wide range of time resolutions, from several seconds down to sub-ms scales. Using the dispersion measures in excess of the Galactic terms, we estimated the upper limits on the redshifts. No convincing signal was found and for each FRB we constrained the gamma-ray isotropic-equivalent luminosity and the released energy as a function of emission timescale. For the nearest FRB source, the periodic repeater FRB180916.J0158+65, we find $L_{\gamma,iso}<5.5\times 10^{47}$ erg/s over 1 s, whereas $L_{\gamma,iso}<10^{49}-10^{51}$ erg/s for the bulk of FRBs. The same values scale up by a factor of ~100 for a ms-long emission. Even on a timescale comparable with that of the radio pulse itself no keV-MeV emission is observed. A systematic association with either long or short GRBs is ruled out with high confidence, except for subluminous events, as is the case for core-collapse of massive stars (long) or binary neutron star mergers (short) viewed off axis. Only giant flares from extra-galactic magnetars at least ten times more energetic than Galactic siblings are ruled out for the nearest FRB.

To trace the unknown particle nature of dark matter (DM) from observations in CMB, we study the co-evolution of DM density perturbation and primordial metric perturbation in the early Universe. By analytically solving the improved Einstein-Boltzmann equation to include all perturbative contributions from the DM productions, we observe an amplified perturbative resonance that takes place between DM density perturbation and scalar modes of metric perturbation. As this perturbative resonance does not affect the tensor modes, it results in a suppression of the tensor-to-scalar ratio of metric perturbation $r$. Then by using a typical reheating model to specify cosmic background, we establish an explicit relation between DM particle mass $m_\chi$ and $r$. As $r$ will be further constrained in the coming round of primordial gravitational wave detections, $m_\chi$ would be constrained accordingly.

The Moon's optical characteristics in visible and long-wavelength infrared (LWIR) have long been observed with our eyes or with instruments. What the mid-infrared (MIR) Moon looks like is still a mystery. For the first time we present detailed appearance of the MIR Moon observed by a high-resolution geostationary satellite and reveal the essence behind its appearance. The appearance of the MIR Moon is opposite to its normal visible appearance. In addition the MIR Moon shows limb darkening. Both the absolute and the relative brightness distribution of the MIR lunar disk changes with the solar incidence angle. The signatures of the MIR Moon are controlled by both the reflection and emission of the lunar surface. We also show first-ever brightness temperature maps of the lunar disk without needing a mosaic, which better show the temperature variation across the lunar disk. They reveal that the relationship between brightness temperature and solar incidence angle i is cos1/bi, and the power parameter is smaller than the Lambertian temperature model of cos1/4i observed for lunar orbit-based measurements. The slower decrease of the brightness temperature when moving away from the sub-solar point than the Lambertian model is due to topographic effects. The brightness temperature is dominated by albedo and the solar incidence angle and influenced by the topography. Our results indicate that the Moon in the MIR exhibits many interesting phenomena which were previously unknown, and contains abundant information about lunar reflection and thermal emission for future study.

We measure the star formation histories (SFH) of a sample of low-mass galaxies with $M_\ast<10^9M_\odot$ from the SDSS-IV MaNGA survey. The large number of IFU spectra for each galaxy are either combined to reach a high signal to noise ratio or used to investigate spatial variations. We use Bayesian inferences based on full spectrum fitting. Our analysis based on Bayesian evidence ratio indicates a strong preference for a model that allows the presence of an old stellar population, and that an improper model for the SFH can significantly underestimate the old population in these galaxies. The addition of NIR photometry to the constraining data can further distinguish between different SFH model families and significantly tighten the constraints on the mass fraction in the old population. On average more than half of the stellar mass in present-day low-mass galaxies formed 8 Gyrs ago, while about 30\% within the past 4 Gyrs. Satellite galaxies on average have formed their stellar mass earlier than central galaxies, and stars in the outer regions of galaxy are younger than those in the central part. Our results suggest that most of the low-mass galaxies have an early episode of active star formation that produces a large fraction of their present stellar mass.

The Galactic plane scanning survey is one of the main scientific objectives of the Hard X-ray Modulation Telescope (known as Insight-HXMT). During the two-year operation of Insight-HXMT, more than 1000 scanning observations have been performed and the whole Galactic plane ($\rm 0^{\circ}<l<360^{\circ}$, $\rm -10^{\circ}<b<10^{\circ}$) has been covered completely. We summarize the Galactic plane scanning survey of Insight-HXMT for two years, including the characteristics of the scanning data, the data analysis process and the preliminary results of the Low-Energy telescope, the Medium-Energy telescope and the High-Energy telescope. With the light curve PSF fitting method, the fluxes of the known sources in the scanned area as well as the flux errors are obtained for each scanning observation. From the relationships of SNRs and fluxes, the $5\sigma$ sensitivities of three telescopes of Insight-HXMT are estimated as $\rm \sim7.6\times10^{-11}~erg cm^{-2}~s^{-1}$ ($\rm 3 mCrab,~1-6 keV$), $\rm \sim4.0\times10^{-10}~erg~cm^{-2}~s^{-1}$ ($\rm 20~mCrab,~7-40 keV$) and $\rm \sim2.6\times10^{-10}~erg cm^{-2}~s^{-1}$ ($\rm 18 mCrab,~25-100 keV$) for an individual scanning observation of $2-3$ hours, respectively. Up to September 2019, more than 800 X-ray sources with various types are monitored by the three telescopes and their long-term light curves with three energy bands are obtained to make further scientific analyses.

We present the calibration of the Insight-Hard X-ray Modulation Telescope (Insight-HXMT) X-ray satellite, which can be used to perform timing and spectral studies of bright X-ray sources. Insight-HXMT carries three main payloads onboard: the High Energy X-ray telescope (HE), the Medium Energy X-ray telescope (ME) and the Low Energy X-ray telescope (LE). In orbit, the radioactive sources, activated lines, the fluorescence lines and celestial sources are used to calibrate the energy scale and energy resolution of the payloads. The Crab nebular is adopted as the primary effective area calibrator and empirical functions are constructed to modify the simulated effective areas of the three payloads respectively. The systematic errors of HE, compared to the model of the Crab nebular, are less than 2% in 28--120 keV and 2%--10% above 120 keV. The systematic errors of ME are less than 1.5% in 10--35 keV. The systematic errors of LE are less than 1% in 1--7 keV except the Si K--edge (1.839 keV, up to 1.5%) and less than 2% in 7--10 keV.

We use the ELUCID N-Body simulation to study the relation of structural properties of dark matter halos to their assembly history and environment. We find that the complexity of halo assembly histories can be well described by a small number of principal components. Using decision trees built with the random ensemble method, we find that the correlation between halo concentration and assembly history is tight: more than 60% of the variance in halo concentration can be explained by assembly history alone. Halo shape and spin parameters are found to depend not only on assembly history, but also on environment and initial condition, and these dependencies are not orthogonal. The structures of more massive halos depend more strongly on higher-order components of the assembly history. The relations between halo intrinsic properties and environment are weak compared to their variances, with the anisotropy of local tidal fields having the strongest correlation with halo properties. Our method of dimension reduction and regression can help simplify the characterization of the halo population and clarify the degeneracy among halo properties.

The WorldWide Telescope(WWT) is a scientific visualization platform which can browse deep space images, star catalogs, and planetary remote sensing data from different observation facilities in a three-dimensional virtual scene. First launched and then open-sourced by Microsoft Research, the WWT is now managed by the American Astronomical Society (AAS). Hierarchical Progressive Survey (HiPS) is an astronomical data release scheme proposed by Centre de Donn\'ees astronomiques de Strasbourg (CDS) and has been accepted as a recommendation by International Virtual Observatory Alliance (IVOA). The HiPS solution has been adopted widely by many astronomical institutions for data release. Since WWT selected Hierarchical Triangular Mesh (HTM) as the standard for data visualization in the early stage of development, data released by HiPS cannot be visualized in WWT, which significantly limits the application of WWT. This paper introduces the implementation method for HiPS dataset visualization in WWT, and introduces HiPS data projection, mesh rendering, and data index implementation in WWT. Taking Chang'E-2 lunar probe data as an example, this paper introduces how to convert planetary remote sensing data into a HiPS dataset and integrate it into WWT. This paper also compares the efficiency and memory consumption of WWT loading its native data and HiPS data, and illustrates the application of HiPS in scientific data visualization and science education in WWT.

Recent studies have revealed that the onset age for the presence of multiple stellar populations (MPs) in star clusters seems to correspond to the disappearance of the extended main-sequence turnoff (eMSTO) in young clusters, a pattern associated with stellar rotations. A speculative suggestion is that MPs might be caused by the magnetic brake, a stellar evolutionary effect linked to the rotation. In this work, we use the young massive cluster NGC 419 as a testbed. We examined if its magnetically baked MS stars would exhibit MPs. Using the deep ultraviolet and visible images observed through the Hubble Space Telescope, combined with a specific color index that is sensitive to the nitrogen (N) abundance, we examined if its late G- and K-type MS stars are affected by N variation. Our analysis reports that the morphology of its GK-type MS is most likely an SSP, and only a negligible probability that indicates a N variation up to 0.4 dex is present. We conclude that there is no significant N variation among its GK-type MS stars. The absence of a significant chemical variation among the late-type MS stars indicates that MPs might not be a specific pattern of magnetically braked stars.

We report on our analysis of the 2019 outburst of the X-ray accreting pulsar 4U 1901+03 observed with Insight-HXMT and NICER. Both spectra and pulse profiles evolve significantly in the decaying phase of the outburst. Dozens of flares are observed throughout the outburst. They are more frequent and brighter at the outburst peak. We find that the flares, which have a duration from tens to hundreds of seconds, are generally brighter than the persistent emission by a factor of $\sim$ 1.5. The pulse profile shape during the flares can be significantly different than that of the persistent emission. In particular, a phase shift is clearly observed in many cases. We interpret these findings as direct evidence of changes of the pulsed beam pattern, due to transitions between the sub- and super-critical accretion regimes on a short time scale. We also observe that at comparable luminosities the flares' pulse profiles are rather similar to those of the persistent emission. This indicates that the accretion on the polar cap of the neutron star is mainly determined by the luminosity, i.e., the mass accretion rate.

We make the in-orbit calibration to the point-spread functions (PSFs) of the collimators of the Hard X-ray Modulation Telescope with the scanning observation of the Crab. We construct the empirical adjustments to the theoretically calculated geometrical PSFs. The adjustments contain two parts: a rotating matrix to adjust the directional deviation of the collimators and a paraboloidal function to correct the inhomogeneity of the real PSFs. The parameters of the adjusting matrices and paraboloidal functions are determined by fitting the scanning data with lower scanning speed and smaller intervals during the calibration observations. After the PSF calibration, the systematic errors in source localization in the Galactic plane scanning survey are 0.010 deg, 0.015 deg, 0.113 deg for the Low-Energy Telescope (LE), the Medium-Energy telescope (ME) and the High-Energy telescope (HE), respectively; meanwhile, the systematic errors in source flux estimation are 1.8%, 1.6%, 2.7% for LE, ME and HE, respectively.

A new black hole X-ray binary (BHXRB) MAXI J1535-571 was discovered by MAXI during its outburst in 2017. Using observations taken by the first Chinese X-ray satellite, the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT), we perform a joint spectral analysis (2-150 keV) in both energy and time domains. The energy spectra provide the essential input for probing the intrinsic Quasi-Periodic Oscillation (QPO) fractional rms spectra (FRS). Our results show that during the intermediate state, the energy spectra are in general consistent with those reported by Swift/XRT and NuSTAR. However, the QPO FRS become harder and the FRS residuals may suggest the presence of either an additional power-law component in the energy spectrum or a turn-over in the intrinsic QPO FRS at high energies.

Efforts to detect low-mass exoplanets using stellar radial velocities (RVs) are currently limited by magnetic photospheric activity. Suppression of convective blueshift is the dominant magnetic contribution to RV variability in low-activity Sun-like stars. Due to convective plasma motions, the magnitude of RV contributions from the suppression of convective blueshift is related to the depth of formation of photospheric spectral lines of a given species used to compute the RV time series. Meunier et al. (2017), used this relation to demonstrate a method for spectroscopic extraction of the suppression of convective blueshift in order to isolate RV contributions, including planetary RVs, that contribute equally to the timeseries for each spectral line. Here, we extract disk-integrated solar RVs from observations over a 2.5 year time span made with the solar telescope integrated with the HARPS-N spectrograph at the Telescopio Nazionale Galileo (La Palma, Canary Islands, Spain). We apply the methods outlined by Meunier et al. (2017). We are not, however, able to isolate physically meaningful contributions of the suppression of convective blueshift from this solar dataset, potentially because our dataset is from solar minimum when the suppression of convective blueshift may not sufficiently dominate activity contributions to RVs. This result indicates that, for low-activity Sun-like stars, one must include additional RV contributions from activity sources not considered in the Meunier et al. (2017) model at different timescales as well as instrumental variation in order to reach the sub-meter per second RV sensitivity necessary to detect low-mass planets in orbit around Sun-like stars.

We study the 2018 outburst of Aql X-1 via the monitor of all sky X-ray image (MAXI) data. We show that the outburst starting in February 2018 is a member of short-low class in the frame of outburst duration and the peak count rate although the outburst morphology is slightly different from the other fast-rise-exponential-decay (FRED) type outbursts with a milder rising stage. We study the partial accretion in the weak propeller stage of Aql X-1 via the MAXI data of the 2018 outburst. We report on the spectral analysis of 3 observations of Aquila X-1 obtained by Insight - hard X-ray modulation telescope (Insight-HXMT) during the late decay stage of the 2018 outburst. We discuss that the data taken by Insight-HXMT is just after the transition to the weak propeller stage. Our analysis shows the necessity of a comptonization component to take into account the existence of an electron cloud resulting photons partly up-scattered.

Standard stellar evolution model predicts a severe depletion of lithium (Li) abundance during the first dredge up process (FDU). Yet a small fraction of giant stars are still found to preserve a considerable amount of Li in their atmospheres after the FDU. Those giants are usually identified as Li-rich by a widely used criterion, A(Li) > 1.5 dex. A large number of works dedicated to searching for investigating this minority of the giant family, and the amount of Li-rich giants, has been largely expanded on, especially in the era of big data. In this paper, we present a catalog of Li-rich giants found from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) survey with Li abundances derived from a template-matching method developed for LAMOST low-resolution spectra.

In this paper we study implications of the possible excess of 21-cm line global signal at the epoch of cosmic dawn on the evolutions of a class of dynamically interacting dark energy (IDE) models. We firstly summarize two dynamical mechanisms in which different background evolutions can exert considerable effects on the 21-cm line global signal. One is the decoupling time of Compton scattering heating, the other stems from the direct change of optical depth due to the different expansion rate of the Universe. After that, we investigate the IDE models to illustrate the tension between the results of Experiment to Detect the Global Epoch of reionization Signature (EDGES) and other experiments. To apply the analyses of these two mechanisms to IDE models, we find that only the optical depth can be significantly changed. Accordingly, in order to relieve the tension by including the effects of the decoupling time of Compton scattering heating, we deduce a possible evolution form for the Hubble parameter within IDE that begins at an early stage around $z\sim 100$ and then smoothly evolves to a value at $z\sim 17$ which is smaller than that obtained in the standard paradigm. Eventually, we fulfill this scenario by adding an early dark energy dominated stage to the cosmological paradigm described by IDE models, which can alleviate the tension between EDGES and other cosmological observations but cannot completely solve it.

We have conducted a large-field simultaneous survey of $^{12}$CO, $^{13}$CO, and C$^{18}$O $J=1-0$ emission toward the Orion A giant molecular cloud (GMC) with a sky coverage of $\sim$ 4.4 deg$^2$ using the PMO-13.7 m millimeter-wavelength telescope. We use the probability distribution function of the column density (N-PDF) to investigate the distribution of molecular hydrogen in the Orion A GMC. The H$_2$ column density, derived from the $^{13}$CO emission, of the GMC is dominated by log-normal distribution in the range from $\sim$4$\times10^{21}$ to $\sim$1.5$\times10^{23}$ cm$^{-2}$ with excesses both at the low-density and high-density ends. The excess of the low-density end is possibly caused by an extended and low-temperature ($\sim$10 K) component with velocities in the range of 5$-$8 km s$^{-1}$. Compared with the northern sub-regions, the southern sub-regions of the Orion A GMC contain less gas with column density in $N_{H_2} > 1.25\times 10^{22}\ \rm{cm}^{-2}$. The dispersions of the N-PDFs of the sub-regions are found to correlate with the evolutionary stages of the clouds across the Orion A GMC. The structure hierarchy of Orion A GMC is explored with the DENDROGRAM algorithm, and it is found that the GMC is composed of two branches. All structures except one in the tree have virial parameters less than 2, indicating self-gravity is important on the spatial scales from $\sim$0.3 to $\sim$4 pc. Although power-laws and departures from log-normal distributions are found at the high-density end of N-PDFs of active star-forming regions, the N-PDFs of structures in the Orion A GMC are predominantly log-normal on scales from R$\sim$0.4 to 4 pc.

Impact craters, as "lunar fossils", are the most dominant lunar surface features and occupy most of the Moon's surface. Their formation and evolution record the history of the Solar System. Sixty years of triumphs in the lunar exploration projects accumulated a large amount of lunar data. Currently, there are 9137 existing recognized craters. However, only 1675 of them have been determined age, which is obviously not satisfactory to reveal the evolution of the Moon. Identifying craters is a challenging task due to their enormous difference in size, large variations in shape and vast presence. Furthermore, estimating the age of craters is extraordinarily difficult due to their complex and different morphologies. Here, in order to effectively identify craters and estimate their age, we convert the crater identification problem into a target detection task and crater age estimation into a taxonomy structure. From an initial small number of available craters, we progressively identify craters and estimate their age from Chang'E data by transfer learning (TL) using deep neural networks. For comprehensive identification of multi-scale craters, a two-stage craters detection approach is developed. Thus 117240 unrecognized lunar craters that range in diameter from 532 km to 1 km are identified. Then, a two-stage classification approach is developed to estimate the age of craters by simultaneously extracting their morphological features and stratigraphic information. The age of 79243 craters larger than 3 km in diameter is estimated. These identified and aged craters throughout the mid and low-latitude regions of the Moon are crucial for reconstructing the dynamic evolution process of the Solar System.

The Askaryan Radio Array is an ultra-high energy (UHE, $>10^{17}$\,eV) neutrino detector designed to observe neutrinos by searching for the radio waves emitted by the relativistic byproducts of neutrino-nucleon interactions in Antarctic ice. In this paper, we present constraints on the diffuse flux of ultra-high energy neutrinos between $10^{16}-10^{21}$\,eV resulting from a search for neutrinos in two complementary analyses, both analyzing four years of data (2013-2016) from the two deep stations (A2, A3) operating at that time. We place a 90\% CL upper limit on the diffuse neutrino flux at $10^{18}$ eV of $EF(E)=3.9\times10^{-16}$~$\textrm{cm}^2$$\textrm{s}^{-1}$$\textrm{sr}^{-1}$. The analysis leverages more than quadruple the livetime of the previous ARA limit, and at this energy is the most sensitive reported by an in-ice radio neutrino detector by more than a factor of four. Looking forward, ARA has yet to analyze more than half of its archived data, and expects that dataset to double in size in the next three years as livetime accumulates in the full five-station array. With this additional livetime, we project ARA will have world-leading sensitivity to neutrinos above $10^{17}$ eV by 2022.

We present CO(J=1-0) and/or CO(J=2-1) spectroscopy for 31 galaxies selected from the ongoing MaNGA survey, obtained with multiple telescopes. This sample is combined with CO observations from the literature to study the correlation of the CO luminosities ($L_{\rm CO(1-0)}$) with the mid-infrared luminosities at 12 ($L_{12 \mu m}$) and 22 $\mu$m ($L_{\rm 22 \mu m}$), as well as the dependence of the residuals on a variety of galaxy properties. The correlation with $L_{\rm 12 \mu m}$ is tighter and more linear, but galaxies with relatively low stellar masses and blue colors fall significantly below the mean $L_{\rm CO(1-0)}-L_{\rm 12\mu m}$ relation. We propose a new estimator of the CO(1-0) luminosity (and thus the total molecular gas mass) that is a linear combination of three parameters: $L_{\rm 12 \mu m}$, $M_\ast$ and $g-r$. We show that, with a scatter of only 0.18 dex in log $(L_{\rm CO(1-0)})$, this estimator provides unbiased estimates for galaxies of different properties and types. An immediate application of this estimator to a compiled sample of galaxies with only CO(J=2-1) observations yields a distribution of the CO(J=2-1) to CO(J=1-0) luminosity ratios ($R21$) that agrees well with the distribution of real observations, in terms of both the median and the shape. Application of our estimator to the current MaNGA sample reveals a gas-poor population of galaxies that are predominantly early-type and show no correlation between molecular gas-to-stellar mass ratio and star formation rate, in contrast to gas-rich galaxies. We also provide alternative estimators with similar scatters, based on $r$ and/or $z$ band luminosities instead of $M_\ast$. These estimators serve as cheap and convenient $M_{\rm mol}$ proxies to be potentially applied to large samples of galaxies, thus allowing statistical studies of gas-related processes of galaxies.

The small-scale crisis is one of the most outstanding puzzles in modern cosmology and astrophysics. It may imply a suppression of matter perturbation at small scale. In this work, by taking into account of the gravitational effects from the non-equilibrium production of DM, we propose a new mechanism, which can realize such desired suppression and alleviate the crisis within the framework of cold dark matter (DM) and simple inflation. Moreover, in this new mechanism, we establish a novel relation between the particle mass of DM, $m_\chi$, and the critical scale of suppression, $k_\star^{-1}$. As $k_\star^{-1}$ will be further constrained in future astrophysical observations, $m_\chi$ can be constrained accordingly with this relation. It thus provides a new method in complementary to other existing strategies of determining $m_\chi$. Furthermore, to illustrate our theoretical prediction, we consider a suppression at $k_\star^{-1}=1~\text{kpc}$ that can partially alleviate the small-scale crisis, and obtain $m_\chi=2.2~\text{MeV}$ for realizing such suppression. Then we plot the power spectrum of linear matter perturbation for this case, and illustrate a salient feature of the suppression that can serves a smoking-gun signature of this new mechanism in future observations.

Low-resolution spectra are proved competitive to high-resolution spectra in determining many stellar labels at comparable precision. It is useful to consider the spectral information content when assessing the capability of a stellar spectrum in deriving precise stellar labels. In this work, we quantify the information content brought by the LAMOST-II medium-resolution spectroscopic survey (MRS) using the gradient spectra and the coefficients-of-dependence (CODs). In general, the wavelength coverage of the MRS well constrains the stellar labels but the sensitivities of different stellar labels vary with spectral types and metallicity of the stars of interest and, therefore, affect the performance of the stellar label determination from the MRS spectra. Applying the SLAM to the synthetic spectra which mimic the MRS data, we find the precision of the fundamental stellar parameters Teff, logg and [M/H] are better when combining both the blue and red bands of the MRS. This is especially important for warm stars since the H$\alpha$ line located in the red part plays a more important role in determining the effective temperature for warm stars. With blue and red parts together, we are able to reach similar performance to the low-resolution spectra except for warm stars. However, at [M/H]$\sim-2.0$ dex, the uncertainties of fundamental stellar labels estimated from MRS are substantially larger than those from low-resolution spectra. We also tested the uncertainties of Teff, logg and [M/H] of from MRS data induced from the radial velocity mismatch and find that a mismatch of about 1 km s$^{-1}$, which is typical for LAMOST MRS data, would not significantly affect the stellar label estimates. At last, reference precision limits are calculated using synthetic gradient spectra, according to which we expect abundances of at least 17 elements to be measured precisely from MRS spectra.

We have used two methods to search for surviving companions of Type Ia supernova progenitors in three Balmer-dominated supernova remnants (SNRs) in the Large Magellanic Cloud: 0519-69.0, 0505-67.9 (DEM L71), and 0548-70.4. In the first method, we use the Hubble Space Telescope photometric measurements of stars to construct color-magnitude diagrams (CMDs), and compare positions of stars in the CMDs with those expected from theoretical post-impact evolution of surviving main sequence or helium star companions. No obvious candidates of surviving companion are identified in this photometric search. Future models for surviving red giant companions or with different explosion mechanisms are needed for thorough comparisons with these observations in order to make more definitive conclusions. In the second method, we use Multi-Unit Spectroscopic Explorer (MUSE) observations of 0519-69.0 and DEM L71 to carry out spectroscopic analyses of stars in order to use large peculiar radial velocities as diagnostics of surviving companions. We find a star in 0519-69.0 and a star in DEM L71 moving at radial velocities of 182 $\pm$ 0 km s$^{-1}$ and 213 $\pm$ 0 km s$^{-1}$, more than 2.5$\sigma$ from the mean radial velocity of the underlying stellar population, 264 km s$^{-1}$ and 270 km s$^{-1}$, respectively. These stars need higher-quality spectra to investigate their abundances and rotation velocities to determine whether they are indeed surviving companions of the SN progenitors.

The 'chromosome map' diagram (ChM) proved a successful tool to identify and characterize multiple populations (MPs) in 59 Galactic Globular Clusters (GCs). Here, we construct ChMs for 11 GCs of both Magellanic Clouds (MCs) and with different ages to compare MPs in Galactic and extra-Galactic environments, and explore whether this phenomenon is universal through 'place' and 'time'. MPs are detected in five clusters. The fractions of 1G stars, ranging from about 50% to more than 80%, are significantly higher than those observed in Galactic GCs with similar present-day masses. By considering both Galactic and MC clusters, the fraction of 1G stars exhibits: (i) a strong anti-correlation with the present-day mass, and (ii) with the present-day mass of 2G stars; (iii) a mild anti-correlation with 1G present-day mass. All Galactic clusters without MPs have initial masses smaller than ~1.5 10^5 solar masses but a mass threshold governing the occurrence of MPs seems challenged by massive simple-population MC GCs; (iv) Milky Way clusters with large perigalactic distances typically host larger fractions of 1G stars, but the difference disappears when we use initial cluster masses. These facts are consistent with a scenario where the stars lost by GCs mostly belong to the 1G. By exploiting recent work based on Gaia, half of the known Type II GCs appear clustered in a distinct region of the integral of motions space, thus suggesting a common progenitor galaxy. Except for these Type II GCs,we do not find any significant difference in the MPs between clusters associated with different progenitors.

As China's first X-ray astronomical satellite, the Hard X-ray Modulation Telescope (HXMT), which was dubbed as Insight-HXMT after the launch on June 15, 2017, is a wide-band (1-250 keV) slat-collimator-based X-ray astronomy satellite with the capability of all-sky monitoring in 0.2-3 MeV. It was designed to perform pointing, scanning and gamma-ray burst (GRB) observations and, based on the Direct Demodulation Method (DDM), the image of the scanned sky region can be reconstructed. Here we give an overview of the mission and its progresses, including payload, core sciences, ground calibration/facility, ground segment, data archive, software, in-orbit performance, calibration, background model, observations and some preliminary results.

We present a detailed timing study of the brightest persistent X-ray source Sco X-1 using the data collected by the Hard X-ray Modulation Telescope ($Insight$-HXMT) from July 2017 to August 2018. A complete $Z$-track hardness-intensity diagram (HID) is obtained. The normal branch oscillations (NBOs) at $\sim$ 6 Hz in the lower part of the normal branch (NB) and the flare branch oscillations (FBOs) at $\sim$ 16 Hz in the beginning part of the flaring branch (FB) are found in observations with the Low Energy X-ray Telescope (LE) and the Medium Energy X-ray Telescope (ME) of $Insight$-HXMT, while the horizontal branch oscillations (HBOs) at $\sim$ 40 Hz and the kilohertz quasi-periodic oscillations (kHz QPOs) at $\sim$ 800 Hz are found simultaneously up to 60 keV for the first time on the horizontal branch (HB) by the High Energy X-ray Telescope (HE) and ME. We find that for all types of the observed QPOs, the centroid frequencies are independent of energy, while the root mean square (rms) increases with energy; the centroid frequencies of both the HBOs and kHz QPOs increase along the $Z$-track from the top to the bottom of the HB; and the NBOs show soft phase lags increasing with energy. A continuous QPO transition from the FB to NB in $\sim$ 200 s are also detected. Our results indicate that the non-thermal emission is the origin of all types of QPOs, the innermost region of the accretion disk is non-thermal in nature, and the corona is nonhomogeneous geometrically.

The evidences for the influence of thermonuclear (type-I) X-ray bursts upon the surrounding environments in neutron star low-mass X-ray binaries (LMXB) were detected previously via spectral and timing analyses. Benefitting from a broad energy coverage of Insight-HXMT, we analyze one photospheric radius expansion (PRE) burst, and find an emission excess at soft X-rays. Our spectral analysis shows that, such an excess is not likely relevant to the disk reflection induced by the burst emission and can be attributed to an enhanced pre-burst/persistent emission. We find that the burst and enhanced persistent emissions sum up to exceed Eddington luminosity by $\sim$ 40 percentages. We speculate that the enhanced emission is from a region beyond the PRE radius, or through the Comptonization of the corona.

The detection of clumps(cores) in molecular clouds is an important issue in sub-millimetre astronomy. However, the completeness of the identification and the accuracy of the returned parameters of the automated clump identification algorithms are still not clear by now. In this work, we test the performance and bias of the GaussClumps, ClumpFind, Fellwalker, Reinhold, and Dendrograms algorithms in identifying simulated clumps. By designing the simulated clumps with various sizes, peak brightness, and crowdedness, we investigate the characteristics of the algorithms and their performance. In the aspect of detection completeness, Fellwalker, Dendrograms, and Gaussclumps are the first, second, and third best algorithms, respectively. The numbers of correct identifications of the six algorithms gradually increase as the size and SNR of the simulated clumps increase and they decrease as the crowdedness increases. In the aspect of the accuracy of retrieved parameters, Fellwalker and Dendrograms exhibit better performance than the other algorithms. The average deviations in clump parameters for all algorithms gradually increase as the size and SNR of clumps increase. Most of the algorithms except Fellwalker exhibit significant deviation in extracting the total flux of clumps. Taken altogether, Fellwalker, Gaussclumps, and Dendrograms exhibit the best performance in detection completeness and extracting parameters. The deviation in virial parameter for the six algorithms is relatively low. When applying the six algorithms to the clump identification for the Rosette molecular cloud, ClumpFind1994, ClumpFind2006, Gaussclumps, Fellwalker, and Reinhold exhibit performance that is consistent with the results from the simulated test.

Corona cooling was detected previously from stacking a series of short type-I bursts occurred during the low/had state of atoll outburst. Type-I bursts are hence regarded as sharp probe to our better understanding on the basic property of the corona. The launch of the first Chinese X-ray satellite Insight-HXMT has large detection area at hard X-rays which provide almost unique chance to move further in this research field. We report the first detection of the corona cooling by Insight-HXMT from single short type-I burst showing up during {\bf flare} of 4U 1636-536. This type-I X-ray burst has a duration of $\sim$13 seconds and hard X-ray shortage is detected with significance 6.2~$\sigma$ in 40-70 keV. A cross-correlation analysis between the lightcurves of soft and hard X-ray band, shows that the corona shortage lag the burst emission by 1.6 $\pm$1.2~s. These results are consistent with those derived previously from stacking a large amount of bursts detected by RXTE/PCA within a series of {\bf flares} of 4U 1636-536. Moreover, the broad bandwidth of Insight-HXMT allows as well for the first time to infer the burst influence upon the continuum spectrum via performing the spectral fitting of the burst, which ends up with the finding that hard X-ray shortage appears at around 40 keV in the continuum spectrum. These results suggest that the evolution of the corona along with the outburst{\bf /flare} of NS XRB may be traced via looking into a series of embedded type-I bursts by using Insight-HXMT.

The Hard X-ray Modulation Telescope HXMT or also dubbed as Insight HXMT is China s first astronomical satellite . It was launched on 15 th June 2017 in JiuQuan, China and is currently in service smoothly. It was designed to perform point ing , scan ning and gamma ray burst (GRB) observations and , based on the Direct Demodulation Method (DDM), the image of the scanned sky region can be reconstructed. Here we introduce the mission and its progresses in aspects of payload, core sciences, ground calibration/facility , ground segment, data archive, software, in orbit performance, calibration, background model, observations and preliminary results .

The Hard X-ray Modulation Telescope (HXMT) named Insight is China's first X-ray astronomical satellite. The Low Energy X-ray Telescope (LE) is one of its main payloads onboard. The detectors of LE adopt swept charge device CCD236 with L-shaped transfer electrodes. Charges in detection area are read out continuously along specific paths, which leads to a time response distribution of photons readout time. We designed a long exposure readout mode to measure the time response distribution. In this mode, CCD236 firstly performs exposure without readout, then all charges generated in preceding exposure phase are read out completely. Through analysis of the photons readout time in this mode, we obtained the probability distribution of photons readout time.

We present a study of timing properties of the accreting pulsar 2S 1417-624 observed during its 2018 outburst, based on Swift/BAT, Fermi/GBM, Insight-HXMT and NICER observations. We report a dramatic change of the pulse profiles with luminosity. The morphology of the profile in the range 0.2-10.0keV switches from double to triple peaks at $\sim2.5$ $\rm \times 10^{37}{\it D}_{10}^2\ erg\ s^{-1}$ and from triple to quadruple peaks at $\sim7$ $\rm \times 10^{37}{\it D}_{10}^2\ erg\ s^{-1}$. The profile at high energies (25-100keV) shows significant evolutions as well. We explain this phenomenon according to existing theoretical models. We argue that the first change is related to the transition from the sub to the super-critical accretion regime, while the second to the transition of the accretion disc from the gas-dominated to the radiation pressure-dominated state. Considering the spin-up as well due to the accretion torque, this interpretation allows to estimate the magnetic field self-consistently at $\sim7\times 10^{12}$G.

The long-term evolution of the centroid energy of the CRSF in Her X-1 is still a mystery. We report a new measurement from a campaign between {\sl Insight}-HXMT and {\sl NuSTAR} performed in February 2018. Generally, the two satellites show well consistent results of timing and spectral properties. The joint spectral analysis confirms that the previously observed long decay phase has ended, and that the line energy instead keeps constant around 37.5 keV after flux correction.

We characterize the kinematic and chemical properties of $\sim$3,000 Sagittarius (Sgr) stream stars, including K-giants, M-giants, and BHBs, select from SEGUE-2, LAMOST, and SDSS separately in Integrals-of-Motion space. The orbit of Sgr stream is quite clear from the velocity vector in $X$-$Z$ plane. Stars traced by K-giants and M-giants present the apogalacticon of trailing steam is $\sim$ 100 kpc. The metallicity distributions of Sgr K-, M-giants, and BHBs present that the M-giants are on average the most metal-rich population, followed by K-giants and BHBs. All of the K-, M-giants, and BHBs indicate that the trailing arm is on average more metal-rich than leading arm, and the K-giants show that the Sgr debris is the most metal-poor part. The $\alpha$-abundance of Sgr stars exhibits a similar trend with the Galactic halo stars at lower metallicity ([Fe/H] $<\sim$ $-$1.0 dex), and then evolve down to lower [$\alpha$/Fe] than disk stars at higher metallicity, which is close to the evolution pattern of $\alpha$-element of Milky Way dwarf galaxies. We find $V_Y$ and metallicity of K-giants have gradients along the direction of line-of-sight from the Galactic center in $X$-$Z$ plane, and the K-giants show that $V_Y$ increases with metallicity at [Fe/H] $>\sim-$1.5 dex. After dividing the Sgr stream into bright and faint stream according to their locations in equatorial coordinate, the K-giants and BHBs show that the bright and faint stream present different $V_Y$ and metallicities, the bright stream is on average higher in $V_Y$ and metallicity than the faint stream.

We report on analysis of observations of the bright transient X-ray pulsar \src obtained during its 2017-2018 giant outburst with Insight-HXMT, \emph{NuSTAR}, and \textit{Swift} observatories. We focus on the discovery of a sharp state transition of the timing and spectral properties of the source at super-Eddington accretion rates, which we associate with the transition of the accretion disk to a radiation pressure dominated (RPD) state, the first ever directly observed for magnetized neutron star. This transition occurs at slightly higher luminosity compared to already reported transition of the source from sub- to super-critical accretion regime associate with onset of an accretion column. We argue that this scenario can only be realized for comparatively weakly magnetized neutron star, not dissimilar to other ultra-luminous X-ray pulsars (ULPs), which accrete at similar rates. Further evidence for this conclusion is provided by the non-detection of the transition to the propeller state in quiescence which strongly implies compact magnetosphere and thus rules out magnetar-like fields.

The origin of the star forming main sequence ( i.e., the relation between star formation rate and stellar mass, globally or on kpc-scales; hereafter SFMS) remains a hotly debated topic in galaxy evolution. Using the ALMA-MaNGA QUEnching and STar formation (ALMaQUEST) survey, we show that for star forming spaxels in the main sequence galaxies, the three local quantities, star-formation rate surface density (\sigsfr), stellar mass surface density (\sigsm), and the \h2~mass surface density (\sigh2), are strongly correlated with one another and form a 3D linear (in log) relation with dispersion. In addition to the two well known scaling relations, the resolved SFMS (\sigsfr~ vs. \sigsm) and the Schmidt-Kennicutt relation (\sigsfr~ vs. \sigh2; SK relation), there is a third scaling relation between \sigh2~ and \sigsm, which we refer to as the `molecular gas main sequence' (MGMS). The latter indicates that either the local gas mass traces the gravitational potential set by the local stellar mass or both quantities follow the underlying total mass distributions. The scatter of the resolved SFMS ($\sigma \sim 0.25$ dex) is the largest compared to those of the SK and MGMS relations ($\sigma \sim$ 0.2 dex). A Pearson correlation test also indicates that the SK and MGMS relations are more strongly correlated than the resolved SFMS. Our result suggests a scenario in which the resolved SFMS is the least physically fundamental and is the consequence of the combination of the SK and the MGMS relations.

The CFHT Large Area U-band Deep Survey (CLAUDS) uses data taken with the MegaCam mosaic imager on CFHT to produce images of 18.60 deg2 with median seeing of FWHM=0.92 arcsec and to a median depth of U = 27.1 AB (5 sigma in 2 arcsec apertures), with selected areas that total 1.36 deg2 reaching a median depth of U=27.7 AB. These are the deepest U-band images assembled to date over this large an area. These data are located in four fields also imaged to comparably faint levels in grizy and several narrowband filters as part of the Hyper Suprime-Cam (HSC) Subaru Strategic Program (HSC-SSP). These CFHT and Subaru datasets will remain unmatched in their combination of area and depth until the advent of the Large Synoptic Survey Telescope (LSST). This paper provides an overview of the scientific motivation for CLAUDS and gives details of the observing strategy, observations, data reduction, and data merging with the HSC-SSP. Three early applications of these deep data are used to illustrate the potential of the dataset: deep U-band galaxy number counts, z~3 Lyman break galaxy (LBG) selection, and photometric redshifts improved by adding CLAUDS U to the Subaru HSC grizy photometry.

Within the standard model of hierarchical galaxy formation in a {\Lambda}CDM Universe, the environment of galaxies is expected to play a key role in driving galaxy formation and evolution. In this paper we investigate whether and how the gas metallicity and the star formation surface density ({\Sigma}_SFR) depend on galaxy environment. To this end we analyse a sample of 1162 local, star-forming galaxies from the galaxy survey Mapping Nearby Galaxies at APO (MaNGA). Generally, both parameters do not show any significant dependence on environment. However, in agreement with previous studies, we find that low-mass satellite galaxies are an exception to this rule. The gas metallicity in these objects increases while their {\Sigma}SFR decreases slightly with environmental density. The present analysis of MaNGA data allows us to extend this to spatially resolved properties. Our study reveals that the gas metallicity gradients of low-mass satellites flatten and their {\Sigma}SFR gradients steepen with increasing environmental density. By extensively exploring a chemical evolution model, we identify two scenarios that are able to explain this pattern: metal-enriched gas accretion or pristine gas inflow with varying accretion timescales. The latter scenario better matches the observed {\Sigma}SFR gradients, and is therefore our preferred solution. In this model, a shorter gas accretion timescale at larger radii is required. This suggests that 'outside-in quenching' governs the star formation processes of low-mass satellite galaxies in dense environments.

We report on analysis of englacial radio-frequency (RF) pulser data received over horizontal baselines of 1--5 km, based on broadcasts from two sets of transmitters deployed to depths of up to 1500 meters at the South Pole. First, we analyze data collected usingtwo RF bicone transmitters 1400 meters below the ice surface, and frozen into boreholes drilled for the IceCube experiment in 2011. Additionally, in Dec., 2018, a fat-dipole antenna, fed by one of three high-voltage (~1 kV), fast (~(1-5 ns)) signal generators was lowered into the 1700-m deep icehole drilled for the South Pole Ice Core Experiment (SPICE), approximately 3 km from the geographic South Pole. Signals from transmitters were recorded on the five englacial multi-receiver ARA stations, with receiver depths between 60--200 m. We confirm the long, >1 km RF electric field attenuation length, test our observed signal arrival timing distributions against models, and measure birefringent asymmetries at the 0.15% level.

The presence of multiple populations (MPs) in almost all globular clusters (GCs) older than ~10 Gyr, has caught lots of attention. Recently, cumulative evidence indicates that extragalactic GCs that are older than 2 Gyr, seem to also harbor MPs, however, those that are younger than 2 Gyr do not. These observations seem to imply that age is a primary property that controls the presence of MPs in star clusters. However, because of the lack of studies of {intermediate-age (~2--6 Gyr-old),} low mass clusters, it is unclear if the cluster mass, in addition to age, also plays a role in the occurrence of MPs. In this work, we studied a $\sim$4 Gyr-old, low mass (~23,000 Msun) cluster, Lindsay-113, in the Small Magellanic Cloud. Using Hubble Space Telescope photometry, we find that the width of the red-giant branch in this cluster, when measured in a specific color index that is sensitive to star-to-star chemical variations, {can be adequately explained by a "simple" stellar population model with some possible noises contributed from measurement uncertainty, photometric artifact, as well as differential reddening. The comparison of observations with predictions from synthetic spectra indicates that the internal chemical spread in nitrogen abundance, which is a signature of MPs, would not exceed 0.2 dex. Since Lindsay 113 is significantly older than other GCs with MPs, we suggest that the onset of MPs is likely determined by the combination of cluster age and mass.

We present a detailed analysis of the projected stellar rotational velocities of the well-separated double main sequence (MS) in the young, $\sim200$Myr-old Milky Way open cluster NGC 2287 and suggest that stellar rotation may drive the split MSs in NGC 2287. We find that the observed distribution of projected stellar rotation velocities could result from a dichotomous distribution of stellar rotation rates. We discuss whether our observations may reflect the effects of tidal locking affecting a fraction of the cluster's member stars in stellar binary systems. The slow rotators are likely stars that initially rotated rapidly but subsequently slowed down through tidal locking induced by low-mass-ratio binary systems. However, the cluster may have a much larger population of short-period binaries than is usually seen in the literature, with relatively low secondary masses.

Stars spend most of their lifetimes on the `main sequence' (MS) in the Hertzsprung--Russell diagram. The obvious double MSs seen in the equivalent color--magnitude diagrams characteristic of Milky Way open clusters pose a fundamental challenge to our traditional understanding of star clusters as `single stellar populations.' The clear MS bifurcation of early-type stars with masses greater than $\sim1.6 M_\odot$ is thought to result from a range in the stellar rotation rates. However, direct evidence connecting double MSs to stellar rotation properties has yet to emerge. Here, we show through analysis of the projected stellar rotational velocities ($v\sin i$, where $i$ represents the star's inclination angle) that the well-separated double MS in the young, $\sim200Myr$-old Milky Way open cluster NGC 2287 is tightly correlated with a dichotomous distribution of stellar rotation rates. We discuss whether our observations may reflect the effects of tidal locking affecting a fraction of the cluster's member stars in stellar binary systems. We show that the slow rotators could potentially be initially rapidly rotating stars that have been slowed down by tidal locking by a low mass-ratio companion in a cluster containing a large fraction of short-period, low-mass-ratio binaries. This demonstrates that stellar rotation drives the split MSs in young, $\lessapprox 300$Myr-old star clusters. However, special conditions, e.g., as regards the mass-ratio distribution, might be required for this scenario to hold.

Resolving spatially-varying exoplanet features from single-point light curves is essential for determining whether Earth-like worlds harbor geological features and/or climate systems that influence habitability. To evaluate the feasibility and requirements of this spatial feature resolving problem, we present an analysis of multi-wavelength single-point light curves of Earth, where it plays the role of a proxy exoplanet. Here, ~10,000 DSCOVR/EPIC frames collected over a two-year period were integrated over the Earth's disk to yield a spectrally-dependent point source and analyzed using singular value decomposition. We found that, between the two dominant principal components (PCs), the second PC contains surface-related features of the planet, while the first PC mainly includes cloud information. We present the first two-dimensional (2D) surface map of Earth reconstructed from light curve observations without any assumptions of its spectral properties. This study serves as a baseline for reconstructing the surface features of Earth-like exoplanets in the future.

In this work, we report the in-orbit demonstration of X-ray pulsar navigation with Insight-Hard X-ray Modulation Telescope (Insight-HXMT), which was launched on Jun. 15th, 2017. The new pulsar navigation method 'Significance Enhancement of Pulse-profile with Orbit-dynamics' (SEPO) is adopted to determine the orbit with observations of only one pulsar. In this test, the Crab pulsar is chosen and observed by Insight-HXMT from Aug. 31th to Sept. 5th in 2017. Using the 5-day-long observation data, the orbit of Insight-HXMT is determined successfully with the three telescopes onboard - High Energy X-ray Telescope (HE), Medium Energy X-ray Telescope (ME) and Low Energy X-ray Telescope (LE) - respectively. Combining all the data, the position and velocity of the Insight-HXMT are pinpointed to within 10 km (3 sigma) and 10 m/s (3 sigma), respectively.

The EDGES collaboration's observation of an anomalously strong 21 cm absorption feature around the cosmic dawn era has energised the cosmological community by suggesting a novel signature of dark matter in the cooling of cosmic hydrogen. In a recent letter we have argued that by virtue of the ability to mediate cooling processes whilst in the condensed phase, a small amount of axion dark matter can explain these observations within the context of standard models of axions and axion-like particles. These axions and axion-like particles (ALPs) can thermalize through gravitational self-interactions and so eventually form a Bose-Einstein condensate (BEC), whereupon large-scale long-range correlation can produce experimentally observable signals such as these. In this context the EDGES best-fit result favours an axion-like-particles mass in the (10, 450) meV range, which can be compressed for the QCD axion to (100, 450) meV in the absence of fine tuning. Future experiments and galaxy surveys, particularly the International Axion Observatory (IAXO) and EUCLID, should have the capability to directly test this scenario. In this paper, we will explore this mechanism in detail and give more thorough computational details of certain key points.

The Askaryan Radio Array (ARA) is an ultra-high energy (UHE) neutrino telescope at the South Pole consisting of an array of radio antennas aimed at detecting the Askaryan radiation produced by neutrino interactions in the ice. Currently, the experiment has five stations in operation that have been deployed in stages since 2012. This contribution focuses on the development of a search for a diffuse flux of neutrinos in two ARA stations (A2 and A3) from 2013-2016. A background of $\sim 0.01-0.02$ events is expected in one station in each of two search channels in horizontal- and vertical-polarizations. The expected new constraints on the flux of ultra-high energy neutrinos based on four years of analysis with two stations improve on the previous limits set by ARA by a factor of about two. The projected sensitivity of ARA's five-station dataset is beginning to be competitive with other neutrino telescopes at high energies near $10^{10.5}\,$GeV.

The recently discovered neutron star transient Swift J0243.6+6124 has been monitored by {\it the Hard X-ray Modulation Telescope} ({\it Insight-\rm HXMT). Based on the obtained data, we investigate the broadband spectrum of the source throughout the outburst. We estimate the broadband flux of the source and search for possible cyclotron line in the broadband spectrum. No evidence of line-like features is, however, found up to $\rm 150~keV$. In the absence of any cyclotron line in its energy spectrum, we estimate the magnetic field of the source based on the observed spin evolution of the neutron star by applying two accretion torque models. In both cases, we get consistent results with $B\rm \sim 10^{13}~G$, $D\rm \sim 6~kpc$ and peak luminosity of $\rm >10^{39}~erg~s^{-1}$ which makes the source the first Galactic ultraluminous X-ray source hosting a neutron star.

We construct a large halo K-giant sample by combining the positions, distances, radial velocities, and metallicities of over 13,000 LAMOST DR5 halo K giants with the Gaia DR2 proper motions, which covers a Galactocentric distance range of 5-120 kpc. Using a position-velocity clustering estimator (the 6Distance), we statistically quantify the presence of position-velocity substructure at high significance: K giants have more close pairs in position-velocity space than a smooth stellar halo. We find that the amount of substructure in the halo increases with increasing distance and metallicity. With a percolation algorithm named friends-of-friends (FoF) to identify groups, we identify members belonging to Sagittarius (Sgr) Streams, Monoceros Ring, Virgo overdensity, Hercules-Aquila Cloud, Orphan Streams and other unknown substructures and find that the Sgr streams account for a large part of grouped stars beyond 20 kpc and enhance the increase of substructure with distance and metallicity. For the first time, we identify spectroscopic members of Monoceros Ring in the south and north Galactic hemisphere, which presents a rotation of about 185 km s^{-1} and mean metallicity is -0.66 dex.

(Abridged) We perform a detailed study of the main optical activity indicators (Ca II H & K, Balmer lines, Na I D$_{\rm 1}$ D$_{\rm 2}$, and He I D$_{\rm 3}$) measured for the Sun using the data provided by the HARPS-N solar-telescope feed at the Telescopio Nazionale Galileo. The value of the solar rotation period is found in all the activity indicators, with the only exception being H$\delta$. The derived values vary from 26.29 days (H$\gamma$ line) to 31.23 days (He I). From an analysis of sliding periodograms we find that in most of the activity indicators the spectral power is split into several "bands" of periods around 26 and 30 days, that might be explained by the migration of active regions between the equator and a latitude of $\sim$ 30$^{\circ}$, spot evolution or a combination of both effects. In agreement with previous works a typical lifetime of active regions of $\sim$ ten rotation periods is inferred from the pooled variance diagrams. We find that H$\alpha$, H$\beta$, H$\gamma$, H$\epsilon$, and He I show a significant correlation with the S index. Significant correlations between the contrast, bisector span, and the heliocentric radial velocity with the activity indexes are also found. We show that the full width at half maximum, the bisector, and the disc-integrated magnetic field correlate with the radial velocity variations. The correlation of the S index and H$\alpha$ changes with time, increasing with larger sun spot numbers and solar irradiance. A similar tendency with the S index - radial velocity correlation is also present in the data. Our results are consistent with a scenario in which higher activity favours the correlation between the S index and the H$\alpha$ activity indicators and between the S index and radial velocity variations.

Cassini/UVIS FUV observations of stellar occultations at Titan are well suited for probing its atmospheric composition and structure. However, due to instrument pointing motion, only five out of tens of observations have been analyzed. We present an innovative retrieval method that corrects for the effect of pointing motion by forward modeling the Cassini/UVIS instrument response function with the pointing motion value obtained from the SPICE C-kernel along the spectral dimension. To illustrate the methodology, an occultation observation made during flyby T52 is analyzed, when the Cassini spacecraft had insufficient attitude control. A high-resolution stellar model and an instrument response simulator that includes the position of the point source on the detector are used for the analysis of the pointing motion. The Markov Chain Monte-Carlo method is used to retrieve the line-of-sight abundance profiles of eleven species (CH4, C2H2, C2H4, C2H6, C4H2, C6H6, HCN, C2N2, HC3N, C6N2 and haze particles) in the spectral vector fitting process. We obtain tight constraints on all of the species aside from C2H6, C2N2 and C6N2, for which we only retrieved upper limits. This is the first time that the T52 occultation was used to derive abundances of major hydrocarbon and nitrile species in Titan's upper and middle atmosphere, as pointing motion prohibited prior analysis. With this new method, nearly all of the occultations obtained over the entire Cassini mission could yield reliable profiles of atmospheric composition, allowing exploration of Titan's upper atmosphere over seasons, latitudes, and longitudes.

It has been suggested that the Z-mode instability driven by energetic electrons with a loss-cone type velocity distribution is one candidate process behind the continuum and zebra pattern of solar type-IV radio bursts. Both the temperature of background plasma ($T_0$) and the energy of energetic electrons ($v_e$) are considered to be important to the variation of the maximum growth rate ($\gamma_{max}$). Here we present a detailed parameter study on the effect of $T_0$ and $v_e$, within a regime of the frequency ratio ($10 \leq \frac{\omega_{pe}}{\Omega_{ce}} \leq 30$). In addition to $\gamma_{max}$, we also analyze the effect on the corresponding wave frequency ($\omega^r_{max}$) and propagation angle ($\theta_{max}$). We find that (1) $\gamma_{max}$ in-general decreases with increasing $v_e$, while its variation with $T_0$ is more complex depending on the exact value of $v_e$; (2) with increasing $T_0$ and $v_e$, $\omega^r_{max}$ presents step-wise profiles with jumps separated by gradual or very-weak variations, and due to the warm-plasma effect on the wave dispersion relation $\omega^r_{max}$ can vary within the hybrid band (the harmonic band containing the upper hybrid frequency) and the band higher; (3) the propagation is either perpendicular or quasi-perpendicular, and $\theta_{max}$ presents variations in line with those of $\omega^r_{max}$, as constrained by the resonance condition. We also examine the profiles of $\gamma_{max}$ with $\frac{\omega_{pe}}{\Omega_{ce}}$ for different combinations of $T_0$ and $v_e$ to clarify some earlier calculations which show inconsistent results.

We have conducted a large-field simultaneous survey of $^{12}$CO, $^{13}$CO, and C$^{18}$O $J=1-0$ emission toward the Cassiopeia A (Cas A) supernova remnant (SNR), which covers a sky area of $3.5^{\circ}\times3.1^{\circ}$. The Cas giant molecular cloud (GMC) mainly consists of three individual clouds with masses on the order of $10^4-10^5\ M_{\odot}$. The total mass derived from the $\rm{^{13}CO}$ emission of the GMC is 2.1$\times10^{5}\ M_{\odot}$ and is 9.5$\times10^5\ M_{\odot}$ from the $\rm{^{12}CO}$ emission. Two regions with broadened (6$-$7 km s$^{-1}$) or asymmetric $^{12}$CO line profiles are found in the vicinity (within a 10$'\times10'$ region) of the Cas A SNR, indicating possible interactions between the SNR and the GMC. Using the GAUSSCLUMPS algorithm, 547 $^{13}$CO clumps are identified in the GMC, 54$\%$ of which are supercritical (i.e. $\alpha_{\rm{vir}}<2$). The mass spectrum of the molecular clumps follows a power-law distribution with an exponent of $-2.20$. The pixel-by-pixel column density of the GMC can be fitted with a log-normal probability distribution function (N-PDF). The median column density of molecular hydrogen in the GMC is $1.6\times10^{21}$ cm$^{-2}$ and half the mass of the GMC is contained in regions with H$_2$ column density lower than $3\times10^{21}$ cm$^{-2}$, which is well below the threshold of star formation. The distribution of the YSO candidates in the region shows no agglomeration.

We present the SCUBA-2 850 ${\mu}m$ component of JINGLE, the new JCMT large survey for dust and gas in nearby galaxies, which with 193 galaxies is the largest targeted survey of nearby galaxies at 850 ${\mu}m$. We provide details of our SCUBA-2 data reduction pipeline, optimised for slightly extended sources, and including a calibration model adjusted to match conventions used in other far-infrared data. We measure total integrated fluxes for the entire JINGLE sample in 10 infrared/submillimetre bands, including all WISE, Herschel-PACS, Herschel-SPIRE and SCUBA-2 850 ${\mu}m$ maps, statistically accounting for the contamination by CO(J=3-2) in the 850 ${\mu}m$ band. Of our initial sample of 193 galaxies, 191 are detected at 250 ${\mu}m$ with a $\geq$ 5${\sigma}$ significance. In the SCUBA-2 850 ${\mu}m$ band we detect 126 galaxies with $\geq$ 3${\sigma}$ significance. The distribution of the JINGLE galaxies in far-infrared/sub-millimetre colour-colour plots reveals that the sample is not well fit by single modified-blackbody models that assume a single dust-emissivity index $(\beta)$. Instead, our new 850 ${\mu}m$ data suggest either that a large fraction of our objects require $\beta < 1.5$, or that a model allowing for an excess of sub-mm emission (e.g., a broken dust emissivity law, or a very cold dust component 10 K) is required. We provide relations to convert far-infrared colours to dust temperature and $\beta$ for JINGLE-like galaxies. For JINGLE the FIR colours correlate more strongly with star-formation rate surface-density rather than the stellar surface-density, suggesting heating of dust is greater due to younger rather than older stellar-populations, consistent with the low proportion of early-type galaxies in the sample.

An increasing number of white dwarf stars show atmospheric chemical composition polluted by heavy elements accreted from debris disk material. The existence of such debris disks strongly suggests the presence of one or more planet(s) whose gravitational interaction with rocky planetesimals is responsible for their disruption by tidal effect. The ZZ Ceti pulsator and polluted DAZ white dwarf GD 133 is a good candidate for searching for such a potential planet. We started in 2011 a photometric follow-up of its pulsations. As a result of this work in progress, we used the data gathered from 2011 to 2015 to make an asteroseismological analysis of GD 133, providing the star parameters from a best fit model with $M$/$M_{\odot}$ = 0.630 $\pm$ 0.002, $T_{\rm eff}$ = 12400 K $\pm$ 70 K, log($M_{\rm He}/M$) = -2.00 $\pm$ 0.02, log($M_{\rm H}/M$) = -4.50 $\pm$ 0.02 and determining a rotation period of $\approx$ 7 days.

Using {\sl Gaia} Data Release 2 photometry, we report the detection of extended main-sequence turnoff (eMSTO) regions in the color--magnitude diagrams (CMDs) of the $\sim 14$ Myr-old double clusters $h$ and $\chi$ Persei (NGC 869 and NGC 884). We find that stars with masses below $\sim$1.3 $M_{\odot}$ in both $h$ and $\chi$ Persei populate narrow main sequences (MSs), while more massive stars define the eMSTO, closely mimicking observations of young Galactic and Magellanic Cloud clusters (with ages older than $\sim$30 Myr). Previous studies based on clusters older than $\sim$30 Myr find that rapidly rotating MS stars are redder than slow rotators of similar luminosity, suggesting that stellar rotation may be the main driver of the eMSTO. By combining photometry and projected rotational velocities from the literature of stars in $h$ and $\chi$ Persei, we find no obvious relation between the rotational velocities and colors of non-emission-line eMSTO stars, in contrast with what is observed in older clusters. Similarly to what is observed in Magellanic Cloud clusters, most of the extremely rapidly rotating stars, identified by their strong H$\alpha$ emission lines, are located in the red part of the eMSTOs. This indicates that stellar rotation plays a role in the color and magnitude distribution of MSTO stars. By comparing the observations with simulated CMDs, we find that a simple population composed of coeval stars that span a wide range of rotation rates is unable to reproduce the color spread of the clusters' MSs. We suggest that variable stars, binary interactions, and stellar rotation affect the eMSTO morphology of these very young clusters.

Star-to-star light-element abundance variations, know as multiple stellar populations (MPs), are common in almost all Galactic globular clusters. Recently, MPs have also been detected in a number of massive clusters with ages in excess of 2 Gyr in the Large Magellanic Cloud (LMC), thus indicating that age is likely a control parameter for the presence of MPs. However, to conclusively confirm this notion, additional studies of star clusters in the LMC's `age gap' of 3--6 Gyr are required. Here, we use {\sl Hubble Space Telescope} observations to study the 3 Gyr-old cluster NGC 2121. Compared with so-called `simple' stellar population models, the cluster's red-giant branch exhibits an apparent spread in a specific color index that is sensitive to intrinsic chemical spreads. The observed spread can be explained by an intrinsic spread in nitrogen abundance of $\sim$0.5--1.0 dex}. NGC 2121 has a comparable mass to its younger counterparts without MPs, indicating that cluster mass might not be the dominant parameter controlling the presence of MPs. The transition phase between the occurrence of clusters with or without MPs seems to occur at ages of 2--3 Gyr, indicating that cluster age may play a dominant role in the establishment of MPs.

Large scale imaging surveys will increase the number of galaxy-scale strong lensing candidates by maybe three orders of magnitudes beyond the number known today. Finding these rare objects will require picking them out of at least tens of millions of images and deriving scientific results from them will require quantifying the efficiency and bias of any search method. To achieve these objectives automated methods must be developed. Because gravitational lenses are rare objects reducing false positives will be particularly important. We present a description and results of an open gravitational lens finding challenge. Participants were asked to classify 100,000 candidate objects as to whether they were gravitational lenses or not with the goal of developing better automated methods for finding lenses in large data sets. A variety of methods were used including visual inspection, arc and ring finders, support vector machines (SVM) and convolutional neural networks (CNN). We find that many of the methods will be easily fast enough to analyse the anticipated data flow. In test data, several methods are able to identify upwards of half the lenses after applying some thresholds on the lens characteristics such as lensed image brightness, size or contrast with the lens galaxy without making a single false-positive identification. This is significantly better than direct inspection by humans was able to do. (abridged)

We investigate the clustering properties and close neighbour counts for galaxies with different types of bulges and stellar masses. We select samples of "classical" and "pseudo" bulges, as well as "bulge-less" disk galaxies, based on the bulge/disk decomposition catalog of SDSS galaxies provided by Simard et al. (2011). For a given galaxy sample we estimate: the projected two-point cross-correlation function with respect to a spectroscopic reference sample, w_p(r_p), and the average background-subtracted neighbour count within a projected separation using a photometric reference sample, N_neighbour(<r_p). We compare the results with the measurements of control samples matched in color, concentration and redshift. We find that, when limited to a certain stellar mass range and matched in color and concentration, all the samples present similar clustering amplitudes and neighbour counts on scales above ~0.1h^{-1}Mpc. This indicates that neither the presence of a central bulge, nor the bulge type is related to intermediate-to-large scale environments. On smaller scales, in contrast, pseudo-bulge and pure-disk galaxies similarly show strong excess in close neighbour count when compared to control galaxies, at all masses probed. For classical bulges, small-scale excess is also observed but only for M_stars < 10^{10} M_sun; at higher masses, their neighbour counts are similar to that of control galaxies at all scales. These results imply strong connections between galactic bulges and galaxy-galaxy interactions in the local Universe, although it is unclear how they are physically linked in the current theory of galaxy formation.

We analyze the stellar initial mass functions (IMF) of a large sample of early type galaxies (ETGs) provided by MaNGA. The large number of IFU spectra of individual galaxies provide high signal-to-noise composite spectra that are essential for constraining IMF and to investigate possible radial gradients of the IMF within individual galaxies. The large sample of ETGs also make it possible to study how the IMF shape depends on various properties of galaxies. We adopt a novel approach to study IMF variations in ETGs, use Bayesian inferences based on full spectrum fitting. The Bayesian method provides a statistically rigorous way to explore potential degeneracy in spectrum fitting, and to distinguish different IMF models with Bayesian evidence. We find that the IMF slope depends systematically on galaxy velocity dispersion, in that galaxies of higher velocity dispersion prefer a more bottom-heavy IMF, but the dependence is almost entirely due to the change of metallicity, $Z$, with velocity dispersion. The IMF shape also depends on stellar age, $A$, but the dependence is completely degenerate with that on metallicity through a combination $AZ^{-1.42}$. Using independent age and metallicity estimates we find that the IMF variation is produced by metallicity instead of age. The IMF near the centers of massive ETGs appears more bottom-heavy than that in the outer parts, while a weak opposite trend is seen for low-mass ETGs. Uncertainties produced by star formation history, dust extinction, $\alpha$-element abundance enhancement and noise in the spectra are tested.

We construct a sample of 70 clusters using data from XMM-Newton and Planck to investigate the $Y_{\rm SZ,Planck}-Y_{\rm SZ, XMM}$ scaling relation and the cool-core influences on the relation. $Y_{\rm SZ,XMM}$ is calculated by accurate de-projected temperature and electron number density profiles derived from XMM-Newton. $Y_{\rm SZ,Planck}$ is the latest Planck data restricted to our precise X-ray size $\theta_{\rm 500}$. To study the cool-core influences on$Y_{\rm SZ,Planck}-Y_{\rm SZ, XMM}$ scaling relation, we apply two criteria, limits of central cooling time and classic mass deposition rate, to distinguish cool-core clusters (CCCs) from non-cool-core clusters (NCCCs). We also use $Y_{\rm SZ,Planck}$ from other papers, which are derived from different methods, to confirm our results. The intercept and slope of the$Y_{\rm SZ,Planck}-Y_{\rm SZ, XMM}$ scaling relation are $A=-0.86 \pm 0.30$, $B=0.83 \pm 0.06$. The intrinsic scatter is $\sigma_{\rm ins}= 0.14 \pm 0.03$. The ratio of \mbox{$Y_{\rm SZ,Planck}/Y_{\rm SZ, XMM}$} is $1.03 \pm 0.05$, which is perfectly agreed with unity. Discrepancies of $Y_{\rm SZ,Planck}-Y_{\rm SZ, XMM}$ scaling relation between CCCs and NCCCs are found in observation. They are independent of cool-core classification criteria and $Y_{\rm SZ,Planck}$ calculation methods, although discrepancies are more significant under the classification criteria of classic mass deposition rate. The intrinsic scatter of CCCs (0.04) is quite small compared to that of NCCCs (0.27). The ratio of $Y_{\rm SZ,Planck}/Y_{\rm SZ, XMM}$ for CCCs is $0.89 \pm 0.05$, suggesting that CCCs' $Y_{\rm SZ,XMM}$ may overestimate SZ signal. By contrast, the ratio of $Y_{\rm SZ,Planck}/Y_{\rm SZ, XMM}$ for NCCCs is $1.14 \pm 0.12$, which indicates that NCCCs' $Y_{\rm SZ,XMM}$ may underestimate SZ signal.

We use a sample of ~3000 galaxies from the MaNGA MPL-7 internal data release to study the alpha abundance distribution within low-redshift galaxies. We use the Lick index ratio Mgb/<Fe> as an alpha abundance indicator to study relationships between the alpha abundance distribution and galaxy properties such as effective stellar velocity dispersion within 0.3 effective radii (sigma_*), galaxy environment, and dark matter halo formation time (z_f). We find that (1) all galaxies show a tight correlation between Mgb/<Fe> and sigma_*; (2) `old' (H_beta < 3) low-sigma_* galaxies in high local density environment and inner regions within galaxy groups are enhanced in Mgb/<Fe>, while `young' (H_beta>3) galaxies and high-mass galaxies show no or less environmental dependence; (3) `old' galaxies with high-z_f show enhanced Mgb/<Fe> over low- and medium-z_f; (4) Mgb/<Fe> gradients are close to zero and show dependence on sigma_* but no obvious dependence on the environment or z_f. Our study indicates that stellar velocity dispersion or galaxy mass is the main parameter driving the Mgb/<Fe> enhancement, although environments appear to have modest effects, particularly for low- and medium-mass galaxies.

State of the art radial-velocity (RV) exoplanet searches are currently limited by RV signals arising from stellar magnetic activity. We analyze solar observations acquired over a 3-year period during the decline of Carrington Cycle 24 to test models of RV variation of Sun-like stars. A purpose-built solar telescope at the High Accuracy Radial velocity Planet Searcher for the Northern hemisphere (HARPS-N) provides disk-integrated solar spectra, from which we extract RVs and $\log{R'_{\rm HK}}$. The Solar Dynamics Observatory (SDO) provides disk-resolved images of magnetic activity. The Solar Radiation and Climate Experiment (SORCE) provides near-continuous solar photometry, analogous to a Kepler light curve. We verify that the SORCE photometry and HARPS-N $\log{R'_{\rm HK}}$ correlate strongly with the SDO-derived magnetic filling factor, while the HARPS-N RV variations do not. To explain this discrepancy, we test existing models of RV variations. We estimate the contributions of the suppression of convective blueshift and the rotational imbalance due to brightness inhomogeneities to the observed HARPS-N RVs. We investigate the time variation of these contributions over several rotation periods, and how these contributions depend on the area of active regions. We find that magnetic active regions smaller than $60 \ \rm Mm^2$ do not significantly suppress convective blueshift. Our area-dependent model reduces the amplitude of activity-induced RV variations by a factor of two. The present study highlights the need to identify a proxy that correlates specifically with large, bright magnetic regions on the surfaces of exoplanet-hosting stars.

We study the spatially resolved star formation history and molecular gas distribution of 58 nearby galaxies, using integral field spectroscopy from the CALIFA survey and CO $J=1\rightarrow 0$ intensity mapping from the CARMA EDGE survey. We use the 4000 \AA\ break (D$_n$4000), the equivalent width of the H$\delta$ absorption line (EW H$\delta_A$), and the equivalent width of the H$\alpha$ emission line (EW H$\alpha$) to measure the recent star formation history (SFH) of these galaxies. We measure radial profiles of the three SFH indicators and molecular gas mass surface density, from which we measure the level of centrally enhanced star formation and the molecular gas concentration. When we separate our galaxies into categories of barred (17 galaxies), unbarred (24 galaxies), and merging/paired (17 galaxies) we find that the galaxies which have centrally-enhanced star formation (19/58) are either barred (13/19) or in mergers/pairs (6/19) with relatively high molecular gas concentrations. A comparison between our barred galaxies and a snapshot of a hydrodynamic $N$-body simulation of a barred galaxy shows that the current theory of bar formation and evolution can qualitatively reproduce the main features of the observed galaxies in our sample, including both the sharp decrease of stellar age in the galactic center and the gradual decrease of age with increasing distance from center. These findings provide substantial evidence for a picture in which cold gas is transported inward by a bar or tidal interaction, which leads to the growth and rejuvenation of star formation in the central region.

Using data collected by the Askaryan Radio Array (ARA) experiment at the South Pole, we have used long-baseline propagation of radio-frequency signals to extract information on the radio-frequency index-of-refraction in South Polar ice. Owing to the increasing ice density over the upper 150--200 meters, rays are observed along two, nearly parallel paths, one of which is direct and a second which refracts through an inflection point, with differences in both arrival time and arrival angle that can be used to constrain the neutrino properties. We also observe indications, for the first time, of radio-frequency ice birefringence for signals propagating along predominantly horizontal trajectories, corresponding to an asymmetry of order 0.1% between the ordinary and extra-ordinary paths, numerically compatible with previous measurements of birefringent asymmetries for vertically-propagating radio-frequency signals at South Pole. Taken together, these effects offer the possibility of redundantly measuring the range from receiver to a neutrino interaction in Antarctic ice, if receiver antennas are deployed at shallow (25 m<z<100 m) depths. Such range information is essential in determining both the neutrino energy, as well as the incident neutrino direction.

A new non-hydrostatic and cloud-resolving atmospheric model is developed for studying moist convection and cloud formation in planetary atmospheres. It is built on top of the Athena++ framework, utilizing its static/adaptive mesh-refinement, parallelization, curvilinear geometry, and dynamic task scheduling. We extend the original hydrodynamic solver to vapors, clouds, and precipitation. Microphysics is formulated generically so that it can be applied to both Earth and Jovian planets. We implemented the Low Mach number Approximate Riemann Solver (LMARS) for simulating low speed atmospheric flows in addition to the usual Roe and HLLC Riemann solvers. Coupled with a fifth-order Weighted Essentially Nonoscillatory (WENO) subgrid-reconstruction method, the sharpness of critical fields such as clouds is well-preserved, and no extra hyperviscosity or spatial filter is needed to stabilize the model. Unlike many atmospheric models, total energy is used as the prognostic variable of the thermodynamic equation. One significant advantage of using total energy as a prognostic variable is that the entropy production due to irreversible mixing process can be properly captured. The model is designed to provide a unified framework for exploring planetary atmospheres across various conditions, both terrestrial and Jovian. First, a series of standard numerical tests for Earth's atmosphere is carried out to demonstrate the performance and robustness of the new model. Second, simulation of an idealized Jovian atmosphere in radiative-convective equilibrium shows that 1) the temperature gradient is superadiabatic near the water condensation level because of the changing of the mean molecular weight, and 2) the mean profile of ammonia gas shows a depletion in the subcloud layer down to nearly 10 bars. Relevance to the recent Juno observations is discussed.

Blue straggler stars have been proposeed as powerful indicators to measure the dynamical state of Galactic globular clusters. Here we examine for the first time if this framework of blue straggler stars as dynamical clocks, which was specifically developed for Galactic globular clusters, may also hold for younger globular clusters in the Large Magellanic Cloud. Using Hubble Space Telescope observations, we study seven Large Magellanic Cloud star clusters with ages between ~700 Myr and 7 Gyr. We confirm that our sample clusters are consistent with the empirical correlation previously derived for Galactic globular clusters, within a modest tolerance. We suggest that to further examine if blue straggler stars can measure the dynamical state of their host clusters over an extended range of dynamical ages, more studies of dynamically older Magellanic Cloud clusters are required. We discuss the physical implications of our results in terms of their central, dimensionless King potential, as well as the initial retention fraction of black holes.

Li et al. (2018) (hereafter L18) detected two apparently distinct populations of blue straggler stars (BSSs) in the young globular cluster NGC 2173, a similar feature as observed in numerous Galactic globular clusters (GCs). Recently, Dalessandro et al. (2018) (D18) compared the observed color--magnitude diagrams (CMDs) of both the cluster and a nearby reference field (although observed with the ACS/WFC instrument). They conclude that the bifurcated pattern of BSSs in NGC 2173 observed by L18 is a field contamination artefact. In this note, we explore the central concentration properties of the removed `field stars' identified by D18. Our purpose is to examine if these `field stars' are spatially homogeneously distributed. Employing a Monte Carlo-based approach, we have carefully studied the probability that any such central concentration may be caused by small number statistics. We find that, in most cases of, the `field stars' removed by D18 exhibit a clear central concentration, which cannot be explained on the basis of small number statistics alone. Therefore, we suggest that D18 may well have overestimated the field contamination level, implying that the bifurcated BSS pattern in NGC 2173 cannot, in fact, be explained by field contamination.

Twenty years have passed since first light for the Sloan Digital Sky Survey (SDSS). Here, we release data taken by the fourth phase of SDSS (SDSS-IV) across its first three years of operation (July 2014-July 2017). This is the third data release for SDSS-IV, and the fifteenth from SDSS (Data Release Fifteen; DR15). New data come from MaNGA - we release 4824 datacubes, as well as the first stellar spectra in the MaNGA Stellar Library (MaStar), the first set of survey-supported analysis products (e.g. stellar and gas kinematics, emission line, and other maps) from the MaNGA Data Analysis Pipeline (DAP), and a new data visualisation and access tool we call "Marvin". The next data release, DR16, will include new data from both APOGEE-2 and eBOSS; those surveys release no new data here, but we document updates and corrections to their data processing pipelines. The release is cumulative; it also includes the most recent reductions and calibrations of all data taken by SDSS since first light. In this paper we describe the location and format of the data and tools and cite technical references describing how it was obtained and processed. The SDSS website (www.sdss.org) has also been updated, providing links to data downloads, tutorials and examples of data use. While SDSS-IV will continue to collect astronomical data until 2020, and will be followed by SDSS-V (2020-2025), we end this paper by describing plans to ensure the sustainability of the SDSS data archive for many years beyond the collection of data.

Although the formation of blue straggler stars (BSSs) is routinely attributed to stellar interactions in binary systems, the relative importance of the direct collision and slow(er) stellar coalescence formation channels is still poorly understood. We selected a sample of 24 Magellanic Cloud star clusters for which multi-passband Hubble Space Telescope images are available to address this outstanding question. We compiled a BSS database, containing both traditional and evolved BSSs. We found a robust correlation between the number of BSSs in a cluster's core and its core mass, $N_{\rm BSS, core} \propto M_{\rm core}^{0.51 \pm 0.07}$, which supports the notion that BSS formation is linked to a population's binary fraction. At low stellar collision rates, the mass-normalised number of BSSs does not appear to depend on the collision rate, which implies that the coalescence-driven BSS formation channel dominates. Comparison with simulations suggests that stellar collisions contribute less than 20% to the total number of BSSs formed.

We investigate the clustering and dark halo properties for the narrow-line active galactic nuclei (AGN) in the SDSS, particularly examining the joint dependence on galaxy mass and color. AGN in galaxies with blue colors or massive red galaxies with M*>~10^{10.5}Msun are found to show almost identical clustering amplitudes at all scales to control galaxies of the same mass, color and structural parameters. This suggests AGN activity in blue galaxies or massive red galaxies is regulated by internal processes, with no correlation with environment. The antibias of AGN at scales between ~100kpc and a few Mpc, as found in Li et al. (2006) for the AGN as a whole, is observed only for the AGN hosted by galaxies with red colors and relatively low masses <10^{10.5}Msun. A simple halo model in which AGN are preferentially found at dark halo centers can reproduce the observational results, but requiring a mass-dependent central fraction which is a factor of ~4 higher than the fraction estimated from the SDSS group catalogue. The same group catalogue reveals that the host groups of AGN in red satellites tend to have lower halo masses than control galaxies, while the host groups of AGN in red centrals tend to form earlier, as indicated by a larger stellar mass gap between the two most massive galaxies in the groups. Our result implies that the mass assembly history of dark halos may play an additional role in the AGN activity in low-mass red galaxies.

A new class of quiescent galaxies harboring possible AGN-driven winds has been discovered using spatially resolved optical spectroscopy from the ongoing SDSS-IV MaNGA survey. These galaxies, termed "red geysers", constitute $5-10\%$ of the local quiescent population and are characterized by narrow bisymmetric patterns in ionized gas emission features. Cheung et al. argued that these galaxies host large-scale AGN-driven winds that may play a role in suppressing star formation at late times. In this work, we test the hypothesis that AGN activity is ultimately responsible for the red geyser phenomenon. We compare the nuclear radio activity of the red geysers to a matched control sample with similar stellar mass, redshift, rest frame $NUV-r$ color, axis ratio and presence of ionized gas. We have used the 1.4 GHz radio continuum data from VLA FIRST survey to stack the radio flux from the red geyser and control samples. In addition to a 3 times higher FIRST detection rate, we find that red geysers have a 5$\sigma$ higher level of average radio flux than control galaxies. After restricting to rest-frame $NUV-r$ color $>$ 5 and checking mid-IR WISE photometry, we rule out star formation contamination and conclude that red geysers are associated with more active AGN. Red geysers and a possibly-related class with disturbed H$\alpha$ emission account for 40\% of all radio-detected red galaxies with $\rm log~(M_\star/M_\odot) < 11$. Our results support a picture in which episodic AGN activity drives large-scale-relatively weak ionized winds that may provide a feedback mechanism for many early-type galaxies.

Source imaging of solar radio bursts can be used to track energetic electrons and associated magnetic structures. Here we present a combined analysis of data at different wavelengths for an eruption associated with a moving type-IV (t-IVm) radio burst. In the inner corona, the sources are correlated with a hot and twisted eruptive EUV structure, while in the outer corona the sources are associated with the top front of the bright core of a white light coronal mass ejection (CME). This reveals the potential of using t-IVm imaging data to continuously track the CME by lighting up the specific component containing radio-emitting electrons. It is found that the t-IVm burst presents a clear spatial dispersion with observing frequencies. The burst manifests broken power-law like spectra in brightness temperature, which is as high as $10^7$-$10^9$ K while the polarization level is in-general weak. In addition, the t-IVm burst starts during the declining phase of the flare with a duration as long as 2.5 hours. From the differential emission measure analysis of AIA data, the density of the T-IVm source is likely at the level of 10$^8$ cm$^{-3}$ at the start of the burst, and the temperature may reach up to several MK. These observations do not favor gyro-synchrotron to be the radiation mechanism, yet in line with a coherent plasma emission excited by energetic electrons trapped within the source. Further studies are demanded to elucidate the emission mechanism and explore the full diagnostic potential of t-IVm bursts.

The EDGES Collaboration has recently announced the detection of the 21-cm spectrum with an absorption profile centred at $78$ megahertz, of which the depth is deeper than that expected by the standard cosmological paradigm. To enrich the heating process of baryons due to scattering with dark matter during dark ages, we in this Letter explore the possibility of extra heat transfer between dark sector compositions and their observational signatures on the 21-cm cosmological spectrum. By parameterizing interaction models of the dark Universe, we find that the observational constraint on the parameter space of dark matter can be slightly relaxed but the discrepancy with the commonly predicted parameter space of weakly interacting massive particles remains. Our analyses also reveal that the interaction between dark compositions may leave observational signatures on the 21-cm spectrum during dark ages and thus would become detectable in the forthcoming 21-cm cosmology.

Ultra-high energy neutrinos are detectable through impulsive radio signals generated through interactions in dense media, such as ice. Subsurface in-ice radio arrays are a promising way to advance the observation and measurement of astrophysical high-energy neutrinos with energies above those discovered by the IceCube detector ($\geq$1 PeV) as well as cosmogenic neutrinos created in the GZK process ($\geq$100 PeV). Here we describe the $\textit{NuPhase}$ detector, which is a compact receiving array of low-gain antennas deployed 185 m deep in glacial ice near the South Pole. Signals from the antennas are digitized and coherently summed into multiple beams to form a low-threshold interferometric phased array trigger for radio impulses. The NuPhase detector was installed at an Askaryan Radio Array (ARA) station during the 2017/18 Austral summer season. $\textit{In situ}$ measurements with an impulsive, point-source calibration instrument show a 50% trigger efficiency on impulses with voltage signal-to-noise ratios (SNR) of $\le$2.0, a factor of $\sim$1.8 improvement in SNR over the standard ARA combinatoric trigger. Hardware-level simulations, validated with $\textit{in situ}$ measurements, predict a trigger threshold of an SNR as low as 1.6 for neutrino interactions that are in the far field of the array. With the already-achieved NuPhase trigger performance included in ARASim, a detector simulation for the ARA experiment, we find the trigger-level effective detector volume is increased by a factor of 1.8 at neutrino energies between 10 and 100 PeV compared to the currently used ARA combinatoric trigger. We also discuss an achievable near term path toward lowering the trigger threshold further to an SNR of 1.0, which would increase the effective single-station volume by more than a factor of 3 in the same range of neutrino energies.

We develop the effective field theory approach to torsional modified gravities, a formalism that allows for the systematic investigation of the background and perturbation levels separately. Starting from the usual effective field theory approach to curvature-based gravity, we suitably generalize it at the background level by including terms of the contracted torsion tensor, and at the perturbation level by including pure torsion perturbative terms and mixed perturbative terms of torsion and curvature. Having constructed the effective field theory action of general torsional modified gravity, amongst others we focus on $f(T)$ gravity and we perform a cosmological application. We investigate the scalar perturbations up to second order, and we derive the expressions of the Newtonian constant and the post Newtonian parameter $\gamma$. Finally, we apply this procedure to two specific and viable $f(T)$ models, namely the power-law and the exponential ones, introducing a new parameter that quantifies the deviation from general relativity and depends on the model parameters. Since this parameter can be expressed in terms of the scalar perturbation mode, a precise measurement of its evolution could be used as an alternative way to impose constraints on $f(T)$ gravity and break possible degeneracies between different $f(T)$ models.

JINGLE is a new JCMT legacy survey designed to systematically study the cold interstellar medium of galaxies in the local Universe. As part of the survey we perform 850um continuum measurements with SCUBA-2 for a representative sample of 193 Herschel-selected galaxies with M*>10^9Msun, as well as integrated CO(2-1) line fluxes with RxA3m for a subset of 90 of these galaxies. The sample is selected from fields covered by the Herschel-ATLAS survey that are also targeted by the MaNGA optical integral-field spectroscopic survey. The new JCMT observations combined with the multi-wavelength ancillary data will allow for the robust characterization of the properties of dust in the nearby Universe, and the benchmarking of scaling relations between dust, gas, and global galaxy properties. In this paper we give an overview of the survey objectives and details about the sample selection and JCMT observations, present a consistent 30 band UV-to-FIR photometric catalog with derived properties, and introduce the JINGLE Main Data Release (MDR). Science highlights include the non-linearity of the relation between 850um luminosity and CO line luminosity, and the serendipitous discovery of candidate z>6 galaxies.

Using the PMO-13.7m millimeter telescope at Delingha in China, we have conducted a large-scale simultaneous survey of $^{12}$CO, $^{13}$CO, and C$^{18}$O J=1-0 emission toward the Rosette molecular cloud (RMC) region with a sky coverage of 3.5 $\times$ 2.5 deg$^2$. The majority of the emission in the region comes from the RMC complex with velocities lying in the range from -2 km s$^{-1}$ to 20.5 km s$^{-1}$. Beyond this velocity range, 73 molecular clumps are identified with kinematic distances from 2.4 kpc to 11 kpc. Based on the spatial and velocity distribution, nine individual clouds, C1-C9, have been identified for the RMC complex. It appears that the C3 cloud is different from other clouds in the RMC complex in view of its characteristic velocity, excitation temperature, and velocity dispersion. Most of the young stellar clusters in the region are located in positions of both high column density and high excitation temperature. Seven new molecular filaments are discovered in the RMC complex. Evidence for cloud-cloud collision is found in the region of young stellar clusters REFL9 and PouF, showing that these young stellar clusters probably result from a cloud-cloud collision. The abundance ratios of $^{13}$CO to C$^{18}$O in the region have a mean value of 13.7 which is 2.5 times larger than the solar system value, showing that UV photons from the nearby OB clusters have strong influence on the chemistry of clouds in the RMC complex.

Halo merger trees are constructed from ELUCID, a constrained $N$-body simulation in the Sloan Digital Sky Survey (SDSS) volume. These merger trees are used to populate dark matter halos with galaxies according to an empirical model of galaxy formation. Mock catalogs in the SDSS sky coverage are constructed, which can be used to study the spatial distribution of galaxies in the low-$z$ Universe. These mock catalogs are used to quantify the cosmic variance in the galaxy stellar mass function (GSMF) measured from the SDSS survey. The GSMF estimated from the SDSS magnitude-limited sample can be affected significantly by the presence of the under-dense region at $z<0.03$, so that the low-mass end of the function can be underestimated significantly. Several existing methods designed to deal with the effects of the cosmic variance in the estimate of GSMF are tested, and none is found to be able to fully account for the cosmic variance. We propose a method based on the conditional stellar mass functions in dark matter halos, which can provide an unbiased estimate of the global GSMF. The application of the method to the SDSS data shows that the GSMF has a significant upturn at $M_*< 10^{9.5} h^{-1}{\rm M}_\odot$, which has been missed in many earlier measurements of the local GSMF.

The EDGES Collaboration has reported an anomalously strong 21cm absorption feature corresponding to the era of first star formation, which may indirectly betray the influence of dark matter during this epoch. We demonstrate that, by virtue of the ability to mediate cooling processes whilst in the condensed phase, a small amount of axion dark matter can explain these observations within the context of standard models of axions and axion-like-particles. The EDGES best-fit result favours an axion-like-particles mass in the (10, 450) meV range, which can be compressed for the QCD axion to (100, 450) meV in the absence of fine tuning. Future experiments and large scale surveys, particularly the International Axion Observatory (IAXO) and EUCLID, should have the capability to directly test this scenario.

We study the radial metallicity gradient $\Delta[M/H]/\Delta R_g$ as a function of [Mg/Fe] and $|Z|$ with the help of a guiding radius based on the Apache Point Observatory Galactic Evolution Experiment and Gaia and then analyze the radial migration effect on the radial metallicity gradient and metallicity-rotation gradient between the Galactic thin and thick disks. The derived trend of gradient $\Delta[M/H]/\Delta R_g$ versus [Mg/Fe] shows a transition at [Mg/Fe] $\sim 0.18$ dex, below which the gradient is negative and varies a little as [Mg/Fe] increases; however, it changes sharply in [Mg/Fe] ranges of 0.16-0.18, above which the gradient increases linearly with increasing [Mg/Fe], being a positive value at [Mg/Fe]$\gtrsim 0.22$ dex. These positive gradients in the high-[Mg/Fe] populations are found at $|Z| < 0.8$ kpc, and there are nearly no gradients toward higher $|Z|$. By comparing the metallicity distributions, the radial metallicity gradients $\Delta[M/H]/\Delta R$ and the metallicity-rotation gradients between the total sample and $|R-R_g|<2$ kpc subsample (or $|R-R_g|>2$ kpc subsample), we find that, for the thick disk, blurring flattens the gradient $\Delta[M/H]/\Delta R$ and favors metal-poor high-eccentricity stars. These stars are responsible for the measured positive metallicity-rotation gradient of the thick disk.

We present the X-ray timing results of the new black hole candidate (BHC) MAXI J1535-571 during its 2017 outburst from Hard X-ray Modulation Telescope (\emph{Insight}-HXMT) observations taken from 2017 September 6 to 23. Following the definitions given by \citet{Belloni2010}, we find that the source exhibits state transitions from Low/Hard state (LHS) to Hard Intermediate state (HIMS) and eventually to Soft Intermediate state (SIMS). Quasi-periodic oscillations (QPOs) are found in the intermediate states, which suggest different types of QPOs. With the large effective area of \emph{Insight}-HXMT at high energies, we are able to present the energy dependence of the QPO amplitude and centroid frequency up to 100 keV which is rarely explored by previous satellites. We also find that the phase lag at the type-C QPOs centroid frequency is negative (soft lags) and strongly correlated with the centroid frequency. By assuming a geometrical origin of type-C QPOs, the source is consistent with being a high inclination system.

Ly$\alpha$ Emitters (LAEs) may represent an important galaxy population in the low mass regime. We present our deep narrowband imaging surveys in the COSMOS and ECDF-S fields and study the properties of LAEs at $z=2.23\pm0.03$. The narrowband surveys conducted at Magellan II telescope allow us to obtain a sample of 452 LAEs reaching a $5\sigma$ limiting magnitude of $\sim26$ mag. Our Ly$\alpha$ luminosity functions extend to $10^{41.8}$ erg s$^{-1}$ with steep faint-end slope. Using multi-wavelength ancillary data, especially the deep Spitzer/IRAC 3.6$\mu$m and 4.5$\mu$m photometric data, we obtained reliable stellar mass estimates for 130 IRAC-detected LAEs, spanning a range of $8 < {\rm log}(M_\star/M_\odot)< 11.5$. For the remaining IRAC-undetected LAEs, the median-stacked spectral energy distribution yields a stellar mass of ${\rm log}(M_\star/M_\odot)=7.97^{+0.05}_{-0.07}$ and the rest-frame ultraviolet emission indicates a median star formation rate of ${\rm log} (SFR/M_\odot$ yr$^{-1})=-0.14\pm0.35$. There are six LAEs detected by the Spitzer/MIPS 24$\mu$m or even Herschel far-infrared observations. Taking into account the six MIR/FIR detected LAEs, our LAEs cover a wide range in the star formation rate (${\rm 1<SFR<2000}$ M$_\odot$ yr$^{-1}$). Although LAEs as a population are diverse in their stellar properties, they are mostly low-mass star-forming galaxies and follow the star formation main sequence relations or their extrapolations to the low-mass end, implying a normal star-forming nature of LAEs. The clustering analysis indicates that our LAEs reside in dark matter halos with ${\rm <\log(M_{h}/M_{\odot})> =10.8^{+0.56}_{-1.1}}$, suggesting that they are progenitors of local Large Magellanic Cloud-like galaxies.

This paper presents a photometric and spectroscopic study of the short-period binary star Cl*~Melotte~111~AV~1224. Measurements in the $B$, $V$, and $R$ passbands obtained during three observing runs between 2014 and 2017 and medium-resolution spectra secured in 2014, are analyzed together with public data from the SuperWASP and LAMOST projects. Our light curves show marked asymmetry with a variable O'Connell effect. The SuperWASP photometry is used to derive a mean binary period of 0.345225 days. The analysis of the $(O-C)$ diagram reveals that the orbital period is decreasing at a rate of $dP/dt = -3.87 \times 10^{-6}$ days yr$^{-1}$, which may be caused by mass transfer from the more-massive component to the less-massive one. The system is found to be a single-lined spectroscopic binary with a systemic velocity, $\gamma = 1 \pm 3$ Km s$^{-1}$, and a semi-amplitude, K$_{1}$ = 21 $\pm$ 5 Km s$^{-1}$. The spectral classification and the effective temperature of the primary component are estimated to be K0V $\pm$ 1 and $5200 \pm 150$ K, respectively. The photometric and spectroscopic solutions reveal that Cl*~Melotte~111~AV~1224 is a low-mass ratio ($q=m_{2}/m_{1} \sim 0.11$), low-inclination ($\sim ~ 38^{\circ}$) near-contact system. The masses, radii and luminosity for the primary and secondary are estimated to be $1.02 \pm 0.06\, M_\odot$, $1.23 \pm 0.05\, R_\odot $, $1.01 \pm 0.06\, L_\odot$ and $0.11 \pm 0.08\, M_\odot$, $0.45 \pm 0.05\, R_\odot$, $0.10 \pm 0.06\, L_\odot$, respectively. The marginal contact, together with the period decrease, suggests that this binary system may be at a key evolutionary stage, as predicted by the theory of thermal relaxation oscillations.

The Large Magellanic Cloud (LMC) has $\sim$60 confirmed supernova remnants (SNRs). Because of the known distance, 50 kpc, the SNRs' angular sizes can be converted to linear sizes, and their X-ray observations can be used to assess X-ray luminosities ($L_X$). We have critically examined the LMC SNRs' sizes reported in the literature to determine the most plausible sizes. These sizes and the $L_X$ determined from \emph{XMM-Newton} observations are used to investigate their relationship in order to explore the environmental and evolutionary effects on the X-ray properties of SNRs. We find: (1) Small LMC SNRs, a few to 10 pc in size, are all of Type Ia with $L_X>10^{36}$ ergs s$^{-1}$. The scarcity of small core-collapse (CC) SNRs is a result of CCSNe exploding in the low-density interiors of interstellar bubbles blown by their massive progenitors during their main sequence phase. (2) Medium-sized (10-30 pc) CC SNRs show bifurcation in $L_X$, with the X-ray-bright SNRs either in an environment associated with molecular clouds or containing pulsars and pulsar wind nebulae and the X-ray-faint SNRs being located in low-density interstellar environments. (3) Large (size$>$30 pc) SNRs show a trend of $L_X$ fading with size, although the scatter is large. The observed relationship between $L_X$ and sizes can help constrain models of SNR evolution.

We analyze {\sl Hubble Space Telescope} observations of nine Large Magellanic Cloud star clusters with ages of 1--2 Gyr to search for evolved counterparts of blue straggler stars. Near the red clump regions in the clusters' color--magnitude diagrams, we find branches of evolved stars that are much brighter than normal evolved stars. We examine the effects of photometric artifacts, differential reddening, and field contamination. We conclude that these bright evolved stars cannot be explained by any of these effects. Our statistical tests show that the contributions of photometric uncertainties and crowding effects, as well as that owing to differential reddening, to these bright evolved stars are insufficient to fully explain the presence of these bright evolved stars. Based on isochrone fitting, we also ruled out the possibility that these bright evolved stars could be reproduced by an internal chemical abundance spread. The spatial distributions of the bright evolved stars exhibit clear concentrations that cannot be explained by homogeneously distributed field stars. This is further confirmed based on Monte Carlo-based tests. By comparing our observations with stellar evolution models, we find that the masses of most of bright evolved stars do not exceed twice the average mass of normal evolved stars. We suggest that these bright evolved stars are, in fact, evolved blue straggler stars.

The Askaryan Radio Array (ARA) reports an observation of radio emission coincident with the "Valentine's Day" solar flare on Feb. 15$^{\rm{th}}$, 2011 in the prototype "Testbed" station. We find $\sim2000$ events that passed our neutrino search criteria during the 70 minute period of the flare, all of which reconstruct to the location of the sun. A signal analysis of the events reveals them to be consistent with that of bright thermal noise correlated across antennas. This is the first natural source of radio emission reported by ARA that is tightly reconstructable on an event-by-event basis. The observation is also the first for ARA to point radio from individual events to an extraterrestrial source on the sky. We comment on how the solar flares, coupled with improved systematic uncertainties in reconstruction algorithms, could aid in a mapping of any above-ice radio emission, such as that from cosmic-ray air showers, to astronomical locations on the sky.

The Tianlai Pathfinder is designed to demonstrate the feasibility of using a wide field of view radio interferometers to map the density of neutral hydrogen in the Universe after the Epoch of Reionizaton. This approach, called 21~cm intensity-mapping, promises an inexpensive means for surveying the large-scale structure of the cosmos. The Tianlai Pathfinder presently consists of an array of three, 15~m $\times$ 40~m cylinder telescopes and an array of sixteen, 6~m diameter dish antennas located in a radio-quiet part of western China. The two types of arrays were chosen to determine the advantages and disadvantages of each approach. The primary goal of the Pathfinder is to make 3D maps by surveying neutral hydrogen over large areas of the sky %$20,000 {\rm deg}^2$ in two different redshift ranges: first at $1.03 > z > 0.78$ ($700 - 800$~MHz) and later at $0.21 > z > 0.12$ ($1170 - 1270$~MHz). The most significant challenge to $21$~cm intensity-mapping is the removal of strong foreground radiation that dwarfs the cosmological signal. It requires exquisite knowledge of the instrumental response, i.e. calibration. In this paper, we provide an overview of the status of the Pathfinder and discuss the details of some of the analysis that we have carried out to measure the beam function of both arrays. We compare electromagnetic simulations of the arrays to measurements, discuss measurements of the gain and phase stability of the instrument, and provide a brief overview of the data processing pipeline.

We investigate the origin of the type III (up-bending) discs based on a sample of 1808 galaxies from SDSS and a sub-sample of 286 galaxies with HI data from ALFALFA. We examine how the type III fraction f3, the fraction of disc galaxies which host up-bending disc breaks, depends on other galactic properties. We confirm that f3 strongly depends on the stellar concentration of galaxies. We find that HI-rich galaxies with low spins tend to have significantly more type III disc breaks than control galaxies, which are matched in concentration and stellar mass. This trend is independent of the existence of strong bars or environment of the galaxies. This result is broadly consistent with predictions from theoretical simulations, and indicates in-situ star formation fueled by gas accretion to be an important mechanism that builds the outer discs of type III galaxies.

Sky survey is one of the most important motivations to improve the astrophysics development, especially when using new photometric bands. We are performing the SAGE (Stellar Abundance and Galactic Evolution) survey with a self-designed SAGE photometric system, which is composed of eight photometric bands. The project mainly aims to study the stellar atmospheric parameters of $\sim$0.5 billion stars in the $\sim12,000$ deg$^2$ of the northern sky, which mainly focuses on the Galactic sciences, as well as some extragalactic sciences. This work introduces the detailed data reduction process of the testing field NGC\,6791, including the data reduction of single-exposure image and stacking multi-exposure images, and properties of the final catalogue.

To investigate in more details of Stellar Abundance and Galactic Evolution (SAGE) and in a huge sample, we are performing a northern sky photometric survey named SAGES with the SAGE photometric system, which consists of 8 filters: Str\"omgren-$u$, SAGE-$v$, SDSS $g$, $r$, $i$, DDO-$51$, $H\alpha_{wide}$, and $H\alpha _{narrow}$, including three Sloan broadband filters, three intermediate-band filters and two narrow-band filters, and one newly-designed narrow-band filter. SAGES covers $\sim$12,000 square degrees of the northern sky with $\delta > -5 ^{\circ}$, excluding the Galactic disk ($|b|<10^{\circ}$) and the sky area of 12 hr $<$ R.A. $<$ 18\,hr. The photometric detection limit depth at signal-to-noise ratio $5\sigma$ can be as deep as $V\sim$20\,mag. The SAGES will produce a depth-uniformed photometric catalogue for $\sim$500 million stars with atmospheric parameters including effective temperature $T_{\rm eff}$, surface gravity log\,g, and metallicity [Fe/H], as well as interstellar extinction to each individual target. In this work, we will briefly introduce the SAGE photometric system, the SAGE survey, and a preliminary test field of the open cluster NGC\,6791 and around.

We have analyzed populations of blue straggler stars (BSSs) in 24 Magellanic Cloud star clusters using multi-passband Hubble Space Telescope images. We compiled a homogeneous BSS database, containing both traditional and evolved BSSs. We uncovered a sub-linear correlation between the number of BSSs in the cluster cores and the clusters' core masses, characterized by a power-law index of $0.51\pm 0.07$. For low stellar collision rates, the mass-normalized number of BSSs depends only weakly (or perhaps not at all) on the collision rate, implying that the binary-driven BSS formation channel dominates. Comparison with simulations suggests that stellar collisions contribute less than 20\% to the total number of BSSs formed. Further tests, including analysis of the BSS specific frequencies and their population numbers at larger cluster radii, suggest that binary interactions may be their main formation channel, hinting at an anti-correlation between a cluster's binary fraction and its core mass.

A High-performance Atmospheric Radiation Package (HARP) is developed for studying multiple-scattering planetary atmospheres. HARP is an open-source program written in C++ that utilizes high-level data structure and parallel-computing algorithms. It is generic in three aspects. First, the construction of the model atmospheric profile is generic. The program can either take in an atmospheric profile or construct an adiabatic thermal and compositional profile, taking into account the clouds and latent heat release due to condensation. Second, the calculation of opacity is generic, based on line-by-line molecular transitions and tabulated continuum data, along with a table of correlated-k opacity provided as an option to speed up the calculation of energy fluxes. Third, the selection of the solver for the radiative transfer equation is generic. The solver is not hardwired in the program. Instead, based on the purpose, a variety of radiative transfer solvers can be chosen to couple with the atmosphere model and the opacity model. We use the program to investigate the radiative heating and cooling rates of all four giant planets in the Solar System. Our Jupiter's result is consistent with previous publications. Saturn has a nearly perfect balance between the heating rate and cooling rate. Uranus has the least radiative fluxes because of the lack of CH4 and its photochemical products. Both Uranus and Neptune suffer from a severe energy deficit in their stratospheres. Possible ways to resolve this issue are discussed. Finally, we recalculate the radiative time constants of all four giant planet atmospheres and find that the traditional values from (Conrath BJ, Gierasch PJ, Leroy SS. Temperature and Circulation in the Stratosphere of the Outer Planets. Icar. 1990;83:255-81) are significantly overestimated.

We report on results of an observing campaign to support the Juno mission. At the beginning of 2016, using TEXES (Texas Echelon cross-dispersed Echelle Spectrograph), mounted on the NASA Infrared Telescope Facility (IRTF), we obtained data cubes of Jupiter in the 1930--1943 cm$^{-1}$ spectral ranges (around 5 $\mu$m), which probe the atmosphere in the 1--4 bar region, with a spectral resolution of $\approx$ 0.15 cm$^{-1}$ and an angular resolution of $\approx$ 1.4". This dataset is analysed by a code that combines a line-by-line radiative transfer model with a non-linear optimal estimation inversion method. The inversion retrieves the vertical abundance profiles of NH$_3$ - which is the main contributor at these wavelengths - with a maximum sensitivity at $\approx$ 1--3 bar, as well as the cloud transmittance. This retrieval is performed on over more than one thousand pixels of our data cubes, producing maps of the disk, where all the major belts are visible. We present our retrieved NH$_3$ abundance maps which can be compared with the distribution observed by Juno's MWR (Bolton et al., 2017; Li et al., 2017) in the 2 bar region and discuss their significance for the understanding of Jupiter's atmospheric dynamics. We are able to show important latitudinal variations - such as in the North Equatorial Belt (NEB), where the NH$_3$ abundance is observed to drop down to 60 ppmv at 2 bar - as well as longitudinal variability. In the zones, we find the NH$_3$ abundance to increase with depth, from 100 $\pm$ 15 ppmv at 1 bar to 500 $\pm$ 30 ppmv at 3 bar. We also display the cloud transmittance--NH$_3$ abundance relationship, and find different behaviour for the NEB, the other belts and the zones. Using a simple cloud model, we are able to fit this relationship, at least in the NEB, including either NH$_3$-ice or NH$_4$SH particles with sizes between 10 and 100 $\mu$m.

Coronal loops interconnecting two active regions, called as interconnecting loops (ILs), are prominent large-scale structures in the solar atmosphere. They carry a significant amount of magnetic flux, therefore are considered to be an important element of the solar dynamo process. Earlier observations show that eruptions of ILs are an important source of CMEs. It is generally believed that ILs are formed through magnetic reconnection in the high corona (>150-200"), and several scenarios have been proposed to explain their brightening in soft X-rays (SXRs). Yet, the detailed IL formation process has not been fully explored and the associated energy release in the corona still remains unresolved. Here we report the complete formation process of a set of ILs connecting two nearby active regions, with successive observations by STEREO-A on the far side of the Sun and SDO and Hinode on the Earth side. We conclude that ILs are formed by gradual reconnection high in the corona, in line with earlier postulations. In addition, we show evidence supporting that ILs become brightened in SXRs and EUVs through heating at or close to the reconnection site in the corona (i.e., through direct heating process of reconnection), a process that has been largely overlooked in earlier studies on ILs

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) has been in operation since July 2014. This paper describes the second data release from this phase, and the fourteenth from SDSS overall (making this, Data Release Fourteen or DR14). This release makes public data taken by SDSS-IV in its first two years of operation (July 2014-2016). Like all previous SDSS releases, DR14 is cumulative, including the most recent reductions and calibrations of all data taken by SDSS since the first phase began operations in 2000. New in DR14 is the first public release of data from the extended Baryon Oscillation Spectroscopic Survey (eBOSS); the first data from the second phase of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE-2), including stellar parameter estimates from an innovative data driven machine learning algorithm known as "The Cannon"; and almost twice as many data cubes from the Mapping Nearby Galaxies at APO (MaNGA) survey as were in the previous release (N = 2812 in total). This paper describes the location and format of the publicly available data from SDSS-IV surveys. We provide references to the important technical papers describing how these data have been taken (both targeting and observation details) and processed for scientific use. The SDSS website (www.sdss.org) has been updated for this release, and provides links to data downloads, as well as tutorials and examples of data use. SDSS-IV is planning to continue to collect astronomical data until 2020, and will be followed by SDSS-V.

In this paper, we use the equation of state based on a modification of 2+1 flavors Nambu-Jona- Lasinio (NJL) model to study the quark matter of hybrid stars. For comparison, we utilize five EOSs of the relativistic mean-field (RMF) model to describe the hadronic phase. With the 3- window crossover interpolation approach, we try to construct relatively soft hybrid EOSs but find the maximum masses of hybrid stars do not differ too much. The results are quite close to two solar mass, which is consistent with the mass constraint of PSR J0348+0432. Furthermore it is noteworthy that the heaviest stable stars have central densities higher than that of the deconfinement transition thus suggesting a pure quark core in the hybrid star.

Recent discoveries of bimodal main sequences (MSs) associated with young clusters (with ages $\lesssim 1$ Gyr) in the Magellanic Clouds have drawn a lot of attention. One of the prevailing formation scenarios attributes these split MSs to a bimodal distribution in stellar rotation rates, with most stars belonging to a rapidly rotating population. In this scenario, only a small fraction of stars populating a secondary blue sequence are slowly or non-rotating stars. Here, we focus on the blue MS in the young cluster NGC 1850. We compare the cumulative number fraction of the observed blue-MS stars to that of the high-mass-ratio binary systems at different radii. The cumulative distributions of both populations exhibit a clear anti-correlation, characterized by a highly significant Pearson coefficient of $-0.97$. Our observations are consistent with the possibility that blue-MS stars are low-mass-ratio binaries, and therefore their dynamical disruption is still ongoing. High-mass-ratio binaries, on the other hand, are more centrally concentrated.

Using the high-resolution observations obtained by the Hubble Space Telescope, we analyzed the blue straggler stars (BSSs) in the Large Magellanic Cloud cluster NGC 2213. We found that the radial distribution of BSSs is consistent with that of the normal giant stars in NGC 2213, showing no evidence of mass segregation. However, an analytic calculation carried out for these BSSs shows that they are already dynamically old, because the estimated half-mass relaxation time for these BSSs is significantly shorter than the isochronal age of the cluster. We also performed direct N-body simulations for a NGC 2213-like cluster to understand the dynamical processes that lead to this none-segregated radial distribution of BSSs. Our numerical simulation shows that the presence of black hole subsystems inside the cluster centre can significantly affect the dynamical evolution of BSSs. The combined effects of the delayed segregation, binary disruption and exchange interactions of BSS progenitor binaries may result in this none-segregated radial distribution of BSSs in NGC 2213.

The split main sequences (MSs) and extended MS turnoffs (eMSTOs) detected in a few young clusters have demonstrated that these stellar systems host multiple populations differing in a number of properties such as rotation and, possibly, age.We analyze Hubble Space Telescope photometry for thirteen clusters with ages between ~40 and ~1000 Myrs and of different masses. Our goal is to investigate for the first time the occurrence of multiple populations in a large sample of young clusters. We find that all the clusters exhibit the eMSTO phenomenon and that MS stars more massive than ~1.6 solar masses define a blue and red MS, with the latter hosting the majority of MS stars. The comparison between the observations and isochrones suggests that the blue MSs are made of slow-rotating stars, while the red MSs host stars with rotational velocities close to the breakup value. About half of the bright MS stars in the youngest clusters are H-alpha emitters. These Be stars populate the red MS and the reddest part of the eMSTO thus supporting the idea that the red MS is made of fast rotators. We conclude that the split MS and the eMSTO are a common feature of young clusters in both Magellanic Clouds. The phenomena of a split MS and an eMSTO occur for stars that are more massive than a specific threshold which is independent of the host-cluster mass. As a by-product, we report the serendipitous discovery of a young SMC cluster, GALFOR1.

We investigate radial gradients in the recent star formation history (SFH) using 1917 galaxies with $0.01<z<0.14$ and integral-field spectroscopy from the ongoing MaNGA survey. For each galaxy, we obtain two-dimensional maps and radial profiles for three spectroscopically-measured parameters that are sensitive to the recent SFH: D$_n$(4000) (the 4000\AA\ break), EW(H$\delta_A$) (equivalent width of the H$\delta$ absorption line), and EW(H$\alpha$) (equivalent width of the H$\alpha$ emission line). We find the majority of the spaxels in these galaxies are consistent with models with continuously declining star formation rate, indicating that starbursts occur rarely in local galaxies. We classify the galaxies into three classes: fully star-forming (SF), partly quenched (PQ) and totally quenched (TQ), according to the fraction of quenched area within 1.5 times the effective radius, $f_Q(1.5R_e)$. We find that galaxies less massive than $10^{10}$M$_{\odot}$ present at most weak radial gradients in all the diagnostic parameters. In contrast, galaxies with stellar mass above $10^{10}$M$_{\odot}$ present significant gradients in all the three diagnostic parameters if classified as SF or PQ, but show weak gradients in D$_n$(4000) and EW(H$\delta_A$) and no gradients in EW(H$\alpha$) if in the TQ class. This implies the existence of a critical stellar mass, at $\sim10^{10}$M$_{\odot}$, above which the star formation in a galaxy gets shutdown from the inside out. Galaxies tend to evolve synchronously from inner to outer regions before their mass reaches the critical value. We have further divided the sample at fixed mass by both bulge-to-total luminosity ratio and morphological type, finding that our conclusions hold regardless of these factors: it appears that the presence of a central dense object is not a driving parameter, but rather a byproduct of the star formation cessation process.

This paper presents the latest observations from the newly-built solar radio spectrograph at the \emph{Chashan Solar Observatory}. On July 18 2016, the spectrograph records a solar spike burst event, which has several episodes showing harmonic structures, with the second, third, and fourth harmonics. The lower harmonic radio spike emissions are observed later than the higher harmonic bands, and the temporal delay of the second (third) harmonic relative to the fourth harmonic is about 30\ --\ 40 (10) ms. Based on the electron cyclotron maser emission mechanism, we analyze possible causes of the temporal delay and further infer relevant coronal parameters, such as the magnetic field strength and the electron density at the radio source.

We perform full spectrum fitting stellar population analysis and Jeans Anisotropic modelling (JAM) of the stellar kinematics for about 2000 early-type galaxies (ETGs) and spiral galaxies from the MaNGA DR14 sample. Galaxies with different morphologies are found to be located on a remarkably tight mass plane which is close to the prediction of the virial theorem, extending previous results for ETGs. By examining an inclined projection (`the mass-size' plane), we find that spiral and early-type galaxies occupy different regions on the plane, and their stellar population properties (i.e. age, metallicity and stellar mass-to-light ratio) vary systematically along roughly the direction of velocity dispersion, which is a proxy for the bulge fraction. Galaxies with higher velocity dispersions have typically older ages, larger stellar mass-to-light ratios and are more metal rich, which indicates that galaxies increase their bulge fractions as their stellar populations age and become enriched chemically. The age and stellar mass-to-light ratio gradients for low-mass galaxies in our sample tend to be positive ($\rm centre<outer$), while the gradients for most massive galaxies are negative. The metallicity gradients show a clear peak around velocity dispersion $\log_{10} \sigma_{\rm e}\approx 2.0$, which corresponds to the critical mass $\sim 3\times 10^{10}M_{\odot}$ of the break in the mass-size relation. Spiral galaxies with large mass and size have the steepest gradients, while the most massive ETGs, especially above the critical mass $M_{\rm crit}\ge 2\times 10^{11} M_{\odot}$, where slow rotator ETGs start dominating, have much flatter gradients. This may be due to differences in their evolution histories, e.g. mergers.

Bifurcated patterns of blue straggler stars in their color--magnitude diagrams have atracted significant attention. This type of special (but rare) pattern of two distinct blue straggler sequences is commonly interpreted as evidence of cluster core-collapse-driven stellar collisions as an efficient formation mechanism. Here, we report the detection of a bifurcated blue straggler distribution in a young Large MagellanicCloud cluster, NGC 2173. Because of the cluster's low central stellar number density and its young age, dynamical analysis shows that stellar collisions alone cannot explain the observed blue straggler stars. Therefore, binary evolution is instead the most viable explanation of the origin of these blue straggler stars. However, the reason why binary evolution would render the color--magnitude distribution of blue straggler stars bifurcated remains unclear.

The New IRAM KID Arrays 2 (NIKA2) consortium has just finished installing and commissioning a millimetre camera on the IRAM 30 m telescope. It is a dual-band camera operating with three frequency multiplexed kilo-pixels arrays of Lumped Element Kinetic Inductance Detectors (LEKID) cooled at 150 mK, designed to observe the intensity and polarisation of the sky at 260 and 150 GHz (1.15 and 2 mm). NIKA2 is today an IRAM resident instrument for millimetre astronomy, such as Intra Cluster Medium from intermediate to distant clusters and so for the follow-up of Planck satellite detected clusters, high redshift sources and quasars, early stages of star formation and nearby galaxies emission. We present an overview of the instrument performance as it has been evaluated at the end of the commissioning phase.

With many large science equipment constructing and putting into use, astronomy has stepped into the big data era. The new method and infrastructure of big data processing has become a new requirement of many astronomers. Cloud computing, Map/Reduce, Hadoop, Spark, etc. many new technology has sprung up in recent years. Comparing to the high performance computing(HPC), Data is the center of these new technology. So, a new computing architecture infrastructure is necessary, which can be shared by both HPC and big data processing. Based on Astronomy Cloud project of Chinese Virtual Observatory (China-VO), we have made much efforts to optimize the designation of the hybrid computing platform. which include the hardware architecture, cluster management, Job and Resource scheduling.

Public observatory project is playing more and more important role in science popularization education and scientific research, and many amateur astronomers also have began to build their own observatories in remote areas. As a result of the limitation of technical condition and construction funds for amateur astronomers, their system often breaks down, and then a stable remote unattended operation system becomes very critical. Hardware connection and control is the basic and core part in observatory design. Here we propose a conception of engineering hardware design for public observatory operation as a bridge between observatory equipment and observation software. It can not only satisfy multiple observation mode requirement, but also save cost.

The combined observation of GW170817 and its electromagnetic counterpart GRB170817A reveals that gravitational waves propagate at the speed of light in high precision. We apply the effective field theory approach to investigate the experimental consequences for the theory of $f(T)$ gravity. We find that the speed of gravitational waves within $f(T)$ gravity is exactly equal to the light speed, and hence the constraints from GW170817 and GRB170817A are trivially satisfied. The results are verified through the standard analysis of cosmological perturbations. Nevertheless, by examining the dispersion relation and the frequency of cosmological gravitational waves, we observe a deviation from the results of General Relativity, quantified by a new parameter. Although its value is relatively small in viable $f(T)$ models, its possible future measurement in advancing gravitational-wave astronomy would be the smoking gun of testing this type of modified gravity.

We derived a new formula for calculating the moist adiabatic temperature profile of an atmosphere consisting of ideal gases with multiple condensing species. This expression unifies various formulas published in the literature and can be generalized to account for chemical reactions. Unlike previous methods, it converges to machine precision independent of mesh size. It accounts for any ratio of condensable vapors to dry gas, from zero to infinity, and for variable heat capacities as a function of temperature. Because the derivation is generic, the new formula is not only applicable to planetary atmosphere in the solar system, but also to hot Jupiters and brown dwarfs in which a variety of alkali metals, silicates and exotic materials condense. We demonstrate that even though the vapors are ideal gases, they interact in their effects on the moist adiabatic lapse rate. Finally, we apply the new thermodynamic model to study the effects of downdrafts on the distribution of minor constituents and thermal profile in the Galileo probe hotspot. We find that the Galileo Probe measurements can be interpreted as a strong downdraft that displaces an air parcel from 1 bar to the 4 bar level.

We have analyzed multi-passband photometric observations, obtained with the {\it Hubble Space Telescope}, of the massive ($1.8 \times 10^5 M_\odot$), intermediate-age (1.8 Gyr-old) Large Magellanic Cloud star cluster NGC 1783. The morphology of the cluster's red giant branch does not exhibit a clear broadening beyond its intrinsic width; the observed width is consistent with that owing to photometric uncertainties alone and independent of our photometric selection boundaries applied to obtain our sample of red-giant stars. The color dispersion of the cluster's red-giant stars around the best-fitting ridgeline is $0.062 \pm 0.009$ mag, which is equivalent to the width of $0.080 \pm 0.001$ mag derived from artificial simple stellar population tests, that is, tests based on single-age, single-metallicity stellar populations. NGC 1783 is comparably massive as other star clusters that show clear evidence of multiple stellar populations. After incorporating mass-loss recipes from its current age of 1.8 Gyr to an age of 6 Gyr, NGC 1783 is expected to remain as massive as some other clusters that host clear multiple populations at these intermediate ages. If we were to assume that mass is an important driver of multiple population formation, then NGC 1783 should have exhibited clear evidence of chemical abundance variations. However, our results support the absence of any chemical abundance variations in NGC 1783.

Further studies on the corner singularity of GS reconstruction are are compiled in this paper. It's focused on solution of the Data Completion (DC) problem with the Extended Hilbert Transform (EHT) over plane rectangular region. Optimal selections of the regular parameters in{\it Tikhonov} solution of corresponding DC problem are developed in this study. The 6-parameter regular solutions and the {\it Jacobian} matrix and one {\it Hessian} tensor to the regular parameters are derived in this work. A concise formula for EHT in the near field of corners, which shows property for EHT near the corner, are also provided. It serves as the additional constraints for our parameter optimization problem (OP). Third, a nonlinear convex function defined by the regular solution and the corner constraints is introduced and is taken as the object function for the OP of the 6 regular parameters on half-space $(\mathbf{p}>0)$. Given an initial guess of $\mathbf{p}_0$, the optimal parameters are solved from the OP through a well known constrained nonlinear optimization method. Last, the benchmark tests to the proposed solution approach are carried out, and detailed results from totally 9 different bench-cases are tabulated. In contrast to solutions with given regular parameters, our bench results demonstrate that an objective way for selection of the optimal $\mathbf{p}$ is successfully laid out here. Robustness and efficiency of the suggested new approach are also highlighted in this study.

Regions of disc galaxies with widespread star formation tend to be both gravitationally unstable and self-shielded against ionizing radiation, whereas extended outer discs with little or no star formation tend to be stable and unshielded on average. We explore what drives the transition between these two regimes, specifically whether discs first meet the conditions for self-shielding (parameterized by dust optical depth, $\tau$) or gravitational instability (parameterized by a modified version of Toomre's instability parameters, $Q_{\rm thermal}$, which quantifies the stability of a gas disc that is thermally supported at $T=10^4$ K). We first introduce a new metric formed by the product of these quantities, $Q_{\rm thermal}\tau$, which indicates whether the conditions for disk instability or self-shielding are easier to meet in a given region of a galaxy, and we discuss how $Q_{\rm thermal}\tau$ can be constrained even in the absence of direct gas information. We then analyse a sample of 13 galaxies with resolved gas measurements and find that on average galaxies will reach the threshold for disk instabilities ($Q_{\rm thermal}<1$) before reaching the threshold for self-shielding ($\tau>1$). Using integral field spectroscopic observations of a sample of 236 galaxies from the MaNGA survey, we find that the value of $Q_{\rm thermal}\tau$ in star-forming discs is consistent with similar behavior. These results support a scenario where disc fragmentation and collapse occurs before self-shielding, suggesting that gravitational instabilities are the primary condition for widespread star formation in galaxy discs. Our results support similar conclusions based on recent galaxy simulations.

Around the turn of the last century, star clusters of all kinds were considered "simple" stellar populations. Over the past decade, this situation has changed dramatically. At the same time, star clusters are among the brightest stellar population components and, as such, they are visible out to much greater distances than individual stars, even the brightest, so that understanding the intricacies of star cluster composition and their evolution is imperative for understanding stellar populations and the evolution of galaxies as a whole. In this review of where the field has moved to in recent years, we place particular emphasis on the properties and importance of binary systems, the effects of rapid stellar rotation, and the presence of multiple populations in Magellanic Cloud star clusters across the full age range. Our most recent results imply a reverse paradigm shift, back to the old simple stellar population picture for at least some intermediate-age (~1--3 Gyr-old) star clusters, opening up exciting avenues for future research efforts.

Dark matter(DM) is the only possible candidate which would be apart from the thermal equilibrium before Big Bang nucleosynthesis(BBN) in accordance with current DM searches. In this work, we report a generic scenario that primordial perturbation of dark matter(PPDM) can be, effectively, generated and encoded with primordial information of very early universe up to the reheating era. We present an analytical solution of the whole evolution of PPDM. A novel and strong constraint on the reheating process imposed by primordial gravitational wave(PGW) is obtained for the first time. It indicates the ratio of PGW to primordial curvature perturbation(PCP) is not only dependent on the slow-roll spectral index but also, strongly, on the decay process of inflaton at reheating. For very generic reheating process, our result provides a natural explanation of the paucity of PGW in current observations.

We study the role of cold gas in quenching star formation in the green valley by analysing ALMA $^{12}$CO (1-0) observations of three galaxies with resolved optical spectroscopy from the MaNGA survey. We present resolution-matched maps of the star formation rate and molecular gas mass. These data are used to calculate the star formation efficiency (SFE) and gas fraction ($f_{\rm~gas}$) for these galaxies separately in the central `bulge' regions and outer disks. We find that, for the two galaxies whose global specific star formation rate (sSFR) deviates most from the star formation main sequence, the gas fraction in the bulges is significantly lower than that in their disks, supporting an `inside-out' model of galaxy quenching. For the two galaxies where SFE can be reliably determined in the central regions, the bulges and disks share similar SFEs. This suggests that a decline in $f_{\rm~gas}$ is the main driver of lowered sSFR in bulges compared to disks in green valley galaxies. Within the disks, there exist common correlations between the sSFR and SFE and between sSFR and $f_{\rm~gas}$ on kpc scales -- the local SFE or $f_{\rm~gas}$ in the disks declines with local sSFR. Our results support a picture in which the sSFR in bulges is primarily controlled by $f_{\rm~gas}$, whereas both SFE and $f_{\rm~gas}$ play a role in lowering the sSFR in disks. A larger sample is required to confirm if the trend established in this work is representative of green valley as a whole.

Finding the electromagnetic (EM) counterpart of binary compact star merger, especially the binary neutron star (BNS) merger, is critically important for gravitational wave (GW) astronomy, cosmology and fundamental physics. On Aug. 17, 2017, Advanced LIGO and \textit{Fermi}/GBM independently triggered the first BNS merger, GW170817, and its high energy EM counterpart, GRB 170817A, respectively, resulting in a global observation campaign covering gamma-ray, X-ray, UV, optical, IR, radio as well as neutrinos. The High Energy X-ray telescope (HE) onboard \textit{Insight}-HXMT (Hard X-ray Modulation Telescope) is the unique high-energy gamma-ray telescope that monitored the entire GW localization area and especially the optical counterpart (SSS17a/AT2017gfo) with very large collection area ($\sim$1000 cm$^2$) and microsecond time resolution in 0.2-5 MeV. In addition, \textit{Insight}-HXMT quickly implemented a Target of Opportunity (ToO) observation to scan the GW localization area for potential X-ray emission from the GW source. Although it did not detect any significant high energy (0.2-5 MeV) radiation from GW170817, its observation helped to confirm the unexpected weak and soft nature of GRB 170817A. Meanwhile, \textit{Insight}-HXMT/HE provides one of the most stringent constraints (~10$^{-7}$ to 10$^{-6}$ erg/cm$^2$/s) for both GRB170817A and any other possible precursor or extended emissions in 0.2-5 MeV, which help us to better understand the properties of EM radiation from this BNS merger. Therefore the observation of \textit{Insight}-HXMT constitutes an important chapter in the full context of multi-wavelength and multi-messenger observation of this historical GW event.

It is now well established that globular clusters (GCs) exhibit star-to-star light-element abundance variations (known as multiple stellar populations, MPs). Such chemical anomalies have been found in (nearly) all the ancient GCs (more than 10 Gyr old) of our Galaxy and its close companions, but so far no model for the origin of MPs is able to reproduce all the relevant observations. To gain new insights into this phenomenon, we have undertaken a photometric Hubble Space Telescope survey to study clusters with masses comparable to that of old GCs, where MPs have been identified, but with significantly younger ages. Nine clusters in the Magellanic Clouds with ages between $\sim$ 1.5-11 Gyr have been targeted in this survey. We confirm the presence of multiple populations in all clusters older than 6 Gyr and we add NGC 1978 to the group of clusters for which MPs have been identified. With an age of $\sim$ 2 Gyr, NGC 1978 is the youngest cluster known to host chemical abundance spreads found to date. We do not detect evident star-to-star variations for slightly younger massive clusters ($\sim$ 1.7 Gyr), thus pointing towards an unexpected age dependence for the onset of multiple populations. This discovery suggests that the formation of MPs is not restricted to the early Universe and that GCs and young massive clusters share common formation and evolutionary processes.

The internal structures of pulsating white dwarfs can be explored only with asteroseismology. Time series photometric observations were made for the pulsating DA white dwarf (ZZ Ceti star) WD~0246+326 during 9 nights in 2014 with a bi-site observation campaign. Eleven frequencies were detected including 1 triplet, 2 doublets, and 4 single modes, which are identified as either $l=1$ or $l=2$ modes with the complementarity of frequencies present in the literature. From the multiplets, the rotation period of \astrobj{WD~0246+326} is derived as $3.78\pm 0.11$ days. The average period spacing of the l=1 modes $\Delta P=29.3\pm 0.2s$, implies that \astrobj{WD~0246+326} may be a massive ZZ Ceti star concerning the $\Delta P-M_*$ relationship for the DAVs. Preliminary analysis derives the stellar parameters of $M_*=0.98\pm0.01$~${\rm M_\odot}$ and $T_{\rm eff}=11700\pm100$~K by fitting the theoretical frequencies of the eigen modes to the observed ones.

Asteroseismology is a unique tool to explore the internal structure of stars through both observational and theoretical research. The internal structure of pulsating hydrogen shell white dwarfs (ZZ Ceti stars) detected by asteroseismology is regarded as the representative of all DA white dwarfs. Observations for KUV~08368+4026, which locates in the middle of the ZZ Ceti instability strip, have been carried out in 1999 and from 2009 to 2012 with either single-site runs or multisite campaigns. Time-series photometric data of about 300 hours were collected in total. Through data reduction and analysis, 30 frequencies were extracted, including four triplets, two doublets, one single mode and further signals. The independent modes are identified as either l=1 or l=2 modes. Hence, a rotation period of $5.52\pm 0.22$ days was deduced from the period spacing in the multiplets. Theoretical static models were built and a best fit model for KUV~08368+4026 was obtained with $0.692\pm0.002$ solar mass, $(2.92\pm0.02)\times 10^{-3}$ solar luminosity and the hydrogen mass fraction of $10^{-4}$ stellar mass.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) began observations in July 2014. It pursues three core programs: APOGEE-2, MaNGA, and eBOSS. In addition, eBOSS contains two major subprograms: TDSS and SPIDERS. This paper describes the first data release from SDSS-IV, Data Release 13 (DR13), which contains new data, reanalysis of existing data sets and, like all SDSS data releases, is inclusive of previously released data. DR13 makes publicly available 1390 spatially resolved integral field unit observations of nearby galaxies from MaNGA, the first data released from this survey. It includes new observations from eBOSS, completing SEQUELS. In addition to targeting galaxies and quasars, SEQUELS also targeted variability-selected objects from TDSS and X-ray selected objects from SPIDERS. DR13 includes new reductions of the SDSS-III BOSS data, improving the spectrophotometric calibration and redshift classification. DR13 releases new reductions of the APOGEE-1 data from SDSS-III, with abundances of elements not previously included and improved stellar parameters for dwarf stars and cooler stars. For the SDSS imaging data, DR13 provides new, more robust and precise photometric calibrations. Several value-added catalogs are being released in tandem with DR13, in particular target catalogs relevant for eBOSS, TDSS, and SPIDERS, and an updated red-clump catalog for APOGEE. This paper describes the location and format of the data now publicly available, as well as providing references to the important technical papers that describe the targeting, observing, and data reduction. The SDSS website, http://www.sdss.org, provides links to the data, tutorials and examples of data access, and extensive documentation of the reduction and analysis procedures. DR13 is the first of a scheduled set that will contain new data and analyses from the planned ~6-year operations of SDSS-IV.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) began observations in July 2014. It pursues three core programs: APOGEE-2, MaNGA, and eBOSS. In addition, eBOSS contains two major subprograms: TDSS and SPIDERS. This paper describes the first data release from SDSS-IV, Data Release 13 (DR13), which contains new data, reanalysis of existing data sets and, like all SDSS data releases, is inclusive of previously released data. DR13 makes publicly available 1390 spatially resolved integral field unit observations of nearby galaxies from MaNGA, the first data released from this survey. It includes new observations from eBOSS, completing SEQUELS. In addition to targeting galaxies and quasars, SEQUELS also targeted variability-selected objects from TDSS and X-ray selected objects from SPIDERS. DR13 includes new reductions of the SDSS-III BOSS data, improving the spectrophotometric calibration and redshift classification. DR13 releases new reductions of the APOGEE-1 data from SDSS-III, with abundances of elements not previously included and improved stellar parameters for dwarf stars and cooler stars. For the SDSS imaging data, DR13 provides new, more robust and precise photometric calibrations. Several value-added catalogs are being released in tandem with DR13, in particular target catalogs relevant for eBOSS, TDSS, and SPIDERS, and an updated red-clump catalog for APOGEE. This paper describes the location and format of the data now publicly available, as well as providing references to the important technical papers that describe the targeting, observing, and data reduction. The SDSS website, http://www.sdss.org, provides links to the data, tutorials and examples of data access, and extensive documentation of the reduction and analysis procedures. DR13 is the first of a scheduled set that will contain new data and analyses from the planned ~6-year operations of SDSS-IV.

We analyze observations of eight quiescent low-mass X-ray binaries in globular clusters and combine them to determine the neutron star mass-radius curve and the equation of state of dense matter. We determine the effect that several uncertainties may have on our results, including uncertainties in the distance, the atmosphere composition, the neutron star maximum mass, the neutron star mass distribution, the possible presence of a hotspot on the neutron star surface, and the prior choice for the equation of state of dense matter. We find that the radius of a 1.4 solar mass neutron star is most likely from 10 to 14 km and that tighter constraints are only possible with stronger assumptions about the nature of the neutron stars, the systematics of the observations, or the nature of dense matter. Strong phase transitions are preferred over other models and interpretations of the data with a Bayes factor of 8 or more, and in this case, the radius is likely smaller than 12 km. However, radii larger than 12 km are preferred if the neutron stars have uneven temperature distributions.

The spectroscopic sky survey carried out by the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) provides the largest stellar spectra library in the world until now. A large number of new DA white dwarfs had been identified based on the LAMOST spectra. The effective temperature ($T_{\rm eff}$) and surface gravity ($\log g$) of most DA white dwarfs were determined and published in the catalogs, e.g. Zhao et al. (2013), Rebassa-Mansergas et al. (2015), Gentile Fusillo et al. (2015) and Guo et al. (2015). We selected ZZ Ceti candidates from the published catalogs by considering whether their $T_{\rm eff}$ are situated in the ZZ Ceti instability strip. The follow-up time-series photometric observations for the candidates were performed in 2015 and 2016. Four stars: LAMOST J004628.31+343319.90, LAMOST J062159.49+252335.9, LAMOST J010302.46+433756.2 and LAMOST J013033.90+273757.9 are finally confirmed to be new ZZ Ceti stars. They show dominant peaks with amplitudes rising above the 99.9% confidence level in the amplitude spectra. As LAMOST J004628.31+343319.90 has an estimated mass of $\sim$ 0.40 $M_{\odot}$ and LAMOST J013033.90+273757.9 has a mass of $\sim$ 0.45 $M_{\odot}$ derived from their $\log g$ values, these two stars are inferred to be potential helium-core white dwarfs.

We study the dependence of galaxy clustering on atomic gas mass using a sample of $\sim$16,000 galaxies with redshift in the range of $0.0025<z<0.05$ and HI mass of $M_{\rm HI}>10^8M_{\odot}$, drawn from the 70% complete sample of the Arecibo Legacy Fast ALFA survey. We construct subsamples of galaxies with $M_{\rm HI}$ above different thresholds, and make volume-limited clustering measurements in terms of three statistics: the projected two-point correlation function, the projected cross-correlation function with respect to a reference sample selected from the Sloan Digital Sky Survey, and the redshift-space monopole moment. In contrast to previous studies, which found no/weak HI-mass dependence, we find both the clustering amplitude on scales above a few Mpc and the bias factors to increase significantly with increasing HI mass for subsamples with HI mass thresholds above $10^9M_{\odot}$. For HI mass thresholds below $10^9M_{\odot}$, while the measurements have large uncertainties caused by the limited survey volume and sample size, the inferred galaxy bias factors are systematically lower than the minimum halo bias factor from mass-selected halo samples. The simple halo model, in which galaxy content is only determined by halo mass, has difficulties in interpreting the clustering measurements of the HI-selected samples. We extend the simple model by including the halo formation time as an additional parameter. A model that puts HI-rich galaxies into halos that formed late can reproduce the clustering measurements reasonably well. We present the implications of our best-fitting model on the correlation of HI mass with halo mass and formation time, as well as the halo occupation distributions and HI mass functions for central and satellite galaxies. These results are compared with the predictions from semi-analytic galaxy formation models and hydrodynamic galaxy formation simulations.

Numerical simulations were carried out to study the origin of multiple stellar populations in the intermediate-age clusters NGC 411 and NGC 1806 in the Magellanic Clouds. We performed NBODY6++ simulations based on two different formation scenarios, an ad hoc formation model where second-generation (SG) stars are formed inside a cluster of first-generation (FG) stars using the gas accumulated from the external intergalactic medium and a minor merger model of unequal mass ($M_{SG}$/$M_{FG}$ ~5-10%) clusters with an age difference of a few hundred million years. We compared our results such as the radial profile of the SG-to-FG number ratio with observations on the assumption that the SG stars in the observations are composed of cluster members, and confirmed that both the ad hoc formation and merger scenarios reproduce the observed radial trend of the SG-to-FG number ratio which shows less centrally concentrated SG than FG stars. It is difficult to constrain the formation scenario for the multiple populations by only using the spatial distribution of the SG stars. SG stars originating from the merger scenario show a significant velocity anisotropy and rotational features compared to those from the ad hoc formation scenario. Thus, observations aimed at kinematic properties like velocity anisotropy or rotational velocities for SG stars should be obtained to better understand the formation of the multiple populations in these clusters. This is, however, beyond current instrumentation capabilities.

We describe the sample design for the SDSS-IV MaNGA survey and present the final properties of the main samples along with important considerations for using these samples for science. Our target selection criteria were developed while simultaneously optimizing the size distribution of the MaNGA integral field units (IFUs), the IFU allocation strategy, and the target density to produce a survey defined in terms of maximizing S/N, spatial resolution, and sample size. Our selection strategy makes use of redshift limits that only depend on i-band absolute magnitude ($M_i$), or, for a small subset of our sample, $M_i$ and color (NUV-i). Such a strategy ensures that all galaxies span the same range in angular size irrespective of luminosity and are therefore covered evenly by the adopted range of IFU sizes. We define three samples: the Primary and Secondary samples are selected to have a flat number density with respect to $M_i$ and are targeted to have spectroscopic coverage to 1.5 and 2.5 effective radii (Re), respectively. The Color-Enhanced supplement increases the number of galaxies in the low-density regions of color-magnitude space by extending the redshift limits of the Primary sample in the appropriate color bins. The samples cover the stellar mass range $5\times10^8 \leq M_* \leq 3\times10^{11} M_{\odot}$ and are sampled at median physical resolutions of 1.37 kpc and 2.5 kpc for the Primary and Secondary samples respectively. We provide weights that will statistically correct for our luminosity and color-dependent selection function and IFU allocation strategy, thus correcting the observed sample to a volume limited sample.

In this paper, we report the peculiar HI morphology of the cluster spiral galaxy NGC 6145, which has a 150 kpc HI filament on one side that is nearly parallel to its major axis. This filament is made up of several HI clouds and the diffuse HI gas between them, with no optical counterparts. We compare its HI distribution with other one-sided HI distributions in the literature, and find that the overall HI distribution is very different from the typical tidal and ram-pressure stripped HI shape, and its morphology is inconsistent with being a pure accretion event. Only about 30% of the total HI gas is anchored on the stellar disk, while most of HI gas forms the filament in the west. At a projected distance of 122 kpc, we find a massive elliptical companion (NGC 6146) with extended radio emission, whose axis points to an HI gap in NGC 6145. The velocity of the HI filament shows an overall light-of- sight motion of 80 to 180 km/s with respect to NGC 6145. Using the long-slit spectra of NGC 6145 along its major stellar axis, we find that some outer regions show enhanced star formation, while in contrast, almost no star formation activities are found in its center (less than 2 kpc). Pure accretion, tidal or ram-pressure stripping is not likely to produce the observed HI filament. An alternative explanation is the jet-stripping from NGC 6146, although direct evidence for a jet-cold gas interaction has not been found.

We describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratio in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spatially-resolved spectroscopy for thousands of nearby galaxies (median redshift of z = 0.03). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between redshifts z = 0.6 and 3.5 to constrain cosmology using baryon acoustic oscillations, redshift space distortions, and the shape of the power spectrum. Within eBOSS, we are conducting two major subprograms: the SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating X-ray AGN and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5-meter Sloan Foundation Telescope at Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5-meter du Pont Telescope at Las Campanas Observatory, with observations beginning in early 2017. Observations at both facilities are scheduled to continue through 2020. In keeping with previous SDSS policy, SDSS-IV provides regularly scheduled public data releases; the first one, Data Release 13, was made available in July 2016.