Recently solar flares exhibiting a double J-shaped ribbons in the lower solar atmosphere have been paid increasing attention in the context of extending the two-dimensional standard flare model to three dimensions, as motivated by the spatial correlation between photospheric current channels and flare ribbons. Here we study the electric currents through the photospheric area swept by flare ribbons (termed synthesized ribbon area or SRA), with a sample of 71 two-ribbon flares, of which 36 are J-shaped. Electric currents flowing through one ribbon are highly correlated with those through the other, therefore belonging to the same current system. The non-neutrality factor of this current system is independent of the flare magnitude, implying that both direct and return currents participate in flares. J-shaped flares are distinct from non-J-shaped flares in the following aspects: 1) Electric current densities within J-shaped SRA are significantly smaller than those within non-J-shaped SRA, but J-shaped SRA and its associated magnetic flux are also significantly larger. 2) Electric currents through SRA are positively correlated with the flare magnitude, but J-shaped flares show stronger correlation than non-J-shaped flares. 3) The majority (75\%) of J-shaped flares are eruptive, while the majority (86\%) of non-J-shaped flares are confined; accordingly, hosting active regions of J-shaped flares are more likely to be sigmoidal than non-J-shaped flares. Thus, J-shaped flares constitute a distinct subset of two-ribbon flares, probably the representative of eruptive ones. Further, we found that combining SRA and its associated magnetic flux has the potential to differentiate eruptive from confined flares.

Since Sep. 2018, LAMOST started the medium-resolution (R $\sim 7500$) spectral survey (MRS). We proposed the spectral survey of Galactic nebulae, including \hii\ regions, HH objects, supernova remnants, planetary nebulae and the special stars with MRS (LAMOST MRS-N). LAMOST MRS-N covers about 1700 square degrees of the northern Galactic plane within $40^\circ < l < 215^\circ$ and $-5^\circ < b < 5^\circ$. In this 5-year survey, we plan to observe about 500 thousand nebulae spectra. According to the commissioning observations, the nebulae spectra can provide precise radial velocity with uncertainty less than 1 km s$^{-1}$. These high-precision spectral data are of great significance to our understanding of star formation and evolution.

We present here the analysis of giant micropulses from the Vela pulsar. A total of 4187 giant micropulses with peak flux density $>$2.5 Jy were detected during almost 4 hours of observations carried out with the Yunnan 40-m radio telescope at 6800 MHz. Nine of the giant micropulses arrived approximately 3 to 4 ms earlier than the peak of average pulse profile, longer than that at lower frequencies. The remaining giant micropulses were clustered into three distributions which correspond to three main emission regions, including four occurring on the trailing edge of averaged profile.

The interstellar oxygen isotopic ratio of 18O/17O can reflect the relative amount of the secular enrichment by ejecta from high-mass versus intermediate-mass stars. Previous observations found a Galactic gradient of 18O/17O, i.e., low ratios in the Galactic center and large values in the Galactic disk, which supports the insideout formation scenario of our Galaxy. However, the observed objects are not many and, in particular, not so many at large galactocentric distances. Thus we started a systematic study on Galactic interstellar 18O/17O, through observations of C18O and C17O multi-transition lines toward a large sample of 286 sources (at least one order of magnitude larger than previous ones), from the Galactic center region to the far outer Galaxy (~22 kpc). Here we present our observations of J=1-0 lines of C18O and C17O, with the ARO12m and the IRAM 30m telescope. We detected successfully both C18O and C17O 1-0 lines for 34 sources among our IRAM30m sample of 50 targets and for 166 sources among our ARO12m sample of 260 targets. The C18O optical depth effect on our ratio results, evaluated by fitting results of C17O spectra with hyperfine components and our RADEX non-LTE model calculation for the strongest source, was found to be insignificant. Beam dilution does not seem to be a problem either, which was supported by the fact of no systematic variation between the isotopic ratio and the heliocentric distance, and consistent measured ratios from two telescopes for most of those detected sources. Our results, though there are still very few detections made for sources in the outer Galaxy, confirm the apparent 18O/17O gradient of 18O/17O = (0.10+-0.03)R_GC+(2.95+-0.30), with a Pearson's rank correlation coefficient R = 0.69. This is supported by the newest Galactic chemical evolution model including the impact of massive stellar rotators and novae.

QPOs are seen as peak features in the X-ray power spectral density of stellar mass black holes and neutron stars, of which frequencies can be measured with high precision. These QPO frequencies are believed to be related to fundamental frequencies of test particles, which are mainly determined by the background metric. We consider the metric introduced in Ghasemi-Nodehi and Bambi (Eur Phys J C 76:290, 2016). The fundamental frequencies in this metric are sensitive to some of the metric parameters but insensitive to other parameters, which means, the differences in fundamental frequencies in this metric and those in the Kerr ones can be significant for small changes of some but not all parameters around the Kerr value. By comparing with the QPO observations of GRO J1655-40, we find that only one parameter of the Ghasemi-Nodehi-Bambi metric can be strongly constrained, but other parameters cannot. We also use nested algorithm to investigate whether better constraints on the Ghasemi-Nodehi-Bambi metric parameters can be obtained from QPO observations of multiple objects by simulations. We find that four parameters can be strongly constrained while other parameters cannot . Our results suggest that QPOs may be important tools for testing the Kerr metric.

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.

Since September 2018, LAMOST starts a new 5-year medium-resolution spectroscopic survey (MRS) using bright/gray nights. We present the scientific goals of LAMOST-MRS and propose a near optimistic strategy of the survey. A complete footprint is also provided. Not only the regular medium-resolution survey, but also a time-domain spectroscopic survey is being conducted since 2018 and will be end in 2023. According to the detailed survey plan, we expect that LAMOST-MRS can observe about 2 million stellar spectra with ~7500 and limiting magnitude of around G=15 mag. Moreover, it will also provide about 200 thousand stars with averagely 60-epoch observations and limiting magnitude of G~14 mag. These high quality spectra will give around 20 elemental abundances, rotational velocities, emission line profiles as well as precise radial velocity with uncertainty less than 1 km/s. With these data, we expect that LAMOST can effectively leverage sciences on stellar physics, e.g. exotic binary stars, detailed observation of many types of variable stars etc., planet host stars, emission nebulae, open clusters, young pre-main-sequence stars etc.

The mode switching phenomenon of PSR J0614+2229 was studied by using the archived observations at 686, 1369 and 3100 MHz with the Parkes radio telescope which have not been published before, and combining existing observations from the literature. Over a wide frequency range from 327 to 3100 MHz, the pulsar switches between one mode occurring earlier in pulse phase (mode A) and the other mode appearing later in phase (mode B), with a generally stable phase offset between their profile peaks. The two modes are found to be different in the following aspects. (1) Mode A has a flatter spectrum than mode B, with a difference in the spectral index of about 0.5. This accounts for the phenomenon that the flux ratio between the modes A and B increases with frequency, and mode A becomes stronger than mode B above $\sim500$ MHz. (2) For mode B, the flux density of the subintegrated profile is anticorrelated with the emission phase, indicating that the emission from earlier phases is relatively stronger than that from later phases; such an anticorrelation is not observed in mode A. (3) The frequency dependence of the FWHM of the two modes are opposite each other; namely, the FWHM of mode A increases with frequency, while that of mode B decreases with frequency. A possible interpretation is suggested that the longitudinal spectral variation across the two beams may be opposite each other.

Electrostatic analyzers of different designs have been used since the earliest days of the space age, beginning with the very earliest solar wind measurements made by Mariner 2 en route to Venus in 1962. The Parker Solar Probe (PSP) mission, NASA's first dedicated mission to study the innermost reaches of the heliosphere, makes its thermal plasma measurements using a suite of instruments called the Solar Wind Electrons, Alphas, and Protons (SWEAP) investigation. SWEAP's electron Parker Solar Probe Analyzer (SPAN-E) instruments are a pair of top-hat electrostatic analyzers on PSP that are capable of measuring the electron distribution function in the solar wind from 2 eV to 30 keV. For the first time, in-situ measurements of thermal electrons provided by SPAN-E will help reveal the heating and acceleration mechanisms driving the evolution of the solar wind at the points of acceleration and heating, closer than ever before to the Sun. This paper details the design of the SPAN-E sensors and their operation, data formats, and measurement caveats from Parker Solar Probe's first two close encounters with the Sun.

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

A second peak in extreme ultraviolet sometimes appears during the gradual phase of solar flares, which is known as EUV late phase (ELP). Stereotypically ELP is associated with two separated sets of flaring loops with distinct sizes, and it has been debated whether ELP is caused by additional heating or extended plasma cooling in the longer loop system. Here we carry out a survey of 55 M-and-above GOES-class flares with ELP during 2010--2014. Based on the flare-ribbon morphology, these flares are categorized as circular-ribbon (19 events), two-ribbon (23 events), and complex-ribbon (13 events) flares. Among them, 22 events (40%) are associated with coronal mass ejections, while the rest are confined. An extreme ELP, with the late-phase peak exceeding the main-phase peak, is found in 48% of two-ribbon flares, 37% of circular-ribbon flares, and 31% of complex-ribbon flares, suggesting that additional heating is more likely present during ELP in two-ribbon than in circular-ribbon flares. Overall, cooling may be the dominant factor causing the delay of the ELP peak relative to the main-phase peak, because the loop system responsible for the ELP emission is generally larger than, and well separated from, that responsible for the main-phase emission. All but one of the circular-ribbon flares can be well explained by a composite "dome-plate" quasi-separatrix layer (QSL). Only half of these show a magnetic null point, with its fan and spine embedded in the dome and plate, respectively. The dome-plate QSL, therefore, is a general and robust structure characterizing circular-ribbon flares.

The fundamental stellar atmospheric parameters T_eff and log g and 13 chemical abundances are derived for medium-resolution spectroscopy from LAMOST Medium-Resolution Survey (MRS) data sets with a deep-learning method. The neural networks we designed, named as SPCANet, precisely map LAMOST MRS spectra to stellar parameters and chemical abundances. The stellar labels derived by SPCANet are with precisions of 119 K for T_eff and 0.17 dex for log g. The abundance precision of 11 elements including [C/H], [N/H], [O/H], [Mg/H], [Al/H], [Si/H], [S/H], [Ca/H], [Ti/H], [Cr/H], [Fe/H], and [Ni/H] are 0.06~0.12 dex, while of [Cu/H] is 0.19 dex. These precisions can be reached even for spectra with signal-to-noise as low as 10. The results of SPCANet are consistent with those from other surveys such as APOGEE, GALAH and RAVE, and are also validated with the previous literature values including clusters and field stars. The catalog of the estimated parameters is available at \url{http://paperdata.china-vo.org/LAMOST/MRS_parameters_elements.csv}.

Galaxy mergers and interactions are expected to play a significant role leading to offsets between gas and stellar motions in galaxies. Herein we crossmatch galaxies in MaNGA MPL-8 with the Dark Energy Spectroscopic Instrument (DESI) Legacy Surveys and identify 538 galaxies with merging/interacting features to investigate their position angle offsets ($\Delta$PA) between gas and stellar rotation. We find that there is a much higher merging/interacting fraction in misaligned galaxies (30$^\circ$ $\leqslant$ $\Delta$PA $<$150$^\circ$) than that in co-rotators ($\Delta$PA $<$ 30$^\circ$). This result corroborates that merging/interacting is one of the sources to produce misaligned galaxies and recent merging events can contribute to maximum 40\% of such misalignment. Furthermore, the marginal merging/interacting fraction in star-forming (SF) counter-rotators ($\Delta$PA $\geqslant$ 150$^\circ$) indicates the minor role of merging in the origin of SF counter-rotators. In addition, the slightly smaller merging/interacting fraction in non star-forming (non-SF) counter-rotators (0.14$\pm$0.039) than that in non-SF misaligned galaxies (0.22$\pm$0.031) agrees with that counter-rotating is a stable state, where tidal features disappear. Finally, the ratio of co-rotators to counter-rotators in non-SF merging/interacting galaxies is about 8:1, much larger than the prediction from the isotropic merging (1:1), which supports the speculation that gas and stars prefer to be aligned during merging, as the orbital angular momentum transfers to gas and stellar spin.

We present and analyze the optical/UV and X-ray observations of a nearby tidal disruption event (TDE) candidate AT2019azh, spanning from 30 d before to ~ 250 d after its early optical peak. The X-rays show a late brightening by a factor of ~ 30-100 around 250 days after discovery, while the UV/opticals continuously decayed. The early X-rays show two flaring episodes of variation, temporally uncorrelated with the early UV/opticals. We found a clear sign of X-ray hardness evolution, i.e., the source is harder at early times, and becomes softer as it brightens later. The drastically different temporal behaviors in X-rays and UV/opticals suggest that the two bands are physically distinct emission components, and probably arise from different locations. These properties argue against the reprocessing of X-rays by any outflow as the origin of the UV/optical peak. The full data are best explained by a two-process scenario, in which the UV/optical peak is produced by the debris stream-stream collisions during the circularization phase; some low angular momentum, shocked gas forms an early, low-mass accretion disk which emits the early X-rays. The major body of the disk is formed after the circularization finishes, whose enhanced accretion rate produces the late X-ray brightening. AT2019azh is a strong case of TDE whose emission signatures of stream-stream collision and delayed accretion are both identified.

In this work, using the Gaussian Process, we explore the potentiality of future gravitational wave (GW) measurement to probe cosmic opacity through comparing its opacity-free luminosity distance (LD) with the opacity-dependent one from type Ia supernovae (SNIa). GW data points are simulated from the third generation Einstein Telescope, and SNIa data are taken from the Joint Light Analysis (JLA) or Pantheon compilation. The advantages of using Gaussian Process are that one may match SNIa data with GW data at the same redshift and use all available data to probe cosmic opacity. We obtain that the error bar of the constraint on cosmic opacity can be reduced to $\sigma_{\epsilon}\sim 0.011$ and $0.006$ at $1\sigma$ confidence level (CL) for JLA and Pantheon respectively in a cosmological-independent way. Thus, the future GW measurements can give competitive results on the cosmic opacity test. Furthermore, we propose a method to probe the spatial homogeneity of the cosmic transparency through comparing the reconstructed LD from the mock GW with the reconstructed one from SNIa data in a flat $\Lambda$CDM with the Gaussian Process. The result shows that a transparent universe is favored at $1\sigma$ CL, although the best-fit value of cosmic opacity is redshift-dependent.

By analysing the data acquired from the Parkes 64-m radio telescope at 1369 MHz, we report on the phase-stationary non-drift amplitude modulation observed in PSR J1048-5832. The high-sensitivity observations revealed that the central and trailing components of the pulse profile of this pulsar switch between a strong mode and a weak mode periodically. However, the leading component remains unchanged. Polarization properties of the strong and weak modes are investigated. Considering the similarity to mode changing, we argue that the periodic amplitude modulation in PSR J1048$-$5832 is periodic mode changing. The fluctuation spectral analysis showed that the modulation period is very short (~2.1 s or 17 P1), where P1 is the rotation period of the pulsar. We find that this periodic amplitude modulation is hard to explain by existing models that account for the periodic phenomena in pulsars like subpulse drifting.

The radial velocity (RV) is a basic physical quantity which can be determined through Doppler shift of the spectrum of a star. The precision of RV measurement depends on the resolution of the spectrum we used and the accuracy of wavelength calibration. In this work, radial velocities of LAMOST-II medium resolution (R ~ 7500) spectra are measured for 1,594,956 spectra (each spectrum has two wavebands) through matching with templates. A set of RV standard stars are used to recalibrate the zero point of the measurement, and some reference sets with RVs derived from medium/high-resolution observations are used to evaluate the accuracy of the measurement. Comparing with reference sets, the accuracy of our measurement can get 0.0227 km s/1 with respect to radial velocities standard stars. The intrinsic precision is estimated with the multiple observations of single stars, which can achieve to 1.36 km s/1,1.08 km s/1, 0.91 km s/1 for the spectra at signal-to-noise levels of 10, 20, 50, respectively.

We present the Global Rapid Advanced Network Devoted to the Multi-messenger Addicts (GRANDMA). The network consists of 21 telescopes with both photometric and spectroscopic facilities. They are connected together thanks to a dedicated infrastructure. The network aims at coordinating the observations of large sky position estimates of transient events to enhance their follow-up and reduce the delay between the initial detection and the optical confirmation. The GRANDMA program mainly focuses on follow-up of gravitational-wave alerts to find and characterise the electromagnetic counterpart during the third observational campaign of the Advanced LIGO and Advanced Virgo detectors. But it allows for any follow-up of transient alerts involving neutrinos or gamma-ray bursts, even with poor spatial localisation. We present the different facilities, tools, and methods we developed for this network, and show its efficiency using observations of LIGO/Virgo S190425z, a binary neutron star merger candidate. We furthermore report on all GRANDMA follow-up observations performed during the first six months of the LIGO-Virgo observational campaign, and we derive constraints on the kilonova properties assuming that the events' locations were imaged by our telescopes.

We obtained 7 spectra of the Be star V423 Aur on Dec. 5th, 2017 using the LAMOST Medium-Resolution Spectrograph with exposures from 600 to 1200 seconds. These spectra show the irregular H{\alpha} emission line profile variations (LPVs). In the seven spectra, from the 4th to 7th, the left part of H{\alpha} profile even show excess. However, no variation can be seen from the follow-up observation of photometry by 1.26-m telescope and High-Resolution spectra by 2.16-m telescope. According to the High-Resolution spectra, we conclude that it is a B7V type star with E(B-V)= 0.709{\pm}0.036 and its vsini is {\sim} 221.8km/s. The short-term H{\alpha} LPVs could be explained as a result of the transient ejection of matter from rotating disk or shell around V423 Aur.

The precise measurement of the spectrum of protons, the most abundant component of the cosmic radiation, is necessary to understand the source and acceleration of cosmic rays in the Milky Way. This work reports the measurement of the cosmic ray proton fluxes with kinetic energies from 40 GeV to 100 TeV, with two and a half years of data recorded by the DArk Matter Particle Explorer (DAMPE). This is the first time an experiment directly measures the cosmic ray protons up to ~100 TeV with a high statistics. The measured spectrum confirms the spectral hardening found by previous experiments and reveals a softening at ~13.6 TeV, with the spectral index changing from ~2.60 to ~2.85. Our result suggests the existence of a new spectral feature of cosmic rays at energies lower than the so-called knee, and sheds new light on the origin of Galactic cosmic rays.

In this work, we revisit the size-luminosity relation of the extended narrow line regions (ENLRs) using a large sample of nearby active galactic nuclei (AGN) from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey. The ENLRs ionized by the AGN are identified through the spatially resolved BPT diagram, which results in a sample of 152 AGN. By combining our AGN with the literature high-luminosity quasars, we found a tight log-linear relation between the size of the ENLR and the AGN [O III]{\lambda}5007{\AA} luminosity over four orders of magnitude of the [O III] luminosity. The slope of this relation is 0.42 $\pm$ 0.02 which can be explained in terms of a distribution of clouds photoionized by the AGN. This relation also indicates the AGN have the potential to ionize and heat the gas clouds at a large distance from the nuclei without the aids of outflows and jets for the low-luminosity Seyferts.

Cosmological tensions can arise within $\Lambda$CDM scenario amongst different observational windows, such as measurements of the $H_0$ and $\sigma_8$ parameters. These tensions, if finally confirmed by measurements, may indicate new physics beyond the standard paradigm. In this Letter, we report how to alleviate both the $H_0$ and $\sigma_8$ tensions simultaneously within torsional gravity from the perspective of effective field theory (EFT). We apply the EFT approach, which allows to investigate the evolution equations at the background and perturbation levels in a systematic way, and examine the conditions followed by the coefficients of various possibly involved operators such that cosmological tensions can be relaxed. Following these observations we construct concrete models of Lagrangians of torsional gravity. Specifically, we consider the parametrization $f(T)=-T -2\Lambda/M_P^2 +\alpha T^\beta$, where two out of the three parameters are independent (in which an additional term of the form $cT^{1/2}$ can be added). This model can efficiently fit observations solving all statistical tensions. To our knowledge, this is the first time where a modified gravity theory can alleviate both $H_0$ and $\sigma_8$ tensions simultaneously, hence offering an additional argument in favor of gravitational modification.

The Analytical Blind Separation (ABS) method is dedicated to extract the CMB signal from its contaminants without any prior information about foregrounds. In this study, we apply the ABS method to simulated polarization sky maps in order to test its ability in the recovery of the CMB $E$- and $B$-mode power spectra. We simulate multi-frequency $Q$ and $U$ microwave sky maps in the presence of two polarized foreground components (synchrotron and thermal dust), together with an uncorrelated Gaussian distributed instrumental noise, measured with a hypothetical PICO-like experiment operating in 10-frequency bands in the range of 30 to 321 GHz. Both full and partial sky surveys are taken into account in our analysis. For the full sky case with the tensor-to-scalar ratios of $r=0$ and $r=0.05$, the ABS can accurately recover the $E$-mode and $B$-mode power spectra with a relative error below $20\%$ with respect to the input CMB power spectra for the full multipole range ($2 \leq \ell \leq 1050$), close to theoretical expectations. For the partial sky analysis with a E/B separation proposed by \cite{ebmixture10} to avoid "$E$-to-$B$" leakage, the results are comparable to the full sky case. We find that the ABS is able to estimate both CMB $E$-mode and $B$-mode power spectra within 1-$\sigma$ at all scales for full sky, and for most of the scales ($50 \leq \ell \leq 1050$) for partial sky. For low-$\ell$ multipoles ($\ell \leq 50$), a noticeable divergence occurs in the partial sky analysis, which might be due to uncertainties originated from both the ABS and the E/B separation. A more detailed analysis for the partial sky case is thus needed for estimating primordial $B$-modes unbiasedly and accurately.

In this paper, we investigate the energy partition of two homologous M1.1 circular-ribbon flares (CRFs) in active region (AR) 12434. They were observed by \textit{SDO}, \textit{GOES}, and \textit{RHESSI} on 2015 October 15 and 16, respectively. The peak thermal energy, nonthermal energy of flare-accelerated electrons, total radiative loss of hot plasma, and radiant energies in 1$-$8 {\AA} and 1$-$70 {\AA} of the flares are calculated. The two flares have similar energetics. The peak thermal energies are (1.94$\pm$0.13)$\times$10$^{30}$ erg. The nonthermal energies in flare-accelerated electrons are (3.9$\pm$0.7)$\times$10$^{30}$ erg. The radiative outputs of the flare loops in 1$-$70 {\AA}, which are $\sim$200 times greater than the outputs in 1$-$8 {\AA}, account for $\sim$62.5\% of the peak thermal energies. The radiative losses of SXR-emitting plasma are one order of magnitude lower than the peak thermal energies. Therefore, the total heating requirements of flare loops including radiative loss are (2.1$\pm$0.1)$\times$10$^{30}$ erg, which could sufficiently be supplied by nonthermal electrons.

The DArk Matter Particle Explorer (DAMPE), a satellite-based cosmic ray and gamma-ray detector, was launched on December 17, 2015, and began its on-orbit operation on December 24, 2015. In this work we document the on-orbit calibration procedures used by DAMPE and report the calibration results of the Plastic Scintillator strip Detector (PSD), the Silicon-Tungsten tracKer-converter (STK), the BGO imaging calorimeter (BGO), and the Neutron Detector (NUD). The results are obtained using Galactic cosmic rays, bright known GeV gamma-ray sources, and charge injection into the front-end electronics of each sub-detector. The determination of the boundary of the South Atlantic Anomaly (SAA), the measurement of the live time, and the alignments of the detectors are also introduced. The calibration results demonstrate the stability of the detectors in almost two years of the on-orbit operation.

Type Ia supernovae (SNe Ia) play key roles in revealing the accelerating expansion of the universe, but our knowledge about their progenitors is still very limited. Here we report the discovery of a rigid dichotomy in circumstellar (CS) environments around two subclasses of type Ia supernovae (SNe Ia) as defined by their distinct photospheric velocities. For the SNe Ia with high photospheric velocities (HV), we found a significant excess flux in blue light during 60-100 days past maximum, while this phenomenon is absent for SNe with normal photospheric velocity (Normal). This blue excess can be attributed to light echoes by circumstellar dust located at a distance of about 1-3x10^{17} cm from the HV subclass. Moreover, we also found that the HV SNe Ia show systematically evolving Na I absorption line by performing a systematic search of variable Na I absorption lines in spectra of all SNe Ia, whereas this evolution is rarely seen in Normal ones. The evolving Na I absorption can be modeled in terms of photoionization model, with the location of the gas clouds at a distance of about 2x10^{17} cm, in striking agreement with the location of CS dust inferred from B-band light curve excess. These observations show clearly that the progenitors of HV and Normal subclasses are systematically different, suggesting that they are likely from single and double degenerate progenitor systems, respectively.

We analyze the intrinsic velocity dispersion properties of 648 star-forming galaxies observed by the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey, to explore the relation of intrinsic gas velocity dispersions with star formation rates (SFRs), SFR surface densities ($\rm{\Sigma_{SFR}}$), stellar masses and stellar mass surface densities ($\rm{\Sigma_{*}}$). By combining with high z galaxies, we found that there is a good correlation between the velocity dispersion and the SFR as well as $\rm{\Sigma_{SFR}}$. But the correlation between the velocity dispersion and the stellar mass as well as $\rm{\Sigma_{*}}$ is moderate. By comparing our results with predictions of theoretical models, we found that the energy feedback from star formation processes alone and the gravitational instability alone can not fully explain simultaneously the observed velocity-dispersion/SFR and velocity-dispersion/$\rm{\Sigma_{SFR}}$ relationships.

Understanding the sources of lunar water is crucial for studying the history of lunar evolution, and also the solar wind interaction with the Moon and other airless bodies. Recent observations revealed lunar hydration is very likely a surficial dynamic process driven by solar wind. Solar wind is shielded over a period of 3-5 days as the Moon passes through the Earth's magnetosphere, during which a significant loss of hydration is expected from previous works.Here we study lunar hydration inside the magnetosphere using orbital spectral data, which unexpectedly found that the polar surficial OH/H2O abundance remains at the same level when in the solar wind and in the magnetosphere. We suggest that particles from the magnetosphere (Earth wind, naturally different from solar wind) contribute to lunar hydration. From lunar orbital plasma observations, we find the existence of optimal energy ranges, other than 1 keV as previously thought, for surface hydration formation. These optimal energy ranges deduced from space observations may provide strong implications for laboratory experiments simulating lunar hydration processes.

We explore the possibility of detecting entangled photon pairs from cosmic microwave background or other cosmological sources coming from two patches of the sky. The measurements use two detectors with different photon polarizer directions. When two photon sources are separated by a large angle relative to the earth, such that each detector has only one photon source in its field of view, a null test of unentangled photons can be performed. The deviation from this unentangled background is, in principle, the signature of photon entanglement. To confirm whether the deviation is consistent with entangled photons, we derive a photon polarization correlation to compare with, similar to that in a Bell inequality measurement. However, since photon coincidence measurement cannot be used to discriminate unentangled cosmic photons, it is unlikely that the correlation expectation value alone can violate Bell inequality to provide the signature for entanglement.

The Near-infrared Sky Brightness Monitor (NIRBM) aims to measure the middle infrared sky background in Antarctica. The NIRBM mainly consists of an InGaAs detector, a chopper, a reflector, a cooler and a black body. The reflector can rotate to scan the sky with a field of view ranging from 0{\deg} to 180{\deg}. Electromechanical control and weak signal readout functions are accomplished by the same circuit, whose core chip is a STM32F407VG microcontroller. Considering the environment is harsh for humans in Antarctica, a multi-level remote control software system is designed and implemented. A set of EPICS IOCs are developed to control each hardware module independently via serial port communication with the STM32 microcontroller. The tornado web framework and PyEpics are introduced as a combination where PyEpics is used to monitor or change the EPICS Process Variables, functioning as a client for the EPICS framework. Tornado is responsible for the specific operation process of inter-device collaboration, and expose a set of interfaces to users to make calls. Considering the high delay and low bandwidth of the network environment, the tornado back-end is designed as a master-and-agent architecture to improve domestic user experience. The master node is deployed in Antarctic while multiple agent nodes can be deployed domestic. The master and agent nodes communicate with each other through the WebSocket protocol to exchange latest information so that bandwidth is saved. The GUI is implemented in the form of single-page application based on the Vue framework which communicates with tornado through WebSocket and AJAX requests. The web page integrates device control, data curve drawing, alarm display, auto observation and other functions together.

Our study aims to recognize M-type stars which are classified as "UNKNOWN" due to bad quality in Large sky Area Multi-Object fibre Spectroscopic Telescope (LAMOST) DR5 V1. A binary nonlinear hashing algorithm based on Multi-Layer Pseudo Inverse Learning (ML-PIL) is proposed to effectively learn spectral features for the M-type star detection, which can overcome the bad fitting problem of template matching, particularly for low S/N spectra. The key steps and the performance of the search scheme are presented. A positive dataset is obtained by clustering the existing M-type spectra to train the ML-PIL networks. By employing this new method, we find 11,410 M-type spectra out of 642,178 "UNKNOWN" spectra, and provide a supplemental catalogue. Both the supplemental objects and released M-type stars in DR5 V1 are composed a whole M type sample, which will be released in the official DR5 to the public in June 2019, All the M-type stars in the dataset are classified to giants and dwarfs by two suggested separators: 1) color diagram of H versus J~K from 2MASS; 2) line indices CaOH versus CaH1, and the separation is validated with HRD derived from Gaia DR2. The magnetic activities and kinematics of M dwarfs are also provided with the EW of H_alpha emission line and the astrometric data from Gaia DR2 respectively.

In this paper, we present extensive optical and ultraviolet (UV) observations of the type Ia supernova (SN Ia) 2013gs discovered during the Tsinghua-NAOC Transient Survey. The photometric observations in the optical show that the light curves of SN 2013gs is similar to that of normal SNe Ia, with an absolute peak magnitude of $M_{B}$ = $-$19.25 $\pm$ 0.15 mag and a post-maximum decline rate $\Delta$m$_{15}$(B) = 1.00 $ \pm $ 0.05 mag. \emph{Gehrels Swift} UVOT observations indicate that SN 2013gs shows unusually strong UV emission (especially in the $uvw1$ band) at around the maximum light (M$_{uvw1}$ $\sim$ $-$18.9 mag). The SN is characterized by relatively weak Fe~{\sc ii} {\sc iii} absorptions at $\sim$ 5000{\AA} in the early spectra and a larger expansion velocity ($v_{Si}$ $\sim$ 13,000 km s$^{-1}$ around the maximum light) than the normal-velocity SNe Ia. We discuss the relation between the $uvw1-v$ color and some observables, including Si~{\sc ii} velocity, line strength of Si~{\sc ii} $\lambda$6355, Fe~{\sc ii}/{\sc iii} lines and $\Delta m_{15}$(B). Compared to other fast-expanding SNe Ia, SN 2013gs exhibits Si and Fe absorption lines with similar strength and bluer $uvw1-v$ color. We briefly discussed the origin of the observed UV dispersion of SNe Ia.

Radial velocity is one of key measurements in understanding the fundamental properties of stars, stellar clusters and the Galaxy. A plate of stars in the Kepler field were observed in May of 2018 with the medium-resolution spectrographs of LAMOST, aiming to test the performance of this new system which is the upgraded equipment of LAMOST after the first five-year regular survey.We present our analysis on the radial velocity measurements (RVs) derived from these data. The results show that slight and significant systematic errors exist among the RVs obtained from the spectra collected by different spectrographs and exposures, respectively. After correcting the systematic errors with different techniques, the precision of RVs reaches ~1.3, ~1.0, ~0.5 and ~0.3 km/s at S/Nr = 10, 20, 50, and 100, respectively. Comparing with the RVs of the standard stars of the APOGEE survey, our RVs are calibrated with a zero-point shift of ~7 km/s. The results indicate that the LAMOST medium-resolution spectroscopic system may provide RVs in a reasonable accuracy and precision for the selected targets.

Radial velocity is one of key measurements in understanding the fundamental properties of stars, stellar clusters and the Galaxy. A plate of stars in the Kepler field were observed in May of 2018 with the medium-resolution spectrographs of LAMOST, aiming to test the performance of this new system which is the upgraded equipment of LAMOST after the first five-year regular survey.We present our analysis on the radial velocity measurements (RVs) derived from these data. The results show that slight and significant systematic errors exist among the RVs obtained from the spectra collected by different spectrographs and exposures, respectively. After correcting the systematic errors with different techniques, the precision of RVs reaches ~1.3, ~1.0, ~0.5 and ~0.3 km/s at S/Nr = 10, 20, 50, and 100, respectively. Comparing with the RVs of the standard stars of the APOGEE survey, our RVs are calibrated with a zero-point shift of ~7 km/s. The results indicate that the LAMOST medium-resolution spectroscopic system may provide RVs in a reasonable accuracy and precision for the selected targets.

We present an empirical stellar spectra library created using spectra from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) DR5. This library represents a uniform data set ranging from 3750 through 8500K in effective temperature (T_eff), from -2.5 through +1.0 dex in metallicity ([Fe/H]), and from 0 to 5.0 dex in gravity (log g). The spectra in the library have resolution R~1800, with well-calibrated fluxes and rest-framed wavelengths. Using a large number of red stars observed by LAMOST, we generated denser K type templates to fill in data missing from current empirical spectral libraries, particularly the late K type. For K giants, we calibrated the spectroscopic surface gravities against the asteroseismic surface gravities. To verify the reliability of the parameters labeled for this library, we performed an internal cross-validation by using a chi^2 minimization method to interpolate parameters of each individual spectrum using the remaining spectra in the library. We obtained precisions of 41 K, 0.11 dex, and 0.05 dex for T_eff, log g, and [Fe/H], respectively, which means the templates are labeled with correct stellar parameters. Through external comparisons, we confirmed that measurements of the stellar parameters through this library can achieve accuracies of approximately 125K in T_eff, 0.1 dex in [Fe/H] and 0.20 dex in log g without systematic offset. This empirical library is useful for stellar parameter measurements because it has large parameter coverage and full wavelength coverage from 3800 to 8900A.

There has been suggestive evidence of intermediate-mass black holes (IMBHs; 10^{3-5} M_sun) existing in some globular clusters (GCs) and dwarf galaxies, but IMBHs as a population remain elusive. As a main-sequence star passes too close by an IMBH it might be tidally captured and disrupted. We study the long-term accretion and observational consequence of such tidal disruption events. The disruption radius is hundreds to thousands of the BH's Schwarzschild radius, so the circularization of the falling-back debris stream is very inefficient due to weak general relativity effects. Due to this and a high mass fallback rate, the bound debris initially goes through a ~10 yr long super-Eddington accretion phase. The photospheric emission of the outflow ejected during this phase dominates the observable radiation and peaks in the UV/optical bands with a luminosity of 10^42 erg/s. After the accretion rate drops below the Eddington rate, the bolometric luminosity follows the conventional t^{-5/3} power-law decay, and X-rays from the inner accretion disk start to be seen. Modeling the newly reported IMBH tidal disruption event candidate 3XMM J2150-0551, we find a general consistency between the data and predictions. The search for these luminous, long-term events in GCs and nearby dwarf galaxies could unveil the IMBH population.

Type Ia supernova cosmology depends on the ability to fit and standardize observations of supernova magnitudes with an empirical model. We present here a series of new models of Type Ia Supernova spectral time series that capture a greater amount of supernova diversity than possible with the models that are currently customary. These are entitled SuperNova Empirical MOdels (\textsc{SNEMO}\footnote{https://snfactory.lbl.gov/snemo}). The models are constructed using spectrophotometric time series from $172$ individual supernovae from the Nearby Supernova Factory, comprising more than $2000$ spectra. Using the available observations, Gaussian Processes are used to predict a full spectral time series for each supernova. A matrix is constructed from the spectral time series of all the supernovae, and Expectation Maximization Factor Analysis is used to calculate the principal components of the data. K-fold cross-validation then determines the selection of model parameters and accounts for color variation in the data. Based on this process, the final models are trained on supernovae that have been dereddened using the Fitzpatrick and Massa extinction relation. Three final models are presented here: \textsc{SNEMO2}, a two-component model for comparison with current Type~Ia models; \textsc{SNEMO7}, a seven component model chosen for standardizing supernova magnitudes which results in a total dispersion of $0.100$~mag for a validation set of supernovae, of which $0.087$~mag is unexplained (a total dispersion of $0.113$~mag with unexplained dispersion of $0.097$~mag is found for the total set of training and validation supernovae); and \textsc{SNEMO15}, a comprehensive $15$ component model that maximizes the amount of spectral time series behavior captured.

Blazar OJ 287 is a candidate nanoHertz (nHz) gravitational wave (GW) source. In this article, we investigate the GWs generated by OJ 287 and their potential detection through a pulsar timing array (PTA). First, we obtain the orbit and the corresponding GW strain of OJ 287. During the time span of the next 10 years (2019 to 2029), the GW of OJ 287 will be active before 2021, with a peak strain amplitude $8 \times 10^{-16}$, and then decay after that. When OJ 287 is silent in the GW channel during 2021 to 2029, the timing residual signals of the PTA will be dominated by the 'pulsar term' of the GW strain and this provides an opportunity to observe this pulsar term. Furthermore, we choose 26 pulsars with white noise below 300 ns to detect the GW signal of OJ 287, evaluating their timing residuals and signal-to-noise ratios (SNRs). The total SNR (with a cadence of 2 weeks in the next 10 years) of the PTA ranges from 1.9 to 2.9, corresponding to a weak GW signal for the current sensitivity level. Subsequently, we investigate the potential measurement of the parameters of OJ 287 using these pulsars. In particular, PSR J0437-4715, with a precisely measured distance, has the potential to constrain the polarization angle with an uncertainty below $8^{\deg}$ and this pulsar will play an important role in future PTA observations.

Magnetohydrodynamic(MHD) turbulence displays velocity anisotropies which reflect the direction of the magnetic field. This anisotropy has led to the development of a number of statistical techniques for studying magnetic fields in the interstellar medium. In this paper, we review and compare three techniques that use radio position-position-velocity data for determining magnetic field strength and morphology : the correlation function anisotropy (CFA), Principal Component Analysis of Anisotropies (PCAA), and the more recent Velocity Gradient Technique (VGT). We compare these three techniques and suggest improvements to the CFA and PCAA techniques to increase their accuracy and versatility. In particular, we suggest and successfully implement a much faster way of calculating non-periodic correlation functions for the CFA. We discuss possible improvements to the current implementation of the PCAA. We show the advantages of the VGT in terms of magnetic field tracing and stress the complementary nature with the other two techniques.

Early-type galaxies (ETGs) are composed of two distinct populations: high-mass and low-mass, which are likely to be built via gas-poor merging and gas-rich merging/accretion, respectively. However, it is difficult to directly associate low-mass ETGs with gas-rich processes, because currently they are gas poor with no signs of ongoing star formation. We report a discovery of an ETG (SDSS J142055.01+400715.7) with Mstellar=10^10 Msun that offers direct evidence for gas-rich merging as the origin of low-mass ETGs. The integrated properties of the galaxy are consistent with a typical low-mass ETG, but the outer and inner regions show distinct dispersion- and rotation-dominated kinematics, respectively. There are some tidal features surrounding the galaxy. These two facts suggest very recent galaxy merging. Furthermore, the inner disk harbors on-going star formation, indicating the merging to be gas rich. This type of galaxy is rare but it may be a demonstration of the role the transient phase of gas-rich merging plays in making a low-mass ETG.

Recent developments of the Velocity Gradient Technique (VGT) show that the velocity gradients provide a reliable tracing of magnetic field direction in turbulent plasmas. In this paper, we explore the ability of velocity gradients to measure the magnetization of interstellar medium. We demonstrate that the distribution of velocity gradient orientations provides a reliable estimation of the magnetization of the media. In particular, we determine the relation between Alfvenic Mach number $M_A$ in the range of $M_A \in [0.2,1.7]$ and properties of the velocity gradient distribution, namely, with the dispersion of velocity gradient orientation as well as with the peak to base ratio of the amplitudes. We apply our technique for a selected GALFA-HI region and find the results consistent with the expected behavior of $M_A$. Using 3D MHD simulations we successfully compare the results with our new measure of magnetization that is based on the dispersion of starlight polarization. We demonstrate that, combined with the velocity dispersion along the line of sight direction, our technique is capable to delivering the magnetic field strength. The new technique opens a way to measure magnetization using other gradient measures such as synchrotron intensity gradients (SIGs) and synchrotron polarization gradients (SPGs).

A long-standing issue in peculiar velocity cosmology is whether the halo/galaxy velocity bias $b_v=1$ at large scale. The resolution of this important issue must resort to high precision cosmological simulations. However, this is hampered by another long-standing `sampling artifact' problem in volume weighted velocity measurement. We circumvent this problem with a hybrid approach. We first measure statistics free of sampling artifact, then link them to volume weighted statistics in theory, finally solve for the velocity bias. $b_v$ determined by our method is not only free of sampling artifact, but also free of cosmic variance. We apply this method to a $\Lambda$CDM N-body simulation of $3072^3$ particles and $1200 Mpc/{\rm h}$ box size. For the first time, we determine the halo velocity bias to $0.1\%$-$1\%$ accuracy. Our major findings are as follows: (1) $b_v\neq 1$ at $k>0.1 h/{\rm Mpc}$. The deviation from unity ($|b_v-1|$) increases with $k$. Depending on halo mass and redshift, it may reach $\mathcal{O}(0.01)$ at $k=0.2 h/{\rm Mpc}$ and $\mathcal{O}(0.05)$ at $k\sim 0.3 h/{\rm Mpc}$. The discovered $b_v\neq 1$ has statistically significant impact on structure growth rate measurement by spectroscopic redshift surveys, including DESI, Euclid and SKA. (2) Both the sign and the amplitude of $b_v-1$ depend on mass and redshift. These results disagree with the peak model prediction in that $b_v$ has much weaker deviation from unity, varies with redshift, and can be bigger than unity. (3) Most of the mass and redshift dependences can be compressed into a single dependence on the halo density bias. Based on this finding, we provide an approximate two-parameter fitting formula.

MARST (Multi-channel Antarctic Solar Telescope) aims to observe the Sun in multiple wavebands in Antarctica and will be China's first solar telescope in Antarctica. The telescope has two tubes, corresponding to Photosphere observation which uses 11 filters and Chromosphere observation in H{\alpha} waveband. The two tubes need to be coordinated to observe at the same time. The telescope will observe the Sun for a long time, so a self-guiding module is needed to improve sun tracking. Besides, performing solar specific flat-field exposure is necessary for analyzing. EPICS is introduced to control each hardware and an autonomous observation system based on RTS2 is designed under such demands. EPICS application modules are implemented for each device: telescope mount & focuser, filter wheel, H{\alpha} filter, dome with webcams, Andor CCD and PI CCD. We also integrate EPICS modules into RTS2 framework with an XML format configuration. To control these applications autonomously we have developed a RTS2 executor module where two plan classes are instantiated to control two sets of filters and CCDs, and to ensure only one could control the mount at the same time. Different types of observation plans are designed to describe different series of processes with different priorities. To improve sun tracking, we calculate the centroid of each image to get the offset, then apply the correction to the telescope during observation process. For frontend users, a GUI based on PyQt5 and QML is implemented and connects to rts2-httpd and rts2-proxy modules so that users can control devices, check images and get logs.

Life appears to have emerged relatively quickly on the Earth, a fact sometimes used to justify a high rate of spontaneous abiogenesis ($\lambda$) among Earth-like worlds. Conditioned upon a single datum - the time of earliest evidence for life ($t_{\mathrm{obs}}$) - previous Bayesian formalisms for the posterior distribution of $\lambda$ have demonstrated how inferences are highly sensitive to the priors. Rather than attempt to infer the true $\lambda$ posterior, we here compute the relative change to $\lambda$ when new experimental/observational evidence is introduced. By simulating posterior distributions and resulting entropic information gains, we compare three experimental pressures on $\lambda$: 1) evidence for an earlier start to life; $t_{\mathrm{obs}}$; 2) constraints on spontaneous abiogenesis from the lab; and 3) an exoplanet survey for biosignatures. First, we find that experiments 1 and 2 can only yield lower limits on $\lambda$, unlike 3. Second, evidence for an earlier start to life can yield negligible information on $\lambda$ if $t_{\mathrm{obs}} \ll \lambda_{\mathrm{max}}^{-1}$. Vice versa, experiment 2 is uninformative when $\lambda_{\mathrm{max}} \gg t_{\mathrm{obs}}^{-1}$. Whilst experiment 3 appears the most direct means of measuring $\lambda$, we highlight that early starts inform us of the conditions of abiogenesis, and that lab experiments could succeed in building new life. Together then, the three experiments are complimentary and we encourage activity in all to solve this grand challenge.

Tibet is known as the third pole of the earth, as high as the South Pole and North Pole. The Ngari (Ali) observatory in Tibet has the advantage of plenty of photometric night, low precipitable water vapor, high transmittance, good seeing. It is a good site, and promising to be one of the best place for infrared and submillimeter observations in the world. However, there is no data available for sky background brightness in such place, which restrict the astronomical development of the sites. In the near infrared band of J, H, Ks, a NIR sky brightness monitor (NISBM) is designed based on InGaAs photoelectric diode. By using the method of chopper modulation and digital lock-in amplifier processing, the SNR (Signal Noise Ratio), detectivity and the data acquisition speed of the device is greatly improved. For each band of J, H, Ks, an independent instrument is designed and calibrated in laboratory. The NISBM has been installed in Ngari observatory in July of 2017 and obtained the first data of NIR sky brightness at Ngari observatory.

Tibet is known as the third pole of the earth, as high as the South Pole and North Pole. The Ngari (Ali) observatory in Tibet has the advantage of plenty of photometric night, low precipitable water vapor, high transmittance, good seeing. It is a good site, and promising to be one of the best place for infrared and submillimeter observations in the world. However, there is no data available for sky background brightness in such place, which restrict the astronomical development of the sites. In the near infrared band of J, H, Ks, a NIR sky brightness monitor (NISBM) is designed based on InGaAs photoelectric diode. By using the method of chopper modulation and digital lock-in amplifier processing, the SNR (Signal Noise Ratio), detectivity and the data acquisition speed of the device is greatly improved. For each band of J, H, Ks, an independent instrument is designed and calibrated in laboratory. The NISBM has been installed in Ngari observatory in July of 2017 and obtained the first data of NIR sky brightness at Ngari observatory.

A near-infrared (NIR) measurement based on digital orthogonal vector lock-in amplifier (LIA) is present in this paper. NIR sky background radiation is very weak. To detect the signals obscured by noise, the best way achieved is to use a chopper to modulate the detected signal and using a LIA to demodulate. The effect of 1/f noise of detector, dark current and other noises can be reduced to get sufficient signal-to-noise ratio (SNR). The orthogonal vector LIA can avoid the phase shift on the accuracy of measurement by two orthogonal components. In order to simplify the system, a digital algorithm is adopted to realize the LIA which is operated in a microcontroller with ARM cortex-M4. Data is obtained through ADC and the signal of detector is amplified and filtered. Then the phase sensitive detection (PSD), low-pass filter (LPF) and amplitude phase calculation are performed. The digital method can greatly simplify the circuit, and conveniently adjust the time constant of the LPF to realize the different equivalent noise bandwidth (ENB). The algorithm has the specification of high precision, flexible usage, simple implementation and low computation resource. By using this method, the weak infrared signal submerged by the noise can be obtained, which extremely improves the detection capability of the system.

Protoplanets may be born into dust-rich environments if planetesimals formed through streaming or gravitational instabilities, or if the protoplanetary disc is undergoing mass loss due to disc winds or photoevaporation. Motivated by this possibility, we explore the interaction between low mass planets and dusty protoplanetary discs with focus on disc-planet torques. We implement Lin & Youdin's newly developed, purely hydrodynamic model of dusty gas into the PLUTO code to simulate dusty protoplanetary discs with an embedded planet. We find that for imperfectly coupled dust and high metallicity, e.g. Stokes number $10^{-3}$ and dust-to-gas ratio $\Sigma_\mathrm{d}/\Sigma_\mathrm{g}=0.5$ , a `bubble' develops inside the planet's co-orbital region, which introduces unsteadiness in the flow. The resulting disc-planet torques sustain large amplitude oscillations that persists well beyond that in simulations with perfectly coupled dust or low dust-loading, where co-rotation torques are always damped. We show that the desaturation of the co-rotation torques by finite-sized particles is related to potential vorticity generation from the misalignment of dust and gas densities. We briefly discuss possible implications for the orbital evolution of protoplanets in dust-rich discs. We also demonstrate Lin & Youdin's dust-free framework reproduces previous results pertaining to dusty protoplanetary discs, including dust-trapping by pressure bumps, dust settling, and the streaming instability.

The interior structure of an exoplanet is hidden from direct view yet likely plays a crucial role in influencing the habitability of the Earth analogues. Inferences of the interior structure are impeded by a fundamental degeneracy that exists between any model comprising more than two layers and observations constraining just two bulk parameters: mass and radius. In this work, we show that although the inverse problem is indeed degenerate, there exists two boundary conditions that enables one to infer the minimum and maximum core radius fraction, CRFmin and CRFmax. These hold true even for planets with light volatile envelopes, but require the planet to be fully differentiated and that layers denser than iron are forbidden. With both bounds in hand, a marginal CRF can also be inferred by sampling in-between. After validating on the Earth, we apply our method to Kepler-36b and measure CRFmin = (0.50 +/- 0.07), CRFmax = (0.78 +/- 0.02), and CRFmarg = (0.64 +/- 0.11), broadly consistent with the Earth's true CRF value of 0.55. We apply our method to a suite of hypothetical measurements of synthetic planets to serve as a sensitivity analysis. We find that CRFmin and CRFmax have recovered uncertainties proportional to the relative error on the planetary density, but CRFmarg saturates to between 0.03 and 0.16 once delta rho/rho drops below 1-2 per cent. This implies that mass and radius alone cannot provide any better constraints on internal composition once bulk density constraints hit around a per cent, providing a clear target for observers.

The disintegration of solar filaments via mass drainage is a frequently observed phenomenon during a variety of filament activities. It is generally considered that the draining of dense filament material is directed by both gravity and magnetic field, yet the detailed process remains elusive. Here we report on a partial filament eruption during which filament material drains downward to the surface not only along the filament's legs, but to a remote flare ribbon through a fan-out curtain-like structure. It is found that the magnetic configuration is characterized by two conjoining dome-like quasi-sepratrix layers (QSLs). The filament is located underneath one QSL dome, whose footprint apparently bounds the major flare ribbons resulting from the filament eruption, whereas the remote flare ribbon matches well with the other QSL dome's far-side footprint. We suggest that the interaction of the filament with the overlying QSLs results in the splitting and disintegration of the filament.

In this work, we present a catalog of 2651 carbon stars from the fourth Data Release (DR4) of the Large Sky Area Multi-Object Fiber Spectroscopy Telescope (LAMOST). Using an efficient machine-learning algorithm, we find out these stars from more than seven million spectra. As a by-product, 17 carbon-enhanced metal-poor (CEMP) turnoff star candidates are also reported in this paper, and they are preliminarily identified by their atmospheric parameters. Except for 176 stars that could not be given spectral types, we classify the other 2475 carbon stars into five subtypes including 864 C-H, 226 C-R, 400 C-J, 266 C-N, and 719 barium stars based on a series of spectral features. Furthermore, we divide the C-J stars into three subtypes of CJ( H), C-J(R), C-J(N), and about 90% of them are cool N-type stars as expected from previous literature. Beside spectroscopic classification, we also match these carbon stars to multiple broadband photometries. Using ultraviolet photometry data, we find that 25 carbon stars have FUV detections and they are likely to be in binary systems with compact white dwarf companions.

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.

The Characterizing Exoplanets Satellite (CHEOPS) mission is planned for launch next year with a major objective being to search for transits of known RV planets, particularly those orbiting bright stars. Since the radial velocity method is only sensitive to planetary mass, the radii, transit depths and transit signal-to-noise values of each RV planet are, a-priori, unknown. Using an empirically calibrated probabilistic mass-radius relation, forecaster (Chen & Kipping 2017a), we address this by predicting a catalog of homogeneous credible intervals for these three keys terms for 468 planets discovered via radial velocities. Of these, we find that the vast majority should be detectable with CHEOPS, including terrestrial bodies, if they have the correct geometric alignment. In particular, we predict that 22 mini-Neptunes and 82 Neptune-sized planets would be suitable for detection and that more than 80% of these will have apparent magnitude of V < 10, making them highly suitable for follow-up characterization work. Our work aims to assist the CHEOPS team in scheduling efforts and highlights the great value of quantifiable, statistically robust estimates for upcoming exoplanetary missions.

In this paper, we present the detailed analysis of recurrent homologous jets originating from an emerging negative magnetic flux at the edge of an Active Region. The observed jets show multi-thermal features. Their evolution shows high consistence with the characteristic parameters of the emerging flux, suggesting that with more free magnetic energy, the eruptions tend to be more violent, frequent and blowout-like. The average temperature, average electron number density and axial speed are found to be similar for different jets, indicating that they should have been formed by plasmas from similar origins. Statistical analysis of the jets and their footpoint region conditions reveals a strong positive relationship between the footpoint-region total 131 {\AA} intensity enhancement and jets' length/width. Stronger linearly positive relationships also exist between the total intensity enhancement/thermal energy of the footpoint regions and jets' mass/kinetic/thermal energy, with higher cross-correlation coefficients. All the above results, together, confirm the direct relationship between the magnetic reconnection and the jets, and validate the important role of magnetic reconnection in transporting large amount of free magnetic energy into jets. It is also suggested that there should be more free energy released during the magnetic reconnection of blowout than of standard jet events.

High energy cosmic ray electrons plus positrons (CREs), which lose energy quickly during their propagation, provide an ideal probe of Galactic high-energy processes and may enable the observation of phenomena such as dark-matter particle annihilation or decay. The CRE spectrum has been directly measured up to $\sim 2$ TeV in previous balloon- or space-borne experiments, and indirectly up to $\sim 5$ TeV by ground-based Cherenkov $\gamma$-ray telescope arrays. Evidence for a spectral break in the TeV energy range has been provided by indirect measurements of H.E.S.S., although the results were qualified by sizeable systematic uncertainties. Here we report a direct measurement of CREs in the energy range $25~{\rm GeV}-4.6~{\rm TeV}$ by the DArk Matter Particle Explorer (DAMPE) with unprecedentedly high energy resolution and low background. The majority of the spectrum can be properly fitted by a smoothly broken power-law model rather than a single power-law model. The direct detection of a spectral break at $E \sim0.9$ TeV confirms the evidence found by H.E.S.S., clarifies the behavior of the CRE spectrum at energies above 1 TeV and sheds light on the physical origin of the sub-TeV CREs.

We present results from the investigation of 5-min umbral oscillations in a single-polarity sunspot of active region NOAA 12132. The spectra of TiO, H$\alpha$, and 304 \AA{} are used for corresponding atmospheric heights from the photosphere to lower corona. Power spectrum analysis at the formation height of H$\alpha$ - 0.6 \AA{} to H$\alpha$ center resulted in the detection of 5-min oscillation signals in intensity interpreted as running waves outside the umbral center, mostly with vertical magnetic field inclination $>15\deg$. A phase-speed filter is used to extract the running wave signals with speed $v_{ph}> 4$ km s$^{-1}$, from the time series of H$\alpha$ - 0.4 \AA{} images, and found twenty-four 3-min umbral oscillatory events in a duration of one hour. Interestingly, the initial emergence of the 3-min umbral oscillatory events are noticed closer to or at umbral boundaries. These 3-min umbral oscillatory events are observed for the first time as propagating from a fraction of preceding Running Penumbral Waves (RPWs). These fractional wavefronts rapidly separates from RPWs and move towards umbral center, wherein they expand radially outwards suggesting the beginning of a new umbral oscillatory event. We found that most of these umbral oscillatory events develop further into RPWs. We speculate that the waveguides of running waves are twisted in spiral structures and hence the wavefronts are first seen at high latitudes of umbral boundaries and later at lower latitudes of the umbral center.

In the previous work, a new kind of inflation model was proposed, which has the interesting property that its perturbation equation of motion gets a correction of k^4, due to the non-linearity of the kinetic term. Nonetheless, the scale-invariance of the power spectrum remains valid, both in large-k and small-k limits. In this paper, we investigate in detail the spectral index, the index running and the tensor/scalar ratio in this model, especially on the potential-driven case, and compare the results to the current PLANCK/BICEP observational data. We also discuss the tensor spectrum in this case, which is expected to be tested by the future observations on primordial gravitational waves.

We present the results obtained from the analysis of dynamic fibrils in NOAA active region (AR) 12132, using high resolution H$\alpha$ observations from New Solar Telescope operating at BBSO. The dynamic fibrils are seen to be moving up and down, and most of these dynamic fibrils are periodic and have jet like appearance. We found from our observations that the fibrils follows perfect parabolic paths at the most in many cases. A statistical measure on the properties of the parabolic paths showing an analysis on deceleration, maximum velocity, duration and kinetic energy of these fibrils is presented here. We found the average maximum velocity to be around 15 km s$^{-1}$ and mean deceleration to be around 100 m s$^{-2}$. The deceleration observed appears to be a fraction of gravity of sun and is not compatible with the path of ballistic at the gravity of sun. We found positive correlation between the deceleration and the maximum velocity. This correlation is consistent with the simulations done earlier on magnetoacoustic shock waves propagating upward.

Recent transit surveys have discovered thousands of planetary candidates with directly measured radii, but only a small fraction have measured masses. Planetary mass is crucial in assessing the feasibility of numerous observational signatures, such as radial velocities (RVs), atmospheres, moons and rings. In the absence of a direct measurement, a data-driven, probabilistic forecast enables observational planning and so here we compute posterior distributions for the forecasted mass of approximately seven thousand Kepler Objects of Interest (KOIs). Our forecasts reveal that the predicted RV amplitudes of Neptunian planets are relatively consistent, as a result of transit survey detection bias, hovering around the few m/s level. We find that mass forecasts are unlikely to improve through more precise planetary radii, with the error budget presently dominated by the intrinsic model uncertainty. Our forecasts identify a couple of dozen KOIs near the Terran-Neptunian divide with particularly large RV semi-amplitudes which could be promising targets to follow-up, particularly in the near-IR. With several more transit surveys planned in the near-future, the need to quickly forecast observational signatures is likely to grow and the work here provides a template example of such calculations.

The DArk Matter Particle Explorer (DAMPE), one of the four scientific space science missions within the framework of the Strategic Pioneer Program on Space Science of the Chinese Academy of Sciences, is a general purpose high energy cosmic-ray and gamma-ray observatory, which was successfully launched on December 17th, 2015 from the Jiuquan Satellite Launch Center. The DAMPE scientific objectives include the study of galactic cosmic rays up to $\sim 10$ TeV and hundreds of TeV for electrons/gammas and nuclei respectively, and the search for dark matter signatures in their spectra. In this paper we illustrate the layout of the DAMPE instrument, and discuss the results of beam tests and calibrations performed on ground. Finally we present the expected performance in space and give an overview of the mission key scientific goals.

The scientific yield of current and future optical surveys is increasingly limited by systematic uncertainties in the flux calibration. This is the case for Type Ia supernova (SN Ia) cosmology programs, where an improved calibration directly translates into improved cosmological constraints. Current methodology rests on models of stars. Here we aim to obtain flux calibration that is traceable to state-of-the-art detector-based calibration. We present the SNIFS Calibration Apparatus (SCALA), a color (relative) flux calibration system developed for the SuperNova Integral Field Spectrograph (SNIFS), operating at the University of Hawaii 2.2 m (UH 88) telescope. By comparing the color trend of the illumination generated by SCALA during two commissioning runs, and to previous laboratory measurements, we show that we can determine the light emitted by SCALA with a long-term repeatability better than 1%. We describe the calibration procedure necessary to control for system aging. We present measurements of the SNIFS throughput as estimated by SCALA observations. The SCALA calibration unit is now fully deployed at the UH\,88 telescope, and with it color-calibration between 4000 {\AA} and 9000 {\AA} is stable at the percent level over a one-year baseline.

We studied the background field for 60 two-ribbon flares of M-and-above classes during 2011--2015. These flares are categorized into two groups, i.e., \emph{eruptive} and \emph{confined} flares, based on whether a flare is associated with a coronal mass ejection or not. The background field of source active regions is approximated by a potential field extrapolated from the $B_z$ component of vector magnetograms provided by the Helioseismic and Magnetic Imager. We calculated the decay index $n$ of the background field above the flaring polarity inversion line, and defined a critical height $h_\mathrm{crit}$ corresponding to the theoretical threshold ($n_\mathrm{crit}=1.5$) of the torus instability. We found that $h_\mathrm{crit}$ is approximately half of the distance between the centroids of opposite polarities in active regions, and that the distribution of $h_\mathrm{crit}$ is bimodal: it is significantly higher for confined flares than for eruptive ones. The decay index increases monotonously with increasing height for 86\% (84\%) of the eruptive (confined) flares but displays a saddle-like profile for the rest 14\% (16\%), which are found exclusively in active regions of multipolar field configuration. Moreover, $n$ at the saddle bottom is significantly smaller in confined flares than that in eruptive ones. These results highlight the critical role of background field in regulating the eruptive behavior of two-ribbon flares.

The angular momentum distribution in dark matter haloes and galaxies is a key ingredient in understanding their formation. Especially, the internal distribution of angular momenta is closely related to the formation of disk galaxies. In this article, we use haloes identified from a high-resolution simulation, the Bolshoi simulation, to study the spatial distribution of specific angular momenta, $j(r,\theta)$. We show that by stacking haloes with similar masses to increase the signal-to-noise ratio, the profile can be fitted as a simple function, $j(r,\theta)=j_s \sin^2(\theta/\theta_s) (r/r_s)^2/(1+r/r_s)^4 $, with three free parameters, $j_s, r_s$, and $\theta_s$. Specifically, $j_s$ correlates with the halo mass $M_\mathrm{vir}$ as $j_s\propto M_\mathrm{vir}^{2/3}$, $r_s$ has a weak dependence on the halo mass as $r_s \propto M_\mathrm{vir}^{0.040}$, and $\theta_s$ is independent of $M_\mathrm{vir}$. This profile agrees with that from a rigid shell model, though its origin is unclear. Our universal specific angular momentum profile $j(r,\theta)$ is useful in modelling haloes' angular momenta. Furthermore, by using an empirical stellar mass - halo mass relation, we can infer the averaged angular momentum distribution of a dark matter halo. The specific angular momentum - stellar mass relation within a halo computed from our profile is shown to share a similar shape as that from the observed disk galaxies.

We in this paper investigate the formation and evolution of primordial black holes (PBHs) in nonsingular bouncing cosmologies. We discuss the formation of PBH in the contracting phase and calculate the PBH abundance as a function of the sound speed and Hubble parameter. Afterwards, by taking into account the subsequent PBH evolution during the bouncing phase, we derive the density of PBHs and their Hawking radiation. Our analysis shows that nonsingular bounce models can be constrained from the backreaction of PBHs.

Sky subtraction is the key technique in data reduction of multi-fiber spectra. Knowledge of the related instrument character is necessary to determine the method adopted in sky subtraction. In this study, we described the sky subtraction method designed for LAMOST(Large sky Area Multi-Object fiber Spectroscopic Telescope) survey. The method has been intergrated into LAMOST 2D Pipeline v2.6 and applied to data of LAMOST DR3 and later. For LAMOST, sky emission line calibration is used to alleviate the position-dependent (thus time-dependent) ~4% fiber throughput uncertainty and the small wavelength instability (0.1\AA ) during observation. PCA (Principal Component Analysis) sky subtraction further reduces 25% of the sky line residual of the OH lines in the red part of the LAMOST spectra after the mater sky spectrum, which is derived from a B-spline fit of 20 sky fibers in each spectrograph, is adjusted by sky emission line and subtracted from each fiber. Further analysis shows that our wavelength calibration accuracy is about 4.5km/s, and the average sky subtraction residuals are about 3% for sky emission lines and 3% for continuum region. The relative sky subtraction residuals vary with the moon light background brightness, could reach as low as 1.5% for the sky emission line regions in the dark night. Tests on the F stars of both similar sky emission line strength and similar object continuum intensity show that the sky emission line residual of LAMOST is smaller than those of SDSS survey.

The requirements-driven OSIRIS-REx Camera Suite (OCAMS) acquires images essential to collecting a sample from the surface of Bennu. During proximity operations, these images document the presence of satellites and plumes, record spin state, enable an accurate digital terrain model of the shape of the asteroid and identify any surface hazards. They confirm the presence of sampleable regolith on the surface, observe the sampling event itself, and image the sample head in order to verify its readiness to be stowed. They document the history of Bennu as an example of early solar system material, as a microgravity body with a planetesimal size-scale, and as a carbonaceous object. OCAMS is fitted with three cameras. The MapCam records point-source color images on approach to the asteroid in order to connect ground-based point-source observations of Bennu to later higher-resolution surface spectral imaging. The SamCam documents the sample site before, during, and after it is disturbed by the sample mechanism. The PolyCam, using its focus mechanism, observes the sample site at sub-centimeter resolutions, revealing surface texture and morphology. While their imaging requirements divide naturally between the three cameras, they preserve a strong degree of functional overlap. OCAMS and the other spacecraft instruments allow the OSIRIS-REx mission to collect a sample from a microgravity body on the same visit during which it was first optically acquired from long range, a useful capability as humanity explores near-Earth, Main-Belt and Jupiter Trojan asteroids.

Jets often occur repeatedly from almost the same location. In this paper, a complex solar jet was observed with two phases to the west of NOAA AR 11513 on July 2nd, 2012. If it had been observed at only moderate resolution, the two phases and their points of origin would have been regarded as identical. However, at high resolution we find the two phases merge into one another and the accompanying footpoint brightenings occur at different locations. The phases originate from different magnetic patches rather than being one phase originating from the same patch. Photospheric line of sight (LOS) magnetograms show that the bases of the two phases lie in two different patches of magnetic flux which decrease in size during the occurrence of the two phases. Based on these observations, we suggest the driving mechanism of the two successive phases is magnetic cancellation of two separate magnetic fragments with an opposite polarity fragment between them.

he DArk Matter Particle Explorer (DAMPE) is a general purposed satellite-borne high energy $\gamma-$ray and cosmic ray detector, and among the scientific objectives of DAMPE are the searches for the origin of cosmic rays and an understanding of Dark Matter particles. As one of the four detectors in DAMPE, the Plastic Scintillator Detector (PSD) plays an important role in the particle charge measurement and the photons/electrons separation. The PSD has 82 modules, each consists of a long organic plastic scintillator bar and two PMTs at both ends for readout, in two layers and covers an overall active area larger than 82 cm $\times$ 82 cm. It can identify the charge states for relativistic ions from H to Fe, and the detector efficiency for Z=1 particles can reach 0.9999. The PSD has been successfully launched with DAMPE on Dec. 17, 2015. In this paper, the design, the assembly, the qualification tests of the PSD and some of the performance measured on the ground have been described in detail.

Virtual Observatory (VO) is a data-intensively online astronomical research and education environment, which takes advantages of advanced information technologies to achieve seamless and global access to astronomical information. AstroCloud is a cyber-infrastructure for astronomy research initiated by Chinese Virtual Observatory (China-VO) project, and also a kind of physical distributed platform which integrates lots of tasks such as telescope access proposal management, data archiving, data quality control, data release and open access, cloud based data processing and analysis. It consists of five application channels, i.e. observation, data, tools, cloud and public and is acting as a full lifecycle management system and gateway for astronomical data and telescopes. Physically, the platform is hosted in six cities currently, i.e. Beijing, Nanjing, Shanghai, Kunming, Lijiang and Urumqi, and serving more than 17 thousand users. Achievements from international Virtual Observatories and Cloud Computing are adopted heavily. In the paper, backgrounds of the project, architecture, Cloud Computing environment, key features of the system, current status and future plans are introduced.

Correction of Type Ia Supernova brightnesses for extinction by dust has proven to be a vexing problem. Here we study the dust foreground to the highly reddened SN 2012cu, which is projected onto a dust lane in the galaxy NGC 4772. The analysis is based on multi-epoch, spectrophotometric observations spanning 3,300 - 9,200 {\AA}, obtained by the Nearby Supernova Factory. Phase-matched comparison of the spectroscopically twinned SN 2012cu and SN 2011fe across 10 epochs results in the best-fit color excess of (E(B-V), RMS) = (1.00, 0.03) and total-to-selective extinction ratio of (RV , RMS) = (2.95, 0.08) toward SN 2012cu within its host galaxy. We further identify several diffuse interstellar bands, and compare the 5780 {\AA} band with the dust-to-band ratio for the Milky Way. Overall, we find the foreground dust-extinction properties for SN 2012cu to be consistent with those of the Milky Way. Furthermore we find no evidence for significant time variation in any of these extinction tracers. We also compare the dust extinction curve models of Cardelli et al. (1989), O'Donnell (1994), and Fitzpatrick (1999), and find the predictions of Fitzpatrick (1999) fit SN 2012cu the best. Finally, the distance to NGC4772, the host of SN 2012cu, at a redshift of z = 0.0035, often assigned to the Virgo Southern Extension, is determined to be 16.6$\pm$1.1 Mpc. We compare this result with distance measurements in the literature.

The analysis of Proxima Centauri's radial velocities recently led Anglada-Escud\'e et al. (2016) to claim the presence of a low mass planet orbiting the Sun's nearest star once every 11.2 days. Although the a-priori probability that Proxima b transits its parent star is just 1.5%, the potential impact of such a discovery would be considerable. Independent of recent radial velocity efforts, we observed Proxima Centauri for 12.5 days in 2014 and 31 days in 2015 with the MOST space telescope. We report here that we cannot make a compelling case that Proxima b transits in our precise photometric time series. Imposing an informative prior on the period and phase, we do detect a candidate signal with the expected depth. However, perturbing the phase prior across 100 evenly spaced intervals reveals one strong false-positive and one weaker instance. We estimate a false-positive rate of at least a few percent and a much higher false-negative rate of 20-40%, likely caused by the very high flare rate of Proxima Centauri. Comparing our candidate signal to HATSouth ground-based photometry reveals that the signal is somewhat, but not conclusively, disfavored (1-2 sigmas) leading us to argue that the signal is most likely spurious. We expect that infrared photometric follow-up could more conclusively test the existence of this candidate signal, owing to the suppression of flare activity and the impressive infrared brightness of the parent star.

The task of flux calibration for LAMOST (Large sky Area Multi-Object Spectroscopic Telescope) spectra is difficult due to many factors. For example, the lack of standard stars, flat fielding for large field of view, and variation of reddening between different stars especially at low galactic latitudes etc. Poor selection, bad spectral quality, or extinction uncertainty of standard stars not only might induce errors to the calculated spectral response curve (SRC), but also might lead to failures in producing final 1D spectra. In this paper, we inspected spectra with Galactic latitude |b|>=60 degree and reliable stellar parameters, determined through the LAMOST Stellar Parameter Pipeline (LASP), to study the stability of the spectrograph. To guarantee the selected stars had been observed by each fiber, we selected 37,931 high quality exposures of 29,000 stars from LAMOST DR2, and more than 7 exposures for each fiber. We calculated the SRCs for each fiber for each exposure, and calculated the statistics of SRCs for spectrographs with both the fiber variations and time variations. The result shows that the average response curve of each spectrograph (henceforth ASPSRC) is relatively stable with statistical errors <= 10%. From the comparison between each ASPSRC and the SRCs for the same spectrograph obtained by 2D pipeline, we find that the ASPSRCs are good enough to use for the calibration. The ASPSRCs have been applied to spectra which were abandoned by LAMOST 2D pipeline due to the lack of standard stars, increasing the number of LAMOST spectra by 52,181 in DR2. Comparing those same targets with SDSS, the relative flux differences between SDSS spectra and that of LAMOST spectra with the ASPSRC method are less than 10%, which underlines that the ASPSRC method is feasible for LAMOST flux calibration.

Mass and radius are two of the most fundamental properties of an astronomical object. Increasingly, new planet discoveries are being announced with a measurement of one of these terms, but not both. This has led to a growing need to forecast the missing quantity using the other, especially when predicting the detectability of certain follow-up observations. We present am unbiased forecasting model built upon a probabilistic mass-radius relation conditioned on a sample of 316 well-constrained objects. Our publicly available code, Forecaster, accounts for observational errors, hyper-parameter uncertainties and the intrinsic dispersions observed in the calibration sample. By conditioning our model upon a sample spanning dwarf planets to late-type stars, Forecaster can predict the mass (or radius) from the radius (or mass) for objects covering nine orders-of-magnitude in mass. Classification is naturally performed by our model, which uses four classes we label as Terran worlds, Neptunian worlds, Jovian worlds and stars. Our classification identifies dwarf planets as merely low-mass Terrans (like the Earth), and brown dwarfs as merely high-mass Jovians (like Jupiter). We detect a transition in the mass-radius relation at $2.0_{-0.6}^{+0.7} M_\oplus$, which we associate with the divide between solid, Terran worlds and Neptunian worlds. This independent analysis adds further weight to the emerging consensus that rocky Super-Earths represent a narrower region of parameter space than originally thought. Effectively, then, the Earth is the Super-Earth we have been looking for.

Typical solar flares display two quasi-parallel, bright ribbons on the chromosphere. In between is the polarity inversion line (PIL) separating concentrated magnetic fluxes of opposite polarity in active regions (ARs). Intriguingly a series of flares exhibiting X-shaped ribbons occurred at the similar location on the outskirts of NOAA AR 11967, where magnetic fluxes were scattered, yet three of them were alarmingly energetic. The X shape, whose center coincided with hard X-ray emission, was similar in UV/EUV, which cannot be accommodated in the standard flare model. Mapping out magnetic connectivities in potential fields, we found that the X morphology was dictated by the intersection of two quasi-separatrix layers, i.e., a hyperbolic flux tube (HFT), within which a separator connecting a double null was embedded. This topology was not purely local but regulated by fluxes and flows over the whole AR. The nonlinear force-free field model suggested the formation of a current layer at the HFT, where the current dissipation can be mapped to the X-shaped ribbons via field-aligned heat conduction. These results highlight the critical role of HFTs in 3D magnetic reconnection and have important implications for astrophysical and laboratory plasmas.

The spectrum of relic gravitational wave (RGW) contains high-frequency divergences, which should be removed. We present a systematic study of the issue, based on the exact RGW solution that covers the five stages, from inflation to the acceleration, each being a power law expansion. We show that the present RGW consists of vacuum dominating at $f>10^{11}$Hz and graviton dominating at $f<10^{11}$Hz, respectively. The gravitons are produced by the four cosmic transitions, mostly by the inflation-reheating one. We perform adiabatic regularization to remove vacuum divergences in three schemes: at present, at the end of inflation, and at horizon-exit, to the 2-nd adiabatic order for the spectrum, and the 4-th order for energy density and pressure. In the first scheme a cutoff is needed to remove graviton divergences. We find that all three schemes yield the spectra of a similar profile, and the primordial spectrum defined far outside horizon during inflation is practically unaffected. We also regularize the gauge-invariant perturbed inflaton and the scalar curvature perturbation by the last two schemes, and find that the scalar spectra, the tensor-to-scalar ratio, and the consistency relation remain unchanged.

Solar active regions (ARs) are the major sources of two kinds of the most violent solar eruptions, namely flares and coronal mass ejections (CMEs). The largest AR in the past 24 years, NOAA AR 12192, crossed the visible disk from 2014 October 17 to 30, unusually produced more than one hundred flares, including 32 M-class and 6 X-class ones, but only one small CME. Flares and CMEs are believed to be two phenomena in the same eruptive process. Why is such a flare-rich AR so CME-poor? We compared this AR with other four ARs; two were productive in both and two were inert. The investigation of the photospheric parameters based on the SDO/HMI vector magnetogram reveals that the flare-rich AR 12192, as the other two productive ARs, has larger magnetic flux, current and free magnetic energy than the two inert ARs, but contrast to the two productive ARs, it has no strong, concentrated current helicity along both sides of the flaring neutral line, indicating the absence of a mature magnetic structure consisting of highly sheared or twisted field lines. Furthermore, the decay index above the AR 12192 is relatively low, showing strong constraint. These results suggest that productive ARs are always large and have enough current and free energy to power flares, but whether or not a flare is accompanied by a CME is seemingly related to (1) if there is mature sheared or twisted core field serving as the seed of the CME, (2) if the constraint of the overlying arcades is weak enough.

This paper reports on polarimetric radiation properties based on the switching modes of normal PSR B2020+28 by analysing the data acquired from the Nanshan 25-m radio telescope at 1556 MHz. With nearly 8 hours quasi-continuous observation, the data presented some striking and updated phenomena. The change of relative intensity between the leading and trailing components is the predominant feature of mode switching. The intensity ratio between the leading and trailing components are measured for the individual profiles averaged over 30 seconds. It is found that there is an excess of high ratios over the normal distribution, which indicates that two modes exist in the pulsar. The distribution of abnormal mode has a narrower width indicating that the abnormal mode is more stable than the normal mode. A total of 76 mode switching events are detected in our data. It spends 89\% in the normal mode and 11\% in the abnormal mode. The intrinsic distributions of mode timescales are constrained with power-law distributions. The significant difference in the index of the duration distribution between normal and abnormal modes possibly indicates that the timescale for the abnormal mode to get stable is shorter than that for the normal mode. The frequent switching between both modes may indicate that the oscillations between different magnetospheric states are rapid.

We perform zoom-in cosmological simulations of a suite of dwarf galaxies, examining the impact of cosmic-rays generated by supernovae, including the effect of diffusion. We first look at the effect of varying the uncertain cosmic ray parameters by repeatedly simulating a single galaxy. Then we fix the comic ray model and simulate five dwarf systems with virial masses range from 8-30 $\times 10^{10}$ Msun. We find that including cosmic ray feedback (with diffusion) consistently leads to disk dominated systems with relatively flat rotation curves and constant star formation rates. In contrast, our purely thermal feedback case results in a hot stellar system and bursty star formation. The CR simulations very well match the observed baryonic Tully-Fisher relation, but have a lower gas fraction than in real systems. We also find that the dark matter cores of the CR feedback galaxies are cuspy, while the purely thermal feedback case results in a substantial core.

The 2011 January 28 M1.4 flare exhibits two side-by-side candle-flame-shaped flare loop systems underneath a larger cusp-shaped structure during the decay phase, as observed at the northwestern solar limb by the Solar Dynamics Observatory (SDO). The northern loop system brightens following the initiation of the flare within the southern loop system, but all three cusp-shaped structures are characterized by ~ 10 MK temperatures, hotter than the arch-shaped loops underneath. The "Ahead" satellite of the Solar Terrestrial Relations Observatory (STEREO) provides a top view, in which the post-flare loops brighten sequentially, with one end fixed while the other apparently slipping eastward. By performing stereoscopic reconstruction of the post-flare loops in EUV and mapping out magnetic connectivities, we found that the footpoints of the post-flare loops are slipping along the footprint of a hyperbolic flux tube (HFT) separating the two loop systems, and that the reconstructed loops share similarity with the magnetic field lines that are traced starting from the same HFT footprint, where the field lines are relatively flexible. These results argue strongly in favor of slipping magnetic reconnection at the HFT. The slipping reconnection was likely triggered by the flare and manifested as propagative dimmings before the loop slippage is observed. It may contribute to the late-phase peak in Fe XVI 33.5 nm, which is even higher than its main-phase counterpart, and may also play a role in the density and temperature asymmetry observed in the northern loop system through heat conduction.

Decadal-long radial velocity surveys have recently started to discover analogs to the most influential planet of our solar system, Jupiter. Detecting and characterizing these worlds is expected to shape our understanding of our uniqueness in the cosmos. Despite the great successes of recent transit surveys, Jupiter analogs represent a terra incognita, owing to the strong intrinsic bias of this method against long orbital periods. We here report on the first validated transiting Jupiter analog, Kepler-167e (KOI-490.02), discovered using Kepler archival photometry orbiting the K4-dwarf KIC-3239945. With a radius of $(0.91\pm0.02)$ $R_{\mathrm{Jup}}$, a low orbital eccentricity ($0.06_{-0.04}^{+0.10}$) and an equilibrium temperature of $(131\pm3)$ K, Kepler-167e bears many of the basic hallmarks of Jupiter. Kepler-167e is accompanied by three Super-Earths on compact orbits, which we also validate, leaving a large cavity of transiting worlds around the habitable-zone. With two transits and continuous photometric coverage, we are able to uniquely and precisely measure the orbital period of this post snow-line planet ($1071.2323\pm0.0006$ d), paving the way for follow-up of this $K=11.8$ mag target.

We present well-sampled optical observations of the bright Type Ia supernova (SN~Ia) SN 2011fe in M101. Our data, starting from $\sim16$ days before maximum light and extending to $\sim463$ days after maximum, provide an unprecedented time series of spectra and photometry for a normal SN~Ia. Fitting the early-time rising light curve, we find that the luminosity evolution of SN 2011fe follows a $t^n$ law, with the index $n$ being close to 2.0 in the $VRI$ bands but slightly larger in the $U$ and $B$ bands. Combining the published ultraviolet (UV) and near-infrared (NIR) photometry, we derive the contribution of UV/NIR emission relative to the optical. SN 2011fe is found to have stronger UV emission and reaches its UV peak a few days earlier than other SNe~Ia with similar $\Delta m_{15}(B)$, suggestive of less trapping of high-energy photons in the ejecta. Moreover, the $U$-band light curve shows a notably faster decline at late phases ($t\approx 100$--300 days), which also suggests that the ejecta may be relatively transparent to UV photons. These results favor the notion that SN 2011fe might have a progenitor system with relatively lower metallicity. On the other hand, the early-phase spectra exhibit prominent high-velocity features (HVFs) of O~I $\lambda$7773 and the Ca~II~NIR triplet, but only barely detectable in Si~II~6355. This difference can be caused either by an ionization/temperature effect or an abundance enhancement scenario for the formation of HVFs; it suggests that the photospheric temperature of SN 2011fe is intrinsically low, perhaps owing to incomplete burning during the explosion of the white dwarf.

We present well-sampled optical observations of the bright Type Ia supernova (SN~Ia) SN 2011fe in M101. Our data, starting from $\sim16$ days before maximum light and extending to $\sim463$ days after maximum, provide an unprecedented time series of spectra and photometry for a normal SN~Ia. Fitting the early-time rising light curve, we find that the luminosity evolution of SN 2011fe follows a $t^n$ law, with the index $n$ being close to 2.0 in the $VRI$ bands but slightly larger in the $U$ and $B$ bands. Combining the published ultraviolet (UV) and near-infrared (NIR) photometry, we derive the contribution of UV/NIR emission relative to the optical. SN 2011fe is found to have stronger UV emission and reaches its UV peak a few days earlier than other SNe~Ia with similar $\Delta m_{15}(B)$, suggestive of less trapping of high-energy photons in the ejecta. Moreover, the $U$-band light curve shows a notably faster decline at late phases ($t\approx 100$--300 days), which also suggests that the ejecta may be relatively transparent to UV photons. These results favor the notion that SN 2011fe might have a progenitor system with relatively lower metallicity. On the other hand, the early-phase spectra exhibit prominent high-velocity features (HVFs) of O~I $\lambda$7773 and the Ca~II~NIR triplet, but only barely detectable in Si~II~6355. This difference can be caused either by an ionization/temperature effect or an abundance enhancement scenario for the formation of HVFs; it suggests that the photospheric temperature of SN 2011fe is intrinsically low, perhaps owing to incomplete burning during the explosion of the white dwarf.

As the cyber-infrastructure for Astronomical research from Chinese Virtual Observatory (China-VO) project, AstroCloud has been archived solid progresses during the last one year. Proposal management system and data access system are re-designed. Several new sub-systems are developed, including China-VO PaperData, AstroCloud Statics and Public channel. More data sets and application environments are integrated into the platform. LAMOST DR1, the largest astronomical spectrum archive was released to the public using the platform. The latest progresses will be introduced.

Photometric and spectroscopic observations of a slowly declining, luminous Type Ia supernova (SN Ia) SN 2011hr in the starburst galaxy NGC 2691 are presented. SN 2011hr is found to peak at $M_{B}=-19.84 \pm 0.40\,\rm{mag}$, with a post-maximum decline rate $\Delta$m$_{15}$(B) = 0.92 $\pm$ 0.03\,$\rm{mag}$. From the maximum-light bolometric luminosity, $L=(2.30 \pm 0.90) \times 10^{43}\,\rm{erg\,s^{-1}}$, we estimate the mass of synthesized \Nifs\ in SN 2011hr to be $M(\rm{^{56}Ni})=1.11 \pm 0.43\,M_{\sun}$. SN 2011hr appears more luminous than SN 1991T at around maximum light, and the absorption features from its intermediate-mass elements (IMEs) are noticeably weaker than the latter at similar phases. Spectral modeling suggests that SN 2011hr has the IMEs of $\sim$\,0.07 M$_{\sun}$ in the outer ejecta, which is much lower than the typical value of normal SNe Ia (i.e., 0.3 -- 0.4 M$_{\sun}$) and is also lower than the value of SN 1991T (i.e., $\sim$\,0.18 M$_{\sun}$). These results indicate that SN 2011hr may arise from a Chandrasekhar-mass white dwarf progenitor that experienced a more efficient burning process in the explosion. Nevertheless, it is still possible that SN 2011hr may serve as a transitional object connecting the SN 1991T-like SNe Ia with the superluminous subclass like SN 2007if given that the latter also shows very weak IMEs at all phases.

We investigate the evolution of NOAA Active Region 11817 during 2013 August 10--12, when it developed a complex field configuration and produced four confined, followed by two eruptive, flares. These C-and-above flares are all associated with a magnetic flux rope (MFR) located along the major polarity inversion line, where shearing and converging photospheric flows are present. Aided by the nonlinear force-free field modeling, we identify the MFR through mapping magnetic connectivities and computing the twist number $\mathcal{T}_w$ for each individual field line. The MFR is moderately twisted ($|\mathcal{T}_w| < 2$) and has a well-defined boundary of high squashing factor $Q$. We found that the field line with the extremum $|\mathcal{T}_w|$ is a reliable proxy of the rope axis, and that the MFR's peak $|\mathcal{T}_w|$ temporarily increases within half an hour before each flare while it decreases after the flare peak for both confined and eruptive flares. This pre-flare increase in $|\mathcal{T}_w|$ has little effect on the active region's free magnetic energy or any other parameters derived for the whole region, due to its moderate amount and the MFR's relatively small volume, while its decrease after flares is clearly associated with the stepwise decrease in free magnetic energy due to the flare. We suggest that $\mathcal{T}_w$ may serve as a useful parameter in forewarning the onset of eruption, and therefore, the consequent space weather effects. The helical kink instability is identified as the prime candidate onset mechanism for the considered flares.

Two major coronal mass ejections (CMEs) observed during the Whole Heliosphere Interval (WHI) are compared with the catastrophe (CA) and eruptive flux rope (EF) models. The objective is to test two distinct mechanisms for CMEs by modeling these well-observed CMEs and comparing predictions of the theories and observed data. The two CMEs selected for this study occurred on 25 March and 5 April 2008, respectively. For the 25 March event, an M 1.7 class flare, a filament eruption, and hard X-ray (HXR) and soft X-ray (SXR) emissions were observed during the CME onset. The observed CME kinematics and SXR light curve of this event are found to be more consistent with the EF model than with the CA model. For the 5 April event, the SXR light curve shows multiple enhancements, some of which temporally coincide with successive side loop brightening and multiple foot points at the source region after the eruption. The physical connection between the side-loop multiple brightenings and the eruption cannot be determined from the data. Both models produced observationally consistent kinematics profiles, and the EF model correctly predicted the first emission enhancement. Neither model includes multiple brightenings in the formulation.

For an unattended telescopes in Antarctic, the remote operation, autonomous observation and control are essential. An EPICS (Experimental Physics and Industrial Control System) and RTS2(Remote Telescope System, 2nd Version) based autonomous observation and control system with remoted operation is introduced in this paper. EPICS is a set of Open Source software tools, libraries and applications developed collaboratively and used worldwide to create distributed soft real-time control systems for scientific instruments while RTS2 is an open source environment for control of a fully autonomous observatory. Using the advantage of EPICS and RTS2 respectively, a combined integrated software framework for autonomous observation and control is established that use RTS2 to fulfill the function of astronomical observation and use EPICS to fulfill the device control of telescope. A command and status interface for EPICS and RTS2 is designed to make the EPICS IOC (Input/Output Controller) components integrate to RTS2 directly. For the specification and requirement of control system of telescope in Antarctic, core components named Executor and Auto-focus for autonomous observation is designed and implemented with remote operation user interface based on Browser-Server mode. The whole system including the telescope is tested in Lijiang Observatory in Yunnan Province for practical observation to complete the autonomous observation and control, including telescope control, camera control, dome control, weather information acquisition with the local and remote operation.

We introduce a method for identifying "twin" Type Ia supernovae, and using them to improve distance measurements. This novel approach to Type Ia supernova standardization is made possible by spectrophotometric time series observations from the Nearby Supernova Factory (SNfactory). We begin with a well-measured set of supernovae, find pairs whose spectra match well across the entire optical window, and then test whether this leads to a smaller dispersion in their absolute brightnesses. This analysis is completed in a blinded fashion, ensuring that decisions made in implementing the method do not inadvertently bias the result. We find that pairs of supernovae with more closely matched spectra indeed have reduced brightness dispersion. We are able to standardize this initial set of SNfactory supernovae to 0.083 +/- 0.012 magnitudes, implying a dispersion of 0.072 +/- 0.010 magnitudes in the absence of peculiar velocities. We estimate that with larger numbers of comparison SNe, e.g, using the final SNfactory spectrophotometric dataset as a reference, this method will be capable of standardizing high-redshift supernovae to within 0.06-0.07 magnitudes. These results imply that at least 3/4 of the variance in Hubble residuals in current supernova cosmology analyses is due to previously unaccounted-for astrophysical differences among the supernovae

In this paper, we report the detections of stellar variabilities from the first 2-year observations of sky area of about 1300 square degrees from the Tsinghua University-NAOC Transient Survey (TNTS). A total of 1237 variable stars (including 299 new ones) were detected with brightness < 18.0 mag and magnitude variation >= 0.1 mag on a timescale from a few hours to few hundred days. Among such detections, we tentatively identified 661 RR Lyrae stars, 431 binaries, 72 Semiregular pulsators, 29 Mira stars, 11 slow irregular variables, 11 RS Canum Venaticorum stars, 7 Gamma Doradus stars, 5 long period variables, 3 W Virginis stars, 3 Delta Scuti stars, 2 Anomalous Cepheids, 1 Cepheid, and 1 nove-like star based on their time-series variability index Js and their phased diagrams. Moreover, we found that 14 RR Lyrae stars show the Blazhko effect and 67 contact eclipsing binaries exhibit the O'Connell effect. Since the period and amplitude of light variations of RR Lyrae variables depend on their chemical compositions, their photometric observations can be used to investigate distribution of metallicity along the direction perpendicular to the Galactic disk. We find that the metallicity of RR Lyrae stars shows large scatter at regions closer to the Galactic plane (e.g., -3.0 < [Fe/H] < 0) but tends to converge at [Fe/H]~ -1.7 at larger Galactic latitudes. This variation may be related to that the RRAB Lyrae stars in the Galactic halo come from globular clusters with different metallicity and vertical distances, i.e. OoI and OoII populations, favoring for the dual-halo model.

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

We present extensive ultraviolet, optical, and near-infrared observations of the type IIP supernova (SN IIP) 2013ej in the nearby spiral galaxy M74. The multicolor light curves, spanning from $\sim$ 8--185 days after explosion, show that it has a higher peak luminosity (i.e., M$_{V}$ $\sim$$-$17.83 mag at maximum light), a faster post-peak decline, and a shorter plateau phase (i.e., $\sim$ 50 days) compared to the normal type IIP SN 1999em. The mass of $^{56}$Ni is estimated as 0.02$\pm$0.01 M$_{\odot}$ from the radioactive tail of the bolometric light curve. The spectral evolution of SN 2013ej is similar to that of SN 2004et and SN 2007od, but shows a larger expansion velocity (i.e., $v_{Fe II} \sim$ 4600 km s$^{-1}$ at t $\sim$ 50 days) and broader line profiles. In the nebular phase, the emission of H$\alpha$ line displays a double-peak structure, perhaps due to the asymmetric distribution of $^{56}$Ni produced in the explosion. With the constraints from the main observables such as bolometric light curve, expansion velocity and photospheric temperature of SN 2013ej, we performed hydrodynamical simulations of the explosion parameters, yielding the total explosion energy as $\sim$0.7$\times$ 10$^{51}$ erg, the radius of the progenitor as $\sim$600 R$_{\odot}$, and the ejected mass as $\sim$10.6 M$_{\odot}$. These results suggest that SN 2013ej likely arose from a red supergiant with a mass of 12--13 M$_{\odot}$ immediately before the explosion.

Previously we used the Nearby Supernova Factory sample to show that SNe~Ia having locally star-forming environments are dimmer than SNe~Ia having locally passive environments.Here we use the \constitution\ sample together with host galaxy data from \GALEX\ to independently confirm that result. The effect is seen using both the SALT2 and MLCS2k2 lightcurve fitting and standardization methods, with brightness differences of $0.094 \pm 0.037\ \mathrm{mag}$ for SALT2 and $0.155 \pm 0.041\ \mathrm{mag}$ for MLCS2k2 with $R_V=2.5$. When combined with our previous measurement the effect is $0.094 \pm 0.025\ \mathrm{mag}$ for SALT2. If the ratio of these local SN~Ia environments changes with redshift or sample selection, this can lead to a bias in cosmological measurements. We explore this issue further, using as an example the direct measurement of $H_0$. \GALEX{} observations show that the SNe~Ia having standardized absolute magnitudes calibrated via the Cepheid period--luminosity relation using {\textit{HST}} originate in predominately star-forming environments, whereas only ~50% of the Hubble-flow comparison sample have locally star-forming environments. As a consequence, the $H_0$ measurement using SNe~Ia is currently overestimated. Correcting for this bias, we find a value of $H_0^{corr}=70.6\pm 2.6\ \mathrm{km\ s^{-1}\ Mpc^{-1}}$ when using the LMC distance, Milky Way parallaxes and the NGC~4258 megamaser as the Cepheid zeropoint, and $68.8\pm 3.3\ \mathrm{km\ s^{-1}\ Mpc^{-1}}$ when only using NGC~4258. Our correction brings the direct measurement of $H_0$ within $\sim 1\,\sigma$ of recent indirect measurements based on the CMB power spectrum.

It is suggested by Chen {\it et al.} that the Earth's surface Ultraviolet irradiance ($280-400$ nm) could influence the Earth's surface temperature variation by "Highly Excited Water Vapor" (HEWV) effect. In this manuscript, we reconstruct the developing history of the ozone layer depth variation from 1860 to 2011 based on the HEWV effect. It is shown that the reconstructed ozone layer depth variation correlates with the observational variation from 1958 to 2005 very well ($R=0.8422$, $P>99.9\%$). From this reconstruction, we may limit the spectra band of the surface Ultraviolet irradiance referred in HEWV effect to Ultraviolet B ($280-320$ nm).

The frequency dependence of pulse width is studied for 150 normal pulsars, mostly selected from the European Pulsar Network, for which the 10\% multifrequency pulse widths can be well fit with the Thorsett relationship W_{10}=A\nu^{\mu}+W_{\rm 10, min}. The relative fraction of pulse width change between 0.4~GHz and 4.85~GHz, \eta=(W_{4.85}-W_{0.4})/W_{0.4}, is calculated in terms of the best-fit relationship for each pulsar. It is found that 81 pulsars (54\%) have $\eta<-10$\% (group A), showing considerable profile narrowing at high frequencies, 40 pulsars (27\%) have $-$10\%$\leq\eta\leq 10$\% (group B), meaning a marginal change in pulse width, and 29 pulsars (19\%) have $\eta>10$\% (group C), showing a remarkable profile broadening at high frequencies. The fractions of the group-A and group-C pulsars suggest that the profile narrowing phenomenon at high frequencies is more common than the profile broadening phenomenon, but a large fraction of the group-B and group-C pulsars (a total of 46\%) is also revealed. The group-C pulsars, together with a portion of group-B pulsars with a slight pulse broadening, can hardly be explained using the conventional radius-to-frequency mapping, which only applies to the profile narrowing phenomenon. Based on a recent version of the fan beam model, a type of broadband emission model, we propose that the diverse frequency dependence of pulse width is a consequence of different types of distribution of emission spectra across the emission region. The geometrical effect predicting a link between the emission beam shrinkage and the spectrum steepening is tested but disfavored.

AstroCloud is a cyber-Infrastructure for Astronomy Research initiated by Chinese Virtual Observatory (China-VO) under funding support from NDRC (National Development and Reform commission) and CAS (Chinese Academy of Sciences){\url{http://astrocloud.china-vo.org}}\citep{O8-5_Cui_adassxxiv}. To archive the astronomical data in China, we present the implementation of the astronomical data archiving system (ADAS). Data archiving and quality control are the infrastructure for the AstroCloud. Throughout the data of the entire life cycle, data archiving system standardized data, transferring data, logging observational data, archiving ambient data, And storing these data and metadata in database. Quality control covers the whole process and all aspects of data archiving.

We develop a new framework for use in exploring Type Ia Supernova (SN Ia) spectra. Combining Principal Component Analysis (PCA) and Partial Least Square analysis (PLS) we are able to establish correlations between the Principal Components (PCs) and spectroscopic/photometric SNe Ia features. The technique was applied to ~120 supernova and ~800 spectra from the Nearby Supernova Factory. The ability of PCA to group together SNe Ia with similar spectral features, already explored in previous studies, is greatly enhanced by two important modifications: (1) the initial data matrix is built using derivatives of spectra over the wavelength, which increases the weight of weak lines and discards extinction, and (2) we extract time evolution information through the use of entire spectral sequences concatenated in each line of the input data matrix. These allow us to define a stable PC parameter space which can be used to characterize synthetic SN Ia spectra by means of real SN features. Using PLS, we demonstrate that the information from important previously known spectral indicators (namely the pseudo-equivalent width (pEW) of Si II 5972 / Si II 6355 and the line velocity of S II 5640 / Si II 6355) at a given epoch, is contained within the PC space and can be determined through a linear combination of the most important PCs. We also show that the PC space encompasses photometric features like B or V magnitudes, B-V color and SALT2 parameters c and x1. The observed colors and magnitudes, that are heavily affected by extinction, cannot be reconstructed using this technique alone. All the above mentioned applications allowed us to construct a metric space for comparing synthetic SN Ia spectra with observations.

The Tsinghua University-National Astronomical Observatories of China (NAOC) Transient Survey (TNTS) is an automatic survey for a systematic exploration of optical transients (OTs), conducted with a 60/90 cm Schmidt telescope at Xinglong station of NAOC. This survey repeatedly covers ~ 1000 square degrees of the north sky with a cadence of 3-4 days. With an exposure of 60 s, the survey reaches a limited unfiltered magnitude of about 19.5 mag. This enables us to discover supernovae at their relatively young stages. In this paper, we describe the overall performance of our survey during the first year and present some preliminary results.

Whether natural factors could interpret the rise of the Earth's surface temperature is still controversial. Though numerous recent researches have reported apparent correlations between solar activity and the Earth's climate, solar activity has encountered a big problem when describing the rapid global warming after 1970s. Our investigation shows the good positive correlations between the Earth's surface Ultraviolet irradiance (280-400 nm) and the Earth's surface temperature both in temporal and spatial variations by analyzing the global surface Ultraviolet irradiance (280-400 nm) and global surface temperature data from 1980-1999. The rise of CO$_2$ cannot interpret the good positive correlations, and we could even get an opposite result to the good correlations when employing the rise of CO$_2$ to describe the relation between them. Based on the good positive correlations, we suggest a new effect, named "Highly Excited Water Vapor" (HEWV) effect, which can interpret how the Sun influences the Earth's surface temperature reasonably, including the rapid warming after 1970s.

We report observations of molecular oxygen (O$_2$) rotational transitions at 487 GHz, 774 GHz, and 1121 GHz toward Orion Peak A. The O2 lines at 487 GHz and 774 GHz are detected at velocities of 10-12 km/s with line widths 3 km/s; however, the transition at 1121 GHz is not detected. The observed line characteristics, combined with the results of earlier observations, suggest that the region responsible for the O$_2$ emission is 9" (6e16 cm) in size, and is located close to the H2 Peak 1position (where vibrationally-excited H$_2$ emission peaks), and not at Peak A, 23" away. The peak O2 column density is 1.1e18/cm2. The line velocity is close to that of 621 GHz water maser emission found in this portion of the Orion Molecular Cloud, and having a shock with velocity vector lying nearly in the plane of the sky is consistent with producing maximum maser gain along the line-of-sight. The enhanced O$_2$ abundance compared to that generally found in dense interstellar clouds can be explained by passage of a low-velocity C-shock through a clump with preshock density 2e4/cm3, if a reasonable flux of UV radiation is present. The postshock O$_2$ can explain the emission from the source if its line of sight dimension is ~10 times larger than its size on the plane of the sky. The special geometry and conditions required may explain why O$_2$ emission has not been detected in the cores of other massive star-forming molecular clouds.

We present extensive optical observations of the normal Type Ic supernova (SN) 2007gr, spanning from about one week before maximum light to more than one year thereafter. The optical light and color curves of SN 2007gr are very similar to those of the broad-lined Type Ic SN 2002ap, but the spectra show remarkable differences. The optical spectra of SN 2007gr are characterized by unusually narrow lines, prominent carbon lines, and slow evolution of the line velocity after maximum light. The earliest spectrum (taken at t=-8 days) shows a possible signature of helium (He~I 5876 at a velocity of ~19,000 km s{-1}). Moreover, the larger intensity ratio of the [O I] 6300 and 6364 lines inferred from the early nebular spectra implies a lower opacity of the ejecta shortly after the explosion. These results indicate that SN 2007gr perhaps underwent a less energetic explosion of a smaller-mass Wolf-Rayet star (~ 8--9 Msun) in a binary system, as favored by an analysis of the progenitor environment through pre-explosion and post-explosion Hubble Space Telescope images. In the nebular spectra, asymmetric double-peaked profiles can be seen in the [O~I] 6300 and Mg~I] 4571 lines. We suggest that the two peaks are contributed by the blueshifted and rest-frame components. The similarity in velocity structure and the different evolution of the strength of the two components favor an aspherical explosion with the ejecta distributed in a torus or disk-like geometry, but inside the ejecta the O and Mg have different distributions.

Precise fiber positioning is crucial to a wide field, multi-fiber spectroscopic survey like LAMOST. Nowadays, most position error measurements are based on CCD photographic and imaging processing techniques. Those methods only work for measuring errors orthogonal to the telescope optical axis, while there also lies errors parallel to the telescope optical axis, like defocusing, and error caused by the existing deviation angle between optical axes of a fiber and the telescope. Directly measuring two latter types of position errors is difficult for individual fiber, especially during observation. Possible sources of fiber position errors are discussed in brief for LAMOST. By constructing a model of magnitude loss due to the fiber position error for a point source, we propose an indirect method to calculate the total and systematic position errors for each individual fiber from spectra data. Restrictions and applications of this method are also discussed.

We present optical photometric and spectroscopic observations of the Type Ibn (SN 2006jc-like) supernova iPTF13beo. Detected by the intermediate Palomar Transient Factory ~3 hours after the estimated first light, iPTF13beo is the youngest and the most distant (~430 Mpc) Type Ibn event ever observed. The iPTF13beo light curve is consistent with light curves of other Type Ibn SNe and with light curves of fast Type Ic events, but with a slightly faster rise-time of two days. In addition, the iPTF13beo R-band light curve exhibits a double-peak structure separated by ~9 days, not observed before in any Type Ibn SN. A low-resolution spectrum taken during the iPTF13beo rising stage is featureless, while a late-time spectrum obtained during the declining stage exhibits narrow and intermediate-width He I and Si II features with FWHM ~ 2000-5000 km/s and is remarkably similar to the prototypical SN Ibn 2006jc spectrum. We suggest that our observations support a model of a massive star exploding in a dense He-rich circumstellar medium (CSM). A shock breakout in a CSM model requires an eruption releasing a total mass of ~0.1 Msun over a time scale of couple of weeks prior to the SN explosion.

A sensitive broadband molecular line survey of the Sagittarius B2(N) star-forming region has been obtained with the HIFI instrument on the Herschel Space Observatory, offering the first high-spectral resolution look at this well-studied source in a wavelength region largely inaccessible from the ground (625-157 um). From the roughly 8,000 spectral features in the survey, a total of 72 isotopologues arising from 44 different molecules have been identified, ranging from light hydrides to complex organics, and arising from a variety of environments from cold and diffuse to hot and dense gas. We present an LTE model to the spectral signatures of each molecule, constraining the source sizes for hot core species with complementary SMA interferometric observations, and assuming that molecules with related functional group composition are cospatial. For each molecule, a single model is given to fit all of the emission and absorption features of that species across the entire 480-1910 GHz spectral range, accounting for multiple temperature and velocity components when needed to describe the spectrum. As with other HIFI surveys toward massive star forming regions, methanol is found to contribute more integrated line intensity to the spectrum than any other species. We discuss the molecular abundances derived for the hot core, where the local thermodynamic equilibrium approximation is generally found to describe the spectrum well, in comparison to abundances derived for the same molecules in the Orion KL region from a similar HIFI survey.

Extensive optical and ultraviolet (UV) observations of the type Ia supernova (SN Ia) 2012fr are presented in this paper. It has a relatively high luminosity, with an absolute $B$-band peak magnitude of about $-19.5$ mag and a smaller post-maximum decline rate than normal SNe Ia [e.g., $\Delta m _{15}$($B$) $= 0.85 \pm 0.05$ mag]. Based on the UV and optical light curves, we derived that a $^{56}$Ni mass of about 0.88 solar masses was synthesized in the explosion. The earlier spectra are characterized by noticeable high-velocity features of \ion{Si}{2} $\lambda$6355 and \ion{Ca}{2} with velocities in the range of $\sim22,000$--$25,000$ km s$^{-1}$. At around the maximum light, these spectral features are dominated by the photospheric components which are noticeably narrower than normal SNe Ia. The post-maximum velocity of the photosphere remains almost constant at $\sim$12,000 km s$^{-1}$ for about one month, reminiscent of the behavior of some luminous SNe Ia like SN 1991T. We propose that SN 2012fr may represent a subset of the SN 1991T-like SNe Ia viewed in a direction with a clumpy or shell-like structure of ejecta, in terms of a significant level of polarization reported in Maund et al. (2013).

Is it possible that the current cosmic accelerating expansion will turn into a decelerating one? Can this transition be realized by some viable theoretical model that is consistent with the standard Big Bang cosmology? We study a class of phenomenological models of a transient acceleration, based on a dynamical dark energy with a very general form of equation of state $p_{de}=\alpha\rho_{de}-\beta\rho_{de}^m$. It mimics the cosmological constant $\rho_{de}\rightarrow$ const for small scale factor $a$, and behaves as a barotropic gas with $\rho_{de}\rightarrow a^{-3(\alpha+1)}$ with $\alpha\ge 0$ for large $a$. The cosmic evolution of four models in the class has been examined in details, and all yields a smooth transient acceleration. Depending on the specific model, the future universe may be dominated either by the dark energy or by the matter. In two models, the dynamical dark energy can be explicitly realized by a scalar field with an analytical potential $V(\phi)$. Moreover, the statistical analysis shows that the models can be as robust as $\Lambda$CDM in confronting the observational data of SN Ia, CMB, and BAO. As improvements over the previous studies, our models overcome the over-abundance problem of dark energy during early eras, and satisfy the constraints on the dark energy from WMAP observations of CMB.

According to traditional gas-phase chemical models, O2 should be abundant in molecular clouds, but until recently, attempts to detect interstellar O2 line emission with ground- and space-based observatories have failed. Following the multi-line detections of O2 with low abundances in the Orion and rho Oph A molecular clouds with Herschel, it is important to investigate other environments, and we here quantify the O2 abundance near a solar-mass protostar. Observations of O2, at 487 GHz toward a deeply embedded low-mass Class 0 protostar, NGC 1333-IRAS 4A, are presented, using the HIFI instrument on the Herschel Space Observatory. Complementary data of the chemically related NO and CO molecules are obtained as well. The high spectral resolution data are analysed using radiative transfer models to infer column densities and abundances, and are tested directly against full gas-grain chemical models. The deep HIFI spectrum fails to show O2 at the velocity of the dense protostellar envelope, implying one of the lowest abundance upper limits of O2/H2 at <6x10^-9 (3 sigma). However, a tentative (4.5 sigma) detection of O2 is seen at the velocity of the surrounding NGC 1333 molecular cloud, shifted by 1 km/s relative to the protostar. For the protostellar envelope, pure gas-phase models and gas-grain chemical models require a long pre-collapse phase (~0.7-1x10^6 years), during which atomic and molecular oxygen are frozen out onto dust grains and fully converted to H2O, to avoid overproduction of O2 in the dense envelope. The same model also reproduces the limits on the chemically related NO molecule. The tentative detection of O2 in the surrounding cloud is consistent with a low-density PDR model with small changes in reaction rates. The low O2 abundance in the collapsing envelope around a low-mass protostar suggests that the gas and ice entering protoplanetary disks is very poor in O2.

We continue our systematic statistical study of various components in gamma-ray burst (GRB) optical lightcurves. We decompose the early onset bump and the late re-brightening bump with empirical fits and analyze their statistical properties. Among the 146 GRBs that have well-sampled optical lightcurves, the onset and re-brightening bumps are observed in 38 and 26 GRBs, respectively. It is found that the typical rising and decaying slopes for both the onset and re-brightening bumps are ~1.5 and -1.15, respectively. No early onset bumps in the X-ray band are detected to be associated with the optical onset bumps, while an X-ray re-brightening bump is detected for half of the re-brightening optical bumps. The peak luminosity is anti-correlated with the peak time, L_p\propto t_{p}^{-1.81+/-0.32} and L_p\propto t_{p}^{-0.83+/-0.17} for the onset and re-brightening bumps, respectively. Both L_p and the isotropic energy release of the onset bumps are correlated with E_{\gamma, iso}, whereas no similar correlation is found for the re-brightening bumps. Taking the onset bumps as probes for the properties of the fireballs and their ambient medium, we find that the typical power-law index of the relativistic electrons is 2.5 and the medium density profile behaves as n\propto r^{-1} within the framework of the synchrotron external shock models. With the medium density profile obtained from our analysis, we also confirm the correlation between initial Lorentz factor (\Gamma_0) and E_{\gamma, iso} in our previous work. The jet component that produces the re-brightening bump seems to be on-axis and independent of the prompt emission jet component. Its typical kinetic energy budget would be about one order of magnitude larger than the prompt emission component, but with a lower \Gamma_0, typically several tens.

We present Herschel-HIFI, SPIRE, and PACS 50-670 {\mu}m imaging and spectroscopy of six FU Orionis-type objects and candidates (FU Orionis, V1735 Cyg, V1515 Cyg, V1057 Cyg, V1331 Cyg, and HBC 722), ranging in outburst date from 1936-2010, from the "FOOSH" (FU Orionis Objects Surveyed with Herschel) program, as well as ancillary results from Spitzer-IRS and the Caltech Submillimeter Observatory. In their system properties (Lbol, Tbol, line emission), we find that FUors are in a variety of evolutionary states. Additionally, some FUors have features of both Class I and II sources: warm continuum consistent with Class II sources, but rotational line emission typical of Class I, far higher than Class II sources of similar mass/luminosity. Combining several classification techniques, we find an evolutionary sequence consistent with previous mid-IR indicators. We detect [O I] in every source at luminosities consistent with Class 0/I protostars, much greater than in Class II disks. We detect transitions of 13CO (J_up of 5 to 8) around two sources (V1735 Cyg and HBC 722) but attribute them to nearby protostars. Of the remaining sources, three (FU Ori, V1515 Cyg, and V1331 Cyg) exhibit only low-lying CO, but one (V1057 Cyg) shows CO up to J = 23 - 22 and evidence for H2O and OH emission, at strengths typical of protostars rather than T Tauri stars. Rotational temperatures for "cool" CO components range from 20-81 K, for ~ 10^50 total CO molecules. We detect [C I] and [N II] primarily as diffuse emission.

The CXOCY J220132.8-320144 system consists of an edge-on spiral galaxy lensing a background quasar into two bright images. Previous efforts to constrain the mass distribution in the galaxy have suggested that at least one additional image must be present (Castander et al. 2006). These extra images may be hidden behind the disk which features a prominent dust lane. We present and analyze Hubble Space Telescope (HST) observations of the system. We do not detect any extra images, but the observations further narrow the observable parameters of the lens system. We explore a range of models to describe the mass distribution in the system and find that a variety of acceptable model fits exist. All plausible models require 2 magnitudes of dust extinction in order to obscure extra images from detection, and some models may require an offset between the center of the galaxy and the center of the dark matter halo of 1 kiloparsec. Currently unobserved images will be detectable by future James Webb Space Telescope (JWST) observations and will provide strict constraints on the fraction of mass in the disk.

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

We present 32 epochs of optical (3300-9700 \AA) spectrophotometric observations of the nearby quintessential "normal" type Ia supernova (SN Ia) SN 2011fe in the galaxy M101, extending from -15 to +97 d with respect to B-band maximum, obtained by the Nearby Supernova Factory collaboration. SN 2011fe is the closest (\mu = 29.04) and brightest (Bmax = 9.94 mag) SN Ia observed since the advent of modern large scale programs for the intensive periodic followup of supernovae. Both synthetic light curve measurements and spectral feature analysis attest to the normality of SN 2011fe. There is very little evidence for reddening in its host galaxy. The homogeneous calibration, intensive time sampling, and high signal-to-noise ratio of the data set make it unique. Thus it is ideal for studying the physics of SN Ia explosions in detail, and for furthering the use of SNe Ia as standardizable candles for cosmology. Several such applications are shown, from the creation of a bolometric light curve and measurement of the 56Ni mass, to the simulation of detection thresholds for unburned carbon, direct comparisons with other SNe Ia, and existing spectral templates.

This work is to study the dynamics of coronal mass ejections and to investigate possible correlation between CME productivity and subsurface structure. Two CMEs and six active regions are selected for the study. The CMEs are examined by comparing observations and theoretical models, and the subsurface structures are probed by local helioseismic inversions. The analysis of the CMEs shows that the eruptive flux-rope model is in good agreement with both events. However, some discrepancies with the observation are also found, indicating that the model can be further improved. The helioseismic investigation results indicate a consistent correlation between the CME productivity and the subsurface temperature structure. The inferred subsurface magnetic structure reveals that the source regions of the two studied CMEs may share similar subsurface structures. Not found in the other CME productive regions selected for consideration, however, the similarity in the subsurface features may indicate the presence of a subsurface structural connection betweenthe two closely located and CME-productive regions.

We study the constraints on the cosmic opacity using the latest BAO and Union2 SNIa data in this paper and find that the best fit values seem to indicate that an opaque universe is preferred in redshift regions $0.20-0.35$, $0.35-0.44$ and $0.60-0.73$, whereas, a transparent universe is favored in redshift regions $0.106-0.20$, $0.44-0.57$and $0.57-0.60$. However, our result is still consistent with a transparent universe at the 1$\sigma$ confidence level, even though the best-fit cosmic opacity oscillates between zero and some nonzero values as the redshift varies.

We present extensive optical observations of a Type IIn supernova (SN) 2010jl for the first 1.5 years after the discovery. The UBVRI light curves demonstrated an interesting two-stage evolution during the nebular phase, which almost flatten out after about 90 days from the optical maximum. SN 2010jl has one of the highest intrinsic H_alpha luminosity ever recorded for a SN IIn, especially at late phase, suggesting a strong interaction of SN ejecta with the dense circumstellar material (CSM) ejected by the progenitor. This is also indicated by the remarkably strong Balmer lines persisting in the optical spectra. One interesting spectral evolution about SN 2010jl is the appearance of asymmetry of the Balmer lines. These lines can be well decomposed into a narrow component and an intermediate-width component. The intermediate-width component showed a steady increase in both strength and blueshift with time until t ~ 400 days after maximum, but it became less blueshifted at t ~ 500 days when the line profile appeared relatively symmetric again. Owing to that a pure reddening effect will lead to a sudden decline of the light curves and a progressive blueshift of the spectral lines, we therefore propose that the asymmetric profiles of H lines seen in SN 2010jl is unlikely due to the extinction by newly formed dust inside the ejecta, contrary to the explanation by some early studies. Based on a simple CSM-interaction model, we speculate that the progenitor of SN 2010jl may suffer a gigantic mass loss (~ 30-50 M_sun) in a few decades before explosion. Considering a slow moving stellar wind (e.g., ~ 28 km/s) inferred for the preexisting, dense CSM shell and the extremely high mass-loss rate (1-2 M_sun per yr), we suggest that the progenitor of SN 2010jl might have experienced a red supergiant stage and explode finally as a post-red supergiant star with an initial mass above 30-40 M_sun.

We present extensive optical observations of a Type IIn supernova (SN) 2010jl for the first 1.5 years after the discovery. The UBVRI light curves demonstrated an interesting two-stage evolution during the nebular phase, which almost flatten out after about 90 days from the optical maximum. SN 2010jl has one of the highest intrinsic H_alpha luminosity ever recorded for a SN IIn, especially at late phase, suggesting a strong interaction of SN ejecta with the dense circumstellar material (CSM) ejected by the progenitor. This is also indicated by the remarkably strong Balmer lines persisting in the optical spectra. One interesting spectral evolution about SN 2010jl is the appearance of asymmetry of the Balmer lines. These lines can be well decomposed into a narrow component and an intermediate-width component. The intermediate-width component showed a steady increase in both strength and blueshift with time until t ~ 400 days after maximum, but it became less blueshifted at t ~ 500 days when the line profile appeared relatively symmetric again. Owing to that a pure reddening effect will lead to a sudden decline of the light curves and a progressive blueshift of the spectral lines, we therefore propose that the asymmetric profiles of H lines seen in SN 2010jl is unlikely due to the extinction by newly formed dust inside the ejecta, contrary to the explanation by some early studies. Based on a simple CSM-interaction model, we speculate that the progenitor of SN 2010jl may suffer a gigantic mass loss (~ 30-50 M_sun) in a few decades before explosion. Considering a slow moving stellar wind (e.g., ~ 28 km/s) inferred for the preexisting, dense CSM shell and the extremely high mass-loss rate (1-2 M_sun per yr), we suggest that the progenitor of SN 2010jl might have experienced a red supergiant stage and explode finally as a post-red supergiant star with an initial mass above 30-40 M_sun.

Well-sampled optical lightcurves of 146 gamma-ray bursts (GRBs) are compiled from the literature. By empirical fitting we identify eight possible emission components and summarize the results in a "synthetic" lightcurve. Both optical flare and early shallow-decay components are likely related to long-term central engine activities. We focus on their statistical properties in this paper. Twenty-four optical flares are obtained from 19 GRBs. The isotropic R-band energy is smaller than 1% of $E_{\gamma, \rm iso}$. The relation between isotropic luminosities of the flares and gamma-rays follows $L^{\rm F}_{\rm R, iso}\propto L_{{\gamma}, \rm iso}^{1.11\pm 0.27}$. Later flares tend to be wider and dimmer, i.e., $w^{\rm F}\sim t^{\rm F}_{\rm p}/2$ and $L^{\rm F}_{\rm R, iso}\propto [t^{\rm F}_{\rm p}/(1+z)]^{-1.15\pm0.15}$. The detection probability of the optical flares is much smaller than that of X-ray flares. An optical shallow decay segment is observed in 39 GRBs. The relation between the break time and break luminosity is a power-law, with an index of $-0.78\pm 0.08$, similar to that derived from X-ray flares. The X-ray and optical breaks are usually chromatic, but a tentative correlation is found. We suggest that similar to the prompt optical emission that tracks $\gamma$-rays, the optical flares are also related to the erratic behavior of the central engine. The shallow decay component is likely related to a long-lasting spinning-down central engine or piling up of flare materials onto the blastwave. Mixing of different emission components may be the reason of the diverse chromatic afterglow behaviors.

We describe the target selection algorithm for the low latitude disk portion of the LAMOST Pilot Survey, which aims to test systems in preparation for the LAMOST spectroscopic survey. We use the PPMXL (Roeser et al. 2010) astrometric catalog, which provides positions, proper motions, B/R/I magnitudes (mostly) from USNO-B (Monet et al. 2003) and J/H/Ks from The Two Micron All Sky Survey (2MASS, see Skrutskie et al. 2006) as well. We chose 8 plates along the Galactic plane, in the region $0^\circ<\alpha<67^\circ$ and $42^\circ<\delta<59^\circ$, that cover 22 known open clusters with a range of ages. Adjacent plates may have small overlapping. Each plate covers an area $2.5^\circ$ in radius,with central star (for Shack-Hartmann guider) brighter than $\sim8^{\rm th}$ magnitude. For each plate, we create an input catalog in the magnitude range $11.3<Imag<16.3$ and $Bmag$ available from PPMXL. The stars are selected to satisfy the requirements of the fiber positioning system and have a uniform distribution in the $I$ vs. $B-I$ color-magnitude diagram. Our final input catalog consists of 12,000 objects on each of 8 plates that are observable during the winter observing season in Xinglong Station of the National Astronomical Observatory of China.

We extend the holographic Ricci dark energy model to include some direct, non-gravitational interaction between dark energy and dark matter. We consider three phenomenological forms for the interaction term $Q$ in the model, namely, $Q$ is taken proportional to the Hubble expansion rate and the energy densities of dark sectors (taken to be $\rho_{\rm de}$, $\rho_{\rm m}$, and $\rho_{\rm de}+\rho_{\rm m}$, respectively). We obtain a uniform analytical solution to the three interacting models. Furthermore, we constrain the models by using the latest observational data, including the 557 Union2 type Ia supernovae data, the cosmic microwave background anisotropy data from the 7-yr WMAP, and the baryon acoustic oscillation data from the SDSS. We show that in the interacting models of the holographic Ricci dark energy, a more reasonable value of $\Omega_{\rm m0}$ will be obtained, and the observations favor a rather strong coupling between dark energy and dark matter.

Well-sampled optical lightcurves of 146 gamma-ray bursts (GRBs) are complied from the literature. Multiple optical emission components are extracted with power-law function fits to these lightcurves. We present a systematical analysis for statistical properties and their relations to prompt gamma-ray emission and X-ray afterglow for each component. We show that peak luminosity in the prompt and late flares are correlated and the evolution of the peak luminosity may signal the evolution of the accretion rate. No tight correlation between the shallow decay phase/plateau and prompt gamma-ray emission is found. Assuming that they are due to a long-lasting wind injected by a compact object, we show that the injected behavior favors the scenarios of a long-lasting wind after the main burst episode. The peak luminosity of the afterglow onset is tightly correlated with Eiso, and it is dimmer as peaking later. Assuming that the onset bump is due to the fireball deceleration by the external medium, we examine the Gamma_0-Eiso relation and find that it is confirmed with the current sample. Optical re-brightening is observed in 30 GRBs in our sample. It shares the same relation between the width and the peak time as found in the onset bump, but no clear correlation between the peak luminosity and Eiso as observed in the onset bumps is found. Although its peak luminosity also decays with time, the slope is much shallower than that of the onset peak. We get L t^{-1}_{p}$, being consistent with off-axis observations to an expanding external fireball in a wind-like circum medium. The late re-brightening may signal another jet component. Mixing of different emission components may be the reason for the observed chromatic breaks in different energy bands.

Models of pure gas-phase chemistry in well-shielded regions of molecular clouds predict relatively high levels of molecular oxygen, O2, and water, H2O. Contrary to expectation, the space missions SWAS and Odin found only very small amounts of water vapour and essentially no O2 in the dense star-forming interstellar medium. Only toward rho Oph A did Odin detect a weak line of O2 at 119 GHz in a beam size of 10 arcmin. A larger telescope aperture such as that of the Herschel Space Observatory is required to resolve the O2 emission and to pinpoint its origin. We use the Heterodyne Instrument for the Far Infrared aboard Herschel to obtain high resolution O2 spectra toward selected positions in rho Oph A. These data are analysed using standard techniques for O2 excitation and compared to recent PDR-like chemical cloud models. The 487.2GHz line was clearly detected toward all three observed positions in rho Oph A. In addition, an oversampled map of the 773.8GHz transition revealed the detection of the line in only half of the observed area. Based on their ratios, the temperature of the O2 emitting gas appears to vary quite substantially, with warm gas (> 50 K) adjacent to a much colder region, where temperatures are below 30 K. The exploited models predict O2 column densities to be sensitive to the prevailing dust temperatures, but rather insensitive to the temperatures of the gas. In agreement with these model, the observationally determined O2 column densities seem not to depend strongly on the derived gas temperatures, but fall into the range N(O2) = (3 to >6)e15/cm^2. Beam averaged O2 abundances are about 5e-8 relative to H2. Combining the HIFI data with earlier Odin observations yields a source size at 119 GHz of about 4 - 5 arcmin, encompassing the entire rho Oph A core.

We explore the halo-to-halo variation of dark matter substructure in galaxy-sized dark matter halos, focusing on its implications for strongly gravitational lensed systems. We find that the median value for projected substructure mass fractions within projected radii of 3% of the host halo virial radius is approximately fsub ~ 0.25%, but that the variance is large with a 95-percentile range of 0 <= fsub <= 1%. We quantify possible effects of substructure on quadruply-imaged lens systems using the cusp relation and the simple statistic, Rcusp. We estimate that the probability of obtaining the large values of the Rcusp which have been observed from substructure effects is roughly ~10^-3 to ~10^-2. We consider a variety of possible correlations between host halo properties and substructure properties in order to probe possible sample biases. In particular, low-concentration host dark matter halos have more large substructures and give rise to large values of Rcusp more often. However, there is no known observational bias that would drive observed quadruply-imaged quasars to be produced by low-concentration lens halos. Finally, we show that the substructure mass fraction is a relatively reliable predictor of the value of Rcusp.

We report observations of three rotational transitions of molecular oxygen (O2) in emission from the H2 Peak 1 position of vibrationally excited molecular hydrogen in Orion. We observed the 487 GHz, 774 GHz, and 1121 GHz lines using HIFI on the Herschel Space Observatory, having velocities of 11 km s-1 to 12 km s-1 and widths of 3 km s-1. The beam-averaged column density is N(O2) = 6.5\times1016 cm-2, and assuming that the source has an equal beam filling factor for all transitions (beam widths 44, 28, and 19"), the relative line intensities imply a kinetic temperature between 65 K and 120 K. The fractional abundance of O2 relative to H2 is 0.3 - 7.3\times10-6. The unusual velocity suggests an association with a ~ 5" diameter source, denoted Peak A, the Western Clump, or MF4. The mass of this source is ~ 10 M\odot and the dust temperature is \geq 150 K. Our preferred explanation of the enhanced O2 abundance is that dust grains in this region are sufficiently warm (T \geq 100 K) to desorb water ice and thus keep a significant fraction of elemental oxygen in the gas phase, with a significant fraction as O2. For this small source, the line ratios require a temperature \geq 180 K. The inferred O2 column density \simeq 5\times1018 cm-2 can be produced in Peak A, having N(H2) \simeq 4\times1024 cm-2. An alternative mechanism is a low-velocity (10 to 15 km s-1) C-shock, which can produce N(O2) up to 1017 cm-2.

The mode switching phenomenon of PSR B0329+54 is investigated based on the long-term monitoring from September 2003 to April 2009 made with the Urumqi 25m radio telescope at 1540 MHz. At that frequency, the change of relative intensity between the leading and trailing components is the predominant feature of mode switching. The intensity ratios between the leading and trailing components are measured for the individual profiles averaged over a few minutes. It is found that the ratios follow normal distributions, where the abnormal mode has a wider typical width than the normal mode, indicating that the abnormal mode is less stable than the normal mode. Our data show that 84.9% of the time for PSR B0329+54 was in the normal mode and 15.1% was in the abnormal mode. From the two passages of eight-day quasi-continuous observations in 2004, and supplemented by the daily data observed with 15 m telescope at 610 MHz at Jodrell Bank Observatory, the intrinsic distributions of mode timescales are constrained with the Bayesian inference method. It is found that the gamma distribution with the shape parameter slightly smaller than 1 is favored over the normal, lognormal and Pareto distributions. The optimal scale parameters of the gamma distribution is 31.5 minutes for the abnormal mode and 154 minutes for the normal mode. The shape parameters have very similar values, i.e. 0.75^{+0.22}_{-0.17} for the normal mode and 0.84^{+0.28}_{-0.22} for the abnormal mode, indicating the physical mechanisms in both modes may be the same. No long-term modulation of the relative intensity ratios was found for both the modes, suggesting that the mode switching was stable. The intrinsic timescale distributions, for the first time constrained for this pulsar, provide valuable information to understand the physics of mode switching.

A number of spectroscopic surveys have been carried out or are planned to study the origin of the Milky Way. Their exploitation requires reliable automated methods and softwares to measure the fundamental parameters of the stars. Adopting the ULySS package, we have tested the effect of different resolutions and signal-to-noise ratios (SNR) on the measurement of the stellar atmospheric parameters (effective temperature Teff, surface gravity log g, and metallicity [Fe/H]). We show that ULySS is reliable to determine these parameters with medium-resolution spectra (R~2000). Then, we applied the method to measure the parameters of 771 stars selected in the commissioning database of the Guoshoujing Telescope (GSJT). The results were compared with the SDSS/SEGUE Stellar Parameter Pipeline (SSPP), and we derived precisions of 167 K, 0.34 dex, and 0.16 dex for Teff, log g and [Fe/H] respectively. Furthermore, 120 of these stars are selected to construct the primary stellar spectra template library (Version 1.0) of GSJT, and will be deployed as basic ingredients for the GSJT automated parametrization pipeline.

The first release of the Chandra Source Catalog (CSC) contains ~95,000 X-ray sources in a total area of ~0.75% of the entire sky, using data from ~3,900 separate ACIS observations of a multitude of different types of X-ray sources. In order to maximize the scientific benefit of such a large, heterogeneous data-set, careful characterization of the statistical properties of the catalog, i.e., completeness, sensitivity, false source rate, and accuracy of source properties, is required. Characterization efforts of other, large Chandra catalogs, such as the ChaMP Point Source Catalog (Kim et al. 2007) or the 2 Mega-second Deep Field Surveys (Alexander et al. 2003), while informative, cannot serve this purpose, since the CSC analysis procedures are significantly different and the range of allowable data is much less restrictive. We describe here the characterization process for the CSC. This process includes both a comparison of real CSC results with those of other, deeper Chandra catalogs of the same targets and extensive simulations of blank-sky and point source populations.

We analyze the submillimeter emission surrounding the new FU Orionis-type object, HBC 722. We present the first epoch of observations of the active environs of HBC 722, with imaging and spectroscopy from PACS, SPIRE, and HIFI aboard the Herschel Space Observatory, as well as CO J= 2-1 and 350 um imaging (SHARC-II) with the Caltech Submillimeter Observatory. The primary source of submillimeter continuum emission in the region -- 2MASS 20581767+4353310 -- is located 16$\arcsec$ south-southeast of the optical flaring source while the optical and near-IR emission is dominated by HBC 722. A bipolar outflow extends over HBC 722; the most likely driver is the submillimeter source. We detect warm (100 K) and hot (246 K) CO emission in the surrounding region, evidence of outflow-driven heating in the vicinity. The region around HBC 722 itself shows little evidence of heating driven by the new outbursting source itself.

In this paper we present a multicolor photometry for A2589 ($z=0.0414$) with 15 intermediate bands in the Beijing-Arizona-Taiwan-Connecticut (BATC) system which covers an optical wavelength range from 3000 \AA\ to 10000 \AA. The spectral energy distributions (SEDs) for more than 5000 sources are achieved down to {\it V} $\sim$ 20 mag in about 1 deg$^{2}$ field. A2589 has been also covered by the Sloan Digital Sky Survey (SDSS) in photometric mode only. A cross-identification of the BATC-detected galaxies with the SDSS photometric catalog achieves 1199 galaxies brighter than $i=19.5$ mag, among which 68 member galaxies with known spectroscopic redshifts are found. After combining the SDSS five-band photometric data and the BATC SEDs, the technique of photometric redshift is applied to these galaxies for selecting faint member galaxies. The color-magnitude relation is taken as a further restriction of early-type cluster galaxies. As a result, 106 galaxies are newly selected as member galaxies. Spatial distribution of member galaxies shows a north-south elongation which agrees with the X-ray brightness profile and the orientation of central cD galaxy, NGC 7647. No substructures are detected on the basis of positions and radial velocities of cluster galaxies, indicating that A2589 is a well-relaxed system. The luminosity function of A2589 exhibits a peak at $M_{R} \sim -20$ mag and a dip at $M_{R} \sim -19 $ mag. The low-density outer regions are the preferred habitat of faint galaxies. With the evolutionary population synthesis model, PEGASE, the environmental effect on the star formation properties for 68 spectroscopically confirmed member galaxies is studied. The outlier faint galaxies tend to have longer time scales of star formation, shorter mean stellar ages, and lower metallicities of interstellar medium, which can be interpreted in the context of hierarchical cosmological scenario.

We present new observations of the CB130 region, composed of three separate cores. Using the \textit{Spitzer Space Telescope} we detected a Class 0 and a Class II object in one of these, CB130-1. The observed photometric data from \textit{Spitzer} and ground-based telescopes are used to establish the physical parameters of the Class 0 object. SED fitting with a radiative transfer model shows that the luminosity of the Class 0 object is 0.14 $-$ 0.16 L$_{\odot}$, which is a low luminosity for a protostellar object. In order to constrain the chemical characteristics of the core having the low luminosity object, we compare our molecular line observations to models of lines including abundance variations. We tested both ad hoc step function abundance models and a series of self-consistent chemical evolution models. In the chemical evolution models, we consider a continuous accretion model and an episodic accretion model to explore how variable luminosity affects the chemistry. The step function abundance models can match observed lines reasonably well. The best fitting chemical evolution model requires episodic accretion and the formation of CO$_2$ ice from CO ice during the low luminosity periods. This process removes C from the gas phase, providing a much improved fit to the observed gas-phase molecular lines and the CO$_2$ ice absorption feature. Based on the chemical model result, the low luminosity of CB130-1 is explained better as a quiescent stage between episodic accretion bursts rather than being at the first hydrostatic core stage.

We summarize and critically evaluate the available data on nuclear fusion cross sections important to energy generation in the Sun and other hydrogen-burning stars and to solar neutrino production. Recommended values and uncertainties are provided for key cross sections, and a recommended spectrum is given for 8B solar neutrinos. We also discuss opportunities for further increasing the precision of key rates, including new facilities, new experimental techniques, and improvements in theory. This review, which summarizes the conclusions of a workshop held at the Institute for Nuclear Theory, Seattle, in January 2009, is intended as a 10-year update and supplement to Reviews of Modern Physics 70 (1998) 1265.

We investigate the feasibility of extracting the gravitational nanolensing signal due to the presence of subsolar mass halos within galaxy-sized dark matter halos. We show that subsolar mass halos in a lensing galaxy can cause strong nanolensing events with shorter durations and smaller amplitudes than microlensing events caused by stars. We develop techniques that can be used in future surveys such as Pan-STARRS, LSST and OMEGA to search for the nanolensing signal from subsolar mass halos.

Properties of proton resonances in $^{18}$Ne have been investigated efficiently by utilizing a technique of proton resonant elastic scattering with a $^{17}$F radioactive ion (RI) beam and a thick proton target. A 4.22~MeV/nucleon $^{17}$F RI beam was produced via a projectile-fragmentation reaction, and subsequently separated by a Radioactive Ion Beam Line in Lanzhou ({\tt RIBLL}). Energy spectra of the recoiled protons were measured by two sets of $\Delta$E-E silicon telescope at center-of-mass scattering angles of $\theta_{c.m.}$$\approx$175${^\circ}$$\pm$5${^\circ}$, $\theta_{c.m.}$$\approx$152${^\circ}$$\pm$8${^\circ}$, respectively. Several proton resonances in $^{18}$Ne were observed, and their resonant parameters have been determined by an $R$-matrix analysis of the differential cross sections in combination with the previous results. The resonant parameters are related to the reaction-rate calculation of the stellar $^{14}$O($\alpha$,$p$)$^{17}$F reaction, which was thought to be the breakout reaction from the hot CNO cycles into the $rp$-process in x-ray bursters. Here, $J^\pi$=(3$^-$, 2$^-$) are tentatively assigned to the 6.15-MeV state which was thought the key 1$^-$ state previously. In addition, a doublet structure at 7.05 MeV are tentatively identified, and its contribution to the resonant reaction rate of $^{14}$O($\alpha$,$p$)$^{17}$F could be enhanced by at least factors of about 4$\sim$6 in comparison with the previous estimation involving only a singlet. The present calculated resonant rates are much larger than those previous values, and it may imply that this breakout reaction could play a crucial role under x-ray bursters conditions.

M31 globular cluster B379 is the first extragalactic cluster, the age of which was determined by main-sequence photometry. In this method, the age of a cluster is obtained by fitting its CMD with stellar evolutionary models. However, different stellar evolutionary models use different parameters of stellar evolution, such as range of stellar masses, different opacities and equations of state, and different recipes, and so on. So, it is interesting to check whether different stellar evolutionary models can give consistent results for the same cluster. Brown et al. (2004a) constrained the age of B379 by comparing its CMD with isochrones of the 2006 VandenBerg models. Using SSP models of BC03 and its multi-photometry, Ma et al. (2007) independently determined the age of B379, which is in good agreement with the determination of Brown et al. (2004a). The BC03 models are calculated based on the Padova evolutionary tracks. It is necessary to check whether the age of B379 which, being determined based on the Padova evolutionary tracks, is in agreement with the determination of Brown et al. (2004a). So, in this paper, we re-determine its age using isochrones of the Padova stellar evolutionary models. In addition, the metal abundance, the distance modulus, and the reddening value for B379 are also determined in this paper. The results obtained in this paper are consistent with the previous determinations, which including the age obtained by Brown et al. (2004a). So, this paper confirms the consistence of the age scale of B379 between the Padova isochrones and the 2006 VandenBerg isochrones, i.e. the results' comparison between Brown et al. (2004a) and Ma et al. (2007) is meaningful. The results obtained in this paper are: the metallicity [M/H]=-0.325, the age $\tau=11.0\pm1.5$ Gyr, the reddening value E(B-V)=0.08, and the distance modulus $(m-M)_{0}=24.44\pm0.10$.

Neutron-rich isotopes with masses near that of iron are produced in type Ia and II supernovae. Traces of such nucleosynthesis are found in primitive meteorites in the form of variations in the isotopic abundance of 54Cr, the most neutron-rich stable isotope of chromium. The hosts of these isotopic anomalies must be presolar grains that condensed in the outflows of supernovae, offering the opportunity to study the nucleosynthesis of iron-peak nuclei in ways that complement spectroscopic observations and can inform models of stellar evolution. However, despite almost two decades of extensive search, the carrier of 54Cr anomalies is still unknown, presumably because it is fine-grained and is chemically labile. Here we identify in the primitive meteorite Orgueil the carrier of 54Cr-anomalies as nanoparticles, most likely spinels that show large enrichments in 54Cr relative to solar composition (54Cr/52Cr ratio >3.6xsolar). Such large enrichments in 54Cr can only be produced in supernovae. The mineralogy of the grains supports condensation in the O/Ne-O/C zones of a type II supernova, although a type Ia origin cannot be excluded. We suggest that planetary materials incorporated different amounts of these nanoparticles, possibly due to late injection by a nearby supernova that also delivered 26Al and 60Fe to the solar system. This idea explains why the relative abundance of 54Cr and other neutron-rich isotopes vary between planets and meteorites. We anticipate that future isotopic studies of the grains identified here will shed new light on the birth of the solar system and the conditions insupernovae.

We report the discovery of eight new quasars in one extragalactic field (five degree centered at RA=$08^h58^m08.2^s$, Dec=$01^o32'29.7''$) with the LAMOST commissioning observations on December 18, 2009. These quasars, with $i$ magnitudes from 16.44 to 19.34 and redshifts from 0.898 to 2.773, were not identified in the SDSS spectroscopic survey, though six of them with redshifts less than 2.5 were selected as quasar targets in SDSS. Except one source without near-IR $Y$-band data, seven of these eight new quasars meet a newly proposed quasar selection criterion involving both near-IR and optical colors. Two of them were found in the 'redshift desert' for quasars ($z$ from 2.2 to 3), indicating that the new criterion is efficient for recovering the missing quasars with similar optical colors as stars. Although LAMOST met some problems during the commissioning observations, we were still able to identify other 38 known SDSS quasars in this field, with $i$ magnitudes from 16.24 to 19.10 and redshifts from 0.297 to 4.512. Our identifications imply that a substantial fraction of quasars may be missing in the previous quasar surveys. The implication of our results to the future LAMOST quasar survey is discussed.

The redshift range from 2.2 to 3, is known as the 'redshift desert' of quasars because quasars with redshift in this range have similar optical colors as normal stars and are thus difficult to be found in optical sky surveys. A quasar candidate, SDSS J085543.40-001517.7, which was selected by a recently proposed criterion involving near-IR $Y-K$ and optical $g-z$ colors, was identified spectroscopically as a new quasar with redshift of 2.427 by the LAMOST commissioning observation in December 2009 and confirmed by the observation made with the NAOC/Xinglong 2.16m telescope in March 2010. This quasar was not targeted in the SDSS spectroscopic survey because it locates in the stellar locus of the optical color-color diagrams, while it is clearly separated from stars in the $Y-K$ vs. $g-z$ diagram. Comparing with other SDSS quasars we found this new quasar with $i$ magnitude of 16.44 is apparently the brightest one in the redshift range from 2.3 to 2.7. From the spectral properties we derived its central black hole mass as $(1.4\sim3.9) \times 10^{10} M_\odot$ and the bolometric luminosity as $3.7\times 10^{48}$ \ergs, which indicates that this new quasar is intrinsically very bright and belongs to the most luminous quasars in the universe. Our identification supports that quasars in the redshift desert can be found by the quasar selection criterion involving the near-IR colors. More missing quasars are expected to be recovered by the future LAMOST spectroscopic surveys, which is important to the study of the cosmological evolution of quasars at redshift higher than 2.2.

The Chandra Source Catalog (CSC) is a general purpose virtual X-ray astrophysics facility that provides access to a carefully selected set of generally useful quantities for individual X-ray sources, and is designed to satisfy the needs of a broad-based group of scientists, including those who may be less familiar with astronomical data analysis in the X-ray regime. The first release of the CSC includes information about 94,676 distinct X-ray sources detected in a subset of public ACIS imaging observations from roughly the first eight years of the Chandra mission. This release of the catalog includes point and compact sources with observed spatial extents <~ 30''. The catalog (1) provides access to the best estimates of the X-ray source properties for detected sources, with good scientific fidelity, and directly supports scientific analysis using the individual source data; (2) facilitates analysis of a wide range of statistical properties for classes of X-ray sources; and (3) provides efficient access to calibrated observational data and ancillary data products for individual X-ray sources, so that users can perform detailed further analysis using existing tools. The catalog includes real X-ray sources detected with flux estimates that are at least 3 times their estimated 1 sigma uncertainties in at least one energy band, while maintaining the number of spurious sources at a level of <~ 1 false source per field for a 100 ks observation. For each detected source, the CSC provides commonly tabulated quantities, including source position, extent, multi-band fluxes, hardness ratios, and variability statistics, derived from the observations in which the source is detected. In addition to these traditional catalog elements, for each X-ray source the CSC includes an extensive set of file-based data products that can be manipulated interactively.

We present far-infrared spectroscopic observations, taken with the Photodetector Array Camera and Spectrometer (PACS) on the Herschel Space Observatory, of the protoplanetary disk around the pre-main-sequence star HD 100546. These observations are the first within the DIGIT Herschel key program, which aims to follow the evolution of dust, ice, and gas from young stellar objects still embedded in their parental molecular cloud core, through the final pre-main-sequence phases when the circumstellar disks are dissipated. Our aim is to improve the constraints on temperature and chemical composition of the crystalline olivines in the disk of HD 100546 and to give an inventory of the gas lines present in its far-infrared spectrum. The 69 \mu\m feature is analyzed in terms of position and shape to derive the dust temperature and composition. Furthermore, we detected 32 emission lines from five gaseous species and measured their line fluxes. The 69 \mu\m emission comes either from dust grains with ~70 K at radii larger than 50 AU, as suggested by blackbody fitting, or it arises from ~200 K dust at ~13 AU, close to the midplane, as supported by radiative transfer models. We also conclude that the forsterite crystals have few defects and contain at most a few percent iron by mass. Forbidden line emission from [CII] at 157 \mu\m and [OI] at 63 and 145 \mu\m, most likely due to photodissociation by stellar photons, is detected. Furthermore, five H2O and several OH lines are detected. We also found high-J rotational transition lines of CO, with rotational temperatures of ~300 K for the transitions up to J=22-21 and T~800 K for higher transitions.

DK Cha is an intermediate-mass star in transition from an embedded configuration to a star plus disk stage. We aim to study the composition and energetics of the circumstellar material during this pivotal stage. Using the Range Scan mode of PACS on the Herschel Space Observatory, we obtained a spectrum of DK Cha from 55 to 210 micron as part of the DIGIT Key Program. Almost 50 molecular and atomic lines were detected, many more than the 7 lines detected in ISO-LWS. Nearly the entire ladder of CO from J=14-13 to 38-37 (E_u/k = 4080 K), water from levels as excited as E_u/k = 843 K, and OH lines up to E_u/k = 290 K were detected. The continuum emission in our PACS SED scan matches the flux expected from a model consisting of a star, a surrounding disk of 0.03 Solar mass, and an envelope of a similar mass, supporting the suggestion that the object is emerging from its main accretion stage. Molecular, atomic, and ionic emission lines in the far-infrared reveal the outflow's influence on the envelope. The inferred hot gas can be photon-heated, but some emission could be due to C-shocks in the walls of the outflow cavity.

In this paper, we present an automated system, which has the capability to catch and track solar limb prominences based on observations from EUV 304 passband. The characteristic parameters and their evolution, including height, position angle, area, length and brightness, are obtained without manual interventions. By applying the system to the STEREO-B/SECCHI/EUVI 304 data during 2007 April -2009 October, we obtain a total of 9477 well-tracked prominences and a catalog of these events available online at http://space.ustc.edu.cn/dreams/slipcat/. A detailed analysis of these prominences suggests that the system has a rather good performance. We have obtained several interesting statistical results based on the catalog. Most prominences appear below the latitude of 60 degrees and at the height of about 26 Mm above the solar surface. Most of them are quite stable during the period they are tracked. Nevertheless, some prominences have an upward speed of more than 100 km/s, and some others show significant downward and/or azimuthal speeds. There are strong correlations among the brightness, area and height. The expansion of a prominence is probably one major cause of its fading during the rising or erupting process.

We present results from the Spitzer Space Telescope and molecular line observations of 9 species toward the dark cloud L43. The Spitzer images and molecular line maps suggest it has a starless core and a Class I protostar evolving in the same environment. CO depletion is seen in both sources, and DCO+ lines are stronger toward the starless core. With a goal of testing the chemical characteristics from pre- to protostellar stages, we adopt an evolutionary chemical model to calculate the molecular abundances and compare with our observations. Among the different model parameters we tested, the best-fit model suggests a longer total timescale at the pre-protostellar stage, but with faster evolution at the later steps with higher densities.

Circumstellar discs are expected to be the nursery of planets. Grain growth within such discs is the first step in the planet formation process. The Bok globule CB 26 harbours such a young disc. We present a detailed model of the edge-on circumstellar disc and its envelope in the Bok globule CB 26. The model is based on HST near-infrared maps in the I, J, H, and K bands, OVRO and SMA radio maps at 1.1mm, 1.3mm and 2.7mm, and the spectral energy distribution (SED) from 0.9 microns to 3mm. New photometric and spectroscopic data from the Spitzer Space Telescope and the Caltech Submilimeter Observatory have been obtained and are part of our analysis. Using the self-consistent radiative transfer code MC3D, the model we construct is able to discriminate parameter sets and dust properties of both its parts, namely envelope and disc. We find that the disc has an inner hole with a radius of 45 +/- 5 AU. Based on a dust model including silicate and graphite the maximum grain size needed to reproduce the spectral millimetre index is 2.5 microns. Features seen in the near-infrared images, dominated by scattered light, can be described as a result of a rotating envelope. Successful employment of ISM dust in both the disc and envelope hint that grain growth may not yet play a significant role for the appearance of this system. A larger inner hole gives rise to the assumption that CB 26 is a circumbinary disc.

We report the results of a time-series CCD photometric survey of variable stars in the field of the open cluster NGC 2126. In about one square degree field covering the cluster, a total number of 21 variable candidates are detected during this survey, of which 16 are newly found. The periods, classifications and spectral types of 14 newly discovered variables are discussed, which consist of six eclipsing binaries systems, three pulsating variable stars, three long period variables, one RS CVn star, and one W UMa or $\delta$ Scuti star. In addition, there are two variable candidates, the properties of which cannot be determined in this paper. By a method based on fitting spectral energy distributions(SEDs) of stars with theoretical ones, the membership probabilities and the fundamental parameters of this cluster are determined. As a result, five variables are probably members of NGC 2126. The fundamental parameters of the this cluster are determined as: the metallicity to be 0.008 $Z_\odot$, the age $\log(t)$= 8.95, the distance modulus $(m-M)_{0}$= 10.34 and the reddening value $E(B-V)$=0.55 mag.

There is a vast menagerie of plausible candidates for the constituents of dark matter, both within and beyond extensions of the Standard Model of particle physics. Each of these candidates may have scattering (and other) cross section properties that are consistent with the dark matter abundance, BBN, and the most scales in the matter power spectrum; but which may have vastly different behavior at sub-galactic "cutoff" scales, below which dark matter density fluctuations are smoothed out. The only way to quantitatively measure the power spectrum behavior at sub-galactic scales at distances beyond the local universe, and indeed over cosmic time, is through probes available in multiply imaged strong gravitational lenses. Gravitational potential perturbations by dark matter substructure encode information in the observed relative magnifications, positions, and time delays in a strong lens. Each of these is sensitive to a different moment of the substructure mass function and to different effective mass ranges of the substructure. The time delay perturbations, in particular, are proving to be largely immune to the degeneracies and systematic uncertainties that have impacted exploitation of strong lenses for such studies. There is great potential for a coordinated theoretical and observational effort to enable a sophisticated exploitation of strong gravitational lenses as direct probes of dark matter properties. This opportunity motivates this white paper, and drives the need for: a) strong support of the theoretical work necessary to understand all astrophysical consequences for different dark matter candidates; and b) tailored observational campaigns, and even a fully dedicated mission, to obtain the requisite data.

The spatial distribution of satellite galaxies around host galaxies can illuminate the relationship between satellites and dark matter subhalos and aid in developing and testing galaxy formation models. Previous efforts to constrain the distribution attempted to eliminate interlopers from the measured projected number density of satellites and found that the distribution is generally consistent with the expected dark matter halo profile of the parent hosts, with a best-fit power-law slope of ~ -1.7 between projected separations of ~30 kpc/h and 0.5 Mpc/h. Here, I use the projected cross-correlation of bright and faint galaxies to analyze contributions from satellites and interlopers together, using a halo occupation distribution (HOD) analytic model for galaxy clustering. This approach is tested on mock catalogs constructed from simulations. I find that analysis of Sloan Digital Sky Survey (SDSS) data gives results generally consistent with interloper subtraction methods between the projected separations of 10 kpc/h and 6.3 Mpc/h, although the errors on the parameters that constrain the radial profile are large, and larger samples of data are required.

The source J 1128+5925 was found recently to show strong intraday variability at radio wavelengths and its radio variability may come from interstellar scintillation. In optical, the object was quiet in our 2007 monitoring session. Here we report the results of our new optical monitoring of this source in 2008. In addition to confirm our 2007 results, that the object did not display any clear variation on timescales from hour--day to month, we provide evidence that the object does not vary on timescale of one year, and it is probably intrinsically quiet in optical domain. Its very different behaviors in optical and radio regimes can be naturally explained if its strong radio variability comes from interstellar scintillation.

We describe a Windows compatible version of the evolutionary synthesis code Starburst99. Starburst99 for Windows was developed from the public UNIX based version at STScI. We converted the original Fortran77 source code into a version for a Win32 environment with an Absoft Fortran Pro x86 compiler. Extensive testing showed no significant numerical differences in comparison with the previous UNIX version. The software application consists of the source code, executable, and a number of auxiliary files. The package installs on any PC running Windows 2000, XP, or Vista and can be obtained as freeware at http://www.stsci.edu/science/starburst/PCStarburst99.html. We give an overview of the different running modes and provide instructions for getting started with the initial set-up.

The amount of mass in small, dark matter clumps within galaxies (substructure) is an important test of cold dark matter. One approach to measuring the substructure mass fraction is to analyze the fluxes of images that have been strongly lensed by a galaxy. Flux ratios between images that are anomalous with respect to smooth (no substructure) models have previously suggested that there is a greater amount of substructure than found in dark matter simulations. One measure of anomalous flux ratios is parity dependence -- that the fluxes of different images of a source are perturbed differently. In this paper, we discuss parity dependence as a probe of dark matter substructure. We find that reproducing the observed parity dependence requires a significant alignment between concentrated dark matter clumps and images. The results may imply a larger fraction of mass in substructures than suggested by some dark matter simulations and that the observed parity dependence is unlikely to be reproduced by luminous satellites of lens galaxies.

The dimensionless universal coefficient $\xi$ defines the ratio of the unitary fermions energy density to that for the ideal non-interacting ones in the non-relativistic limit with T=0. The classical Thomson Problem is taken as a nonperturbative quantum many-body arm to address the ground state energy including the low energy nonlinear quantum fluctuation/correlation effects. With the relativistic Dirac continuum field theory formalism, the concise expression for the energy density functional of the strongly interacting limit fermions at both finite temperature and density is obtained. Analytically, the universal factor is calculated to be $\xi={4/9}$. The energy gap is $\Delta=\frac{{5}{18}{k_f^2}/(2m)$.

The carrier of the 2175 Angstrom interstellar extinction feature remains unidentified since its first detection over 40 years ago. In recent years carbon buckyonions have been proposed as a carrier of this feature, based on the close similarity between the electronic transition spectra of buckyonions and the 2175 Angstrom interstellar feature. We examine this hypothesis by modeling the interstellar extinction with buckyonions as a dust component. It is found that dust models containing buckyonions (in addition to amorphous silicates, polycyclic aromatic hydrocarbon molecules, graphite) can closely reproduce the observed interstellar extinction curve. To further test this hypothesis, we call for experimental measurements and/or theoretical calculations of the infrared vibrational spectra of hydrogenated buckyonions. By comparing the infrared emission spectra predicted for buckyonions vibrationally excited by the interstellar radiation with the observed emission spectra of the diffuse interstellar medium, we will be able to derive (or place an upper limit on) the abundance of interstellar buckyonions.

We explore the color dependence of the radial profile of satellite galaxies around isolated parent galaxies. Samples of potential satellites selected from large galaxy redshift surveys are significantly contaminated by interlopers -- objects not bound to the parent galaxy. We use the Sloan Digital Sky Survey to estimate the interloper fraction in samples of candidate satellite galaxies. We show that samples of red and blue satellites have different interloper populations: a larger fraction of blue galaxies are likely to be interlopers compared to red galaxies. Both with and without interloper subtraction, the radial profile of blue satellites is significantly shallower than that of red satellites. In addition, while red and blue primaries have different interloper fractions, the slope of the corrected radial profiles are consistent after interloper correction. We discuss the implications of these results for galaxy formation models.

The lensing cross section of triaxial halos depends on the relative orientation between a halo's principal axes and its line of sight. Consequently, a lensing subsample of randomly oriented halos is not, in general,randomly oriented. Using an isothermal mass model for the lensing galaxies and their host halos, we show that the lensing subsample of halos that produces doubles is preferentially aligned along the lines of sight, whereas halos that produce quads tend to be projected along their middle axes. These preferred orientations result in different projected ellipticity distributions for quad, doubles, and random galaxies. We show that ~300 lens systems must be discovered to detect this effect at the 95% confidence level. We also investigate the importance of halo shape for predicting the quad-to-double ratio and find that the latter depends quite sensitively on the distribution of the short-to-long axis ratio, but is otherwise nearly independent of halo shape. Finally, we estimate the impact of the preferred orientation of lensing galaxies on their projected substructure mass fraction, and find that the observed alignment between the substructure distribution and the mass distribution of halos result in a negligible bias.

Very recently, J 1128+5925 was found to show strong intraday variability at radio wavelengths and may be a new source with annual modulation of the timescale of its radio variability. Therefore, its radio variability can be best explained via interstellar scintillation. Here we present the properties of its optical variability for the first time after a monitoring program in 2007 May. Our observations indicate that in this period J 1128+5925 only showed trivial optical variability on internight timescale, and did not show any clear intranight variability. This behavior is quite different from its strong radio intraday variability. Either this object was in a quiescent state in optical in this period, or it is intrinsically not so active in optical as it is in radio regimes.

This paper is the third in the series of papers on M81 globular clusters. In this paper, we present spatial and metal abundance properties of 95 M81 globular clusters, which comprise nearly half of all the M81 globular cluster system. These globular clusters are divided into two M81 metallicity groups by a KMM test. Our results show that, the metal-rich clusters did not demonstrate a centrally concentrated spatial distribution as ones in M31, and metal-poor clusters tend to be less spatially concentrated. In other words, the distribution of the metal-rich clusters in M81 is not very similar to that of M31. Most of the metal-rich clusters distribute at projected radii of 4-8 kpc. It is also noted that the metal-rich clusters distribute within the inner 20 kpc, and the metal-poor ones do out to radii of ~40 kpc. Like our Galaxy and M31, the metallicity distribution of globular clusters in M81 along galactocentric radius suggests that some dissipation occurred during the formation of the globular cluster system, i.e. smooth, pressure-supported collapse models of galaxies are unlikely to produce such radial distribution of metallicity presented in this paper. There is not evident correlation between globular cluster luminosity and metallicity in M81 globular clusters. The overwhelming conclusion of this paper seems to be that a more complete and thorough cluster search is needed in M81.

This paper has been withdrawn by the second author since he is suspected by the first author (the corresponding author) of commiting plagiarism.

A hidden local symmetry formalism with a two-photon counterterm approach is performed based on the relativistic continuum quantum many-body theory. The underlying electromagnetic under-screening as well as screening effects between the electric charged point-like electrons and composite protons are discussed by analyzing the in-medium isospin mixing of Lorentz vector with scalar due to electromagnetic photon. Besides the usual screening results, the main conclusion is that an effective oscillatory instantaneous Coulomb potential between the like-charged collective electrons contributes a very large negative term to the equation of state. This counterintuitive like-charged attraction results from the modulation factor of the opposite charged baryon superfluid background. The anomalous long distance quantum dragging effects between the collective electrons can be induced in a compact Coulomb confinement environment. The physics is of the strongly coupling characteristic in a specific dilute regime.

We present 13-band CCD intermediate-band spectrophotometry of a field centered on the open cluster NGC 7789 from 400 to nearly 1000 nm, taken with Beijing-Arizona-Taiwan-Connecticut (BATC) Multi-Color Survey photometric system. By comparing observed spectral energy distributions of NGC 7789 stars with theoretical ones, the fundamental parameters of this cluster are derived: an age of $1.4\pm0.1$ Gyr, a distance modulus $(m-M)_{0}=11.27\pm0.04$, a reddening $E(B-V)=0.28\pm0.02$, and a metallicity with the solar composition $Z=0.019$. When the surface density profile for member stars with limiting magnitudes of 19.0 in the BATC $e$ band ($\lambda_{\textrm{eff}}=4925$ \AA) is fitted by King model, the core radius $R_{c}=7.52\arcmin$ and the tidal radius $R_{t}=28.84\arcmin$ are derived for NGC 7789. The observed mass function (MF) for main sequence stars of NGC 7789 with masses from 0.95 to 1.85 $M_{\odot}$ is fitted with a power-law function $\phi(m)\propto m^{\alpha}$ and the slope $\alpha=-0.96$ is derived. Strong mass segregation in NGC 7789 is reflected in the significant variation of the concentration parameters $C_{0}=\log (R_{t}/R_{c})$ for member stars of NGC 7789 within different mass ranges: $C_{0}=1.02$ for most of massive stars; $C_{0}=0.37$ for the lowest-mass MS stars. Strong mass segregation in NGC 7789 is also indicated in the significant variation of the slopes $\alpha$ in different spatial regions of the cluster: the MF for stars within the core region has $\alpha=-0.71$, much flatter than that for stars in external regions of the cluster ($\alpha=-1.20$).

By analyzing the Chandra ACIS S3 data we studied the X-ray properties of the low-mass and high-mass X-ray binary populations in the nearby spiral galaxy NGC 5055. A total of 43 X-ray point sources were detected within the 2 effective radii, with 31 sources located on the disk and the rest 12 sources in the bulge. The resolved point sources dominate the total X-ray emission of the galaxy by accounting for about 80% of the total counts in 0.3--10 keV. By carrying out the spectral fittings we calculated the 0.3--10.0 keV luminosities of all the detected X-ray point sources and found that they span a wide range from a few 10^{37} erg s^{-1} to over 10^{39} erg s^{-1}. After compensating for the incompleteness at the low luminosity end, we find that the corrected XLF of the bulge population is well fitted with a broken power-law model with a break at 1.57^{+0.21}_{-0.20}\times 10^{38} erg s^{-1}, while the profile of the disk population's XLF agrees with a single power-law distribution with a slope of 0.93^{+0.07}_{-0.06}. The disk population is significantly richer at {}^{>}_{\sim}2\times10^{38} erg s^{-1} than the bulge population, inferring that the disk may have undergone relatively recent, strong starbursts that significantly increased the HMXB population, although ongoing starbursts are also observed in the nuclear region. Similar XLF profiles of the bulge and disk populations were found in M81. However, in most other spiral galaxies different patterns of the spatial variation of the XLF profiles from the bulge to the disk have been observed, indicating that the star formation and evolution history may be more complex than we have expected.

In recent years, gravitational lensing has been used as a means to detect substructure in galaxy-sized halos, via anomalous flux ratios in quadruply-imaged lenses. In addition to causing anomalous flux ratios, substructure may also perturb the positions of lensed images at observable levels. In this paper, we numerically investigate the scale of such astrometric perturbations using realistic models of substructure distributions. Substructure distributions that project clumps near the Einstein radius of the lens result in perturbations that are the least degenerate with the best-fit smooth macromodel, with residuals at the milliarcsecond scale. Degeneracies between the center of the lens potential and astrometric perturbations suggest that milliarcsecond constraints on the center of the lensing potential boost the observed astrometric perturbations by an order of magnitude compared to leaving the center of the lens as a free parameter. In addition, we discuss methods of substructure detection via astrometric perturbations that avoid full lens modeling in favor of local image observables and also discuss modeling of systems with luminous satellites to constrain the masses of those satellites.

We at first introduce a novel photometric system, which consists of a Schmidt telescope, an objective prism, a CCD camera, and, especially, a multi-peak interference filter. The multi-peak interference filter enables light in multi passbands to pass through it simultaneously. The light in different passbands is differentially refracted by the objective prism and is focused on the CCD separately, so we have multi "images" for each object on the CCD frames. This system enables us to monitor blazars exactly simultaneously in multi wavebands on a single telescope, and to accurately trace the color change during the variation. We used this novel system to monitor the BL Lacertae object S5 0716+714 during 2006 January and February and achieved a very high temporal resolution. The object was very bright and very active during this period. Two strong flares were observed, with variation amplitudes of about 0.8 and 0.6 mags in the $V'$ band, respectively. Strong bluer-when-brighter correlations were found for both internight and intranight variations. No apparent time lag was observed between the $V'$- and $R'$-band variations, and the observed bluer-when-brighter chromatism may be mainly attributed to the larger variation amplitude at shorter wavelength. In addition to the bluer-when-brighter trend, the object also showed a bluer color when it was more active. The observed variability and its color behaviors are consistent with the shock-in-jet model.

An effective Proca Lagrangian action is used to address the vector condensation Lorentz violation effects on the equation of state of the strongly interacting fermions system. The interior quantum fluctuation effects are incorporated as an external field approximation indirectly through a fictive generalized Thomson Problem counterterm background. The general analytical formulas for the $d$-dimensions thermodynamics are given near the unitary limit region. In the non-relativistic limit for $d=3$, the universal dimensionless coefficient $\xi ={4}/{9}$ and energy gap $\Delta/\epsilon_f ={5}/{18}$ are reasonably consistent with the existed theoretical and experimental results. In the unitary limit for $d=2$ and T=0, the universal coefficient can even approach the extreme occasion $\xi=0$ corresponding to the infinite effective fermion mass $m^*=\infty$ which can be mapped to the strongly coupled two-dimensions electrons and is quite similar to the three-dimensions Bose-Einstein Condensation of ideal boson gas. Instead, for $d=1$, the universal coefficient $\xi$ is negative, implying the non-existence of phase transition from superfluidity to normal state. The solutions manifest the quantum Ising universal class characteristic of the strongly coupled unitary fermions gas.

Beijing-Arizona-Taiwan-Connecticut (BATC) multi-band photometric data in the field of open cluster M48 are used to determine its membership. By comparing observed spectral energy distributions (SEDs) of stars with theoretical ones, membership probabilities of 750 stars with limiting magnitude of 15.0 in BATC $c$ band ($\lambda_{eff}=4194$ \AA) are determined. 323 stars with membership probabilities higher than 30% are considered as candidate members of M48. Comparing membership probabilities of 229 common stars obtained by the present method and the proper-motion based methods, a 80% agreement among these methods is obtained.

We present the results of our optical monitoring of the BL Lacertae object OJ 287 during the first half of 2005. The source did not show large-amplitude variations during this period and was in a relatively quiescent state. A possible period of 40 days was derived from its light curves in three BATC wavebands. A bluer-when-brighter chromatism was discovered, which is different from the extremely stable color during the outburst in 1994--96. The different color behaviors imply different variation mechanisms in the two states. We then re-visited the optical data on OJ 287 from the OJ-94 project and found as well a probable period of 40 days in its optical variability during the late-1994 outburst. The results suggest that two components contribute to the variability of OJ 287 during its outburst state. The first component is the normal {\sl blazar} variation. This component has an amplitude similar to that of the quiescent state and also may share a similar periodicity. The second component can be taken as a `low-frequency modulation' to the first component. It may be induced by the interaction of the assumed binary black holes in the center of this object. The 40-day period may be related to the helical structure of the magnetic field at the base of the jet, or to the orbital motion close to the central primary black hole.

We use the Sloan Digital Sky Survey (SDSS) spectroscopic sample to constrain the projected radial distribution of satellites around isolated ~ L* galaxies. We employ mock galaxy catalogs derived from high-resolution cosmological simulations to investigate the effects of interloper contamination and show that interlopers significantly bias the estimated slope of the projected radial distribution of satellites. We also show that the distribution of interlopers around galaxies is expected to be non-uniform in velocity space because galaxies are clustered and reside in crowded environments. Successful methods of interloper contamination correction should therefore take into account environments of the host galaxies. Two such new methods are presented and the most reliable of them is used to correct for interloper contamination in analyses of the SDSS galaxy sample. The best fit power-law slope of the interloper-corrected surface density distribution of satellites, Sigma(R) ~ R^alpha, in the volume-limited SDSS sample is alpha ~ -1.7 +/- 0.1, independent of the galaxy and satellite luminosities. Comparison with Lambda-CDM simulations shows that the radial distribution of the SDSS satellites is more concentrated than that of subhalos around galaxy-sized halos, especially at R < 100 kpc/h. The predicted dark matter radial distribution is somewhat more concentrated than the profile of the SDSS satellites, but the difference is not statistically significant for our sample.

We independently redetermine the reddening and age of the globular cluster 037-B327 in M31 by comparing independently obtained multicolour photometry with theoretical stellar population synthesis models. 037-B327 has long been known to have a very large reddening value, which we confirm to be E(B-V)=1.360+-0.013, in good agreement with the previous results. We redetermine its most likely age at 12.4+- 3.2 Gyr. 037-B327 is a prime example of an unusually bright early counterpart to the ubiquitous ``super'' star clusters presently observed in most high-intensity star-forming regions in the local Universe. In order to have survived for a Hubble time, we conclude that its stellar IMF cannot have been top-heavy. Using this constraint, and a variety of SSP models, we determine a photometric mass of M_{GC}=(3.0+-0.5)*10^7, somewhat depending on the SSP models used, the metallicity and age adopted and the IMF representation. This mass, and its relatively small uncertainties, make this object the most massive star cluster of any age in the Local Group. Assuming that the photometric mass estimate thus derived is fairly close to its dynamical mass, we predict that this GC has a (one-dimensional) velocity dispersion of order (72+-13)km/s. As a surviving "super" star cluster, this object is of prime importance for theories aimed at describing massive star cluster evolution.

We present photometric data of the type II-P supernova (SN) 2004dj in NGC 2403. The multicolor light curves cover the SN from $\sim$ 60 to 200 days after explosion, and are measured with a set of intermediate-band filters that have the advantage of tracing the strength variations of some spectral features. The light curves show a flat evolution in the middle of the plateau phase, then decline exponentially at the late times, with a rate of 0.10$\pm$0.03 mag (10 days)$^{-1}$ in most of the filters. In the nebular phase, the spectral energy distribution (SED) of SN 2004dj shows a steady increase in the flux near 6600 \AA and 8500 \AA, which may correspond to the emission lines of H$\alpha$ and Ca II near-IR triplet, respectively. The photometric behavior suggests that SN 2004dj is a normal SN II-P. Compared with the light curves of another typical SN II-P 1999em, we estimate the explosion date to be June 10$\pm$21 UT, 2004 (JD 2453167$\pm$21) for SN 2004dj. We also estimate the ejected nickel mass during the explosion to be $M(^{56}\rm{Ni})$ = 0.023 $\pm$ 0.005 $M_{\odot}$ from two different methods, which is typical for a SN II-P. We derive the explosion energy $E \approx 0.75^{+0.56}_{-0.38}\times10^{51}$ erg, the ejecta mass $M \approx 10.0^{+7.4}_{-5.2}$ $M_{\odot}$, and the initial radius $R \approx 282^{+253}_{-122}$ $R_{\odot}$ for the presupernova star of SN 2004dj, which are consistent with other typical SNe II-P.

In this paper, we give the spectral energy distributions of 42 M81 globular clusters in 13 intermediate-band filters from 4000 to 10000 A, using the CCD images of M81 observed as part of the BATC multicolor survey of the Sky. The BATC multicolor filter system is specifically designed to exclude most of the bright and variable night-sky emission lines including the OH forest. Hence, it can present accurate SEDs of the observed objects. These spectral energy distributions are low-resolution spectra, and can reflect the stellar populations of the globular clusters. This paper confirms the conclusions of Schroder et al. that, M81 contains clusters as young as a few Gyrs, which also were observed in both M31 and M33

This paper presents accurate spectral energy distributions (SEDs) of 16 M31 globular clusters (GCs) confirmed by spectroscopy and/or high spatial-resolution imaging, as well as 30 M31 globular cluster candidates detected by Mochejska et al. Most of these candidates have m_V > 18, deeper than previous searches, and these candidates have not yet been confirmed to be globular clusters. The SEDs of these clusters and candidates are obtained as part of the BATC Multicolor Survey of the Sky, in which the spectrophotometrically-calibrated CCD images of M31 in 13 intermediate-band filters from 4000 to 10000 A were observed. These filters are specifically designed to exclude most of the bright and variable night-sky emission lines including the OH forest. In comparison to the SEDs of true GCs, we find that some of the candidate objects are not GCs in M31. SED fits show that theoretical simple stellar population (SSP) models can fit the true GCs very well. We estimate the ages of these GCs by comparing with SSP models. We find that, the M31 clusters range in age from a few ten Myr to a few Gyr old, as well as old GCs, confirming the conclusion that has been found by Barmby et a, Williams & Hodge, Beasley et al., Burstein et al. and Puzia et al. in their investigations of the SEDs of M31 globular clusters.

037-B327 is of interest because it is both the most luminous and the most highly reddened cluster known in M31. Deep observations with the Advanced Camera for Surveys on the Hubble Space Telescope provide photometric data in the F606W band, and also show that this cluster is crossed by a dust lane. We determined the structural parameters of 037-B327 by fitting the observed surface brightness distribution to a King model with r_c=0.72" (=2.69 pc), and r_t=5.87" (=21.93 pc), and a concentration index c=log(r_t/r_c)=0.91. The surface brightness profile appears to be essentially flat within 0.25" of the center and shows no signs of core collapse. Although the dust lane affects the photometry, the King model fits the surface brightness profile well except for the regions badly affected by the dust lane. We also calculate the half-light radius r_h=1.11" (=4.15 pc). Combined with previous photometry, we find that this object falls in the same region of the M_V versus log R_h diagram as do Omega Centauri, M54 and NGC 2419 in the Milky Way and the massive cluster G1 in M31. All four of these objects have been claimed to be the stripped cores of former dwarf galaxies. This suggests that 037-B327 may also be the stripped core of a former dwarf companion to M31.

We study the statistical properties of magnification perturbations by substructures in strong lensed systems using linear perturbation theory and an analytical substructure model including tidal truncation and a continuous substructure mass spectrum. We demonstrate that magnification perturbations are dominated by perturbers found within a tidal radius of an image, and that sizable magnification perturbations may arise from small, coherent contributions from several substructures within the lens halo. We find that the root-mean-square (rms) fluctuation of the magnification perturbation is 10% to 20% and both the average and rms perturbations are sensitive to the mass spectrum and density profile of the perturbers. Interestingly, we find that relative to a smooth model of the same mass, the average magnification in clumpy models is lower (higher) than that in smooth models for positive (negative) parity images. This is opposite from what is observed if one assumes that the image magnification predicted by the best-fit smooth model of a lens is a good proxy for what the observed magnification would have been if substructures were absent. While it is possible for this discrepancy to be resolved via nonlinear perturbers, we argue that a more likely explanation is that the assumption that the best-fit lens model is a good proxy for the magnification in the absence of substructure is not correct. We conclude that a better theoretical understanding of the predicted statistical properties of magnification perturbations by CDM substructure is needed in order to affirm that CDM substructures have been unambiguously detected.

Using the High Dispersion Spectrograph at the Subaru Telescope, echelle spectra of two giant arcs, i.e. nebulosities Cw and L associated with HH 131 in Orion are presented. Typical emission lines of Herbig-Haro objects have been detected towards Cw. With the 2.16 m telescope of National Astronomical Observatories, spectra of Nebu. C, L and K are obtained, which also show strong [SII]6717/6731, H$\alpha$ and [NII]6583 emission lines. Position-velocity distributions of Cw and L are analyzed. The fastest radial velocity of Cw is V_r ~ -18.0 km/s. When the flow at L goes to the south, it slows down. The fastest radial velocity of L has been observed of -45.0 km/s and the slowest value is about -18.3 km/s. The similarity of the velocities and their positional connection indicate that Cw and L are physically associated. The entire flow tends to become less excited and less ionized when going further to the south (i.e., from Nebu. K, L to C). The electron densities of all the observed nebulosities are low (n_e ~ 10^2 cm^-3). Double kinematic signatures have been found in Cw from its [NII]6583 profiles while the observed H$\alpha$ profiles of Cw are almost symmetric. Bow shock models appear to agree with the observed position-velocity diagrams of the [NII spectra better than H$\alpha$ spectra. With the suggestion that these arcs are HH shocks possibly ejected out of the Orion A molecular cloud by an uncertain source, their spectra show low to intermediate excitation from their diagnostic line ratios.

This paper includes two parts. The first is to present the spectral energy distributions (SEDs) of 49 globular cluster (GC) X-ray sources in the BATC 13 intermediate-band filters from 3800 to 10000 A, and identify 8 unidentified X-ray sources in M31. Using the X-ray data of Einstein observation from 1979 to 1980, ROSAT HRI observation in 1990, Chandra HRC and ACIS-I observations from 1999 to 2001, and the BATC optical survey from 1995 to 1999, we find 49 GC X-ray sources and 8 new unidentified X-ray sources in the BATC M31 field. By analyzing SEDs and FWHMs, 4 of the 8 X-ray sources may be GC candidates. The second is to present some statistical relationships about 62 GC X-ray sources, of which 58 are already known, and 4 are identified in this paper. The distribution of M31 GC X-ray sources' V mags is bimodal, with peaks at m_v = 15.65 and m_v = 17.89, which is different from the distribution of GC candidates. The distribution of B-V color shows that,the GC X-ray sources seem to be associated preferentially with the redder GCs, in agreement with the previous results. Kolmogorov-Smirnov test shows that the maximum value of the absolute difference of B-V distributions of GC X-ray sources and GCs is D_{max}=0.181, and the probability P=0.068 which means we can reject the hypothesis that the two distributions are the same at the 90.0% confidence level. In the end, we study the correlation between X-ray luminosity (0.3-10 keV) and the optical luminosity (in V band) of the GC X-ray sources in M31, and find that there exits a weak relationship with the linear correlation coefficient r = 0.36 at the confidence level of 98.0%.

If our universe underwent inflation, its entropy during the inflationary phase was substantially lower than it is today. Because a low-entropy state is less likely to be chosen randomly than a high-entropy one, inflation is unlikely to arise through randomly-chosen initial conditions. To resolve this puzzle, we examine the notion of a natural state for the universe, and argue that it is a nearly-empty spacetime. If empty space has a small vacuum energy, however, inflation can begin spontaneously in this background. This scenario explains why a universe like ours is likely to have begun via a period of inflation, and also provides an origin for the cosmological arrow of time.

Three deep flux minima were observed with nearly the same time-scales and intervals for the blazar PKS 1510-089 in the past few years. A binary black hole system was proposed to be at the nucleus of this object, and a new minimum was predicted to occur in 2002 March. We monitored this source with a 60/90 cm Schmidt telescope from 2002 February to April. In combination with the data obtained by Xie et al. (2004) in the same period, we presented for the 2002 minimum a nearly symmetric light curve, which would be required by an eclipsing model of a binary black hole system. We also constrained the time-scale of the minimum to be 35 min, which is more consistent with the time-scales ~42 min of the three previous minima than the 89 min time-scale given by the same authors. The wiggling miniarcsecond radio jet observed in this object is taken as a further evidence for the binary black hole system. The `coupling' of the periodicity in light curve and the helicity in radio jet is discussed in the framework of a binary black hole system.

Using the theoretical stellar population synthesis models of BC96, Kong et al. (2003) showed that some BATC colors and color indices could be used to disentangle the age and metallicity effect. They found that there is a very good relation between the flux ratio of L_{8510}/L_{9170} and the metallicity for stellar populations older than 1 Gyr. In this paper, based on the Kong et al. results and on the multicolor spectrophotometry of Ma et al. (2001, 2002a,b,c), we estimate the metallicities of 31 old star clusters in the nearby spiral galaxy M33, 23 of which are ``true'' globular clusters. The results show that most of these old clusters are metal poor. We also find that the ages and metal abundance for these old star clusters of M33 do not vary with deprojected radial position.

In this paper we present catalogs of photometric and spectroscopic data for M81 globular clusters (GCs). The catalogs include B- and V- photometric and reddening data of 95 GCs, and spectroscopic metallicities of 40 GCs in M81. Using these data, we make some statistical correlations. The results show that the distributions of intrinsic B and V colors and metallicities are bimodal, with metallicity peaks at [Fe/H]=-1.45 and -0.53, respectively as has been demonstrated for our Milky Way and M31. The relation between spectroscopic metallicity and intrinsic B and V color also exists as it does for the Milky Way and M31.

The electromagnetic(EM) interactions between charged protons on the correlations of nucleons are discussed by introducing the Anderson-Higgs mechanism of broken U(1) EM symmetry into the relativistic nuclear theory with a parametric photon mass. The non-saturating Coulomb force contribution is emphasized on the equation of state of nuclear matter with charge symmetry breaking(CSB) at finite temperature and the breached $^1S_0$ pairing correlations of proton-proton and neutron-neutron. The universal properties given by an order parameter field with a non-zero vacuum expectation value (VEV) nearby phase transition are explored within the mean field theory(MFT) level. This mechanism can be extended to the charged or charge neutralized strongly coupling multi-components system for the discussion of binding or pairing issues.

Optical monitoring of S5 0716+714 was performed with a 60/90 Schmidt telescope in 2003 November and December and 2004 January for studying the variability of the object on short timescales. Due to the high brightness of the source we could carry out quasi-simultaneous measurements in three bands with a temporal resolution of about 20 minutes by using one single telescope. Intraday and intranight variations were observed with an overall change of ~0.9 mag during the whole campaign. Two outbursts were recorded on JD 2453005 and JD 2453009. Minimum timescales of a few hours were derived from the light curves of individual nights but were different from night to night. A bluer-when-brighter chromatism was present when the object was in fast flare, but was absent when it was in a relatively quiescent state. Our results are basically consistent with the shock-in-jet model and demonstrate that the geometrical effects can sometimes play an important role in the variability of blazars.

This paper presents 13 color CCD intermediate-band spectrophotometry of a field centered on the open cluster M48 (NGC 2548), from 400nm to nearly 1000nm, taken with Beijing-Arizona-Taiwan-Connecticut (BATC) Multi-Color Survey Photometric system. The fundamental parameters of this cluster are derived with a new method which based on the comparison between the spectral energy distributions (SEDs) of cluster stars and the theoretical SEDs of Padova models. We find that the best-fitting age of M48 is 0.32 Gyr, the distance is 780 pc and a reddening E(B-V)=0.04 with a solar metallicity Z=0.019

We suggest that spontaneous eternal inflation can provide a natural explanation for the thermodynamic arrow of time, and discuss the underlying assumptions and consequences of this view. In the absence of inflation, we argue that systems coupled to gravity usually evolve asymptotically to the vacuum, which is the only natural state in a thermodynamic sense. In the presence of a small positive vacuum energy and an appropriate inflaton field, the de Sitter vacuum is unstable to the spontaneous onset of inflation at a higher energy scale. Starting from de Sitter, inflation can increase the total entropy of the universe without bound, creating universes similar to ours in the process. An important consequence of this picture is that inflation occurs asymptotically both forwards and backwards in time, implying a universe that is (statistically) time-symmetric on ultra-large scales.

We present Chandra X-ray observations of galaxy cluster A168 (z=0.045). Two X-ray peaks with a projected distance of 676 kpc are found to be located close to two dominant galaxies, respectively. Both peaks are significantly offset from the peak of the number density distribution of galaxies. This suggests that A168 consists of two subclusters, a northern subcluster (A168N) and a southern subcluster (A168S). Further X-ray imaging analysis reveals that (1) the X-ray isophotes surrounding the two X-ray peaks are heavily distorted, (2) an elongated and ontinuous filament connects the two X-ray peaks. These suggest that strong interactions have occurred between the two subclusters. Spectral analysis shows that A168 has a mean temperature of 2.53 +/- 0.09 keV and a mean metallicity of 0.31 +/- 0.04 Z_{solar}. The metallicity is roughly a constant across the cluster but the temperature shows some systematic variations. Most X-ray, optical and radio properties of A168 are consistent with it being an off-axis merger several Gyrs after a core passage, although detailed numerical simulations are required to see whether the observed properties, in particular the significant offset between the optical and X-ray centers, can be reproduced in such a scenario.

We present optical candidates for 75 X-ray sources in a $\sim 1$ deg$^2$ overlapping region with the medium deep ROSAT survey. These candidates are selected using the multi-color CCD imaging observations made for the T329 field of the Beijing-Arizona-Taipei-Connecticut (BATC) Sky Survey. These X-ray sources are relatively faint (CR $<< 0.2 s^{-1}$) and thus mostly are not included in the RBS catalog, they also remain as X-ray sources without optical candidates in a previous identification program carried out by the Hamburg Quasar Survey. Within their position-error circles, almost all the X-ray sources are observed to have one or more spatially associated optical candidates within them down to the magnitude $m_V \sim 23.1$. We have classified 149 of 156 detected optical candidates with 73 of the 75 X-ray sources with a SED-based Object Classification Approach (SOCA). These optical candidates include: 31 QSOs, 39 stars, 37 starburst galaxies, 42 galaxies, and 7 "just" visible objects. We have also cross-correlated the positions of these optical objects with NED, the FIRST radio source catalog and the 2MASS catalog. Separately, we have also SED-classified the remaining 6011 objects in our field of view. Optical objects are found at the $6.5\sigma$ level above what one would expect from a random distribution, only QSOs are over-represented in these error circles at greater than 4$\sigma$ frequency. We estimate redshifts for all extragalactic objects, and find a good correspondence of our predicted redshift with the measured redshift (a mean error of 0.04 in $\Delta z$. There appears to be a supercluster at z $\sim$ 0.3-0.35 in this direction, including many of the galaxies in the X-ray error circles are found in this redshift range.

The $^1 S_0$ pairing of symmetric nuclear matter is discussed in the frame work of relativistic nuclear theory with Dyson-Schwinger equations (DSEs). The in-medium nucleon and meson propagators are treated in a more self-consistent way through meson polarizations. The screening effects on mesons due to in-medium nucleon excitation are found to reduce the $^1S_0$ pairing gap and shift remarkably the gap peak to low density region.

The influence of the in-medium mesons on the effective nucleon mass and in turn on the equation of state of hot/dense nuclear matter is discussed in the Walecka model. Due to the self-consistent treatment of couplings between nucleons and $\sigma $ and $\omega$ mesons, the temperature and density dependence of the effective hadron masses approaches more towards the Brown-Rho scaling law, and the compression modulus $K$ is reduced from $550 MeV$ in mean field theory to an accepted value $318.2 MeV$.

This paper presents a multicolor study of the nearby cluster of galaxies Abell 168 (z=0.045) with 13 intermediate-band filters in the Beijing-Arizona-Taiwan-Connecticut (BATC) filter system. After a cross-identification between the photometric data obtained from the BATC and the Sloan Digital Sky Survey (SDSS), a catalog containing 1553 galaxies down to r'<20.0 mag is achieved, which includes 121 spectroscopically confirmed member galaxies. The technique of photometric redshift has been applied to all these galaxies with combined 18-band (13 from BATC and 5 from SDSS) spectral energy distributions (SEDs), in order to perform a faint membership selection in Abell 168. As a result, 255 galaxies are newly selected as the member candidates. Based on the enlarged sample of cluster galaxies, the spatial distribution and dynamics of Abell 168 are investigated. In the light of the spatial distribution of the member galaxies and the 0.2-3.5 keV X-ray image by Einstein observatory, it seems that Abell 168 consists of two merging subclusters with a relative radial velocity of 264 +/- 142 km/s. With the help of ROSTAT software, a detailed investigation of the dynamics shows the intrinsic difference in the velocity distributions for these two subclusters. Under a linear two-body model, they are found to be a gravitationally bound system with 92% probability. The slight deviation of the local velocity distribution from the overall distribution in the central region may suggest a picture of head-on collision in the projection plane. By using a double Schechter function, the close examination of the luminosity functions reveals the different galaxy content of these two subclusters, which implies the merging process is possibly at the early stage.

In a first of a series of studies of the H-alpha + [NII] emission from nearby spiral galaxies, we present measurements of H-alpha + [NII] emission from HII regions in M81. Our method uses large-field-CCD images and long-slit spectra, and is part of the ongoing Beijing-Arizona-Taipei-Connecticut Sky Survey. The CCD images are taken with the NAOC 0.6/0.9m f/3 Schmidt telescope at the Xinglong Observing Station, using a multicolor filter set. Spectra of 10 of the brightest HII regions are obtained using the NAOC 2.16m telescope with a Tek 1024 X 1024 CCD. The continua of the spectra are calibrated by flux-calibrated images taken from the Schmidt observations. We determine the continuum component of our H-alpha + [NII] image via interpolation from the more accurately-measured backgrounds (M81 starlight) obtained from the two neighboring (in wavelength) BATC filter images. We use the calibrated fluxes of H-alpha + [NII] emission from the spectra to normalize this interpolated, continuum-subtracted H-alpha + [NII] image. We estimate the zero point uncertainty of the measured H-alpha + [NII] emission flux to be $\sim$ 8%. A catalogue of H-alpha + [NII] fluxes for 456 HII regions is provided, with those fluxes are on a more consistent linear scale than previously available. The logarithmically-binned H-alpha + [NII] luminosity function of HII regions is found to have slope $\alpha$ = -0.70, consistent with previous results (which allowed $\alpha=-0.5 \sim -0.8$). From the overall H-alpha + [NII] luminosity of the HII regions, the star formation rate of M81 is found to be $\sim 0.68 M_{\odot} {\rm yr}^{-1}$, modulo uncertainty with extinction corrections.

We report the first results of star counts carried out with the National Astronomical Observatories (NAOC) 60/90 cm Schmidt Telescope in 15 intermediate-band filters from 3000 to 10000 {\AA} in the BATC survey. We analyze a sample of over 1400 main sequence stars ($14\le$V$\le21$), which lie in the field with central coordinates R.A.=$09^h53^m13^s.30$ and DEC=47$^\circ49^{\prime}00^{\prime\prime}.0$ (J2000). The field of view is 0.95 deg$^{2}$, and the spatial scale was $1\arcsec.67$. In our model, the distribution of stars perpendicular to the plane of the Galaxy is given by two exponential disks (thin disk plus thick disk) and a de Vaucouleurs halo. Based on star counts, we derive the scale heights of the thin disk to be $320^{+14}_{-15}$ pc and of the thick disk to be $640^{+30}_{-32}$ pc, respectively, with a local density of $7.0\pm1%$ of the thin disk. We find that the observed counts support an axial ratio of $c/a\le0.6$ for a de Vaucouleurs $r^{1/4}$ law, implying a more flattened halo. We also derive the stellar luminosity function (SLF) for the thin disk, and it partly agrees with the Hipparcos luminosity function.

Multiply-imaged quasar lenses can be used to constrain the substructure mass fraction in galaxy-sized dark matter halos via anomalous flux ratios in lensed images. The flux ratios, however, can be affected by both the substructure in the lens halo and by isolated small-mass halos along the entire line-of-sight to the lensed source. While lensing by dark matter clumps near the lens galaxy is more efficient than elsewhere, the cumulative effect of all objects along the line-of-sight could be significant. Here we estimate the potential contribution of isolated clumps to the substructure lensing signal using a simple model motivated by cosmological simulations. We find that the contribution of isolated clumps to the total lensing optical depth ranges from a few to tens percent, depending on assumptions and the particular configuration of the lens. Therefore, although the contribution of isolated clumps to the lensing signal is not dominant, it should not be neglected in detailed analyses of substructure lensing. For the currently favored Lambda-CDM model, the total calculated optical depth for lensing is high, tau ~ 0.2-20 and could, therefore, naturally explain the high frequency of anomalous flux ratios in observed lenses. The prediction, however, is highly sensitive to the spatial distribution of substructure halos in the innermost regions of the lens halo, which is still very uncertain. Therefore, constraints on the properties of the substructure population or accurate cosmological constraints, such as the mass of the warm dark matter particle, are difficult if not impossible to derive at this point.

This paper is the fourth in a series presenting spectrophotometry of 51 globular cluster candidates, that were detected by Mochejska et al. in the nearby galaxy M33 using the data collected by the DIRECT project. The frames of M33 in this study were taken as part of the BATC Multicolor Sky Survey. We obtained the spectral energy distributions of these candidates in 13 intermediate-band filters. By comparing the integrated photometric measurements with theoretical stellar population synthesis models of Bruzual & Charlot, we estimated their ages. The BC96 models provide the evolution in time of the spectrophotometric properties of simple stellar populations for a wide range of stellar metallicity. Our results show that half of the candidates are younger than 10^8 years, whose age degeneracy is not pronounced. We also find that globular clusters formed continuously in M33 from ~4*10^6 -- 10^{10} years. Our results are in agreement with Chandar et al., who estimated ages for 35 globular clusters candidates in common by comparing the photometric measurements to integrated colors from theoretical models by Bertelli et al. The Kolmogorov-Smirnov Test shows that the maximum value of the absolute difference of estimated ages between Chandar et al. and us is 0.48, and the significance level probability is 100.00 per cent.

In this paper we present CCD multicolor photometry for 172 globular clusters (GCs), taken from the Bologna catalog (Battistini et al. 1987), in the nearby spiral galaxy M31. The observations were carried out by using the National Astronomical Observatories 60/90 cm Schmidt Telescope in 13 intermediate-band filters, which covered a range of wavelength from 3800 to 10000A. This provides a multicolor map of M31 in pixels of 1.7"*1.7". By aperture photometry, we obtain the spectral energy distributions (SEDs) for these GCs. Using the relationship between the BATC intermediate-band system used for the observations and the UBVRI broad-band system, the magnitudes in the B and V bands are derived. The computed V and B-V are in agreement with the values given by Battistini et al. (1987) and Barmby et al. (2000). Finally, by comparing the photometry of each GC with theoretical stellar population synthesis models, we estimate ages of the sample GCs for different metallicities. The results show that nearly all our sample GCs have ages more than 10^{9} years, and most of them are around 10^{10} years old. Also, we find that GCs fitted by the metal-poor model are generally older than ones fitted by the metal-rich model.

We discuss how to reduce theoretical uncertainties in the neutrino-deuteron breakup cross-sections crucial to the Sudbury Neutrino Observatory's efforts to measure the solar neutrino flux. In effective field theory, the dominant uncertainties in all neutrino-deuteron reactions can be expressed through a single, common, isovector axial two-body current parameterized by $L_{1,A}$. After briefly reviewing the status of fixing $L_{1,A}$ experimentally, we present a constraint on $L_{1,A}$ imposed by existing reactor antineutrino-deuteron breakup data. This constraint alone leads to an uncertainty of 6-7% at 7 MeV neutrino energy in the cross-sections relevant to the Sudbury Neutrino Observatory. However, more significantly for the Sudbury experiment, the constraint implies an uncertainty of only 0.7% in the ratio of charged to neutral current cross-sections used to verify the existence of neutrino oscillations, at the same energy. This is the only direct experimental constraint from the two-body system, to date, of the uncertainty in these cross-sections.

In this paper, we present the methodology of photometric redshift determination with the BATC 15-color system by using hyperz program. Both simulated galaxies and real galaxies with known redshifts were used to estimate the accuracy of redshifts inferred from the multicolor photometry. From the test with simulated galaxies, the uncertainty in the inferred redshifts is about $0.02\sim0.03$ for a given range of photometric uncertainty of $0.05 \sim 0.10$. The results with the 27 real galaxies are in good agreement with the simulated ones. The advantage of using BATC intermediate-band system to derive redshift is clear through the comparison with the UBVRI broad-band system. The accuracy in redshift determination with BATC system is mainly affected by the selection of filters and the photometric uncertainties in the observation. When we take the limiting magnitudes of the 15 filters into account, we find that redshift can be determined with good accuracy for galaxies with redshifts less than 0.5, using only filters with central wavelengths shorter than 6000 A.

In this paper, we present a set of photometric observations in 15 colors of stars in the Landolt covered by the system range from 300 nm to 1000 nm. Visual magnitudes of the stars being studied in the field are from 10th to 20th mag. The observational methodology and the data reduction procedures are described. The relationships between the BATC intermediate-band system and the Landolt UBVRI broad band system are obtained. A catalogue of the photometry has been produced which contains the SEDs of 3613 stars. The electronic form of this catalogue can be accessed at the CDS via anonymous ftp to cdsarc.u-strasbg.fr.

This paper is the first in a series presenting CCD multicolor photometry for 145 HII regions, selected from 369 candidate regions from Boulesteix et al. (1974), in the nearby spiral galaxy M33. The observations, which covered the whole area of M33, were carried out by the Beijing Astronomical Observatory 60/90 cm Schmidt Telescope, in 13 intermediate-band filters, covering a range of wavelength from 3800 to 10000 A. This provides a series of maps which can be converted to a multicolor map of M33, in pixels of 1.7''*1.7''. Using aperture photometry we obtain the spectral energy distributions (SEDs) for these HII regions. We also give their identification charts. Using the relationship between the BATC intermediate-band system used for the observations and the UBVRI broad-band system, the magnitudes in the B and V bands are then derived. Histograms of the magnitudes in V and in B-V are plotted, and the color-magnitude diagram is also given. The distribution of magnitudes in the V band shows that the apparent magnitude of almost all the regions is brighter than 18, corresponding to an absolute magnitude of -6.62 for an assumed distance modulus of 24.62, which corresponds to a single main sequence O5 star, while the distribution of color shows that the sample is blue, with a mode close to -0.05 as would be expected from a range of typical young clusters.

In this third paper of our series, we present CCD spectrophotometry of 78 star clusters that were detected by Chandar, Bianchi, & Ford in the nearby spiral galaxy M33. CCD images of M33 were obtained as a part of the BATC Color Survey of the sky in 13 intermediate-band filters from 3800 to 10000{\AA}. By aperture photometry, we obtain the spectral energy distributions of these 78 star clusters. As Chandar, Bianchi, & Ford did, we estimate the ages of our sample clusters by comparing the photometry of each object with theoretical stellar population synthesis models for different values of metallicity. We find that the sample clusters formed continuously in M33 from $\sim 3\times10^6$ -- $10^{10}$ years. This conclusion is consistent with Chandar, Bianchi, & Ford. The results also show that, there are two peaks in cluster formation, at $\sim 8\times10^6$ and $\sim 10^9$ years in these clusters.

In this paper we study the properties of ten halo globular clusters in the nearby spiral galaxy M33. CCD images of M33 were obtained as a part of the BATC Colour Survey of the sky in 13 intermediate-band filters from 3800 to 10000{\AA}. By aperture photometry, we obtain the spectral energy distributions (SEDs) of these globular clusters. We estimate the ages of our sample clusters by comparing the photometry of each object with theoretical stellar population synthesis models for different values of metallicity. Our results suggest that eight of the ten sample halo globular clusters have ``intermediate'' ages, i.e. between 1 and 8 Gyrs.

An optical photometric observation with the BATC multicolor system is carried out for the central 56' x 56' region of the nearby cluster of galaxies, Abell 2634. We achieved the spectral energy distributions (SEDs) of 5572 sources detected down to $V \sim 20$ mag, including 178 previously known galaxies, with fourteen filters covering a wavelength range from 3600 \AA\ to 10000 \AA. After excluding the foreground and background galaxies, a sample of 124 known members is formed for an investigation of the SED properties. Based on the knowledge of SED properties of member galaxies, we performed the selection of faint galaxies belonging to Abell 2634. The color-color diagrams are powerful in the star/galaxy separation, and 359 faint galaxies are selected by their color features. The technique of photometric redshift and color-magnitude correlation for the early-type galaxies are applied for these faint galaxies, and a list of 74 faint member galaxies is achieved. Basis on the enlarged sample of member galaxies, the spatial distribution and color-magnitude relation of the galaxies in core region of Abell 2634 are discussed. We find a tendency that the color index dispersion of the early-type members is larger for the outer region, which might reflect some clues about the environmental effect on the evolution of galaxies in a cluster.

We present new interferometric HI and CO (1-0) observations of the luminous infrared source, Arp 256. Arp 256 consists of two spiral galaxies in an early stage of merging, with a projected nuclear separation of 29 kpc (54") and an infrared luminosity of 2.0E11 L_sun. Despite the large separation of the galaxies' nuclei and mildly disrupted stellar components, the HI disks are found to be strongly disrupted, and the southern galaxy in Arp 256 shows an elevated star formation efficiency, which is consistent with a nuclear starburst. Both of these results run contrary to expectations, posing interesting questions on the physical mechanisms involved in stimulating star formation during an interaction.

In this paper we present CCD spectrophotometry of the star clusters that were detected by Chandar, Bianchi, & Ford in the nearby spiral galaxy M33, using the images obtained with the Beijing Astronomical Observatory 60/90 cm Schmidt Telescope in 13 intermediate-band filters from 3800 to 10000{\AA}. The observations cover the whole area of M33 with a total integration of 32.75 hours from September 23, 1995 to August 28, 2000. This provides a multi-color map of M33 in pixels of $1.7'' \times 1.7''$. By aperture photometry, we obtain the spetral energy distributions (SEDs) of these star clusters. Using theoretical stellar population synthesis models, we also obtain the distributions of age and metallicity of of these star clusters. These clusters formed continuously from $\sim 3\times10^6$ -- $10^{10}$ years, and have a large span of metallicity from ${\rm {Z}=0.0004}$ to ${\rm {Z}=0.05}$.

We present the analysis of archival data from ROSAT and ASCA of a serendipitous BL Lac object PKS 2316-423. Because of its featureless non-thermal radio/optical continuum, PKS 2316-423 has been called as a BL Lac candidate in the literature. PKS 2316-423 was evidently variable over the multiple X-ray observations, in particular, a variable high-energy tail of the synchrotron radiation is revealed. The X-ray spectral analysis provides further evidence of the synchrotron nature of its broad-band spectrum: a steep and downward curving spectrum between 0.1--10 keV, typical of high-energy peaked BL Lacs (HBL). The spectral energy distribution (SED) through radio-to-X-ray yields the synchrotron radiation peak at frequency Vp=7.3 10^{15} Hz, with integrated luminosity of Lsyn=2.1 10^{44} ergs s^{-1}. The averaged SED properties of PKS 2316-423 are very similar to those ``intermediate'' BL Lac objects (IBL) found recently in several deep surveys, such as Deep X-ray Radio Blazar, Radio-Emitting X-ray, and ROSAT-Green Bank surveys. We suggest that PKS 2316-423 is an IBL though it also shows some general features of a HBL. Actually, this double attribute of PKS 2316-423 provides a good test of the prediction that an IBL object can show either synchrotron or inverse-Compton characteristics in different variability states.

In this paper, we present a multi-color photometric study of the nearby spiral galaxy M81, using images obtained with the Beijing Astronomical Observatory 60/90 cm Schmidt Telescope in 13 intermediate-band filters from 3800 to 10000{\AA}. The observations cover the whole area of M81 with a total integration of 51 hours from February 1995 to February 1997. This provides a multi-color map of M81 in pixels of $1\arcsec.7 \times 1\arcsec.7$. Using theoretical stellar population synthesis models, we demonstrate that some BATC colors and color indices can be used to disentangle the age and metallicity effect. We compare in detail the observed properties of M81 with the predictions from population synthesis models and quantify the relative chemical abundance, age and reddening distributions for different components of M81. We find that the metallicity of M81 is about $Z=0.03$ with no significant difference over the whole galaxy. In contrast, an age gradient is found between stellar populations of the central regions and of the bulge and disk regions of M81: the stellar population in its central regions is older than 8 Gyr while the disk stars are considerably younger, $\sim 2$ Gyr. We also give the reddening distribution in M81. Some dust lanes are found in the galaxy bulge region and the reddening in the outer disk is higher than that in the central regions.

The cross section for n+p -> d+gamma is calculated at energies relevant to big bang nucleosynthesis using the recently developed effective field theory that describes the two-nucleon sector. The E1 amplitude is computed up to NNNLO and depends only upon nucleon-nucleon phase shift data. In contrast, the M1 contribution is computed up to NLO, and the four-nucleon-one-magnetic-photon counterterm that enters is determined by the cross section for cold neutron capture. The uncertainty in the calculation for nucleon energies up to E ~ 1 MeV is estimated to be < 4%.

The purpose of this paper is to present a method to self-calibrate the spectral energy distribution (SED) of objects in a survey based on the fitting of an SED library to the observed multi-color photometry. We adopt for illustrative purposes the Vilnius (Strizyz and Sviderskiene 1972) and Gunn & Stryker (1983) SED libraries. The self-calibration technique can improve the quality of observations which are not taken under perfectly photometric conditions. The more passbands used for the photometry, the better the results. This technique has been applied to the BATC 15-passband CCD survey.

We obtained very deep exposures of NGC 5907 with a Schmidt telescope, large-format CCD, and intermediate-band filters centered at 6660A and 8020A. These two filters, part of a 15-filter set, are custom-designed to avoid the brightest (and most variable) night sky lines, with the result that our images go deeper, with lower sky noise, than those taken with broad-band filters. Our errors of observation reach 1 mag/arcsec^2 at 29.00 mag/arcsec-2 in the 6660A image (= 28.7 in R-band), and 27.4 mag/arcsec^ in the 8020A image (same zero point as I-band). Fainter than R = 27 mag arcsec^, the surface brightness around NGC 5907 is strongly asymmetric, being mostly brighter the NW side of the galaxy midplane. This asymmetry rules out a halo for the origin of the faint surface brightness we see. Rather, this asymmetry is likely an artifact due to a combination of light from a faint ring around this galaxy (as described in Shang et al. 1998, ApJL 504, 23) and residual surface brightness at faint levels from stars that our star-masking procedure cannot completely eliminate. Good agreement with the surface photometry of NGC 5907 by Morrison et al. and other workers lead us to conclude that their data are similarly affected at faint levels by ring light and residual effects from their star masking procedures. Inspection of the images published by Morrison et al. and Sackett et al. confirm this to be the case. Thus, we conclude that NGC 5907 does not have a faint, extended halo.

The edge-on, nearby spiral galaxy NGC5907 has long been used as the prototype of a ``non-interacting'' warped galaxy. We report here the discovery of two interactions with companion dwarf galaxies that substantially change this picture. First, a faint ring structure is discovered around this galaxy that is likely due to the tidal disruption of a companion dwarf spheroidal galaxy. The ring is elliptical in shape with the center of NGC5907 close to one of the ring's foci. This suggests the ring material is in orbit around NGC5907. No gaseous component to the ring has been detected either with deep H$\alpha$ images or in Very Large Array (VLA) HI 21-cm line maps. The visible material in the ring has an integrated luminosity <= 10^8 L_sun and its brightest part has a color R-I ~ 0.9. All of these properties are consistent with the ring being a tidally-disrupted dwarf spheroidal galaxy. Second, we find that NGC5907 has a dwarf companion galaxy, PGC54419, projected to be only 36.9 kpc from the center of NGC5907, close in radial velocity (\Delta V = 45 km s^{-1}) to the giant spiral galaxy. This dwarf is seen at the tip of the HI warp and in the direction of the warp. Hence, NGC5907 can no longer be considered ``non-interacting,'' but is obviously interacting with its dwarf companions much as the Milky Way interacts with its dwarf galaxies. These results, coupled with the finding by others that dwarf galaxies tend to be found around giant galaxies, suggest that tidal interaction with companions, even if containing a mere 1% of the mass of the parent galaxy, might be sufficient to excite the warps found in the disks of many large spiral galaxies.

We present a new, but simple, procedure for estimating the total magnitudes of galaxies. This procedure involves the out-focusing of digital galaxy images numerically, the fitting of the resulting surface-brightness profiles with a single generalised profile model and the extrapolation of the fitted profiles to infinite radial distances. This new system, which we denote `t', differs fundamentally from the `T' system (of the `Reference Catalog of Bright Galaxies' series) in that: (1) it enables a galaxy's luminosity profile to be extrapolated without the need for any prior morphological classification, and (2) it is applicable to images of widely different spatial resolutions (including unresolved ones) because it takes into account systematic effects due to differential image resolution. It also differs fundamentally from the Kron system in that: (1) it can be derived directly from surface photometry without the need to go back to the plate scans or CCD frames (unless the surface photometry is of high resolution and/or the galaxies being measured are very bright), and (2) it can cope with merged images (provided they are separable by image-segmentation software). Through worked examples, we demonstrate the stability of `t'-system total magnitudes with respect to morphological type, the seeing conditions at the time of observation, degree of smoothing and limiting isophote. We also compare and contrast the new system with both the `T' system and the Kron system, and investigate the advantages and limitations of each of the three systems.

We report our new results from wide-field surveys of Herbig-Haro (HH) objects in the nearby star forming regions. The surveys covered approximately 56 degree in Perseus, Taurus, Orion, Monoceros, and other regions. Using refined techniques, we discovered in total 68 HH objects, of which 32 were in Perseus, 4 in Taurus, 13 in Orion, 18 in Monoceros, and 1 in S287 regions. The newly discovered HH objects demonstrate a great variety of morphological structures, including 5 jets, 7 arcs, 12 cirri or cirrus groups, 13 patches, and many knots. These objects provide a new comprehensive database for the study of HH objects in the regions of recent star formation.

We report on a global CCD time-series photometric campaign to decode the pulsations of the nucleus of the planetary nebula NGC1501. The star is hot and hydrogen-deficient, similar to the pre-white-dwarf PG 1159 stars. NGC1501 shows pulsational brightness variations of a few percent with periods ranging from 19 to 87 minutes. The variations are very complex, suggesting a pulsation spectrum that requires a long unbroken time series to resolve. Our CCD photometry of the star covers a two-week period in 1991 November, and used a global network of observatories. We obtained nearly continuous coverage over an interval of one week in the middle of the run. We have identified 10 pulsation periods, ranging from 5235 s down to 1154 s. We find strong evidence that the modes are indeed nonradial g-modes. The ratios of the frequencies of the largest-amplitude modes agree with those expected for modes that are trapped by a density discontinuity in the outer layers. We offer a model for the pulsation spectrum that includes a common period spacing of 22.3 s and a rotation period of 1.17 days; the period spacing allows us to assign a seismological mass of 0.55+/-0.03 Msun.

We present nine color CCD intermediate-band spectrophotometry of a two square degree field centered on the old open cluster M67, from 3890$\rm \AA$ to nearly 1$\mu$. These observations are taken as a part of the BATC (Beijing-Arizona-Taipei-Connecticut) Color Survey of the Sky, for both scientific and calibration reasons. With these data we show that the BATC survey can reach its goal of obtaining spectrophotometry to a zero point accuracy of 0.01 mag, and down to V = 21 with 0.3 mag random error. We fit the color-magnitude diagrams (CMDs) with Worthey's theoretical models. The net result is the excellent fit of the 4.0 Gyr, [Fe/H] = $-0.10$ model to our data, including a good fit to the main sequence (MS) turn-off. Our data are consistent with a toy model with 50\% of the stars in M67 being binaries and a random distribution of binary mass-ratios, although other models with different mass-ratio distributions cannot be ruled out. The spatial distribution and mass function (MF) of stars in M67 show marked effects of dynamical evolution and evaporation of stars from the cluster. Blue stragglers and binary stars are the most condensed within the cluster, with degree of condensation depending on mass.We find M67 to have an elongated shape, oriented at an angle of $15^{\circ}$ relative to the galactic plane. Within its tidal radius, the observed MF of M67 between 1.2 $\rm M_\odot$ and $\rm 0.8 M_\odot$ has a Salpeter slope $\rm \eta = -1.93 \pm 0.66$. For stars of mass below 0.8 $\rm M_\odot$, $\rm \eta \sim 0$. It is plausible that the leveling-off of the MF at lower masses is a result of evaporation of lower mass stars in this mass range at a rate of one every $\sim 10^7$ years. If so, it is plausible that the IMF of M67 has the canonical field value of $\rm \eta = -2.0$.