Strong gravitational lensing has been a powerful probe of cosmological models and gravity. To date, constraints in either domain have been obtained separately. We propose a new methodology through which the cosmological model, specifically the Hubble constant, and post-Newtonian parameter can be simultaneously constrained. Using the time-delay cosmography from strong lensing combined with the stellar kinematics of the deflector lens, we demonstrate the Hubble constant and post-Newtonian parameter are incorporated in two distance ratios which reflect the lensing mass and dynamical mass, respectively. Through the reanalysis of the four publicly released lenses distance posteriors from the H0LiCOW collaboration, the simultaneous constraints of Hubble constant and post-Newtonian parameter are obtained. Our results suggests no deviation from the General Relativity, $\gamma_{\texttt{PPN}}=0.87^{+0.19}_{-0.17}$ with a Hubble constant favors the local Universe value, $H_0=73.65^{+1.95}_{-2.26}$ km s$^{-1}$ Mpc$^{-1}$. Finally, we forecast the robustness of gravity tests by using the time-delay strong lensing for constraints we expect in the next few years. We find that the joint constraint from 40 lenses are able to reach the order of $7.7\%$ for the post-Newtonian parameter and $1.4\%$ for Hubble constant.

We present a general model-independent approach to study the coupled dark energy and the string Swampland criteria. We show how the dark sector interaction is degenerated with the equation of state of dark energy in the context of the expansion of the Universe. With priors for either of them, the dynamics of dark energy and the dark sector interactions can be reconstructed together with the bounds of the Swampland criteria. Combining cosmic chronometers, baryon acoustic oscillation (BAO), and Type Ia supernovae our results suggest a mild $1 \sigma$ significance of dark sector interactions at low redshift for the coupled quintessence. The Lyman-$\alpha$ BAO at $z=2.34$ leads a $2 \sigma$ signal of non-zero interactions at high redshift. The implications for coupled quintessence are discussed.

Local determinations of the Hubble constant $H_0$ favour a higher value than Planck based on CMB and $\Lambda$CDM. Through a model-independent expansion, we show that low redshift ($z \lesssim 0.7$) data comprising baryon acoustic oscillations (BAO), cosmic chronometers and Type Ia supernovae has a preference for Quintessence models that lower $H_0$ relative to $\Lambda$CDM. In addition, we confirm that an exponential coupling to dark matter cannot alter this conclusion in the same redshift range. Our results leave open the possibility that a coupling in the matter-dominated epoch, potentially even in the dark ages, may yet save $H_0$ from sinking in the string theory Swampland.

We report pulsar timing observations carried out in L-band with NTSC's 40-meter Haoping Radio Telescope (HRT), which was constructed in 2014. The observations were carried out using the pulsar machine we developed. Timing observations toward millisecond pulsar J0437-4715 obtains a timing residual (r.m.s) of 397ns in the time span of 284 days. And our observations successfully detected Crab pulsar's glitch that happened on July 23rd, 2019.

Strong gravitational lensing by galaxies provides us with a unique opportunity to understand the nature of gravity on galactic and extra-galactic scales. In this paper, we propose a new multimessenger approach using data from both gravitational wave (GW) and the corresponding electromagnetic (EM) counterpart to infer the constraint of the modified gravity (MG) theory denoted by the scale dependent phenomenological parameter. To demonstrate the robustness of this approach, we calculate the time-delay predictions by choosing various values of the phenomenological parameters and then compare them with that from the general relativity (GR). For the third generation ground-based GW observatory, with one typical strongly lensed GW+EM event, and assuming that the dominated error from the stellar velocity dispersions is 5\%, the GW time-delay data can distinguish an 18\% MG effect on a scale of tens of kiloparsecs with a $68\%$ confidence level. Assuming GR and a Singular Isothermal Sphere mass model, there exists a simplified consistency relationship between time-delay and imaging data. This relationship does not require for the velocity dispersion measurement, and hence can avoid major uncertainties. By using this relationship, the multimessenger approach is able to distinguish an $8\%$ MG effect. Our results show that the GW multimessenger approach can play an important role in revealing the nature of gravity on the galactic and extra-galactic scales.

We consider a low redshift $(z<0.7)$ cosmological dataset comprising megamasers, cosmic chronometers, type Ia SNe and BAO, which we bin according to their redshift. For each bin, we read the value of $H_0$ by fitting directly to the flat $\Lambda$CDM model. Doing so, we find that $H_0$ descends with redshift, allowing one to fit a line with a \textit{non-zero} slope of statistical significance $2.1 \, \sigma$. Our results are in accord with a similar descending trend reported by the H0LiCOW collaboration. If substantiated going forward, early Universe solutions to the Hubble tension will struggle explaining this trend.

In this paper, we analyse the light variation of KIC 8840638 using high-precision time-series data delivered by the $Kepler$ mission. The analysis reveals that the target KIC 8840638 is an eclipsing binary with a $\delta$ Scuti component, not a pure single $\delta$ Scuti star previously known. The frequency analysis of the high-precision short cadence light curve reveals 95 significant frequencies, most of which lie in the frequency range of 23$-$32 d$^{-1}$. Among them, seven independent frequencies are detected in the typical frequency range of the $\delta$ Scuti variables. Besides, the orbital frequency $f_{orb}$ (=0.320008 d$^{-1}$) and its harmonics are also detected directly in the frequency spectrum. The binary modelings derived from the Wilson-Devinney code indicate the binary parameters are not affected by the pulsations, and the binary system is in semi-detached configure with a mass ratio of 1.6, which is higher than that of all other $\delta$ Scuti eclipsing binaries. The locations of the components of KIC 8840638 indicate the primary star lies in the $\delta$ Scuti instability region below the terminal-age main-sequence, and the secondary component might be an evolved star.

The cold dark matter model predicts that dark matter haloes are connected by filaments. Direct measurements of the masses and structure of these filaments are difficult, but recently several studies have detected these dark-matter-dominated filaments using weak lensing. Here we study the efficiency of galaxy formation within the filaments by measuring their total mass-to-light ratios and stellar mass fractions. Specifically, we stack pairs of Luminous Red Galaxies (LRGs) with a typical separation on the sky of $8 h^{-1}$ Mpc. We stack background galaxy shapes around pairs to obtain mass maps through weak lensing, and we stack galaxies from the Sloan Digital Sky Survey (SDSS) to obtain maps of light and stellar mass. To isolate the signal from the filament, we construct two matched catalogues of physical and non-physical (projected) LRG pairs, with the same distributions of redshift and separation. We then subtract the two stacked maps. Using LRG pair samples from the BOSS survey at two different redshifts, we find that the evolution of the mass in filament is consistent with the predictions from perturbation theory. The filaments are not entirely dark: their mass-to-light ratios ($M/L = 351\pm87$ in solar units in the $r$-band) and stellar mass fractions ($M_{\rm stellar}/M = 0.0073\pm0.0020$) are consistent with the cosmic values (and with their redshift evolutions)

Recently constructed high-redshift Hubble diagrams of supernovae, quasars and gamma-ray bursts show $\sim4 \, \sigma$ deviations from flat $\Lambda$CDM, where the discrepancies are observed in cosmographic expansions. In this note, we show that existing approaches to high-redshift cosmography can lead to "tensions" that are simply an artifact of the expansions, even if the underlying model is fully consistent with flat $\Lambda$CDM. Thus, in this note we revisit the $\sim 4 \, \sigma$ tensions to ascertain if results have been impacted by the cosmography. By directly fitting respective datasets to flat $\Lambda$CDM we find that the tension is real, but the quasar data, which has been calibrated by external Type Ia supernovae has a 7:3 preference for a Universe with no dark energy over one described by $\Lambda$CDM with canonical (Planck) values. Results based on cosmography mask this feature.

Recent work has shown that modified gravitational wave (GW) propagation can be a powerful probe of dark energy and modified gravity, specific to GW observations. We use the technique of Gaussian processes, that allows the reconstruction of a function from the data without assuming any parametrization, to measurements of the GW luminosity distance from simulated joint GW-GRB detections, combined with measurements of the electromagnetic luminosity distance by simulated DES data. For the GW events we consider both a second-generation LIGO/Virgo/Kagra (HVLKI) network, and a third-generation detector such as the Einstein Telescope. We find that the HVLKI network at target sensitivity, with $O(15)$ neutron star binaries with electromagnetic counterpart, could already detect deviations from GR at a level predicted by some modified gravity models, and a third-generation detector such as ET would have a remarkable discovery potential. We discuss the complementarity of the Gaussian processes technique to the $(\Xi_0,n)$ parametrization of modified GW propagation.

Most measurements of mass in Astronomy that use kinematics of stars or gas rely on assumptions of equilibrium that are often hard to verify. Instead, we develop a novel idea that uses the clustering in action space, as a probe of underlying gravitational potential: the correct potential should maximize small-scale clustering in the action space. We provide a first-principle derivation of likelihood using the two-point correlation function in action space, and test it against simulations of stellar streams. We then apply this method to the 2nd data release of Gaia, and use it to measure the radial force fraction $f_h$ and logarithmic slope $\alpha$ of dark matter halo profile. We investigate stars within 9-11 kpc and 11.5-15 kpc from Galactic centre, and find $(f_h,\alpha)= (0.391\pm 0.009, 1.835\pm 0.092) $ and $(0.351\pm 0.012,1.687\pm 0.079)$, respectively. We also confirm that the set of parameters that maximize the likelihood function do correspond to the most clustering in the action space. The best-fit circular velocity curve for Milky Way potential is consistent with past measurements (although it is $\sim$ 5-10\% lower than previous methods that use masers or globular clusters). Our work provides a clear demonstration of the full statistical power that lies in the full phase space information, relieving the need for {\it ad hoc} assumptions such as virial equilibrium, circular motion, or steam-finding algorithms.

In this paper, we present an analysis of the pulsating behavior of $Kepler$ target KIC 10284901. The Fourier transform of the high-precision light curve reveals seven independent frequencies for its light variations. Among them, the first two frequencies are main pulsation modes: F0 = 18.994054(1) $\rm{day^{-1}}$ and F1 = 24.335804(4) $\rm{day^{-1}}$, the ratio of F0/F1 = 0.7805 classifies this star as a double-mode high-amplitude $\delta$ Scuti (HADS) star; another two frequencies, $f_{m1}$ = 0.4407 day$^{-1}$ and $f_{m2}$ = 0.8125 day$^{-1}$, are detected directly, and the modulations of $f_{m1}$ and $f_{m2}$ to F0 and F1 modes (seen as quintuplet structures centered on these two modes in the frequency spectrum) are also found. This is the first detection of a double-modulation effect in the HADS stars. The features of the frequency pattern and the ratio ($f_{m1}$/$f_{m2}$ $\approx$ 1:2), as well as the cyclic variation of amplitude of the two dominant pulsation modes, which seem to be similar to that in Blazhko RR Lyrae stars, indicate this modulation might be related to the Blazhko effect. A preliminary analysis suggests that KIC 10284901 is in the bottom of the HADS instability strip and situated in the main sequence.

We carried out photometric observations of the SX Phe star BL Cam in 2014, 2017 and 2018 using Nanshan 1-m telescope. In addition to the dominated frequency of 25.5790(3) cd$^{-1}$ and its two harmonics, an independent frequency of 25.247 (2) cd$^{-1}$, which is a nonradial mode frequency, was detected from the data in 2014. A total of 123 new times of light maxima were determined from our light curves in the three years, which, together with that published in the literature, were used to analyze the $O$$-$$C$ diagram. The change rate of the main period was derived as (1/P)(dP/dt) = -2.39 (8)$\times$10$^{-8}$ yr$^{-1}$, which is lower than that published in previous literature. A periodical change with a period of 14.01 (9) yr was found in the residuals of the $O$$-$$C$ curve fitting. If it was caused by the light-time effect, BL Cam should be a binary system. The mass of the companion was restricted as low as that of a brown dwarf. No evidence of the triple system suggested by previous authors was shown in our analysis.

Employing multiple pulsars and using an appropriate algorithm to establish ensemble pulsar timescale can reduce the influences of various noises on the long-term stability of pulsar timescale, compared to a single pulsar. However, due to the low timing precision and the significant red noises of some pulsars, their participation in the construction of ensemble pulsar timescale is often limited. Inspired by the principle of solving non-stationary sequence modeling using co-integration theory, we puts forward an algorithm based on the co-integration theory to establish ensemble pulsar timescale. It is found that this algorithm can effectively suppress some noise sources if a co-integration relationship between different pulsar data exist. Different from the classical weighted average algorithm, the co-integration method provides the chances of the pulsar with significant red noises to attend the establishment of ensemble pulsar timescale. Based on the data from the North American Nanohertz Observatory for Gravitational Waves, we found that the co-integration algorithm can successfully reduce several timing noises and improve the long-term stability of the ensemble pulsar timescale.

The cosmic distance duality relation (CDDR) has been test through several astronomical observations in the last years. This relation establishes a simple equation relating the angular diameter ($D_A$) and luminosity ($D_L$) distances at a redshift $z$, $D_LD_A^{-1}(1+z)^{-2}=\eta=1$. However, only very recently this relation has been observationally tested at high redshifts ($z \approx 3.6$) by using luminosity distances from type Ia supernovae (SNe Ia) and gamma ray bursts (GRBs) plus angular diameter distances from strong gravitational lensing (SGL) observations. The results show that no significant deviation from the CDDR validity has been verified. In this work, we test the potentialities of future luminosity distances from gravitational waves (GWs) sources to impose limit on possible departures of CDDR jointly with current SGL observations. The basic advantage of $D_L$ from GWs is being insensitive to non-conservation of the number of photons. By simulating 600, 900 and 1200 data of GWs using the Einstein Telescope (ET) as reference, we derive limits on $\eta(z)$ function and obtain that the results will be at least competitive with current limits from the SNe Ia $+$ GRBs $+$ SGLs analyses.

This paper presents the Charge-Coupled Device time-series photometric observations of the open cluster NGC 1582 and its surrounding field with Johnson B, V and R filters by using Nanshan 1m telescope of Xinjiang Astronomical Observatory. 19 variable stars and 3 variable candidates were detected in a 45 {\times} 48.75 field around the cluster. 12 of the variable stars are newly-discovered variable objects. The physical properties, classifications, and membership of the 22 objects are studied through their light curves, their positions on the color-magnitude diagram, and with the archival data from the Naval Observatory Merged Astrometric Dataset. Among these objects, 5 are eclipsing binary systems, 6 are pulsating variable stars including one known Delta Scuti star and one newly-discovered RR Lyrae star. The distance of the RR Lyrae star is estimated to be 7.9 {\pm} 0.3 kpc, indicates that the star locates far behind the cluster. 4 variable stars are probable members of the cluster, and 13 of the 22 objects are confirmed to be field stars.

In this paper, the high-precision light curves of $Kepler$ target KIC 5950759 are analyzed. The Fourier analysis of the long cadence light curve reveals 3 independent frequencies. Two of them are main pulsation modes: F0 = 14.221373(21) $\rm{d^{-1}}$ and F1 = 18.337249(44) $\rm{d^{-1}}$. The third independent frequency $f_m$ = 0.3193 d$^{-1}$ is found in long cadence data with a signal-to-noise ratio of 6.2. A weak modulation of $f_m$ to F0 and F1 modes (triplet structures centred on F0 and F1) are detected both in long and short cadence data. This is the first detection of the modulation effect in a double-mode HADS star. The most possible cause of the modulation effect in the light curves is amplitude modulation with the star's rotation frequency of 0.3193 d$^{-1}$. The preliminary analysis suggests that KIC 5950759 is in the bottom of the HADS instability strip and likely situated in the main sequence. Spectroscopic observations are necessary to verify the true nature of the modulation terms.

In this paper, we simulate Super-Eddington accreting massive black holes (SEAMBHs) as the candles to probe cosmology for the first time. SEAMBHs have been demonstrated to be able to provide a new tool for estimating cosmological distance. Thus, we create a series of mock data sets of SEAMBHs, especially in the high redshift region, to check their abilities to probe the cosmology. To fulfill the potential of the SEAMBHs on the cosmology, we apply the simulated data to three projects. The first is the exploration of their abilities to constrain the cosmological parameters, in which we combine different data sets of current observations such as the cosmic microwave background from {\it Planck} and type Ia supernovae from Joint Light-curve Analysis (JLA). We find that the high redshift SEAMBHs can help to break the degeneracies of the background cosmological parameters constrained by {\it Planck} and JLA, thus giving much tighter constraints of the cosmological parameters. The second uses the high redshift SEAMBHs as the complements of the low redshift JLA to constrain the early expansion rate and the dark energy density evolution in the cold dark matter frame. Our results show that these high redshift SEAMBHs are very powerful on constraining the early Hubble rate and the evolution of the dark energy density; thus they can give us more information about the expansion history of our Universe, which is also crucial for testing the $\Lambda$CDM model in the high redshift region. Finally, we check the SEAMBH candles' abilities to reconstruct the equation of state of dark energy at high redshift. In summary, our results show that the SEAMBHs, as the rare candles in the high redshift region, can provide us a new and independent observation to probe cosmology in the future.

The gravitational wave (GW) as a standard siren directly determines the luminosity distance from the gravitational waveform without reference to the specific cosmological model, of which the redshift can be obtained separately by means of the electromagnetic counterpart like GW events from binary neutron stars and massive black hole binaries (MBHBs). To see to what extent the standard siren can reproduce the presumed dipole anisotropy written in the simulated data of standard siren events from typical configurations of GW detectors, we find that (1) for the Laser Interferometer Space Antenna with different MBHB models during five-year observations, the cosmic isotropy can be ruled out at $3\sigma$ confidence level (C.L.) and the dipole direction can be constrained roughly around $20\%$ at $2\sigma$ C.L., as long as the dipole amplitude is larger than $0.03$, $0.06$ and $0.025$ for MBHB models Q3d, pop III and Q3nod with increasing constraining ability, respectively; (2) for Einstein Telescope with no less than $200$ standard siren events, the cosmic isotropy can be ruled out at $3\sigma$ C.L. if the dipole amplitude is larger than $0.06$, and the dipole direction can be constrained within $20\%$ at $3\sigma$ C.L. if the dipole amplitude is near $0.1$; (3) for the Deci-Hertz Interferometer Gravitational wave Observatory with no less than $100$ standard siren events, the cosmic isotropy can be ruled out at $3\sigma$ C.L. for dipole amplitude larger than $0.03$ , and the dipole direction can even be constrained within $10\%$ at $3\sigma$ C.L. if the dipole amplitude is larger than $0.07$. Our work manifests the promising perspective of the constraint ability on the cosmic anisotropy from the standard siren approach.

The Kepler-9 system harbors three known transiting planets. The system holds significant interest for several reasons. First, the outer two planets exhibit a period ratio that is close to a 2:1 orbital commensurability, with attendant dynamical consequences. Second, both planets lie in the planetary mass "desert" that is generally associated with the rapid gas agglomeration phase of the core accretion process. Third, there exist attractive prospects for accurately measuring both the sky-projected stellar spin-orbit angles as well as the mutual orbital inclination between the planets in the system. Following the original \textit{Kepler} detection announcement in 2010, the initially reported orbital ephemerides for Kepler-9~b and c have degraded significantly, due to the limited time base-line of observations on which the discovery of the system rested. Here, we report new ground-based photometric observations and extensive dynamical modeling of the system. These efforts allow us to photometrically recover the transit of Kepler-9~b, and thereby greatly improve the predictions for upcoming transit mid-times. Accurate ephemerides of this system are important in order to confidently schedule follow-up observations of this system, for both in-transit Doppler measurements as well as for atmospheric transmission spectra taken during transit.

We present a study on the pulsations of the high-amplitude $\delta$ Scuti star YZ Boo based on photometric observations in Johnson $V$ and $R$ bands with both the Nanshan 1-m telescope of Xinjiang Astronomical Observatory (XAO) and the Xinglong 85-cm telescope of National Astronomical Observatories, Chinese Academy of Sciences (NAOC). The Fourier analysis of the light curves reveals the fundamental radial mode and its five harmonics, with the fourth and the fifth newly detected. 39 new times of light maximum are determined from the light curves, and combined with those in the literature, we construct the $O - C$ diagram and derive a new ephemeris and the determination of a new value of the updated period 0.104091579(2). Besides, the $O - C$ diagram reveals an increasing period change rate of YZ Boo. Theoretical models are calculated and constrained with the observationally determined parameters of YZ Boo. The mass and age of YZ Boo are hence derived as $M$ = 1.61 $\pm$ 0.05 $M_\odot$ and $age$ = (1.44 $\pm$ 0.14) $\times$ 10$^9$ yr, respectively. With both the frequency of the fundamental radial mode and the period change rate, YZ Boo is located at the post-main-sequence stage.

A recent surge in the discovery of new ultrafaint dwarf satellites of the Milky Way has inspired the idea of searching for faint satellites, $10^3\, \mathrm{M_{\odot}} < M_* < 10^6\, \mathrm{M_{\odot}}$, around less massive field galaxies in the Local Group. Such satellites would be subject to weaker environmental influences than Milky Way satellites, and could lead to new insights on low mass galaxy formation. In this paper, we predict the number of luminous satellites expected around field dwarf galaxies by applying several abundance matching models and a reionization model to the dark-matter only \textit{Caterpillar} simulation suite. For three of the four abundance matching models used, we find a $>99\%$ chance that at least one satellite with stellar mass $M_*> 10^5\, \mathrm{M_{\odot}}$ exists around the combined five Local Group field dwarf galaxies with the largest stellar mass. When considering satellites with $M_*> 10^4\, \mathrm{M_{\odot}}$, we predict a combined $5-25$ satellites for the five largest field dwarfs, and $10-50$ for the whole Local Group field dwarf population. Because of the relatively small number of predicted dwarfs, and their extended spatial distribution, a large fraction each Local Group dwarf's virial volume will need to be surveyed to guarantee discoveries. We compute the predicted number of satellites in a given field of view of specific Local Group galaxies, as a function of minimum satellite luminosity, and explicitly obtain such values for the Solitary Local dwarfs survey. Uncertainties in abundance matching and reionization models are large, implying that comprehensive searches could lead to refinements of both models.

The gravitational waves from compact binary systems are viewed as a standard siren to probe the evolution of the universe. This paper summarizes the potential and ability to use the gravitational waves to constrain the cosmological parameters and the dark sector interaction in the Gaussian process methodology. After briefly introducing the method to reconstruct the dark sector interaction by the Gaussian process, the concept of standard sirens and the analysis of reconstructing the dark sector interaction with LISA are outlined. Furthermore, we estimate the constraint ability of the gravitational waves on cosmological parameters with ET. The numerical methods we use are Gaussian process and the Markov-Chain Monte-Carlo. Finally, we also forecast the improvements of the abilities to constrain the cosmological parameters with ET and LISA combined with the \it Planck.

The direct detection of gravitational wave by Laser Interferometer Gravitational-Wave Observatory indicates the coming of the era of gravitational-wave astronomy and gravitational-wave cosmology. It is expected that more and more gravitational-wave events will be detected by currently existing and planned gravitational-wave detectors. The gravitational waves open a new window to explore the Universe and various mysteries will be disclosed through the gravitational-wave detection, combined with other cosmological probes. The gravitational-wave physics is not only related to gravitation theory, but also is closely tied to fundamental physics, cosmology and astrophysics. In this review article, three kinds of sources of gravitational waves and relevant physics will be discussed, namely gravitational waves produced during the inflation and preheating phases of the Universe, the gravitational waves produced during the first-order phase transition as the Universe cools down and the gravitational waves from the three phases: inspiral, merger and ringdown of a compact binary system, respectively. We will also discuss the gravitational waves as a standard siren to explore the evolution of the Universe.

We perform a forecast analysis of the ability of the LISA space-based interferometer to reconstruct the dark sector interaction using gravitational wave (GW) standard sirens at high redshift. We employ Gaussian process methods to reconstruct the distance-redshift relation in a model independent way. We adopt simulated catalogues of standard sirens given by merging massive black hole binaries (MBHBs) visible by LISA, with an electromagnetic (EM) counterpart detectable by future telescopes. The catalogues are constructed considering three different astrophysical scenarios for the evolution of MBHB mergers based on the semi-analytic model of E. Barausse (2012). We first use these standard siren datasets to assess the potential of LISA in reconstructing a possible interaction between vacuum dark energy and dark matter. Then we combine the LISA cosmological data with supernovae data simulated for the Dark Energy Survey (DES). We consider two scenarios distinguished by the time duration of the LISA mission: 5 and 10 years. Using only LISA standard siren data, the dark sector interaction can be well reconstructed from redshift $z\sim1$ to $z\sim3$ (5 yr) and $z\sim1$ up to $z\sim5$ (10 yr), though the reconstruction is inefficient at lower redshift. When combined with the DES datasets, the interaction is well reconstructed in the whole redshift region from $z\sim0$ to $z\sim3$ (5 yr) and $z\sim0$ to $z\sim5$ (10 yr). MBHB standard sirens can thus be used to constrain the dark sector interaction at redshift ranges not reachable by usual supernovae datasets which probe only the $z\lesssim 1.5$ range. GW standard sirens will not only constitute a complementary and alternative way, with respect to familiar EM observations, to probe the cosmic expansion, but will also provide new tests to constrain possible deviations from the standard $\Lambda$CDM dynamics, especially at high redshift.

We investigate the constraint ability of the gravitational wave (GW) as the standard siren on the cosmological parameters by using the third-generation gravitational wave detector: the Einstein Telescope. We simulate the luminosity distances and redshift measurements from 100 to 1000 GW events. We use two different algorithms to constrain the cosmological parameters. For the Hubble constant $H_0$ and dark matter density parameter $\Omega_m$, we adopt the Markov chain Monte Carlo approach. We find that with about 500-600 GW events we can constrain the Hubble constant with an accuracy comparable to \textit{Planck} temperature data and \textit{Planck} lensing combined results, while for the dark matter density, GWs alone seem not able to provide the constraints as good as for the Hubble constant; the sensitivity of 1000 GW events is a little lower than that of \textit{Planck} data. It should require more than 1000 events to match the \textit{Planck} sensitivity. Yet, for analyzing the more complex dynamical property of dark energy, i.e., the equation of state $w$, we adopt a new powerful nonparametric method: the Gaussian process. We can reconstruct $w$ directly from the observational luminosity distance at every redshift. In the low redshift region, we find that about 700 GW events can give the constraints of $w(z)$ comparable to the constraints of a constant $w$ by \textit{Planck} data with type Ia supernovae. Those results show that GWs as the standard sirens to probe the cosmological parameters can provide an independent and complementary alternative to current experiments.

The wide-field synoptic sky surveys, known as the Palomar Transient Factory (PTF) and the intermediate Palomar Transient Factory (iPTF), will accumulate a large number of known and new RR Lyrae. These RR Lyrae are good tracers to study the substructure of the Galactic halo if their distance, metallicity, and galactocentric velocity can be measured. Candidates of halo RR Lyrae can be identified from their distance and metallicity before requesting spectroscopic observations for confirmation. This is because both quantities can be obtained via their photometric light curves, because the absolute V-band magnitude for RR Lyrae is correlated with metallicity, and the metallicity can be estimated using a metallicity-light curve relation. To fully utilize the PTF and iPTF light-curve data in related future work, it is necessary to derive the metallicity-light curve relation in the native PTF/iPTF R-band photometric system. In this work, we derived such a relation using the known ab-type RR Lyrae located in the Kepler field, and it is found to be $[Fe/H]_{PTF} = -4.089 - 7.346 P + 1.280 \phi_{31}$ (where $P$ is pulsational period and $\phi_{31}$ is one of the Fourier parameters describing the shape of the light curve), with a dispersion of 0.118 dex. We tested our metallicity-light curve relation with new spectroscopic observations of a few RR Lyrae in the Kepler field, as well as several data sets available in the literature. Our tests demonstrated that the derived metallicity-light curve relation could be used to estimate metallicities for the majority of the RR Lyrae, which are in agreement with the published values.

We report our measurements for orbital and spin parameters of X 1822-371 using its X-ray partial eclipsing profile and pulsar timing from data collected by the Rossi X-ray Timing Explorer (RXTE). Four more X-ray eclipse times obtained by the RXTE 2011 observations were combined with historical records to trace evolution of orbital period. We found that a cubic ephemeris likely better describes evolution of the X-ray eclipse times during a time span of about 34 years with a marginal second order derivative of $\ddot{P}_{orb}=(-1.05 \pm 0.59) \times 10^{-19}$ s$^{-1}$. Using the pulse arrival time delay technique, the orbital and spin parameters were obtained from RXTE observations from 1998 to 2011. The detected pulse periods show that the neutron star in X 1822-371 is continuously spun-up with a rate of $\dot{P}_{s}=(-2.6288 \pm 0.0095) \times 10^{-12}$ s s$^{-1}$. Evolution of the epoch of the mean longitude $l=\pi /2$ (i.e. $T_{\pi / 2}$) gives an orbital period derivative value consistent with that obtained from the quadratic ephemeris evaluated by the X-ray eclipse but the detected $T_{\pi / 2}$ values are significantly and systematically earlier than the corresponding expected X-ray eclipse times by $90 \pm 11$ s. This deviation is probably caused by asymmetric X-ray emissions. We also attempted to constrain the mass and radius of the neutron star using the spin period change rate and concluded that the intrinsic luminosity of X 1822-371 is likely more than $10^{38}$ ergs s$^{-1}$.

We develop a new model-independent method to probe the constancy of the speed of light $c$. In our method, the degeneracy between the cosmic curvature and the speed of light can be eliminated, which makes the test more natural and general. Combining the independent observations of Hubble parameter $H(z)$ and luminosity distance $d_L(z)$, we use the model-independent smoothing technique, Gaussian processes, to reconstruct them and then detect variation of the speed of light. We find no signal of deviation from the present value of the speed of light $c_0$. Moreover, to demonstrate the improvement in probing the constancy of the speed of light from future experiments, we produce a series of simulated data. The Dark Energy Survey will be able to detect $\Delta c /c_0 \sim 1\%$ at $\sim 1.5\sigma$ confidence level and $\Delta c /c_0 \sim 2\%$ at $\sim 3\sigma$ confidence level. If the errors are reduced to one-tenth of the expected DES ones, it can detect a $\Delta c /c_0 \sim 0.1\%$ variation at $\sim 2\sigma$ confidence level.

We report the detection of orbital modulation, a model solution, and X-ray properties of a newly discovered contact binary, 2MASS J11201034$-$2201340. We serendipitously found this X-ray point source outside the error ellipse when searching for possible X-ray counterparts of $\gamma$-ray millisecond pulsars among the unidentified objects detected by the {\it Fermi Gamma-ray Space Telescope}. The optical counterpart of the X-ray source (unrelated to the $\gamma$-ray source) was then identified using archival databases. The long-term CRTS survey detected a precise signal with a period of $P=0.28876208(56)$ days. A follow-up observation made by the SLT telescope of Lulin Observatory revealed the binary nature of the object. Utilizing archived photometric data of multi-band surveys, we construct the spectral energy distribution, which is well fitted by a K2V spectral template. The fitting result of the orbital profile using the Wilson--Devinney code suggests that 2MASS J11201034-2201340 is a short-period A-type contact binary and the more massive component has a cool spot. The X-ray emission was first noted in observations made by \swift\, then further confirmed and characterized by an \xmm\ observation. The X-ray spectrum can be described by a power law or thermal Bremsstrahlung. Unfortunately, we could not observe significant X-ray orbital modulation. Finally, according to the spectral energy distribution, this system is estimated to be 690 pc from Earth with a calculated X-ray intensity of $(0.7-1.5)\times 10^{30}$ erg s$^{-1}$, which is in the expected range of an X-ray emitting contact binary.

We introduce a model-independent approach to the null test of the cosmic curvature which is geometrically related to the Hubble parameter $H(z)$ and luminosity distance $d_L(z)$. Combining the independent observations of $H(z)$ and $d_L(z)$, we use the model-independent smoothing technique, Gaussian processes, to reconstruct them and determine the cosmic curvature $\Omega_K^{(0)}$ in the null test relation. The null test is totally geometrical and without assuming any cosmological model. We show that the cosmic curvature $\Omega_K^{(0)}=0$ is consistent with current observational data sets, falling within the $1\sigma$ limit. To demonstrate the effect on the precision of the null test, we produce a series of simulated data of the models with different $\Omega_K^{(0)}$. Future observations in better quality can provide a greater improvement to constrain or refute the flat universe with $\Omega_K^{(0)}=0$.

Relic gravitational waves (RGWs) , a background originated during inflation, would give imprints on the pulsar timing residuals. This makes RGWs be one of important sources for detection using the method of pulsar timing. In this paper, we discuss the effects of RGWs on the single pulsar timing, and give quantitively the timing residuals caused by RGWs with different model parameters. In principle, if the RGWs are strong enough today, they can be detected by timing a single millisecond pulsar with high precision after the intrinsic red noise in pulsar timing residuals were understood, even though observing simultaneously multiple millisecond pulsars is a more powerful technique in extracting gravitational wave signals. We corrected the normalization of RGWs using observations of the cosmic microwave background (CMB), which leads to the amplitudes of RGWs being reduced by two orders of magnitude or so compared to our previous works. We made new constraints on RGWs using the recent observations from the Parkes Pulsar Timing Array, employing the tensor-to-scalar ratio $r=0.2$ due to the tensor-type polarization observations of CMB by BICEP2 as a referenced value even though it has been denied. Moreover, the constraints on RGWs from CMB and BBN (Big Bang nucleosynthesis) will also be discussed for comparison.

We present the time-frequency analysis results based on the Hilbert-Huang transform (HHT) for the evolution of a 4-Hz low-frequency quasi-periodic oscillation (LFQPO) around the black hole X-ray binary XTE J1550-564. The origin of LFQPOs is still debated. To understand the cause of the peak broadening, we utilized a recently developed time-frequency analysis, HHT, for tracking the evolution of the 4-Hz LFQPO from XTE J1550 564. By adaptively decomposing the ~4-Hz oscillatory component from the light curve and acquiring its instantaneous frequency, the Hilbert spectrum illustrates that the LFQPO is composed of a series of intermittent oscillations appearing occasionally between 3 Hz and 5 Hz. We further characterized this intermittency by computing the confidence limits of the instantaneous amplitudes of the intermittent oscillations, and constructed both the distributions of the QPO's high and low amplitude durations, which are the time intervals with and without significant ~4-Hz oscillations, respectively. The mean high amplitude duration is 1.45 s and 90% of the oscillation segments have lifetimes below 3.1 s. The mean low amplitude duration is 0.42 s and 90% of these segments are shorter than 0.73 s. In addition, these intermittent oscillations exhibit a correlation between the oscillation's rms amplitude and mean count rate. This correlation could be analogous to the linear rms-flux relation found in the 4-Hz LFQPO through Fourier analysis. We conclude that the LFQPO peak in the power spectrum is broadened owing to intermittent oscillations with varying frequencies, which could be explained by using the Lense-Thirring precession model.

This article investigates the problem of estimating stellar atmospheric parameters from spectra. Feature extraction is a key procedure in estimating stellar parameters automatically. We propose a scheme for spectral feature extraction and atmospheric parameter estimation using the following three procedures: firstly, learn a set of basic structure elements (BSE) from stellar spectra using an autoencoder; secondly, extract representative features from stellar spectra based on the learned BSEs through some procedures of convolution and pooling; thirdly, estimate stellar parameters ($T_{eff}$, log$~g$, [Fe/H]) using a back-propagation (BP) network. The proposed scheme has been evaluated on both real spectra from Sloan Digital Sky Survey (SDSS)/Sloan Extension for Galactic Understanding and Exploration (SEGUE) and synthetic spectra calculated from Kurucz's new opacity distribution function (NEWODF) models. The best mean absolute errors (MAEs) are 0.0060 dex for log$~T_{eff}$, 0.1978 dex for log$~g$ and 0.1770 dex for [Fe/H] for the real spectra and 0.0004 dex for log$~T_{eff}$, 0.0145 dex for log$~g$ and 0.0070 dex for [Fe/H] for the synthetic spectra.

Two dedicated asteroid rotation-period surveys have been carried out using data taken on January 6-9 and February 20-23 of 2014 by the Intermediate Palomar Transient Factory (iPTF) in the $R$~band with $\sim 20$-min cadence. The total survey area covered 174~deg$^2$ in the ecliptic plane. Reliable rotation periods for 1,438 asteroids are obtained from a larger data set of 6,551 mostly main-belt asteroids, each with $\geq 10$~detections. Analysis of 1751, PTF based, reliable rotation periods clearly shows the "spin barrier" at $\sim 2$~hours for "rubble-pile" asteroids. We also found a new large-sized super-fast rotator, 2005 UW163 (Chang et al., 2014), and other five candidates as well. Our spin-rate distributions of asteroids with $3 < D < 15$~km shows number decrease when frequency greater than 5 rev/day, which is consistent to that of the Asteroid Light Curve Database (LCDB, Warner et al., 2009) and the result of (Masiero et al., 2009). We found the discrepancy in the spin-rate distribution between our result and (Pravec et al., 2008, update 2014-04-20) is mainly from asteroids with $\Delta m < 0.2$ mag that might be primarily due to different survey strategies. For asteroids with $D \leq 3$~km, we found a significant number drop at $f = 6$ rev/day. The YORP effect timescale for small-sized asteroid is shorter that makes more elongate objets spun up to reach their spin-rate limit and results in break-up. The K-S test suggests a possible difference in the spin-rate distributions of C- and S-type asteroids. We also find that C-type asteroids have a smaller spin-rate limit than the S-type, which agrees with the general sense that the C-type has lower bulk density than the S-type.

We present a nonparametric approach to reconstruct the interaction between dark energy and dark matter directly from SNIa Union 2.1 data using Gaussian processes, which is a fully Bayesian approach for smoothing data. In this method, once the equation of state ($w$) of dark energy is specified, the interaction can be reconstructed as a function of redshift. For the decaying vacuum energy case with $w=-1$, the reconstructed interaction is consistent with the standard $\Lambda$CDM model, namely, there is no evidence for the interaction. This also holds for the constant $w$ cases from $-0.9$ to $-1.1$ and for the Chevallier-Polarski-Linder (CPL) parametrization case. If the equation of state deviates obviously from $-1$, the reconstructed interaction exists at $95\%$ confidence level. This shows the degeneracy between the interaction and the equation of state of dark energy when they get constraints from the observational data.

Large-scale and deep sky survey missions are rapidly collecting a large amount of stellar spectra, which necessitate the estimation of atmospheric parameters directly from spectra and makes it feasible to statistically investigate latent principles in a large dataset. We present a technique for estimating parameters $T_{eff}$, log$~g$ and [Fe/H] from stellar spectra. With this technique, we first extract features from stellar spectra using the LASSO algorithm; then, the parameters are estimated from the extracted features using the SVR. On a subsample of 20~000 stellar spectra from SDSS with reference parameters provided by SDSS/SEGUE Pipeline SSPP, estimation consistency are 0.007458 dex for log$~T_{eff}$ (101.609921 K for $T_{eff}$), 0.189557 dex for log$~g$ and 0.182060 for [Fe/H], where the consistency is evaluated by mean absolute error. Prominent characteristics of the proposed scheme are sparseness, locality, and physical interpretability. In this work, every spectrum consists of 3821 fluxes, and 10, 19, and 14 typical wavelength positions are detected respectively for estimating $T_{eff}$, log$~g$ and [Fe/H]. It is shown that the positions are related to typical lines of stellar spectra. This characteristic is important in investigating physical indications from analysis results. Then, stellar spectra can be described by the individual fluxes on the detected positions (PD) or local integration of fluxes near them (LI). The abovementioned consistency is the result based on features described by LI. If features are described by PD, consistency are 0.009092 dex for log$~T_{eff}$ (124.545075 K for $T_{eff}$), 0.198928 dex for log$~g$, and 0.206814 dex for [Fe/H].

A new asteroid rotation period survey have been carried out by using the Palomar Transient Factory (PTF). Twelve consecutive PTF fields, which covered an area of 87 deg$^2$ in the ecliptic plane, were observed in $R$ band with a cadence of $\sim$20 min during February 15--18, 2013. We detected 2500 known asteroids with a diameter range of 0.5 km $\leq D \leq$ 200 km. Of these, 313 objects had highly reliable rotation periods and exhibited the "spin barrier" at $\sim2$ hours. In contrast to the flat spin rate distribution of the asteroids with 3 km $\leq D \leq$ 15 km shown by Pravec et al. (2008), our results deviated somewhat from a Maxwellian distribution and showed a decrease at the spin rate greater than 5 rev/day. One super-fast-rotator candidate and two possible binary asteroids were also found in this work.

The pulsar timing residuals induced by gravitational waves from non-evolving single binary sources are affected by many parameters related to the relative positions of the pulsar and the gravitational wave sources. We will fully analyze the effects due to different parameters one by one. The standard deviations of the timing residuals will be calculated with a variable parameter fixing a set of other parameters. The orbits of the binary sources will be generally assumed to be elliptical. The influences of different eccentricities on the pulsar timing residuals will also studied in detail. We find that effects of the related parameters are quite different, and some of them present certain regularities.

RE J1034+396, a narrow-line Seyfert-1 active galactic nucleus (AGN), is the first example of AGNs that exhibited a nearly coherent quasi-periodic oscillation (QPO) for the data collected by XMM-Newton in 2007. The spectral behaviors and timing properties of the QPO have been studied since its discovery. We present an analysis of the QPO in RE J1034+396 based on the Hilbert-Huang transform (HHT). Comparing other time-frequency analysis methods, the Hilbert spectrum reveals the variation of the QPO period in great detail. Furthermore, the empirical mode decomposition provides band-pass filtered data that can be used in the O -- C and correlation analysis. We suggest that it is better to divide the evolution of the QPO in this observation into three epochs according to their different periodicities. In addition to the periodicities, the correlations between the QPO periods and corresponding mean count rates are also different in these three epochs. Further examining the phase lags in these epochs, we found no significant phase lags between the soft and hard X-ray bands, which is also confirmed in the QPO phase-resolved spectral analysis. Finally, we discuss the indications of current models including a spotted accretion disk, diskoseismology, and oscillation of shock, according to the observed time-frequency and spectral behaviors.

We present 855 cataclysmic variable candidates detected by the Catalina Real-time Transient Survey (CRTS) of which at least 137 have been spectroscopically confirmed and 705 are new discoveries. The sources were identified from the analysis of five years of data, and come from an area covering three quarters of the sky. We study the amplitude distribution of the dwarf novae CVs discovered by CRTS during outburst, and find that in quiescence they are typically two magnitudes fainter compared to the spectroscopic CV sample identified by SDSS. However, almost all CRTS CVs in the SDSS footprint have ugriz photometry. We analyse the spatial distribution of the CVs and find evidence that many of the systems lie at scale heights beyond those expected for a Galactic thin disc population. We compare the outburst rates of newly discovered CRTS CVs with the previously known CV population, and find no evidence for a difference between them. However, we find that significant evidence for a systematic difference in orbital period distribution. We discuss the CVs found below the orbital period minimum and argue that many more are yet to be identified among the full CRTS CV sample. We cross-match the CVs with archival X-ray catalogs and find that most of the systems are dwarf novae rather than magnetic CVs.

In the non-standard model of relic gravitational waves (RGWs) generated in the early universe, the theoretical spectrum of is mainly described by an amplitude $r$ and a spectral index $\beta$, the latter usually being determined by the slope of the inflation potential. Pulsar timing arrays (PTAs) data have imposed constraints on the amplitude of strain spectrum for a power-law form as a phenomenological model. Applying these constraints to a generic, theoretical spectrum with $r$ and $\beta$ as independent parameters, we convert the PTAs constraint into an upper bound on the index $\beta$, which turns out to be less stringent than those upper bounds from BBN, CMB, and LIGO/VIRGO, respectively. Moreover, it is found that PTAs constrain the non-standard RGWs more stringent than the standard RGWs. If the condition of the quantum normalization is imposed upon a theoretical spectrum of RGWs, $r$ and $\beta$ become related. With this condition, a minimum requirement of the horizon size during inflation is greater than the Planck length results in an upper bound on $\beta$, which is comparable in magnitude to that by PTAs. When both PTAs and the quantum normalization are applied to a theoretical spectrum of RGWs, constraints can be obtained for other cosmic processes of the early universe, such as the reheating, a process less understood observationally so far. The resulting constraint is consistent with the slow-roll, massive scalar inflation model. The future SKA will be able to constrain RGWs further and might even detect RGWs, rendering an important probe to the very early universe.

The pulsar timing residuals induced by gravitational waves from non-evolving single binary sources with general elliptical orbits will be analyzed. For different orbital eccentricities, the timing residuals present different properties. The standard deviations of the timing residuals induced by a fixed gravitational wave source will be calculated for different values of the eccentricity. We will also analyze the timing residuals of PSR J0437-4715 induced by one of the best known single gravitational wave sources, the supermassive black hole binary in the blazar OJ287.

The high-mass X-ray binary SMC X-1 is an eclipsing binary with an orbital period of 3.89 d. This system exhibits a superorbital modulation with a period varying between ~40 d and ~65 d. The instantaneous frequency and the corresponding phase of the superorbital modulation can be obtained by a recently developed time-frequency analysis technique, the Hilbert-Huang transform (HHT). We present a phase-resolved analysis of both the spectra and the orbital profiles with the superorbital phase derived from the HHT. The X-ray spectra observed by the Proportional Counter Array onboard the Rossi X-ray Timing Explorer are fitted well by a blackbody plus a Comptonized component. The plasma optical depth, which is a good indicator of the distribution of material along the line of sight, is significantly anti-correlated with the flux detected at 2.5-25 keV. However, the relationship between the plasma optical depth and the equivalent width of the iron line is not monotonic: there is no significant correlation for fluxes higher than ~35 mCrab but clear positive correlation when the intensity is lower than ~20 mCrab. This indicates that the iron line production is dominated by different regions of this binary system in different superorbital phases. To study the dependence of the orbital profile on the superorbital phase, we obtained the eclipse profiles by folding the All Sky Monitor light curve with the orbital period for different superorbital states. A dip feature, similar to the pre-eclipse dip in Her X-1, lying at orbital phase ~0.6-0.85, was discovered during the superorbital transition state. This indicates that the accretion disk has a bulge that absorbs considerable X-ray emission in the stream-disk interaction region. The dip width is anti-correlated with the flux, and this relation can be interpreted by the precessing tilted accretion disk scenario.

The high-mass X-ray binary (HMXB) SMC X-1 exhibits a superorbital modulation with a dramatically varying period ranging between ~40 d and ~60 d. This research studies the time-frequency properties of the superorbital modulation of SMC X-1 based on the observations made by the All-Sky Monitor (ASM) onboard the Rossi X-ray Timing Explorer (RXTE).We analyzed the entire ASM database collected since 1996. The Hilbert-Huang Transform (HHT), developed for non-stationary and nonlinear time series analysis, was adopted to derive the instantaneous superorbital frequency. The resultant Hilbert spectrum is consistent with the dynamic power spectrum while it shows more detailed information in both the time and frequency domains. The RXTE observations manifest that the superorbital modulation period was mostly betweenn ~50 d and ~65 d, whenas it changed to ~45 d around MJD 50,800 and MJD 54,000. Our analysis further indicates that the instantaneous frequency changed in a time scale of hundreds of days between ~MJD 51,500 and ~MJD 53,500. Based on the instantaneous phase defined by HHT, we folded the ASM light curve to derive a superorbital profile, from which an asymmetric feature and a low state with barely any X-ray emissions (lasting for ~0.3 cycles) were observed. We also calculated the correlation between the mean period and the amplitude of the superorbital modulation. The result is similar to the recently discovered relationship between the superorbital cycle length and the mean X-ray flux for Her X-1.

This study reports pulse variation analysis results for the forth discovered accretion-powered millisecond pulsar XTE J1807-294 during its 2003 outburst observed by {\it Rossi X-ray Timing Explorer}. The pulsation is significantly detected only in the first $\sim$90d out of $\sim$150d observations. The pulse phase variation is too complex to be described as an orbital motion plus a simple polynomial model. The precise orbital parameters with $P_{orb}=40.073601(8)$ min and ${\it a_x}\sin {\it i}=4.823(5)$ lt-ms were obtained after applying the trend removal to the daily observed 150s segments pulse phases folded with a constant spin frequency without Keplerian orbit included. The binary barycenter corrected pulse phases show smooth evolution and clear negative phase shifts coincident with the flares seen on the light curve and the enhancements of fractional pulse amplitude. The non-flare pulse phases for the first $\sim$60d data are well described as a fourth order polynomial implying that the neutron star was spun-up during the first $\sim$60d with a rate $\dot \nu=(1.7\pm0.3) \times 10^{-13}$ Hz/s at the beginning of the outburst. Significant soft phase lags up to $\sim$500 $\mu s$ ($\sim$10% cycle) between 2 to 20 keV were detected for the nonflare pulse phases. We conclude that the anomalous phase shifts are unlikely due to the accretion torque but could result from the ``hot spot'' moving on the surface of neutron star.

Evidence of TeV emission from GRB970417a has been previously reported using data from the Milagrito detector. Constraints on the TeV fluence and the energy spectrum are now derived using additional data from a scaler system that recorded the rate of signals from the Milagrito photomultipliers. This analysis shows that if emission from GRB970417a has been observed, it must contain photons with energies above 650 GeV. Some consequences of this observation are discussed.

Milagrito, a detector sensitive to very high energy gamma rays, monitored the northern sky from February 1997 through May 1998. With a large field of view and a high duty cycle, this instrument was well suited to perform a search for TeV gamma-ray bursts (GRBs). We report on a search made for TeV counterparts to GRBs observed by BATSE. BATSE detected 54 GRBs within the field of view of Milagrito during this period. An excess of events coincident in time and space with one of these bursts, GRB 970417a, was observed by Milagrito. The excess has a chance probability of $2.8 \times 10^{-5}$ of being a fluctuation of the background. The probability for observing an excess at least this large from any of the 54 bursts is $1.5 \times 10^{-3}$. No significant correlations were detected from the other bursts.

Milagrito, a detector sensitive to gamma-rays at TeV energies, monitored the northern sky during the period February 1997 through May 1998. With a large field of view and high duty cycle, this instrument was used to perform a search for TeV counterparts to gamma-ray bursts. Within the Milagrito field of view 54 gamma-ray bursts at keV energies were observed by the Burst And Transient Satellite Experiment (BATSE) aboard the Compton Gamma-Ray Observatory. This paper describes the results of a preliminary analysis to search for TeV emission correlated with BATSE detected bursts. Milagrito detected an excess of events coincident both spatially and temporally with GRB 970417a, with chance probability $2.8 \times 10^{-5}$ within the BATSE error radius. No other significant correlations were detected. Since 54 bursts were examined the chance probability of observing an excess with this significance in any of these bursts is $1.5 \times 10^{-3}$. The statistical aspects and physical implications of this result are discussed.

Previous works have shown that water Cherenkov detectors have superior sensitivity to those of scintillation counters as applied to detecting extensive air showers (EAS). This is in large part due to their much higher sensitivity to EAS photons which are more than five times more numerous than EAS electrons. Large area water Cherenkov detectors can be constructed relatively cheaply and operated reliably. A sparse detector array has been designed which uses these types of detectors to substantially increase the area over which the Milagro Gamma Ray Observatory collects EAS information. Improvements to the Milagro detector's performance characteristics and sensitivity derived from this array and preliminary results from a prototype array currently installed near the Milagro detector will be presented.

The Milagro detector is an air shower array which uses the water Cherenkov technique and is capable of continuously monitoring the sky at energies near 1 TeV. The detector consists of 20000 metric tons of pure water instrumented with 723 photo-multiplier tubes (PMTs). The PMTs are arranged in a two-layer structure on a lattice of 3 m spacing covering 5000 $m^2$ area. The direction of the shower is determined from the relative timing of the PMT signals, necessitating a common time reference and amplitude slewing corrections to improve the time resolution. The calibration system to provide these consists of a pulsed laser driving 30 diffusing light sources deployed in the pond to allow cross-calibration of the PMTs. The system is capable of calibrating times and the pulse-heights from the PMTs using the time-over-threshold technique. The absolute energy scale is provided using single muons passing through the detector. The description of the calibration system of the Milagro detector and its prototype Milagrito will be presented.

The Milagro Gamma Ray Observatory is the world's first large-area water Cherenkov detector capable of continuously monitoring the sky at TeV energies. Located in northern New Mexico, Milagro will perform an all sky survey of the Northern Hemisphere at energies between ~250 GeV and 50 TeV. With a high duty cycle, large detector area (~5000 square meters), and a wide field-of-view (~1 sr), Milagro is uniquely capable of searching for transient and DC sources of high-energy gamma-ray emission. Milagro has been operating since February, 1999. The current status of the Milagro Observatory and initial results will be discussed.

Observing gamma ray bursts (GRBs) in the TeV energy range can be extremely valuable in providing insight to GRB radiation mechanisms and in constraining source distances. The Milagrito detector was an air shower array which used the water Cherenkov technique to search for TeV sources. Data from this detector was analyzed to look for a TeV component of GRBs coincident with low energy -rays detected by the BATSE instrument on the Compton Gamma Ray Observatory. A sample of 54 BATSE GRBs which were in the field of view of the Milagrito detector during its lifetime (February 1997 to May 1998) was used.

The Milagro water Cherenkov detector near Los Alamos, New Mexico, has been operated as a sky monitor at energies of a few TeV between February 1997 and April 1998. Serving as a test run for the full Milagro detector, Milagrito has taken data during the strong and long-lasting 1997 flare of Mrk 501. We present results from the analysis of Mrk 501 and compare the excess and background rates with expectations from the detector simulations.

The Milagrito water Cherenkov telescope operated for over a year. The most probable gamma-ray energy was ~1 TeV and the trigger rate was as high as 400 Hz. We have developed an efficient technique for searching the entire sky for short duration bursts of TeV photons. Such bursts may result from "traditional" gamma-ray bursts that were not in the field-of-view of any other instruments, the evaporation of primordial black holes, or some as yet undiscovered phenomenon. We have begun to search the Milagrito data set for bursts of duration 10 seconds. Here we will present the technique and the expected results. Final results will be presented at the conference.

Milagrito, a prototype for the Milagro detector, operated for 15 months in 1997-8 and collected 8.9 billion events. It was the first extensive air shower (EAS) array sensitive to showers intiated by primaries with energy below 1 TeV. The shadows of the sun and moon observed with cosmic rays can be used to study systematic pointing shifts and measure the angular resolution of EAS arrays. Below a few TeV, the paths of cosmic rays coming toward the earth are bent by the helio- and geo-magnetic fields. This is expected to distort and displace the shadows of the sun and the moon. The moon shadow, offset from the nominal (undeflected) position, has been observed with high statistical significance in Milagrito. This can be used to establish energy calibrations, as well as to search for the anti-matter content of the VHE cosmic ray flux. The shadow of the sun has also been observed with high significance.

Solar Energetic Particles from the 6 November 1997 solar flare/CME(coronal mass ejection) with energies exceeding 10 GeV have been detected by Milagrito, a prototype of the Milagro Gamma Ray Observatory. While particle acceleration beyond 1 GeV at the Sun is well established, few data exist for protons or ions beyond 10 GeV. The Milagro observatory, a ground based water Cherenkov detector designed for observing very high energy gamma ray sources, can also be used to study the Sun. Milagrito, which operated for approximately one year in 1997/98, was sensitive to solar proton and neutron fluxes above ~5- 10 GeV. Milagrito operated in a scaler mode, which was primarily sensitive to muons, low energy photons, and electrons, and the detector operated in a mode sensitive to showers and high zenith angle muons. In its scaler mode, Milagrito registered a rate increase coincident with the 6 November 1997 ground level event observed by Climax and other neutron monitors. A preliminary analysis suggests the presence of >10 GeV particles.

The Milagrito water Cherenkov detector in the Jemez Mountains near Los Alamos, New Mexico took data from February 1997 to April 1998. Milagrito served as a prototype for the larger Milagro detector, which has just begun operations. Milagrito was the first large-aperture gamma-ray detector with sensitivity to gamma rays below 1 TeV. We report here on a search for steady emission from point sources over most of the northern sky using data from Milagrito.