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.

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.

Gravitational wave (GW) data can be used to test the parity symmetry of gravity by investigating the difference between left-hand and right-hand circular polarization modes. In this article, we develop a method to decompose the circular polarizations of GWs produced during the inspiralling stage of compact binaries, with the help of stationary phase approximation. The foremost advantage is that this method is simple, clean, independent of GW waveform, and is applicable to the existing detector network. Applying it to the mock data, we test the parity symmetry of gravity by constraining the velocity birefringence of GWs. If a nearly edge-on binary neutron-stars with observed electromagnetic counterparts at 40 Mpc is detected by the second-generation detector network, one could derive the model-independent test on the parity symmetry in gravity: the lower limit of the energy scale of parity violation can be constrained within $\mathcal{O}(10^4{\rm eV})$.

We propose a compact binary model with an eccentric orbit to explain periodically active fast radio burst (FRB) sources, where the system consists of a neutron star (NS) with strong dipolar magnetic fields and a magnetic white dwarf (WD). In our model, the WD fills its Roche lobe at periastron, and mass transfer occurs from the WD to the NS around this point. The accreted material may be fragmented into a number of parts, which arrive at the NS at different times. The fragmented magnetized material may trigger magnetic reconnection near the NS surface. The electrons can be accelerated to an ultra-relativistic speed, and therefore the curvature radiation of the electrons can account for the burst activity. In this scenario, the duty cycle of burst activity is related to the orbital period of the binary. We show that such a model may work for duty cycles roughly from ten minutes to two days. For the recently reported 16.35-day periodicity of FRB 180916.J0158+65, our model does not naturally explain such a long duty cycle, since an extremely high eccentricity ($e>0.95$) is required.

We present new observations of the active star-formation region NGC 1333 in the Perseus molecular cloud complex from the James Clerk Maxwell Telescope (JCMT) B-Fields In Star-forming Region Observations (BISTRO) survey with the POL-2 instrument. The BISTRO data cover the entire NGC 1333 complex (~ 1.5 pc x 2 pc) at 0.02 pc resolution and spatially resolve the polarized emission from individual filamentary structures for the first time. The inferred magnetic field structure is complex as a whole, with each individual filament aligned at different position angles relative to the local field orientation. We combine the BISTRO data with low- and high- resolution data derived from Planck and interferometers to study the multi-scale magnetic field structure in this region. The magnetic field morphology drastically changes below a scale of ~ 1 pc and remains continuous from the scales of filaments (~ 0.1 pc) to that of protostellar envelopes (~ 0.005 pc or ~ 1000 au). Finally, we construct simple models in which we assume that the magnetic field is always perpendicular to the long axis of the filaments. We demonstrate that the observed variation of the relative orientation between the filament axes and the magnetic field angles are well reproduced by this model, taking into account the projection effects of the magnetic field and filaments relative to the plane of the sky. These projection effects may explain the apparent complexity of the magnetic field structure observed at the resolution of BISTRO data towards the filament network.

We report the detection of four new hot corino sources, G211.47-19.27S, G208.68-19.20N1, G210.49-19.79W and G192.12-11.10 from a survey study of Planck Galactic Cold Clumps in the Orion Molecular Cloud Complex with the Atacama Compact Array (ACA). Three sources had been identified as low mass Class 0 protostars in the Herschel Orion Protostar Survey (HOPS). One source in the lambda Orionis region is firstly reported as a protostellar core. We have observed abundant complex organic molecules (COMs), primarily methanol but also other oxygen-bearing COMs (in G211.47-19.27S and G208.68-19.20N1) and the molecule of prebiotic interest NH2CHO (in G211.47-19.27S), signifying the presence of hot corinos. While our spatial resolution is not sufficient for resolving most of the molecular emission structure, the large linewidth and high rotational temperature of COMs suggest that they likely reside in the hotter and innermost region immediately surrounding the protostar. In G211.47-19.27S, the D/H ratio of methanol ([CH2DOH]/[CH3OH]) and the 12C/13C ratio of methanol ([CH3OH]/[13CH3OH]) are comparable to those of other hot corinos. Hydrocarbons and long carbon-chain molecules such as c-C3H2 and HCCCN are also detected in the four sources, likely tracing the outer and cooler molecular envelopes.

In this paper, we present a cosmological model-independent test for the cosmic opacity at high redshifts ($z\sim5$). We achieve this with the opacity-dependent luminosity distances derived from nonlinear relation between X-ray and UV emissions of quasars, combined with two types of opacity-independent luminosity distances derived from the Hubble parameter measurements and simulated gravitational wave (GW) events achievable with the Einstein Telescope (ET). In the framework of two phenomenological parameterizations adopted to describe cosmic opacity at high redshifts, our main results show that a transparent universe is supported by the current observational data at 2$\sigma$ confidence level. However, the derived value of the cosmic opacity is slightly sensitive to the parametrization of $\tau(z)$, which highlights the importance of choosing a reliable parametrization to describe the optical depth $\tau(z)$ in the early universe. Compared with the previous works, the combination of the quasar data and the $H(z)$/GW observations in similar redshift ranges provides a novel way to confirm a transparent universe ($\epsilon=0$ at higher redshifts $z\sim 5$), with an accuracy of $\Delta \epsilon\sim 10^{-2}$. More importantly, our findings indicate that a strong degeneracy between the cosmic opacity parameter and the parameters characterizing the $L_{UV}-L_X$ relation of quasars, which reinforces the necessity of proper calibration for such new type of high-redshift standard candle (in a cosmological model-independent way).

The Tianlai Cylinder Pathfinder is a radio interferometer array designed to test techniques for 21 cm intensity mapping in the post-reionization Universe, with the ultimate aim of mapping the large scale structure and measuring cosmological parameters such as the dark energy equation of state. Each of its three parallel cylinder reflectors is oriented in the north-south direction, and the array has a large field of view. As the Earth rotates, the northern sky is observed by drift scanning. The array is located in Hongliuxia, a radio-quiet site in Xinjiang, and saw its first light in September 2016. In this first data analysis paper for the Tianlai cylinder array, we discuss the sub-system qualification tests, and present basic system performance obtained from preliminary analysis of the commissioning observations during 2016-2018. We show typical interferometric visibility data, from which we derive the actual beam profile in the east-west direction and the frequency band-pass response. We describe also the calibration process to determine the complex gains for the array elements, either using bright astronomical point sources, or an artificial on site calibrator source, and discuss the instrument response stability, crucial for transit interferometry. Based on this analysis, we find a system temperature of about 90 K, and we also estimate the sensitivity of the array.

Scaling laws of dust, HI gas and metal mass with stellar mass, specific star formation rate and metallicity are crucial to our understanding of the buildup of galaxies through their enrichment with metals and dust. In this work, we analyse how the dust and metal content varies with specific gas mass ($M_{\text{HI}}$/$M_{\star}$) across a diverse sample of 423 nearby galaxies. The observed trends are interpreted with a set of Dust and Element evolUtion modelS (DEUS) - incluidng stellar dust production, grain growth, and dust destruction - within a Bayesian framework to enable a rigorous search of the multi-dimensional parameter space. We find that these scaling laws for galaxies with $-1.0\lesssim \log M_{\text{HI}}$/$M_{\star}\lesssim0$ can be reproduced using closed-box models with high fractions (37-89$\%$) of supernova dust surviving a reverse shock, relatively low grain growth efficiencies ($\epsilon$=30-40), and long dus lifetimes (1-2\,Gyr). The models have present-day dust masses with similar contributions from stellar sources (50-80\,$\%$) and grain growth (20-50\,$\%$). Over the entire lifetime of these galaxies, the contribution from stardust ($>$90\,$\%$) outweighs the fraction of dust grown in the interstellar medium ($<$10$\%$). Our results provide an alternative for the chemical evolution models that require extremely low supernova dust production efficiencies and short grain growth timescales to reproduce local scaling laws, and could help solving the conundrum on whether or not grains can grow efficiently in the interstellar medium.

The cosmic curvature, a fundamental parameter for cosmology could hold deep clues to inflation and cosmic origins. We propose an improved model-independent method to constrain the cosmic curvature by combining the constructed Hubble diagram of high-redshift quasars with galactic-scale strong lensing systems expected to be seen by the forthcoming LSST survey. More specifically, the most recent quasar data are used as a new type of standard candles in the range $0.036<z<5.100$, whose luminosity distances can be directly derived from the non-linear relation between X-ray and UV luminosities. Compared with other methods, the proposed one involving the quasar data achieves constraints with higher precision ($\Delta \Omega_k\sim 10^{-2}$) at high redshifts ($z\sim 5.0$). We also investigate the influence of lens mass distribution in the framework of three types of lens models extensively used in strong lensing studies (SIS model, power-law spherical model, and extended power-law lens model), finding the strong correlation between the cosmic curvature and the lens model parameters. When the power-law mass density profile is assumed, the most stringent constraint on the cosmic curvature $\Omega_k$ can be obtained. Therefore, the issue of mass density profile in the early-type galaxies is still a critical one that needs to be investigated further.

Theory predicts that a supermassive black hole binary (SMBHB) could be observed as a luminous active galactic nucleus (AGN) that periodically varies on the order of its orbital timescale. In X-rays, periodic variations could be caused by mechanisms including relativistic Doppler boosting and shocks. Here we present the first systematic search for periodic AGNs using $941$ hard X-ray light curves (14-195 keV) from the first 105 months of the Swift Burst Alert Telescope (BAT) survey (2004-2013). We do not find evidence for periodic AGNs in Swift-BAT, including the previously reported SMBHB candidate MCG+11$-$11$-$032. We find that the null detection is consistent with the combination of the upper-limit binary population in AGNs in our adopted model, their expected periodic variability amplitudes, and the BAT survey characteristics. We have also investigated the detectability of SMBHBs against normal AGN X-ray variability in the context of the eROSITA survey. Under our assumptions of a binary population and the periodic signals they produce which have long periods of hundreds of days, up to $13$% true periodic binaries can be robustly distinguished from normal variable AGNs with the ideal uniform sampling. However, we demonstrate that realistic eROSITA sampling is likely to be insensitive to long-period binaries because longer observing gaps reduce their detectability. In contrast, large observing gaps do not diminish the prospect of detecting binaries of short, few-day periods, as 19% can be successfully recovered, the vast majority of which can be identified by the first half of the survey.

In direct analogy to the 4-body decay of a heavy scalar particle, the 4-point correlation function of primordial fluctuations carries P and CP information. The CP violation appears as a P-odd angular dependence in the imaginary part of the trispectrum in momentum space. We construct a model with axion-like couplings which leads to observably large CP-violating trispectrum for future surveys. Furthermore, we show the importance of on-shell particle production in observing P- and CP-violating signals. It is impossible to observe these signals from local 4-scalar EFT operators that respect dS isometries, and thus any such observation can rule out single-field EFT with sufficiently small slow-roll parameters. This calculation opens a new frontier of studying P and CP at very high energy scales.

We present results of continuum and spectral line observations with ALMA and 22 GHz water (H$_2$O) maser observations using KaVA and VERA toward a high-mass star-forming region, G25.82-0.17. Multiple 1.3 mm continuum sources are revealed, indicating the presence of young stellar objects (YSOs) at different evolutionary stages, namely an ultra-compact HII region, G25.82-E, a high-mass young stellar object (HM-YSO), G25.82-W1, and starless cores, G25.82-W2 and G25.82-W3. Two SiO outflows, at N-S and SE-NW orientations, are identified. The CH$_3$OH 8$_{-1}$-7$_{0}$ E line, known to be a class I CH$_3$OH maser at 229 GHz is also detected showing a mixture of thermal and maser emission. Moreover, the H$_2$O masers are distributed in a region ~0.25" shifted from G25.82-W1. The CH$_3$OH 22$_{4}$-21$_{5}$ E line shows a compact ring-like structure at the position of G25.82-W1 with a velocity gradient, indicating a rotating disk or envelope. Assuming Keplerian rotation, the dynamical mass of G25.82-W1 is estimated to be $>$25 M$_{\odot}$ and the total mass of 20 M$_\odot$-84 M$_\odot$ is derived from the 1.3 mm continuum emission. The driving source of the N-S SiO outflow is G25.82-W1 while that of the SE-NW SiO outflow is uncertain. Detection of multiple high-mass starless$/$protostellar cores and candidates without low-mass cores implies that HM-YSOs could form in individual high-mass cores as predicted by the turbulent core accretion model. If this is the case, the high-mass star formation process in G25.82 would be consistent with a scaled-up version of low-mass star formation.

We mapped two molecular cloud cores in the Orion A cloud with the ALMA ACA 7-m Array and with the Nobeyama 45-m radio telescope. These cores have bright N$_2$D$^+$ emission in single-pointing observations with the Nobeyama 45-m radio telescope, have relatively high deuterium fraction, and are thought to be close to the onset of star formation. One is a star-forming core, and the other is starless. These cores are located along filaments observed in N$_2$H$^+$, and show narrow linewidths of 0.41 km s$^{-1}$ and 0.45 km s$^{-1}$ in N$_2$D$^+$, respectively, with the Nobeyama 45-m telescope. Both cores were detected with the ALMA ACA 7m Array in the continuum and molecular lines at Band 6. The starless core G211 shows clumpy structure with several sub-cores, which in turn show chemical differences. Also, the sub-cores in G211 have internal motions that are almost purely thermal. The starless sub-core G211D, in particular, shows a hint of the inverse P Cygni profile, suggesting infall motion. The star-forming core G210 shows an interesting spatial feature of two N$_2$D$^+$ peaks of similar intensity and radial velocity located symmetrically with respect to the single dust continuum peak. One interpretation is that the two N$_2$D$^+$ peaks represent an edge-on pseudo-disk. The CO outflow lobes, however, are not directed perpendicular to the line connecting both N$_2$D$^+$ peaks.

The proposed ARIANNA-200 neutrino detector, located at sea-level on the Ross Ice Shelf, Antarctica, consists of 200 autonomous and independent detector stations separated by 1 kilometer in a uniform triangular mesh, and serves as a pathfinder mission for the future IceCube-Gen2 project. The primary science mission of ARIANNA-200 is to search for sources of neutrinos with energies greater than 10^17 eV, complementing the reach of IceCube. An ARIANNA observation of a neutrino source would provide strong insight into the enigmatic sources of cosmic rays. ARIANNA observes the radio emission from high energy neutrino interactions in the Antarctic ice. Among radio based concepts under current investigation, ARIANNA-200 would uniquely survey the vast majority of the southern sky at any instant in time, and an important region of the northern sky, by virtue of its location on the surface of the Ross Ice Shelf in Antarctica. The broad sky coverage is specific to the Moore's Bay site, and makes ARIANNA-200 ideally suited to contribute to the multi-messenger thrust by the US National Science Foundation, Windows on the Universe - Multi-Messenger Astrophysics, providing capabilities to observe explosive sources from unknown directions. The ARIANNA architecture is designed to measure the angular direction to within 3 degrees for every neutrino candidate, which too plays an important role in the pursuit of multi-messenger observations of astrophysical sources.

Because the binary neutron star (BNS) merger GW170817 and its electromagnetic (EM) counterparts were successfully detected, there is a chance to explore the joint effect of the host galaxy and our Milky Way on the weak equivalence principle (WEP) test. In this paper, through using the Navarro-Frenk-White profile and the Herquist profile, we present a new analytic model to calculate the galactic gravitational potential, in which the possible locations of the source by the observed angle offset and the second SN (SN2) kick is accounted for. We show that the upper limit of $\Delta \gamma$ is $10^{-9}$ for the comparison of GW170817/GRB 170817A, and is $10^{-4}$ for the comparison of GW170817/AT2017gfo, which are tighter by one or two orders of magnitude than these of only by the measured Milky Way gravitational potential in literatures. We demonstrate that the WEP test could be strengthened by the contribution of the host galaxy NGC 4993 to the Shapiro time delay. Meanwhile, we also find the large natal kick can produce a maximum deviation of about $20\%$ on the results with a typical kick velocity $400\sim 500 km s^{-1}$. Finally, we study and analyse the impact from the halo mass of NGC 4993 with a typical 0.2 dex uncertainty, and find the upper limit of $\Delta \gamma$ of the maximum mass $10^{12.4}h^{-1} M_{\odot}$ is nearly two times tighter than that of the minimum mass $10^{12.0}h^{-1} M_{\odot}$. It suggests the constraint on WEP test would be also enhanced when the halo mass increases in the host galaxy NGC 4993.

We report results of a project to map HCN and HCO+ J = 1-0 emission toward a sample of molecular clouds in the inner Galaxy, all containing dense clumps that are actively engaged in star formation. We compare these two molecular line tracers with millimeter continuum emission and extinction, as inferred from 13CO, as tracers of dense gas in molecular clouds. The fraction of the line luminosity from each tracer that comes from the dense gas, as measured by AV > 8 mag, varies substantially from cloud to cloud. In all cases, a substantial fraction (in most cases, the majority) of the total luminosity arises in gas below the AV > 8 mag threshold and outside the region of strong mm continuum emission. Measurements of the luminosity of HCN toward other galaxies will likely be dominated by such gas at lower surface density. Substantial, even dominant, contributions to the total line luminosity can arise in gas with densities typical of the cloud as a whole (densities about 100 per cubic cm). Defining the dense clump from the HCN or HCO+ emission itself, similarly to previous studies, leads to a wide range of clump properties, with some being considerably larger and less dense than in previous studies. HCN and HCO+ have similar ability to trace dense gas for the clouds in this sample. For the two clouds with low virial parameters, the 13CO is definitely a worse tracer of the dense gas, but for the other four, it is equally good (or bad) at tracing dense gas.

Filament S242 is 25 pc long with massive clumps and YSO clusters concentrated in its end regions; it is considered a good example of edge collapse. We mapped this filament in the $^{12}$CO(1-0) and $^{13}$CO(1-0) lines. A large-scale velocity gradient along filament S242 has been detected; the relative velocity between the two end-clumps is $\sim$ 3 km s$^{-1}$, indicating an approaching motion between them. These signatures are consistent with the filament S242 being formed through the collapse of a single elongated entity, where an effect known as "gravitational focusing" drives the ends of the filament to collapse (edge collapse). Based on this picture, we estimate a collapse timescale of $\sim$ 4.2 Myr, which is the time needed for a finite and elongated entity evolving to the observed filament S242. For the whole filament, we find that increases in surface densities lead to increases in velocity dispersion, which can be consistently explained as the result of self-gravity. We also calculated the contribution of longitudinal collapse to the observed velocity dispersion and found it to be the dominant effect in driving the gas motion near the end-clumps. We propose that our filament S242 is formed through a two-stage collapse model, where the edge collapse of a truncated filament is followed by a stage of longitudinal accretion toward the dense end-clumps.

Gamma-ray bursts (GRBs) might be powered by a black hole (BH) hyperaccretion systems via the Blandford-Znajek (BZ) mechanism or neutrino annihilation from neutrino-dominated accretion flows (NDAFs). Magnetic coupling (MC) between the inner disc and BH can transfer angular momentum and energy from the fast-rotating BH to the disc. The neutrino luminosity and neutrino annihilation luminosity are both efficiently enhanced by the MC process. In this paper, we study the structure, luminosity, MeV neutrinos, and gravitational waves (GWs) of magnetized NDAFs (MNDAFs) under the assumption that both the BZ and MC mechanisms are present. The results indict that the BZ mechanism will compete with the neutrino annihilation luminosity to trigger jets under the different partitions of the two magnetic mechanisms. The typical neutrino luminosity and annihilation luminosity of MNDAFs are definitely higher than those of NDAFs. The typical peak energy of neutrino spectra of MNDAFs is higher than that of NDAFs, but similar to those of core-collapse supernovae. Moreover, if the MC process is dominant, then the GWs originating from the anisotropic neutrino emission will be stronger particularly for discs with high accretion rates.

We present a cosmic ray muon tomographic experiment for measuring the muon flux and imaging the tunnel overburden structures in Changshu, China. The device used in this study is a tracking detector based on the plastic scintillator with SiPM technology, which can be conveniently operated in field works. The compact system with sensitive area of $6400 cm^2$ can measure the angular distribution of cosmic muons. It's able to image the overburden density length from the surface of overburden to the detector along the muon tracks. The open sky muon flux measurement outside the tunnel has a good agreement with the modified Gassier Formula model. The distributions of muon flux at three positions inside the tunnel are very similar to that of open sky. Assuming the average density of overburden compact sandstone is $2.65 g/cm^3$, the overburden thickness can be obtained from the density length derived from the difference of muon flux inside and outside the tunnel. Moreover, for known penetrated lengths (i.e., topography of overburden), the density anomalies of the overburden can also been obtained. This study suggests a potential application for imaging and detecting subsurface structures in civil engineering, tunnels or caverns with the cosmic ray muon telescope.

We present the properties of magnetic fields (B-fields) in two clumps (clump 1 and clump 2), located at the waist of the bipolar H II region Sh2-201, based on JCMT SCUBA-2/POL-2 observations of 850 $\mu$m polarized dust emission. We find that B-fields in the direction of the clumps are bent and compressed, showing bow-like morphologies, which we attribute to the feedback effect of the H II region on the surface of the clumps. Using the modified Davis-Chandrasekhar-Fermi method we estimate B-fields strengths of 266 $\mu$G and 65 $\mu$G for clump 1 and clump 2, respectively. From virial analyses and critical mass ratio estimates, we argue that clump 1 is gravitationally bound and could be undergoing collapse, whereas clump 2 is unbound and stable. We hypothesize that the interplay between thermal pressure imparted by the H II region, B-field morphologies, and the various internal pressures of the clumps (such as magnetic, turbulent, and gas thermal pressure), has the following consequences: (a) formation of clumps at the waist of the H II region; (b) progressive compression and enhancement of the B-fields in the clumps; (c) stronger B-fields will shield the clumps from erosion by the H II region and cause pressure equilibrium between the clumps and the H II region, thereby allowing expanding I-fronts to blow away from the filament ridge, forming bipolar H II regions; and (d) stronger B-fields and turbulence will be able to stabilize the clumps. A study of a larger sample of bipolar H II regions would help to determine whether our hypotheses are widely applicable.

Testing gravity by binary pulsars nowadays becomes a key issue. Screened modified gravity is a kind of scalar-tensor theory with screening mechanism in order to satisfy the tight Solar System tests. In this paper, we investigate how the screening mechanism affects the orbital dynamics of binary pulsars, and calculate in detail the five post-Keplerian (PK) parameters in this theory. These parameters differ from those of general relativity (GR), and the differences are quantified by the scalar charges, which lead to the dipole radiation in this theory. We combine the observables of PK parameters for the ten binary pulsars, respectively, to place the constraints on the scalar charges and possible deviations from GR. The dipole radiation in the neutron star (NS) - white dwarf (WD) binaries leads to more stringent constraints on deviations from GR. The most constraining systems for the scalar charges of NS and WD are PSR~B1913$+$16 and PSR~J1738$+$0333, respectively. The results of all tests exclude significant strong-field deviations and show good agreement with GR.

General relativity reproduces main current cosmological observations, assuming the validity of cosmic distance duality relation (CDDR) at all scales and epochs. However, CDDR is poorly tested in the redshift interval between the farthest observed Type Ia supernovae (SN Ia) and that of the Cosmic Microwave background (CMB). We present a new idea of testing the validity of CDDR, through the multiple measurements of high-redshift quasars. Luminosity distances are derived from the relation between the UV and X-ray luminosities of quasars, while angular diameter distances are obtained from the compact structure in radio quasars. This will create a valuable opportunity where two different cosmological distances from the same kind of objects at high redshifts are compared. Our constraints are more stringent than other currently available results based on different observational data and show no evidence for the deviation from CDDR at $z\sim 3$. Such accurate model-independent test of fundamental cosmological principles can become a milestone in precision cosmology.

Active galactic nuclei (AGNs) have long been observed to "twinkle" (i.e., their brightness varies with time) on timescales from days to years in the UV/optical bands. Such AGN UV/optical variability is essential for probing the physics of supermassive black holes (SMBHs), the accretion disk, and the broad-line region. Here we show that the temperature fluctuations of an AGN accretion disk, which is magnetically coupled with the corona, can account for observed high-quality AGN optical light curves. We calculate the temperature fluctuations by considering the gas physics of the accreted matter near the SMBH. We find that the resulting simulated AGN UV/optical light curves share the same statistical properties as the observed ones as long as the dimensionless viscosity parameter $\alpha$, which is widely believed to be controlled by magnetohydrodynamic (MHD) turbulence in the accretion disk, is about $0.01$---$0.2$. Moreover, our model can simultaneously explain the larger-than-expected accretion disk sizes and the dependence of UV/optical variability upon wavelength for NGC 5548. Our model also has the potential to explain some other observational facts of AGN UV/optical variability, including the timescale-dependent bluer-when-brighter color variability and the dependence of UV/optical variability on AGN luminosity and black hole mass. Our results also demonstrate a promising way to infer the black-hole mass, the accretion rate, and the radiative efficiency, thereby facilitating understanding of the gas physics and MHD turbulence near the SMBH and its cosmic mass growth history by fitting the AGN UV/optical light curves in the era of time-domain astronomy.

We present a catalog of 316 trans-Neptunian bodies detected by the Dark Energy Survey (DES). These objects include 245 discoveries by DES (139 not previously published) detected in $\approx 60,000$ exposures from the first four seasons of the survey ("Y4" data). The survey covers a contiguous 5000 deg$^2$ of the southern sky in the $grizY$ optical/NIR filter set, with a typical TNO in this part of the sky being targeted by $25-30$ Y4 exposures. We describe the processes for detection of transient sources and the linkage into TNO orbits, which are made challenging by the absence of the few-hour repeat observations employed by TNO-optimized surveys. We also describe the procedures for determining detection efficiencies vs. magnitude and estimating rates of false-positive linkages. This work presents all TNOs which were detected on $\ge 6$ unique nights in the Y4 data and pass a "sub-threshold confirmation" test wherein we demand the the object be detectable in a stack of the individual images in which the orbit indicates an object should be present, but was not detected. This eliminates false positives and yields TNO detections complete to $r\lesssim 23.3$ mag with virtually no dependence on orbital properties for bound TNOs at distance $30\,{\rm AU}<d<2500\,{\rm AU}.$ The final DES TNO catalog is expected to yield $>0.3$ mag more depth, and arcs of $>4$ years for nearly all detections.

Offsets of molecular line emission peaks from continuum peaks are very common but frequently difficult to explain with a single spherical cloud chemical model. We propose that the spatial projection effects of an irregular three dimensional (3D) cloud structure can be a solution. This work shows that the idea can be successfully applied to the Planck cold clump G224.4-0.6 by approximating it with four individual spherically symmetric cloud cores whose chemical patterns overlap with each other to produce observable line maps. With the empirical physical structures inferred from the observation data of this clump and a gas-grain chemical model, the four cores can satisfactorily reproduce its 850 $\mu$m continuum map and the diverse peak offsets of CCS, HC$_3$N and N$_2$H$^+$ simultaneously at chemical ages of about $8\times 10^5\sim 3\times 10^6$ yrs. The 3D projection effects on chemistry has the potential to explain such asymmetrical distributions of chemicals in many other molecular clouds.

We look to provide a detailed description of the SPectroscopic IDentification of ERosita Sources (SPIDERS) survey, an SDSS-IV programme aimed at obtaining spectroscopic classification and redshift measurements for complete samples of sufficiently bright X-ray sources. We describe the SPIDERS X-ray Point Source Spectroscopic Catalogue, considering its store of 11,092 observed spectra drawn from a parent sample of 14,759 ROSAT and XMM sources over an area of 5,129 deg$^2$ covered in SDSS-IV by the eBOSS survey. This programme represents the largest systematic spectroscopic observation of an X-ray selected sample. A total of 10,970 (98.9\%) of the observed objects are classified and 10,849 (97.8\%) have secure redshifts. The majority of the spectra (10,070 objects) are active galactic nuclei (AGN), 522 are cluster galaxies, and 294 are stars. The observed AGN redshift distribution is in good agreement with simulations based on empirical models for AGN activation and duty cycle. Forming composite spectra of type 1 AGN as a function of the mass and accretion rate of their black holes reveals systematic differences in the H-beta emission line profiles. This study paves the way for systematic spectroscopic observations of sources that are potentially to be discovered in the upcoming eROSITA survey over a large section of the sky.

Astrophysical observations provide a unique opportunity to test possible signatures of Lorentz Invariance Violation (LIV), due to the high energies and long distances involved. In quantum theory of gravity, one may expect the modification of the dispersion relation between energy and momentum for photons, which can be probed with the time-lag (the arrival time delay between light curves in different energy bands) of Gamma-ray bursts (GRBs). In this paper, by using the detailed time-delay measurements of GRB 160625B at different energy bands, as well as 23 time-delay GRBs covering the redshifts range of $z=0.168-2.5$ (which were measured at different energy channels from the light curves), we propose an improved model-independent method (based on the newly-compiled sample of $H(z)$ measurements) to probe the energy-dependent velocity due to the modified dispersion relation for photons. In the framework of a more complex and reasonable theoretical expression to describe the time delays, our results imply that the intrinsic time lags can be better described with more GRBs time delay data. More importantly, through direct fitting of the time-delay measurements of a sample of GRBs, our limit on the LIV energy scale is comparable to that with unknown constant for the intrinsic time lag, much lower than the Planck energy scale in both linear LIV and quadratic LIV cases.

We present a statistical study on the orientation of outflows with respect to large-scale filaments and the magnetic fields. Although filaments are widely observed toward Galactic star-forming regions, the exact role of filaments in star formation is unclear. Studies toward low-mass star-forming regions revealed both preferred and random orientation of outflows respective to the filament long-axes, while outflows in massive star-forming regions mostly oriented perpendicular to the host filaments, and parallel to the magnetic fields at similar physical scales. Here, we explore outflows in a sample of 11 protoclusters in HII regions, a more evolved stage compared to IRDCs, using ALMA CO (3-2) line observations. We identify a total of 105 outflow lobes in these protoclusters. Among the 11 targets, 7 are embedded within parsec-scale filamentary structures detected in $^{13}$CO line and 870 $\mu m$ continuum emissions. The angles between outflow axes and corresponding filaments ($\gamma_\mathrm{Fil}$) do not show any hint of preferred orientations (i.e., orthogonal or parallel as inferred in numerical models) with respect to the position angle of the filaments. Identified outflow lobes are also not correlated with the magnetic fields and Galactic plane position angles. Outflows associated with filaments aligned along the large-scale magnetic fields are also randomly orientated. Our study presents the first statistical results of outflow orientation respective to large-scale filaments and magnetic fields in evolved massive star-forming regions. The random distribution suggests a lack of alignment of outflows with filaments, which may be a result of the evolutionary stage of the clusters.

In this paper, we study in detail the evolution of a mini-sigmiod originating in a cross-equatorial coronal hole, where the magnetic field is mostly open and seriously distinct from the closed background field above active-region sigmoids. The source region first appeared as a bipole, which subsequently experienced a rapid emergence followed by a long-term decay. Correspondingly, the coronal structure initially appeared as arc-like loops, then gradually sheared and transformed into continuously sigmoidal loops, mainly owing to flux cancellation near the polarity inversion line. The temperature of J-shaped and sigmoidal loops is estimated to be about $2.0\times10^{6}$ K, greater than that of the background coronal hole. Using the flux-rope insertion method, we further reconstruct the nonlinear force-free fields that well reproduces the transformation of the potential field into a sigmoidal field. The fact that the sheared and sigmoidal loops are mainly concentrated at around the high-Q region implies that the reconnection most likely takes place there to form the sigmoidal field and heat the plasma. Moreover, the twist of sigmoidal field lines is estimated to be around 0.8, less than the values derived for the sigmoids from active regions. However, the sigmoidal flux may quickly enter an unstable regime at the very low corona ($<10\ Mm$) due to the open background field. The results suggest that the mini-sigmoid, at least the one in our study, has the same formation and eruption process as the large-scale one, but is significantly influenced by the overlying flux.

As progenitors of gamma-ray bursts (GRBs), core collapse of massive stars and coalescence of compact object binaries are believed to be powerful sources of gravitational waves (GWs). In the collapsar scenario, a rotating stellar-mass black hole (BH) surrounded by a hyperaccretion disk might be running in the center of a massive collapsar, which is one of the plausible central engines of long GRBs. Such a BH hyperaccretion disk would be in a state of a neutrino-dominated accretion flow (NDAF) at the initial stage of the accretion process; meanwhile, the jets attempt to break out from the envelope and circumstellar medium to power GRBs. In addition to collapsars, the BH hyperaccretion systems are important sources of neutrinos and GWs. In this paper, we investigated the GW emission generated by the anisotropic neutrino emission from NDAFs in the collapsar scenarios. As the results indicate, the typical frequency of GWs is $\sim$ 1-100 Hz, and the masses and metallicities of the progenitor stars have slight effects on the GW strains. The GWs from NDAFs might be detected by operational or planned detectors at the distance of 10 kpc. Moreover, comparisons of the detectable GWs from collapsars, NDAFs, and GRB jets (internal shocks) are displayed. By combining the electromagnetic counterparts, neutrinos, and GWs, one may constrain the characteristics of collapsars and central BH accretion systems.

To deepen our understanding of the chemical properties of the Planck Galactic Cold Clump (PGCC) G168.72-15.48, we performed observations of nine molecular species, namely, \ce{c-C3H}, \ce{H2CO}, \ce{HC5N}, \ce{HC7N}, \ce{SO}, \ce{CCH}, \ce{N2H+}, \ce{CH3OH}, and \ce{CH3CCH}, toward two dense cores in PGCC G168.72-15.48 using the Tianma Radio Telescope and Purple Mountain Observatory Telescope. We detected \ce{c-C3H}, \ce{H2CO}, \ce{HC5N}, \ce{N2H+}, \ce{CCH}, and \ce{CH3OH} in both G168-H1 and G168-H2 cores, whereas \ce{HC7N} and \ce{CH3CCH} were detected only in G168-H1 and SO was detected only in G168-H2. Mapping observations reveal that the \ce{CCH}, \ce{N2H+}, \ce{CH3OH}, and \ce{CH3CCH} emissions are well coupled with the dust emission in G168-H1. Additionally, \ce{N2H+} exhibits an exceptionally weak emission in the denser and more evolved G168-H2 core, which may be attributed to the \ce{N2H+} depletion. We suggest that the \ce{N2H+} depletion in G168-H2 is dominated by \ce{N2} depletion, rather than the destruction by CO. The local thermodynamic equilibrium calculations indicate that the carbon-chain molecules of \ce{CCH}, \ce{HC5N}, \ce{HC7N}, and \ce{CH3CCH} are more abundant in the younger G168-H1 core. We found that starless core G168-H1 may have the properties of cold dark clouds based on its abundances of carbon-chain molecules. While, the prestellar core G168-H2 exhibits lower carbon-chain molecular abundances than the general cold dark clouds. With our gas-grain astrochemical model calculations, we attribute the observed chemical differences between G168-H1 and G168-H2 to their different gas densities and different evolutionary stages.

We report on the HST detection of the Lyman-continuum (LyC) radiation emitted by a galaxy at redshift z=3.794, dubbed Ion1 (Vanzella et al. 2012). The LyC from Ion1 is detected at rest-frame wavelength 820$\sim$890 \AA with HST WFC3/UVIS in the F410M band ($m_{410}=27.60\pm0.36$ magnitude (AB), peak SNR = 4.17 in a circular aperture with radius r = 0.12'') and at 700$\sim$830 \AA with the VLT/VIMOS in the U-band ($m_U = 27.84\pm0.19$ magnitude (AB), peak SNR = 6.7 with a r = 0.6'' aperture). A 20-hr VLT/VIMOS spectrum shows low- and high-ionization interstellar metal absorption lines, the P-Cygni profile of CIV and Ly$\alpha$ in absorption. The latter spectral feature differs from what observed in known LyC emitters, which show strong Ly$\alpha$ emission. An HST far-UV color map reveals that the LyC emission escapes from a region of the galaxy that is bluer than the rest, presumably because of lower dust obscuration. The F410M image shows that the centroid of the LyC emission is offset from the centroid of the non-ionizing UV emission by 0.12''$\pm$0.03'', corresponding to 0.85$\pm$0.21 kpc (physical), and that its morphology is likely moderately resolved. These morphological characteristics favor a scenario where the LyC photons produced by massive stars escape from low HI column-density "cavities" in the ISM, possibly carved by stellar winds and/or supernova. We also collect the VIMOS U-band images of a sample of 107 Lyman-break galaxies with spectroscopic redshifts at $3.40<z<3.95$, i.e. sampling the LyC, and stack them with inverse-variance weights. No LyC emission is detected in the stacked image, resulting in a 32.5 magnitude (AB) flux limit (1$\sigma$) and an upper limit of absolute LyC escape fraction $f_{esc}^{abs} < 0.63\%$. LyC emitters like Ion1 are very likely at the bright-end of the LyC luminosity function.

The first detected gravitational wave GW170817 from a binary neutron star merger is associated with an important optical transient AT 2017gfo, which is a direct observation of kilonova. Recent observations suggest that the remnant compact object of the binary neutron star merger associated with GW170817/GRB 170817A may be a stable long-lived magnetized neutron star. In this situation, there would be a pulsar wind nebula (PWN) embedded inside the dynamic ejecta. The PWN emission may be absorbed by the ejecta or leak out of the system. We study the effect of the PWN emission on the observed light curves and radiation spectra. Different from previous works, the absorption and leakage of the PWN emission are all involved in our model, where the absorption of the PWN emission heats up the ejecta and alters its radiation. It is found that the characteristic emission of the embedded PWN quickly evolves. For the multiband and long-term observations of AT 2017gfo, we find that the dynamic ejecta with a PWN emission can fit the observational data very well, especially for the light curves at $t\sim 5$ days and those in the late phase. In addition, our model can naturally generate the thermal to nonthermal spectrum evolution of AT 2017gfo. Our fitting result suggests that a PWN is embedded in the AT 2017gfo.

The Askaryan Radio Array is an ultra-high energy (UHE, $>10^{17}$\,eV) neutrino detector designed to observe neutrinos by searching for the radio waves emitted by the relativistic byproducts of neutrino-nucleon interactions in Antarctic ice. In this paper, we present constraints on the diffuse flux of ultra-high energy neutrinos between $10^{16}-10^{21}$\,eV resulting from a search for neutrinos in two complementary analyses, both analyzing four years of data (2013-2016) from the two deep stations (A2, A3) operating at that time. We place a 90\% CL upper limit on the diffuse neutrino flux at $10^{18}$ eV of $EF(E)=3.9\times10^{-16}$~$\textrm{cm}^2$$\textrm{s}^{-1}$$\textrm{sr}^{-1}$. The analysis leverages more than quadruple the livetime of the previous ARA limit, and at this energy is the most sensitive reported by an in-ice radio neutrino detector by more than a factor of four. Looking forward, ARA has yet to analyze more than half of its archived data, and expects that dataset to double in size in the next three years as livetime accumulates in the full five-station array. With this additional livetime, we project ARA will have world-leading sensitivity to neutrinos above $10^{17}$ eV by 2022.

We present a pilot HI survey of 17 Planck Galactic Cold Clumps (PGCCs) with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). HI Narrow Self-Absorption (HINSA) is an effective method to detect cold HI being mixed with molecular hydrogen H$_2$ and improves our understanding of the atomic to molecular transition in the interstellar medium. HINSA was found in 58\% PGCCs that we observed. The column density of HINSA was found to have an intermediate correlation with that of $^{13}$CO, following $\rm log( N(HINSA)) = (0.52\pm 0.26) log(N_{^{13}CO}) + (10 \pm 4.1) $. HI abundance relative to total hydrogen [HI]/[H] has an average value of $4.4\times 10^{-3}$, which is about 2.8 times of the average value of previous HINSA surveys toward molecular clouds. For clouds with total column density N$\rm_H >5 \times 10^{20}$ cm$^{-2}$, an inverse correlation between HINSA abundance and total hydrogen column density is found, confirming the depletion of cold HI gas during molecular gas formation in more massive clouds. Nonthermal line width of $^{13}$CO is about 0-0.5 km s$^{-1}$ larger than that of HINSA. One possible explanation of narrower nonthermal width of HINSA is that HINSA region is smaller than that of $^{13}$CO. Based on an analytic model of H$_2$ formation and H$_2$ dissociation by cosmic ray, we found the cloud ages to be within 10$^{6.7}$-10$^{7.0}$ yr for five sources.

Cosmological applications of HII galaxies (HIIGx) and giant extragalactic HII regions (GEHR) to construct the Hubble diagram at higher redshifts require knowledge of the "$L$--$\sigma$" relation of the standard candles used. In this paper, we study the properties of a large sample of 156 sources (25 high-$z$ HII galaxies, 107 local HII galaxies, and 24 giant extragalactic HII regions) compiled by Terlevich et al.(2015). Using the the cosmological distances reconstructed through two new cosmology-independent methods, we investigate the correlation between the H$\beta$ emission-line luminosity $L$ and ionized-gas velocity dispersion $\sigma$. The method is based on non-parametric reconstruction using the measurements of Hubble parameters from cosmic clocks, as well as the simulated data of gravitational waves from the third-generation gravitational wave detector (the Einstein Telescope, ET), which can be considered as standard sirens. Assuming the emission-line luminosity versus ionized gas velocity dispersion relation, $\log L ($H$\beta) = \alpha \log \sigma($H$\beta)+\kappa$, we find the full sample provides a tight constraint on the correlation parameters. However, similar analysis done on three different sub-samples seems to support the scheme of treating HII galaxies and giant extragalactic HII regions with distinct strategies. Using the corrected "$L$--$\sigma$" relation for the HII observational sample beyond the current reach of Type Ia supernovae, we obtain a value of the matter density parameter, $\Omega_{m}=0.314\pm0.054$ (calibrated with standard clocks) and $\Omega_{m}=0.311\pm0.049$ (calibrated with standard sirens), in the spatially flat $\Lambda$CDM cosmology.

In this paper, we present a scheme to investigate the opacity of the Universe in a cosmological-model-independent way, with the combination of current and future measurements of type Ia supernova sample and galactic-scale strong gravitational lensing systems with SNe Ia acting as background sources. The observational data include the current newly-compiled SNe Ia data (Pantheon sample) and simulated sample of SNe Ia observed by the forthcoming Large Synoptic Survey Telescope (LSST) survey, which are taken for luminosity distances ($D_L$) possibly affected by the cosmic opacity, as well as strongly lensed SNe Ia observed by the LSST, which are responsible for providing the observed time-delay distance ($D_{\Delta t}$) unaffected by the cosmic opacity. Two parameterizations, $\tau(z)=2\beta z$ and $\tau(z)=(1+z)^{2\beta}-1$ are adopted for the optical depth associated to the cosmic absorption. Focusing on only one specific type of standard cosmological probe, this provides an original method to measure cosmic opacity at high precision. Working on the simulated sample of strongly lensed SNe Ia observed by the LSST in 10 year $z$-band search, our results show that, with the combination of the current newly-compiled SNe Ia data (Pantheon sample), there is no significant deviation from the transparency of the Universe at the current observational data level. Moreover, strongly lensed SNe Ia in a 10 year LSST $z$-band search would produce more robust constraints on the validity of cosmic transparency (at the precision of $\Delta\beta=10^{-2}$), with a larger sample of unlensed SNe Ia detected in future LSST survey. We have also discussed the ways in which our methodology could be improved, with the combination of current and future available data in gravitational wave (GW) and electromagnetic (EM) domain.

This paper was withdrawed from the ApJ after the comments from the referee, please be carefully.

Strongly gravitational lensing systems (SGL) encodes cosmology information in source/lens distance ratios $\mathcal{D}_{\rm obs}=\mathcal{D}_{\rm ls}/\mathcal{D}_{\rm s}$, which can be used to precisely constrain cosmological parameters. In this paper, based on future measurements of 390 strong lensing systems from the forthcoming LSST survey, we have successfully reconstructed the distance ratio $\mathcal{D}_{\rm obs}$ (with the source redshift $z_s\sim 4.0$), directly from the data without assuming any parametric form. A recently developed method based on model-independent reconstruction approach, Gaussian Processes (GP) is used in our study of these strong lensing systems. Our results show that independent measurement of the matter density parameter ($\Omega_m$) could be expected from such strong lensing statistics. More specifically, one can expect $\Omega_m$ to be estimated at the precision of $\Delta\Omega_m\sim0.015$ in the concordance $\Lambda$CDM model, which provides comparable constraint on $\Omega_m$ with Planck 2015 results. In the framework of the modified gravity theory (DGP), 390 detectable galactic lenses from future LSST survey would lead to stringent fits of $\Delta\Omega_m\sim0.030$. \textbf{Finally, we have discussed three possible sources of systematic errors (sample incompleteness, the determination of length of lens redshift bin, and the choice of lens redshift shells), and quantified their effects on the final cosmological constraints. Our results strongly indicate that future strong lensing surveys, with the accumulation of a larger and more accurate sample of detectable galactic lenses, will considerably benefit from the methodology described in this analysis. }

In this paper, we study the gravitational waveforms, polarizations and radiation powers of the first relativistic triple systems PSR J0337 + 1715, observed in 2014, by using the post-Newtonian approximations to their lowest order. Although they cannot be observed either by current or next generation of the detectors, they do provide useful information to test different theories of gravity. In particular, we carry out the studies in three different theories, general relativity (GR), Brans-Dicke (BD) gravity, and Einstein-aether ($\ae$) theory. The tensor modes $h_{+}$ and $h_{\times}$ exist in all three theories and have almost equal amplitudes. Their frequencies are all peaked at two locations, $ f^{(+, \times)}_1 \simeq 0.0686656 \mu$Hz and $ f^{(+, \times)}_2 \simeq 14.2138 \mu$Hz, which are about twice of the outer and inner orbital frequencies of the triple system, as predicted in GR. In $\ae$-theory, all the six polarization modes are different from zero, but the breathing $h_b$ and longitudinal $h_L$ modes are not independent and also peaked at two frequencies, but at the frequencies, $ f^{(b, L, \ae)}_{1} \simeq 0.0457771 \mu$Hz and $f^{(b, L, \ae)}_{ 2} \simeq 7.09545 \mu$Hz. A similar phenomenon is also observed in BD gravity, in which only the three modes $h_{+},\; h_{\times}$ and $h_{b}$ exit. We also study the radiation powers, and find that the quadrupole emission in each of the three theories has almost the same amplitude, but the dipole emission can be as big as the quadrupole emission in $\ae$-theory. This can provide a very promising window to obtain severe constraints on $\ae$-theory by the multi-band gravitational wave astronomy.

We present a systematic search for periodically varying quasar and supermassive black hole binary (SMBHB) candidates in the Pan-STARRS1 Medium Deep Survey. From $\sim9,000$ color-selected quasars in a $\sim50$ deg$^{2}$ sky area, we initially identify $26$ candidates with more than $1.5$ cycles of variation. We extend the baseline of observations via our imaging campaign with the Discovery Channel Telescope and the Las Cumbres Observatory network and reevaluate the candidates using a more rigorous, maximum likelihood method. Using a range of statistical criteria and assuming the Damped Random Walk model for normal quasar variability, we identify one statistically significant periodic candidate. We also investigate the capabilities of detecting SMBHBs by the Large Synoptic Survey Telescope using our study with MDS as a benchmark and explore any complementary, multiwavelength evidence for SMBHBs in our sample.

(Abridged) ASHES has been designed to systematically characterize the earliest stages and to constrain theories of high-mass star formation. A total of 12 massive (>500 $M_{\odot}$), cold (<15 K), 3.6-70 $\mu$m dark prestellar clump candidates, embedded in IRDCs, were carefully selected in the pilot survey to be observed with ALMA. We mosaiced each clump (~1 arcmin^2) in dust and line emission with the 12m/7m/TP arrays at 224 GHz, resulting in ~1.2" resolution (~4800 AU). As the first paper of the series, we concentrate on the dust emission to reveal the clump fragmentation. We detect 294 cores, from which 84 (29%) are categorized as protostellar based on outflow activity or 'warm core' line emission. The remaining 210 (71%) are considered prestellar core candidates. The number of detected cores is independent of the mass sensitivity range of the observations. On average, more massive clumps tend to form more cores. We find a large population of low-mass (<1 M) cores and no high-mass (>30 $M_{\odot}$) prestellar cores. The most massive prestellar core has a mass of 11 $M_{\odot}$. From the prestellar CMF, we derive a power law index of 1.17+-0.1, slightly shallower than Salpeter (1.35). We use the MST technique to characterize the separation between cores and their spatial distribution, and derive mass segregation ratios. While there is a range of core masses and separations detected in the sample, the mean separation and mass of cores are well explained by thermal fragmentation and are inconsistent with turbulent Jeans fragmentation. The core spatial distribution is well described by hierarchical subclustering rather than centrally peaked clustering. There is no conclusive evidence of mass segregation. We test several theoretical conditions, and conclude that overall, competitive accretion and global hierarchical collapse scenarios are favored over the turbulent core accretion scenario.

We started a multi-scale analysis of G202.3+2.5, an intertwined filamentary region of Monoceros OB1. In Paper I, we examined the distributions of dense cores and protostars and found enhanced star formation (SF) activity in the junction region of the filaments. In this second paper, we aim to unveil the connections between the core and filament evolutions, and between the filament dynamics and the global evolution of the cloud. We characterise the gas dynamics and energy balance using Herschel and WISE observations and molecular tracers observed with the IRAM 30m and TRAO 14m telescopes. The velocity field of the cloud is examined and velocity-coherent structures are put in perspective with the cloud environment. Two main velocity components (VCs) are revealed, well separated in the north and merged around the location of intense N2H+ emission where Paper I found the peak of SF activity. The relative position of the two VCs along the sightline, and the velocity gradient in N2H+ emission imply that the VCs have been undergoing collision for ~10^5 yrs. The dense gas where N2H+ is detected is interpreted as the compressed region between the two filaments, which corresponds to a high mass inflow rate of ~1e-3 Msun/yr and possibly leads to an increase in its SF efficiency. We identify a protostar in the junction region that possibly powers two crossed intermittent outflows. We show that the HII region around the nearby cluster NCG 2264 is still expanding and its role in the collision is examined. However, we cannot rule out the idea that the collision arises mostly from the global collapse of the cloud. The (sub-)filament-scale observables examined in this paper reveal a collision between G202.3+2.5 sub-structures and its probable role in feeding the cores in the junction region. One must now characterise the cloud morphology, its fragmentation, and magnetic field, all at high resolution.

RNO is the mid-scale discovery instrument designed to make the first observation of neutrinos from the cosmos at extreme energies, with sensitivity well beyond current instrument capabilities. This new observatory will be the largest ground-based neutrino telescope to date, enabling the measurement of neutrinos above $10^{16}$ eV, determining the nature of the astrophysical neutrino flux that has been measured by IceCube at higher energies, similarly extending the reach of multi-messenger astrophysics to the highest energies, and enabling investigations of fundamental physics at energies unreachable by particle accelerators on Earth.

Current theories and models attempt to explain star formation globally, from core scales to giant molecular cloud scales. A multi-scale observational characterisation of an entire molecular complex is necessary to constrain them. We investigate star formation in G202.3+2.5, a ~10x3 pc sub-region of the Monoceros OB1 cloud with a complex morphology harbouring interconnected filamentary structures. We aim to connect the evolution of cores and filaments in G202.3+2.5 with the global evolution of the cloud and to identify the engines of the cloud dynamics. In this first paper, the star formation activity is evaluated by surveying the distributions of dense cores and protostars, and their evolutionary state, as characterised using both infrared observations from the Herschel and WISE telescopes and molecular line observations with the IRAM 30-m telescope. We find ongoing star formation in the whole cloud, with a local peak in star formation activity around the centre of G202.3+2.5 where a chain of massive cores (10-50 Msun) forms a massive ridge (>150 Msun). All evolutionary stages from starless cores to Class II protostars are found in G202.3+2.5, including a possibly starless, large column density (8x10^{22} cm-2), and massive (52 Msun) core. All the core-scale observables examined in this paper point to an enhanced star formation activity centred on the junction between the three main branches of the ramified structure of G202.3+2.5. This suggests that the increased star-formation activity results from the convergence of these branches. To further investigate the origin of this enhancement, it is now necessary to extend the analysis to larger scales, in order to examine the relationship between cores, filaments and their environment. We address these points through the analysis of the dynamics of G202.3+2.5 in a joint paper.

We report on analysis of englacial radio-frequency (RF) pulser data received over horizontal baselines of 1--5 km, based on broadcasts from two sets of transmitters deployed to depths of up to 1500 meters at the South Pole. First, we analyze data collected usingtwo RF bicone transmitters 1400 meters below the ice surface, and frozen into boreholes drilled for the IceCube experiment in 2011. Additionally, in Dec., 2018, a fat-dipole antenna, fed by one of three high-voltage (~1 kV), fast (~(1-5 ns)) signal generators was lowered into the 1700-m deep icehole drilled for the South Pole Ice Core Experiment (SPICE), approximately 3 km from the geographic South Pole. Signals from transmitters were recorded on the five englacial multi-receiver ARA stations, with receiver depths between 60--200 m. We confirm the long, >1 km RF electric field attenuation length, test our observed signal arrival timing distributions against models, and measure birefringent asymmetries at the 0.15% level.

Recently it was shown that cosmic anisotropy can be well tested using either standard siren measurement of luminosity distance $d_\mathrm{L}(z)$ from gravitational-wave (GW) observation or dispersion measure ($\mathrm{DM}(z)$) from fast radio burst (FRB). It was also observed that the combined measurement of $d_\mathrm{L}(z)\cdot\mathrm{DM}(z)$ from the GW/FRB association system as suggested in some of FRB models is more effective to constrain cosmological parameters than $d_\mathrm{L}(z)$ or $\mathrm{DM}(z)$ separately due to its independence from Hubble constant. In this paper, we will show both theoretically and with simulation that, this upgraded sirens from combined GW/FRB observations could test cosmic anisotropy with a double relative sensitivity compared to the usual standard siren from GW observation alone.

We present the highest resolution single-dish submillimetre observations of the detached shell source U Antliae to date. The observations were obtained at $450~\micron$ and $850~\micron$ with SCUBA-2 instrument on the James Clerk Maxwell Telescope as part of the Nearby Evolved Stars Survey. The emission at $850~\micron$ peaks at $40\arcsec$ with hints of a second peak seen at $\sim 20\arcsec$. The emission can be traced out to a radius of $56\arcsec$ at a $3\sigma$ level. The outer peak observed at $850~\micron$ aligns well with the peak observed at Herschel/PACS wavelengths. With the help of spectral energy distribution fitting and radiative transfer calculations of multiple-shell models for the circumstellar envelope, we explore the various shell structures and the variation of grain sizes along the in the circumstellar envelope. We determine a total shell dust mass of $(2.0 \pm 0.3) \times 10^{-5}$ M$_{\odot}$ and established that the thermal pulse which gave rise to the detached shell occurred 3500 $\pm$ 500 years ago.

The Lorentz factor of a relativistic jet and its evolution during the jet expansion are difficult to estimate, especially for the jets in gamma-ray bursts (GRBs). However, it is related to the understanding of jet physics. Owing to the absorption of two-photon pair production ($\gamma\gamma{\leftrightarrow}e^+e^-$), a high-energy spectral cutoff may appear in the radiation spectrum of GRBs. We search such kind of high-energy cutoff in GRB 160625B, which is one of the brightest bursts in recent years. It is found that the high-energy spectral cutoff is obvious for the first pulse in the second emission episode of GRB 160625B (i.e., $\sim$186-192 s after the burst first trigger), which is smooth and well-shaped. Then, we estimate the Lorentz factor and radiation location of the jet shell associated with the first pulse in the second emission episode of GRB 160625B. It is found that the radiation location increases with time. In addition, the Lorentz factor remains almost constant during the expansion of the jet shell. This reveals that the magnetization of the jet is low or intermediate in the emission region, event though the jet could be still Poynting flux dominated at smaller radii to avoid a bright thermal component in the emission episode.

A Monte Carlo simulation program for the radio detection of Ultra High Energy (UHE) neutrino interactions in the Antarctic ice as viewed by the Antarctic Impulsive Transient Antenna (ANITA) is described in this article. The program, icemc, provides an input spectrum of UHE neutrinos, the parametrization of the Askaryan radiation generated by their interaction in the ice, and the propagation of the radiation through ice and air to a simulated model of the third and fourth ANITA flights. This paper provides an overview of the icemc simulation, descriptions of the physics models used and of the ANITA electronics processing chain, data/simulation comparisons to validate the predicted performance, and a summary of the impact of published results.

We present the B-fields mapped in IRDC G34.43+0.24 using 850\,$\mu$m polarized dust emission observed with the POL-2 instrument at JCMT. We examine the magnetic field geometries and strengths in the northern, central, and southern regions of the filament. The overall field geometry is ordered and aligned closely perpendicular to the filament's main axis, particularly in regions containing the central clumps MM1 and MM2, whereas MM3 in the north has field orientations aligned with its major axis. The overall field orientations are uniform at large (POL-2 at 14$\arcsec$ and SHARP at 10$\arcsec$) to small scales (TADPOL at 2.5$\arcsec$ and SMA at 1.5$\arcsec$) in the MM1 and MM2 regions. SHARP/CSO observations in MM3 at 350\,$\mu$m from Tang et al. show a similar trend as seen in our POL-2 observations. TADPOL observations demonstrate a well-defined field geometry in MM1/MM2 consistent with MHD simulations of accreting filaments. We obtained a plane-of-sky magnetic field strength of 470$\pm$190\,$\mu$G, 100$\pm$40\,$\mu$G, and 60$\pm$34\,$\mu$G in the central, northern and southern regions of G34, respectively, using the updated Davis-Chandrasekhar-Fermi relation. The estimated value of field strength, combined with column density and velocity dispersion values available in the literature, suggests G34 to be marginally critical with criticality parameter $\rm \lambda$ values 0.8$\pm$0.4, 1.1$\pm$0.8, and 0.9$\pm$0.5 in the central, northern, and southern regions, respectively. The turbulent motions in G34 are sub-Alfv\'{e}nic with Alfv\'{e}nic Mach numbers of 0.34$\pm$0.13, 0.53$\pm$0.30, and 0.49$\pm$0.26 in the three regions. The observed aligned B-fields in G34.43+0.24 are consistent with theoretical models suggesting that B-fields play an important role in guiding the contraction of the cloud driven by gravity.

L1521F is found to be forming multiple cores and it is cited as an example of the densest core with an embedded VeLLO in a highly dynamical environment. We present the core-scale magnetic fields (B-fields) in the near vicinity of the VeLLO L1521F-IRS using submm polarization measurements at 850$~\mu$m using JCMT POL-2. This is the first attempt to use high-sensitivity observations to map the sub-parsec scale B-fields in a core with a VeLLO. The B-fields are ordered and very well connected to the parsec-scale field geometry seen in our earlier optical polarization observations and the large-scale structure seen in Planck dust polarization. The core scale B-field strength estimated using Davis-Chandrasekhar-Fermi relation is $\rm 330\pm100~\mu$G which is more than ten times of the value we obtained in the envelope (envelope in this paper is "core envelope"). This indicates that B-fields are getting stronger on smaller scales. The magnetic energies are found to be 1 to 2 orders of magnitude higher than non-thermal kinetic energies in the envelope and core. This suggests that magnetic fields are more important than turbulence in the energy budget of L1521F. The mass-to-flux ratio of 2.3$\pm$0.7 suggests that the core is magnetically-supercritical. The degree of polarization is steadily decreasing towards the denser part of the core with a power law slope of -0.86.

The Askaryan Radio Array (ARA) is an ultra-high energy (UHE) neutrino telescope at the South Pole consisting of an array of radio antennas aimed at detecting the Askaryan radiation produced by neutrino interactions in the ice. Currently, the experiment has five stations in operation that have been deployed in stages since 2012. This contribution focuses on the development of a search for a diffuse flux of neutrinos in two ARA stations (A2 and A3) from 2013-2016. A background of $\sim 0.01-0.02$ events is expected in one station in each of two search channels in horizontal- and vertical-polarizations. The expected new constraints on the flux of ultra-high energy neutrinos based on four years of analysis with two stations improve on the previous limits set by ARA by a factor of about two. The projected sensitivity of ARA's five-station dataset is beginning to be competitive with other neutrino telescopes at high energies near $10^{10.5}\,$GeV.

This white paper gives a brief summary of Galactic continuum surveys with the next generation 850 micron camera at the James Clerk Maxwell Telescope (JCMT) in the next decade. This new camera will have mapping speeds at least 10-20 times faster than the present SCUBA-2 camera, and will enable deep (<10 mJy/beam) and extensive continuum surveys such as a wider (-5 deg.<l<240 deg., |b|<2-5 deg.) Galactic Plane survey, a larger Gould Belt survey, and a follow-up survey of all Planck compact objects visible from the northern hemisphere. These surveys will provide a complete census of molecular clouds, filaments, and dense cores across the Galaxy, vital for studying star formation in various environments.

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

The gas depletion factor $\gamma(z)$, i.e., the average ratio of the gas mass fraction to the cosmic mean baryon fraction of galaxy clusters, plays a very important role in the cosmological application of the gas mass fraction measurements. In this paper, using the newest catalog of 182 galaxy clusters detected by the Atacama Cosmology Telescope (ACT) Polarization experiment, we investigate the possible redshift evolution of $\gamma(z)$ through a new cosmology-independent method. The method is based on non-parametric reconstruction using the measurements of Hubble parameters from cosmic chronometers. Unlike hydrodynamical simulations suggesting constant depletion factor, our results reveal the trend of $\gamma(z)$ decreasing with redshift. This result is supported by a parametric model fit as well as by calculations on the reduced ACTPol sample and on the alternative sample of 91 SZ clusters reported earlier in ACT compilation. Discussion of possible systematic effects leaves an open question about validity of the empirical relation $M_{tot}$-$f_{gas}$ obtained on very close clusters. These results might pave the way to explore the hot gas fraction within large radii of galaxy clusters as well as its possible evolution with redshift, which should be studied further on larger galaxy cluster samples in the upcoming X-ray/SZ cluster surveys.

We analyze extensive spectroscopic and photometric data of the hypervariable quasar SDSS J131424+530527 (RMID 017) at z=0.456, an optical "changing look" quasar from the Sloan Digital Sky Survey Reverberation Mapping project that increased in optical luminosity by a factor of 10 between 2014 and 2017. The observed broad emission lines all respond in luminosity and width to the changing optical continuum, as expected for photoionization in a stratified, virialized broad emission line region. The luminosity changes therefore result from intrinsic changes in accretion power rather than variable obscuration. The variability is continuous and apparently stochastic, disfavoring an origin as a discrete event such as a tidal disruption flare or microlensing event. It is coordinated on day timescales with blue leading red, consistent with reprocessing powering the entire optical SED. We show that this process cannot work in a standard thin disk geometry on energetic grounds, and would instead require a large covering factor reprocessor. Disk instability models could potentially also explain the data, provided that the instability sets in near the inner radius of a geometrically thick accretion disk.

In the early stages of star formation, a protostar is deeply embedded in an optically thick envelope such that it is not directly observable. Variations in the protostellar accretion rate, however, will cause luminosity changes that are reprocessed by the surrounding envelope and are observable at submillimeter wavelengths. We searched for submillimeter flux variability toward 12 Planck Galactic Cold Clumps detected by the James Clerk Maxwell Telescope (JCMT)-SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE) survey. These observations were conducted at 850 um using the JCMT/SCUBA-2. Each field was observed three times over about 14 months between 2016 April and 2017 June. We applied a relative flux calibration and achieved a calibration uncertainty of ~ 3.6% on average. We identified 136 clumps across 12 fields and detected four sources with flux variations of ~ 30%. For three of these sources, the variations appear to be primarily due to large-scale contamination, leaving one plausible candidate. The flux change of the candidate may be associated with low- or intermediate-mass star formation assuming a distance of 1.5 kpc, although we cannot completely rule out the possibility that it is a random deviation. Further studies with dedicated monitoring would provide a better understanding of the detailed relationship between submillimeter flux and accretion rate variabilities while enhancing the search for variability in star-forming clumps farther away than the Gould Belt.

The dependence of polarization fraction $p$ on total intensity $I$ in polarized submillimeter emission measurements is typically parameterized as $p\propto I^{-\alpha}$ $(\alpha \leq 1)$, and used to infer dust grain alignment efficiency in star-forming regions, with an index $\alpha=1$ indicating near-total lack of alignment of grains with the magnetic field. In this work we demonstrate that the non-Gaussian noise characteristics of polarization fraction may produce apparent measurements of $\alpha \sim 1$ even in data with significant signal-to-noise in Stokes $Q$, $U$ and $I$ emission, and so with robust measurements of polarization angle. We present a simple model demonstrating this behavior, and propose a criterion by which well-characterized measurements of polarization fraction may be identified. We demonstrate that where our model is applicable, $\alpha$ can be recovered by fitting the $p-I$ relationship with the mean of the Rice distribution, without statistical debiasing of polarization fraction. We apply our model to JCMT BISTRO Survey POL-2 850$\mu$m observations of three clumps in the Ophiuchus Molecular Cloud, finding that in the externally-illuminated Oph A region, $\alpha\approx 0.34$, while in the more isolated Oph B and C, despite their differing star formation histories, $\alpha \sim 0.6-0.7$. Our results thus suggest that dust grain alignment in dense gas is more strongly influenced by incident interstellar radiation field than by star formation history. We further find that grains may remain aligned with the magnetic field at significantly higher gas densities than has previously been believed, thus allowing investigation of magnetic field properties within star-forming clumps and cores.

In the context of the upcoming SRG/eROSITA survey, we present an N-body simulation-based mock catalogue for X-ray selected AGN samples. The model reproduces the observed hard X-ray AGN luminosity function (XLF) and the soft X-ray logN-logS from redshift 0 to 6. The XLF is reproduced to within $\pm5\%$ and the logN-logS to within $\pm20\%$. We develop a joint X-ray -- optical extinction and classification model. We adopt a set of empirical spectral energy distributions to predict observed magnitudes in the UV, optical and NIR. With the latest eROSITA all sky survey sensitivity model, we create a high-fidelity full-sky mock catalogue of X-ray AGN. It predicts their distributions in right ascension, declination, redshift and fluxes. Using empirical medium resolution optical spectral templates and an exposure time calculator, we find that $1.1\times10^6$ ($4\times10^5$) fiber-hours are needed to follow-up spectroscopically from the ground the detected X-ray AGN with an optical magnitude $21<r<22.8$ ($22.8<r<25$) with a 4-m (8-m) class multi-object spectroscopic facility. We find that future clustering studies will measure the AGN bias to the percent level at redshift $z<1.2$ and should discriminate possible scenarios of galaxy-AGN co-evolution. We predict the accuracy to which the baryon acoustic oscillation standard ruler will be measured using X-ray AGN: better than 3\% for AGN between redshift 0.5 to 3 and better than 1\% using the Ly$\alpha$ forest of X-ray QSOs discovered between redshift 2 and 3. eROSITA will provide an outstanding set of targets for future galaxy evolution and cosmological studies.

We present a multi-wavelength study of the Planck cold clump G181.84+0.31, which is located at the northern end of the extended filamentary structure S242. We have extracted 9 compact dense cores from the SCUBA-2 850 um map, and we have identified 18 young stellar objects (YSOs, 4 Class I and 14 Class II) based on their Spitzer, Wide-field Infrared Survey Explorer (WISE) and Two-Micron All-Sky Survey (2MASS) near- and mid-infrared colours. The dense cores and YSOs are mainly distributed along the filamentary structures of G181.84 and are well traced by HCO$^{+}$(1-0) and N$_{2}$H$^{+}$(1-0) spectral-line emission. We find signatures of sequential star formation activities in G181.84: dense cores and YSOs located in the northern and southern sub-structures are younger than those in the central region. We also detect global velocity gradients of about 0.8$\pm$0.05 km s$^{-1}$pc$^{-1}$ and 1.0$\pm$0.05 km s$^{-1}$pc$^{-1}$ along the northern and southern sub-structures, respectively, and local velocity gradients of 1.2$\pm$0.1 km s$^{-1}$pc$^{-1}$ in the central substructure. These results may be due to the fact that the global collapse of the extended filamentary structure S242 is driven by an edge effect, for which the filament edges collapse first and then further trigger star formation activities inward. We identify three substructures in G181.84 and estimate their critical masses per unit length, which are $\sim$ 101$\pm$15 M$_{\odot}$ pc$^{-1}$, 56$\pm$8 M$_{\odot}$ pc$^{-1}$ and 28$\pm$4 M$_{\odot}$ pc$^{-1}$, respectively. These values are all lower than the observed values ($\sim$ 200 M$_{\odot}$ pc$^{-1}$), suggesting that these sub-structures are gravitationally unstable.

Using the new equipment of the Shanghai Tian Ma Radio Telescope, we have searched for carbon-chain molecules (CCMs) towards five outflow sources and six Lupus I starless dust cores, including one region known to be characterized by warm carbon-chain chemistry (WCCC), Lupus I-1 (IRAS 15398-3359), and one TMC-1 like cloud, Lupus I-6 (Lupus-1A). Lines of HC3N J=2-1, HC5N J=6-5, HC7N J=14-13, 15-14, 16-15 and C3S J=3-2 were detected in all the targets except in the outflow source L1660 and the starless dust core Lupus I-3/4. The column densities of nitrogen-bearing species range from 10$^{12}$ to 10$^{14}$ cm$^{-2}$ and those of C$_3$S are about 10$^{12}$ cm$^{-2}$. Two outflow sources, I20582+7724 and L1221, could be identified as new carbon-chain--producing regions. Four of the Lupus I dust cores are newly identified as early quiescent and dark carbon-chain--producing regions similar to Lup I-6, which together with the WCCC source, Lup I-1, indicate that carbon-chain-producing regions are popular in Lupus I which can be regard as a Taurus like molecular cloud complex in our Galaxy. The column densities of C3S are larger than those of HC7N in the three outflow sources I20582, L1221 and L1251A. Shocked carbon-chain chemistry (SCCC) is proposed to explain the abnormal high abundances of C3S compared with those of nitrogen-bearing CCMs. Gas-grain chemical models support the idea that shocks can fuel the environment of those sources with enough $S^+$ thus driving the generation of S-bearing CCMs.

We detected carbon-chain molecules (CCMs) HC2n+1N (n=1-3) and C3S in Ku band as well as {high-energy excitation lines including C4H N=9-8, J=17/2-15/2, 19/2-17/2, and CH3CCH J=5-4, K=2 in the 3 mm band toward a starless core called the eastern molecular core (EMC) of L1489 IRS. Maps of all the observed lines were also obtained. Comparisons with a number of early starless cores and WCCC source L1527 show that the column densities of C4H and CH3CCH are close to those of L1527, and the CH3CCH column densities of the EMC and L1527 are slightly higher than those of TMC-1. The EMC and L1527 have similar C3S column densities, but they are much lower than those of all the starless cores, with only 6.5% and 10% of the TMC-1 value, respectively. The emissions of the N-bearing species of the EMC and L1527 are at the medium level of the starless cores. These comparisons show that the CCM emissions in the EMC are similar to those of L1527, though L1527 contains a protostar. Although dark and quiescent, the EMC is warmer and at a later evolutionary stage than classical carbon-chain--producing regions in the cold, dark, quiescent early phase. The PACS, SPIRE, and SCUBA maps evidently show that the L1489 IRS seems to be the heating source of the EMC. Although it is located at the margins of the EMC, its bolometric luminosity and bolometric temperature are relatively high. Above all, the EMC is a rather particular carbon-chain-producing region and is quite significant for CCM science.

We present magnetohydrodynamic simulation of the evolution from quasi-equilibrium to onset of eruption of a twisted, prominence-forming coronal magnetic flux rope underlying a corona streamer. The flux rope is built up by an imposed flux emergence at the lower boundary. During the quasi-static phase of the evolution, we find the formation of a prominence-cavity system with qualitative features resembling observations, as shown by the synthetic SDO/AIA EUV images with the flux rope observed above the limb viewed nearly along its axis. The cavity contains substructures including ``U''-shaped or horn-liked features extending from the prominence enclosing a central ``cavity'' on top of the prominence. The prominence condensations form in the dips of the highly twisted field lines due to runaway radiative cooling and the cavity is formed by the density depleted portions of the prominence-carrying field lines extending up from the dips. The prominence ``horns'' are threaded by twisted field lines containing shallow dips, where the prominence condensations have evaporated to coronal temperatures. The central ``cavity'' enclosed by the horns is found to correspond to a central hot and dense core containing twisted field lines that do not have dips. The flux rope eventually erupts as its central part rises quasi-statically to a critical height consistent with the onset of the torus instability. The erupting flux rope accelerates to a fast speed of nearly 900 km/s and the associated prominence eruption shows significant rotational motion and a kinked morphology.

The structured jet is involved to explain the afterglows and even the prompt emission of GRB 170817A. In this paper, we stress that for a precessing jet, the jet structure in the prompt emission phase and that in the afterglow phase may be different. The jet structure in the afterglow phase can be non-uniform even if a narrow-uniform jet is presented in the prompt emission phase. We estimate the jet structure in the afterglow phase under the situation that a narrow-uniform-precessing jet is launched from the central engine of gamma-ray burst. With different precession angles, it is found that the structured jet can be roughly described as follows: a narrow uniform core with power-law wings and sharp cut-off edges, a Gaussian profile, a ring shape, or other complex profile in energy per solid angle. Correspondingly, the afterglows for our obtained structured jets are also estimated. We find that the estimates of the intrinsic kinetic energy, the electron index, and the jet opening angle based on the afterglows formed in a precessing system may be incorrect. Our obtained structured jet is likely to be revealed by future observations for a fraction of gravitational wave detected merging compact binary systems (e.g., black hole-neutron star mergers).

As the plausible central engine of gamma-ray bursts (GRBs), a black hole (BH) hyperaccretion disk should be in a state of neutrino-dominated accretion flow (NDAF) if the accretion rate is larger than the ignition rate of an NDAF. A rotating stellar-mass BH surrounded by a hyperaccretion disk might be born in the center of a massive core collapsar. In the initial hundreds of seconds of the accretion process, the mass supply rate of the massive progenitor is generally higher than the ignition accretion rate, but the jets are generally choked in the envelope. Thus, neutrinos can be emitted from the center of a core collapsar. In this paper, we study the effects of the masses and metallicities of progenitor stars on the time-integrated spectra of electron neutrinos from NDAFs. The peak energies of the calculated spectra are approximately 10-20 MeV. The mass of a collapsar has little influence on the neutrino spectrum, and a low metallicity is beneficial to the production of low-energy ($\lesssim$ 1 MeV) neutrinos. We also investigate the differences in the electron neutrino spectra between NDAFs and proto-neutron stars. Combining with the electromagnetic counterparts and multi-messenger astronomy, one may verify the possible remnants of the core collapse of massive stars with future neutrino detectors.

"TRAO FUNS" is a project to survey Gould Belt's clouds in molecular lines. This paper presents its first results on the central region of the California molecular cloud, L1478. We performed On-The-Fly mapping observations using the Taedeok Radio Astronomy Observatory (TRAO) 14m single dish telescope equipped with a 16 multi-beam array covering $\sim$1.0 square degree area of this region using C$^{18}$O (1-0) mainly tracing low density cloud and about 460 square arcminute area using N$_{2}$H$^{+}$ (1-0) mainly tracing dense cores. CS (2-1) and SO $(3_{2}-2_{1})$ were also used simultaneously to map $\sim$440 square arcminute area of this region. We identified 10 filaments by applying the dendrogram technique to the C$^{18}$O data-cube and 8 dense N$_{2}$H$^{+}$ cores by using {\sc FellWalker}. Basic physical properties of filaments such as mass, length, width, velocity field, and velocity dispersion are derived. It is found that L1478 consists of several filaments with slightly different velocities. Especially the filaments which are supercritical are found to contain dense cores detected in N$_{2}$H$^{+}$. Comparison of non-thermal velocity dispersions derived from C$^{18}$O and N$_{2}$H$^{+}$ for the filaments and dense cores indicates that some of dense cores share similar kinematics with those of the surrounding filaments while several dense cores have different kinematics with those of their filaments. This suggests that the formation mechanism of dense cores and filaments can be different in individual filaments depending on their morphologies and environments.

We present the POL-2 850 $\mu$m linear polarization map of the Barnard 1 clump in the Perseus molecular cloud complex from the B-fields In STar-forming Region Observations (BISTRO) survey at the James Clerk Maxwell Telescope. We find a trend of decreasing polarization fraction as a function of total intensity, which we link to depolarization effects towards higher density regions of the cloud. We then use the polarization data at 850 $\mu$m to infer the plane-of-sky orientation of the large-scale magnetic field in Barnard 1. This magnetic field runs North-South across most of the cloud, with the exception of B1-c where it turns more East-West. From the dispersion of polarization angles, we calculate a turbulence correlation length of $5.0 \pm 2.5$ arcsec ($1500$ au), and a turbulent-to-total magnetic energy ratio of $0.5 \pm 0.3$ inside the cloud. We combine this turbulent-to-total magnetic energy ratio with observations of NH$_3$ molecular lines from the Green Bank Ammonia Survey (GAS) to estimate the strength of the plane-of-sky component of the magnetic field through the Davis-Chandrasekhar-Fermi method. With a plane-of-sky amplitude of $120 \pm 60$ $\mu$G and a criticality criterion $\lambda_c = 3.0 \pm 1.5$, we find that Barnard 1 is a supercritical molecular cloud with a magnetic field nearly dominated by its turbulent component.

We present a detailed X-ray spectral analysis of 1152 AGNs selected in the Chandra Deep Fields (CDFs), in order to identify highly obscured AGNs ($N_{\rm H} > 10^{23}\ \rm cm^{-2}$). By fitting spectra with physical models, 436 (38%) sources with $L_{\rm X} > 10^{42}\ \rm erg\ s^{-1}$ are confirmed to be highly obscured, including 102 Compton-thick (CT) candidates. We propose a new hardness-ratio measure of the obscuration level which can be used to select highly obscured AGN candidates. The completeness and accuracy of applying this method to our AGNs are 88% and 80%, respectively. The observed logN-logS relation favors cosmic X-ray background models that predict moderate (i.e., between optimistic and pessimistic) CT number counts. 19% (6/31) of our highly obscured AGNs that have optical classifications are labeled as broad-line AGNs, suggesting that, at least for part of the AGN population, the heavy X-ray obscuration is largely a line-of-sight effect, i.e., some high-column-density clouds on various scales (but not necessarily a dust-enshrouded torus) along our sightline may obscure the compact X-ray emitter. After correcting for several observational biases, we obtain the intrinsic NH distribution and its evolution. The CT-to-highly-obscured fraction is roughly 52% and is consistent with no evident redshift evolution. We also perform long-term (~17 years in the observed frame) variability analyses for 31 sources with the largest number of counts available. Among them, 17 sources show flux variabilities: 31% (5/17) are caused by the change of NH, 53% (9/17) are caused by the intrinsic luminosity variability, 6% (1/17) are driven by both effects, and 2 are not classified due to large spectral fitting errors.

We present the 850 $\mu$m polarization observations toward the IC5146 filamentary cloud taken using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) and its associated polarimeter (POL-2), mounted on the James Clerk Maxwell Telescope (JCMT), as part of the B-fields In STar forming Regions Observations (BISTRO). This work is aimed at revealing the magnetic field morphology within a core-scale ($\lesssim 1.0$ pc) hub-filament structure (HFS) located at the end of a parsec-scale filament. To investigate whether or not the observed polarization traces the magnetic field in the HFS, we analyze the dependence between the observed polarization fraction and total intensity using a Bayesian approach with the polarization fraction described by the Rice likelihood function, which can correctly describe the probability density function (PDF) of the observed polarization fraction for low signal-to-noise ratio (SNR) data. We find a power-law dependence between the polarization fraction and total intensity with an index of 0.56 in $A_V\sim$ 20--300 mag regions, suggesting that the dust grains in these dense regions can still be aligned with magnetic fields in the IC5146 regions. Our polarization maps reveal a curved magnetic field, possibly dragged by the contraction along the parsec-scale filament. We further obtain a magnetic field strength of 0.5$\pm$0.2 mG toward the central hub using the Davis-Chandrasekhar-Fermi method, corresponding to a mass-to-flux criticality of $\sim$ $1.3\pm0.4$ and an Alfv\'{e}nic Mach number of $<$0.6. These results suggest that gravity and magnetic field is currently of comparable importance in the HFS, and turbulence is less important.

We aim to reveal the physical properties and chemical composition of the cores in the California molecular cloud (CMC), so as to better understand the initial conditions of star formation. We made a high-resolution column density map (18.2") with Herschel data, and extracted a complete sample of the cores in the CMC with the \textsl{fellwalker} algorithm. We performed new single-pointing observations of molecular lines near 90 GHz with the IRAM 30m telescope along the main filament of the CMC. In addition, we also performed a numerical modeling of chemical evolution for the cores under the physical conditions. We extracted 300 cores, of which 33 are protostellar and 267 are starless cores. About 51\% (137 of 267) of the starless cores are prestellar cores. Three cores have the potential to evolve into high-mass stars. The prestellar core mass function (CMF) can be well fit by a log-normal form. The high-mass end of the prestellar CMF shows a power-law form with an index $\alpha=-0.9\pm 0.1$ that is shallower than that of the Galactic field stellar mass function. Combining the mass transformation efficiency ($\varepsilon$) from the prestellar core to the star of $15\pm 1\%$ and the core formation efficiency (CFE) of 5.5\%, we suggest an overall star formation efficiency of about 1\% in the CMC. In the single-pointing observations with the IRAM 30m telescope, we find that 6 cores show blue-skewed profile, while 4 cores show red-skewed profile. [$\rm {HCO}^{+}$]/[HNC] and [$\rm {HCO}^{+}$]/$\rm [N_{2}H^{+}]$ in protostellar cores are higher than those in prestellar cores; this can be used as chemical clocks. The best-fit chemical age of the cores with line observations is $\sim 5\times 10^4$~years.

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.

Non-relativistic QCD axions or axion-like particles are among the most popular candidates for cold Dark Matter (DM) in the universe. We proposed to detect axion-like DM, using linearly polarized pulsar light as a probe. Because of birefringence effect potentially caused by an oscillating galactic axion DM background, when pulsar light travels across the galaxy, its linear polarization angle may vary with time. With a soliton+NFW galactic DM density profile, we show that this strategy can potentially probe an axion-photon coupling as small as $\sim 10^{-13}$ GeV$^{-1}$ for axion mass $m_a \sim 10^{-22}-10^{-20}$ eV, given the current measurement accuracy. An exclusion limit stronger than CAST ($ \sim 10^{-10}$ GeV$^{-1}$) and SN1987A ($ \sim 10^{-11}$ GeV$^{-1}$) could be extended up to $m_a \sim 10^{-18}$ eV and $\sim 10^{-19}$ eV, respectively.

We report 850~$\mu$m dust polarization observations of a low-mass ($\sim$12 $M_{\odot}$) starless core in the $\rho$ Ophiuchus cloud, Ophiuchus C, made with the POL-2 instrument on the James Clerk Maxwell Telescope (JCMT) as part of the JCMT B-fields In STar-forming Region Observations (BISTRO) survey. We detect an ordered magnetic field projected on the plane of sky in the starless core. The magnetic field across the $\sim$0.1~pc core shows a predominant northeast-southwest orientation centering between $\sim$40$^\circ$ to $\sim$100$^\circ$, indicating that the field in the core is well aligned with the magnetic field in lower-density regions of the cloud probed by near-infrared observations and also the cloud-scale magnetic field traced by Planck observations. The polarization percentage ($P$) decreases with an increasing total intensity ($I$) with a power-law index of $-$1.03 $\pm$ 0.05. We estimate the plane-of-sky field strength ($B_{\mathrm{pos}}$) using modified Davis-Chandrasekhar-Fermi (DCF) methods based on structure function (SF), auto-correlation (ACF), and unsharp masking (UM) analyses. We find that the estimates from the SF, ACF, and UM methods yield strengths of 103 $\pm$ 46 $\mu$G, 136 $\pm$ 69 $\mu$G, and 213 $\pm$ 115 $\mu$G, respectively. Our calculations suggest that the Ophiuchus C core is near magnetically critical or slightly magnetically supercritical (i.e. unstable to collapse). The total magnetic energy calculated from the SF method is comparable to the turbulent energy in Ophiuchus C, while the ACF method and the UM method only set upper limits for the total magnetic energy because of large uncertainties.

We present a multiwavelength study of the extended green object, G12.42+0.50 in this paper. The associated ionized, dust, and molecular components of this source are studied in detail employing various observations at near-, mid- and far-infrared, submillimeter and radio wavelengths. Radio continuum emission mapped at 610 and 1390 MHz, using the Giant Meterwave Radio Telescope, India, advocates for a scenario of coexistence of an UC H II region and an ionized thermal jet possibly powered by the massive young stellar object, IRAS 18079-1756 with an estimated spectral type of B1 - B0.5. Shock-excited lines of H2 and [FeII], as seen in the near-infrared spectra obtained with UKIRT-UIST, lend support to this picture. Cold dust emission shows a massive clump of mass 1375 M{\sun} enveloping G12.42+0.50. Study of the molecular gas kinematics using the MALT90 and JCMT archival data unravels the presence of both infall activity and large-scale outflow suggesting an early stage of massive star formation in G12.42+0.50. A network of filamentary features are also revealed merging with the massive clump mimicking a hub-filament layout. Velocity structure along these indicate bulk inflow motion.

The ANtarctic Impulsive Transient Antenna (ANITA) NASA long-duration balloon payload completed its fourth flight in December 2016, after 28 days of flight time. ANITA is sensitive to impulsive broadband radio emission from interactions of ultra-high-energy neutrinos in polar ice (Askaryan emission). We present the results of two separate blind analyses searching for signals from Askaryan emission in the data from the fourth flight of ANITA. The more sensitive analysis, with a better expected limit, has a background estimate of $0.64^{+0.69}_{-0.45}$ and an analysis efficiency of $82\pm2\%$. The second analysis has a background estimate of $0.34^{+0.66}_{-0.16}$ and an analysis efficiency of $71\pm6\%$. Each analysis found one event in the signal region, consistent with the background estimate for each analysis. The resulting limit further tightens the constraints on the diffuse flux of ultra-high-energy neutrinos at energies above $10^{19.5}$ eV.

iPTF14hls is an unusually bright, long-lived II-P supernova (SN), whose light curve has at least five peaks. We propose that the outflows from the black hole hyperaccretion systems in the center of the collapsars should continuously inject into the envelope. For a jet-driven core-collapsar model, the outflow feedback results in prolonging the accretion timescale and fluctuating accretion rates in our analytic solutions. Thus, the long period of luminous, varying SN iPTF14hls might originate from the choked jets, which are regulated by the feedback of the strong disk outflows in a massive core-collapsar. One can expect that jet-driven iPTF14hls may last no more than approximately 3,000 days, and the luminosity may quickly decrease in the later stages. Moreover, the double-peak light curves in some SNe might be explained by the outflow feedback mechanism.

Recently, the ANITA collaboration reported on two upward-going extensive air shower events consistent with a primary particle that emerges from the surface of the ice. These events may be of $\nu_\tau$ origin, in which the neutrino interacts within the Earth to produce a $\tau$ lepton that emerges from the Earth, decays in the atmosphere, and initiates an extensive air shower. In this paper we estimate an upper bound on the ANITA acceptance to a diffuse $\nu_\tau$ flux detected via $\tau$-lepton-induced air showers within the bounds of Standard Model (SM) uncertainties. By comparing this estimate with the acceptance of Pierre Auger Observatory and IceCube and assuming SM interactions, we conclude that a $\nu_\tau$ origin of these events would imply a neutrino flux at least two orders of magnitude above current bounds.

In this paper we analyse the implications of the latest cosmological data sets to test the Etherington's distance duality relation (DDR), which connects the luminosity distance $D_L$ and angular diameter distance $D_A$ at the same redshift. For $D_L$ we consider the simulated data of gravitational waves from the third-generation gravitational wave detector (the Einstein Telescope, ET), which can be considered as standard candles (or standard siren), while the angular diameter distances $D_A$ are derived from the newly-compiled sample of compact radio quasars observed by very-long-baseline interferometry (VLBI), which represents a type of new cosmological standard ruler. Alleviating the absorption and scattering effects of dust in the Universe, this will create a valuable opportunity to directly test DDR at much higher precision with the combination of gravitational wave (GW) and electromagnetic (EM) signals. Our results show that, with the combination of the current radio quasar observations, the duality-distance relation can be verified at the precision of $10^{-2}$. Moreover, the Einstein Telescope ET would produce more robust constraints on the validity of such distance duality relation (at the precision of $10^{-3}$), with a larger sample of compact milliarcsecond radio quasars detected in future VLBI surveys.

In this paper, we assemble a well-defined sample of early-type gravitational lenses extracted from a large collection of 158 systems, and use the redshift distribution of galactic-scale lenses to test the standard cosmological model ($\Lambda$CDM) and the modified gravity theory (DGP). Two additional sub-samples are also included to account for possible selection effect introduced by the detectability of lens galaxies. Our results show that independent measurement of the matter density parameter ($\Omega_m$) could be expected from such strong lensing statistics. Based on future measurements of strong lensing systems from the forthcoming LSST survey, one can expect $\Omega_m$ to be estimated at the precision of $\Delta\Omega_m\sim 0.006$, which provides a better constraint on $\Omega_m$ than \textit{Planck} 2015 results. Moreover, use the lens redshift test is also used to constrain the characteristic velocity dispersion of the lensing galaxies, which is well consistent with that derived from the optical spectroscopic observations. A parameter $f_E$ is adopted to quantify the relation between the lensing-based velocity dispersion and the corresponding stellar value. Finally, the accumulation of detectable galactic lenses from future LSST survey would lead to more stringent fits of $\Delta f_E\sim10^{-3}$, which encourages us to test the global properties of early-type galaxies at much higher accuracy.

Gravitationally bound hierarchies containing three or more components are very common in our Universe. In this paper we study {\em periodic} gravitational wave (GW) form, their polarizations, response function, its Fourier transform, and energy loss rate of a triple system through three different channels of radiation, the scalar, vector and tensor modes, in Einstein-aether theory of gravity. In the weak-field approximations and with the recently obtained constraints of the theory, we first analyze the energy loss rate of a binary system, and find that the dipole contributions from the scalar and vector modes could be of the order of ${\cal{O}}\left(c_{14}\right){\cal{O}}\left(G_Nm/d\right)^2$, where $c_{14} \; (\equiv c_{1} + c_{4})$ is constrained to $c_{14} \lesssim {\cal{O}}\left(10^{-5}\right)$ by current observations, where $c_i$'s are the four coupling constants of the theory. On the other hand, the "strong-field" effects for a binary system of neutron stars are about six orders lower than that of GR. So, in this paper we ignore these "strong-field" effects and first develop the general formulas to the lowest post-Newtonian order, by taking the coupling of the aether field with matter into account. Within this approximation, we find that the scalar breather mode and the scalar longitudinal mode are all suppressed by a factor of ${\cal{O}}\left(c_{14}\right)$ with respect to the transverse-traceless modes ($h_{+}$ and $h_{\times}$), while the vectorial modes $(h_{X}$ and $h_{Y}$) are suppressed by a factor of $c_{13} \lesssim {\cal{O}}\left(10^{-15}\right)$. Applying the general formulas to a triple system with periodic orbits, we find that the corresponding GW form, response function, and its Fourier transform depend sensitively on both the configuration of the triple system and their orientations with respect to the detectors.

Using data collected by the Askaryan Radio Array (ARA) experiment at the South Pole, we have used long-baseline propagation of radio-frequency signals to extract information on the radio-frequency index-of-refraction in South Polar ice. Owing to the increasing ice density over the upper 150--200 meters, rays are observed along two, nearly parallel paths, one of which is direct and a second which refracts through an inflection point, with differences in both arrival time and arrival angle that can be used to constrain the neutrino properties. We also observe indications, for the first time, of radio-frequency ice birefringence for signals propagating along predominantly horizontal trajectories, corresponding to an asymmetry of order 0.1% between the ordinary and extra-ordinary paths, numerically compatible with previous measurements of birefringent asymmetries for vertically-propagating radio-frequency signals at South Pole. Taken together, these effects offer the possibility of redundantly measuring the range from receiver to a neutrino interaction in Antarctic ice, if receiver antennas are deployed at shallow (25 m<z<100 m) depths. Such range information is essential in determining both the neutrino energy, as well as the incident neutrino direction.

A survey of C2H N=1-0 and N2H+ J=1-0 toward Planck Galactic cold clumps (PGCCs) was performed using the Purple Mountain Observatory's 13.7 m telescope. C2H and N2H+ were chosen to study the chemical evolutionary states of PGCCs. Among 121 observed molecular cores associated with PGCCs, 71 and 58 are detected with C2H N=1-0 and N2H+ J=1-0, respectively. The detected lines of most sources can be fitted with a single component with compatible Vlsr and line widths, which confirms that these PGCC cores are very cold (with gas temperatures 9-21 K) and quiescent while still dominanted by turbulence. The ratio between the column densities of C2H and N2H+ (N(C2H)/N(N2H+)) is found to be a good tracer for the evolutionary states of PGCC cores. Gas-grain chemical model can reproduce the decreasing trend of N(C2H)/N(N2H+) as a function of time. The cores with the lowest abundances of N2H+ (X[N2H+] < 10^{-10}) are the youngest, and have nearly constant abundances of C2H. In evolved cores with X[N2H+] ~ 1E-9, abundances of C2H drop quickly as the exhaustion of carbon atoms. Although these PGCC cores are in different evolutionary states, they are all quite young (<5E5 yr) with N(C2H) > N(N2H+). Mapping observations are carried out toward 20 PGCC cores. The PGCC cores in Cepheus have lower N(C2H)/N(N2H+) and larger line widths compared with those in Taurus. This implies that PGCC cores in Taurus are less chemically evolved than those in Cepheus.

Long-duration gamma-ray bursts (LGRBs) associated with supernovae (SNe) are possibly born out of the death of a massive star. After the star collapses, a stellar-mass black hole (BH) is formed, surrounded by a hyperaccretion disk with outflows. Blandford-Znajek jets can be launched and then break out from the envelope to power LGRBs. The jet luminosity depends on the net inflow accretion rate at the inner radius of the disk. Furthermore, $\rm ^{56}Ni$ synthesis should occur in the strong outflows from the accretion disk. The decay of $\rm ^{56}Ni$ is considered to be the possible origin of SN bumps in the subsequent optical afterglows of LGRBs. If $\rm ^{56}Ni$ originates entirely from the outflows, there is competition between the luminosities of LGRBs and those of the corresponding $\rm ^{56}Ni$ bumps because of the material distribution between the disk inflows and outflows. In this paper, we investigated these two luminosities based on 15 cases of LGRB-SN in the framework of the BH hyperaccretion inflow-outflow model. Then, one can constrain the characteristics of the progenitor stars of these LGRBs. The results indicate that these LGRBs may originate from the low-metallicity ($Z\lesssim 10^{-2}Z_{\odot}$, where $Z$ and $Z_{\odot}$ are the metallicities of the stars and the Sun, respectively) stars or some massive solar-metallicity stars. For ultra-LGRBs (ULGRBs), such as GRB 111209A, most of the massive low-metallicity stars with $Z \lesssim 10^{-2}Z_{\odot}$ could be progenitors only if very strong outflows are launched from the disks. When the contributions of nucleosynthesis in the disk outflows are considered, there is no shortage of $\rm ^{56}Ni$ mass for luminous SNe associated with ULGRBs.

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.

Galaxy interactions are often accompanied by an enhanced star formation rate (SFR). Since molecular gas is essential for star formation, it is vital to establish whether, and by how much, galaxy interactions affect the molecular gas properties. We investigate the effect of interactions on global molecular gas properties by studying a sample of 58 galaxies in pairs and 154 control galaxies. Molecular gas properties are determined from observations with the JCMT, PMO, CSO telescopes, and supplemented with data from the xCOLD GASS and JINGLE surveys at $^{12}$CO(1-0) and $^{12}$CO(2-1). The SFR, gas mass ($M_\mathrm{H_{2}}$), and gas fraction ($f_{gas}$) are all enhanced in galaxies in pairs by $\sim$ 2.5 times compared to the controls matched in redshift, mass, and effective radius, while the enhancement of star formation efficiency (SFE $\equiv$ SFR/$M_{H_{2}}$) is less than a factor of 2. We also find that the enhancements in SFR, $M_{H_{2}}$ and $f_{gas}$ increase with decreasing pair separation and are larger in systems with smaller stellar mass ratio. Conversely, the SFE is only enhanced in close pairs (separation $<$ 20 kpc) and equal-mass systems; therefore most galaxies in pairs lie in the same parameter space on the SFR-$M_{H_{2}}$ plane as controls. This is the first time that the dependence of molecular gas properties on merger configurations is probed statistically with a relatively large sample and with a carefully-selected control sample for individual galaxies. We conclude that galaxy interactions do modify the molecular gas properties, although the strength of the effect is merger configuration dependent.

Ultra-high energy neutrinos are detectable through impulsive radio signals generated through interactions in dense media, such as ice. Subsurface in-ice radio arrays are a promising way to advance the observation and measurement of astrophysical high-energy neutrinos with energies above those discovered by the IceCube detector ($\geq$1 PeV) as well as cosmogenic neutrinos created in the GZK process ($\geq$100 PeV). Here we describe the $\textit{NuPhase}$ detector, which is a compact receiving array of low-gain antennas deployed 185 m deep in glacial ice near the South Pole. Signals from the antennas are digitized and coherently summed into multiple beams to form a low-threshold interferometric phased array trigger for radio impulses. The NuPhase detector was installed at an Askaryan Radio Array (ARA) station during the 2017/18 Austral summer season. $\textit{In situ}$ measurements with an impulsive, point-source calibration instrument show a 50% trigger efficiency on impulses with voltage signal-to-noise ratios (SNR) of $\le$2.0, a factor of $\sim$1.8 improvement in SNR over the standard ARA combinatoric trigger. Hardware-level simulations, validated with $\textit{in situ}$ measurements, predict a trigger threshold of an SNR as low as 1.6 for neutrino interactions that are in the far field of the array. With the already-achieved NuPhase trigger performance included in ARASim, a detector simulation for the ARA experiment, we find the trigger-level effective detector volume is increased by a factor of 1.8 at neutrino energies between 10 and 100 PeV compared to the currently used ARA combinatoric trigger. We also discuss an achievable near term path toward lowering the trigger threshold further to an SNR of 1.0, which would increase the effective single-station volume by more than a factor of 3 in the same range of neutrino energies.

We investigate the formation and magnetic topology of four flare/CME events with filament-sigmoid systems, in which the sigmoidal hot channels are located above the filaments, and they appear in pairs prior to eruption.} The formation of hot channels usually takes several to dozens of hours during which two J-shape sheared arcades gradually evolve into sigmoidal hot channels, then they keep stable for tens of minutes or hours and erupt. While the low-lying filaments show no significant change. \textbf{We construct a series of magnetic field models and find that the best-fit preflare models contain magnetic flux ropes with hyperbolic flux tubes (HFTs). The field lines above the HFT correspond to the high-lying hot channel, while those below the HFT surround the underlying filaments. In particular, the continuous and long field lines representing the flux rope located above the HFT match the observed hot channels well in three events. While for SOL2014-04-18 event, the flux bundle that mimics the observed hot channel is located above the flux rope. The flux rope axis lies in a height range of 19.8 Mm and 46 Mm above the photosphere for the four events, among which the flux rope axis in SOL2012-07-12 event has a maximum height, which probably explains why it is often considered as a double-decker structure. Our modeling suggests that the high-lying hot channel may be formed by magnetic reconnections between sheared field lines occurring above the filament prior to eruption.

Previous observations have revealed an accretion disk and outflow motion in high-mass star-forming region G192.16-3.84. While collapse have not been reported before. We present here molecular line and continuum observations toward massive core G192.16-3.84 with the Submillimeter Array. C$^{18}$O(2-1) and HCO$^{+}$(3-2) lines show pronounced blue profiles, indicating gas infalling in this region. This is the first time that the infall motion has been reported in G192.16-3.84 core. Two-layer model fitting gave infall velocities of 2.0$\pm$0.2 and 2.8$\pm$0.1 km s$^{-1}$. Assuming that the cloud core follows a power-law density profile ($\rho$$\propto$$r^{1.5}$), the corresponding mass infall rates are (4.7$\pm$1.7)$\times10^{-3}$ and (6.6$\pm$2.1)$\times10^{-3}$ M$_{\sun}$ yr$^{-1}$ for C$^{18}$O(2-1) and HCO$^{+}$(3-2), respectively. The derived infall rates are in agreement with the turbulent core model and those in other high-mass star-forming regions, suggesting that high accretion rate is a general requirement to form a massive star.

These proceedings address a recent publication by the ANITA collaboration of four upward- pointing cosmic-ray-like events observed in the first flight of ANITA. Three of these events were consistent with stratospheric cosmic-ray air showers where the axis of propagation does not inter- sect the surface of the Earth. The fourth event was consistent with a primary particle that emerges from the surface of the ice suggesting a possible {\tau}-lepton decay as the origin of this event. These proceedings follow-up on the modeling and testing of the hypothesis that this event was of {\tau} neutrino origin.

The dust sub-millimetre polarisation of star-forming clouds carries information on dust and the role of magnetic fields in cloud evolution. With observations of a dense filamentary cloud G035.39-00.33, we aim to characterise the dust emission properties and the variations of the polarisation fraction. JCMT SCUBA-2/POL-2 data at 850um are combined with Planck 850um (353GHz) data to map polarisation fractions. With previous SCUBA-2 observations (450um and 850um) and Herschel data, the column densities are determined via modified blackbody fits and via radiative transfer modelling. Models are constructed to examine how the polarisation angles and fractions depend on potential magnetic field geometries and grain alignment. POL-2 data show clear changes in the magnetic field orientation. The filament has a peak column density of N(H2)~7 10^22 cm-2, a minimum dust temperature of T~12 K, and a mass of some 4300Msun for the area N(H2)> 5 10^21 cm-2. The estimated average value of the dust opacity spectral index is beta ~ 1.9. The ratio of sub-millimetre and J band optical depths is tau(250 um)/tau(J) ~ 2.5 10^-3, more than four times the typical values for diffuse medium. The polarisation fraction decreases as a function of column density to p ~ 1% in the central filament. Because of noise, the observed decrease of p(N) is significant only at N(H2)>2 10^22 cm-2. The observations suggest that the grain alignment is not constant. Although the data can be explained with a complete loss of alignment at densities above ~ 10^4 cm-3 or using the predictions of radiative torques alignment, the uncertainty of the field geometry and the spatial filtering of the SCUBA-2 data prevent strong conclusions. G035.39-00.33 shows strong signs of dust evolution and the low polarisation fraction is suggestive of a loss of polarised emission from its densest parts.

We present a multiwavelength study of an atypical submillimeter galaxy in the GOODS-North field, with the aim to understand its physical properties of stellar and dust emission, as well as the central AGN activity. Although it is shown that the source is likely an extremely dusty galaxy at high redshift, its exact position of submillimeter emission is unknown. With the new NOEMA interferometric imaging, we confirm that the source is a unique dusty galaxy. It has no obvious counterpart in the optical and even NIR images observed with HST at lambda~<1.4um. Photometric-redshift analyses from both stellar and dust SED suggest it to likely be at z~>4, though a lower redshift at z~>3.1 cannot be fully ruled out (at 90% confidence interval). Explaining its unusual optical-to-NIR properties requires an old stellar population (~0.67 Gyr), coexisting with a very dusty ongoing starburst component. The latter is contributing to the FIR emission, with its rest-frame UV and optical light being largely obscured along our line of sight. If the observed fluxes at the rest-frame optical/NIR wavelengths were mainly contributed by old stars, a total stellar mass of ~3.5x10^11Msun would be obtained. An X-ray spectral analysis suggests that this galaxy harbors a heavily obscured AGN with N_H=3.3x10^23 cm^-2 and an intrinsic 2-10 keV luminosity of L_X~2.6x10^44 erg/s, which places this object among distant type 2 quasars. The radio emission of the source is extremely bright, which is an order of magnitude higher than the star-formation-powered emission, making it one of the most distant radio-luminous dusty galaxies. The combined characteristics of the galaxy suggest that the source appears to have been caught in a rare but critical transition stage in the evolution of submillimeter galaxies, where we are witnessing the birth of a young AGN and possibly the earliest stage of its jet formation and feedback.

We investigate the sensitivity of reheating history on the abundance of primordial black holes (PBHs). Contrary to the monochromatic case of mass fraction of PBH, reheating era with different e-folding and equation-of-state can have a substantial impact on the abundance of PBH with an extended mass fraction. We demonstrate explicitly this reheating sensitivity in an illustrative model of single field inflation with a quasi-inflection point, and find that both the peak position and amplitude of the extended mass fraction as well as the abundance of PBH in DM can vary by many orders of magnitude, which adds another layer of uncertainty on the PBH scenarios as dark matter.

The radiation mechanism of the prompt emission of gamma-ray bursts (GRBs) remains an open question. Although their spectra are usually well fitted with the empirical Band function, which is widely believed to be fully non-thermal and interpreted as an optically thin synchrotron emission, accumulating evidence shows that a thermal component actually exists. In this paper, a multi-color blackbody (mBB) model is proposed for the time-integrated spectrum of GRB 081221 by assuming a power-law distribution of the thermal luminosities with temperature, which manifests photospheric emissions from a different radius and/or angle. The effects of the minimum temperature $kT_{\min }$, the maximum temperature $kT_{\max }$ and the power law index $m$ of the luminosity distribution of an mBB are discussed. The fitting to the time-integrated spectrum during the bright phase (from 20s to 30s since the trigger) of GRB 081221 by the mBB model yields $kT_{\min }$ = $4.4\pm 0.3$ keV, $kT_{\max }$ = 57.0 $_{-1.4}^{+1.6}$ keV, and $m=-0.46_{-0.06}^{+0.05}$. When the time bin is small enough, the time-resolved spectra of GRB 081221 are well fitted with a series of single-temperature blackbodies. Our results imply the prompt emission of GRB 081221 is dominated by the photosphere emission and its time-integrated spectrum is a superposition of pure blackbody components at different times, indicating that some empirical Band spectra may be interpreted as mBB if the temperature is widely distributed.

The Askaryan Radio Array (ARA) reports an observation of radio emission coincident with the "Valentine's Day" solar flare on Feb. 15$^{\rm{th}}$, 2011 in the prototype "Testbed" station. We find $\sim2000$ events that passed our neutrino search criteria during the 70 minute period of the flare, all of which reconstruct to the location of the sun. A signal analysis of the events reveals them to be consistent with that of bright thermal noise correlated across antennas. This is the first natural source of radio emission reported by ARA that is tightly reconstructable on an event-by-event basis. The observation is also the first for ARA to point radio from individual events to an extraterrestrial source on the sky. We comment on how the solar flares, coupled with improved systematic uncertainties in reconstruction algorithms, could aid in a mapping of any above-ice radio emission, such as that from cosmic-ray air showers, to astronomical locations on the sky.

Using the X-ray-selected active galactic nuclei (AGN) from the XMM-XXL north survey and the SDSS Baryon Oscillation Spectroscopic Survey (BOSS) spectroscopic follow-up of them, we compare the properties of X-ray unobscured and obscured broad-line AGN (BLAGN1 and BLAGN2; $N_\textrm{H}$below and above $10^{21.5}$ cm$^{-2}$), including their X-ray luminosity $L_X$, black hole mass, Eddington ratio $\lambda_{\textrm{Edd}}$, optical continuum and line features. We find that BLAGN2 have systematically larger broad line widths and hence apparently higher (lower) $M_{\textrm{BH}}$ ($\lambda_{\textrm{Edd}}$) than BLAGN1. We also find that the X-ray obscuration in BLAGN tends to coincide with optical dust extinction, which is optically thinner than that in narrow-line AGN (NLAGN) and likely partial-covering to the broad line region. All the results can be explained in the framework of a multi-component, clumpy torus model by interpreting BLAGN2 as an intermediate type between BLAGN1 and NLAGN in terms of an intermediate inclination angle.

The Tianlai Pathfinder is designed to demonstrate the feasibility of using a wide field of view radio interferometers to map the density of neutral hydrogen in the Universe after the Epoch of Reionizaton. This approach, called 21~cm intensity-mapping, promises an inexpensive means for surveying the large-scale structure of the cosmos. The Tianlai Pathfinder presently consists of an array of three, 15~m $\times$ 40~m cylinder telescopes and an array of sixteen, 6~m diameter dish antennas located in a radio-quiet part of western China. The two types of arrays were chosen to determine the advantages and disadvantages of each approach. The primary goal of the Pathfinder is to make 3D maps by surveying neutral hydrogen over large areas of the sky %$20,000 {\rm deg}^2$ in two different redshift ranges: first at $1.03 > z > 0.78$ ($700 - 800$~MHz) and later at $0.21 > z > 0.12$ ($1170 - 1270$~MHz). The most significant challenge to $21$~cm intensity-mapping is the removal of strong foreground radiation that dwarfs the cosmological signal. It requires exquisite knowledge of the instrumental response, i.e. calibration. In this paper, we provide an overview of the status of the Pathfinder and discuss the details of some of the analysis that we have carried out to measure the beam function of both arrays. We compare electromagnetic simulations of the arrays to measurements, discuss measurements of the gain and phase stability of the instrument, and provide a brief overview of the data processing pipeline.

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.

The Antarctic Impulsive Transient Antenna (ANITA), a NASA long-duration balloon payload, searches for radio emission from interactions of ultra-high-energy neutrinos in polar ice. The third flight of ANITA (ANITA-III) was launched in December 2014 and completed a 22-day flight. We present the results of three analyses searching for Askaryan radio emission of neutrino origin. In the most sensitive of the analyses, we find one event in the signal region on an expected a priori background of $0.7^{+0.5}_{-0.3}$. Though consistent with the background estimate, the candidate event remains compatible with a neutrino hypothesis even after additional post-unblinding scrutiny.

We report evidence of environmental quenching among galaxies at redshift ~ 2, namely the probability that a galaxy quenches its star formation activity is enhanced in the regions of space in proximity of other quenched, more massive galaxies. The effect is observed as strong clustering of quiescent galaxies around quiescent galaxies on angular scales \theta < 20 arcsec, corresponding to a proper(comoving) scale of 168 (502) kpc at z = 2. The effect is observed only for quiescent galaxies around other quiescent galaxies; the probability to find star-forming galaxies around quiescent or around star-forming ones is consistent with the clustering strength of galaxies of the same mass and at the same redshift, as observed in dedicated studies of galaxy clustering. The effect is mass dependent in the sense that the quenching probability is stronger for galaxies of smaller mass ($\rm{M_*<10^{10} Msun}$) than for more massive ones, i.e. it follows the opposite trend with mass relative to gravitational galaxy clustering. The spatial scale where the effect is observed suggests these environments are massive halos, in which case the observed effect would likely be satellite quenching. The effect is also redshift dependent in that the clustering strength of quiescent galaxies around other quiescent galaxies at z = 1.6 is ~ 1.7 times larger than that of the galaxies with the same stellar mass at z = 2.6. This redshift dependence allows for a crude estimate of the time scale of environmental quenching of low-mass galaxies, which is in the range 1.5 - 4 Gyr, in broad agreement with other estimates and with our ideas on satellite quenching.

We present the results of dust emission polarization measurements of Ophiuchus-B (Oph-B) carried out using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) camera with its associated polarimeter (POL-2) on the James Clerk Maxwell Telescope (JCMT) in Hawaii. This work is part of the B-fields In Star-forming Region Observations (BISTRO) survey initiated to understand the role of magnetic fields in star formation for nearby star-forming molecular clouds. We present a first look at the geometry and strength of magnetic fields in Oph-B. The field geometry is traced over $\sim$0.2 pc, with clear detection of both of the sub-clumps of Oph-B. The field pattern appears significantly disordered in sub-clump Oph-B1. The field geometry in Oph-B2 is more ordered, with a tendency to be along the major axis of the clump, parallel to the filamentary structure within which it lies. The degree of polarization decreases systematically towards the dense core material in the two sub-clumps. The field lines in the lower density material along the periphery are smoothly joined to the large scale magnetic fields probed by NIR polarization observations. We estimated a magnetic field strength of 630$\pm$410 $\mu$G in the Oph-B2 sub-clump using a Davis-Chandeasekhar-Fermi analysis. With this magnetic field strength, we find a mass-to-flux ratio $\lambda$= 1.6$\pm$1.1, which suggests that the Oph-B2 clump is slightly magnetically supercritical.

The occurrence rates of bright X-ray flares in z<1 gamma-ray bursts (GRBs) with or without observed supernovae (SNe) association were compared. Our Sample I: the z<1 long GRBs (LGRBs) with SNe association (SN-GRBs) and with early Swift/X-Ray Telescope (XRT) observations, consists of 18 GRBs, among which only two GRBs have bright X-ray flares. Our Sample II: for comparison, all the z<1 LGRBs without observed SNe association and with early Swift/XRT observations, consists of 45 GRBs, among which 16 GRBs present bright X-ray flares. Thus, the study indicates a lower occurrence rate of bright X-ray flares in Sample I (11.1%) than in Sample II (35.6%). In addition, if dim X-ray fluctuations are included as flares, then 16.7% of Sample I and 55.6% of Sample II are found to have flares, again showing the discrepancy between these two samples. We examined the physical origin of these bright X-ray flares and found that most of them are probably related to the central engine reactivity. To understand the discrepancy, we propose that such a lower occurrence rate of flares in the SN-GRB sample may hint at an energy partition among the GRB, SNe, and X-ray flares under a saturated energy budget of massive star explosion.

Based on the 850 $\mu$m dust continuum data from SCUBA-2 at James Clerk Maxwell Telescope (JCMT), we compare overall properties of Planck Galactic Cold Clumps (PGCCs) in the $\lambda$ Orionis cloud to those of PGCCs in the Orion A and B clouds. The Orion A and B clouds are well known active star-forming regions, while the $\lambda$ Orionis cloud has a different environment as a consequence of the interaction with a prominent OB association and a giant Hii region. PGCCs in the $\lambda$ Orionis cloud have higher dust temperatures ($Td=16.13\pm0.15$ K) and lower values of dust emissivity spectral index ($ \beta=1.65\pm0.02$) than PGCCs in the Orion A (Td=13.79$\pm 0.21$K, $\beta=2.07\pm0.03$) and Orion B ($Td=13.82\pm0.19$K, $\beta=1.96\pm0.02$) clouds. We find 119 sub-structures within the 40 detected PGCCs and identify them as cores. Of total 119 cores, 15 cores are discovered in the $\lambda$ Orionis cloud, while 74 and 30 cores are found in the Orion A and B clouds, respectively. The cores in the $\lambda$ Orionis cloud show much lower mean values of size R=0.08 pc, column density N(H2)=$(9.5\pm1.2) \times 10^{22}$ cm$^{-2}$, number density n(H2)=$(2.9 \pm 0.4)\times10^{5}$ cm$^{-3}$, and mass $M_{core}$=$1.0\pm0.3$ M$_{\odot}$ compared to the cores in the Orion A (R=0.11pc, $N(H2)=(2.3\pm0.3) \times 10^{23}$ cm$^{-2}$, n(H2)=$(3.8\pm0.5) \times 10^{5}$cm$^{-3}$, and $M_{core}$=$2.4 \pm 0.3$ M$_{\odot}$) and Orion B (R=0.16pc, N(H2)=$(3.8 \pm 0.4) \times 10^{23}$cm$^{-2}$, n(H2)=$(15.6\pm1.8)\times10^{5}$ cm$^{-3}$, and $M_{core}$= $2.7\pm0.3$ M$_{\odot}$) clouds. These core properties in the $\lambda$ Orionis cloud can be attributed to the photodissociation and external heating by the nearby Hii region, which may prevent the PGCCs from forming gravitationally bound structures and eventually disperse them. These results support the idea of negative stellar feedback on core formation.

Graham et al. (2015a) reported a periodically varying quasar and supermassive black hole binary candidate, PG1302-102 (hereafter PG1302), which was discovered in the Catalina Real-Time Transient Survey (CRTS). Its combined Lincoln Near-Earth Asteroid Research (LINEAR) and CRTS optical light curve is well fitted to a sinusoid of an observed period of $\approx 1,884$ days and well modeled by the relativistic Doppler boosting of the secondary mini-disk (D'Orazio et al. 2015). However, the LINEAR+CRTS light curve from MJD $\approx 52700$ to MJD $\approx 56400$ covers only $\sim 2$ cycles of periodic variation, which is a short baseline that can be highly susceptible to normal, stochastic quasar variability (Vaughan et al. 2016). In this Letter, we present a re-analysis of PG1302, using the latest light curve from the All-Sky Automated Survey for Supernovae (ASAS-SN), which extends the observational baseline to the present day (MJD $\approx 58200$), and adopting a maximum likelihood method which searches for a periodic component in addition to stochastic quasar variability. When the ASAS-SN data are combined with the previous LINEAR+CRTS data, the evidence for periodicity decreases. For genuine periodicity one would expect that additional data would strengthen the evidence, so the decrease in significance may be an indication that the binary model is disfavored.

In order to understand the initial conditions and early evolution of star formation in a wide range of Galactic environments, we carried out an investigation of 64 \textit{Planck} Galactic Cold Clumps (PGCCs) in the second quadrant of the Milky Way. Using the $^{13}$CO and C$^{18}$O $J = 1 - 0$ lines, and 850\,$\mu$m continuum observations, we investigated cloud fragmentation and evolution associated with star formation. We extracted 468 clumps and 117 cores from the $^{13}$CO line and 850\,$\mu$m continuum maps, respectively. We make use of the Bayesian Distance Calculator and derived the distances of all 64 PGCCs. We found that in general, the mass-size plane follows a relation of $m\sim r^{1.67}$. At a given scale, the masses of our objects are around 1/10 of that of typical Galactic massive star-forming regions. Analysis of the clump and core masses, virial parameters, densities, and mass-size relation suggests that the PGCCs in our sample have a low core formation efficiency ($\sim$3.0\%), and most PGCCs are likely low-mass star-forming candidates. Statistical study indicates that the 850\,$\mu$m cores are more turbulent, more optically thick, and denser than the $^{13}$CO clumps for star formation candidates, suggesting that the 850\,$\mu$m cores are likely more appropriate future star-formation candidates than the $^{13}$CO clumps.

Using Galactic Plane surveys, we have selected a massive (1200 M$_\odot$), cold (14 K) 3.6-70 $\mu$m dark IRDC G331.372-00.116. This IRDC has the potential to form high-mass stars and, given the absence of current star formation signatures, it seems to represent the earliest stages of high-mass star formation. We have mapped the whole IRDC with the Atacama Large Millimeter/submillimeter Array (ALMA) at 1.1 and 1.3 mm in dust continuum and line emission. The dust continuum reveals 22 cores distributed across the IRDC. In this work, we analyze the physical properties of the most massive core, ALMA1, which has no molecular outflows detected in the CO (2-1), SiO (5-4), and H$_2$CO (3-2) lines. This core is relatively massive ($M$ = 17.6 M$_\odot$), subvirialized (virial parameter $\alpha_{vir}=M_{vir}/M=0.14$), and is barely affected by turbulence (transonic Mach number of 1.2). Using the HCO$^+$ (3-2) line, we find the first detection of infall signatures in a relatively massive, prestellar core (ALMA1) with the potential to form a high-mass star. We estimate an infall speed of 1.54 km s$^{-1}$ and a high accretion rate of 1.96 $\times$ 10$^{-3}$ M$_\odot$ yr$^{-1}$. ALMA1 is rapidly collapsing, out of virial equilibrium, more consistent with competitive accretion scenarios rather than the turbulent core accretion model. On the other hand, ALMA1 has a mass $\sim$6 times larger than the clumps Jeans mass, being in an intermediate mass regime ($M_{J}=2.7<M\lesssim$ 30 M$_\odot$), contrary to what both the competitive accretion and turbulent core accretion theories predict.

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

We present 850 $\mu$m imaging polarimetry data of the $\rho$ Oph-A core taken with the Submillimeter Common-User Bolometer Array-2 (SCUBA-2) and its polarimeter (POL-2), as part of our ongoing survey project, BISTRO (B-fields In STar forming RegiOns). The polarization vectors are used to identify the orientation of the magnetic field projected on the plane of the sky at a resolution of 0.01 pc. We identify 10 subregions with distinct polarization fractions and angles in the 0.2 pc $\rho$ Oph A core; some of them can be part of a coherent magnetic field structure in the $\rho$ Oph region. The results are consistent with previous observations of the brightest regions of $\rho$ Oph-A, where the degrees of polarization are at a level of a few percents, but our data reveal for the first time the magnetic field structures in the fainter regions surrounding the core where the degree of polarization is much higher ($> 5 \%$). A comparison with previous near-infrared polarimetric data shows that there are several magnetic field components which are consistent at near-infrared and submillimeter wavelengths. Using the Davis-Chandrasekhar-Fermi method, we also derive magnetic field strengths in several sub-core regions, which range from approximately 0.2 to 5 mG. We also find a correlation between the magnetic field orientations projected on the sky with the core centroid velocity components.

The detections of some long gamma-ray bursts (LGRBs) relevant to mergers of neutron star (NS)-NS or black hole (BH)-NS, as well as some short gamma-ray bursts (SGRBs) probably produced by collapsars, muddle the boundary of two categories of gamma-ray bursts (GRBs). In both cases, a plausible candidate of central engine is a BH surrounded by a hyperaccretion disc with strong outflows, launching relativistic jets driven by Blandford-Znajek mechanism. In the framework of compact binary mergers, we test the applicability of the BH hyperaccretion inflow-outflow model on powering observed GRBs. We find that, for a low outflow ratio, $\sim 50\%$, postmerger hyperaccretion processes could power not only all SGRBs but also most of LGRBs. Some LGRBs might do originate from merger events in the BH hyperaccretion scenario, at least on the energy requirement. Moreover, kilonovae might be produced by neutron-rich outflows, and their luminosities and timescales significantly depend on the outflow strengths. GRBs and their associated kilonovae are competitive with each other on the disc mass and total energy budgets. The stronger the outflow, the more similar the characteristics of kilonovae to supernovae (SNe). This kind of `nova' might be called `quasi-SN'.

Short gamma-ray bursts may originate from the merger of double neutron stars (NS) or that of a black hole (BH) and an NS. We propose that the bright X-ray flare related to the central engine reactivity may hint a BH-NS merger, since such a merger can provide more fall-back materials and therefore a more massive accretion disk than the NS-NS merger. Based on the observed 49 short bursts with Swift/X-ray Telescope follow-up observations, we find that three bursts have bright X-ray flares, among which three flares from two bursts are probably related to the central engine reactivity. We argue that these two bursts may originate from the BH-NS merger rather than the NS-NS merger. Our suggested link between the central engine-powered bright X-ray flare and the BH-NS merger event can be checked by the future gravitational wave detections from advanced LIGO and Virgo.

We report on an upward traveling, radio-detected cosmic-ray-like impulsive event with characteristics closely matching an extensive air shower. This event, observed in the third flight of the Antarctic Impulsive Transient Antenna (ANITA), a NASA-sponsored long-duration balloon payload, is consistent with a similar event reported in a previous flight. These events may be produced by the atmospheric decay of an upward-propagating $\tau$-lepton produced by a $\nu_{\tau}$ interaction, although their relatively steep arrival angles create tension with the standard model (SM) neutrino cross section. Each of the two events have $a~posteriori$ background estimates of $\lesssim 10^{-2}$ events. If these are generated by $\tau$-lepton decay, then either the charged-current $\nu_{\tau}$ cross section is suppressed at EeV energies, or the events arise at moments when the peak flux of a transient neutrino source was much larger than the typical expected cosmogenic background neutrinos.

We study the first electromagnetic pulse after the gravitational wave chirp signal, focusing on the profile and duration. It is found that the light curve, especially the steep decay (SD) phase, can be very different by adopting different viewing angle $\theta_{\rm view}$ on the jet shell. For an on-axis jet with a power-law radiation spectrum, the observed flux in the SD is proportional to $t_{\rm{obs}}^{-2-\beta}$ with $\beta$ being the spectral index and $t_{\rm{obs}}$ being the observer time. Here, $t_{\rm{obs}}=0$ is set at the observed time of the jet ejected from the central engine. The SD may become steep by increasing $\theta_{\rm view}$. We also study the bolometric luminosity $L$ from a jet shell with a non-power-law radiation spectrum. For an on-axis jet, $L{\propto}t_{\rm{obs}}^{-3}$ is found in the SD. However, the SD is steeper than $L{\propto}t_{\rm{obs}}^{-3}$ for the radiation from an off-axis jet. The higher value of $\theta_{\rm view}$ is, the steeper of SD would be. Then, we suggest that the SD phase can be used to discriminate an off-axis jet from an on-axis jet. The reason for above behaviors is discussed. In addition, we find that the duration of first electromagnetic pulse is close to its peak time, especially for $\theta_{\rm{view}}\sim20^\circ$. This result is consistent with that found in GW~170817/GRB~170817A. Thus, the jet corresponding to the prompt emission of GRB~170817A should be immediately ejected after the merger. Our results also reveal that the duration of the first electromagnetic pulse can provide the information of the time to search gravitational waves.

Planck Galactic Cold Clumps (PGCCs) possibly represent the early stages of star formation. To understand better the properties of PGCCs, we studied 16 PGCCs in the L1495 cloud with molecular lines and continuum data from Herschel, JCMT/SCUBA-2 and the PMO 13.7 m telescope. Thirty dense cores were identified in 16 PGCCs from 2-D Gaussian fitting. The dense cores have dust temperatures of $T_{\rm d}$ = 11-14 K, and H$_{2}$ column densities of $N_{\rm H_{2}}$ = 0.36-2.5$\times10^{22}$ cm$^{-2}$. We found that not all PGCCs contain prestellar objects. In general, the dense cores in PGCCs are usually at their earliest evolutionary stages. All the dense cores have non-thermal velocity dispersions larger than the thermal velocity dispersions from molecular line data, suggesting that the dense cores may be turbulence-dominated. We have calculated the virial parameter $\alpha$ and found that 14 of the dense cores have $\alpha$ $<$ 2, while 16 of the dense cores have $\alpha$ $>$ 2. This suggests that some of the dense cores are not bound in the absence of external pressure and magnetic fields. The column density profiles of dense cores were fitted. The sizes of the flat regions and core radii decrease with the evolution of dense cores. CO depletion was found to occur in all the dense cores, but is more significant in prestellar core candidates than in protostellar or starless cores. The protostellar cores inside the PGCCs are still at a very early evolutionary stage, sharing similar physical and chemical properties with the prestellar core candidates.

In recent years, superfluid dark matter (SfDM) has become a competitive model of emergent modified Newtonian dynamics (MOND) scenario: MOND phenomenons naturally emerge as a derived concept due to an extra force mediated between baryons by phonons as a result of axionlike particles condensed as superfluid at galactic scales; Beyond galactic scales, these axionlike particles behave as normal fluid without phonon-mediated MOND-like force between baryons, therefore SfDM also maintains the usual success of $\Lambda$CDM at cosmological scales. In this paper, we use gravitational waves (GWs) to probe the relevant parameter space of SfDM. GWs through Bose-Einstein condensate (BEC) could propagate with a speed slightly deviation from the speed-of-light due to the change in the effective refractive index, which depends on the SfDM parameters and GW-source properties. We find that Five hundred meter Aperture Spherical Telescope (FAST), Square Kilometre Array (SKA) and International Pulsar Timing Array (IPTA) are the most promising means as GW probe of relevant parameter space of SfDM. Future space-based GW detectors are also capable of probing SfDM if a multimessenger approach is adopted.

Outflows launched from a remnant disc of compact binary merger may have essential contribution to the kilonova emission. Numerical calculations are conducted in this work to study the structure of accretion flows and outflows. By the incorporation of limited-energy advection in the hyper-accretion discs, outflows occur naturally from accretion flows due to imbalance between the viscous heating and the sum of the advective and radiative cooling. Following this spirit, we revisit the properties of the merger outflow ejecta. Our results show that around $10^{-3} \sim 10^{-1} M_\odot$ of the disc mass can be launched as powerful outflows. The amount of unbound mass varies with the disc mass and the viscosity. The outflow-contributed peak luminosity is around $10^{40} \sim 10^{41}$ erg/s. Such a scenario can account for the observed kilonovae associated with short gamma-ray bursts, including the recent event AT2017gfo (GW170817).

Two plateaus and one following bump in the X-ray light curve of GRB 170714A have been detected by the \textit{Swift}/X-Ray Telescope, which could be very meaningful for the central engine of gamma-ray bursts (GRBs), implying that the origin of this burst might be different from that of other ultra-long GRBs. We propose that merging two neutron stars into a hyper-massive quark star (QS) and then collapsing into a black hole (BH), with a delay time around $10^4$~s, could be responsible for those X-ray components. The hyper-massive QS is initially in a fluid state, being turbulent and differentially rotating, but would be solidified and release its latent heat injected into the GRB fireball (lasting about $10^3$~s during the liquid-solid phase transition). Magnetic field as high as $\sim 10^{15}$~G could be created by dynamo action of the newborn liquid QS, and a magnetar-like central engine (after solidification) supplies significant energy for the second plateau. More energy could be released during a fall-back accretion after the post-merger QS collapses to a BH, and the X-ray bump forms. This post-merger QS model might be tested by future observations, with either advanced gravitational wave detectors (e.g., advanced LIGO and VIRGO) or X-ray/optical telescopes.

The balloon-borne HiCal radio-frequency (RF) transmitter, in concert with the ANITA radio-frequency receiver array, is designed to measure the Antarctic surface reflectivity in the RF wavelength regime. The amplitude of surface-reflected transmissions from HiCal, registered as triggered events by ANITA, can be compared with the direct transmissions preceding them by O(10) microseconds, to infer the surface power reflection coefficient $\cal{R}$. The first HiCal mission (HiCal-1, Jan. 2015) yielded a sample of 100 such pairs, resulting in estimates of $\cal{R}$ at highly-glancing angles (i.e., zenith angles approaching $90^\circ$), with measured reflectivity for those events which exceeded extant calculations. The HiCal-2 experiment, flying from Dec., 2016-Jan., 2017, provided an improvement by nearly two orders of magnitude in our event statistics, allowing a considerably more precise mapping of the reflectivity over a wider range of incidence angles. We find general agreement between the HiCal-2 reflectivity results and those obtained with the earlier HiCal-1 mission, as well as estimates from Solar reflections in the radio-frequency regime. In parallel, our calculations of expected reflectivity have matured; herein, we use a plane-wave expansion to estimate the reflectivity R from both a flat, smooth surface (and, in so doing, recover the Fresnel reflectivity equations) and also a curved surface. Multiplying our flat-smooth reflectivity by improved Earth curvature and surface roughness corrections now provides significantly better agreement between theory and the HiCal 2a/2b measurements.

Gamma-ray bursts (GRBs) are luminous and violent phenomena in the universe. Traditionally, long GRBs are expected to be produced by the collapse of massive stars and associated with supernovae. However, some low-redshift long GRBs have no detection of supernova association, such as GRBs 060505, 060614 and 111005A. It is hard to classify these events convincingly according to usual classifications, and the lack of the supernova implies a non-massive star origin. We propose a new path to produce long GRBs without supernova association, the unstable and extremely violent accretion in a contact binary system consisting of a stellar-mass black hole and a white dwarf, which fills an important gap in compact binary evolution.

We make the first data release (DR1) of the ongoing ESO Public Survey SAMPLING (http://dx.doi.org/10.7910/DVN/0L8NHX). DR1 comprises of 124 fields distributed in $70^\circ < l < 216^\circ$, $-35^\circ < b < 25^\circ$. The 12CO and 13CO (2-1) cubes are gridded in $8"$ pixels, with an effective resolution of $36"$. The channel width is 0.33 km/s and the RMS noise is $T_{\rm mb}<0.2$ K. Once completed, SAMPLING and complementary surveys will initiate the first major step forward to characterize molecular clouds and star formation on truly Galactic scales.

The $f(R)$ gravity can be cast into the form of a scalar-tensor theory, and scalar degree of freedom can be suppressed in high-density regions by the chameleon mechanism. In this article, for the general $f(R)$ gravity, using a scalar-tensor representation with the chameleon mechanism, we calculate the parameterized post-Newtonian parameters $\gamma$ and $\beta$, the effective gravitational constant $G_{\rm eff}$, and the effective cosmological constant $\Lambda_{\rm eff}$. In addition, for the general $f(R)$ gravity, we also calculate the rate of orbital period decay of the binary system due to gravitational radiation. Then we apply these results to specific $f(R)$ models (Hu-Sawicki model, Tsujikawa model and Starobinsky model) and derive the constraints on the model parameters by combining the observations in solar system, cosmological scales and the binary systems.

Analysis of all-sky Planck submillimetre observations and the IRAS 100um data has led to the detection of a population of Galactic cold clumps. The clumps can be used to study star formation and dust properties in a wide range of Galactic environments. Our aim is to measure dust spectral energy distribution (SED) variations as a function of the spatial scale and the wavelength. We examine the SEDs at large scales using IRAS, Planck, and Herschel data. At smaller scales, we compare with JCMT/SCUBA-2 850um maps with Herschel data that are filtered using the SCUBA-2 pipeline. Clumps are extracted using the Fellwalker method and their spectra are modelled as modified blackbody functions. According to IRAS and Planck data, most fields have dust colour temperatures T_C ~ 14-18K and opacity spectral index values of beta=1.5-1.9. The clumps/cores identified in SCUBA-2 maps have T~ 13K and similar beta values. There are some indications of the dust emission spectrum becoming flatter at wavelengths longer than 500um. In fits involving Planck data, the significance is limited by the uncertainty of the corrections for CO line contamination. The fits to the SPIRE data give a median beta value slightly above 1.8. In the joint SPIRE and SCUBA-2 850um fits the value decreases to beta ~1.6. Most of the observed T-beta anticorrelation can be explained by noise. The typical submillimetre opacity spectral index beta of cold clumps is found to be ~1.7. This is above the values of diffuse clouds but lower than in some previous studies of dense clumps. There is only tentative evidence of T-beta anticorrelation and beta decreasing at millimetre wavelengths.

The low dust temperatures (<14 K) of Planck Galactic Cold Clumps (PGCCs) make them ideal targets to probe the initial conditions and very early phase of star formation. TOP-SCOPE is a joint survey program targeting ~2000 PGCCs in J=1-0 transitions of CO isotopologues and ~1000 PGCCs in 850 micron continuum emisison. The objective of the TOP-SCOPE survey and the joint surveys (SMT 10-m, KVN 21-m and NRO 45-m) is to statistically study the initial conditions occurring during star formation and the evolution of molecular clouds, across a wide range of environments. The observations, data analysis and example science cases for these surveys are introduced with an exemplar source, PGCC G26.53+0.17 (G26), which is a filamentary infrared dark cloud (IRDC). The total mass, the length and the mean line-mass (M/L) of the G26 filament are ~6200 Msun, ~12 pc and ~500 Msun/pc, respectively. Ten massive clumps including eight starless ones are found along the filament. The most massive Clump as a whole may be still in global collapse while its denser part seems to be undergoing expansion due to outflow feedback. The fragmentation in G26 filament from cloud scale to clump scale is in agreement with gravitational fragmentation of an isothermal, non-magnetized, and turbulent supported cylinder. A bimodal behavior in dust emissivity spectral index ($\beta$) distribution is found in G26, suggesting grain growth along the filament. The G26 filament may be formed due to large-scale compression flows evidenced by the temperature and velocity gradients across its natal cloud.

Central engine of gamma-ray bursts (GRBs) may be intermittent and launch several episodes of ejecta separated by a long quiescent interval. In this scenario, an external shock is formed due to the propagation of the first launched ejecta into the circum-burst medium and the later launched ejecta may interact with the external shock at later period. Owing to the internal dissipation, the later launched ejecta may be observed at a later time ($t_{\rm{jet}}$). In this paper, we study the relation of $t_{\rm{b}}$ and $t_{\rm{jet}}$, where $t_{\rm{b}}$ is the collision time of the later launched ejecta with the formed external shock. It is found that the relation of $t_{\rm{b}}$ and $t_{\rm{jet}}$ depends on the bulk Lorentz factor ($\Gamma_{\rm{jet}}$) of the later launched ejecta and the density ($\rho$) of the circum-burst medium. If the value of $\Gamma_{\rm{jet}}$ or $\rho$ is low, the $t_{\rm{b}}$ would be significantly larger than $t_{\rm{jet}}$. However, the $t_{\rm{b}}\sim t_{\rm{jet}}$ can be found if the value of $\Gamma_{\rm{jet}}$ or $\rho$ is significantly large. Our results can explain the large lag of the optical emission relative to the $\gamma$-ray/X-ray emission in GRBs, e.g., GRB~111209A. For GRBs with a precursor, our results suggest that the energy injection into the external shock and thus more than one external-reverse shock may appear in the main prompt emission phase. According to our model, we estimate the Lorentz factor of the second launched ejecta in GRB~160625B.

Long-duration gamma-ray bursts (LGRBs) and ultra-LGRBs (ULGRBs) originate from collapsars, in the center of which a newborn rotating stellar-mass black hole (BH) surrounded by a massive accretion disk may form. In the scenario of BH hyperaccretion inflow-outflow model and Blandford-Znajek (BZ) mechanism to trigger gamma-ray bursts (GRBs), the real accretion rate to power a BZ jet is far lower than the mass supply rate from the progenitor star. The characteristics of the progenitor stars can be constrained by GRB luminosity observations, and the results exceed usual expectations. LGRBs lasting from several seconds to tens of seconds in the rest frame may originate from solar-metallicity ($Z \sim 1~ Z_\odot$, where $Z$ and $Z_\odot$ are the metallicities of progenitor stars and the Sun), massive ($M \gtrsim 34 ~M_\odot$, where $M$ and $M_\odot$ are the masses of progenitor stars and the Sun) stars or some zero-metallicity ($Z \sim 0$) stars. A fraction of low-metallicity ($Z \lesssim 10^{-2}~Z_\odot$) stars, including Population III stars, can produce ULGRBs such as GRB 111209A. The fraction of LGRBs lasting less than tens of seconds in the rest frame is more than 40$\%$, which cannot conform to the fraction of the demanded type of progenitor star. It possibly implies that the activity timescale of central engine may be much longer than the observed timescale of prompt emission phase, as indicated by X-ray late-time activities. Alternatively, LGRBs and ULGRBs may be powered by a millisecond magnetar central engine.

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.

We present analysis of the first Chandra observation of PSO J334.2028+01.4075 (PSO J334), targeted as a binary-AGN candidate based on periodic variations of the optical flux. With no prior targeted X-ray coverage for PSO J334, our new 40 ksec Chandra observation allows for the opportunity to differentiate between a single or binary-AGN system, and if a binary, can characterize the mode of accretion. Simulations show that the two expected accretion disk morphologies for binary-AGN systems are (i) a "cavity", where the inner region of the accretion disk is mostly empty and emission is truncated blueward of the wavelength associated with the temperature of the innermost ring, or (ii) "minidisks", where there is substantial accretion from the cirumbinary disk onto one or both of the members of the binary, each with their own shock-heated thin-disk accretion system. We find the X-ray emission to be well-fit with an absorbed power-law, incompatible with the simple cavity scenario. Further, we construct an SED of PSO J334 by combining radio through X-ray observations and find that the SED agrees well with that of a normal AGN, most likely incompatible with the minidisk scenario. Other analyses, such as locating the quasar on IR color-color diagrams and analyzing the quasar mass predicted by the fundamental plane of black hole activity, further highlight the similarity of PSO J334 with respect to normal AGN. On the multi-wavelength fronts we investigated, we find no evidence supporting PSO J334 as a binary-AGN system, though our analysis remains insensitive to some binary configurations.

Aims. We present a multiwavelength study of two southern Galactic H II regions G346.056-0.021 and G346.077-0.056 which are located at a distance of 10.9 kpc. The distribution of ionized gas, cold and warm dust and the stellar population associated with the two H II regions are studied in detail using measurements at near-infrared, mid-infrared, far-infrared, submillimeter and radio wavelengths. Methods. The radio continuum maps at 1280 and 610 MHz were obtained using the Giant Metrewave Radio Telescope to probe the ionized gas. The dust temperature, column density and dust emissivity maps were generated by using modified blackbody fits in the far-infrared wavelength range 160 - 500 {\mu}m. Various near- and mid-infrared colour and magnitude criteria were adopted to identify candidate ionizing star(s) and the population of young stellar objects in the associated field. Results. The radio maps reveal the presence diffuse ionized emission displaying distinct cometary morphologies. The 1280 MHz flux densities translate to ZAMS spectral types in the range O7.5V - O7V and O8.5V - O8V for the ionizing stars of G346.056-0.021 and G346.077-0.056, respectively. A few promising candidate ionizing star(s) are identified using near-infrared photometric data. The column density map shows the presence of a large, dense dust clump enveloping G346.077-0.056. The dust temperature map shows peaks towards the two H II regions. The submillimetre image shows the presence of two additional clumps one being associated with G346.056-0.021. The masses of the clumps are estimated to range between {\sim} 1400 to 15250 M{\sun}. Based on simple analytic calculations and the correlation seen between the ionized gas distribution and the local density structure, the observed cometary morphology in the radio maps is better explained invoking the champagne-flow model.

We systematically analyze X-ray variability of active galactic nuclei (AGNs) in the 7~Ms \textit{Chandra} Deep Field-South survey. On the longest timescale ($\approx~17$ years), we find only weak (if any) dependence of X-ray variability amplitudes on energy bands or obscuration. We use four different power spectral density (PSD) models to fit the anti-correlation between normalized excess variance ($\sigma^2_{\rm nxv}$) and luminosity, and obtain a best-fit power law index $\beta=1.16^{+0.05}_{-0.05}$ for the low-frequency part of AGN PSD. We also divide the whole light curves into 4 epochs in order to inspect the dependence of $\sigma^2_{\rm nxv}$ on these timescales, finding an overall increasing trend. The analysis of these shorter light curves also infers a $\beta$ of $\sim 1.3$ that is consistent with the above-derived $\beta$, which is larger than the frequently-assumed value of $\beta=1$. We then investigate the evolution of $\sigma^2_{\rm nxv}$. No definitive conclusion is reached due to limited source statistics but, if present, the observed trend goes in the direction of decreasing AGN variability at fixed luminosity toward large redshifts. We also search for transient events and find 6 notable candidate events with our considered criteria. Two of them may be a new type of fast transient events, one of which is reported here for the first time. We therefore estimate a rate of fast outbursts $\langle\dot{N}\rangle = 1.0^{+1.1}_{-0.7}\times 10^{-3}~\rm galaxy^{-1}~yr^{-1}$ and a tidal disruption event~(TDE) rate $\langle\dot{N}_{\rm TDE}\rangle=8.6^{+8.5}_{-4.9}\times 10^{-5}~\rm galaxy^{-1}~yr^{-1}$ assuming the other four long outbursts to be TDEs.

Neutrino-dominated accretion flow (NDAF) around a rotating stellar-mass black hole (BH) is one of the plausible candidates for the central engines of gamma-ray bursts (GRBs). Because the time-variant and anisotropic emission of neutrinos from NDAFs leads to GRB variability, NDAFs can be regarded as the sources of the strong gravitational waves (GWs). We calculate the dependences of the GW strains on both the BH spin and the accretion rate. We demonstrate that for typical GRBs with either single pulse or multiple pulses, the GWs from NDAFs might be detected at a distance of $\sim 100$ kpc/$\sim 1$ Mpc by advanced LIGO/Einstein Telescope with a typical frequency of $\sim 10-100$ Hz. Besides NDAFs, the other two competitive candidates for GRB central engine are Blandford-Znajek (BZ) mechanism and millisecond magnetars. We explore the GW signals from these two as well, and compare the corresponding results with NDAFs'. We find that for a certain GRB, the possible detected distance from NDAFs is about two orders of magnitude higher than that from BZ mechanism, but at least two orders of magnitude lower than that from magnetars. The typical GW frequency for BZ mechanism is the same with that of NDAFs, $\sim 10-100$ Hz, while the typical frequency for magnetars is $\sim 2000$ Hz. Therefore, the GWs released by the central engines of adjacent GRBs might be used to determine whether there is an NDAF, a BZ jet or a magnetar in GRB center.

We investigate the population of high-redshift ($3\leq z < 6$) AGN selected in the two deepest X-ray surveys, the 7 Ms \textit{Chandra} Deep Field-South and 2 Ms \textit{Chandra} Deep Field-North. Their outstanding sensitivity and spectral characterization of faint sources allow us to focus on the sub-$L_*$ regime (log$L_{\mathrm{X}}\lesssim44$), poorly sampled by previous works using shallower data, and the obscured population. Taking fully into account the individual photometric-redshift probability distribution functions, the final sample consists of $\approx102$ X-ray selected AGN at $3\leq z < 6$. The fraction of AGN obscured by column densities log$N_{\mathrm{H}}>23$ is $\sim0.6-0.8$, once incompleteness effects are taken into account, with no strong dependence on redshift or luminosity. We derived the high-redshift AGN number counts down to $F_{\mathrm{0.5-2\,keV}}=7\times10^{-18}\,\mathrm{erg\,cm^{-2}\,s^{-1}}$, extending previous results to fainter fluxes, especially at $z>4$. We put the tightest constraints to date on the low-luminosity end of AGN luminosity function at high redshift. The space-density, in particular, declines at $z>3$ at all luminosities, with only a marginally steeper slope for low-luminosity AGN. By comparing the evolution of the AGN and galaxy densities, we suggest that such a decline at high luminosities is mainly driven by the underlying galaxy population, while at low luminosities there are hints of an intrinsic evolution of the parameters driving nuclear activity. Also, the black-hole accretion rate density and star-formation rate density, which are usually found to evolve similarly at $z\lesssim3$, appear to diverge at higher redshifts.

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.

The Antarctic Impulsive Transient Antenna (ANITA) is a NASA long-duration balloon experiment with the primary goal of detecting ultra-high-energy ($>10^{18}\,\mbox{eV}$) neutrinos via the Askaryan Effect. The fourth ANITA mission, ANITA-IV, recently flew from Dec 2 to Dec 29, 2016. For the first time, the Tunable Universal Filter Frontend (TUFF) boards were deployed for mitigation of narrow-band, anthropogenic noise with tunable, switchable notch filters. The TUFF boards also performed second-stage amplification by approximately 45 dB to boost the $\sim\,\mu\mbox{V-level}$ radio frequency (RF) signals to $\sim$ mV-level for digitization, and supplied power via bias tees to the first-stage, antenna-mounted amplifiers. The other major change in signal processing in ANITA-IV is the resurrection of the $90^{\circ}$ hybrids deployed previously in ANITA-I, in the trigger system, although in this paper we focus on the TUFF boards. During the ANITA-IV mission, the TUFF boards were successfully operated throughout the flight. They contributed to a factor of 2.8 higher total instrument livetime on average in ANITA-IV compared to ANITA-III due to reduction of narrow-band, anthropogenic noise before a trigger decision is made.

A multiwavelength investigation of the southern infrared dust bubble CS51 is presented in this paper. We probe the associated ionized, cold dust, molecular and stellar components. Radio continuum emission mapped at 610 and 1300 MHz, using the Giant Metrewave Radio Telescope, India, reveal the presence of three compact emission components (A, B, and C) apart from large-scale diffuse emission within the bubble interior. Radio spectral index map show the coexistence of thermal and non-thermal emission components. Modified blackbody fits to the thermal dust emission using Herschel PACS and SPIRE data is performed to generate dust temperature and column density maps. We identify five dust clumps associated with CS51 with masses and radius in the range 810 - 4600 M{\sun} and 1.0 - 1.9 pc, respectively. We further construct the column density probability distribution functions of the surrounding cold dust which display the impact of ionization feedback from high-mass stars. The estimated dynamical and fragmentation timescales indicate the possibility of collect and collapse mechanism in play at the bubble border. Molecular line emission from the MALT90 survey is used to understand the nature of two clumps which show signatures of expansion of CS51.

We study the X-ray variability properties of distant AGNs in the Chandra Deep Field-South region over 17 years, up to $z\sim 4$, and compare them with those predicted by models based on local samples. We use the results of Monte Carlo simulations to account for the biases introduced by the discontinuous sampling and the low-count regime. We confirm that variability is an ubiquitous property of AGNs, with no clear dependence on the density of the environment. The variability properties of high-z AGNs, over different temporal timescales, are most consistent with a Power Spectral Density (PSD) described by a broken (or bending) power-law, similar to nearby AGNs. We confirm the presence of an anti-correlation between luminosity and variability, resulting from the dependence of variability on BH mass and accretion rate. We explore different models, finding that our acceptable solutions predict that BH mass influences the value of the PSD break frequency, while the Eddington ratio $\lambda_{Edd}$ affects the PSD break frequency and, possibly, the PSD amplitude as well. We derive the evolution of the average $\lambda_{Edd}$ as a function of redshift, finding results in agreement with measurements based on different estimators. The large statistical uncertainties make our results consistent with a constant Eddington ratio, although one of our models suggest a possible increase of $\lambda_{Edd}$ with lookback time up to $z\sim 2-3$. We conclude that variability is a viable mean to trace the accretion history of supermassive BHs, whose usefulness will increase with future, wide-field/large effective area X-ray missions.

Neutrino-dominated accretion flows (NDAFs) around rotating stellar-mass black holes (BHs) are plausible candidates for the central engines of gamma-ray bursts (GRBs). NDAFs are hyperaccretion disks with accretion rates in the range of around 0.001-10 $M_\odot ~\rm s^{-1}$, which have high density and temperature and therefore are extremely optically thick and geometrically slim or even thick. We review the theoretical progresses in studying the properties of NDAFs as well as their applications to the GRB phenomenology. The topics include: the steady radial and vertical structure of NDAFs and the implications for calculating neutrino luminosity and annihilation luminosity, jet power due to neutrino-antineutrino annihilation and Blandford-Znajek mechanism and their dependences on parameters such as BH mass, spin, and accretion rate, time evolution of NDAFs, effect of magnetic fields, applications of NDAF theories to the GRB phenomenology such as lightcurve variability, extended emission, X-ray flares, kilonovae, etc., as well as probing NDAFs using multi-messenger signals such as MeV neutrinos and gravitational waves.

We present the results from a series of ground-based radio observations toward a Planck Galactic Cold Clump (PGCC), PGCC G108.84-00.81, which is located in one curved filamentary cloud in the vicinity of an extended HII region Sh2-152 and SNR G109.1-1.0. PGCC G108.84-00.81 is mainly composed of two clumps, "G108-N" and "G108-S". In the 850 micron dust continuum emission map, G108-N is shown as one component while G108-S is fragmented into four components. There is no infrared source associated with G108-N while there are two infrared sources (IRS 1 and IRS 2) associated with G108-S. The total mass of G108-N is larger than the jeans mass, suggesting that G108-N is gravitationally unstable and a potential place for a future star formation. The clump properties of G108-N and G108-S such as the gas temperature and the column density, are not distinctly different. However, G108-S is slightly more evolved than G108-N, in the consideration of the CO depletion factor, molecular abundances, and association with infrared sources. G108-S seems to be affected by the compression from Sh2-152, while G108-N is relatively protected from the external effect

The primary science goal of the NASA-sponsored ANITA project is measurement of ultra-high energy neutrinos and cosmic rays, observed via radio-frequency signals resulting from a neutrino- or cosmic ray- interaction with terrestrial matter (atmospheric or ice molecules, e.g.). Accurate inference of the energies of these cosmic rays requires understanding the transmission/reflection of radio wave signals across the ice-air boundary. Satellite-based measurements of Antarctic surface reflectivity, using a co-located transmitter and receiver, have been performed more-or-less continuously for the last few decades. Satellite-based reflectivity surveys, at frequencies ranging from 2--45 GHz and at near-normal incidence, yield generally consistent reflectivity maps across Antarctica. Using the Sun as an RF source, and the ANITA-3 balloon borne radio-frequency antenna array as the RF receiver, we have also measured the surface reflectivity over the interval 200-1000 MHz, at elevation angles of 12-30 degrees, finding agreement with the Fresnel equations within systematic errors. To probe low incidence angles, inaccessible to the Antarctic Solar technique and not probed by previous satellite surveys, a novel experimental approach ("HiCal-1") was devised. Unlike previous measurements, HiCal-ANITA constitute a bi-static transmitter-receiver pair separated by hundreds of kilometers. Data taken with HiCal, between 200--600 MHz shows a significant departure from the Fresnel equations, constant with frequency over that band, with the deficit increasing with obliquity of incidence, which we attribute to the combined effects of possible surface roughness, surface grain effects, radar clutter and/or shadowing of the reflection zone due to Earth curvature effects.

We present a detailed spectral analysis of the brightest Active Galactic Nuclei (AGN) identified in the 7Ms Chandra Deep Field South (CDF-S) survey over a time span of 16 years. Using a model of an intrinsically absorbed power-law plus reflection, with possible soft excess and narrow Fe K$\alpha$ line, we perform a systematic X-ray spectral analysis, both on the total 7Ms exposure and in four different periods with lengths of 2-21 months. With this approach, we not only present the power-law slopes, column densities $N_H$, observed fluxes, and absorption-corrected 2-10~keV luminosities $L_X$ for our sample of AGNs, but also identify significant spectral variabilities among them on time scales of years. We find that the $N_H$ variabilities can be ascribed to two different types of mechanisms, either flux-driven or flux-independent. We also find that the correlation between the narrow Fe line EW and $N_H$ can be well explained by the continuum suppression with increasing $N_H$. Accounting for the sample incompleteness and bias, we measure the intrinsic distribution of $N_H$ for the CDF-S AGN population and present re-selected subsamples which are complete with respect to $N_H$. The $N_H$-complete subsamples enable us to decouple the dependences of $N_H$ on $L_X$ and on redshift. Combining our data with that from C-COSMOS, we confirm the anti-correlation between the average $N_H$ and $L_X$ of AGN, and find a significant increase of the AGN obscured fraction with redshift at any luminosity. The obscured fraction can be described as $f_{obscured}\thickapprox 0.42\ (1+z)^{0.60}$.

Gravitational radiation is an excellent field for testing theories of gravity in strong gravitational fields. The current observations on the gravitational-wave (GW) bursts by LIGO have already placed various constraints on the alternative theories of gravity. In this paper, we investigate the possible bounds which could be placed on the Brans-Dicke gravity using GW detection from inspiralling compact binaries with the proposed Einstein Telescope, a third-generation GW detector. We first calculate in details the waveforms of gravitational radiation in the lowest post-Newtonian approximation, including the tensor and scalar fields, which can be divided into the three polarization modes, i.e. "plus mode", "cross mode" and "breathing mode". Applying the stationary phase approximation, we obtain their Fourier transforms, and derive the correction terms in amplitude, phase and polarization of GWs, relative to the corresponding results in General Relativity. Imposing the noise level of Einstein Telescope, we find that the GW detection from inspiralling compact binaries, composed of a neutron star and a black hole, can place stringent constraints on the Brans-Dicke gravity. The bound on the coupling constant $\omega_{\rm BD}$ depends on the mass, sky-position, orbital angle, polarization angle, luminosity distance, redshift distribution and total observed number $N_{\rm GW}$ of the binary systems. Taking into account all the burst events up to redshift $z=5$, we find that the bound could be $\omega_{\rm BD}\gtrsim 10^{6}\times(N_{\rm GW}/10^4)^{1/2}$. Even for the conservative estimation with $10^{4}$ observed events, the bound is still more than one order tighter than the current limit from Solar System experiments. So, we conclude that Einstein Telescope will provide a powerful platform to test alternative theories of gravity.

We investigate the dependence of black-hole accretion rate (BHAR) on host-galaxy star formation rate (SFR) and stellar mass ($M_*$) in the CANDELS/GOODS-South field in the redshift range of $0.5\leq z < 2.0$. Our sample consists of $\approx 18000$ galaxies, allowing us to probe galaxies with $0.1 \lesssim \mathrm{SFR} \lesssim 100\ M_\odot\ \mathrm{yr}^{-1}$ and/or $10^8 \lesssim M_* \lesssim 10^{11}\ M_{\odot}$. We use sample-mean BHAR to approximate long-term average BHAR. Our sample-mean BHARs are derived from the $Chandra$ Deep Field-South 7 Ms observations, while the SFRs and $M_*$ have been estimated by the CANDELS team through SED fitting. The average BHAR is correlated positively with both SFR and $M_*$, and the BHAR-SFR and BHAR-$M_*$ relations can both be described acceptably by linear models with a slope of unity. However, BHAR appears to be correlated more strongly with $M_*$ than SFR. This result indicates that $M_*$ is the primary host-galaxy property related to black-hole growth, and the apparent BHAR-SFR relation is largely a secondary effect due to the star-forming main sequence. Among our sources, massive galaxies ($M_* \gtrsim 10^{10} M_{\odot}$) have significantly higher BHAR/SFR ratios than less-massive galaxies, indicating the former have higher black-hole fueling efficiency and/or higher SMBH occupation fraction than the latter. Our results can naturally explain the observed proportionality between $M_{\rm BH}$ and $M_*$ for local giant ellipticals, and suggest their $M_{\rm BH}/M_*$ is higher than that of local star-forming galaxies. Among local star-forming galaxies, massive systems might have higher $M_{\rm BH}/M_*$ compared to dwarfs.

Screened modified gravity (SMG) is a kind of scalar-tensor theory with screening mechanisms, which can suppress the fifth force in dense regions and allow theories to evade the solar system and laboratory tests. In this paper, we investigate how the screening mechanisms in SMG affect the gravitational radiation damping effects, calculate in detail the rate of the energy loss due to the emission of tensor and scalar gravitational radiations, and derive their contributions to the change in the orbital period of the binary system. We find that the scalar radiation depends on the screened parameters and the propagation speed of scalar waves, and the scalar dipole radiation dominates the orbital decay of the binary system. For strongly self-gravitating bodies, all effects of scalar sector are strongly suppressed by the screening mechanisms in SMG. By comparing our results to observations of binary system PSR J1738+0333, we place the stringent constraints on the screening mechanisms in SMG. As an application of these results, we focus on three specific models of SMG (chameleon, symmetron, and dilaton), and derive the constraints on the model parameters, respectively.

A wide bandwidth, dual polarized, modified four-square antenna is presented as a feed antenna for radio astronomical measurements. A linear array of these antennas is used as a line-feed for cylindrical reflectors for Tianlai, a radio interferometer designed for 21~cm intensity mapping. Simulations of the feed antenna beam patterns and scattering parameters are compared to experimental results at multiple frequencies across the 650 - 1420 MHz range. Simulations of the beam patterns of the combined feed array/reflector are presented as well.

We report a sample of 463 high-mass starless clump (HMSC) candidates within $-60\deg<l<60\deg$ and $-1\deg<b<1\deg$. This sample has been singled out from 10861 ATLASGAL clumps. All of these sources are not associated with any known star-forming activities collected in SIMBAD and young stellar objects identified using color-based criteria. We also make sure that the HMSC candidates have neither point sources at 24 and 70 \micron~nor strong extended emission at 24 $\mu$m. Most of the identified HMSCs are infrared ($\le24$ $\mu$m) dark and some are even dark at 70 $\mu$m. Their distribution shows crowding in Galactic spiral arms and toward the Galactic center and some well-known star-forming complexes. Many HMSCs are associated with large-scale filaments. Some basic parameters were attained from column density and dust temperature maps constructed via fitting far-infrared and submillimeter continuum data to modified blackbodies. The HMSC candidates have sizes, masses, and densities similar to clumps associated with Class II methanol masers and HII regions, suggesting they will evolve into star-forming clumps. More than 90% of the HMSC candidates have densities above some proposed thresholds for forming high-mass stars. With dust temperatures and luminosity-to-mass ratios significantly lower than that for star-forming sources, the HMSC candidates are externally heated and genuinely at very early stages of high-mass star formation. Twenty sources with equivalent radius $r_\mathrm{eq}<0.15$ pc and mass surface density $\Sigma>0.08$ g cm$^{-2}$ could be possible high-mass starless cores. Further investigations toward these HMSCs would undoubtedly shed light on comprehensively understanding the birth of high-mass stars.

We present the first results from the B-fields In STar-forming Region Observations (BISTRO) survey, using the Sub-millimetre Common-User Bolometer Array 2 (SCUBA-2) camera, with its associated polarimeter (POL-2), on the James Clerk Maxwell Telescope (JCMT) in Hawaii. We discuss the survey's aims and objectives. We describe the rationale behind the survey, and the questions which the survey will aim to answer. The most important of these is the role of magnetic fields in the star formation process on the scale of individual filaments and cores in dense regions. We describe the data acquisition and reduction processes for POL-2, demonstrating both repeatability and consistency with previous data. We present a first-look analysis of the first results from the BISTRO survey in the OMC 1 region. We see that the magnetic field lies approximately perpendicular to the famous 'integral filament' in the densest regions of that filament. Furthermore, we see an 'hour-glass' magnetic field morphology extending beyond the densest region of the integral filament into the less-dense surrounding material, and discuss possible causes for this. We also discuss the more complex morphology seen along the Orion Bar region. We examine the morphology of the field along the lower-density north-eastern filament. We find consistency with previous theoretical models that predict magnetic fields lying parallel to low-density, non-self-gravitating filaments, and perpendicular to higher-density, self-gravitating filaments.

We derive a simple analytical formula to describe the evolution of spectral index $\beta$ in the steep decay phase shaped by the curvature effect with assumption that the spectral parameters and Lorentz factor of jet shell is the same for different latitude. Here, the value of $\beta$ is estimated in 0.3$-$10keV energy band. For a spherical thin shell with a cutoff power law (CPL) intrinsic radiation spectrum, the spectral evolution can be read as a linear function of observer time. For the situation with Band function intrinsic radiation spectrum, the spectral evolution may be complex. If the observed break energy of radiation spectrum is larger than 10keV, the spectral evolution is the same as that shaped by jet shells with a CPL spectrum. If the observed break energy is less than 0.3keV, the value of $\beta$ would be a constant. Others, the spectral evolution can be approximated as a logarithmal function of the observer time in generally.

The curvature effect may be responsible for the steep decay phase observed in gamma-ray bursts. For testing the curvature effect with observations, the zero time point $t_0$ adopted to plot observer time and flux on a logarithmic scale should be appropriately selected. In practice, however, the true $t_0$ cannot be directly constrained from the data. Then, we move $t_0$ to a certain time in the steep decay phase, which can be easily identified. In this situation, we derive an analytical formula to describe the flux evolution of the steep decay phase. The analytical formula is read as $F_\nu\propto (1+\tilde t_{\rm obs}/{\tilde t_c})^{-\alpha}$ with $\alpha(\tilde{t}_{\rm obs})=2+{\int_{0}^{\log (1+\tilde{t}_{\rm obs}/{\tilde{t}_c})} {\beta(\tau)d[\log(1+\tau/{\tilde{t}_c})]}}/{\log (1 + {\tilde t}_{\rm obs}/{{\tilde t}_c})}$, where $F_\nu$ is the flux observed at frequency $\nu$, $\tilde t_{\rm obs}$ is the observer time by setting zero time point $t_0$ at a certain time in the steep decay phase, $\beta$ is the spectral index estimated around $\nu$, and ${\tilde t}_c$ is the decay timescale of the phase with $\tilde{t}_{\rm obs}{\geqslant}0$. We test the analytical formula with the data from numerical calculations. It is found that the analytical formula presents a well estimation about the evolution of flux shaped by the curvature effect. Our analytical formula can be used to confront the curvature effect with observations and estimate the decay timescale of the steep decay phase.

We report the detection of morphology dependent stellar age in massive quenched galaxies (QGs) at z~1.2. The sense of the dependence is that compact QGs are 0.5-2 Gyr older than normal-sized ones. The evidence comes from three different age indicators, Dn4000, Hdelta and fits to spectral synthesis models, applied to their stacked optical spectra. All age indicators consistently show that the stellar populations of compact QGs are older than their normally-sized counterparts. We detect weak [OII] emission in a fraction of QGs, and the strength of the line, when present, is similar between the two samples; however, compact galaxies exhibit significantly lower frequency of [OII] emission than normal ones. A fraction of both samples are individually detected in 7 Ms Chandra X-ray images (luminosities$\sim10^{40}$-$10^{41}$ erg/sec). 7 Ms stacks of non-detected galaxies show similarly low luminosities in the soft band only, consistent with a hot gas origin for the X-ray emission. While both [OII] emitters and non-emitters are also X-ray sources among normal galaxies, no compact galaxy with [OII] emission is an X-ray source, arguing against an AGN powering the line in compact galaxies. We interpret the [OII] properties as further evidence that compact galaxies are older and further along into the process of quenching star-formation and suppressing gas accretion. Finally, we argue that the older age of compact QGs is evidence of progenitor bias: compact QGs simply reflect the smaller sizes of galaxies at their earlier quenching epoch, with stellar density most likely having nothing directly to do with cessation of star-formation.

Recent simulations on super-Eddington accretion flows have shown that, apart from the diffusion process, the vertical advection based on magnetic buoyancy can be a more efficient process to release the trapped photons in the optically thick disk. As a consequence, the radiative luminosity from the accretion disk can be far beyond the Eddington value. Following this spirit, we revisit the structure and radiation of hyper-accretion disks with mass accretion rates in the range $10^{-3}\sim 10~M_{\sun}~{\rm s}^{-1}$. Our results show that, due to the strong cooling through the vertical advection, the disk temperature becomes lower than that in the classic model without the vertical advection process, and therefore the neutrino luminosity from the disk is lower. On the other hand, the gamma-ray photons released through the vertical advection can be extremely super-Eddington. We argue that the large amount of escaped gamma-ray photons may have more significant contribution to the primordial fireball than the neutrino annihilation, and may hint a link between gamma-ray bursts and kilonovae in the black hole hyper-accretion scenario.

We report on a search for ultra-high-energy (UHE) neutrinos from gamma-ray bursts (GRBs) in the data set collected by the Testbed station of the Askaryan Radio Array (ARA) in 2011 and 2012. From 57 selected GRBs, we observed no events that survive our cuts, which is consistent with 0.12 expected background events. Using NeuCosmA as a numerical GRB reference emission model, we estimate upper limits on the prompt UHE GRB neutrino fluence and quasi-diffuse flux from $10^{7}$ to $10^{10}$ GeV. This is the first limit on the prompt UHE GRB neutrino quasi-diffuse flux above $10^{7}$ GeV.

Cosmological phase transitions (CPTs), such as the Grand Unified Theory (GUT) and the electroweak (EW) ones, play a significant role in both particle physics and cosmology. In this letter, we propose to probe the first-order CPTs, by detecting gravitational waves (GWs) which are generated during the phase transitions through the cosmic microwave background (CMB). If happened around the inflation era, the first-order CPTs may yield low-frequency GWs due to bubble dynamics, leaving imprints on the CMB. In contrast to the nearly scale-invariant primordial GWs caused by vacuum fluctuation, these bubble-generated GWs are scale dependent and have non-trivial B-mode spectra. If decoupled from inflaton, the EWPT during inflation may serve as a probe for the one after reheating where the baryon asymmetry could be generated via EW baryogenesis (EWBG). The CMB thus provides a potential way to test the feasibility of the EWBG, complementary to the collider measurements of Higgs potential and the direct detection of GWs generated during EWPT.

We present X-ray source catalogs for the $\approx7$ Ms exposure of the Chandra Deep Field-South (CDF-S), which covers a total area of 484.2 arcmin$^2$. Utilizing WAVDETECT for initial source detection and ACIS Extract for photometric extraction and significance assessment, we create a main source catalog containing 1008 sources that are detected in up to three X-ray bands: 0.5-7.0 keV, 0.5-2.0 keV, and 2-7 keV. A supplementary source catalog is also provided including 47 lower-significance sources that have bright ($K_s\le23$) near-infrared counterparts. We identify multiwavelength counterparts for 992 (98.4%) of the main-catalog sources, and we collect redshifts for 986 of these sources, including 653 spectroscopic redshifts and 333 photometric redshifts. Based on the X-ray and multiwavelength properties, we identify 711 active galactic nuclei (AGNs) from the main-catalog sources. Compared to the previous $\approx4$ Ms CDF-S catalogs, 291 of the main-catalog sources are new detections. We have achieved unprecedented X-ray sensitivity with average flux limits over the central $\approx1$ arcmin$^2$ region of $\approx1.9\times10^{-17}$, $6.4\times10^{-18}$, and $2.7\times10^{-17}$ erg cm$^{-2}$ s$^{-1}$ in the three X-ray bands, respectively. We provide cumulative number-count measurements observing, for the first time, that normal galaxies start to dominate the X-ray source population at the faintest 0.5-2.0 keV flux levels. The highest X-ray source density reaches $\approx50\,500$ deg$^{-2}$, and $47\%\pm4\%$ of these sources are AGNs ($\approx23\,900$ deg$^{-2}$).

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.

We observed thirteen Planck cold clumps with the James Clerk Maxwell Telescope/SCUBA-2 and with the Nobeyama 45 m radio telescope. The N$_2$H$^+$ distribution obtained with the Nobeyama telescope is quite similar to SCUBA-2 dust distribution. The 82 GHz HC$_3$N, 82 GHz CCS, and 94 GHz CCS emission are often distributed differently with respect to the N$_2$H$^+$ emission. The CCS emission, which is known to be abundant in starless molecular cloud cores, is often very clumpy in the observed targets. We made deep single-pointing observations in DNC, HN$^{13}$C, N$_2$D$^+$, cyclic-C$_3$H$_2$ toward nine clumps. The detection rate of N$_2$D$^+$ is 50\%. Furthermore, we observed the NH$_3$ emission toward 15 Planck cold clumps to estimate the kinetic temperature, and confirmed that most of targets are cold ($\lesssim$ 20 K). In two of the starless clumps observe, the CCS emission is distributed as it surrounds the N$_2$H$^+$ core (chemically evolved gas), which resembles the case of L1544, a prestellar core showing collapse. In addition, we detected both DNC and N$_2$D$^+$. These two clumps are most likely on the verge of star formation. We introduce the Chemical Evolution Factor (CEF) for starless cores to describe the chemical evolutionary stage, and analyze the observed Planck cold clumps.

We present a systematic search for periodically varying quasars and supermassive black hole binary (SMBHB) candidates in the Pan-STARRS1 (PS1) Medium Deep Survey's MD09 field. From a color-selected sample of 670 quasars extracted from a multi-band deep-stack catalog of point sources, we locally select variable quasars and look for coherent periods with the Lomb-Scargle periodogram. Three candidates from our sample demonstrate strong variability for more than ~3 cycles, and their PS1 light curves are well fitted to sinusoidal functions. We test the persistence of the candidates' apparent periodic variations detected during the 4.2 years of the PS1 survey with archival photometric data from the SDSS Stripe 82 survey or new monitoring with the Large Monolithic Imager at the Discovery Channel Telescope. None of the three periodic candidates (including PSO J334.2028+1.4075) remain persistent over the extended baseline of 7 - 14 years, corresponding to a detection rate of < 1 in 670 quasars in a search area of 5 deg^2. Even though SMBHBs should be a common product of the hierarchal growth of galaxies, and periodic variability in SMBHBs has been theoretically predicted, a systematic search for such signatures in a large optical survey is strongly limited by its temporal baseline and the "red noise" associated with normal quasar variability. We show that follow-up long-term monitoring (>5 cycles) is crucial to our search for these systems.

A supramassive, strongly-magnetized millisecond neutron star (NS) has been proposed to be the candidate central engine of at least some short gamma-ray bursts (SGRBs), based on the "internal plateau" commonly observed in the early X-ray afterglow. While a previous analysis shows a qualitative consistency between this suggestion and the Swift SGRB data, the distribution of observed break time $t_b$ is much narrower than the distribution of the collapse time of supramassive NSs for the several NS equations-of-state (EoSs) investigated. In this paper, we study four recently-constructed "unified" NS EoSs, as well as three developed strange quark star (QS) EoSs within the new confinement density-dependent mass model. All the EoSs chosen here satisfy the recent observational constraints of the two massive pulsars whose masses are precisely measured. We construct sequences of rigidly rotating NS/QS configurations with increasing spinning frequency $f$, from non-rotating ($f = 0$) to the Keplerian frequency ($f = f_{\rm K}$), and provide convenient analytical parametrizations of the results. Assuming that the cosmological NS-NS merger systems have the same mass distribution as the Galactic NS-NS systems, we demonstrate that all except the BCPM NS EoS can reproduce the current $22\%$ supramassive NS/QS fraction constraint as derived from the SGRB data. We simultaneously simulate the observed quantities (the break time $t_b$, the break time luminosity $L_b$ and the total energy in the electromagnetic channel $E_{\rm total}$) of SGRBs, and find that while equally well reproducing other observational constraints, QS EoSs predict a much narrower $t_b$ distribution than that of the NS EoSs, better matching the data. We therefore suggest that the post-merger product of NS-NS mergers might be fast-rotating supramassive QSs rather than NSs.

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

We perform long-term ($\approx 15$ yr, observed-frame) X-ray variability analyses of the 68 brightest radio-quiet active galactic nuclei (AGNs) in the 6 Ms $Chandra$ Deep Field-South (CDF-S) survey; the majority are in the redshift range of $0.6-3.1$, providing access to penetrating rest-frame X-rays up to $\approx 10-30$ keV. Twenty-four of the 68 sources are optical spectral type I AGNs, and the rest (44) are type II AGNs. The time scales probed in this work are among the longest for X-ray variability studies of distant AGNs. Photometric analyses reveal widespread photon-flux variability: $90\%$ of AGNs are variable above a 95% confidence level, including many X-ray obscured AGNs and several optically classified type II quasars. We characterize the intrinsic X-ray luminosity ($L_{\rm{X}}$) and absorption ($N_{\rm{H}}$) variability via spectral fitting. Most (74%) sources show $L_{\rm{X}}$ variability; the variability amplitudes are generally smaller for quasars. A Compton-thick candidate AGN shows variability of its high-energy X-ray flux, indicating the size of reflecting material to be $\lesssim 0.3$ pc. $L_{\rm{X}}$ variability is also detected in a broad absorption line (BAL) quasar. The $N_{\rm{H}}$ variability amplitude for our sample appears to rise as time separation increases. About 16% of sources show $N_{\rm{H}}$ variability. One source transitions from an X-ray unobscured to obscured state while its optical classification remains type I; this behavior indicates the X-ray eclipsing material is not large enough to obscure the whole broad-line region.

We observed 146 Galactic clumps in HCN (4-3) and CS (7-6) with the Atacama Submillimeter Telescope Experiment (ASTE) 10-m telescope. A tight linear relationship between star formation rate and gas mass traced by dust continuum emission was found for both Galactic clumps and the high redshift (z>1) star forming galaxies (SFGs), indicating a constant gas depletion time of ~100 Myr for molecular gas in both Galactic clumps and high z SFGs. However, low z galaxies do not follow this relation and seem to have a longer global gas depletion time. The correlations between total infrared luminosities (L_TIR) and molecular line luminosities (L'_mol) of HCN (4-3) and CS (7-6) are tight and sublinear extending down to clumps with LTIR~10^{3} L_sun. These correlations become linear when extended to external galaxies. A bimodal behavior in the LTIR--L'mol correlations was found for clumps with different dust temperature, luminosity-to-mass ratio, and sigmaline/sigmavir. Such bimodal behavior may be due to evolutionary effects. The slopes of LTIR--L'mol correlations become more shallow as clumps evolve. We compared our results with lower J transition lines in wu et al. (2010). The correlations between clump masses and line luminosities are close to linear for low effective excitation density tracers but become sublinear for high effective excitation density tracers for clumps with LTIR larger than LTIR~10^4.5 Lsun. High effective excitation density tracers cannot linearly trace the total clump masses, leading to a sublinear correlations for both Mclump-L'mol and LTIR-L'mol relations.

HLX-1, currently the best intermediate-mass black hole candidate, has undergone seven violent outbursts, each with a peak X-ray luminosity of $L_{\mathrm{peak},\mathrm{X}}\sim 10^{42}\ \rm{erg\ s^{-1}}$. Interestingly, the properties of the HLX-1 outbursts evolve with time. In this work, we aim to constrain the physical parameters of the central engine of the HLX-1 outbursts in the framework of the black hole accretion. We find that the physical properties of the HLX-1 outbursts are consistent with being driven by the radiation pressure instability. This scenario can explain the evolution of the recurrent timescales of the HLX-1 outbursts as a function of the durations.

The emission from black hole binaries (BHBs) and active galactic nuclei (AGNs) displays significant aperiodic variabilities. The most promising explanation for these variabilities is the propagating fluctuations in the accretion flow. It is natural to expect that the mechanism driving variabilities in BHBs and AGNs may operate in a black hole hyper-accretion disk, which is believed to power gamma-ray bursts (GRBs). We study the variabilities of jet power in GRBs based on the model of propagating fluctuations. It is found that the variabilities of jet power and the temporal profile of erratic spikes in this scenario are similar to those in observed light curves of prompt gamma-ray emission of GRBs. Our results show that the mechanism driving X-ray variabilities in BHBs and AGNs may operate in the central engine to drive the variabilities of GRBs.

We report on four radio-detected cosmic-ray (CR) or CR-like events observed with the Antarctic Impulsive Transient Antenna (ANITA), a NASA-sponsored long-duration balloon payload. Two of the four were previously identified as stratospheric CR air showers during the ANITA-I flight. A third stratospheric CR was detected during the ANITA-II flight. Here we report on characteristics these three unusual CR events, which develop nearly horizontally, 20-30~km above the surface of the Earth. In addition, we report on a fourth steeply upward-pointing ANITA-I CR-like radio event which has characteristics consistent with a primary that emerged from the surface of the ice. This suggests a possible $\tau$-lepton decay as the origin of this event, but such an interpretation would require significant suppression of the Standard Model ${\tau}$-neutrino cross section.

The neutrino mass hierarchy is one of the neutrino fundamental properties yet to be determined. We introduce a method to determine neutrino mass hierarchy by comparing the interaction rate of neutral current (NC) interactions, $\nu(\hat{\nu}) + p\rightarrow\nu(\hat{\nu}) + p$, and inverse beta decays (IBD), $\bar{\nu}_e + p\rightarrow n + e^+$, of supernova neutrinos in scintillation detectors. Neutrino flavor conversions inside the supernova are sensitive to neutrino mass hierarchy. Due to Mikheyev-Smirnov-Wolfenstein effects, the full swapping of $\bar{\nu}_e$ flux with the $\bar{\nu}_x$ ($x=\mu,~\tau$) one occurs in the inverted hierarchy, while such a swapping does not occur in the normal hierarchy. As a result, more high energy IBD events occur in the detector for the inverted hierarchy than the high energy IBD events in the normal hierarchy. By comparing IBD interaction rate with the mass hierarchy independent NC interaction rate, one can determine the neutrino mass hierarchy.

Neutrino-dominated accretion flows (NDAFs) around rotating stellar-mass black holes (BHs) have been theorized as the central engine of relativistic jets launched in massive star core collapse events or compact star mergers. In this work, we calculate the electron neutrino/anti-neutrino spectra of NDAFs by fully taking into account the general relativistic effects, and investigate the effects of viewing angle, BH spin, and mass accretion rate on the results. We show that even though a typical NDAF has a neutrino luminosity lower than that of a typical supernova (SN), it can reach $10^{50}-10^{51}~{\rm erg~s^{-1}}$ peaking at $\sim 10$ MeV, making them potentially detectable with the upcoming sensitive MeV neutrino detectors if they are close enough to Earth. Based on the observed GRB event rate in the local universe and requiring that at least 3 neutrinos are detected to claim a detection, we estimate a detection rate up to $\sim$ (0.10-0.25) per century for GRB-related NDAFs by the Hyper-Kamiokande (Hyper-K) detector if one neglects neutrino oscillation. If one assumes that all Type Ib/c SNe have an engine-driven NDAF, the Hyper-K detection rate would be $\sim$ (1-3) per century. By considering neutrino oscillations, the detection rate may decrease by a factor of 2-3. Detecting one such event would establish the observational evidence of NDAFs in the universe.

We present a multi-wavelength analysis of the host galaxy of short-duration gamma-ray burst (GRB) 150101B. Follow-up optical and X-ray observations suggested that the host galaxy, 2MASX J12320498-1056010, likely harbors a low-luminosity active galactic nuclei (AGN). Our modeling of the spectral energy distribution (SED) has confirmed the nature of the AGN, making it the first reported GRB host that contains an AGN. We have also found the host galaxy is a massive elliptical galaxy with stellar population of $\sim 5.7\ Gyr$, one of the oldest among the short-duration GRB hosts. Our analysis suggests that the host galaxy can be classified as an X-ray bright, optically normal galaxy (XBONG), and the central AGN is likely dominated by a radiatively inefficient accretion flow (RIAF). Our work explores interesting connection that may exist between GRB and AGN activities of the host galaxy, which can help understand the host environment of the GRB events and the roles of AGN feedback.

The maximal $U(1)_L$ supersymmetric inverse seesaw mechanism (M$L$SIS) provides a natural way to relate asymmetric dark matter (ADM) with neutrino physics. In this paper we point out that, M$L$SIS is a natural outcome if one dynamically realizes the inverse seesaw mechanism in the next-to minimal supersymmetric standard model (NMSSM) via the dimension-five operator $(N)^2S^2/M_*$, with $S$ the NMSSM singlet developing TeV scale VEV; it slightly violates lepton number due to the suppression by the fundamental scale $M_*$, thus preserving $U(1)_L$ maximally. The resulting sneutrino is a distinguishable ADM candidate, oscillating and favored to have weak scale mass. A fairly large annihilating cross section of such a heavy ADM is available due to the presence of singlet.

Fast radio bursts (FRBs) are radio transients lasting only about a few milliseconds. They seem to occur at cosmological distances. We propose that these events can be originated in the collapse of the magnetosphere of Kerr-Newman black holes (KNBHs). We show that the closed orbits of charged particles in the magnetosphere of these objects are unstable. After examining the dependencies on the specific charge of the particle and the spin and charge of the KNBH, we conclude that the resulting timescale and radiation mechanism fit well with the extant observations of FRBs. Furthermore, we argue that the merger of a KNBH binary is one of the plausible central engines for potential gamma-ray or radio afterglow following a certain FRBs, and can also account for gravitational wave (GW) events like GW 150914. Our model leads to predictions that can be tested by combined multi-wavelength electromagnetic and GW observations.

Ultra-high energy neutrinos are interesting messenger particles since, if detected, they can transmit exclusive information about ultra-high energy processes in the Universe. These particles, with energies above $10^{16}\mathrm{eV}$, interact very rarely. Therefore, detectors that instrument several gigatons of matter are needed to discover them. The ARA detector is currently being constructed at South Pole. It is designed to use the Askaryan effect, the emission of radio waves from neutrino-induced cascades in the South Pole ice, to detect neutrino interactions at very high energies. With antennas distributed among 37 widely-separated stations in the ice, such interactions can be observed in a volume of several hundred cubic kilometers. Currently 3 deep ARA stations are deployed in the ice of which two have been taking data since the beginning of the year 2013. In this publication, the ARA detector "as-built" and calibrations are described. Furthermore, the data reduction methods used to distinguish the rare radio signals from overwhelming backgrounds of thermal and anthropogenic origin are presented. Using data from only two stations over a short exposure time of 10 months, a neutrino flux limit of $3 \cdot 10^{-6} \mathrm{GeV} / (\mathrm{cm^2 \ s \ sr})$ is calculated for a particle energy of 10^{18}eV, which offers promise for the full ARA detector.

We propose a compact binary model for the fast radio burst (FRB) repeaters, where the system consists of a magnetic white dwarf (WD) and a neutron star (NS) with strong bipolar magnetic fields. When the WD fills its Roche lobe, mass transfer will occur from the WD to the NS through the inner Lagrange point. The accreted magnetized materials may trigger magnetic reconnection when they approach the NS surface, and therefore the electrons can be accelerated to an ultra-relativistic speed. In this scenario, the curvature radiation of the electrons moving along the NS magnetic field lines can account for the characteristic frequency and the timescale of an FRB. Owing to the conservation of angular momentum, the WD may be kicked away after a burst, and the next burst may appear when the system becomes semi-detached again through the gravitational radiation. By comparing our analyses with the observations, we show that such an intermittent Roche lobe overflow mechanism can be responsible for the observed repeating behavior of FRB 121102.

For fifty years, cosmic-ray air showers have been detected by their radio emission. We present the first laboratory measurements that validate electrodynamics simulations used in air shower modeling. An experiment at SLAC provides a beam test of radio-frequency (RF) radiation from charged particle cascades in the presence of a magnetic field, a model system of a cosmic-ray air shower. This experiment provides a suite of controlled laboratory measurements to compare to particle-level simulations of RF emission, which are relied upon in ultra-high-energy cosmic-ray air shower detection. We compare simulations to data for intensity, linearity with magnetic field, angular distribution, polarization, and spectral content. In particular, we confirm modern predictions that the magnetically induced emission in a dielectric forms a cone that peaks at the Cherenkov angle and show that the simulations reproduce the data within systematic uncertainties.

Most bipolar outflows are associated with individual young stellar objects and have small opening angles. Here we report the discovery of an extremely wide-angle ($\sim$180$\arcdeg$) bipolar outflow ("EWBO") in a cluster forming region AFGL 5142 from low-velocity emission of the HCN (3-2) and HCO$^{+}$ (3-2) lines. This bipolar outflow is along a north-west to south-east direction with a line-of-sight flow velocity of about 3 km~s$^{-1}$ and is spatially connected to the high-velocity jet-like outflows. It seems to be a collection of low-velocity material entrained by the high-velocity outflows due to momentum feedback. The total ejected mass and mass loss rate due to both high velocity jet-like outflows and the "EWBO" are $\sim$24.5 M$_{\sun}$ and $\sim1.7\times10^{-3}$ M$_{\sun}$~yr$^{-1}$, respectively. Global collapse of the clump is revealed by the "blue profile" in the HCO$^{+}$ (1-0) line. A hierarchical network of filaments was identified in NH$_{3}$ (1,1) emission. Clear velocity gradients of the order of 10 km~s$^{-1}$~pc$^{-1}$ are found along filaments, indicating gas inflow along the filaments. The sum of the accretion rate along filaments and mass infall rate along the line of sight is $\sim$3.1$\times10^{-3}$ M$_{\sun}$~yr$^{-1}$, which exceeds the total mass loss rate, indicating that the central cluster is probably still gaining mass. The central cluster is highly fragmented and 22 condensations are identified in 1.1 mm continuum emission. The fragmentation process seems to be determined by thermal pressure and turbulence. The magnetic field may not play an important role in fragmentation.

We study the effects of an asymmetric radiation field on the properties of a molecular cloud envelope. We employ observations of carbon monoxide (12CO and 13CO), atomic carbon, ionized carbon, and atomic hydrogen to analyze the chemical and physical properties of the core and envelope of L1599B, a molecular cloud forming a portion of the ring at approximately 27 pc from the star Lambda Ori. The O III star provides an asymmetric radiation field that produces a moderate enhancement of the external radiation field. Observations of the [CII] fine structure line with the GREAT instrument on SOFIA indicate a significant enhanced emission on the side of the cloud facing the star, while the [Ci], 12CO and 13CO J = 1-0 and 2-1, and 12CO J = 3-2 data from the PMO and APEX telescopes suggest a relatively typical cloud interior. The atomic, ionic, and molecular line centroid velocities track each other very closely, and indicate that the cloud may be undergoing differential radial motion. The HI data from the Arecibo GALFA survey and the SOFIA/GREAT [CII] data do not suggest any systematic motion of the halo gas, relative to the dense central portion of the cloud traced by 12CO and 13CO.

Black holes (BHs) hide themselves behind various astronomical phenomena, and their properties, i.e., mass and spin, are usually difficult to constrain. One leading candidate for the central engine model of gamma-ray bursts (GRBs) invokes a stellar mass BH and a neutrino-dominated accretion flow (NDAF), with the relativistic jet launched due to neutrino-anti-neutrino annihilations. Such a model gives rise to a matter-dominated fireball, and is suitable to interpret GRBs with a dominant thermal component with a photospheric origin. We propose a method to constrain BH mass and spin within the framework of this model, and apply the method to a thermally-dominant GRB 101219B whose initial jet launching radius $r_0$ is constrained from the data. Using our numerical model of NDAF jets, we estimate the following constraints on the central BH: mass $M_{\rm BH} \sim 5-9~M_\odot$, spin parameter $a_* \gtrsim 0.6$, and disk mass $3~M_\odot \lesssim M_{\rm disk} \lesssim 4~M_\odot$. Our results also suggest that the NDAF model is a competitive candidate for the central engine of GRBs with a strong thermal component.

Long-duration gamma-ray bursts (LGRBs) are generally considered to originate from the massive collapsars. It is believed that the central engine of gamma-ray bursts (GRBs) is a neutrino-dominated accretion flow (NDAF) around a rotating stellar-mass black hole (BH). The neutrino annihilation above the NDAF is a feasible mechanism to power GRB. In this work, we analyse the distributions of the isotropic gamma-ray radiated energy and jet kinetic energy of 48 LGRBs. According to the NDAF and fireball models, we estimate the mean accreted masses of LGRBs in our sample to investigate whether the NDAFs can power LGRBs with the reasonable BH parameters and conversion efficiency of neutrino annihilation. The results indicate that most of the values of the accreted masses are less than $5~M_\odot$ for the extreme Kerr BHs and high conversion efficiency. It suggests that the NDAFs may be suitable for most of LGRBs except for some extremely high energy sources.

The physical mechanisms that induce the transformation of a certain mass of gas in new stars are far from being well understood. Infrared bubbles associated with HII regions have been considered to be good samples of investigating triggered star formation. In this paper we report on the investigation of the dust properties of the infrared bubble N4 around the HII region G11.898+0.747, analyzing its interaction with its surroundings and star formation histories therein, with the aim of determining the possibility of star formation triggered by the expansion of the bubble. Using Herschel PACS and SPIRE images with a wide wavelength coverage, we reveal the dust properties over the entire bubble. Meanwhile, we are able to identify six dust clumps surrounding the bubble, with a mean size of 0.50 pc, temperature of about 22 K, mean column density of 1.7 $\times10^{22}$ cm$^{-2}$, mean volume density of about 4.4 $\times10^{4}$ cm$^{-3}$, and a mean mass of 320 $M_{\odot}$. In addition, from PAH emission seen at 8 $\mu$m, free-free emission detected at 20 cm and a probability density function in special regions, we could identify clear signatures of the influence of the HII region on the surroundings. There are hints of star formation, though further investigation is required to demonstrate that N4 is the triggering source.

We present the first survey of dense gas towards Planck Galactic Cold Clumps (PGCCs). Observations in the J=1-0 transitions of HCO+ and HCN towards 621 molecular cores associated with PGCCs were performed using the Purple Mountain Observatory 13.7-m telescope. Among them, 250 sources have detection, including 230 cores detected in HCO+ and 158 in HCN. Spectra of the J=1-0 transitions from CO, 13CO, and C18O at the centers of the 250 cores were extracted from previous mapping observations to construct a multi-line data set. The significantly low detection rate of asymmetric double-peaked profiles, together with the well consistence among central velocities of CO, HCO+, and HCN spectra, suggests that the CO-selected Planck cores are more quiescent compared to classical star-forming regions. The small difference between line widths of C18O and HCN indicates that the inner regions of CO-selected Planck cores are not more turbulent than the exterior. The velocity-integrated intensities and abundances of HCO+ are positively correlated with those of HCN, suggesting these two species are well coupled and chemically connected. The detected abundances of both HCO+ and HCN are significantly lower than values in other low- to high-mass star-forming regions. The low abundances may be due to beam dilution. On the basis of the inspection of the parameters given in the PGCC catalog, we suggest that there may be about 1 000 PGCC objects having sufficient reservoir of dense gas to form stars.

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

We are performing a series of observations with ground-based telescopes toward Planck Galactic cold clumps (PGCCs) in the $\lambda$ Orionis complex in order to systematically investigate the effects of stellar feedback. In the particular case of PGCC G192.32-11.88, we discovered an extremely young Class 0 protostellar object (G192N) and a proto-brown dwarf candidate (G192S). G192N and G192S are located in a gravitationally bound bright-rimmed clump. The velocity and temperature gradients seen in line emission of CO isotopologues indicate that PGCC G192.32-11.88 is externally heated and compressed. G192N probably has the lowest bolometric luminosity ($\sim0.8$ L$_{\sun}$) and accretion rate (6.3$\times10^{-7}$ M$_{\sun}$~yr$^{-1}$) when compared with other young Class 0 sources (e.g. PACS Bright Red sources (PBRs)) in the Orion complex. It has slightly larger internal luminosity ($0.21\pm0.01$ L$_{\sun}$) and outflow velocity ($\sim$14 km~s$^{-1}$) than the predictions of first hydrostatic cores (FHSCs). G192N might be among the youngest Class 0 sources, which are slightly more evolved than a FHSC. Considering its low internal luminosity ($0.08\pm0.01$ L$_{\odot}$) and accretion rate (2.8$\times10^{-8}$ M$_{\sun}$~yr$^{-1}$), G192S is an ideal proto-brown dwarf candidate. The star formation efficiency ($\sim$0.3\%-0.4\%) and core formation efficiency ($\sim$1\%) in PGCC G192.32-11.88 are significantly smaller than in other giant molecular clouds or filaments, indicating that the star formation therein is greatly suppressed due to stellar feedback.

We present Submillimeter Array (SMA) observations of the HCN\,(3--2) and HCO$^+$\,(3--2) molecular lines toward the W3(H$_2$O) and W3(OH) star-forming complexes. Infall and outflow motions in the W3(H$_2$O) have been characterized by observing HCN and HCO$^+$ transitions. High-velocity blue/red-shifted emission, tracing the outflow, show multiple knots, which might originate in episodic and precessing outflows. `Blue-peaked' line profiles indicate that gas is infalling onto the W3(H$_2$O) dust core. The measured large mass accretion rate, 2.3$\times$10$^{-3}$~M$_{\odot}$~yr$^{-1}$, together with the small free-fall time scale, 5$\times$10$^{3}$~yr, suggest W3(H$_2$O) is in an early evolutionary stage of the process of formation of high-mass stars. For the W3(OH), a two-layer model fit to the HCN and HCO$^+$ spectral lines and Spizter/IRAC images support that the W3(OH) H{\sc ii} region is expanding and interacting with the ambient gas, with the shocked neutral gas being expanding with an expansion timescale of 6.4$\times$10$^{3}$~yr. The observations suggest different kinematical timescales and dynamical states for the W3(H$_2$O) and W3(OH).

We present Submillimeter Array (SMA) molecular line observations in two 2 GHz-wide bands centered at 217.5 and 227.5 GHz, toward the massive star forming region W51 North. We identified 84 molecular line transitions from 17 species and their isotopologues. The molecular gas distribution of these lines mainly peaks in the continuum position of W51 North, and has a small tail extending to the west, probably associated with W51 d2. In addition to the commonly detected nitrogen and oxygen-bearing species, we detected a large amount of transitions of the Acetone (CH$_3$COCH$_3$) and Methyl Formate (CH$_3$OCHO), which may suggest that these molecules are present in an early evolutionary stage of the massive stars. We also found that W51 North is an ethanol-rich source. There is no obvious difference in the molecular gas distributions between the oxygen-bearing and nitrogen-bearing molecules. Under the assumption of Local Thermodynamic Equilibrium (LTE), with the XCLASS tool, the molecular column densities, and rotation temperatures are estimated. We found that the oxygen-bearing molecules have considerable higher column densities and fractional abundances than the nitrogen-bearing molecules. The rotation temperatures range from 100 to 200 K, suggesting that the molecular emission could be originated from a warm environment. Finally, based on the gas distributions, fractional abundances and the rotation temperatures, we conclude that CH$_3$OH, C$_2$H$_5$OH, CH$_3$COCH$_3$ and CH$_3$CH$_2$CN might be synthesized on the grain surface, while gas phase chemistry is responsible for the production of CH$_3$OCH$_3$, CH$_3$OCHO and CH$_2$CHCN.

CIT 6 is a carbon star in the transitional phase from the asymptotic giant branch (AGB) to the protoplanetary nebulae (pPN). Observational evidences of two point sources in the optical, circumstellar arc segments in an HC$_3$N line emission, and a bipolar nebula in near-infrared provide strong support for the presence of a binary companion. Hence, CIT 6 is very attractive for studying the role of companions in the AGB-pPN transition. We have carried out high resolution $^{12}$CO $J=2-1$ and $^{13}$CO $J=2-1$ observations of CIT 6 with the Submillimeter Array combined with the Submillimeter Telescope (single-dish) data. The $^{12}$CO channel maps reveal a spiral-shell pattern connecting the HC$_3$N segments in a continuous form, and an asymmetric outflow corresponding to the near-infrared bipolar nebula. Rotation of the $^{12}$CO channel peak position may be related to the inner spiral winding and/or the bipolar outflow. An eccentric orbit binary is suggested for the presences of an anisotropic mass loss to the west and a double spiral pattern. The lack of interarm emission to the west may indicate a feature corresponding to the periastron passage of a highly eccentric orbit of the binary. Spatially-averaged radial and spectral profiles of $^{12}$CO $J=2-1$ and $^{13}$CO $J=2-1$ are compared with simple spherical radiative transfer models, suggesting a change of $^{12}$CO/$^{13}$CO abundance ratio from $\sim30$ to $\sim50$ inward in the CSE of CIT 6. The millimeter continuum emission is decomposed into extended dust thermal emission (spectral index $\sim-2.4$) and compact emission from radio photosphere (spectral index $\sim-2.0$).

W49N is a mini-starburst in the Milky Way and thus an ideal laboratory for high-mass star formation studies. Due to its large distance (11.1$_{-0.7}^{+0.9}$ kpc), the kinematics inside and between the dense molecular clumps in W49N are far from well understood. The SMA observations resolved the continuum emission into two clumps. The molecular line observation of SO$_{2}$ (28$_{4,24}$-28$_{3,25}$) suggests that the two clumps have a velocity difference of $\sim$7 km~s$^{-1}$. The eastern clump is very close to two radio sources "G1" and "G2", and the western clump coincides with a radio source "B". The HCN (3-2) line reveals an extremely energetic outflow, which is among the most energetic molecular outflows in the Milky Way. This is the first report of high-velocity molecular outflow detection in W49N. The outflow jet might be in precession, which could account for the distribution, velocity and rotation of water maser spots. Three absorption systems are identified in HCO$^{+}$ (3-2) spectra. The absorption features are blueshifted with respect to the emission of SO$_{2}$ (28$_{4,24}$-28$_{3,25}$) lines, indicating that a cold layer is expanding in front of the warm gas. Further analysis indicates that the expansion is decelerated from the geometric expansion centers.

The Ultra-Fast Flash Observatory pathfinder (UFFO-p) is a new space mission dedicated to detect Gamma-Ray Bursts (GRBs) and rapidly follow their afterglows in order to provide early optical/ultraviolet measurements. A GRB location is determined in a few seconds by the UFFO Burst Alert & Trigger telescope (UBAT) employing the coded mask imaging technique and the detector combination of Yttrium Oxyorthosilicate (YSO) scintillating crystals and multi-anode photomultiplier tubes. The results of the laboratory tests of UBAT's functionality and performance are described in this article. The detector setting, the pixel-to-pixel response to X-rays of different energies, the imaging capability for <50 keV X-rays, the localization accuracy measurements, and the combined test with the Block for X-ray and Gamma-Radiation Detection (BDRG) scintillator detector to check the efficiency of UBAT are all described. The UBAT instrument has been assembled and integrated with other equipment on UFFO-p and should be launched on board the Lomonosov satellite in late-2015.

Fe and Al are two of the most important rock-forming elements other than Mg, Si, and O. Their presence in the lower mantle's most abundant minerals, MgSiO_3 bridgmanite, MgSiO_3 post-perovskite and MgO periclase, alters their elastic properties. However, knowledge on the thermoelasticity of Fe- and Al-bearing MgSiO_3 bridgmanite, and post-perovskite is scarce. In this study, we perform ab initio molecular dynamics to calculate the elastic and seismic properties of pure, Fe^{3+}- and Fe^{2+}-, and Al^{3+}-bearing MgSiO_3 perovskite and post-perovskite, over a wide range of pressures, temperatures, and Fe/Al compositions. Our results show that a mineral assemblage resembling pyrolite fits a 1D seismological model well, down to, at least, a few hundred kilometers above the core-mantle boundary, i.e. the top of the D'' region. In D'', a similar composition is still an excellent fit to the average velocities and fairly approximate to the density. We also implement polycrystal plasticity with a geodynamic model to predict resulting seismic anisotropy, and find post-perovskite with predominant (001) slip across all compositions agrees best with seismic observations in the D''.

The Askaryan Radio Array (ARA) is an ultra-high energy ($>10^{17}$ eV) cosmic neutrino detector in phased construction near the South Pole. ARA searches for radio Cherenkov emission from particle cascades induced by neutrino interactions in the ice using radio frequency antennas ($\sim150-800$ MHz) deployed at a design depth of 200 m in the Antarctic ice. A prototype ARA Testbed station was deployed at $\sim30$ m depth in the 2010-2011 season and the first three full ARA stations were deployed in the 2011-2012 and 2012-2013 seasons. We present the first neutrino search with ARA using data taken in 2011 and 2012 with the ARA Testbed and the resulting constraints on the neutrino flux from $10^{17}-10^{21}$ eV.

Supermassive black hole binaries (SMBHBs) should be an inevitable consequence of the hierarchical growth of massive galaxies through mergers, and the strongest sirens of gravitational waves (GWs) in the cosmos. And yet, their direct detection has remained elusive due to the compact (sub-parsec) orbital separations of gravitationally bound SMBHBs. Here we exploit a theoretically predicted signature of a SMBHB in the time domain: periodic variability caused by a mass accretion rate that is modulated by the binary's orbital motion. We report our first significant periodically varying quasar detection from the systematic search in the Pan-STARRS1 (PS1) Medium Deep Survey. Our SMBHB candidate, PSO J334.2028+01.4075, is a luminous radio-loud quasar at $z=2.060$, with extended baseline photometry from the Catalina Real-Time Transient Survey, as well as archival spectroscopy from the FIRST Bright Quasar Survey. The observed period ($542 \pm 15$ days) and estimated black hole mass ($\log (M_{\rm BH}/M_\odot) = 9.97 \pm 0.50$), correspond to an orbital separation of $7^{+8}_{-4}$ Schwarzschild radii ($\sim 0.006^{+0.007}_{-0.003}$ pc), assuming the rest-frame period of the quasar variability traces the orbital period of the binary. This SMBHB candidate, discovered at the peak redshift for SMBH mergers, is in a physically stable configuration for a circumbinary accretion disk, and within the regime of GW-driven orbital decay. Our search with PS1 is a benchmark study for the exciting capabilities of LSST, which will have orders of magnitude larger survey power, and will potentially pinpoint the locations of thousands of SMBHBs in the variable night sky.

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

A neutrino-dominated accretion flow (NDAF) around a rotating stellar-mass black hole (BH) is one of the plausible candidates for the central engine of gamma-ray bursts (GRBs). Two mechanisms, i.e., Blandford-Znajek (BZ) mechanism and neutrino annihilation process, are generally considered to power GRBs. Using the analytic solutions from Xue et al. (2013) and ignoring the effects of the magnetic field configuration, we estimate the BZ and neutrino annihilation luminosities as the functions of the disk masses and BH spin parameters to contrast the observational jet luminosities of GRBs. The results show that, although the neutrino annihilation processes could account for most of GRBs, the BZ mechanism is more effective, especially for long-duration GRBs. Actually, if the energy of afterglows and flares of GRBs is included, the distinction between these two mechanisms is more significant. Furthermore, massive disk mass and high BH spin are beneficial to power high luminosities of GRBs. Finally, we discuss possible physical mechanisms to enhance the disk mass or the neutrino emission rate of NDAFs and relevant difference between these two mechanisms.

Stellar-mass black holes (BHs) surrounded by neutrino-dominated accretion flows (NDAFs) are the plausible candidates to power gamma-ray bursts (GRBs) via neutrinos emission and their annihilation. The progenitors of short-duration GRBs (SGRBs) are generally considered to be compact binaries mergers. According to the simulation results, the disk mass of the NDAF has been limited after merger events. We can estimate such disk mass by using the current SGRB observational data and fireball model. The results show that the disk mass of a certain SGRB mainly depends on its output energy, jet opening angle, and central BH characteristics. Even for the extreme BH parameters, some SGRBs require massive disks, which approach or exceed the limits in simulations. We suggest that there may exist alternative magnetohydrodynamic processes or some mechanisms increasing the neutrino emission to produce SGRBs with the reasonable BH parameters and disk mass.

We present the effects of the vertical convection on the structure and luminosity of the neutrino-dominated accretion flow (NDAF) around a stellar-mass black hole in spherical coordinates. We found that the convective energy transfer can suppress the radial advection in the NDAF, and that the density, temperature and opening angle are slightly changed. As a result, the neutrino luminosity and annihilation luminosity are increased, which is conducive to achieve the energy requirement of gamma-ray bursts.

We present the Swift X-ray Cluster Survey (SWXCS) catalog obtained using archival data from the X-ray telescope (XRT) on board the Swift satellite acquired from 2005 to 2012, extending the first release of the SWXCS. The catalog provides positions, soft fluxes, and, when possible, optical counterparts for a flux-limited sample of X-ray group and cluster candidates. We consider the fields with Galactic latitude |b| > 20 degree to avoid high HI column densities. We discard all of the observations targeted at groups or clusters of galaxies, as well as particular extragalactic fields not suitable to search for faint extended sources. We finally select ~3000 useful fields covering a total solid angle of ~400 degree^2. We identify extended source candidates in the soft-band (0.5-2keV) images of these fields using the software EXSdetect, which is specifically calibrated for the XRT data. Extensive simulations are used to evaluate contamination and completeness as a function of the source signal, allowing us to minimize the number of spurious detections and to robustly assess the selection function. Our catalog includes 263 candidate galaxy clusters and groups down to a flux limit of 7E-15 erg/cm^2/s in the soft band, and the logN-logS is in very good agreement with previous deep X-ray surveys. The final list of sources is cross-correlated with published optical, X-ray, and SZ catalogs of clusters. We find that 137 sources have been previously identified as clusters, while 126 are new detections. Currently, we have collected redshift information for 158 sources (60% of the entire sample). Once the optical follow-up and the X-ray spectral analysis of the sources are complete, the SWXCS will provide a large and well-defined catalog of groups and clusters of galaxies to perform statistical studies of cluster properties and tests of cosmological models.

We report on the Submillimeter Array (SMA) observations of molecular lines at 270 GHz toward W3(OH) and W3(H$_2$O) complex. Although previous observations already resolved the W3(H$_2$O) into two or three sub-components, the physical and chemical properties of the two sources are not well constrained. Our SMA observations clearly resolved W3(OH) and W3(H$_2$O) continuum cores. Taking the advantage of the line fitting tool XCLASS, we identified and modeled a rich molecular spectrum in this complex, including multiple CH$_3$CN and CH$_3$OH transitions in both cores. HDO, C$_2$H$_5$CN, O$^{13}$CS, and vibrationally excited lines of HCN, CH$_3$CN, and CH$_3$OCHO were only detected in W3(H$_2$O). We calculate gas temperatures and column densities for both cores. The results show that W3(H$_{2}$O) has higher gas temperatures, and larger column densities than W3(OH) as previously observed, suggesting physical and chemical differences between the two cores. We compare the molecular abundances in W3(H$_2$O) to those in the Sgr B2(N) hot core, the Orion KL hot core and the Orion Compact Ridge, and discuss the chemical origin of specific species. An east-west velocity gradient is seen in W3(H$_2$O), and the extension is consistent with the bipolar outflow orientation traced by water masers and radio jets. A north-south velocity gradient across W3(OH) is also observed. However, with current observations we can not assure if the velocity gradients are caused by rotation, outflow or radial velocity differences of the sub-components in W3(OH).

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

In the center of active galactic nuclei (AGN), the dusty torus absorb the radiation from the central engine and re-emit in mid-infrared (MIR). Observations have detected moderate anisotropy in the dust MIR emission, in the way that type 1 AGNs (type1s) are mildly brighter in MIR comparing with type 2 sources (type2s). However, type1s and type2s were found to follow statistically the same tight MIR -- hard X-ray correlation, suggesting the MIR emission is highly isotropic assuming the hard X-ray radiation is inclination independent. We argue this discrepancy could be solved considering the hard X-ray emission in AGN is also mildly anisotropic as we recently discovered. To verify this diagram, we compare the sub-arcsecond 12\mu m flux densities of type1s and type2s using [OIV]$\lambda$25.89\mu m emission line as an isotropic luminosity indicator. We find that on average type1s are brighter in nuclei 12\mu m radiation by a factor of $2.6 \pm 0.6$ than type2s at given [OIV]$\lambda$25.89\mu m luminosities, confirming the mild anisotropy of the nuclei 12\mu m emission. We show that the anisotropy of the 12\mu m emission we detected is in good agreement with radiative transfer models of clumpy torus. The fact that type 1 and type 2 AGNs follow the same tight MIR -- hard X-ray correlation instead supports that both the MIR and hard X-ray emission in AGNs are mildly anisotropic.

We present a multi-wavelength study of the IR bubble G24.136+00.436. The J=1-0 observations of $^{12}$CO, $^{13}$CO and C$^{18}$O were carried out with the Purple Mountain Observatory 13.7 m telescope. Molecular gas with a velocity of 94.8 km s$^{-1}$ is found prominently in the southeast of the bubble, shaping as a shell with a total mass of $\sim2\times10^{4}$ $M_{\odot}$. It is likely assembled during the expansion of the bubble. The expanding shell consists of six dense cores. Their dense (a few of $10^{3}$ cm$^{-3}$) and massive (a few of $10^{3}$ $M_{\odot}$) characteristics coupled with the broad linewidths ($>$ 2.5 km s$^{-1}$) suggest they are promising sites of forming high-mass stars or clusters. This could be further consolidated by the detection of compact HII regions in Cores A and E. We tentatively identified and classified 63 candidate YSOs based on the \emph{Spitzer} and UKIDSS data. They are found to be dominantly distributed in regions with strong emission of molecular gas, indicative of active star formation especially in the shell. The HII region inside the bubble is mainly ionized by a $\sim$O8V star(s), of the dynamical age $\sim$1.6 Myr. The enhanced number of candidate YSOs and secondary star formation in the shell as well as time scales involved, indicate a possible scenario of triggering star formation, signified by the "collect and collapse" process.

How high-mass stars form remains unclear currently. Calculation suggests that the radiation pressure of a forming star can halt spherical infall, preventing its further growth when it reaches 10 M$_{\odot}$. Two major theoretical models on the further growth of stellar mass were proposed. One model suggests the mergence of less massive stellar objects, and the other is still through accretion but with the help of disk. Inflow motions are the key evidence of how forming stars further gain mass to build up massive stars. Recent development in technology has boosted the search of inflow motion. A number of high-mass collapse candidates were obtained with single dish observations, mostly showed blue profile. The infalling signatures seem to be more common in regions with developed radiation pressure than in younger cores, which opposes the theoretical prediction and is also very different from that of low mass star formation. Interferometer studies so far confirm such tendency with more obvious blue profile or inverse P Cygni profile. Results seem to favor the accretion model. However, the evolution tendency of the infall motion in massive star forming cores needs to be further explored. Direct evidence for monolithic or competitive collapse processes is still lack. ALMA will enable us to probe more detail of gravity process.

The physical and chemical properties of prestellar cores, especially massive ones, are still far from being well understood due to the lack of a large sample. The low dust temperature ($<$14 K) of Planck cold clumps makes them promising candidates for prestellar objects or for sources at the very initial stages of protostellar collapse. We have been conducting a series of observations toward Planck cold clumps (PCCs) with ground-based radio telescopes. In general, when compared with other star forming samples (e.g. infrared dark clouds), PCCs are more quiescent, suggesting that most of them may be in the earliest phase of star formation. However, some PCCs are associated with protostars and molecular outflows, indicating that not all PCCs are in a prestellar phase. We have identified hundreds of starless dense clumps from the mapping survey with the Purple Mountain Observatory (PMO) 13.7-m telescope. Follow-up observations suggest that these dense clumps are ideal targets to search for prestellar objects.

The boiling of possible quark nuggets during the quark-hadron phase transition of the Universe at nonzero chemical potential is revisited within the microscopic Brueckner-Hartree-Fock approach employed for the hadron phase, using two kinds of baryon interactions as fundamental inputs. To describe the deconfined phase of quark matter, we use a recently developed quark mass density-dependent model with a fully self-consistent thermodynamic treatment of confinement. We study the baryon number limit $A_{\rm boil}$ (above which boiling may be important) with three typical values for the confinement parameter $D$. It is firstly found that the baryon interaction with a softer equation of state for the hadron phase would only lead to a small increase of $A_{\rm boil}$. However, results depend sensitively on the confinement parameter in the quark model. Specifically, boiling might be important during the Universe cooling for a limited parameter range around $D^{1/2} = 170$ MeV, a value satisfying recent lattice QCD calculations of the vacuum chiral condensate, while for other choices of this parameter, boiling might not happen and cosmological quark nuggets of $10^2 < A < 10^{50}$ could survive.

This work reports a high spatial resolution observations toward Orion KL region with high critical density lines of CH$_{3}$CN (12$_{4}$-11$_{4}$) and CH$_{3}$OH (8$_{-1, 8}$-7$_{0, 7}$) as well as continuum at $\sim$1.3 mm band. The observations were made using the Atacama Large Millimeter/Submillimeter Array with a spatial resolution of $\sim$1.5$^{\prime\prime}$ and sensitives about 0.07 K and $\sim$0.18 K for continuum and line, respectively. The observational results showed that the gas in the Orion KL region consists of jet-propelled cores at the ridge and dense cores at east and south of the region, shaped like a wedge ring. The outflow has multiple lobes, which may originate from an explosive ejection and is not driven by young stellar objects. Four infrared bubbles were found in the Spitzer/IRAC emissions. These bubbles, the distributions of the previously found H$_2$ jets, the young stellar objects and molecular gas suggested that BN is the explosive center. The burst time was estimated to be $\leq$ 1300 years. In the mean time, signatures of gravitational collapse toward Source I and hot core were detected with material infall velocities of 1.5 km~s$^{-1}$ and $\sim$ 0.6 km~s$^{-1}$, corresponding to mass accretion rates of 1.2$\times$10$^{-3}$M$_{\sun}$/Yr and 8.0$\times$10$^{-5}$M$_{\sun}$/Yr, respectively. These observations may support that high-mass stars form via accretion model, like their low-mass counterparts.

(Abridged) We present a spectral analysis of a new, flux-limited sample of 72 X-ray selected clusters of galaxies identified with the X-ray Telescope (XRT) on board the Swift satellite down to a flux limit of ~10-14 erg/s/cm2 (SWXCS, Tundo et al. 2012). We carry out a detailed X-ray spectral analysis with the twofold aim of measuring redshifts and characterizing the properties of the Intra-Cluster Medium (ICM). Optical counterparts and spectroscopic or photometric redshifts are obtained with a cross-correlation with NED. Additional photometric redshifts are computed with a dedicated follow-up program with the TNG and a cross-correlation with the SDSS. We also detect the iron emission lines in 35% of the sample, and hence obtain a robust measure of the X-ray redshift zX. We use zX whenever the optical redshift is not available. Finally, for all the sources with measured redshift, background-subtracted spectra are fitted with a mekal model. We perform extensive spectral simulations to derive an empirical formula to account for fitting bias. The bias-corrected values are then used to investigate the scaling properties of the X-ray observables. Overall, we are able to characterize the ICM of 46 sources. The sample is mostly constituted by clusters with temperatures between 3 and 10 keV, plus 14 low-mass clusters and groups with temperatures below 3 keV. The redshift distribution peaks around z~0.25 and extends up to z~1, with 60% of the sample at 0.1<z<0.4. We derive the Luminosity-Temperature relation for these 46 sources, finding good agreement with previous studies. The quality of the SWXCS sample is comparable to other samples available in the literature and obtained with much larger X-ray telescopes. Our results have interesting implications for the design of future X-ray survey telescopes, characterised by good-quality PSF over the entire field of view and low background.

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

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.

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

Hard X-ray emission in radio-quiet active galactic nuclei (AGNs) is believed to be produced via inverse Compton scattering by hot and compact coronae near the super massive black hole. However the origin and physical properties of the coronae, including geometry, kinematics and dynamics, yet remain poorly known. In this work, taking [OIV] 25.89um emission line as an isotropic indicator of AGN's intrinsic luminosity, we compare the intrinsic corona X-ray emission between Seyfert 1 and Compton-thin Seyfert 2 galaxies, which are viewed at different inclinations according to the unification scheme. We compile a sample of 130 Compton-thin Seyfert galaxies with both [OIV] 25.89um line luminosities measured with Spitzer-IRS and X-ray spectra observed by XMM-Newton, Chandra, Suzaku or Swift. Known radio-loud sources are excluded. We fit the X-ray spectra to obtain the absorption-corrected 2-10 keV continuum luminosities. We find that Seyfert 1 galaxies are intrinsically brighter in intrinsic 2-10 keV emission by a factor of 2.8+0.5-0.4 (2.2+0.9-0.3 in Swift-BAT 14-195 keV emission), comparing with Compton-thin Seyfert 2 galaxies. The Seyfert 1 and Compton-thin Seyfert 2 galaxies follow a statistically identical correlation between the absorption-corrected 2-10 keV luminosity and the 14-195 keV luminosity, indicating that our absorption correction to the 2-10 keV flux is sufficient. The difference in X-ray emission between the two populations is thus unlikely due to X-ray absorption, and instead implies an intrinsic anisotropy in the corona X-ray emission. This striking anisotropy of X-ray emission can be explained by a bipolar outflowing corona with a bulk velocity of ~0.3-0.5c. This would provide a natural link between the so-called coronae and weak jets in these systems. Other consequences of outflowing coronae are also discussed.

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

Gas accretion is a key process in star formation. However, the gas infall detections in high-mass star forming regions with high-spatial resolution observations are rare. Here we report the detection of gas infall towards a cometary ultracompact H{\sc ii} region "C" in G34.26+0.15 complex. The hot core associated with "C" has a mass of $\sim$76 M$_{\sun}$ and a volume density of 1.1$\times10^{8}$ cm$^{-3}$. The HCN (3--2), HCO$^{+}$ (1--0) lines observed by single-dishes and the CN (2--1) lines observed by the SMA show redshifted absorption features, indicating gas infall. We found a linear relationship between the line width and optical depth of the CN (2--1) lines. Those transitions with larger optical depth and line width have larger absorption area. However, the infall velocities measured from different lines seem to be constant, indicating the gas infall is uniform. We also investigated the evolution of gas infall in high-mass star forming regions. At stages prior to hot core phase, the typical infall velocity and mass infall rate are $\sim$ 1 km s$^{-1}$ and $\sim10^{-4}$ M$_{\sun}\cdot$yr$^{-1}$, respectively. While in more evolved regions, the infall velocity and mass infall rates can reach as high as serval km s$^{-1}$ and $\sim10^{-3}-10^{-2}$ M$_{\sun}\cdot$yr$^{-1}$, respectively. Accelerated infall has been detected towards some hypercompact H{\sc ii} and ultracompact H{\sc ii} regions. However, the acceleration phenomenon becomes inapparent in more evolved ultracompact H{\sc ii} regions (e.g. G34.26+0.15).

Planck cold clumps are among the most promising objects to investigate the initial conditions of the evolution of molecular clouds. In this work, by combing the dust emission data from the survey of Planck satellite with the molecular data of $^{12}$CO/$^{13}$CO (1-0) lines from observations with the Purple Mountain Observatory (PMO) 13.7 m telescope, we investigate the CO abundance, CO depletion and CO-to-H$_{2}$ conversion factor of 674 clumps in the early cold cores (ECC) sample. The median and mean values of the CO abundance are 0.89$\times10^{-4}$ and 1.28$\times10^{-4}$, respectively. The mean and median of CO depletion factor are 1.7 and 0.9, respectively. The median value of $X_{CO-to-H_{2}}$ for the whole sample is $2.8\times10^{20}$ cm$^{-2}$K$^{-1}$km$^{-1}$ s. The CO abundance, CO depletion factor and CO-to-H$_{2}$ conversion factor are strongly (anti-)correlated to other physical parameters (e.g. dust temperature, dust emissivity spectral index, column density, volume density and luminosity-to-mass ratio). To conclude, the gaseous CO abundance can be used as an evolutionary tracer for molecular clouds.

Observation of ultra-high energy astrophysical neutrinos and identification of their flavors have been proposed for future neutrino telescopes. The flavor ratio of astrophysical neutrinos observed on the Earth depends on both the initial flavor ratio at the source and flavor transitions taking place during propagations of these neutrinos. The flavor transition mechanisms are well-classified with our model-independent parametrization. We find a new parameter R={\phi}_e/({\phi}_{\mu} + {\phi}_{\tau}) can probe directly the flavor transition in the framework of our model-independent parametrization, without the assumption of the {\nu}_{\mu}-{\nu}_{\tau} symmetry. A few flavor transition models are employed to test our parametrization with this new observable. The observational constraints on flavor transition mechanisms by the new observable is discussed through our model-independent parametrization.

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

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

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

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

Hot matter with nucleons can be produced in the inner region of the neutrino-dominated accretion flow in gamma-ray bursts or during the proto-neutron star birth in successful supernova. The composition and equation of state of the matter depend on the dynamic $\beta$ equilibrium under various neutrino opacities. The strong interaction between nucleons may also play an important role. We plan to extend the previous studies by incorporating these two aspects in our model. The modification of the $\beta$-equilibrium condition from neutrino optically thin to thick has been modeled by an equilibrium factor $\chi$ ranging between the neutrino-freely-escaping case and the neutrino-trapped case. We employ the microscopic Brueckner-Hartree-Fock approach extended to the finite temperature regime to study the interacting nucleons. We show that the composition and chemical potentials of the hot nuclear matter for different densities and temperatures at each stage of $\beta$ equilibrium. We also compare our realistic equation of states with those of the free gas model. We find the neutrino opacity and the strong interaction between nucleons are important for the description and should be taken into account in model calculations.

We present the results of high spatial resolution observations at 1.1 mm waveband, with the Submillimetre Array (SMA), towards the protocluster G10.6-0.4. The 1.1 mm continuum emission reveals seven dense cores, towards which infall motions are all detected with the red-shifted absorption dips in HCN (3--2) line. This is the first time that infall is seen towards multiple sources in a protocluster. We also identified four infrared point sources in this region, which are most likely Class 0/I protostars. Two jet-like structures are also identified from Spitzer/IRAC image. The dense core located in the centre has much larger mass than the off-centre cores. The clump is in overall collapse and the infall motion is supersonic. The standard deviation of core velocities and the velocity differences between the cores and the cloud/clump are all larger than the thermal velocity dispersion. The picture of G10.6-0.4 seems to favor the "competitive accretion" model but needs to be tested by further observations.

The Swift Gamma-ray Burst (GRB) observatory responds to GRB triggers with optical observations in ~ 100 s, but cannot respond faster than ~ 60 s. While some ground-based telescopes respond quickly, the number of sub-60 s detections remains small. In mid- to late-2013, the Ultra-Fast Flash Observatory-Pathfinder is to be launched on the Lomonosov spacecraft to investigate early optical GRB emission. This pathfinder mission is necessarily limited in sensitivity and event rate; here we discuss a next generation rapid-response space observatory. We list science topics motivating our instruments, those that require rapid optical-IR GRB response, including: A survey of GRB rise shapes/times, measurements of optical bulk Lorentz factors, investigation of magnetic dominated (vs. non-magnetic) jet models, internal vs. external shock origin of prompt optical emission, the use of GRBs for cosmology, and dust evaporation in the GRB environment. We also address the impacts of the characteristics of GRB observing on our instrument and observatory design. We describe our instrument designs and choices for a next generation observatory as a second instrument on a low-earth orbit spacecraft, with a 120 kg instrument mass budget. Restricted to relatively modest mass and power, we find that a coded mask X-ray camera with 1024 cm2 of detector area could rapidly locate about 64 GRB triggers/year. Responding to the locations from the X-ray camera, a 30 cm aperture telescope with a beam-steering system for rapid (~ 1 s) response and a near-IR camera should detect ~ 29 GRB, given Swift GRB properties. Am additional optical camera would give a broadband optical-IR slope, allowing dynamic measurement of dust extinction at the source, for the first time.

Aims. We present a stand-alone software (named EXSdetect) for the detection of extended sources in X-ray images. Our goal is to provide a flexible tool capable of detecting extended sources down to the lowest flux levels attainable within instrumental limitations, while maintaining robust photometry, high completeness, and low contamination, regardless of source morphology. EXSdetect was developed mainly to exploit the ever-increasing wealth of archival X-ray data, but is also ideally suited to explore the scientific capabilities of future X-ray facilities, with a strong focus on investigations of distant groups and clusters of galaxies. Methods. EXSdetect combines a fast Voronoi tessellation code with a friends-of-friends algorithm and an automated deblending procedure. The values of key parameters are matched to fundamental telescope properties such as angular resolution and instrumental background. In addition, the software is designed to permit extensive tests of its performance via simulations of a wide range of observational scenarios. Results. We applied EXSdetect to simulated data fields modeled to realistically represent the Swift X-ray Cluster Survey (SXCS), which is based on archival data obtained by the X-ray telescope onboard the Swift satellite. We achieve more than 90% completeness for extended sources comprising at least 80 photons in the 0.5-2 keV band, a limit that corresponds to 10^-14 erg cm^-2 s^-1 for the deepest SXCS fields. This detection limit is comparable to the one attained by the most sensitive cluster surveys conducted with much larger X-ray telescopes. While evaluating the performance of EXSdetect, we also explored the impact of improved angular resolution and discuss the ideal properties of the next generation of X-ray survey missions.

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

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

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

The Planck Early Cold Cores Catalog (ECC) provides an unbiased list of Galactic cold clumps, which form an ideal sample for studying the early phases of star formation (\cite[Planck Collabrators et al. 2011]{Planck_etal11}). To study their properties, we have carried out a molecular line ($^{12}$CO/$^{13}$CO/C$^{18}$O) survey towards 674 Planck cold clumps in the ECC with the PMO 13.7 m telescope.

Soft extended emission (EE) following initial hard spikes up to 100 seconds was observed with {\em Swift}/BAT for about half of short-type gamma-ray bursts (SGRBs). This challenges the conversional central engine models of SGRBs, i.e., compact star merger models. In the framework of the black hole-neutron star merger models, we study the roles of the radial angular momentum transfer in the disk and the magnetic barrier around the black hole for the activity of SGRB central engines. We show that the radial angular momentum transfer may significantly prolong the lifetime of the accretion process and multiple episodes may be switched by the magnetic barrier. Our numerical calculations based on the models of the neutrino-dominated accretion flows suggest that the disk mass is critical for producing the observed EE. In case of the mass being $\sim 0.8M_{\odot}$, our model can reproduce the observed timescale and luminosity of both the main and EE episodes in a reasonable parameter set. The predicted luminosity of the EE component is lower than the observed EE with about one order of magnitude and the timescale is shorter than 20 seconds if the disk mass being $\sim 0.2M_{\odot}$. {\em Swift}/BAT-like instruments may be not sensitive enough to detect the EE component in this case. We argue that the EE component would be a probe for merger process and disk formation for compact star mergers.

The environments surrounding nine Wolf-Rayet stars were studied in molecular emission. Expanding shells were detected surrounding these WR stars (see left panels of Figure 1). The average masses and radii of the molecular cores surrounding these WR stars anti-correlate with the WR stellar wind velocities (middle panels of Figure 1), indicating the WR stars has great impact on their environments. The number density of Young Stellar Objects (YSOs) is enhanced in the molecular shells at $\sim$5 arcmin from the central WR star (lower-right panel of Figure 1). Through detailed studies of the molecular shells and YSOs, we find strong evidences of triggered star formation in the fragmented molecular shells (\cite[Liu et al. 2010]{liu_etal12}

A mapping survey towards 51 Planck cold clumps projected on Orion complex was performed with J=1-0 lines of $^{12}$CO and $^{13}$CO at the 13.7 m telescope of Purple Mountain Observatory. The mean column densities of the Planck gas clumps range from 0.5 to 9.5$\times10^{21}$ cm$^{-2}$, with an average value of (2.9$\pm$1.9)$\times10^{21}$ cm$^{-2}$. While the mean excitation temperatures of these clumps range from 7.4 to 21.1 K, with an average value of 12.1$\pm$3.0 K. The averaged three-dimensional velocity dispersion $\sigma_{3D}$ in these molecular clumps is 0.66$\pm$0.24 km s$^{-1}$. Most of the clumps have $\sigma_{NT}$ larger than or comparable with $\sigma_{Therm}$. The H$_{2}$ column density of the molecular clumps calculated from molecular lines correlates with the aperture flux at 857 GHz of the dust emission. Through analyzing the distributions of the physical parameters, we suggest turbulent flows can shape the clump structure and dominate their density distribution in large scale, but not affect in small scale due to the local fluctuations. Eighty two dense cores are identified in the molecular clumps. The dense cores have an averaged radius and LTE mass of 0.34$\pm$0.14 pc and 38$_{-30}^{+5}$ M$_{\sun}$, respectively. And structures of low column density cores are more affected by turbulence, while those of high column density cores are more concerned by other factors, especially by gravity. The correlation of the velocity dispersion versus core size is very weak for the dense cores. The dense cores are found most likely gravitationally bounded rather than pressure confined. The relationship between $M_{vir}$ and $M_{LTE}$ can be well fitted with a power law. The core mass function here is much more flatten than the stellar initial mass function. The lognormal behavior of the core mass distribution is most likely determined by the internal turbulence.

A survey toward 674 Planck cold clumps of the Early Cold Core Catalogue (ECC) in the J=1-0 transitions of $^{12}$CO, $^{13}$CO and C$^{18}$O has been carried out using the PMO 13.7 m telescope. 673 clumps were detected with the $^{12}$CO and $^{13}$CO, and 68% of the samples have C$^{18}$O emission. Additional velocity components were also identified.A close consistency of the three line peak velocities was revealed for the first time. Kinematic distances are given out for all the velocity components and half of the clumps are located within 0.5 and 1.5 kpc. Excitation temperatures range from 4 to 27 K, slightly larger than those of $T_d$. Line width analysis shows that the majority of ECC clumps are low mass clumps. Column densities N$_{H_{2}}$ span from 10$^{20}$ to 4.5$\times10^{22}$ cm$^{-2}$ with an average value of (4.4$\pm$3.6)$\times10^{21}$ cm$^{-2}$. N$_{H_{2}}$ cumulative fraction distribution deviates from the lognormal distribution, which is attributed to optical depth. The average abundance ratio of the $^{13}$CO to C$^{18}$O in these clumps is 7.0$\pm$3.8, higher than the terrestrial value. Dust and gas are well coupled in 95% of the clumps. Blue profile, red profile and line asymmetry in total was found in less than 10% of the clumps, generally indicating star formation is not developed yet. Ten clumps were mapped. Twelve velocity components and 22 cores were obtained. Their morphologies include extended diffuse, dense isolated, cometary and filament, of which the last is the majority. 20 cores are starless.Only 7 cores seem to be in gravitationally bound state. Planck cold clumps are the most quiescent among the samples of weak-red IRAS, infrared dark clouds, UC H{\sc ii} region candidates, EGOs and methanol maser sources, suggesting that Planck cold clumps have expanded the horizon of cold Astronomy.

The environment surrounding Wolf-Rayet star HD 211853 is studied in molecular emission, infrared emission, as well as radio and HI emission. The molecular ring consists of well-separated cores, which have a volume density of 10$^{3}$ cm$^{-3}$ and kinematic temperature $\sim$20 K. Most of the cores are under gravitational collapse due to external pressure from the surrounding ionized gas. From SED modeling towards the young stellar objects (YSOs), sequential star formation is revealed on a large scale in space spreading from the Wolf-Rayet star to the molecular ring. A small scale sequential star formation is revealed towards core A, which harbors a very young star cluster. Triggered star formations is thus suggested. The presence of PDR, the fragmentation of the molecular ring, the collapse of the cores, the large scale sequential star formation indicate the "Collect and Collapse" process functions in this region. The star forming activities in core A seem to be affected by the "Radiation-Driven Implosion" (RDI) process.

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

We present the result of the Chandra high-resolution observation of the Seyfert~2 galaxy NGC 7590. This object was reported to show no X-ray absorption in the low-spatial resolution ASCA data. The XMM observations show that the X-ray emission of NGC 7590 is dominated by an off-nuclear ultra-luminous X-ray source (ULX) and an extended emission from the host galaxy, and the nucleus is rather weak, likely hosting a Compton-thick AGN. Our recent Chandra observation of NGC 7590 enables to remove the X-ray contamination from the ULX and the extended component effectively. The nuclear source remains undetected at ~4x10^{-15} erg/s/cm^-2 flux level. Although not detected, Chandra data gives a 2--10 keV flux upper limit of ~6.1x10^{-15} erg/s/cm^-2 (at 3 sigma level), a factor of 3 less than the XMM value, strongly supporting the Compton-thick nature of the nucleus. In addition, we detected five off-nuclear X-ray point sources within the galaxy D25 ellipse, all with 2 -- 10 keV luminosity above 2x10^{38} erg/s (assuming the distance of NGC 7590). Particularly, the ULX previously identified by ROSAT data was resolved by Chandra into two distinct X-ray sources. Our analysis highlights the importance of high spatial resolution images in discovering and studying ULXs.

Using the 13.7 m telescope of Purple Mountain Observatory (PMO), a survey of J=1-0 lines of CO and its isotopes was carried out towards 98 methanol maser sources in January 2008. Eighty-five sources have infrared counterparts within one arcmin. In the survey, except 43 sources showing complex or multiple-peak profiles, almost all the $^{13}$CO line profiles of the other 55 sources have large line widths of 4.5 km s$^{-1}$ on average and are usually asymmetric. Fifty corresponding Infrared Astronomical Satellite (IRAS) sources of these 55 sources are with $L_{bol}$ larger than $10^{3}L_{\odot}$, which can be identified as possible high-mass young stellar sources. Statistics show that the $^{13}$CO line widths correlate with the bolometric luminosity of the associated IRAS sources. We also report the mapping results of two sources: IRAS 06117+1350 and IRAS 07299-1651 here. Two cores were found in IRAS 06117+1350 and one core was detected in IRAS 07299-1651. The northwest core in IRAS 06117+1350 and the core in IRAS 07299-1651 can be identified as precursors of UC\simH{\sc ii} regions or high-mass protostellar objects (HMPOs). The southeast core of IRAS 06117+1350 has no infrared counterpart, seeming to be on younger stages than pre-UC\simH{\sc ii} phase.

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.

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

We address the possibility that CP-violation in $B_s-\bar B_s$ mixing may help explain the origin of the cosmic baryon asymmetry. We propose a new baryogenesis mechanism - "Electroweak Beautygenesis" - explicitly showing that these two CP-violating phenomena can be sourced by a common CP-phase. As an illustration, we work in the Two-Higgs-Doublet model. Because the relevant CP-phase is flavor off-diagonal, this mechanism is less severely constrained by null results of electric dipole moment searches than other scenarios. We show how measurements of flavor observables by the D0, CDF, and LHCb collaborations test this scenario.

We study a limit of the nearly-Peccei-Quinn-symmetric Next-to-Minimal Supersymmetric Standard Model possessing novel Higgs and dark matter (DM) properties. In this scenario, there naturally co-exist three light singlet-like particles: a scalar, a pseudoscalar, and a singlino-like DM candidate, all with masses of order 0.1-10 GeV. The decay of a Standard Model-like Higgs boson to pairs of the light scalars or pseudoscalars is generically suppressed, avoiding constraints from collider searches for these channels. For a certain parameter window annihilation into the light pseudoscalar and exchange of the light scalar with nucleons allow the singlino to achieve the correct relic density and a large direct detection cross section consistent with the CoGeNT and DAMA/LIBRA preferred region simultaneously. This parameter space is consistent with experimental constraints from LEP, the Tevatron, and Upsilon- and flavor physics.

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

We report on studies of the viability and sensitivity of the Askaryan Radio Array (ARA), a new initiative to develop a Teraton-scale ultra-high energy neutrino detector in deep, radio-transparent ice near Amundsen-Scott station at the South Pole. An initial prototype ARA detector system was installed in January 2011, and has been operating continuously since then. We report on studies of the background radio noise levels, the radio clarity of the ice, and the estimated sensitivity of the planned ARA array given these results, based on the first five months of operation. Anthropogenic radio interference in the vicinity of the South Pole currently leads to a few-percent loss of data, but no overall effect on the background noise levels, which are dominated by the thermal noise floor of the cold polar ice, and galactic noise at lower frequencies. We have also successfully detected signals originating from a 2.5 km deep impulse generator at a distance of over 3 km from our prototype detector, confirming prior estimates of kilometer-scale attenuation lengths for cold polar ice. These are also the first such measurements for propagation over such large slant distances in ice. Based on these data, ARA-37, the 200 km^2 array now under construction, will achieve the highest sensitivity of any planned or existing neutrino detector in the 10^{16}-10^{19} eV energy range.

We present an observational study towards the young high-mass star forming region G23.44-0.18 using the Submillimeter Array. Two massive, radio-quiet dusty cores MM1 and MM2 are observed in 1.3 mm continuum emission and dense molecular gas tracers including thermal CH$_3$OH, CH$_3$CN, HNCO, SO, and OCS lines. The $^{12}$CO (2--1) line reveals a strong bipolar outflow originated from MM2. The outflow consists of a low-velocity component with wide-angle quasi-parabolic shape and a more compact and collimated high-velocity component. The overall geometry resembles the outflow system observed in the low-mass protostar which has a jet-driven fast flow and entrained gas shell. The outflow has a dynamical age of $6\times10^3$ years and a mass ejection rate $\sim10^{-3} M_{\odot}$ year$^{-1}$. A prominent shock emission in the outflow is observed in SO and OCS, and also detected in CH$_3$OH and HNCO. We investigated the chemistry of MM1, MM2 and the shocked region. The dense core MM2 have molecular abundances of 3 to 4 times higher than those in MM1. The abundance excess, we suggest, can be a net effect of the stellar evolution and embedded shocks in MM2 that calls for further inspection.

The Ultra-Fast Flash Observatory (UFFO) pathfinder is a space mission devoted to the measurement of Gamma-Ray Bursts (GRBs), especially their early light curves which will give crucial information on the progenitor stars and central engines of the GRBs. It consists of two instruments: the UFFO Burst Alert & Trigger telescope (UBAT) for the detection of GRB locations and the Slewing Mirror Telescope (SMT) for the UV/optical afterglow observations, upon triggering by UBAT. The UBAT employs a coded-mask {\gamma}/X-ray camera with a wide field of view (FOV), and is comprised of three parts: a coded mask, a hopper, and a detector module (DM). The UBAT DM consists of a LYSO scintillator crystal array, multi-anode photo multipliers, and analog and digital readout electronics. We present here the design and fabrication of the UBAT DM, as well as its preliminary test results.

Hundreds of gamma-ray burst (GRB) optical light curves have been measured since the discovery of optical afterglows. However, even after nearly 7 years of operation of the Swift Observatory, only a handful of measurements have been made soon (within a minute) after the gamma ray signal. This lack of early observations fails to address burst physics at short time scales associated with prompt emissions and progenitors. Because of this lack of sub-minute data, the characteristics of the rise phase of optical light curve of short-hard type GRB and rapid-rising GRB, which may account for ~30% of all GRB, remain practically unknown. We have developed methods for reaching sub-minute and sub-second timescales in a small spacecraft observatory. Rather than slewing the entire spacecraft to aim the optical instrument at the GRB position, we use rapidly moving mirror to redirect our optical beam. As a first step, we employ motorized slewing mirror telescope (SMT), which can point to the event within 1s, in the UFFO Pathfinder GRB Telescope onboard the Lomonosov satellite to be launched in Nov. 2011. UFFO's sub-minute measurements of the optical emission of dozens of GRB each year will result in a more rigorous test of current internal shock models, probe the extremes of bulk Lorentz factors, provide the first early and detailed measurements of fast-rise GRB optical light curves, and help verify the prospect of GRB as a new standard candle. We will describe the science and the mission of the current UFFO Pathfinder project, and our plan of a full-scale UFFO-100 as the next step.

With the KOSMA 3-m telescope, 54 Herbig Ae/Be stars were surveyed in CO and $^{13}$CO emission lines. The properties of the stars and their circumstellar environments are studied by fitting the SEDs. The mean line width of $^{13}$CO (2-1) lines of this sample is 1.87 km s$^{-1}$. The average column density of H$_{2}$ is found to be $4.9\times10^{21}$ cm$^{-2}$ for the stars younger than $10^{6}$ yr, while drops to $2.5\times10^{21}$ cm$^{-2}$ for those older than $10^{6}$ yr. No significant difference is found among the SEDs of Herbig Ae stars and Herbig Be stars at the same age. The infrared excess decreases with age. The envelope masses and the envelope accretion rates decease with age after $10^{5}$ yr. The average disk mass of the sample is $3.3\times10^{-2} M_{\sun}$. The disk accretion rate decreases more slowly than the envelope accretion rate. A strong correlation between the CO line intensity and the envelope mass is found.

[abridged] We present point-source catalogs for the 4Ms Chandra Deep Field-South (CDF-S), which is the deepest Chandra survey to date and covers an area of 464.5 arcmin^2. We provide a main source catalog, which contains 740 X-ray point sources that are detected with wavdetect at a false-positive probability threshold of 1E-5 and also satisfy a binomial-probability source-selection criterion of P<0.004; this approach is designed to maximize the number of reliable sources detected. A total of 300 main-catalog sources are new compared to the previous 2Ms CDF-S main-catalog sources. We also provide a supplementary catalog, which consists of 36 sources that are detected with wavdetect at 1E-5, satisfy 0.004< P<0.1, and have an optical counterpart with R<24. Multiwavelength identifications, basic optical/infrared/radio photometry, and spectroscopic/photometric redshifts are provided for the X-ray sources. Basic analyses of the X-ray and multiwavelength properties of the sources indicate that >75% of the main-catalog sources are AGNs; of the 300 new main-catalog sources, about 35% are likely normal and starburst galaxies, reflecting the rise of normal and starburst galaxies at the very faint flux levels uniquely accessible to the 4Ms CDF-S. Near the center of the 4Ms CDF-S, the observed AGN and galaxy source densities have reached ~9800 and 6900 per square degree, respectively. The 4 Ms CDF-S reaches on-axis flux limits of ~9.1E-18 and 5.5E-17 erg/cm^2/s for the soft and hard bands, respectively. An increase in the CDF-S exposure by a factor of ~2-2.5 would provide further significant gains and probe key unexplored discovery space.

We argue that an effective flavor discrimination in neutrino telescopes is the key to probe the flavor ratios of astrophysical neutrinos at the source and flavor transition mechanisms of these neutrinos during their propagations from the source to the Earth. We first discuss how well one can reconstruct the flavor ratios of astrophysical neutrinos at the source. We then discuss how to probe flavor transition mechanisms of propagating astrophysical neutrinos. In this regard, we propose a model independent parametrization for neutrino flavor transitions, with the neutrino oscillation as a special case.

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.

One of the important goals for future neutrino telescopes is to identify the flavors of astrophysical neutrinos and therefore determine the flavor ratio. The flavor ratio of astrophysical neutrinos observed on the Earth depends on both the initial flavor ratio at the source and flavor transitions taking place during propagations of these neutrinos. We propose a model independent parametrization for describing the above flavor transitions. A few flavor transition models are employed to test our parametrization. The observational test for flavor transition mechanisms through our parametrization is discussed.

We compile a sample of 89 Seyfert galaxies with both [OIV] 25.89 um line luminosities observed by Spitzer IRS and X-ray spectra observed by XMM EPIC. Using [OIV] emission as proxy of AGN intrinsic luminosity, we find that although type 2 AGNs have higher line equivalent width, the narrow Fe Kalpha line in Compton-Thin and Compton-Thick Seyfert 2 galaxies are 2.9 +0.8 -0.6 and 5.6 +1.9 -1.4 times weaker in terms of luminosity than Seyfert 1 galaxies respectively. This indicates different correction factors need to be applied for various types of AGNs before the narrow Fe Kalpha line luminosity could serve as intrinsic AGN luminosity indicator. We also find Seyfert 1 galaxies in our sample have on average marginally larger line width and higher line centroid energy, suggesting contamination from highly ionized Fe line or broader line emission from much smaller radius, but this effect is too weak to explain the large difference in narrow Fe Kalpha line luminosity between type 1 and type 2 AGNs. This is the first observational evidence showing the narrow Fe Kalpha line emission in AGNs is anisotropic. The observed difference is consistent with theoretical calculations assuming a smoothly distributed obscuring torus, and could provide independent constraints on the clumpiness of the torus.

We present the analysis of three XMM observations of the Seyfert 2 galaxy NGC 7590. The source was found to have no X-ray absorption in the low spatial resolution ASCA data. The XMM observations provide a factor of 10 better spatial resolution than previous ASCA data. We find that the X-ray emission of NGC 7590 is dominated by an off-nuclear ULX and extended emission from the host galaxy. The nuclear X-ray emission is rather weak comparing with the host galaxy. Based on its very low X-ray luminosity as well as the small ratio between the 2-10 keV and the [O III] fluxes, we interpret NGC 7590 as Compton-thick rather than being an "unobscured" Seyfert 2 galaxy. Future higher resolution observations such as Chandra are crucial to shed light on the nature of NGC 7590 nucleus.

Cosmic ray anomalies observed by PAMELA and Fermi-LAT experiments may be interpreted by heavy (TeV-scale) dark matter annihilation enhanced by Sommerfeld effects mediated by a very light (sub-GeV) U(1)_X gauge boson, while the recent direct searches from CoGeNT and DAMA/LIBRA experiments may indicate a rather light (\sim 7 GeV) dark matter with weak interaction. Motivated by these apparently different scales, we consider a gauge mediated next-to-the minimal supersymmetric standard model (NMSSM) entended with a light U(1)_X sector plus a heavy sector (\bar H_h,H_h), which can provide both a light (\sim 7 GeV) and a heavy (TeV-scale) dark matter without introducing any ad hoc new scale. Through the Yukawa coupling between H_h and the messager fields, the U(1)_X gauge symmetry is broken around the GeV scale radiatively and a large negative m_S^2 is generated for the NMSSM singlet S. Furthermore, the small kinetic mixing parameter between U(1)_X and U(1)_Y is predicted to be \theta\sim 10^{-5}-10^{-6} after integrating out the messengers. Such a light dark matter, which can have a normal relic density from the late decay of the right-handed sneutrino (assumed to be the ordinary next-to-the lightest supersymmetric particle and thermally produced in the early universe), can serve a good candidate to explain the recent CoGeNT and DAMA/LIBRA results.

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

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

We discuss the reconstruction of neutrino flavor ratios at astrophysical sources through the future neutrino-telescope measurements. Taking the ranges of neutrino mixing parameters $\theta_{ij}$ as those given by the current global fit, we demonstrate by a statistical method that the accuracies in the measurements of energy-independent ratios $R\equiv\phi (\nu_{\mu})/(\phi (\nu_{e})+\phi (\nu_{\tau}))$ and $S\equiv\phi (\nu_e)/\phi (\nu_{\tau})$ among integrated neutrino flux should both be better than 10% in order to distinguish between the pion source and the muon-damped source at the $3 \sigma$ level. The 10% accuracy needed for measuring $R$ and $S$ requires an improved understanding on the background atmospheric neutrino flux to a better than 10% level in the future. We discuss the applicability of our analysis to practical situations that the diffuse astrophysical neutrino flux arises from different types of sources and each point source has a neutrino flavor ratio varying with energies. We also discuss the effect of leptonic CP phase on the flavor-ratio reconstruction.

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