We present SPISEA (Stellar Population Interface for Stellar Evolution and Atmospheres), an open-source Python package that simulates simple stellar populations. The strength of SPISEA is its modular interface which offers the user control of 13 input properties including (but not limited to) the Initial Mass Function, stellar multiplicity, extinction law, and the metallicity-dependent stellar evolution and atmosphere model grids used. The user also has control over the Initial-Final Mass Relation in order to produce compact stellar remnants (black holes, neutron stars, and white dwarfs). We demonstrate several outputs produced by the code, including color-magnitude diagrams, HR-diagrams, luminosity functions, and mass functions. SPISEA is object-oriented and extensible, and we welcome contributions from the community. The code and documentation are available on GitHub and ReadtheDocs, respectively.

Providing a possible connection between neutrino emission and gravitational-wave (GW) bursts is important to our understanding of the physical processes that occur when black holes or neutron stars merge. In the Daya Bay experiment, using data collected from December 2011 to August 2017, a search has been performed for electron-antineutrino signals coinciding with detected GW events, including GW150914, GW151012, GW151226, GW170104, GW170608, GW170814, and GW170817. We used three time windows of $\mathrm{\pm 10~s}$, $\mathrm{\pm 500~s}$, and $\mathrm{\pm 1000~s}$ relative to the occurrence of the GW events, and a neutrino energy range of 1.8 to 100 MeV to search for correlated neutrino candidates. The detected electron-antineutrino candidates are consistent with the expected background rates for all the three time windows. Assuming monochromatic spectra, we found upper limits (90% confidence level) on electron-antineutrino fluence of $(1.13~-~2.44) \times 10^{11}~\rm{cm^{-2}}$ at 5 MeV to $8.0 \times 10^{7}~\rm{cm^{-2}}$ at 100 MeV for the three time windows. Under the assumption of a Fermi-Dirac spectrum, the upper limits were found to be $(5.4~-~7.0)\times 10^{9}~\rm{cm^{-2}}$ for the three time windows.

Based on the dynamics of single scalar field slow-roll inflation and the theory of reheating, we investigate the generalized natural inflationary (GNI) model. Concretely, we give constraints on the scalar spectral index $n_{s}$ and tensor-to scalar ratio $r$ for $\Lambda$CDM $+r$ model according to the latest data from Plack 2018 TT,TE,EE+lowE+lensing (P18) and BICEP2/Keck 2015 season (BK15), i.e., $n_{s}=0.9659\pm0.0044$ at $68\%$ confidence level (CL) and $r<0.0623$ at $95\%$CL. We find that the GNI model is favored by P18 plus BK15 in the ranges of $\log_{10}(f/M_{p})=0.62^{+0.17}_{-0.18}$ and $m=0.35^{+0.13}_{-0.23}$ at $68\%$CL. In addition, the corresponding predictions of the general and two-phase reheating are respectively discussed. It follows that the parameter $m$ has the significant effect on the model behaviors.

It is conjectured that coherent re-emission of cyclotron resonance absorption could result in pulsar giant pulses. This conjecture seems reasonable as it can naturally explain the distribution of pulsars with giant pulses on the $P$-$\dot{P}$ diagram.

Microlensing surveys have discovered thousands of events with almost all events discovered within the Galactic bulge or toward the Magellanic clouds. The Zwicky Transient Facility (ZTF), while not designed to be a microlensing campaign, is an optical time-domain survey that observes the entire northern sky every few nights including the Galactic plane. ZTF observes $\sim10^9$ stars in g-band and r-band and can significantly contribute to the observed microlensing population. We predict that ZTF will observe $\sim$1100 microlensing events in three years of observing within $10^\circ$ degrees latitude of the Galactic plane, with $\sim$500 events in the outer Galaxy ($\ell \geq 10^\circ$). This yield increases to $\sim$1400 ($\sim$800) events by combining every three ZTF exposures, $\sim$1800 ($\sim$900) events if ZTF observes for a total of five years, and $\sim$2400 ($\sim$1300) events for a five year survey with post-processing image stacking. Using the microlensing modeling software PopSyCLE, we compare the microlensing populations in the Galactic bulge and the outer Galaxy. We also present an analysis of the microlensing event ZTF18abhxjmj to demonstrate how to leverage these population statistics in event modeling. ZTF will constrain Galactic structure, stellar populations, and primordial black holes through photometric microlensing.

We present a new method to photometrically delineate between various sub-types of type Ia supernovae (SNe Ia). Using the color-stretch parameters, $s_{BV}$ or $s_{gr}$, and the time of i-band primary maximum relative to the B-band or g-band maximum it is demonstrated that 2003fg-like, 1991bg-like, and 2002cx-like SNe Ia can readily be identified. In the cases of these extreme SNe Ia, their primary i-band maximum occurs after the time of the B or g band maxima. We suggest that the timing of the i-band maximum can reveal the physical state of the SN Ia explosion as it traces: i) the speed of the recombination front of iron group elements in the ejecta, ii) the temperature evolution and rate of adiabatic cooling in the ejecta and, iii) the presence of interaction with a stellar envelope. This photometric sub-typing can be used in conjunction with other SNe analysis, such as the Branch diagram, to examine the physics and diversity of SNe Ia. The results here can also be used to screen out non-Ia SNe from cosmological samples that do not have complete spectroscopic typing. Finally, as future surveys like LSST create large databases of light curves of many objects this photometric identification can be used to readily identify and study the rates and bulk properties of peculiar SNe Ia.

Searching for space-time variations of the constants of Nature is a promising way to search for new physics beyond General Relativity and the standard model motivated by unification theories and models of dark matter and dark energy. We propose a new way to search for a variation of the fine-structure constant using measurements of late-type evolved giant stars from the S-star cluster orbiting the supermassive black hole in our Galactic Center. A measurement of the difference between distinct absorption lines (with different sensitivity to the fine structure constant) from a star leads to a direct estimate of a variation of the fine structure constant between the star's location and Earth. Using spectroscopic measurements of 5 stars, we obtain a constraint on the relative variation of the fine structure constant below $10^{-5}$. This is the first time a varying constant of Nature is searched for around a black hole and in a high gravitational potential. This analysis shows new ways the monitoring of stars in the Galactic Center can be used to probe fundamental physics.

We report the discovery of an eclipsing binary millisecond pulsar in the globular cluster M92 (NGC6341) with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). PSR J1717+4308A, or M92A, has a pulse frequency of 316.5~Hz (3.16~ms) and a dispersion measure of 35.45 pc cm$^{-3}$. The pulsar is a member of a binary system with an orbital period of 0.20~days around a low-mass companion which has a median mass of $\sim$0.18~\Ms. From observations so far, at least two eclipsing events have been observed in each orbit. The longer one lasted for ~5000~s in the orbital phase range 0.1--0.5. The other lasted for ~500~s and occurred between 1000--2000~s before or after the longer eclipsing event. The lengths of these two eclipsing events also change. These properties suggest that J1717+4308A is a ``red-back'' system with a low-mass main sequence or sub-giant companion. Timing observations of the pulsar and further searches of the data for additional pulsars are ongoing.

We report the discovery of a binary millisecond pulsar (namely PSR J1641+3627F or M13F) in the globular cluster M13 (NGC 6205) and timing solutions of M13A to F using observations made with the Five-hundred-metre Aperture Spherical radio Telescope (FAST). PSR J1641+3627F has a spin period of 3.00 ms and an orbital period of 1.4 days. The most likely companion mass is 0.16 M$_{\odot}$. M13A to E all have short spin periods and small period derivatives. We also confirm that the binary millisecond pulsar PSR J1641$+$3627E (also M13E) is a black widow with a companion mass around 0.02 M$_{\odot}$. We find that all the binary systems have low eccentricities compared to those typical for globular cluster pulsars and that they decrease with distance from the cluster core. This is consistent with what is expected as this cluster has a very low encounter rate per binary.

We present adaptive optics imaging from the NIRC2 instrument on the Keck-2 telescope that resolves the exoplanet host (and lens) star as it separates from the brighter source star. These observations yield the $K$-band brightness of the lens and planetary host star, as well as the lens-source relative proper motion, $\mu_{\rm rel,H}$. in the heliocentric reference frame. The $\mu_{\rm rel,H}$ measurement allows determination of the microlensing parallax vector, $\pi_E$, which had only a single component determined by the microlensing light curve. The combined measurements of $\mu_{\rm rel,H}$ and $K_L$ provide the masses of the host stat, $M_{\rm host} = 0.426\pm 0.037 M_\odot$, and planet, $m_p = 3.27 \pm 0.32 M_{\rm Jup}$ with a projected separation of $3.4\pm 0.5\,$AU. This confirms the tentative conclusion of a previous paper (Dong et al. 2009) that this super-Jupiter mass planet, OGLE-2005-BLG-071Lb, orbits an M-dwarf. Such planets are predicted to be rare by the core accretion theory and have been difficult to find with other methods, but there are two such planets with firm mass measurements from microlensing, and an additional 11 planetary microlens events with host mass estimates $< 0.5M_\odot$ and planet mass estimates $> 2$ Jupiter masses that could be confirmed by high angular follow-up observations. We also point out that OGLE-2005-BLG-071L has separated far enough from its host star that it should be possible to measure the host star metallicity withspectra from a high angular resolution telescope such as Keck, the VLT, the Hubble Space Telescope or the James Webb Space Telescope.

We apply the first law of thermodynamics to the apparent horizon of the universe with the power-law corrected and non-extensive Tsallis entropies rather than the Bekenstein-Hawking one. We examine the cosmological properties in the two entropy models by using the CosmoMC package. In particular, the first numerical study for the cosmological observables with the power-law corrected entropy is performed. We also show that the neutrino mass sum has a non-zero central value with a relaxed upper bound in the Tsallis entropy model comparing with that in the $\Lambda$CDM one.

We present a new Milky Way microlensing simulation code, dubbed PopSyCLE (Population Synthesis for Compact object Lensing Events). PopSyCLE is the first resolved microlensing simulation to include a compact object distribution derived from numerical supernovae explosion models and both astrometric and photometric microlensing effects. We demonstrate the capabilities of PopSyCLE by investigating the optimal way to find black holes (BHs) with microlensing. Candidate BHs have typically been selected from wide-field photometric microlensing surveys, such as OGLE, by selecting events with long Einstein crossing times ($t_E>120$ days). These events can be selected at closest approach and monitored astrometrically in order to constrain the mass of each lens; PopSyCLE predicts a BH detection rate of ~40% for such a program. We find that the detection rate can be enhanced to ~85% by selecting events with both $t_E>120$ days and a microlensing parallax of $\pi_E<0.08$. Unfortunately, such a selection criterion cannot be applied during the event as $\pi_E$ requires both pre- and post-peak photometry. However, historical microlensing events from photometric surveys can be revisited using this new selection criteria in order to statistically constrain the abundance of BHs in the Milky Way. The future WFIRST microlensing survey provides both precise photometry and astrometry and will yield individual masses of $\mathcal{O}(100-1000) black holes, which is at least an order of magnitude more than is possible with individual candidate follow-up with current facilities. The resulting sample of BH masses from WFIRST will begin to constrain the shape of the black hole present-day mass function, BH multiplicity, and BH kick velocity distributions.

The Wide-Field InfraRed Space Telescope (WFIRST) will be capable of delivering precise astrometry for faint sources over the enormous field of view of its main camera, the Wide-Field Imager (WFI). This unprecedented combination will be transformative for the many scientific questions that require precise positions, distances, and velocities of stars. We describe the expectations for the astrometric precision of the WFIRST WFI in different scenarios, illustrate how a broad range of science cases will see significant advances with such data, and identify aspects of WFIRST's design where small adjustments could greatly improve its power as an astrometric instrument.

Gaia16aye was a binary microlensing event discovered in the direction towards the northern Galactic disc and was one of the first microlensing events detected and alerted to by the Gaia space mission. Its light curve exhibited five distinct brightening episodes, reaching up to I=12 mag, and it was covered in great detail with almost 25,000 data points gathered by a network of telescopes. We present the photometric and spectroscopic follow-up covering 500 days of the event evolution. We employed a full Keplerian binary orbit microlensing model combined with the motion of Earth and Gaia around the Sun to reproduce the complex light curve. The photometric data allowed us to solve the microlensing event entirely and to derive the complete and unique set of orbital parameters of the binary lensing system. We also report on the detection of the first-ever microlensing space-parallax between the Earth and Gaia located at L2. The properties of the binary system were derived from microlensing parameters, and we found that the system is composed of two main-sequence stars with masses 0.57$\pm$0.05 $M_\odot$ and 0.36$\pm$0.03 $M_\odot$ at 780 pc, with an orbital period of 2.88 years and an eccentricity of 0.30. We also predict the astrometric microlensing signal for this binary lens as it will be seen by Gaia as well as the radial velocity curve for the binary system. Events such as Gaia16aye indicate the potential for the microlensing method of probing the mass function of dark objects, including black holes, in directions other than that of the Galactic bulge. This case also emphasises the importance of long-term time-domain coordinated observations that can be made with a network of heterogeneous telescopes.

In a spatially flat Friedmann-Lema\^{i}tre-Robertson-Walker universe, we investigate a unified cosmic fluid scenario endowed with bulk viscosity in which the coefficient of the bulk viscosity has a power law evolution. The power law in the bulk viscous coefficient is a general case in this study which naturally includes several choices as special cases. Considering such a general bulk viscous scenario, in the present work we have extracted the observational constraints using the latest cosmological datasets and examine their behaviour at the level of both background and perturbations. From the observational analyses, we find that a non-zero bulk viscous coefficient is always favored and some of the models in this series are able to weaken the current tension on $H_0$ for some dataset. However, from the Bayesian evidence analysis, $\Lambda$CDM is favored over the bulk viscous model.

We present $81$ near-infrared (NIR) spectra of $30$ Type II supernovae (SNe II) from the Carnegie Supernova Project-II (CSP-II), the largest such dataset published to date. We identify a number of NIR features and characterize their evolution over time. The NIR spectroscopic properties of SNe II fall into two distinct groups. This classification is first based on the strength of the He I $\lambda1.083\,\mu$m absorption during the plateau phase; SNe II are either significantly above (spectroscopically strong) or below $50$ angstroms (spectroscopically weak) in pseudo equivalent width. However between the two groups, other properties, such as the timing of CO formation and the presence of Sr II, are also observed. Most surprisingly, the distinct weak and strong NIR spectroscopic classes correspond to SNe II with slow and fast declining light curves, respectively. These two photometric groups match the modern nomenclature of SNe IIP and IIL. Including NIR spectra previously published, 18 out of 19 SNe II follow this slow declining-spectroscopically weak and fast declining-spectroscopically strong correspondence. This is in apparent contradiction to the recent findings in the optical that slow and fast decliners show a continuous distribution of properties. The weak SNe II show a high-velocity component of helium that may be caused by a thermal excitation from a reverse-shock created by the outer ejecta interacting with the red supergiant wind, but the origin of the observed dichotomy is not understood. Further studies are crucial in determining whether the apparent differences in the NIR are due to distinct physical processes or a gap in the current data set.

The authors discuss potential remarkable achievements for pulsar science with the FAST (pulsar monitoring, timing and searching, as well as others related), and expect a FAST era of pulsar science to come.

We report new infrared measurements of the supermassive black hole at the Galactic Center, Sgr A*, over a decade that was previously inaccessible at these wavelengths. This enables a variability study that addresses variability timescales that are ten times longer than earlier published studies. Sgr A* was initially detected in the near-infrared with adaptive optics observations in 2002. While earlier data exists in form of speckle imaging (1995 - 2005), Sgr A* was not detected in the initial analysis. Here, we improved our speckle holography analysis techniques. This has improved the sensitivity of the resulting speckle images by up to a factor of three. Sgr A* is now detectable in the majority of epochs covering 7 years. The brightness of Sgr A* in the speckle data has an average observed K magnitude of 16.0, which corresponds to a dereddened flux density of $3.4$ mJy. Furthermore, the flat power spectral density (PSD) of Sgr A* between $\sim$80 days and 7 years shows its uncorrelation in time beyond the proposed single power-law break of $\sim$245 minutes. We report that the brightness and its variability is consistent over 22 years. This analysis is based on simulations using Witzel et al. (2018) model to characterize infrared variability from 2006 to 2016. Finally, we note that the 2001 periapse of the extended, dusty object G1 had no apparent effect on the near-infrared emission from accretion flow onto Sgr A*. The result is consistent with G1 being a self-gravitating object rather than a disrupting gas cloud.

Adaptive optics (AO) systems deliver high-resolution images that may be ideal for precisely measuring positions of stars (i.e. astrometry) if the system has stable and well-calibrated geometric optical distortions. A calibration unit, equipped with back-illuminated pinhole mask, can be utilized to measure instrumental optical distortions. AO systems on the largest ground-based telescopes, such as the W. M. Keck Observatory and the Thirty Meter Telescope require pinhole positions known to 20 nm to achieve an astrometric precision of 0.001 of a resolution element. We characterize a photo-lithographic pinhole mask and explore the systematic errors that result from different experimental setups. We characterized the nonlinear geometric distortion of a simple imaging system using the mask; and we measured 857 nm RMS of optical distortion with a final residual of 39 nm (equivalent to 20 {\mu}as for TMT). We use a sixth order bivariate Legendre polynomial to model the optical distortion and allow the reference positions of the individual pinholes to vary. The nonlinear deviations in the pinhole pattern with respect to the manufacturing design of a square pattern are 47.2 nm +/- 4.5 nm (random) +/- 10.8 nm (systematic) over an area of 1788 mm$^2$. These deviations reflect the additional error induced when assuming the pinhole mask is manufactured perfectly square. We also find that ordered mask distortions are significantly more difficult to characterize than random mask distortions as the ordered distortions can alias into optical camera distortion. Future design simulations for astrometric calibration units should include ordered mask distortions. We conclude that photo-lithographic pinhole masks are >10 times better than the pinhole masks deployed in first generation AO systems and are sufficient to meet the distortion calibration requirements for the upcoming thirty meter class telescopes.

WFIRST is NASA's first flagship mission with pre-defined core science programs to study dark energy and perform a statistical census of wide orbit exoplanets with a gravitational microlensing survey. Together, these programs are expected to use more than half of the prime mission observing time. Previously, only smaller, PI-led missions have had core programs that used such a large fraction of the observing time, and in many cases, the data from these PI-led missions was reserved for the PI's science team for a proprietary period that allowed the PI's team to make most of the major discoveries from the data. Such a procedure is not appropriate for a flagship mission, which should provide science opportunities to the entire astronomy community. For this reason, there will be no proprietary period for WFIRST data, but we argue that a larger effort to make WFIRST science accessible to the astronomy community is needed. We propose a plan to enhance community involvement in the WFIRST exoplanet microlensing survey in two different ways. First, we propose a set of high level data products that will enable astronomers without detailed microlensing expertise access to the statistical implications of the WFIRST exoplanet microlensing survey data. And second, we propose the formation of a WFIRST Exoplanet Microlensing Community Science Team that will open up participation in the development of the WFIRST exoplanet microlensing survey to the general astronomy community in collaboration for the NASA selected science team, which will have the responsibility to provide most of the high level data products. This community science team will be open to volunteers, but members should also have the opportunity to apply for funding.

General Relativity predicts that a star passing close to a supermassive black hole should exhibit a relativistic redshift. We test this using observations of the Galactic center star S0-2. We combine existing spectroscopic and astrometric measurements from 1995-2017, which cover S0-2's 16-year orbit, with measurements in 2018 March to September which cover three events during its closest approach to the black hole. We detect the combination of special relativistic- and gravitational-redshift, quantified using a redshift parameter, $\Upsilon$. Our result, $\Upsilon=0.88 \pm 0.17$, is consistent with General Relativity ($\Upsilon=1$) and excludes a Newtonian model ($\Upsilon=0$ ) with a statistical significance of 5 $\sigma$.

We study the properties of hot beta-stable nuclear matter using equations of state derived within the Brueckner-Hartree-Fock approach at finite temperature including consistent three-body forces. Simple and accurate parametrizations of the finite-temperature equations of state are provided. The properties of hot neutron stars are then investigated within this framework, in particular the temperature dependence of the maximum mass. We find very small temperature effects and analyze the interplay of the different contributions.

We propose a new approach to Bayesian prior probability distributions (priors) that can improve orbital solutions for low-phase-coverage orbits, where data cover less than approximately 40% of an orbit. In instances of low phase coverage such as with stellar orbits in the Galactic center or with directly-imaged exoplanets, data have low constraining power and thus priors can bias parameter estimates and produce under-estimated confidence intervals. Uniform priors, which are commonly assumed in orbit fitting, are notorious for this. We propose a new observable-based prior paradigm that is based on uniformity in observables. We compare performance of this observable-based prior and of commonly assumed uniform priors using Galactic center and directly-imaged exoplanet (HR 8799) data. The observable-based prior can reduce biases in model parameters by a factor of two and helps avoid under-estimation of confidence intervals for simulations with less than about 40% phase coverage. Above this threshold, orbital solutions for objects with sufficient phase coverage such as S0-2, a short-period star at the Galactic center with full phase coverage, are consistent with previously published results. Below this threshold, the observable-based prior limits prior influence in regions of prior dominance and increases data influence. Using the observable-based prior, HR 8799 orbital analyses favor lower eccentricity orbits and provide stronger evidence that the four planets have a consistent inclination around 30 degrees to within 1-sigma. This analysis also allows for the possibility of coplanarity. We present metrics to quantify improvements in orbital estimates with different priors so that observable-based prior frameworks can be tested and implemented for other low-phase-coverage orbits.

In various fields of physics and astronomy, access to experimental facilities or to telescopes is becoming more and more competitive and limited. It becomes therefore important to optimize the type of measurements and their scheduling to reach a given scientific objective and to increase the chances of success of a scientific project. In this communication, extending the work of Ford (2008) and of Loredo et al. (2012), we present an efficient adaptive scheduling tool aimed at prioritzing measurements in order to reach a scientific goal. The algorithm, based on the Fisher matrix, can be applied to a wide class of measurements. We present this algorithm in detail and discuss some practicalities such as systematic errors or measurements losses due to contigencies (such as weather, experimental failure, ...). As an illustration, we consider measurements of the short-period star S0-2 in our Galactic Center. We show that the radial velocity measurements at the two turning points of the radial velocity curve are more powerful for detecting the gravitational redshift than measurements at the maximal relativistic signal. We also explicitly present the methodology that was used to plan measurements in order to detect the relativistic redshift considering systematics and possible measurements losses. For the future, we identify the astrometric turning points to be highly sensitive to the relativistic advance of the periastron. Finally, we also identify measurements particularly sensitive to the distance to our Galactic Center: the radial velocities around periastron and the astrometric measurements just before closest approach and at the maximal right ascension astrometric turning point.

Over the past three decades instruments on the ground and in space have discovered thousands of planets outside the solar system. These observations have given rise to an astonishingly detailed picture of the demographics of short-period planets, but are incomplete at longer periods where both the sensitivity of transit surveys and radial velocity signals plummet. Even more glaring is that the spectra of planets discovered with these indirect methods are either inaccessible (radial velocity detections) or only available for a small subclass of transiting planets with thick, clear atmospheres. Direct detection can be used to discover and characterize the atmospheres of planets at intermediate and wide separations, including non-transiting exoplanets. Today, a small number of exoplanets have been directly imaged, but they represent only a rare class of young, self-luminous super-Jovian-mass objects orbiting tens to hundreds of AU from their host stars. Atmospheric characterization of planets in the <5 AU regime, where radial velocity (RV) surveys have revealed an abundance of other worlds, is technically feasible with 30-m class apertures in combination with an advanced AO system, coronagraph, and suite of spectrometers and imagers. There is a vast range of unexplored science accessible through astrometry, photometry, and spectroscopy of rocky planets, ice giants, and gas giants. In this whitepaper we will focus on one of the most ambitious science goals --- detecting for the first time habitable-zone rocky (<1.6 R_Earth) exoplanets in reflected light around nearby M-dwarfs

Astrophysical and cosmological observations currently provide the only robust, empirical measurements of dark matter. Future observations with Large Synoptic Survey Telescope (LSST) will provide necessary guidance for the experimental dark matter program. This white paper represents a community effort to summarize the science case for studying the fundamental physics of dark matter with LSST. We discuss how LSST will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational couplings to the Standard Model, and compact object abundances. Additionally, we discuss the ways that LSST will complement other experiments to strengthen our understanding of the fundamental characteristics of dark matter. More information on the LSST dark matter effort can be found at https://lsstdarkmatter.github.io/ .

The Quintuplet star cluster is one of only three known young ($<10$ Myr) massive (M $>10^4$ M$_\odot$) clusters within $\sim100$ pc of the Galactic Center. In order to explore star cluster formation and evolution in this extreme environment, we analyze the Quintuplet's dynamical structure. Using the HST WFC3-IR instrument, we take astrometric and photometric observations of the Quintuplet covering a $120''\times120''$ field-of-view, which is $19$ times larger than those of previous proper motion studies of the Quintuplet. We generate a catalog of the Quintuplet region with multi-band, near-infrared photometry, proper motions, and cluster membership probabilities for $10,543$ stars. We present the radial density profile of $715$ candidate Quintuplet cluster members with $M\gtrsim4.7$ M$_\odot$ out to $3.2$ pc from the cluster center. A $3\sigma$ lower limit of $3$ pc is placed on the tidal radius, indicating the lack of a tidal truncation within this radius range. Only weak evidence for mass segregation is found, in contrast to the strong mass segregation found in the Arches cluster, a second and slightly younger massive cluster near the Galactic Center. It is possible that tidal stripping hampers a mass segregation signature, though we find no evidence of spatial asymmetry. Assuming that the Arches and Quintuplet formed with comparable extent, our measurement of the Quintuplet's comparatively large core radius of $0.62^{+0.10}_{-0.10}$ pc provides strong empirical evidence that young massive clusters in the Galactic Center dissolve on a several Myr timescale.

Gravitational microlensing is one of the few means of finding primordial black holes (PBHs), if they exist. Recent LIGO detections of 30 Msun black holes have re-invigorated the search for PBHs in the 10-100 Msun mass regime. Unfortunately, individual PBH microlensing events cannot easily be distinguished from stellar lensing events from photometry alone. However, the distribution of microlensing timescales (tE, the Einstein radius crossing time) can be analyzed in a statistical sense using models of the Milky Way with and without PBHs. While previous works have presented both theoretical models and observational constrains for PBHs (e.g. Calcino et al. 2018; Niikura et al. 2019), surprisingly, they rarely show the observed quantity -- the tE distribution -- for different abundances of PBHs relative to the total dark matter mass (fPBH). We present a simple calculation of how the tE distribution changes between models with and without PBHs.

One of the key phases of stellar evolution that remains poorly understood is stellar death. We lack a predictive model for how a star of a given mass explodes and what kind of remnant it leaves behind (i.e. the initial-final mass relation, IFMR). Progress has been limited due to the difficulty in finding and weighing black holes and neutron stars in large numbers. Technological advances that allow for sub-milliarcsecond astrometry in crowded fields have opened a new window for finding black holes and neutron stars: astrometric gravitational lensing. Finding and weighing a sample of compact objects with astrometric microlensing will allow us to place some of the first constraints on the present-day mass function of isolated black holes and neutron stars, their multiplicity, and their kick velocities. All of these are fundamental inputs into understanding the death phase of stellar evolution, improving supernovae models, and interpreting LIGO detections in an astrophysical context. To achieve these goals, we require large area surveys, such as the WFIRST exoplanet microlensing survey, to photometrically identify long-duration (>120 day), un-blended microlensing events as candidate compact objects. We also require high-precision astrometric follow-up monitoring using extremely large telescopes, equipped with adaptive optics, such as TMT and GMT.

The Wide Field Infrared Survey Telescope (WFIRST) will monitor $\sim 2$ deg$^2$ toward the Galactic bulge in a wide ($\sim 1-2~\mu$m) W149 filter at 15-minute cadence with exposure times of $\sim$50s for 6 seasons of 72 days each, for a total $\sim$41,000 exposures taken over $\sim$432 days, spread over the 5-year prime mission. This will be one of the deepest exposures of the sky ever taken, reaching a photon-noise photometric precision of 0.01 mag per exposure and collecting a total of $\sim 10^9$ photons over the course of the survey for a W149$_{\rm AB}\sim 21$ star. Of order $4 \times 10^7$ stars will be monitored with W149$_{\rm AB}$<21, and 10$^8$ stars with W145$_{\rm AB}$<23. The WFIRST microlensing survey will detect $\sim$54,000 microlensing events, of which roughly 1% ($\sim$500) will be due to isolated black holes, and $\sim$3% ($\sim$1600) will be due to isolated neutron stars. It will be sensitive to (effectively) isolated compact objects with masses as low as the mass of Pluto, thereby enabling a measurement of the compact object mass function over 10 orders of magnitude. Assuming photon-noise limited precision, it will detect $\sim 10^5$ transiting planets with sizes as small as $\sim 2~R_\oplus$, perform asteroseismology of $\sim 10^6$ giant stars, measure the proper motions to $\sim 0.3\%$ and parallaxes to $\sim 10\%$ for the $\sim 6 \times 10^6$ disk and bulge stars in the survey area, and directly detect $\sim 5 \times 10^3$ Trans-Neptunian objects (TNOs) with diameters down to $\sim 10$ km, as well as detect $\sim 10^3$ occulations of stars by TNOs during the survey. All of this science will completely serendipitous, i.e., it will not require modifications of the WFIRST optimal microlensing survey design. Allowing for some minor deviation from the optimal design, such as monitoring the Galactic center, would enable an even broader range of transformational science.

We have compelling evidence for stellar-mass black holes (BHs) of ~5-80 M_sun that form through the death of massive stars. We also have compelling evidence for so-called supermassive BHs (10^5-10^10 M_sun) that are predominantly found in the centers of galaxies. We have very good reason to believe there must be BHs with masses in the gap between these ranges: the first ~10^9 M_sun BHs are observed only hundreds of millions of years after the Big Bang, and all theoretically viable paths to making supermassive BHs require a stage of "intermediate" mass. However, no BHs have yet been reliably detected in the 100-10}^5 M_sun mass range. Uncovering these intermediate-mass BHs of 10^3-10^5 M_sun is within reach in the coming decade. In this white paper we highlight the crucial role that 30-m class telescopes will play in dynamically detecting intermediate-mass black holes, should they exist.

PSR B0919+06 is known for its abnormal emission phenomenon, where the pulse emission window occasionally shifts progressively in longitude and returns afterwards. The physical mechanism behind this phenomenon is still under investigation. In this paper, we present our ultra-wideband observation of this pulsar using the Five-hundred-meter Aperture Spherical radio Telescope (FAST), with simultaneous measurements in the frequency ranges 280-780 MHz and 1250-1550 MHz. We have identified three abnormal events, each of which becomes less apparent as the frequency decreases. At 1400 MHz, the averaged profile slightly shifted after the first and third abnormal events, implying a relationship between abnormal event and profile variation. We also found a linear trend in the left-edge position of the averaged profiles between the first and third event as well as after the third event, suggesting the existence of a slow-drifting mode between the two major events. The second event has a comparatively small shift in phase and is thus categorized as a `small flare state'. During the third event, a sequence of approximately nine pulses was seen to significantly weaken in all frequency bands, likely associated with the pseudo-nulling observed at 150 MHz. A three-component de-composition analysis of the normal averaged profiles shows that the trailing component is dominant at our observing frequencies, while the centre component has a comparatively steeper spectrum. We found the overall flux density in an abnormal event to slightly differ from that in an ordinary state, and the difference shows a frequency dependence. A comparison of the normal, abnormal and dimmed averaged profile indicates that the leading component is likely to be stable in all states.

Rotating radio transients (RRATs) are peculiar astronomical objects whose emission mechanism remains under investigation. In this paper, we present observations of three RRATs, J1538+2345, J1854+0306 and J1913+1330, observed with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Specifically, we analyze the mean pulse profiles and temporal flux density evolutions of the RRATs. Owing to the high sensitivity of FAST, the derived burst rates of the three RRATs are higher than those in previous reports. RRAT J1854+0306 exhibited a time-dynamic mean pulse profile, whereas RRAT J1913+1330 showed distinct radiation and nulling segments on its pulse intensity trains. The mean pulse profile variation with frequency is also studied for RRAT J1538+2345 and RRAT J1913+1330, and the profiles at different frequencies could be well fitted with a cone-core model and a conal-beam model, respectively.

With the largest dish Five-hundred-meter Aperture Spherical radio Telescope (FAST), both the mean and single pulses of PSR B2016$+$28, especially including the single-pulse structure, are investigated in detail in this study. The mean pulse profiles at different frequencies can be well fitted in a conal model, and the peak separation of intensity-dependent pulse profiles increases with intensity. The integrated pulses are obviously frequency dependent (pulse width decreases by $\sim\,20\%$ as frequency increases from 300 MHz to 750 MHz), but the structure of single pulses changes slightly (the corresponding correlation scale decreases by only $\sim\,1\%$). This disparity between mean and single pulses provides independent evidence for the existence of the RS-type vacuum inner gap, indicating a strong bond between particles on the pulsar surface. Diffused drifting sub-pulses are analyzed. The results show that the modulation period along pulse series ($P_3$) is positively correlated to the separation between two adjacent sub-pulses ($P_2$). This correlation may hint a rough surface on the pulsar, eventually resulting in the irregular drift of sparks. All the observational results may have significant implications in the dynamics of pulsar magnetosphere and are discussed extensively in this paper.

The Five-hundred-meter Aperture Spherical radio Telescope (FAST) was completed with its main structure installed on September 25, 2016, after which it entered the commissioning phase. This paper aims to introduce the commissioning progress of the FAST over the past two years. To improve its operational reliability and ensure effective observation time, FAST has been equipped with a real-time information system for the active reflector system and hierarchical commissioning scheme for the feed support system, which ultimately achieves safe operation of the two systems. For meeting the high-performance indices, a high-precision measurement system was set up based on the effective control methods that were implemented for the active reflector system and feed support system. Since the commissioning of the FAST, a low-frequency ultra-wideband receiver and 19-beam 1.05-1.45 GHz receiver have been mainly used. Telescope efficiency, pointing accuracy, and system noise temperature were completely tested and ultimately achieved the acceptance indices of the telescope. The FAST has been in the process of national acceptance preparations and has begun to search for pulsars. In the future, it will still strive to improve its capabilities and expand its application prospects.

As the closest example of a galactic nucleus, the Galactic center (GC) presents an exquisite laboratory for learning about supermassive black holes (SMBH) and their environment. We describe several exciting new research directions that, over the next 10 years, hold the potential to answer some of the biggest scientific questions raised in recent decades: Is General Relativity (GR) the correct description for supermassive black holes? What is the nature of star formation in extreme environments? How do stars and compact objects dynamically interact with the supermassive black hole? What physical processes drive gas accretion in low-luminosity black holes? We describe how the high sensitivity, angular resolution, and astrometric precision offered by the next generation of large ground-based telescopes with adaptive optics will help us answer these questions. First, it will be possible to obtain precision measurements of stellar orbits in the Galaxy's central potential, providing both tests of GR in the unexplored regime near a SMBH and measurements of the extended dark matter distribution that is predicted to exist at the GC. Second, we will probe stellar populations at the GC to significantly lower masses than are possible today, down to brown dwarfs. Their structure and dynamics will provide an unprecedented view of the stellar cusp around the SMBH and will distinguish between models of star formation in this extreme environment. This increase in depth will also allow us to measure the currently unknown population of compact remnants at the GC by observing their effects on luminous sources. Third, uncertainties on the mass of and distance to the SMBH can be improved by a factor of $\sim$10. Finally, we can also study the near-infrared accretion onto the black hole at unprecedented sensitivity and time resolution, which can reveal the underlying physics of black hole accretion.

The stellar initial mass function (IMF) is a fundamental property of star formation, offering key insight into the physics driving the process as well as informing our understanding of stellar populations, their by-products, and their impact on the surrounding medium. While the IMF appears to be fairly uniform in the Milky Way disk, it is not yet known how the IMF might behave across a wide range of environments, such as those with extreme gas temperatures and densities, high pressures, and low metallicities. We discuss new opportunities for measuring the IMF in such environments in the coming decade with JWST, WFIRST, and thirty-meter class telescopes. For the first time, we will be able to measure the high-mass slope and peak of the IMF via direct star counts for massive star clusters across the Milky Way and Local Group, providing stringent constraints for star formation theory and laying the groundwork for understanding distant and unresolved stellar systems.

The outer regions of globular clusters can enable us to answer many fundamental questions concerning issues ranging from the formation and evolution of clusters and their multiple stellar populations to the study of stars near and beyond the hydrogen-burning limit and to the dynamics of the Milky Way. The outskirts of globular clusters are still uncharted territories observationally. A very efficient way to explore them is through high-precision proper motions of low-mass stars over a large field of view. The Wide Field InfraRed Survey Telescope (WFIRST) combines all these characteristics in a single telescope, making it the best observational tool to uncover the wealth of information contained in the clusters' outermost regions.

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.

We present improved relative astrometry for stars within the central half parsec of our Galactic Center based on data obtained with the 10 m W. M. Keck Observatory from 1995 to 2017. The new methods used to improve the astrometric precision and accuracy include correcting for local astrometric distortions, applying a magnitude dependent additive error, and more carefully removing instances of stellar confusion. Additionally, we adopt jackknife methods to calculate velocity and acceleration uncertainties. The resulting median proper motion uncertainty is 0.05 mas/yr for our complete sample of 1184 stars in the central 10'' (0.4 pc). We have detected 24 accelerating sources, 2.6 times more than the number of previously published accelerating sources, which extend out to 4'' (0.16 pc) from the black hole. Based on S0-2's orbit, our new astrometric analysis has reduced the systematic error of the supermassive black hole (SMBH) by a factor of 2. The linear drift in our astrometric reference frame is also reduced in the North-South direction by a factor of 4. We also find the first potential astrometric binary candidate S0-27 in the Galactic center. These astrometric improvements provide a foundation for future studies of the origin and dynamics of the young stars around the SMBH, the structure and dynamics of the old nuclear star cluster, the SMBH's properties derived from orbits, and tests of General Relativity (GR) in a strong gravitational field.

The Orion Nebula Cluster (ONC) is the nearest site of ongoing massive star formation, which allows us to study the kinematics and dynamics of the region in detail and constrain star formation theories. Using HST ACS/WFPC2/WFC3IR and Keck II NIRC2 data, we have measured the proper motions of 701 stars within an $\sim6'\times6'$ field of view around the center of the ONC. We have found more than 10 escaping star candidates, concentrated predominantly at the core of the cluster. The proper motions of the bound stars are consistent with a normal distribution, albeit elongated North-South along the Orion filament, with proper motion dispersions of $(\sigma_{\mu,\alpha^*}, \sigma_{\mu,\delta}) = (0.83\pm0.02,\,1.12\pm0.03)$ mas yr$^{-1}$ or intrinsic velocity dispersions of $(\sigma_{v,\alpha^*}, \sigma_{v,\delta}) = (1.57\pm0.04,\,2.12\pm0.06)$ km s$^{-1}$ assuming a distance of 400 pc to the ONC. The cluster shows no evidence for tangential-to-radial anisotropy. Our velocity dispersion profile agrees with the prediction from the observed stellar + gas density profile from Da Rio et al. (2014), indicating that the ONC is in virial equilibrium. This finding suggests that the cluster was formed with a low star formation efficiency per dynamical timescale based on comparisons with current star formation theories. Our survey also recovered high-velocity IR sources BN, x, and n in the BN/KL region. The estimated location of the first two sources $\sim500$ years ago agrees with that of the radio source I, consistent with their proposed common origin from a multi-stellar disintegration. However, source n appears to have a small proper motion and is unlikely to have been involved in the event.

Precision measurements of the stars in short-period orbits around the supermassive black hole at the Galactic Center are now being used to constrain general relativistic effects, such as the gravitational redshift and periapse precession. One of the largest systematic uncertainties in the measured orbits has been errors in the astrometric reference frame, which is derived from seven infrared-bright stars associated with SiO masers that have extremely accurate radio positions, measured in the Sgr A*-rest frame. We have improved the astrometric reference frame within 14'' of the Galactic Center by a factor of 2.5 in position and a factor of 5 in proper motion. In the new reference frame, Sgr A* is localized to within a position of 0.645 mas and proper motion of 0.03 mas/yr. We have removed a substantial rotation (2.25 degrees per decade), that was present in the previous less-accurate reference frame used to measure stellar orbits in the field. With our improved methods and continued monitoring of the masers, we predict that orbital precession predicted by General Relativity will become detectable in the next ~5 years.

OSIRIS is a near-infrared (1.0--2.4 $\mu$m) integral field spectrograph operating behind the adaptive optics system at Keck Observatory, and is one of the first lenslet-based integral field spectrographs. Since its commissioning in 2005, it has been a productive instrument, producing nearly half the laser guide star adaptive optics (LGS AO) papers on Keck. The complexity of its raw data format necessitated a custom data reduction pipeline (DRP) delivered with the instrument in order to iteratively assign flux in overlapping spectra to the proper spatial and spectral locations in a data cube. Other than bug fixes and updates required for hardware upgrades, the bulk of the DRP has not been updated since initial instrument commissioning. We report on the first major comprehensive characterization of the DRP using on-sky and calibration data. We also detail improvements to the DRP including characterization of the flux assignment algorithm; exploration of spatial rippling in the reduced data cubes; and improvements to several calibration files, including the rectification matrix, the bad pixel mask, and the wavelength solution. We present lessons learned from over a decade of OSIRIS data reduction that are relevant to the next generation of integral field spectrograph hardware and data reduction software design.

Shifting the focus of Type Ia supernova (SN Ia) cosmology to the near-infrared (NIR) is a promising way to significantly reduce the systematic errors, as the strategy minimizes our reliance on the empirical width-luminosity relation and uncertain dust laws. Observations in the NIR are also crucial for our understanding of the origins and evolution of these events, further improving their cosmological utility. Any future experiments in the rest-frame NIR will require knowledge of the SN Ia NIR spectroscopic diversity, which is currently based on a small sample of observed spectra. Along with the accompanying paper, Phillips et al. (2018), we introduce the Carnegie Supernova Project-II (CSP-II), to follow up nearby SNe Ia in both the optical and the NIR. In particular, this paper focuses on the CSP-II NIR spectroscopy program, describing the survey strategy, instrumental setups, data reduction, sample characteristics, and future analyses on the data set. In collaboration with the Harvard-Smithsonian Center for Astrophysics (CfA) Supernova Group, we obtained 661 NIR spectra of 157 SNe Ia. Within this sample, 451 NIR spectra of 90 SNe Ia have corresponding CSP-II follow-up light curves. Such a sample will allow detailed studies of the NIR spectroscopic properties of SNe Ia, providing a different perspective on the properties of the unburned material, radioactive and stable nickel produced, progenitor magnetic fields, and searches for possible signatures of companion stars.

Giant pulses (GPs) are extremely bright individual pulses of radio pulsar. In microbursts of Crab pulsar, which is an active GP emitter, zebra-pattern-like spectral structures are observed, which are reminiscent of the `zebra bands' that are observed in type IV solar radio flares. However, band spacing linearly increases with the band center frequency of $\sim5-30$\,GHz. In this study, we propose that the Crab pulsar GP can originate from the coherent instability of plasma near a light cylinder. Further, the growth of coherent instability can be attributed to the resonance observed between the cyclotron-resonant-excited wave and the background plasma oscillation. The particles can be injected into the closed-field line regions owing to magnetic reconnection near a light cylinder. These particles introduce a large amount of free energy that further causes cyclotron-resonant instability, which grows and amplifies radiative waves at frequencies close to the electron cyclotron harmonics that exhibit zebra-pattern-like spectral band structures. Further, these structures can be modulated by the resonance between the cyclotron-resonant-excited wave and the background plasma oscillation. In this scenario, the band structures of the Crab pulsar can be well fitted by a coherent instability model, where the plasma density of a light cylinder should be $\sim10^{13-15}\,\rm{cm^{-3}}$, with an estimated gradient of $>5.5\times10^5\,\rm{cm^{-4}}$. This process may be accompanied by high-energy emissions. Similar phenomena are expected to be detected in other types of GP sources that have magnetic fields of $\simeq10^6$\,G in a light cylinder.

We present a $\approx 11.5$ year adaptive optics (AO) study of stellar variability and search for eclipsing binaries in the central $\sim 0.4$ pc ($\sim 10''$) of the Milky Way nuclear star cluster. We measure the photometry of 563 stars using the Keck II NIRC2 imager ($K'$-band, $\lambda_0 = 2.124 \text{ } \mu \text{m}$). We achieve a photometric uncertainty floor of $\Delta m_{K'} \sim 0.03$ ($\approx 3\%$), comparable to the highest precision achieved in other AO studies. Approximately half of our sample ($50 \pm 2 \%$) shows variability. $52 \pm 5\%$ of known early-type young stars and $43 \pm 4 \%$ of known late-type giants are variable. These variability fractions are higher than those of other young, massive star populations or late-type giants in globular clusters, and can be largely explained by two factors. First, our experiment time baseline is sensitive to long-term intrinsic stellar variability. Second, the proper motion of stars behind spatial inhomogeneities in the foreground extinction screen can lead to variability. We recover the two known Galactic center eclipsing binary systems: IRS 16SW and S4-258 (E60). We constrain the Galactic center eclipsing binary fraction of known early-type stars to be at least $2.4 \pm 1.7\%$. We find no evidence of an eclipsing binary among the young S-stars nor among the young stellar disk members. These results are consistent with the local OB eclipsing binary fraction. We identify a new periodic variable, S2-36, with a 39.43 day period. Further observations are necessary to determine the nature of this source.

Major galaxy mergers are thought to play an important part in fuelling the growth of supermassive black holes. However, observational support for this hypothesis is mixed, with some studies showing a correlation between merging galaxies and luminous quasars and others showing no such association. Recent observations have shown that a black hole is likely to become heavily obscured behind merger-driven gas and dust, even in the early stages of the merger, when the galaxies are well separated (5 to 40 kiloparsecs). Merger simulations further suggest that such obscuration and black-hole accretion peaks in the final merger stage, when the two galactic nuclei are closely separated (less than 3 kiloparsecs). Resolving this final stage requires a combination of high-spatial-resolution infrared imaging and high-sensitivity hard-X-ray observations to detect highly obscured sources. However, large numbers of obscured luminous accreting supermassive black holes have been recently detected nearby (distances below 250 megaparsecs) in X-ray observations. Here we report high-resolution infrared observations of hard-X-ray-selected black holes and the discovery of obscured nuclear mergers, the parent populations of supermassive-black-hole mergers. We find that obscured luminous black holes (bolometric luminosity higher than 2x10^44 ergs per second) show a significant (P<0.001) excess of late-stage nuclear mergers (17.6 per cent) compared to a sample of inactive galaxies with matching stellar masses and star formation rates (1.1 per cent), in agreement with theoretical predictions. Using hydrodynamic simulations, we confirm that the excess of nuclear mergers is indeed strongest for gas-rich major-merger hosts of obscured luminous black holes in this final stage.

As a young massive cluster in the Central Molecular Zone, the Arches cluster is a valuable probe of the stellar Initial Mass Function (IMF) in the extreme Galactic Center environment. We use multi-epoch Hubble Space Telescope observations to obtain high-precision proper motion and photometric measurements of the cluster, calculating cluster membership probabilities for stars down to 1.8 M$_{\odot}$ between cluster radii of 0.25 pc -- 3.0 pc. We achieve a cluster sample with just ~8% field contamination, a significant improvement over photometrically-selected samples which are severely compromised by the differential extinction across the field. Combining this sample with K-band spectroscopy of 5 cluster members, we forward model the Arches cluster to simultaneously constrain its IMF and other properties (such as age and total mass) while accounting for observational uncertainties, completeness, mass segregation, and stellar multiplicity. We find that the Arches IMF is best described by a 1-segment power law that is significantly top-heavy: $\alpha$ = 1.80 $\pm$ 0.05 (stat) $\pm$ 0.06 (sys), where dN/dm $\propto$ m$^{-\alpha}$, though we cannot discount a 2-segment power law model with a high-mass slope only slightly shallower than local star forming regions ($\alpha$ = 2.04$^{+0.14}_{-0.19}$ $\pm$ 0.04) but with a break at 5.8$^{+3.2}_{-1.2}$ $\pm$ 0.02 M$_{\odot}$. In either case, the Arches IMF is significantly different than the standard IMF. Comparing the Arches to other young massive clusters in the Milky Way, we find tentative evidence for a systematically top-heavy IMF at the Galactic Center.

We present the analysis of the simultaneous high resolution images from the {\it Hubble Space Telescope} and Keck Adaptive Optics system of the planetary event OGLE-2012-BLG-0950 that determine that the system consists of a $0.58 \pm 0.04 \rm{M}_\odot$ host star orbited by a $39\pm 8 \rm{M}_\oplus$ planet of at projected separation of $2.54 \pm 0.23\,$AU. The planetary system is located at a distance of $2.19\pm 0.23$ kpc from Earth. This is the second microlens planet beyond the snow line with a mass measured to be in the mass range $20$--$80 \rm{M}_\oplus$. The runaway gas accretion process of the core accretion model predicts few planets in this mass range, because giant planets are thought to be growing rapidly at these masses and they rarely complete growth at this mass. So, this result suggests that the core accretion theory may need revision. This analysis also demonstrates the techniques that will be used to measure the masses of planets and their host stars by the WFIRST exoplanet microlensing survey: one-dimensional microlensing parallax combined with the separation and brightness measurement of the unresolved source and host stars to yield multiple redundant constraints on the masses and distance of the planetary system.

We summarize the red channel (2-5 micron) of the Planetary Systems Imager (PSI), a proposed second-generation instrument for the TMT. Cold exoplanets emit the majority of their light in the thermal infrared, which means these exoplanets can be detected at a more modest contrast than at other wavelengths. PSI-Red will be able to detect and characterize a wide variety of exoplanets, including radial-velocity planets on wide orbits, accreting protoplanets in nearby star-forming regions, and reflected-light planets around the nearest stars. PSI-Red will feature an imager, a low-resolution lenslet integral field spectrograph, a medium-resolution lenslet+slicer integral field spectrograph, and a fiber-fed high-resolution spectrograph.

Ground-layer adaptive optics (GLAO) systems offer the possibility of improving the "seeing" of large ground-based telescopes and increasing the efficiency and sensitivity of observations over a wide field-of-view. We explore the utility and feasibility of deploying a GLAO system at the W. M. Keck Observatory in order to feed existing and future multi-object spectrographs and wide-field imagers. We also briefly summarize science cases spanning exoplanets to high-redshift galaxy evolution that would benefit from a Keck GLAO system. Initial simulations indicate that a Keck GLAO system would deliver a 1.5x and 2x improvement in FWHM at optical (500 nm) and infrared (1.5 micron), respectively. The infrared instrument, MOSFIRE, is ideally suited for a Keck GLAO feed in that it has excellent image quality and is on the telescope's optical axis. However, it lacks an atmospheric dispersion compensator, which would limit the minimum usable slit size for long-exposure science cases. Similarly, while LRIS and DEIMOS may be able to accept a GLAO feed based on their internal image quality, they lack either an atmospheric dispersion compensator (DEIMOS) or flexure compensation (LRIS) to utilize narrower slits matched to the GLAO image quality. However, some science cases needing shorter exposures may still benefit from Keck GLAO and we will investigate the possibility of installing an ADC.

`Imaka is a ground layer adaptive optics (GLAO) demonstrator on the University of Hawaii 2.2m telescope with a 24'x18' field-of-view, nearly an order of magnitude larger than previous AO instruments. In 15 nights of observing with natural guide star asterisms ~16' in diameter, we measure median AO-off and AO-on empirical full-widths at half-maximum (FWHM) of 0''95 and 0''64 in R-band, 0''81 and 0''48 in I-band, and 0''76 and 0''44 at 1 micron. This factor of 1.5-1.7 reduction in the size of the point spread function (PSF) results from correcting both the atmosphere and telescope tracking errors. The AO-on PSF is uniform out to field positions ~5' off-axis, with a typical standard deviation in the FWHM of 0''018. Images exhibit variation in FWMM by 4.5% across the field, which has been applied as a correction to the aforementioned quantities. The AO-on PSF is also 10x more stable in time compared to the AO-off PSF. In comparing the delivered image quality to proxy measurements, we find that in both AO-off and AO-on data, delivered image quality is correlated with `imaka's telemetry, with R-band correlation coefficients of 0.68 and 0.70, respectively. At the same wavelength, the data are correlated to DIMM and MASS seeing with coefficients of 0.45 and 0.55. Our results are an essential first step to implementing facility-class, wide-field GLAO on Maunakea telescopes, enabling new opportunities to study extended astronomical sources, such as deep galaxy fields, nearby galaxies or star clusters, at high angular resolution.

Adaptive optics (AO) restore the angular resolution of ground-based telescopes, but at the cost of delivering a time- and space-varying point spread function (PSF) with a complex shape. PSF knowledge is crucial for breaking existing limits on the measured accuracy of photometry and astrometry in science observations. In this paper, we concentrate our analyses on anisoplanatism signature only onto PSF: for large-field observations (20") with single- conjugated AO, PSFs are strongly elongated due to anisoplanatism that manifests itself as three different terms for Laser-guide star (LGS) systems: angular, focal and tilt. We propose a generalized model that relies on a point-wise decomposition of the phase and encompasses the non-stationarity of LGS systems. We demonstrate it is more accurate and less computationally demanding than existing models: it agrees with end-to-end physical-optics simulations to within 0.1% of PSF measurables, such as Strehl-ratio, FWHM and fraction of variance unexplained.Secondly, we study off-axis PSF modelling is with respect to $C_n^2(h)$ profile (heights and fractional weights). For 10m class telescope, PSF morphology is estimated at 1% level as long as we model the atmosphere with at least 7 layers whose heights and weights are known respectively with 200m and 10% precision. As a verification test we used the Canada's NRC-Herzberg HeNOS testbed data, featuring four lasers. We highlight capability of retrieving off-axis PSF characteristics within 10% of fraction of variance unexplained, which complies with the expected range from the sensitivity analysis. Our new off-axis PSF modelling method lays the ground-work for testing on-sky in the near future.

X-ray-dim isolated neutron stars (XDINSs), also known as the Magnificent Seven, exhibits a Planck-like soft X-ray spectrum. In the optical/ultraviolet(UV) band, there is an excess of radiation compared to the extrapolation from the X-ray spectrum. A model of bremsstrahlung emission from a nonuniform plasma atmosphere is proposed in the regime of a strangeon star to explain the optical/UV excess and its spectral deviation. The nonuniform atmosphere is formed by the accretion of interstellar matter which moves along the magnetic field lines to polar caps, and these particles may spread out of the pole regions. The modeled electron temperatures are $\sim100-200$ eV with radiation radii $R_{\rm opt}^{\infty}\sim5-14$ km. The spectra of five sources (RX J0720.4--3125, RX J0806.4--4123, RX J1308.6+2127, RX J1605.3+3249, RX J1856.5--3754) from optical/UV to X-ray bands could be well fitted by the radiative model, and exhibit gaussian absorption lines at $\sim 100-500$ eV as would be expected. Furthermore, the surroundings (i.e., fallback disks or dusty belts) of XDINSs could be tested by future observations with the Large Optical Telescope (LOT).

Pulsar-like compact stars provide us a unique laboratory to explore properties of dense matter at supra-nuclear densities. One of the models for pulsar-like stars is that they are totally composed of "strangeons", and in this paper we studied the pulsar glitches in a strangeon star model. Strangeon stars would be solidified during cooling, and the solid stars would be natural to have glitches as the result of starquakes. Based on the starquake model established before, we proposed that when the starquake occurs, the inner motion of the star which changes the moment of inertia and has impact on the glitch sizes, is divided into plastic flow and elastic motion. The plastic flow which is induced in the fractured part of the outer layer, would move tangentially to redistribute the matter of the star and would be hard to recover. The elastic motion, on the other hand, changes its shape and would recover significantly. Under this scenario, we could understand the behaviors of glitches without significant energy releasing, including the Crab and the Vela pulsars, in an uniform model. We derive the recovery coefficient as a function of glitch size, as well as the time interval between two successive glitches as the function of the released stress. Our results show consistency with observational data under reasonable ranges of parameters. The implications on the oblateness of the Crab and the Vela pulsars are discussed.

We present adaptive-optics assisted near-infrared high-spectral resolution observations of late-type giants in the nuclear star cluster of the Milky Way. The metallicity and elemental abundance measurements of these stars offer us an opportunity to understand the formation and evolution of the nuclear star cluster. In addition, their proximity to the supermassive black hole ($\sim 0.5$ pc) offers a unique probe of the star formation and chemical enrichment in this extreme environment. We observed two stars identified by medium spectral-resolution observations as potentially having very high metallicities. We use spectral-template fitting with the PHOENIX grid and Bayesian inference to simultaneously constrain the overall metallicity, [M/H], alpha-element abundance [$\alpha$/Fe], effective temperature, and surface gravity of these stars. We find that one of the stars has very high metallicity ([M/H] $> 0.6$) and the other is slightly above solar metallicity. Both Galactic center stars have lines from scandium (Sc), vanadium (V), and yttrium (Y) that are much stronger than allowed by the PHOENIX grid. We find, using the spectral synthesis code Spectroscopy Made Easy, that [Sc/Fe] may be an order of magnitude above solar. For comparison, we also observed an empirical calibrator in NGC6791, the highest metallicity cluster known ([M/H] $\sim 0.4$). Most lines are well matched between the calibrator and the Galactic center stars, except for Sc, V, and Y, which confirms that their abundances must be anomalously high in these stars. These unusual abundances, which may be a unique signature of nuclear star clusters, offer an opportunity to test models of chemical enrichment in this region.

The state equation of strangeon matter is very stiff due to the non-relativistic nature of and the repulsing interaction between the particles, and pulsar masses as high as $\sim 3M_\odot$ would be expected. However, an adiabatic sound speed, $c_s=\sqrt{\partial P/\partial \rho}$, is usually superluminal for strangeon matter, and dynamic response of strangeon star (e.g., binary merger) could not be tractable in calculation. We examine signal propagation in strangeon matter, and calculate the propagation speed, $c_{\rm signal}$, in reality. It is found that as the causality condition is satisfied, i.e., $c_{\rm signal}<c$, and the signal speed as a function of stellar radius is presented.

A precise extinction law is a critical input when interpreting observations of highly reddened sources such as young star clusters and the Galactic Center (GC). We use Hubble Space Telescope observations of a region of moderate extinction and a region of high extinction to measure the optical and near-infrared extinction law (0.8 $\mu$m -- 2.2 $\mu$m). The moderate extinction region is the young massive cluster Westerlund 1 (Wd1; A$_{Ks} \sim$ 0.6 mag), where 453 proper motion-selected main-sequence stars are used to measure the shape of the extinction law. To quantify the shape we define the parameter $\mathcal{S}_{1/\lambda}$, which behaves similarly to a color excess ratio but is continuous as a function of wavelength. The high extinction region is the GC (A$_{Ks} \sim$ 2.5 mag), where 819 red clump stars are used to determine the normalization of the law. The best-fit extinction law is able to reproduce the Wd1 main sequence colors, which previous laws misestimate by 10%-30%. The law is inconsistent with a single power law, even when only the near-infrared filters are considered, and has A$_{F125W}$/A$_{Ks}$ and A$_{F814W}$/A$_{Ks}$ values that are 18% and 24% larger than the commonly used \citet{Nishiyama:2009fc} law, respectively. Using the law we recalculate the Wd1 distance to be 3896 $\pm$ 328 pc from published observations of eclipsing binary W13. This new extinction law should be used for highly reddened populations in the Milky Way, such as the Quintuplet cluster and Young Nuclear Cluster. A python code is provided to generate the law for future use.

We present a kinematic study of the nuclear stellar disk in M31 at infrared wavelengths using high spatial resolution integral field spectroscopy. The spatial resolution achieved, FWHM = 0."12 (0.45 pc at the distance of M31), has only previously been equaled in spectroscopic studies by space-based long-slit observations. Using adaptive optics-corrected integral field spectroscopy from the OSIRIS instrument at the W. M. Keck Observatory, we map the line-of-sight kinematics over the entire old stellar eccentric disk orbiting the supermassive black hole (SMBH) at a distance of r<4 pc. The peak velocity dispersion is 381+/-55 km/s , offset by 0.13 +/- 0.03 from the SMBH, consistent with previous high-resolution long-slit observations. There is a lack of near-infrared (NIR) emission at the position of the SMBH and young nuclear cluster, suggesting a spatial separation between the young and old stellar populations within the nucleus. We compare the observed kinematics with dynamical models from Peiris & Tremaine (2003). The best-fit disk orientation to the NIR flux is [$\theta_l$, $\theta_i$, $\theta_a$] = [-33 +/- 4$^{\circ}$, 44 +/- 2$^{\circ}$, -15 +/- 15$^{\circ}$], which is tilted with respect to both the larger-scale galactic disk and the best-fit orientation derived from optical observations. The precession rate of the old disk is $\Omega_P$ = 0.0 +/- 3.9 km/s/pc, lower than the majority of previous observations. This slow precession rate suggests that stellar winds from the disk will collide and shock, driving rapid gas inflows and fueling an episodic central starburst as suggested in Chang et al. (2007).

The population of young stars near the Supermassive black hole (SMBH) in the Galactic Center (GC) has presented an unexpected challenge to theories of star formation. Kinematics measurements of these stars have revealed a stellar disk structure (with an apparent 20% disk membership) that has provided important clues to the origin of these mysterious young stars. However many of the apparent disk properties are difficult to explain, including the low disk membership fraction and the high eccentricities, given the youth of this population. Thus far, all efforts to derive the properties of this disk have made the simplifying assumption that stars at the GC are single stars. Nevertheless, stellar binaries are prevalent in our Galaxy, and recent investigations suggested that they may also be abundant in the Galactic Center. Here we show that binaries in the disk can largely alter the apparent orbital properties of the disk. The motion of binary members around each other adds a velocity component, which can be comparable to the magnitude of the velocity around the SMBH in the GC. Thus neglecting the contribution of binaries can significantly vary the inferred stars' orbital properties. While the disk orientation is unaffected the apparent disk's 2D width is increased to about 11.2deg, similar to the observed width. For a population of stars orbiting the SMBH with zero eccentricity, unaccounted for binaries will create a wide apparent eccentricity distribution with an average of 0.23.This is consistent with the observed average eccentricity of the stars' in the disk. We suggest that this high eccentricity value, which poses a theoretical challenge, may be an artifact of binary stars. Finally our results suggest that the actual disk membership might be significantly higher than the one inferred by observations that ignore the contribution of binaries, alleviating another theoretical challenge.

The star S0-2, which orbits the supermassive black hole (SMBH) in our Galaxy with a period of 16 years, provides the strongest constraint on both the mass of the SMBH and the distance to the Galactic center. S0-2 will soon provide the first measurement of relativistic effects near a SMBH. We report the first limits on the binarity of S0-2 from radial velocity monitoring, which has implications for both understanding its origin and robustness as a probe of the central gravitational field. With 87 radial velocity measurements, which include 12 new observations presented, we have the data set to look for radial velocity variations from S0-2's orbital model. Using a Lomb-Scargle analysis and orbit fitting for potential binaries, we detect no radial velocity variation beyond S0-2's orbital motion and do not find any significant periodic signal. The lack of a binary companion does not currently distinguish between different formation scenarios for S0-2. The upper limit on the mass of a companion star ($M_{\text{comp}}$) still allowed by our results has a median upper limit of $M_{\text{comp}}$ $\sin i \leq$ 1.6 M$_{\odot}$ for periods between 1 and 150 days, the longest period to avoid tidal break up of the binary. We also investigate the impact of the remaining allowed binary system on the measurement of the relativistic redshift at S0-2's closest approach in 2018. While binary star systems are important to consider for this experiment, we find plausible binaries for S0-2 will not alter a 5$\sigma$ detection of the relativistic redshift.

This paper studies acceleration processes of background thermal electrons in X-ray binary jets via turbulent stochastic interactions and shock collisions. By considering turbulent magnetized jets mixed with fluctuation magnetic fields and ordered, large-scale one, and numerically solving the transport equation along the jet axis, we explore the influence of such as magnetic turbulence, electron injections, location of an acceleration region, and various cooling rates on acceleration efficiency. The results show that (1) the existence of the dominant turbulent magnetic fields in the jets is necessary to accelerate background thermal electrons to relativistic energies. (2) Acceleration rates of electrons depend on magnetohydrodynamic turbulence types, from which the turbulence type with a hard slope can accelerate electrons more effectively. (3) An effective acceleration region should be located at the distance $>10^3R_{\rm g}$ away from the central black hole ($R_{\rm g}$ being a gravitational radius). As a result of acceleration rates competing with various cooling rates, background thermal electrons obtain not only an increase in their energies but also their spectra are broadened beyond the given initial distribution to form a thermal-like distribution. (4) The acceleration mechanisms explored in this work can reasonably provide the electron maximum energy required for interpreting high-energy $\gamma$-ray observations from microquasars, but it needs to adopt some extreme parameters in order to predict a possible very high-energy $\gamma$-ray signal.

By the end of 2014, a cosmic ray muon telescope was installed at Zhongshan Station in Antarctic and has been continuously collecting data since then. It is the first surface muon telescope to be built in Antarctic. In June 2015, five CMEs were ejected towards the Earth initiating a big large Forbush decrease (FD) event. We conduct a comprehensive study of the galactic cosmic ray intensity fluctuations during the FD using the data from cosmic ray detectors of multiple stations (Zhongshan, McMurdo, South Polar and Nagoya) and he solar wind measurements from ACE and WIND. A pre-increase before the shock arrival was observed. Distinct differences exist in the timelines of the galactic cosmic ray recorded by the neutron monitors and the muon telescopes. FD onset for Zhongshan muon telescope is delayed (2.5h) with respect to SSC onset. This FD had a profile of four-step decrease. The traditional one- or two-step classification of FDs was inadequate to explain this FD.

Stable micro-nucleus is 2-flavored (u and d), whereas stable macro-nucleus could be 3-flavored (u, d and s) if the light flavor symmetry restores there. Nucleons are the constituent of a nucleus, while strangeons are named as the constituent of 3-flavored baryonic matter. Gravity-compressed baryonic object created after core-collapse supernova could be strangeon star if the energy scale (~0.5 GeV) cannot be high enough for quark deconfinement and if there occurs 3-flavor symmetry restoration. Strangeon stars are explained here, including their formation and manifestation/identification. Much work, coupled with effective {micro-model} of strangeon matter, is needed to take advantage of the unique opportunities advanced facilities will provide.

The liquid drop model of 2-flavored ($u$ and $d$) nucleus is well known and successful, analogically, a similar drop model for 3-flavored ($u$, $d$ and $s$) nucleus is developed. A 3-flavored nucleus conjectured could be stable only if its baryon number is lager than a critical one, $A_{\rm c}$, in which strangeons are the constituent as an analogy of nucleons for nucleus. We try to model strangeon matter in a sense of phenomenological liquid drop, with two free parameters: the mass per bayron of a strangeon in vacuum, $M$, and potential deep between strangeons, $\epsilon$. It is found that, for $M\sim$ GeV and $\epsilon\sim 100$ MeV, strangeon matter could be stable and its critical number could be as low as $A_{\rm c}=300$.

At a distance of 2~pc, our nearest brown dwarf neighbor, Luhman 16 AB, has been extensively studied since its discovery 3 years ago, yet its most fundamental parameter -- the masses of the individual dwarfs -- has not been constrained with precision. In this work we present the full astrometric orbit and barycentric motion of Luhman 16 AB and the first precision measurements of the individual component masses. We draw upon archival observations spanning 31 years from the European Southern Observatory (ESO) Schmidt Telescope, the Deep Near-Infrared Survey of the Southern Sky (DENIS), public FORS2 data on the Very Large Telescope (VLT), and new astrometry from the Gemini South Multiconjugate Adaptive Optics System (GeMS). Finally, we include three radial velocity measurements of the two components from VLT/CRIRES, spanning one year. With this new data sampling a full period of the orbit, we use a Markov Chain Monte Carlo algorithm to fit a 16-parameter model incorporating mutual orbit and barycentric motion parameters and constrain the individual masses to be~$27.9^{+1.1}_{-1.0}$~$M_{J}$ for the T dwarf and~$34.2^{+1.3}_{-1.1}$~$M_{J}$ for the L dwarf. Our measurements of Luhman 16 AB's mass ratio and barycentric motion parameters are consistent with previous estimates in the literature utilizing recent astrometry only. The GeMS-derived measurements of the Luhman 16 AB separation in 2014-2015 agree closely with Hubble Space Telescope (HST) measurements made during the same epoch Bedin et al. 2017, and the derived mutual orbit agrees with those measurements to within the HST uncertainties of $0.3 - 0.4$ milliarcseconds.

Integrated pulse profiles at 8.6~GHz obtained with the Shanghai Tian Ma Radio Telescope (TMRT) are presented for a sample of 26 pulsars. Mean flux densities and pulse width parameters of these pulsars are estimated. For eleven pulsars these are the first high-frequency observations and for a further four, our observations have a better signal-to-noise ratio than previous observations. For one (PSR J0742-2822) the 8.6~GHz profiles differs from previously observed profiles. A comparison of 19 profiles with those at other frequencies shows that in nine cases the separation between the outmost leading and trailing components decreases with frequency, roughly in agreement with radius-to-frequency mapping, whereas in the other ten the separation is nearly constant. Different spectral indices of profile components lead to the variation of integrated pulse profile shapes with frequency. In seven pulsars with multi-component profiles, the spectral indices of the central components are steeper than those of the outer components. For the 12 pulsars with multi-component profiles in the high-frequency sample, we estimate the core width using gaussian fitting and discuss the width-period relationship.

We present new Adaptive Optics (AO) imaging and spectroscopic measurements of Galactic Center source G1 from W. M. Keck Observatory. Our goal is to understand its nature and relationship to G2, which is the first example of a spatially-resolved object interacting with the supermassive black hole (SMBH). Both objects have been monitored with AO for the past decade (2003 - 2014) and are comparatively close to the black hole ($a_{\rm{min}} \sim$200-300 AU) on very eccentric orbits ($e_{\rm{G1}}\sim$0.99; $e_{\rm{G2}}\sim$0.96). While G2 has been tracked before and during periapse passage ($T_{0} \sim$ 2014.2), G1 has been followed since soon after emerging from periapse ($T_{0} \sim$ 2001.3). Our observations of G1 double the previously reported observational time baseline, which improves its orbital parameter determinations. G1's orbital trajectory appears to be in the same plane as that of G2, but with a significantly different argument of periapse ($\Delta\omega$ = 21$\pm$4 degrees). This suggests that G1 is an independent object and not part of a gas stream containing G2 as has been proposed. Furthermore, we show for the first time that: (1) G1 is extended in the epochs closest to periapse along the direction of orbital motion and (2) G1 becomes significantly smaller over time, (450 AU in 2004 to less than 170 AU in 2009). Based on these observations, G1 appears to be the second example of an object tidally interacting with a SMBH. G1's existence 14 years after periapse, along with its compactness in epochs further from the time of periapse, suggest that this source is stellar in nature.

Motion of short-period stars orbiting the supermassive black hole in our Galactic Center has been monitored for more than 20 years. These observations are currently offering a new way to test the gravitational theory in an unexplored regime: in a strong gravitational field, around a supermassive black hole. In this proceeding, we present three results: (i) a constraint on a hypothetical fifth force obtained by using 19 years of observations of the two best measured short-period stars S0-2 and S0-38 ; (ii) an upper limit on the secular advance of the argument of the periastron for the star S0-2 ; (iii) a sensitivity analysis showing that the relativistic redshift of S0-2 will be measured after its closest approach to the black hole in 2018.

The neutrino burst detected during supernova SN1987A is explained in a strangeon star model, in which it is proposed that a pulsar-like compact object is composed of strangeons (strangeon: an abbreviation of "strange nucleon"). A nascent strangeon star's initial internal energy is calculated, with the inclusion of pion excitation (energy around 10^53 erg, comparable to the gravitational binding energy of a collapsed core). A liquid-solid phase transition at temperature ~ 1-2 MeV may occur only a few ten-seconds after core-collapse, and the thermal evolution of strangeon star is then modeled. It is found that the neutrino burst observed from SN 1987A could be re-produced in such a cooling model.

In this Letter, we demonstrate that short-period stars orbiting around the supermassive black hole in our Galactic Center can successfully be used to probe the gravitational theory in a strong regime. We use 19 years of observations of the two best measured short-period stars orbiting our Galactic Center to constrain a hypothetical fifth force that arises in various scenarios motivated by the development of a unification theory or in some models of dark matter and dark energy. No deviation from General Relativity is reported and the fifth force strength is restricted to an upper 95% confidence limit of $\left|\alpha\right| < 0.016$ at a length scale of $\lambda=$ 150 astronomical units. We also derive a 95% confidence upper limit on a linear drift of the argument of periastron of the short-period star S0-2 of $\left|\dot \omega_\textrm{S0-2} \right|< 1.6 \times 10^{-3}$ rad/yr, which can be used to constrain various gravitational and astrophysical theories. This analysis provides the first fully self-consistent test of the gravitational theory using orbital dynamic in a strong gravitational regime, that of a supermassive black hole. A sensitivity analysis for future measurements is also presented.

About 0.4s after the Laser Interferometer Gravitational-Wave Observatory (LIGO) detected a transient gravitational-wave (GW) signal GW150914, the Fermi Gamma-ray Burst Monitor (GBM) also found a weak electromagnetic transient (GBM transient 150914). Time and location coincidences favor a possible association between GW150904 and GBM transient 150914. Under this possible association, we adopt Fermi's electromagnetic (EM) localization and derive constraints on possible violations of the Weak Equivalence Principle (WEP) from the observations of two events. Our calculations are based on four comparisons: (1)The first is the comparison of the initial GWs detected at the two LIGO sites. From the different polarizations of these initial GWs, we obtain a limit on any difference in the parametrized post-Newtonian (PPN) parameter $\Delta\gamma\lesssim 10^{-10}$. (2) The second is a comparison of GWs and possible EM waves. Using a traditional super-Eddington accretion model for GBM transient 150914, we again obtain an upper limit $\Delta\gamma\lesssim 10^{-10}$. Compared with previous results for photons and neutrinos, our limits are five orders of magnitude stronger than those from PeV neutrinos in blazar flares, and seven orders stronger than those from MeV neutrinos in SN1987A. (3) The third is a comparison of GWs with different frequencies in the range [35 Hz, 250 Hz]. (4) The fourth is a comparison of EM waves with different energies in the range [1 keV, 10 MeV]. These last two comparisons lead to an even stronger limit, $\Delta\gamma\lesssim 10^{-8}$. Our results highlight the potential of multi-messenger signals exploiting different emission channels to strengthen existing tests of the WEP.

The emission geometries, e.g. the emission region height, the beam shape, and radius-to-frequency mapping, are important predictions of pulsar radiation model. The multi-band radio observations carry such valuable information. In this paper, we study two bright pulsars, (PSRs B0329+54 and B1642-03) and observe them in high frequency (2.5 GHz, 5 GHz, and 8 GHz). The newly acquired data together with historical archive provide an atlas of multi-frequency profiles spanning from 100 MHz to 10 GHz. We study the frequency evolution of pulse profiles and the radiation regions with the these data. We firstly fit the pulse profiles with Gaussian functions to determine the phase of each component, and then calculate the radiation altitudes of different emission components and the radiation regions. We find that the inverse Compton scattering (ICS) model can reproduce the radiation geometry of these two pulsars. But for PSR B0329+54 the radiation can be generated in either annular gap (AG) or core gap (CG), while the radiation of PSR B1642-03 can only be generated in the CG. This difference is caused by the inclination angle and the impact angle of these two pulsars. The relation of beaming angle (the angle between the radiation direction and the magnetic axis) and the radiation altitudes versus frequency is also presented by modelling the beam-frequency evolution in the ICS model. The multi-band pulse profiles of these two pulsars can be described well by the ICS model combined with the CG and AG.

X-ray dim isolated neutron stars (XDINSs) are characterized by Planckian spectra in X-ray bands, but show optical/ultraviolet(UV) excesses which are the measured photometry exceeding that is extrapolated from X-ray spectra. To solve this problem, a radiative model of bremsstrahlung emission from a plasma atmosphere is established in the regime of strangeon star. A strangeon star atmosphere could simply be regarded as the upper layer of a normal neutron star. This plasma atmosphere, formed and maintained by the ISM-accreted matter due to the so-called strangeness barrier, is supposed to be of two-temperature. All the seven XDINS spectra could be well fitted by the radiative model, from optical/UV to X-ray bands. The fitted radiation radii of XDINSs are from 7 to 13km, while the modelled electron temperatures are between 50 and 250eV, except RX J0806.4-4123 with a radiation radius ~3.5km, indicating that this source could be a low-mass strangeon star candidate. This strangeon star model could further be tested by soft X-ray polarimetry, such as the Lightweight Asymmetry and Magnetism Probe which is expected to work on Chinese space station around 2020.

We study the properties of polarized radiation in turbulent magnetic fields from X-ray binary jets. These turbulent magnetic fields are composed of large- and small-scale configurations, which result in the polarized jitter radiation when the characteristic length of turbulence is less than the non-relativistic Larmor radius. On the contrary, the polarized synchrotron emission occurs, corresponding to a large-scale turbulent environment. We calculate the spectral energy distributions and the degree of polarization for a general microquasar. Numerical results show that turbulent magnetic field configurations can indeed provide a high degree of polarization, which does not mean that a uniform, large-scale magnetic field structure exists. The model is applied to investigate the properties of polarized radiation of black hole X-ray binary Cygnus X--1. Under the constraint of multiband observations of this source, our studies demonstrate that the model can explain the high polarization degree at MeV tail and predict the highly polarized properties at high-energy $\gamma$-ray region, and that the dominant small-scale turbulent magnetic field plays an important role for explaining the highly polarized observation at hard X-ray/soft $\gamma$-ray bands. This model can be tested by polarization observations of upcoming polarimeters at high-energy $\gamma$-ray bands.

Tensions between several cosmic observations were found recently, such as the inconsistent values of $H_{0}$ (or $\sigma_{8}$) were indicated by the different cosmic observations. Introducing the massive neutrinos in $\Lambda$CDM could potentially solve the tensions. Viable $f(R)$ gravity producing $\Lambda$CDM background expansion with massive neutrinos is investigated in this paper. We fit the current observational data: Planck-2015 CMB, RSD, BAO and SNIa to constrain the mass of neutrinos in viable $f(R)$ theory. The constraint results at 95\% confidence level are: $\Sigma m_\nu<0.202$ eV for the active neutrino case, $m_{\nu, sterile}^{eff}<0.757$ eV with $N_{eff}<3.22$ for the sterile neutrino case. For the effects by the mass of neutrinos, the constraint results on model parameter at 95\% confidence level become $f_{R0}\times 10^{-6}> -1.89$ and $f_{R0}\times 10^{-6}> -2.02$ for two cases, respectively. It is also shown that the fitting values of several parameters much depend on the neutrino properties, such as the cold dark matter density, the cosmological quantities at matter-radiation equality, the neutrino density and the fraction of baryonic mass in helium. At last, the constraint result shows that the tension between direct and CMB measurements of $H_0$ gets slightly weaker in the viable $f(R)$ model than that in the base $\Lambda$CDM model.

We report the discovery of an extremely red planetary-mass companion to 2MASS J22362452+4751425, a $\approx$0.6 $M_{\odot}$ late-K dwarf likely belonging to the $\sim$120 Myr AB Doradus moving group. 2M2236+4751 b was identified in multi-epoch NIRC2 adaptive optics imaging at Keck Observatory at a separation of 3.7$"$, or 230 $\pm$ 20 AU in projection at the kinematic distance of 63 $\pm$ 5 pc to its host star. Assuming membership in the AB Dor group, as suggested from its kinematics, the inferred mass of 2M2236+4751 b is 11-14 $M_\mathrm{Jup}$. Follow-up Keck/OSIRIS $K$-band spectroscopy of the companion reveals strong CO absorption similar to other faint red L dwarfs and lacks signs of methane absorption despite having an effective temperature of $\approx$900-1200 K. With a ($J$-$K$)$_\mathrm{MKO}$ color of 2.69 $\pm$ 0.12 mag, the near-infrared slope of 2M2236+4751 b is redder than all of the HR 8799 planets and instead resembles the $\approx$23 Myr isolated planetary-mass object PSO J318.5-22, implying that similarly thick photospheric clouds can persist in the atmospheres of giant planets at ages beyond 100 Myr. In near-infrared color-magnitude diagrams, 2M2236+4751 b is located at the tip of the red L dwarf sequence and appears to define the "elbow" of the AB Dor substellar isochrone separating low-gravity L dwarfs from the cooler young T dwarf track. 2M2236+4751 b is the reddest substellar companion to a star and will be a valuable benchmark to study the shared atmospheric properties of young low-mass brown dwarfs and extrasolar giant planets.

Due to its proximity, the Milky Way nuclear star cluster provides us with a wealth of data not available in other galactic nuclei. In particular, with adaptive optics, we can observe the detailed properties of individual stars, which can offer insight into the origin and evolution of the cluster. We summarize work on the central parsec of the Galactic center based on imaging and spectroscopic observations at the Keck and Gemini telescopes. These observations include stellar positions in two dimension and the velocity in three dimensions. Spectroscopic observations also enable measurements of the physical properties of individual stars, such as the spectral type and in some cases the effective temperature, metallicity, and surface gravity. We present a review of our latest measurements of the density profiles and luminosity functions of the young and old stars in this region. These observations show a complex stellar population with a young (4-6 Myr) compact star cluster in the central 0.5 pc embedded in an older and much more massive nuclear star cluster. Surprisingly, the old late-type giants do not show a cusp profile as long been expected from theoretical work. The solution to the missing cusp problem may offer us insight into the dynamical evolution of the cluster. Finally, we also discuss recent work on the metallicity of stars in this region and how they might be used to trace their origin. The nuclear star cluster shows a large variation in metallicity ([M/Fe]). The majority of the stars have higher than solar metallicity, with about 6\% having [M/Fe] $< -0.5$. These observations indicate that the NSC was not built from the globular clusters that we see today. The formation of the nuclear star cluster is more likely from the inward migration of gas originating in the disk of the Milky Way.

The coupled dark energy model provides a possible approach to mitigate the coincidence problem of cosmological standard model. Here, the coupling term is assumed as $\bar{Q}=3H\xi_x\bar{\rho}_x$, which is related to the interaction rate and energy density of dark energy. We derive the background and perturbation evolution equations for several coupled models. Then, we test these models by currently available cosmic observations which include cosmic microwave background radiation from Planck 2015, baryon acoustic oscillation, type Ia supernovae, $f\sigma_8(z)$ data points from redshift-space distortions, and weak gravitational lensing. The constraint results tell us there is no evidence of interaction at 2$\sigma$ level, it is very hard to distinguish different coupled models from other ones.

We present constraints on the variability and binarity of young stars in the central 10 arcseconds (~0.4 pc) of the Milky Way Galactic Center (GC) using Keck Adaptive Optics data over a 12 year baseline. Given our experiment's photometric uncertainties, at least 36% of our sample's known early-type stars are variable. We identified eclipsing binary systems by searching for periodic variability. In our sample of spectroscopically confirmed and likely early-type stars, we detected the two previously discovered GC eclipsing binary systems. We derived the likely binary fraction of main sequence, early-type stars at the GC via Monte Carlo simulations of eclipsing binary systems, and find that it is at least 32% with 90% confidence.

The late-type stellar population in the Galactic Center was first predicted to reside in a dynamically relaxed cusp (power law slope ranging from 3/2 to 7/4). However, observations - which rely on models to correct for projection effects - have suggested a flat distribution instead. The need for this correction is due to the lack of information regarding the line-of-sight distances. With a two decade long baseline in astrometric measurements, we are now able to measure significant projected radial accelerations, six of which are newly reported here, that directly constrain line-of-sight distances. Here we present a statistical approach to take advantage of this information and more accurately constrain the shape of the radial density profile of the late-type stellar population in the Galactic Center.

This work focuses on a sample of seven extremely late-time X-ray flares with peak time $t_{\rm p} > 10^4 {\rm s}$, among which two flares can be confirmed as the late-time activity of central engine. The main purpose is to investigate the mechanism of such late-time flares based on the internal origin assumption. In the hyper-accreting black hole (BH) scenario, we study the possibility of two well-known mechanisms as the central engine to power such X-ray flares, i.e., the neutrino-antineutrino annihilation and the Blandford-Znajek (BZ) process. Our results show that the annihilation luminosity is far below the observational data. Thus, the annihilation mechanism cannot account for such late-time flares. For the BZ process, if the role of outflows is taken into consideration, the inflow mass rate near the horizon will be quite low such that the magnetic field will probably be too weak to power the observed X-ray flares. We therefore argue that, for the late-time flares with internal origin, the central engine is unlikely to be associated with BHs. On the contrary, a fast rotating neutron star with strong bipolar magnetic fields may be responsible for such flares.

We present the integration status for `imaka, the ground-layer adaptive optics (GLAO) system on the University of Hawaii 2.2-meter telescope on Maunakea, Hawaii. This wide-field GLAO pathfinder system exploits Maunakea's highly confined ground layer and weak free-atmosphere to push the corrected field of view to ~1/3 of a degree, an areal field approaching an order of magnitude larger than any existing or planned GLAO system, with a FWHM ~ 0.33 arcseconds in the visible and near infrared. We discuss the unique design aspects of the instrument, the driving science cases and how they impact the system, and how we will demonstrate these cases on the sky.

The Thirty Meter Telescope (TMT) first light instrument IRIS (Infrared Imaging Spectrograph) will complete its preliminary design phase in 2016. The IRIS instrument design includes a near-infrared (0.85 - 2.4 micron) integral field spectrograph (IFS) and imager that are able to conduct simultaneous diffraction-limited observations behind the advanced adaptive optics system NFIRAOS. The IRIS science cases have continued to be developed and new science studies have been investigated to aid in technical performance and design requirements. In this development phase, the IRIS science team has paid particular attention to the selection of filters, gratings, sensitivities of the entire system, and science cases that will benefit from the parallel mode of the IFS and imaging camera. We present new science cases for IRIS using the latest end-to-end data simulator on the following topics: Solar System bodies, the Galactic center, active galactic nuclei (AGN), and distant gravitationally-lensed galaxies. We then briefly discuss the necessity of an advanced data management system and data reduction pipeline.

While dozens of stellar mass black holes have been discovered in binary systems, isolated black holes have eluded detection. Their presence can be inferred when they lens light from a background star. We attempt to detect the astrometric lensing signatures of three photometrically identified microlensing events, OGLE-2011-BLG-0022, OGLE-2011-BLG-0125, and OGLE-2012-BLG-0169 (OB110022, OB110125, and OB120169), located toward the Galactic Bulge. These events were selected because of their long durations, which statistically favors more massive lenses. Astrometric measurements were made over 1-2 years using laser-guided adaptive optics observations from the W. M. Keck Observatory. Lens model parameters were first constrained by the photometric light curves. The OB120169 light curve is well-fit by a single-lens model, while both OB110022 and OB110125 light curves favor binary-lens models. Using the photometric fits as prior information, no significant astrometric lensing signal was detected and all targets were consistent with linear motion. The significant lack of astrometric signal constrains the lens mass of OB110022 to 0.05-1.79 Msun in a 99.7% confidence interval, which disfavors a black hole lens. Fits to OB110125 yielded a reduced Einstein crossing time and insufficient observations during the peak, so no mass limits were obtained. Two degenerate solutions exist for OB120169, which have a lens mass between 0.2-38.8 Msun and 0.4-39.8 Msun for a 99.7% confidence interval. Follow-up observations of OB120169 will further constrain the lens mass. Based on our experience, we use simulations to design optimal astrometric observing strategies and show that, with more typical observing conditions, detection of black holes is feasible.

We present new, more precise measurements of the mass and distance of our Galaxy's central supermassive black hole, Sgr A*. These results stem from a new analysis that more than doubles the time baseline for astrometry of faint stars orbiting Sgr A*, combining two decades of speckle imaging and adaptive optics data. Specifically, we improve our analysis of the speckle images by using information about a star's orbit from the deep adaptive optics data (2005 - 2013) to inform the search for the star in the speckle years (1995 - 2005). When this new analysis technique is combined with the first complete re-reduction of Keck Galactic Center speckle images using speckle holography, we are able to track the short-period star S0-38 (K-band magnitude = 17, orbital period = 19 years) through the speckle years. We use the kinematic measurements from speckle holography and adaptive optics to estimate the orbits of S0-38 and S0-2 and thereby improve our constraints of the mass ($M_{bh}$) and distance ($R_o$) of Sgr A*: $M_{bh} = 4.02\pm0.16\pm0.04\times10^6~M_{\odot}$ and $7.86\pm0.14\pm0.04$ kpc. The uncertainties in $M_{bh}$ and $R_o$ as determined by the combined orbital fit of S0-2 and S0-38 are improved by a factor of 2 and 2.5, respectively, compared to an orbital fit of S0-2 alone and a factor of $\sim$2.5 compared to previous results from stellar orbits. This analysis also limits the extended dark mass within 0.01 pc to less than $0.13\times10^{6}~M_{\odot}$ at 99.7% confidence, a factor of 3 lower compared to prior work.

The Rapid Transient Surveyor (RTS) is a proposed rapid-response, high-cadence adaptive optics (AO) facility for the UH 2.2-m telescope on Maunakea. RTS will uniquely address the need for high-acuity and sensitive near-infrared spectral follow-up observations of tens of thousands of objects in mere months by combining an excellent observing site, unmatched robotic observational efficiency, and an AO system that significantly increases both sensitivity and spatial resolving power. We will initially use RTS to obtain the infrared spectra of ~4,000 Type Ia supernovae identified by the Asteroid Terrestrial-Impact Last Alert System over a two year period that will be crucial to precisely measuring distances and mapping the distribution of dark matter in the z < 0.1 universe. RTS will comprise an upgraded version of the Robo-AO laser AO system and will respond quickly to target-of-opportunity events, minimizing the time between discovery and characterization. RTS will acquire simultaneous-multicolor images with an acuity of 0.07-0.10" across the entire visible spectrum (20% i'-band Strehl in median conditions) and <0.16" in the near infrared, and will detect companions at 0.5" at contrast ratio of ~500. The system will include a high-efficiency prism integral field unit spectrograph: R = 70-140 over a total bandpass of 840-1830 nm with an 8.7" by 6.0" field of view (0.15" spaxels). The AO correction boosts the infrared point-source sensitivity of the spectrograph against the sky background by a factor of seven for faint targets, giving the UH 2.2-m the H-band sensitivity of a 5.7-m telescope without AO.

A cosmology of Poincare gauge theory is developed, where several properties of universe corresponding to the cosmological equations with the pseudoscalar torsion function are investigated. The cosmological constant is found to be the intrinsic torsion and curvature of the vacuum universe and is derived from the theory naturally rather than added artificially, i.e. the dark energy originates from geometry and includes the cosmological constant but differs from it. The cosmological constant puzzle, the coincidence and fine tuning problem are relieved naturally at the same time. By solving the cosmological equations, the analytic cosmological solution is obtained and can be compared with the LCDM model. In addition, the expressions of density parameters of the matter and the geometric dark energy are derived, and it is shown that the evolution of state equations for the geometric dark energy agrees with the current observational data. At last, the full equations of linear cosmological perturbations and the solutions are obtained.

Samples of early-type galaxies show a correlation between stellar velocity dispersion and the stellar initial mass function (IMF) as inferred from gravity-sensitive absorption lines in the galaxies' central regions. To search for spatial variations in the IMF, we have observed two early-type galaxies with Keck/LRIS and measured radial gradients in the strengths of absorption features from 4000-5500 \AA $\,$ and 8000-10,000 \AA. We present spatially resolved measurements of the dwarf-sensitive spectral indices NaI (8190 \AA) and Wing-Ford FeH (9915 \AA), as well as indices for species of H, C$_2$, CN, Mg, Ca, TiO, and Fe. Our measurements show a metallicity gradient in both objects, and Mg/Fe consistent with a shallow gradient in \alpha-enhancement, matching widely observed trends for massive early-type galaxies. The NaI index and the CN$_1$ index at 4160 \AA $\,$ exhibit significantly steeper gradients, with a break at $r \sim 0.1 r_{\rm eff}$ ($r \sim 300$ pc). Inside this radius NaI strength increases sharply toward the galaxy center, consistent with a rapid central rise in [Na/Fe]. In contrast, the ratio of FeH to Fe index strength decreases toward the galaxy center. This behavior cannot be reproduced by a steepening IMF inside $0.1 r_{\rm eff}$ if the IMF is a single power law. While gradients in the mass function above $\sim 0.4 M_\odot$ may occur, exceptional care is required to disentangle these IMF variations from the extreme variations in individual element abundances near the galaxies' centers.

The pulse profile of PSR B1133+16 is usually regarded as a conal-double structure. However, its multifrequency profiles cannot simply be fitted with two Gaussian functions, and a third component is always needed to fit the bridge region (between two peaks). This would introduce additional, redundant parameters. In this paper, through a comparison of five fitting functions (Gaussian, von Mises, hyperbolic secant, square hyperbolic secant, and Lorentz), it is found that the square hyperbolic secant function can best reproduce the profile, yielding an improved fit. Moreover, a symmetric 2D radiation beam function, instead of a simple 1D Gaussian function, is used to fit the profile. Each profile with either well-resolved or not-so-well-resolved peaks could be fitted adequately using this beam function, and the bridge emission between the two peaks does not need to be a new component. Adopting inclination and impact angles based on polarization measurements, the opening angle ({\theta}_{\mu}0) of the radiation beam in a certain frequency band is derived from beam-function fitting. The corresponding radiation altitudes are then calculated. Based on multi-frequency profiles, we also computed the Lorentz factors of the particles and their dispersion at those locations in both the curvature-radiation (CR) and inverse-Compton-scattering (ICS) models. We found that the Lorentz factors of the particles decrease rapidly as the radiation altitude increases. Besides, the radiation prefers to be generated in annular region rather than core region, and this needs further validation.

A neutrino-dominated accretion disk around a stellar-mass black hole (BH) can power a gamma-ray burst (GRB) via annihilation of neutrinos launched from the disk. For the BH hyperaccretion system, high accretion rate should trigger the violent evolution of the BH's characteristics, which further leads to the evolution of the neutrino annihilation luminosity. In this paper, we consider the evolution of the accretion system to analyze the mean time-dependent neutrino annihilation luminosity with the different mean accretion rates and initial BH parameters. By time-integrating the luminosity, the total neutrino annihilation energy with the reasonable initial disk mass can satisfy the most of short-duration GRBs and about half of long-duration GRBs. Moreover, the extreme Kerr BH should exist in the cental engines of some high-luminosity GRBs. GRBs with higher energy have to request the alternative magnetohydrodynamics processes in the centers, such as the Blandford-Znajek jet from the accretion system or the millisecond magnetar.

TThe glitch of anomalous X-ray pulsars \& soft gamma repeaters (AXP/SGRs) usually accompanied with detectable energy releases manifesting as X-ray bursts or outbursts, while the glitch of some pulsars like Vela release negligible energy. We find that these two types of glitches can naturally correspond to two types of starquake of solid strange stars. By applying the EoS of quark cluster star and some realistic pulsar parameters, we can reproduce consistent results compared with previous constraints and observations.

Terahertz imaging systems have received substantial attention from the scientific community for their use in astronomy, spectroscopy, plasma diagnostics and security. One approach to designing such systems is to use focal plane arrays. Although the principle of these systems is straightforward, realizing practical architectures has proven deceptively difficult. A different approach to imaging consists of spatially encoding the incoming flux of electromagnetic energy prior to detection using a reconfigurable mask. This technique is referred to as coded aperture or Hadamard imaging. This paper details the design, fabrication and testing of a prototype coded aperture mask operating at WR 1.5 (500 to 750 GHz) that uses the switching properties of vanadium dioxide (VO2). The reconfigurable mask consists of bowtie antennas with vanadium dioxide VO2 elements at the feed points. From the symmetry, a unit cell of the array can be represented by an equivalent waveguide whose dimensions limit the maximum operating frequency. In this design, the cutoff frequency of the unit cell is 640 GHz. The VO2 devices are grown using reactive-biased target ion beam deposition. A reflection coefficient (S11) measurement of the mask in the WR 1.5 (500 to 750 GHz) band is conducted. The results are compared with circuit models and found to be in good agreement. A simulation of the transmission response of the mask is conducted and shows a transmission modulation of up to 28 dB. This project is a first step towards the development of a full coded aperture imaging system operating at WR 1.5 with VO2 as the mask switching element.

At a projected distance of ~26 pc from Sgr A*, the Arches cluster provides insight to star formation in the extreme Galactic Center (GC) environment. Despite its importance, many key properties such as the cluster's internal structure and orbital history are not well known. We present an astrometric and photometric study of the outer region of the Arches cluster (R > 6.25") using HST WFC3IR. Using proper motions we calculate membership probabilities for stars down to F153M = 20 mag (~2.5 M_sun) over a 120" x 120" field of view, an area 144 times larger than previous astrometric studies of the cluster. We construct the radial profile of the Arches to a radius of 75" (~3 pc at 8 kpc), which can be well described by a single power law. From this profile we place a 3-sigma lower limit of 2.8 pc on the observed tidal radius, which is larger than the predicted tidal radius (1 - 2.5 pc). Evidence of mass segregation is observed throughout the cluster and no tidal tail structures are apparent along the orbital path. The absence of breaks in the profile suggests that the Arches has not likely experienced its closest approach to the GC between ~0.2 - 1 Myr ago. If accurate, this constraint indicates that the cluster is on a prograde orbit and is located front of the sky plane that intersects Sgr A*. However, further simulations of clusters in the GC potential are required to interpret the observed profile with more confidence.

Observations indicate that most of the universal matter are invisible and the gravitational constant $G(t)$ maybe depends on the time. A theory of the variational $G$ (VG) is explored in this paper, with naturally producing the useful dark components in universe. We utilize the observational data: lookback time data, model-independent gamma ray bursts, growth function of matter linear perturbations, type Ia supernovae data with systematic errors, CMB and BAO to restrict the unified model (UM) of dark components in VG theory. Using the best-fit values of parameters with the covariance matrix, constraints on the variation of $G$ are $(\frac{G}{G_{0}})_{z=3.5}\simeq 1.0015^{+0.0071}_{-0.0075}$ and $(\frac{\dot{G}}{G})_{today}\simeq -0.7252^{+2.3645}_{-2.3645}\times 10^{-13} yr^{-1}$, the small uncertainties around constants. Limit on the equation of state of dark matter is $w_{0dm}=0.0072^{+0.0170}_{-0.0170}$ with assuming $w_{0de}=-1$ in unified model, and dark energy is $w_{0de}=-0.9986^{+0.0011}_{-0.0011}$ with assuming $w_{0dm}=0$ at prior. Restriction on UM parameters are $B_{s}=0.7442^{+0.0137+0.0262}_{-0.0132-0.0292}$ and $\alpha=0.0002^{+0.0206+0.0441}_{-0.0209-0.0422}$ with $1\sigma$ and $2\sigma$ confidence level. In addition, the effect of a cosmic string fluid on unified model in VG theory are investigated. In this case it is found that the $\Lambda$CDM ($\Omega_{s}=0$, $\beta=0$ and $\alpha=0$) is included in this VG-UM model at $1\sigma$ confidence level, and the larger errors are given: $\Omega_{s}=-0.0106^{+0.0312+0.0582}_{-0.0305-0.0509}$ (dimensionless energy density of cosmic string), $(\frac{G}{G_{0}})_{z=3.5}\simeq 1.0008^{+0.0620}_{-0.0584}$ and $(\frac{\dot{G}}{G})_{today}\simeq -0.3496^{+26.3135}_{-26.3135}\times 10^{-13}yr^{-1}$.

Cheung identified a sample of 100 candidate X-shaped radio galaxies using the NRAO FIRST survey; these are small-axial-ratio extended radio sources with off-axis emission. Here we present radio images of 52 of these sources that have been made from archival Very Large Array data with resolution of about 1 arcsec. Fifty-one of the 52 were observed at 1.4 GHz, seven were observed at 1.4 GHz and 5 GHz, and one was observed only at 5 GHz. We also present overlays of the SDSS red images for 48 of the sources, and DSS II overlays for the remainder. Optical counterparts have been identified for most sources, but there remain a few empty fields. Our higher resolution VLA images along with FIRST survey images of the sources in the sample reveal that extended extragalactic radio sources with small axial ratios are largely (60%) cases of double radio sources with twin lobes that have off-axis extensions, usually with inversion-symmetric structure. The available radio images indicate that at most 20% sources might be genuine X-shaped radio sources that could have formed by a restarting of beams in a new direction following an interruption and axis flip. The remaining 20% are in neither of these categories. The implications of this result for the gravitational wave background are discussed in Roberts, Saripalli, and Subrahmanyan.

We present a metallicity analysis of 83 late-type giants within the central 1 pc of the Milky Way. K-band spectroscopy of these stars were obtained with the medium-spectral resolution integral-field spectrograph NIFS on Gemini North using laser-guide star adaptive optics. Using spectral template fitting with the MARCS synthetic spectral grid, we find that there is large variation in metallicity, with stars ranging from [M/H] $<$ -1.0 to above solar metallicity. About 6\% of the stars have [M/H] $<$ -0.5. This result is in contrast to previous observations, with smaller samples, that show stars at the Galactic center have approximately solar metallicity with only small variations. Our current measurement uncertainties are dominated by systematics in the model, especially at [M/H] $>$ 0, where there are stellar lines not represented in the model. However, the conclusion that there are low metallicity stars, as well as large variations in metallicity is robust. The metallicity may be an indicator of the origin of these stars. The low-metallicity population is consistent with that of globular clusters in the Milky Way, but their small fraction likely means that globular cluster infall is not the dominant mechanism for forming the Milky Way nuclear star cluster. The majority of stars are at or above solar metallicity, which suggests they were formed closer to the Galactic center or from the disk. In addition, our results indicate that it will be important for star formation history analyses using red giants at the Galactic center to consider the effect of varying metallicity.

We present a high spatial resolution optical and infrared study of the circumnuclear region in Arp 220, a late-stage galaxy merger. Narrowband imaging using HST/WFC3 has resolved the previously observed peak in H$\alpha$+[NII] emission into a bubble-shaped feature. This feature measures 1.6" in diameter, or 600 pc, and is only 1" northwest of the western nucleus. The bubble is aligned with the western nucleus and the large-scale outflow axis seen in X-rays. We explore several possibilities for the bubble origin, including a jet or outflow from a hidden active galactic nucleus (AGN), outflows from high levels of star formation within the few hundred pc nuclear gas disk, or an ultraluminous X-ray source. An obscured AGN or high levels of star formation within the inner $\sim$100 pc of the nuclei are favored based on the alignment of the bubble and energetics arguments.

The TMT Detailed Science Case describes the transformational science that the Thirty Meter Telescope will enable. Planned to begin science operations in 2024, TMT will open up opportunities for revolutionary discoveries in essentially every field of astronomy, astrophysics and cosmology, seeing much fainter objects much more clearly than existing telescopes. Per this capability, TMT's science agenda fills all of space and time, from nearby comets and asteroids, to exoplanets, to the most distant galaxies, and all the way back to the very first sources of light in the Universe. More than 150 astronomers from within the TMT partnership and beyond offered input in compiling the new 2015 Detailed Science Case. The contributing astronomers represent the entire TMT partnership, including the California Institute of Technology (Caltech), the Indian Institute of Astrophysics (IIA), the National Astronomical Observatories of the Chinese Academy of Sciences (NAOC), the National Astronomical Observatory of Japan (NAOJ), the University of California, the Association of Canadian Universities for Research in Astronomy (ACURA) and US associate partner, the Association of Universities for Research in Astronomy (AURA).

By proposing a pure leptonic radiation model, we study the potential gamma-ray emissions from jets of the low-mass X-ray binaries. In this model, the relativistic electrons that are accelerated in the jets are responsible for radiative outputs. Nevertheless, dynamics of jets are dominated by the magnetic and proton-matter kinetic energies. The model involves all kinds of related radiative processes and considers the evolution of relativistic electrons along the jet by numerically solving the kinetic equation. Numerical results show that the spectral energy distributions can extend up to TeV bands, in which synchrotron radiation and synchrotron self-Compton scattering are dominant components. As an example, we apply the model to the low-mass X-ray binary GX 339-4. The results can not only reproduce the currently available observations from GX 339-4, but also predict detectable radiation at GeV and TeV bands by Fermi and CTA telescopes. The future observations with Fermi and CTA can be used to test our model, which could be employed to distinguish the origin of X-ray emissions.

We present new near-infrared spectroscopic observations of the outer edges of the young stellar cluster around the supermassive black hole at the Galactic center. The observations show a break in the surface-density profile of young stars at approximately 13 arcsec (0.52 pc). These observations spectroscopically confirm previous suggestions of a break based on photometry. Using Gemini North's Near-Infrared Integral Field Spectrometer (NIFS) we are able to detect and separate early- and late-type stars with a 75% completeness at Ks = 15.5. We sample a region with radii between 7" to 23" (0.28 pc to 0.92 pc) from Sgr A*, and present new spectral classifications of 144 stars brighter than Ks = 15.5, where 140 stars are late-type (> 1 Gyr) and only four stars are early-type (young, 4-6 Myr). A broken power-law fit of the early-type surface-density matches well with our data and previously published values. The projected surface-density of late-type stars is also measured and found to be consistent with previous results. We find that the observed early-type surface-density profile is inconsistent with the theory of the young stars originating from a tightly bound infalling cluster, as no significant trail of young stars is found at radii above 13". We also note that either a simple disk instability criterion or a cloud-cloud collision could explain the location of the outer edge, though we lack information to make conclusive remarks on either alternative. If this break in surface-density represents an edge to the young stellar cluster it would set an important scale for the most recent episode of star formation at the Galactic center.

In this paper, we study cosmologies containing the reduced modified Chaplygin gas (RMCG) fluid which is reduced from the modified Chaplygin gas $p=A\rho-B\rho^{-\alpha}$ for the value of $\alpha=-1/2$. In this special case, dark cosmological models can be realized for different values of model parameter $A$. We investigate the viabilities of these dark cosmological models by discussing the evolutions of cosmological quantities and using the currently available cosmic observations. It is shown that the special RMCG model ($A=0$ or $A=1$) which unifies the dark matter and dark energy should be abandoned. For $A=1/3$, RMCG which unifies the dark energy and dark radiation is the favorite model according to the objective Akaike information criteria. In the case of $A<0$, RMCG can achieve the features of the dynamical quintessence and phantom models, where the evolution of the universe is not sensitive to the variation of model parameters.

We investigate the circumstellar disc fraction as determined from L-band excess observations of the young, massive Arches and Quintuplet clusters residing in the central molecular zone of the Milky Way. The Quintuplet cluster was searched for L-band excess sources for the first time. We find a total of 26 excess sources in the Quintuplet cluster and 21 in the Arches cluster, of which 13 are new detections. With the aid of proper motion membership samples, the disc fraction of the Quintuplet cluster was derived for the first time to be 4.0 +/- 0.7%. There is no evidence for a radially varying disc fraction in this cluster. In the case of the Arches cluster, a disc fraction of 9.2 +/- 1.2% approximately out to the cluster's predicted tidal radius, r < 1.5 pc, is observed. This excess fraction is consistent with our previously found disc fraction in the cluster in the radial range 0.3 < r < 0.8 pc. In both clusters, the host star mass range covers late A- to early B-type stars, 2 < M < 15 Msun, as derived from J-band photospheric magnitudes. We discuss the unexpected finding of dusty circumstellar discs in these UV intense environments in the context of primordial disc survival and formation scenarios of secondary discs. We consider the possibility that the L-band excess sources in the Arches and Quintuplet clusters could be the high-mass counterparts to T Tauri pre-transitional discs. As such a scenario requires a long pre-transitional disc lifetime in a UV intense environment, we suggest that mass transfer discs in binary systems are a likely formation mechanism for the B-star discs observed in these starburst clusters.

We numerically investigate the thermally unstable accretion disks around black holes. We adopt an evolutionary viscous stress equation to replace the standard alpha-prescription based on the results of two MHD simulations. We find a kind of interesting oscillations on some running models in limit-cycle outburst state. The oscillations arise near the inner boundary and propagate radially outwards. We deem that they are the trapped $p$-mode oscillations excited by sonic-point instability. We directly integrate the local radiation cooling fluxes to construct the mimic bolometric light-curve. We find a series of overtones beside the fundamental harmonic on the power spectra of mimic light-curves. The frequency of the fundamental harmonic is very close to the maximum epicyclic frequency of the disk and the frequency ratio of the fundamental harmonic and overtones is a regular integer series. We suggest that the code for ray-tracing calculation must be time-dependent in virtual observation and point out the robustness of the black hole spin measurement with high frequency QPOs.

We investigate the thermal stability of optically thin, two-temperature, radiative cooling-dominated accretion disks. Our linear analysis shows that the disk is thermally unstable without magnetic fields, which agrees with previous stability analysis on the Shapiro-Lightman-Eardley disk. By taking into account the effects of magnetic fields, however, we find that the disk can be or partly be thermally stable. Our results may be helpful to understand the outflows in optically thin flows. Moreover, such radiative cooling-dominated disks may provide a new explanation of the different behaviors between black hole and neutron star X-ray binaries on the radio/X-ray correlation.

In order to understand the main observational characteristics of the Galactic X-ray binary Cygnus X-3, we propose a radiation model in which high-energy electrons accelerated in the dissipation zone of a jet produce the non-thermal broadband emissions. Broadband spectral energy distributions are computed to confront with the AGILE and\emph{ Fermi} LAT data together with the multi-band data, during soft X-ray spectral states. By fitting observations at different locations of the jet, we find that the emission region is rather compact and should be located at the distance of about one orbital radius. Our results can explain the current multi-frequency observations and also predict the TeV band emission. The model could be tested by a polarization measurement at IR band, and/or by a correlation study between the GeV and TeV bands once very high-energy observations are available.

We report new observations of the Galactic Center source G2 from the W. M. Keck Observatory. G2 is a dusty red object associated with gas that shows tidal interactions as it nears closest approach with the Galaxy's central black hole. Our observations, conducted as G2 passed through periapse, were designed to test the proposal that G2 is a 3 earth mass gas cloud. Such a cloud should be tidally disrupted during periapse passage. The data were obtained using the Keck II laser guide star adaptive optics system (LGSAO) and the facility near-infrared camera (NIRC2) through the K' [2.1 $\mu$m] and L' [3.8 $\mu$m] broadband filters. Several results emerge from these observations: 1) G2 has survived its closest approach to the black hole as a compact, unresolved source at L'; 2) G2's L' brightness measurements are consistent with those over the last decade; 3) G2's motion continues to be consistent with a Keplerian model. These results rule out G2 as a pure gas cloud and imply that G2 has a central star. This star has a luminosity of $\sim$30 $L_{\odot} $ and is surrounded by a large ($\sim$2.6 AU) optically thick dust shell. The differences between the L' and Br-$\gamma$ observations can be understood with a model in which L' and Br-$\gamma$ emission arises primarily from internal and external heating, respectively. We suggest that G2 is a binary star merger product and will ultimately appear similar to the B-stars that are tightly clustered around the black hole (the so-called S-star cluster).

Linear properties of kinetic Alfv\'{e}n waves (KAWs) and kinetic slow waves (KSWs) are studied in the framework of two-fluid magnetohydrodynamics. We obtain the wave dispersion relations that are valid in a wide range of the wave frequency $\omega $ and plasma-to-magnetic pressure ratio $\beta $. The KAW frequency can reach and exceed the ion cyclotron frequency at ion kinetic scales, whereas the KSW frequency remains sub-cyclotron. At $\beta \sim 1$, the plasma and magnetic pressure perturbations of both modes are in anti-phase, so that there is nearly no total pressure perturbations. However, these modes exhibit also several opposite properties. At high $% \beta $, the electric polarization ratios of KAWs and KSWs are opposite at the ion gyroradius scale, where KAWs are polarized in sense of electron gyration (right-hand polarized) and KSWs are left-hand polarized. The magnetic helicity $\sigma \sim 1$ for KAWs and $\sigma \sim -1$ for KSWs, and the ion Alfv\'{e}n ratio $R_{Ai}\ll 1$ for KAWs and $R_{Ai}\gg 1$ for KSWs. We also found transition wavenumbers where KAWs change their polarization from left- to right-hand. These new properties can be used to discriminate KAWs and KSWs when interpreting kinetic-scale electromagnetic fluctuations observed in various solar-terrestrial plasmas. This concerns, in particular, identification of modes responsible for kinetic-scale pressure-balanced fluctuations and turbulence in the solar wind.

The Gemini Multi-conjugate adaptive optics System (GeMS) is a facility instrument for the Gemini-South telescope. It delivers uniform, near-diffraction-limited image quality at near-infrared wavelengths over a 2 arcminute field of view. Together with the Gemini South Adaptive Optics Imager (GSAOI), a near-infrared wide field camera, GeMS/GSAOI's combination of high spatial resolution and a large field of view will make it a premier facility for precision astrometry. Potential astrometric science cases cover a broad range of topics including exo-planets, star formation, stellar evolution, star clusters, nearby galaxies, black holes and neutron stars, and the Galactic center. In this paper, we assess the astrometric performance and limitations of GeMS/GSAOI. In particular, we analyze deep, mono-epoch images, multi-epoch data and distortion calibration. We find that for single-epoch, un-dithered data, an astrometric error below 0.2 mas can be achieved for exposure times exceeding one minute, provided enough stars are available to remove high-order distortions. We show however that such performance is not reproducible for multi-epoch observations, and an additional systematic error of ~0.4 mas is evidenced. This systematic multi-epoch error is the dominant error term in the GeMS/GSAOI astrometric error budget, and it is thought to be due to time-variable distortion induced by gravity flexure.

The High Energy cosmic-Radiation Detection (HERD) facility onboard China's Space Station is planned for operation starting around 2020 for about 10 years. It is designed as a next generation space facility focused on indirect dark matter search, precise cosmic ray spectrum and composition measurements up to the knee energy, and high energy gamma-ray monitoring and survey. The calorimeter plays an essential role in the main scientific objectives of HERD. A 3-D cubic calorimeter filled with high granularity crystals as active material is a very promising choice for the calorimeter. HERD is mainly composed of a 3-D calorimeter (CALO) surrounded by silicon trackers (TK) from all five sides except the bottom. CALO is made of 9261 cubes of LYSO crystals, corresponding to about 55 radiation lengths and 3 nuclear interaction lengths, respectively. Here the simulation results of the performance of CALO with GEANT4 and FLUKA are presented: 1) the total absorption CALO and its absorption depth for precise energy measurements (energy resolution: 1\% for electrons and gamma-rays beyond 100 GeV, 20\% for protons from 100 GeV to 1 PeV); 2) its granularity for particle identification (electron/proton separation power better than $10^{-5}$); 3) the homogenous geometry for detecting particles arriving from every unblocked direction for large effective geometrical factor ($>$3 ${\rm m}^{2}{\rm sr}$ for electron and diffuse gamma-rays, $>$2 $ {\rm m}^{2}{\rm sr}$ for cosmic ray nuclei); 4) expected observational results such as gamma-ray line spectrum from dark matter annihilation and spectrum measurement of various cosmic ray chemical components.

IRIS (InfraRed Imaging Spectrograph) is a first light near-infrared diffraction limited imager and integral field spectrograph being designed for the future Thirty Meter Telescope (TMT). IRIS is optimized to perform astronomical studies across a significant fraction of cosmic time, from our Solar System to distant newly formed galaxies (Barton et al. [1]). We present a selection of the innovative science cases that are unique to IRIS in the era of upcoming space and ground-based telescopes. We focus on integral field spectroscopy of directly imaged exoplanet atmospheres, probing fundamental physics in the Galactic Center, measuring 10^4 to 10^10 Msun supermassive black hole masses, resolved spectroscopy of young star-forming galaxies (1 < z < 5) and first light galaxies (6 < z < 12), and resolved spectroscopy of strong gravitational lensed sources to measure dark matter substructure. For each of these science cases we use the IRIS simulator (Wright et al. [2], Do et al. [3]) to explore IRIS capabilities. To highlight the unique IRIS capabilities, we also update the point and resolved source sensitivities for the integral field spectrograph (IFS) in all five broadband filters (Z, Y, J, H, K) for the finest spatial scale of 0.004" per spaxel. We briefly discuss future development plans for the data reduction pipeline and quicklook software for the IRIS instrument suite.

Glitch (sudden spinup) is a common phenomenon in pulsar observations. However, the physical mechanism of glitch is still a matter of debate because it depends on the puzzle of pulsar's inner structure, i.e., the equation of state of dense matter. Some pulsars (e.g., Vela-like) show large glitches ({\Delta}{\nu}/{\nu}~10^-6) but release negligible energy, whereas the large glitches of AXPs/SGRs (anomalous X-ray pulsars/soft gamma repeaters) are usually (but not always) accompanied with detectable energy releases manifesting as X-ray bursts or outbursts. We try to understand this aspect of glitches in a starquake model of solid quark stars. There are actually two kinds of glitches in this scenario: bulk-invariable (Type I) and bulk-variable (Type II) ones. The total stellar volume changes (and then energy releases) significantly for the latter but not for the former. Therefore, glitches accompanied with X-ray bursts (e.g., that of AXP/SGRs) could originate from Type II starquakes induced probably by accretion, while the others without evident energy release (e.g., that of Vela pulsar) would be the result of Type I starquakes due to, simply, a change of stellar ellipticity.

We investigate the orbital motion of the Quintuplet cluster near the Galactic center with the aim of constraining formation scenarios of young, massive star clusters in nuclear environments. Three epochs of adaptive optics high-angular resolution imaging with Keck/NIRC2 and VLT/NACO were obtained over a time baseline of 5.8 years, delivering an astrometric accuracy of 0.5-1 mas/yr. Proper motions were derived in the cluster reference frame and were used to distinguish cluster members from the majority of field stars. Fitting the cluster and field proper motion distributions with 2D gaussian models, we derive the orbital motion of the cluster for the first time. The Quintuplet is moving with a 2D velocity of 132 +/- 15 km/s with respect to the field along the Galactic plane, which yields a 3D orbital velocity of 167 +/- 15 km/s when combined with the previously known radial velocity. From a sample of 119 stars measured in three epochs, we derive an upper limit to the velocity dispersion in the core of the Quintuplet cluster of sigma_1D < 10 km/s. Knowledge of the three velocity components of the Quintuplet allows us to model the cluster orbit in the potential of the inner Galaxy. Comparing the Quintuplet's orbit with the Arches orbit, we discuss the possibility that both clusters originated in the same area of the central molecular zone. [abridged]

We present the effects of the self-gravity on the vertical structure and neutrino luminosity of the neutrino-dominated accretion disks in cylindrical coordinates. It is found that significant changes of the structure appear in the outer region of the disk, especially for high accretion rates (e.g., $\ga 1 M_\odot~\rm s^{-1}$), and thus causes the slight increase of the neutrino luminosity. Furthermore, the gravitational instability of the disk is reviewed by the vertical distribution of the Toomre parameter \citep{Toomre1964}, which may account for the late-time flares in gamma-ray bursts and the extended emission in short-duration gamma-ray bursts.

The prototype Robo-AO system at the Palomar Observatory 1.5-m telescope is the world's first fully automated laser adaptive optics instrument. Scientific operations commenced in June 2012 and more than 12,000 observations have since been performed at the ~0.12" visible-light diffraction limit. Two new infrared cameras providing high-speed tip-tilt sensing and a 2' field-of-view will be integrated in 2014. In addition to a Robo-AO clone for the 2-m IGO and the natural guide star variant KAPAO at the 1-m Table Mountain telescope, a second generation of facility-class Robo-AO systems are in development for the 2.2-m University of Hawai'i and 3-m IRTF telescopes which will provide higher Strehl ratios, sharper imaging, ~0.07", and correction to {\lambda} = 400 nm.

We study the origin of non-thermal emissions from the Galactic black hole X-ray binary Cygnus X--1 which is a confirmed high mass microquasar. By analogy with those methods used in studies of active galactic nuclei, we propose a two-dimensional, time-dependent radiation model from the microquasar Cygnus X--1. In this model, the evolution equation for relativistic electrons in a conical jet are numerically solved by including escape, adiabatic and various radiative losses. The radiative processes involved are synchrotron emission, its self-Compton scattering, and inverse Compton scatterings of an accretion disk and its surrounding stellar companion. This model also includes an electromagnetic cascade process of an anisotropic $\gamma$-$\gamma$ interaction. We study the spectral properties of electron evolution and its emission spectral characteristic at different heights of the emission region located in the jet. We find that radio data from Cygnus X--1 are reproduced by the synchrotron emission, the \emph{Fermi} LAT measurements by the synchrotron emission and Comptonization of photons of the stellar companion, the TeV band emission fluxes by the Comptonization of the stellar photons. Our results show that: (1) Radio emission region extends from the binary system scales to the termination of the jet. (2) The GeV band emissions should originate from the distance close to the binary system scales. (3) The TeV band emissions could be inside the binary system, and these emissions could be probed by the upcoming CTA telescope. (4) The MeV tail emissions, which produces a strongly linearly polarized signal, are emitted inside the binary system. The location of the emissions is very close to the inner region of the jet.

The particular giant X-ray bump of GRB 121027A triggered by \emph{Swift} is quite different from the typical X-ray flares in gamma-ray bursts. There exhibit four parts of the observed structural variabilities in the rise and decay phase of the bump. Considering the quality of four parts of the data, we can only analyze the data from about 5300 s to about 6100 s in the bump using the stepwise filter correlation method (Gao et al. 2012), and find that the $86^{+5.9}_{-9.4}~\rm s$ periodic oscillation may exist, which is confirmed by the Lomb-Scargle method (Scargle 1982). Furthermore, a jet precession model (Liu et al. 2010) is proposed to account for such a variability.

GRB 120326A is an unusual gamma-ray burst (GRB) which has a quite long plateau and a very late rebrightening both in X-ray and optical bands. The similar behavior of the optical and X-ray light curves suggests that they maybe have a common origin. The long plateau starts from several hundred seconds and ends at tens of thousands seconds. The peak time of the late rebrightening is about 30000 s. We analyze the energy injection model by means of numerical and analytical solutions, considering both the wind environment and ISM environment for GRB afterglows. We especially study the influence of the injection starting time, ending time, stellar wind density (or density of the circumburst environment), and injection luminosity on the shape of the afterglow light curves, respectively. We find that the light curve is largely affected by the parameters in the wind model. There is a "bump" at the late time only in the wind model too. In the wind case, it is interesting that the longer the energy injected, the more obvious the rebrightening will be. We also find the peak time of bump is determined by the stellar wind density. We use the late continuous injection model to interpret the unusual afterglow of GRB 120326A. The model can well fit the observational data, however, we find that the time scale of the injection must be larger than ten thousands seconds. This implies that the time scale of the central engine activity must be more than ten thousands seconds. This can give useful constraints on the central engine of GRBs. We consider a new born millisecond pulsar with strong magnetic field as the central engine. On the other hand, our results suggest that the circumburst environment of GRB 120326A is very likely a stellar wind.

We present new kinematic measurements and modeling of a sample of 116 young stars in the central parsec of the Galaxy in order to investigate the properties of the young stellar disk. The measurements were derived from a combination of speckle and laser guide star adaptive optics imaging and integral field spectroscopy from the Keck telescopes. Compared to earlier disk studies, the most important kinematic measurement improvement is in the precision of the accelerations in the plane of the sky, which have a factor of six smaller uncertainties (~10 uas/yr/yr). We have also added the first radial velocity measurements for 8 young stars, increasing the sample at the largest radii (6"-12") by 25%. We derive the ensemble properties of the observed stars using Monte-Carlo simulations of mock data. There is one highly significant kinematic feature (~20 sigma), corresponding to the well-known clockwise disk, and no significant feature is detected at the location of the previously claimed counterclockwise disk. The true disk fraction is estimated to be ~20%, a factor of ~2.5 lower than previous claims, suggesting that we may be observing the remnant of what used to be a more densely populated stellar disk. The similarity in the kinematic properties of the B stars and the O/WR stars suggests a common star formation event. The intrinsic eccentricity distribution of the disk stars is unimodal, with an average value of <e> = 0.27 +/- 0.07, which we show can be achieved through dynamical relaxation in an initially circular disk with a moderately top-heavy mass function.

GRB 130925A, composed of three gamma-ray emission episodes and a series of orderly flares, has been detected by $Swift$, Fermi, Konus-$Wind$, and $INTEGRAL$. If the third weakest gamma-ray episode can be considered as a giant flare, we find that after the second gamma-ray episode observed by $INTEGRAL$ located at about 2000 s, there exists a positive relation between the time intervals of the adjacent flares and the time since the episode. We suggest that the second gamma-ray episode and its flares originate from the resumption of the accretion process due to the fragments from the collapsar falling back, such a relation may be related to a hyperaccretion disk around a precessed black hole (BH). We propose that the origin and time evolution of the flares, and the approximately symmetrical temporal structure and spectral evolution of the single flare, can be explained well by the jet precession model \citep{Liu2010}. In addition, the mass and spin of the BH can be constrained, which indicates a stellar-mass, fast rotating BH located in the center of GRB 130925A.

We give an update of the observations and analysis of G2 - the gaseous red emission-line object that is on a very eccentric orbit around the Galaxy's central black hole and predicted to come within 2400 Rs in early 2014. During 2013, the laser guide star adaptive optics systems on the W. M. Keck I and II telescopes were used to obtain three epochs of spectroscopy and imaging at the highest spatial resolution currently possible in the near-IR. The updated orbital solution derived from radial velocities in addition to Br-Gamma line astrometry is consistent with our earlier estimates. Strikingly, even ~6 months before pericenter passage there is no perceptible deviation from a Keplerian orbit. We furthermore show that a proposed "tail" of G2 is likely not associated with it but is rather an independent gas structure. We also show that G2 does not seem to be unique, since several red emission-line objects can be found in the central arcsecond. Taken together, it seems more likely that G2 is ultimately stellar in nature, although there is clearly gas associated with it.

We present 3D kinematic observations of stars within the central 0.5 pc of the Milky Way nuclear star cluster using adaptive optics imaging and spectroscopy from the Keck telescopes. Recent observations have shown that the cluster has a shallower surface density profile than expected for a dynamically relaxed cusp, leading to important implications for its formation and evolution. However, the true three dimensional profile of the cluster is unknown due to the difficulty in de-projecting the stellar number counts. Here, we use spherical Jeans modeling of individual proper motions and radial velocities to constrain for the first time, the de-projected spatial density profile, cluster velocity anisotropy, black hole mass ($M_\mathrm{BH}$), and distance to the Galactic center ($R_0$) simultaneously. We find that the inner stellar density profile of the late-type stars, $\rho(r)\propto r^{-\gamma}$ to have a power law slope $\gamma=0.05_{-0.60}^{+0.29}$, much more shallow than the frequently assumed Bahcall $\&$ Wolf slope of $\gamma=7/4$. The measured slope will significantly affect dynamical predictions involving the cluster, such as the dynamical friction time scale. The cluster core must be larger than 0.5 pc, which disfavors some scenarios for its origin. Our measurement of $M_\mathrm{BH}=5.76_{-1.26}^{+1.76}\times10^6$ $M_\odot$ and $R_0=8.92_{-0.55}^{+0.58}$ kpc is consistent with that derived from stellar orbits within 1$^{\prime\prime}$ of Sgr A*. When combined with the orbit of S0-2, the uncertainty on $R_0$ is reduced by 30% ($8.46_{-0.38}^{+0.42}$ kpc). We suggest that the MW NSC can be used in the future in combination with stellar orbits to significantly improve constraints on $R_0$.

The existence of intermediate-mass black holes (IMBHs) in globular clusters (GCs) remains a crucial problem. Searching IMBHs in GCs reveals a discrepancy between radio observations and dynamical modelings: the upper mass limits constrained by radio observations are systematically lower than that of dynamical modelings. One possibility for such a discrepancy is that, as we suggest in this work, there exist outflows in accretion flows. Our results indicate that, for most sources, current radio observations cannot rule out the possibility that IMBHs may exist in GCs. In addition, we adopt an \dot{M}-L_{R} relation to revisit this issue, which confirms the results obtained by the Fundamental Plane relation.

Rapid temporal variability has been widely observed in the light curves of gamma-ray bursts (GRBs). One possible mechanism for such variability is related to the relativistic eddies in the jet. In this paper, we include the contribution of the inter-eddy medium together with the eddies to the gamma-ray emission. We show that the gamma-ray emission can either lead or lag behind the observed synchrotron emission, where the latter originates in the inter-eddy medium and provides most of seed photons for producing gamma-ray emission through the inverse-Compton scattering. As a consequence, we argue that the lead/lag found in non-stationary short-lived light curves may not reveal the intrinsic lead/lag of different emission components. In addition, our results may explain the lead of gamma-ray emission with respect to optical emission observed in GRB 080319B.

The neutrino-cooled accretion disk, which was proposed to work as the central engine of gamma-ray bursts, encounters difficulty in interpreting the X-ray flares after the prompt gamma-ray emission. In this paper, the magnetic coupling between the inner disk and the central black hole is taken into consideration. For mass accretion rates around $0.001 \sim 0.1 M_{\sun} s^{-1}$, our results show that the luminosity of neutrino annihilation can be significantly enhanced due to the coupling effects. As a consequence, after the gamma-ray emission, a remnant disk with mass $M_{disk} \la 0.5 M_{\sun}$ may power most of the observed X-ray flares with the rest frame duration less than 100 seconds. In addition, a comparison between the magnetic coupling process and the Blandford-Znajek mechanism is shown on the extraction of black hole rotational energy.

Neutrino-dominated accretion flows (NDAFs) around rotating stellar-mass black holes are plausible candidates for the central engines of gamma-ray bursts (GRBs). We investigate one-dimensional global solutions of NDAFs, taking account of general relativity in Kerr metric, neutrino physics and nucleosynthesis more precisely than previous works. We calculate sixteen solutions with different characterized accretion rates and black hole spins to exhibit the radial distributions of various physical properties in NDAFs. \iffalse We find that the gas pressure and the neutrino cooling always become dominated in the inner region for large accretion rate, $\gtrsim 0.1 M_\odot$ $\rm s^{-1}$. \fi We confirm that the electron degeneracy has important effects in NDAFs and we find that the electron fraction is about 0.46 in the outer region for all the sixteen solutions. From the perspective of the mass fraction, free nucleons, $\rm ^4He$, and $\rm ^{56}Fe$ dominate in the inner, middle, and outer region, respectively. The influence of neutrino trapping on the annihilation is of importance for the superhigh accretion ($\dot M=10M_{\odot} \rm s^{-1}$) and most of the sixteen solutions have an adequate annihilation luminosity for GRBs.

We present new observations and analysis of G2 - the intriguing red emission-line object which is quickly approaching the Galaxy's central black hole. The observations were obtained with the laser guide star adaptive optics systems on the W. M. Keck I and II telescopes and include spectroscopy (R ~ 3600) centered on the Hydrogen Br-gamma line as well as K' (2.1 micrometer) and L' (3.8 micrometer) imaging. Analysis of these observations shows the Br-gamma line emission has a positional offset from the L' continuum. This offset is likely due to background source confusion at L'. We therefore present the first orbital solution derived from Br-gamma line astrometry, which when coupled with radial velocity measurements, results in a later time of closest approach (2014.21 +/- 0.14), closer periastron (130 AU, 1900Rs), and higher eccentricity (0.9814 +/- 0.0060) compared to a solution using L' astrometry. The new orbit casts doubt on previous associations of G2 and a low surface brightness "tail". It is shown that G2 has no K' counterpart down to K' ~ 20 mag. G2's L' continuum and the Br-gamma line-emission is unresolved in almost all epochs; however it is marginally extended in our highest quality Br-gamma data set from 2006 and exhibits a clear velocity gradient at that time. While the observations altogether suggest that G2 has a gaseous component which is tidally interacting with the central black hole, there is likely a central star providing the self-gravity necessary to sustain the compact nature of this object.

We have made the first detection of a near-infrared counterpart associated with the disk around Radio Source "I," a massive protostar in the Kleinmann-Low Nebula in Orion using imaging with laser guide star adaptive optics on the Keck II telescope. The infrared emission is evident in images acquired using L' (3.8 microns) and Ms (4.7 microns) filters and is not detectable at K' (2.1 microns). The observed morphology strongly suggests that we are seeing some combination of scattered and thermal light emanating from the disk. The disk is also manifest in the L'/Ms flux ratio image. We interpret the near-infrared emission as the illuminated surface of a nearly edge-on disk, oriented so that only the northern face is visible; the opposite surface remains hidden by the disk. We do not see infrared radiation associated directly with the star proposed to be associated with Source "I." The data also suggest that there is a cavity above and below the disk that is oriented perpendicular to the disk, and is sculpted by the known, strong outflow from the inner disk of Source I. We compare our data to models of a protostar with a surrounding disk, envelope, and wind-blown cavity in order to elucidate the nature of the disk around Radio Source I.

The nature of pulsar is still unknown because of non-perturbative effects of the fundamental strong interaction, and different models of pulsar inner structures are then suggested, either conventional neutron stars or quark stars. Additionally, a state of quark-cluster matter is conjectured for cold matter at supranuclear density, as a result pulsars could thus be quark-cluster stars. Besides understanding different manifestations, the most important issue is to find an effective way to observationally differentiate those models. X-ray polarimetry would play an important role here. In this letter, we focus on the thermal X-ray polarization of quark/quark-cluster stars. While the thermal X-ray linear polarization percentage is typically higher than ~10% in normal neutron star models, the percentage of quark/quark-cluster stars is almost zero. It could then be an effective method to identify quark/quark-cluster stars by soft X-ray polarimetry. We are therefore expecting to detect thermal X-ray polarization in the coming decades.

The supermassive black hole at the center of the Milky Way plays host to a massive, young cluster that may have formed in one of the most inhospitable environments in the Galaxy. We present new measurements of the global properties of this cluster, including the initial mass function (IMF), age, and cluster mass. These results are based on Keck laser-guide-star adaptive optics observations used to identify the young stars and measure their Kp-band luminosity function as presented in Do et al. 2013. A Bayesian inference methodology is developed to simultaneously fit the global properties of the cluster utilizing the observations and extensive simulations of synthetic star clusters. We find that the slope of the mass function for this cluster is \alpha = 1.7 +/- 0.2, which is steeper than previously reported, but still flatter than the traditional Salpeter slope of 2.35. The age of the cluster is between 2.5-5.8 Myr with 95% confidence, which is a younger age than typically adopted but consistent within the uncertainties of past measurements. The exact age of the cluster is difficult to determine since our results show two distinct age solutions (3.9 Myr and 2.8 Myr) due to model degeneracies in the relative number of Wolf-Rayet and OB stars. The total cluster mass is between 14,000 - 37,000 \msun above 1 \msun and it is necessary to include multiple star systems in order to fit the observed luminosity function and the number of observed Wolf-Rayet stars. The new IMF slope measurement is now consistent with X-ray observations indicating a factor of 10 fewer X-ray emitting pre-main-sequence stars than expected when compared with a Salpeter IMF. The young cluster at the Galactic center is one of the few definitive examples of an IMF that deviates significantly from the near-universal IMFs found in the solar neighborhood.

We present new high angular resolution near-infrared spectroscopic observations of the nuclear star cluster surrounding the Milky Way's central supermassive black hole. Using the integral-field spectrograph OSIRIS on Keck II behind the laser-guide-star adaptive optics system, this spectroscopic survey enables us to separate early-type (young, 4-6 Myr) and late-type (old, >1 Gyr) stars with a completeness of 50% down to K' = 15.5 mag, which corresponds to ~10 \msun for the early-type stars. This work increases the radial extent of reported OSIRIS/Keck measurements by more than a factor of 3 from 4" to 14" (0.16 pc to 0.56 pc), along the projected disk of young stars. For our analysis, we implement a new method of completeness correction using a combination of star-planting simulations and Bayesian inference. We assign probabilities for the spectral type of every source detected in deep imaging down to K' = 15.5 mag using information from spectra, simulations, number counts, and the distribution of stars. The inferred radial surface-density profiles, $\Sigma(R) \propto R^{-\Gamma}$, for the young stars and late-type giants are consistent with earlier results ($\Gamma_{early} = 0.93 \pm 0.09$, $\Gamma_{late} = 0.16 \pm 0.07$). The late-type surface-density profile is approximately flat out to the edge of the survey. While the late-type stellar luminosity function is consistent with the Galactic bulge, the completeness-corrected luminosity function of the early-type stars has significantly more young stars at faint magnitudes compared to previous surveys with similar depth. This luminosity function indicates that the corresponding mass function of the young stars is likely less top-heavy than that inferred from previous surveys.

We investigate the vertical structure and elements distribution of neutrino-dominated accretion flows around black holes in spherical coordinates with the reasonable nuclear statistical equilibrium. According our calculations, heavy nuclei tend to be produced in a thin region near the disk surface, whose mass fractions are primarily determined by the accretion rate and the vertical distribution of temperature and density. In this thin region, we find that $^{56}\rm Ni$ is dominant for the flow with low accretion rate (e.g., $0.05 $$M_{\odot}$ $\rm s^{-1}$) but $^{56}\rm Fe$ is dominant for the high counterpart (e.g., $1 M_{\odot}$ $\rm s^{-1}$). The dominant $^{56}\rm Ni$ in the special region may provide a clue to understand the bumps in the optical light curve of core-collapse supernovae.

We study the fluctuations of standard thin accretion disks by linear analysis of the time-dependent energy equation together with the vertical hydrostatic equilibrium and the equation of state. We show that some of the simulation results in Hirose et al. (2009b), such as the time delay, the relationship of power spectra, and the correlation between magnetic energy and radiation energy, can be well understood by our analytic results.

The system of accretion disk and black hole is usually considered as the central engine of Gamma-ray Bursts (GRBs). It is usually thought that the disk in the central engine of GRBs is the advection-dominated accretion disk, which is developed from a massive (mass M_disk) torus at radius r_disk. We find a positive correlation between the isotropic gamma-ray energy E_iso and duration (so-called T_90) for GRBs. We interpret this correlation within the advection-dominated accretion disk model, associating E_iso and T_90 with M_ disk and viscous timescale respectively.

Stars with short orbital periods at the center of our galaxy offer a powerful and unique probe of a supermassive black hole. Over the past 17 years, the W. M. Keck Observatory has been used to image the Galactic center at the highest angular resolution possible today. By adding to this data set and advancing methodologies, we have detected S0-102, a star orbiting our galaxy's supermassive black hole with a period of just 11.5 years. S0-102 doubles the number of stars with full phase coverage and periods less than 20 years. It thereby provides the opportunity with future measurements to resolve degeneracies in the parameters describing the central gravitational potential and to test Einstein's theory of General Relativity in an unexplored regime.

Adaptive Optics observations have dramatically improved the quality and versatility of high angular resolution measurements of the center of our Galaxy. In this paper, we quantify the quality of our Adaptive Optics observations and report on the astrometric precision for the young stellar population that appears to reside in a stellar disk structure in the central parsec. We show that with our improved astrometry and a 16 year baseline, including 10 years of speckle and 6 years of laser guide star AO imaging, we reliably detect accelerations in the plane of the sky as small as 70 microarcsec/yr/yr (~2.5 km/s/yr) and out to a projected radius from the supermassive black hole of 1.5" (~0.06 pc). With an increase in sensitivity to accelerations by a factor of ~6 over our previous efforts, we are able to directly probe the kinematic structure of the young stellar disk, which appears to have an inner radius of 0.8". We find that candidate disk members are on eccentric orbits, with a mean eccentricity of < e > = 0.30 +/- 0.07. Such eccentricities cannot be explained by the relaxation of a circular disk with a normal initial mass function, which suggests the existence of a top-heavy IMF or formation in an initially eccentric disk.

In this paper, we revisit generalized Chaplygin gas (GCG) model as a unified dark matter and dark energy model. The energy density of GCG model is given as $\rho_{GCG}/\rho_{GCG0}=[B_{s}+(1-B_{s})a^{-3(1+\alpha)}]^{1/(1+\alpha)}$, where $\alpha$ and $B_s$ are two model parameters which will be constrained by type Ia supernova as standard candles, baryon acoustic oscillation as standard rulers and the seventh year full WMAP data points. In this paper, we will not separate GCG into dark matter and dark energy parts any more as adopted in the literatures. By using Markov Chain Monte Carlo method, we find the result: $\alpha=0.00126_{- 0.00126- 0.00126}^{+ 0.000970+ 0.00268}$ and $B_s= 0.775_{- 0.0161- 0.0338}^{+ 0.0161+ 0.0307}$.

The physical nature of gamma-ray bursts (GRBs) are believed to involve an ultra-relativistic jet. The observed complex structure of light curves motivate the idea of jet precession. In this work, we study the gravitational waves of jet precession based on neutrino-dominated accretion disks around black holes, which may account for the central engine of GRBs. In our model, the jet and the inner part of the disk may precess along with the black hole, which is driven by the outer part of the disk. Gravitational waves are therefore expected to be significant from this black hole-inner disk precession system. By comparing our numerical results with the sensitivity of some detectors, we find that it is possible for DECIGO and BBO to detect such gravitational waves, particularly for GRBs in the Local Group.

A modified f(G) gravity model with coupling between matter and geometry is proposed, which is described by the product of the Lagrange density of the matter and an arbitrary function of the Gauss-Bonnet term. The field equations and the equations of motion corresponding to this model show the non-conservation of the energy-momentum tensor, the presence of an extra-force acting on test particles and the non-geodesic motion. Moreover, the energy conditions and the stability criterion at de Sitter point in the modified f(G) gravity models with curvature-matter coupling are derived, which can degenerate to the well-known energy conditions in general relativity. Furthermore, in order to get some insight on the meaning of these energy conditions, we apply them to the specific models of f(G) gravity and the corresponding constraints on the models are given. In addition, the conditions and the candidate for late-time cosmic accelerated expansion in the modified f(G) gravity are studied by means of conditions of power-law expansion and the equation of state of matter less than -1/ 3 .

We use the Markov ChainMonte Carlo method to investigate global constraints on the generalized holographic (GH) dark energy with flat and non-flat universe from the current observed data: the Union2 dataset of type supernovae Ia (SNIa), high-redshift Gamma-Ray Bursts (GRBs), the observational Hubble data (OHD), the cluster X-ray gas mass fraction, the baryon acoustic oscillation (BAO), and the cosmic microwave background (CMB) data. The most stringent constraints on the GH model parameter are obtained. In addition, it is found that the equation of state for this generalized holographic dark energy can cross over the phantom boundary wde =-1.

Time variable cosmological constant (TVCC) of holographic origin with interaction in Brans-Dicke theory is discussed in this paper. We investigate some characters for this model, and show the evolutions of deceleration parameter and equation of state (EOS) for dark energy. It is shown that in this scenario an accelerating universe can be obtained and the evolution of EOS for dark energy can cross over the boundary of phantom divide. In addition, a geometrical diagnostic method, jerk parameter is applied to this model to distinguish it with cosmological constant.

In this paper we investigate the interaction between dark matter and dark energy by considering two different interacting scenarios, i.e. the cases of constant interaction function and variable interaction function. By fitting the current observational data to constrain the interacting models, it is found that the interacting strength is non-vanishing, but weak for the case of constant interaction function, and the interaction is not obvious for the case of variable interaction function. In addition, for seeing the influence from interaction we also investigate the evolutions of interaction function, effective state parameter for dark energy and energy density of dark matter. At last some geometrical quantities in the interacting scenarios are discussed.

In this paper, a periodicity analysis of the radio light curves of the blazar NRAO 530 at 14.5, 8.0, and 4.8 GHz is presented employing an improved Phase Dispersion Minimization (PDM) technique. The result, which shows two persistent periodic components of $ \sim 6$ and $ \sim 10$ years at all three frequencies, is consistent with the results obtained with the Lomb-Scargle periodogram and weighted wavelet Z-transform algorithms. The reliability of the derived periodicities is confirmed by the Monte Carlo numerical simulations which show a high statistical confidence. (Quasi-)Periodic fluctuations of the radio luminosity of NRAO 530 might be associated with the oscillations of the accretion disk triggered by hydrodynamic instabilities of the accreted flow. \keywords{methods: statistical -- galaxies: active -- galaxies: quasar: individual: NRAO 530}

We report the first detection of the intrinsic velocity dispersion of the Arches cluster - a young (~2 Myr), massive (~10,000 Solar Mass) starburst cluster located near the Galactic center. This was accomplished using proper motion measurements within the central region of the cluster, obtained with the laser guide star adaptive optics system at Keck Observatory over a 3 year time baseline (2006-2009). This uniform dataset results in proper motion measurements that are improved by a factor ~5 over previous measurements from heterogeneous instruments, yielding internal velocity dispersion estimates 0.15 +/- 0.01 mas/yr, which corresponds to 5.4 +/- 0.4 km/s at a distance of 8.4 kpc. Projecting a simple model for the cluster onto the sky to compare with our proper motion dataset, in conjunction with surface density data, we estimate the total present-day mass of the cluster to be 15,000 (+7400 -6000) Solar masses. The mass in stars observed within a cylinder of radius R=0.4 pc is found to be 9000 (+4000 -3500) Solar Masses at formal 3-sigma confidence. This mass measurement is free from assumptions about the mass function of the cluster, and thus may be used to check mass estimates from photometry and simulation. When we conduct this check, we find that the present-day mass function of the Arches cluster is likely either top-heavy or truncated at low-mass, or both. Collateral benefits of our data and analysis include: 1. cluster membership probabilities, which may be used to extract a clean cluster sample for future photometric work; 2. a refined estimate of the bulk motion of the Arches cluster with respect to the field, which we find to be 172 +/- 15 km/s, which is slightly slower than suggested by previous VLT-Keck measurements; and 3. a velocity dispersion estimate for the field itself, which is likely dominated by the inner galactic bulge and the nuclear disk.

In accretion systems, outflows may have significant influence on the luminosity fluctuations. In this paper, following the Lyubarskii's general scheme, we revisit the power spectral density of luminosity fluctuations by taking into account the role of outflows. Our analysis is based on the assumption that the coupling between the local outflow and inflow is weak on the accretion rate fluctuations. We find that, for the inflow mass accretion rate $\dot M \propto r^{s}$, the power spectrum of flicker noise component will present a power-law distribution $p(f) \propto f^{-(1+4s/3)}$ for advection-dominated flows. We also obtain descriptions of $p(f)$ for both standard thin discs and neutrino-cooled discs, which show that the power-law index of a neutrino-cooled disc is generally larger than that of a photon-cooled disc. Furthermore, the obtained relationship between $p(f)$ and $s$ indicates the possibility of evaluating the strength of outflows by the power spectrum in X-ray binaries and gamma-ray bursts. In addition, we discuss the possible influence of the outflow-inflow coupling on our results.

We systematically study black holes in the Horava-Lifshitz (HL) theory by following the kinematic approach, in which a horizon is defined as the surface at which massless test particles are infinitely redshifted. Because of the nonrelativistic dispersion relations, the speed of light is unlimited, and test particles do not follow geodesics. As a result, there are significant differences in causal structures and black holes between general relativity (GR) and the HL theory. In particular, the horizon radii generically depend on the energies of test particles. Applying them to the spherical static vacuum solutions found recently in the nonrelativistic general covariant theory of gravity, we find that, for test particles with sufficiently high energy, the radius of the horizon can be made as small as desired, although the singularities can be seen in principle only by observers with infinitely high energy. In these studies, we pay particular attention to the global structure of the solutions, and find that, because of the foliation-preserving-diffeomorphism symmetry, ${Diff}(M,{\cal{F}})$, they are quite different from the corresponding ones given in GR, even though the solutions are the same. In particular, the ${Diff}(M,{\cal{F}})$ does not allow Penrose diagrams. Among the vacuum solutions, some give rise to the structure of the Einstein-Rosen bridge, in which two asymptotically flat regions are connected by a throat with a finite non-zero radius. We also study slowly rotating solutions in such a setup, and obtain all the solutions characterized by an arbitrary function $A_{0}(r)$. The case $A_{0} = 0$ reduces to the slowly rotating Kerr solution obtained in GR.

The evolution of universe in Brans-Dicke (BD) theory is discussed in this paper. Considering a parameterized scenario for BD scalar field $\phi=\phi_{0}a^{\alpha}$ which plays the role of gravitational "constant" $G$, we apply the Markov Chain Monte Carlo method to investigate a global constraints on BD theory with a self-interacting potential according to the current observational data: Union2 dataset of type supernovae Ia (SNIa), high-redshift Gamma-Ray Bursts (GRBs) data, observational Hubble data (OHD), the cluster X-ray gas mass fraction, the baryon acoustic oscillation (BAO), and the cosmic microwave background (CMB) data. It is shown that an expanded universe from deceleration to acceleration is given in this theory, and the constraint results of dimensionless matter density $\Omega_{0m}$ and parameter $\alpha$ are, $\Omega_{0m}=0.286^{+0.037+0.050}_{-0.039-0.047}$ and $\alpha=0.0046^{+0.0149+0.0171}_{-0.0171-0.0206}$ which is consistent with the result of current experiment exploration, $\mid\alpha\mid \leq 0.132124$. In addition, we use the geometrical diagnostic method, jerk parameter $j$, to distinguish the BD theory and cosmological constant model in Einstein's theory of general relativity.

We revisit the vertical structure of neutrino dominated accretion flows (NDAFs) in spherical coordinates with a new boundary condition based on the mechanical equilibrium. The solutions show that NDAF is significantly thick. The Bernoulli parameter and neutrino trapping are determined by the mass accretion rate and the viscosity parameter. According to the distribution of the Bernoulli parameter, the possible outflow may appear in the outer region of the disk. The neutrino trapping can essentially affect the neutrino radiation luminosity. The vertical structure of NDAF is like a "sandwich", and the multilayer accretion may account for the flares in gamma-ray bursts.

For the first time ever, we derive equations governing the time-evolution of fully relativistic slim accretion disks in the Kerr metric, and numerically construct their detailed non-stationary models. We discuss applications of these general results to a possible limit-cycle behavior of thermally unstable disks. Our equations and numerical method are applicable in a wide class of possible viscosity prescriptions, but in this paper we use a diffusive form of the "standard alpha prescription" that assumes the viscous torque is proportional to the total pressure. In this particular case, we find that the parameters which dominate the limit-cycle properties are the mass-supply rate and the value of the alpha-viscosity parameter. Although the duration of the cycle (or the outburst) does not exhibit any clear dependence on the black hole spin, the maximal outburst luminosity (in the Eddington units) is positively correlated with the spin value. We suggest a simple method for a rough estimate of the black hole spin based on the maximal luminosity and the ratio of outburst to cycle durations. We also discuss a temperature-luminosity relation for the Kerr black hole accretion discs limit-cycle. Based on these results we discuss the limit-cycle behavior observed in microquasar GRS 1915+105. We also extend this study to several non-standard viscosity prescriptions, including a "delayed heating" prescription recently stimulated by the recent MHD simulations of accretion disks.

We use the Markov Chain Monte Carlo method to investigate a global constraints on the modified Chaplygin gas (MCG) model as the unification of dark matter and dark energy from the latest observational data: the Union2 dataset of type supernovae Ia (SNIa), the observational Hubble data (OHD), the cluster X-ray gas mass fraction, the baryon acoustic oscillation (BAO), and the cosmic microwave background (CMB) data. In a flat universe, the constraint results for MCG model are, $\Omega_{b}h^{2}=0.02263^{+0.00184}_{-0.00162}$ ($1\sigma$) $^{+0.00213}_{-0.00195}$ $(2\sigma)$, $B_{s}=0.7788^{+0.0736}_{-0.0723}$ ($1\sigma$) $^{+0.0918}_{-0.0904}$ $(2\sigma)$, $\alpha=0.1079^{+0.3397}_{-0.2539}$ ($1\sigma$) $^{+0.4678}_{-0.2911}$ $(2\sigma)$, $B=0.00189^{+0.00583}_{-0.00756}$ ($1\sigma$) $^{+0.00660}_{-0.00915}$ $(2\sigma)$, and $H_{0}=70.711^{+4.188}_{-3.142}$ ($1\sigma$) $^{+5.281}_{-4.149}$ $(2\sigma)$.

Kinematical models are constrained by the latest observational data from geometry-distance measurements, which include 557 type Ia supernovae (SNIa) Union2 data and 15 observational Hubble data. Considering two parameterized deceleration parameter, the values of current deceleration parameter $q_{0}$, jerk parameter $j_{0}$ and transition redshift $z_{T}$, are obtained. Furthermore, we show the departures for two parameterized kinematical models from $\Lambda$CDM model according to the evolutions of jerk parameter $j(z)$. Also, it is shown that the constraint on jerk parameter $j(z)$ is weak by the current geometrical observed data.

The stress evolution process is taken into account in the linear stability analysis of standard thin accretion discs. We find that the growth rate of thermally unstable modes can decrease significantly owing to the stress delay, which may help to understand the quasi-periodic variability of GRS 1915+105. We also discuss possible application of stress evolution to the stability of Shapiro-Lightman-Eardley disc.

The standard thin disk model predicts that when the accretion rate is over a small fraction of the Eddington rate, which corresponds to $L \ga 0.06 L_{Edd}$, the inner region of the disk is radiation-pressure-dominated and thermally unstable. However, observations of the high/soft state of black hole X-ray binaries with luminosity well within this regime ($0.01L_{Edd} \la L \la 0.5L_{Edd}$) indicate that the disk has very little variability, i.e., quite stable. Recent radiation magnetohydrodynamic simulations of a vertically stratified shearing box have confirmed the absence of the thermal instability. In this paper, we revisit the thermal stability by linear analysis, taking into account the role of magnetic field in the accretion flow. By assuming that the field responses negatively to a positive temperature perturbation, we find that the threshold of accretion rate above which the disk becomes thermally unstable increases significantly compared with the case of not considering the role of magnetic field. This accounts for the stability of the observed sources with high luminosities. Our model also presents a possible explanation as to why only GRS 1915+105 seems to show thermally unstable behavior. This peculiar source holds the highest accretion rate (or luminosity) among the known high state sources, which is well above the accretion rate threshold of the instability.

The luminescent properties of CsI(Na) crystals are studied in this report. By using a TDS3054C oscilloscope with a sampling frequency of 5 GS/s, we find out that nuclear recoil signals are dominated by very fast light pulse with a decay time of ~20 ns, while {\gamma}-ray signals have a decay time of ~600 ns. The wavelength of nuclear recoil and {\gamma}-ray signals are also different. The study of n/{\gamma} separation shows that the rejection factor can reach an order of 10-7 with signal efficiency more than 80% at an equivalent electron recoil energy of 20 keV or more. Such a property makes CsI(Na) an ideal candidate for dark matter searches.

In this paper, we study the motion of photons around a Kehagias-Sfetsos (KS) black hole and obtain constraints on IR modified Ho$\check{r}$ava gravity without cosmological constant ($\sim \Lambda_{W}$). An analytic formula for the light deflection angle is obtained. For a propagating photon, the deflection angle $\delta \phi$ increases with large values of the Ho$\check{r}$ava gravity parameter $\omega$. Under the UV limit $\omega \longrightarrow \infty$, deflection angle reduces to the result of usual Schwarzschild case, $4GM/R$. It is also found that with increasing scale of astronomical observation system the Ho$\check{r}$ava-Lifshitz gravity should satisfy $|\omega M^2|>1.1725 \times10^{-16}$ with 12% precision for Earth system, $|\omega M^2| > 8.27649 \times 10^{-17}$ with 17% precision for Jupiter system and $|\omega M^2| > 8.27650\times 10^{-15}$ with 0.17% precision for solar system.

In this paper, we study the motion of photons around a Kehagias-Sfetsos (KS) black hole and obtain constraints on IR modified Ho$\check{r}$ava gravity without cosmological constant ($\sim \Lambda_{W}$). An analytic formula for the light deflection angle is obtained. For a propagating photon, the deflection angle $\delta \phi$ increases with large values of the Ho$\check{r}$ava gravity parameter $\omega$. Under the UV limit $\omega \longrightarrow \infty$, deflection angle reduces to the result of usual Schwarzschild case, $4GM/R$. It is also found that with increasing scale of astronomical observation system the Ho$\check{r}$ava-Lifshitz gravity should satisfy $|\omega M^2|>1.1725 \times10^{-16}$ with 12% precision for Earth system, $|\omega M^2| > 8.27649 \times 10^{-17}$ with 17% precision for Jupiter system and $|\omega M^2| > 8.27650\times 10^{-15}$ with 0.17% precision for solar system.

We present significantly improved proper motion measurements of the Milky Way's central stellar cluster. These improvements are made possible by refining our astrometric reference frame with a new geometric optical distortion model for the W. M. Keck II 10 m telescope's Adaptive Optics camera (NIRC2) in its narrow field mode. For the first time, this distortion model is constructed from on-sky measurements, and is made available to the public. When applied to widely dithered images, it produces residuals in the separations of stars that are a factor of ~3 smaller compared to the outcome using previous models. By applying this new model, along with corrections for differential atmospheric refraction, to widely dithered images of SiO masers at the Galactic center, we improve our ability to tie into the precisely measured radio Sgr A*-rest frame. The resulting infrared reference frame is ~2-3 times more accurate and stable than earlier published efforts. In this reference frame, Sgr A* is localized to within a position of 0.6 mas and a velocity of 0.09 mas/yr, or ~3.4 km/s at 8 kpc (1 sigma). While earlier proper motion studies defined a reference frame by assuming no net motion of the stellar cluster, this approach is fundamentally limited by the cluster's intrinsic dispersion and therefore will not improve with time. We define a reference frame with SiO masers and this reference frame's stability should improve steadily with future measurements of the SiO masers in this region. This is essential for achieving the necessary reference frame stability required to detect the effects of general relativity and extended mass on short-period stars at the Galactic center.

In this paper, motivated from the fact that a de Sitter cosmological boundary corresponds to a positive cosmological constant, we consider time variable cosmological constants, dubbed {\it horizon cosmological constants}. The horizon cosmological constants correspond to Hubble horizon, future event horizon and particle horizon are discussed respectively. When the Hubble horizon is taken as a cosmological length scale, the effective equation of state of horizon cosmological constant is quintessence-like. The values of model parameter $c$ will determine the current status of our universe. When particle horizon is taken as the cosmological length scale, non viable cosmological model can be obtained for the requirement of $\Omega_{\Lambda}<1/3$ which conflicts with current comic observations. When the future event horizon is taken as the role of cosmological length scale, the forms of effective equation of state of horizon cosmological constants are the the same as the holographic ones. But, their evolutions are different because of the effective interaction with cold dark matter.

Large ground-based telescopes equipped with adaptive optics (AO) systems have ushered in a new era of high-resolution infrared photometry and astrometry. Relative astrometric accuracies of <0.2 mas have already been demonstrated from infrared images with spatial resolutions of 55-95 mas resolution over 10-20'' fields of view. Relative photometric accuracies of 3% and absolute photometric accuracies of 5%-20% are also possible. I will review improvements and current limitations in astrometry and photometry with adaptive optics of crowded stellar fields. These capabilities enable experiments such as measuring orbits for brown dwarfs and exoplanets, studying our Galaxy's supermassive black hole and its environment, and identifying individual stars in young star clusters, which can be used test the universality of the initial mass function.

The InfraRed Imaging Spectrograph (IRIS) is a first-light instrument being designed for the Thirty Meter Telescope (TMT). IRIS is a combination of an imager that will cover a 16.4" field of view at the diffraction limit of TMT (4 mas sampling), and an integral field unit spectrograph that will sample objects at 4-50 mas scales. IRIS will open up new areas of observational parameter space, allowing major progress in diverse fields of astronomy. We present the science case and resulting requirements for the performance of IRIS. Ultimately, the spectrograph will enable very well-resolved and sensitive studies of the kinematics and internal chemical abundances of high-redshift galaxies, shedding light on many scenarios for the evolution of galaxies at early times. With unprecedented imaging and spectroscopy of exoplanets, IRIS will allow detailed exploration of a range of planetary systems that are inaccessible with current technology. By revealing details about resolved stellar populations in nearby galaxies, it will directly probe the formation of systems like our own Milky Way. Because it will be possible to directly characterize the stellar initial mass function in many environments and in galaxies outside of the the Milky Way, IRIS will enable a greater understanding of whether stars form differently in diverse conditions. IRIS will reveal detailed kinematics in the centers of low-mass galaxies, allowing a test of black hole formation scenarios. Finally, it will revolutionize the characterization of reionization and the first galaxies to form in the universe.

Deep Keck/NIRC2 HK'L' observations of the Arches cluster near the Galactic center reveal a significant population of near-infrared excess sources. We combine the L'-band excess observations with K'-band proper motions, to confirm cluster membership of excess sources in a starburst cluster for the first time. The robust removal of field contamination provides a reliable disk fraction down to our completeness limit of H=19 mag, or about 5 Msun at the distance of the Arches. Of the 24 identified sources with K'-L' > 2.0 mag, 21 have reliable proper motion measurements, all of which are proper motion members of the Arches cluster. VLT/SINFONI K'-band spectroscopy of three excess sources reveals strong CO bandhead emission, which we interpret as the signature of dense circumstellar disks. The detection of strong disk emission from the Arches stars is surprising in view of the high mass of the B-type main sequence host stars of the disks and the intense starburst environment. We find a disk fraction of 6 +/- 2% among B-type stars in the Arches cluster. A radial increase in the disk fraction from 3 to 10% suggests rapid disk destruction in the immediate vicinity of numerous O-type stars in the cluster core. A comparison between the Arches and other high- and low-mass star-forming regions provides strong indication that disk depletion is significantly more rapid in compact starburst clusters than in moderate star-forming environments.

So far, there have been no theories or observational data that deny the presence of interaction between dark energy and dark matter. We extend naturally the holographic dark energy (HDE) model, proposed by Granda and Oliveros, in which the dark energy density includes not only the square of the Hubble scale, but also the time derivative of the Hubble scale to the case with interaction and the analytic forms for the cosmic parameters are obtained under the specific boundary conditions. The various behaviors concerning the cosmic expansion depend on the introduced numerical parameters which are also constrained. The more general interacting model inherits the features of the previous ones of HDE, keeping the consistency of the theory.

We investigate observational constraints on the generalized Chaplygin gas (GCG) model as the unification of dark matter and dark energy from the latest observational data: the Union SNe Ia data, the observational Hubble data, the SDSS baryon acoustic peak and the five-year WMAP shift parameter. It is obtained that the best fit values of the GCG model parameters with their confidence level are $A_{s}=0.73^{+0.06}_{-0.06}$ ($1\sigma$) $^{+0.09}_{-0.09}$ $(2\sigma)$, $\alpha=-0.09^{+0.15}_{-0.12}$ ($1\sigma$) $^{+0.26}_{-0.19}$ $(2\sigma)$. Furthermore in this model, we can see that the evolution of equation of state (EOS) for dark energy is similar to quiessence, and its current best-fit value is $w_{0de}=-0.96$ with the $1\sigma$ confidence level $-0.91\geq w_{0de}\geq-1.00$.

In this paper, the properties of dark energy are investigated according to the parameterized deceleration parameter $q(z)$, which is used to describe the extent of the accelerating expansion of the universe. The potential of dark energy $V(\phi)$ and the cosmological parameters, such as the dimensionless energy density $\Omega_{\phi}$, $\Omega_{m}$, and the state parameter $w_\phi$, are connected to it. Concretely, by giving two kinds of parameterized deceleration parameters $q(z)=a+\frac{bz}{1+z}$ and $q(z)=1/2+\frac{az+b}{(1+z)^2}$, the evolution of these parameters and the reconstructed potentials $V(\phi)$ are plotted and analyzed. It's found that the potentials run away with the evolution of universe.

The integrated Sachs-Wolfe (ISW) effect is an important implication for dark energy. In this paper, we have calculated the power spectrum of the ISW effect in the time varying vacuum cosmological model, where the model parameter $\beta=4.407$ is obtained by the observational constraint of the growth rate. It's found that the source of the ISW effect is not only affected by the different evolutions of the Hubble function $H(a)$ and the dimensionless matter density $\Omega_m(a)$, but also by the different growth function $D_+(a)$, all of which are changed due to the presence of matter production term in the time varying vacuum model. However, the difference of the ISW effect in $\Lambda(t)\textmd{CDM}$ model and $\Lambda \textmd{CDM}$ model is lessened to a certain extent due to the integration from the time of last scattering to the present. It's implied that the observations of the galaxies with high redshift are required to distinguish the two models.

A new diagnostic method, $Om$ is applied to generalized Chaplygin gas (GCG) model as the unification of dark matter and dark energy. On the basis of the recently observed data: the Union supernovae, the observational Hubble data, the SDSS baryon acoustic peak and the five-year WMAP shift parameter, we show the discriminations between GCG and $\Lambda$CDM model. Furthermore, it is calculated that the current equation of state of dark energy $w_{0de}=-0.964$ according to GCG model.

In this Letter, a modified Chaplygin gas (MCG) model of unifying dark energy and dark matter with the exotic equation of state $p_{MCG}=B\rho_{MCG} -\frac A{\rho_{MCG}^\alpha}$ is constrained from recently observed data: the 182 Gold SNe Ia, the 3-year WMAP and the SDSS baryon acoustic peak. It is shown that the best fit value of the three parameters ($B$,$B_{s}$,$\alpha$) in MCG model are (-0.085,0.822,1.724). Furthermore, we find the best fit $w(z)$ crosses -1 in the past and the present best fit value $w(0)=-1.114<-1$, and the $1\sigma$ confidence level of $w(0)$ is $-0.946\leq w(0)\leq-1.282$. Finally, we find that the MCG model has the smallest $\chi^{2}_{min}$ value in all eight given models. According to the Alaike Information Criterion (AIC) of model selection, we conclude that recent observational data support the MCG model as well as other popular models.

We use the Markov Chain Monte Carlo method to investigate a global constraints on the generalized Chaplygin gas (GCG) model as the unification of dark matter and dark energy from the latest observational data: the Constitution dataset of type supernovae Ia (SNIa), the observational Hubble data (OHD), the cluster X-ray gas mass fraction, the baryon acoustic oscillation (BAO), and the cosmic microwave background (CMB) data. In a non-flat universe, the constraint results for GCG model are, $\Omega_{b}h^{2}=0.0235^{+0.0021}_{-0.0018}$ ($1\sigma$) $^{+0.0028}_{-0.0022}$ $(2\sigma)$, $\Omega_{k}=0.0035^{+0.0172}_{-0.0182}$ ($1\sigma$) $^{+0.0226}_{-0.0204}$ $(2\sigma)$, $A_{s}=0.753^{+0.037}_{-0.035}$ ($1\sigma$) $^{+0.045}_{-0.044}$ $(2\sigma)$, $\alpha=0.043^{+0.102}_{-0.106}$ ($1\sigma$) $^{+0.134}_{-0.117}$ $(2\sigma)$, and $H_{0}=70.00^{+3.25}_{-2.92}$ ($1\sigma$) $^{+3.77}_{-3.67}$ $(2\sigma)$, which is more stringent than the previous results for constraint on GCG model parameters. Furthermore, according to the information criterion, it seems that the current observations much support $\Lambda$CDM model relative to the GCG model.

Motivated by the recent work about a new physical interpretation of quasinormal modes by Maggiore, we investigate the quantization of near-extremal Schwarzschild-de Sitter black holes in the four dimensional spacetime. Following Kunstatter's method, we derive the area and entropy spectrum of near-extremal Schwarzschild-de Sitter black holes which differs from Setare's result. Furthermore, we find that the derived a universal area spectrum is $2\pi n$ which is equally spaced.

We use observational data from Type Ia Supernovae (SN), Baryon Acoustic Oscillations (BAO), Cosmic Microwave Background (CMB) and observational Hubble data (OHD), and the Markov Chain Monte Carlo (MCMC) method, to constrain the cosmological scenario of holographic dark energy with varying gravitational constant. We consider both flat and non-flat background geometry, and we present the corresponding constraints and contour-plots of the model parameters. We conclude that the scenario is compatible with observations. In 1$\sigma$ we find $\Omega_{\Lambda0}=0.72^{+0.03}_{-0.03}$, $\Omega_{k0}=-0.0013^{+0.0130}_{-0.0040}$, $c=0.80^{+0.19}_{-0.14}$ and $\Delta_G\equiv G'/G=-0.0025^{+0.0080}_{-0.0050}$, while for the present value of the dark energy equation-of-state parameter we obtain $w_0=-1.04^{+0.15}_{-0.20}$.

A model of jet precession driven by a neutrino-cooled disc around a spinning black hole is present in order to explain the temporal structure and spectral evolution of gamma-ray bursts (GRBs). The differential rotation of the outer part of a neutrino dominated accretion disc may result in precession of the inner part of the disc and the central black hole, hence drives a precessed jet via neutrino annihilation around the inner part of the disc. Both analytic and numeric results for our model are present. Our calculations show that a black hole-accretion disk system with black hole mass $M \simeq 3.66 M_\odot$, accretion rate $\dot{M} \simeq 0.54 M_\odot \rm s^{-1}$, spin parameter $a=0.9$ and viscosity parameter $\alpha=0.01$ may drive a precessed jet with period P=1 s and luminosity $L=10^{51}$ erg s$^{-1}$, corresponding to the scenario for long GRBs. A precessed jet with $P=0.1$s and $L=10^{50}$ erg s$^{-1}$ may be powered by a system with $M \simeq 5.59 M_\odot$, $\dot{M} \simeq 0.74 M_\odot \rm s^{-1}$, $a=0.1$, and $\alpha=0.01$, possibly being responsible for the short GRBs. Both the temporal and spectral evolution in GRB pulse may explained with our model. GRB central engines likely power a precessed jet driven by a neutrino-cooled disc. The global GRB lightcurves thus could be modulated by the jet precession during the accretion timescale of the GRB central engine. Both the temporal and spectral evolution in GRB pulse may be due to an viewing effect due to the jet precession.

(Abridged) We report on the structure of the nuclear star cluster in the innermost 0.16 pc of the Galaxy as measured by the number density profile of late-type giants. Using laser guide star adaptive optics in conjunction with the integral field spectrograph, OSIRIS, at the Keck II telescope, we are able to differentiate between the older, late-type (~ 1 Gyr) stars, which are presumed to be dynamically relaxed, and the unrelaxed young (~ 6 Myr) population. This distinction is crucial for testing models of stellar cusp formation in the vicinity of a black hole, as the models assume that the cusp stars are in dynamical equilibrium in the black hole potential. We find that contamination from young stars is significant, with more than twice as many young stars as old stars in our sensitivity range (K < 15.5) within the central arcsecond. Based on the late-type stars alone, the surface stellar number density profile, is flat, with a projected power law slope of -0.26+-0.24. These results are consistent with the nuclear star cluster having no cusp, with a core profile that is significantly flatter than predicted by most cusp formation theories. Here, we also review the methods for further constraining the true three-dimensional radial profile using kinematic measurements. Precise acceleration measurements in the plane of the sky as well as along the line of sight has the potential to directly measure the density profile to establish whether there is a "hole" in the distribution of late-type stars in the inner 0.1 pc.

We report a factor of $\sim$3 improvement in Keck Laser Guide Star Adaptive Optics (LGSAO) astrometric measurements of stars near the Galaxy's supermassive black hole (SMBH). By carrying out a systematic study of M92, we have produced a more accurate model for the camera's optical distortion. Updating our measurements with this model, and accounting for differential atmospheric refraction, we obtain estimates of the SMBH properties that are a factor of $\sim$2 more precise, and most notably, increase the likelihood that the black hole is at rest with respect to the nuclear star cluster. These improvements have also allowed us to extend the radius to which stellar orbital parameter estimates are possible by a factor of 2.

We present progenitor-star detections, light curves, and optical spectra of SN2009ip and the 2009 optical transient in UGC2773 (U2773-OT), which were not genuine SNe. Precursor variability in the decade before outburst indicates that both of the progenitor stars were LBVs. Their pre-outburst light curves resemble the S Doradus phases that preceded giant eruptions of eta Carinae and SN1954J (V12 in NGC2403), with intermediate progenitor luminosities. HST detections a decade before discovery indicate that the SN2009ip and U2773-OT progenitors were supergiants with likely initial masses of 50-80 Msun and $\ga$20 Msun, respectively. Both outbursts had spectra befitting known LBVs, although in different physical states. SN 2009ip exhibited a hot LBV spectrum with characteristic speeds of 550 km/s, plus faster material up to 5000 km/s, resembling the slow Homunculus and fast blast wave of eta Carinae. U2773-OT shows a forest of narrow absorption and emission lines comparable to that of S Dor in its cool state, plus [CaII] emission and an IR excess indicative of dust, similar to SN2008S and N300-OT. [CaII] emission is probably tied to a dusty pre-outburst environment, and not the outburst mechanism. SN2009ip and U2773-OT may provide a critical link between historical LBV eruptions, while U2773-OT may provide a link between LBVs and SN2008S and N300-OT. Future searches will uncover more examples of precursor LBV variability of this kind, providing key clues that may help unravel the instability driving LBVs.

We revisit the vertical structure of neutrino-dominated accretion flows in spherical coordinates. We stress that the flow should be geometrically thick when advection becomes dominant. In our calculation, the luminosity of neutrino annihilation is enhanced by one or two orders of magnitude. The empty funnel along the rotation axis can naturally explain the neutrino annihilable ejection.

We investigate the stability of strange quark matter and the properties of the corresponding strange stars, within a wide range of quark mass scaling. The calculation shows that the resulting maximum mass always lies between 1.5 solor mass and 1.8 solor mass for all the scalings chosen here. Strange star sequences with a linear scaling would support less gravitational mass, and a change (increase or decrease) of the scaling around the linear scaling would lead to a larger maximum mass. Radii invariably decrease with the mass scaling. Then the larger the scaling, the faster the star might spin. In addition, the variation of the scaling would cause an order of magnitude change of the strong electric field on quark surface, which is essential to support possible crusts of strange stars against gravity and may then have some astrophysical implications.

We revisit the vertical structure of black hole accretion disks in spherical coordinates. By comparing the advective cooling with the viscous heating, we show that advection-dominated disks are geometrically thick, i.e., with the half-opening angle > 2\pi/5, rather than slim as supposed previously in the literature.

We provide a detailed analysis of 12 XMM observations of the nearest persistent extragalactic ultraluminous X-ray source (ULX), M33 X-8. No significant spectral evolution is detected between the observations, therefore we combine the individual observations to increase the signal-to-noise ratio for spectral fitting. The combined spectra are best fitted by the self-consistent p-free disk plus power-law component model with p = 0.571_{-0.030}^{+0.032}, kT_{in} = 1.38_{-0.08}^{+0.09} keV, and the flux ratio of the p-free disk component to the power-law component being 0.63:0.37 in the 0.3 -- 10 keV band. The fitting indicates that the black hole in M33 X-8 is of \sim 10 M_{\odot} and accretes at a super-Eddington rate (\sim 1.5 L_{Edd}), and the phase of the accretion disk is close to a slim disk (p = 0.5). We report, for the first time, that an extra power-law component is required in addition to the p-free disk model for ULXs. In super-Eddington cases, the power-law component may possibly result from the optically thin inner region f the disk or a comptonized corona similar to that of a standard thin disk.

We report on the structure of the nuclear star cluster in the innermost 0.16 pc of the Galaxy as measured by the number density profile of late-type giants. Using laser guide star adaptive optics in conjunction with the integral field spectrograph, OSIRIS, at the Keck II telescope, we are able to differentiate between the older, late-type ($\sim$ 1 Gyr) stars, which are presumed to be dynamically relaxed, and the unrelaxed young ($\sim$ 6 Myr) population. This distinction is crucial for testing models of stellar cusp formation in the vicinity of a black hole, as the models assume that the cusp stars are in dynamical equilibrium in the black hole potential. Based on the late-type stars alone, the surface stellar number density profile, $\Sigma(R) \propto R^{-\Gamma}$, is flat, with $\Gamma = -0.27\pm0.19$. Monte Carlo simulations of the possible de-projected volume density profile, n(r) $\propto r^{-\gamma}$, show that $\gamma$ is less than 1.0 at the 99.73 % confidence level. These results are consistent with the nuclear star cluster having no cusp, with a core profile that is significantly flatter than predicted by most cusp formation theories, and even allows for the presence of a central hole in the stellar distribution. Of the possible dynamical interactions that can lead to the depletion of the red giants observable in this survey -- stellar collisions, mass segregation from stellar remnants, or a recent merger event -- mass segregation is the only one that can be ruled out as the dominant depletion mechanism. The lack of a stellar cusp around a supermassive black hole would have important implications for black hole growth models and inferences on the presence of a black hole based upon stellar distributions.

We take the vertical distribution of the radial and azimuthal velocity into account in spherical coordinates, and find that the analytic relation c_{s0}/(v_K \Theta) = [(\gamma -1)/(2\gamma)]^{1/2} is valid for both geometrically thin and thick accretion flows, where c_{s0} is the sound speed on the equatorial plane, v_K is the Keplerian velocity, \Theta is the half-opening angle of the flow, and \gamma is the adiabatic index.

In this paper, a holographic dark energy model, dubbed Ricci dark energy, is confronted with cosmological observational data from type Ia supernovae (SN Ia), baryon acoustic oscillations (BAO) and cosmic microwave background (CMB). By using maximum likelihood method, it is found out that Ricci dark energy model is a viable candidate of dark energy model with the best fit parameters: $\Omega_{m0}=0.34\pm 0.04$, $\alpha=0.38\pm 0.03$ with $1\sigma$ error. Here, $\alpha$ is a dimensionless parameter related with Ricci dark energy $\rho_{R}$ and Ricci scalar $R$, i.e., $\rho_{R}\propto \alpha R$.

In this paper, we consider generalized holographic and Ricci dark energy models where the energy densities are given as $\rho_{R}=3c^2M^{2}_{pl}Rf(H^2/R)$ and $\rho_{h}=3c^2M^{2}_{pl}H^2g(R/H^2)$ respectively, here $f(x),g(y)$ are positive defined functions of dimensionless variables $H^2/R$ or $R/H^2$. It is interesting that holographic and Ricci dark energy densities are recovered or recovered interchangeably when the function $f(x)=g(y)\equiv 1$ or$f=g\equiv Id$ is taken respectively (for example $f(x),g(x)=1-\epsilon(1-x)$, $\epsilon=0 \text{or} 1$ respectively). Also, when $f(x)\equiv xg(1/x)$ is taken, the Ricci and holographic dark energy models are equivalents to a generalized one. When the simple forms $f(x)=1-\epsilon(1-x)$ and $g(y)=1-\eta(1-y)$ are taken as examples, by using current cosmic observational data, generalized dark energy models are researched. As expected, in these cases, the results show that they are equivalent ($\epsilon=1-\eta=1.312$) and Ricci-like dark energy is more favored relative to the holographic one where the Hubble horizon was taken as an IR cut-off. And, the suggestive combination of holographic and Ricci dark energy components would be $1.312 R-0.312H^2$ which is $2.312H^2+1.312\dot{H}$ in terms of $H^2$ and $\dot{H}$.

In this paper, time variable cosmological constant, dubbed {\it age cosmological constant}, is investigated motivated by the fact: any cosmological length scale and time scale can introduce a cosmological constant or vacuum energy density into Einstein's theory. The age cosmological constant takes the form $\rho_{\Lambda}=3c^2M^2_P/t_{\Lambda}^2$, where $t_{\Lambda}$ is the age of our universe or conformal time. The effective equation of state of age cosmological constant are $w^{eff}_{\Lambda}=-1+{2/3}\frac{\sqrt{\Omega_{\Lambda}}}{c}$ and $w^{eff}_{\Lambda}=-1+{2/3}\frac{\sqrt{\Omega_{\Lambda}}}{c}(1+z)$ when the age of universe and conformal time are taken as the role of cosmological time scales respectively. They are the same as the so-called agegraphic dark energy models. However, the evolution history are different from the agegraphic ones for their different evolution equations.

In this paper, linear first order expansion of deceleration parameter $q(z)=q_0+q_1(1-a)$ ($M_1$), constant jerk $j=j_0$ ($M_2$) and third order expansion of luminosity distance ($M_3$) are confronted with cosmic observations: SCP 307 SN Ia, BAO and observational Hubble data (OHD). Likelihood is implemented to find the best fit model parameters. All these models give the same prediction of the evolution of the universe which is undergoing accelerated expansion currently and experiences a transition from decelerated expansion to accelerated expansion. But, the transition redshift depends on the concrete parameterized form of the model assumed. $M_1$ and $M_2$ give value of transition redshift about $z_t\sim 0.6$. $M_3$ gives a larger one, say $z_t\sim 1$. The $\chi^2/dof$ implies almost the same goodness of the models. But, for its badness of evolution of deceleration parameter at high redshift $z>1$, $M_3$ can not be reliable. $M_1$ and $M_2$ are compatible with $\Lambda$CDM model at the $2\sigma$ and $1\sigma$ confidence levels respectively. $M_3$ is not compatible with $\Lambda$CDM model at $2\sigma$ confidence level. From $M_1$ and $M_2$ models, one can conclude that the cosmic data favor a cosmological model having $j_0<-1$.

In this paper, the holographic dark energy in Brans-Dicke theory is confronted by cosmic observations from SN Ia, BAO and CMB shift parameter. The best fit parameters are found in $1\sigma$ region: $\Omega_{h0}=0.683^{+0.035}_{-0.038}$, $c=0.605^{+0.138}_{-0.107}$ and $\alpha=0.00662^{+0.00477}_{-0.00467}$ (equivalently $\omega=905.690^{+637.906}_{-651.471}$ which is less the solar system bound and consistent with other constraint results). With these best fit values of the parameters, it is found the universe is undergoing accelerated expansion, and the current value of equation of state of holographic dark energy $w_{h0}=-1.246^{+0.191}_{-0.144}$ which is phantom like in Brans-Dicke theory. The evolution of effective Newton's constant is also explored.

In this paper, the holographic dark energy model is considered in Brans-Dicke theory where the holographic dark energy density $\rho_{\Lambda} =3c^2 M^{2}_{pl} L^{-2}$ is replaced with $\rho_{h}=3c^2 \Phi(t)L^{-2}$. Here $\Phi(t)=\frac{1}{8\pi G}$ is a time variable Newton constant. With this replacement, it is found that no accelerated expansion universe will be achieved when the Hubble horizon is taken as the role of IR cut-off. When the event horizon is adopted as the IR cut-off, an accelerated expansion universe is obtained. In this case, the equation of state of holographic dark energy $w_h$ takes a modified form $w_h=-{1/3}(1+\alpha+\frac{2}{c}\sqrt{\Omega_{h}})$. In the limit $\alpha\to 0$, the 'standard' holographic dark energy is recovered. In the holographic dark energy dominated epoch, power-law and de Sitter time-space solutions are obtained.

As the closest example of a galactic nucleus, the Galactic center presents an exquisite laboratory for learning about supermassive black holes (SMBH) and their environs. Detailed studies of stellar dynamics deep in the potential well of a galaxy, with exisiting and future large ground-based telescopes, offer several exciting directions in the coming decade. First, it will be possible to obtain precision measurements of the Galaxy's central potential, providing both a unique test of General Relativity (GR) and a detection of the extended dark matter distribution that is predicted to exist around the SMBH. Tests of gravity have not previously been possible on the mass scale of a SMBH. Similarly, only upper limits on the extended matter distribution on small scales currently exist; detection of dark matter on these scales is an important test of Lambda-CDM and the detection of stellar remnants would reveal a population that may dominate the stellar dynamics on the smallest scales. Second, our detailed view of the SMBH and its local gas and stellar environment provides insight into how SMBHs at the centers of galaxies form, grow and interact with their environs as well as on the exotic processes at work in the densest stellar clusters in the Universe. The key questions, still unanswered, of when and how SMBHs formed in the early universe, and the myriad ways in which feedback from SMBHs can affect structure formation, can be informed by directly observing the physical processes operating at the SMBH.

We use recently observed data: the 192 ESSENCE type Ia supernovae (SNe Ia), the 182 Gold SNe Ia, the 3-year WMAP, the SDSS baryon acoustic peak, the X-ray gas mass fraction in clusters and the observational $H(z)$ data to constrain models of the accelerating universe. Combining the 192 ESSENCE data with the observational $H(z)$ data to constrain a parameterized deceleration parameter, we obtain the best fit values of transition redshift and current deceleration parameter $z_{T}=0.632^{+0.256}_{-0.127}$, $q_{0}=-0.788^{+0.182}_{-0.182}$. Furthermore, using $\Lambda$CDM model and two model-independent equation of state of dark energy, we find that the combined constraint from the 192 ESSENCE data and other four cosmological observations gives smaller values of $\Omega_{0m}$ and $q_{0}$, but a larger value of $z_{T}$ than the combined constraint from the 182 Gold data with other four observations. Finally, according to the Akaike information criterion it is shown that the recently observed data equally supports three dark energy models: $\Lambda$CDM, $w_{de}(z)=w_{0}$ and $w_{de}(z)=w_{0}+w_{1}\ln(1+z)$.

We show that there exists a maximal possible accretion rate, beyond which global slim disk solutions cannot be constructed because in the vertical direction the gravitational force would be unable to balance the pressure force to gather the accreted matter. The principle for this restriction is the same as that for the Eddington luminosity and the corresponding critical accretion rate, which were derived for spherical accretion by considering the same force balance in the radial direction. If the assumption of hydrostatic equilibrium is waived and vertical motion is included, this restriction may become even more serious as the value of the maximal possible accretion rate becomes smaller. Previous understanding in the literature that global slim disk solutions could stand for any large accretion rates is due to the overestimation of the vertical gravitational force by using an approximate potential. For accretion flows with large accretion rates at large radii, outflows seem unavoidable in order for the accretion flow to reduce the accretion rate and follow a global solution till the central black hole.

We present the results of near-infrared (2 and 3 microns) monitoring of Sgr A*-IR with 1 min time sampling using the natural and laser guide star adaptive optics (LGS AO) system at the Keck II telescope. Sgr A*-IR was observed continuously for up to three hours on each of seven nights, between 2005 July and 2007 August. Sgr A*-IR is detected at all times and is continuously variable, with a median observed 2 micron flux density of 0.192 mJy, corresponding to 16.3 magnitude at K'. These observations allow us to investigate Nyquist sampled periods ranging from about 2 minutes to an hour. Using Monte Carlo simulations, we find that the variability of Sgr A* in this data set is consistent with models based on correlated noise with power spectra having frequency dependent power law slopes between 2.0 to 3.0, consistent with those reported for AGN light curves. Of particular interest are periods of ~20 min, corresponding to a quasi-periodic signal claimed based upon previous near-infrared observations and interpreted as the orbit of a 'hot spot' at or near the last stable orbit of a spinning black hole. We find no significant periodicity at any time scale probed in these new observations for periodic signals. This study is sensitive to periodic signals with amplitudes greater than 20% of the maximum amplitude of the underlying red noise component for light curves with duration greater than ~2 hours at a 98% confidence limit.

We present new proper motions from the 10 m Keck telescopes for a puzzling population of massive, young stars located within 3.5" (0.14 pc) of the supermassive black hole at the Galactic Center. Our proper motion measurements have uncertainties of only 0.07 mas/yr (3 km/s), which is ~7 times better than previous proper motion measurements for these stars, and enables us to measure accelerations as low as 0.2 mas/yr^2 (7 km/s/yr). Using these measurements, line-of-sight velocities from the literature, and 3D velocities for additional young stars in the central parsec, we constrain the true orbit of each individual star and directly test the hypothesis that the massive stars reside in two stellar disks as has been previously proposed. Analysis of the stellar orbits reveals only one of the previously proposed disks of young stars using a method that is capable of detecting disks containing at least 7 stars. The detected disk contains 50% of the young stars, is inclined by ~115 deg from the plane of the sky, and is oriented at a position angle of ~100 deg East of North. Additionally, the on-disk and off-disk populations have similar K-band luminosity functions and radial distributions that decrease at larger projected radii as \propto r^-2. The disk has an out-of-the-disk velocity dispersion of 28 +/- 6 km/s, which corresponds to a half-opening angle of 7 +/- 2 deg, and several candidate disk members have eccentricities greater than 0.2. Our findings suggest that the young stars may have formed in situ but in a more complex geometry than a simple, thin circular disk.

We report new precision measurements of the properties of our Galaxy's supermassive black hole. Based on astrometric (1995-2007) and radial velocity (2000-2007) measurements from the W. M. Keck 10-meter telescopes, a fully unconstrained Keplerian orbit for the short period star S0-2 provides values for Ro of 8.0+-0.6 kpc, M_bh of 4.1+-0.6x10^6 Mo, and the black hole's radial velocity, which is consistent with zero with 30 km/s uncertainty. If the black hole is assumed to be at rest with respect to the Galaxy, we can further constrain the fit and obtain Ro = 8.4+-0.4 kpc and M_bh = 4.5+-0.4x10^6 Mo. More complex models constrain the extended dark mass distribution to be less than 3-4x10^5 Mo within 0.01 pc, ~100x higher than predictions from stellar and stellar remnant models. For all models, we identify transient astrometric shifts from source confusion and the assumptions regarding the black hole's radial motion as previously unrecognized limitations on orbital accuracy and the usefulness of fainter stars. Future astrometric and RV observations will remedy these effects. Our estimates of Ro and the Galaxy's local rotation speed, which it is derived from combining Ro with the apparent proper motion of Sgr A*, (theta0 = 229+-18 km/s), are compatible with measurements made using other methods. The increased black hole mass found in this study, compared to that determined using projected mass estimators, implies a longer period for the innermost stable orbit, longer resonant relaxation timescales for stars in the vicinity of the black hole and a better agreement with the M_bh-sigma relation.

Energetic flares are observed in the Galactic supermassive black hole Sagittarius A* from radio to X-ray wavelengths. On a few occasions, simultaneous flares have been detected in IR and X-ray observations, but clear counterparts at longer wavelengths have not been seen. We present a flare observed over several hours on 2006 July 17 with the Chandra X-Ray Observatory, the Keck II telescope, the Caltech Submillimeter Observatory, and the Submillimeter Array. All telescopes observed strong flare events, but the submillimeter peak is found to occur nearly 100 minutes after the X-ray peak. Submillimeter polarization data show linear polarization in the excess flare emission, increasing from 9% to 17% as the flare passes through its peak, consistent with a transition from optically thick to thin synchrotron emission. The temporal and spectral behavior of the flare require that the energetic electrons responsible for the emission cool faster than expected from their radiative output. This is consistent with adiabatic cooling in an expanding emission region, with X-rays produced through self-Compton scattering, although not consistent with the simplest model of such expansion. We also present a submillimeter flare that followed a bright IR flare on 2005 July 31. Compared to 2006, this event had a larger peak IR flux and similar submillimeter flux, but it lacked measurable X-ray emission. It also showed a shorter delay between the IR and submillimeter peaks. Based on these events we propose a synchrotron and self-Compton model to relate the submillimeter lag and the variable IR/X-ray luminosity ratio.

We present a unified description of all the three known classes of optically thick accretion disks around black holes, namely Shakura-Sunyaev disks, slim disks, and neutrino-dominated accretion flows (NDAFs). It is found that NDAFs have both a maximal and a minimal possible mass accretion rate at their each radius. This may be suggestive of an interpretation for the origin of X-ray flares observed in gamma-ray bursts.

We present the first measurement of the proper motion of the young, compact Arches cluster near the Galactic center from near-infrared adaptive optics (AO) data taken with the recently commissioned laser-guide star (LGS) at the Keck 10-m telescope. The excellent astrometric accuracy achieved with LGS-AO provides the basis for a detailed comparison with VLT/NAOS-CONICA data taken 4.3 years earlier. Over the 4.3 year baseline, a spatial displacement of the Arches cluster with respect to the field population is measured to be 24.0 +/- 2.2 mas, corresponding to a proper motion of 5.6 +/- 0.5 mas/yr or 212 +/- 29 km/s at a distance of 8 kpc. In combination with the known line-of-sight velocity of the cluster, we derive a 3D space motion of 232 +/- 30 km/s of the Arches relative to the field. The large proper motion of the Arches cannot be explained with any of the closed orbital families observed in gas clouds in the bar potential of the inner Galaxy, but would be consistent with the Arches being on a transitional trajectory from x1 to x2 orbits. We investigate a cloud-cloud collision as the possible origin for the Arches cluster. The integration of the cluster orbit in the potential of the inner Galaxy suggests that the cluster passes within 10 pc of the supermassive black hole only if its true GC distance is very close to its projected distance. A contribution of young stars from the Arches cluster to the young stellar population in the inner few parsecs of the GC thus appears increasingly unlikely. The measurement of the 3D velocity and orbital analysis provides the first observational evidence that Arches-like clusters do not spiral into the GC. This confirms that no progenitor clusters to the nuclear cluster are observed at the present epoch.

We report the first time-series of broadband infrared (IR) color measurements of Sgr A*, the variable emission source associated with the supermassive black hole at the Galactic Center. Using the laser and natural guide star AO systems on the Keck II telescope, we imaged Sgr A* in multiple near-infrared broadband filters with a typical cycle time of ~3 min during 4 observing runs (2005-2006), two of which were simultaneous with Chandra X-ray measurements. In spite of the large range of dereddened flux densities for Sgr A* (2-30 mJy), all of our near-IR measurements are consistent with a constant spectral index of alpha = -0.6+-0.2. Furthermore, this value is consistent with the spectral indices observed at X-ray wavelengths during nearly all outbursts; which is consistent with the synchrotron self-Compton model for the production of the X-ray emission. During the coordinated observations, one IR outburst occurs <36 min after a possibly associated X-ray outburst, while several similar IR outbursts show no elevated X-ray emission. A variable X-ray to IR ratio and constant infrared spectral index challenge the notion that the IR and X-ray emission are connected to the same electrons. We, therefore, posit that the population of electrons responsible for both the IR and X-ray emission are generated by an acceleration mechanism that leaves the bulk of the electron energy distribution responsible for the IR emission unchanged, but has a variable high-energy cutoff. Occasionally a tail of electrons >1 GeV is generated, and it is this high-energy tail that gives rise to the X-ray outbursts. One possible explanation for this type of variation is from the turbulence induced by a magnetorotational instability, in which the outer scale length of the turbulence varies and changes the high-energy cutoff.

We propose in this paper a quintom model of dark energy with a single scalar field $\phi$ given by the lagrangian ${\cal L}=-V(\phi)\sqrt{1-\alpha^\prime\nabla_{\mu}\phi\nabla^{\mu}\phi +\beta^\prime \phi\Box\phi}$. In the limit of $\beta^\prime\to$0 our model reduces to the effective low energy lagrangian of tachyon considered in the literature. We study the cosmological evolution of this model, and show explicitly the behaviors of the equation of state crossing the cosmological constant boundary.

We present a numerical method for spatially 1.5-dimensional and time-dependent studies of accretion disks around black holes, that is originated from a combination of the standard pseudo-spectral method and the adaptive domain decomposition method existing in the literature, but with a number of improvements in both the numerical and physical senses. In particular, we introduce a new treatment for the connection at the interfaces of decomposed subdomains, construct an adaptive function for the mapping between the Chebyshev-Gauss-Lobatto collocation points and the physical collocation points in each subdomain, and modify the over-simplified 1-dimensional basic equations of accretion flows to account for the effects of viscous stresses in both the azimuthal and radial directions. Our method is verified by reproducing the best results obtained previously by Szuszkiewicz & Miller on the limit-cycle behavior of thermally unstable accretion disks with moderate viscosity. A new finding is that, according to our computations, the Bernoulli function of the matter in such disks is always and everywhere negative, so that outflows are unlikely to originate from these disks. We are encouraged to study the more difficult case of thermally unstable accretion disks with strong viscosity, and wish to report our results in a subsequent paper.

We show that when the gravitational force is correctly calculated in dealing with the vertical hydrostatic equilibrium of black hole accretion disks, the relationship that is valid for geometrically thin disks, i.e., $c_s/\Omega_K H =$ constant, where $c_s$ is the sound speed, $\Omega_K$ is the Keplerian angular velocity, and $H$ is the half-thickness of the disk, does not hold for slim disks. More importantly, by adopting the correct vertical gravitational force in studies of thermal equilibrium solutions, we find that there exists a maximally possible accretion rate for each radius in the outer region of optically thick accretion flows, so that only the inner region of these flows can possibly take the form of slim disks, and strong outflows from the outer region are required to reduce the accretion rate in order for slim disks to be realized.

Neutrino-cooled hyperaccretion disks around stellar mass black holes are plausible candidates for the central engine of gamma-ray bursts. We calculate the one-dimensional structure and the annihilation luminosity of such disks. The neutrino optical depth is of crucial importance in determining the neutrino cooling rate and is in turn dependent on the electron fraction, the free nucleon fraction, and the electron degeneracy, with given density and temperature of the disk matter. We construct a bridging formula for the electron fraction that works for various neutrino optical depths, and give exact definitions for the free proton fraction and free neutron fraction. We show that the electron degeneracy has important effects in the sense that it enlarges the absorption optical depth for neutrinos, and it along with the neutronization processes favored by high temperature cause the electron fraction to drop to be below 0.1 in the inner region of the disk. The resulting neutrino annihilation luminosity is considerably reduced comparing with that obtained in previous works where the electron degeneracy was not considered and the electron fraction was simply taken to be 0.5, but it is still likely to be adequate for gamma-ray bursts, and it is ejected mainly from the inner region of the disk and has an anisotropic distribution.

We report the results of a diffraction-limited, photometric variability study of the central 5"x5" of the Galaxy conducted over the past 10 years using speckle imaging techniques on the W. M. Keck I 10 m telescope. Within our limiting magnitude of mK < 16 mag for images made from a single night of data, we find a minimum of 15 K[2.2 micron]-band variable stars out of 131 monitored stars. The only periodic source in our sample is the previously identified variable IRS 16SW, for which we measure an orbital period of 19.448 +- 0.002 days. In contrast to recent results, our data on IRS 16SW show an asymmetric phased light curve with a much steeper fall-time than rise-time, which may be due to tidal deformations caused by the proximity of the stars in their orbits. We also identify a possible wind colliding binary (IRS 29N) based on its photometric variation over a few year time-scale which is likely due to episodic dust production. None of the 4 LBV candidates in our sample show the characteristic large increase or decrease in luminosity, however, our time baseline is too short to rule them out as LBVs. Among the remaining variable stars, the majority are early-type stars and three are possibly variable due to line of sight extinction variations. For the 7 OB stars at the center of our field of view that have well-determined 3-dimensional orbits, we see no evidence of flares or dimming of their light, which limits the possibility of a cold, geometrically-thin inactive accretion disk around the supermassive black hole, Sgr A*.

We first refine the fixed concept in the literature that the usage of the Newtonian potential in studies of black hole accretion is invalid and the general relativistic effect must be considered. As our main results, we then show that the energy released by neutrino annihilation in neutrino-dominated accretion flows is sufficient for gamma-ray bursts when the contribution from the optically thick region of the flow is included, and that in the optically thick region advection does not necessarily dominate over neutrino cooling because the advection factor is relevant to the geometrical depth rather than the optical depth of the flow.

We show that the observational data of extragalactic radio sources tend to support the theoretical relationship between the jet precession period and the optical luminosity of the sources, as predicted by the model in which an accretion disk causes the central black hole to precess.

We investigate the linear stability of a shocked accretion flow onto a black hole in the adiabatic limit. Our linear analyses and numerical calculations show that, despite the post-shock deceleration, the shock is generally unstable to non-axisymmetric perturbations. The simulation results of Molteni, T\'oth & Kuznetsov can be well explained by our linear eigenmodes. The mechanism of this instability is confirmed to be based on the cycle of acoustic waves between the corotation radius and the shock. We obtain an analytical formula to calculate the oscillation period from the physical parameters of the flow. We argue that the quasi-periodic oscillation should be a common phenomenon in accretion flows with angular momentum.

(Abridged) We present the first Laser Guide Star Adaptive Optics (LGS-AO) observations of the Galactic center. LGS-AO has dramatically improved the quality, robustness, and versatility with which high angular resolution infrared images of the Galactic center can be obtained with the W. M. Keck II 10-meter telescope. Specifically, Strehl ratios of 0.7 and 0.3 at L'[3.8 micron] and K'[2.1 micron], respectively, are achieved in these LGS-AO images. During our observations, the infrared counterpart to the central supermassive black hole, Sgr A*-IR, showed significant infrared intensity variations, with observed L' magnitudes ranging from 12.6 to 14.5 mag. The faintest end of our L' detections, 1.3 mJy (dereddened), is the lowest level of emission yet observed for this source by a factor of 3. No significant variation in the location of SgrA*-IR is detected as a function of either wavelength or intensity. Near a peak in its intensity, we obtained the first measurement of SgrA*-IR's K'-L' color (3.0 +- 0.2 mag, observed), which corresponds to an intrinsic spectral index of -0.5 +- 0.3. This is significantly bluer than other recent infrared measurements. Because our measurement was taken at a time when Sgr A* was ~6 times brighter in the infrared than the other measurements, we posit that the spectral index of the emission arising from the vicinity of our Galaxy's central black hole may depend on the strength of the flare, with stronger flares giving rise to a higher fraction of high energy electrons in the emitting region.

One of the most perplexing problems associated with the supermassive black hole at the center of our Galaxy is the origin of the young stars in its close vicinity. Using proper motion measurements and stellar number density counts based on 9 years of diffraction-limited K(2.2 micron)-band speckle imaging at the W. M. Keck 10-meter telescopes, we have identified a new comoving group of stars, which we call the IRS 16SW comoving group, located 1.9" (0.08 pc, in projection) from the central black hole. Four of the five members of this comoving group have been spectroscopically identified as massive young stars, specifically He I emission-line stars and OBN stars. This is the second young comoving group within the central parsec of the Milky Way to be recognized and is the closest, by a factor of 2, in projection to the central black hole. These comoving groups may be the surviving cores of massive infalling star clusters that are undergoing disruption in the strong tidal field of the central supermassive black hole.

Recent X-ray and radio observations have identified a transient low-mass X-ray binary (LMXB) located only 0.1 pc in projection from the Galactic center, CXOGC J174540.0-290031. In this paper, we report the detailed analysis of X-ray and infrared observations of the transient and its surroundings. Chandra bservations detect the source at a flux of F_X = 2e-12 erg cm^-2 s^-1 (2-8 keV). After accounting for absorption both in the interstellar medium and in material local to the source, the implied luminosity of the source is only L_X = 4e34 erg/s (2-8 keV; D=8 kpc). However, the diffuse X-ray emission near the source also brightened by a factor of 2. The enhanced diffuse X-ray emission lies on top of a known ridge of dust and ionized gas that is visible infrared images. We interpret the X-ray emission as scattered flux from the outburst, and determine that the peak luminosity of CXOGC J174540.0-290031 was >2e36 erg/s. We suggest that the relatively small observed flux results from the fact that the system is observed nearly edge-on, so that the accretion disk intercepts most of the flux emitted along our line of sight. We compare the inferred peak X-ray luminosity to that of the radio jet. The ratio of the X-ray to radio luminosities, L_X/L_R<1e4, is considerably smaller than in other known LMXBs (> 1e5). This is probably because the jets are radiating with unusually high efficiency at the point where they impact the surrounding interstellar medium. This hypothesis is supported by a comparison with mid-infrared images of the surrounding dust. Finally, we find that the minimum power required to produce the jet, L_jet~1e37 erg/s, is comparable to the inferred peak X-ray luminosity. This is the most direct evidence yet obtained that LMXBs accreting at low rates release about half of their energy as jets.

During five years of Chandra observations, we have identified seven X-ray transients located within 23 pc of Sgr A*. These sources each vary in luminosity by more than a factor of 10, and have peak X-ray luminosities greater than 5e33 erg/s, which strongly suggests that they are accreting black holes or neutron stars. The peak luminosities of the transients are intermediate between those typically considered outburst and quiescence for X-ray binaries. Remarkably four of these transients lie within only 1 pc of Sgr A*. This implies that, compared to the numbers of similar systems located between 1 and 23 pc, transients are over-abundant by a factor of 20 per unit stellar mass within 1 pc of Sgr A*. It is likely that the excess transient X-ray sources are low-mass X-ray binaries that were produced, as in the cores of globular clusters, by three-body interactions between binary star systems and either black holes or neutron stars that have been concentrated in the central parsec through dynamical friction. Alternatively, they could be high-mass X-ray binaries that formed among the young stars that are present in the central parsec.

We present new proper motion measurements and simultaneous orbital solutions for three newly identified (S0-16, S0-19, and S0-20) and four previously known (S0-1, S0-2, S0-4, and S0-5) stars at the Galactic Center. This analysis pinpoints the Galaxy's central dark mass to within +-1 milli-arcsec and, for the first time from orbital dynamics, limits its proper motion to 1.5+-0.5 mas/y, which is consistent with our derivation of the position of Sgr A* in the infrared reference frame (+-10 mas). The estimated central dark mass from orbital motions is 3.7 (+-0.2) x 10^6 (Ro/8kpc)^3 Mo; this is a more direct measure of mass than those obtained from velocity dispersion measurements, which are as much as a factor of two smaller. The smallest closest approach is achieved by S0-16, which confines the mass to within a radius of a mere 45 AU and increases the inferred dark mass density by four orders of magnitude compared to earlier analyses based on velocity and acceleration vectors, making the Milky Way the strongest existing case by far for a supermassive black hole at the center of any normal type galaxy. The stellar orbital properties suggest that the distributions of eccentricities and angular momentum vector and apoapse directions are consistent with those of an isotropic system. Therefore many of the mechanisms proposed for the formation of young stars in the vicinity of a supermassive black hole, such as formation from a pre-existing disk, are unlikely solutions for the Sgr A* cluster stars. Unfortunately, all existing alternative theories are also somewhat problematic. Understanding the apparent youth of stars in the Sgr A* cluster, as well as the more distant He I emission line stars, has now become one of the major outstanding issues in the study of the Galactic Center.

We obtain global solutions of radiatively inefficiently accretion flows around black holes. Whether and where convection develops in a flow are self-consistently determined with the mixing-length theory. The solutions can be divided into three types according to the strength of normal viscosity. Type I solution corresponds to large viscosity parameter $\alpha \ga 0.1$, which is purely advection-dominated and with no convection, and has been extensively studied in the literature. Type II solution is for moderate $\alpha \sim 0.01$, which has a three-zone structure. The inner zone is advection-dominated, the middle zone is convection-dominated and ranges from a few tens to a few thousands of gravitational radii, and the outer zone is convectively stable and matches outward a Keplerian disc. The net energy flux throughout the flow is inward as in type I solution. Type III solution which is for small $\alpha \la 0.001$ consists of two zones as Abramowicz et al. suggested previously: an inner advection-dominated zone and an outer convection-dominated zone, separated at a radius of a few tens of gravitational radii. This type of solution has an outward net energy flux. In both type II and III solutions the radial density profile is between the 1/2 law of self-similar convection-dominated accretion flow model and the 3/2 law of self-similar advection-dominated accretion flow model, and the efficiency of energy release is found to be extremely low. Our results are in good agreement with those of recent numerical simulations.

We use the standard Runge-Kutta method to solve the set of basic equations describing black hole accretion flows composed of two-temperature plasma. We do not invoke any extra energy transport mechanism such as thermal conduction and do not specify any ad hoc outer boundary condition for the advection-dominated accretion flow (ADAF) solution. We find that in the case of high viscosity and non-zero radiative cooling, the ADAF solution can have an asymptotic approach to the Shakura-Sunyaev disk (SSD) solution, and the SSD-ADAF transition radius is close to the central black hole. Our results further prove the mechanism of thermal instability-triggered SSD-ADAF transition suggested previously by Takeuchi & Mineshige and Gu & Lu.

We show that, unlike the results presented previously in the literature, the transition from an outer Shakura-Sunyaev disk (SSD) to an advection-dominated accretion flow (ADAF) is possible for large values of the viscosity parameter $\alpha >0.5$. The transition is triggered by thermal instability of a radiation-pressure-supported SSD. The transition radius is close to the central black hole. We confirm our qualitative prediction by actually constructing global bimodal SSD-ADAF solutions.

The influence of the outer boundary condition (OBC) on the dynamics and radiation of optically thin accretion flow is investigated. Bremsstrahlung and synchrotron radiations amplified by Comptonization are taken into account and two-temperature plasma assumption is adopted. The three OBCs we adopted are the temperatures of the electrons and ions and the specific angular momentum of the accretion flow at a certain outer boundary. We find that when the general parameters such as the mass accretion rate and the viscous parameter are fixed, the peak flux at various bands such as radio, IR and X-ray, can differ by as large as several orders of magnitude under different OBCs in our example. Our results indicate that OBC is both dynamically and radiatively important therefore should be regarded as a new ``parameter'' in accretion disk models. We apply the above results to the compact radio source Sgr A* and find that the discrepancy between the mass accretion rate favored by ADAF models in the literature and that favored by the three dimensional hydrodynamical simulation can be naturally resolved by seriously considering the outer boundary condition of the accretion flow.

We numerically solve the set of dynamical equations describing advection-dominated accretion flows (ADAF) around black holes, using a method similar to that of Chakrabarti (1996a). We choose the sonic radius of the flow $R_s$ and the integration constant in angular momentum equation $j$ as free parameters, and integrate the equations from the sonic point inwards to see if the solution can extend supersonically to the black hole horizon, and outwards to see if and where an acceptable outer boundary of the flow can be found. We recover the ADAF-thin disk solution constructed in Narayan, Kato, & Honma (1997, NKH97), an representative paper of the previous works on global ADAF solutions, although in that paper an apparently very different procedure was adopted. We obtain a complete picture in the form of $R_s-j$ parameter space which sums up the situation of ADAF solution at a glance. For comparison we also present the distribution of global solutions for inviscid flows in the $R_s-j$ space, which supports the view that there should be some similarities between the dynamical behavior of ADAF and that of adiabatic flows, and that there should be a continuous change from the properties of viscous flows to those of inviscid ones.