Ammonia is a classic interstellar thermometer. The estimation of the rotational and kinetic temperatures can be affected by the blended Hyperfine Components (HFCs). We developed a new recipe, referred to as the HyperFine Group Ratio (HFGR), which utilizes only direct observables, namely the intensity ratios between the grouped HFCs. As tested on the model spectra, the empirical formulae in HFGR can derive the rotational temperature ($T_{\rm rot}$) from the HFC group ratios in an unambiguous manner. We compared HFGR to the two conventional methods, hyperfine fitting and line intensity ratio, based on both real data and simulated spectra. The HFGR has three major improvements. First, HFGR does not require modeling the HFC or fitting the line profiles, thus is more robust against the effect of HFC blending. Second, the simulation-enabled empirical formulae are much faster than fitting the spectra over the parameter space, so the computer time and human time can be both largely saved. Third, the statistical uncertainty of the temperature $\Delta T_{\rm rot}$ as a function of the signal-to-noise ratio (SNR) are also provided. HFGR can keep an internal error of $\Delta T_{\rm rot}\leq0.5$ K over a broad parameter space of temperature (10 to 70 K), line width (0.3 to 4 \kms), and optical depth (0 to 5). When applied to the noisy spectra, HFGR can keep an uncertainty of $\Delta T_{\rm rot}\leq 1.0$ K (1 $\sigma$) when SNR$>4$.

The "ATOMS," standing for {\it ALMA Three-millimeter Observations of Massive Star-forming regions}, survey has observed 146 active star forming regions with ALMA Band 3, aiming to systematically investigate the spatial distribution of various dense gas tracers in a large sample of Galactic massive clumps, to study the roles of stellar feedback in star formation, and to characterize filamentary structures inside massive clumps. In this work, the observations, data analysis, and example science of the "ATOMS" survey are presented, using a case study for the G9.62+0.19 complex. Toward this source, some transitions, commonly assumed to trace dense gas, including CS $J = 2-1$, HCO$^+$ $J = 1-0$ and HCN $J = 1-0$, are found to show extended gas emission in low density regions within the clump; less than 25\% of their emission is from dense cores. SO, CH$_3$OH, H$^{13}$CN and HC$_3$N show similar morphologies in their spatial distributions and reveal well the dense cores. Widespread narrow SiO emission is present (over $\sim$1 pc), which may be caused by slow shocks from large--scale colliding flows or H{\sc ii} regions. Stellar feedback from an expanding H{\sc ii} region has greatly reshaped the natal clump, significantly changed the spatial distribution of gas, and may also account for the sequential high-mass star formation in the G9.62+0.19 complex. The ATOMS survey data can be jointly analyzed with other survey data, e.g., "MALT90", "Orion B", "EMPIRE", "ALMA\_IMF", and "ALMAGAL", to deepen our understandings of "dense gas" star formation scaling relations and massive proto-cluster formation.

We report studies of the relationships between the total bolometric luminosity ($L_{\rm bol}$ or $L_{\rm TIR}$) and the molecular line luminosities of $J=1-0$ transitions of H$^{13}$CN, H$^{13}$CO$^+$, HCN, and HCO$^+$ with data obtained from ACA observations in the "ATOMS" survey of 146 active Galactic star forming regions. The correlations between $L_{\rm bol}$ and molecular line luminosities $L'_{\rm mol}$ of the four transitions all appear to be approximately linear. Line emission of isotopologues shows as large scatters in $L_{\rm bol}$-$L'_{\rm mol}$ relations as their main line emission. The log($L_{\rm bol}$/$L'_{\rm mol}$) for different molecular line tracers have similar distributions. The $L_{\rm bol}$-to-$L'_{\rm mol}$ ratios do not change with galactocentric distances ($R_{\rm GC}$) and clump masses ($M_{\rm clump}$). The molecular line luminosity ratios (HCN-to-HCO$^+$, H$^{13}$CN-to-H$^{13}$CO$^+$, HCN-to-H$^{13}$CN and HCO$^+$-to-H$^{13}$CO$^+$) all appear constant against $L_{\rm bol}$, dust temperature ($T_{\rm d}$), $M_{\rm clump}$ and $R_{\rm GC}$. Our studies suggest that both the main lines and isotopologue lines are good tracers of the total masses of dense gas in Galactic molecular clumps. The large optical depths of main lines do not affect the interpretation of the slopes in star formation relations. We find that the mean star formation efficiency (SFE) of massive Galactic clumps in the "ATOMS" survey is reasonably consistent with other measures of the SFE for dense gas, even those using very different tracers or examining very different spatial scales.

We studied the filament structures and dense cores in OMC-2,3 region in Orion A North molecular cloud using the high-resolution N2H+ (1-0) spectral cube observed with the Atacama Large Millimeter/Submillimeter Array (ALMA). The filament network over a total length of 2 pc is found to contain 170 intersections and 128 candidate dense cores. The dense cores are all displaced from the infrared point sources (possible young stars), and the major fraction of cores (103) are located around the intersections. Towards the intersections, there is also an increasing trend for the total column density Ntot as well as the the power-law index of the column-density Probability Distribution Function (N-PDF), suggesting that the intersections would in general have more significant gas assembly than the other part of the filament paths. The virial analysis shows that the dense cores mostly have virial mass ratio of alpha_vir=M_vir/M_gas<1.0, suggesting them to be bounded by the self gravity. In the mean time, only about 23 percent of the cores have critical mass ratio of alpha_crit=M_crit/M_gas<1.0, suggesting them to be unstable against core collapse. Combining these results, it shows that the major fraction of the cold starless and possible prestellar cores in OMC-2,3 are being assembled around the intersections, and currently in a gravitationally bound state. But more extensive core collapse and star formation may still require continuous core-mass growth or other perturbatio

We report here Atacama Large Millimeter/submillimeter Array (ALMA) N$_2$H$^+$ (1-0) images of the Orion Molecular Cloud 2 and 3 (OMC-2/3) with high angular resolution (3'' or 1200 au) and high spatial dynamic range. Combining dataset from the ALMA main array, ALMA Compact Array (ACA), the Nobeyama 45m Telescope, and the JVLA (providing temperature measurement on matching scales), we find that most of the dense gas in OMC-2/3 is subsonic ($\rm \sigma_{NT}/c_{s}$ = 0.62) with a mean line width ($\Delta\upsilon$) of 0.39 km s$^{-1}$ FWHM. This is markedly different from the majority of previous observations of massive star-forming regions. In contrast, line widths from the Nobeyama Telescope are transonic at 0.69 km s$^{-1}$ ($\rm \sigma_{NT}/c_{s}$ = 1.08). We demonstrated that the larger line widths obtained by the single-dish telescope arose from unresolved sub-structures within their respective beams. The dispersions from larger scales $\sigma_{ls}$ (as traced by the Nobeyama Telescope) can be decomposed into three components $\rm \sigma_{ls}^2 = \sigma_{ss}^2+ \sigma_{bm}^2+ \sigma_{rd}^2$, where small-scale $\sigma_{ss}$ is the line dispersion of each ALMA beam, bulk motion $\sigma_{bm}$ is dispersion between peak velocity of each ALMA beam, and $\sigma_{rd}$ is the residual dispersion. Such decomposition, though purely empirical, appears to be robust throughout our data cubes. Apparent supersonic line widths, commonly found in massive molecular clouds, are thus likely due to the effect of poor spatial resolution. The observed non-thermal line dispersion (sometimes referred to as 'turbulence') transits from supersonic to subsonic at $\sim 0.05$ pc scales in OMC-2/3 region. Such transition could be commonly found with sufficient spatial (not just angular) resolution, even in regions with massive young clusters, such as Orion molecular clouds studied here.

To investigate the use of saliency-map analysis to aid in searches for transient signals, such as fast radio bursts and individual pulses from radio pulsars. We aim to demonstrate that saliency maps provide the means to understand predictions from machine learning algorithms and can be implemented in piplines used to search for transient events. We have implemented a new deep learning methodology to predict whether or not any segment of the data contains a transient event. The algorithm has been trained using real and simulated data sets. We demonstrate that the algorithm is able to identify such events. The output results are visually analysed via the use of saliency maps. We find that saliency maps can produce an enhanced image of any transient feature without the need for de-dispersion or removal of radio frequency interference. Such maps can be used to understand which features in the image were used in making the machine learning decision and to visualise the transient event. Even though the algorithm reported here was developed to demonstrate saliency-map analysis, we have detected, in archival data, a single burst event with dispersion measure of $41$\,cm$^{-3}$pc that is not associated with any currently known pulsar.

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.

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.

We present an observational study of the sulfur (S)-bearing species towards Orion KL at 1.3 mm by combining ALMA and IRAM-30\,m single-dish data. At a linear resolution of $\sim$800 au and a velocity resolution of 1 $\mathrm{km\, s^{-1}\, }$, we have identified 79 molecular lines from 6 S-bearing species. In these S-bearing species, we found a clear dichotomy between carbon-sulfur compounds and carbon-free S-bearing species in various characteristics, e.g., line profiles, spatial morphology, and molecular abundances with respect to $\rm H_2$. Lines from the carbon-sulfur compounds (i.e., OCS, $^{13}$CS, H$_2$CS) exhibit spatial distributions concentrated around the continuum peaks and extended to the south ridge. The full width at half maximum (FWHM) linewidth of these molecular lines is in the range of 2 $\sim$ 11 $\mathrm{km\, s^{-1}\, }$. The molecular abundances of OCS and H$_2$CS decrease slightly from the cold ($\sim$68 K) to the hot ($\sim$176 K) regions. In contrast, lines from the carbon-free S-bearing species (i.e., SO$_2$, $^{34}$SO, H$_2$S) are spatially more extended to the northeast of mm4, exhibiting broader FWHM linewidths (15 $\sim$ 26 $\mathrm{km\, s^{-1}\, }$). The molecular abundances of carbon-free S-bearing species increase by over an order of magnitude as the temperature increase from 50 K to 100 K. In particular, $\mathrm{^{34}SO/^{34}SO_2}$ and $\mathrm{OCS/SO_2}$ are enhanced from the warmer regions ($>$100 K) to the colder regions ($\sim$50 K). Such enhancements are consistent with the transformation of SO$_2$ at warmer regions and the influence of shocks.

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.

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

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.

We performed a systematic 6.7 GHz Class \uppercase\expandafter{\romannumeral2} methanol maser survey using the Shanghai Tianma Radio Telescope toward targets selected from the all-sky \emph{Wide-Field Infrared Survey Explorer (WISE)} point catalog. In this paper, we report the results from the survey of those at high Galactic latitudes, i.e. $|b|>$ 2$^\circ$. Of 1473 selected \emph{WISE} point sources at high latitude, 17 point positions that were actually associated with 12 sources were detected with maser emission, reflecting the rarity (1$-$2\%) of methanol masers in the region away from the Galactic plane. Out of the 12 sources, 3 are detected for the first time. The spectral energy distribution (SED) at infrared bands shows that these new detected masers occur in the massive star forming regions. Compared to previous detections, the methanol maser changes significantly in both spectral profiles and flux densities. The infrared \emph{WISE} images show that almost all of these masers are located in the positions of the bright \emph{WISE} point sources. Compared with the methanol masers at the Galactic plane, these high-latitude methanol masers provide good tracers for investigating the physics and kinematics around massive young stellar objects, because they are believed to be less affected by the surrounding cluster environment.

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 present the first release of the data and compact-source catalogue for the JCMT Large Program SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE). SCOPE consists of 850-um continuum observations of 1235 Planck Galactic Cold Clumps (PGCCs) made with the Submillimetre Common-User Bolometer Array 2 on the James Clerk Maxwell Telescope. These data are at an angular resolution of 14.4 arcsec, significantly improving upon the 353-GHz resolution of Planck at 5 arcmin, and allowing for a catalogue of 3528 compact sources in 558 PGCCs. We find that the detected PGCCs have significant sub-structure, with 61 per cent of detected PGCCs having 3 or more compact sources, with filamentary structure also prevalent within the sample. A detection rate of 45 per cent is found across the survey, which is 95 per cent complete to Planck column densities of $N_{H_{2}}$ $>$ 5 $\times$ 10$^{21}$ cm$^{-2}$. By positionally associating the SCOPE compact sources with YSOs, the star formation efficiency, as measured by the ratio of luminosity to mass, in nearby clouds is found to be similar to that in the more distant Galactic Plane, with the column density distributions also indistinguishable from each other.

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.

How stellar feedback from high-mass stars (e.g., H{\sc ii} regions) influences the surrounding interstellar medium and regulates new star formation is still unclear. To address this question, we observed the G9.62+0.19 complex in 850 $\mu$m continuum with the JCMT/POL-2 polarimeter. An ordered magnetic field has been discovered in its youngest clump, the G9.62 clump. The magnetic field strength is determined to be $\sim$1 mG. Magnetic field plays a larger role than turbulence in supporting the clump. However, the G9.62 clump is still unstable against gravitational collapse even if thermal, turbulent, and magnetic field support are taken into account all together. The magnetic field segments in the outskirts of the G9.62 clump seem to point toward the clump center, resembling a dragged-in morphology, indicating that the clump is likely undergoing magnetically-regulated global collapse. However, The magnetic field in its central region is aligned with the shells of the photodissociation regions (PDRs) and is approximately parallel to the ionization (or shock) front, indicating that the magnetic field therein is likely compressed by the expanding H{\sc ii} regions that formed in the same complex.

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.

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.

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.

Magnetic field is one of the key agents that play a crucial role in shaping molecular clouds and regulating star formation, yet the complete information on the magnetic field is not well constrained due to the limitations in observations. We study the magnetic field in the massive infrared dark cloud G035.39-00.33 from dust continuum polarization observations at 850 $\micron$ with SCUBA-2/POL-2 at JCMT. The magnetic field tends to be perpendicular to the densest part of the main filament (F$_{M}$), whereas it has a less defined relative orientation in the rest of the structure, where it tends to be parallel to some diffuse regions. A mean plane-of-the-sky magnetic field strength of $\sim$50 $\mu$G for F$_{M}$ is obtained using Davis-Chandrasekhar-Fermi method. Based on $^{13}$CO (1-0) line observations, we suggest a formation scenario of F$_{M}$ due to large-scale ($\sim$10 pc) cloud-cloud collision. Using additional NH$_3$ line data, we estimate that F$_{M}$ will be gravitationally unstable if it is only supported by thermal pressure and turbulence. The northern part of F$_{M}$, however, can be stabilized by a modest additional support from the local magnetic field. The middle and southern parts of F$_{M}$ are likely unstable even if the magnetic field support is taken into account. We claim that the clumps in F$_{M}$ may be supported by turbulence and magnetic fields against gravitational collapse. Finally, we identified for the first time a massive ($\sim$200 M$_{\sun}$), collapsing starless clump candidate, "c8", in G035.39-00.33. The magnetic field surrounding "c8" is likely pinched, hinting at an accretion flow along the filament.

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.

Continuum emissions from dust grains are used as a general probe to constrain the initial physical conditions of molecular dense cores where new stars may born. To get as much information as possible from dust emissions, we have developed a tool, named as $COREGA$, which is capable of identifying positions of dense cores, optimizing a three-dimensional model for the dense cores with well characterized uncertainties. $COREGA$ can also estimate the physical properties of dense cores, such as density, temperature, and dust emissivity, through analyzing multi-wavelength dust continuum data sets. In the numerical tests on $COREGA$, the results of fitting simulated data are consistent with initial built-in parameters. We also demonstrate $COREGA$ by adding random gaussian noises with Monte Carlo methods and show that the results are stable against varying observational noise intensities within certain levels. A beam size $<$ 3 arcsec and rms $<$ 0.2mJy/pixel (1 pixel = 0.1") is needed for ALMA to distinguish different collapse models, such as power law and Bonner-Ebert sphere, during continuum observations of massive dense cores in Orion molecular cloud. Based on its advanced algorithm, $COREGA$ is capable of giving a quick and deep analysis on dust cores.

The axion-like particle (ALP)-photon mixing in the magnetic field around $\gamma$-ray sources or along the line-of-sight could induce oscillation between photons and ALPs, which then causes irregularities in the $\gamma$-ray spectra. In this work we try to search for such spectral irregularities in the spectrum of PKS $2155-304$ using 8.6 years of the Fermi Large Area Telescope (Fermi-LAT) data. No significant evidence for the presence of ALP-photon oscillation is obtained, and the parameter space of ALPs is constrained. The exclusion region sensitively depends on the poorly known magnetic field of host galaxy cluster of PKS $2155-304$. If the magnetic field is as high as $\sim 10~{\rm \mu G}$, the "hole"-like parameter region allowed in Ref.~\cite{Fermi} can be ruled out.

A correlation between the secondary cosmic ray flux and the near-earth electric field intensity, measured during thunderstorms, has been found by analyzing the data of the ARGO-YBJ experiment, a full coverage air shower array located at the Yangbajing Cosmic Ray Laboratory (4300 m a. s. l., Tibet, China). The counting rates of showers with different particle multiplicities, have been found to be strongly dependent upon the intensity and polarity of the electric field measured during the course of 15 thunderstorms. In negative electric fields (i.e. accelerating negative charges downwards), the counting rates increase with increasing electric field strength. In positive fields, the rates decrease with field intensity until a certain value of the field EFmin (whose value depends on the event multiplicity), above which the rates begin increasing. By using Monte Carlo simulations, we found that this peculiar behavior can be well described by the presence of an electric field in a layer of thickness of a few hundred meters in the atmosphere above the detector, which accelerates/decelerates the secondary shower particles of opposite charge, modifying the number of particles with energy exceeding the detector threshold. These results, for the first time, give a consistent explanation for the origin of the variation of the electron/positron flux observed for decades by high altitude cosmic ray detectors during thunderstorms.

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.

How mass is accumulated from cloud-scale down to individual stars is a key open question in understanding high-mass star formation. Here, we present the mass accumulation process in a hub-filament cloud G22 which is composed of four supercritical filaments. Velocity gradients detected along three filaments indicate that they are collapsing with a total mass infall rate of about 440 $M_\odot$ Myr$^{-1}$, suggesting the hub mass would be doubled in six free-fall times, adding up to $ \sim2 $ Myr. A fraction of the masses in the central clumps C1 and C2 can be accounted for through large-scale filamentary collapse. Ubiquitous blue profiles in HCO$^+$ $ (3-2) $ and $^{13}$CO $ (3-2) $ spectra suggest a clump-scale collapse scenario in the most massive and densest clump C1. The estimated infall velocity and mass infall rate are 0.31 km s$^{-1}$ and $ 7.2 \times10^{-4} $ $M_\odot$ yr$^{-1}$, respectively. In clump C1, a hot molecular core (SMA1) is revealed by the SMA observations and an outflow-driving high-mass protostar is located at the center of SMA1. The mass of the protostar is estimated to be $ 11-15 $ $M_\odot$ and it is still growing with an accretion rate of $ 7\times10^{-5} $ $M_\odot$ yr$^{-1}$. The coexistent infall in filaments, clump C1, and the central hot core in G22 suggests that pre-assembled mass reservoirs (i.e., high-mass starless cores) may not be required to form high-mass stars. In the course of high-mass star formation, the central protostar, the core, and the clump can simultaneously grow in mass via core-fed/disk accretion, clump-fed accretion, and filamentary/cloud collapse.

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.

Stellar feedback from high-mass stars (e.g., H{\sc ii} regions) can strongly influence the surrounding interstellar medium and regulate star formation. Our new ALMA observations reveal sequential high-mass star formation taking place within one sub-virial filamentary clump (the G9.62 clump) in the G9.62+0.19 complex. The 12 dense cores (MM 1-12) detected by ALMA are at very different evolutionary stages, from starless core phase to UC H{\sc ii} region phase. Three dense cores (MM6, MM7/G, MM8/F) are associated with outflows. The mass-velocity diagrams of outflows associated with MM7/G and MM8/F can be well fitted with broken power laws. The mass-velocity diagram of SiO outflow associated with MM8/F breaks much earlier than other outflow tracers (e.g., CO, SO, CS, HCN), suggesting that SiO traces newly shocked gas, while the other molecular lines (e.g., CO, SO, CS, HCN) mainly trace the ambient gas continuously entrained by outflow jets. Five cores (MM1, MM3, MM5, MM9, MM10) are massive starless core candidates whose masses are estimated to be larger than 25 M$_{\sun}$, assuming a dust temperature of $\leq$ 20 K. The shocks from the expanding H{\sc ii} regions ("B" \& "C") to the west may have great impact on the G9.62 clump through compressing it into a filament and inducing core collapse successively, leading to sequential star formation. Our findings suggest that stellar feedback from H{\sc ii} regions may enhance the star formation efficiency and suppress the low-mass star formation in adjacent pre-existing massive clumps.

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

The High Altitude Water Cherenkov Observatory (HAWC) has a wide field-of-view (FOV, $\sim$2sr) and a high duty cycle ($\sim$95\%), which make it a powerful survey and monitoring experiment for sources of TeV gamma rays. We present a systematic survey of gamma-ray sources based on the Fermi 3FHL catalog. Sources are restricted to HAWC's FOV (Declination 19$^\circ$ $\pm$ 40$^\circ$) and to extragalactic sources with redshift: 0.001 $<$ z $<$ 0.3. Extragalactic gamma-ray sources are dominated by active galactic nuclei (AGN) and TeV gamma-ray sources are mostly BL Lac-type blazars. The study of AGNs through high energy gamma rays has opened a new window into the extreme processes of particle acceleration in the jets of these objects and provides a way to study the photon propagation and extra-galactic background light. We have improved the HAWC sensitivity at low energies (100 GeV to 1 TeV) based on the Crab pulsar, which is an excellent calibration source for TeV gamma rays. We will present the results of searching for and monitoring nearby AGNs with the improved analysis.

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.

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 ARGO-YBJ detector, located at the Yangbajing Cosmic Ray Laboratory (4300 m a. s. l., Tibet, China), was a full coverage air shower array dedicated to gamma ray astronomy and cosmic ray studies. The wide field of view (~ 2 sr) and high duty cycle (> 86%), made ARGO-YBJ suitable to search for short and unexpected gamma ray emissions like gamma ray bursts (GRBs). Between 2007 November 6 and 2013 February 7, 156 satellite-triggered GRBs (24 of them with known redshift) occurred within the ARGO-YBJ field of view. A search for possible emission associated to these GRBs has been made in the two energy ranges 10-100 GeV and 10-1000 GeV. No significant excess has been found in time coincidence with the satellite detections nor in a time window of one hour after the bursts. Taking into account the EBL absorption, upper limits to the energy fluence at 99% of confidence level have been evaluated,with values ranging from ~ 10-5 erg cm-2 to ~10-1 erg cm-2.

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.

We present a Submillimeter Array (SMA) observation towards the young massive double-core system G350.69-0.49. This system consists of a northeast (NE) diffuse gas Bubble and a southwest (SW) massive young stellar object (MYSO), both clearly seen in the Spitzer images. The SMA observations reveal a gas flow between the NE Bubble and the SW MYSO in a broad velocity range from 5 to 30 km/s with respect to the system velocity. The gas flow is well confined within the interval between the two objects, and traces a significant mass transfer from the NE gas Bubble to the SW massive core. The transfer flow can supply the material accreted onto the SW MYSO at a rate of $4.2\times10^{-4} year^{-1}$. The whole system therefore suggests a mode for the mass growth in MYSO from a gas transfer flow launched from its companion gas clump, despite that the driving mechanism of the transfer flow is not yet fully determined from the current data.

Dark matter may interact with Standard Model (SM) particle through the exchange of a massive spin-2 graviton producing signals that can be detected. In this work we examine the $\gamma$-ray emission signals, including the line emission and the continuous spectrum component in such a massive graviton-mediated dark matter model. The constraints of LHC data, dark matter relic density as well as the dark matter indirect detection data have been applied to narrow down the parameter space. We focus on the vector dark matter model which could produce detectable $\gamma$-ray line signal. It is found that the $\gamma$-ray line data is effective on constraining the model parameters and the ongoing and upcoming space or ground-based $\gamma$-ray experiments can constrain the model further. As for the continuous $\gamma$-ray emission, the total effective annihilation cross section is $\sim 10^{-26} {\rm cm^{3}~s^{-1}}$ except at the region where dark matter mass is around the graviton mass or half of it, which is consistent with current observational data and will be reliably probed by the upcoming CTA.

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.

In this work the relativistic mean field (RMF) FSUGold model extended to include hyperons is employed to study the properties of neutron stars with strong magnetic fields. The chaotic magnetic field approximation is utilized. The effect of the anomalous magnetic moments (AMMs) is also investigated. It is shown that the equation of state (EOS)of neutron star matter is stiffened by the presence of the magnetic field, which increases the maximum mass of neutron star by around 6%. The AMMs only have a small influence on the EOS of neutron star matter, and increase the maximum mass of neutron star by 0.02 M$_{sun}$. Neutral particles are spin polarized due to the presence of the AMMs.

We present a multi-wavelength observational study of the NGC 2024 filament using infrared to sub-millimeter continuum and the NH$_3$ $(1,1)$ and $(2,2)$ inversion transitions centered on FIR-3, the most massive core therein. FIR-3 is found to have no significant infrared point sources in the Spitzer/IRAC bands. But the NH$_3$ kinetic temperature map shows a peak value at the core center with $T_{\rm k}=25$ K which is significantly higher than the surrounding level ($T_{\rm k}=15-19$ K). Such internal heating signature without an infrared source suggests an ongoing core collapse possibly at a transition stage from first hydrostatic core (FHSC) to protostar. The eight dense cores in the filament have dust temperatures between 17.5 and 22 K. They are much cooler than the hot ridge ($T_{\rm d}=55$ K) around the central heating star IRS-2b. Comparison with a dust heating model suggests that the filament should have a distance of $3-5$ pc from IRS-2b. This value is much larger than the spatial extent of the hot ridge, suggesting that the filament is spatially separated from the hot region along the line of sight.

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.

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.

We present a high-resolution $1.5"$ observational study towards two massive dust-and-gas cores, ORI8nw\_2 and ORI2\_6 in Orion Molecular Cloud using the Combined Array for Research in Millimeter-wave Astronomy (CARMA). In both regions the 3.2 mm continuum emission reveals a dense and compact dust core of 1 to 3 solar masses. The cores are estimated to have number densities more than $10^9$ cm$^{-3}$, which are among the highest volume density ever published for the star-forming cores. In both regions the ${\rm N_2H^+}$ shows multiple gas clumps which are spatially displaced from the dense gas traced by the 3.2 mm continuum as well as the ${\rm HCO^+}$. The ${\rm N_2H^+}$ clumps in ORI8nw\_2 are gravitationally unbound. The clumps have masses comparable to the thermal Jeans mass, and are separated with intervals comparable with the Jeans length, thus likely result from thermal Jeans fragmentation. The fragmentation may have prohibited the high-mass star formation in the two cores which originally have larger total masses. In both regions the abundance ratio of ${\rm [N_2H^+]/[HCO^+]}$ appears to be higher than those in infrared dark clouds as a result of CO depletion. The extreme high central density in these Orion cores implies a very small Jeans scale ($<500$ AU). To further reveal the evolution of such dense cores requires higher resolution and sensitivity.

Context: We investigated the molecular gas associated with 6.7 GHz methanol masers throughout the Galaxy using a J=1-0 transition of the CO isotopologues. Methods:Using the 13.7-meter telescope at the Purple Mountain Observatory (PMO), we have obtained ^{12}CO and ^{13}CO (1-0) lines for 160 methanol masers sources from the first to the third Galactic quadrants. We made efforts to resolve the distance ambiguity by careful comparison with the radio continuum and HI 21 cm observations. Results: First, the maser sources show increased ^{13}CO line widths toward the Galactic center, suggesting that the molecular gas are more turbulent toward the Galactic center. This trend can be noticeably traced by the ^{13}CO line width. Second, the ^{12}CO excitation temperature shows a noticeable correlation with the H_2 column density. A possible explanation consistent with the collapse model is that the higher surface-density gas is more efficient to the stellar heating and/or has a higher formation rate of high-mass stars. Third, comparing the IRDCs, the maser sources on average have significantly lower H_2 column densities, moderately higher temperatures, and similar line widths. Fourth, in the mapped regions, 51 ^{13}CO cores have been revealed. Only 17 coincide with the radio continuum emission (F_{cm}>6 mJy), while a larger fraction (30 cores) coincide with the infrared emissions. The IR-bright and radio-bright sources exhibit significantly higher $^{12}$CO excitation temperatures than the IR-faint and radio-faint sources, respectively. Conclusions: The 6.7 GHz masers show a moderately low ionization rate but have a common-existing stellar heating that generates the IR emissions. The relevant properties can be characterized by the ^{12}CO and ^{13}CO (1-0) emissions in several aspects as described above.

A physical model of Geant4-based simulation of galactic cosmic ray (GCR) particles interaction with the lunar surface matter has been developed to investigate the production rate of cosmogenic nuclei. In this model the GCRs, mainly very high energy protons and $\alpha$ particles, bombard the surface of the Moon and produce many secondary particles such as protons and neutrons. The energies of proton and neutron at different depths are recorded and saved into ROOT files, and the analytical expressions for the differential proton and neutron fluxes are obtained through the best-fit procedure under the ROOT software. To test the validity of this model, we calculate the production rates of long-lived nuclei $^{10}$Be and $^{26}$Al in the Apollo 15 long drill core by combining the above differential fluxes and the newly evaluated spallation reaction cross sections. Numerical results show that the theoretical production rates agree quite well with the measured data. It means that this model works well. Therefore, it can be expected that this model can be used to investigate the cosmogenic nuclei in lunar samples returned by Chinese lunar exploration program and can be extended to study other objects, such as the meteorites and the Earth's atmosphere.

We present a multi-wavelength observational study towards the high-mass young stellar object G8.68-0.37. A single massive gas-and-dust core is observed in the (sub)millimeter continuum and molecular line emissions. We fitted the spectral energy distribution (SED) from the dust continuum emission. The best-fit SED suggests the presence of two components with temperature of $T_{\rm d}=20$ K and 120 K, respectively. The core has a total mass of up to $1.5\times10^3$ $M_{\odot}$ and bolometric luminosity of $2.3\times10^4 L_{\odot}$. Both the mass and luminosity are dominated by the cold component ($T_{\rm d}=20$ K). The molecular lines of C$^{18}$O, C$^{34}$S, DCN, and thermally excited CH$_3$OH are detected in this core. Prominent infall signatures are observed in the $^{12}$CO $(1-0)$ and $(2-1)$. We estimated an infall velocity of 0.45 km s$^{-1}$ and mass infall rate of $7\times10^{-4} M_{\odot}$ year$^{-1}$. From the molecular lines, we have found a high DCN abundance and relative abundance ratio to HCN. The overabundant DCN may originate from a significant deuteration in the previous cold pre-protostellar phase. And the DCN should now be rapidly sublimated from the grain mantles to maintain the overabundance in the gas phase.

We present the results of a high resolution study with the Submillimeter Array towards the massive star forming complex G9.62+0.19. Three sub-mm cores are detected in this region. The masses are 13, 30 and 165 M$_{\sun}$ for the northern, middle and southern dust cores, respectively. Infall motions are found with HCN (4-3) and CS (7-6) lines at the middle core (G9.62+0.19 E). The infall rate is $4.3\times10^{-3}~M_{\odot}\cdot$yr$^{-1}$. In the southern core, a bipolar-outflow with a total mass about 26 M$_{\sun}$ and a mass-loss rate of $3.6\times10^{-5}~M_{\odot}\cdot$yr$^{-1}$ is revealed in SO ($8_{7}-7_{7}$) line wing emission. CS (7-6) and HCN (4-3) lines trace higher velocity gas than SO ($8_{7}-7_{7}$). G9.62+0.19 F is confirmed to be the driving source of the outflow. We also analyze the abundances of CS, SO and HCN along the redshifted outflow lobes. The mass-velocity diagrams of the outflow lobes can be well fitted by a single power law. The evolutionary sequence of the cm/mm cores in this region are also analyzed. The results support that UC~H{\sc ii} regions have a higher blue excess than their precursors.

We present a high-resolution study of a massive dense core JCMT 18354-0649S with the Submillimeter Array. The core is mapped with continuum emission at 1.3 mm, and molecular lines including CH$_{3}$OH ($5_{23}$-$4_{13}$) and HCN (3-2). The dust core detected in the compact configuration has a mass of $47 M_{\odot}$ and a diameter of $2\arcsec$ (0.06 pc), which is further resolved into three condensations with a total mass of $42 M_{\odot}$ under higher spatial resolution. The HCN (3-2) line exhibits asymmetric profile consistent with infall signature. The infall rate is estimated to be $2.0\times10^{-3} M_{\odot}\cdot$yr$^{-1}$. The high velocity HCN (3-2) line wings present an outflow with three lobes. Their total mass is $12 M_{\odot}$ and total momentum is $121 M_{\odot}\cdot$km s$^{-1}$, respectively. Analysis shows that the N-bearing molecules especially HCN can trace both inflow and outflow.