M51 and NGC 5195 is an interacting system that can be explored in great details with ground-based telescopes. The H II regions in M51 were observed using the 2.16 m telescope of the National Astronomical Observatories of the Chinese Academy of Sciences and the 6.5 m Multiple Mirror Telescope with spatial resolution of less than $\sim100$ pc. We obtain a total of 113 spectra across the galaxy and combine the literature data of Croxall et al. to derive a series of physical properties, including the gas-phase extinction, stellar population age, star formation rate (SFR) surface density, and oxygen abundance. The spatial distributions and radial profiles of these properties are investigated in order to study the characteristics of M51 and the clues to the formation and evolution of this galaxy. M51 presents a mild radial extinction gradient. The lower gas-phase extinction in the north spiral arms compared to the south arms are possibly caused by the past encounters with the companion galaxy of NGC 5195. A number of H II regions have the stellar age between 50 and 500 Myr, consistent with the recent interaction history by simulations in the literatures. The SFR surface density presents a mild radial gradient, which is ubiquitous in spiral galaxies. There is a negative metallicity gradient of $-0.08$ dex $R_{e}^{-1}$ in the disk region, which is also commonly found in many spiral galaxies. It is supported by the "inside-out" scenario of galaxy formation. We find a positive abundance gradient of 0.26 dex $R_{e}^{-1}$ in the inner region. There are possible reasons causing the positive gradient, including the freezing of the chemical enrichment due to the star-forming quenching in the bulge and the gas infall and dilution due to the pseudobulge growth and/or galactic interaction.

Providing a possible connection between neutrino emission and gravitational-wave (GW) bursts is important to our understanding of the physical processes that occur when black holes or neutron stars merge. In the Daya Bay experiment, using data collected from December 2011 to August 2017, a search has been performed for electron-antineutrino signals coinciding with detected GW events, including GW150914, GW151012, GW151226, GW170104, GW170608, GW170814, and GW170817. We used three time windows of $\mathrm{\pm 10~s}$, $\mathrm{\pm 500~s}$, and $\mathrm{\pm 1000~s}$ relative to the occurrence of the GW events, and a neutrino energy range of 1.8 to 100 MeV to search for correlated neutrino candidates. The detected electron-antineutrino candidates are consistent with the expected background rates for all the three time windows. Assuming monochromatic spectra, we found upper limits (90% confidence level) on electron-antineutrino fluence of $(1.13~-~2.44) \times 10^{11}~\rm{cm^{-2}}$ at 5 MeV to $8.0 \times 10^{7}~\rm{cm^{-2}}$ at 100 MeV for the three time windows. Under the assumption of a Fermi-Dirac spectrum, the upper limits were found to be $(5.4~-~7.0)\times 10^{9}~\rm{cm^{-2}}$ for the three time windows.

The density structure of the interstellar medium (ISM) determines where stars form and release energy, momentum, and heavy elements, driving galaxy evolution. Density variations are seeded and amplified by gas motion, but the exact nature of this motion is unknown across spatial scale and galactic environment. Although dense star-forming gas likely emerges from a combination of instabilities, convergent flows, and turbulence, establishing the precise origin is challenging because it requires quantifying gas motion over many orders of magnitude in spatial scale. Here we measure the motion of molecular gas in the Milky Way and in nearby galaxy NGC 4321, assembling observations that span an unprecedented spatial dynamic range ($10^{-1}{-}10^3$ pc). We detect ubiquitous velocity fluctuations across all spatial scales and galactic environments. Statistical analysis of these fluctuations indicates how star-forming gas is assembled. We discover oscillatory gas flows with wavelengths ranging from $0.3{-}400$ pc. These flows are coupled to regularly-spaced density enhancements that likely form via gravitational instabilities. We also identify stochastic and scale-free velocity and density fluctuations, consistent with the structure generated in turbulent flows. Our results demonstrate that ISM structure cannot be considered in isolation. Instead, its formation and evolution is controlled by nested, interdependent flows of matter covering many orders of magnitude in spatial scale.

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

We present kinematic orientations and high resolution (150 pc) rotation curves for 67 main sequence star-forming galaxies surveyed in CO (2-1) emission by PHANGS-ALMA. Our measurements are based on the application of a new fitting method tailored to CO velocity fields. Our approach identifies an optimal global orientation as a way to reduce the impact of non-axisymmetric (bar and spiral) features and the uneven spatial sampling characteristic of CO emission in the inner regions of nearby galaxies. The method performs especially well when applied to the large number of independent lines-of-sight contained in the PHANGS CO velocity fields mapped at 1'' resolution. The high resolution rotation curves fitted to these data are sensitive probes of mass distribution in the inner regions of these galaxies. We use the inner slope as well as the amplitude of our fitted rotation curves to demonstrate that CO is a reliable global dynamical mass tracer. From the consistency between photometric orientations from the literature and kinematic orientations determined with our method, we infer that the shapes of stellar disks in the mass range of log($\rm M_{\star}(M_{\odot})$)=9.0-10.9 probed by our sample are very close to circular and have uniform thickness.

This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

We report the discovery of a highly dispersed fast radio burst, FRB~181123, from an analysis of $\sim$1500~hr of drift-scan survey data taken using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). The pulse has three distinct emission components, which vary with frequency across our 1.0--1.5~GHz observing band. We measure the peak flux density to be $>0.065$~Jy and the corresponding fluence $>0.2$~Jy~ms. Based on the observed dispersion measure of 1812~cm$^{-3}$~pc, we infer a redshift of $\sim 1.9$. From this, we estimate the peak luminosity and isotropic energy to be $\lesssim 2\times10^{43}$~erg~s$^{-1}$ and $\lesssim 2\times10^{40}$~erg, respectively. With only one FRB from the survey detected so far, our constraints on the event rate are limited. We derive a 95\% confidence lower limit for the event rate of 900 FRBs per day for FRBs with fluences $>0.025$~Jy~ms. We performed follow-up observations of the source with FAST for four hours and have not found a repeated burst. We discuss the implications of this discovery for our understanding of the physical mechanisms of FRBs.

In this paper, a new measurement system based on integration method is presented,which can provide all-weather dependability and higher precision for the measurement of FAST's feed support system. The measurement system consists of three types of measuring equipments, and a processing software with the core data fusion algorithm. The Strapdown Inertial Navigation System(SINS) can autonomously measure the position, speed and attitude of the carrier. Its own shortcoming is the measurement data diverges rapidly over time. SINS must combine the Global Positioning System(GPS) and the Total Station(TS) to obtain high-precision measurement data. Kalman filtering algorithm is adopted for the integration measurement system, which is an optimal algorithm to estimate the measurement errors. To evaluate the performance, series of tests are carried out. For the feed cabin, the maximum RMS of the position is 14.56mm, the maximum RMS of the attitude is 0.095{\deg}, these value are less than 15mm and 0.1{\deg} as the precision for measuring the feed cabin. For the Stewart manipulator, the maximum RMS of the position is 2.99mm, the maximum RMS of the at titude is 0.093{\deg}, these value are less than 3mm and 0.1{\deg} as the precision for measuring the Stewart manipulator. As a result, the new measurement meets the requirement of measurement precision for FAST's feed support system.

We compare the observed turbulent pressure in molecular gas, $P_\mathrm{turb}$, to the required pressure for the interstellar gas to stay in equilibrium in the gravitational potential of a galaxy, $P_\mathrm{DE}$. To do this, we combine arcsecond resolution CO data from PHANGS-ALMA with multi-wavelength data that traces the atomic gas, stellar structure, and star formation rate (SFR) for 28 nearby star-forming galaxies. We find that $P_\mathrm{turb}$ correlates with, but almost always exceeds the estimated $P_\mathrm{DE}$ on kiloparsec scales. This indicates that the molecular gas is over-pressurized relative to the large-scale environment. We show that this over-pressurization can be explained by the clumpy nature of molecular gas; a revised estimate of $P_\mathrm{DE}$ on cloud scales, which accounts for molecular gas self-gravity, external gravity, and ambient pressure, agrees well with the observed $P_\mathrm{turb}$ in galaxy disks. We also find that molecular gas with cloud-scale ${P_\mathrm{turb}}\approx{P_\mathrm{DE}}\gtrsim{10^5\,k_\mathrm{B}\,\mathrm{K\,cm^{-3}}}$ in our sample is more likely to be self-gravitating, whereas gas at lower pressure appears more influenced by ambient pressure and/or external gravity. Furthermore, we show that the ratio between $P_\mathrm{turb}$ and the observed SFR surface density, $\Sigma_\mathrm{SFR}$, is compatible with stellar feedback-driven momentum injection in most cases, while a subset of the regions may show evidence of turbulence driven by additional sources. The correlation between $\Sigma_\mathrm{SFR}$ and kpc-scale $P_\mathrm{DE}$ in galaxy disks is consistent with the expectation from self-regulated star formation models. Finally, we confirm the empirical correlation between molecular-to-atomic gas ratio and kpc-scale $P_\mathrm{DE}$ reported in previous works.

In order to measure the Universal Time (UT1) in real time, National Time Service Center (NTSC) has built a VGOS-like (VLBI Global Observing System) broadband VLBI network, which includes three 13-m radio telescopes located in Jilin, Sanya and Kashi, and a data analysis center in Xi'an. Each station is equipped with a highly stable hydrogen atomic clock and a self-developed VLBI backend, and is co-located with two GPS receivers. This VGOS-like VLBI network may play an important role in improving the Chinese broadband VLBI technology and making valuable contributions to domestic VLBI measurements of UT1. In this paper, we introduce the specifications of this VLBI network, and present the UT1 measurements at C-band conducted in 2018 using the Jilin-Kashi baseline of this network. The comparisons between our UT1 estimates and those provided by IERS suggest that the NTSC VLBI network is capable to determine UT1 accurate at the level of 58.8 microseconds.

The Five-hundred-meter Aperture Spherical radio Telescope (FAST) passed national acceptance and is taking pilot cycle of 'Shared-Risk' observations. The 19-beam receiver covering 1.05-1.45 GHz was used for most of these observations. The electronics gain fluctuation of the system is better than 1\% over 3.5 hours, enabling enough stability for observations. Pointing accuracy, aperture efficiency and system temperature are three key parameters of FAST. The measured standard deviation of pointing accuracy is 7.9$''$, which satisfies the initial design of FAST. When zenith angle is less than 26.4$^\circ$, the aperture efficiency and system temperature around 1.4 GHz are $\sim$ 0.63 and less than 24 K for central beam, respectively. The measured value of these two parameters are better than designed value of 0.6 and 25 K, respectively. The sensitivity and stability of the 19-beam backend are confirmed to satisfy expectation by spectral HI observations toward N672 and polarization observations toward 3C286. The performance allows FAST to take sensitive observations in various scientific goals, from studies of pulsar to galaxy evolution.

The chemical abundances of spiral galaxies, as probed by HII regions across their disks, are key to understanding the evolution of galaxies over a wide range of environments. We present LBT/MODS spectra of 52 HII regions in NGC3184 as part of the CHemical Abundances Of Spirals (CHAOS) project. We explore the direct-method gas-phase abundance trends for the first four CHAOS galaxies, using temperature measurements from one or more auroral line detections in 190 individual HII regions. We find the dispersion in Te-Te relationships is dependent on ionization, as characterized by F_5007/F_3727, and so recommend ionization-based temperature priorities for abundance calculations. We confirm our previous results that [NII] and [SIII] provide the most robust measures of electron temperature in low-ionization zones, while [OIII] provides reliable electron temperatures in high-ionization nebula. We measure relative and absolute abundances for O, N, S, Ar, and Ne. The four CHAOS galaxies marginally conform with a universal O/H gradient, as found by empirical IFU studies when plotted relative to effective radius. However, after adjusting for vertical offsets, we find a tight universal N/O gradient of alpha_N/O = -0.33 dex/Re with sigma_tot. = 0.08 for Rg/Re < 2.0, where N is dominated by secondary production. Despite this tight universal N/O gradient, the scatter in the N/O-O/H relationship is significant. Interestingly, the scatter is similar when N/O is plotted relative to O/H or S/H. The observable ionic states of S probe lower ionization and excitation energies than O, which might be more appropriate for characterizing abundances in metal-rich HII regions.

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

The expansion of hot electrons in flaring magnetic loops is crucial to understanding the dynamics of solar flares. In this paper we investigate, for the first time, the transport of hot electrons in a magnetic mirror field based on a 1-D particle-in-cell (PIC) simulation. The hot electrons with small pitch angle transport into the cold plasma, which leads to the generation of Langmuir waves in the cold plasma and ion acoustic waves in the hot plasma. The large pitch angle electrons can be confined by the magnetic mirror, resulting in the different evolution time scale between electron parallel and perpendicular temperature. This will cause the formation of electron temperature anisotropy, which can generate the whistler waves near the interface between hot electrons and cold electrons. The whistler waves can scatter the large pitch angle electrons to smaller value through the cyclotron resonance, leading to electrons escaping from the hot region. These results indicate that the whistler waves may play an important role in the transport of electrons in flaring magnetic loops. The findings from this study provide some new insights to understand the electron dynamics of solar flares.

We present CO(J=1-0) and/or CO(J=2-1) spectroscopy for 31 galaxies selected from the ongoing MaNGA survey, obtained with multiple telescopes. This sample is combined with CO observations from the literature to study the correlation of the CO luminosities ($L_{\rm CO(1-0)}$) with the mid-infrared luminosities at 12 ($L_{12 \mu m}$) and 22 $\mu$m ($L_{\rm 22 \mu m}$), as well as the dependence of the residuals on a variety of galaxy properties. The correlation with $L_{\rm 12 \mu m}$ is tighter and more linear, but galaxies with relatively low stellar masses and blue colors fall significantly below the mean $L_{\rm CO(1-0)}-L_{\rm 12\mu m}$ relation. We propose a new estimator of the CO(1-0) luminosity (and thus the total molecular gas mass) that is a linear combination of three parameters: $L_{\rm 12 \mu m}$, $M_\ast$ and $g-r$. We show that, with a scatter of only 0.18 dex in log $(L_{\rm CO(1-0)})$, this estimator provides unbiased estimates for galaxies of different properties and types. An immediate application of this estimator to a compiled sample of galaxies with only CO(J=2-1) observations yields a distribution of the CO(J=2-1) to CO(J=1-0) luminosity ratios ($R21$) that agrees well with the distribution of real observations, in terms of both the median and the shape. Application of our estimator to the current MaNGA sample reveals a gas-poor population of galaxies that are predominantly early-type and show no correlation between molecular gas-to-stellar mass ratio and star formation rate, in contrast to gas-rich galaxies. We also provide alternative estimators with similar scatters, based on $r$ and/or $z$ band luminosities instead of $M_\ast$. These estimators serve as cheap and convenient $M_{\rm mol}$ proxies to be potentially applied to large samples of galaxies, thus allowing statistical studies of gas-related processes of galaxies.

It remains a major challenge to derive a theory of cloud-scale ($\lesssim100$ pc) star formation and feedback, describing how galaxies convert gas into stars as a function of the galactic environment. Progress has been hampered by a lack of robust empirical constraints on the giant molecular cloud (GMC) lifecycle. We address this problem by systematically applying a new statistical method for measuring the evolutionary timeline of the GMC lifecycle, star formation, and feedback to a sample of nine nearby disc galaxies, observed as part of the PHANGS-ALMA survey. We measure the spatially-resolved ($\sim100$ pc) CO-to-H$\alpha$ flux ratio and find a universal de-correlation between molecular gas and young stars on GMC scales, allowing us to quantify the underlying evolutionary timeline. GMC lifetimes are short, typically 10-30 Myr, and exhibit environmental variation, between and within galaxies. At kpc-scale molecular gas surface densities $\Sigma_{\rm H_2}\geqslant8$M$_{\odot}$pc$^{-2}$, the GMC lifetime correlates with time-scales for galactic dynamical processes, whereas at $\Sigma_{\rm H_2}\leqslant8$M$_{\odot}$pc$^{-2}$ GMCs decouple from galactic dynamics and live for an internal dynamical time-scale. After a long inert phase without massive star formation traced by H$\alpha$ (75-90% of the cloud lifetime), GMCs disperse within just 1-5 Myr once massive stars emerge. The dispersal is most likely due to early stellar feedback, causing GMCs to achieve integrated star formation efficiencies of 4-10% These results show that galactic star formation is governed by cloud-scale, environmentally-dependent, dynamical processes driving rapid evolutionary cycling. GMCs and HII regions are the fundamental units undergoing these lifecycles, with mean separations of 100-300 pc in star-forming discs. Future work should characterise the multi-scale physics and mass flows driving these lifecycles.

Cloud-scale surveys of molecular gas reveal the link between molecular clouds properties and star formation (SF) across a range of galactic environments. Cloud populations in galaxy disks are considered to be representative of the `normal' SF. At high resolution, however, clouds with exceptional gas properties and SF activity may also be observed in normal disk environments. In this paper, we study the brightest cloud traced in CO emission in the disk of NGC628. The cloud is spatially coincident with an extremely bright HII region. We characterize its molecular gas properties and investigate how feedback and large-scale processes influence the properties of the molecular gas. High resolution CO ALMA observations are used to characterize its mass and dynamical state, which are compared to other clouds in NGC628. A LVG analysis is used to constrain the beam-diluted density and temperature of the molecular gas. We analyze the MUSE spectrum using Starburst99 to characterize the young stellar population associated with the HII region. The cloud is massive ($1-2\times10^7$M$_{\odot}$), with a beam-diluted density of $n_{\rm H_2}=5\times10^4$ cm$^{-3}$. It has a low virial parameter, suggesting that its CO emission may be overluminous due to heating by the HII region. A young ($2-4$ Myr), massive $3\times10^{5}$ M$_{\odot}$ stellar population is associated. We argue that the cloud is currently being destroyed by feedback from young massive stars. Due to the cloud's large mass, this phase of the cloud's evolution is long enough for the impact of feedback on the excitation of the gas to be observed. Its high mass may be related to its location at a spiral co-rotation radius, where gas experiences reduced galactic shear compared to other regions of the disk, and receives a sustained inflow of gas that can promote the cloud's mass growth.

The processes regulating star formation in galaxies are thought to act across a hierarchy of spatial scales. To connect extragalactic star formation relations from global and kpc-scale measurements to recent cloud-scale resolution studies, we have developed a simple, robust method that quantifies the scale dependence of the relative spatial distributions of molecular gas and recent star formation. In this paper, we apply this method to eight galaxies with roughly 1 arcsec resolution molecular gas imaging from the PHANGS-ALMA and PAWS surveys that have matched resolution, high quality narrowband Halpha imaging. At a common scale of 140pc, our massive (log(Mstar/Msun)=9.3-10.7), normally star-forming (SFR/Msun/yr=0.3-5.9) galaxies exhibit a significant reservoir of quiescent molecular gas not associated with star formation as traced by Halpha emission. Galactic structures act as backbones for both molecular and HII region distributions. As we degrade the spatial resolution, the quiescent molecular gas disappears, with the most rapid changes occurring for resolutions up to about 0.5kpc. As the resolution becomes poorer, the morphological features become indistinct for spatial scales larger than about 1kpc. The method is a promising tool to search for relationships between the quiescent or star-forming molecular reservoir and galaxy properties, but requires a larger sample size to identify robust correlations between the star-forming molecular gas fraction and global galaxy parameters.

The POLAR detector is a space based Gamma-Ray Burst (GRB) polarimeter sensitive in the 15-500 keV energy range. Apart from its main scientific goal as a Gamma-Ray Burst polarimeter it is also able to detect photons from pulsars in orbit. By using the six-months in-orbit observation data, significant pulsation from the PSR B0531+21 (Crab pulsar) was obtained. In this work, we present the precise timing analysis of the Crab pulsar, together with a phase-resolved spectroscopic study using a joint-fitting method adapted for wide field of view instruments like POLAR. By using single power law fitting over the pulsed phase, we obtained spectral indices ranging from 1.718 to 2.315, and confirmed the spectral evolution in a reverse S shape which is homogenous with results from other missions over broadband. We will also show, based on the POLAR in-orbit performance and Geant4 Monte-Carlo simulation, the inferred capabilities of POLAR-2, the proposed follow-up mission of POLAR on board the China Space Station (CSS), for pulsars studies.

Simultaneous $\gamma$-ray measurements of gamma-ray burst (GRB) spectra and polarization offer a unique way to determine the underlying emission mechanism(s) in these objects as well as probing the particle acceleration mechanism(s) that lead to the observed $\gamma$-ray emission. Herein we examine the jointly-observed data from POLAR and GBM of GRB 170114A to determine its spectral and polarization properties and seek to understand the emission processes that generate these observations. We aim to develop an extensible and statistically sound framework for these types of measurements applicable to other instruments. We leverage the existing 3ML analysis framework to develop a new analysis pipeline for simultaneously modeling the spectral and polarization data. We derive the proper Poisson likelihood for $\gamma$-ray polarization measurements in the presence of background. The developed framework is publicly available for similar measurements with other $\gamma$-ray polarimeters. The data are analyzed within a Bayesian probabilistic context and the spectral data from both instruments are simultaneously modeled with a physical, numerical synchrotron code. The spectral modeling of the data is consistent with a synchrotron photon model as has been found in a majority of similarly analyzed single-pulse GRBs. The polarization results reveal a slight trend of growing polarization in time reaching values of ~30% at the temporal peak of the emission. Additionally, it is observed that the polarization angle evolves with time throughout the emission. These results hint at a synchrotron origin of the emission but further observations of many GRBs are required to verify these evolutionary trends. Furthermore, we encourage the development of time-resolved polarization models for the prompt emission of GRBs as the current models are not predictive enough to enable a full modeling of our current data.

POLAR is a compact space-borne detector initially designed to measure the polarization of hard X-rays emitted from Gamma-Ray Bursts in the energy range 50-500keV. This instrument was launched successfully onboard the Chinese space laboratory Tiangong-2 (TG-2) on 2016 September 15. After being switched on a few days later, tens of gigabytes of raw detection data were produced in-orbit by POLAR and transferred to the ground every day. Before the launch date, a full pipeline and related software were designed and developed for the purpose of quickly pre-processing all the raw data from POLAR, which include both science data and engineering data, then to generate the high level scientific data products that are suitable for later science analysis. This pipeline has been successfully applied for use by the POLAR Science Data Center in the Institute of High Energy Physics (IHEP) after POLAR was launched and switched on. A detailed introduction to the pipeline and some of the core relevant algorithms are presented in this paper.

We describe PSR J1926-0652, a pulsar recently discovered with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Using sensitive single-pulse detections from FAST and long-term timing observations from the Parkes 64-m radio telescope, we probed phenomena on both long and short time scales. The FAST observations covered a wide frequency range from 270 to 800 MHz, enabling individual pulses to be studied in detail. The pulsar exhibits at least four profile components, short-term nulling lasting from 4 to 450 pulses, complex subpulse drifting behaviours and intermittency on scales of tens of minutes. While the average band spacing P3 is relatively constant across different bursts and components, significant variations in the separation of adjacent bands are seen, especially near the beginning and end of a burst. Band shapes and slopes are quite variable, especially for the trailing components and for the shorter bursts. We show that for each burst the last detectable pulse prior to emission ceasing has different properties compared to other pulses. These complexities pose challenges for the classic carousel-type models.

We describe a next major frontier in observational studies of galaxy evolution and star formation: linking the physical conditions in the cold, star-forming interstellar medium to host galaxy and local environment. The integrated gas content of galaxies has been surveyed extensively over the last decades. The link between environment and cold gas density, turbulence, excitation, dynamical state, and chemical makeup remain far less well understood. We know that these properties do vary dramatically and theoretical work posits a strong connection between the state of the gas, its ability to form stars, and the impact of stellar feedback. A next major step in the field will be to use sensitive cm-, mm-, and submm-wave spectroscopy and high resolution spectroscopic imaging to survey the state of cold gas across the whole local galaxy population. Such observations have pushed the capabilities of the current generation of telescopes. We highlight three critical elements for progress in the next decade: (1) robust support and aggressive development of ALMA, (2) the deployment of very large heterodyne receiver arrays on single dish telescopes, and (3) development of a new interferometric array that dramatically improves on current capabilities at cm- and mm-wavelengths (~ 1-115 GHz).

The Five-hundred-meter Aperture Spherical radio Telescope(FAST) is expected to complete its commissioning in 2019. FAST will soon begin the Commensal Radio Astronomy FasT Survey(CRAFTS), a novel and unprecedented commensal drift scan survey of the entire sky visible from FAST. The goal of CRAFTS is to cover more than 20000 $deg^{2}$ and reach redshift up to about 0.35. We provide empirical measurements of the beam size and sensitivity of FAST across the 1.05 to 1.45 GHz frequency range of the FAST L-band Array of 19-beams(FLAN). Using a simulated HI-galaxy catalogue based on the HI Mass Function(HIMF), we estimate the number of galaxies that CRAFTS may detect. At redshifts below 0.35, over $6\, \times \, 10^{5}$ HI galaxies may be detected. Below the redshift of 0.07, the CRAFTS HIMF will be complete above a mass threshold of $10^{9.5}\,M_{\odot}$. FAST will be able to investigate the environmental and redshift dependence of the HIMF to an unprecedented depth, shedding light onto the missing baryon and missing satellite problems.

To assist with the commissioning (Jiang et al. 2019) of the Five-hundred-meter Aperture Spherical radio Telescope (FAST), we performed a pulsar search, with the primary goal of developing and testing the pulsar data acquisition and processing pipelines. We tested and used three pipelines, two (P1 and P2 hereafter) searched for the periodic signature of pulsars whereas the other one was used to search for bright single pulses (P3 hereafter). A pulsar candidate was discovered in the observation on the 22nd August, 2017, and later confirmed by the Parkes radio telescope on the 10th September, 2017. This pulsar, named PSR J1900-0134, was the first pulsar discovered by FAST. The pulsar has a pulse period of 1.8 s and a dispersion measure (DM) of 188\,pc\,cm$^{-3}$.

The Five-hundred-meter Aperture Spherical radio Telescope (FAST) was completed with its main structure installed on September 25, 2016, after which it entered the commissioning phase. This paper aims to introduce the commissioning progress of the FAST over the past two years. To improve its operational reliability and ensure effective observation time, FAST has been equipped with a real-time information system for the active reflector system and hierarchical commissioning scheme for the feed support system, which ultimately achieves safe operation of the two systems. For meeting the high-performance indices, a high-precision measurement system was set up based on the effective control methods that were implemented for the active reflector system and feed support system. Since the commissioning of the FAST, a low-frequency ultra-wideband receiver and 19-beam 1.05-1.45 GHz receiver have been mainly used. Telescope efficiency, pointing accuracy, and system noise temperature were completely tested and ultimately achieved the acceptance indices of the telescope. The FAST has been in the process of national acceptance preparations and has begun to search for pulsars. In the future, it will still strive to improve its capabilities and expand its application prospects.

We compare the properties of clouds in simulated M33 galaxies to those observed in the real M33. We apply a friends of friends algorithm and CPROPS to identify clouds, as well as a pixel by pixel analysis. We obtain very good agreement between the number of clouds, and maximum mass of clouds. Both are lower than occurs for a Milky Way-type galaxy and thus are a function of the surface density, size and galactic potential of M33. We reproduce the observed dependence of molecular cloud properties on radius in the simulations, and find this is due to the variation in gas surface density with radius. The cloud spectra also show good agreement between the simulations and observations, but the exact slope and shape of the spectra depends on the algorithm used to find clouds, and the range of cloud masses included when fitting the slope. Properties such as cloud angular momentum, velocity dispersions and virial relation are also in good agreement between the simulations and observations, but do not necessarily distinguish between simulations of M33 and other galaxy simulations. Our results are not strongly dependent on the level of feedback used here (10 and 20%) although they suggest that 15% feedback efficiency may be optimal. Overall our results suggest that the molecular cloud properties are primarily dependent on the gas and mass surface density, and less dependent on the localised physics such as the details of stellar feedback, or the numerical code used.

Gamma-ray bursts are the strongest explosions in the Universe since the Big Bang, believed to be produced either in forming black holes at the end of massive star evolution or merging of compact objects. Spectral and timing properties of gamma-ray bursts suggest that the observed bright gamma-rays are produced in the most relativistic jets in the Universe; however, the physical properties, especially the structure and magnetic topologies in the jets are still not well known, despite several decades of studies. It is widely believed that precise measurements of the polarization properties of gamma-ray bursts should provide crucial information on the highly relativistic jets. As a result there have been many reports of gamma-ray burst polarization measurements with diverse results, see, however many such measurements suffered from substantial uncertainties, mostly systematic. After the first successful measurements by the GAP and COSI instruments, here we report a statistically meaningful sample of precise polarization measurements, obtained with the dedicated gamma-ray burst polarimeter, POLAR onboard China's Tiangong-2 spacelab. Our results suggest that the gamma-ray emission is at most polarized at a level lower than some popular models have predicted; although our results also show intrapulse evolution of the polarization angle. This indicates that the low polarization degrees could be due to an evolving polarization angle during a gamma-ray burst.

Twenty years have passed since first light for the Sloan Digital Sky Survey (SDSS). Here, we release data taken by the fourth phase of SDSS (SDSS-IV) across its first three years of operation (July 2014-July 2017). This is the third data release for SDSS-IV, and the fifteenth from SDSS (Data Release Fifteen; DR15). New data come from MaNGA - we release 4824 datacubes, as well as the first stellar spectra in the MaNGA Stellar Library (MaStar), the first set of survey-supported analysis products (e.g. stellar and gas kinematics, emission line, and other maps) from the MaNGA Data Analysis Pipeline (DAP), and a new data visualisation and access tool we call "Marvin". The next data release, DR16, will include new data from both APOGEE-2 and eBOSS; those surveys release no new data here, but we document updates and corrections to their data processing pipelines. The release is cumulative; it also includes the most recent reductions and calibrations of all data taken by SDSS since first light. In this paper we describe the location and format of the data and tools and cite technical references describing how it was obtained and processed. The SDSS website (www.sdss.org) has also been updated, providing links to data downloads, tutorials and examples of data use. While SDSS-IV will continue to collect astronomical data until 2020, and will be followed by SDSS-V (2020-2025), we end this paper by describing plans to ensure the sustainability of the SDSS data archive for many years beyond the collection of data.

We use ALMA and IRAM 30-m telescope data to investigate the relationship between the spectroscopically-traced dense gas fraction and the cloud-scale (120 pc) molecular gas surface density in five nearby, star-forming galaxies. We estimate the dense gas mass fraction at 650 pc and 2800 pc scales using the ratio of HCN (1-0) to CO (1-0) emission. We then use high resolution (120 pc) CO (2-1) maps to calculate the mass-weighted average molecular gas surface density within 650 pc or 2770 pc beam where the dense gas fraction is estimated. On average, the dense gas fraction correlates with the mass-weighted average molecular gas surface density. Thus, parts of a galaxy with higher mean cloud-scale gas surface density also appear to have a larger fraction of dense gas. The normalization and slope of the correlation do vary from galaxy to galaxy and with the size of the regions studied. This correlation is consistent with a scenario where the large-scale environment sets the gas volume density distribution, and this distribution manifests in both the cloud-scale surface density and the dense gas mass fraction.

The Serpens filament, as one of the nearest infrared dark clouds, is regarded as a pristine filament at a very early evolutionary stage of star formation. In order to study its molecular content and dynamical state, we mapped this filament in seven species. Among them, HCO$^{+}$, HNC, HCN, and CS show self-absorption, while C$^{18}$O is most sensitive to the filamentary structure. A kinematic analysis demonstrates that this filament forms a velocity-coherent (trans-)sonic structure, a large part of which is one of the most quiescent regions in the Serpens cloud. Widespread C$^{18}$O depletion is found throughout the Serpens filament. Based on the Herschel dust-derived H$_{2}$ column density map, the line mass of the filament is 36--41~M$_{\odot}$~pc$^{-1}$, and its full width at half maximum is 0.17$\pm$0.01~pc, while its length is ~1.6~pc. The inner radial column density profile of this filament can be well fitted with a Plummer profile with an exponent of 2.2$\pm$0.1, a scale radius of $0.018\pm 0.003$ pc, and a central density of $(4.0\pm 0.8)\times 10^{4}$~cm$^{-3}$. The Serpens filament appears to be slightly supercritical. The widespread blue-skewed HNC and CS line profiles and HCN hyperfine line anomalies across this filament indicate radial infall in parts of the Serpens filament. C$^{18}$O velocity gradients also indicate accretion flows along the filament. The velocity and density structures suggest that such accretion flows are likely due to a longitudinal collapse parallel to the filament's long axis. Both the radial infall rate and the longitudinal accretion rate along the Serpens filament are lower than all previously reported values in other filaments. This indicates that the Serpens filament lies at an early evolutionary stage when collapse has just begun, or that thermal and non-thermal support are effective in providing support against gravity.

Star formation is a multi-scale process that requires tracing cloud formation and stellar feedback within the local (<kpc) and global galaxy environment. We present first results from two large observing programs on ALMA and VLT/MUSE, mapping cloud scales (1arcsec = 47pc) in both molecular gas and star forming tracers across 90 kpc^2 of the central disk of NGC 628 to probe the physics of star formation. Systematic spatial offsets between molecular clouds and HII regions illustrate the time evolution of star-forming regions. Using uniform sampling of both maps on 50-500 pc scales, we infer molecular gas depletion times of 1-3 Gyr, but also find the increase of scatter in the star formation relation on small scales is consistent with gas and HII regions being only weakly correlated at the cloud (50 pc) scale. This implies a short overlap phase for molecular clouds and HII regions, which we test by directly matching our catalog of 1502 HII regions and 738 GMCs. We uncover only 74 objects in the overlap phase, and we find depletion times >1 Gyr, significantly longer than previously reported for individual star-forming clouds in the Milky Way. Finally, we find no clear trends that relate variations in the depletion time observed on 500 pc scales to physical drivers (metallicity, molecular and stellar mass surface density, molecular gas boundedness) on 50 pc scales.

We estimate the star formation efficiency per gravitational free fall time, $\epsilon_{\rm ff}$, from observations of nearby galaxies with resolution matched to the typical size of a Giant Molecular Cloud. This quantity, $\epsilon_{\rm ff}$, is theoretically important but so far has only been measured for Milky Way clouds or inferred indirectly in a few other galaxies. Using new, high resolution CO imaging from the PHANGS-ALMA survey, we estimate the gravitational free-fall time at 60 to 120 pc resolution, and contrast this with the local molecular gas depletion time to estimate $\epsilon_{\rm ff}$. Assuming a constant thickness of the molecular gas layer ($H = 100$ pc) across the whole sample, the median value of $\epsilon_{\rm ff}$ in our sample is $0.7\%$. We find a mild scale-dependence, with higher $\epsilon_{\rm ff}$ measured at coarser resolution. Individual galaxies show different values of $\epsilon_{\rm ff}$, with the median $\epsilon_{\rm ff}$ ranging from $0.3\%$ to $2.6\%$. We find the highest $\epsilon_{\rm ff}$ in our lowest mass targets, reflecting both long free-fall times and short depletion times, though we caution that both measurements are subject to biases in low mass galaxies. We estimate the key systematic uncertainties, and show the dominant uncertainty to be the estimated line-of-sight depth through the molecular gas layer and the choice of star formation tracers.

POLAR is a compact wide-field space-borne detector for precise measurements of the linear polarisation of hard X-rays emitted by transient sources in the energy range from 50 keV to 500 keV. It consists of a 40$\times$40 array of plastic scintillator bars used as a detection material. The bars are grouped in 25 detector modules. The energy range sensitivity of POLAR is optimized to match with the prompt emission photons from the gamma-ray bursts (GRBs). Polarization measurements of the prompt emission would probe source geometries, emission mechanisms and magnetic structures in GRB jets. The instrument can also detect hard X-rays from solar flares and be used for precise measurement of their polarisation. POLAR was launched into a low Earth orbit on-board the Chinese space-lab TG-2 on September 15th, 2016. To achieve high accuracies in polarisation measurements it is essential to assure both before and during the flight a precise energy calibration. Such calibrations are performed with four low activity $^{22}$Na radioactive sources placed inside the instrument. Energy conversion factors are related to Compton edge positions from the collinear annihilation photons from the sources. This paper presents main principles of the in-flight calibration, describes studies of the method based on Monte Carlo simulations and its laboratory verification and finally provides some observation results based on the in-flight data analysis.

POLAR is a compact space-borne detector designed to perform reliable measurements of the polarization for transient sources like Gamma-Ray Bursts in the energy range 50-500keV. The instrument works based on the Compton Scattering principle with the plastic scintillators as the main detection material along with the multi-anode photomultiplier tube. POLAR has been launched successfully onboard the Chinese space laboratory TG-2 on 15th September, 2016. In order to reliably reconstruct the polarization information a highly detailed understanding of the instrument is required for both data analysis and Monte Carlo studies. For this purpose a full study of the in-orbit performance was performed in order to obtain the instrument calibration parameters such as noise, pedestal, gain nonlinearity of the electronics, threshold, crosstalk and gain, as well as the effect of temperature on the above parameters. Furthermore the relationship between gain and high voltage of the multi-anode photomultiplier tube has been studied and the errors on all measurement values are presented. Finally the typical systematic error on polarization measurements of Gamma-Ray Bursts due to the measurement error of the calibration parameters are estimated using Monte Carlo simulations.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) has been in operation since July 2014. This paper describes the second data release from this phase, and the fourteenth from SDSS overall (making this, Data Release Fourteen or DR14). This release makes public data taken by SDSS-IV in its first two years of operation (July 2014-2016). Like all previous SDSS releases, DR14 is cumulative, including the most recent reductions and calibrations of all data taken by SDSS since the first phase began operations in 2000. New in DR14 is the first public release of data from the extended Baryon Oscillation Spectroscopic Survey (eBOSS); the first data from the second phase of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE-2), including stellar parameter estimates from an innovative data driven machine learning algorithm known as "The Cannon"; and almost twice as many data cubes from the Mapping Nearby Galaxies at APO (MaNGA) survey as were in the previous release (N = 2812 in total). This paper describes the location and format of the publicly available data from SDSS-IV surveys. We provide references to the important technical papers describing how these data have been taken (both targeting and observation details) and processed for scientific use. The SDSS website (www.sdss.org) has been updated for this release, and provides links to data downloads, as well as tutorials and examples of data use. SDSS-IV is planning to continue to collect astronomical data until 2020, and will be followed by SDSS-V.

We measure the velocity dispersion, $\sigma$, and surface density, $\Sigma$, of the molecular gas in nearby galaxies from CO spectral line cubes with spatial resolution $45$-$120$ pc, matched to the size of individual giant molecular clouds. Combining $11$ galaxies from the PHANGS-ALMA survey with $4$ targets from the literature, we characterize ${\sim}30,000$ independent sightlines where CO is detected at good significance. $\Sigma$ and $\sigma$ show a strong positive correlation, with the best-fit power law slope close to the expected value for resolved, self-gravitating clouds. This indicates only weak variation in the virial parameter $\alpha_\mathrm{vir}\propto\sigma^2/\Sigma$, which is ${\sim}1.5$-$3.0$ for most galaxies. We do, however, observe enormous variation in the internal turbulent pressure $P_\mathrm{turb}\propto\Sigma\,\sigma^2$, which spans ${\sim}5\rm\;dex$ across our sample. We find $\Sigma$, $\sigma$, and $P_\mathrm{turb}$ to be systematically larger in more massive galaxies. The same quantities appear enhanced in the central kpc of strongly barred galaxies relative to their disks. Based on sensitive maps of M31 and M33, the slope of the $\sigma$-$\Sigma$ relation flattens at $\Sigma\lesssim10\rm\;M_\odot\,pc^{-2}$, leading to high $\sigma$ for a given $\Sigma$ and high apparent $\alpha_\mathrm{vir}$. This echoes results found in the Milky Way, and likely originates from a combination of lower beam filling factors and a stronger influence of local environment on the dynamical state of molecular gas in the low density regime.

Rapid response and short time latency are very important for Time Domain Astronomy, such as the observations of Gamma-ray Bursts (GRBs) and electromagnetic (EM) counterparts of gravitational waves (GWs). Based on the near real-time Fermi/GBM data, we developed a low-latency pipeline to automatically calculate the temporal and spectral properties of GRBs. With this pipeline, some important parameters can be obtained, such as T90 and fluence, within ~20 minutes after the GRB trigger. For ~90% GRBs, T90 and fluence are consistent with the GBM catalog results within 2 sigma errors. This pipeline has been used by the Gamma-ray Bursts Polarimeter (POLAR) and the Insight Hard X-ray Modulation Telescope (Insight-HXMT) to follow up the bursts of interest. For GRB 170817A, the first EM counterpart of GW events detected by Fermi/GBM and INTEGRAL/SPI-ACS, the pipeline gave T90 and spectral information in 21 minutes after the GBM trigger, providing important information for POLAR and Insight-HXMT observations.

The POLAR detector is a space based Gamma Ray Burst (GRB) polarimeter with a wide field of view, which covers almost half the sky. The instrument uses Compton scattering of gamma rays on a plastic scintillator hodoscope to measure the polarization of the incoming photons. The instrument has been successfully launched on board of the Chinese space laboratory Tiangong~2 on September 15, 2016. The construction of the instrument components is described in this article. Details are provided on problems encountered during the construction phase and their solutions. Initial performance of the instrument in orbit is as expected from ground tests and Monte Carlo simulation.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) began observations in July 2014. It pursues three core programs: APOGEE-2, MaNGA, and eBOSS. In addition, eBOSS contains two major subprograms: TDSS and SPIDERS. This paper describes the first data release from SDSS-IV, Data Release 13 (DR13), which contains new data, reanalysis of existing data sets and, like all SDSS data releases, is inclusive of previously released data. DR13 makes publicly available 1390 spatially resolved integral field unit observations of nearby galaxies from MaNGA, the first data released from this survey. It includes new observations from eBOSS, completing SEQUELS. In addition to targeting galaxies and quasars, SEQUELS also targeted variability-selected objects from TDSS and X-ray selected objects from SPIDERS. DR13 includes new reductions of the SDSS-III BOSS data, improving the spectrophotometric calibration and redshift classification. DR13 releases new reductions of the APOGEE-1 data from SDSS-III, with abundances of elements not previously included and improved stellar parameters for dwarf stars and cooler stars. For the SDSS imaging data, DR13 provides new, more robust and precise photometric calibrations. Several value-added catalogs are being released in tandem with DR13, in particular target catalogs relevant for eBOSS, TDSS, and SPIDERS, and an updated red-clump catalog for APOGEE. This paper describes the location and format of the data now publicly available, as well as providing references to the important technical papers that describe the targeting, observing, and data reduction. The SDSS website, http://www.sdss.org, provides links to the data, tutorials and examples of data access, and extensive documentation of the reduction and analysis procedures. DR13 is the first of a scheduled set that will contain new data and analyses from the planned ~6-year operations of SDSS-IV.

POLAR is a new satellite-born detector aiming to measure the polarization of an unprecedented number of Gamma-Ray Bursts in the 50-500 keV energy range. The instrument, launched on-board the Tiangong-2 Chinese Space lab on the 15th of September 2016, is designed to measure the polarization of the hard X-ray flux by measuring the distribution of the azimuthal scattering angles of the incoming photons. A detailed understanding of the polarimeter and specifically of the systematic effects induced by the instrument's non-uniformity are required for this purpose. In order to study the instrument's response to polarization, POLAR underwent a beam test at the European Synchrotron Radiation Facility in France. In this paper both the beam test and the instrument performance will be described. This is followed by an overview of the Monte Carlo simulation tools developed for the instrument. Finally a comparison of the measured and simulated instrument performance will be provided and the instrument response to polarization will be presented.

We describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratio in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spatially-resolved spectroscopy for thousands of nearby galaxies (median redshift of z = 0.03). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between redshifts z = 0.6 and 3.5 to constrain cosmology using baryon acoustic oscillations, redshift space distortions, and the shape of the power spectrum. Within eBOSS, we are conducting two major subprograms: the SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating X-ray AGN and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5-meter Sloan Foundation Telescope at Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5-meter du Pont Telescope at Las Campanas Observatory, with observations beginning in early 2017. Observations at both facilities are scheduled to continue through 2020. In keeping with previous SDSS policy, SDSS-IV provides regularly scheduled public data releases; the first one, Data Release 13, was made available in July 2016.

We compare the structure of molecular gas at $40$ pc resolution to the ability of gas to form stars across the disk of the spiral galaxy M51. We break the PAWS survey into $370$ pc and $1.1$ kpc resolution elements, and within each we estimate the molecular gas depletion time ($\tau_{\rm Dep}^{\rm mol}$), the star formation efficiency per free fall time ($\epsilon_{\rm ff}$), and the mass-weighted cloud-scale (40 pc) properties of the molecular gas: surface density, $\Sigma$, line width, $\sigma$, and $b\equiv\Sigma/\sigma^2\propto\alpha_{\rm vir}^{-1}$, a parameter that traces the boundedness of the gas. We show that the cloud-scale surface density appears to be a reasonable proxy for mean volume density. Applying this, we find a typical star formation efficiency per free-fall time, $\epsilon_{ff} \left( \left< \Sigma_{40pc} \right> \right) \sim 0.3{-}0.36\%$, lower than adopted in many models and found for local clouds. More, the efficiency per free fall time anti-correlates with both $\Sigma$ and $\sigma$, in some tension with turbulent star formation models. The best predictor of the rate of star formation per unit gas mass in our analysis is $b \equiv \Sigma / \sigma^2$, tracing the strength of self gravity, with $\tau_{\rm Dep}^{\rm mol} \propto b^{-0.9}$. The sense of the correlation is that gas with stronger self-gravity (higher $b$) forms stars at a higher rate (low $\tau_{\rm Dep}^{\rm mol}$). The different regions of the galaxy mostly overlap in $\tau_{\rm Dep}^{\rm mol}$ as a function of $b$, so that low $b$ explains the surprisingly high $\tau_{\rm Dep}^{\rm mol}$ found towards the inner spiral arms found by by Meidt et al. (2013).

POLAR is space-borne detector designed for a precise measurement of gamma-ray polarization of the prompt emissions of Gamma-Ray Bursts in the energy range 50 keV - 500 keV. POLAR is a compact Compton polarimeter consisting of 40$\times$ 40 plastic scintillator bars read out by 25 multi-anode PMTs. In May 2015, we performed a series of tests of the POLAR flight model with 100\% polarized x-rays beams at the European Synchrotron Radiation Facility beam-line ID11 aming to study thresholds, crosstalk between channels and responses of POLAR flight model to polarized X-ray beams. In this paper we present the data analysis method and some analysis results. According the results, POLAR FM has good polarimetric capabilities.

As a space-borne detector POLAR is designed to conduct hard X-ray polarization measurements of gamma-ray bursts on the statistically significant sample of events and with an unprecedented accuracy. During its development phase a number of tests, calibrations runs and verification measurements were carried out in order to validate instrument functionality and optimize operational parameters. In this article we present results on gain optimization togeter with verification data obtained in the course of broad laboratory and environmental tests. In particular we focus on exposures to the $^{137}$Cs radioactive source and determination of the gain dependence on the high voltage for all 1600 detection channels of the polarimeter. Performance of the instrument is described in detail with respect to the dynamic range, energy resolution and temperature dependence. Gain optimization algorithms and response non-uniformity studies are also broadly discussed. Results presented below constitute important parts for development of the POLAR calibration and operation database.

In the holographic description of cosmic preheating proposed in an accompanied Letter [1], the energy transfer between the inflaton and matter field at strong coupling is suggested to be mimicked by superfluid and normal components of a superconductor on Friedmann-Robertson-Walker (FRW) boundary in an asymptotically Anti-de Sitter (AdS) spacetime. In this paper we investigated two aspects of the scenario of holographic preheating that are not included in the accompanied work. Firstly, we study in detail the evolution of the quasi-normal modes (QNMs) surrounding a metastable hairy black hole. This analysis can quantitatively describe the preheating process of the matter field that is produced continuously in the case of strong coupling. Secondly, we present a detailed analysis of the holographic renormalization for the AdS-FRW background, which allows us to extract operator expectation values for studying cosmological implications.

We propose a holographic description of cosmic preheating at strong coupling. In this scenario the energy transfer between the inflaton and matter field is mimicked by a model of holographic superconductor. An exponential amplification of the matter field during preheating can be described by the quasi-normal modes of a metastable "black hole" in the bulk spacetime with an expanding boundary. Our results reveal that the matter field can be produced continuously at strong coupling in contrast to the case of weak coupling with a discontinuous matter growth as inflaton oscillates. Furthermore, the amplification of matter field has an enhanced dependence on the vacuum expectation value of the inflaton at strong coupling. By virtue of the proposed mechanism, physics of the very early universe at an extremely high temperature right after inflation may become accessible.

Gamma-ray polarimetry is a new powerful tool to study the processes responsible for the emission from astrophysical sources and the environments in which this emission takes place. Few successful polarimetric measurements have however been performed thus far in the gamma-ray energy band due to the difficulties involved. POLAR is a dedicated polarimeter designed to perform high precision measurements of the polarization of the emission from gamma-ray burst in the 50-500 keV energy range. This new polarimeter is expected to detect approximately 50 gamma-ray bursts per year while performing high precision polarization measurements on approximately 10 bursts per year. The instrument was launched into lower earth orbit as part of the second Chinese space lab, the Tiangong-2, on September 15th 2016 and has been taking data successfully since being switched on one week after. The instrument uses a segmented scintillator array consisting of 1600 plastic scintillator bars, read out by 25 flat-panel multi-anode photomultipliers, to measure the Compton scattering angles of incoming photons. The small segmentation and relatively large uniform effective area allow the instrument to measure the polarization of a large number of transient events, such as gamma-ray bursts, with an unprecedented precision during its two year life-time. The final flight model underwent detailed calibration prior to launch as well as intensive space qualification tests, a summary of which will be presented in this paper. The instrument design will be discussed first followed by an overview of the on-ground tests, finally the in-orbit behavior as measured during the first weeks of the mission will be presented.

Magnetic null has long been recognized as a special structure serving as a preferential site for magnetic reconnection (MR). However, the direct observational study of MR at null-points is largely lacking. Here, we show the observations of MR around a magnetic null associated with an eruption that resulted in an M1.7 flare and a coronal mass ejection. The GOES X- ray profile of the flare exhibited two peaks at 02:23 UT and 02:40 UT on 2012 November 8, respectively. Based on the imaging observations, we find that the first and also primary X- ray peak was originated from MR in the current sheet underneath the erupting magnetic flux rope (MFR). On the other hand, the second and also weaker X-ray peak was caused by MR around a null-point located above the pre-eruption MFR. The interaction of the null-point and the erupting MFR can be described as a two-step process. During the first step, the erupting and fast expanding MFR passed through the null-point, resulting in a significant displacement of the magnetic field surrounding the null. During the second step, the displaced magnetic field started to move back, resulting in a converging inflow and subsequently the MR around the null. The null-point reconnection is a different process from the current sheet reconnection in this flare; the latter is the cause of the main peak of the flare, while the former is the cause of the secondary peak of the flare and the conspicuous high-lying cusp structure.

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.

High-redshift Gamma-Ray Bursts (GRBs) beyond redshift $\sim6$ are potentially powerful tools to probe the distant early Universe. Their detections in large numbers and at truly high redshifts call for the next generation of high-energy wide-field instruments with unprecedented sensitivity at least one order of magnitude higher than the ones currently in orbit. On the other hand, follow-up observations of the afterglows of high-redshift GRBs and identification of their host galaxies, which would be difficult for the currently operating telescopes, require new, extremely large facilities of at multi-wavelengths. This chapter describes future experiments that are expected to advance this exciting field, both being currently built and being proposed. The legacy of Swift will be continued by SVOM, which is equipped with a set of space-based multi-wavelength instruments as well as and a ground segment including a wide angle camera and two follow-up telescopes. The established Lobster-eye X-ray focusing optics provides a promising technology for the detection of faint GRBs at very large distances, based on which the {THESEUS}, {Einstein Probe} and other mission concepts have been proposed. Follow-up observations and exploration of the reionization era will be enabled by large facilities such as {SKA} in the radio, the 30m class telescopes in the optical/near-IR, and the space-borne {WFIRST} and {JWST} in the optical/near-IR/mid-IR. In addition, the X-ray and $\gamma$-ray polarization experiment POLAR is also introduced.

As a part of the Milky Way Imaging Scroll Painting (MWISP) survey, we performed a simultaneous 12CO(1-0), 13CO(1-0), C18O(1-0) mapping toward molecular clouds in a region encompassing 3.75 square degrees. We reveal three molecular clouds, the 15 km/s cloud, the 27 km/s cloud, and the 50 km/s cloud, in the surveyed region. The 50 km/s cloud is resolved with an angular resolution of around 1' for the first time. Investigating their morphology and velocity structures, we find that the 27 km/s cloud is likely affected by feedback from the stellar association Mon OB3 and the 50 km/s cloud is characterised by three large expanding molecular shells. The region is mapped in C18O (1-0) for the first time. We discover seven C18O clumps, which are likely to form massive stars, and 15 dust clumps based on the BGPS archival data. Using infrared color-color diagrams, we find 56 Class I and 107 Class II young stellar object (YSO) candidates. According to the distribution of YSO candidates, an overdensity is found around the HII region S287 and the intersection of two shells, indicative of triggering. The star formation efficiency and rate of the 27~km/s cloud are discussed. Comparing the observed values of the filament S287-main with the models of fragmentation, we suggest that turbulence controls the large scale of fragmentation in the filament while gravitational fragmentation plays an important role in the formation of YSOs on the small scale. We find that star-forming gas tends to have higher excitation temperature, higher 13CO opacity, and higher column density than non-star-forming gas, which is consistent with the point that star formation happens in denser gas and star-forming gas is heated by YSOs. Using the 1.1 mm dust emission to trace dense gas, we obtain a dense gas fraction of 2.7%-10.4% for the 27 km/s cloud.

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.

Gamma Ray Bursts (GRBs) are the strongest explosions in the universe which might be associated with creation of black holes. Magnetic field structure and burst dynamics may influence polarization of the emitted gamma-rays. Precise polarization detection can be an ultimate tool to unveil the true GRB mechanism. POLAR is a space-borne Compton scattering detector for precise measurements of the GRB polarization. It consists of a 40$\times$40 array of plastic scintillator bars read out by 25 multi-anode PMTs (MaPMTs). It is scheduled to be launched into space in 2016 onboard of the Chinese space laboratory TG2. We present a dedicated methodology for POLAR calibration and some calibration results based on the combined use of the laboratory radioactive sources and polarized X-ray beams from the European Synchrotron Radiation Facility. They include calibration of the energy response, computation of the energy conversion factor vs. high voltage as well as determination of the threshold values, crosstalk contributions and polarization modulation factors.

We have performed an unbiased deep near-infrared survey toward the Aquila molecular cloud with a sky coverage of ~1 deg2. We identified 45 molecular hydrogen emission-line objects(MHOs), of which only 11 were previously known. Using the Spitzer archival data we also identified 802 young stellar objects (YSOs) in this region. Based on the morphology and the location of MHOs and YSO candidates, we associate 43 MHOs with 40 YSO candidates. The distribution of jet length shows an exponential decrease in the number of outflows with increasing length and the molecular hydrogen outflows seem to be oriented randomly. Moreover, there is no obvious correlation between jet lengths, jet opening angles, or jet H2 1-0 S(1) luminosities and spectral indices of the possible driving sources in this region. We also suggest that molecular hydrogen outflows in the Aquila molecular cloud are rather weak sources of turbulence, unlikely to generate the observed velocity dispersion in the region of survey.

Magnetic reconnection, a change of magnetic field connectivity, is a fundamental physical process in which magnetic energy is released explosively. It is responsible for various eruptive phenomena in the universe. However, this process is difficult to observe directly. Here, the magnetic topology associated with a solar reconnection event is studied in three dimensions (3D) using the combined perspectives of two spacecraft. The sequence of extreme ultraviolet (EUV) images clearly shows that two groups of oppositely directed and non-coplanar magnetic loops gradually approach each other, forming a separator or quasi-separator and then reconnecting. The plasma near the reconnection site is subsequently heated from $\sim$1 to $\ge$5 MK. Shortly afterwards, warm flare loops ($\sim$3 MK) appear underneath the hot plasma. Other observational signatures of reconnection, including plasma inflows and downflows, are unambiguously revealed and quantitatively measured. These observations provide direct evidence of magnetic reconnection in a 3D configuration and reveal its origin.

The diverse properties of broad-line quasars appear to follow a well-defined main sequence along which the optical FeII strength increases. It has been suggested that this sequence is mainly driven by the Eddington ratio (L/L_Edd) of the black hole (BH) accretion. Shen & Ho demonstrated with quasar clustering analysis that the average BH mass decreases with increasing FeII strength when quasar luminosity is fixed, consistent with this suggestion. Here we perform an independent test by measuring the stellar velocity dispersion sigma* (hence the BH mass via the M-sigma* relation) from decomposed host spectra in low-redshift Sloan Digital Sky Survey quasars. We found that at fixed quasar luminosity, sigma* systematically decreases with increasing FeII strength, confirming that Eddington ratio increases with FeII strength. We also found that at fixed luminosity and FeII strength, there is little dependence of sigma* on the broad Hbeta FWHM. These new results reinforce the framework put forward by Shen & Ho that Eddington ratio and orientation govern most of the diversity seen in broad-line quasar properties.

Type II radio bursts are thought to be a signature of coronal shocks. In this paper, we analyze a short-lived type II burst that started at 07:40 UT on 2011 February 28. By carefully checking white-light images, we find that the type II radio burst is not accompanied by a coronal mass ejection, only with a C2.4 class flare and narrow jet. However, in the extreme-ultraviolet (EUV) images provided by the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory (SDO), we find a wave-like structure that propagated at a speed of $\sim$ 600 km s$^{-1}$ during the burst. The relationship between the type II radio burst and the wave-like structure is in particular explored. For this purpose, we first derive the density distribution under the wave by the differential emission measure (DEM) method, which is used to restrict the empirical density model. We then use the restricted density model to invert the speed of the shock that produces the observed frequency drift rate in the dynamic spectrum. The inverted shock speed is similar to the speed of the wave-like structure. This implies that the wave-like structure is most likely a coronal shock that produces the type II radio burst. We also examine the evolution of the magnetic field in the flare-associated active region and find continuous flux emergence and cancellation taking place near the flare site. Based on these facts, we propose a new mechanism for the formation of the type II radio burst, i.e., the expansion of the strongly-inclined magnetic loops after reconnected with nearby emerging flux acts as a piston to generate the shock wave.

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.

Many synoptic surveys are observing large parts of the sky multiple times. The resulting lightcurves provide a wonderful window to the dynamic nature of the universe. However, there are many significant challenges in analyzing these light curves. These include heterogeneity of the data, irregularly sampled data, missing data, censored data, known but variable measurement errors, and most importantly, the need to classify in astronomical objects in real time using these imperfect light curves. We describe a modeling-based approach using Gaussian process regression for generating critical measures representing features for the classification of such lightcurves. We demonstrate that our approach performs better by comparing it with past methods. Finally, we provide future directions for use in sky-surveys that are getting even bigger by the day.

Magnetic flux rope (MFR) is a coherent and helical magnetic field structure that is recently found probably to appear as an elongated hot-channel prior to a solar eruption. In this paper, we investigate the relationship between the hot-channel and associated prominence through analyzing a limb event on 2011 September 12. In the early rise phase, the hot-channel was cospatial with the prominence initially. It then quickly expanded, resulting in a separation of the top of the hot-channel from that of the prominence. Meanwhile, both of them experienced an instantaneous morphology transformation from a $\Lambda$ shape to a reversed-Y shape and the top of these two structures showed an exponential increase in height. These features are a good indication for the occurrence of the kink instability. Moreover, the onset of the kink instability is found to coincide in time with the impulsive enhancement of the flare emission underneath the hot-channel, suggesting that the ideal kink instability likely also plays an important role in triggering the fast flare reconnection besides initiating the impulsive acceleration of the hot-channel and distorting its morphology. We conclude that the hot-channel is most likely the MFR system and the prominence only corresponds to the cool materials that are collected in the bottom of the helical field lines of the MFR against the gravity.

In the solar corona, magnetic reconnection usually takes place at the singular configuration of magnetic field, in particular near a magnetic null owing to its high susceptibility to perturbations. In this Letter, we report a rare X-shaped structure, encompassing a magnetic null, above a trans-equatorial quadrupole active region group that is well observed by the Atmospheric Imaging Assembly (AIA). The observations show that this X-shaped structure is visible in all AIA EUV passbands and stably exists for days. However, possibly induced by flare activities at the northern part of the quadrupole active region group, the X-shaped structure starts to destabilize and meanwhile a jet erupted near its center at 15:05 UT on 2013 October 7. Through the non-linear force-free field modeling, we identify a magnetic null, which is above the quadrupole polarities and well corresponds to the X-shaped structure. After the jet eruption, the temperature and emission measure of the plasma near the X-shaped structure rise from 2.3 MK and 1.2x10^27 cm-5 at 15:01 UT to 5.4 MK and 3.7x10^27 cm-5 at 15:36 UT, respectively, revealed by the differential emission measure analysis, indicating that magnetic reconnection most likely takes place there to heat the plasma. Moreover, the height of the null has an increase of 10 Mm, which is most likely due to the partial opening of the field lines near the fan surface that makes the null underneath rise to seek for a new equilibrium.

We perform Differential Emission Measure (DEM) analysis of an M7.7 flare that occurred on 2012 July 19 and was well observed by the Atmospheric Imaging Assembly (AIA) aboard the Solar Dynamic Observatory. Using the observational data with unprecedented high temporal and spatial resolution from six AIA coronal passbands, we calculate the DEM of the flare and derive the time series of maps of DEM-weighted temperature and emission measure (EM). It is found that, during the flare, the highest EM region is located in the flare loop top with a value varying between 8.4 10^28 cm^-5 and 2.5 10^30 cm^-5. The temperature there rises from 8 MK at about 04:40 UT (the initial rise phase) to a maximum value of 13 MK at about 05:20 UT (the hard X-ray peak). Moreover, we find a hot region that is above the flare loop top with a temperature even up to 16 MK. We also analyze the DEM properties of the reconnection site. The temperature and density there are not as high as that in the loop top and the flux rope, indicating that the main heating may not take place inside the reconnection site. In the end, we examine the dynamic behavior of the flare loops. Along the flare loop, both the temperature and the EM are the highest in the loop top and gradually decrease towards the footpoints. In the northern footpoint, an upward force appears with a biggest value in the impulsive phase, which we conjecture originates from chromospheric evaporation.

The phase diagram and equation of state of dense nitrogen are of interest in understanding the fundamental physics and chemistry under extreme conditions, including planetary processes, and in discovering new materials. We predict several stable phases of nitrogen at multi-TPa pressures, including a P4/nbm structure consisting of partially charged N2 pairs and N5 tetrahedra, which is stable in the range 2.5-6.8 TPa. This is followed by a modulated layered structure between 6.8 and 12.6 TPa, which also exhibits a significant charge transfer. The P4/nbm metallic nitrogen salt and the modulated structure are stable at high pressures and temperatures, and they exhibit strongly ionic features and charge density distortions, which is unexpected in an element under such extreme conditions and could represent a new class of nitrogen materials. The P-T phase diagram of nitrogen at TPa pressures is investigated using quasiharmonic phonon calculations and ab initio molecular dynamics simulations.

The magnetic flux rope (MFR) is believed to be the underlying magnetic structure of coronal mass ejections (CMEs). However, it remains unclear how an MFR evolves into and forms the multi-component structure of a CME. In this paper, we perform a comprehensive study of an extreme-ultraviolet (EUV) MFR eruption on 2013 May 22 by tracking its morphological evolution, studying its kinematics, and quantifying its thermal property. As EUV brightenings begin, the MFR starts to rise slowly and shows helical threads winding around an axis. Meanwhile, cool filamentary materials descend spirally down to the chromosphere. These features provide direct observational evidence of intrinsically helical structure of the MFR. Through detailed kinematical analysis, we find that the MFR evolution experiences two distinct phases: a slow rise phase and an impulsive acceleration phase. We attribute the first phase to the magnetic reconnection within the quasi-separatrix-layers surrounding the MFR, and the much more energetic second phase to the fast magnetic reconnection underneath the MFR. We suggest that the transition between these two phases be caused by the torus instability. Moreover, we identify that the MFR evolves smoothly into the outer corona and appears as a coherent structure within the white light CME volume. The MFR in the outer corona was enveloped by bright fronts that originated from plasma pile-up in front of the expanding MFR. The fronts are also associated with the preceding sheath region followed the outmost MFR-driven shock.

We investigate the acceleration source of the impulsive solar energetic particle (SEP) events on 2007 January 24. Combining the in situ electron measurements and remote-sensing solar observations, as well as the calculated magnetic fields obtained from a potential-field source-surface model, we demonstrate that the jets associated with the hard X-ray flares and type-III radio bursts, rather than the slow and partial coronal mass ejections, are closely related to the production of interplanetary electron streams. The jets, originated from the well-connected active region (AR 10939) whose magnetic polarity structure favors the eruption, are observed to be forming in a coronal site, extending to a few solar radii, and having a good temporal correlation with the electron solar release. The open-field lines near the jet site are rooted in a negative polarity, along which energetic particles escape from the flaring AR to the near-Earth space, consistent with the in situ electron pitch angle distribution. The analysis enables us to propose a coronal magnetic topology relating the impulsive SEP events to their solar source.

Outflows of plasma at the edges of active regions surrounded by quiet Sun are now a common observation with the Hinode satellite. While there is observational evidence to suggest that the outflows are originating in the magnetic field surrounding the active regions, there is no conclusive evidence that reveals how they are driven. Motivated by observations of outflows at the periphery of a mature active region embedded in a coronal hole, we have used a three-dimensional simulation to emulate the active region's development in order to investigate the origin and driver of these outflows. We find outflows are accelerated from a site in the coronal hole magnetic field immediately surrounding the active region and are channelled along the coronal hole field as they rise through the atmosphere. The plasma is accelerated simply as a result of the active region expanding horizontally as it develops. Many of the characteristics of the outflows generated in the simulation are consistent with those of observed outflows: velocities up to 45 km per sec, properties akin to the coronal hole, proximity to the active region's draining loops, expansion with height, and projection over monopolar photospheric magnetic concentrations. Although the horizontal expansion occurs as a consequence of the active region's development in the simulation, expansion is also a general feature of established active regions. Hence, it is entirely possible and plausible that the expansion acceleration mechanism displayed in the simulation is occurring in active regions on the Sun and, in addition to reconnection, is driving the outflows observed at their edges.

Using a sample of 248 metal-poor stars (RR Lyraes, red giants and RHB stars) which is remarkable for the accuracy of its 6-D kinematical data, we find a new component for the local halo which has an axial ratio c/a ~ 0.2, a similar flattening to the thick disk. It has a small prograde rotation but is supported by velocity anisotropy, and contains more intermediate-metallicity stars (with -1.5 < [Fe/H] < -1.0) than the rest of our sample. We suggest that this component was formed quite late, during or after the formation of the disk. It formed either from the gas that was accreted by the last major mergers experienced by the Galaxy, or by dynamical friction of massive infalling satellite(s) with the halo and possibly the stellar disk or thick disk. The remainder of the stars in our sample exhibit a clumpy distribution in energy and angular momentum, suggesting that the early, chaotic conditions under which the inner halo formed were not violent enough to erase the record of their origins. The clumpy structure suggests that a relatively small number of progenitors were responsible for building up the inner halo, in line with theoretical expectations. We find a difference in mean binding energy between the RR Lyrae variables and the red giants in our sample, suggesting that more of the RR Lyraes in the sample belong to the outer halo, and that the outer halo may be somewhat younger, as first suggested by Searle and Zinn (1978). We also find that the RR Lyrae mean rotation is more negative than the red giants, which is consistent with the recent result of Carollo et al.(2007) that the outer halo has a retrograde rotation and with the difference in kinematics seen between RR Lyraes and BHB stars by Kinman et al.(2007).

We develop a mixture-based approach to robust density modeling and outlier detection for experimental multivariate data that includes measurement error information. Our model is designed to infer atypical measurements that are not due to errors, aiming to retrieve potentially interesting peculiar objects. Since exact inference is not possible in this model, we develop a tree-structured variational EM solution. This compares favorably against a fully factorial approximation scheme, approaching the accuracy of a Markov-Chain-EM, while maintaining computational simplicity. We demonstrate the benefits of including measurement errors in the model, in terms of improved outlier detection rates in varying measurement uncertainty conditions. We then use this approach in detecting peculiar quasars from an astrophysical survey, given photometric measurements with errors.

Parametric Embedding (PE) has recently been proposed as a general-purpose algorithm for class visualisation. It takes class posteriors produced by a mixture-based clustering algorithm and projects them in 2D for visualisation. However, although this fully modularised combination of objectives (clustering and projection) is attractive for its conceptual simplicity, in the case of high dimensional data, we show that a more optimal combination of these objectives can be achieved by integrating them both into a consistent probabilistic model. In this way, the projection step will fulfil a role of regularisation, guarding against the curse of dimensionality. As a result, the tradeoff between clustering and visualisation turns out to enhance the predictive abilities of the overall model. We present results on both synthetic data and two real-world high-dimensional data sets: observed spectra of early-type galaxies and gene expression arrays.

We present a study of molecular outflows using high-resolution mapping of the CO (1-0) line emission toward eight relatively nearby 6.7 GHz methanol masers, which are associated with massive star forming regions. Outflows were detected in seven out of eight sources, and five of them clearly show bipolar or multiple outflow morphologies. These outflows have typical masses of a few solar masses, momenta of tens of $M_{\odot}$ km s$^{-1}$, kinetic energies of $\sim$ 10$^{45}$ ergs, and the mass entrainment rate of a few $10^{-5} M_{\odot}$ yr$^{-1}$. They have significantly more mass and kinetic energy than their low-mass counterparts. In some of the sources, the massive outflow is obviously associated with a particular massive star in the cluster, while in others the origin remains uncertain. The high detection rate of outflows toward methanol masers suggests that the outflow phase of massive protostars encompasses the methanol maser phase.