Comparing chemical abundances of a planet and the host star reveals the origin and formation path. Stellar abundance is measured with high-resolution spectroscopy. Planet abundance, on the other hand, is usually inferred from low-resolution data. For directly imaged exoplanets, the data are available from a slew of high-contrast imaging/spectroscopy instruments. Here, we study the chemical abundance of HR 8799 and its planet c. We measure stellar abundance using LBT/PEPSI (R=120,000) and archival HARPS data: stellar [C/H], [O/H], and C/O are 0.11$\pm$0.12, 0.12$\pm$0.14, and 0.54$^{+0.12}_{-0.09}$, all consistent with solar values. We conduct atmospheric retrieval using newly obtained Subaru/CHARIS data together with archival Gemini/GPI and Keck/OSIRIS data. We model the planet spectrum with petitRADTRANS and conduct retrieval using PyMultiNest. Retrieved planetary abundance can vary by $\sim$0.5 dex, from sub-stellar to stellar C and O abundances. The variation depends on whether strong priors are chosen to ensure a reasonable planet mass. Moreover, comparison with previous works also reveals inconsistency in abundance measurements. We discuss potential issues that can cause the inconsistency, e.g., systematics in individual data sets and different assumptions in the physics and chemistry in retrieval. We conclude that no robust retrieval can be obtained unless the issues are fully resolved.

General relativity is a fully conservative theory, but there exist other possible metric theories of gravity. We consider non-conservative ones with a parameterized post-Newtonian (PPN) parameter, $\zeta_2$. A non-zero $\zeta_2$ induces a self-acceleration for the center of mass of an eccentric binary pulsar system, which contributes to the second time derivative of the pulsar spin frequency, $\ddot{\nu}$. In our work, using the method in Will (1992), we provide an improved analysis with four well-timed, carefully-chosen binary pulsars. In addition, we extend Will's method and derive $\zeta_2$'s effect on the third time derivative of the spin frequency, $\dddot{\nu}$. For PSR B1913+16, the constraint from $\dddot{\nu}$ is even tighter than that from $\ddot{\nu}$. We combine multiple pulsars with Bayesian inference, and obtain an upper limit, $\left|\zeta_{2}\right|<1.3\times10^{-5}$ at 95% confidence level, assuming a flat prior in $\log_{10} \left| \zeta_{2}\right|$. It improves the existing bound by a factor of three. Moreover, we propose an analytical timing formalism for $\zeta_2$. Our simulated times of arrival with simplified assumptions show binary pulsars' capability in limiting $\zeta_{2}$, and useful clues are extracted for real data analysis in future. In particular, we discover that for PSRs B1913+16 and J0737$-$3039A, $\dddot{\nu}$ can yield more constraining limits than $\ddot{\nu}$.

Time series data of celestial objects are commonly used to study valuable and unexpected objects such as extrasolar planets and supernova in time domain astronomy. Due to the rapid growth of data volume, traditional manual methods are becoming extremely hard and infeasible for continuously analyzing accumulated observation data. To meet such demands, we designed and implemented a special tool named AstroCatR that can efficiently and flexibly reconstruct time series data from large-scale astronomical catalogues. AstroCatR can load original catalogue data from Flexible Image Transport System (FITS) files or databases, match each item to determine which object it belongs to, and finally produce time series datasets. To support the high-performance parallel processing of large-scale datasets, AstroCatR uses the extract-transform-load (ETL) preprocessing module to create sky zone files and balance the workload. The matching module uses the overlapped indexing method and an in-memory reference table to improve accuracy and performance. The output of AstroCatR can be stored in CSV files or be transformed other into formats as needed. Simultaneously, the module-based software architecture ensures the flexibility and scalability of AstroCatR. We evaluated AstroCatR with actual observation data from The three Antarctic Survey Telescopes (AST3). The experiments demonstrate that AstroCatR can efficiently and flexibly reconstruct all time series data by setting relevant parameters and configuration files. Furthermore, the tool is approximately 3X faster than methods using relational database management systems at matching massive catalogues.

The 0.5\,m Antarctic Survey Telescopes (AST3) were designed for time-domain optical/infrared astronomy. They are located in Dome~A, Antarctica, where they can take advantage of the continuous dark time during winter. Since the site is unattended in winter, everything for the operation, from observing to data reduction, had to be fully automated. Here, we present a brief overview of the AST3 project and some of its unique characteristics due to its location in Antarctica. We summarise the various components of the survey, including the customized hardware and software, that make complete automation possible.

The location of an astronomical observatory is a key factor that affects its scientific productivity. The best astronomical sites are generally those found at high altitudes. Several such sites in western China and the Tibetan plateau are presently under development for astronomy. One of these is Ali, which at over 5000 m is one of the highest astronomical sites in the world. In order to further investigate the astronomical potential of Ali, we have installed a lunar scintillometer, for the primary purpose of profiling atmospheric turbulence. This paper describes the instrument and technique, and reports results from the first year of observations. We find that ground layer (GL) turbulence at Ali is remarkably weak and relatively thin. The median seeing, from turbulence in the range 11- 500 m above ground is 0.34 arcsec, with seeing better than 0.26 arcsec occurring 25 per cent of the time. Under median conditions, half of the GL turbulence lies below a height of 62 m. These initial results, and the high altitude and relatively low temperatures, suggest that Ali could prove to be an outstanding site for ground-based astronomy.

Rich complexes of associated absorption lines (AALs) in quasar spectra provide unique information about gaseous infall, outflows, and feedback processes in quasar environments. We study five quasars at redshifts 3.1 to 4.4 with AAL complexes containing from 7 to 18 CIV 1548, 1551 systems in high-resolution spectra. These complexes span velocity ranges $\lesssim$3600 km/s within $\lesssim$8200 km/s of the quasar redshifts. All are highly ionised with no measurable low-ionisation ions like SiII or CII, and all appear to form in the quasar/host galaxy environments based on evidence for line locking, partial covering of the background light source, strong NV absorption, and/or roughly solar metallicities, and on the implausibility of such complexes forming in unrelated intervening galaxies. Most of the lines in all five complexes identify high-speed quasar-driven outflows at velocity shifts $v\lesssim -1000$ km/s. Four of the complexes also have lines at smaller blueshifted velocities that might form in ambient interstellar clouds, low-speed outflows or at feedback interfaces in the host galaxies where high-speed winds impact and shred interstellar clouds. The partial covering we measure in some of the high-speed outflow lines require small absorbing clouds with characteristic sizes $\lesssim$1 pc or $\lesssim$0.01 pc. The short survival times of these clouds require locations very close to the quasars, or cloud creation in situ at larger distances perhaps via feedback/cloud-shredding processes. The AAL complex in one quasar, J1008+3623, includes unusually narrow CIV systems at redshifted velocities $350\lesssim v\lesssim640$ km/s that are excellent candidates for gaseous infall towards the quasar, e.g., ''cold-mode" accretion or a gravitationally-bound galactic fountain.

Gamma-ray bursts (GRBs), as extremely energetic explosions in the universe, are widely believed to consist of two stages: the prompt phase and the subsequent afterglow. Recent studies indicate that some high-energy photons are emitted earlier at source than the prompt phase. Due to the light speed variation, these high-energy photons travel slowly than the low-energy photons, so that they are observed after the prompt low-energy photons at the detector. Based on the data from the Fermi Gamma-ray Space Telescope (FGST), we analyse the photon distribution before the prompt emission in detail and propose the existence of a hitherto unknown pre-burst stage of GRBs by adopting a classification method of machine learning. Analysis on the photons automatically selected by machine learning also produce a light speed variation at $E_{\mathrm{LV}}=\mathrm{3.55\times 10^{17} GeV}$.

The strength and vertical distribution of atmospheric turbulence is a key factor determining the performance of optical and infrared telescopes, with and without adaptive optics. Yet, this remains challenging to measure. We describe a new technique using a sequence of short-exposure images of a star field, obtained with a small telescope. Differential motion between all pairs of star images is used to compute the structure functions of longitudinal and transverse wavefront tilt for a range of angular separations. These are compared with theoretical predictions of simple turbulence models by means of a Markov-Chain Monte-Carlo optimization. The method is able to estimate the turbulence profile in the lower atmosphere, the total and free-atmosphere seeing, and the outer scale. We present results of Monte-Carlo simulations used to verify the technique, and show some examples using data from the second AST3 telescope at Dome A in Antarctica.

We investigate the scalar-tensor gravity of Damour and Esposito-Far\`ese (DEF), which predicts non-trivial phenomena in the nonperturbative strong-field regime for neutron stars (NSs). Instead of solving the modified Tolman-Oppenheimer-Volkoff equations, we construct reduced-order surrogate models, coded in the pySTGROM package, to predict the relations of a NS radius, mass, and effective scalar coupling to its central density. Our models are accurate at $\sim1\%$ level and speed up large-scale calculations by two orders of magnitude. As an application, we use pySTGROM and Markov-chain Monte Carlo techniques to constrain parameters in the DEF theory, with five well-timed binary pulsars, the binary NS (BNS) inspiral GW170817, and a hypothetical BNS inspiral in the Advanced LIGO and next-generation GW detectors. In the future, as more binary pulsars and BNS mergers are detected, our surrogate models will be helpful in constraining strong-field gravity with essential speed and accuracy.

A recent proposal to associate 60~TeV to 2~PeV IceCube neutrino events with gamma-ray bursts~(GRBs) indicates the Lorentz violation of cosmic neutrinos and leads further to the $CPT$ symmetry violation between neutrinos and antineutrinos. Here we find that another 12 northern hemisphere track events possibly correlated with GRBs from three-year IceCube data satisfy the same regularity at a lower energy scale around 1~TeV. The combined fitting indicates a Lorentz violation scale ${E}_{\rm LV}=(6.4\pm 1.5)\times10^{17}~{ \rm GeV}$ and an intrinsic time difference ${\Delta {t}_{\rm in}=(-2.8\pm 0.7)\times10^2~{\rm s}}$, from which we find an earlier emission of neutrinos than photons at the GRB source. We also suggest analyzing neutrino events detected a few minutes before the GRB trigger time to test the $CPT$ violation of ultrahigh-energy neutrinos.

We report a ground X-ray calibration of two X-ray telescope prototypes at the PANTER X-ray Test Facility, of the Max-Planck-Institute for Extraterrestrial Physics, in Neuried, Germany. The X-ray telescope prototypes were developed by the Institute of Precision Optical Engineering (IPOE) of Tongji University, in a conical Wolter-I configuration, using thermal glass slumping technology. Prototype #1 with 3 layers and Prototype #2 with 21 layers were tested to assess the prototypes' on-axis imaging performance. The measurement of Prototype #1 indicates a Half Power Diameter (HPD) of 82" at 1.49 keV. As for Prototype #2, we performed more comprehensive measurements of on-axis angular resolution and effective area at several energies ranging from 0.5-10 keV. The HPD and effective area are 111" and 39 cm^2 at 1.49 keV, respectively, at which energy the on-axis performance of the prototypes is our greatest concern.

In Einstein's general relativity, gravity is mediated by a massless spin-2 metric field, and its extension to include a mass for the graviton has profound implication for gravitation and cosmology. In 2002, Finn and Sutton used the gravitational-wave (GW) back-reaction in binary pulsars, and provided the first bound on the mass of graviton. Here we provide an improved analysis using 9 well-timed binary pulsars with a phenomenological treatment. First, individual mass bounds from each pulsar are obtained in the frequentist approach with the help of an ordering principle. The best upper limit on the graviton mass, $m_{g}<3.5\times10^{-20} \, {\rm eV}/c^{2}$ (90% C.L.), comes from the Hulse-Taylor pulsar PSR B1913+16. Then, we combine individual pulsars using the Bayesian theorem, and get $m_{g}<5.2\times10^{-21} \, {\rm eV}/c^{2}$ (90% C.L.) with a uniform prior for $\ln m_g$. This limit improves the Finn-Sutton limit by a factor of more than 10. Though it is not as tight as those from GWs and the Solar System, it provides an independent and complementary bound from a dynamic regime.

TOU is an extremely high resolution optical spectrograph (R=$100,000$, 380-900~nm), which is designed to detect low mass exoplanets using the radial velocity technique.We describe an IDL-based radial velocity (RV) data reduction pipeline for the TOU spectrograph and its performance with stable stars. This pipeline uses a least-squares fitting algorithm to match observed stellar spectra to a high signal-to-noise template created for each star. By carefully controlling all of the error contributions to RV measurements in both the hardware and data pipeline, we have achieved $\sim$0.9\ms long-term RV precision with one of the most RV stable stars, Tau Ceti, similar to what has been achieved with HARPS. This paper presents steps and details in our data pipeline on how to reach the \sms RV precision and also all major error sources which contribute to the final RV measurement uncertainties. The lessons learned in this pipeline development can be applied to other environmentally controlled, very high resolution optical spectrographs to improve RV precision.

We report first results from the CHinese Exoplanet Searching Program from Antarctica (CHESPA)---a wide-field high-resolution photometric survey for transiting exoplanets carried out using telescopes of the AST3 (Antarctic Survey Telescopes times 3) project. There are now three telescopes (AST3-I, AST3-II, and CSTAR-II) operating at Dome A---the highest point on the Antarctic Plateau---in a fully automatic and remote mode to exploit the superb observing conditions of the site, and its long and uninterrupted polar nights. The search for transiting exoplanets is one of the key projects for AST3. During the Austral winters of 2016 and 2017 we used the AST3-II telescope to survey a set of target fields near the southern ecliptic pole, falling within the continuous viewing zone of the TESS mission \citep{Ricker10}. The first data release of the 2016 data, including images, catalogs and lightcurves of 26578 bright stars ($7.5\le i \le15$) was presented in \citet{Zhang18}. The best precision, as measured by the RMS of the lightcurves at the optimum magnitude of the survey ($i=10$), is around 2\,mmag. We detect 222 objects with plausible transit signals from these data, 116 of which are plausible transiting exoplanet candidates according to their stellar properties as given by the TESS Input Catalog \citep{Stassun17}, Gaia DR2 \citep{Gaia18} and TESS-HERMES spectroscopy \citep{Sharma18}. With the first data release from TESS expected in late 2018, this candidate list will be a timely for improving the rejection of potential false positives.

Located at Dome A, the highest point of the Antarctic plateau, the Chinese Kunlun station is considered to be one of the best ground-based photometric sites because of its extremely cold, dry, and stable atmosphere(Saunders et al. 2009). A target can be monitored from there for over 40 days without diurnal interruption during a polar winter. This makes Kunlun station a perfect site to search for short-period transiting exoplanets. Since 2008, an observatory has been built at Kunlun station and three telescopes are working there. Using these telescopes, the AST3 project has been carried out over the last six years with a search for transiting exoplanets as one of its key programs (CHESPA). In the austral winters of 2016 and 2017, a set of target fields in the Southern CVZ of TESS (Ricker et al. 2009) were monitored by the AST3-II telescope. In this paper, we introduce the CHESPA and present the first data release containing photometry of 26,578 bright stars (m_i < 15). The best photometric precision at the optimum magnitude for the survey is around 2 mmag. To demonstrate the data quality, we also present a catalog of 221 variables with a brightness variation greater than 5 mmag from the 2016 data. Among these variables, 179 are newly identified periodic variables not listed in the AAVSO databasea), and 67 are listed in the Candidate Target List(Stassun et al. 2017). These variables will require careful attention to avoid false-positive signals when searching for transiting exoplanets. Dozens of new transiting exoplanet candidates will be also released in a subsequent paper(Zhang et al. 2018b).

We examine a method to detect the light speed variation from gamma ray burst data observed by the Fermi Gamma-ray Space Telescope (FGST). We suggest new criteria to determine the characteristic time for low energy photons by the energy curve and the average energy curve, and obtain similar results compared with those from the light curve. We offer a new criterion with both the light curve and the average energy curve to determine the characteristic time for low energy photons. We then apply the new criteria to the GBM NaI data, the GBM BGO data, and the LAT LLE data, and obtain consistent results for three different sets of low energy photons from different FERMI detectors.

We present an analysis of meteorological data from the second generation of the Kunlun Automated Weather Station (KLAWS-2G) at Dome A, Antarctica during 2015 and 2016. We find that a strong temperature inversion exists for all the elevations up to 14 m that KLAWS-2G can reach, and lasts for more than 10 hours for 50% or more of the time when temperature inversion occurs. The average wind speeds at 4 m elevation are 4.2 m/s and 3.8 m/s during 2015 and 2016, respectively. The strong temperature inversion and moderate wind speed lead to a shallow turbulent boundary layer height at Dome A. By analyzing the temperature and wind shear profiles, we note telescopes should be elevated by at least 8 m above the ice. We also find that the duration of temperature inversions, and the wind speed, vary considerably from year to year. Therefore, long-term and continuous data are still needed for the site survey at Dome A.

The Lorentz violation~(LV) effect of ultra-relativistic particles can be tested by gamma-ray burst~(GRB) neutrinos and photons. The IceCube Collaboration has observed plenty of ultra-high energy neutrinos, including four events of PeV scale neutrinos. Recent studies suggested a possible energy dependent speed variation of GRB neutrinos in a similar way to that of GRB photons. Here we find that all four events of PeV neutrinos with associated GRB candidates can satisfy a regularity found from TeV neutrinos about a linear form correlation between the observed time difference and the LV factor. Such regularity indicates a Lorentz violation scale $E_{\rm LV}=(6.5\pm 0.4)\times10^{17}~{\rm GeV}$, which is comparable with that determined by GRB photons. We also suggest that neutrinos and anti-neutrinos can be superluminal and subluminal respectively due to opposite signs of LV correction.

The Dharma Planet Survey (DPS) aims to monitor about 150 nearby very bright FGKM dwarfs (within 50 pc) during 2016$-$2020 for low-mass planet detection and characterization using the TOU very high resolution optical spectrograph (R$\approx$100,000, 380-900nm). TOU was initially mounted to the 2-m Automatic Spectroscopic Telescope at Fairborn Observatory in 2013-2015 to conduct a pilot survey, then moved to the dedicated 50-inch automatic telescope on Mt. Lemmon in 2016 to launch the survey. Here we report the first planet detection from DPS, a super-Earth candidate orbiting a bright K dwarf star, HD 26965. It is the second brightest star ($V=4.4$ mag) on the sky with a super-Earth candidate. The planet candidate has a mass of 8.47$\pm0.47M_{\rm Earth}$, period of $42.38\pm0.01$ d, and eccentricity of $0.04^{+0.05}_{-0.03}$. This RV signal was independently detected by Diaz et al. (2018), but they could not confirm if the signal is from a planet or from stellar activity. The orbital period of the planet is close to the rotation period of the star (39$-$44.5 d) measured from stellar activity indicators. Our high precision photometric campaign and line bisector analysis of this star do not find any significant variations at the orbital period. Stellar RV jitters modeled from star spots and convection inhibition are also not strong enough to explain the RV signal detected. After further comparing RV data from the star's active magnetic phase and quiet magnetic phase, we conclude that the RV signal is due to planetary-reflex motion and not stellar activity.

The three Antarctic Survey Telescopes (AST3) aim to carry out time domain imaging survey at Dome A, Antarctica. The first of the three telescopes (AST3-1) was successfully deployed on January 2012. AST3-1 is a 500\,mm aperture modified Schmidt telescope with a 680\,mm diameter primary mirror. AST3-1 is equipped with a SDSS $i$ filter and a 10k $\times$ 10k frame transfer CCD camera, reduced to 5k $\times$ 10k by electronic shuttering, resulting in a 4.3 deg$^2$ field-of-view. To verify the capability of AST3-1 for a variety of science goals, extensive commissioning was carried out between March and May 2012. The commissioning included a survey covering 2000 deg$^2$ as well as the entire Large and Small Magellanic Clouds. Frequent repeated images were made of the center of the Large Magellanic Cloud, a selected exoplanet transit field, and fields including some Wolf-Rayet stars. Here we present the data reduction and photometric measurements of the point sources observed by AST3-1. We have achieved a survey depth of 19.3\,mag in 60 s exposures with 5\,mmag precision in the light curves of bright stars. The facility achieves sub-mmag photometric precision under stable survey conditions, approaching its photon noise limit. These results demonstrate that AST3-1 at Dome A is extraordinarily competitive in time-domain astronomy, including both quick searches for faint transients and the detection of tiny transit signals.

Combining stellar atmospheric parameters, such as effective temperature, surface gravity, and metallicity, with barycentric radial velocity data provides insight into the chemo-dynamics of the MilkyWay and our local Galactic environment. We analyze 3075 stars with spectroscopic data from the Sloan Digital Sky Survey III (SDSS-III) MARVELS radial velocity survey and present atmospheric parameters for 2343 dwarf stars using the spectral indices method, a modified version of the equivalent width method. We present barycentric radial velocities for 2610 stars with uncertainties of $\sim$0.3 km s$^{-1}$. We determine stellar ages using two independent methods and calculate ages for 2343 stars with a maximum-likelihood isochronal age-dating method and for 1694 stars with a Bayesian age-dating method. Using previously published parallax data we compute Galactic orbits and space velocities for 1995 stars to explore stellar populations based on kinematic and age parameters. This study combines good ages and exquisite velocities to explore local chemo-kinematics of the Milky Way, which complements many of the recent studies of giant stars with the APOGEE survey, and we find our results to be in agreement with current chemo-dynamical models of the Milky Way. Particularly, we find from our metallicity distributions and velocity-age relations of a kinematically-defined thin disk that the metal rich end has stars of all ages, even after we clean the sample of highly eccentric stars, suggesting that radial migration plays a key role in the metallicity scatter of the thin disk. All stellar parameters and kinematic data derived in this work are catalogued and published online in machine-readable form.

The effect of Quantum Gravity (QG) may bring a tiny light speed variation as $v(E)=c(1-E/E_{\rm LV})$, where $E$ is the photon energy and $E_{\rm LV}$ is a Lorentz violation scale. A remarkable regularity was suggested in previous studies to look for the light speed variation from high energy photon events of Gamma Ray Bursts (GRBs). We provide a general analysis on the data of 25 bright GRBs observed by the Fermi Gamma-ray Space Telescope (FGST). Such method allows a completed scan over all possibilities in a more clean and impartial way without any bias compared to previous intuitive analysis. The results show that with the increase in the intrinsic energies of photons, such regularity truly emerges and gradually becomes significant. For photons with intrinsic energies higher than 40~GeV, the regularity exists at a significance of 3-5~$\sigma$ with $E_{\rm LV}=3.6\times 10^{17}~\rm GeV$ determined by the GRB data.

The LIGO detection of gravitational waves (GW) from merging black holes in 2015 marked the beginning of a new era in observational astronomy. The detection of an electromagnetic signal from a GW source is the critical next step to explore in detail the physics involved. The Antarctic Survey Telescopes (AST3), located at Dome A, Antarctica, is uniquely situated for rapid response time-domain astronomy with its continuous night-time coverage during the austral winter. We report optical observations of the GW source (GW~170817) in the nearby galaxy NGC 4993 using AST3. The data show a rapidly fading transient at around 1 day after the GW trigger, with the $i$-band magnitude declining from $17.23\pm0.13$ magnitude to $17.72\pm0.09$ magnitude in $\sim 1.8$ hour. The brightness and time evolution of the optical transient associated with GW~170817 are broadly consistent with the predictions of models involving merging binary neutron stars. We infer from our data that the merging process ejected about $\sim 10^{-2}$ solar mass of radioactive material at a speed of up to $30\%$ the speed of light.

The discovery of the first electromagnetic counterpart to a gravitational wave signal has generated follow-up observations by over 50 facilities world-wide, ushering in the new era of multi-messenger astronomy. In this paper, we present follow-up observations of the gravitational wave event GW170817 and its electromagnetic counterpart SSS17a/DLT17ck (IAU label AT2017gfo) by 14 Australian telescopes and partner observatories as part of Australian-based and Australian-led research programs. We report early- to late-time multi-wavelength observations, including optical imaging and spectroscopy, mid-infrared imaging, radio imaging, and searches for fast radio bursts. Our optical spectra reveal that the transient source afterglow cooled from approximately 6400K to 2100K over a 7-day period and produced no significant optical emission lines. The spectral profiles, cooling rate, and photometric light curves are consistent with the expected outburst and subsequent processes of a binary neutron star merger. Star formation in the host galaxy probably ceased at least a Gyr ago, although there is evidence for a galaxy merger. Binary pulsars with short (100 Myr) decay times are therefore unlikely progenitors, but pulsars like PSR B1534+12 with its 2.7 Gyr coalescence time could produce such a merger. The displacement (about 2.2 kpc) of the binary star system from the centre of the main galaxy is not unusual for stars in the host galaxy or stars originating in the merging galaxy, and therefore any constraints on the kick velocity imparted to the progenitor are poor.

Planet searches using the radial velocity technique show a paucity of companions to solar-type stars within ~5 AU in the mass range of ~10 - 80 M$_{\text{Jup}}$. This deficit, known as the brown dwarf desert, currently has no conclusive explanation. New substellar companions in this region help asses the reality of the desert and provide insight to the formation and evolution of these objects. Here we present 10 new brown dwarf and two low-mass stellar companion candidates around solar-type stars from the Multi-object APO Radial-Velocity Exoplanet Large-Area Survey (MARVELS) of the Sloan Digital Sky Survey III (SDSS-III). These companions were selected from processed MARVELS data using the latest University of Florida Two Dimensional (UF2D) pipeline, which shows significant improvement and reduction of systematic errors over previous pipelines. The 10 brown dwarf companions range in mass from ~13 to 76 M$_{\text{Jup}}$ and have orbital radii of less than 1 AU. The two stellar companions have minimum masses of ~98 and 100 M$_{\text{Jup}}$. The host stars of the MARVELS brown dwarf sample have a mean metallicity of [Fe/H] = 0.03 $\pm$ 0.08 dex. Given our stellar sample we estimate the brown dwarf occurrence rate around solar-type stars with periods less than ~300 days to be ~0.56%.

AST3-1 is the second-generation wide-field optical photometric telescope dedicated to time domain astronomy at Dome A, Antarctica. Here we present the results of $i$ band images survey from AST3-1 towards one Galactic disk field. Based on time-series photometry of 92,583 stars, 560 variable stars were detected with $i$ magnitude $\leq$ 16.5 mag during eight days of observations; 339 of these are previously unknown variables. We tentatively classify the 560 variables as 285 eclipsing binaries (EW, EB, EA), 27 pulsating variable stars ($\delta$~Scuti, $\gamma$~Doradus, $\delta$~Cephei variable and RR Lyrae stars) and 248 other types of variables (unclassified periodic, multi-periodic and aperiodic variable stars). Among the eclipsing binaries, 34 show O'Connell effects. One of the aperiodic variables shows a plateau light curve and another one shows a secondary maximum after peak brightness. We also detected a complex binary system with RS CVn-like light curve morphology; this object is being followed-up spectroscopically using the Gemini South telescope.

Our Universe is described jointly by Einstein's relativity principle and the quantum principle; there the existence of the smallest units of matter, such as electrons, neutrinos, quarks, and photons, is well established and the smallest units of matter are described by the Standard Model. Based on these, we have a clear picture for the propagation of ultra high energy cosmic rays (UHECR's), for energies greater than $E > 10^{13}\, eV$ but less than $E < 10^{26}\, eV$, in the cosmic medium of the Cosmic Microwave Background (CMB) and of the cosmic background neutrinos (CB$\nu$'s). We find that the CMB plays a pivot role in this energy range, so that the observed "knee(s)" and the "ankle" could be understood in reasonable terms. What we might observe at the energy near $10^{25}\,eV$', the so-called "$W^\pm$ bursts" or "$Z^0$ bursts", is briefly discussed. Meanwhile, in the early Universe, the CB$\nu$'s, in view of the tiny neutrino mass, gave rise to the clustering phenomenon, resulting in a neutrino halo attached to each visual ordinary-matter object, such as the Earth, the Sun, and stars - with the $PeV$ neutrinos near the surface of the Fermi-Dirac sphere of the neutrino halo of the Earth playing the most fundamental role.

We study the probability distribution function (PDF) of mass surface densities, $\Sigma$, of infrared dark cloud (IRDC) G028.37+00.07 and its surrounding giant molecular cloud. This PDF constrains the physical processes, such as turbulence, magnetic fields and self-gravity, that are expected to be controlling cloud structure and star formation activity. The chosen IRDC is of particular interest since it has almost 100,000 solar masses within a radius of 8 parsecs, making it one of the most massive, dense molecular structures known and is thus a potential site for the formation of a "super star cluster." We study $\Sigma$ in two ways. First, we use a combination of NIR and MIR extinction maps that are able to probe the bulk of the cloud structure up to $\Sigma\sim1\:{\rm{g\:cm}^{-2}}\:$($A_V\simeq200$ mag). Second, we study the FIR and sub-mm dust continuum emission from the cloud utilizing Herschel PACS and SPIRE images and paying careful attention to the effects of foreground and background contamination. We find that the PDFs from both methods, applied over a $\sim20^\prime$(30 pc)-sized region that contains $\simeq1.5\times10^5\:M_\odot$ and encloses a minimum closed contour with $\Sigma\simeq0.013\:{\rm{g\:cm}^{-2}}\:$($A_V\simeq3$ mag), shows a log-normal shape with the peak measured at $\Sigma\simeq0.021\:{\rm{g\:cm}^{-2}}\:$($A_V\simeq4.7$ mag). There is tentative evidence for the presence of a high-$\Sigma$ power law tail that contains from $\sim3\%$ to 8\% of the mass of the cloud material. We discuss the implications of these results for the physical processes occurring in this cloud.

The effect of quantum gravity can bring a tiny light speed variation which is detectable through energetic photons propagating from gamma ray bursts (GRBs) to an observer such as the space observatory. Through an analysis of the energetic photon data of the GRBs observed by the Fermi Gamma-ray Space Telescope (FGST), we reveal a surprising regularity of the observed time lags between photons of different energies with respect to the Lorentz violation factor due to the light speed energy dependence. Such regularity suggests a linear form correction of the light speed $v(E)=c(1-E/E_{\rm LV})$, where $E$ is the photon energy and $E_{\rm LV}=(3.60 \pm 0.26) \times 10^{17}~ \rm GeV$ is the Lorentz violation scale measured by the energetic photon data of GRBs. The results support an energy dependence of the light speed in cosmological space.

We report the detections of a giant planet (MARVELS-7b) and a brown dwarf candidate (MARVELS-7c) around the primary star in the close binary system, HD 87646. It is the first close binary system with more than one substellar circum-primary companion discovered to the best of our knowledge. The detection of this giant planet was accomplished using the first multi-object Doppler instrument (KeckET) at the Sloan Digital Sky Survey (SDSS) telescope. Subsequent radial velocity observations using ET at Kitt Peak National Observatory, HRS at HET, the "Classic" spectrograph at the Automatic Spectroscopic Telescope at Fairborn Observatory, and MARVELS from SDSS-III confirmed this giant planet discovery and revealed the existence of a long-period brown dwarf in this binary. HD 87646 is a close binary with a separation of $\sim22$ AU between the two stars, estimated using the Hipparcos catalogue and our newly acquired AO image from PALAO on the 200-inch Hale Telescope at Palomar. The primary star in the binary, HD 87646A, has Teff = 5770$\pm$80K, log(g)=4.1$\pm$0.1 and [Fe/H] = $-0.17\pm0.08$. The derived minimum masses of the two substellar companions of HD 87646A are 12.4$\pm$0.7M$_{\rm Jup}$ and 57.0$\pm3.7$M$_{\rm Jup}$. The periods are 13.481$\pm$0.001 days and 674$\pm$4 days and the measured eccentricities are 0.05$\pm$0.02 and 0.50$\pm$0.02 respectively. Our dynamical simulations show the system is stable if the binary orbit has a large semi-major axis and a low eccentricity, which can be verified with future astrometry observations.

eXTP is a science mission designed to study the state of matter under extreme conditions of density, gravity and magnetism. Primary targets include isolated and binary neutron stars, strong magnetic field systems like magnetars, and stellar-mass and supermassive black holes. The mission carries a unique and unprecedented suite of state-of-the-art scientific instruments enabling for the first time ever the simultaneous spectral-timing-polarimetry studies of cosmic sources in the energy range from 0.5-30 keV (and beyond). Key elements of the payload are: the Spectroscopic Focusing Array (SFA) - a set of 11 X-ray optics for a total effective area of about 0.9 m^2 and 0.6 m^2 at 2 keV and 6 keV respectively, equipped with Silicon Drift Detectors offering <180 eV spectral resolution; the Large Area Detector (LAD) - a deployable set of 640 Silicon Drift Detectors, for a total effective area of about 3.4 m^2, between 6 and 10 keV, and spectral resolution <250 eV; the Polarimetry Focusing Array (PFA) - a set of 2 X-ray telescope, for a total effective area of 250 cm^2 at 2 keV, equipped with imaging gas pixel photoelectric polarimeters; the Wide Field Monitor (WFM) - a set of 3 coded mask wide field units, equipped with position-sensitive Silicon Drift Detectors, each covering a 90 degrees x 90 degrees FoV. The eXTP international consortium includes mostly major institutions of the Chinese Academy of Sciences and Universities in China, as well as major institutions in several European countries and the United States. The predecessor of eXTP, the XTP mission concept, has been selected and funded as one of the so-called background missions in the Strategic Priority Space Science Program of the Chinese Academy of Sciences since 2011. The strong European participation has significantly enhanced the scientific capabilities of eXTP. The planned launch date of the mission is earlier than 2025.

It is postulated in Einstein's relativity that the speed of light in vacuum is a constant for all observers. However, the effect of quantum gravity could bring an energy dependence of light speed. Even a tiny speed variation, when amplified by the cosmological distance, may be revealed by the observed time lags between photons with different energies from astrophysical sources. From the newly detected long gamma ray burst GRB~160509A, we find evidence to support the prediction for a linear form modification of light speed in cosmological space.

The census of exoplanets is incomplete for orbital distances larger than 1 AU. Here, we present 41 long-period planet candidates in 38 systems identified by Planet Hunters based on Kepler archival data (Q0-Q17). Among them, 17 exhibit only one transit, 14 have two visible transits and 10 have more than three visible transits. For planet candidates with only one visible transit, we estimate their orbital periods based on transit duration and host star properties. The majority of the planet candidates in this work (75%) have orbital periods that correspond to distances of 1-3 AU from their host stars. We conduct follow-up imaging and spectroscopic observations to validate and characterize planet host stars. In total, we obtain adaptive optics images for 33 stars to search for possible blending sources. Six stars have stellar companions within 4". We obtain high-resolution spectra for 6 stars to determine their physical properties. Stellar properties for other stars are obtained from the NASA Exoplanet Archive and the Kepler Stellar Catalog by Huber et al. (2014). We validate 7 planet candidates that have planet confidence over 0.997 (3-{\sigma} level). These validated planets include 3 single-transit planets (KIC-3558849b, KIC-5951458b, and KIC-8540376c), 3 planets with double transits (KIC-8540376b, KIC-9663113b, and KIC-10525077b), and 1 planet with 4 transits (KIC-5437945b). This work provides assessment regarding the existence of planets at wide separations and the associated false positive rate for transiting observation (17%-33%). More than half of the long-period planets with at least three transits in this paper exhibit transit timing variations up to 41 hours, which suggest additional components that dynamically interact with the transiting planet candidates. The nature of these components can be determined by follow-up radial velocity and transit observations.

We report the discovery of KELT-4Ab, an inflated, transiting Hot Jupiter orbiting the brightest component of a hierarchical triple stellar system. The host star is an F star with $T_{\rm eff}=6206\pm75$ K, $\log g=4.108\pm0.014$, $\left[{\rm Fe}/{\rm H}\right]=-0.116_{-0.069}^{+0.065}$, ${\rm M_*}=1.201_{-0.061}^{+0.067} \ {\rm M}_{\odot}$, and ${\rm R_*}=1.610_{-0.068}^{+0.078} \ {\rm R}_{\odot}$. The best-fit linear ephemeris is $\rm {BJD_{TDB}} = 2456193.29157 \pm 0.00021 + E\left(2.9895936 \pm 0.0000048\right)$. With a magnitude of $V\sim10$, a planetary radius of $1.699_{-0.045}^{+0.046} \ {\rm R_J}$, and a mass of $0.902_{-0.059}^{+0.060} \ {\rm M_J}$, it is the brightest host among the population of inflated Hot Jupiters ($R_P > 1.5R_J$), making it a valuable discovery for probing the nature of inflated planets. In addition, its existence within a hierarchical triple and its proximity to Earth ($210$ pc) provides a unique opportunity for dynamical studies with continued monitoring with high resolution imaging and precision radial velocities. In particular, the motion of the binary stars around each other and of both stars around the primary star relative to the measured epoch in this work should be detectable when it rises in October 2015.

We describe the search for Lyman-break galaxies (LBGs) near the sub-millimeter bright starburst galaxy HFLS3 at $z$$=$6.34 and a study on the environment of this massive galaxy during the end of reionization.We performed two independent selections of LBGs on images obtained with the \textit{Gran Telescopio Canarias} (GTC) and the \textit{Hubble Space Telescope} (HST) by combining non-detections in bands blueward of the Lyman-break and color selection. A total of 10 objects fulfilling the LBG selection criteria at $z$$>$5.5 were selected over the 4.54 and 55.5 arcmin$^2$ covered by our HST and GTC images, respectively. The photometric redshift, UV luminosity, and the star-formation rate of these sources were estimated with models of their spectral energy distribution. These $z$$\sim$6 candidates have physical properties and number densities in agreement with previous results. The UV luminosity function at $z$$\sim$6 and a Voronoi tessellation analysis of this field shows no strong evidence for an overdensity of relatively bright objects (m$_{F105W}$$<$25.9) associated with \textit{HFLS3}. However, the over-density parameter deduced from this field and the surface density of objects can not excluded definitively the LBG over-density hypothesis. Moreover we identified three faint objects at less than three arcseconds from \textit{HFLS3} with color consistent with those expected for $z$$\sim$6 galaxies. Deeper data are needed to confirm their redshifts and to study their association with \textit{HFLS3} and the galaxy merger that may be responsible for the massive starburst.

We present the rest-frame optical spectral energy distribution and stellar masses of six Herschel- selected gravitationally lensed dusty, star-forming galaxies (DSFGs) at 1 < z < 3. These galaxies were first identified with Herschel/SPIRE imaging data from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) and the Herschel Multi-tiered Extragalactic Survey (HerMES). The targets were observed with Spitzer/IRAC at 3.6 and 4.5um. Due to the spatial resolution of the IRAC observations at the level of 2 arcseconds, the lensing features of a background DSFG in the near-infrared are blended with the flux from the foreground lensing galaxy in the IRAC imaging data. We make use of higher resolution Hubble/WFC3 or Keck/NIRC2 Adaptive Optics imaging data to fit light profiles of the foreground lensing galaxy (or galaxies) as a way to model the foreground components, in order to successfully disentangle the foreground lens and background source flux densities in the IRAC images. The flux density measurements at 3.6 and 4.5um, once combined with Hubble/WFC3 and Keck/NIRC2 data, provide important constraints on the rest-frame optical spectral energy distribution of the Herschel-selected lensed DSFGs. We model the combined UV- to millimeter-wavelength SEDs to establish the stellar mass, dust mass, star-formation rate, visual extinction, and other parameters for each of these Herschel-selected DSFGs. These systems have inferred stellar masses in the range 8 x 10^10 to 4 x 10^11 Msun and star-formation rates of around 100 Msun yr-1. This puts these lensed sub-millimeter systems well above the SFR-M* relation observed for normal star-forming galaxies at similar redshifts. The high values of SFR inferred for these systems are consistent with a major merger-driven scenario for star formation.

We present Hubble Space Telescope (HST) WFC3 imaging and grism spectroscopy observations of the {\it Herschel}-selected gravitationally-lensed starburst galaxy HATLASJ1429-0028. The lensing system consists of an edge-on foreground disk galaxy at $z=0.218$ with a nearly complete Einstein ring of the infrared luminous galaxy at $z=1.027$. The WFC3 spectroscopy with G102 and G141 grisms, covering the wavelength range of 0.8 to 1.7 $\mu$m, resulted in detections of H$\alpha$+[NII], H$\beta$, [SII], and [OIII] for the background galaxy from which we measure line fluxes and ratios. The Balmer line ratio H$\alpha$/H$\beta$ of $7.5 \pm 4.4$, when corrected for [NII], results in an extinction for the starburst galaxy of $E(B-V)=0.8 \pm 0.5$. The $H\alpha$ based star-formation rate, when corrected for extinction, is $60 \pm 50$ M$_{\odot}$ yr$^{-1}$, lower than the instantaneous star-formation rate of 390 $\pm$ 90 M$_{\odot}$ yr$^{-1}$ from the total IR luminosity. We also compare the nebular line ratios of HATLASJ1429-0028 with other star-forming and sub-mm bright galaxies. The nebular line ratios are consistent with an intrinsic ultra-luminous infrared galaxy with no evidence for excitation by an active galactic nucleus (AGN). We estimate the metallicity, $12 + log(O/H)$, of HATLASJ1429-0028 to be 8.49 $\pm$ 0.16. Such a low value is below the average relations for stellar mass vs. metallicity of galaxies at $z \sim 1$ for a galaxy with stellar mass of $\sim 2 \times 10^{11}$ M$_{\odot}$. The combination of high stellar mass, lack of AGN indicators, low metallicity, and the high star-formation rate of HATLASJ1429-0028 suggest that this galaxy is currently undergoing a rapid formation.

The Kepler mission has yielded a large number of planet candidates from among the Kepler Objects of Interest (KOIs), but spectroscopic follow-up of these relatively faint stars is a serious bottleneck in confirming and characterizing these systems. We present motivation and survey design for an ongoing project with the SDSS-III multiplexed APOGEE near-infrared spectrograph to monitor hundreds of KOI host stars. We report some of our first results using representative targets from our sample, which include current planet candidates that we find to be false positives, as well as candidates listed as false positives that we do not find to be spectroscopic binaries. With this survey, KOI hosts are observed over ~20 epochs at a radial velocity precision of 100-200 m/s. These observations can easily identify a majority of false positives caused by physically-associated stellar or substellar binaries, and in many cases, fully characterize their orbits. We demonstrate that APOGEE is capable of achieving RV precision at the 100-200 m/s level over long time baselines, and that APOGEE's multiplexing capability makes it substantially more efficient at identifying false positives due to binaries than other single-object spectrographs working to confirm KOIs as planets. These APOGEE RVs enable ancillary science projects, such as studies of fundamental stellar astrophysics or intrinsically rare substellar companions. The coadded APOGEE spectra can be used to derive stellar properties (T_eff, log(g)) and chemical abundances of over a dozen elements to probe correlations of planet properties with individual elemental abundances.

We present Keck-Adaptive Optics and Hubble Space Telescope high resolution near-infrared (IR) imaging for 500 um-bright candidate lensing systems identified by the Herschel Multi-tiered Extra-galactic Survey (HerMES) and Herschel Astrophysical Terahertz Survey (H-ATLAS). Out of 87 candidates with near-IR imaging, 15 (~17%) display clear near-IR lensing morphologies. We present near-IR lens models to reconstruct and recover basic rest-frame optical morphological properties of the background galaxies from 12 new systems. Sources with the largest near-IR magnification factors also tend to be the most compact, consistent with the size bias predicted from simulations and pre- vious lensing models for sub-millimeter galaxies. For four new sources that also have high-resolution sub-mm maps, we test for differential lensing between the stellar and dust components and find that the 880 um magnification factor (u_880) is ~1.5 times higher than the near-IR magnification factor (u_NIR), on average. We also find that the stellar emission is ~2 times more extended in size than dust. The rest-frame optical properties of our sample of Herschel-selected lensed SMGs are consistent with those of unlensed SMGs, which suggests that the two populations are similar.

Twilight/night sky images are often used for flat-fielding CCD images, but the brightness gradient in twilight/night sky causes problems of accurate flat-field correction in astronomical images for wide-field telescopes. Using data from the Antarctic Survey Telescope (AST3), we found that when the sky brightness gradient is minimum and stable, there is still a gradient of 1% across AST3's field-of-view of 4.3 square degrees. We tested various approaches to remove the varying gradients in individual flat-field images. Our final optimal method can reduce the spatially dependent errors caused by the gradient to the negligible level. We also suggest a guideline of flat-fielding using twilight/night sky images for wide-field robotic autonomous telescopes.

The photon transfer curve (PTC, variance vs. signal level) is a commonly used and effective tool in characterizing CCD performance. It is theoretically linear in the range where photon shot noise dominates, and its slope is utilized to derive the gain of the CCD. However, recent researches on different CCDs have revealed that the variance progressively drops at high signal levels, while the linearity shown by signal versus exposure time is still excellent and unaffected. On the other hand, bright stars are found to exhibit fatter point spread function (PSF). Both nonlinear PTC and the brighter-fatter effect are regarded as the result of spreading of charges between pixels, an interaction progress increasing with signal level. In this work we investigate the nonlinear PTC based on the images with a STA1600FT CCD camera, whose PTC starts to become nonlinear at about 1/3 full well. To explain the phenomenon, we present a model to characterize the charge-sharing PSF. This signal-dependent PSF can be derived from flat-field frames, and allow us to quantify the effects on photometry and measured shape of stars. This effect is essentially critical for projects requiring accurate photometry and shape parameters.

We have developed a new method to correct dark current at relatively high temperatures for Charge-Coupled Device (CCD) images when dark frames cannot be obtained on the telescope. For images taken with the Antarctic Survey Telescopes (AST3) in 2012, due to the low cooling efficiency, the median CCD temperature was -46$^\circ$C, resulting in a high dark current level of about 3$e^-$/pix/sec, even comparable to the sky brightness (10$e^-$/pix/sec). If not corrected, the nonuniformity of the dark current could even overweight the photon noise of the sky background. However, dark frames could not be obtained during the observing season because the camera was operated in frame-transfer mode without a shutter, and the telescope was unattended in winter. Here we present an alternative, but simple and effective method to derive the dark current frame from the scientific images. Then we can scale this dark frame to the temperature at which the scientific images were taken, and apply the dark frame corrections to the scientific images. We have applied this method to the AST3 data, and demonstrated that it can reduce the noise to a level roughly as low as the photon noise of the sky brightness, solving the high noise problem and improving the photometric precision. This method will also be helpful for other projects that suffer from similar issues.

We present an analysis of the meteorological data collected at Dome A, Antarctica by the Kunlun Automated Weather Station, including temperatures and wind speeds at eight elevations above the snow surface between 0m and 14.5m. The average temperatures at 2m and 14.5m are $-54^{\circ}$C and $-46^{\circ}$C, respectively. We find that a strong temperature inversion existed at all heights for more than 70% of the time, and the temperature inversion typically lasts longer than 25 hours, indicating an extremely stable atmosphere. The temperature gradient is larger at lower elevations than higher elevations. The average wind speed was 1.5m/s at 4m elevation. We find that the temperature inversion is stronger when the wind speed is lower and the temperature gradient decreases sharply at a specific wind speed for each elevation. The strong temperature inversion and low wind speed results in a shallow and stable boundary layer with weak atmospheric turbulence above it, suggesting that Dome A should be an excellent site for astronomical observations. All the data from the weather station are available for download.

The constancy of light speed is a basic assumption in Einstein's special relativity, and consequently the Lorentz invariance is a fundamental symmetry of space-time in modern physics. However, it is speculated that the speed of light becomes energy-dependent due to the Lorentz invariance violation~(LV) in various new physics theories. We analyse the data of the energetic photons from the gamma-ray bursts (GRBs) by the Fermi Gamma-Ray Space Telescope, and find more events to support the energy dependence in the light speed with both linear and quadratic form corrections. We provide two scenarios to understand all the new-released Pass~8 data of bright GRBs by the Fermi-LAT Collaboration, with predictions from such scenarios being testable by future detected GRBs.

We discuss the restframe UV emission from the starbursting galaxy HFLS3 at z=6.34, discovered in Herschel/SPIRE data due to its red color in the submm wavelengths from 250-500 um. The apparent inst. SFR of HFLS3 inferred from the total FIR luminosity measured with over 15 photometric data points between 100 to 1000 um is 2900 Msun/yr. Keck/NIRC2 Ks band adaptive optics imaging data showed two potential NIR counterparts near HFLS3. Previously, the northern galaxy was taken to be in the foreground at z=2.1 while the southern galaxy was assumed to HFLS3's NIR counterpart. New HST/WFC3 and ACS imaging data show both optically bright galaxies are in the foreground at z<6. A new lensing model based on HST data and mm-wave continuum emission yields a magnification of 2.2+/-0.3. The lack of multiple imaging constrains the lensing magnification to be lower than either 2.7 or 3.5 at the 95% confidence level for the two scenarios, which attribute one or two components to HFLS3 in the source plane. Correcting for gravitational lensing, the inst. SFR is 1320 Msun/yr with the 95% confidence lower limit around 830 Msun/yr. Using models for the restframe UV to FIR SED, the ave. SFR over the last 100 Myr is around 660 Msun/yr. The dust and stellar masses of HFLS3 from the same SED models are 3x10^8 Msun and ~5x10^10 Msun, respectively, with large systematic uncertainties on assumptions related to the SED model. With HST/WFC3 images we also find diffuse NIR emission about 0.5" (~3 kpc) SW of HFLS3 that remains undetected in the ACS data. The emission has a photometric redshift consistent with either z~6 or a dusty galaxy template at z~2. If at the same redshift as HFLS3 the detected diffuse emission could be part of the complex merger system that could be triggering the starburst. Alternatively, it could be part of the foreground structure at z~2.1 that is responsible for lensing of HFLS3.

The Sloan Digital Sky Survey (SDSS) has been in operation since 2000 April. This paper presents the tenth public data release (DR10) from its current incarnation, SDSS-III. This data release includes the first spectroscopic data from the Apache Point Observatory Galaxy Evolution Experiment (APOGEE), along with spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS) taken through 2012 July. The APOGEE instrument is a near-infrared R~22,500 300-fiber spectrograph covering 1.514--1.696 microns. The APOGEE survey is studying the chemical abundances and radial velocities of roughly 100,000 red giant star candidates in the bulge, bar, disk, and halo of the Milky Way. DR10 includes 178,397 spectra of 57,454 stars, each typically observed three or more times, from APOGEE. Derived quantities from these spectra (radial velocities, effective temperatures, surface gravities, and metallicities) are also included.DR10 also roughly doubles the number of BOSS spectra over those included in the ninth data release. DR10 includes a total of 1,507,954 BOSS spectra, comprising 927,844 galaxy spectra; 182,009 quasar spectra; and 159,327 stellar spectra, selected over 6373.2 square degrees.

(Abridged) The Census of High- and Medium-mass Protostars (CHaMP) is the first large-scale (280 degree<l<300 degree, -4 degree<b<2 degree), unbiased, sub-parsec resolution survey of Galactic molecular clumps and their embedded stars. Barnes et al. (2011) presented the source catalog of ~300 clumps based on HCO+(1-0) emission, used to estimate masses M. Here we use archival mid-infrared to mm continuum data to construct spectral energy distributions. Fitting two-temperature grey-body models, we derive bolometric luminosities, L. We find the clumps have 10Lsun<L<1E6.5Lsun and 0.1<L/M<1E3, consistent with theoretical expectations of a clump population that spans a range of instantaneous star formation efficiencies from 0 to ~50%. We thus expect L/M to be a useful, strongly-varying indicator of clump evolution during the star cluster formation process. We find correlations of the ratio of warm to cold component fluxes and of cold component temperature with L/M. We also find a near linear relation between L/M and Spitzer-IRAC specific intensity (surface brightness), which may thus also be useful as a star formation efficiency indicator. The lower bound of the clump $L/M$ distribution suggests the star formation efficiency per free-fall time is epsilon<0.2. We do not find strong correlations of L/M with mass surface density, velocity dispersion or virial parameter. We find a linear relation between L and L_{HCO+(1-0}}, although with large scatter for any given individual clump. Fitting together with extragalactic systems, the linear relation still holds, extending over 10 orders of magnitude in luminosity. The complete nature of the CHaMP survey over a several kiloparsec-scale region allows us to derive a measurement at an intermediate scale bridging those of individual clumps and whole galaxies.

We report the discovery of a candidate brown dwarf or a very low mass stellar companion (MARVELS-5b) to the star HIP 67526 from the Multi-object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The radial velocity curve for this object contains 31 epochs spread over 2.5 years. Our Keplerian fit using a Markov Chain Monte Carlo approach, reveals that the companion has an orbital period of $90.2695^{+0.0188}_{-0.0187}$ days, an eccentricity of $0.4375 \pm 0.0040$ and a semi-amplitude of $2948.14^{+16.65}_{-16.55}$ m s$^{-1}$. Using additional high-resolution spectroscopy, we find the host star has an effective temperature $T_{\rm{eff}}=6004 \pm 34$ K, a surface gravity $\log g$ [cgs] $=4.55 \pm 0.17$ and a metallicity [Fe/H] $=+0.04 \pm 0.06$. The stellar mass and radius determined through the empirical relationship of Torres et al. (2010), yields 1.10$\pm$0.09 $M_{\sun}$ and 0.92$\pm$0.19 $R_{\sun}$. The minimum mass of MARVELS-5b is $65.0 \pm 2.9 M_{Jup}$, indicating that it is likely to be either a brown dwarf or a very low mass star, thus occupying a relatively sparsely-populated region of the mass function of companions to solar-type stars. The distance to this system is 101$\pm$10 pc from the astrometric measurements of Hipparcos. No stellar tertiary is detected in the high-contrast images taken by either FastCam lucky imaging or Keck adaptive optics imaging, ruling out any star with mass greater than 0.2$M_{\sun}$ at a separation larger than 40 AU.

We report the discovery of KELT-3b, a moderately inflated transiting hot Jupiter with a mass of 1.477 (-0.067, +0.066) M_J, and radius of 1.345 +/- 0.072 R_J, with an orbital period of 2.7033904 +/- 0.000010 days. The host star, KELT-3, is a V=9.8 late F star with M_* = 1.278 (-0.061, +0.063) M_sun, R_* = 1.472 (-0.067, +0.065) R_sun, T_eff = 6306 (-49, +50) K, log(g) = 4.209 (-0.031, +0.033), and [Fe/H] = 0.044 (-0.082, +0.080), and has a likely proper motion companion. KELT-3b is the third transiting exoplanet discovered by the KELT survey, and is orbiting one of the 20 brightest known transiting planet host stars, making it a promising candidate for detailed characterization studies. Although we infer that KELT-3 is significantly evolved, a preliminary analysis of the stellar and orbital evolution of the system suggests that the planet has likely always received a level of incident flux above the empirically-identified threshold for radius inflation suggested by Demory & Seager (2011).

We make a reanalysis on the issue of neutrino speed anomaly by taking into account the newly reported data from the ICARUS experiment and other CNGS collaborations. We examine the consequence of the Lorentz violation on the neutrino speed in a new framework of standard model supplement (SMS), and find that the Lorentz violating parameters are constrained at least one order stronger than that of the earlier OPERA report. The combination with other phenomenological considerations puts more stringent constraints on the Lorentz violation of neutrinos.

We describe the discovery of a likely brown dwarf (BD) companion with a minimum mass of 31.7 +/- 2.0 M_Jup to GSC 03546-01452 from the MARVELS radial velocity survey, which we designate as MARVELS-6b. For reasonable priors, our analysis gives a probability of 72% that MARVELS-6b has a mass below the hydrogen-burning limit of 0.072 M_Sun, and thus it is a high-confidence BD companion. It has a moderately long orbital period of 47.8929 +0.0063/-0.0062 days with a low eccentricty of 0.1442 +0.0078/-0.0073, and a semi-amplitude of 1644 +12/-13 m/s. Moderate resolution spectroscopy of the host star has determined the following parameters: T_eff = 5598 +/- 63, log g = 4.44 +/- 0.17, and [Fe/H] = +0.40 +/- 0.09. Based upon these measurements, GSC 03546-01452 has a probable mass and radius of M_star = 1.11 +/- 0.11 M_Sun and R_star = 1.06 +/- 0.23 R_Sun with an age consistent with less than ~6 Gyr at a distance of 219 +/- 21 pc from the Sun. Although MARVELS-6b is not observed to transit, we cannot definitively rule out a transiting configuration based on our observations. There is a visual companion detected with Lucky Imaging at 7.7 arcsec from the host star, but our analysis shows that it is not bound to this system. The minimum mass of MARVELS-6b exists at the minimum of the mass functions for both stars and planets, making this a rare object even compared to other BDs.

The mass domain where massive extrasolar planets and brown dwarfs overlap is still poorly understood due to the paucity of brown dwarfs orbiting close to solar-type stars, the so-called brown dwarf desert. In this paper we collect all of available data about close brown dwarfs around solar type stars and their host stars from literature and study the demographics of the brown dwarf desert. The data clearly show a short period and a medium mass gap in the brown dwarf period-mass distribution diagram ($ 35<m\sin i<55 M_{\rm Jup}$ and $P<100$ days), representing the "driest land" in the brown dwarf desert. Observation biases are highly unlikely to cause this gap due to its short period and medium mass, of which brown dwarfs can be easily detected by previous RV surveys. Brown dwarfs above and below this gap have significantly different eccentricity distribution, which not only confirms that this gap is real, but also implies that they may have different origins. Our further statistical study indicates that brown dwarfs below this gap may primarily form in the protoplanetary disk through disk gravitational instability, while brown dwarfs above this gap may dominantly form like a stellar binary through molecular cloud fragmentation. Our discoveries have offered important insights about brown dwarf formation mechanisms and their possible relationships with planet and star formation.

We present an eccentric, short-period brown dwarf candidate orbiting the active, slightly evolved subgiant star TYC 2087-00255-1, which has effective temperature T_eff = 5903+/-42 K, surface gravity log (g) = 4.07+/-0.16 (cgs), and metallicity [Fe/H] = -0.23+/-0.07. This candidate was discovered using data from the first two years of the Multi-object APO Radial Velocity Exoplanets Large-area Survey (MARVELS), which is part of the third phase of Sloan Digital Sky Survey. From our 38 radial velocity measurements spread over a two-year time baseline, we derive a Keplerian orbital fit with semi-amplitude K=3.571+/-0.041 km/s, period P=9.0090+/-0.0004 days, and eccentricity e=0.226+/-0.011. Adopting a mass of 1.16+/-0.11 Msun for the subgiant host star, we infer that the companion has a minimum mass of 40.0+/-2.5 M_Jup. Assuming an edge-on orbit, the semimajor axis is 0.090+/-0.003 AU. The host star is photometrically variable at the \sim1% level with a period of \sim13.16+/-0.01 days, indicating that the host star spin and companion orbit are not synchronized. Through adaptive optics imaging we also found a point source 643+/-10 mas away from TYC 2087-00255-1, which would have a mass of 0.13 Msun if it is physically associated with TYC 2087-00255-1 and has the same age. Future proper motion observation should be able to resolve if this tertiary object is physically associated with TYC 2087-00255-1 and make TYC 2087-00255-1 a triple body system. Core Ca II H and K line emission indicate that the host is chromospherically active, at a level that is consistent with the inferred spin period and measured v_{rot}*sin i, but unusual for a subgiant of this T_eff. This activity could be explained by ongoing tidal spin-up of the host star by the companion.

The Sloan Digital Sky Survey III (SDSS-III) presents the first spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS). This ninth data release (DR9) of the SDSS project includes 535,995 new galaxy spectra (median z=0.52), 102,100 new quasar spectra (median z=2.32), and 90,897 new stellar spectra, along with the data presented in previous data releases. These spectra were obtained with the new BOSS spectrograph and were taken between 2009 December and 2011 July. In addition, the stellar parameters pipeline, which determines radial velocities, surface temperatures, surface gravities, and metallicities of stars, has been updated and refined with improvements in temperature estimates for stars with T_eff<5000 K and in metallicity estimates for stars with [Fe/H]>-0.5. DR9 includes new stellar parameters for all stars presented in DR8, including stars from SDSS-I and II, as well as those observed as part of the SDSS-III Sloan Extension for Galactic Understanding and Exploration-2 (SEGUE-2). The astrometry error introduced in the DR8 imaging catalogs has been corrected in the DR9 data products. The next data release for SDSS-III will be in Summer 2013, which will present the first data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) along with another year of data from BOSS, followed by the final SDSS-III data release in December 2014.

We report the discovery via radial velocity of a short-period (P = 2.430420 \pm 0.000006 days) companion to the F-type main sequence star TYC 2930-00872-1. A long-term trend in the radial velocities indicates the presence of a tertiary stellar companion with $P > 2000$ days. High-resolution spectroscopy of the host star yields T_eff = 6427 +/- 33 K, log(g) = 4.52 +/- 0.14, and [Fe/H]=-0.04 +/- 0.05. These parameters, combined with the broad-band spectral energy distribution and parallax, allow us to infer a mass and radius of the host star of M_1=1.21 +/- 0.08 M_\odot and R_1=1.09_{-0.13}^{+0.15} R_\odot. We are able to exclude transits of the inner companion with high confidence. The host star's spectrum exhibits clear Ca H and K core emission indicating stellar activity, but a lack of photometric variability and small v*sin(I) suggest the primary's spin axis is oriented in a pole-on configuration. The rotational period of the primary from an activity-rotation relation matches the orbital period of the inner companion to within 1.5 \sigma, suggesting they are tidally locked. If the inner companion's orbital angular momentum vector is aligned with the stellar spin axis, as expected through tidal evolution, then it has a stellar mass of M_2 ~ 0.3-0.4 M_\odot. Direct imaging limits the existence of stellar companions to projected separations < 30 AU. No set of spectral lines and no significant flux contribution to the spectral energy distribution from either companion are detected, which places individual upper mass limits of M < 1.0 M_\odot, provided they are not stellar remnants. If the tertiary is not a stellar remnant, then it likely has a mass of ~0.5-0.6 M_\odot, and its orbit is likely significantly inclined from that of the secondary, suggesting that the Kozai-Lidov mechanism may have driven the dynamical evolution of this system.

TYC 4110-01037-1 has a low-mass stellar companion, whose small mass ratio and short orbital period are atypical amongst solar-like (Teff ~< 6000 K) binary systems. Our analysis of TYC 4110-01037-1 reveals it to be a moderately aged (~<5 Gyr) solar-like star having a mass of 1.07 +/- 0.08 MSun and radius of 0.99 +/- 0.18 RSun. We analyze 32 radial velocity measurements from the SDSS-III MARVELS survey as well as 6 supporting radial velocity measurements from the SARG spectrograph on the 3.6m TNG telescope obtained over a period of ~2 years. The best Keplerian orbital fit parameters were found to have a period of 78.994 +/- 0.012 days, an eccentricity of 0.1095 +/- 0.0023, and a semi-amplitude of 4199 +/- 11 m/s. We determine the minimum companion mass (if sin i = 1) to be 97.7 +/- 5.8 MJup. The system's companion to host star mass ratio, >0.087 +/- 0.003, places it at the lowest end of observed values for short period stellar companions to solar-like (Teff ~< 6000 K) stars. One possible way to create such a system would be if a triple-component stellar multiple broke up into a short period, low q binary during the cluster dispersal phase of its lifetime. A candidate tertiary body has been identified in the system via single-epoch, high contrast imagery. If this object is confirmed to be co-moving, we estimate it would be a dM4 star. We present these results in the context of our larger-scale effort to constrain the statistics of low mass stellar and brown dwarf companions to FGK-type stars via the MARVELS survey.

Building on the legacy of the Sloan Digital Sky Survey (SDSS-I and II), SDSS-III is a program of four spectroscopic surveys on three scientific themes: dark energy and cosmological parameters, the history and structure of the Milky Way, and the population of giant planets around other stars. In keeping with SDSS tradition, SDSS-III will provide regular public releases of all its data, beginning with SDSS DR8 (which occurred in Jan 2011). This paper presents an overview of the four SDSS-III surveys. BOSS will measure redshifts of 1.5 million massive galaxies and Lya forest spectra of 150,000 quasars, using the BAO feature of large scale structure to obtain percent-level determinations of the distance scale and Hubble expansion rate at z<0.7 and at z~2.5. SEGUE-2, which is now completed, measured medium-resolution (R=1800) optical spectra of 118,000 stars in a variety of target categories, probing chemical evolution, stellar kinematics and substructure, and the mass profile of the dark matter halo from the solar neighborhood to distances of 100 kpc. APOGEE will obtain high-resolution (R~30,000), high signal-to-noise (S/N>100 per resolution element), H-band (1.51-1.70 micron) spectra of 10^5 evolved, late-type stars, measuring separate abundances for ~15 elements per star and creating the first high-precision spectroscopic survey of all Galactic stellar populations (bulge, bar, disks, halo) with a uniform set of stellar tracers and spectral diagnostics. MARVELS will monitor radial velocities of more than 8000 FGK stars with the sensitivity and cadence (10-40 m/s, ~24 visits per star) needed to detect giant planets with periods up to two years, providing an unprecedented data set for understanding the formation and dynamical evolution of giant planet systems. (Abridged)

In the electromagnetism of loop quantum gravity, two helicities of a photon have different phase velocities and group velocities, termed as "vacuum birefringence". Two novel phenomenons, "peak doubling" and "de-polarization", are expected to appear for a linearly polarized light from astrophysical sources. We show that the criteria to observe these two phenomenons are the same. Further, from recently observed $\gamma$-ray polarization from Cygnus X-1, we obtain an upper limit $\sim 8.7\times10^{-12}$ for Lorentz-violating parameter $\chi$, which is the most firm constraint from well-known systems. We also suggest to analyze possible existence of "peak doubling" through Fermi LAT GRBs.

Doubly special relativity (DSR), with both an invariant velocity and an invariant length scale, elegantly preserves the principle of relativity between moving observers, and appears as a promising candidate of the quantum theory of gravity. We study the modifications of photon gas thermodynamics in the framework of DSR with an invariant length $|\lambda|$, after properly taking into account the effects of modified dispersion relation, upper bounded energy-momentum space, and deformed integration measure. We show that with a positive $\lambda$, the grand partition function, the energy density, the specific heat, the entropy, and the pressure are smaller than those of special relativity (SR), while the velocity of photons and the ratio of pressure to energy are larger. In contrast, with a negative $\lambda$, the quantum gravity effects show up in the opposite direction. However, these effects only manifest themselves significantly when the temperature is larger than $10^{-3} E_{\rm P}$. Thus, DSR can have considerable influence on the early universe in cosmological study.

We study the self-gravitational effect on the equation of state (EoS) of Bose and Fermi gases in thermal equilibrium at the end of reheating, the period after quark-hadron transition and before Big Bang Nucleosynthesis (BBN). After introducing new grand canonical partition functions based on the work of Uhlenbeck and Gropper, we notice some interesting features of the newly developed EoSs with distinct behaviors of relativistic and non-relativistic gases under self-gravity. The usual negligence of the self-gravitational effect when solving the background expansion of the early universe is justified with numerical results, showing the magnitude of the self-gravitational modification of the state constant to be less than $O(10^{-78})$. This helps us to clarify the background thermal evolution of the primordial patch. Such clarification is crucial in testing gravity theories, evaluating inflation models and determining element abundances in BBN.

Flux-limited and volume-limited galaxy samples are constructed from SDSS data releases DR4, DR6 and DR7 for statistical analysis. The two-point correlation functions $\xi(s)$, monopole of three-point correlation functions $\zeta_0$, projected two-point correlation function $w_p$ and pairwise velocity dispersion $\sigma_{12}$ are measured to test if galaxy samples are fair for these statistics. We find that with increment of sky coverage of SDSS, $\xi(s)$ of flux-limited sample is extremely robust and insensitive to local structures at low redshift. But for volume-limited samples fainter than $L^*$ at large scales $s>\sim 10\hmpc$, deviation of $\xi(s)$ and $\zeta_0$ of DR7 to those of DR4 and DR6 increases with larger absolute magnitude. In the weakly nonlinear regime, there is no agreement between $\zeta_0$ of different data releases in all luminosity bins. Furthermore, $w_p$ of volume-limited samples of DR7 in luminosity bins fainter than $-M_{r,0.1}=[18.5,19.5]$ are significantly larger, and $\sigma_{12}$ of the two faintest volume-limited samples of DR7 display very different scale dependence than results of DR4 and DR6. Our findings call for cautions in understanding clustering analysis results of SDSS faint galaxy samples, and higher order statistics of SDSS volume-limited samples in the weakly nonlinear regime. The first zero-crossing points of $\xi(s)$ of volume-limited samples are also investigated and discussed.

The occurrence of digits one through nine as the leftmost nonzero digit of numbers from real world sources is often not uniformly distributed, but instead, is distributed according to a logarithmic law, known as Benford's law. Here, we investigate systematically the mantissa distributions of some pulsar quantities, and find that for most quantities their first digits conform to this law. However, the barycentric period shows significant deviation from the usual distribution, but satisfies a generalized Benford's law roughly. Therefore pulsars can serve as an ideal assemblage to study the first digit distributions of real world data, and the observations can be used to constrain theoretical models of pulsar behavior.

We show that black holes can be quantized in an intuitive and elegant way with results in agreement with conventional knowledge of black holes by using Bohr's idea of quantizing the motion of an electron inside the atom in quantum mechanics. We find that properties of black holes can be also derived from an Ansatz of quantized entropy $\Delta S=4\pi k {\Delta R / \lambdabar}$, which was suggested in a previous work to unify the black hole entropy formula and Verlinde's conjecture to explain gravity as an entropic force. Such an Ansatz also explains gravity as an entropic force from quantum effect. This suggests a way to unify gravity with quantum theory. Several interesting and surprising results of black holes are given from which we predict the existence of primordial black holes ranging from Planck scale both in size and energy to big ones in size but with low energy behaviors.

Lorentz violation (LV) is predicted by some quantum gravity theories, where photon dispersion relation is modified, and the speed of light becomes energy-dependent. Consequently, it results in a tiny time delay between high energy photons and low energy ones. Very high energy (VHE) photon emissions from cosmological distance can amplify these tiny LV effects into observable quantities. Here we analyze four VHE $\gamma$-ray bursts (GRBs) from Fermi observations, and briefly review the constraints from three TeV flares of active galactic nuclei (AGNs) as well. One step further, we present a first robust analysis of VHE GRBs taking the intrinsic time lag caused by sources into account, and give an estimate to quantum gravity energy $\sim 2 \times 10^{17}$ GeV for linear energy dependence, and $\sim 5 \times 10^9$ GeV for quadratic dependence. However, the statistics is not sufficient due to the lack of data, and further observational results are desired to constrain LV effects better.

We study entropic force effects on black holes and photons. We find that application of an entropic analysis restricts the radial change $\Delta R$ of a black hole of radius $R_{\mathrm{H}}$, due to a test particle of a Schwartzchild radius $R_{h}$ moving towards the black hole by $\Delta x$ near black body surface, to be given by a relation $R_{\mathrm{H}} \Delta R= R_h \Delta x/2$, or ${\Delta R}/{\lambdabar_M} = {\Delta x}/{2 \lambdabar_m}$. We suggest a new rule regarding entropy changes in different dimensions, $\Delta S= 2\pi k D \Delta l /\lambdabar,$ which unifies Verlinde's conjecture and the black hole entropy formula. We also propose to extend the entropic force idea to massless particles such as a photon. We find that there is an entropic force on a photon of energy $E_\gamma$, with $F=G M m_{\gamma}/R^2$, and therefore the photon has an effective gravitational mass $m_\gamma = E_\gamma/c^2$.

We obtain modified dispersion relations by requiring the vanishing of determinant of inverse of modified photon propagators in Lorentz invariance violation (LIV) theory. Inspired by these dispersion relations, we give a more general dispersion relation with less assumption and apply it to the recent observed gamma-ray burst GRB090510 to extract various constraints on LIV parameters. We find that the constraint on quantum gravity mass is slightly larger than the Planck mass but is consistent with the other recent observations, so the corresponding LIV coefficient $\xi_1$ has reached the natural order ($o(1)$) as one expects. From our analysis, the linear LIV corrections to photon group velocity might be not excluded yet.

We present detailed clustering analysis of a large K-band selected local galaxy sample, which is constructed from the 2MASS and the SDSS and consists of $82,486$ galaxies with $10 < K < 13.5$ and $0.01 < z < 0.1$. The two-point correlation function of the magnitude-limited sample in real space at small scales is well described by a power law $\xi(r)=(r/6.44\pm0.23)^{-1.81\pm0.02}$. The pairwise velocity dispersion is derived from the anisotropic two-point correlation function and find the dispersion $\sigma_{12}=685\pm 17\kms$ if its scale invariance is assumed, which is larger than values measured in optical bands selected galaxy samples. We further investigate the dependence of the two-point correlation function and the $\sigma_{12}$ on the $g-r$ color and the $K$-band luminosity, obtain similar results to previous works in optical bands. Comparing a mock galaxy sample with our real data indicates that the semi-analytical model can not mimic the $\sigma_{12}$ in observation albeit it can approximate the two-point correlation function within measurement uncertainties.

The formation of ultra-compact X-ray binaries (UCXBs) has not been well understood. Previous works show that ultra-short orbital periods ($<1$ hr) may be reached through mass transfer driven by magnetic braking in normal low/intermediate-mass X-ray binaries (L/IMXBs) only for an extremely small range of initial binary parameters, which makes it difficult to account for the rather large population of UCXBs. In this paper we report the calculated results on mass transfer processes in L/IMXBs with a circum-binary disk. We show that when the orbital angular momentum loss due to a circum-binary disk is included, ultra-short orbital periods could be reached for a relatively wide range of initial binary parameters. The results of our binary models suggest an alternative formation channel for UCXBs.

The bifurcation period in low-mass X-ray binaries is the initial orbital pe- riod which separates the formation of converging systems (which evolve with decreasing orbital periods until the donor becomes degenerate) from the diverging systems (which evolve with increasing orbital periods until the donor star loses its envelope and a wide detached binary is formed). We calculate systematically the bifurcation periods of binary systems with a 1.4M_\sun neutron star and a 0.5-2M_\sun donor star, taking into account different kinds of magnetic braking and mass loss mechanisms. Our results show that the saturated magnetic braking can considerably decrease the values of bifurcation period compared to the traditional magnetic braking, while the influence of mass loss mechanisms on bifurcation periods is quite weak. We also develop a semi-analytical method to compute the bifurcation period, the result of which agrees well with the numerical method in the leading order.

The muon charge ratio of ultrahigh energy cosmic rays may provide information to detect the composition of the primary cosmic rays. We propose to extract the charge information of high energy muons in very inclined extensive air showers by analyzing their relative lateral positions in the shower transverse plane.

Giant flares from soft gamma-ray repeaters (SGRs) are one of the most violent phenomena in neutron stars. Quasi-periodic oscillations (QPOs) with frequencies ranging from 18 to 1840 Hz have been discovered in the tails of giant flares from two SGRs, and were ascribed to be seismic vibrations or torsional oscillations of magnetars. Here we propose an alternative explanation for the QPOs in terms of standing sausage mode oscillations of flux tubes in the magnetar coronae. We show that most of the QPOs observed in SGR giant flares could be well accounted for except for those with very high frequencies (625 and 1840 Hz).

We study a five-dimensional perfect fluid coupled with Kaluza-Klein (KK) gravity. By dimensional reduction, a modified form of Maxwell's equation is obtained, which is relevant to the equation of state of the source. Since the relativistic magnetohydrodynamics (MHD) and the 3-dimensional formulation are widely used to study space matter, we derive the modified Maxwell's equations and relativistic MHD in 3+1 form. We then take an ideal Fermi gas as an example to study the modified effect, which can be visible under high density or high energy condition, while the traditional Maxwell's equation can be regarded as a result in the low density and low temperature limit. We also indicate the possibility to test the state-relevant effect of KK theory in a telluric laboratory.

A new model for the description of the behaviour of strangelets in the Earth's atmosphere is presented. Strangelet fission induced by collision with air nuclei is included. It is shown that strangelets with certain parameters of initial mass and energy may reach depths near sea level, which can be examined by ground-based experiments.

We propose to extract the charge information of high energy muons in very inclined extensive air showers by analyzing their relative lateral positions in the shower transverse plane. We calculate the muon lateral deviation under the geomagnetic field and compare it to dispersive deviations from other causes. By our criterion of resolvability, positive and negative muons with energies above $10^4$ GeV will be clearly separated into two lobes if the shower zenith angle is larger than $70^\circ$. Thus we suggest a possible approach to measure the $\mu^+ / \mu^-$ ratio for high energy muons.

The neutron fraction in the very high energy cosmic rays near the Greisen-Zatsepin-Kuzmin (GZK) cutoff energy is analyzed by taking into account the time dilation effect of the neutron decays and also the pion photoproduction behaviors above the GZK cutoff. We predict a non-trivial neutron fraction above the GZK cutoff and a negligibly small neutron fraction below. However, there should be a large antiproton fraction in the high energy cosmic rays below the GZK cutoff in several existing models for the observed cosmic-ray events above and near the GZK cutoff. Such a large antiproton fraction can manifest itself by the muon charge ratio $\mu^+/\mu^-$ in the collisions of the primary nucleon cosmic rays with the atmosphere, if there is no neutron contribution. We suggest to use the muon charge ratio as one of the information to detect the composition of the primary cosmic rays near or below the GZK cutoff.

We discuss the relations between the mixing angles and the mixing matrices of quarks and leptons. With Raidal's numerical relations, we parametrize the lepton mixing (PMNS) matrix with the parameters of the quark mixing (CKM) matrix, and calculate the products of $V_{\mathrm{CKM}}U_{\mathrm{PMNS}}$ and $U_{\mathrm{PMNS}}V_{\mathrm{CKM}}$. Also, under the conjectures $V_{\mathrm{CKM}}U_{\mathrm{PMNS}}=U_{\mathrm{bimax}}$ or $U_{\mathrm{PMNS}}V_{\mathrm{CKM}}=U_{\mathrm{bimax}}$, we get the PMNS matrix naturally, and test Raidal's relations in these two different versions. The similarities and the differences between the different versions are discussed in detail.

We apply Koide's mass relation of charged leptons to neutrinos and quarks, with both the normal and inverted mass schemes of neutrinos discussed. We introduce the parameters $k_{\nu}$, $k_u$ and $k_d$ to describe the deviations of neutrinos and quarks from Koide's relation, and suggest a quark-lepton complementarity of masses such as $ k_{l}+k_{d} \approx k_{\nu}+k_{u} \approx 2$. The masses of neutrinos are determined from the improved relation, and they are strongly hierarchical (with the different orders of magnitude of $10^{-5} eV$, $10^{-3} eV$, and $10^{-2} eV$).

We present a unified parametrization of quark and lepton mixing matrices. By using some simple relations between the mixing angles of quarks and leptons, i.e., the quark-lepton complementarity, we parametrize the lepton mixing matrix with the Wolfenstein parameters $\lambda$ and $A$ of the quark mixing matrix. It is shown that the Wolfenstein parameter $\lambda$ can measure both the deviation of the quark mixing matrix from the unit matrix, and the deviation of the lepton mixing matrix from the exactly bimaximal mixing pattern.

We propose a method to investigate the scenario that cosmic relic neutrinos are highly clustered around stars and galaxies, or dark-matter clusters, rather than uniformly distributed in the universe. Such a scenario can be detected or constrained by the interaction of high energy cosmic ray protons and nuclei with the cosmic neutrinos. There should be observable signature in the energy spectra of cosmic ray protons and nuclei for a neutrino clustering factor beyond $10^{13}$. We provide a relation on the signature onset positions between proton and nuclei spectra, and discuss possible support from existing experiments. It is also suggested that the relative abundance of cosmic ray nuclei may detect or constrain the cosmic neutrinos with smaller clustering.

The NuTeV Collaboration reported a value of $\sin^{2}\theta_{w}$ measured in neutrino-nucleon deep inelastic scattering, and found that the value is three standard deviations from the world average value of other electroweak measurements. If this result cannot be explained within conventional physics, it must imply some novel physics beyond the standard model. We report the correction from the asymmetric strange-antistrange sea by using both the light-cone baryon-meson fluctuation model and the chiral quark model, and show that a significant part of the NuTeV anomaly can be explained by the strange-antistrange asymmetry.

We calculate the strange quark and antiquark distributions of the nucleon by using the effective chiral quark model, and find that the strange-antistrange asymmetry can bring a contribution of about 60--100% to the NuTeV deviation of $\sin^{2}\theta_{w}$ from the standard value measured in other electroweak processes. The results are insensitive to different inputs. The light-flavor quark asymmetry of $\bar{d}-\bar{u}$ is also investigated and found to be consistent with the experimental measurements. Therefore the chiral quark model provides a successful picture to understand the NuTeV anomaly, as well as the light-flavor quark asymmetry and the proton spin problem in previous studies.

The neutrino mixing matrix is expanded in powers of a small parameter $\lambda$, which approximately equals to 0.1. The meaning of every order of the expansion is discussed respectively, and the range of $\lambda$ is carefully calculated. We also present some applications of this new parametrization, such as to the expression of the Jarlskog parameter $J$, in which the simplicities and advantages of this parametrization are shown.

We confirmed the following idea by numerical calculation: The extended structure of baryons may make the nuclear equation of state stiffer at higher density while keeping the compression modulus for normal nuclear matter around its empirical value $K=240MeV$.The model built in this wayy fits all empirical data for normal nuclear matter and gives the mass limit of neutron stars compatible with the observation.