We report the discovery of a mid-infrared (MIR) flare using WISE data in the center of the nearby Seyfert 1.9 galaxy MCG-02-04-026. The MIR flare began in the first half of 2014, peaked around the end of 2015, and faded in 2017. During these years, energy more than $7\times10^{50}$ erg was released in the infrared, and the flare's MIR color was generally turning red. We detected neither optical nor ultraviolet (UV) variation corresponding to the MIR flare based on available data. We explained the MIR flare using a dust echo model in which the radiative transfer is involved. The MIR flare can be well explained as thermal reradiation from dust heated by UV-optical photons of a primary nuclear transient event. Although the transient event was not seen directly due to dust obscuration, we can infer that it may produce a total energy of at least $\sim10^{51}$ erg, most of which was released in less than $\sim$3 years. The nature of the transient event could be a stellar tidal disruption event by the central supermassive black hole (SMBH), or a sudden enhancement of the existing accretion flow onto the SMBH, or a supernova which was particularly bright.

The buildup of small-scale magnetic helicity which accompanies the oppositely signed growth on large scales is central to conventional dynamical quenching theories of mean-field dynamos. However, the conventional formalism presumes isotropic turbulence and thereby excludes part of the magnetic Lorentz back-reaction. This renders it insufficient to predict the full quenching for general anisotropic flows. To overcome this deficiency, we derive a new generalized quenching formalism that includes the full back-reacting Lorentz force, and a new "selective-damping-$\tau$" closure which conserves magnetic helicity. We apply the formalism to examples of $\bm\alpha^2$ dynamos and show its predicted quenching for different cases of turbulence---isotropic helical, anisotropic helical, and anisotropic non-helical. It predicts stronger-than-conventional quenching in general, but reduces to the conventional case in the helical isotropic limit.

We present the detection of VHE gamma-ray emission above 100 TeV from HAWC J2227+610 with the HAWC observatory. Combining our observations with previously published results by VERITAS, we interpret the gamma-ray emission from HAWC J2227+610 as emission from protons with a lower limit in their cutoff energy of 800 TeV. The most likely source of the protons is the associated supernova remnant G106.3+2.7, making it a good candidate for a Galactic PeVatron. However, a purely leptonic origin of the observed emission cannot be excluded at this time.

The co-evolution between black holes and galaxies suggests that feedback of active galactic nuclei influence host galaxies through ejecting radiative and kinetic energies to surroundings. Larger scale outflow in local universe are frequently observed by spatially resolved spectroscopy, while smaller scale outflow cannot be directly resolved by current observations. At the scale of the dusty torus, radiative and kinetic energies ejected from the central active nucleus interact with the materials. However, observations of such outflow are rarely reported due to the lack detection of unambiguously gas emission. Here we report the detection of clear and rich emission lines origin from the scale of dusty tours in an partially obscured quasar. The lines share a common intermediate width with full width at half maximum about 1900 \kmps\ and are shown in two systems: a major system is unshifted and a minor system has a blue-shifts of about 2600 \kmps. The line intensity ratios, combining photo-ionization simulations, indicates an ultradense line-emitting region with the density as high as $\sim$ $10^{13}~\rm cm^{-3}$. We interpret this as the lines being excited by a shock induced by the high-density and high-temperature gases at the scale of dusty torus, rather than photo-ionized by the central accretion disk. We speculate that the outflow, launched from the accretion disk, collides onto the inner wall of the dusty torus and shock-heat the gases to cause the major emission lines. The outflowing gases may also collide onto surrounding isolated clouds, and give rise to blue-shifted minor emission lines.

Due to the high energies and long distances to the sources, astrophysical observations provide a unique opportunity to test possible signatures of Lorentz invariance violation (LIV). Superluminal LIV enables the decay of photons at high energy. The High Altitude Water Cherenkov (HAWC) Observatory is among the most sensitive gamma-ray instruments currently operating above 10 TeV. HAWC finds evidence of 100 TeV photon emission from at least four astrophysical sources. These observations exclude, for the strongest of the limits set, the LIV energy scale to $2.2\times10^{31}$ eV, over 1800 times the Planck energy and an improvement of 1 to 2 orders of magnitude over previous limits.

We present multiwavelength studies of a transient X-ray source, XMMSL1 J131952.3+225958, associated with the galaxy NGC 5092 at $z=0.023$ detected in the XMM-Newton SLew survey (XMMSL). The source brightened in the 0.2--2 keV band by a factor of $>20$ in 2005 as compared with previous flux limits and then faded by a factor of $>200$ as observed with it XMM-Newton in 2013 and with it Swift in 2018. At the flaring state, the X-ray spectrum can be modeled with a blackbody at a temperature of $\sim$ 60 eV and an overall luminosity of $\sim$ $1.5 \times 10^{43}$ erg s$^{-1}$. A UV flare and optical flare were also detected with the Galaxy Evolution Explorer and the Sloan Digital Sky Survey, respectively, within several months of the X-ray flare, whose nonstellar UV--optical spectrum can be described with a blackbody at a temperature of $\sim$ $(1-2) \times 10^4$ K and a luminosity of $\sim$ $(2-6) \times 10^{43}$ erg s$^{-1}$. Interestingly, mid-infrared monitoring observations of NGC 5092 with the Wide-field Infrared Survey Explorer 5--13 yr later show a continuous flux decline. These dramatic variability properties, from the X-ray through UV and optical to infrared, appear to be orderly, suggestive of a stellar tidal disruption event (TDE) by a massive black hole, confirming the postulation by Kanner et al.(2013). This TDE candidate belongs to a rare sample with contemporaneous bright emission detected in the X-ray, UV, and optical, which are later echoed by dust-reprocessed light in the mid-infrared. The black hole has a mass of $\sim$ $5 \times 10^{7} \rm M_{\odot}$, residing in a galaxy that is dominated by a middle-aged stellar population of 2.5 Gyr.

We examined the solar gravitational lens (SGL) as the means to produce direct high-resolution, multipixel images of exoplanets. The properties of the SGL are remarkable: it offers maximum light amplification of ~1e11 and angular resolution of ~1e-10 arcsec. A probe with a 1-m telescope in the SGL focal region can image an exoplanet at 30 pc with 10-kilometer resolution on its surface, sufficient to observe seasonal changes, oceans, continents, surface topography. We reached and exceeded all objectives set for our study: We developed a new wave-optical approach to study the imaging of exoplanets while treating them as extended, resolved, faint sources at large but finite distances. We properly accounted for the solar corona brightness. We developed deconvolution algorithms and demonstrated the feasibility of high-quality image reconstruction under realistic conditions. We have proven that multipixel imaging and spectroscopy of exoplanets with the SGL are feasible. We have developed a new mission concept that delivers an array of optical telescopes to the SGL focal region relying on three innovations: i) a new way to enable direct exoplanet imaging, ii) use of smallsats solar sails fast transit through the solar system and beyond, iii) an open architecture to take advantage of swarm technology. This approach enables entirely new missions, providing a great leap in capabilities for NASA and the greater aerospace community. Our results are encouraging as they lead to a realistic design for a mission that will be able to make direct resolved images of exoplanets in our stellar neighborhood. It could allow exploration of exoplanets relying on the SGL capabilities decades, if not centuries, earlier than possible with other extant technologies. The architecture and mission concepts for a mission to the strong interference region of the SGL are promising and should be explored further.

Cosmic rays, along with stellar radiation and magnetic fields, are known to make up a significant fraction of the energy density of galaxies such as the Milky Way. When cosmic rays interact in the interstellar medium, they produce gamma-ray emission which provides an important indication of how the cosmic rays propagate. Gamma rays from the Andromeda Galaxy (M31), located 785 kpc away, provide a unique opportunity to study cosmic-ray acceleration and diffusion in a galaxy with a structure and evolution very similar to the Milky Way. Using 33 months of data from the High Altitude Water Cherenkov Observatory, we search for TeV gamma rays from the galactic plane of M31. We also investigate past and present evidence of galactic activity in M31 by searching for Fermi Bubble-like structures above and below the galactic nucleus. No significant gamma-ray emission is observed, so we use the null result to compute upper limits on the energy density of cosmic rays $>10$ TeV in M31. The computed upper limits are approximately ten times higher than expected from the extrapolation of the Fermi LAT results.

Primordial Black Holes (PBHs) may have been created by density fluctuations in the early Universe and could be as massive as $> 10^9$ solar masses or as small as the Planck mass. It has been postulated that a black hole has a temperature inversely-proportional to its mass and will thermally emit all species of fundamental particles via Hawking Radiation. PBHs with initial masses of $\sim 5 \times 10^{14}$ g (approximately one gigaton) should be expiring today with bursts of high-energy gamma radiation in the GeV--TeV energy range. The High Altitude Water Cherenkov (HAWC) Observatory is sensitive to gamma rays with energies of $\sim$300 GeV to past 100 TeV, which corresponds to the high end of the PBH burst spectrum. With its large instantaneous field-of-view of $\sim 2$ sr and a duty cycle over 95%, the HAWC Observatory is well suited to perform an all-sky search for PBH bursts. We conducted a search using 959 days of HAWC data and exclude the local PBH burst rate density above $3400~\mathrm{pc^{-3}~yr^{-1}}$ at 99% confidence, the strongest limit on the local PBH burst rate density from any existing electromagnetic measurement.

We present the first catalog of gamma-ray sources emitting above 56 and 100 TeV with data from the High Altitude Water Cherenkov (HAWC) Observatory, a wide field-of-view observatory capable of detecting gamma rays up to a few hundred TeV. Nine sources are observed above 56 TeV, all of which are likely Galactic in origin. Three sources continue emitting past 100 TeV, making this the highest-energy gamma-ray source catalog to date. We report the integral flux of each of these objects. We also report spectra for three highest-energy sources and discuss the possibility that they are PeVatrons.

Detection of non-gravitational interactions of massive dark matter with visible sector so far have given null results. The DM may communicate with the ordinary matter only through gravitational interactions. Besides, the majority of traditional direct detections have poor sensitivity of searching for light DM particle because of the small recoil energy. Thanks to the high energy cosmic rays (CRs), the light DM can be boosted by scattering with CRs and thus may be detected in those existing experiments. In this work, we first derive the exclusion limits of the cosmic ray boosted sub-GeV Dark Matter with gravitational mediator from the Xenon1T data.

Model-independent measurements for the cosmic spatial curvature, which is related to the nature of cosmic space-time geometry, plays an important role in cosmology. On the basis of the Distance Sum Rule in the Friedmann-Lema{\^i}tre-Robertson-Walker metric, (distance ratio) measurements of strong gravitational lensing (SGL) systems together with distances from type Ia supernovae observations have been proposed to directly estimate the spatial curvature without any assumptions for the theories of gravity and contents of the universe. However, previous studies indicated that a spatially closed universe was strongly preferred. In this paper, we re-estimate the cosmic curvature with the latest SGL data which includes 163 well-measured systems. In addition, possible factors, e.g. combination of SGL data from different surveys and stellar mass of the lens galaxy, which might affect estimations for the spatial curvature, are considered in our analysis. We find that, except the case where only SGL systems from the Sloan Lens ACS Survey are considered, a spatially flat universe is consistently favored at very high confidence level by the latest observations. It is suggested that the increasing number of well-measured strong lensing events might significantly reduce the bias of estimation for the cosmic curvature.

IRAS F11119+3257 is a quasar-dominated Ultra-Luminous InfraRed Galaxy, with a partially obscured narrow-line seyfert 1 nucleus. In this paper, we present the NIR spectroscopy of F11119+3257, in which we find unusual Paschen emission lines, and metastable He I* $\lambda$10830 absorption associated with the previously reported atomic sodium and molecular OH mini-BAL (Broad Absorption Line) outflow. Photo-ionization diagnosis confirms previous findings that the outflows are at kilo-parsec scales. Such large-scale outflows should produce emission lines. We indeed find that high-ionization emission lines ([O III], [Ne III], and [Ne V]) are dominated by blueshifted components at similar speeds to the mini-BALs. The blueshifted components are also detected in some low-ionization emission lines, such as [O II] $\lambda$3727 and some Balmer lines (H$\alpha$, H$\beta$, and H$\gamma$), even though their cores are dominated by narrow ($FWHM_{\rm NEL} = 570\pm40$km s$^{-1}$) or broad components at the systemic redshift of $z=0.18966\pm0.00006$. The mass flow rate (230-730$~M_\odot \rm yr^{-1}$) and the kinetic luminosity ($\dot{E}_k \sim 10^{43.6-44.8} $erg s$^{-1}$) are then inferred jointly from the blueshifted emission and absorption lines. In the NIR spectrum of F11119+3257, we also find that the Paschen emission lines are unique, in which a very narrow ($FWHM=260\pm20~$km s$^{-1}$) component is shown in only Pa$\alpha$. This narrow component most probably comes from heavily obscured star formation. Based on the Pa$\alpha$ and Pa$\beta$ emissions, we obtain an extinction at the $H$ band, $A_H~>~2.1$ (or a reddenning of $E_{B-V}~>~$3.7), and a star formation rate of $SFR~>~130\rm M_\odot yr^{-1}$ that resembles the estimates inferred from the FIR emissions ($SFR_{\rm FIR} = 190\pm90$ $M_\odot$ yr$^{-1}$).

We present spectropolarimetric and spectrophotometric observations of the peculiar active galactic nucleus (AGN) SDSS J120300.19+162443.7 (hereafter J1203+1624) at $z=0.1656$. Its optical total flux spectra clearly show broad emission lines (BELs) in H$\alpha$ and H$\beta$. After removal of narrow emission lines (NELs), the full width at half maximum of the lorentzian BEL is $FWHM_{\rm BEL} \approx 1,270$ km s$^{-1}$, fulfilling the conventional definition of a narrow line Seyfert 1 (NLS1) galaxy. However, its NELs are unprecedentedly strong when compared to type 1AGNs. This, together with its large MIR excess ($g-W_4 = 13.172$ mag), implies that the observer and the narrow emission line region (NELR) might see a different ionization continuum. Our optical spectropolarimetry confirms its type 2 nature by detecting a polarized blue continuum and Balmer BELs ($FWHM_{\rm Polarized~BEL} \approx 1,183$ km s$^{-1}$), with a high polarization degree of $> 20$% in the blue wing. The mass and Eddington rate of the central black hole is estimated based on both transmitted and scattered AGN radiation, which is $M_{\bullet} < 2.9 \times 10^7 M_\odot$ and $L_{\rm bol}/L_{\rm Edd} > 1.5$. Severe extinction of the AGN emission also enables a clear view of the compact host galaxy. Discovery of J1203+1624 suggests that NLS1 follows AGN unfication scheme, studying of its analogs could blaze a new trail for exploring the connection between black hole growth and star formation in the host galaxy. Interesting features of J1203+1624 like the peculiar NELs, inflowing scattering clouds within sublimation radius are worth detailed follow-ups in the future.

We report the detection of a large amplitude MIR outburst in the white dwarf (WD) 0145+234 in the NEOWISE Survey data. The source had a stable MIR flux before 2018, and was brightened by about 1.0 magnitude in the W1 and W2 bands within half a year and has been continuously brightening since then. No significant variations are found in the optical photometry data during the same period. This suggests that this MIR outburst is caused by recent replenishing or redistribution of dust, rather than intrinsic variations of the WD. SED modeling of 0145+234 suggests that there was already a dust disk around the WD in the quiescent state, and both of the temperature and surface area of the disk evolved rapidly since the outburst. The dust temperature was about 1770K in the initial rising phase, close to the sublimation temperature of silicate grains, and gradually cooled down to around 1150K, while the surface area increased by a factor of about 6 during the same period. The inferred closest distance of dust to the WD is within the tidal disruption radius of a gravitationally bounded asteroid. We estimated the dust mass to be between $3\times10^{15}$ to $3\times10^{17}$ $\rho/(1 \mathrm g~cm^{-3})$ kg for silicate grains of a power-law size distribution with a high cutoff size from 0.1 to 1000 microns. We interpret this as a possible tidal breakup of an exo-asteroid by the WD. Further follow-up observations of this rare event may provide insights on the origin of dust disk and metal pollution in some white dwarfs.

We extract the multiple-epoch Balmer-line profiles of the heavily obscured quasar 3C 445 from the spectral curves in the literature, and analyze the emission-line profiles of the H$\alpha$ and H$\beta$ lines and the profile variability in the H$\alpha$ line in the large time interval of more than three decades. The profile comparison between the H$\alpha$ and H$\beta$ lines shows that both Balmer lines share the profile with the same form, while the blue system of the H$\beta$ line is seriously weaker than that of the H$\alpha$ line. Moreover, the blue system of the H$\alpha$ line suddenly disappeared completely and then did not appear again, however the other two components did not exhibit significant variation in the velocity or the amplitude. These findings suggest that the blue system of 3C 445, as with SDSS J153636.22+044127.0 and its analogs, is probable the result of the shock-heated outflowing gases. The observation angle of almost edge-on which the previous studies suggested can easily produce the high-speed and high-temperature shock in the collision between the massive outflow and the inner surface of the dusty torus.

We present TeV gamma-ray observations of the Crab Nebula, the standard reference source in ground-based gamma-ray astronomy, using data from the High Altitude Water Cherenkov (HAWC) Gamma-Ray Observatory. In this analysis we use two independent energy-estimation methods that utilize extensive air shower variables such as the core position, shower angle, and shower lateral energy distribution. In contrast, the previously published HAWC energy spectrum roughly estimated the shower energy with only the number of photomultipliers triggered. This new methodology yields a much improved energy resolution over the previous analysis and extends HAWC's ability to accurately measure gamma-ray energies well beyond 100 TeV. The energy spectrum of the Crab Nebula is well fit to a log parabola shape $\left(\frac{dN}{dE} = \phi_0 \left(E/\textrm{7 TeV}\right)^{-\alpha-\beta\ln\left(E/\textrm{7 TeV}\right)}\right)$ with emission up to at least 100 TeV. For the first estimator, a ground parameter that utilizes fits to the lateral distribution function to measure the charge density 40 meters from the shower axis, the best-fit values are $\phi_o$=(2.35$\pm$0.04$^{+0.20}_{-0.21}$)$\times$10$^{-13}$ (TeV cm$^2$ s)$^{-1}$, $\alpha$=2.79$\pm$0.02$^{+0.01}_{-0.03}$, and $\beta$=0.10$\pm$0.01$^{+0.01}_{-0.03}$. For the second estimator, a neural network which uses the charge distribution in annuli around the core and other variables, these values are $\phi_o$=(2.31$\pm$0.02$^{+0.32}_{-0.17}$)$\times$10$^{-13}$ (TeV cm$^2$ s)$^{-1}$, $\alpha$=2.73$\pm$0.02$^{+0.03}_{-0.02}$, and $\beta$=0.06$\pm$0.01$\pm$0.02. The first set of uncertainties are statistical; the second set are systematic. Both methods yield compatible results. These measurements are the highest-energy observation of a gamma-ray source to date.

Quasars are high-luminosity active galactic nuclei believed to be powered by accretion of interstellar matter onto a super-massive black hole (SMBH) therein. Most of the observed energy is released in an accretion disk of inspiralling gas surrounding the SMBH. An enormous amount of fueling material is expected to be transported inwards. However, basic questions remain unanswered as to whether and how the accretion disks are supplied with external gas, since no disk-feeding inflow has hitherto been observed clearly. Here we report the discovery of highly redshifted broad absorption lines arising from neutral hydrogen and helium atoms in a small sample of quasars. Their absorption troughs show a broad range of Doppler velocities from zero extending continuously inward up to as high as $\sim 5,000$ km s$^{-1}$. We thus see through streams of cold gas moving with a radially inward velocity component that spans an immense gradient---a result of gravitational acceleration by the central SMBH. Extensive photo-ionization modeling for the archetypical object SDSS J103516.20+142200.6 indicates the inflowing gas to be dense, thick and moderately ionized, with a characteristic distance to the SMBH of $\sim 1,000$ gravitational radii, possibly overlapping or close to the outer accretion disk. Our results present the first compelling evidence for the long-sought inflow directly feeding quasars' accretion disks with external materials, likely originating from the dusty torus at a parsec scale. Our approach provides a new tool to probe the bulk of the so far elusive fueling inflows in quasars. Their studies may help address some of the fundamental questions concerning accretion physics, the onset and sustainment of quasar activity, and the SMBH growth at centers of most galaxies.

List of proceedings from the HAWC Collaboration presented at the 36th International Cosmic Ray Conference, 24 July - 1 August 2019, Madison, Wisconsin, USA.

The High Altitude Water Cherenkov (HAWC) observatory recently published the discovery of SS 433 as a TeV source, reporting the observation of multi-TeV gamma-ray emission from the jet interaction regions e1 and w1, suggesting in-situ particle acceleration. This showed the first direct evidence of acceleration in jets at energies greater than a few TeV. SS 433 was the first microquasar to be discovered and is still considered special in that the accretion is supercritical and the luminosity of the system is very high ($\sim10^{40}$ erg s$^{-2}$). The lobes of the supernova remnant W 50 in which the jets terminate, about 40 parsecs from the central binary, are expected to accelerate charged particles, and indeed radio and X-ray emission consistent with electron synchrotron emission in a magnetic field have been observed. SS 433 has also been a strong candidate for hadronic acceleration due to spectroscopic evidence of ionized nuclei in the inner jets. However, multiwavelength fits including the HAWC measurements favor the leptonic production of the observed gamma rays. Here, we present new follow-up measurements of the jet interaction regions of SS 433 using the most recent data from HAWC.

We describe plans for the development of the Southern Wide-field Gamma-ray Observatory (SWGO), a next-generation instrument with sensitivity to the very-high-energy (VHE) band to be constructed in the Southern Hemisphere. SWGO will provide wide-field coverage of a large portion of the southern sky, effectively complementing current and future instruments in the global multi-messenger effort to understand extreme astrophysical phenomena throughout the universe. A detailed description of science topics addressed by SWGO is available in the science case white paper [1]. The development of SWGO will draw on extensive experience within the community in designing, constructing, and successfully operating wide-field instruments using observations of extensive air showers. The detector will consist of a compact inner array of particle detection units surrounded by a sparser outer array. A key advantage of the design of SWGO is that it can be constructed using current, already proven technology. We estimate a construction cost of 54M USD and a cost of 7.5M USD for 5 years of operation, with an anticipated US contribution of 20M USD ensuring that the US will be a driving force for the SWGO effort. The recently formed SWGO collaboration will conduct site selection and detector optimization studies prior to construction, with full operations foreseen to begin in 2026. Throughout this document, references to science white papers submitted to the Astro2020 Decadal Survey with particular relevance to the key science goals of SWGO, which include unveiling Galactic particle accelerators [2-10], exploring the dynamic universe [11-21], and probing physics beyond the Standard Model [22-25], are highlighted in red boldface.

A new, complete sample of 14,584 broad-line AGNs at $z<0.35$ is presented, which are uncovered homogeneously from the complete database of galaxies and quasars observed spectroscopically in the Sloan Digital Sky Survey Seventh Data Release. The stellar continuum is properly removed for each spectrum with significant host absorption line features, and careful analyses of the emission-line spectra, particularly in the H$\rm \alpha$ and H$\rm \beta$ wavebands, are carried out. The broad Balmer emission line, particularly, H$\rm \alpha$, is used to indicate the presence of an AGN. The broad H$\rm \alpha$ lines have luminosities in a range of $10^{38.5}$-$10^{44.3}$ erg s$^{-1}$, and line widths (FWHMs) of 500-34,000 km s$^{-1}$. The virial black hole masses, estimated from the broad line measurements, span a range of $10^{5.1}$-$10^{10.3}$ $M_\odot$, and the Eddington ratios vary from $-3.3$ to $1.3$ in logarithmic scale. Other quantities such as multi-wavelength photometric properties and flags denoting peculiar line profiles are also included in this catalog. We describe the construction of this catalog and briefly discuss its properties. The catalog is publicly available online. This homogeneously selected AGN catalog, along with the accurately measured spectral parameters, provide the most updated, largest AGN sample data, which will enable further comprehensive investigations of the properties of the AGN population in the low-redshift universe.

Simulations of dark matter show a discrepancy between the expected number of Galactic dark matter sub-halos and how many have been optically observed. Some of these unseen satellites may exist as dark dwarf galaxies: sub-halos like dwarf galaxies with no luminous counterpart. Assuming WIMP dark matter, it may be possible to detect these unseen sub-halos from gamma-ray signals originating from dark matter annihilation. The High Altitude Water Cherenkov Observatory (HAWC) is a very high energy (500 GeV to 100 TeV) gamma ray detector with a wide field-of-view and near continuous duty cycle, making HAWC ideal for unbiased sky surveys. We perform such a search for gamma ray signals from dark dwarfs in the Milky Way halo. We perform a targeted search of HAWC gamma-ray sources which have no known association with lower-energy counterparts, based on an unbiased search of the entire sky. With no sources found to strongly prefer dark matter models, we calculate the ability of HAWC to observe dark dwarfs. We also compute the HAWC sensitivity to potential future detections for a given model of dark matter substructure.

As hybrids of narrow-line Seyfert 1 (NLS1) galaxies and blazars, {\gamma}-ray emitting NLS1s are important probes of jet physics in the high Eddington-ratio regime. Only very few of them are known to date; the majority of them below redshift z = 0.5. Here we present the identification of the {\gamma}-ray emitting AGN TXS 0943+105 (SDSS J094635.06+101706.1) as a high-redshift NLS1 galaxy. It turns out to be one of the radio-loudest NLS1s known, highly variable at all wavelengths, and shows widely extended radio emission at a (projected) > 100 kpc scale. It is a known strong {\gamma}-ray emitter with a luminous flare reported previously. At redshift z=1.004, this is the most distant {\gamma}-NLS1 known to date.

Pulsar Wind Nebula (PWN) DA 495 (G65.7+1.2) was detected in TeV gamma-rays by the High Altitude Water Cherenkov Observatory (HAWC) in 2017 (2HWC J1953+294). Follow-up observations by the Very Energetic Radiation Imaging Telescope Array System (VERITAS) confirmed the association between 2HWC J1953+294 and DA 495 and found the TeV emission to be spatially coincident with the radio emission first reported in 1968. The detection of TeV gamma-rays from DA 495, along with past X-ray detection up to 10 keV, prompted high energy X-ray observations as part of the NuSTAR Galactic Legacy Survey. We present the results of these NuSTAR observations, combined with archival Chandra and XMM-Newton observations, and confirm the previous X-ray photon index of $\Gamma_{2-20 \rm\ keV} = 2.0 \pm 0.1$. We find no spectral cutoff up to 20 keV. With the spectral information for DA 495 extended to TeV gamma-rays, we were able to perform analytical modeling to test leptonic and hadronic emission scenarios. The leptonic models can explain the broadband emission, but also imply a diffuse X-ray nebula of similar extent to the radio and TeV nebulae, which cannot be confirmed by our observations. The hadronic models can simultaneously explain the spectrum and the spatial extent in all wavelengths; however, we need a very high magnetic field strength pervading the radio and TeV nebulae and a surprisingly high particle kinetic energy. These requirements deepen the mystery of the physical nature of DA 495. Future observations in radio to infrared bands and spatially resolved $\gamma$-rays can further constrain the physical conditions and radiation mechanisms in DA 495.

The binary emission-line system, particularly the binary broad-line emission system, is considered the most effective indicator of the black hole binary. A plausible example of such a system, SDSS J153636.22+044127.0, was reported as the first known object with two hydrogen Balmer broad-line systems, which are interpreted to be the result of broad-line regions around a pair of black holes (Boroson \& Lauer 2009). Here, we show the follow-up optical and near-infrared spectral observations of SDSS J153636.22+044127.0 and its analogues. In these objects, the broad hydrogen Balmer and Paschen, He I and Mg II lines share the same peculiar emission-line profile (including a blue system, a red system and a double-peaked disk-line component); however, the invariance in the large time interval, the absence of the blue system in He I $\lambda$10830 profile and the abnormally strong emission of the hydrogen Pa$\beta$ blue system oppose the binary proposal. We suggest that these unique broad lines arise from the AGN emission-line region and the shock-heated outflowing gases rather than a binary system of two active black holes.

Mean field electrodynamics (MFE) facilitates practical modeling of secular, large scale properties of astrophysical or laboratory systems with fluctuations.Practitioners commonly assume wide scale separation between mean and fluctuating quantities, to justify equality of ensemble and spatial or temporal averages.Often however, real systems do not exhibit such scale separation. This raises two questions: (I) what are the appropriate generalized equations of MFE in the presence of mesoscale fluctuations? (II) how precise are theoretical predictions from MFE? We address both by first deriving the equations of MFE for different types of averaging, along with mesoscale correction terms that depend on the ratio of averaging scale to variation scale of the mean. We then show that even if these terms are small, predictions of MFE can still have a significant precision error. This error has an intrinsic contribution from the dynamo input parameters and a filtering contribution from differences in the way observations and theory are projected through the measurement kernel.Minimizing the sum of these contributions can produce an optimal scale of averaging that makes the theory maximally precise.The precision error is important to quantify when comparing to observations because it quantifies the resolution of predictive power. We exemplify these principles for galactic dynamos, comment on broader implications, and identify possibilities for further work.

TeV Halos, extended regions of TeV gamma-ray emission around middle-aged pulsars, have recently been established as a new source class in gamma-ray astronomy. These halos have been attributed to relativistic electrons and positrons that have left the acceleration region close to the pulsar and are diffusing in the surrounding medium. Measuring the morphology of TeV Halos enables for the first time a direct measurement of the electron diffusion on scales of tens of parsecs. There are hints that the presence of relativistic particles affects the diffusion rate in the pulsars' surroundings. Understanding electron diffusion is necessary to constrain the origins of the apparent `excess' of cosmic-ray positrons at tens of GeV. TeV Halos can also be used to find mis-aligned pulsars, as well as study certain properties of the Galaxy's pulsar population. Future VHE gamma-ray instruments will detect more of those TeV Halos and determine how much pulsars contribute to the observed cosmic-ray electron and positron fluxes, and how they affect diffusion in their environments.

The Andromeda Galaxy (M31) is a nearby ($\sim$780 kpc) galaxy similar to our own Milky Way. Observational evidence suggests that it resides in a large halo of dark matter (DM), making it a good target for DM searches. We present a search for gamma rays from M31 using 1017 days of data from the High Altitude Water Cherenkov (HAWC) Observatory. With its wide field of view and constant monitoring, HAWC is well-suited to search for DM in extended targets like M31. No DM annihilation or decay signal was detected for DM masses from 1 to 100 TeV in the $b\bar{b}$, $t\bar{t}$, $\tau^{+}\tau^{-}$, $\mu^{+}\mu^{-}$, and $W^{+}W^{-}$ channels. Therefore we present limits on those processes. Our limits nicely complement the existing body of DM limits from other targets and instruments. Specifically the DM decay limits from our benchmark model are the most constraining for DM masses from 25 TeV to 100 TeV in the $b\bar{b}, t\bar{t}$ and $\mu^{+}\mu{-}$ channels. In addition to DM-specific limits, we also calculate general gamma-ray flux limits for M31 in 5 energy bins from 1 TeV to 100 TeV.

We outline the science motivation for SGSO, the Southern Gamma-Ray Survey Observatory. SGSO will be a next-generation wide field-of-view gamma-ray survey instrument, sensitive to gamma-rays in the energy range from 100 GeV to hundreds of TeV. Its science topics include unveiling galactic and extragalactic particle accelerators, monitoring the transient sky at very high energies, probing particle physics beyond the Standard Model, and the characterization of the cosmic ray flux. SGSO will consist of an air shower detector array, located in South America. Due to its location and large field of view, SGSO will be complementary to other current and planned gamma-ray observatories such as HAWC, LHAASO, and CTA.

Quasar outflows carry mass, momentum and energy into the surrounding environment, and have long been considered a potential key factor in regulating the growth of supermassive black holes and the evolution of their host galaxies. A crucial parameter for understanding the origin of these outflows and measuring their influence on their host galaxies is the distance (R) between the outflow gas and the galaxy center. While R has been measured in a number of individual galaxies, its distribution remains unknown. Here we report the distributions of R and the kinetic luminosities of quasars outflows, using the statistical properties of broad absorption line variability in a sample of 915 quasars from the Sloan Digital Sky Surveys. The mean and standard deviation of the distribution of R are 10^{1.4+/-0.5} parsecs. The typical outflow distance in this sample is tens of parsec, which is beyond the theoretically predicted location (0.01 ~ 0.1 parsecs) where the accretion disc line-driven wind is launched, but is smaller than the scales of most outflows that are derived using the excited state absorption lines. The typical value of the mass-flow rate is of tens to a hundred solar masses per year, or several times the accretion rate. The typical kinetic-to-bolometric luminosity ratio is a few per cent, indicating that outflows are energetic enough to influence the evolution of their host galaxies.

We present the first full-sky analysis of the cosmic ray arrival direction distribution with data collected by the HAWC and IceCube observatories in the Northern and Southern hemispheres at the same median primary particle energy of 10 TeV. The combined sky map and angular power spectrum largely eliminate biases that result from partial sky coverage and holds a key to probe into the propagation properties of TeV cosmic rays through our local interstellar medium and the interaction between the interstellar and heliospheric magnetic fields. From the map we determine the horizontal dipole components of the anisotropy $\delta_{0h} = 9.16 \times 10^{-4}$ and $\delta_{6h} = 7.25 \times 10^{-4}~(\pm0.04 \times 10^{-4})$. In addition, we infer the direction ($229.2\pm 3.5^\circ$ RA , $11.4\pm 3.0^\circ$ Dec.) of the interstellar magnetic field from the boundary between large scale excess and deficit regions from which we estimate the missing corresponding vertical dipole component of the large scale anisotropy to be $\delta_N \sim -3.97 ^{+1.0}_{-2.0} \times 10^{-4}$.

The HAWC Collaboration released the 2HWC catalog of TeV sources, in which 19 show no association with any known high-energy (HE; E > 10 GeV) or very-high-energy (VHE; E > 300 GeV) sources. This catalog motivated follow-up studies by both the MAGIC and Fermi-LAT observatories with the aim of investigating gamma-ray emission over a broad energy band. In this paper, we report the results from the first joint work between HAWC, MAGIC and Fermi-LAT on three unassociated HAWC sources: 2HWC J2006+341, 2HWC J1907+084* and 2HWC J1852+013*. Although no significant detection was found in the HE and VHE regimes, this investigation shows that a minimum 1 degree extension (at 95% confidence level) and harder spectrum in the GeV than the one extrapolated from HAWC results are required in the case of 2HWC J1852+013*, while a simply minimum extension of 0.16 degrees (at 95% confidence level) can already explain the scenario proposed by HAWC for the remaining sources. Moreover, the hypothesis that these sources are pulsar wind nebulae is also investigated in detail.

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

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

Transport coefficients in turbulence are comprised of correlation functions between turbulent fluctuations and efficient methods to calculate them are desirable. For example, in mean field dynamo theories used to model the growth of large scale magnetic fields of stars and galaxies, the turbulent electromotive force is commonly approximated by a series of tensor products of turbulent transport coefficients with successively higher order spatial derivatives of the mean magnetic field. One ingredient of standard models is the kinematic coefficient of the zeroth order term, namely the averaged kinetic pseudotensor $\bm{\alpha}$, that converts toroidal to poloidal fields. Here we demonstrate an efficient way to calculate this quantity for rotating stratified turbulence, whereby the pre-averaged quantity is calculated for the motion of a single plume, and the average is then taken over an ensemble of plumes of different orientations. We calculate the plume dynamics in the most convenient frame, before transforming back to the lab frame and averaging. Our concise configuration space calculation gives essentially identical results to previous lengthier approaches. The present application exemplifies what is a broadly applicable method.

We analyze the Sun as a source for the indirect detection of dark matter through a search for gamma rays from the solar disk. Capture of dark matter by elastic interactions with the solar nuclei followed by annihilation to long-lived mediators can produce a detectable gamma-ray flux. We search three years of data from the High Altitude Water Cherenkov Observatory and find no statistically significant detection of TeV gamma-ray emission from the Sun. Using this, we constrain the spin-dependent elastic scattering cross section of dark matter with protons for dark matter masses above 1 TeV, assuming an unstable mediator with a favorable lifetime. The results complement constraints obtained from Fermi-LAT observations of the Sun and together cover WIMP masses between 4 GeV and $10^6$ GeV. The cross section constraints for mediator decays to gamma rays can be as strong as $\sim10^{-45}$ cm$^{-2}$, which is more than four orders of magnitude stronger than current direct-detection experiments for 1 TeV dark matter mass. The cross-section constraints at higher masses are even better, nearly 7 orders of magnitude better than the current direct-detection constraints for 100 TeV dark matter mass. This demonstration of sensitivity encourages detailed development of theoretical models in light of these powerful new constraints.

Steady gamma-ray emission up to at least 200 GeV has been detected from the solar disk in the Fermi-LAT data, with the brightest, hardest emission occurring during solar minimum. The likely cause is hadronic cosmic rays undergoing collisions in the Sun's atmosphere after being redirected from ingoing to outgoing in magnetic fields, though the exact mechanism is not understood. An important new test of the gamma-ray production mechanism will follow from observations at higher energies. Only the High Altitude Water Cherenkov (HAWC) Observatory has the required sensitivity to effectively probe the Sun in the TeV range. Using three years of HAWC data from November 2014 to December 2017, just prior to the solar minimum, we search for 1--100 TeV gamma rays from the solar disk. No evidence of a signal is observed, and we set strong upper limits on the flux at a few $10^{-12}$ TeV$^{-1}$ cm$^{-2}$ s$^{-1}$ at 1 TeV. Our limit, which is the most constraining result on TeV gamma rays from the Sun, is $\sim10\%$ of the theoretical maximum flux (based on a model where all incoming cosmic rays produce outgoing photons), which in turn is comparable to the Fermi-LAT data near 100 GeV. The prospects for a first TeV detection of the Sun by HAWC are especially high during solar minimum, which began in early 2018.

We report the discovery of an active intermediate-mass black hole (IMBH) candidate in the center of nearby barred bulgeless galaxy $\rm NGC~3319$. The point X-ray source revealed by archival Chandra and XMM-Newton observations is spatially coincident with the optical and UV galactic nuclei from Hubble Space Telescope observations. The spectral energy distribution derived from the unresolved X-ray and UV-optical flux is comparable with active galactic nuclei (AGNs) rather than ultra-luminous X-ray sources, although its bolometric luminosity is only $3.6\times10^{40}~\rm erg~s^{-1}$. Assuming an Eddington ratio range between 0.001 and 1, the black hole mass ($M_\rm{BH}$) will be located at $3\times10^2-3\times10^5~M_{\odot}$, placing it in the so-called IMBH regime and could be the one of the lowest reported so far. Estimates from other approaches (e.g., fundamental plane, X-ray variability) also suggest $M_\rm{BH}\lesssim10^5~M_{\odot}$.Similar to other BHs in bulgeless galaxies, the discovered IMBH resides in a nuclear star cluster with mass of $\sim6\times10^6~M_{\odot}$. The detection of such a low-mass BH offers us an ideal chance to study the formation and early growth of SMBH seeds, which may result from the bar-driven inflow in late-type galaxies with a prominent bar such as $\rm NGC~3319$.

SS 433 is a binary system containing a supergiant star that is overflowing its Roche lobe with matter accreting onto a compact object (either a black hole or neutron star). Two jets of ionized matter with a bulk velocity of $\sim0.26c$ extend from the binary, perpendicular to the line of sight, and terminate inside W50, a supernova remnant that is being distorted by the jets. SS 433 differs from other microquasars in that the accretion is believed to be super-Eddington, and the luminosity of the system is $\sim10^{40}$ erg s$^{-1}$. The lobes of W50 in which the jets terminate, about 40 pc from the central source, are expected to accelerate charged particles, and indeed radio and X-ray emission consistent with electron synchrotron emission in a magnetic field have been observed. At higher energies (>100 GeV), the particle fluxes of $\gamma$ rays from X-ray hotspots around SS 433 have been reported as flux upper limits. In this energy regime, it has been unclear whether the emission is dominated by electrons that are interacting with photons from the cosmic microwave background through inverse-Compton scattering or by protons interacting with the ambient gas. Here we report TeV $\gamma$-ray observations of the SS 433/W50 system where the lobes are spatially resolved. The TeV emission is localized to structures in the lobes, far from the center of the system where the jets are formed. We have measured photon energies of at least 25 TeV, and these are certainly not Doppler boosted, because of the viewing geometry. We conclude that the emission from radio to TeV energies is consistent with a single population of electrons with energies extending to at least hundreds of TeV in a magnetic field of $\sim16$~micro-Gauss.

The HAWC (High Altitude Water Cherenkov) collaboration recently published their 2HWC catalog, listing 39 very high energy (VHE; >100~GeV) gamma-ray sources based on 507 days of observation. Among these, there are nineteen sources that are not associated with previously known TeV sources. We have studied fourteen of these sources without known counterparts with VERITAS and Fermi-LAT. VERITAS detected weak gamma-ray emission in the 1~TeV-30~TeV band in the region of DA 495, a pulsar wind nebula coinciding with 2HWC J1953+294, confirming the discovery of the source by HAWC. We did not find any counterpart for the selected fourteen new HAWC sources from our analysis of Fermi-LAT data for energies higher than 10 GeV. During the search, we detected GeV gamma-ray emission coincident with a known TeV pulsar wind nebula, SNR G54.1+0.3 (VER J1930+188), and a 2HWC source, 2HWC J1930+188. The fluxes for isolated, steady sources in the 2HWC catalog are generally in good agreement with those measured by imaging atmospheric Cherenkov telescopes. However, the VERITAS fluxes for SNR G54.1+0.3, DA 495, and TeV J2032+4130 are lower than those measured by HAWC and several new HAWC sources are not detected by VERITAS. This is likely due to a change in spectral shape, source extension, or the influence of diffuse emission in the source region.

Large or even medium sized asteroids impacting the Earth can cause damage on a global scale. Existing and planned concepts for finding near-Earth objects (NEOs) with diameter of 140 m or larger would take ~15-20 years of observation to find ~90% of them. This includes both ground and space based projects. For smaller NEOs (~50-70 m in diameter), the time scale is many decades. The reason it takes so long to detect these objects is because most of the NEOs have highly elliptical orbits that bring them into the inner solar system once per orbit. If these objects cross the Earth's orbit when the Earth is on the other side of the Sun, they will not be detected by facilities on or around the Earth. A constellation of MicroSats in orbit around the Sun can dramatically reduce the time needed to find 90% of NEOs ~100-140 m in diameter.

After two years of operation, the High-Altitude Water Cherenkov (HAWC) Observatory has analyzed the TeV cosmic-ray sky over an energy range between $2.0$ and $72.8$ TeV. The HAWC detector is a ground-based air-shower array located at high altitude in the state of Puebla, Mexico. Using 300 light-tight water tanks, it collects the Cherenkov light from the particles of extensive air showers from primary gamma rays and cosmic rays. This detection method allows for uninterrupted observation of the entire overhead sky (2~sr instantaneous, 8.5~sr integrated) in the energy range from a few TeV to hundreds of TeV. Like other detectors in the northern and southern hemisphere, HAWC observes an energy-dependent anisotropy in the arrival direction distribution of cosmic rays. The observed cosmic-ray anisotropy is dominated by a dipole moment with phase $\alpha\approx40^{\circ}$ and amplitude that slowly rises in relative intensity from $8\times10^{-4}$ at 2 TeV to $14\times10^{-4}$ around 30.3 TeV, above which the dipole decreases in strength. A significant large-scale ($>60^{\circ}$ in angular extent) signal is also observed in the quadrupole and octupole moments, and significant small-scale features are also present, with locations and shapes consistent with previous observations. Compared to previous measurements in this energy range, the HAWC cosmic-ray sky maps improve on the energy resolution and fit precision of the anisotropy. These data can be used in an effort to better constrain local cosmic-ray accelerators and the intervening magnetic fields.

We identified a peculiar proximate sub-damped Lya absorption system (sub-DLA) at z=3.234 with a neutral hydrogen column density of $N_{H} =10^{19.70\pm0.15}~\rm cm^{-2}$ toward the quasar SDSS J125903.26+621211.5 in two epoch optical spectra of the Sloan Digital Sky Survey. We detected Lya residue in the proximate sub-DLA trough at $> 8\sigma$ level. To our surprise, extreme ultraviolet (EUV) continuum emission is significantly ($> 4\sigma$) detected in the corresponding Lymann limit absorption edge at both of the FUV and NUV bands by the Galaxy Evolution Explorer. The high neutral hydrogen column density should allow a negligible transmission of both the Lya line photons and EUV continuum photons due to the high optical depth of the gas. The possible scenarios of foreground galaxy contamination, partial coverage, emission from the quasar host galaxy, and extended Lya emission are excluded in turn, and we speculate that the residual Lya and EUV emission is due to photons scattering (broad Lya and the continuum emission) of electrons residing at a spatial scale larger than that of the proximate sub-DLA. Electron scattering is wavelength-independent, the scattered light is therefore a copy of the incident spectrum that might originate in the accretion disk. With the assistances from the neutral hydrogen absorbers as the "natural coronagraph" and the scatterers as the "natural mirror", we yielded a very hard EUV spectral index of $\alpha = 0.65\pm 0.25$ ($F_{\nu}\propto \nu^{\alpha}$), consistent with the standard picture of the locally heated accretion disk in the inner EUV-emitting radii, as well as in the outer near-infrared-emitting radii suggested by Kishimoto et al. (2008).

We report the discovery of highly reddening and hydrogen Balmer and metastable helium broad absorption lines in the quasar SDSS J145057.28+530007.6, based on the optical and near-infrared spectra taken from the SDSS-III/BOSS and the TripleSpec observations. The nuclear continuum, Balmer decrement and absorption-line depth analyses suggest that (1) the accretion disk is completely obscured and the covering factor of the broad-line region (BLR) is only $0.39\pm0.03$, (2) the power-law continuum is reddened by the SMC extinction law of $E(B-V )=0.72\pm 0.01$ mag and the dusty materials are mainly associated with \ion{Ca}{2} H and K rather than the Balmer and \ion{He}{1}* absorption-line system, (3) the unsaturated Balmer (H$\beta$, H$\gamma$, and H$\delta$) and \ion{He}{1}* $\lambda3889$ absorption lines have same two-Gaussian profiles with the shifts of $-931\pm33$ and $-499\pm39$ km s$^{-1}$ and the widths of $121\pm28$ and $196\pm37$ km s$^{-1}$, respectively. Constrained mutually by the Balmer, \ion{He}{1}* absorption lines and undetected \ion{Fe}{2}* $\lambda5169$ in the photoionization simulations, the physical properties of the outflow gases are derived as follows: ionization parameter $10^{-1.4} \lesssim U \lesssim 10^{-0.8}$, density $10^{8.2\pm0.4} \lesssim n_{\rm H} \lesssim 10^{9.0\pm0.4}$ cm$^{-3}$, and column density $10^{22.0\pm0.2} \lesssim N_{\rm H} \lesssim 10^{22.2-22.3}$ cm$^{-2}$. We propose that the ultra-dense outflow gases appear in the vicinity of the surface of the BLR or are located at most 3.12 pc away from the engine. That probably implies that the outflow originates from the BLR, and this kind of ultra-dense BLR scale outflow gases can effectively test the physical properties of the outer gases of the BLR.

We report on the identification of a new {\gamma}-ray-emitting narrow-line Seyfert 1 (NLS1) galaxy, SDSS J211852.96-073227.5 (hereinafter J2118-0732). The galaxy, at a redshift of 0.26, is associated with a radio source of flat/inverted spectrum at high radio frequencies. The analysis of its optical spectrum obtained in the Sloan Digital Sky Survey (SDSS) revealed a small linewidth of the broad component of the H{\beta} line (full width at half-maximum = 1585 km s^{-1}), making it a radio-loud NLS1 galaxy -- an intriguing class of active galactic nuclei with exceptional multiwavelength properties. A new {\gamma}-ray source centred at J2118-0732 was sporadically detected during 2009--2013 in form of flares by the Fermi-LAT. Our XMM-Newton observations revealed a flat X-ray spectrum described by a simple power law, and a flux variation by a factor of ~2.5 in five months. The source also shows intraday variability in the infrared band. Its broad-band spectral energy distribution can be modelled by emission from a simple one-zone leptonic jet model, and the flux drop from infrared to X-rays in five months can be explained by changes of the jet parameters, though the exact values may be subject to relatively large uncertainties. With the NLS1-blazar composite nucleus, the clear detection of the host galaxy, and the synchronous variations in the multiwavelength fluxes, J2118-0732 provides a new perspective on the formation and evolution of relativistic jets under the regime of relatively small black hole masses and high accretion rates.

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

A new sample of 204 low-mass black holes (LMBHs) in active galactic nuclei (AGNs) is presented with black hole masses in the range of (1-20) * 10^5 M_sun. The AGNs are selected from a systematic search among galaxies in the Seventh Data Release (DR 7) of the Sloan Digital Sky Survey (SDSS), and careful analyses of their optical spectra and precise measurement of spectral parameters. Combining them with our previous sample selected from the SDSS DR 4 makes it the largest LMBH sample so far, totaling over 500 objects. Some of the statistical properties of the combined LMBH AGN sample are briefly discussed, in the context of exploring the low-mass end of the AGN population. Their X-ray luminosities follow the extension of the previously known correlation with the [O III] luminosity. The effective optical-to-X-ray spectral indices \alpha_OX, albeit with a large scatter, are broadly consistent with the extension of the relation with the near-UV luminosity L_2500\AA. Interestingly, a correlation of \alpha_OX with black hole mass is also found in the sense that \alpha_OX is statistically flatter (stronger X-ray relative to optical) for lower black hole mass. Only 26 objects, mostly radio loud, were detected in radio at 20 cm in the FIRST survey, giving a radio loud fraction of 4%. The host galaxies of LMBHs have stellar masses in the range of 10^8.8-10^12.4 M_sun and optical colors typical of Sbc spirals. They are dominated by young stellar populations that seem to have undergone a continuous star formation history.

In this paper, we present an intercomparison study of two quasars, SDSS J145618.32+340037.2 and SDSS J215331.50-025514.1, which have proximate damped Lya systems (PDLAs) with residual flux upon the Lya absorption trough. Though they both have residual flux as luminous as $10^{43}erg/s$, their PDLAs are quite different in, e.g., HI column density, metal line absorption strength, high-ionization absorption lines as well as residual flux strength. For J1456, the HI column density is $20.6\pm0.2$, with za=2.3138, nearly identical to the quasar redshift (z=2.3142) determined from the [OIII] emission line. The metallicity of this system is typical of DLAs and there is high ionization therein, suggesting that the PDLA system is multiphase, putting it in the quasar environment. For J2153, we measure the HI column density to be $21.5\pm0.1$ at za=3.511, slightly redshifted with respect to the quasar (z=3.490) measured from CIII]. The metallicity of this system is quite low and there is a lack of significant high-ionization absorption lines therein, suggesting that the system is beyond the quasar host galaxy. The residual flux is wide (1000 km/s) in J1456, with a significance of $8\sigma$, while also wide (1500 km/s) but with a smaller significance of $3\sigma$ in J2153. Among many explanations, we find that Lya fuzz or resonant scattering can be used to explain the residual flux in the two sources while partial coverage cannot be excluded for J1456. By comparing these two cases, together with similar cases reported previously, we suggest that the strength of the residual flux is related to properties such as metallicity and high-ionization absorption lines of PDLAs. The residual flux recorded upon the PDLA absorption trough opens a window for us to see the physical conditions and processes of the quasar environment, and their profile and strength further remind us of their spatial scales.

An indirect measurement of the antiproton flux in cosmic rays is possible as the particles undergo deflection by the geomagnetic field. This effect can be measured by studying the deficit in the flux, or shadow, created by the Moon as it absorbs cosmic rays that are headed towards the Earth. The shadow is displaced from the actual position of the Moon due to geomagnetic deflection, which is a function of the energy and charge of the cosmic rays. The displacement provides a natural tool for momentum/charge discrimination that can be used to study the composition of cosmic rays. Using 33 months of data comprising more than 80 billion cosmic rays measured by the High Altitude Water Cherenkov (HAWC) observatory, we have analyzed the Moon shadow to search for TeV antiprotons in cosmic rays. We present our first upper limits on the $\bar{p}/p$ fraction, which in the absence of any direct measurements, provide the tightest available constraints of $\sim1\%$ on the antiproton fraction for energies between 1 and 10 TeV.

The Solar Gravitational Lens (SGL) allows for major brightness amplification ($\sim 10^{11}$ at wavelength of $1~\mu$m) and extreme angular resolution ($\sim10^{-10}$ arcsec) within a narrow field of view. A meter-class telescope, with a modest coronagraph to block solar light with 1e-6 suppression placed in the focal area of the SGL, can image an exoplanet at a distance of 30 parsec with few kilometer-scale resolution on its surface. Notably, spectroscopic broadband SNR is $\sim 10^{-6}$ in two weeks of integration time, providing this instrument with incredible remote sensing capabilities. A mission capable of exploiting the remarkable optical properties of the SGL allows for direct high-resolution imaging/spectroscopy of a potentially habitable exoplanet. Such missions could allow exploration of exoplanets relying on the SGL capabilities decades, if not centuries, earlier than possible with other extant technologies.

The remarkable optical properties of the solar gravitational lens (SGL) include major brightness amplification (~1e11 at wavelength of 1 um) and extreme angular resolution (~1e-10 arcsec) in a narrow field of view. A mission to the SGL carrying a modest telescope and coronagraph opens up a possibility for direct megapixel imaging and high-resolution spectroscopy of a habitable Earth-like exoplanet at a distance of up to 100 light years. The entire image of such a planet is compressed by the SGL into a region with a diameter of ~1.3 km in the vicinity of the focal line. The telescope, acting as a single pixel detector while traversing this region, can build an image of the exoplanet with kilometer-scale resolution of its surface, enough to see its surface features and signs of habitability. We report here on the results of our initial study of a mission to the deep outer regions of our solar system, with the primary mission objective of conducting direct megapixel high-resolution imaging and spectroscopy of a potentially habitable exoplanet by exploiting the remarkable optical properties of the SGL. Our main goal was to investigate what it takes to operate spacecraft at such enormous distances with the needed precision. Specifically, we studied i) how a space mission to the focal region of the SGL may be used to obtain high-resolution direct imaging and spectroscopy of an exoplanet by detecting, tracking, and studying the Einstein ring around the Sun, and ii) how such information could be used to detect signs of life on another planet. Our results indicate that a mission to the SGL with an objective of direct imaging and spectroscopy of a distant exoplanet is challenging, but possible. We composed a list of recommendations on the mission architectures with risk and return tradeoffs and discuss an enabling technology development program.

Proxima Centauri is known as the closest star from the Sun. Recently, radial velocity observations revealed the existence of an Earth-mass planet around it. With an orbital period of ~11 days, the surface of Proxima Centauri b is temperate and might be habitable. We took a photometric monitoring campaign to search for its transit, using the Bright Star Survey Telescope at the Zhongshan Station in Antarctica. A transit-like signal appearing on 2016 September 8th, is identified tentatively. Its midtime, $T_{C}=2,457,640.1990\pm0.0017$ HJD, is consistent with the predicted ephemeris based on RV orbit in a 1$\sigma$ confidence interval. Time-correlated noise is pronounced in the light curve of Proxima Centauri, affecting detection of transits. We develop a technique, in a Gaussian process framework, to gauge the statistical significance of potential transit detection. The tentative transit signal reported here, has a confidence level of $2.5\sigma$. Further detection of its periodic signals is necessary to confirm the planetary transit of Proxima Centauri b. We plan to monitor Proxima Centauri in next Polar night at Dome A in Antarctica, taking the advantage of continuous darkness. \citet{Kipping17} reported two tentative transit-like signals of Proxima Centauri b, observed by the Microvariability and Oscillation of Stars space Telescope in 2014 and 2015, respectively. The midtransit time of our detection is 138 minutes later than that predicted by their transit ephemeris. If all the signals are real transits, the misalignment of the epochs plausibly suggests transit timing variations of Proxima Centauri b induced by an outer planet in this system.

The High Altitude Water Cherenkov observatory (HAWC) is an air shower array devised for TeV gamma-ray astronomy. HAWC is located at an altitude of 4100 m a.s.l. in Sierra Negra, Mexico. HAWC consists of 300 Water Cherenkov Detectors, each instrumented with 4 photomultiplier tubes (PMTs). HAWC re-uses the Front-End Boards from the Milagro experiment to receive the PMT signals. These boards are used in combination with Time to Digital Converters (TDCs) to record the time and the amount of light in each PMT hit (light flash). A set of VME TDC modules (128 channels each) is operated in a continuous (dead time free) mode. The TDCs are read out via the VME bus by Single-Board Computers (SBCs), which in turn are connected to a gigabit Ethernet network. The complete system produces ~ 500 MB/s of raw data. A high-throughput data processing system has been designed and built to enable real-time data analysis. The system relies on off-the-shelf hardware components, an open-source software technology for data transfers (ZeroMQ) and a custom software framework for data analysis (AERIE). Multiple trigger and reconstruction algorithms can be combined and run on blocks of data in a parallel fashion, producing a set of output data streams which can be analyzed in real time with minimal latency (< 5 s). This paper provides an overview of the hardware set-up and an in-depth description of the software design, covering both the TDC data acquisition system and the real-time data processing system. The performance of these systems is also discussed.

On February 17 2016, the IceCube real-time neutrino search identified, for the first time, three muon neutrino candidates arriving within 100 s of one another, consistent with coming from the same point in the sky. Such a triplet is expected once every 13.7 years as a random coincidence of background events. However, considering the lifetime of the follow-up program the probability of detecting at least one triplet from atmospheric background is 32%. Follow-up observatories were notified in order to search for an electromagnetic counterpart. Observations were obtained by Swift's X-ray telescope, by ASAS-SN, LCO and MASTER at optical wavelengths, and by VERITAS in the very-high-energy gamma-ray regime. Moreover, the Swift BAT serendipitously observed the location 100 s after the first neutrino was detected, and data from the Fermi LAT and HAWC observatory were analyzed. We present details of the neutrino triplet and the follow-up observations. No likely electromagnetic counterpart was detected, and we discuss the implications of these constraints on candidate neutrino sources such as gamma-ray bursts, core-collapse supernovae and active galactic nucleus flares. This study illustrates the potential of and challenges for future follow-up campaigns.

We present the cold neutral content (H I and C I gas) of 13 quasar 2175 \AA$ $ dust absorbers (2DAs) at $z$ = 1.6 - 2.5 to investigate the correlation between the presence of the UV extinction bump with other physical characteristics. These 2DAs were initially selected from the Sloan Digital Sky Surveys I - III and followed up with the Keck-II telescope and the Multiple Mirror Telescope as detailed in our Paper I. We perform a correlation analysis between metallicity, redshift, depletion level, velocity width, and explore relationships between 2DAs and other absorption line systems. The 2DAs on average have higher metallicity, higher depletion levels, and larger velocity widths than Damped Lyman-$\alpha$ absorbers (DLAs) or subDLAs. The correlation between [Zn/H] and [Fe/Zn] or [Zn/H] and log$\Delta$V$_{90}$ can be used as alternative stellar mass estimators based on the well-established mass-metallicity relation. The estimated stellar masses of the 2DAs in this sample are in the range of $\sim$ 10$^9$ to $\sim$2 $\times$ 10$^{11}$ $M_{\odot}$ with a median value of $\sim$2 $\times$ 10$^{10}$ $M_{\odot}$. The relationship with other quasar absorption line systems can be described as (1) 2DAs are a subset of Mg II and Fe II absorbers, (2) 2DAs are preferentially metal-strong DLAs/subDLAs, (3) More importantly, all of the 2DAs show C I detections with logN(C I) $>$ 14.0 cm$^{-2}$, (4) 2DAs can be used as molecular gas tracers. Their host galaxies are likely to be chemically enriched, evolved, massive (more massive than typical DLA/subDLA galaxies), and presumably star-forming galaxies.

We investigate the hypothesized tidal disruption of the Hercules ultra-faint dwarf galaxy (UFD). Previous tidal disruption studies of the Hercules UFD have been hindered by the high degree of foreground contamination in the direction of the dwarf. We bypass this issue by using RR Lyrae stars, which are standard candles with a very low field-volume density at the distance of Hercules. We use wide-field imaging from the Dark Energy Camera on CTIO to identify candidate RR Lyrae stars, supplemented with observations taken in coordination with the Beijing-Arizona Sky Survey on the Bok Telescope. Combining color, magnitude, and light-curve information, we identify three new RR Lyrae stars associated with Hercules. All three of these new RR Lyrae stars lie outside its published tidal radius. When considered with the nine RR Lyrae stars already known within the tidal radius, these results suggest that a substantial fraction of Hercules' stellar content has been stripped. With this degree of tidal disruption, Hercules is an interesting case between a visibly disrupted dwarf (such as the Sagittarius dwarf spheroidal galaxy) and one in dynamic equilibrium. The degree of disruption also shows that we must be more careful with the ways we determine object membership when estimating dwarf masses in the future. One of the three discovered RR Lyrae stars sits along the minor axis of Hercules, but over two tidal radii away. This type of debris is consistent with recent models that suggest Hercules' orbit is aligned with its minor axis.

The unexpectedly high flux of cosmic ray positrons detected at Earth may originate from nearby astrophysical sources, dark matter, or unknown processes of cosmic-ray secondary production. We report the detection, using the HighAltitude Water Cherenkov Observatory (HAWC), of extended tera-electron volt gamma-ray emission coincident with the locations of two nearby middle-aged pulsars (Geminga and PSR B0656+14). The HAWC observations demonstrate that these pulsars are indeed local sources of accelerated leptons, but the measured tera-electron volt emission profile constrains the diffusion of particles away from these sources to be much slower than previously assumed. We demonstrate that the leptons emitted by these objects are therefore unlikely to be the origin of the excess positrons, which may have a more exotic origin.

The High Altitude Water Cherenkov (HAWC) gamma-ray observatory is a wide field-of-view observatory sensitive to 500 GeV - 100 TeV gamma rays and cosmic rays. With its observations over 2/3 of the sky every day, the HAWC observatory is sensitive to a wide variety of astrophysical sources, including possible gamma rays from dark matter. Dark matter annihilation and decay in the Milky Way Galaxy should produce gamma-ray signals across many degrees on the sky. The HAWC instantaneous field-of-view of 2 sr enables observations of extended regions on the sky, such as those from dark matter in the Galactic halo. Here we show limits on the dark matter annihilation cross-section and decay lifetime from HAWC observations of the Galactic halo with 15 months of data. These are some of the most robust limits on TeV and PeV dark matter, largely insensitive to the dark matter morphology. These limits begin to constrain models in which PeV IceCube neutrinos are explained by dark matter which primarily decays into hadrons.

We report on the measurement of the all-particle cosmic ray energy spectrum with the High Altitude Water Cherenkov (HAWC) Observatory in the energy range 10 to 500 TeV. HAWC is a ground based air-shower array deployed on the slopes of Volcan Sierra Negra in the state of Puebla, Mexico, and is sensitive to gamma rays and cosmic rays at TeV energies. The data used in this work were taken from 234 days between June 2016 to February 2017. The primary cosmic-ray energy is determined with a maximum likelihood approach using the particle density as a function of distance to the shower core. Introducing quality cuts to isolate events with shower cores landing on the array, the reconstructed energy distribution is unfolded iteratively. The measured all-particle spectrum is consistent with a broken power law with an index of $-2.49\pm0.01$ prior to a break at $(45.7\pm0.1$) TeV, followed by an index of $-2.71\pm0.01$. The spectrum also respresents a single measurement that spans the energy range between direct detection and ground based experiments. As a verification of the detector response, the energy scale and angular resolution are validated by observation of the cosmic ray Moon shadow's dependence on energy.

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

The positron excess measured by PAMELA and AMS can only be explained if there is one or several sources injecting them. Moreover, at the highest energies, it requires the presence of nearby ($\sim$hundreds of parsecs) and middle age (maximum of $\sim$hundreds of kyr) sources. Pulsars, as factories of electrons and positrons, are one of the proposed candidates to explain the origin of this excess. To calculate the contribution of these sources to the electron and positron flux at the Earth, we developed EDGE (Electron Diffusion and Gamma rays to the Earth), a code to treat the propagation of electrons and compute their diffusion from a central source with a flexible injection spectrum. Using this code, we can derive the source's gamma-ray spectrum, spatial extension, the all-electron density in space, the electron and positron flux reaching the Earth and the positron fraction measured at the Earth. We present in this paper the foundations of the code and study how different parameters affect the gamma-ray spectrum of a source and the electron flux measured at the Earth. We also studied the effect of several approximations usually performed in these studies.

Galactic diffuse TeV gamma-ray emission is produced by the interaction of high-energy cosmic-ray particles with matter and radiation in our Galaxy. The measurement of Galactic diffuse TeV gamma-ray emission would provide strong constraints on the acceleration and propagation of Galactic cosmic rays. The High Altitude Water Cherenkov (HAWC) Observatory, located in central Mexico at 4100 m above sea level, is sensitive to gamma rays between a few hundreds GeV and about 100 TeV. Thanks to its large field of view of 2 steradians and excellent background estimation, HAWC has a unique capability to observe large extended sources such as Galactic diffuse emission. We present the analysis strategy to detect the Galactic diffuse emission with HAWC, including techniques to remove contaminations from localized sources.

Nearby electron/positron accelerators, mostly Pulsar Wind Nebulae (PWNe), have been proposed as potential origins of the positron excess above 10 GeV. The HAWC Observatory reveals two very extended sources spatially coincident with two nearby middle-aged pulsars: Geminga and PSR B0656+14, suggesting ultrarelativistic electrons/positrons accelerated in our backyard. Morphological studies on these two PWNe provide a constraint on the diffusion coefficient at HAWC energies. In this poster, we will present the model development and morphological studies on these PWNe, and the derived diffusion coefficient that best fits the data.

The HAWC Gamma-Ray Observatory has reported the discovery of TeV gamma-ray emission extending several degrees around the positions of Geminga and B0656+14 pulsars. Assuming these gamma rays are produced by inverse Compton scattering off low-energy photons in electron halos around the pulsars, we determine the diffusion of electrons and positrons in the local interstellar medium. We will present the morphological and spectral studies of these two VHE gamma-ray sources and the derived positron spectrum at Earth.

List of proceedings from the HAWC Collaboration presented at the 35th International Cosmic Ray Conference, 12 July - 20 July 2017, Bexco, Busan, Korea.

The rest-frame UV spectra of three recent tidal disruption events (TDEs), ASASSN-14li, PTF15af and iPTF16fnl display strong nitrogen emission lines but weak or undetectable carbon lines. In these three objects, the upper limits of the C III] $\lambda 1908$/N III] $\lambda 1750$ ratio are about two orders of magnitude lower than those of quasars, suggesting a high abundance ratio of [N/C]. With detailed photoionization simulations, we demonstrate that $C^{2+}$ and $N^{2+}$ are formed in the same zone, so the C III]/N III] ratio depends only moderately on the physical conditions in the gas and weakly on the shape of the ionizing continuum. There are smaller than $0.5$ dex variations in the line ratio over wide ranges of gas densities and ionization parameters at a given metallicity. This allows a robust estimate of the relative abundance ratio nitrogen to carbon. We derive a relative abundance ratio of [N/C]$>1.5$ for ASASSN-14li, and even higher for PTF15af and iPTF16fnl. This suggests that the broad line region in those TDE sources is made of nitrogen-enhanced core material that falls back at later times. Based on stellar evolution models, the lower limit of the disrupted star should be larger than 0.6M$_{\odot}$. The chemical abundance of the line emitting gas gives a convincing evidence that the flares origin from stellar tidal disruptions. The coincidence of the weakness of the X-ray emission with the strong broad absorption lines in PTF15af, iPTF16fnl and the strong X-ray emission without such lines in ASASSN-li14 are in analogy to quasars with and without broad absorption lines.

The positron excess measured by PAMELA and AMS can only be explained if there is one or several sources injecting them. Moreover, at the highest energies, it requires the presence of nearby ($\sim$hundreds of parsecs) and middle age (maximum of $\sim$hundreds of kyr) source. Pulsars, as factories of electrons and positrons, are one of the proposed candidates to explain the origin of this excess. To calculate the contribution of these sources to the electron and positron flux at the Earth, we developed EDGE (Electron Diffusion and Gamma rays to the Earth), a code to treat diffusion of electrons and compute their diffusion from a central source with a flexible injection spectrum. We can derive the source's gamma-ray spectrum, spatial extension, the all-electron density in space and the electron and positron flux reaching the Earth. We present in this contribution the fundamentals of the code and study how different parameters affect the gamma-ray spectrum of a source and the electron flux measured at the Earth.

The High Altitude Water Cherenkov (HAWC) Gamma-ray Observatory is an extensive air shower detector operating in central Mexico, which has recently completed its first two years of full operations. If for a burst like GRB 130427A at a redshift of 0.34 and a high-energy component following a power law with index -1.66, the high-energy component is extended to higher energies with no cut-off other than from extragalactic background light attenuation, HAWC would observe gamma rays with a peak energy of $\sim$300 GeV. This paper reports the results of HAWC observations of 64 gamma-ray bursts (GRBs) detected by $\mathit{Swift}$ and $\mathit{Fermi}$, including three GRBs that were also detected by the Large Area Telescope ($\mathit{Fermi}$-LAT). An ON/OFF analysis method is employed, searching on the time scale given by the observed light curve at keV-MeV energies and also on extended time scales. For all GRBs and time scales, no statistically significant excess of counts is found and upper limits on the number of gamma rays and the gamma-ray flux are calculated. GRB 170206A, the third brightest short GRB detected by the Gamma-ray Burst Monitor on board the $\mathit{Fermi}$ satellite ($\mathit{Fermi}$-GBM) and also detected by the LAT, occurred very close to zenith. The LAT measurements can neither exclude the presence of a synchrotron self-Compton (SSC) component nor constrain its spectrum. Instead, the HAWC upper limits constrain the expected cut-off in an additional high-energy component to be less than $100~\rm{GeV}$ for reasonable assumptions about the energetics and redshift of the burst.

We report the detection of heavy reddening and the He I* $\lambda10830$ absorption lines at the AGNs' redshift in three Narrow-line Seyfert 1 galaxies: SDSS J091848.61+211717.0, SDSS J111354.66+124439.0, and SDSS J122749.13+321458.9. They exhibit very red optical to near-infrared colors, narrow Balmer/Paschen broad emission lines and He I* $\lambda10830$ absorption lines. The ultraviolet-optical-infrared nucleus continua are reddened by the SMC extinction law of $E(B-V)\sim 0.74$, 1.17, and 1.24 mag for three objects, which are highly consistent with the values obtained from the broad-line Balmer decrements, but larger than those of narrow emission lines. The reddening analysis suggests the extinction dust simultaneously obscure the accretion disk, the broad emission line region, and the hot dust from the inner edge of the torus. It is possible that the dust obscuring the AGN structures are the dusty torus itself. Furthermore, the Cloudy analysis of the He I* $\lambda10830$ absorption lines propose the distance of the absorption materials to be the extend scale of the torus. That greatly increases probabilities of the obscure and absorption materials being the dusty torus.

The High Altitude Water Cherenkov (HAWC) gamma-ray observatory is a wide field of view observatory sensitive to 500 GeV - 100 TeV gamma rays and cosmic rays. It can also perform diverse indirect searches for dark matter (DM) annihilation and decay. Among the most promising targets for the indirect detection of dark matter are dwarf spheroidal galaxies. These objects are expected to have few astrophysical sources of gamma rays but high dark matter content, making them ideal candidates for an indirect dark matter detection with gamma rays. Here we present individual limits on the annihilation cross section and decay lifetime for 15 dwarf spheroidal galaxies within the HAWC field-of-view, as well as their combined limit. These are the first limits on the annihilation cross section and decay lifetime using data collected with HAWC.

We present the development of a real-time flare monitor for the High Altitude Water Cherenkov (HAWC) observatory. The flare monitor has been fully operational since 2017 January and is designed to detect very high energy (VHE; $E\gtrsim100$ GeV) transient events from blazars on time scales lasting from 2 minutes to 10 hours in order to facilitate multiwavelength and multimessenger studies. These flares provide information for investigations into the mechanisms that power the blazars' relativistic jets and accelerate particles within them, and they may also serve as probes of the populations of particles and fields in intergalactic space. To date, the detection of blazar flares in the VHE range has relied primarily on pointed observations by imaging atmospheric Cherenkov telescopes. The recently completed HAWC observatory offers the opportunity to study VHE flares in survey mode, scanning 2/3 of the entire sky every day with a field of view of $\sim$1.8 steradians. In this work, we report on the sensitivity of the HAWC real-time flare monitor and demonstrate its capabilities via the detection of three high-confidence VHE events in the blazars Markarian 421 and Markarian 501.

We present a search of very high energy gamma-ray emission from the Northern $\textit{Fermi}$ Bubble region using data collected with the High Altitude Water Cherenkov (HAWC) gamma-ray observatory. The size of the data set is 290 days. No significant excess is observed in the Northern $\textit{Fermi}$ Bubble region, hence upper limits above $1\,\text{TeV}$ are calculated. The upper limits are between $3\times 10^{-7}\,\text{GeV}\, \text{cm}^{-2}\, \text{s}^{-1}\,\text{sr}^{-1}$ and $4\times 10^{-8}\,\text{GeV}\,\text{cm}^{-2}\,\text{s}^{-1}\,\text{sr}^{-1}$. The upper limits disfavor a proton injection spectrum that extends beyond $100\,\text{TeV}$ without being suppressed. They also disfavor a hadronic injection spectrum derived from neutrino measurements.

We present results from daily monitoring of gamma rays in the energy range $\sim0.5$ to $\sim100$ TeV with the first 17 months of data from the High Altitude Water Cherenkov (HAWC) Observatory. Its wide field of view of 2 steradians and duty cycle of $>95$% are unique features compared to other TeV observatories that allow us to observe every source that transits over HAWC for up to $\sim6$ hours each sidereal day. This regular sampling yields unprecedented light curves from unbiased measurements that are independent of seasons or weather conditions. For the Crab Nebula as a reference source we find no variability in the TeV band. Our main focus is the study of the TeV blazars Markarian (Mrk) 421 and Mrk 501. A spectral fit for Mrk 421 yields a power law index $\Gamma=2.21 \pm0.14_{\mathrm{stat}}\pm0.20_{\mathrm{sys}}$ and an exponential cut-off $E_0=5.4 \pm 1.1_{\mathrm{stat}}\pm 1.0_{\mathrm{sys}}$ TeV. For Mrk 501, we find an index $\Gamma=1.60\pm 0.30_{\mathrm{stat}} \pm 0.20_{\mathrm{sys}}$ and exponential cut-off $E_0=5.7\pm 1.6_{\mathrm{stat}} \pm 1.0_{\mathrm{sys}}$ TeV. The light curves for both sources show clear variability and a Bayesian analysis is applied to identify changes between flux states. The highest per-transit fluxes observed from Mrk 421 exceed the Crab Nebula flux by a factor of approximately five. For Mrk 501, several transits show fluxes in excess of three times the Crab Nebula flux. In a comparison to lower energy gamma-ray and X-ray monitoring data with comparable sampling we cannot identify clear counterparts for the most significant flaring features observed by HAWC.

The first limits on the prompt emission from the long gamma-ray burst (GRB) 130427A in the $>100\nobreakspace\rm{GeV}$ energy band are reported. GRB 130427A was the most powerful burst ever detected with a redshift $z\lesssim0.5$ and featured the longest lasting emission above $100\nobreakspace\rm{MeV}$. The energy spectrum extends at least up to $95\nobreakspace\rm{GeV}$, clearly in the range observable by the High Altitude Water Cherenkov (HAWC) Gamma-ray Observatory, a new extensive air shower detector currently under construction in central Mexico. The burst occurred under unfavourable observation conditions, low in the sky and when HAWC was running 10% of the final detector. Based on the observed light curve at MeV-GeV energies, eight different time periods have been searched for prompt and delayed emission from this GRB. In all cases, no statistically significant excess of counts has been found and upper limits have been placed. It is shown that a similar GRB close to zenith would be easily detected by the full HAWC detector, which will be completed soon. The detection rate of the full HAWC detector may be as high as one to two GRBs per year. A detection could provide important information regarding the high energy processes at work and the observation of a possible cut-off beyond the $\mathit{Fermi}$-LAT energy range could be the signature of gamma-ray absorption, either in the GRB or along the line of sight due to the extragalactic background light.

Galactic dynamo models sustained by supernova (SN) driven turbulence and differential rotation have revealed that the sustenance of large scale fields requires a flux of small scale magnetic helicity to be viable. Here we generalize a minimalist analytic version of such galactic dynamos to explore some heretofore unincluded contributions from shear on the total turbulent energy and turbulent correlation time, with the helicity fluxes maintained by either winds, diffusion, or magnetic buoyancy. We construct an analytic framework for modeling the turbulent energy and correlation time as functions of SN rate and shear. We compare our prescription with previous approaches that only include rotation. The solutions depend separately on the rotation period and the eddy turnover time and not just on their ratio (the Rossby number). We consider models in which these two time scales are allowed to be independent and also a case in which they are mutually dependent on radius when a radial dependent SN rate model is invoked. For the case of a fixed rotation period (or fixed radius) we show that the influence of shear is dramatic for low Rossby numbers, reducing the correlation time of the turbulence, which in turn, strongly reduces the saturation value of the dynamo compared to the case when the shear is ignored. We also show that even in the absence of winds or diffusive fluxes, magnetic buoyancy may be able to sustain sufficient helicity fluxes to avoid quenching.

RXJ 1301.9+2747 is an optically identified very low mass AGN candidate with M_BH~1x10^6M_sun, which shows extremely soft X-ray emission and unusual X-ray variability in the form of short-lived flares. We present an analysis of multiwavelength observations of RXJ 1301.9+2747 in order to study the properties of the active nucleus and its host galaxy. The UV-to-X-ray spectrum in the quiescent state can be well and self-consistently described by a thermal and a Comptonized emission from accretion disk, with blackbody dominating ~70% of the X-rays in the 0.2-2 keV. The same model can describe the X-ray spectrum in the flare state but the Comptonized component becomes dominant (~80%). The best-fit implies an Eddington ratio of ~0.14, and a black hole mass of 1.7-2.8x10^6M_sun, in agreement with the estimation from the optical data within errors. However, the best-fitting model under-predicts the optical flux for the HST point source by a factor of ~2. The excess of nuclear optical emission could be attributed to a nuclear stellar cluster which is frequently seen in low mass AGNs. The X-ray to optical spectral slope (alpha_ox) is lower than in most other active galaxies, which may be attributed to intrinsically X-ray weakness due to very little hot and optically thin coronal emission. We performed a pilot search for weak or hidden broad emission lines using optical spectropolarimetry observations, but no any polarized broad lines are detected. The host galaxy appears to be a disk galaxy with a boxy pseudobulge or nuclear bar accounting for ~15% of the total starlight, which is consistent with the general characteristics of the host of low mass AGNs.

Very few low-ionization broad absorption line (LoBAL) QSOs have been found at high redshifts to date. One high-redshift LoBAL QSO, J0122+1216, was recently discovered at the Lijiang 2.4-m Telescope with an initial redshift determination of 4.76. Aiming to investigate its physical properties, we carried out follow-up observations in the optical and near-IR spectroscopy. Near-IR spectra from UKIRT and P200 confirms that it is a LoBAL, with a new redshift determination of $4.82\pm0.01$ based on the \mgii~ emission-line. The new \mgii~ redshift determination reveals strong blueshifts and asymmetry of the high-ionization emission lines. We estimated a black hole mass of $\sim 2.3\times 10^9 M_\odot$ and Eddington ratio of $\sim 1.0$ according to the empirical \mgii-based single-epoch relation and bolometric correction factor. It is possible that strong outflows are the result of an extreme quasar environment driven by the high Eddington ratio. A lower limit on the outflowing kinetic power ($>0.9\% L_{Edd}$) was derived from both emission and absorption lines, indicating these outflows play a significant role in the feedback process to regulate the growth of its black hole as well as host galaxy evolution.

In this paper we present an analysis of absorption line variability in mini-BAL quasar LBQS 1206+1052. The SDSS spectrum demonstrates that the absorption troughs can be divided into two components of blueshift velocities of $\sim$700 km s$^{-1}$ and $\sim$1400 km s$^{-1}$ relative to the quasar rest-frame. The former component shows rare Balmer absorption, which is an indicator of high density absorbing gas, thus the quasar is worth follow-up spectroscopic observations. Our follow-up optical and near-infrared spectra using MMT, YFOSC, TripleSpec and DBSP reveal that the strengths of the absorption lines vary for both of the two components, while the velocities do not change. We reproduce all of the spectral data by assuming that only the ionization state of the absorbing gas is variable and that all other physical properties are invariable. The variation of ionization is consistent with the variation of optical continuum from the V-band light-curve. Additionally, we can not interpret the data by assuming that the variability is due to a movement of the absorbing gas. Therefore, our analysis strongly indicates that the absorption line variability in LBQS 1206+1052 is photoionization-driven. As shown from photo-ionization simulations, the absorbing gas with blueshift velocity of $\sim$700 km s$^{-1}$ has a density in the range of $10^9$ to $10^{10}$ cm$^{-3}$ and a distance of $\sim$1 pc, and the gas with blueshift velocity of $\sim$1400 km s$^{-1}$ has a density of $10^3$ cm$^{-3}$ and a distance of $\sim$1 kpc.

We present a statistical analysis of the variability of broad absorption lines (BALs) in quasars using the large multi-epoch spectroscopic dataset of the Sloan Digital Sky Survey Data Release 12 (SDSS DR12). We divide the sample into two groups according to the pattern of the variation of C iv BAL with respect to that of continuum: the equivalent widths (EW) of the BAL decreases (increases) when the continuum brightens (dims) as group T1; and the variation of EW and continuum in the opposite relation as group T2. We find that T2 has significantly (P_T<10-6 , Students T Test) higher EW ratios (R) of Si iv to C iv BAL than T1. Our result agrees with the prediction of photoionization models that C +3 column density increases (decreases) if there is a (or no) C +3 ionization front while R decreases with the incident continuum. We show that BAL variabilities in at least 80% quasars are driven by the variation of ionizing continuum while other models that predict uncorrelated BAL and continuum variability contribute less than 20%. Considering large uncertainty in the continuum flux calibration, the latter fraction may be much smaller. When the sample is binned into different time interval between the two observations, we find significant difference in the distribution of R between T1 and T2 in all time-bins down to a deltaT < 6 days, suggesting that BAL outflow in a fraction of quasars has a recombination time scale of only a few days.

We present the first catalog of TeV gamma-ray sources realized with the recently completed High Altitude Water Cherenkov Observatory (HAWC). It is the most sensitive wide field-of-view TeV telescope currently in operation, with a 1-year survey sensitivity of ~5-10% of the flux of the Crab Nebula. With an instantaneous field of view >1.5 sr and >90% duty cycle, it continuously surveys and monitors the sky for gamma ray energies between hundreds GeV and tens of TeV. HAWC is located in Mexico at a latitude of 19 degree North and was completed in March 2015. Here, we present the 2HWC catalog, which is the result of the first source search realized with the complete HAWC detector. Realized with 507 days of data and represents the most sensitive TeV survey to date for such a large fraction of the sky. A total of 39 sources were detected, with an expected contamination of 0.5 due to background fluctuation. Out of these sources, 16 are more than one degree away from any previously reported TeV source. The source list, including the position measurement, spectrum measurement, and uncertainties, is reported. Seven of the detected sources may be associated with pulsar wind nebulae, two with supernova remnants, two with blazars, and the remaining 23 have no firm identification yet.

Broad emission-line outflows of active galactic nuclei (AGNs) have been proposed for many years but are very difficult to quantitatively study because of the coexistence of the gravitationally-bound and outflow emission. We present detailed analysis of a heavily reddened quasar, SDSS J000610.67+121501.2, whose normal ultraviolet (UV) broad emission lines (BELs) are heavily suppressed by the Dusty Torus as a natural "Coronagraph", thus the blueshifted BELs (BBELs) can be reliably measured. The physical properties of the emission-line outflows are derived as follows: ionization parameter $U \sim 10^{-0.5}$, column density $N_{\rm H}\sim 10^{22.0}$ cm$^{-2}$, covering fraction of $\sim 0.1$ and upper limit density of $n_{\rm H}\sim 10^{5.8}$ cm$^{-3}$. The outflow gases are located at least 41 pc away from the central engine, which suggests that they have expanded to the scale of the dust torus or beyond. Besides, Lya shows a narrow symmetric component, to our surprise, which is undetected in any other lines. After inspecting the narrow emission-line region and the starforming region as the origin of the Lya narrow line, we propose the end-result of outflows, diffusing gases in the larger region, acts as the screen of Lya photons. Future high spatial resolution spectrometry and/or spectropolarimetric observation are needed to make a final clarification.

The Crab Nebula is the brightest TeV gamma-ray source in the sky and has been used for the past 25 years as a reference source in TeV astronomy, for calibration and verification of new TeV instruments. The High Altitude Water Cherenkov Observatory (HAWC), completed in early 2015, has been used to observe the Crab Nebula at high significance across nearly the full spectrum of energies to which HAWC is sensitive. HAWC is unique for its wide field-of-view, nearly 2 sr at any instant, and its high-energy reach, up to 100 TeV. HAWC's sensitivity improves with the gamma-ray energy. Above $\sim$1 TeV the sensitivity is driven by the best background rejection and angular resolution ever achieved for a wide-field ground array. We present a time-integrated analysis of the Crab using 507 live days of HAWC data from 2014 November to 2016 June. The spectrum of the Crab is fit to a function of the form $\phi(E) = \phi_0 (E/E_{0})^{-\alpha -\beta\cdot{\rm{ln}}(E/E_{0})}$. The data is well-fit with values of $\alpha=2.63\pm0.03$, $\beta=0.15\pm0.03$, and log$_{10}(\phi_0~{\rm{cm}^2}~{\rm{s}}~{\rm{TeV}})=-12.60\pm0.02$ when $E_{0}$ is fixed at 7 TeV and the fit applies between 1 and 37 TeV. Study of the systematic errors in this HAWC measurement is discussed and estimated to be $\pm$50\% in the photon flux between 1 and 37 TeV. Confirmation of the Crab flux serves to establish the HAWC instrument's sensitivity for surveys of the sky. The HAWC survey will exceed sensitivity of current-generation observatories and open a new view of 2/3 of the sky above 10 TeV.

SDSS J072910.34+333634.3 is a reddened quasar at z=0.96. The archivel Keck/ESI spectrum and our new P200/TripleSpec spectrum reveal an absorption line system in He I*, Ca II H\&K and Na ID. The absorption line system has a width of $\sim$600 km/s and a blueshift velocity of $\sim$800 km/s relative to the core of narrow emission lines, indicating an outflow. Using the Ca II doublet, we determined that the outflowing gas covers $\sim$70\% of the continuum. On the other hand, the HST/ACS image which taken in rest-frame 4130 \AA\ show that the fraction of the quasar emission in ESI aperture was $<$40\%. We thus conclude that the absorbing gas covers a significant fraction of extended starlight emission, and the best-estimated fraction of $\sim$50\% yields a lower limit of the crosssectional area of the outflowing gas to be $>$8 kpc$^2$. The strong Na\&Ca absorption suggests that the absorbing gas is thick and mostly neutral, which is also supported by dust extinction $A_V\sim3$. Using the best-estimated hydrogen column density $N_H\sim3\times10^{22}$ cm$^{-2}$, the total mass of the outflowing gas is $>3\times10^9 M_\odot$. The outflow is likely to be driven by AGN because of the $\sim$800 km/s blueshift velocity, suggesting SDSS J072910.34+333634.3 is undergoing one of the most violent AGN feedback we have seen. In the future, one can find more such massive neutral gas outflows in other reddened quasars using similar method, and this can shed new light on the study of AGN feedback.

The sporadic accretion following the tidal disruption of a star by a super-massive black hole (TDE) leads to a bright UV and soft X-ray flare in the galactic nucleus. The gas and dust surrounding the black hole responses to such a flare with an echo in emission lines and infrared emission. In this paper, we report the detection of long fading mid-IR emission lasting up to 14 years after the flare in four TDE candidates with transient coronal lines using the WISE public data release. We estimate that the reprocessed mid-IR luminosities are in the range between $4\times 10^{42}$ and $2\times 10^{43}$ erg~s$^{-1}$ and dust temperature in the range of 570-800K when WISE first detected these sources three to five years after the flare. Both luminosity and dust temperature decreases with time. We interpret the mid-IR emission as the infrared echo of the tidal disruption flare. We estimate the UV luminosity at the peak flare to be 1 to 30 times $10^{44}$ erg s$^{-1}$ and for warm dust masses to be in the range of 0.05-1.3 Msun within a few parsecs. Our results suggest that the mid-infrared echo is a general signature of TDE in the gas-rich environment.

Using low-redshift (z<0.09) samples of AGNs, normal galaxies and groups of galaxies selected from the Sloan Digital Sky Survey (SDSS), we study the environments of type 1 and type 2 AGNs both on small and large scales. Comparisons are made for galaxy samples matched in redshift, $r$-band luminosity, [OIII] luminosity, and also the position in groups (central or satellite). We find that type 2 AGNs and normal galaxies reside in similar environments. Type 1 and type 2 AGNs have similar clustering properties on large scales ($\gtrsim1$Mpc), but at scales smaller than 100 kpc, type 2s have significant more neighbors than type 1s ($3.09\pm0.69$ times more for central AGNs at $\lesssim30$kpc). These results suggest that type 1 and type 2 AGNs are hosted by halos of similar masses, as is also seen directly from the mass distributions of their host groups ($\sim10^{12}h^{-1} M_{\odot}$ for centrals and $\sim10^{13}h^{-1} M_{\odot}$ for satellites). Type~2s have significantly more satellites around them, and the distribution of their satellites is also more centrally concentrated. The host galaxies of both types of AGNs have similar optical properties, but their infrared colors are significantly different. Our results suggest that the simple unified model based solely on torus orientation is not sufficient, but that galaxy interactions in dark matter halos must have played an important role in the formation of the dust structure that obscures AGNs.

We report the detection of a significant infrared variability of the nearest tidal disruption event (TDE) ASASSN-14li using Wide-field Infrared Survey Explorer} and newly released Near-Earth Object WISE Reactivation} data. In comparison with the quiescent state, the infrared flux is brightened by 0.12 and 0.16 magnitude in the W1 ($3.4\mu$m) and W2 ($4.6\mu$m) bands at 36 days after the optical discovery (or $\sim110$ days after the peak disruption date). The flux excess is still detectable $\sim170$ more days later. Assuming that the flare-like infrared emission is from the dust around the black hole, its blackbody temperature is estimated to be $\sim2.1\times10^3$~K, slightly higher than the dust sublimation temperature, indicating that the dust is likely located close to the dust sublimation radius. The equilibrium between the heating and radiation of the dust claims a bolometric luminosity of $\sim10^{43}-10^{45}$~\lum, comparable with the observed peak luminosity. This result has for the first time confirmed the detection of infrared emission from the dust echoes of TDEs.

The accretion of interstellar medium onto the central super massive black holes is widely accepted as the source of the gigantic energy released by the active galactic nuclei. But few pieces of observational evidence have been confirmed directly demonstrating the existence of the inflows. The absorption line system in the spectra of quasar SDSS J112526.12+002901.3 presents an interesting example, in which the rarely detected hydrogen Balmer and metastable He I absorption lines are found redshifted to the quasar's rest frame along with the low-ionization metal absorption lines Mg II, Fe II, etc. The repeated SDSS spectroscopic observations suggest a transverse velocity smaller than the radial velocity. The motion of the absorbing medium is thus dominated by infall. The He I* lines present a powerful probe to the strength of ionizing flux, while the Balmer lines imply a dense environment. With the help of photoionization simulations, we find the absorbing medium is exposed to the radiation with ionization parameter $U\approx 10^{-1.8}$, and the density is $n(\mathrm{H})\approx 10^9\ \mathrm{cm}^{-3}$. Thus the absorbing medium is located $\sim 4\ \mathrm{pc}$ away from the central engine. According to the similarity in the distance and physical conditions between the absorbing medium and the torus, we strongly propose the absorption line system as a candidate for the accretion inflow which originates from the inner surface of the torus.

Emission lines from the broad emission line region (BELR) and the narrow emission line region (NELR) of active galactic nuclei (AGNs) are extensively studied. However, between these two regions emission lines are rarely detected. We present a detailed analysis of a quasar SDSS J232444.80-094600.3 (SDSS J2324$-$0946), which is remarkable for its strong intermediate-width emission lines (IELs) with FWHM $\approx$ 1800 \kmps. The IEL component is presented in different emission lines, including the permitted lines \lya\ $\lambda$1216, \civ\ $\lambda$1549, semiforbidden line \ciii\ $\lambda$1909, and forbidden lines \oiii\ $\lambda\lambda$4959, 5007. With the aid of photo-ionization models, we found that the IELs are produced by gas with a hydrogen density of $n_{\rm H} \sim 10^{6.2}-10^{6.3}~\rm cm^{-3}$, a distance to the central ionizing source of $R \sim 35-50$ pc, a covering factor of CF $\sim$ 6\%, and a dust-to-gas ratio of $\leq 4\%$ times of SMC. We suggest that the strong IELs of this quasar are produced by nearly dust-free and intermediate-density gas located at the skin of the dusty torus. Such strong IELs, served as a useful diagnose, can provide an avenue to study the properties of gas between the BELR and the NELR.

SDSS J1634+2049 is a local (z = 0.1293) infrared-luminous quasar with LIR= 10^11.91 Lsun. We present a detailed multiwavelength study of both the host galaxy and the nucleus. The host galaxy demonstrates violent, obscured star formation activities with SFR ~ 140 Msun yr^-1, estimated from either the PAH emission or IR luminosity. The optical to NIR spectra exhibit a blueshifted narrow cuspy component in Hb, HeI5876,10830 and other emission lines consistently with an offset velocity of ~900 km/s, as well as additional blueshifting phenomena in high-ionization lines , while there exist blueshifted broad absorption lines (BALs) in NaID and HeI*3889,10830, indicative of the AGN outflows producing BALs and emission lines. Constrained mutually by the several BALs with CLOUDY, the physical properties of the absorption-line outflow are derived as follows: 10^4 < n_H <= 10^5 cm^-3, 10^-1.3 <= U <= 10^-0.7 and 10^22.5<= N_H <= 10^22.9 cm^-2 , similar to those derived for the emission-line outflows. The similarity suggests a common origin. Taking advantages of both the absorption lines and outflowing emission lines, we find that the outflow gas is located at a distance of 48 - 65 pc from the nucleus, and that the kinetic luminosity of the outflow is 10^44-10^46 erg s^-1. J1634+2049 has a off-centered galactic ring on the scale of ~ 30 kpc that is proved to be formed by a recent head-on collision by a nearby galaxy. Thus this quasar is a valuable object in the transitional phase emerging out of dust enshrouding as depicted by the co-evolution scenario.

SDSS J095253.83+011422.0 (SDSS J0952+0114) was reported by Hall et al. (2004) as an exotic quasar at $z_{\rm em}=3.020$. In contrast to prominent broad metal--line emissions with FWHM~9000 km/s, only a narrow Ly \alpha emission line is present with FWHM~1000 km/s. The absence of broad Ly alpha emission line has been a mystery for more than a decade. In this paper, we demonstrate that this is due to dark Proximate Damped Ly alpha Absorption (PDLA) at $z_{\rm abs}=3.010$ by identifying associated Lyman absorption line series from the damped Ly beta up to Ly9, as well as the Lyman limit absorption edge. The PDLA cloud has a column density of $\log N_{\rm H\,I}({\rm cm}^{-2})=21.8\pm0.2$, a metallicity of [Zn/H]$>-1.0$, and a spatial extent exceeding the Narrow Emission Line Region (NELR) of the quasar. With a luminosity of $L_{{\rm Ly}\alpha}\sim10^{45}$ erg s$^{-1}$, the residual Ly alpha emission superposed on the PDLA trough is of two orders of magnitude stronger than previous reports. This is best explained as re-radiated photons arising from the quasar outflowing gas at a scale larger than the NELR. The PDLA here, acting like a natural coronagraph, provides us with a good insight into the illuminated gas in the vicinity of the quasar, which are usually hard to resolve due to their small size and "seeing fuzz" of bright quasars. Notably, SDSS J0952+0114 analogs might be easily omitted in the spectroscopic surveys of DLAs and PDLAs, as their damped Ly alpha troughs can be fully filled by additional strong Ly alpha emissions. Our preliminary survey shows that such systems are not very rare. They are potentially a unique sample for probing strong quasar feedback phenomena in the early universe.

We present an analysis of spectrum and variability of the bright reddened narrow line Seyfert 1 galaxy Was~61 using 90 ks archival {\it XMM-Newton} data. The X-ray spectrum in 0.2-10 keV can be characterized by an absorbed power-law plus soft excess and an Fe K$\alpha$ emission line. The power-law spectral index remains constant during the flux variation. The absorbing material is mildly ionized, with a column density of 3.2$\times$10$^{21}$ cm$^{-2}$, and does not appear to vary during the period of the X-ray observation. If the same material causes the optical reddening (E(B-V)$\simeq$0.6 mag), it must be located outside the narrow line region with a dust-to-gas ratio similar to the average Galactic value. We detect significant variations of the Fe K$\alpha$ line during the observational period. A broad Fe K$\alpha$ line at $\simeq$6.7 keV with a width of $\sim$0.6 keV is detected in the low flux segment of the first 40 ks exposure, and is absent in the spectra of other segments; a narrow Fe K$\alpha$ emission line $\sim$6.4 keV with a width of $\sim$0.1 keV is observed in the subsequent 20 ks segment, which has a count rate of 35% higher and is in the next day. We believe this is due to the change in geometry and kinematics of the X-ray emitting corona. The temperature and flux of soft X-ray excess appear to correlate with the flux of the hard power-law component. Comptonization of disc photons by a warm and optically thick inner disk is preferred to interpret the soft excess, rather than the ionized reflection.

We report on the discovery of broad Balmer absorption lines variability in the QSO SDSS J125942.80+121312.6, based on the optical and near-infrared spectra taken from the SDSS-I, SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS), and TripleSpec observations over a timescale of 5.8 years in the QSO's rest-frame. The blueshifted absorption profile of H$\beta$ shows a variation of more than 5$\sigma$ at a high velocity portion ($>3000\ \mathrm{km\ s}^{-1}$) of the trough. We perform a detailed analysis for the physical conditions of the absorber using Balmer lines as well as metastable He I and optical Fe II absorptions ($\lambda 4233$ from b$^4$P$_{5/2}$ level and $\lambda 5169$ from a$^6$S$_{5/2}$) at the same velocity. These Fe II lines are identified in the QSO spectra for the first time. According to the photoionization simulations, we estimate a gas density of $n(\mathrm{H})\approx 10^{9.1}\ \mathrm{cm}^{-3}$ and a column density of $N_{\mathrm{col}}(\mathrm{H})\approx 10^{23}\ \mathrm{cm}^{-2}$ for the BOSS data, but the model fails to predict the variations of ionic column densities between the SDSS and BOSS observations if changes in ionizing flux are assumed. We thus propose transverse motion of the absorbing gas being the cause of the observed broad Balmer absorption line variability. In fact, we find that the changes in covering factors of the absorber can well-reproduce all of the observed variations. The absorber is estimated $\sim 0.94$ pc away from the central engine, which is where the outflow likely experiences deceleration due to the collision with the surrounding medium. This scheme is consistent with the argument that LoBAL QSOs may represent the transition from obscured star-forming galaxies to classic QSOs.

Quasi-periodic oscillation (QPO) detected in the X-ray radiation of black hole X-ray binaries (BHXBs) is thought to originate from dynamical processes in the close vicinity of the black holes (BHs), and thus carries important physical information therein. Such a feature is extremely rare in active galactic nuclei (AGNs) with supermassive BHs. Here we report on the detection of a possible X-ray QPO signal with a period of 3800\,s at a confidence level $>99.99\%$ in the narrow-line Seyfert 1 galaxy (NLS1) 1H~0707-495 in one data set in 0.2-10\,keV taken with {\it XMM-Newton}. The statistical significance is higher than that of most previously reported QPOs in AGNs. The QPO is highly coherent (quality factor $Q=\nu/\Delta\nu \geqslant 15$) with a high rms fractional variability ($\sim15\%$). A comprehensive analysis of the optical spectra of this AGN is also performed, yielding a central BH mass $5.2\times10^6\,M_{\odot}$ from the broad emission lines based on the scaling relation. The QPO follows closely the known frequency-BH mass relation, which spans from stellar-mass to supermassive BHs. The absence of the QPO in other observations of the object suggests it a transient phenomenon. We suggest that the (high-frequency) QPOs tend to occur in highly accreting BH systems, from BHXBs to supermassive BHs. Future precise estimation of the BH mass may be used to infer the BH spin from the QPO frequency.

We present Keck/ESI long-slit spectroscopy of SBS 1421+511, a system consisting of a quasar at z = 0.276 and an extended source 3" northern to the quasar. The quasar shows a blue-skewed profile of Balmer broad emission lines, which can be well modeled as emissions from a circular disk with a blueshift velocity of ~1400 km/s. The blueshift is better interpreted as resulting from a recoiling active black hole than from a super-massive black hole binary, since the line profile almost kept steady over one decade in the quasar rest-frame. Alternative interpretations are possible as well, such as emissions from a bipolar outflow or a circular disk with spiral emissivity perturbations. The extended source shows Seyfert-like narrow line ratios and a [OIII] luminosity of >1.4\times10^8L_\odot, with almost the same redshift as the quasar and a projected distance of 12.5 kpc at the redshift. SBS 1421+511 is thus likely to be an interacting galaxy pair with dual AGN. Alternatively, the quasar companion only appears to be active but not necessarily so: the gas before/in/behind the companion galaxy is illuminated by the quasar as an extended emission line region is detected at a similar distance in the opposite direction southern to the quasar, which may be generated either by tidal interactions between the galaxy pair or large-scale outflows from the quasar.

The existence of intermediate-width emission line regions (IELRs) in active galactic nuclei has been discussed for over two decades. A consensus, however, is yet to be arrived at due to the lack of convincing evidence for their detection. We present a detailed analysis of the broadband spectrophotometry of the partially obscured quasar OI 287. The ultraviolet intermediate-width emission lines (IELs) are very prominent, in high contrast to the corresponding broad emission lines (BELs) which are heavily suppressed by dust reddening. Assuming that the IELR is virialized, we estimated its distance to the central black hole of $\sim 2.9$ pc, similar to the dust sublimation radius of $\sim 1.3$ pc. Photo-ionization calculations suggest that the IELR has a hydrogen density of $\sim 10^{8.8}-10^{9.4} ~ \rm cm^{-3}$, within the range of values quoted for the dusty torus near the sublimation radius. Both its inferred location and physical conditions suggest that the IELR originates from the inner surface of the dusty torus. In the spectrum of this quasar, we identified only one narrow absorption-line system associated with the dusty material. With the aid of photo-ionization model calculations, we found that the obscuring material might originate from an outer region of the dusty torus. We speculate that the dusty torus, which is exposed to the central ionizing source, may produce IELs through photo-ionization processes, while also obscure BELs as a natural "coronagraph". Such a "coronagraph" could be found in a large number of partially obscured quasars and be a useful tool to study IELRs.

We apply a method recently introduced to the statistical literature to directly estimate the precision matrix from an ensemble of samples drawn from a corresponding Gaussian distribution. Motivated by the observation that cosmological precision matrices are often approximately sparse, the method allows one to exploit this sparsity of the precision matrix to more quickly converge to an asymptotic 1/sqrt(Nsim) rate while simultaneously providing an error model for all of the terms. Such an estimate can be used as the starting point for further regularization efforts which can improve upon the 1/sqrt(Nsim) limit above, and incorporating such additional steps is straightforward within this framework. We demonstrate the technique with toy models and with an example motivated by large-scale structure two-point analysis, showing significant improvements in the rate of convergence.For the large-scale structure example we find errors on the precision matrix which are factors of 5 smaller than for the sample precision matrix for thousands of simulations or, alternatively, convergence to the same error level with more than an order of magnitude fewer simulations.

We present the discovery of Balmer line absorption from H$\alpha$ to H$\gamma$ in an iron low-ionization broad absorption line (FeLoBAL) quasar SDSS J152350.42+391405.2 (hereafter J1523), by the quasi-simultaneous optical and near-infrared spectroscopy. The Balmer line absorption is at $z_{absor}$ = 0.6039 +/-0.0021 and blueshifted by v=10,353 km/s with respect to the Balmer emission lines. All Balmer BALs have uniform absorption profile with the widths of $\Delta$ v ~12,000 km/s. We also found the absorption trough in He 1* $\lambda$10830 with the same velocity and width in the H-band TripleSpec spectrum of J1523. This object is only the tenth active galactic nucleus known to exhibit non-stellar Balmer absorption, and also the case with the highest velocity and broadest Balmer absorption lines which have ever been found. A CLOUDY analysis shows that the absorbers require an gas density of $log_{10} n_ e (cm^{-3})=9$ and an ionization parameter of $log_{10} U=-1.0$. They locate at a distance of ~0.2 pc from the central ionizing source which is slightly farther than that of BELRs. Furthermore, J1523 is one of the brightest Balmer-BAL quasar ever reported, with unique iron absorption variations, making it as the most promising candidate for follow up high-resolution spectroscopy, multi-band observations, and long-term monitoring.

We compared optical spectroscopic and photometric data for 18 AGN galaxies over 2 to 3 epochs, with time intervals of typically 5 to 10 years. We used the Multi-Object Double Spectrograph (MODS) at the Large Binocular Telescope (LBT) and compared the spectra to data taken from the SDSS database and the literature. We find variations in the forbidden oxygen lines as well as in the hydrogen recombination lines of these sources. For 4 of the sources we find that, within the calibration uncertainties, the variations in continuum and line spectra of the sources are very small. We argue that it is mainly the difference in black hole mass between the samples that is responsible for the different degree of continuum variability. In addition we find that for an otherwise constant accretion rate the total line variability (dominated by the narrow line contributions) reverberates the continuum variability with a dependency $\Delta L_{line} \propto (\Delta L_{cont.})^{\frac{3}{2}}$. Since this dependency is prominently expressed in the narrow line emission it implies that the luminosity dominating part of the narrow line region must be very compact with a size of the order of at least 10 light years. A comparison to literature data shows that these findings describe the variability characteristics of a total of 61 broad and narrow line sources.

List of proceedings from the HAWC Collaboration presented at the 34th International Cosmic Ray Conference, 30 July - 6 August 2015, The Hague, The Netherlands.

A survey of the inner Galaxy region of Galactic longitude l in [+15, +50] degree and latitude b in [-4,+4] degree is performed using one-third of the High Altitude Water Cherenkov (HAWC) Observatory operated during its construction phase. To address the ambiguities arising from unresolved sources in the data, we use a maximum likelihood technique to identify point source candidates. Ten sources and candidate sources are identified in this analysis. Eight of these are associated with known TeV sources but not all have differential fluxes compatible with previous measurements. Three sources are detected with significances $>5\,\sigma$ after accounting for statistical trials, and are associated with known TeV sources.

We report on the identification of a new $\gamma$-ray-emitting narrow-line Seyfert 1 (NLS1) galaxy, SDSS J122222.55+041315.7, which increases the number of known objects of this remarkable but rare type of active galactic nuclei (AGN) to seven. Its optical spectrum, obtained in the Sloan Digital Sky Survey-Baryon Oscillation Spectroscopic Survey, reveals a broad H $\beta$ emission line with a width (FWHM) of 1734$\pm$104 km s$^{-1}$. This, along with strong optical Fe II multiplets [$R_{4570}=0.9$] and a weak [O III] $\lambda 5007$ emission line, makes the object a typical NLS1. On the other hand, the source exhibits a high radio brightness temperature, rapid infrared variability, and a flat X-ray spectrum extending up to $\sim$200 keV. It is associated with a luminous $\gamma$-ray source detected significantly with {\it Fermi}/LAT. Correlated variability with other wavebands has not yet been tested. The spectral energy distribution can be well modelled by a one-zone leptonic jet model. This new member is by far the most distant $\gamma$-ray-emitting NLS1, at a redshift of $z=0.966$.

A number of Galactic sources emit GeV-TeV gamma rays that are produced through leptonic and/or hadronic mechanisms. Spectral analysis in this energy range is crucial in order to understand the emission mechanisms. The HAWC Gamma-Ray Observatory, with a large field of view and location at $19^\circ$ N latitude, is surveying the Galactic Plane from high Galactic longitudes down to near the Galactic Center. Data taken with partially-constructed HAWC array in 2013-2014 exhibit TeV gamma-ray emission along the Galactic Plane. A high-level analysis likelihood framework for HAWC, also presented at this meeting, has been developed concurrently with the Multi-Mission Maximum Likelihood (3ML) architecture to deconvolve the Galactic sources and to perform multi-instrument analysis. It has been tested on early HAWC data and the same method will be applied on HAWC data with the full array. I will present preliminary results on Galactic sources from TeV observations with HAWC and from joint analysis on Fermi and HAWC data in GeV-TeV energy range.

We present the detection of compact radio structures of fourteen radio-loud narrow line Seyfert 1 (NLS1) galaxies from Very Long Baseline Array observations at 5 GHz, which were performed in 2013. While 50\% of the sources of our sample show a compact core only, the remaining 50\% exhibit a core-jet structure. The measured brightness temperatures of the cores range from $10^{8.4}$ to $10^{11.4}$ K with a median value of $10^{10.1}$ K, indicating that the radio emission is from non-thermal jets, and that, likely, most sources are not strongly beamed, then implying a lower jet speed in these radio-loud NLS1 galaxies. In combination with archival data taken at multiple frequencies, we find that seven sources show flat or even inverted radio spectra, while steep spectra are revealed in the remaining seven objects. Although all these sources are very radio-loud with $R > 100$, their jet properties are diverse, in terms of their milli-arcsecond (mas) scale (pc scale) morphology and their overall radio spectral shape. The evidence for slow jet speeds (i.e., less relativistic jets), in combination with the low kinetic/radio power, may offer an explanation for the compact VLBA radio structure in most sources. The mildly relativistic jets in these high accretion rate systems are consistent with a scenario, where jets are accelerated from the hot corona above the disk by the magnetic field and the radiation force of the accretion disk. Alternatively, a low jet bulk velocity can be explained by low spin in the Blandford-Znajek mechanism.

The High Altitude Water Cherenkov (HAWC) Observatory continuously observes gamma-rays between 100 GeV to 100 TeV in an instantaneous field of view of about 2 steradians above the array. The large amount of raw data, the importance of small number statistics, the large dynamic range of gamma-ray signals in time (1 - $10^8$ sec) and angular extent (0.1 - 100 degrees), and the growing need to directly compare results from different observatories pose some special challenges for the analysis of HAWC data. To address these needs, we have designed and implemented a modular analysis framework based on the method of maximum likelihood. The framework facilitates the calculation of a binned Poisson Log-likelihood value for a given physics model (i.e., source model), data set, and detector response. The parameters of the physics model (sky position, spectrum, angular extent, etc.) can be optimized through a likelihood maximization routine to obtain a best match to the data. In a similar way, the parameters of the detector response (absolute pointing, angular resolution, etc.) can be optimized using a well-known source such as the Crab Nebula. The framework was designed concurrently with the Multi-Mission Maximum Likelihood (3ML) architecture, and allows for the definition of a general collection of sources with individually varying spectral and spatial morphologies. Compatibility with the 3ML architecture allows to easily perform powerful joint fits with other observatories. In this contribution, we overview the design and capabilities of the HAWC analysis framework, stressing the overarching design points that have applicability to other astronomical and cosmic-ray observatories.

The majority of Galactic TeV gamma-ray sources are pulsar wind nebulae (PWNe) and supernova remnants (SNRs), and the most common association for unidentified sources is PWN. Many of these sources were discovered in TeV by imaging air Cherenkov telescopes using overlapping pointed observations over sections of the Galactic plane. The HAWC observatory is a survey type instrument in the Northern hemisphere with an energy range of 100 GeV to 100 TeV. Preliminary analysis of data recorded with the partially completed HAWC array taken since 2013 shows extended detections that are coincident with known TeV SNRs and PWNe. The full array became operational in early 2015 and has been steadily surveying the Northern sky since. I will discuss detections in HAWC data taken since 2013 associated with PWNe and SNRs.

The HAWC collaboration has recently completed the construction of a gamma-ray observatory at an altitude of 4100 meters on the slope of the Sierra Negra volcano in the state of Puebla, Mexico. In order to achieve an optimal angular resolution, energy reconstruction, and cosmic-ray background suppression for the air showers observed by HAWC, it is crucial to obtain good timing and charge calibrations of the photosensors in the detector. The HAWC calibration is based on a laser system which is able to deliver short light pulses to all the tanks in the array. The light intensity can range over 7 orders of magnitude, broad enough to cover all the dynamic range of the PMT readout electronics. In this contribution we will present the HAWC calibration system, together with the methods used to calibrate the detector.

Astrophysical sources are now observed by many different instruments at different wavelengths, from radio to high-energy gamma-rays, with an unprecedented quality. Putting all these data together to form a coherent view, however, is a very difficult task. Each instrument has its own data format, software and analysis procedure, which are difficult to combine. It is for example very challenging to perform a broadband fit of the energy spectrum of the source. The Multi-Mission Maximum Likelihood framework (3ML) aims to solve this issue, providing a common framework which allows for a coherent modeling of sources using all the available data, independent of their origin. At the same time, thanks to its architecture based on plug-ins, 3ML uses the existing official software of each instrument for the corresponding data in a way which is transparent to the user. 3ML is based on the likelihood formalism, in which a model summarizing our knowledge about a particular region of the sky is convolved with the instrument response and compared to the corresponding data. The user can choose between a frequentist analysis, and a Bayesian analysis. In the former, parameters of the model are optimized in order to obtain the best match to the data (i.e., the maximum of the likelihood). In the latter, the priors specified by the user are used to build the posterior distribution, which is then sampled with Markov Chain Monte Carlo or Multinest. Our implementation of this idea is very flexible, allowing the study of point sources as well as extended sources with arbitrary spectra. We will review the problem we aim to solve, the 3ML concepts and its innovative potential.

We present the variability study of broad absorption lines (BALs) in a uniformly radio-selected sample of 28 BAL quasars using the archival data from the FIRST Bright Quasar Survey (FBQS) and the Sloan Digital Sky Survey (SDSS), as well as those obtained by ourselves, covering time scales $\sim 1-10$ years in the quasar's rest-frame. To our surprise, 5 quasars showing strong variations are all belong to a special subclass of overlapping iron low ionization BAL (OFeLoBAL) quasars, however, other 4 non-overlapping FeLoBALs (non-OFeLoBALs) are invariable except one case with weak optical depth change. Meanwhile, we also find 6 typical variations of high-ionization and low-ionization BALs in this BAL quasar sample. Photoionization models suggest that OFeLoBALs are formed in a relative dense ($n_e>10^6$ cm$^{-3}$) outflows at a distance from the subparsec to the dozens of parsecs from the continuum source. They differ from those of non-OFeLoBALs, which are likely produced by low-density gas, locating at a distance of hundreds to thousands parsecs. Thus, OFeLoBALs and non-OFeLoBALs, i.e., FeLoBALs with/without strong BAL variations, perhaps represent the bimodality of Fe II absorption, the former is located in the active galactic nucleus environment rather than the host galaxy. We suggest that high density and small distance are the necessary conditions that cause OFeLoBALs. As suggested in the literature, strong BAL variability is possibly due to variability of the covering factor of BAL regions caused by clouds transiting across the line of sight rather than ionization variations.

Neutral Helium multiplets, HeI*3189,3889,10830 are very useful diagnostics to the geometry and physical conditions of the absorbing gas in quasars. So far only a handful of HeI* detections have been reported. Using a newly developed method, we detected HeI*3889 absorption line in 101 sources of a well-defined sample of 285 MgII BAL quasars selected from the SDSS DR5. This has increased the number of HeI* BAL quasars by more than one order of magnitude. We further detected HeI*3189 in 50% (52/101) quasars in the sample. The detection fraction of HeI* BALs in MgII BAL quasars is about 35% as a whole, and increases dramatically with increasing spectral signal-to-noise ratios, from 18% at S/N <= 10 to 93% at S/N >= 35. This suggests that HeI* BALs could be detected in most MgII LoBAL quasars, provided spectra S/N is high enough. Such a surprisingly high HeI* BAL fraction is actually predicted from photo-ionization calculations based on a simple BAL model. The result indicates that HeI* absorption lines can be used to search for BAL quasars at low-z, which cannot be identified by ground-based optical spectroscopic survey with commonly seen UV absorption lines. Using HeI*3889, we discovered 19 BAL quasars at z<0.3 from available SDSS spectral database. The fraction of HeI* BAL quasars is similar to that of LoBAL objects.

We report the discovery of excess broad band absorption near 2250 A (EBBA) in the spectra of seven broad absorption line (BAL) quasars. By comparing with the statistical results from the control quasar sample, the significance for the detections are all above the > 4{\sigma} level, with five above > 5{\sigma}. The detections have also been verified by several other independent methods. The EBBAs present broader and weaker bumps at smaller wavenumbers than the Milky Way, and similar to the Large Magellanic Cloud. The EBBA bump may be related to the 2175 A bump seen in the Local Group and may be a counterpart of the 2175 A bump under different conditions in the early Universe. Furthermore, five objects in this sample show low-ionization broad absorption lines (LoBALs), such as Mg II and Al III, in addition to the high-ionization broad absorption lines (HiBALs) of C IV and Si IV. The fraction of LoBALs in our sample, ~70%, is surprisingly high compared to that of general BAL quasars, ~10%. Although the origin of the bump is still not clear, the coexistence of both BALs and bumps and the significantly high fraction of LoBALs may indicate the bump carriers is closely related to the early evolution phase of quasars.

The High Altitude Water Cherenkov (HAWC) observatory is a wide field-of-view detector sensitive to gamma rays of 100 GeV to a few hundred TeV. Located in central Mexico at 19 degrees North latitude and 4100 m above sea level, HAWC will observe gamma rays and cosmic rays with an array of water Cherenkov detectors. The full HAWC array is scheduled to be operational in Spring 2015. In this paper, we study the HAWC sensitivity to the gamma-ray signatures of high-mass (multi- TeV) dark matter annihilation. The HAWC observatory will be sensitive to diverse searches for dark matter annihilation, including annihilation from extended dark matter sources, the diffuse gamma-ray emission from dark matter annihilation, and gamma-ray emission from non-luminous dark matter subhalos. Here we consider the HAWC sensitivity to a subset of these sources, including dwarf galaxies, the M31 galaxy, the Virgo cluster, and the Galactic center. We simulate the HAWC response to gamma rays from these sources in several well-motivated dark matter annihilation channels. If no gamma-ray excess is observed, we show the limits HAWC can place on the dark matter cross-section from these sources. In particular, in the case of dark matter annihilation into gauge bosons, HAWC will be able to detect a narrow range of dark matter masses to cross-sections below thermal. HAWC should also be sensitive to non-thermal cross-sections for masses up to nearly 1000 TeV. The constraints placed by HAWC on the dark matter cross-section from known sources should be competitive with current limits in the mass range where HAWC has similar sensitivity. HAWC can additionally explore higher dark matter masses than are currently constrained.

We report the identification of an unusual absorption line system in the quasar SDSS J080248.18$+$551328.9 and present a detailed study of the system, incorporating follow-up optical and NIR spectroscopy. A few tens of absorption lines are detected, including He I*, Fe II* and Ni II* that arise from metastable or excited levels, as well as resonant lines in Mg I, Mg II, Fe II, Mn II, and Ca II. All of the isolated absorption lines show the same profile of width $\Delta v\sim 1,500$km s$^{-1}$ centered at a common redshift as that of the quasar emission lines, such as [O II], [S II], and hydrogen Paschen and Balmer series. With narrow Balmer lines, strong optical Fe II multiplets, and weak [O III] doublets, its emission line spectrum is typical for that of a narrow-line Seyfert 1 galaxy (NLS1). We have derived reliable measurements of the gas-phase column densities of the absorbing ions/levels. Photoionization modeling indicates that the absorber has a density of $n_{\rm H} \sim (1.0-2.5)\times 10^5~ {\rm cm}^{-3}$ and a column density of $N_{\rm H} \sim (1.0-3.2)\times 10^{21} \sim {\rm cm}^{-2}$, and is located at $R\sim100-250$ pc from the central super-massive black hole. The location of the absorber, the symmetric profile of the absorption lines, and the coincidence of the absorption and emission line centroid jointly suggest that the absorption gas is originated from the host galaxy and is plausibly accelerated by stellar processes, such as stellar winds \zhy{and/or} supernova explosions. The implications for the detection of such a peculiar absorption line system in an NLS1 are discussed in the context of co-evolution between super-massive black hole growth and host galaxy build-up.

We report the finding of an unobscured type II Active Galactic Nuclei (AGN) candidate, SDSS J012032.19-005501.9 at a relatively high redshift of 0.601,which shows a number of unusual properties. It varies significantly on timescales of years as typical type I AGNs and marginally on timescales of weeks. The color-magnitude relation and the structure function are also consistent with that of type I AGNs, which imply that its variability likely originates from the black hole accretion system .However, no broad emission line is detected in the SDSS spectrum, and the upper limit of the equivalent width of the H$\rm \beta$ broad emission line is much less than that of type I AGNs. These properties suggest that SDSS J012032.19-005501.9 may be an unobscured quasar without broad emission lines intrinsically, namely an unobscured type II AGN or "true" type II AGN. Furthermore, its continuum luminosity is at least one order of magnitude fainter than the average value of thepast century expected from the [OIII] emission line. It indicates that SDSS J012032.19-005501.9 may be switching off. Additional possible scenarios to explain this intriguing source are also discussed. Future deep observations at multi-wavelengths are needed to reveal the nature of this peculiar and intriguing AGN.

We derive photometric redshifts (\zp) for sources in the entire ($\sim0.4$ deg$^2$) Hawaii-Hubble Deep Field-North (\hdfn) field with the EAzY code, based on point spread function-matched photometry of 15 broad bands from the ultraviolet (\bandu~band) to mid-infrared (IRAC 4.5 $\mu$m). Our catalog consists of a total of 131,678 sources. We evaluate the \zp~quality by comparing \zp~with spectroscopic redshifts (\zs) when available, and find a value of normalized median absolute deviation \sigm$=$0.029 and an outlier fraction of 5.5\% (outliers are defined as sources having $\rm |\zp - \zs|/(1+\zs) > 0.15$) for non-X-ray sources. More specifically, we obtain \sigm$=0.024$ with 2.7\% outliers for sources brighter than $R=23$~mag, \sigm$=0.035$ with 7.4\% outliers for sources fainter than $R=23$~mag, \sigm$=$0.026 with 3.9\% outliers for sources having $z<1$, and \sigm$=$0.034 with 9.0\% outliers for sources having $z>1$. Our \zp\ quality shows an overall improvement over an earlier \zp\ work that focused only on the central \hdfn\ area. We also classify each object as star or galaxy through template spectral energy distribution fitting and complementary morphological parametrization, resulting in 4959 stars and 126,719 galaxies. Furthermore, we match our catalog with the 2~Ms {\it Chandra} Deep Field-North main \xray~catalog. For the 462 matched non-stellar \xray~sources (281 having \zs), we improve their \zp~quality by adding three additional AGN templates, achieving \sigm$=0.035$ and an outlier fraction of 12.5\%. We make our catalog publicly available presenting both photometry and \zp, and provide guidance on how to make use of our catalog.

The High-Altitude Water Cherenkov (HAWC) Observatory is sensitive to gamma rays and charged cosmic rays at TeV energies. The detector is still under construction, but data acquisition with the partially deployed detector started in 2013. An analysis of the cosmic-ray arrival direction distribution based on $4.9\times 10^{10}$ events recorded between June 2013 and February 2014 shows anisotropy at the $10^{-4}$ level on angular scales of about $10^\circ$. The HAWC cosmic-ray sky map exhibits three regions of significantly enhanced cosmic-ray flux; two of these regions were first reported by the Milagro experiment. A third region coincides with an excess recently reported by the ARGO-YBJ experiment. An angular power spectrum analysis of the sky shows that all terms up to $\ell=15$ contribute significantly to the excesses.

Primordial Black Holes (PBHs) are gravitationally collapsed objects that may have been created by density fluctuations in the early universe and could have arbitrarily small masses down to the Planck scale. Hawking showed that due to quantum effects, a black hole has a temperature inversely proportional to its mass and will emit all species of fundamental particles thermally. PBHs with initial masses of ~5.0 x 10^14 g should be expiring in the present epoch with bursts of high-energy particles, including gamma radiation in the GeV - TeV energy range. The Milagro high energy observatory, which operated from 2000 to 2008, is sensitive to the high end of the PBH evaporation gamma-ray spectrum. Due to its large field-of-view, more than 90% duty cycle and sensitivity up to 100 TeV gamma rays, the Milagro observatory is well suited to perform a search for PBH bursts. Based on a search on the Milagro data, we report new PBH burst rate density upper limits over a range of PBH observation times. In addition, we report the sensitivity of the Milagro successor, the High Altitude Water Cherenkov (HAWC) observatory, to PBH evaporation events.

VAMOS was a prototype detector built in 2011 at an altitude of 4100m a.s.l. in the state of Puebla, Mexico. The aim of VAMOS was to finalize the design, construction techniques and data acquisition system of the HAWC observatory. HAWC is an air-shower array currently under construction at the same site of VAMOS with the purpose to study the TeV sky. The VAMOS setup included six water Cherenkov detectors and two different data acquisition systems. It was in operation between October 2011 and May 2012 with an average live time of 30%. Besides the scientific verification purposes, the eight months of data were used to obtain the results presented in this paper: the detector response to the Forbush decrease of March 2012, and the analysis of possible emission, at energies above 30 GeV, for long gamma-ray bursts GRB111016B and GRB120328B.

We compute the intrinsic isotropic peak luminosity function (LF) and formation rate of long gamma-ray bursts (LGRBs) using a novel approach. We complement a standard log\,$N$\,--\,log\,$P$ brightness distribution and $V_{\mathrm{max}}$ estimations with two observation-time relations: a redshift--observation-time relation (log\,$z$\,--\,log\,$T$) and a new luminosity--observation-time relation (log\,$L$\,--\,log\,$T$). We show that this approach reduces degeneracies that exist between the rate and LF of a brightness distribution. To account for the complex triggering algorithm employed by \emph{Swift} we use recent results of \citet{Lien_2014ApJ} to produce a suite of efficiency functions. Using these functions with the above methods, we show that a log\,$L$\,--\,log\,$T$ method can provide good constraints on the form of the LF, particularly the high end. Using a sample of 175 peak luminosities determined from redshifts with well defined selection criteria our results suggest that LGRBs occur at a local rate (without beaming corrections) of $[\,0.7 < \rho_{0} < 0.8\,]\,\mathrm{Gpc}^{-3}\mathrm{yr}^{-1}$. Within this range, assuming a broken-power-law LF, we find best estimates for the low and high energy indices of $-0.95 \pm 0.09$ and $-2.59 \pm0.93$ respectively, separated by a break luminosity $0.80 \pm0.43 \times 10^{52}$\,erg\,s$^{-1}$.

We report a detection of a faint near-Earth asteroid (NEA), which was done using our synthetic tracking technique and the CHIMERA instrument on the Palomar 200-inch telescope. This asteroid, with apparent magnitude of 23, was moving at 5.97 degrees per day and was detected at a signal-to-noise ratio (SNR) of 15 using 30 sec of data taken at a 16.7 Hz frame rate. The detection was confirmed by a second observation one hour later at the same SNR. The asteroid moved 7 arcseconds in sky over the 30 sec of integration time because of its high proper motion. The synthetic tracking using 16.7 Hz frames avoided the trailing loss suffered by conventional techniques relying on 30-sec exposure, which would degrade the surface brightness of image on CCD to an approximate magnitude of 25. This detection was a result of our 12-hour blind search conducted on the Palomar 200-inch telescope over two nights on September 11 and 12, 2013 scanning twice over six 5.0 deg x 0.043 deg fields. The fact that we detected only one NEA, is consistent with Harris's estimation of the asteroid population distribution, which was used to predict the detection of 1--2 asteroids of absolute magnitude H=28--31 per night. The design of experiment, data analysis method, and algorithms for estimating astrometry are presented. We also demonstrate a milli-arcsecond astrometry using observations of two bright asteroids with the same system on Apr 3, 2013. Strategies of scheduling observations to detect small and fast-moving NEAs with the synthetic tracking technique are discussed.

We have investigated a sample of 2099 broad absorption line (BAL) quasars with z=1.7-2.2 built from the Sloan Digital Sky Survey Data Release Seven and the Wide-field Infrared Survey. This sample is collected from two BAL quasar samples in the literature, and refined by our new algorithm. Correlations of outflow velocity and strength with hot dust indicator (beta_NIR) and other quasar physical parameters, such as Eddington ratio, luminosity and UV continuum slope, are explored in order to figure out which parameters drive outflows. Here beta_NIR is the near-infrared continuum slope, a good indicator of the amount of hot dust emission relative to accretion disk emission. We confirm previous findings that outflow properties moderately or weakly depends on Eddington ratio, UV slope and luminosity. For the first time, we report moderate and significant correlations of outflow strength and velocity with beta_NIR in BAL quasars. It is consistent with the behavior of blueshifted broad emission lines in non-BAL quasars. The statistical analysis and composite spectra study both reveal that outflow strength and velocity are more strongly correlated with beta_NIR than Eddington ratio, luminosity and UV slope. In particular, the composites show that the entire C IV absorption profile shifts blueward and broadens as beta_NIR increases, while Eddington ratio and UV slope only affect the high and low velocity part of outflows, respectively. We discuss several potential processes and suggest that dusty outflow scenario, i.e. dust is intrinsic to outflows and may contribute to the outflow acceleration, is most likely. The BAL quasar catalog is available from the authors upon request.

We report on Chandra X-ray observations of four candidate low-mass black hole (<10^6Msun) active galactic nuclei (AGNs) that have the estimated Eddington ratios among the lowest (~10^(-2)) found for this class. The aims are to validate the nature of their AGNs and to confirm the low Eddington ratios that are derived from the broad H_alpha line, and to explore this poorly studied regime in the AGN parameter space. Among them, two objects with the lowest significance of the broad lines are also observed with Multi-Mirror Telescope, and the high-quality optical spectra taken confirm them as Seyfert 1 AGNs and as having small black hole masses. X-ray emission is detected from the nuclei of two of the galaxies, which is variable on timescales of 10^3s, whereas no significant (or only marginal at best) detection is found for the remaining two. The X-ray luminosities are on the order of 10^(41) ergs/s or even lower, on the order of 10^(40) ergs/s for non-detections, which are among the lowest regimes ever probed for Seyfert galaxies. The low X-ray luminosities, compared to their black hole masses derived from H_alpha, confirm their low accretion rates assuming typical bolometric corrections. Our results hint at the existence of a possibly large population of under-luminous low-mass black holes in the local universe. An off-nucleus ultra-luminous X-ray source (ULX) in one of the dwarf galaxies is detected serendipitously, with a luminosity (6-9)x10^(39) ergs/s in 2-10 keV.

The High-Altitude Water Cherenkov Gamma Ray Observatory (HAWC) is under construction 4100 meters above sea level at Sierra Negra, Mexico. We describe the design and cabling of the detector, the characterization of the photomultipliers, and the timing calibration system. We also outline a next-generation detector based on the water Cherenkov technique.

The High-Altitude Water Cherenkov (HAWC) Gamma-Ray Observatory is designed to record air showers produced by cosmic rays and gamma rays between 100 GeV and 100 TeV. Because of its large field of view and high livetime, HAWC is well-suited to measure gamma rays from extended sources, diffuse emission, and transient sources. We describe the sensitivity of HAWC to emission from the extended Cygnus region as well as other types of galactic diffuse emission; searches for flares from gamma-ray bursts and active galactic nuclei; and the first measurement of the Crab Nebula with HAWC-30.

We describe measurements of GeV and TeV cosmic rays with the High-Altitude Water Cherenkov Gamma-Ray Observatory, or HAWC. The measurements include the observation of the shadow of the moon; the observation of small-scale and large-scale angular clustering of the TeV cosmic rays; the prospects for measurement of transient solar events with HAWC; and the observation of Forbush decreases with the HAWC engineering array and HAWC-30.

The High-Altitude Water Cherenkov Gamma Ray Observatory (HAWC) is designed to perform a synoptic survey of the TeV sky. The high energy coverage of the experiment will enable studies of fundamental physics beyond the Standard Model, and the large field of view of the detector will enable detailed studies of cosmologically significant backgrounds and magnetic fields. We describe the sensitivity of the full HAWC array to these phenomena in five contributions shown at the 33rd International Cosmic Ray Conference in Rio de Janeiro, Brazil (July 2013).

Correlations of hot dust emission with outflow properties are investigated, based on a large z~2 non-broad absorption lines quasar sample built from the Wide-field Infrared Survey and the Sloan Digital Sky Survey data releases. We use the near infrared slope and the infrared to UV luminosity ratio to indicate the hot dust emission relative to the emission from the accretion disk. In our luminous quasars, these hot dust emission indicators are almost independent of the fundamental parameters, such as luminosity, Eddington ratio and black hole mass, but moderately dependent on the blueshift and asymmetry index (BAI) and full width at half-maximum (FWHM) of CIV lines. Interestingly, the latter two correlations dramatically strengthen with increasing Eddington ratio. We suggest that, in high Eddington ratio quasars, CIV regions are dominated by outflows so the BAI and FWHM(CIV) can reliably reflect the general properties and velocity of outflows, respectively. While in low Eddington ratio quasars, CIV lines are primarily emitted by virialized gas so the BAI and FWHM(CIV) become less sensitive to outflows. Therefore, the correlations for the highest Eddington ratio quasars are more likely to represent the true dependence of hot dust emission on outflows and the correlations for the entire sample are significantly diluted by the low Eddington ratio quasars. Our results show that an outflow with a large BAI or velocity can double the hot dust emission on average. We suggest that outflows either contain hot dust in themselves or interact with the dusty interstellar medium or torus.

A fraction of the heavily reddened quasars require a reddening curve which is even steeper than that of the Small Magellanic Cloud. In this paper, we thoroughly characterize the anomalously steep reddening law in quasars, via an exceptional example observed in IRAS 14026+4341. By comparing the observed spectrum to the quasar composite spectrum, we derive a reddening curve in the rest-frame wavelength range of 1200 {\AA}--10000 {\AA}. It is featured with a steep rise at wavelengths shorter than 3000 {\AA}, but no significant reddening at longer wavelengths. The absence of dust reddening in optical continuum is confirmed by the normal broad-line Balmer decrement (the H$\alpha$/H$\beta$ ratio) in IRAS 14026+4341. The anomalous reddening curve can be satisfactorily reproduced by a dust model containing silicate grains in a power-law size distribution, $dn(a)/da \propto a^{-1.4}$, truncated at a maximum size $a_{max}=70 {\rm nm}$. The unusual size distribution may be caused by the destruction of large "stardust" grains by quasar activities or a different dust formation mechanism (i.e., the in situ formation of dust grains in quasar outflows). It is also possible that the analogies of the dust grains observed toward the Galactic center is responsible for the steep reddening curve. In addition, we find that IRAS 14026+4341 is a weak emission-line quasar (i.e., PHL 1811 analogies) with heavy dust reddening and blueshifted broad absorption lines.