The unique combination of superb spatial resolution, wide-area coverage, and deep depth of the optical imaging from the Hyper Suprime-Cam (HSC) Subaru Strategic Program is utilized to search for dual quasar candidates. Using an automated image analysis routine on 34,476 known SDSS quasars, we identify those with two (or more) distinct optical point sources in HSC images covering 796 deg^2. We find 421 candidates out to a redshift of 4.5 of which one hundred or so are more likely after filtering out contaminating stars. Angular separations of 0.6 - 4.0" correspond to projected separations of 3 - 30 kpc, a range relatively unexplored for population studies of luminous dual quasars. Using Keck-I/LRIS and Gemini-N/NIFS, we spectroscopically confirm three dual quasar systems at z < 1, two of which are previously unknown out of eight observed, based on the presence of characteristic broad emission lines in each component, while highlighting that the continuum of one object in one of the pairs is reddened. In all cases, the [OIII]5007 emission lines have mild velocity offsets, thus the joint [OIII] line profile is not double-peaked. We find a dual quasar fraction of 0.26+/-0.18% and no evidence for evolution. A comparison with the Horizon-AGN simulation seems to support the case of no evolution in the dual quasar fraction when broadly matching the quasar selection. These results may indicate a scenario in which the frequency of the simultaneous triggering of luminous quasars is not as sensitive as expected to the cosmic evolution of the merger rate or gas content of galaxies.

Astronomical observation data require long-term preservation, and the rapid accumulation of observation data makes it necessary to consider the cost of long-term archive storage. In addition to low-speed disk-based online storage, optical disk or tape-based offline storage can be used to save costs. However, for astronomical research that requires historical data (particularly time-domain astronomy), the performance and energy consumption of data-accessing techniques cause problems because the requested data (which are organized according to observation time) may be located across multiple storage devices. In this study, we design and develop a tool referred to as AstroLayout to redistribute the observation data using spatial aggregation. The core algorithm uses graph partitioning to generate an optimized data placement according to the original observation data statistics and the target storage system. For the given observation data, AstroLayout can copy the long-term archive in the target storage system in accordance with this placement. An efficiency evaluation shows that AstroLayout can reduce the number of devices activated when responding to data-access requests in time-domain astronomy research. In addition to improving the performance of data-accessing techniques, AstroLayout can also reduce the storage systems power consumption. For enhanced adaptability, it supports storage systems of any media, including optical disks, tapes, and hard disks.

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

The post-merger gravitational wave radiation of the remnant formed in the double neutron star coalescence has not been directly measured, yet. We show in this work that the properties of the binary neutron stars involved in GW170817, additionally constrained by PSR J0030+0451, the lower limit on the maximum gravitational mass of non-rotating neutron star and some nuclear data, are in favor of strong post-merger gravitational wave radiation. This conclusion applies to the mergers of Galactic binary neutron star systems as well. Significant post-merger gravitational wave radiation is also preferred to improve the consistency between the maximum gravitational mass of the non-rotating neutron star inferred from GW170817/GRB170817A/AT2017gfo and the latest mass measurements of pulsars. The prominent post-merger gravitational radiation of GW170817-like events are expected to be detectable by advanced LIGO/Virgo detectors in the next decade and then shed valuable lights on the properties of the matter in the extremely high density.

Very recently the NICER collaboration has published the first-ever accurate measurement of mass and radius together for PSR J0030+0451, a nearby isolated quickly-rotating neutron star (NS). In this work we set the joint constraints on the equation of state (EoS) and some bulk properties of NSs with the data of PSR J0030+0451, GW170817 and some nuclear experiments. The piecewise polytropic expansion method and the spectral decomposition method have been adopted to parameterize the EoS. The resulting constraints are consistent with each other. Assuming the maximal gravitational mass of non-rotating NS $M_{\rm TOV}$ lies between $2.04 \rm M_{\odot}$ and $2.4 \rm M_{\odot}$, with the piecewise method the pressure at twice nuclear saturation density is measured to be $3.38^{+2.43}_{-1.50}\times 10^{34}~{\rm dyn~cm^{-2}}$ at the $90\%$ level. For a NS with canonical mass of $1.4 \rm M_\odot$, we have the moment of inertia $I_{1.4} = {1.43}^{+0.28}_{-0.13} \times 10^{38}~{\rm kg \cdot m^2}$, tidal deformability $\Lambda_{1.4} = 390_{-140}^{+320}$, radius $R_{1.4} = 12.2_{-0.9}^{+1.0}~{\rm km}$, and binding energy $BE_{1.4} = {0.16}^{+0.01}_{-0.02} \rm M_{\odot}$ at the $90\%$ level, which are improved in comparison to the constraints with the sole data of GW170817.

Assuming that the differential rotation of the massive neutron star (NS) formed in the double NS (DNS) mergers has been effectively terminated by the magnetic braking and a uniform rotation has been subsequently established (i.e., a supramassive NS is formed), we analytically derive in this work an approximated expression for the critical total gravitational mass ($M_{\rm tot,c}$) to form supramassive NS (SMNS) in the DNS mergers, benefited from some equation of state (EoS) insensitive relationships. The maximum gravitational mass of the nonrotating NSs ($M_{\rm TOV}$) as well as the dimensionless angular momentum of the remnant ($j$) play the dominant roles in modifying $M_{\rm tot,c}$, while the radius and mass differences of the premerger NSs do not. The GW170817/GRB 170817A/AT2017gfo observations have provided so far the best opportunity to quantitatively evaluate $M_{\rm TOV}$. Supposing the central engine for GRB 170817A is a black hole quickly formed in the collapse of an SMNS, we find $M_{\rm TOV}=2.13^{+0.09}_{-0.08}M_\odot$ (68.3% credibility interval, including also the uncertainties of the EoS insensitive relationships), which is consistent with the constraints set by current NS mass measurements.

The discovery of gravitational waves from compact objects coalescence opens a brand-new window to observe the universe. With more events being detected in the future, statistical examinations would be essential to better understand the underlying astrophysical processes. In this work we investigate the prospect of measuring the mass function of black holes that are merging with the neutron stars. Applying Bayesian parameter estimation for hundreds of simulated neutron star$-$black hole (NSBH) mergers, we find that the parameters for most of the injected events can be well recovered. We also take a Bayesian hierarchical model to reconstruct the population properties of the masses of black holes, in the presence of a low mass gap, both the mass gap and power-law index ($\alpha$) of black hole mass function can be well measured, thus we can reveal where the $\alpha$ is different for binary black hole (BBH) and NSBH systems. In the absence of a low mass gap, the gravitational wave data as well as the electromagnetic data can be used to pin down the nature of the merger event and then measure the mass of these very light black holes. However, as a result of the misclassification of BBH into NSBH, the measurement of $\alpha$ is more challenging and further dedicated efforts are needed.

GW190425 is the second neutron star merger event detected by the Advanced LIGO/Virgo detectors. If interpreted as a double neutron star merger, the total gravitational mass is substantially larger than that of the binary systems identified in the Galaxy. In this work we analyze the gravitational-wave data within the neutron star$-$black hole merger scenario. For the black hole, we yield a mass of $2.40^{+0.36}_{-0.32}M_\odot$ and an aligned spin of $0.141^{+0.067}_{-0.064}$. As for the neutron star we find a mass of $1.15^{+0.15}_{-0.13}M_\odot$ and the dimensionless tidal deformability of $1.4^{+3.8}_{-1.2}\times 10^{3}$. These parameter ranges are for 90\% credibility. The inferred masses of the neutron star and the black hole are not in tension with current observations and we suggest that GW190425 is a viable candidate of a neutron star$-$black hole merger event. Benefitting from the continual enhancement of the sensitivities of the advanced gravitational detectors and the increase of the number of the observatories, similar events are anticipated to be much more precisely measured in the future and the presence of black holes below the so-called mass gap will be unambiguously clarified. If confirmed, the mergers of neutron stars with (quickly rotating) low-mass black holes are likely important production sites of the heaviest r-process elements.

For some neutron stars (NSs) in the binary systems, the masses have been accurately measured. While for the isolated neutron stars (INSs), no mass measurement has been reported yet. The situation will change soon thanks to the successful performance of the Neutron Star Interior Composition Explorer (NICER), with which the radius and mass of the isolated PSR J0030+0451 can be simultaneously measured. For most INSs, no mass measurements are possible for NICER because of observational limitations. Benefiting from recent significant progress made on constraining the equation of state of NSs, in this work we propose a way to estimate the masses of the INSs with the measured gravitational redshifts. We apply our method to RX J1856.5-3754, RX J0720.4-3125, and RBS 1223, three members of "The Magnificent Seven" (M7), and estimate their masses to be $1.24_{-0.29}^{+0.29}M_{\odot}$, $1.23_{-0.05}^{+0.10}M_{\odot}$, and $1.08_{-0.11}^{+0.20}M_{\odot}$, respectively. These masses are consistent with that of binary NS systems, suggesting no evidence for experiencing significant accretion of these isolated objects.

In this work we re-analyze the data of GW170817, the first binary neutron star (BNS) merger event, in two ways, including the parameterized equation of state (EoS) method and gravitational wave (GW) parameter estimation analysis. Three sets of data, including the combined LIGO/Virgo (HLV) data, the LIGO-Hanford (H1) data and the LIGO-Livingston (L1) data, have been analyzed, respectively. We confirm the bimodal probability distribution of the tidal deformability of GW170817 with the HLV data and resolve the different contributions of H1 and L1 detectors. Our simulation reveals a tendency of ``overestimating" the tidal parameter with the decreasing signal-to-noise ratio (SNR). Such a tendency is likely caused by the dependence of the result on the prior (which is not well understood and has been assumed to be widely-distributed) in the low SNR case. Though this effect is interesting, the slightly higher SNR of H1 than L1 at frequencies above $\sim 100$ Hz is most-likely not the main reason for the lower tidal parameter region favored by H1. In view of the large fluctuation of the expected tide parameter in the case of low SNR, the different probability distributions favored by the L1 and H1 data of GW170817 are reasonable. We have also explored the dependence of expected difference of SNRs of H1 and L1 detectors on the source locations in the O3 and full sensitivity runs. Such maps are helpful in evaluating the relative powers of individual detectors on measuring the tidal deformabilities for the new double neutron star merger events.

We present the stellar population, using {\it Gaia}\,DR2 parallax, kinematics, and photometry, of the young ($\sim 100$~Myr), nearby ($\sim 230$~pc) open cluster, Blanco1. A total of 644 member candidates are identified via the unsupervised machine learning method \textsc{StarGO} to find the clustering in the 5-dimensional position and proper motion parameter ($X$, $Y$, $Z$, $\mu_\alpha \cos\delta$, $\mu_\delta$) space. Within the tidal radius of $10.0 \pm 0.3$~pc, there are 488 member candidates, 3 times more than those outside. A leading tail and a trailing tail, each of 50--60~pc in the Galactic plane, are found for the first time for this cluster, with stars further from the cluster center streaming away faster, manifest stellar stripping. Blanco1 has a total detected mass of $285\pm32$~M$_\odot $ with a mass function consistent with a slope of $\alpha=1.35\pm0.2$ in the sense of $dN/dm \propto m^{-\alpha}$, in the mass range of 0.25--2.51~M$_\odot $, where $N$ is the number of members and $m$ is stellar mass. A Minimum Spanning Tree ($\Lambda_{\rm MSR}$) analysis shows the cluster to be moderately mass segregated among the most massive members ($\gtrsim 1.4$~M$_\odot$), suggesting an early stage of dynamical disintegration.

In this work we parameterize the Equation of State of dense neutron star (NS) matter with four pressure parameters of $\{\hat{p}_1, \hat{p}_2, \hat{p}_3, \hat{p}_4\}$ and then set the combined constraints with the data of GW 170817 and the data of 6 Low Mass X-ray Binaries (LMXBs) with thermonuclear burst or alternatively the symmetry energy of the nuclear interaction. We find that the nuclear data effectively narrow down the possible range of $\hat{p}_1$, the gravitational wave data plays the leading role in bounding $\hat{p}_2$, and the LMXB data as well as the lower bound on maximal gravitational mass of non-rotating NSs govern the constraints on $\hat{p}_3$ and $\hat{p}_4$. Using posterior samples of pressure parameters and some universal relations, we further investigate how the current data sets can advance our understanding of tidal deformability ($\Lambda$), moment of inertia ($I$) and binding energy ($BE$) of NSs. For a canonical mass of $1.4M_\odot$, we have $I_{1.4} = {1.43}^{+0.30}_{-0.13} \times 10^{38}~{\rm kg \cdot m^2}$, $\Lambda_{1.4} = 390_{-210}^{+280}$ , $R_{1.4} = 11.8_{-0.7}^{+1.2}~{\rm km}$ and $BE_{1.4} = {0.16}^{+0.01}_{-0.02} M_{\odot}$ if the constraints from the nuclear data and the gravitational wave data have been jointly applied. For the joint analysis of gravitational wave data and the LMXB data, we have $I_{1.4} = {1.28}^{+0.15}_{-0.08} \times 10^{38}~{\rm kg \cdot m^2}$, $\Lambda_{1.4} = 220_{-90}^{+90}$, $R_{1.4} = 11.1_{-0.6}^{+0.7}~{\rm km}$ and $BE_{1.4} = {0.18}^{+0.01}_{-0.01} M_{\odot}$. These results suggest that the current constraints on $\Lambda$ and $R$ still suffer from significant systematic uncertainties, while $I_{1.4}$ and $BE_{1.4}$ are better constrained.

The Gamma-Ray Integrated Detectors (GRID) is a space mission concept dedicated to monitoring the transient gamma-ray sky in the energy range from 10 keV to 2 MeV using scintillation detectors onboard CubeSats in low Earth orbits. The primary targets of GRID are the gamma-ray bursts (GRBs) in the local universe. The scientific goal of GRID is, in synergy with ground-based gravitational wave (GW) detectors such as LIGO and VIRGO, to accumulate a sample of GRBs associated with the merger of two compact stars and study jets and related physics of those objects. It also involves observing and studying other gamma-ray transients such as long GRBs, soft gamma-ray repeaters, terrestrial gamma-ray flashes, and solar flares. With multiple CubeSats in various orbits, GRID is unaffected by the Earth occultation and serves as a full-time and all-sky monitor. Assuming a horizon of 200 Mpc for ground-based GW detectors, we expect to see a few associated GW-GRB events per year. With about 10 CubeSats in operation, GRID is capable of localizing a faint GRB like 170817A with a 90% error radius of about 10 degrees, through triangulation and flux modulation. GRID is proposed and developed by students, with considerable contribution from undergraduate students, and will remain operated as a student project in the future. The current GRID collaboration involves more than 20 institutes and keeps growing. On August 29th, the first GRID detector onboard a CubeSat was launched into a Sun-synchronous orbit and is currently under test.

As revealed recently by the modeling of the multi-wavelength data of the emission following GW170817/GRB 170817A, there was an off-axis energetic relativistic outflow component launched by this historic double neutron star merger event. In this work we use the results of these modeling to examine the energy extraction process of the central engine. We show that the magnetic process (i.e., the Blandford-Znajek mechanism) is favored, while the neutrino process usually requires a too massive accretion disk if the duration of the central engine activity is comparable to the observed $T_{90}$ of GRB 170817A, unless the timescale of the central engine activity is less than $\sim$ 0.2s. We propose that the GRB observations are helpful to constrain the combined tidal parameter $\tilde{\Lambda}$, and by adopting the accretion disk mass distribution estimated in BZ mechanism, the $90\%$ credible interval of $\tilde{\Lambda}$ for the progenitor of GW170817 is inferred as $309-954$.

We report the discovery of tidal structures around the intermediate-aged ($\sim$ 700--800~Myr), nearby ($\sim85$~pc) star cluster Coma Berenices. The spatial and kinematic grouping of stars is determined with the {\it Gaia} DR2 parallax and proper motion data, by a clustering analysis tool, \textsc{StarGO}, to map 5D parameters ($X, Y, Z$, $\mu_\alpha \cos\delta, \mu_\delta$) onto a 2D neural network. A leading and a trailing tails, each with an extension of $\sim50$~pc are revealed for the first time around this disrupting star cluster. The cluster members, totaling to $\sim115^{+5}_{-3}\,\rm {M_\odot}$, are clearly mass segregated, and exhibit a flat mass function with $\alpha \sim 0.79\pm0.16$, in the sense of $dN/dm \propto m^{-\alpha}$, where $N$ is the number of member stars and $m$ is stellar mass, in the mass range of $m=0.25$--$2.51~{\rm M_\odot}$. Within the tidal radius of $\sim$6.9~pc, there are 77 member candidates with an average position, i.e., as the cluster center, of R.A.= 186.8110~deg, and decl.= 25.8112~deg, and an average distance of 85.8~pc. Additional 120 member candidates reside in the tidal structures, i.e., outnumbering those in the cluster core. The expansion of escaping members lead to an anisotropy in the velocity field of the tidal tails. Our analysis also serendipitously uncovers an adjacent stellar group, part of which has been cataloged in the literature. We identify 218 member candidates, 10 times more than previously known. This star group is some 65~pc away from, and $\sim400$~Myr younger than, Coma Ber, but is already at the final stage of disruption.

We report the outbursting source PTF 14jg, which prior to the onset of its late 2013 eruption, was a faint, unstudied, and virtually uncatalogued star. The salient features of the PTF 14jg outburst are: (i) projected location near the W4 HII region and radial velocity consistent with physical association; (ii) a lightcurve that underwent a $\sim$6-7 mag optical (R-band) through mid-infrared (L-band) brightening on a few month time scale, that peaked and then faded by $\sim$3 mag, but plateaued still $>$3.5 mag above quiescence by $\sim$8 months post-peak, lasting to at least five years after eruption; (iii) strong outflow signatures, with velocities reaching -530 km/s; (iv) a low gravity and broad ($\sim$100-150 km/s FWHM) optical absorption line spectrum that systematically changes its spectral type with wavelength; (v) lithium; and (vi) ultraviolet and infrared excess. We tentatively identify the outburst as exhibiting characteristics of a young star FU Ori event. However, the burst would be unusually hot, with an absorption spectrum exhibiting high-excitation ($\sim$11,000-15,000 K) lines in the optical, and no evidence of CO in the near-infrared, in addition to exhibiting an unusual lightcurve. We thus also consider alternative scenarios including various forms of novae, nuclear burning instabilities, massive star events, and mergers -- finding them all inferior to the atypically hot FU Ori star classification. The source eventually may be interpreted as a new category of young star outburst with larger amplitude and shorter rise time than most FU Ori-like events. Continued monitoring of the lightcurve and spectral evolution will reveal its true nature.

We have identified stellar and substellar members in the nearby star cluster Coma Berenices, using photometry, proper motions, and distances of a combination of 2MASS, UKIDSS, URAT1, and {\it Gaia}/DR2 data. Those with {\it Gaia}/DR2 parallax measurements provide the most reliable sample to constrain the distance, averaging 86.7~pc with a dispersion 7.1~pc, and the age $\sim800$~Myr, of the cluster. This age is older than the 400--600~Myr commonly adopted in the literature. Our analysis, complete within 5deg of the cluster radius, leads to identification of 192 candidates, among which, after field contamination is considered, about 148 are true members. The members have $J\sim3$~mag to $\sim17.5$~mag, corresponding to stellar masses 2.3--0.06~$M_\odot$. The mass function of the cluster peaks around 0.3~$M_\odot$ and, in the sense of $dN/dm = m^{-\alpha}$, where $N$ is the number of members and $m$ is stellar mass, has a slope $\alpha\approx 0.49\pm0.03$ in the mass range 0.3--2.3~$M_\odot$. This is much shallower than that of the field population in the solar neighborhood. The slope $\alpha=-1.69\pm0.14$ from 0.3~$M_\odot$ to 0.06~$M_\odot$, the lowest mass in our sample. The cluster is mass segregated and has a shape elongated toward the Galactic plane. Our list contains nine substellar members, including three new discoveries of an M8, an L1 and an L4 brown dwarfs, extending from the previously known coolest members of late-M types to even cooler types.

We present a sample of 59 periodic variables from the Palomar Transient Factory, selected from published catalogues of white dwarf (WD) candidates. The variability can likely be attributed to ellipsoidal variation of the tidally distorted companion induced by the gravity of the primary (WD or hot subdwarf) or to the reflection of hot emission by a cooler companion. We searched 11311 spectroscopically or photometrically selected WD candidates from three hot star/WD catalogues, using the Lomb-Scargle periodogram to single out promising sources. We present period estimates for the candidates, 45 of which were not previously identified as periodic variables, and find that most have a period shorter than a few days. Additionally, we discuss the eclipsing systems in our sample and present spectroscopic data on selected sources.

We present TYC 2505-672-1 as a newly discovered and remarkable eclipsing system comprising an M-type red giant that undergoes a ~3.45 year long, near-total eclipse (depth of ~4.5 mag) with a very long period of ~69.1 yr. TYC 2505-672-1 is now the longest-period eclipsing binary system yet discovered, more than twice as long as that of the currently longest-period system, $\epsilon$ Aurigae. We show from analysis of the light curve including both our own data and historical data spanning more than 120 yr and from modeling of the spectral energy distribution, both before and during eclipse, that the red giant primary is orbited by a moderately hot source (T$_{eff}$~8000 K) that is itself surrounded by an extended, opaque circumstellar disk. From the measured ratio of luminosities, the radius of the hot companion must be in the range 0.1-0.5 Rsun (depending on the assumed radius of the red giant primary), which is an order of magnitude smaller than that for a main sequence A star and 1-2 orders of magnitude larger than that for a white dwarf. The companion is therefore most likely a "stripped red giant" subdwarf-B type star destined to become a He white dwarf. It is however somewhat cooler than most sdB stars, implying a very low mass for this "pre-He-WD" star. The opaque disk surrounding this hot source may be a remnant of the stripping of its former hydrogen envelope. However, it is puzzling how this object became stripped, given that it is at present so distant (orbital semi-major axis of ~24 AU) from the current red giant primary star. Extrapolating from our calculated ephemeris, the next eclipse should begin in early UT 2080 April and end in mid UT 2083 September (eclipse center UT 2081 December 24). This system is poised to become an exemplar of a very rare class of systems, even more extreme in several respects than the well studied archetype $\epsilon$ Aurigae.

PSR J2129-0429 is a "redback" eclipsing millisecond pulsar binary with an unusually long 15.2 hour orbit. It was discovered by the Green Bank Telescope in a targeted search of unidentified Fermi gamma-ray sources. The pulsar companion is optically bright (mean $m_R = 16.6$ mag), allowing us to construct the longest baseline photometric dataset available for such a system. We present ten years of archival and new photometry of the companion from LINEAR, CRTS, PTF, the Palomar 60-inch, and LCOGT. Radial velocity spectroscopy using the Double-Beam Spectrograph on the Palomar 200-inch indicates that the pulsar is massive: $1.74\pm0.18 M_\odot$. The G-type pulsar companion has mass $0.44\pm0.04 M_\odot$, one of the heaviest known redback companions. It is currently 95\% Roche-lobe filling and only mildly irradiated by the pulsar. We identify a clear 13.1 mmag yr$^{-1}$ secular decline in the mean magnitude of the companion as well as smaller-scale variations in the optical lightcurve shape. This behavior may indicate that the companion is cooling. Binary evolution calculations indicate that PSR J2129-0429 has an orbital period almost exactly at the bifurcation period between systems that converge into tighter orbits as black widows and redbacks and those that diverge into wider pulsar--white dwarf binaries. Its eventual fate may depend on whether it undergoes future episodes of mass transfer and increased irradiation.

We present the results of a Palomar Transient Factory (PTF) archival search for blue transients which lie in the magnitude range between "normal" core-collapse and superluminous supernovae (i.e. with $-21\,{\leq}M_{R\,(peak)}\,{\leq}-19$). Of the six events found after excluding all interacting Type~IIn and Ia-CSM supernovae, three (PTF09ge, 09axc and 09djl) are coincident with the centers of their hosts, one (10iam) is offset from the center, and for two (10nuj and 11glr) a precise offset can not be determined. All the central events have similar rise times to the He-rich tidal disruption candidate PS1-10jh, and the event with the best-sampled light curve also has similar colors and power-law decay. Spectroscopically, PTF09ge is He-rich, while PTF09axc and 09djl display broad hydrogen features around peak magnitude. All three central events are in low star-formation hosts, two of which are E+A galaxies. Our spectrum of the host of PS1-10jh displays similar properties. PTF10iam, the one offset event, is different photometrically and spectroscopically from the central events and its host displays a higher star formation rate. Finding no obvious evidence for ongoing galactic nuclei activity or recent star formation, we conclude that the three central transients likely arise from the tidal disruption of a star by a super-massive black hole. We compare the spectra of these events to tidal disruption candidates from the literature and find that all of these objects can be unified on a continuous scale of spectral properties. The accumulated evidence of this expanded sample strongly supports a tidal disruption origin for this class of nuclear transients.

We report the optical identification of the companion to the {\it Fermi} black widow millisecond pulsar PSR J1544+4937. We find a highly variable source on Keck LRIS images at the nominal pulsar position, with 2 magnitude variations over orbital period in the B, g, R, and I bands. The nearly achromatic light curves are difficult to explain with a simply irradiated hemisphere model, and suggest that the optical emission is dominated by a nearly isothermal hot patch on the surface of the companion facing the pulsar. We roughly constrain the distance to PSR J1544+4937 to be between 2 and 5 kpc. A more reliable distance measurement is needed in order to constrain the composition of the companion.

The iPTF detection of the most recent outburst of the recurrent nova system RX J0045.4+4154 in the Andromeda Galaxy has enabled the unprecedented study of a massive ($M>1.3\ M_\odot$) accreting white dwarf (WD). We detected this nova as part of the near daily iPTF monitoring of M31 to a depth of $R\approx 21$\,mag and triggered optical photometry, spectroscopy and soft X-ray monitoring of the outburst. Peaking at an absolute magnitude of $M_R=-6.6$ mag, and with a decay time of 1 mag per day, it is a faint and very fast nova. It shows optical emission lines of He/N and expansion velocities of 1900 to 2600 km s$^{-1}$ 1--4 days after the optical peak. The {\it Swift} monitoring of the X-ray evolution revealed a supersoft source (SSS) with $kT_{\rm eff}\approx 90-110\ {\rm eV}$ that appeared within 5 days after the optical peak, and lasted only 12 days. Most remarkably, this is not the first event from this system, rather it is a recurrent nova with a time between outbursts of approximately 1 year, the shortest known. Recurrent X-ray emission from this binary was detected by ROSAT in 1992 and 1993, and the source was well characterized as a $M>1.3\ M_\odot$ WD SSS. Based on the observed recurrence time between different outbursts, the duration and effective temperature of the SS phase, MESA models of accreting WDs allow us to constrain the accretion rate to $\dot M>1.7\times10^{-7}\ {M_{\odot}\ {\rm yr}}^{-1}$ and WD mass $>1.30\ M_{\odot}$. If the WD keeps $30\%$ of the accreted material, it will take less than a Myr to reach core densities high enough for carbon ignition (if made of C/O) or electron capture (if made of O/Ne) to end the binary evolution.

Recent observations strongly suggest that the millisecond pulsar binary PSR J1023+0038 has developed an accretion disk since 2013 June. We present a multi-wavelength analysis of PSR J1023+0038, which reveals that 1) its gamma-rays suddenly brightened within a few days in June/July 2013 and has remained at a high gamma-ray state for several months; 2) both UV and X-ray fluxes have increased by roughly an order of magnitude, and 3) the spectral energy distribution has changed significantly after the gamma-ray sudden flux change. Time variabilities associated with UV and X-rays are on the order of 100-500 seconds and 50-100 seconds, respectively. Our model suggests that a newly formed accretion disk due to the sudden increase of the stellar wind could explain the changes of all these observed features. The increase of UV is emitted from the disk, and a new component in gamma-rays is produced by inverse Compton scattering between the new UV component and pulsar wind. The increase of X-rays results from the enhancement of injection pulsar wind energy into the intra-binary shock due to the increase of the stellar wind. We also predict that the radio pulses may be blocked by the evaporated winds from the disk and the pulsar is still powered by rotation.

We report positions, velocities and metallicities of 50 ab-type RR Lyrae (RRab) stars observed in the vicinity of the Orphan stellar stream. Using about 30 RRab stars classified as being likely members of the Orphan stream, we study the metallicity and the spatial extent of the stream. We find that RRab stars in the Orphan stream have a wide range of metallicities, from -1.5 dex to -2.7 dex. The average metallicity of the stream is -2.1 dex, identical to the value obtained by Newberg et al. (2010) using blue horizontal branch stars. We find that the most distant parts of the stream (40-50 kpc from the Sun) are about 0.3 dex more metal-poor than the closer parts (within ~30 kpc), suggesting a possible metallicity gradient along the stream's length. We have extended the previous studies and have mapped the stream up to 55 kpc from the Sun. Even after a careful search, we did not identify any more distant RRab stars that could plausibly be members of the Orphan stream. If confirmed with other tracers, this result would indicate a detection of the end of the leading arm of the stream. We have compared the distances of Orphan stream RRab stars with the best-fit orbits obtained by Newberg et al. (2010). We find that model 6 of Newberg et al. (2010) cannot explain the distances of the most remote Orphan stream RRab stars, and conclude that the best fit to distances of Orphan stream RRab stars and to the local circular velocity is provided by potentials where the total mass of the Galaxy within 60 kpc is M_{60}~2.7x10^{11} Msun, or about 60% of the mass found by previous studies. More extensive modelling that would consider non-spherical potentials and the possibility of misalignment between the stream and the orbit, is highly encouraged.

The Digital Access to a Sky Century@Harvard (DASCH) project is digitizing the ~500,000 glass plate images obtained (full sky) by the Harvard College Observatory from 1885-1992. Astrometry and photometry for each resolved object are derived with photometric rms values of ~0.15mag for the initial photometry analysis pipeline. Here we describe new developments for DASCH photometry, applied to the Kepler field, that has yielded further improvements, including better identification of image blends and plate defects by measuring image profiles and astrometric deviations. A local calibration procedure using nearby stars in a similar magnitude range as the program star (similar to what has been done for visual photometry from the plates) yields additional improvement for a net photometric rms ~0.1mag. We also describe statistical measures of light curves that are now used in the DASCH pipeline processing to autonomously identify new variables. The DASCH photometry methods described here are used in the pipeline processing for the Data Releases of DASCH data (Grindlay et al. 2012; http://dasch.rc.fas.harvard.edu), as well as for the long-term variables discovered by DASCH in the Kepler field (Tang et al, in preparation).

We present 100 year light curves of Kepler planet-candidate host stars from the Digital Access to a Sky Century at Harvard (DASCH) project. 261 out of 997 host stars have at least 10 good measurements on DASCH scans of the Harvard plates. 109 of them have at least 100 good measurements, including 70% (73 out of 104) of all host stars with g<=13 mag, and 44% (100 out of 228) of all host stars with g<=14 mag. Our typical photometric uncertainty is ~0.1-0.15 mag. No variation is found at 3-sigma level for these host stars, including 21 confirmed or candidate hot Jupiter systems which might be expected to show enhanced flares from magnetic interactions between dwarf primaries and their close and relatively massive planet companions.

We report the discovery of R Coronae Borealis (RCB) stars in the Andromeda galaxy (M31) using the Palomar Transient Factory (PTF). RCB stars are rare hydrogen-deficient, carbon-rich supergiant variables, most likely the merger products of two white dwarfs. These new RCBs, including two confirmed ones and two candidates, are the first to be found beyond the Milky Way and the Magellanic Clouds. All of M31 RCBs showed >1.5 mag irregular declines over timescales of weeks to months. Due to the limiting magnitude of our data (R~21-22 mag), these RCB stars have R~19.5 to 20.5 mag at maximum light, corresponding to M_R= -4 to -5, making them some of the most luminous RCBs known. Spectra of two objects show that they are warm RCBs, similar to the Milky Way RCBs RY Sgr and V854 Cen. We consider these results, derived from a pilot study of M31 variables, as an important proof-of-concept for the study of rare bright variables in nearby galaxies with the PTF or other synoptic surveys.

We analyzed the 100-yr light curves of Galactic high-mass X-ray binaries using the Harvard photographic plate collection, made accessible through the DASCH project (Digital Access to a Sky Century at Harvard). As scanning is still in progress, we focus on the four objects that are currently well covered: the supergiant X-ray binary Cyg X-1 (V1357 Cyg), and the Be X-ray binaries 1H 1936+541 (BD+53 2262), RX J1744.7-2713 (HD 161103), and RX J2030.5+4751 (SAO 49725), the latter two objects being similar to gamma Cas. The star associated with Cyg X-1 does not show evidence for variability with an amplitude higher than 0.3 magnitude over a hundred years. We found significant variability of one magnitude with timescales of more than 10 years for SAO 49725, as well as a possible period of 500-600 days and an amplitude of 0.05 magnitude that might be the orbital, or super-orbital period of the system. The data is insufficient to conclude for HD 161103 but suggests a similar long-term variability. We thus observe an additional characteristic of gamma Cas-like objects: their long-term variability. This variability seems to be due to the slow evolution of a decretion disk around the Be star, but may be triggered by the presence of a compact object in the system, possibly a white dwarf. This characteristic could be used to identify further similar objects otherwise difficult to detect.

The large-scale surveys such as PTF, CRTS and Pan-STARRS-1 that have emerged within the past 5 years or so employ digital databases and modern analysis tools to accentuate research into Time Domain Astronomy (TDA). Preparations are underway for LSST which, in another 6 years, will usher in the second decade of modern TDA. By that time the Digital Access to a Sky Century @ Harvard (DASCH) project will have made available to the community the full sky Historical TDA database and digitized images for a century (1890--1990) of coverage. We describe the current DASCH development and some initial results, and outline plans for the "production scanning" phase and data distribution which is to begin in 2012. That will open a 100-year window into temporal astrophysics, revealing rare transients and (especially) astrophysical phenomena that vary on time-scales of a decade. It will also provide context and archival comparisons for the deeper modern surveys

We present multi-band optical photometry of 94 spectroscopically-confirmed Type Ia supernovae (SN Ia) in the redshift range 0.0055 to 0.073, obtained between 2006 and 2011. There are a total of 5522 light curve points. We show that our natural system SN photometry has a precision of roughly 0.03 mag or better in BVr'i', 0.06 mag in u', and 0.07 mag in U for points brighter than 17.5 mag and estimate that it has a systematic uncertainty of 0.014, 0.010, 0.012, 0.014, 0.046, and 0.073 mag in BVr'i'u'U, respectively. Comparisons of our standard system photometry with published SN Ia light curves and comparison stars reveal mean agreement across samples in the range of ~0.00-0.03 mag. We discuss the recent measurements of our telescope-plus-detector throughput by direct monochromatic illumination by Cramer et al (in prep.). This technique measures the whole optical path through the telescope, auxiliary optics, filters, and detector under the same conditions used to make SN measurements. Extremely well-characterized natural-system passbands (both in wavelength and over time) are crucial for the next generation of SN Ia photometry to reach the 0.01 mag accuracy level. The current sample of low-z SN Ia is now sufficiently large to remove most of the statistical sampling error from the dark energy error budget. But pursuing the dark-energy systematic errors by determining highly-accurate detector passbands, combining optical and near-infrared (NIR) photometry and spectra, using the nearby sample to illuminate the population properties of SN Ia, and measuring the local departures from the Hubble flow will benefit from larger, carefully measured nearby samples.

We present photometric and spectroscopic observations of a peculiar variable (designated DASCH J075731.1+201735 or J0757) discovered from our DASCH project using the digitized Harvard College Observatory archival photographic plates. It brightened by about 1.5 magnitudes in B within a year starting in 1942, and then slowly faded back to its pre-outburst brightness from 1943 to the 1950s. The mean brightness level was stable before and after the outburst, and ellipsoidal variations with a period of $P=119.18\pm0.07$ days are seen, suggesting that the star is tidally distorted. Radial-velocity measurements indicate that the orbit is nearly circular ($e=0.02\pm0.01$) with a spectroscopic period that is the same as the photometric period. The binary consists of a $1.1\pm0.3 M_\odot$ M0III star, and a $0.6\pm0.2 M_\odot$ companion, very likely a white dwarf (WD). Unlike other symbiotic binaries, there is no sign of emission lines or a stellar wind in the spectra. With an outburst timescale of ~10 years and estimated B band peak luminosity M_B~0.7, J0757 is different from any other known classic or symbiotic novae. The most probable explanation of the outburst is Hydrogen shell-burning on the WD, although an accretion-powered flare cannot be ruled out.

We report the detection of a planet whose orbit surrounds a pair of low-mass stars. Data from the Kepler spacecraft reveal transits of the planet across both stars, in addition to the mutual eclipses of the stars, giving precise constraints on the absolute dimensions of all three bodies. The planet is comparable to Saturn in mass and size, and is on a nearly circular 229-day orbit around its two parent stars. The eclipsing stars are 20% and 69% as massive as the sun, and have an eccentric 41-day orbit. The motions of all three bodies are confined to within 0.5 degree of a single plane, suggesting that the planet formed within a circumbinary disk.

NSV 11154 has been confirmed as a new member of the rare hydrogen deficient R Coronae Borealis (RCB) stars based on new photometric and spectroscopic data. Using new photometry, as well as archival plates from the Harvard archive, we have constructed the historical lightcurve of NSV 11154 from 1896 to the present. The lightcurve shows the sudden, deep, irregularly spaced declines characteristic of RCB stars. The visible spectrum is typical of a cool (Teff < 5000 K) RCB star showing no hydrogen lines, strong C2 Swan bands, and no evidence of 13C. In addition, the star shows small pulsations typical of an RCB star, and an infrared excess due to circumstellar dust with a temperature of ~800 K. The distance to NSV 11154 is estimated to be ~14.5 kpc. RCB stars are very rare in the Galaxy so each additional star is important to population studies leading to a better understanding the origins of these mysterious stars. Among the known sample of RCB stars, NSV 11154 is unusual in that it lies well above the Galactic plane (5 kpc) and away from the Galactic Center which suggests that its parent population is neither thick disk nor bulge.

KU Cyg is an eclipsing binary consisting of a F-type star accreting through a large accretion disk from a K5III red giant. Here we present the discovery of a 5-yr dip around 1900 found from its 100 yr DASCH (Digital Access to a Sky Century at Harvard) light curve. It showed a ~0.5 mag slow fading from 1899 to 1903, and brightened back around 1904 on a relatively shorter timescale. The light curve shape of the 1899-1904 fading-brightening event differs from the dust production and dispersion process observed in R Coronae Borealis (RCB) stars, which usually has a faster fading and slower recovery, and for KU Cyg is probably related to the accretion disk surrounding the F star. The slow fading in KU Cyg is probably caused by increases in dust extinction in the disk, and the subsequent quick brightening may be due to the evaporation of dust transported inwards through the disk. The extinction excess which caused the fading may arise from increased mass transfer rate in the system, or from dust clump ejections from the K giant.

We describe the process implemented in the DASCH pipeline which applies a reliable astrometric correction to each scanned plate. Our first blind astrometric fit resolves the pointing, scale and orientation of the plate in the sky using astrometry. net code. Then we iteratively improve this solution with WCSTools imwcs. Finally, we apply a 6th order polynomial fit with SCAMP to correct the image for distortions. During a test on 140 plates, this process has allowed us to successfully correct 79% of the plates. With further refinements of the process we now reach a 95% success rate after reprocessing all our scanned plates (~11 000 in Nov. 2010). We could extract a lightcurve for 2.85 times more objects than with the previous Pipeline, down to magnitude 17. The resulting median RMS error is 0.13" for objects with mag. 8 to 17.

Digital Access to a Sky Century @ Harvard (DASCH) is a project to digitize the collection of approximately 525,000 astronomical plates held at the Harvard College Observatory. This paper presents an overview of the DASCH data processing pipeline, with special emphasis on the processing of multiple-exposure plates. Such plates extended the dynamic range of photograph emulsions and improved photometric accuracy by minimizing variations in plate development procedures. Two approaches are explored in this paper: The repetitive use of astrometry.net (Lang et al. 2010) and local correlation searches. Both procedures have yielded additional quality control checks useful to the pipeline.

Based on new and archival Chandra observations of the Sombrero galaxy (M 104), we study the diffuse X-ray emission in and around its massive stellar bulge. The 2-6 keV unresolved emission from the bulge region closely follows the K-band star light and most likely arises from unresolved stellar sources. At lower energies, however, the unresolved emission reaches a galactocentric radius of at least 23 kpc, significantly beyond the extent of the starlight, clearly indicating the presence of diffuse hot gas. We isolate the emission of the gas by properly accounting for the emission from unresolved stellar sources, predominantly cataclysmic variables and coronally active binaries, whose quasi-universal X-ray emissivity is recently established. We find a gas temperature of ~0.6 keV with little variation across the field of view, except for a lower temperature of ~0.3 keV along the stellar disk. We measure a total intrinsic 0.3-2 keV luminosity of ~2e39 erg/s, which is comparable to the prediction by the latest galaxy formation models for disk galaxies as massive as Sombrero. However, such numerical models do not fully account for internal feedback processes, such as nuclear feedback and stellar feedback, against accretion from the intergalactic medium. On the other hand, we find no evidence for either the nucleus or the very modest star-forming activities in the disk to be a dominant heating source for the diffuse gas. We also show that neither the expected energy released by Type Ia supernovae nor the expected mass returned by evolved stars is recovered by observations. We argue that in Sombrero a galactic-scale subsonic outflow of hot gas continuously removes much of the "missing" energy and mass input from the bulge region. The observed density and temperature distributions of such an outflow, however, continues to pose challenges to theoretical studies.

Gamma-ray spectroscopy provides diagnostics of particle acceleration in solar flares, but care must be taken when interpreting the spectra due to effects of the angular distribution of the accelerated particles (such as relativistic beaming) and Compton reprocessing of the radiation in the solar atmosphere. In this paper, we use the GEANT4 Monte Carlo package to simulate the interactions of accelerated electrons and protons and study these effects on the gamma-rays resulting from electron bremsstrahlung and pion decay. We consider the ratio of the 511~keV annihilation-line flux to the continuum at 200~keV and in the energy band just above the nuclear de-excitation lines (8--15~MeV) as a diagnostic of the accelerated particles and a point of comparison with data from the X17 flare of 2003 October 28. We also find that pion secondaries from accelerated protons produce a positron annihilation line component at a depth of $\sim$ 10 g cm$^{-2}$, and that the subsequent Compton scattering of the 511~keV photons produces a continuum that can mimic the spectrum expected from the 3$\gamma$ decay of orthopositronium.

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

We present a study of discrete X-ray sources in and around the bulge-dominated, massive Sa galaxy, Sombrero (M104), based on new and archival Chandra observations with a total exposure of ~200 ks. With a detection limit of L_X = 1E37 erg/s and a field of view covering a galactocentric radius of ~30 kpc (11.5 arcminute), 383 sources are detected. Cross-correlation with Spitler et al.'s catalogue of Sombrero globular clusters (GCs) identified from HST/ACS observations reveals 41 X-rays sources in GCs, presumably low-mass X-ray binaries (LMXBs). We quantify the differential luminosity functions (LFs) for both the detected GC and field LMXBs, whose power-low indices (~1.1 for the GC-LF and ~1.6 for field-LF) are consistent with previous studies for elliptical galaxies. With precise sky positions of the GCs without a detected X-ray source, we further quantify, through a fluctuation analysis, the GC LF at fainter luminosities down to 1E35 erg/s. The derived index rules out a faint-end slope flatter than 1.1 at a 2 sigma significance, contrary to recent findings in several elliptical galaxies and the bulge of M31. On the other hand, the 2-6 keV unresolved emission places a tight constraint on the field LF, implying a flattened index of ~1.0 below 1E37 erg/s. We also detect 101 sources in the halo of Sombrero. The presence of these sources cannot be interpreted as galactic LMXBs whose spatial distribution empirically follows the starlight. Their number is also higher than the expected number of cosmic AGNs (52+/-11 [1 sigma]) whose surface density is constrained by deep X-ray surveys. We suggest that either the cosmic X-ray background is unusually high in the direction of Sombrero, or a distinct population of X-ray sources is present in the halo of Sombrero.

Type Ia supernovae (Ia SNe) provide a rich source of iron for hot gas in galactic stellar spheroids. However, the expected super-solar iron abundance of the hot gas is not observed. Instead, X-ray observations often show decreasing iron abundance toward galactic central regions, where the Ia SN enrichment is expected to be the highest. We examine the cause of this missing iron problem by studying the enrichment process and its effect on X-ray abundance measurements of the hot gas. The evolution of Ia SN iron ejecta is simulated in the context of galaxy-wide hot gas outflows, in both supersonic and subsonic cases, as may be expected for hot gas in galactic bulges or elliptical galaxies of intermediate masses. SN reverse-shock heated iron ejecta is typically found to have a very high temperature and low density, hence producing little X-ray emission. Such hot ejecta, driven by its large buoyancy, can quickly reach a substantially higher outward velocity than the ambient medium, which is dominated by mass loss from evolved stars. The ejecta is gradually and dynamically mixed with the medium at large galactic radii. The ejecta is also slowly diluted and cooled by {\sl insitu} mass injection from evolved stars. These processes together naturally result in the observed positive gradient in the average radial iron abundance distribution of the hot gas, even if mass-weighted. This trend is in addition to the X-ray measurement bias that tends to underestimate the iron abundance for the hot gas with a temperature distribution.

Here we present the discovery of three unusual long-term variables found in the Digital Access to a Sky Century at Harvard (DASCH) project, with ~1 magnitude variations in their lightcurves on ~10-100 yr timescales. They are all spectroscopically identified as K2III giant stars, probably in the thick disk. Their lightcurves do not match any previously measured for known types of variable stars, or any theoretical model reported for red giants, and instead suggest a new dust formation mechanism or the first direct observation of "short" timescale evolution-driven variability. More theoretical work on the lithium flash near the Red Giant Branch (RGB) bump and the helium shell ignition in the lower Asymptotic Giant Branch (AGB), as well as long term monitoring of K2III thick disk stars is needed.

Double-peaked emission lines are believed to be originated from accretion disks around supermassive black holes (SMBHs), and about 3% of z<0.33 AGNs are found to be double-peaked emitters. The quasar SDSS J153636.22+044127.0 has recently been identified with peculiar broad-line emission systems exhibiting multiple redshifts. We decompose the H$\alpha$ and H$\beta$ profiles into a circular Keplerian disk line component and other Gaussian components. We propose that the system is both a double-peaked emitter and a binary SMBH system, where the extra-flux in the blue peaks of the broad lines comes from the region around the secondary black hole. We suggest that such black hole binary systems might also exist in many known double-peaked emitters, where the tidal torques from the secondary black hole clear the outer region of the disk around the primary black hole, similar to the gap in a protostellar disk due to the process of planetary migration, and might also stimulate the formation of a vertical extended source in the inner region around the primary which illuminates the disk. However, most secondary SMBHs in such systems might be too small to maintain a detectable broad line region (BLR), so that the disk line from the primary dominates.

A two-components model for Broad Line Region (BLR) of Active Galactic Nuclei (AGN) has been suggested for many years but not widely accepted (e.g., Hu et al. 2008; Sulentic et al. 2000; Brotherton et al. 1996; Mason et al. 1996). This model indicates that the broad line can be described with superposition of two Gaussian components (Very Broad Gaussian Component (VBGC) and InterMediate Gaussian Component (IMGC)) which are from two physically distinct regions; i.e., Very Broad Line Region (VBLR) and InterMediate Line Region (IMLR). We select a SDSS sample to further confirm this model and give detailed analysis to the geometry, density and evolution of these two regions. Micro-lensing result of BLR in J1131-1231 and some unexplained phenomena in Reverberation Mapping (RM) experiment provide supportive evidence for this model. Our results indicate that the radius obtained from the emission line RM normally corresponds to the radius of the VBLR, and the existence of the IMGC may affect the measurement of the black hole masses in AGNs. The deviation of NLS1s from the M-sigma relation and the Type II AGN fraction as a function of luminosity can be explained in this model in a coherent way. The evolution of the two emission regions may be related to the evolutionary stages of the broad line regions of AGNs from NLS1s to BLS1s. Based on the results presented here, a unified picture of hierarchical evolution of black hole, dust torus and galaxy is proposed.

Using a cross-correlation method, we study the X-ray halo of Cyg X-3. Two components of dust distributions are needed to explain the time lags derived by the cross-correlation method. Assuming the distance as 1.7 kpc for Cygnus OB2 association (a richest OB association in the local Galaxy) and another uniform dust distribution, we get a distance of $7.2^{+0.3}_{-0.5}$ kpc (68$%$ confidence level) for Cyg X-3. When using the distance estimation of Cygnus OB2 as 1.38 kpc or 1.82 kpc, the inferred distance for Cyg X-3 is $3.4^{+0.2}_{-0.2}$ kpc or $9.3^{+0.6}_{-0.4}$ kpc respectively. The distance estimation uncertainty of Cyg X-3 is mainly related to the distance of the Cygnus OB2, which may be improved in the future with high precision astrometric measurements. The advantage of this method is that the result depends weakly on the photon energy, dust grain radius, scattering cross-section and so on.

Stellar feedback in galactic bulges plays an essential role in shaping the evolution of galaxies. To quantify this role and facilitate comparisons with X-ray observations, we conduct 3D hydrodynamical simulations with the adaptive mesh refinement code, FLASH, to investigate the physical properties of hot gas inside a galactic bulge, similar to that of our Galaxy or M31. We assume that the dynamical and thermal properties of the hot gas are dominated by mechanical energy input from SNe, primarily Type Ia, and mass injection from evolved stars as well as iron enrichment from SNe. We study the bulge-wide outflow as well as the SN heating on scales down to ~4 pc. An embedding scheme that is devised to plant individual SNR seeds, allows to examine, for the first time, the effect of sporadic SNe on the density, temperature, and iron ejecta distribution of the hot gas as well as the resultant X-ray morphology and spectrum. We find that the SNe produce a bulge wind with highly filamentary density structures and patchy ejecta. Compared with a 1D spherical wind model, the non-uniformity of simulated gas density, temperature, and metallicity substantially alters the spectral shape and increases the diffuse X-ray luminosity. The differential emission measure as a function of temperature of the simulated gas exhibits a log-normal distribution, with a peak value much lower than that of the corresponding 1D model. The bulk of the X-ray emission comes from the relatively low temperature and low abundance gas shells associated with SN blastwaves. SN ejecta are not well mixed with the ambient medium, at least in the bulge region. These results, at least partly, account for the apparent lack of evidence for iron enrichment in the soft X-ray-emitting gas in galactic bulges and intermediate-mass elliptical galaxies.[...]

Many hydrodynamic processes can be studied in a way that is scalable over a vastly relevant physical parameter space. We systematically examine this scalability, which has so far only briefly discussed in astrophysical literature. We show how the scalability is limited by various constraints imposed by physical processes and initial conditions. Using supernova remnants in different environments and evolutionary phases as application examples, we demonstrate the use of the scaling as a powerful tool to explore the interdependence among relevant parameters, based on a minimum set of simulations. In particular, we devise a scaling scheme that can be used to adaptively generate numerous seed remnants and plant them into 3D hydrodynamic simulations of the supernova-dominated interstellar medium.

We have decomposed the broad H-alpha, H-beta and H-gamma lines of 90 Active Galactic Nuclei (AGNs) into a superposition of a very broad and an intermediate Gaussian components (VBGC and IMGC) and discovered that the two Gaussian components evolve with FWHM of the whole emission lines. We suggest that the VBGC and the IMGC are produced in different emission regions, namely, Very Broad Line Region (VBLR) and Intermediate Line Region (IMLR). The details of the two components of H-alpha, H-beta and H-gamma lines indicate that the radius obtained from the emission line reverberation mapping normally corresponds to the radius of the VBLR, but the radius obtained from the infrared reverberation mapping corresponding to IMLR, i.e., the inner boundary of the dusty torus. The existence of the IMGC may affect the measurement of the black hole mass in AGNs. Therefore, the deviation of NLS1s from the M-sigma relation may be explained naturally in this way. The evolution of the two emission line regions may be related to the evolutionary stages of the broad line regions of AGNs from NLS1s to BLS1s. Other evidences for the existence of the IMLR are also presented.

X-ray photons scattered by the interstellar medium carry information about dust distribution, dust grain model, scattering cross section, and the distance of the source; they also take longer time than unscattered photons to reach the observer. Using a cross-correlation method, we study the light curves of the X-ray dust scattering halo of Cyg X-1, observed with the \textit{Chandra X-ray Observatory}. Significant time lags are found between the light curves of the point source and its halo. This time lag increases with the angular distance from Cyg X-1, implying a dust concentration at a distance along the line of sight (LOS) of 2.0 kpc $\times$ (0.876 $\pm$ 0.002) from the Earth. By fitting the observed light curves of the halo at different radii with simulated light curves, we obtain a width of $\mathit{\Delta L}=33_{-13}^{+18}$ pc of this dust concentration. The origin of this dust concentration is still not clearly known. The advantage of our method is that we need no assumption of scattering cross section, dust grain model, or dust distribution along the LOS. Combining the derived dust distribution from the cross-correlation study with the surface brightness distribution of the halo, we conclude that the two commonly accepted models of dust grain size distribution need to be modified significantly.

Digital Access to a Sky Century @ Harvard (DASCH) is a project to digitize the collection of ~500,000 glass photographic plates held at Harvard College Observatory. The collection spans the time period from 1880 to 1985, during which time every point on the sky has been observed approximately 500 to 1000 times. In this paper we describe the results of the DASCH commissioning run, during which we developed the data-reduction pipeline and fine-tuned the digitzer's performance and operation. This initial run consisted of 500 plates taken from a variety of different plate-series, all containing the open cluster Praeseppe (M44). We report that accurate photometry at the 0.1mag level is possible on the majority of plates, and demonstrate century-long light-curves of various types of variable stars in and around M44.

We demonstrate that the feedback from stellar bulges can play an essential role in shaping the halo gas of galaxies with substantial bulge components by conducting 1-D hydrodynamical simulations. The feedback model we consider consists of two distinct phases: 1) an early starburst during the bulge formation and 2) a subsequent long-lasting mass and energy injection from stellar winds of low-mass stars and Type Ia SNe. An energetic outward blastwave is initiated by the starburst and is maintained and enhanced by the long-lasting stellar feedback. For a MW-like galactic bulge, this blastwave sweeps up the halo gas in the proto-galaxy and heats up the surrounding medium to a scale much beyond the virial radius of the halo, thus the accretion of the halo hot gas can be completely stopped. In addition, the long-lasting feedback in the later phase powers a galactic bulge wind that is reverse-shocked at a large radius in the presence of surrounding intergalactic medium and hence maintains a hot gaseous halo. As the mass and energy injection decreases with time, the feedback evolves to a subsonic and quasi-stable outflow, which is enough to prevent halo gas from cooling. The two phases of the feedback thus re-enforce each-other's impact on the gas dynamics. The simulation results demonstrate that the stellar bulge feedback may provide a plausible solution to the long-standing problems in understanding the MW type galaxies, such as the "missing stellar feedback" problem and the "over-cooling" problem. The simulations also show that the properties of the hot gas in the subsonic outflow state depend sensitively on the environment and the formation history of the bulge. This dependence and variance may explain the large dispersion in the X-ray to B-band luminosity ratio of the low $L_X/L_B$ Es.

In the unusual intrinsic QSO redshift models, QSOs are ejected by active galaxies with periodic non-cosmological reshifts, thus QSOs are generally associated with active galaxies, and certain structures will be revealed in the QSO redshift distribution. As the largest homogeneous sample of QSOs and galaxies, SDSS data provide the best opportunity to examine this issue. We review the debates on this issue, focused on those based on SDSS and 2dF data, and conclude that there is no strong connection between foreground active galaxies and high-redshift QSOs. The existence of two dips in the SDSS QSO redshift distribution at z=2.7 and 3.5 has recently re-ignited those controversial debates on the origin of QSO redshift. It also turned out that both dips are totally caused by selection effects and after selection effects have been corrected, the two dips disappear and no structure in the redshift distribution of SDSS DR5 sample. These results support that the reshifts of QSOs are cosmological.

This paper has been withdrawn temporarily by the authors, because we are waiting for referee report of the paper submitted to ApJ.

We critically examine the dependence of spectral index on luminosity in optically selected AGN samples. An analysis of optically selected high-redshift quasars showed an anti-correlation of $\alpha_{OX}$, the spectral index between the rest-frame 2500 A and 2 keV, with optical luminosity (Miyaji et al. 2006). We examine this relationship by means of Monte Carlo simulations and conclude that a constant spectral index independent of optical luminosity is still consistent with this high-z sample. We further find that that contributions of large dispersions and narrow range of optical luminosity are most important for the apparent, yet artificial, $\alpha_{OX} - l_o$ correlation reported. We also examine another, but more complete low-z optical selected AGN sub-sample from Steffen et al. (2006), and our analysis shows that a constant spectral index independent of optical luminosity is also consistent with the data. By comparing X-ray and optical luminosity functions, we find that a luminosity independent spectral index is in fact more preferred than the luminosity dependent spectral index model. We also discuss the selection effects caused by flux limits, which might systematically bias the $l_x - l_o$ relation and cause discrepancy in optically selected and X-ray selected AGN samples. To correctly establish a dependence of spectral index of AGNs on their luminosity, a larger and more complete sample is needed and consequences of luminosity dispersions and selection effects in flux limited samples must be taken into account properly.

The prompt optical emission contemporaneous with the $\gamma$-rays from $\gamma$-ray bursts (GRBs) carries important information on the central engine and explosion mechanism. We study the time lag between prompt optical emission and $\gamma$-rays in GRB 990123 and GRB 041219a, which are the only two GRBs had been detected at optical wavelengths during the ascending burst phase. Assuming profiles of prompt optical light curves are the same as the prompt $\gamma$-rays, we simulate optical light curves with different time lags and compare them with the observed optical flux. Then the best fit time lag and its error are determined by chi-squared values. We find that time lags between prompt optical emission and $\gamma$-rays in GRB host galaxy rest-frames are consistent in the two GRBs, which is $5\sim7$ s for GRB 990123 and $1\sim5$ s for GRB 041219a. This result is consistent with a common origin of prompt optical and $\gamma$-ray emissions in the two GRBs. Based on synchrotron cooling model, we also derive the parameters for the two GRBs.

We have used the publicly available data from the Sloan Digital Sky Survey and 2dF QSO Redshift Survey to test the hypothesis that QSOs are ejected from active galaxies with periodic non-cosmological redshifts. For two different intrinsic redshift models, namely the Karlsson $\log(1+z)$ model and Bell's decreasing intrinsic redshift (DIR) model, we do two tests respectively. First, using different criteria, we generate four sets of QSO-galaxy pairs and find there is no evidence for a periodicity at the predicted frequency in $\log(1+z)$, or at any other frequency. We then check the relationship between high redshift QSOs and nearby active galaxies, and we find the distribution of projected distance between high redshift QSOs and nearby active galaxies and the distribution of redshifts of those active galaxies are consistent with a distribution of simulated random pairs, completely different from Bell's previous conclusion. We also analyze the periodicity in redshifts of QSOs, and no periodicity is found in high completeness samples, contrary to the DIR model. These results support that QSOs are not ejected from active galaxies.

Most supernovae are expected to explode in low-density hot media, particularly in galactic bulges and elliptical galaxies. The remnants of such supernovae, though difficult to detect individually, can be profoundly important in heating the media on large scales. We characterize the evolution of this kind of supernova remnants, based on analytical approximations and hydrodynamic simulations. We generalize the standard Sedov solution to account for both temperature and density effects of the ambient media. Although cooling can be neglected, the expansion of such a remnant deviates quickly from the standard Sedov solution and asymptotically approaches the ambient sound speed as the swept-up thermal energy becomes important. The relatively steady and fast expansion of the remnants over large volumes provides an ideal mechanism for spatially distributed heating, which may help to alleviate the over-cooling problem of hot gas in groups and clusters of galaxies as well as in galaxies themselves. The simulations were performed with the FLASH code.

The amplification caused by gravitational light bending by compact objects in a foreground galaxy can affect the apparent number density of background QSOs, as well as their distribution in the fields of galactic halos. In this work we investigate the number enhancement of QSOs in the fields of galactic halos caused by point mass lensing effect and singular isothermal lensing effect, and apply the microlensing effect due to dark compact objects in the halo to NGC 3628. NGC 3628 is a well-studied nearby edge-on Sbc peculiar galaxy, where QSOs are shown to be concentrated around the galaxy with a density much higher than background. We show that if present understanding of the luminosity function of QSOs is right, such concentration could not be caused by gravitational lensing.