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.

Matching the number counts of high-$z$ sub-millimetre-selected galaxies (SMGs) has been a long standing problem for galaxy formation models. In this paper, we use 3D dust radiative transfer to model the sub-mm emission from galaxies in the SIMBA cosmological hydrodynamic simulations, and compare predictions to the latest single-dish observational constraints on the abundance of 850$\mathrm{\mu m}$-selected sources. We find unprecedented agreement with the integrated 850$\mathrm{\mu m}$ luminosity function, along with good agreement in the redshift distribution of bright SMGs. The excellent agreement is driven primarily by SIMBA's good match to infrared measures of the star formation rate (SFR) function between $z = 2-4$ at high SFRs. Also important is the self-consistent on-the-fly dust model in SIMBA, which predicts, on average, higher dust masses (by up to a factor of 7) compared to using a fixed dust-to-metals ratio of 0.3. We construct a lightcone to investigate the effect of far-field blending, and find minimal contribution to the shape and normalisation of the luminosity function. We provide new fits to the 850$\mathrm{\mu m}$ luminosity as a function of SFR and dust mass. Our results demonstrate that exotic solutions to the discrepancy between sub-mm counts in simulations and observations, such as a top-heavy IMF, are unnecessary, and that sub-millimetre-bright phases are a natural consequence of massive galaxy evolution.

The authors of this white paper met on 16-17 January 2020 at the New Jersey Institute of Technology, Newark, NJ, for a 2-day workshop that brought together a group of heliophysicists, data providers, expert modelers, and computer/data scientists. Their objective was to discuss critical developments and prospects of the application of machine and/or deep learning techniques for data analysis, modeling and forecasting in Heliophysics, and to shape a strategy for further developments in the field. The workshop combined a set of plenary sessions featuring invited introductory talks interleaved with a set of open discussion sessions. The outcome of the discussion is encapsulated in this white paper that also features a top-level list of recommendations agreed by participants.

We present Powderday, a flexible, fast, open-source dust radiative transfer package designed to interface with galaxy formation simulations. Powderday builds on FSPS population synthesis models, Hyperion dust radiative transfer, and employs yt to interface between different software packages. We include our stellar population synthesis modeling on the fly, which allows for significant run-time flexibility in the assumed stellar physics. We include a model for nebular line emission that can employ either precomputed Cloudy lookup tables (for efficiency), or direct photoionization calculations for all young stars (for flexibility). The dust content follows either observationally-motivated prescriptions, direct modeling from galaxy formation simulations, or a novel approach that includes the dust content via learning-based algorithms from the SIMBA cosmological galaxy formation simulation. AGN can additionally be included via a range of prescriptions. The output of these models are broadband SEDs, as well as filter-convolved images. Powderday is designed to eliminate last-mile efforts by researchers that employ different hydrodynamic galaxy formation models, and seamlessly interfaces with GIZMO, AREPO, GASOLINE, CHANGA, and ENZO. We demonstrate the capabilities of the code via three applications: a model for the star formation rate (SFR) - infrared luminosity relation in galaxies (including the impact of AGN); the impact of circumstellar dust around AGB stars on the mid-infrared emission from galaxy SEDs; and the impact of galaxy inclination angle on dust attenuation laws.

Combining the LAMOST radial velocities with Gaia parallaxes and proper motions, we presented 3D Galactic space motions and the orbits of 182 single-lined hot subdwarf stars. These stars have been identified by Lei et al. (2020) in Gaia DR2 with LAMOST DR6 and DR7 spectra. He-rich hot subdwarf stars with log(y)>0 show the largest standard deviations of the Galactic velocity components and orbital parameters, while those with -1<log(y)<0 exhibit the second largest standard deviations. The two groups of He-deficient stars with log(y)<-1 show similar standard deviations, which is systematically lower compared to He-rich stars. We also presented a kinematic population classification of the four hot subdwarf helium groups based on their positions in the U-V velocity diagram, J_z-eccentricity diagram and their Galactic orbits. The overall tendency of the fractional distributions of the four hot subdwarf helium groups in the halo, thin disk and thick disk is largely consistent with the findings reported by Luo et al.(2019) based on LAMOST DR5, which appears to support the predictions of binary population synthesis (Han et al. 2003; 2008). He-deficient stars with -2.2<log(y)<-1 likely origin from stable the Roche lobe overflow channel, He-deficient stars with log(y)<-2.2 from the common envelope ejection channel, and He-rich stars with log(y)>0 from the merger channel of double He white dwarf stars. The fraction of He-rich hot subdwarf stars with -1<log(y)<0 in the thin disk and the halo are far higher than in the thick disk, which implies that these stars have different formation channels in the thin disk and in the halo.

Using the catalogues of galaxy clusters from The Three Hundred project, modelled with both hydrodynamic simulations, (Gadget-X and Gadget-MUSIC), and semi-analytic models (SAMs), we study the scatter and self-similarity of the profiles and distributions of the baryonic components of the clusters: the stellar and gas mass, metallicity, the stellar age, gas temperature, and the (specific) star formation rate. Through comparisons with observational results, we find that the shape and the scatter of the gas density profiles matches well the observed trends including the reduced scatter at large radii which is a signature of self-similarity suggested in previous studies. One of our simulated sets, Gadget-X, reproduces well the shape of the observed temperature profile, while Gadget-MUSIC has a higher and flatter profile in the cluster centre and a lower and steeper profile at large radii. The gas metallicity profiles from both simulation sets, despite following the observed trend, have a relatively lower normalisation. The cumulative stellar density profiles from SAMs are in better agreement with the observed result than both hydrodynamic simulations which show relatively higher profiles. The scatter in these physical profiles, especially in the cluster centre region, shows a dependence on the cluster dynamical state and on the cool-core/non-cool-core dichotomy. The stellar age, metallicity and (s)SFR show very large scatter, which are then presented in 2D maps. We also do not find any clear radial dependence of these properties. However, the brightest central galaxies have distinguishable features compared to the properties of the satellite galaxies.

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.

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.

The first repeating fast radio burst (FRB), FRB 121102, was found to be associated with a spatially coincident, persistent nonthermal radio source, but the origin of the persistent emission remains unknown. In this paper, we propose that the persistent emission is produced via synchrotron-heating process by multiple bursts of FRB 121102 in a self-absorbed synchrotron nebula. As a population of bursts of the repeating FRB absorbed by the synchrotron nebula, the energy distribution of electrons in the nebula will change significantly. As a result, the spectrum of the nebula will show a hump steadily. For the persistent emission of FRB 121102, the total energy of bursts injecting into the nebula is required to be about $3.3\times10^{49}\,\unit{erg}$, the burst injection age is over $6.7\times 10^4\,\unit{yr}$, the nebula size is $\sim0.02\,\unit{pc}$, and the electron number is about $3.2\times10^{55}$. We predict that as more bursts inject, the brightness of the nebula would be brighter than the current observation, and meanwhile, the peak frequency would become higher. Due to the synchrotron absorption of the nebula, some low-frequency bursts would be absorbed, which may explain why most bursts were detected above $\sim1~\unit{GHz}$.

With the launch of SDO, many EUV waves were observed during solar eruptions. However, the joint observations of Moreton and EUV waves are still relatively rare. We present an event that a small-scale filament eruption simultaneously results in a Moreton wave, an EUV wave and a Coronal Mass Ejection in active region NOAA 12740. Firstly, we find that some dark elongate lanes or filamentary structures in the photosphere existed under the small-scale filament and drifted downward, which manifests that the small-scale filament was emerging and lifting from subsurface. Secondly, combining the simultaneous observations in different Extreme UltraViolet (EUV) and H$\alpha$ passbands, we study the kinematic characteristics of the Moreton and EUV waves. The comparable propagating velocities and the similar morphology of Moreton and different passbands EUV wavefronts were obtained. We deduce that Moreton and different passbands EUV waves were the perturbations in different temperature-associated layers induced by the coronal magneto-hydrodynamic shock wave. We also find the refracted, reflected and diffracted phenomena during the propagation of the EUV wave. By using power-law fittings, the kinematic characteristics of unaffected, refracted and diffracted waves were obtained. The extrapolation field derived by the potential field source surface (PFSS) model manifests that the existence of an interface of different magnetic system (magnetic separatrix) result in refraction, reflection and deviation of the EUV wave.

The physical origin of fast radio bursts (FRBs) is still unknown. Multiwavelength and polarization observations of an FRB source would be helpful to diagnose its progenitor and environment. So far only the first repeating source FRB 121102 appears to be spatially coincident with a persistent radio emission. Its bursts also have very large values of the Faraday rotation measure (RM) i.e., $|{\rm RM}|\sim10^5~{\rm rad~m^{-2}}$. We show that theoretically there should be a simple relation between RM and the luminosity of the persistent source of an FRB source if the observed RM mostly arises from the persistent emission region. FRB 121102 follows this relation given that the magnetic field in the persistent emission region is highly ordered and that the number of relativistic electrons powering the persistent emission is comparable to that of non-relativistic electrons that contribute to RM. The non-detections of persistent emission sources from all other localized FRB sources are consistent with their relatively small RMs ($\left|{\rm RM}\right|\lesssim{\rm a~few}\times100~{\rm rad~m^{-2}}$) according to this relation. Based on this picture, the majority of FRBs without a large RM are not supposed to be associated with bright persistent sources.

The effective vacuum energy density contributed by the non-trivial contortion distribution and the bare vacuum energy density can be viewed as the energy density of the auxiliary quintessence field potential. We find that the negative bare vacuum energy density from string landscape leads to a monotonically decreasing quintessence potential while the positive one from swampland leads to the meta stable or stable de Sitter like potential. Moreover, the non-trivial Brans-Dicke like coupling between quintessence field and gravitation field is necessary in the latter case.

We report the discovery of a binary millisecond pulsar (namely PSR J1641+3627F or M13F) in the globular cluster M13 (NGC 6205) and timing solutions of M13A to F using observations made with the Five-hundred-metre Aperture Spherical radio Telescope (FAST). PSR J1641+3627F has a spin period of 3.00 ms and an orbital period of 1.4 days. The most likely companion mass is 0.16 M$_{\odot}$. M13A to E all have short spin periods and small period derivatives. We also confirm that the binary millisecond pulsar PSR J1641$+$3627E (also M13E) is a black widow with a companion mass around 0.02 M$_{\odot}$. We find that all the binary systems have low eccentricities compared to those typical for globular cluster pulsars and that they decrease with distance from the cluster core. This is consistent with what is expected as this cluster has a very low encounter rate per binary.

Since their discovery in 1991, carbon nanotubes (CNTs) -- a novel one-dimensional carbon allotrope -- have attracted considerable interest worldwide because of their potential technological applications such as electric and optical devices. In the astrophysical context, CNTs may be present in the interstellar space since many of the other allotropes of carbon (e.g., amorphous carbon, fullerenes, nanodiamonds, graphite, polycyclic aromatic hydrocarbons, and possibly graphene as well) are known to be widespread in the Universe, as revealed by presolar grains in carbonaceous primitive meteorites and/or by their fingerprint spectral features in astronomical spectra. In addition, there are also experimental and theoretical pathways to the formation of CNTs in the interstellar medium (ISM). In this work, we examine their possible presence in the ISM by comparing the observed interstellar extinction curve with the ultraviolet/optical absorption spectra experimentally obtained for single-walled CNTs of a wide range of diameters and chiralities. Based on the absence in the interstellar extinction curve of the ~4.5 and 5.25 eV $\pi$-plasmon absorption bands which are pronounced in the experimental spectra of CNTs, we place an upper limit of ~10 ppm of C/H (i.e., ~4% of the total interstellar C) on the interstellar CNT abundance.

The MAMMOTH-I Nebula at redshift 2.3 is one of the largest known Ly-alpha nebulae in the Universe, spanning ~440 kpc. Enormous Ly-alpha nebulae like MAMMOTH-I typically trace the densest and most active regions of galaxy formation. Using sensitive low-surface-brightness observations of CO(1-0) with the Very Large Array, we trace the cold molecular gas in the inner 150 kpc of the MAMMOTH-I Nebula. CO is found in four regions that are associated with either galaxies or groups of galaxies that lie inside the nebula. In three of the regions, the CO stretches up to ~30 kpc into the circum-galactic medium (CGM). In the centermost region, the CO has a very low velocity dispersion (FWHM$_{\rm CO}$ ~ 85 km/s), indicating that this gas is dynamically cold. This dynamically cold gas coincides with diffuse restframe optical light in the CGM around a central group of galaxies, as discovered with the Hubble Space Telescope. We argue that this likely represents cooling of settled and enriched gas in the center of MAMMOTH-I. This implies that the dynamically cold gas in the CGM, rather than the obscured AGN, marks the core of the potential well of this Ly-alpha nebula. In total, the CO in the MAMMOTH-I Nebula traces a molecular gas mass of M$_{\rm H2}$ ~ 1.4 ($\alpha_{\rm CO}$/3.6) $\times$ 10$^{11}$ M$_{\odot}$, with roughly 50% of the CO(1-0) emission found in the CGM. Our results add to the increasing evidence that extended reservoirs of molecular gas exist in the CGM of massive high-z galaxies and proto-clusters.

To better understand the formation and decay of sunspot penumbra, we studied the evolution of sunspots in the three regions of the active-region NOAA 12673 in detail. The evolution of sunspots in the three regions were all involved in the interaction of two magnetic field systems: the pre-existing magnetic field system and the later emerging magnetic field system. Through analyzing the photospheric magnetic field properties, it is found that the formation of the penumbra originated from the newly emerging magnetic bipole that were trapped in the photosphere. The change of magnetic fields in penumbra from horizontal to vertical can cause the disappearance of penumbra. The transformation of the magnetic field between the umbra and the penumbra is found and the outward moat flow around sunspot gradually decreased and vanished during the decay of sunspot. In addition, we found that the mean longitudinal magnetic strength in penumbra decreased and the mean transverse magnetic strength in penumbra increased with the increasing penumbral area during the formation of sunspots. However, during the decay of sunspots, the mean longitudinal magnetic strength in penumbra increased and the mean transverse magnetic strength in penumbra decreased with the decreasing penumbral area. Comparatively, the dependence of the area and the mean transverse/longitudinal magnetic field strength in umbra is not remarkable. These results reveal that the formation and decay process of umbra are different with penumbra.

The massive galaxy population above the characteristic Schechter mass M*~10^10.6 Msun contributes to about half of the total stellar mass in the local Universe. These massive galaxies usually reside in hot dark matter haloes above the critical shock-heating mass ~10^12 Msun, where the external cold gas supply to these galaxies is expected to be suppressed. When galaxies run out of their cold gas reservoir, they become dead and quiescent. Therefore, massive quiescent galaxies living in hot haloes are commonly believed to be gas-poor. Based on the data from SDSS, ALFALFA, GASS and COLD GASS surveys, here we show that the vast majority of the massive, quiescent, central disk galaxies in the nearby Universe have a remarkably large amount of cold atomic hydrogen gas, surprisingly similar to star-forming galaxies. Both star-forming and quiescent disk galaxies show identical symmetric double-horn HI spectra, indicating similar regularly rotating HI disks. Relative to their star-forming counterparts, massive quiescent central disk galaxies are quenched because of their significantly reduced molecular gas content, lower dust content, and lower star formation efficiency. Our findings reveal a new picture, which clearly demonstrates the detailed star-formation quenching process in massive galaxies and provides a stringent constraint on the physical mechanism of quenching.

The possible presence of graphene in the interstellar medium (ISM) is examined by comparing the interstellar extinction curve with the ultraviolet absorption of graphene calculated from its dielectric functions experimentally obtained with the electron energy loss spectroscopy (EELS) method. Based on the absence in the interstellar extinction curve of the ~2755 Angstrom $\pi$--$\pi^{\ast}$ electronic interband transition of graphene, we place an upper limit of ~20ppm of C/H on the interstellar graphene abundance, exceeding the previous estimate by a factor of ~3 which made use of the dielectric functions measured with the spectroscopic ellipsometry (SE) method. Compared with the SE method which measures graphene in air (and hence its surface is contaminated) in a limited energy range of ~0.7--5 eV, the EELS probes a much wider energy range of ~0--50 eV and is free of contamination. The fact that the EELS dielectric functions are substantially smaller than that of SE naturally explains why a higher upper limit on the graphene abundance is derived with EELS. Inspired by the possible detection of C$_{24}$, a planar graphene sheet, in several Galactic and extragalactic planetary nebulae, we also examine the possible presence of C$_{24}$ in the diffuse ISM by comparing the model IR emission of C$_{24}$ with the observed IR emission of the Galactic cirrus and the diffuse ISM toward $l = 44^{o}20^\prime$ and $b=-0^{o}20^\prime$. An upper limit of ~20 ppm on C$_{24}$ is also derived from the absence of the characteristic vibrational bands of C$_{24}$ at ~6.6, 9.8 and 20 $\mu$m in the observed IR emission.

Motivated by the recent discovery of the near-ubiquity of Ly$\alpha$ emission around $z \gtrsim 3$ QSOs, we performed a systematic study of QSO circumgalactic Ly$\alpha$ emission at $z\approx2$, utilizing the unique capability of the Keck Cosmic Web Imager (KCWI) -- a new wide-field, blue sensitive integral-field spectrograph (IFU). In this paper, we present KCWI observations on a sample of 16 ultraluminous Type-I QSOs at $z=2.1-2.3$ with ionizing luminosities of $L_{\rm{\nu_{\rm{LL}}}}=10^{31.1-32.3}$ erg s$^{-1}$ Hz$^{-1}$. We found that 14 out of 16 QSOs are associated with Ly$\alpha$ nebulae with projected linear-sizes larger than 50 physical kpc (pkpc). Among them, four nebulae have enormous Ly$\alpha$ emission with the Ly$\alpha$ surface brightness $SB_{\rm{Ly\alpha}} >10^{-17}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$ on the $>100$ kpc scale, extending beyond the field of view of KCWI. Our KCWI observations reveal that most $z\approx2$ QSO nebulae have a more irregular morphology compared to those at $z\gtrsim3$. In turn, we measure that the circularly-averaged surface brightness (SB) at $z\approx2$ is 0.4 dex fainter than the redshift-corrected, median SB at $z\gtrsim 3$. The Ly$\alpha$ SB profile (SB$_{\rm{Ly\alpha}}$) of QSOs at $z\approx 2$ can be described by a power law of SB$_{\rm{Ly\alpha},z\approx2.3}= 3.7\times10^{-17}\times(r/40)^{-1.8}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$, with the slope similar to that at $z\gtrsim3$. The observed lower redshift-corrected, circularly-averaged SB may be mainly due to the lower covering factor of cool gas clouds in massive halos at $z\approx2$.

We present predictions for the evolution of the galaxy dust-to-gas (DGR) and dust-to-metal (DTM) ratios from z=0 to 6, using a model for the production, growth, and destruction of dust grains implemented into the \simba\ cosmological hydrodynamic galaxy formation simulation. In our model, dust forms in stellar ejecta, grows by the accretion of metals, and is destroyed by thermal sputtering and supernovae. Our simulation reproduces the observed dust mass function at z=0, but modestly under-predicts the mass function by ~x3 at z ~ 1-2. The z=0 DGR vs metallicity relationship shows a tight positive correlation for star-forming galaxies, while it is uncorrelated for quenched systems. There is little evolution in the DGR-metallicity relationship between z=0-6. We use machine learning techniques to search for the galaxy physical properties that best correlate with the DGR and DTM. We find that the DGR is primarily correlated with the gas-phase metallicity, though correlations with the depletion timescale, stellar mass and gas fraction are non-negligible. We provide a crude fitting relationship for DGR and DTM vs. the gas-phase metallicity, along with a public code package that estimates the DGR and DTM given a set of galaxy physical properties.

The Solar System harbors deep unresolved mysteries despite centuries-long study. A highly intriguing case concerns anomalous non-dipolar and non-axisymmetric magnetic fields of Uranus and Neptune that have long eluded explanation by the prevailing theory. A thin-shell dynamo conjecture captures observed phenomena but leaves unexplained fundamental material basis and underlying mechanism. Here, we report the discovery of trihydrogen oxide (H3O) in metallic liquid state stabilized at extreme pressure and temperature conditions inside these icy planets. Calculated stability pressure field compared to known pressure-radius relation for Uranus and Neptune places metallic liquid H3O in a thin-shell zone near planetary cores. These findings from accurate quantum mechanical calculations rationalize the empirically conjectured thin-shell dynamo model and establish key physical benchmarks that are essential to elucidating the enigmatic magnetic-field anomaly of Uranus and Neptune, resolving a major mystery in planetary science.

We selected 90 massive star-forming clumps with strong N2H+, HCO+, HCN, and HNC emission from the Millimetre Astronomy Legacy Team 90 GHz survey. We obtained Herschel data for all 90 sources and NRAO VLA Sky Survey data for 51 of them. We convolved and regridded all images to the same resolution and pixel size and derived the temperature, H2 column density, molecules' abundances and abundance, and ratios of each pixel. Our analysis yields three main conclusions. First, the abundances of N2H+, HCO+, HCN, and HNC increase when the column density decreases and the temperature increases, with spatial variations in their abundances dominated by changes in the H2 column density. Second, the abundance ratios between N2H+, HCO+, HCN, and HNC also display systemic variations as a function of the column density due to the chemical properties of these molecules. Third, the sources associated with the 20 cm continuum emission can be classified into four types based on the behavior of the abundances of the four molecules considered here as a function of this emission. The variations of the first three types could also be attributed to the variation of the H2 column density.

A total of 188 high-mass outflows have been identified from a sample of 694 clumps from the Millimetre Astronomy Legacy Team 90 GHz survey, representing a detection rate of approximately 27%. The detection rate of outflows increases from the protostellar stage to the H II stage, but decreases again at the photodissociation (PDR) stage suggesting that outflows are being switched off during the PDR stage. An intimate relationship is found between outflow action and the presence of masers, and water masers appear together with 6.7 GHz methanol masers. Comparing the infall detection rate of clumps with and without outflows, we find that outflow candidates have a lower infall detection rate. Finally, we find that outflow action has some influence on the local environment and the clump itself, and this influence decreases with increasing evolutionary time as the outflow action ceases.

A hot corona is suggested to be above the standard thin disk. The anisotropy of hard X-ray emission in radio-quiet active galactic nuclei implies that the corona is not static and probably moves outwards like winds. We perform two-dimensional magnetohydrodynamical simulations to study the outflowing corona driven by magnetic field and radiation force. In our simulations, as the initial state and the boundary condition at the disk surface, the corona temperature is set to $10^9$ K inside 10 Schwarzschild radius ($r_{\rm s}$) while the corona temperature is set to $10^7$ K a 10 $r_{\rm s}$. We employ a weak poloidal magnetic field as the initial magnetic field. A collimated outflow and a wide-angle ordered outflow are observed in our simulations. The collimated outflow is around the rotational axis and has a bulk velocity of $\sim$0.03--0.3 $c$ ( $c$ is speed of light) at 90 $r_{\rm s}$, while their mass outflow rate is very low. The collimated outflow is a weak jet. The wide-angle ordered outflow is distributed at the middle and high latitudes and moves outwards with a velocity of $10^2$--$10^4$ km s$^{-1}$. The outflow velocity depends on the disk luminosity. The gas around the disk surface is turbulent, especially outside of 10 $r_{\rm s}$. The other properties of outflows are discussed in detail.

With the observations of SDO/AIA 304 {\AA} and NVST H{\alpha} bands, we present the formation process of an active-region filament in active region NOAA 11903 during the period from 02:00 UT to 10:00 UT on November 25, 2013. A series of jets occurring in the vicinity of the south-western footpoint of the filament directly ejected cool and hot plasma to filament height and supplied material for the filament. Some newly emerging flux is found in the vicinity of the filament south-western footpoint during these jets. In this paper, we mainly focus on the material supply of the filament. The plasma mass uploaded by the jets and the mass of the filament are estimated, which manifests the fact that the mass carried by jets can supply sufficient material for the formation of the filament. We found two types of jets. one is H-alpha jet, and the other is EUV jet. The significant finding is that some cool jets seen in the H{\alpha} band but not in the SDO/AIA bands also could eject the cool material for the filament. These results suggest that cool plasma in the low atmosphere can be directly injected into the upper atmosphere and become the filament material by two types of jets. Moreover, the newly emerging flux with non-potential field plays an important role in the appearances of the jets and the magnetic structure of the filament.

We show quantitatively that an exact solution of Einstein's conformal gravity can explain very well the galactic rotation curves for a sample of 104 galaxies without need for dark matter or other exotic modification of gravity. The metric is an overall rescaling of the Schwarzschild-de Sitter spacetime as required by Weyl conformal invariance, which has to be spontaneously broken, and the velocity of the stars depends only on two universal parameters determined on the base of the observational data.

A multi-wavelength analysis of the large Galactic infrared bubble N 24 is has been presented in this paper in order to investigate the molecular and star formation environment around expanding HII regions. Using archival data from Herschel and ATLASGAL, the distribution and physical properties of the dust over the entire bubble are studied. Twenty three dense clumps are identified using the Clumpfind2d algorithm with sizes and masses in the range 0.65-1.73 pc and 600-16300 Msun, respectively. To analyse the molecular environment in N 24, observations of NH3 (1,1) and (2,2) were carried out using the Nanshan 26m radio telescope. Analysis of the kinetic temperature and gravitational stability of these clumps suggests gravitational collapse in several of them. The mass-size distributions of the clumps and the presence of massive young protostars indicate that the shell of N 24 is a region of ongoing massive star formation. The compatibility of the dynamical and fragmentation time-scales and the overabundance of YSOs and clumps on the rim suggest that the collect and collapse mechanism is in play at the boundary of the bubble, but the existence of the IRDC at the edge of bubble indicates that radiation-driven implosion mechanism may also have played a role there.

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 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 introduce the Simba simulations, the next generation of the Mufasa cosmological galaxy formation simulations run with Gizmo's meshless finite mass hydrodynamics. Simba includes updates to Mufasa's sub-resolution star formation and feedback prescriptions, and introduces black hole growth via the torque-limited accretion model of Angl\'es-Alc\'azar et al. (2017) from cold gas and Bondi accretion from hot gas, along with black hole feedback via kinetic bipolar outflows and X-ray energy. Ejection velocities are taken to be ~10^3 km/s at high Eddington ratios, increasing to ~8000 km/s at Eddington ratios below 2%, with a constant momentum input of 20L/c. Simba further includes an on-the-fly dust production, growth, and destruction model. Our Simba run with (100 Mpc/h)^3 and 1024^3 gas elements reproduces numerous observables, including galaxy stellar mass functions at z=0-6, the stellar mass--star formation rate main sequence, HI and H2 fractions, the mass-metallicity relation at z=0 and z=2, star-forming galaxy sizes, hot gas fractions in massive halos, and z=0 galaxy dust properties. However, Simba also yields an insufficiently sharp truncation of the z=0 mass function, and too-large sizes for low-mass quenched galaxies. We show that Simba's jet feedback is primarily responsible for quenching massive galaxies.

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 present Keck/KCWI integral field spectrograph observations of the complex system surrounding SMM J02399$-$0136 (a lensed $z=2.8$ sub-mm galaxy), including an associated Ly$\rm \alpha$ nebula, a dust-obscured, broad-absorption-line quasar, and neighboring galaxies. At a 3$\sigma$ surface brightness contour of 1.6$\times$ 10$^{-17}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$, the Ly$\rm \alpha$ nebula extends over 17 arcsec ($\gtrsim 140$ physical kpc) and has a total Ly$\rm \alpha$ luminosity of $2.5 \times 10^{44} \, \rm erg \, s^{-1}$ (uncorrected for lensing). The nebula exhibits a kinematic shear of $\sim$ 1000 km s$^{-1}$ over 100 pkpc with lowest velocities east of SMM J02399$-$0136 and increasing to the southwest. We also discover a bright, Ly$\rm \alpha$ emitter, separated spatially and kinematically from the nebula, at a projected separation of $\approx$60 kpc from the quasar. This source has no clear central counterpart in deep $Hubble \, Space \, Telescope$ imaging, giving an intrinsic Ly$\rm \alpha$ rest-frame equivalent width greater than 312 $ \rm \mathring{A}$ (5-$\sigma$). We argue that this 'dark cloud' is illuminated by the quasar with a UV flux that is orders of magnitude brighter than the emission along our sightline. This result confirms statistical inferences that luminous quasars at $z>2$ emit UV radiation anisotropically. Future KCWI observations of other lines, e.g. Ly$\rm \beta$, He$\rm _{II}$, C$\rm _{IV}$, etc, and with polarimetry will further reveal the origin of the Ly$\rm \alpha$ nebula and nature of the dark cloud.

In this study, the fundamental stellar atmospheric parameters (Teff, log g, [Fe/H] and [{\alpha}/Fe]) were derived for low-resolution spectroscopy from LAMOST DR5 with Generative Spectrum Networks (GSN). This follows the same scheme as a normal artificial neural network with stellar parameters as the input and spectra as the output. The GSN model was effective in producing synthetic spectra after training on the PHOENIX theoretical spectra. In combination with Bayes framework, the application for analysis of LAMOST observed spectra exhibited improved efficiency on the distributed computing platform, Spark. In addition, the results were examined and validated by a comparison with reference parameters from high-resolution surveys and asteroseismic results. Our results show good consistency with the results from other survey and catalogs. Our proposed method is reliable with a precision of 80 K for Teff, 0.14 dex for log g, 0.07 dex for [Fe/H] and 0.168 dex for [{\alpha}/Fe], for spectra with a signal-to-noise in g bands (SNRg) higher than 50. The parameters estimated as a part of this work are available at http://paperdata.china-vo.org/GSN_parameters/GSN_parameters.csv.

The $\rm H_2$ mass of molecular clouds has traditionally been traced by the CO(J=1-0) rotational transition line. This said, CO is relatively easily photodissociated, and can also be destroyed by cosmic rays, thus rendering some fraction of molecular gas to be "CO-dark". We investigate the amount and physical properties of CO-dark gas in two $z \sim 0$ disc galaxies, and develop predictions for the expected intensities of promising alternative tracers ([CI 609 $\mu$m and [CII] 158 $\mu$m emission). We do this by combining cosmological zoom simulations of disc galaxies with thermal-radiative-chemical equilibrium interstellar medium (ISM) calculations to model the predicted H~\textsc{i} and $\rm H_2$ abundances and CO(J=1-0), [CI] 609 $\mu$m and [CII] 158 $\mu$m emission properties. Our model treats the ISM as a collection of radially stratified clouds whose properties are dictated by their volume and column densities, the gas-phase metallicity, and the interstellar radiation field and cosmic ray ionization rates. Our main results follow. Adopting an observationally motivated definition of CO-dark gas, i.e. $\rm H_2$ gas with $W_{\rm CO} < 0.1 $ K-km/s, we find that a significant amount ($> 50\%$) of the total $\rm H_2$ mass lies in CO-dark gas, most of which is diffuse gas, poorly shielded due to low dust column density. The CO-dark molecular gas tends to be dominated by [CII], though [CI] also serves as a bright tracer of the dark gas in many instances. At the same time, [CII] also tends to trace neutral atomic gas. As a result, when we quantify the conversion factors for the three carbon-based tracers of molecular gas, we find that [CI] suffers the least contamination from diffuse atomic gas, and is relatively insensitive to secondary parameters.

Observations show that winds can be driven from the innermost region (inside a 50 Schwarschild radius) of a thin disk. It is interesting to study the winds launched from the innermost region. A hot corona above the black hole (BH) thin disk is irradiated by the disk. We perform two-dimensional hydrodynamical simulations to study the winds driven by radiation force from the corona in the innermost regions. The hard X-ray spectrum from active galactic nuclei (AGNs) suggests that the corona temperature is about $10^9$ K, so that we mainly analyze the properties of winds (or outflows) from the $10^9$ K corona. The disk luminosity plays an important role in driving the outflows. The more luminous the disk, the stronger the outflows. Mass outflow rate ($\dot{M}_{\rm out}$) at a 90 Schwarschild radius depends on disk luminosity, which can be described as $\dot{M}_{\rm out}\propto 10^{3.3 \Gamma}$ ($\Gamma$ is the ratio of the disk luminosity to the Eddington luminosity). In the case of high luminosity (e.g. $\Gamma=0.75$), the supersonic outflows with maximum speed $1.0 \times 10^4$ Km s$^{-1}$ are launched at $\sim17^{o}$ --$30^{o}$ and $\sim50^{o}$ --$80^{o}$ away from the pole axis. The Bernoulli parameter keeps increasing with the outward propagation of outflows. The radiation force keeps accelerating the outflows when outflows move outward. Therefore, we can expect the outflows to escape from the BH gravity and go to the galactic scale. The interaction between outflows and interstellar medium may be an important AGN feedback process.

An unbiased search of molecular outflows within the region of the COHRS survey has identified 157 high-mass outflows from a sample of 770 ATLASGAL clumps with a detection rate of 20\%. The detection rate of outflows increases for clumps with higher M$_{clump}$, L$_{bol}$, L$_{bol}$/M$_{clump}$, N$_{H_{2}}$, and T$_{dust}$ compared to the clumps with no outflow. The detection rates of the outflow increases from protostellar (8\%) to YSO clump (17\%) and to MSF clump (29\%). The detection rate 26\% for quiescent clump is preliminary, because the sample of quiescent clumps is small. A statistical relation between the outflow and clump masses for our sample is $\log(M_{out}/M_{\bigodot}) = (-1.1\pm0.21) + (0.9\pm0.07)\log(M_{clump}/M_{\bigodot})$. The detection rate of outflows and the outflow mass-loss rate show an increase with increasing M$_{clump}$, L$_{bol}$, N$_{H_{2}}$, and T$_{dust}$, which indicates that clumps with outflow with higher parameter values are at a more advanced evolutionary stage. The outflow mechanical force increases with increasing bolometric luminosities. No clear evidence has yet been found that higher mass outflows have different launching conditions than low-mass outflows.

Stars are commonly formed in binary systems, which provide a natural laboratory for studying planet formation in extreme conditions. In our first paper (Paper I) of a series Xie et al. (2011), we have shown that the intermediate stage - from planetesimals to planetary embryos/cores - of planet formation can proceed even in highly inclined binaries. Following Paper I, here we numerically study the late stage of terrestrial planet formation, i.e., from embryos to full planets, in binary systems of various orbital configurations. We identify an orbital alignment or anti-alignment effect; namely, although an inclined binary generally misaligns the planetary orbits with respect to the spin axis of the primary host star (i.e., causing large obliquity), it could align or anti-align the planetary orbits with respect to the binary orbit. Such an orbital (anti-)alignment effect is caused by the combination of orbital differential precession and self-damping, and it is mostly significant in cases of intermediate binary separations, i.e., a_B~40-200 AU for terrestrial planet formation around 1 AU from the primary stars. In such intermediate separation binaries, somewhat contrary to intuition, the binary companion can aid planet growth by having increased the rate of collisions, forming significantly more massive but fewer planets. In the other two ends, the companion is ether too close thus plays a violently disruptive role or too wide to have significant effect on planet formation. Future observations, which can discover more planet-bearing binary star systems and constrain their masses and 3-D orbital motions will test our numerical findings.

Enormous Ly$\alpha$ nebulae (ELANe) represent the extrema of Ly$\alpha$ nebulosities. They have detected extents of $>200$ kpc in Ly$\alpha$ and Ly$\alpha$ luminosities $>10^{44}$ erg s$^{-1}$. The ELAN population is an ideal laboratory to study the interactions between galaxies and the intergalactic/circumgalactic medium (IGM/CGM) given their brightness and sizes. The current sample size of ELANe is still very small, and the few $z\approx2$ ELANe discovered to date are all associated with local overdensities of active galactic nuclei (AGNs). Inspired by these results, we have initiated a survey of ELANe associated with QSO pairs using the Palomar and Keck Cosmic Web Imagers (PCWI/KCWI). In this letter, we present our first result: the discovery of ELAN0101+0201 associated with a QSO pair at $z=2.45$. Our PCWI discovery data shows that, above a 2-$\sigma$ surface brightness of $1.2\times10^{-17}$ \sbunit, the end-to-end size of ELAN0101+0201 is $\gtrsim 232$ kpc. We have conducted follow-up observations using KCWI, resolving multiple Ly$\alpha$ emitting sources within the rectangular field-of-view of $\approx 130\times165$ projected kpc$^2$, and obtaining their emission line profiles at high signal-to-noise ratios. Combining both KCWI and PCWI, our observations confirm that ELAN0101+0201 resides in an extremely overdense environment. Our observations further support that a large amount of cool ($T\sim10^4$K) gas could exist in massive halos (M$\gtrsim10^{13}$M$_\odot$) at $z\approx2$. Future observations on a larger sample of similar systems will provide statistics of how cool gas is distributed in massive overdensities at high-redshift and strongly constrain the evolution of the intracluster medium (ICM).

We present a clear case study on the occurrence of two successive X-class flares including a decade-class flare (X9.3) and two coronal mass ejections (CMEs) triggered by shearing motion and sunspot rotation in active region NOAA 12673 on 2017 September 6. A shearing motion between the main sunspots with opposite polarities started on September 5 and even lasted after the second X-class flare on September 6. Moreover, the main sunspot with negative polarity rotated around its umbral center and another main sunspot with positive polarity also exhibited a slow rotation. The sunspot with negative polarity at the northwest of active region also began to rotate counter-clockwise before the onset of the first X-class flare. The successive formation and eruption of two S-shaped structures were closely related to the counter-clockwise rotation of three sunspots. It is also found that the rotation of sunspots is faster during four hours prior to the onset of the flares than the period before. The existence of a flux rope is found prior to the onset of two flares by using non-linear force free field extrapolation based on the vector magnetograms observed by SDO/HMI. These results suggest that shearing motion and sunspot rotation play an important role in the buildup of the free energy and the formation of flux ropes in the corona which produces solar flares and CMEs.

In this letter, we present a spectacular eruptive flare (X8.2) associated with a coronal mass ejection (CME) on 2017 September 10 at the west limb of the Sun. A flux rope eruption is followed by the inflow, the formation of a current sheet and a cusp structure, which were simultaneously observed during the occurrence of this flare. The hierarchical layers of the cusp-shaped structure are well observed in 131 \AA\ observation. The scenario that can be created from these observations is very consistent with the predictions of some eruptive models. Except for the characteristics mentioned above in the process of the flare predicted by classical eruption models, the current sheet separating into several small current sheets is also observed at the final stage of the flux rope eruption. The quantitative calculation of the velocities and accelerations of the inflow, hot cusp structure, and post-flare loops is presented. The width of the current sheet is estimated to be about 3 $\times$ $10^{3}$ km. These observations are very useful to understand the process of solar eruptions.

We analyze the high-resolution observations from the 1.6m telescope at Big Bear Solar Observatory that cover an active region filament. Counter-streaming motions are clearly observed in the filament. The northern end of the counter-streaming motions extends to a light bridge, forming a spectacular circulation pattern around a sunspot, with clockwise motion in the blue wing and counterclockwise motion in the red wing as observed in Halpha off-bands. The apparent speed of the flow is around 10 to 60 km/s in the filament, decreasing to 5 to 20 km/s in the light bridge. The most intriguing results are the magnetic structure and the counter-streaming motions in the light bridge. Similar to those in the filament, magnetic fields show a dominant transverse component in the light bridge. However, the filament is located between opposite magnetic polarities, while the light bridge is between strong fields of the same polarity. We analyze the power of oscillations with the image sequences of constructed Dopplergrams, and find that the filament's counter-streaming motion is due to physical mass motion along fibrils, while the light bridge's counter-streaming motion is due to oscillation in the direction along the line-of-sight. The oscillation power peaks around 4 minutes. However, the section of the light bridge next to the filament also contains a component of the extension of the filament in combination with the oscillation, indicating that some strands of the filament are extended to and rooted in that part of the light bridge.

We present a series of adaptive mesh refinement (AMR) hydrodynamic simulations of flat rotation curve galactic gas disks with a detailed treatment of the interstellar medium (ISM) physics of the atomic to molecular phase transition under the influence of diffuse FUV radiation fields and cosmic ray backgrounds. We explore the effects of different FUV intensities, including a model with a radial gradient designed to mimic the Milky Way. The effects of cosmic rays, including radial gradients in their heating and ionization rates, are also explored. The final simulations in this series achieve $4\:$pc resolution across the $\sim20\:$kpc global disk diameter, with heating and cooling followed down to temperatures of $\sim10\:$K. The disks are evolved for $300\:$Myr, which is enough time for the ISM to achieve a quasi-statistical equilibrium. In particular, the mass fraction of molecular gas stabilizes by $\sim$200 Myr. Additional global ISM properties are analysed. Giant molecular clouds (GMCs) are also identified and the statistical properties of their populations examined. GMCs are tracked as the disks evolve. GMC collisions, which may be a means of triggering star cluster formation, are counted and the rates compared with analytic models. Relatively frequent GMC collision rates are seen in these simulations and their implications for understanding GMC properties, including the driving of internal turbulence, are discussed.

We investigate collisions between giant molecular clouds (GMCs) as potential generators of their internal turbulence. Using magnetohydrodynamic (MHD) simulations of self-gravitating, magnetized, turbulent, GMCs, we compare kinematic and dynamic properties of dense gas structures formed when such clouds collide compared to those that form in non-colliding clouds as self-gravity overwhelms decaying turbulence. We explore the nature of turbulence in these structures via distribution functions of density, velocity dispersions, virial parameters, and momentum injection. We find that the dense clumps formed from GMC collisions have higher effective Mach number, greater overall velocity dispersions, sustain near-virial equilibrium states for longer times, and are the conduit for injection of turbulent momentum into high density gas at high rates.

Studying the dynamics of filaments at pre-eruption phase can shed light on the precursor of eruptive events. Such studies in high-resolution (in the order of 0.1") are highly desirable yet very rare so far. In this work, we present a detailed observation of a pre-eruption evolution of a filament obtained by the 1.6 m New Solar Telescope (NST) at Big Bear Solar Observatory (BBSO). One end of the filament is anchored at the sunspot in NOAA active region (AR) 11515, which is well observed by NST H$\alpha$ off-bands four hours before till one hour after the filament eruption. A M1.6 flare is associated with the eruption. We observed persistent downflowing materials along the H$\alpha$ multi-threaded component of the loop towards the AR end during the pre-eruption phase. We traced the trajectories of plasma blobs along the H$\alpha$ threads and obtained the plane-of-sky velocity of 45 km s$^{-1}$ on average. We further estimated the real velocities of the downflows and the altitude of the filament by matching the observed H$\alpha$ threads with magnetic field lines extrapolated from a nonlinear force-free field (NLFFF) model. Observation of chromospheric brightenings (BZs) at the footpoints of the falling plasma blobs is also presented in the paper. The lower limit of the kinetic energy per second of the downflows through the BZs is found to be $\sim$ 10$^{21}$ erg. Larger FOV observations from BBSO full disk H$\alpha$ images show that the AR end of the filament started ascending four hours before the flare. We attribute the observed downflows at the AR end of the filament to the draining effect of the filament rising prior to its eruption. During the slow-rise phase, the downflows continuously drained away $\sim$ 10$^{15}$g mass from the filament over a few hours, which is believed to be essential for the instability at last, and could be an important precursor of eruptive events.

We report new IRAM/PdBI, JCMT/SCUBA-2, and VLA observations of the ultraluminous quasar SDSSJ010013.02+280225.8 (hereafter, J0100+2802) at z=6.3, which hosts the most massive supermassive black hole (SMBH) of 1.24x10^10 Msun known at z>6. We detect the [C II] 158 $\mu$m fine structure line and molecular CO(6-5) line and continuum emission at 353 GHz, 260 GHz, and 3 GHz from this quasar. The CO(2-1) line and the underlying continuum at 32 GHz are also marginally detected. The [C II] and CO detections suggest active star formation and highly excited molecular gas in the quasar host galaxy. The redshift determined with the [C II] and CO lines shows a velocity offset of ~1000 km/s from that measured with the quasar Mg II line. The CO (2-1) line luminosity provides direct constraint on the molecular gas mass which is about (1.0+/-0.3)x10^10 Msun. We estimate the FIR luminosity to be (3.5+/-0.7)x10^12 Lsun, and the UV-to-FIR spectral energy distribution of J0100+2802 is consistent with the templates of the local optically luminous quasars. The derived [C II]-to-FIR luminosity ratio of J0100+2802 is 0.0010+/-0.0002, which is slightly higher than the values of the most FIR luminous quasars at z~6. We investigate the constraint on the host galaxy dynamical mass of J0100+2802 based on the [C II] line spectrum. It is likely that this ultraluminous quasar lies above the local SMBH-galaxy mass relationship, unless we are viewing the system at a small inclination angle.

The aim of this paper is to present an analytical relationship between the power spectral density of GRACE-like mission measurements and the accuracies of the gravity field coefficients mainly from the point of view of theory of signal and system, which indicates the one-to-one correspondence between spherical harmonic error degree variances and frequencies of the measurement noise. In order to establish this relationship, the average power of the errors due to gravitational acceleration difference and the relationship between perturbing forces and range-rate perturbations are derived, based on the orthogonality property of associated Legendre functions and the linear orbit perturbation theory, respectively. This method provides a physical insight into the relation between mission parameters and scientific requirements. By taking GRACE-FO as the object of research, the effects of sensor noises and time variable gravity signals are analyzed. If LRI measurements are applied, a mission goal with a geoid accuracy of 7.4 cm at a spatial resolution of 101 km is reachable, whereas if the KBR measurement error model is applied, a mission goal with a geoid accuracy of 10.2 cm at a spatial resolution of 125 km is reachable. Based on the discussion of the spectral matching of instrument accuracies, an improvement in accuracy of accelerometers is necessary for the match between the range errors and accelerometer noises in the future mission. Temporal aliasing caused by the time variable gravity signals is also discussed by this method.

G-band bright points (GBPs) are regarded as good manifestations of magnetic flux concentrations. We aim to investigate the relationship between the dynamic properties of GBPs and their longitudinal magnetic field strengths. High spatial and temporal resolution observations were recorded simultaneously with G-band filtergrams and Narrow-band Filter Imager (NFI) Stokes I and V images with Hinode /Solar Optical Telescope. The GBPs are identified and tracked in the G-band images automatically, and the corresponding longitudinal magnetic field strength of each GBP is extracted from the calibrated NFI magnetograms by a point-to-point method. After categorizing the GBPs into five groups by their longitudinal magnetic field strengths, we analyze the dynamics of GBPs of each group. The results suggest that with increasing longitudinal magnetic field strengths of GBPs correspond to a decrease in their horizontal velocities and motion ranges as well as by showing more complicated motion paths. This suggests that magnetic elements showing weaker magnetic field strengths prefer to move faster and farther along straighter paths, while stronger ones move more slowly in more erratic paths within a smaller region. The dynamic behaviors of GBPs with different longitudinal magnetic field strengths can be explained by that the stronger flux concentrations withstand the convective flows much better than weaker ones.

The Large sky Area Multi-Object Spectroscopic Telescope (LAMOST) General Survey is a spectroscopic survey that will eventually cover approximately half of the celestial sphere and collect 10 million spectra of stars, galaxies and QSOs. Objects both in the pilot survey and the first year general survey are included in the LAMOST First Data Release (DR1). The pilot survey started in October 2011 and ended in June 2012, and the data have been released to the public as the LAMOST Pilot Data Release in August 2012. The general survey started in September 2012, and completed its first year of operation in June 2013. The LAMOST DR1 includes a total of 1202 plates containing 2,955,336 spectra, of which 1,790,879 spectra have observed signal-to-noise S/N >10. All data with S/N>2 are formally released as LAMOST DR1 under the LAMOST data policy. This data release contains a total of 2,204,696 spectra, of which 1,944,329 are stellar spectra, 12,082 are galaxy spectra and 5,017 are quasars. The DR1 includes not only spectra, but also three stellar catalogues with measured parameters: AFGK-type stars with high quality spectra (1,061,918 entries), A-type stars (100,073 entries), and M stars (121,522 entries). This paper introduces the survey design, the observational and instrumental limitations, data reduction and analysis, and some caveats. Description of the FITS structure of spectral files and parameter catalogues is also provided.

While it is well recognized that interstellar grains are made of amorphous silicates and some form of carbonaceous materials, it remains debated regarding what exact chemical and physical form the carbonaceous component takes. Contemporary grain models assume that the silicate and carbon components are either physically separated, or they form a core-mantle structure, or they agglomerate to form porous composites. The core-mantle model posits that the mantle is made of some sort of aliphatic hydrocarbon materials and is responsible for the 3.4 micrometer absorption feature ubiquitously seen in the diffuse interstellar medium (ISM) of the Milky Way and external galaxies. This model is challenged by the nondetection of polarization in the 3.4 micrometer absorption feature as the 9.7 micrometer silicate feature is observed to be polarized. To alleviate this challenge, we calculate the degree of polarization of the 3.4 micrometer feature for spheroidal silicate dust coated by a layer of spherical aliphatic hydrocarbon. It is found that the 3.4 micrometer feature polarization still exceeds the observational upper limit, even though spherical aliphatic hydrocarbon mantles are expected to cause much less polarization than nonspherical (e.g., spheroidal) mantles. We have also shown that the composite grain model which consists of amorphous silicate, aliphatic hydrocarbon, and vacuum also predicts the 3.4 micrometer feature polarization to well exceed what is observed. These results support the earlier arguments that the aliphatic hydrocarbon component is physically separated from the silicate component unless the 3.4 micrometer absorption feature is just a minor carbon sink in the ISM.

We discuss the criteria to associate PSRs with SNRs, and summary 78 pairs of possible PSR-SNR associations which is the most complete sample so far. We refine them into three categories according to degree of reliability. Statistic study on PSR-SNR associations helps us understand massive star evolution and constrain pulsar's theory models.

We carried out intra-day variability (IDV) observations from August 2005 to January 2010 with the Urumqi 25m radio telescope for a dozen IDV sources including the quasar 0917+624. This target exhibited pronounced centimeter-band, intra-day variability during the 1980s--1990s, but its strong IDV phase ceased in 2000. The source showed no IDV in the majority of the Urumqi observing sessions, although weak IDV activity was detected in some. Multifrequency UMRAO data for 0917+624 show that the spectral index is steeper during the weak and non-IDV phases than during the strong IDV phase, supporting the idea that the size of the scintillating component may be enlarged in the weak/non IDV phases.

We aim to investigate the relation between the long-term flux density and the position angle (PA) evolution of inner-jet in blazars. We have carried out the elliptic Gaussian model-fit to the `core' of 50 blazars from 15 GHz VLBA data, and analyzed the variability properties of three blazars from the model-fit results. Diverse correlations between the long-term peak flux density and the PA evolution of the major axis of the `core' have been found in $\sim$ 20% of the 50 sources. Of them, three typical blazars have been analyzed, which also show quasi-periodic flux variations of a few years (T). The correlation between the peak flux density and the PA of inner-jet is positive for S5~0716+714, and negative for S4~1807+698. The two sources cannot be explained with the ballistic jet models, the non-ballistic models have been analyzed to explain the two sub-luminal blazars. A correlation between the peak flux density and the PA (with a T/4 time lag) of inner-jet is found in [HB89]~1823+568, this correlation can be explained with a ballistic precession jet model. All the explanations are based mainly on the geometric beaming effect; physical flux density variations from the jet base would be considered for more complicated situations in future, which could account for the no or less significance of the correlation between the peak flux density and the PA of inner-jet in the majority blazars of our sample.

X-ray loud M dwarfs are a major source of by-product (contamination) in the X-ray band of the multiwavelength quasar survey (MWQS). As a by-product, the low dispersion spectra of 22 M dwarfs are obtained in which the spectra of 16 sources are taken for the first time. The spectral types and distance of the sample are given based on spectral indices CaH2, CaH3, and TiO5. The parameter {\zeta}TiO/CaH is calculated to make the metallicity class separation among dwarfs, subdwarfs and extreme subdwarfs. We also discuss the distributions in the diagrams of Log(Lx/Lbol) versus spectral type and infrared colors.

CdZnTe (CZT) is a wide bandgap II-VI semiconductor developed for the spectroscopic detection of X-rays and {\gamma}-rays at room temperature. The Swift Burst Alert Telescope is using an 5240 cm2 array of 2 mm thick CZT detectors for the detection of 15-150 keV X-rays from Gamma-Ray Bursts. We report on the systematic tests of thicker (\geq 0.5 cm) CZT detectors with volumes between 2 cm3 and 4 cm3 which are potential detector choices for a number of future X-ray telescopes that operate in the 10 keV to a few MeV energy range. The detectors contacted in our laboratory achieve Full Width Half Maximum energy resolutions of 2.7 keV (4.5%) at 59 keV, 3 keV (2.5%) at 122 keV and 4 keV (0.6%) at 662 keV. The 59 keV and 122 keV energy resolutions are among the world-best results for \geq 0.5 cm thick CZT detectors. We use the data set to study trends of how the energy resolution depends on the detector thickness and on the pixel pitch. Unfortunately, we do not find clear trends, indicating that even for the extremely good energy resolutions reported here, the achievable energy resolutions are largely determined by the properties of individual crystals. Somewhat surprisingly, we achieve the reported results without applying a correction of the anode signals for the depth of the interaction. Measuring the interaction depths thus does not seem to be a pre-requisite for achieving sub-1% energy resolutions at 662 keV.

We present new imaging data and archival multiwavelength observations of the little studied emission nebula K 1-6 and its central star. Narrow-band images in H-alpha (+ [NII]) and [OIII] taken with the Faulkes Telescope North reveal a stratified, asymmetric, elliptical nebula surrounding a central star which has the colours of a late G- or early K-type subgiant or giant. GALEX ultraviolet images reveal a very hot subdwarf or white dwarf coincident in position with this star. The cooler, optically dominant star is strongly variable with a period of 21.312 +/- 0.008 days, and is possibly a high amplitude member of the RS CVn class, although an FK Com classification is also possible. Archival ROSAT data provide good evidence that the cool star has an active corona. We conclude that K 1-6 is most likely an old bona fide planetary nebula at a distance of ~1.0 kpc, interacting with the interstellar medium, and containing a binary or ternary central star. The observations and data analyses reported in this paper were conducted in conjunction with Year 11 high school students as part of an Australian Research Council Linkage Grant science education project, denoted Space To Grow, conducted jointly by professional astronomers, educational researchers, teachers, and high-school students.

We report on the design and estimated performance of a balloon-borne hard X-ray polarimeter called HX-POL. The experiment uses a combination of Si and Cadmium Zinc Telluride detectors to measure the polarization of 50 keV-400 keV X-rays from cosmic sources through the dependence of the angular distribution of Compton scattered photons on the polarization direction. On a one-day balloon flight, HX-POL would allow us to measure the polarization of bright Crab-like sources for polarization degrees well below 10%. On a longer (15-30 day) flight from Australia or Antarctica, HX-POL would be be able to measure the polarization of bright galactic X-ray sources down to polarization degrees of a few percent. Hard X-ray polarization measurements provide unique venues for the study of particle acceleration processes by compact objects and relativistic outflows. In this paper, we discuss the overall instrument design and performance. Furthermore, we present results from laboratory tests of the Si and CZT detectors.

The Chinese Small Telescope ARray (hereinafter CSTAR) is the first Chinese astronomical instrument on the Antarctic ice cap. The low temperature and low pressure testing of the data acquisition system was carried out in a laboratory refrigerator and on the 4500m Pamirs high plateau, respectively. The results from the final four nights of test observations demonstrated that CSTAR was ready for operation at Dome A, Antarctica. In this paper we present a description of CSTAR and the performance derived from the test observations.

The quality of Cadmium Zinc Telluride (CZT) detectors is steadily improving. For state of the art detectors, readout noise is thus becoming an increasingly important factor for the overall energy resolution. In this contribution, we present measurements and calculations of the dark currents and capacitances of 0.5 cm-thick CZT detectors contacted with a monolithic cathode and 8x8 anode pixels on a surface of 2 cm x 2 cm. Using the NCI ASIC from Brookhaven National Laboratory as an example, we estimate the readout noise caused by the dark currents and capacitances. Furthermore, we discuss possible additional readout noise caused by pixel-pixel and pixel-cathode noise coupling.

The polarimetric and photometric variability of Wolf-Rayet (WR) stars as caused by clumps in the winds, is revisited. In the model which is improved from Li et al. 2000, the radial expansion of the thickness is accounted for, but we retain the dependence on the beta velocity law, stellar occultation effects. We again search for parameters that can yield results consistent with observations in regards to the mean polarization, the ratio of polarimetric to photometric variability, and the volume filling factor. Clump generation and spatial distribution are randomized by the Monte Carlo method so as to produce clumps which are, in the mean, distributed uniformly in space and have time intervals with a Gaussian distribution. The generated clumps move radially outward with a velocity law determined by a beta index, and the angular size of the clumps is assumed to keep fixed. By fitting the observational results and the volume filling factor, the clump velocity law index beta and clump ejection rate are inferred, and are found to be well constrained. In addition, the subpeak features on broad emission lines seem to support the clump ejection rate. Meanwhile, the fraction of the total mass loss rate that is contained in the clumps is obtained by fitting the observed polarization. We conclude that this picture for the clump properties produces a valuable diagnostic of WR wind structure.

The main methods grown Cadmium Zinc Telluride (CZT) crystals with high yield and excellent homogeneity are Modified Horizontal Bridgman (MHB) and High Pressure Bridgman (HPB) processes, respectively. In this contribution, the readout system based on two 32-channel NCI-ASICs for pixellated CZT detector arrays has been developed and tested. The CZT detectors supplied by Orbotech (MHB) and eV products (HPB) are tested by NCI-ASIC readout system. The CZT detectors have an array of 8x8 or 11x11 pixel anodes fabricated on the anode surface with the area up to 2 cm x2 cm and the thickness of CZT detectors ranges from 0.5 cm to 1 cm. Energy spectra resolution and electron mobility-lifetime products of 8x8 pixels CZT detector with different thicknesses have been investigated.

The Modified Horizontal Bridgman (MHB) process produces Cadmium Zinc Telluride (CZT) crystals with high yield and excellent homogeneity. Various groups,including our own, previously reported on the test of 2x2x0.5 cm3 MHB CZT detectors grown by the company Orbotech and read out with 8x8 pixels. In this contribution, we describe the optimization of the photolithographic process used for contacting the CZT detector with pixel contacts. The optimized process gives a high yield of good pixels down to pixel diameters/pitches of 50 microns. Furthermore, we discuss the performance of 0.5 cm and 0.75 cm thick detectors contacted with 64 and 225 pixel read out with the RENA-3 ASICs from the company NOVA R&D.

We focus attention on the Oe/Be stars to test the concept that the disks of these stars form by magnetic channeling of wind material toward the equator. Calculations are made of the X-rays expected from the Magnetically Torqued Disk (MTD) model for Be stars discussed by Cassinelli et al. (2002), by Maheswaran (2003), and by Brown et al. (2004). The dominant parameters in the model are the $\beta$ value of the velocity law, the rotation rate of the star, $S_o$, and the ratio of the magnetic field energy density to the disk gravitational energy density, $\gamma$. The model predictions are compared with the $ROSAT$ observations obtained for an O9.5 star $\zeta$ Oph from \Berghofer et al. (1996) and for 7 Be stars from Cohen et al. (1997). Extra considerations are also given here to the well studied Oe star $\zeta$ Oph for which we have $Chandra$ observations of the X-ray line profiles of the triad of He-like lines from the ion Mg XI.

In the magnetically torqued disk (MTD) model for hot star disks, as proposed and formulated by \citet{cassi02}, stellar wind mass loss was taken to be uniform over the stellar surface. Here account is taken of the fact that as stellar spin rate S_o (=\sqrt {\Omega_o^2 R^3/GM}) is increased, and the stellar equator is gravity darkened, the equatorial mass flux and terminal speed are reduced, compared to the poles, for a given total \mdot. As a result, the distribution of equatorial disk density, determined by the impact of north and southbound flows, is shifted further out from the star. This results, for high S_o (\gtrsim 0.5), in a fall in the disk mass and emission measure, and hence in the observed emission line EW, scattering polarization, and IR emission. Consequently, contrary to expectations, critical rotation S_o \to 1 is not the optimum for creation of hot star disks which, in terms of EM for example, is found to occur in a broad peak around S_o\approx 0.5-0.6 depending slightly on the wind velocity law. The relationship of this analytic quasi-steady parametric MTD model to other work on magnetically guided winds is discussed. In particular the failures of the MTD model for Be-star disks alleged by \citet{owo03} are shown to revolve largely around open observational tests, rather in the basic MTD physics, and around their use of insufficiently strong fields.

The hot star wind momentum problem \eta = mdot*vinfty/(L/c) >> 1 is revisited, and it is shown that the conventional belief, that it can be solved by a combination of clumping of the wind and multiple scattering of photons, is not self-consistent for optically thick clumps. Clumping does reduce the mass loss rate mdot, and hence the momentum supply, required to generate a specified radio emission measure \varepsilon, while multiple scattering increases the delivery of momentum from a specified stellar luminosity L. However, in the case of thick clumps, when combined the two effects act in opposition rather than in unison since clumping reduces multiple scattering. Thus the ratio of momentum deliverable to momentum required is maximal for a smooth wind and the momentum problem remains for the thick clump case. In the case of thin clumps, all of the benefit of clumping in reducing \eta lies in reducing \mdot for a given \varepsilon so that extremely small filling factors f = 10e-4 are needed. It is also shown that clumping affects the inference of \mdot from radio \varepsilon not only by changing the emission measure per unit mass but also by changing the radio optical depth unity radius R_{rad}, and hence the observed wind volume, at radio wavelengths.

The low dispersion (400$\AA$/mm) spectrum of the optical counterpart of a flat-spectrum radio source 87GB 080315.5+512613, which is one of two possible radio counterparts of 2EG J0809+5117, was obtained recently. The optical counterpart, which is $2.02''$ away from 87GB 080315.5+512613 and $19.3'$ away from 2EG J0809+5117, was identified as a quasar with redshift of 1.14. We noted that Mattox et al. (1997) suggested the other radio counterpart 87GB 080459.4+495915 (OJ 508), which is $87.1'$ away from 2EG J0809+5117, is the more potential identification, though it was previously suggested to be the identification (with low confidence) of another nearby EGRET source 2EG J0807+4849. Our observation suggests that it is quite possible that 87GB 080315.5+512613 is the identification of 2EG J0809+5117 rather than 87GB 080459.4+495915. But we still can not exclude the possibility of 87GB 080459.4+495915 at present. Moreover, in order to determine whether or not 87GB 080315.5+512613 is a blazar type quasar, the optical polarization and variability measures of its optical counterpart are strongly encouraged.

The local stability of accretion disks with advection is studied together with the considerations of radial viscous force and thermal diffusion. For a geometrically thin, radiative cooling dominated disk, the thermal diffusion has nearly no effects on the thermal and viscous modes. The including of thermal diffusion, however, tends to stabilize the acoustic modes which, if without advection, are unstable if the disk is optically thick, radiation pressure dominated or optically thin, and are stable if the disk is optically thick, gas pressure dominated. The including of very little advection has significant effects on two acoustic modes. Independent on the optical depth, the instability of the outward propagating mode (O-mode) is enhanced and that of the inward propagating mode (I-mode) is damped if the disk is gas pressure dominated, while the instability of O-mode is damped and that of I-mode is enhanced if the disk is radiation pressure dominated. For a geometrically slim, advection-dominated disk, both the thermal and viscous modes, as well as I-mode, are always stable if the disk is optically thin. The including of thermal diffusion tends to make these modes more stable. However, the O-mode can become unstable when $q/m$ is very large ($q$ is the ratio of advective to viscous dissipated energy and $m$ the Mach number), even if the thermal diffusion is considered. On the other hand, if the advection-dominated disk is optically thick, we found there is no self-consistent acoustic modes in our local analyses. The thermal diffusion has no effect on the stable viscous mode but has a significant contribution to enhance the thermal instability.