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.

Dense molecular filaments are central to the star formation process, but the detailed manner in which they fragment into prestellar cores is not yet well understood. Here, we investigate the fragmentation properties and dynamical state of several star-forming filaments in the X-shape Nebula region of the California MC, in an effort to shed some light on this issue. We used multi-wavelength far-infrared images from {\it Herschel} and the \textsl{getsources} and \textsl{getfilaments} extraction methods to identify dense cores and filaments and derive their basic properties. We also used a map of $\rm ^{13}CO (2-1)$ emission from SMT 10m submillimeter telescope to constrain the dynamical state of the filaments. We identified 10 filaments, as well as 57 dense cores. Two star-forming filaments (\# 8 and \# 10) stand out in that they harbor quasi-periodic chains of dense cores with a typical projected core spacing of $\sim$0.15 pc. These two filaments have thermally supercritical line masses and are not static. Filament~8 exhibits a prominent transverse velocity gradient, suggesting that it is accreting gas from the parent cloud gas reservoir. In both cases, the observed (projected) core spacing is similar to the filament width and significantly shorter than the canonical separation of $\sim \,$4 times the filament width predicted by classical cylinder fragmentation theory. We suggest that continuous accretion of gas onto the two star-forming filaments, as well as geometrical bending of the filaments, may account for the observed core spacing. Our findings suggest that the characteristic fragmentation lengthscale of molecular filaments is quite sensitive to external perturbations from the parent cloud, such as gravitational accretion of ambient material.

It is generally believed that the variability of photospheric emission in gamma-ray bursts (GRBs) traces that of the jet power. This work further investigates the variability of photospheric emission in a variable jet. By setting a constant $\eta$ (dimensionless entropy of the jet), we find that the light curve of the photospheric emission shows a ``tracking'' pattern on the time profile of jet power. However, the relative variability is significantly low in the photospheric emission compared with that in the jet power. If the $\eta$ is genetic variable, the variability of the photospheric emission is not only limited by the jet power but also affected by $\eta$ strongly. It becomes complex and is generally different from that of the jet power. Moreover, the opposite phase may stand in the variabilities of the photospheric emission at different photon energies. We also find that the relative variability does not remain constant over the photon energies with an obvious reduction at a certain energy. This is consistent with the analysis of GRB 090902B in which an appreciable thermal component has been detected in a wide energy range. For several other GRBs coupling with the thermal component, we conservatively evaluate the variability of the thermal and non-thermal emission, respectively. Our results show that the relative variability of the thermal emission is likely comparable to that of the non-thermal emission for these bursts. In addition, the analysis of GRB~120323A reveals that the variability of the photospheric emission may be of the opposite phase from that of the non-thermal emission.

Dark energy (DE) plays an important role in the expansion history of our universe. In the earlier works, universe including DE with a constant or time-varying equation of state (EOS) has been well studied. Different from predecessors, we investigate the cosmic expansion with multicomponent DE. Today's total energy density of DE is separated into several parts equally and every part has a constant EOS $w_i$. We modify the Friedmann equation and the parameterized post-Friedmann description of dark energy, then put constraints on $w_i$ from Planck 2018 TT,TE,EE$+$lowE$+$lensing, BAO data, PANTHEON samples and SH0ES measurement. We find that the largest and the smallest values of $w_i$ have no overlap at $95\%$C.L.. Our results indicate that DE is composed of more than one candidate. Besides, the density evolution of multicomponent DE is tick-like and different from the other DE models'. That is to say, multicomponent DE makes a contribution to late and early universe.

Cosmic Microwave Background (CMB) anisotropy encodes a lot of information about our Universe. In this paper we take the ground-based CMB observations (GCMB), including the South Pole Telescope (SPT), SPTpol and the Atacama Cosmology Telescope Polarimeter (ACTPol), as a new probe to the CMB anisotropy independent of two satellite observations, i.e. Wilkinson Microwave Anisotropy Probe (WMAP) and Planck. The combination of current GCMB data is consistent with WMAP and Planck. In the spatially flat $\Lambda$CDM model, the Hubble constant is $H_0=69.72\pm 1.63$ km/s/Mpc at $68\%$ confidence level (CL). Combining with baryon acoustic oscillation (BAO) and the Pantheon sample of Type Ia supernovae (SN), we find that $H_0=68.40\pm 0.58$ km/s/Mpc ($68\%$ CL) in the spatially flat $\Lambda$CDM cosmology which has a tension with local measurement given by Riess et al. in 2019 at $3.7\sigma$ level, and $\Omega_k=-0.0013\pm 0.0039$ and $N_{\rm{eff}}=2.90\pm 0.41$ ($68\%$ CL) in the extended cosmological models.

We simulate the evolution of a dust universe from $z=1089$ to $z=0$ by numerically integrating the Einstein's equation for a spatially flat Friedmann-Lemaire-Robertson-Walker (FLRW) background spacetime with scalar perturbations which are derived from the matter power spectrum produced with the Code for Anisotropies in the Microwave Background (CAMB). To investigate the effects of primordial gravitational waves (GWs) on the inhomogeneity of the universe, we add an additional decaying, divergenceless and traceless primordial tensor perturbation with its initial amplitude being $3\times 10^{-4}$ to the above metric. We find that this primordial tensor perturbation suppresses the matter power spectrum by about $0.01\%$ at $z=0$ for modes with wave number similar to its. This suppression may be a possible probe of a GWs background in the future.

Identifying galaxy groups from redshift surveys of galaxies plays an important role in connecting galaxies with the underlying dark matter distribution. Current and future high-$z$ spectroscopic surveys, usually incomplete in redshift sampling, present both opportunities and challenges to identifying groups in the high-$z$ Universe. We develop a group finder that is based on incomplete redshift samples combined with photometric data, using a machine learning method to assign halo masses to identified groups. Test using realistic mock catalogs shows that $\gtrsim 90\%$ of true groups with halo masses $\rm M_h \gtrsim 10^{12} M_{\odot}/h$ are successfully identified, and that the fraction of contaminants is smaller than $10\%$. The standard deviation in the halo mass estimation is smaller than 0.25 dex at all masses. We apply our group finder to zCOSMOS-bright and describe basic properties of the group catalog obtained.

Newton's constant has observational effects on both the CMB power spectra and the light curves of SNIa. We use Planck data, BAO data and the SNIa measurement to constrain the varying Newton's constant $G$ during the CMB epoch and the redshift ranges of PANTHEON samples, and find no evidence indicating that $G$ is varying with redshift. By extending the $\Lambda$CDM model with one free parameter $G$, we get $G =(6.65635_{-0.18560}^{+0.18766} ) \times 10^{-11} \rm m^3kg^{-1}s^{-2}$ and $H_0=67.62^{+1.24}_{-1.25} $ km s$^{-1}$ Mpc$^{-1}$ at 68$\%$ C.L. from Planck$+$BAO$+$uncalibrated PANTHEON. The results show the value of $G$ is consistent with CODATA 2018, but the $H_0$ tension can't be solved in this way.

We present 1.3 mm ALMA dust polarization observations at a resolution of $\sim$0.02 pc of three massive molecular clumps, MM1, MM4, and MM9, in the infrared dark cloud G28.34+0.06. With the sensitive and high-resolution continuum data, MM1 is resolved into a cluster of condensations. The magnetic field structure in each clump is revealed by the polarized emission. We found a trend of decreasing polarized emission fraction with increasing Stokes $I$ intensities in MM1 and MM4. Using the angular dispersion function method (a modified Davis-Chandrasekhar-Fermi method), the plane-of-sky magnetic field strength in two massive dense cores, MM1-Core1 and MM4-Core4, are estimated to be $\sim$1.6 mG and $\sim$0.32 mG, respectively. \textbf{The ordered magnetic energy is found to be smaller than the turbulent energy in the two cores, while the total magnetic energy is found to be comparable to the turbulent energy.} The total virial parameters in MM1-Core1 and MM4-Core4 are calculated to be $\sim$0.76 and $\sim$0.37, respectively, suggesting that massive star formation does not start in equilibrium. Using the polarization-intensity gradient-local gravity method, we found that the local gravity is closely aligned with intensity gradient in the three clumps, and the magnetic field tends to be aligned with the local gravity in MM1 and MM4 except for regions near the emission peak, which suggests that the gravity plays a dominant role in regulating the gas collapse. Half of the outflows in MM4 and MM9 are found to be aligned within 10$^{\circ}$ of the condensation-scale ($<$0.05 pc) magnetic field, indicating that the magnetic field could play an important role from condensation to disk scale in the early stage of massive star formation. We also found that the fragmentation in MM1-Core1 cannot be solely explained by thermal Jeans fragmentation or turbulent Jeans fragmentation.

The concentration parameter is a key characteristic of a dark matter halo that conveniently connects the halo's present-day structure with its assembly history. Using 'Dark Sky', a suite of cosmological $N$-body simulations, we investigate how halo concentration evolves with time and emerges from the mass assembly history. We also explore the origin of the scatter in the relation between concentration and assembly history. We show that the evolution of halo concentration has two primary modes: (1) smooth increase due to pseudo-evolution; and (2) intense responses to physical merger events. Merger events induce lasting and substantial changes in halo structures, and we observe a universal response in the concentration parameter. We argue that merger events are a major contributor to the uncertainty in halo concentration at fixed halo mass and formation time. In fact, even haloes that are typically classified as having quiescent formation histories experience multiple minor mergers. These minor mergers drive small deviations from pseudo-evolution, which cause fluctuations in the concentration parameters and result in effectively irreducible scatter in the relation between concentration and assembly history. Hence, caution should be taken when using present-day halo concentration parameter as a proxy for the halo assembly history, especially if the recent merger history is unknown.

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.

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.

Although the physical origin of prompt emission in gamma-ray bursts (GRBs) remains inconclusive, previous studies have considered the synchrotron radiation of relativistic electrons as a promising mechanism. These works usually adopted a invariable injection rate of electrons ($Q$) which may be discordant with that in a Poynting-flux dominated jet. In a Poynting-flux dominated jet (e.g., ICMART model, Zhang & Yan 2011), the number of magnetic reconnections occurred simultaneously may grow rapidly with time and results in an increase of $Q$ with time. This paper is dedicated to study the synchrotron radiation spectrum in this scenario. It is found that the radiation spectrum would obviously get harder if an increasing $Q$ is adopted and a Band-like radiation spectrum can be obtained if the increase of $Q$ is fast enough. The latter is related to the fact that a bump-shape rather than a power-law spectrum appears in the low-energy regime of the obtained electron spectrum. This effect can strongly harden the low-energy radiation spectrum. It indicates that an increasing $Q$ can help to alleviate the "fast-cooling problem" of synchrotron radiation for GRBs. Our studies also reveal that a Poynting-flux dominated jet with a large emission radius, a small length of the magnetic reconnection region, or a low-minimum energy of injected electron would prefer to form a Band-like radiation spectrum. We suggest that the Band spectrum found in GRBs may be the synchrotron emission of the electrons with a bump-shape distribution in its low-energy regime.

Solar system materials are variably depleted in moderately volatile elements (MVEs) relative to the proto-solar composition. To address the origin of this MVE depletion, we conducted a systematic study of high-precision K isotopic composition on 16 carbonaceous chondrites (CCs) of types CM1-2, CO3, CV3, CR2, CK4-5 and CH3 and 28 ordinary chondrites (OCs) covering petrological types 3 to 6 and chemical groups H, L, and LL. We observed significant overall K isotope (delta41K) variations (-1.54 to 0.70 permil). The K isotope compositions of CCs are largely higher than the Bulk Silicate Earth (BSE) value, whereas OCs show typically lower values than BSE. Neither CCs nor OCs show resolvable correlations between K isotopes and chemical groups, petrological types, shock levels, exposure ages, fall or find occurrence, or terrestrial weathering. The lack of a clear trend between K isotopes and K content indicates that the K isotope fractionations were decoupled from the relative elemental K depletions. The range of K isotope variations in the CCs is consistent with a four-component (chondrule, refractory inclusion, matrix and water) mixing model that is able to explain the bulk elemental and isotopic compositions of the main CC groups, but requires a fractionation in K isotopic compositions in chondrules. We propose that the major control of the isotopic compositions of group averages is condensation or vaporization in nebular environments that is preserved in the compositional variation of chondrules. Parent-body processes (aqueous alteration, thermal metamorphism, and metasomatism) can mobilize K and affect the K isotopes in individual samples. In the case of the OCs, the full range of K isotopic variations can only be explained by the combined effects of the size and relative abundances of chondrules, parent-body aqueous and thermal alteration.

Enstatite chondrites and aubrites are meteorites that show the closest similarities to the Earth in many isotope systems that undergo mass-independent and mass-dependent isotope fractionations. Due to the analytical challenges to obtain high-precision K isotopic compositions in the past, potential differences in K isotopic compositions between enstatite meteorites and the Earth remained uncertain. We report the first high-precision K isotopic compositions of eight enstatite chondrites and four aubrites and find that there is a significant variation of K isotopic compositions among enstatite meteorites (from -2.34 permil to -0.18 permil). However, K isotopic compositions of nearly all enstatite meteorites scatter around the Bulk Silicate Earth (BSE) value. The average K isotopic composition of the eight enstatite chondrites (-0.47 +/- 0.57 permil) is indistinguishable from the BSE value (-0.48 +/- 0.03 permil), thus further corroborating the isotopic similarity between Earth' building blocks and enstatite meteorite precursors. We found no correlation of K isotopic compositions with the chemical groups, petrological types, shock degrees, and terrestrial weathering conditions; however, the variation of K isotopes among enstatite meteorite can be attributed to the parent body processing. Our sample of the main group aubrite MIL 13004 is exceptional and has an extremely light K isotopic composition (delta 41K= -2.34 +/- 0.12 permil). We attribute this unique K isotopic feature to the presence of abundant djerfisherite inclusions in our sample because this K-bearing sulfide mineral is predicted to be enriched in 39K during equilibrium exchange with silicates.

We use the ELUCID N-Body simulation to study the relation of structural properties of dark matter halos to their assembly history and environment. We find that the complexity of halo assembly histories can be well described by a small number of principal components. Using decision trees built with the random ensemble method, we find that the correlation between halo concentration and assembly history is tight: more than 60% of the variance in halo concentration can be explained by assembly history alone. Halo shape and spin parameters are found to depend not only on assembly history, but also on environment and initial condition, and these dependencies are not orthogonal. The structures of more massive halos depend more strongly on higher-order components of the assembly history. The relations between halo intrinsic properties and environment are weak compared to their variances, with the anisotropy of local tidal fields having the strongest correlation with halo properties. Our method of dimension reduction and regression can help simplify the characterization of the halo population and clarify the degeneracy among halo properties.

Using the multiconfiguration Dirac-Hartree-Fock and the relativistic configuration interaction methods, a consistent set of transition energies and radiative transition data for the lowest 546 (623, 701, 745) states of the $3p^4 3d$, $3s 3p^2 3d^2$, $3s 3p^3 4p$, $3s 3p^4$, $3s^2 3d^3$, $3s^2 3p^2 3d$, $3s^2 3p^2 4d$, $3s^2 3p^2 4s$, $3p^3 3d^2$, $3p^5$, $3s 3p 3d^3$, $3s 3p^3 3d$, $3s 3p^3 4s$, $3s^2 3p 3d^2$, %$3s^2 3p^2 4f$, $3s^2 3p^2 4p$, $3s^2 3p^3$ configurations in Mn~XI (Fe~XII, Co~XIII, Ni~XIV) is provided. The comparison between calculated excitation energies for the $n=4$ states and available experimental values for Fe XII indicate that the calculations are highly accurate, with uncertainties of only a few hundred cm$^{-1}$. Lines from these states are prominent in the soft X-rays. With the present calculations, several recent new identifications are confirmed. Other identifications involving $3p^2 4d$ levels in Fe~XII that were found questionable are discussed and a few new assignments are recommended. As some $n=4$ states of the other ions also show large discrepancies between experimental and calculated energies, we reassess their identification. The present study provides highly accurate atomic data for the $n=4$ states of P-like ions of astrophysical interest, for which experimental data are scarce.

It has long been speculated that many starburst or compact dwarf galaxies are resulted from dwarf-dwarf galaxy merging, but unequivocal evidence for this possibility has rarely been reported in the literature. We present the first study of deep optical broadband images of a gas-dominated blue compact dwarf galaxy (BCD) VCC848 (Mstar=2e8Msun) which hosts extended stellar shells and thus is confirmed to be a dwarf-dwarf merger. VCC848 is located in the outskirts of the Virgo Cluster. By analyzing the stellar light distribution, we found that VCC848 is the result of a merging between two dwarf galaxies with a primary-to-secondary mass ratio < ~ 5 for the stellar components and < ~ 2 for the presumed dark matter halos. The secondary progenitor galaxy has been almost entirely disrupted. The age-mass distribution of photometrically selected star cluster candidates in VCC848 implies that the cluster formation rate (CFR, proportional to star formation rate) was enhanced by a factor of ~ 7 - 10 during the past 1 Gyr. The merging-induced enhancement of CFR peaked near the galactic center a few hundred Myr ago and has started declining in the last few tens of Myr. The current star formation activities, as traced by the youngest clusters, mainly occur at large galactocentric distances (> ~ 1 kpc). The fact that VCC848 is still (atomic) gas-dominated after the period of most violent collision suggests that gas-rich dwarf galaxy merging can result in BCD-like remnants with extended atomic gas distribution surrounding a blue compact center, in general agreement with previous numerical simulations.

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.

CASA, the Common Astronomy Software Applications package, is the primary data processing software for the Atacama Large Millimeter/submillimeter Array (ALMA) and NSF's Karl G. Jansky Very Large Array (VLA), and is frequently used also for other radio telescopes. The CASA software can process data from both single-dish and aperture-synthesis telescopes, and one of its core functionalities is to support the data reduction and imaging pipelines for ALMA, VLA and the VLA Sky Survey (VLASS). CASA has recently undergone several exciting new developments, including an increased flexibility in Python (CASA 6), support of Very Long Baseline Interferometry (VLBI), performance gains through parallel imaging, data visualization with the new Cube Analysis Rendering Tool for Astronomy (CARTA), enhanced reliability and testing, and modernized documentation. These proceedings of the 2019 Astronomical Data Analysis Software & Systems (ADASS) conference give an update of the CASA project, and detail how these new developments will enhance user experience of CASA.

The Wide-Field InfraRed Survey Telescope (WFIRST) mission is expected to launch in the mid-2020s. Its weak lensing program is designed to enable unprecedented systematics control in photometric measurements, including shear recovery, point-spread function (PSF) correction, and photometric calibration. This will enable exquisite weak lensing science and allow us to adjust to and reliably contribute to the cosmological landscape after the initial years of observations from other concurrent Stage IV dark energy experiments. This potential requires equally careful planning and requirements validation as the mission prepares to enter its construction phase. We present a suite of image simulations based on GalSim that are used to construct a complex, synthetic WFIRST weak lensing survey that incorporates realistic input galaxies and stars, relevant detector non-idealities, and the current reference five-year WFIRST survey strategy. We present a first study to empirically validate the existing WFIRST weak lensing requirements flowdown using a suite of 12 matched image simulations, each representing a different perturbation to the wavefront or image motion model. These are chosen to induce a range of potential static and low- and high-frequency time-dependent PSF model errors. We analyze the measured shapes of galaxies from each of these simulations and compare them to a reference, fiducial simulation to infer the response of the shape measurement to each of these modes in the wavefront model. We then compare this to existing analytic flowdown requirements, and find general agreement between the empirically derived response and that predicted by the analytic model.

We investigate at a high angular resolution the spatial and kinematic structure of the S255IR high mass star-forming region, which demonstrated recently the first disk-mediated accretion burst in the massive young stellar object. The observations were performed with ALMA in Band 7 at an angular resolution $ \sim 0.1^{\prime\prime}$, which corresponds to $ \sim 180 $ AU. The 0.9 mm continuum, C$^{34}$S(7-6) and CCH $N=4-3$ data show a presence of very narrow ($ \sim 1000 $ AU), very dense ($n\sim 10^7$ cm$^{-3}$) and warm filamentary structures in this area. At least some of them represent apparently dense walls around the high velocity molecular outflow with a wide opening angle from the S255IR-SMA1 core, which is associated with the NIRS3 YSO. This wide-angle outflow surrounds a narrow jet. At the ends of the molecular outflow there are shocks, traced in the SiO(8-7) emission. The SiO abundance there is enhanced by at least 3 orders of magnitude. The CO(3-2) and SiO(8-7) data show a collimated and extended high velocity outflow from another dense core in this area, SMA2. The outflow is bent and consists of a chain of knots, which may indicate periodic ejections possibly arising from a binary system consisting of low or intermediate mass protostars. The C$^{34}$S emission shows evidence of rotation of the parent core. Finally, we detected two new low mass compact cores in this area (designated as SMM1 and SMM2), which may represent prestellar objects.

We use the multiconfiguration Dirac-Hartree-Fock (MCDHF) method combined with the relativistic configuration interaction (RCI) approach (GRASP2K) to provide a consistent set of transition energies and radiative transition data for the lower $n =3$ states in all Cl-like ions of astrophysical importance, from \ion{Cr}{8} to \ion{Zn}{14}. We also provide excitation energies calculated for \mbox{Fe X} using the many-body perturbation theory (MBPT, implemented within FAC). The comparison of the present MCDHF results with MBPT and with the available experimental energies indicates that the theoretical excitation energies are highly accurate, with uncertainties of only a few hundred cm$^{-1}$. Detailed comparisons for Fe~X and Ni~XII highlight discrepancies in the experimental energies found in the literature. Several new identifications are proposed.

We study the long-term variability in the optical monitoring database of Ark~120, a nearby radio-quiet active galactic nucleus (AGN) at a distance of 143 Mpc (z=0.03271). We compiled the historical archival photometric and spectroscopic data since 1974 and conducted a new two-year monitoring campaign in 2015-2017, resulting in a total temporal baseline over four decades. The long-term variations in the optical continuum exhibit a wave-like pattern and the Hbeta integrated flux series varies with a similar behavior. The broad Hbeta profiles have asymmetric double peaks, which change strongly with time and tend to merge into a single peak during some epochs. The period in the optical continuum determined from various period-search methods is about 20 yr and the estimated false alarm probability with null hypothesis simulations is about 1*10^-3. The overall variations of the broad Hbeta profiles also follow the same period. However, the present database only covers two cycles of the suggested period, which strongly encourages continued monitoring to track more cycles and confirm the periodicity. Nevertheless, in light of the possible periodicity and the complicated Hbeta profile, Ark~120 is one candidate of the nearest radio-quiet AGNs with possible periodic variability, and it is thereby a potential candidate host for a sub-parsec supermassive black hole binary.

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.

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.

Extracting accurate cosmological information from galaxy-galaxy and galaxy-matter correlation functions on non-linear scales ($\lesssim 10 h^{-1} \mathrm{Mpc}$) requires cosmological simulations. Additionally, one has to marginalise over several nuisance parameters of the galaxy-halo connection. However, the computational cost of such simulations prohibits naive implementations of stochastic posterior sampling methods like Markov chain Monte Carlo (MCMC) that would require of order $\mathcal{O}(10^6)$ samples in cosmological parameter space. Several groups have proposed surrogate models as a solution: a so-called emulator is trained to reproduce observables for a limited number of realisations in parameter space. Afterwards, this emulator is used as a surrogate model in an MCMC analysis. Here, we demonstrate a different method called Cosmological Evidence Modelling (CEM). First, for each simulation, we calculate the Bayesian evidence marginalised over the galaxy-halo connection by repeatedly populating the simulation with galaxies. We show that this Bayesian evidence is directly related to the posterior probability of cosmological parameters. Finally, we build a physically motivated model for how the evidence depends on cosmological parameters as sampled by the simulations. We demonstrate the feasibility of CEM by using simulations from the Aemulus simulation suite and forecasting cosmological constraints from BOSS CMASS measurements of redshift-space distortions. Our analysis includes an exploration of how galaxy assembly bias affects cosmological inference. Overall, CEM has several potential advantages over the more common approach of emulating summary statistics, including the ability to easily marginalise over highly complex models of the galaxy-halo connection and greater accuracy, thereby reducing the number of simulations required.

Most models for the connection between galaxies and their haloes ignore the possibility that galaxy properties may be correlated with halo properties other than mass, a phenomenon known as galaxy assembly bias. Yet, it is known that such correlations can lead to systematic errors in the interpretation of survey data. At present, the degree to which galaxy assembly bias may be present in the real Universe, and the best strategies for constraining it remain uncertain. We study the ability of several observables to constrain galaxy assembly bias from redshift survey data using the decorated halo occupation distribution (dHOD), an empirical model of the galaxy--halo connection that incorporates assembly bias. We cover an expansive set of observables, including the projected two-point correlation function $w_{\mathrm{p}}(r_{\mathrm{p}})$, the galaxy--galaxy lensing signal $\Delta \Sigma(r_{\mathrm{p}})$, the void probability function $\mathrm{VPF}(r)$, the distributions of counts-in-cylinders $P(N_{\mathrm{CIC}})$, and counts-in-annuli $P(N_{\mathrm{CIA}})$, and the distribution of the ratio of counts in cylinders of different sizes $P(N_2/N_5)$. We find that despite the frequent use of the combination $w_{\mathrm{p}}(r_{\mathrm{p}})+\Delta \Sigma(r_{\mathrm{p}})$ in interpreting galaxy data, the count statistics, $P(N_{\mathrm{CIC}})$ and $P(N_{\mathrm{CIA}})$, are generally more efficient in constraining galaxy assembly bias when combined with $w_{\mathrm{p}}(r_{\mathrm{p}})$. Constraints based upon $w_{\mathrm{p}}(r_{\mathrm{p}})$ and $\Delta \Sigma(r_{\mathrm{p}})$ share common degeneracy directions in the parameter space, while combinations of $w_{\mathrm{p}}(r_{\mathrm{p}})$ with the count statistics are more complementary. Therefore, we strongly suggest that count statistics should be used to complement the canonical observables in future studies of the galaxy--halo connection.

We present new observations of C-band continuum emission and masers to assess high-mass ($>$8 $M_\odot$) star formation at early evolutionary phases in the inner 200 pc of the Central Molecular Zone (CMZ) of the Galaxy. The continuum observation is complete to free-free emission from stars above 10-11 $M_\odot$ in 91% of the covered area. We identify 104 compact sources in the continuum emission, among which five are confirmed ultracompact H II regions, 12 are candidates of ultra-compact H II regions, and the remaining 87 sources are mostly massive stars in clusters, field stars, evolved stars, pulsars, extragalactic sources, or of unknown nature that is to be investigated. We detect class II CH$_3$OH masers at 23 positions, among which six are new detections. We confirm six known H$_2$CO masers in two high-mass star forming regions, and detect two new H$_2$CO masers toward the Sgr C cloud, making it the ninth region in the Galaxy that contains masers of this type. In spite of these detections, we find that current high-mass star formation in the inner CMZ is only taking place in seven isolated clouds. The results suggest that star formation at early evolutionary phases in the CMZ is about 10 times less efficient than expected by the dense gas star formation relation, which is in line with previous studies that focus on more evolved phases of star formation. This means that if there will be any impending, next burst of star formation in the CMZ, it has not yet begun.

We investigate the influence of Wolf-Rayet (W-R) stars on their surrounding star-forming molecular clouds. We study five regions containing W-R stars in the inner Galactic plane ($l\sim$[14$^\circ$-52$^\circ$]), using multi-wavelength data from near-infrared to radio wavelengths. Analysis of $^{13}$CO line data reveals that these W-R stars have developed gas-deficient cavities in addition to molecular shells with expansion velocities of a few km s$^{-1}$. The pressure owing to stellar winds primarily drives these expanding shells and sweeps up the surrounding matter to distances of a few pc. The column densities of shells are enhanced by a minimum of 14% for one region to a maximum of 88% for another region with respect to the column densities within their central cavities. No active star formation - including molecular condensations, protostars, or ionized gas - is found inside the cavities, whereas such features are observed around the molecular shells. Although the expansion of ionized gas is considered an effective mechanism to trigger star formation, the dynamical ages of the HII regions in our sample are generally not sufficiently long to do so efficiently. Overall, our results hint at the possible importance of negative W-R wind-driven feedback on the gas-deficient cavities, where star formation is quenched as a consequence. In addition, the presence of active star formation around the molecular shells indicates that W-R stars may also assist in accumulating molecular gas, and that they could initiate star formation around those shells.

The chemical structure of high-mass star nurseries is important for a general understanding of star formation. Deuteration is a key chemical process in the earliest stages of star formation because its efficiency is sensitive to the environment. Using the IRAM-30 m telescope at 1.3--4.3 mm wavelengths, we have imaged two parsec-scale high-mass protostellar clumps (P1 and S) that show different evolutionary stages but are located in the same giant filamentary {infrared dark cloud} G28.34+0.06. Deep spectral images at subparsec resolution reveal the dust and gas physical structures of both clumps. We find that (1) the low-$J$ lines of $\rm N_2H^+$, HCN, HNC, and $\rm HCO^+$ isotopologues are subthermally excited; and (2) the deuteration of $\rm N_2H^+$ is more efficient than that of $\rm HCO^+$, HCN, and HNC by an order of magnitude. The deuterations of these species are enriched toward the chemically younger clump S compared with P1, indicating that this process favors the colder and denser environment ($\rm T_{kin} \sim14 K$, $\rm N(NH_3) \sim 9\times 10^{15}\,cm^{-2}$). In contrast, single deuteration of $\rm NH_3$ is insensitive to the environmental difference between P1 and S; and (3) single deuteration of $\rm CH_3OH$ ($\rm > 10\%$) is detected toward the location where CO shows a depletion of $\sim10$. This comparative chemical study between P1 and S links the chemical variations to the environmental differences and shows chemical similarities between the early phases of high- and low-mass star-forming regions.

A high-energy muon neutrino event, IceCube-170922A, was recently discovered in both spatial and temporal coincidence with a gamma-ray flare of the blazar TXS 0506+056. It has been shown, with standard one-zone models, that neutrinos can be produced in the blazar jet via hadronic interactions, but with a flux which is mostly limited by the X-ray data. In this work, we explore the neutrino production from TXS 0506+056 by invoking two physically distinct emission zones in the jet, separated by the broad line region (BLR). Using the Doppler-boosted radiation of the BLR as the target photon field, the inner zone accounts for the neutrino and gamma-ray emission via $p\gamma$ interactions and inverse Compton scattering respectively, while the outer zone produces the optical and X-ray emission via synchrotron and synchrotron self-Compton processes. The different conditions of the two zones allow us to suppress the X-ray emission from the electromagnetic cascade, and set a much higher upper limit on the muon neutrino flux (i.e., $\sim 10^{-11}\rm erg~cm^{-2}s^{-1}$) than in one-zone models. We compare, in detail, our scenario with one-zone models discussed in the literature, and argue that differentiating between such scenarios will become possible with next generation neutrino telescopes, such as IceCube-Gen2.

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.

We present molecular line and dust continuum observations of a Planck-detected cold cloud, G074.11+00.11. The cloud consists of a system of curved filaments and a central star-forming clump. The clump is associated with several infrared sources and H2O maser emission. We aim to determine the mass distribution and gas dynamics within the clump, to investigate if the filamentary structure seen around the clump repeats itself on a smaller scale, and to estimate the fractions of mass contained in dense cores and filaments. The velocity distribution of pristine dense gas can be used to investigate the global dynamical state of the clump, the role of filamentary inflows, filament fragmentation and core accretion. We use molecular line and continuum observations from single dish observatories and interferometric facilities to study the kinematics of the region. The molecular line observations show that the central clump may have formed as a result of a large-scale filament collision. The central clump contains three compact cores. Assuming a distance of 2.3 kpc, based on Gaia observations and a three-dimensional extinction method of background stars, the mass of the central clump exceeds 700 solar masses, which is roughly 25% of the total mass of the cloud. Our virial analysis suggests that the central clump and all identified substructures are collapsing. We find no evidence for small-scale filaments associated with the cores. Our observations indicate that the clump is fragmented into three cores with masses in the range [10,50] solar masses and that all three are collapsing. The presence of an H2O maser emission suggests active star formation. However, the CO lines show only weak signs of outflows. We suggest that the region is young and any processes leading to star formation have just recently begun.

We report a ground X-ray calibration of two X-ray telescope prototypes at the PANTER X-ray Test Facility, of the Max-Planck-Institute for Extraterrestrial Physics, in Neuried, Germany. The X-ray telescope prototypes were developed by the Institute of Precision Optical Engineering (IPOE) of Tongji University, in a conical Wolter-I configuration, using thermal glass slumping technology. Prototype #1 with 3 layers and Prototype #2 with 21 layers were tested to assess the prototypes' on-axis imaging performance. The measurement of Prototype #1 indicates a Half Power Diameter (HPD) of 82" at 1.49 keV. As for Prototype #2, we performed more comprehensive measurements of on-axis angular resolution and effective area at several energies ranging from 0.5-10 keV. The HPD and effective area are 111" and 39 cm^2 at 1.49 keV, respectively, at which energy the on-axis performance of the prototypes is our greatest concern.

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.

We compiled a sample of 73 active galactic nuclei (AGNs) with reverberation mapping (RM) observations from RM campaigns including our ongoing campaign of monitoring super-Eddington accreting massive black holes (SEAMBHs). This sample covers a large range of black hole (BH) mass $(M_{\bullet}=10^{6-9}~M_\odot)$, dimensionless accretion rates $(\dot{\mathscr{M}}=10^{-2.7}-10^{2.7})$ and 5100~\AA~luminosity $(L_{5100}=10^{42-46}~\rm erg~s^{-1})$, allowing us to systematically study the AGN variability and their relations with BH mass, accretion rates, and optical luminosity. We employed the damped random walk (DRW) model to delineate the optical variability of continuum at 5100~\AA~and obtained damped variability timescale ($\tau_{\rm d}$) and amplitude ($\sigma_{\rm d}$) using a Markov Chain Monte Carlo (MCMC) method. We also estimated the traditional variability amplitudes ($F_{\rm var}$), which provide a model-independent measure and therefore are used to test the DRW results. We found that AGN variability characteristics are generally correlated with $(M_{\bullet},\dot{\mathscr{M}},L_{5100})$. These correlations are smooth from sub-Eddington to super-Eddington accretion AGNs, probably implying that the AGN variability may be caused by the same physical mechanism.

We report results of the first reverberation mapping campaign of I Zwicky 1 during $2014$-$2016$, which showed unambiguous reverberations of the broad H$\beta$ line emission to the varying optical continuum. From analysis using several methods, we obtain a reverberation lag of $\tau_{\rm H\beta}=37.2^{+4.5}_{-4.9}\,$ days. Taking a virial factor of $f_{_{\rm BLR}}=1$, we find a black hole mass of $M_{\bullet}=9.30_{-1.38}^{+1.26}\times 10^6 M_{\odot}$ from the mean spectra. The accretion rate is estimated to be $203.9_{-65.8}^{+61.0}\,L_{\rm Edd}c^{-2}$, suggesting a super-Eddington accretor, where $L_{\rm Edd}$ is the Eddington luminosity and $c$ is the speed of light. By decomposing {\it Hubble Space Telescope} images, we find that the stellar mass of the bulge of its host galaxy is $\log (M_{\rm bulge}/M_{\odot}) = \rm 10.92\pm 0.07$. This leads to a black hole to bulge mass ratio of $\sim 10^{-4}$, which is significantly smaller than that of classical bulges and elliptical galaxies. After subtracting the host contamination from the observed luminosity, we find that I Zw 1 follows the empirical $R_{\rm BLR}\propto L_{5100}^{1/2}$ relation.

Using {\sl Gaia} Data Release 2 photometry, we report the detection of extended main-sequence turnoff (eMSTO) regions in the color--magnitude diagrams (CMDs) of the $\sim 14$ Myr-old double clusters $h$ and $\chi$ Persei (NGC 869 and NGC 884). We find that stars with masses below $\sim$1.3 $M_{\odot}$ in both $h$ and $\chi$ Persei populate narrow main sequences (MSs), while more massive stars define the eMSTO, closely mimicking observations of young Galactic and Magellanic Cloud clusters (with ages older than $\sim$30 Myr). Previous studies based on clusters older than $\sim$30 Myr find that rapidly rotating MS stars are redder than slow rotators of similar luminosity, suggesting that stellar rotation may be the main driver of the eMSTO. By combining photometry and projected rotational velocities from the literature of stars in $h$ and $\chi$ Persei, we find no obvious relation between the rotational velocities and colors of non-emission-line eMSTO stars, in contrast with what is observed in older clusters. Similarly to what is observed in Magellanic Cloud clusters, most of the extremely rapidly rotating stars, identified by their strong H$\alpha$ emission lines, are located in the red part of the eMSTOs. This indicates that stellar rotation plays a role in the color and magnitude distribution of MSTO stars. By comparing the observations with simulated CMDs, we find that a simple population composed of coeval stars that span a wide range of rotation rates is unable to reproduce the color spread of the clusters' MSs. We suggest that variable stars, binary interactions, and stellar rotation affect the eMSTO morphology of these very young clusters.

Despite many decades of study, the kinematics of the broad-line region of 3C~273 are still poorly understood. We report a new, high signal-to-noise, reverberation mapping campaign carried out from November 2008 to March 2018 that allows the determination of time lags between emission lines and the variable continuum with high precision. The time lag of variations in H$\beta$ relative to those of the 5100 Angstrom continuum is $146.8_{-12.1}^{+8.3}$ days in the rest frame, which agrees very well with the Paschen-$\alpha$ region measured by the GRAVITY at The Very Large Telescope Interferometer. The time lag of the H$\gamma$ emission line is found to be nearly the same as for H$\beta$. The lag of the Fe II emission is $322.0_{-57.9}^{+55.5}$ days, longer by a factor of $\sim$2 than that of the Balmer lines. The velocity-resolved lag measurements of the H$\beta$ line show a complex structure which can be possibly explained by a rotation-dominated disk with some inflowing radial velocity in the H$\beta$-emitting region. Taking the virial factor of $f_{\rm BLR} = 1.3$, we derive a BH mass of $M_{\bullet} = 4.1_{-0.4}^{+0.3} \times 10^8 M_{\odot}$ and an accretion rate of $9.3\,L_{\rm Edd}\,c^{-2}$ from the H$\beta$ line. The decomposition of its $HST$ images yields a host stellar mass of $M_* = 10^{11.3 \pm 0.7} M_\odot$, and a ratio of $M_{\bullet}/M_*\approx 2.0\times 10^{-3}$ in agreement with the Magorrian relation. In the near future, it is expected to compare the geometrically-thick BLR discovered by the GRAVITY in 3C 273 with its spatially-resolved torus in order to understand the potential connection between the BLR and the torus.

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.

Context. We present the first high-cadence multi-wavelength radial-velocity observations of the Sun-as-a-star, carried out during 57 consecutive days using the stellar \'echelle spectrograph at the Hertzsprung SONG Telescope operating at the Teide Observatory. Aims. The aim was to produce a high-quality data set and reference values for the global helioseismic parameters {\nu_{max}}, and {\Delta \nu} of the solar p-modes using the SONG instrument. The obtained data set or the inferred values should then be used when the scaling relations are applied to other stars showing solar-like oscillations which are observed with SONG or similar instruments. Methods. We used different approaches to analyse the power spectrum of the time series to determine {\nu_{max}}; simple Gaussian fitting and heavy smoothing of the power spectrum. {\Delta\nu} was determined using the method of autocorrelation of the power spectrum. The amplitude per radial mode was determined using the method described in Kjeldsen et al. (2008). Results. We found the following values for the solar oscillations using the SONG spectrograph: {\nu_{max}} = 3141 {\pm} 12 {\mu}Hz, {\Delta\nu} = 134.98 {\pm} 0.04 {\mu}Hz and an average amplitude of the strongest radial modes of 16.6 {\pm} 0.4 cm/s. These values are consistent with previous measurements with other techniques.

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.

Although many high-energy neutrinos detected by the IceCube telescope are believed to have anextraterrestrial origin, their astrophysical sources remain a mystery. Recently, an unprecedenteddiscovery of a high-energy muon neutrino event coincident with a multiwavelength flare from ablazar, TXS 0506+056, shed some light on the origin of the neutrinos. It is usually believed that ablazar is produced by a relativistic jet launched from an accreting supermassive black hole (SMBH).Here we show that the high-energy neutrino event can be interpreted by the inelastic hadronuclearinteractions between the accelerated cosmic-ray protons in the relativistic jet and the dense gasclouds in the vicinity of the SMBH. Such a scenario only requires a moderate proton power in thejet, which could be much smaller than that required in the conventional hadronic model whichinstead calls upon the photomeson process. Meanwhile, the flux of the multiwavelength flare fromthe optical to gamma-ray band can be well explained by invoking a second radiation zone in thejet at a larger distance to the SMBH. In our model, the neutrino emission lasts a shorter time thanthe multiwavelength flare so the neutrino event is not necessarily correlated with the flare but it is probably accompanied by a spectrum hardening above a few GeV.

The infrared dark clouds (IRDCs) G11.11$-$0.12 and G28.34$+$0.06 are two of the best-studied IRDCs in our Galaxy. These two clouds host clumps at different stages of evolution, including a massive dense clump in both clouds that is dark even at 70 and 100$\mu$m. Such seemingly quiescent massive dense clumps have been speculated to harbor cores that are precursors of high-mass stars and clusters. We observed these two "prestellar" regions at 1mm with the Submillimeter Array (SMA) with the aim of characterizing the nature of such cores. We show that the clumps fragment into several low- to high-mass cores within the filamentary structure of the enveloping cloud. However, while the overall physical properties of the clump may indicate a starless phase, we find that both regions host multiple outflows. The most massive core though 70 $\mu$m dark in both clumps is clearly associated with compact outflows. Such low-luminosity, massive cores are potentially the earliest stage in the evolution of a massive protostar. We also identify several outflow features distributed in the large environment around the most massive core. We infer that these outflows are being powered by young, low-mass protostars whose core mass is below our detection limit. These findings suggest that low-mass protostars have already formed or are coevally formed at the earliest phase of high-mass star formation.

We present long-term, multi-color photometric monitoring and spectroscopic observations of the short-period M-type eclipsing binary BX Tri. Six flare events were recorded in 4 years from 2014 to 2017. Three of them were detected on one night within an orbital cycle. The strongest one was identified on December 23, 2014. With the amplitudes of $\Delta$B=0.48 mag, $\Delta$V=0.28 mag, $\Delta$R=0.10 mag, $\Delta$I=0.02 mag, the total energy due to the flare event was measured to be 4.08 ($\pm$0.24) $\times$$10^{34}$ erg, exceeding the superflare level ($10^{34}$). Based on the observations, the evolutionary status of the binary system as well as the long-term orbital period changes were analyzed. It reveals that BX Tri is probably a semi-detached system with the primary component filling its Roche lobe. The extremely high occurrence of flare events of the binary could be related to the rapid mass transfer between components.

The dense wind environment (or circumstellar medium) may be ubiquitous for the regular Type II supernovae (SNe) before the explosion, the interaction of which with the SN ejecta could result in a wind breakout event. The shock generated by the interaction of the SN ejecta and the wind can accelerate the protons and subsequently the high-energy gamma-rays and neutrinos could arise from the inelastic pp collisions. In this work, we present the detailed calculations of gamma-ray and neutrino production for the regular Type II SNe. The calculation is executed by applying time-dependent evolutions of dynamic and proton distribution so that the emission could be shown at different times. Our results show, for the SN 2013fs-like wind environment, the multi-GeV and ~ few-100 TeV gamma-rays are detectable with a time window of several days at <~ 2-3 Mpc by Fermi/LAT and CTA during the ejecta-wind interaction, respectively, and can be detected at a further distance if the wind environment is denser. Besides, we found the contribution of the wind breakouts of regular Type II SNe to diffuse neutrino flux is subdominant by assuming all Type II SNe are SN 2013fs-like, whereas for a denser wind environment the contribution could be conspicuous above 300TeV.

To probe the initial conditions for high-mass star and cluster formation, we investigate the properties of dense filaments within the infrared dark cloud G035.39-00.33 (IRDC G035.39) in a combined Very Large Array (VLA) and the Green Bank Telescope (GBT) mosaic tracing the NH3 (1,1) and (2,2) emission down to 0.08 pc scales. Using agglomerative hierarchical clustering on multiple line-of-sight velocity component fitting results, we identify seven extended velocity-coherent components in our data, likely representing spatially coherent physical structures, some exhibiting complex gas motions. The velocity gradient magnitude distribution peaks at its mode of 0.35 km/s/pc and has a long tail extending into higher values of 1.5 - 2 km/s/pc, and is generally consistent with those found toward the same cloud in other molecular tracers and with the values found towards nearby low-mass dense cloud cores at the same scales. Contrary to observational and theoretical expectations, we find the non-thermal ammonia line widths to be systematically narrower (by about 20%) than those of N2H+ (1-0) line transition observed with similar resolution. If the observed ordered velocity gradients represent the core envelope solid-body rotation, we estimate the specific angular momentum to be about 2 x 10^21 cm^2/s, similar to the low-mass star-forming cores. Together with the previous finding of subsonic motions in G035.39, our results demonstrate high levels of similarity between kinematics of a high-mass star-forming IRDC and the low-mass star formation regime.

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

We present new results on the relationship between central galaxies and dark matter haloes inferred from observations of satellite kinematics in the Sloan Digital Sky Survey (SDSS) DR7. We employ an updated analysis framework that includes detailed mock catalogues to model observational effects in SDSS. Our results constrain the colour-dependent conditional luminosity function (CLF) of dark matter haloes, as well as the radial profile of satellite galaxies. Confirming previous results, we find that red central galaxies live in more massive haloes than blue galaxies at fixed luminosity. Additionally, our results suggest that satellite galaxies have a radial profile less centrally concentrated than dark matter but not as cored as resolved subhaloes in dark matter-only simulations. Compared to previous works using satellite kinematics by More et al., we find much more competitive constraints on the galaxy-halo connection, on par with those derived from a combination of galaxy clustering and galaxy-galaxy lensing. We compare our results on the galaxy-halo connection to other studies using galaxy clustering and group catalogues, showing very good agreement between these different techniques. We discuss future applications of satellite kinematics in the context of constraining cosmology and the relationship between galaxies and dark matter haloes.

We have started a long-term reverberation mapping project using the Wyoming Infrared Observatory 2.3 meter telescope titled "Monitoring AGNs with H\beta\ Asymmetry" (MAHA). The motivations of the project are to explore the geometry and kinematics of the gas responsible for complex H\beta\ emission-line profiles, ideally leading to an understanding of the structures and origins of the broad-line region (BLR). Furthermore, such a project provides the opportunity to search for evidence of close binary supermassive black holes. We describe MAHA and report initial results from our first campaign, from December 2016 to May 2017, highlighting velocity-resolved time lags for four AGNs with asymmetric H\beta\ lines. We find that 3C 120, Ark 120, and Mrk 6 display complex features different from the simple signatures expected for pure outflow, inflow, or a Keplerian disk. While three of the objects have been previously reverberation mapped, including velocity-resolved time lags in the cases of 3C 120 and Mrk 6, we report a time lag and corresponding black hole mass measurement for SBS 1518+593 for the first time. Furthermore, SBS 1518+593, the least asymmetric of the four, does show velocity-resolved time lags characteristic of a Keplerian disk or virialized motion more generally. Also, the velocity-resolved time lags of 3C 120 have significantly changed since previously observed, indicating an evolution of its BLR structure. Future analyses of the data for these objects and others in MAHA will explore the full diversity of H\beta\ lines and the physics of AGN BLRs.

The kinematics of satellite galaxies moving in a dark matter halo are a direct probe of the underlying gravitational potential. Thus, the phase-space distributions of satellites represent a powerful tool to determine the galaxy-halo connection from observations. By stacking the signal of a large number of satellite galaxies this potential can be unlocked even for haloes hosting a few satellites on average. In this work, we test the impact of various modelling assumptions on constraints derived from analysing satellite phase-space distributions in the non-linear, 1-halo regime. We discuss their potential to explain the discrepancy between average halo masses derived from satellite kinematics and gravitational lensing previously reported. Furthermore, we develop an updated, more robust analysis to extract constraints on the galaxy-halo relation from satellite properties in spectroscopic galaxy surveys such as the SDSS. We test the accuracy of this approach using a large number of realistic mock catalogues. Furthermore, we find that constraints derived from such an analysis are complementary and competitive with respect to the commonly used galaxy clustering and galaxy-galaxy lensing observables.

NGC 5548 has been well spectroscopically monitored for reverberation mapping of the central kinematics by 19 campaigns. Using the maximum entropy method in this Letter, we build up a high-quality velocity-delay map of the H$\beta$ emission line in the light curves of the continuum and the line variations observed between 2015-2016. The map shows the response strength and lags of the velocity fields of the H$\beta$ emitting regions. The velocity-delay structure of the map is generally symmetric, with strong red and blue wings at time lag $\tau \leq 15$ days, a narrower velocity distribution at $\tau \geq 15$ days, and a deficit of response in the core. This is suggestive of a disk geometry of the broad-line region (BLR). The relatively weaker H$\beta$ response at the longer lags in the red side indicates anisotropic emission from the outer part of the BLR. We also recover the velocity-delay maps of NGC 5548 from the historical data of 13 years to investigate the long-term variability of its BLR. In general, the BLR of NGC 5548 was switching between the inflow and virialized phases in the past years. The resultant maps of seven years reveal inflow signatures and show decreasing lags, indicating that the changes in the BLR size are related to the infalling BLR gas. The other four maps show potential disk signatures which are similar to our map.

Recently, a high-energy neutrino event IceCube-170922A in the spatial and temporal coincidence with the flaring gamma-ray blazar TXS 0506+056 was reported. A neutrino outburst between September 2014 and March 2015 was discovered in the same direction by a further investigation of $9.5$ years of IceCube data, while the blazar is in a quiescent state during the outburst with a gamma-ray flux only about one-fifth of the neutrino flux. In this letter, we propose the neutrino outburst originates from the interaction between a relativistic jet and a dense gas cloud which may be formed via the tidally disrupted envelope of a red giant being blown by the impact of the jet. Gamma-ray photons and electron/positron pairs that are produced correspondingly will induce electromagnetic cascades. Comptonization of the cascade emission inside the cloud forms an X-ray photon field with Wien distribution. GeV flux is suppressed due to the absorption by the Comptonized photon field and, as a result, a hard spectrum above 10 GeV is formed. The gamma-ray spectrum predicted in our model is consistent with the Fermi-LAT data of TXS 0506+056.

We report the discovery of a new eclipsing EL CVn-type binary, consisting of a $\delta$ Sct-type pulsator and a thermally bloated low-mass pre-He white dwarf (WD). Spectroscopy and time-series $BV$ photometry of V1224 Cas were carried out. The spectroscopy reveals a spectral type of A3 for the star. Light-curve modelling indicate that V1224 Cas is a short-period detached system containing a possible low-mass WD with an effective temperature of about 9516 K. Based on the effective temperature and the surface gravity of the A-type primary star from the spectroscopic results, the absolute parameters of the components were estimated as: $M_{P}$=2.16$\pm$0.22$M_{\sun}$, $R_{P}$=3.54$\pm$0.12$R_{\sun}$, $L_{P}$=55.9$\pm$6.9$L_{\sun}$, $M_{S}$=0.19$\pm$0.02$M_{\sun}$, $R_{S}$=0.97$\pm$0.04$R_{\sun}$, and $L_{S}$=6.9$\pm$0.6$L_{\sun}$. We therefore introduce V1224 Cas as a new EL CVn-type binary candidate. The light curves in both filters all show multi-periodic pulsations, superimposed on binary effects. We performed a preliminary frequency analysis of the light residuals after removing the synthetic eclipsing curve from the original observational data. The results suggest that the rapid light variations among the light curves could be attributed to the $\delta$ Sct-type primary component. We therefore conclude that V1224 Cas is very likely a WD+$\delta$ Sct binary.

By using a non-local and time-dependent convection theory, we have calculated radial and low-degree non-radial oscillations for stellar evolutionary models with $M=1.4$--3.0\,$\mathrm{M}_\odot$. The results of our study predict theoretical instability strips for $\delta$ Scuti and $\gamma$ Doradus stars, which overlap with each other. The strip of $\gamma$ Doradus is slightly redder in colour than that of $\delta$ Scuti. We have paid great attention to the excitation and stabilization mechanisms for these two types of oscillations, and we conclude that radiative $\kappa$ mechanism plays a major role in the excitation of warm $\delta$ Scuti and $\gamma$ Doradus stars, while the coupling between convection and oscillations is responsible for excitation and stabilization in cool stars. Generally speaking, turbulent pressure is an excitation of oscillations, especially in cool $\delta$ Scuti and $\gamma$ Doradus stars and all cool Cepheid- and Mira-like stars. Turbulent thermal convection, on the other hand, is a damping mechanism against oscillations that actually plays the major role in giving rise to the red edge of the instability strip. Our study shows that oscillations of $\delta$ Scuti and $\gamma$ Doradus stars are both due to the combination of $\kappa$ mechanism and the coupling between convection and oscillations, and they belong to the same class of variables at the low-luminosity part of the Cepheid instability strip. Within the $\delta$ Scuti--$\gamma$ Doradus instability strip, most of the pulsating variables are very likely hybrids that are excited in both p and g modes.

We present the distance priors from the finally released $Planck~ \textrm{TT,TE,EE}+\textrm{lowE}$ data in 2018. The uncertainties are around $40\%$ smaller than those from $Planck$ 2015 TT$+$lowP. In order to check the validity of these new distance priors, we adopt the distance priors to constrain the cosmological parameters in different dark energy models, including the $\Lambda$CDM model, the $w$CDM model and the CPL model, and conclude that the distance priors provide consistent constraints on the relevant cosmological parameters compared to those from the full $Planck$ 2018 data release.

We present high resolution (0.2", 1000 AU) 1.3 mm ALMA observations of massive infrared dark cloud clump, G028.37+00.07-C1, thought to harbor the early stages of massive star formation. Using $\rm N_2D^+$(3-2) we resolve the previously identified C1-S core, separating the bulk of its emission from two nearby protostellar sources. C1-S is thus identified as a massive ($\sim50\:M_\odot$), compact ($\sim0.1\:$pc diameter) starless core, e.g., with no signs of outflow activity. Being highly deuterated, this is a promising candidate for a pre-stellar core on the verge of collapse. An analysis of its dynamical state indicates a sub-virial velocity dispersion compared to a trans-Alfv\'enic turbulent core model. However, virial equilibrium could be achieved with sub-Alfv\'enic conditions involving $\sim2\:$mG magnetic field strengths.

Temporal photometric variations at near infrared to submillimeter wavelengths have been found in low-mass young stellar objects. These phenomena are generally interpreted as accretion events of star-disk systems with varying accretion rates. There is growing evidence suggesting that similar luminosity flaring also occurs in high-mass star/cluster-forming regions. We report in this Letter the rise and fall of the 900 ${\mu}$m continuum emission and the newly found 349.1 GHz methanol maser emission in the massive star forming region S255IR~SMA1 observed with the Submillimeter Array and the Atacama Large Millimeter/submillimeter Array. The level of flux variation at a factor of $\sim$ 2 at the submillimeter band and the relatively short 2-year duration of this burst suggest that the event is probably similar to those milder and more frequent minor bursts seen in 3D numerical simulations.

We numerically integrate the Einstein's equations for a spatially flat Friedmann-Lemaire-Robertson-Walker (FLRW) background spacetime with a spatial curvature perturbation and evolving primordial tensor perturbations using the Einstein Toolkit. We find that although the primordial tensor perturbation doesn't play an important role in the evolution of the overdensity produced by the scalar perturbation, there is an obvious imprint left by the primordial tensor perturbation on the distribution of the fractional density perturbation in the nonlinear region. This imprint may be a possible probe of a gravitational waves background in the future.

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.

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.

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.

The high-energy (> 100MeV) emission observed by Fermi-LAT during the prompt phase of some luminous gamma-ray bursts (GRBs) could arise from the cascade induced by interactions between accelerated protons and the radiation field of GRBs. The photomeson process, which is usually suggested to operate in such a hadronic explanation, requires a rather high proton energy (> 10^17eV) for an efficient interaction. However, whether GRBs can accelerate protons to such a high energy is far from guaranteed, although they have been suggested as the candidate source for ultrahigh-energy cosmic rays. In this work, we revisit the hadronic model for the prompt high-energy emission of GRBs with a smaller maximum proton energy than the usually adopted value estimated from the Bohm condition. In this case, the Bethe-Heitler pair production process becomes comparably important or even dominates over the photomeson process. We show that with a relatively low maximum proton energy with a Lorentz factor of 10^5 in the comoving frame, the cascade emission can still reproduce various types of high-energy spectrum of GRBs. For most GRBs without high-energy emission detected, the maximum proton energy could be even lower and relax the constraints on the parameters of GRB jet resulting from the fact of non-detection of GRB neutrinos by IceCube.

High-mass star forming regions are typically thought to be dominated by supersonic motions. We present combined Very Large Array and Green Bank Telescope (VLA+GBT) observations of NH$_3$ (1,1) and (2,2) in the infrared dark cloud (IRDC) G035.39-00.33, tracing cold and dense gas down to scales of 0.07 pc. We find that, in contrast to previous similar studies of IRDCs, more than a third of the fitted ammonia spectra show subsonic non-thermal motions (mean line width of 0.71 $\mathrm{km~s^{-1}}$), and the sonic Mach number distribution peaks around $\mathcal{M} = 1$. As possible observational and instrumental biases would only broaden the line profiles, our results provide strong upper limits to the actual value of $\mathcal{M}$, further strengthening our findings of narrow line widths. This finding calls for a reevaluation of the role of turbulent dissipation and subsonic regions in massive-star and cluster formation. Based on our findings in G035.39, we further speculate that the coarser spectral resolution used in the previous VLA NH$_3$ studies may have inhibited the detection of subsonic turbulence in IRDCs. The reduced turbulent support suggests that dynamically important magnetic fields of the 1 mG order would be required to support against possible gravitational collapse. Our results offer valuable input into the theories and simulations that aim to recreate the initial conditions of high-mass star and cluster formation.

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.

As one of the series of papers reporting on a large reverberation mapping campaign of super-Eddington accreting massive black holes (SEAMBHs) in active galactic nuclei (AGNs), we present the results of 10 SEAMBHs monitored spectroscopically during 2015-2017. Six of them are observed for the first time, and have generally higher 5100 \AA\ luminosities than the SEAMBHs monitored in our campaign from 2012 to 2015; the remaining four are repeat observations to check if their previous lags change. Similar to the previous SEAMBHs, the H$\beta$ time lags of the newly observed objects are shorter than the values predicted by the canonical $R_{\mathrm{H\beta}}$-$L_{5100}$ relation of sub-Eddington AGNs, by factors of $\sim2-6$, depending on the accretion rate. The four previously observed objects have lags consistent with previous measurements. We provide linear regressions for the $R_{\mathrm{H\beta}}$-$L_{5100}$ relation, solely for the SEAMBH sample and for low-accretion AGNs. We find that the relative strength of Fe II and the profile of the H$\beta$ emission line can be used as proxies of accretion rate, showing that the shortening of H$\beta$ lags depends on accretion rates. The recent SDSS-RM discovery of shortened H$\beta$ lags in AGNs with low accretion rates provides compelling evidence for retrograde accretion onto the black hole. These evidences show that the canonical $R_{\mathrm{H\beta}}$-$L_{5100}$ relation holds only in AGNs with moderate accretion rates. At low accretion rates, it should be revised to include the effects of black hole spin, whereas the accretion rate itself becomes a key factor in the regime of high accretion rates.

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.

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.

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.

In this paper we calculate the radio burst signals from three kinds of structures of superconducting cosmic strings. By taking into account the observational factors including scattering and relativistic effects, we derive the event rate of radio bursts as a function of redshift with the theoretical parameters $G\mu$ and $\mathcal{I}$ of superconducting strings. Our analyses show that cusps and kinks may have noticeable contributions to the event rate and in most cases cusps would dominate the contribution, while the kink-kink collisions tend to have secondary effects. By fitting theoretical predictions with the normalized data of fast radio bursts, we for the first time constrain the parameter space of superconducting strings and report that the parameter space of $G\mu \sim [10^{-14}, 10^{-12}]$ and $\mathcal{I} \sim [10^{-1}, 10^{2}] ~ \rm{GeV}$ fit the observation well although the statistic significance is low due to the lack of observational data. Moreover, we derive two types of best fittings, with one being dominated by cusps with a redshift $z = 1.3$, and the other dominated by kinks at the range of the maximal event rate.

We update the constraints on the cosmological parameters by adopting the Planck data released in 2015 and Baryon Acoustic Oscillation (BAO) measurements including the new DR14 quasar sample measurement at redshift $z=1.52$, and we conclude that the based six-parameter $\Lambda$CDM model is preferred. Exploring some extensions to the $\Lambda$CDM models, we find that the equation of state of dark energy reads $w=-1.036\pm 0.056$ in the $w$CDM model, the effective relativistic degrees of freedom in the Universe is $N_\text{eff}=3.09_{-0.20}^{+0.18}$ in the $N_\text{eff}+\Lambda$CDM model and the spatial curvature parameter is $\Omega_k=(1.8\pm1.9)\times 10^{-3}$ in the $\Omega_k+\Lambda$CDM model at $68\%$ confidence level (C.L.), and the $95\%$ C.L. upper bounds on the sum of three active neutrinos masses are $\sum m_\nu<0.16$ eV for the normal hierarchy (NH) and $\sum m_\nu<0.19$ eV for the inverted hierarchy (IH) with $\Delta\chi^2\equiv \chi^2_\text{NH}-\chi^2_\text{IH}=-1.25$.

Dark matter halo clustering depends not only on halo mass, but also on other properties such as concentration and shape. This phenomenon is known broadly as assembly bias. We explore the dependence of assembly bias on halo definition, parametrized by spherical overdensity parameter, $\Delta$. We summarize the strength of concentration-, shape-, and spin-dependent halo clustering as a function of halo mass and halo definition. Concentration-dependent clustering depends strongly on mass at all $\Delta$. For conventional halo definitions ($\Delta \sim 200\mathrm{m}-600\mathrm{m}$), concentration-dependent clustering at low mass is driven by a population of haloes that is altered through interactions with neighbouring haloes. Concentration-dependent clustering can be greatly reduced through a mass-dependent halo definition with $\Delta \sim 20\mathrm{m}-40\mathrm{m}$ for haloes with $M_{200\mathrm{m}} \lesssim 10^{12}\, h^{-1}\mathrm{M}_{\odot}$. Smaller $\Delta$ implies larger radii and mitigates assembly bias at low mass by subsuming altered, so-called backsplash haloes into now larger host haloes. At higher masses ($M_{200\mathrm{m}} \gtrsim 10^{13}\, h^{-1}\mathrm{M}_{\odot}$) larger overdensities, $\Delta \gtrsim 600\mathrm{m}$, are necessary. Shape- and spin-dependent clustering are significant for all halo definitions that we explore and exhibit a relatively weaker mass dependence. Generally, both the strength and the sense of assembly bias depend on halo definition, varying significantly even among common definitions. We identify no halo definition that mitigates all manifestations of assembly bias. A halo definition that mitigates assembly bias based on one halo property (e.g., concentration) must be mass dependent. The halo definitions that best mitigate concentration-dependent halo clustering do not coincide with the expected average splashback radii at fixed halo mass.

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.

We present historical VLA and Jansky Very Large Array multi-frequency multi-array radio continuum imaging of a unique sample of 100 radio sources that have been selected on the basis of low axial ratios. These observations allow us the opportunity to study radio sources with synchrotron plasma that is significantly offset from the main radio axis and therefore open a window into investigations of physical mechanisms responsible for depositing the plasma in off-axis regions. These images are discussed in detail in subsequent papers in this series (Saripalli & Roberts 2017; Roberts & Saripalli 2017).

Aims. Infrared dark clouds represent the earliest stages of high-mass star formation. Detailed observations of their physical conditions on all physical scales are required to improve our understanding of their role in fueling star formation. Methods. We investigate the large-scale structure of the IRDC G035.39-00.33, probing the dense gas with the classical ammonia thermometer. This allows us to put reliable constraints on the temperature of the extended, pc-scale dense gas reservoir and to probe the magnitude of its non-thermal motions. Available far-infrared observations can be used in tandem with the observed ammonia emission to estimate the total gas mass contained in G035.39-00.33. Results. We identify a main velocity component as a prominent filament, manifested as an ammonia emission intensity ridge spanning more than 6 pc, consistent with the previous studies on the Northern part of the cloud. A number of additional line-of-sight components are found, and a large scale, linear velocity gradient of ~0.2 km s$^{-1}$ pc$^{-1}$ is found along the ridge of the IRDC. In contrast to the dust temperature map, an ammonia-derived kinetic temperature map, presented for the entirety of the cloud, reveals local temperature enhancements towards the massive protostellar cores. We show that without properly accounting for the line of sight contamination, the dust temperature is 2-3 K larger than the gas temperature measured with NH$_3$. Conclusions. While both the large scale kinematics and temperature structure are consistent with that of starless dark filaments, the kinetic gas temperature profile on smaller scales is suggestive of tracing the heating mechanism coincident with the locations of massive protostellar cores.

High-mass stars have generally been assumed to accrete most of their mass while already contracted onto the main sequence, but this hypothesis has not been observationally tested. We present ALMA observations of a 3 x 1.5 pc area in the W51 high-mass star-forming complex. We identify dust continuum sources and measure the gas and dust temperature through both rotational diagram modeling of CH3OH and brightness-temperature-based limits. The observed region contains three high-mass YSOs that appear to be at the earliest stages of their formation, with no signs of ionizing radiation from their central sources. The data reveal high gas and dust temperatures (T > 100 K) extending out to about 5000 AU from each of these sources. There are no clear signs of disks or rotating structures down to our 1000 AU resolution. The extended warm gas provides evidence that, during the process of forming, these high-mass stars heat a large volume and correspondingly large mass of gas in their surroundings, inhibiting fragmentation and therefore keeping a large reservoir available to feed from. By contrast, the more mature massive stars that illuminate compact H ii regions have little effect on their surrounding dense gas, suggesting that these main sequence stars have completed most or all of their accretion. The high luminosity of the massive protostars (L > 10^4 solar luminosities), combined with a lack of centimeter continuum emission from these sources, implies that they are not on the main sequence while they accrete the majority of their mass; instead, they may be bloated and cool.

The early reionization (ERE) is supposed to be a physical process which happens after recombination, but before the instantaneous reionization caused by the first generation of stars. We investigate the effect of the ERE on the temperature and polarization power spectra of cosmic microwave background (CMB), and adopt principal components analysis (PCA) to model-independently reconstruct the ionization history during the ERE. In addition, we also discuss how the ERE affects the cosmological parameter estimates, and find that the ERE does not impose any significant influences on the tensor-to-scalar ratio $r$ and the neutrino mass at the sensitivities of current experiments. The better CMB polarization data can be used to give a tighter constraint on the ERE and might be important for more precisely constraining cosmological parameters in the future.

The formation of deuterated molecules is favoured at low temperatures and high densities. Therefore, the deuteration fraction D$_{frac}$ is expected to be enhanced in cold, dense prestellar cores and to decrease after protostellar birth. Previous studies have shown that the deuterated forms of species such as N2H+ (formed in the gas phase) and CH3OH (formed on grain surfaces) can be used as evolutionary indicators and to constrain their dominant formation processes and time-scales. Formaldehyde (H2CO) and its deuterated forms can be produced both in the gas phase and on grain surfaces. However, the relative importance of these two chemical pathways is unclear. Comparison of the deuteration fraction of H2CO with respect to that of N2H+, NH3 and CH3OH can help us to understand its formation processes and time-scales. With the new SEPIA Band 5 receiver on APEX, we have observed the J=3-2 rotational lines of HDCO and D2CO at 193 GHz and 175 GHz toward three massive star forming regions hosting objects at different evolutionary stages: two High-mass Starless Cores (HMSC), two High-mass Protostellar Objects (HMPOs), and one Ultracompact HII region (UCHII). By using previously obtained H2CO J=3-2 data, the deuteration fractions HDCO/H2CO and D2CO/HDCO are estimated. Our observations show that singly-deuterated H2CO is detected toward all sources and that the deuteration fraction of H2CO increases from the HMSC to the HMPO phase and then sharply decreases in the latest evolutionary stage (UCHII). The doubly-deuterated form of H2CO is detected only in the earlier evolutionary stages with D2CO/H2CO showing a pattern that is qualitatively consistent with that of HDCO/H2CO, within current uncertainties. Our initial results show that H2CO may display a similar D$_{frac}$ pattern as that of CH3OH in massive young stellar objects. This finding suggests that solid state reactions dominate its formation.

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.

The Infrared Dark Cloud (IRDC) G028.23-00.19 hosts a massive (1,500 Msun), cold (12 K), and 3.6-70 um IR dark clump (MM1) that has the potential to form high-mass stars. We observed this prestellar clump candidate with the SMA (~3.5" resolution) and JVLA (~2.1" resolution) in order to characterize the early stages of high-mass star formation and to constrain theoretical models. Dust emission at 1.3 mm wavelength reveals 5 cores with masses <15 Msun. None of the cores currently have the mass reservoir to form a high-mass star in the prestellar phase. If the MM1 clump will ultimately form high-mass stars, its embedded cores must gather a significant amount of additional mass over time. No molecular outflows are detected in the CO (2-1) and SiO (5-4) transitions, suggesting that the SMA cores are starless. By using the NH3 (1,1) line, the velocity dispersion of the gas is determined to be transonic or mildly supersonic (DeltaV_nt}/DeltaV_th ~1.1-1.8). The cores are not highly supersonic as some theories of high-mass star formation predict. The embedded cores are 4 to 7 times more massive than the clump thermal Jeans mass and the most massive core (SMA1) is 9 times less massive than the clump turbulent Jeans mass. These values indicate that neither thermal pressure nor turbulent pressure dominates the fragmentation of MM1. The low virial parameters of the cores (0.1-0.5) suggest that they are not in virial equilibrium, unless strong magnetic fields of ~1-2 mG are present. We discuss high-mass star formation scenarios in a context based on IRDC G028.23-00.19, a study case believed to represent the initial fragmentation of molecular clouds that will form high-mass stars.

We have modified the iterative procedure introduced by Lin et al. (2016), to systematically combine the submm images taken from ground based (e.g., CSO, JCMT, APEX) and space (e.g., Herschel, Planck) telescopes. We applied the updated procedure to observations of three well studied Infrared Dark Clouds (IRDCs): G11.11-0.12, G14.225-0.506 and G28.34+0.06, and then performed single-component, modified black-body fits to derive $\sim$10$"$ resolution dust temperature and column density maps. The derived column density maps show that these three IRDCs exhibit complex filamentary structures embedding with rich clumps/cores. We compared the column density probability distribution functions (N-PDFs) and two-point correlation (2PT) functions of the column density field between these IRDCs with several OB cluster-forming regions. Based on the observed correlation and measurements, and complementary hydrodynamical simulations for a 10$^{4}$ $\rm M_{\odot}$ molecular cloud, we hypothesize that cloud evolution can be better characterized by the evolution of the (column) density distribution function and the relative power of dense structures as a function of spatial scales, rather than merely based on the presence of star-forming activity. Based on the small analyzed sample, we propose four evolutionary stages, namely: {\it cloud integration, stellar assembly, cloud pre-dispersal and dispersed-cloud.} The initial {\it cloud integration} stage and the final {\it dispersed cloud} stage may be distinguished from the two intermediate stages by a steeper than $-$4 power-law index of the N-PDF. The {\it cloud integration} stage and the subsequent {\it stellar assembly} stage are further distinguished from each other by the larger luminosity-to-mass ratio ($>$40 $\rm L_{\odot}/M_{\odot}$) of the latter.

Muon pairs can be produced in the annihilation of ultra-high energy (UHE, $E \gtrsim 10^{18} \,\mathrm{eV}$) photons with low energy cosmic background radiation in the intergalactic space, giving birth to neutrinos. Although the branching ratio of muon pair production is low, products of other channels, which are mainly electron/positron pairs, will probably transfer most of their energies into the new generated UHE photon in the subsequent interaction with the cosmic background radiation via Compton scattering in deep Klein-Nishina regime. The regeneration of these new UHE photons then provides a second chance to produce the muon pairs, enhancing the neutrino flux. We investigate the neutrino production in the propagation of UHE photons in the intergalactic space at different redshifts, considering various competing processes such as pair production, double pair production for UHE photons, and triplet production and synchrotron radiation for UHE electrons. Regarding the least energetic outgoing particles as energy loss, we obtain the effective penetration length of the leading particle, as well as energy loss rate including the neutrino emission rate in the cascade process. Finally, we find that an extra component of UHE neutrinos will arise from the propagation of UHE cosmic rays due to the generated UHE photons and electron/positrons. However, the flux of this component is quite small, with a flux of at most $10\%$ of that of the conventional cosmogenic neutrino at a few EeV, in the absence of a strong intergalactic magnetic field and a strong cosmic radio background. The precise contribution of extra component depends on several factors, e.g., cosmic radio background, intergalactic magnetic field, the spectrum of proton, which will be discussed in this work.

The origin and life-cycle of molecular clouds are still poorly constrained, despite their importance for understanding the evolution of the interstellar medium. We have carried out a systematic, homogeneous, spectroscopic survey of the inner Galactic plane, in order to complement the many continuum Galactic surveys available with crucial distance and gas-kinematic information. Our aim is to combine this data set with recent infrared to sub-millimetre surveys at similar angular resolutions. The SEDIGISM survey covers 78 deg^2 of the inner Galaxy (-60 deg < l < +18 deg, |b| < 0.5 deg) in the J=2-1 rotational transition of 13CO. This isotopologue of CO is less abundant than 12CO by factors up to 100. Therefore, its emission has low to moderate optical depths, and higher critical density, making it an ideal tracer of the cold, dense interstellar medium. The data have been observed with the SHFI single-pixel instrument at APEX. The observational setup covers the 13CO(2-1) and C18O(2-1) lines, plus several transitions from other molecules. The observations have been completed. Data reduction is in progress, and the final data products will be made available in the near future. Here we give a detailed description of the survey and the dedicated data reduction pipeline. Preliminary results based on a science demonstration field covering -20 deg < l < -18.5 deg are presented. Analysis of the 13CO(2-1) data in this field reveals compact clumps, diffuse clouds, and filamentary structures at a range of heliocentric distances. By combining our data with data in the (1-0) transition of CO isotopologues from the ThrUMMS survey, we are able to compute a 3D realization of the excitation temperature and optical depth in the interstellar medium. Ultimately, this survey will provide a detailed, global view of the inner Galactic interstellar medium at an unprecedented angular resolution of ~30".

We investigate the primordial gravitational waves (PGWs) in the general scenario where the inflation is preceded by a pre-inflationary stage with the effective equation of state $w$. Comparing with the results in the usual inflationary models, the power spectrum of PGWs is modified in two aspects: One is the mixture of the perturbation modes caused by he presence of the pre-inflationary period, and the other is the thermal initial state formed at the Planck era of the early Universe. By investigating the observational imprints of these modifications on the B-mode polarization of cosmic microwave background (CMB) radiation, we obtain the constraints on the conformal temperature of the thermal gravitational-wave background $T<5.01\times 10^{-4}$Mpc$^{-1}$ and a tensor-to-scalar ratio $r<0.084$ ($95\%$ confident level), which follows the bounds on total number of e-folds $N>63.5$ for the model with $w=1/3$, and $N>65.7$ for that with $w=1$. By taking into account various noises and the foreground radiations, we forecast the detection possibility of the thermal gravitational-wave background by the future CMBPol mission, and find that if $r>0.01$, the detection is possible as long as $T>1.5\times 10^{-4}$Mpc$^{-1}$. However, the effect of different $w$ is quite small, and it seems impossible to determine its value from the potential observations of CMBPol mission.

We use Planck data released in 2015 to constrain the abundance of primordial black holes (PBHs) in dark matter in two different reionization models (one is the instantaneous reionization and the other is the asymmetric reionization), and significantly improve the upper limits on the abundance of PBHs from WMAP 3-year data by around two orders of magnitude. Furthermore, these new limits imply that the event rates of mergers of PBH binaries (Gpc$^{-3}$ yr$^{-1}$) are less than $0.002$ for $M_\text{pbh}=30M_\odot$, $5$ for $M_\text{pbh}=10M_\odot$ and $2000$ for $M_\text{pbh}=2M_\odot$ at $95\%$ confidence level (C.L.), and thus GW150914 seems very unlikely produced by the merger of a PBH binary.

One of the main goals of modern cosmology is to search for primordial gravitational waves by looking on their imprints in the B-type polarization in the cosmic microwave background radiation. However, this signal is contaminated by various sources, including cosmic weak lensing, foreground radiations, instrumental noises, as well as the E-to-B leakage caused by the partial sky surveys, which should be well understood to avoid the misinterpretation of the observed data. In this paper, we adopt the E/B decomposition method suggested by Smith in 2006, and study the imprints of E-to-B leakage residuals in the constructed B-type polarization maps, $\mathcal{B}(\hat{n})$, by employing various statistical tools. We find that the effects of E-to-B leakage are negligible for the $\mathcal{B}$-mode power spectrum, as well as the skewness and kurtosis analyses of $\mathcal{B}$-maps. However, if employing the morphological statistical tools, including Minkowski functionals and/or Betti numbers, we find the effect of leakage can be detected at very high confidence level, which shows that in the morphological analysis, the leakage can play a significant role as a contaminant for measuring the primordial B-mode signal and must be taken into account for a correct explanation of the data.

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 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.

eXTP is a science mission designed to study the state of matter under extreme conditions of density, gravity and magnetism. Primary targets include isolated and binary neutron stars, strong magnetic field systems like magnetars, and stellar-mass and supermassive black holes. The mission carries a unique and unprecedented suite of state-of-the-art scientific instruments enabling for the first time ever the simultaneous spectral-timing-polarimetry studies of cosmic sources in the energy range from 0.5-30 keV (and beyond). Key elements of the payload are: the Spectroscopic Focusing Array (SFA) - a set of 11 X-ray optics for a total effective area of about 0.9 m^2 and 0.6 m^2 at 2 keV and 6 keV respectively, equipped with Silicon Drift Detectors offering <180 eV spectral resolution; the Large Area Detector (LAD) - a deployable set of 640 Silicon Drift Detectors, for a total effective area of about 3.4 m^2, between 6 and 10 keV, and spectral resolution <250 eV; the Polarimetry Focusing Array (PFA) - a set of 2 X-ray telescope, for a total effective area of 250 cm^2 at 2 keV, equipped with imaging gas pixel photoelectric polarimeters; the Wide Field Monitor (WFM) - a set of 3 coded mask wide field units, equipped with position-sensitive Silicon Drift Detectors, each covering a 90 degrees x 90 degrees FoV. The eXTP international consortium includes mostly major institutions of the Chinese Academy of Sciences and Universities in China, as well as major institutions in several European countries and the United States. The predecessor of eXTP, the XTP mission concept, has been selected and funded as one of the so-called background missions in the Strategic Priority Space Science Program of the Chinese Academy of Sciences since 2011. The strong European participation has significantly enhanced the scientific capabilities of eXTP. The planned launch date of the mission is earlier than 2025.

NGC 5548 is the best-observed reverberation-mapped active galactic nucleus with long-term, intensive monitoring. Here we report results from a new observational campaign between January and July, 2015. We measure the centroid time lag of the broad H$\beta$ emission line with respect to the 5100 \AA continuum and obtain $\tau_{\rm cent} = 7.20^{+1.33}_{-0.35}$ days in the rest frame. This yields a black hole mass of $M_{\bullet}=8.71^{+3.21}_{-2.61} $x$ 10^{7}M_{\odot}$ using a broad H$\beta$ line dispersion of $3124\pm302$ km s$^{-1}$ and a virial factor of $f_{_{\rm BLR}}=6.3\pm1.5$ for the broad-line region (BLR), consistent with the mass measurements from previous H$\beta$ campaigns. The high-quality data allow us to construct a velocity-binned delay map for the broad H$\beta$ line, which shows a symmetric response pattern around the line center, a plausible kinematic signature of virialized motion of the BLR. Combining all the available measurements of H$\beta$ time lags and the associated mean 5100 {\AA} luminosities over 18 campaigns between 1989 and 2015, we find that the H$\beta$ BLR size varies with the mean optical luminosity, but, interestingly, with a possible delay of $2.35_{-1.25}^{+3.47}$ yrs. This delay coincides with the typical BLR dynamical timescale of NGC 5548, indicating that the BLR undergoes dynamical changes, possibly driven by radiation pressure.

Large-scale gaseous filaments with length up to the order of 100 pc are on the upper end of the filamentary hierarchy of the Galactic interstellar medium. Their association with respect to the Galactic structure and their role in Galactic star formation are of great interest from both observational and theoretical point of view. Previous "by-eye" searches, combined together, have started to uncover the Galactic distribution of large filaments, yet inherent bias and small sample size limit conclusive statistical results to be drawn. Here, we present (1) a new, automated method to identify large-scale velocity-coherent dense filaments, and (2) the first statistics and the Galactic distribution of these filaments. We use a customized minimum spanning tree algorithm to identify filaments by connecting voxels in the position-position-velocity space, using the Bolocam Galactic Plane Survey spectroscopic catalog. In the range of $7.^{\circ}5 \le l \le 194^{\circ}$, we have identified 54 large-scale filaments and derived mass ($\sim 10^3 - 10^5 \, M_\odot$), length (10-276 pc), linear mass density (54-8625 $M_\odot \, \rm{pc}^{-1}$), aspect ratio, linearity, velocity gradient, temperature, fragmentation, Galactic location and orientation angle. The filaments concentrate along major spiral arms. They are widely distributed across the Galactic disk, with 50% located within $\pm$20 pc from the Galactic mid-plane and 27% run in the center of spiral arms (aka "bones"). An order of 1% of the molecular ISM is confined in large filaments. Massive star formation is more favorable in large filaments compared to elsewhere. This is the first comprehensive catalog of large filaments useful for a quantitative comparison with spiral structures and numerical simulations.

We used the newly commissioned 50 cm Binocular Network (50BiN) telescope at Qinghai Station of Purple Mountain Observatory (Chinese Academy of Sciences) to observe the old open cluster NGC 188 in V and R as part of a search for variable objects. Our time-series data span a total of 36 days. Radial velocity and proper-motion selection resulted in a sample of 532 genuine cluster members. Isochrone fitting was applied to the cleaned cluster sequence, yielding a distance modulus of (m - M)0V = 11.35 \pm 0.10 mag and a total foreground reddening of E(V - R) = 0.062 \pm 0.002 mag. Light-curve solutions were obtained for eight W Ursae Majoris eclipsing-binary systems (W UMas) and their orbital parameters were estimated. Using the latter parameters, we estimate a distance to the W UMas which is independent of the host cluster's physical properties. Based on combined fits to six of the W UMas (EP Cep, EQ Cep, ES Cep, V369 Cep, and--for the first time--V370 Cep and V782 Cep), we obtain an average distance modulus of (m - M)0V = 11.31 \pm 0.08 mag, which is comparable with that resulting from our isochrone fits. These six W UMas exhibit an obvious period-luminosity relation. We derive more accurate physical parameters for the W UMa systems and discuss their initial masses and ages. The former show that these W UMa systems have likely undergone angular-momentum evolution within a convective envelope (W-type evolution). The ages of the W UMa systems agree well with the cluster's age.

The detection of the magnetic type $B$-mode polarization is the main goal of future cosmic microwave background (CMB) experiments. In the standard model, the $B$-mode map is a strongly non-gaussian field due to the lensed component. Besides the two-point correlation function, the other statistics are also very important to dig the information of the polarization map. In this paper, we employ the Minkowski functionals to study the morphological properties of the lensed $B$-mode maps. We find that the deviations from Gaussianity are very significant for both full and partial-sky surveys. As an application of the analysis, we investigate the morphological imprints of the foreground residuals in the $B$-mode map. We find that even for very tiny foreground residuals, the effects on the map can be detected by the Minkowski functional analysis. Therefore, it provides a complementary way to investigate the foreground contaminations in the CMB studies.

We try to reconcile the tension between the local 2.4\% determination of Hubble constant and its global determination by $\textit{Planck}$ CMB data and BAO data through modeling the dark energy variously. We find that the chi-square is significantly reduced by $\Delta\chi^2_\text{all}=-6.76$ in the redshift-binned dark energy model where the $68\%$ limits of the equation of state of dark energy read $w(0\leq z\leq 0.1)=-1.958_{-0.508}^{+0.509}$, $w(0.1< z\leq 1.5)=-1.006_{-0.082}^{+0.092}$, and here $w(z>1.5)$ is fixed to $-1$.

We present observations towards a high-mass ($\rm >40\,M_{\odot}$), low luminosity ($\rm <10\,L_{\odot}$) $\rm 70\,\mu$m dark molecular core G 28.34 S-A at 3.4 mm, using the IRAM 30 m telescope and the NOEMA interferometer. We report the detection of $\rm SiO$ $J=\rm 2\rightarrow1$ line emission, which is spatially resolved in this source at a linear resolution of $\sim$0.1 pc, while the 3.4 mm continuum image does not resolve any internal sub-structures. The SiO emission exhibits two W-E oriented lobes centring on the continuum peak. Corresponding to the red-shifted and blue-shifted gas with velocities up to $\rm 40\,km\,s^{-1}$ relative to the quiescent cloud, these lobes clearly indicate the presence of a strong bipolar outflow from this $\rm 70\,\mu$m dark core, a source previously considered as one of the best candidates of "starless" core. Our SiO detection is consistent with ALMA archival data of $\rm SiO$ $J=\rm 5\rightarrow4$, whose high-velocity blue-shifted gas reveals a more compact lobe spatially closer to the dust center. This outflow indicates that the central source may be in an early evolutionary stage of forming a high-mass protostar. We also find that the low-velocity components (in the range of $\rm V_{lsr}$$\rm_{-5}^{+3}\,km\,s^{-1}$) have an extended, NW-SE oriented distribution. Discussing the possible accretion scenarios of the outflow-powering young stellar object, we argue that the molecular line emission and the molecular outflows may provide a better indication of the accretion history when forming young stellar object, than that from a snapshot observations of the present bolometric luminosity. This is particularly significant for the cases of episodic accretion, which may occur during the collapse of the parent molecular core.

We provide the latest constraints on the power spectra of both scalar and tensor perturbations from the CMB data (including \textit{Planck}~2015, BICEP2 \& \textit{Keck Array} experiments) and the new BAO scales from SDSS-III BOSS observation. We find that the inflation model with a convex potential is not favored and both the inflation model with a monomial potential and the natural inflation model are marginally disfavored at around $95\%$ confidence level. But both the Brane inflation model and the Starobinsky inflation model fit the data quite well.

The strikingly low luminosity of Uranus (Teff ~ Teq) constitutes a long-standing challenge to our understanding of Ice Giant planets. Here we present the first Uranus structure and evolution models that are constructed to agree with both the observed low luminosity and the gravity field data. Our models make use of modern ab initio equations of state at high pressures for the icy components water, methane, and ammonia. Proceeding step by step, we confirm that adiabatic models yield cooling times that are too long, even when uncertainties in the ice:rock ratio (I:R) are taken into account. We then argue that the transition between the ice/rock-rich interior and the H/He-rich outer envelope should be stably stratified. Therefore, we introduce a simple thermal boundary and adjust it to reproduce the low luminosity. Due to this thermal boundary, the deep interior of the Uranus models are up to 2--3 warmer than adiabatic models, necessitating the presence of rocks in the deep interior with a possible I:R of $1\times$ solar. Finally, we allow for an equilibrium evolution (Teff ~ Teq) that begun prior to the present day, which would therefore no longer require the current era to be a "special time" in Uranus' evolution. In this scenario, the thermal boundary leads to more rapid cooling of the outer envelope. When Teff ~ Teq is reached, a shallow, subadiabatic zone in the atmosphere begins to develop. Its depth is adjusted to meet the luminosity constraint. This work provides a simple foundation for future Ice Giant structure and evolution models, that can be improved by properly treating the heat and particle fluxes in the diffusive zones.

In this paper, we investigate numerically the stochastic ABC model, a toy model in the theory of astrophysical kinematic dynamos, within the recently proposed supersymmetric theory of stochastics (STS). STS characterises stochastic differential equations (SDEs) by the spectrum of the stochastic evolution operator (SEO) on elements of the exterior algebra or differentials forms over the system's phase space, X. STS can thereby classify SDEs as chaotic or non-chaotic by identifying the phenomenon of stochastic chaos with the spontaneously broken topological supersymmetry that all SDEs possess. We demonstrate the following three properties of the SEO, deduced previously analytically and from physical arguments: the SEO spectra for zeroth and top degree forms never break topological supersymmetry, all SDEs possesses pseudo-time-reversal symmetry, and each de Rahm cohomology class provides one supersymmetric eigenstate. Our results also suggests that the SEO spectra for forms of complementary degrees, i.e., k and dim X -k, may be isospectral.

This paper reports results of the third-year campaign of monitoring super-Eddington accreting massive black holes (SEAMBHs) in active galactic nuclei (AGNs) between 2014-2015. Ten new targets were selected from quasar sample of Sloan Digital Sky Survey (SDSS), which are generally more luminous than the SEAMBH candidates in last two years. H$\beta$ lags ($\tau_{_{\rm H\beta}}$) in five of the 10 quasars have been successfully measured in this monitoring season. We find that the lags are generally shorter, by large factors, than those of objects with same optical luminosity, in light of the well-known $R_{_{\rm H\beta}}-L_{5100}$ relation. The five quasars have dimensionless accretion rates of $\dot{\mathscr{M}}=10-10^3$. Combining measurements of the previous SEAMBHs, we find that the reduction of H$\beta$ lags tightly depends on accretion rates, $\tau_{_{\rm H\beta}}/\tau_{_{R-L}}\propto\dot{\mathscr{M}}^{-0.42}$, where $\tau_{_{R-L}}$ is the H$\beta$ lag from the normal $R_{_{\rm H\beta}}-L_{5100}$ relation. Fitting 63 mapped AGNs, we present a new scaling relation for the broad-line region: $R_{_{\rm H\beta}}=\alpha_1\ell_{44}^{\beta_1}\,\min\left[1,\left(\dot{\mathscr{M}}/\dot{\mathscr{M}}_c\right)^{-\gamma_1}\right]$, where $\ell_{44}=L_{5100}/10^{44}\,\rm erg~s^{-1}$ is 5100 \AA\ continuum luminosity, and coefficients of $\alpha_1=(29.6_{-2.8}^{+2.7})$ lt-d, $\beta_1=0.56_{-0.03}^{+0.03}$, $\gamma_1=0.52_{-0.16}^{+0.33}$ and $\dot{\mathscr{M}}_c=11.19_{-6.22}^{+2.29}$. This relation is applicable to AGNs over a wide range of accretion rates, from $10^{-3}$ to $10^3$. Implications of this new relation are briefly discussed.

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.

Considering the mass splitting between three active neutrinos, we represent the new constraints on the sum of neutrino mass $\sum m_\nu$ by updating the anisotropic analysis of Baryon Acoustic Oscillation (BAO) scale in the CMASS and LOWZ galaxy samples from Data Release 12 of the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS DR12). Combining the BAO data of 6dFGS, MGS, LOWZ and CMASS with $\textit{Planck}$~2015 data of temperature anisotropy and polarizations of Cosmic Microwave Background (CMB), we find that the $95\%$ C.L. upper bounds on $\sum m_\nu$ refer to $\sum m_{\nu,\rm NH}<0.18$ eV for normal hierarchy (NH), $\sum m_{\nu,\rm IH}<0.20$ eV for inverted hierarchy (IH) and $\sum m_{\nu,\rm DH}<0.15$ eV for degenerate hierarchy (DH) respectively, and the normal hierarchy is slightly preferred than the inverted one ($\Delta \chi^2\equiv \chi^2_{\rm NH}-\chi^2_{\rm IH} \simeq -3.4$). In addition, the additional relativistic degrees of freedom and massive sterile neutrinos are neither favored at present.

Previous studies of the initial conditions of massive star formation have mainly targeted Infrared-Dark Clouds (IRDCs) toward the inner Galaxy. This is due to the fact that IRDCs were first detected in absorption against the bright mid-IR background, requiring a favourable location to be observed. By selection, IRDCs represent only a fraction of the Galactic clouds capable of forming massive stars and star clusters. Due to their low dust temperatures, IRDCs are bright in the far-IR and millimeter and thus, observations at these wavelengths have the potential to provide a complete sample of star-forming massive clouds across the Galaxy. Our aim is to identify the clouds at the initial conditions of massive star formation across the Galaxy and compare their physical properties as a function of their Galactic location. We have examined the physical properties of a homogeneous galactic cold core sample obtained with the Planck satellite across the Galactic Plane. With the use of Herschel Hi-GAL observations, we have characterized the internal structure of them. By using background-subtracted Herschel images, we have derived the H2 column density and dust temperature maps for 48 Planck clumps. Their basic physical parameters have been calculated and analyzed as a function of location within the Galaxy. These properties have also been compared with the empirical relation for massive star formation derived by Kauffmann & Pillai (2010). Most of the Planck clumps contain signs of star formation. About 25% of them are massive enough to form high mass stars. Planck clumps toward the Galactic center region show higher peak column densities and higher average dust temperatures than those of the clumps in the outer Galaxy. Although we only have seven clumps without associated YSOs, the Hi-GAL data show no apparent differences in the properties of Planck cold clumps with and without star formation.

As a natural consequence of cosmological hierarchical structure formation, sub-parsec supermassive black hole binaries (SMBHBs) should be common in galaxies but thus far have eluded spectroscopic identification. Based on four decades of optical spectroscopic monitoring, we report that the nucleus of NGC 5548, a nearby Seyfert galaxy long suspected to have experienced a major merger about one billion years ago, exhibits long-term variability with a period of 14 years in the optical continuum and broad Hbeta emission line. Remarkably, the double-peaked profile of Hbeta shows systematic velocity changes with a similar period. These pieces of observations plausibly indicate that a SMBHB resides in the center of NGC 5548. The complex, secular variations in the line profiles can be explained by orbital motion of a binary with equal mass and a semi-major axis of 22 light-days (corresponding to 18 milli-parsec). At a distance of 75 Mpc, NGC 5548 is one of the nearest sub-parsec SMBHB candidates that offers an ideal laboratory for gravitational wave detection.

SONG is a global ground based network of 1 meter telescopes for stellar time-domain science, an international collaboration involving many countries across the world. In order to enable a favourable duty cycle, the SONG network plans to create a homogeneous distribution of 4 nodes in each of the northern and southern hemispheres. A natural possibility was building one of the northern nodes in East Asia, preferably on the Qinghai-Tibetan Plateau. During the last decade, a great deal of effort has been invested in searching for high a quality site for ground based astronomy in China, since this has been one of the major concerns for the development of Chinese astronomy. A number of sites on the plateau have been in operation for many years, but most of them are used only for radio astronomy, as well as small optical telescopes for applied astronomy. Several potential sites for large optical instruments have been identified by the plateau site survey, but as yet none of them have been adequately quantitatively characterised. Here we present results from a detailed multi-year study of the Delingha site, which was eventually selected for the SONG-China node. We also describe the site monitoring system that will allow an isolated SONG and 50BiN node to operate safely in an automated mode.

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.

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 update the distance priors by adopting $Planck~ \textrm{TT,TE,EE}+\textrm{lowP}$ data released in 2015, and our results impose at least $30\%$ tighter constraints than those from $Planck~ \textrm{TT}+\textrm{lowP}$. Combining the distance priors with the combination of supernova Union~2.1 compilation of 580 SNe (Union~2.1) and low redshift Baryon Acoustic Oscillation (BAO) data, we constrain the cosmological parameters in the freely binned dark energy (FBDE) and FBDE$+\Omega_k$ models respectively, and find that the equations of state of dark energy in both models are consistent with $w=-1$. Furthermore, we show that the tension with the BAO data at $z=2.34$ from Ly$\alpha$ forest (Ly$\alpha$F) auto-correlation and Combined Ly$\alpha$F cannot be relaxed in the FBDE and FBDE$+\Omega_k$ models.

Fomalhaut plays an important role in the study of debris disks and small bodies in other planetary systems. The proximity and luminosity of the star make key features of its debris, like the water ice-line, accessible. Here we present ALMA cycle 1, 870 \mu m (345 GHz) observations targeted at the inner part of the Fomalhaut system with a synthesized beam of 0.45"x0.37" (~3 AU linear resolution at the distance of Fomalhaut) and a rms of 26 \mu Jy/beam. The high angular resolution and sensitivity of the ALMA data enable us to place strong constraints on the nature of the warm excess revealed by Spitzer and Herschel observations. We detect a point source at the star position with a total flux consistent with thermal emission from the stellar photosphere. No structures that are brighter than 3\sigma\ are detected in the central 15 AU x 15 AU region. Modeling the spectral energy distribution using parameters expected for a dust-producing planetesimal belt indicates a radial location in the range ~8-15 AU. This is consistent with the location where ice sublimates in Fomalhaut, i.e., an asteroid-belt analog. The 3\sigma\ upper limit for such a belt is <1.3 mJy at 870 \mu m. We also interpret the 2 and 8-13 \mu m interferometric measurements to reveal the structure in the inner 10 AU region as dust naturally connected to this proposed asteroid belt by Poynting-Robertson drag, dust sublimation, and magnetically trapped nano grains.

The anisotropies of the B-mode polarization in the cosmic microwave background radiation play a crucial role for the study of the very early Universe. However, in the real observation, the mixture of the E-mode and B-mode can be caused by the partial sky surveys, which must be separated before applied to the cosmological explanation. The separation method developed by Smith (\citealt{PhysRevD.74.083002}) has been widely adopted, where the edge of the top-hat mask should be smoothed to avoid the numerical errors. In this paper, we compare three different smoothing methods, and investigate the leakage residuals of the E-B mixture. We find that, if the less information loss is needed and the smaller region is smoothed in the analysis, the \textit{sin}- and \textit{cos}-smoothing methods are better. However, if we need a clean constructed B-mode map, the larger region around the mask edge should be smoothed. In this case, the \textit{Gaussian}-smoothing method becomes much better. In addition, we find that the leakage caused by the numerical errors in the \textit{Gaussian}-smoothing method mostly concentrates on two bands, which is quite easy to be reduced for the further E-B separations.

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 ubiquity of filamentary structure at various scales through out the Galaxy has triggered a renewed interest in their formation, evolution, and role in star formation. The largest filaments can reach up to Galactic scale as part of the spiral arm structure. However, such large scale filaments are hard to identify systematically due to limitations in identifying methodology (i.e., as extinction features). We present a new approach to directly search for the largest, coldest, and densest filaments in the Galaxy, making use of sensitive Herschel Hi-GAL data complemented by spectral line cubes. We present a sample of the 9 most prominent Herschel filaments, including 6 identified from a pilot search field plus 3 from outside the field. These filaments measure 37-99 pc long and 0.6-3.0 pc wide with masses (0.5-8.3)$\times10^4 \, M_\odot$, and beam-averaged ($28"$, or 0.4-0.7 pc) peak H$_2$ column densities of (1.7-9.3)$\times 10^{22} \, \rm{cm^{-2}}$. The bulk of the filaments are relatively cold (17-21 K), while some local clumps have a dust temperature up to 25-47 K. All the filaments are located within <~60 pc from the Galactic mid-plane. Comparing the filaments to a recent spiral arm model incorporating the latest parallax measurements, we find that 7/9 of them reside within arms, but most are close to arm edges. These filaments are comparable in length to the Galactic scale height and therefore are not simply part of a grander turbulent cascade.

We present multi-color CCD photometry of the neglected contact binary XZ Leo. Completely covered VRI band light curves and four times of minimum light were obtained. Combining the photometric and previously published radial velocity data, a revised photometric analysis was carried out for the binary system by applying the Wilson-Devinney code. With a hot spot placed on the massive primary component near the neck region of the common envelope, the light curves were satisfactorily modeled. The photometric solution combined with the radial velocity solution reveals that XZ Leo is an A-type contact binary with a degree of contact of 24($\pm1)\%$. The absolute parameters of the components were determined as M_1 = 1.74($\pm$0.06)M_\odot, M_2 = 0.61($\pm$0.02)M_\odot, R_1 = 1.69($\pm$0.01)R_\odot, R_2 = 1.07($\pm0.01$)R_\odot, L_1 = 6.73($\pm0.08$) L_\odot, L_2 = 2.40($\pm$0.04)L_\odot. Based on all the available data, the long-term orbital period behavior of the system was investigated. It indicates that the binary system was undergoing continuous orbital period increase in the past three decades with a rate of dP/dt = + 6.12 \times {10^{-8}} days yr^{-1}, which suggests a probable mass transfer from the secondary to the primary component at a rate of dM/dt= 3.92\times 10^{-8} M_\odot yr^{-1}. The binary system is expected to evolve into the broken-contact stage in $1.56 \times 10^6$ years. This could be evidence supporting the Thermal Relaxation Oscillation theory.

A major goal of the Atacama Large Millimeter/submillimeter Array (ALMA) is to make accurate images with resolutions of tens of milliarcseconds, which at submillimeter (submm) wavelengths requires baselines up to ~15 km. To develop and test this capability, a Long Baseline Campaign (LBC) was carried out from September to late November 2014, culminating in end-to-end observations, calibrations, and imaging of selected Science Verification (SV) targets. This paper presents an overview of the campaign and its main results, including an investigation of the short-term coherence properties and systematic phase errors over the long baselines at the ALMA site, a summary of the SV targets and observations, and recommendations for science observing strategies at long baselines. Deep ALMA images of the quasar 3C138 at 97 and 241 GHz are also compared to VLA 43 GHz results, demonstrating an agreement at a level of a few percent. As a result of the extensive program of LBC testing, the highly successful SV imaging at long baselines achieved angular resolutions as fine as 19 mas at ~350 GHz. Observing with ALMA on baselines of up to 15 km is now possible, and opens up new parameter space for submm astronomy.

A major goal of the Atacama Large Millimeter/submillimeter Array (ALMA) is to make accurate images with resolutions of tens of milliarcseconds, which at submillimeter (submm) wavelengths requires baselines up to ~15 km. To develop and test this capability, a Long Baseline Campaign (LBC) was carried out from September to late November 2014, culminating in end-to-end observations, calibrations, and imaging of selected Science Verification (SV) targets. This paper presents an overview of the campaign and its main results, including an investigation of the short-term coherence properties and systematic phase errors over the long baselines at the ALMA site, a summary of the SV targets and observations, and recommendations for science observing strategies at long baselines. Deep ALMA images of the quasar 3C138 at 97 and 241 GHz are also compared to VLA 43 GHz results, demonstrating an agreement at a level of a few percent. As a result of the extensive program of LBC testing, the highly successful SV imaging at long baselines achieved angular resolutions as fine as 19 mas at ~350 GHz. Observing with ALMA on baselines of up to 15 km is now possible, and opens up new parameter space for submm astronomy.

To study the fragmentation and gravitational collapse of dense cores in infrared dark clouds (IRDCs), we have obtained submillimeter continuum and spectral line data as well as multiple inversion transitions of NH3 and H2O maser data of four massive clumps in an IRDC G28.53-0.25. Combining single dish and interferometer NH3 data, we derive the rotation temperature of G28.53. We identity 12 dense cores at 0.1 pc scale based on submillimeter continuum, and obtain their physical properties using NH3 and continuum data. By comparing the Jeans masses of cores with the core masses, we find that turbulent pressure is important in supporting the gas when 1 pc scale clumps fragment into 0.1 pc scale cores. All cores have a virial parameter smaller than 1 assuming a inverse squared radial density profile, suggesting they are gravitationally bound, and the three most promising star forming cores have a virial parameter smaller than 1 even taking magnetic field into account. We also associate the cores with star formation activities revealed by outflows, masers, or infrared sources. Unlike what previous studies suggested, MM1 turns out to harbor a few star forming cores and is likely a progenitor of high-mass star cluster. MM5 is intermediate while MM7/8 are quiescent in terms of star formation, but they also harbor gravitationally bound dense cores and have the potential of forming stars as in MM1.

Sufficiently massive clumps of molecular gas collapse under self-gravity and fragment to spawn a cluster of stars that have a range of masses. We investigate observationally the early stages of formation of a stellar cluster in a massive filamentary infrared dark cloud, G28.34+0.06 P1, in the 1.3mm continuum and spectral line emission using the ALMA. Sensitive continuum data reveal further fragmentation in five dusty cores at a resolution of several 10^3 AU. Spectral line emission from C18O, CH3OH, 13CS, H2CO and N2D+ are detected for the first time toward these dense cores. We found that three cores are chemically more evolved as compared with the other two; interestingly though, all of them are associated with collimated outflows as suggested by evidence from the CO, SiO, CH3OH, H2CO and SO emissions. The parsec-scale kinematics in NH3 exhibit velocity gradients along the filament, consistent with accretion flows toward the clumps and cores. The moderate luminosity and the chemical signatures indicate that the five cores harbor low- to intermediate-mass protostars that likely become massive ones at the end of the accretion. Despite the fact that the mass limit reached by the 1\sigma dust continuum sensitivity is 30 times lower than the thermal Jeans mass, there is a lack of a distributed low-mass protostellar population in the clump. Our observations indicate that in a protocluster, low-mass stars form at a later stage after the birth of more massive protostars.

We report the first spectroscopic detection of ethyl cyanide (C$_2$H$_5$CN) in Titan's atmosphere, obtained using spectrally and spatially resolved observations of multiple emission lines with the Atacama Large Millimeter/submillimeter array (ALMA). The presence of C$_2$H$_5$CN in Titan's ionosphere was previously inferred from Cassini ion mass spectrometry measurements of C$_2$H$_5$CNH$^+$. Here we report the detection of 27 rotational lines from C$_2$H$_5$CN (in 19 separate emission features detected at $>3\sigma$ confidence), in the frequency range 222-241 GHz. Simultaneous detections of multiple emission lines from HC$_3$N, CH$_3$CN and CH$_3$CCH were also obtained. In contrast to HC$_3$N, CH$_3$CN and CH$_3$CCH, which peak in Titan's northern (spring) hemisphere, the emission from C$_2$H$_5$CN is found to be concentrated in the southern (autumn) hemisphere, suggesting a distinctly different chemistry for this species, consistent with a relatively short chemical lifetime for C$_2$H$_5$CN. Radiative transfer models show that most of the C$_2$H$_5$CN is concentrated at altitudes 300-600 km, suggesting production predominantly in the mesosphere and above. Vertical column densities are found to be in the range (2-5)$\times10^{14}$ cm$^{-2}$.

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.

Supersonic radiation diffusion approximation is a useful way to study the radiation transportation. Considering the bent Marshak wave theory in 2-dimensions, and an invariable source temperature, we get the supersonic radiation diffusion conditions which are about the Mach number $M>8(1+\sqrt{\ep})/3$, and the optical depth $\tau>1$. A large Mach number requires a high temperature, while a large optical depth requires a low temperature. Only when the source temperature is in a proper region these conditions can be satisfied. Assuming the material opacity and the specific internal energy depend on the temperature and the density as a form of power law, for a given density, these conditions correspond to a region about source temperature and the length of the sample. This supersonic diffusion region involves both lower and upper limit of source temperature, while that in 1-dimension only gives a lower limit. Taking $\rm SiO_2$ and the Au for example, we show the supersonic region numerically.

We present a comprehensive photometric study of the pulsating, eclipsing binary OO Dra. Simultaneous B- and V-band photometry of the star was carried out on 14 nights. Revised orbital period and a new ephemeris were derived from the data. The first photometric solution of the binary system and the physical parameters of the component stars are determined. It reveals that OO Dra could be a detached system with the less-massive secondary component nearly filling in its Roche lobe. By subtracting the eclipsing light changes from the data, we obtained the intrinsic pulsating light curves of the hotter and massive primary component. Frequency analysis of the residuals light yields two confident pulsation modes in both B- and V-band data with the dominant frequency detected at 41.865 c/d. A brief discussion concerning the evolutionary status and the pulsation nature of the binary system is finally given.

The present study aims at characterizing the massive star forming region G35.20N, which is found associated with at least one massive outflow and contains multiple dense cores, one of them recently found associated with a Keplerian rotating disk. We used ALMA to observe the G35.20N region in the continuum and line emission at 350 GHz. The observed frequency range covers tracers of dense gas (e.g. H13CO+, C17O), molecular outflows (e.g. SiO), and hot cores (e.g. CH3CN, CH3OH). The ALMA 870 um continuum emission map reveals an elongated dust structure (0.15 pc long and 0.013 pc wide) perpendicular to the large-scale molecular outflow detected in the region, and fragmented into a number of cores with masses 1-10 Msun and sizes 1600 AU. The cores appear regularly spaced with a separation of 0.023 pc. The emission of dense gas tracers such as H13CO+ or C17O is extended and coincident with the dust elongated structure. The three strongest dust cores show emission of complex organic molecules characteristic of hot cores, with temperatures around 200 K, and relative abundances 0.2-2x10^(-8) for CH3CN and 0.6-5x10^(-6) for CH3OH. The two cores with highest mass (cores A and B) show coherent velocity fields, with gradients almost aligned with the dust elongated structure. Those velocity gradients are consistent with Keplerian disks rotating about central masses of 4-18 Msun. Perpendicular to the velocity gradients we have identified a large-scale precessing jet/outflow associated with core B, and hints of an east-west jet/outflow associated with core A. The elongated dust structure in G35.20N is fragmented into a number of dense cores that may form massive stars. Based on the velocity field of the dense gas, the orientation of the magnetic field, and the regularly spaced fragmentation, we interpret this elongated structure as the densest part of a 1D filament fragmenting and forming massive stars.

This paper presents the richness of submillimeter spectral features in the high-mass star forming region AFGL 2591. As part of the CHESS (Chemical Herschel Survey of Star Forming Regions) Key Programme, AFGL 2591 was observed by the Herschel/HIFI instrument. The spectral survey covered a frequency range from 480 up to 1240 GHz as well as single lines from 1267 to 1901 GHz (i.e. CO, HCl, NH3, OH and [CII]). Rotational and population diagram methods were used to calculate column densities, excitation temperatures and the emission extents of the observed molecules associated with AFGL 2591. The analysis was supplemented with several lines from ground-based JCMT spectra. From the HIFI spectral survey analysis a total of 32 species were identified (including isotopologues). In spite of the fact that lines are mostly quite week, 268 emission and 16 absorption lines were found (excluding blends). Molecular column densities range from 6e11 to 1e19 cm-2 and excitation temperatures range from 19 to 175 K. One can distinguish cold (e.g. HCN, H2S, NH3 with temperatures below 70 K) and warm species (e.g. CH3OH, SO2) in the protostellar envelope.

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

We present Submillimeter Array (SMA) $\lambda =$ 0.88 and 1.3 mm broad band observations, and the Jansky Very Large Array (VLA) observations in $\rm{NH_3}$ $(J,K) = (1,1)$ up to $(5,5)$, $\rm{H_2O}$ and $\rm{CH_3OH}$ maser lines toward the two most massive molecular clumps in infrared dark cloud (IRDC) G11.11-0.12. Sensitive high-resolution images reveal hierarchical fragmentation in dense molecular gas from the $\sim 1$ pc clump scale down to $\sim 0.01$ pc condensation scale. At each scale, the mass of the fragments is orders of magnitude larger than the Jeans mass. This is common to all four IRDC clumps we studied, suggesting that turbulence plays an important role in the early stages of clustered star formation. Masers, shock heated $\rm{NH_3}$ gas, and outflows indicate intense ongoing star formation in some cores while no such signatures are found in others. Furthermore, chemical differentiation may reflect the difference in evolutionary stages among these star formation seeds. We find $\rm{NH_3}$ ortho/para ratios of $1.1\pm0.4$, $2.0\pm0.4$, and $3.0\pm0.7$ associated with three outflows, and the ratio tends to increase along the outflows downstream. Our combined SMA and VLA observations of several IRDC clumps present the most in depth view so far of the early stages prior to the hot core phase, revealing snapshots of physical and chemical properties at various stages along an apparent evolutionary sequence.

[Context] Two competing models describe the formation of massive stars in objects like the Orion Trapezium. In the turbulent core accretion model, the resulting stellar masses are directly related to the mass distribution of the cloud condensations. In the competitive accretion model, the gravitational potential of the protocluster captures gas from the surrounding cloud for which the individual cluster members compete. [Aims] With high resolution submillimeter observations of the structure, kinematics, and chemistry of the proto-Trapezium cluster W3 IRS5, we aim to determine which mode of star formation dominates. [Methods] We present 354 GHz Submillimeter Array observations at resolutions of 1"-3" (1800-5400 AU) of W3 IRS5. ...... [Results] The observations show five emission peaks (SMM1-5). SMM1 and SMM2 contain massive embedded stars (~20 Msun); SMM3-5 are starless or contain low-mass stars (<8 Msun). The inferred densities are high, >= 10^7 cm^-3, but the core masses are small, 0.2-0.6 Msun. The detected molecular emission reveals four different chemical zones. ...... [Conclusions] The proto-Trapezium cluster W3 IRS5 is an ideal test case to discriminate between models of massive star formation. Either the massive stars accrete locally from their local cores; in this case the small core masses imply that W3 IRS5 is at the very end stages (1000 yr) of infall and accretion, or the stars are accreting from the global collapse of a massive, cluster forming core. We find that the observed masses, densities and line widths observed toward W3 IRS 5 and the surrounding cluster forming core are consistent with the competitive accretion of gas at rates of Macc~10^-4 Msun yr^-1 by the massive young forming stars. ......

Recently, the GeV radiation during the X-ray flare activity in GRB 100728A was detected by Feimi/LAT. Here we study the dynamics and emission properties of collision between two homogeneous shells based on the late internal shock model. The GeV photons can be produced from X-ray flare photons up-scattered by relativistic electrons accelerated by forward-reverse shocks, where involved radiative processes include synchrotron self-Compton and crossing inverse-Compton scattering. By analytical and numerical calculations, the observed spectral properties in GRB 100728A can be well explained.

The LAMOST spectroscopic survey for Galactic structure and evolution has been in operation since October 2012, following a one-year pilot survey. The pilot survey produced a data release containing over 600,000 stellar spectra. By cross-checking with a large time series photometric database of RR Lyrae stars in high Galactic latitude regions, we found a total number of 157 stars that have been observed with LAMOST. In this sample, we successfully captured three RR Lyrae stars in the fast expansion phase, all of them showing significant hypersonic shock wave features in the spectra. By fitting the H{\alpha} line shape, we determine that the emission features seen within the broader H{\alpha} absorption line suggest hypersonic relative motions in the atmospheres of these three objects. With a further LAMOST survey of millions of stars, we will capture a large sample of RR Lyrae stars in their hypersonic expansion phase, and therefore provide a large database for the study of the internal structure and pulsation mechanism of RR Lyrae stars.

We report on ALMA observations of continuum and molecular line emission with 0.4" resolution towards the high-mass star forming region G35.20-0.74 N. Two dense cores are detected in typical hot-core tracers, such as CH3CN, which reveal velocity gradients. In one of these cores, the velocity field can be fitted with an almost edge-on Keplerian disk rotating about a central mass of 18 Msun. This finding is consistent with the results of a recent study of the CO first overtone bandhead emission at 2.3mum towards G35.20-0.74 N. The disk radius and mass are >2500 au and 3 Msun. To reconcile the observed bolometric luminosity (3x10^4 Lsun) with the estimated stellar mass of 18 Msun, we propose that the latter is the total mass of a binary system.

[Context] Recent detections of disks around young high-mass stars support the idea of massive star formation through accretion rather than coalescence, but the detailed kinematics in the equatorial region of the disk candidates is not well known, which limits our understanding of the accretion process. [Aims] This paper explores the kinematics of the gas around a young massive star with millimeter-wave interferometry to improve our understanding of the formation of massive stars though accretion. [Methods] We use Plateau de Bure interferometric images to probe the environment of the nearby (~1 kpc) and luminous (~20000 Lsun) high-mass (10-16 Msun) young star AFGL 2591-VLA3 in continuum and in lines of HDO, H218O and SO2 in the 115 and 230 GHz bands. Radiative transfer calculations are employed to investigate the kinematics of the source. [Results] At ~0.5" (500 AU) resolution, the line images clearly resolve the velocity field of the central compact source (diameter of ~ 800 AU) and show linear velocity gradients in the northeast-southwest direction. Judging from the disk-outflow geometry, the observed velocity gradient results from rotation and radial expansion in the equatorial region of VLA3. Radiative transfer calculations suggest that the velocity field is consistent with sub-Keplerian rotation plus Hubble-law like expansion. The line profiles of the observed molecules suggest a layered structure, with HDO emission arising from the disk mid-plane, H218O from the warm mid-layer, and SO2 from the upper disk. [Conclusions] We propose AFGL 2591-VLA3 as a new massive disk candidate, with peculiar kinematics. The rotation of this disk is sub-Keplerian, probably due to magnetic braking, while the stellar wind may be responsible for the expansion of the disk. The expansion motion [...]

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.

The dense molecular clump P1 in the infrared dark cloud (IRDC) complex G28.34+0.06 harbors a massive protostellar cluster at its extreme youth. Our previous Submillimeter Array (SMA) observations revealed several jet-like CO outflows emanating from the protostars, indicative of intense accretion and potential interaction with ambient natal materials. Here we present the Expanded Very Large Array (EVLA) spectral line observations toward P1 in the NH3 (J,K) = (1,1), (2,2), (3,3) lines, as well as H2O and class I CH3OH masers. Multiple NH3 transitions reveal the heated gas widely spread in the 1 pc clump. The temperature distribution is highly structured; the heated gas is offset from the protostars, and morphologically matches the outflows very well. Hot spots of spatially compact, spectrally broad NH3 (3,3) emission are also found coincident with the outflows. A weak NH3 (3,3) maser is discovered at the interface between an outflow jet and the ambient gas. These findings suggest that protostellar heating may not be effective in suppressing fragmentation during the formation of massive cores.

We observe the 1.2 mm continuum emission around the OB cluster forming region G10.6-0.4, using the IRAM 30m telescope MAMBO-2 bolometer array and the Submillimeter array. Comparison of the Spitzer 24 $\mu$m and 8 $\mu$m images with our 1.2 mm continuum maps reveals the ionization front of an HII region, the photon-dominated layer, and several 5 pc scale filaments following the outer edge of the photon-dominated layer. The filaments, which are resolved in the MAMBO-2 observations, show regularly spaced parsec-scale molecular clumps, embedded with a cluster of submillimeter molecular cores as shown in the SMA 0.87 mm observations. Toward the center of the G10.6-0.4 region, the combined SMA+IRAM 30m continuum image reveals several, parsec-scale protrusions. They may continue down to within 0.1 pc of the geometric center of a dense 3 pc size structure, where a 200 M$_{\odot}$ OB cluster resides. The observed filaments may facilitate mass accretion onto the central cluster--forming region in the presence of strong radiative and mechanical stellar feedbacks. Their filamentary geometry may also facilitate fragmentation. We did not detect any significant polarized emission at 0.87 mm in the inner 1 pc region with the SMA.

We present Submillimeter Array (SMA) \lambda = 0.88mm observations of an infrared dark cloud (IRDC) G28.34+0.06. Located in the quiescent southern part of the G28.34 cloud, the region of interest is a massive ($>10^3$\,\msun) molecular clump P1 with a luminosity of $\sim 10^3$ \lsun, where our previous SMA observations at 1.3mm have revealed a string of five dust cores of 22-64 \msun\ along the 1 pc IR-dark filament. The cores are well aligned at a position angle of 48 degrees and regularly spaced at an average projected separation of 0.16 pc. The new high-resolution, high-sensitivity 0.88\,mm image further resolves the five cores into ten compact condensations of 1.4-10.6 \msun, with sizes a few thousands AU. The spatial structure at clump ($\sim 1$ pc) and core ($\sim 0.1$ pc) scales indicates a hierarchical fragmentation. While the clump fragmentation is consistent with a cylindrical collapse, the observed fragment masses are much larger than the expected thermal Jeans masses. All the cores are driving CO(3-2) outflows up to 38 km/s, majority of which are bipolar, jet-like outflows. The moderate luminosity of the P1 clump sets a limit on the mass of protostars of 3-7 \msun. Because of the large reservoir of dense molecular gas in the immediate medium and ongoing accretion as evident by the jet-like outflows, we speculate that P1 will grow and eventually form a massive star cluster. This study provides a first glimpse of massive, clustered star formation that currently undergoes through an intermediate-mass stage.

Massive stars (M $\gsim 10$ \msun) form from collapse of parsec-scale molecular clumps. How molecular clumps fragment to give rise to massive stars in a cluster with a distribution of masses is unclear. We search for cold cores that may lead to future formation of massive stars in a massive ($> 10^3$ \msun), low luminosity ($4.6 \times 10^2$ \lsun) infrared dark cloud (IRDC) G030.88+00.13. The \nh3 data from VLA and GBT reveal that the extinction feature seen in the infrared consists of two distinctive clumps along the same line of sight: The C1 clump at 97 \kms-1 coincides with the extinction in the Spitzer 8 and 24 $\mu$m. Therefore, it is responsible for the majority of the IRDC. The C2 clump at 107 \kms-1 is more compact and has a peak temperature of 45 K. Compact dust cores and \h2O masers revealed in the SMA and VLA observations are mostly associated with C2, and none is within the IRDC in C1. The luminosity indicates that neither the C1 nor C2 clump has yet to form massive protostars. But C1 might be at a precluster forming stage. The simulated observations rule out 0.1pc cold cores with masses above 8 \msun\ within the IRDC. The core masses in C1 and C2, and those in high-mass protostellar objects suggest an evolutionary trend that the mass of cold cores increases over time. Based on our findings, we propose an empirical picture of massive star formation that protostellar cores and the embedded protostars undergo simultaneous mass growth during the protostellar evolution.

We search for evidence of a relation between properties of young stellar objects (YSOs) and their parent molecular clouds to understand the initial conditions of high-mass star formation. A sample of 135 sources was selected from the Infrared Astronomical Satellite (IRAS) Point Source Catalog, on the basis of their red color to enhance the possibility of discovering young sources. Using the Kolner Observatorium fur SubMillimeter Astronomie (KOSMA) 3-m telescope, a single-point survey in 13CO(2-1) was carried out for the entire sample, and 14 sources were mapped further. Archival mid-infrared (MIR) data were compared with the 13CO emissions to identify evolutionary stages of the sources. A 13CO observed sample was assembled to investigate the correlation between 13CO line width of the clouds and the luminosity of the associated YSOs. We identified 98 sources suitable for star formation analyses for which relevant parameters were calculated. We detected 18 cores from 14 mapped sources, which were identified with eight pre-UC HII regions and one UC HII region, two high-mass cores earlier than pre-UC HII phase, four possible star forming clusters, and three sourceless cores. By compiling a large (360 sources) 13CO observed sample, a good correlation was found between the 13CO line width of the clouds and the bolometric luminosity of the associated YSOs, which can be fitted as a power law: lg(dV13/km/s)=-0.023+0.135lg(Lbol/Lsolar). Results show that luminous (>10^3Lsolar) YSOs tend to be associated with both more massive and more turbulent (dV13>2km/s) molecular cloud structures.

We present observations of CH3CN (12-11) emission at a resolution of 2" toward the shell-like ultracompact HII region G5.89-0.39 with the Submillimeter Array. The integrated CH3CN emission reveals dense and hot molecular cocoon in the periphery of the HII region G5.89-0.39, with a CH3CN deficient region roughly centered at G5.89-0.39. By analyzing the CH3CN emission using population diagram analysis, we find, for the first time, a decreasing temperature structure from 150 to 40 K with the projected distance from Feldt's star, which is thought to be responsible for powering the HII region. Our results further indicate that the majority of the heating energy in the observed dense gas is supplied by the Feldt's star. From the derived CH3CN column density profile, we conclude that the dense gas is not uniformly-distributed but centrally-concentrated, with a power-law exponent of 5.5 for r < 8000 AU, and 2.0 for 8000 AU < r < 20000 AU, where r is the distance to Feldt's star. The estimated large power index of 5.5 can be attributed to an enhancement of CH3CN abundance in the close vicinity of Feldt's star.

This paper has been withdraw by the author due to adding more modeling results.

Evidence of TeV emission from GRB970417a has been previously reported using data from the Milagrito detector. Constraints on the TeV fluence and the energy spectrum are now derived using additional data from a scaler system that recorded the rate of signals from the Milagrito photomultipliers. This analysis shows that if emission from GRB970417a has been observed, it must contain photons with energies above 650 GeV. Some consequences of this observation are discussed.

A search for steady TeV point sources anywhere in the northern sky has been made with data from the Milagrito air-shower-particle detector. Over 3 x 10**9 events collected from 1997 February to 1998 May have been used in this study. No statistically significant excess above the background from the isotropic flux of cosmic rays was found for any direction of the sky with declination between -5 degrees and 71.7 degrees. Upper limits are derived for the photon flux above 1 TeV from any steady point source in the northern sky.

Milagrito, a detector sensitive to very high energy gamma rays, monitored the northern sky from February 1997 through May 1998. With a large field of view and a high duty cycle, this instrument was well suited to perform a search for TeV gamma-ray bursts (GRBs). We report on a search made for TeV counterparts to GRBs observed by BATSE. BATSE detected 54 GRBs within the field of view of Milagrito during this period. An excess of events coincident in time and space with one of these bursts, GRB 970417a, was observed by Milagrito. The excess has a chance probability of $2.8 \times 10^{-5}$ of being a fluctuation of the background. The probability for observing an excess at least this large from any of the 54 bursts is $1.5 \times 10^{-3}$. No significant correlations were detected from the other bursts.

Milagrito, a detector sensitive to gamma-rays at TeV energies, monitored the northern sky during the period February 1997 through May 1998. With a large field of view and high duty cycle, this instrument was used to perform a search for TeV counterparts to gamma-ray bursts. Within the Milagrito field of view 54 gamma-ray bursts at keV energies were observed by the Burst And Transient Satellite Experiment (BATSE) aboard the Compton Gamma-Ray Observatory. This paper describes the results of a preliminary analysis to search for TeV emission correlated with BATSE detected bursts. Milagrito detected an excess of events coincident both spatially and temporally with GRB 970417a, with chance probability $2.8 \times 10^{-5}$ within the BATSE error radius. No other significant correlations were detected. Since 54 bursts were examined the chance probability of observing an excess with this significance in any of these bursts is $1.5 \times 10^{-3}$. The statistical aspects and physical implications of this result are discussed.

Previous works have shown that water Cherenkov detectors have superior sensitivity to those of scintillation counters as applied to detecting extensive air showers (EAS). This is in large part due to their much higher sensitivity to EAS photons which are more than five times more numerous than EAS electrons. Large area water Cherenkov detectors can be constructed relatively cheaply and operated reliably. A sparse detector array has been designed which uses these types of detectors to substantially increase the area over which the Milagro Gamma Ray Observatory collects EAS information. Improvements to the Milagro detector's performance characteristics and sensitivity derived from this array and preliminary results from a prototype array currently installed near the Milagro detector will be presented.

The Milagro detector is an air shower array which uses the water Cherenkov technique and is capable of continuously monitoring the sky at energies near 1 TeV. The detector consists of 20000 metric tons of pure water instrumented with 723 photo-multiplier tubes (PMTs). The PMTs are arranged in a two-layer structure on a lattice of 3 m spacing covering 5000 $m^2$ area. The direction of the shower is determined from the relative timing of the PMT signals, necessitating a common time reference and amplitude slewing corrections to improve the time resolution. The calibration system to provide these consists of a pulsed laser driving 30 diffusing light sources deployed in the pond to allow cross-calibration of the PMTs. The system is capable of calibrating times and the pulse-heights from the PMTs using the time-over-threshold technique. The absolute energy scale is provided using single muons passing through the detector. The description of the calibration system of the Milagro detector and its prototype Milagrito will be presented.

The Milagro Gamma Ray Observatory is the world's first large-area water Cherenkov detector capable of continuously monitoring the sky at TeV energies. Located in northern New Mexico, Milagro will perform an all sky survey of the Northern Hemisphere at energies between ~250 GeV and 50 TeV. With a high duty cycle, large detector area (~5000 square meters), and a wide field-of-view (~1 sr), Milagro is uniquely capable of searching for transient and DC sources of high-energy gamma-ray emission. Milagro has been operating since February, 1999. The current status of the Milagro Observatory and initial results will be discussed.

Observing gamma ray bursts (GRBs) in the TeV energy range can be extremely valuable in providing insight to GRB radiation mechanisms and in constraining source distances. The Milagrito detector was an air shower array which used the water Cherenkov technique to search for TeV sources. Data from this detector was analyzed to look for a TeV component of GRBs coincident with low energy -rays detected by the BATSE instrument on the Compton Gamma Ray Observatory. A sample of 54 BATSE GRBs which were in the field of view of the Milagrito detector during its lifetime (February 1997 to May 1998) was used.

Solar Energetic Particles from the 6 November 1997 solar flare/CME(coronal mass ejection) with energies exceeding 10 GeV have been detected by Milagrito, a prototype of the Milagro Gamma Ray Observatory. While particle acceleration beyond 1 GeV at the Sun is well established, few data exist for protons or ions beyond 10 GeV. The Milagro observatory, a ground based water Cherenkov detector designed for observing very high energy gamma ray sources, can also be used to study the Sun. Milagrito, which operated for approximately one year in 1997/98, was sensitive to solar proton and neutron fluxes above ~5- 10 GeV. Milagrito operated in a scaler mode, which was primarily sensitive to muons, low energy photons, and electrons, and the detector operated in a mode sensitive to showers and high zenith angle muons. In its scaler mode, Milagrito registered a rate increase coincident with the 6 November 1997 ground level event observed by Climax and other neutron monitors. A preliminary analysis suggests the presence of >10 GeV particles.

Milagrito, a prototype for the Milagro detector, operated for 15 months in 1997-8 and collected 8.9 billion events. It was the first extensive air shower (EAS) array sensitive to showers intiated by primaries with energy below 1 TeV. The shadows of the sun and moon observed with cosmic rays can be used to study systematic pointing shifts and measure the angular resolution of EAS arrays. Below a few TeV, the paths of cosmic rays coming toward the earth are bent by the helio- and geo-magnetic fields. This is expected to distort and displace the shadows of the sun and the moon. The moon shadow, offset from the nominal (undeflected) position, has been observed with high statistical significance in Milagrito. This can be used to establish energy calibrations, as well as to search for the anti-matter content of the VHE cosmic ray flux. The shadow of the sun has also been observed with high significance.

The Milagrito water Cherenkov telescope operated for over a year. The most probable gamma-ray energy was ~1 TeV and the trigger rate was as high as 400 Hz. We have developed an efficient technique for searching the entire sky for short duration bursts of TeV photons. Such bursts may result from "traditional" gamma-ray bursts that were not in the field-of-view of any other instruments, the evaporation of primordial black holes, or some as yet undiscovered phenomenon. We have begun to search the Milagrito data set for bursts of duration 10 seconds. Here we will present the technique and the expected results. Final results will be presented at the conference.

The Milagro water Cherenkov detector near Los Alamos, New Mexico, has been operated as a sky monitor at energies of a few TeV between February 1997 and April 1998. Serving as a test run for the full Milagro detector, Milagrito has taken data during the strong and long-lasting 1997 flare of Mrk 501. We present results from the analysis of Mrk 501 and compare the excess and background rates with expectations from the detector simulations.

The Milagrito water Cherenkov detector in the Jemez Mountains near Los Alamos, New Mexico took data from February 1997 to April 1998. Milagrito served as a prototype for the larger Milagro detector, which has just begun operations. Milagrito was the first large-aperture gamma-ray detector with sensitivity to gamma rays below 1 TeV. We report here on a search for steady emission from point sources over most of the northern sky using data from Milagrito.